王培濤，于福江,2，趙聯(lián)大,2，侯京明

（1.國家海洋環(huán)境預(yù)報(bào)中心,北京 100081；2.國家海洋環(huán)境預(yù)報(bào)中心海洋災(zāi)害預(yù)報(bào)技術(shù)研究國家海洋局重點(diǎn)實(shí)驗(yàn)室，北京 100081）

1 引言

2004年12月26日蘇門答臘-安達(dá)曼Mw9.1地震引發(fā)的海嘯，不僅造成了印度洋周邊十余個國家的重大經(jīng)濟(jì)損失（超過100億美元），還奪去了近30萬人的生命，超過100萬人無家可歸[1-7]。2010年2月27日智利中部近海發(fā)生Mw8.8地震，并引發(fā)越洋海嘯，此次地震及海嘯在智利當(dāng)?shù)卦斐山?000人死亡，超過150萬所房屋建筑倒塌，至少200萬人受災(zāi)，遠(yuǎn)在大洋彼岸的中國在地震發(fā)生25個小時后監(jiān)測到海嘯，最大海嘯波為32 cm[8-9]。同年10月25日西蘇門答臘省明打威群島附近海域發(fā)生的Mw 7.7地震再次引發(fā)了海嘯，并造成近400人死亡[10]。

當(dāng)人類還未擺脫印度洋海嘯、智利海嘯的陰影時，又一次“惡夢”震驚世界—2011年3月11日14時46 min（當(dāng)?shù)貢r間）日本東北近海發(fā)生Mw9.0級地震，為日本歷史上最強(qiáng)地震，引發(fā)的越洋海嘯不僅對日本東北部近岸造成災(zāi)難性影響，還對環(huán)太平洋多個國家造成災(zāi)害[11-13]；海嘯近場的福島縣相馬海嘯波幅9.3 m、宮城縣石卷8.6 m以上、巖手縣宮古8.5 m以上,局部海嘯爬高達(dá)40 m[14];美國西海岸也監(jiān)測到了2.5 m的海嘯波。地震、海嘯、火災(zāi)和核泄漏構(gòu)成了災(zāi)害的主體，形成并發(fā)的復(fù)合型災(zāi)害，共造成15240人死亡死亡、8173人失蹤、112528棟房屋被毀、40多萬人無家可歸和100萬人缺少飲用水，綜合多家國際知名金融公司的預(yù)測，此次災(zāi)難造成的經(jīng)濟(jì)損失超過18萬億日元[15]。受此影響，我國發(fā)布了首個海嘯藍(lán)色警報(bào)，海嘯波震6 h后到達(dá)我國大陸沿海，我國沿海儀測最大海嘯波幅為55 cm,數(shù)值模擬研究發(fā)現(xiàn)此次地震海嘯在我國沿海局部地區(qū)海嘯波幅可達(dá)1.5 m[12]。

