周依 曹筠 馬旭東 單連君

摘要：為分析唐山老震區(qū)地震活動背景，對該區(qū)域現(xiàn)今地震活動特征進行研究。老震區(qū)現(xiàn)今地震活動與1976年唐山MS7.8地震早期余震密集區(qū)域存在一定差異，結合區(qū)域背景地震活動統(tǒng)計結果，認為該區(qū)域地震活動受唐山地震序列影響逐漸減弱，逐步趨近于背景地震。地震分布呈現(xiàn)出明顯的3分區(qū)差異，主要集中于中區(qū)，東區(qū)次之，西區(qū)最弱。中等地震的空間分布特征顯示，地震活動存在東、中、西區(qū)交替發(fā)震的規(guī)律。目前，唐山老震區(qū)構造應力積累加速，存在地震活動水平由平靜向活躍轉(zhuǎn)變的可能，需要重點關注唐山斷裂帶北段以及灤縣—樂亭斷裂和盧龍斷裂的交匯區(qū)域。

關鍵詞：唐山地震；地震活動特征；背景地震；遷移規(guī)律

中圖分類號：P315.72文獻標識碼：A文章編號：1000-0666（2024）03-0430-11

doi：10.20015/j.cnki.ISSN1000-0666.2024.0046

0引言

1976年7月28日3時42分，唐山發(fā)生MS7.8地震，15個小時之后在其東部的灤縣（現(xiàn)灤州市）發(fā)生MS7.1地震，3個半月后在其西部的寧河再次發(fā)生MS6.9地震，主震及兩次強余震之后形成3個明顯的余震區(qū)，震中最高烈度XI度，造成了嚴重的災害損失（國家地震局《一九七六年唐山地震》編輯組，1982；孫麗娜等，2021）。唐山地震序列復雜且豐富，由7.8級主震、兩次7級左右強余震和一系列5級以上余震構成。隨著時間推移，序列總體呈衰減特征，衰減過程中雖有階段性起伏，但是符合強震的衰減規(guī)律（張躍剛等，2013；郭蕾等，2021）。

截至2022年，唐山地震已過去46年，該地區(qū)現(xiàn)今小震活動仍十分頻繁。唐山老震區(qū)巖石圈黏彈性松弛模擬結果顯示，1976年唐山地震發(fā)生15年內(nèi)庫倫應力變化顯著，而15年之后庫倫應力呈現(xiàn)緩慢的穩(wěn)定增加狀態(tài)，顯示區(qū)域應力調(diào)整已基本穩(wěn)定，因此15年左右的時間尺度可能是唐山余震活動的一個時間分界點（朱琳等，2021）。王想等（2012）認為唐山地區(qū)1992年進入余震活動晚期階段之后，受唐山MS7.8主震和強余震的控制減弱，地震活動趨于穩(wěn)定狀態(tài)。1996年以來唐山老震區(qū)中等地震活動顯示，唐山MS7.8地震震中附近區(qū)域基本沒有4級以上地震，余震活動可能已經(jīng)不受主震引起的局部應力狀態(tài)變化所制約（張素欣等，2020）。鄭建常等（2014，2022）通過對唐山地震序列不同時段進行叢集地震識別，認為2011年之后唐山地震序列受主震和較大余震影響的事件越來越少，更多小震已表現(xiàn)為彼此獨立的隨機事件。以上研究均表明，唐山余震活動進入晚期階段后，地震活動受主震影響逐漸減小，目前該區(qū)域發(fā)生的地震可看作背景地震。

唐山MS7.8地震屬于大陸板塊內(nèi)部地震（Robinson，Zhou，2005）。[HJ2mm]研究顯示，板內(nèi)地震衰減較慢，余震持續(xù)時間可達數(shù)百年（Stein，Liu，2009）。仲秋和史保平（2012）根據(jù)對數(shù)線性回歸和余震觸發(fā)理論計算，得到唐山地震的余震序列還將持續(xù)近100年的結論。唐山老震區(qū)的GPS測量結果顯示，唐山地區(qū)的地殼應變率高于華北地區(qū)的背景值，王輝等（2013）認為這一現(xiàn)象表明該地區(qū)現(xiàn)今仍受1976年唐山地震序列影響。關于唐山老震區(qū)現(xiàn)今地震活動是否為唐山MS7.8地震余震的討論始終沒有明確的結論。唐山地區(qū)如今經(jīng)過40多年的恢復重建，經(jīng)濟發(fā)達、人口稠密，其所處的首都圈地區(qū)，一直是我國地震重點監(jiān)視防御區(qū)（孫啟凱等，2017）。因此，開展唐山老震區(qū)現(xiàn)今地震活動特征研究，分析其是否為該區(qū)域的背景地震活動，對于進行唐山地區(qū)地震風險評估以及助力當?shù)胤勒饻p災建設具有一定的現(xiàn)實意義。

1區(qū)域地質(zhì)構造背景與資料選取

