发明名称 一种为电力系统实现电力平衡的电储能设备建模方法
摘要 本发明涉及一种为电力系统实现电力平衡的电储能设备建模方法,属于电力系统运行和控制技术领域。首先对电储能设备的充电、放电功率和充电过程与放电过程不能同时进行的特点进行数学建模。然后对电储能设备为电力系统提供的上旋转备用、下旋转备用与电储能设备充放电功率之间的关系、电储能设备上下旋转备用与电储能设备为电力系统提供的上旋转备用、下旋转备用容量之间的关系进行数学建模,最终得到为电力系统提供旋转备用的电储能设备调度模型。本方法解决了电储能设备应用于考虑旋转备用的电力系统调度的建模问题,为电储能设备应用于考虑旋转备用的电力系统调度问题提供模型基础。
申请公布号 CN105762821A 申请公布日期 2016.07.13
申请号 CN201610210652.7 申请日期 2016.04.06
申请人 清华大学 发明人 孙宏斌;郭庆来;王彬;张伯明;吴文传;晋宏杨
分类号 H02J3/28(2006.01)I;G06Q50/06(2012.01)I 主分类号 H02J3/28(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 一种为电力系统实现电力平衡的电储能设备建模方法,其特征在于,该方法包括以下步骤:(1)分别对电力系统电储能设备的充放电功率以及充电过程与放电过程不能同时进行的逻辑关系进行建模,具体过程如下:(1‑1)设电储能设备在时刻t的状态为<img file="FDA0000959160840000011.GIF" wi="91" he="63" />t=1,2,....,T,T为设定的电力系统调度时间尺度,若电储能设备在时刻t处于充电状态,则<img file="FDA0000959160840000012.GIF" wi="163" he="63" />若电储能设备在时刻t处于放电状态,则<img file="FDA0000959160840000013.GIF" wi="171" he="70" />(1‑2)设模型中电储能设备在时刻t的充电功率为<img file="FDA0000959160840000014.GIF" wi="126" he="63" />最大充电功率为<img file="FDA0000959160840000015.GIF" wi="118" he="63" />t=1,2,....,T,则电储能设备的充电功率的范围表述为调度模型如下:<maths num="0001"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>&CenterDot;</mo><msubsup><mi>x</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000959160840000016.GIF" wi="406" he="71" /></maths>(1‑3)设模型中电储能设备在时刻t的放电功率为<img file="FDA0000959160840000017.GIF" wi="126" he="62" />最大充电功率为<img file="FDA0000959160840000018.GIF" wi="130" he="63" />t=1,2,....,T,则电储能设备的放电功率的范围表述为调度模型如下:<maths num="0002"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>x</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000959160840000019.GIF" wi="502" he="63" /></maths>(1‑4)设模型中电储能设备在时刻t的存储能量值为<img file="FDA00009591608400000110.GIF" wi="163" he="63" />t=1,2,....,T,T为电力系统的调度时间尺度,则电储能设备的存储能量值<img file="FDA00009591608400000111.GIF" wi="141" he="63" />与步骤(1‑2)中电储能设备充电功率<img file="FDA00009591608400000112.GIF" wi="99" he="63" />和步骤(1‑3)中放电功率<img file="FDA00009591608400000113.GIF" wi="101" he="70" />的关系表述为模型如下:<maths num="0003"><math><![CDATA[<mrow><mo>-</mo><msub><mi>&eta;</mi><mi>c</mi></msub><mo>&CenterDot;</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>+</mo><mfrac><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><msub><mi>&eta;</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub></mfrac><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>+</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>SOC</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA00009591608400000114.GIF" wi="1054" he="135" /></maths>其中,η<sub>c</sub>为电储能设备充电损耗,η<sub>dc</sub>为电储能设备放电损耗,电储能设备自放电系数为γ,由电储能设备厂商给出,Δt为两个相邻调度时刻之间的时间间隔,取值为1;(2)对电力系统中电储能设备为电力系统提供的上旋转备用、下旋转备用及电储能设备为电力系统提供的上旋转备用、下旋转备用容量进行建模,具体步骤如下:(2‑1)设电力系统中电储能设备向电力系统提供的上旋转备用为<img file="FDA00009591608400000115.GIF" wi="116" he="62" />t=1,2,....,T,T为设定的电力系统的调度时间尺度,根据上述步骤(1)中电储能设备状态<img file="FDA00009591608400000116.GIF" wi="93" he="70" />电储能设备充电功率<img file="FDA00009591608400000117.GIF" wi="99" he="63" />和电储能设备放电功率<img file="FDA00009591608400000118.GIF" wi="129" he="63" />电储能设备向电力系统提供的上旋转备用<img file="FDA00009591608400000119.GIF" wi="88" he="63" />的范围为:<maths num="0004"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>Rh</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&le;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>x</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>)</mo></mrow><msubsup><mi>P</mi><mi>max</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000959160840000021.GIF" wi="683" he="63" /></maths>(2‑2)设电力系统中电储能设备向电力系统提供的下旋转备用为<img file="FDA0000959160840000022.GIF" wi="107" he="62" />t=1,2,....,T,T为电力系统的调度时间尺度,由人工事前设定,根据上述步骤(1)中的电储能设备状态<img file="FDA0000959160840000023.GIF" wi="94" he="63" />电储能设备充电功率<img file="FDA0000959160840000024.GIF" wi="96" he="62" />和电储能设备放电功率<img file="FDA0000959160840000025.GIF" wi="129" he="70" />电储能设备向电力系统提供的下旋转备用<img file="FDA0000959160840000026.GIF" wi="91" he="62" />的范围为:<maths num="0005"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>Rl</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>x</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><msubsup><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000959160840000027.GIF" wi="574" he="62" /></maths>(2‑3)设电力系统中电储能设备用于向电力系统提供上旋转备用的容量为<img file="FDA0000959160840000028.GIF" wi="187" he="61" />t=1,2,....,T,T为设定的电力系统的调度时间尺度,根据步骤(2‑1)中电储能设备向电力系统提供的上旋转备用<img file="FDA00009591608400000221.GIF" wi="91" he="63" />和上述步骤(1‑4)中所述的电储能设备放电损耗η<sub>dc</sub>,则电储能设备用于向电力系统提供下旋转备用的容量<img file="FDA0000959160840000029.GIF" wi="162" he="62" />范围表述为调度模型如下:<maths num="0006"><math><![CDATA[<mrow><mfrac><mrow><msubsup><mi>Rh</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup></mrow><msub><mi>&eta;</mi><mi>d</mi></msub></mfrac><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>+</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>h</mi></mrow></msubsup><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>SOC</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>h</mi></mrow></msubsup><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA00009591608400000210.GIF" wi="782" he="127" /></maths>其中,γ为电储能设备自放电系数,由电储能设备厂商给出,Δt为相邻两个调度时刻之间的时间间隔,取值为1;(2‑4)设电力系统中电储能设备用于向电力系统提供下旋转备用的容量为<img file="FDA00009591608400000211.GIF" wi="186" he="62" />t=1,2,....,T,T为设定的电力系统的调度时间尺度,则根据步骤(2‑2)中所述的电储能设备向电力系统提供的下旋转备用<img file="FDA00009591608400000220.GIF" wi="83" he="63" />和步骤(1‑4)中所述的电储能设备充电损耗η<sub>c</sub>,则电储能设备用于向电力系统提供下旋转备用的容量<img file="FDA00009591608400000212.GIF" wi="161" he="63" />范围表述为调度模型如下:<maths num="0007"><math><![CDATA[<mrow><mo>-</mo><msub><mi>&eta;</mi><mi>c</mi></msub><mo>&CenterDot;</mo><msubsup><mi>Rl</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>+</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>l</mi></mrow></msubsup><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>SOC</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>l</mi></mrow></msubsup><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA00009591608400000213.GIF" wi="869" he="63" /></maths>(2‑5)设电储能设备内部存储能量允许的最小能量值为<img file="FDA00009591608400000214.GIF" wi="174" he="70" />电储能设备内部存储能量允许的最大能量值为<img file="FDA00009591608400000215.GIF" wi="178" he="69" />将步骤(1‑4)中电储能设备在时刻t时内部存储能量值<img file="FDA00009591608400000216.GIF" wi="169" he="67" />步骤(2‑2)中电储能设备用于向电力系统提供上旋转备用的容量<img file="FDA00009591608400000217.GIF" wi="163" he="61" />以及步骤(2‑3)中电储能设备用于向电力系统提供下旋转备用的容量<img file="FDA00009591608400000218.GIF" wi="164" he="62" />三者之间的关系表述为调度模型如下:<maths num="0008"><math><![CDATA[<mrow><msubsup><mi>SOC</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>SOC</mi><mi>max</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>-</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>l</mi></mrow></msubsup><mo>-</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>h</mi></mrow></msubsup></mrow>]]></math><img file="FDA00009591608400000219.GIF" wi="902" he="63" /></maths>(3)综合上述步骤(1)所述的电储能设备的充电功率、放电功率和充电过程与放电过程不能同时进行的调度模型和步骤(2)所述的电储能设备为电力系统提供的上旋转备用、下旋转备用及电储能设备为电力系统提供的上旋转备用、下旋转备用容量调度模型,得到为电力系统提供旋转备用的电储能设备调度模型如下:<maths num="0009"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>&CenterDot;</mo><msubsup><mi>x</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000959160840000031.GIF" wi="405" he="63" /></maths><maths num="0010"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>x</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000959160840000032.GIF" wi="502" he="63" /></maths><maths num="0011"><math><![CDATA[<mrow><mo>-</mo><msub><mi>&eta;</mi><mi>c</mi></msub><mo>&CenterDot;</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>+</mo><mfrac><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><msub><mi>&eta;</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub></mfrac><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>+</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>SOC</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0000959160840000033.GIF" wi="1052" he="143" /></maths><maths num="0012"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>Rh</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&le;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>x</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>)</mo></mrow><msubsup><mi>P</mi><mi>max</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000959160840000034.GIF" wi="686" he="63" /></maths><maths num="0013"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>Rl</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>x</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><msubsup><mi>P</mi><mi>max</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>c</mi></mrow></msubsup><mo>+</mo><msubsup><mi>P</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>d</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000959160840000035.GIF" wi="581" he="63" /></maths><maths num="0014"><math><![CDATA[<mrow><mfrac><mrow><msubsup><mi>Rh</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup></mrow><msub><mi>&eta;</mi><mi>d</mi></msub></mfrac><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>+</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>h</mi></mrow></msubsup><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>SOC</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>h</mi></mrow></msubsup><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0000959160840000036.GIF" wi="782" he="135" /></maths><maths num="0015"><math><![CDATA[<mrow><mo>-</mo><msub><mi>&eta;</mi><mi>c</mi></msub><mo>&CenterDot;</mo><msubsup><mi>Rl</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi><mo>+</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>l</mi></mrow></msubsup><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>SOC</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>l</mi></mrow></msubsup><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0000959160840000037.GIF" wi="846" he="63" /></maths><maths num="0016"><math><![CDATA[<mrow><msubsup><mi>SOC</mi><mi>min</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>SOC</mi><mi>max</mi><mrow><mi>E</mi><mi>S</mi></mrow></msubsup><mo>-</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>l</mi></mrow></msubsup><mo>-</mo><msubsup><mi>SOC</mi><mi>t</mi><mrow><mi>E</mi><mi>S</mi><mo>,</mo><mi>r</mi><mi>e</mi><mi>h</mi></mrow></msubsup><mo>.</mo></mrow>]]></math><img file="FDA0000959160840000038.GIF" wi="933" he="79" /></maths>
地址 100084 北京市海淀区清华园1号