京津冀地區(qū)地殼應(yīng)力場特征

王曉山 馮向東 趙英萍

摘要：基于標量斷層類型值，對京津冀地區(qū)及鄰區(qū)2 187個中小地震震源機制解進行分類，統(tǒng)計結(jié)果顯示研究區(qū)震源機制類型以走滑斷層和正斷層為主，P軸優(yōu)勢方位為NEE—EW和SWW—EW向;采用MSATSI軟件包反演該區(qū)1°×1°網(wǎng)格的精細地殼應(yīng)力場，結(jié)果表明：最大主壓應(yīng)力軸最優(yōu)解的優(yōu)勢方向為NEE—EW向，與P軸優(yōu)勢方位一致;所有網(wǎng)格的相對應(yīng)力大小R值均小于0.5，表明京津冀地區(qū)應(yīng)力狀態(tài)偏拉張性質(zhì)，而且最小主壓應(yīng)力軸的不確定度變化范圍相對穩(wěn)定，表明現(xiàn)今京津冀地區(qū)地殼應(yīng)力場處于一個相對統(tǒng)一的NNW—SSE向的拉張作用控制下。39°N以北地區(qū)最大主壓應(yīng)力軸方位最優(yōu)解顯示一定角度的偏轉(zhuǎn)，同時最大、中等、最小主壓應(yīng)力軸最優(yōu)解推斷的應(yīng)力狀態(tài)由西向東存在一個正斷層—走滑斷層—正斷層的轉(zhuǎn)換過程;而39°N以南地區(qū)的現(xiàn)今構(gòu)造應(yīng)力場保持穩(wěn)定，最優(yōu)主壓應(yīng)力軸呈NEE—SWW向，大部分網(wǎng)格應(yīng)力狀態(tài)顯示走滑型。構(gòu)造應(yīng)力場的反演結(jié)果與活動構(gòu)造、GPS主應(yīng)變方向和剪切波分裂的快波偏振方向等相關(guān)研究結(jié)果基本一致。

關(guān)鍵詞：京津冀地區(qū);震源機制;應(yīng)力張量阻尼反演;地殼應(yīng)力場

中圖分類號：P315.727文獻標識碼：A文章編號：1000-0666（2020）04-0610-10

0引言

地殼應(yīng)力場的時空分布特征，特別是地殼深部應(yīng)力狀態(tài)，是解決地震發(fā)生機理、區(qū)域地殼穩(wěn)定性、震源物理過程等地球動力學有關(guān)科學問題的基礎(chǔ)。地質(zhì)調(diào)查與大地測量結(jié)果表明，華北平原現(xiàn)今不存在大規(guī)模的拉張變形，區(qū)內(nèi)主要斷裂帶的滑移速率低（Shen et al，1996;郭良遷等，2008），是一個地殼變形相對較弱的地區(qū)。在這種緩慢的變形區(qū)域中，沉積蓋層或淺部巖石中的直接應(yīng)力測量通常不能代表深部的應(yīng)力場，而震源機制解是可用于估計緩慢變形區(qū)域中震源深度處應(yīng)力場的少數(shù)信息之一。

