李穎, 雷熠, 莫爾兵, 王其君, 蘭杰
(東方電氣風(fēng)電有限公司, 四川德陽, 618000)
近年來, 隨著激光雷達技術(shù)的發(fā)展, 其在工業(yè)領(lǐng)域的應(yīng)用程度得到廣泛提高, 包括航天航空、汽車行業(yè)、 氣象領(lǐng)域、 災(zāi)難預(yù)警、 城市交通安全等等。 在風(fēng)能行業(yè), 激光雷達在陸上和海上得到了越來越廣泛的應(yīng)用。 激光測風(fēng)雷達主要運用在風(fēng)場投資前期的風(fēng)資源勘測階段、 進行單機或風(fēng)場級控制的運營階段、 以及風(fēng)機全生命周期內(nèi)的功率特性測試等等。 然而, 由于激光雷達測風(fēng)技術(shù)的數(shù)據(jù)有效性、 完整性以及成本較高等問題,并沒有得到風(fēng)電行業(yè)的普遍認(rèn)可。
激光雷達是通過發(fā)出激光光束并檢測碰到空氣流場中的微粒或者氣溶膠產(chǎn)生反向散射信號的一種激光多普勒效應(yīng)風(fēng)速計。 根據(jù)多普勒效應(yīng),在發(fā)出和返回的光束方向上, 給出了反射粒子速度分量的直接測量, 而這些測量信號的頻率會發(fā)生變化。 激光雷達根據(jù)發(fā)出信號的類型又分為脈沖激光雷達和連續(xù)波激光雷達。 激光雷達測風(fēng)原理存在一些前提性的假設(shè), 而這些假設(shè)往往就是導(dǎo)致激光雷達測風(fēng)數(shù)據(jù)完整性不高、 有效性不高的原因。 這些假設(shè)包括測量點分布空間內(nèi)流場速度均勻, 空間中氣溶膠顆粒均勻分布, 且空氣流動主導(dǎo)粒子運動, 其運動速度為空氣流動速度。因此, 空氣中氣溶膠顆粒的含量、 風(fēng)場的天氣狀況會對激光測風(fēng)雷達的測量造成較大影響。
選取1 臺安裝機艙式激光測風(fēng)雷達的風(fēng)機,進行長期數(shù)據(jù)采集和分析。 測試風(fēng)機位于內(nèi)蒙古風(fēng)場, 風(fēng)場地貌屬于壩上高平原地區(qū)的丘陵地段,測試時間為3 個月。 雷達選用某機艙式激光雷達,光束結(jié)構(gòu)為4 光束、 采樣頻率為4 Hz。 此機艙式激光雷達提供記錄風(fēng)力發(fā)電機組機艙位置射出前方50~200 m 共10 個測量距離的風(fēng)場信息數(shù)據(jù)。
徑向風(fēng)速表示空間風(fēng)矢量在雷達4 個激光束的投影風(fēng)速, 是雷達的原始測量數(shù)據(jù)。 輪轂高度處合成風(fēng)速為4 個激光束在某個距離點測得風(fēng)速信息的合成風(fēng)速, 因此, 輪轂高度處合成風(fēng)速的有效率能夠反應(yīng)4 光束同時被遮擋的情況和信噪比隨距離變化的情況, 也能反應(yīng)激光雷達測量風(fēng)速數(shù)據(jù)的可用程度。 定義nνhub為輪轂高度處合成風(fēng)速的有效數(shù)據(jù)(排除采集周期內(nèi)NaN 的數(shù)據(jù)), n為采集數(shù)據(jù)的總數(shù)。 (影響激光雷達系統(tǒng)數(shù)據(jù)可用性的任何條件不包括在內(nèi): 電源中斷、 自然行為導(dǎo)致的系統(tǒng)損壞、 通訊中斷、 維護等等。 )則激光雷達測量風(fēng)速數(shù)據(jù)的有效率計算公式為:
數(shù)據(jù)可利用率定義為在整個時間戳范圍內(nèi)的可用的合理性數(shù)據(jù)nνhub_R(排除采集周期內(nèi)偏差過大數(shù)據(jù), 排除NaN 數(shù)據(jù)) 占有效數(shù)據(jù)的百分比。(影響激光雷達系統(tǒng)數(shù)據(jù)可用性的任何條件不包括在內(nèi): 電源中斷、 自然行為導(dǎo)致的系統(tǒng)損壞、 通訊中斷、 維護等等。 )則激光雷達測量風(fēng)速數(shù)據(jù)的可利用率計算公式為:
通過對激光雷達數(shù)據(jù)進行有效性和可利用率計算, 如表1 所示, 說明激光雷達風(fēng)速信息的可利用率滿足需求, 但在整個觀測期內(nèi)的有效率不高, 超出可接受范圍, 不建議作為實時控制的輸入信息。
由于測試?yán)走_的測量范圍為50~200 m, 對相近測量距離的數(shù)據(jù)進行對比, 分析雷達自身數(shù)據(jù)的關(guān)聯(lián)性, 從而確認(rèn)雷達本身數(shù)據(jù)的可靠性和正確性。 測試數(shù)據(jù)選取2019 年 12 月至2020 年 2 月共3 個月的數(shù)據(jù), 在保證數(shù)據(jù)為相同時間段內(nèi)的滑動平均值、 數(shù)據(jù)無時移的前提下, 篩選掉無數(shù)據(jù)和不合理數(shù)據(jù), 進行激光雷達射出距離分別為50 m、 60 m、 80 m、 100 m、 120 m、 140 m、 160 m、 180 m、 190 m、 200 m 共 10 個距離段兩兩相近距離(共9 個比較結(jié)果)的風(fēng)速/風(fēng)向時序?qū)Ρ群拖嚓P(guān)性分析。 相關(guān)性分析采用線性回歸算法模型進行曲線擬合,并用最小二乘法求解損失函數(shù)的最小化問題,擬合程度的評估值用判定系數(shù)R2表示。
此線性回歸預(yù)測模型可定義為假設(shè)存在一條最佳擬合的直線方程: y=ax+b, 對于每一個樣本點Xi代入直線方程, 可得到預(yù)測值Y^i=aXi+b。 且樣本真實值Yi與預(yù)測值Y^i的距離要盡可能的小。則此問題就是最小化損失函數(shù)即可以定義為每一個樣本點和預(yù)測值的距離總和最小, 即:
式中n 為樣本數(shù)量。
通過此模型擬合方法將每個10 個距離分段數(shù)據(jù)共9 組對比結(jié)果計算得出的模型擬合函數(shù)和判定系數(shù)如表2~3 所示。
表2 10 min 平均風(fēng)速相關(guān)性分析表
表3 10 min 平均風(fēng)向相關(guān)性分析表
10 分鐘平均風(fēng)速/風(fēng)向時序圖不同距離進行了對比。
對比數(shù)據(jù)分析顯示, 激光雷達測得不同距離處風(fēng)速/風(fēng)向的相關(guān)性程度高, 縱向說明了雷達設(shè)備本身測量數(shù)據(jù)的可靠性與真實性。 后續(xù)還需要將雷達數(shù)據(jù)與機艙風(fēng)速計、 機艙風(fēng)向標(biāo)進行橫向?qū)Ρ龋?確認(rèn)不同測量設(shè)備測量數(shù)據(jù)的相關(guān)程度。由于雷達不同距離相關(guān)性程度高, 只需要將距離風(fēng)輪面最近的50 m 雷達風(fēng)速數(shù)據(jù)與機艙風(fēng)速計、機艙風(fēng)向標(biāo)進行對比即可。
現(xiàn)階段, 單臺風(fēng)力發(fā)電機組測量風(fēng)速風(fēng)向,往往選取經(jīng)濟性高的機械式或超聲波風(fēng)速計。 那么機艙式激光測風(fēng)雷達的測量數(shù)據(jù)需要跟同類產(chǎn)品的測量數(shù)據(jù)進行橫向?qū)Ρ龋?確認(rèn)其關(guān)聯(lián)程度,從而確認(rèn)雷達測得風(fēng)速風(fēng)向的真實性。 同樣地,選取相同時間段雷達數(shù)據(jù)與機艙風(fēng)速計和風(fēng)向標(biāo)數(shù)據(jù)進行采集, 剔除掉無效數(shù)據(jù), 采取線性回歸模型對其相關(guān)性進行分析可以得到風(fēng)速風(fēng)向時序圖對比和相關(guān)性分析圖, 如圖1~2 所示。
以上數(shù)據(jù)分析表明, 激光雷達與機艙風(fēng)速計測得風(fēng)速相關(guān)程度較高, 但激光雷達與風(fēng)向標(biāo)測得風(fēng)向數(shù)據(jù)相關(guān)程度較差,擬合程度只有0.525 2。
激光雷達根據(jù)安裝位置又可分為機艙式、 地面式。 根據(jù)不同的應(yīng)用場景, 選取不同類型的激光雷達。 比如, 地基式激光雷達可作為測風(fēng)塔的替代方案, 進行微觀場址的風(fēng)資源評估。 而機艙式激光雷達可進行單臺風(fēng)機的風(fēng)速預(yù)測, 進行前饋控制、 尾流控制等。
單臺機組可以通過離線或在線的方式, 收集一段時間的數(shù)據(jù)量進行單臺機組的偏航校正。 根據(jù)雷達測得的水平距離50 m 處的風(fēng)速和風(fēng)向分布關(guān)系, 分布圖如圖3 所示。
圖3 雷達測量風(fēng)速/風(fēng)向分布圖 (10 min 風(fēng)速)
根據(jù)風(fēng)速風(fēng)向分布, 定義偏航角度誤差計算公式為:
式中Di為每個有效風(fēng)向,νi為每個有效風(fēng)速。
通過采集3 個月10 min 平均風(fēng)速風(fēng)向數(shù)據(jù),可計算出該風(fēng)機的平均偏航角度誤差為+0.846 7°(雷達風(fēng)向數(shù)據(jù)以雷達射線方向左邊為正)。 此誤差通過一段時間數(shù)據(jù)采集統(tǒng)計, 由修改主控參數(shù)實現(xiàn)。 此方法可以與設(shè)備定檢同步進行, 間隔一段時間校正偏航誤差, 提高全風(fēng)場對風(fēng)精確度,從而提升機組的發(fā)電效率。
采集1 個月的雷達風(fēng)速、 機艙風(fēng)速計風(fēng)速、有功功率的數(shù)據(jù)進行分析。 篩選掉無效數(shù)據(jù)的前提下, 通過對測試風(fēng)力發(fā)電機組雷達測得10 min平均風(fēng)速與機艙風(fēng)速計測得10 min 平均風(fēng)速進行功率曲線對比數(shù)據(jù), 可以看到, 在同等功率下,機艙風(fēng)速計統(tǒng)計的功率曲線比雷達統(tǒng)計的功率曲線有1 m/s 左右的過高估計, 如圖4 所示。 而在激光雷達風(fēng)速、 風(fēng)速計風(fēng)速的時序?qū)Ρ戎胁浑y看出,風(fēng)速計風(fēng)速時序圖略高于激光雷達測得風(fēng)速時序圖。 由此可見, 風(fēng)速計測得風(fēng)速無論作為秒級數(shù)據(jù)查看機組運行曲線, 或者用作10 min 平均數(shù)據(jù)作為機組功率曲線統(tǒng)計, 過高的風(fēng)速估計會導(dǎo)致功率曲線過低。 而激光雷達數(shù)據(jù)在數(shù)據(jù)準(zhǔn)確度方面優(yōu)勢較為明顯, 統(tǒng)計功率曲線更具優(yōu)勢。
圖4 風(fēng)速計與雷達測風(fēng)功率曲線對比圖
本文通過對機艙式測風(fēng)雷達數(shù)據(jù)進行多維度多角度的初步分析, 并對雷達數(shù)據(jù)在風(fēng)力發(fā)電機組的應(yīng)用進行舉例分析, 評價雷達在風(fēng)能行業(yè)的可用性程度。 就3 個月試驗數(shù)據(jù)而言, 可得出以下結(jié)論:
(1)機艙式激光雷達在測量范圍內(nèi)的相關(guān)程度較高, 風(fēng)速信息的判定系數(shù)高達0.99 以上, 風(fēng)向信息的判定系數(shù)>0.95, 不同距離的測量數(shù)據(jù)準(zhǔn)確程度較高。
(2)激光雷達與機艙風(fēng)速計作對比, 從時序圖看, 機艙風(fēng)速計測量風(fēng)速普遍略高于激光雷達測量風(fēng)速, 從擬合曲線看, 風(fēng)速相關(guān)性較好, 判定系數(shù)為0.945 5。 但風(fēng)向數(shù)據(jù)的關(guān)聯(lián)程度較差, 判定系數(shù)只有0.525 2。 可以看出, 輪轂高度處的來流風(fēng)向經(jīng)過風(fēng)輪面掃風(fēng)后, 已經(jīng)發(fā)生了較大改變。
(3)激光雷達測量風(fēng)速風(fēng)向的數(shù)據(jù)有效性、 完整性較低, 其中輪轂高度處合成風(fēng)速的有效率并未達到統(tǒng)計指標(biāo)的可接受范圍, 而機艙風(fēng)速計的測量數(shù)據(jù)完整性高。 激光雷達測得風(fēng)速是風(fēng)輪面前方的風(fēng)速信息, 但由于風(fēng)輪遮擋和空氣氣溶膠顆粒不足等原因?qū)е吕走_測得風(fēng)速風(fēng)向存在無效數(shù)據(jù)。 因此, 雷達數(shù)據(jù)作為風(fēng)機單機控制的風(fēng)險較大。 而機艙式風(fēng)速計測得風(fēng)速風(fēng)向為風(fēng)輪面后的尾流, 雖然數(shù)據(jù)完整性較高, 但精準(zhǔn)程度存在較大誤差, 并非實際風(fēng)輪面的來流風(fēng)速。
(4)根據(jù)實際雷達數(shù)據(jù)的應(yīng)用與雷達有效性相結(jié)合分析,雷達數(shù)據(jù)完整性程度不高, 不建議作為實時控制的測量信號。 但由于其數(shù)據(jù)的準(zhǔn)確性較高, 對于有一定時間間隔的離線控制或數(shù)據(jù)分析具有一定指導(dǎo)意義。
綜上所述, 在風(fēng)能行業(yè), 激光雷達具有測量范圍大、 運用場景多、 數(shù)據(jù)準(zhǔn)確性高等優(yōu)點, 但由于數(shù)據(jù)完整性、 有效性不高, 市場價格高、 設(shè)備壽命短等特點, 許多整機廠商認(rèn)為其設(shè)備作為風(fēng)機控制信號還需觀望。 但在科技高速發(fā)展的今天, 隨著產(chǎn)業(yè)配套的加速發(fā)展, 相信不久的將來,雷達測風(fēng)設(shè)備可以像車載雷達一樣, 成為風(fēng)電機組的標(biāo)準(zhǔn)配置, 為智慧風(fēng)電做出必要的貢獻。