梁連松，李瑞杰，2，豐青，孔從穎，宋曉波

（1.河海大學(xué)海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室，南京210098；2.河海大學(xué)環(huán)境海洋實(shí)驗(yàn)室，南京210098；3.中交天津港灣工程研究院有限公司，天津300000）

舟山海域臺風(fēng)浪數(shù)值模擬

梁連松1，李瑞杰1，2，豐青1，孔從穎3，宋曉波1

（1.河海大學(xué)海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室，南京210098；2.河海大學(xué)環(huán)境海洋實(shí)驗(yàn)室，南京210098；3.中交天津港灣工程研究院有限公司，天津300000）

為提高臺風(fēng)浪模擬的精度，將Holland臺風(fēng)風(fēng)場模型與CCMP背景風(fēng)場相疊加，構(gòu)造合成風(fēng)場來驅(qū)動SWAN模型。運(yùn)用Jason?2衛(wèi)星數(shù)據(jù)進(jìn)行驗(yàn)證，比較了不同的最大風(fēng)速半徑計(jì)算公式和Holland B參數(shù)組合構(gòu)造出的合成風(fēng)場對模擬結(jié)果影響。選取最優(yōu)組合，運(yùn)用自嵌套，模擬了臺風(fēng)“米雷”通過舟山海域的波浪場。結(jié)果表明舟山群島對臺風(fēng)浪的阻擋效果明顯，臺風(fēng)期間舟山東部海域波高大、周期長且涌浪影響明顯，西部海域波高較小且以風(fēng)浪影響為主。

臺風(fēng)浪；SWAN；Holland；Jason?2；舟山

在臺風(fēng)浪的數(shù)值模擬中，波浪場的模擬精度很大程度上依賴于海面風(fēng)場的計(jì)算精度。通過在經(jīng)驗(yàn)臺風(fēng)風(fēng)場中加入再分析風(fēng)場資料，能有效提高海面風(fēng)場的精度。常見的臺風(fēng)風(fēng)場模型有藤田模型、Myers模型、Jelenianski模型、Fujita模型等，其中Holland臺風(fēng)風(fēng)場模型能有效提高臺風(fēng)過境期間大風(fēng)、大浪的模擬精度［1］。

舟山群島港灣眾多，臺風(fēng)過境期間，臺風(fēng)浪往往對其海岸工程、海洋交通造成較大威脅［1-2］。本文以臺風(fēng)“米雷”為例，利用SWAN模型，通過與Jason?2衛(wèi)星波高數(shù)據(jù)的比較，比較不同Holland風(fēng)場模型中最大風(fēng)速半徑和Holland B參數(shù)計(jì)算組合對臺風(fēng)浪模擬精度的影響，并分析了臺風(fēng)“米雷”在舟山海域的臺風(fēng)浪分布特征。

1 風(fēng)場和波浪模型

1.1 波浪模型

SWAN模型在笛卡爾坐標(biāo)系下控制方程為

式中：N為作用量密度；σ為頻率；θ為傳播方向；Stot為源匯項(xiàng)；Cx、Cy、Cσ、Cθ分別為在x、y、σ、θ空間上的群速度。SWAN是Delft大學(xué)提出的考慮地形變化影響下的波浪淺化、破碎、波-波非線性相互作用、風(fēng)能輸入和波浪耗散等效應(yīng)的波浪模型，能用于近岸海域等淺水區(qū)域的波浪計(jì)算［1］，模型采用全隱格式，具有無條件穩(wěn)定的特點(diǎn)。

1.2Holland臺風(fēng)風(fēng)場模型

臺風(fēng)風(fēng)場通常利用梯度風(fēng)原理，由臺風(fēng)氣壓場計(jì)算得到。1980年，Holland在Schloemer指數(shù)型氣壓分布模型的基礎(chǔ)上，引入Holland B參數(shù)，用梯度風(fēng)方程得到Holland臺風(fēng)風(fēng)場模型［3］

表1 最大風(fēng)速半徑Rmax和Holland B參數(shù)計(jì)算公式Tab.1Summary of Rmax and Holland B estimation methods

