发明名称 一种考虑交流潮流约束的不确定性机组组合模型及求解方法
摘要 本发明涉及一种考虑交流潮流约束的不确定性机组组合模型及求解方法,本发明以一天24时段中调火电机组的燃料费用总和最小为优化目标,采用机会约束方法描述风电出力的不确定性,构建基于交流潮流模型的网络安全约束条件,从而提出一种考虑交流潮流安全约束的不确定性机组组合模型。针对模型求解困难的问题,提出了一种随机约束序优化方法,成功实现对模型的快速求解。相比与传统基于直流潮流约束的不确定性机组组合模型,本发明所提模型有效降低了风电大规模接入后电网电压越限的风险,提升了日前发电计划决策的有效性,所建模型还可以精细化的计算系统网络损耗,从而为调度人员提供数据参考。同时本发明提求解算法较传统算法有了显著提升。
申请公布号 CN105790265A 申请公布日期 2016.07.20
申请号 CN201610254532.7 申请日期 2016.04.21
申请人 三峡大学 发明人 杨楠;周峥;崔家展;李宏圣;王璇;黎索亚
分类号 H02J3/00(2006.01)I;G06Q50/06(2012.01)I;G06F17/50(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 宜昌市三峡专利事务所 42103 代理人 吴思高
主权项 一种考虑交流潮流安全约束的不确定性机组组合建模方法,其特征在于,考虑风电出力的不确定性,并且以交流潮流约束对其进行安全校核,其目标函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>minF</mi><mrow><mi>G</mi><mi>t</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mo>&lsqb;</mo><msub><mi>U</mi><mrow><mi>G</mi><mi>i</mi><mi>t</mi></mrow></msub><msub><mi>Y</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>G</mi><mi>i</mi><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>U</mi><mrow><mi>G</mi><mi>i</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>U</mi><mrow><mi>G</mi><mi>i</mi><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><msub><mi>S</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&tau;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971327580000011.GIF" wi="1350" he="142" /></maths>式中:F<sub>Gt</sub>为系统的总运行成本,系统中的发电机组个数为M,P<sub>Git</sub>和U<sub>Git</sub>分别表示第i号机组在t时段的有功出力和启停状态,U<sub>Git</sub>=0表示发电机组处于停机状态,U<sub>Git</sub>=1表示发电机组处于开机状态。Y<sub>it</sub>(P<sub>Git</sub>)为发电机组的运行成本,S<sub>it</sub>(τ<sub>i</sub>)为发电机的启停成本;其中,发电机组的运行成本、启停成本的具体数学表达形式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>Y</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>G</mi><mi>i</mi><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>+</mo><msub><mi>&beta;</mi><mi>i</mi></msub><msub><mi>P</mi><mrow><mi>G</mi><mi>i</mi><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><msubsup><mi>P</mi><mrow><mi>G</mi><mi>i</mi><mi>t</mi></mrow><mn>2</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971327580000012.GIF" wi="1198" he="79" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>i</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&tau;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>S</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><msub><mi>&tau;</mi><mi>i</mi></msub><mo>/</mo><msub><mi>&omega;</mi><mi>i</mi></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="FDA0000971327580000013.GIF" wi="1187" he="71" /></maths>式中:α<sub>i</sub>、β<sub>i</sub>、γ<sub>i</sub>为机组运行成本参数,S<sub>0i</sub>、S<sub>1i</sub>、τ<sub>i</sub>为i号机组的停机时间,ω<sub>i</sub>为启停成本参数;决策变量通常满足以下常规约束条件:系统有功功率平衡约束;发电机组有功、无功出力上下限约束;最小启停次数约束;机组爬坡约束;除上述常规约束条件外,该方法提出的网络安全约束条件如下:‑P<sub>l,max</sub>≤P<sub>l,t</sub>≤P<sub>l,max</sub>               (4)‑Q<sub>l,max</sub>≤Q<sub>l,t</sub>≤Q<sub>l,max</sub>              (5)V<sub>b,min</sub>≤V<sub>bt</sub>≤V<sub>b,max</sub>              (6)式中:P<sub>l,t</sub>、Q<sub>l,t</sub>分别为线路l在t时段的有功功率和无功功率,V<sub>bt</sub>为节点b在t时段的电压,P<sub>l,max</sub>、Q<sub>l,max</sub>为分别线路l允许的最大有功、无功容量;V<sub>b,min</sub>,V<sub>b,max</sub>分别为节点b允许的最小最大电压;公式(4‑6)中的变量P<sub>l,max</sub>、Q<sub>l,max</sub>、V<sub>bt</sub>需要经过交流潮流方程求出,其具体形式如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>&Delta;</mi><mi>P</mi><mo>=</mo><msub><mi>V</mi><mi>b</mi></msub><munder><mo>&Sigma;</mo><mrow><mi>c</mi><mo>&Element;</mo><mi>b</mi></mrow></munder><msub><mi>V</mi><mi>c</mi></msub><mrow><mo>(</mo><msub><mi>G</mi><mrow><mi>b</mi><mi>c</mi></mrow></msub><msub><mi>cos&theta;</mi><mrow><mi>b</mi><mi>c</mi></mrow></msub><mo>+</mo><msub><mi>B</mi><mrow><mi>b</mi><mi>c</mi></mrow></msub><msub><mi>sin&theta;</mi><mrow><mi>b</mi><mi>c</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>&Delta;</mi><mi>Q</mi><mo>=</mo><msub><mi>V</mi><mi>b</mi></msub><munder><mo>&Sigma;</mo><mrow><mi>c</mi><mo>&Element;</mo><mi>b</mi></mrow></munder><msub><mi>V</mi><mi>c</mi></msub><mrow><mo>(</mo><msub><mi>G</mi><mrow><mi>b</mi><mi>c</mi></mrow></msub><msub><mi>sin&theta;</mi><mrow><mi>b</mi><mi>c</mi></mrow></msub><mo>-</mo><msub><mi>B</mi><mrow><mi>b</mi><mi>c</mi></mrow></msub><msub><mi>cos&theta;</mi><mrow><mi>b</mi><mi>c</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971327580000021.GIF" wi="1300" he="206" /></maths>式中:△P、△Q表示节点注入有功和无功;b=1,2,L K为节点数;c∈b表示与节点b相连的节点c;V表示节点电压;G<sub>bc</sub>、B<sub>bc</sub>分别是节点b和节点c之间的网络导纳的实部和虚部;采用机会约束方法描述机组组合问题中的风电出力不确定性,分别构建系统正负旋转备用风险指标:Q<sub>d</sub>≤λ                     (13)Q<sub>u</sub>≤λ              (14)式中:Q<sub>d</sub>、Q<sub>u</sub>分别为正、负旋转备用风险指数;λ为旋转备用风险槛值,系统调度部门可利用年总费用最小法获取可靠性标准后换算而得,通常取0~10%之间。
地址 443002 湖北省宜昌市大学路8号