王 印 劉飛虹 吳 迪 苗得勝 陳俊鵬

明陽智慧能源集團股份公司

0 引言

在實現(xiàn)“雙碳”目標的過程中，清潔能源產(chǎn)業(yè)的發(fā)展得到了大力支持，海上風電作為清潔能源的重要組成部分，對于實現(xiàn)能源結(jié)構(gòu)轉(zhuǎn)型和應(yīng)對氣候變化具有重要意義。我國具有綿長的海岸線、廣闊的可利用海域面積和豐富的海上風能資源，海上風電產(chǎn)業(yè)呈快速發(fā)展態(tài)勢，但我國沿海地區(qū)臺風高發(fā)。臺風特性參數(shù)與海上風電場的設(shè)計與控制息息相關(guān)，特別是風向參數(shù)，直接參與風機的運行與調(diào)控[1］，臺風期間風向驟變是風力機倒塌的主要原因之一[2］，因而對臺風期間風向變化特性研究極為重要。

已有的臺風特性研究多集中在臺風影響期間風速的脈動特性、湍流強度、風向角等方面[3-5］，對臺風期間風向特性的研究較少。常蕊等[6］利用數(shù)值模擬的方法研究30 min 平均下，相鄰時刻風向變化角度沿臺風半徑方向的變化規(guī)律及臺風眼周圍最大風向變化角度分布；朱云輝等[7］對廈門地區(qū)150 m高空的四次臺風過境時，2 h 數(shù)據(jù)內(nèi)10 min 平均下風向角的整體變化范圍進行對比，發(fā)現(xiàn)四次臺風的風向角變化范圍差異較大；陳秋陽等[8］對比分析了3 座海上測風塔測量同一場臺風過境時的數(shù)據(jù)，發(fā)現(xiàn)10 min 風向偏轉(zhuǎn)角度隨著距臺風中心距離的增加呈逐漸減小的趨勢。何如等[9］利用超聲波測風儀測量臺風影響期間55 m 高度處的風速風向數(shù)據(jù)，并針對風向偏轉(zhuǎn)變化進行了分析。麥梁詩[10］通過沿海橋梁上記錄的3 場完整的臺風過境期間的數(shù)據(jù)，分析了日常平均風向、標準差，與臺風期間的風向、標準差的差異，但缺少對相鄰時刻風向的瞬變分析，不能很好地代表臺風期間海上風電場內(nèi)的風況。

以上有關(guān)臺風期間風向的研究多集中在臺風的整體風向變化，極少有涉及臺風期間的風電場內(nèi)的風向變化，且臺風期間的風向研究仍基于10 min時間精度以上的數(shù)據(jù)分析，不能體現(xiàn)出風向的驟變過程。因此，基于以上學者的研究成果，本文將依據(jù)多個臺風過境期間，海上風電場內(nèi)所記錄到的高精度風速風向數(shù)據(jù)，分析1 s 平均及30 s 平均下海上風電場的風向偏轉(zhuǎn)特性，為海上風電場的設(shè)計及臺風期間風電機組的安全運行提供參考。

1 數(shù)據(jù)說明

1.1 數(shù)據(jù)來源

統(tǒng)計數(shù)據(jù)源于臺風“山竹”“海高斯”“圓規(guī)”及“查帕卡”登陸我國廣東沿海地區(qū)期間，臺風路徑沿途風電場內(nèi)測風設(shè)備記錄下的數(shù)據(jù)，測風設(shè)備安裝于風電機組機艙的頂部，記錄的主要要素包括：風速、風向、機組運行模式、機艙對北角度、風向?qū)Ρ苯嵌鹊?，?shù)據(jù)記錄頻率為1 Hz。本文分析的多場臺風及其所影響的風電場信息見表1，臺風移動路徑與受影響風電場的位置關(guān)系如圖1所示，圖2為風電場A在臺風“山竹”影響下所記錄的原始時程數(shù)據(jù)。

圖2 “山竹”影響下風電場A風速風向時程數(shù)據(jù)

表1 臺風及其影響風電場的信息統(tǒng)計

1.2 數(shù)據(jù)的篩選與統(tǒng)計方法

對于臺風過境期間風電場所記錄到的風速風向數(shù)據(jù)，首先識別出其中的無效數(shù)據(jù)，并對其進行剔除處理，處理完畢后有效數(shù)據(jù)完整率如表2所示。

表2 數(shù)據(jù)統(tǒng)計信息

風向偏轉(zhuǎn)角表示的是當前時間點與前一個時間點的風向差值，本文統(tǒng)計時間點的尺度分別為1 s平均與30 s 平均。在統(tǒng)計風向偏轉(zhuǎn)角時，忽略數(shù)據(jù)斷點間的風向變化。

為了便于下文分析，根據(jù)臺風中心風力等級[11］，確定風速特征參量（見表3），并根據(jù)統(tǒng)計數(shù)據(jù)的工程應(yīng)用可參考性，確定數(shù)據(jù)的統(tǒng)計特征參量（見表4）。風電機組葉輪的旋轉(zhuǎn)會使位于機艙頂部的測風設(shè)備記錄到的風速風向數(shù)據(jù)產(chǎn)生較大偏差，因此重點統(tǒng)計分析強風速段所對應(yīng)的風向偏轉(zhuǎn)特性，此時風電機組處于停機狀態(tài)，不僅能降低葉輪旋轉(zhuǎn)對分析結(jié)果的影響，還能分析獲取到機組在強風速下承受極限載荷時的風向偏轉(zhuǎn)特性。

表3 風力強度分級

表4 統(tǒng)計特征參量

2 風向偏轉(zhuǎn)特性分析

本節(jié)主要對臺風過境期間，風電場內(nèi)不同風速段內(nèi)風向偏轉(zhuǎn)角的大小進行分析，探索具有典型特征的臺風過境時，風電場的風向偏轉(zhuǎn)特征，并對多場臺風下風向偏轉(zhuǎn)統(tǒng)計特征及其隨風速的變化規(guī)律進行總結(jié)。

2.1 風向偏轉(zhuǎn)數(shù)據(jù)的統(tǒng)計特點

對臺風“山竹”過境期間風電場A的記錄數(shù)據(jù)進行讀取，風速范圍為2.5~57.5 m/s，數(shù)據(jù)的時間精度為1 s，按5 m/s 劃分風速區(qū)間，統(tǒng)計各風速區(qū)間內(nèi)，處于不同風向偏轉(zhuǎn)角范圍的數(shù)據(jù)量。

不同風速強度下的風向偏轉(zhuǎn)量如圖3 所示，風向偏轉(zhuǎn)角主要分布在2°以內(nèi)，且隨著風速的增加，風向偏轉(zhuǎn)角呈減小趨勢。各風速區(qū)間內(nèi)風向偏轉(zhuǎn)角的統(tǒng)計結(jié)果表明，在全風速段下，絕大部分風向偏轉(zhuǎn)角小于18°，強風速段的（風速大于32.5 m/s）主體風向偏轉(zhuǎn)角小于10°，而超強風速段（風速大于47.5 m/s）的主體風向偏轉(zhuǎn)角小于8°，圖3（b）和（c）清楚地展現(xiàn)出上述變化。

圖3 不同風速強度下的風向偏轉(zhuǎn)量統(tǒng)計

2.2 典型風向偏轉(zhuǎn)變化分析

風電機組的偏航控制需同時關(guān)注1 s 平均風向及30 s 平均風向，因此，分別獲取臺風“山竹”過境期間風電場A 中#1 機組、臺風“查帕卡”過境期間風電場C 中#1 機組所記錄的風速風向數(shù)據(jù)，對臺風過境期間風電場中1 s平均及30 s平均的典型風向偏轉(zhuǎn)變化特性進行分析。

2.2.1 1 s平均統(tǒng)計分析

臺風過境時1 s 平均下的風速、風向偏轉(zhuǎn)隨時序變化如圖4 和圖5 所示，在此基礎(chǔ)上，統(tǒng)計可對風電機組安全性造成影響的強風（≥35 m/s）下風向偏轉(zhuǎn)角特征參數(shù)。

圖4 “山竹”過境1 s平均下的風速、風向偏轉(zhuǎn)變化

圖5 “查帕卡”過境1 s平均下的風速、風向偏轉(zhuǎn)變化

臺風“山竹”過境期間，風電場A中風速風向1 s平均下的特征參量統(tǒng)計結(jié)果見表5。臺風過境期間所有1 s平均下大于35 m/s風速下風向偏轉(zhuǎn)角統(tǒng)計數(shù)據(jù)表明，偏轉(zhuǎn)角的最大值為26.10°，對應(yīng)的風速為40.01 m/s；偏轉(zhuǎn)角的99%分位數(shù)對應(yīng)的值為7.55°，對應(yīng)的風速為36.05 m/s；偏轉(zhuǎn)角的95%分位數(shù)對應(yīng)的值為4.07°，對應(yīng)的風速為42.40 m/s。

表5 “山竹”過境1 s平均下的特征參量

臺風“查帕卡”過境期間，風電場C 中風速風向1 s 平均下的特征參量統(tǒng)計結(jié)果見表6。臺風過境期間所有1 s 平均下大于35 m/s 風速下風向偏轉(zhuǎn)角統(tǒng)計數(shù)據(jù)表明，偏轉(zhuǎn)角的最大值為23.14°，對應(yīng)的風速為42.03 m/s；偏轉(zhuǎn)角的99%分位數(shù)對應(yīng)的值為11.52°，對應(yīng)的風速為37.40 m/s；偏轉(zhuǎn)角的95% 分位數(shù)對應(yīng)的值為7.92°，對應(yīng)的風速為63.12 m/s。

表6 “查帕卡”過境1 s平均下的統(tǒng)計參數(shù)

2.2.2 30 s平均統(tǒng)計分析

臺風過境時30 s 平均下的風速、風向偏轉(zhuǎn)隨時序變化如圖6和圖7所示，對于30 s平均分析，統(tǒng)計在30 m/s 以上的強風速所對應(yīng)的風向偏轉(zhuǎn)角特征參量。

圖6 “山竹”過境30 s平均下的風速、風向偏轉(zhuǎn)變化

圖7 “查帕卡”過境30 s平均下的風速、風向偏轉(zhuǎn)變化

臺風“山竹”過境期間，風電場A中風速風向30 s平均下的特征參量統(tǒng)計結(jié)果見表7。臺風過境期間所有30 s 平均下大于30 m/s 風速下風向偏轉(zhuǎn)角統(tǒng)計結(jié)果表明，偏轉(zhuǎn)角的最大值為18.45°，對應(yīng)的風速值為30.24 m/s；偏轉(zhuǎn)角的99%分位數(shù)對應(yīng)的值為11.16°，對應(yīng)的風速值為36.29 m/s；偏轉(zhuǎn)角的95%分位數(shù)對應(yīng)的值為5.57°，對應(yīng)的風速值為34.40 m/s。

表7 “山竹”過境30 s平均統(tǒng)計參數(shù)

臺風“查帕卡”過境期間，風電場C 中風速風向30 s 平均下的特征參量統(tǒng)計結(jié)果見表8，臺風過境期間所有30 s 平均下大于30 m/s 風速下風向偏轉(zhuǎn)角統(tǒng)計結(jié)果表明，偏轉(zhuǎn)角的最大值為28.12°，對應(yīng)的風速值為38.63 m/s；偏轉(zhuǎn)角的99%分位數(shù)對應(yīng)的值為13.03°，對應(yīng)的風速值為43.60 m/s；偏轉(zhuǎn)角的95%分位數(shù)對應(yīng)的值為9.24°，對應(yīng)的風速值為42.98 m/s。

表8 “查帕卡”過境30 s平均下統(tǒng)計參數(shù)

將30 s 平均與1 s 平均所統(tǒng)計出的風向偏轉(zhuǎn)特征參數(shù)進行對比，發(fā)現(xiàn)30 s 平均下99%、95%分位數(shù)的風向偏轉(zhuǎn)值均大于1 s 平均的統(tǒng)計結(jié)果，說明在臺風影響期間，風電場的風向處于一個持續(xù)偏轉(zhuǎn)的過程。

2.3 風向偏轉(zhuǎn)特征參數(shù)統(tǒng)計分析

基于多個臺風過境時沿途風電場的記錄數(shù)據(jù)，分別統(tǒng)計風向偏轉(zhuǎn)角的最大值、99%分位數(shù)、95%分位數(shù)及其所對應(yīng)的風速，1 s 平均及30 s 平均下的風向偏轉(zhuǎn)角與對應(yīng)風速的統(tǒng)計值，統(tǒng)計結(jié)果見表9~表12，基于統(tǒng)計數(shù)據(jù)得到偏轉(zhuǎn)角與對應(yīng)風速的關(guān)系如圖8和圖9所示。

圖8 偏轉(zhuǎn)角與對應(yīng)風速（1 s平均）

圖9 偏轉(zhuǎn)角與對應(yīng)風速（30 s平均）

表9 “山竹”影響下風電場A風向偏轉(zhuǎn)角與對應(yīng)風速統(tǒng)計

表10 “海高斯”影響下風電場A風向偏轉(zhuǎn)角與對應(yīng)風速統(tǒng)計

表11 “圓規(guī)”影響下風電場B風向偏轉(zhuǎn)角與對應(yīng)風速統(tǒng)計

表12 “查帕卡”影響下風電場C風向偏轉(zhuǎn)角與對應(yīng)風速統(tǒng)計

圖8 為臺風過境期間風速風向數(shù)據(jù)的1 s 平均統(tǒng)計結(jié)果，由圖8（a）可知風向偏轉(zhuǎn)角最大值與其所對應(yīng)的風速呈正相關(guān)，而99%分位數(shù)、95%分位數(shù)所對應(yīng)的風向偏轉(zhuǎn)角與風速呈負相關(guān)，30 s 平均下的風向偏轉(zhuǎn)角與對應(yīng)風速的統(tǒng)計結(jié)果也發(fā)現(xiàn)同樣的規(guī)律，如圖9（a）所示。結(jié)合不同風速強度下的風向偏轉(zhuǎn)量統(tǒng)計結(jié)果（圖3）可知，風向偏轉(zhuǎn)角有隨風速增大而減小的趨勢，即風向偏轉(zhuǎn)角最大值不具有風向隨風速偏轉(zhuǎn)變化的代表性，而99%、95%分位數(shù)統(tǒng)計結(jié)果顯示風向偏轉(zhuǎn)角隨風速的變化趨勢與散點數(shù)據(jù)分布結(jié)果具有一致性趨勢。因此，對風速的偏轉(zhuǎn)進行參量統(tǒng)計時，99%分位數(shù)統(tǒng)計結(jié)果具有代表性。將99%分位數(shù)對應(yīng)的風向偏轉(zhuǎn)角與風速值的散點數(shù)據(jù)進行線性擬合，1 s 平均和30 s 平均下擬合直線斜率分別為-0.060、-0.058，斜率相近。

以上所述數(shù)據(jù)是針對多場臺風、多個風場的統(tǒng)計參數(shù)結(jié)果，該統(tǒng)計結(jié)果能反映臺風過境期間臺風路徑沿途風電場的風向偏轉(zhuǎn)變化特征，因此，對風向偏轉(zhuǎn)角與對應(yīng)風速進行分位數(shù)統(tǒng)計，結(jié)果如圖8（b）和圖9（b）所示。統(tǒng)計結(jié)果表明，1 s 平均下，99%分位數(shù)偏轉(zhuǎn)角度最大值、中位值分別為16.38°、10.5°，對應(yīng)風速分別為27.41 m/s、57.2 m/s；30 s平均下，99%分位數(shù)偏轉(zhuǎn)角最大值、中位值分別為17.06°、13.19°，對應(yīng)風速分別為36.66 m/s、43.82 m/s。

3 結(jié)論

本文基于登陸廣東沿海的多場臺風過境期間風電場內(nèi)風電機組的記錄數(shù)據(jù)，分析臺風過境時風電場內(nèi)輪轂高度處的風向偏轉(zhuǎn)角變化規(guī)律，得到以下主要結(jié)論：

1）風向偏轉(zhuǎn)角隨著風速增大呈減小趨勢；

2）最大風向偏轉(zhuǎn)角與風速大小呈正相關(guān)，99%、95%分位數(shù)風向偏轉(zhuǎn)角與風速大小呈負相關(guān)，而99%分位數(shù)值遵循風向偏轉(zhuǎn)角變化規(guī)律，對反映風向偏轉(zhuǎn)的變化特征具有代表性；

3）根據(jù)風速風向數(shù)據(jù)的1 s 平均統(tǒng)計結(jié)果，99%分位數(shù)偏轉(zhuǎn)角最大值、中位值分別為16.38°、10.5°，對應(yīng)風速分別為27.41 m/s、57.2 m/s；

4）根據(jù)風速風向數(shù)據(jù)的30 s 平均統(tǒng)計結(jié)果，99%分位數(shù)偏轉(zhuǎn)角最大值、中位值分別為17.06°、13.19°，對應(yīng)風速分別為36.66 m/s、43.82 m/s。