基于震源機(jī)制解的分區(qū)地震烈度衰減關(guān)系研究1

鄭　韻姜立新楊天青劉　杰1）福建省地震局，福州 350003中國地震臺網(wǎng)中心，北京 100045

基于震源機(jī)制解的分區(qū)地震烈度衰減關(guān)系研究1

鄭韻1）姜立新2）楊天青2）劉杰2）
1）福建省地震局，福州 350003
2）中國地震臺網(wǎng)中心，北京 100045

鄭韻，姜立新，楊天青，劉杰，2016.基于震源機(jī)制解的分區(qū)地震烈度衰減關(guān)系研究.震災(zāi)防御技術(shù)，11（2）：349—359. doi：10.11899/zzfy20160218

根據(jù)中國大陸構(gòu)造應(yīng)力場分區(qū)，本文擬合不同震源機(jī)制解的烈度衰減關(guān)系，收集了我國西部198次5.0級以上震例，共419條的等震線記錄；東部47次5.0級以上震例，共84條等震線記錄，采用橢圓模型進(jìn)行衰減關(guān)系擬合，得到了地震烈度隨震級和長短軸長度的衰減關(guān)系。文章初步分析了幾種震源機(jī)制地震烈度衰減規(guī)律的差異，以西部地區(qū)為例，對比了西部走滑型、逆斷型烈度衰減關(guān)系和汪素云西部衰減關(guān)系與實(shí)際值之間的差距。本文研究結(jié)果表明，區(qū)分發(fā)震方式的烈度衰減關(guān)系長短軸可以作為應(yīng)急地震影響場的修正因子，在震后幾小時(shí)得到震源機(jī)制解時(shí)，對烈度圈長短軸進(jìn)行初步修正。

烈度衰減關(guān)系震源機(jī)制解構(gòu)造應(yīng)力場分區(qū)

引言

地震烈度衰減關(guān)系在抗震設(shè)防和震后烈度快速評估等工作中，具有重要的實(shí)際意義。我國幅員遼闊，地質(zhì)結(jié)構(gòu)復(fù)雜，很多研究者已對不同地區(qū)、不同地質(zhì)構(gòu)造環(huán)境的地震烈度衰減關(guān)系進(jìn)行了研究，得出了適用于各個(gè)地區(qū)的統(tǒng)計(jì)關(guān)系。汪素云等（2000）將中國分為東、西部兩個(gè)區(qū)，建立了分區(qū)的烈度衰減關(guān)系；有些研究人員考慮了地質(zhì)構(gòu)造環(huán)境和震害分布特點(diǎn)，將川滇地區(qū)分為西南地區(qū)和盆地地區(qū)，分別確定了烈度衰減規(guī)律（雷建成等，2007；孫繼浩等，2011）；譚明等（2011）分不同時(shí)段和區(qū)域，得到了符合新疆地區(qū)的衰減關(guān)系。但是，除近地表環(huán)境對地表強(qiáng)震動有很大影響外，震源破裂方式也會對烈度的衰減特性造成影響。以往的研究中，缺乏考慮因地震破裂方式不同造成的烈度衰減規(guī)律的差別。本文將根據(jù)中國大陸構(gòu)造應(yīng)力場分區(qū)（謝富仁等，2004），回歸擬合不同震源機(jī)制解的烈度衰減關(guān)系，為地震應(yīng)急救援提供參考。

1　大陸構(gòu)造應(yīng)力場分區(qū)和震例選取

在區(qū)域構(gòu)造應(yīng)力作用下地殼應(yīng)變能積累到一定程度時(shí)造成的斷層錯(cuò)動，就是地震。區(qū)域應(yīng)力的強(qiáng)度提供了地震可能釋放能量的信息，應(yīng)力的特性顯示了可能的地震類型。地震釋放的能量和地震類型是影響烈度衰減關(guān)系和震害分布的重要因素。謝富仁等（2004）根據(jù)中國大陸構(gòu)造應(yīng)力場的基本特征，劃分了4個(gè)二級應(yīng)力區(qū)，分別為東北—華北應(yīng)力區(qū)，華南應(yīng)力區(qū)，新疆應(yīng)力區(qū)和青藏高原應(yīng)力區(qū)。考慮到地震資料的多少和研究區(qū)域范圍對地震烈度衰減關(guān)系合理性的影響，本文將青藏高原應(yīng)力區(qū)進(jìn)一步分為川滇地區(qū)和甘青寧地區(qū)進(jìn)行研究，并將東北—華北應(yīng)力區(qū)和華南應(yīng)力區(qū)合并為東部應(yīng)力區(qū)進(jìn)行研究。