如果說2004年印度洋海嘯事件改變了人們對印度洋海嘯危險性的認(rèn)識[16]，那么，1960年、2010年智利海嘯、2011年日本海嘯則應(yīng)該改變?nèi)藗儗ξ覈粫茉窖蠛[影響的觀點(diǎn)[8-9,12,16]。加之，近年來頻發(fā)的地震海嘯災(zāi)害再次敦促各濱海國家進(jìn)一步完善海嘯預(yù)警系統(tǒng)，加強(qiáng)海嘯基礎(chǔ)理論研究、加大海嘯預(yù)警監(jiān)測投入，適時開展海嘯災(zāi)害風(fēng)險評估及區(qū)劃工作。早在上世紀(jì)90年代于福江等就利用自主開發(fā)的地震海嘯模型對沿海的多個核電站進(jìn)行可能最大海嘯增水計(jì)算；日本3.11地震海嘯后，國家海洋環(huán)境預(yù)報(bào)中心海嘯組針對區(qū)域海嘯評估了我國沿海所有在建和已建核電站海嘯危險性，但是上述研究成果均是針對區(qū)域海嘯和局地海嘯所開展的論述。本文主要關(guān)注越洋海嘯對研究區(qū)域危險性的影響評估。

2 我國東南沿海的地震海嘯風(fēng)險

我國是世界上最早記錄到地震海嘯的國家，2000多年來我國沿海共發(fā)生海嘯大約60多次[17]，其中大多數(shù)集中在我國臺灣和南海沿岸。我國東南沿海處于環(huán)太平洋地震帶上，既面臨著局地海嘯的威脅，也受區(qū)域和越洋海嘯的影響。影響我國的海嘯最有可能發(fā)生在南海東部馬尼拉海溝、臺灣島周邊海域、琉球群島。尤其是南海馬尼拉海溝，是國際公認(rèn)的海嘯潛在發(fā)生源地[18-19]。局地海嘯一般發(fā)生在近海和岸邊，海嘯波到達(dá)岸邊的時間很短，有時僅幾分鐘或幾十分鐘，通常猝不及防，往往造成嚴(yán)重危害。我國歷史上1867年臺灣基隆和1992年海南島的海嘯事件均屬于局地海嘯。越洋海嘯主要來自于太平洋海域，盡管有外圍島鏈的阻擋和大陸架地形的影響，越洋海嘯對我國東海沿岸和臺灣島的影響不容忽視。1960年、2010年智利海嘯和2011年日本海嘯事件均屬于越洋海嘯，我國東部沿海儀器監(jiān)測到55 cm的最大海嘯波幅。這并不能說明，在這幾次事件中我國東南沿海最大海嘯波為55 cm，在沒有儀器觀測的局部地區(qū)可能會有更大的海嘯波存在.我國沿海越洋海嘯風(fēng)險不容忽視[8-9,12]。

根據(jù)美國哈佛大學(xué)提供的震源機(jī)制解資料分析，環(huán)太平洋地震帶7級以上地震主要為逆沖型破裂，由于引發(fā)海嘯的地震震源機(jī)制類型是傾滑型或帶傾滑分量的走滑型的地震，從地震錯動方式上，該區(qū)的地震多數(shù)具備引發(fā)海嘯的條件。因此研究溫州甌江口淺灘地區(qū)海嘯風(fēng)險除應(yīng)考慮區(qū)域海嘯源（琉球海溝、馬尼拉海溝）的影響，也應(yīng)該考慮越洋海嘯對該地區(qū)的威脅。本文將針對該區(qū)域的越洋海嘯危險進(jìn)行計(jì)算評估。

3 越洋海嘯傳播特征及海嘯源選取

3.1 越洋海嘯傳播特征

海嘯通常是由海底地震、海底火山噴發(fā)、海底或海岸滑坡、崩塌、海底核爆炸或者小行星撞擊所產(chǎn)生的系列具有超長波長和周期（與風(fēng)浪相比）的大洋行波，其中海底地震是誘發(fā)海嘯的主要原因，在所有的海嘯事件中86%是由地震觸發(fā)，5%由火山活動引起，4%由海底滑坡產(chǎn)生，5%是由這些因素的綜合過程引起[20]，因此海嘯通常被稱為“地震海嘯”。按照海嘯影響范圍海嘯又可以分為越洋海嘯、區(qū)域海嘯和局地海嘯[21].越洋海嘯是指海嘯源距離受災(zāi)地區(qū)超過1000 km。越洋海嘯首波波長在量級相當(dāng)于震源長軸尺度，通常可達(dá)幾百千米，此時大洋水深和海嘯波長的關(guān)系為：<<1，因此海嘯波可視為淺水長波，即海嘯波在大洋中傳播時不發(fā)生頻散，傳播速度為，可達(dá)700 km/h，但波高卻只有1 m左右，甚至更小；加之海嘯波在大洋中受到的摩擦力較小，所以越洋海嘯可以在大洋中傳播數(shù)千海里而能量衰減甚少，對數(shù)千海里外的沿岸地區(qū)也會造成海嘯災(zāi)害。當(dāng)海嘯波傳播到近岸區(qū)域后由于水深急劇變淺，水體波動特性因受到地形影響而改變，主要表現(xiàn)為其傳播方向、海嘯高度以及剖面形狀都要隨之變化；特別當(dāng)水深沿波向線逐漸變淺使得原先的海嘯波波長急劇變短，波高驟增。如1960年、2010年智利海嘯、2011年日本海嘯，上述3次越洋海嘯奔襲24 h，海嘯波行進(jìn)了大約1.8×104km，海嘯波平均波速達(dá)到750 km/h（見圖1 a）,周期為60 min（見圖1b）。海嘯到達(dá)彼岸時海嘯波幅依舊高達(dá)數(shù)米。

3.2 海嘯源選取

我國東南沿海越洋海嘯威脅主要來自環(huán)太平洋地震帶，因此本文進(jìn)行評估計(jì)算的海嘯源取自歷次海嘯演習(xí)中UNESCO/IOC推薦的環(huán)太平洋地震帶上的潛在地震海嘯源。包括日本東北部潛在地震源、菲律賓東部潛在地震源、智利中南部潛在地震源、厄瓜多爾西部潛在地震源、中美洲潛在地震源、卡斯凱迪亞潛在地震源、阿留申群島潛在地震源、堪察加半島潛在地震源、湯加以及瓦努阿圖潛在地震源（見圖2，表1）。

圖2 UNESCO/IOC推薦環(huán)太平洋地震帶可能最大越洋海嘯震源位置分布圖

表1 環(huán)太平洋地震帶潛在震源參數(shù)

4 地震海嘯模式建立及驗(yàn)證

4.1 地震海嘯模式建立

地震海嘯模型是探索海嘯生成機(jī)制、評估海嘯災(zāi)害和建立實(shí)時海嘯預(yù)警系統(tǒng)的重要工具和技術(shù)手段。根據(jù)海嘯波的物理特性，海嘯數(shù)值模型通常包括：計(jì)算海嘯初始位移場的斷層模型(Genera?tion)；基于非線性淺水方程或類Boussinesq方程，計(jì)算海嘯波傳播的傳播模型(Propagation)；考慮沿岸地形、地物、采用動態(tài)移動邊界技術(shù)來實(shí)現(xiàn)海嘯爬高及與海岸相互作用的淹沒模型(Inundation)[2,8,22-23]。

4.1.1 海嘯源模型(斷層模型)

海嘯源模型是海嘯數(shù)值模型研究的基礎(chǔ)，它的適用性直接關(guān)系到海嘯波的傳播及海嘯與海岸的相互作用。斷層模型通過利用海床位移量來估算地震引起的初始水面高度，為海嘯數(shù)值模型計(jì)算提供初始條件。目前，Mansinha&Smylie(1971)[24]以及Okada(1985)[25]基于彈性錯移理論發(fā)展的兩套斷層模型被廣泛應(yīng)用，大量的研究和應(yīng)用實(shí)例表明此類模型對大部分海嘯源計(jì)算具有較好的適用性[8,26-27]。

4.1.2 海嘯傳播數(shù)值模型

海嘯在大洋和大陸架傳播時，水質(zhì)點(diǎn)的垂直加速度與重力加速度相比可看做一小量；這時水質(zhì)點(diǎn)的垂直運(yùn)動對壓力分布的影響可以忽略，所以不考慮頻散項(xiàng)的非線性淺水方程通常被用來作為海嘯的傳播模型。在近岸此類模型對海嘯波到達(dá)時間和海嘯波首波波高均能給出較精確的模擬結(jié)果，但對首波后的系列波形的預(yù)報(bào)卻有較大的偏差。Imamura F和Liu P.L認(rèn)同頻散作用的重要性，認(rèn)為模型物理頻散效應(yīng)可以通過選擇適當(dāng)?shù)臅r間和空間步長，用差分方程的數(shù)值頻散代替微分方程的物理頻散(physical dispersion)既保持了海嘯傳播過程的物理本質(zhì),又降低了求解物理頻散項(xiàng)的機(jī)時代價。

