基于連續(xù)潮流的靜態(tài)分岔點追蹤方法

李林劉玲

(德陽電業(yè)局變電檢修中心，四川德陽 618000)

電力市場化改革以后，經(jīng)濟的快速增長，投入系統(tǒng)的設(shè)備增長迅猛，同時電力網(wǎng)絡(luò)的發(fā)展相對滯后，很多地方的電網(wǎng)都經(jīng)常運行在極限狀態(tài)，可能遭受電壓穩(wěn)定的威脅［1］。靜態(tài)電壓穩(wěn)定分析能判斷電力系統(tǒng)有多大的靜態(tài)電壓穩(wěn)定裕度并指出采取何種合適的控制措施使電力系統(tǒng)運行在安全裕度之內(nèi)［2］。

在電力系統(tǒng)靜態(tài)電壓穩(wěn)定性分析中，電壓崩潰點是與靜態(tài)分岔點相聯(lián)系的。電力系統(tǒng)中主要有兩種靜態(tài)分岔點:鞍結(jié)分岔點(saddle－node bifurcation point，SNBP)和極限誘導分岔點(limit－induced bifurcation point，LIBP)。前者對應(yīng)潮流雅可比矩陣的奇異點;后者則是由于某發(fā)電機節(jié)點的無功出力達到極限而導致系統(tǒng)穩(wěn)定平衡點忽然消失［3］。因此，只有準確搜索和識別這兩類分岔點，才能準確地計算系統(tǒng)的電壓穩(wěn)定裕度，進而準確地計算在該分岔點處各種控制變量對于穩(wěn)定裕度的靈敏度信息，從而進行預(yù)防控制［4］。求取電壓靜態(tài)分岔點的算法主要有連續(xù)潮流法［5－7］、崩潰點直接計算法［8－10］、最優(yōu)潮流計算法［11－13］。

文獻［4］利用連續(xù)潮流法和崩潰點法的混合方法對極限誘導分岔點和鞍結(jié)分岔點進行了識別和計算。文獻［14］利用預(yù)測－校正方法快速穿越電壓崩潰臨界點的方法來搜索。文獻［15］采用觸發(fā)條件和兩個檢測條件的判斷準則，對結(jié)構(gòu)變化點進行檢測，從而判斷極限誘導分岔點的存在。這3篇文獻中都提到利用最優(yōu)步長的方法來計算，但是都沒有告訴具體的方法。這里采用發(fā)電機無功極限值引導變步長的連續(xù)潮流法來搜索靜態(tài)分岔點，繼而利用雅可比矩陣特征值的符號變化進行識別。該方法應(yīng)用于IEEE 39節(jié)點測試系統(tǒng)，顯示了該方法的正確性和有效性。

1 連續(xù)潮流的數(shù)學模型

一般，在靜態(tài)輸電計算中，極坐標系統(tǒng)下的潮流方程可用式(1)表示為

式中，PGi、QGi為節(jié)點 i的發(fā)電機有功和無功出力;Qmin，i、Qmax，i為發(fā)電機的最小和最大無功極限值;PLi、QLi為節(jié)點i負荷的有功和無功;Vi、Vsi是節(jié)點i的電壓幅值和節(jié)點i的標準電壓幅值;θij是節(jié)點i和節(jié)點j的電壓相角差值;Gij、Bij為節(jié)點i與節(jié)點j之間的網(wǎng)絡(luò)導納矩陣的實部和虛部。

若用λ表示發(fā)電機和負荷的增長參數(shù)，即為負荷因子;nGi、nPLi分別表示發(fā)電機和負荷有功增長的方向向量;nQLi表示負荷無功增長的方向向量。則有

由于引入變量λ，使潮流方程的未知數(shù)比方程數(shù)多一個，為了使未知數(shù)的個數(shù)和方程數(shù)相同，避免潮流雅可比矩陣奇異，此時需要增加一維參數(shù)化方程。

連續(xù)潮流主要由4部分組成:預(yù)估、校正、參數(shù)化和步長控制。具體步驟不再敘述，詳情請參考文獻［16］、［17］。

2 發(fā)電機無功受限后現(xiàn)象分析