由式（2）可知，當(dāng)臺風(fēng)中心氣壓確定后，Holland B參數(shù)和最大風(fēng)速半徑即可確定Holland臺風(fēng)風(fēng)場模型。最大風(fēng)速半徑能夠反映出臺風(fēng)氣壓分布的陡度，Holland B參數(shù)則確定臺風(fēng)氣壓剖面和風(fēng)速剖面的形狀。表1列舉了國內(nèi)外學(xué)者提出的最大風(fēng)速半徑和B參數(shù)的計(jì)算公式。表1中，φ是臺風(fēng)中心緯度，V為臺風(fēng)中心移動速度，為飛行層最大風(fēng)速，此處飛行層氣壓為700 hPa，高度約3 km，其值為［4］

臺風(fēng)風(fēng)場往往具有不對稱性。然而Holland理想臺風(fēng)模型是圓對稱風(fēng)場，造成了臺風(fēng)模型應(yīng)用的局限性。為此，在臺風(fēng)外圍采用背景風(fēng)場，臺風(fēng)中心附近仍用臺風(fēng)模型風(fēng)場，構(gòu)造出合成的臺風(fēng)風(fēng)場［10］

式中：Vi為合成風(fēng)場；Vh為Holland風(fēng)場；Vc是背景風(fēng)場；e為權(quán)重系數(shù)系數(shù)n通常取9或10。背景風(fēng)場采用NASA研發(fā)的CCMP風(fēng)場，其空間分辨率為0.25°×0.25°；時(shí)間分辨率為6 h。

2 數(shù)值模擬

2.1 模型設(shè)置

文中用SWAN模型模擬“米雷”經(jīng)過舟山海域時(shí)的臺風(fēng)浪過程。“米雷”于2011年6月22日生成，24日下午加強(qiáng)為強(qiáng)熱帶風(fēng)暴，26日晚在山東登陸，27日上午減弱后于朝鮮再次登陸。臺風(fēng)期間在東海海面形成4～8 m高的臺風(fēng)浪，國家海洋局QF207浮標(biāo)實(shí)測最大有效波高達(dá)6.7 m［11］。

大范圍模型采用非結(jié)構(gòu)三角網(wǎng)格，網(wǎng)格數(shù)72 096，網(wǎng)格點(diǎn)數(shù)36 877（圖1）。初始化方法為基于初始輸入有限風(fēng)場的JONSWAP譜，方向分為36段，分辨率為10°。計(jì)算時(shí)間從6月21日～6月29日，時(shí)間步長取600 s。模型考慮海底摩擦、波浪破碎、三波和四波相互作用。自嵌套的小范圍模擬區(qū)域覆蓋舟山海域，網(wǎng)格數(shù)85 928，網(wǎng)格點(diǎn)43 931。由于舟山海域島嶼眾多，小范圍模型還考慮了繞射項(xiàng)，其余設(shè)置與大范圍模型一致。模型所使用的地形資料來自于ETOPO1全球地形數(shù)據(jù)集。臺風(fēng)路徑以及中心氣壓來源于中國氣象局熱帶氣旋資料中心的CMA最佳路徑數(shù)據(jù)集。

圖11105 臺風(fēng)“米雷”移動路徑、Jason-2衛(wèi)星軌道及SWAN模型大、小范圍計(jì)算區(qū)域Fig.1Computational area and 1105 typhoon moving path and Jason?2 track

2.2 風(fēng)場方案

為比較不同的Holland B參數(shù)和最大風(fēng)速半徑Rmax計(jì)算公式構(gòu)造出的臺風(fēng)風(fēng)場對臺風(fēng)浪模擬精度的影響，文中設(shè)置了6組實(shí)驗(yàn)方案。方案設(shè)置如表2所示，其中T2、T3、T4的B參數(shù)均選用Hubbert的經(jīng)驗(yàn)公式，但三者選取的最大風(fēng)速半徑不同；T4、T5、T6選用Willoughby的最大風(fēng)速半徑公式，B參數(shù)選取不同。

