李 穎, 方偉華, 林 偉, 葉妍婷

(1. 北京師范大學環(huán)境演變與自然災(zāi)害教育部重點實驗室, 北京 100875; 2. 北京師范大學地表過程與資源生態(tài)國家重點實驗室, 北京 100875; 3. 北京師范大學民政部-教育部減災(zāi)與應(yīng)急管理研究院, 北京 100875)

可能最大風暴潮風險評估中各等級熱帶氣旋設(shè)定方法

李 穎1,2,3, 方偉華1,2,3, 林 偉1,2,3, 葉妍婷1,2,3

(1. 北京師范大學環(huán)境演變與自然災(zāi)害教育部重點實驗室, 北京 100875; 2. 北京師范大學地表過程與資源生態(tài)國家重點實驗室, 北京 100875; 3. 北京師范大學民政部-教育部減災(zāi)與應(yīng)急管理研究院, 北京 100875)

可能最大熱帶氣旋的設(shè)定是可能最大風暴潮計算的基礎(chǔ), 對風暴潮災(zāi)害應(yīng)急疏散具有重要意義。利用1949～2011年中國氣象局(CMA)西北太平洋熱帶氣旋最佳路徑數(shù)據(jù)集、美國聯(lián)合臺風預(yù)警中心(JTWC)以及美國國家海洋和大氣管理局(NOAA)最大風速半徑數(shù)據(jù)集, 基于各等級熱帶氣旋參數(shù)之間的定量關(guān)系, 建立了各等級可能最大熱帶氣旋最大風速、中心氣壓、最大風速半徑、移動速度、移動方向等參數(shù)設(shè)定及路徑合成的方法。以福建省連江縣為例, 按照臺風、強臺風及超強臺風強度等級,分強度衰減和不衰減2種情況, 設(shè)定3種移動方向, 合成了共216場熱帶氣旋作為可能最大風暴潮的計算輸入。另外, 對參數(shù)敏感性、風場參數(shù)設(shè)定、參數(shù)設(shè)定與計算量的關(guān)系、疊加天文潮以及潰堤等問題進行了討論。

熱帶氣旋; 可能最大風暴潮; 參數(shù)設(shè)定; 風險評估

風暴潮災(zāi)害風險評估是風暴潮災(zāi)害管理的基礎(chǔ)。我國瀕臨西北太平洋, 登陸我國的熱帶氣旋頻率高、影響范圍廣、突發(fā)性強, 其引發(fā)的熱帶風暴潮災(zāi)害經(jīng)常造成嚴重的人員傷亡和財產(chǎn)損失[1-3]。根據(jù)評估對象、方法及結(jié)果的不同, 風暴潮風險評估包括[4]: (1)某點各多年一遇水平下風暴潮增水或水深評估, (2)某區(qū)域各多年一遇水平下風暴潮淹沒水深及范圍評估, (3)不同等級熱帶氣旋下可能最大風暴潮淹沒范圍及水深評估, (4)風險評估與區(qū)劃等。其中, 不同多年一遇水平評估產(chǎn)品可為工程設(shè)計、各類規(guī)劃及保險產(chǎn)品設(shè)計等服務(wù)[5], 而可能最大風暴潮則主要面向于熱帶風暴潮災(zāi)害的應(yīng)急指揮、疏散避難、社區(qū)綜合減災(zāi)等應(yīng)用[6]。

可能最大熱帶氣旋(Probable Maximum Tropical Cyclone, PMTC)是一種假想的平穩(wěn)狀態(tài)的熱帶氣旋,它是根據(jù)可以在特定海岸地區(qū)發(fā)生最大持續(xù)風速所選擇的氣象參數(shù)值的組合[7]。其中, 平穩(wěn)狀態(tài)是指熱帶氣旋各氣象參數(shù)相互約束且相互適應(yīng)的狀態(tài)。可能最大風暴潮(Probable Maximum Storm Surge, PMSS)是 PMTC沿最不利途徑逼近海岸的情景下產(chǎn)生的沿海各點的最大風暴潮增水及淹沒。目前, 對于如何設(shè)定PMTC以計算PMSS, 已開展了一系列研究。

美國對可能最大颶風(Probable Maximum Hurricane, PMH)及其引起的 PMSS的研究開展早、理論強、應(yīng)用廣。1959年, 美國陸軍工程兵團(USACE)在美國納拉干海灣與新奧爾良地區(qū)開展了 PMH研究[8],但僅重點考慮了中心氣壓; 隨后, 出現(xiàn)了越來越多的 PMH理論方法探索及區(qū)域應(yīng)用研究, 例如, Wall采用概率性方法設(shè)定可能最大颶風參數(shù), 并依此評估了紐約杰弗森港口的PMSS增水[9]; 美國核安全導則規(guī)定, 設(shè)定 PMSS時應(yīng)在可能最大颶風引起的增水基礎(chǔ)上疊加10%的超高潮位并考慮海平面異常[10];美國國家海洋和大氣管理局(NOAA)技術(shù)報告NWS23中提供了可能最大颶風的參數(shù)范圍, 通過設(shè)定多種情景來組合颶風參數(shù), 從而評估PMSS[11]。

2007年, 美國國家颶風中心(National Hurricane Center, NHC)在SLOSH Display Program(SDP)項目的支持下建立了最大可能增水產(chǎn)品(Maximum Envelope of Water, MEOW)[6]。目前NHC已經(jīng)在美國37個海區(qū)完成了 80多個 MEOW 產(chǎn)品的計算及應(yīng)用, MEOW結(jié)果中某個點和其相鄰點的最大增水可能是由不同合成路徑的颶風產(chǎn)生的。用戶雖然不知道是哪場颶風計算的最大增水, 但是MEOW提供了1種在已知颶風強度等級前提下, 綜合考慮颶風路徑預(yù)報不確定性、天文潮位疊加不確定性的最壞情景下的最大可能增水。

