发明名称 一种微电网的能量优化管理方法
摘要 一种微电网的能量优化管理方法,该方法基于微网内分布式电源和储能元件的类别、个数以及储能元件的充电状态及电压百分比建立微网内约束矩阵,根据微网内各发电单元和储能单元的发电时刻、荷电状态等计算各发电单元储能单元的发电功率上下限,根据不同分布式能源特性确定相应惩罚因子,以便于优化这些系统运行;将优化条件作用于优化算法确定微网能量优化管理系统的目标函数。本发明充分考虑了分布式能源发电能力的动态变化,使分布式发电系统能量优化管理更有效,在能量优化管理时根据不同分布式能源特性提出相应惩罚因子,可以优化这些系统运行,延长系统寿命。
申请公布号 CN104978609A 申请公布日期 2015.10.14
申请号 CN201510366748.8 申请日期 2015.06.27
申请人 云南电网有限责任公司电力科学研究院 发明人 苏适;杨洋;严玉廷;李萍;杨家全;陆海;戴海琪
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 昆明大百科专利事务所 53106 代理人 何健
主权项 一种微电网能量优化管理方法,本发明特征是,包括以下步骤:步骤一:基于微网内分布式电源和储能元件的类别、个数以及储能元件的充电状态及电压百分比建立微网内约束矩阵如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>l</mi><mrow><mn>1</mn><mi>min</mi></mrow></msub></mtd><mtd><msub><mi>l</mi><mrow><mn>1</mn><mi>max</mi></mrow></msub></mtd><mtd><msub><mi>p</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>l</mi><mrow><mn>2</mn><mi>min</mi></mrow></msub></mtd><mtd><msub><mi>l</mi><mrow><mn>2</mn><mi>max</mi></mrow></msub></mtd><mtd><msub><mi>p</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>l</mi><mrow><mn>3</mn><mi>min</mi></mrow></msub></mtd><mtd><msub><mi>l</mi><mrow><mn>3</mn><mi>max</mi></mrow></msub></mtd><mtd><msub><mi>p</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>l</mi><mrow><mn>4</mn><mi>min</mi></mrow></msub></mtd><mtd><msub><mi>l</mi><mrow><mn>4</mn><mi>max</mi></mrow></msub></mtd><mtd><msub><mi>p</mi><mn>4</mn></msub></mtd></mtr><mtr><mtd><msub><mi>l</mi><mrow><mn>5</mn><mi>min</mi></mrow></msub></mtd><mtd><msub><mi>l</mi><mrow><mn>5</mn><mi>max</mi></mrow></msub></mtd><mtd><msub><mi>p</mi><mn>5</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>l</mi><mrow><mi>n</mi><mi>min</mi></mrow></msub></mtd><mtd><msub><mi>l</mi><mrow><mi>n</mi><mi>max</mi></mrow></msub></mtd><mtd><msub><mi>p</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000747126010000011.GIF" wi="392" he="391" /></maths>步骤二:根据微网内各发电单元和储能单元的发电时刻、荷电状态等计算各发电单元储能单元的发电功率上下限l<sub>min</sub>、l<sub>max</sub>;对于几种微网常见的发电、储能单元如光伏发电、风电、燃气轮机、蓄电池、超级电容的考虑方法:1)对于光伏发电,a.对于有光伏发电预测的微网能量优化管理系统,l<sub>max</sub>=P<sub>n</sub>,l<sub>min</sub>=0;b.对于没有光伏发电预测的微网能量优化管理系统:由于光伏与一天内的发电时刻有关,则当:<img file="FDA0000747126010000012.GIF" wi="1161" he="211" />2)对于风力发电l<sub>max</sub>=P<sub>n</sub> l<sub>min</sub>=0;3)对于燃气轮机l<sub>max</sub>=P<sub>n</sub> l<sub>min</sub>=0;4)对于蓄电池,基于其荷电状态确定其发电功率上下限:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mn>0</mn><mi>%</mi><mo>&le;</mo><msub><mi>S</mi><mrow><mi>o</mi><mi>c</mi></mrow></msub><mo>&lt;</mo><mn>20</mn><mi>%</mi></mrow></mtd><mtd><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>=</mo><mn>20</mn><msub><mi>%P</mi><mi>n</mi></msub><mo>,</mo><msub><mi>l</mi><mi>min</mi></msub><mo>=</mo><mo>-</mo><mn>20</mn><msub><mi>%P</mi><mi>n</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>20</mn><mi>%</mi><mo>&le;</mo><msub><mi>S</mi><mrow><mi>o</mi><mi>c</mi></mrow></msub><mo>&lt;</mo><mn>30</mn><mi>%</mi></mrow></mtd><mtd><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>=</mo><mn>50</mn><msub><mi>%P</mi><mi>n</mi></msub><mo>,</mo><msub><mi>l</mi><mi>min</mi></msub><mo>=</mo><mo>-</mo><mn>50</mn><msub><mi>%P</mi><mi>n</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>30</mn><mi>%</mi><mo>&le;</mo><msub><mi>S</mi><mrow><mi>o</mi><mi>c</mi></mrow></msub><mo>&le;</mo><mn>90</mn><mi>%</mi></mrow></mtd><mtd><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>=</mo><msub><mi>P</mi><mi>n</mi></msub><mo>,</mo><msub><mi>l</mi><mi>min</mi></msub><mo>=</mo><mo>-</mo><msub><mi>P</mi><mi>n</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>90</mn><mi>%</mi><mo>&lt;</mo><msub><mi>S</mi><mrow><mi>o</mi><mi>c</mi></mrow></msub><mo>&lt;</mo><mn>100</mn><mi>%</mi></mrow></mtd><mtd><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>=</mo><mn>20</mn><msub><mi>%P</mi><mi>n</mi></msub><mo>,</mo><msub><mi>l</mi><mi>min</mi></msub><mo>=</mo><mo>-</mo><mn>20</mn><msub><mi>%P</mi><mi>n</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000747126010000013.GIF" wi="1132" he="268" /></maths>5)对于超级电容,基于加在电容两端电压确定其发电功率上下限:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mrow><mn>0</mn><mi>%</mi><mo>&le;</mo><msub><mi>U</mi><mi>c</mi></msub><mo>&lt;</mo><mn>75</mn><msub><mi>%U</mi><mi>n</mi></msub></mrow></mtd><mtd><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>=</mo><mn>50</mn><msub><mi>%P</mi><mi>n</mi></msub><mo>,</mo><msub><mi>l</mi><mi>min</mi></msub><mo>=</mo><mo>-</mo><mn>50</mn><msub><mi>%P</mi><mi>n</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>75</mn><msub><mi>%U</mi><mi>n</mi></msub><mo>&le;</mo><msub><mi>U</mi><mi>c</mi></msub><mo>&lt;</mo><msub><mi>U</mi><mi>n</mi></msub></mrow></mtd><mtd><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>=</mo><msub><mi>P</mi><mi>n</mi></msub><mo>,</mo><msub><mi>l</mi><mi>min</mi></msub><mo>=</mo><mo>-</mo><msub><mi>P</mi><mi>n</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>U</mi><mi>n</mi></msub><mo>&le;</mo><msub><mi>U</mi><mi>c</mi></msub><mo>&le;</mo><msub><mi>U</mi><mi>max</mi></msub></mrow></mtd><mtd><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>=</mo><mn>50</mn><msub><mi>%P</mi><mi>n</mi></msub><mo>,</mo><msub><mi>l</mi><mi>min</mi></msub><mo>=</mo><mo>-</mo><mn>50</mn><msub><mi>%P</mi><mi>n</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000747126010000021.GIF" wi="1162" he="205" /></maths>步骤三:根据不同分布式能源特性确定相应惩罚因子,以便于优化这些系统运行;1)对于光伏发电和风力发电,根据其运行时间确定惩罚因子:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>p</mi><mo>=</mo><mfrac><msub><mi>t</mi><mi>u</mi></msub><mn>24</mn></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000747126010000022.GIF" wi="172" he="144" /></maths>式中为光伏或风力发电的运行时间;2)对于燃气轮机,其惩罚因子与发出的功率成正比,取p=b;式中b为一常数,取值:0∠b∠1;3)对于蓄电池,根据其荷电状态、发电功率上下限l<sub>min</sub>、及当前功率x确定其惩罚因子p:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>p</mi><mo>=</mo><mi>f</mi><mi>l</mi><mi>a</mi><mi>g</mi><mo>*</mo><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>l</mi><mi>min</mi></msub></mrow><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>-</mo><msub><mi>l</mi><mi>min</mi></msub></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000747126010000024.GIF" wi="455" he="155" /></maths>其中:flag为惩罚方向参数,当0%<S<sub>oc</sub><3时,flag=1;当30%<S<sub>oc</sub><10时,flag=‑1;x为当前所发出的功率;l<sub>min</sub>、与之前同;4)对于超级电容,根据其加在电容两端电压、发电功率上下限l<sub>min</sub>、及当前功率x确定其惩罚因子p:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>p</mi><mo>=</mo><mi>f</mi><mi>l</mi><mi>a</mi><mi>g</mi><mo>*</mo><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>l</mi><mi>min</mi></msub></mrow><mrow><msub><mi>l</mi><mi>max</mi></msub><mo>-</mo><msub><mi>l</mi><mi>min</mi></msub></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000747126010000023.GIF" wi="444" he="147" /></maths>其中:flag为惩罚方向参数,当0%≤U<sub>c</sub><75%时,flag=1;当75%U<sub>n</sub>≤U<sub>c</sub>≤U时,flag=‑1;x为当前所发出的功率;l<sub>min</sub>、与之前同;步骤四:将优化条件作用于优化算法确定微网能量优化管理系统的目标函数:1)目标函数a.孤网运行模型下:保证本地负载供电,并使微网总发电成本最小;<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>min</mi><mi>C</mi><mo>=</mo><mi>k</mi><mn>1</mn><msub><mi>C</mi><mi>G</mi></msub><mo>+</mo><mi>k</mi><mn>2</mn><msub><mi>C</mi><mrow><mi>O</mi><mi>M</mi></mrow></msub><mo>+</mo><mi>k</mi><mn>3</mn><msub><mi>C</mi><mrow><mi>D</mi><mi>P</mi></mrow></msub><mo>+</mo><mi>k</mi><mn>4</mn><msub><mi>C</mi><mi>L</mi></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msub><mi>x</mi><mi>i</mi></msub><msub><mi>p</mi><mi>i</mi></msub></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000747126010000031.GIF" wi="972" he="197" /></maths>其中:<img file="FDA0000747126010000032.GIF" wi="399" he="85" />为发电单元燃料消耗成本,i为发电单元编号,t为运行时刻;<img file="FDA0000747126010000033.GIF" wi="656" he="72" />发电单元运行管理成本,K<sub>oMi</sub>为运行管理系数;<img file="FDA0000747126010000034.GIF" wi="401" he="114" />为发电单元折旧成本,P<sub>fc,i</sub>为最大输出功率,C<sub>fi</sub>为容量因素;<img file="FDA0000747126010000035.GIF" wi="388" he="79" />为负荷停运补偿成本;k1、k2、k3、k4为费用考虑系数,取0或1,(0表示不考虑该项费用,1表示考虑该项费用);<img file="FDA0000747126010000036.GIF" wi="105" he="77" />为微网系统内对应的惩罚项,xi表示发电或储能单元当前吸收或者发出的功率,p<sub>i</sub>为相应惩罚因子;b.并网运行模式下:保证本地负载供电,并使微网经济效益最大;<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>min</mi><mi>C</mi><mo>=</mo><mi>k</mi><mn>1</mn><msub><mi>C</mi><mi>G</mi></msub><mo>+</mo><mi>k</mi><mn>2</mn><msub><mi>C</mi><mrow><mi>O</mi><mi>M</mi></mrow></msub><mo>+</mo><mi>k</mi><mn>3</mn><msub><mi>C</mi><mrow><mi>D</mi><mi>P</mi></mrow></msub><mo>+</mo><mi>k</mi><mn>4</mn><msub><mi>C</mi><mi>L</mi></msub><mo>+</mo><mi>k</mi><mn>5</mn><msub><mi>C</mi><mrow><mi>g</mi><mi>r</mi><mi>i</mi><mi>d</mi></mrow></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msub><mi>x</mi><mi>i</mi></msub><msub><mi>p</mi><mi>i</mi></msub></mrow></mrow>]]></math><img file="FDA0000747126010000037.GIF" wi="1150" he="167" /></maths>其中:C<sub>grid</sub>=C<sub>bt</sub>C<sub>bp</sub>‑C<sub>st</sub>C<sub>sp</sub>为微网与配电网交互成本,C<sub>bt</sub>,C<sub>st</sub>为向电网买入、卖出的电量,C<sub>bp</sub>,C<sub>sp</sub>为向电网买入、卖出的电价;C<sub>G</sub>、C<sub>OM</sub>、C<sub>grid</sub>、C<sub>DP</sub>、C<sub>L</sub>、k1、k2、k3、k4、k5与上同;2)约束条件a.功率平衡约束:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></munderover><msub><mi>P</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><mrow><msub><mi>P</mi><mrow><mi>S</mi><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>g</mi><mi>r</mi><mi>i</mi><mi>d</mi></mrow></msub></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>P</mi><mrow><mi>L</mi><mi>i</mi></mrow></msub></mrow>]]></math><img file="FDA0000747126010000038.GIF" wi="642" he="178" /></maths>其中:为系统中第i个分布式发电单元发出的功率,I为系统中发电单元的个数;P<sub>Sj</sub>为第j个储能单元吸收或发出的功率,J为系统中储能单元的个数;P<sub>Li</sub>为第i个负荷需求的功率,K为系统中负荷的个数;p<sub>grid</sub>为配电网的交换功率,当微网向电网输出功率时,P<sub>grid</sub>≥0;当微网从电网吸收功率时,P<sub>grid</sub><0;P<sub>Li</sub>为负荷单元所需功率;b.发电单元输出功率(发电和储能)<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msubsup><mi>P</mi><mrow><mi>G</mi><mi>i</mi></mrow><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>P</mi><mrow><mi>G</mi><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>G</mi><mi>i</mi></mrow><mi>max</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>P</mi><mrow><mi>S</mi><mi>j</mi></mrow><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>P</mi><mrow><mi>S</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>S</mi><mi>j</mi></mrow><mi>max</mi></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000747126010000041.GIF" wi="475" he="157" /></maths>式中<img file="FDA0000747126010000042.GIF" wi="514" he="75" />分别为系统中第i个分布式电源发电单元和第j个储能单元的输出功率下限和输出功率上限;c.微网与配电网交换容量约束<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>g</mi><mi>r</mi><mi>i</mi><mi>d</mi></mrow><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>P</mi><mrow><mi>g</mi><mi>r</mi><mi>i</mi><mi>d</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>g</mi><mi>r</mi><mi>i</mi><mi>d</mi></mrow><mi>max</mi></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000747126010000043.GIF" wi="484" he="76" /></maths>式中,<img file="FDA0000747126010000044.GIF" wi="258" he="86" />分别为微网系统与大电网能量交换的功率下限和功率上限。
地址 650217 云南省昆明市经济技术开发区云大西路105号