電力系統暫態穩定最小負荷裕度的計算方法

常乃超，陶洪鑄，郭建成，辛耀中

(國家電網公司，北京100032)

電力系統暫態穩定最小負荷裕度的計算方法

常乃超，陶洪鑄，郭建成，辛耀中

(國家電網公司，北京100032)

為考慮不同的負荷—發電增長方向對電力系統暫態穩定負荷裕度的影響，提出了計算電力系統暫態穩定最小負荷裕度的最優化模型。該模型使用改進PEBS法提供的暫態穩定量化指標構成穩定約束條件，具有與暫態時域仿真相同的模型適應性。根據點到曲面最短距離的基本定理，給出了求解所提出最優化模型的數值方法。該方法的基本思想是以暫態穩定臨界點處暫態穩定邊界曲面的法向量作為求解迭代過程中的負荷—發電增長方向，反復迭代求出暫態穩定最小負荷裕度。理論上，該模型不但可以計算各種發電—負荷增長方向上的暫態穩定最小負荷裕度，也可用于計算暫態穩定最大負荷裕度。所提出的計算電力系統暫態穩定最小負荷裕度的最優化模型及其求解方法可以考慮各種系統運行參數的約束。數值算例驗證了所提出的模型及求解方法。

暫態穩定;負荷裕度;電力系統

1 引言

在線穩定分析及預警系統已成為智能電網調度技術支持系統的重要應用。目前，在線穩定分析及預警系統的主要功能是在給定的故障及過渡方案(即任意系統運行參數的變化方向，通常是考察負荷—發電增長)下，逐步惡化系統運行方式，找到系統的穩定極限，給出系統的穩定裕度［1-8］。顯然，這樣得到的穩定裕度依賴于選取的過渡方案，以負荷—發電增長方向為例，不同的負荷—發電增長方向下得到的負荷功率裕度相差可能很大。為了給調度運行人員提供更為全面深入的系統穩定裕度信息，本文探討如何求取最小負荷裕度的問題。最小負荷裕度是各種可能的負荷—發電增長方向上負荷裕度的最小值。以暫態穩定負荷裕度為例，某個故障下最小負荷裕度的意義是只要負荷增加量小于此負荷裕度，則無論所增加的負荷在各負荷母線上如何分配，系統都是暫態穩定的。

本文提出了計算暫態穩定約束下最小負荷裕度的最優化模型，并根據點到曲面最短距離的基本定理，給出了求解該模型的數值方法。求解方法的基本過程是給定初始負荷—發電增長方向后，求出該增長方向上的暫態穩定臨界點(位于注入空間中暫態穩定邊界曲面上)，以該臨界點處注入空間中暫態穩定邊界曲面的法向量作為下一次計算的負荷—發電增長方向，如此反復迭代，直至當前運行點與暫態穩定臨界點之間連線矢量的方向與暫態穩定臨界點處暫態穩定邊界曲面法向量的方向相同。

2 暫態穩定最小負荷裕度的最優化模型

為求出暫態穩定限制下的最小負荷裕度，構造如下最優化模型:

其中，x、y分別為電力系統暫態過程狀態變量和代數變量;f(·)、g(·)為描述暫態過程的微分代數方程;PL∈Rm為負荷注入向量;PG∈Rn為發電注入向量;PL，0為當前運行點負荷注入向量;PG，0為當前運行點發電注入向量;η(PL，PG)為暫態穩定量化指標;ε為暫態穩定裕度門檻值;g(·)≤0為運行約束;r∈R為負荷裕度;為負荷增長方向向量;為發電增長方向向量;=1。

暫態穩定約束下的最小負荷裕度r*表示不論負荷增量在各負荷母線上如何分配，系統至少能夠承受總和為r*的負荷增量而不發生暫態失穩。

3 求解方法

3.1 暫態穩定量化指標

電力系統暫態穩定量化指標η表征了系統運行點與暫態安全域邊界的接近程度，η接近于零意味著系統接近暫態不穩定。目前，已有很多能夠提供暫態穩定量化指標的分析方法，例如BCU、PEBS、Hybrid Method、EEAC、半 張 量 積 等 方法［8-15］。這里選擇改進PEBS方法提供的暫態穩定量化指標［12］。

改進PEBS方法的基本思想是基于時域仿真得到的系統受擾軌跡，結合PEBS方法的基本思想提取暫態穩定量化指標，綜合了時域仿真法及PEBS方法的優點，能夠考慮各種復雜的系統模型且不需要計算故障后的不穩定平衡點。

慣性中心坐標下多機系統轉子搖擺方程為

慣性中心坐標系定義為

其中，δi、ωi、Pmis、Pei、Mi分別為第i臺發電機的轉子角度、角速度、輸入機械功率、輸出電磁功率、轉子慣量。

改進PEBS方法中系統的能量函數定義為

對于暫態不穩定情形，改進PEBS方法提供的穩定量化指標為

