電力系統負荷節點在線電壓穩定指標研究

高峰, 王彥文，郭林，徐睿

(1.中國礦業大學(北京)，北京市 100083；2.國網山東省電力公司德州供電公司，山東省德州市 253000)

電力系統負荷節點在線電壓穩定指標研究

高峰1, 王彥文1，郭林1，徐睿2

(1.中國礦業大學(北京)，北京市 100083；2.國網山東省電力公司德州供電公司，山東省德州市 253000)

利用同步相量測量單元(phasor measurement unit，PMU)測得任意(非PV)節點及其相鄰母線的電壓、線路潮流等實時數據，建立一種新的網絡等值模型。同時，基于該模型提出了綜合考慮負荷節點有功功率和無功功率的在線電壓穩定指標；提出在線電壓穩定有功功率指標和在線電壓穩定無功功率指標分別用作表征負荷節點的有功功率裕度和無功功率裕度，作為對節點(系統)在線電壓穩定指標的補充。如果負荷節點在線電壓穩定指標值接近臨界值1，則表明系統臨近電壓崩潰。EPRI-36節點仿真分析表明該方法可有效判斷系統電壓穩定性，適用于在線電壓穩定監測和預估。與其他指標的對比結果體現了本文模型和指標的合理性與優越性。最后通過對山東電網500 kV主網負荷節點電壓穩定裕度評估簡要說明本文所提方法的可行性。

電壓穩定；在線電壓穩定指標；同步相量測量單元(PMU)；網絡等值模型

0 引 言

近年來世界范圍內發生了多起大停電事故，給社會和經濟帶來了巨大的損失。停電事故的原因是多方面的，但最終主要表現為電壓失穩和電壓崩潰[1-4]。為了預防電壓失穩和電壓崩潰，系統電壓穩定裕度，即當前電力系統運行狀態離崩潰點的距離，成為電力調度工作人員最關心的問題。因此，制定一個準確衡量電壓穩定程度的指標，對于有效體現電壓穩定裕度，預警電壓失穩和電壓崩潰至關重要。現有的靜態電壓穩定性指標主要有靈敏度指標[5-6]、最小特征值(奇異值)指標[7]、試驗函數指標[8]、負荷裕度指標[9-10]、L指標[11]等。但這些指標普遍存在計算量大、準確性低、線性度差的問題，都不能在線預警電壓失穩[12-13]。

為了實用、簡捷地快速判斷系統的電壓穩定裕度，基于上述靜態電壓穩定性指標，許多新的電壓穩定指標已被提出。文獻[14]提出的算法關鍵在于求取電網絡戴維南等效參數，并假設等效參數在2個或多個系統狀態間是不變的，計算誤差較大；文獻[15]基于雅可比矩陣，從負載母線臨界電壓崩潰點的電壓相角和無功損耗特征，推導出的指標算法不適用于電壓穩定的在線監測；文獻[16]提出了一種利用母線電氣參數求取電壓穩定裕度的新算法，其算法精度嚴重依賴于母線電氣參數的連續求取；文獻[17]提出一種新的L指標，但該推導方法采用的線路模型忽略了支路對地電容的影響，不適用于高壓大電網系統；文獻[18]采用的П型線路模型雖然考慮了支路對地電容，但未考慮系統對被測線路的影響，結果誤差較大；文獻[19]提出的電壓穩定指標建立在等值電壓源和等效負荷功率因數不變的基礎上，指標無法跟蹤體現系統影響下的節點電壓變化。

針對以上問題，本文提出一種新的電力系統在線電壓穩定指標。該指標基于同步相量測量單元(phasor measurement unit，PMU)測得的實時數據，建立被測負荷節點的相關支路網絡等值模型，考慮相關支路對被測線路的影響。同時，相對于求取網絡戴維南等效參數，需要求取的等值阻抗只是該模型阻抗的一部分，減小了網絡等值模型誤差，提高了模型的準確度。更重要的是，新指標綜合考慮負荷節點的有功功率和無功功率極限裕度，能適用于線路傳輸功率因數和負荷節點功率因數變化的情況；能在系統發生電壓崩潰時接近于1，且具有較高的精度和線性度；能直接估算出節點負荷的電壓穩定裕度，無須計算系統潮流；計算速度快，便于在線實用。

