姜廣緒，潘晶雯，田景奎

(1.中國電力工程顧問集團(tuán)華北電力設(shè)計院有限公司,北京市 100120；2.中國人民解放軍裝甲兵工程學(xué)院，北京市 100072)

?

雙參數(shù)威布爾分布風(fēng)況中基于k值分析的能量分布研究

姜廣緒1，潘晶雯2，田景奎1

(1.中國電力工程顧問集團(tuán)華北電力設(shè)計院有限公司,北京市 100120；2.中國人民解放軍裝甲兵工程學(xué)院，北京市 100072)

在風(fēng)力發(fā)電場風(fēng)能資源分析中，僅采用年平均風(fēng)速和年平均風(fēng)功率密度來判斷風(fēng)電場風(fēng)能資源的好壞是不全面的，其無法體現(xiàn)風(fēng)電場風(fēng)能資源的能量輸出能力。通過對單峰雙參數(shù)威布爾分布風(fēng)況中形狀參數(shù)k值的研究，找出其與風(fēng)速、風(fēng)能頻率分布及風(fēng)中蘊(yùn)含能量的關(guān)系，并進(jìn)一步分析對于不同風(fēng)力發(fā)電機(jī)組Cp值曲線，k值與總體有效風(fēng)能利用率的關(guān)系，從而在風(fēng)電場建設(shè)的咨詢、決策、設(shè)計等階段對風(fēng)電場風(fēng)能資源評估和風(fēng)機(jī)預(yù)選型作出科學(xué)的指導(dǎo)。

威布爾分布；風(fēng)力發(fā)電；能量分布；形狀參數(shù)

0 引 言

風(fēng)力發(fā)電場是將風(fēng)能轉(zhuǎn)化為機(jī)械能，再將機(jī)械能轉(zhuǎn)化成電能的裝置總成在一起的系統(tǒng)[1]。在風(fēng)電場開發(fā)設(shè)計過程中，對場區(qū)內(nèi)風(fēng)能資源的評估研究是前期工作最重要的環(huán)節(jié)，這將從根本上決定風(fēng)電場的優(yōu)劣[2-3]。

實際工程中經(jīng)常用年平均風(fēng)速和年平均風(fēng)功率密度來直觀地評價某地的風(fēng)能資源情況，這種直觀的評價為風(fēng)電場選址提供了重要依據(jù)[4-5]。但在實際工程中，僅利用年平均風(fēng)速和年平均風(fēng)功率密度評價風(fēng)能資源有一定的偏差甚至評價結(jié)果與實際相反。例如，山西省北部和中部A、B這2個風(fēng)場測風(fēng)塔處80 m高年平均風(fēng)速均為6.9 m/s，但A風(fēng)場測風(fēng)塔80 m高年平均風(fēng)功率密度為393 W/m2，B風(fēng)場的為278 W/m2。直觀來看，顯然是A風(fēng)場較好，但對這2個測風(fēng)塔處選用同種機(jī)型進(jìn)行發(fā)電量計算發(fā)現(xiàn)，B的發(fā)電量要遠(yuǎn)高于A的發(fā)電量。

隨著我國風(fēng)電場的大規(guī)模建設(shè)，如何準(zhǔn)確評估風(fēng)能資源成為研究熱點(diǎn)。文獻(xiàn)[6]指出風(fēng)功率密度等級的劃分僅用風(fēng)功率密度、年平均風(fēng)速來要求存在一定的缺陷，并建議根據(jù)風(fēng)場測風(fēng)數(shù)據(jù)統(tǒng)計出的A、k值草擬《風(fēng)功率密度等級和風(fēng)能區(qū)等級劃分標(biāo)準(zhǔn)》確定等級。 一些學(xué)者從風(fēng)速頻率分布的參數(shù)(如威布爾分布的k值)影響入手，研究了其對風(fēng)速頻率分布自身的影響[7]，這從理論上對風(fēng)況自身的好壞評價可作出較好的分析。而實際上風(fēng)力發(fā)電場畢竟是要通過能量轉(zhuǎn)化裝置即風(fēng)力發(fā)電機(jī)組來將風(fēng)能轉(zhuǎn)化為電能，因此，上述研究很難結(jié)合風(fēng)力發(fā)電機(jī)組的性能對風(fēng)能資源的能量輸出能力進(jìn)行分析。文獻(xiàn)[8]針對上述情況提出了相當(dāng)風(fēng)速和有功風(fēng)功率密度這2個概念，結(jié)合風(fēng)力發(fā)電機(jī)組的性能對風(fēng)能資源情況進(jìn)行評估。但上述研究僅限定在理論概念上，沒有結(jié)合風(fēng)速頻率分布及分布參數(shù)的影響進(jìn)行分析。

