張致偉，喬慧珍，吳 朋，宮 悅，趙小艷

0 引言

地震是地殼構(gòu)造應(yīng)力積累超過了巖層的受力極限，應(yīng)變能快速釋放，斷裂帶破裂失穩(wěn)的結(jié)果，大地震的發(fā)生必然受到構(gòu)造應(yīng)力場的制約，而震源機(jī)制是研究構(gòu)造應(yīng)力場的基礎(chǔ)資料，不同階段的震源機(jī)制解在一定程度上反映了大震孕育過程中構(gòu)造應(yīng)力場的變化。國內(nèi)外諸多學(xué)者研究認(rèn)為基于震源機(jī)制解反演應(yīng)力場的時空變化能夠捕捉到強(qiáng)震前的震兆信息，震源機(jī)制一致性參數(shù)是判斷強(qiáng)震危險性的一個有用判據(jù) (Michael et al．，1990;Hauksson，1994;Wiemer et al．，2002;陳颙，1978;刁桂苓等，1994;程萬正等，2006b;刁桂苓等，2011)。

震源機(jī)制一致性參數(shù)雖然物理含義明確，但由于受地震臺網(wǎng)布局和小震波形信噪比的制約，往往難以給出比較精確的小震震源機(jī)制解。鑒于上述問題，Lund和B?evarsson(2002)提出了采用微震體波譜振幅相關(guān)系數(shù)來研究地震序列震源機(jī)制是否相似的問題。國內(nèi)學(xué)者隨后也將該方法運(yùn)用于強(qiáng)震震例研究中 (朱航等，2006;崔子健等，2012)，并取得了較好的成果。研究表明強(qiáng)震前的小震群序列與一般性小震群序列相比，譜振幅相關(guān)系數(shù)接近于1．0，說明震源區(qū)應(yīng)力場增強(qiáng)，震源機(jī)制相似性更為顯著。且譜振幅相關(guān)系數(shù)的測定和分析，只需地震序列周圍若干數(shù)字地震臺站即可，具有一定的可操作性。

四川東南部存在多處采鹽、采氣及廢水回注井，地震活動普遍表現(xiàn)為強(qiáng)度低、頻次高、持續(xù)時間長的特征，且與當(dāng)?shù)丶訅鹤⑺安蓺?(鹵)等工業(yè)活動具有很好的相關(guān)性 (阮祥等，2008;張致偉等，2012)。本文選取四川自貢和長寧兩個注水地區(qū)作為研究區(qū)域，基于測震臺網(wǎng)記錄的數(shù)字地震波形資料和觀測報告，分析注水誘發(fā)地震的譜振幅相關(guān)系數(shù)特征及其與當(dāng)?shù)丶訅鹤⑺當(dāng)?shù)據(jù)的關(guān)系，并結(jié)合表征絕對應(yīng)力水平的視應(yīng)力綜合分析研究區(qū)域的應(yīng)力高低。

1 資料和方法

1． 1 研究區(qū)地震活動概況

本文以自貢和長寧兩個注水地區(qū)為例 (圖1a中圓形區(qū)域)。2000年以來研究區(qū)域ML≥2．0地震時序圖 (圖1b)顯示，自貢地區(qū)在2000～2008年地震活動較弱，未發(fā)生ML≥3．0地震，2009年1月上旬該區(qū)小震活動明顯增強(qiáng)，并于2009年2月16日、5月22日先后發(fā)生了ML4．4和 ML4．2地震，且在2011、2012年又發(fā)生了兩次比較顯著的4級地震震群活動 (圖1中圓形區(qū)域b)，張致偉等 (2012)研究表明該區(qū)地震活動主要受當(dāng)?shù)貜U水回注井的注水壓力調(diào)控。長寧地區(qū)在2000～2006年3月地震活動同樣呈現(xiàn)相對較弱狀態(tài)，2006年4月以來該區(qū)域地震活動頻次和強(qiáng)度明顯升高，并于2008年2月1日發(fā)生了ML4．8地震，隨后該區(qū)域地震頻次略有下降(圖1c)，但ML≥4．0地震依然持續(xù)發(fā)生，并于2013年4月25日發(fā)生了ML5．2地震 (圖1a中圓形區(qū)域c)，阮祥等 (2008)指出對該區(qū)域的地震活動阮祥等 (2008)研究表明該區(qū)域的地震活動主要受長寧雙河鹽礦注 (出)水量的影響。

