郝正航， 王映康， 李健波， 樊 容， 魏 猛

(1.天津大學(xué)電氣與自動(dòng)化學(xué)院， 天津 300072； 2.貴州大學(xué)電氣工程學(xué)院， 貴陽(yáng) 550003； 3.重慶渝電信息通信有限公司， 重慶 400014； 4.重慶銅梁縣供電有限責(zé)任公司， 重慶 402500)

電力系統(tǒng)電磁轉(zhuǎn)矩系數(shù)的性質(zhì)及其計(jì)算方法①

郝正航1，2， 王映康3， 李健波4， 樊 容3， 魏 猛4

(1.天津大學(xué)電氣與自動(dòng)化學(xué)院， 天津 300072； 2.貴州大學(xué)電氣工程學(xué)院， 貴陽(yáng) 550003； 3.重慶渝電信息通信有限公司， 重慶 400014； 4.重慶銅梁縣供電有限責(zé)任公司， 重慶 402500)

為了揭示電磁轉(zhuǎn)矩系數(shù)的實(shí)質(zhì)，證明了用電磁轉(zhuǎn)矩系數(shù)表達(dá)電磁轉(zhuǎn)矩只是某種近似，還證明了將單機(jī)系統(tǒng)的電磁轉(zhuǎn)矩定義推廣于多機(jī)系統(tǒng)，其前提條件是所研究的機(jī)組相對(duì)于慣性中心COI(center of inertia)只有1個(gè)運(yùn)動(dòng)模式；基于狀態(tài)空間模型，提出了單機(jī)、2機(jī)和2機(jī)群的電磁轉(zhuǎn)矩系數(shù)的定義，給出了理論計(jì)算方法。在多機(jī)情形下，針對(duì)某一模式將系統(tǒng)視為2個(gè)機(jī)群，通過(guò)廣域測(cè)量系統(tǒng)得到局部系統(tǒng)COI的轉(zhuǎn)速和功角、各發(fā)電機(jī)電磁功率，從而不需等值運(yùn)算而轉(zhuǎn)化為等值2機(jī)，并測(cè)辨其電磁轉(zhuǎn)矩系數(shù)。算例分析表明，所提出的電磁轉(zhuǎn)矩系數(shù)的理論計(jì)算和辨識(shí)方法是有效的，二者能夠相互印證。

電力系統(tǒng)； 電磁轉(zhuǎn)矩； 在線辨識(shí)； 廣域測(cè)量系統(tǒng)； 小擾動(dòng)穩(wěn)定

電力系統(tǒng)機(jī)電振蕩分析有兩種基本方法[3]：其一是根據(jù)線性化狀態(tài)空間模型進(jìn)行特征分析，即計(jì)算系統(tǒng)矩陣的特征值和特征向量，從而得到振蕩頻率、阻尼和振型；其二是基于物理意義的電磁轉(zhuǎn)矩分析方法，將作用在轉(zhuǎn)子上的電磁轉(zhuǎn)矩分解為阻尼轉(zhuǎn)矩和同步轉(zhuǎn)矩，根據(jù)兩種轉(zhuǎn)矩的數(shù)值大小判斷同步發(fā)電機(jī)的小擾動(dòng)穩(wěn)定性。一般認(rèn)為，阻尼轉(zhuǎn)矩過(guò)小或?yàn)樨?fù)值將導(dǎo)致低頻振蕩。

電磁轉(zhuǎn)矩的概念不但在早期文獻(xiàn)[4～6]中有較多論述，而且至今仍是阻尼分析與控制的理論基礎(chǔ)[7～13]。文獻(xiàn)[4]在單機(jī)無(wú)窮大系統(tǒng)中，給出了電磁轉(zhuǎn)矩的定義并使用最小二乘法計(jì)算。文獻(xiàn)[5]基于派克方程推導(dǎo)了電磁轉(zhuǎn)矩系數(shù)的計(jì)算公式，但沒(méi)有考慮勵(lì)磁調(diào)節(jié)器。文獻(xiàn)[6]分別計(jì)算勵(lì)磁繞組、阻尼繞組和PSS提供的電磁轉(zhuǎn)矩，再用疊加原理得到總的電磁轉(zhuǎn)矩。文獻(xiàn)[7，8]提出了基于神經(jīng)網(wǎng)絡(luò)和卡爾曼濾波的辨識(shí)方法，該算法的抗噪性能好。文獻(xiàn)[9，10]將電磁轉(zhuǎn)矩的概念推廣到多機(jī)系統(tǒng)并提出計(jì)算方法和用于阻尼控制。文獻(xiàn)[11]通過(guò)在線辨識(shí)阻尼轉(zhuǎn)矩，預(yù)測(cè)振蕩穩(wěn)定極限。文獻(xiàn)[12]基于Heffron-Phillips模型討論了電磁轉(zhuǎn)矩系數(shù)，并導(dǎo)出與機(jī)電特征值的關(guān)系。