2003年, 日本內(nèi)閣府聯(lián)合其他風暴潮災(zāi)害管理部門, 在多個研究機構(gòu)支持下, 編制發(fā)布了《海嘯-風暴潮危險制圖手冊》, 對 PMTC參數(shù)的設(shè)定方法進行了說明[12]。具體方法可分為2類, 一是參照影響研究區(qū)的最嚴重風暴潮災(zāi)害的歷史熱帶氣旋; 二是參照整個日本范圍內(nèi)造成過嚴重風暴潮災(zāi)害的熱帶氣旋, 一般選定1959年的伊勢灣臺風作為參照。另外, 對較嚴重風暴潮的潰堤情景進行了粗略規(guī)范。

1991年, 中國國家核安全局發(fā)布了《核安全導則匯編》, 在核電站設(shè)計基準熱帶氣旋相關(guān)部分給出了PMTC參數(shù)的設(shè)定方法[6], 隨后, 其他研究針對廣東陽江、江蘇連云港等核電站附近或其他熱帶風暴潮易發(fā)地區(qū)PMTC相關(guān)水文氣象參數(shù)進行了設(shè)定和計算[13-15]。

雖然美、日、中三國均有一些PMTC方面的研究, 但也存在著較大的差別: (1)美國在颶風路徑參數(shù)設(shè)定上考慮的情景比較全面, 主要有三個特點:一是考慮了不同颶風等級分別設(shè)定; 二是考慮不同情景進行參數(shù)組合, 并能計算出各種組合的概率;三是采用風暴潮數(shù)值模式對颶風參數(shù)進行敏感性分析, 提高了設(shè)定結(jié)果的可靠性。(2)日本的PMTC設(shè)定方法有2個特點: 一是對于最大風速、最低氣壓、最大風速半徑、登陸方向、移動速度等參數(shù)的具體設(shè)定方法沒有做出明確的規(guī)定, 而是由地方政府依據(jù)本地區(qū)情況進行靈活設(shè)定; 二是沒有按照不同熱帶氣旋等級分別計算制圖, 只是選擇本地或者全國最壞的歷史情景進行設(shè)定, 部分原因是日本地處太平洋地震帶, 在全國很多地區(qū)海嘯造成的淹沒范圍和水位要明顯高于熱帶風暴潮。(3)中國的《核安全導則匯編》中的PMTC參數(shù)設(shè)定, 具有2大特點: 一是不考慮不同熱帶氣旋等級分別進行設(shè)定; 二是參數(shù)設(shè)定為最不利情況的組合, 出現(xiàn)這種極端組合的概率非常小。

可以看出, 我國以往的PMTC和PMSS研究, 目的是為了核工程的絕對安全, 其參數(shù)設(shè)定和計算沒有針對不同等級的熱帶氣旋。目前, 我國地方政府越來越注重以人為本科學應(yīng)對災(zāi)害, 需要大范圍地疏散轉(zhuǎn)移沿岸居民并強制漁船進港、漁民上岸, 但在確定疏散規(guī)模時一些地區(qū)存在過度疏散情況, 易導致災(zāi)害應(yīng)對社會成本過高和行政資源浪費。2011年“3.11”日本地震海嘯后, 國家海洋局開始組織開展國家、省、市、縣等尺度的風暴潮災(zāi)害風險評估和區(qū)劃工作, 現(xiàn)已編制了《風暴潮災(zāi)害風險評估與區(qū)劃技術(shù)導則(試行)》[5], 其中提出要對我國沿海各縣進行不同熱帶氣旋等級的 PMSS淹沒范圍和水深評估,為地方科學合理進行風暴潮災(zāi)害應(yīng)急避難、轉(zhuǎn)移疏散提供依據(jù)。

風暴潮增水大小通常與熱帶氣旋強度、范圍、移速、移向及下墊面情況有關(guān)。一般而言, 最大風速越大、最大風速半徑越廣、登陸時移動方向越垂直于海岸, 則風暴潮增水越大; 中心最低氣壓對增水的直接影響則較小, 氣壓越低增水越大; 通常移動速度越小, 海灣和河口地區(qū)增水越高、淹沒范圍越廣,而移動速度越大則會使開闊海岸地區(qū)增水越高[16-17]。另外, 局地的岸線形狀、海灣、河流、海岬及島嶼等地形作用、天文潮位、風場計算中的Holland B系數(shù)[18]、潰堤情景等也會影響PMSS的計算。

據(jù)此, 本文將結(jié)合中國沿海熱帶氣旋及熱帶風暴潮的特點, 綜合考慮縣尺度熱帶風暴潮風險評估所需的數(shù)據(jù)及技術(shù)要求, 以福建省連江縣為例, 探索適合中國沿?？h尺度可能最大熱帶風暴潮風險評估中的熱帶氣旋熱帶參數(shù)設(shè)定方法, 并給出具體的設(shè)定結(jié)果, 為開展不同熱帶氣旋等級的可能最大風暴潮淹沒范圍和水深評估工作提供參考。

1 資料處理

