用于風力發(fā)電系統(tǒng)中的仿真技術

0 引言

出于保護環(huán)境的考慮以及全球面臨的能源短缺現(xiàn)狀，風力發(fā)電在世界范圍內(nèi)得到了快速發(fā)展。隨著風電行業(yè)的技術進步，風力發(fā)電成本逐步降低，在經(jīng)濟性上已經(jīng)能夠與核能發(fā)電、水力發(fā)電展開競爭。當前，我國面臨電力短缺局面，在煤電占主導地位的我國電力行業(yè)，因環(huán)境承載力限制以及各種因素導致的煤炭短缺局面，煤電發(fā)展受到制約。而我國風能資源豐富，風能利用得到了政府的政策支持，風力發(fā)電產(chǎn)業(yè)面臨前所未有的發(fā)展機遇。近幾年來，我國風電產(chǎn)業(yè)發(fā)展勢頭良好，多個大型風電場處于建設或規(guī)劃階段。

我國風電產(chǎn)業(yè)發(fā)展尚存在諸多制約因素。就技術層面而言，大型風力發(fā)電設備生產(chǎn)技術不過關，大多從國外采購或引進技術生產(chǎn)，研發(fā)基礎薄弱;對風電機組或風電場的運行特性的研究不足，設備運行管理水平還有待提高。將仿真技術廣泛應用于風力發(fā)電設備的設計、試驗測試、運行分析等各個方面，將有助于加快我國風力發(fā)電技術的發(fā)展步伐，是縮小與發(fā)達國家技術差距的捷徑。

1 仿真技術在風力發(fā)電系統(tǒng)的應用概況

隨著風力發(fā)電在世界的廣泛應用，為降低風力發(fā)電成本、提高風能利用效率，風力發(fā)電設備單機容量越來越大，同時為風機的設計制造、控制系統(tǒng)設計和運行等各個方面提出了更多的研究課題。傳統(tǒng)的實物測試研究方法已不能滿足發(fā)展需要，仿真技術因不受氣象條件的限制，且投入低等優(yōu)點，逐漸在風力發(fā)電機組的研究和測試領域得到了越來越廣泛的應用。

仿真即選取一個物理的或抽象的系統(tǒng)的某些行為特征，用數(shù)學模型來表示它們的過程，若用計算機求解數(shù)學模型，稱為計算機仿真。通俗說來，仿真是指使用儀器設備、模型、多媒體技術，以及利用場地、環(huán)境的布置，模仿出真實系統(tǒng)的工作特性和環(huán)境，進而用于科學研究、工業(yè)設計、預測預報或教學訓練等目的的一項綜合技術。仿真若僅限于設計研究目的，則勿需仿真對象系統(tǒng)的環(huán)境，亦無實時仿真的必要，借助一臺主流微型計算機和商業(yè)仿真軟件即可開展仿真研究工作。

1980 年代初，國外學者開始將仿真方法用于風電機組的的性能研究[1] ，其后，仿真技術在風力發(fā)電系統(tǒng)的應用范圍逐漸拓展。目前，從風電關鍵設備和控制系統(tǒng)的設計、制造、性能測試與研究，風電機組或風電場運行分析等各個方面均有仿真技術的應用。仿真技術的應用在很大程度上替代了傳統(tǒng)的利用實際設備開展的設計檢驗等手段。主要的研究方向整理如下。

1)風能特性仿真，反映風能的位置分布和時間變化特性。風特性仿真結(jié)果將用于風力發(fā)電機組或風電場的仿真分析中，是風電仿真研究的基礎。

2)風力發(fā)電機組仿真，仿真特定風力機組在風能變化下輸出電能的變化規(guī)律，分析其特性，尋找設備本身存在的不足，提供改進建議。表征電能特性的參數(shù)主要包括有功功率、無功功率、電壓和頻率。