1.1區(qū)域地質(zhì)構造背景

唐山地區(qū)位于張渤地震帶和華北平原地震帶交匯處，受NEE向豐臺—野雞坨斷裂、寧河—昌黎斷裂和NW向灤縣—樂亭斷裂、薊運河斷裂的切割，形成了一個NEE向菱形塊體（虢順民等，1977）。菱形塊體內(nèi)部蓋層發(fā)育有唐山斷裂帶，西南界的薊運河斷裂與塊體西南發(fā)育的NNE向滄東斷裂相交，東邊界的灤縣—樂亭斷裂與NNE向盧龍斷裂相交（圖1）。塊體內(nèi)部的唐山斷裂帶是一條走向NE、具右旋走滑特征的褶皺斷裂帶，由陡河斷裂、巍山—長山南坡斷裂、唐山—古冶斷裂等數(shù)條次級斷裂組成，1976年唐山MS7.8地震序列主要沿唐山斷裂帶展布（徐杰等，2012）。本文研究區(qū)唐山老震區(qū)主要指唐山地震序列的小震密集區(qū)（張素欣等，2017），其深部地球物理探測顯示，在唐山斷裂帶的正下方存在一條傾角陡立，具有走滑性質(zhì)的深大斷裂，被稱為“唐山地殼深斷裂”。該斷裂錯斷了上、下地殼分界和莫霍面，為幔源物質(zhì)的侵入提供了通道（劉保金等，2011）。據(jù)此推測，唐山地區(qū)深、淺部耦合的地殼構造系統(tǒng)可能是該區(qū)域地震活動頻繁的原因（李赫等，2020）。

1.2資料選取

1992年以來，唐山余震活動進入晚期階段（王想等，2012）。據(jù)此，本文基于河北地震臺網(wǎng)

資料，整理了1992—2021年唐山地區(qū)（38.8°～40.5°N，117°～120°E）地震目錄，對該區(qū)域現(xiàn)今地震活動特征進行分析。為保證分析結果的可靠性，采用震級-序號法和多統(tǒng)計方法（蔣長勝等，2013）對最小完整性震級MC進行評估。其中，多統(tǒng)計方法是使用最大曲率法（MAXC）和擬合優(yōu)度法（GFT-90%、GFT-95%）的一種綜合評估方法。如圖2所示，兩種方法計算得到的結果較為一致，1992年以來，唐山老震區(qū)最小完整性震級MC為ML2.0左右；2008年以來，隨著“十五”數(shù)字地震觀測網(wǎng)絡項目投入運行，研究區(qū)的地震監(jiān)測能力顯著提升，MC已低于ML1.5；2015年以后，得益于首都圈測震臺站升級改造項目的完成，研究區(qū)地震監(jiān)測能力進一步提升，MC達到ML1.0。因此，本文分析唐山老震區(qū)現(xiàn)今地震活動特征時，將地震起始震級設定為ML2.0。

2小震密集分布特征分析

根據(jù)1992年以來唐山老震區(qū)最小完整性震級MC，設置起始震級為ML2.0，繪制區(qū)域地震震中分布圖。如圖3a所示，地震表現(xiàn)出明顯的密集分布特征，主要沿唐山斷裂帶和盧龍斷裂呈NE向密集分布，沿著薊運河斷裂也存在少量地震呈NW向帶狀分布。在盧龍斷裂和灤縣—樂亭斷裂交匯區(qū)域附近，地震密集活動出現(xiàn)“丁”字形轉(zhuǎn)折。為進一步分析地震活動的時空分布特征，沿地震密集區(qū)繪制了AB軸，長軸方向NE。在分析小震時空分布圖像時，設定A點為研究區(qū)軸線原點，根據(jù)距離軸線20 km范圍內(nèi)地震相對原點的距離，得到唐山地區(qū)小震時空演化圖像（圖3b）。地震密集區(qū)集中在距離A點50～110 km處，該區(qū)域主要為唐山斷裂帶北段和盧龍斷裂南段，大部分ML4.0以上地震都發(fā)生在該區(qū)域，其中包括所有ML5.0以上地震，其余的ML4.0以上地震主要分布在薊運河斷裂與唐山斷裂帶交匯附近。2004年1月20日昌黎ML5.0地震后，雖然小震活動仍然集中于該區(qū)域，但是地震密集范圍明顯收縮，更為集中在唐山斷裂帶北端和盧龍斷裂與灤縣—樂亭斷裂的交匯區(qū)域，使得在唐山斷裂帶和灤縣—樂亭斷裂交而不匯的位置出現(xiàn)了一條清晰的地震分界線。

