陽(yáng)育德， 李 雨,2， 袁 輝， 吳忠標(biāo)， 楊 健

(1. 廣西電力系統(tǒng)最優(yōu)化與節(jié)能技術(shù)重點(diǎn)實(shí)驗(yàn)室， 廣西大學(xué)， 廣西 南寧 530004;2. 廣西電網(wǎng)有限責(zé)任公司枊州供電局， 廣西 柳州 545005)

電力系統(tǒng)小干擾穩(wěn)定的時(shí)域計(jì)算及波形分析

陽(yáng)育德1， 李 雨1,2， 袁 輝1， 吳忠標(biāo)1， 楊 健1

(1. 廣西電力系統(tǒng)最優(yōu)化與節(jié)能技術(shù)重點(diǎn)實(shí)驗(yàn)室， 廣西大學(xué)， 廣西 南寧 530004;2. 廣西電網(wǎng)有限責(zé)任公司枊州供電局， 廣西 柳州 545005)

電力系統(tǒng)小干擾穩(wěn)定分析的實(shí)質(zhì)是系統(tǒng)近似線性化后所得的定常時(shí)不變系統(tǒng)的分析，系統(tǒng)的時(shí)域解特性使得仿真曲線的阻尼表征了系統(tǒng)的振蕩模態(tài)?；谶@一特性，本文提出一種基于時(shí)域仿真的小干擾穩(wěn)定分析方法。該方法通過(guò)直觀分析系統(tǒng)的時(shí)域仿真曲線的收斂程度，快速判定系統(tǒng)的小干擾穩(wěn)定性，無(wú)需對(duì)系統(tǒng)狀態(tài)矩陣進(jìn)行特征值求解，進(jìn)而避免了特征值求解方法存在的各種不足，為電力系統(tǒng)小干擾穩(wěn)定在線分析提供了一種新的依據(jù)。并對(duì)系統(tǒng)的仿真曲線采用HHT方法進(jìn)行波形分析，得到系統(tǒng)的最小阻尼比來(lái)體現(xiàn)系統(tǒng)穩(wěn)定裕度。最后，分別對(duì)WSCC3機(jī)9節(jié)點(diǎn)小系統(tǒng)及某地區(qū)實(shí)際大電網(wǎng)的某個(gè)時(shí)間點(diǎn)的典型運(yùn)行方式下的算例進(jìn)行仿真實(shí)驗(yàn)，驗(yàn)證所提方法的有效性。

電力系統(tǒng)； 振蕩模態(tài)； 時(shí)域仿真； 小干擾穩(wěn)定； 阻尼比； HHT

1 引言

隨著現(xiàn)代電力系統(tǒng)大區(qū)域互聯(lián)[1,2]、特高壓交直流電網(wǎng)跨區(qū)網(wǎng)架的形成[3,4]，電力系統(tǒng)小干擾穩(wěn)定成為影響互聯(lián)系統(tǒng)運(yùn)行、限制區(qū)域間傳輸功率的主要因素。同時(shí)清潔能源的接入和電動(dòng)汽車新市場(chǎng)的開(kāi)發(fā)也給系統(tǒng)的小干擾穩(wěn)定帶來(lái)很大的影響，大量風(fēng)電并網(wǎng)后的風(fēng)機(jī)類型、并網(wǎng)接入點(diǎn)數(shù)目、接入方式、并網(wǎng)容量比例等因素都對(duì)系統(tǒng)小干擾穩(wěn)定特性產(chǎn)生不同的影響[5]；電動(dòng)汽車需要大規(guī)模充電站，充電站的充電容量會(huì)對(duì)小干擾穩(wěn)定的阻尼比產(chǎn)生影響[6]。電力系統(tǒng)小干擾穩(wěn)定問(wèn)題[7]的研究越來(lái)越具有重要意義，也備受研究人員的關(guān)注。