表2 實(shí)驗(yàn)方案設(shè)置Tab.2Numbers and content of the experiments

Jason?2衛(wèi)星于2008年發(fā)射，其海面測高精度可達(dá)3.4 cm。文中選用Jason?2衛(wèi)星高度計(jì)Ku波段的波浪資料作為驗(yàn)證。臺風(fēng)期間，衛(wèi)星正處于cycle109，有3條軌道經(jīng)過模擬區(qū)域，分別為P214、P229和P240（圖1）。

圖2為模擬波高和衛(wèi)星波高的對比圖?？傮w上，散點(diǎn)較為均勻地分布在直線x=y的兩側(cè)，表明模擬結(jié)果較好，且當(dāng)衛(wèi)星波高大于4 m時(shí)，合成風(fēng)場的散點(diǎn)分布更為接近直線x=y。為定量分析各實(shí)驗(yàn)方案模擬效果，采用平均相對誤差EMAE、均方根誤差Rmse和相關(guān)系數(shù)r對臺風(fēng)浪模擬結(jié)果和衛(wèi)星波高資料進(jìn)行比較分析。

式中：yi為模擬值；xi為衛(wèi)星數(shù)據(jù)；xˉ、yˉ分別為衛(wèi)星數(shù)據(jù)和模擬值的平均值；N為樣本總數(shù)，此處共371個(gè)數(shù)據(jù)點(diǎn)。

由表3模擬波高和衛(wèi)星數(shù)據(jù)對比的統(tǒng)計(jì)結(jié)果可知，各方案中相關(guān)系數(shù)r均大于0.9，表明模擬結(jié)果與衛(wèi)星波高數(shù)據(jù)相關(guān)性較高。合成風(fēng)場的平均相對誤差均不超過0.15，除T2中均方根誤差要稍大于T1外，合成風(fēng)場的計(jì)算結(jié)果較T1均有一定程度的改善，可見合成風(fēng)場能夠提高CCMP風(fēng)場模擬有效波高的精度。對比T2、T3、T4，發(fā)現(xiàn)不同最大風(fēng)速半徑計(jì)算公式的選取對臺風(fēng)浪的模擬結(jié)果影響較大，且T4的模擬結(jié)果要優(yōu)于T2、T3，表明Willoughby提出的最大風(fēng)速半徑的計(jì)算公式模擬效果較好。比較T4，T5，T6，三者的平均相對誤差相差不大，這可能是衛(wèi)星數(shù)據(jù)所處的軌道離臺風(fēng)中心有一定距離及風(fēng)場分辨率的限制導(dǎo)致的?？傮w上T6組合的模擬結(jié)果最優(yōu)，T4次之。用T6組合構(gòu)造的合成風(fēng)場作為SWAN模型的輸入風(fēng)場。圖3為T6組合的合成風(fēng)場下輸出大范圍風(fēng)場矢量和有效波高分布圖。

表3 實(shí)驗(yàn)?zāi)M有效波高與衛(wèi)星數(shù)據(jù)的對比Tab.3Comparisons of significant wave height with satellite data

圖2 實(shí)驗(yàn)方案的衛(wèi)星數(shù)據(jù)和模式模擬波高對比的散點(diǎn)圖（圖中實(shí)線表示x=y）Fig.2Comparison between model results and satellite data

圖3 大范圍風(fēng)場矢量和有效波高分布圖Fig.3Distribution of significant wind speed，wave height at large computation area

2.3 結(jié)果分析

圖4給出了從24日21時(shí)～26日11時(shí)臺風(fēng)“米雷”經(jīng)過舟山海域時(shí)波高分布圖。當(dāng)臺風(fēng)逼近時(shí)，舟山東南部海域的波高最先增加，隨著臺風(fēng)中心的北移，最大波高位置也隨之北移，同時(shí)，波浪也沿著群島之間的水道向舟山西部海域傳播。從總體上看，群島大致為南北走向，越靠近群島波高等值線分布越密，東部海面波高明顯要高于群島西部海域，可見舟山群島對臺風(fēng)浪起到了阻擋作用。