本文主要采用了中國氣象局(CMA)熱帶氣旋最佳路徑數(shù)據(jù)集(1949～2011年), 其記錄了西北太平洋熱帶氣旋每6 h的時間(T)、中心經(jīng)度(Dlon)、中心緯度(Dlat)、中心最低氣壓(Pmin)以及平均近中心最大風速(Vmax)等指標[19], 但是該數(shù)據(jù)集中未包含對熱帶氣旋設(shè)定具有重要作用的登陸點相關(guān)參數(shù)。為此, 本文利用登陸前后的數(shù)據(jù)進行了插值。首先, 建立亞歐大陸較為光滑的海岸線作為熱帶氣旋是否登陸的分界線, 避免計算過程中由于實際海岸線過于復(fù)雜造成短時間內(nèi)多次登陸的問題;其次, 利用登陸前后的參數(shù), 從線性、拋物線、樣條、三次樣條及分形等函數(shù)形式中優(yōu)選出不同參數(shù)的最佳插值方法, 建立了自適應(yīng)插值方法, 對登陸點參數(shù)T、Dlon、Dlat、Pmin、Vmax等進行插值, 并計算了移動速度(Vt)和移動方向(θ), 詳見文獻[20]。

由于最大風速半徑Rm實測數(shù)據(jù)較少, 本文還使用了兩組包含Rm的數(shù)據(jù)集。一是美國聯(lián)合臺風預(yù)警中心(JTWC)最佳路徑數(shù)據(jù)集(2001～2011年)[21], 二是NOAA西北太平洋飛機探測臺風資料(1960～1974年), 記錄了173個熱帶氣旋的Pmin及Rm數(shù)據(jù)[11]。所有數(shù)據(jù)集在使用前均將單位轉(zhuǎn)換為與CMA一致。其中, 由于CMA數(shù)據(jù)集中Vmax為2 min最大平均風速, 而JTWC中為1 min最大平均風速。為此, 參照世界氣象組織標準及已有研究[22-23], 對2 min與1 min的平均風速進行了轉(zhuǎn)換,轉(zhuǎn)換關(guān)系為2 min最大平均風速等于1 min最大平均風速減去 0.91 m/s。除特別指出外, 下文Vmax均指2 min最大平均風速。

CMA、JTWC和NOAA是西北太平洋熱帶氣旋歷史數(shù)據(jù)收集最為權(quán)威的機構(gòu), 在整編上述數(shù)據(jù)集中均進行了數(shù)據(jù)質(zhì)量及一致性檢驗。

2 參數(shù)設(shè)定方法

根據(jù)PMSS與PMTC參數(shù)之間的關(guān)系, 為了達到PMSS要求, 確定PMTC參數(shù)設(shè)計步驟如下: (1)基于熱帶氣旋等級標準確定各等級PMTC的最大風速; (2)根據(jù)中心氣壓、最大風速半徑與最大風速的經(jīng)驗統(tǒng)計關(guān)系, 結(jié)合已設(shè)定的最大風速, 在約束條件下盡可能提高中心氣壓差、最大風速半徑的設(shè)定值; (3)鑒于移動速度對 PMSS的影響與熱帶氣旋路徑所經(jīng)歷的局地地形有很大關(guān)系, 分情況給出參考建議; (4)設(shè)定與計算量關(guān)系較為密切的指標大小,包括移動方向多少、路徑間距的大小、是否考慮衰減等, 在計算能力滿足的情況下, 盡量增加多種移動方向、縮小路徑間距及考慮衰減情況。基于歷史觀測數(shù)據(jù), 分析熱帶氣旋主要參數(shù)(Dlon、Dlat、Vmax、Pmin、Rm、Vt、θ)(圖 1)之間的定量關(guān)系, 結(jié)合其對風暴潮增水的影響[16-17], 進行PMTC參數(shù)組合設(shè)定。

2.1 登陸參考點

為了對研究區(qū)內(nèi)及附近的歷史熱帶氣旋參數(shù)進行特征統(tǒng)計, 或者對合成路徑的參數(shù)進行設(shè)定, 需要確定一個基準點作為熱帶氣旋登陸參考點(Central Reference Point, CREF)。CREF一般選擇位于研究區(qū)光滑后的海岸線的中間位置, 并根據(jù)沿海海岸線及島嶼分布情況進行適當調(diào)整。圖 1給出了中國沿海光滑處理后的海岸線、沿海240個市(縣)行政邊界以及各市(縣)CREF分布。其中, 大陸沿海203個市(縣), CREF共209個; 海南省沿海13個市(縣), CREF共13個; 臺灣省沿海24個市(縣), CREF共24個。

圖 1 中國沿海各市(縣)CREF分布、福建連江縣位置及PMTC參數(shù)示意圖Fig. 1 CREF for China's coastal cities (counties), location of Lianjiang County, Fujian Province and schematic diagram of PMTC parameters

2.2 最大風速

熱帶氣旋其他參數(shù)不變時,Vmax越大引起的風暴潮增水越大, 因而, 可以將各等級 PMTC的Vmax設(shè)為該等級Vmax的上限值, 即熱帶低壓17.1 m/s, 熱帶風暴24.4 m/s, 強熱帶風暴32.6 m/s, 臺風41.4 m/s,強臺風50.9 m/s[24]。

由于超強臺風Vmax沒有上限, 其設(shè)定應(yīng)考慮PMSS風險評估成果的應(yīng)用特點, 并結(jié)合研究區(qū)及附近歷史熱帶氣旋的Vmax來綜合判定。例如, 一是考慮西北太平洋所有歷史熱帶氣旋登陸點的Vmax的最大值(1949～2011年為54.03 m/s); 二是考慮研究區(qū)CREF附近范圍內(nèi)歷史熱帶氣旋登陸點的Vmax的最大值; 三是考慮CREF附近不同重現(xiàn)期(如200 a一遇)下的Vmax。

