发明名称 一种面向系统调控需求的储能系统最优配置方法
摘要 本发明涉及一种面向系统调控需求的储能系统最优配置方法,其特点是:从风电对整个电网功率的不平衡性角度出发,调控风电场群输出功率不超出整个电网功率平衡目标的同时,在剩余的电网可利用空间下,设定储能系统放电区间,基于SOC分层控制策略,对储能系统能量进行优化管理,实时修正储能系统充放电功率,优化储能系统工作性能,并建立风储运行经济性评估模型,以综合收益最大为目标,确定最佳储能系统放电区间,优化储能系统容量,延长其使用寿命,提高新能源电量的入网规模,有效利用电网可接纳风电空间。
申请公布号 CN104538983B 申请公布日期 2016.08.17
申请号 CN201510038095.0 申请日期 2015.01.26
申请人 东北电力大学 发明人 李军徽;严干贵;葛延峰;王月;冯凯翔
分类号 H02J3/32(2006.01)I;H02J3/38(2006.01)I 主分类号 H02J3/32(2006.01)I
代理机构 吉林市达利专利事务所 22102 代理人 陈传林
主权项 一种面向系统调控需求的储能系统最优配置方法,其特征是,调控风电场群输出功率不超出整个电网功率平衡目标的同时,在剩余的电网可利用空间下,设定储能系统放电区间,基于SOC分层控制策略,对储能系统能量进行优化管理,实时修正储能系统充放电功率,优化储能系统工作性能,建立风储运行经济性评估模型,以综合收益最大为目标,确定最佳所需储能系统放电区间,优化储能系统容量,延长其使用寿命,提高新能源电量的入网规模,有效利用电网可接纳风电空间,它包括以下步骤:1)储能系统的放电区间设定在接纳风电功率后未突破电网可利用空间极限值的时段,设定储能系统的放电区间α,0≤α&lt;100%;若系统无剩余可利用空间时α=1,若储能系统不放电α=0;基于放电区间α的储能系统充放电功率如下:<maths num="0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>E</mi><mi>S</mi><mi>S</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mrow><mi>w</mi><mi>d</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>P</mi><mrow><mi>lim</mi><mi>i</mi><mi>t</mi></mrow><mrow><mi>s</mi><mi>p</mi><mi>a</mi><mi>c</mi><mi>e</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>P</mi><mrow><mi>w</mi><mi>d</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&gt;</mo><msubsup><mi>P</mi><mrow><mi>lim</mi><mi>i</mi><mi>t</mi></mrow><mrow><mi>s</mi><mi>p</mi><mi>a</mi><mi>c</mi><mi>e</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>E</mi><mi>S</mi><mi>S</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>&alpha;</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>w</mi><mi>d</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><msubsup><mi>P</mi><mrow><mi>lim</mi><mi>i</mi><mi>t</mi></mrow><mrow><mi>s</mi><mi>p</mi><mi>a</mi><mi>c</mi><mi>e</mi></mrow></msubsup><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msub><mi>P</mi><mrow><mi>w</mi><mi>d</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&lt;</mo><msubsup><mi>P</mi><mrow><mi>lim</mi><mi>i</mi><mi>t</mi></mrow><mrow><mi>s</mi><mi>p</mi><mi>a</mi><mi>c</mi><mi>e</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001006615350000011.GIF" wi="1774" he="195" /></maths>其中P<sub>ESS</sub>(t)为t时刻储能系统充放电功率;P<sub>wd</sub>(t)、<img file="FDA0001006615350000012.GIF" wi="161" he="63" />分别为t时刻风电场群实际输出功率之和以及风电可运行域极值;面向系统需求区间的储能系统充放电能量E<sub>t</sub>以及储能系统在各调度时段结束后充放电累积容量W<sub>t</sub>如下所示:<maths num="0002"><math><![CDATA[<mrow><msub><mi>E</mi><mi>t</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mo>&Integral;</mo><mrow><mi>t</mi><mn>1</mn></mrow><mrow><mi>t</mi><mn>2</mn></mrow></msubsup><msub><mi>P</mi><mrow><mi>E</mi><mi>S</mi><mi>S</mi></mrow></msub><msub><mi>&eta;</mi><mrow><mi>c</mi><mi>h</mi><mi>arg</mi><mi>e</mi></mrow></msub><mi>d</mi><mi>t</mi></mrow></mtd><mtd><mrow><msub><mi>P</mi><mrow><mi>E</mi><mi>S</mi><mi>S</mi></mrow></msub><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mo>&Integral;</mo><mrow><mi>t</mi><mn>1</mn></mrow><mrow><mi>t</mi><mn>2</mn></mrow></msubsup><msub><mi>P</mi><mrow><mi>E</mi><mi>S</mi><mi>S</mi></mrow></msub><mo>/</mo><msub><mi>&eta;</mi><mrow><mi>d</mi><mi>i</mi><mi>s</mi><mi>c</mi><mi>h</mi><mi>arg</mi><mi>e</mi></mrow></msub><mi>d</mi><mi>t</mi></mrow></mtd><mtd><mrow><msub><mi>P</mi><mrow><mi>E</mi><mi>S</mi><mi>S</mi></mrow></msub><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001006615350000013.GIF" wi="1582" he="262" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>W</mi><mi>t</mi></msub><mo>=</mo><msub><mi>E</mi><mn>0</mn></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>t</mi></munderover><msub><mi>E</mi><mi>t</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001006615350000014.GIF" wi="1366" he="166" /></maths>其中t<sub>1</sub>,t<sub>2</sub>分别为充放电的起始与结束时刻;η<sub>charge</sub>,η<sub>discharge</sub>分别为储能系统的充放电效率;P<sub>ESS</sub>为储能系统充放电功率;E<sub>0</sub>为储能系统初始能量;风电场群所需配置储能系统容量W<sub>opt</sub>如下所示:W<sub>opt</sub>=max{W<sub>t</sub>}‑min{W<sub>t</sub>},t=1,2,…,N        (4)其中max{W<sub>t</sub>}为各调度时段累计容量最大值,min{W<sub>t</sub>}为各调度时段累计容量最小值,N为储能系统总控制时段;2)SOC分层控制策略构建SOC分层控制策略,对储能系统充放电功率进行实时修正,确保储能系统具有良好工作性能;将储能系统SOC按照充放电能力分为以下五个层次:不充电紧急层、少充电预防层、正常充放电安全层、少放电预防层、不放电紧急层;储能系统充放电功率P<sub>ESS</sub>,经储能能量管理系统确定的修正系数K<sub>SOC</sub>进行动态调整,得到储能系统实际充放电指令P<sub>SOC_ESS</sub>;K<sub>SOC</sub>值与Sigmoid函数特性类似,因此利用Sigmoid函数对其进行修正,具体表达如下所示:储能系统处于充电状态下,P<sub>ESS</sub>(t)&gt;0<maths num="0004"><math><![CDATA[<mrow><msub><mi>K</mi><mrow><mi>S</mi><mi>O</mi><mi>C</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>S</mi><mi>max</mi></msub><mo>&le;</mo><mi>S</mi><mo>&le;</mo><mn>100</mn><mi>%</mi></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><mn>10</mn><mrow><mo>(</mo><msub><mi>x</mi><mi>c</mi></msub><mo>-</mo><mn>0.5</mn><mo>)</mo></mrow></mrow></msup></mrow></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>S</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi><mo>_</mo><mi>max</mi></mrow></msub><mo>&lt;</mo><mi>S</mi><mo>&lt;</mo><msub><mi>S</mi><mi>max</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mn>0</mn><mo>&le;</mo><mi>S</mi><mo>&le;</mo><msub><mi>S</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi><mo>_</mo><mi>max</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001006615350000021.GIF" wi="1605" he="302" /></maths>x<sub>c</sub>=(S‑S<sub>max</sub>)/(S<sub>pre_max</sub>‑S<sub>max</sub>)                         (6)储能系统处于放电状态下,P<sub>ESS</sub>(t)&lt;0<maths num="0005"><math><![CDATA[<mrow><msub><mi>K</mi><mrow><mi>S</mi><mi>O</mi><mi>C</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mn>0</mn><mi>%</mi><mo>&le;</mo><mi>S</mi><mo>&le;</mo><msub><mi>S</mi><mi>min</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><mn>10</mn><mrow><mo>(</mo><msub><mi>x</mi><mi>f</mi></msub><mo>-</mo><mn>0.5</mn><mo>)</mo></mrow></mrow></msup></mrow></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>S</mi><mi>min</mi></msub><mo>&lt;</mo><mi>S</mi><mo>&lt;</mo><msub><mi>S</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi><mo>_</mo><mi>min</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>S</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi><mo>_</mo><mi>min</mi></mrow></msub><mo>&le;</mo><mi>S</mi><mo>&le;</mo><mn>100</mn><mi>%</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001006615350000022.GIF" wi="1606" he="336" /></maths>x<sub>f</sub>=(S‑S<sub>min</sub>)/(S<sub>pre_min</sub>‑S<sub>min</sub>)                           (8)经调整系数K<sub>SOC</sub>修正确定储能系统实际充放电功率P<sub>SOC_ESS</sub>(t)为:P<sub>SOC_ESS</sub>(t)=K<sub>SOC</sub>P<sub>ESS</sub>(t)                             (9)其中S为储能系统的荷电状态;S<sub>max</sub>为不充电紧急层的下限;S<sub>max</sub>、S<sub>pre_max</sub>为少充电预防层的上下限;S<sub>pre_max</sub>、S<sub>pre_min</sub>为正常充放电安全层的上下限;S<sub>min</sub>为少放电预防层的下限;X<sub>c</sub>为储能系统充电状态下计算K<sub>SOC</sub>的系数;X<sub>f</sub>为储能系统放电状态下计算K<sub>SOC</sub>的系数;3)风储系统经济性评估模型评估模型以风储运行最终收益最大为目标,其目标函数表示如下:max(R<sub>income</sub>‑T<sub>invest</sub>)                               (10)其中R<sub>income</sub>为储能系统总收益费用,T<sub>invest</sub>为储能系统投资费用;模型的投资主要由储能系统容量投资与其运行维护费用两部分组成,表示为:T<sub>invest</sub>=W<sub>opt</sub>P+W<sub>opt</sub>Mn                            (11)其中,P为储能系统的单位容量价格;M为储能系统年运行维护费用;n为储能系统的运行年限;模型的收益主要分为配置储能系统多接纳的电量收益,储能系统存储能量改变量收益以及减少化石燃料发电的气体排放量所带来的环境保护收益三部分;采用储能系统面向系统需求的控制策略,多接纳电量收益R<sub>sw</sub>为:R<sub>sw</sub>=E<sub>sw</sub>I<sub>sw</sub>                               (12)<maths num="0006"><math><![CDATA[<mrow><msub><mi>E</mi><mrow><mi>s</mi><mi>w</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>365</mn><mi>n</mi></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>t</mi></munderover><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>s</mi><mi>w</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mi>w</mi></msub><mo>)</mo></mrow><mi>&Delta;</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001006615350000031.GIF" wi="1462" he="167" /></maths>其中E<sub>sw</sub>为利用储能装置面向系统需求控制后电网多接纳的风电电量;P<sub>sw</sub>为利用储能装置面向系统需求调控后风场实际发电功率;I<sub>sw</sub>为风电上网电价;P<sub>w</sub>为未安装储能装置风电场输出功率;t为每个调度日内的调度时段,Δt为调度时间,储能系统在调度日内累计的存储电量的变化量E<sub>s</sub>,其带来存储电量收益R<sub>s</sub>:R<sub>s</sub>=E<sub>s</sub>I<sub>s</sub>                                 (14)<maths num="0007"><math><![CDATA[<mrow><msub><mi>E</mi><mi>s</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>365</mn><mi>n</mi></mrow></munderover><mrow><mo>(</mo><msub><mi>E</mi><mi>T</mi></msub><mo>-</mo><msub><mi>E</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001006615350000032.GIF" wi="1406" he="174" /></maths>其中I<sub>s</sub>为储能电量上网电价,E<sub>T</sub>表示储能系统1个调度日内剩余电量;风储电量的增加促进新能源开发利用,减少化石燃料用量与减排费用,为环境保护带来收益;具体环境收益R<sub>CO2</sub>为:<maths num="0008"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><msub><mi>CO</mi><mn>2</mn></msub></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>E</mi><mrow><mi>s</mi><mi>w</mi></mrow></msub><mo>+</mo><msub><mi>E</mi><mi>s</mi></msub><mo>)</mo></mrow><msub><mi>I</mi><mrow><msub><mi>CO</mi><mn>2</mn></msub></mrow></msub><msub><mi>m</mi><mrow><msub><mi>CO</mi><mn>2</mn></msub></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001006615350000033.GIF" wi="1574" he="103" /></maths>其中I<sub>CO2</sub>为传统发电机组产生单位电能所排放CO<sub>2</sub>的减排价格,m<sub>co2</sub>为1kW·h的电能向大气排放的CO<sub>2</sub>的质量。
地址 132012 吉林省吉林市船营区长春路169号