京津冀地區(qū)是我國東部強震活動頻度和強度最強的地區(qū)，同時也是新構(gòu)造活動最強烈的地區(qū)之一，區(qū)內(nèi)分布有張家口—渤海地震構(gòu)造帶（張渤地震帶）和華北平原地震構(gòu)造帶，中強地震絕大多數(shù)都發(fā)生在這2條構(gòu)造帶或其交會部位。該區(qū)是國內(nèi)使用地震記錄開展地殼應(yīng)力場較早的地區(qū)，李欽祖等（1973）利用河北紅山臺和沙城臺記錄的小地震P波初動資料確定了臺站周邊區(qū)域的應(yīng)力場;許忠淮等（1979，1983）對上述方法進行改進，形成了“綜合震源機制解法”，求解了京津唐張和華北地區(qū)的構(gòu)造應(yīng)力場方向。首都圈數(shù)字地震臺網(wǎng)于2002 年開始正式運行（劉瑞豐等，2008），為地震精確定位和震源機制解反演提供了有利條件，研究人員利用P波初動極性、振幅比和波形反演求解了大量的中小震震源機制解，且其質(zhì)量也更加可靠（胡新亮等，2004）?；诖?，不同研究人員使用不同方法開展了包括首都圈在內(nèi)的京津冀地區(qū)的地殼應(yīng)力場研究（蘭從欣等，2005;李瑞莎等，2008;張紅艷等，2009;武敏捷等，2011，2012;胡幸平，崔效鋒，2013;黃驥超，萬永革，2015;劉靜等，2016;劉麗等，2017;樊文杰等，2019）。其中大部分研究根據(jù)構(gòu)造特征和地震分布將研究區(qū)域分成幾個應(yīng)力小區(qū)，每個應(yīng)力小區(qū)的范圍都和許忠淮等（1983）對京津冀地區(qū)的分區(qū)基本一致，所不同的是張紅艷等（2009）在此基礎(chǔ)上進行了細分，但是所得的構(gòu)造應(yīng)力場方向基本相同;劉靜等（2016）和劉麗等（2017）對河北地區(qū)整體構(gòu)造應(yīng)力場的研究結(jié)果也驗證了該地區(qū)處于華北地區(qū)統(tǒng)一地殼應(yīng)力場的控制下（李欽祖，1980;魏光興等，1982;許忠淮等，1983）。

在京津冀地區(qū)震源機制解數(shù)據(jù)日益豐富的當下，武敏捷等（2011）采用細分網(wǎng)格加滑動步長的方式初步討論了華北北部地區(qū)的構(gòu)造應(yīng)力場，由于每個網(wǎng)格都是獨立反演，分區(qū)方式不同和分區(qū)大小的變化都可能對反演的應(yīng)力張量結(jié)果產(chǎn)生影響（Maury et al，2013）。相鄰2個網(wǎng)格主應(yīng)力方向的不同或許反映了應(yīng)力場真實的變化，但也可能是震源機制數(shù)據(jù)的誤差或反演的約束較差造成的假象。地殼中的真實應(yīng)力場是連續(xù)分布的，為了解決應(yīng)力反演模型依賴于分區(qū)的問題，更準確地反映研究區(qū)域的應(yīng)力場空間變化特征，Hardebeck和Michael（2006）構(gòu)建了一組可調(diào)整的阻尼參數(shù)，引入平滑約束來抑制相鄰網(wǎng)格應(yīng)力模型的差異，提出了區(qū)域尺度的應(yīng)力張量阻尼反演方法。本文基于京津冀地區(qū)及鄰區(qū)的中小地震震源機制解，將該方法應(yīng)用于對該區(qū)精細地殼應(yīng)力場特征的研究。

1數(shù)據(jù)和方法

1.1震源機制解數(shù)據(jù)

本文所用震源機制解共2 187個，震級區(qū)間2.0≤ML≤5.5，覆蓋范圍包括京津冀地區(qū)以及山西、遼寧、內(nèi)蒙部分地區(qū)。數(shù)據(jù)來源分為2個部分：第一部分為2002—2008年京津冀地區(qū)ML≥2.0地震的震源機制解，使用梁尚鴻等（1984）提出的利用區(qū)域地震臺網(wǎng)垂直向P，S最大振幅比資料測定小震震源參數(shù)的方法，該方法反演結(jié)果的可靠性在前人的對比論證中得到了證實（胡新亮等，2004;付虹等，2009），且已經(jīng)在首都圈的應(yīng)力場研究中得到應(yīng)用（李瑞莎等，2008;張紅艷等，2009;武敏捷等，2011，2012）;第二部分為2009—2019年10月京津冀地區(qū)ML≥2.0地震的震源機制解，使用格點嘗試法（許忠淮等，1983）利用P波初動極性計算震源機制解，P波初動極性資料主要來自“地震編目系統(tǒng)”下載的震相報告和從事件波形中人工讀取的P波初動。為了保證震源機制解的可靠性，參與計算的P波初動極性數(shù)據(jù)一般要求10個以上，且在震中周圍均勻分布，矛盾比小于0.20。胡新亮等（2004）對首都圈數(shù)字地震臺網(wǎng)內(nèi)的地震分別使用格點嘗試法和垂直向最大振幅比方法測定震源機制，二者結(jié)果基本一致，表明2種方法計算結(jié)果可以聯(lián)合起來進行京津冀地區(qū)地殼應(yīng)力場特征分析。

