江蘇沿海波浪多重嵌套模擬研究

儲錫君，徐福敏

(河海大學 海岸災害及防護教育部重點實驗室，港口海岸與近海工程學院，江蘇 南京 210098)

江蘇沿海波浪多重嵌套模擬研究

儲錫君，徐福敏

(河海大學 海岸災害及防護教育部重點實驗室，港口海岸與近海工程學院，江蘇 南京 210098)

根據24年CCMP風場資料和江蘇沿海4個方向(N、NE、E和SE)百年一遇風速，構建西北太平洋、東中國海和江蘇沿海上述4個方向的百年一遇風場。首次建立一個基于第三代海浪模型SWAN的自西北太平洋、東中國海至江蘇沿海的三重嵌套數值模型，以AVISO衛星觀測數據和江蘇沿海定點實測數據進行驗證。以三個計算域4個方向百年一遇風場為驅動風場，驅動該多重嵌套模型，高精度數值模擬江蘇沿海4個方向百年一遇有效波高分布并進行分析。結果表明，江蘇沿海輻射沙洲地形對有效波高分布影響顯著；E向百年一遇風場作用下海域有效波高最大，NE向次之，N向和SE較小。

CCMP風場；SWAN模型；多重嵌套模擬；有效波高分布；江蘇沿海

江蘇地理位置優越、經濟發展迅速，綜合實力在全國名列前茅，其海岸線北起蘇魯交界的繡針河口，南抵長江口，全長約888.95 km[1]，沿海地區自然資源豐富。同時其地理位置和復雜自然條件，沿海地區常受到臺風和寒潮的影響，造成重大經濟損失，嚴重威脅人員生命財產安全。確定江蘇沿海不同重現期的波浪要素，合理確立海岸工程使用期間可能出現的極值波高，對江蘇沿海地區的安全、海岸工程建筑物的設計和建設、防災減災具有極其重要作用。

江蘇沿海的大浪主要是由極端天氣引起的，由于海域缺乏長期連續的波浪實測資料，本研究首次以海浪數值模型多重嵌套模擬方法得到該海域百年一遇波浪要素。研究首次建立一個基于海浪模型SWAN的自西北太平洋、東中國海至江蘇沿海的三重嵌套模擬模型并以AVISO衛星觀測數據和江蘇沿海響水站實測數據進行驗證。以共計24年(1988年1月至2011年12月)的CCMP(Cross Calibrated，Multi Platform)風場數據按泊松-耿貝爾復合極值分布推算西北太平洋和東中國海4個方向(N、NE、E和SE)的百年一遇風場，以江蘇沿海上述4個方向的百年一遇風速構建江蘇沿海的百年一遇風場。以三個海域4個方向的百年一遇風場驅動該多重嵌套數值模型模擬江蘇沿海的上述4個方向百年一遇有效波高分布，并對有效波高分布特征進行分析。

1 西北太平洋至江蘇沿海的百年一遇風場分布特征

目前常用的風場有NCEP風場、QN(QuickSCAT/NCEP)混合風場和CCMP風場等，其中NCEP風場的空間分辨率較低(1.875°×1.905°)，QN混合風場具有較高空間分辨率(0.5°×0.5°)，但時間序列較短(1999年10月至2009年7月)，CCMP風場相比則具有分辨率高、時間序列長等優點[2]。CCMP風場資料來自美國航空航天局(NASA)，數據自1987年7月至2011年12月，覆蓋全球，空間分辨率為0.25°×0.25°，時間分辨率為6 h。

Atlas[3]通過實驗證實了CCMP較其他單個的衛星平臺測量的風場數據在精度方面有很大地提高，能夠滿足很多海洋和大氣環境應用研究的需要。張鵬等[4]利用渤海海區一個點位11天的風場數據，對CCMP的風場數據進行驗證，結果表明風場的風速和風向與相應的實測值之間的相關系數分別為0.89和0.91。故CCMP風場數據能較好的代表實際風場狀況。

本研究使用24年(1988年1月至2011年12月)的CCMP風場數據，根據泊松-耿貝爾復合極值分布方法推算西北太平洋至東中國海的4個方向(N、NE、E和SE)百年一遇風場。

1.1泊松-耿貝爾復合極值分布

