南通沿海臺風(fēng)風(fēng)暴潮分析及其經(jīng)驗預(yù)報初探

高清清，曹兵，高鑫鑫，徐常三

(國家海洋局南通海洋環(huán)境監(jiān)測中心站，江蘇南通226005)

1 引言

風(fēng)暴潮是來自海上的一種海洋災(zāi)害，是指由強烈大氣擾動如熱帶氣旋、溫帶氣旋等天氣系統(tǒng)所帶來的海面異常升降的現(xiàn)象。如果風(fēng)暴潮恰好與天文高潮相遇尤其是與天文大潮期間的高潮疊加，水位驟漲，將對沿海地區(qū)造成嚴(yán)重的財產(chǎn)損失，甚至對人身安全帶來嚴(yán)重威脅。因此，風(fēng)暴潮的預(yù)測與分析一直是相關(guān)領(lǐng)域研究人員的重要工作之一。按照產(chǎn)生風(fēng)暴潮現(xiàn)象的大氣擾動特征，將風(fēng)暴潮分為由熱帶氣旋和溫帶氣旋及寒潮大風(fēng)造成的兩大類，前者稱為臺風(fēng)風(fēng)暴潮，后者統(tǒng)稱為溫帶風(fēng)暴潮[1-2]。對于風(fēng)暴潮的研究早在20世紀(jì)20年代就已經(jīng)開展，至目前已經(jīng)取得了較多的成果[3-8]。

風(fēng)暴潮預(yù)報方法分為兩大類：一種是經(jīng)驗預(yù)報方法；另一種是數(shù)值預(yù)報方法。前者包括預(yù)報員的主觀經(jīng)驗和經(jīng)驗統(tǒng)計預(yù)報方法，后者包括諾模圖方法和數(shù)值模式預(yù)報方法[9]。葉鳳娟等[10]研究了單站資料在風(fēng)暴潮預(yù)報中的應(yīng)用，并指出單站資料可補充數(shù)值預(yù)報的不足，但也有一定的局限，與其他方法配合使用，才能更好地發(fā)揮其作用。曾德美[11]分析了青島港風(fēng)暴潮概況，并求取了青島港極值增減水的預(yù)報公式。謝亞力，黃世昌[12]將動力線性模型應(yīng)用于錢塘江河口臺風(fēng)風(fēng)暴潮經(jīng)驗預(yù)報中，通過合理選取預(yù)報量及預(yù)報因子，建立經(jīng)驗預(yù)報方程，取得較好的效果。

隨著計算機的發(fā)展和觀測技術(shù)的進步，風(fēng)暴潮的數(shù)值預(yù)報方法發(fā)展迅速。20世紀(jì)80年代以來，我國的風(fēng)暴潮數(shù)值模擬取得了很大的進展和應(yīng)用。在之后的風(fēng)暴潮研究工作中，大量的研究者將風(fēng)暴潮的數(shù)值預(yù)報作為風(fēng)暴潮的重要研究方向之一。朱建榮等[13]基于原正交ECOM模式，改進為非正交曲線坐標(biāo)系，應(yīng)用于長江河口、杭州灣及鄰近海區(qū)，改進后的模式可更好地應(yīng)用于河口海岸的研究中。端義宏等[14]利用三維河口、海岸模式（ECOMSi）建立了一個適用于長江口區(qū)風(fēng)暴潮的數(shù)值預(yù)報模式。

江蘇南通位于江蘇省東南部，居長江入海口北岸，瀕江臨海。該地區(qū)基本的氣候特點是氣候溫和，日照充分，四季分明，無霜期較長，海洋性氣候顯著。由于其獨特的地理位置，常受到西風(fēng)帶、副熱帶等天氣系統(tǒng)影響，氣象和氣象誘發(fā)災(zāi)害時有發(fā)生，其中主要的災(zāi)害之一就是風(fēng)暴潮災(zāi)害。南通幾乎每年夏秋季節(jié)都要受到臺風(fēng)影響，而由臺風(fēng)引起的風(fēng)暴潮災(zāi)害更是得到了南通市各級相關(guān)部門的極大關(guān)注。

本文選取南通沿海呂泗站的水文氣象資料，分析臺風(fēng)影響期間，該站的風(fēng)暴增水特征，并對該站臺風(fēng)風(fēng)暴潮經(jīng)驗預(yù)報公式進行初步探討，以期應(yīng)用于南通沿海臺風(fēng)風(fēng)暴潮預(yù)報業(yè)務(wù)工作中，更好地開展風(fēng)暴潮預(yù)報工作。

