陈 梦 黄 翔 屈悦滢
(西安工程大学城市规划与市政工程学院 西安 710048)
蒸发冷却空调技术以水作为制冷剂,最大限度地利用了我国干燥地区丰富的干空气能源,能有效减少温室气体排放、节约能源[1]。作为我国暖通空调行业绿色、环保、节能的关键技术之一,目前蒸发冷却空调技术已日益成熟并实现了产业化,已用于工业建筑、公共建筑、住宅建筑、轨道交通、农业建筑的空调系统中,可以大大减少峰值的电力消耗和电力需求[2]。
露点间接蒸发冷却技术利用一次空气干球温度和露点温度的差值为驱动势,能有效提高蒸发冷却空调设备的温降。由于其显著的优势,成为目前国际上许多专家和学者研究的热点。关于露点间接蒸发冷却性能预测方面还鲜有研究,如何寻找一种简单、快速、高效的性能预测模型已成为研究重点之一。采用神经网络方法的优势在于可以充分模拟预测、规避产品缺陷;降低机组的研发成本、缩短机组研发的周期、增强机组的核心竞争力[3]。
本项目用的露点间接蒸发冷却空调机组位于陕西省榆林市锦界镇,机组实物图如图1所示,机组测点布置图如图3所示,测点分别布置在机组的进风口、排风口、送风口处。其中1 台机组应用于中控室,面积约50m2/间;2 台应用于电子间,面积约60m2/间,机组风量均为20000m3/h。室内还装有三台机械制冷柜机,作为辅助冷源使用。
图1 露点间接蒸发冷却机组实物图Fig.1 Physical map of dew point indirect evaporative cooling air conditioning unit
图2 露点间接蒸发冷却空调系统平面图Fig.2 Dew point indirect evaporative cooling air conditioning unit system plan
图3 露点间接蒸发冷却机组测点布置示意图Fig.3 Dew point indirect evaporative cooling unit measurement point arrangement diagram
本项目用的是在福州地区某公司实验室安装的露点机组[4],机组实物图如图4所示,机组测点布置图同样如图3所示,测点分别布置在机组的进风口、排风口、送风口处。该机组额定风量15000m3/h,对上述机组进行测试,机组的风口截面积不变,实测的机组进风迎面平均风速为5.5m/s,二/一次风量比为0.6,一次风量为6375m3/h,二次风量为3825m3/h。部分测试数据如表2所示。
图4 机组实物图Fig.4 Physical map of the unit
表2 15000m3/h 风量机组测试数据(部分数据)Table 2 Test data of 15000m3/h air volume unit(partial data)
本工程为机场扩建工程子项目,主要范围为敦煌机场航空配餐中心。总建筑面积2880m2,共2层,层高4.2m,机组实物图如图5所示,共有6台机组[5,6],机组测点布置图同样如图3所示,测点分别布置在机组的进风口、排风口、送风口处。该机组额定风量分别为20000m3/h、10000m3/h,对上述机组进行测试,机组的风口截面积不变,实测的机组进风迎面平均风速分别为6.1m/s、4.7m/s,二/一次风量比分别为0.56、0.62,一次风量分别为11804m3/h、4320m3/h,二次风量分别为6610m3/h、2680m3/h。部分测试数据如表3所示。
图5 机组实物图Fig.5 Physical map of the unit
表3 10000m3/h 和20000m3/h 风量机组测试数据(部分数据)Table 3 Test data of 10000m3/h and 20000m3/h air volume unit(partial data)
续表3 10000m3/h 和20000m3/h 风量机组测试数据(部分数据)
2.1.1 BP 神经网络算法介绍
BP 算法包括样本正向、误差反向传播两部分。如图6所示,算法结构分为三部分[7]。当网络的结构被确定后,根据输入、输出数据对网络进行训练,找到两者之间的关系,修正隐含层的阈值。
图6 BP 神经网络拓扑结构图Fig.6 BP Neural network topology structure diagram
2.1.2 BP 神经网络研究方法
(1)训练样本的选取
训练样本:本文所研究的机组功能段核心为露点间接段。由于影响机组冷却效率的因素对机组性能产生的影响是非线性的,同时需要较多的测试数据进行训练,才能更好地寻找两者之间的规律。因此,采用机组实测数据作为网络的训练样本,利用Matlab 软件建立预测露点蒸发冷却空调机组性能的神经网络模型,从而利用训练好的神经网络实现对露点蒸发冷却空调机组的蒸发效率进行预测。
(2)输入/输出样本设计
影响露点间接蒸发冷却效率的影响因素很多,本课题主要研究在机组其他条件不变的情况下,主要考虑机组进风口干球温度、相对湿度、总风量、迎面风速、一次风量、二次风量、二/一次风量比对机组性能(蒸发效率)的影响。为了更有效地进行网络训练,对训练数据归一化[8],对网络输出数据进行反归一化。
(3)神经网络的结构设计
输入层包括机组进风口的空气状态参数(干球温度、相对湿度)、总风量、迎面风速、一次风量、二次风量、二/一次风量比7 个重要的参数,由于考虑到多输出参数导致的神经网络参数权重以及参数之间共用最优解的问题,因此输出层有蒸发效率1 个参数。网络采用一个隐含层,隐含层的神经元个数按Mr.Kawashima 提出的方法选择,即神经元个数为7×2+2=16,即隐含层神经元个数初始值取16 个,通过程序运行过程中采用试凑法得出最佳的隐含层神经元个数为14。
(4)网络传递及训练函数
隐含层神经元采用对数S 型正切函数Tansig,输出层神经元采用purelin 函数。为避免BP 神经网络算法的预测结果陷入局部最优的问题,学习和训练函数采用LM 算法的trainlm 函数。同时还需要设置训练次数1000 次,训练误差10-9,学习率0.1。
本小节对中等湿度地区(陕西榆林)露点间接蒸发冷却空调机组性能预测,测试时间为2019年8月1—6日、每天9:00—18:00 的实际机组运行的测试数据,每小时记录一次,每天10 组,共有60组测试数据。以2019年8月1—5日每天9:00—18:00 的测试数据为训练数据,设置2019年8月6日9:00—18:00 的数据为预测数据。
由图7 和图8,BP 网络对机组蒸发效率的预测值与实测值的变化趋势整体上基本一致。但存在个别数据误差较大的问题,其中,在测试样本号为7 的相对误差较大,接近10%。通过对预测结果进行检验与验证,得到该预测模型蒸发效率的RMSE为3.1367,决定系数r2为0.9659。
图7 BP 神经网络对机组蒸发效率的预测结果Fig.7 BP Neural Network's Forecast Results of Unit Evaporation Efficiency
图8 BP 神经网络对机组蒸发效率的预测相对误差Fig.8 BP neural network predicts relative error of unit evaporation efficiency
本小节对高湿度地区(福建福州)机组的性能预测方法同上,并且机组结构相同,机组测试时间为2017年6月23日、27—30日,通过对数据的整理,选择每天9:00—18:00 的机组运行的测试数据,每30 分钟记录一次,每天20 组,共有100 组测试数据。以2017年6月23日、27—29日每天9:00—18:00 的测试数据为训练数据,设置2017年6月30日9:00—18:00 的数据为预测数据。
由图9 和图10,BP 网络对机组蒸发效率的预测值与实测值的变化趋势整体上基本一致。但存在个别数据误差较大的问题,其中,在测试样本号为3、10、16 的相对误差较大,最大的相对误差接近10%。通过对预测结果进行检验与验证,得到该预测模型蒸发效率的RMSE 为4.6170,决定系数r2为0.9581。因此可认为,BP 网络能够预测出蒸发效率及其变化趋势。
图9 BP 神经网络对机组蒸发效率的预测结果Fig.9 BP Neural Network's Forecast Results of Unit Evaporation Efficiency
图10 BP 神经网络对机组蒸发效率的预测相对误差Fig.10 BP neural network predicts relative error of unit evaporation efficiency
本小节对干燥地区(敦煌机场)机组性能预测方法同上,并且机组结构相同,机组测试时间为2018年8月19-23日,通过对数据的整理,选择每天9:00—18:00 的机组运行的测试数据,每30分钟记录一次,每天20 组,共有100 组测试数据。以2018年8月19—22日每天9:00—18:00 的测试数据为训练数据,设置2018年8月23日9:00—18:00 的数据为预测数据。
由图11 和图12,BP 网络对机组蒸发效率的预测值与实测值的变化趋势整体上也基本一致。但同时存在个别数据误差较大的问题,其中,在测试样本号为4 的相对误差较大,接近14%。通过对预测结果进行检验与验证,得到该预测模型蒸发效率的RMSE 为3.4934,决定系数r2为0.9673。因此可认为,BP 网络能够准确的预测出蒸发效率及其变化趋势。
图11 BP 神经网络对机组蒸发效率的预测结果Fig.11 BP Neural Network's Forecast Results of Unit Evaporation Efficiency
图12 BP 神经网络对机组蒸发效率的预测相对误差Fig.12 BP neural network predicts relative error of unit evaporation efficiency
本文2.2-2.4 节分别是对不同地区露点间接蒸发冷却空调机组的性能预测,本节神经网络输入输出样本为三个地区的260 组数据,设置陕西地区2019年8月1—2日9:00—18:00 的数据为预测值,其他数据均为模型的训练值。
由图13 和图14,BP 网络对机组蒸发效率的预测值与实测值的变化趋势整体上是基本一致的。同样存在个别数据误差较大的问题,其中,在测试样本号为11、12 的相对误差较大,接近13%。通过对预测结果进行检验与验证,得到该预测模型蒸发效率的RMSE 为6.4157,决定系数r2为0.9381。因此可认为,BP 网络能够准确的预测出机组蒸发效率及其变化趋势。
图13 BP 神经网络对机组蒸发效率的预测结果Fig.13 BP Neural Network's Forecast Results of Unit Evaporation Efficiency
图14 BP 神经网络对机组蒸发效率的预测相对误差Fig.14 BP neural network predicts relative error of unit evaporation efficiency
由表4 可得,对于BP 神经网络在不同气象条件下露点间接蒸发冷却空调机组的预测,决定系数r2均在0.9 以上,表明所搭建的模型具有一定的稳定性。均方根误差RMSE 平均为4.4157。
表4 BP 神经网络对不同气象条件露点间接蒸发冷却空调机组蒸发效率的预测结果对比Table 4 Comparison of prediction results of BP neural network on the evaporation efficiency of indirect evaporative cooling air conditioning units with dew point in different meteorological conditions
2.6.1 灰色关联分析方法理论基础
灰色关联分析是一种多因素统计方法,衡量各个因素间的关联度,适用于动态过程分析[9],可以弥补数理统计方法的不足,具有很强的适用性。该方法通过对数据进行量化,计算量小,非常方便,因此灰色关联分析的结果不会与定性结果矛盾。本文采用该方法来实现多个因素对结果影响程度的分析。
2.6.2 输入变量对输出变量影响程度结果分析
采用Matlab 来实现输入变量与输出变量之间的关联度分析。得到的关联度如表5 和图15所示,通过比较7 个输入变量与输出变量关联度的大小可以看出:二/一次风量比>迎面风速>进风温度>总风量>相对湿度>二次风量>一次风量。
图15 关联度结果图Fig.15 Relevance result graph
表5 关联度的计算结果Table 5 Calculation result of relevance
因此得出结论:输出变量蒸发效率受到输入变量二/一次风量比的影响最大。二次风量的影响虽然是最小的,但是关联度为0.5657>0.5,因此也是重要的参数。
由于第2 节所建立的BP 神经网络模型预测效果具有较好的稳定性,在不同气象条件下(陕西、福建、甘肃)的预测效果差别不大,因此本节将建立的PSO-BP 神经网络预测模型仅应用在陕西榆林,即本节采用的数据为第1.2 节露点间接蒸发冷却空调机组的实测数据。训练数据和预测数据的设置同2.2 节。
PSO 优化BP 神经网络模型是对网络的阈值、权值进行不断更新的一个过程,通过构建合适的网络结构来对机组的蒸发效率进行预测。网络的结构设计同2.1 节。优化算法的流程图如图16所示。
图16 PSO 优化BP 神经网络模型的流程图Fig.16 Flow chart of PSO optimized BP neural network model
(1)确定粒子的维度
式中:m为输入层神经元数;n为隐含层神经元数;l为输出层神经元数。
根据公式(1),由2.1 节可知BP 神经网络模型的输入层神经元数为7、隐含层神经元数为14、输出层神经元数为1,因此可得D=127。
(2)确定粒子种群规模N
N与该算法的计算量和全局寻优能力有很大关系,N越大,参与搜索的粒子数越多,所以算法全局寻优能力越强。N过大也会导致其计算时间呈现指数式增长。N过少会使全局信息减少从而易陷入局部最优。N的增长不与搜索能力成正比[10]。一般来说,经验参数值多在[20,80]。综上,N取50。
(3)确定学习因子C1和C2
一般情况下C1和C2为非负常数。PSO 算法的两个学习因子代表的是单个粒子的认知能力与整个群体的信息交换能力,C1和C2可以确定粒子在整体信息扩展和自我认知学习之间的关系,决定粒子的搜索轨迹。C1和C2越大越容易让算法快速达到最优解范围,若C1和C2过大,解会在最优解中徘徊,不能获取最优;过小则会导致算法缓慢,得不到最优解。一般经验取值在[0,4]。
因此本节在建立PSO 优化BP 神经网络算法时,设定C1=C2=2。
(4)确定惯性权重ω
ω表示粒子在每次迭代前的惯性,不断寻找最优解。ω越大粒子移动速度越大,搜索速度变快,在全局的寻优能力越强,有时候会使算法无法收敛,因此就需要进行合理的取值。ω越小粒子的移动速度越慢,可能会导致算法的局部收敛。可以分为固定和时变两类,前者的值是固定的,而后者可以针对实际情况自动进行及时的调整,在一定程度上能够提高算法准确性。由于其弊端是在高维度粒子群里,运行速度减慢。因此本节选用固定权重。一般经验取值ω范围为[0.8,1.2],所以取ω=1。
(5)确定粒子的搜索范围和速度范围
粒子的搜索能力与速度范围、搜索范围相关。搜索范围越大,算法的精度越高,但是运行速度会越低。搜索范围越小,运行速度加快,陷入局部最小解。由于BP 神经网络的权值、阈值的范围为[0,1],所以设置PSO 的搜索范围和速度范围也为[0,1]。
(6)PSO-BP 模型输入参数的确定
根据对影响蒸发效率的影响因素考虑,本节依然选取机组进风口干球温度、相对湿度、总风量、迎面风速、一次风量、二次风量、二/一次风量比这7 个参数作为PSO-BP 神经网络的输入变量,神经网络的输出变量为机组的蒸发效率。
3.2.1 PSO-BP 神经网络模型预测结果
利用采集并处理过的测试数据对优化后的网络进行训练和测试分析,以此来验证优化后网络的预测效果。图17、图18 分别为两者的预测结果、相对误差,可以看出,PSO-BP 神经网络对机组蒸发效率的预测值与测试变化趋势基本一致,且两者的数据较为接近,不存在个别数据误差较大的问题,最大相对误差2.48%。
图18 PSO-BP 神经网络预测相对误差Fig.18 PSO-BP Neural network predicts relative error
通过对该PSO-BP神经网络的预测结果进行检验与验证,得到该预测模型蒸发效率的RMSE 为1.102,决定系数r2为0.9745。因此可认为,本节所建立的PSO-BP 神经网络模型精度相对较高,该PSO-BP 神经网络能够准确预测蒸发效率及其变化趋势。
3.2.2 BP、PSO-BP 神经网络模型预测结果对比分析
图19 为BP 神经网络和优化后的PSO-BP 神经网络预测露点间接蒸发冷却空调机组蒸发效率的结果对比图,可以看出,优化后的模型预测曲线与实测值更接近,没有出现太大的误差波动点,预测精度明显要比BP 神经网络好。
图19 BP、PSO-BP 神经网络的预测结果对比图Fig.19 Comparison chart of prediction results of BP and PSO-BP neural network
从图20 和表6 可以明显看出,优化后模型的相对误差整体较优化前的BP 神经网络小,因此可以用于预测露点间接蒸发冷却空调机组的蒸发效率。且在同一时刻,优化前的BP 神经网络相对误差为9.6%,优化后仅为2.21%,准确度有效提高了7.39%。
图20 BP、PSO-BP 神经网络预测相对误差对比Fig.20 Comparison of relative error of BP and PSO-BP neural network prediction
表6 BP、PSO-BP 神经网络预测结果对比Table 6 Comparison of prediction results of BP and PSO-BP neural network
(1)数据的采集:从三个不同地区整理的测试数据来看,机组送风参数均满足室内送风要求,说明应用露点间接蒸发冷却空调机组能保证室内空气调节要求。采集到了可供分析应用的基础数据,为神经网络预测和对比奠定基础。
(2)参数的选取:在机组其他条件不变的情况下,仅考虑机组进风参数干球温度、进风相对湿度、总风量、迎面风速、一次风量、二次风量、二/一次风量比7 个重要参数对机组蒸发效率的影响。采用灰色关联分析方法,得出蒸发效率受到二/一次风量比的影响最大,受到二次风量的影响最小,但关联度为0.5657>0.5,因此也是重要的参数。
(3)模型的研究:虽然BP 神经网络对机组蒸发效率的预测结果在可接受的范围内,预测结果与实测数据虽然规律相近,但是有个别点存在误差过大的问题,因此必须对该模型进行改进。
(4)模型的优化:实验表明PSO 算法优化后的BP 神经网络预测模型在预测准确度和精度上都优于BP 神经网络预测模型。在同一时刻,BP 神经网络预测的相对误差为9.6%时,PSO-BP 神经网络预测的相对误差为2.21%,精度提高了7.39%。