圖4 舟山海域波高分布圖Fig.4Distribution of wave height in Zhoushan

為更好地分析臺風(fēng)“米雷”對舟山海域的影響，依據(jù)舟山海域的地形特征選取了3個(gè)特征站位，東面海域?yàn)锳、B站位，水深大約為40 m，分別代表舟山東南部外海和東北部外海，C站位則處于群島西部海域（圖5）。

觀測點(diǎn)在臺風(fēng)期間的波高變化如圖5～圖6所示。初始時(shí)，A、B、C的波高均處在0.6～1.1 m，是當(dāng)?shù)氐娘L(fēng)場并非“米雷”臺風(fēng)系統(tǒng)的風(fēng)場所造成。24日6時(shí)，A站位涌浪波高開始增大，B站位增幅較小，此時(shí)臺風(fēng)中心距離A站位900 km左右。隨后，A、B兩站位海域的混合浪波高迅速增加，當(dāng)臺風(fēng)中心距A站位120 km時(shí)，A站位于25日13時(shí)到達(dá)混合浪波高最大值5.36 m，Tm01達(dá)到10 s，涌浪波高在短短30 h內(nèi)就增加到了3.9 m，表明A站位在臺風(fēng)逼近的過程中受涌浪影響較大。25日13時(shí)后，A站位波高下降明顯，B站位則于17時(shí)波高升到峰值，達(dá)到4.2 m，Tm01達(dá)到了10 s。C站位于19時(shí)波高達(dá)到最大。25日19時(shí)后，臺風(fēng)中心開始遠(yuǎn)離，3處海域波高均下降。

相比較而言，位于舟山西部海域的C站位的混合浪波高則一直處于1.5 m以下，涌浪波高最大時(shí)僅0.18 m，Tm01維持在1.8～4 s。由于舟山群島的阻擋作用，僅部分涌浪得以從島嶼之間的水道傳入，當(dāng)臺風(fēng)中心開始遠(yuǎn)離舟山海域時(shí)，C站位海域波高下降較A、B站位明顯，表明在此過程中舟山西部海域的波浪以風(fēng)浪為主。

圖5 小區(qū)域計(jì)算范圍及觀測點(diǎn)布置Fig.5Computed area by SWAN at small grid and distribution of characteristic stations

圖6 特征站點(diǎn)的波高過程線Fig.6 Variation of mixed wave height and swell with time in A，B，C

3 結(jié)語

文中采用Holland臺風(fēng)模型與CCMP背景風(fēng)場疊加后的合成風(fēng)場模擬臺風(fēng)風(fēng)場，通過SWAN模型計(jì)算臺風(fēng)浪，利用Jason?2衛(wèi)星波浪資料作為驗(yàn)證，探討了不同的最大風(fēng)速半徑計(jì)算公式和B參數(shù)組合構(gòu)造出的風(fēng)場對臺風(fēng)浪模擬精度的影響。表明T6組合即Willoughby的最大風(fēng)速半徑計(jì)算公式和Vickery的Holland B參數(shù)經(jīng)驗(yàn)公式組合計(jì)算效果最優(yōu)。通過自嵌套，模擬了臺風(fēng)經(jīng)過舟山海域附近時(shí)的臺風(fēng)浪。模型中未考慮島嶼反射現(xiàn)象，可能對島嶼附近波高有一定影響。

結(jié)果表明舟山群島對臺風(fēng)浪的阻擋效果明顯。群島大致為南北走向，臺風(fēng)期間，越靠近群島波高等值線分布越密，群島東部海面波高明顯要高于西部海域的波高；臺風(fēng)接近舟山海域500～900 km時(shí)，涌浪開始出現(xiàn)在舟山海域，東南部外海上升最快且最大，由于群島的阻擋作用，舟山西部海域最晚；舟山東部海域在臺風(fēng)期間涌浪作用明顯，波高大，波浪周期較長。西部海域則以風(fēng)浪作用為主，波浪周期小。

