楊振勇，呂迎雪，孔叢穎，李茂元

（1.中交第一航務工程局有限公司，天津300461；2.中交天津港灣工程研究院有限公司，天津300222）

海上工程施工期波浪（涌浪）預報系統(tǒng)

楊振勇1，呂迎雪2，孔叢穎2，李茂元1

（1.中交第一航務工程局有限公司，天津300461；2.中交天津港灣工程研究院有限公司，天津300222）

針對海上施工涌浪問題開發(fā)了波浪（涌浪）預報系統(tǒng)GWFS。首先針對涌浪的特點，包括涌浪的形成與傳播、涌浪周期的特點以及涌浪在全球的分布進行了分析。然后，對涌浪預報程序、風場數(shù)據(jù)、外海及近岸波浪模型和模型的驗證進行了詳細說明。模型對太平洋兩個實測數(shù)據(jù)點的波高和周期進行了詳細對比驗證，結(jié)果表明預報結(jié)果可靠。該系統(tǒng)的應用可有效提高海上工程可作業(yè)時間的利用率，避免或減少突變海況造成的損失。

海上工程；涌浪；波浪預報系統(tǒng)

0 引言

對于海上工程，包括海洋工程、海岸及港口工程，可靠的波浪數(shù)據(jù)是工程設計、施工以及項目運營的重要技術(shù)保障。波浪的推算可分為后報與預報，其中，波浪后報首先復演以往多年的波浪數(shù)據(jù)，對波浪各種分布進行分析，推算不同重現(xiàn)期極值波高等，其結(jié)果多用于設計階段，基于可靠的波浪數(shù)據(jù)進行的設計可保證工程抵御設計年限內(nèi)可能出現(xiàn)的各種極端海況。而波浪預報，則可使海上工程在施工期或運營期間事先預知未來幾天的海況變化，據(jù)此可合理地安排作業(yè)。對于海上項目施工，波浪預報使得船舶作業(yè)的安排可有效地利用小波高海況窗口期，避免可能出現(xiàn)較大波浪或長周期波浪對海上作業(yè)的影響。因此，波浪預報可提高海上有限可作業(yè)時間的利用率，并可避免或減少突變海況造成的損失。

對于風暴天氣產(chǎn)生的波浪，比如臺風波浪，由于目前的氣象預報對此可以預報，海上作業(yè)隊伍很容易事先獲知，因此這類風暴波浪不是本研究關(guān)注的重點。本文所介紹的波浪預報系統(tǒng)可預報各類波浪變化，但最突出的是對施工隊伍比較關(guān)注的施工期涌浪海況的預報。該系統(tǒng)是在作者多年波浪后報研究工作基礎上開發(fā)的，所采用的波浪數(shù)學模型WAVEWATCHIII[1]已在國際海洋波浪研究領域得到廣泛應用。由于波浪預報系統(tǒng)的驗證經(jīng)常借用后報研究中的一些內(nèi)容，本預報系統(tǒng)在世界各地的應用借鑒了不同海區(qū)波浪后報的許多經(jīng)驗。

由于涌浪具有和風浪顯著不同的特點，對其特點的理解非常有助于預報模型的建立和結(jié)果分析。因此，在介紹波浪預報系統(tǒng)之前首先對涌浪的特點做簡單分析。

1 涌浪的特性

1.1 涌浪的起源與傳播

涌浪是指離開風暴中心區(qū)域繼續(xù)傳播的波浪。在風暴中心區(qū)域波浪具有典型的風浪特征，其波高大小、周期及方向的分布相當分散。當波浪傳播離開風暴中心之后，波浪將不再從風中接受能量，不同的波浪成分將依據(jù)自身的特性傳播。周期較長者，傳播速度較快，在大洋環(huán)境中衰減相當緩慢。而周期較小者，傳播速度較慢。隨著距離的不斷增加，長波與短波逐漸分離，波浪傳播變化的這一過程即為頻率分散。具體表現(xiàn)為在某一海岸地區(qū)，某一時段內(nèi)所出現(xiàn)的涌浪周期相差不多。相似地，波浪傳播也存在方向分散過程，即在海岸地區(qū)觀察到某一時段的涌浪基本來自一個方向。由于涌浪的上述特點，當采用數(shù)學模型模擬涌浪的產(chǎn)生和傳播變形時，模型的范圍需要充分考慮原始風暴的區(qū)域，比如，涌浪可能是3 000 km之外的風暴產(chǎn)生的，那么模型范圍只考慮了工程附近1 000 km范圍的區(qū)域就不可能得到可靠的計算結(jié)果，這就是為什么涌浪研究一般要考慮一個大洋的范圍。