為了更直觀地反映唐山老震區(qū)小震分布情況，區(qū)分地震活動相對密集地區(qū)與稀疏地區(qū)，本文基于小震震中分布繪制唐山地區(qū)地震密集分布圖。對研究區(qū)域進行0.1°×0.1°的網(wǎng)格劃分，統(tǒng)計網(wǎng)格節(jié)點附近10 km范圍內(nèi)的地震次數(shù)并對其進行歸一化處理，由此得到的密集程度值能夠直觀描述研究區(qū)域的地震密集程度和密集區(qū)域。從小震密集分布特征來看，唐山老震區(qū)現(xiàn)今地震活動存在明顯差異，呈現(xiàn)空間分布不均勻特征。結合區(qū)域構造環(huán)境并參考前人的分區(qū)規(guī)則（楊理華，趙喜柱，1985；張宏志等，2008），本文將唐山老震區(qū)劃分為東區(qū)、中區(qū)和西區(qū)，其中灤縣—樂亭斷裂及以東為東區(qū)，唐山斷裂帶所在的菱形塊體內(nèi)部為中區(qū)，薊運河斷裂及以西為西區(qū)。如圖4a所示，1992年以來，中區(qū)地震密集程度尤為突出，地震密集區(qū)域主要集中在唐山斷裂帶北段，其南段地震分布相對稀疏；東區(qū)地震密集程度次之，主要集中在灤縣—樂亭斷裂與盧龍斷裂交匯區(qū)域，盧龍斷裂北段有少量地震聚集；西區(qū)地震活動沿著薊運河斷裂呈NW向分布，地震密集程度較弱。將1976年唐山MS7.8主震及兩次強余震的震中位置投影到地震密集分布圖上發(fā)現(xiàn)，唐山MS7.8和灤縣MS7.1地震震中均位于現(xiàn)今地震密集區(qū)域的邊緣附近，寧河MS6.9地震震中位于現(xiàn)今地震活動較弱區(qū)域。而1992年以前的余震密集分布圖顯示（圖4b），唐山MS7.8和灤縣MS7.1地震震中均處于地震密集程度較高區(qū)域，且寧河MS6.9地震震中也位于地震密集區(qū)域的邊緣。與1992年以前相比，現(xiàn)今地震密集區(qū)域收縮，偏離3次強震震中，可能預示區(qū)域地震受唐山MS7.8強震序列影響減弱。

3中等地震時空分布特征分析

3.1時間叢集特征

區(qū)域地震活動時間分布特征研究是在一定統(tǒng)計區(qū)內(nèi)（地震帶或地震區(qū)內(nèi)）進行的，描述在時間進程中不同時間段的地震活動特征。1992年以來唐山老震區(qū)ML4.0以上地震累積應變能量釋放曲線整體表現(xiàn)比較平緩，未出現(xiàn)較大波動，說明整體地震活動比較平穩(wěn)。從細節(jié)上看，累積應變能量釋放曲線在2004年和2019年出現(xiàn)了小的臺階，參考1992年以來唐山老震區(qū)ML4.0以上地震的M-t圖，可以將該區(qū)域現(xiàn)今中等地震活動大致分為3個階段（圖5）。第1階段1992—2004年，累積應變能量釋放用J來表示，其速率為11.2 J/a，ML4.0以上地震平均年頻次為1.7次。期間發(fā)生了3次ML5.0以上地震，最大地震為1995年10月6日唐山ML5.4。該階段地震發(fā)震時間間隔多集中在0.5～1 a，最長發(fā)震間隔為2 a，期間多次出現(xiàn)短期（3個月內(nèi)）地震成對發(fā)生。第2階段2005—2018年，累積應變能量釋放速率和地震平均年頻次較第1階段有所降低，分別為6.0 J/a和0.9次。期間發(fā)生ML5.0以上地震1次，為2012年5月28日唐山ML5.1。該階段地震存在2 a左右的準周期特征，并且“雙震”特征更加顯著。但2015年9月14日昌黎ML4.8地震擾亂了2 a左右的準周期特征，此后發(fā)震時間間隔縮短，此時可能處于這1階段末期的剩余應變釋放調(diào)整階段。2016年9月10日唐山ML4.6地震后，唐山老震區(qū)進入了長達3.2 a的ML4地震平靜，直至2019年12月5日豐南ML4.9地震的發(fā)生，預示著進入了地震活動的第3階段。目前處于第3階段早期，共計發(fā)生3次ML4.0以上地震，最大地震為2020年7月12日古冶ML5.4地震，該階段地震發(fā)震時間間隔保持在0.5～1.0 a，目前累積應變能量釋放速率為10.5 J/a。

前2個階段的持續(xù)時間均在10 a左右，第2階段累積應變能量釋放速率、地震平均年頻次和最大地震震級均低于第1階段，符合老震區(qū)地震衰減規(guī)律。但是從第3階段唐山老震區(qū)目前發(fā)生的中等地震活動情況來看，累積應變能量釋放速率和最大地震震級均有所回升，ML4.0以上地震的發(fā)震時間間隔縮短，存在地震活動水平由平靜向活躍轉(zhuǎn)變的可能，與衰減趨勢相悖，可能意味著該區(qū)域地震的發(fā)生受唐山MS7.8強震序列約束進一步減弱。

