楊雅瓊 王曉山 萬永革 盛書中 陳 婷

1) 中國石家莊050021河北省地震局 2) 中國北京100081中國地震局地球物理研究所3) 中國河北三河065201防災(zāi)科技學(xué)院

?

由震源機(jī)制解推斷唐山地震序列發(fā)震斷層的分段特征*

楊雅瓊1)王曉山1,2),*萬永革3,2)盛書中3)陳 婷1)

1) 中國石家莊050021河北省地震局 2) 中國北京100081中國地震局地球物理研究所3) 中國河北三河065201防災(zāi)科技學(xué)院

使用唐山地區(qū)2002年1月—2015年11月ML≥2.5地震的255個震源機(jī)制解，采用構(gòu)造應(yīng)力場均勻性的分段方法對唐山地震序列的發(fā)震斷層進(jìn)行分段．在已有唐山地區(qū)震源機(jī)制解分區(qū)特征的基礎(chǔ)上，給出了5個參考應(yīng)力張量，并通過差異顯著性的z值檢驗(yàn)計算將唐山地震序列的發(fā)震斷層分為寧河、唐山、灤縣和盧龍等4個子段，進(jìn)而分別對4個子段的應(yīng)力場進(jìn)行反演．結(jié)果顯示：4個子段的最大主壓應(yīng)力方向均呈近EW向，且唐山、寧河和盧龍子段的應(yīng)力場均表現(xiàn)出較大的拉張分量；唐山、寧河子段的最佳應(yīng)力張量與唐山主震對唐山斷裂帶兩端點(diǎn)所產(chǎn)生的引張應(yīng)力場的作用方式一致．此外，唐山子段的應(yīng)力場符合基于接收函數(shù)給出的上地幔物質(zhì)隆升模型，灤縣子段走滑型的應(yīng)力狀態(tài)反映了該區(qū)的共軛構(gòu)造運(yùn)動，盧龍子段的最佳應(yīng)力張量為正斷兼右旋走滑．從當(dāng)前唐山地震序列發(fā)震斷層分段的應(yīng)力場特征可以推斷，現(xiàn)今唐山地區(qū)的地震活動具有繼承性，主要受區(qū)域構(gòu)造應(yīng)力場和該區(qū)深、淺共存的斷裂構(gòu)造體系控制．

唐山地震序列 震源機(jī)制解 構(gòu)造應(yīng)力場 發(fā)震斷層分段

引言

1976年7月28日河北唐山MS7.8地震發(fā)生在唐山菱形塊體內(nèi)部NE走向的唐山斷裂帶上(虢順民等，1977)．該斷裂帶可劃分為南、北兩個段落，北段自西向東依次分布有陡河斷裂、唐山—巍山—長山南坡斷裂和唐山—古冶斷裂，南段自西向東分布著唐山—豐南斷裂和唐山—古冶斷裂(尤惠川等，2002)；與NW向的薊運(yùn)河斷裂、灤縣—樂亭斷裂和EW向的豐臺—豐南斷裂交匯，形成南北兩段構(gòu)造斷裂閉鎖和唐山、寧河、灤縣等3個應(yīng)力集中區(qū)(馬瑾等，1980)．?dāng)嗔哑屏训挠邢拊M結(jié)果(Xie，Yao，1991；Huang，Yeh，1997)顯示，唐山地震發(fā)震斷層的破裂起始于唐山斷裂帶南段的北端(唐山斷裂帶南北兩段的交接點(diǎn)在唐山市附近)，然后沿走向?yàn)镹E30°的南段擴(kuò)展，5.9 s之后破裂開始沿走向約為NE50°—60°的北段傳播(Náběleketal，1987)．郭慧等(2011a，b)根據(jù)唐山地震的最新地表破裂帶的位移特征，將唐山地震地表破裂帶分為南北兩段：北段以右旋走滑為主，兼有西升東降的垂直位移；南段以東升西降的傾滑活動為主，兼有右旋走滑活動；且兩段地表破裂帶的重疊區(qū)與MS7.8地震的微觀震中相近．Shedlock等(1987)關(guān)于唐山地震序列早期余震分布的研究也表明南北兩段的交接點(diǎn)在唐山市附近．唐山主震發(fā)生后3.5小時和15小時左右，分別發(fā)生了寧河MS6.2和灤縣MS7.1地震；主震后大約3.5個月，在主震破裂的SW延伸斷裂部位上發(fā)生寧河MS6.9地震(國家地震局《一九七六年唐山地震》編輯組，1982)．顯然，這3次地震是因唐山主震產(chǎn)生破裂而發(fā)生的觸發(fā)型破裂事件，與主震共同構(gòu)成了一次復(fù)雜的板內(nèi)地震序列(Butleretal，1979；劉桂萍，傅征祥，2000；Robinson，Zhou，2005；萬永革等，2008)．