1.2 涌浪的周期分布

對于影響某一工程地點的涌浪，其特性很大程度上取決于原始風暴。除原始風暴的強度之外，原始風暴與工程地點的距離等因素對抵達工程地點的涌浪影響很大。因此，拋開原始風暴而孤立地分析某一工程地點的涌浪特點有可能導致不準確的結(jié)論。由于涌浪起源的這些特點，一般涌浪的周期分布并不像風浪那樣與波高的相關(guān)性很強，相反地，涌浪波高越大相對應的周期可能變小，主要原因是產(chǎn)生大波高的風暴與工程地點距離一般較近一些。圖1為非洲西海岸剛果黑角港（Pointe Noire）實測波浪周期與波高的聯(lián)合分布，該海區(qū)波浪屬于典型的涌浪。圖2則是北美安大略湖實測波浪周期與波高分布?？梢钥闯霾ㄖ芷诤筒ǜ呦嚓P(guān)性很好，這里屬于典型的風浪特征。

圖1 剛果黑角外海實測波高Hs與周期Tp分布（涌浪為主）Fig.1Measured wave height Hsand wave period TpinPointe Noire,Congo(mainly in swell)

圖2 安大略湖實測波高Hs與周期Tp分布（風浪）Fig.2Measured wave height Hsand wave period Tpin Lake Ontario(wind wave)

1.3 全球各地涌浪的分布及影響區(qū)域

從全球范圍來看，經(jīng)常出現(xiàn)大浪頻率較高的地區(qū)主要是南北半球中高緯度的西風帶（35°~ 60°）。西風帶，又稱暴風圈，該區(qū)域極易形成溫帶氣旋，常常是一個氣旋未完，另一個氣旋生成。在此區(qū)域，西風一般比較強勁持久，因而在該海洋區(qū)域經(jīng)常出現(xiàn)大浪。一般地，南半球波浪更強，且冬季大浪幾率明顯高于夏季。在北半球的冬季，位于北半球的西風帶頻繁出現(xiàn)大風或風暴。在南半球的冬季，南半球西風帶風暴頻繁風速強勁。這就是全球很多海岸地區(qū)經(jīng)常或出現(xiàn)季節(jié)性涌浪的原因，其中包括我國參與建設的一些港口海岸工程項目所在國家和地區(qū)，如印尼印度洋海岸、斯里蘭卡、阿拉伯海、非洲大西洋海岸、南美太平洋海岸等。我國海岸的涌浪主要由季風和西太平洋的風暴（臺風）造成，與上述國家和地區(qū)相比，我國海岸涌浪的頻率和強度要小許多。

2 波浪（涌浪）預報系統(tǒng)GWFS

2.1 預報模型系統(tǒng)GWFS

波浪預報系統(tǒng)GWFS由外海波浪模型、近岸波浪模型以及多個數(shù)據(jù)處理程序組成。對于不同大洋分別建立了不同的模型。該模型每日定時自動啟動，每日計算1次（預報1次），每次預報未來1周的波況，預報數(shù)據(jù)間隔為1 h，輸出結(jié)果包括有效波高、峰值周期和平均波向等數(shù)據(jù)。輸出結(jié)果包括波浪參數(shù)變化圖和結(jié)果文本文件，根據(jù)需要可自動發(fā)至網(wǎng)上或通過郵件發(fā)送給有關(guān)人員。

圖3為波浪預報系統(tǒng)流程圖。

圖3 波浪預報系統(tǒng)（程序）流程圖Fig.3Flow-chart of the wave forecasting system

