易桂喜，聞學(xué)澤，辛 華，喬慧珍，王思維，宮 悅

1 四川省地震局，成都 610041

2 成都理工大學(xué)地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，成都 610059

3 中國(guó)地震局地震預(yù)測(cè)研究所，北京 100036

1 引 言

2008年5月12日汶川M8．0級(jí)大地震發(fā)生在青藏高原東緣龍門(mén)山斷裂帶中－北段，破裂沿中－北段呈北東向單側(cè)擴(kuò)展［1－3］，龍門(mén)山斷裂帶南段未參與汶川地震的破裂與余震活動(dòng)［4］（見(jiàn)圖1）．南段無(wú)6．5級(jí)以上歷史地震記錄，近代僅發(fā)生了1941 年6 月12日寶興與康定間M6．0、1970年2 月24 日大邑M6．2 兩次6 級(jí)以上地震，屬于長(zhǎng)期缺少大震級(jí)（M≥7．0）地震 事 件 的 破 裂 空 區(qū)［4－5］．1970 年 大 邑M6．2級(jí)地震后，南段未發(fā)生5級(jí)以上地震．汶川地震導(dǎo)致周邊斷層的應(yīng)力狀態(tài)發(fā)生不同程度的變化，龍門(mén)山斷裂帶南段庫(kù)侖破裂應(yīng)力增強(qiáng)［6－13］，有利于強(qiáng)震的發(fā)生．然而，汶川地震至今，南段及近鄰發(fā)生的最大震級(jí)地震為2012年10月19日天全—蘆山間M4．2（或ML4．6），震中位于龍門(mén)山斷裂帶南段前陸地帶，南段處于較低水平的地震活動(dòng)狀態(tài)．屬于長(zhǎng)期的大地震空區(qū)、近年來(lái)地震活動(dòng)水平偏低的龍門(mén)山斷裂帶南段目前的應(yīng)力狀態(tài)及其強(qiáng)震危險(xiǎn)性成為重點(diǎn)關(guān)注的問(wèn)題之一．

圖1 龍門(mén)山及鄰區(qū)活動(dòng)構(gòu)造圖［4］黑色粗虛線以南為龍門(mén)山斷裂帶南段．彩色圓為汶川5級(jí)以上余震震中．

強(qiáng)震通常發(fā)生于活動(dòng)斷裂帶上具有高應(yīng)力積累的凹凸體［14］或 閉 鎖 斷 裂 段［15－16］．實(shí) 驗(yàn) 室 巖 石 破 裂實(shí)驗(yàn)［17］、與流體抽取有關(guān)的地震活動(dòng)研究［18］以及礦山巖石破裂觀測(cè)［19］等均表明，應(yīng)力與古登堡－里克特震級(jí)（M）－頻度（N）關(guān)系式［20］

中的b值呈反比，b值越低，應(yīng)力水平越高．式中a、b為回歸常數(shù)．b值作為估計(jì)斷裂帶相對(duì)應(yīng)力水平的主要地震活動(dòng)性參數(shù)之一，已廣泛用于中－長(zhǎng)期地震危 險(xiǎn) 性 評(píng) 價(jià)［15－16，21－30］．近年 來(lái)，地震 視 應(yīng) 力 也 逐漸用于斷裂帶應(yīng)力狀態(tài)分析與地震趨勢(shì)判定［31－36］．本文利用b值與地震視應(yīng)力值，評(píng)價(jià)龍門(mén)山斷裂帶南段當(dāng)前應(yīng)力狀態(tài)和強(qiáng)震危險(xiǎn)性．

2 資料和方法

2．1 b值

利用四川區(qū)域地震臺(tái)網(wǎng)1977年1月至2012年12月的地震目錄資料計(jì)算b值．龍門(mén)山斷裂帶南段位于該臺(tái)網(wǎng)監(jiān)測(cè)能力相對(duì)較高的區(qū)域［37］．ML≥2．0級(jí)地震的震級(jí)－頻度線性擬合（圖2）相關(guān)系數(shù)R＝0．999，直線斜率b＝1．02±0．007，與全球范圍的區(qū)域平均b值一致．因此，在計(jì)算b值時(shí)，取最小完整性震級(jí)MC為2．0．圖3給出了1977—2012年龍門(mén)山斷裂帶南段ML≥2．0級(jí)地震事件的震中分布，圖中用黑色和綠色區(qū)別汶川地震前、后的地震事件．

圖2 龍門(mén)山斷裂帶南段震級(jí)－頻度分布圖（1977－01－01—2012－12－31）

沿?cái)嗔褞值計(jì)算方法：沿龍門(mén)山斷裂帶南段，以0．1°×0．1°間距進(jìn)行網(wǎng)格化，挑選出以每個(gè)網(wǎng)格節(jié)點(diǎn)為圓心、半經(jīng)為r的圓形統(tǒng)計(jì)單元內(nèi)的地震．確定統(tǒng)計(jì)單元內(nèi)能滿足整個(gè)研究時(shí)段的最小完整性震級(jí)MC；然后，利用最小二乘法由各單元震級(jí)M≥MC的地震資料計(jì)算出（1）式中的b值，將其作為相應(yīng)單元中心點(diǎn)（即網(wǎng)格節(jié)點(diǎn)）的計(jì)算值，進(jìn)而獲得b值空間分布．計(jì)算時(shí)震級(jí)分檔間隔取0．1，每個(gè)統(tǒng)計(jì)單元內(nèi)的地震樣本數(shù)不少于50，參與擬合的有效震級(jí)分檔數(shù)不低于5檔．統(tǒng)計(jì)單元的半徑r值通常取20km；對(duì)于地震分布較稀疏的局部區(qū)域，擴(kuò)大r值，最大不超過(guò)40km．

