代元軍 ， 汪建文 ， 劉雄飛 ， 溫彩鳳 ，4， 車(chē) 飛

（1.內(nèi)蒙古工業(yè)大學(xué)能源動(dòng)力學(xué)院，內(nèi)蒙古 呼和浩特 010051；2.內(nèi)蒙古可再生能源實(shí)驗(yàn)室，內(nèi)蒙古 呼和浩特 010051；3.新疆工業(yè)高等專(zhuān)科學(xué)校，新疆 烏魯木齊 830091；4.內(nèi)蒙古呼和浩特職業(yè)學(xué)院，內(nèi)蒙古 呼和浩特 010051）

0 引言

目前在石油、煤炭等傳統(tǒng)能源緊張的環(huán)境下，世界各國(guó)均將目光投向新能源領(lǐng)域，其中風(fēng)能倍受關(guān)注。風(fēng)力機(jī)葉片作為將風(fēng)能轉(zhuǎn)換為動(dòng)能的“心臟”，其設(shè)計(jì)直接決定著風(fēng)力機(jī)組的效率[1]，而翼型又是葉片的重要組成部分，因此，對(duì)風(fēng)力機(jī)翼型的設(shè)計(jì)是當(dāng)前風(fēng)力機(jī)空氣動(dòng)力學(xué)研究的一個(gè)熱點(diǎn)，也是風(fēng)力機(jī)研究的基礎(chǔ)[2-5]。該文對(duì)采用不同翼型設(shè)計(jì)的風(fēng)力機(jī)在風(fēng)洞進(jìn)行了氣動(dòng)特性實(shí)驗(yàn)，得到了風(fēng)力機(jī)風(fēng)速功率特性曲線和風(fēng)速風(fēng)能利用系數(shù)曲線，然后利用中國(guó)氣象科學(xué)數(shù)據(jù)共享網(wǎng)提供的當(dāng)?shù)仫L(fēng)資源數(shù)據(jù)和WAsP軟件在當(dāng)?shù)貙?shí)際風(fēng)場(chǎng)的數(shù)字地圖上，建立模擬風(fēng)力機(jī)站[6]，求解風(fēng)力機(jī)的年發(fā)電量和單位面積年發(fā)電量，用以分析和評(píng)價(jià)不同翼型風(fēng)力機(jī)的性能好壞。

1 風(fēng)力機(jī)氣動(dòng)特性實(shí)驗(yàn)

1.1 實(shí)驗(yàn)設(shè)備

實(shí)驗(yàn)設(shè)備主要由低速開(kāi)口直流式風(fēng)洞、9033G測(cè)試模塊和風(fēng)力機(jī)等組成，如圖1所示。實(shí)驗(yàn)風(fēng)洞為B1/K2直流式低速風(fēng)洞，分為開(kāi)口段和閉口段，采用數(shù)字變頻技術(shù)調(diào)節(jié)風(fēng)速，變頻器的變頻范圍為0~60.0 Hz。風(fēng)洞開(kāi)口實(shí)驗(yàn)段的圓形橫截面直徑為2.04 m，風(fēng)速調(diào)節(jié)范圍為0~18m/s。測(cè)試樣機(jī)1#風(fēng)力機(jī)是某公司生產(chǎn)的400W風(fēng)力機(jī)，2#風(fēng)力機(jī)是該課題組自行設(shè)計(jì)制造的采用新翼型S809的300 W風(fēng)力機(jī)，進(jìn)行氣動(dòng)特性實(shí)驗(yàn)，其參數(shù)見(jiàn)表1。

圖1 實(shí)驗(yàn)設(shè)備

表1 風(fēng)力機(jī)參數(shù)

1.2 實(shí)驗(yàn)方法

