中國大陸活動地塊邊界帶主要斷層的強震震間晚期綜合判定

邵志剛，武艷強，季靈運，刁法啟，石富強，李玉江，龍鋒，張輝，朱良玉，陳長云，王武星，魏文薪 ，王芃，劉曉霞，劉琦，潘正洋，尹曉菲，劉月，馮蔚，鄒鎮宇，曹建玲， 徐晶，韓立波，程佳，魯人齊0，徐岳仁，李西，孫鑫喆

1 中國地震局地震預測研究所，地震預測重點實驗室，北京 100036 2 中國地震局第一監測中心，天津 300180 3 中國地震局第二監測中心，西安 710054 4 中國地質大學（武漢）地球物理與空間信息學院，武漢 430074 5 陜西省地震局，西安 710068 6 國家自然災害防治研究院，北京 100085 7 四川省地震局，成都 610044 8 甘肅省地震局，蘭州 730013 9 中國地震局地球物理研究所，北京 100081 10 中國地震局地質研究所，北京 100029 11 云南大學地球科學學院，昆明 650500

0 引言

強震因強烈破壞性和災害的不可逆轉性，成為影響社會發展和百姓福祉的嚴重威脅，因此強震預報有著強烈的社會需求.當然，對于地震綜合預測結論必須給出確定的“時、空、強”三要素，而且最好能用定量的方式表述預測結果（陳運泰，2009）.對于這三個要素，目前給出準確的預測結果仍比較困難，但是否可以針對單個要素適當放寬要求并通過不斷改進和完善，以期達到提高地震預報準確性的效果.

其中對于時間預測來講，針對強震孕育發生過程不同階段，時間可分為：同震、震后、震間和震前四個時間段，而震間又分為：早期、中期和晚期（力武常次，1978；Scholz,1998）.對于目標斷層來講，震前的時間尺度、物理機制、有效監測等方面仍存在諸多問題.比較現實的是，進一步放寬對時間精度的要求，能否利用現有技術方法來確定目標斷層是否處于強震孕育晚期.

近20年來基于較密集觀測資料的積累、理論計算方法的發展、實際地震活動的檢驗，在很多強震前地震地質、大地測量、地震活動、數值模擬等不同學科均有研究方法用于強震長期危險性分析，而且取得較好的長期預測效果.其中有些用于強震孕育晚期判定的有效研究方法，對2010年智利8.8級、2011年日本9.0級、2014年智利8.1級、2015年尼泊爾8.1級等地震做出了較好的強震背景預測（Madariaga et al.,2010；Hashimoto et al.,2009;Loveless and Meade,2010；Schurr et al.,2014；Avouac et al.,2015）；例如，2010年智利8.8級地震發生后就有學者在Science上發表文章稱，早在十幾年前就知道要發生大地震（Madariaga et al.,2010）.因此，2011年日本9級地震后，著名地震學家金森博雄呼吁更好地利用已有技術手段來改進地震預報研究（Kanamori,2012）.

可用于震間晚期判定的方法，概括起來最具代表性的有，地震地質的強震破裂空段、大地測量的斷層運動閉鎖段、地震活動的中小地震稀疏段、數值模擬的庫侖應力增強顯著段等，從震例角度來講，這些方法在很多強震之前都有很好的研究結果.但對于震間晚期來講，這些方法給出的強震孕育晚期大多數判定依據是必要不充分條件，即斷層處于強震孕育晚期時多呈現為強震破裂空段、大地測量的斷層運動閉鎖段、數值模擬的庫侖應力增強顯著段，部分出現地震活動的中小地震稀疏段，但出現上述特征時，目標斷層并不一定處于強震孕育晚期.因此，本文嘗試根據以往震例表現，采用上述方法綜合給出強震孕育晚期判定.

中國大陸是典型的大陸動力學環境，“垂向分層，橫向分塊”是大陸活動地塊發生相對運動和變形是“陸內變形”的重要方式之一.而活動地塊是指被形成于晚新生代、晚第四紀（10—12萬年）至今仍強烈活動的構造所分割和圍限、具有相對統一運動方式的地質單元（張培震等，2003），中國大陸地區共有6個I級活動地塊，22個II級活動地塊；而介于活動地塊之間的便是活動地塊邊界帶，目前確認的活動地塊邊界帶有26個（圖1）.中國大陸主要構造變形和淺源強震大都發生在邊界帶上，有歷史記載以來幾乎全部8 級和超過80%的7 級以上強震發生在活動地塊邊界帶上（張國民等，2005；張培震等，2013；圖1）.

