摘要
基于冷负荷时间序列固有的复杂性和不规则性,针对预测过程中容易出现梯度消失、模态混叠和过拟合等问题,提出一种集成变分模态分解(variational mode decomposition,VMD)和门控循环单元网络(gated recurrent unit,GRU)的VMD-GRU模型。对原始数据进行相关性分析,挑选出相关性高的进行预测;使用VMD将原始数据序列分解为独立固有模式函数;使用GRU对每个分量进行预测;将分量预测结果相加得出冷负荷预测值。为验证模型的有效性,以西安某大型公共建筑为例进行能耗分析,并与BP、 GRU、EMD-BP、VMD-BP、EMD-GRU等其他预测模型进行对比。实验结果表明,提出的VMD-GRU模型可有效解决梯度消失、模态混叠和过拟合等问题,预测精度显著提高,预测效果优于其它预测模型,符合大型公共建筑冷负荷的变化规律,为节能优化提供有力数据支撑。
全球能源消耗量伴随经济发展迅速增长,其中建筑能耗占比高达40
建筑能耗受多种因素影响,例如气候状况、室内条件和人员流
有学者将经验模式分解( empirical mode decomposition,EMD) 应用于预测领域并取得初步成
笔者结合VMD与GRU优势,提出VMD-GRU冷负荷预测模型,有效解决特征信息的提取,用于大型公共建筑冷负荷预测,并展开以下研究:
1)构建VMD-GRU冷负荷预测模型,使用VMD将原始数据序列分解为独立固有模式函数,用GRU对每个分量进行预测,将分量预测结果相加得出冷负荷预测值。
2)以西安某大型公共建筑相关数据为例进行实例分析,对模型的输入变量进行相关性分析,选取对冷负荷影响较大的输入变量。
3)采用VMD-GRU模型进行冷负荷预测,与其他预测模型进行比较,验证模型的有效性。
VMD模型依据序列数据的特点逐级进行平稳处理,GRU网络具有较强的非线性拟合能力, VMD-GRU模型处理冷负荷的非线性拟合及预测,降低复杂度的同时提升模型预测精度。
变分模态分解是一种新型非平稳信号自适应分解估计方法,目的是将原始复杂信号分解为K个调幅调频子信号。将大型公建冷负荷相关变量进行Pearson相关性分析,选取相关性高的数据重组为新的相关性序列。利用VMD进行分解及平稳化处理,分解为表示原始数据特征的多个分量。VMD分解的具体步骤如下
步骤一:初始化各模态、中心频率和算子。
步骤二:根据
, | (1) |
, | (2) |
其中:为原始信号;是信号的第k个分量;、、和分别代表、、和的傅里叶变换;n代表迭代次数。
步骤三:根据
。 | (3) |
步骤四:对于给定的判别精度e>0,若,则停止迭代,否则返回步骤二。
GRU网络是长短期记忆网络(long short-term memory,LSTM)的一种变体。通过使用“门”结构,极大避免梯度消失现象。GRU网络包括更新门、重置门和输出门3个部分。它将LSTM网络中的遗忘门和输入门合并成,原有重置门, LSTM门控网络的优点更新细胞状态和隐藏状态。
Step1:更新门和重置门。
更新门是过去时刻特征信息对现在特征信息的影响程度,阈值越大说明前时刻特征信息对现在影响越大,如
, | (4) |
。 | (5) |
Step2:候选隐藏状态
为候选隐藏状态,与ht-1、rt有关,当rt的取值接近于0,上一刻信息则被选择遗忘,见
。 | (6) |
Step3:隐藏状态
ht为当前隐藏状态,取决于ht-1和。如果zt趋近于0,表示上一时刻信息被遗忘; 如果 zt趋近于1,表示当前输入信息被遗忘,见
。 | (7) |
Step4:输出
。 | (8) |
过去的冷负荷状态会对当前状态产生长期影响,GRU通过控制与更新门限层控制模型的记忆能力,在不断迭代中,对历史数据的特征信息进行记忆与更新,历史数据会被赋予不同的权重值,已经训练的模型会对接下来的数据进行预测。

