李金磊，萬永革.希臘克里特MW6.0地震后的應(yīng)力方向變化與強(qiáng)余震發(fā)生［J］.地震工程學(xué)報(bào)，2024，46（2）：491500.DOI：10.20000j.10000844.20220929001

摘要：

地震的震源機(jī)制是地殼應(yīng)力變化的指示器，而地殼應(yīng)力變化與強(qiáng)震的發(fā)生直接相關(guān)。前人研究了地震震源機(jī)制變化在視應(yīng)力較高的走滑型大震前的應(yīng)力變化過程，而未見到震源機(jī)制變化對視應(yīng)力較低的正斷型大震發(fā)生的指示作用的研究。文章以2021年希臘克里特MW6.0正斷型地震序列為例，通過計(jì)算地震序列震源機(jī)制解與區(qū)域應(yīng)力場方向之間最小空間旋轉(zhuǎn)角的變化，揭示應(yīng)力變化與強(qiáng)震發(fā)生的關(guān)系。為保證震源機(jī)制解的準(zhǔn)確性，采用多家機(jī)構(gòu)確定的震源機(jī)制得到中心震源機(jī)制作為該地震的震源機(jī)制，而后采用該地震序列精確的震源機(jī)制求解當(dāng)?shù)貞?yīng)力場，最后計(jì)算地震震源機(jī)制與主震震源機(jī)制及與所估計(jì)的地殼應(yīng)力場方向的空間旋轉(zhuǎn)角隨時(shí)間的變化，探索強(qiáng)震發(fā)生與應(yīng)力場變化的關(guān)系。結(jié)果表明：在主震發(fā)生的短期內(nèi)，余震震源機(jī)制與該區(qū)域應(yīng)力場方向的空間旋轉(zhuǎn)角較大，與其后小震級(jí)的弱地震活動(dòng)對應(yīng);隨后余震震源機(jī)制與應(yīng)力場方向的空間旋轉(zhuǎn)角減小，對應(yīng)后面發(fā)生的3次MWgt;5.0的強(qiáng)余震，在此之后的長時(shí)間內(nèi)余震震源機(jī)制和應(yīng)力場方向的空間旋轉(zhuǎn)角再次增大，對應(yīng)的余震震級(jí)及頻度皆明顯下降。文章以2021年希臘克里特MW6.0地震序列為例，發(fā)現(xiàn)視應(yīng)力較低的正斷型地震前也存在應(yīng)力方向集中現(xiàn)象，為探索地震應(yīng)力前兆提供了范例。

關(guān)鍵詞：

希臘克里特地震序列;震源機(jī)制中心解;構(gòu)造應(yīng)力場;最小空間旋轉(zhuǎn)角

中圖分類號(hào)：P319.56文獻(xiàn)標(biāo)志碼：A文章編號(hào)：10000844（2024）02-0491-10

DOI：10.20000j.10000844.20220929001

0引言

地震的發(fā)生與地殼應(yīng)力狀態(tài)有著直接的因果關(guān)系，而震源機(jī)制解直觀地反映了地震破裂幾何特征和運(yùn)動(dòng)學(xué)特征，是研究區(qū)域構(gòu)造應(yīng)力的基礎(chǔ)［14］。從20世紀(jì)90年代開始，地震學(xué)家們就嘗試從地震序列、震源機(jī)制解之間尋找一定的聯(lián)系，來促進(jìn)人們對發(fā)震機(jī)理和地震預(yù)測等研究領(lǐng)域的了解［57，10］。在唐山地震序列應(yīng)力釋放過程研究中，李欽祖等［5］根據(jù)唐山地震序列震源機(jī)制數(shù)據(jù)對其應(yīng)力釋放調(diào)整機(jī)制進(jìn)行了分析，但對于這一過程中強(qiáng)余震震源機(jī)制與應(yīng)力場的分布關(guān)系并未進(jìn)行詳細(xì)說明。王俊國等［6］利用Harvarduniversity測定的千島島弧地區(qū)地震的矩心矩張量（CMT）解，分析該地區(qū)震源機(jī)制一致性特征，結(jié)果表明大震前一致性參數(shù)重復(fù)出現(xiàn)低值。在研究1992年美國蘭德斯MW7.3地震事件和1999年美國赫克托礦MW7.1地震事件時(shí)，萬永革［7］利用Kagan［8］提出的雙力偶地震震源機(jī)制空間旋轉(zhuǎn)以及Hauksson［9］對美國南加州地區(qū)小震震源機(jī)制的詳細(xì)數(shù)據(jù)，通過對比分析獲得了該區(qū)域震源機(jī)制和主震震源機(jī)制之間10年內(nèi)的平均空間旋轉(zhuǎn)角，并發(fā)現(xiàn)2次大震前小震均傾向于接近主震震源機(jī)制。韓曉明等［10］利用美國南加州地區(qū)1981—2011年間的地震震源機(jī)制解資料，研究發(fā)現(xiàn)小震發(fā)震應(yīng)力場P軸向區(qū)域構(gòu)造應(yīng)力場主壓力的趨近現(xiàn)象。

