亢川川，李建亮，孔 軍，劉 韶

(1.四川省地震局工程地震研究院，四川 成都 610041；2.四川賽思特科技有限責任公司，四川 成都 610041)

2015年4月25日尼泊爾Ms8.1級地震震動圖的不確定性分析

亢川川1,2，李建亮1，孔 軍1，劉 韶2

(1.四川省地震局工程地震研究院，四川 成都 610041；2.四川賽思特科技有限責任公司，四川 成都 610041)

尼泊爾Ms8.1級地震造成了重大的人員傷亡和財產損失。中國地震局在震后2小時就及時發(fā)布了此次地震的地震動強度預測圖(即震動圖)。作為地震應急救援的重要圖件之一，震動圖對地震救援和震后災害損失評估等工作提供了科學有效的參考。但是在震后短時間內快速生成的震動圖難以避免地存在一些不確定性因素，需要根據震后信息的不斷積累及時加以修正，同時還要結合考慮局部地質構造背景、地震波的衰減關系特性以及土層對地震動參數放大效應等方面的不確定性因素來對于強度預測圖進行修正。

震動圖；不確定性；分析

2015年4月25日14時11分(北京時間)在尼泊爾發(fā)生了8.1級地震。據中國地震臺網中心測定，震中位置為28.2°N、84.7°E，震源深度20 km[1]。地震震中位于尼泊爾的博克拉，最大烈度為X度，等震線長軸總體呈北西西走向。截至2015年5月13日，地震至少造成8 216人死亡，近2萬人受傷。地震發(fā)生后又發(fā)生了多次強余震，其中7級以上強余震就有3次。當地時間2015年5月12日發(fā)生的尼泊爾加德滿都7.5強烈地震導致65人遇難，另有約1 200人受傷[2]。尼泊爾發(fā)生了8.1級地震后不久，中國地震局先后發(fā)布了相應的地震動強度預測圖、地震烈度圖。相比地震烈度圖，震動圖的優(yōu)勢在于快速產出，此次地震震動圖對地震動強度的預測基本可靠，但仍存在一定的誤差，如何盡可能地改進圖的精度是本文的目的。

1 地震動強度預測圖和地震烈度圖基本情況

中國地震局在地震發(fā)生后2小時迅速發(fā)布了“尼泊爾8.1地震地震動強度預測圖”(參見圖1)[3]，以下簡稱“震動圖”。綜合考慮了現有的局部地質構造背景、地震波的衰減關系特性以及土層對地震動參數放大效應的影響，給出了此次地震的震動圖?！罢饎訄D”是根據場地效應在很短的時間里面快速生成的[4]，其中地震波的衰減采用了汪素云2000的中國西部衰減模型，以點橢圓衰減模型計算了此次地震的影響場，并且考慮了土層對地震動參數放大效應。而這里的地震動參數放大效應是采用了Allen等于2007年和Wald等于2006和2007分別提出的一種直接利用地形數據確定場地類別的方法計算的，依據這些計算參數對此次地震的場地條件進行了估計[5-7]。在此基礎上，還利用了Borcherdt的場地放大系數成果參加計算[8]，所有這些工作量化了土層對于地震動參數的放大效應。計算結果表明，這次地震預測極震區(qū)的地震烈度可能達到Ⅹ度，可能的受災面積達到20萬平方千米，中國西藏吉隆鎮(zhèn)受影響的烈度可能達到Ⅷ度。隨后基于收集到的震源破裂過程信息，利用斷層投影面源的衰減模型修正了最初版本的震動圖，大大優(yōu)化了地震震動圖的精度。