在實際潮流計算中，發(fā)電機的無功并不是不變的。當發(fā)電機的無功始終維持在最小值與最大值之間范圍內(nèi)時，潮流雅可比矩陣在最大功率傳輸點奇異，其靜態(tài)分岔類型為鞍結(jié)分岔(saddle－node bifurcation，SNB);當發(fā)電機的無功越限時，系統(tǒng)中PV節(jié)點的無功功率達到極限將保持不變，而節(jié)點電壓發(fā)生跌落現(xiàn)象，在潮流方程中發(fā)電機節(jié)點由PV節(jié)點類型轉(zhuǎn)化為PQ節(jié)點類型，潮流方程中QGi取值如下。

若 QGi≤ Qmin，i時，QGi=Qmin，i，Vi≤Vsi若 QGi≥ Qmax，i時，QGi=Qmax，i，Vi≤Vsi

此時，系統(tǒng)可能發(fā)生兩種現(xiàn)象:無功/電壓約束轉(zhuǎn)換(constraint exchange pomt，CEP)和極限誘導分岔(Gmit－induead bifurcation，LIB)。

2.1 無功/電壓約束轉(zhuǎn)換

CEP對應(yīng)的分岔圖如圖1所示。圖中，點A為系統(tǒng)初始運行點，此時發(fā)電機無功沒有越限。隨著負荷的增加，當發(fā)電機無功達到極限(B點)時，系統(tǒng)的結(jié)構(gòu)發(fā)生突然變化，運行點從QGi≤QGmaxi曲線的上半分支轉(zhuǎn)化到Q=QGmaxi曲線的上半分支運行，此時系統(tǒng)是穩(wěn)定的。如果再增加負荷，系統(tǒng)將運行到C點，就會發(fā)生SNB而導致電壓崩潰。對于B點來講，是一個

圖1 無功/電壓約束轉(zhuǎn)換

2.2 極限誘導分岔

LIB對應(yīng)的分岔圖如圖2所示。圖中，點A為系統(tǒng)初始運行點，此時發(fā)電機無功沒有越限。隨著系統(tǒng)負荷的增加，當發(fā)電機無功達到極限(B點)時，運行點突然從QGi≤QGmaxi曲線的上半分支轉(zhuǎn)化到QGi=QGmaxi曲線的下半分支運行，形成一個不同于C點(SNBP)的尖點形狀，系統(tǒng)的結(jié)構(gòu)突然改變導致了電壓崩潰。對于B點來講，是一個不穩(wěn)定的運行點，這就是LIBP。

圖2 極限誘導分岔

3 靜態(tài)分岔點的追蹤

3.1 靜態(tài)分岔點的搜索

搜索方法通常用常規(guī)連續(xù)潮流法，從起始點開始逐個搜索。這樣的話，一方面降低了算法的效率，另一方面，如果連續(xù)潮流法的步長過大的話，多臺發(fā)電機都穿越了分岔點。因此，要用盡可能少的步數(shù)達到和穿越分岔點，從而可以把分岔點限制在一個較小的范圍內(nèi)，為精確識別做好準備。

于是采用發(fā)電機無功極限值引導變步長連續(xù)潮流法來搜索分岔點。當發(fā)電機無功沒有越限時，采用較大步長的非線性預(yù)估，當發(fā)電機快要達到無功極限時，采用較小的步長，避免多臺發(fā)電機在相鄰短時間內(nèi)先后越限。具體步長控制策略如下。

式(1)中發(fā)電機的無功功率QGi可以如下描述。

式中，Yij和αij是導納矩陣在極坐標下的元素。對于一個特定的電力系統(tǒng)，導納矩陣是一定的，所以Yij和αij是個常數(shù)。對于發(fā)電機節(jié)點來說Vi也是個常數(shù)，如果電力系統(tǒng)在正常條件下運行，相角差θij通常很小，所以，式(3)中的正弦函數(shù)接近于常數(shù)，描述如下。

從式(3)中可以看出電壓幅值Vj是唯一的變量，即發(fā)電機的無功功率QGi是關(guān)于電壓幅值Vj的函數(shù)，這樣可以簡化式(3)如下。

式中，K=［k1，k2，k3，…kn］是線性系數(shù)矩陣，V=［V1，V2，V3，…Vn］T是系統(tǒng)的電壓幅值矩陣。從式(5)可以看出Vi與QGi成線性關(guān)系。負荷因子λ與電壓幅值Vi接近于平方關(guān)系［18］，所以負荷因子λ與發(fā)電機的無功功率QGi也接近于平方關(guān)系。關(guān)系示意圖如圖3所示，表達式描述如下。

圖3 負荷因子λ與發(fā)電機無功QGi的關(guān)系