前述研究的強(qiáng)震前震源機(jī)制所揭示的地震孕育過程均對應(yīng)視應(yīng)力較高的走滑型地震，而對于視應(yīng)力較低的正斷型強(qiáng)震前是否有類似的現(xiàn)象，目前未見報(bào)道。2021年9月27日6時(shí)17分21秒（當(dāng)?shù)貢r(shí)間）在希臘克里特發(fā)生MW6.0正斷型地震。震中位于35.24°N，25.27°E，震源深度6km。隨后發(fā)生多次強(qiáng)余震，并且均可得到多個(gè)機(jī)構(gòu)給出的震源機(jī)制解。這為得到該地震發(fā)生后震源機(jī)制變化和余震發(fā)生的關(guān)系創(chuàng)造了良好的契機(jī)。本文以該正斷型地震序列為實(shí)例進(jìn)行研究，首先，為降低各機(jī)構(gòu)震源機(jī)制解對于分析結(jié)論不確定性的影響，采用萬永革［11］提出的震源機(jī)制中心解的方法，運(yùn)用LevenbergMarquardt方案獲得主震和余震震源機(jī)制中心解；然后根據(jù)得到的震源機(jī)制中心解，求得地震序列的局部應(yīng)力場，以備求解余震震源機(jī)制解和局部應(yīng)力場最優(yōu)滑動(dòng)節(jié)面的最小空間旋轉(zhuǎn)角；最后，采用震源機(jī)制最小空間旋轉(zhuǎn)角計(jì)算方法，得出該區(qū)域各時(shí)間段內(nèi)震源機(jī)制中心解與該地區(qū)應(yīng)力場節(jié)面最小空間旋轉(zhuǎn)角的變化，以及主震震源機(jī)制和余震震源機(jī)制的最小空間旋轉(zhuǎn)角隨時(shí)間的變化，以此來分析地震發(fā)生與應(yīng)力場方向集中的關(guān)系。

1震源機(jī)制中心解的求解

克里特島位于希臘島弧前緣，是歐亞板塊與非洲板塊相互作用的結(jié)果，構(gòu)造變形、地震活動(dòng)極其頻繁［1213］。該區(qū)域地質(zhì)構(gòu)造環(huán)境復(fù)雜，含有許多微板塊運(yùn)動(dòng)及區(qū)域性尺度結(jié)構(gòu)。本次地震余震較多，有網(wǎng)站（如EMSC網(wǎng)站，https：www.emsccsem.org#2）提供了這次地震的主震震源機(jī)制，以及16組余震的震源機(jī)制?？紤]到該地區(qū)地震完整性，以MW3.0為所采用的地震序列目錄震級(jí)下限。文中所用的地震序列目錄如表1所列，包括1個(gè)MW6.0以上的地震事件、3個(gè)MW5.0～6.0地震事件和13個(gè)MW5.0以下的地震事件（圖1）。

1.1震源機(jī)制中心解求解

震源機(jī)制是地殼應(yīng)力場研究的關(guān)鍵基礎(chǔ)資料，但由于各個(gè)機(jī)構(gòu)采用的計(jì)算方法及參考資料各不相同，導(dǎo)致同一次地震震源機(jī)制離散性較大。為了找到最可靠的震源機(jī)制解，采用萬永革［11］提出的震源機(jī)制中心解模型來分析地震震源機(jī)制。首先設(shè)定某震源機(jī)制為起始值，對該值域進(jìn)行一階泰勒擴(kuò)展;然后采用LevenbergMarquardt迭代法求出，同時(shí)逐步調(diào)整初始值;最后，將與各備選震源機(jī)制最小空間旋轉(zhuǎn)角平方和最小的解答作為此次地震的中心解。

1.2主震震源機(jī)制中心解的求解