目前，特征值計(jì)算是最為成熟的小干擾穩(wěn)定分析方法，其包括全部特征值算法和部分特征值算法。全部特征值算法是一類以實(shí)現(xiàn)實(shí)Schur分解為任務(wù)的GR算法，包括基于LU分解的LR迭代算法[8]、含帶狀上Hessenberg約化的BR迭代的BR算法[9]和QR算法等。QR法具有最穩(wěn)定的數(shù)值性能，魯棒性強(qiáng)，收斂速度快，但不能稀疏實(shí)現(xiàn)，容易產(chǎn)生維數(shù)災(zāi)，不能在大規(guī)模電力系統(tǒng)中得到廣泛應(yīng)用[10]。而針對(duì)大規(guī)模系統(tǒng)的特征值計(jì)算多采用部分特征值算法，文獻(xiàn)[11,12]應(yīng)用非Krylov子空間Jacobi-Davidson方法獲取子空間的最佳正交修正來(lái)逼近關(guān)鍵特征子集所對(duì)應(yīng)的特征空間，實(shí)現(xiàn)對(duì)特征值的自動(dòng)搜索，求取系統(tǒng)狀態(tài)矩陣的關(guān)鍵特征子集。另一類Krylov子空間的Arnoldi方法也是有效的特征值求解方法，隱式重啟Arnoldi 方法已被嵌入實(shí)用的小干擾穩(wěn)定分析工具中，但其收縮困難、收斂性差。文獻(xiàn)[13]應(yīng)用了擴(kuò)維重啟技術(shù)，既能鎖定已收斂的特征值，又可在此基礎(chǔ)上把特征值個(gè)數(shù)轉(zhuǎn)化成搜索半徑逐步擴(kuò)大搜索范圍，實(shí)現(xiàn)了不依賴設(shè)置參數(shù)自動(dòng)不漏根搜索，但是該方法在多次擴(kuò)維后殘差有增大的可能，需要多次非擴(kuò)維重啟才能在設(shè)定殘差之內(nèi)，增加了計(jì)算量，同時(shí)可能獲得多余特征值。文獻(xiàn)[11-13]的方法只求解影響穩(wěn)定性的關(guān)鍵特征值，以確保滿足大規(guī)模電力系統(tǒng)的計(jì)算效率的要求，但不能保證準(zhǔn)確地獲得系統(tǒng)所有的負(fù)阻尼和弱阻尼。

針對(duì)特征值計(jì)算存在的不足，文獻(xiàn)[14]避開(kāi)了對(duì)特征值的直接求解，通過(guò)狀態(tài)矩陣指數(shù)函數(shù)的跡曲線的振蕩信息反映系統(tǒng)的穩(wěn)定性。其應(yīng)用精細(xì)積分法計(jì)算矩陣指數(shù)函數(shù)，獲得系統(tǒng)線性化高階狀態(tài)矩陣的跡的數(shù)值曲線，通過(guò)西伯爾特-黃變換 (Hilbert-Huang Transform，HHT)方法分析曲線，獲取分析小干擾穩(wěn)定的阻尼比等模態(tài)參數(shù)。但精細(xì)積分方法在所求矩陣是滿陣、維數(shù)非常大的情況下，矩陣相乘計(jì)算量及存儲(chǔ)量隨維數(shù)迅速增長(zhǎng)；而且它的迭代次數(shù)對(duì)計(jì)算效率有著很大的影響，次數(shù)越高效率越低。

為了便于獲取表征系統(tǒng)振蕩的曲線，本文采用了時(shí)域仿真的方法。本文基于數(shù)值分析求解系統(tǒng)在小的擾動(dòng)后各狀態(tài)變量的時(shí)間響應(yīng)曲線，通過(guò)HHT方法對(duì)曲線分析計(jì)算出系統(tǒng)關(guān)鍵模態(tài)參數(shù)，如阻尼比、振蕩頻率等。

2 小干擾穩(wěn)定時(shí)域仿真模型

電力系統(tǒng)動(dòng)態(tài)特性都被描述為如下的微分-代數(shù)方程[10]：

(1)

式中，x為描述系統(tǒng)動(dòng)態(tài)特性的狀態(tài)變量；y為系統(tǒng)的運(yùn)行參量。本文選擇發(fā)電機(jī)六階模型，根據(jù)李雅普諾夫線性化理論，在系統(tǒng)平衡點(diǎn)(x0,y0)進(jìn)行線性化：