對于唐山地震序列的應(yīng)力場研究，由于受震源機(jī)制解豐富程度和求解精度所限(許忠淮等，1994)，早期僅能在時間進(jìn)程和空間分布上對震源區(qū)的應(yīng)力場進(jìn)行粗略的分析．李欽祖等(1983)關(guān)于唐山地震序列早期震源機(jī)制解的研究揭示了震源區(qū)應(yīng)力場由主震前的緊張平衡態(tài)經(jīng)主震后的極不穩(wěn)定態(tài)逐漸向松弛平衡態(tài)演變；刁桂苓等(1995)關(guān)于唐山地震序列震源機(jī)制解的系統(tǒng)聚類結(jié)果顯示，該序列的震源機(jī)制解類型多，且取向分散；Xu和Wang(1986)利用滑動方向擬合法得到的唐山老震區(qū)平均應(yīng)力場結(jié)果顯示，震后唐山地區(qū)的主壓應(yīng)力軸發(fā)生了順時針旋轉(zhuǎn)．2002年投入運(yùn)行的首都圈數(shù)字地震臺網(wǎng)極大地改善了唐山地區(qū)地震震源的定位精度和震源機(jī)制的求解質(zhì)量，為求取唐山地區(qū)更加精細(xì)的應(yīng)力場空間分布特征奠定了數(shù)據(jù)基礎(chǔ)．李軼群和王健(2008)采用網(wǎng)格密集值方法將唐山余震區(qū)分為5個區(qū)域并分別求解各分區(qū)的小震綜合斷層面解；張宏志等(2008)利用雙差定位結(jié)果對唐山震區(qū)分區(qū)，并給出各分區(qū)震源機(jī)制解的平均應(yīng)力軸取向；萬永革等(2008)則利用模擬退火算法和高斯-牛頓算法反演了唐山地震序列的5段地震斷層面的走向、傾角和滑動角．上述研究均從地震活動的叢集性出發(fā)，對唐山地震序列的發(fā)震斷層進(jìn)行分段，然后分別求取各子段的應(yīng)力場，但均并未將構(gòu)造應(yīng)力場均勻性對地震活動的影響考慮在內(nèi)．

唐山地區(qū)的斷裂構(gòu)造和主震的破裂過程極為復(fù)雜，各子段在區(qū)域應(yīng)力場作用下的地震活動存在明顯差異，本文將通過基于構(gòu)造應(yīng)力場均勻性對斷層或板塊邊界進(jìn)行分段(Wyss，Lu，1995)的方法，利用唐山地震序列近期的255個震源機(jī)制解，從構(gòu)造應(yīng)力場角度對唐山地震序列的發(fā)震斷層進(jìn)行分段，并給出經(jīng)過差異顯著性檢驗(yàn)的各子段的應(yīng)力場特征，以期為深入理解唐山地震序列的發(fā)震機(jī)理、發(fā)震構(gòu)造以及該地區(qū)的地震危險性評估提供詳盡的基礎(chǔ)資料．

1 構(gòu)造背景與數(shù)據(jù)

唐山地震發(fā)生在近EW向的陰山—燕山南緣活動構(gòu)造帶與NNE向的華北平原活動構(gòu)造帶的交匯處．如圖1所示，唐山地區(qū)主要發(fā)育有3組淺層活動斷裂，即NE向的豐臺—野雞坨斷裂(F1)、寧河—昌黎斷裂(F3)、唐山斷裂帶(F5)，NW向的灤縣—樂亭斷裂(F2)、薊運(yùn)河斷裂(F4)和EW向的豐臺—豐南斷裂(F9)，唐山地震即發(fā)生在唐山斷裂帶的深部．

本文基于首都圈數(shù)字地震臺網(wǎng)記錄的2002年1月—2015年11月唐山地區(qū)ML≥2.5地震的波形資料，使用梁尚鴻等(1984)提出的垂直向直達(dá)P波、S波最大振幅比方法計算得到255個震源機(jī)制解，如圖1中灰色和紅色沙灘球所示，其震源位置利用hypo2000(Klein，1978)經(jīng)重新絕對定位得到．圖1中壓縮區(qū)用黑色填充的震源機(jī)制解為矩心矩張量解(Dziewonskietal，1981；Ekstr?metal，2012)，震中位置引自《中國地震目錄(公元1970—1979年)》(顧功敘等，1983). 由圖1可以看出，震源機(jī)制解相對集中在唐山斷裂帶的北段和灤縣地區(qū)，在盧龍和寧河地區(qū)呈叢集性分布，基本上與馬瑾等(1980)給出的震前3個應(yīng)力集中區(qū)分布一致．

圖1 唐山地區(qū)地質(zhì)構(gòu)造和2002年1月—2015年11月ML≥2.5地震的震源機(jī)制解分布F1: 豐臺—野雞坨斷裂；F2: 灤縣—樂亭斷裂；F3: 寧河—昌黎斷裂；F4: 薊運(yùn)河斷裂；F5: 唐山斷裂帶；F6: 盧龍斷裂；F7: 漢沽斷裂；F8: 滄東斷裂；F9: 豐臺—豐南斷裂

根據(jù)震源機(jī)制解的3個應(yīng)力軸傾角大小，Zoback(1992)將震源機(jī)制解分為6種類型，即正斷型、正走滑型、走滑型、逆走滑型、逆斷型和無法確定型．對本文所獲取的255個震源機(jī)制解進(jìn)行分類統(tǒng)計，結(jié)果為：走滑型地震為93次，占總數(shù)的36.5%；正斷型和正斷兼走滑分量類型為91次，占35.9%；逆斷型和逆斷兼走滑分量類型為23次，占9.0%；無法確定型為48次，占18.4%．可以看出，走滑型與正斷型所占的比例基本相同，這可能是由于本文將具有走滑分量的正斷層和逆斷層分別劃歸至正斷型和逆斷型所致，從整體上看研究區(qū)的應(yīng)力狀態(tài)以剪切拉張為主．

2 方法

Wyss和Lu(1995)首次提出了利用構(gòu)造應(yīng)力場均勻性對斷層或板塊邊界進(jìn)行分段的方法，并將該方法應(yīng)用于圣安德列斯斷層的分段研究，結(jié)果顯示，該斷層的5個子斷層中，4個子斷層的位置與基于非定量的構(gòu)造劃分法得到的分段結(jié)果是一致的．因此，應(yīng)用該方法對斷層進(jìn)行分段，不但使斷層或板塊邊界的劃分更量化和客觀化，而且還能找出構(gòu)造劃分法難以找出的子段．之后，該方法被應(yīng)用于阿留申板塊邊界(Lu，Wyss，1996)和阿拉斯加俯沖帶(Luetal，1997)的分段；這些地區(qū)的構(gòu)造應(yīng)力場相對簡單，因此應(yīng)用效果很好．Ratchkovski(2003)將其用于板緣地區(qū)的阿拉斯加德納里(Denali)斷層的分段，盛書中和萬永革(2012)將其用于板內(nèi)地區(qū)的汶川地震斷層的分段研究，二者均給出了較好的分段結(jié)果．這些研究結(jié)果表明利用構(gòu)造應(yīng)力場的均勻性對板塊邊界或斷層進(jìn)行分段的方法適用于各種類型的構(gòu)造環(huán)境，能夠?qū)Π鍓K邊界或斷層給出經(jīng)過統(tǒng)計顯著性檢驗(yàn)的分段結(jié)果．

