发明名称 基于能量裕度约束的电压在线协调控制方法
摘要 本发明涉及电力系统暂态稳定分析技术领域。本发明公开了一种基于能量裕度约束的电压在线协调控制方法,包括以下步骤:a、建立系统势能函数模型E<sub>p</sub>,基于李导数计算能量裕度E<sub>m</sub>;b、以能量裕度E<sub>m</sub>作为目标函数,无功裕度和电压允许变化范围作为约束条件,形成近似线性规划问题;c、获得协调控制策略,完成多种电压控制方式的协调。本发明将多种电压控制方式纳入统一的能量函数框架下,引入李导数法来计算能量裕度变化,以能量裕度不小于零作为目标函数,以无功充裕度和电压变化范围作为约束来形成近似线性规划问题,其解即为协调控制策略。本发明可以提高计算速度,有助于多种电压控制方式的在线协调控制策略的生成,提高电网的安全运行水平。
申请公布号 CN103887800A 申请公布日期 2014.06.25
申请号 CN201410103209.0 申请日期 2014.03.19
申请人 电子科技大学 发明人 刘群英;刘起方;田贵云;田云翔;孟劲松;任文伟;高斌
分类号 H02J3/12(2006.01)I 主分类号 H02J3/12(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 李顺德
主权项 1.基于能量裕度约束的电压在线协调控制方法,包括以下步骤:a、建立系统势能函数模型E<sub>p</sub>,基于李导数计算能量裕度E<sub>m</sub>;所述能量函数模型E<sub>p</sub>为:<maths num="0001"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>E</mi><mi>p</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>-</mo><mo>[</mo><mfrac><mn>1</mn><msub><mi>D</mi><mi>Li</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>Lj</mi><mn>0</mn></mrow></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>U</mi><mi>j</mi></msub><msub><mi>U</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub></mfrac><mo>)</mo></mrow><mi>&alpha;</mi></msup></mrow><mo>-</mo><msup><mi>k</mi><mn>2</mn></msup><msubsup><mi>U</mi><mi>j</mi><mn>2</mn></msubsup><msub><mi>G</mi><mi>ij</mi></msub><mo>+</mo><msubsup><mi>kU</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><mi>cos</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>-</mo><msub><mi>B</mi><mi>ij</mi></msub><mi>sin</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>[</mo><mfrac><mn>1</mn><msub><mi>M</mi><mi>i</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>P</mi><mi>mi</mi></msub><mo>-</mo><msup><mi>k</mi><mn>2</mn></msup><msubsup><mi>U</mi><mi>j</mi><mn>2</mn></msubsup><msub><mi>G</mi><mi>ij</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>-</mo><msub><mi>&delta;</mi><mrow><mi>ij</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>kU</mi><mi>j</mi><mn>2</mn></msubsup><mi>G</mi><mrow><mo>(</mo><mi>sin</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>-</mo><mi>sin</mi><msub><mi>&delta;</mi><mrow><mi>ij</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>[</mo><msubsup><mi>kU</mi><mi>j</mi><mn>2</mn></msubsup><msub><mi>B</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>cos</mi><msub><mi>&delta;</mi><mi>ij</mi></msub><mo>-</mo><mi>cos</mi><msub><mi>&delta;</mi><mrow><mi>ij</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mfrac><mn>1</mn><msub><mi>T</mi><mi>j</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>Q</mi><mrow><mi>Lj</mi><mn>0</mn></mrow></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>U</mi><mi>j</mi></msub><msub><mi>U</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub></mfrac><mo>)</mo></mrow><mi>&beta;</mi></msup><mo>-</mo><msub><mi>Q</mi><mi>ei</mi></msub><mo>)</mo></mrow><mi>ln</mi><msub><mi>U</mi><mi>j</mi></msub><msubsup><mo>|</mo><msub><mi>U</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><msub><mi>U</mi><mi>j</mi></msub></msubsup><mo>]</mo></mtd></mtr><mtr><mtd><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></munderover><mfrac><msub><mi>B</mi><mi>cj</mi></msub><mn>2</mn></mfrac><mo>[</mo><msup><mrow><mo>(</mo><msub><mi>U</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>U</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mtd></mtr></mtable></mfenced></math>]]></maths>其中:D<sub>Li</sub>是节点i的阻尼系数;U<sub>i</sub>、U<sub>j</sub>是节点i、j的电压;U<sub>i0</sub>、U<sub>j0</sub>是节点i、j的初始电压;α、β分别是有功和无功的电压比例系数;k是有载调压分接头的调节率;Gi<sub>j</sub>、Bi<sub>j</sub>是节点i、j之间的电导和电纳;ω<sub>i</sub>是发电机i的角速度;ω<sub>i0</sub>是发电机的初始角速度;M<sub>i</sub>是发电机i的惯性常数;T<sub>j</sub>是负荷节点j的负荷恢复系数;Bc<sub>j</sub>是节点j并联电容器组的补偿电纳;P<sub>mi</sub>是节点i的机械功率;P<sub>Lj</sub>、Q<sub>Lj</sub>是节点j的有功和无功负荷;Q<sub>ei</sub>是发电机节点i的无功注入;δ<sub>ij0</sub>为节点i、j之间的初始相位差;δ<sub>ij</sub>为节点i、j之间的相位差;所述能量裕度E<sub>m</sub>为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>E</mi><mi>m</mi></msub><mo>=</mo><msubsup><mi>E</mi><mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>&delta;</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>,</mo><msub><mi>U</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><mo>(</mo><msub><mi>&delta;</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>,</mo><msub><mi>U</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow></msubsup><mo>-</mo><msubsup><mi>E</mi><mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>&delta;</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>U</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><mo>(</mo><mi>&delta;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>U</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></msubsup><mo>=</mo><msubsup><mi>E</mi><mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>02</mn></msub><mo>)</mo></mrow></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>u</mi></msub><mo>)</mo></mrow></msubsup><mo>-</mo><msubsup><mi>E</mi><mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>01</mn></msub><mo>)</mo></mrow></mrow><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></msubsup></mrow></math>]]></maths>其中,x<sub>02</sub>:(δ<sub>s2</sub>,U<sub>s2</sub>)和x<sub>u</sub>:(δ<sub>u2</sub>,U<sub>u2</sub>)是故障后的稳定平衡点和不稳定平衡点;x<sub>01</sub>:(δ<sub>s1</sub>,U<sub>s1</sub>)是故障前的稳定平衡点;x:(δ(t),U(t))是在时间t时刻的当前运行点;<img file="FDA0000479189800000013.GIF" wi="537" he="99" />是故障后的能量裕度;<img file="FDA0000479189800000014.GIF" wi="482" he="101" />是在时间t的当前势能值;b、以能量裕度E<sub>m</sub>作为目标函数,无功裕度和电压允许变化范围作为约束条件,形成近似线性规划问题;c、获得协调控制策略,完成多种电压控制方式的协调。
地址 611731 四川省成都市高新西区西源大道2006号