溫燕林, 于海英, 朱艾斕, 宋治平, 李春峰

(1.同濟大學(xué)海洋與地球科學(xué)學(xué)院,上海 200092; 2.上海市地震局,上海 201203;3.中國地震臺網(wǎng)中心,北京 100045)

?

設(shè)定琉球海溝發(fā)生罕遇地震評估我國東南沿海地區(qū)的海嘯風(fēng)險①

溫燕林1,2, 于海英2, 朱艾斕2, 宋治平3, 李春峰1

(1.同濟大學(xué)海洋與地球科學(xué)學(xué)院,上海 200092; 2.上海市地震局,上海 201203;3.中國地震臺網(wǎng)中心,北京 100045)

摘要:根據(jù)構(gòu)造相似條件分析,琉球海溝與日本海溝、智利海溝、印尼巽他海溝一樣具備發(fā)生9級罕遇超巨大地震的可能。在對近幾年來全球發(fā)生的超巨大地震參數(shù)及構(gòu)造對比分析的基礎(chǔ)上,設(shè)定琉球海溝9.0級地震參數(shù),并將其引發(fā)的海嘯進行數(shù)值模擬研究。結(jié)果表明,該地震可引發(fā)初始波高為8 m的海嘯,臺灣東北部半小時后遭受10 m以上海嘯,3～4小時左右傳至浙南、閩北沿岸,近岸各處波高在1～2 m；5小時左右傳至浙北、粵北沿岸,浙江近岸各處波高在2 m左右,廣東沿海、臺灣海峽由于臺灣島的正面阻擋,海嘯波高低于50 cm；8小時后靠近上海海岸線,最大波高約1 m。海嘯的上岸高度與海岸附近的海深和海岸線的形態(tài)密切相關(guān),我國東南海域地形變化復(fù)雜、海灣眾多,對海嘯波有放大作用,模擬結(jié)果可能比實際海嘯偏小。我國沿海地區(qū)分布著不少已建和在建的核電廠,在核電設(shè)計時未考慮海嘯,一旦發(fā)生這種罕遇地震海嘯則影響不可忽視,尤其是若與風(fēng)暴潮、天文大潮疊加則可能出現(xiàn)嚴重后果。由于核電安全要求萬無一失,故須制訂有效預(yù)警和應(yīng)對措施。

關(guān)鍵詞:琉球海溝; 罕遇地震; 海嘯模擬; 海嘯風(fēng)險圖; 核電廠安全

0引言

20世紀(jì)以來觀測到的全球9級以上超巨大地震,即1952年俄羅斯堪察加MW9.0、1960年智利MW9.5、1964年美國阿拉斯加MW9.2及2004年印尼蘇門答臘—安達曼MW9.1地震都發(fā)生在大洋邊緣海溝,引發(fā)了巨大海嘯,造成重大人員傷亡和財產(chǎn)損失。特別是2011年3月11日日本東北方太平洋海域發(fā)生的罕遇MW9.0地震還造成震中附近的福島核電廠遭受超設(shè)計水平的海嘯襲擊,導(dǎo)致嚴重核泄漏等次生災(zāi)害,引發(fā)全球?qū)穗姀S安全問題的關(guān)注。這次地震是發(fā)生在太平洋邊緣日本海溝的特大俯沖型地震,是有記錄以來日本最大的地震。從板塊構(gòu)造理論來看,海溝是兩大板塊匯聚碰撞之處,能量積聚很快,最易發(fā)生俯沖型超巨大地震,而俯沖帶發(fā)生的大地震又最易引發(fā)大海嘯。西太平洋邊緣的琉球海溝俯沖帶地震活動同樣強烈,有資料記載的巨大地震有1911年琉球群島北部的奄美大島外海8.2級地震和1920年琉球海溝南端臺灣花蓮?fù)夂?.0級地震。我國華東沿海地區(qū)正對琉球海溝,若是該處發(fā)生罕見的9級地震并引發(fā)海嘯,將對我國造成多大沖擊？上述幾次太平洋越洋大海嘯傳播過來時受摩擦衰減和琉球島鏈阻擋,對我國沿海基本沒有造成災(zāi)害影響,而琉球海溝發(fā)生罕遇大地震海嘯的情況雖未發(fā)現(xiàn)過,但這種可能性不能不考慮。我國東南沿海各省都有已建和在建核電廠,鑒于日本核電廠的教訓(xùn),本文擬設(shè)定琉球海溝發(fā)生9級罕遇地震并引發(fā)海嘯,進行數(shù)值模擬,定量評估海嘯對我國東南沿海地區(qū)的影響,以期為核電廠海嘯風(fēng)險評估提供參考。