利用構(gòu)造應(yīng)力場的均勻性對板塊邊界或斷層進(jìn)行分段的方法假設(shè)板塊邊界或斷層的各子段內(nèi)的應(yīng)力張量方向相同，且構(gòu)造應(yīng)力場方向在不同子段呈顯著變化，用單個震源機(jī)制解與參考應(yīng)力張量之間的殘差作為板塊邊界或斷層的分段參數(shù)(Wyss，Lu，1995)．基于上述假設(shè)，將發(fā)生在三維空間里的地震事件視為沿板塊邊界或斷層走向分布的一維事件，計算出沿斷層或板塊邊界走向的每個震源機(jī)制解相對于參考應(yīng)力張量的最小三維空間旋轉(zhuǎn)角(Kagan，2007)，即本文的分段參數(shù)．相鄰子段間的差異顯著性采用z檢驗(yàn)的統(tǒng)計方法(Meyer，1975)來計算，z檢驗(yàn)是基于中心極限定理，即從相同母體中抽取的兩個樣本組的平均值之差服從正態(tài)分布．z值是用來檢驗(yàn)來自同一序列的兩組獨(dú)立樣本平均值之間的差異，具體檢驗(yàn)統(tǒng)計公式為

(1)

式中，μ1為子段1的平均殘差，μ2為子段2的平均殘差，s1和s2分別為兩個子段內(nèi)的殘差標(biāo)準(zhǔn)差，n1和n2分別為兩個子段內(nèi)的震源機(jī)制解個數(shù)．計算不同參考應(yīng)力張量相對于唐山地震序列震源機(jī)制解的殘差累積的z值，當(dāng)z值超過z檢驗(yàn)顯著性水平為0.05的臨界值(1.96)時，判定其通過顯著性檢驗(yàn)．本文要求當(dāng)2個或2個以上的參考應(yīng)力張量通過差異顯著性檢驗(yàn)時，才可將該震源機(jī)制解所處位置確定為唐山地震序列發(fā)震斷層的分段點(diǎn)．

3 唐山地震序列分段

3.1 參考應(yīng)力張量反演

利用構(gòu)造應(yīng)力場均勻性對斷層進(jìn)行分段，可通過3種方法求取所需要的參考應(yīng)力張量．其一，不同參考應(yīng)力張量與不同子段中震源機(jī)制解的殘差通常不等，當(dāng)所選的參考應(yīng)力張量與某一子段的應(yīng)力場相近時，該子段的殘差累積曲線的斜率變化不明顯，因而易于分辨；其二，根據(jù)斷裂帶的產(chǎn)狀變化進(jìn)行分段，類似于地質(zhì)上的活斷層分段方法；其三，根據(jù)地震活動叢集性進(jìn)行分段．對復(fù)雜的斷裂帶進(jìn)行分段時，需全面考慮各種因素，盡可能選用多個應(yīng)力張量作為參考應(yīng)力張量(Luetal，1997)．關(guān)于唐山地區(qū)震源機(jī)制解的分區(qū)特征和發(fā)震斷層的分段性，目前已有很多研究成果，本文將基于前人根據(jù)震源機(jī)制解得到的分段結(jié)果，定量地對其進(jìn)行差異顯著性檢驗(yàn)，并在此基礎(chǔ)上對唐山地震序列進(jìn)行更加合理的分段．

根據(jù)唐山地區(qū)震源機(jī)制解的分區(qū)特征(許忠淮，1985；李軼群，王健，2008；張宏志等，2008)和唐山地震序列地震斷層面的擬合結(jié)果(萬永革等，2008)，本文將研究區(qū)分為5個子段：寧河段、唐山斷裂帶南段、唐山斷裂帶北段、灤縣段、盧龍斷裂段．將基于震源機(jī)制應(yīng)力反演程序(focal mechanism stress inversion package，簡寫為FMSI)(Gephart，F(xiàn)orsyth，1984)獲取的最佳應(yīng)力模型平均殘差作為表征區(qū)域內(nèi)應(yīng)力場均勻程度的一個指標(biāo)，當(dāng)平均殘差介于6°—9°之間時，反演所得到的最佳應(yīng)力模型視為可接受解，但區(qū)域應(yīng)力場仍存在一定程度的不均勻性(Wyss，Lu，1995；Luetal，1997)．利用FMSI方法分別對研究區(qū)的5個子段進(jìn)行應(yīng)力張量反演，結(jié)果列于表1．可以看出：除唐山斷裂帶北段外，其它4段的應(yīng)力張量均可視為可接受解；寧河、唐山斷裂帶南段和盧龍斷裂這3個子段的最大主壓應(yīng)力方向?yàn)镋SE--WNW向；唐山斷裂帶北段和灤縣段這兩個子段的最大主壓應(yīng)力方向近EW向，且表現(xiàn)出一定的拉張分量．表1還顯示出：唐山斷裂帶北段的平均殘差大于9°，表明該段的應(yīng)力場具有較高的非均勻性；唐山斷裂帶北段的震源機(jī)制解個數(shù)相對于其它4段最多，各種類型的震源機(jī)制解沿唐山斷裂帶北段均勻分布，并且在該斷裂帶東北端形成一條近NS向的地震條帶，表明此處可能存在隱伏于地殼深處的活動斷裂(張四昌，刁桂苓，1992)，復(fù)雜的斷裂結(jié)構(gòu)互相作用可能導(dǎo)致應(yīng)力場的非均勻性．