2．2 地震視應(yīng)力

地震視應(yīng)力σapp定義為［38］：

式中，ES為地震波輻射能量；M0為地震矩；μ為震源區(qū)介質(zhì)剪切模量，通常取μ＝3．0×104MPa［38－40］．ES與M0之比表示單位地震矩輻射出的地震波能量．地震視應(yīng)力為震源區(qū)絕對(duì)應(yīng)力水平的下限估計(jì)［31，40］，σapp值越高，震源區(qū)應(yīng)力水平越高［31－36］．

在近震源情況下，基于Brune圓盤(pán)模型［41］，忽略非彈性衰減，震源譜可表示為［42］：

根據(jù)地震波形三分向速度記錄，可得合成位移譜：

在低頻段f1—f2，Ω（f）≈Ω0，則合成位移譜平方的積分SD 為：

獲得地震波零頻極限Ω0：

進(jìn)而求得地震矩M0［42］、考慮低頻和高頻補(bǔ)償后 的地震波輻射能量ES［43］：

式中，f為頻率；fc為拐角頻率，求取方法參見(jiàn)相關(guān)文獻(xiàn)［34］；V（f）為速度譜，下標(biāo)Z、NS、EW 分別表示垂直、南北、東西向分量；d為震源距；β為S波速，取β＝3．5km／s；ρ為介質(zhì)密度，取ρ＝2．71g／cm3；Rθφ為輻射圖型因子，取S 波的均方根輻射圖型因子

本文選取四川區(qū)域地震臺(tái)網(wǎng)震中距≤200km的臺(tái)站記錄的具有較高信噪比的寬頻帶數(shù)字化地震波形資料，利用S波1．0～20Hz頻帶范圍的波形計(jì)算震源譜，由（2）式逐臺(tái)求取地震視應(yīng)力σapp，取各臺(tái)計(jì)算結(jié)果的平均值作為相應(yīng)地震的視應(yīng)力值．為消除個(gè)別臺(tái)站異常高值對(duì)結(jié)果的影響，采用Archuleta等［45］的方法求取平均值ˉx與標(biāo)準(zhǔn)差Δx：

其中，N為計(jì)算使用的臺(tái)站數(shù)，xi為第i個(gè)臺(tái)站的計(jì)算值．

地震視應(yīng)力與震級(jí)的關(guān)系已有不少研究．例如，李艷娥等先后計(jì)算了2007年寧洱M6．4級(jí)地震前滇西南地區(qū)98 次ML≥2．5 級(jí)地震的視應(yīng)力、2008年汶川M8．0級(jí)地震前四川地區(qū)439次ML≥3．0級(jí)地震的視應(yīng)力，發(fā)現(xiàn)：滇西南地區(qū)ML2．5～5．0級(jí)地震的視應(yīng)力與震級(jí)具有較好的線性關(guān)系，相關(guān)系數(shù)0．77［34］，四川地區(qū)視應(yīng)力σapp的對(duì)數(shù)與震級(jí)ML呈正比，相關(guān)系數(shù)R＝0．634［35］，數(shù)據(jù)離散性高于滇西地區(qū)的．王瓊等［36］對(duì)1999年岫巖M5．4級(jí)地震序列的視應(yīng)力研究結(jié)果顯示，視應(yīng)力與震級(jí)關(guān)系具有顯著的分段性：ML2．0～4．5級(jí)地震的視應(yīng)力與震級(jí)具有較好的線性關(guān)系，對(duì)于ML≥4．5級(jí)地震，兩者非線性相關(guān)．上述研究也表明，地震視應(yīng)力與震級(jí)的相關(guān)性存在較顯著的地區(qū)差異．基于中小地震的視應(yīng)力具有隨震級(jí)增大而升高的特征［34－36］，本文在分析比較不同地震的視應(yīng)力值與區(qū)域應(yīng)力的關(guān)系時(shí)，采用如下原則：震級(jí)相當(dāng)時(shí)，視應(yīng)力值相對(duì)高的震源區(qū)應(yīng)力水平高；視應(yīng)力值相等時(shí)，震級(jí)較低的震源區(qū)應(yīng)力水平高．

3 結(jié) 果

圖4為龍門(mén)山斷裂帶南段b值圖像，除個(gè)別統(tǒng)計(jì)單元b值標(biāo)準(zhǔn)誤差為0．02外，其余統(tǒng)計(jì)單元b值標(biāo)準(zhǔn)誤差≤0．01．其中，圖4a為當(dāng)前b值分布，利用1977年1月1日至2012年12月31日ML≥2．0級(jí)地震獲得．圖4b為Δb值分布圖，是當(dāng)前b值與利用1977年1月1日至2008年5月11日時(shí)段ML≥2．0級(jí)地震獲得的汶川地震前的b值兩者之間的差異，用于推測(cè)大地震后斷裂帶應(yīng)力水平的變化．Δb正值表示b值升高，應(yīng)力水平下降；負(fù)值則說(shuō)明b值降低，應(yīng)力水平升高．

結(jié)果表明，當(dāng)前龍門(mén)山斷裂帶南段整體平均b值約1．02，與由1977年1月1日至2008年5月11日資料獲得的汶川地震前南段的平均值1．03相近．

