電力系統(tǒng)暫態(tài)穩(wěn)定直接法的前沿追蹤與述評

劉楊帆，楊金成，劉艷芳

(三峽大學(xué)電氣與新能源學(xué)院，湖北 宜昌 443000)

電力系統(tǒng)暫態(tài)穩(wěn)定直接法的前沿追蹤與述評

劉楊帆，楊金成，劉艷芳

(三峽大學(xué)電氣與新能源學(xué)院，湖北 宜昌 443000)

穩(wěn)定破壞是電網(wǎng)中較為嚴(yán)重的事故之一。為了防止穩(wěn)定破壞，直接法作為一種快速、精確的穩(wěn)定算法，已經(jīng)在電力系統(tǒng)穩(wěn)定分析中取得了豐碩的成果?；仡櫫酥苯臃☉?yīng)用于暫態(tài)穩(wěn)定分析中的歷史和現(xiàn)狀。在介紹直接法的一些主要方法的基礎(chǔ)上，對直接法用于電力系統(tǒng)暫態(tài)穩(wěn)定分析作了簡要的評價，最后對直接法的發(fā)展進行了展望。

暫態(tài)穩(wěn)定分析；直接法；歷史及現(xiàn)狀；述評與展望

1 引言

電力系統(tǒng)暫態(tài)穩(wěn)定性問題是指電力系統(tǒng)在某個運行情況下突然受到大的干擾(各種短路故障，切除大容量發(fā)電機或輸電設(shè)備以及某些負荷的突然變化等)后，能否經(jīng)過暫態(tài)過程達到新的穩(wěn)態(tài)運行狀態(tài)或者恢復(fù)到原來運行狀態(tài)的能力[1]。電力系統(tǒng)規(guī)模的迅速發(fā)展使得互聯(lián)電力網(wǎng)絡(luò)變得越來越復(fù)雜，隨之而帶來的后果是暫態(tài)穩(wěn)定問題變得日趨嚴(yán)重。統(tǒng)計表明，在電力系統(tǒng)穩(wěn)定破壞事故中，暫態(tài)穩(wěn)定破壞的事故是最嚴(yán)重的。因此，研究電力系統(tǒng)暫態(tài)穩(wěn)定分析的方法，對于提高系統(tǒng)運行的安全和穩(wěn)定性具有十分重要的意義。

求解電力系統(tǒng)暫態(tài)穩(wěn)定問題的方法主要有三類：時域仿真法、直接法、其它派生方法(例如：數(shù)據(jù)挖掘技術(shù)、向量分析法、人工智能法等)。本文重點就電力系統(tǒng)暫態(tài)穩(wěn)定分析中的直接法進行簡要介紹與評述。

2 直接法暫態(tài)穩(wěn)定分析的研究與發(fā)展

李雅普諾夫直接法(簡稱直接法)是從一個古典的力學(xué)概念發(fā)展而來的。該概念指出：“對于一個自由的(無外力作用)動態(tài)系統(tǒng)，若系統(tǒng)的總能量隨時間的變化率恒為負，則系統(tǒng)總能量不斷減少直至最終達到一個最小值，即平衡狀態(tài)，此時系統(tǒng)是穩(wěn)定的”[2]。直接法從系統(tǒng)能量及其轉(zhuǎn)化的角度分析穩(wěn)定問題，基本方法是確定故障后的暫態(tài)穩(wěn)定域，通過判斷故障后穩(wěn)定平衡點與故障清除瞬間系統(tǒng)的狀態(tài)是否處于該穩(wěn)定域中，可得到系統(tǒng)的運動軌跡是否是穩(wěn)定的。

2.1 直接法的發(fā)展歷史

