舒 進1, 都勁松1, 王靖程1, 張保會2

?

基于奇異攝動降階的風(fēng)電接入系統(tǒng)阻尼分析

舒 進, 都勁松, 王靖程, 張保會

（1.西安熱工研究院有限公司，陜西 西安710043；2.西安交通大學(xué)電力設(shè)備電氣絕緣國家重點實驗室，陜西 西安 710049）

為分析風(fēng)電機組接入對同步機主系統(tǒng)低頻振蕩的影響，首先建立FSIG及DFIG風(fēng)電機組小干擾分析模型，其次應(yīng)用基于奇異攝動系統(tǒng)動態(tài)降階技術(shù)揭示兩類風(fēng)電機組接入后同步機主系統(tǒng)阻尼變化機理。將奇異攝動動態(tài)降階與特征值分析相結(jié)合，分析不同控制參數(shù)的FSIG及DFIG風(fēng)電機組接入對系統(tǒng)低頻振蕩模式及阻尼的影響?；陲L(fēng)電接入IEEE測試系統(tǒng)仿真表明，所提出的分析方案能夠直觀、有效地對兩類風(fēng)電接入系統(tǒng)低頻振蕩進行分析。

風(fēng)電機組；奇異攝動系統(tǒng)；動態(tài)降階；低頻振蕩；阻尼

0 引言

作為可再生能源的主要利用形式，風(fēng)能發(fā)電一直保持高速增長。2012年我國新增風(fēng)電裝機12 960 MW，同比增長20.8%，總裝機達75 324 MW，裝機總量居全球首位。在華能集團加快裝機結(jié)構(gòu)轉(zhuǎn)型升級，積極發(fā)展清潔能源的宏觀形勢下，2012年底，集團風(fēng)電裝機達9 356.9 MW，已成為僅次于水電的第三大發(fā)電形式。

由欠阻尼導(dǎo)致的低頻振蕩是制約電力系統(tǒng)安全運行的重要因素。大規(guī)模風(fēng)電將對其接入的同步機主系統(tǒng)低頻振蕩特性產(chǎn)生顯著影響。文獻[4]基于數(shù)值仿真，分析了恒速與變速風(fēng)電機組及不同控制方式的雙饋式(Doubly Fed Induction Generator，DFIG)風(fēng)電機組對系統(tǒng)穩(wěn)定性的影響，指出恒速風(fēng)電機組接入系統(tǒng)阻尼特性更好而DFIG機組對系統(tǒng)小干擾穩(wěn)定性能影響不大；若風(fēng)電機組裝設(shè)穩(wěn)定控制器，其接入可改善系統(tǒng)阻尼特性。文獻[5]基于內(nèi)蒙風(fēng)電系統(tǒng)數(shù)值仿真，指出風(fēng)電機組出力改變使網(wǎng)內(nèi)常規(guī)機組運行方式變化，導(dǎo)致一些局部振蕩模態(tài)出現(xiàn)或消失。各大機群間振蕩特性發(fā)生變化，增加了與風(fēng)電場強相關(guān)的振蕩模式。文獻[6]采用DFIG機電暫態(tài)模型，指出風(fēng)電接入對系統(tǒng)小干擾穩(wěn)定性的作用與其并網(wǎng)位置有關(guān)，而風(fēng)電場不參與振蕩。文獻[7]指出定速感應(yīng)型(Fixed-Speed Induction Generator，F(xiàn)SIG)風(fēng)電場接入后增加了與風(fēng)場強相關(guān)的振蕩模式，但這些模式具有很好的阻尼。文獻[8]分析了不同無功控制方案風(fēng)電機組對接入系統(tǒng)穩(wěn)定性影響，指出按電壓偏差進行比例式控制的DFIG電壓調(diào)節(jié)有助于其他節(jié)點電壓恢復(fù)，但過大的反饋增益將造成系統(tǒng)穩(wěn)定性下降。文獻[9-10]針對風(fēng)電接入系統(tǒng)，設(shè)計了附加控制器改善接入系統(tǒng)穩(wěn)定性。數(shù)值分析雖能夠定量計算風(fēng)電對接入系統(tǒng)小干擾穩(wěn)定性的影響，但難以從本質(zhì)上揭示風(fēng)電場與系統(tǒng)小干擾穩(wěn)定性間的相互作用規(guī)律。

