王曉山， 呂堅(jiān)， 謝祖軍， 龍鋒， 趙小艷， 鄭勇

1 河北省地震局， 石家莊 050021 2 江西省地震局， 南昌 330039 3 大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，中國(guó)科學(xué)院測(cè)量與地球物理研究所， 武漢 430077 4 四川省地震局， 成都 610041 5 云南省地震局， 昆明 650224

?

南北地震帶震源機(jī)制解與構(gòu)造應(yīng)力場(chǎng)特征

王曉山1， 呂堅(jiān)2, 3*， 謝祖軍3， 龍鋒4， 趙小艷5， 鄭勇3

1 河北省地震局， 石家莊 050021 2 江西省地震局， 南昌 330039 3 大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，中國(guó)科學(xué)院測(cè)量與地球物理研究所， 武漢 430077 4 四川省地震局， 成都 610041 5 云南省地震局， 昆明 650224

南北地震帶作為中國(guó)大陸地應(yīng)力場(chǎng)一級(jí)分區(qū)的邊界，其構(gòu)造應(yīng)力場(chǎng)的研究對(duì)理解大陸強(qiáng)震機(jī)理、構(gòu)造變形和地震應(yīng)力的相互作用具有重要意義.本文收集南北地震帶1970—2014年的震源機(jī)制解819條，按照全球應(yīng)力圖的分類標(biāo)準(zhǔn)對(duì)震源機(jī)制解進(jìn)行分類，發(fā)現(xiàn)其空間分布特征與地質(zhì)構(gòu)造活動(dòng)性質(zhì)比較吻合.P軸水平投影指示了活動(dòng)塊體的運(yùn)動(dòng)方向，T軸水平投影在川滇塊體及鄰近地區(qū)空間差異特征最為突出，存在順時(shí)針旋轉(zhuǎn)的趨勢(shì).南北地震帶的最大水平主應(yīng)力方向具有明顯的分區(qū)特征，北段為NE向走滑類型的應(yīng)力狀態(tài)，中段為NEE—EW—NWW向的逆沖類型，南段為SE—SSE—NS—NNE向走滑和正斷類型，在川滇塊體的北部和西邊界應(yīng)力狀態(tài)為EW—SE—SSE的正斷層類型，表明來(lái)自印度板塊的NNE或NE向的水平擠壓應(yīng)力和青藏高原物質(zhì)東向滑移沿大型走滑斷裂帶向SE向平移的復(fù)合作用控制了南北地震帶的巖石圈應(yīng)力場(chǎng).川滇塊體西邊界正斷層類型應(yīng)力狀態(tài)范圍與高分辨率地震學(xué)觀測(cè)得到的中下地殼低速帶范圍基本吻合，青藏高原向東擴(kuò)張的塑性物質(zhì)流與橫向邊界(麗江—小金河斷裂帶)的弱化易于應(yīng)變能的釋放，在局部地區(qū)使NS向拉張的正斷層向EW向拉張正斷層轉(zhuǎn)變.反演得到的應(yīng)力狀態(tài)基本上與各種類型地震的破裂方式比較吻合，也進(jìn)一步驗(yàn)證反演結(jié)果的可靠性，可為地球動(dòng)力學(xué)過(guò)程的模擬和活動(dòng)斷層滑動(dòng)性質(zhì)的厘定提供參考.

南北地震帶； 震源機(jī)制解； 構(gòu)造應(yīng)力場(chǎng)； 最大水平主應(yīng)力方向

1 引言

南北地震帶北起鄂爾多斯地塊西緣，跨越秦嶺，穿過(guò)龍門(mén)山，再沿安寧河—小江斷裂帶向南延至緬甸境內(nèi)，成為分隔中國(guó)大陸東部相對(duì)穩(wěn)定的鄂爾多斯高原、四川盆地和華南地塊與西部強(qiáng)烈隆升的青藏高原之間的邊界活動(dòng)構(gòu)造帶(Zhang，2013).深部地球物理探測(cè)表明，南北地震帶本身為非常明顯的地殼厚度梯度帶和重力梯度帶(馬杏垣，1989；Huang et al.，2003).在不同地球動(dòng)力學(xué)過(guò)程共同作用下，晚新生代和現(xiàn)今構(gòu)造變形在南北帶兩側(cè)的差異非常明顯，形成了一條由北向南分段的、不同性質(zhì)斷裂和褶皺構(gòu)成近SN向的復(fù)雜構(gòu)造帶和地震活動(dòng)帶，集中了中國(guó)有歷史記錄以來(lái)一半的8級(jí)以上大地震，是中國(guó)大陸淺源地震最活躍、地震分布帶狀明顯的巨型地震帶，亦是中國(guó)大陸地應(yīng)力場(chǎng)一級(jí)分區(qū)的邊界(謝富仁等，2004).2008年5月12日汶川8.0級(jí)及其后續(xù)一系列6級(jí)以上地震即是南北地震帶構(gòu)造活動(dòng)的直接反映.因此，獲取整個(gè)南北地震帶的構(gòu)造應(yīng)力場(chǎng)，對(duì)探討驅(qū)動(dòng)南北地震帶構(gòu)造變形和強(qiáng)震發(fā)生的動(dòng)力學(xué)機(jī)制具有十分重要的意義.

構(gòu)造應(yīng)力場(chǎng)在地球動(dòng)力學(xué)研究中起著十分重要的作用，區(qū)域地殼應(yīng)力場(chǎng)及其時(shí)空變化特征為深入了解地震機(jī)理、構(gòu)造變形及地震應(yīng)力的相互作用提供了可靠的信息.對(duì)地殼內(nèi)構(gòu)造應(yīng)力進(jìn)行直接測(cè)量非常困難，但可以利用其他的一些地球物理觀測(cè)資料反演出構(gòu)造應(yīng)力的某些特征，例如地震震源機(jī)制和地應(yīng)力測(cè)量資料等.基于這些資料，許多學(xué)者就南北地震帶的構(gòu)造應(yīng)力場(chǎng)進(jìn)行了研究.通過(guò)地震P波初動(dòng)極性(許忠淮等，1987，1989)、震源機(jī)制解(闞榮舉等，1977；鄢家全等，1979；成爾林，1981；秦保燕等，1986；徐紀(jì)人和尾池和夫，1995；杜興信和邵輝成，1999；程萬(wàn)正等，2003；謝富仁等，2004；徐紀(jì)人等，2008；Wan，2010；郝平等，2012；郭祥云等，2014)、斷層擦痕數(shù)據(jù)(謝富仁等，2001)、實(shí)測(cè)地應(yīng)力(楊樹(shù)新等，2012)等進(jìn)行整體和分段的研究，目前對(duì)于南北地震帶的基本一致結(jié)論是該地震帶為中國(guó)地震構(gòu)造應(yīng)力場(chǎng)的一條分界線，應(yīng)力場(chǎng)方向變化很劇烈；北段成北北東向，中段成近東西向，南段成北北西向.同時(shí)，南北地震帶在不同區(qū)域的應(yīng)力大小存在著差異，杜興信和邵輝成(1999)利用FMSI軟件(Gephart and Forsyth，1984)給出了相對(duì)應(yīng)力大小R值，北段R值高于0.60，南段低于0.50，大多數(shù)斷層為走滑斷層.

此外，對(duì)于南北地震帶的構(gòu)造應(yīng)力背景也得到了一定的認(rèn)識(shí).例如，楊樹(shù)新等(2012)得出南北地震帶北段擠壓作用強(qiáng)度僅次于青藏板塊，中段擠壓作用強(qiáng)度弱于北段強(qiáng)于南段.闞榮舉等(1977)指出，川滇菱形塊體由于受到印度洋板塊對(duì)歐亞板塊的擠壓的影響，向東南方向運(yùn)動(dòng)并控制著那里的應(yīng)力場(chǎng).根據(jù)在亞洲東部和東南部設(shè)為自由邊界的條件下的實(shí)驗(yàn)?zāi)M結(jié)果，Peltzer和Tapponnier(1988)認(rèn)為，由北向運(yùn)動(dòng)的印度—澳大利亞板塊所產(chǎn)生的強(qiáng)擠壓力控制著南北地震帶南段的構(gòu)造應(yīng)力場(chǎng).徐紀(jì)人和尾池和夫(1995)發(fā)現(xiàn)在南北地震帶南段P軸的方位在空間組成一個(gè)倒“V”字形，東部和西部地區(qū)的邊界與青藏高原和揚(yáng)子塊體之間的邊界是一致的.王振聲等(1976)分析南北地震帶的范圍及其地震活動(dòng)特征，結(jié)果顯示104°E線為一天然的東西分界線，震中的南北呼應(yīng)關(guān)系在這一線上較為顯著，這些跡象暗示地下深處可能存在一斷續(xù)的南北構(gòu)造線，由此推測(cè)，沿此構(gòu)造線可能存在一南北走向的應(yīng)力集中帶.

