劉書成

(東北電力大學(xué)電氣工程學(xué)院，吉林吉林132012)

0 引言

隨著能源危機的不斷加深，新能源的利用和開發(fā)也越來越受到人們的重視和青睞，風(fēng)能作為一種清潔能源被納入中國的能源長期規(guī)劃范圍之內(nèi)［1］。由于風(fēng)能具有一定的隨機性，并入電力系統(tǒng)之后會給系統(tǒng)帶來一定的沖擊。因此，研究合理的風(fēng)機接入容量的問題成為一個亟待解決的現(xiàn)實問題。而其中的重要一方面就是需分析風(fēng)機接入后系統(tǒng)的穩(wěn)定性的問題。對此，許多學(xué)者曾做出一些深入的研究［2－4］。如文獻(xiàn)［5］分析了在不同運行狀態(tài)下，風(fēng)電并網(wǎng)前后互聯(lián)系統(tǒng)的小干擾穩(wěn)定性變化，對實際的穩(wěn)定性分析提供了理論依據(jù);文獻(xiàn)［6］推導(dǎo)了雙饋機組小干擾穩(wěn)定的數(shù)學(xué)模型的基礎(chǔ)上，以通遼外送型電網(wǎng)接入大規(guī)模風(fēng)電為例，系統(tǒng)研究大規(guī)模風(fēng)電外送對電力系統(tǒng)小干擾穩(wěn)定的影響;文獻(xiàn)［7］利用DIgSILENT/PowerFactory進(jìn)行頻域分析和時域仿真，分析異步風(fēng)電機組對電力系統(tǒng)小干擾穩(wěn)定性及阻尼特性的影響;文獻(xiàn)［8］考慮同步發(fā)電機勵磁調(diào)節(jié)系統(tǒng)動態(tài)行為的全系統(tǒng)狀態(tài)矩陣公式;利用特征值分析法，分析了系統(tǒng)的小干擾穩(wěn)定性?；诖?，本文先介紹小干擾穩(wěn)定分析的線性化理論，給出小干擾穩(wěn)定域的條件，然后推導(dǎo)風(fēng)機與同步機并聯(lián)的系統(tǒng)的小干擾分析的雅克比矩陣，通過改變風(fēng)機的輸入功率和原動機的機械功率，得到對應(yīng)的小干擾穩(wěn)定域。

1 小擾動穩(wěn)定域的邊界

對同時含有微分方程和代數(shù)方程的非線性系統(tǒng)，數(shù)學(xué)表達(dá)式為

式中:x為系統(tǒng)狀態(tài)變量;y為代數(shù)變量;p為控制變量。對系統(tǒng)式(1)在平衡點處線性化，可得:

式中:fx、fy、gx、gy分別為函數(shù)對狀態(tài)變量和代數(shù)變量的偏導(dǎo)數(shù)。

若矩陣gy非奇異，則式(2)可寫為

式中:A(p)=fx－fygy－1gx。

當(dāng)gy非奇異且A(p)的特征根具有負(fù)實部時，得出的小擾動穩(wěn)定域是穩(wěn)定的。

根據(jù)式(3)的方程，選擇適當(dāng)?shù)乃惴ǖ玫椒€(wěn)定域的步驟如下:

1)選定需要變化的參數(shù)值，將其它參數(shù)設(shè)定為常數(shù)。

2)確定參數(shù)的初始點，一般選擇系統(tǒng)穩(wěn)態(tài)的參數(shù)值作為初始值。

3)在對應(yīng)的參數(shù)空間中，從初始點沿某個特定的方向，并以適當(dāng)?shù)牟介L改變系統(tǒng)的參數(shù)值，得到給1個點對應(yīng)的雅可比矩陣特征值。

4)當(dāng)參數(shù)出現(xiàn)1對共軛純虛特征值且其余特征值均有負(fù)實部時，系統(tǒng)發(fā)生HB，記錄此時的參數(shù)，并停止該方向的搜索。

5)改變3)步驟中搜索邊方向，重復(fù)3)和4)步驟，得到新的邊界點。該方法在參數(shù)空間(p1，p2)中占用較多的計算時間，但具有較好的準(zhǔn)確性，可根據(jù)上個步驟的特征值自適應(yīng)地調(diào)整步長，進(jìn)而增加該方法的靈活性。

在具有多個發(fā)電機的電力系統(tǒng)中，若考慮發(fā)電機注入功率為分岔參數(shù)時，對應(yīng)的小干擾穩(wěn)定域比發(fā)電機臺數(shù)少一維的超平面;當(dāng)SSSR的維數(shù)大于三維時，已經(jīng)失去了直觀形象的優(yōu)點。常見的方法是將其投影到低維平面中，使其具有一定的實用價值，但也喪失很大一部分的有用信息。1個有價值的方案是按照二機等值的方法，將n－1維的曲面投影到n－1個2維空間上，節(jié)省了計算量，具有直觀的物理意義［9］。

2 含有異步風(fēng)機的電力系統(tǒng)

2.1 異步風(fēng)力發(fā)電機組數(shù)學(xué)模型