3)控制系統(tǒng)仿真，建立待檢驗的控制系統(tǒng)的仿真模型和控制對象的仿真模型，建立模型間的相互聯(lián)系。改變仿真風電機組的風能參數(shù)或工作狀態(tài)，測試在各種不同運行方式下控制系統(tǒng)的動作特性和工作效果，尋找控制系統(tǒng)設計中存在的問題，改進設計后修正仿真模型進一步驗證，直到控制系統(tǒng)滿足設計和運行要求。

4)風電場仿真，針對特定風電場的具體風能特性和實際(或規(guī)劃設計)安裝的風電機組情況，建立整個風電場的仿真模型。研究風能變化、風機介入或退出系統(tǒng)對風電場電能特性的影響，進而分析風電場建設的可行性，分析風電場不同運行方式下對電力系統(tǒng)的影響，或用于運行人員培訓，提高風電場運行管理水平。

國內(nèi)利用仿真技術開展風電系統(tǒng)研究的起步較晚，公開發(fā)表的仿真研究成果不多，尚未形成氣候。近幾年情況在發(fā)生變化，更多的研究人員已經(jīng)將仿真技術引入風電系統(tǒng)的研究工作中，相信近期內(nèi)將會有更多的高水平仿真研究成果發(fā)表，并能有力促進風電產(chǎn)業(yè)的技術進步。

2 風能及風電機組仿真模型

與常規(guī)發(fā)電機組(如火電、水電、核電)相比，風力發(fā)電機組的突出特點是輸入能量不受控制，這一特點導致風力發(fā)電機組在構(gòu)成上與常規(guī)發(fā)電機組有著很大的不同且呈現(xiàn)出多樣化特點。我們知道，常規(guī)發(fā)電機組的機械能-電能轉(zhuǎn)換裝置普遍采用同步發(fā)電機，而并網(wǎng)型風電機組采用的發(fā)電機則形式多樣，如恒速恒頻同步/異步發(fā)電機、交/直/交發(fā)電機、磁場調(diào)制發(fā)電機、交流勵磁雙饋發(fā)電機等。因采用的發(fā)電機類型不同，相應的控制系統(tǒng)區(qū)別很大，電能參數(shù)隨風能變化的特性也有很大的不同。

仿真研究人員需要根據(jù)風力發(fā)電機組的特點開發(fā)針對性的仿真模型軟件。限于篇幅，本文主要介紹共性部分的仿真。

2.1 典型風力發(fā)電機組的仿真模型總體結(jié)構(gòu)

在風電場中得到廣泛應用的恒速風力機如圖1 所示[2]，異步發(fā)電機將風輪吸收的機械能轉(zhuǎn)化成電能，發(fā)電機轉(zhuǎn)速隨發(fā)電量的變化而在一定范圍內(nèi)變化，因轉(zhuǎn)速變化范圍很小(1% 左右)，通常稱為恒速系統(tǒng)。恒速系統(tǒng)通常選用失速型調(diào)節(jié)方式。

2012-10-6 09:37:50 上傳

下載附件 (13.24 KB) 圖 1 恒速風力機系統(tǒng)示意圖

一種典型的變速風力發(fā)電機組見圖2，它采用雙饋異步發(fā)電機(DFIG)。發(fā)電機的定子線圈直接與電網(wǎng)相連，轉(zhuǎn)子線圈則通過滑環(huán)和電力電子逆變器與電網(wǎng)連接。因此，當風速變化引起發(fā)電機轉(zhuǎn)速變化時，通過控制轉(zhuǎn)子電流的頻率，可保持定子頻率的恒定，進而實現(xiàn)風力發(fā)電機組的變速運行。在高風速條件下，通過調(diào)整葉片槳距限制風力機的輸出功率。

上述兩種風力發(fā)電機組的仿真模型的總體結(jié)構(gòu)分別表示在圖3 和圖4 中[3] ，變速風力機的控制系統(tǒng)要比恒速系統(tǒng)復雜得多，其仿真模型相應增加了槳距角、轉(zhuǎn)速、端電壓等控制器子模型和變頻器的仿真模型。

