用于風(fēng)力發(fā)電系統中的仿真技術(shù)

0 引言

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

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

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

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

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

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

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

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

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

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

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

2 風(fēng)能及風(fēng)電機組仿真模型

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

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

2.1 典型風(fēng)力發(fā)電機組的仿真模型總體結構

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

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

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

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

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

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

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

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

下載附件 (16.06 KB) 圖 3 恒速風(fēng)力機仿真模型的總體結構

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

下載附件 (29.28 KB) 圖 4 變速風(fēng)力機仿真模型的總體結構

2.2 風(fēng)能特性模型

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

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

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

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

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

2.3 風(fēng)力機模型[4]

風(fēng)力機實(shí)際能夠獲得的機械功輸出為：

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

下載附件 (3.35 KB) (1)

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

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