1． 2 譜振幅相關(guān)系數(shù)測定方法

地震波能量在傳播過程中不但隨著傳播距離出現(xiàn)幾何擴(kuò)散，還受到傳播路徑上介質(zhì)的吸收和散射等影響，在到達(dá)臺站接收儀器前地表下方淺層的介質(zhì)也會對地震波產(chǎn)生影響，因此觀測得到的地震數(shù)據(jù)是震源激發(fā)的信息經(jīng)過上述各種過程的產(chǎn)物。在頻率域，臺站j記錄到的地震i的觀測位移譜Uij(f)可表示為

式中，S'i(f)為地震i的震源譜;φij為地震i的震源輻射圖型因子;Pij(f)為傳播路徑效應(yīng)，包括地震波的幾何擴(kuò)散和非彈性衰減;L'j(f)為第j個臺站的局部場地效應(yīng);N'j(f)為第j個臺站附近的地面運(yùn)動噪聲;I'j(f)為第j個臺站儀器響應(yīng);Surj為臺站附近地表自由表面效應(yīng)。

根據(jù)Brune模型 (1970)，S'i(f)可表示為

式中，Ω'0i為地震i的震源譜低頻漸近線值，即零頻極限值;fc為低頻漸近線與高頻漸近線交點(diǎn)處的頻率，稱為拐角頻率。

定義譜振幅Ω0ij為臺站j記錄到的地震i的包含輻射圖型因子的震源譜零頻極限值，可表示為

式 (1)中的Pij(f)可表示為

式中，Gij為幾何擴(kuò)散因子;Rij為震源距;Q(f)為介質(zhì)的品質(zhì)因子;v為地震波 P波或 S波傳播速度。

從地震記錄中扣除噪聲和儀器響應(yīng)，并令Lj(f)=L'j(f)Surj)(Lj稱為廣義的臺站場地效應(yīng))，地震i在臺站j的地面運(yùn)動位移譜Uij(f)為

可見，若臺站j記錄的兩次地震x、y震源位置足夠接近，即其間距比震源距Rij小得多，且震源機(jī)制 (輻射圖型因子)相同，則臺站j記錄的這兩次地震的地面運(yùn)動位移譜Uxj(f)與Uyj(f)的相對大小只與這兩次地震的Ω'0相對大小相關(guān);若震源機(jī)制不同，則還應(yīng)與這兩次地震震源機(jī)制的差異有關(guān)。反過來說，當(dāng)由地震記錄來反演譜振幅時，若這兩次地震震源機(jī)制相同，則反演得到的Ω0xj與Ω0yj的相對大小只與Ω'0x和Ω'0y的相對大小有關(guān);如果震源機(jī)制不同，則與震源機(jī)制有關(guān)。因此可由臺站的地震記錄反演地震x、y的振幅譜Ω0xj與Ω0yj，通過計(jì)算其相關(guān)系數(shù)rxy來描述震源機(jī)制是否相似的問題，若震源機(jī)制相似，譜振幅相關(guān)系數(shù)rxy應(yīng)較大，接近1;若震源機(jī)制不相似，則相關(guān)系數(shù)較小。地震x和y的譜振幅相關(guān)系數(shù) rxy可表示為(Lund，B?evarsson，2002)

式中，xjl、yjl分別表示臺站j記錄的地震x、y的l分量波形資料所反演得到的譜振幅Ω0xj與Ω0yj的對數(shù)值;x、y分別為xjl與tjl的平均值;n為使用的臺站數(shù)目;l取值為1到5，依次為垂向和徑向P波，垂向、徑向和切向S波。

本研究考慮到持續(xù)地震震源機(jī)制是否相似是針對多數(shù)地震而言的，因此采用崔子健等 (2012)提出的計(jì)算方法，將地震按發(fā)生時間順序排列，地震m與前面m-1個地震為一組，計(jì)算組內(nèi)每兩個地震的rxy，得出N=m(m-1)/2個相關(guān)系數(shù)rxy，對rxy求算數(shù)平均值，其結(jié)果代表了地震m時刻的組內(nèi)譜振幅的相關(guān)程度，以步長為1進(jìn)行滑動，計(jì)算每個組內(nèi)的譜振幅相關(guān)系數(shù)的算數(shù)平均值，最后可以得到隨時間變化的譜振幅相關(guān)系數(shù)。