图1 预测模型结构
Fig. 1 Forecast model structure
分别是系统n个不同输入变量的观测值,是模型输出的观测值,GRU神经网络的输入向量。将历史数据依次通过多层GRU全连接层,并在接下来通过完全连接层合并特征,生成预测当前时刻后m个采样周期的输出。将GRU神经网络的预测值与真实值之间的平均平方差定义为损失函数,随时间反向传播。
选用Pearson分析法对冷负荷和输入变量间的相关性进行分析。公式为
, | (9) |
式中:rxy为2个变量的相关系数;xi、yi分别为2个变量的第i个数据点;分别为2个变量的均值;n为该变量中数据个数。
VMD-GRU模型的构建分为4部分:VMD分解、分解分量重构、GRU网络训练及最终预测结果输出。使用VMD对冷负荷序列和变量进行相关性分析,舍弃相关性低的变量,对相关性高的变量和冷负荷序列进行VMD分解,将分解后的分量进行重组。选取所有相关性高的变量第一个分量作为输入,冷负荷分解后的第一分量作为输出,在GRU网络中对预测模型进行训练,依次类推,对第二分量、第三分量等分别进行训练。对冷负荷预测数据相关性高的变量进行VMD分解,输入GRU网络中进行预测,得到冷负荷的预测分量,将分量预测值叠加得到预测值并输出结果,见
。 | (10) |

图2 方法流程图
Fig. 2 Flow chart of the method
实验数据来源于某大型公共建筑,建筑物高40.6 m,总建筑面积258×1
大型公共建筑冷负荷预测训练样本的输入层节点考虑到太阳辐射、室外温度导致建筑冷负荷存在滞后,预测模型中加了(T-1)h时刻室外空气温
实验以6月、7月60天中前50天每天14 h运行数据做为训练数据,采用7月份的25、26、27、28数据作为验证分析数据,使用Pearson相关性分析法对其相关性进行分析。
如

图3 变量间相关性热点图
Fig. 3 Hot spot map of Inter-variable correlation
相关性分析前 | 相关性分析后 | ||||
---|---|---|---|---|---|
变量 | 变量解释 | 变量 | 变量解释 | ||
相关变量
| O | T时刻冷负荷 |
输入向量 | U1 | T-1 h冷负荷 |
X1 | T-1 h时刻冷负荷 | U2 | T-2h冷负荷 | ||
X2 | T-2 h时刻冷负荷 | U3 | T h室外温度 | ||
X3 | T 时刻室外温度 | U4 | T-1 h室外温度 | ||
X4 | T-1 h时刻室外温度 | U5 | T时刻太阳辐射强度 | ||
X5 | T时刻太阳辐射强度 |
输出向量 | |||
X6 | T-1 h时刻太阳辐射强度 | O | T时刻冷负荷 | ||
X7 | T时刻相对湿度 | ||||
X8 | T时刻室外风速 |
经过相关性分析后,提取5个相关性较高的变量反映原始数据的主特征,进行VMD分解。参数设置,确定模态数量,按照实验选择模态数量对冷负荷序列进行VMD分解,获取VMD结果。
参数设置为:惩罚参数α=1 000;初始中心频率ω=0;收敛判据 r=1
模态个数 | 中心频率 | |||||
---|---|---|---|---|---|---|
2 | 6 | 1 452 | - | - | - | - |
3 | 2 | 675 | 2 548 | - | - | - |
4 | 1 | 662 | 1 444 | 3 278 | - | - |
5 | 1 | 641 | 1 381 | 2 420 | 3 750 | - |
6 | 1 | 640 | 1 368 | 2 288 | 2 336 | 3 827 |
对冷负荷序列进行VMD分解,分解结果如

图4 冷负荷VMD分解图
Fig. 4 The decomposed results of cooling load by VMD
基于VMD-GRU的大型公共建筑冷负荷预测,是将VMD分解的每一个分量分别采用GRU网络进行预测,最终获得所有冷负荷预测分量相加求和,得到冷负荷预测结果。为验证所提出模型的有效性,采用研究模型分别与BP、GRU、EMD-BP、VMD-BP、EMD-GRU模型进行实验对比分析。经过实验参数调试比较,GRU网络设置2个隐含层为最佳,激活函数为sigmoid,学习率为0.05,迭代次数为1 500,时间步的长度为28,最终设置GRU网络为5-3-7-1。
由于很多输入向量与预测结果相关性差、信息冗余度高等原因,导致多数预测模型在预测精度上未达到理想预测效
实验模型 | 训练迭代次数 | 相关性分析前 | 相关性分析后 | ||||
---|---|---|---|---|---|---|---|
训练时间/s | 测试时间/s | MAE | 训练时间/s | 测试时间/s | MAE | ||
BP模型 |
1 500 | 56.79 | 0.60 | 1 982.473 | 37.86 | 0.54 | 1 924.206 |
GRU模型 | 44.95 | 0.58 | 1 881.445 | 34.58 | 0.42 | 1 813.956 | |
EMD-BP模型 | 84.69 | 1.00 | 1 243.514 | 70.58 | 0.63 | 1 157.865 | |
VMD-BP模型 | 80.16 | 0.88 | 1 305.724 | 60.83 | 0.59 | 1 126.459 | |
EMD-GRU模型 | 77.26 | 0.79 | 1 149.968 | 60.79 | 0.57 | 1 082.471 | |
VMD-GRU模型 | 58.45 | 0.63 | 563.766 | 42.66 | 0.49 | 495.532 |