表1 唐山地震序列發(fā)震斷層分段的參考應(yīng)力張量反演結(jié)果

注：R表示相對應(yīng)力大小，R=(σ2-σ1)/(σ3-σ1).

3.2 基于應(yīng)力場均勻性的斷層分段

將所選用的震源機(jī)制解沿斷層走向從SW至NE進(jìn)行排列，分別計算5個參考應(yīng)力張量和排序后的每個震源機(jī)制解的最小空間旋轉(zhuǎn)角，并將其作為每個地震震源機(jī)制解與參考應(yīng)力張量之間的殘差，繪制殘差累積曲線和z檢驗(yàn)絕對值曲線，如圖2所示．由殘差累積曲線(圖2a)可以看出，在第31，145和219次地震處，寧河段、唐山斷裂帶南段、唐山斷裂帶北段和盧龍斷裂段的參考應(yīng)力張量的殘差累積曲線斜率發(fā)生改變，而灤縣段的參考應(yīng)力張量殘差累積曲線的整體趨勢變化則不明顯．由相應(yīng)殘差累積曲線的z檢驗(yàn)值曲線(圖2b)可以看出：在第31次地震處，寧河段、唐山斷裂帶北段和盧龍斷裂段參考應(yīng)力張量的殘差累積曲線z檢驗(yàn)值均超過臨界值；在第145次地震處，灤縣段和盧龍斷裂段參考應(yīng)力張量的殘差累積曲線z檢驗(yàn)值均超過臨界值；在第219次地震處，唐山斷裂帶北段和盧龍斷裂段參考應(yīng)力張量的殘差累積曲線z檢驗(yàn)值均超過臨界值．對于唐山斷裂帶南段參考應(yīng)力張量而言，各子段間的殘差累積曲線z檢驗(yàn)值未通過統(tǒng)計檢驗(yàn)，且殘差累積z檢驗(yàn)值曲線的趨勢變化與寧河段參考應(yīng)力張量的殘差累積z檢驗(yàn)值曲線大體一致．按照差異顯著性的選擇標(biāo)準(zhǔn)，本文選擇第31，145和219次地震的位置(圖1中壓縮區(qū)用紅色填充的震源機(jī)制解所示)作為唐山地震序列發(fā)震斷層分段的分界點(diǎn)，即由z檢驗(yàn)值結(jié)果可見，在0.05的顯著性水平上，各子段的差異較為顯著．如圖1所示，第31次地震位于唐山斷裂帶與豐臺—豐南斷裂的交匯處，第145次地震位于唐山斷裂帶與一條隱伏斷裂的交匯處，第219次地震位于盧龍斷裂上，因此從構(gòu)造應(yīng)力場均勻性的角度出發(fā)，將唐山地震序列的發(fā)震斷層分為4個子斷層段，第31，145和219次地震的震源位置即為唐山地震序列發(fā)震斷層在地殼深部介質(zhì)力學(xué)性質(zhì)的分界點(diǎn)．

圖2 累積殘差曲線(a)及其相應(yīng)的z檢驗(yàn)絕對值曲線(b)垂直虛線為分段界限，圖(b)中的水平線表示z檢驗(yàn)顯著水平為0.05的臨界值1.96

大地震由于破裂尺度太大對應(yīng)力場方向的局部變化不敏感，因而只能反映大尺度的區(qū)域應(yīng)力場，而小地震既能反映大尺度的應(yīng)力場，又能對局部擾動的應(yīng)力場有所顯示．現(xiàn)今唐山地區(qū)以中小地震(2.5≤ML≤5.1)為主，對應(yīng)的破裂尺度約為幾百米，能夠反映出其應(yīng)力場在10—100 km尺度上是否具有均勻性(Luetal，1997)，因此該震級范圍適于研究唐山地震序列的斷層分段性，即根據(jù)應(yīng)力場的均勻性可以確定各段分界點(diǎn)的位置．根據(jù)上述分段方法，本文將唐山地震序列的發(fā)震斷層分為寧河、唐山、灤縣和盧龍等4個子段，并利用FMSI方法反演得到各子段的最佳應(yīng)力張量，具體列于表2，圖3給出了重新分段后的4個子段的應(yīng)力場3個主應(yīng)力軸的分布圖．

表2 基于應(yīng)力場均勻性得到的唐山地震序列發(fā)震斷層4個子段的應(yīng)力張量反演結(jié)果

圖3 重新分段后4個子段的應(yīng)力場3個主應(yīng)力軸分布圖

從表2所示的反演結(jié)果可以看出：唐山子段的擬合殘差較大，應(yīng)力場非均勻程度較高，這可能與唐山斷裂帶NE端存在一近NS向的地震條帶有關(guān)，將該條帶與唐山斷裂帶北段分別進(jìn)行應(yīng)力場反演，結(jié)果也列于表2，可以看到: 二者的應(yīng)力場都較為均勻；灤縣子段的相對應(yīng)力大小R值為0.5，表明灤縣地區(qū)的應(yīng)力場比較均勻，分段較為合理．從唐山子段分割出來的隱伏斷裂的最佳應(yīng)力張量與灤縣MS7.1地震(圖1中發(fā)震時間為1976年7月28日18時45分的震源機(jī)制解)的矩心矩張量解比較接近，呈純正斷層錯動．由圖3可以看出，4個子段的最大主壓應(yīng)力方向均為近EW向，與李瑞莎等(2008)和盛書中等(2013)對該區(qū)應(yīng)力場的研究結(jié)果基本一致．

