郭祥云 陳學(xué)忠

(中國(guó)地震局地球物理研究所,北京 100081)

2009年中國(guó)大陸地區(qū)部分地震的震源機(jī)制解＊

郭祥云 陳學(xué)忠

(中國(guó)地震局地球物理研究所,北京 100081)

選用國(guó)家地震臺(tái)網(wǎng)和區(qū)域臺(tái)網(wǎng)的初至波P波初動(dòng)符號(hào)資料,采用下半球等面積投影,用格點(diǎn)嘗試法計(jì)算了2009年中國(guó)大陸27次MS4.5級(jí)以上地震的震源機(jī)制解。根據(jù)Zoback[1]研究全球應(yīng)力場(chǎng)時(shí)的分類(lèi)標(biāo)準(zhǔn),對(duì)斷層面解類(lèi)型進(jìn)行了分類(lèi)。結(jié)果表明:走滑型(SS)地震13次,逆斷層為主兼走滑(TS)8次,逆斷層(TF)2次,正斷層為主兼走滑(NS)3次,正斷層(NF)1次。主壓應(yīng)力軸P軸方位的優(yōu)勢(shì)分布具有明顯的區(qū)域特征,除東北地區(qū)有稍微的偏移,其余地區(qū)都與當(dāng)?shù)氐谋尘皹?gòu)造應(yīng)力場(chǎng)基本一致。

初動(dòng)符號(hào);中國(guó)大陸地區(qū);震源機(jī)制

引言

震源機(jī)制解可以用于描述震源的性質(zhì)和受力情況,是了解地震應(yīng)力狀態(tài)的一個(gè)有效途徑,對(duì)監(jiān)視地殼應(yīng)力動(dòng)態(tài)變化和認(rèn)識(shí)地震孕育發(fā)生過(guò)程具有一定的意義。常用的求解震源機(jī)制的方法主要有P波初動(dòng)符號(hào)法、波形矩張量反演方法和振幅比方法[2-4],國(guó)內(nèi)外專(zhuān)家們已經(jīng)在這方面做了大量的工作[5-21]。

P波初動(dòng)符號(hào)法因其物理圖像簡(jiǎn)單,地震波信息穩(wěn)定,在確定震源機(jī)制解方面的獨(dú)特優(yōu)勢(shì),被國(guó)內(nèi)外學(xué)者廣泛關(guān)注。至今,已發(fā)展了多個(gè)利用P波初動(dòng)符號(hào)求解震源機(jī)制解的方法,大致可歸為兩類(lèi):迭代法[22-23]和網(wǎng)格搜索法。因迭代算法對(duì)資料質(zhì)量的依賴(lài)性更大,所以網(wǎng)格搜索方法更為流行[24]。許忠淮等[25]最早提出了求P波震源機(jī)制解的格點(diǎn)嘗試法,得到解區(qū)(可選解的集合)及平均解。Reasenberg和Oppenheimer[26]編寫(xiě)了一套求解P波初動(dòng)震源機(jī)制解的程序FPFIT,求最佳解及解的置信區(qū)間,并對(duì)解的質(zhì)量進(jìn)行評(píng)價(jià)。Hardebeek和Shearer[27]改進(jìn)了程序,并討論了P波初動(dòng)極性、震源位置以及速度模型的不確定性對(duì)震源機(jī)制解的影響。俞春泉等[28]對(duì)許忠淮的格點(diǎn)嘗試法進(jìn)行了改進(jìn),提出了數(shù)據(jù)點(diǎn)依賴(lài)于質(zhì)量以及密集程度的權(quán)重計(jì)算方法;提出了通過(guò)聚類(lèi)分析給出單個(gè)或多個(gè)可能平均解的方法;借鑒了杰克刀技術(shù)對(duì)解穩(wěn)定性的測(cè)試;改進(jìn)了解的搜索方式;并提出了評(píng)價(jià)震源機(jī)制解質(zhì)量的新方案。