所以，本文以中國大陸活動地塊邊界帶為研究對象，根據活動斷層研究結果，重點關注中國大陸地區的391個斷層段（圖2），嘗試綜合應用強震破裂空段、斷層運動閉鎖段、庫侖應力顯著增強段和中小地震稀疏段等方法來判定處于強震孕育晚期的斷層段，期望能給長期地點預測提供可參考的科學依據.

1 強震破裂空段

強震破裂空段是指板塊邊界帶或活動塊體邊界帶等大型構造帶內，長時間未發生強震破裂的斷層段（Sykes,1971），這些段落相對于發生強震破裂不久的斷層段所積累的斷層應力應變更高，因而，強震破裂空段成為強震預測的重要方法（M7專項工作組，2012）.

在上述強震空間和時間上并非完全隨機的認識下，從長期預測角度來看，強震破裂空段取得較好的地點預測效果，成為高震級地震地點預測的重要地震地質學標志（徐錫偉等，2017）.例如，2004年蘇門答臘9.2級地震（Kagan and Jackson,1991）、2014年智利8.1級地震（Kagan and Jackson,1991;Nishenko,1991）、2010年智利8.8級地震（Nishenko,1991;Madariaga et al.,2010）、2015年尼泊爾8.1級地震（Bilham et al.,2001；Bollinger et al.,2014）等均發生在強震破裂空段.在中國大陸長期重點危險區判定中同樣發揮重要作用，例如，2008年汶川8.0級地震、2010年玉樹7.1級地震、2013年蘆山7.0級地震等均發生在強震破裂空段（M7專項工作組，2012；《2006～2020年中國大陸地震危險區與地震災害損失預測研究》項目組，2007）.

圖1 中國大陸地區活動地塊與歷史強震空間分布圖圖中名稱為二級活動地塊名稱，編號為活動地塊邊界帶編號（張培震等，2003；張國民等，2004）.Fig.1 Active tectonic blocks in Chinese mainland and its spatial distribution of historical strong earthquakesThe font is the name of secondary active blocks,the number represents the boundary zone of active tectonic block （Zhang et al.,2003;Zhang et al.,2004）.

圖2 中國大陸地區活動地塊邊界帶及其主要斷層空間分布圖（張培震等，2003；張國民等，2004）Fig.2 The distribution of the boundary zone of active tectonic block and main faults in Chinese mainland （Zhang et al.,2003;Zhang et al.,2004）

對于活動構造帶內強震破裂空段的工作基礎是古地震、歷史地震和現代地震的同震破裂信息，具體包括強震發生時間、同震破裂空間展布、同震破裂平均位錯矢量或破裂過程等.對于保存于地質記錄中的史前和歷史沒有明確記載的古地震事件，一般通過強震破裂地表遺跡探測技術、古地震年代測定技術等確定古地震同震破裂信息（冉勇康和鄧起東，1999；冉勇康等，2014）.對于歷史地震，還有地震災害的歷史文字記載可以用于判定極震區、地表破裂等空間展布（Wen et al.,2008）；而對于現代地震，可以用于同震破裂過程或同震位錯反演的資料越來越多，包括大地測量儀器、地震儀、強震儀等各種資料.依據古地震、歷史地震和現代地震的同震破裂，結合活動構造相關研究，確定具體斷層段的斷層分段、滑動速率、強震離逝時間、復發周期等（鄧起東和聞學澤，2008），這些地震地質研究結果是長期地震預測的重要基礎，是判定強震破裂空段的主要依據.