式中，tPEBS表示系統軌跡穿越PEBS的時刻。

對于暫態穩定情形，改進PEBS方法提供的穩定量化指標為

其中，tcl為故障切除時刻;分別為故障切除時刻的勢能、動能;為持續故障下的勢能最大值。所謂持續故障是指，在時刻tins(tins＞tcl)重新插入故障勢能所能達到的最大值。tins滿足

3.2 求解過程

注入空間中系統的暫態穩定邊界曲面方程為

假定當前運行點 (PL0，PG0)與曲面上的點之間的連線長度為當前運行點(PL0，PG0)到暫態穩定邊界曲面的最短距離，則根據點到曲面距離的基本定理，點(PL0，P)G0與點之間連線的方向必與暫態穩定邊界曲面在點處的法向量方向重合。暫態穩定邊界曲面在點處的法向量為

基于暫態穩定邊界曲面的法向量，可構造如下求解過程:

(1)置k=0;

(2)給定負荷及發電增長方向dL(k)，dG(k)，在該方向上以一定步長h逐步增加負荷及發電，即PL=PL0+mh dL，PG=PG0+mh dG，m=1，2，…，Nk，當m=Nk時系統暫態失穩。dL(k)，dG(k)方向上的負荷裕度為r(k)=Nkh，暫態穩定邊界曲面上的臨界點為，其中;

(4)將負荷及發電增長方向修正為

并歸一化

(5)置k∶=k+1，返回第(2)步。

4 數值仿真

數值仿真采用如圖1所示的5機兩區域系統［16］。發電機均采用6階模型，每臺發電機均配有簡單的比例勵磁控制器和PSS模型，母線1、2、12、22上的發電機配有簡化的調速器模型，發電機1設為平衡機。母線4、14上的負荷采用恒阻抗模型，母線21上的負荷采用恒阻抗和電動機組合模型，其中電動機模型負荷占40%?；鶞使β蕿?00MVA，以下功率均為標幺值，表1給出了初始發電和負荷分布。

圖1 兩區域系統Fig．1 Two area power system

表1 發電和負荷初始分布(pu)Tab．1 Initial distribution of generation and load(pu)

負荷及發電的初始增長方向為按比例增長dL= ［0.3333，0.5，0.1667］T， dG= ［0.3043，0.3043，0.3261，0.0652］T。所考察的故障為0.1s時線路3-25在母線25側發生三相短路，0.15s時母線25側開關跳開，0.2s時母線3側開關跳開。初始發電和負荷狀態下各發電機相對于母線1上的發電機的轉子角度曲線如圖2所示。

沿負荷及發電初始增長方向按照rstep=0.05的步長增加負荷和發電，當r=2.7時，系統處于臨界穩定狀態，響應如圖3所示。

按照3.2節給出的求解過程進行計算，迭代8次后收斂，最終的負荷及發電增長方向為 dL= ［0.0564，0.0806，0.863］T，dG=［0.0546，0.0546，0.0526，0.8381］T，暫態穩定最小負荷功率裕度r=1.09，即最小負荷裕度為109MW。這意味著無論所增加的109MW負荷是以何種比例分配到各個負荷節點上，所增加的109MW負荷由哪些發電機承擔，系統都是暫態穩定的。所得到的最小負荷裕度可作為暫態穩定監控的重要指標。另外，可以看到，由于負荷21、發電機22離故障地點距離較近，因此其負荷和發電的增長對暫態穩定影響較大，表現在最終的負荷及發電增長方向中負荷21、發電機22對應的比例較大。為驗證所得到的最小負荷裕度，將118MW全部加在負荷21及發電機22上，仿真表明系統暫態穩定，如圖4所示。

圖2 初始狀態下的發電機轉子角度響應Fig．2 Generator angle response at initial state

圖3 初始增長方向下r=2.7時的角度響應Fig．3 Generator angle response under initial increasing direction(r=2.7)

5 結論

本文提出了計算電力系統暫態穩定最小負荷裕度的最優化模型，該模型中的暫態穩定約束量化指標由改進PEBS方法給出，具有與時域仿真方法相同的模型適應性，可以考慮系統復雜動態模型。

圖4 負荷21及發電機增加109MW時的角度響應Fig．4 Generator angle response when load 21 and generation 22 increasing 109MW

選定某個負荷及發電初始增長方向后，逐步增加負荷及發電得到該方向上的最小負荷裕度及暫態穩定邊界曲面上的臨界穩定點，計算該點處暫態穩定邊界曲面的法向量，并將該法向量歸一化后作為下一次迭代的負荷及發電增長方向，反復迭代得到最小負荷功率裕度。數值算例表明，所提出的模型及求解方法是有效的。本文提出的電力系統暫態穩定最小負荷裕度實際上還可看作注入空間中電力系統暫態穩定域的最短半徑。