2.2 風場

預報模型是通過輸入風場變化數(shù)據(jù)計算波浪。本系統(tǒng)輸入風場采用美國大氣與海洋研究署環(huán)境預報中心（NOAA/NCEP）全球風場預報數(shù)據(jù)GFS（也可采用歐洲中期天氣預報中心ECMWF的數(shù)據(jù)）。該數(shù)據(jù)可較好地預報未來6~7 d的全球風場變化。NOAA/NCEP同時提供超過7 d以上的天氣預報數(shù)據(jù)，根據(jù)經(jīng)驗，超過7 d的海洋天氣預報數(shù)據(jù)準確性較低，因此，本模型僅給出未來7 d的結(jié)果。風場數(shù)據(jù)間隔為1 h，空間網(wǎng)格為0.5°。經(jīng)試算，模型可以在各大洋不同海岸給出可靠的波浪推算結(jié)果。

2.3 外海波浪模型

波浪數(shù)學模型采用Tolman的WAVEWATCHIII（Tolman，2009）。該模型可較好地模擬風浪生成及衰減過程、涌浪的傳播、波浪的非線性相互作用、波浪的衰減等。該模型經(jīng)海洋波浪工作者大量檢驗和應用，現(xiàn)已成為海洋波浪模型最可靠的模型之一，特別適合大洋波浪包括涌浪的傳播變形的計算[2-4]。該模型所依據(jù)的控制方程可表示如下：

式中：N（k，θ；x，t）為波浪作用譜（wave action density spectrum）；k、θ為波數(shù)及方向；S為源項，包括風影響、波浪非線性作用、波能的各種損耗等。

2.4 近岸波浪傳播模型

外海模型范圍很大，因此一般模型網(wǎng)格較大（0.3°~1.0°）。此外，由于近岸區(qū)域地形和岸線走向的變化，近岸波浪的傳播變形，比如港口附近，甚至港內(nèi)水域，一般采用近岸模型對波浪自外海至近岸的傳播變形進行推算。

在本文介紹的預報系統(tǒng)中，近岸波浪的傳播變形采用MIKE21SW模型或SWAN模型[5]。對于某一項目地點，首先采用該模型對不同的波浪分布（波高周期波向水位的組合）進行計算，從而可得到各種波浪的傳播變形系數(shù)。之后，結(jié)合外海波浪數(shù)據(jù)即可預報近岸甚至港內(nèi)區(qū)域的波況變化。

2.5 預報模型驗證

由于實測數(shù)據(jù)有限，本模型采用太平洋北美海岸部分測波儀實測數(shù)據(jù)對預報結(jié)果進行了驗證，此外，對該模型在安哥拉項目的測試進行了簡單介紹。需要說明的是，以下模型驗證中各點的波浪均以涌浪為主。

驗證實例1：太平洋測波儀數(shù)據(jù)驗證。

圖4所示兩個測波儀分別為美國大氣與海洋管理署NOAA下屬浮標數(shù)據(jù)中心（NDBC）所有（測站編號分別為32012，46006）。圖5為2016年10月19—29日預報數(shù)據(jù)和實測數(shù)據(jù)對比結(jié)果。對比結(jié)果表明，模型預報在6~7 d以內(nèi)的結(jié)果可靠程度很高，而超過1周的結(jié)果則有很大的不確定性，這主要是受海洋天氣預報數(shù)據(jù)準確程度的影響。換句話說，波浪模型結(jié)果可靠與否，取決于海洋區(qū)域天氣預報結(jié)果的準確程度。因此，目前海洋天氣預報水平基本上在6~7 d以內(nèi)比較可靠，超過此期限可靠程度明顯下降。當然，對于陸地天氣預報，由于地形影響明顯，陸地地形變化較大的區(qū)域準確程度則更低一些。

圖4 太平洋波浪預報模型實測數(shù)據(jù)測波儀位置示意圖Fig.4Location of wave buoys used in the Pacific Ocean wave forecasting model validation

46006測站

32012測站