李雅普諾夫在1892年提出了直接法穩(wěn)定性理論，但是直到1947年才由Magnusson提出應(yīng)用李雅普諾夫能量函數(shù)研究電力系統(tǒng)暫態(tài)穩(wěn)定。早期對直接法的研究主要集中在運用數(shù)學(xué)方法構(gòu)造李雅普諾夫能量函數(shù)和如何求取不穩(wěn)定平衡點(UEP)，在此基礎(chǔ)上的一些研究成果對系統(tǒng)中的結(jié)構(gòu)大量簡化，導(dǎo)致結(jié)果保守性和誤差均較大，從而直接法的發(fā)展一度處于低潮。

直到20世紀(jì)70年代之后，Athay和Kakimoto等人通過構(gòu)造暫態(tài)能量函數(shù)應(yīng)用于穩(wěn)定域的確定，為直接法的發(fā)展打開了一個新的局面。隨后，不斷有學(xué)者對直接法的應(yīng)用進行補充和完善，包括單機能量函數(shù)法的提出；應(yīng)用穩(wěn)定域概念對勢能界面法進行分析；以及在穩(wěn)定域基礎(chǔ)上BCU法的提出等等。

我國電力科學(xué)界對應(yīng)用直接法分析暫態(tài)穩(wěn)定性的研究大致始于20世紀(jì)80年代，在1983年夏道止與Heydt等人對分解-聚合法在線穩(wěn)定理論進行了研究。1987年清華大學(xué)教授倪以信與美國Fouad等人對不穩(wěn)定平衡點法的勵磁模型與直流輸電模型進行了研究；其后，關(guān)于PEBS法、EEAC法、高階Taylor級數(shù)研究快速暫態(tài)穩(wěn)定計算問題、PEBS法復(fù)雜模型的研究等，使得直接法的應(yīng)用越來越走向?qū)嵱没?/p>

從直接法被提出用于電力系統(tǒng)暫態(tài)穩(wěn)定分析以來，已經(jīng)有了一百多年的研究歷史，在不斷發(fā)展的同時產(chǎn)生了大量的研究成果?？梢钥闯觯簢鴥?nèi)外學(xué)者們對直接法的不斷改進和完善(例如：改進現(xiàn)有模型，尋求更精確、快速和實用化的方法；充分集成各種方法之長，使其完美的有機結(jié)合[3])，目的是使得直接法更好地服務(wù)于電力系統(tǒng)暫態(tài)穩(wěn)定的應(yīng)用。

2.2 直接法的分支結(jié)構(gòu)

直接法的分支結(jié)構(gòu)如圖1所示。

圖1

3 經(jīng)典模型的能量函數(shù)法

能量函數(shù)法是在直接法的基礎(chǔ)上發(fā)展起來的，它利用暫態(tài)能量判斷系統(tǒng)的暫態(tài)穩(wěn)定性，思路與直接法一脈相承。在研究人員對借助能量函數(shù)分析暫態(tài)穩(wěn)定進行大量研究后，使得直接法慢慢轉(zhuǎn)變?yōu)闀簯B(tài)能量函數(shù)法。從系統(tǒng)失穩(wěn)時對暫態(tài)能量認識不同的角度出發(fā)，能量函數(shù)法可以分為全部能量函數(shù)法和局部能量函數(shù)法。

利用暫態(tài)能量函數(shù)法進行電力系統(tǒng)暫態(tài)穩(wěn)定分析，就是解決如下兩個基本問題：(1)構(gòu)造系統(tǒng)的暫態(tài)能量函數(shù)；(2)求取系統(tǒng)暫態(tài)能量函數(shù)值。

因此，一般包括如下四個基本步驟：

(1)首先，定義適當(dāng)?shù)臅簯B(tài)能量函數(shù)V(X)；

(2)計算故障清除時刻系統(tǒng)暫態(tài)能量函數(shù)V(X)的值Vcr；

(3)尋找對應(yīng)故障類型的暫態(tài)能量函數(shù)V(X)的臨界值Vcl；

(4)對故障后的微分方程從0～tcr進行定積分，使積分結(jié)果與暫態(tài)能量臨界值Vcr相等，從而求得故障臨界切除時間tcr。

當(dāng)Vcl

規(guī)格化的暫態(tài)穩(wěn)定域度ΔVn定義為：