整理主震數(shù)據(jù)并列于表2。震源機(jī)制中心解（表2第3列）和相應(yīng)的起始解標(biāo)準(zhǔn)差（表2第4列）是將每個(gè)機(jī)構(gòu)確定的希臘克里特MW6.0地震震源機(jī)制解（表2第2列）作為起始解計(jì)算出來的。經(jīng)過比較后得出：無論使用哪家機(jī)構(gòu)得到的震源機(jī)制結(jié)果，差別都很小，對應(yīng)標(biāo)準(zhǔn)差只發(fā)生在小數(shù)點(diǎn)之后5位處。結(jié)果表明，用這種方法所獲得的震源機(jī)制的中心解是穩(wěn)定可靠的。即便如此，我們依然選擇標(biāo)準(zhǔn)差最小的震源機(jī)制解作為最終解。對于這次地震，我們發(fā)現(xiàn)用雅典國家天文臺(tái)（NationalObservatoryofAthens，NOA）測量的震源機(jī)制解為起始解所得到的標(biāo)準(zhǔn)偏離值最小，將此中心解（第一個(gè)節(jié)面走向、傾角和滑動(dòng)角的值分別是：32.05°，40.91°，-92.12°;另一個(gè)節(jié)面走向、傾角和滑動(dòng)角的值分別是：214.85°，49.13°，-88.16°）當(dāng)做最終結(jié)果（表2第8行）。根據(jù)走向、傾角和滑動(dòng)角與壓應(yīng)力軸、張應(yīng)力軸、中間軸間的關(guān)系［14］，可以得到壓應(yīng)力軸走向和傾伏角分別為142.28°，85.66°，不確定范圍為46.08°～296.04°和82.87°～103.97°;張應(yīng)力軸走向和傾伏角分別為303.55°，4.11°，不確定范圍為285.02°～322.02°和-14.41°～22.57°;中間軸走向和傾伏角分別為33.65°，1.39°，不確定范圍為15.12°～52.12°和-10.58°～13.40°。

由于震源機(jī)制P、T、B軸的二義性，通過4種坐標(biāo)軸旋轉(zhuǎn)方式可以完全由一個(gè)震源機(jī)制的P、T、B軸旋轉(zhuǎn)到另外一個(gè)震源機(jī)制的P、T、B軸上。這4種旋轉(zhuǎn)角度中最小的一個(gè)可以代表這兩個(gè)震源機(jī)制的差別，通常稱為最小空間旋轉(zhuǎn)角［11］。表2第5列給出了不同機(jī)構(gòu)發(fā)布的震源機(jī)制和作者測定的中心解之間的最小空間旋轉(zhuǎn)角。經(jīng)比較可知，此震源機(jī)制和中心解的最小空間旋轉(zhuǎn)角的最小值為7.43°，最大值為31.36°，最小三維空間旋轉(zhuǎn)角的標(biāo)準(zhǔn)差最小值為19.27°。必須指出，盡管由NOA提出的震源機(jī)制解作為起始解得到的實(shí)際中心解的標(biāo)準(zhǔn)差為最小值，但這并不代表此解離實(shí)際震源機(jī)制中心解之間的空間距離最近，而本次地震離實(shí)際中心解較近者應(yīng)是GCMT（GlobalCentroidMomentTensor）。這一現(xiàn)象的根源在于，求解中心解是一個(gè)具有非線性性質(zhì)的問題，用雅克比矩陣求解各參數(shù)偏量時(shí)有可能超過或達(dá)不到最優(yōu)解（表2第3列），這是本文利用各機(jī)構(gòu)實(shí)測震源機(jī)制解反演初值的基礎(chǔ)。將最終確定的中心震源機(jī)制解的不確定性以及空間三維輻射花樣繪制于圖2。圖2（a）中黑色曲線為震源機(jī)制中心解兩個(gè)不同節(jié)面，綠色曲線所含區(qū)域范圍為其不確定區(qū)間;中心震源機(jī)制解壓縮軸、拉張軸與中間軸則分別以紅、藍(lán)、黃點(diǎn)表示，四周相應(yīng)顏色的圓圈是其不確定度區(qū)間，綠色、黑色和藍(lán)色的小點(diǎn)代表每個(gè)作者獲得的震源機(jī)制解的壓縮軸、張力軸和中間軸的分布;紫色弧是每個(gè)作者獲得的震源機(jī)制解節(jié)面。圖2（b）中的壓縮范圍為藍(lán)色，拉張范圍為紅色。

1.3余震震源機(jī)制中心解的求解、分類及總體特征分析