電力系統暫態穩定最小負荷裕度應該成為在線穩定分析及預警系統的一個重要監視指標，對于調度員掌握電力系統實時運行工況穩定性的深層信息具有重要作用，有利于避免對特殊運行方式下系統穩定裕度的誤判。如何減少最小負荷裕度指標的保守性是一個值得進一步研究的課題，例如，可根據超短期負荷預測結果將負荷增長方向限制在一定范圍之內進行尋優，從而得到更貼近實際運行經驗的負荷裕度指標。

［1］Vittal V．Consequence and impact of electric utility industry restructuring on transient stability and small-signal stability analysis［J］．Proceedings of IEEE，2000，88 (2):196-207．

［2］Wu F F，Kumagai S．Steady-state security regions of power systems［J］．IEEE Transactions on Circuits and Systems，1982，29(11):703-711．

［3］Wu F F，Tsai Y-K，Yu Y．Probabilistic dynamic security assessment of power systems，Part I:Basic model［J］．IEEE Transactions on Circuits and Systems，1983，30 (3):148-159．

［4］Kaye R J，Wu F F．Dynamic security regions for power systems［J］．IEEE Transactions on Circuits and Systems，1982，29(9):612-623．

［5］Kundur P．Power system stability and control［M］．New York，USA:MacGraw-Hill，1993．

［6］倪以信，陳壽孫，張寶霖 (Ni Yixin，Chen Shousun，Zhang Baolin)．動態電力系統的理論和分析 (Theory and analysis of dynamic power systems)［M］．北京:清華大學出版社 (Beijing:Tsinghua University Press)，2002．

［7］Vittal V，Zhou E Z，Hwang C，et al．Derivation of stability limits using analytical sensitivity of the transient energy margin［J］．IEEE Transactions on Power Systems，1989，4(4):1363-1372．

［8］Chiang H D，Wu F F，Varaiya P P．A BCU method for direct analysis of power system transient stability［J］．IEEE Transactions on Power Systems，1994，9(3): 1194-1208．

［9］Chiang H D，Wu F F，Varaiya P P．Foundations of the potential energy boundary surface method for power system transient stability analysis［J］．IEEE Transactions on Circuits and Systems，1988，35(6):160-172．

［10］Sauer P W，Behera A K，Pai M A，et al．Trajectory approximations for direct energy methods that use sustained faults with detailed power system models［J］．IEEE Transactions on Power Systems，1989，4(2):499-506．

［11］Maria G A，Tang C，Kim J．Hybrid transient stability analysis［J］．IEEE Transactions on Power Systems，1990，5(2):384-391．

［12］Tang C，Graham C E，El-Kady M，et al．Transient stability index from conventional time domain simulation ［J］．IEEE Transactions on Power Systems，1994，(3): 1524-1530．

［13］Fouad A A，Kruempel K C，Vittal V，et al．Transient stability program output analysis［J］．IEEE Transactions on Power Systems，1986，1(1):2-9．

［14］Haque M H．Hybrid method of determining the transient stability margin of a power system［J］．IEE Proc．-Gener．Transm．Distrib．，1996，143(1):27-32．

［15］梅生偉，劉鋒，薛安成 (Mei Shengwei，Liu Feng，Xue Ancheng)．電力系統暫態分析中的半張量積方法(The semi-tensor product method for power systems transient stability analysis) ［M］．北京:清華大學出版社(Beijing:Tsinghua University Press)，2010．

［16］Power System Toolbox(PST)［EB/OL］．http://www．ecse．rpi．edu/pst/Downloads．html．

Method for calculating minimal load margin constrained by power system transient stability

CHANG Nai-chao，TAO Hong-zhu，GUO Jian-cheng，XIN Yao-zhong (State Grid Corporation of China，Beijing 100032，China)

In order to consider the effect of different load-generation increasing directions to the load margin constrained by power system transient stability，an optimization model for calculating the minimal load margin constrained by power system transient stability is proposed．The stability quantitative index proposed by the improved PEBS method is used to construct the transient stability constraint，and it has the same model adaptation as the time domain transient simulation method．A numerical method for solving the proposed optimization model is given out based on the basic theorem about the distance from a point to a curved surface．The basic idea of the proposed solution method is that the normal vector at the critical point of the power systems transient stability boundary surface is taken as the load-generation increasing direction in the iterative solution procedure，and the minimal load margin constrained by power systems transient stability is finally acquired．Theoretically，the proposed model and method can also be applied to calculating the maximal load margin constrained by power system transient stability，and all kinds of power system operating constraints can also be considered．A numerical example verifies the proposed model and method．

transient stability;load margin;power systems

TM712

A

1003-3076(2014)06-0047-05

2013-07-13

國家高技術研究發展計劃(863計劃)研究資助項目(2011AA05A118)

常乃超(1977-)，男，河南籍，高級工程師，博士，研究方向為電力系統分析與調度自動化;陶洪鑄(1973-)，男，安徽籍，高級工程師，碩士，研究方向為電網調度自動化。