2.3 中心氣壓差

首先, 找出Vmax與Pmin之間的經(jīng)驗關(guān)系。熱帶氣旋中心壓差PD是外圍氣壓Pn與Pmin的差值。PD與Vmax有著很強的相關(guān)性, Atkinson和Holliday[25]于1977年對西北太平洋熱帶氣旋的風壓關(guān)系進行了擬合(公式(1)), 其中A、B為經(jīng)驗參數(shù)。另外, 在探討風壓關(guān)系時也有研究考慮了緯度、Vt等因素的影響[23,26]。

本文利用 CMA數(shù)據(jù)集所有樣本對風壓關(guān)系進行了擬合, 如公式(2)所示, 其中,R為相關(guān)系數(shù), 通過Vmax推算出PD(圖2)。

其次, 比較登陸點PD的歷史觀測值與平均經(jīng)驗值之間的關(guān)系。分析1949～2011年登陸點風壓資料,可以看出, 平均經(jīng)驗值加上歷史觀測值的 1倍標準差(20.77)后, 結(jié)果不僅一般均大于歷史登陸點觀測值的最高值, 而且正好位于其附近(圖2)。

圖2 西北太平洋熱帶氣旋PD與Vmax關(guān)系圖Fig. 2 Relationship betweenPD and Vmax for cyclones in the Western North Pacific

最后, 利用歷史觀測值與平均經(jīng)驗值之間的關(guān)系, 設(shè)定Pmin, 如公式(3)所示, 其中SD為 1倍方差,Pn取為1 008 hPa。

2.4 最大風速半徑

首先, 繪制Rm與Vmax之間的散點圖。通常采用經(jīng)驗擬合法利用Vmax估算Rm, 也有一些研究采用逐步回歸法利用PD及緯度等參數(shù)估算Rm[27-28]。本文對Vmax和Rm之間的關(guān)系進行了探索, 總的規(guī)律為Vmax越大時Rm越小(圖3), 但同時可以看出相關(guān)性非常弱。

其次, 找出Rm的最大值與Vmax之間的經(jīng)驗關(guān)系。設(shè)定PMTC的Rm時, 重點研究最大值而非平均值。由圖4可以看出, 熱帶氣旋Rm的最大值與Vmax具有較好的線性關(guān)系(公式(4)):

然后, 比較Rm最大值的歷史觀測值與平均經(jīng)驗值之間的關(guān)系。利用擬合曲線加1倍標準差(21.37 km),同DP類似, 其值也均高于歷史觀測值且位于其附近(圖4)。

最后, 利用公式(5)設(shè)定PMTC登陸點的Rm。

圖3 西北太平洋熱帶氣旋Rm與Vmax散點圖Fig. 3 Relationship between Rm and Vmax for tropical cyclones in the Western North Pacific

圖4 西北太平洋熱帶氣旋Rm的最大值與Vmax關(guān)系圖Fig. 4 Relationship between maximum Rm and Vmax for tropical cyclones in the Western North Pacific

2.5 移動速度

首先, 探索Vt與強度、緯度等參數(shù)的關(guān)系。根據(jù)CMA資料統(tǒng)計,Vt與Vmax、Pmin等強度變量關(guān)系不大?？紤]到Vt受地轉(zhuǎn)偏向力和大氣環(huán)流影響較大,與其本身所處緯度關(guān)系密切[29-30]。由此, 本文對以5°為間隔的各緯度帶內(nèi)登陸點的平均移速進行了統(tǒng)計(表1)。

表1 西北太平洋熱帶氣旋登陸點Vt分緯度統(tǒng)計結(jié)果Tab. 1 Statistical results of translating velocity Vt of landing tropical cyclones in the Western North Pacific indifferent latitudes

其次, 參考各緯度帶內(nèi)Vt的統(tǒng)計結(jié)果, 設(shè)定登陸點。移動速度大, 可造成較大的風暴增水極值, 而移動速度小持續(xù)時間較長的熱帶氣旋, 也可能導致較大的淹沒范圍?？紤]到Vt本身波動較大, 為減少PMSS估算的不確定性, 特定區(qū)域PMTC的Vt應(yīng)通過敏感性分析來確定。但為便于實際應(yīng)用, 一些研究根據(jù)PD設(shè)定了Vt上下限[7,11], 或?qū)t設(shè)為區(qū)域樣本最大值[13]或者平均值[15,31]。本文案例中將PMTC的Vt設(shè)為研究區(qū)歷史記錄的平均值。

2.6 移動方向

首先, 探索登陸時θ與強度等參數(shù)間的關(guān)系。θ與Vmax、Rm等參數(shù)沒有顯著相關(guān)關(guān)系, 而更多受引導氣流的影響[32]。因此設(shè)定θ時, 不考慮其與熱帶氣旋其他參數(shù)之間的關(guān)系。

其次, 對熱帶氣旋登陸θ進行頻率分析。由于歷史觀測樣本較少, 本文采用 Chouinard方法分析[33],實質(zhì)為以高斯核函數(shù)對在距 CREF點一定空間范圍內(nèi)的歷史樣本進行插值, 從而得到各方向熱帶氣旋登陸頻次(公式(6)、(7)、(8))。其中,λ為熱帶氣旋頻次, 表示一定方向(如1°)、一定距離(如1 km)、一定時間(如圖1 a)范圍內(nèi)熱帶氣旋登陸的次數(shù)。w為權(quán)重, 是距離熱帶氣旋中心的距離和角度的函數(shù)。θi為第i個樣本點的移動方向,hd為最優(yōu)距離窗寬,hθ為最優(yōu)角度窗寬。

