彭利媚，魏婭玲

（1.海南省地震局，?？?72900；2.四川省地震局，成都610041）

由ML≥4.0地震震源機制解推斷四川理縣-黑水斷裂構(gòu)造特征

彭利媚1，魏婭玲2

（1.海南省地震局，?？?72900；2.四川省地震局，成都610041）

基于四川臺網(wǎng)寬頻帶數(shù)字地震波形資料和CAP波形反演方法，提取了四川理縣-黑水地區(qū)2008年1月至2014年12月10個ML≥4.0中等強度地震的震源機制解。結(jié)果顯示：10個中等強度地震震源機制解節(jié)面分布主要集中在北西-南東和北東-南西兩個優(yōu)勢分布方向。據(jù)眾多學者研究結(jié)果推斷，北西-南東向展布的節(jié)面為發(fā)震斷層面；理縣-黑水斷裂大致沿北47°西（近北西向）展布，屬高傾角陡峭斷裂，主傾向南西；壓應力（P軸）呈近東西水平擠壓特征，與斷層面斜交，以左旋走滑運動為主；張應力（T軸）呈近南北水平撕裂特征；中等強度地震的震源主壓應力方向與青藏地塊東部地區(qū)主壓應變率方向比較一致，說明理縣-黑水斷裂主要是受區(qū)域應力場控制的一個區(qū)域性構(gòu)造左旋走滑斷裂。

發(fā)震斷層；震源機制解；CAP方法；波形互相關(guān)系數(shù)；斷層類型

0 引言

地震震源機制解直觀地反映了地震破裂的幾何特征和運動學特征，是研究區(qū)域構(gòu)造應力的基礎(chǔ)，因此震源機制解一直是地震工作者們關(guān)心的重要參數(shù)[1]。區(qū)域性寬頻帶地震波形資料中包含了豐富的地震震源和區(qū)域地殼信息，利用寬頻帶波形資料研究區(qū)域性地震震源機制解有助于了解震源區(qū)應力狀態(tài)和發(fā)震斷層構(gòu)造特征，也是理解地震孕育過程的重要途徑。

表1 四川理縣-黑水ML≥4.0地震目錄Table 1 The catalog of the earthquakes with ML≥4.0 in Lixian-Heishui

2008年汶川8.0級地震發(fā)生后，國內(nèi)眾多學者研究了汶川地震余震的精確定位、空間分布、震源機制解等[2-8]，認為余震序列主要分布在龍門山斷裂帶的中、北段，沿三條主斷裂帶展布；同時，小魚洞-理縣-黑水的北西向余震分支也是目前各個研究結(jié)果中普遍存在的現(xiàn)象。據(jù)地質(zhì)調(diào)查[9-10]已證實小魚洞斷裂的存在，該斷裂南東起于彭州市通濟場西南側(cè)附近，北西止于彭州市小魚洞鎮(zhèn)后壩村附近，大體呈北西-南東向展布。繼而向西的理縣-黑水地震分支條帶的存在[4，7-8]也預示著有一條沿該地震帶延伸的隱伏斷裂，我們將其稱為理縣-黑水斷裂；而汶川地震前、后，研究區(qū)域并沒有發(fā)現(xiàn)出露地表的斷裂帶（或地震帶）存在，因此很有必要對其分支條帶上的中等強度地震震源機制解進行研究分析。為了深入了解這一斷裂帶地震震源的發(fā)震機制，本文選用目前使用較為廣泛的CAP[11-12]方法和四川臺網(wǎng)寬頻帶地震波形資料，反演計算了四川理縣-黑水地區(qū)ML≥4.0中等強度地震震源機制解和最佳震源深度，并根據(jù)地震震源機制解的空間分布，對理縣-黑水斷裂帶的發(fā)震斷層面展布特征進行了初步分析。

1 數(shù)據(jù)和模型

1.1 臺站和資料的選取

