大型風(fēng)電場(chǎng)并網(wǎng)運行系統暫態(tài)穩定性分析

隨著(zhù)風(fēng)力發(fā)電技術(shù)的迅速發(fā)展，大容量風(fēng)電場(chǎng)并網(wǎng)運行勢在必行。由于大型風(fēng)力發(fā)電機并網(wǎng)時(shí)具有和其他常規能源電廠(chǎng)不同的特點(diǎn)，可能會(huì )對電力系統穩定性 產(chǎn)生一定影響[1～8]。建立風(fēng)力發(fā)電機組及風(fēng)電場(chǎng)的數學(xué)模型并進(jìn)行仿真分析是重要的手段之一。近年來(lái)，大多側重于風(fēng)力發(fā)電系統本身的穩定性及其模型對穩定性影響的研究，而對風(fēng)電并入后電網(wǎng)側發(fā)生大擾動(dòng)的系統暫態(tài)穩定性未作深入研究。文獻[9，10]對包含風(fēng)電場(chǎng)的電力系統小干擾穩定性進(jìn)行了分析，但建模中忽略了風(fēng)力機的傳動(dòng)部分。傳動(dòng)軸連接風(fēng)輪和發(fā)電機兩個(gè)慣量相差很大的器件必然是柔性的，在風(fēng)力機模型中考慮傳動(dòng)軸的柔性可更準確地反映擾動(dòng)后的動(dòng)態(tài)過(guò)程。

本文以恒速恒頻風(fēng)力發(fā)電系統為研究對象，在風(fēng)力發(fā)電機建模中考慮傳動(dòng)軸的柔性，采用風(fēng)力機傳動(dòng)部分的兩質(zhì)量塊模型，在PSAT軟件中搭建仿真模型，分析大型風(fēng)電場(chǎng)并入電網(wǎng)側發(fā)生大擾動(dòng)時(shí)的系統暫態(tài)穩定性，并通過(guò)加入TCSC改善了包含風(fēng)電場(chǎng)的系統暫態(tài)穩定性。

1 風(fēng)力發(fā)電機的數學(xué)模型

風(fēng)電場(chǎng)由多臺風(fēng)力發(fā)電機按一定規則排列構成。風(fēng)力發(fā)電機主要由風(fēng)力機和發(fā)電機等主要元件組成，自然界的風(fēng)推動(dòng)風(fēng)輪葉片，將風(fēng)能轉化為機械能。風(fēng)力機的機械傳動(dòng)再將機械能傳遞到發(fā)電機的轉子上。恒速恒頻風(fēng)力發(fā)電機由風(fēng)輪、低速軸、增速齒輪箱、高速軸、發(fā)電機和無(wú)功補償電容器組組成。

1.1風(fēng)力機的數學(xué)模型

風(fēng)力機機械傳動(dòng)部分如圖1所示。

齒輪箱兩側的傳動(dòng)軸以?xún)煞N轉速旋轉，連接了三個(gè)質(zhì)量塊，轉速比為1：n。因可忽略齒輪箱的慣量，故本文將風(fēng)輪和低速軸的慣量、低速軸的剛性系數和阻尼系數轉換到高速側(分別除以n2)就可將圖1的機械傳動(dòng)部分轉化為兩質(zhì)量塊模型，如圖2所示。

兩質(zhì)量塊風(fēng)機傳動(dòng)部分用s函數表示為：

θk=S(ωt-ωg)

Tw-Tk=(JtS+Dt)ωt

Tk-Tem=(JgenS+Dg)ωg

Tk=(k/S+Dtg)(ωt-ωg)

式中，Jt為風(fēng)輪轉換到高速側的慣量;Jgen為發(fā)電機轉子慣量;θk為風(fēng)輪和發(fā)電機轉子間的轉角差;ωt為風(fēng)輪的角速度;ωg為發(fā)電機的角速度;k為 傳動(dòng)軸剛性系數;Dtg為傳動(dòng)軸阻尼系數;Dt為風(fēng)輪的阻尼系數;Dg為發(fā)電機的阻尼系數;Tw為風(fēng)推動(dòng)風(fēng)機產(chǎn)生的機械轉矩，轉換到高速側應除以n。

其中