龍門(mén)山斷裂帶南段及鄰區(qū)b值橫向差異顯著（圖4a）．瀘定附近、西側(cè)的鮮水河斷裂帶康定以南段以明顯高于區(qū)域平均值（1．02）的高b值為特征，天全北部區(qū)域b值高于1．3，反映龍門(mén)山斷裂帶最西南瀘定、天全北部及西側(cè)的鮮水河斷裂帶康定以南斷裂段目前應(yīng)力水平不高．天全與寶興之間及寶興以東，大部分區(qū)域b值＜0．9，低于區(qū)域平均值，表明天全—寶興及其以東目前的應(yīng)力水平高于最西南的瀘定及以西地區(qū)．

南段當(dāng)前b值與汶川地震前b值之間的差異Δb圖像（圖4b）顯示，2008年汶川8．0級(jí)地震后，瀘定、天全—蘆山以及寶興北部地區(qū)b值有較明顯降低，瀘定附近區(qū)域b值下降尤為突出，局部b值最大降幅達(dá)0．39，由汶川地震前的1．66下降到目前的1．27，下降幅度超過(guò)23%，顯示該區(qū)域在汶川地震后應(yīng)力水平有明顯上升．相反，靠近汶川地震主破裂區(qū)南端的大邑地區(qū)、南段西側(cè)的鮮水河斷裂帶康定以南段b值升高，應(yīng)力水平下降．

2007—2012年，龍門(mén)山斷裂帶南段及近鄰區(qū)域共發(fā)生ML≥3．8級(jí)地震14次，其震中依據(jù)發(fā)震時(shí)間先后由小到大編號(hào)繪于圖4中．本文計(jì)算了其中13次地震的視應(yīng)力σapp，1次事件（圖4中4號(hào)地震）因波形重疊未給出計(jì)算結(jié)果．13 次地震的視應(yīng)力值在0．12～0．82 MPa之間（圖5）．為便于與b值結(jié)果進(jìn)行對(duì)比分析，上述地震的震級(jí)與視應(yīng)力值一并標(biāo)示于圖4a中．

視應(yīng)力計(jì)算結(jié)果表明，汶川地震后，發(fā)生在龍門(mén)山斷裂帶南段的ML≥3．8級(jí)地震的視應(yīng)力水平高于其西側(cè)鮮水河斷裂帶康定段，如，2008年6 月26日大邑ML4．2、2009年9月19日寶興ML4．1級(jí)地震的視應(yīng)力值（圖4a與圖5中6、8號(hào)）明顯高于2012年3月5日康定ML4．5級(jí)地震的視應(yīng)力值（圖4a與圖5中12號(hào)）．汶川地震后，原地應(yīng)力測(cè)量結(jié)果顯示，龍門(mén)山斷裂帶南段寶興測(cè)點(diǎn)的水平最大主應(yīng)力為9．8 MPa，遠(yuǎn)高于鮮水河斷裂帶上康定測(cè)點(diǎn)的水平最大主應(yīng)力2．6 MPa［46］，反映龍門(mén)山斷裂帶南段應(yīng)力水平高于鮮水河斷裂帶康定段，與本文結(jié)果相吻合．2012年3月5日康定ML4．5級(jí)地震的視應(yīng)力值（圖4a與圖5中12 號(hào)）遠(yuǎn)低于2008 年2月16日瀘定ML4．4和27日康定ML5．0級(jí)地震的視應(yīng)力值（圖4a與圖5中2、3號(hào)），表明汶川地震后鮮水河斷裂帶康定以南段應(yīng)力水平下降．發(fā)生在龍門(mén)山斷裂帶南段南側(cè)四川盆地及盆地邊緣的2007年7月31日滎經(jīng)ML4．3、2011年7月25日雅安—名山間ML4．0、10月27日洪雅ML4．1、2012年10月19日天全—蘆山間ML4．6 級(jí)地震的視應(yīng)力值（圖4a與圖5中1、10、11和14號(hào)）均處于擬合趨勢(shì)線上方，明顯高于龍門(mén)山斷裂帶南段大邑地震的視應(yīng)力值（圖4a與圖5中6、9號(hào)），說(shuō)明四川盆地西北部的背景應(yīng)力水平高于龍門(mén)山斷裂帶南段東端．

圖5 13次ML≥3．8級(jí)地震的視應(yīng)力與震級(jí)分布地震序號(hào)同圖4．

4 討 論

