发明名称 一种微电网储能裕度的检测方法
摘要 本发明公开了一种微电网储能裕度检测方法,其特点是根据引入储能设备前后,微电网所能承载负荷的能力的变化,对微电网储能裕度进行检测或是计算微电网要达到某一储能裕度时,所需的储能装置的最小容量。不投入储能设备,建立基于成本目标、可再生能源波动抑制目标和功率需求匹配目标的模型,将非间歇性微电源作为决策变量带入,采用粒子群算法对模型求解,记录所求得的非间歇性微电源的出力和单位发电负荷所对应的综合成本。投入储能设备,将非间歇性微电源和储能设备出力作为决策变量带入模型并采用粒子群算法求解,记录此时所的储能设备的出力和单位发电负荷所对应的综合成本。采用二分法,不断调整投入储能设备后微电网的负荷水平,直到达到终止条件。
申请公布号 CN104392394A 申请公布日期 2015.03.04
申请号 CN201410673998.1 申请日期 2014.11.20
申请人 四川大学;国网上海市电力公司 发明人 刘天琪;胡晓通;李兴源;时珊珊;刘舒;柳劲松;张程翔;周宜广
分类号 G06Q50/06(2012.01)I 主分类号 G06Q50/06(2012.01)I
代理机构 成都科海专利事务有限责任公司 51202 代理人 邓继轩
主权项 一种微电网储能裕度的检测方法,其特征在于该方法包括以下步骤:1)储能裕度的检测目标函数本发明建立了含有风力发电机、光伏电池、微型燃气轮机等多种微电源、负荷及储能设备的微电网储能裕度检测模型;考虑了成本目标、可再生能源波动抑制目标、功率需求匹配目标在内的目标函数,考虑了微电网系统中的功率平衡、各个微电源容量限制、微型燃气轮机爬坡率、联络线传输功率、蓄电池容量和蓄电池充放电各类约束;minf=f<sub>1</sub>+f<sub>2</sub>+f<sub>3</sub>                    (1)其中<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msub><mi>C</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mi>DP</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mi>OM</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>+</mo><msub><mi>C</mi><mi>grid</mi></msub><mo>+</mo><msub><mi>C</mi><mi>punish</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000011.GIF" wi="1375" he="151" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>f</mi><mn>2</mn></msub><mo>=</mo><mi>E</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>1</mn></msub></mrow></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000012.GIF" wi="1252" he="77" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>f</mi><mn>3</mn></msub><mo>=</mo><mi>E</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>2</mn></msub></mrow></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000013.GIF" wi="1254" he="74" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>C</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>F</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000014.GIF" wi="1253" he="137" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>C</mi><mi>DP</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>[</mo><mfrac><msub><mi>C</mi><mrow><mi>az</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mn>8760</mn><msub><mi>k</mi><mi>i</mi></msub></mrow></mfrac><mo>&times;</mo><mfrac><mrow><mi>r</mi><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>r</mi><mo>)</mo></mrow><msub><mi>n</mi><mi>i</mi></msub></msup></mrow><mrow><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>r</mi><mo>)</mo></mrow><msub><mi>n</mi><mi>i</mi></msub></msup><mo>-</mo><mn>1</mn></mrow></mfrac><mo>&times;</mo><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000015.GIF" wi="1259" he="147" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>C</mi><mi>OM</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>[</mo><msub><mi>M</mi><mrow><mi>OM</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>&times;</mo><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000016.GIF" wi="1199" he="134" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>C</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mo>[</mo><msub><mi>V</mi><mi>ej</mi></msub><msub><mi>Q</mi><mi>ij</mi></msub><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000017.GIF" wi="1202" he="142" /></maths>C<sub>grid</sub>(t)=E(t)×P<sub>grid</sub>(t)            (9)<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>23</mn></munderover><msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>DG</mi><mo>,</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>DG</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000018.GIF" wi="1193" he="141" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><msub><mi>P</mi><mrow><mi>L</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>DG</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>|</mo></mrow><msub><mi>P</mi><mi>N</mi></msub></mfrac><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000019.GIF" wi="1198" he="167" /></maths>式中,f<sub>1</sub>、f<sub>2</sub>、f<sub>3</sub>分别为成本目标、可再生能源波动抑制目标、功率需求匹配目标,α<sub>1</sub>、α<sub>2</sub>分别为各个时刻分布式电源的波动之和与功率需求不匹配之和,C<sub>f</sub>(t)、C<sub>DP</sub>(t)、C<sub>OM</sub>(t)、C<sub>e</sub>(t)分别为t时刻各微电源的燃料成本、投资折旧成本、运行维护成本、排污处理成本;C<sub>grid</sub>、P<sub>grid</sub>分别为微电网与配电网的交互成本与交换功率;C<sub>punish</sub>为不满足负荷需求所引入的惩罚函数;E(t)为t时刻的实时购售电价;C<sub>az,i</sub>、k<sub>i</sub>、r、n<sub>i</sub>、K<sub>OM,i</sub>、V<sub>ej</sub>、Q<sub>ij</sub>分别为单位容量的微电源的安装成本、容量因素、年利率、投资偿还期、单位电量运行维护成本系数、污染物的环境价值和污染物的排放量;其约束条件为:A、功率平衡约束<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>P</mi><mi>load</mi></msub><mo>=</mo><msub><mi>P</mi><mi>grid</mi></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>P</mi><mrow><mi>DG</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000612826920000021.GIF" wi="1372" he="135" /></maths>式中,P<sub>load</sub>为系统总的负荷,P<sub>grid</sub>为微电网与配网的交换功率,P<sub>DG,i</sub>为各微电源的出力;B、微电源有功出力约束P<sub>i,min</sub>≤P<sub>i</sub>≤P<sub>i,max</sub>                      (13)C、微型燃气轮机爬坡率约束增加出力时为:P<sub>MT</sub>(t)‑P<sub>MT</sub>(t‑1)≤R<sub>up</sub>                   (14)减少出力时为:P<sub>MT</sub>(t‑1)‑P<sub>MT</sub>(t)≤R<sub>down</sub>                    (15)D、微电网与配电网交换功率约束P<sub>grid,min</sub>≤P≤P<sub>grid,max</sub>                    (16)E、蓄电池运行约束P<sub>SB,min</sub>≤P≤P<sub>SB,max</sub>                       (17)SOC<sub>min</sub>≤SOC≤SOC<sub>max</sub>                         (18)对微电网中的风、光等出力按照分布式电源模型进行求解,将求得的风、光出力以及负荷作为已知量代入粒子群算法中;将其它微电源作为未知量,进行编码操作,生成初始种群,分别求取出投入与未投入单位发电负荷所对应的综合成本,采用二分法,不断调整投入储能设备后微电网的负荷水平,直到达到终止条件(即投入储能设备前后单位发电负荷的综合成本之差的绝对值小于设定值),从而求得微电网某运行方式下的储能裕度,即在维持单位发电负荷所对应的综合成本不变的前提下,新增的可承载负荷的能力即为微电网的储能裕度;2)储能裕度的检测算法A、对微电网储能检测模型中的各微电源、负荷参数、目标函数参数以及约束条件参数进行设置,对粒子群算法中的控制参数进行设置,其中包括仿真代数、粒子规模、学习因子和惯性权重等;B、不投入储能设备,对微电网中的非间歇性微电源出力进行编码,随机生成初始种群,通过解码,将个体代入适应度函数中,对个体进行适应度值评价,选出种群中的最优个体;C、对种群进行位置和速度行进更新,同时在操作过程中采用最优保留政策,生成新的子代种群;D、新生成种群重新计算群体中各个体的适应度值,进入判别条件,满足条件的个体作为最优解输出,不满足条件则返回步骤C;E、输出的最优解即为未投入储能设备的一天24小时非间歇性微电源和单位发电负荷所对应的综合成本。F、投入储能设备,对微电网中的储能设备和非间歇性微电源的出力进行编码,随机生成初始种群,通过解码,将个体代入适应度函数中,对个体进行适应度值评价,选出种群中的最优个体;G、对种群进行位置和速度行进更新,同时在操作过程中采用最优保留政策,生成新的子代种群;H、新生成种群重新计算群体中各个体的适应度值,进入判别条件,满足条件的个体作为最优解输出,不满足条件则返回步骤G;I、输出的最优解即为投入储能设备的一天24小时非间歇性微电源和储能设备的出力和单位发电负荷所对应的综合成本;J、将投入储能装置前后的单位发电负荷所对应的综合成本作差;K、进入判别条件,满足条件转到步骤M,不满足条件则转到步骤L;L、采用二分法,不断调整投入储能设备后微电网的负荷水平;M、最终输出结构为即在维持系统评价指标不变的前提下,由于储能设备的投入,微电网新增的可承载负荷的能力,即为储能裕度。
地址 610065 四川省成都市武侯区一环路南一段24号