對(duì)電磁轉(zhuǎn)矩分析理論的研究，目前有3個(gè)關(guān)鍵問(wèn)題尚未得到很好的解決：

①電磁轉(zhuǎn)矩系數(shù)的定義和計(jì)算方法是基于單機(jī)無(wú)窮大系統(tǒng)提出的，可否直接推廣于多機(jī)電力系統(tǒng)仍然存在疑問(wèn)；

②電磁轉(zhuǎn)矩系數(shù)的定義是從辨識(shí)角度提出的，定義不夠嚴(yán)謹(jǐn),需要更深入的科學(xué)論證和完善；

③電磁轉(zhuǎn)矩法與特征分析法的聯(lián)系沒(méi)有建立；電磁轉(zhuǎn)矩系數(shù)的理論計(jì)算以及在多機(jī)電力系統(tǒng)中的應(yīng)用尚無(wú)深入研究。

本文將針對(duì)上述問(wèn)題進(jìn)行深入討論并提出解決途徑。

1 電力系統(tǒng)電磁轉(zhuǎn)矩的定義

1.1 單機(jī)無(wú)窮大系統(tǒng)的電磁轉(zhuǎn)矩定義

電力系統(tǒng)電磁轉(zhuǎn)矩(ΔTe)的最早定義是基于單機(jī)無(wú)窮大系統(tǒng)的，后來(lái)Demello和Concordia用這一概念分析低頻振蕩[4]。電磁轉(zhuǎn)矩的定義式為

ΔTe(t)=KdΔω(t)+KsΔδ(t)

(1)

式中：Δω和Δδ分別是發(fā)電機(jī)的轉(zhuǎn)速和功角增量；Kd和Ks分別稱為阻尼轉(zhuǎn)矩系數(shù)和同步轉(zhuǎn)矩系數(shù)，它們的數(shù)值大小和符號(hào)可以決定系統(tǒng)的小擾動(dòng)穩(wěn)定性。為了后文表達(dá)方便，本文補(bǔ)充電磁轉(zhuǎn)矩系數(shù)的規(guī)格化定義，即KdN=Kd/H，KsN=Ks/H。H是發(fā)電機(jī)慣性時(shí)間常數(shù)。

1.2 多機(jī)電力系統(tǒng)的電磁轉(zhuǎn)矩定義

電磁轉(zhuǎn)矩在單機(jī)無(wú)窮大系統(tǒng)中提出定義后，提出者并沒(méi)有進(jìn)一步推廣至多機(jī)系統(tǒng)，而后來(lái)的不少研究者試圖在多機(jī)電力系統(tǒng)中建立其定義。其形式[9,10]為

ΔTe,i(t)=Kd,iΔωi(t)+Ks,iΔδi(t)

(2)

式中：Δωi(t)和Δδi(t)是第i臺(tái)機(jī)相對(duì)于慣性中心COI(center of inertia)的轉(zhuǎn)速和功角?？梢?jiàn)，式(2)只是式(1)的簡(jiǎn)單推廣。

俄國(guó)學(xué)者還提出過(guò)另一種形式[12]，即

(3)

式(3)對(duì)電磁轉(zhuǎn)矩的表達(dá)與式(2)不同，它將電磁轉(zhuǎn)矩分解為n-1個(gè)分量，認(rèn)為第i臺(tái)機(jī)的電磁轉(zhuǎn)矩和其他n-1臺(tái)發(fā)電機(jī)相對(duì)轉(zhuǎn)速、相對(duì)功角有關(guān)。

2 電磁轉(zhuǎn)矩系數(shù)的數(shù)學(xué)本質(zhì)與物理本質(zhì)

2.1 電磁轉(zhuǎn)矩系數(shù)的存在性分析