2 譜振幅相關(guān)系數(shù)結(jié)果分析

2． 1 場地響應(yīng)和區(qū)域Q值反演

地震儀器記錄到的地震波形資料包含震源信息、傳播路徑及臺站場地響應(yīng)等，要想獲得趨于更真實(shí)的震源參數(shù)，必須對觀測信號做地震波的幾何擴(kuò)散、非彈性衰減以及臺站場地響應(yīng)校正 (劉杰等，2003)?；?009-1-1～2012-12-31四川東南地區(qū)20次ML≥3．5地震波形資料及觀測報告，結(jié)合Atkinson和Mereu(1992)，Moya等 (2000)提出的方法，反演獲得了川東南地區(qū)P、S波Q值和7個測震臺站的場地響應(yīng)。

Q值是描述地球介質(zhì)特性的主要參數(shù)，與地質(zhì)構(gòu)造、地震活動性及地域熱流有密切的關(guān)系 (周連慶等，2008)。Atkison方法是考慮頻率依賴的Q模型，通常用頻率的冪函數(shù)，即Q(f)=Q0fη的形式來擬合Q值與頻率f的關(guān)系，參數(shù)η反映了Q值對頻率的依賴程度。圖2給出了川東南地區(qū)P、S波Q值隨頻率的變化及其擬合曲線，P波Q=53．8f0．879，S 波 Q=101．9f1．025，與阮祥等 (2008)基于早期資料計(jì)算獲得長寧地區(qū)的尾波Q值比較一致，而與川西北地區(qū)的Q值結(jié)果 (喬慧珍等，2006)差別較大，表明Q值具有地區(qū)差異特征。

圖3給出了川東南地區(qū)7個測震臺站的場地響應(yīng)結(jié)果，圖中藍(lán)線是由每個地震記錄得到的臺站場地響應(yīng)，紅線則是該臺站場地響應(yīng)的平均值。臺站場地響應(yīng)總體表現(xiàn)為高頻衰減特征，JLI、MBI、LBO和LZH臺在1～4 Hz的放大倍數(shù)在1～2之間，而其余臺站的放大效應(yīng)都大于2。臺站場地響應(yīng)形態(tài)之所以存在差異，主要受臺基、周邊地形等因素的影響。臺站場地響應(yīng)的計(jì)算不僅對研究區(qū)臺站場地選點(diǎn)有一定的參考意義，而且也為下一步計(jì)算譜振幅相關(guān)系數(shù)奠定基礎(chǔ)。

2． 2 注水誘發(fā)地震的譜振幅相關(guān)系數(shù)

在獲得川東南地區(qū)P、S波的Q值及臺站場地響應(yīng)的基礎(chǔ)上，以自貢和長寧兩個注水地區(qū)為例，采用上述方法研究注水誘發(fā)地震的譜振幅相關(guān)系數(shù)特征及其與當(dāng)?shù)丶訅鹤⑺當(dāng)?shù)據(jù)的關(guān)系。

選取自貢地區(qū)2009-1-1～2013-4-30共32次 ML≥3．0地震，基于馬邊 (MBI)、漢王山(HWS)、花馬寺 (HMS)、五馬坪 (WMP)和瀘州臺 (LZH)記錄的地震波資料和觀測報告，計(jì)算獲得該區(qū)域地震譜振幅相關(guān)系數(shù)在0．84～0．92之間波動，其平均值為0．89。結(jié)合研究區(qū)地震活動及注水井的加壓數(shù)據(jù)，將研究時段分為兩個階段進(jìn)行分析 (圖4a)，第一階段為2009-1～2010-7，該時段注水壓力持續(xù)升高，地震活動顯著，譜振幅相關(guān)系數(shù)由0．84逐漸上升至0．92，呈現(xiàn)隨加壓數(shù)據(jù)的升高而增大的特征;第二階段為2010-8～2013-4，該時段注水井停止加壓，地震活動明顯減弱，譜振幅相關(guān)系數(shù)也有逐漸降低的趨勢。