4 討論與結(jié)論

本文使用唐山地震序列近期的255個震源機(jī)制解，將該序列的發(fā)震斷層劃分為寧河、唐山、灤縣、盧龍等4個子段．灤縣、盧龍子段的范圍大致上與李軼群和王健(2008)給出的灤縣、盧龍區(qū)一致，本文獲取的灤縣、盧龍子段的最佳應(yīng)力張量與李軼群和王健(2008)的綜合斷層面解的類型基本一致，只是盧龍子段的最大主壓應(yīng)力方向與李軼群和王健(2008)所得的結(jié)果相比有所偏轉(zhuǎn)；寧河、唐山子段的范圍與許忠淮(1985)給出的寧河、唐山區(qū)基本一致，其最大主壓應(yīng)力方向相對于該序列早期有所偏轉(zhuǎn)，最大主應(yīng)力軸傾角大于45°，具有較大的拉張分量．唐山子段與寧河子段的分界點(diǎn)(第31次地震)在EW向的豐臺—豐南斷裂附近，與震源破裂過程反演結(jié)果(Náběleketal，1987；Xie，Yao，1991；Huang，Yeh，1997)中的南北兩段分界點(diǎn)、地表破裂帶南北兩段分界點(diǎn)(郭慧等，2011a，b)的位置大致吻合，唐山子段與唐山斷裂帶北段展布范圍基本一致；而震源破裂過程反演結(jié)果中的南段和地表破裂帶的南段近期地震活動較弱，由于震級較小地震的震源機(jī)制解可信度不高，所以南段可用的震源機(jī)制解個數(shù)較少，且大部分集中在薊運(yùn)河斷裂與昌黎—寧河斷裂交匯處，寧河子段與唐山斷裂帶南段的最佳應(yīng)力張量基本相同，故將唐山斷裂帶南段合并到寧河子段，這樣反演的最佳應(yīng)力張量與前人的研究結(jié)果(許忠淮，1985；李軼群，王健，2008；萬永革等，2008)較為一致．在唐山斷裂帶NE端，一條近NS向的小震條帶(張四昌，刁桂苓，1992；劉亢等，2015)將唐山子段與灤縣子段隔開，唐山子段的北端分界點(diǎn)(第145次地震)恰好就在近NS向的小震條帶上；灤縣子段處于盧龍斷裂與灤縣—樂亭斷裂的共軛交匯處，地震活動叢集，張宏志等(2008)和萬永革等(2008)使用雙差定位后的震源位置結(jié)合斷裂構(gòu)造細(xì)分出盧龍斷裂段和灤縣地震段，但本文認(rèn)為不能將地震活動簡單地歸并至某條斷裂上，因此本文將交匯處的震源機(jī)制解看成一個整體，劃分出灤縣子段，其最佳應(yīng)力張量表現(xiàn)為走滑型的應(yīng)力結(jié)構(gòu)．綜上所述，本文根據(jù)構(gòu)造應(yīng)力場的均勻性所得到的唐山地震序列發(fā)震斷層的分段與震源破裂過程反演結(jié)果、地表破裂分段特征、地震活動叢集以及構(gòu)造特征大體一致，在一定程度上反映了唐山地震震源區(qū)應(yīng)力場的現(xiàn)今狀態(tài)．

唐山地震序列發(fā)震斷層各子段的最佳應(yīng)力張量與該區(qū)較大地震的矩心矩張量解類似，也從一定程度上驗(yàn)證了本文給出的斷層分段較為合理．灤縣MS7.1和寧河MS6.9地震以及絕大多數(shù)余震均發(fā)生在唐山斷裂帶的SW端和NE端，而這兩個區(qū)域恰好是唐山MS7.8地震所產(chǎn)生的引張應(yīng)力區(qū)；本文中寧河子段和唐山子段的最佳應(yīng)力張量為拉張狀態(tài)，與主震對寧河、唐山兩個子段的應(yīng)力場作用方式一致．近期震源機(jī)制解中正斷型地震較多，唐山、寧河和盧龍子段的應(yīng)力場表現(xiàn)出的拉張分量，使我們不能僅限于二維平面討論唐山地區(qū)的應(yīng)力場作用方式，而忽略掉地殼深部構(gòu)造環(huán)境所產(chǎn)生的垂直力作用．劉啟元等(2007)給出的唐山及其鄰區(qū)60 km 深度范圍內(nèi)的三維地殼上地幔S波速度結(jié)構(gòu)顯示，沿唐山斷裂帶布設(shè)的臺站下方的中上地殼存在明顯的向上擠出的S波低速體，流體和熱作用導(dǎo)致并加速了地殼物質(zhì)的弱化和物質(zhì)向兩側(cè)的擠出，同時在地殼上部產(chǎn)生張應(yīng)力，這可能是唐山子段近期走滑型和正斷型事件活躍的深部動力學(xué)成因；在剖面NE端存在的局部低速體可能與盧龍子段的應(yīng)力場作用方式有關(guān)，而灤縣子段在剖面上則未顯示出較強(qiáng)的橫向不均勻性．唐山地震序列的形成與地殼內(nèi)部垂直構(gòu)造差異運(yùn)動和張剪性局部應(yīng)力場的觸發(fā)作用密切相關(guān)(Shaoetal，1986；劉啟元等，2007)．

