張國強，田宏哲，楊政厚

（北京華能新銳控制技術(shù)有限公司，北京 102209）

0 引言

風(fēng)能作為一種綠色清潔能源，風(fēng)力發(fā)電發(fā)展迅速。隨著技術(shù)進(jìn)步，對風(fēng)電機組控制系統(tǒng)要求越來越高，控制趨向多元化、精細(xì)化、智能化方向。為了提高機組發(fā)電效率，提升在役風(fēng)電場發(fā)電量，最大風(fēng)能捕獲控制算法成為了研究熱點[1]～[3]。

風(fēng)電機組發(fā)電功率直接受到風(fēng)速和環(huán)境參數(shù)的影響，許多專家學(xué)者針對此問題提出了解決方法。文獻(xiàn)[4]考慮環(huán)境溫度和濕度，提出了功率反饋的最大風(fēng)能捕獲控制PSF 算法，得到較好的仿真測試結(jié)果。文獻(xiàn)[5]提出了在空氣密度降低情況下的功率自適應(yīng)的優(yōu)化方法，并在MATLAB 平臺上搭建直驅(qū)機組模型進(jìn)行仿真，優(yōu)化方案在風(fēng)速階躍階段效果較好。但是僅做了恒定風(fēng)速的仿真，未考慮機組的非線性氣動特性。文獻(xiàn)[6]分析了影響空氣密度的環(huán)境因素，提出了自適應(yīng)空氣密度變化的控制算法，并結(jié)合實際機組進(jìn)行驗證。但測試過程中未考慮兩臺機組自身的性能差異導(dǎo)致的發(fā)電功率不同。

本文從風(fēng)電機組基本發(fā)電原理出發(fā)，分析額定風(fēng)速以下最優(yōu)轉(zhuǎn)矩控制與空氣密度的關(guān)系，提出了空氣密度變化修正的最優(yōu)轉(zhuǎn)矩控制方法。結(jié)合吉林省某風(fēng)電場1年的實際測量數(shù)據(jù)，分析了風(fēng)電場環(huán)境下空氣密度變化規(guī)律。最后，基于Bladed 軟件平臺，在平均風(fēng)速6 m/s 湍流風(fēng)況下，進(jìn)行了5 組仿真試驗，驗證了所提出的修正算法的功率提升效果。

1 空氣密度與風(fēng)電機組轉(zhuǎn)矩控制的關(guān)系

在實際運行中，風(fēng)電機組輸出功率的大小直接受到風(fēng)的影響。通常只關(guān)注風(fēng)速和風(fēng)向，忽略了空氣密度變化的影響，本文從理論角度，結(jié)合兆瓦級雙饋型風(fēng)電機組，分析了空氣密度對風(fēng)電機組功率的影響。

風(fēng)電機組的輸出功率與風(fēng)速的關(guān)系為

式中：Cp為功率系數(shù)，表征風(fēng)電機組的風(fēng)能捕獲效率；ρ 為空氣密度；A 為風(fēng)輪的掃略面積 A=πR2；V 為來流風(fēng)速。

風(fēng)電機組的輸出功率與空氣密度大小直接相關(guān)，圖1 為同一臺風(fēng)電機組在不同的空氣密度下的風(fēng)功率曲線。

風(fēng)電廠商一般提供標(biāo)準(zhǔn)空氣密度（1.225 kg/m3）下的功率曲線，而風(fēng)電場的實際空氣密度并不是標(biāo)準(zhǔn)條件下的密度，因此，在進(jìn)行運行數(shù)據(jù)擬合功率曲線時須要依據(jù)現(xiàn)場的空氣密度進(jìn)行修正[7]。

圖1 不同空氣密度下的風(fēng)功率曲線Fig.1 Wind-power curve in different air density

式中：ωr為風(fēng)輪轉(zhuǎn)速。