1.2震源機制解分類方法

對震源機制解分類有多種方法，如使用P，T軸參數(shù)進行三角形分類（Frohlich，1992，2001）、根據(jù)震源機制解3個應(yīng)力軸傾角的世界應(yīng)力圖劃分標準（Zoback，1992）、以3個應(yīng)力軸的傾角最大值判斷震源機制破裂類型（刁桂苓等，2011）等，以上都是以應(yīng)力軸參數(shù)對震源機制解分類。本文采用Shearer 等（2006）提出的依據(jù)節(jié)面滑動角ri對震源機制解進行分類的方法。該方法的優(yōu)勢在于提供單個標量斷層類型值用于表征震源錯動類型，標量斷層類型值fptype為：fptype=r90（1）式中：當r190ri，ri≤90（i=1，2）（2）標量斷層類型值從-1（正斷層）→0（走滑斷層）→1（逆斷層）變化。Lin 和Okubo（2016）取一個介于0～1的Xf值來區(qū)分斷層類型：-1≤fptype<-Xf，正斷層-Xf≤fptype≤Xf，走滑斷層Xf

1.3應(yīng)力場反演方法

目前，使用震源機制解進行應(yīng)力場反演應(yīng)用最廣泛的2種方法為：震源機制應(yīng)力反演（FMSI）的網(wǎng)格搜索方法（Gephart，F(xiàn)orsyth，1984;Gephart，1990）和自助線性應(yīng)力反演（LSIB）方法（Michael，1984，1987）。Hardebeck和Hauksson（2001）使用含噪聲的合成數(shù)據(jù)對這2種方法進行測試，2種方法都能精確地確定應(yīng)力方向;FMSI方法通常對于高質(zhì)量數(shù)據(jù)更準確，而LSIB方法對于離散度大的數(shù)據(jù)更準確;LSIB方法產(chǎn)生的置信區(qū)間通常較為合理，而FMSI方法的置信區(qū)間通常過大。合成數(shù)據(jù)測試結(jié)果表明，只要滿足應(yīng)力場均勻和斷層的滑動方向與剪應(yīng)力方向一致的模型假設(shè)，由LSIB方法產(chǎn)生的置信區(qū)間也應(yīng)適用于真實數(shù)據(jù)。Hardebeck和Michael（2006）引入阻尼最小二乘反演方法對LSIB方法加以改進，提出了區(qū)域尺度的應(yīng)力張量阻尼反演方法SATSI（Spatial And Temporal Stress Inversion）算法;Martínez-Garzón等（2014）將SATSI算法移植到Matlab環(huán)境中，形成一款結(jié)合可靠經(jīng)典方法的、新簡化用戶處理及可視化工具的應(yīng)力反演Matlab軟件包。由于本文所用震源機制解數(shù)據(jù)量較大，且2種不同的方法計算得到的數(shù)據(jù)離散度較大，故選用MSATSI軟件包反演京津冀地區(qū)地殼應(yīng)力場，以期可以準確地約束應(yīng)力張量的方向（Hardebeck，Michael，2006）。

2震源機制解統(tǒng)計特征

本文參考Lin和Okubo（2016）的做法，對京津冀地區(qū)中小地震震源機制解的破裂類型比例進行定量估計，取Xf=0.33，將[-1，1]分成3個相等的部分，震源機制解也相應(yīng)分為走滑斷層、正斷層和逆斷層（圖1）。對上述地震進行統(tǒng)計后發(fā)現(xiàn)，京津冀地區(qū)有1 315個走滑型震源機制解，占整體數(shù)量的60.1%，628個正斷型地震，占28.7%，244個逆斷型地震，占11.2%，該區(qū)震源錯動類型主要以走滑型和正斷型地震為主。不同類型的震源機制解空間分布（圖1）顯示，走滑型地震在整個區(qū)域彌散分布，而正斷型和逆斷型地震基本上沿著張渤地震帶和華北平原地震帶分布，沒有形成明顯的單一類型的叢集區(qū)。