而實際判定過程中強震破裂空段分為兩類：（1）高離逝率空段（離逝率等于上次強震距離現在的時間與復發周期的比值），上次強震離逝時間接近或超過復發周期的斷層段；（2）歷史缺震空段，沒有離逝率相關研究結果，且沒有歷史強震記錄的斷層段（徐錫偉等，2017）.本文收集整理了已有的研究結果，主要包括“八五”期間中國大陸中長期強震預測研究結果（時振梁等，1997）、第二版中國大陸地震重點監視防御區研究結果（《2006～2020年中國大陸地震危險區與地震災害損失預測研究》項目組，2007）、M7工作研究結果（M7專項工作組，2012）、活動斷層探察等地震地質研究結果（徐錫偉等，2016，2017），中國大陸目前能識別出的強震破裂空段共有66個，其中高離逝率空段29個，缺震空段37個（圖3）；當然，現有結果的完備程度需要更全面的活動斷層探察工作，研究工作需要更準確的測試技術.

圖3 中國大陸活動地塊邊界帶主要斷層與強震破裂空段分布圖（張培震等，2003；張國民等，2004；徐錫偉等，2017）Fig.3 The distribution of main active faults and some strong seismic rupture gaps in the boundary zone of active tectonic blocks in Chinese mainland （Zhang et al.,2003;Zhang et al.,2004;Xu et al.,2017）

近年來隨著地貌觀測技術和測年技術的發展，古地震和歷史地震的同震破裂的相關信息越來越準確，強震破裂空段的判定會越來越可靠、越來越完備.但是對于地震地質的研究方法來講，強震破裂空段除了測年的精度、破裂空間展布等問題外，還有兩個難題，一是沒有地表破裂或地表破裂不明顯的古地震事件的判定，另外一個是復雜構造地區強震空區的確定.

另外，強震離逝率是定量地描述特定斷層是否處于強震孕育晚期的有效參數，離逝率越大的強震破裂空段時間上越逼近下次強震，一般認為接近1的強震破裂空段存在強震危險性，但實際上有些強震破裂空段的離逝率超過，甚至達到4，仍沒發生地震（徐錫偉等，2017）.由中國大陸強震原地復發時離逝率的概率密度函數可知（聞學澤，1995），在離逝率為1左右的強震發生個數比例最高，但在離逝率為0.5～1.5范圍之間，也即離逝時間達到復發周期的一半后，強震的危險性會越來越高（圖4b）.

圖4 中國大陸強震發生時離逝率的概率密度函數（a）及積累概率函數（b）藍線和紅線分別代表對數正態函數與布朗時間過程函數.Fig.4 Probability density function （a） and accumulation probability function （b） of lapse rate during strong earthquakes in Chinese mainlandBlue line and red line represent logarithm-normal function and Brownian time process function respectively.

2 斷層運動閉鎖段

彈性回跳理論的基礎觀測資料是1906年舊金山大地震震前、同震三角測量資料給出的地震前后地殼變形.該理論認為，由于震前圣安德烈斯斷層兩側地塊存在相對運動，而斷層處于閉鎖狀態，導致了斷層附近地殼中彈性應變能的積累，當應變積累到臨界點后斷層回跳發生錯動，產生地震波動（Reid,1910）.此后，隨著大地測量觀測儀器的不斷改進和震例的不斷豐富，逐漸認識到地震過程具體應分為震間（孕育）、同震（發生）和震后（后效）（力武常次，1978）.而不同孕震階段跨斷層地殼形變具有不同的表現形式（圖5；Savage and Burford，1973；Savage and Prescott，1973；Meade and Hager,2005）；而且斷層上地震應變的積累以及震后變形調整不僅與上地殼脆性層有關，還與下地殼或上地幔深部韌性層長期穩態相對運動有關（Nur and Mavko,1974；Deng et al.,1998）；在大量研究基礎上，Scholz（1998）給出了斷層運動在強震孕育發生不同階段隨深度變化的理論模式.處于震間不同階段的斷層運動與地表變形速率空間分布特征之間的理論關系，是利用大地測量觀測資料開展孕震階段判定的理論基礎.

圖5 斷層運動不同階段地殼地表變形速率空間分布示意圖（Meade and Hager,2005）（a） 斷層處于同震（紅色）、震間期（橙色）地表變形速率空間分布及和斷層兩側地塊平均相對運動速率（黑色）空間分布示意圖；（b） 斷層處于強震孕育階段不同時期的地表變形速率空間分布示意圖，黃色為孕震早期、橙色為孕震中期、紅色為孕震晚期.圖中v是地表變形速度,v0是斷層長期滑動速率，x是點位與斷層的距離，D是斷層閉鎖深度.Fig.5 The map of spatial distribution of crustal surface deformation rate at different stages of fault movement （Meade and Hager,2005）（a） Spatial distribution diagram of surface deformation rate during coseismic （red） and interseismic period （orange） of the fault and the spatial distribution of average relative movement rate （black） of blocks on both sides of the fault;（b） Spatial distribution diagram of surface deformation rate at different stages of the fault in the strong earthquake preparation stage,yellow is the early seismogenic stage,orange is the epicenter stage,and red is the late seismogenic stage.In the figure,v is the surface deformation velocity,v0 is the long-term slip rate of the fault,x is the distance between the point and the fault,and D is the fault locking depth.

