发明名称 一种电力系统短期电压稳定控制方法
摘要 本发明公开了一种电力系统短期电压稳定控制方法,该方法首先通过现场测量互感器测量电力系统运行状态,经过状态估计计算后,得到系统运行的稳态运行数据,构建基于最优潮流的无功优化非线性规划模型,然后从电力系统稳定数据库中提取系统动态数据,构建描述系统动态特性的微分代数方程组及其初值方程,依据上述稳态数据和动态数据,对动态无功补偿装置进行建模,然后构造动态优化模型,最后求解基于稳定约束最优潮流的动态无功优化模型。将无功备用的优化结果编码为控制指令,通过电力系统遥控系统发送至现场发电机和无功补偿装置;本发明能够高效、准确地计算出系统运行中最优的动态无功备用调度方案,提升系统运行的短期电压稳定性。
申请公布号 CN104037775A 申请公布日期 2014.09.10
申请号 CN201410204127.5 申请日期 2014.05.14
申请人 浙江大学 发明人 耿光超;江全元;周姝灿
分类号 H02J3/16(2006.01)I 主分类号 H02J3/16(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 邱启旺
主权项 一种电力系统短期电压稳定控制方法,其特征在于,包括如下步骤:第一步:通过现场测量互感器测量电力系统运行状态,经过状态估计计算后,得到系统运行的稳态运行数据,构建基于最优潮流的无功优化非线性规划模型P<sub>1</sub>,形如:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mn>1</mn></msub><mo>:</mo><mfenced open='' close=''><mtable><mtr><mtd><munder><mrow><mi>min</mi><mi>&Phi;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mrow><mi>u</mi></munder></mtd></mtr><mtr><mtd><mi>s</mi><mo>,</mo><mi>t</mi><mfenced open='{' close=''><mtable><mtr><mtd><mi>H</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><munder><mi>G</mi><mo>&OverBar;</mo></munder><mo>&le;</mo><mi>G</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>&le;</mo><mover><mi>G</mi><mo>&OverBar;</mo></mover></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000505386400000011.GIF" wi="581" he="270" /></maths>其中,u为待优化的变量,包括发电机出力、节点电压和无功补偿方案;Φ为目标函数,通常设置为最小化燃料成本或系统网损;H为等式约束,包括潮流方程;G为不等式约束,包括发电机出力约束、节点电压约束和线路潮流约束;<img file="FDA0000505386400000012.GIF" wi="58" he="79" />和<img file="FDA0000505386400000013.GIF" wi="58" he="75" />分别为不等式约束的上下限。第二步:从电力系统稳定数据库中提取系统动态数据,构建描述系统动态特性的微分代数方程组F<sub>1</sub>及其初值方程E<sub>1</sub>,形如:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>F</mi><mn>1</mn></msub><mrow><mo>(</mo><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0000505386400000014.GIF" wi="466" he="78" /></maths>E<sub>1</sub>(x(0),u)=0其中t为仿真时间,x和<img file="FDA0000505386400000015.GIF" wi="54" he="62" />分别为状态变量及其导数,x(0)代表状态变量的初值。第三步:依据上述稳态数据和动态数据,对动态无功补偿装置进行建模。将其稳态特性加入非线性规划模型P<sub>1</sub>中,即在潮流方程中考虑动态无功补偿装置的无功输出,将其无功出力作为待优化变量,同时在不等式约束中加入其无功出力容量限制。同时,将其动态特性加入微分代数方程组F<sub>1</sub>及其初值方程E<sub>1</sub>中。从而,得到考虑动态无功补偿装置的非线性规划模型P<sub>2</sub>、微分代数方程F<sub>2</sub>及其初值方程E<sub>2</sub>,形如:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mn>2</mn></msub><mo>:</mo><mfenced open='' close=''><mtable><mtr><mtd><munder><mrow><mi>min</mi><mi>&Phi;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mrow><msup><mi>u</mi><mo>'</mo></msup></munder></mtd></mtr><mtr><mtd><mi>s</mi><mo>,</mo><mi>t</mi><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>H</mi><mo>'</mo></msup><mrow><mo>(</mo><msup><mi>u</mi><mo>'</mo></msup><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><munder><msup><mi>G</mi><mo>'</mo></msup><mo>&OverBar;</mo></munder><mo>&le;</mo><msup><mi>G</mi><mo>'</mo></msup><mrow><mo>(</mo><msup><mi>u</mi><mo>'</mo></msup><mo>)</mo></mrow><mo>&le;</mo><mover><msup><mi>G</mi><mo>'</mo></msup><mo>&OverBar;</mo></mover></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000505386400000016.GIF" wi="637" he="271" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>F</mi><mn>2</mn></msub><mrow><mo>(</mo><mover><msup><mi>x</mi><mo>'</mo></msup><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msup><mi>x</mi><mo>'</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msup><mi>u</mi><mo>'</mo></msup><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0000505386400000017.GIF" wi="510" he="73" /></maths>E<sub>2</sub>(x′(0),u′)=0其中u′为考虑动态无功补偿装置后的待优化变量,H′和G′分别为考虑动态无功补偿装置后的等式和不等式约束,<img file="FDA0000505386400000021.GIF" wi="81" he="79" />和<img file="FDA0000505386400000022.GIF" wi="72" he="77" />分别为考虑动态无功补偿装置后的不等式约束上限和下限,x′和<img file="FDA0000505386400000026.GIF" wi="56" he="59" />分别为考虑动态无功补偿装置后的状态变量及其导数。第四步:基于上述非线性规划模型P<sub>2</sub>、微分代数方程组F<sub>2</sub>和初值方程E<sub>2</sub>,构造动态优化模型P<sub>3</sub>,即将微分代数方程组F<sub>2</sub>和初值方程E<sub>2</sub>作为等式约束集成到原有非线性规划模型P<sub>2</sub>中,同时在不等式约束中加入稳定约束S。从而得到基于稳定约束最优潮流的动态无功优化模型P<sub>3</sub>。<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>P</mi><mn>3</mn></msub><mo>:</mo><mfenced open='' close=''><mtable><mtr><mtd><munder><mi>min</mi><msup><mi>u</mi><mo>'</mo></msup></munder><mi>&Phi;</mi><mrow><mo>(</mo><msup><mi>u</mi><mo>'</mo></msup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mtable></mtable><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>F</mi><mn>2</mn></msub><mrow><mo>(</mo><msup><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>'</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msup><mi>x</mi><mo>'</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msup><mi>u</mi><mo>'</mo></msup><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>E</mi><mn>2</mn></msub><mrow><mo>(</mo><msup><mi>x</mi><mo>'</mo></msup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>,</mo><msup><mi>u</mi><mo>'</mo></msup><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msup><mi>H</mi><mo>'</mo></msup><mrow><mo>(</mo><msup><mi>u</mi><mo>'</mo></msup><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msup><munder><mi>G</mi><mo>&OverBar;</mo></munder><mo>'</mo></msup><mo>&le;</mo><msup><mi>G</mi><mo>'</mo></msup><mrow><mo>(</mo><msup><mi>u</mi><mo>'</mo></msup><mo>)</mo></mrow><mo>&le;</mo><msup><mover><mi>G</mi><mo>&OverBar;</mo></mover><mo>'</mo></msup></mtd></mtr><mtr><mtd><munder><mi>S</mi><mo>&OverBar;</mo></munder><mo>&le;</mo><mi>S</mi><mrow><mo>(</mo><msup><mi>x</mi><mo>'</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&le;</mo><mover><mi>S</mi><mo>&OverBar;</mo></mover></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000505386400000023.GIF" wi="755" he="531" /></maths>其中,<img file="FDA0000505386400000024.GIF" wi="58" he="79" />和<img file="FDA0000505386400000025.GIF" wi="42" he="75" />分别为稳定约束的S的上限和下限。第五步:使用动态优化算法求解基于稳定约束最优潮流的动态无功优化模型P<sub>3</sub>。将无功备用的优化结果编码为控制指令,通过电力系统遥控系统发送至现场发电机和无功补偿装置。该控制指令得到的系统运行点可以有效规避预想故障中可能存在的短期电压稳定问题,提升系统运行的安全性。
地址 310058 浙江省杭州市西湖区余杭塘路866号