文 | 葉林,馮標(biāo),鐘德華,趙重陽,蔡曉峰
功率曲線特性測試儀有效風(fēng)速補償算法的研究
文 | 葉林,馮標(biāo),鐘德華,趙重陽,蔡曉峰
隨著我國風(fēng)力發(fā)電的迅速發(fā)展,以及大型風(fēng)電機組制造本地化的不斷深入,在機組研制和運行過程中,對機組功率曲線特性進行評估顯得尤為重要。通過評估機組實際功率曲線可以反映出機組多方面特性的優(yōu)劣,對機組的設(shè)計和制造過程中的質(zhì)量控制以及后期運行過程中的技改提升有著指導(dǎo)作用。但是,生成準(zhǔn)確的功率曲線需要機組的有效風(fēng)速和有功功率,由于測量風(fēng)速的功率曲線測試儀安裝于機艙尾部,葉輪轉(zhuǎn)動所引起的尾流效應(yīng)等影響,使得實測風(fēng)速與有效風(fēng)速之間存在一定的偏差。因此,需要研究一種利用實測風(fēng)速的補償算法,計算得到有效風(fēng)速。經(jīng)過大量的數(shù)據(jù)分析,能夠證明此有效風(fēng)速能夠生成較為準(zhǔn)確的功率曲線。目前,我國機組的功率、風(fēng)速等指標(biāo)的檢測主要通過機組主控系統(tǒng)的檢測設(shè)備和風(fēng)電場測風(fēng)塔,很多單位對功率曲線的檢測已取得了一些成果,但對實測風(fēng)速的補償算法和功率曲線測試儀產(chǎn)品的研究尚處于起步階段。
本文提出了一種基于大量功率曲線測試儀和激光雷達測風(fēng)儀所采集風(fēng)速數(shù)據(jù)的相關(guān)性分析,研究該功率曲線測試儀有效風(fēng)速的補償算法,最終生成風(fēng)電機組有功功率與有效風(fēng)速的功率曲線,從而實現(xiàn)準(zhǔn)確評估風(fēng)電機組風(fēng)能利用率和實際發(fā)電效率的目的。
一、數(shù)據(jù)測量對象
數(shù)據(jù)測試對象為3臺雙饋型1.5MW機組,分別為機組#1、機組#2和機組#3,均為山地高原型機組。每臺機組的機艙尾部裝有2個機械式風(fēng)杯、一套功率曲線測試儀,機組#1與機組#2的側(cè)面平行位置裝有一臺激光雷達測風(fēng)儀。
二、數(shù)據(jù)采集工具
(一)功率曲線特性測試儀。它能獨立地測量風(fēng)電機組的有功功率與風(fēng)速等數(shù)據(jù),精度:CT 0.2級,電能檢測模塊精度:0.5級,超聲波風(fēng)速儀偏差:±3%,功率曲線特性測試儀結(jié)構(gòu)圖如圖1所示。
(二)激光雷達測風(fēng)儀。Zephir公司的ZP3000測風(fēng)儀,可采集機組附近各度梯度的自由風(fēng)速、風(fēng)向、氣壓、湍流強度等數(shù)據(jù),風(fēng)速測量精度:0.5%。
(三)風(fēng)電機組主控系統(tǒng)。機組主控系統(tǒng)采集機組的有功功率和風(fēng)速等數(shù)據(jù)。
三、數(shù)據(jù)測量方法
(一)數(shù)據(jù)收集。同步采集數(shù)據(jù)的基準(zhǔn)時鐘,匯總功率曲線測試儀和風(fēng)電機組主控系統(tǒng)采集的有功功率和風(fēng)速,以及激光雷達測風(fēng)儀采集的風(fēng)速,采樣時間上具有分段連續(xù)性。測風(fēng)儀安裝位置如圖2,滿足IEC61400和GL標(biāo)準(zhǔn)。
(二)數(shù)據(jù)平均值。由于風(fēng)速具有較大的隨機性,考慮風(fēng)速儀的風(fēng)速響應(yīng)速度等因素,對采集數(shù)據(jù)作了濾波處理,計算10min平均值。
(三)風(fēng)速切變。考慮了不同高度梯度下風(fēng)速變化的影響,選擇與功率曲線測試儀處于同一高度層的測風(fēng)儀數(shù)據(jù)進行比對分析。
(四)湍流強度。數(shù)據(jù)分析中已充分考慮湍流強度等因素的影響。
(五)數(shù)據(jù)處理。根據(jù)風(fēng)電機組、測試設(shè)備的運行狀況對采集數(shù)據(jù)進行篩選,剔除無效數(shù)據(jù)。
一、數(shù)據(jù)分析工具
數(shù)據(jù)分析工具采用Matlab擬合工具箱cftool。
二、功率測試儀與風(fēng)電機組主控系統(tǒng)所采集數(shù)據(jù)相關(guān)性分析
數(shù)據(jù)相關(guān)性分析,旨在說明功率曲線測試儀和風(fēng)電機組主控系統(tǒng)所采集數(shù)據(jù)的相關(guān)性較強(相關(guān)性指兩個變量之間的密切程度),確保采集數(shù)據(jù)切實可用。
(一)風(fēng)速的相關(guān)性分析
通過分析工具,生成功率曲線測試儀與主控系統(tǒng)所采集風(fēng)速的散點圖。X1為功率曲線測試儀采集風(fēng)速,Y為風(fēng)電機組主控系統(tǒng)采集風(fēng)速。圖3為風(fēng)電機組#1的風(fēng)速散點圖,大量散點分布在對角線區(qū)域。
cftool工具線性擬合方法可以得到相關(guān)性參數(shù);風(fēng)電機組#1主控系統(tǒng)采集風(fēng)速與功率曲線特性測試儀采集風(fēng)速之間關(guān)系為Y=1.022 X1,它的風(fēng)速線性關(guān)系圖如圖4所示。
對采集風(fēng)速做相關(guān)性分析,得到相關(guān)系數(shù)為ρX1Y=0.964,根據(jù)相關(guān)系數(shù)范圍:0.8<ρX1Y≤1,具有極強相關(guān)性。
(二)有功功率數(shù)據(jù)相關(guān)性分析
對功率曲線測試儀和風(fēng)電機組主控系統(tǒng)采集的有功功率進行比對分析,P1為功率曲線測試儀采集的有功功率,P2為機組主控系統(tǒng)采集的有功功率。
機組#1的有功功率散點圖如下圖5所示,散點集中在對角線區(qū)域。通過線性關(guān)系擬合,兩者之間的關(guān)系為P2=1.062 P1,有功功率線性關(guān)系圖如圖6所示。
對采集數(shù)據(jù)點做相關(guān)性分析,得到相關(guān)系數(shù)為ρP1P2=0.957,根據(jù)相關(guān)系數(shù)范圍:0.8<ρP1P2≤1,具有極強相關(guān)性。
三、功率測試儀與激光雷達測風(fēng)儀采集風(fēng)速相關(guān)性分析
圖7為機組#1的功率曲線測試儀與激光雷達測風(fēng)儀采集風(fēng)速進行對比分析,顯示激光雷達測風(fēng)儀采集風(fēng)速比功率曲線測試儀采集風(fēng)速值大,擬合線性函數(shù):Z=1.197X1,X1為功率曲線測試儀采集風(fēng)速,Z為激光雷達測風(fēng)儀采集風(fēng)速。對兩組數(shù)據(jù)做相關(guān)性分析,相關(guān)系數(shù)為ρX1Z=0.684。
從兩組數(shù)據(jù)做的相關(guān)性分析和圖7的數(shù)據(jù)顯示可以看出,功率測試儀與激光雷達測風(fēng)儀兩個設(shè)備所采集到的風(fēng)速并沒有極強的相關(guān)性。
根據(jù)IEC61400-12標(biāo)準(zhǔn)對風(fēng)速測量的具體要求,考慮到IEC61400-12標(biāo)準(zhǔn)是對相對嚴格的型式試驗而作的規(guī)定,因此,為提高功率曲線測試儀的便利性,本文采用對機組尾部實測風(fēng)速數(shù)據(jù)進行補償,來達到接近IEC61400-12標(biāo)準(zhǔn)的方法。
激光雷達測風(fēng)儀采集風(fēng)電場的自由風(fēng)速,功率曲線測試儀采集機組尾流風(fēng)速,通過風(fēng)速補償算法得到風(fēng)電場自由風(fēng)速與功率曲線測試儀采集風(fēng)速之間的擬合關(guān)系,以達到補償功率曲線測試儀采集風(fēng)速、計算有效風(fēng)速的目的。
風(fēng)速補償算法以功率曲線測試儀采集風(fēng)速為輸入,激光雷達測風(fēng)儀采集風(fēng)速為輸出。擬合時使用樣本點分布在機組運行時的風(fēng)速范圍,數(shù)據(jù)來源于被檢測機組#1和機組#2。
根據(jù)散點分布情況做分段擬合,得到如下函數(shù)模型,X2為補償后的功率曲線特性測試儀風(fēng)速:
根據(jù)風(fēng)速補償算法的擬合函數(shù)計算得到機組#2某時間段補償后的有效風(fēng)速,補償后的有效風(fēng)速與激光雷達測風(fēng)儀采集風(fēng)速曲線如圖9所示。圖9中顯示為補償后的有效風(fēng)速與激光雷達測風(fēng)儀采集風(fēng)速較為接近,擬合線性函數(shù):Z=1.0565 X2,計算相關(guān)系數(shù)為ρX2Z=0.836,相關(guān)性有所提高。
由圖7和圖9兩圖對照可以看出,有效風(fēng)速補償算法提高了功率測試儀與激光雷達測風(fēng)儀兩個設(shè)備所采集風(fēng)速的相關(guān)性,增強了實測風(fēng)速的精確度,最終實現(xiàn)準(zhǔn)確評估風(fēng)電機組風(fēng)能利用率和實際發(fā)電效率的目的。
根據(jù)功率曲線測試儀的風(fēng)速補償算法,通過擬合函數(shù)得到補償后的有效風(fēng)速,生成基于有效風(fēng)速的功率曲線,此功率曲線即為功率曲線特性測試儀所生成的功率曲線。
圖10和圖11分別為云南某風(fēng)電場兩臺1.5MW機組(下文分別稱機組#1和風(fēng)電機組#2)的功率曲線特性測試儀所測功率曲線、風(fēng)電機組主控功率曲線(備注:風(fēng)電機組主控功率曲線為整機廠商SCADA軟件所生成的功率曲線。)和機型保證功率曲線。
下文分別以風(fēng)電機組主控功率曲線低于保證功率曲線
和主控功率曲線高于保證功率曲線的兩種情況為例,說明使用功率曲線特性測試儀可以準(zhǔn)確地檢測此兩臺風(fēng)電機組的實際功率曲線特性。
圖10為機組#1的主控功率曲線低于保證功率曲線,主控功率曲線對比保證功率曲線,發(fā)現(xiàn)主控功率曲線在風(fēng)速<10m/s時,接近合同的保證功率曲線,在風(fēng)速>10m/s時,主控功率曲線明顯低于保證功率曲線(大概50kW)。功率曲線特性測試儀所測功率曲線明顯低于保證功率曲線和主控功率曲線,爬坡段的功率曲線偏差較多,說明經(jīng)功率曲線特性測試儀檢測后證明,機組#1的主控功率曲線不合格。分析原因后發(fā)現(xiàn)該機組采用轉(zhuǎn)速-轉(zhuǎn)矩查表法控制,改進主控程序控制算法后機組功率曲線在爬坡段有了大幅提升,同時對該機組風(fēng)速儀進行了維護和程序升級。
圖11為風(fēng)電機組#2的主控功率曲線高于保證功率曲線,主控功率曲線對比保證功率曲線,發(fā)現(xiàn)主控功率曲線在整個風(fēng)速段內(nèi),主控功率曲線都高于保證功率曲線,說明該風(fēng)電機組的功率曲線合格。當(dāng)風(fēng)速<10m/s時,功率曲線特性測試儀所測功率曲線高于保證功率曲線;當(dāng)風(fēng)速>10m/s時,功率曲線特性測試儀所測功率低于保證功率曲線,說明經(jīng)功率曲線特性測試儀檢測后證明,風(fēng)電機組#2的主控功率曲線不合格。分析原因后發(fā)現(xiàn)該機組主控功率測量值為發(fā)電機機端功率,而功率曲線特性測試儀測量的是機組送出功率(變流器網(wǎng)側(cè)功率),將該機組額定功率限值少量增加至1.55MW。
功率曲線測試儀能夠獨立地檢測且直接反映風(fēng)電機組的風(fēng)能利用率和實際發(fā)電效率,對風(fēng)電機組研制和技術(shù)改進有著十分重要的意義。根據(jù)風(fēng)電機組功率曲線測試儀研究的要求,本文提出了一種有效風(fēng)速的補償算法。從仿真結(jié)果可以看出,有效風(fēng)速的補償算法可以較準(zhǔn)確地補償實測風(fēng)速的偏差,并生成更符合實際工況的功率曲線,對已投運風(fēng)電機組功率曲線的測定和機組開展技改工作都有著十分重要的作用。
(作者單位:葉林,馮標(biāo):華能洱源風(fēng)力發(fā)電有限公司;鐘德華:福氏新能源技術(shù)(上海)有限公司;趙重陽,蔡曉峰:沈陽華創(chuàng)風(fēng)能有限公司)