周少輝 蔣海昆 曲均浩 李健

1）山東省地震局，濟南市文化東路20號 250014

2）中國地震臺網(wǎng)中心，北京 100045

3）海南省地震局，海口 570203

0 引言

地震是指震源區(qū)介質(zhì)在應(yīng)力作用下產(chǎn)生的破裂錯動或震源區(qū)原有斷層的破裂擴展。換言之，地震的發(fā)生與巖石發(fā)生錯動處的應(yīng)力狀態(tài)有關(guān)。從余震預(yù)測的角度，就整個地震序列而言，地震發(fā)生后震源區(qū)應(yīng)力狀態(tài)的判定，對震后趨勢估計及強余震預(yù)測具有重要意義（鐘羽云等，2004）。由于地震發(fā)生在地下深處，目前技術(shù)條件下無法直接測量震源區(qū)應(yīng)力場的強弱，因此主要根據(jù)震源機制、應(yīng)力降、視應(yīng)力等來綜合研究局部應(yīng)力場強弱、方向等作用方式（陳學(xué)忠，2005）。

近年來，隨著數(shù)字地震觀測的逐漸完善，基于數(shù)字地震資料計算震源及介質(zhì)參數(shù)，進而開展余震預(yù)測的方法逐漸受到重視，將應(yīng)力降用于震源區(qū)地殼應(yīng)力動態(tài)變化研究和地震序列趨勢判斷也取得了很多重要結(jié)果，這方面的研究亦是今后余震物理預(yù)測的一個重要發(fā)展方向（蔣海昆等，2015）。不少研究指出，震前、震后應(yīng)力降的變化與地震活動性有關(guān)（Chinnery，1964；Hanks，1976；Kanamori et al，1975；Richardson et al，1977）。應(yīng)力降表征地震發(fā)生瞬間錯動時位錯面上的應(yīng)力變化，反映地震過程中構(gòu)造應(yīng)力場大小和釋放水平，是直接與地震的發(fā)生、震源介質(zhì)和構(gòu)造應(yīng)力密切相關(guān)的參數(shù)。中小地震應(yīng)力降隨時間的變化可能反映了應(yīng)力狀態(tài)隨時間的變化，對震后趨勢判定具有重要意義（華衛(wèi)，2007）。大震后，若余震區(qū)中小地震應(yīng)力降明顯偏高或有逐步增大的趨勢，可能說明震源區(qū)仍處于較高應(yīng)力狀態(tài)；反之，則可能表示余震區(qū)應(yīng)力狀態(tài)基本恢復(fù)正常（蔣海昆等，2015）。通過計算強震序列中余震的震源參數(shù)，分析其應(yīng)力降隨時間、空間的變化，可研究地震序列的震源物理過程，進而研究震源區(qū)地殼介質(zhì)的物理力學(xué)性質(zhì)、應(yīng)力狀態(tài)動態(tài)變化過程等，而了解中小地震孕育、發(fā)生的應(yīng)力背景，對序列后續(xù)地震活動判定及分析局部地區(qū)的地震危險性有著重要意義（華衛(wèi)等，2009）。因此，從地震預(yù)測的角度深入研究應(yīng)力降的主要特征，進而開展利用應(yīng)力降對序列后續(xù)地震活動及強余震的判定研究，將有利于地震短臨預(yù)測。在查閱大量文獻的基礎(chǔ)上，本文將簡述國內(nèi)外應(yīng)力降研究領(lǐng)域的一些進展，著重整理應(yīng)力降的主要特征，并對將應(yīng)力降應(yīng)用于震源區(qū)地殼應(yīng)力動態(tài)變化研究和地震序列趨勢判斷等方面所取得的最新研究進展進行了簡要綜述。

1 應(yīng)力降基本概念

1.1 應(yīng)力降定義

應(yīng)力降Δσ表征地震發(fā)生瞬間錯動時位錯面上的應(yīng)力變化，定義為地震前、后初始應(yīng)力σ0與終止應(yīng)力σ1之差

Δσ稱為靜態(tài)應(yīng)力降，與靜態(tài)應(yīng)力降相對應(yīng)的還有有效應(yīng)力降（也稱動態(tài)應(yīng)力降），有效應(yīng)力降定義為初始應(yīng)力σ0與動摩擦應(yīng)力σf之差（圖1）

斷層錯動期間斷層面上的平均應(yīng)力為

圖1 應(yīng)力降示意圖

1.2 應(yīng)力降計算方法

主要通過2種方法來計算應(yīng)力降，一種是在比較地震震級與其輻射能或斷層破裂面積估計值的基礎(chǔ)上，運用標(biāo)度關(guān)系估算應(yīng)力降（kanamori et al，1975；Andrews，1986）；另一種是在假定震源模型的基礎(chǔ)上，通過震源譜分析求取模型中的相關(guān)參數(shù)，進而估計應(yīng)力降（Brune，1970）。

在利用標(biāo)度關(guān)系來估計應(yīng)力降的方法中，若假定地震時尺度為L的斷層錯動的位移量為D，則依據(jù)Hooke定律可估計整個斷層上的應(yīng)力降（Mohammadioun，2002）

式中，μ為摩擦系數(shù)；平均滑移量D可由地震矩 M0進行表達，即 D≈cM0／（μL2），其中，斷層形狀參數(shù)c＝L／W，L、W分別為斷層的長與寬。由此，Δσ可表達為

由式（4）、（5）可見，地震應(yīng)力降Δσ與地震滑移量D、斷層尺度L及所釋放的地震矩M0等有關(guān)。

通過震源譜分析的應(yīng)力降計算方法中，對于L×W矩形走滑型斷層模型，應(yīng)力降 Δσ為（蔣海昆等，2015）