［1］唐建.基于臺風(fēng)風(fēng)場模型的臺風(fēng)浪數(shù)值模擬［J］.海洋湖沼通報(bào)，2013（2）:24-30. TANG J.Numerical simulation of typhoon waves with typhoon wind model［J］.Transactions of oceanology and Limnology，2013（2）:24-30.

［2］姬厚德，藍(lán)尹余，趙東波.SWAN在臺灣海峽臺風(fēng)浪數(shù)值模擬中的應(yīng)用研究［J］.水道港口，2013，34（4）:113-117. JI H D，LAN Y Y，ZHAO D B.Application of SWAN in Taiwan Strait typhoon waves numerical simulation［J］.Journal of Water?way and Harbor，2013，34（4）:113-117.

［3］Holland G J.An Analytic Model of the Wind and Pressure Profiles in Hurricanes［J］.Monthly Weather Review，1980，108（8）: 1 212-1 218.

［4］林偉，方偉華.西北太平洋臺風(fēng)風(fēng)場模型中Holland B系數(shù)區(qū)域特征研究［J］.熱帶地理.2013，33:124-132. LIN W，F(xiàn)ANG W H.Rigional Characteristics of Holland B Parameter in Typhoon Wind Field Model for Northwest Pacific［J］. Tropical Geography，2013，33:124-132.

［5］Graham H E，Nunn D E.National Hurricane Research Project，Report No.33:Meteorological conditions pertinent to standard project hurricane Atlantic and Gulf Coasts of the United States［R］.US:US Department of Commerce and Weather Bureau，1959.

［6］Willoughby H E，Rahn M E.Parametric representation of the primary hurricane vortex.Part I:observations and evaluation of the holland（1980）model［J］.Monthly Weather Review，2004，132（12）:3 033-3 048.

［7］江志輝，華鋒.一個(gè)新的熱帶氣旋參數(shù)調(diào)整方案［J］.海洋科學(xué)進(jìn)展，2008，26（1）:1-7. JIANG Z H，HUA F.A New Scheme for Adjusting the Tropical Cyclone Parameters［J］.Adwancs in Marine Science，2008，26（1）: 1-7.

［8］Hubbert G D，Holland G J，Leslie L M，et al.A real?time system for forecasting tropical cyclone storm surges［J］.Weather and Forecasting，1991，6（1）:86-97.

［9］Vickery P J，Skerlj P F，Steckley A C，et al.Hurricane wind field model for use in hurricane simulations［J］.Journal of Structural Engineering，2000，126（10）:1 203-1 221.

［10］聞斌，汪鵬，萬雷，等.中國近海海域臺風(fēng)浪模擬實(shí)驗(yàn)［J］.海洋通報(bào)，2008，27（3）:1-6. WEN B，WANG P，WAN L，et al.Experiments on the simulation of typhoon waves in the inshore area of China sea［J］.MarineScience Bullein，2008，27（3）:1-6.

［11］王宏.中國海洋年鑒［J］.海洋災(zāi)害與海洋環(huán)境預(yù)報(bào)服務(wù)，2012:296-298. WANG H.China ocean yearbook［J］.Marine Disasters and Marine Environmental Forecasting Service，2012:296-298.

長興島恒力石化碼頭工程使用港口岸線獲批復(fù)

本刊從交通運(yùn)輸部獲悉，長興島恒力石化精對苯二甲酸項(xiàng)目配套液體化工碼頭工程、配套散雜貨碼頭工程和2萬t級散貨碼頭工程使用港口岸線近日已獲得交通運(yùn)輸部批復(fù)。3個(gè)碼頭工程位于大連港長興島港區(qū)北岸作業(yè)區(qū)，其中，恒力石化精對苯二甲酸項(xiàng)目配套液體化工碼頭工程將建設(shè)2個(gè)10萬t級液體化工泊位，設(shè)計(jì)年通過能力638萬t。交通運(yùn)輸部同意工程可行性研究報(bào)告推薦的總平面布置方案，按640 m泊位長度使用所對應(yīng)的港口岸線。恒力石化項(xiàng)目配套散雜貨碼頭和2萬t級散貨碼頭工程將建設(shè)3個(gè)5萬t級散雜貨泊位，設(shè)計(jì)年通過能力404萬t，泊位長度816 m；建設(shè)一個(gè)2萬t級散貨泊位(水工結(jié)構(gòu)按照靠泊5萬t級船舶設(shè)計(jì)和建設(shè))，設(shè)計(jì)年通過能力128萬t，泊位長度242 m。交通運(yùn)輸部同意工程可行性研究報(bào)告推薦的總平面布置方案，按1 058 m泊位長度使用所對應(yīng)的港口岸線。（殷缶，梅深）