式中，ρ為空氣密度，kg/m3;ωw為風(fēng)機的機械轉速，r/s;R為風(fēng)機的葉輪半徑，m;Cp(λ)為風(fēng)力機的風(fēng)能利用系數，即單位時(shí)間內風(fēng)力機所吸收的風(fēng)能與通過(guò)葉片旋轉面的全部風(fēng)能之比。按貝茲理論[11]，Cp，max(λ)取0.593，與風(fēng)力機的葉尖速比λ(風(fēng)力機葉片頂端線(xiàn)速度與風(fēng)速之比)有關(guān)，λ=ωwR/V;Vin、Vout分別為風(fēng)力機的切入風(fēng)速和切出風(fēng)速。風(fēng)速變化的時(shí)空模型由基本風(fēng)Vwa、陣風(fēng)Vwt、漸變風(fēng)Vwr、噪聲風(fēng)Vwg四部分組成[10]。

1.2 發(fā)電機的數學(xué)模型

異步發(fā)電機采用考慮轉子機電暫態(tài)模型：

發(fā)電機機電轉矩為：

Tem= E′qiq + E′did

定子電壓為：

U= - ( r1+ jx′ ) I + E′

式中，s為異步發(fā)電機滑差，s=(ωr-ω0)/ω0;ωr為發(fā)電機轉子電角速度，ω0為同步轉速;x′為發(fā)電機暫態(tài)電抗，x′=x1+x2xm/(x2+xm) ，x=x1+xm，x為發(fā)電機電抗，x1為定子漏抗，xm為激磁電抗，x2為轉子漏抗;r1、r2分別為定子電阻、轉子電阻;T′d0為定子開(kāi)路時(shí)轉子回路的時(shí)間常數，T′d0=(x2+xm)/(2πf0r);E′= E′d+jE′q，E′為暫態(tài)電勢，E′d為d軸暫態(tài)電勢，E′q為q軸暫態(tài)電勢。

2、PSAT軟件暫態(tài)仿真

PSAT為基于MATLAB環(huán)境的電力系統分析工具箱，可進(jìn)行潮流、連續潮流、最優(yōu)潮流的計算及分岔分析、小干擾穩定分析、時(shí)域仿真及PMU配置，并支持用戶(hù)自定義模型。模型數據的輸入有兩種方式：①*.m文件格式輸 入;②PSAT支持在SIMULINK環(huán)境下從模型庫中拖拽元件搭接仿真算例(*.mdl文件)。SIMULINK環(huán)境下編譯的文件最終是解釋①的文本文 件，由MATLAB函數完成計算。因此，搭建暫態(tài)仿真模型時(shí)要符合元件使用的可行性規則[12]，確保*.mdl文件編譯成功并生成相應的*.m文件。

3、算例分析

在PSAT中搭建IEEE39 節點(diǎn)仿真模型。IEEE39節點(diǎn)系統接線(xiàn)如圖3所示，具體參數見(jiàn)文獻[13]。發(fā)電機總有功/無(wú)功為6192.9MW/923.72Mvar，負荷總有功 /無(wú)功為6098.1MW/1408.9Mvar。母線(xiàn)31取為松弛節點(diǎn)。同步發(fā)電機采用四階模型。每臺風(fēng)力發(fā)電機出力2MW，Jt=5.0，Jgen=1.0，k=0.3，Dtg=0.02，Dt=0.02，Dg=0.02，R=37.5，n=89，r1=0.01，x1=0.1，r2=0.01，x2=0.08，xm=3.0，機端電容補償無(wú)功使其功率因數達到0.98。

系統在1s時(shí)線(xiàn)路2—3靠近3處發(fā)生三相接地短路，故障在6個(gè)周波(1.1s)時(shí)由斷開(kāi)線(xiàn)路2—3而被消除。線(xiàn)路2—3在4s時(shí)重合閘。以下分系統無(wú)風(fēng)電場(chǎng)、接入風(fēng)電場(chǎng)、有風(fēng)電場(chǎng)并接入TCSC三種情況，比較系統發(fā)生上述故障并斷開(kāi)線(xiàn)路、重合閘動(dòng)作時(shí)的暫態(tài)穩定性。

3.1 原IEEE39節點(diǎn)系統(無(wú)風(fēng)電場(chǎng))

圖4虛線(xiàn)為無(wú)風(fēng)電場(chǎng)的原系統發(fā)生故障后，位于母線(xiàn)30、39的#10、#9兩臺同步發(fā)電機與松弛節點(diǎn)31的相對搖擺角隨時(shí)間變化曲線(xiàn)。圖5虛線(xiàn)為#9、#10發(fā)電機角速度與時(shí)間的關(guān)系曲線(xiàn)。

3.2 434MW風(fēng)電場(chǎng)