如果知道式(6)的3個解，根據(jù)二階拉格朗日插值就可以得出系數(shù)ai、bi和ci。當發(fā)電機的無功到達極限，把上限QGimax和下限QGimin代入式(6)，可以得到發(fā)電機無功越限時相應(yīng)的負荷因子。

3.2 靜態(tài)分岔點的識別

分岔點可根據(jù)靜態(tài)分岔類型的特點進行識別。SNB的特點是，雅可比矩陣有零特征值出現(xiàn);CEP的特點是，在無功/電壓約束轉(zhuǎn)換過程的前后，雅可比矩陣的所有特征值仍在左開半平面;LIB的特點是，在發(fā)生分岔之前，雅可比矩陣的所有特征值仍在左開半平面，而在發(fā)生分岔之后，雅可比矩陣的最小特征值跳躍到右開半平面。

3.3 追蹤分岔點的步驟

1)采用連續(xù)潮流法計算出3個運行點，獲得λ與QGi值;

2)利用二階拉格朗日插值法獲取式(6)的系數(shù)ai、bi和 ci;

3)利用式(7)計算出發(fā)電機無功越限后的λq值，從而計算出連續(xù)潮流的步長h;

4)利用連續(xù)潮流法進行預(yù)估、校正，得到潮流解;

5)對潮流解進行識別，若不是分岔點就采用固定小步長進行搜索，執(zhí)行上一步;若是分岔點則搜索成功，把PV節(jié)點轉(zhuǎn)化為PQ節(jié)點，執(zhí)行第二步，進行下一個點的搜索。

4 算例分析

下面以IEEE 39節(jié)點測試系統(tǒng)對采用的方法進行了驗證。仿真中所有的負荷和發(fā)電機的有功按同一比例λ增長，負荷為恒功率模型，并且保持功率因數(shù)不變。發(fā)電機無功功率取值范圍如表1所示，其中Bus_38發(fā)電機的最大無功功率Qmax有兩個值，對應(yīng)兩種不同的模式。

表1 發(fā)電機無功功率取值范圍

模式1:Bus_38的 Qmax=1.5。Bus_6的 PV曲線如圖4所示。隨著負荷因子的增加，系統(tǒng)的初始運行曲線為曲線1，某一發(fā)電機到達無功受限點A后，系統(tǒng)的潮流方程發(fā)生了突變，運行曲線變化到曲線2。接著發(fā)電機無功受限點依次為B、C、D，運行曲線對應(yīng)依次變化為曲線3、4、5，隨著負荷因子的進一步增加，系統(tǒng)運行到E點，發(fā)生電壓崩潰。

圖4 模式1下Bus_6的PV曲線圖

圖4中，A、B、C、D、E 對應(yīng)的識別結(jié)果如表2所示。從表中可以看出在受限點A、B、C、D潮流雅可比矩陣特征值符號都是從負變化到負，即潮流雅可比矩陣的所有特征值仍在左開半平面，系統(tǒng)仍然穩(wěn)定。在受限點E負荷達到最大，潮流雅可比矩陣特征值從符號為負變化為零，此時對應(yīng)電壓穩(wěn)定臨界點。所以模式1下的電壓崩潰是由于鞍結(jié)分岔所致。

圖5 模式2下Bus_6的PV曲線圖

表2 模式1下測試系統(tǒng)的識別結(jié)果

模式2:Bus_38的Qmax=0.32。Bus_6的PV曲線如圖5所示。隨著負荷因子的增加，發(fā)電機的無功功率受限點依次為A、B、C、D、E，各點對應(yīng)的識別結(jié)果如表3所示。

表3 模式2下測試系統(tǒng)的識別結(jié)果

從圖5和表3可以看出在受限點A、B、C、D潮流雅可比矩陣特征值符號都是從負變化到負，即潮流雅可比矩陣的所有特征值仍在左開半平面，系統(tǒng)仍然穩(wěn)定。在受限點E負荷達到最大，潮流雅可比矩陣特征值符號從負變化到正，即在發(fā)生分岔之前，潮流雅可比矩陣的所有特征值仍在左開半平面，而在發(fā)生分岔之后，潮流雅可比矩陣的特征值跳躍到右開半平面，此時發(fā)生電壓崩潰。所以模式2下的電壓崩潰是由于極限誘導分岔所致。

5 結(jié)論