據(jù)我國(guó)地震臺(tái)網(wǎng)測(cè)定,2009年1月至12月中國(guó)大陸共發(fā)生MS4.5級(jí)以上中強(qiáng)地震28次,最大地震為2009年8月28日青海省西蒙自治州MS6.4級(jí)地震。汶川震區(qū)發(fā)生MS4.5級(jí)以上余震5次,最大為2009年6月30日綿竹MS5.6級(jí)地震。為了給廣大的地震工作者及時(shí)提供中國(guó)大陸中強(qiáng)地震的震源機(jī)制結(jié)果,本文作者在中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司《數(shù)字地震觀測(cè)臺(tái)網(wǎng)產(chǎn)出強(qiáng)化工作》專(zhuān)項(xiàng)的大力支持下,選用國(guó)家地震臺(tái)網(wǎng)和區(qū)域臺(tái)網(wǎng)(震中距500km之內(nèi))的初至波P波初動(dòng)符號(hào)資料,采用下半球等面積投影,用格點(diǎn)嘗試法計(jì)算了2009年中國(guó)大陸MS4.5級(jí)以上地震的震源機(jī)制解共27個(gè),其中MS5.0級(jí)以上地震23次,MS4.5級(jí)以上地震4次。

2 方法及資料

許忠淮等[25]提出了利用P波初動(dòng)方向資料求解震源機(jī)制解的格點(diǎn)嘗試法,其主要技術(shù)思路是:在整個(gè)解的可能空間按一定的步長(zhǎng)進(jìn)行全面搜索嘗試,找出矛盾符號(hào)比φ在一定范圍內(nèi)的所有可能解,計(jì)算其平均解作為震源機(jī)制解的結(jié)果。郭祥云等[29]運(yùn)用本方法計(jì)算汶川地震的震源機(jī)制結(jié)果,并與USGS和哈佛大學(xué)給出的結(jié)果做了一定的對(duì)比,結(jié)果表明該方法比較可靠。

由初動(dòng)符號(hào)求解地震的震源機(jī)制要求地震被盡量多的臺(tái)站記錄到,且臺(tái)站的分布盡量均勻、合理。在本研究中,震源位置參數(shù)采用中國(guó)地震臺(tái)網(wǎng)中心(http:∥www.csndmc.ac.cn/newweb/index.jsp)提供的地震目錄,選用國(guó)家地震臺(tái)網(wǎng)和區(qū)域臺(tái)網(wǎng)(震中距500km之內(nèi))的臺(tái)站記錄。隨著我國(guó)數(shù)字地震臺(tái)網(wǎng)的不斷完善,尤其是“十五”數(shù)字地震觀測(cè)網(wǎng)絡(luò)項(xiàng)目完成之后,提供更好的地震觀測(cè)資料已成為可能。研究所用的地震事件和臺(tái)站分別如圖1和圖2所示。

依據(jù)讀取的P波初動(dòng)數(shù)據(jù),采用格點(diǎn)嘗試法對(duì)每個(gè)地震所有可能的震源機(jī)制解進(jìn)行了全空間掃描計(jì)算,掃描步長(zhǎng)為5°×5°×5°,從中選取矛盾符號(hào)比φ介于最小值φmin至φmin+2%之間的解作為該地震震源機(jī)制解的離散區(qū)(即可能解的分布范圍)。對(duì)離散區(qū)內(nèi)所有可能解平均后作為求解結(jié)果。2009年MS4.5級(jí)以上27次地震,每個(gè)地震的初動(dòng)符號(hào)都在20個(gè)以上,2/3地震的初動(dòng)符號(hào)在30個(gè)以上。所得到的震源機(jī)制解結(jié)果的矛盾符號(hào)比最大為20%,矛盾比≤10%的占88%。所有MS5.0級(jí)以上地震的矛盾比最大為12%,矛盾比≤10%的占95%。每個(gè)地震的震源機(jī)制解、所用初動(dòng)符號(hào)臺(tái)站數(shù)以及矛盾符號(hào)比見(jiàn)表1。

圖1 2009年中國(guó)大陸4.5級(jí)以上地震震中空間分布圖

圖2 求解震源機(jī)制解所用到的臺(tái)站空間分布圖

3 震源機(jī)制解分析

圖3 2009年中國(guó)大陸MS4.5級(jí)以上地震的震源機(jī)制解分布圖(圖中序號(hào)與表1一致)

