发明名称 基于数据特征提取的风功率预测误差估计方法
摘要 本发明公开了一种基于数据特征提取的风功率预测误差估计方法,包括以下步骤:(1)统计分析近一年风电运行历史数据,分析计算风功率幅值、日前预测出力的波动程度、近3天风电出力波动程度和预测精度对风功率预测误差的影响权重系数;(2)根据(1)得到的权重系数,利用日前预测出力和近3天风电运行数据,计算日前风功率预测幅值和波动程度、近3天风电出力波动程度和预测精度,然后估计风功率预测误差。本发明统计分析并提取风电运行历史数据特征,利用近3天风电运行数据,估计风功率预测误差,该方法具有在线计算强度低、数据来源可靠且容易获得的特点,具有很高的工程实用价值。
申请公布号 CN103473461B 申请公布日期 2016.06.15
申请号 CN201310422579.6 申请日期 2013.09.16
申请人 东南大学;国电南瑞科技股份有限公司 发明人 张凯锋;丁恰;杨国强;王颖;陈汉一
分类号 G06F19/00(2011.01)I;H02J3/38(2006.01)I 主分类号 G06F19/00(2011.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 一种基于数据特征提取的风功率预测误差估计方法,其特征在于,具体包括以下步骤:1)统计分析近一年风电运行历史数据,计算风功率幅值、日前预测出力的波动程度、近3天风电出力波动程度和预测精度对风功率预测误差的影响权重系数;统计分析近一年风电运行历史数据,计算风功率幅值对风功率预测误差的影响权重系数α的公式如(1)所示,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>j</mi></mrow></msub></mfrac></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&alpha;</mi><mi>min</mi></msub><mo>=</mo><mover><msub><mi>K</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>/</mo><mover><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover></mrow></mtd><mtd><mrow><mo>&ForAll;</mo><msub><mi>K</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&alpha;</mi><mi>max</mi></msub><mo>=</mo><mover><msub><mi>K</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>/</mo><mover><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover></mrow></mtd><mtd><mrow><mo>&ForAll;</mo><msub><mi>K</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>&GreaterEqual;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000939169750000011.GIF" wi="1461" he="478" /></maths>其中,j表示历史数据的粒度,每15分钟为一个时间点,P<sub>act,j</sub>表示第j时间点风电的实际出力,P<sub>fore,j</sub>表示第j时间点风电的预测出力,<img file="FDA0000939169750000012.GIF" wi="85" he="85" />表示所有符合约束条件的K<sub>1,j</sub>的平均值,<img file="FDA0000939169750000013.GIF" wi="123" he="87" />表示所有符合约束条件K<sub>1,j</sub>对应的预测出力P<sub>fore,j</sub>的平均值,α<sub>min</sub>表示风功率幅值对风功率预测误差的影响权重系数的下限值,α<sub>max</sub>表示风功率幅值对风功率预测误差的影响权重系数的上限值;统计分析近一年风电运行历史数据,按出力相近的原则对日前预测出力数据进行分段处理,计算日前预测出力的波动程度对风功率预测误差的影响权重系数β的公式如(2)所示,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mn>2</mn><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>)</mo><mo>|</mo><msub><mi>t</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>|</mo></mrow><mrow><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mi>&beta;</mi><mo>=</mo><mover><msub><mi>K</mi><mrow><mn>2</mn><mo>,</mo><mi>i</mi></mrow></msub><mo>&OverBar;</mo></mover></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000939169750000014.GIF" wi="1318" he="270" /></maths>其中,i表示对日前预测出力数据所分的段数,P<sub>act,i,max</sub>表示第i个所分的时间段中风电实际出力的最大值;P<sub>act,i,min</sub>表示第i个所分的时间段中风电实际出力的最小值;P<sub>fore,i,max</sub>表示第i个所分的时间段中风电预测出力的最大值,其对应的时间为t<sub>fore,i,max</sub>;P<sub>fore,i,min</sub>表示第i个所分的时间段中风电预测出力的最小值,其对应的时间为t<sub>fore,i,min</sub>,<img file="FDA0000939169750000015.GIF" wi="90" he="86" />表示所有K<sub>2,i</sub>的平均值;统计分析近一年风电运行历史数据,计算近3天风电出力波动程度和预测精度对风功率预测误差的影响权重系数K<sub>3</sub>和K<sub>4</sub>的公式如(3)所示,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mn>3</mn></msub><mo>=</mo><mn>1</mn><mo>/</mo><mover><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>K</mi><mn>4</mn></msub><mo>=</mo><mover><msub><mi>K</mi><mrow><mn>4</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>=</mo><mover><mrow><mo>(</mo><mfrac><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>j</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>&OverBar;</mo></mover></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000939169750000021.GIF" wi="1305" he="318" /></maths>其中,j表示历史数据的粒度,每15分钟为一个时间点,P<sub>act,j</sub>表示第j时间点的风电实际出力,P<sub>fore,j</sub>为第j时间点的风电预测出力,<img file="FDA0000939169750000022.GIF" wi="191" he="176" />表示<img file="FDA0000939169750000023.GIF" wi="130" he="143" />的平均值,<img file="FDA0000939169750000024.GIF" wi="90" he="84" />表示P<sub>act,j</sub>的平均值;2)根据步骤1)得到的权重系数,利用日前风功率预测出力和近3天风电数据,计算日前风功率预测幅值和波动程度、近3天风电出力波动程度和近3天预测精度,然后估计风功率预测误差;分析风电日前预测出力数据,按出力相近的原则对日前预测出力数据进行分段处理,计算第i个所分的时间段数据的波动程度σ<sub>i</sub>的公式如(4)所示:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub></mrow><mrow><mo>|</mo><msub><mi>t</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000939169750000025.GIF" wi="1311" he="158" /></maths>其中,i表示对日前预测出力数据所分的段数,P<sub>fore,i,max</sub>表示第i个所分的时间段中风电预测出力的最大值,其对应的时间为t<sub>fore,i,max</sub>;P<sub>fore,i,min</sub>表示第i个所分的时间段中风电预测出力的最小值,其对应的时间为t<sub>fore,i,min</sub>;分析风电运行近3天数据,计算近3天风电出力波动程度S和近3天预测精度ε的公式如(5)所示<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>S</mi><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><mover><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mtd></mtr><mtr><mtd><mrow><mi>&epsiv;</mi><mo>=</mo><mover><msub><mi>&epsiv;</mi><mi>j</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mover><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>|</mo></mrow><msub><mi>P</mi><mrow><mi>a</mi><mi>c</mi><mi>t</mi><mo>,</mo><mi>j</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>&OverBar;</mo></mover></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000939169750000026.GIF" wi="1326" he="365" /></maths>其中,j表示历史数据的粒度,每15分钟为一个时间点,P<sub>act,j</sub>和P<sub>fore,j</sub>为第j时间点风电的实际出力和预测出力,<img file="FDA0000939169750000031.GIF" wi="334" he="182" />表示<img file="FDA0000939169750000032.GIF" wi="271" he="165" />的平均值,<img file="FDA0000939169750000033.GIF" wi="96" he="87" />表示P<sub>act,j</sub>的平均值;估计风功率预测误差的公式如(6)所示:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>&lsqb;</mo><msub><mi>e</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>,</mo><msub><mi>e</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&rsqb;</mo><mo>=</mo><mo>&lsqb;</mo><msub><mi>&alpha;</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mover><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>&beta;&sigma;</mi><mi>i</mi></msub><mo>+</mo><msub><mi>K</mi><mn>3</mn></msub><mi>S</mi><mo>+</mo><msub><mi>K</mi><mn>4</mn></msub><mi>&epsiv;</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>&alpha;</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mover><msub><mi>P</mi><mrow><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>&beta;&sigma;</mi><mi>i</mi></msub><mo>+</mo><msub><mi>K</mi><mn>3</mn></msub><mi>S</mi><mo>+</mo><msub><mi>K</mi><mn>4</mn></msub><mi>&epsiv;</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000939169750000034.GIF" wi="1795" he="95" /></maths>其中,i表示对日前预测出力数据所分的段数,<img file="FDA0000939169750000035.GIF" wi="123" he="85" />表示第i个所分的时间段中日前预测风功率P<sub>fore,i</sub>的平均值,e<sub>max,i</sub>和e<sub>min,i</sub>分别表示第i个所分的时间段中风功率预测误差的上限值和下限值,其他参数同上所述。
地址 210096 江苏省南京市四牌楼2号