3.2空間遷移特征

1976年唐山MS7.8主震及兩次強余震依次發(fā)生于3分區(qū)中的中、東和西區(qū)，存在空間遷移特征（圖4）。眾多專家學者已對其成因及動力機制進行了闡述（Butler et al，1979；國家地震局《一九七六年唐山地震》編輯組，1982；聞學澤，馬勝利，2006）。本文對唐山老震區(qū)現(xiàn)今中等地震空間分布特征進行分析時，發(fā)現(xiàn)這一特征同樣顯著。

將1992年以來唐山老震區(qū)ML4.0以上地震分組討論，若以發(fā)震時間間隔超過2.0 a作為分組條件，則大致可以分為5組，3個階段分別存在2組、2組和1組地震活動。從地震遷移曲線上看（圖6），存在東、中、西3個分區(qū)交替發(fā)震的特點。5組地震活動中有3組的始發(fā)地震出現(xiàn)在中區(qū)，即遷移過程大多從中區(qū)開始，而結束地震多發(fā)生在東區(qū)。從遷移特征來看，雖未發(fā)現(xiàn)較明顯的優(yōu)勢遷移方向，但遷移過程中存在一定的有序性，當?shù)卣鸹顒影唇?jīng)度從東向西或從西向東活動時，會在3個分區(qū)依次進行遷移，每個遷移過程中區(qū)都會參加。在19次不同分區(qū)的地震遷移過程中，僅有3次沒有經(jīng)過中區(qū)，其中有2次從東區(qū)直接跳躍到西區(qū)（第1、4組），還有1次則從西區(qū)跳躍到東區(qū)（第1組）。由此可以將中區(qū)作為唐山老震區(qū)地震活動的窗口。研究表明，在東、中、西3個分區(qū)均存在應力集中點（馬瑾等，1980；劉蒲雄，呂曉健，2011），當區(qū)域應力積累到一定程度后，在其中一個應力集中點進行應力釋放。1976年唐山7.8級地震后，由于老震區(qū)地質(zhì)介質(zhì)相對破碎，介質(zhì)強度不易積累更大能量（Li et al，2007），較難阻止地震波能量的傳播，初始破裂的沖擊波誘發(fā)了其余兩個應力集中點的應力釋放，導致發(fā)震。3個分區(qū)中等地震空間遷移的現(xiàn)象，就是在3個應力集中點依次應力釋放的過程。

若以遍歷東、中、西區(qū)為一重遷移，再一次遍歷為二重遷移，則前2組地震活動存在二重遷移，后3組地震活動僅存在一重遷移。5組地震活動存在的7次遷移過程中，遷移路徑為“西→中→東”模式的情況最多（圖7a、b、e），“東→中→西”模式次之（圖7b、c）。還有2次遷移過程為跨區(qū)跳躍的模式：第1組地震活動的第一重遷移過程中，天津發(fā)生ML4.9地震后，直接在東區(qū)發(fā)生了灤縣ML4.0地震（圖7a）；第4組地震活動遷移過程中，東區(qū)發(fā)生昌黎ML4.7地震后，接著西區(qū)發(fā)生豐南ML4.0地震（圖7d）。但是從統(tǒng)計來看，跨區(qū)跳躍的情況較少發(fā)生，遷移過程主要遵循3個分區(qū)依次遷移的規(guī)律。

研究時段內(nèi)最后一次ML4.0以上地震為2021年4月16日灤州ML4.8，位于東區(qū)。自此，第5組地震活動已經(jīng)完成了一重遷移，遷移模式為“西→中→東”。參考前幾次該模式下的遷移過程，在第1組地震活動的第二重遷移過程中，中區(qū)依次發(fā)生了唐山ML4.0和ML5.4地震后，再往東區(qū)遷移；第2組地震活動的第二重遷移過程，在東區(qū)接連發(fā)生了昌黎ML4.1和ML5.0兩次結束地震?？紤]到“西→中→東”遷移模式下“雙震”現(xiàn)象的存在，第5組地震遷移過程有可能還沒有結束，即下一次ML4.0以上地震存在發(fā)生在東區(qū)的可能。

4討論

4.1現(xiàn)今地震活動性質(zhì)判定

1976年唐山MS7.8地震早期余震主要分布在主震及兩次強余震附近（Shedlock et al，1987；張素欣等，2017，2020），且1992年前、后的余震密集分布也表現(xiàn)出不同特征（圖4）。許忠淮等（1983，1985）利用唐山地震前10年的小震P波初動數(shù)據(jù)，得到震前區(qū)域主壓應力方位約為N60°E，而震后早期表現(xiàn)為近EW方向，順時針旋轉(zhuǎn)了約30°。之后，唐山老震區(qū)的主壓應力方向長期處在近EW向（李瑞莎等，2008）。但是，2002年以來的地殼應力場反演結果和GPS應變場顯示，該區(qū)域目前主要受NEE向的主壓應變控制（陳長云，2016；王曉山等，2020）。表明唐山老震區(qū)現(xiàn)今主壓應力方向較1976年唐山地震序列早期可能發(fā)生了調(diào)整，趨于恢復到震前水平。而唐山老震區(qū)現(xiàn)今地震活動密集區(qū)域的變化，也是區(qū)域應力場調(diào)整下的一種間接表現(xiàn)。