1 靜態電壓穩定分析模型

圖1 被測節點及其相連線路Fig.1 Measured node with network around

圖2 被測節點及其相連線路等值模型Fig.2 Equivalent model of measured node with network around

2 在線電壓穩定指標

(1)

(2)

分解Pn、Qn得：

(3)

(4)

式中δ=δeq-δn。將Yeq= 1/Zeq代入式(3)、式(4)可分別得：

(5)

ZeqQn=0

(6)

2.1 在線電壓穩定指標

(7)

由式(7)可知，為保證Un有2個正實根，要求其判別式大于等于0，即

(8)

式(8)等價于

(9)

將節點在線電壓穩定指標(on-linevoltagestabilityindex)定義為

(10)

整個系統的在線電壓穩定指標定義為

(11)

式中：LnVSI為時間斷面t時，節點n的電壓穩定指標；S代表系統的節點集。用LnVSI的最大值LsVSI表征系統的電壓穩定指標，系統首先在這里發生電壓崩潰。LsVSI的值越接近1，則節點(系統)越接近電壓崩潰點。

式(11)是考慮了負荷節點有功功率Pn和無功功率Qn而提出的電壓穩定指標，體現了負荷節點在(P，Q，U)空間曲面[20]上運行時，對曲面臨界點的接

近程度。但是，式(11)僅能反映負荷節點電壓穩定的“綜合指標”，無法分別體現負荷節點的有功功率和無功功率裕度。

而PV曲線是負荷節點Pn、Qn、Un三者關系曲面(屬于(P，Q，U)空間)在(P，U)平面上的投影，反映了一定PQ關系下(比如一定功率因數)，節點的有功功率裕度；VQ曲線是曲面在(U，Q)平面上的投影，反映了P值一定時，節點的無功功率裕度[20-22]。由此，本文在前述等值模型的基礎上，利用文獻[18]和文獻[23]所提的電壓穩定判據提出電壓穩定有功功率指標和無功功率指標，作為對式(10)和式(11)的補充，也作為與本文所提指標的對比。

2.2 在線電壓穩定有功功率指標

當式(5)的根判別式大于等于0時，即

(12)

式(5)有2個實數解：

(13)

假設0 <θeq-δ< 90°(一般情況下該命題亦是成立的)，且0 <θeq< 90°，為保證式(5)有2個正實根，式(12)成立即可。式(12)等價于

(14)

將式(14)定義為節點n的在線電壓穩定有功功率指標(on-linevoltagestabilityindexofactivepower)[18]，即

(15)

式中LnVSIP表征節點n在時間斷面t時，在PV曲線上對曲線拐點的接近程度，反映節點的有功功率裕度。LnVSIP的值越接近1，則節點越接近電壓崩潰點(有功功率負荷能力極限)。

2.3 在線電壓穩定無功功率指標

當式(6)的根判別式大于等于0時，即

ΔQn=[Ueqsin(θeq-δ)]2-

4QnZeq(sinθeq-ZeqBeq/2)≥0

(16)

式(6)有2個實數解：

(17)

而0 <θeq-δ< 90°，0 <θeq< 90°，sinθeq≥ZeqBeq/2。為保證式(6)有2個正實根，式(16)成立即可。式(16)等價于

(18)

將式(18)定義為節點n的在線電壓穩定無功功

率指標(on-linevoltagestabilityindexofreactivepower)[23]，即

(19)

式中LnVSIQ表征節點n在時間斷面t時，在VQ曲線上對曲線拐點的接近程度，反映節點的無功功率裕度。LnVSIQ的值越接近1，則節點越接近電壓崩潰點(無功功率負荷能力極限)。

具有LnVSIP最大值的負荷節點，表示該節點是系統內有功功率裕度最小的負荷節點；具有LnVSIQ最大值的負荷節點，表示該節點是系統內無功功率裕度最小的負荷節點；具有LnVSI最大值的負荷節點，則表示該節點電壓穩定裕度最小，最有可能引起整個系統的電壓崩潰。

3 仿真結果與分析