隨著時間的推移，尼泊爾8.1級地震發(fā)生后6天，中國地震局于2015年5月1日發(fā)布了“尼泊爾8.1級地震烈度圖”[9]，以下簡稱“烈度圖”。這次地震波及了尼泊爾、中國、印度、孟加拉等國。中國地震局又根據災情、實地調查數據、觀測數據、衛(wèi)星遙感震害解譯和全球有關機構數據，并參照《中國地震烈度表》(GB/T 17742-2008)，進一步對此次地震烈度分布進行了評估。此次地震災區(qū)最高烈度為Ⅸ度，等震線長軸總體走向呈北西西方向，Ⅵ度區(qū)及以上的總面積約為214 700平方千米；Ⅸ度區(qū)及以上面積約8 300平方千米，長軸長155千米，短軸長63千米；Ⅷ度區(qū)面積約為20 500平方千米，長軸長260千米，短軸長135千米；Ⅶ度區(qū)面積約為45 000平方千米，長軸長383千米，短軸長236千米；Ⅵ度區(qū)面積約為140 900平方千米，長軸長588千米，短軸長470千米。其中，地震烈度圖中國境內是根據中國地震局現場工作隊實地調查數據進行評定和繪制的；中國境外地震烈度圖是根據震源機制解、震源破裂過程、地質構造、建筑物破壞特征等地球物理、地質、工程結構等分析資料，并結合衛(wèi)星遙感震害解譯和中國國際救援隊調查數據進行估計和繪制的。

圖1 2015年4月25日尼泊爾8.1級地震動強度預測圖

圖2 尼泊爾8.1地震地烈度圖

2 討論分析

從尼泊爾8.1級地震于地震后2小時發(fā)布的“震動圖”和于震后6天發(fā)布的“烈度圖”的對比分析來看，震動圖在極震區(qū)最大可能烈度、受災面積以及對于中國受災最嚴重的西藏吉隆鎮(zhèn)的估計是準確的。說明震動圖生成的方法科學可行，這些圖件的及時發(fā)布對地震救援和震后災害損失評估等工作提供了科學有效的參考依據。但是在Ⅶ、Ⅷ、Ⅸ度的面積、長短軸長度等參數方面，震動圖的估計在發(fā)布初期仍存在一定的誤差，尤其對極震區(qū)的估計，誤差更大。其原因在于震后短時間內可獲取的“震動圖”相應的資料十分有限，導致問題的邊界條件過少，這給精準估計地震動強度帶來了很大的不確定性，從而導致“震動圖”可能出現較大的誤差?！罢饎訄D”的精度受到局部地質構造背景、地震波的衰減關系特性以及土層對地震動參數放大效應等多方面因素的影響，而在尼泊爾地區(qū)相關基礎資料缺乏的情況下，就可能導致地震動強度估計上的誤差。在局部地質構造背景下，中國地震局和美國地質調查局[10]都給出了這次地震的震源機制解，兩個結果大體一致，所以大家認為此次地震是一次以逆沖型為主的斷層活動所產生。一般逆沖型的斷層活動產生的地震，上盤地震動強度一般要高于下盤，這點在這次發(fā)布的“烈度圖”上得到了體現，也就是說斷層南側的地震烈度高于北側。對此次以逆沖為主的地震事件在局部地質構造背景上估計不足，可能會導致估計結果存在一定的不確定性。

對于地震波的衰減問題，這次“震動圖”采用了點橢圓模型來表達。地震發(fā)生以后，得到了更多的資料和信息，有關機構隨后發(fā)布了這次地震的震源破裂過程信息，由此認為采用線源或者面源衰減模型更符合實際。另外，“震動圖”采用的衰減關系在某種程度上不一定適用于尼泊爾地區(qū)。在衰減關系方面，我國許多研究者做過大量的研究，根據情況的不同給出了相應區(qū)域的衰減關系[11-14]。然而，這些衰減關系各有特點，從資料到模型都有差別，必然就有一定的區(qū)域性特征。目前，中國的地震學研究者對此次地震影響區(qū)或尼泊爾地區(qū)的地震波衰減關系特性研究較少，選擇適合該地區(qū)的衰減關系是十分重要的，而選擇的本身就存在一定的不確定性，還有很多值得研究的課題。

