干旱气象, 2021, 39(5): 864-870 doi: 10.11755/j.issn.1006-7639(2021)-05-0864

业务技术应用

基于EMD与BP神经网络的北京市采暖季天然气消耗量预测

闵晶晶,, 王华,, 董颜

北京市气象服务中心,北京 100097

Forecast of Natural Gas Consumption in Heating Season Based on EMD and BP Neural Network Methods in Beijing

MIN Jingjing,, WANG Hua,, DONG Yan

Beijing Meteorological Service Center, Beijing 100097, China

通讯作者: 王华(1975— ),女,河南邓州人,高级工程师,主要从事应用气象研究工作. E-mail: huawang0707@sina.com

收稿日期: 2021-06-2   修回日期: 2021-08-27  

基金资助: 国家重点研发计划“科技冬奥”重点专项项目“冬奥会气象条件预测保障关键技术”(2018YFF0300100)
第5课题(2018YFF0300105)

Received: 2021-06-2   Revised: 2021-08-27  

作者简介 About authors

闵晶晶(1984— ),女,湖北大冶人,高级工程师,博士,主要从事应用气象研究工作. E-mail: minjj06@163.com

摘要

基于2002—2018年北京市采暖季天然气消耗量和地面常规气象观测资料以及社会统计年度资料,采用经验模态分解、相关分析等方法分析北京地区采暖季天然气消耗量的年际变化特征及影响要素。在此基础上,利用BP神经网络方法构建采暖季天然气消耗量预测模型,并对模型进行评估检验。结果表明:(1)近17 a北京市采暖季天然气消耗量呈现持续增加趋势,经验模态分解方法能够较好地分离出天然气的社会消耗量和气象消耗量,分别反映了天然气消耗量的长期变化趋势和短期波动特征。(2)采暖季天然气的社会消耗量与GDP、集中供热面积和常住人口数量呈显著正相关;气象消耗量与气温和负积温呈显著负相关,而与降水量和持续低温日数呈显著正相关,当采暖季气温明显偏低或出现较强降雪、持续低温等天气过程时,天然气的气象消耗量将大幅增加。(3)北京市采暖季天然气消耗量EMD_BP预测模型具有较好的预测效果,平均相对误差为5.6%,能够准确预测天然气气象消耗量的峰谷变化。

关键词: 经验模态分解; BP神经网络; 采暖季; 天然气消耗量; 预测模型

Abstract

Based on natural gas consumption, ground conventional meteorological observation data in Beijing in heating season from 2002 to 2018, as well as yearly social statistical information, the inter-annual variation characteristics of natural gas consumption in heating season in Beijing and its impact factors were analyzed by using empirical mode decomposition (EMD) and correlation analysis methods. And on this basis the forecast model of natural gas consumption in heating season was established by using back propagation (BP) neural network method, further the model was tested and evaluated. The results are as follows: (1) The natural gas consumption increased persistently in heating season from 2002 to 2018 in Beijing, and it was decomposed better into social and meteorological consumptions by EMD, which reflected long-term variation trend and short-term fluctuation of natural gas consumption, respectively. (2) The social consumption of natural gas in heating season had significantly positive correlation with GDP, intensive heating supply area and resident population number in Beijing. The meteorological consumption had significantly negative correlation with air temperature and negative accumulative temperature, while it was significantly positive correlated with precipitation and persistent low-temperature days. In heating season, when the air temperature was obviously lower or the continuous low-temperature and strong snowfall processes appeared, the meteorological consumption of natural gas would increase sharply. (3) The forecast model of gas consumption in heating season based on EMD_BP method had better prediction effect in Beijing, the average relative error was 5.6%, especially the model could predict accurately the peak and valley change of meteorological consumption of gas, which could provide scientific reference to a certain extent for energy planning and regulating.

Keywords: empirical mode decomposition; BP neural network; heating season; natural gas consumption; forecast model