16組余震震源機(jī)制中心解匯總于表3，其各余震最后確定的中心震源機(jī)制解及其不確定性繪于圖3中，空間三維輻射花樣繪于圖4。地震震源機(jī)制分類對于地震動(dòng)力學(xué)分析有著重要作用。采用萬永革［11，15］提出的震源機(jī)制水平應(yīng)變花面應(yīng)變的地震震源機(jī)制分類方法。圖5表示了2021年希臘克里特地震序列的余震震源機(jī)制分類情況，其中正斷型地震有11個(gè)，占總體數(shù)量的68.75%;正走滑型地震有3個(gè)，占總體數(shù)量的18.75%;走滑型地震有1個(gè)，占總體數(shù)量的6.25%;逆走滑型地震有1個(gè)，占總體數(shù)量的6.25%。再考慮到主震為正斷型地震，總體來講該地震序列的震源機(jī)制以正斷型和正走滑型為主，應(yīng)屬于正斷型的構(gòu)造應(yīng)力背景。而研究區(qū)東北部的Kastelli斷層（圖1）也是NNESSW為走向的正斷層體系［16］，佐證了該地區(qū)為正斷型應(yīng)力體系。

2應(yīng)力場方向和相對大小的求解

地震的孕育和發(fā)生是地殼內(nèi)應(yīng)力、應(yīng)變長期積累與釋放的結(jié)果，通過運(yùn)用數(shù)值模擬技術(shù)對現(xiàn)代應(yīng)力場的產(chǎn)生和發(fā)展過程進(jìn)行系統(tǒng)分析研究，也是探究地震產(chǎn)生機(jī)理的一種有效手段。為解決研究區(qū)域應(yīng)力場的方向及其相對大小問題，本文采用萬永革等［14，1719］基于斷層面上應(yīng)力張量剪應(yīng)力方向符合斷層面上滑動(dòng)方向以及斷層滑動(dòng)的性質(zhì)，所提出的擦痕數(shù)據(jù)結(jié)合定性斷層滑動(dòng)數(shù)據(jù)及相應(yīng)位置在應(yīng)力場按照搜索1°的方向間隔和0.1的應(yīng)力比值間隔網(wǎng)格搜索的方法，為反演參數(shù)提供一定置信度的置信區(qū)三角形的三個(gè)邊分別為震源機(jī)制P，T和B軸傾伏角的刻度，三角形中的點(diǎn)劃線表示震源機(jī)制類型分類的界限。相對面應(yīng)變（As）以圖下方的色棒為標(biāo)準(zhǔn)采用背景顏色繪制。白線為網(wǎng)格線;SS：走滑型；NS：正走滑型;N：正斷型;RS：逆走滑型;R：逆斷型間。利用先前求解得到的主震及余震序列的震源機(jī)制中心解（表2第8行第3列、表3第3列），采用網(wǎng)格搜索法對應(yīng)力場進(jìn)行求解。應(yīng)力場反演如圖6所示，從中可以看出在95%置信度下，最優(yōu)解的置信范圍都較小，說明以圖中結(jié)果來反映此區(qū)域應(yīng)力場是較為準(zhǔn)確的。最終求得該地區(qū)應(yīng)力場最優(yōu)節(jié)面參數(shù)：走向293°，傾角67°，滑動(dòng)角-57°。根據(jù)走向、傾角、滑動(dòng)角與壓縮軸、張力軸、中間軸間的關(guān)系［19］，

（a）圖中黑弧代表選定的“可能斷層面”施密特投影，球外藍(lán)色箭頭即S1軸的水平方向，球內(nèi)藍(lán)色小箭頭即“可能斷層面\"觀測錯(cuò)動(dòng)方位。球體外側(cè)面紅色箭頭為S3軸水平方向，內(nèi)側(cè)面紅色小箭頭為“可能斷層面”理論錯(cuò)動(dòng)方向。綠色曲線表示置信度95%應(yīng)力場中最大剪應(yīng)力節(jié)面，黃色箭頭表示該節(jié)面的錯(cuò)動(dòng)方向，紅色、黃色和藍(lán)色封閉曲線分別表示主壓應(yīng)力軸、中應(yīng)力軸和主張應(yīng)力軸95%置信度下的置信區(qū)間。（b）圖為R=0.9情況下的應(yīng)力輻射花樣圖，紅色代表壓縮軸，藍(lán)色代表拉伸軸得到壓縮軸走向和傾伏角分別為245.98°，55.45°，不確定范圍分別為242.95°～246.98°和54.47°～55.95°;張力軸走向和傾伏角分別為359.77°，15.52°，不確定范圍分別為358.77°～360.77°和15.02°～16.02°;中間軸走向和傾伏角分別為99.00°，30.00°，不確定范圍分別為98.00°～100.00°和29.50°～30.50°。其中，應(yīng)力形因子R=（S2-S3）（S1-S3），這里S1、S2和S3是在主軸坐標(biāo)系下的三個(gè)本征值，分別代表主張應(yīng)力、中間應(yīng)力和主壓應(yīng)力的相對大小，拉張為正。