本文通過對雙參數(shù)威布爾分布風(fēng)況中形狀參數(shù)k值的研究，找出其與風(fēng)速頻率分布、風(fēng)能頻率分布以及風(fēng)電場能量輸出能力之間的關(guān)系和相應(yīng)規(guī)律，為風(fēng)電場風(fēng)能資源的評估提供科學(xué)依據(jù)。

1 理論基礎(chǔ)

雙參數(shù)威布爾分布被普遍認(rèn)為是適用于描述風(fēng)速分布統(tǒng)計的概率密度函數(shù)，其是一種單峰的、雙參數(shù)分布函數(shù)簇[9]。風(fēng)速概率密度函數(shù)可表達(dá)為

(1)

式中：k為形狀參數(shù)，無量綱；c為尺度參數(shù)，m/s；v為風(fēng)速，m/s。

風(fēng)能公式為

(2)

式中ρ為空氣密度，kg/m3。

由此可以推導(dǎo)平均風(fēng)功率密度的數(shù)學(xué)期望E(w)[10-11]，假定空氣密度一定,則

(3)

(4)

式中Г為伽馬函數(shù)。

同理可以推導(dǎo)平均風(fēng)速的數(shù)學(xué)期望：

(5)

即

(6)

將式(6)代入式(4)中可以得到：

(7)

(8)

2 k值與風(fēng)中蘊(yùn)含能量的關(guān)系

由公式(7)可知，某一符合雙參數(shù)威布爾分布的風(fēng)況中，其平均風(fēng)功率密度與空氣密度、平均風(fēng)速和形狀參數(shù)(即k值)有關(guān)。假定空氣密度為定值不變，當(dāng)平均風(fēng)速一定時，平均風(fēng)功率密度只與k值有關(guān)。

根據(jù)工程經(jīng)驗，中國的風(fēng)電場k值一般為1.2～3[12]，本次選取1～5作為研究區(qū)段，以0.1作為步長對θ進(jìn)行定量計算，結(jié)果如圖1所示。

由圖1可知，θ值隨著k值的增大逐漸減小，當(dāng)k趨近于∞時，θ值趨近于1，即對于服從雙參數(shù)威布爾分布的風(fēng)況，平均風(fēng)速一定時，k值越小，風(fēng)中蘊(yùn)含的能量越多。

3 k值與風(fēng)速、風(fēng)能頻率分布的關(guān)系

以平均風(fēng)速7 m/s的風(fēng)況為例，分別對k=1.5，2.0，2.5和3的風(fēng)速頻率分布情況進(jìn)行分析，如圖2所示。

圖2 風(fēng)速頻率分布隨k值變化情況Fig.2 Wind speed frequency distribution with k value changes

由圖2可知，在研究區(qū)域內(nèi)，隨著k值的增大，分布的峰值逐漸增大，且出現(xiàn)的風(fēng)速區(qū)域越來越大并接近平均風(fēng)速，同時小風(fēng)速區(qū)域(0～5 m/s)頻率逐漸減小，平均風(fēng)速附近區(qū)域(5～11 m/s)頻率逐漸增大，大風(fēng)速區(qū)域(11 m/s以上)頻率逐漸減小。也就是說k值越大，分布越集中，k值越小，分布越分散。這也可以從k值的估算公式上加以驗證：

(9)

式中s為分布的標(biāo)準(zhǔn)差[13]。

可見k值越小，標(biāo)準(zhǔn)差s越大，說明數(shù)據(jù)分布越分散。

對上述風(fēng)速頻率分布分別作出相應(yīng)的風(fēng)能頻率分布曲線，如圖3所示。