目前常用的分布函數有皮爾遜Ⅲ型累計頻率曲線、對數正態概率分布、第一極值分布(耿貝爾分布)和泊松-耿貝爾復合極值分布等。在統計我國沿海臺風引起的大浪，常采用泊松-耿貝爾復合極值分布，因為該分布不僅考慮了年極值分布，還考慮了年臺風出現的次數，因而比較合理[5]。

由于每年的臺風路徑、強度、頻率都是隨機的，構成離散型分布，而TC影響下的最大風速構成連續型分布。由于每年臺風出現的頻次n符合泊松分布，假設臺風影響下海域風速服從耿貝爾分布，可推得泊松-耿貝爾復合極值分布的分布函數，并可進一步推算得概率為P的大風極值VP為：

1.2西北太平洋和東中國海的百年一遇風場

基于自1988年1月至2011年12月共24年的CCMP風場數據，計算西北太平洋和東中國海每隔6小時的風場，并劃分為8個方位(N、NE、E、SE、S、SW、W和NW)。由于本研究主要考慮百年一遇風場作用下江蘇沿海大浪，只考慮西北太平洋和東中國海范圍內對江蘇沿海影響較大的4個方向(N、E、NE和SE)。

由中國臺風網資料統計表明西北太平洋上述24年間共計有熱帶氣旋692個，根據泊松-耿貝爾復合極值分布，推算得到西北太平洋和東中國海計算域的上述4個方向百年一遇風速分布情況，圖1為西北太平洋和江蘇沿海E向百年一遇風速等值線分布。

可見西北太平洋百年一遇風速的最大值可達65 m/s，位于北緯40°～50°、東經152°～165°范圍內；西北太平洋高緯海域百年一遇風速比低緯普遍偏大，靠近赤道時則不足20 m/s；東中國海百年一遇風速的最大值可達50 m/s，位于外海北緯23.3°～28.5°、東經126°～128°范圍內，百年一遇風速從外海向海岸逐漸減小。計算研究發現其余3個方向的西北太平洋和東中國海的百年一遇風場具有與E向相似的分布特征。

圖1 E向百年一遇風速等值線分布(單位：m/s)Fig. 1 The contour distribution of East wind speed with a 100-year return period

1.3江蘇沿海的百年一遇風場

江蘇沿海范圍內CCMP風場網格數據較少，故不采用CCMP風場數據推算江蘇沿海的百年一遇風場。本研究使用贛榆、西連島、燕尾島、射陽、大豐、如東、呂四和啟東8個測站1981年至1998年風速實測資料推算出的江蘇沿海N、NE、E和SE向百年一遇風速和西連島1974年至2003年風速實測資料推算出的江蘇沿海NE和E向百年一遇風速。表1為基于上述資料推算得到的江蘇沿海百年一遇風速表。

可見基于1981—1998年風速實測資料推求得到的江蘇沿海4個方向海上百年一遇風速中，SE向比其它三向明顯偏小，E向最大，接近40 m/s，NE向其次；基于1974—2003年風速實測資料推求得到的NE向和E向百年一遇風速均大于40 m/s，E向略大些。

表1 江蘇沿海百年一遇風速表Tab. 1 Wind speed along Jiangsu Coast with a 100-year return period

2 SWAN模型簡介

海浪數值預報始于20世紀50年代[6]，數值模擬發展到20世紀末已達到比較成熟的階段，數值模式被廣泛應用于模擬海浪的發展過程和規律，形成了以WAM模型[7]、WAVEWATCH模型[8]和SWAN(simulating waves nearshore)模型[9]為代表的第三代海浪數值模型。Booij等以波浪線性理論及荷蘭、德國的一些現場觀測資料進行驗證，表明SWAN模型對復雜地形下的波浪模擬能力好。Zijlema等[10]就SWAN模型收斂性及數值計算的準確性作了分析及改進。Gorman等[11]利用SWAN模型模擬了潮間帶河口的波浪生成及傳播過程。Rogers和Lin等的研究表明，SWAN模型能很好地模擬近岸涌浪、風浪和混合浪。