靜態(tài)庫(kù)侖破裂應(yīng)力增加0．01 MPa可導(dǎo)致斷層破裂［47－48］，使其地震活動(dòng)產(chǎn)生顯著變化．Parsons等［6］計(jì)算了汶川地震周邊斷層的靜態(tài)庫(kù)侖應(yīng)力，2008年汶川地震加載到龍門(mén)山斷裂帶南段的靜態(tài)庫(kù)侖破裂應(yīng)力值達(dá)0．1 MPa；Toda等［7］的研究結(jié)果顯示，汶川地震使龍門(mén)山斷裂帶彭－灌斷裂南西段的靜態(tài)庫(kù)侖破裂應(yīng)力增加0．419 MPa；單斌等［9］計(jì)算出汶川地震使龍門(mén)山斷裂帶南段庫(kù)侖破裂應(yīng)力增加0．005～0．015 MPa，認(rèn)為南段存在發(fā)生強(qiáng)余震的可能；萬(wàn)永革等［10］基于GPS 和InSAR 資料反演的 汶川地震破裂模型，計(jì)算出汶川地震加載到龍門(mén)山斷裂帶南段靜態(tài)庫(kù)侖破裂應(yīng)力最大值達(dá)0．01 MPa；Xie等［11］的研究結(jié)果表明，汶川地震使北川—映秀斷裂和彭縣—灌縣斷裂上的應(yīng)力明顯增加，前者的應(yīng)力增加量為0．004～0．045 MPa，后者為0～0．05 MPa；華衛(wèi)等［49］的研究則顯示，寶興至映秀間庫(kù)侖應(yīng)力增加量與多次發(fā)生6級(jí)強(qiáng)余震的北川至青川段相同，均達(dá)0．05 MPa．上述研究表明，汶川8．0級(jí)地震加載到龍門(mén)山斷裂帶南段的靜態(tài)庫(kù)侖破裂應(yīng)力值已達(dá)到甚至遠(yuǎn)高于靜態(tài)應(yīng)力觸發(fā)閾值0．01 MPa．然而，與1999 年8月17日土耳其伊茲米特Mw7．4 級(jí)地震觸 發(fā)1999 年11 月12 日 迪 茲 杰Mw7．2 級(jí) 地震［50－51］、2004年12 月26日 蘇 門(mén) 答 臘Mw9．2 級(jí) 地震觸發(fā)2005年3月28日Mw8．7級(jí)地震［52］不同，汶川地震沒(méi)有在龍門(mén)山斷裂帶南段快速觸發(fā)出與汶川地震震級(jí)相似大小的地震．圖3顯示，汶川地震后，龍門(mén)山斷裂帶南段地震活動(dòng)的空間分布特征與汶川地震前基本一致，表明汶川地震的發(fā)生未明顯改變南段地震活動(dòng)的空間格局．可見(jiàn)，庫(kù)侖應(yīng)力增強(qiáng)并非引發(fā)強(qiáng)震的決定性因素，是否發(fā)生強(qiáng)震取決于增強(qiáng)區(qū)的應(yīng)力積累水平及是否接近破裂極限．

靠近汶川余震區(qū)南端的大邑附近區(qū)域b值（圖4a）低于區(qū)域平均值，但未達(dá)到異常低值（b＜0．7），汶川地震后該區(qū)域b值升高（圖4b），應(yīng)力水平降低，汶川地震后發(fā)生在該區(qū)域的兩次地震（圖4a與圖5中6和9號(hào)）的視應(yīng)力不高，反映其長(zhǎng)期應(yīng)力積累水平尚未達(dá)到發(fā)生強(qiáng)震的應(yīng)力極限，這或許是汶川8．0級(jí)地震的余震帶向南西止于大邑的重要原因之一．該區(qū)域僅在2008 年6 月28 日發(fā)生了1 次ML4．2級(jí)地震（圖4a編號(hào)6），之后再無(wú)4級(jí)地震活動(dòng)．大邑附近區(qū)域的b值升高，也說(shuō)明汶川地震后并非整個(gè)龍門(mén)山斷裂帶南段都處于應(yīng)力增強(qiáng)區(qū)．

天全—蘆山、瀘定以及寶興北部等區(qū)域b值降低（圖4b），應(yīng)力水平升高，但整個(gè)南段目前沒(méi)有出現(xiàn)如同2008年汶川8．0級(jí)地震前龍門(mén)山斷裂帶中北段那樣的異常低b值［4，27］．原地應(yīng)力測(cè)量結(jié)果顯示，汶川地震前，龍門(mén)山斷裂帶北段青川測(cè)點(diǎn)的水平最大主應(yīng)力為21～22 MPa［53］，汶川地震的最大余震即發(fā)生在青川，汶川地震后，青川測(cè)點(diǎn)復(fù)測(cè)的水平最大主應(yīng)力值仍高達(dá)15～16 MPa［53］，遠(yuǎn)高于寶興測(cè)點(diǎn)的9．8 MPa［46］．因此，即使龍門(mén)山斷裂帶南段位于汶川地震引起的庫(kù)侖應(yīng)力增強(qiáng)區(qū)［6－13，49］，其當(dāng)前的應(yīng)力積累水平也不具備發(fā)生類似汶川8．0級(jí)地震的應(yīng)力條件，這可能是汶川地震未能在南段快速觸發(fā)出相似震級(jí)大震的主要原因．趙靜等［54］利用1999—2007期GPS水平速度場(chǎng)觀測(cè)資料反演了汶川地震前龍門(mén)山斷裂帶的閉鎖程度，認(rèn)為南段閉鎖程度低于中－北段．南段短期內(nèi)不具備發(fā)生類似汶川8．0級(jí)大震的危險(xiǎn)性，但天全至寶興段低于區(qū)域平均值的b值以及相對(duì)較高的地震視應(yīng)力值（圖4a）反映該斷裂段已有一定的應(yīng)力積累，具備發(fā)生中－強(qiáng)地震的應(yīng)力條件．

龍門(mén)山斷裂帶南端西側(cè)的鮮水河斷裂帶康定以南段屬于1786年73／4級(jí)地震的破裂段［55］（圖1）．正的Δb表明汶川地震后該段b值升高，汶川地震后發(fā)生在該段的地震的視應(yīng)力值遠(yuǎn)低于汶川地震前，均揭示該段在汶川地震后應(yīng)力水平下降．該段當(dāng)前的高b值、低地震視應(yīng)力值以及相對(duì)較低的原地應(yīng)力測(cè)量值［46］反映該斷裂段在經(jīng)歷上一次強(qiáng)震破裂后斷面尚未耦合，應(yīng)力水平較低，短期內(nèi)發(fā)生強(qiáng)震的可能性較?。?/p>