而基于斷層運動特征分析目標斷層所處的孕震階段的技術基礎，則是負位錯、位錯解析解、位錯數值模擬等.Matsu′ura等（1986）提出負位錯理論模型，該模型認為受震間構造應力的加載作用引起的發生在深部“延性”區域相對運動，可認為是無應變積累的剛性塊體運動加上淺部“脆性”閉鎖層的反方向均勻位錯；與此同時，關于斷層運動引起地殼變形的計算理論逐漸成熟，尤其是Okada的解析解（Okada,1985,1992）為相關地震變形的模擬和反演提供非常堅實的理論基礎.

近年來利用大地測量給出的地殼變形速率為約束，基于Okada理論計算和負位錯理論模型，通過斷層運動反演確定斷層震間強閉鎖區，進而判定目標斷層是否處于震間晚期階段，累積了越來越多的震例，例如，2010年智利8.8級、2011年日本9.0級、2014年智利8.1級、2015年尼泊爾8.1級等地震（Hashimoto et al.,2009;Madariaga et al.,2010；Loveless and Meade,2010；Schurr et al.,2014；Avouac et al.,2015）.其中對于觀測比較好的地方，應用比較多的軟件是DEFNODE（McCaffrey et al.,2000；McCaffrey,2002），該軟件基于Okada模型，以地表變形觀測為約束，反演斷層運動速率，結合地震地質研究給出的活動斷層最新地質年代的斷層平均運動速率，給出斷層運動閉鎖程度，一般認為處于震間晚期的斷層運動閉鎖程度較強.而對于觀測基礎不太好的地區，基于負位錯理論利用反正切的方法反演斷層閉鎖深度，一般認為震間晚期的斷層運動閉鎖深度較深.

而依據現有中國大陸大地測量觀測資料開展斷層孕震階段的判定，對中國大陸活動地塊邊界帶主要斷層的391個斷層段，本文將閉鎖程度結果分為4類，即強閉鎖、部分閉鎖、不閉鎖及無資料.強閉鎖是指斷層段的閉鎖深度接近該區域孕震層深度，且閉鎖率不低于0.8，例如，川滇地區小江斷裂帶南段、安寧河斷裂帶等.部分閉鎖是指斷層段的閉鎖深度小于該區域孕震層深度，但大于10 km，且閉鎖率不低于0.8，例如，青藏高原北部地區阿爾金斷裂（Li et al.,2018）；或者斷層段的閉鎖深度接近該區域孕震層深度，閉鎖率在0.6至0.8之間，例如，川滇地區鮮水河斷裂帶中段等（Jiang et al.,2015）.不閉鎖是指斷層段的閉鎖深度小于5 km，例如，川滇地區龍門山斷裂帶北段、小江斷裂帶中段等.無資料是指斷層段附近區域的大地測量資料不充分，無法得到可靠的反演結果，例如，川滇地區瀾滄江斷裂、新疆地區興地斷裂等.對于無法開展二維反演的部分地區采用反正切的方法確定閉鎖深度，如果閉鎖深度達到或接近彈性層上地殼厚度則認為是強閉鎖，而閉鎖深度較淺的認為不閉鎖，例如，海原斷裂帶閉鎖深度僅1 km左右（Cavalié et al.,2008），因此被認為有可能仍處于震后調整階段（M7專項工作組，2012）.

利用現有大地測量觀測資料為約束，反演中國大陸活動地塊邊界帶主要斷層391個斷層段的斷層運動閉鎖程度，目前能識別出強閉鎖斷層段103個，部分閉鎖斷層段58個，不閉鎖斷層段64個，無資料斷層段166個（圖6）.