SWAN海浪模型在國內也得到了廣泛研究應用。陳希等[12]采用SWAN模型嵌套方法，對一次臺風浪過程進行了模擬。徐福敏等[13-15]采用SWAN模型對不同水深時流場對于波浪傳播的影響進行了探討，并將模型應用于海安灣的淺水波浪數值計算，模擬灣內波高及周期場分布，與觀測結果對比驗證表明符合良好。張洪生等[16]基于CCMP衛星遙感風場數據，通過將WAVEWATCH和SWAN模型嵌套的方法，數值模擬了珠江口附近海域的風浪場，定量分析表明數值模型模擬結果良好。

SWAN是Delft大學開發的基于動譜平衡方程的第三代譜波浪數值模型。動譜密度N(σ,θ)與能譜密度E(σ,θ)的關系為：N(σ,θ)=E(σ,θ)/σ，σ為相對波頻，θ為波向。在球坐標系下，動譜平衡方程可表示為：

SWAN模型采用全隱式有限差分格式，無條件的穩定[17]，具有計算精度高的特點，發展至今已成功應用于河口、近海水域和大尺度海域的海浪預報[13,15]。

3 西北太平洋至江蘇沿海三重嵌套模型驗證

本文建立的三重嵌套模型區域自外向內分別：西北太平洋計算域(100°E～180°E， 0°N～50°N)，如圖1(a)所示；東中國海計算域(116.9°E～130.8°E， 23.3°N～40.8°N)，如圖1(b)所示；江蘇沿海計算域(119°10' E～122°24' E，31°40' N～35°6' N)，如圖2所示。外兩層計算域水深數據是將從NOAA(National Oceanic and Atmospheric Administration)下載的1＇×1＇水深數據插值處理得到，西北太平洋計算域水深插值后的分辨率為0.25°×0.25°，東中國海計算域插值后的分辨率為4＇×4＇。江蘇沿海水深地形條件復雜，近海存在輻射沙洲群，海域使用海圖水深，分辨率為1＇×1＇。

圖2 江蘇沿海的水深分布(單位：m)Fig. 2 The depth distribution of Jiangsu coast

西北太平洋計算域模型計算網格的空間分辨率為0.25°×0.25°；東中國海模型計算網格分辨率為0.1°×0.1°；江蘇沿海模型計算網格經度方向上的分辨率為0.776＇，緯度方向上為0.687＇。外層計算域為內層提供波浪邊界條件。模型采用球形坐標系，計算頻率為0.04～1.0 Hz，以fi+1=λfi劃分為20個頻域，λ為常數；譜方向網格數為24個；底摩擦造成的能量損耗采用Collins非線性表達式，系數Cf=0.006；破波參數γ=0.73；其余物理過程如白浪耗散、非線性相互作用等參數采用模型默認值，模型時間步長為30分鐘。

本研究分別對西北太平洋至東中國海海浪模型和江蘇沿海海浪模型進行驗證。驗證時段自2011年5月15日0時至5月24日18時(UTC，世界標準時間)，驅動風場為該時段的CCMP風場。

3.1西北太平洋至東中國海海浪模型的驗證

AVISO(Arching Validation and Interpretation of Satellite Oceanographic data)是法國國家空間研究中心開發的衛星觀測數據，本研究使用其有效波高觀測數據。有效波高觀測數據的時間范圍為2009年9月至今，覆蓋全球，空間分辨率為1°×1°，時間分辨率為24小時。圖3(a)為SWAN模型數值模擬得到的西北太平洋2011年5月20日0時的有效波高分布；圖3(b)為該時刻的AVISO衛星觀測有效波高分布。

根據表7顯示：標準舞練習組女生在上肢快速反應測試中使用時間平均減少了7.35s，拉丁舞練習組女生平均減少了8.17s，兩組之間沒有顯著性差異。觸桿反應測試和十字繞桿跑測試中，標準舞練習組女生使用時間分別平均減少了3.06s和3.09s，拉丁舞練習組女生分別平均減少了2.88s和2.67s，兩組之間差異不顯著。全身協調測試中兩組之間差異顯著，因標準舞練習組女生比拉丁舞練習組女生前測成績相對較好，且無訓練基礎但協調能力較好的人群在經歷訓練后協調能力能夠得到提升，但相對協調能力較差的人群，其提高程度一般較低[19]。

圖3 西北太平洋5月20日0時有效波高分布(單位：m)Fig. 3 The distribution of significant wave height at Northwest Pacific at zero on May 20th

