臺風影響下西沙海域與島礁內外的波浪分布特征梁海萍;馮偉忠;梁海燕;石軍【摘要】通過利用1991-2018年西沙海洋站波浪觀測資料與SWAN模式反演計算結果,初步探究不同路徑熱帶氣旋影響下西沙海域、西沙礁盤與島嶼波浪掩護區海域的波浪分布特征?結果表明:西北型路徑是影響西沙群島的TC最大可能路徑，通常以臺風強度居多,波高以大到巨浪級別為主;西行型路徑是影響西沙群島的TC第二可能路徑，通常以強臺風居多，巨浪為主，海況極差；TC影響下各島嶼波浪掩護區波高基本大于2.0m,甚至接近5.0m,礁盤內波高基本在1.0m以下?本研究結果為西沙海域漁船避風與海洋防災減災工作提供基礎性參考.期刊名稱】《海洋預報》年(卷),期】2019(036)004【總頁數】7頁(P28-34)【關鍵詞】西沙；熱帶氣旋；波浪；不同路徑；分布特征【作者】梁海萍；馮偉忠；梁海燕；石軍【作者單位】國家海洋局海口海洋環境監測中心站,海南海口570100；中國海洋大學,山東青島2661002；國家海洋局南海預報中心,廣東廣州510000；國家海洋局海口海洋環境監測中心站,海南海口570100；國家海洋局海口海洋環境監測中心站,海南海口570100【正文語種】中文中圖分類】P731.22引言西沙群島是連接太平洋與印度洋的重要通道，戰略位置極為重要。西沙海域有著名的西沙漁場，同時又是南下南海船只的避風、補給基地。西沙群島易受熱帶氣旋(TropicalCyclone,TC)侵襲，每年受影響個數約3—4個，由于西沙群島遠離大陸，避風、補給基礎設施遠未滿足要求，利用天然島嶼淺灘避風是目前的主要途徑。西沙海域特殊地理位置導致人們長期對西沙海域波浪分布特征缺乏全面認識。即便2008年以來，海南省安全管理部門相繼投入大量財力物力，利用全球定位系統、AIS(船舶自動識別系統，AutomaticIdentificationSystem)基站、短波通信基站、手機等實現船舶位置動態監控，建成了全省的漁船安全應急指揮系統，但在獲知TC惡劣海洋環境預警信息條件下如何在廣闊的西沙海域就近選擇地點避險，仍然是當前面臨的難題。海浪特征分析的實際應用比較廣泛，劉桂海等［1］利用長期波浪觀測資料，尋求合適的海岸工程波浪要素長期統計分布概率模式；楊憲章［2］通過對波浪和系泊船舶運動的研究和分析，對毛里塔尼亞友誼港修建港口碼頭提出看法與意見；耿寶磊等［3］收集了2012年13號TC“啟德”和23號TC“山神”發生期間的波浪數據，并結合波浪方向譜和頻譜的變化，分析出了瓊州海峽海域在TC發生期間的波浪特征，同時也對風浪和涌浪的特征進行了判別；丁千龍等［4］根據西沙海洋站實測波浪及風要素資料，對西沙海區波浪分布特征進行了分析，通過多元回歸計算，得出波要素的預報方程，后報檢驗結果良好；楊永增等［5］運用MASNUM(KeyLaboratoryofMarineScienceandNumericalModeling)海浪數值預報系統里的全球10a后報數據庫資料，對北印度洋區域波浪分布特征進行了統計分析，發現了夏季的波浪高于冬季，在空間分布上，西部風、浪比東部大；楊永增等［6］建立了球坐標系下的全球海浪數值模式，重點導出了球坐標系下的海浪能量譜平衡方程及其復雜特征線方程；韓淑宗等［7］運用SWAN模式(SimulatingWavesNearshore，近岸海浪數值計算模式)計算出有效波高，采用實測資料驗證，通過模式計算的結果分析出東中國海海域的臺風浪有效波高分布特征；劉永玲等［8-11］參考國外學者對風暴潮模式與波浪之間的相關關系及相關作用，建立了考慮海浪因素的三維風暴潮模式［12］。船只避風方面的研究比較缺乏，丁慧彥等［13］通過比較不同種類船只的不同避風方式，設計出可以給出船只回港優化方案的模型，可有效指導船只安全回港；王傳聰等［14］從避風錨地的底質、防波堤設計標準與錨泊方式等3個方面分別討論了避風型漁港的避風因素，為建設避風型漁港與災害決策提供了科學依據，目前關于西沙海域波浪分布特征的研究較少，大多數學者只停留在研究大范圍波浪分布，未細化到礁盤與島嶼波浪掩護區，而關于西沙海域漁船避風方面尚未有相關研究，本研究探討不同路徑熱帶氣旋影響下西沙海域、礁盤與島嶼波浪掩護區的波浪分布特征，旨在為西沙海域漁船避風與海洋防災減災工作提供基礎性參考。資料來源本研究使用的資料來自于國家海洋局西沙海洋站(16°8'N,112°3'E)的實測波高、風等，統計選取資料時間為1991—2018年。定義在TC活動期間監測到7級及7級以上風速過程稱為受TC影響過程。年鑒資料上將西北太平洋TC共分成13種路徑，結合實際情況，本研究在此基礎上，經過適當歸并，將對影響西沙群島的TC路徑分為以下6種，各類路徑定義［15］如下：西行型：TC移動的全過程中，平均移向大致在250。~289。之間。西北型：TC移動的全過程中，平均移向大致在290。~339。之間。右折型：TC移動路徑中至少有一個折向點，并且前后12h平均移向相比較，右折角度＞45°O左折型：TC移動路徑中至少有一個折向點，并且前后12h平均移向相比較，左折角度＞30°。偏北型：TC移動的全過程中，平均移向是在340。~20。之間。回旋型：TC移動路徑變幻多端，擺動、打轉或停滯少動，無法判定。1991—2018年期間影響西沙群島的TC過程共計89次［16］，其中西北型路徑占的比率最大，占比為43.8%，其次為西行型路徑，占比為22.5%（見表1）。由于其他路徑TC樣本空間少，此處不做研究。結果與討論3.1不同路徑TC波浪分布特征西北型路徑1991—2018年期間影響西沙群島的西北型TC數量為39個，占比約為43.8%（見表1），其中熱帶風暴占10.3%，強熱帶風暴占25.6%，臺風占30.8%，強臺風占25.6%，超強臺風所占比率最小，為7.7%，其統計如表2所示。基于西沙海洋站實測波高條件，圖1表明1991—2018年期間影響西沙群島的西北型TC波高分布主要集中在2.5~3.9m之間，占比為41%;其次是4.0~5.9m之間，占比為28%;2.5m以下波高占比為18%;波高在6.0~8.9m之間占比為13%。從表3可以看出：同一次TC過程中最大波高發生在最大風速之前的占比為12.8%,發生在最大風速之后的占比為61.5%，發生在最大風速發生時的占比為25.7%，由此可見，在西北型TC過程中，最大波高一般出現在最大風速之后，所以要注意提防大風過后的大浪災害。西行型路經1991—2018年期間影響西沙群島的西行型TC數量為20個，占總數的22.5%,其中強臺風所占比最高，為30%;其次為熱帶風暴與強熱帶風暴，占比均為20%，超強臺風與臺風占比均為15%（見表4）。基于西沙海洋站實測波高條件，從圖1可以看出，1991—2018年期間影響西沙群島的TC波高分布主要集中在4.0~5.9m之間，占比為40%;其次是2.5~3.9m之間，占比為35%;2.5m以下波高占比為20%;波高分布在6.0~8.9m之間占比較小，僅為5%。從表5可以看出，同一次TC過程中最大波高發生在最大風速之前的占比為25%，發生在最大風速之后的占比為60%，發生在最大風速發生時的占比為15%。由此可見，在西行型TC過程中，最大波高一般出現在最大風速之后，所以要注意提防大風過后的大浪災害。表11991—2018西沙群島的TC路徑分布路徑西行型西北型右折型左折型偏北型回旋型編號占比/%93239508952196129718001600110123040505180615061608160916100211191321 1325 1519 1723 22.47910391069111920492059207930394199515 9520 9606 9713 9803990301030307032014151621170443.820602 0606 0623 0814 0817100310051108111712231306130913101311 1330 0516 0521 14091719931895169810060106215.62931294249621991503080312071511181329172411.25950696109710080115085.62940394229511951496189902021405080917180411.25圖1各類路徑TCH1/10波高分布表21991—2018年西北型TC各等級強度統計TC強度等級熱帶風暴強熱帶風暴臺風強臺風超強臺風37.69TC數量/個占比/%410.261025.641230.801025.64表31991—2018年西北型路徑最大波高與最大風速關系統計出現時間最大風速之前最大風速之后最大風、浪同時出現TC數量/個占比/%512.822461.541025.65表41991—2018年西行型TC各等級強度統計TC強度等級TC數量/個占比/%熱帶風暴420強熱帶風暴420臺風315強臺風630超強臺風315表51991—2018年西行型路徑最大波高與最大風速關系統計出現時間TC數量/個占比/%最大風速之前525最大風速之后1260最大風、浪同時出現315島礁內外波浪特征西沙海洋站波浪觀測資料多為一天4次，且未能較好代表西沙礁盤與島嶼波浪掩護區狀況，為更全面、具體了解西沙波浪分布特征，本研究采用SWAN波浪模式做進一步研究。3.2.1SWAN模式介紹主要對特定波浪條件下西沙海域中環礁與島嶼附近海域波浪進行研究，主要考慮地形對波浪的影響。采用基于能量平衡方程，適用于近岸、河口、海灣、湖泊的SWAN模式對研究區域進行了模擬。SWAN41.01版模式考慮了波浪的淺變、折射、底摩擦、破碎、白帽、風能輸入、三波相互作用和四波相互非線性效應，因而具有以下4個方面的功能［17］：（1）適用于海岸、河口水域風浪、涌浪及混合浪的預報，在直角坐標和球坐標下以矩形網格進行嵌套計算、以曲線網格進行數值計算，適用于大、中、小水域；（2）模擬水底地形和流場的變化引起的波浪折射、淺化，逆流時波浪的反射和破碎，波浪遇到障礙物的透射及阻礙，波浪增水；（3）能預報計算域內波高、波周期、波長、波陡、波浪行進方向、近底水質點的運動速度、波能傳播方向、能量耗散、單位水面所受波力等海岸工程所需的重要參數。SWAN模式的基本控制方程為動譜能量平衡方程：3.2.1模式推算范圍本研究利用SWAN模式兩重嵌套方法研究西沙海域TC過程波浪，南海區域為：99°~129.5°E,2.0°~30.75°N,空間分辨率為1/20°:西沙海域范圍為：111.35°~112.85°E,15.75°~17.25°N(見圖2),空間分辨率為1/400。。南海區域SWAN模式推算結果為西沙海域模式提供嵌套的波譜邊界條件，風場為WRF模式提供的逐時風場。波浪所需要的地形數據分別取自ETOPO1與海圖融合，西沙海域由航保部海圖(17100,1:25000)讀取，限于西沙詳細地形水深難以獲得，環礁內與邊緣水深概化處理，環礁中部水深取4.0~10.0m，環礁邊緣水深取0.5m。在西沙群島已經形成環礁的礁盤中、或島嶼波浪掩護區選取9個代表點，研究不同波向,不同波級條件下各島礁的波浪條件。SWAN模式采用SWAN41.01版，風能輸入和白帽耗散項采用Komen參數化方案、底摩擦采用Hasselmann的JONSWAP實驗常數理論，考慮Hasselmann四波相互和三相波相互作用，波浪破碎能量耗散采用Battjes和Janssen的borebased模式。利用1321號臺風“蝴蝶”和0518號臺風“達維”在西沙海洋站的波浪觀測資料與模式結果進行比對,模式計算結果絕對平均誤差小于30%,因此可作為研究西沙海域的基礎資料。圖2小區推算范圍和島礁代表點3.2.3模式計算結果TC影響下西沙海域波浪方向隨TC風向改變，本研究按TC期西沙海域各代表點累年各向波浪分布頻率將主要波浪方向確定為NE、S、E、N、W5個方向，每個方向設定兩種情景，其中Hs是按照3.1節中西北型及西行型路徑均在2.5~3.9m與4.0~5.9m波高范圍出現頻率最高結論確定，每個波高范圍取平均值，分別確定為3.0m、5.0m，以研究10種情景下礁盤與島嶼的波浪特征，其中H1/10=1.18Hs轉換。各情景小區邊界條件設置見表6,計算結果見表7。表7表明，在NE方向波浪作用下，無論情景1還是情景2，礁盤內H1/10都小于1.0m，其中北礁、羚羊礁、盤石嶼海況好，島嶼波浪掩護區內H1/10都大于1.0m，情景2時深航島口門H1/10有2.45m，永興島西南區域H1/10高達2.76m；S波浪方向影響下，礁盤內與島嶼波浪掩護區的H1/10都比NE向波浪高，其中七連嶼西南區域與永興島西南區域波浪較大，情景4時永興島西南區域波浪高達3.32m；在E或N方向波浪作用下，除了情景6時華光礁H1/10略大于1.0m外，其他的礁盤內H1/10都小于1.0m，海況較好，E方向波浪作用時各島嶼掩護區H1/10相差小于0.5m，而N方向波浪作用時七連嶼西南區域H1/10明顯小于其他兩個島嶼波浪掩護區；在W方向的情景10中，除了羚羊礁、玉琢礁與琛航島口門外，礁盤內與島嶼掩護區H1/10相對其他海況都明顯偏高，其中永興島西南向高達4.7m，位居第一。表6各情景小區邊界條件情景波浪方向Hs/mTz/s風速/(m/s)NENE12345678910SSEENNWW35353535356.98.96.98.96.98.96.98.96.98.912161216121612161216西沙波浪分布特征成因分析TC大多數發生在南、北緯度的5。~20。，其中10。~20。居多［18］。影響到西沙群島的TC基本都是在菲律賓以東洋面生成后進入南海的，少部分TC是在南海內形成。西沙群島位于南海的西北部海域，所以主要受西北行路徑與西行型路徑影響。圖3顯示臺風主要分布在7—10月份，強臺風主要分布在9—11月份；圖4顯示西北型路徑主要分布在7—8月份，西行型路徑主要分布在8、9、11月份，所以影響西沙的西北型路徑的TC多數以臺風等級出現，而西行型路徑的TC多數以強臺風等級出現，另外在菲律賓以東洋面生成的TC強度一般較大，以臺風與強臺風為主。風浪的成長取決于風速、風時、風區，在傳播過程中風浪不斷從風中獲得能量和上游風浪傳播能量而發展［18］，所以最大波高一般出現在最大風速之后。根據風浪原理，大風浪的產生需同時具備風速、風時、風區3個良好條件，而西沙礁盤基本四面都有礁石遮擋且區域較小，風區條件較差，所以無法形成大浪，礁盤內海況較好；而島嶼波浪掩護區內，一般具備風速、風時、風區3個良好條件，所以H1/10基本都在2m以上，甚至高達4.7m。表7各礁盤和島嶼代表點H1/10計算結果（單位：m）情景北礁華光礁浪花礁羚羊礁盤石嶼七連嶼西南深航島口門永興島西南玉琢礁0.060.130.640.710.170.250.130.240.591.23123456789100.750.980.860.880.741.010.490.820.711.030.640.830.570.70.550.810.190.410.761.070.150.180.750.850.130.20.150.260.30.430.080.140.860.870.050.280.110.150.480.921.211.762.182.771.642.221.231.62.323.871.622.451.071.371.632.451.732.571.622.491.462.762.493.321.472.661.752.852.824.70.570.770.670.780.540.740.210.330.590.66圖3影響西沙群島的TC強度分布特征（1991—2018年，TS熱帶風暴、STS強熱帶風暴、TY臺風、STY強臺風、SuperTY超強臺風）圖4影響西沙群島的TC路徑月分布特征（1991—2018年）影響西沙海域波浪分布特征的因素還有很多，特別是西沙群島不利于海水輻聚、輻散的海岸線，四面環海同時島嶼之間相互遮擋，與中沙群島毗鄰的特殊且復雜的地理位置產生的局部海水震蕩［15］。結論（1）西北型路徑是影響西沙群島的TC的最大可能路徑，該路徑TC最大可能等級是臺風，引起的波高主要集中在2.5~5.9m之間；西行型路徑是影響西沙群島的TC第二可能路徑，該路徑TC最可能等級是強臺風，引起的波高主要集中在4.0~5.9m之間；（2）西沙礁盤內風浪較小，基本小于1.0m，而西沙島嶼波浪掩護區內風浪較大，基本大于2.0m，甚至接近5m，一般永興島西南向掩護區波浪高于其他島嶼波浪掩護區；(3)無論西北型路徑還是西行型TC路徑，最大波高一般出現在最大風速之后，所以要注意提防大風過后的大浪災害。參考文獻：【相關文獻】劉桂海，劉松?我國近海波浪要素的長期分布[幾港口工程,1989,(6):25-37.楊憲章?長周期波的特性及對系泊船舶動態特性的影響[J]?港口工程,1989,(6):37-43.耿寶磊，文先華?臺風作用下瓊州海峽海域波浪特征分析[J].海洋工程,2013,31(6):59-67.丁千龍，曾永明,王江偉?西沙海區波浪特征分析與預報[J].海洋通報,2000,19(3):19-23.楊永增，孫玉娟,王關鎖，等?基于MASNUM海浪預報系統的北印度洋波浪特征模擬與預報分析[J].海洋科學進展,2011,29(1):1-9.楊永增，喬方利，趙偉，等?球坐標系下MASNUM海浪數值模式的建立及其應用[J].海洋學報,2005,27(2):1-7.韓淑宗，史玉姣?東中國海臺風浪分布特征研究[J].中國海洋大學學報,2013,43(10):1-7.BooijN,RisRC,HolthuijsenLH.Athird-generationwavemodelforcoastalregions:1.Modeldescriptionandvalidation[J].JournalofGeophysicalResearch:Oceans,1999,104(C4):7649-7666.Tolman