首先給出常見FSIG及DFIG風(fēng)電機組接入電力系統(tǒng)簡化小干擾分析模型。其次，采用所建立的分析模型，提出基于奇異攝動的系統(tǒng)動態(tài)降階技術(shù)，通過消除與風(fēng)電機組相關(guān)的快動態(tài)狀態(tài)量，揭示兩類風(fēng)電機組接入后，經(jīng)典同步機系統(tǒng)低頻振蕩阻尼變化的機理。進一步，將奇異攝動動態(tài)降階與特征值分析技術(shù)相結(jié)合，對DFIG風(fēng)電機組控制參數(shù)對系統(tǒng)低頻振蕩模式及阻尼的影響進行研究?；陲L(fēng)電接入IEEE 9節(jié)點系統(tǒng)分析表明，所提出的方案可直觀、有效地分析風(fēng)電接入系統(tǒng)低頻振蕩。

1 風(fēng)電機組小干擾分析模型

1.1 DFIG風(fēng)電機組模型

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

式中：為機械轉(zhuǎn)矩；為轉(zhuǎn)子慣性常數(shù)；為轉(zhuǎn)速；電磁轉(zhuǎn)矩，為內(nèi)電勢，為定子電流向量的共軛。

定子磁場定向DFIG定子有功、無功功率為

式中：為定子電壓向量幅值；、、、為定、轉(zhuǎn)子電流向量、軸分量；為定子磁鏈；定、轉(zhuǎn)子互感=1.5；定、轉(zhuǎn)子自感=+、=+；、為定、轉(zhuǎn)子漏感；為定、轉(zhuǎn)子互感。

由DFIG轉(zhuǎn)子磁鏈、電壓方程及式(2)，可知

式中：為轉(zhuǎn)子電阻；、為轉(zhuǎn)子電壓向量的M、T軸分量；Δ、Δ表達式參見附錄。

對補償后轉(zhuǎn)子電壓控制量為

式中，U，U為轉(zhuǎn)子電壓參考值。

其與轉(zhuǎn)子電流為一階關(guān)系。

由式(5)知，通過配置閉環(huán)系統(tǒng)零極點，可實現(xiàn)有功與無功功率解耦控制。為突出風(fēng)電機組對同步機的影響，簡化風(fēng)電機組動態(tài)模型，小干擾分析時認為DFIG功率動態(tài)為一階慣性環(huán)節(jié)。

由式(2)，單位功率因數(shù)運行的DFIG風(fēng)電機組，其總功率為

式中，_w為風(fēng)電機組的轉(zhuǎn)速。

DFIG小干擾定子有功、無功動態(tài)過程可表示為一階慣性環(huán)節(jié)。

式中：_ws、_ws為定子側(cè)有功；無功功率；、為狀態(tài)變量；、為積分增益。

電網(wǎng)接口方程為

由式(1)、式(6)~式(8)，DFIG線性化模型如式(9)所示。

1.2 FSIG風(fēng)電機組模型

忽略轉(zhuǎn)子暫態(tài)過程，F(xiàn)SIG電磁轉(zhuǎn)矩為

采用一階機械暫態(tài)模型，F(xiàn)SIG線性化模型為

式中：

1.3 同步機組模型

內(nèi)電勢恒定的同步機線性化模型為

式中：；；；；、；、表達式參見附錄。

1.4電力網(wǎng)絡(luò)模型

采用恒阻抗負荷，僅保留同步機與風(fēng)電機組并網(wǎng)節(jié)點，收縮后系統(tǒng)線性化模型為

2 基于動態(tài)降階的風(fēng)電接入系統(tǒng)阻尼分析

2.1奇異攝動系統(tǒng)的動態(tài)降階

風(fēng)電機組具有良好的小干擾穩(wěn)定性，分析中關(guān)注的是其對同步機主系統(tǒng)阻尼的影響。本節(jié)通過動態(tài)降階，消去風(fēng)電機組狀態(tài)量，僅保留同步機狀態(tài)變量，直接分析同步機主系統(tǒng)阻尼變化。