將震源機制解按節(jié)面走向、傾角、滑動角和P軸方位角繪制成直方圖進行統(tǒng)計（圖2）。由于絕大多數(shù)地震震級太小，無法區(qū)分真正的斷層面和輔助斷層面，故將2 個節(jié)面等同看待，參數(shù)合并統(tǒng)計。節(jié)面走向大體呈各個方位均勻分布，在NNE—SSW和NWW—SEE向稍顯優(yōu)勢（圖2a）;節(jié)面傾角以接近直立的高傾角最多，表明以走向滑動為主，幾乎沒有近水平的節(jié)面（圖2b）;節(jié)面滑動角以走向滑動為主，左旋、右旋走滑數(shù)量相仿，正傾滑明顯多于逆傾滑（圖2c）。P軸方位角在NEE—EW和SWW—EW向存在顯著優(yōu)勢取向，可以看作是京津冀地區(qū)最大水平壓應(yīng)力方向（圖2d）。上述震源機制解的統(tǒng)計特征與前人研究結(jié)果（李欽祖，1980;許忠淮等，1983;武敏捷等，2012;劉靜等，2016;劉麗等，2017）基本一致。

3地殼應(yīng)力場反演

采用1°×1°的步長，將京津冀地區(qū)（36.0°～43.0°N，113°～120°E）劃分網(wǎng)格，使用MSATSI軟件包反演該地區(qū)精細的地殼應(yīng)力場。該方法通過在反演過程中加上一個合理的阻尼，使反演的相鄰區(qū)域的應(yīng)力張量變化最小，消除人為劃分反演區(qū)域?qū)е碌膽?yīng)力偏轉(zhuǎn)假象，更好地凸顯整體應(yīng)力場特征（Hardebeck，Michael，2006）。在反演之前要先確定阻尼系數(shù)e，它控制著理論值與觀測數(shù)據(jù)之間擬合殘差和應(yīng)力反演模型長度的相對權(quán)重，對反演結(jié)果起著至關(guān)重要的作用。圖3為模型長度—數(shù)據(jù)擬合殘差的折中曲線，其拐點在1.2，因此最佳阻尼系數(shù)取e=1.2。反演中每個網(wǎng)格采用的震源機制解至少為6個，隨機選取其中一個節(jié)面進行構(gòu)造應(yīng)力場反演。反演應(yīng)力場參數(shù)的置信水平的可設(shè)定范圍為68%～95%，本文設(shè)定為95%。重采樣的迭代次數(shù)對于解的不確定度評價很重要，重采樣次數(shù)過少導致評估結(jié)果沒有統(tǒng)計意義，過多會導致計算量大，反演中重采樣次數(shù)的默認范圍是1 000～5 000次（Martínez-Garzón et al，2014），本文對所有數(shù)據(jù)進行2 000次抽樣估計。經(jīng)過計算得到最優(yōu)狀態(tài)下的3個應(yīng)力主軸的方向和相對應(yīng)力大?。≧值）（Gephart，F(xiàn)orsyth，1984）為：R=（σ2-σ1）/（σ3-σ1）（4）式中：σ1，σ2，σ3分別表示最大、中間、最小主壓應(yīng)力。

從表1和圖4可以看出，最優(yōu)主壓應(yīng)力軸σ1的優(yōu)勢方向為NEE—EW向，且有大約半數(shù)的網(wǎng)格傾角接近直立，表明應(yīng)力狀態(tài)為正斷層類型，而網(wǎng)格編號05，09，23的最優(yōu)主壓應(yīng)力軸σ1方向為NNW—NWW向，傾角都接近直立，表明最大主應(yīng)力軸垂直，這與前人用中強震震源機制解和綜合震源機制解研究華北地區(qū)地殼應(yīng)力場的結(jié)果有所差異（李欽祖，1980，1982;魏光興等，1982;許忠淮等，1983），主要表現(xiàn)為最優(yōu)主應(yīng)力軸的傾角直立或水平，這可能與本文所用的震源機制解多為小震、易受局部構(gòu)造運動控制、且具有隨機性特征有關(guān)，反映了整體背景應(yīng)力場下的非均勻性和復雜性;最優(yōu)主張應(yīng)力軸σ3的優(yōu)勢方向為SSE—NWW向，且傾角均接近水平。