實(shí)驗(yàn)按照GB/T 19068.3-2003《離網(wǎng)型風(fēng)力發(fā)電機(jī)組第3部分：風(fēng)洞實(shí)驗(yàn)方法》[7]和GB/T 19068.2-2003《離網(wǎng)型風(fēng)力發(fā)電機(jī)組第2部分：試驗(yàn)方法》[8]，風(fēng)洞試驗(yàn)在內(nèi)蒙古工業(yè)大學(xué)的B1/K2低速風(fēng)洞出口段進(jìn)行。該實(shí)驗(yàn)具體采用的實(shí)驗(yàn)方法是對(duì)兩種不同翼型的風(fēng)力機(jī)，改變來(lái)流風(fēng)速，改變電阻大小，得出不同尖速比下風(fēng)力機(jī)的輸出功率和功率系數(shù)，繪制風(fēng)速功率特性曲線和風(fēng)速功率系數(shù)曲線。

1.3 實(shí)驗(yàn)結(jié)果及分析

由實(shí)驗(yàn)得到的2個(gè)風(fēng)力機(jī)的風(fēng)速功率特性曲線和風(fēng)速風(fēng)能利用系數(shù)曲線，如圖2和圖3。

圖2 風(fēng)速功率特性曲線

圖3 風(fēng)速功率系數(shù)曲線

從圖2可以看出，1#風(fēng)力機(jī)在額定風(fēng)速12 m/s的輸出功率是405.3W，2#風(fēng)力機(jī)在額定風(fēng)速10m/s的輸出功率是303W，都達(dá)到設(shè)計(jì)工況要求。2#風(fēng)力機(jī)在風(fēng)速3～14 m/s的范圍內(nèi)輸出功率明顯大于1#風(fēng)力機(jī)輸出功率。

從圖3可以看出，1#風(fēng)力機(jī)和2#風(fēng)力機(jī)的功率系數(shù)最大值都出現(xiàn)在5 m/s，分別為0.248和0.395。2#風(fēng)力機(jī)的功率系數(shù)和1#風(fēng)力機(jī)功率系數(shù)風(fēng)速3～14 m/s的范圍變化趨勢(shì)相同，2#風(fēng)力機(jī)的功率系數(shù)明顯大于1#風(fēng)力機(jī)功率系數(shù)。

以上分析說(shuō)明2#風(fēng)力機(jī)所采用的翼型能夠滿足風(fēng)力機(jī)的出力，同時(shí)提高了風(fēng)力機(jī)對(duì)風(fēng)能的利用率，達(dá)到了既追求輸出功率的增大，又追求在各風(fēng)速段較高的功率系數(shù)的設(shè)計(jì)要求。

2 風(fēng)譜圖建立和分析

風(fēng)譜圖是指某個(gè)區(qū)域獨(dú)立于地形地貌障礙物等因素之外的風(fēng)況特征。根據(jù)中國(guó)氣象科學(xué)數(shù)據(jù)共享網(wǎng)提供的2008年內(nèi)蒙古24個(gè)地區(qū)的風(fēng)資源數(shù)據(jù)，利用WAsP軟件進(jìn)行分析統(tǒng)計(jì)，選出具有代表特性的4個(gè)地區(qū)，分別是阿巴嘎旗、巴彥毛道、朱日和、海力素，對(duì)以上4個(gè)地區(qū)進(jìn)行分析，得到風(fēng)譜圖的相關(guān)參數(shù)，為計(jì)算風(fēng)力機(jī)年發(fā)電量提供數(shù)據(jù)。

利用WAsP OWC Wizard工具，首先導(dǎo)入欲安裝風(fēng)力機(jī)地區(qū)1年的原始的按時(shí)間序列的風(fēng)速風(fēng)向數(shù)據(jù)；其次輸入取得風(fēng)資源數(shù)據(jù)的地點(diǎn)及設(shè)備距地面高度，主要是指實(shí)測(cè)點(diǎn)的經(jīng)度、緯度、當(dāng)?shù)貥?biāo)準(zhǔn)氣壓、溫度及海拔高度；然后設(shè)置單位換算系數(shù)和風(fēng)速的上下限以及16風(fēng)向方位；最后計(jì)算得到風(fēng)譜圖[9]。按照以上步驟分析得到內(nèi)蒙古阿巴嘎旗、巴彥毛道、朱日和、海力素4個(gè)地區(qū)2008年的風(fēng)譜圖，見(jiàn)圖4～圖7。從圖4～圖7可以看出4個(gè)地區(qū)的2008年平均風(fēng)速分別是阿巴嘎旗3.5 m/s，巴彥毛道3.98 m/s，朱日和4.19m/s，海力素4.94m/s。