本文選用美國Cornell大學(xué)開發(fā)COMCOT(Cor?nell Multi-grid Coupled Tsunami Model)海嘯模型作為本文的研究模型。該模型采用基于多層網(wǎng)格嵌套的有限差分法。針對海嘯波不同物理特性，模式可以靈活配置所需坐標(biāo)系(直角/球面)和控制方程類型(線性/非線性)[28-30]。本次模擬采用三層嵌套網(wǎng)格（見圖3）。

本文中深水模塊采用球坐標(biāo)系下線性方程:

淺水模塊采用球坐標(biāo)系下非線性方程:

圖3 模型計(jì)算區(qū)域網(wǎng)格配置

圖4 2010年智利海嘯（a）、2011年日本海嘯（b）在近場及遠(yuǎn)場傳播的數(shù)值模擬與觀測數(shù)據(jù)對比

式中η為相對于平均海平面的自由表面位移；?為緯度；ψ為經(jīng)度；R為地球半徑；h為凈水深，H=h+η為總水深；P為沿經(jīng)度單位寬度的通量；Q為沿緯度單位寬度的通量；f為科氏力系數(shù)；g為重力加速度。Fx，F(xiàn)y分別為經(jīng)度和緯度方向的底摩擦力。

4.2 模型驗(yàn)證