本文震例資料主要取自《中國近代地震目錄（公元1912年—1990年）》、《中國震例》（1966年—2007年）和《中國大陸地震災(zāi)害損失評估》（1990年—2014年）等。為避免I和M之間相關(guān)性的影響，本文采用的烈度等震線資料的I和M均為獨(dú)立測定。

震例選取時(shí)，若資料中給出等震線長短軸半徑，則直接采用；若沒有給出，數(shù)字化后采用長軸可轉(zhuǎn)向方法進(jìn)行量取。本文震源機(jī)制解結(jié)果來自美國哈佛大學(xué)的CMT解、中國地震局及其下屬機(jī)構(gòu)解析處理結(jié)果（《中國震例》，1966年—2007年；《中國大陸地震災(zāi)害損失評估》，1990年—2014年）。若存在不同震源機(jī)制解，以哈佛大學(xué)CMT解為準(zhǔn)。本文將震源機(jī)制按照走滑型、逆斷-逆斷走滑型（以下簡稱逆斷型）、正斷-正斷走滑型（以下簡稱正斷型）和未確定型四種類型進(jìn)行烈度衰減關(guān)系分類統(tǒng)計(jì)。依據(jù)上述原則，在已知震源機(jī)制解的條件下，本文選取了我國西部198次5.0級以上震例，共419條的等震線記錄；我國東部47次5.0級以上震例，共84條等震線記錄。表1、表2給出了所選用地震資料的震源機(jī)制解分布。

表1　中國大陸震源機(jī)制解分布Table 1　The geographical distribution of focal mechanism solution in Chinese mainland

表2　西部分區(qū)震源機(jī)制解分布Table 2　The geographical distribution of focal mechanism solution in western China

2　地震烈度衰減關(guān)系

由于受到發(fā)震構(gòu)造的影響，近場等震線一般呈橢圓分布，遠(yuǎn)場逐漸衰減趨于圓形。本文采用橢圓模型進(jìn)行衰減關(guān)系擬合。

式中，c1、c2、c3、a0為回歸系數(shù)，I表示烈度，a是烈度I的等震線的長軸或短軸長度，Ms表示面波震級?；貧w系數(shù)a0的選取原則是將a0取一系列初值進(jìn)行計(jì)算，選取標(biāo)準(zhǔn)差最小的值為最終結(jié)果。

極震區(qū)內(nèi)任意一點(diǎn)具有相同的烈度值，筆者在極震區(qū)不同距離進(jìn)行適當(dāng)?shù)慕鼒鲅a(bǔ)點(diǎn)，補(bǔ)點(diǎn)僅在震中烈度7度以上及最內(nèi)圈等震線半徑r＞5km時(shí)進(jìn)行。另外，遠(yuǎn)場發(fā)震構(gòu)造的影響逐漸消失，等震線趨于圓形，取有感半徑作為遠(yuǎn)場控制點(diǎn)，有感烈度值為4度時(shí)進(jìn)行遠(yuǎn)場補(bǔ)點(diǎn)。有感半徑R和震級M的關(guān)系利用下式計(jì)算（中國地震烈度區(qū)劃圖編委會，1992）：

本文對選取資料基于最小二乘法，對震例頻次≥15的資料進(jìn)行擬合，得到地震烈度隨震級和長短軸長度的衰減關(guān)系。表3至表7給出了中國大陸東、西部以及西部各分區(qū)的地震烈度衰減關(guān)系回歸系數(shù)。

表3　中國東部地震烈度衰減系數(shù)Table 3　Seismic intensity attenuation coefficients in eastern China

表4　中國西部地震烈度衰減系數(shù)Table 4　Seismic intensity attenuation coefficients in western China

表5　新疆應(yīng)力區(qū)地震烈度衰減系數(shù)Table 5　Seismic intensity attenuation coefficients in Xinjiang stress zone

表6　甘青寧區(qū)地震烈度衰減系數(shù)Table 6　Seismic intensity attenuation coefficients in the Gansu-Qinghai-Ningxia areas

表7　川滇區(qū)地震烈度衰減系數(shù)Table 7　Seismic intensity attenuation coefficients in the Sichuan-Yunnan stress zone

2.1各分區(qū)回歸結(jié)果與分析