3 風(fēng)機(jī)站的建立

圖4 2008年阿巴嘎旗地區(qū)風(fēng)譜圖

圖6 2008年朱日和地區(qū)風(fēng)譜圖

圖7 2008年海力素地區(qū)風(fēng)譜圖

首先，利用地圖編輯WAsP Map Editor工具在數(shù)字化圖板或地圖掃描圖形文件上，用鼠標(biāo)選取需要的地圖直接生成WAsP可以接受地圖文件；然后提供風(fēng)力機(jī)站的位置坐標(biāo)：阿巴嘎旗地區(qū)是北緯44.01°，東經(jīng) 114.57°，海拔高度 112.61 m；巴彥毛道地區(qū)是北緯40.10°，東經(jīng)104.48°，海拔高度132.39m；朱日和地區(qū)是北緯42.24°，東經(jīng)112.54°，海拔高度115.08 m；海力素地區(qū)北緯 41.24°，東經(jīng) 106.24°，海拔高度150.96m，按照給定坐標(biāo)在地圖上標(biāo)定風(fēng)力機(jī)站的位置；最后在風(fēng)力機(jī)站周?chē)鶕?jù)實(shí)際情況選擇設(shè)置障礙物參數(shù)和地面粗糙度參數(shù)，由于風(fēng)力機(jī)安裝內(nèi)蒙古4個(gè)地區(qū)的戈壁、草原等非常開(kāi)闊平坦的區(qū)域，阻風(fēng)物很少，障礙物參數(shù)的取值從最小的房屋建筑為0到最大的防風(fēng)樹(shù)林為0.5，故該文選取障礙物參數(shù)房屋建筑0，粗糙度1級(jí)（Z=0.003）[10]。

4 風(fēng)力機(jī)發(fā)電量計(jì)算及分析

利用2008年4個(gè)地區(qū)的年風(fēng)譜圖和兩種風(fēng)力機(jī)的功率特性曲線，通過(guò)計(jì)算得到兩種風(fēng)力機(jī)在4個(gè)地區(qū)的年發(fā)電量和單位面積年發(fā)電量，見(jiàn)圖8和圖9。

從圖8和圖9可以看出，2#風(fēng)力機(jī)無(wú)論在年平均風(fēng)速較小的阿巴嘎旗風(fēng)場(chǎng)，還是在年平均風(fēng)速較大的海力素風(fēng)場(chǎng)的單位面積年發(fā)電量明顯大于1#風(fēng)力機(jī)單位面積年發(fā)電量，說(shuō)明2#風(fēng)力機(jī)采用新翼型的優(yōu)化設(shè)計(jì)提高了風(fēng)力機(jī)的輸出功率，也就是在實(shí)際風(fēng)場(chǎng)中提高了風(fēng)力機(jī)的單位面積的年發(fā)電量。通過(guò)計(jì)算分析，2#風(fēng)力機(jī)單位面積的發(fā)電量比1#風(fēng)力機(jī)單位面積年發(fā)電量在4個(gè)平均風(fēng)速由低到高的風(fēng)場(chǎng)分別提高了30.6%、29.1%、27.7%、26.5%，平均提高28.4%，說(shuō)明2#風(fēng)力機(jī)的性能明顯優(yōu)于1#風(fēng)力機(jī)[11-12]。

圖8 風(fēng)力機(jī)年發(fā)電量

圖9 風(fēng)力機(jī)單位面積年發(fā)電量

5 結(jié)束語(yǔ)

該文對(duì)采用兩種不同翼型的水平軸風(fēng)力機(jī)利用風(fēng)洞進(jìn)行了氣動(dòng)特性實(shí)驗(yàn)，得到了2個(gè)風(fēng)力機(jī)的風(fēng)速功率特性曲線和風(fēng)速風(fēng)能利用系數(shù)曲線，2#風(fēng)力機(jī)所采用的翼型能夠滿足風(fēng)力機(jī)的出力，同時(shí)提高了風(fēng)力機(jī)對(duì)風(fēng)能的利用率，達(dá)到了既追求輸出功率的增大，又追求在各風(fēng)速段較高的功率系數(shù)的設(shè)計(jì)要求。

