发明名称 一种基于相似数据选取的光伏电站输出功率预测方法
摘要 本发明涉及一种基于相似数据选取的光伏电站输出功率预测方法,属于光伏发电技术领域。该方法包括以下步骤:步骤1:采集历史天的辐照强度值、温度值和实际光伏输出功率值,以及气象预报中预测天辐照强度值和温度值;步骤2:确定每天6点到18点各整点辐照强度对应的权值w<sub>1</sub>(i),确定每天6点到18点各整点温度对应的权值w<sub>2</sub>(i);步骤3:进行相似天的选取;步骤4:根据所选相似天与预测天的相关程度,确定功率预测时各相似天功率所占权重;步骤5:计算得到预测天进行光伏输出功率预测值,并对预测结果进行评估;本发明较好的挖掘了预测天和历史数据的相关性,易于实现,提高了光伏输出功率预测的准确性。
申请公布号 CN103390202B 申请公布日期 2016.01.20
申请号 CN201310303922.5 申请日期 2013.07.18
申请人 华北电力大学 发明人 杨锡运;刘欢;姜飞飞;武晓宁;肖运启
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 黄家俊
主权项 一种基于相似数据选取的光伏电站输出功率预测方法,其特征在于,该方法包括以下步骤:步骤1:采集历史天的辐照强度值、温度值和实际光伏输出功率值,以及气象预报中预测天辐照强度值和温度值;(1)采集与预测天临近的过去N天的历史数据每天6点到18点各整点的辐照强度值,组成一组辐照强度向量:I<sub>1n</sub>[I<sub>1n</sub>(1),I<sub>1n</sub>(2),…I<sub>1n</sub>(12),I<sub>1n</sub>(13)],(n=1…N);其中,N为设定值;则预测天临近的过去N天中的第n天6点到18点各整点的辐照强度值表示为I<sub>1n</sub>(i),i=1,2,......,13;(2)采集与预测天临近的过去N天的历史数据每天6点到18点各整点的温度值,组成一组温度向量:I<sub>2n</sub>[I<sub>2n</sub>(1),I<sub>2n</sub>(2),…I<sub>2n</sub>(12),I<sub>2n</sub>(13)],(n=1…N);则预测天临近的过去N天中的第n天6点到18点各整点的温度值表示为I<sub>2n</sub>(i),i=1,2,......,13;(3)采集与预测天临近N天的实际光伏输出功率,第n天实际光伏输出功率向量为P<sub>n</sub>(n=1…N),每隔15分钟采集一次,每天采集96个点,则第n天每个采集时刻的功率值表示为P<sub>n</sub>(j),j=1,2,......,96;(4)从气象预报中提取预测天6点到18点各整点辐照强度值,组成一组预测天辐照强度向量:I<sub>1o</sub>[I<sub>1o</sub>(1),I<sub>1o</sub>(2),…I<sub>1o</sub>(12),I<sub>1o</sub>(13)];预测天6点到18点各整点辐照强度值表示为I<sub>1o</sub>(i),i=1,2,......,13,并定义预测天辐照强度误差系数为ε<sub>1</sub>,即代表天气预报给出的辐照强度的可信程度;(5)从气象预报中提取预测天6点到18点各整点温度值,组成一组预测天温度向量:I<sub>2o</sub>[I<sub>2o</sub>(1),I<sub>2o</sub>(2),…I<sub>2o</sub>(12),I<sub>2o</sub>(13)];预测天6点到18点各整点温度值表示为I<sub>2o</sub>(i),i=1,2,......,13,并定义预测天温度误差系数为ε<sub>2</sub>,即代表天气预报中给出的温度的可信程度;步骤2:基于光伏输出功率值与辐照强度变化的关系,确定每天6点到18点各整点辐照强度对应的权值w<sub>1</sub>=[0.001,0.001,0.028,0.06,0.09,0.13,0.15,0.15,0.15,0.15,0.05,0.039,0.001]',i=1,2,......,13;基于光伏输出功率值与温度变化的关系,确定每天6点到18点各整点温度对应的权值w<sub>2</sub>=[0.01,0.015,0.02,0.05,0.09,0.13,0.15,0.15,0.15,0.145,0.05,0.03,0.01]',i=1,2,......,13;步骤3:利用与预测天临近的过去N天的辐照强度向量I<sub>1n</sub>、与预测天临近的过去N天的温度向量I<sub>2n</sub>、预测天辐照强度向量I<sub>1o</sub>、预测天温度向量I<sub>2o</sub>、预测天辐照强度误差系数ε<sub>1</sub>和预测天温度误差系数ε<sub>2</sub>,以及各整点辐照强度对应的权值w<sub>1</sub>(i)和各整点温度对应的权值w<sub>2</sub>(i)进行相似天的选取;(1)求取预测天辐照强度向量I<sub>1o</sub>和预测天临近的N个历史天辐照强度向量I<sub>1n</sub>(n=1…N)的差值绝对值向量,得到每个历史天6点到18点各整点所对应的辐照强度差值绝对值向量:δ<sub>1n</sub>=abs[I<sub>1o</sub>(1)‑I<sub>1n</sub>(1),I<sub>1o</sub>(2)‑I<sub>1n</sub>(2),…I<sub>1o</sub>(12)‑I<sub>1n</sub>(12),I<sub>1o</sub>(13)‑I<sub>1n</sub>(13)];其中,δ<sub>1n</sub>代表第n个历史天辐照强度差值绝对值向量;abs[......]代表对预测天辐照强度向量I<sub>1o</sub>和与预测天临近的N个历史天中的第n天的辐照强度向量I<sub>1n</sub>的差值取绝对值;则预测天临近的N个历史天中第n天6点到18点各整点所对应的辐照强度差值绝对值表示为:δ<sub>1n</sub>(i),i=1,2,......,13;(2)求取预测天温度向量I<sub>2o</sub>和预测天临近的N个历史天温度向量I<sub>2n</sub>(n=1…N)的差值绝对值向量,得到每个历史天6点到18点各整点所对应的温度差值绝对值向量:δ<sub>2n</sub>=abs[I<sub>2o</sub>(1)‑I<sub>2n</sub>(1),I<sub>2o</sub>(2)‑I<sub>2n</sub>(2),…I<sub>2o</sub>(12)‑I<sub>2n</sub>(12),I<sub>2o</sub>(13)‑I<sub>2n</sub>(13)];其中,δ<sub>2n</sub>代表第n个历史天温度差值绝对值向量;abs[......]代表对预测天温度向量I<sub>2o</sub>和预测天临近的N个历史天温度向量I<sub>2n</sub>的差值取绝对值;则预测天临近的N个历史天中第n天6点到18点各整点所对应的温度差值绝对值表示为:δ<sub>2n</sub>(i),i=1,2,......,13;(3)利用各整点辐照强度对应的权值w<sub>1</sub>(i)、各整点温度对应的权值w<sub>2</sub>(i),每个历史天6点到18点各整点所对应的辐照强度差值绝对值δ<sub>1n</sub>(i)和每个历史天6点到18点各整点所对应的温度差值绝对值δ<sub>2n</sub>(i),计算得到第n个历史天6点到18点各整点所对应的差值绝对值修正值:a<sub>n</sub>(i)=δ<sub>1n</sub>(i)×w<sub>1</sub>(i)×ε<sub>1</sub>×50%+δ<sub>2n</sub>(i)×w<sub>2</sub>(i)×ε<sub>2</sub>×50%,i=1…13;再计算第n个历史天的各整点差值绝对值修正值的和:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>S</mi><mi>n</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>13</mn></munderover><msub><mi>a</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000850654100000031.GIF" wi="305" he="141" /></maths>(4)将S<sub>n</sub>(n=1…N)从小到大排序,取排好序列的前三个S<sub>n</sub>的值对应的天作为预测天的相似天;步骤4:根据所选相似天与预测天的相关程度,确定功率预测时各相似天功率所占权重;(1)假设所选相似天各整点差值绝对值修正值的和分别为S<sub>a</sub>、S<sub>b</sub>、S<sub>c</sub>,S<sub>a</sub>&lt;S<sub>b</sub>&lt;S<sub>c</sub>;构造向量<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>p</mi><mo>=</mo><mrow><mo>(</mo><mfrac><msub><mi>S</mi><mi>c</mi></msub><msub><mi>S</mi><mi>a</mi></msub></mfrac><mo>,</mo><mfrac><msub><mi>S</mi><mi>c</mi></msub><msub><mi>S</mi><mi>b</mi></msub></mfrac><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000850654100000032.GIF" wi="328" he="141" /></maths>其中<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>p</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mi>S</mi><mi>c</mi></msub><msub><mi>S</mi><mi>a</mi></msub></mfrac><mo>,</mo><msub><mi>p</mi><mn>2</mn></msub><mo>=</mo><mfrac><msub><mi>S</mi><mi>c</mi></msub><msub><mi>S</mi><mi>b</mi></msub></mfrac><mo>,</mo><msub><mi>p</mi><mn>3</mn></msub><mo>=</mo><mn>1</mn><mo>;</mo></mrow>]]></math><img file="FDA0000850654100000033.GIF" wi="607" he="143" /></maths>(2)计算构造向量p各元素的和<img file="FDA0000850654100000034.GIF" wi="389" he="148" />(3)计算各相似天的权重向量<img file="FDA0000850654100000041.GIF" wi="423" he="119" />其中,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>r</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mi>p</mi><mn>1</mn></msub><mrow><mi>s</mi><mi>u</mi><mi>m</mi></mrow></mfrac><mo>,</mo><msub><mi>r</mi><mn>2</mn></msub><mo>=</mo><mfrac><msub><mi>p</mi><mn>2</mn></msub><mrow><mi>s</mi><mi>u</mi><mi>m</mi></mrow></mfrac><mo>,</mo><msub><mi>r</mi><mn>3</mn></msub><mo>=</mo><mfrac><msub><mi>p</mi><mn>3</mn></msub><mrow><mi>s</mi><mi>u</mi><mi>m</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000850654100000042.GIF" wi="673" he="119" /></maths>其中,p<sub>1</sub>、p<sub>2</sub>和p<sub>3</sub>是差值绝对值修正和S<sub>a</sub>、S<sub>b</sub>和S<sub>c</sub>相对应的构造向量值;r<sub>1</sub>、r<sub>2</sub>和r<sub>3</sub>代表与差值绝对值修正和S<sub>a</sub>、S<sub>b</sub>和S<sub>c</sub>相对应的权重值;步骤5:计算得到预测天进行光伏输出功率预测值,并对预测结果进行评估;(1)预测光伏输出功率值的计算公式为:P<sub>f</sub>=r<sub>1</sub>P<sub>1</sub>+r<sub>2</sub>P<sub>2</sub>+r<sub>3</sub>P<sub>3</sub>;其中,P<sub>f</sub>表示预测光伏输出功率值;P<sub>1</sub>为选出的最相似天的功率,P<sub>2</sub>为选出的次相似天的功率,P<sub>3</sub>为选出的第3个相似天的功率;(2)计算预测值的均方根误差:<img file="FDA0000850654100000043.GIF" wi="460" he="235" />其中,e<sub>j</sub>=P<sub>f</sub>(j)‑P<sub>n</sub>(j);其中,e<sub>j</sub>代表每一个预测时刻预测功率与实际功率的差值,用P<sub>f</sub>(j)表示第j个的预测光伏输出功率值,P<sub>总</sub>表示总装机的容量;rmse表示预测值的均方根误差,P<sub>n</sub>(j)表示第n天第j个预测时刻的光伏实际输出功率值。
地址 102206 北京市昌平区朱辛庄北农路2号