本文以WEPRI-36系統為例，采用綜合電力系統分析綜合程序(power system analysis software package，PSASP)進行數值仿真。本文首先以2種負荷波動形式進行仿真分析：(1)所有節點負荷以初始功率的10%增加，主要通過該負荷增加方式分析系統最弱支路；(2)單獨增加某一節點負荷(增加量保持恒功率因數)，主要通過該負荷增加方式分析負荷節點間的相互影響。最后與其他指標進行對比，并利用本文所提指標對山東電網500 kV主網負荷節點進行電壓穩定裕度評估。

3.1 所有節點負荷以初始功率的10%增加

系統在原負荷的基礎上，各節點負荷以初始功率的10%增加(保持負荷原功率因數不變)，直到系統臨近電壓崩潰。PSASP中，采用這種負荷增加方式，系統發生電壓崩潰時，不僅負荷節點的功率增加倍數相同，而且某個負荷節點有功功率增加倍數和無功功率增加倍數也相同。

本文對具有代表性的負荷節點9、16、18、20、21和29的電壓穩定性進行判斷分析。

(1)LnVSI曲線分析

綜合考慮節點有功功率和無功功率，判斷負荷節點電壓穩定裕度。圖3為系統電壓崩潰前節點LnVSI曲線圖。

從圖3中可以看出，隨著功率的增加，L20VSI到達最大值0.997 6 pu，較其他LnVSI首先接近于臨界值1，據此判定是節點20發生電壓失穩導致系統電壓崩潰。從圖3中可以看出，最弱負荷節點(綜合考慮節點有功功率和無功功率時，電壓穩定裕度最低的負荷節點)從節點21轉移到了節點20，直至系統電壓崩潰。

圖3 各負荷節點以其初始值的10%增加時的 電壓穩定指標Fig.3 Each load node and voltage stability index with load increasing by 10% of initial value

(2)LnVSIP和LnVSIQ分析

系統各負荷節點以其初始功率的10%增加，負荷節點電壓Un隨其無功功率Qn與初始無功功率Qno的比值變化曲線(Un-Qn/Qno曲線)與Un-Pn/Pno曲線相同，只是橫坐標代表Qn/Qno，這里僅以Un-Pn/Pno曲線為例說明。但是節點LnVSIQ曲線和LnVSIP曲線是不同的，后文將進行分析說明。

圖4為負荷節點電壓Un隨其有功功率Pn與初始有功功率Pno的比值變化曲線(Un-Pn/Pno曲線)。從圖4可以看出，BUS16一直保持著較高的電壓水平，這是因為BUS16連接并聯電容器，具有充足的無功功率，同時BUS16和6條支路相連接，是連接支路最多的節點，能夠很好地得到周邊節點的功率/電壓支撐。BUS16的電壓落差值(初始電壓幅值與電壓崩潰時刻的電壓幅值之差)是最小的，但電壓穩定指標L16VSI并不是最大的(見圖3)；BUS9的電壓落差值最大，但其電壓穩定指標L9VSI也不總是最小的(見圖3中L9VSI前半段)。這說明不能簡單地以負荷節點電壓降落大小來判斷系統電壓穩定水平。

圖5為負荷節點在系統電壓崩潰前的電壓穩定有功功率指標變化曲線。隨著負荷節點有功功率的增加(Pn/Pno值增大)，負荷節點20的LnVSIP曲線L20VSIP達到最大值0.976 1 pu，較其他LnVSIP首先接近于臨界值1。圖5中節點20和節點21的電壓穩定有功功率指標曲線L20VSIP和L21VSIP的交點，與圖4中其Un-Pn/Pno曲線BUS20和BUS21的交點很明顯并不在同一時刻(相同的Pn/Pno值)。這說明不能以有功功率、電壓幅值是否相同來判斷2個負荷節點是否具有相同的電壓穩定有功功率裕度。

負荷節點在系統電壓崩潰前的電壓穩定無功功率指標變化曲線如圖6所示。從圖中可以看出，隨著負荷節點無功功率的增加(Qn/Qno值增大)，負荷節點29的LnVSIQ曲線L29VSIQ達到最大值0.968 9 pu，較其他LnVSIQ首先接近于臨界值1。需要注意的是：在圖5中節點29的電壓穩定有功功率指標L29VSIP比較低，但圖6中節點29的電壓穩定無功功率指標L29VSIQ不但初始值較高，而且最終發展為系統最大值。可見，同一負荷節點對有功功率和無功功率的敏感程度是不一樣的，需要對節點的有功功率裕度和無功功率裕度區別對待。