以往的這些研究給出了南北地震帶總體的構(gòu)造應(yīng)力場(chǎng)分布，以及一些局部斷層的應(yīng)力狀態(tài)，但對(duì)于南北地震帶整體上的各斷層性質(zhì)，以及其對(duì)應(yīng)的孕震環(huán)境研究方面仍存在一定的不足.由于地震震源機(jī)制解提供了斷層的走向、傾角、滑動(dòng)角以及應(yīng)力釋放的定量信息，因此，通過(guò)對(duì)震源機(jī)制解的分析，將能夠給出整個(gè)區(qū)域的發(fā)震環(huán)境和構(gòu)造應(yīng)力場(chǎng)特征.特別是南北地震帶地震活動(dòng)性很高，從1970年以來(lái)，該區(qū)域發(fā)生了近千次4級(jí)以上的地震，空間分布基本上覆蓋了整個(gè)南北地震帶，加上臺(tái)網(wǎng)布設(shè)較為密集，為研究地震位置和震源機(jī)制解提供了良好的數(shù)據(jù)基礎(chǔ).因此，獲取該區(qū)域近幾十年來(lái)的地震震源位置和震源機(jī)制解，并由此分析南北地震帶的斷層性質(zhì)和構(gòu)造應(yīng)力場(chǎng)特征成為可能.本文將通過(guò)反演和收集1970年以來(lái)的4.0級(jí)以上地震震源機(jī)制解，并綜合Global Centroid Moment Tensor (GCMT) database提供的部分震源機(jī)制解，反演其應(yīng)力張量和水平最大主壓應(yīng)力方向，進(jìn)而給出南北地震帶整體的構(gòu)造應(yīng)力特征.

2 數(shù)據(jù)和方法

2.1 地震數(shù)據(jù)

選取南北地震帶(21°N—43°N，97°E—107°E區(qū)域)作為預(yù)訂研究范圍，使用的數(shù)據(jù)包括：來(lái)自中國(guó)及其鄰區(qū)地殼應(yīng)力場(chǎng)數(shù)據(jù)庫(kù)(謝富仁等，2003)的發(fā)生在1970—1983年的275個(gè)震源機(jī)制解數(shù)據(jù)；下載自GCMT database的1976—2014年的211個(gè)矩心矩張量數(shù)據(jù)(Dziewonski et al.，1981；Ekstr?m et al.，2012)；其他數(shù)據(jù)為筆者等(將另擬文發(fā)表)利用國(guó)家或者區(qū)域數(shù)字化測(cè)震臺(tái)站的波形數(shù)據(jù)(鄭秀芬等，2009)反演的發(fā)生在2007年10月—2013年12月的333個(gè)4級(jí)以上地震震源機(jī)制解數(shù)據(jù).以上數(shù)據(jù)如果出現(xiàn)重合，則以GCMT解數(shù)據(jù)為準(zhǔn)，最終參與反演的數(shù)據(jù)共計(jì)819個(gè)震源機(jī)制解(圖1).

2.2 研究方法

通過(guò)震源機(jī)制解可以計(jì)算地震的應(yīng)力張量.應(yīng)力張量既可以由震源機(jī)制解反演得到，也可以直接用P波初動(dòng)極性得到(Abers and Grphart，2001；萬(wàn)永革等，2011a).大多數(shù)成熟的法向應(yīng)力反演技術(shù)有兩個(gè)主要的假設(shè)：(1)研究區(qū)域地殼范圍內(nèi)應(yīng)力場(chǎng)是均勻的；(2)斷層的滑動(dòng)方向與應(yīng)力張量投影在斷層面上的剪切應(yīng)力方向一致(Wallace，1951；Bott，1959).估計(jì)應(yīng)力場(chǎng)方向是非線性反演問(wèn)題，既可以直接求解，也可以采用線性化求解.在求解非線性反演問(wèn)題時(shí)，使用網(wǎng)格搜索法(Gephart and Forsyth，1984；Gephart，1990；Arnold and Townend，2007)或基于采樣優(yōu)化的蒙特卡羅方法(Angelier，1984；Xu，2004)得到一組與震源機(jī)制達(dá)到最優(yōu)擬合的應(yīng)力張量.線性反演方法是通過(guò)計(jì)算量小、被廣泛應(yīng)用的最小二乘法求得(Michael，1984；Hardebeck and Michael，2006).本研究使用MSATSI軟件包(Martínez-Garzón et al.，2014)，根據(jù)SATSI(Spatial And Temporal Stress Inversion)算法(Hardebeck and Michael，2006)重新設(shè)計(jì)，采用各種圖表示0D—3D的應(yīng)力反演結(jié)果，使其可視化.SATSI算法程序中，輸入的震源機(jī)制解在反演之前被分到按一定維度分布的許多子區(qū)(網(wǎng)格).維度從0D(單個(gè)網(wǎng)格的應(yīng)力反演)到1D(應(yīng)力場(chǎng)隨時(shí)間變化)再上升到4D(時(shí)空分布).4D是把時(shí)間作為一個(gè)額外的坐標(biāo).然后，采用阻尼最小二乘反演(Menke，1989)同時(shí)得到每個(gè)網(wǎng)格點(diǎn)的應(yīng)力張量并進(jìn)行平滑.SATSI算法采用了線性反演技術(shù)，相對(duì)于網(wǎng)格搜索法反演速度更快，能夠最大程度上擬合數(shù)據(jù)，消除人為影響，在反映應(yīng)力場(chǎng)真實(shí)變化的同時(shí)，又去掉了過(guò)度平滑所帶來(lái)的虛假變化.對(duì)于SATSI算法的詳細(xì)描述，請(qǐng)參考Hardebeck和Michael(2006)和Michael(1984，1987).

圖1 南北地震帶震源機(jī)制解類型空間分布

Lund和Townend(2007)從數(shù)學(xué)上推導(dǎo)給出了由應(yīng)力張量嚴(yán)格計(jì)算最大水平主應(yīng)力方向的公式和方法，可以得到“真實(shí)的”最大水平主應(yīng)力方向.識(shí)別最大水平應(yīng)力的常規(guī)方法是簡(jiǎn)單地采用最大近水平主應(yīng)力的水平投影方向，但該情況只有當(dāng)應(yīng)力張量中的其中一個(gè)主應(yīng)力軸嚴(yán)格垂直的才是適用的.Anderson應(yīng)力模型可以很好地解釋不同構(gòu)造條件下觀測(cè)到的應(yīng)力狀態(tài)，當(dāng)與最大主應(yīng)力軸平行時(shí)，即該軸垂直于地表水平面.一般情況下由震源機(jī)制解估計(jì)得到的應(yīng)力張量，幾乎沒(méi)有應(yīng)力軸是嚴(yán)格垂直于地表水平面的，則不能簡(jiǎn)單地采用最大近水平主應(yīng)力的水平投影，二者之間在方位上可能存在數(shù)十度的偏差.已知構(gòu)造應(yīng)力張量S，由此可以得到3個(gè)主應(yīng)力及其對(duì)應(yīng)的主應(yīng)力矢量s1，s2，s3，在Aki和Richards(1980)定義的地理坐標(biāo)系下(北、東、下為正)

si=[siN,siE,siD]T,

(1)

作用在我們感興趣的垂直平面上法向正應(yīng)力Sn與該法向應(yīng)力的方位角α之間存在如下關(guān)系(Lund and Townend，2007，式(5))：

(2)

式中，α為法向正應(yīng)力Sn的方位角.一般情況下，Sn在區(qū)間[0，π)具有一個(gè)最大值和一個(gè)最小值，使得式(2)為零，分別對(duì)應(yīng)了水平最大主應(yīng)力和最小主應(yīng)力.

反演結(jié)果不確定度的評(píng)定對(duì)理解應(yīng)力場(chǎng)的時(shí)間、空間變化具有重要意義.通常震源機(jī)制解參數(shù)具有較大的不確定性(25°或更大)，需要大量的(30個(gè)震源機(jī)制解)、離散分布的數(shù)據(jù)來(lái)準(zhǔn)確地約束應(yīng)力張量的方向.考察應(yīng)力狀態(tài)的空間分布特征時(shí)，通常的方法是將區(qū)域分成很多小區(qū)，對(duì)每個(gè)小區(qū)的震源機(jī)制擬合應(yīng)力張量，分區(qū)方式不同和分區(qū)大小的變化都可能對(duì)反演的應(yīng)力張量結(jié)果產(chǎn)生影響(Maury et al., 2013).近年來(lái)使用震源機(jī)制解推斷構(gòu)造應(yīng)力場(chǎng)特征的方法又有了新的發(fā)展，一方面基于大量的資料進(jìn)行小區(qū)域精細(xì)構(gòu)造應(yīng)力場(chǎng)的研究(Hardebeck and Hauksson，2001；萬(wàn)永革等，2011b)；另一方面通過(guò)增加阻尼約束應(yīng)用概率信息準(zhǔn)則等數(shù)學(xué)手段來(lái)更好地約束反演結(jié)果(Hardebeck and Michalel，2006；Arnold and Townend，2007).