本研究得到了2009年MS4.5級(jí)以上27次地震的震源機(jī)制解,空間分布見(jiàn)圖3。根據(jù)表2中Zoback[1]研究全球應(yīng)力場(chǎng)時(shí)的分類(lèi)標(biāo)準(zhǔn),對(duì)這些MS4.5級(jí)以上地震的斷層面解類(lèi)型進(jìn)行了分類(lèi)(表1)。結(jié)果表明2009年發(fā)生在中國(guó)大陸的MS4.5級(jí)以上地震以走滑型地震(SS)為主的共13次,逆斷層為主兼走滑的(TS)8次,逆斷層(TF)2次,正斷層為主兼走滑的(NS)3次,正斷層(NF)1次(表2)。這與謝富仁等[30]得到中國(guó)及鄰區(qū)的震源機(jī)制資料50.3%為走滑型,32%為逆斷型,13.9%為正斷型基本一致。其中新疆地區(qū)及青藏塊體的震源機(jī)制解類(lèi)型以走滑(SS)和逆斷層為主兼走滑(TS)的為主。這與徐紀(jì)人等[31]取得青藏高原周緣地區(qū),南部的喜馬拉雅山前沿以及青藏高原的北部、西部邊緣地區(qū)所發(fā)生的地震絕大多數(shù)都屬于逆斷層型或走滑逆斷層型地震結(jié)果相一致。

新疆天山地區(qū)活動(dòng)斷層非常發(fā)育,主要有近EW向、NW-NWW向和NE向3組。新疆的6次地震,主要發(fā)生在3個(gè)地震帶上:北天山地震帶、南天山地震帶、帕米爾-西昆侖地震帶。它們的震源機(jī)制解的P軸方位優(yōu)勢(shì)分布以近SN向?yàn)橹?傾角較小,T軸傾角較大,顯示出新疆天山地區(qū)主要受到近SN向水平擠壓的作用。這與新疆近SN向的區(qū)域構(gòu)造應(yīng)力場(chǎng)方向基本一致[9]。

表2 Zoback[1]研究全球應(yīng)力場(chǎng)時(shí)的分類(lèi)標(biāo)準(zhǔn)

汶川震區(qū)發(fā)生5次地震,位于龍門(mén)山斷裂帶南段的地震4次,其P軸的方位以近NW-SE向?yàn)橹?龍門(mén)山斷裂帶北段的地震1次,P軸方位為NWW-SEE向。龍門(mén)山斷裂帶南段地震的P軸方位與2008年汶川8.0級(jí)地震的P軸的方位和當(dāng)?shù)氐膮^(qū)域應(yīng)力場(chǎng)方位近乎一致,而北段地震的P軸方位則明顯與之前發(fā)生的當(dāng)?shù)氐挠嗾鸬腜軸方位的優(yōu)勢(shì)分布一致。因此,汶川震區(qū)發(fā)生5次地震仍然屬于汶川8.0級(jí)地震的余震。

云南地區(qū)3次地震均為走滑型地震,其P軸方位近SN向,傾角較小,T軸的傾角也較小,顯示出云南地區(qū)主要受近SN向近水平的壓力作用。節(jié)面的滑動(dòng)角接近0°和-180°,表明地震以走滑錯(cuò)動(dòng)為主。

青海地區(qū)6次MS5.0級(jí)以上地震均發(fā)生在走向NWW-SEE、傾向NE的大柴旦-宗務(wù)隆山現(xiàn)代活動(dòng)斷裂帶上。最大地震為2009年8月28日1時(shí)52分青海省西蒙自治州MS6.4級(jí)地震,其節(jié)面Ⅰ的走向、傾角和滑動(dòng)角分別為152°、82°、119°,節(jié)面Ⅱ的走向、傾角和滑動(dòng)角分別為256°、30°、16°。其中節(jié)面Ⅱ的走向與大柴旦-宗務(wù)隆山斷裂帶的走向一致,節(jié)面Ⅱ?yàn)榘l(fā)震斷層。主壓應(yīng)力軸的方位角為SW向,此次地震是受近NE-SW擠壓應(yīng)力場(chǎng)所致。P軸方位角優(yōu)勢(shì)分布(圖4),為近NE-SW方向,與青藏高原整體擠壓應(yīng)力場(chǎng)的方向大體相同。