1構(gòu)造環(huán)境與地震海嘯

琉球海溝是歐亞板塊與菲律賓海板塊匯聚邊界,總體走向NE-SW向,向東南方向凸出。琉球海溝地區(qū)在地質(zhì)構(gòu)造上為典型的太平洋溝-弧-盆體系,即琉球海溝-琉球島弧-沖繩海槽弧后盆地構(gòu)造體系[1]。琉球島弧是雙列島弧,內(nèi)弧是水下火山或島嶼,外弧是琉球群島。內(nèi)弧從日本九州向西南延伸,經(jīng)吐噶喇群島至沖繩島,呈楔狀插入沖繩海槽與外弧之間,是從第四紀(jì)初發(fā)展至今的一條活火山帶。琉球群島為新生的火山島,是東海陸架海與西北太平洋的天然分界。群島由100多個島嶼組成,北起日本九州南端,南至我國臺灣島附近,是一個向東南突出的弧形島鏈,延伸 1 200 km,自北而南分別由大隅群島和吐噶喇列島、奄美群島和沖繩群島、宮古群島和八重山群島組成。琉球群島東側(cè)是琉球海溝,水深超過6 km。

琉球地震火山帶屬于環(huán)太平洋地震帶的一部分,是至今仍在活動的強地震帶。該處地震活動頻度高、震級大,僅發(fā)生在20世紀(jì)7級以上的大地震就有四十余次,特別引人注目的有1911年奄美大島外海的8.2級大地震、1920年臺灣花蓮?fù)夂５?.0級大地震和1938年宮古島西北的7.7級地震(圖1)。據(jù)統(tǒng)計,琉球島弧處90%的中源地震集中在沖繩海槽,琉球群島處多為淺震[2]。震源機制解表明該地區(qū)主壓應(yīng)力方向通常與琉球島弧正交,反映菲律賓海板塊向西的俯沖與擠壓。

菲律賓海板塊在日本東海—南海的俯沖帶也常常發(fā)生巨大地震。據(jù)統(tǒng)計,日本駿河灣—南海海槽每隔100～150年左右就會發(fā)生8級以上巨大地震,如日本1498年和1707年南海8.6級、1854年南海和東南海接連8.4級、1944東南海8.0級與1946南海8.1級地震(圖1),每次都引發(fā)了大海嘯[3-4],而且影響到我國沿海。1498年9月20日日本廣大地區(qū)發(fā)生地震(震中推測在日本南海,震級8.6),京都、三河、熊野震感最為強烈,同日我國江浙多處發(fā)生水溢。1707年10月28日日本南海發(fā)生8.6級地震(震后富士山大規(guī)模噴發(fā)),浙江吳興縣雙林地震水涌,烏青鎮(zhèn)河水暴漲,海鹽縣地震水沸；地震不僅波及浙江,錢塘江口亦有海嘯反映,但強度甚低。1854年12月24日日本南海8.4級地震引發(fā)東京以東的房總半島至九州太平洋沿岸海嘯,最高波高16.1 m。這次日本地震海嘯對我國的影響最為明顯[5],如江蘇如皋:海嘯,淹斃多人；上海:黃浦水沸,有高二三尺者；海寧:河水無風(fēng)自涌,如潮漲落,池沼皆然;慈溪:河水驟騰三四尺；臨海:海潮泛濫,城鄉(xiāng)溝池積潦,同時俱沸,歷二時止,沿海廬舍多被淹沒?；『笈璧匾矔a(chǎn)生地震海嘯,如日本西海岸一側(cè)的日本海發(fā)生過一些局地海嘯,但規(guī)模均小于其東海岸太平洋一側(cè)的越洋海嘯。琉球島弧后的沖繩海槽盆地與日本海類似,多形成局地海嘯。1771年4月24日石垣島西北沖繩海槽發(fā)生了7.4級地震,引發(fā)局地大海嘯,由于地形效應(yīng),波高在石垣島和宮古島海岸上沖超過10 m[6];1938年沖繩海槽中部(25.5°N,125°E)發(fā)生7.7級地震,也引發(fā)了海嘯;但均未見有對我國沿海造成影響的記載。