因此，我們使用MSATSI計(jì)算應(yīng)力場(chǎng)方向的主要步驟如下：(1)確定反演中的最佳阻尼參數(shù)e，阻尼系數(shù)e控制著理論值與觀測(cè)數(shù)據(jù)之間錯(cuò)配值和應(yīng)力反演模型長(zhǎng)度的相對(duì)權(quán)重，它的取值對(duì)反演結(jié)果起著至關(guān)重要的作用；(2)利用SATSI對(duì)所有網(wǎng)格同時(shí)進(jìn)行反演，確定每個(gè)網(wǎng)格的應(yīng)力張量方向和相對(duì)應(yīng)力大小值R；(3)利用bootstrap對(duì)每個(gè)網(wǎng)格的數(shù)據(jù)進(jìn)行重采樣，在特定置信區(qū)間內(nèi)選擇重采樣數(shù)據(jù)進(jìn)行不確定性評(píng)價(jià)；(4)利用msatsi_plot對(duì)應(yīng)力場(chǎng)反演結(jié)果根據(jù)給定參數(shù)進(jìn)行可視化，水平最大主應(yīng)力方向的求解在畫(huà)圖程序中給出，并在特定置信區(qū)間內(nèi)選擇重采樣數(shù)據(jù)進(jìn)行不確定評(píng)價(jià).

3 震源機(jī)制解的總體特征

參照世界應(yīng)力圖的劃分原則(Zoback，1992)，根據(jù)震源機(jī)制解3個(gè)應(yīng)力軸傾角的大小，將震源機(jī)制解類型(即其反映的應(yīng)力狀態(tài))分為6種：正斷型(NF)、正走滑型(NS)、走滑型(SS)、逆走滑型(TS)、逆斷型(TF)和無(wú)法確定型(U)，具體分類標(biāo)準(zhǔn)見(jiàn)表1.結(jié)果表明：走滑斷層型地震為413個(gè)，占50.4%，具有一定走滑分量的正斷層和正斷層類型為91個(gè)，占11.1%，具有一定走滑分量的逆斷層和逆斷層為257個(gè)，占31.4%，無(wú)法確定型為58個(gè)，占7.1%.

從圖1可以看出，大多數(shù)地震都發(fā)生在活動(dòng)地塊邊界帶上，這與張國(guó)民等(2005)得到中國(guó)大陸活動(dòng)地塊邊界帶與強(qiáng)震分布的關(guān)系一致.活動(dòng)地塊具有分級(jí)性，高級(jí)別地塊內(nèi)部可能存在次級(jí)地塊，但不同地塊之間或不同級(jí)別地塊之間的構(gòu)造變形在更大區(qū)域框架下具有協(xié)調(diào)性.地塊內(nèi)部的構(gòu)造活動(dòng)有兩種形式：一種是相對(duì)穩(wěn)定，不發(fā)生大幅度構(gòu)造變形；另一種是內(nèi)部次級(jí)地塊之間發(fā)生相對(duì)運(yùn)動(dòng)，具有一定的構(gòu)造活動(dòng)性，但不論是其活動(dòng)強(qiáng)度還是頻度都遠(yuǎn)小于邊界活動(dòng)構(gòu)造帶.

表1 震源機(jī)制解類型劃分表Table 1 Categories of tectonic stress regime for focal mechanism

走滑類型(圖1a)遍布整個(gè)南北地震帶區(qū)域，尤以南段最多，主要分布在川滇塊體的邊界和內(nèi)部，印度板塊向北東推擠以及青藏高原物質(zhì)向東擴(kuò)張導(dǎo)致川滇菱形塊體上部地殼沿?cái)鄬拥幕?Liu et al., 2014)，引起塊體內(nèi)部次級(jí)斷裂的左旋走滑以及塊體內(nèi)部各構(gòu)造單元之間的差異滑動(dòng)造成區(qū)域內(nèi)呈現(xiàn)走滑的應(yīng)力狀態(tài)，川滇地塊本身相對(duì)于鄰區(qū)向東南的滑脫也促成了其應(yīng)力場(chǎng)的走滑性質(zhì)，造成這一地區(qū)走滑類型地震的彌散分布.龍門(mén)山斷裂帶上的走滑類型地震大多發(fā)生在汶川地震的破裂區(qū)，屬于汶川地震的余震序列，受控于主破裂的滑動(dòng)，表明余震區(qū)應(yīng)力仍處于調(diào)整狀態(tài).正斷層類型(圖1b)主要分布在川滇塊體北部和西南部、龍門(mén)山斷裂帶以及銀川地塹.青藏高原中部羌塘地塊的東向擴(kuò)張?jiān)谒拇ǖ嘏_(tái)強(qiáng)烈受阻引起被動(dòng)式的轉(zhuǎn)向，數(shù)值模擬(朱守彪和石耀霖，2004)和高分辨率地震學(xué)觀測(cè)(Liu et al., 2014)均表明青藏高原東部下地殼易于流動(dòng)，從而對(duì)上地殼有拖曳作用，多種力源的組合決定了川滇塊體北部的應(yīng)力狀態(tài)；川滇塊體的西南部位于一種特殊構(gòu)造背景下的擴(kuò)展裂陷區(qū)，既是紅河斷裂帶的尾端拉伸區(qū)(虢順民等，1996)，又有橫切川滇菱形塊體的麗江—小金河(逆左旋走滑)斷裂在區(qū)內(nèi)與之交會(huì)，青藏高原向東擴(kuò)張引起的中下地殼物質(zhì)流動(dòng)對(duì)上地殼拖曳作用與構(gòu)造部位的特殊性疊加形成了現(xiàn)今顯著的拉張應(yīng)力狀態(tài).逆斷層類型(圖1c)主要集中在龍門(mén)山斷裂帶、馬邊—鹽津斷裂帶、青藏高原東北緣等區(qū)域.汶川地震和蘆山地震發(fā)生在青藏高原東邊界的龍門(mén)山逆沖推覆斷裂帶上，呈現(xiàn)地殼明顯縮短的擠壓破裂帶，其余震序列的震源機(jī)制大多為逆沖類型；馬邊—鹽津斷裂帶表現(xiàn)出左旋走滑-擠壓逆沖的晚第四紀(jì)活動(dòng)特征，NE向昭通、蓮峰斷裂帶表現(xiàn)出右旋走滑兼逆沖、或者以逆沖為主要錯(cuò)動(dòng)方式，震源機(jī)制解類型與構(gòu)造背景一致；青藏高原東北緣的逆斷層型地震顯示了在印度板塊NE向歐亞大陸擠壓碰撞作用下，青藏高原在不斷隆升的構(gòu)造運(yùn)動(dòng)特征.無(wú)法確定型震源機(jī)制解(圖1d)主要集中分布在龍門(mén)山斷裂帶，基本上屬于汶川地震的余震活動(dòng)，反映了不同來(lái)源的動(dòng)力作用使主應(yīng)力場(chǎng)方向傾斜偏離了水平面和垂直面而引起的變形.綜上所述，南北地震帶上的活動(dòng)構(gòu)造和塊體的運(yùn)動(dòng)特性、構(gòu)造變動(dòng)樣式與震源機(jī)制類型之間存在較好的一致性，表明區(qū)域構(gòu)造應(yīng)力場(chǎng)和構(gòu)造運(yùn)動(dòng)特征具有很好的關(guān)聯(lián)性.

為了研究南北地震帶構(gòu)造應(yīng)力場(chǎng)分布的方向特征，圖2給出了研究區(qū)內(nèi)819個(gè)震源機(jī)制解的主壓應(yīng)力P軸和張應(yīng)力T軸的水平投影分布.圖中線段的長(zhǎng)短與P、T軸的傾角有關(guān)，當(dāng)P、T軸平行于水平面時(shí)線段最長(zhǎng)，反之當(dāng)P、T軸垂直于地表時(shí)，其投影成為一個(gè)點(diǎn).P、T軸的水平投影反映了震源應(yīng)力場(chǎng)主壓應(yīng)力軸和張應(yīng)力軸的地表水平投影方位.P軸方位(圖2a)從北向南呈現(xiàn)規(guī)律的變化特征：NNE向—NE向—近EW向—NW向—近NS向—NNE向，表明來(lái)自印度板塊的北北東或北東向的水平擠壓應(yīng)力和青藏高原物質(zhì)向東擴(kuò)張引起的下部地殼物質(zhì)流動(dòng)和上部地殼沿大型走滑斷裂南東向塊體滑移(Liu et al., 2014)的復(fù)合作用控制了南北地震帶的巖石圈應(yīng)力場(chǎng).T軸方位(圖2b)在川滇塊體及鄰近地區(qū)空間差異特征最為突出，存在順時(shí)針旋轉(zhuǎn)的趨勢(shì)，T軸方位從塊體北部的NS向越過(guò)麗江—小金河斷裂后突然轉(zhuǎn)為近EW向拉張，穿過(guò)紅河斷裂后在滇西塊體轉(zhuǎn)為近EW—NWW向.T軸方位的突然轉(zhuǎn)向可能與橫切川滇塊體的麗江—小金河斷裂帶有關(guān)，該斷裂帶在重力、地殼厚度和人工地震剖面上均顯示為一明顯的線性異常帶，這種深位凹凸不平的梯度轉(zhuǎn)折正是阻擋中上地殼物質(zhì)塑性流動(dòng)或彈塑性應(yīng)變傳播的主要屏障，自然會(huì)屏蔽和吸收一部分高原物質(zhì)沿川滇塊體自北向南的擠出和流動(dòng)(向宏發(fā)等，2002).同時(shí)，緬甸弧后擴(kuò)張的遠(yuǎn)場(chǎng)效應(yīng)影響也不能忽視.