統(tǒng)計(jì)本文收集的震源機(jī)制解數(shù)據(jù)和烈度資料可知，中國大陸5級以上地震中，走滑斷層活動占58%，逆斷型占31%，僅有10%為正斷型斷層活動的結(jié)果。圖1至圖5分別為“東部”、“西部”、“新疆”、“甘青寧”、“川滇”不同震源機(jī)制解的烈度衰減曲線。從圖中可以看出，不同震源錯(cuò)動方式的烈度衰減曲線總體趨勢一致，但不同地區(qū)不同錯(cuò)動方式具有不同的衰減特征。

東部地區(qū)：近場范圍內(nèi)，走滑型地震烈度略低于全部震例，中遠(yuǎn)場的烈度值幾乎重合，這與東部地區(qū)走滑型選例占72%有很大關(guān)系。

西部地區(qū)：近場范圍內(nèi)，M=7、8時(shí)，地震烈度正斷型＞逆斷型＞全部震例＞走滑型；M=5、6時(shí)，逆斷型地震烈度仍然大于走滑型，但正斷型烈度隨震級衰減快于其它類型，5級正斷型地震衰減曲線略低于其它幾條曲線。隨著距離的增加，不同類型地震烈度差異縮小。遠(yuǎn)場范圍內(nèi)，長軸方向走滑型地震影響范圍最大，正斷型影響范圍最小，逆斷型居中；短軸方向不同類型地震影響范圍趨于一致。

圖1　東部不同震源機(jī)制解的烈度衰減關(guān)系對比Fig.1　Comparison of seismic intensity attenuation relationships with different focal mechanism solution in eastern China

圖2　西部不同震源機(jī)制解的烈度衰減關(guān)系對比Fig.2　Comparison of seismic intensity attenuation relationships with different focal mechanism solution in western China

圖3　新疆不同震源機(jī)制解的烈度衰減關(guān)系對比Fig.3　Comparison of seismic intensity attenuation relationships with different focal mechanism solution in Xinjiang

圖4　甘青寧不同震源機(jī)制解的烈度衰減關(guān)系對比Fig.4　Comparison of seismic intensity attenuation relationships with different focal mechanism solution in Gansu-Qinghai-Ningxia

圖5　川滇不同震源機(jī)制解的烈度衰減關(guān)系對比Fig.5　Comparison of seismic intensity attenuation relationships with different focal mechanism solution in Sichuan-Yunnan

新疆地區(qū)：近場范圍內(nèi)，逆斷型地震烈度大于走滑型，隨距離的增加兩者烈度逐漸接近。遠(yuǎn)場范圍內(nèi)，長軸方向走滑型震例影響范圍略大于逆斷型，即逆斷型烈度衰減曲線衰減速度快于走滑型。全部震例衰減曲線位于兩者之間。

甘青寧地區(qū)：烈度衰減曲線整體情況和新疆地區(qū)接近。近場范圍內(nèi)，逆斷型地震烈度大于走滑型，隨距離的增加兩者烈度逐漸接近。遠(yuǎn)場范圍內(nèi)，長軸方向走滑型震例影響范圍略大于逆斷型，即逆斷型烈度衰減曲線衰減速度快于走滑型。全部震例衰減曲線位于兩者之間。

川滇地區(qū)：走滑型地震在震級大于7時(shí)，衰減曲線整體低于全部震例擬合曲線，6級左右時(shí)兩者衰減曲線幾乎重合，5級走滑型地震烈度高于全部震例，也就是說，走滑型大地震破壞小一些，隨震級下降烈度衰減速率慢于全部震例。

2.2與前人結(jié)果的對比

我國地震等震線資料豐富，許多學(xué)者都已對不同區(qū)域的烈度衰減關(guān)系進(jìn)行過研究，得到了符合區(qū)域特點(diǎn)的結(jié)果。圖6至圖10分別為本文所得到的各地區(qū)的擬合曲線與汪素云等（2000）、肖亮等（2001）、譚明等（2001）、周中紅等（2010）、雷建成等（2007）和孫繼浩等（2011）所給出的幾種擬合曲線的對比。