PDF (1796KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

闵晶晶, 王华, 董颜. 基于EMD与BP神经网络的北京市采暖季天然气消耗量预测[J]. 干旱气象, 2021, 39(5): 864-870 doi:10.11755/j.issn.1006-7639(2021)-05-0864

MIN Jingjing, WANG Hua, DONG Yan. Forecast of Natural Gas Consumption in Heating Season Based on EMD and BP Neural Network Methods in Beijing[J]. Journal of Arid Meteorology, 2021, 39(5): 864-870 doi:10.11755/j.issn.1006-7639(2021)-05-0864

引言

近年来,我国经济的迅速发展导致能源消费量快速增加,能源安全与环境问题日益显著。天然气作为优质高效的清洁能源极大改善了空气环境质量。伴随北京煤改气的全面推进,天然气消费量呈现爆发式增长[1,2,3,4]。合理预测天然气消耗量,不仅有利于城市管网的优化调度,保证管网供气量,而且在实现节能减排、国家政策分析与制定以及国民经济核算等方面具有深远意义。据统计,天然气作为北京市主要的供暖燃料,采暖季的消耗量占全年消耗总量的80%,具有明显的季节性特征。如2009—2010、2012—2013年冬季北京遭遇了罕见的大雪、低温天气,天然气出现供需不平衡情况,直接威胁到城市居民的冬季供暖。因此,灾害性天气影响下冬季采暖季天然气预测显得尤为重要。

天然气消耗量常用的预测方法主要分为统计方法、经济数学模型和智能算法三类[5]。统计学方法以时间序列分析为主,广泛应用于能源预测[6,7,8],但统计模型需要长期数据积累,且模型中影响因子的贡献率较为固定,无法反映随时间推移某些特征因子贡献率的衰减。经济数学模型如灰色系统理论,可量化不确定因素和小样本信息,常用于构建短期天然气消耗量预测模型[9,10,11]。近年来,为了解决非线性决策问题,陆续采用神经网络、遗传算法等智能算法构建天然气消耗量预测模型[12,13,14,15]。经验模态分解是处理非平稳信号的一种时频分析方法,可依据输入信号的自身特点,将原始信号分解成若干个本征模态函数与残差,该方法是对傅里叶分析和小波变换的重大突破[16]。BP神经网络方法在处理非线性问题上具有极大的灵活性和自适应性,可通过训练数据集合不断刺激系统来调整网络权值和阈值,最终达到网络误差平方和最小的标准[17]

从天然气消耗量的影响因素来看,一方面受国民经济、城市人口、能源政策等社会经济发展要素的影响,燃气负荷表现出长期变化趋势[18];另一方面受天气影响,燃气负荷呈现季节变化和波动特征[19,20,21]。然而,以往研究天然气消耗量的影响因素大多局限在温度、湿度、天气状况、降雨量、风等基本要素,且多以温度作为唯一气象参考因子[22,23],没有过多考虑持续低温或强降雪等灾害性天气过程造成天然气消耗量大幅增加的情况,同时也缺乏社会经济发展要素影响的考量。为此,本文以北京市为例,采用经验模态分解、BP神经网络和相关分析等方法,根据天然气消耗量的非线性和非平稳的特点,综合分析经济、人口、供热面积和气象等多种影响要素,构建北京市采暖季天然气消耗量预测模型,并进行评估检验,以期为北京能源规划和调控提供一定参考。

1 资料和方法

1.1 资料

天然气资料为北京市2002—2018年采暖季(11月至次年3月)天然气月消耗量,累加后得到采暖季天然气总消耗量,其中2002—2015年作为预测模型构建样本,2016—2018年作为模型评估样本。

北京市有20个国家级气象观测站,其中佛爷顶、汤河口、上甸子、斋堂、霞云岭5站海拔较高且地处远郊,人类活动较少,不具代表性,因此选用采暖季其他15站最高、最低气温及平均气温和降水量、风速等气象要素逐日观测资料,并用上述15站平均值代表北京市气候特征,同时还引入了表征冬季冷暖程度的负积温和低温日数2个参数[24]。其中,负积温是指日平均气温低于0 ℃的累积值;低温日数是指日最低气温小于某个阈值的累积天数,本文将北京供暖部门关注的-5 ℃和-8 ℃作为阈值。

社会因素资料选用了2002—2018年GDP、集中供热面积、常住人口。

1.2 经验模态分解方法

经验模态分解(empirical mode decomposition, EMD)作为数据驱动的自适应非线性时变信号分解方法,是以傅里叶分析和小波分析为基础,适用于对非线性、非平稳信号逐级进行平稳化处理。EMD算法的筛选过程是将复杂时序数据分解为有限个本征模函数(intrinsic mode function, IMF),分解得到的IMF分量包含了原数据在不同时间尺度上的波动信息[25]。具体处理过程如下:

(1)对任一处理信号x(t),确定其局部最大值、最小值,将信号数据x(t)与上、下包络线均值m1(t)的差值记为h1(t),其表达式如下:

h1(t)=x(t)-m1(t)

(2)重复以上过程,从原始信号中筛选出的第一阶段IMF分量h1(t)通常包含信号的最高频分量。将h1(t)从x(t)中分离出来,得到去除高频分量的差值信号r1(t),把r1(t)作为新信号重复以上筛分步骤,直到第n阶段的残余信号为单调函数不再筛分出IMF分量的rn(t),其表达式如下:

r1(t)=x(t)-h1(t)r2(t)=r1(t)-h2(t)   rn(t)=rn-1(t)-hn(t)

根据分解算法,x(t)可表示为n个IMF分量和1个残差的和,表达式如下:

x(t)=j=1nhj(t)+rn(t)

式中:rn(t)为残差,表示信号中的平均趋势;hj(t)为第j个IMF分量,j=1,2,…,n,分别表示信号从高到低不同频率的成分。

筛分过程的结束主要依据类似Cauchy收敛准则[26],标准偏差SD通常设置为0.2~0.3,表达式如下:

SD=1n1n|hj(t)-hj-1(t)|2|hj-1(t)|2dt

1.3 BP神经网络方法

采用BP(back propagation)神经网络方法构建预测模型。BP神经网络适用于一些高度复杂和非线性问题的建模,主要包含输入层、输出层和隐含层。一般采用梯度下降法作为权值修正,常利用Sigmoid函数计算各节点神经元的激励函数,并归一化到-1~1之间。设W(k+1)和W(k)为BP网络各层节点之间相邻训练数的连接权值[27],△W(k)为权值的修正值,k为训练次数,权值修正公式如下:

W(k+1)=W(k)+W(k)

设有N组训练样本和M组检测样本,TiYi(k)分别为第i组训练样本的目标值和网络实际输出值,而TjYj(k)分别为第j组检验样本的目标值和网络实际输出值,则训练误差函数Etrain(k)和检验误差函数Etest(k)可表示为:

Etrain(k)=1Ni=1N[Ti-Yi(k)]2
Etest(k)=1Mj=1M[Tj-Yj(k)]2

2 天然气消耗量特征及其影响要素

2.1 采暖季天然气消耗量年际变化特征

2002—2018年,北京地区采暖季天然气消耗量呈波动增加趋势,2008年以前增加缓慢,之后迅速增加[图1(a)];采暖季天然气消耗量的IMF1高频分量呈波动变化,无明显变化趋势,2008年以前变化较平稳,之后波动幅度明显增大,最高峰值和次高峰值分别出现在2012年和2009年,谷值出现在2013年和2016年[图1(b)]。与采暖季天然气消耗量的变化趋势相似,天然气消耗量的IMF2趋势项呈持续增加趋势,2002—2012年消耗量平稳增加,2014—2015年急剧增加[图1(c)],增量超过20%。

图1

图1   北京地区2002—2018年采暖季天然气消耗量(a)及分解量IMF1(b)和IMF2(c)的年际变化

Fig.1   The inter-annual changes of natural gas consumption (a) and its decomposed IMF1 (b) and IMF2 (c) components in Beijing in heating season from 2002 to 2018


2.2 影响要素及年际变化特征

天然气消耗量受社会因素和气象条件的影响,具有典型的非平稳特征,其中气象条件是造成天然气负荷短期波动的主要原因。利用EMD方法分析北京市2002—2018年采暖季天然气消耗量,得到不同频率的本征模态分量IMF1和趋势分量IMF2,并与各影响因子作相关分析(表1)。可以看出,IMF1与平均气温、最低气温、最高气温、负积温均呈显著负相关,而与持续低温日数和降水量呈显著正相关;IMF2与GDP、集中供热面积、常住人口数量呈显著正相关,而与气象因子相关性不大,表明IMF2分量综合反映了GDP、人口、城市化、能源结构等社会因素的影响。因此,设定高频分量IMF1为气象消耗量,趋势分量IMF2为社会消耗量。

表1   北京市2002—2018年采暖季天然气消耗量的IMF1和IMF2分量与各因子的相关系数

Tab.1  Correlation coefficients between IMF1, IMF2 components of natural gas consumption and each factor in Beijing in heating season from 2002 to 2018

要 素IMF1IMF2
GDP-0.130.99*
风速-0.010.58
集中供热面积-0.090.94*
降水量0.69*0.19
常住人口数量-0.090.92*
负积温-0.82*0.11
平均气温-0.84*0.18
-5 ℃以下低温日数0.79*0.03
最低气温-0.79*0.06
-8 ℃以下低温日数0.76*0.02
最高气温-0.80*0.31

注:*表示通过α=0.05的显著性检验。

新窗口打开| 下载CSV


2.2.1 气象消耗量年际变化特征

近17 a北京地区采暖季平均气温大多在0 ℃以上,2009年和2012年采暖季因出现多次强降雪过程导致平均气温较常年同期(0.7 ℃)明显偏低,分别为-1.3、-0.9 ℃[图2(a)],而降水量较常年同期(27.7 mm)异常偏多,分别为74.6、106 mm[图2(b)],致使低温持续时间较长,低温日数多于常年,负积温低于-400 ℃(常年平均为-295.6 ℃);2013、2014、2016年采暖季平均气温分别为2.4、2.3、2.0 ℃,均明显高于常年同期,而降水量相对较少,低于-8 ℃的低温日数在20 d左右,较常年同期(42 d)明显偏少,负积温仅在-150 ℃左右。因此,当采暖季气温明显偏低或出现较强降雪、持续低温等天气过程时,天然气的气象消耗量会大幅增加,而当采暖季气温明显偏高、负积温偏少时,天然气的气象消耗量明显减少。

图2

图2   北京地区2002—2018年采暖季平均气温(a)、降水量(b)年际变化

Fig.2   The inter-annual changes of average temperature (a) and precipitation (b) in Beijing in heating season from 2002 to 2018


2.2.2 社会消耗量年际变化特征

近17 a北京市GDP、集中供热面积持续平稳增长,而常住人口数量增长相对缓慢,特别是2015年以后人口略有减少(图3),其中2013—2015年GDP、集中供热面积及人口数量均线性平稳增长,未出现与天然气社会消耗量IMF2相对应的明显增长现象。2013年9月北京市政府发布了《北京市2013—2017年清洁空气行动计划》,这可能是导致2013—2015年采暖季天然气用量明显增加的主要原因。

图3

图3   北京地区2002—2018年采暖季GDP、集中供热面积及常住人口数量的年际变化

Fig.3   The inter-annual changes of GDP, intensive heating supply area and permanent population number in Beijing in heating season from 2002 to 2018


3 天然气消耗量预测模型与检验

3.1 预测模型的建立

基于EMD方法在处理非平稳信号和BP神经网络方法在处理非线性问题的优势,提出了基于EMD_BP方法的采暖季天然气消耗量预测模型(图4)。选取平均气温Tave、最低气温Tmin、最高气温Tmax、风速V、降水量R、负积温NAT、小于-8 ℃的低温日数LTD-8等作为IMF1预测模型的输入因子,GDP、集中供热面积、常住人口数量、时间(年份)等作为IMF2预测模型的输入因子,分别建立天然气气象消耗量IMF1和社会消耗量IMF2的预测模型,两分量预测值求和得到天然气总消耗量预测值。其中,2002—2015年采暖季天然气消耗量数据作为训练样本,2016—2018年数据作为检验样本。

图4

图4   基于EMD_BP方法的采暖季天然气消耗量预测模型构建流程

Fig.4   The building flow of prediction model of natural gas consumption in heating season based on EMD_BP method


3.2 预测结果检验

利用采暖季天然气消耗量EMD_BP预测模型,分别对2002—2015年建模样本进行拟合、2016—2018年检验样本进行预测。可以看出,EMD_BP预测模型对天然气消耗量IMF1、IMF2分量和总消耗量均具有很好的拟合与预测效果(图5),尤其是气象消耗量IMF1预测模型能够准确预测出陡增点(2009、2012、2015、2017年)和陡减点(2010、2013、2016年)。

图5

图5   北京地区2002—2018年采暖季天然气消耗量IMF1(a)、IMF2(b)分量和总消耗量(c)的拟合值、预测值与真实值对比

Fig.5   Comparison of fitted and forecasted IMF1 (a), IMF2 (b) components of natural gas consumption and total consumption (c) with actual value in Beijing in heating season from 2002 to 2018


通过拟合、预测与真实值的结果进行检验,发现EMD_BP预测模型对2002—2015年数据样本拟合的平均相对误差为4.5%,对2016—2018年检验样本预测的平均相对误差为5.6%。总体来看,采暖季天然气消耗量EMD_BP预测模型具有较好的预测效果,能够准确预测出天气变化带来的影响。

4 结论和讨论

本文引入对非线性非平稳信号进行平稳化处理的EMD方法,对北京市采暖季天然气消耗量的时间序列进行分解及特征分析,探讨天然气消耗量的社会和气象影响因子,采用BP神经网络方法构建预测模型并检验,主要结论如下:

(1)EMD方法能够较好地提取北京市采暖季天然气消耗量线性增加的社会消耗量和波动变化的气象消耗量。

(2)北京市采暖季天然气的社会消耗量与GDP、集中供热面积、常住人口数量呈显著正相关,但受政策影响,2014—2015年天然气的社会消耗量出现剧增,此后增量有所减缓。

(3)北京市采暖季天然气的气象消耗量与气温、负积温呈显著负相关,而与降水量、小于-8 ℃的低温日数呈显著正相关。当采暖季气温明显偏低或出现较强降雪、持续低温等天气过程时,天然气的气象消耗量会大幅增加,而在气温明显偏高、负积温偏少时,气象消耗量明显减少。

(4)采暖季天然气消耗量EMD_BP预测模型具有较好的预测效果,2002—2015年数据样本拟合的平均相对误差为4.5%,2016—2018年检验样本预测的平均相对误差为5.6%,能够准确预测因天气变化而引起的天然气气象消耗量的峰谷变化。

参考文献

胡奥林.

新版《天然气利用政策》解读

[J]. 天然气工业, 2013, 33(2):110-114.

[本文引用: 1]

张涛, 杨娇敏.

我国区域天然气消费量的趋势预测

[J]. 天然气工业, 2016, 36(9):135-140.

[本文引用: 1]

甄仟, 郭晓茜, 闫强.

基于行业视角的中国天然气消费因素分解

[J]. 中国矿业, 2018, 27(2):50-57.

[本文引用: 1]

孙相博, 王岳.

基于改进灰色GM(1,1)模型的天然气负荷预测

[J]. 辽宁石油化工大学学报, 2019, 39(3):52-57.

[本文引用: 1]

苏欣, 段康, 袁宗明, .

城市天然气负荷特点及其预测研究

[J]. 油气储运, 2007, 26(1):5-9.

[本文引用: 1]

刘蓬军, 赵媛.

西安市天然气冬季短期负荷预测研究

[J]. 煤气与热力, 2009, 29(10):37-42.

[本文引用: 1]

BOROOJENI K G, AMINI M H, BAHRAMI S, et al.

A novel multi-time-scale modeling for electric power demand forecasting: From short-term to medium-term horizon

[J]. Electric Power Systems Research, 2017, 142:58-73.

DOI:10.1016/j.epsr.2016.08.031      URL     [本文引用: 1]

段旭, 王志强, 马军, .

河北省供暖季天然气负荷与气温关系分析研究

[J]. 石油与天然气化工, 2019, 48(5):42-48.

[本文引用: 1]

王兵, 田丹, 刘书文, .

天然气消费量的灰色模型预测

[J]. 管道技术与设备, 2008(5):3-5.

[本文引用: 1]

SHAIKH F, JI Q, SHAIKH P H, et al.

Forecasting China’s natural gas demand based on optimized nonlinear grey models

[J]. Energy, 2017, 140:941-951.

DOI:10.1016/j.energy.2017.09.037      URL     [本文引用: 1]

叶志宏, 陈二龙.

基于残差灰色模型的天然气消费预测研究

[J]. 重庆理工大学学报(自然科学), 2018, 32(9):99-102.

[本文引用: 1]

代丽娴.

天然气负荷预测的建模与应用研究

[J]. 计算机仿真, 2011, 28(10):180-183.

[本文引用: 1]

YU F, XU X Z.

A short-term load forecasting model of natural gas based on optimized genetic algorithm and improved BP neural network

[J]. Applied Energy, 2014, 134:102-113.

DOI:10.1016/j.apenergy.2014.07.104      URL     [本文引用: 1]

宋超, 宋娟, 任军.

基于Elman神经网络的短期天然气负荷预测建模与仿真

[J]. 工业控制计算机, 2016, 29(4):87-88.

[本文引用: 1]

PANAPAKIDIS I P, DAGOUMAS A S.

Day-ahead natural gas demand forecasting based on the combination of wavelet transform and ANFIS/genetic algorithm/neural network model

[J]. Energy, 2017, 118:231-245.

DOI:10.1016/j.energy.2016.12.033      URL     [本文引用: 1]

徐晓刚, 徐冠雷, 王孝通, .

经验模式分解(EMD)及其应用

[J]. 电子学报, 2009, 37(3):581-585.

[本文引用: 1]

经验模式分解(Empirical Mode Decomposition,EMD)是一种数据驱动的自适应非线性时变信号分解方法,可以把数据分解成具有物理意义的少数几个模式函数分量.本文总结归纳了一维EMD、二维EMD方面的主要工作,比较了不同方法存在的优点与不足,指出了EMD研究存在的难题和瓶颈,并给出了EMD研究与应用的发展趋势.

马学款, 普布次仁, 唐叔乙, .

人工神经网络在西藏中短期温度预报中的应用

[J]. 高原气象, 2007, 26(3):491-495.

[本文引用: 1]

王耀庭, 孟春雷, 苗世光, .

基于人口和能耗数据估算北京人为热排放

[J]. 干旱气象, 2020, 38(1):89-99.

[本文引用: 1]

井元元, 李宏宇, 许启慧.

河北省采暖期低空大气逆温特征时空变化

[J]. 干旱气象, 2018, 36(4):624-635.

[本文引用: 1]

顾婷婷, 骆月珍, 潘娅英.

杭州市燃气负荷与气象条件的响应关系及其预测模型

[J]. 气象科技, 2014, 42(6):1154-1158.

[本文引用: 1]

舒漫, 刘夏兰, 徐婷, .

城市天然气短期日需求量预测新模型

[J]. 天然气工业, 2018, 38(6):128-132.

[本文引用: 1]

李谦益, 吴渊.

基于气温累积效应的冬季燃气日负荷预测——以西安市为例

[J]. 油气储运, 2016, 35(9):1014-1017.

[本文引用: 1]

何春蕾, 段言志, 邬宗婧, .

基于气温的城市燃气短期日负荷预测模型——以四川省成都地区为例

[J]. 天然气工业, 2013, 33(4):131-134.

[本文引用: 1]

高庆九, 张荣, 管兆勇.

近50 a来华北地区负积温变化特征

[J]. 大气科学学报, 2012, 35(4):448-457.

[本文引用: 1]

陈川, 陈冬林, 何李凯.

基于BPNN-EMD-LSTM组合模型的城市短期燃气负荷预测

[J]. 安全与环境工程, 2019, 26(1):149-155.

[本文引用: 1]

HUANG N E, SHEN Z, LONG S R, et al.

The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis

[J]. Proceedings of the Royal Society A: Mathematical Physical and Engineering Sciences, 1998, 454:903-995.

DOI:10.1098/rspa.1998.0193      URL     [本文引用: 1]

闵晶晶, 孙景荣, 刘还珠, .

一种改进的BP算法及在降水预报中的应用

[J]. 应用气象学报, 2010, 21(1):55-62.

[本文引用: 1]

/