圖3 風(fēng)能頻率分布隨k值變化情況Fig.3 Wind energy frequency distribution with k value changes

從圖3中可以看出，在研究區(qū)域內(nèi)，隨著k值的增大，分布的峰值逐漸增大且出現(xiàn)的風(fēng)速區(qū)域越來越小并接近平均風(fēng)速，同時小風(fēng)速區(qū)域頻率變化很小，平均風(fēng)速附近區(qū)域頻率逐漸增大，大風(fēng)速區(qū)域頻率逐漸減小。

由上述分析可知，在威布爾單峰雙參數(shù)風(fēng)況分布中，隨著k值的減小，風(fēng)中蘊(yùn)含的能量逐漸增加，但同時能量的分布也越來越分散。

4 k值與總體有效風(fēng)能利用率的關(guān)系

風(fēng)力發(fā)電機(jī)組在變轉(zhuǎn)速區(qū)(額定風(fēng)速以下)需要通過變槳對轉(zhuǎn)速進(jìn)行相應(yīng)控制，以達(dá)到較好的葉尖速比，從而提高相應(yīng)的風(fēng)能利用系數(shù)，即Cp值[14]。而在定轉(zhuǎn)速區(qū)(額定風(fēng)速以上)，為了穩(wěn)定輸出和安全性問題，又需要通過變槳使葉輪轉(zhuǎn)速恒定，這就導(dǎo)致了葉尖速比隨著來流風(fēng)速的增加顯著降低，導(dǎo)致了大風(fēng)速區(qū)域風(fēng)能利用系數(shù)較低[15]。圖4給出了某風(fēng)力發(fā)電機(jī)組的設(shè)計Cp曲線。

圖4 某風(fēng)力發(fā)電機(jī)組的Cp值曲線Fig.4 Cp curve of a kind of wind turbine

表1k值與總體有效風(fēng)能利用率的對應(yīng)表

Table 1 Correspondence table ofkand overall effective wind energy utilization

結(jié)合前述分析，當(dāng)平均風(fēng)速一定時，隨著k值的減小，風(fēng)力發(fā)電機(jī)組風(fēng)能利用系數(shù)較大的小風(fēng)速區(qū)域風(fēng)能頻率明顯降低，而在風(fēng)能頻率顯著增大的大風(fēng)速區(qū)域，風(fēng)能利用系數(shù)又很小，這就導(dǎo)致了表中μCp減小。因此隨著k值的減小，雖然風(fēng)中蘊(yùn)含的風(fēng)能逐漸增加(θ值變大)，但平均有效利用風(fēng)功率密度卻不是一直在增加，而是在達(dá)到某個最優(yōu)值后逐漸減小。因此會出現(xiàn)相同的年平均風(fēng)速，年平均風(fēng)功率密度較高的風(fēng)況下發(fā)電量反而較低的情況。

通過分析，發(fā)現(xiàn)針對特定的風(fēng)力發(fā)電機(jī)組和平均風(fēng)速，平均有效利用風(fēng)功率密度存在一個最優(yōu)值，相應(yīng)也就存在一個最優(yōu)k值。例如采用文中提到的機(jī)型，當(dāng)平均風(fēng)速為7 m/s時，對應(yīng)的最優(yōu)k值約為2.5；當(dāng)平均風(fēng)速為6.5 m/s時，對應(yīng)的最優(yōu)k值約為2.0；當(dāng)平均風(fēng)速為7.5 m/s時，對應(yīng)的最優(yōu)k值約為3.4。一般來說，平均風(fēng)速越大，對應(yīng)的最優(yōu)k值也越大。在風(fēng)力發(fā)電機(jī)組選型時，機(jī)組在實際平均風(fēng)速下對應(yīng)的最優(yōu)k值越接近風(fēng)場風(fēng)況的實際k值，發(fā)電量也就越高。