近年發(fā)生在龍門(mén)山斷裂帶南段南側(cè)四川盆地西北部的幾次ML≥3．8級(jí)地震（圖4a與圖5中1、10、11、14號(hào)）的視應(yīng)力值明顯偏高，尤其是2012年10月19日天全—蘆山間ML4．6級(jí)地震的視應(yīng)力值異常高（圖4a與圖5中14號(hào)），反映毗鄰龍門(mén)山斷裂帶南段的四川盆地區(qū)域應(yīng)力水平相對(duì)較高．Parsons等［6］的庫(kù)侖應(yīng)力計(jì)算結(jié)果顯示，2008年汶川地震使位于龍門(mén)山斷裂帶南側(cè)的四川盆地區(qū)域應(yīng)力顯著升高，這或許是汶川地震后四川盆地內(nèi)4級(jí)地震活動(dòng)水平明顯高于汶川地震前的重要原因，四川盆地具有發(fā)生中－強(qiáng)地震的危險(xiǎn)性．2013年2月19日在四川盆地歷史少震與弱震區(qū)三臺(tái)縣發(fā)生M4．7級(jí)地震，支持本文的研究結(jié)果．

5 結(jié) 論

本文利用1977年1月至2012年12月的地震資料，通過(guò)精細(xì)b值空間掃描，對(duì)比汶川地震前、后b值的空間分布差異，結(jié)合中小地震視應(yīng)力值，分析了龍門(mén)山斷裂帶南段及附近區(qū)域目前的應(yīng)力狀態(tài)，進(jìn)而分析其強(qiáng)震危險(xiǎn)性，并初步探討了南段未參與汶川地震破裂與余震活動(dòng)的原因．

2008年汶川地震后，龍門(mén)山斷裂帶南段僅部分區(qū)域應(yīng)力增強(qiáng)．天全—蘆山、瀘定以及寶興北部等區(qū)域b值下降、應(yīng)力水平明顯升高．靠近汶川余震區(qū)南端的大邑附近區(qū)域b值上升、應(yīng)力水平下降，這或許是汶川余震活動(dòng)止于大邑的重要原因．

龍門(mén)山斷裂帶南段目前整體平均b值約1．02．盡管天全—蘆山、瀘定、以及寶興北部等區(qū)域在汶川地震后b值下降，但南段目前仍未出現(xiàn)類似汶川8．0級(jí)地震前中－北段那樣的異常低b值，南段整體的當(dāng)前應(yīng)力積累水平低于汶川地震前中－北段應(yīng)力水平，可能是汶川地震沒(méi)能在南段快速觸發(fā)強(qiáng)震的重要原因．短期內(nèi)南段發(fā)生類似汶川8．0級(jí)大震的可能性不大，但天全至寶興段已具有一定的應(yīng)力積累，存在發(fā)生中－強(qiáng)地震的危險(xiǎn)性．南段西側(cè)的鮮水河斷裂帶康定以南段應(yīng)力水平低，短期內(nèi)發(fā)生強(qiáng)震的可能性較?。?/p>

龍門(mén)山斷裂帶南側(cè)四川盆地西北部地區(qū)的地震視應(yīng)力水平明顯高于龍門(mén)山斷裂帶南段東端及西側(cè)的鮮水河斷裂帶康定以南段，因此，對(duì)四川盆地的中－強(qiáng)地震危險(xiǎn)性應(yīng)引起足夠重視．

由于震中距、方位角、臺(tái)站場(chǎng)地響應(yīng)、傳播路徑和震源輻射花瓣效應(yīng)的差異，不同臺(tái)站記錄到的地震波形具有較大差異，使得震源參數(shù)的計(jì)算存在種種不確定性［56］，導(dǎo)致地震視應(yīng)力結(jié)果存在一定的不確定性．但對(duì)于地震活動(dòng)不活躍、不滿足地震活動(dòng)性參數(shù)空間掃描條件的區(qū)域，地震視應(yīng)力在中、強(qiáng)地震危險(xiǎn)區(qū)判定中仍可發(fā)揮重要作用．

致 謝兩位審稿專家提出了具有建設(shè)性的修改建議．視應(yīng)力計(jì)算得到中國(guó)地震局地球物理研究所李艷娥同志的大力幫助．特此致謝！

（

）

［1］ Xu X W，Wen X Z，Yu G H，et al．Coseismic reverse－and oblique－slip surface faulting generated by the 2008Mw7．9 Wenchuan earthquake，China．Geology，2009，37（6）：515－518．

［2］ Zhang P Z，Wen X Z，Shen Z K，et al．Oblique，high－angle，Listric－reverse faulting and associated development of strain：The Wenchuan Earthquake of May 12，2008，Sichuan，China．Annu．Rev．EarthPlanet．Sci．，2010，38（1）：353－382．

［3］ 王衛(wèi)民，趙連鋒，李娟等．四川汶川8．0級(jí)地震震源過(guò)程．地球物理學(xué)報(bào)，2008，51（5）：1403－1410．

Wang W M，Zhao L F，Li J，et al．Rupture process of theMs8．0 Wenchuan earthquake of Sichuan，China．ChineseJ．Geophys．（in Chinese），2008，51（5）：1403－1410．