圖1　 琉球海溝—日本南海海槽俯沖帶歷史上8級以上地震破裂范圍示意圖Fig.1　Rupture ranges of historical M≥8.0 earthquakes on Ryukyu Trench-Nankai Trough subduction zone

琉球海溝與日本南海海槽相接,同屬于菲律賓海板塊西北面與歐亞板塊交界的俯沖帶,二者間被九州—帕勞海嶺脊分割,向東北與日本海溝、伊豆—小笠原海溝交接,構(gòu)成三聯(lián)點,向南西與臺灣島交接(圖1),與已發(fā)生過超巨大地震的印尼巽他海溝、日本海溝的俯沖構(gòu)造相似。因此琉球海溝同樣具備8級以上乃至9級地震發(fā)生的能力,而海溝處發(fā)生的板間特大逆沖地震常常又會引發(fā)巨大海嘯。已有的實例表明,日本東北地區(qū)的大地震海嘯不會影響到我國,其東海、南海地區(qū)的大海嘯則對我國沿海有一定影響。由此分析,若琉球海溝發(fā)生巨大地震海嘯,會對我國東南沿海有更明顯的影響。

2震源及海嘯模型

據(jù)日本氣象廳等進行的震源破裂時空過程反演,2011年日本東北地區(qū)MW9.0地震破裂長度約500 km,寬度約200 km,最大滑動量約18 m,觀測到的震中附近海底隆起5～8 m,震中北邊100多 km外的兩個深海驗潮站觀測到的海嘯最大波高達5 m。日本海嘯波源區(qū)550 km×200 km,最大沖進陸地10 km,海嘯強度為IV級(http://outreach.eri.u-tokyo.ac.jp/eqvolc/201103_tohoku/)。2004年印尼蘇門答臘—安達曼海溝MW9.1地震破裂長度約1 200 km,寬度100 km,最大滑動量約25 m,初始波高超8 m[7-8]。與同具有溝-弧-盆海底構(gòu)造體系的安達曼海溝相比,日本海溝后的本州島則是比較寬的大島(島弧寬度較大,溝弧間距離較寬),故俯沖破裂寬度較大,破裂長度較短。每次俯沖地震都會使溝弧間距離縮短,此次地震后日本島向海溝方向移動了24 m就是例證。而琉球島弧狹窄(島弧寬度窄,溝弧間距離較短),說明琉球島弧地殼滑入海溝多,只剩下一個個分割的小島,是溝-弧-盆體系演化的中期形態(tài),日本列島則是早期形態(tài),晚期形態(tài)是智利海溝后那樣火山弧與大陸板塊拼合到一起。因此琉球、安達曼俯沖帶類型破裂寬度較窄?；谝陨险J識,本文設(shè)定琉球海溝MW9.0地震震源區(qū)參數(shù)為:長度1 000 km,寬度100 km,平均滑動量18 m。

地質(zhì)和地球物理學(xué)家們對琉球地區(qū)的構(gòu)造特征進行過深入研究,并取得了不少成果[9-12]。臧紹先等討論了琉球海溝不同地段貝尼奧夫帶的形態(tài)和應(yīng)力狀態(tài),認為位于29°N附近的吐噶喇海峽是應(yīng)力場特征的一個分界線。進一步研究發(fā)現(xiàn),吐噶喇海峽北側(cè)俯沖帶傾角陡,屬于拉張型；以南傾角較緩,屬于擠壓型。再考慮到20世紀(jì)琉球海溝發(fā)生的2次8級以上地震位置,下一次大地震最有可能發(fā)生在未破裂的地震空區(qū),因此本文假設(shè)琉球海溝9級地震破裂段位于八重山群島到奄美群島之間的俯沖帶。根據(jù)琉球地區(qū)地震層析成像和貝尼奧夫帶分布反映的俯沖板片構(gòu)造特征,琉球海溝至島弧方向100 km寬度內(nèi)俯沖角度十分平緩(小于30°),震源深度都小于30 km,表明菲律賓海板塊在琉球海溝處低傾角俯沖至琉球群島地殼下面[13]。震源機制解表明,印尼與日本兩次地震的斷層面傾角均小于30°、滑動角接近90°,都是低角度逆沖運動[14](http://www.geol.tsukuba.ac.jp/%7Eyagi-y/EQ/Tohoku/)。因此本文設(shè)定琉球海溝9級地震為傾角30°、滑動角90°的逆沖運動(按該構(gòu)造約束下可產(chǎn)生的最大位錯設(shè)置),震源深度為20 km。

海底板塊交界處斷層的錯動是引發(fā)地震海嘯的主要原因[15]。地震引發(fā)大海嘯的條件:震級要大,水要深,水體被抬升幅度(即海底變形)要大。由斷層錯動產(chǎn)生的海底變形可以應(yīng)用彈性位錯理論計算[16]。此外地震破裂速度(破裂速度小于1 km/s難以產(chǎn)生大海嘯)、震源區(qū)上覆的沉積層密度也影響生成海嘯的大小[17-18]。本文暫不考慮沉積層因素,并且設(shè)定琉球9級地震破裂速度滿足形成大海嘯條件,因此可以假設(shè)海底變形是瞬態(tài)過程,并且流-固耦合較好,即海面初始位移與海底變形是一致的(圖2)。

海嘯是一種淺水長波,其模型采用球坐標(biāo)系下的非線性淺水波方程,同時考慮了地球自轉(zhuǎn)時的科氏力和海底摩擦作用,表達式如下[19]:

圖2　地震海嘯模型Fig.2　 The earthquake-induced tsunami model

?η?t+1Rcosφ?P?Ψ+??φ(Qcosφ)é?êêù?úú=-?h?t=0?P?t+ 1Rcosφ??ΨP2H?è???÷+1R??φPQH?è???÷+ gHRcosφ?η?Ψ-fQ+Fx=0?Q?t+1Rcosφ??ΨPQH?è???÷+1R??φQ2H?è???÷+ gHR?η?φ+fP+Fy=0ì?í??????????????

(1)

式中:η為相對于平均海平面的自由表面位移;h為海水靜深度;H=h+η為總水深;R為地球半徑;f為科氏力系數(shù);P、Q分別為水流沿經(jīng)度Ψ(東西方向)、緯度φ(南北方向)單位長度的流量;Fx、Fy分別表示經(jīng)度、緯度方向的海底摩擦作用。采用Manning底部摩擦公式:

(2)

式中:n是Manning粗糙系數(shù)。海底摩擦影響著淺水區(qū)海嘯傳播過程的水動力特性和爬高過程。

3數(shù)值模擬及分析

海嘯數(shù)值模擬采用了康奈爾大學(xué)P.L.F.Liu等開發(fā)的COMCOT 模式,該模式已被證明是一個成熟、有效的海嘯計算模式[19]。本文模擬區(qū)域為113°～135°E,20°～35°N,空間網(wǎng)格尺度采用規(guī)則的1′弧度網(wǎng)格(約1 852 m),網(wǎng)格數(shù)為1 321×901,時間步長1 s。計算區(qū)域水深及陸地地形數(shù)據(jù)來源于美國國家地球物理數(shù)據(jù)中心(NGDC)的ETOPO1數(shù)據(jù)庫。該模式在陸地邊界采用反射邊界條件計算,在水區(qū)邊界則采用輻射開邊界條件以“透浪”消波。

考慮到琉球海溝俯沖帶呈弧形,計算時將1 000 km 破裂帶分成相接的3段同時破裂,將這3個子事件組合成一次9級事件[圖3(a)],然后模擬9級地震海嘯在我國海域及鄰區(qū)傳播10小時的過程,海嘯波場見圖3。9級地震產(chǎn)生了初始波高超過8 m的海嘯。

圖3　海嘯傳播示意圖Fig.3　Snapshots of simulated tsunami propagation

我們共設(shè)置12個虛擬海嘯監(jiān)測站,其中大陸沿海9個,臺灣東部3個[編號及分布見圖3(a)中“△”標(biāo)示]。虛擬監(jiān)測站記錄了本次模擬9級地震海嘯的時程曲線(圖4),這些曲線記錄了海嘯到達各地的時間和波高變化信息。

由海嘯傳播過程可以看出,9級地震可產(chǎn)生初始波高8 m多的海嘯,向我國東海和太平洋兩個方向傳播,日本和韓國亦將受到海嘯沖擊。琉球諸島幾分鐘后首先遭到超過10 m的巨大海嘯掃蕩,最大溯上到20 m左右；十幾分鐘后海嘯到達臺灣東部,臺灣東北海岸受到5～10 m的海嘯襲擊,半小時后海嘯繞過臺灣最南端進入南海,海嘯約一小時后進入東海大陸架傳播。當(dāng)由深水傳到大陸架時,由于水深很快變淺,速度受阻,波高會驟然升高,但隨后繼續(xù)傳播,由于水深越來越淺,速度也越來越慢,波高就會受到越來越大的摩擦而逐漸衰減。3 h后海嘯最先接近我國大陸的浙閩交界和閩粵交界地帶海岸(廣東汕頭以南海域和福州—溫州海域),4 h左右從臺灣島南北兩端繞過來的海嘯在臺灣海嘯相遇,臺灣西部海岸和福建大部份海岸也將受到海嘯影響。廈門地區(qū)要比南邊的汕頭和北邊的福州晚半個多小時遭到海嘯,不過由于臺灣島的阻擋,海嘯不是正面?zhèn)鞑ミ^來的,其威力比廈門南北兩邊地區(qū)的小,到達廈門海岸的海嘯波高約0.5 m。3.5～7 h內(nèi)海嘯將陸續(xù)抵達浙江和廣東沿海各地,浙江沿海普遍波高2 m左右,最大可達3 m；而傳到廣東沿海的海嘯由于先繞過了臺灣南端(海嘯傳到臺灣南端波高僅為1 m左右),再加上南海大陸架衰減,抵岸波高都在0.5 m以下。海嘯傳至上海海岸線約在8 h后,最大波高接近1 m(圖4)。從圖4的海嘯時程曲線可以看出,虛擬臺從1～12臺記錄的海嘯波周期總體上是由長變短,這跟各臺水深條件有關(guān)。1～12臺分布由北向南,海岸地形由平緩到復(fù)雜,水深變化也是南部臺比北部臺劇烈,造成海嘯波形在近岸傳播時南部臺比北部臺變化更劇烈,頻散現(xiàn)象嚴重,因此南部臺記錄的高頻成分較多,北部臺記錄的海嘯頻率較單一。由于海嘯的上岸高度與海岸附近的海深和海岸線的形態(tài)、陸地的地形以及海灣的朝向密切相關(guān),一般海灣內(nèi)部的高度大于海灣開口處,因此海嘯進入杭州灣和長江口后波高會有所增大,考慮到網(wǎng)格的分辨率和計算誤差,預(yù)計最大波高仍小于2 m。也就是說,罕遇情況下上海地區(qū)遭受的最大海嘯為Ⅰ級小型海嘯。根據(jù)以上模擬結(jié)果繪制了我國大陸沿海及臺灣部分地區(qū)可能遭遇的最大海嘯風(fēng)險圖(圖5)。

圖4　12個虛擬監(jiān)測站記錄的海嘯時程曲線Fig.4　Time history curves of simulated tsunami in twelve virtual tidal stations

4結(jié)論與討論

海嘯是從海底到海面的海水整體擾動以波的形式傳播出去,是一種淺水波。海嘯傳播速度只與海水深度有關(guān),傳播到大陸架后隨水深變淺,速度變緩。海水只要從海底到海面受到整體擾動都能形成海嘯,只是深淺不同能量不同,發(fā)生在深海激起的海嘯能量大,發(fā)生在淺海其海嘯能量小,且波長比深海短,衰減快,到岸破壞力小,而我國東海淺水大陸架寬廣,故通常認為我國東部沿海不會遭受海嘯災(zāi)害。但我們的研究表明,琉球海溝發(fā)生罕遇9級地震條件下,波源區(qū)可達上千米,琉球群島都將遭受滅頂之災(zāi),臺灣臺北部分地區(qū)可能遭受最大10 m左右的海嘯,閩北及浙江沿海即使有東海淺水大陸架摩擦衰減,大部分地區(qū)可能遭受的最大海嘯仍超過2 m,上海及鄰區(qū)海岸最大海嘯在1 m左右,而海嘯傳到臺灣海峽及廣東沿海時因臺灣島阻擋作用會變得很小。這說明琉球罕遇海嘯對我國沿海一些地方的影響也是不可忽視的。一旦琉球海溝發(fā)生海嘯,有關(guān)部門可參考本文的研究結(jié)果對海嘯影響程度迅速估算,為制定及時、合理的防災(zāi)減災(zāi)措施提供依據(jù)。

圖5　我國大陸沿海和臺灣地區(qū)的海嘯風(fēng)險圖Fig.5　The tsunami risk map in China mainland coast and Taiwan region

我國核電廠設(shè)計中,認為我國近海不易產(chǎn)生大海嘯,遠洋海嘯傳過來也很小,主要考慮的是最大風(fēng)暴潮與最大天文大潮組合條件的影響[20-21]。但是對于核電廠這種安全要求“萬無一失”的工程,在選址與設(shè)計時還是應(yīng)該考慮罕遇海嘯和各種海浪(風(fēng)暴潮,天文大潮)組合條件下的安全問題。我國在沿海核設(shè)施都興建有防波堤,是按最大風(fēng)暴潮和最高天文潮組合條件設(shè)計的,從上面的模擬結(jié)果看這些防波堤用以防范單獨的海嘯應(yīng)該沒問題,但若遇罕遇海嘯疊加上最大風(fēng)暴潮和天文大潮組合則面臨超設(shè)計水平的危險。即使在設(shè)計中認為這種組合概率很小而沒有考慮,那么在運營中也應(yīng)對此可能的風(fēng)險情況制訂有效的防范應(yīng)對措施予以補救。此外還需注意的是,這樣的大地震極易造成沖繩海槽兩側(cè)海底斜坡出現(xiàn)滑坡,也有誘發(fā)琉球火山爆發(fā)的可能,這兩種情況都會引發(fā)海嘯,與地震海嘯疊加則會大大增加海嘯能量,使海嘯成倍放大,核電廠設(shè)計的極限地震海嘯安全應(yīng)該考慮這種極端情況。

參考文獻(References)

[1]吳時國,喻普之.海底構(gòu)造學(xué)導(dǎo)論[M].北京:科學(xué)出版社,2006:247-268.

WU Shi-guo,YU Pu-zhi. Introduction to Marine Tectonics[M].Beijing:Science Press,2006:247-268.(in Chinese)

[2]秦蘊珊.東海地質(zhì)[M].北京:科學(xué)出版社,1987.

QIN Yun-shan.Geology of the East China Sea[M].Beijing:Science Press,1987.(in Chinese)

[3]全球重大災(zāi)害性地震目錄(2150 B.C.-1991 A.D.)[M].北京:地震出版社,1996.

The Catalogue of Global Fetal Calamity Earthquakes (2150 B.C.-1991A.D.)[M].Beijing:Seismology Press,1996.(in Chinese)

[4][日]金井清,著.工程地震學(xué)[M].常寶琦,張虎男,譯.北京:地震出版社,1987:46-47.

KIYOSHI Kanai.Engineering Seismology[M].Translator:CHANG Bao-qi,ZHANG Hu-nan.Beijing:Seismology Press,1987:46-47.(in Chinese)

[5]王鋒,劉昌森,章振銓.中國古籍中的地震海嘯記錄[J].中國地震,2005,21(33):437-443.

WANG Feng,LIU Chang-sen,ZHANG Zhen-quan.Earthquake Tsunami Record in Chinese Ancient Books[J].Earthquake Research in China,2005,21(3):437-443.(in Chinese)

[6]Nakamura M.Source Fault Model of the 1771 Yaeyama Tsunami,Southern Ryukyu Lslands,Japan,Inferred from Numerical Simulation[J].Pure appl Geophys,2006,163:41-54.

[7]Chlieh M,Avouac J P,Hjorleifsdottir V,et al.Coseismic Slip and Afterslip of the GreatMW9.15 Sumatra—Andaman Earthquake of 2004[J].Bull Seism Soc Amer,2007,97(1):152-173.

[8]Geist E L,Titov V V,Arcas D,et al.Implications of the 26 December 2004 Sumatra—Andaman Earthquake on Tsunami Forecast and Assessment Models for Great Subduction-Zone Earthquakes[J].Bull Seism Soc Amer,2007,97(1):249-270.

[9]高祥林.琉球海溝的構(gòu)造和運動特征[J].地球物理學(xué)進展,2003,18(2):293-301.

GAO Xiang-lin.Tectonics and Motion of the Ryukyu Trench[J].Progress in Geophysics,2003,18(2):293-301.(in Chinese).

[10]臧紹先,寧杰遠,許立忠.琉球島弧地區(qū)的地震分布、Beniof帶形態(tài)及應(yīng)力狀態(tài)[J].地震學(xué)報,1989,11(2):113-123.

ZANG Shao-xian,NING Jie-yuan,XU Li-zhong.The Distribution of Earthquakes,Behavior of the Subduction Zone and Stress State in the Ryukyu Island Arc[J].Acta Seismologica Sinica,1989,11(2):113-123.(in Chinese)

[11]臧紹先,寧杰遠.西太平洋俯沖帶的研究及其動力學(xué)意義[J].地球物理學(xué)報,1996,39(2):188-202.

ZANG Shao-xian,NING Jie-yuan.Study on the Subduction Zone in Western Pacific and Its Implication for the Geodynamics[J].Chinese Journal of Geophysics,1996,39(2):188-202.(in Chinese)

[12]皇甫海泉,魏東平.沖繩板塊應(yīng)力場數(shù)值模擬及其動力學(xué)特征[J].地震學(xué)報,2010,32(5):529-543.

HUANGFU Hai-quan,WEI Dong-ping.Numerical Simulation of the Stress Field in Okinawa Plate and Its Dynamic Characteristics[J].Acta Seismologica Sinica,2010,32(5):529-543.(in Chinese)

[13]李家彪.中國邊緣海形成演化與資源效應(yīng)[M].北京:海洋出版社,2005.

LI Jia-biao.Evolution and Resources Effect of Chinese Marginal Seas[M].Beijing:China Ocean Press,2005.(in Chinese)

[14]Subarya C,Chlieh M,Prawtrodirdjo L,et al.Plate-boundary Deformation Associated with the Great Sumatra—Andaman Earthquake[J].Nature,2006,440(2):46-51.

[15]Satake K,Tanioka Y.Sources of Tsunami and Tsunamigenic Earthquakes in Subduction Zones[J].Pure appl Geophys,1999,154:467-483.

[16]Mansinha L,Smylie D E.The Displacement Fields of Inclined Faults[J].Bull Seism Soc Amer,1971,61:1433-1440.

[17]Geist E L,Bilek S L,Arcas D,et al.Differences in Tsunami Generation between the December 26,2004 and March 28,2005 Sumatra Earthquakes[J].Earth Planets Space,2006,58:185-193.

[18]Simon M D,Mcneill L C,Henstock T J,et al.Contrasting Décollement and Prism Properties over the Sumatra 2004—2005 Earthquake Rupture Boundary[J].Science,2010,329:207-210.

[19]Wang X M,Liu P L-F.An Analysis of 2004 Sumatra Earthquake Fault Plane Mechanisms and Lndian Ocean Tsunami[J].Journal of Hydraulic Research,2006,44(2):147-154.

[20]陳鋒.我國濱海核電廠廠坪標(biāo)高的確定[J].核安全,2006(2):44-48.

CHEN Feng.How to Determine the Platform Elevation of Nuclear Power Plants on Coastal Sites in China[J].Nuclear Safety,2006(2):44-48.(in Chinese)

[21]張愛玲.對濱海核電廠防洪評價中海嘯影響的一些認識[J].核安全,2008(4):28-32.

ZHANG Ai-ling.Discussion of Tsunami Impact in Flooding Assessment of Nuclear Power Plant on Coastal Site[J].Nuclear Safety,2008(4):28-32.(in Chinese)

Numerical Simulation of Risk in the Southeast Coastal Region of China Owing to a Rare Earthquake-induced Tsunami in the Ryukyu Trench

WEN Yan-ling1, LIU Shuang-qing2, ZHU Ai-lan2, SONG Zhi-ping3, LI Chun-feng1

(1.SchoolofOceanandEarthScience,TongjiUniversity,Shanghai200092,China;2.EarthquakeAdministrationofShanghaiMunicipality,Shanghai201203,China;3.ChinaEarthquakeNetworksCenter,Beijing100045,China)

Abstract:Based on their tectonic similarities,the Ryukyu trench is as likely to experience a rare mega-thrust earthquake as the Japanese,Chilean,and Sunda trenches.Based on an analysis of the mega-thrust hypocenter parameters along trenches in recent years,this study proposed a scenario of an earthquake of MW9.0 in the Ryukyu trench and then simulated the resulting earthquake-induced tsunami.Our study results indicated that the initial tsunami reached a height of up to 8 m,then propagated to northeast Taiwan half an hour later and increased in height to over 10 m.Three to four hours later,the tsunami reached the south coast of Zhejiang province and the north coast of Fujian province,and its height was then 1～2 m.Roughly five hours later,the tsunami propagated to the north coasts of Zhejiang and Guangdong provinces,at heights of approximately 2 m and less than 0.5 m,respectively.The tsunami reached the shoreline of Shanghai eight hours later,and its maximum height was then approximately 1 m.These simulated tsunami height data tend to be lower than would occur in reality because of the shoaling effect in the bays along the southeast China shoreline.Owing to the presence of a number of nuclear power plants along the coast of southeast China,the potential for tsunami impact on these plants should not be ignored in the event of a rare earthquake in the Ryukyu trench.We recommend that nuclear plants take measures to prevent potential damage from tsunami hazards.

Key words:Ryukyu trench; rare earthquake; tsunami simulation; tsunami risk map; nuclear power plant safety

DOI:10.3969/j.issn.1000-0844.2016.02.0285

中圖分類號:TU47

文獻標(biāo)志碼:A

文章編號:1000-0844(2016)02-0285-08

作者簡介:溫燕林(1981-),男,在讀博士生,主要研究方向為地球動力學(xué)及海嘯模擬。E-mail:wenzhengyi@sina.com。

基金項目:國家自然科學(xué)基金項目(41372222);上海市科委科研計劃項目(14231202600)

收稿日期:①2015-02-22