綜觀以往文獻(xiàn)，電磁轉(zhuǎn)矩表達(dá)形式多樣化。若不對(duì)其深入分析而應(yīng)用，則可能導(dǎo)致錯(cuò)誤。為了從數(shù)學(xué)上認(rèn)識(shí)它，給出命題1。

命題1在實(shí)數(shù)范圍內(nèi)，嚴(yán)格滿足式(1)～式(3)的電磁轉(zhuǎn)矩系數(shù)Kd和Ks、Kd，i和Ks,i，Kd,ij和Ks,ij都不存在。

證明：使用反證法。首先考察單機(jī)無(wú)窮大系統(tǒng)。假設(shè)存在Kd和Ks使式(1)成立。根據(jù)電力系統(tǒng)暫態(tài)分析理論，同步電機(jī)電磁轉(zhuǎn)矩是橫、縱軸次暫態(tài)電勢(shì)和功角的多元函數(shù)，表示為

(4)

根據(jù)式(4)，小擾動(dòng)下的偏差方程為

(5)

正常運(yùn)行條件下，C1和C2不全為零。將式(5)代入式(1)可得

(6)

由于不存在滿足式(1)～式(3)的電磁轉(zhuǎn)矩系數(shù)(陣)，所以用電磁轉(zhuǎn)矩系數(shù)表達(dá)電磁轉(zhuǎn)矩只是某種近似。以上三式本質(zhì)上都是矛盾方程組，求解電磁轉(zhuǎn)矩系數(shù)就是求解矛盾方程的近似解。

2.2 電磁轉(zhuǎn)矩系數(shù)與系統(tǒng)動(dòng)態(tài)等值的關(guān)系

電磁轉(zhuǎn)矩系數(shù)除了原始的物理意義外，本文給出另一物理解釋。首先討論單機(jī)情形，假設(shè)原動(dòng)機(jī)機(jī)械轉(zhuǎn)矩Tm恒定，轉(zhuǎn)子運(yùn)動(dòng)方程可表達(dá)為

(7)

式(1)代入式(7)，并考慮規(guī)格化定義，得

(8)

由此可見(jiàn)，當(dāng)電磁轉(zhuǎn)矩用式(1)表達(dá)時(shí)，轉(zhuǎn)子運(yùn)動(dòng)方程中不再包含Δδ及Δω之外的其他狀態(tài)變量，相當(dāng)于將詳細(xì)模型X=AX等值為二階模型。所以，求解Kd和Ks，本質(zhì)上就是進(jìn)行動(dòng)態(tài)等值。

等值后的系統(tǒng)不再包含轉(zhuǎn)子各繞組磁鏈模式、勵(lì)磁機(jī)模式及調(diào)節(jié)器模式，僅剩機(jī)械模式。

依據(jù)系統(tǒng)等值觀點(diǎn)，式(8)表達(dá)的降階系統(tǒng)特征值應(yīng)該等于原系統(tǒng)的機(jī)電特征值(λ=σ±jω)。因此，問(wèn)題轉(zhuǎn)化為：當(dāng)Kd和Ks取何值時(shí)，降階系統(tǒng)(式(8))特征值等于期望(原系統(tǒng))特征值。

可見(jiàn)，這是給定特征值求微分方程系數(shù)的問(wèn)題，屬代數(shù)特征值反問(wèn)題的范疇[13]。對(duì)于單機(jī)(或二機(jī))系統(tǒng)，該反問(wèn)題的解存在且唯一。對(duì)于式(8)，反問(wèn)題的解為

(9)

對(duì)于多機(jī)情形(式(2))，同樣表達(dá)的是系統(tǒng)等值關(guān)系。仿照前面不難分析，式(2)將多機(jī)系統(tǒng)等值為

(10)

式(10)表明，第i臺(tái)機(jī)相對(duì)于COI的運(yùn)動(dòng)是獨(dú)立的，以二階微分方程表達(dá)，只能體現(xiàn)一個(gè)振蕩模式的動(dòng)態(tài)行為?；蛘哒f(shuō)，第i臺(tái)機(jī)在COI坐標(biāo)下的轉(zhuǎn)子搖擺曲線只包含一個(gè)頻率。據(jù)此可以提出如下命題。

命題2在多機(jī)系統(tǒng)中，將電磁轉(zhuǎn)矩表示為式(2)的前提條件是，所研究的機(jī)組相對(duì)于COI只有一個(gè)主導(dǎo)模式。該條件是苛刻的。