從唐山MS7.8地震序列時序變化來看（圖8），地震頻度隨時間呈指數(shù)衰減，地震強度逐漸減弱但存在起伏變化。在序列早期，余震衰減較快，2年內(nèi)最大余震衰減至ML5.0左右，85%左右的余震活動集中在該時間段內(nèi)。其后余震年頻次減少，地震強度較為穩(wěn)定，直至1991年5月30日發(fā)生唐山ML5.6地震，地震強度出現(xiàn)起伏。此后老震區(qū)ML4.0以上地震年頻次進一步減少，地震活動呈現(xiàn)階段性變化，與余震衰減趨勢相悖，可能意味著唐山老震區(qū)現(xiàn)今地震活動受1976年唐山MS7.8地震序列影響較弱。通常定義余震活動持續(xù)時間為余震活動減弱到與背景地震活動水平相當?shù)臅r間（蔣海昆等，2006），1992年以來區(qū)域ML4.0以上地震平均年頻次為1.2次，與震前區(qū)域地震平均年頻次（1.3次）基本相當（圖8b）。推測唐山老震區(qū)現(xiàn)今地震活動趨近于震前活動水平，發(fā)生的地震可能為背景地震。

4.2地震活動趨勢分析

唐山老震區(qū)現(xiàn)今中等地震活動存在3個階段，目前處于第3階段早期。從累積應變能量釋放速率、地震平均年頻次和最大地震震級的計算及統(tǒng)計結果來看，第2階段均低于第1階段，地震活動由活躍轉(zhuǎn)為平靜，符合老震區(qū)地震衰減規(guī)律。小震密集分布顯示，第2階段之后（即2004年之后），唐山老震區(qū)密集活動范圍收縮。因此，大致可將2004年看作該區(qū)域地震活動由活躍轉(zhuǎn)為平靜的一個時間節(jié)點。2004年以來唐山老震區(qū)中等地震活動形態(tài)呈現(xiàn)出“平靜—叢集”的周期性變化。2019年以后，累積應變能量釋放曲線顯示區(qū)域構造應力積累加速（圖5b），存在地震活動水平由平靜向活躍轉(zhuǎn)變的可能。結合小震密集分布特征，唐山斷裂帶北段以及灤縣—樂亭斷裂和盧龍斷裂的交匯區(qū)域可能是未來中強地震活躍的場所。

唐山老震區(qū)3個分區(qū)的地震活動差異顯著，并且存在中等地震東、中、西區(qū)交替發(fā)震的規(guī)律。依據(jù)對中等地震遷移過程的統(tǒng)計，遷移過程大多從中區(qū)開始，即構造應力多數(shù)從中區(qū)積累釋放。并且中區(qū)也是小震活動最為密集的區(qū)域，因此可以將中區(qū)作為區(qū)域地震活動的窗口，中區(qū)開始活動，意味著唐山老震區(qū)可能開啟活動。

本文對唐山老震區(qū)現(xiàn)今中等地震遷移現(xiàn)象的分析，結果可以嘗試應用于區(qū)域地震遷移趨勢分析，對未來中強地震發(fā)震位置的判定提供一定參考。根據(jù)遷移模式分析，認為2021年4月16日灤州ML4.8地震是一重遷移過程的結束地震，但考慮到以往震例中雙震情況的出現(xiàn)，認為東區(qū)還有可能發(fā)生一次ML4.0以上地震。2022年5月21日盧龍ML4.1地震的發(fā)生也驗證了這一預測結論。

5結論

本文根據(jù)唐山老震區(qū)現(xiàn)今地震的最小完整性震級，利用1992—2021年的小震目錄，對該區(qū)域現(xiàn)今地震活動特征進行了分析，得到以下結論：

（1）1992年以來唐山老震區(qū)小震密集分布于唐山斷裂帶北段以及灤縣—樂亭斷裂和盧龍斷裂的交匯區(qū)域，與1976年唐山MS7.8地震早期余震密集區(qū)域存在一定差異，可能是區(qū)域地震活動受強震序列影響減弱的表現(xiàn)。唐山MS7.8地震前、后老震區(qū)中等地震頻次統(tǒng)計結果顯示，現(xiàn)今中等地震年頻次與震前基本相當，區(qū)域地震活動趨近于震前活動水平，由此認為該區(qū)域所發(fā)生地震逐步趨近于背景地震。