［4］ 易桂喜，聞學(xué)澤，辛華等．2008年汶川MS8．0地震前龍門(mén)山－岷山構(gòu)造帶的地震活動(dòng)性參數(shù)與地震視應(yīng)力分布．地球物理學(xué)報(bào)，2011，54（6）：1490－1500．

Yi G X，Wen X Z，Xin H，et al．Distributions of seismicity parameters and seismic apparent stresses on the Longmenshan－Minshan tectonic zone before the 2008MS8．0 Wenchuan earthquake．ChineseJ．Geophys．（in Chinese），2011，54（6）：1490－1500．

［5］ 聞學(xué)澤，張培震，杜方等．2008年汶川8．0級(jí)地震發(fā)生的歷史與現(xiàn)今地震活動(dòng)背景．地球物理學(xué)報(bào)，2009，52（2）：444－454．

Wen X Z，Zhang P Z，Du F，et al．The background of historical and modern seismic activities of the occurrence of the 2008Ms8．0Wenchuan，Sichuan，earthquake．ChineseJ．Geophys．（in Chinese），2009，52（2）：444－454．

［6］ Parsons T，Ji C，Kirby E．Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin．Nature，2008，454（7023）：509－510．

［7］ Toda S，Lin J，Meghraoui M，et al．12 May 2008Ms7．9 Wenchuan，China，earthquake calculated to increase failure stress and seismicity rate on three major fault systems．Geophys．Res．Lett．，2008，35：L17305，doi：10．1029／2008GL034903．

［8］ 邵志剛，周龍泉，蔣長(zhǎng)勝等．2008年汶川Ms8．0地震對(duì)周邊斷層地震活動(dòng)的影響．地球物理學(xué)報(bào)，2010，53（8）：1784－1795．

Shao Z G，Zhou L Q，Jiang C S，et al．The impact of WenchuanMs8．0 earthquake on the seismic activity of surrounding faults．ChineseJ．Geophys．（in Chinese），2010，53（8）：1784－1795．

［9］ 單斌，熊熊，鄭勇等．2008年5月12日Mw7．9汶川地震導(dǎo)致的周邊斷層應(yīng)力變化．中國(guó)科學(xué)D 輯：地球科學(xué)，2009，39（5）：537－545．

Shan B，Xiong X，Zheng Y，et al．Stress changes on major faults caused byMw7．9 Wenchuan earthquake，May 12，2008．Sci．ChinaSer．D－EarthSci．，2009，52（5）：593－601．

［10］ 萬(wàn)永革，沈正康，盛書(shū)中等．2008年汶川大地震對(duì)周?chē)鷶鄬拥挠绊懀卣饘W(xué)報(bào)，2009，31（2）：128－139．

Wang Y G，Shen Z K，Sheng S Z，et al．The influence of 2008 Wenchuan earthquake on surrounding faults．Acta SeismologicaSinica（in Chinese），2009，31（2）：128－139．

［11］ Xie C D，Zhu Y Q，Lei X L，et al．Pattern of stress change and its effect on seismicity rate caused byMs8．0 Wenchuan earthquake．Sci．ChinaEarthSci．，2010，53（9）：1260－1270．

［12］ Nalbant S S，McCloskey J．Stress evolution before and after the 2008 Wenchuan，China earthquake．EarthPlanet．Sci．Lett．，2011，307（1－2）：222－232．

［13］ Zhu H，Wen X Z．Static stress triggering effects related withMs8．0 Wenchuan earthquake．J．EarthSci．，2010，21（1）：32－41．

［14］ Aki K．Asperities，barriers，characteristic earthquakes and strong motion prediction．J．GeophysRes．，1984，89（B7）：5867－5872．

［15］ Wiemer S， Wyss M． Mapping the frequency－magnitude distribution in asperities：An improved technique to calculate recurrence times？J．Geophys．Res．，1997，102（B7）：15115－15128．

［16］ Wyss M，Schorlemmer D，Wiemer S．Mapping asperities by minima of local recurrence time：San Jacinto－Elsinore fault zones．J．Geophys．Res．，2000，105（B4）：7829－7844．

［17］ Scholz C H． The frequency－magnitude relation of microfracturing in rock and its relation to earthquakes．Bull．Seism．Soc．Am．，1968，58（1）：399－415．

［18］ Wyss M． Towards a physical understanding of the earthquake frequency distribution．Geophys．J．Int．，1973，31（4）：341－359．

［19］ Urbancic T I，Trifu C I，Long J M，et al．Space－time correlation ofbvalues with stress release．PureAppl．Geophys．，1992，139（3－4）：449－462．

［20］ Gutenberg R，Richter C F．Frequency of earthquakes in California．Bull．Seism．Soc．Am．，1944，34（4）：185－188．

［21］ Wyss M， Wiemer S． Change in the probability for earthquakes in southern California due to the Landers magnitude 7．3earthquake．Science，2000，290（5495）：1334－1338．

［22］ Zu?iga F R，Wyss M．Most－and least－likely locations of large to great earthquakes along the Pacific coast of Mexico estimated from local recurrence times based onb－values．Bull．Seism．Soc．Am．，2001，91（6）：1717－1728．

［23］ Schorlemmer D， Wiemer S， Wyss M． Variations in earthquake－size distribution across different stress regimes．Nature，2005，437（7058）：539－542．

［24］ Wyss M， Stefansson R． Nucleation points of recent mainshocks in Southern Iceland，mapped byb－values．Bull．Seism．Soc．Am．，2006，96（2）：599－608．