3最小空間旋轉(zhuǎn)角計(jì)算結(jié)果與分析

通過比較各震源機(jī)制中心解與該地區(qū)應(yīng)力場節(jié)面的最小空間旋轉(zhuǎn)角，可以定量得出震源機(jī)制中心解間的差別——其值越小表明此地震與該地區(qū)應(yīng)力場節(jié)面的差別越小，反之亦然。

按照每一個(gè)地震的發(fā)震時(shí)間，將各余震震源機(jī)制中心解（表3第3列）、主震震源機(jī)制中心解（表2第3列）與應(yīng)力場最優(yōu)節(jié)面解（走向293°，傾角67°，滑動(dòng)角-57°）之間的最小空間旋轉(zhuǎn)角整合成表4。圖7（a）顯示了每個(gè)余震震源機(jī)制中心解與主震震源機(jī)制中心解之間最小空間旋轉(zhuǎn)角的關(guān)系變化。

圖7（b）為各震源機(jī)制中心解與該地區(qū)應(yīng)力場節(jié)面的最小空間旋轉(zhuǎn)角關(guān)系變化。圖7（c）為希臘克里特島MW6.0主震和余震地震序列的MT圖。

通過對數(shù)據(jù)的比較，將所有余震活動(dòng)分成三個(gè)階段：在主震發(fā)生后的48h內(nèi)，各個(gè)余震震源機(jī)制的中心解與該區(qū)域應(yīng)力場節(jié)面的最小空間旋轉(zhuǎn)角明顯較大，且明顯靠近主震震源機(jī)制中心解，在達(dá)到1次震級(jí)為MW5.3（表1第8行）的強(qiáng)余震后，余震數(shù)量及震級(jí)明顯降低，其震源機(jī)制中心解與該區(qū)域應(yīng)力場節(jié)面的最小空間旋轉(zhuǎn)角出現(xiàn)減小的趨勢；主震發(fā)生一個(gè)月后，余震又呈集群式出現(xiàn)，此時(shí)余震的數(shù)量及出現(xiàn)的最大震級(jí)均小于第一階段，且各個(gè)余震震源機(jī)制的中心解與該區(qū)域應(yīng)力場節(jié)面的最小空間旋轉(zhuǎn)角明顯減小;主震發(fā)生三個(gè)月后，余震再次集群出現(xiàn)，但此次僅在數(shù)量上少于前兩次，震級(jí)整體較前兩次高，出現(xiàn)2次MW5.0（表1第16、17行）以上的強(qiáng)余震，最高達(dá)MW5.7，其震源機(jī)制最小空間旋轉(zhuǎn)角與該區(qū)域應(yīng)力場節(jié)面的最小空間旋轉(zhuǎn)角明顯增大;直到2022年9月該地區(qū)未再發(fā)生震級(jí)達(dá)MW5.0以上的強(qiáng)余震。而這三個(gè)階段可以按照強(qiáng)震前的地震發(fā)生與地殼應(yīng)力場的一致性相關(guān)來解釋。主震后48h內(nèi)發(fā)生的地震震源機(jī)制與地殼應(yīng)力場相差較大（最小空間旋轉(zhuǎn)角較大），預(yù)示著震源機(jī)制發(fā)生的應(yīng)力場較為混亂，沒有主流的方向，說明地殼應(yīng)力得到了較大釋放，使應(yīng)力場的方向的主流方向不明顯，意味著后面強(qiáng)震發(fā)生的可能性較小，而主震一個(gè)月后的余震叢中的確沒有大的地震發(fā)生。這些地震震源機(jī)制與應(yīng)力場的差別（最小空間旋轉(zhuǎn)角）明顯變小，意味著震源機(jī)制的發(fā)生與應(yīng)力場較為一致，呈現(xiàn)了地震發(fā)生的震源機(jī)制的主流方向，預(yù)示著將來會(huì)有較強(qiáng)地震發(fā)生，結(jié)果在主震發(fā)生后的12月29日，發(fā)生了本次余震序列的最大地震。在本次地震叢中，震源機(jī)制與應(yīng)力場的最小空間旋轉(zhuǎn)角再次增大，意味著后面沒有強(qiáng)余震發(fā)生。到目前為止，還未見到較大地震的發(fā)生，這說明余震的震源機(jī)制與局部應(yīng)力場的最小空間旋轉(zhuǎn)角確實(shí)跟強(qiáng)震的發(fā)生具有時(shí)空對應(yīng)關(guān)系。