SWAN數值模擬結果與衛星觀測有效波高數據對比分析，西北太平洋的海域有效波高最大值均位于北緯35°～50°，東經160°～180°范圍內，等值線分布較為相似，大部分海域的有效波高相近。計算西北太平洋所有網格節點的SWAN模擬有效波高與衛星觀測數據的均方根誤差為0.44 m，表明對西北太平洋海浪模型驗證良好。計算研究發現東中國海SWAN數值模擬的有效波高等值線分布與衛星觀測相似，海域有效波高值相近。

3.2江蘇沿海海浪模型的驗證

響水站位于江蘇沿海響水縣34.437°N、120.1°E，如圖2所示。研究使用海圖水深插值得到響水站的水深為8.36 m，測站由河海大學維護，使用波浪浮標遙感系統獲取波浪要素。測站數據包括有效波高、平均波周期、波向和最大波高等，時間范圍自2010年11月1日至2012年7月3日，分辨率為1小時。研究對江蘇沿海海浪模型進行驗證時，對模型參數選取進行數值模擬研究和分析。將底摩擦系數Cf分別取值為0.006、0.010和0.015 ，SWAN三重嵌套模型模擬得到響水站在2011年5月15日0時至5月24日18時的有效波高，并進行敏感性分析。結果表明當Cf=0.006時，數值模擬有效波高最大值與實測最為接近，故模型底摩擦系數取為0.006。圖4(a)為SWAN三重嵌套模型模擬得到的響水站在2011年5月15日0時至5月24日18時的有效波高與測站實測數據的對比；圖4(b)為模型模擬得到的該時間范圍內的平均波周期與測站實測數據的對比。

可見SWAN數值模擬得到的響水站有效波高、平均波周期與實測值變化趨勢均相近，兩者的有效波高均方根誤差為0.15 m，平均波周期均方根誤差為1.42 s，表明結果吻合較好。

對西北太平洋至江蘇沿海的模型驗證表明SWAN模型三重嵌套進行數值模擬可較好得模擬計算海域的有效波高分布和海浪動態變化過程，證明了SWAN三重嵌套數值模擬模型的合理性。

4 多重嵌套方法模擬江蘇沿海百年一遇波浪

以第1節的三海域的百年一遇風場為驅動風場，基于上述西北太平洋至江蘇沿海SWAN三重嵌套數值模型，模擬得到江蘇沿海4個方向(N、NE、E和SE)百年一遇有效波高分布。

4.1江蘇沿海百年一遇波浪分布模擬結果

基于1974-2003年風速資料推算得到的江蘇沿海E、NE向百年一遇風速比基于1981—1998年風速資料的推算結果偏大(見表1)，模擬得到的江蘇沿海外海百年一遇有效波高偏大，近海輻射沙洲有效波高分布則相似。分別作圖5中(a)-(d)(基于1981—1998年風速資料，見表1)和(e)-(f)(基于1974—2003年風速資料，見表1)。

可見，江蘇沿海外海處有效波高等值線較為順滑，有和岸線平行的趨勢，在圖5(d)SE向百年一遇有效波高分布體現的最為明顯。隨著波浪向近岸傳播，由于底摩擦、淺化破碎等損耗，有效波高迅速減小，北部海域有效波高變化比較平緩，在南部輻射沙洲海域(32.5°～33.5°N)近岸百年一遇有效波高的變化較為復雜，對比江蘇沿海的水深分布(見圖2)，可明顯觀察到輻射沙洲處百年一遇有效波高偏小，有效波高分布與水深分布相關性很大。表明輻射沙洲對江蘇近岸有效波高分布的影響顯著。

圖5 江蘇沿海百年一遇有效波高分布Fig. 5 The distribution of significant wave height a with 100-year return period along Jiangsu coast

4.2江蘇沿海百年一遇有效波高分布特征

表2為SWAN多重嵌套模擬得到江蘇沿海4個方向(N、E、NE和SE)的百年一遇有效波高的最大值。

表2 江蘇沿海百年一遇有效波高的最大值Tab. 2 The maximum significant wave height with a 100-year return period along Jiangsu coast