從近期唐山地區(qū)應(yīng)力場的基本特征可以看出，唐山地震的震源斷層并沒有愈合，而是處于持續(xù)活動狀態(tài)，且現(xiàn)今的中小地震具有繼承性斷裂活動特征．寧河、唐山、灤縣和盧龍這4個子段相對之前研究結(jié)果所表現(xiàn)出來的最大主壓應(yīng)力方向偏轉(zhuǎn)和傾伏角增大等現(xiàn)象需要基于更豐富、更精確的資料進(jìn)行精細(xì)研究．隨著數(shù)字地震觀測技術(shù)的進(jìn)步和數(shù)字地震臺網(wǎng)的廣泛布設(shè)，高精度地震觀測資料為我們研究歷史大震震源構(gòu)造的精細(xì)結(jié)構(gòu)提供了數(shù)據(jù)基礎(chǔ)，使得重新考慮唐山、灤縣、寧河地震的破裂模型、機(jī)制及其震后的動力學(xué)演化過程成為可能．本文僅從構(gòu)造應(yīng)力場均勻性的角度對唐山地震序列的發(fā)震斷層進(jìn)行分段，期望下一步能夠加入更小震級地震的震源機(jī)制解資料，對唐山地震序列的發(fā)震斷層再進(jìn)行下一級子段的細(xì)分．

刁桂苓，于利民，李欽祖，王紹晉. 1995. 唐山和瀾滄地震序列震源區(qū)應(yīng)力場的對比分析[J]. 地震學(xué)報，17(3)：305--311.

Diao G L，Yu L M，Li Q Z，Wang S J. 1995. Comparison of the stress field between Tangshan and Lancang earthquake sequence[J].ActaSeismologicaSinica，17(3): 305--311 (in Chinese).

顧功敘，林庭煌，時振梁. 1983. 中國地震目錄(公元1970—1979年)[M]. 北京：地震出版社: 110-190.

Gu G X，Lin T H，Shi Z L. 1983.CatalogueofChineseEarthquakes(AD 1970--1979)[M]. Beijing: Seismological Press: 110--190 (in Chinese).

郭慧，江娃利，謝新生. 2011a. 鉆孔與探槽揭示1976年河北唐山MS7.8地震發(fā)震構(gòu)造晚第四紀(jì)強(qiáng)震活動[J]. 中國科學(xué)：地球科學(xué)，41(7)：1009--1028.

Guo H，Jiang W L，Xie X S. 2011a. Late-Quaternary strong earthquakes on the seismogenic fault of the 1976MS7.8 Tangshan earthquake，Hebei，as revealed by drilling and trenching[J].ScienceChinaEarthSciences，54(11): 1696--1715.

郭慧，江娃利，謝新生. 2011b. 對1976 年河北唐山MS7.8地震地表破裂帶展布及位移特征的新認(rèn)識[J]. 地震地質(zhì)，33(3)：506--524.

Guo H，Jiang W L，Xie X S. 2011b. New evidence for the distribution of surface rupture zone of the 1976MS7.8 Tang-shan earthquake[J].SeismologyandGeology，33(3): 506--524 (in Chinese).

國家地震局《一九七六年唐山地震》編輯組. 1982. 一九七六年唐山地震[M]. 北京：地震出版社：1--32.

Editorial Group ofTheTangshanEarthquakeof1976，State Seismological Bureau. 1982.TangshanEarthquakeof1976[M]. Beijing: Seismological Press: 1--32 (in Chinese).

虢順民，李志義，程紹平，陳獻(xiàn)程，陳孝德，楊主恩，李如成. 1977. 唐山地震區(qū)域構(gòu)造背景和發(fā)震模式的討論[J]. 地質(zhì)科學(xué)，12(4)：305--321.

Guo S M，Li Z Y，Cheng S P，Chen X C，Chen X D，Yang Z E，Li R C. 1977. Discussion on the regional structural background and the seismogenic model of the Tangshan earthquake[J].ScientiaGeologicaSinica，12(4): 305--321 (in Chinese).

李欽祖，刁桂苓，戴英華. 1983. 唐山地震序列的應(yīng)力釋放調(diào)整過程[J]. 地球物理學(xué)報，26(3)：224--236.

Li Q Z，Diao G L，Dai Y H. 1983. The stress releases and adjustment process during Tangshan earthquake sequence[J].ChineseJournalofGeophysics，26(3): 224--236 (in Chinese).

李瑞沙，崔效鋒，刁桂苓，張紅艷. 2008. 華北北部地區(qū)現(xiàn)今應(yīng)力場時空變化特征研究[J]. 地震學(xué)報，30(6)：570--580.

Li R S，Cui X F，Diao G L，Zhang H Y. 2008. Temporal and spatial variation of the present crustal stress in northern part of North China[J].ActaSeismologicaSinica，30(6): 570--580 (in Chinese).

李軼群，王健. 2008. 唐山余震區(qū)中小地震震源機(jī)制解分區(qū)特征的初步研究[J]. 中國地震，24(2)：150--158.

Li Y Q，Wang J. 2008. The preliminary study on focal mechanism of small and moderate earthquakes in divided zones around Tangshan[J].EarthquakeResearchinChina，24(2): 150--158 (in Chinese).

梁尚鴻，李幼銘，束沛鎰，朱碚定．1984．利用區(qū)域地震臺網(wǎng)P、S振幅比資料測定小震震源參數(shù)[J]．地球物理學(xué)報，27(3)：249--257.

劉桂萍，傅征祥. 2000. 1976年7月28日唐山7.8級地震觸發(fā)的區(qū)域地震活動和靜應(yīng)力場變化[J]. 地震學(xué)報，22(1)：17--26.

Liu G P，F(xiàn)u Z X. 2000. The seismicity and static stress change of the trigger area of 1976 Tangshan earthquake[J].ActaSeismologicaSinica，22(1): 17--26 (in Chinese).

劉亢，曲國勝，房立華，李紅光. 2015. 唐山古冶、灤縣地區(qū)中小地震活動與構(gòu)造關(guān)系研究[J]. 地震，35(2)：111--120.

Liu K，Qu G S，F(xiàn)ang L H，Li H G. 2015. Relationship between tectonics and small to moderate earthquakes in Guye and Luanxian area，Hebei Province[J].Earthquake，35(2): 111--120 (in Chinese).