東北地區(qū)3個(gè)MS4.5級(jí)以上地震,其P軸方位近EW向,相對(duì)于東北地區(qū)的NEESWW方向的主壓應(yīng)力背景應(yīng)力場(chǎng)方向,發(fā)生了一定的偏轉(zhuǎn)。

4 對(duì)比不同計(jì)算方法的震源機(jī)制解

圖4 2009年青海大柴旦地震序列P軸方位角隨時(shí)間的變化

為進(jìn)一步考察本研究結(jié)果的可靠性,將本研究結(jié)果與目前國(guó)內(nèi)外使用不同方法、不同資料來(lái)源中強(qiáng)以上地震的常規(guī)計(jì)算結(jié)果進(jìn)行對(duì)比分析,其中包括全球矩心矩張量解(GCMT,http:∥www.globalcmt.org/CMT-search.html),中國(guó)地震局地球物理研究所提供的采用體波面波聯(lián)合反演矩張量的“剪切-粘貼(CAP)法結(jié)果,以及P波波形快速矩張量反演結(jié)果(本研究中用PCMT表示)。由于P波初動(dòng)法得到的震源機(jī)制解,反映的是地震的初始破裂方式,矩張量解的最佳雙力偶解反映的是整個(gè)地震破裂過(guò)程的平均結(jié)果。一般來(lái)說(shuō),地震過(guò)程的初始破裂方式與平均破裂方式總體性質(zhì)是一致的,但也有不完全相同的時(shí)候。為了對(duì)不同方法得到的震源機(jī)制解結(jié)果進(jìn)行對(duì)比分析,選取了上述4種方法得到的2009年MS5.6級(jí)以上地震的4次地震的震源機(jī)制解的結(jié)果(表3)。

由表3中可見(jiàn),本文所得到的結(jié)果,對(duì)于斷層類(lèi)型來(lái)說(shuō),與GCMT和PCMT的結(jié)果一致,對(duì)于P軸方位角來(lái)說(shuō),本文得到的結(jié)果似乎與GCMT的結(jié)果更接近。

表3 用不同方法得到的中國(guó)大陸地區(qū)2009年MS5.6級(jí)以上地震的震源機(jī)制解

5 結(jié)論與分析

本文采用P波初動(dòng)格點(diǎn)嘗試法得到了2009年中國(guó)大陸MS4.5級(jí)以上地震的震源機(jī)制解,并將結(jié)果與國(guó)內(nèi)外不同計(jì)算方法的結(jié)果進(jìn)行了比較。結(jié)果顯示:

(1)2009年中國(guó)大陸地震的震源機(jī)制解以走滑型地震(SS)和逆斷層為主兼走滑型地震(TS)為主。主壓應(yīng)力軸方位優(yōu)勢(shì)分布具有明顯的區(qū)域特征。表明中國(guó)構(gòu)造應(yīng)力場(chǎng)的格局明顯受制于周邊板塊的動(dòng)力學(xué)作用。我國(guó)新疆及云南地區(qū)主壓應(yīng)力方向主體表現(xiàn)為近SN向;青海省所在區(qū)域主壓應(yīng)力方位主體表現(xiàn)為近NE-SW;中國(guó)東北部主壓應(yīng)力方位主體表現(xiàn)為近EW向,與主壓應(yīng)力方向以NEE-SWW方向?yàn)橹鲗?dǎo)的背景構(gòu)造場(chǎng)略有偏轉(zhuǎn)。2008年發(fā)生8.0級(jí)地震的汶川震區(qū),汶川龍門(mén)山南段地震的P軸方位與汶川8.0級(jí)地震的主應(yīng)力軸的方位和該區(qū)應(yīng)力場(chǎng)方向近乎一致;而龍門(mén)山北段地震的主壓應(yīng)力軸方位相對(duì)于南段的主壓應(yīng)力軸向西發(fā)生了明顯偏轉(zhuǎn)。