(2)

應(yīng)用隱式梯形積分公式：

(3)

對(duì)微分方程式(2)進(jìn)行差分化：

(4)

式中，dt表示仿真步長(zhǎng)；下標(biāo)n表示時(shí)段，n為當(dāng)前時(shí)段，則n+1則表示下一個(gè)時(shí)段。

3 HHT在小干擾穩(wěn)定分析中的應(yīng)用

HHT方法具有很強(qiáng)的自適應(yīng)性[15]，由兩部分組成[16，17]：經(jīng)驗(yàn)?zāi)Ｊ椒纸?Empirical Mode Decomposition，EMD)和Hilbert變換。對(duì)輸入信號(hào)進(jìn)行EMD分解獲得固有模式函數(shù)分量(Intrinsic Mode Function，IMF)，對(duì)每個(gè)IMF經(jīng)過(guò)Hilbert變化得到瞬時(shí)幅值、相角和頻率。

由EMD分解的具體步驟可知，對(duì)原始輸入信號(hào)s(t)可以分解為第一階IMF，記作c1(t)，第二階IMFc2(t)，第三階IMFc3(t)，……，繼而得到原始信號(hào)的IMF組成表達(dá)式：

(5)

式中，m為分解得到的IMF個(gè)數(shù)；r(t)為剩余量。

對(duì)任一經(jīng)EMD分解得到的IMF信號(hào)ci(t)，Hilbert變換為：

(6)

即可獲得解析信號(hào)：

(7)

從而得到信號(hào)的瞬時(shí)幅值a(t)、瞬時(shí)相位φ(t)和瞬時(shí)頻率f(t)的計(jì)算公式：

(8)

當(dāng)電力系統(tǒng)中的一個(gè)振蕩模式λ=σ±jω具有式(9)的表達(dá)形式：

(9)

則，一個(gè)振蕩信號(hào)是多個(gè)(如m個(gè))不同頻率的振蕩模式的線性和，系統(tǒng)振蕩曲線即可表達(dá)為：

(10)

當(dāng)對(duì)該曲線進(jìn)行EMD分解，則可獲得由多個(gè)IMF線性和組成的表達(dá)形式：

(11)

由式(10)和式(11)的對(duì)應(yīng)關(guān)系，可得到每個(gè)IMF對(duì)應(yīng)的振蕩模態(tài)參數(shù)的計(jì)算公式：

(12)

即

(13)

式中，fi,σi,ζi分別為原信號(hào)的第i個(gè)模態(tài)分量的頻率、衰減系數(shù)和阻尼比；fj,aj分別為對(duì)應(yīng)的IMF的瞬時(shí)頻率和瞬時(shí)幅值；Δt為時(shí)間間隔；p為選擇的時(shí)間間隔數(shù)。

4 算例分析

本文采用WSCC3機(jī)9節(jié)點(diǎn)系統(tǒng)[18]算例和某電網(wǎng)某個(gè)時(shí)刻的實(shí)際運(yùn)行方式(95臺(tái)發(fā)電機(jī)、615節(jié)點(diǎn)大系統(tǒng))算例進(jìn)行時(shí)域仿真。

4.1 WSCC3機(jī)9節(jié)點(diǎn)系統(tǒng)

該系統(tǒng)的結(jié)構(gòu)框架如圖1所示。電力系統(tǒng)運(yùn)行中，當(dāng)前的平衡點(diǎn)對(duì)于下一個(gè)運(yùn)行點(diǎn)來(lái)說(shuō)算是一個(gè)小的擾動(dòng)，本文設(shè)置這一擾動(dòng)為初始擾動(dòng)。仿真時(shí)間設(shè)為10s，針對(duì)系統(tǒng)三種情況：特征值都為負(fù)、有兩對(duì)實(shí)部為正的共軛特征值和兩個(gè)正特征值，分別進(jìn)行仿真計(jì)算。

圖1 WSCC3機(jī)9節(jié)點(diǎn)系統(tǒng)框架Fig.1 Framework of WSCC system with 3 machines and 9 buses

