劉自鳳 龍鋒 彭關(guān)靈 趙小艷

摘要：利用1965年1月—2018年11月云南地震臺網(wǎng)提供的月報目錄和1999年1月—2018年10月中小地震震源機制解，對滇西北地區(qū)主要斷裂的b值和震源機制一致性參數(shù)進(jìn)行了空間掃描，結(jié)合歷史強震空區(qū)，分析了滇西北地區(qū)的強震危險性。結(jié)果表明：①2000年1月以來程海斷裂、龍蟠—喬后斷裂及維西—喬后斷裂南段b值顯著降低，而中甸斷裂地震空區(qū)兩端、程海斷裂及鄰區(qū)則為低b值（b<07）主體區(qū)域。②震源機制一致性較好的區(qū)域分布于程海斷裂中部、龍蟠—喬后斷裂及維西—喬后斷裂南段及鄰區(qū)。分析認(rèn)為，低b值和震源機制一致性好可能是區(qū)域應(yīng)力水平增強所致，據(jù)此判斷程海斷裂、龍蟠—喬后斷裂及維西—喬后斷裂南段及鄰區(qū)可能是滇西北地區(qū)未來發(fā)生強震的潛在危險區(qū)。

關(guān)鍵詞：b值;震源機制一致性;應(yīng)力張量方差;強震危險性;滇西北

中圖分類號：P31572?? 文獻(xiàn)標(biāo)識碼：A?? 文章編號：1000-0666（2019）03-0330-08

0 引言

強震活動往往發(fā)生在活動斷裂帶上具有高應(yīng)力積累的凹凸體段或者閉鎖段（Aki，1984;Wiemer，Wyss，1997;Wyss et al，2000）。地震一旦成核，破裂將會在高應(yīng)力環(huán)境下擴展，并穿過凹凸體及其附近的非凹凸體的斷裂段（Wyss et al，2004）。古登堡-里克特所提出的震級-頻度關(guān)系（Gutenber，Richter，1944）中的b值與介質(zhì)的非均勻度、有效剪應(yīng)力等參數(shù)有關(guān)，并與應(yīng)力呈反比關(guān)系，這一現(xiàn)象早已在巖石破裂實驗（Scholz，1968）、與流體抽取有關(guān)的地震活動（Wyss，1973）以及地下礦山巖石破裂（Urbancic et al，1992）中被觀測到。特定巖體內(nèi)構(gòu)造應(yīng)力的大小與b值成反比，因此低b值可作為評判高有效剪應(yīng)力的一個指標(biāo)（Wyss et al，2000）。目前，利用b值的空間分布揭示和推斷一個區(qū)域內(nèi)活動斷裂帶不同段落的相對應(yīng)力水平，已成為分析不同斷裂段現(xiàn)今活動習(xí)性與強震（大震）潛在危險地段的重要手段之一（Wyss et al，2000;Wyss，Stefansson，2006;錢曉東，秦嘉政，2008;易桂喜等，2008，2010，2013;邵延秀等，2015;吳萍萍等，2015;毛燕等，2016;韓曉明等，2016;張廣偉，2016）。

震源機制一致性參數(shù)（陳颙，1978;刁桂苓等，2004;王俊國，刁桂苓，2005;程萬正等，2006;澤仁志瑪?shù)龋?010）在震情跟蹤工作中已經(jīng)成為判斷地震危險性的一個有效判據(jù)，基于震源機制解反演應(yīng)力場的時空變化捕捉強震前兆信息的方法也被廣泛運用。

滇西北地區(qū)處于青藏高原東南緣、川滇菱形塊體西南部，受印度板塊與歐亞板塊碰撞作用（Dewey，Sun，1988;Rowley，1996）的影響，新構(gòu)造活動顯著，發(fā)育了眾多受近南北向斷裂控制的地塹或半地塹盆地，構(gòu)成了“Z”字形的滇西北裂陷帶（吳中海等，2008;黃小巾等，2014;羅睿潔等，2015）。該區(qū)域歷史上M≥60地震頻發(fā)，且大多發(fā)生在主要斷裂及其附近區(qū)域，但M7專項工作組（2012）研究表明，在M≥60地震主體活動斷裂上仍存在3個地震空區(qū)，分別為中甸斷裂地震空區(qū)、小金河斷裂地震空區(qū)和賓川地震空區(qū)。1996年2月3日麗江M70地震后，滇西北地區(qū)M≥60地震活動較弱，僅發(fā)生2組6級地震，最大震級MS62，可見該區(qū)自1996年麗江地震后強震進(jìn)入弱活動狀態(tài)，地震空區(qū)持續(xù)發(fā)展。2013年3月以來，滇西北地區(qū)ML≥40地震顯著增強，且主要集中分布在維西—喬后斷裂、程海斷裂及其附近地區(qū)。2013年3月—2018年11月共計發(fā)生ML≥40地震20次，其中40～49級16次、50級以上4次，最大為2013年3月3日洱源55級地震。因此，探討滇西北地區(qū)的強震危險性具有重要意義。本文結(jié)合滇西北地區(qū)歷史強震空區(qū)，計算b值與震源機制一致性參數(shù)，分析該區(qū)域主要斷裂和歷史強震空區(qū)現(xiàn)今應(yīng)力狀態(tài)和強震危險性。