圖4 各負荷節點以其初始值的10%增加時的 電壓值(有功功率部分)Fig.4 Each load node and voltage (of active power) with load increasing by 10% of initial value

圖5 各負荷節點以其初始值的10%增加時的 電壓穩定有功功率指標Fig.5 Each load node and voltage stability index of active power with load increasing by 10% of initial value

(3)綜合分析

1)Un-Pn/Pno(和Un-Qn/Qno)曲線位置和變化趨勢十分接近的負荷節點，其電壓穩定指標曲線，包括LnVSI、LnVSIP和LnVSIQ，可能相差較大(如圖4中的BUS18和BUS21圖形十分相似且位置接近，但是圖3、圖5、圖6中的兩者的指標曲線都相差較大)。這是因為負荷節點各自周邊網絡線路及其負載情況不同，導致節點電壓崩潰點、電壓穩定指標曲線走勢等會有較大差別。采用網絡等效模型，計及被測負荷節點周邊線路的影響體現于此。

圖6 各負荷節點以其初始值的10%增加時的 電壓穩定無功功率指標Fig.6 Each load node and voltage stability index of reactive power with load increasing by 10% of initial value

2)在圖5中，隨著有功功率的增加，L20VSIP首先達到電壓崩潰點；而在圖6中，隨著無功功率的增加，L29VSIQ首先達到電壓崩潰點。

3)以電壓穩定有功功率指標判斷系統最弱支路情況是從節點21向節點20轉移，如圖5所示；而以電壓穩定無功功率指標判斷系統最弱支路情況是節點20向節點29轉移，如圖6所示。

系統通過不同的功率指標(有功功率指標或無功功率指標)得出不同的最弱節點和不同的電壓穩定裕度值。這是因為負荷節點的Un-Pn/Pno和LnVSIP曲線(或Un-Qn/Qno和LnVSIQ曲線)只是從有功功率(或無功功率)側面體現節點的功率極限和裕度，沒有綜合考慮有功功率和無功功率同時影響。因此，如文獻[18]、[24]和[25]等僅僅依靠有功功率或者無功功率作出的電壓穩定指標來判斷系統的電壓穩定裕度是不準確的，甚至是矛盾的。這為本文判據式(10)的提出提供了實驗依據。

3.2 恒功率因數增加節點負荷

以0.95的恒功率因數逐漸增加節點21的負荷，其他節點的負荷保持不變，直到系統臨近電壓崩潰。

(1)LnVSI曲線分析

圖7為系統電壓崩潰前節點在線電壓穩定指標曲線圖。圖7中可以看出，隨著節點21負荷功率的增加，LnVSI到達最大值(0.995 0 pu)，較其他LnVSI首先接近于臨界值1，據此判定是節點21發生電壓失穩導致系統電壓崩潰。

在臨近電壓崩潰時，L16VSI快速上升且具有較大值，最終僅次于L21VSI，而在圖3中，L16VSI一直保持較低值。節點16與無功功率補償器相連并且連接6條支路，可以說是負荷節點中可以獲得來自周邊節點的功率/電壓支撐最多的節點，但是因受節點21牽連，電壓穩定裕度仍舊急速減少。可見，一旦某一負荷節點發生電壓失穩，則會大大降低附近節點(特別是為系統提供無功功率輸出的節點)的電壓穩定裕度，并波及附近區域乃至整個網絡。

圖7 以恒功率因數0.95增加節點21負荷時的電壓穩定指標Fig.7 Voltage stability index with the 21st bus’s load increasing by the power factor 0.95

(2)LnVSIP和LnVSIQ分析

圖8為負荷節點電壓Un隨節點21負荷有功功率P21與其初始有功功率P21o的比值變化曲線(Un-P21/P21o曲線)。圖9為負荷節點電壓Un隨節點21負荷無功功率Q21與其初始無功功率Q21o的比值變化曲線(Un-Q21/Q21o曲線)。從圖8、圖9看出，系統節點的Un-P21/P21o曲線和Un-Q21/Q21o曲線十分相似，僅僅是橫坐標最大值不一樣，即節點有功功率和無功功率的增加倍數不一樣。