對于L×W矩形傾滑型斷層模型，應(yīng)力降Δσ為（蔣海昆等，2015）

對于中小地震，一般使用圓盤斷層模型（Brune，1970、1971），即將地震斷層面與一個半徑為R的圓盤等效，假設(shè)剪應(yīng)力同時作用于整個斷層面，則應(yīng)力降 Δσ為（Brune，1970、1971）

2 應(yīng)力降主要特征

2.1 應(yīng)力降數(shù)值分布特征

地震應(yīng)力降數(shù)值大多為零點幾兆帕至數(shù)十兆帕。Kanamori等（1975）對1923～1968年全球6級以上地震的研究表明，應(yīng)力降為0.1～10.0MPa，平均應(yīng)力降為 6MPa；Purcaru等（1978）對1857～1976年全球 7級以上大地震的研究表明，應(yīng)力降為 2～13MPa；臧紹先（1984）對這2組地震事件進一步分析指出，對于大多數(shù)淺源地震，應(yīng)力降均為幾兆帕，而2～6MPa為主要分布范圍。Allmann等（2009）對1900～2007年全球約2000個5級以上地震的研究表明，應(yīng)力降為0.3～50.0MPa，中位數(shù)為4MPa；Oth等（2010）研究顯示，震級為 2.7～8.0的日本淺源地震，應(yīng)力降為0.1～100.0MPa，平均應(yīng)力降為1.1MPa。中國大陸地區(qū)中小地震應(yīng)力降為0.1～20.0MPa，絕大多數(shù)小于10MPa（趙翠萍等，2011）。

2.2 應(yīng)力降與震級間的關(guān)系

地震釋放的應(yīng)力降與地震震級之間的關(guān)系是震源參數(shù)研究中被廣泛討論的重要問題之一，至今仍未有完全統(tǒng)一的認識。部分研究者認為，對于矩震級MW1～8的地震，其地震矩與破裂尺度間呈簡單的冪律關(guān)系，其應(yīng)力降則近似于常數(shù)，不隨地震震級的改變而改變（Hanks，1976；McGarr，1984；Abercrombie，1995；Ide et al，2001；Allmann et al，2009），此即常數(shù)應(yīng)力降模型的觀測基礎(chǔ)。常數(shù)應(yīng)力降模型認為，地震震級與凹凸體破裂尺度有關(guān)，應(yīng)力降不隨地震矩的變化而明顯改變（Allmann et al，2007、2009；Baltay et al，2011）。但也有許多觀測表明，小地震應(yīng)力降隨震級的增大而增加（Izutani et al，2001；吳忠良等，1999；Hardebeck et al，2009；趙翠萍等，2011；華衛(wèi)等，2012），并不完全與常應(yīng)力降模型的預(yù)期相一致。如Mayeda等（1996）的結(jié)果表明，小地震應(yīng)力降進一步也有研究者指出，應(yīng)力降隨震級的變化存在多重標(biāo)度特征，陳運泰等（2000）曾指出，“大”地震與“小”地震的情況有別：對于不同地震矩的大地震，應(yīng)力降接近于常數(shù)，地震的大小是通過破裂面積和錯距來區(qū)分的；對于 “小”地震，應(yīng)力降大小隨震級（或地震矩）的增大而增大。Atkinson等（2004）對加拿大東南和美國東北地區(qū)的研究表明，MW≥4.0地震的應(yīng)力降與震級間正相關(guān)，MW≤4.0地震則顯示常應(yīng)力降特征。王培玲等（2013）的研究也顯示，應(yīng)力降隨震級的變化可能存在2種標(biāo)度特征：①對3.0級以上地震，應(yīng)力降有很大的變化范圍，應(yīng)力降有隨震級增大而增大的趨向；②對于震級小于3.0的地震，應(yīng)力降大多分布在10MPa附近，隨震級的變化不明顯。值得注意的是，以上討論均是從多樣本平均的角度來進行的，事實上大地震釋放的應(yīng)力降往往也存在較大起伏（臧紹先，1984）。另外，Nuttli（1983）研究認為，板間地震滿足常數(shù)應(yīng)力降模型，即應(yīng)力降不隨地震震級而變化；而板內(nèi)地震則滿足增加應(yīng)力降模型，即應(yīng)力降隨地震震級增大而增加?？傊?，應(yīng)力降與震級之間關(guān)系復(fù)雜，認識不一。因而，在討論應(yīng)力降變化時，首先應(yīng)盡可能消除或回避震級的影響。

2.3 應(yīng)力降與震源深度間的關(guān)系