當(dāng)ΔVn>0時，說明受擾后系統(tǒng)是暫態(tài)穩(wěn)定的。

3.1 全部能量函數(shù)法

全部能量函數(shù)法的觀點考慮了故障發(fā)生時系統(tǒng)中所有發(fā)電機的能量，即認為：機組失穩(wěn)與此時系統(tǒng)中包含的全部暫態(tài)能量有關(guān)，不僅包括分離出去的失穩(wěn)機組的能量，還有非失穩(wěn)機組之間相互搖擺的能量。

3.1.1 RUEP法

相關(guān)不穩(wěn)定平衡點法(RUEP)是在不穩(wěn)定平衡點法(UEP)的研究過程中發(fā)展起來的，學(xué)者通過對動能修正、相關(guān)求解以及能量裕度等方面的進一步發(fā)展和豐富，使人們對系統(tǒng)失穩(wěn)模式，故障軌跡及RUEP等一些概念之間的相互關(guān)系有了更清晰的認識。自從直接法借助能量的觀點分析系統(tǒng)暫態(tài)穩(wěn)定性以來，大量的研究無不建立在RUEP法的基礎(chǔ)之上。Athay[5]等人在研究中引入了慣性中心坐標(biāo)，認為機組失穩(wěn)與否與慣性中心的動能無關(guān)，在計算時需要將其扣除。同時在研究中，作者使用一種優(yōu)化搜索方法，即用系統(tǒng)近似故障軌跡代替實際受擾軌跡，從勢能最大點出發(fā)，尋找到相關(guān)不穩(wěn)定平衡點處的勢能就可以作為臨界能量。

隨后，F(xiàn)ouad對RUEP法進一步研究后指出，若故障穩(wěn)定平衡點與故障前系統(tǒng)穩(wěn)定時的平衡點相同，就需要對暫態(tài)能量的起始點做出一定的修正。但其實，由于一般在求解故障的臨界軌跡時，并不能順利的確定。因此，無論是Athay的優(yōu)化搜尋法，還是Fouad在進一步發(fā)展之后提出的確定主導(dǎo)不穩(wěn)定平衡點準(zhǔn)則都只是得到與首擺失穩(wěn)相關(guān)的不穩(wěn)定平衡點。

用RUEP法作暫態(tài)分析的步驟為：

(1)輸入原始數(shù)據(jù)，形成穩(wěn)態(tài)導(dǎo)納矩陣Y，計算故障前平衡點θS0；

(2)有故障或操作，修改Y；否則按時間步長仿真，用數(shù)值解法求解轉(zhuǎn)子運動方程，得到故障切除時間(wC，δC)后作為初值，求取故障后的平衡點δS；

(3)判斷失穩(wěn)模式并以(2)中求得的平衡點δS為初值，在慣性中心坐標(biāo)下用以下功率平衡方程計算θu。

(1)

3.1.2 勢能界面法(PEBS)

勢能界面法(PEBS)是暫態(tài)能量函數(shù)法的一個分支，全稱為位能邊界曲面方法。它于1978年由日本學(xué)者Kakwi Oto等人提出，后來Athay及其他一些學(xué)者對PEBS法的發(fā)展做出了重要貢獻。文獻[6]對PEBS法的理論基礎(chǔ)給出了詳細介紹，指出PEBS是電力系統(tǒng)中的相關(guān)梯度系統(tǒng)的穩(wěn)定邊界。在系統(tǒng)遭受大擾動時，PEBS法通過搜索持續(xù)故障軌跡和勢能邊界的交點，將其作為臨界能量估計值，確定系統(tǒng)的穩(wěn)定邊界和穩(wěn)定性?；舅枷胧怯霉收蠒r持續(xù)故障軌線來代替臨界故障軌跡，這也同時避免了計算主導(dǎo)不穩(wěn)定平衡點(CUEP)的計算。