(1)系統(tǒng)的特征值都為負(fù)數(shù)的情況下，對(duì)系統(tǒng)的仿真結(jié)果與時(shí)域解的波形進(jìn)行比較，步長(zhǎng)為0.01s時(shí)，結(jié)果如圖2所示。三臺(tái)發(fā)電機(jī)的兩類波形的最大誤差點(diǎn)的誤差分別只有0.092°、0.244°、0.270°，都比較小，換算成誤差比分別為1.577%、0.3855%、0.4661%。

圖2 步長(zhǎng)為0.01s，仿真曲線與時(shí)域解波形比較Fig.2 Comparison of simulation curve and time domain solution wave in step of 0.01s

當(dāng)步長(zhǎng)取值為0.05s時(shí)，系統(tǒng)響應(yīng)曲線如圖3所示。可以看出，兩類波形有著非同步性，但曲線的整體趨勢(shì)一致。

圖3 步長(zhǎng)為0.05s，仿真曲線與時(shí)域解波形比較Fig.3 Comparison of simulation curve and time domain solution wave in step of 0.05s

具體誤差值分析表明，隨著仿真步長(zhǎng)的增大，誤差疊加效果比較明顯，兩類曲線的誤差也隨之增加。通過(guò)比較分析可知，在步長(zhǎng)為0.01s時(shí)產(chǎn)生的誤差只有2%，比較小，在可以接受范圍內(nèi)。因此后文中采用0.01s的步長(zhǎng)，對(duì)0.05s步長(zhǎng)的波形不再給予展示。

(2)對(duì)發(fā)電機(jī)參數(shù)進(jìn)行修改，三臺(tái)發(fā)電機(jī)的x″d參數(shù)的數(shù)值分別改為0.4pu、0.89pu、1.07pu，此時(shí)的特征值計(jì)算結(jié)果存在兩對(duì)實(shí)部為正的共軛特征值：λ1,2=0.7762±i10.183，λ3,4=0.905±i7.1001，此時(shí)的系統(tǒng)時(shí)域曲線如圖4所示。可以看出，兩類曲線都是振蕩發(fā)散的。

圖4 振蕩發(fā)散時(shí)仿真曲線與時(shí)域解波形比較Fig.4 Comparison of simulation curve and waveform of time domain solutions for oscillating divergence

(3)分析x″d參數(shù)設(shè)置過(guò)大的情況。x″d數(shù)值分別改為4pu、89pu、107pu，系統(tǒng)會(huì)出現(xiàn)兩個(gè)正的特征值：λ1=4.999，λ2=0.6801。存在兩個(gè)正特征值時(shí)系統(tǒng)仿真曲線及時(shí)域解波形如圖5所示。可以看出，仿真波形能較好地接近時(shí)域解波形，同時(shí)可以反映系統(tǒng)的不穩(wěn)定狀態(tài)。

圖5 直接發(fā)散時(shí)仿真曲線與時(shí)域解波形比較Fig.5 Comparison between simulation curve and waveform of time domain solutions when system is divergent

通過(guò)比較仿真曲線與時(shí)域解波形可知，曲線發(fā)散的趨勢(shì)是一致，同時(shí)仿真步長(zhǎng)足夠小時(shí)，兩收斂曲線能夠基本重合，誤差比較小，可以忽略不計(jì)，因此時(shí)域仿真可以應(yīng)用于電力系統(tǒng)的小干擾穩(wěn)定分析。同時(shí)表明，對(duì)系統(tǒng)穩(wěn)定性判據(jù)起主要作用的特征值與仿真曲線的收斂程度是相互對(duì)應(yīng)的，存在正的共軛特征值時(shí)曲線振蕩失穩(wěn)；存在正數(shù)特征值時(shí)曲線持續(xù)增大而發(fā)散。反之，通過(guò)曲線的收斂程度亦可以判定系統(tǒng)的小干擾穩(wěn)定，當(dāng)曲線發(fā)散，系統(tǒng)是小干擾不穩(wěn)定的；當(dāng)曲線收斂，系統(tǒng)是穩(wěn)定的。