為了驗(yàn)證模型的可靠性和適用性，本文利用該模型對2010年智利海嘯、2011年日本地震海嘯的近場、遠(yuǎn)場傳播進(jìn)行了數(shù)值模擬（模擬參數(shù)配置見文獻(xiàn)[8,12]），將模擬結(jié)果與觀測數(shù)據(jù)進(jìn)行對比。從對比的結(jié)果（見圖4）可以看出所采用模型對海嘯波到達(dá)時間及最大波幅的模擬計(jì)算具有較高的精度。但是模型對首波以后的高頻波動模擬仍存在一定的偏差，這與該模型所采用的海嘯源模型及基礎(chǔ)的水深數(shù)據(jù)精度有直接關(guān)系。但就海嘯評估工作，此模型精度滿足需求。

5 越洋海嘯對甌江口地區(qū)的影響評估計(jì)算及危險性分析

5.1 計(jì)算網(wǎng)格和地形數(shù)據(jù)

圖5 第三層嵌套網(wǎng)格范圍及評估點(diǎn)設(shè)置

越洋海嘯的計(jì)算分3層嵌套網(wǎng)格進(jìn)行（見表2），第一層太平洋范圍網(wǎng)格分辨率4'，第二層陸架范圍網(wǎng)格分辨率0.5'，第三層評估區(qū)所在的近岸區(qū)域（見圖5）網(wǎng)格分辨率0.1'。第一層的水深數(shù)據(jù)主要取自ETOPO1數(shù)據(jù)；第二層陸架范圍數(shù)據(jù)采用GEBCO（General Bathymetric Chart of the Oceans）數(shù)據(jù)與部分海圖數(shù)據(jù)融合；第三層數(shù)據(jù)由大比例尺海圖數(shù)字化獲取，其中近岸海陸交接面的分離采用strm 30衛(wèi)星數(shù)據(jù)訂正?？紤]海嘯波進(jìn)入大陸架，特別是行進(jìn)到研究區(qū)域時海嘯波的傳播特性，圍繞甌江口淺灘近岸區(qū)域選擇12個控制點(diǎn)作為該區(qū)域海嘯評估的主要參考依據(jù)。其中，1#—5#位置主要考慮海嘯波直接進(jìn)入該區(qū)域后對研究區(qū)域的影響，7#—10#位置主要考慮繞射波對評估區(qū)域的影響作用，11#、12#位置主要考查開闊海域海嘯波特征參數(shù)及模式評估驗(yàn)證，6#、7#綜合考慮波浪反射、繞射、散射等共同作用后海嘯波幅的高度。

5.2 計(jì)算結(jié)果及危險性分析

