文 | 付炳瑞,王海斌,胡威
在風電場基礎建設過程中,為了節(jié)省成本,減少挖平臺的土方量,經(jīng)常將山頭挖至可以澆筑風電機組基礎時就停止施工,在風電機組基礎旁留下10m-20m高未削去的半個山頭,機位微地形的改變對前期風資源評估會有一定的偏差,從而對機組運行安全造成隱患。本文將以山西某風電場為例,利用機組SCADA實測數(shù)據(jù)和CFD流體仿真手段,研究上述微地形的改變對輪轂高度處湍流和機組振動的影響。
山西某風電場共吊裝29臺機組,2014年年底完成調試并網(wǎng)運行,運行過程中部分機組出現(xiàn)機艙加速度超限故障,其中大部分機組存在基礎平臺半削平的情況。
本文選用機組運行的SCADA數(shù)據(jù),時間長度為2015年1月1日到2015年12月31日,選取場區(qū)內A、B兩臺典型故障機組和C、D兩臺正常機組,四臺機組所處的地形情況如圖1所示。故障機組與正常機組所處的地形相似,A機組和C機組與B機組和D機組各位于一道山梁,大地形下環(huán)境大體相同,因此對于研究微地形的影響具有較好的參考性。
通過處理機組SCADA數(shù)據(jù),本文將從湍流標準偏差、陣風因子、機艙加速度有效值和功率曲線四個方面討論機組局地風資源情況和運行狀態(tài)。
通過分析問題機組的故障文件可知,該項目機組發(fā)生機艙加速度超限故障的時間段主要集中在四月和十二月,即春冬季節(jié),一年中風速最大的月份也恰好在四月和十二月。
圖1 四臺機組所處的地形情況
在北方春冬季節(jié),植被稀疏,影響機位處湍流的因素以地形為主。湍流標準偏差是衡量10min樣本內數(shù)據(jù)強弱的量度,反映數(shù)據(jù)的脈動程度,用公式(1)表示:
其中ui為瞬時風速,u為樣本平均風速,n為樣本內數(shù)據(jù)個數(shù)。
在正常湍流模型中,湍流標準偏差代表值σ1應為輪轂高度處風速的90%,對于標準風力發(fā)電機等級,這個值由式(2)給出:
將實測數(shù)據(jù)計算湍流標準偏差值與代表值曲線進行對比,來驗證實際湍流特征與IEC標準適應性。以2015年4月為例,A和B機組湍流標準偏差隨時間的變化如圖3(a)和(b)所示。從圖中可知機位湍流標準偏差在大于10m/s的風速段內有較多的數(shù)據(jù)點,超出IEC標準規(guī)定的正常湍流模型中A類湍流標準偏差曲線,且這些數(shù)據(jù)點在風速維度上分布較為分散。反觀C和D正常運行的機組,其湍流標準偏差隨時間的變化如圖3(c)和(d)所示,機位湍流標準偏差在大于10m/s的風速段內的數(shù)據(jù)點超出在IEC標準規(guī)定的正常范圍。模型中超出A類湍流標準偏差曲線的數(shù)據(jù)點較少,且這些數(shù)據(jù)點在風速維度上相對集中,均在風速18m/s左右分布。四臺機組湍流標準偏差與IEC標準值的對比見表1,從表中可知,四臺機組實測湍流標準偏差在IEC湍流標準偏差A類代表值曲線以下的涵蓋比均超過90%,但A和B兩臺故障機組實測湍流標準偏差在IEC湍流標準偏差B類代表值曲線以下的占比卻低于90%。
圖2 正常湍流模型(NTM)的湍流標準偏差
表1 2015年4月四臺風電機組湍流標準偏差與IEC標準值對比
圖3 2015年4月四臺機組湍流標準偏差與風速關系圖
2015年4月1日,A和B機組因發(fā)生機艙加速度超限故障,兩臺故障機組當天的機艙加速度有效值如圖4(a)和(b)所示。由圖可知兩臺機組在部分時段內機艙加速度有效值超過0.1g,甚至達到0.16g,觸發(fā)機艙加速度超限故障的最低條件。然而同一時段內這兩臺故障機組附近的C和D機組運行正常,其機艙加速度有效值見圖4(c)和(d),C和D機組機艙加速度有效值均未超過0.1g。
通常定義陣風因子為陣風持續(xù)期tg內平均脈動風速的最大值與基本時段內的平均風速之比,即:
2015年4月1日四臺機組陣風因子隨時間變化的散點圖如圖5所示,從圖中可知A和B故障機組的陣風因子存在超過2的時段,而且這一時段恰好與故障發(fā)生時間相吻合;而在同一時間段內,C和D機組陣風因子均未超過2。經(jīng)統(tǒng)計風電機組實測數(shù)據(jù)樣本中陣風因子超過1.4的頻率見表2,由表可知,A和B故障機組陣風因子超過1.4的頻率占比均超過10%,且日平均陣風因子偏大。
表2 風電機組實測數(shù)據(jù)樣本中陣風因子超過1.4的頻率統(tǒng)計
圖4 2015年4月1日四臺機組機艙加速度有效值隨時間變化圖
圖5 2014年4日1日四臺機組陣風因子隨時間變化散點圖
故障機組不僅在湍流標準偏差、機艙加速度、陣風因子等方面存在異常,而且對功率輸出方面也有較大的影響,其功率曲線如圖6所示,2015年4月A和B機組的功率曲線在額定風速段內出現(xiàn)較大幅度的下降。
從現(xiàn)場地形分析,距離A機位的東側約22m處有一高約10m的半削平山包,同樣在距離B機位的東部約20m也有一高約11m的半削平山包,兩臺機組所處的地形相似。另外,兩臺機組與測風塔的相關性較好(見圖7),說明測風塔對這兩臺機組有一定的代表性,因此用測風塔數(shù)據(jù)進行分析是可行的。
經(jīng)分析2015年4月1日測風塔的風向為東偏東南,由此可知兩臺機組旁邊的半削平山頭恰好位于機組的主風向。機組葉片最底端距離地面約30m,其與附近的半削平山頭在垂直高度上有不足20m的高差,當風從山頭一側吹向風電機組時,氣流受到阻擋而發(fā)生畸變分離,不僅產(chǎn)生較大的湍流,而且導致風速沿塔筒高度方向變化很大。
借助機組SCADA數(shù)據(jù),只能監(jiān)測輪轂高度處來流變化情況,而無法獲知其他高度層的風參數(shù)以及流場變化。為了更好地解釋半削平山頭對下風向風電機組的影響作用,下文將采用CFD流體模擬仿真技術對場區(qū)地形進行流場模擬。
圖6 2015年4月A和B機組網(wǎng)側有功功率與標準功率曲線的對比
圖7 測風塔與A和B機組風速同期相關性驗證
CFD模擬仿真技術可較好地應用于風電場復雜地形的流動模擬, 再現(xiàn)復雜地形的大氣邊界層內的流動。為了更加直觀地描述氣流流經(jīng)半削平山頭后的流場變化情況,本文首先對半削平山頭實現(xiàn)3D建模,然后采用流場仿真模型對半削平山頭的流場進行仿真模擬,最后從機組平臺的垂直剖面和水平剖面分析半削平平臺對機組的影響。
圖8展示了模擬方案設置及數(shù)據(jù)傳遞流程。半削平處的小山包與機位點平臺的海拔相差10m,模擬區(qū)域設置網(wǎng)格分辨率為1m×1m×1m,提取測風塔風速、風向、氣溫、濕度數(shù)據(jù)生成邊界廓線作為模擬區(qū)域的初始條件,對半削平風電機組平臺處的流場進行模擬。
假設平臺處高度為0m,根據(jù)前述方法模擬半削平風電機組平臺垂直方向上0m-80m各高度處的環(huán)境湍流變化情況,具體模擬結果見圖9和表3。
根據(jù)CFD流體模擬仿真結果,風電機組平臺處的湍流強度在半削平小山包的10m高度附近達到最大,隨著高度的上升,環(huán)境湍流值逐漸減小,且環(huán)境湍流值在30m高度處有較大的變化,該結論與障礙物對湍流的影響結論一致,即受障礙物的影響,渦流擾動區(qū)的影響高度范圍約為障礙物高度的3倍。風電機組葉片掃風最低點所處的高度一般距地約為20m-30m,若氣流從半削平風電機組平臺處的小山包流經(jīng)下風向的風電機組平臺會發(fā)生氣流分離,山包與平臺間的海拔差越大,則氣流分離就越嚴重,影響的葉片掃風面的面積越大,進而給葉片的受力、機組的載荷帶來一定影響,最終影響機組正常發(fā)電。
假設半削平的小山包與機位點平臺處海拔相差10m,如圖10所示,在水平剖面上(其高度與平臺高度一致),由于受半削平平臺的影響,半削平的小山包對風有阻擋作用,在小山包前后均形成回流區(qū),對機組的發(fā)電量和載荷安全產(chǎn)生影響。
原山頭處垂直剖面上環(huán)境湍流見圖11,半削平風電機組平臺垂直剖面上環(huán)境湍流見表4。
由表4可知,從原山頭與半削平風電機組平臺環(huán)境湍流值變化情況來看,環(huán)境湍流值突變的高度層發(fā)生在風與地表接觸的那一層高度,即山地地形中貼近地表處的湍流值最大,貼近地表處的湍流對機組有一定的影響,從而增加機組振動的可能性。
圖8 地形模擬仿真流程
圖9 半削平風電機組平臺垂直剖面上環(huán)境湍流
表3 半削平風電機組平臺垂直方向上0m-80m處的環(huán)境湍流變化情況
圖10 半削平風電機組平臺水平剖面上環(huán)境湍流情況(風流場平面高度與平臺高度一致)
圖11 原來山頭處垂直剖面上環(huán)境湍流
通過對機組SCADA數(shù)據(jù)和CFD仿真模擬的分析,位于機組上風向的半削平山頭前后形成較大湍流的回流區(qū),特別容易對機組下旋葉片的氣流造成擾動。從發(fā)電量角度來看,氣流擾動會影響葉片正常變槳,使其無法達到理論最優(yōu)Cp值,最終導致機組出力較差,發(fā)電量降低;從機組載荷適應性角度來看,氣流擾動會引起機組振動、葉片壽命減少等問題。風速越大,氣流流經(jīng)半削平山頭產(chǎn)生的脫流效應會更加強烈,最終會導致機組機艙加速度超限故障。
表4 半削平山頭與原山頭處環(huán)境湍流對比
建議在今后的風電場建設中,盡量保證機組平臺周邊平整,避免半削平山頭情況出現(xiàn),雖然在前期投入成本增加,但在后期風電場運維和發(fā)電量收益中,可以減少不必要的損失。
攝影:黃斌