4.2 某電網(wǎng)615節(jié)點(diǎn)系統(tǒng)

獲得某實(shí)際電網(wǎng)在某個(gè)時(shí)間點(diǎn)的運(yùn)行方式下的運(yùn)行數(shù)據(jù)及發(fā)電機(jī)參數(shù)，在發(fā)電機(jī)六階模型下的仿真結(jié)果如下：

(1)在特征值都為負(fù)值的情況下，仿真曲線與時(shí)域解波形比較如圖6所示。曲線收斂未發(fā)生發(fā)散現(xiàn)象，系統(tǒng)是穩(wěn)定的。

(2)當(dāng)改變運(yùn)行參數(shù)或發(fā)電機(jī)參數(shù)時(shí)，會(huì)使得特征值的正負(fù)情況發(fā)生改變。改變發(fā)電機(jī)參數(shù)后使得系統(tǒng)有一對(duì)實(shí)部為正的共軛特征值：λ1,2=0.3804±i14.388，獲得的仿真曲線與時(shí)域解波形如圖7所示?？梢钥闯?，波形振蕩失穩(wěn)。

圖6 系統(tǒng)特征值為負(fù)的仿真曲線 與時(shí)域解波形圖Fig.6 Time domain simulation curve of angle and waveforms of time domain solution with negative system’s eigenvalues

圖7 系統(tǒng)存在一對(duì)實(shí)部為正的共軛特征值 時(shí)仿真及時(shí)域解波形Fig.7 Time domain simulation curve of angle and waveforms of time domain solution when there are one couple of conjugate eigenvalues with positive real part

(3)改變發(fā)電機(jī)參數(shù)后使得系統(tǒng)存在兩個(gè)正特征值，分別是17.057，0.29062，獲得的仿真曲線與時(shí)域解波形如圖8所示。可以看出，在1s內(nèi)系統(tǒng)出現(xiàn)機(jī)組功角超過(guò)±180°，說(shuō)明存在正的特征值時(shí)，仿真曲線是發(fā)散的，系統(tǒng)不穩(wěn)定。

圖8 系統(tǒng)存在正特征值的時(shí)域仿真及時(shí)域解波形Fig.8 Time domain simulation curve of angle and waveforms of time domain solution with positive system’s eigenvalues

仿真結(jié)果表明，不論小系統(tǒng)還是實(shí)際大電網(wǎng)，當(dāng)步長(zhǎng)足夠小時(shí)，時(shí)域仿真波形都能很好地逼近時(shí)域解的波形，通過(guò)時(shí)域仿真曲線與特征值計(jì)算結(jié)果判定系統(tǒng)穩(wěn)定性所得結(jié)論一致。

4.3 算例的波形分析

4.1節(jié)和4.2節(jié)仿真了9節(jié)點(diǎn)和615節(jié)點(diǎn)系統(tǒng)穩(wěn)定狀態(tài)下的功角曲線，但是應(yīng)用曲線的收斂性來(lái)表達(dá)系統(tǒng)的穩(wěn)定時(shí)，并不能給出系統(tǒng)的穩(wěn)定裕度。為了克服這一不足，本文應(yīng)用HHT方法對(duì)曲線進(jìn)行分析獲得表征系統(tǒng)穩(wěn)定裕度的曲線阻尼參數(shù)。

(1)9節(jié)點(diǎn)系統(tǒng)

對(duì)9節(jié)點(diǎn)系統(tǒng)中發(fā)電機(jī)G1的功角仿真曲線進(jìn)行EMD分解，分解結(jié)果如圖9所示。

圖9 發(fā)電機(jī)G1的功角仿真曲線的經(jīng)驗(yàn)?zāi)Ｊ椒纸釬ig.9 EMD of simulation curves of power angle of generator1

圖9顯示發(fā)電機(jī)G1功角曲線被分解為1個(gè)IMF分量IMF1和剩余分量r。應(yīng)用式(8)求解IMF1的瞬時(shí)幅值、瞬時(shí)相位、瞬時(shí)頻率，并根據(jù)式(13)計(jì)算曲線的關(guān)鍵振蕩模態(tài)的頻率f、衰減系數(shù)σ和阻尼比ξ，結(jié)果見(jiàn)表1。