地震釋放的應(yīng)力降隨震中位置及震源深度的不同而有所差異。主震斷層之外區(qū)域地震的應(yīng)力降高于斷層帶內(nèi)，反映出斷裂區(qū)內(nèi)此前較充分的應(yīng)力釋放水平（Jones et al，1998）。從應(yīng)力降隨震源深度的變化來看，通常情況下較深處地震的應(yīng)力降一般較高，如Kanamori（2004）指出，深源地震應(yīng)力降比淺源地震約高10倍。同時，對于淺源地震，應(yīng)力降隨深度增加而增大。Hardebeck等（2009）研究認為，在震源深度為1～7km時，中值應(yīng)力降約為5MPa；7～13km時，約為10MPa；大于13km時，約為50MPa，認為這主要表征脆性地殼內(nèi)部較深處相對較高的應(yīng)力環(huán)境這一特征，即理論上在更深處，斷層承受的正應(yīng)力更大，因此斷層在破裂前可以支撐更大的剪切載荷，當(dāng)?shù)卣鸢l(fā)生時，會有更大的剪應(yīng)力釋放，因而在較深處可以產(chǎn)生更大的應(yīng)力降（Hardebeck et al，1997）。也有研究指出，某些區(qū)域震源深度大于40km的地震的平均應(yīng)力降隨深度的增加有更明顯的增大（Allmann et al，2007），這可用俯沖板塊內(nèi)的剛度變化來進行解釋（Bilek et al，1998），即這些地震的震源主要分布于俯沖帶。也有研究發(fā)現(xiàn)，應(yīng)力降的變化對震源深度的依賴程度隨區(qū)域不同而有差異，如Hardebeck等（1997）對1994年北嶺MW6.7地震序列的279個ML2.5～4.0余震應(yīng)力降的研究表明，當(dāng)震源深度小于15km時，應(yīng)力降上限隨震源深度的增加而增大，應(yīng)力降分布范圍下限則隨深度的變化不明顯；而當(dāng)震源深度大于15km時，應(yīng)力降下限隨震源深度的增加而增大，而上限則不隨深度而變化；他們認為這種現(xiàn)象是由于該區(qū)域存在異常高速層所致。Fletcher等（2006）對幾個地震序列的研究表明，對同一地震序列中震源深度小于主震震源深度的地震事件，其應(yīng)力降隨深度的增加而增大；而對于震源深度大于主震震源深度的地震事件，應(yīng)力降則隨深度的增加而減少。他們認為這是由于應(yīng)力降和介質(zhì)剛度分布主要受破裂帶內(nèi)高速滑移部分的尺度控制所致。就中國大陸地震而言，震源深度小于20km時，中小地震釋放的應(yīng)力降對深度的依賴不顯著，但震源深度大于20km時，應(yīng)力降顯示出隨深度增加而增大的現(xiàn)象（趙翠萍等，2011）。

也有研究認為，應(yīng)力降與震源深度之間不存在統(tǒng)計可信的依賴關(guān)系（華衛(wèi)，2007）。Mori等（2003）對北嶺地震序列余震的研究也未發(fā)現(xiàn)應(yīng)力降對深度明顯的依賴關(guān)系。Allmann等（2009）對1900～2007年全球約2000個5級以上地震的研究顯示，這些中強地震釋放的應(yīng)力降不隨深度的增加而增大。他們認為應(yīng)力降隨深度增加而增大的結(jié)果，可能是由于研究中對地殼淺部采用了常剪切波速這一做法而造成的，由于應(yīng)力降隨剪切波速度的3次方變化，因而在上地殼假設(shè)一個常剪切波速可能會對應(yīng)力降的增加產(chǎn)生明顯的誤差影響（A llmann et al，2007）。

2.4 應(yīng)力降與震源機制間的關(guān)系

地震震源機制一般分為逆斷層、正斷層和走滑斷層等3種類型（惠少興等，2018）。在不同的斷裂方式中，構(gòu)造負荷的速率和摩擦強度不同，因而不同震源機制地震產(chǎn)生的應(yīng)力降可能有所差異（Campbell，1984）。一般而言，逆斷層事件應(yīng)力降大于走滑斷層事件，而走滑斷層事件應(yīng)力降又大于正斷層事件（Mohammadioun，2002），這主要是由于逆斷層事件摩擦強度最大，其剪應(yīng)力最高而正斷層事件最低所致（McGarr，2002，曲均浩等，2014）。但也有完全相反的認識，有研究認為走滑斷層事件應(yīng)力降顯著高于其他類型斷層事件，而正斷層事件應(yīng)力降也高于逆斷層事件（Pérez-Campos et al，2001；Houston，2001），如 Choy等（1995）研究發(fā)現(xiàn)走滑斷層地震應(yīng)力降比其他震源機制地震的應(yīng)力降高1個數(shù)量級，他們用走滑斷層地震比傾滑斷層地震具有更大的破裂速度對此進行解釋，也有觀點認為該現(xiàn)象可能是由于在矩震級相等的情況下，走滑斷層地震的平均拐角頻率比其他斷裂類型的地震更高所致（Allmann et al，2007）。還有研究者認為，基于現(xiàn)有資料，應(yīng)力降與斷層類型之間的依賴關(guān)系缺乏統(tǒng)計顯著性（Scholz et al，1985）；對于中國大陸地區(qū)的研究結(jié)果也表明，應(yīng)力降與斷層類型之間的相關(guān)性不明顯（趙翠萍等，2011）。由此可見，應(yīng)力降與震源機制之間的關(guān)系仍存在較大爭議，有待更系統(tǒng)、深入的分析和研究。

2.5 應(yīng)力降區(qū)域特征及其與構(gòu)造環(huán)境間的關(guān)系

應(yīng)力降有一定的空間分布特征，這是由于不同地區(qū)積累應(yīng)力的能力有所差異（肖蔚文等，1992）。中美洲科科斯俯沖帶是全球地震應(yīng)力降最低的區(qū)域，平均應(yīng)力降小于1MPa；在太平洋蘇拉西威西島北部附近也存在低應(yīng)力降事件。應(yīng)力降顯著偏高的地震事件在興都庫什地區(qū)出現(xiàn)較多，這可能與該地區(qū)多發(fā)生中深源逆斷型地震有關(guān)；在南太平洋貝尼奧夫俯沖帶北端（湯加-克馬德克地區(qū)）也有應(yīng)力降系統(tǒng)地偏高的地震事件分布（Chung et al，1980；Allmann et al，2009）。就中國大陸而言，應(yīng)力降相對較高的區(qū)域主要包括龍門山斷裂帶北東段汶川地震余震區(qū)、青藏高原東北緣以及新疆天山中段地區(qū)，這幾個區(qū)域地震震源機制以逆沖型為主，空間上主要沿斷裂帶集中展布，而其他地區(qū)的應(yīng)力降則相對較小（趙翠萍等，2011）。在更小的空間尺度上，對一個地震序列而言，不同區(qū)域應(yīng)力降同樣有一定差異，華衛(wèi)等（2009）發(fā)現(xiàn)汶川地震后地震應(yīng)力降的空間分布沿龍門山斷裂帶呈一定的分段特征，應(yīng)力水平較低的地段短期內(nèi)發(fā)生強震的可能性較小。周少輝等（2017）對發(fā)生時間相近、震級相當(dāng)、破裂方式相同、大區(qū)域動力環(huán)境相似，而余震活動差別巨大的2014年魯?shù)?.5級、景谷6.6級地震序列震源參數(shù)的研究表明，后續(xù)未發(fā)生強余震的魯?shù)榈卣鹦蛄兄薪^大多數(shù)高應(yīng)力降余震都圍繞主震發(fā)生；而后續(xù)有6級左右強余震發(fā)生的景谷地震序列，在5.8級強余震發(fā)生前，大多數(shù)高應(yīng)力降余震均發(fā)生在未來5.8級強余震震中周圍。