［25］ Ghosh A，Newman A V，Thomas A M，et al．Interface locking along the subduction megathrust fromb－value mapping near Nicoya Peninsula，Costa Rica．Geophys．Res．Lett．，2008，35：L01301，doi：10．1029／2007GL031617．

［26］ 易桂喜，聞學(xué)澤，范軍等．由地震活動(dòng)參數(shù)分析安寧河—?jiǎng)t木河斷裂帶的現(xiàn)今活動(dòng)習(xí)性及地震危險(xiǎn)性．地震學(xué)報(bào)，2004，26（3）：294－303．

Yi G X，Wen X Z，F(xiàn)an J，et al．Assessing current faulting behaviors and seismic risk of the Anninghe－Zemuhe fault zone from seismicity parameters．ActaSeismologicaSinica（in Chinese），2004，17（3）：322－333．

［27］ 易桂喜，聞學(xué)澤，王思維等．由地震活動(dòng)參數(shù)分析龍門(mén)山－岷山斷裂帶的現(xiàn)今活動(dòng)習(xí)性與強(qiáng)震危險(xiǎn)性．中國(guó)地震，2006，22（2）：117－125．

Yi G X，Wen X Z，Wang S W，et al．Study on fault sliding behaviors and strong－earthquake risk of the Longmenshan－Minshan fault zones from current seismicity parameters．EarthquakeResearchinChina（in Chinese），2006，22（2）：117－125．

［28］ Zhao Y Z， Wu Z L． Mapping theb－values along the Longmenshan fault zone before and after the 12 May 2008，Wenchuan，China，MS8．0earthquake．Nat．HazardsEarth Syst．Sci．，2008，8：1375－1385．

［29］ M7專項(xiàng)工作組．中國(guó)大陸大地震中－長(zhǎng)期危險(xiǎn)性研究．北京：地震出版社，2012：336．

Working Group of M7．Study on the Mid－to Long－term Potential of Large Earthquakes on the Chinese Continent（in Chinese）．Beijing：Seismological Press，2012：336．

［30］ 龍鋒，蔣長(zhǎng)勝，馮建剛等．歷史大地震破裂區(qū)地震危險(xiǎn)性的地震活動(dòng)性定量分析——以南北地震帶中北段為例．地震，2012，3（3）：98－108．

Long F，Jiang C S，F(xiàn)eng J G，et al．Quantitative seismicity analysis for the risk of historical large earthquake rupture zones：applied to the mid－north segment of the South－North Seismic Belt．Earthquake（in Chinese），2012，3（3）：98－108．

［31］ 吳忠良，黃靜，林碧蒼．中國(guó)西部地震視應(yīng)力的空間分布．地震學(xué)報(bào)，2002，24（3）：293－301．

Wu Z L，Huang J，Lin B C．Distribution of apparent stress in western China．ActaSeismologicaSinica（in Chinese），2002，24（3）：293－301．

［32］ 陳學(xué)忠，王小平，王林瑛等．地震視應(yīng)力用于震后趨勢(shì)快速判定的可能性．國(guó)際地震動(dòng)態(tài)，2003，（7）：1－4．

Chen X Z，Wang X P，Wang L Y，et al．Possibility of earlier judgement on seismic tendency after strong earthquakes upon their apparent stresses．RecentDevelopmentsinWorld Seismology（in Chinese），2003，（7）：1－4．

［33］ 喬慧珍，程萬(wàn)正，陳學(xué)忠．安寧河—?jiǎng)t木河斷裂帶地震視應(yīng)力研究．地震研究，2006，29（2）：125－130．

Qiao H Z，Chen W Z，Chen X Z．Study on the apparent stress of the earthquakes on Anninghe－Zemuhe fault．J．Seism．Res．（in Chinese），2006，29（2）：125－130．

［34］ 李艷娥，陳學(xué)忠，付虹．2007年云南濘洱Ms6．4地震前滇西南地區(qū)震源動(dòng)力學(xué)參數(shù)時(shí)空變化特征．地震，2012，32（1）：28－39．

Li Y E，Chen X Z，F(xiàn)u H．Temporal and spatial variation of focal dynamic parameters in southwest Yunnan before the 2007Ms6．4Ning′er earthquake．Earthquake（in Chinese），2012，32（1）：28－39．

［35］ 李艷娥，陳學(xué)忠，王恒信．汶川8．0級(jí)地震前四川地區(qū)地震視應(yīng)力時(shí)空變化特征．地震，2012，32（4）：113－122．

Li Y E，Chen X Z，Wang H X．Temporal and spatial variation of apparent stress in Sichuan area before theMs8．0 Wenchuan earthquake．Earthquake（in Chinese），2012，32（4）：113－122．

［36］ 王瓊，陳學(xué)忠，王琳瑛．岫巖地震序列視應(yīng)力變化及其預(yù)測(cè)意義．地震，2005，25（2）：91－97．

Wang Q，Chen X Z，Wang L Y．Change of apparent stress of Xiuyan earthquake sequence and its prediction significance．Earthquake（in Chinese），2005，25（2）：91－97．

［37］ 易桂喜，龍鋒，張致偉．汶川Ms8．0地震余震震源機(jī)制時(shí)空分布特征．地球物理學(xué)報(bào)，2012，55（4）：1213－1227．

