李健華 黃亞玨

摘? ? 要：湛江海域是廣東省海上風電發(fā)展的重點規(guī)劃區(qū)域，數(shù)值模擬是為規(guī)模化海上風電場建設(shè)提供波浪設(shè)計參數(shù)的重要手段。本文以歐洲氣象中心數(shù)據(jù)為輸入風場，結(jié)合C-MAP和實測地形得到計算水深，基于MIKE21-SW建立了湛江海域的波浪數(shù)學模型，模擬海域2010～2020年的逐時波浪過程。基于實測數(shù)據(jù)對模型進行了驗證，波高的相關(guān)系數(shù)為0.94；統(tǒng)計了研究海域多年波浪要素特征，得到有效波高分布地圖；分析了波浪的季節(jié)變化、年際變化趨勢及地形對波浪的影響。波高同時受風和地形影響，2010～2020年間的年平均波高無明顯變化，冬季波高最大、夏季最小。本文可為當前項目及后續(xù)廣東海區(qū)繼續(xù)發(fā)展海上風電提供借鑒。

關(guān)鍵詞：波浪特征；數(shù)值模擬；海上風電場；地形變化

中圖分類號：P731.22 ?? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標識碼：A

Study on Wave Characteristics of Zhanjiang Offshore Wind Farm Sea Area

LI Jianhua1，2，? HUANG Yajue3

（ 1. China Energy Engineering Group Guangdong Electric Power Design Institute Co.， Ltd.，? ?Guangzhou 510663;? 2. Guangdong Kenuo Surveying Engineering Co.， Ltd.，? Guangzhou? 510663;? 3. Pearl River Water Resources Commission of the Ministry of Water Resources，? Guangzhou? 510610 ）

Abstract： Zhanjiang sea area is a key plan area for the development of offshore wind power in Guangdong Province. And it is important to provide wave design parameters for the scale construction of offshore wind farms based on numerical simulation. By taking the data from the European Meteorological Center as the input wind field， and obtaining the calculated water depth? in combination with C-MAP and measured terrain， the wave mathematical model of Zhanjiang sea area was established upon MIKE21-SW model in this study. The hourly wave process in the sea area from 2010 to 2020 was simulated. The model was validated based on measured data， and the correlation coefficient of wave height was 0.94. The characteristics of wave parameters at the study area were statistically analyzed， and the distribution of significant wave heights was obtained. The seasonal and inter-annual variations of wave were analyzed， as well as the impact of terrain on wave. The wave height was influenced by both wind and topographic change， and there was no significant trend in the annual average wave height during 2010-2020. However， the wave height was the highest in winter and the lowest in summer. The results could provide scientific reference for the present and future development of offshore wind farm in Guangdong coastal area.

Key words： wave characteristics; numerical simulation; offshore wind farm; topographic change

1? ? 前言

海上風力發(fā)電是清潔能源的重要組成部分[1]，其有多種類型，如固定式、漂浮式等。由于要長期經(jīng)受風、浪、流等載荷的共同作用，工作環(huán)境復(fù)雜，因此海上風電在選址時必須基于長期的海洋水文資料[2]。湛江海域是廣東省海上風電發(fā)展的重點規(guī)劃區(qū)域，已建成外羅海上風電場、徐聞海上風電場、新寮海上風電場、中國海裝湛江漂浮式示范項目等，規(guī)劃建設(shè)徐聞東海上風電場等。湛江海域具有水下地形變化大的特點，研究湛江海上風電場海域的波浪時空變化特征，是海上風電場規(guī)模化發(fā)展的迫切需求。

現(xiàn)場實測是獲取波浪資料以及分析波浪特征的主要手段[3]。由于海上風電開發(fā)逐漸向深遠海發(fā)展，對于儀器安全、數(shù)據(jù)安全、作業(yè)窗口期等方面均提出了更加嚴峻的考驗；同時長時間的觀測需要更高的勘察成本，因此短期的現(xiàn)場觀測配合數(shù)值模擬的手段，能夠彌補現(xiàn)場觀測存在的問題，保證項目的順利推進，已逐漸成為研究海洋特征的常用方法。

盡管目前已有對南海波浪場的模擬研究，但大多是對南海波浪場整體的數(shù)值模擬，湛江海上風電場海域地形變化較大，波浪要素的空間分布變化也較大，導(dǎo)致波浪場具有其獨特的特征，這些特征會對風電場的設(shè)計和運行造成影響，而目前對于湛江風電場海域的波浪場特征研究還不夠深入。