據(jù)中國地震臺網(wǎng)中心統(tǒng)一地震目錄統(tǒng)計，2008年5月至2014年12月四川理縣-黑水地區(qū)共發(fā)生ML≥4.0地震10個（見表1）。理縣地區(qū)發(fā)生8次，最大震級為MS5.9；黑水地區(qū)發(fā)生2次，最大震級為ML4.8；這些中等強度地震為研究該地區(qū)地震活動構(gòu)造和應力場提供了重要信息。這些地震發(fā)生前,震中區(qū)域附近250 km范圍內(nèi)正常運行的固定臺站有30個，2010年災后重建8個，在0.05～20 Hz的頻率范圍內(nèi)，這些臺站儀器的幅頻特性曲線都是平直的，獲得了較好的寬頻帶儀器記錄波形，為研究地震震源性質(zhì)提供了重要的資料。經(jīng)過臺站選擇，最終選用了40 km≤Δ≤250 km、P波初動明顯、波形連續(xù)且信噪比相對較高的臺站進行波形反演計算。地震震中與所選臺站分布如圖1所示，從圖中可以看出，所選臺站對研究區(qū)域地震震中具有較強的控制能力。

1.2 方法和速度模型

圖1 震中及參加反演的臺站分布Fig.1 The distribution of epicenters and inversion stations

本文利用CAP波形反演方法求取理縣-黑水ML≥4.0地震震源機制及震源深度，具體方法參見相關(guān)文獻[5，11-14]。該方法將寬頻帶數(shù)字地震波形記錄分為體波Pnl與面波兩部分，分別對Pnl波、面波進行帶通濾波，并計算其理論地震波形與觀測波形之間的誤差函數(shù)，通過網(wǎng)格搜索，獲取給定參數(shù)空間中誤差函數(shù)達到最小的最佳解。本文計算時，體波Pnl與面波截取波形窗長分別設(shè)置為30 s與60 s，其帶通濾波頻帶寬度分別為0.05～0.2 Hz、0.05～0.1 Hz。與其他求解震源機制解的方法相比，CAP方法彌補了利用P波初動求解震源機制解受臺站數(shù)量限制的缺陷，其優(yōu)點體現(xiàn)在反演震源機制解時使用的是近震寬頻帶波形資料，反演結(jié)果對地殼橫向變化不敏感、對速度模型依賴性相對較小等，保證了反演結(jié)果的穩(wěn)定性。本研究中，我們選用了林向東等的速度模型[15]進行地震波反演，這一模型是根據(jù)川西地區(qū)人工地震測深剖面和波形擬合結(jié)果，并做適當調(diào)整后建立的速度模型（見圖2）。

圖2 地震波反演的速度模型Fig.2 The velocity models inverted by seismic waveform

2 反演結(jié)果及分析

2.1 震源機制解

我們采用CAP方法和圖2速度模型進行地震波反演計算，獲得了研究區(qū)域內(nèi)的10個ML≥4.0地震震源機制解和最佳震源深度。為了保證反演結(jié)果的穩(wěn)定，我們盡量保證參加反演計算的臺站四方位包圍震中、且各臺站的波形擬合互相關(guān)系數(shù)要大于60%。通過反演計算獲得的地震波形擬合圖和最佳深度誤差擬合圖得知，參加反演計算的10個ML≥4.0地震的波形擬合效果均比較理想，震源機制解結(jié)果是穩(wěn)定、可靠的，表2列出了地震震源機制解詳細參數(shù)和最佳震源深度值，地震序號與表1的相同。

圖3（a）和圖4分別是2008年5月16日理縣MS5.9地震震源深度擬合誤差和地震波形擬合圖。從圖3（a）可以看出，震源機制解反演結(jié)果在2～20 km震源深度范圍內(nèi)均穩(wěn)定，震源深度在13 km處獲得最小深度擬合誤差；以13 km為最佳震源深度計算出的理論波形均能較好地與實際記錄波形匹配（圖中紅線表示理論波形，黑線表示實際記錄波形），10個臺站的波形擬合互相關(guān)系數(shù)均大于62%，波形擬合效果理想，且結(jié)果穩(wěn)定、可靠（見圖4）。其中，節(jié)面I的解：走向為147°，傾角為81°，滑動角為-31°，節(jié)面Ⅱ的解：走向為242°，傾角為59°，滑動角為-170°，矩震級為MW=5.4。該地震震源機制解反演結(jié)果與鄭勇等[5]的研究結(jié)果一致。