然后, 分析最有利于造成研究區(qū)內(nèi)風暴潮增水的θ。根據(jù)歷史風暴潮資料, 或者采用數(shù)值模式計算,對熱帶風暴潮增水對登陸時θ的敏感性進行分析,確定最有利于風暴潮增水的θ。

最后, 確定θ間隔并設(shè)定登陸θ列表。綜合考慮登陸θ頻率分布、最有利風暴潮增水的θ后, 既可以采用等間距設(shè)定θ, 也可以在增水對登陸θ比較敏感的角度多設(shè)定一些方向。

2.7 路徑間距

為保證最大風暴潮淹沒水深空間分布的合理性,首先, 所設(shè)定的路徑間距要足夠小, 才能保證熱帶氣旋路徑對模擬區(qū)域的的平滑覆蓋, 若路徑間隔太遠, 則可能導致風暴潮增水估算結(jié)果波動較大, 從而低估模擬風暴潮的峰值增水出現(xiàn)的可能性。其次,若路徑間距設(shè)定過小則會導致模擬次數(shù)過多, 帶來較大的計算成本, 影響PMSS風險評估效率。已有研究表明, 當路徑間距設(shè)定為Rm時, 已能得到較好的評估結(jié)果[32]。另外, 由于熱帶氣旋的右側(cè)風速更大,所引起的風暴潮更強, 因而主路徑方向左側(cè)設(shè)定的路徑數(shù)量一般應(yīng)略多于右側(cè), 若研究區(qū)附近熱帶氣旋不對稱性越高, 則兩側(cè)路徑數(shù)量差異越大。

因此, 應(yīng)綜合考慮敏感性及計算能力進行路徑間距設(shè)定, 路徑間距設(shè)定值應(yīng)不大于Rm, 在計算能力允許的條件下, 可適當減小路徑間距設(shè)定值。

2.8 強度衰減

若假設(shè)熱帶氣旋登陸前后無衰減, 對 PMTC各路徑點參數(shù)進行簡化設(shè)定, 則各路徑點Pmin、Vmax、Rm等特征參數(shù)可設(shè)定為與CREF一致。

實際熱帶氣旋登陸前后受陸地影響, 其強度會從登陸前一定距離開始衰減。如果考慮衰減, 可根據(jù)研究區(qū)歷史熱帶氣旋實際情況, 分洋面衰減和陸地衰減兩部分進行設(shè)定:

首先, 設(shè)定洋面強度衰減系數(shù)。熱帶氣旋在登陸前一個時段,PD、Rm均呈較明顯的線性變化趨勢[32]。例如, 美國墨西哥灣的研究表明, 距離海岸約為90 n mile時的Rm約為登陸時Rm的1/1.3倍, 而PD的衰減系數(shù)與登陸點Rm有關(guān)[32]。由于離岸較遠的強度變化與當時的環(huán)境有關(guān), 具有一定的隨機性,可設(shè)定為最多增強至某一強度(如Vmax約70 m/s)處即可。

其次, 設(shè)定登陸強度衰減系數(shù)。熱帶氣旋登陸后,PD隨時間呈指數(shù)衰減[34-35], 其衰減規(guī)律近似符合以下關(guān)系(公式(9)), 其中,PD0為登陸點中心氣壓差;T為登陸后時間;a、c為常參數(shù)。

最后, 設(shè)定各路徑點強度參數(shù)。根據(jù)洋面衰減系數(shù)及登陸后衰減系數(shù), 計算每個時間間隔的PD衰減量, 從登陸點開始推算各路徑點的Pmin、Vmax、Rm等。

2.9 路徑合成

在熱帶氣旋參數(shù)確定后, 再通過設(shè)定起始點及終止點位置、時間分辨率等, 即可合成整條路徑。具體設(shè)定步驟如下:

首先, 設(shè)定起始點位置。起始點設(shè)定于登陸點前移路徑的相反方向, 其與登陸點的距離由選用的風暴潮模式的計算能力決定, 通常應(yīng)不小于 90 n mile[32],且設(shè)定時要覆蓋風暴潮數(shù)值模式計算的網(wǎng)格范圍。

其次, 設(shè)定終止點位置。熱帶氣旋登陸后對沿岸風暴潮的影響逐漸減小, 終止點處的熱帶氣旋所引起的風暴增水應(yīng)足夠小, 例如可以設(shè)定為登陸后300 km終止。

然后, 設(shè)定路徑點時間間隔。一般情況下, 可與最佳路徑集中一致設(shè)定為6 h, 也可以設(shè)定為更密的時間間隔。

最后, 組合所有參數(shù)。根據(jù)熱帶氣旋Vt和θ, 計算每個時間間隔的移動距離, 根據(jù)各參數(shù)隨與登陸點距離的變化規(guī)律, 從登陸點組合參數(shù)推算出整條合成路徑的所有組合參數(shù)。

為了自動組合參數(shù)形成合成路徑, 本文開發(fā)了一個自動化程序, 輸入CREF、登陸參數(shù)、路徑點時間間隔、起始點及終止點位置, 即可以GIS格式或者文本格式輸出所有合成路徑[36]。

3 參數(shù)設(shè)定示例