進(jìn)一步分析在低風(fēng)速情況下，如平均風(fēng)速為5.5 m/s時，其最優(yōu)k值小于1.3，當(dāng)k值為2.0時，其電量折合滿發(fā)小時數(shù)約為2 202 h，如果有20%的發(fā)電量折減，則折合滿發(fā)小時數(shù)僅為1 762 h；如平均風(fēng)速為5 m/s時，其最優(yōu)k值小于1.1，當(dāng)k值為2.0時，其電量折合滿發(fā)小時數(shù)約為1 765 h，如果有20%的發(fā)電量折減，則折合滿發(fā)小時數(shù)僅為1 412 h。在服從雙參數(shù)單峰威布爾分布風(fēng)況下，平均風(fēng)速小于5.5 m/s，k值又較大時，該處可能不適宜布置本文所選的風(fēng)力發(fā)電機(jī)組，或應(yīng)在滿足安全性的前提下，選擇更大葉輪直徑機(jī)組，并在低風(fēng)速區(qū)域盡量地提高風(fēng)能利用系數(shù)，然而這種提高是有限的[16]。

5 結(jié) 論

(1)在研究范圍內(nèi)，平均風(fēng)速一定的情況下，平均有效利用風(fēng)功率密度系θ3μCp對于k值是一個減函數(shù)與增函數(shù)的乘積，其存在一個最大值。這個最大值對應(yīng)的k值隨著平均風(fēng)速的變化而不同。

(2)平均風(fēng)速較低的風(fēng)場在資源評估時容易出現(xiàn)決策失誤，故建議低風(fēng)速風(fēng)場的區(qū)分與決策不應(yīng)簡單通過平均風(fēng)速或平均風(fēng)功率密度來決定，而應(yīng)同時結(jié)合k值來分析。

(3)在風(fēng)電場建設(shè)的咨詢、決策、設(shè)計等階段，對擬建風(fēng)電場內(nèi)風(fēng)能資源的平均能量風(fēng)速系數(shù)θ和平均能量風(fēng)速v能進(jìn)行研究分析，可更準(zhǔn)確地指導(dǎo)風(fēng)電場的宏觀選址、微觀選址以及風(fēng)力發(fā)電機(jī)組選型等工作。因此建議風(fēng)力發(fā)電機(jī)組廠家在機(jī)組技術(shù)參數(shù)文件中增加針對該機(jī)組的k值與總體有效風(fēng)能利用率的對應(yīng)表。

[1]姚興佳，王士榮，董麗萍.風(fēng)力發(fā)電技術(shù)講座(一)風(fēng)力發(fā)電機(jī)組(1)[J].可再生能源，2006，127(3)：102-104.Yao Xingjia,Wang Shirong,Dong Liping.Lessons of status of wind power technology(Ⅲ) Wind power generator(1)[J].Renewable Energy,2006(3):102-104.

[2]張洋，劉雪楓，李強(qiáng).風(fēng)資源評估綜述[J].水利電力機(jī)械，2007，29(5)：72-75.Zhang Yang, LiuXuefeng, LiQiang.Evaluation outline of wind resources[J].Water Conservancy & Electric Power Machinery,2007,29(5):72-75 .

[3]華澤嘉，閆少軍，楊振華.風(fēng)資源評估方法的研究[J].技術(shù)與市場，2012，19(11)：33-35.

[4]GB/T 18709—2002 風(fēng)電場風(fēng)能資源測量方法[S].

[5]GB/T 18710—2002風(fēng)電場風(fēng)能資源評估方法[S].

[6]王曉林.風(fēng)功率密度等級劃分和風(fēng)能區(qū)劃的修訂建議[J].電力技術(shù)，2010，19(8)：78-82.Wang Xiaolin.Proposals of amendments of wind power density grade and wind energy area classification[J].Electric Power Technology, 2010, 19(8):78-82.

[7]朱德臣，汪建文.風(fēng)工況雙參數(shù)威布爾分布k值影響研究[J].太陽能，2007(6)：34-36.Zhu Dechen, Wang Jianwen.Study on the effects of k for two-parameter Weibull Distribution in wind resource assessment[J].Solar Energy,2007(6):34-36.

[8]楊振斌，朱瑞兆，薛桁.風(fēng)電場風(fēng)能資源評價兩個新參數(shù)：相當(dāng)風(fēng)速、有功風(fēng)功率密度[J].太陽能學(xué)報，2007，28(3)：248-251.Yang Zhenbin,Zhu Ruizhao,XueHeng.Two New Concepts on wind energy assessment in wind farm: equivalent wind speed, available wind power density[J].ActaEnergiae Solaris Sinica,2007,28(3):248-251.