最后考察式(3)，它表達(dá)的也是一種動(dòng)態(tài)等值方案。等值系統(tǒng)為

從系統(tǒng)等值觀點(diǎn)，式(11)的n-1對(duì)特征值應(yīng)和原系統(tǒng)n-1對(duì)特征值相等，也屬于一類矩陣特征值的反問(wèn)題。然而遺憾的是，該類特征值反問(wèn)題的解是不確定的，即可能無(wú)解，也可能多解[13]。即，無(wú)法保證電磁轉(zhuǎn)矩系數(shù)陣的存在性和唯一性。故式(3)定義的科學(xué)性還需商榷。

(11)

綜上所述，電磁轉(zhuǎn)矩系數(shù)的數(shù)學(xué)本質(zhì)是線性矛盾方程組的近似解，或者是特征值反問(wèn)題的解；物理本質(zhì)是僅保留機(jī)電模式條件下，對(duì)原系統(tǒng)的降階等值表達(dá)。

3 電磁轉(zhuǎn)矩系數(shù)的理論計(jì)算

3.1 基于狀態(tài)空間模型的電磁轉(zhuǎn)矩系數(shù)的定義

應(yīng)當(dāng)說(shuō)明，式(1)和式(2)的電磁轉(zhuǎn)矩定義是從辨識(shí)角度提出的，求解Kd、Ks(Kd,i、Ks,i)只能以“曲線擬合”的手段得到最小二乘解[4，8](即辨識(shí)值)。為了評(píng)價(jià)辨識(shí)值的有效性，也為了溝通電磁轉(zhuǎn)矩法與狀態(tài)空間法之間的聯(lián)系。下面提出一種定義及計(jì)算方法。

基于狀態(tài)空間模型，從物理意義出發(fā)，給出2機(jī)(2區(qū))系統(tǒng)電磁轉(zhuǎn)矩系數(shù)的定義。

定義1在二機(jī)系統(tǒng)中，規(guī)格化的電磁轉(zhuǎn)矩系數(shù)定義如下：

KdN=Re[G(s)]s=jωd，KsN=-Im[G(s)]s=jωd，

定義2在A、B二區(qū)域多機(jī)系統(tǒng)中，規(guī)格化的電磁轉(zhuǎn)矩系數(shù)定義如下：

KdN=Re[G(s)]s=jωd，KsN=-Im[G(s)]s=jωd，其中

對(duì)定義2的說(shuō)明：定義2是對(duì)定義1的推廣，用子系統(tǒng)(局部)慣性中心的轉(zhuǎn)速，代替定義1中單臺(tái)電機(jī)轉(zhuǎn)速；用子系統(tǒng)各電機(jī)功率和，代替單臺(tái)電機(jī)功率。所以，定義2的物理意義仍然是明確的。

3.2 電磁轉(zhuǎn)矩定義在多機(jī)系統(tǒng)的推廣問(wèn)題

4 電磁轉(zhuǎn)矩系數(shù)的在線辨識(shí)

由第3節(jié)的定義可見(jiàn)，電磁轉(zhuǎn)矩系數(shù)的理論計(jì)算運(yùn)算量很大，超過(guò)了特征分析法的計(jì)算量，但作為理論研究的基本環(huán)節(jié)是必要的。電磁轉(zhuǎn)矩系數(shù)的在線辨識(shí)往往更具實(shí)用性。單機(jī)情形已有較多研究，以下著重討論2機(jī)及多機(jī)系統(tǒng)的情況。

4.1 二機(jī)系統(tǒng)電磁轉(zhuǎn)矩系數(shù)在線辨識(shí)

仿照單機(jī)系統(tǒng)Kd、Ks的辨識(shí)方法[4]，提出2機(jī)系統(tǒng)的辨識(shí)方案：設(shè)2機(jī)系統(tǒng)的ΔPe,1，ΔPe,2、Δω12和Δδ12均通過(guò)廣域測(cè)量系統(tǒng)(WAMS)獲得N個(gè)采集數(shù)據(jù)。假設(shè)

ΔPe,1(i)/H1-ΔPe,21(i)/H2=

KsNΔδ12(i)+KdNΔω12(i)

i=1，2，…，N

(12)

令：ΔP=ΔPe,1/H1-ΔPe,2/H2

