发明名称 基于自适应增强算法的短期电力负荷在线预测方法
摘要 本发明公开了电力负荷预测技术领域中的一种基于自适应增强算法的短期电力负荷在线预测方法。包括选取影响气象数据的M个因素,并提取过去L天每个影响气象数据的因素的实测值,形成气象数据矩阵S<sub>L×M</sub>;提取过去L天中每天n个时刻的电力负荷数据,形成电力负荷数据矩阵D<sub>L×n</sub>;在影响气象数据的因素中,选取与电力负荷数据关联度最大的m个因素,将其作为有效成分,并根据过去L天有效成分的实测值,形成有效气象数据矩阵T<sub>L×m</sub>;利用有效气象数据矩阵T<sub>L×m</sub>和电力负荷数据矩阵D<sub>L×n</sub>求解短期电力负荷预测模型;利用短期电力负荷预测模型进行电力负荷预测。本发明可以有效地消除数据噪音对模型预测精度的影响,得到较为精确且稳定的预测结果。
申请公布号 CN103793887B 申请公布日期 2016.05.18
申请号 CN201410053462.X 申请日期 2014.02.17
申请人 华北电力大学 发明人 许刚;谈元鹏;马爽
分类号 G06T5/00(2006.01)I;G06N3/12(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 张文宝
主权项 一种基于自适应增强算法的短期电力负荷在线预测方法,其特征是所述方法包括:步骤1:选取影响气象数据的M个因素,并提取过去L天每个影响气象数据的因素的实测值,形成气象数据矩阵S<sub>L×M</sub>;其中,M和L为设定值;步骤2:提取过去L天中每天n个时刻的电力负荷数据,形成电力负荷数据矩阵D<sub>L×n</sub>;其中,n为设定值;步骤3:在影响气象数据的因素中,选取与电力负荷数据关联度最大的m个因素,将所述m个因素作为有效成分,根据过去L天有效成分的实测值,形成有效气象数据矩阵T<sub>L×m</sub>;其中,m为设定值;步骤4:利用有效气象数据矩阵T<sub>L×m</sub>和电力负荷数据矩阵D<sub>L×n</sub>求解短期电力负荷预测模型;步骤5:利用短期电力负荷预测模型进行电力负荷预测;所述步骤4包括:子步骤B1:以第k天有效成分的实测值T<sub>k</sub>、与第k天日期类型相同的前一天的有效成分实测值T<sub>k‑1</sub>、与第k天日期类型相同的前一天的电力负荷数据D<sub>k‑1,j</sub>和第k天的日期类型w作为输入样本,以第k天的电力负荷数据D<sub>k,j</sub>作为输出样本,建立样本集合{[T<sub>k‑1</sub>,T<sub>k</sub>,D<sub>k‑1,j</sub>,w],D<sub>k,j</sub>},记为{x<sub>p</sub>,t<sub>p</sub>};其中,k=2,3,...,L,j=1,2,...,n,p=1,2,...,Ln;子步骤B2:随机选取单隐层人工神经网络,设定权重步长a和收敛阈值ε,令初始迭代次数iter=1,权重分布<img file="FDA0000929660670000021.GIF" wi="231" he="142" />子步骤B3:在样本集合{x<sub>p</sub>,t<sub>p</sub>}的权重分布<img file="FDA0000929660670000022.GIF" wi="82" he="79" />下,训练所述单隐层人工神经网络,利用极端学习机ELM算法计算得到所述单隐层人工神经网络的参数<img file="FDA0000929660670000023.GIF" wi="83" he="74" />和<img file="FDA0000929660670000024.GIF" wi="117" he="92" />其中,参数<img file="FDA0000929660670000025.GIF" wi="81" he="69" />为单隐层人工神经网络中第l个神经元和输出节点之间的权重;参数<img file="FDA0000929660670000026.GIF" wi="76" he="77" />为单隐层人工神经网络中第l个神经元的阈值;子步骤B4:判断是否满足收敛条件<img file="FDA0000929660670000027.GIF" wi="485" he="238" />当<img file="FDA0000929660670000028.GIF" wi="454" he="215" />时,执行子步骤B6;否则,执行子步骤B5;其中,<maths num="0001"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>&sigma;</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></munderover><msub><mi>&alpha;</mi><mi>&sigma;</mi></msub><msub><mi>h</mi><mi>&sigma;</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000929660670000029.GIF" wi="469" he="151" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mi>&sigma;</mi></msub><mo>=</mo><mi>l</mi><mi>n</mi><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>&beta;</mi><mi>&sigma;</mi></msub></mfrac><mo>)</mo></mrow><mo>/</mo><munderover><mo>&Sigma;</mo><mrow><mi>&sigma;</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></munderover><mi>l</mi><mi>n</mi><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>&beta;</mi><mi>&sigma;</mi></msub></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00009296606700000210.GIF" wi="516" he="157" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mi>&sigma;</mi></msub><mo>=</mo><mfrac><mrow><mo>|</mo><msub><mi>h</mi><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>&sigma;</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>t</mi><mi>&sigma;</mi></msub><mo>|</mo><mo>-</mo><mi>e</mi><mi>r</mi><mi>r</mi></mrow><mrow><mi>e</mi><mi>r</mi><mi>r</mi></mrow></mfrac><mo>&times;</mo><mi>a</mi><mo>;</mo></mrow>]]></math><img file="FDA00009296606700000211.GIF" wi="582" he="155" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>h</mi><mi>&sigma;</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mover><mi>N</mi><mo>&OverBar;</mo></mover></munderover><msubsup><mi>&gamma;</mi><mi>l</mi><mi>&sigma;</mi></msubsup><mi>G</mi><mrow><mo>(</mo><msup><mrow><mo>(</mo><msubsup><mi>&omega;</mi><mi>l</mi><mi>&sigma;</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><msub><mi>x</mi><mi>p</mi></msub><mo>+</mo><msubsup><mi>b</mi><mi>l</mi><mi>&sigma;</mi></msubsup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00009296606700000212.GIF" wi="678" he="159" /></maths>G(·)为单隐层人工神经网络的激活函数且<img file="FDA00009296606700000213.GIF" wi="319" he="142" /><maths num="0005"><math><![CDATA[<mrow><mi>e</mi><mi>r</mi><mi>r</mi><mo>=</mo><msqrt><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>L</mi><mi>n</mi></mrow></munderover><mrow><mo>(</mo><msubsup><mi>d</mi><mi>p</mi><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>&times;</mo><msup><mrow><mo>&lsqb;</mo><msub><mi>h</mi><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>t</mi><mi>p</mi></msub><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>L</mi><mi>n</mi></mrow></munderover><msubsup><mi>d</mi><mi>p</mi><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup></mrow></mfrac></msqrt><mo>;</mo></mrow>]]></math><img file="FDA00009296606700000214.GIF" wi="700" he="311" /></maths><img file="FDA00009296606700000215.GIF" wi="61" he="61" />为单隐层人工神经网络神经元的个数;<img file="FDA0000929660670000031.GIF" wi="70" he="71" />为第l个输入节点和神经元之间的权重;子步骤B5:令<img file="FDA0000929660670000032.GIF" wi="474" he="103" />iter=iter+1,转至子步骤B3;其中,<maths num="0006"><math><![CDATA[<mrow><msub><mi>e</mi><mi>p</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>0</mn><mo>,</mo><mo>|</mo><msub><mi>h</mi><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>)</mo><mo>-</mo><msub><mi>t</mi><mi>p</mi></msub><mo>|</mo><mo>&gt;</mo><mi>e</mi><mi>r</mi><mi>r</mi></mtd></mtr><mtr><mtd><mn>1</mn><mo>,</mo><mo>|</mo><msub><mi>h</mi><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>)</mo><mo>-</mo><msub><mi>t</mi><mi>p</mi></msub><mo>|</mo><mo>&le;</mo><mi>e</mi><mi>r</mi><mi>r</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000929660670000033.GIF" wi="574" he="207" /></maths><maths num="0007"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mi>p</mi></msub><mo>=</mo><mfrac><mrow><mo>|</mo><msub><mi>h</mi><mrow><mi>i</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>t</mi><mi>p</mi></msub><mo>|</mo><mo>-</mo><mi>e</mi><mi>r</mi><mi>r</mi></mrow><mrow><mi>e</mi><mi>r</mi><mi>r</mi></mrow></mfrac><mo>&times;</mo><mi>a</mi><mo>;</mo></mrow>]]></math><img file="FDA0000929660670000034.GIF" wi="581" he="173" /></maths>子步骤B6:H(x<sub>p</sub>)即为所求的短期电力负荷预测模型。
地址 102206 北京市昌平区回龙观朱辛庄2号
您可能感兴趣的专利