利用WAsP軟件對(duì)內(nèi)蒙古阿巴嘎旗，巴彥毛道、朱日和、海力素4個(gè)地區(qū)2008年風(fēng)能資源進(jìn)行了評(píng)估，得到了風(fēng)譜圖。通過(guò)分析計(jì)算得到1#風(fēng)力機(jī)和2#風(fēng)力機(jī)的年發(fā)電量、單位面積年發(fā)電量，2#風(fēng)力機(jī)單位面積的發(fā)電量比1#風(fēng)力機(jī)單位面積年發(fā)電量平均提高28.4%，2#風(fēng)力機(jī)的性能明顯優(yōu)于1#風(fēng)力機(jī)，說(shuō)明新翼型在氣動(dòng)性能方面是優(yōu)于傳統(tǒng)翼型的。

但是對(duì)于采用新翼型的風(fēng)力機(jī)的可靠性能，需要對(duì)其結(jié)構(gòu)動(dòng)態(tài)進(jìn)行測(cè)試和分析，為設(shè)計(jì)輸出功率高、結(jié)構(gòu)穩(wěn)定的風(fēng)力機(jī)提供一定參考。

[1]王軍，周丙超.基于Matlab的小型風(fēng)力機(jī)葉片設(shè)計(jì)[J].水電能源科學(xué)，2007，25（5）：142-144.

[2]竇秀榮.水平軸風(fēng)力機(jī)起動(dòng)性能及結(jié)構(gòu)動(dòng)力學(xué)特性研究[D].濟(jì)南：山東工業(yè)大學(xué)，1995.

[3]巫發(fā)明，王立鵬，楊從新，等.基于遺傳算法的直葉片垂直軸風(fēng)力機(jī)風(fēng)輪優(yōu)化設(shè)計(jì)[J].西華大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，28（5）：36-39.

[4]George S.Dulikravich.Aerodynamic shape design and optimization：status and trends[J].Journal of Aircraft，1992，29（6）：1020-1026.

[5]Jamshid A.Generation Samareh.Status and future of geometry modeling and grid for design[J].Journal of Aircraft，1999，36（1）：97-104.

[6]Dai Y J，Wang J W，Wen C F.Determination of characteristic wind speed for wind power under conditions of low cost and constant power[C]∥Proceedings of the first IEEE International Conference on Mechanic Automation and Control Engineering，2010：4054-4057.

[7]GB/T 19068.3—2003離網(wǎng)型風(fēng)力發(fā)電機(jī)組第3部分：風(fēng)洞實(shí)驗(yàn)方法[S].北京：中國(guó)標(biāo)準(zhǔn)出版社，2003.

[8]GB/T 190682.2—2003離網(wǎng)型風(fēng)力發(fā)電機(jī)組第2部分：試驗(yàn)方法[S].北京：中國(guó)標(biāo)準(zhǔn)出版社，2003.

[9]Justus C G，Hargraves W R，Mikhail，et al.Methods for Estimating Wind Speed Frequency Distributions[J].Journal of Applied Meteorology，1978，17（3）：350-353.

[10]代元軍，汪建文，溫彩鳳，等.利用WasP軟件對(duì)風(fēng)力機(jī)發(fā)電量的預(yù)測(cè)[J].可再生能源，2011，29（1）：103-106.

[11]Alexiadis M C，Dokopoulos P S，Sahsamanoglou H S，et al.Wind speed and powerforecasting based on spatial correlation models[C]∥Institute of Electrical and Electronics Engineers Transactions on Energy Conversion，1999，14（3）：836-842.

[12]Jaramillo O A， Borja，BimOda M A.Weibull wind speed distributions：an analysis of wind energy potential in laventa，mexico[J].Wind Engineering，2004（28）：225-234.