需要說明，實測與預報波浪的周期均為峰值周期。由于大洋中的波浪譜經(jīng)常為多峰波譜，不同譜峰的變化可造成最大峰值周期的跳躍性變化。

圖5 太平洋不同位置波浪預報與實測波高對比結(jié)果Fig.5Comparison between measured wave height and forecasted results at different locations in the Pacific Ocean

驗證實例2：安哥拉卡賓達海岸波浪預報。

安格拉卡賓達海岸波浪以來自南大西洋的涌浪為主，時間上和南大西洋西風帶風暴強度變化一致，即南半球的冬季5—10月期間波浪較大，而11—翌年4月波浪偏小。本項目位于南緯6°，目前尚無波浪測量設施，但自2016年7月以來每日對岸邊波浪進行目測，表明預報結(jié)果相當可靠。圖6為2016年8—9月期間部分預報結(jié)果?，F(xiàn)場觀測顯示，模型對9月3日之前的小浪期和9月3—4日出現(xiàn)的大浪期預報與現(xiàn)場觀測十分吻合。

圖6 安哥拉卡賓達海岸波浪預報結(jié)果（有效波高與峰值周期）Fig.6Forecasted wave height and period in Cabinda,Angola(effective wave height and peak period)

3 結(jié)語

本文首先對涌浪的一些特點進行了分析，包括涌浪的起源與傳播、涌浪周期的特點以及全球不同海域的涌浪分布等。之后對作者開發(fā)的波浪（涌浪）預報系統(tǒng)GWFS進行了詳細介紹。該系統(tǒng)可定時自動啟動，自動下載風場數(shù)據(jù)然后進行計算和結(jié)果處理，計算結(jié)束后自動發(fā)布預報結(jié)果。模型驗證和現(xiàn)場測試結(jié)果表明，該系統(tǒng)可較好地預報未來6~7 d的海況，可在全球各大洋和海域應用。本預報的應用將有助于提高海上項目有限可作業(yè)時間的利用率，避免或減少突變海況造成的損失。

[1]TOLMAN H.User manual and system documentation of WAVEWATCH III version 3.14[Z].NOAA/NWS/NCEP/MMAB Technical Note,2009.

[2]HANSON J,TRACY B,TOLMAN H.Pacific hindcast performance evaluation of three numerical wave models[J].Journal of Atmospheric and Oceanic Technology,2009,26(26):1 614-1 633.

[3]LEE B,FAN Y,CHEN J.Verification on swell simulation of WAVEWATCH III wave modeling[C]//International conference on computational science&engineering,2011:65-74.

[4]PADILLA-HERNANDEZ R,SCOTT D.Inter-comparison of modern operational wave models[C]//8th international workshop on wave hindcasting and forecasting.Hawaii,2004:14-19.

[5]DHI.MIKE21 SW user manual[Z].Denmark:DHIGROUP,2012.

A wave(swell)forecasting system for marine works

YANG Zhen-yong1,Lü Ying-xue2,KONG Cong-ying2,LI Mao-yuan1
（1.CCCC First Harbour Engineering Co.,Ltd.,Tianjin 300461,China; 2.CCCC Tianjin Port Engineering Institute Co.,Ltd.,Tianjin 300222,China）

Based on the swell problems during marine construction,we developed a wave（swell)forecasting system.First, swell properties,including swell generation and propagation,characteristics of swell period and global distribution of swells have been reviewed.Next,the wave forecasting system,input wind data,ocean wave model and nearshore wave model of the system,and detailed verification of the model system have been discussed.Detailed comparisons have been conducted to verify the model using two measured wave buoy datasets,and the result shows that the forecasting model is reliable.With the wave forecasting system,marine works can be reasonably scheduled,as a result,the working efficiency can be improved,and some unpredictable downtime or damage could be avoided.

marine works;swell;wave forecasting system

U655.11；P731.22

A

2095-7874（2017）05-0006-04

10.7640/zggwjs201705002

2016-09-28

2017-02-24

楊振勇（1963—），男，河北欒城人，高級工程師，主要從事海岸動力研究及國際工程管理。

E-mail：yangzhenyong2013@163.com