何 佳

(自然資源部第三海洋研究所，福建 廈門(mén) 361005)

閩東地區(qū)位于中國(guó)東南沿海、福建省東部，包括福州、寧德兩市，北接浙江溫州，東臨臺(tái)灣海峽，其海岸線漫長(zhǎng)曲折、港灣眾多、海島棋布，海域遼闊，是發(fā)展港口航運(yùn)、海洋工程、海產(chǎn)捕撈及養(yǎng)殖、濱海旅游的理想?yún)^(qū)域，亦是西北太平洋風(fēng)暴登陸最頻繁的區(qū)域[1].熱帶風(fēng)暴是一種極具破壞性的災(zāi)害性海洋天氣系統(tǒng)，當(dāng)其移近和登陸時(shí)帶來(lái)的狂風(fēng)、巨浪、風(fēng)暴潮等災(zāi)害會(huì)極大危及國(guó)民經(jīng)濟(jì)建設(shè)和生命財(cái)產(chǎn)安全.在全球變暖、風(fēng)暴生成年總數(shù)量和登陸我國(guó)風(fēng)暴的年頻數(shù)58a來(lái)呈線性增長(zhǎng)趨勢(shì)的背景下[1]，對(duì)影響閩東海域典型臺(tái)風(fēng)過(guò)程進(jìn)行模擬研究，可為海岸帶規(guī)劃及海上防災(zāi)減災(zāi)等提供科學(xué)依據(jù).

近年來(lái)，第三代淺水波浪數(shù)值模型(SWAN)已廣泛應(yīng)用于中國(guó)東南沿海臺(tái)風(fēng)浪的研究.韓樹(shù)宗等(2013)采用SWAN模式對(duì)1981—2010年92次臺(tái)風(fēng)進(jìn)行模擬，據(jù)此分析了其研究海域臺(tái)風(fēng)浪有效波高的分布特征[2].袁凱瑞等(2014)、姬厚德等(2013)、鄭祥靖等(2013)、唐艷平等(2017)采用常用的臺(tái)風(fēng)風(fēng)場(chǎng)模型驅(qū)動(dòng)SWAN模式成功地模擬了臺(tái)灣海峽鄰近海域的臺(tái)風(fēng)浪[3-6].此外，史劍等(2011)、宗芳伊等(2014)將非結(jié)構(gòu)三角網(wǎng)格下的SWAN模式應(yīng)用于黃、渤海及嵊泗海域表面重力波模擬，提升了波浪模擬的精確度[7-8].

學(xué)者們致力于常見(jiàn)臺(tái)風(fēng)風(fēng)場(chǎng)模型的研究，以提升臺(tái)風(fēng)風(fēng)場(chǎng)的計(jì)算精度.張亮等(2015)比較了5種梯度風(fēng)的臺(tái)風(fēng)場(chǎng)模型，發(fā)現(xiàn)Holland風(fēng)場(chǎng)模型通過(guò)調(diào)整和臺(tái)風(fēng)結(jié)構(gòu)相關(guān)的參數(shù)B可適應(yīng)不同的臺(tái)風(fēng)場(chǎng)結(jié)構(gòu)[9].唐建等(2013)采用4種常用的臺(tái)風(fēng)風(fēng)場(chǎng)模型與交叉校準(zhǔn)多平臺(tái)(CCMP)風(fēng)場(chǎng)資料構(gòu)造的合成風(fēng)場(chǎng)對(duì)臺(tái)風(fēng)“米雷”進(jìn)行模擬，結(jié)果表明Holland臺(tái)風(fēng)風(fēng)場(chǎng)模型的模擬效果最好[10].楊萬(wàn)康等(2017)基于Holland風(fēng)場(chǎng)模型，通過(guò)臺(tái)風(fēng)“菲特”的后報(bào)試驗(yàn)，研究得到最大風(fēng)速半徑公式和HollandB參數(shù)的最優(yōu)組合[11].蔣小平等(2007)采用Holland模型得到的風(fēng)場(chǎng)疊加美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NCEP)背景風(fēng)場(chǎng)，較好地描述了臺(tái)風(fēng)內(nèi)外風(fēng)區(qū)的分布[12].