福建省受熱帶氣旋災(zāi)害嚴重, 2001～2010年間因熱帶氣旋災(zāi)害造成的直接經(jīng)濟損失達606.23億元[37]。福州市連江縣是福建東部沿海受熱帶氣旋及風暴潮影響較大的縣(圖1)。本文將以連江縣為例, 給出PMTC的參考設(shè)定。

1) 確定設(shè)定強度等級及研究區(qū) CREF。由于連江縣熱帶低壓、熱帶風暴和強熱帶風暴所引起的風暴潮增水相對較小, 本文僅以臺風、強臺風、超強臺風三個等級為例設(shè)定參數(shù)。連江縣海岸線長、島嶼眾多, 結(jié)合熱帶氣旋主要登陸方向, 近似取其東南部海岸線中點(119.7°E, 26.28°N)作為CREF。

2) 設(shè)定PMTC登陸時的參數(shù)。

臺風、強臺風的Vmax、Pmin、Rm可依據(jù)前文進行設(shè)定, 另外, 綜合考慮西北太平洋登陸點Vmax歷史最大值、連江縣附近歷史最大值(50.13 m/s), 以及連江縣不同重現(xiàn)期下的Vmax(圖5), 連江縣超強臺風的Vmax設(shè)定為60 m/s。由Vmax進而設(shè)定Pmin及Rm(表2)。

圖5 連江縣熱帶氣旋不同重現(xiàn)期下Vmax曲線Fig. 5 Vmax of tropical cyclones at different return periods for Lianjiang County

考慮連江縣特殊的地理位置, 一方面由于熱帶氣旋主方向為-65°(圖 6), 南岸易受風暴潮侵襲, 另一方面, 一旦熱帶氣旋從東北方向來襲, 恰使潮水涌入羅源灣, 則會造成較大的增水和淹沒。綜合考慮, 設(shè)定PMTC的登陸θ為以歷史熱帶氣旋主方向為基準, 以 60°為間隔, 即θ可取為–125°, –65°和–5°。另外, 參照25°～30°N緯度帶內(nèi)登陸熱帶氣旋的Vt歷史平均值, 設(shè)定Vt為21.29 km/h。

圖6 1949～2011年連江縣CREF附近140 km內(nèi)熱帶氣旋登陸θ頻次分布圖Fig. 6 Directional rate of translating direction θ for tropical cyclones landing near the CREF(with in 140 km) in Lianjiang County

3) 分別針對無強度衰減和強度衰減兩種情況,確定各路徑點Pmin、Vmax和Rm。

不考慮強度衰減情況下, 各路徑點參數(shù)設(shè)定與登陸時設(shè)定一致。

考慮強度衰減情況進行參數(shù)設(shè)定。首先, 設(shè)定Pmin。登陸前, 結(jié)合登陸點Rm, 各等級下Pmin變化速率(即登陸前靠近CREF1 km所增加的氣壓值)分別為:臺風0.2603 hPa/km, 強臺風0.2688 hPa/km, 超強臺風0.2769 hPa/km, 設(shè)定由CREF向海上路徑點按此速率線性增強至884.56 hPa(相當于Vmax達70 m/s), 然后強度保持不變。登陸后, 假設(shè)強度依照歷史平均衰減速率衰減, 具體地, 基于連江縣附近歷史登陸點樣本,利用公式(9)擬合每場熱帶氣旋參數(shù), 進而計算a、c平均值, 分別為0.022 172 8和0.689 49。其次, 結(jié)合公式(3), 根據(jù)Pmin反算Vmax。最后, 設(shè)定Rm。各等級下Rm線性變化速率表示登陸前靠近或登陸后遠離登陸點1 km時, 熱帶氣旋最大風速半徑所減小的量值(km)的結(jié)果分別為臺風為0.548 7, 強臺風為0.564 4, 超強臺風為0.579 5。

4) 設(shè)定合成路徑的時間精度和空間密度, 合成路徑。設(shè)定沿移動方向的反方向距CREF約600 km處為起始點, 登陸后約300 km處路徑終止; 考慮強度無衰減情況和強度衰減兩種情況; 設(shè)定時間間隔為1 h、路徑間距為23 km(約為平均Rm一半); 對每一種參數(shù)組合, 設(shè)定路徑12條, 其中主路徑1條、主路徑左側(cè)6條、右側(cè)5條。對臺風、強臺風及超強臺風等級下共216條PMTC路徑進行合成(圖7)。以臺風等級為例, 在不考慮衰減的情況下, 計算圖 7中無衰減情況下的 36場熱帶氣旋的風暴潮淹沒范圍及水深, 則PMSS淹沒范圍取值為36場風暴潮的所有淹沒范圍, 水深取值為36場風暴潮中的最大水深。

4 結(jié)論及討論

4.1 結(jié)論

PMTC的設(shè)定是PMSS計算的基礎(chǔ), 但目前國內(nèi)不同等級PMTC設(shè)定方法研究較少。基于1949～2011年間影響中國的熱帶氣旋關(guān)鍵參數(shù)統(tǒng)計, 探討了可能最大熱帶氣旋參數(shù)設(shè)定方法, 主要結(jié)論如下:

1) 基于各等級熱帶氣旋參數(shù)之間的定量關(guān)系,建立了各等級可能最大熱帶氣旋的最大風速、中心氣壓、最大風速半徑、移動速度、移動方向等參數(shù)設(shè)定及路徑合成方法。

2) 設(shè)定方法在福建省連江縣臺風、強臺風及超強臺風等級下合成路徑結(jié)果表明, 該方法可行有效。

4.2 討論

PMSS的計算還需要考慮一些其他問題, 討論如下:

1) 主要熱帶氣旋參數(shù)敏感性分析。在各等級PMTC給定Vmax、Pmin等強度變量下, 不同熱帶氣旋的Vt及θ對風暴潮增水影響差異較大[17], 同時, 不同研究區(qū)風暴潮增水也會受當?shù)睾０毒€分布、地形等因素影響而呈現(xiàn)異質(zhì)性。因此, 采用風暴潮數(shù)值模式, 針對研究區(qū)開展風暴潮增水對 PMTC的參數(shù)敏感性試驗, 確定區(qū)域最適參數(shù)是下一步 PMSS研究重點。

2) 參數(shù)設(shè)定與計算量關(guān)系。對熱帶氣旋強度等級及各種相關(guān)參數(shù)設(shè)定時應(yīng)考慮計算量, 例如, 減小路徑間隔及θ間隔、設(shè)定更遠的起算和終止距離會有助于 PMSS的計算結(jié)果精度的提高, 但是一方面會增加熱帶氣旋場次, 另一方面會增加熱帶氣旋持續(xù)時間, 從而增加風暴潮數(shù)值模式的總計算量,可綜合衡量評估精度需求及計算能力設(shè)定參數(shù)。

3) 熱帶氣旋風場參數(shù)設(shè)定。表征風速徑向分布形狀的Holland B系數(shù)等對熱帶氣旋風場以及風暴潮增水的計算有重要影響, 在計算 PMSS時盡量選擇能反映Holland B系數(shù)特點的經(jīng)驗風場模型。

4) PMSS計算中的其他參數(shù)設(shè)定。實際潮位不僅與 PMTC引起的增水有關(guān), 還與天文潮位相關(guān), 在實際計算時需考慮天文潮, 另外還應(yīng)考慮沿海防御能力并進行潰堤情景設(shè)定等。

致謝: 國志興、譚駿、吳維登、陳延盛等參加了可能最大風暴潮計算的相關(guān)討論, 王喜年提供了部分最大風速半徑數(shù)據(jù), 吳小寧協(xié)助了部分制圖工作, 兩位匿名評審人給出了有益的修改建議, 特表感謝。

[1] 王喜年. 風暴潮災(zāi)害及其地理分布[J]. 海洋預(yù)報, 2001, 18(2): 70-77.

[2] 王喜年. 風暴潮風險分析與計算[J]. 海洋預(yù)報, 2002, 19(4): 73-76.

[3] 中國氣象局. 中國氣象災(zāi)害年鑒(2010)[M]. 北京:氣象出版社, 2010.

[4] Resio D T, Irish J, Cialone M. A surge response function approach to coastal hazard assessment – Part 1: Basic concepts[J]. Nat Hazards, 2009, 51: 163-182.

[5] 國家海洋局, 風暴潮災(zāi)害風險評估和區(qū)劃技術(shù)導則(試行)[R]. 北京: 國家海洋局, 2012.

[6] Glahn B, Taylor A, Kurkowski N, et al. The role of the slosh model in national weather service storm surge forecasting[J]. National Weather Digest, 2009, 33(1): 3-14.

[7] HAD101/11. 核電廠設(shè)計基準熱帶氣旋[S]. 國家核安全局, 1991.

[8] Harris D L. An interim hurricane storm surge guide[R]. Washington: Weather Bureau Department of Commerce, 1959.

[9] Wall. Probabilistic assessment of flooding hazard for nuclear power plants[J]. Nuclear Safe, 1974, 15(4): 399-408.

[10] Regulatory Guide 1.59-1977, Design basis flood for nuclear power plants (revision 2)[S]. Washington D.C.: US Nuclear Regulatory Commission.

[11] Schwerdt R W. Meteorological criteria for standard project hurricane and probable maximum hurricane windfields, Gulf and East coasts of the United States[R]. Washington D.C.: US Department of Commerce, 1979.

[12] 津波? 高潮ハザ ー ドマップ研究會事 務(wù)局. 津波?高潮ハザ ードマップマニュアル(案)[EB/OL]. [2003-12-16].http: //www.mlit.go.jp/kowan/hazard_ map/5/shiryou2.pdf.

[13] 尹慶江, 王喜年, 吳少華. 鎮(zhèn)?？赡茏畲笈_風增水的計算[J]. 海洋學報, 1995, 17(6): 21-27.

[14] 吳少華, 王喜年, 于福江, 等. 連云港溫帶風暴潮及可能最大溫帶風暴潮的計算[J]. 海洋學報, 2002, 24(5): 8-18.

[15] 趙鑫. 錢塘江北岸可能最大熱帶氣旋引起的風暴高潮位研究[D]. 杭州: 浙江大學, 2009.

[16] Rego J L, Li C. On the importance of the forward speed of hurricanes in storm surge forecasting: A numerical study[J]. Geophsical Research Letters, 2009, 36: L07609, doi: 10.1029/2008GL036953.

[17] Peng M S, Jeng B F, Williams R T. A numerical study on tropical cyclone intensification. Part I: Beta effect and mean flow effect[J]. Journal of the Atmospheric Sciences, 1999, 56: 1404-1423.

[18] Niedoroda a W, Resio D T, Toro G R, et al. Analysis of the coastal Mississippi storm surge hazard[J]. Ocean Engineering, 2010, 37: 82-90.

[19] 中國氣象局. CMA-STI熱帶氣旋最佳路徑數(shù)據(jù)集(1949-2011)[DB/OL].[2011-07-29]. http: //www.typhoon.gov.cn/index.php?controller=spage&pid=169.