表1 各IMF分量參數(shù)Tab.1 Parameters of each IMF

對(duì)9節(jié)點(diǎn)系統(tǒng)中G2、G3兩臺(tái)機(jī)組的功角仿真曲線進(jìn)行EMD分解，對(duì)每個(gè)IMF分量進(jìn)行幅值、頻率、相角的計(jì)算，并換算為曲線的阻尼參數(shù)，結(jié)果如圖10和圖11所示，其中展示了表征曲線振蕩模式的參數(shù)。

圖10 發(fā)電機(jī)G2的功角仿真曲線及其參數(shù)Fig.10 Modal parameters of simulation curves of power angle of generator 2

圖11 發(fā)電機(jī)G3的功角仿真曲線及其參數(shù)Fig.11 Modal parameters of simulation curves of power angle of generator 3

分析結(jié)果表明，發(fā)電機(jī)G1、G2的仿真功角曲線只分解出一個(gè)IMF，該IMF的阻尼比則表征了仿真曲線的阻尼比。發(fā)電機(jī)G3的曲線分解得到IMF1和IMF2，說(shuō)明G3機(jī)組的功角曲線有兩個(gè)關(guān)鍵模態(tài)，其中最小阻尼比為0.0430。G2機(jī)組功角仿真曲線振蕩最為厲害，計(jì)算得到的阻尼比最小為0.0391，即為9節(jié)點(diǎn)系統(tǒng)的最小阻尼比，與穩(wěn)定狀態(tài)特征值計(jì)算得到的最小阻尼比0.04221比較接近。

(2)615節(jié)點(diǎn)系統(tǒng)

對(duì)于615節(jié)點(diǎn)大系統(tǒng)來(lái)說(shuō)，對(duì)每一臺(tái)發(fā)電機(jī)的功角曲線進(jìn)行分析的計(jì)算量很大，并且沒(méi)有必要，只需對(duì)其中振蕩最厲害的曲線進(jìn)行分析就可以分析系統(tǒng)的穩(wěn)定裕度情況。因此，對(duì)615節(jié)點(diǎn)系統(tǒng)的仿真曲線進(jìn)行判斷計(jì)算，得到振蕩最大的一條曲線為46號(hào)機(jī)組對(duì)應(yīng)的功角曲線，如圖12所示。對(duì)該曲線進(jìn)行分析和阻尼計(jì)算，結(jié)果如圖13所示，最小阻尼比為0.04380。

圖12 46號(hào)機(jī)組的時(shí)域仿真曲線Fig.12 Time domain simulation curve of generator 46

圖13 46號(hào)機(jī)組的時(shí)域仿真曲線的阻尼計(jì)算Fig.13 Damping calculation of time domain simulation curve of generator 46

兩個(gè)算例的波形分析說(shuō)明，通過(guò)曲線分析可獲得系統(tǒng)仿真曲線的關(guān)鍵振蕩模態(tài)的阻尼參數(shù)，即得到系統(tǒng)的最小阻尼比，從而表征系統(tǒng)的穩(wěn)定裕度。

5 結(jié)論

本文基于時(shí)域仿真方法獲取系統(tǒng)在初始擾動(dòng)下的狀態(tài)變量的時(shí)域曲線，通過(guò)曲線的收斂程度，判斷系統(tǒng)的小干擾穩(wěn)定性，并通過(guò)HHT方法分析曲線獲得系統(tǒng)的穩(wěn)定裕度。該方法的優(yōu)點(diǎn)如下：

(1)無(wú)需進(jìn)行特征值計(jì)算，通過(guò)波形直觀判斷系統(tǒng)的小干擾穩(wěn)定性；通過(guò)對(duì)波形分析可獲得系統(tǒng)的穩(wěn)定裕度。

(2)對(duì)模型及規(guī)模沒(méi)有要求，適合任意模型下的大小系統(tǒng)的小干擾穩(wěn)定計(jì)算。