對比圖4和圖8、圖9可以發現：BUS9的電壓跌落幅值明顯減小；BUS29和BUS18的交點后移；各節點在臨近電壓崩潰點處的電壓幅值也發生不同變化。

圖8 以恒功率因數0.95增加節點21負荷時的 電壓值(有功功率部分)Fig.8 Bus voltage against the 21st bus’s load increasing by the power factor 0.95 (active part)

圖9 以恒功率因數0.95增加節點21 負荷時的電壓值(無功功率部分)Fig.9 Bus voltage against the 21st bus’s load increasing by the power factor 0.95 (reactive part)

這說明不同的負荷增加方式對整個系統電壓的影響是不同的，負荷節點的電壓穩定裕度模型需要從整個網絡的層面去考慮。

圖10和圖11分別為負荷節點在系統電壓崩潰前的電壓穩定有功功率指標和無功功率指標。隨著節點21負荷有功功率和無功功率的增加(P21/P21o和Q21/Q21o值增大)，負荷節點21的電壓穩定指標值L21VSIP和L21VSIQ較其他節點指標值首先接近于臨界值1。與圖6中的L29VSIQ不同，圖11中L29VSIQ沒有隨負荷的增加快速上升而超過其他曲線，這是因為負荷增加的節點21離節點29較遠，還不足以引起L29VSIQ超過L21VSIQ。但是，可以明顯看出L29VSIQ快速上升，臨近電壓崩潰時甚至超過了L20VSIQ，僅次于L21VSIQ。可見節點29對系統無功功率變化十分敏感，需要引起重視。

圖10 以恒功率因數0.95增加節點 21負荷時的電壓穩定有功功率指標Fig.10 Voltage stability index of active power with the 21st bus’s load increasing by the power factor 0.95

(3)綜合分析

1)節點18的電壓穩定有功功率指標和無功功率

圖11 以恒功率因數0.95增加節點 21負荷時的電壓穩定無功功率指標Fig.11 Voltage stability index of reactive power with the 21st bus’s load increasing by the power factor 0.95

指標曲線(L18VSIP和L18VSIQ)都較低，甚至其電壓穩定指標曲線(L18VSI)表現為系統最低，但是其電壓幅值下降值卻不是最小的，再次說明不能簡單地以負荷節點電壓降落大小來判斷系統電壓穩定水平。

2)圖11中，隨著無功功率的增加，系統無功功率最弱節點從節點20向節點21轉移，較好地反映了系統的無功功率最弱節點變化情況。而在這方面，LnVSI指標是做不到的，因此把節點的LnVSIQ指標(和LnVSIP指標)作為對LnVSI指標的補充是有必要的。

通過對比2個仿真算例還可以發現：節點18在兩次仿真中都表現出較高的電壓穩定性(L18VSI、L18VSIP和L18VSIQ曲線都較低)，這是因為節點18通過變壓器與電源相連(電氣距離為j 0.037 5 pu)，還與節點16相連(電氣距離為0.003 3 + j 0.033 3 pu)，節點16處安裝有無功功率補償裝置，電源和無功功率補償裝置對節點電壓具有支撐作用，使得節點18的電壓始終保持在較高水平。

系統負荷節點的Un-Pn/Pno曲線和Un-Qn/Qno曲線非常相似，這是因為PSASP采用節點有功功率和無功功率同步增長的方式進行仿真的；而節點的電壓穩定有功功率指標曲線和無功功率指標曲線相差較大，這是因為它們由不同的變量因子決定的(具體表達式見式(15)和式(19))。

3.3 與其他指標對比分析

式(20)中的4個電壓穩定指標[18,23-25]在電壓崩潰點處臨界值也為1。

(20)

為使指標之間具有可對比性，上述4個指標也采用本文提供的網絡等值模型。在上述2個仿真事例中同時計算5個指標，系統臨近電壓崩潰點時的計算結果如表1所示。

表1 臨近系統崩潰點時的各指標值

Table 1 Each inex on system critical point