風(fēng)電機組 Cp是 V，ωr和槳距角 β 3 個變量的函數(shù)。在額定風(fēng)速以下，通常β 保持在最佳角度不改變，此時Cp和λ 的關(guān)系曲線如圖2 所示。由圖2 可知，在額定風(fēng)速以下功率系數(shù)為λ 的單值凸函數(shù)，即在運行區(qū)間內(nèi)存在一個最大的Cp值，此時，對應(yīng)葉尖速比λopt稱為最優(yōu)λ。最大風(fēng)能捕獲原理就是通過控制發(fā)電機轉(zhuǎn)速，在不同的風(fēng)速下使 λ 保持在 λopt，從而保證 Cp最大。

圖2 β=0 時功率系數(shù)與葉尖速比的關(guān)系Fig.2 Relationship between power factor and tip speed ratio when β=0

風(fēng)電機組發(fā)電系統(tǒng)屬于旋轉(zhuǎn)機械，其輸出功率也可以表示為轉(zhuǎn)速和轉(zhuǎn)矩的乘積，忽略機械損失和轉(zhuǎn)化效率，可得：

式中：Tr為風(fēng)輪系統(tǒng)的轉(zhuǎn)矩；Tg，ωg分別為發(fā)電機的轉(zhuǎn)矩和轉(zhuǎn)速值，其中 ωg=Gωr，G 為齒輪箱的傳動比。

聯(lián)立式（2）和式（3），得到：

式中：Kopc為控制參數(shù)。

在額定風(fēng)速以下，雙饋發(fā)電機轉(zhuǎn)速控制是通過調(diào)節(jié)電機電磁轉(zhuǎn)矩實現(xiàn)的，因此，只要控制發(fā)電機轉(zhuǎn)矩和轉(zhuǎn)速使其滿足式（4）的關(guān)系，就可以保證風(fēng)電機組獲得最大的風(fēng)能捕獲效率。Kopc與ρ 呈正相關(guān)關(guān)系，現(xiàn)場的空氣密度實時變化，為保證獲得最大功率，須要根據(jù)實際空氣密度對控制系統(tǒng)進(jìn)行優(yōu)化。

2 現(xiàn)場數(shù)據(jù)分析

根據(jù)式（1）可知，空氣密度變化對風(fēng)電機組輸出功率有很大的影響。本文選取吉林省某風(fēng)電場2017-12-01-2018-11-30日測風(fēng)塔處的空氣密度數(shù)據(jù)，分析了風(fēng)電場環(huán)境下空氣密度的變化規(guī)律。數(shù)據(jù)間隔15 min，得到的原始數(shù)據(jù)如圖3 所示。

圖3 空氣密度原始數(shù)據(jù)曲線Fig.3 Curve of air density raw data

由于原始數(shù)據(jù)存在很多錯誤點，須對原始數(shù)據(jù)進(jìn)行篩選修正及分析統(tǒng)計。

（1）數(shù)據(jù)篩選修正

風(fēng)電場現(xiàn)場環(huán)境惡劣，由于傳感器故障、數(shù)據(jù)通訊存儲錯誤等原因造成采集的數(shù)據(jù)不真實。在進(jìn)行數(shù)據(jù)分析前，須要進(jìn)行篩選及修正。

限值判斷：結(jié)合空氣密度實際物理意義，將圖3 中數(shù)值為0 和小于0.8 的數(shù)據(jù)判斷為錯誤點。

變化率判斷： 空氣密度為大氣環(huán)境的一項表征參數(shù)，是連續(xù)緩變的過程，若相鄰采樣點數(shù)據(jù)發(fā)生突變（本文設(shè)置為0.12），判斷為變化率超限錯誤點。

修正方法：對原始錯誤點進(jìn)行標(biāo)記。本文采用錯誤數(shù)據(jù)點前一天（96 點）平均值替換錯誤數(shù)據(jù)。

（2）數(shù)據(jù)處理

根據(jù)修正后的數(shù)據(jù)，選取其中無錯誤的3 d時間段的數(shù)據(jù)，得到圖4 所示的空氣密度變化曲線。分析空氣密度的每日變化規(guī)律，分別計算空氣密度的日平均值、月平均值和年平均值，并統(tǒng)計365 個日平均值的最大值和最小值，計算結(jié)果如圖5 所示。

圖4 空氣密度3 d 數(shù)據(jù)Fig.4 Air density data in 3 days

