发明名称 一种基于coif5小波实时分解的风电功率实时预测计算方法
摘要 一种基于coif5小波实时分解的风电功率实时预测计算方法,其特点是:包括以下步骤:数据采集及预处理,采集所需预测风电场的整场连续6天的实际历史功率值,并剔除不良数据;数据分解及预测,利用coif5小波进行三层塔式分解,对各分量采用时间序列法分别进行单步预测并求和,将得到的新序列利用coif5小波实时分解并滚动预测,得到整场未来4个小时的风电功率预测值;仿真计算,输入仿真输入量,采用上述方法进行未来4个小时的风电功率预测;误差分析,误差评价标准执行国家能源局于2011年《风电场功率预测预报管理暂行办法》中的准确率和合格率的要求,给出该计算方法的精确度。具有科学、合理、实用,精度更高,能够满足在线使用要求等优点。
申请公布号 CN103065049B 申请公布日期 2015.08.19
申请号 CN201210582482.7 申请日期 2012.12.28
申请人 东北电力大学 发明人 杨茂;穆钢;严干贵;熊昊;王东;宋薇
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 吉林市达利专利事务所 22102 代理人 陈传林
主权项 一种基于coif5小波实时分解的风电功率实时预测计算方法,其特征在于:它包括以下步骤:(1)数据采集及预处理采集所需预测风电场的整场连续6天的实际历史功率值,采样时间间隔15分钟,剔除功率数值小于零和大于开机容量的数据,得到所需的输入数据序列:X<sub>0</sub>={x(t‑n),x(t‑n+1),x(t‑n+2),…,x(t)},n=576;(2)数据分解及预测按步骤(1)采集的数据序列X<sub>0</sub>={x(t‑n),x(t‑n+1),x(t‑n+2),…,x(t)},将序列X<sub>0</sub>用小波基coif5函数进行三层小波塔式分解,得到数据集a<sub>3</sub>,d<sub>1</sub>,d<sub>2</sub>,d<sub>3</sub>,分别对应于第三层低频分量和三层高频分量,其中:x(t)为整个风电场t时刻的实际输出功率,a<sub>3</sub>(t)为t时刻整个风电场实际输出功率小波分解后的第三层低频分量a<sub>3</sub>的值,d<sub>i</sub>(t)为t时刻整个风电场实际输出功率小波分解后的高频分量d<sub>i</sub>的值;i=1,2,3,风电功率实时预测是一种超短期4个小时的多步预测,采用时间序列法中的ARMA模型对小波分解后的4个分量分别进行预测,将得到的4个预测值作和,得到t+1时刻整场的预测值;即<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><msub><mover><mi>&alpha;</mi><mo>^</mo></mover><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000732766060000011.GIF" wi="1318" he="130" /></maths>其中,<img file="FDA0000732766060000012.GIF" wi="221" he="92" />为低频分量a<sub>3</sub>在t+1时刻的预测值;<img file="FDA0000732766060000013.GIF" wi="225" he="102" />为高频分量d<sub>i</sub>在t+1时刻的预测值,i=1,2,3,进而得到数据序列<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mrow><msub><mi>X</mi><mn>1</mn></msub><mo>=</mo><mo>{</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>}</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000732766060000014.GIF" wi="985" he="82" /></maths>重复上述方法类推16步,得到<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>X</mi><mn>16</mn></msub><mo>=</mo><mo>{</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>n</mi><mo>+</mo><mn>16</mn><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>16</mn><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000732766060000015.GIF" wi="1047" he="87" /></maths>其中<img file="FDA0000732766060000016.GIF" wi="434" he="78" />为一次多步预测的预测值;(3)仿真计算输入仿真输入量:风电场实测开机状态每15分钟采样间隔下6天的整场风电功率数据,根据步骤(2)数据分解及预测过程,得到整场风电功率实时预测结果;(4)误差分析误差评价标准执行中国华人民共和国国家能源局2011年《风电场功率预测预报管理暂行办法》中规定的准确率和合格率的要求,设P<sub>Mk</sub>为k时段的实际平均功率,P<sub>Pk</sub>为k时段的预测平均功率,N为日考核总时段,Cap为风电场开机容量,那么,实时预测准确率定义为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>r</mi><mn>1</mn></msub><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msqrt><mfrac><mn>1</mn><mn>16</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>16</mn></munderover><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>P</mi><mi>Mk</mi></msub><mo>-</mo><msub><mi>P</mi><mi>Pk</mi></msub></mrow><mi>Cap</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>)</mo></mrow><mo>&times;</mo><mn>100</mn><mo>%</mo></mrow>]]></math><img file="FDA0000732766060000021.GIF" wi="889" he="164" /></maths>合格率定义为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>r</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>16</mn></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>16</mn></munderover><msub><mi>B</mi><mi>k</mi></msub><mo>&times;</mo><mn>100</mn><mo>%</mo></mrow>]]></math><img file="FDA0000732766060000022.GIF" wi="470" he="148" /></maths>其中若<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msub><mi>P</mi><mi>Mk</mi></msub><mo>-</mo><msub><mi>P</mi><mi>Pk</mi></msub></mrow><mi>Cap</mi></mfrac><mo>)</mo></mrow><mo>&times;</mo><mn>100</mn><mo>%</mo><mo>&GreaterEqual;</mo><mn>85</mn><mo>%</mo><mo>,</mo></mrow>]]></math><img file="FDA0000732766060000023.GIF" wi="619" he="142" /></maths>则B<sub>k</sub>=1;若<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msub><mi>P</mi><mi>Mk</mi></msub><mo>-</mo><msub><mi>P</mi><mi>Pk</mi></msub></mrow><mi>Cap</mi></mfrac><mo>)</mo></mrow><mo>&times;</mo><mn>100</mn><mo>%</mo><mo>&lt;</mo><mn>85</mn><mo>%</mo><mo>,</mo></mrow>]]></math><img file="FDA0000732766060000024.GIF" wi="596" he="141" /></maths>则B<sub>k</sub>=0给出风电场记录的实测开机状态15分钟采样间隔下6天的整场风电功率数据作为输入,把模型计算的预测功率与实测功率通过上述误差评价标准进行误差计算,给出该计算方法的预测准确率。
地址 132012 吉林省吉林市船营区长春路169号