在土層對地震動參數放大效應上，這次地震的“震動圖”采用了Allen等和Wald等分別于2006、2007年提出的地形坡度與Vs30(地下30m的平均剪切波速數據)相關關系來確定場地條件的方法計算相應參數。該方法提出后，被廣泛地應用到震后的地震動強度分布的快速預測工作中[15-20]。這次地震的“震動圖”根據USGS全球地形數據計算了坡度等參數，進一步推導后得到了此次地震受災區(qū)Vs30參數，從而獲得了場地條件圖。之后，運用Borcherdt的場地放大系數[8]，得到考慮局部場地效應影響后的圖件，并將基巖上的地震動參數分布修正到土層地表上，最終得到了考慮土層對地震動參數放大效應后的“震動圖”。然而，地形坡度與Vs30關系也有一定的區(qū)域性特征，Allen等統(tǒng)計的是美國加利福尼亞、尤他州，意大利和中國臺灣等活動構造區(qū)，而美國孟菲斯、澳大利亞等穩(wěn)定構造區(qū)得到的統(tǒng)計關系哪類更適合尼泊爾地區(qū)是需要進一步探討的。當然，通過尼泊爾地區(qū)的Vs30實測數據所得到的關系肯定更合適該地區(qū)。同樣的，Borcherdt等的場地放大系數同樣存在區(qū)域特性，選取其參數也會對于相應的地震烈度結果造成相應誤差。

3 結束語

尼泊爾Ms8.1級地震后2小時，中國地震局發(fā)布了此次地震的地震動強度預測圖，對地震救援與震后災害損失評估等工作提供了科學有效的參考。其估計結果基本可靠，準確性可通過修正不斷提高。震后短時間內可獲取的震動圖估計所需的相關資料十分有限，在追求實效性的同時，客觀降低了準確性。比如在衰減關系方面，對地震周邊研究程度比較欠缺的地區(qū)，參考研究程度較高地區(qū)的衰減相關特征，通過對模型和衰減特點的修正預先給出適合該地區(qū)的衰減關系，作為數據儲備以備今后地震應急救援時使用；再者，土層對地震動參數放大效應不容忽視，應大力開展區(qū)域性的場地條件研究工作，以便更為準確地量化場地效應。另外，可利用震源破裂過程、強震記錄等及時修正地震動強度預測圖，提高準確性。最后，此次地震主震發(fā)生后，又連續(xù)發(fā)生多次高震級余震，又造成了人員傷亡，同時對建筑物等的破壞有一定累加作用。建議疊加主震和余震“震動圖”，給地震烈度圖以參考，以更好地服務于震后災害損失評估工作。

[1] 中國地震局，2015，中國地震局發(fā)布尼泊爾發(fā)生8.1級地震.[2015-4-25].http://www.cea.gov.cn/publish/dizhenj/468/553/101803/101809/20150425152150837617854/index.html.

[2] 中國新聞網，2015，尼泊爾8.1級強地震遇難人數升至8019人[EB/OL],http://www.hi.chinanews.com/hnnew/2015-05-11/4_48694.html.

[3] 中國地震局，2015，中國地震局發(fā)布2015年4月25日尼泊爾8.1地震地震動強度預測圖.http://www.csi.ac.cn/manage/eqDown/31RollingTop/shakemap.html.

[4] 陳鯤，俞言祥，高孟潭.考慮場地效應的Shake Map系統(tǒng)研究[J].中國地震，2010，26(1)：92-102.

[5] Allen, T.I. and Wald, D.J., 2007, Topographic Slope as a Proxy for Seismic Site-Conditions (Vs30) and Amplification around the Globe: U.S. Geological Survey Open-File Report 2007-1357, 69p.

[6] Wald D J, Worden B C, Quitoriano V, and Pankow K L. 2006, Shake Map Manual: Technical manual, user’s guide, and software guide: U.S. Geological Survey, 156p.