自從PEBS法提出之后，能量函數(shù)法在理論和實用化方面都有了很大的發(fā)展。但該方法在應(yīng)用中也避免不了其局限性，例如：如果系統(tǒng)故障時與故障切除后結(jié)構(gòu)差別比較大，導(dǎo)致故障軌跡不能代替臨界軌跡，這時顯然單獨使用PEBS法是解決不了問題的；另外，計算過程中轉(zhuǎn)子角度會時間變化，這就需要把握計算速度問題。因此，長期以來，對PEBS法的研究不是孤立的，它常常與其它暫態(tài)分析方法相互補充，相互滲透。3.1.3 基于穩(wěn)定域邊界的主導(dǎo)不穩(wěn)定平衡點法(BCU)

BCU法是建立在現(xiàn)代動力學(xué)穩(wěn)定理論基礎(chǔ)之上的一種新的暫態(tài)穩(wěn)定域計算方法[3]。它借助梯度系統(tǒng)的CUEP尋找故障后的不穩(wěn)定平衡點，并以該點處的暫態(tài)能量作為近似系統(tǒng)的局部穩(wěn)定域邊界。BCU法搜索主導(dǎo)不穩(wěn)定平衡點的示意圖見圖2。

圖2

該方法基于梯度系統(tǒng)關(guān)于穩(wěn)定域邊界關(guān)系的理論，求解步驟基本上可以歸結(jié)如下：(1)消去系統(tǒng)中除發(fā)電機外的所有內(nèi)節(jié)點，建立梯度系統(tǒng)中發(fā)電機內(nèi)節(jié)點的模型；(2)通過逐步仿真，由積分方程得到故障軌跡，檢測出口點EP(故障軌跡投影與收縮系統(tǒng)邊界的交點)；(3)從出口點出發(fā)，對事故梯度系統(tǒng)方程積分，并且沿積分軌跡尋找最小梯度點GMP；(4)將最小梯度點作為初值，求解梯度系統(tǒng)的功率偏差方程，進而得到主導(dǎo)不穩(wěn)定平衡點CUEP；(5)最終求解系統(tǒng)臨界暫態(tài)能量。

3.2 局部能量函數(shù)法

局部能量函數(shù)法提出在系統(tǒng)發(fā)生故障的暫態(tài)過程中，對暫態(tài)能量起決定作用的只有失穩(wěn)機組的的能量，并非系統(tǒng)的全部能量，因為如果計入非失穩(wěn)機組相互搖擺的能量，不僅計算量大，而且計算結(jié)果也相對不夠準(zhǔn)確。

3.2.1 單機能量函數(shù)法(IMEF)

Michel等人在20世紀(jì)80年代之后提出了單機能量法。IMEF法從系統(tǒng)結(jié)構(gòu)內(nèi)部進行了分析，它通過求解在暫態(tài)過程中對穩(wěn)定性起決定作用的失穩(wěn)機組能量變化來確定機組在暫態(tài)過程中的行為。文獻[7]作者采用慣性中心坐標(biāo)，借助首次積分法推導(dǎo)出單臺機組的能量函數(shù)。

(2)

在系統(tǒng)是單機失穩(wěn)模式時應(yīng)用IMEF法計算結(jié)果較為準(zhǔn)確，但當(dāng)系統(tǒng)出現(xiàn)多臺發(fā)電機失穩(wěn)時，需要對失穩(wěn)電機逐個進行失穩(wěn)計算，不僅耗時多，而且此時的能量又包括了失穩(wěn)機群內(nèi)部的相對動能，故計算結(jié)果偏于保守。

3.2.2 擴展等面積法(EEAC)

擴展等面積法(EEAC)是由我國學(xué)者薛禹勝教授于1986年提出來的，此方法將現(xiàn)代電力系統(tǒng)暫態(tài)穩(wěn)定分析的思想與古典的等面積定則結(jié)合，具有快速、直觀的特點[9]。EEAC法的發(fā)展經(jīng)歷了SEEAC法、DEEAC法和IEEAC法三個階段。

