发明名称 基于广义动态模糊神经网络的燃气日负荷组合预报方法
摘要 本发明提供的是一种基于广义动态模糊神经网络的燃气日负荷组合预报方法。(1)采集城市燃气历史记录数据作为历史时间序列数据;(2)对历史时间序列数据进行异常数据判断和处理;(3)利用广义回归神经网络,对历史负荷时间序列即样本数据进行差分处理再通过网络进行预测;(4)利用灰色神经网络,将输入历史负荷时间序列的一次累加生成数据作为网络的输入,输出对应预测日负荷的一次累加生成数据,训练网络,最后将输出值进行一次累减逆生成处理;(5)将步骤(3)、(4)得到的预测值作为广义动态模糊神经网络的输入,并对数据进行分组。本发明针对燃气日负荷的随机性,不稳定性,周期性等特点,采用了组合预报的方法,预测精度更高。
申请公布号 CN102073785B 申请公布日期 2013.02.27
申请号 CN201010561217.1 申请日期 2010.11.26
申请人 哈尔滨工程大学 发明人 陈虹丽;王辉;齐红芳;郑薇;李少阳;王岩
分类号 G06F19/00(2006.01)I;G06N3/02(2006.01)I 主分类号 G06F19/00(2006.01)I
代理机构 代理人
主权项 1.一种基于广义动态模糊神经网络的燃气日负荷组合预报方法,其特征在于具体包括以下步骤:(1)由城市燃气调度中心采集城市燃气历史记录数据作为历史时间序列数据;(2)对历史时间序列数据进行异常数据判断和处理,异常数据判断和处理包括:参考相似日数据,进行平滑、修正、归一化、组成历史负荷时间序列;所述相似日是指离其最近的具有相近天气、温度、星期几的日期;(3)利用广义回归神经网络,对历史负荷时间序列即样本数据进行差分处理再通过网络进行预测,根据样本数据逼近隐含的预测映射关系,得到预测日负荷值;具体为:广义回归神经网络的输出为:<maths num="0001"><![CDATA[<math><mrow><msup><mi>y</mi><mi>q</mi></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>r</mi><mi>i</mi><mi>q</mi></msubsup><mo>&times;</mo><mi>&omega;</mi><msub><mn>2</mn><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>而<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>r</mi><mi>i</mi><mi>q</mi></msubsup><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msup><mrow><mo>(</mo><msqrt><munder><mi>&Sigma;</mi><mi>j</mi></munder><mrow><mo>(</mo><msub><mrow><mi>&omega;</mi><mn>1</mn></mrow><mi>ji</mi></msub><mo>-</mo><msubsup><mi>x</mi><mi>j</mi><mi>q</mi></msubsup><msup><mo>)</mo><mn>2</mn></msup><mo>)</mo></mrow></msqrt><mo>&times;</mo><msub><mrow><mi>b</mi><mn>1</mn></mrow><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中ω2<sub>i</sub>,ω1<sub>ji</sub>为网络权系数,b1<sub>i</sub>为网络阈值,通过历史数据利用梯度法,通过改变步长,增加动量项自适应确定,y<sup>q</sup>为网络输出,<img file="FDA00002168683100013.GIF" wi="44" he="56" />为网络输入值,i=1,2,…,n,j=1,2,…,m;通过对时间序列的差分来去除非平稳性,具体步骤为:1)将预处理后的时间序列进行差分得到的一个新的信号,χ<sup>(1)</sup>(k)=χ<sup>(0)</sup>(k)-χ<sup>(0)</sup>(k-1)        (3)作为网络训练时间序列;2)建立和训练广义回归神经网络;3)进行预测时,网络输出值即为预测日负荷量的一阶差分,经还原即得到预测的日负荷量;(4)利用灰色神经网络,将输入历史负荷时间序列的一次累加生成数据作为网络的输入,输出对应预测日负荷的一次累加生成数据,训练网络,最后将输出值进行一次累减逆生成处理,得到预测日负荷值;具体步骤为:1)将预处理后的时间序列{χ<sup>(0)</sup>(1),χ<sup>(0)</sup>(2),…χ<sup>(0)</sup>(n)}进行一次累加生成处理,得到χ<sup>(1)</sup>={χ<sup>(1)</sup>(1),χ<sup>(1)</sup>(2),…χ<sup>(1)</sup>(n)},作为网络训练时间序列,其中;<maths num="0003"><![CDATA[<math><mrow><msup><mi>&chi;</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msup><mi>&chi;</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>2)建立和训练广义回归神经网络;3)进行预测时,将输出值<img file="FDA00002168683100022.GIF" wi="67" he="54" />做累减还原:<maths num="0004"><![CDATA[<math><mrow><msup><mover><mi>&chi;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>o</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mover><mi>&chi;</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msup><mover><mi>&chi;</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>得到预测日负荷值;(5)将步骤(3)、(4)得到的预测值作为广义动态模糊神经网络的输入,并对数据进行分组,最终得到预测日负荷值;具体方法为:1)初始化系统预定义参数ε<sub>min</sub>,ε<sub>max</sub>,k<sub>s</sub>,k<sub>err</sub>;2)当前第一个观察到来时产生第一个规则,确定第一个规则的参数,其高斯函数宽度根据ε-完备性确定,将当前输入训练样本作为新规则的高斯函数中心,其结论参数的确定采用最小二乘法;3)从第二组样本数据开始,每来一组样本,分别计算实际的输出误差e<sub>k</sub>和马氏距离md<sub>k</sub>(j),找到<maths num="0005"><![CDATA[<math><mrow><mi>J</mi><mo>=</mo><msub><mi>md</mi><mrow><mi>k</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mrow><mn>1</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>u</mi></mrow></munder><mrow><mo>(</mo><msub><mi>md</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>(ⅰ)如果md<sub>k,min</sub>&gt;k<sub>d</sub>且e<sub>k</sub>&gt;k<sub>e</sub>,则产生一个新规则;(a)规则参数的确定同步骤(5)中的第2)步;(b)计算出误差减小率ρ<sub>l</sub>和η<sub>l</sub>,l=1,2,…,u,若η<sub>l</sub>&lt;k<sub>err</sub>,则删除第l个规则,否则调整结论参数;η<sub>l</sub>:第l条规则的重要性;(ⅱ)如果md<sub>k,min</sub>&gt;k<sub>d</sub>且e<sub>k</sub>&lt;k<sub>e</sub>,只调整结论参数;(ⅲ)如果md<sub>k,min</sub>&lt;k<sub>d</sub>且e<sub>k</sub>&gt;k<sub>e</sub>,则调整已存在规则高斯函数宽度及调整结论参数;(ⅳ)如果md<sub>k,min</sub>&lt;k<sub>d</sub>且e<sub>k</sub>&lt;k<sub>e</sub>,调整结论参数。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室