應(yīng)力降與地震破裂過程及構(gòu)造發(fā)育程度有關(guān)，閉鎖、不連續(xù)或年輕構(gòu)造帶上地震的應(yīng)力降一般較高，低應(yīng)力降地震則常出現(xiàn)在曾經(jīng)經(jīng)歷過區(qū)域性地震矩釋放的地區(qū)（相對破碎區(qū)域）（Molnar et al，1972；蔣海昆等，2015）。高應(yīng)力降地震集中分布的區(qū)域有可能是中強地震的潛在震源成核區(qū)（Hardebeck et al，2009）。從已有的研究結(jié)果來看，大洋轉(zhuǎn)換斷層上的地震具有最高的應(yīng)力降，板內(nèi)地震釋放的應(yīng)力降普遍比板間地震高（表1）（Kanamori et al，1975、1986；陳運 泰等，2000；Allmann et al， 2007、 2009；Seno，2014），板間地震應(yīng)力降相對較低可能與俯沖區(qū)域剛度相對較低有關(guān)（Bilek et al，1999）。Molnar等（1972）也指出，印度板塊內(nèi)陸地震應(yīng)力降遠高于板塊邊界的淺源地震事件。Scholz等（1985）進一步指出，板內(nèi)斷層與板間邊界斷層在以下幾個方面明顯不同：①板內(nèi)斷層的滑移速率通常比板間斷層低1～2個數(shù)量級；②板內(nèi)地震總滑移尺度通常為1～10km，而板塊邊界斷層往往為數(shù)百千米；③板內(nèi)斷層在有限的長度上并不連續(xù)，板內(nèi)地震斷層寬度通常只有1～10m，而板塊邊界地區(qū)往往有100～1000m的破碎帶。同時，板內(nèi)地震普遍比板間地震有更長的重復(fù)周期（Kanamori et al，1986），這與板內(nèi)地震較低的滑移速率有關(guān)。實驗研究也表明，巖石摩擦強度隨滑移速率的增加而減少，隨地震重復(fù)周期的增加而增大（Dieterich，1972；Kanamori et al，1986；Goldsby et al，2011），因此，板塊內(nèi)部斷層比板塊邊界有更高的摩擦強度。也有研究指出俯沖區(qū)域剛度較低是產(chǎn)生低應(yīng)力降的重要原因（Bilek et al，1999）。還需要指出的是，流體作用會導(dǎo)致小震應(yīng)力降發(fā)生變化，如對水庫地震與天然地震應(yīng)力降的對比研究顯示，震級較低級時水庫地震應(yīng)力降明顯低于天然地震，但隨震級的增大，兩者逐漸接近，當(dāng)震級大于5時兩者趨于一致（華衛(wèi)等，2012），這是由于流體滲透使得介質(zhì)強度降低、介質(zhì)孔隙壓力增大，進而引起破裂面剪應(yīng)力降低，即在一個較低的應(yīng)力環(huán)境下發(fā)生水庫誘發(fā)地震。與此類似，工業(yè)注水活動也會引起小震應(yīng)力降的變化，如Wu等（2018）對美國俄克拉荷馬州由注水誘發(fā)的4個地震序列的研究表明，每個序列的絕大部分前震和余震應(yīng)力降明顯偏低；對瑞士巴塞爾由于注水采礦誘發(fā)的1000個地震應(yīng)力降的研究結(jié)果顯示，隨著地震空間位置與注水點間的距離由10m增加到300m，應(yīng)力降增高約5倍，這進一步證實了注水可引起孔隙壓力增大，導(dǎo)致低應(yīng)力條件下發(fā)生地震（Goertz-Allmann et al，2011）。

表1 不同構(gòu)造系統(tǒng)地震的中值應(yīng)力降表（據(jù)A llm ann等（2009））

3 地震序列應(yīng)力降在震后地震判定中的應(yīng)用研究

已有研究表明，大震前后中小地震應(yīng)力降變化與地震活動強弱有關(guān)（Chinnery，1964；Richardson et al，1977；Baltay et al，2011）。近年來，隨著數(shù)字地震觀測技術(shù)的發(fā)展，在應(yīng)力降用于震源區(qū)地殼應(yīng)力變化監(jiān)測和序列趨勢判定方面也開展了許多工作。

3.1 應(yīng)力降與地震序列類型判定