(3)線性化的計(jì)算能夠保證計(jì)算量和計(jì)算時(shí)間達(dá)到在線計(jì)算的要求，為小干擾穩(wěn)定在線分析提供了新的思路。

[1] 朱方，趙紅光，劉增煌，等(Zhu Fang，Zhao Hongguang，Liu Zenghuang，et al.)．大區(qū)電網(wǎng)互聯(lián)對(duì)電力系統(tǒng)動(dòng)態(tài)穩(wěn)定性的影響(The influence of large power grid interconnected on power system dynamic stability)[J]．中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2007，27(1)：1-7.

[2] 余貽鑫，李鵬(Yu Yixin，Li Peng)．大區(qū)電網(wǎng)弱互聯(lián)對(duì)互聯(lián)系統(tǒng)阻尼和動(dòng)態(tài)穩(wěn)定性的影響(The impact of weak internection of bulk power grids to damping and dynamic stability of power systems)[J]．中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2005，25(11)：6-11．

[3] 于之虹，李芳，孫璐，等(Yu Zhihong，Li Fang，Sun Lu，et al.)．小干擾穩(wěn)定調(diào)度控制策略在線計(jì)算方法(An enhanced algorithm of online small signal stability control strategy generation)[J]，中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2014，34(34)：6191-6198．

[4] 顧麗鴻，周孝信，嚴(yán)劍峰，等(Gu Lihong，Zhou Xiaoxin，Yan Jianfeng，et al.)．特高壓聯(lián)網(wǎng)區(qū)域?qū)崟r(shí)小干擾穩(wěn)定分析策略(Real time small signal stability analysis of power systems interconnected)[J]．中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2010，30(13)：1-7．

[5] 和萍，文福拴，薛禹勝，等(He Ping，Wen Fushuan，Xue Yusheng，et al.)．風(fēng)力發(fā)電對(duì)電力系統(tǒng)小干擾穩(wěn)定性影響述評(píng)(Survey on impact of wind power on small signal stability in power system)[J]．電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings of the CSU-EPSA)，2014，26(1)：1-8．

[6] 蔡暉，黃俊輝，王海潛，等(Cai Hui，Huang Junhui，Wang Haiqian，et al.)．大規(guī)模電動(dòng)汽車智能充電站對(duì)電力系統(tǒng)小干擾穩(wěn)定性的影響(Influences of large-scale smart charging stations on small-signal stability of power systems)[J]．中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2014，34(16)：2561-2573．

[7] 葛潤(rùn)東，劉文穎，郭鵬，等(Ge Rundong，Liu Wenying，Guo Peng，et al.)．電力系統(tǒng)低頻振蕩預(yù)警及動(dòng)態(tài)阻尼控制策略研究(Research of low-frequency oscillation warning and dynamic damping control strategy in power system)[J]．電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy)，2015，34(2)：7-12．

[8] G H Golub，C F Van Loan．Matrix computations [M]．Baltimore, Maryland, USA：The John Hopkins University Press，1996. 321-322．

[9] 趙文愷，嚴(yán)正，房鑫炎，等(Zhao Wenkai，Yan Zheng，F(xiàn)ang Xinyan，et al.)．電力系統(tǒng)小信號(hào)穩(wěn)定全特征值分析的改進(jìn)BR算法(Complete eigenanalysis of power system small signal stability using the improved BR algorithm)[J]．中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2014，34(10)：1599-1608．

[10] 王錫凡，方萬(wàn)良，杜正春(Wang Xifan, Fang Wanliang, Du Zhengchun)．現(xiàn)代電力系統(tǒng)分析(Modern Power System Analysis) [M]．北京：科學(xué)出版社(Beijing: Science Press)，2003．

[11] 杜正春，劉偉，方萬(wàn)良，等(Du Zhengchun，Liu Wei，F(xiàn)ang Wanliang，et al.)．基于JACOBI-DAVIDSON方法的小干擾穩(wěn)定性分析中關(guān)鍵特征值計(jì)算(The application of the JACOBI-DAVIDSON method to the calculation of critical eigenvalues in the small signal stability analysis)[J]．中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2005，25(14)：19-24．