Yi G X，Long F，Zhang Z W．Spatial and temporal variation of focal mechanisms for aftershocks of the 2008Ms8．0 Wenchuan earthquake．ChineseJ．Geophys．（in Chinese），2012，55（4）：1213－1227．

［38］ Wyss M，Brune J N．Seismic moment，stress，and source dimensions for earthquakes in the California－Nevada region．J．Geophys．Res．，1968，73（14）：4681－4694．

［39］ Choy G L，Boatwright J L．Global patterns of radiated seismic energy and apparent stress．J．Geophys．Res．，1995，100（B9）：18205－18228．

［40］ Wu Z L．Scaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence．EarthPlanetsSpace，2001，53：943－948．

［41］ Brune J N．Tectonic stress and the spectra of seismic shear waves from earthquakes．J．Geophys．Res．，1970，75（26）：4997－5009．

［42］ Jiménez A，García J M，Romacho M D．Simultaneous inversion of source parameters and attenuation factor using genetic algorithms．Bull．Seism．Soc．Am．，2005，95（4）：1401－1411．

［43］ Snoke J A．Stable determination of（Brune）stress drops．Bull．Seism．Soc．Am．，1987，77（2）：530－538．

［44］ Aki K，Richards P G．Quantitative Seismology，Theory and Methods．San Francisco：Freeman ＆Company，1980．

［45］ Archuleta R J，Cranswick E，Mueller C，et al．Source parameters of the 1980Mammoth Lakes，California，earthquakes sequence．J．Geophys．Res．，1982，87（B6）：4595－4607．

［46］ 吳滿路，張?jiān)罉?，廖椿庭等．汶川地震后沿龍門(mén)山裂斷帶原地應(yīng)力測(cè)量初步結(jié)果．地質(zhì)學(xué)報(bào)，2010，84（9）：1291－1299．

Wu M L，Zhang Y Q，Liao C T，et al．Preliminary results of in situ stress measurements along the Longmenshan Fault Zone after the WenchuanMs8．0earthquake．ActaGeologica Sinica（in Chinese），2010，84（9）：1291－1299．

［47］ King G C P，Stein R S，Lin J．Static stress changes and the triggering of earthquakes．Bull．Seism．Soc．Am．，1994，84（3）：935－953．

［48］ Harris R A．Introduction to special section：stress triggers，stress shadows，and implications for seismic hazard．J．Geophys．Res．，1998，103（B10）：24347－24358．

［49］ 華 衛(wèi)，陳章立，李志雄等．汶川8．0級(jí)地震觸發(fā)與余震活動(dòng)空間分布研究．地震，2009，29（1）：33－39．

Hua W，Chen Z L，Li Z X，et al．Seismic triggering and the aftershock distribution of the WenchuanM8．0earthquake．Earthquake（in Chinese），2009，29（1）：33－39．

［50］ Utkucu M， Nalbant S S， McCloskey J，et al．Slip distribution and stress changes associated with the 1999 November 12，Düzce （Turkey）earthquake （Mw＝7．1）．Geophys．J．Int．，2003，153（1）：229－241．

［51］ Daniel G，Marsan D，Bouchon M．Perturbation of the Izmit earthquake aftershock decaying activity following the 1999Mw7．2 Düzce，Turkey，earthquake．J．Geophys．Res．，2006，111，B05310，doi：10．1029／2005JB003978．

［52］ Pollitz F F，Banerjee P，Bürgmann R，et al．Stress changes along the Sunda trench following the 26 December 2004 Sumatra－Andaman and 28 March 2005 Nias earthquakes．Geophys．Res．Lett．，2006，33（6），doi：10．1029／2005GL024558．

［53］ 郭啟良，王成虎，馬洪生等．汶川Ms8．0級(jí)大震前后的水壓致裂原地應(yīng)力測(cè)量．地球物理學(xué)報(bào)，2009，52（5）：1395－1401．

Guo Q L，Wang C H，Ma H S，et al．In－situ hydro－fracture stress measurement before and after the WenchuanMs8．0 earthquake of China．ChineseJ．Geophys．（in Chinese），2009，52（5）：1395－1401．

［54］ 趙靜，江在森，武艷強(qiáng)等．汶川地震前龍門(mén)山斷裂帶閉鎖程度和滑動(dòng)虧損分布研究．地球物理學(xué)報(bào)，2012，55（9）：2963－2972．

Zao J，Jiang Z S，Wu Y Q，et al．Study on fault locking and fault slip deficit of the Longmenshan fault zone before the Wenchuan earthquake．ChineseJ．Geophys．（in Chinese），2012，55（9）：2963－2972．

［55］ 易桂喜，范軍，聞學(xué)澤．由現(xiàn)今地震活動(dòng)分析鮮水河斷裂帶中－南段活動(dòng)習(xí)性與強(qiáng)震危險(xiǎn)地段．地震，2005，25（1）：58－66．

Yi G X，F(xiàn)an J，Wen X Z．Study on faulting behavior and fault－segments for pot－ential strong earthquake risk along the central－southern segment of Xianshuihe fault zone based on current seismicity．Earthquake（in Chinese），2005，25（1）：58－66．

［56］ 楊志高，張曉東．紫坪鋪水庫(kù)地區(qū)蓄水前后視應(yīng)力標(biāo)度率變化研究．地球物理學(xué)報(bào)，2010，53（12）：2861－2868．

Yang Z G，Zhang X D．The transition of apparent stress scaling law before and after impounding in Zipingpu reservoir region．ChineseJ．Geophys．（in Chinese），2010，53（12）：2861－2868．