圖2 南北地震帶震源機(jī)制解P(a)、T(b)軸水平分量投影分布(圖例同圖1)Fig.2 Horizontal projections of principal compressive (P) axes and extensional (T) axes of focal mechanism solutions in North-South Seismic Belt (Legend same as Fig.1)

龍門(mén)山斷裂帶NWW—SEE向P軸緩傾角而T軸高傾角分布，揭示了巴顏喀拉塊體與華南地塊的相互作用方式——強(qiáng)烈逆沖.川滇塊體東邊界的安寧河—?jiǎng)t木河—小江斷裂帶、大涼山次級(jí)塊體的馬邊—鹽津斷裂帶、蓮峰—昭通斷裂帶上的P軸均為NW向緩傾角分布，震源機(jī)制表明川滇塊體SE向滑移矢量部分被分解到大涼山次級(jí)塊體，繼續(xù)向SE向滑移受到堅(jiān)固的華南地塊的阻擋而產(chǎn)生逆沖兼有少量走滑分量的地震.2012年9月7日云南彝良M5.7地震為蓮峰—昭通斷裂帶有文字記載以來(lái)的最大地震，震源機(jī)制解反映出SEE向水平擠壓的應(yīng)力場(chǎng)(呂堅(jiān)等，2013)，活動(dòng)構(gòu)造野外調(diào)查和高分辨率衛(wèi)星影像分析得到的該斷裂帶的活動(dòng)方式與震源機(jī)制解的錯(cuò)動(dòng)方式基本一致，由此可推知上述斷裂活動(dòng)和變形的動(dòng)力源來(lái)自川滇塊體的SSE向運(yùn)動(dòng)(聞學(xué)澤等，2013).

4 應(yīng)力張量反演

我們將整個(gè)南北帶劃分成2°×2°的網(wǎng)格，利用MSATSI軟件對(duì)每個(gè)網(wǎng)格反演最佳擬合應(yīng)力張量.Hardebeck和Michael(2006)對(duì)不同的數(shù)據(jù)選取方法對(duì)反演結(jié)果的影響做過(guò)對(duì)比研究，發(fā)現(xiàn)加入阻尼的應(yīng)力場(chǎng)反演能消除傳統(tǒng)應(yīng)力場(chǎng)反演方法中應(yīng)力旋轉(zhuǎn)的人為因素，同時(shí)消除簡(jiǎn)單平滑模型和滑動(dòng)窗反演方法產(chǎn)生的較為明顯的真實(shí)應(yīng)力旋轉(zhuǎn).因此，本研究?jī)H對(duì)2°×2°網(wǎng)格內(nèi)的震源機(jī)制解進(jìn)行應(yīng)力張量反演，不進(jìn)行滑動(dòng)和增加周圍網(wǎng)格的數(shù)據(jù).

阻尼系數(shù)e控制著理論值與觀測(cè)數(shù)據(jù)之間錯(cuò)配值和應(yīng)力反演模型長(zhǎng)度的相對(duì)權(quán)重，簡(jiǎn)化模型，則錯(cuò)配值升高，反演誤差增大；反之提高錯(cuò)配的相對(duì)權(quán)重，反演誤差減小，模型逐漸變得復(fù)雜，甚至失去了阻尼約束的意義.通常根據(jù)圖3所示的折中曲線選取最佳的阻尼系數(shù)，折中曲線的拐點(diǎn)在e=1.6處，表明低于該相對(duì)權(quán)重值，提高模型復(fù)雜程序?qū)Ψ囱輾埐畹母纳谱饔貌淮?；而提高相?duì)權(quán)重，隨著模型的簡(jiǎn)化，反演誤差會(huì)急劇增加.本研究中我們采用e=1.6作為最佳阻尼系數(shù).

圖3 模型長(zhǎng)度與數(shù)據(jù)擬合誤差之間的折中曲線圖

在震源機(jī)制應(yīng)力反演程序中有4個(gè)待求參數(shù)，最佳應(yīng)力模型3個(gè)主應(yīng)力軸σ1，σ2和σ3(σ1≥σ2≥σ3)的方位角和傾角，同時(shí)還可以得到相對(duì)應(yīng)力大小值R，R=(σ1-σ2)/(σ1-σ3).相對(duì)應(yīng)力大小是中等主應(yīng)力軸σ2接近最大主應(yīng)力軸σ1或者接近最小主應(yīng)力軸σ3的度量，有助于我們區(qū)分應(yīng)力場(chǎng)的類型.地殼和地幔內(nèi)的應(yīng)力場(chǎng)或許存在差異，所以我們僅選取深度在60 km以內(nèi)的數(shù)據(jù)，對(duì)于震源機(jī)制解個(gè)數(shù)少于6個(gè)的網(wǎng)格，不反演其應(yīng)力張量.在95%的置信區(qū)間對(duì)原始數(shù)據(jù)進(jìn)行600次bootstrap重采樣利用2D阻尼法向應(yīng)力反演技術(shù)獲得每個(gè)網(wǎng)格的最佳擬合構(gòu)造應(yīng)力張量，然后計(jì)算每個(gè)網(wǎng)格的水平最大主應(yīng)力方向(表2和圖4).為了平滑地震空白的區(qū)域，沒(méi)有數(shù)據(jù)的格點(diǎn)也包括在反演中，但這些點(diǎn)上的結(jié)果沒(méi)有顯示.

圖4給出以全球應(yīng)力圖WSM(Heidbach et al.，2010)分類標(biāo)準(zhǔn)顯示的從鄂爾多斯塊體西緣經(jīng)龍門(mén)山斷裂帶、川滇塊體到滇緬塊體最大水平主應(yīng)力的方向變化與應(yīng)力狀態(tài)，直觀地反映了南北地震帶背景應(yīng)力場(chǎng)的作用方式.南北地震帶的最大水平主應(yīng)力方向具有明顯的分區(qū)特征，與P軸方位有著基本一致的變化規(guī)律.大體分為三段，北段為NE向走滑類型的應(yīng)力狀態(tài)，中段為NEE—EW—NWW向的逆沖類型，南段為SE—SSE—NS—NNE向走滑和正斷類型，在川滇塊體的北部和西邊界應(yīng)力狀態(tài)為EW—SE—SSE的正斷層類型.反演得到的應(yīng)力狀態(tài)基本上與各種類型地震的發(fā)生比較吻合，也進(jìn)一步證明結(jié)果的可靠性.水平最大主壓應(yīng)力方向與前人對(duì)南北地震帶研究得到的結(jié)果一致(許忠淮等，1987；闞榮舉等，1977；成爾林，1981；秦保燕等，1986；杜興信和邵輝成，1999；Wan，2010；楊樹(shù)新等，2012)，在南北地震帶南段水平最大主壓應(yīng)力方向存在順時(shí)針旋轉(zhuǎn)的趨勢(shì)，與GPS觀測(cè)得到的地表水平速度場(chǎng)基本一致(Shen et al.，2005；Zhang，2013)；在空間組成一個(gè)倒“V”字形，與徐紀(jì)人和尾池和夫(1995)發(fā)現(xiàn)P軸方位空間分布特征類似.

鄂爾多斯西緣的銀川—吉蘭泰斷陷帶內(nèi)發(fā)育著第四紀(jì)以來(lái)活動(dòng)性最強(qiáng)的NNE至NE向斷裂，控制著地塹的形成，在NE向主壓應(yīng)力作用下近NS向斷陷帶做右旋剪切運(yùn)動(dòng)，致使盆地內(nèi)部和邊緣發(fā)生走滑和正斷層類型地震.1920年海原地震(Selva and Marzocchi，2004)發(fā)生在海原斷裂上, 是左旋走滑型破裂，基本與所屬區(qū)域的應(yīng)力狀態(tài)吻合.祁連山、柴達(dá)木盆地和六盤(pán)山構(gòu)成了青藏高原的東北邊緣，該地區(qū)地震活動(dòng)強(qiáng)烈，震源機(jī)制解表明地震破裂以逆沖和左旋走滑為特征，與地表活動(dòng)斷裂研究結(jié)果一致(Meyer et al.，1998).鄂爾多斯塊體西南緣NEE向的逆沖類型應(yīng)力狀態(tài)反映了青藏塊體對(duì)鄂爾多斯塊體的直接作用，六盤(pán)山斷裂帶第四紀(jì)以來(lái)以強(qiáng)烈的擠壓性質(zhì)及左旋走滑運(yùn)動(dòng)區(qū)別于鄂爾多斯邊緣的其他活動(dòng)斷裂(國(guó)家地震局“鄂爾多斯周緣活動(dòng)斷裂系”課題組，1988).印度板塊對(duì)歐亞板塊NNE向的推擠在青藏高原的邊界遇到剛性塊體的阻擋，逆沖類型地震主要發(fā)生在高原周邊，反映了高原向周邊的擠壓和縮短作用；走滑斷裂發(fā)育在高原的不同海拔不同部位，北部是左旋走滑運(yùn)動(dòng)，南部是右旋走滑運(yùn)動(dòng)，反映了高原內(nèi)部物質(zhì)東流所導(dǎo)致的巨型左旋和右旋剪切作用(Gan et al.，2007； Zhang，2013)，在巴顏喀拉塊體的南邊界鮮水河斷裂上走滑類型地震十分發(fā)育；青藏高原的物質(zhì)向東擴(kuò)張由于受到穩(wěn)定的四川盆地的阻擋，一部分東流物質(zhì)使得高角度(50°～70°)、犁狀的龍門(mén)山斷層面上的正應(yīng)力不斷增大，導(dǎo)致該斷層的閉鎖性逐步加強(qiáng)，另一方面青藏高原較柔軟的下地殼物質(zhì)的不斷向東運(yùn)動(dòng)對(duì)龍門(mén)山斷裂帶上盤(pán)的推擠作用不斷加強(qiáng)，導(dǎo)致斷裂帶上剪應(yīng)力越來(lái)越大；當(dāng)剪應(yīng)力超過(guò)摩擦強(qiáng)度時(shí)，斷層解鎖產(chǎn)生滑動(dòng)，發(fā)生汶川地震(朱守彪和張培震，2009).