仿照前人采用地震和主震震源機(jī)制的最小三維空間旋轉(zhuǎn)角來解釋強(qiáng)震發(fā)生的研究［7，10］，本研究也計(jì)算了本次地震序列的余震和主震震源機(jī)制的最小空間旋轉(zhuǎn)角隨時(shí)間的變化［圖7（a）］。我們發(fā)現(xiàn)，主震后48h內(nèi)的余震震源機(jī)制與主震震源機(jī)制的三維空間旋轉(zhuǎn)角較小，且后面的兩叢余震表現(xiàn)與應(yīng)力場最優(yōu)節(jié)面、余震震源機(jī)制的最小三維空間旋轉(zhuǎn)角是變化一致的，這一方面說明前人采用震源機(jī)制和主震震源機(jī)制的最小空間旋轉(zhuǎn)角代表應(yīng)力場是可行的;另一方面也表明，主震震源機(jī)制的表現(xiàn)不一定代表應(yīng)力場的應(yīng)力變化情況，而采用局部應(yīng)力場進(jìn)行最小空間旋轉(zhuǎn)角才是最合適的。

4結(jié)論與討論

本文從EMSC網(wǎng)站收集到2021年希臘克里特島地區(qū)MW6.0地震后多家機(jī)構(gòu)提供的多次余震資料，采用萬永革［11］提出的震源機(jī)制中心解程序求解了主震和余震序列的震源機(jī)制中心解，按照搜索1°的方向間隔和0.1的應(yīng)力比值間隔，通過網(wǎng)格搜索法對應(yīng)力場展開求解，計(jì)算了各震源機(jī)制中心解與該地區(qū)應(yīng)力場節(jié)面的最小空間旋轉(zhuǎn)角。通過對比研究結(jié)果，證實(shí)了地震震源機(jī)制與局部應(yīng)力場最優(yōu)節(jié)面震源機(jī)制的最小空間旋轉(zhuǎn)角與強(qiáng)震的發(fā)生有明顯的對應(yīng)關(guān)系，最小空間旋轉(zhuǎn)角的增大意味著其后的地震危險(xiǎn)性降低，而最小空間旋轉(zhuǎn)角的減小意味著后續(xù)該地區(qū)的地震危險(xiǎn)性增高。

需要指出的是，這里的最小空間旋轉(zhuǎn)角究竟多低才是后續(xù)強(qiáng)震發(fā)生的閾值，后續(xù)的強(qiáng)震的震級(jí)可能是多大，這些問題目前還沒有一致的標(biāo)準(zhǔn)。本文只是從短時(shí)間較為密集的地震震源機(jī)制數(shù)據(jù)的角度給出了一個(gè)實(shí)例。針對這些問題進(jìn)行進(jìn)一步研究也許是將來的一個(gè)方向，對地震預(yù)測研究有一定幫助。

地震釋放能量和地震矩關(guān)系的研究表明：走滑型地震多具有較高的視應(yīng)力水平，為傾滑型地震的數(shù)倍［2022］。對全球強(qiáng)震之間的觸發(fā)研究認(rèn)為具有較低視應(yīng)力水平的逆沖地震更容易被觸發(fā)。前人研究的很多地球物理現(xiàn)象與震源機(jī)制類型有關(guān)，如潮汐對地震的觸發(fā)［23］、前震發(fā)生率［24］、震級(jí)頻度關(guān)系［25］。前人討論了視應(yīng)力較高的走滑型強(qiáng)震發(fā)生前呈現(xiàn)明顯的應(yīng)力軸與主震震源機(jī)制的趨近［7］，對于傾滑型強(qiáng)震是否也可以觀測到應(yīng)力方向集中現(xiàn)象，目前尚不清楚。本文僅對正斷型地震給出了強(qiáng)震前震源機(jī)制與局部應(yīng)力場方向一致的實(shí)例，而中國大陸有很多正斷型的地塹帶，如藏南亞東—谷露裂谷帶、山西地塹帶等，這些地區(qū)發(fā)生正斷型的地震的可能性較大。本研究可為正斷型地震之前的應(yīng)力方向集中提供參考。將來也許從震源機(jī)制類型分類討論強(qiáng)震之前的震源機(jī)制的變化趨勢是地震預(yù)測研究的一個(gè)方向。