本文所得的東部衰減曲線在震級為7級以上時(shí)略低于曲線a，中強(qiáng)震級地震烈度較為一致。西部地區(qū)的衰減曲線和曲線b、曲線c的結(jié)果總體趨勢一致，曲線b影響范圍稍大于本文的結(jié)果和曲線c，曲線c的衰減速率快于本文的結(jié)果和曲線b。新疆地區(qū)長軸方向，本文所得的結(jié)果在近場范圍內(nèi)略低于曲線e和曲線d，中遠(yuǎn)場位于兩者之間；短軸方向近場范圍內(nèi)本文曲線和曲線d接近，均高于曲線e，中遠(yuǎn)場烈度值高于其它兩者。甘青寧全震例衰減曲線的長軸方向在近場范圍時(shí)，本文曲線位于曲線g和曲線f之間，中遠(yuǎn)場烈度值低于曲線g和曲線f的結(jié)果；短軸方向近場處甘青寧全震例和曲線g的烈度大于曲線f的結(jié)果，中遠(yuǎn)場范圍內(nèi)甘青寧全震例的烈度值略低于其它兩條曲線。川滇地區(qū)本文擬合曲線的衰減速率快于曲線j，慢于曲線i的計(jì)算結(jié)果，曲線j的地震影響范圍最大，曲線i擬合結(jié)果的地震影響范圍較小，本文結(jié)果影響范圍居于二者之間。這主要由于孫繼浩沒有采用遠(yuǎn)場補(bǔ)點(diǎn)，且選用資料為6級以上震例，而雷建成遠(yuǎn)場補(bǔ)點(diǎn)依據(jù)四川有感半徑與震級關(guān)系，且選用資料為4級以上震例所造成的差異。同一地區(qū)衰減曲線差異主要是震級選取、烈度資料分布以及是否進(jìn)行遠(yuǎn)場補(bǔ)點(diǎn)等因素制造而成。

圖6　中國東部烈度衰減關(guān)系對比Fig.6　Comparison of seismic intensity attenuation relationships in eastern China

圖7　中國西部烈度衰減關(guān)系對比Fig.7 Comparison of seismic intensity attenuation relationships in western China

圖8　新疆地區(qū)烈度衰減關(guān)系對比Fig.8 Comparison of seismic intensity attenuation relationships in Xinjiang

圖9　甘青寧地區(qū)烈度衰減關(guān)系對比Fig.9 Comparison of seismic intensity attenuation relationships in Gansu-Qinghai-Ningxia

圖10　川滇地區(qū)烈度衰減關(guān)系對比Fig.10　Comparison of seismic intensity attenuation relationships in Sichuan-Yunnan

為進(jìn)一步探討地震破裂方式對烈度衰減的影響，本文比較了西部地區(qū)走滑型、逆斷型地震的烈度衰減關(guān)系，以及汪素云西部衰減關(guān)系計(jì)算的Ⅵ至Ⅸ度影響場半徑與實(shí)際結(jié)果的差距，長度差均值如表8、表9所示：

表8　走滑型長短軸與實(shí)際值長度差均值計(jì)算結(jié)果Table 8　The mean deviation value of major and minor axis length of strike-slip faults between calculated and practical value

表9　逆斷型長短軸與實(shí)際值長度差均值計(jì)算結(jié)果Table 9　The mean deviation value of major and minor axis length of thrust faults between calculated and practical value

從表8、表9可以看出，本文所得結(jié)果除走滑型Ⅸ度長軸和逆斷型Ⅵ度短軸略差于汪素云西部之外，其它情況劃分震源機(jī)制的地震烈度衰減關(guān)系計(jì)算平均結(jié)果均優(yōu)于汪素云西部，長軸提升率在-6.26%—9.96%之間，短軸提升率在-0.87%—16.63%之間。總體來說，區(qū)分發(fā)震方式的烈度衰減關(guān)系比未區(qū)分的更接近真實(shí)情況，進(jìn)一步佐證了地震破裂方式對烈度衰減的影響，這將有助于我們更加準(zhǔn)確地估計(jì)影響場的大致范圍，對于震后災(zāi)害評估，震后應(yīng)急救援具有一定指導(dǎo)意義。

3　討論和結(jié)論

發(fā)震方式對地震烈度衰減關(guān)系的影響是較為顯著的?？傮w來看，震級相同時(shí)走滑型地震的烈度衰減曲線值是最小的，逆斷型以及正斷型衰減曲線烈度值是最大的，全部震例的擬合曲線介于二者之間。長軸方向走滑型地震影響范圍似乎略大于其它類型，但是地震的不同類型在短軸方向的影響范圍差異不大?？梢钥闯?，同一區(qū)域不同震源機(jī)制的地震烈度衰減特征有差別，不同區(qū)域相同震源機(jī)制的烈度衰減關(guān)系并不完全一致，震源機(jī)制對不同震級的烈度影響也不同。此外，地震烈度衰減規(guī)律的影響因素是復(fù)雜的，除震源特性之外，地震波傳播過程中非彈性和非均勻性導(dǎo)致的衰減主要取決于品質(zhì)因子Q值，而Q值是區(qū)域性變量。另外，局部場地條件也會對烈度衰減造成顯著影響。

