发明名称 基于卡尔曼滤波的风电场中储能系统充放电控制方法
摘要 本发明涉及一种基于卡尔曼滤波的风电场中储能系统充放电控制方法,属于电力系统控制技术领域。首先计算k时刻风电场中风力发电系统与储能系统的联合输出功率的先验估计值和计算误差协方差先验估计值,计算当前风力发电系统的功率波动量,并根据其确定观测噪声协方差值,根据k时刻风力发电系统功率的实测值,对联合输出功率的先验估计值进行修正,得到后验估计值和计算误差协方差的后验估计值;根据储能系统的剩余能量和输出功率的限值,计算k时刻风电场中储能系统输出功率的上、下限。本方法综合考虑风电波动情况、储能系统能量和功率情况,自适应地调整滤波器参数以获取更好的滤波效果同时减小所需的储能系统容量。
申请公布号 CN103618326B 申请公布日期 2015.08.19
申请号 CN201310571570.1 申请日期 2013.11.13
申请人 清华大学;中国南方电网有限责任公司调峰调频发电公司 发明人 袁志昌;张百华;陆超;郭海峰;徐罕晨;李勇琦;陈满
分类号 H02J3/28(2006.01)I 主分类号 H02J3/28(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 一种基于卡尔曼滤波的风电场中储能系统充放电控制方法,其特征在于该方法包括以下步骤:(1)在k‑1时刻,计算k时刻风电场中风力发电系统与储能系统的联合输出功率的先验估计值(P<sub>O</sub>)<sub>k|k‑1</sub>以及联合输出功率的计算误差协方差P<sub>k|k‑1</sub>的先验估计值:(P<sub>O</sub>)<sub>k|k‑1</sub>=(P<sub>O</sub>)<sub>k‑1</sub>,P<sub>k|k‑1</sub>=P<sub>k‑1</sub>+Q,其中,P<sub>O</sub>是风电场中风力发电系统与储能系统的联合输出功率,(P<sub>O</sub>)<sub>k‑1</sub>是k‑1时刻风电场中风力发电系统与储能系统的联合输出功率的后验估计值,P是联合输出功率的计算误差协方差,P<sub>k‑1</sub>是k‑1时刻联合输出功率的计算误差协方差的后验估计值,Q是过程激励噪声协方差,Q=1×10<sup>‑4</sup>,初始时刻的联合输出功率(P<sub>O</sub>)<sub>0</sub>为初始时刻的风电功率测量值(P<sub>W</sub>)<sub>0</sub>,并设初始计算误差协方差P<sub>0</sub>为0;(2)在k时刻,计算当前风力发电系统的功率波动量F(k):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>max</mi><mrow><mo>(</mo><msub><mi>P</mi><mi>W</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>P</mi><mi>O</mi><mi>min</mi></msubsup><mo>,</mo><msubsup><mi>P</mi><mi>O</mi><mi>max</mi></msubsup><mo>-</mo><msub><mi>P</mi><mi>W</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000415637620000015.GIF" wi="815" he="67" /></maths>其中,P<sub>W</sub>为风力发电系统功率,P<sub>W</sub>(k)为k时刻风力发电系统功率,<img file="FDA0000415637620000011.GIF" wi="104" he="72" />和<img file="FDA0000415637620000012.GIF" wi="92" he="79" />分别为k时刻的前T时段内联合输出功率的最大值和最小值;根据上述计算得到的F(k)的值,由下表获取观测噪声协方差R的值:<tables num="0001" id="ctbl0001"><img file="FDA0000415637620000013.GIF" wi="1987" he="227" /></tables>其中,F<sub>max</sub>为T时段内风力发电系统功率的最大波动量,根据上述观测噪声协方差R,计算k时刻的卡尔曼增益Kk:<img file="FDA0000415637620000014.GIF" wi="329" he="151" />其中,P<sub>k|k‑1</sub>为步骤(1)中计算得到的k时刻计算误差协方差的先验估计值;(3)根据k时刻风力发电系统功率的实测值(P<sub>W</sub>)<sub>k</sub>,对上述联合输出功率的先验估计值(P<sub>O</sub>)<sub>k|k‑1</sub>进行修正,得到联合输出功率的后验估计值(P<sub>O</sub>)<sub>k</sub>和k时刻计算误差协方差的后验估计值P<sub>k</sub>:(P<sub>O</sub>)<sub>k</sub>=(P<sub>O</sub>)<sub>k|k‑1</sub>+K<sub>k</sub>((P<sub>W</sub>)<sub>k</sub>‑(P<sub>O</sub>)<sub>k|k‑1</sub>),P<sub>k</sub>=(1‑K<sub>k</sub>)P<sub>k|k‑1</sub>,其中,(P<sub>O</sub>)<sub>k|k‑1</sub>是步骤(1)中联合输出功率的先验估计值,(P<sub>W</sub>)<sub>k</sub>是k时刻的风力发电系统功率实测值,K<sub>k</sub>是步骤(2)中计算得到的卡尔曼增益,P<sub>k|k‑1</sub>为步骤(1)中计算得到的计算误差协方差的先验估计值;(4)根据测量得到的(P<sub>W</sub>)<sub>k</sub>和上述计算得到的(P<sub>O</sub>)<sub>k</sub>,计算k时刻储能系统的输出功率(P<sub>ES</sub>)<sub>k</sub>:(P<sub>ES</sub>)<sub>k</sub>=(P<sub>O</sub>)<sub>k</sub>‑(P<sub>W</sub>)<sub>k</sub>;(5)根据储能系统的剩余能量和输出功率的限值,计算k时刻风电场中储能系统输出功率的上、下限:其中,上限U为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>U</mi><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msubsup><mi>P</mi><mrow><mi>ES</mi><mo>,</mo><mi>disch</mi></mrow><mi>max</mi></msubsup><mo>,</mo><mfrac><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>-</mo><msub><mi>E</mi><mi>min</mi></msub></mrow><mi>&Delta;t</mi></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000415637620000021.GIF" wi="587" he="129" /></maths>下限L为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>L</mi><mo>=</mo><mi>max</mi><mrow><mo>(</mo><msubsup><mi>P</mi><mrow><mi>ES</mi><mo>,</mo><mi>disch</mi></mrow><mi>min</mi></msubsup><mo>,</mo><mfrac><mrow><msub><mi>E</mi><mi>k</mi></msub><mo>-</mo><msub><mi>E</mi><mi>max</mi></msub></mrow><mi>&Delta;t</mi></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000415637620000022.GIF" wi="593" he="129" /></maths>当上述计算得到的(P<sub>ES</sub>)<sub>k</sub>大于U时,将U作为k时刻风电场中储能系统的输出功率,当上述计算得到的(P<sub>ES</sub>)<sub>k</sub>小于L时,将L作为k时刻风电场中储能系统的输出功率,其中,Δt为风电场中储能系统的控制周期,Δt为(k‑1)~k,E<sub>k</sub>为k时刻储能系统的剩余能量,<img file="FDA0000415637620000023.GIF" wi="117" he="77" />和<img file="FDA0000415637620000024.GIF" wi="144" he="77" />分别为储能系统的最大充电功率和最大放电功率,E<sub>min</sub>和E<sub>max</sub>分别为储能系统的最小和最大剩余能量;(6)重复步骤(1)‑(5),实现对风电场中储能系统的充放电控制。
地址 100084 北京市海淀区清华园1号