由表1可知，本文研究提出的指標能較好地反映系統臨近電壓崩潰時的穩定狀態，其他指標則相對較差。這是因為這些指標不能綜合考慮被測節點的有功功率和無功功率變化，與實際系統存在較大差別。

3.4 實例分析

本小節以山東電網500 kV主網為例(如圖12所示)，按照前述2種負荷波動形式對山東省500 kV主網進行計算機仿真，簡要說明本文所提指標在實際系統中的適用性。系統臨近電壓崩潰點時的計算結果如表2所示。表2中，事例1為所有節點負荷以初始功率的10%增加的情況；事例2為單獨增加淄博節點負荷(增加量保持恒功率因數)的情況。

從表2中可以看出，當全網節點負荷同步增加時，德州節點率先達到電壓不穩定點。同時，濰坊、魯中節點也具有較小的電壓穩定裕度。單獨增加淄博節點負荷，使得該節點電壓失穩，濰坊節點LnVSI指標受其影響急劇上升，雖未超越淄博節點的LnVSI指標，但從中可見節點電壓失穩對其周邊節點，尤其是對周邊重負荷節點影響甚大(濰坊節點在事例1中就表現出較低的電壓穩定裕度)。

在仿真中發現，按山東電網500 kV主網各節點規劃負荷進行電壓穩定裕度評估，其3項指標均較低，體現出山東500 kV主網具有較高的電壓穩定裕度。

圖12 山東電網500 kV主網架接線圖(2015年)Fig.12 500 kV main network of Shandong grid (2015)表2 山東電網500 kV主網臨近系統崩潰點時各節點指標值Table 2 Nodal indices on system critical point of 500 kV Shandong grid

綜上所述，在系統負荷逐漸變化的過程中，指標LnVSI能較好地反映系統靜態電壓穩定裕度的變化情況。在系統趨于電壓崩潰的過程中，總存在LnVSI值最大的負荷節點，其值跟隨電壓變化并趨近于1。

4 結 論

本文提出的節點電壓穩定指標采用了網絡等值模型，考慮了被測負荷節點相關支路網絡的影響，更接近輸電網絡的實際情況；指標參數可以由PMU實時測得，可實現在線電壓穩定監測和預估。本文提出使用在線電壓穩定有功功率指標LnVSIP和在線電壓穩定無功功率指標LnVSIQ分別表征節點的有功功率裕度和無功功率裕度，并由綜合考慮了有功功率和無功功率的電壓穩定指標LnVSI作為負荷節點的在線電壓穩定裕度指標，計算簡單、快速，具有一定的工程實用價值。

數值仿真分析結果表明：在線電壓穩定有功功率指標和在線電壓穩定無功功率指標反映了節點電壓分別關于其有功功率和無功功率的穩定狀況；在線電壓穩定指標綜合反映了電力系統的電壓穩定狀況。

[1]薛禹勝.綜合防御由偶然故障演化為電力災難—北美“8.14”大停電的警示[J].電力系統自動化，2003，27(18)：1-5. XUE Yusheng.The way from a simple contingency to system-wide disaster: lessons from the eastern interconnection blackout in 2003[J]. Automation of Electric Power Systems，2003，27(18)：1-5.

[2]唐桃波, 夏云非, 魯文, 等.美國近年的停電事故及對我國電力系統安全穩定運行的啟示[J].電力建設，2003，24(11)：2-4. TANG Taobo, XIA Yunfei, LU Wen, et al. Power failure of this year in USA and its inspiration to safe stable operation of power system in China[J]. Electric Power Construction, 2003, 24(11): 2-4.

[3]阮前途，馬若新，肖飛，等. 菲律賓呂宋電網系統崩潰及黑啟動措施分析[J]. 電力系統自動化，2011, 35(4): 82-91. RUAN Qiantu, MA Ruoxin, XIAO Fei, et al.Blackout and restoration experience in Luzon Grid[J]. Automation of Electric Power Systems，2011，35(4)：82-91.

[4]湯涌，卜廣全，易俊.印度“7.30”、“7.31”大停電事故分析及啟示[J] .中國電機工程學報，2012，32(25)：167-174. TANG Yong，BU Guangquan，YI Jun.Analysis and lessons of the blackout in Indian power grid on July 30 and 31, 2012[J].Proceedings of the CSEE，2012，32(25)：167-174.