本研究基于加密的非結(jié)構(gòu)三角網(wǎng)格，對(duì)水深急劇變化且岸線復(fù)雜的近海區(qū)地形進(jìn)行刻畫(huà),采用Holland理論模型得到的風(fēng)場(chǎng)疊加NCEP海面風(fēng)場(chǎng)資料構(gòu)造的臺(tái)風(fēng)合成風(fēng)場(chǎng)驅(qū)動(dòng)第三代海浪模式SWAN對(duì)2017年影響閩東沿海的“納沙”和“泰利”臺(tái)風(fēng)浪進(jìn)行模擬，并對(duì)臺(tái)風(fēng)期間巨浪的時(shí)空分布特征進(jìn)行分析，以期為閩東沿海臺(tái)風(fēng)浪災(zāi)害預(yù)警和應(yīng)急管理提供技術(shù)支撐和參考依據(jù).

1 模型簡(jiǎn)介與數(shù)據(jù)資料

1.1 海浪模式

由荷蘭Delft大學(xué)(Delft University of Technology)開(kāi)發(fā)的第三代近岸海浪模式SWAN采用動(dòng)譜密度N(σ,θ)為控制變量，能夠準(zhǔn)確合理的模擬復(fù)雜的潮流、地形及風(fēng)場(chǎng)環(huán)境下的波浪場(chǎng).Rogers等(2002)運(yùn)用新的數(shù)值計(jì)算方法，減少了計(jì)算尺度比較大的海域時(shí)產(chǎn)生的數(shù)值耗散，使之適用于淺水到大洋各種尺度的波浪數(shù)值模擬[13].

在球坐標(biāo)系下，動(dòng)譜平衡方程表示如下:

(1)

式(1)中：N為動(dòng)譜密度(m2/Hz)，t代表時(shí)間(s)，λ為經(jīng)度(°)，φ為緯度(°)，σ為相對(duì)波頻(Hz)，θ為波向(°)，cλ、cφ、cσ和cθ分別表示動(dòng)譜密度在λ、φ、σ和θ空間上的傳播速度(m2/s)，Stot代表以譜密度表示的源匯項(xiàng)，包含風(fēng)能輸入、白浪、破碎、海底摩擦、波—波非線性相互作用等物理過(guò)程.

1.2 臺(tái)風(fēng)風(fēng)場(chǎng)

臺(tái)風(fēng)過(guò)程中風(fēng)速的觀測(cè)不容易進(jìn)行，風(fēng)速遙感資料在臺(tái)風(fēng)中心附近明顯偏小.目前國(guó)內(nèi)外普遍根據(jù)梯度風(fēng)原理，采用臺(tái)風(fēng)氣壓場(chǎng)計(jì)算臺(tái)風(fēng)風(fēng)場(chǎng)，包括理論氣壓模型和經(jīng)驗(yàn)?zāi)Ｐ?將臺(tái)風(fēng)氣壓場(chǎng)計(jì)算的風(fēng)場(chǎng)與遙感實(shí)測(cè)的風(fēng)場(chǎng)合成，除臺(tái)風(fēng)過(guò)境期間的風(fēng)場(chǎng)外，還可以較好的考慮臺(tái)風(fēng)生成和衰減后遠(yuǎn)離臺(tái)風(fēng)中心區(qū)域的風(fēng)場(chǎng)[11-12].

1.2.1Holland理論模型風(fēng)場(chǎng) 在Schloemer(1954)指數(shù)氣壓分布模型基礎(chǔ)上，澳大利亞學(xué)者Holland[14](1980)引入B參數(shù)，利用梯度風(fēng)風(fēng)場(chǎng)得到如下的臺(tái)風(fēng)風(fēng)場(chǎng)和氣壓方程：

Vg(r)=

(2)

(3)

式(2、3)中：Vg(r)為距離臺(tái)風(fēng)中心r處的風(fēng)速(m/s)，Pc為臺(tái)風(fēng)中心氣壓(hPa)，Pn為外圍氣壓(取值為1013hPa)，Rmax為最大風(fēng)速半徑(km)，ρa(bǔ)為空氣密度(取值為1.29kg/m3)，f為科氏力參數(shù)，r為計(jì)算點(diǎn)與臺(tái)風(fēng)中心的距離(km)，B參數(shù)決定了臺(tái)風(fēng)的峰度和強(qiáng)度，P(r)為距離臺(tái)風(fēng)中心r處的氣壓(hPa).