（2）唐山老震區(qū)現(xiàn)今地震分布3個分區(qū)的差異顯著，主要集中于中區(qū)，東區(qū)次之，西區(qū)最弱。中等地震活動存在較為明顯的遷移規(guī)律，即地震在東、中、西區(qū)交替發(fā)震。將這一規(guī)律應用于區(qū)域地震遷移趨勢預測，可以對未來中強地震發(fā)震位置的判定提供一定的參考。

（3）唐山老震區(qū)現(xiàn)今地震活動呈現(xiàn)階段性變化，2019年以來，區(qū)域構造應力積累加速，地震活動水平存在由平靜向活躍轉(zhuǎn)變的可能?，F(xiàn)今小震活動密集分布于唐山斷裂帶北段以及灤縣—樂亭斷裂和盧龍斷裂的交匯區(qū)域，上述地區(qū)可能成為未來中強地震活躍的場所。[HTK]

文中圖件使用GMT軟件繪制，在此表示感謝。

參考文獻：

陳長云.2016.張家口—渤海斷裂帶分段運動變形特征分析［J］.地震，36（1）：1-11.Chen C Y.2016.Characteristics of segmentary motion and deformation along the Zhangjiakou-Bohai fault［J］.Earthquake，36（1）：1-11.（in Chinese）

郭蕾，孫麗娜，周依，等.2021.唐山地區(qū)近十年地震的重定位研究［J］.地震學報，43（6）：706-715.Guo L，Sun L N，Zhou Y，et al.2021.Relocation of earthquakes in Tangshan area in the recent decade［J］.Acta Seismologica Sinica，43（6）：706-715.（in Chinese）

國家地震局《一九七六年唐山地震》編輯組.1982.一九七六年唐山地震［M］.北京：地震出版社.Editorial Group of The Tangshan Earthquake of 1976，State Seismological Bureau.1982.Tangshan Earthquake of 1976［M］.Beijing：Seismological Press.（in Chinese）

虢順民，李志義，程紹平，等.1977.唐山地震區(qū)域構造背景和發(fā)震模式的討論［J］.地質(zhì)科學，12（4）：305-321.Guo S M，Li Z Y，Cheng S P，et al.1977.Discussion on the regional structural background and the seismogenic model of the Tangshan earthquake［J］.Scientia Geologica Sinica，12（4）：305-321.（in Chinese）

蔣長勝，吳忠良，莊建倉.2013.地震的“序列歸屬”問題與ETAS模型——以唐山序列為例［J］.地球物理學報，56（9）：2971-2981.Jiang C S，Wu Z L，Zhuang J C.2013.ETAS model applied to the Earthquake-Sequence Association（ESA）problem：the Tangshan sequence［J］.Chinese Journal of Geophysics，56（9）：2971-2981.（in Chinese）

蔣海昆，曲延軍，李永莉，等.2006.中國大陸中強地震余震序列的部分統(tǒng)計特征［J］.地球物理學報，49（4）：1110-1117.Jiang H K，Qu Y J，Li Y L，et al.2006.Some statistic features of aftershock sequences in Chinese mainland［J］.Chinese Journal of Geophysics，49（4）：1110-1117.（in Chinese）

李赫，董一兵，王熠熙，等.2020.唐山地區(qū)地震成因初探——以豐南2010年M4.1和2019年M4.5地震序列為例［J］.地球物理學進展，35（4）：1276-1283.Li H，Dong Y B，Wang Y X，et al.2020.Discussion on the seismogenesis of Tangshan area：take the sequence Fengnan M4.1 in 2010 and M4.5 in 2019 as an example［J］.Progress in Geophysics，35（4）：1276-1283.（in Chinese）

李瑞莎，崔效鋒，刁桂苓，等.2008.華北北部地區(qū)現(xiàn)今應力場時空變化特征研究［J］.地震學報，30（6）：570-580.Li R S，Cui X F，Diao G L，et al.2008.Temporal and spatial variation of the present crustal stress in northern part of North China［J］.Acta Seismologica Sinica，30（6）：570-580.（in Chinese）

劉保金，曲國勝，孫銘心，等.2011.唐山地震區(qū)地殼結構和構造：深地震反射剖面結果［J］.地震地質(zhì)，33（4）：901-912.Liu B J，Qu G S，Sun M X，et al.2011.Crustal structures and tectonics of Tangshan earthquake area：results from deep seismic reflection profiling［J］.Seismology and Geology，33（4）：901-912.（in Chinese）

劉蒲雄，呂曉健.2011.唐山地震序列空間分布復雜性原因研究［J］.地震，3（4）：1-14.Liu P，Lyu X J.2011.On the causes of the spatial distribution complexity of the 1976 Tangshan earthquake sequence［J］.Earthquake，31（4）：1-14.（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］.Seismology and Geology，2（2）：43-54.（in Chinese）

孫麗娜，齊玉妍，陳婷，等.2021.震源機制和場地條件在河北地震影響場判定中的應用研究［J］.地震學報，43（4）：508-520.Sun L N，Qi Y Y，Chen T，et al.2021.Research on application of focal mechanism and site conditions in judgment of Hebei earthquake influence field［J］.Acta Seismologica Sinica，43（4）：508-520.（in Chinese）