表2 南北地震帶應(yīng)力場(chǎng)反演結(jié)果(只列出有結(jié)果的網(wǎng)格)Table 2 Results of stress field inverted in the North-South Seismic Belt (only give the grids with results)

闞榮舉等(1977)指出，由鮮水河斷裂帶、安寧河斷裂帶、小江斷裂帶、通?！瓟嗔褞Ш图t河斷裂帶所圍成的一個(gè)菱形塊體，由于受到印度洋板塊對(duì)歐亞板塊的推擠影響，向SE方向運(yùn)動(dòng)并控制著那里的應(yīng)力場(chǎng).川滇塊體及鄰區(qū)的水平最大主壓應(yīng)力方向呈現(xiàn)有序的變化特征，北部EW向的主壓應(yīng)力方向顯示出羌塘塊體曾發(fā)生過(guò)或正經(jīng)歷著東向滑移，與劉平江等(2007)的分析結(jié)果一致；中部和南部水平最大主壓應(yīng)力方向?yàn)镾E—SSE向，青藏高原物質(zhì)東向滑移遇到華南塊體的強(qiáng)烈阻擋，使塊體滑動(dòng)方向被迫向南偏轉(zhuǎn)，引起川西北次級(jí)塊體SE向的滑移，而塊體東邊界斷裂走向從北西偏轉(zhuǎn)成近南北，滇中次級(jí)塊體被迫SSE向平移.滇緬地塊水平最大主壓應(yīng)力方向?yàn)镹NE向，應(yīng)力狀態(tài)為走滑類型，與謝富仁等(2001)的分析結(jié)果一致，該地區(qū)遭受著來(lái)自印度板塊側(cè)壓、川滇塊體楔入和東南亞板塊向北的聯(lián)合作用(謝富仁等，1999).川滇塊體及鄰區(qū)的運(yùn)動(dòng)圖像與構(gòu)造地質(zhì)學(xué)觀測(cè)結(jié)果相一致(徐錫偉等，2003).

圖4 南北地震帶應(yīng)力場(chǎng)反演結(jié)果

值得一提的是，川滇塊體的西邊界的應(yīng)力結(jié)構(gòu)為EW—SE—SSE向的正斷層類型，該區(qū)域與高分辨率地震學(xué)觀測(cè)得到的中下地殼低速帶范圍基本吻合，青藏高原向東擴(kuò)張的塑性物質(zhì)流與橫向邊界(麗江—小金河斷裂帶)的弱化易于應(yīng)變能的釋放，在局部地區(qū)使NS向拉張的正斷層向EW向拉張正斷層轉(zhuǎn)變(Liu et al., 2014).

5 結(jié)論與討論

震源機(jī)制解是研究構(gòu)造應(yīng)力場(chǎng)特征的基礎(chǔ)資料，本研究收集和計(jì)算了南北地震帶上819條震源機(jī)制解.通過(guò)MSATSI軟件(Martínez-Garzón et al.，2014)分析了南北地震帶現(xiàn)今地殼應(yīng)力狀態(tài)總體特征與不同分段地殼應(yīng)力狀態(tài)特征.主要結(jié)論如下：

(1) 南北地震帶上的震源機(jī)制解類型主要以走滑類型和逆沖類型為主.走滑類型地震主要分布在塊體邊界的大型走滑斷裂帶上，反映了青藏高原內(nèi)部物質(zhì)向東擴(kuò)張所導(dǎo)致的巨型左旋和右旋剪切作用.南北地震帶南段走滑類型最多，青藏高原向東擴(kuò)張引起川滇菱形塊體的上地殼沿邊界斷裂南南東向的滑移導(dǎo)致其內(nèi)部次級(jí)斷裂的左旋走滑，塊體本身相對(duì)于鄰區(qū)向東南的滑脫也促成了其應(yīng)力場(chǎng)的走滑性質(zhì)，再加上印度板塊的側(cè)向剪切擠壓，造成這一地區(qū)走滑類型地震的彌散分布.逆沖類型主要分布在柴達(dá)木盆地、六盤(pán)山斷裂帶和龍門(mén)山斷裂帶上，反映了青藏高原向周邊的擠壓隆升和縮短作用.

(2) 南北地震帶P軸水平投影從北向南呈現(xiàn)規(guī)律的變化特征：NNE向—NE向—近EW向—NW向—近NS向—NNE向，表明來(lái)自印度板塊的NNE或NE向的水平擠壓應(yīng)力和青藏高原物質(zhì)東向滑移沿大型走滑斷裂帶向SE向平移的復(fù)合作用控制了南北地震帶的巖石圈應(yīng)力場(chǎng).

(3) 南北地震帶T軸水平投影在川滇塊體及鄰近地區(qū)空間差異特征最為突出，T軸方位從塊體北部的NS向越過(guò)麗江—小金河斷裂后突然轉(zhuǎn)為近EW向拉張，而水平最大主應(yīng)力方向只是存在漸變順時(shí)針旋轉(zhuǎn)的趨勢(shì)，并且與GPS觀測(cè)得到的地表水平速度場(chǎng)基本一致.T軸方位的突然轉(zhuǎn)向可能與橫切川滇塊體的麗江—小金河斷裂帶有關(guān)，青藏高原向東擴(kuò)張引起的中下地殼塑性物質(zhì)流與橫向邊界(麗江—小金河斷裂帶)的弱化易于應(yīng)變能的積累，擠壓強(qiáng)度的減弱(楊樹(shù)新等，2012)和下地殼對(duì)上地殼的拖曳作用，在局部地區(qū)使NS向拉張的正斷層向EW向拉張正斷層轉(zhuǎn)變.同時(shí)，緬甸弧后擴(kuò)張的遠(yuǎn)場(chǎng)效應(yīng)影響也不能忽視.

(4) 南北地震帶應(yīng)力狀態(tài)的分布特征：北段為NE向走滑類型，中段為NEE—EW—SEE向逆沖類型，南段為SE—SSE—NNE向走滑和正斷類型；由北向南分別對(duì)應(yīng)于鄂爾多斯塊體西緣的吉蘭泰—銀川斷陷盆地，六盤(pán)山斷裂帶，龍門(mén)山斷裂帶和川滇塊體等，震源機(jī)制解類型比較符合反演得到的應(yīng)力狀態(tài).川滇塊體的西邊界的應(yīng)力結(jié)構(gòu)為EW—SE—SSE向的正斷層類型，在深部結(jié)構(gòu)上對(duì)應(yīng)中下地殼低速帶，淺部構(gòu)造對(duì)應(yīng)紅河斷裂帶尾端拉伸區(qū)，深淺構(gòu)造的聯(lián)合作用導(dǎo)致了該區(qū)域獨(dú)特的應(yīng)力狀態(tài).

本研究得到的水平最大主壓應(yīng)力方向與張培震等(2013)綜合考慮斷裂錯(cuò)動(dòng)、地塊旋轉(zhuǎn)和內(nèi)部變形所獲得的青藏高原東部的構(gòu)造變形速度場(chǎng)基本一致，表明印度板塊對(duì)中國(guó)大陸的碰撞和推擠作用是中國(guó)大陸晚新生代構(gòu)造變形最重要的動(dòng)力作用，而不同活動(dòng)地塊本身的性質(zhì)決定著地塊的整體性和變形方式，活動(dòng)地塊邊界帶由于其差異運(yùn)動(dòng)強(qiáng)烈而構(gòu)造變形非連續(xù)性最強(qiáng)，最有利于應(yīng)力的高度積累而孕育強(qiáng)震.先前研究中利用初動(dòng)和震源機(jī)制資料采用不同反演方法也得到相同的應(yīng)力場(chǎng)特征(許忠淮等，1987；杜興信和邵輝成等，1999；Wan，2010).與先前研究給出同樣的應(yīng)力場(chǎng)分布特征，證實(shí)了本研究所使用的方法和數(shù)據(jù).本研究給出的應(yīng)力場(chǎng)結(jié)果與其他方法獲得的結(jié)果也是一致的，如GPS觀測(cè)結(jié)果(Wang et al.，2001)及其應(yīng)變率場(chǎng)(Zhang，2013)以及GPS、地質(zhì)和剪切波分裂數(shù)據(jù)的綜合研究結(jié)果(Flesch et al.，2005).以上由大量震源機(jī)制解反演得到的南北地震帶地殼應(yīng)力場(chǎng)宏觀分布的一些規(guī)律性的認(rèn)識(shí)，可為地球動(dòng)力學(xué)過(guò)程的模擬和活動(dòng)斷層滑動(dòng)性質(zhì)的厘定提供參考.鑒于本研究只是使用了4級(jí)以上地震的震源機(jī)制，后續(xù)仍需要填加質(zhì)量可靠的中小地震震源機(jī)制解來(lái)豐富數(shù)據(jù)的多樣性，增加數(shù)據(jù)分布的均勻性，使我們可以更加精細(xì)地討論南北地震帶的構(gòu)造應(yīng)力狀態(tài).