[7] Wald D J，and Allen T I.Topographic Slope as a proxy for Seismic Site Conditions and Amplification，Bul1．Seism．Soc．Amer．2007，97(5):1379-1395．

[8] Borcherdt R D，1994，Estimates of site-dependent response spectra for design (methodology and justification)，Earthquake Spectra，10，617-654．

[9] 中國地震局，2015，中國地震局發(fā)布尼泊爾8.1級地震烈度圖.http://www.cea.gov.cn/publish/dizhenj/464/478/20150501221233264365835/index.html

[10] USGS,2015a, M7.8 -34 km ESE of Lamjung, Nepal[EB/OL].[2015-4-25],http://earthquake.usgs.gov/earthquakes/eventpage/us20002926#scientific_tensor:us_us_20002926_mwc.

[11] 俞言祥，李山有，肖亮.為新區(qū)劃圖編制所建立的地震動衰減關系[J]．震災防御技術，2013，8(1)：24-33.

[12] 王玉石，李小軍，周正華.川滇地區(qū)水平向強地震動衰減關系研究[J]．地震學報，2013，35(2)：238-249.

[13] 肖亮，俞言祥. 中國西部地區(qū)地震烈度衰減關系[J].震災防御技術，2011，6(4)：358-371.

[14] 雷建成，高孟潭，俞言祥.四川及鄰區(qū)地震動衰減關系[J].地震學報，2007，5(29)：500-511.

[15] 陳鯤，俞言祥，高孟潭.2010年4月14日青海玉樹地震震動圖[J].中國地震，2011，27(1)：99-102.

[16] 陳鯤，俞言祥，高孟潭，等.考慮震源破裂過程的青海玉樹地震震動圖研究[J].中國地震，2011，27(1)：56-64.

[17] 陳鯤，俞言祥，高孟潭，等.用有限強地震動記錄校正等震線的估計研究[J]．地震學報，2012，34(5)：633-645.

[18] 陳鯤，俞言祥，高孟潭.基于地震記錄的震動圖校正方法研究[J].應用基礎與工程科學學報, 2013，21(4)：679-691.

[19] 陳鯤，高孟潭.“抗大地震”與低超越概率水準地震動關系的討論[J].震災防御技術,2013，8(1)：34-39.

[20] 陳鯤，俞言祥，高孟潭.強震記錄校正的2013年4月月20日四川蘆山地震峰值加速度震動圖[J].地震地質，2013，35(3)：1-7.

Uncertainty analysis for the ShakeMaps ofMs8.1 Nepal Earthquake on 25 April 2015

Kang Chuanchuan1,2，Li Jianliang1，Kong Jun1，Liu Shao3

(1.Earthquake Administration of Sichuan Province, Sichuan Chengdu 610041；2.Sichuan Seistech Corporation, Ltd, Sichuan Chengdu 610041， China)

TheMs8.1 Nepal Earthquake caused heavy casualties and property losses. Two hours after the earthquake, China Earthquake Administration released the seismic intensity prediction map (i.e. ShakeMaps). The intensity map can provide a scientific and effective reference to earthquake rescue and earthquake damage assessment. But the rapid generation of ShakeMaps has many uncertain factors and we need to modify it according to the information we get after the earthquake. This paper discussed several uncertainty factors such as the local geological structure background, seismic wave attenuation relationship and soil amplification effect.

ShakeMaps；uncertainty；analysis

2015-05-18；

2015-07-30

四川省地震局地震科技專項LY1309“通過地形坡度與Vs30關系確定場地類別的研究”、四川賽思特科技有限責任公司董事長基金共同資助.

亢川川(1982-)，男，蒙古族，吉林省松原市人，工程師.E-mail:kangchuan@126.com.

P315.32

B

1001-8115(2015)03-0001-04

10.13716/j.cnki.1001-8115.2015.03.001