有研究顯示，應(yīng)力降與地震序列類型間有一定的關(guān)系，即主-余型地震序列的余震應(yīng)力降相對較高，而震群型地震序列余震的應(yīng)力降則普遍較低（梁秀英，1982）。鐘羽云等（2004）對我國大陸地區(qū)18例強震的靜態(tài)應(yīng)力降與動態(tài)應(yīng)力降之比的研究表明，當(dāng)比值大于0.2時，強震可能為孤立型或近孤立型的主-余型地震，強震后再次發(fā)生較大或更大震級地震的可能性較??；當(dāng)比值小于0.2時，強震可能為多震型或強余震豐富的主-余型地震，強震后再次發(fā)生較大或更大震級地震的可能性較大。此外，陳章立等（2008）也指出，對震群型序列而言，地震應(yīng)力降隨時間而增大可能是前震序列的重要標(biāo)志。還有研究顯示，前震序列的平均應(yīng)力降比余震序列的平均應(yīng)力降大約高1倍以上（秦嘉政等，2005）。

3.2 應(yīng)力降隨時間的變化與序列后續(xù)活動趨勢

從地震序列應(yīng)力降隨時間的變化來看，應(yīng)力降變化趨勢與地震序列的衰減趨勢基本吻合，余震序列應(yīng)力降在主震后一段時間內(nèi)的變化特征為有起伏的逐漸衰減。如華衛(wèi)等（2009）對汶川8.0級地震序列2008年5月12日～6月8日間136次ML3.0地震的研究結(jié)果顯示，汶川地震余震序列中相同震級余震的應(yīng)力降隨時間逐漸減小，這意味著主震后余震區(qū)應(yīng)力水平逐漸降低。王培玲等（2013）計算了玉樹2006年5級震群序列和2010年7.1級地震序列的應(yīng)力降發(fā)現(xiàn)，隨著序列的衰減，應(yīng)力降也逐漸減小。當(dāng)然也有不一致的認識，鄧菲等（2014）對2011年盈江5.8級地震序列中92個2.5級以上地震的應(yīng)力降研究結(jié)果顯示，早期余震密集活動階段的應(yīng)力降比后期余震活動稀疏階段的應(yīng)力降要小，他們認為這種差異可能是由于主震較充分地釋放應(yīng)力，主震后早期余震密集活動階段震源區(qū)應(yīng)力水平不高，而隨著余震活動趨向結(jié)束，后期應(yīng)力水平逐漸恢復(fù)所致。對小震群序列的研究也有類似的認識，如對2003年青島嶗山ML4.1地震序列的研究表明，嶗山序列4.1級主震后，其余震序列的應(yīng)力降遠低于主震，普遍小于5MPa，且約有80%的余震應(yīng)力降小于1MPa，顯示出類似的余震低應(yīng)力降現(xiàn)象（鄭建常等，2008）。鄭建常等（2008）結(jié)合嶗山地震序列活動持續(xù)時間較長（近2年）、主余震震級差較小等特點認為，嶗山序列余震活動屬于低能量的破裂過程，這種余震低應(yīng)力降現(xiàn)象與嶗山序列衰減強度較弱的現(xiàn)象相一致。

部分震例研究顯示，強余震發(fā)生前中小地震應(yīng)力降有一定的上升趨勢，這是介質(zhì)狀態(tài)在高應(yīng)力條件下的反映，而高應(yīng)力背景往往與強余震發(fā)生之間有較密切的關(guān)系。Baltay等（2011）對日本4個地震序列應(yīng)力降的研究表明，在地震以及強余震孕育過程中，孕震區(qū)中小地震應(yīng)力降呈現(xiàn)出逐漸增大的趨勢，震后應(yīng)力降逐漸趨于平穩(wěn)。鄧菲等（2014）對2011年盈江5.8級地震序列應(yīng)力降的研究也顯示，在主震和強余震發(fā)生前應(yīng)力降表現(xiàn)出一定程度的上升趨勢；劉建明等（2016）對2012年新源、和靜交界6.5級地震序列余震的研究表明，總體而言，應(yīng)力降在序列中余震平靜活動階段變化平穩(wěn)，而在強余震前呈現(xiàn)升高—回落的變化過程。王培玲等（2013）計算了玉樹2006年5級震群序列和2010年7.1級地震序列的應(yīng)力降，發(fā)現(xiàn)2次玉樹地震序列的應(yīng)力降在強余震前都出現(xiàn)應(yīng)力降的上升過程。周少輝等（2017）對2014年魯?shù)?.5級、景谷6.6級地震進行對比研究發(fā)現(xiàn)，魯?shù)榈卣鸷罄m(xù)沒有強余震發(fā)生，最大余震震級僅4.8；而景谷6.6級地震后2個月，震源區(qū)先后發(fā)生5.8級、5.9級強余震。周少輝等（2017）的研究還顯示，魯?shù)榈卣鸷笥嗾饝?yīng)力降隨時間幾乎不變，他們認為這由震源區(qū)應(yīng)力釋放已較為充分、主震后震源區(qū)應(yīng)力調(diào)整不顯著所致；而景谷地震序列5.8級強余震發(fā)生前，應(yīng)力降存在明顯的升高過程。還有研究顯示，對于雙震型地震序列，首個強震的應(yīng)力降往往明顯偏低，而第2個強震的應(yīng)力降一般偏大（葉建慶等，2014）。綜上所述，序列地震應(yīng)力降隨時間轉(zhuǎn)折升高的變化過程，可能對后續(xù)強余震的發(fā)生有一定的時間指示意義。

3.3 序列地震應(yīng)力降空間差異與強余震地點判定