(1)接入母線(xiàn)8，同時(shí)原IEEE系統中10臺同步發(fā)電機出力各減小7%。圖4實(shí)線(xiàn)為接有風(fēng)電場(chǎng)的系統 發(fā)生故障后，#9、#10同步發(fā)電機與松弛節點(diǎn)31的相對搖擺角隨時(shí)間變化曲線(xiàn)。圖5實(shí)線(xiàn)為#9、#10發(fā)電機角速度與時(shí)間的關(guān)系曲線(xiàn)?？梢?jiàn)，較之原 IEEE系統，總出力7%的風(fēng)電場(chǎng)接入電網(wǎng)后系統發(fā)生故障時(shí)發(fā)電機的最大相對搖擺角增大，衰減較緩，系統的暫態(tài)穩定性降低，出現差異主要原因尚待進(jìn)一步探討。

(2)系統母線(xiàn)8接入7%出力的風(fēng)電場(chǎng)并在線(xiàn)路5—6安裝TCSC。系統發(fā)生上述故障。TCSC模型見(jiàn)圖6，補償度為30%，Tr=0.5，αmax=0.5，αmin=-0.5，KP=5，KI=1，XL=0.2，XC=0.1，Kr=10。

圖7為#9、#10發(fā)電機的相對搖擺角隨時(shí)間變化的曲線(xiàn)，其中虛線(xiàn)為未裝TCSC含風(fēng)電場(chǎng)的電力系統，實(shí)線(xiàn)為安裝TCSC含風(fēng)電場(chǎng)的電力系統?？梢?jiàn)，安裝TCSC后系統的暫態(tài)穩定性好轉，振蕩持續時(shí)間縮短。TCSC在暫態(tài)過(guò)程中，快速地改變了線(xiàn)路5—6的電氣距離，對包含風(fēng)電場(chǎng)的電力系統暫態(tài)穩定性有所 提高。

4、結語(yǔ)

a 建立了考慮風(fēng)力機傳動(dòng)部分兩質(zhì)量塊的風(fēng)力發(fā)電機數學(xué)模型，搭建PSAT仿真模型，對比分析了電網(wǎng)側發(fā)生故障情況下大型風(fēng)電場(chǎng)接入前后的系統暫態(tài)穩定性。

b 通過(guò)仿真分析，得出風(fēng)電場(chǎng)并入電力系統后系統的暫態(tài)穩定性降低，系統中同步發(fā)電機的最大相對搖擺角增大，衰減減緩;加入TCSC后，包含風(fēng)電場(chǎng)的電力系統暫態(tài)穩定性有所提高。

參考文獻：

[1]王承煦，張源. 風(fēng)力發(fā)電[M].北京：中國電力出版社，2003.

[2]Slootweg J G.Wind Power：Modeling and Impacton Power System Dynamic [D]. Delft University of Technology ，2003.

[3]李東東.風(fēng)力發(fā)電機組并網(wǎng)控制與仿真分析[J]. 水電能源科學(xué)，2006，24(1)：9211

[4]雷亞洲.與風(fēng)電并網(wǎng)相關(guān)的研究課題[J].電力系統自動(dòng)化，2003，27(8)：84289

[5]馬幼捷，張繼東.風(fēng)電場(chǎng)的穩定問(wèn)題[J].可再生能源，2006(3)：37239

[6]李鋒，魯一川.大規模風(fēng)力發(fā)電對電力系統的影響[J].中國電力，2006，39(11)：80284

[7]遲永寧，王偉勝，戴慧珠.大型風(fēng)電場(chǎng)對電力系統暫態(tài)穩定性的影響[J].電力系統自動(dòng)化，2006，30(5)：10214

[8]Ledesma P，Usaola J，Rodriguez J L.Transient Stability of a Fixed Speed Wind Farm [J]. Renewable Energy，2003，28(9)：38241

[9]湯宏，吳俊玲，周雙喜.包含風(fēng)電場(chǎng)電力系統的小干擾穩定分析建模和仿真[J].電網(wǎng)技術(shù)，2004，28(1)：38241

[10]吳學(xué)光，張學(xué)成.異步風(fēng)力發(fā)電系統動(dòng)態(tài)穩定性分析的數學(xué)模型及其應用[J].電網(wǎng)技術(shù)，1998，22(6)：69270

[11]Wiik J，Gjerde J O，Gj Engedal T，et al.Steady State Power System Issues Wind Farms[J] .IEEE Power Engineering Society Winter Meeting，2002(17)：3662371

[12]Power System Analysis Toolbox-Documentation for PSAT-version 2.0.0β2[EB/OL]. http：//www.power.uwaterloo.ca/~fmilano/downloads.htm，2007203208.

[13]Pai M A.Energy Function Analysis for Power System Stability [M]. Boston：BostonMA，1989.