[9]王承煦，張源.風(fēng)力發(fā)電[M].北京：中國電力出版社，2002：31-47.

[10]宮靖遠(yuǎn)，賀德馨，孫如林，等.風(fēng)電場工程技術(shù)手冊[M].北京：機(jī)械工業(yè)出版社，2004：19-24.

[11]伯頓 T，等.風(fēng)能技術(shù)[M].武鑫，等，譯.北京：科學(xué)出版社，2007：12-34.

[12]李軍，胡非，王丙蘭，等.風(fēng)速的Weibull分布參數(shù)[J].太陽能學(xué)報，2012，33(10)：1667-1671.Li Jun, Hu Fei, Wang Binglan, et al.Weibull parameters for wind speed[J].ActaEnergiae Solaris Sinica, 2012, 33(10):1667-1671.

[13]徐衛(wèi)民，張星琳，孔新紅.三種計算風(fēng)速威布爾分布參數(shù)的比較[J].江西電力，2007,31(1):1-3.XuWeiming,Zhang Xinglin,Kong Xinhong.Comparison for three methods to calculate Weibull distribution parameters of wind speed[J].Jiangxi Electric Power,2007,31(1):1-3.

[14]胡燕平,甄海華,戴巨川.變槳距風(fēng)力發(fā)電機(jī)額定風(fēng)速的確定方法[J].太陽能學(xué)報，2011，32(3)：307-310.Hu Yanping,Zhen Haihua,Dai Juchuan.The study on the determining method for the nominal wind speed of pitch wind turbies[J].ActaEnergiae Solaris Sinica, 2011,32(3):307-310.

[15]高俊瑩，徐建軍，許愛華，等.風(fēng)能利用系數(shù)模糊控制策略的研究[J].電氣技術(shù)，2009(8)：36-42.GaoJunying,XuJianjun,XuAihui, et al.Study on fuzzy control strategy of wind energy utilization factor[J].Electrical Engineering,2009(8):36-42.

[16]鄧英，張丹，周峰，等.提高風(fēng)電機(jī)組風(fēng)能轉(zhuǎn)化效率的變速控制技術(shù)研究[J].電力與能源進(jìn)展，2013，1(2)：64-71.Deng Ying, Zhang Dan, Zhou Feng, et al.Study to improve the wind energy conversion efficiency with the variable speed control[J].Advances in Energy and Power Engineering, 2013, 1(2):64-71.

(編輯：張小飛)

Energy Distribution Research Based on the ParameterkAnalysis of Two-Parameter Weibull Distribution Wind Conditions

JIANG Guangxu1, PAN Jingwen2, TIAN Jingkui1

(1.North China Power Engineering Co., Ltd.of China Power Engineering Consulting Group, Beijing 100120, China;2.Academy of Armored Forces Engineering of PLA, Beijing 100072, China)

Only by using the annual average wind speed and wind power density to determine the quality of wind energy resources is not comprehensive enough in wind resource analysis, which can not reflect the output capability of the wind resource.Investigating thekvalue，which is the shape parameter of two-parameter Weibull distribution wind conditions,is to find out its inherent relationship with the wind energy, wind speed distribution and wind energy distribution.And further research is put forward to find out the relationship between thekvalue and the effective wind energy utilization for differentCpcurves of wind turbine generators.The above research will supply a scientific guidance for wind resource assessment and wind turbine generator pre-selection during the consultation, decision-making, design and other stages of wind farm construction.

Weibull distribution; wind power; energy distribution; shape parameter

TM 614

A

1000-7229(2015)03-0105-04

10.3969/j.issn.1000-7229.2015.03.018

2014-09-04

2014-09-24

姜廣緒(1983)，男，碩士，工程師，主要從事風(fēng)力發(fā)電、太陽能光伏發(fā)電設(shè)計方面的研究工作；

潘晶雯(1984)，女，碩士，講師，主要從事工程力學(xué)、材料力學(xué)方面的教學(xué)和研究工作；

田景奎(1960)，男，工學(xué)學(xué)士，教授級高級工程師，主要從事風(fēng)力發(fā)電、太陽能光伏發(fā)電設(shè)計方面的研究和管理工作。