劉啟元，王峻，陳九輝，李順成，郭彪. 2007. 1976 年唐山大地震的孕震環(huán)境: 密集地震臺陣觀測得到的結(jié)果[J]. 地學(xué)前緣，14(6)：205--213.

Liu Q Y，Wang J，Chen J H，Li S C，Guo B. 2007. Seismogenic tectonic environment of 1976 great Tangshan earthquake: Results given by dense seismic array observations[J].EarthScienceFrontiers，14(6): 205--213 (in Chinese).

馬瑾，張渤濤，袁淑榮. 1980. 唐山地震與地震危險區(qū)[J]. 地震地質(zhì)，2(2)：43--54.

Ma J，Zhang B T，Yuan S R. 1980. Tangshan earthquake and the earthquake risk area[J].SeismologyandGeology，2(2): 43--54 (in Chinese).

盛書中，萬永革．2012．由構(gòu)造應(yīng)力場研究汶川地震斷層的分段性[J]．地震學(xué)報，34(6)：741--753.

Sheng S Z，Wan Y G. 2012. Segmentation of the Wenchuan earthquake fault derived from tectonic stress analysis[J].ActaSeismologicaSinica，34(6): 741--753 (in Chinese).

盛書中，萬永革，徐志國，卜玉菲，武曄，李迎秋. 2013. 由地震釋放的地震矩疊加推導(dǎo)平均應(yīng)力場[J]. 地震地質(zhì)，35(1)：92--100.

Sheng S Z，Wan Y G，Xu Z G，Bu Y F，Wu Y，Li Y Q. 2013. Mean stress field inferred from the total seismic moment released by earthquakes[J].SeismologyandGeology，35(1): 92--100 (in Chinese).

萬永革，沈正康，刁桂苓，王福昌，胡新亮，盛書中. 2008. 利用小震分布和區(qū)域應(yīng)力場確定大震斷層面參數(shù)方法及其在唐山地震序列中的應(yīng)用[J]. 地球物理學(xué)報，51(3)：793--804.

Wan Y G，Shen Z K，Diao G L，Wang F C，Hu X L，Sheng S Z. 2008. An algorithm of fault parameter determination using distribution of small earthquakes and parameters of regional stress field and its application to Tangshan earthquake sequence[J].ChineseJournalofGeophysics，51(3): 793--804 (in Chinese).

萬永革，沈正康，曾躍華，盛書中，徐曉楓. 2008. 唐山地震序列應(yīng)力場觸發(fā)的粘彈性力學(xué)模型研究[J]. 地震學(xué)報，30(6)：581--593.

Wan Y G，Shen Z K，Zeng Y H，Sheng S Z，Xu X F. 2008. Study on visco-elastic stress triggering model of the 1976 Tangshan earthquake sequence[J].ActaSeismologicaSinica，30(6): 581--593 (in Chinese).

許忠淮. 1985. 用滑動方向擬合法反演唐山余震區(qū)的平均應(yīng)力場[J]. 地震學(xué)報，7(4)：349--362.

Xu Z H. 1985. Mean stress field in Tangshan aftershock area obtained from focal mechanism data by fitting slip directions[J].ActaSeismologicaSinica，7(4): 349--362 (in Chinese).

許忠淮，汪素云，高阿甲，郭瑛. 1994. 我國部分早期震源機(jī)制解的重新測定[J]. 地震地磁觀測與研究，15(5)：1--9.

Xu Z H，Wang S Y，Gao A J，Guo Y. 1994. Redetermination of some early focal mechanism solutions of Chinese earthquakes[J].SeismologicalandGeomagneticObservationandResearch，15(5): 1--9 (in Chinese).

尤惠川，徐錫偉，吳建平，何正勤. 2002. 唐山地震深淺構(gòu)造關(guān)系研究[J]. 地震地質(zhì)，24(4)：571--582.

You H C，Xu X W，Wu J P，He Z Q. 2002. Study on the relationship between shallow and deep structures in the 1976 Tangshan earthquake area[J].SeismologyandGeology，24(4): 571--582 (in Chinese).

張宏志，刁桂苓，陳棋福，胡新亮，王曉山，張躍剛. 2008. 1976年唐山7.8級地震震區(qū)現(xiàn)今地震震源機(jī)制分析[J]. 地震研究，31(1)：1--6.

Zhang H Z，Diao G L，Chen Q F，Hu X L，Wang X S，Zhang Y G. 2008. Focal mechanism analysis of the recent earthquakes in Tangshan seismic region ofM7.8 in 1976[J].JournalofSeismologicalResearch，31(1): 1--6 (in Chinese).

張四昌，刁桂苓. 1992. 唐山地震序列的構(gòu)造過程[J]. 中國地震，8(2)：73--80.

Zhang S C，Diao G L. 1992. The tectonic process of the Tangshan earthquake sequence[J].EarthquakeResearchinChina，8(2): 73--80 (in Chinese).

Butler R，Stewart G S，Kanamori H. 1979. The July 27，1976 Tangshan，China earthquake：A complex sequence of intraplate events[J].BullSeismolSocAm，69(1): 207-220.

Dziewonski A M，Chou T A，Woodhouse J H. 1981. Determination of earthquake source parameters from waveform data for studies of global and regional seismicity[J].JGeophysRes，86(B4): 2825--2852.doi:10.1029/JB086iB04p02825.

Ekstr?m G，Nettles M，Dziewonski A M. 2012. The global CMT project 2004--2010: Centroid-moment tensors for 13017 earthquakes[J].PhysEarthPlanetInter，200/201: 1--9.doi:10.1016/j.pepi.2012.04.002.

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].JGeophysRes，89(B11): 9305--9320.doi:10.1029/JB089iB11p09305.