從各越洋海嘯源計(jì)算結(jié)果（見表3）可以看出：智利海嘯源產(chǎn)生的海嘯在評估對象近岸的海嘯波最大，其中06點(diǎn)處的最大海嘯波幅為100 cm；其次是菲律賓東部、日本海嘯源在06點(diǎn)海嘯波幅分別達(dá)到97 cm、94 cm，厄瓜多爾、堪察加海嘯源、瓦努阿圖、湯加、中美洲、卡斯凱迪亞海嘯源同樣最大海嘯波幅出現(xiàn)在06點(diǎn)位，分別為34 cm、64 cm、60 cm、29 cm、85 cm和62 cm。06點(diǎn)位處的海嘯波幅最大及淺灘南側(cè)海嘯波幅偏大分析可能的原因：越洋海嘯傳播至琉球海溝后，經(jīng)過繞射及折射調(diào)整后，波傳播方向與06點(diǎn)、05點(diǎn)位連線的方向接近一致，06點(diǎn)位又處于最外側(cè)；而淺灘北側(cè)由于受到群島環(huán)繞，海嘯波繞射后能量受到損失。按照太平洋海嘯警報(bào)中心海嘯危險性等級分類標(biāo)準(zhǔn)（見表4）劃分：智利海嘯在所評估的12個點(diǎn)位處淺灘附近的10個點(diǎn)的海嘯波幅都達(dá)到Ⅱ級；菲律賓東部、日本東北、阿留申群島、中美洲以及卡斯凱迪亞海嘯源只在01#—08#點(diǎn)達(dá)到Ⅱ級，厄瓜多爾、堪察加、瓦努阿圖海嘯源在05點(diǎn)、06點(diǎn)處的海嘯波達(dá)Ⅱ級，湯加海嘯源各點(diǎn)處的海嘯波幅均不會對評估區(qū)近海、近岸造成影響；同時可知：各海嘯源對評估區(qū)近岸均不構(gòu)成淹沒風(fēng)險，但是可以看出菲律賓東部、南美智利、日本東北部海嘯以及北美洲中部潛在海嘯源產(chǎn)生的海嘯對評估區(qū)近海工程及作業(yè)船會有一定的影響。

表2 網(wǎng)格及模型參數(shù)設(shè)置

表3 太平洋潛在地震海嘯在評估點(diǎn)處海嘯波幅計(jì)算結(jié)果（單位/cm）

表4 海嘯危險性等級劃分

6 結(jié)語與討論

（1）從太平洋潛在地震海嘯過程數(shù)值計(jì)算結(jié)果可以看到，評估區(qū)受越洋海嘯過程的影響比較小，普遍小于100 cm。但由于工況位置不同，評估區(qū)所受海嘯影響也各不相同，主要表現(xiàn)為：評估區(qū)南側(cè)比北側(cè)海嘯危險性大、東側(cè)比西側(cè)的危險性大。這主要與海嘯波傳播方向及局部的地形變化有關(guān)；

（2）從10個越洋海嘯評估計(jì)算結(jié)果分析可知，南美中南部海嘯源、日本東北部海嘯源、菲律賓東部、北美洲中部海嘯源比較容易對我國東南沿海造成影響。但厄瓜多爾、湯加、瓦努阿圖、堪察加等海嘯源由于其海嘯能量傳播主軸方向不直接面向我國東南沿海方向，因而上述幾處海嘯源地發(fā)生的海嘯不易對評估區(qū)造成影響；

（3）值得注意的是如果日本3.11海嘯源的位置發(fā)生在日本“南海海槽”，那么對評估區(qū)將造成較大影響。另外，沖繩海槽地震潛在海嘯對評估區(qū)域有較大威脅，在充分認(rèn)識越洋海嘯對該研究區(qū)域的危險性后，加強(qiáng)該區(qū)域可能最大海嘯地震分析研究，評估計(jì)算研究區(qū)可能最大的海嘯將顯得尤為重要。

[1]陳運(yùn)泰,楊智嫻,許力生.海嘯、地震海嘯與海嘯地震[J].物理,2005,34(12):864-872.

[2]姚遠(yuǎn),蔡樹群,王盛安.海嘯波數(shù)值模擬的研究現(xiàn)狀[J].海洋科學(xué)進(jìn)展,2007,25(4):487-494.

[3]王培濤,趙聯(lián)大,于福江,等.海嘯災(zāi)害數(shù)值預(yù)報(bào)技術(shù)研究現(xiàn)狀[J].海洋預(yù)報(bào),2011,28(3):74-79.

[4]陳虹,李成日.印尼8.7級地震海嘯災(zāi)害及應(yīng)急救援[J].國際地震動態(tài),2005,316(4):22-26.

[5]Kenji S,Than TA.Report On Post Tsunami Survey Along The Myanmar Coast For The December 2004 Sumatra-Andam An Earthquake [R].Active Faults-Paleoseismic Study,2005(5):161-188.

[6]任葉飛.基于數(shù)值模擬的我國地震海嘯危險性分析研究[D].哈爾濱：中國地震局工程力學(xué)研究所,2007.

[7]褚芹芹.印尼海嘯過程數(shù)值模擬研究及其風(fēng)險評估[D].青島:中國海洋大學(xué),2011.