图5 冷负荷预测结果
Fig. 5 Prediction results of cooling load

图6 6种模型的预测结果对比
Fig. 6 Comparison of prediction results of six models

图7 6种模型的预测相对误差的对比
Fig. 7 Comparison of prediction relative errors of the six models
传统BP神经网络单独使用预测效果不能很好挖掘时间序列中的长时序关系。GRU引入门控单元可很好挖掘其中的长短期时序关
实验模型 | 拟合曲线 | 拟合 | ||
---|---|---|---|---|
BP模型 | -3.636% | 18.764 | 0.951 | |
GRU模型 | 5.196% | 19.593 | 0.871 | |
EMD-BP模型 | 3.457% | 15.371 | 0.926 | |
VMD-BP模型 | 3.654% | 16.176 | 0.942 | |
EMD-GRU模型 | 3.136% | 12.136 | 0.968 | |
VMD-GRU模型 | -3.022% | 9.306 | 0.98 5 |
在正态分布中值决定了曲线位置,越接近0,说明位置越靠近0,即大多数的预测相对误差值接近于0,模型预测效果越好。决定了曲线的尖峭程度,越小图形越尖峭,说明在该误差范围内的包含数值越多。通过对6个模型的相对误差分布图形对比可知,VMD-GRU模型,最接近于0,并且最小,说明该预测模型的相对误差在0附近数量值最大,模型最精确。