表2 四川理縣-黑水ML≥4.0地震震源機制反演結(jié)果Table 2 The inversion results of the earthquakes with ML≥4.0 in Lixian-Heishui

圖3 地震震源深度擬合誤差Fig.3 The fitting error of earthquake focal depth

圖4 理縣MS5.9地震波形擬合圖Fig.4 The fitting figure of seismic waveform of the MS5.9 earthquake in Lixian

圖5黑水ML4.8地震波擬形合圖Fig.5 The fitting figure of seismic waveform of the ML4.8 earthquake in Heishui

圖3 （b）和圖5分別是我們得到的2008年11月14日黑水ML4.8地震震源深度擬合誤差和地震波形擬合圖。從圖3（b）可以看出，震源機制反演結(jié)果在4～20 km震源深度范圍內(nèi)均穩(wěn)定，震源深度在11 km處獲得最小深度擬合誤差；以11 km為最佳震源深度計算出的理論波形均能較好地與實際記錄波形匹配，8個臺站的波形擬合互相關(guān)系數(shù)均大于60%，波形擬合效果理想，且結(jié)果穩(wěn)定、可靠（見圖5）。其中，節(jié)面I的解為：走向311°，傾角72°，滑動角7°，節(jié)面Ⅱ的解為：走向219°，傾角83°，滑動角162°，矩震級MW=4.6。

2.2 結(jié)果分析

在地震破裂面未達到地表的情況下，地震震中分布是判斷斷裂帶分布的主要依據(jù)。上述國內(nèi)眾多學者研究結(jié)果推斷，理縣-黑水斷裂斷層面是沿北西向展布的；且從圖1和圖6也可以看出其斷裂帶的大致分布方向，圖6顯出的地震震源機制解節(jié)面分布主要集中在北西-南東和北東-南西兩個優(yōu)勢分布方向。由此可知，表2中10個地震震源機制解的節(jié)面I應為發(fā)震斷層面，其走向空間分布見圖7（a），斷層面走向的優(yōu)勢分布方向為南東，主要傾向南西；節(jié)面II應為輔助面,其走向空間分布見圖7（b），輔助面走向的優(yōu)勢分布方向為南西西。根據(jù)Zoback等的斷層分類原則[16]對表2中地震震源機制解結(jié)果進行分類，理縣-黑水10個ML≥4.0地震的震源機制解類型均屬于走滑斷層。

圖6 理縣-黑水ML≥4.0地震震源機制解空間分布Fig.6 The spatial distribution of the focal mechanism solution of the earthquakes with ML≥4.0 in Lixian-Heishui

根據(jù)發(fā)震斷層面機制解結(jié)果及其空間分布進行分析，總體上說，斷層面分布范圍沿北33°～58°西（平均約北47°西，近北西向）展布；傾角范圍為61°～82°（平均約75°），屬高傾角陡峭斷裂，傾向南西，南西盤為上盤，北東盤為下盤；滑動角范圍為-31°～10°（平均約-7°），呈現(xiàn)出一定的正傾滑分量，西南盤為下降盤；壓應力（P軸）方位為80°～100°（平均約89°），呈現(xiàn)近東西水平擠壓特征，與斷層面斜交（見圖8），其仰角范圍為5°～28°（平均約為17°），水平滑動分量大于垂直滑動分量，以左旋走滑運動為主；張應力（T軸）方位平均約為北2°西，呈現(xiàn)近南北水平撕裂特征。因此說明理縣-黑水斷裂帶發(fā)震斷層面大致沿近北西方向展布，以高角度的左旋走滑運動為主，略具一定正傾滑分量作用的性質(zhì)。