地震的發(fā)生與巖石發(fā)生錯動所處的應(yīng)力環(huán)境有關(guān)（杜迎春，2000），震源區(qū)中小地震的應(yīng)力降主要反映震源區(qū)域的應(yīng)力水平，高應(yīng)力降表明該區(qū)域的應(yīng)力水平相對較高，后續(xù)發(fā)生中強震的可能性較大。就強余震而言，可供研究的震例不多，但從已有的研究來看，應(yīng)力降空間差異與強震發(fā)生地點間似乎也有一定關(guān)系，即高應(yīng)力降地震大多緊鄰主震震中或與未來強余震震中相距較近。華衛(wèi)等（2009）研究發(fā)現(xiàn)，汶川地震余震應(yīng)力降空間分布沿龍門山斷裂帶有一定的分段特征，應(yīng)力水平較低的地段（余震釋放的應(yīng)力降相對較低的地段）短期內(nèi)發(fā)生強余震的可能性較小。秦嘉政等（2001）曾對3個月內(nèi)發(fā)生10次5級以上地震的武定6.5級地震序列進行研究，結(jié)果顯示高應(yīng)力降地震主要分布在緊鄰主震與強余震震中的較小范圍內(nèi)，多數(shù)5級以上強余震發(fā)生在應(yīng)力降等值線變化明顯的區(qū)域。對后續(xù)沒有強余震發(fā)生的魯?shù)榈卣鹦蛄信c先后發(fā)生5.8級、5.9級強余震的景谷地震序列的對比研究表明，魯?shù)榈卣鹦蛄兄薪^大部分高應(yīng)力降余震圍繞主震震中分布；而景谷地震序列中，在5.8級強余震發(fā)生前，絕大多數(shù)高應(yīng)力降余震發(fā)生在5.8級強余震震中周圍（周少輝等，2017）。

4 應(yīng)力降計算結(jié)果不確定性討論

應(yīng)力降計算涉及多種可能的影響因素，如震源模型選擇、震源輻射因子、研究區(qū)域介質(zhì)密度、零頻極限和拐角頻率等。此外，計算方法及模型參數(shù)選擇對最終結(jié)果也有影響（Ide et al，2001；蔣海昆等，2015）。

對中小地震，一般選擇點源模型，選取同一計算模型，對各個單一地震應(yīng)力降相對大小的比較并無影響。點源模型條件下應(yīng)力降計算公式為

式中，ρ為研究區(qū)域介質(zhì)密度；Rθφ為 SH波的震源輻射因子；Ω0、fc分別為零頻極限、拐角頻率。

地震震源譜一般由多臺平均得到，這在一定程度上減弱了震源輻射因子對計算結(jié)果的影響，且震源輻射因子的取值只有當(dāng)臺站非常接近S波節(jié)線或臺站恰好處于S波節(jié)線位置上時才會產(chǎn)生嚴重的誤差（康英等，2004；Stork，2004），因此震源輻射因子對應(yīng)力降計算結(jié)果的影響可忽略不計。

關(guān)于介質(zhì)密度ρ對應(yīng)力降的可能影響，根據(jù)對滇西地殼三維密度結(jié)構(gòu)（楊文采等，2015）及魯?shù)檎饏^(qū)地殼平均密度的研究（楊光亮等，2014）可知，云南魯?shù)椤⒕肮鹊貐^(qū)平均密度最大不確定性約為10%，基于式（9），當(dāng)不考慮其他因素的影響時，其對應(yīng)力降不確定性的影響也約為10%，但對魯?shù)榈卣?、景谷地震序列各震級段余震?yīng)力降隨時間的變化，即使考慮最大10%的應(yīng)力降不確定性，也遠小于實際應(yīng)力降隨時間的變化，因此，區(qū)域密度不確定性所導(dǎo)致的應(yīng)力降不確定性遠小于實際應(yīng)力降隨時間的變化，據(jù)此，亦可粗略忽略區(qū)域密度不確定性對應(yīng)力降的影響（周少輝等，2017）。

