发明名称 一种径流式小水电集群发电功率短期预测方法及预测系统
摘要 一种径流式小水电集群发电功率短期预测方法及预测系统,属于水力发电功率预测技术领域。该方法将接入同一变电站的若干个小水电视为一个集群,分别对集群内所有单个小水电功率数据和集群整体功率数据进行预测,然后将两方面预测结果相融合获得最终的集群发电预测数值。本发明还提供一种径流式小水电集群发电功率短期预测系统。本发明将分散的多个径流式小水电归为一个集群,充分考虑径流式小水电的累积效应和滞后效应,计算未来时刻功率预测数值时将单点功率变化率与平均功率变化率相结合,集群功率预测结果融合了整体功率预测结果和所有单个小水电功率预测结果,针对径流式小水电功率输出规律性差而无法准确预测技术问题,本发明预测准确性较好。
申请公布号 CN103679288A 申请公布日期 2014.03.26
申请号 CN201310648773.6 申请日期 2013.12.04
申请人 华北电力大学 发明人 彭文;刘文霞;辜庭帅;赵天阳;李鹤
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京麟保德和知识产权代理事务所(普通合伙) 11428 代理人 韩建功
主权项 1.一种径流式小水电集群发电功率短期预测方法,其特征在于,该方法,(I)将两个长度为n的数据序列X={X<sub>1</sub>,...,X<sub>k</sub>,...,X<sub>n</sub>}和Y={Y<sub>1</sub>,...,Y<sub>k</sub>,...,Y<sub>n</sub>}的归一化形系数S定义为:<![CDATA[<math><mrow><mi>S</mi><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>w</mi><mi>k</mi></msub><mo>|</mo><mover><msub><mi>X</mi><mi>k</mi></msub><mo>&OverBar;</mo></mover><mo>-</mo><mover><msub><mi>Y</mi><mi>k</mi></msub><mo>&OverBar;</mo></mover><mo>|</mo><mo>,</mo></mrow></math>]]></maths>   公式(1);其中,<![CDATA[<math><mrow><mover><msub><mi>X</mi><mi>k</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mrow><mo>(</mo><msub><mi>X</mi><mi>max</mi></msub><mo>-</mo><msub><mi>X</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><msub><mi>X</mi><mi>max</mi></msub><mo>-</mo><msub><mi>X</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><![CDATA[<math><mrow><mover><msub><mi>Y</mi><mi>k</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mrow><mo>(</mo><msub><mi>Y</mi><mi>max</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><msub><mi>Y</mi><mi>max</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>w<sub>k</sub>=T<sup>(n-k)</sup>(0<T<1),<img file="FDA0000430202310000014.GIF" wi="166" he="75" />分别是X和Y归一化后的序列,若X<sub>max</sub>=X<sub>min</sub>或者Y<sub>max</sub>=Y<sub>min</sub>,则<img file="FDA0000430202310000015.GIF" wi="165" he="74" />为全0序列;w<sub>k</sub>为日期权重,其中的T<sup>(n-k)</sup>随着k值增加而逐渐变小,反映了功率预测中“近大远小”的原则;S能够很好地反映了X和Y的形状差异,如X和Y形状完全一致,即使在垂直方向上有位移偏差,S始终为0;如X和Y形状存在差异,则S不为0且差异越大S越大;(II)基于相似日区间的功率短期预测方法该方法既适用于单个小水电功率预测又适用于集群整体功率预测,只需替换相应的预测数据,因此在方法描述过程中不对数据进行区分;设P={P<sub>k</sub>}(k=1,...,Pn)为预测日前n天的功率序列,其中Pn=n*DC,DC为每天采集的功率数值个数,Q={Q<sub>j</sub>}(j=1,...,PN)为历史功率数据,其中PN=N*DC,N为历史样本天数;功率预测方法就是在序列Q中检索出与序列P形状最相近的若干相似日区间,将各相似日区间之后时间点功率变化趋势进行累加,并扩展到序列P之后以获得预测日的功率序列;方法具体如下:(1)相似日区间选取功率数据采集的频率非常密集,如果直接在这一层面寻找相似日区间,则检索过程将非常耗时;本发明提出一种基于日电量数据和功率数据两层搜索的相似日区间选取算法;首先,通过功率计算每天的日电量,在日电量数据层面上,寻找相似度较高的若干个候选相似日区间;然后在功率数据层面上再次对这个候选相似日区间进行过滤,最终得到最佳相似日区间;(1.1)基于日电量数据的相似日区间选取方法设<img file="FDA0000430202310000021.GIF" wi="452" he="78" />为预测日前n天的日电量序列,其中P<sub>i</sub><sup>E</sup>的计算公式为<![CDATA[<math><mrow><msup><msub><mi>P</mi><mi>i</mi></msub><mi>E</mi></msup><mo>=</mo><mfrac><mn>24</mn><mi>DC</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mi>DC</mi><mo>+</mo><mn>1</mn></mrow><mrow><mi>i</mi><mo>*</mo><mi>DC</mi></mrow></munderover><msub><mi>P</mi><mi>k</mi></msub><mo>,</mo><msup><mi>Q</mi><mi>E</mi></msup><mo>=</mo><mo>{</mo><msubsup><mi>Q</mi><mi>j</mi><mi>E</mi></msubsup><mo>}</mo><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>)</mo></mrow></mrow></math>]]></maths>为日电量历史序列,其中<img file="FDA0000430202310000024.GIF" wi="73" he="83" />的计算公式为<img file="FDA0000430202310000025.GIF" wi="490" he="149" />基于日电量数据的相似日区间选取算法是在序列Q<sup>E</sup>中找出长度为n且与序列P<sup>E</sup>形状最为相似的m个子序列,方法步骤为:(a)设最相似序列集合<img file="FDA00004302023100000212.GIF" wi="163" he="57" />C<sub>worst</sub>为C中与P<sup>E</sup>相似度最差的序列,初始化C<sub>worst</sub>=NULL,S<sub>worst</sub>为相应的归一化形系数,初始化为,S<sub>worst</sub>=-1,h=1;(b)采用公式(1)计算序列P<sup>E</sup>和<img file="FDA0000430202310000026.GIF" wi="529" he="79" />的归一化形系数<img file="FDA0000430202310000027.GIF" wi="118" he="83" />(c)如果|C|<m,将<img file="FDA0000430202310000028.GIF" wi="71" he="79" />添加到C中并转到步骤(e);(d)将<img file="FDA0000430202310000029.GIF" wi="70" he="79" />与序列C<sub>worst</sub>进行比较,如果<img file="FDA00004302023100000210.GIF" wi="92" he="83" />小于其归一化形系数S<sub>worst</sub>,用<img file="FDA00004302023100000211.GIF" wi="60" he="79" />替换C<sub>worst</sub>;(e)对C中所有的序列按归一化形系数从小到大排序,并更新C<sub>worst</sub>和S<sub>worst</sub>;(f)h=h+1,如果h&lt;N-n,转向步骤(b),否则退出;(1.2)基于功率数据的相似日区间选取方法根据集合C中的m个日电量序列中对应的日期序列,获得m个长度为DC*n的功率序列,然后筛选出其中与功率序列P最相似的M个功率序列(M&lt;m),仍然采用归一化形系数作为两个序列形状相似程度的评价函数,本方法步骤如下:(a)根据序列Q<sup>E</sup>,根据集合C中的所有日电量序列对应的日期序列,得到m个长度为DC*n的功率序列<img file="FDA0000430202310000031.GIF" wi="371" he="72" />(l=1,...,m),其中<img file="FDA0000430202310000032.GIF" wi="126" he="72" />是序列Q<sup>E</sup>中的子序列;(b)采用公式(1)计算P和<img file="FDA0000430202310000033.GIF" wi="122" he="70" />的归一化形系数S<sub>l</sub>;(c)将C<sup>power</sup>中的序列按照S<sub>l</sub>从小到大排序,并选取前M个作为最终的相似日区间序列集合C<sup>final</sup>;(2)基于相似日区间的功率预测计算方法(2.1)相似日区间的权重计算集合C<sup>fin</sup>中各相似日区间的功率序列与序列P的归一化形系数为S<sub>i</sub>(i=1,2,...,M),其权重为:<![CDATA[<math><mrow><msub><mi>W</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>/</mo><msub><mi>S</mi><mi>i</mi></msub><mo>/</mo><mover><mi>S</mi><mo>&OverBar;</mo></mover><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>M</mi><mo>)</mo></mrow><mtext>,</mtext></mrow></math>]]></maths>其中,<![CDATA[<math><mrow><mover><mi>S</mi><mo>&OverBar;</mo></mover><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mn>1</mn><mo>/</mo><msub><mi>S</mi><mi>i</mi></msub><mo>;</mo></mrow></math>]]></maths>(2.2)设集合C<sup>final</sup>中的元素为功率序列<![CDATA[<math><mrow><msubsup><mi>C</mi><mi>i</mi><mi>final</mi></msubsup><mo>=</mo><mo>{</mo><msub><mi>Q</mi><mi>ik</mi></msub><mo>}</mo><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>M</mi><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>Pn</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>将序列<img file="FDA0000430202310000036.GIF" wi="130" he="78" />在时间轴上之后时刻的功率序列定义为<img file="FDA0000430202310000037.GIF" wi="214" he="79" />(i=1,...,M,t=1,...,TC),即<img file="FDA0000430202310000038.GIF" wi="111" he="76" />与<img file="FDA0000430202310000039.GIF" wi="56" he="76" />在时间上是连续的,其中TC是预测的数据个数;计算预测日各时间点的功率数值R<sub>j</sub>(j=1,2,..,TC),公式如下:R<sub>j</sub>=αE<sub>j1</sub>+βE<sub>j2</sub>,   公式(2);<![CDATA[<math><mrow><msub><mi>E</mi><mrow><mi>j</mi><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>P</mi><mi>Pn</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><mo>(</mo><msub><mi>W</mi><mi>s</mi></msub><mrow><mo>(</mo><msubsup><mi>Q</mi><mi>sj</mi><mo>*</mo></msubsup><mo>-</mo><msub><mi>Q</mi><mi>sPn</mi></msub><mo>)</mo></mrow><mo>/</mo><msub><mi>Q</mi><mi>sPn</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths><![CDATA[<math><mrow><msub><mi>E</mi><mrow><mi>j</mi><mn>2</mn></mrow></msub><mo>=</mo><mrow><mo>(</mo><mi>Pn</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mn>1</mn><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><mo>(</mo><msub><mi>W</mi><mi>s</mi></msub><mrow><mo>(</mo><mover><msubsup><mi>Q</mi><mi>sj</mi><mo>*</mo></msubsup><mo>&OverBar;</mo></mover><mo>-</mo><mover><msub><mi>Q</mi><mi>s</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>/</mo><mover><msub><mi>Q</mi><mi>s</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>Pn</mi><mover><mi>P</mi><mo>&OverBar;</mo></mover></mrow></math>]]></maths><![CDATA[<math><mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>Pn</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>u</mi><mo>=</mo><mn>1</mn></mrow><mi>Pn</mi></munderover><msub><mi>P</mi><mi>u</mi></msub></mrow></math>]]></maths><![CDATA[<math><mrow><mover><msub><mi>Q</mi><mi>s</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>Pn</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>u</mi><mo>=</mo><mn>1</mn></mrow><mi>Pn</mi></munderover><msub><mi>Q</mi><mi>su</mi></msub></mrow></math>]]></maths><![CDATA[<math><mrow><mover><msubsup><mi>Q</mi><mi>sj</mi><mo>*</mo></msubsup><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mrow><mi>Pn</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>u</mi><mo>=</mo><mn>1</mn></mrow><mi>Pn</mi></munderover><msub><mi>Q</mi><mi>su</mi></msub><mo>+</mo><msubsup><mi>Q</mi><mi>sj</mi><mo>*</mo></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>其中,E<sub>j1</sub>表示在j时间点各相似日区间功率变化率加权结果对预测值的贡献,E<sub>j2</sub>表示在j时间点各相似日区间平均功率变化率加权结果对预测值的贡献,α和β是调节因子,可通过训练过程确定数值大小;对于一个有M<sub>D</sub>个小水电站的集群D={D<sub>i</sub>}(i=1,2,...,M<sub>D</sub>),本方法步骤为:步骤1:对所有小水电应用上述(II)中的基于相似日区间的功率预测方法,分别获得其功率预测数据T<sub>ij</sub>(i=1,2,...,M<sub>D</sub>,j=1,2,...,TC);步骤2:对小水电集群整体数据应用上述(II)中的基于相似日区间的预测模型,获得集群整体功率预测数据T<sub>j</sub>(j=1,2,...,TC);步骤3:将基于两种数据的预测结果融合,获得集群最终的功率预测数值<img file="FDA0000430202310000044.GIF" wi="83" he="85" />公式如下:<img file="FDA0000430202310000045.GIF" wi="540" he="175" />公式(3);其中,<img file="FDA0000430202310000046.GIF" wi="145" he="80" />是权重因子,可通过训练过程确定数值大小。
地址 102206 北京市昌平区回龙观朱辛庄2号