對發(fā)電機無功受限后出現(xiàn)的無功/電壓約束轉(zhuǎn)換和極限誘導分岔現(xiàn)象進行了分析，并采用了發(fā)電機無功極限值引導變步長的連續(xù)潮流法，對靜態(tài)電壓穩(wěn)定分析中兩類主要的分岔點進行了搜索，繼而利用潮流雅可比矩陣特征值的符號變化進行識別。通過IEEE 39節(jié)點測試系統(tǒng)的算例分析，取得了理想的效果，表明了該方法的正確性和有效性。本方法可以提供系統(tǒng)可能的電壓崩潰方式，這對電力系統(tǒng)運行人員進行預(yù)防控制具有非常重要的指導意義。

［1］王冰，游振華，韓學軍，等.考慮負荷靜態(tài)電壓特性的重負荷節(jié)點靜態(tài)電壓穩(wěn)定分析［J］.繼電器，2007，35(11):21－25.

［2］李娟，陳繼軍，司雙.連續(xù)潮流與免疫遺傳算法結(jié)合的靜態(tài)電壓穩(wěn)定裕度計算［J］.電力系統(tǒng)保護與控制，2010，38(18):24 －27.

［3］Canizares C.A.，Mithulananthan N.，Berizzi A.，et al.On the Linear Profile of Indices for the Prediction of Saddle－node and Limit－induced Bifurcation Points in Power Systems［J］.IEEE Trans on Circuits and Systems I:Fundamental Theory and Applications，2003，50(12):1588－1595.

［4］馬冠雄，劉明波，王奇.一種識別靜態(tài)電壓穩(wěn)定分岔點的混合方法［J］.電力系統(tǒng)自動化，2006，30(24):17 －20.

［5］Ajjarapu V.，Christy C.The Continuation Power Flow:A Tool for Steady State Voltage Stability Analysis［J］.IEEE Transactions on Power Systems，1992，7(1):416 –423.

［6］Feng Z.，Ajjarapu V.，Long B..Identification of Voltage Collapse Through Direct Equilibrium Tracing［J］.IEEE Transactions on Power Systems，2000，15(1):342 －349.

［7］Yorino N.，LI Huaqiang，，Sasaki H..A Predictor/corrector Scheme for Obtaining Q－limit Points for Power Flow Studies［J］.IEEE Trans on Power Systems，2005，20(1):130－137.

［8］Canizares C.A.，Alvarado F.L..Point of Collapse and Continuation Methods forLlarge AC/DC Systems［J］.IEEE Trans on Power Systems，1993，8(1):1 －8.

［9］郭瑞鵬，韓禎祥.電壓崩潰臨界點計算的改進零特征根法［J］.中國電機工程學報，2000，20(5):63－66.

［10］Rosehart W.D.，Canizares C.A.，Quintana V.H.Multiob jectiveOptimal Power Flows to Evaluate Voltage Security Cost in Power Networks［J］.IEEE Trans on Power Systems，2003，18(2):578 －587.

［11］Irisarri G.D.，Wang X.，Tong J.，et al.Maximum Load A-bility of Power Systems Using Interior Point Non－linear Optimization Method［J］.IEEE Trans on Power Systems，1997，12(1):162 －172.

［12］韋化，丁曉鶯.基于現(xiàn)代內(nèi)點理論的電壓穩(wěn)定臨界點算法［J］.中國電機工程學報，2002，22(3):27－31.

［13］李華強，劉亞梅，YORINO N.鞍結(jié)分岔與極限誘導分岔的電壓穩(wěn)定性評估［J］.中國電機工程學報，2005，25(24):56－60.

［14］趙晉泉，江曉東，張伯明.一種靜態(tài)電壓穩(wěn)定臨界點的識別和計算方法［J］.電力系統(tǒng)自動化，2004，28(23):28－32.

［15］熊寧，程浩忠，馬則良，等.基于連續(xù)潮流的極限誘導分岔檢測方法［J］.電力系統(tǒng)自動化，2008，32(18):35－38.

［16］徐曉春，萬秋蘭，顧偉，等.在線電壓穩(wěn)定分析中快速連續(xù)潮流的應(yīng)用［J］.電力系統(tǒng)保護與控制，2008，36(14):37－41.

［17］周雙喜，朱凌志，郭錫玖，等.電力系統(tǒng)電壓穩(wěn)定性及其控制［M］.北京:中國電力出版社，2004.

［18］Antonio C.Z.de Souza，Claudio A.Canizzares，Victor H.Quintana.New Techniques to Speed up Voltage Collapse Computations Using Tangent Vectors［J］.IEEE Trans on Power Systems，1997，12(3):1380 – 1387.