发明名称 基于时间序列模型的光伏发电功率超短期预测方法
摘要 本发明公开了一种基于时间序列模型的光伏发电功率超短期预测方法,其特征在于,包括以下步骤:步骤一:采集光伏电站历史功率数据,并对历史功率数据进行归一化处理。步骤二:根据归一化后历史功率数据建立拟合方程,根据建立的拟合方程和残差方差确定模型的阶数即p和q的值;步骤三:确定<img file="2013104168327100004DEST_PATH_IMAGE002.GIF" wi="114" he="22" />;<img file="2013104168327100004DEST_PATH_IMAGE004.GIF" wi="93" he="24" />的值;步骤四:建立自回归滑动平均模型。通过历史功率数据、以及拟合方程和自回归滑动平均模型方程对光伏发电功率超短期预测建立预测模型,根据预测模型和现有的数据即可达到对光伏发电功率进行短期准确预测的目的。
申请公布号 CN103473322A 申请公布日期 2013.12.25
申请号 CN201310416832.7 申请日期 2013.09.13
申请人 国家电网公司;甘肃省电力公司;甘肃省电力公司风电技术中心 发明人 路亮;汪宁渤;李照荣;赵龙;王有生;刘光途
分类号 G06F17/30(2006.01)I 主分类号 G06F17/30(2006.01)I
代理机构 北京中恒高博知识产权代理有限公司 11249 代理人 姜万林
主权项 1.一种基于时间序列模型的光伏发电功率超短期预测方法,其特征在于,包括以下步骤:步骤一:采集光伏电站历史功率数据,并对历史功率数据进行归一化处理,归一化处理公式如下:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>x</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>min</mi></msub></mrow><mrow><msub><mi>x</mi><mi>max</mi></msub><mo>-</mo><msub><mi>x</mi><mi>min</mi></msub></mrow></mfrac></mrow></math>]]></maths>其中,x<sub>max</sub>是历史功率数据中的最大值,x<sub>min</sub>是历史功率数据中的最小值,x<sub>i</sub>是第i个历史功率数据,x′<sub>i</sub>是归一化后的第i个历史功率数据;步骤二:根据上述归一化后历史功率数据建立拟合方程,根据建立的拟合方程和残差方差确定模型的阶数即p和q的值;步骤三:如步骤一中采集的光伏电站历史功率数据利用数据序列x<sub>1</sub>,x<sub>2</sub>,…,x<sub>t</sub>表示,其样本自协方差定义为<maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>&gamma;</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>t</mi></msub><msub><mi>x</mi><mrow><mi>t</mi><mo>-</mo><mi>k</mi></mrow></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo></mrow></math>]]></maths>其中<img file="FDA0000381692740000013.GIF" wi="411" he="199" />x<sub>t</sub>和x<sub>t-k</sub>均为数据序列x<sub>1</sub>,x<sub>2</sub>,…,x<sub>t</sub>中的数值;则历史功率数据样本自相关函数为:<maths num="0003"><![CDATA[<math><mrow><msub><mover><mi>&rho;</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><mfrac><msub><mover><mi>&gamma;</mi><mo>^</mo></mover><mi>k</mi></msub><msub><mover><mi>&gamma;</mi><mo>^</mo></mover><mn>0</mn></msub></mfrac><mo>=</mo><mfrac><mrow><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>t</mi></msub><msub><mi>x</mi><mrow><mi>t</mi><mo>-</mo><mi>k</mi></mrow></msub></mrow><mrow><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>x</mi><mi>t</mi><mn>2</mn></msubsup></mrow></mfrac><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>t</mi></msub><msub><mi>x</mi><mrow><mi>t</mi><mo>-</mo><mi>k</mi></mrow></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>x</mi><mi>t</mi><mn>2</mn></msubsup></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>k</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></math>]]></maths>则自回归部分的矩估计为<img file="FDA0000381692740000023.GIF" wi="1049" he="351" />令<img file="FDA0000381692740000024.GIF" wi="925" he="79" />协方差函数为<img file="FDA0000381692740000025.GIF" wi="1662" he="537" />用<img file="FDA0000381692740000026.GIF" wi="42" he="56" />的估计代替γ<sub>k</sub>,则<img file="FDA0000381692740000027.GIF" wi="720" he="201" />即得出,<img file="FDA0000381692740000028.GIF" wi="393" he="93" />对滑动平均模型系数θ<sub>1</sub>,θ<sub>2</sub>,...,θ<sub>q</sub>采用矩估计,有<maths num="0005"><![CDATA[<math><mrow><msub><mi>&gamma;</mi><mn>0</mn></msub><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msubsup><mi>&theta;</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&theta;</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msubsup><mi>&theta;</mi><mi>q</mi><mn>2</mn></msubsup><mo>)</mo></mrow><msubsup><mi>&sigma;</mi><mi>a</mi><mn>2</mn></msubsup></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>&gamma;</mi><mi>k</mi></msub><mo>=</mo><mrow><mo>(</mo><mo>-</mo><msub><mi>&theta;</mi><mi>k</mi></msub><mo>+</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&theta;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>&theta;</mi><mrow><mi>q</mi><mo>-</mo><mi>k</mi></mrow></msub><msub><mi>&theta;</mi><mi>q</mi></msub><mo>)</mo></mrow><msubsup><mi>&sigma;</mi><mi>a</mi><mn>2</mn></msubsup><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi></mrow></math>]]></maths>以上共包含m+1个方程,对其参数而言,方程为非线性,采用迭代法进行求解即得出θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>q</sub>的值;步骤四:根据上述步骤二和步骤三得出的p和q值以及<img file="FDA0000381692740000034.GIF" wi="368" he="92" />、θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>q</sub>的值建立自回归滑动平均模型;所述自回归滑动平均模型的公式如下:<img file="FDA0000381692740000035.GIF" wi="1110" he="220" />其中,<img file="FDA0000381692740000036.GIF" wi="405" he="86" />和θ<sub>j</sub>(1≤j≤q)是系数,α<sub>t</sub>是白噪声序列。
地址 100031 北京市西城区西长安街86号