(2)國(guó)內(nèi)外不同計(jì)算方法得到的2009年中國(guó)大陸4個(gè)MS5.6級(jí)以上地震的震源機(jī)制結(jié)果有一定的差異。本文得到的結(jié)果與GCMT和PCMT的結(jié)果一致性較好,而與CAP方法得到的震源機(jī)制解有一定的差異。

致謝

感謝中國(guó)地震局地球物理研究所“國(guó)家地震臺(tái)網(wǎng)數(shù)據(jù)備份中心”為本研究提供地震波形數(shù)據(jù)。

(作者電子信箱,郭祥云:gxy-1@sohu.com)

[1]ZobackML.First-andsecond-orderpatternsofstressinthelithosphere:Theworldstressmapproject.GeophysicsRes,1992,97(B8):11703-11728

[2]王曰風(fēng),刁桂苓,張秀萍,等.2001年云南永勝6.0級(jí)地震余震序列震源機(jī)制解與震源應(yīng)力場(chǎng)分析.地震研究,2008,31(2):119-122

[3]許向彤,許忠淮,張東寧,等.求震源機(jī)制解P波初動(dòng)解的格點(diǎn)嘗試概率法.地震地磁觀測(cè)與研究,1995,16(4):34-42

[4]許忠淮,汪素云,高阿甲,等.我國(guó)部分早期震源機(jī)制解的重新測(cè)定.地震地磁觀測(cè)與研究,1994,15(5):1-9

[5]程萬(wàn)正,阮祥,張永久.川滇次極地塊震源機(jī)制解類(lèi)型與一致性參數(shù).地震學(xué)報(bào),2006,28(6):561-573

[6]胡幸平,俞春泉,陶開(kāi),等.利用P波初動(dòng)資料求解汶川地震及其余震震源機(jī)制解.地球物理學(xué)報(bào),2008,51(6):1711-1718

[7]焦青,楊選輝,許麗卿,等.汶川8.0級(jí)地震前后龍門(mén)山斷裂活動(dòng)特征淺析.大地測(cè)量與地球動(dòng)力學(xué),2008,28(4):7-11

[8]劉超,張勇,許力生,等.一種矩張量反演新方法及其對(duì)2008年汶川MS8.0地震序列的應(yīng)用.地震學(xué)報(bào),2008,30(4):329-339

[9]龍海英,聶曉紅,高國(guó)英.新疆天山地區(qū)與川滇地區(qū)中強(qiáng)地震震源機(jī)制解對(duì)比分析.地震研究,2007,30(3):229-234

[10]屠泓為,都昌庭,馬震,等.利用P波初動(dòng)計(jì)算青海中強(qiáng)地震的震源機(jī)制解.高原地震,2008,20(1):22-26

[11]許力生,蔣長(zhǎng)勝,陳運(yùn)泰,等.2004年首都圈地區(qū)中小地震的矩張量反演.地震學(xué)報(bào),2007,29(3):229-239

[12]許忠淮.地震參數(shù)—數(shù)字地震學(xué)在地震預(yù)測(cè)中的應(yīng)用.中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司.北京:地震出版社,2003:16-25

[13]楊中書(shū),李超.利用P波初動(dòng)及振幅比初步計(jì)算2005年11月26日九江—瑞昌地震余震序列震源機(jī)制.華南地震,2006,26(3):77-82

[14]周翠英,華愛(ài)軍,蔣海昆,等.以格點(diǎn)嘗試法求取山東地區(qū)現(xiàn)代中小地震震源機(jī)制解.東北地震研究,2003,19(1):1-11

[15]張萍,谷光裕,高艷玲.岫巖-海城MS5.4級(jí)地震序列震源機(jī)制解.地震,2001,21(1):98-102

[16]張萍,蔣秀琴.岫巖-海城MS5.4級(jí)地震序列的震源機(jī)制解及應(yīng)力場(chǎng)特征.地震地磁觀測(cè)與研究,2001,22(2):76-82

[17]張萍,于龍偉,李涯,等.岫巖-海城MS5.4級(jí)地震前小震震源機(jī)制解與記錄特征分析.地震地磁觀測(cè)與研究,2003,24(1):29-38.