致謝 中國(guó)地震局地球物理研究所國(guó)家數(shù)字測(cè)震臺(tái)網(wǎng)數(shù)據(jù)備份中心(doi:10.7914/SN/CB)為本研究提供地震波形數(shù)據(jù)，文中圖件利用Generic Mapping Tools (GMT)(Wessel and Smith，1995)繪制，審稿專家給出寶貴修改意見(jiàn)，在此一并表示衷心感謝！

Abers G A, Gephart J W. 2001. Direct inversion of earthquake first motions for both the stress tensor and focal mechanisms and application to southern California.J.Geophys.Res., 106(B11): 26523-26540, doi: 10.1029/2001JB000437.

Aki K, Richards P G. 1980. Quantitative Seismology: Theory and Methods, Vol. Ⅰ. San Francisco: W H Freeman, 1-557.

Angelier J. 1984. Tectonic analysis of fault slip data sets.J.Geophys.Res., 89(B7): 5835-5848, doi: 10.1029/JB089iB07p05835. Arnold R, Townend J. 2007.

A Bayesian approach to estimating tectonic stress from seismological data.Geophys.J.Int., 170(3): 1336-1356, doi: 10.1111/j.1365-246X.2007.03485.x.

Bott M H P. 1959. The mechanics of oblique slip faulting.Geol.Mag., 96(2): 109-117.

Cheng E L. 1981. Recent tectonic stress field and tectonic movement of the Sichuan province and its vicinity.ActaSeismologicaSinica(in Chinese), 3(3): 231-241.

Cheng W Z, Diao G L, Lü G P, et al. 2003. Focal mechanisms, displacement rate and mode of motion of the Sichuan-Yunnan block.SeismologyandGeology(in Chinese), 25(1): 71-87.

Du X X, Shao H C. 1999. Modern tectonic stress field in the Chinese mainland inversed from focal mechanism solutions.ActaSeismologicaSinica(in Chinese), 21(4): 354-360.

Dziewonski A M, Chou T A, Woodhouse J H. 1981. Determination of earthquake source parameters from waveform data for studies of global and regional seismicity.J.Geophys.Res., 86(B4): 2825-2852, doi: 10.1029/JB086iB04p02825.

Ekstr?m G, Nettles M, Dziewoński A M. 2012. The global CMT project 2004—2010: Centroid-moment tensors for 13017 earthquakes.Phys.EarthPlanet.Inter., 200-201: 1-9, doi: 10.1016/j.pepi.2012.04.002.

Flesch L M, Holt W E, Silver P G, et al. 2005. Constraining the extent of crust-mantle coupling in central Asia using GPS, geologic, and shear wave splitting data.EarthPlanet.Sci.Lett., 238(1-2): 248-268.

Gan W J, Zhang P Z, Shen Z K, et al. 2007. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements.J.Geophys.Res., 112: B08416, doi: 10.1029/2005JB004120.

Gephart J W, Forsyth D W. 1984. An improved method for determining the regional stress tensor using earthquake focal mechanism data: Application to the San Fernando earthquake sequence.J.Geophys.Res., 89(B11): 9305-9320, doi: 10.1029/JB089iB11p09305.

Gephart J W. 1990. Stress and the direction of slip on fault planes.Tectonics, 9(4): 845-858, doi: 10.1029/TC009i004p00845.

Guo S M, Xiang H F, Ji F J, et al. 1996. A study on the relation between Quaternary right-lateral slip and tip extension along the Honghe Fault.SeismologyandGeology(in Chinese), 18(4): 301-309.

Guo X Y, Chen X Z, Wang S W, et al. 2014. Focal mechanism of small and moderate earthquakes and tectonic stress field in Sichuan-Yunnan Areas.ChinaEarthquakeEngineeringJournal(in Chinese), 36(3): 599-607.

Hao P, Lü X J, Tian Q J, et al. 2012. Contemporary tectonic stress field on boundaries of active tectonic blocks in and around western China.ActaSeismologicaSinica(in Chinese), 34(4): 439-450.

Hardebeck J L, Hauksson E. 2001. Crustal stress field in southern California and its implications for fault mechanics.J.Geophys.Res., 106(B10): 21859-21882, doi: 10.1029/2001JB000292.

Hardebeck J L, Michael A J. 2006. Damped regional-scale stress inversions: Methodology and examples for southern California and the Coalinga aftershock sequence.J.Geophys.Res., 111(B11): B11310, doi: 10.1029/2005JB004144.

Heidbach O, Tingay M, Barth A, et al. 2010. Global crustal stress pattern based on the World Stress Map database release 2008.Tectonophysics, 482(1-4): 3-15, doi: 10.1016/j.tecto.2009.07.023.

Huang Z W, Su Y, Peng Y J, et al. 2003. Rayleigh wave tomography of China and adjacent regions.J.Geophys.Res., 108(B2): 2073, doi: 10.1029/2001JB001696.

Kan R J, Zhang S C, Yan F T, et al. 1977. Present tectonic stress field and its relation to the characteristics of recent tectonic activity in southwestern China.ChineseJ.Geophys. (ActaGeophysicaSinica) (in Chinese), 20(2): 96-109.

Liu P J, Diao G L, Ning J Y. 2007. Fault plane solutions in Sichuan Yunnan rhombic block and their dynamic implications.ActaSeismologicaSinica(in Chinese), 29(5): 449-458.

Liu Q Y, van der Hilst R D, Li Y, et al. 2014. Eastward expansion of the Tibetan Plateau by crustal flow and strain partitioning across faults.NatureGeoscience, 7(5): 361-365, doi: 10.1038/ngeo2130.

Lund B, Townend J. 2007. Calculating horizontal stress orientations with full or partial knowledge of the tectonic stress tensor.Geophys.J.Int., 170(3): 1328-1335.

Lü J, Zheng X F, Xiao J, et al. 2013. Rupture characteristics and seismogenic structures of theMS5.7 andMS5.6 Yiliang earthquakes of Sep. 7, 2012.ChineseJ.Geophys. (in Chinese), 56(8): 2645-2654, doi: 10.6038/cjg20130814.

Ma X Y. 1989. Lithospheric Dynamics Atlas of China (in Chinese). Beijing: China Cartographic Publishing House.

Martínez-Garzón P, Kwiatek G, Ickrath M, et al. 2014. MSATSI: A MATLAB package for stress inversion combining solid classic methodology, a new simplified user-handling, and a visualization tool.Seismol.Res.Lett., 85(4): 896-904.

Maury J, Cornet F H, Dorbath L. 2013. A review of methods for determining stress fields from earthquakes focal mechanisms; Application to the Sierentz 1980 seismic crisis (Upper Rhine graben).Bull.Soc.Géol.France, 184(4-5): 319-334.

Menke W. 1989. Geophysical Data Analysis: Discrete Inverse Theory. San Diego, California: Academic Press.

Meyer B, Tapponnier P, Bourjot L, et al. 1998. Crustal thickening in Gansu-Qinghai, lithospheric mantle subduction, and oblique, strike-slip controlled growth of the Tibet Plateau.Geophys.J.Int., 135: 1-47.

Michael A J. 1984. Determination of stress from slip data: Faults and folds.J.Geophys.Res., 89(B13): 11517-11526, doi: 10.1029/JB089iB13p11517.

Michael A J. 1987. Use of focal mechanisms to determine stress: A control study.J.Geophys.Res., 92(B1): 357-368.

Peltzer G, Tapponnier P. 1988. Formation and evolution of strike-slip collision: an experimental approach.J.Geophys.Res., 93(B12): 15085-15117, doi: 10.1029/JB093iB12p15085.

Qin B Y, Xu J R, Zhai W J. 1986. Discussion on the current crustal stress field in the north segment of the South-North Seismic Belt.NorthChinaEarthquakeSciences(in Chinese), 4(1): 1-15.

Selva J, Marzocchi W. 2004. Focal parameters, depth estimation, and plane selection of the worldwide shallow seismicity withMS≥7.0 for the period 1900～1976.Geochem.Geophys.Geosyst., 5(5): Q05005, doi: 10.1029/2003GC000669.

Shen Z K, Lü J N, Wang M, et al. 2005. Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau.J.Geophys.Res., 110(B11), doi: 10.1029/2004JB003421.

The Research Group on Active Fault System around Ordos Massif, SSB. 1988. Active Fault System around Ordos Massif (in Chinese). Beijing: Seismological Press.

Wallace R E. 1951. Geometry of shearing stress and relation to faulting.J.Geol., 59(2): 118-130.

Wan Y G. 2010. Contemporary tectonic stress field in China.EarthquakeScience, 23(4): 377-386, doi: 10.1007/s11589-010-0735-5.