參考文獻(xiàn)（References）

［1］萬永革，沈正康，盛書中，等.2008年汶川大地震對周圍斷層的影響［J］.地震學(xué)報(bào)，2009，31（2）：128139.

WANYongge，SHENZhengkang，SHENGShuzhong，etal.Theinfluenceof2008Wenchuanearthquakeonsurroundingfaults［J］.ActaSeismologicaSinica，2009，31（2）：128139.

［2］解朝娣，王英楠，閆如玉，等.2014年四川康定MS6.3和MS5.8地震的應(yīng)力觸發(fā)研究［J］.地震研究，2021，44（1）：914.

XIEChaodi，WANGYingnan，YANRuyu，etal.Studyonstresstriggeringof2014KangdingMS6.3andMS5.8earthquakes［J］.JournalofSeismologicalResearch，2021，44（1）：914.

［3］郝平，傅征祥，田勤儉，等.昆侖山口西8.1級(jí)地震強(qiáng)余震庫侖破裂應(yīng)力觸發(fā)研究［J］.地震學(xué)報(bào)，2004，26（1）：3037.

HAOPing，F(xiàn)UZhengxiang，TIANQinjian，etal.Largeaftershockstriggeringbycoulombfailurestressfollowingthe2001MS8.1greatKunlunearthquake［J］.ActaSeismologicaSinica，2004，26（1）：3037.

［4］陳颙.用震源機(jī)制一致性作為描述地震活動(dòng)性的新參數(shù)［J］.地球物理學(xué)報(bào)，1978，21（2）：142159.

CHENRong.Consistencyoffocalmechanismasanewparameterindescribingseismicactivity［J］.ChineseJournalofSinica，1978，21（2）：142159.

［5］李欽祖，刁桂苓，戴英華.唐山地震序列的應(yīng)力釋放調(diào)整過程［J］.地球物理學(xué)報(bào)，1983，26（3）：224236.

LIQinzu，DIAOGuiling，DAIYinghua.ThestressreleaseandadjustmentprocessduringTangshanearthquakesequence［J］.ChineseJournalofGeophysics，1983，26（3）：224236.

［6］王俊國，刁桂苓.千島島弧大震前哈佛大學(xué)矩心矩張量（CMT）解一致性的預(yù)測意義［J］.地震學(xué)報(bào)，2005，27（2）：178183.

WANGJunguo，DIAOGuiling.ConsistentcmtsolutionsfromHarvarduniversitybeforethegreatearthquakesinKurileIslandsanditssignificanceforearthquakeprediction［J］.ActaSeismologicaSinica，2005，27（2）：178183.

［7］萬永革.美國Landers地震和HectorMine地震前震源機(jī)制與主震機(jī)制一致現(xiàn)象的研究［J］.中國地震，2008，24（3）：216225.

WANYongge.StudyonconsistencyoffocalmechanismofmainshockandthatofpreshocksinLandersandHectorMineearthquakeinUnitedStates［J］.EarthquakeResearchinChina，2008，24（3）：216225.

［8］KAGANYY.Simplifiedalgorithmsforcalculatingdoublecouplerotation［J］.GeophysicalJournalInternational，2007，171（1）：411418.

［9］HAUKSSONE.CrustalstructureandseismicitydistributionadjacenttothePacificandNorthAmericaplateboundaryinsouthernCalifornia［J］.JournalofGeophysicalResearch：SolidEarth，2000，105（B6）：1387513903.

［10］韓曉明，榮代潞.美國南加州地區(qū)1981—2011年MW≥6.0地震前發(fā)震應(yīng)力場與構(gòu)造應(yīng)力場趨于一致現(xiàn)象研究［J］.地震學(xué)報(bào)，2015，37（6）：948958，1.

HANXiaoming，RONGDailu.Consistencyofseismogenicstressfieldofpreshockstothetectonicstressfieldbeforeeightearthquakes（MW≥6.0）insouthernCaliforniaofUnitedStatesfrom1981to2011［J］.ActaSeismologicaSinica，2015，37（6）：948958，1

［11］萬永革.同一地震多個(gè)震源機(jī)制中心解的確定［J］.地球物理學(xué)報(bào)，2019，62（12）：47184728.

WANYongge.Determinationofcenterofseveralfocalmechanismsofthesameearthquake［J］.ChineseJournalofGeophysics，2019，62（12）：47184728.