图8 实际冷负荷与模型预测冷负荷的回归拟合曲线
Fig. 8 Cooling load regression fitting curve between actual value and model prediction value
实验模型 | 拟合曲线 | 拟合 | 截距值 标准误差 | ||
---|---|---|---|---|---|
BP模型 | 4.917 | 0.928 | 1231.902 | 194.558 | |
GRU模型 | 2.891 | 0.961 | 587.624 | 149.179 | |
EMD-BP模型 | 2.014 | 0.971 | 765.991 | 124.519 | |
VMD-BP模型 | 1.601 | 0.977 | 649.514 | 111.027 | |
EMD-GRU模型 | 1.114 | 0.985 | 300.081 | 92.600 | |
VMD-GRU模型 | 1.045 | 0.992 | 56.098 | 68.774 |
针对大型公共建筑冷负荷预测,提出VMD-GRU预测模型,对输入输出数据间的相关性进行分析,利用VMD算法对冷负荷序列进行初步分解,利用GRU网络对分解序列分别进行预测,得出以下结论:
1)该大型公共建筑的 T 时刻输出冷负荷与T时刻室外温度、T-1时刻室外温度、T时刻太阳辐射量、T-1 时刻冷负荷、T-2 时刻冷负荷这5个输入量相关性较高。对变量进行相关性分析,避免人工经验选取输入变量的不足,缩短预测时间。
2)以MAE、MAPE和r作为3种预测模型评价标准,实验结果表明,BP、GRU、EMD-BP、VMD-BP、EMD-GRU和VMD-GRU模型的MAE分别为1 924.206、1 813.956、1 157.865、1 126.459、1082.47和495.532;MAPE分别为0.181 6、0.171 8、0.102 5、0.122 5、0.092 7和0.041 9;r分别为0.515 64、0.525 22、0.803 61、0.793 24、0.781 76和0.954 9。
3) 冷负荷时间序列往往是非平稳、非线性。VMD- GRU模型更容易掌握建筑冷负荷时间序列的特征,实现原始序列平稳化,提高建筑冷负荷预测的精度,更适合工程实际应用。
参考文献
Liu T, Xu C, Guo Y, Chen H, A novel deep reinforcement learning based methodology for short-term HVAC system energy[J].International Journal of Refrigeration,2019,107(3): 39-51. [百度学术]
王烨, 白玉洁, 闫若文. 不同通风模式下变温度送风的节能特性研究[J]. 重庆大学学报, 2020, 43(5): 40-49. [百度学术]
Wang Y, Bai Y J, Yan R W. Study on energy saving characteristics of variable temperature air supply under different ventilation modes[J]. Journal of Chongqing University, 2020, 43(5): 40-49.(in Chinese) [百度学术]
Dong B, Li Z X, Mahbobur Rahman S M, et al. A hybrid model approach for forecasting future residential electricity consumption[J]. Energy and Buildings, 2016, 117: 341-351. [百度学术]
Yildiz B, Bilbao J I, Sproul A B. A review and analysis of regression and machine learning models on commercial building electricity load forecasting[J]. Renewable and Sustainable Energy Reviews, 2017, 73: 1104-1122. [百度学术]
Fumo N. A review on the basics of building energy estimation[J]. Renewable and Sustainable Energy Reviews, 2014, 31: 53-60. [百度学术]
钱志. 基于改进型SVR的电网短期负荷预测[J]. 中国电力, 2016, 49(8): 54-58. [百度学术]
Qian Z. Short-term load forecasting of power grid based on improved SVR[J]. Electric Power, 2016, 49(8): 54-58.(in Chinese) [百度学术]
李军, 黄杰. 基于自组织映射神经网络的局部自回归方法在网络流量预测中的应用[J]. 信息与控制, 2016, 45(1): 120-128. [百度学术]
Li J, Huang J. Prediction of network traffic using local auto-regressive methods based on self-organizing map neural network[J]. Information and Control, 2016, 45(1): 120-128.(in Chinese) [百度学术]
杨祖莨, 丁洁, 刘晋峰. 一种新的结合卷积神经网络的隧道内停车检测方法[J]. 重庆大学学报, 2021, 44(6): 49-59. [百度学术]
Yang Z L, Ding J, Liu J F. A new tunnel vehicle stopping detection methodology combined with convolutional neural network[J]. Journal of Chongqing University, 2021, 44(6): 49-59.(in Chinese) [百度学术]
彭飞, 马煜, 张晓华, 等. 基于改进深度森林的短期电力负荷预测模型[J]. 重庆大学学报, 2022, 45(5): 1-8. [百度学术]
Peng F, Ma Y, Zhang X H, et al. Short term power load forecasting model based on improved deep forest[J]. Journal of Chongqing University, 2022, 45(5): 1-8.(in Chinese) [百度学术]
Lauster M, Teichmann J, Fuchs M, et al. Low order thermal network models for dynamic simulations of buildings on city district scale[J]. Building and Environment, 2014, 73: 223-231. [百度学术]
李军, 卓泽赢. 基于经验小波变换和多核学习的风电功率短期预测[J]. 信息与控制, 2018, 47(4): 437-447. [百度学术]
Li J, Zhuo Z Y. Wind power short-term forecasting based on empirical wavelet transform and multiple kernel learning[J]. Information and Control, 2018, 47(4): 437-447.(in Chinese) [百度学术]
罗勇, 郑金, 宁美凤. 基于相似日搜索的PSO-WNN组合模型在短期电力负荷预测中的应用[J]. 信息与控制, 2013, 42(3): 371-376. [百度学术]
Luo Y, Zheng J, Ning M F. Application of PSO-WNN combined model based on similar days searching to short-term electric power load forecasting[J]. Information and Control, 2013, 42(3): 371-376.(in Chinese) [百度学术]
杨照坤, 宋万清, 曹琨. 基于量子遗传算法的FARIMA模型电力负荷短期预测[J]. 传感器与微系统, 2019, 38(10): 143-145. [百度学术]
Yang Z K, Song W Q, Cao K. Short-term prediction of power load based on FARIMA model for quantum genetic algorithm[J]. Transducer and Microsystem Technologies, 2019, 38(10): 143-145.(in Chinese) [百度学术]
Wang L, Lv S X, Zeng Y R. Effective sparse adaboost method with ESN and FOA for industrial electricity consumption forecasting in China[J]. Energy, 2018, 155: 1013-1031. [百度学术]
吴云, 雷建文, 鲍丽山, 等. 基于改进灰色关联分析与蝙蝠优化神经网络的短期负荷预测[J]. 电力系统自动化, 2018, 42(20): 67-72. [百度学术]
Wu Y, Lei J W, Bao L S, et al. Short-term load forecasting based on improved grey relational analysis and neural network optimized by bat algorithm[J]. Automation of Electric Power Systems, 2018, 42(20): 67-72.(in Chinese) [百度学术]
李伟, 韩力. 组合灰色预测模型在电力负荷预测中的应用[J]. 重庆大学学报(自然科学版), 2004, 27(1): 36-39. [百度学术]
Li W, Han L. Application of combination grey model in power load forecasting[J]. Journal of Chongqing University (Natural Science Edition), 2004, 27(1): 36-39.(in Chinese) [百度学术]
Wang J J, Li Y N. Multi-step ahead wind speed prediction based on optimal feature extraction, long short term memory neural network and error correction strategy[J]. Applied Energy, 2018, 230: 429-443. [百度学术]
赵兵, 王增平, 纪维佳, 等. 基于注意力机制的CNN-GRU短期电力负荷预测方法[J]. 电网技术, 2019, 43(12): 4370-4376. [百度学术]
Zhao B, Wang Z P, Ji W J, et al. A short-term power load forecasting method based on attention mechanism of CNN-GRU[J]. Power System Technology, 2019, 43(12): 4370-4376.(in Chinese) [百度学术]
钟波, 赵渊, 周家启. 基于粗神经网络的电力系统可靠性评估方法[J]. 重庆大学学报(自然科学版), 2005, 28(7): 38-42. [百度学术]
Zhong B, Zhao Y, Zhou J Q. Methods of power systenms remiability assessment based on rough nerual network[J]. Journal of Chongqing University (Natural Science Edition), 2005, 28(7): 38-42.(in Chinese) [百度学术]
Paudel S, Elmitri M, Couturier S, et al. A relevant data selection method for energy consumption prediction of low energy building based on support vector machine[J]. Energy and Buildings, 2017, 138: 240-256. [百度学术]
Escrivá-Escrivá G, Álvarez-Bel C, Roldán-Blay C, et al. New artificial neural network prediction method for electrical consumption forecasting based on building end-uses[J]. Energy and Buildings, 2011, 43(11): 3112-3119. [百度学术]
Yau Y H, Hasbi S. A comprehensive case study of climate change impacts on the cooling load in an air-conditioned office building in Malaysia[J]. Energy Procedia, 2017, 143: 295-300. [百度学术]
王蓓蓓, 朱峰, 嵇文路, 等. 中央空调降负荷潜力建模及影响因素分析[J]. 电力系统自动化, 2016, 40(19): 44-52. [百度学术]
Wang B B, Zhu F, Ji W L, et al. Load cutting potential modeling of central air-conditioning and analysis on influencing factors[J]. Automation of Electric Power Systems, 2016, 40(19): 44-52.(in Chinese) [百度学术]
Zhang Y G, Chen B, Pan G F, et al. A novel hybrid model based on VMD-WT and PCA-BP-RBF neural network for short-term wind speed forecasting[J]. Energy Conversion and Management, 2019, 195: 180-197. [百度学术]
杨熊, 于军琪, 郭晨露, 等. 基于改进PSO-BP神经网络的冰蓄冷空调冷负荷动态预测模型[J]. 土木与环境工程学报(中英文), 2019, 41(1): 168-174. [百度学术]
Yang X, Yu J Q, Guo C L, et al. Dynamic load forecasting model of ice storage air conditioning based on improved PSO-BP neural network[J]. Journal of Civil and Environmental Engineering, 2019, 41(1): 168-174.(in Chinese) [百度学术]
Maier M, Brückmann A, Schleußner E, et al. Using critical flicker frequency in the evaluation of visual impairment in preeclamptic women[J]. European Journal of Obstetrics & Gynecology and Reproductive Biology, 2017, 211: 188-193. [百度学术]
Gao S A, Huang Y F, Zhang S, et al. Short-term runoff prediction with GRU and LSTM networks without requiring time step optimization during sample generation[J]. Journal of Hydrology, 2020, 589: 125188. [百度学术]
Zhang Y G, Zhao Y, Kong C H, et al. A new prediction method based on VMD-PRBF-ARMA-E model considering wind speed characteristic[J]. Energy Conversion and Management, 2020, 203: 112254. [百度学术]