孫啟凱，池國民，徐東卓.2017.首都圈地區(qū)地殼垂直形變特征及剖面分析［J］.大地測量與地球動力學，37（5）：497-501.Sun Q K，Chi G M，Xu D Z.2017.Characteristics of vertical deformation and profile analysis for the capital circle［J］.Journal of Geodesy and Geodynamics，37（5）：497-501.（in Chinese）

王輝，劉勉，李振，等.2013.華北地區(qū)大地震的長余震序列［J］.地震，33（3）：1-12.Wang H，Liu M，Li Z，et al.2013.Long aftershock sequences of major earthquakes in North China［J］.Earthquake，33（3）：1-12.（in Chinese）

王想，馮向東，宮猛，等.2012.唐山老震區(qū)地震活動特征［J］.地震地磁觀測與研究，33（S2）：44-48.Wang X，F(xiàn)eng X D，Gong M，et al.2012.Study on seismicity features of Tangshan seismic region［J］.Seismological and Geomagnetic Observation and Research，33（S2）：44-48.（in Chinese）

王曉山，馮向東，趙英萍.2020.京津冀地區(qū)地殼應力場特征［J］.地震研究，43（4）：610-619.Wang X S，F(xiàn)eng X D，Zhao Y P.2020.Characteristics of crustal stress field in Beijing-Tianjin-Hebei region［J］.Journal of Seismological Research，43（4）：610-619.（in Chinese）

聞學澤，馬勝利.2006.唐山大地震對相鄰斷裂段地震復發(fā)的影響［J］.自然科學進展，16（10）：1346-1350.Wen X Z，Ma S L.2006.Effect of strong Tangshan earthquake on seismic recurrence in adjacent faults［J］.Progress in Natural Science，16（10）：1346-1350.（in Chinese）

徐杰，周本剛，計鳳桔，等.2012.華北渤海灣盆地區(qū)大震發(fā)震構造的基本特征［J］.地震地質(zhì)，34（4）：618-636.Xu J，Zhou B G，Ji F J，et al.2012.Features of seismogenic structures of great earthquakes in the Bohai Bay Basin area，North China［J］.Seismology and Geology，34（4）：618-636.（in Chinese）

許忠淮，閻明，趙仲和，等.1983.由多個小地震推斷的華北地區(qū)構造應力場的方向［J］.地震學報，5（3）：268-279.Xu Z H，Yan M，Zhao Z H，et al.1983.Evaluation of the direction of tectonic stress in North China from recorded data of a large number of small earthquake［J］.Acta Seismologica Sinica，5（3）：268-279.（in Chinese）

許忠淮.1985.用滑動方向擬合法反演唐山余震區(qū)的平均應力場［J］.地震學報，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］.Acta Seismologica Sinica，7（4）：349-362.（in Chinese）

楊理華，趙喜柱.1984.1976年唐山地震為三主震型地震［J］.地震研究，7（6）：639-648.Yang L H，Zhao X Z.1984.The Tangshan earthquake of 1976-an event of triple main shock pattern［J］.Journal of Seismological Research，7（6）：639-648.（in Chinese）

張宏志，刁桂苓，陳祺福，等.2008.1976年唐山7.8級地震震區(qū)現(xiàn)今地震震源機制分析［J］.地震研究，31（1）：1-6.Zhang H Z，Diao G L，Chen Q F，et al.2008.Focal mechanism analysis of the recent earthquakes in Tangshan seismic region of M7.8 in 1976［J］.Journal of Seismological Research，31（1）：1-6.（in Chinese）

張素欣，孫麗娜，李永慶，等.2020.2019年12月5日河北唐山M4.5地震的發(fā)震背景［J］.華北地震科學，38（S1）：76-79.Zhang S X，Sun L N，Li Y Q，et al.2020.The seismogenic background of the M4.5 earthquake in Hebei Tangshan on December 5，2019［J］.North China Earthquake Sciences，38（S1）：76-79.（in Chinese）

張素欣，王曉山，陳婷，等.2017.唐山老震區(qū)40年地震時空演化特征分析［J］.華北地震科學，35（1）：32-37.Zhang S X，Wang X S，Chen T，et al.2017.Seismic space-time evolution characteristics in Tangshan old seismic region during the past 40 years［J］.North China Earthquake Sciences，35（1）：32-37.（in Chinese）

張躍剛，王玉珍，尹寶軍，等.2013.2012年5月28日河北省唐山4.8級地震［J］.中國地震，29（2）：219-229.Zhang Y G，Wang Y Z，Yin B J，et al.2013.The Tangshan MS4.8 earthquake on May 28，2012 in Hebei Province［J］.Earthquake Research in China，29（2）：219-229.（in Chinese）