從圖4中最大主壓應(yīng)力σ1軸的指示方向可以看出，在39°N以北的首都圈地區(qū)最大主壓應(yīng)力軸自西向東存在由NE向近EW向的偏轉(zhuǎn)，這與前人對首都圈構(gòu)造應(yīng)力場的研究結(jié)果基本一致（張紅艷等，2009;武敏捷等，2012;黃驥超，萬永革，2015;樊文杰等，2019），而且本文的分區(qū)更加精細，應(yīng)力場的變化連續(xù)可靠。由最大、中間、最小主壓應(yīng)力軸最優(yōu)解推斷的張渤地震帶陸地段的應(yīng)力狀態(tài)表現(xiàn)為，自西向東存在一個正斷層—走滑斷層—正斷層的轉(zhuǎn)換過程，這與張紅艷等（2009）、武敏捷等（2012）、樊文杰等（2019）以走滑型為主的應(yīng)力場特征存在一定差異;同時，在黃驥超和萬永革（2015）對首都圈的構(gòu)造應(yīng)力場反演結(jié)果中，北京—張家口—大同地區(qū)和唐山及鄰區(qū)都表現(xiàn)出少量正斷層性質(zhì)。39°N以南地區(qū)現(xiàn)今構(gòu)造應(yīng)力場保持穩(wěn)定，最優(yōu)主壓應(yīng)力軸σ1為NEE—SWW向，大部分網(wǎng)格應(yīng)力狀態(tài)顯示走滑型。

相對應(yīng)力大小R值是中間主應(yīng)力σ2軸的值，是更接近于最大主壓應(yīng)力軸σ1還是最小主壓應(yīng)力軸σ3的大小的度量。當R值趨于中間值0.5時，表示3個應(yīng)力主軸方位可完全分辨;而當R值趨于取值端點0或者1時，表示壓縮或拉張應(yīng)力狀態(tài)完全不能分辨，只能確定它在垂直于拉張或壓縮應(yīng)力主軸的平面內(nèi)。一般認為當R值小于0.5時，獲得的應(yīng)力狀態(tài)偏拉張性質(zhì)，反之則為偏壓縮性質(zhì)（Guiraud et al，1989）。從表1和圖5可以看出，所有網(wǎng)格的R值均小于0.5，表明京津冀地區(qū)的應(yīng)力狀態(tài)呈現(xiàn)偏拉張性質(zhì)。同時從圖4可知，最大主應(yīng)力軸σ1的不確定范圍最大，最大主壓應(yīng)力軸σ1和中間應(yīng)力軸σ2表現(xiàn)的壓應(yīng)力狀態(tài)是一致的，兩軸無法區(qū)分（Guiraud et al，1989;萬永革等，2011）;而最小主壓應(yīng)力軸σ3的不確定度變化范圍相對穩(wěn)定，我們可以推斷京津冀地區(qū)現(xiàn)今地殼應(yīng)力場處于一個相對統(tǒng)一的NNW—SSE向的拉張作用控制下。4結(jié)論和討論

基于京津冀地區(qū)及鄰區(qū)2 187個中小地震震源機制解，通過計算標量斷層類型值將其劃分為3種震源錯動類型，采用MSATSI軟件包反演該地區(qū)1°×1°網(wǎng)格的精細地殼應(yīng)力場，研究結(jié)果表明：研究區(qū)震源機制類型以走滑斷層和正斷層為主，占比約90%;不同類型的震源機制解沒有形成明顯的單一類型的叢集區(qū)，反映出小地震的發(fā)生除了受背景應(yīng)力場的控制，還受局部構(gòu)造條件控制且具有一定的隨機性。這為我們利用小地震的震源機制解來推斷構(gòu)造應(yīng)力主軸的方向提供了更大的優(yōu)越性，震源機制解的一對節(jié)面解中既含有構(gòu)造應(yīng)力場的信息，也可能含有局部介質(zhì)不均勻結(jié)構(gòu)的影響（許忠淮等，1983）。