風(fēng)電機組慣性小且控制迅速，其動態(tài)比同步機快的多，含風(fēng)電接入系統(tǒng)為奇異攝動系統(tǒng)。

若存在可微解(,)，滿足

且對小參數(shù)，=(,)，(,(,))在=0處可展開為級數(shù)

(16)

式中：

由式(15)、式(16)有

(17)

由對應(yīng)項系數(shù)，得到代數(shù)方程組

對奇異攝動系統(tǒng)，依次求解代數(shù)方程組(18)，可得(,)的任意階近似。

則奇異攝動系統(tǒng)(14)可降階為

(20)

2.2 DFIG接入系統(tǒng)分析

由式(9)、式(12)及式(13)，消去Δ、Δ、Δ、Δ，含DFIG接入系統(tǒng)方程如式(20)所示。式中各參數(shù)表達式參見附錄。

若選小參數(shù)陣為

式(20)中，對DFIG機組子系統(tǒng)，有

()的一階近似為

(23)

滿足方程

對式(24)，求解方程組(18)，有

故，降階同步機主系統(tǒng)為

(26)

式(26)表明，DFIG對系統(tǒng)同步機阻尼影響與其控制時間常數(shù)及慣性時間常數(shù)均相關(guān)，若控制增益、很大，則風(fēng)電機組不影響同步機阻尼陣。上述結(jié)論解釋如下：對大積分增益，DFIG功率調(diào)節(jié)十分迅速，可忽略調(diào)節(jié)動態(tài)，近似認為動態(tài)過程中其功率恒定。此時，其線性化模型為

由式(12)、式(13)及式(27)，有

式中，表達式參見附錄。

式(29)表明，采用恒功率控制，DFIG電功率調(diào)節(jié)迅速，其接入影響系統(tǒng)同步轉(zhuǎn)矩陣而不影響阻尼陣。DFIG接入不影響系統(tǒng)總阻尼，僅影響總阻尼在各模態(tài)間的分配。DFIG接入后，由于總阻尼守恒，隨著部分模態(tài)阻尼增加，其余模態(tài)阻尼將減小。

2.3 FSIG接入系統(tǒng)分析

由式(11)~式(13)，含F(xiàn)SIG機組接入的電力系統(tǒng)動態(tài)方程為

式中，各參數(shù)參見附錄。

同上節(jié)分析，得到降階同步機主系統(tǒng)方程為

式(31)表明，F(xiàn)SIG對同步機主系統(tǒng)同步轉(zhuǎn)矩陣及阻尼陣均產(chǎn)生影響。機組對系統(tǒng)阻尼影響與風(fēng)電機組慣性與接入系統(tǒng)方式均有關(guān)系，風(fēng)電機組容量越大，對系統(tǒng)影響越顯著。

3 大規(guī)模風(fēng)電接入系統(tǒng)算例分析

本節(jié)通過風(fēng)電機組接入IEEE 3機9節(jié)點系統(tǒng)考查不同類型風(fēng)電機組對系統(tǒng)小干擾穩(wěn)定性的影響。仿真系統(tǒng)如圖1所示，系統(tǒng)運行方式參見附錄。母線8經(jīng)線路L8-10接入不同類型電源。G1為平衡機。風(fēng)電場采用等值風(fēng)電場，單機參數(shù)參見附錄。

不同電源接入后系統(tǒng)特征值如表1與表2所示。易知，原系統(tǒng)存在兩個振蕩模態(tài)，由相應(yīng)特征向量知：模態(tài)1為G2相對G3的局部模態(tài)，模態(tài)2為G1相對G2、G3的區(qū)間模態(tài)，兩個模態(tài)均有良好的阻尼。母線11接入同步機G4后，系統(tǒng)增加一個G3相對G11的局部振蕩模態(tài)，原有各模態(tài)阻尼略有降低，而新增模態(tài)具有良好阻尼。

母線11接入DFIG后，原各模態(tài)頻振蕩率均降低，但幅度不大。模態(tài)1阻尼降低而模態(tài)2阻尼增加。同步機振蕩模態(tài)與無風(fēng)電接入時相同，而風(fēng)電機相應(yīng)右特征向量幅值很小。主要原因是DFIG與系統(tǒng)為柔性連接，低頻振蕩模態(tài)與DFIG轉(zhuǎn)速無關(guān)。

