劉勇
(大連供電公司,遼寧省大連市 116001)
隨著能源危機(jī)和溫室效應(yīng)等問題的日趨凸顯,利用清潔無污染的風(fēng)能發(fā)電成為一種有效的解決措施。從2006年的《可再生能源法》頒布以來,我國的風(fēng)電事業(yè)得到了快速發(fā)展[1-2],大型風(fēng)電場接連涌現(xiàn)。研究整個(gè)風(fēng)電場的并網(wǎng)運(yùn)行與電網(wǎng)的相互影響,以及系統(tǒng)發(fā)生短路故障情況下風(fēng)電場的暫態(tài)運(yùn)行特性等,對(duì)整個(gè)電力系統(tǒng)的安全穩(wěn)定運(yùn)行非常重要[3-4]。本文從單臺(tái)風(fēng)電機(jī)組模型入手,研究風(fēng)電場的等值方法,搭建包含風(fēng)電場的電力系統(tǒng)模型,分析電網(wǎng)發(fā)生三相短路和單相短路故障時(shí)并網(wǎng)風(fēng)電場的暫態(tài)特性。
風(fēng)電場等值方法有參數(shù)變換單機(jī)等值法[5]、變尺度降階多機(jī)等值法[6]、容量加權(quán)單機(jī)等值法[7]以及在此基礎(chǔ)上的改進(jìn)容量加權(quán)單機(jī)等值法[8-9]等。
假定某風(fēng)電場由n臺(tái)雙饋風(fēng)電機(jī)組構(gòu)成,且滿足以下條件[8]:(1)風(fēng)電場中的風(fēng)電機(jī)組類型單一且出口變壓器接在同一條母線上;(2)忽略連接相鄰2臺(tái)風(fēng)電機(jī)組的線路阻抗。由此可對(duì)包含這n臺(tái)風(fēng)電機(jī)組的風(fēng)電場進(jìn)行等效處理[10-11]。由于雙饋電機(jī)也是一種異步電機(jī),所以可以借鑒異步電動(dòng)機(jī)作為負(fù)載時(shí)的等值方法,研究風(fēng)電場的并網(wǎng)運(yùn)行。接在同一母線的n臺(tái)雙饋風(fēng)電機(jī)組的等效過程如圖1所示。
圖1(a)為單臺(tái)異步電動(dòng)機(jī)的T型等值電路,根據(jù)電路等效變換原理,引入阻抗Zf,將圖1(a)等效變換為圖1(b),有
圖1 參數(shù)變換法Fig.1Parameter transformation method
式中:Zs為定子阻抗;Zr為轉(zhuǎn)子阻抗;Zm為勵(lì)磁阻抗。當(dāng)風(fēng)電場的n臺(tái)風(fēng)電機(jī)組都接在同一母線上時(shí),每臺(tái)風(fēng)電機(jī)組經(jīng)過等效變換后的電路如圖1(c)所示,化簡后如圖1(d)所示,可將圖1(d)的電路變換成如圖1(e)所示的T型等效電路,于是得到了并聯(lián)在同一母線上的n臺(tái)風(fēng)電機(jī)組的單機(jī)等值模型。圖1(a)~(e)的參數(shù)分別為
式中:Zs、Zr、Zm分別為定子阻抗、轉(zhuǎn)子阻抗和勵(lì)磁阻抗;Z's、Z'r、Z'm分別為等效模型的定子阻抗、轉(zhuǎn)子阻抗和勵(lì)磁阻抗。
風(fēng)電場的容量、風(fēng)力機(jī)的風(fēng)能利用系數(shù)計(jì)算則采用改進(jìn)的容量加權(quán)風(fēng)電場單臺(tái)等值方法[8]。設(shè)單臺(tái)風(fēng)機(jī)的等值容量為S'和風(fēng)能利用系數(shù)為C',則有
式中Si、Ci分別為第i臺(tái)風(fēng)機(jī)容量和風(fēng)能利用系數(shù)。
對(duì)于風(fēng)電機(jī)組其他參數(shù)X'的等值方法,則要引入加權(quán)系數(shù)ξi,為
考慮到風(fēng)電場的尾流效應(yīng),假定風(fēng)電場布置在地域開闊、地形變化比較平緩的區(qū)域,選用典型的風(fēng)速尾流效應(yīng)Jensen模型,為
式中:V1為第1排(列)的風(fēng)速;Vn為第n排(列)的風(fēng)速;Ct為風(fēng)力機(jī)推力系數(shù);D為風(fēng)力機(jī)直徑;X為每排(列)相鄰風(fēng)力機(jī)間的軸向距離;K為尾流衰減常數(shù),K=0.5/ln(h/z0),其中h為風(fēng)力機(jī)的輪轂高度,z0為粗糙度,z0=2 mm。
當(dāng)風(fēng)力發(fā)電機(jī)組呈N列M行規(guī)則分布,即可以通過(i,j)進(jìn)行標(biāo)記,其中i代表風(fēng)電機(jī)組所在的列,j代表風(fēng)電機(jī)組所在的行。在風(fēng)電場內(nèi)某行或某列機(jī)組參數(shù)完全相同時(shí),可采用參數(shù)變換法將其等值為單機(jī)模型,為
式中:下標(biāo)i和j為風(fēng)電場內(nèi)第j行第i臺(tái)風(fēng)力發(fā)電機(jī)參數(shù);標(biāo)示為第j行其對(duì)應(yīng)的等值機(jī)組參數(shù)。
風(fēng)電場包括同一型號(hào)的16臺(tái)1.5 MW和10臺(tái)2 MW雙饋式風(fēng)力發(fā)電機(jī)組。搭建的仿真模型如圖2所示。
圖2 仿真模型Fig.2Simulation Model
圖2中G1為帶有自動(dòng)勵(lì)磁調(diào)節(jié)裝置的同步發(fā)電機(jī),額定容量為10 MVA,額定電壓為10.8 kV;G2為理想電壓源,額定電壓為110 kV,等效阻抗Xs=4 Ω;G3為等值風(fēng)電場,額定容量為44 MVA;S1為RLC構(gòu)成的三相對(duì)稱負(fù)荷,負(fù)荷功率為15+j3 MVA;S2為包含電動(dòng)機(jī)的小型工業(yè)負(fù)荷,負(fù)荷功率為8.8+j1 MVA。
母線B5在5.0 s時(shí)發(fā)生三相接地短路,接地電阻為1 Ω,持續(xù)時(shí)間為0.1 s,風(fēng)電場輸出有功功率和無功功率如圖3所示。
圖3 風(fēng)電場輸出有功功率和無功功率Fig.3Active power and reactive power of wind farm output
圖3中Pwind、Qwind和分別為正常運(yùn)行和發(fā)生故障時(shí)風(fēng)電場輸出的有功功率和無功功率。在5 s時(shí)發(fā)生三相短路故障,風(fēng)電場的有功輸出降低,從16 MW降到10 MW;5.1 s時(shí)故障消除,有功輸出增加至20 MW,比非故障運(yùn)行時(shí)高出約4 MW;故障消除后有功輸出出現(xiàn)短暫的波動(dòng)。發(fā)生三相故障時(shí),風(fēng)電場從電網(wǎng)中吸收的無功功率大幅增加,在短路瞬間,吸收的無功功率從5 Mvar驟升至10 Mvar,然后逐漸減少;到5.1 s時(shí)故障消除,吸收的無功功率減少至 0;故障消除后,風(fēng)電場從電網(wǎng)中吸收的無功上升,且吸收的無功出現(xiàn)波動(dòng),在故障消除后的0.4 s時(shí),無功吸收最多,再過0.4 s后無功吸收最少。線路L1的電流波形如圖4所示。
圖4 線路L1電流Fig.4Current of line L1
圖4中IL1、I(3)L1為分別為正常運(yùn)行和發(fā)生三相短路故障時(shí)的線路L1電流。發(fā)生三相接地時(shí)線路L1的電流急劇增加,故障期間,線路L1電流從350 A上升到500 A,最大值為非故障時(shí)的2倍。故障消除后線路L1電流有所波動(dòng)。母線B1的電壓波形如圖5所示。
圖5中UB1、U(3)B1分別為正常運(yùn)行和發(fā)生三相短路故障時(shí)的母線B1電壓。發(fā)生三相接地時(shí)母線B1的電壓從1.05 pu驟降至0.55 pu,在5.1 s時(shí)跌落至0.4 pu。故障消除后電壓回升并出現(xiàn)波動(dòng),在5.6 s時(shí)電壓最低至0.85 pu,在6 s時(shí)電壓最高至1.1 pu。正常運(yùn)行和發(fā)生三相短路時(shí)的系統(tǒng)頻率波形如圖6所示。
圖5 母線B1電壓Fig.5Voltage of bus B1
圖6 頻率Fig.6Frequency
圖6中f、f(3)分別為正常運(yùn)行和發(fā)生三相短路時(shí)的系統(tǒng)頻率。在5 s時(shí)發(fā)生三相接地故障,系統(tǒng)頻率從50.1 Hz開始下降;在5.1 s時(shí)頻率減少至48.4 Hz;在5.1 s時(shí)故障消除,頻率上升并進(jìn)入暫態(tài)波動(dòng)過程;在5.2 s時(shí)出現(xiàn)波峰,頻率為52.4 Hz;在5.3 s時(shí)出現(xiàn)波谷,頻率跌至49.7 Hz。
母線B5在5.0 s時(shí)發(fā)生A相接地,接地電阻為1 Ω,持續(xù)時(shí)間為0.1 s,風(fēng)電場輸出的有功功率和無功功率如圖7所示。
圖7 風(fēng)電場輸出有功功率和無功功率Fig.7Active power and reactive power of wind farm output
圖7中Pwind、Qwind和P(1)wind、Q(1)wind分別為正常運(yùn)行和單相接地時(shí)風(fēng)電場輸出的有功功率和無功功率。在5 s時(shí)發(fā)生單相接地,風(fēng)電場的有功輸出略有減少;在5.1 s時(shí)故障消除,有功輸出迅速上升;在5.2 s時(shí)有功功率曲線和非故障運(yùn)行時(shí)的風(fēng)電場有功輸出曲線基本重合。故障時(shí)風(fēng)電場從電網(wǎng)中吸收的無功功率略有增加,增加幅度很小;在5.1 s時(shí)故障消除,吸收的無功功率減少至4 Mvar;故障消除后,風(fēng)電場從電網(wǎng)中吸收的無功功率出現(xiàn)很小的波動(dòng)。線路L1的電流波形如圖8所示。
圖8 線路L1電流Fig.8Current of line L1
圖8中IA為正常運(yùn)行時(shí)線路L1的A相電流分別為發(fā)生單相接地時(shí)線路L1的A相、B相、C相電流。A相發(fā)生接地短路故障時(shí),線路L1的C相電流升高到400 A,A相電流上升到340 A,B相電流變?yōu)?40 A。與非故障時(shí)的電流相比,A相、C相電流上升,B相電流略有減少。母線B1的電壓波形如圖9所示。
圖9 母線B1電壓Fig.9Voltage of bus B1
圖9中U為正常運(yùn)行時(shí)的母線B1的電壓分別為單相接地時(shí)母線B1出口的連接線路A相、B相、C相電壓。在5 s時(shí)發(fā)生單相接地故障,35 kV母線B1的故障相(A相)電壓從1.05 pu驟降至0.75 pu;故障期間電壓繼續(xù)下降,在5.1 s時(shí)跌落到0.7 pu;B相電壓和C相電壓也發(fā)生跌落,在5.1 s時(shí)分別降至0.9和0.75 pu。故障消除后,三相電壓回升并出現(xiàn)較小的波動(dòng)。風(fēng)電場系統(tǒng)頻率如圖10所示。
圖10中f、f(1)分別為正常運(yùn)行和發(fā)生單相接地時(shí)的系統(tǒng)頻率。在5 s時(shí)發(fā)生單相接地,系統(tǒng)頻率從50.1 Hz開始下降;在5.1 s時(shí)頻率減少至49.7 Hz;在5.1 s時(shí)故障消除,頻率上升并進(jìn)入暫態(tài)波動(dòng)過程;在5.2 s時(shí)出現(xiàn)頻率的最大值,為50.6 Hz。
圖10 系統(tǒng)頻率Fig.10System frequency
當(dāng)電網(wǎng)發(fā)生短路時(shí),風(fēng)電場的功率輸出、并網(wǎng)點(diǎn)母線電壓、故障線路電流以及系統(tǒng)頻率都會(huì)發(fā)生變化。短路時(shí)風(fēng)電場的有功輸出驟減,從電網(wǎng)吸收的無功驟增,并網(wǎng)點(diǎn)母線電壓驟減;故障線路的電流變化和故障類型有關(guān),并非故障線路的短路電流最大;短路時(shí)系統(tǒng)頻率驟減,故障消除后會(huì)發(fā)生波動(dòng),短路程度越大、波動(dòng)幅值越大。
[1]李俊峰,施鵬飛,高虎.中國風(fēng)電發(fā)展報(bào)告[M].???海南出版社,2010.
[2]肖湘寧.電能質(zhì)量分析與控制[M].北京:中國電力出版社,2004:71-72.
[3]李建林,許洪華.風(fēng)力發(fā)電系統(tǒng)低電壓運(yùn)行技術(shù)[M].北京:機(jī)械工業(yè)出版社,2008:130-135.
[4]Shafiu A,Anaya-Lara O,Bathurst G,et al.Aggregated wind turbine models for power system dynamic studies[J].Wind Engineering,2006,30(3):171-186.
[5]孫建峰,焦連偉,吳俊玲,等.風(fēng)電場發(fā)電機(jī)動(dòng)態(tài)等值問題研究[J].電網(wǎng)技術(shù),2004,28(7):56-61.
[6]吳學(xué)光,張學(xué)成,印永華,等.異步風(fēng)力發(fā)電系統(tǒng)動(dòng)態(tài)穩(wěn)定性分析的數(shù)學(xué)模型及其應(yīng)用[J].電網(wǎng)技術(shù),1998,22(6):68-72.
[7]Morren J,Haan S W H d.Short-circuit current of wind turbines with turbines with doubly fed induction generator[J].IEEE Transactions on Energy conversion,2007,22(1):174-180.
[8]湯宏,吳俊玲,周雙喜.包含風(fēng)電場電力系統(tǒng)的小干擾穩(wěn)定分析建模和仿真[J].電網(wǎng)技術(shù),2004,28(1):38-41.
[9]魏巍,王渝紅,李興源,等.基于PSASP的雙饋風(fēng)電場建模及接入電網(wǎng)仿真[J].電力自動(dòng)化設(shè)備,2009,29(12):68-73.
[10]袁貴川,倪林,陳穎.大連風(fēng)電場等值建模與仿真分析[J].廣東電力,2010,23(11):6-9.
[11]卞松江,呂曉美,相會(huì)杰,等.交流勵(lì)磁變速恒頻風(fēng)力發(fā)電系統(tǒng)控制策略的仿真研究[J].中國電機(jī)工程學(xué)報(bào),2005,25(16): 57-62.
(編輯:張磊)