Wan Y G, Wu Y M, Sheng S Z, et al. 2011a. Preliminary result of Taiwan 3-D stress field from P wave polarity data.ChineseJ.Geophys. (in Chinese), 54(11): 2809-2818, doi: 10.3969/j.issn.0001-5733.2011.11.011.

Wan Y G, Sheng S Z, Hsu Y J, et al. 2011b. Effect of stress ratio and friction coefficient on composite P wave radiation patterns.ChineseJ.Geophys. (in Chinese), 54(4): 994-1001, doi: 10.3969/j.issn.0001-5733.2011.04.014.

Wang Q, Zhang P Z, Freymueller J, et al. 2001. Present-day crustal deformation in China constrained by Global Positioning System measurements.Science, 294(5542): 574-577.

Wang Z S, Wang Z Y, Gu J P, et al. 1976. A preliminary investigation of the limits and certain features of the North-South Seismic Zone of China.ChineseJ.Geophys. (ActaGeophysicaSinica) (in Chinese), 19(2): 110-117.

Wen X Z, Du F, Yi G X, et al. 2013. Earthquake potential of the Zhaotong and Lianfeng fault zones of the eastern Sichuan-Yunnan border region.ChineseJ.Geophys., 56(10): 3361-3372, doi: 10.6038/cjg20131012.

Wessel P, Smith W H F. 1995. New version of the generic mapping tools.Eos,TransactionsAmericanGeophysicalUnion, 76(33): 329.

Xiang H F, Xu X W, Guo S M, et al. 2002. Sinistral thrusting along the Lijiang-Xiaojinhe Fault since Quaternary and its geologic-tectonic significance: Shielding effect of transverse structure of intracontinental active block.SeismologyandGeology(in Chinese), 24(2): 188-198.

Xie F R, Zhang S M, Dou S Q, et al. 1999. Evolution characteristics of Quaternary tectonic stress field in the north and east margin of Qinghai-Xizang plateau.ActaSeismologicaSinica(in Chinese), 21(5): 501-512.

Xie F R, Su G, Cui X F, et al. 2001. Modern tectonic stress field in southwestern Yunnan, China.ActaSeismologicaSinica(in Chinese), 23(1): 17-23.

Xie F R, Chen Q C, Cui X F, et al. 2003. Research on Crustal Stress State in China and Adjacent Area (in Chinese). Beijing: Geological Publishing House.

Xie F R, Cui X F, Zhao J T, et al. 2004. Regional division of the recent tectonic stress field in China and adjacent areas.ChineseJ.Geophys. (in Chinese), 47(4): 654-662.

Xu J R, Oike Kazuo. 1995. Earthquake mechanisms and its implication for tectonic stress field in the southern part of the North-South Seismic Belt in China.ActaSeismologicaSinica(in Chinese), 17(1): 31-40.

Xu J R, Zhao Z X, Ishikawa Yuzo. 2008. Regional characteristics of crustal stress field and tectonic motions in and around Chinese mainland.ChineseJ.Geophys. (in Chinese), 51(3): 770-781.

Xu P L. 2004. Determination of regional stress tensors from fault-slip data.Geophys.J.Int., 157(3): 1316-1330, doi: 10.1111/j.1365-246X.2004.02271.x.

Xu X W, Wen X Z, Zheng R Z, et al. 2003. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region, China.ScienceinChina(Ser.D), 46(Suppl.2): 210-226.

Xu Z H, Wang S Y, Huang Y R, et al. 1987. Directions of mean stress axes in southwestern China deduced from microearthquake data.ChineseJ.Geophys. (ActaGeophysicaSinica) (in Chinese), 30(5): 476-486.

Xu Z H, Wang S Y, Huang Y R, et al. 1989. The tectonic stress field of Chinese continent deduced from a great number of earthquakes.ChineseJ.Geophys. (ActaGeophysicaSinica) (in Chinese), 32(6): 636-647.

Yan J Q, Shi Z L, Wang S Y, et al. 1979. Some features of the recent tectonic stress field of China and environs.ActaSeismologicaSinica(in Chinese), 1(1): 9-24.

Yang S X, Yao R, Cui X F, et al. 2012. Analysis of the characteristics of measured stress in Chinese mainland and its active blocks and North-South seismic belt.ChineseJ.Geophys. (in Chinese), 55(12): 4207-4217, doi: 10.6038/j.issn.0001-5733.2012.12.032.

Zhang G M, Ma H S, Wang H, et al. 2005. Boundaries between active-tectonic blocks and strong earthquakes in the China mainland.ChineseJ.Geophys. (in Chinese), 48(3): 602-610.

Zhang P Z. 2013. A review on active tectonics and deep crustal processes of the Western Sichuan region, eastern margin of the Tibetan Plateau.Tectonophysics, 584: 7-22, doi: 10.1016/j.tecto.2012.02.021.

Zhang P Z, Deng Q D, Zhang Z Q, et al. 2013. Active faults, earthquake hazards and associated geodynamic processes in continental China.ScienceinChina(EarthScience) (in Chinese), 43(10): 1607-1620.

Zheng X F, Ouyang B, Zhang D N, et al. 2009. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake.ChineseJ.Geophys. (in Chinese), 52(5): 1412-1417, doi: 10.3969/j.issn.0001-5733.2009.05.031.

Zhu S B, Shi Y L. 2004. Genetic algorithm-finite element inversion of drag forces exerted by the lower crust on the upper crust in the Sichuan-Yunnan area.ChineseJ.Geophys. (in Chinese), 47(2): 232-239.

Zhu S B, Zhang P Z. 2009. A study on the dynamical mechanisms of the WenchuanMS8.0 earthquake, 2008.ChineseJ.Geophys. (in Chinese), 52(2): 418-427.

Zoback M L. 1992. First- and second-order patterns of stress in the lithosphere: The World Stress Map Project.J.Geophys.Res., 97(B8): 11703-11728, doi: 10.1029/92JB00132.

附中文參考文獻(xiàn)

成爾林. 1981. 四川及其鄰區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)和現(xiàn)代構(gòu)造運(yùn)動(dòng)特征. 地震學(xué)報(bào), 3(3): 231-241.

程萬(wàn)正, 刁桂苓, 呂弋培等. 2003. 川滇地塊的震源力學(xué)機(jī)制、運(yùn)動(dòng)速率和活動(dòng)方式. 地震地質(zhì), 25(1): 71-87.

杜興信, 邵輝成. 1999. 由震源機(jī)制解反演中國(guó)大陸現(xiàn)代構(gòu)造應(yīng)力場(chǎng). 地震學(xué)報(bào), 21(4): 354-360.

國(guó)家地震局“鄂爾多斯周緣活動(dòng)斷裂系”課題組. 1988. 鄂爾多斯周緣活動(dòng)斷裂系. 北京: 地震出版社.

虢順民, 向宏發(fā), 計(jì)鳳桔等. 1996. 紅河斷裂帶第四紀(jì)右旋走滑與尾端拉張轉(zhuǎn)換關(guān)系研究. 地震地質(zhì), 18(4): 301-309.

郭祥云, 陳學(xué)忠, 王生文等. 2014. 川滇地區(qū)中小地震震源機(jī)制解及構(gòu)造應(yīng)力場(chǎng)的研究. 地震工程學(xué)報(bào), 36(3): 599-607.

郝平, 呂曉健, 田勤儉等. 2012. 中國(guó)西部及鄰區(qū)活動(dòng)地塊邊界帶現(xiàn)代構(gòu)造應(yīng)力場(chǎng). 地震學(xué)報(bào), 34(4): 439-450.

闞榮舉, 張四昌, 晏鳳桐等. 1977. 我國(guó)西南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)與現(xiàn)代構(gòu)造活動(dòng)特征的探討. 地球物理學(xué)報(bào), 20(2): 96-109.

劉平江, 刁桂苓, 寧杰遠(yuǎn). 2007. 川滇地塊的震源機(jī)制解特征及其地球動(dòng)力學(xué)解釋. 地震學(xué)報(bào), 29(5): 449-458.

呂堅(jiān), 鄭秀芬, 肖健等. 2013. 2012年9月7日云南彝良MS5.7、MS5.6地震震源破裂特征與發(fā)震構(gòu)造研究. 地球物理學(xué)報(bào), 56(8): 2645-2654, doi: 10.6038/cjg20130814.

馬杏垣. 1989. 中國(guó)巖石圈動(dòng)力學(xué)地圖集. 北京: 中國(guó)地圖出版社.

秦保燕, 徐紀(jì)人, 翟文杰. 1986. 南北地震帶北段現(xiàn)今地殼應(yīng)力場(chǎng)討論. 華北地震科學(xué), 4(1): 1-15.

萬(wàn)永革, 吳逸民, 盛書(shū)中等. 2011a. P波極性數(shù)據(jù)所揭示的臺(tái)灣地區(qū)三維應(yīng)力結(jié)構(gòu)的初步結(jié)果. 地球物理學(xué)報(bào), 54(11): 2809-2818, doi: 10.3969/j.issn.0001-5733.2011.11.011.

萬(wàn)永革, 盛書(shū)中, 許雅儒等. 2011b. 不同應(yīng)力狀態(tài)和摩擦系數(shù)對(duì)綜合P波輻射花樣影響的模擬研究. 地球物理學(xué)報(bào), 54(4): 994-1001, doi: 10.3969/j.issn.0001-5733.2011.04.014.