表2的江蘇沿海各向百年一遇有效波高最大值與謝冬梅等[1]計算得到江蘇外海深水邊界點分方向統計百年一遇波高結果相比較為接近。由圖5和表2可見，基于1981—1998年風速資料模擬得到的4個方向(N、E、NE和SE)江蘇沿海百年一遇有效波高分布，外海N向和SE向的有效波高比其它向明顯偏小，NE向偏大，近海4個方向的有效波高分布相似；海域E向有效波高最大值最大，為11.09 m，NE向其次，N向最小。江蘇沿海N向百年一遇風速較SE明顯偏大，但兩向百年一遇風場影響下有效波高分布相似，最大值相近，可知若在相同風速條件下，SE向風作用下產生的江蘇沿海的有效波高比N向大，對江蘇沿海的影響更大。SE向與N向百年一遇風場作用下江蘇沿海百年一遇有效波高偏小，對江蘇沿海影響也最小，E向最大，NE向其次。

5 結 語

本研究首次建立一個基于海浪模型SWAN的自西北太平洋、東中國海至江蘇沿海的三重嵌套模型，經驗證模型的合理性，該模型對江蘇沿海重現期波浪的數值模擬提供了一個有效途徑。以西北太平洋、東中國海和江蘇沿海三個計算域4個方向百年一遇風場驅動三重嵌套模型模擬得到江蘇沿海百年一遇有效波高分布。研究表明，輻射沙洲對有效波高分布有顯著影響，E向風場作用下的江蘇沿海百年一遇有效波高分布最大，NE向次之，N向和SE較小。

本研究根據泊松-耿貝爾復合極值分布推算西北太平洋和東中國海4個方向百年一遇風速時只統計考慮了臺風影響總次數。下一步工作將綜合考慮臺風與寒潮的影響，推算西北太平洋至江蘇沿海4個方向不同重現期的風速和多重嵌套模擬江蘇沿海不同重現期的有效波高分布。

[1] 謝冬梅，陳永平，張長寬. 江蘇外海深水波要素推算[J]. 海洋工程，2013，31(2)：94-99. (XIE Dongmei, CHEN Yongping, ZHANG Changkuan. Calculation of deepwater waves along Jiangsu coast[J]. The Ocean Engineering, 2013, 31(2):94-99. (in Chinese))

[2] 鄭崇偉. 基于CCMP風場的近22年中國海海表風場特征分析[J]. 氣象與減災研究，2011，34(3)：41-46. (ZHENG Congwei. Sea surface wind field analysis the China Sea during the last 22 years with CCMP wind field[J]. Meteorology and Disaster Reduction Research, 2011, 24(3):41-46. (in Chinese))

[3] ATLAS R, JOSEPH A, HOFFMAN R N. Application of satellite surface wind data to ocean wind analysis[C]//Processing of the SPIE―The International Society for Optical Engineering. 2008.

[4] 張鵬，陳曉玲，陸建忠，等. 基于CCMP衛星遙感海面風場數據的渤海風浪模擬研究[J]. 海洋通報，2011，30(3)：266-271. (ZHANG Peng, CHEN Xiaoling, LU Jianzhong, et al. Research on wave simulation of Bohai Sea based on the CCMP remotely sensed sea winds[J]. Marine Science Bulletin, 2011, 30(3):266-271. (in Chinese))

[5] 邱大洪. 工程水文學[M]. 北京：人民交通出版社，2011. (QIU Dahong. Engineering hydrology[M]. Beijing: China Communications Press, 2011. (in Chinese))

[6] 趙紅軍，宋志堯，徐福敏，等. 南中國海臺風浪數值模擬研究—以臺風“珍珠”為例[J]. 海洋工程，2010，28(3)： 128-134. (ZHAO Hongjun, SONG Zhiyao, XU Fumin, et al. Numerical simulation of typhoon waves in the South China Sea―a case study of typhoon Chanchu[J]. The Ocean Engineering, 2010, 28(3):128-134. (in Chinese))

[7] WAMDI Group. The WAM model-a third generation ocean wave prediction model[J]. Journal of Physical Oceanography, 1988, 18:1775-1810.

[8] TOLMAN. User manual and system documentation of WAVEWATCH-Ⅲ version 2.22[M/OL]. Technical Note 222, NCEP/NOAA/NWS. Washington D C: National Center for Environment Prediction, 2002.