地殼應(yīng)力場反演結(jié)果顯示，所有網(wǎng)格的相對應(yīng)力大小R值均小于0.5，表明京津冀地區(qū)應(yīng)力狀態(tài)偏拉張性質(zhì)，而且最小主壓應(yīng)力軸σ3的不確定度變化范圍相對穩(wěn)定，表明現(xiàn)今京津冀地區(qū)地殼應(yīng)力場處于一個相對統(tǒng)一的NNW—SSE向的拉張作用控制下。這一結(jié)果與前人的研究結(jié)果比較一致，高名修（1979）認為華北地區(qū)處于引張應(yīng)力場環(huán)境，有限元模擬顯示主張應(yīng)力大約為主壓應(yīng)力2～8倍（劉峽等，2010）。利用反演結(jié)果中最大、中間、最小主壓應(yīng)力軸最優(yōu)解推斷的應(yīng)力狀態(tài)顯示，39°N以北地區(qū)的最大主壓應(yīng)力軸方向存在一定角度的偏轉(zhuǎn)，同時應(yīng)力狀態(tài)由西向東存在一個正斷層—走滑斷層—正斷層的轉(zhuǎn)換過程;而39°N以南地區(qū)的現(xiàn)今構(gòu)造應(yīng)力場保持穩(wěn)定，大部分網(wǎng)格應(yīng)力狀態(tài)顯示走滑型。

39°N以北地區(qū)與前人研究中的張渤帶陸地段和首都圈范圍基本一致，大致分為3段：西段、中段和東段。西段為張家口—大同地區(qū)，應(yīng)力狀態(tài)表現(xiàn)為正斷層類型，該區(qū)為山西斷陷帶右旋走滑在NE端形成的一個張性構(gòu)造區(qū)（徐錫偉等，2002），活動構(gòu)造、遙感和垂直形變研究結(jié)果都表明山西斷陷帶中北段至張家口地區(qū)總體以張性正斷層活動為主（孫啟凱等，2018;吳玉濤等，2018;申星等，2019;彭遠黔等，2019），該區(qū)的中小地震震源機制解以正斷型和走滑型為主（孫小入等，2019）。中段為北京—天津地區(qū)，應(yīng)力狀態(tài)為走滑類型和正斷層，該區(qū)歷史上曾發(fā)生過1679年三河—平谷地震，使用現(xiàn)今小震資料推斷該地震的錯動類型為走滑兼少量正斷（王曉山，2017），2018年發(fā)生在該區(qū)的河北永清地震的震源機制為走滑兼正斷類型（王曉山等，2018;李赫等，2020），GPS觀測資料亦顯示該區(qū)存在高剪切變形區(qū)（陳長云，2016）。東段為唐山地區(qū)，應(yīng)力狀態(tài)為走滑類型和正斷層，楊雅瓊等（2016）對唐山地震序列的現(xiàn)今小震震源機制分段擬合最佳應(yīng)力張量，得到的結(jié)果為正斷層和走滑類型兼有。綜合前人多種資料的論證表明，本文的研究結(jié)果與前人結(jié)果具有較好的一致性，但仍存在少許差異，可能由所使用的反演方法、網(wǎng)格劃分和數(shù)據(jù)量等導致，有待進一步研究分析。

GPS主應(yīng)變率結(jié)果顯示：京津冀大部分地區(qū)地殼變形特征主要受NEE的主壓應(yīng)變控制（陳長云，2016），京津冀地區(qū)的剪切波分裂的快波偏振方向為NEE—EW向（吳晶等，2008a，b;趙博等，2011;吳鵬等，2017）。該結(jié)果與該區(qū)使用中小地震震源機制解反演的構(gòu)造應(yīng)力場最優(yōu)解基本一致，同時也和前人使用中強地震震源機制解和綜合震源機制解得到的華北地區(qū)一致性良好的應(yīng)力場吻合，表明中小地震震源機制的結(jié)果也可以較好地用于區(qū)域構(gòu)造應(yīng)力場的研究。京津冀地區(qū)是中國大陸較早實現(xiàn)數(shù)字化地震觀測網(wǎng)絡(luò)的地區(qū)，目前積累了豐富的中小地震波形資料，開展中小地震震源機制解及構(gòu)造應(yīng)力場的動態(tài)演化研究可以為該區(qū)地震動力學研究提供具有物理意義的信息，更好地服務(wù)于該地區(qū)的防震減災(zāi)綜合研究。