王振聲, 王周元, 顧僅萍等. 1976. 中國(guó)南北地震帶的范圍及其活動(dòng)特征初步探討. 地球物理學(xué)報(bào), 19(2): 110-117.

聞學(xué)澤, 杜方, 易桂喜等. 2013. 川滇交界東段昭通、蓮峰斷裂帶的地震危險(xiǎn)背景. 地球物理學(xué)報(bào), 56(10): 3361-3372, doi: 10.6038/cjg20131012.

向宏發(fā), 徐錫偉, 虢順民等. 2002. 麗江—小金河斷裂第四紀(jì)以來(lái)的左旋逆推運(yùn)動(dòng)及其構(gòu)造地質(zhì)意義——陸內(nèi)活動(dòng)地塊橫向構(gòu)造的屏蔽作用. 地震地質(zhì), 24(2): 188-198.

謝富仁, 張世民, 竇素芹等. 1999. 青藏高原北、東邊緣第四紀(jì)構(gòu)造應(yīng)力場(chǎng)演化特征. 地震學(xué)報(bào), 21(5): 501-512.

謝富仁, 蘇剛, 崔效鋒等. 2001. 滇西南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)分析. 地震學(xué)報(bào), 23(1): 17-23.

謝富仁, 陳群策, 崔效鋒等. 2003. 中國(guó)大陸地殼應(yīng)力環(huán)境研究. 北京: 地質(zhì)出版社.

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

徐紀(jì)人, 尾池和夫. 1995. 南北地震帶南段應(yīng)力場(chǎng)特征及其與板塊

運(yùn)動(dòng)的關(guān)系. 地震學(xué)報(bào), 17(1): 31-40.

徐紀(jì)人, 趙志新, 石川有三. 2008. 中國(guó)大陸地殼應(yīng)力場(chǎng)與構(gòu)造運(yùn)動(dòng)區(qū)域特征研究. 地球物理學(xué)報(bào), 51(3): 770-781.

徐錫偉, 聞學(xué)澤, 鄭榮章等. 2003. 川滇地區(qū)活動(dòng)塊體最新構(gòu)造變動(dòng)樣式及其動(dòng)力來(lái)源. 中國(guó)科學(xué)(D輯), 33(增刊): 151-162.

許忠淮, 汪素云, 黃雨蕊等. 1987. 由多個(gè)小震推斷的青甘和川滇地區(qū)地殼應(yīng)力場(chǎng)的方向特征. 地球物理學(xué)報(bào), 30(5): 476-486.

許忠淮, 汪素云, 黃雨蕊等. 1989. 由大量的地震資料推斷的我國(guó)大陸構(gòu)造應(yīng)力場(chǎng). 地球物理學(xué)報(bào), 32(6): 636-647.

鄢家全, 時(shí)振梁, 汪素云等. 1979. 中國(guó)及鄰區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)的區(qū)域特征. 地震學(xué)報(bào), 1(1): 9-24.

楊樹(shù)新, 姚瑞, 崔效鋒等. 2012. 中國(guó)大陸與各活動(dòng)地塊、南北地震帶實(shí)測(cè)應(yīng)力特征分析. 地球物理學(xué)報(bào), 55(12): 4207-4217, doi: 10.6038/j.issn.0001-5733.2012.12.032.

張國(guó)民, 馬宏生, 王輝等. 2005. 中國(guó)大陸活動(dòng)地塊邊界帶與強(qiáng)震活動(dòng). 地球物理學(xué)報(bào), 48(3): 602-610.

張培震, 鄧起東, 張竹琪等. 2013. 中國(guó)大陸的活動(dòng)斷裂、地震災(zāi)害及其動(dòng)力過(guò)程. 中國(guó)科學(xué): 地球科學(xué), 43(10): 1607-1620.

鄭秀芬, 歐陽(yáng)飚, 張東寧等. 2009. “國(guó)家數(shù)字測(cè)震臺(tái)網(wǎng)數(shù)據(jù)備份中心”技術(shù)系統(tǒng)建設(shè)及其對(duì)汶川大地震研究的數(shù)據(jù)支撐. 地球物理學(xué)報(bào), 52(5): 1412-1417, doi: 10.3969/j.issn.0001-5733.2009.05.031.

朱守彪, 石耀霖. 2004. 用遺傳有限單元法反演川滇下地殼流動(dòng)對(duì)上地殼的拖曳作用. 地球物理學(xué)報(bào), 47(2): 232-239.

朱守彪, 張培震. 2009. 2008年汶川MS8.0地震發(fā)生過(guò)程的動(dòng)力學(xué)機(jī)制研究. 地球物理學(xué)報(bào), 52(2): 418-427.

(本文編輯 何燕)

Focal mechanisms and tectonic stress field in the North-South Seismic Belt of China

WANG Xiao-Shan1, Lü Jian2, 3*, XIE Zu-Jun3, LONG Feng4, ZHAO Xiao-Yan5, ZHENG Yong3

1EarthquakeAdministrationofHebeiProvince，Shijiazhuang050021，China2EarthquakeAdministrationofJiangxiProvince，Nanchang330039，China3StateKeyLaboratoryofGeodesyandEarth'sDynamics,InstituteofGeodesyandGeophysics,ChineseAcademyofSciences,Wuhan430077,China4EarthquakeAdministrationofSichuanProvince，Chengdu610041，China5EarthquakeAdministrationofYunnanProvince，Kunming650224，China

The North-South Seismic Belt is the boundary between the first-order blocks of the tectonic stress field in China. The focal mechanisms of 819 events that occurred in this belt from 1970 to 2014 are used to study its tectonic stress field, which is important for understanding the mechanism of continental earthquakes, tectonic deformation and seismic stress interaction.

According to the classification criteria of World Stress Map (WSM), focal mechanism solutions are divided into six types, in which the horizontal projections ofPaxes mark the movement direction of the block. The stress tensors and the orientations of the maximum horizontal stress of the North-South Seismic Belt were determined by the MATLAB Package for Spatial And Temporal Stress Inversion (MSATSI) using 2°×2° grids with the damping parameter 1.6 at the 95% confidence level.

The spatial distribution of different types of focal mechanisms is consistent with the nature of tectonic activity in the North-South Seismic Belt. It demonstrates that the regional tectonic stress field is controlled by tectonic movement. The horizontal projections ofPaxes show the segmentation patterns from north to south: NNE—NE—nearly EW—NW—nearly NS—NNE direction, the horizontal projections ofTaxes exhibit the most prominent characteristic of spatial differences in the Sichuan-Yunnan block and neighboring regions, a clockwise rotational pattern. The spatial distribution of stress regime and the orientation of the maximum horizontal stress has obvious regional characteristics, i.e. NE direction strike-slip type stress state in the northern section, NEE—EW—NWW direction thrusting type in the middle section, SE—SSE—NS—NNE direction strike-slip and normal faulting type in the southern section, and the stress regime in the northern part and west boundary of the Sichuan-Yunnan block is EW—SE—SSE direction normal faulting. The normal fault stress regime in the western boundary of the Sichuan-Yunnan block is consistent with the scope of the low velocity zone in mid-lower crust observed by the high resolution seismic array, the ductile flow of eastern expansion of the Tibetan Plateau and the weakness of lateral boundaries (Lijiang-Xiaojinhe fault) facilitate energy release, limiting earthquake magnitude. It also localizes the transition from normal fault earthquakes with N—S extension to normal faulting with E—W extension.

The results suggest that the lithospheric stress field of the North-South Seismic Belt is controlled jointly by the NE or NNE directed horizontal compressive stress from the Indian plate and eastward motion of the material of the Tibetan Plateau that slip towards SE along large strike-slip faults. Stress regime is almost consistent with the rupture mechanism of all types of earthquakes, verifying the reliability of the inversion. It can offer a reference for simulation of geodynamic processes and determination of the nature of fault slip.

North-South Seismic Belt; Focal mechanism; Tectonic stress field; Orientation of the maximum horizontal stress

10.6038/cjg20151122.

中國(guó)地震局地震科技星火計(jì)劃(XH12027)，大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(SKLGED2013-5-4-E)，國(guó)家自然科學(xué)基金重大研究計(jì)劃項(xiàng)目(91214201)資助.

王曉山，男，1980年生，博士研究生，主要從事數(shù)字地震學(xué)與地震危險(xiǎn)性研究.E-mail：wxs@eq-he.ac.cn

*通訊作者 呂堅(jiān)，男，1975年生，高級(jí)工程師，主要從事地震活動(dòng)性與地震危險(xiǎn)性研究.E-mail：lvjian8438@163.com

10.6038/cjg20151122

P315

2015-02-02，2015-06-29收修定稿

王曉山， 呂堅(jiān)， 謝祖軍等. 2015. 南北地震帶震源機(jī)制解與構(gòu)造應(yīng)力場(chǎng)特征.地球物理學(xué)報(bào)，58(11):4149-4162,

Wang X S, Lü J, Xie Z J, et al. 2015. Focal mechanisms and tectonic stress field in the North-South Seismic Belt of China.ChineseJ.Geophys. (in Chinese),58(11):4149-4162,doi:10.6038/cjg20151122.