關(guān)于拐角頻率fc和零頻極限Ω0對應(yīng)力降計算結(jié)果的影響，由于中小地震的拐角頻率往往落在儀器帶寬的邊緣附近（Hough et al，1995），而獲得拐角頻率和零頻極限在大多數(shù)情況下依賴于人工判斷，有很強的主觀性（劉杰等，2003），同時在地震波傳播過程中，介質(zhì)對高頻成分存在較大的衰減作用，對同一地震，若接收點位于地震波輻射的不同方位，則接收到的頻率也會不同（傅承義等，1985），因而，具體到頻譜圖上，拐角頻率和零頻極限的測定往往就有很大的不確定性（吳忠良等，1999）。更為復(fù)雜的是，在實際研究中，很少能夠判定拐角頻率的變化是由測量誤差引起的，還是由震源性質(zhì)變化引起的（Sonley et al，2006；Prieto et al，2007）。因此，對震源譜拐角頻率fc和零頻極限Ω0的準確識別本質(zhì)上非常困難。此外，由震源譜擬合計算拐角頻率fc和零頻極限Ω0的過程中，還存在噪聲項、儀器項、傳播路徑效應(yīng)和臺站的局部場地效應(yīng)等諸多影響因素，若考慮誤差傳遞的話，消除這些影響因素過程中所導(dǎo)致的對應(yīng)力降計算結(jié)果的影響究竟如何也尚不明晰。一般認為：① 由于所用到的臺站地震儀一般為同類儀器，因而儀器響應(yīng)對記錄信號的影響是固定、一致的，短期內(nèi)不隨時間而改變；②由于臺站環(huán)境短期內(nèi)不會有大的改變，因而臺站周圍環(huán)境噪聲對結(jié)果無太大影響；③一般采用多臺多震聯(lián)合反演方法計算介質(zhì)品質(zhì)因子Q值，以此消除傳播路徑的影響（Atkinson et al，1992、1995；黃玉龍等，2003）。由于不同構(gòu)造區(qū)具有不同的 Q值，因此應(yīng)保證臺網(wǎng)射線覆蓋區(qū)屬于同一個構(gòu)造區(qū)域，否則，結(jié)果誤差可能很大。對地震序列應(yīng)力降研究而言，相對于較大的震中距，所有地震均局限在一個較小的空間范圍內(nèi)，這可以確保在同一個構(gòu)造區(qū)域里計算Q值；場地響應(yīng)與近表面地層密度及地震波速度等有關(guān)，與地震波傳播距離無關(guān)（Shearer，1999）。目前大多采用多臺多震聯(lián)合反演的方法求取臺站場地響應(yīng)；而在此之前采用多臺多震聯(lián)合反演方法計算非彈性衰減系數(shù)并采用三段幾何衰減模型計算幾何衰減（Atkinson et al，1992、1995；黃玉龍等，2003），最終場地響應(yīng)取多震平均結(jié)果，因而也可認為對結(jié)果無顯著影響（劉杰等，2003）?？梢姡瑥目赡苡绊懝战穷l率、零頻極限測量的各項因素來看，每個單項因素的影響在一定的假設(shè)前提下似乎都可忽略。然而，周少輝等（2017）以魯?shù)?、景谷地震余震序列為例并結(jié)合數(shù)值實驗所進行的研究則得到不一致的認識。他們選擇某一較小震級范圍（如2.5≤ML≤2.9）內(nèi)的地震進行研究，以盡可能消除震級對拐角頻率的影響（Savage，1972；華衛(wèi)，2007），并將該震級范圍內(nèi)所有地震事件拐角頻率的平均標(biāo)準偏差粗略地作為每次地震拐角頻率的測量誤差，由式（9）求Δσ對fc的偏微分，當(dāng)不考慮其他因素的影響時，即可獲得與拐角頻率測定誤差相對應(yīng)的應(yīng)力降計算的不確定性。結(jié)果顯示，拐角頻率測量誤差對應(yīng)力降計算結(jié)果產(chǎn)生了顯著的影響，誤差棒幾乎完全覆蓋相鄰地震事件應(yīng)力降隨時間的變化。對零頻極限開展類似的研究，也得到同樣的認識。從這個角度來看，前述3.2節(jié)根據(jù)應(yīng)力降隨時間變化判定震源區(qū)后續(xù)應(yīng)力狀態(tài)、開展強余震預(yù)測研究的認識似乎也存在較大的不確定性。

綜上所述，輻射因子與區(qū)域介質(zhì)密度對應(yīng)力降隨時間變化趨勢的誤差影響可暫且忽略。對于獲取拐角頻率和零頻極限過程中可能產(chǎn)生誤差的各項影響因素，雖然分開來看，在一定假設(shè)前提下各項因素似乎都可忽略，但數(shù)值實驗結(jié)果卻顯示，拐角頻率、零頻極限測量的綜合誤差對應(yīng)力降計算結(jié)果會產(chǎn)生顯著影響，誤差棒幾乎可以完全覆蓋余震應(yīng)力降隨時間的變化趨勢。從這個角度來看，根據(jù)應(yīng)力降的時空變化判定震后地震活動趨勢的認識似乎也存在較大的不確定性。

5 討論及結(jié)論

（1）地震應(yīng)力降一般為零點幾兆帕至數(shù)十兆帕，1～3MPa為優(yōu)勢分布范圍。中國大陸地區(qū)的中小地震應(yīng)力降為0.1～20.0MPa，絕大多數(shù)小于10MPa。關(guān)于地震應(yīng)力降與地震震級之間的關(guān)系至今仍未有統(tǒng)一的認識，常應(yīng)力降模型認為，應(yīng)力降不隨地震震級的變化而變化，但事實上“大”地震與“小”地震情況有別，對于“大”地震，應(yīng)力降宏觀來看不隨地震震級變化而改變；而對于“小”地震，應(yīng)力降隨震級（或地震矩）的增加而增大。從應(yīng)力降隨深度的變化來看，理論上在更深處，斷層承受的正應(yīng)力更大，因此斷層在破裂前可以支撐更大的剪切載荷，地震發(fā)生時會有更大的剪應(yīng)力釋放，因而在較深處可以產(chǎn)生更大的應(yīng)力降；然而有研究認為，應(yīng)力降對深度不存在明顯的依賴關(guān)系，應(yīng)力降隨深度增加而增大的研究結(jié)果可能是由于研究中對地殼淺部采用常剪切波速所致。對于不同的斷裂方式，有研究認為，由于逆斷層事件摩擦強度最大，其剪應(yīng)力最高而正斷層事件最低，因此逆斷層事件應(yīng)力降大于走滑事件，而走滑事件應(yīng)力降又大于正斷事件；但也有研究認為，走滑斷層事件應(yīng)力降顯著高于其他類型斷層，而正斷層事件應(yīng)力降也高于逆斷層事件應(yīng)力降，他們認為這種現(xiàn)象可能是由于矩震級相等時，走滑型地震的平均拐角頻率比其他斷裂類型地震更高所致。

應(yīng)力降空間分布存在一定的區(qū)域特征。沿中美洲科科斯俯沖帶與太平洋蘇拉西威西島北部附近為全球地震應(yīng)力降最低的區(qū)域，全球最高應(yīng)力降出現(xiàn)在興都庫什地區(qū)，南太平洋貝尼奧夫帶的俯沖帶北端（湯加-克馬德克地區(qū)）有系統(tǒng)性較高的應(yīng)力降分布。就更小空間尺度而言，對同一地震序列，高應(yīng)力降地震主要在靠近主震與強余震震中的范圍內(nèi)分布。應(yīng)力降與地震破裂過程及構(gòu)造發(fā)育程度有關(guān)，閉鎖、不連續(xù)或年輕構(gòu)造帶上地震的應(yīng)力降一般較高，低應(yīng)力降地震則常出現(xiàn)在曾經(jīng)歷過區(qū)域性地震矩釋放的地區(qū)（相對破碎區(qū)域）。水庫地震的應(yīng)力降明顯小于天然地震。板內(nèi)地震釋放的應(yīng)力降普遍比板間地震高，板間地震應(yīng)力降相對較低可能與俯沖區(qū)域剛度相對較低有關(guān)。