[9] BOOIJ N, RIS R C, HOLTHUIJSEN L H. A third-generation wave model for coastal regions, PartⅠ: Model description and validation[J]. Journal Geophysical Research, 1999, 104(C4):7649-7666.

[10] ZIJLEMA M, ANDRE J. van der Weathuysen. On convergence behavior and numerical accuracy in stationary SWAN simulations of nearshore wind wave spectra[J]. Coastal Engineering, 2005,52:237-256.

[11] GORMAN R M, NEILSON C G. Modeling shallow water wave and transformation in an intertidal estuary[J]. Coastal Engineering, 1999, 36:197-217.

[12] 陳希，閩錦忠，沙文鈺，等. 近岸海浪模式在中國東海臺風浪模擬中的應用—數值模擬及物理過程研究[J]. 海洋通報，2003，22(2)：9-16. (CHEN Xi, MIN Jinzhong, SHA Wenyu, et al. A simulation of nearshore typhoon waves[J]. Marine Science Bulletin, 2003, 22(2):9-16. (in Chinese))

[13] 徐福敏，張長寬，茅麗華，等. 一種淺水波浪數值模擬的應用研究[J]. 水動力學研究進展，2000，12(4)：429-434. (XU Fumin, ZHANG Changkuan, MAO Lihua, et al. Application of a numerical model for shallow water waves[J]. Journal of Hydrodynamics, 2000, 12(4):429-434. (in Chinese))

[14] XU Fumin, YAN Yixin, ZHANG Changkuan. Wave numerical model for shallow water[J]. China Ocean Engineering, 2002(2):193-202.

[15] CAVALERI L, RIZZOLI P M. Wind wave prediction in shallow water: Theory and applications[J]. Journal of Geophysical Research: Oceans(1978-2012), 1981, 86(C11):10961-10973.

[16] 張洪生，辜俊波，王海龍，等. 利用WAVEWATCH和SWAN嵌套計算珠江口附近海域的風浪場[J]. 熱帶海洋學報，2013，32(1)：8-17. (ZHANG Hongsheng, GU Junbo, WANG Hailong, et al. Simulating wind wave field near the Pearl River Estuary with SWAN nested in WAVEWATCH[J]. Journal of Tropical Oceanography, 2013, 32(1):8-17. (in Chinese))

[17] 李紹武，梁超，莊茜. SWAN風浪成長模型在近海設計波浪要素推算中的應用[J]. 港工技術，2012，49(2)：5-7. (LI Shaowu, LIANG Chao, ZHUANG Qian. Application of SWAN model in onshore design wave factors calculation[J]. Port Engineering Technology, 2012, 49(2):5-7. (in Chinese))

Study on multiple-nested simulation of waves along Jiangsu coast

CHU Xijun， XU Fumin

(Key Laboratory of Coastal Disasters and Defence, Ministry of Education, College of Harbor, Coastal and Offshore Engineering, Hohai University, Nanjing 210098, China)

Wind field with a 100-year return period in Northwest Pacific, East China Sea and Jiangsu Coast in four directions of North, Northeast, East and Southeast is built with CCMP wind data of 24 years and wind speed with a 100-year return period along Jiangsu coast in each direction as mentioned above. Based on the third-generation wave model SWAN, a multiple-nested simulation model is built in three levels of nested domains from Northwest Pacific, East China Sea to Jiangsu coast for the first time and verified with AVISO satellite data and measured data at fixed points along Jiangsu coast. Using wind field with a 100-year return period in these three domains in four directions as the driving wind field of this multiple-nested simulation model, the distribution of significant wave height with a 100-year return period along Jiangsu coast can be simulated with high accuracy. It is concluded that the radial sandbanks along Jiangsu coast affect the distribution of significant wave height obviously; significant wave height is generally the largest under the effect of East wind field with a 100-year return period, Northeast the second, North and Southeast the least.

CCMP wind field; SWAN model; multiple-nested simulation; distribution of significant wave height; Jiangsu coast

P731.22

A

10.16483/j.issn.1005-9865.2015.06.016

1005-9865(2015)06-112-07

2015-01-22

“十二五”國家科技支撐計劃(2012BAB03B01)

儲錫君(1990-)，男，浙江寧波人，碩士生，主要從事災害性海洋動力過程研究。E-mail:chu293546@126.com