中國大陸活動地塊邊界帶的391個斷層段，目前大約有40%無法開展斷層運動反演工作.而監測能力的提升，并非單靠簡單的空間加密就能改變的，因為地表變形觀測對深度約束能力畢竟有限，例如，對于監測能力較好的美國圣安德列斯斷層中段，斷層運動速率反演空間分辨率由淺層的2～5 km，到20 km深度會達到20 km左右（Jolivet et al.,2015）.我們通過理論反演實驗，來討論最優觀測方案，針對不同閉鎖深度的斷層運動模型，通過理論計算地表運動速率，對不同閉鎖深度模型對應的跨斷層地表運動速率剖面，重點分析刻畫地殼運動速率梯度所需的測點間距；然后以理論計算地表運動速率為約束，通過反演分析不同深度空間分辨率.理論分析認為，每個斷層分段布設跨斷層GPS觀測剖面，斷層附近50 km范圍內臺間距2～5 km，超過50 km臺間距20 km，其單個剖面斷層運動一維反演空間分辨率和精度如圖7；剖面間沿斷層兩側布設約20 km間距，經過理論測試該綜合加密觀測方案的空間分辨率應該能滿足識別出7級地震凹凸體的監測需求.

3 斷層庫侖應力顯著增強段

斷層應力累積水平計算的常用方法是斷層面上庫侖應力變化數值模擬，而很多震例表明，強震引起的庫侖應力變化可以用于強震地點的長期預測.尤其是1999年土耳其伊茲米特7.4級地震之后該研究引起地震學界普遍關注，在此地震前Nalbant等（1998）和Stein等（1997）利用歷史強震引起的同震變化分析認為，土耳其北安納托利亞斷層存在強震危險.震后不同物理機制也在地震活動之間庫侖應力的影響中發揮了關鍵作用，例如，1992年美國蘭德斯7.3級地震對1999年赫可托礦地震庫侖應力影響中震后黏滯松弛作用起到關鍵作用（Freed and Lin,2001）.而震后黏滯松弛作用可能會持續數十年甚至上百年（Lorenzo-Martín et al.,2006），有些震例表明周邊強震引起的震后黏滯松弛作用在某些時段有可能達到構造加載作用的水平（徐晶等，2013），例如，2008年汶川8.0級地震震后黏滯松弛在鮮水河斷裂帶中段引起的震后庫侖應力變化隨著時間的累積量可能會達到甚至超過同震引起的庫侖應力變化（圖8）.而在震后黏滯松弛作用引起的庫侖應力變化實際計算過程中，區域殼幔黏滯模型的確定非常重要（Freed and Lin,2001），其主要方法有：（1）基于震后形變觀測做約束反演得到區域殼幔黏滯結構（Deng et al.,1998；Diao et al.,2014）；（2）基于地表運動速率和殼幔溫度約束反演得到區域殼幔黏滯結構（Copley,2008；石耀霖和曹建玲，2008）.

圖6 中國大陸活動地塊邊界帶主要斷層與各斷層閉鎖程度結果分布圖Fig.6 The distribution of main active faults and their locking degree in the boundary zone of active tectonic blocks in Chinese mainland

圖7 不同斷層閉鎖深度反演分析（a），垂向反演空間分辨率和精度（b）（a） 依據不同閉鎖深度斷層運動模型給出的地表測點數值模擬結果做約束，進行斷層運動反演測試，給出不同斷層閉鎖深度模型的空間分辨率（km）和反演精度（mm）；（b） 為基于15—20測點跨斷層GPS剖面觀測方案的一維垂向反演的空間分辨率和精度.Fig.7 Inversion of different fault locking depths （a）,Spatial resolution and accuracy of vertical inversion （b） （a） Inversion test of fault motion was conducted based on the numerical simulation results of surface measurement points given by fault motion models with different locked depths,and the spatial resolution （km） and inversion accuracy （mm） of models with different locked depths were given.（b） shows the spatial resolution and accuracy of one-dimensional vertical inversion based on GPS profile observation scheme across faults at 15—20 survey points.