根據(jù)上述研究結(jié)果和地質(zhì)調(diào)查[9-10]得知，理縣-黑水斷裂位于小魚洞-理縣-黑水北西向地震條帶的北西段，小魚洞斷裂位于小魚洞-理縣-黑水北西向地震條帶的東南段，且兩者的斷層面產(chǎn)狀較為一致，可以推斷出兩者為同一斷裂。另外，理縣-黑水10個ML≥4.0地震震源主壓應力呈近東西擠壓的特征，與武艷強等[17]用GPS數(shù)據(jù)得到的近年來青藏地塊東部地區(qū)主壓應變率為近東西向的結(jié)論一致，從而說明小魚洞-理縣-黑水斷裂主要是受區(qū)域應力場控制的一個區(qū)域性構(gòu)造左旋走滑斷裂。

3 結(jié)語

本文基于四川臺網(wǎng)寬頻帶數(shù)字地震波形資料和CAP波形反演方法，求取了四川理縣-黑水地區(qū)2008年1月至2014年12月10個ML≥4.0地震的震源機制解，分析了震源機制解空間分布特征，并對發(fā)震斷層的幾何產(chǎn)狀及斷裂與區(qū)域應力場的關(guān)系進行了初步探討。綜上所述，得出以下幾點結(jié)論：

圖7 理縣-黑水ML≥4.0地震震源機制解節(jié)面及其P軸空間分布Fig.7 The nodal plane and P axis distribution of the focal mechanism solutions of the earthquakes with ML≥4.0 in Lixian-Heishui

圖8 理縣-黑水ML≥4.0地震震源機制解應力空間分布（左圖P軸，右圖T軸）Fig.8 The stress distribution of the focal mechanism solutions of the earthquakes with ML≥4.0 in Lixian-Heishui

（1）10個中等強度地震震源機制解節(jié)面分布主要集中在北西-南東和北東-南西兩個優(yōu)勢分布方向。據(jù)眾多學者研究結(jié)果推斷，北西-南東向展布的節(jié)面為發(fā)震斷層面。

（2）理縣-黑水斷裂大致沿北47°西（近北西向）展布，屬高傾角陡峭斷裂，主傾向南西。

（3）壓應力（P軸）呈近東西水平擠壓特征，與斷層面斜交，以左旋走滑運動為主；張應力（T軸）呈近南北水平撕裂特征。

（4）中等強度地震的震源主壓應力方向與青藏地塊東部地區(qū)主壓應變率方向比較一致，說明理縣-黑水斷裂主要是受區(qū)域應力場控制的一個區(qū)域性構(gòu)造左旋走滑斷裂。

致謝：四川省地震局魏婭玲高級工程師給予大力的幫助和指導。

[1]鄧起東，張培震，冉勇康，等．中國活動構(gòu)造基本特征[J].中國科學D輯：地球科學，2002，32（12）：1 020-1 030.

[2]呂堅，蘇金蓉，靳玉科，等，汶川MS8.0級地震余震序列重新定位及其發(fā)震構(gòu)造初探[J].地震地質(zhì)，2008，30（4）：917-925.

[3]王勤彩，陳章立，鄭斯華.汶川大地震余震序列震源機制的空間分段特征[J].科學通報，2009，54（16）：2 348-2 354.

[4]陳九輝，劉啟元，李順成，等，汶川8.0地震余震序列重新定位及其地震構(gòu)造研究[J].地球物理學報，2009，52（2）:390-397.

[5]鄭勇，馬宏生，呂堅，等.汶川地震強余震（MS≥5.6）的震源機制解及其與發(fā)震構(gòu)造的關(guān)系.中國科學D輯：地球科學，2009，39（4）：413-426.

[6]張勇，許力生，陳運泰．2008年汶川大地震震源機制的時空變化．地球物理學報，2009，52（2）：379-389.

[7]劉巧霞，朱介壽，曹俊興，等，汶川MS8.0級地震余震序列重新定位及其空間分段特征[J].第四紀研究，2010，30（4）：736-744.

[8]易桂喜，龍鋒，張致偉．汶川MS8．0地震余震震源機制時空分布特征[J]．地球物理學報，2012，55（4）：1213-1227.