選取長寧地區(qū)2006-1-1～2013-4-30共88次 ML≥3．0地震，基于馬邊 (MBI)、雷波(LBO)、花馬寺 (HMS)、五馬坪 (WMP)和瀘州臺 (LZH)記錄的地震波資料和觀測報告計(jì)算獲得了譜振幅相關(guān)系數(shù)，該區(qū)域地震相關(guān)系數(shù)總體上在0．87～0．97之間波動，其平均值為 0．90。結(jié)合研究區(qū)地震活動及鹽礦注 (出)水量數(shù)據(jù)，同樣將研究時段分為兩個階段 (圖4b)，第一階段為2006-4～2007-11，相比前期注 (出)水量，該時段內(nèi)注水量大幅增加，但出水量的增加并不明顯，造成注水量和出水量之差大幅上升，小震活動明顯增強(qiáng)，相應(yīng)的譜振幅相關(guān)系數(shù)也在高值(0．90～0．97)波動，而且在注 (出)水量差值變化最大時相關(guān)系數(shù)出現(xiàn)了明顯的階躍;第二階段為2007-12～2013-4，該時段注水量和出水量相當(dāng)，且有逐漸下降的趨勢，但地震活動依然持續(xù)，譜振幅相關(guān)系數(shù)隨著注 (出)水量的下降也有所降低，但仍然大于0．87，且比較穩(wěn)定。

綜上所述，自貢和長寧地區(qū)注水誘發(fā)地震的譜振幅相關(guān)系數(shù)均在高值波動，其均值分別為0．89和0．90，且變化過程與加壓注水?dāng)?shù)據(jù)具有很好的相關(guān)性。Seeber和Armbruster(2000)認(rèn)為地震可作為應(yīng)力變化的標(biāo)志，雖然小地震對總的變形貢獻(xiàn)很小，但它們在空間、時間和運(yùn)動上的分布對應(yīng)力變化非常敏感。影響地殼力學(xué)狀態(tài)的各種現(xiàn)象，包括諸如斷層蠕動或巖漿注入的自然現(xiàn)象，如水庫蓄水或油田注水的大型工程活動。自貢和長寧地區(qū)的介質(zhì)因加壓注水、采鹽等工業(yè)活動已發(fā)生了變化，注水對局部區(qū)域應(yīng)力場產(chǎn)生擾動影響，使地下淺層裂隙呈優(yōu)勢取向排列，引發(fā)的中小地震震源機(jī)制表現(xiàn)出較好的相似性，體現(xiàn)在較高的譜振幅相關(guān)系數(shù)上。

3 地震視應(yīng)力結(jié)果分析

地震視應(yīng)力是在有關(guān)震源介質(zhì)的均勻彈性、地震破裂的脆性性質(zhì)假定下 (Kanamori，Anderson，1975;Kanamori，1977;Wu，2001)，可由地震波信息獲得的、與區(qū)域平均應(yīng)力呈正比的物理量。對一個地區(qū)中引起地震滑動的應(yīng)力水平進(jìn)行區(qū)域平均，則可作為當(dāng)?shù)亟^對應(yīng)力水平的一個間接估計(jì) (吳忠良等，2002)。

地震視應(yīng)力σapp被定義為單位地震矩輻射出的地震波能量與震源區(qū)介質(zhì)的剪切模型的積 (Wyss，Brune，1968)，即

式中，ES為地震波輻射能量;M0為地震矩;μ為震源區(qū)介質(zhì)剪切模量，通常取 μ=3．5×1011dyn/cm2。

基于Brune(1970)圓盤模型，地震震源譜可由震源譜的低頻水平Ω'oi和拐角頻率fc表示 (式(2))。而地震波輻射能量ES和地震矩M0分別表示為

式中，β為S波波速，取β=3．5 km/s;ρ為介質(zhì)密度，取ρ=2．71 g/cm3;Rθφ為輻射圖型因子，取S波的均方根輻射圖形因子■2/5≈0．63。