圖8 2008年汶川8.0級地震在鮮水河斷裂帶色拉哈—康定段引起的庫侖應力變化（a） 2008年汶川8.0級地震引起的同震庫侖應力變化，接收斷層為色拉哈—康定段的斷層，深度10 km；（b） 汶川地震引起的同震和震后庫侖應力累積變化.Fig.8 Coulomb stress changes caused by the 2008 MS8.0 Wenchuan earthquake in the Salah-Kangding segment of the Xianshuihe fault zone（a） shows the coseismic Coulomb stress change caused by the 2008 Wenchuan MS8.0 earthquake.The geometry of the receiving fault is the fault geometry of the Salah-Kangding section,with a depth of 10 km.（b） shows the coseismic and post-seismic Coulomb stress accumulation changes caused by Wenchuan earthquake.

庫侖應力計算基礎資料包括歷史強震同震位錯信息、接收斷層信息和殼幔黏滯系數.因此，在進行庫侖應力計算前開展了大量基礎資料收集工作.一是整理了中國大陸及周邊179次歷史強震同震位錯，以及391個斷層段的地震地質研究結果（徐錫偉等，2016）.二是通過基于震源機制解的應力張量反演，確定缺乏明確地震地質研究結果的斷層段和活動地塊內部的接收斷層參數（Toda and Enescu,2011），共計得到中國大陸及周邊地區12424次地震事件的震源機制解.三是系統梳理了利用中國大陸強震震后位移、地殼變形速率和溫度模型等方法，給出的不同區域殼幔短期黏滯系數和長期黏滯系數（石耀霖和曹建玲，2008）.

基于收集的中國大陸及周邊179次強震震源參數和不同分區殼幔黏滯系數，利用德國地學研究中心汪榮江研究員開發的庫侖應力計算程序PSGRN/PSCMP（Wang et al.,2006），分區計算了區域歷史強震對活動地塊邊界帶主要斷層各斷層段和活動地塊內部的庫侖應力變化.

中國大陸活動地塊邊界帶主要斷層391個斷層段中，庫侖應力顯著增強（庫侖應力增強不少于0.1 MPa）的斷層段有107個，庫侖應力增強（庫侖應力增強處于0.01～0.1 MPa）的斷層段有84個（圖9）.

4 中小地震稀疏段

根據以往震例，1999年中國臺灣集集7.6級（Lin,2001）、2004年蘇門答臘9.0級（Katsumata,2015）、2011年日本9.0級（Katsumata,2011）等地震震前震源區存在中小地震活動較弱區域.針對中國大陸活動地塊邊界帶主要斷層的各斷層段，基于1970年以來中小地震目錄，以斷層節點為中心，統計該剖面內3級以上地震的數量，按照各斷層段所在邊界帶進行活動程度歸一化處理，若歸一化地震活動程度小于0.4，判定該斷層段為中小地震稀疏段.例如，川滇地區龍日壩斷裂帶中段、安寧河斷裂帶南段.

另外，1989年美國加州7.1級地震震前震源區中小地震呈現立體圍空現象（Olson,1990），僅通過地震事件經緯度可能識別不出中小地震稀疏段.因此，針對南北地震帶、天山地震帶和華北地區等重點區域，采取中小地震重新定位方法（Long et al.,2015），對區域剖面的定位結果進行分析，進而確定是否為中小地震稀疏段.例如，天橋溝—黃羊川斷裂、小金河斷裂南段.

利用上述方法，逐一對中國大陸活動地塊邊界帶主要斷層391個斷層段進行分析，共識別出中小地震稀疏段83個（圖10），需要說明的是，由于監測能力的原因，并未對青藏高原內部進行中小地震稀疏段的分析.

5 強震孕育晚期綜合判定與結果

從以往震例來看，強震破裂空段、斷層運動閉鎖段、斷層庫侖應力增強段是強震孕育晚期判定的必要不充分條件，而中小地震稀疏段的震例相對較少，基于國內外強震預測效果實際情況（表1），分別賦予強震破裂空段、斷層運動閉鎖段、中小地震稀疏段、斷層庫侖應力增強段等4個震間晚期判定方法3.0、3.0、1.0、3.0的權重，各類方法具體權重分配規則（圖11）：（1）高離逝率強震破裂空段為3.0，歷史缺震空段為2.0；（2）斷層運動強閉鎖段為3.0，部分閉鎖為1.5；（3）中小地震稀疏段為1.0；（4）庫侖應力增強≥0.1 MPa為3.0，（0.1～0.01） MPa為0.15，綜合加權給出目標斷層的震間晚期判定系數.這種策略與美國加州地區震級判定中間部分過程類似（Field et al.,2014），只不過本文更多地參考方法在實際震例中的印證情況.