[12] 李崇濤，崔勇，任大偉，等(Li Chongtao，Cui Yong，Ren Dawei，et al.)．使用雙邊Jacobi-Davidson計(jì)算大規(guī)模電力系統(tǒng)關(guān)鍵模態(tài)(Calculation of critical modes of large-scale power systems using two-sided Jacobi-Davidson method)[J]．電力系統(tǒng)自動(dòng)化(Automation of Electric Power Systems)，2012，36(12)：7-11．

[13] 田鵬飛，劉崇茹，贠飛龍，等(Tian Pengfei, Liu Chongru, Yun Feilong，et al.)．一種動(dòng)態(tài)擴(kuò)維的IRA 算法及其在電力系統(tǒng)關(guān)鍵特征值計(jì)算中的應(yīng)用(An improved IRA algorithm and its application in computation of critical eigenvalues in small signal stability analysis)[J]．中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2015，35(17)：4318-4325．

[14] Guowei Cai，Zhenxin Li，D Y Yang, et al．Large-scale power system small signal stability analysis based on matrix exponential [A]．2011 4th International Conference on Electric Utility Deregulation and Restructuring and Power Technologies (DRPT) [C]．2011. 1021-1025．

[15] 郭謀發(fā)，王鵬，楊耿杰(Guo Moufa，Wang Peng，Yang Gengjie)．基于HHT和SVM的諧振接地系統(tǒng)故障選線新方法(New approach to detect fault line in resonant grounded system based on HHT and SVM)[J]．電工電能新技術(shù)(Advanced Technology of Electrical Engineering and Energy)，2014，33(9)：68-80．

[16] 李天云，高磊，趙妍(Li Tianyun，Gao Lei，Zhao Yan)．基于HHT的電力系統(tǒng)低頻振蕩分析(Analysis of low frequency oscillations using HHT method)[J]，中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE)，2006，26(14)：24-30．

[17] G Rilling, P Flandrin, P Goncalves. On empirical mode decomposition and its algorithms [A]. IEEE-EURASIP workshop on nonlinear signal and image processing [C]. 2003. 3: 8-11.

[18] P W Sauer，M Pai．Power system dynamics and stability [M]．Urbana: John Wiley & Sons, 1997.

Time domain computation and waveform analysis of small signal stability of power system

YANG Yu-de1, LI Yu1,2, YUAN Hui1, WU Zhong-biao1, YANG Jian1

(1. Guangxi Key Laboratory of Power System Optimization and Energy Technology,Guangxi University, Nanning 530004, China;2. Liuzhou Power Supply Bureau, Guangxi Power Grid Co. Ltd., Liuzhou 545005, China)

The small signal stability analysis is the analysis of a linear time-invariant system that is based on the approximate linearization of the original system. According to the characteristic of time-domain solutions which illustrates the embodiment of the system oscillation mode by the damping characteristics of the simulation curve, this paper proposes a method of small signal stability analysis based on time domain simulation. The proposed method rapidly determines the small signal stability by analysis of the degree of convergence of the time domain simulation curve of the system without the need of system state matrix eigenvalue calculation. So this method can avoid all the shortcomings of the method of eigenvalue solution, which provides a new way for the online analysis of power system small signal stability. Then, the HHT method is employed to analyze the waveforms of the simulation to obtain the minimum damping ratio which can reflect the system’s stability margin. Finally, the validity of the proposed method is illustrated by simulation of a WSCC3 system with 9 buses and a calculation example of a typical operating mode on a certain time point of a large grid.

power system; oscillation mode; time domain simulation; small signal stability; damping ratio; HHT

2016-04-29

國(guó)家自然科學(xué)基金項(xiàng)目(51541707)

陽(yáng)育德(1971-)， 男， 廣西籍， 副教授， 博士， 研究方向?yàn)殡娏ο到y(tǒng)穩(wěn)定運(yùn)行與控制； 李 雨(1989-)， 女， 廣西籍， 助理工程師， 碩士， 研究方向?yàn)殡娏ο到y(tǒng)小干擾穩(wěn)定分析與控制。

TM71

A

1003-3076(2017)01-0044-08