SEEAC法是最早被提出來的，其采用經(jīng)典模型的多機系統(tǒng)，并借助互補群模型聚合技術(shù)將多機系統(tǒng)聚合到單機無窮大系統(tǒng)中，以得到解析的動態(tài)方程，最終利用等面積法則對系統(tǒng)的穩(wěn)定性進行判別。該方法在故障掃查、排序及參數(shù)極限值的預(yù)估方面非常有效。因此，SEEAC法對小型經(jīng)典模型的研究產(chǎn)生了大量成果。但是，它在應(yīng)用于規(guī)模較大的模型時，強烈的非同調(diào)性使得此時SEEAC法的結(jié)果無法滿足精度要求。也正是為了滿足工程精度要求，國內(nèi)和國際上的學(xué)者相繼提出了DEEAC法。這種方法巧妙地計及了能量在機群內(nèi)部和機群之間的轉(zhuǎn)換，通過計算精度和魯棒性的改善，保持了EEAC法的快速性。1993年7月薛禹勝教授再次綜合了SBS法和EEAC法，提出IEEAC法，利用SBS法得到的發(fā)電機動態(tài)軌跡代替DEEAC中分段大步長的故障切除前后的軌跡，使得積分空間和觀察空間得以分離。至此，可以有效地用EEAC法來處理多擺穩(wěn)定問題。

EEAC法的突出優(yōu)點是計算速度快，而且重新定義了穩(wěn)定裕度的概念，使用該方法前提是需要正確識別臨界機群。經(jīng)過長時間的探索，它已成為穩(wěn)定研究的理論武器和系統(tǒng)穩(wěn)定控制決策的實用工具。

3.2.3 時間尺度解耦法(TSD)

TSD法的提出是因為系統(tǒng)有如下特征：當(dāng)系統(tǒng)受到擾動時，就會分解為多組機群，而每群內(nèi)機組之間的功角差只是在一個微小范圍內(nèi)變化，因而此時不考慮群內(nèi)機組之間的振蕩對系統(tǒng)失穩(wěn)的影響，只是認為群間機組慣性中心的運動對系統(tǒng)臨界失穩(wěn)有較大影響?；诖耍涂梢园严到y(tǒng)分為兩個分量：快速和慢速分量。其中，快與慢分別對應(yīng)于同一機群和不同機群內(nèi)機組間的相對運動。因為這兩個分量之間耦合很小，可以解耦化為局部模式和區(qū)域模式。

TSD經(jīng)過進不斷的發(fā)展，形成了基于模態(tài)的能量函數(shù)法，該方法也已被應(yīng)用于大規(guī)模重載系統(tǒng)中，并且取得了很好的效果。

4 直接法分析電力系統(tǒng)暫態(tài)穩(wěn)定性的評價

4.1 直接法是基于穩(wěn)定理論支持的分析方法

以穩(wěn)定理論為基礎(chǔ)的直接法為電力系統(tǒng)暫態(tài)穩(wěn)定分析開辟了一個新的領(lǐng)域。直接法從能量的角度對系統(tǒng)暫態(tài)穩(wěn)定性進行分析，它具備了用域的方法定量分析穩(wěn)定程度的獨特優(yōu)勢。

4.2 直接法是一般穩(wěn)定分析的輔助方法

直接法可以給出穩(wěn)定度的定量值，并且提供穩(wěn)定裕度對運行條件變化的靈敏度分析，它不需要求解大規(guī)模非線性方程組，其計算速度非常快，從而可以滿足在線穩(wěn)定分析的要求。但考慮到運用直接法時模型的適應(yīng)性相對較差，精度和可靠性方面也只是用于第一擺暫態(tài)穩(wěn)定分析，因此直接法一般是作為暫態(tài)穩(wěn)定分析中的分類器，成為暫態(tài)穩(wěn)定分析中的一種輔助方法。

5 展望