母線11接入FSIG后，系統(tǒng)原有各模態(tài)頻振蕩率與阻尼均有改善，且改善程度較DFIG顯著。同步機振蕩模態(tài)1、2與無風(fēng)電接入時相同。

表1各模態(tài)振蕩頻率與阻尼

Table 1 Frequency and damping ratio of each mode

表2各模態(tài)特征向量

Table 2 Eigenvector of each mode

不同類型電源接入后系統(tǒng)特征值如圖2所示。由該圖知，DFIG通過定子磁場定向?qū)崿F(xiàn)恒功率控制，其與系統(tǒng)為柔性連接，對系統(tǒng)各振蕩模態(tài)頻率及阻尼影響均不顯著，接入后部分模態(tài)阻尼略有增加而其余模態(tài)阻尼降低。而采用經(jīng)典同步機與一階FSIG模型，系統(tǒng)交換功率僅由同步機功角及FSIG轉(zhuǎn)速確定，F(xiàn)SIG功率僅與同步機功角相關(guān)而與同步機轉(zhuǎn)速無關(guān)。FSIG為一隨自身轉(zhuǎn)速變化的可變電阻，故同步機電功率與風(fēng)電機轉(zhuǎn)速相關(guān)，這正是式(36)的矩陣中各非零元意義所在。通過系統(tǒng)電氣聯(lián)系，轉(zhuǎn)速與同步機轉(zhuǎn)速產(chǎn)生關(guān)聯(lián)，從而對系統(tǒng)阻尼產(chǎn)生影響。機組接入后，系統(tǒng)原有各振蕩模態(tài)阻尼略有增加，幅度較機組顯著，其中模態(tài)2阻尼有顯著增加。

圖2不同電源接入系統(tǒng)振蕩模態(tài)

圖3不同控制參數(shù)系統(tǒng)振蕩模態(tài)

DFIG風(fēng)電機組控制增益、變化時，系統(tǒng)特征值情況如圖3。由該圖知，隨著、增大，系統(tǒng)各模態(tài)振蕩頻率變化較小。模態(tài)1阻尼先增大然后減小，而模態(tài)2阻尼先減小后增大，均隨控制參數(shù)非線性變化。系統(tǒng)總阻尼亦呈現(xiàn)先增大后減小的趨勢，最終趨于風(fēng)電機組接入前系統(tǒng)總阻尼值。

對經(jīng)典同步機模型接入系統(tǒng)，系統(tǒng)總阻尼（含非周期振蕩模態(tài)）由同步機機械阻尼確定。使用恒功率控制，若DFIG控制迅速，小干擾條件下其等效于一阻抗，該阻抗影響同步機主系統(tǒng)網(wǎng)絡(luò)阻抗陣，DFIG機組接入將影響系統(tǒng)振蕩頻率，但不改變系統(tǒng)總阻尼，僅影響阻尼在各模態(tài)間的分配?？紤]到DFIG風(fēng)電機功率調(diào)節(jié)十分迅速，其在影響系統(tǒng)振蕩頻率的同時，對系統(tǒng)總阻尼影響不大。

4 結(jié)論

應(yīng)用奇異攝動系統(tǒng)降階技術(shù)，揭示了FSIG與DFIG風(fēng)電機組接入對同步機系統(tǒng)低頻振蕩模態(tài)影響的機理，分析表明：

與同步發(fā)電機組相比，F(xiàn)SIG風(fēng)電機組能夠改善系統(tǒng)小干擾穩(wěn)定性，且對系統(tǒng)小干擾穩(wěn)定性改善效果較DFIG更為顯著。

僅使用最大功率追蹤控制的DFIG不能改善接入系統(tǒng)總阻尼，僅影響系統(tǒng)總阻尼在各振蕩模態(tài)間的分配。

附錄

2 式(9)中表達式

3 式(11)中各參數(shù)表達式

4 式(12)中參數(shù)表達式

其中：Δδ為同步機功角增量；Δω為同步機轉(zhuǎn)速增量；δ為同步機與同步參考坐標系間角度差；Δ、Δ、Δ、Δ為同步旋轉(zhuǎn)坐標系x、y軸電流、電壓分量增量，為內(nèi)電勢，_G為發(fā)電機內(nèi)導(dǎo)納。

