发明名称 一种电网中新能源容量配比分层优化方法
摘要 本发明公开了一种电网中新能源容量配比分层优化方法。本发明分为内外两层进行迭代计算,计算模型如下:内层在考虑该地区新能源特性的基础上对其出力进行时间序列建模,以电网节能减排效益最佳为目标,综合考虑负荷特性、机组调峰特性、不同种类供热机组热电耦合特性等因素,建立了计及新能源发电的年度时序生产模拟仿真模型;外层为容量配比优化模型,以内层模型的节能减排效益为适应度函数更新个体寻优方向,确定新能源发电配比容量,减少了随机生成新能源装机容量的盲目性,提高了优化效率和精度。本发明可用于省级电网新能源容量优化,对省(区)电网的新能源装机容量规划、低碳电力要求下的网源规划和实际电力系统调度都具有重要的指导意义。
申请公布号 CN103762589B 申请公布日期 2015.10.21
申请号 CN201410007489.5 申请日期 2014.01.08
申请人 河海大学 发明人 曹阳;袁越;孙承晨;郭思琪
分类号 H02J3/00(2006.01)I;H02J3/38(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 杨楠
主权项 一种电网中新能源容量配比分层优化方法,所述电网中的电能来源包括火电以及至少一种新能源,其特征在于,所述优化方法包括以下步骤:步骤1、初始化电网中各类新能源的装机容量,作为外层优化模型的初始值;步骤2、将外层优化模型得到的当前各类新能源的装机容量输入内层优化模型;内层优化模型通过求解以下模型,得到在当前各类新能源的装机容量下,使得二氧化碳排放量F最小的火电机组的全年运行状态:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>min</mi><mi>F</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>j</mi></msub></munderover><mrow><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mrow><mo>{</mo><msub><mi>&alpha;</mi><mi>j</mi></msub><msubsup><mi>Y</mi><mi>j</mi><mi>t</mi></msubsup><mo>+</mo><msub><mi>&beta;</mi><mi>j</mi></msub><msubsup><mi>Z</mi><mi>j</mi><mi>t</mi></msubsup><mo>+</mo><msub><mi>a</mi><mi>j</mi></msub><msubsup><mi>P</mi><mi>j</mi><mi>t</mi></msubsup><mo>+</mo><msub><mi>b</mi><mi>j</mi></msub><mo>}</mo><mo>&CenterDot;</mo><mi>&gamma;</mi></mrow></mrow></mrow>]]></math><img file="FDA0000771904230000011.GIF" wi="914" he="165" /></maths>受约束于:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>j</mi><mi>t</mi></msubsup><mo>&le;</mo><msubsup><mi>&Delta;P</mi><mi>j</mi><mrow><mi>u</mi><mi>p</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000771904230000012.GIF" wi="367" he="93" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>j</mi><mi>t</mi></msubsup><mo>-</mo><msubsup><mi>P</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>&le;</mo><msubsup><mi>&Delta;P</mi><mi>j</mi><mrow><mi>d</mi><mi>o</mi><mi>w</mi><mi>n</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000771904230000013.GIF" wi="387" he="86" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>X</mi><mi>j</mi><mi>t</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>P</mi><mi>j</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mi>j</mi><mi>t</mi></msubsup><mo>&le;</mo><msubsup><mi>X</mi><mi>j</mi><mi>t</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>P</mi><mi>j</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000771904230000014.GIF" wi="558" he="88" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>Y</mi><mi>j</mi><mi>t</mi></msubsup><mo>+</mo><msubsup><mi>Z</mi><mi>j</mi><mi>t</mi></msubsup><mo>&le;</mo><mn>1</mn></mrow>]]></math><img file="FDA0000771904230000015.GIF" wi="348" he="90" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>Y</mi><mi>j</mi><mi>t</mi></msubsup><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><msubsup><mi>Z</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mi>i</mi></mrow></msubsup><mo>&le;</mo><mn>1</mn></mrow></mrow>]]></math><img file="FDA0000771904230000016.GIF" wi="348" he="159" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>Z</mi><mi>j</mi><mi>t</mi></msubsup><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><msubsup><mi>Y</mi><mi>j</mi><mrow><mi>t</mi><mo>+</mo><mi>i</mi></mrow></msubsup><mo>&le;</mo><mn>1</mn></mrow></mrow>]]></math><img file="FDA0000771904230000017.GIF" wi="363" he="159" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>B</mi><mi>Y</mi></mrow><mi>t</mi></msubsup><mo>=</mo><msubsup><mi>C</mi><mi>j</mi><mi>b</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>H</mi><mi>j</mi><mi>t</mi></msubsup></mrow>]]></math><img file="FDA0000771904230000018.GIF" wi="343" he="107" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msubsup><mi>H</mi><mi>j</mi><mi>t</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>C</mi><mi>j</mi><mi>b</mi></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mrow><mi>C</mi><mi>Q</mi></mrow><mi>t</mi></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mi>j</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup><mo>-</mo><msubsup><mi>H</mi><mi>j</mi><mi>t</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>C</mi><mi>j</mi><mi>v</mi></msubsup></mrow>]]></math><img file="FDA0000771904230000019.GIF" wi="639" he="93" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>j</mi></msub></munderover><msubsup><mi>P</mi><mi>j</mi><mi>t</mi></msubsup><mo>+</mo><msubsup><mi>P</mi><mn>0</mn><mi>t</mi></msubsup><mo>=</mo><msubsup><mi>P</mi><mn>1</mn><mi>t</mi></msubsup></mrow>]]></math><img file="FDA00007719042300000110.GIF" wi="336" he="169" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>j</mi></msub></munderover><msubsup><mi>p</mi><mi>j</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup><mo>-</mo><msubsup><mi>C</mi><mi>N</mi><mi>t</mi></msubsup><mo>&le;</mo><mo>-</mo><msubsup><mi>P</mi><mn>1</mn><mi>t</mi></msubsup><mo>-</mo><msub><mi>S</mi><mi>p</mi></msub></mrow>]]></math><img file="FDA00007719042300000111.GIF" wi="501" he="149" /></maths><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>j</mi></msub></munderover><msubsup><mi>P</mi><mi>j</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msubsup><mo>+</mo><msubsup><mi>C</mi><mi>N</mi><mi>t</mi></msubsup><mo>&le;</mo><msubsup><mi>P</mi><mn>1</mn><mi>t</mi></msubsup><mo>-</mo><msub><mi>S</mi><mi>N</mi></msub></mrow>]]></math><img file="FDA00007719042300000112.GIF" wi="478" he="154" /></maths><maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>P</mi><mi>i</mi><mi>t</mi></msubsup><mo>&le;</mo><msub><mi>N</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msubsup><mi>P</mi><mi>i</mi><mrow><mi>t</mi><mo>*</mo></mrow></msubsup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>C</mi><mi>L</mi></mrow>]]></math><img file="FDA00007719042300000113.GIF" wi="764" he="84" /></maths>其中,<img file="FDA0000771904230000021.GIF" wi="57" he="77" />为二进制变量,表示第j台火电机组t时刻的启动状态,1表示机组正在启动,0表示机组不在启动状态;<img file="FDA0000771904230000022.GIF" wi="55" he="77" />也为二进制变量,表示第j台火电机组t时刻停机状态,1表示机组正在停机,0表示机组不在停机状态;<img file="FDA0000771904230000023.GIF" wi="64" he="78" />为t时刻的第j台火电机组出力;<img file="FDA0000771904230000024.GIF" wi="89" he="78" /><img file="FDA0000771904230000025.GIF" wi="62" he="77" />以及<img file="FDA0000771904230000026.GIF" wi="60" he="75" />均为自变量;N<sub>j</sub>表示参与优化火电机组的总台数;T表示内层优化模型一次仿真时间长度;α<sub>j</sub>为第j台火电机组启机煤耗;β<sub>j</sub>为第j台火电机组停机煤耗;a<sub>j</sub>为单台火电机组煤耗随功率变化斜率;b<sub>j</sub>为单台火电机组煤耗常数;γ为二氧化碳排放系数;<img file="FDA0000771904230000027.GIF" wi="275" he="82" />分别为第j台火电机组的上爬坡率和下爬坡率;<img file="FDA0000771904230000028.GIF" wi="234" he="82" />分别为第j台火电机组的最小出力值和最大出力值;<img file="FDA0000771904230000029.GIF" wi="66" he="76" />为二进制变量,表示第j台火电机组t时刻运行状态,1表示机组正在运行,0表示机组没有运行;k为预设的最小启机或停机的时间步长;<img file="FDA00007719042300000210.GIF" wi="199" he="81" />分别为供热期内,火电机组中背压式供热机组和抽气式供热机组t时刻的出力;<img file="FDA00007719042300000211.GIF" wi="71" he="78" />为t时刻热负荷大小;<img file="FDA00007719042300000212.GIF" wi="162" he="82" />为供热机组热电耦合系数;<img file="FDA00007719042300000213.GIF" wi="68" he="69" />为t时刻电网的总用电负荷;<img file="FDA00007719042300000214.GIF" wi="59" he="69" />为t时刻电网接纳的各类新能源产生的电力之和;S<sub>p</sub>、S<sub>N</sub>分别为电网正/负旋转备用;<img file="FDA00007719042300000215.GIF" wi="75" he="69" />为新能源发电在各时段的可信容量;<img file="FDA00007719042300000216.GIF" wi="76" he="71" />为t时刻电网接纳的第i类新能源产生的电力;N<sub>i</sub>为外层优化模型输入的电网中第i类新能源的装机容量;<img file="FDA00007719042300000217.GIF" wi="68" he="76" />为电网中第i类新能源长时间尺度出力时间序列的归一化值;CL为电网中新能源的类别总数;步骤3、将得到的最小二氧化碳排放量输出至外层优化模型;步骤4、外层优化模型判断是否满足预设的迭代终止条件,如是,则将二氧化碳排放量最小的各新能源的装机容量N<sub>i</sub>以及火电机组的启动状态<img file="FDA00007719042300000218.GIF" wi="89" he="76" />停机状态<img file="FDA00007719042300000219.GIF" wi="62" he="78" />和火电机组出力<img file="FDA00007719042300000220.GIF" wi="66" he="83" />作为最终的优化结果输出,优化结束;如否,则以内层优化模型输出的最小二氧化碳排放量作为适应度函数值,对各类新能源的装机容量进行更新,然后转至步骤2。
地址 210098 江苏省南京市鼓楼区西康路1号