[5]李佳，劉天琪，陳亮，等.基于理想點法的多準則綜合靈敏度電壓穩定評估指標[J].電力自動化設備，2014，34(3)：108-112. LI Jia，LIU Tianqi，CHEN Liang，et al.Multi-criterion integrated-sensitivity voltage stability evaluation index based on ideal point method[J]. Electric Power Automation Equipment，2014，34(3)：108-112.

[6]李帥虎，曹一家，劉光曄，等.基于電壓穩定在線監測指標的預防控制方法[J].中國電機工程學報, 2015, 35(18): 4598-4606. LI Shuaihu, CAO Yijia, LIU Guangye, et al.Preventive control method of electric power system based on on-line voltage stability monitoring index[J]. Proceedings of the CSEE, 2015, 35(18): 4598-4606.

[7]SMON I, PANTOS M, GUBINA F. An improved voltage-collapse protection algorithm based on local phasors[J].Electric Power Systems Research，2008，78(3): 434-440.

[8]王慶紅，CHEN K. 基于試驗函數的霍普夫分岔搜索新算法與電壓穩定指標分析[J].電工技術學報，2008，23(1)：120-128. WANG Qinghong，CHEN K. Novel algorithms based on test function for tracing hopf bifurcation and analysis of power system voltage stability index[J].Transactions of China Electrotechnical Society，2008，23(1)：120-128.

[9]周任軍，吳潘，童小嬌.極限誘導分岔最小負荷裕度計算方法[J].電力自動化設備，2010，30(7)：19-23. ZHOU Renjun，WU Pan，TONG Xiaojiao.Calculation of minimum load power margin for limit induced bifurcation [J]. Electric Power Automation Equipment，2010，30(7)：19-23.

[10]張曉英, 蔣拯, 丁寧, 等. 負荷裕度及其與典型參數的靈敏度計算[J].電機與控制學報, 2015, 19(8): 28-35. ZHANG Xiaoying, JIANG Zheng, DING Ning, et al. Load margin and its sensitivity calculation with typical parameters[J]. Electric Machines and Control, 2015, 19(8): 28-35.

[11]姜濤，李國慶，賈宏杰，等.電壓穩定在線監控的簡化L指標及其靈敏度分析方法[J].電力系統自動化，2012，36(21)：13-18. JIANG Tao，LI Guoqiang，JIA Hongjie，et al.Simplified L-index and its sensitivity analysis method for on-line monitoring of voltage stability control[J].Automation of Electric Power Systems，2012，36(21)：13-18.

[12]蘇永春，程時杰，文勁宇，等. 電力系統電壓穩定性及其研究現狀(一)[J].電力自動化設備，2006，26(6)：97 - 101. SU Yongchun，CHENG Shijie，WEN Jinyu，et al. Power system voltage stability and its present investigation (I)[J]. Electric Power Automation Equipment，2006，26(6)：97-101.

[13]吳政球, 李日波, 鐘浩, 等. 電力系統靜態電壓穩定極限及裕度計算綜述[J]. 電力系統及其自動化學報, 2010, 22(1): 126-132. WU Zhengqiu, LI Ribo, ZHONG Hao, et al. Summary of power system’s static voltage stability limitation and load margin calculation[J]. Proceedings of the CSU-EPSA, 2010, 22(1): 126-132.

[14]LIU J H, CHU C C. Wide-area measurement-based voltage stability indicators by modified coupled single-port models[J]. IEEE Transactions on Power Systems, 2014, 29(2): 756-764.

[15]GHIOCEL S G, CHOW J H. A power flow method using a new bus type for computing steady-state voltage stability margins[J].IEEE Transactions on Power Systems, 2014, 29(2): 958-965.

[16]HAZARIKA D. New method for monitoring voltage stability condition of a bus of an interconnected power system using measurements of the bus variables[J]. IET Generation, Transmission & Distribution, 2012, 6(10): 977-985.

[17]ALTHOWIBI F A, MUSTAFA M W. Voltage stability calculations in power transmission lines：Indications and allocations[C]// 2010 IEEE International Conference on Power and Energy (PECon). Kuala Lumpur: IEEE, 2010: 390-395.

