发明名称 一种基于动态电价的微电网运行优化方法
摘要 本发明公开了一种基于动态电价的微电网运行优化方法,其特点是将配电网对微电网下达的参考计划交换功率曲线与微电网的计划交换功率曲线进行模糊化,计算二者的欧几里得贴近度来实现电价的动态化,建立考虑动态电价、运行维护成本和排污处理成本目标函数,采用粒子群算法求取微电网经济优化运行方案。使得微电网在运行优化过程中能够协调配电网运行,减轻微电网接入对配电网带来的影响。
申请公布号 CN104065072A 申请公布日期 2014.09.24
申请号 CN201410270976.0 申请日期 2014.06.16
申请人 四川大学;国网上海市电力公司 发明人 刘天琪;胡晓通;李兴源;时珊珊;李妮;刘舒;柳劲松;周宜广
分类号 H02J3/06(2006.01)I 主分类号 H02J3/06(2006.01)I
代理机构 成都科海专利事务有限责任公司 51202 代理人 邓继轩
主权项 一种基于动态电价的微电网运行优化方法,其特征在于:1)考虑微电网参与配电网协调运行优化模型将微电网当作一个可正可负的负荷接入配电网中,配电网根据运行优化目标确定出一条24h的微电网参考计划交换功率控制曲线,将参考计划交换功率和计划交换功率模糊化,分别求出隶属度函数值,通过计算二者的欧几里得贴近度来实现电价的动态化,曲线贴近度越高,售电价越高,购电价越低;2)计及动态电价的微电网运行优化目标函数本发明建立了含有风力发电机、光伏电池、微型燃气轮机、蓄电池多种微电源的微电网运行优化模型;考虑了包括动态电价、运行维护成本、排污处理成本在内的目标函数,考虑了微电网系统中的功率平衡、各个微电源容量限制、微型燃气轮机爬坡率、联络线传输功率、蓄电池容量和蓄电池充放电各类约束;同时模型通过计算参考计划交换功率曲线和计划交换功率曲线的欧几里得贴近度来确定动态电价,来计及配电网对微电网的运行;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>min</mi><mi> C</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><mo>[</mo><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>sc</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000521738690000011.GIF" wi="1446" he="141" /></maths>其中<maths num="0002" id="cmaths0002"><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>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000521738690000012.GIF" wi="1471" he="134" /></maths><maths num="0003" id="cmaths0003"><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>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000521738690000013.GIF" wi="1479" he="145" /></maths><maths num="0004" id="cmaths0004"><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>K</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>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000521738690000014.GIF" wi="1500" he="134" /></maths><maths num="0005" id="cmaths0005"><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>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000521738690000015.GIF" wi="1500" he="140" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>N</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><msqrt><mn>24</mn></msqrt></mfrac><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msup><mrow><mo>(</mo><mi>A</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mi>B</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000521738690000016.GIF" wi="1504" he="142" /></maths>C'<sub>P</sub>(t)=C<sub>P</sub>(t)×(0.5+N(A,B))                     (7)C'<sub>S</sub>(t)=C<sub>S</sub>(t)×(1.5+N(A,B))                     (8)C<sub>grid</sub>(t)=C'<sub>P</sub>(t)C<sub>GP</sub>(t)‑C'<sub>S</sub>(t)C<sub>SP</sub>(t)               (9)C<sub>sc</sub>(t)=Q<sub>ce</sub>(i)K<sub>pc</sub>                                (10)式中,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>、C<sub>sc</sub>分别为微电网与配电网的交互成本和微型燃气轮机制冷收益,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(i)、B(i)分别为参考计划交换功率曲线和计划交换功率曲线模糊化后的隶属度函数值;C<sub>P</sub>(t)、C<sub>P</sub>(t)为t时刻的实时购售电价,C'<sub>P</sub>(t)、C<sub>s</sub>'(t)为引入欧几里得贴近度后的动态实时购售电价;Q<sub>ce</sub>(i)为微型燃气轮机各时刻制冷量;K<sub>pc</sub>为单位制冷的售价;其约束条件为:A、功率平衡约束<maths num="0007" id="cmaths0007"><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>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000521738690000021.GIF" wi="1372" he="134" /></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>                         (12)C、微型燃气轮机爬坡率约束增加出力时为:P<sub>MT</sub>(t)‑P<sub>MT</sub>(t‑1)≤R<sub>up</sub>                       (13)减少出力时为:P<sub>MT</sub>(t‑1)‑P<sub>MT</sub>(t)≤R<sub>down</sub>                     (14)D、微网与配电网交换功率约束P<sub>grid,min</sub>≤P≤P<sub>grid,max</sub>                     (15)E、蓄电池运行约束P<sub>SB,min</sub>≤P≤P<sub>SB,max</sub>                         (16)SOC<sub>min</sub>≤SOC≤SOC<sub>max</sub>                        (17)对微电网中的风、光等出力按照分布式电源模型进行求解,将求得的风、光出力以及负荷作为已知量代入粒子群算法中;将其它微电源作为未知量,进行编码操作,生成初始种群,可以求取出一条计划交换功率曲线;3)计及动态电价的粒子群优化算法A、对微电网优化调度模型中的各微电源、负荷参数、目标函数参数以及约束条件参数进行设置,电价通过计算参考计划交换功率曲线和计划交换功率曲线的欧几里得贴近度来实现动态化,即曲线贴近度越高,售电价越高,购电价越低;B、对粒子群算法中的控制参数进行设置,其中包括仿真代数、粒子规模、学习因子和惯性权重等;C、对微电网中的燃气轮机、燃料电池和交换功率进行编码,随机生成初始种群;D、通过解码,将个体代入适应度函数中,对个体进行适应度值评价,选出种群中的最优个体;E、对种群进行位置和速度行进更新,同时在操作过程中采用最优保留政策,生成新的子代种群;F、新生成种群重新计算群体中各个体的适应度值,进入判别条件,满足条件的个体作为最优解输出,不满足条件则返回步骤4);G、最终输出的最优解即为一天24小时各微电源出力以及交换功率。
地址 610065 四川省成都市武侯区一环路南一段24号