发明名称 一种光伏输出功率值的分频预测方法
摘要 本发明公开了一种光伏输出功率值的分频预测方法,采用了基于灰色系统校正的小波神经网络,能从历史光伏数据中提取出预测日的出力趋势及特征,提高预测精度;采用了基于小波分解对历史数据进行分解,能对历史光伏数据进行多尺度分析,有效提取历史光伏数据的局部信息;采用了基于神经网络构建预测模型,有助于提取各种外界条件下光伏电站的出力趋势,提高了非晴天条件下光伏输出功率的预测精度;采用了基于灰色系统模型对历史数据进行建模,通过对历史光伏数据的误差分析,找出预测误差分布范围,为灰色系统校正提供了依据;采用了基于历史数据样本进行灰色建模,并对神经网络的预测结果进行进一步的校正,提高了神经网络预测结果的预测精度。
申请公布号 CN105117975A 申请公布日期 2015.12.02
申请号 CN201510478986.8 申请日期 2015.08.06
申请人 国家电网公司;江苏省电力公司;江苏省电力公司电力科学研究院 发明人 朱卫平;袁晓冬;陈兵;孙健;戴强晟;李强;史明明;贾萌萌;柳丹;吕振华;罗珊珊
分类号 G06Q50/06(2012.01)I 主分类号 G06Q50/06(2012.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 一种光伏输出功率值的分频预测方法,其特征在于:包括以下步骤,步骤A,获取相似日数据;A1)选取光伏发电设备的历史光伏数据,构造历史样本集S={s<sub>1</sub>,s<sub>2</sub>,…,s<sub>i</sub>,…,s<sub>n</sub>};其中,s<sub>i</sub>表示历史样本集中的第i个样本;1≤i≤n;并有s<sub>i</sub>={w<sub>i</sub>,x<sub>i,max</sub>,x<sub>i,min</sub>,p<sub>i,1</sub>,p<sub>i,2</sub>,…,p<sub>i,l</sub>,…,p<sub>i,m</sub>};w<sub>i</sub>表示历史样本s<sub>i</sub>的天气类型;x<sub>i,max</sub>表示历史样本s<sub>i</sub>的最高温度;x<sub>i,min</sub>表示历史样本s<sub>i</sub>的最低温度;p<sub>i,1</sub>,p<sub>i,2</sub>,…,p<sub>i,l</sub>,…,p<sub>i,m</sub>表示历史样本s<sub>i</sub>各时刻的光伏输出功率值;p<sub>i,l</sub>表示历史样本s<sub>i</sub>在第l个时刻的光伏输出功率值;1≤l≤m;A2)假设光伏发电设备的预测日光伏数据为{w<sub>p</sub>,x<sub>p,max</sub>,x<sub>p,min</sub>};其中,w<sub>p</sub>表示预测日的天气类型;x<sub>p,max</sub>表示预测日的最高温度;x<sub>p,min</sub>表示预测日的最低温度;A3)计算历史样本集S与预测日光伏数据的温度欧式距离,构造温度欧式距离集D={d<sub>1</sub>,d<sub>2</sub>,…,d<sub>i</sub>,…,d<sub>n</sub>};其中,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>d</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>&lsqb;</mo><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>p</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>p</mi><mo>,</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>,</mo><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mrow>]]></math><img file="FDA0000775743550000021.GIF" wi="1425" he="205" /></maths>d<sub>i</sub>表示历史样本s<sub>i</sub>与预测日光伏数据的温度欧式距离;A4)对温度欧式距离集D进行降序排序,选择前r组温度欧式距离对应的历史样本,构造与所述预测日光伏数据相对应的神经网络训练样本集S<sub>c</sub>;S<sub>c</sub>={s<sub>c,1</sub>,s<sub>c,2</sub>,…,s<sub>c,j</sub>,…,s<sub>c,r</sub>},其中s<sub>c,j</sub>表示第j个神经网络训练样本,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>s</mi><mrow><mi>c</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mo>&lsqb;</mo><msup><mi>w</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msup><mo>,</mo><msubsup><mi>x</mi><mrow><mi>c</mi><mo>,</mo><mi>max</mi></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>x</mi><mrow><mi>c</mi><mo>,</mo><mi>min</mi></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>p</mi><mrow><mi>c</mi><mo>,</mo><mn>1</mn></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>p</mi><mrow><mi>c</mi><mo>,</mo><mn>2</mn></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>p</mi><mrow><mi>c</mi><mo>,</mo><mi>l</mi></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>p</mi><mrow><mi>c</mi><mo>,</mo><mi>m</mi></mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>&rsqb;</mo><mo>,</mo></mrow>]]></math><img file="FDA0000775743550000022.GIF" wi="1692" he="171" /></maths>1≤j≤r,w<sup>(j)</sup>表示神经网络训练样本s<sub>c,j</sub>的天气类型;<img file="FDA0000775743550000023.GIF" wi="199" he="147" />表示神经网络训练样本s<sub>c,j</sub>的最高温度;<img file="FDA0000775743550000024.GIF" wi="189" he="148" />表示神经网络训练样本s<sub>c,j</sub>的最低温度;<img file="FDA0000775743550000025.GIF" wi="983" he="165" />表示神经网络训练样本s<sub>c,j</sub>各时刻的光伏输出功率值;<img file="FDA0000775743550000026.GIF" wi="167" he="164" />表示神经网络训练样本s<sub>c,j</sub>在第l个时刻的光伏输出功率值;步骤B,生成初始预测值;B1)对神经网络训练样本集S<sub>c</sub>中的每个元素进行λ次小波分解,获得神经网络训练样本集S<sub>c</sub>的分解小波系数序列集;对s<sub>c,j</sub>进行λ次小波分解,获得第j个分解小波系数序列<img file="FDA0000775743550000027.GIF" wi="909" he="149" />其中,<img file="FDA0000775743550000028.GIF" wi="143" he="137" />表示对s<sub>c,j</sub>进行第λ次小波分解所获得的小波系数序列;<img file="FDA0000775743550000031.GIF" wi="128" he="124" />表示对s<sub>c,j</sub>进行第λ次小波分解所获得的平均系数序列;B2)以第j‑1个分解小波系数序列<img file="FDA0000775743550000032.GIF" wi="969" he="141" />作为神经网络的输入信号,以第j个分解小波系数序列<img file="FDA0000775743550000033.GIF" wi="704" he="129" />作为神经网络的输出信号,以所述输入信号和所述输出信号作为神经网络的学习样本,通过神经网络对所述学习样本进行训练,从而获得光伏输出功率值分频预测模型;B3)以第r个分解小波系数序列<img file="FDA0000775743550000034.GIF" wi="808" he="143" />作为光伏输出功率值分频预测模型的输入信号,通过光伏输出功率值分频预测模型,获得初始预测小波序列<img file="FDA0000775743550000035.GIF" wi="968" he="166" />B4)对初始预测小波序列<img file="FDA0000775743550000036.GIF" wi="933" he="147" />进行小波重组,获得初始光伏输出功率预测序列{p<sub>p,1</sub>,p<sub>p,2</sub>,...,p<sub>p,l</sub>,...,p<sub>p,m</sub>},p<sub>p,l</sub>表示预测日在第l个时刻的光伏输出功率预测值;步骤C,确立灰色系统预测值的判断区间;C1)根据历史样本集获得预测序列X(0)=[x<sup>(0)</sup>(1),x<sup>(0)</sup>(2),…,x<sup>(0)</sup>(k),…,x<sup>(0)</sup>(t)];其中,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>x</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>t</mi><mo>,</mo></mrow>]]></math><img file="FDA0000775743550000037.GIF" wi="975" he="226" /></maths>t表示与预测日相邻的前t天;C2)利用预测序列X(0)构建灰色系统DGM模型;C3)对n个历史样本,利用灰色系统DGM模型分别进行灰色预测,获得n个预测值W,并根据n个预测值W计算n个历史样本的误差值;C4)根据n个历史样本的误差值选取误差区间E<sub>g</sub>=(δ,ε);其中δ,ε是由n个历史样本的误差值集中分布范围决定的,δ是误差区间下限,ε是误差区间上限。C5)将t+1的值赋给k,从而利用灰色系统DGM模型获得预测日的预测序列值W<sub>g</sub>=x<sup>(0)</sup>(t+1),获得判断区间J<sub>g</sub>=(W<sub>g</sub>‑W<sub>g</sub>ε,W<sub>g</sub>‑W<sub>g</sub>δ)∪(W<sub>g</sub>+W<sub>g</sub>δ,W<sub>g</sub>+W<sub>g</sub>ε);步骤D,判断<img file="FDA0000775743550000041.GIF" wi="386" he="196" />是否成立,若成立,则不进行校正,从而以预测日在第l个时刻的光伏输出功率预测值p<sub>p,l</sub>作为最终预测值;否则,利用校正公式进行校正,从而获得预测日在第l个时刻的光伏输出功率预测校正值p<sub>p,lnew</sub>,并以预测校正值p<sub>p,lnew</sub>作为最终预测值;校正公式为,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>p</mi><mrow><mi>p</mi><mo>,</mo><mi>l</mi><mi>n</mi><mi>e</mi><mi>w</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>W</mi><mi>g</mi></msub><mrow><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>p</mi><mrow><mi>p</mi><mo>,</mo><mi>l</mi></mrow></msub></mrow></mfrac><mo>&times;</mo><msub><mi>p</mi><mrow><mi>p</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>.</mo></mrow>]]></math><img file="FDA0000775743550000042.GIF" wi="1255" he="444" /></maths>
地址 211103 江苏省南京市江宁区帕威尔路1号
您可能感兴趣的专利