國(guó)內(nèi)外學(xué)者對(duì)于Rmax和B參數(shù)進(jìn)行了大量的研究.根據(jù)文獻(xiàn)[11]的研究，本研究Rmax采用Willoughby(2004)風(fēng)速半徑公式[15]，B參數(shù)采用Vickery(2000)的HollandB參數(shù)經(jīng)驗(yàn)公式[16].計(jì)算式如下:

Rmax=51.6exp(-0.0223Vmax+0.0281φ)

(4)

B=1.881-0.00557Rmax-0.01295φ

(5)

式(4、5)中：Vmax為最大風(fēng)速(m/s)，φ表示緯度(°).

1.2.2合成風(fēng)場(chǎng) 本研究采用Holland理論模型得到的風(fēng)場(chǎng)疊加NCEP海面風(fēng)場(chǎng)資料構(gòu)造的臺(tái)風(fēng)合成風(fēng)場(chǎng).NCEP海面風(fēng)場(chǎng)資料源于美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心提供的全球大氣再分析產(chǎn)品CFSv2(https://rda.ucar.edu/datasets/ds094.1/#description)，該風(fēng)場(chǎng)為1h時(shí)間間隔，經(jīng)緯度坐標(biāo)下0.2°水平分辨率，距海面10m高處的風(fēng)速資料.

Holland理論模型風(fēng)場(chǎng)與NCEP海面風(fēng)場(chǎng)資料以一定權(quán)重系數(shù)相疊加構(gòu)造成相應(yīng)的SWAN模型風(fēng)場(chǎng).其計(jì)算公式如下:

Vc=VH×(1-e)+e·VNCEP

(6)

式(6)中：Vc為合成風(fēng)場(chǎng)，VH為Holland模型風(fēng)場(chǎng)，VNCEP為NCEP海面風(fēng)場(chǎng)，e為權(quán)重系數(shù)，e=C4/(1+C4)，C為臺(tái)風(fēng)影響范圍系數(shù)，取C=r/(9Rmax)，r是計(jì)算點(diǎn)與臺(tái)風(fēng)中心的距離(km)，Rmax為最大風(fēng)速半徑(km)，若r>5Rmax，則取NCEP海面風(fēng)場(chǎng).

1.3 數(shù)據(jù)資料

1.3.1岸線及水深數(shù)據(jù) 岸線采用美國(guó)NOAA網(wǎng)站提供的GSHHG(A Global Self-consistent, Hierarchical, High-resolution Geography Database)數(shù)據(jù)(https://www.ngdc.noaa.gov/mgg/shorelines/gshhs.html)和海圖提取的福建近海岸線數(shù)據(jù).

水深采用美國(guó)NOAA網(wǎng)站提供的全球陸地海洋高程數(shù)據(jù)ETOPO1(https://www.ngdc.noaa.gov/mgg/global/global.html)和海圖提取的閩浙近海水深數(shù)據(jù).

1.3.2臺(tái)風(fēng)數(shù)據(jù) 熱帶氣旋包括除熱帶擾動(dòng)之外的所有等級(jí)氣旋，即熱帶低氣壓、熱帶風(fēng)暴、強(qiáng)熱帶風(fēng)暴、臺(tái)風(fēng)、強(qiáng)臺(tái)風(fēng)和超強(qiáng)臺(tái)風(fēng).2017年太平洋上共生成38個(gè)熱帶氣旋，其中的1709號(hào)臺(tái)風(fēng)“納沙”和1718號(hào)超強(qiáng)臺(tái)風(fēng)“泰利”對(duì)閩東海域造成的影響最為嚴(yán)重.本研究主要對(duì)這兩個(gè)臺(tái)風(fēng)期間的臺(tái)風(fēng)浪進(jìn)行模擬研究.

1709號(hào)臺(tái)風(fēng)“納沙”于2017年7月21日在帕勞附近海面上生成，7月25日升格為熱帶低壓，7月26日升格為熱帶風(fēng)暴，7月27日升格為強(qiáng)熱帶風(fēng)暴，7月28日升格為臺(tái)風(fēng)，其最強(qiáng)達(dá)到臺(tái)風(fēng)級(jí)(40m/s，13級(jí))，7月29日19：40在臺(tái)灣省宜蘭縣東部沿海登陸，7月30日06：00在福建省福清市沿海再次登陸(33m/s，12級(jí)).

1718號(hào)超強(qiáng)臺(tái)風(fēng)“泰利”于2017年9月6日在科斯雷島東北方海面上生成，9月8日升格為熱帶低壓，9月9日升格為熱帶風(fēng)暴，9月11日升格為強(qiáng)熱帶風(fēng)暴，9月12日升格為臺(tái)風(fēng)，9月13日升格為強(qiáng)臺(tái)風(fēng).爾后，“泰利”進(jìn)入東海中部的低海溫海域，對(duì)流趨于衰弱，強(qiáng)度快速下降，并轉(zhuǎn)向NE向移動(dòng)，9月16日降格為臺(tái)風(fēng)，9月17日上午以臺(tái)風(fēng)級(jí)在日本九州島第一次登陸，下午以強(qiáng)熱帶風(fēng)暴級(jí)第二次登陸四國(guó)島，晚上以強(qiáng)熱帶風(fēng)暴級(jí)第三次登陸本州島，9月18日上午以強(qiáng)熱帶風(fēng)暴級(jí)第四次登陸北海道.

臺(tái)風(fēng)路徑和各基本資料來(lái)自浙江省水利廳公布的數(shù)據(jù)(http://typhoon.zjwater.gov.cn/default.aspx)和日本氣象廳整編的路徑數(shù)據(jù)(http://www.jma.go.jp/jma/jma-eng/jma-center/rsmc-hp-pub-eg/trackarchives.html).臺(tái)風(fēng)“納沙”和“泰利”在研究區(qū)域內(nèi)的路徑如圖1所示.

圖1 區(qū)域水深、臺(tái)風(fēng)路徑及浮標(biāo)站位置Fig.1 Bathymetry map with typhoon tracks and positions of mooring buoys

1.4 驗(yàn)證資料

本研究使用的驗(yàn)證資料為閩東沿海錨系浮標(biāo)測(cè)站TW03(25.5005°N，120.3080°E)、TW04(26.2830°N，120.7166°E)和TW05(27.0094°N，121.0008°E)的實(shí)測(cè)風(fēng)及波浪數(shù)據(jù).3個(gè)浮標(biāo)站均位于與岸線平行的50m等深線附近(圖1).

1.5 模型試驗(yàn)

根據(jù)臺(tái)風(fēng)“納沙”和“泰利”的移行路徑和時(shí)間，風(fēng)場(chǎng)的計(jì)算區(qū)域?yàn)?6°～31°N，116°～131°E.臺(tái)風(fēng)“納沙”和“泰利”的模擬時(shí)間分別為2017年7月27—30日和2017年9月12—17日.

采用第三代海浪譜模型SWAN的最新版本41.20對(duì)海浪進(jìn)行模擬.模式計(jì)算區(qū)域?yàn)轱L(fēng)場(chǎng)覆蓋下的海域，采用球坐標(biāo)系下的非靜態(tài)模式，頻率的計(jì)算范圍為0.05～1.00Hz，劃分為24個(gè)分布；方向分為36段，分辨率為10°，計(jì)算時(shí)間步長(zhǎng)取600s.本次模擬所采用水深的基面為平均海平面，且不考慮潮汐影響.

使用SMS軟件對(duì)岸線數(shù)據(jù)進(jìn)行處理，構(gòu)建出研究海域精細(xì)化的非結(jié)構(gòu)三角網(wǎng)格(圖2，圖中藍(lán)框范圍內(nèi)為重點(diǎn)區(qū)域，該區(qū)域內(nèi)的網(wǎng)格進(jìn)行了加密處理).然后，再將水深數(shù)據(jù)插值到三角網(wǎng)格點(diǎn)上，生成SWAN模式所需的地形(圖1).

圖2 研究使用非結(jié)構(gòu)三角網(wǎng)格Fig.2 Unstructured triangular grids used in this study

2 結(jié)果與討論

2.1 風(fēng)場(chǎng)模擬結(jié)果

由臺(tái)風(fēng)“納沙”模型計(jì)算風(fēng)速與浮標(biāo)觀測(cè)值的時(shí)間序列比較(圖3)可見(jiàn): 臺(tái)風(fēng)在臺(tái)灣省宜蘭縣第一次登陸前，其中心逐漸向臺(tái)灣島靠近且強(qiáng)度不斷加強(qiáng)，期間3個(gè)浮標(biāo)處的實(shí)測(cè)及模擬風(fēng)速逐漸增大至15m/s；隨著臺(tái)風(fēng)登陸宜蘭并西移進(jìn)入臺(tái)灣海峽，3個(gè)浮標(biāo)站的風(fēng)速迅速攀升至峰值(TW03站約為30m/s，TW04和TW05站約為20m/s)；當(dāng)臺(tái)風(fēng)第二次在福建沿海登陸以后，實(shí)測(cè)及模擬風(fēng)速迅速減小.

由臺(tái)風(fēng)“泰利”模型計(jì)算風(fēng)速與浮標(biāo)觀測(cè)值的時(shí)間序列比較(圖3)可見(jiàn):9月12日期間臺(tái)風(fēng)中心向西北迅速移動(dòng)，逐漸靠近臺(tái)灣島，3個(gè)浮標(biāo)站的實(shí)測(cè)及模擬風(fēng)速亦迅速增大；13日期間3個(gè)浮標(biāo)站的風(fēng)速接近15m/s；15日后隨著臺(tái)風(fēng)中心向東北移動(dòng)，逐漸遠(yuǎn)離臺(tái)灣島，3個(gè)站的風(fēng)速逐漸減小.

對(duì)臺(tái)風(fēng)“納沙”、“泰利”期間模型計(jì)算風(fēng)速與實(shí)測(cè)值進(jìn)行比較(圖4)，可見(jiàn)3個(gè)站模型計(jì)算的風(fēng)速值與實(shí)測(cè)值符合較好，相關(guān)系數(shù)均在0.9以上，表明合成風(fēng)場(chǎng)能較好地模擬臺(tái)風(fēng)期間的風(fēng)速變化過(guò)程，可為臺(tái)風(fēng)浪的模擬提供較高精度的風(fēng)場(chǎng)數(shù)據(jù).

圖3 2017年臺(tái)風(fēng)“納沙”和 “泰利”期間各觀測(cè)站模型計(jì)算風(fēng)速與浮標(biāo)觀測(cè)值的時(shí)間序列比較Fig.3 Time series comparison of simulated wind speed with the observed data from mooring buoys during typhoon Nesat and Talim in 2017

圖4 2017年臺(tái)風(fēng)“納沙”和 “泰利”期間各觀測(cè)站模型計(jì)算風(fēng)速與實(shí)測(cè)值比較Fig.4 Comparison between simulated and measured wind speed during typhoon Nesat and Talim in 2017

2.2 臺(tái)風(fēng)浪模擬結(jié)果

將SWAN模式模擬的“納沙”和“泰利”臺(tái)風(fēng)浪場(chǎng)數(shù)據(jù)插值到浮標(biāo)所在位置，并對(duì)模式計(jì)算有效波高與實(shí)測(cè)值進(jìn)行比較.結(jié)果表明：模型計(jì)算的有效波高與實(shí)際觀測(cè)值的時(shí)序變化一致性較好(圖5)，且相關(guān)系數(shù)均高于0.9(圖6)，模型能較好地模擬波浪對(duì)臺(tái)風(fēng)的響應(yīng)過(guò)程，其模擬結(jié)果可為分析臺(tái)風(fēng)浪的時(shí)空分布特征提供較高精度的波浪場(chǎng)數(shù)據(jù).

圖5 2017年臺(tái)風(fēng)“納沙”和 “泰利”期間各觀測(cè)站模型計(jì)算有效波高與浮標(biāo)觀測(cè)值的時(shí)間序列比較Fig.5 Time series of simulated significant wave height comparison against the observed data from mooring buoys during typhoon Nesat and Talim in 2017

圖6 2017年臺(tái)風(fēng)“納沙”和 “泰利”期間各觀測(cè)站模式計(jì)算有效波高與實(shí)測(cè)值比較Fig.6 Comparison between simulated and measured significant wave height during typhoon Nesat and Talim in 2017

1709號(hào)臺(tái)風(fēng)“納沙”為西北行路徑的臺(tái)風(fēng)，其中心在15°～25°N之間向西北移動(dòng)，橫穿臺(tái)灣島和臺(tái)灣海峽，在福建中東部第二次登陸.臺(tái)風(fēng)第二次登陸前，閩東沿海海域主要受東北風(fēng)控制，臺(tái)風(fēng)登陸后，風(fēng)向由NE向往NW—S向偏轉(zhuǎn)，由臺(tái)風(fēng)“納沙”典型時(shí)刻有效波高分布(圖7)可見(jiàn)，7月29日20：00(臺(tái)風(fēng)剛在臺(tái)灣省東北角的宜蘭縣登陸)宜蘭以東洋面分布一大范圍的狂濤區(qū)(浪高9～12m)；臺(tái)灣海峽北部海面在東北風(fēng)的控制下，出現(xiàn)NE—SW向舌狀分布的巨浪區(qū)(浪高4～6m)；30日02：00臺(tái)風(fēng)中心越過(guò)臺(tái)灣島，西行進(jìn)入臺(tái)灣海峽北部海面，由于受到海峽地形的約束，波浪場(chǎng)呈NE—SW向橢圓狀分布，北部海域浪高大于南部，閩東沿海遍布大范圍的巨浪到狂浪(浪高4～8m).

圖7 臺(tái)風(fēng)“納沙”典型時(shí)刻有效波高分布Fig.7 Significant wave height distribution of typhoon Nesat at representative time

1718號(hào)超強(qiáng)臺(tái)風(fēng)“泰利”路徑在15°～27°N之間先往NW向移動(dòng)，后在閩浙外海海域轉(zhuǎn)向NE向移動(dòng)，途經(jīng)琉球群島北部海面登陸日本.受其影響，閩東沿海海域臺(tái)風(fēng)過(guò)程中始終受NE—N向風(fēng)控制，由臺(tái)風(fēng)“泰利”典型時(shí)刻有效波高分布(圖8)可見(jiàn)，9月14日13：00—19：00臺(tái)風(fēng)強(qiáng)度升級(jí)至超強(qiáng)臺(tái)風(fēng)，臺(tái)風(fēng)中心偏右側(cè)附近海域極值有效波高為14m以上，121°E以東海面為大范圍的狂浪到狂濤區(qū)，臺(tái)灣海峽北部的閩東沿海(50m等深線以西)出現(xiàn)NE—SW向舌狀分布的巨浪區(qū)(浪高4～6m)；14日20：00臺(tái)風(fēng)往NE向移動(dòng)，強(qiáng)度略微減弱，波浪場(chǎng)的分布與13：00的分布情況相近，由于能量的持續(xù)輸入，大浪區(qū)域范圍更廣，臺(tái)風(fēng)中心鄰近左側(cè)海域的極值波高達(dá)到過(guò)程最大值.

圖8 臺(tái)風(fēng)“泰利”典型時(shí)刻有效波高分布Fig.8 Significant wave height distribution of typhoon Talim at representative time

3 結(jié)論

(1)合成風(fēng)場(chǎng)能較好地模擬臺(tái)風(fēng)期間的風(fēng)速變化過(guò)程，計(jì)算的風(fēng)速值與實(shí)測(cè)值符合較好，模擬的風(fēng)場(chǎng)數(shù)據(jù)精度較高.數(shù)值模擬的有效波高與浮標(biāo)實(shí)測(cè)值的時(shí)序變化一致性較好，表明SWAN模式能較好地模擬近岸波浪對(duì)臺(tái)風(fēng)的響應(yīng)過(guò)程.

(2) 臺(tái)風(fēng)“納沙”期間，當(dāng)臺(tái)風(fēng)中心穿過(guò)臺(tái)灣島進(jìn)入臺(tái)灣海峽北部海面時(shí)，在海峽地形的約束下，波浪場(chǎng)呈橢圓狀分布，閩東沿海遍布大范圍的巨浪到狂浪(浪高4～8m).“泰利”臺(tái)風(fēng)于9月14日13：00—19：00期間強(qiáng)度升至超強(qiáng)臺(tái)風(fēng)級(jí)，臺(tái)風(fēng)中心偏右側(cè)海域極值有效波高達(dá)到14m以上，121°E以東海面為大范圍的狂浪到狂濤，閩東沿海海域受東北風(fēng)影響，巨浪區(qū)呈NE—SW舌狀分布(浪高4～6m).

(3) 較之臺(tái)風(fēng)“泰利”(臺(tái)風(fēng)中心遠(yuǎn)離臺(tái)灣島)，西北行路徑的臺(tái)風(fēng)“納沙”在其臺(tái)風(fēng)中心進(jìn)入臺(tái)灣海峽期間，在閩東沿海造成更為惡劣的海況.

本次模擬暫未考慮潮位變化、波流相互作用對(duì)波浪的影響，這會(huì)致使計(jì)算結(jié)果與實(shí)際存在一定誤差，特別是在近岸淺水區(qū)(水深20m以內(nèi)).在后續(xù)近岸波浪的模擬工作中將會(huì)對(duì)此做進(jìn)一步的研究.