[18]劉道偉, 謝小榮, 穆鋼, 等. 基于同步相量測量的電力系統在線電壓穩定指標[J]. 中國電機工程學報, 2005, 25(1): 13-17. LIU Daowei, XIE Xiaorong, MU Gang, et al. An on-line voltage stability index of power system based on synchronized phasor measurement [J]. Proceedings of the CSEE, 2005, 25(1): 13-17.

[19]汪洋，盧繼平，李文沅，等.基于局部網絡電壓相量的等值模型及其電壓穩定性指標[J].中國電機工程學報，2008，28(34)：52-58. WANG Yang, LU Jiping, LI Wenyuan, et al. An equivalent model and voltage stability index based on local network voltage phasors [J].Proceedings of the CSEE，2008，28(34)：52-58.

[20]湯涌.電力系統電壓穩定性分析[M].北京：科學出版社，2011：107-110.

[21]周雙喜，朱凌志，郭錫玖，等.電力系統電壓穩定性及其控制[M].北京：中國電力出版社，2003：162-168.

[22]JOHNSON M, VAZIRI M, VADHVA S. Understanding slow voltage instability [C]//2012 IEEE 13th International Conference on Information Reuse and Integration (IRI). Las Vegas, NV: IEEE, 2012: 502-508.

[23]MOGHAVEMMI M, OMAR F M. Technique for contingency monitoring and voltage collapse prediction [J]. IEE Proceedings-Generation, Transmission and Distribution. 1998, 145(6): 634-640.

[24]MUSIRIN I， RAHMAN T K A.Novel fast voltage stability index (FVSI) for voltage stability analysis in power transmission system[C]// Research and Development, 2002. SCOReD2002. Student Conference on. Shah Alam，Malaysia: IEEE, 2002: 265-268.

[25]MUSIRIN I, RAHMAN T K A.Estimating maximum load ability for weak bus identification using FVSI [J].IEEE Power Engineering Review, 2002, 22(11): 50-52.

(實習編輯 景賀峰)

Load-Nodal On-Line Voltage Stability Index for Power System

GAO Feng1, WANG Yanwen1, GUO Lin1, XU Rui2

(1. China University of Mining & Technology (Beijing), Beijing 100083, China;2. Dezhou Power Supply Company,State Grid Shandong Electric Power Company, Dezhou 253000, Shandong Province, China)

This paper applies phasor measurement unit (PMU) to measure the real-time data including the voltage and power flow of any (non PV) node and adjacent busbar, and constructs a network equivalent model. Meanwhile, we present the on-line voltage stability index with comprehensive consideration of the active power and reactive power of load node based on this model. And we suggest that the on-line voltage stability indexes of active power and reactive power can be used as active power or reactive power margin for the characterization of load node respectively, which can be the supplement to the on-line voltage stability index of node (system). If the on-line voltage stability index of load note approaches to 1, it shows that the power system draws near voltage collapse. Simulation results of EPRI-36 bus system show that the proposed method can effectively judge the system voltage stability, and is suitable for the monitoring and prediction of on-line voltage stability. Compared with other indexes, the rationality and superiority of the proposed model and index are reflected. Finally, we briefly explain the feasibility of the proposed method through the voltage stability margin evaluation of load node in Shandong 500 kV main network.

voltage stability; on-line voltage stability index; phasor measurement unit(PMU); network equivalent model

“十二五”國家科技支撐計劃重大項目(2012BAK04B00)

TM 712

A

1000-7229(2016)02-0069-09

10.3969/j.issn.1000-7229.2016.02.010

2015-10-27

高峰(1988)，男，博士研究生，主要研究方向為廣域同步測量技術和電力系統分析；

王彥文(1962)，男，博士，教授，主要研究方向為廣域同步測量系統、網絡型繼電保護技術等；

郭林(1986)，男，博士研究生，主要研究方向為網絡型繼電保護技術；

徐睿(1989)，男，碩士，工程師，主要研究方向為電網安全運行及電網調度控制。

Project supported by key Project of the National Twelfth-Five Year Research Program of China((2012BAK04B00))