圖12為中國大陸活動地塊邊界帶主要斷層各斷層段的震間晚期判定系數空間分布圖.一般情況，強震發生后，發震斷層段的離逝率低，處于震間晚期的可能性相對較小；中國大陸7級淺源地震（西部1900年以來，東部1484年以來）的發震斷層，約90%的震間晚期判定系數<5.0.因此，基于現有地震觀測資料分析認為，震間晚期判定系數≥5.0斷層處于震間晚期的可能性較大.

由圖12可見，中國大陸地區可能處于強震震間晚期的斷層段主要有：

中國大陸東部的郯廬斷裂帶的宿遷—嘉山段、萊州灣段、下遼河段、五常段，河北平原地震帶磁縣—鶴壁段，環鄂爾多斯的桌子山西緣斷裂、狼山山前斷裂、黃河斷裂臨武段、六盤山東麓斷裂、小關山斷裂、渭河斷裂、隴縣—岐山—馬昭斷裂、大青山山前斷裂、懷安鎮盆地北緣斷裂、陽原盆地南緣斷裂、交城斷裂、中條山山前斷裂，東南沿海的濱海斷裂東北段；

青藏高原北部地區的托萊山斷裂、冷龍嶺斷裂、日月山斷裂、天橋溝—黃羊川斷裂、長嶺山北麓斷裂、金強河斷裂、毛毛山斷裂、老虎山斷裂、香山天景山斷裂東南段，西秦嶺北緣斷裂的鍋麻灘段、黃香溝段、漳縣段與臨潭—宕昌斷裂，東昆侖斷裂帶西大灘—東大灘段、瑪沁—瑪曲段、瑪曲—荷葉段、塔藏段，阿爾金斷裂阿克塞段、疏勒河段、瑪依色列克—茫崖段、索爾庫里段，康西瓦斷裂；青藏高原中南部及川滇地區的喜馬拉雅主斷裂中西段、中東段，嘉黎斷裂，甘孜玉樹斷裂當江段，班公怒江縫合帶東段，龍日壩斷裂中段、北段，鮮水河斷裂中南段，龍門山斷裂帶南段，小金河—麗江斷裂中段，理塘斷裂大毛埡壩段、沙灣段，安寧河斷裂、大涼山斷裂、馬邊斷裂、昭通斷裂、小江斷裂巧家—東川段、小江斷裂建水段，南汀河斷裂南段、龍陵—瑞麗斷裂、孟連斷裂，楚雄—建水斷裂、紅河斷裂中南段、程海斷裂南段、元謀—昔格達斷裂羅川以南段；

圖11 4類方法的權重分配規則確定震間晚期判定系數Fig.11 The judgment coefficients of the late stage of interseismic period determined by the rules of the weight allocation in four kinds of methods

南天山地區的邁丹斷裂、柯坪斷裂、霍拉山斷裂、庫木格列木斷裂，北天山地區的博阿斷裂、阜康斷裂、達爾布特斷裂.

本文結果只是基于現有給出的判定結果，有些判定為處于震間晚期可能性較小的斷層段是由于資料不完備的原因，例如，青藏高原西部地區的地震地質探察、大地測量觀測等基礎資料比較弱.

6 結論與討論

發震斷層的震間晚期判定是強震時間預測的基礎工作，雖然本文嘗試給出中國大陸活動地塊邊界帶主要斷層的震間晚期判定結果，但結果的可靠程度、精細程度等均存在巨大的改善空間.其主要原因在于，目前有些區域的地震監測能力較弱，例如，青藏高原中西部地區的大地測量、測震以及地震地質探察均比較薄弱，另外，現今對于大陸型強震原地復發的動力學過程更多地停留在宏觀概念理論模型.

圖12 中國大陸活動地塊邊界帶主要斷層及其震間晚期判定系數空間分布圖Fig.12 The spatial distribution of major faults in the boundary zone of active tectonic blocks in Chinese mainland and their judgment coefficients of the late stage of interseismic period