為了給湛江海上風電場海域多個風電場提供設(shè)計輸入資料，基于MIKE21-SW建立了波浪數(shù)學模型，在通過實測資料對模型進行充分驗證的基礎(chǔ)上，對研究海域進行長時序逐時波浪后報，模擬了研究海域2010～2020的長時序逐時波浪過程，統(tǒng)計了波浪場的空間分布以及年內(nèi)、年際變化規(guī)律，分析了地形對波浪場分布的影響以及波浪的季節(jié)性變化規(guī)律。本文的結(jié)論既能夠了解研究海域的波浪分布特征，又能夠為該海域的規(guī)?；L場設(shè)計提供時序更精細的基礎(chǔ)水文數(shù)據(jù)。

2? ? ?數(shù)值模型

本文基于MIKE21-SW建立波浪數(shù)學模型，MIKE21-SW是基于非結(jié)構(gòu)網(wǎng)格的新一代波浪譜模型， 主要適用于模擬河口海岸地區(qū)大范圍的風浪和涌浪的成長、衰減及變形。波浪數(shù)學模型的基本控制方程為[9]：

式中：N為動譜能量密度；σ為相對波浪頻率；θ為波向；Cx、Cy為波浪沿x、y方向傳播的速度；Cσ、Cθ為波浪在σ、θ坐標下的傳播速度；S為源匯項，由下式計算得到

S=Sin+Snl+Sds+Sbot+Ssurf ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （2）

式中：Sin為風能輸入項；Snl為非線性波-波相互作用的能量傳輸；Sds為波浪白帽耗散造成的能量損失；Sbot為波浪底部摩阻所造成的能量損失；Ssurf為波浪破碎所導(dǎo)致的能量損失。

3? ? 數(shù)值模型驗證

風場數(shù)據(jù)的選擇對模擬結(jié)果的精度有顯著影響，目前主流的風場數(shù)據(jù)有歐洲的ERA5風場、美國的NCEP風場、日本的JRA25等。其中， ERA5數(shù)據(jù)集具有較高的分辨率[10]，并在中國海域經(jīng)過廣泛驗證，具有較高的精度[11]，故選用ERA5海面10 m風場逐時再分析數(shù)據(jù)作為模型計算的風場條件，其數(shù)據(jù)空間分辨率為0.25°×0.25°，時間分辨率為逐時。

3.1? ?模型范圍

為方便描述，選擇如下5個代表點：外海點P2、海上風電場核心區(qū)域位置V1、漂浮式風電場址P1、靠近海峽內(nèi)部的P3、P4點，如圖 1所示。

3.2? ?模型設(shè)置

研究海域地形變化復(fù)雜，采用北南海大范圍區(qū)域和近岸區(qū)域兩套網(wǎng)格進行波浪數(shù)值模擬；圖 2（a）為大范圍區(qū)域網(wǎng)格劃分示意圖，最大網(wǎng)格尺度為1.5°，最小為0.3′×0.3′；圖 2（b）為近岸區(qū)域網(wǎng)格劃分示意圖，最大網(wǎng)格尺度為5 km，最小為300 m；海域水深采用實際測量水深和MIKE模型中C-MAP地形數(shù)據(jù)庫水深插值而成。計算過程中，波浪破碎系數(shù)取為0.9，底摩阻系數(shù)Kn取為0.01 m，白帽耗散系數(shù)Cdis和DELTAdis分別取值為4.5和0.5。

3.3? ?模型驗證

采用研究海域V1驗證點的實測數(shù)據(jù)對計算結(jié)果進行驗證，如圖 3所示：有效波高計算值與實測值平均誤差、均方根誤差、相關(guān)系數(shù)誤差分別為-0.14 m、0.22 m、0.94；有效波高最大值和平均值分別為3.42 m、1.40 m；實測有效波高最大值為3.10 m、1.27 m，計算值略大于實測值；譜峰周期計算值與實測值平均誤差、均方根誤差、相關(guān)系數(shù)，分別為-0.8 s、1.99 s、0.73；譜峰周期最大值和平均值分別為15.9 s和8.2 s；實測譜峰周期最大值為15.6 s、7.4 s ；波浪常浪向為E向，與研究海域的常浪向一致。