2 資料的選取

2000—2011年CMA-STI熱帶氣旋最佳路徑數(shù)據(jù)集；2000—2011年呂泗站實測潮位資料，天文潮資料；呂泗站地面氣象觀測資料中的風(fēng)速、風(fēng)向、氣壓等。

3 影響南通沿海的熱帶氣旋

3.1 基本概況

根據(jù)2000—2011年熱帶氣旋資料，共有13 個熱帶氣旋影響或者外圍影響南通沿海（見表1），近十年間沒有直接登陸南通的熱帶氣旋。

從表1 中還可看出，熱帶氣旋影響南通沿海多發(fā)生在7、8、9三個月份，6月、10月份也有出現(xiàn)。且近十年間，每年一般最多有2 個熱帶氣旋影響南通沿海。

3.2 熱帶氣旋分類

根據(jù)熱帶氣旋的登陸情況以及移動路徑，對上述進入南通或者影響南通的13 個熱帶氣旋進行分類，共得兩大類，分別為不在我國登陸的（Ⅰ、Ⅱ、Ⅲ型）以及在江浙一帶登陸的（Ⅳ、Ⅴ、Ⅵ型），見表2及圖1—6。

4 風(fēng)暴潮分析及預(yù)報

4.1 各類熱帶氣旋風(fēng)暴潮概況

表3列出了熱帶氣旋影響期間呂泗站的最大增水值以及影響時間段，可以看出，在浙江一帶登陸后由內(nèi)陸北上進入江蘇的熱帶氣旋對南通沿海造成的風(fēng)暴增水較其他類別的熱帶氣旋造成的增水要小。

結(jié)合圖1—6，發(fā)現(xiàn)Ⅰ型和Ⅱ型熱帶氣旋路徑較為相似，對南通沿海造成的風(fēng)暴增水也較為接近，除了1007 圓規(guī)以及1105 米雷造成的影響不大，較嚴(yán)重的風(fēng)暴潮基本上在125 cm 以上。而Ⅳ型和Ⅴ型這兩種臺風(fēng)路徑，雖然出海位置不同，但是差別并不明顯，且這兩種類型的臺風(fēng)造成的增水幅度較為接近，基本在50—90 cm。因此，將以上六種路徑的熱帶氣旋分為四類，分別為：第一類南通外海北上；第二類浙江外海轉(zhuǎn)向；第三類內(nèi)陸北上進入江蘇并從江蘇沿海出海；第四類在浙江一帶登陸后未在內(nèi)陸北上，很快又轉(zhuǎn)出海。

表1 2000—2011年影響或者外圍影響南通沿海的熱帶氣旋

表2 熱帶氣旋路徑分類

圖1 Ⅰ型.北上右轉(zhuǎn)向

圖2 Ⅱ型.外海北上

圖3 Ⅲ型.浙江外海轉(zhuǎn)向

圖4 Ⅳ型.登陸后進入江蘇，蘇北出海

圖5 Ⅴ型.登陸后進入江蘇，蘇南出海

圖6 Ⅵ型.登陸后不在內(nèi)陸北上，再轉(zhuǎn)入海上

表3 臺風(fēng)風(fēng)暴潮要素

表4 第一類熱帶氣旋影響期間，呂泗站風(fēng)暴潮情況

4.2 各類臺風(fēng)風(fēng)暴潮分析及預(yù)報公式

4.2.1 第一類南通外海北上

表4列出了南通外海北上的熱帶氣旋引起的風(fēng)暴增水情況，可以看出，此類熱帶氣旋影響期間，呂泗站的風(fēng)向前期均為東北東，隨著臺風(fēng)進一步北上，風(fēng)向轉(zhuǎn)為偏北或者西北西至偏北，最大增水出現(xiàn)時刻與最低氣壓出現(xiàn)的時刻較為接近。

選擇表4 中前四個熱帶氣旋作為研究分析對象，1109梅花過程用以對計算結(jié)果的檢驗。選取呂泗站增水值為預(yù)報量，臺風(fēng)過程該站最低氣壓與該月多年平均氣壓差、該站風(fēng)速、風(fēng)向以及臺風(fēng)中心距該站距離為預(yù)報因子。分別將臺風(fēng)影響過程的增水△H 與氣壓差△P、臺風(fēng)中心距呂泗站的直線距離L、風(fēng)速V 、風(fēng)向θ 之間做相關(guān)分析，發(fā)現(xiàn)風(fēng)暴增水△H 與上述物理量之間的相關(guān)性較好（通過95%信度檢驗）。