早在美國南加州地震中心成立時，地震學家們便期望通過努力構建包括孕震斷層在內的巖石圈應力應變時空演化模型，基于該模型能像氣象預報一樣開展物理預報（Aki，2002）；數十年來，全世界地震學家積累了大量地震學、大地測量、地震地質等野外觀測資料，就斷層運動、斷層應力、破裂模型等等涉及震源物理的眾多基礎研究逐漸深入，使得地震學家在強震物理預測方向有了更多期待（Avouac,2015）.斷層是否處于震間晚期，只是時間預測最初步工作，但對于這個問題，我們基本監測和基礎研究仍有很大改善空間，另一方面，從震間晚期到震前階段的更進一步的時間預測判定，尚需更多的震例積累和理論研究.

基于板塊俯沖百年尺度強震輪回過程給出斷層運動模式，與斷層活動和強震孕育密切相關的地殼變形在時空上是不均勻的，這導致了震間不同階段的判定與強震原地復發的動力學過程密切相關，而中國大陸多數地區7級強震原地復發周期長達數千年（鄧起東等，1994），相關的斷層運動模式是否適用于中國大陸地區，也即描述不同孕震階段斷層運動時空演化特征的中國大陸地區強震斷層運動模式尚未明確.雖然與強震動力學過程有關的現代觀測涉及的頻段越來越寬，但均無法彌補數十年現代地球物理和大地測量觀測時間太短的缺陷.數百年至數千年的強震原地復發間隔正處于上述兩者間隙中，地震位錯、非地震斷層運動和區域地殼變形所形成的永久形變在強震輪回中的時空演化特征很難用現代觀測完整覆蓋（Tectonics Program,2003），因此，對中國大陸地區處于不同孕震階段的斷層段，開展地震地質、大地測量、地震學和地球化學等加密觀測，以“空間換時間”的思路分析其整個強震過程不同階段的斷層運動、斷層應力等震源物理特征，以此構建大陸型震源物理動力學模型，然后逐漸積累更多的資料不斷完善模型，為大陸強震物理預測不斷提供更好的研究基礎.

另外，2021年瑪多7.4級地震發生在巴顏喀拉塊體內部，是我國大陸記錄到為數不多的活動地塊內部的淺源7級地震，而活動地塊內部的地殼變形機制、活動斷層形成與演化、強震孕育發生過程等與活動地塊邊界帶相關問題是不同的.而對于不同區域、不同學科，震源物理模型具體內容存在較大差異，而且隨著監測條件和科學認識不斷改進，既定震源物理模型也不斷調整與完善（Aki，2002；Jordan,2006;2017SCEC）.根據本文斷層震間階段判定過程和相關基礎研究，震源物理模型至少包括：斷層破裂模型、斷層運動模型、斷層應力模型、地震活動模型等，為做好上述模型還需要更基礎的大地測量模型、流變模型、熱力學模型、結構模型、介質模型、地塊模型等（圖13）.只有針對中國大陸活動地塊邊界帶的震源物理模型不斷細化和完善，我國的主要斷層孕震階段的判定才更科學、更可靠，才能有效加強強震預測的物理基礎.

圖13 中國大陸活動地塊邊界帶主要斷層震間晚期判定所需震源物理模型Fig.13 Physics of earthquake source model required for determining late stage of interseismic period of major faults in the boundary zone of active tectonic blocks in Chinese mainland

致謝衷心感謝中國科學院地質與地球物理研究所滕吉文院士、中山大學張培震院士、國家自然災害防治研究院徐錫偉研究員和崔效峰研究員、北京大學沈正康教授、GFZ中心汪榮江研究員、中國地質大學（武漢）熊熊教授和單斌教授、四川地震局聞學澤研究員、中國地震局地震預測研究所江在森研究員和田勤儉研究員、廣東工業大學王華教授、新疆地震局李杰研究員等對本工作的指導，資料的分享、軟件的幫助以及結果的討論；衷心感謝中國地震局科技與國際合作司車時司長、監測預報司馬宏生副司長、張浪平處長、中國地震局地震預測所吳忠良所長、張曉東書記、中國地震局地球物理研究所王曙光處長等的指導、協調與幫助；還有國內外關注地震研究的近百位老師、領導和同事為本次工作提供諸多無私幫助，由于篇幅所限，不能一一列出，在此一并表示感謝.