由驗證結(jié)果可以得到：所建立的波浪模型計算結(jié)果與實測波浪過程吻合良好，滿足研究海域波浪要素后報的精度要求。

4? ? 計算結(jié)果分析

4.1? ?有效波高的年平均值和最大值

采用上述波浪數(shù)值模型，模擬了研究海域2010～2020年間的逐時波浪過程計算結(jié)果。如圖 4 所示：代表點P1的逐時最大有效波高為4.13 m，平均有效波高為1.02 m；P2點的逐時最大有效波高為7.97 m；從P2到P1降低了接近50%；平均有效波高的分布基本符合折射定律，其等值線與等深線基本一致；最大有效波高的分布特點與平均有效波高存在明顯區(qū)別，最大波高出現(xiàn)在臺風“浪卡”過境期間，故波向受臺風風向影響顯著。

4.2? ?波浪受地形影響大

本文研究海域波況受地形影響大。從圖 4可以看出：波高從30 m等深線開始就迅速衰減，尤其是瓊州海峽受地形和海南島的遮掩，越往海峽內(nèi)波高越?。籔1點東側(cè)的地形丘壑相間，海底坡度較大，水深在10～30 m間。從外海傳過來的波浪經(jīng)過此區(qū)域波能發(fā)生快速衰減，盡管P1位置水深較大，但受地形影響，該位置以局部風成浪為主，波高顯著小于外海；另外，在V1點南側(cè)，有成片的有效波高變化較大的區(qū)域，大浪區(qū)正好對應(yīng)于水深較小的區(qū)域，這說明在此地形起伏較大的位置，局部的淺化效應(yīng)引起的波高增大不可忽略，在海上風電場設(shè)計和布置時必須考慮。

通過對比典型位置的風、浪玫瑰圖，進一步分析地形的影響。圖 5為P1、P2、V1點的風、浪玫瑰圖。從圖中可以看出，地形對波浪的影響體現(xiàn)在兩個方面，一是波高的顯著減小，二是浪向的顯著變化：

1）從P2到V1再到P1，水深依次減小，三點的風、浪玫瑰圖形狀接近，風向和風級組成變化不大；

2）三點的多年平均風速分別為6.60 m/s、6.07 m/s和5.61 m/s，V1、P1點的風速分別降低了7.9%、14.9%；

3）三點的多年平均有效波高分別為1.74 m、1.35 m、1.05 m，V1、P1點的有效波高分別降低22.6%、40.0%，波高的降幅顯著大于風速的降幅，這是由于地形快速變化，底摩擦和波浪破碎等引起的波能耗散；

4）在波向方面，對于外海的P2點浪向分布與風向分布較為一致，當?shù)竭_水深較淺的V1點處時，由于波浪的折射作用E向浪的比重顯著增加；當波浪傳播到P1點時，波浪進一步受到水下通道的影響，波向更為集中，E向浪占比高達75%，完全占據(jù)主導(dǎo)地位。

4.2? ?風和浪的變化趨勢

隨著全球氣候的變化，極端天氣的發(fā)生頻次和強度都逐漸增強，因此有必要研究海域的風浪變化趨勢，做好長遠規(guī)劃。統(tǒng)計各代表點2010-2020年的年平均風速和年平均有效波高，如圖 6和圖 7所示?？梢钥闯觯瑥耐夂５胶{內(nèi)部風速逐漸減小，波高也逐漸減小，風、浪的年際變化均較為平穩(wěn)，沒有明顯增大趨勢，越往海峽內(nèi)部，波浪的年際變化越小。

對比風和浪的年際變化可以發(fā)現(xiàn)，在靠近外海的P1、V1點，浪高和風速基本呈正相關(guān)關(guān)系，但到了海峽內(nèi)部的P1、P3、P4點，兩者的相關(guān)關(guān)系變?nèi)?，?018～2019年期間，甚至出現(xiàn)了風速減小、波高增大的情況，這可能是因為受到海峽地形的影響。

4.3? ?波高變化特征

整個南海海域的波浪存在顯著的季節(jié)性，圖8為外海P2點的波高隨月份的變化圖，可以看出：5～8月份波高較小，10～1月份波高較大。

圖9為P1、P2和V1點的月平均有效波高隨月份的變化圖，可以看出：各點平均有效波高隨月份的變化規(guī)律完全一致，外海波高的變化（誤差棒的長度）更為顯著。