圖5 空氣密度日均值、月均值和年均值分布Fig.5 Distribution of air density about daily average，monthly mean and annual mean

（3）結(jié)果分析

由圖4 可知，空氣密度整體趨勢呈現(xiàn)以日為周期的變化規(guī)律，夜晚空氣密度大于白天，每個時刻的值在一定范圍內(nèi)進(jìn)行波動，每日的最大值和最小值約相差0.1 kg/m3。

由圖5 可知，12 個月空氣密度均值在年度區(qū)間內(nèi)呈季節(jié)性變化，冬季大于夏季，整體趨勢近似按照正弦規(guī)律變化。365 個日均值不是一條平滑的曲線，在一定范圍內(nèi)有較大波動。所研究風(fēng)電場空氣密度年平均值為1.253 3 kg/m3，日均值的最大值為1.425 4 kg/m3，日均值最小值為1.140 4 kg/m3，最大值與最小值相差25%，日均空氣密度在冬季和夏季變化很大。

3 考慮空氣密度變化的轉(zhuǎn)矩控制修正

由式（5）可知 Kopt與 ρ 成正比關(guān)系，實際機組控制系統(tǒng)中采用恒定值作為控制參考，未考慮現(xiàn)場空氣密度的季節(jié)性變化，考慮空氣密度修正后式（5）為

式中：ρ0和Kopt-0分別為標(biāo)準(zhǔn)空氣密度及其對應(yīng)的最優(yōu)轉(zhuǎn)矩系數(shù)；ρact和Kopt-act分別為現(xiàn)場實際空氣密度和修正后的最優(yōu)轉(zhuǎn)矩系數(shù)。

根據(jù)式（6），在傳統(tǒng)最優(yōu)轉(zhuǎn)矩PID 控制算法基礎(chǔ)上，增加空氣密度修正系數(shù)，控制流程如圖6所示。

圖6 帶有空氣密度修正的轉(zhuǎn)矩控制算法Fig.6 Torque control algorithm with air density correction

本文基于Bladed 軟件平臺，在1.5 MW 雙饋風(fēng)電機組模型基礎(chǔ)上，進(jìn)行了空氣密度修正轉(zhuǎn)矩控制算法仿真，設(shè)計了5 組對照實驗。其中：T1，T2 兩組空氣密度設(shè)置為風(fēng)電場日均最大值，T1組轉(zhuǎn)矩控制算法未修正，T2 組采用轉(zhuǎn)矩修正算法；T3，T4 兩組空氣密度設(shè)置為日均最小值，T3組轉(zhuǎn)矩控制算法未修正，T4 組采用轉(zhuǎn)矩修正算法；T5 組為對照組，空氣密度和轉(zhuǎn)矩系數(shù)皆采用標(biāo)準(zhǔn)值。

仿真具體參數(shù)見表1。風(fēng)速設(shè)置為平均值6 m/s 的Kalman 湍流風(fēng)模型，仿真時間為600 s，其他模型參數(shù)和控制參數(shù)均不作改變。

表1 仿真實驗組設(shè)置參數(shù)表Table 1 Parameter of simulation experiment group setting

4 仿真結(jié)果分析

對比5 組仿真結(jié)果，分別計算發(fā)電功率、發(fā)電機轉(zhuǎn)矩和發(fā)電機轉(zhuǎn)速的均值和標(biāo)準(zhǔn)差，計算結(jié)果如表2 所示。

根據(jù) IEC 標(biāo)準(zhǔn)的建議[8]，采用 Bins 法對功率和風(fēng)速數(shù)據(jù)進(jìn)行分析，在運行范圍內(nèi)將風(fēng)速分為間隔0.5 m/s 的區(qū)間，在每個區(qū)間內(nèi)分別統(tǒng)計風(fēng)速及功率的平均值，得到如圖7 所示風(fēng)功率曲線。

表2 5 組湍流風(fēng)況仿真結(jié)果Table 2 Result of five simulation experiment in turbulent wind

圖7 仿真結(jié)果的統(tǒng)計風(fēng)功率曲線Fig.7 Wind-power curve of the simulation result

由圖7 可知：