B=[ΔP(1)，ΔP(2)，…，ΔP(N)]T；

A1=[Δδ12(1)，Δδ12(2)，…，Δδ12(N)]T；

A2=[Δω12(1)，Δω12(2)，…，Δω12(N)]T；

A=[A1，A2]；K=[KsN，KdN]T

則式(12)的矩陣表達(dá)式為

AK=B

(13)

由命題1知，式(13)是關(guān)于K的矛盾方程組。根據(jù)數(shù)值代數(shù)理論，該矛盾方程的法方程為ATAK=ATB，法方程的解K=(ATA)-1ATB是式(13)的最小二乘解，由此得出KdN、KsN的估計(jì)值。

4.2 多機(jī)系統(tǒng)電磁轉(zhuǎn)矩系數(shù)在線辨識(shí)

5 算例分析

5.1 2種理論計(jì)算結(jié)果比較示例

為了檢驗(yàn)本文提出的基于狀態(tài)空間模型的電磁轉(zhuǎn)矩定義與文獻(xiàn)[5]定義的等價(jià)性，以單機(jī)無(wú)窮大系統(tǒng)為例(系統(tǒng)參數(shù)見(jiàn)文獻(xiàn)[17])，采用2種定義分別計(jì)算Kd和Ks。計(jì)算結(jié)果如表1所示。

從表1看出，在各種負(fù)荷水平下，2種定義的相對(duì)誤差均小于3×10-3?？梢?jiàn)，2種定義是等價(jià)的，微小誤差是由計(jì)算工具有限精度所致。定義1是建立在狀態(tài)空間模型上，適用于任何復(fù)雜程度的模型。而文獻(xiàn)[5]定義不能考慮調(diào)節(jié)器作用。

表1 單機(jī)無(wú)窮大系統(tǒng)電磁轉(zhuǎn)矩系數(shù)2種理論計(jì)算結(jié)果的比較(p.u.)

5.2 多機(jī)系統(tǒng)功角軌跡功率譜分析

4機(jī)2區(qū)域電力系統(tǒng)如圖1所示，系統(tǒng)參數(shù)見(jiàn)文獻(xiàn)[6]。為了說(shuō)明多機(jī)系統(tǒng)中電磁轉(zhuǎn)矩法的局限性，首先考察功角軌跡功率譜。圖2是受擾后發(fā)電機(jī)功角搖擺軌跡的功率譜。

圖1 2區(qū)域4機(jī)示例系統(tǒng)

功率譜波峰分布直觀地顯示了時(shí)域信號(hào)所包含的頻率成分。圖2(a)是發(fā)電機(jī)G2在COI坐標(biāo)下的功角搖擺曲線的功率譜。可以看出，它包含3.422 rad/s和6.076 rad/s兩個(gè)頻率成分。如4節(jié)所述，這正是不能用式(2)表達(dá)其電磁轉(zhuǎn)矩的原因。圖2(b)是區(qū)域1和區(qū)域2局部COI的相對(duì)搖擺軌跡的功率譜，該譜只有一個(gè)顯著模式，頻率為3.422 rad/s，即區(qū)間搖擺模式，故等值為2機(jī)。

5.3 2區(qū)域系統(tǒng)電磁轉(zhuǎn)矩計(jì)算與辨識(shí)

為檢驗(yàn)定義2的合理性，按此定義計(jì)算圖1示例系統(tǒng)的電磁轉(zhuǎn)矩系數(shù)。先采用格式化的方法形成系統(tǒng)矩陣A，然后依定義2求出傳遞函數(shù)并得到不同負(fù)荷水平下的KdN和KsN，即理論計(jì)算結(jié)果。同時(shí)驗(yàn)證第4.2節(jié)的在線辨識(shí)方案，基于電力系統(tǒng)分析綜合程序(PSASP)，建立仿真系統(tǒng)并設(shè)計(jì)擾動(dòng)試驗(yàn)，采集受擾系統(tǒng)軌跡數(shù)據(jù)，按第4.2節(jié)方法和步驟，得KdN和KsN，即辨別值。 基于定義2和第4.2節(jié)辨識(shí)方法的兩種計(jì)算結(jié)果均列入表2。

圖2 發(fā)電機(jī)功角搖擺軌跡的功率譜分析