[9]徐錫偉，聞學澤，葉建青，等.汶川MS8.0地震地表破裂帶及其發(fā)震構(gòu)造[J].地質(zhì)科學，2008，30（3）：598-629.

[10]李勇，黃潤秋，A.L.Densmone，等.龍門山小魚洞破裂在汶川地震中的地表破裂及地質(zhì)意義[J].第四紀研究，2009，29（3）：502-511.

[11]Zhao L S，Helmberger D V．Source estimation from broadband regional seismograms[J].Bull Seis Soc Amer，1994，84：91-104.

[12]Zhu L P，Helmberger D V．Advancement in source estimation techniques using broadband regional seismograms [J].Bull Seis Soc Amer，1996，86：1 634-1 64l.

[13]陳長云，賀建明.西秦嶺北緣斷裂分段變形特征分析[J].大地測量與動力學，2016，36（9）：784-788.

[14]陳長云.張家口-渤海斷裂帶分段運動變形特征分析[J].地震，2016，36（1）：1-11.

[15]林向東，葛洪魁，徐平，等.近場全波形反演：蘆山7.0級地震及余震矩張量解[J].地球物理學報，2013，56（12）：4 037-4 047.

[16]Zoback M L.Firest-and second-order patterns of stress in the lithosphere:The world stress map project[J].J.Geophys. Res.，1992，97（B8）：11703-11728.

[17]武艷強，江在森，閆偉，等.中國大陸應力、應變率場方向特征分析[J].大地測量與地球動力學，2012，32（1）：5-9.

Tectonic Characteristics of Lixian-Heishui Fault Inferred by the Focal Mechanism Solution of the Earthquakes with ML≥4.0 in Sichuan

PENG Limei1，WEI Yaling2
（1.Earthquake Administration of Hainan Province，Haikou 572900，China；2.Earthquake Administration of Sichuan Province，Chengdu 610041，China）

Based on broadband digital seismic waveform recorded by Sichuan earthquake network and CAP waveform inversion method,the paper obtained the focal mechanism solutions of the ten moderately strong earthquakes with ML≥4.0 from January 2008 to December 2014 in region of sichuan Lixian-Heishui.The results show that:the nodal plane of the focal mechanism solutions of ten moderately strong earthquakes mainly concentrated in two dominance distribution directions（NW-SE and NE-SW），and according to many scholars'research results inference,NW-SE trending section is earthquake fault plane;Lixian-Heishui fault is mainly distributed along N47°W（near north west），which belongs to high dip angle steep fault with dipping to SW；compressive stress（P axis）presents nearly EW horizontal extrusion features,and inclines with fault plane with takeing the left lateral strike slip movement;tensile stress（T axis）presents nearly NS horizontal tearing features.The orientation of principal compressive stress of moderately strong earthquakes is consistent with the orientation of the principal compressive strain rate of eastern part of the Qinghai-Tibet block,which indicates that Lixian-Heishui fault is mainly controlled by regional stress field and is left lateral strike slip fault.

Seismogenic fault；Focal mechanism solution；CAP method；Waveform correlation coefficient to each otherl；Fault type

P315.33

A

1001-8662（2016）03-0009-08

10.13512/j.hndz.2016.03.002

2015-10-02

中國地震局測震臺網(wǎng)青年骨干培養(yǎng)專項（20140320）

彭利媚(1978-)，女，工程師，從事地震監(jiān)測及編目工作.

E-mail：623294290@qq.com

魏婭玲(1979-)，女，高級工程師，從事大震速報工作.

E-mail：sctwsbs@163.com

彭利媚，魏婭玲.由ML≥4.0地震震源機制解推斷四川理縣-黑水斷裂構(gòu)造特征[J].華南地震，2016，36（3）：9-16.[PENG Limei，WEI Yaling. Tectonic Characteristics of Lixian-Heishui Fault Inferred by the Focal Mechanism Solution of the Earthquakes with ML≥4.0 in Sichuan[J].South china journal of seismology，2016，36（3）：9-16.]