鄭建常，李東梅，王鵬，等.2014.基于最鄰近事件距離算法的叢集地震識別方法及應用［J］.地震，34（4）：100-109.Zhang J C，Li D M，Wang P，et al.2014.Method and application of clustering seismicity identification based on nearest-neighbor distance［J］.Earthquake，34（4）：100-109.（in Chinese）

鄭建常，李冬梅，冀東普，等.2022.日本與華北地區(qū)強震活動相關性再討論［J］.華南地震，42（4）：7-19.Zheng J C，Li D M，Ji D P，et al.2022.Rediscussion on the correlation of strong earthquake activities between Japan and North China［J］.South China Journal of Seismology，42（4）：7-19.（in Chinese）仲秋，史保平.2012.1976年唐山MS7.8地震余震序列持續(xù)時間及對地震危險性分析的意義［J］.地震學報，34（4）：494-508.Zhong Q，Shi B P.2012.Aftershock duration of the 1976 MS7.8 Tangshan earthquake and implication for seismic hazard estimation［J］.Acta Seismologica Sinica，34（4）：494-508.（in Chinese）

朱琳，李騰飛，石富強，等.2021.1976年唐山強震群震后庫侖應力演化及其與2020年古冶5.1級地震的關系［J］.地震研究，44（1）：1-8.Zhu L，Li T F，Shi F Q，et al.2021.The coulomb stress evolution after the 1976 Tangshan strong earthquake swarm and its relationship with the 2020 Guye M5.2 earthquake［J］.Journal of Seismological Research，44（1）：1-8.（in Chinese）

Butler R，Stewart G S，Kanamori H.1979.The July 27，1976 Tangshan，China，earthquake：a complex sequence of intraplate events［J］.Bulletin of the Seismological Society of America，69（1）：207-220.

Li L，Chen Q F，Cheng X，et al.2007.Spatial clustering and repeating of seismic events observed along the 1976 Tangshan fault，north China［J］.Geophysical Research Letters，34（23）：L23309.

Robinson R，Zhou S.2005.Stress interactions within the Tangshan，China，earthquake sequence of 1976［J］.Bulletin of the Seismological Society of America，95（6）：2501-2505.

Shedlock K M，Baranowski J，Xiao W W，et al.1987.The Tangshan aftershock sequence［J］.Journal of Geophysical Research，92（B3）：2791-2803.

Stein S，Liu M.2009.Long aftershock sequences within continents and implications for earthquake hazard assessment［J］.Nature，462（5）：87-89.

Characteristics of the Recent Seismicity in the Tangshan Old Seismic Region

ZHOU Yi1，2，CAO Jun2，MA Xudong2，SHAN Lianjun2

（1.Hebei Hongshan National Observatory on Thick Sediments and Seismic Hazards，Xingtai 054000，Hebei，China）?（2.Hebei Earthquake Agency，Shijiazhuang 050021，Hebei，China）

Abstract

To analyze the seismic background in the Tangshan old seismic region and provide earthquake trend predictions，we describe the characteristics of the recent seismicity in this region.Results show that the recent seismicity in the Tangshan old seismic region is different from the early aftershock intensive area of the 1976 Tangshan MS7.8 earthquake.Considering the statistical results of regional background seismicity，the influence of the Tangshan earthquake is gradual weakening，and the recent seismicity in the Tangshan old seismic region can be considered as background seismicity.The characteristics of seismicity are different in the three sub-areas of the Tangshan old seismic region，the earthquakes mainly concentrate in the central area，followed by the eastern area，and the western area was the weakest，and the moderate earthquake following a law of alternate occurrence in the east，central and west sub-areas.This seismic migration law constitutes a reference for the determination of the locations of future moderate-to-strong earthquakes.At present，the accumulation of tectonic stress is accelerating in the old seismic region of Tangshan.Seismicity changing from quiet to active is possible，and moderate-to-strong earthquakes may occur at the northern segment of the Tangshan Fault and the intersection of the Luanxian-Laoting Fault and the Lulong Fault in the future.

Keywords：?the 1976 Tangshan MS7.8 earthquake；characteristics of seismicity；background seismicity；seismic migration law

收稿日期：2022-10-27.

基金項目：河北省重點研發(fā)項目（20375404D）；河北省地震科技星火計劃項目（DZ2021121700014、DZ2024112100001）；河北省地震局創(chuàng)新團隊項目（DZ202109080114）.

第一作者簡介：周依（1991-），工程師，主要從事數(shù)字地震學和地震活動性研究.E-mail：zhouyimonday@163.com.

周依，曹筠，馬旭東，等.2024.唐山老震區(qū)現(xiàn)今地震活動特征［J］.地震研究，47（3）：430-440，doi：10.20015/j.cnki.ISSN1000-0666.2024.0046.

Zhou Y， Cao J， Ma X D，et al.2024.Characteristics of the recent seismicity in the Tangshan old seismic region［J］.Journal of Seismological Research，47（3）：430-440，doi：10.20015/j.cnki.ISSN1000-0666.2024.0046.