區(qū)域1有功出力P機(jī)電特征值定義2KdN測(cè)辨值KdN定義2KsN測(cè)辨值KsN12．2-0．1336±j3．42200．28310．270811．93011．46212．4-0．1275±j3．41100．27220．259411．86411．39912．6-0．1209±j3．40000．26050．248111．79511．32212．8-0．1140±j3．38910．24810．234711．72311．24413．0-0．1066±j3．37700．23480．222111．64711．160

由表2可見(jiàn)，定義2的理論值和辨識(shí)值是接近的，同時(shí)還看到，電磁轉(zhuǎn)矩系數(shù)和特征值的關(guān)系與式(9)相符，達(dá)到三者相互統(tǒng)一和印證。

6 結(jié)語(yǔ)

本文證明了電力系統(tǒng)的電磁轉(zhuǎn)矩表達(dá)式不能恒等，只是近似相等。這種近似關(guān)系和動(dòng)態(tài)等值存在內(nèi)在聯(lián)系，求解電磁轉(zhuǎn)矩系數(shù)過(guò)程就是動(dòng)態(tài)等值的過(guò)程，從而揭示了電磁轉(zhuǎn)矩系數(shù)的數(shù)學(xué)本質(zhì)和物理本質(zhì)。在單機(jī)或2機(jī)系統(tǒng)中，等值系統(tǒng)保留了原系統(tǒng)完整的機(jī)電動(dòng)態(tài)，但在多機(jī)系統(tǒng)中，等值系統(tǒng)只能保留一個(gè)主要模式，由此在多機(jī)系統(tǒng)中推廣電磁轉(zhuǎn)矩的定義需要苛刻的假設(shè)條件。為了使電磁轉(zhuǎn)矩分析法在理論上更加完整，并為在線辨識(shí)提供參考標(biāo)準(zhǔn)，本文還基于狀態(tài)空間模型，提出了電磁轉(zhuǎn)矩系數(shù)理論計(jì)算公式。在多機(jī)系統(tǒng)中，電磁轉(zhuǎn)矩系數(shù)的在線辨識(shí)應(yīng)在分群基礎(chǔ)上進(jìn)行，針對(duì)某一模式劃分為2機(jī)群，按等值2機(jī)進(jìn)行電磁轉(zhuǎn)矩系數(shù)辨識(shí)是可行的。

[1] 徐政. 交直流電力系統(tǒng)動(dòng)態(tài)特性行為分析[M]. 北京: 機(jī)械工業(yè)出版社, 2005.

[2] 周雙喜, 蘇小林(Zhou Shuangxi, Su Xiaolin). 電力系統(tǒng)小干擾穩(wěn)定性研究的新進(jìn)展(Advance of study on small disturbance stability of power systems)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings of the CSU-EPSA), 2007, 19(2) : 1- 8, 34.

[3] 姜齊榮, 陳湘, 王天,等(Jiang Qirong, Chen Xiang, Wang Tian,etal). 多機(jī)系統(tǒng)復(fù)轉(zhuǎn)矩系數(shù)分析及PSS參數(shù)計(jì)算(Analysis of complex torque coefficient in multi-machines system and PSS parameters assignment)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings of the CSU-EPSA), 2010, 22(1): 20-25.

[4] DeMello F P,Concordia C. Concept of synchronous machine stability as affected by excitation control[J].IEEE Trans on Power Systems, 1969, 88(4): 316-329.

[5] 高景德, 王祥珩, 李發(fā)海. 交流電機(jī)及其系統(tǒng)的分析[M].北京: 清華大學(xué)出版社,1992.

[6] Kundur P. Power System Stability and Control[M].New York: McGraw-Hill, 1994.

[7] Abu-Al-Feilat E, Bettayeb M, Al-Duwaish H,etal.Neural network-based approach for on-line dynamic stability assessment using synchronizing and damping torque coefficients[J].Electric Power Systems Research, 1996, 39(2): 103-110.

[8] Abu-Al-Feilat E A, Younan N, Grzybowski S. Estimating the synchronizing and damping torque coefficients using Kalman filtering[J].Electric Power Systems Research, 1999, 52(2): 145-149.

[9] Shaltout A A, Abu Al-Feilat K A. Damping and synchronizing torque computation in multi-machine power systems[J].IEEE Trans on Power Systems, 1992, 7(1): 280-286.