采用回歸分析方法，計算過程中發(fā)現(xiàn)L 對整個回歸方程的貢獻不大，因此得出最大增水△Hmax的預(yù)報方程為：

式（1）中：△P 為過程最低氣壓與該月多年平均氣壓之差的絕對值，單位hPa；Vmax、θ 分別為最大平均風(fēng)速以及對應(yīng)風(fēng)向；以下同。

對1109梅花影響期間的增水值與氣壓差、臺風(fēng)中心距呂泗站的直線距離、風(fēng)速分別進行相關(guān)分析，發(fā)現(xiàn)在該站風(fēng)向轉(zhuǎn)為西北向之前，相關(guān)性均很好，通過了99%信度檢驗。將梅花影響期間呂泗站測得的最低氣壓、最大平均風(fēng)速以及對應(yīng)風(fēng)向代入上式，計算得最大增水為152 cm，比實際最大增水值（139 cm）僅高出14 cm，誤差較小。

表5 第二類熱帶氣旋影響期間，呂泗站風(fēng)暴潮情況

4.2.2 第二類浙江外海轉(zhuǎn)向

表5 給出了浙江外海轉(zhuǎn)向的0014 臺風(fēng)桑美引起的風(fēng)暴增水情況。

根據(jù)其路徑，0014桑美于9月14日18時起轉(zhuǎn)為東北向，且移動速度明顯加快，隨后強度逐漸減弱。呂泗站最大增水發(fā)生在其轉(zhuǎn)向前8 小時左右。呂泗站增水△H 與氣壓差△P、距離L、風(fēng)速V 、風(fēng)向θ 分別進行相關(guān)分析，但是得到的相關(guān)性均不高，同時由于數(shù)據(jù)序列較少，回歸分析效果亦不明顯。

圖7 0014桑美影響期間呂泗站增水曲線

圖7 給出了“桑美”影響期間呂泗站增水曲線。可以看出，呂泗站增水過程有明顯的半日周期波動；天文潮和風(fēng)暴潮耦合作用[15]較強，周期內(nèi)增水最大值一般發(fā)生在天文高潮前3—4 h，因此建立此類臺風(fēng)暴潮與其他物理量之間的相互關(guān)系還需進一步深入探討，將在以后的工作中開展這方面的研究。

4.2.3 第三類內(nèi)陸北上進入江蘇并從江蘇沿海出海

表6 為第三類熱帶氣旋引起的風(fēng)暴增水情況。可以看出，此類熱帶氣旋引起的最大增水均出現(xiàn)在登陸之后。

選擇表6 中0509 麥莎、0515 卡努、0807 海鷗和0908莫拉克這四個臺風(fēng)過程為分析對象，0713韋帕用以檢驗計算的結(jié)果。

采用相關(guān)分析以及回歸分析方法，計算發(fā)現(xiàn)L對整個回歸方程的貢獻不大，因此在公式中未列出這一因子，第三類熱帶氣旋增水極值預(yù)報公式如下：

根據(jù)式（2）計算0713 韋帕引起的最大增水△Hmax=82 cm，實際最大增水為78 cm，結(jié)果較為理想。

4.2.4 第四類在浙江一帶登陸后未在內(nèi)陸北上，很快又轉(zhuǎn)出海

表6 第三類熱帶氣旋影響期間，呂泗站風(fēng)暴潮情況

表7 第四類熱帶氣旋影響期間，呂泗站風(fēng)暴潮情況

表7為登陸后又轉(zhuǎn)入海上的熱帶氣旋所引起的風(fēng)暴增水情況。此類熱帶氣旋影響過程中，呂泗風(fēng)向基本為東北向，平均風(fēng)暴增水較其他類型的熱帶氣旋大。

分析結(jié)果顯示，增水△H 與測站風(fēng)速V12、臺風(fēng)中心到呂泗站的距離L 相關(guān)性均較好。回歸分析得到增水極值預(yù)報形式為（3）式，由于此類臺風(fēng)個例較少，故有待進一步對其進行驗證。

式（3）中，V12為前12小時風(fēng)速。

5 小結(jié)和討論

綜上分析，近十年影響南通沿海的臺風(fēng)路徑可分為四類，分別為：南通外海北上；浙江外海轉(zhuǎn)向；內(nèi)陸北上進入江蘇并從江蘇沿海出海；在浙江一帶登陸后未在內(nèi)陸北上，很快又轉(zhuǎn)出海。臺風(fēng)暴潮由于不同臺風(fēng)路徑其變化亦有所差異。登陸后由內(nèi)陸進入江蘇的熱帶氣旋對南通沿海造成的風(fēng)暴增水較其他類別的熱帶氣旋造成的增水要小。

