发明名称 一种用于风电场实时功率预报修正的电池能量管理方法
摘要 本发明公开了一种用于风电场实时功率预报修正的电池能量管理方法,所述方法不仅从功率预报误差修正的主动需求出发,在考虑变换器效率及自放电因素下,对基本处于连续不完全充放电交替工作复杂状态的电池荷电状态进行动态预测与管理。而且结合电池实际充放电电流时间积分监测量,进一步修正实际SOC变化,同时修正输入输出电能与电池荷电量当量关系参数,修正变换器效率及自放电效率,判断电池健康使用情况,统计使用周期数,方便后续的SOC预测。本发明带来的有益效果是:用电池充放电过程外围电路输入输出参数来估计电池SOC动态变化,解决了电池模型构建难的问题,实现电池科学管理,有效延长电池寿命,达到风电场预报精准要求。
申请公布号 CN105207243A 申请公布日期 2015.12.30
申请号 CN201510651416.4 申请日期 2015.10.10
申请人 河海大学 发明人 朱建红;潘文霞
分类号 H02J3/32(2006.01)I;H02J3/38(2006.01)I 主分类号 H02J3/32(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 一种用于风电场实时功率预报修正的电池能量管理方法,其特征在于,包括如下步骤:1)通过实际风电功率采样周期内功率预报误差修正所需消纳的电能值来预测此采样周期T内电池所需充放的荷电量,估计本采样周期T末的电池荷电状态,并结合监测的电能及荷电量变化实际值,动态修正电池的荷电状态;1.1)离线对新进电池参数检测,根据新电池出产时的数据,在正式投入使用前通过1到2次完整充放电实验获取;1.2)系统实际投入运行时,由风功率预报采样周期T点的预测值与实际采样值的误差数据,对采样周期T内需输入输出做功电能值进行计算;W<sub>e</sub>(nT)表示预报功率与实际功率差值在第n个预报周期内的做功情况;P<sub>e</sub>(nT)为第n个预报采样点预报功率与实际功率的误差;每个预报采样点内有m个实际功率采样周期,实际采样周期<img file="FDA0000818549250000011.GIF" wi="210" he="131" />l为m个采样周期的任一点,<img file="FDA0000818549250000012.GIF" wi="222" he="122" />表示T周期内的第l点的预报功率;每两个预报点之间的功率采用直线拟合,用两点直线拟合P<sub>(n‑1)T</sub>到P<sub>nT</sub>的变化轨迹,并求取两预报点之间的动态预报数据;1.2.1)由预报误差预测输入输出电能变化,<img file="FDA0000818549250000013.GIF" wi="241" he="117" />为正,表示实际风力发电功率小于预报设定值,需要储能系统提供的能量,此时<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>W</mi><mrow><mi>e</mi><mrow><mo>{</mo><mrow><mfrac><mi>l</mi><mi>m</mi></mfrac><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>e</mi><mrow><mo>{</mo><mrow><mfrac><mi>l</mi><mi>m</mi></mfrac><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></mrow></msub><mo>*</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><mi>T</mi><mo>/</mo><mi>&eta;</mi><mo>;</mo></mrow>]]></math><img file="FDA0000818549250000014.GIF" wi="718" he="154" /></maths><img file="FDA0000818549250000015.GIF" wi="238" he="117" />为负,表示实际风力发电功率大于预报设定值,需要储能系统储存的能量,此时<img file="FDA0000818549250000016.GIF" wi="733" he="150" />对采样周期T内电池荷电情况变化进行预测,并根据实际电流时间积分,修正当量系数k;如果<img file="FDA0000818549250000021.GIF" wi="236" he="132" />为正,则充电电流i<sub>c</sub>为零,储能系统需要放电,存在放电电流i<sub>d</sub>,按放电过程荷电量变化量Q<sub>d</sub>=k<sub>d</sub>W<sub>out</sub>估计荷电状态变化量<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>&Delta;SOC</mi><mo>&prime;</mo></msup><mo>=</mo><mfrac><mrow><msup><mi>&Delta;Q</mi><mo>&prime;</mo></msup></mrow><msub><mi>Q</mi><mrow><mi>k</mi><mi>t</mi></mrow></msub></mfrac><mo>=</mo><mfrac><mrow><mo>-</mo><msub><mi>k</mi><mi>d</mi></msub><msub><mi>W</mi><mrow><mi>e</mi><mrow><mo>{</mo><mrow><mfrac><mi>l</mi><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></mrow></msub></mrow><msub><mi>Q</mi><mrow><mi>k</mi><mi>t</mi></mrow></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000818549250000022.GIF" wi="704" he="196" /></maths>如果<img file="FDA0000818549250000023.GIF" wi="230" he="125" />为负,则放电电流i<sub>d</sub>为零,储能系统需要充电,存在充电电流i<sub>c</sub>,按充电过程荷电量变化量Q<sub>c</sub>=k<sub>c</sub>W<sub>in</sub>估计荷电状态变化量<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>&Delta;SOC</mi><mo>&prime;</mo></msup><mo>=</mo><mfrac><mrow><msup><mi>&Delta;Q</mi><mo>&prime;</mo></msup></mrow><msub><mi>Q</mi><mrow><mi>k</mi><mi>t</mi></mrow></msub></mfrac><mo>=</mo><mfrac><mrow><msub><mi>k</mi><mi>c</mi></msub><msub><mi>W</mi><mrow><mi>e</mi><mrow><mo>{</mo><mrow><mfrac><mi>l</mi><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></mrow></msub></mrow><msub><mi>Q</mi><mrow><mi>k</mi><mi>t</mi></mrow></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000818549250000024.GIF" wi="633" he="190" /></maths>荷电状态预测值<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msubsup><mi>SOC</mi><mrow><mo>{</mo><mrow><mfrac><mi>l</mi><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>SOC</mi><mrow><mo>{</mo><mrow><mfrac><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></msub><mo>+</mo><msup><mi>&Delta;SOC</mi><mo>&prime;</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mrow><msub><mi>Q</mi><mn>0</mn></msub><mo>+</mo><mi>&sigma;</mi><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mo>{</mo><mrow><mfrac><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></msubsup><msub><mi>i</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>d</mi><mi>t</mi><mo>-</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mo>{</mo><mrow><mfrac><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></msubsup><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>d</mi><mi>t</mi><mo>-</mo><mi>S</mi><mrow><mo>(</mo><mrow><msub><mi>T</mi><mi>a</mi></msub><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow></mrow><msub><mi>Q</mi><mrow><mi>k</mi><mi>t</mi></mrow></msub></mfrac><mo>+</mo><mfrac><mrow><msub><mi>kW</mi><mrow><mi>e</mi><mrow><mo>{</mo><mrow><mfrac><mi>l</mi><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></mrow></msub></mrow><msub><mi>Q</mi><mrow><mi>k</mi><mi>t</mi></mrow></msub></mfrac></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000818549250000025.GIF" wi="1498" he="325" /></maths>Q<sub>kt</sub>是k次完整充放电循环统计修正后的蓄电池容量;<img file="FDA0000818549250000026.GIF" wi="935" he="181" />其中,Q<sub>0</sub>表示电池初始电量;σ表示充电效率;S(T<sub>a</sub>,t)表示自放电电量;1.2.2)按实际监测数值修正相关参数:第n个预报采样周期内实际监测该周期间隔内第l个<img file="FDA0000818549250000027.GIF" wi="62" he="127" />内的电池荷电变化量<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>Q</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mo>{</mo><mrow><mfrac><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow><mrow><mo>{</mo><mrow><mfrac><mi>l</mi><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></msubsup><msub><mi>i</mi><mi>c</mi></msub><mi>d</mi><mi>t</mi><mo>-</mo><msubsup><mo>&Integral;</mo><mrow><mo>{</mo><mrow><mfrac><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow><mrow><mo>{</mo><mrow><mfrac><mi>l</mi><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>T</mi></mrow><mo>}</mo></mrow></msubsup><msub><mi>i</mi><mi>d</mi></msub><mi>d</mi><mi>t</mi><mo>,</mo></mrow>]]></math><img file="FDA0000818549250000028.GIF" wi="778" he="172" /></maths>对照先前荷电变化量预测<img file="FDA0000818549250000029.GIF" wi="422" he="123" />由ΔQ‑ΔQ′来进一步修正当量常数k<sub>d</sub>或k<sub>c</sub>;当|W<sub>int</sub>‑W<sub>outt</sub>/(η<sub>bc1</sub>*η<sub>bd1</sub>*σ<sub>1</sub>)|≈0.5W<sub>in1</sub>时,修正<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>Q</mi><mrow><mi>k</mi><mi>t</mi></mrow></msub><mo>=</mo><mn>0.5</mn><msub><mi>Q</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo><mi>t</mi></mrow></msub><mo>+</mo><mi>&Delta;</mi><mi>Q</mi><mo>=</mo><mn>0.5</mn><msub><mi>Q</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo><mi>t</mi></mrow></msub><mo>+</mo><mo>|</mo><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mi>k</mi><mi>t</mi></mrow></msubsup><mrow><msub><mi>i</mi><mi>c</mi></msub><mi>d</mi><mi>t</mi></mrow></mrow><mo>-</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mi>k</mi><mi>t</mi></mrow></msubsup><msub><mi>i</mi><mi>d</mi></msub><mi>d</mi><mi>t</mi><mo>|</mo><mo>;</mo></mrow>]]></math><img file="FDA00008185492500000210.GIF" wi="1011" he="127" /></maths>Q<sub>kt</sub>的修正在风电场故障检修时经过一次完整的充放电循环进行,k为完成循环的次数;1.3)对充放电累计值ΔQ<sub>zong</sub>进行过充过放管理:储能系统平繁充放电过程中,对应充放电电量累加统计,充电为正,放电为负,充电或放电累计值始终不能超越值Q<sub>kt</sub>,控制<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mo>|</mo><msub><mi>&Delta;Q</mi><mrow><mi>z</mi><mi>o</mi><mi>n</mi><mi>g</mi></mrow></msub><mo>|</mo><mo>=</mo><mo>|</mo><msubsup><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mrow><mo>{</mo><mfrac><mi>l</mi><mi>m</mi></mfrac><mi>T</mi><mo>+</mo><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>T</mi><mo>}</mo></mrow></msubsup><msub><mi>&Delta;Q</mi><mfrac><mi>T</mi><mi>m</mi></mfrac></msub><mo>|</mo><mo>&le;</mo><mn>0.5</mn><msub><mi>Q</mi><mrow><mi>d</mi><mn>1</mn></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA00008185492500000211.GIF" wi="759" he="193" /></maths>当电池充电过程剩余可充荷电量Q<sub>cs</sub>或放电过程剩余可放荷电量Q<sub>ds</sub>小于采样周期T内风电预报误差修正要求对应荷电量时,以电池的荷电状态SOC的充放电保护极限对应的可用荷电量作为标准,提前调整修正量,通过输入输出电能限制调整充放电量;2)对电池生命周期进行预测:2.1)根据风电输出每次采样得到的周期T内的数据,如果<img file="FDA0000818549250000031.GIF" wi="235" he="120" />为正,修正误差输出的电能<img file="FDA0000818549250000032.GIF" wi="244" he="127" />累计入W<sub>outm</sub>;如果<img file="FDA0000818549250000033.GIF" wi="232" he="119" />为负,修正误差输入的电能累计入W<sub>inn</sub>;2.2)当满足|W<sub>outm</sub>‑0.5W<sub>out1</sub>|≤γ且W<sub>inm</sub>≥0.5W<sub>in1</sub>;或<img file="FDA0000818549250000034.GIF" wi="367" he="99" />且W<sub>outm</sub>≥0.5W<sub>out1</sub>时,符合条件的m或n值即为半个循环周期数;2.3)根据系统修正需求,输出总和与输入总和分别与基准值0.5W<sub>in1</sub>,0.5W<sub>out1</sub>相比,其中较小值就是目前完成半个循环的次数;3)电池健康状态的诊断:根据上述步骤所得的效率系数及相关修正参数,利用循环计量法,采用如下办法进行电池健康状态监测:用某时间段输入输出储能装置过程的电能及电池本身的荷电量变化绝对损耗或此过程的相对损耗来监测;或在计及电池自放电效率及充放电变换器效率情况下,对于同等的充放荷电量,监测输入输出电能的实际数量,由电能损耗在多次循环中的数值变化趋势,判断电池健康状态;或对于计及传输效率的同等的输入输出电能,监测充放电的实际荷电量变化,由荷电量自放电损耗在多次循环中的数值变化趋势,判断电池健康状态;或根据充放电传输效率与自放电效率变化趋势来判断电池健康状态;选定上述某一循环计量过程,对某一时间段t,当<img file="FDA0000818549250000041.GIF" wi="518" he="108" />时,利用此过程中W<sub>outt</sub>与W<sub>int</sub>关系修正<img file="FDA0000818549250000042.GIF" wi="258" he="150" />如果η<sub>bk</sub><30%,表明电池寿命已到。
地址 210098 江苏省南京市鼓楼区西康路1号