5 式(20)中各參數(shù)表達式

式中：×維矩陣，，，；×維矩陣，，，是的分塊。

6式(24)中表達式

7 式(28)中表達式

8 式(30)中各參數(shù)表達式

式中：×維矩陣，，，；×維矩陣，，，是的分塊。

9 仿真系統(tǒng)方式( 1p.u.=100MVA)

功率G1G2G3W4PL5PL6PL8 P0.331.630.850.4-2.3-1.5-1.0 Q0.320.1500-0.5-0.3-0.35

10 DFIG仿真風(fēng)電場機組參數(shù)

設(shè)計數(shù)據(jù)額定風(fēng)速/(m·s-1)切斷風(fēng)速/(m·s-1)空氣密度/(kg·m-3)額定轉(zhuǎn)速/rpm 13.5251.2251 025 DFIG額定功率/MW額定電壓/kVJg/kg·m2rr/p.u. 2.00.69800.01 xr/p.u.xs/p.u.xm/p.u.轉(zhuǎn)速范圍/rpm 0.080.13.0650~1 350 風(fēng)力機風(fēng)輪半徑/m轉(zhuǎn)速范圍/rpmJw/kg·m2 429~161×104 傳動軸及槳距角調(diào)節(jié)Ks/(Nm·rad-1)齒輪箱速比τβ/s變槳范圍/deg 8.1×10777.440.10~18

11 FSIG仿真風(fēng)電場機組參數(shù)

設(shè)計數(shù)據(jù)額定風(fēng)速/(m·s-1)切斷風(fēng)速/(m·s-1)空氣密度/(kg·m-3)額定轉(zhuǎn)速/rpm 12201.2251 520 發(fā)電機額定功率/MW額定電壓/kVJg/(kg·m2)rr/p.u. 1.00.69500.01 xr/p.u.xs/p.u.xm/p.u. 0.050.082.0 風(fēng)力機風(fēng)輪半徑/m轉(zhuǎn)速范圍/rpmJw/kg·m2 30.31200.7×104 傳動軸Ks/(Nm·rad-1)齒輪箱速比 9.1×10776

[1] 郝正航, 余貽鑫. 雙饋風(fēng)電場對電力系統(tǒng)阻尼影響的轉(zhuǎn)矩分析[J]. 電工技術(shù)學(xué)報, 2011, 26(5): 152-157.

HAO Zheng-hang, YU Yi-xin. Torque analysis for the influence of DFIG-based wind farm on damping of power system[J]. Transactions of China Electrotechnical Society, 2011, 26(5): 152-157.

[2] 錢少鋒, 林俐, 沈輝, 等. 基于PSS/E Wind大型風(fēng)電場并入輸電網(wǎng)的動態(tài)特性研究[J]. 電力系統(tǒng)保護與控制, 2009, 37(6): 11-16.

QIAN Shao-feng, LIN Li, SHEN Hui, et al. Dynamic charateristic analysis of transmission grid including large wind farm by PSS/E wind package[J]. Power System Protection and Control, 2009, 37(6): 11-16.

[3] 郭東杰, 王靈梅, 韓西貴, 等. 基于PSS和SSSC的風(fēng) 火打捆輸電系統(tǒng)暫態(tài)穩(wěn)定性研究[J]. 電力系統(tǒng)保護與控制, 2012, 40(19): 61-65.

GUO Dong-jie, WANG Ling-mei, HAN Xi-gui, et a1. Study on transient stability of wind & thermal power transmission system based on PSS and SSSC[J]. Power System Protection and Control, 2012, 40(19): 61-65.

[4] ANAYA-LARA O, HUGHES F M, JENKINS N, et al. Influence of wind farms on power system dynamic and transient stability[J]. Wind Engineering, 2009, 30(2): 107-127.

[5] 張紅光, 張粒子, 陳樹勇, 等. 大容量風(fēng)電場對電力系統(tǒng)小干擾穩(wěn)定和阻尼特性的影響[J]. 電網(wǎng)技術(shù), 2007, 31(13): 75-80.