計(jì)算過程中，一般選取區(qū)域地震臺網(wǎng)震中距≤200 km的臺站記錄的具有較高信噪比的寬頻帶數(shù)字化地震波形資料，利用S波頻帶范圍在1．0～20 Hz的波形計(jì)算震源譜，在獲得地震波輻射能量和地震矩的基礎(chǔ)上，再由 (7)式逐臺求取地震視應(yīng)力，最終取各臺計(jì)算結(jié)果的平均值作為相應(yīng)地震的視應(yīng)力。

自貢和長寧地區(qū)注水誘發(fā)地震的譜振幅相關(guān)系數(shù)較高，可能反映了中小地震震源機(jī)制具有很好的相似性，為了配合譜振幅相關(guān)系數(shù)分析研究區(qū)域的應(yīng)力水平，基于上述方法計(jì)算獲得了2000～2013年川東南地區(qū)ML≥3．0地震的視應(yīng)力結(jié)果。圖5分別給出了獲得地震視應(yīng)力的ML≥3．0地震及其視應(yīng)力空間分布圖，由圖5a可見，3級地震空間分布比較分散，ML≥4．0地震則主要集中分布在四川自貢、長寧、重慶榮昌和云南東北地區(qū)。視應(yīng)力相對高值區(qū)均在發(fā)生4級地震的地區(qū)，唯獨(dú)四川自貢和長寧4級地震集中區(qū) (圓圈)的視應(yīng)力較低 (圖5b)，表明川東南局部地區(qū)地震視應(yīng)力存在差異，自貢和長寧地區(qū)的應(yīng)力水平相對偏低。

為了進(jìn)一步研究自貢和長寧注水地區(qū)地震視應(yīng)力的特征，圖6給出了自貢、長寧及川東南地區(qū)ML≥3．5地震的視應(yīng)力對比圖。由圖可見，川東南地區(qū)地震視應(yīng)力總體呈現(xiàn)隨震級增大而增大的特征，但與程萬正等 (2006a)給出的四川其它地區(qū)視應(yīng)力相比，川東南地區(qū)同等震級的地震視應(yīng)力相對偏低。自貢地區(qū)不同震級檔的地震視應(yīng)力均在川東南地震視應(yīng)力擬合線下方，長寧地區(qū)不同震級檔的地震視應(yīng)力分布在擬合線附近，與川東南其它區(qū)域地震視應(yīng)力相比，自貢和長寧地區(qū)同等震級的地震視應(yīng)力相對偏低。地殼介質(zhì)應(yīng)力水平有隨深度增加而變大的特征 (丁健民，高莉青，1981;李方全，祁英男，1988)，自貢和長寧地區(qū)小震精定位均顯示其震源深度較淺 (阮祥等，2008;張致偉等，2013)，可作為研究區(qū)域應(yīng)力水平較低的佐證。

結(jié)合自貢和長寧地區(qū)地震活動及加壓、注水?dāng)?shù)據(jù)，分析研究區(qū)域地震視應(yīng)力在不同階段的變化特征 (圖7)。自貢地區(qū)除發(fā)生在第Ⅱ時段內(nèi)的14、15和17日ML≥4．5地震視應(yīng)力比較顯著外，加壓時段 (Ⅰ)與停止加壓時段 (Ⅱ)內(nèi)同等震級檔的地震視應(yīng)力相當(dāng) (圖7a);長寧地區(qū)與上述現(xiàn)象相似，除第Ⅱ時段內(nèi)的13和17日ML≥4．5地震外，整個時段內(nèi)同等震級的地震視應(yīng)力相差不大 (圖7b)。可見自貢和長寧地區(qū)地震視應(yīng)力同樣依賴于震級大小，加壓注水等工業(yè)活動對該地區(qū)的地震視應(yīng)力影響不大。

4 結(jié)論與討論

本文以四川自貢和長寧兩個注水地區(qū)為例，基于區(qū)域臺網(wǎng)記錄的數(shù)字地震波形資料和觀測報告，分析注水誘發(fā)地震的譜振幅相關(guān)系數(shù)特征及其與當(dāng)?shù)丶訅鹤⑺當(dāng)?shù)據(jù)的關(guān)系，并結(jié)合表征絕對應(yīng)力水平的視應(yīng)力綜合分析研究區(qū)域的應(yīng)力高低。