1 資料選取與預(yù)處理

11 地震目錄及最小完整性震級

本文利用云南地震臺網(wǎng)提供的1965年1月—2018年11月月報目錄計算滇西北地區(qū)（250°～285°N，990°～1015°）b值。為避免統(tǒng)計結(jié)果 受資料選取的影響，余震的刪除、小震群去叢和地震目錄的完整性分析非常重要。本文在利用k-k法初步去除余震的基礎(chǔ)上，人工進(jìn)行了二次核對，最終獲得0級以上地震94 148次。再采用EMR方法（基于G-R關(guān)系基礎(chǔ)上的完整性震級范圍法）計算研究區(qū)地震目錄最小完整性震級MC值。由圖1a可見，MC值隨時間呈減小趨勢，這與地震臺網(wǎng)的擴展、地震分析方法的不斷改進(jìn)和監(jiān)測能力的不斷提高有關(guān);2000年前后MC值大幅度降低，這是因為該時段云南地區(qū)臺站實現(xiàn)了從模擬記錄向數(shù)字化記錄的轉(zhuǎn)變（皇甫崗，李忠華，2010）。為了減弱MC的大幅變化和研究區(qū)內(nèi)發(fā)生的1996年2月3日麗江M70地震對b值計算結(jié)果的影響，分2個時段分別計算b值，即1965年1月—1995年12月（Ⅰ時段）、2000年1月—2018年11月（Ⅱ時段）。圖1b，c顯示，MC在Ⅰ時段為22±003，在Ⅱ時段為13±003。

12 震源機制解

本文根據(jù)2002年立項的國家重大基礎(chǔ)研究前期研究專項“新一代地震參數(shù)目錄測定”、專著（徐彥，2013）、2016年立項的地震監(jiān)測預(yù)報專項和中國地震局滇西地震預(yù)報實驗場小孔徑臺網(wǎng)給出的相關(guān)結(jié)果整理了研究區(qū)1999年1月—2018年10月絕大部分中小地震震源機制參數(shù)。為了排除序列的余震對震源機制一致性參數(shù)的影響，剔除了余震的震源機制解，最終獲得1 158次ML≥25地震的震源機制解，其中25～29級地震456次、30～39級地震576次、40～49級地震107次、50～59級地震15次、60～69級地震4次。

2 研究方法

21 b值計算

古登堡-里克特提出的震級與對應(yīng)頻度的對數(shù)線性關(guān)系為：

按照一定間隔將研究區(qū)進(jìn)行網(wǎng)格化處理，采用最大似然法計算出ML≥MC地震的各統(tǒng)計單元的b值;再根據(jù)每一單元的b值得到相應(yīng)區(qū)域的b值;最后基于應(yīng)力與b值呈反比關(guān)系（Wyss et al，2000）分析滇西北主要斷裂帶和M7專項工作組（2012）提出的歷史強震空區(qū)的地震危險性。易桂喜等（2008）研究指出，b>10的區(qū)域表明區(qū)域內(nèi)斷層具有以頻繁小震活動為主要特征的蠕滑性質(zhì)，未來發(fā)生強震的可能性相對較小，反之，則區(qū)域整體應(yīng)力水平相對較高，而Wyss等（2004）所獲得的凹凸體b<07，因此一般b<07的區(qū)域為低b值異常區(qū)。

22 震源機制一致性計算