利用呂泗站水文氣象資料，分析了各類臺風(fēng)暴潮的特征，在此基礎(chǔ)上，運用相關(guān)分析、回歸分析等方法探討了風(fēng)暴潮極值預(yù)報公式，得出第一類、第三類和第四類的臺風(fēng)暴潮極值公式，其中第一類和第三類風(fēng)暴潮極值公式經(jīng)驗證后，效果較好；第四類臺風(fēng)暴潮公式由于個例有限，需要進一步加以驗證；結(jié)果還發(fā)現(xiàn)，第二類浙江外海轉(zhuǎn)向的臺風(fēng)對呂泗帶來的風(fēng)暴增水具有顯著的半日振蕩周期，研究需考慮風(fēng)暴潮與天文潮的耦合作用。

臺風(fēng)暴潮的過程值、極值與臺風(fēng)結(jié)構(gòu)、路徑、強度、大陸架形態(tài)和臺風(fēng)移動速度等物理量均有關(guān)。這些因子綜合作用于海水上，如何在預(yù)報中體現(xiàn)出來是較為困難的。本研究僅根據(jù)臺風(fēng)路徑來分類討論，建立了相應(yīng)風(fēng)暴增水極值預(yù)報方程。由于臺站資料的連續(xù)性不夠，本次研究工作中統(tǒng)計樣本總體偏少，在今后的工作中將把時間范圍擴大至30年，繼續(xù)改進預(yù)報方程，從而更好地做好臺風(fēng)暴潮預(yù)報工作。此外，對風(fēng)暴增水極值出現(xiàn)時刻的預(yù)報、風(fēng)暴潮與天文潮的耦合作用也將在后續(xù)的工作中進一步討論。

[1]馮士筰.風(fēng)暴潮導(dǎo)論[M].北京:科學(xué)出版社,1982.

[2]王喜年.關(guān)于溫帶風(fēng)暴潮[J].海洋預(yù)報,2005(S1):17-23.

[3]楊桂山.中國沿海風(fēng)暴潮災(zāi)害的歷史變化及未來趨向[J].自然災(zāi)害學(xué)報,2000,9(3):23-30.

[4]王華,姚圣康,龔茂珣,等.江蘇省洋口港風(fēng)暴潮數(shù)據(jù)分析[J].海洋通報,2007,26(3):26-32.

[5]胡德寶,龔茂珣,孔亞珍.強風(fēng)暴潮對上海地區(qū)影響研究[J].華東師范大學(xué)學(xué)報(自然科學(xué)版),2005(5/6):177-182.

[6]Mcinnes K L,Walsh K J E,Hubbert G D,et al.Impact of Sea-level Rise and Storm Surges on a Coastal Community[J]. Natural Hazards,2003,30(2):187-207.

[7]Zweers N C,Makin V K,de Vries J W,et al.On the influence of changes in the drag relation on surface wind speeds and storm surge forecasts[J].Natural Hazards,2012,62(2):207-219.

[8]Claudio T. Numerical study of non-linear storm-surges[J]. Pure and Applied Geophysics PAGEOPH,1973,111(1):2346-2358.

[9]王喜年. 風(fēng)暴潮預(yù)報知識講座第四講[J]. 海洋預(yù)報, 2001, 18(4):63-69.

[10]葉鳳娟,戴雪亮,鄒濤.單站資料在風(fēng)暴潮預(yù)報中的應(yīng)用[J].海洋預(yù)報,2001,18(4):38-42.

[11]曾德美.青島港風(fēng)暴潮經(jīng)驗統(tǒng)計預(yù)報[J].海洋預(yù)報,1992,9(3):66-73.

[12]謝亞力, 黃世昌. 錢塘江河口風(fēng)暴潮經(jīng)驗預(yù)報[J]. 海洋預(yù)報,2006,23(1):54-58.

[13]朱建榮,朱首賢.ECOM模式的改進及在長江河口、杭州灣及鄰近海區(qū)的應(yīng)用[J].海洋與湖沼,2003,34(4):364-373.

[14]端義宏,朱建榮,秦曾灝,等.一個高分辨率的長江口臺風(fēng)風(fēng)暴潮數(shù)值預(yù)報模式及其應(yīng)用[J].海洋學(xué)報(中文版),2005,27(03):11-19.

[15]趙永良,張延延,陳則實.黃海風(fēng)暴潮和天文潮非線性耦合作用的數(shù)值研究[J].海洋學(xué)報,1992,14(3):37-46.