2012-10-6 09:37:49 上傳

下載附件 (16.94 KB) 圖 2 典型的雙饋發(fā)電機組系統(tǒng)示意圖

2012-10-6 09:37:48 上傳

下載附件 (16.06 KB) 圖 3 恒速風力機仿真模型的總體結(jié)構(gòu)

2012-10-6 09:37:47 上傳

下載附件 (29.28 KB) 圖 4 變速風力機仿真模型的總體結(jié)構(gòu)

2.2 風能特性模型

描述風能特性的參數(shù)主要有風速、風向和風密度。風的密度主要取決于風機所處的地理位置，氣候變化也會產(chǎn)生一定影響，對于特定風機而言，風密度可以直接取自測量數(shù)據(jù)，并可以忽略密度的變化;針對研究型的仿真應用，風向的變動可不予考慮，即假定風力機一直跟蹤風向的變化。因此，我們主要關注風速的變化特性。

風因大氣環(huán)流形成，風速是一個典型的隨機變量。若不考慮風的方向性，風速是其空間坐標位置和時間的函數(shù)，即v=f(x,y,z,t) 。我們將描述某一區(qū)域風速的空域、時域分布變化特性的模型又稱為風場模型(Wind Field Model)。嚴格說來，各空間位置上的風速因風的隨機性、風場地形等因素影響而各不相同，因此，要建立一個準確的風場模型幾乎是不可能的，需要進行一定的簡化處理。

若風場的地形相對平坦，周邊空曠，則基本可以認為在同一高度層上整個風場內(nèi)各點的風速是相同的，這樣可以將風場風速的三維空間模型簡化為沿高度方向變化的一維模型。對于空間分布廣，且地形復雜的大型風電場，可以將整個風場劃分成幾個區(qū)域，針對不同區(qū)域的風能特點建立簡化的一維空間模型，形成分段集總式一維模型。風速空域模型轉(zhuǎn)化為研究風速沿地平面高度方向的變化規(guī)律，借助空氣動力學理論和風場測量數(shù)據(jù)，模型不難建立。

在時間維度上，大時期尺度(小時、天)的風速變化范圍很大，且沒有規(guī)律可循，只能根據(jù)風場監(jiān)測記錄數(shù)據(jù)擬合出風速變化模型。對于絕大多數(shù)仿真應用而言，我們不太關心大時間尺度的風速變化，而重點關注小時間尺度上的風速變化特性。在小時間尺度上觀察，風速隨時間的變化呈現(xiàn)出脈動變化的特點，即風速均值在一段時間內(nèi)基本不變，風速在均值附近波動，國內(nèi)外學者據(jù)此提出了各種描述風頻分布的方法，如概率分布模型、瑞利分布模型、對數(shù)正態(tài)分布模型等[1] 。

需要說明的是，在建立風電場內(nèi)各風力發(fā)電機組的仿真模型時，需要考慮到風力機的尾流效應，即上游風力機對下游風力機流入風速的影響，影響關系和程度取決于風向、風速和風機安裝位置關系，此時，風力機的輸出機械能通常由尾流系數(shù)予以校正。

2.3 風力機模型[4]

風力機實際能夠獲得的機械功輸出為：

2012-10-6 09:37:46 上傳

下載附件 (3.35 KB) (1)

式中， R 、ρ、 v 和Cp分別為風輪半徑、空氣密度、風速和風能利用系數(shù)。

Cp代表風力機能夠從風能中提取出機械能的程度，它取決于風力機葉片的結(jié)構(gòu)和運行狀態(tài)，其數(shù)值由風機廠家提供。Cp主要是葉尖速比λ和槳距角β的函數(shù)，即 ：Cp=f(λ，β)，對于投入運行的風力機，葉片的潔凈程度對Cp的影響很大，譬如葉片結(jié)冰、污物聚集等會改變?nèi)~片的氣動外型，進而降低風能利用系數(shù)的數(shù)值。