自貢和長寧地區(qū)注水誘發(fā)地震的譜振幅相關(guān)系數(shù)均在高值波動，且變化過程與加壓注水?dāng)?shù)據(jù)具有很好的相關(guān)性。自貢地區(qū)地震譜振幅相關(guān)系數(shù)在0．84～0．92之間波動，其平均值為0．89，在注水壓力持續(xù)升高的背景下，譜振幅相關(guān)系數(shù)呈現(xiàn)隨加壓數(shù)據(jù)的升高而增大的特征;長寧地區(qū)地震譜振幅相關(guān)系數(shù)在0．87～0．97之間波動，其平均值為0．90，在注水量和出水量之差大幅上升時段，相應(yīng)的譜振幅相關(guān)系數(shù)也維持在高值。而在研究區(qū)域注 (出)水量相當(dāng)或停止加壓時，譜振幅相關(guān)系數(shù)也有相應(yīng)降低，但依然在高值波動。分析認(rèn)為加壓注水對局部區(qū)域應(yīng)力場產(chǎn)生擾動影響，使地下淺層裂隙呈優(yōu)勢取向排列，引發(fā)的中小地震震源機(jī)制表現(xiàn)出較好的相似性，從而表現(xiàn)為較高的譜振幅相關(guān)系數(shù)。

川東南地區(qū)地震視應(yīng)力存在局部區(qū)域差異，自貢和長寧地區(qū)同等震級的地震視應(yīng)力與其它區(qū)域結(jié)果相比相對偏低，可能反映了上述區(qū)域的應(yīng)力水平相對偏低，研究區(qū)域的地震視應(yīng)力具有隨震級增大而增大的特征，加壓注水等工業(yè)活動對所在區(qū)域地震視應(yīng)力影響不大。

自貢和長寧兩個注水地區(qū)地震的譜振幅相關(guān)系數(shù)整體偏高，反映了中小地震震源機(jī)制具有很好的相似性，而視應(yīng)力結(jié)果顯示研究區(qū)域的應(yīng)力水平并不高，且歷史上該區(qū)域的地震活動水平也不高。崔子健等 (2012)研究發(fā)現(xiàn)強(qiáng)震前的小震群同樣具有較高的譜振幅相關(guān)系數(shù)，其反映了大震孕育過程中震源區(qū)介質(zhì)應(yīng)力水平增強(qiáng)。注水地區(qū)及其它強(qiáng)震孕育區(qū)雖然譜振幅相關(guān)系數(shù)具有相同的特征，但其孕震機(jī)理不同，局部區(qū)域應(yīng)力水平也就存在差異。

鑒于小震震源機(jī)制不易求解，理論上譜振幅相關(guān)系數(shù)能夠反映震源機(jī)制的相似性。但小震活動具有隨機(jī)性，其受局部應(yīng)力場及局部構(gòu)造等因素的影響，應(yīng)結(jié)合表征應(yīng)力水平的參數(shù)進(jìn)行綜合判定。

本項(xiàng)工作得到了中國地震局“地震分析預(yù)測研究青年工作組”的大力支持;中國地震臺網(wǎng)中心周龍泉博士和四川省地震局程萬正研究員給予了重要指導(dǎo)和有益討論;中國地震局預(yù)測研究所崔子健和云南省地震局鄔成棟在程序上給予了幫助和指導(dǎo)，在此一并致謝!

陳颙．1978．用震源機(jī)制一致性作為描述地震活動性的新參數(shù)［J］．地球物理學(xué)報，21(2):146-159．

程萬正，陳學(xué)忠，喬慧珍．2006a．四川地震輻射能量和視應(yīng)力的研究［J］．地球物理學(xué)進(jìn)展，21(3):692-699．

程萬正，阮祥，張永久．2006b．川滇次級地塊震源機(jī)制類型與一致性參數(shù)［J］．地震學(xué)報，28(6):561-573．

崔子健，李志雄，陳章立，等．2012．判別小震群序列類型的新方法研究——譜振幅相關(guān)分析法［J］．地球物理學(xué)報，55(5):1718-1724．

刁桂苓，徐錫偉，陳于高，等．2011．汶川 MW7．9和集集MW7．6地震前應(yīng)力場轉(zhuǎn)換現(xiàn)象及其可能的前兆意義［J］．地球物理學(xué)報，54(1):128-136．