甌江航道整治工程麗水段工程可行性研究報(bào)告獲省發(fā)改委批復(fù)

本刊從交通運(yùn)輸部獲悉，甌江航道整治工程麗水段（船寮鎮(zhèn)黃言村—溫溪鎮(zhèn)馱灘上游）工程可行性研究報(bào)告已于2014年11月3日獲浙江省發(fā)改委批復(fù)。工程起自船寮鎮(zhèn)黃言村，終于溫溪鎮(zhèn)馱灘上游麗水、溫州交界處，全場約30 km。其中船寮鎮(zhèn)黃言村至青田溫溪大橋段約27.5 km，按內(nèi)河四級天然河道和渠化河流航道標(biāo)準(zhǔn)整治；青田溫溪大橋至溫溪鎮(zhèn)馱灘上游段約2.5 km，按可乘潮通航1 000 t級海輪的內(nèi)河三級天然河流航道標(biāo)準(zhǔn)整治；三溪口和青田樞紐船閘按四級標(biāo)準(zhǔn)建設(shè)。本項(xiàng)目航道疏浚約125萬m3，新建三溪口和青田樞紐船閘2座，新建護(hù)岸約1.2 km，新建橋梁防撞設(shè)施14座，錨地、服務(wù)區(qū)、船閘遠(yuǎn)方調(diào)度站各1座，配套建設(shè)導(dǎo)助航設(shè)施和信息化配套設(shè)施。項(xiàng)目總投資估算為67 851萬元，建設(shè)資金除國家和省級內(nèi)河建設(shè)資金補(bǔ)助外，其余由麗水市政府財(cái)政負(fù)責(zé)落實(shí)。（殷缶，梅深）

Numerical simulation of typhoon wave in Zhoushan

LIANG Lian?song1,LI Rui?jie1,2,FENG Qing1,KONG Cong?ying3,SONG Xiao?bo1
（1.Key Laboratory of Coastal Disaster and Defence,Ministry of Education,Hohai University,Nanjing 210098, China；2.Laboratory of Ocean Environment,Hohai University,Nanjing 210098,China;3.Tianjin Port Engineering Institute Ltd.of CCCC，Tianjin 300000,China）

To improve the accuracy of near?shore typhoon wave simulation,numerical simulations were made for typhoon during the period of typhoon Meari,based on SWAN,with combination of the Holland model wind field and CCMP wind field.The Jason?2 satellite data were used to verify the simulation results.The merged wind fields constructed by different maximum wind speed radius formula and parameter B were compared.Then,the typhoon Meari was simulated by adopted nested grid in SWAN.The results demonstrate that the retarding effect of Zhoushan Archipelago for typhoon wave is significant.During typhoon period,the significant wave height in the offshore areas of eastern Zhoushan Archipelago is higher than the one of the west.The offshore areas of eastern Zhoushan Archipel?ago are greatly affected by surge while the western Zhoushan Archipelago areas are dominated by wind wave.

typhoon wave；SWAN；Holland；Jason?2；Zhoushan Archipelago

TV 143；O 242.1

A

1005-8443（2014）06-0582-07

2014-06-09；

2014-08-04

國家自然科學(xué)基金項(xiàng)目（41276017）；國家海洋局海洋公益性行業(yè)科研專項(xiàng)項(xiàng)目（201205005）

梁連松（1990-），男，浙江省蒼南人，碩士研究生，主要從事港口海岸及近海工程方向研究。

Biography：LIANG Lian?song(1990-)，male，master student.