2.1.1 風(fēng)力機數(shù)學(xué)模型

風(fēng)力機機械功率為

風(fēng)力機機械轉(zhuǎn)矩方程為

式中:PW為風(fēng)力機機械功率;TW為風(fēng)力機機械轉(zhuǎn)矩;r為風(fēng)力機葉輪半徑;ρ為空氣密度;β為槳距角;λ=rω/v為葉尖速比;ω為風(fēng)力機轉(zhuǎn)速;Cp為風(fēng)能轉(zhuǎn)換效率系數(shù);v為風(fēng)速。

連接葉片和齒輪箱的輪轂，有較大慣性，其兩邊的轉(zhuǎn)矩可用一階慣性環(huán)節(jié)模擬，其表達(dá)式為

式中:Tt為輸入齒輪箱的機械轉(zhuǎn)矩;τh為輪轂的慣性時間常數(shù)。

齒輪箱和聯(lián)軸器傳遞風(fēng)力機和異步發(fā)電機間的轉(zhuǎn)矩，其方程表達(dá)式為

式中:Tm為齒輪箱輸出轉(zhuǎn)矩;Tt為異步發(fā)電機輸入轉(zhuǎn)矩;ττ為齒輪箱慣性時間常數(shù)。通常認(rèn)為風(fēng)力機轉(zhuǎn)速基本保持不變，有 Tm≈Tt。

傳動部分模型方程表達(dá)式為

式中:τw為風(fēng)力機慣性時間常數(shù)。

2.1.2 槳距角控制系統(tǒng)數(shù)學(xué)模型

通過調(diào)節(jié)葉片槳距角，使風(fēng)力機跟蹤 Cp達(dá)到最大值，提高風(fēng)能利用效率。槳距角控制系統(tǒng)通常表示為

式中:β0為槳距角初始值;τβ為槳距角控制系統(tǒng)的慣性時間常數(shù)。

2.1.3 異步發(fā)電機組數(shù)學(xué)模型

簡單的風(fēng)機系統(tǒng)的方程主要包括定子電壓方程、機電暫態(tài)方程以及轉(zhuǎn)子電壓方程等。

定子電壓方程形式為

機電暫態(tài)方程形式為

而發(fā)電機的暫態(tài)電抗為

式中:Ud、Uq分別為發(fā)電機定子電壓 d、q軸分量;Id、Iq分別為定子電流 d、q 軸分量;E′d、E′q分別為暫態(tài)電勢 d、q軸分量;T′d0為定子開路時間常數(shù);r1、r2、xl、x2、xm分別為異步發(fā)電機定子電阻、轉(zhuǎn)子電阻、定子電抗、轉(zhuǎn)子電抗和激磁電抗標(biāo)么值。

轉(zhuǎn)子運動方程為

式中:τj為異步發(fā)電機的慣性時間常數(shù);s為轉(zhuǎn)率差;Te為電磁轉(zhuǎn)矩。

將式(13)線性化方程為

整理上述微分方程，并對其進(jìn)行在平衡點處進(jìn)行線性化分析，得到以 Δβ、ΔTt、Δs、ΔE′d、ΔE′q為狀態(tài)變量的異步風(fēng)電機組狀態(tài)方程為

式中:ids0、iqs0為定子電流的 d、q 軸分量;E′d0、E′q0為暫態(tài)電勢的d、q軸分量。

將式(15)簡寫為

2.2 同步發(fā)電機組數(shù)學(xué)模型

同步發(fā)電機的三階模型為

式中:δ為發(fā)電機功角;ω為發(fā)電機轉(zhuǎn)速;ω0為發(fā)電機的穩(wěn)態(tài)轉(zhuǎn)速;TJ為機械轉(zhuǎn)動慣量;D為機械阻尼系數(shù);Pm為原動機機械功率;Pe為發(fā)電機電磁功率;Eq為q軸暫態(tài)電勢;Ef為勵磁的控制輸入;Td0′為發(fā)電機定子開路時勵磁繞組的時間常數(shù);Ut為發(fā)電機機端電壓;Us為無窮大母線電壓;xd為發(fā)電機的d軸同步電抗;xd∑為計入了輸電系統(tǒng)總電抗后的d軸總同步電抗;xd′為發(fā)電機d軸暫態(tài)電抗;xd∑′為計入了輸電系統(tǒng)總電抗后的d軸暫態(tài)總電抗。

對其在平衡點處進(jìn)行線性化，可得到類似于式(15)的矩陣，將其寫為

式中:(19)式的參數(shù)的形式和含義見文獻(xiàn)［11］。

式(19)同樣可以簡寫為

ΔY˙=A2ΔY。

將式(15)和式(19)聯(lián)立，可得

3 算例分析

3.1 系統(tǒng)的主要參數(shù)

為了驗證本文算法的有效性，特引入1組風(fēng)機與單臺同步機并聯(lián)的系統(tǒng)，然后與輸電線路域無窮大系統(tǒng)相連。算例系統(tǒng)如圖1所示。其中，1組風(fēng)機可以簡化為單臺的異步發(fā)電機，具體的等值方法見文獻(xiàn)［10］。同步機的主要參數(shù)如下:

線路及變壓器的參數(shù)為

單臺異步風(fēng)力發(fā)電機的主要參數(shù)如下:

圖1 算例系統(tǒng)

分別令風(fēng)機的輸入機械功率PW和發(fā)電機的輸入機械功率Pm為自變量，根據(jù)式(3)可得到對應(yīng)的小干擾穩(wěn)定域，如圖2所示。

圖2 本文所示系統(tǒng)的小干擾穩(wěn)定域

從圖2可以看出，陰影部分對應(yīng)的為系統(tǒng)的小干擾穩(wěn)定域，異步風(fēng)力發(fā)電機維持在較高出力水平，同步發(fā)電機維持在正常出力水平下就可維持系統(tǒng)的小干擾穩(wěn)定。若不滿足上述條件，則有可能出現(xiàn)系統(tǒng)的振蕩，需要其它輔助措施才能維持系統(tǒng)的小干擾穩(wěn)定性。

3.2 勵磁放大倍數(shù)對小干擾穩(wěn)定域的影響

勵磁放大倍數(shù)作為電力系統(tǒng)的重要參數(shù)，對電力系統(tǒng)的穩(wěn)定性具有一定的影響。因此，在這里通過改變勵磁放大倍數(shù)，可以得出對應(yīng)的小干擾穩(wěn)定域。通過對比可知放大倍數(shù)對小干擾穩(wěn)定域的影響。Kv=5時的小干擾穩(wěn)定域如圖3所示，Kv=20時的小干擾穩(wěn)定域如圖4所示。

圖3 Kv=5時的小干擾穩(wěn)定域

圖4 Kv=20時的小干擾穩(wěn)定域

通過對比圖2、圖3、圖4(其中圖2為勵磁放大倍數(shù)為10的穩(wěn)定域)陰影部分，可以看出勵磁放大倍數(shù)對系統(tǒng)的小干擾穩(wěn)定域有相當(dāng)大的影響，主要是對異步風(fēng)力發(fā)電機出力影響較大，對同步機出力基本上沒有影響，說明實際系統(tǒng)中合理配置系統(tǒng)的勵磁放大倍數(shù)，對系統(tǒng)的穩(wěn)定性具有至關(guān)重要的作用。

4 結(jié)論

通過線性穩(wěn)定性理論研究了含有風(fēng)機系統(tǒng)的小干擾穩(wěn)定域問題，得到如下結(jié)論:

1)系統(tǒng)的穩(wěn)定性是由異步風(fēng)機和同步發(fā)電機相互協(xié)調(diào)決定的，兩者的出力水平需要保持在適當(dāng)水平才能保證系統(tǒng)的穩(wěn)定性。

2)勵磁放大倍數(shù)對系統(tǒng)的穩(wěn)定性具有重要的影響。但相對而言，對異步風(fēng)力發(fā)電機的出力影響很大，對同步發(fā)電機沒有很大變化。

［1］雷亞洲．與風(fēng)電并網(wǎng)相關(guān)的研究課題［J］．電力系統(tǒng)自動化，2003，27(8):84 －89．

［2］張衛(wèi)紅，于吉慶，王坤．新能源發(fā)電并網(wǎng)標(biāo)準(zhǔn)比較［J］．黑龍江電力，2011，33(6):403－408．

［3］申洪．變速恒頻風(fēng)電機組并網(wǎng)運行模型及其應(yīng)用［D］．北京:中國電力科學(xué)研究院，2003．

［4］吳義純，丁明，張輝．含風(fēng)電場的電力系統(tǒng)潮流計算［J］．中國電機工程學(xué)報，2005，25(4):36－39．

［5］符楊，霍淑珍，邊曉燕，等．基于NEVA的風(fēng)電并網(wǎng)小干擾穩(wěn)定研究［J］．華東電力，2011，39(3):443－445．

［6］楊濤，鄭濤，遲永寧，等．大規(guī)模風(fēng)電外送對電力系統(tǒng)小干擾穩(wěn)定性影響［J］．中國電力 ．2010，43(6):20－25．

［7］關(guān)宏亮，遲永寧，戴慧珠，等．異步風(fēng)電機組接入系統(tǒng)的小干擾穩(wěn)定及控制［J］．電力系統(tǒng)自動化，2008，32(4):54 －58．

［8］范偉，趙書強．考慮風(fēng)力發(fā)電的電力系統(tǒng)小干擾穩(wěn)定性分析［J］．華北電力大學(xué)學(xué)報，2009，36(2):23 －27，32．

［9］郝正航，陳卓．注入功率空間上電力系統(tǒng)小擾動穩(wěn)定域的實用邊界［J］．電工電能新技術(shù)，2011，30(3):38－42．

［10］李桂紅．含風(fēng)電場電力系統(tǒng)小干擾穩(wěn)定分析［D］．太原:太原理工大學(xué)，2010:34－37．

［11］李光琦．電力系統(tǒng)暫態(tài)分析［M］．北京:中國電力出版社，2007:184－187．