[发明专利]一种温室能耗的模型优化预测方法

摘要

一种温室能耗的模型优化预测方法，包括以下步骤：1)建立温室能量物理模型，温室热平衡定义为温度变化速率，同时温度变化速率转换来计算温室内部供热系统需要提供的能量，供热系统采用风机盘管；2)确定温室中需要辨识的温室因子，目标函数被表示为实际能量消耗和模拟能量消耗的均方根误差RMSE；3)模型参数IPSO‑GA优化。本发明提供一种有效降低计算时间、提高精度的一种温室能耗的模型优化预测方法。

主权项

1.一种温室能耗的模型优化预测方法，其特征在于：所述预测方法包括以下步骤：1)建立温室能量物理模型温室热平衡定义为温度变化速率，同时通过温度变化速率转换来计算温室内部供热系统需要提供的能量，供热系统采用风机盘管，风机盘管需要提供的能量表述如下：式中，Q_g为风机盘管提供的能量，W；ρ_air为空气密度，kg/m³；V为温室的体积，m³；C_air为空气的比热容，J/(kg·k)；T_i为温室室内的空气温度，℃，为单位时间内温度的变化率；Q_solar为进入到温室的净太阳光照辐射能量，W；Q_long为进入到温室的长波辐射能量，W；Q_cover为温室与外界的传导热量，W；Q_vent为温室与外界的通风换热能量，W；Q_trans为温室中的作物蒸腾所需要的能量，W；Q_crop为温室中的作物冠层与空气热传导的能量，W；根据热辐射定律，进入到温室的净太阳光照辐射能量Qsolar表示如下：Qsolar＝AsIaTa·TcTnXt               (3)式中，As为温室地表面积，m2；Ia为室外辐射通量密度，W/m2；Ta为覆盖材料的透光率，Tc为保温幕透光率；Tn为内遮阳透光率；Xt为净太阳光照辐射的修正系数；温室的长波辐射能量Qlong表示如下：Qlong＝ε12Agσ(Ti4‑Tsky4)·Xl        (7)式中，ε12为覆盖材料和空气之间的发射率，分别由空气和覆盖材料各自的发射率系数ε1和ε2决定；Ag为温室覆盖材料的表面积，m2；σ为斯蒂芬‑玻耳兹曼常数；Tsky为天空温度，K；Xl为长波辐射的修正系数；温室与外界的传导热量Qcover表示如下：Qcover＝AtopGaKcKn·(Ti‑To)+S1(Ti‑To)·Ga+S2(Ti‑Tog1)·Kg+S3(Ti‑Tog2)·Kg+S4(Ti‑Tog3)·Kg      (9)式中，Kg为覆盖材料的传热系数，W/(m2·K)；Kc为保温幕阻碍系数；Kn为内遮阳阻碍系数；Atop为顶窗玻璃的面积，m2；Ga为减反射玻璃的传热系数，W/(m2·K)；S1、S2、S3、S4分别是温室4个侧面的面积，m2，其中S1对应减反射玻璃，S2、S3、S4分别对应普通玻璃；Tog1、Tog2、Tog3分别是对应的普通玻璃紧邻的温室的室内温度，℃，To为室外空气温度，℃；温室与外界的通风换热能量Qvent表示如下：式中，Aw为开窗面积，m2；Cd为平均通风排气系数；ΔT为温室内外的温度差，K；Cw为与风速相关的通风率系数；U0为室外风速，m/s；Xw为通风换热修正系数；植物蒸腾作用所需要的能量Qtrans表示如下：Qtrans＝geL(χcrop‑χair)As               (12)式中，ge为蒸腾的电导系数，m/s；L为单位质量的叶片表面蒸发水所需要的能量，J/g；χcrop为作物水平高度的绝对湿度，g/m3；χair为温室室内空气的绝对湿度，g/m3；As为温室地表面积；温室中的作物冠层与空气热传导的能量Qcrop表示如下：式中，LAI为叶面积指数；rb表示边界层阻抗系数，s/m；Tleaf为叶表面的温度，℃；2)确定温室中需要辨识的温室因子，目标函数可表示为实际能量消耗和模拟能量消耗的均方根误差RMSE：上式中，X是物理模型参数向量灵敏度分析结果，l为时间序列的最大数，qreal为温室实际消耗能量，W；3)模型参数优化，过程如下：3.1)：按照模型参数，设置初始参数，包括群体大小M，收敛适应和最大代max_gen，然后所有粒子生成随机的位置和速度向量；3.2)：所有粒子按适应值f被分为两个部分，适应度小于等于平均适应度的粒子会被保存，然后会被PSO算法过程加强并且直接被传递到下一代的候选列表中，而那些适应度大于平均值的粒子将会通过遗传算法的交叉和变异算子被改进；3.3)：粒子通过群优化算法被增强，1≤j≤N，位置向量Xi和速度向量Vi第j维度更新为：vij(n+1)＝w(n)·vij(n)+c1·r1·(pij(n)‑xij(n))+c2·r2·(pbj(n)‑xij(n))  (27)xij(n+1)＝xij(n)+vij(n)                                     (28)上式中，vij是速度向量Vi的第j个维度，w是指惯性权重，xij是位置向量Xi的第j个维度，c1和c2是正的常数被称为加速率，r1和r2是在[0,1]范围内的随机变量，n是指当前这代的数字；在公式(27)中，新速度vij(n+1)由三部分组成，包括了惯性部分w(n)·vij(n)，个体部分c1·r1·(pij(n)‑xij(n))以及社交部分c2·r2·(pbj(n)‑xij(n))；惯性权重w在当前产生的数字表达为，w(n)＝wmax‑n·(wmax‑wmin)/max_n                              (29)上式中，max_n为最大迭代数，wmax是1.4，wmin是0.35；3.4)：假设位置向量Xi第j维，且xij∈[xij_min,xij_max]，这个位置向量Xi应该被编码而且在xij_min到xij_max区间相等的时间间隔内插入2R‑2个点，作为R的分辨率，两点之间的距离δ计算为：其中这个第j维的位置向量xij被转换为二进制，表达为：在IPSO‑GA算法中，在每一代中自适应交叉概率Pc和变异概率Pm是会改变的，如下所示Pc(n)＝Pc(1)×(1‑n/max_gen)                           (32)Pm(n)＝Pm(1)×(1‑n/max_gen)                           (33)上式中，Pc(n)和Pm(n)分别为第n代的交叉概率和变异概率，Pc(1)和Pm(1)分别为第一代的交叉概率和变异概率；3.5)：由于新粒子被遗传算法加强后便缺少了速度向量Vi和先前的最佳位置Pi，随机指定一个速度向量，然后把X″作为先前的最佳位置Pi，并作为PSO算法的算子；3.6)：比较先前每个粒子的适应度和每个粒子的最佳适应度，得到个体最佳适应度fi和总体最佳适应度Pbest；3.7)：如果n＜max_gen且RMSE＜min_rmse，min_rmse为设定的最小均方根偏差，就令n＝n+1返回3.3)；否则输出最佳值RMSE和优化参数向量X。