[8]于福江,王培濤,趙聯(lián)大,等.2010年智利地震海嘯數(shù)值模擬及其對我國沿海的影響分析[J].地球物理學(xué)報(bào),2011,54(4):918-925.

[9]于福江,原野,趙聯(lián)大,等.2010年2月27日智利8.8級地震海嘯對我國的影響分析[J].科學(xué)通報(bào),2011,56(3):239-246.

[10]AndrewVN,GavinH,WeiY,etal.The25October 2010 Mentawai tsunami earthquake from real-time discriminants,finite-fault rupture and tsunami excitation[J].Geophysical Research Letters,2011,38,L05302,doi:10.1029/2010GL046198.

[11]溫瑞智,任葉飛,李小軍.日本Mw9.0級地震海嘯數(shù)值模擬與啟示[J].國際地震動態(tài),2011,388(4):22-27.

[12]王培濤,于福江,趙聯(lián)大,等.2011年3月11日日本地震海嘯越洋傳播及對中國影響的數(shù)值分析[J].地球物理學(xué)報(bào),2012,55(9):3088-3096.

[13]吳佳曄.日本東北大地震對我們的警示[J].四川理工學(xué)院學(xué)報(bào)(自然科學(xué)版),2011,24(2):125-128.

[14]NobuhitoM,TomoyukiT,TomohiroY,etal.Survey of 2011Tohoku earthqua ketsunamiinun dationandrun-up[J].Geophysical Research Letters,Vol.38,L00G14,6PP,2011.doi:10.1029/2011GL049210.

[15]吳新燕,苗崇剛,顧建華,等.日本東海岸9.0級地震應(yīng)對給我國地震應(yīng)急的啟示[J].國際地震動態(tài),2011,388(4):38-43.

[16]王培濤,于福江,趙聯(lián)大,等.越洋海嘯的數(shù)值模擬及其對我國的影響分析[J].海洋學(xué)報(bào),2012,34(2):39-47.

[17]趙聯(lián)大,于福江,王培濤.我國的海嘯風(fēng)險與預(yù)警.國家綜合防災(zāi)減災(zāi)與可持續(xù)發(fā)展論壇文集(上)[C].國家減災(zāi)委專家委員會.中國·北京：2010年5月:363-367.

[18]Liu L F,Wang X M,Salisbury A J.Tsunami hazard and early warning system in South China Sea[J].Journal of Asian Earth Sciences,2009,36(1):2-12.

[19]李林燕,毛獻(xiàn)忠.南海海嘯反問題預(yù)報(bào)模式[J].水動力學(xué)研究與進(jìn)展（A輯）,2012,27(1):62-67.

[20]Rhodes B,Tuttle M,Horton B,et al.Paleotsunami research[J].EOS TransAGU,2006,87(21):205-209.

[21]Intergovernmental Oceanographic Commission.Tsunami Glossary[K].Paris,UNESCO.IOC Technical Series,2008.

[22]祝會兵.海嘯模型研究進(jìn)展(I)[J].寧波大學(xué)學(xué)報(bào)(理工版),2006,19(4):536-541.

[23]張超凡,石耀霖.海嘯災(zāi)害的數(shù)值模擬研究[J].中國科學(xué)院研究生院學(xué)報(bào),2008,25(3):289-296.

[24]Mansinha L,Smylie D E.The displacement fields of inclined faults[J].Bulletin of the Seismological Society of America,1971,61(5):1433-1440.

[25]Okada,Y,Surface deformation due to shear and tensile faults in a half-space[J].Bulletin of the Seismological Society of America,1985,75(4),113-1154.

[26]TitovVV,Gonzalez FI.Implementation and testing of the method of Splitting Tsunami(MOST)model[R].NOAA Technical Memorandum ERL,1997.

[27]TitovVV,Numericalmodeling of longwaverun-up[D].LosAngeles:University of Southern California,1997.

[28]WangXM.USERMANUALFORCOMCOTVERSION1.7[K].Cornell University,2009.

[29]Wang X M,Liu LF.An analysis of 2004 sumatra earthquake fault plane mechanisms and indian ocean tsunami[J].J Hydraulic Res,2006,44(2):147-154.

[30]Wang X M,Liu L F.Numerical simulations of the 2004 indian ocean tsunamis coastal effects[J].Journal of Earthquake and Tsunami,2007,1(3):273-297.