（1）對比實驗組 T1，T3 和 T5，空氣密度是直接影響機組輸出功率的因素，空氣密度越大，相同風(fēng)速下機組的發(fā)電功率越大；

（2）對比 T1 和 T2，空氣密度相差 16.36%，平均發(fā)電功率經(jīng)過空氣密度修正后降低了1.00%；對比T3 和T4，空氣密度相差6.91%，修正后平均功率提升了0.59%。在同一空氣密度下，經(jīng)過空氣密度修正轉(zhuǎn)矩控制算法，對機組功率輸出影響不大；

（3）對比 T1 和 T2 的結(jié)果，發(fā)電機轉(zhuǎn)矩和轉(zhuǎn)速變化分別為4.20%和5.56%；對比T3 和T4 的結(jié)果，發(fā)電機轉(zhuǎn)矩和轉(zhuǎn)速變化分別為1.73%和2.39%?？諝饷芏刃拚蟀l(fā)電機轉(zhuǎn)矩和轉(zhuǎn)速將會產(chǎn)生較大的變化，影響風(fēng)電機組載荷。

針對空氣密度修正后，風(fēng)電機組功率改變較小的結(jié)果，進(jìn)一步進(jìn)行仿真實驗研究。設(shè)計兩組平均風(fēng)速為6 m/s 恒定風(fēng)，空氣密度設(shè)定為1.425 4 kg/m3仿真實驗，其中仿真C1 不進(jìn)行修正，仿真C2 根據(jù)空氣密度對Kopt 進(jìn)行修正，仿真得到功率曲線如圖8 所示，仿真結(jié)果見表3。

圖8 風(fēng)速6 m/s 空氣密度修正前后的功率Fig.8 Electrical power before and after air density correction at wind speed 6 m/s

表3 恒定風(fēng)況仿真實驗結(jié)果Table 3 Result of simulation experiment in constant wind

根據(jù)圖8 的仿真結(jié)果，空氣密度修正前后，機組實際輸出功率未得到明顯提升。

由表3 可知，由于修正后的最優(yōu)轉(zhuǎn)矩控制系數(shù)變大，導(dǎo)致同等風(fēng)速下風(fēng)輪轉(zhuǎn)速降低，計算得到葉尖速比相應(yīng)地減小。根據(jù)圖2 所示的葉尖速比和功率系數(shù)Cp的關(guān)系，可以計算得到相應(yīng)的Cp值。修正后的功率系數(shù)與未修正的功率系數(shù)差別不大，僅提高了0.46%。與上述湍流風(fēng)況下的仿真結(jié)果一致。因此，空氣密度修正算法改變機組的風(fēng)輪轉(zhuǎn)速進(jìn)而改變?nèi)~尖速比，但是由于Cp-λ 曲線頂部平滑，對Cp影響效果小。

5 結(jié)論

本文在分析空氣密度對風(fēng)電機組功率影響的基礎(chǔ)上，結(jié)合風(fēng)電場的實際環(huán)境數(shù)據(jù)，設(shè)計了基于空氣密度修正的轉(zhuǎn)矩控制算法。在不同的空氣密度下進(jìn)行仿真實驗，并對實驗結(jié)果進(jìn)行總結(jié)分析，得到以下結(jié)論。

①風(fēng)電場環(huán)境中，空氣密度每天呈現(xiàn)周期性變化；并且，季節(jié)性變化為冬季大于夏季；整年呈近似正弦變化?？諝饷芏仁軠囟扔绊懨黠@，在風(fēng)電場沒有氣壓傳感器的條件下，可以采用實時溫度進(jìn)行近似估計。

②空氣密度修正后的轉(zhuǎn)矩控制算法對風(fēng)電機組輸出功率影響不明顯。

③空氣密度修正算法導(dǎo)致相同風(fēng)速下的轉(zhuǎn)矩、轉(zhuǎn)速值顯著改變，優(yōu)化時須考慮對機組載荷的影響。

④空氣密度的變化導(dǎo)致同一風(fēng)速下機組的發(fā)電功率不同，在風(fēng)電場功率曲線統(tǒng)計中應(yīng)注意空氣密度的修正。