通過以上對直接法的介紹，可以看出直接法的研究不是孤立的，對屬于直接法的各類方法，在應(yīng)用中都發(fā)現(xiàn)了其缺點(例如計算精度、計算速度、模型適應(yīng)性等方面)。因此，在直接法的發(fā)展過程中，一些方法常常相互滲透，互相補充。另外，時域仿真法和直接法結(jié)合的“混合法”也已經(jīng)被應(yīng)用于電力系統(tǒng)暫態(tài)穩(wěn)定分析中，并且效果良好。無論如何，減少對模型的依賴性，充分應(yīng)用穩(wěn)定理論給出所要求的定量判據(jù)并對系統(tǒng)的穩(wěn)定性做出正確評估，同時減少計算時間是人們對直接法的期望。相信隨著直接法的逐步成熟，將會為電力系統(tǒng)快速在線安全穩(wěn)定分析提供一個強有力的工具。

[1] 潘艷霞.電力系統(tǒng)暫態(tài)穩(wěn)定直接法的分析[D].太原理工大學(xué),2004.

[2] 傅書逷,倪以信,薛禹勝.直接法穩(wěn)定分析[M].北京：中國電力出版社,1999:39-77.

[3] 阮青松.基于能量函數(shù)的電力系統(tǒng)暫態(tài)穩(wěn)定性分析[D].西南交通大學(xué),2006.

[4] Chiang H D.Wu FE,Veraiya P.Foundation of Direct Methods for Power System Transient Stability Analysis[J].IEEE Trans.On Circuits and System.1987,CAS-34(2):160-173,712,712-728.

[5] T.Athay,R.Podmore,Virmani.A Practical Method for Direct Analysis of Transient Stability.IEEE Trans.on Power Apparatus and Systems,Vol.PAS-98,NO.2,1979.

[6] Chiang.H.D.Wu.F.F.Varaiya.P.P.Foundation of the Potential Energy Boundary Surface Method for Power System Transient Stability Analysis[J].IEEE Trans.on Circuits and System.1998,CAS-35(6);712-718.

[7] Y.Xue,et al.A Simple Direct Method for Fast Transient Stability Assessment of Large Power System,IEEETrans.On Power System,1983,3(5):pp266-276

[8] Michel A N.Power System Transient Stability Using Individual Machine Energy Function[J].Trans.On Circuits and Systems,1983,30(5):266-276.

[9] 薛禹勝.EEAC與直接法的機理比較(一)——受擾程度函數(shù)[J].電力系統(tǒng)自動化,2001,11:6-11.

[10] Han Xiaoyan,Han Zhenxiang,Research on intrinsic parallel algorithm of power system transient analysis[J].Proceedings of CSEE;1997;03:2-5(in Chinese).

[11] Xue Y,Rousseaux P,Gao Z,et al.Dynamic Extended Equal Area Criterion.Part 1.Basic Formulation in:Athens Power Tech,1993.

[12] 尹孔殷.電力系統(tǒng)暫態(tài)穩(wěn)定性分析的新方法與直接法[J].電網(wǎng)技術(shù),1995,1:58-61.

[13] 傅書逖.勢能界面法(PEBS)暫態(tài)穩(wěn)定分析的綜述及展望[J].電力系統(tǒng)自動化,1998,9:16-18.

[14] 張妍,王治,周寧.對于電力系統(tǒng)暫態(tài)穩(wěn)定分析的幾點認識[J].江西電力,2011,4:18-20.

Frontier Tracing and Review of the Direct Method for Power System Transient Stability

LIUYang-fan，YANGJin-cheng,LIUYan-fang

(School of Electrical and New Energy,China Three Gorges University,Yichang,443000 Hubei)

Stability failure is one of the more serious accidents in the power grid.In order to prevent the destruction of stability,direct method has been used to analyze the stability of power system,because it is a fast and accurate algorithm.The paper reviews the history and current situation of direct method in transient stability analysis.On the basis of some main methods of the direct method,it briefly evaluates the direct method for transient stability analysis of the power system.Finally,the development of the direct method is prospected.

transient stability;direct method;history and present situation;review and prospect

1004-289X(2016)03-0011-05

TM71

B

2015-11-13