已有研究提出多種表征震源機制一致性的參數(shù)，如單個地震震源機制解的P，T，B軸和構(gòu)造應(yīng)力場的3個正交應(yīng)力主軸在三維空間的夾角之和（刁桂苓等，2004;趙英萍等，2004）、單個地震的滑動矢量與在平均剪切應(yīng)力作用下產(chǎn)生的滑動矢量之間的夾角misfit角（Michael，1984，1987）、應(yīng)力張量方差（Michael et al，1990;Lu et al，1997）和體波譜振幅相關(guān)系數(shù)（朱航等，2006;段云歌等，2018）等。本文采用Michael等（1990）反演研究區(qū)的應(yīng)力張量方差，并用此表征震源機制一致性。Michael等（1990）的方法是按照剪應(yīng)力方向與斷層滑動方向相一致的原則提出的將非線性問題轉(zhuǎn)化為線性反演應(yīng)力張量的算法。計算過程中，由于無法區(qū)分震源機制解的斷層面和輔助面，因此在反演應(yīng)力張量時將2個節(jié)面等同看待，通過節(jié)面與假定應(yīng)力張量的擬合函數(shù)關(guān)系選取其中擬合誤差較小的節(jié)面作為震源機制斷層面?；谘芯繀^(qū)每個網(wǎng)格節(jié)點及其周圍一定范圍內(nèi)的多個震源機制解，用每個震源機制解給激發(fā)斷層運動的應(yīng)力張量1個限制，通過若干個解的逼近則可以得到作用于該節(jié)點的應(yīng)力張量解的方向，在應(yīng)力張量結(jié)果的基礎(chǔ)上求取應(yīng)力張量方差和misfit角。

Michael等（1990）認(rèn)為misfit角大小或應(yīng)力張量方差高低均能表征震源機制一致性程度，其低值反映震源機制一致性好，應(yīng)力場是均勻的;高值則反映震源機制一致性差，應(yīng)力場是非均勻的。Lu等（1997）研究指出當(dāng)應(yīng)力張量方差大于02時，說明該區(qū)域的應(yīng)力場在時空分布上都是非均勻的，震源機制一致性程度低，可解釋為研究區(qū)區(qū)域應(yīng)力水平不高;當(dāng)應(yīng)力張量方差小于01時，說明震源機制一致性程度高，且可以用統(tǒng)一的應(yīng)力張量解釋在該區(qū)域觀測到的震源機制解，也可以理解為研究區(qū)域統(tǒng)一受控于較高的區(qū)域應(yīng)力場。通過反演震源機制一致性參數(shù)的空間分布特征，可確定相對高應(yīng)力區(qū)，再對相對高應(yīng)力區(qū)的一致性參數(shù)進(jìn)行時間掃描，根據(jù)變化形態(tài)可分析區(qū)域內(nèi)的強震危險程度和緊迫性。綜上，misfit角或應(yīng)力張量方差的時空變化特征可以用來分析研究區(qū)未來強震的可能發(fā)生地點和時間。

3 結(jié)果分析

31 b值空間變化特征

本文研究區(qū)構(gòu)造復(fù)雜，展布有龍蟠—喬后斷裂、維西—喬后斷裂、中甸斷裂、小金河斷裂、程海斷裂和紅河斷裂等（圖2a），歷史強震幾乎都發(fā)生在這些斷裂帶上，而賓川地震空區(qū)、小金河斷裂地震空區(qū)和中甸斷裂地震空區(qū)分別位于SN向的程海斷裂中南段、NE向的小金河斷裂帶南段和NW向中甸斷裂帶上。因此將斷裂帶附近的中小地震密集區(qū)以003°的間隔進(jìn)行網(wǎng)格化，以每個網(wǎng)格節(jié)點為圓心，在半徑20 km的范圍內(nèi)進(jìn)行地震統(tǒng)計并估算b值。計算條件要求每個網(wǎng)格的地震數(shù)量不少于30個，參與擬合的有效震級分檔數(shù)不低于5檔，不滿足條件的網(wǎng)格在圖中用空白顯示，即沒有b值計算結(jié)果。

分別以22，13為起始震級掃描計算研究區(qū)Ⅰ，Ⅱ時段的b值，空間分布如圖2b，c所示。由圖2b明顯可見，Ⅰ時段的b值在中甸斷裂地震空區(qū)、小金河斷裂地震空區(qū)和程海斷裂中段普遍低于10，表明該時段這些區(qū)域應(yīng)力水平相對于寧蒗—鹽源（b>10）和彌渡—祥云地區(qū)（b>10）要高，而在中甸斷裂地震空區(qū)東南端和永勝—鶴慶一帶的低b值異常尤為顯著（b<07）;圖2c所示，Ⅱ時段b值在中甸斷裂地震空區(qū)、程海斷裂中北段和劍川—洱源—漾濞地區(qū)低于10，而中甸斷裂地震空區(qū)南端、賓川地震空區(qū)北端與永勝地區(qū)、程海斷裂北段與小金河斷裂地震空區(qū)東南角交匯區(qū)域b值低于07，幾個局部地區(qū)b值甚至低于06。