[20] 李心怡, 方偉華, 林偉.西北太平洋熱帶氣旋路徑及強度插值方法比較研究[J]. 北京師范大學學報(自然科學版), 2014, 50(2): 111-116.

[21] JTWC. JTWC tropical cyclone best track data(2001-2011)[DB/OL]. [2013-05-26]. http: //www.usno.navy. mil/NOOC/nmfc-ph/RSS/JTWC/best_tracks/wpindex. html.

[22] Tropical Cyclone Programme (WMO/TD-No. 196). Typhoon committee operational manual[R]. Geneva: World Meteorological Organization, 2011.

[23] 鄒燕, 趙平. 幾種臺風資料的對比及臺風年鑒數(shù)據(jù)的訂正[J]. 熱帶氣象學報, 2009, 25(3): 295-299.

[24] GB/T 19201-2006, 熱帶氣旋等級[S].

[25] Atkinson G D, Holliday C R. Tropical cyclone minimum sea level pressure maximum sustained wind relationship for the western north pacific[J]. Monthly Weather Review, 1977, 105: 421-427.

[26] Knaff J A, Harper B A. Tropical cyclone surface wind structure and wind-pressure relationships[R]. La Reunion : WMO International Workshop on Tropical Cyclones-VII, 2010.

[27] Vickery P, Skerlj P, Twisdale L. Simulation of hurricane risk in the U.S. using empirical track model[J]. Journal of Structure Engineering, 2000, 26(10): 1222-1237.

[28] Willoughby H, Rahn M. Parametric representation of the primary hurricane vortex. Part I: Observations and evaluation of the holland (1980) model[J]. Monthly Weather Review, 2004, 132: 3033-3048.

[29] Chan J C-L. An observational study of the physical processes responsible for tropical cyclone motion[J]. Journal of the Atmospheric Sciences, 1984, 41: 1036-1048.

[30] Liang X, Chan J C L. The effects of the full coriolis force on the structure and motion of a tropical cyclone. Part I: Effects due to vertical motion[J]. Journal of the Atmospheric Sciences, 2005, 62: 3825-3830.

[31] 董劍希, 仉天宇, 付祥, 等. 福建省沙埕港百年一遇臺風風暴潮的計算[J]. 海洋通報, 2008, 27(1): 9-16.

[32] Toro G R. Joint probability analysis of hurricane flood hazards for Mississippi[R]. Boulder: Risk Engineering, Inc, 2008.

[33] Chouinard L E, Liu C, Cooper C K. Model for severity of hurricanes in Gulf of Mexico[J]. Journal of Waterway, Port and Ocean Engineering, 1997, 123(3): 120-129.

[34] Vickery P J. Simple empirical models for estimating the increase in the central pressure of tropical cyclones after landfall along the coastline of the United States[J]. Journal of Applied Meteorology, 2005, 44: 1807-1826.

[35] Kaplan J, Demaria M. On the decay of tropical cyclone winds after landfall in the New England area[J]. Journal of Appled Meteorology, 2001, 40: 280-286.

[36] Opencyclone(臺風災(zāi)害風險模型)[DB/OL].[2003-12-16]. https: //github.com/OpenCyclone.

[37] 中國氣象局國家氣候中心. 全國氣候影響評價[M].北京: 氣象出版社, 2001-2010.

(本文編輯: 劉珊珊 李曉燕)

Parameterization of synthetic tropical cyclones at various scales for probable maximum storm surge risk modeling

LI Ying1,2,3, FANG Wei-hua1,2,3, LIN Wei1,2,3, YE Yan-ting1,2,3
(1. Key Laboratory of Environmental Change and Natural Disaster, Ministry of Education of China, Beijing Normal University, Beijing 100875, China; 2. State Key Laboratory of Earth Surface Processes and Resources Ecology, Beijing Normal University, Beijing 100875, China; 3. Academy of Disaster Reduction and Emergency Management, Ministry of Civil Affairs & Ministry of Education, the People’s Republic of China, Beijing Normal University, Beijing 100875, China)

Aug., 29, 2012

Tropical Cyclone; Probable Maximum Storm Surge; Parameterization; Risk Assessment

The parameterization of Probable Maximum Tropical Cyclones (PMTC) is of great importance to the computation of Probable Maximum Storm Surge (PMSS) for emergency evacuation planning at different scales of tropical cyclones. Based on the Best Track Dataset (1949-2011) of the western north pacific from CMA and other two databases, the radius of maximum wind from JTWC and NOAA, the quantitative relationships between the maximum winds and other parameters of various scales wereanalyzed and modeled. A method wasdeveloped to set maximum sustained wind, central pressure, radius of maximum winds, forward velocity and direction, and then to synthetize the full track with the above parameters. Taking Liangjiang County, Fujian province as an example area, 216 tropical cyclones were synthetized for tropical cyclones scaled at typhoon, severe typhoon and super typhoon, with decay and no-decay at 3 major directions. At last, parameter sensitivity, Holland B parameter setting, wind field parameter setting and computation resource requirement, astronomical tide, and dam failure, etc. were discussed.

P463.1; P47

A

1000-3096(2014)04-0071-010

10.11759/hykx20120829002

2012-08-29;

2013-04-11

公益性行業(yè)科研專項(201305020); 國家科技支撐項目(2012 BAK10B03)

李穎(1986-), 女, 山東淄博人, 博士生, 主要從事臺風災(zāi)害風險評估研究, E-mail: ying.li@mail.bnu.edu.cn; 方偉華, 通信作者,副教授, 主要從事自然災(zāi)害風險評估研究, E-mail: weihua.fang@ bnu.edu.cn