Huang B S，Yeh Y T. 1997. The fault ruptures of the 1976 Tangshan earthquake sequence inferred from coseismic crustal deformation[J].BullSeismolSocAm，87(4): 1046--1057.

Kagan Y Y. 2007. Simplified algorithms for calculating double-couple rotation[J].GeophysJInt，171(1): 411--418.

Klein F W. 1978.HypocenterLocationProgramHYPOINVERSEPartI:UsersGuidetoVersions1，2，3and4[R]. Menlo Park: U.S. Geological Survey Open-File Rept：78--694.

Lu Z，Wyss M. 1996. Segmentation of the Aleutian plate boundary derived from stress direction estimates based on fault plane solutions[J].JGeophysRes，101(B1): 803--816.

Lu Z，Wyss M，Pulpan H. 1997. Details of stress directions in the Alaska subduction zone from fault plane solutions[J].JGeophysRes，102(B3): 5385--5402.

Meyer S L. 1975.DataAnalysisforScientistsandEngineers[M]. New York: John Wiley & Sons: 235.

Nábělek J，Chen W P，Ye H. 1987. The Tangshan earthquake sequence and its implications for the evolution of the North China basin[J].JGeophysRes，92(B12): 12615--12628.doi:10.1029/JB092iB12p12615.

Ratchkovski N A. 2003. Change in stress directions along the central Denali fault，Alaska after the 2002 earthquake sequence[J].GeophysResLett，30(19): 2017.doi: 10.1029/2003GL017905.

Robinson R，Zhou S Y. 2005. Stress interactions within the Tangshan，China，earthquake sequence of 1976[J].BullSeismolSocAm，95(6): 2501--2505.

Shao X Z，Zhang J R，Liu Q Y，Zhang S Y. 1986. The common features of deep structures in some large earthquake areas of the North China plain and their implication for earthquake site prediction[J].JPhysEarth，34(Suppl): 223--239.

Shedlock K M，Baranowski J，Xiao W W，Hu X L. 1987. The Tangshan aftershock sequence[J].JGeophysRes，92(B3): 2791--2803.

Wyss M，Lu Z. 1995. Plate boundary segmentation by stress directions: Southern San Andreas fault，California[J].GeophysResLett，22(5): 547--550.

Xie X B，Yao Z X. 1991. The faulting process of Tangshan earthquake inverted simultaneously from the teleseismic waveforms and geodesic deformation data[J].PhysEarthPlanetInter，66(3/4): 265--277.

Xu Z H，Wang S Y. 1986. A possible change in stress field orientation due to the 1976 Tangshan earthquake[J].PureApplGeophys，124(4/5): 941--955.

Zoback M L. 1992. First- and second-order patterns of stress in the lithosphere: The World Stress Map Project[J].JGeophysRes，97(B8): 11703--11728.doi:10.1029/92JB00132.

Seismogenic fault segmentation of Tangshan earthquake sequence derived from focal mechanism solutions

Yang Yaqiong1)Wang Xiaoshan1,2),*Wan Yongge3,2)Sheng Shuzhong3)Chen Ting1)

1)EarthquakeAdministrationofHebeiProvince,Shijiazhuang050021,China2)InstituteofGeophysics,ChinaEarthquakeAdministration,Beijing100081,China3)InstituteofDisasterPrevention,HebeiSanhe065201,China

Based on the method of tectonic stress field homogeneity, this paper determined the segmentation of the seismogenic fault of Tangshan earthquake sequence by using the focal mechanisms of earthquake withML≥2.5 occurred in Tangshan area from January of 2002 to November of 2015. According to the previous subarea results related to focal mechanisms in Tangshan area, this paper obtained five reference stress tensors, and the seismogenic fault of Tang-shan earthquake sequence was divided into four segments,i.e., Ninghe, Tang-shan, Luanxian and Lulong, by the significant difference analysis ofz-test of cumulative misfit. The inversion for stress fields of the four segments shows the orientation of maximum principal stress axes is nearly in EW, and the stress filed of the Tangshan, Ninghe and Lulong segments showed a greater tensile component. The best stress tensors of Tangshan and Ninghe segments are consistent with the tensile stress field caused by the Tangshan main shock’s pulling the northeastern and southwestern ends of Tangshan fault zone. The uplift model of upper mantle material derived from receiver function may also conform to the stress field of Tangshan segment. The strike-slip stress field of Luanxian segment is closely related to the conjugated structure in this area. The best stress tensor of Lulong segment is normal fault with right-lateral strike-slip. From the present characteristics of stress field for the seismogenic fault segmentation of Tangshan earthquake sequence, it is concluded that the seismicity in Tangshan area has inherited the development mode of Tangshan earthquake sequence, and mainly can be controlled by the regional tectonic stress field and the coexistence of deep and shallow faults.

Tangshan earthquake sequence; focal mechanism; tectonic stress field; segmentation of the seismogenic fault

10.11939/jass.2016.04.009.

國家自然科學(xué)基金(91214201)和河北省地震局地震科技星火計劃重點(diǎn)項(xiàng)目(DZ20150428102, DZ20140101002)聯(lián)合資助.

2016-04-11收到初稿，2016-06-10決定采用修改稿.

10.11939/jass.2016.04.009

P315.2

A

楊雅瓊，王曉山，萬永革，盛書中，陳婷. 2016. 由震源機(jī)制解推斷唐山地震序列發(fā)震斷層的分段特征. 地震學(xué)報, 38(4): 632--643.

Yang Y Q, Wang X S, Wan Y G, Sheng S Z, Chen T. 2016. Seismogenic fault segmentation of Tangshan earthquake sequence derived from focal mechanism solutions.ActaSeismologicaSinica, 38(4): 632--643.doi:10.11939/jass.2016.04.009.

*通訊作者 e-mail: wxs@eq-he.ac.cn