[18]周云好,許力生,陳運(yùn)泰.2000年6月4日印度尼西亞蘇門(mén)答臘南部MS8.0地震的震源機(jī)制.地震學(xué)報(bào),2002,24(5):462-469

[19]Banghar A R.Focalmechanism of Indian earthquakes.Bull.Seismol.Soc.Amer.,1972,62(2):603-608

[20]Crosson R S,Eedo E T.Focalmechanisms of earthquakes related to the29November1975Kalapana,Hawaii,earthquake:The effect of structure models.Bull.Seismol.Soc.Amer.,1981,71(3):713-729

[21]Wu Yih-Min,Zhao Li,Chang Chien-Hsin,et al.Focal-mechanism determination in Taiwan by genetic algorithm.Bull.Seismol.Soc.Amer.,2008,98(2):651-661

[22]Kasahara K.Computer program for afault-p lane solution.Bull.Seismol.Soc.Amer.,1963,53(1):1-13

[23]Brillinger DR,Udias A,Bolt BA.Aprobability model for regional focalmechanism solutions.Bull.Seismol.Soc.Amer.,1980,70(1):149-170

[24]Udias A,Buforn E.Single and joint fault-plane solutions from first-motion data.In:Doornbos DJ ed.Seismological Algorithms.Academic Press,London,1988:443-453

[25]許忠淮,閻明,趙仲和.由多個(gè)小地震推斷的華北地區(qū)構(gòu)造應(yīng)力場(chǎng)的方向.地震學(xué)報(bào),1983,5(3):268-279

[26]ReasenbergP,OppenheimerD.FPFIT,FPPLOT,and FPPAGE:FORTRANcomputer programs for calculating and displaying earthquake fault-plane solutions.U.S.Geo1.Surv.Open-File Rept,85-739,109

[27]Hardebeek J L,Shearer P M.A new method for determining first-motion focal mechanisms.Bull.Seismol.Soc.Amer.,2002,92(6):2264-2276

[28]俞春泉,陶開(kāi),崔效鋒,等.用格點(diǎn)嘗試法求解P波初動(dòng)震源機(jī)制解及解的質(zhì)量評(píng)價(jià).地球物理學(xué)報(bào),2009,52(5):1402-1411

[29]郭祥云,陳學(xué)忠,李艷娥.2008年5月12日四川汶川8.0級(jí)地震與部分余震的震源機(jī)制解.地震,2010,30(1):50-60

[30]謝富仁,崔效鋒,趙建濤.全球應(yīng)力場(chǎng)與構(gòu)造分析.地學(xué)前緣,2003,10(8):22-30

[31]徐紀(jì)人,趙志新.青藏高原及其周?chē)貐^(qū)區(qū)域應(yīng)力場(chǎng)與構(gòu)造運(yùn)動(dòng)特征.中國(guó)地質(zhì),2006,33(2):275-285

Focal Mechanism Solutions for Some Earthquakes in the Chinese Mainland in2009

Guo Xiangyun and Chen Xuezhong
(Institute of Geophysics,CEA,Beijing100081,China)

In this paper,the focal mechanism solutions of some earthquakes of MS≥4.5occurring in the Chinese mainland in2009are determined with the lower hemisphere equal-area projection and first motion direction data of PnandPgwaves recorded by the regional and distant seismic stations.In the calculation,the grid point test method is employed.The types of the focal fault plane are classified according to the classificatory criteria proposed by Zoback(1992)through the study of the global stress field.The results show that focal mechanism solutions of13earthquakes show strike-slip(SS)faulting,8are predominately thrust-faulting with strike-slip component(TS),2thrust faulting(TF),3predominately normal faulting with strike-slip component(NS),and only 1 normal faulting(NF).The predominant orientation of principal compressive stressP axis is different for different regions,and coincident with the local tectonic stress field,but the predom inant orientation of P axis in the northeast of China deflected from the local tectonic stress field.

firstmotion direction;the Chinese mainland;focal mechanism

P315.3+3;

A;

10.3969/j.issn.0235-4975.2010.02.006

2009-12-29;

2010-01-13。

中國(guó)地震局一般項(xiàng)目(項(xiàng)目編號(hào):049214)資助。