刁桂苓，于利民，李欽祖．1994．強(qiáng)震前后震源區(qū)應(yīng)力場變化一例［J］．地震學(xué)報，16(1):64-69．

丁健民，高莉青．1981．地殼水平應(yīng)力與垂直應(yīng)力隨深度的變化［J］．地震，2:46-48．

李方全，祁英男．1988．地殼應(yīng)力隨深度的變化規(guī)律［J］．巖石力學(xué)與工程學(xué)報，7(4):301-309．

劉杰，鄭斯華，黃玉龍，等．2003．利用遺傳算法反演非彈性衰減系數(shù)、震源參數(shù)和場地響應(yīng)［J］．地震學(xué)報，25(2):211-218．

喬慧珍，張永久，程萬正，等．2006．川西北地區(qū)介質(zhì)衰減特征研究［J］．地震地磁觀測與研究，27(4):1-7．

阮祥，程萬正，張永久，等．2008．四川長寧鹽礦井注水誘發(fā)地震研究［J］．中國地震，24(3):226-234．

吳忠良，黃靜，林碧蒼．2002．中國西部地震視應(yīng)力的分布［J］．地震學(xué)報，24(3):293-301．

張致偉，程萬正，梁明劍，等．2012．四川自貢—隆昌地區(qū)注水誘發(fā)地震研究［J］．地球物理學(xué)報，55(5):1635-1645．

張致偉，程萬正，吳朋，等．2013．自貢—隆昌地區(qū)地震重新定位及P波速度結(jié)構(gòu)研究［J］．中國地震，29(1):37-47．

周連慶，趙翠萍，修濟(jì)剛，等．2008．利用天然地震研究地殼Q值的方法和進(jìn)展［J］．國際地震動態(tài)，(2):1-11．

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

Atkinson G．M．，Mereu R．F．．1992．The Shape of Ground Motion Attenuation Curves in Southeastern Canada［J］．Bull．Seism．Soc．Amer．，82(5):2014-2031．

Brune J．N．．1970．Tectonic Stress and the Spectra of Seismic Shear Waves from Earthquakes［J］．J．Geophys．Res．，75(26):4997-5009．

Hauksson E．．1994．State of Stress from Focal Mechanisms before and after the 1992 Landers Earthquake Sequence［J］．Bull．Seism．Soc．Am，84(3):917-934．

Kanamori H．，Anderson D．L．．1975．Theoretical Basis of Some Empirical Relations in Seismology［J］．Bull．Seism．Soc．Am．，65(5):1073-1095．

Kanamori H．．1977．The Energy Release in Great Earthquakes［J］．J．Geophys．Res．，82(20):2981-2987．

Lund B．，B?evarsson R．．2002．Correlation of Microearthquake Bodywave Spectral Amplitudes［J］．Bull．Seism．Soc．Am．，92(6):2419-2433．

Michael A．J．，Ellsworth W．L．，Oppenheimer D．H．．1990．Coseismic Stress Changes Induced by the 1989 Loma Prieta，California Earthquake［J］．Geophys．Res．Lett．，17(9):1441-1444．

Moya C．A．，Aguirre J．，Irikura K．．2000．Inversion of Source Parameters and Site Effects from Strong Ground Motion Records using Genetic Algorithms［J］．Bull．Seism．Soc．Amer．，90(4):977-992．

Seeber L．，Armbruster J．G．．2000．Earthquakes as Beacons of Stress Change［J］．Nature，407(6800):69-72．

Wiemer S．，Gerstenberger M．，Hauksson E．．2002．Properties of the Aftershock Sequence of the 1999 MW7．1 Hector Mine earthquake:Implications for Aftershock Hazard［J］．Bull．Seism．S．Am．，92(4):1227-1240．

Wu Z．L．．2001．Scaling of Apparent Stress from Broadband Radiated Energy Catalogue and Seismic Moment Catalogue and its Focal Mechanism Dependence［J］．Earth．Planets．Space．，53:943-948．

Wyss M．，Brune J．N．．1968．Seismic Moment，Stress，and Source Dimensions forEarthquakes in theCalifornia-NevadaRegion［J］．J．Gephys．Res．，73(14):4681-4694．