ZHANG Hong-guang, ZHANG Li-zi, CHEN Shu-yong, et a1. Influence of large scale wind farms on power system small signal stability and damping characteristic[J]. Power System Technology, 2007, 31(13): 75-80.

[6] 王忱, 石立寶, 姚良忠, 等. 大規(guī)模雙饋型風(fēng)電場的小擾動穩(wěn)定分析[J]. 中國電機工程學(xué)報, 2010, 30(4): 63-70.

WANG Chen, SHI Li-bao, YAO Liang-zhong, et al. Small signal stability analysisi of the large scale wind farm with DFIGs[J]. Proceedings of the CSEE, 2010, 30(4): 63-70.

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

GUAN Hong-liang, CHI Yong-ning, DAI Hui-zhu, et al. Small signal stability and control of wind tubine with asynchronous generator intergration into power system[J]. Automation of Eletric Power Systems, 2008, 32(4): 54-58.

[8] FERNANDEZ R D, BATTAIOTTO P E, MANTZ R J. Impact of wind farms voltage regulation on the stability of the network frequency[J]. International Journal of Hydrogen Energy, 2008, 33(13): 3543-3548.

[9] 張子泳, 胡志堅, 李勇匯, 等. 大型雙饋風(fēng)力發(fā)電系統(tǒng)小信號動態(tài)建模及附加阻尼控制器設(shè)計[J]. 電力系統(tǒng)保護與控制, 2011, 39(18): 127-133.

ZHANG Zi-yong, HU Zhi-jian, LI Yong-hui, et a1. Small signal dynamic modelling and additionaller damping controller designing for large wind generationsystem based on DFIG[J]. Power System Protection and Control, 2011, 39(18): 127-133.

[10]范高鋒, 遲永寧, 趙海翔, 等. 用STATCOM提高風(fēng)電 場暫態(tài)電壓穩(wěn)定性[J]. 電工技術(shù)學(xué)報, 2007, 22(11): 158-162.

FAN Gao-feng, CHI Yong-ning, ZHAO Hai-xiang, et al. Transient voltage stability enhancement of wind farm using STATCOM[J]. Transactions of China Electrotechnical Society, 2007, 22(11): 158-162.

[11] 王毅, 張祥宇, 李和明, 等. 永磁直驅(qū)風(fēng)電機組對系統(tǒng)功率振蕩的阻尼控制[J]. 電工技術(shù)學(xué)報, 2012, 27(12): 162-170.

WANG Yi, ZHANG Xiang-yu, LI He-ming, et al. Damping control of PMSG base wind trubines for power system osillations[J]. Transactions of China Electrotechnical Society, 2012, 27(12): 162-170.

[12]徐光虎. 交直流輸電系統(tǒng)低頻振蕩相關(guān)問題分析及阻尼控制策略[D]. 上海: 上海交通大學(xué), 2005.

XU Guang-hu. Special issues of low frequency osillation and damping control strategies study for AC/DC power system[D]. Shanghai: Shanghai Jiaotong University, 2005.

System damping analysis after wind power integration using singular perturbation system reduction

SHU Jin, DU Jin-song, WANG Jing-cheng, ZHANG Bao-hui

(1. Thermal Power Research Institute, Co., Ltd. Xi’an 710043, China; 2. State Key Laboratory of Electrical Insulation and Power Equipment, Xi’an Jiaotong University, Xi’an 710049, China)

In order to analyze the effect of wind power integration on low frequency oscillation of synchronous machine main system, the small-signal analysis models of FSIG and DFIG are built firstly, then the damping change mechanism after the above two generators integration is revealed by using the singularly perturbed system reduction. The influence of FSIG and DFIG integration with different parameters on system oscillation modes and damping is analyzed combined with the singularly perturbed system reduction and eigenvalue analysis. The IEEE testing system simulation shows that the proposed scheme can analyze the system low-frequency oscillation after wind power integration intuitively and effectively.

wind turbine generator; singularly perturbed system; system reduction; low frequency oscillation; damping

TM73

A

1674-3415(2014)19-0018-08

2013-12-06；

2014-04-25

舒 進（1983-），男，工學(xué)博士，工程師，主要從事風(fēng)力發(fā)電相關(guān)研究工作。E-mail：shujin@tpri.com.cn