本文根據(jù)中國大陸構(gòu)造應(yīng)力場分區(qū)，分別擬合不同震源機(jī)制解的烈度衰減關(guān)系，得出以下結(jié)論：

（1）建立了中國大陸不同震源機(jī)制解的分區(qū)地震烈度衰減關(guān)系。

（2）震源機(jī)制對地震烈度衰減關(guān)系的影響顯著?？傮w來看，逆斷型、正斷型地震的烈度衰減曲線烈度值大于走滑型。

（3）震源機(jī)制對不同震級的地震烈度影響程度不同。

（4）區(qū)分發(fā)震方式的烈度衰減關(guān)系長短軸可以作為應(yīng)急影響場的修正因子，在震后幾小時(shí)得到震源機(jī)制解時(shí)，對烈度圈長短軸進(jìn)行初步修正。

雷建成，高孟潭，俞言祥，2007.四川及鄰區(qū)地震動衰減關(guān)系.地震學(xué)報(bào)，29（5）：500—511.

孫繼浩，帥向華，2011.川滇及其鄰區(qū)中強(qiáng)地震烈度衰減關(guān)系適用性研究.地震工程與工程振動，31（1）：11—18.

譚明，李帥，孫靜等，2011.新疆地震烈度衰減關(guān)系模型參數(shù)擬合.內(nèi)陸地震，25（1）：29—35.

汪素云，俞言祥，高阿甲等，2000.中國分區(qū)地震動衰減關(guān)系的確定.中國地震，16（2）：5—12.

肖亮，俞言祥，2011.中國西部地區(qū)地震烈度衰減關(guān)系.震災(zāi)防御技術(shù)，6（4）：358—371.

謝富仁，崔效鋒，趙建濤等，2004.中國大陸及鄰區(qū)現(xiàn)代構(gòu)造應(yīng)力場分區(qū).地球物理學(xué)報(bào)，47（4）：654—662.

中國地震烈度區(qū)劃圖編委會，1992.中國地震烈度區(qū)劃圖（1990）及其說明書.中國地震，8（4）：1—11.

周中紅，何少林，陳文凱，2010.甘肅省地震烈度衰減關(guān)系研究.西北地震學(xué)報(bào)，31（1）：72—75.

Study on Seismic Intensity Attenuation Relationship with Regions via Focal Mechanism Solutions

Zheng Yun1），Jiang Lixin2），Yang Tianqing2）and Liu Jie2）
1）Earthquak Administration of Fujian province，F(xiàn)uzhou 350003，China
2）China Earthquake Networks Center，Beijing 100045，China

On the basis of regional division of tectonic stress field in China，we fit the seismic intensity attention relationships with different focal mechanism solution.The data we have collected includes 198 earthquakes and 419 isoseismal lines with MS≥5.0 in western China，altogether；47 earthquakes and 84 isoseismal lines with MS≥5.0 in eastern China.By using the ellipse model fitting attenuation relation，seismic intensity attenuation relationship between the magnitude and the length of axis is established.In this paper，we preliminarily analyze the discrepancy between seismic intensity attenuation relationships of several focal mechanism solutions.Then，taking western China as an example，we compare western seismic intensity attenuation relationship of thrust fault and strike-slip fault with Wang Suyun’s results.Our research results indicate that the length of axis of intensity attenuation relationship which have distinguished focal mechanism solutions could be a modifying factor for seismic influence field in case of emergency.When we obtain focal mechanism solution several hours after the earthquake，we can make a preliminary correction to the length of axis.

Intensity attenuation relationship；Focal mechanism solution；Regional division of tectonic stress field

地震科技星火計(jì)劃項(xiàng)目（XH14056）

2015-06-21

鄭韻，女，生于1989年。碩士研究生。主要研究方向?yàn)榈卣饝?yīng)急技術(shù)研究。E-mail：kayu1989@sohu.com

姜立新，男，生于1966年。研究員。主要從事震害預(yù)測、地震應(yīng)急技術(shù)等研究。E-mail：jlx@seis.ac.cn