[10]Swift F J, Wang H F. The connection between modal analysis and electric torque analysis in studying the oscillation stability of multi-machine power systems[J].International Journal of Electrical Power and Energy System, 1997, 19(5): 321-330.

[11]Ghasemi Hassan, Canizares Claudio. On-line damping torque estimating and oscillatory stability margin prediction[J].IEEE Trans on Power Systems, 2007, 22(2): 667-674.

[12]王青,閔勇,張毅威(Wang Qing, Min Yong, Zhang Yiwei).多機(jī)電力系統(tǒng)電磁轉(zhuǎn)矩分析方法(Electric torque analysis in multi-machine power systems)[J].清華大學(xué)學(xué)報(bào)：自然科學(xué)版(Journal of Tsinghua University：Sci & Tech), 2008,48(1):9-12.

[13]周樹(shù)荃,戴華. 代數(shù)特征值反問(wèn)題[M]. 鄭州: 河南科學(xué)技術(shù)出版社, 1990.

[14]余貽鑫. 線性系統(tǒng)[M]. 北京: 科學(xué)出版社, 2009.

[15]李鵬,余貽鑫,孫強(qiáng),等(Li Peng, Yu Yixin,Sun Qiang,etal).基于Prony分析的多機(jī)系統(tǒng)電磁轉(zhuǎn)矩系數(shù)計(jì)算(Computation of electric torque in multi-machine power systems based on Prony analysis)[J].電網(wǎng)技術(shù)(Power System Technology),2006,30(10):39-44.

[16]趙書強(qiáng), 常鮮戎, 潘云江,等(Zhao Shuqiang, Chang Xianrong, Pan Yunjiang,etal). 基于部分慣量中心等值的多機(jī)系統(tǒng)特征根計(jì)算(A partial center of inertia transformation based method for calculating eigenvalues of multi-machine power system)[J].電力系統(tǒng)自動(dòng)化(Automation of Electric Power Systems), 1998,22(9):43-45,76.

[17]趙書強(qiáng),常鮮戎,賀仁睦,等(Zhao Shuqiang, Chang Xianrong, He Renmu,etal). PSS控制過(guò)程中的借阻尼現(xiàn)象與負(fù)阻尼效應(yīng)(Borrow damping phenomena and negative damping effect of PSS control)[J].中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings of the CSEE),2004,24(5):7-11.

CharacterofElectricTorqueCoefficientsinPowerSystemandItsComputation

HAO Zheng-hang1，2， WANG Ying-kang3， LI Jian-bo4， FAN Rong3， WEI Meng4

(1. School of Electrical Engineering and Automation, Tianjin University, Tianjin 300072, China; 2.School of Electrical Engineering, Guizhou University, Guiyang 550003, China; 3.Chongqing Yudian Information & Communication Company Limited, Chongqing 400014, China; 4.Chongqing Tongliang Power Supply Company Limited, Chongqing 402500, China)

In order to investigate the essence of electrical torque coefficients, it is proved that the equation between active power and electric torque based on two coefficients is not satisfied completely and electric torque method lead to a simplified path of accurate power system model essentially. If the electric torque definition of single machine is applied to multi-machine system, the given machine must be dominated by only one oscillation mode on the center of inertia(COI). This paper also presents electric toque definition and calculating method which can be used for single machine system, two-machine system, and two-area system. In large power system, a two-area system is defined in term of certain mode and its electric torque coefficient can be computed with measured signals (e.g.speed, angle based on partial COI and active power) from wide area measurement system(WAMS). A four-machine system case showed that the methodology proposed can be utilized to identify the coefficient of electric torques.

power system； electric torque； on-line identification； wide area measurement system(WAMS)； small-signal stabilit

2010-03-05

2010-11-16

TM743

A

1003-8930(2011)01-0086-06

郝正航(1972-)，男，博士研究生，副教授，研究方向?yàn)殡娏ο到y(tǒng)穩(wěn)定和風(fēng)力發(fā)電。Email:haozhenghang@163.com 王映康(1965-)，男，生產(chǎn)技術(shù)部主任，研究方向?yàn)殡娏ο到y(tǒng)自動(dòng)化。Email:wyk@cq.sgcc.com.cn 李健波(1970-)，男，副總經(jīng)理，研究方向?yàn)殡娏ο到y(tǒng)自動(dòng)化。Email:li.jb@cq.sgcc.com.cn