參考文獻：

陳長云.2016.張家口—渤海斷裂帶分段運動變形特征分析[J].地震，36（1）：1-11.

刁桂苓，徐錫偉，陳于高，等.2011.汶川MW7.9和集集MW7.6地震前應(yīng)力場轉(zhuǎn)換現(xiàn)象及其可能的前兆意義[J].地球物理學報，54（1）：128-136，doi：10.3969/j.issn.0001-5733.2011.01.014.

樊文杰，崔效鋒，胡幸平，等.2019.首都圈地區(qū)震源機制解及現(xiàn)今構(gòu)造應(yīng)力場時空變化特征研究[J].地震學報，41（1）：33-45，doi：10.11939/jass.20180064.

付虹，王紹晉，龍曉帆，等.2009.利用數(shù)字化速度波形和仿真位移波形資料求解震源機制解的比較[J].地震研究，32（1）：18-24.

高名修.1979.華北斷塊構(gòu)造區(qū)的現(xiàn)代引張應(yīng)力場[J].地震地質(zhì)，1（2）：1-12.

郭良遷，馬青，杜雪松，等.2008.華北地區(qū)斷層形變與地震的關(guān)系[J].大地測量與地球動力學，28（3）：14-20.

胡新亮，刁桂苓，馬瑾，等.2004.利用數(shù)字地震記錄的P，S振幅比資料測定小震震源機制解的可靠性分析[J].地震地質(zhì)，26（2）：347-354.

Gephart J W，F(xiàn)orsyth D W.1984.An improved method for determining the regional stress tensor using earthquake focal mechanism data：application to the San Fernando Earthquake sequence[J].J Geophys Res，89（B11）：9305-9320.

Gephart J W.1990.FMSI：A FORTRAN program for inverting fault/slikenside and earthquake focal mechanism data to obtain the regional stress tensor[J].Computers and Geosciences，16（7）：953-989.

Guiraud M，Laborde O，Philip H.1989.Characterization of various types of deformation and their corresponding deviatoric stress tensors using microfault analysis[J].Tectonophysics，170（3-4）：289-316.

Hardebeck J L，Hauksson E.2001.Stress orientations obtained from earthquake focal mechanisms：what are appropriate uncertainty estimates？[J].Bull Seis Soc Amer，91（2）：250-262.

Hardebeck L L，Michael A J.2006.Damped regional-scale stress inversions：methodology and examples for southern california and the coalinga aftershock sequence[J].J Geophys Res，111（B11）：B11310，doi：10.1029/2005JB004144.

Lin G，Okubo P G.2016.A large refined catalog of earthquake relocations and focal mechanisms for the island of Hawaii and Its seismotectonic implications[J].J Geophys Res，121（7）：5031-5048，doi：10.1002/2016JB013042.

Martínez-Garzón P，Kwiatek G，Ickrath M，et al.2014.MSATSI：A MATLAB package for stress inversion combining solid classic methodology，a new simplified user-handling，and a visualization tool[J].Seismol Res Lett，85（4）：896-904.

Maury J，Cornet F H，Dorbath L.2013.A review of methods for determining stress fields from earthquakes focal mechanisms;application to the sierentz 1980 seismic crisis（upper Rhine graben）[J].Bull Soc geol，184（4-5）：319-334.

Michael A J.1984.Determination of stress from slip data：faults and Folds[J].J Geophys Res，89（B13）：11517-11526.

Michael A J.1987.Use of focal mechanisms to determine stress：a control study[J].J Geophys Res，92（B1）：357-368.

Shearer P M，Prieto G A，Hauksson E.2006.Comprehensive analysis of earthquake source spectra in Southern California[J].J Geophys Res，111（B6）：303，doi：10.1029/2005JB003979.

Shen Z K，Jackson D D，Ge B X.1996.Crustal deformation across and beyond the Los Angeles Basin from geodetic measurements[J].J Geophys Res，101（B12）：27957-27980.

Zoback M L.1992.First- and secondorder patterns of stress in the lithosphere：the world stress map project[J].J Geophys Res，97（B8）：11703-11728，doi：10.1029/92JB00132.

Characteristics of Crustal Stress Field in Beijing-Tianjin-Hebei Region