［12］DANM.ActivetectonicsoftheMediterraneanregion［J］.GeophysicalJournalInternational，1972，30（2）：109185.

［13］LEPICHONX，ANGELIERJ.TheHellenicarcandtrenchsystem：akeytotheneotectonicevolutionoftheeasternMediterraneanarea［J］.Tectonophysics，1979，60（12）：142.

［14］萬永革.地震學(xué)導(dǎo)論［M］.北京：科學(xué)出版社，2016：395402.

WANYongge.Introductiontoseismology［M］.Beijing：SciencePress，2016：395402.

［15］萬永革.震源機(jī)制水平應(yīng)變花面應(yīng)變的地震震源機(jī)制分類方法及序列震源機(jī)制總體特征分析［JOL］.地球科學(xué)，2022：116（20220929）［20230402］.http：kns.cnki.netkcmsdetail42.1874.p.20220715.1532.014.html.

WANYongge.Focalmechanismclassificationbasedonarealstrainofthehorizontalstrainrosetteoffocalmechanismandcharacteristicanalysisofoverallfocalmechanismoftheearthquakesequence［JOL］.EarthScience，2022：116（20220929）［20230402］.http：kns.cnki.netkcmsdetail42.1874.p.20220715.1532.014.html.

［16］CAPUTOR，CATALANOS，MONACOC，etal.MiddlelateQuaternarygeodynamicsofCrete，southernAegean，andseismotectonicimplications［J］.BulletinoftheGeologicalSocietyofGreece，2017，43（1）：400.

［17］WANYG，SHENGSZ，HUANGJC，etal.ThegridsearchalgorithmoftectonicstresstensorbasedonfocalmechanismdataanditsapplicationintheboundaryzoneofChina，VietnamandLaos［J］.JournalofEarthScience，2016，27（5）：777785.

［18］萬永革.聯(lián)合采用定性和定量斷層資料的應(yīng)力張量反演方法及在烏魯木齊地區(qū)的應(yīng)用［J］.地球物理學(xué)報(bào)，2015，58（9）：31443156.

WANYongge.AgridsearchmethodfordeterminationoftectonicstresstensorusingqualitativeandquantitativedataofactivefaultsanditsapplicationtotheUrumqiarea［J］.ChineseJournalofGeophysics，2015，58（9）：31443156.

［19］萬永革，盛書中，許雅儒，等.不同應(yīng)力狀態(tài)和摩擦系數(shù)對綜合P波輻射花樣影響的模擬研究［J］.地球物理學(xué)報(bào)，2011，54（4）：9941001.

WANYongge，SHENGShuzhong，XUYaru，etal.EffectofstressratioandfrictioncoefficientoncompositePwaveradiationpatterns［J］.ChineseJournalofGeophysics，2011，54（4）：9941001.

［20］CHOYGL，BOATWRIGHTJL.Globalpatternsofradiatedseismicenergyandapparentstress［J］.JournalofGeophysicalResearch：SolidEarth，1995，100（B9）：1820518228.

［21］PREZCAMPOSX，BEROZAGC.Anapparentmechanismdependenceofradiatedseismicenergy［J］.JournalofGeophysicalResearch（SolidEarth），2001，106（B6）：1112711136.

［22］WUZL.Scalingofapparentstressfrombroadbandradiatedenergycatalogueandseismicmomentcatalogueanditsfocalmechanismdependence［J］.Earth，PlanetsandSpace，2001，53（10）：943948.

［23］TSURUOKAH，OHTAKEM，SATOH.Statisticaltestofthetidaltriggeringofearthquakes：contributionoftheoceantideloadingeffect［J］.GeophysicalJournalInternational，1995，122（1）：183194.

［24］REASENBERGPA.Foreshockoccurrenceratesbeforelargeearthquakesworldwide［J］.PureandAppliedGeophysics，1999，155（2）：355379.

［25］FROHLICHC，DAVISSD.Teleseismicbvalues;or，muchadoabout1.0［J］.JournalofGeophysicalResearch（SolidEarth），1993，98（B1）：631644.

（本文編輯：賈源源）

收稿日期：20220929

基金項(xiàng)目：國家自然科學(xué)基金（42174074，41674055）

第一作者簡介：李金磊（1999-），男，碩士研究生，主要從事地震學(xué)、勘探地球物理等方面的研究工作。Email：jllee0620@163.com。

通信作者：萬永革（1967-），男，研究員，主要從事構(gòu)造應(yīng)力場、地震應(yīng)力觸發(fā)等方面的研究工作。

Email：wanyg217217@vip.sina.com.cn。