圖10為各代表點平均波浪周期隨月份變化圖，可以看出：周期與波高變化的規(guī)律一致，從外海到海峽內(nèi)部周期逐漸減小，冬季平均周期顯著大于夏季。

P1位置位于海峽入口，受地形影響顯著，又是漂浮式海上風電場的規(guī)劃場址。統(tǒng)計表明：P1點多年平均H1/3為1.05 m，最大為2.30 m；多年平均Hmax為2.05 m，最大為4.46 m；多年平均的Tmean為4.92 s；平均波高的季節(jié)性變化特征顯著，月均有效波高和月均最大波高均從10月份開始顯著增加，至來年一月份波高又逐漸減小，說明該區(qū)域冬季波高較大；最大波高的變化規(guī)律和平均有效波高完全一致，月最大有效波高和月最大波高在10～1月最大。

5? ? 結(jié)論

針對廣東省海上風電的重點發(fā)展區(qū)域湛江，基于MIKE21-SW建立了波浪數(shù)學模型，模擬了研究海域2010～2020年的長時序逐時波浪過程，分析研究海域的波浪參數(shù)分布特征，統(tǒng)計研究海域多年波浪要素特征值，主要成果如下：

1）得到了多年有效波高平均值和極值的分布地圖，由于受地形影響，波浪能量從30 m等深線開始發(fā)生快速衰減；在規(guī)劃電場的P1和V1時位置，P1點大浪過程以局部風成浪為主，在地形變化較大的V1點南側(cè)，淺化效應(yīng)導(dǎo)致平均有效波高分布極不均勻，形成成片的大浪區(qū)，在風電場設(shè)計和布置時需加以考慮；地形還會影響波向，受折射作用和水下通道的綜合影響，波向越往海峽內(nèi)部越集中，到P1位置是E向浪接近75%；

2）2010～2020年間波浪沒有明顯增大趨勢，波高較為平穩(wěn)，越往海峽內(nèi)部波浪的年際變化越??；外海風速和有效波高正相關(guān)，海峽內(nèi)部兩者相關(guān)關(guān)系變?nèi)酰?/p>

3）有效波高呈現(xiàn)顯著的季節(jié)性，冬季（10～1月）有效波高顯著大于夏季（5～8月）；周期的季節(jié)變化和波高一致，冬季周期顯著大于夏季周期。

參考文獻

[1]李志川， 胡鵬， 馬佳星， et al. 中國海上風電發(fā)展現(xiàn)狀分析及展望 [J]. 中國海?上油氣， 2022， 34（05）.

[2]秦天慶， 袁長豐， 袁銘徽， et al. 基于2010—2020年東海風浪統(tǒng)計數(shù)據(jù)的海上風電吸力桶基礎(chǔ)動荷載承載特性 [J]. 海洋工程， 2023， 41（01）.

[3]楊斌， 楊萬康， 謝華偉， et al. 杭州灣北部實測波浪特征 [J]. 海洋與湖沼， 2022，?53（04）.

[4]夏瑞彬， 路超越， 梁楚進， et al. 基于WW3的南海北部有效波高時空變化及其極值重現(xiàn)期估算方法分析 [J]. 海洋學報， 2023， 45（02）.

[5]耿浩博， 張洪生， 洪楊彬， et al. 長江口海域風浪場的計算與特征分析 [J]. 水利水運工程學報， 2020， （05）.

[6]劉功鵬， 孫曉娟. 江蘇如東海域波浪數(shù)值模擬研究 [J]. 水電與新能源， 2020，?34（06）.

[7]韓樹宗， 董楊楊， 張水平， et al. 南海波浪時空變化特征研究 [J]. 海洋湖沼通報，?2020， （02）.

[8] KIM T R， LEE J H. Comparison of high wave hindcasts during typhoon Bolaven （1215）?using SWAN and WAVEWATCH III Model [J]. Journal of Coastal Research， 2018， （85?（10085））.

[9] DHI. Mike21 spectral wave module scientific documentation [M].Denmark： DHI Water?& Environment， 2009.

[10]吉會峰， 劉吉堂， 宋心剛， et al. 基于ERA5數(shù)據(jù)的江蘇海域風能資源評估 [J].?太陽能學報， 2023， 44（01）.

[11]譚海燕， 邵珠曉， 梁丙臣， et al. ERA5風場與NCEP風場在黃海、東海波浪模擬的適用性對比研究 [J]. 海洋通報， 2021， 40（05）.