對比Ⅱ時段與Ⅰ時段b值空間分布可看出，中甸地區(qū)低b值（b<07）區(qū)域有所擴展，永勝—鶴慶地區(qū)低b值（b<07）區(qū)域明顯向程海斷裂方向收縮，在程海斷裂北段新增一個低b值（b=05～07）區(qū)域。而麗江地區(qū)b值從Ⅰ時段低于10顯著升高為11～15，賓川地震空區(qū)南部的b???值也從Ⅰ時段的08～10顯著升高為13以上。因此，統(tǒng)一以22（Ⅰ時段的最小完整性震級）為起始震級計算Ⅱ時段相對于Ⅰ時段的b值變化量Δb（圖2d）。圖2d顯示，Ⅱ時段寧蒗南部—永勝—賓川一帶和鶴慶—劍川—洱源—漾濞一帶b值降低（Δb<0）。由此可見，2000年以來，b值顯著降低的區(qū)域幾乎沿著程海斷裂、龍蟠—喬后斷裂和維西—喬后斷裂南段展布，這與2000年以來滇西地區(qū)4級以上地震空間分布格局相吻合。

32 震源機制一致性參數(shù)空間變化特征

本文收集的1999年1月—2018年10月滇西北地區(qū)的1 158次ML≥25地震震源機制空間分布如圖3a所示，為了便于與b值計算結(jié)果進(jìn)行對比，震源機制一致性參數(shù)掃描區(qū)域與b值掃描區(qū)域一致。將擬掃描區(qū)域按照005°×005°進(jìn)行網(wǎng)格化，選取每個網(wǎng)格節(jié)點及其周圍至少15個地震的震源機制解反演每個節(jié)點的應(yīng)力張量方差。研究區(qū)應(yīng)力張量方差空間分布特征如圖3b所示，從圖中可見，永勝—鶴慶—劍川—洱源—賓川一帶、小金河斷裂

4 結(jié)論與討論

本文利用1965年1月—2018年11月云南地震臺網(wǎng)提供的月報目錄和1999年1月—2018年10月中小地震震源機制解，分析了滇西北主要斷裂帶和3個歷史地震空區(qū)的強震危險性，主要得出以下結(jié)論：

（1）2000年1月以來，研究區(qū)的低b值（<07）主體區(qū)域為中甸斷裂地震空區(qū)南端、賓川地震空區(qū)北端與永勝地區(qū)、程海斷裂中北段與小金河斷裂地震空區(qū)交匯地區(qū)，上述區(qū)域應(yīng)力水平較高，具備發(fā)生強震的條件，而麗江地區(qū)b值較高，發(fā)生強震的可能性不大。

（2）對比Ⅱ時段與Ⅰ時段b值空間分布可知：中甸地區(qū)低b值區(qū)域有所擴展，這可能與Ⅱ時段MC降低，該地區(qū)參與計算的地震相對增多有關(guān);永勝—鶴慶地區(qū)低b值區(qū)域往程海斷裂方向大幅收縮，且在程海斷裂北段新增一個低b值區(qū)域，而在其南段b值升高且空間上明顯出現(xiàn)擴展，這表明程海斷裂中北段的強震危險性在不斷增加;麗江中東部和鶴慶北部地區(qū)b值顯著升高，這可能與1996年麗江7級地震能量強釋放有一定關(guān)系。Δb計算結(jié)果顯示，Ⅱ時段永勝—賓川一帶和鶴慶—劍川—洱源—漾濞一帶b值相對于Ⅰ時段降低。由此可見2000年以來，b值顯著降低的區(qū)域幾乎沿著程海斷裂、龍蟠—喬后斷裂和維西—喬后斷裂南段分布，這與2000年以來滇西地區(qū)4級以上地震空間分布格局基本一致。

（3）應(yīng)力張量方差低值主體區(qū)域為永勝仁和地區(qū)、賓川—鶴慶—劍川—洱源地區(qū)，其值為010～015。這表明上述地區(qū)應(yīng)力是相對均勻的，震源機制一致性較好，應(yīng)力水平較高，發(fā)生強震的可能性較大。

綜上，滇西北地區(qū)b值和應(yīng)力張量方差的空間掃描結(jié)果顯示程海斷裂、龍蟠—喬后斷裂和維西—喬后斷裂南段及鄰區(qū)低值變化特征相對顯著，這可能是區(qū)域應(yīng)力水平增強的表現(xiàn)。因此分析認(rèn)為，未來滇西北地區(qū)的強震發(fā)生在程海斷裂、龍蟠—喬后斷裂和維西—喬后斷裂南段及其鄰區(qū)的危險性較大。