總之，應(yīng)力降與震級、震源深度、構(gòu)造系統(tǒng)、震源機制、序列類型等之間的關(guān)系復(fù)雜，認識不一。這種差異可能主要是因為大部分研究結(jié)果都是針對某個特定區(qū)域的部分地震序列的，而地區(qū)差異、序列類型差異及同一區(qū)域地震事件發(fā)生時間的不同都可能對結(jié)果產(chǎn)生很大影響。因而基于現(xiàn)有資料，應(yīng)力降與上述各因素之間的關(guān)系缺乏統(tǒng)計顯著性檢驗，仍待更系統(tǒng)、深入地研究應(yīng)力降與上述因子之間的關(guān)系。

（2）應(yīng)力降在震后趨勢判定中可以發(fā)揮一定的作用。從現(xiàn)有有限的研究結(jié)果來看，主-余型地震序列的余震應(yīng)力降相對較高，而震群型序列余震的應(yīng)力降則普遍較低，前震序列的平均應(yīng)力降比余震序列大約高1倍以上。當(dāng)強震的靜態(tài)應(yīng)力降與動態(tài)應(yīng)力降之比大于0.2時，可能為孤立型或近孤立型的主-余型地震；當(dāng)比值小于0.2時，強震可能為多震型或強余震豐富的主-余型地震。

應(yīng)力降的時空變化對于后續(xù)地震活動趨勢判斷有一定指示意義，盡管研究樣本有限，但現(xiàn)有研究還是顯示，消除震級影響之后，小震應(yīng)力降在強余震發(fā)生前大多存在一定的上升趨勢，高應(yīng)力降背景往往與強余震發(fā)生之間有較密切的關(guān)系。因而，震后應(yīng)力降隨時間持續(xù)升高—快速轉(zhuǎn)折的變化對后續(xù)中強震可能有一定的時間指示意義。同時，應(yīng)力降空間變化與強震發(fā)生地點間似乎也有一定的相關(guān)性，高應(yīng)力降地震主要在靠近主震與強余震震中的范圍內(nèi)分布。震源區(qū)中小地震應(yīng)力降主要反映震源區(qū)域的背景應(yīng)力水平，高值應(yīng)力降表明該區(qū)域的應(yīng)力水平相對較高，后續(xù)發(fā)生中強震的可能性較大。

（3）諸多因素會影響應(yīng)力降的計算。針對中小地震，選擇布倫圓盤模型之后，破裂模型對所有地震均具有一致的影響，可不再考慮。從圓盤模型應(yīng)力降計算公式出發(fā)，簡化條件下，當(dāng)不考慮其他因素的影響時，1%的拐角頻率變化可引起3%的應(yīng)力降變化；而1%的區(qū)域介質(zhì)密度、零頻極限和輻射因子變化，均可引起1%的應(yīng)力降變化。具體來看，震源輻射因子對應(yīng)力降計算結(jié)果的影響可忽略不計。對研究區(qū)域密度不確定性所導(dǎo)致的應(yīng)力降誤差而言，即使考慮區(qū)域密度最大不確定性，其導(dǎo)致的應(yīng)力降誤差也并未能覆蓋實際應(yīng)力降隨時間的變化，因此區(qū)域密度不確定性對應(yīng)力降的影響也可暫不考慮。拐角頻率和零頻極限由震源譜擬合計算得到，在計算過程中存在諸多影響因素，如噪聲項、儀器項、傳播路徑效應(yīng)和臺站的局部場地效應(yīng)等。從單項來看，每項因素的影響在一定的假設(shè)前提下似乎均可忽略，但事實上，數(shù)值實驗結(jié)果卻顯示，拐角頻率、零頻極限的綜合誤差對應(yīng)力降計算結(jié)果有顯著影響，誤差棒幾乎完全覆蓋余震應(yīng)力降隨時間的變化。從這個角度來看，根據(jù)應(yīng)力降時空變化判定震后趨勢的認識似乎也存在較大的不確定性。因此，應(yīng)力降計算過程中各項因素對結(jié)果的影響，是基于應(yīng)力降變化的各項后續(xù)研究的基礎(chǔ)。

（4）需要指出的是，從現(xiàn)狀來看，對于應(yīng)力降應(yīng)用于地震預(yù)測或序列趨勢判定的工作，已有的震例研究不多，且絕大多數(shù)缺乏必要的統(tǒng)計檢驗和震例驗證，從預(yù)測的角度，這不能不說是一個缺憾。許多研究更是從物理概念或模型出發(fā)，認為什么樣的應(yīng)力降變化“應(yīng)該”與高的或低的應(yīng)力狀態(tài)相對應(yīng)，其中包含了諸多的人為假設(shè)，但事實是否如此，尚無定論甚至還有爭議。簡言之，由于受研究深入程度有限、統(tǒng)計樣本數(shù)偏少、模型及參數(shù)差異等眾多因素的影響，所得到的認識在外推運用時需要慎重對待，對同一地區(qū)、相同震級范圍、不同時期的地震的應(yīng)力降進行對比，對實際序列趨勢判定工作可能更為重要（蔣海昆等，2015）。此外，震源參數(shù)計算的誤差分析，也是一個需要認真討論的基礎(chǔ)問題。