根據(jù)b值和震源機制一致性參數(shù)結(jié)果可判別出強震的相對危險地點，但在時間尺度上屬于中-長期類，如果要在時間判斷方面給出更明確的結(jié)論，需要結(jié)合其它地震學(xué)參數(shù)和地球物理觀測資料開展進(jìn)一步的研究后再進(jìn)行綜合判定。

參考文獻(xiàn)：

陳陽，呂悅軍，謝卓娟，等2013地震活動性參數(shù)b值的研究[C]北京：地殼構(gòu)造與地殼應(yīng)力所，38-47

陳颙1978用震源機制一致性作為描述地震活動性的新參數(shù)[J].地球物理學(xué)報，21（2）：142-159

程萬正，阮祥，張永久2006川滇次級地塊震源機制類型與一致性參數(shù)[J].地震學(xué)報，28（6）：561-573

刁桂苓，趙英萍，啜永清，等2004大同晚期強余震前震源機制解的一致性特征[J].內(nèi)陸地震，18（4）：202-206

段云歌，陳天長，蘇金蓉，等20182013年蘆山MS70地震余震序列震源機制一致性的時空變化[J].地震研究，41（2）：252-257

韓曉明，張文韜，王樹波2016河套地震帶的b值時空變化特征分析[J].中國地震，32（3）：522-532

皇甫崗，李忠華201020世紀(jì)云南地區(qū)地震記錄完整性評價[J].地震研究，33（1）：1-6

黃小巾，吳中海，李家存，等2014滇西北裂陷帶的構(gòu)造地貌特征與第四紀(jì)構(gòu)造活動性[J].地質(zhì)通報，33（4）：578-593

李永莉，蔡靜觀，曹刻2002云南地區(qū)強震活動過程中的調(diào)制比b值[J].地震研究，25（1）：25-30

羅睿潔，吳中海，黃小龍，等2015滇西北賓川地區(qū)主要活動斷裂及其活動構(gòu)造體系[J].地質(zhì)通報，34（1）：155-170

毛燕，劉自鳳，葉建慶2016小江斷裂帶強震危險性分析[J].地震研究，39（2）：213-217

錢曉東，秦嘉政2008小江斷裂帶及周邊地區(qū)強震危險性分析[J]地震研究，31（4）：354-361

邵延秀，袁道陽，梁明劍2015滇西南地區(qū)龍陵—瀾滄斷裂帶地震危險性評價[J].地震學(xué)報，37（6）：1011-1023

王輝，曹建玲，荊燕，等2012川滇地區(qū)強震活動前b值的時空分布特征[J].地震地質(zhì)，34（3）：531-543

王俊國，刁桂苓2005千島島弧大震前哈佛大學(xué)矩心矩張量（CMT）解一致性的預(yù)測意義[J].地震學(xué)報，27（2）：178-183

吳萍萍，李振，葉慶東2015郯廬斷裂帶南段及鄰區(qū)地震b值的空間分布特征[J].中國地震，31（2）：372-381

吳中海，張永雙，胡道功，等2008滇西北哈巴玉龍雪山東麓斷裂的晚第四紀(jì)正斷層作用及其動力學(xué)機制探討[J].中國科學(xué)：地球科學(xué)，38（11）：1361-1375

徐彥2013云南地區(qū)ML30以上中小地震震源機制解匯編[M].昆明：云南科技出版社，168-186

易桂喜，聞學(xué)澤，蘇有錦2008川滇活動地塊東邊界強震危險性研究[J].地球物理學(xué)報，51（6）：1719-1725

易桂喜，聞學(xué)澤，辛華等2013龍門山斷裂帶南段應(yīng)力狀態(tài)與強震危險性研究[J].地球物理學(xué)報，56（4）：1112-1120

易桂喜，聞學(xué)澤，張致偉2010川南馬邊地區(qū)強震危險性分析[J].地震地質(zhì)，32（2）：282-293

澤仁志瑪，刁桂苓，李志雄，等2010大震前顯示的地震震源機制趨于一致的變化[J].地震，30（1）：108-114

張廣偉2016云南地區(qū)地震的重新定位及b值研究[J].中國地震，32（1）：54-62

趙英萍，刁桂苓，高景春，等2004張北強余震前震源機制解的一致性特征[J].華北地震學(xué)報，22（1）：1-4

朱航，劉杰，陳天長2006采用體波譜振幅相關(guān)系數(shù)方法研究地震序列的震源機制變化過程[J].地震，26（2）：1-11

M7專項工作組2012中國大陸大地震中-長期危險性研究[M].北京：地震出版社，234-235