2022年9月5日四川瀘定6.8級地震及其余震前后b值變化

摘要：

收集2012年1月1日—2023年2月8日四川瀘定地區（29.28°～30.04°N，101.82°～102.28°E）記錄的地震資料，利用最大似然法進行b值的空間掃描，揭示瀘定主震區及周邊區域b值的時空變化特征。計算結果顯示，2022年瀘定6.8級地震前5個月，主震區及其東南區b值出現加速下降變化，可能預示著局部巖石受力變形接近臨界狀態;5.6級強余震前，震中區b值仍然維持低值異常，顯示出區域內的應力水平仍然偏高，預示著可能會發生強余震;5.6級余震發生后，異常區域明顯縮小但未消失，預示著余震可能還將持續，但再次發生強余震的可能性較小。

關鍵詞：

瀘定地震; b值; 極大似然法; 地震危險性

中圖分類號： P315""""" 文獻標志碼：A"" 文章編號： 1000-0844（2024）05-1214-09

DOI：10.20000/j.1000-0844.20230302001

Variation of b-value before and after the Luding MS6.8

earthquake and its aftershocks in Sichuan

Province on September 5， 2022

LI Mengya1，2， ZENG Xianwei1， YAO Huajian2， LI Xinyan1， ZHANG Zhi3， YANG Shuo1

（1. Ningxia Earthquake Agency， Yinchuan 750001， Ningxia， China;

2. School of Earth and Space Sciences， University of Science and Technology of China， Hefei 230026， Anhui， China;

3. No. 327 Geological Team， Bureau of Geology and Mineral Exploration of Anhui Province， Hefei 230026， Anhui， China）

Abstract：

This study presents a collection of seismic data recorded in the Luding area of Sichuan Province， China （29.28°-30.04°N， 101.82°-102.28°E）， covering the period from January 1， 2012， to February 8， 2023. The spatial scanning of the b-value using the maximum likelihood method reveals the spatiotemporal variation characteristics of the b-value in the Luding source area and its surroundings. The results indicate that five months prior to the Luding MS6.8 earthquake in 2022， the b-value in both the source area and the southeastern region of Luding experienced a rapid decrease， suggesting that local rocks were approaching a critical state of stress deformation. Before the MS5.6 strong aftershock， the b-value in the epicentral region remained anomalously low， indicating that the stress level in the area was still high and suggesting the potential for further strong aftershocks. Following the MS5.6 aftershock， the anomalous area noticeably shrank but did not disappear entirely， implying the likelihood of continued aftershocks， though the possibility of another strong aftershock is low.

Keywords：

Luding earthquake; b-value; maximum likelihood method; seismic hazard

0 引言

據中國地震臺網中心測定，北京時間2022年9月5日12時52分，四川省甘孜州瀘定縣發生MS6.8地震（下文簡稱瀘定6.8級地震），震中位于瀘定縣磨西鎮（29.59°N，102.08°E），震源深度16 km，最高烈度達到Ⅸ度。震中區域水系發育，地形高差大，且構造應力較強［1-2］。截至2023年2月8日，共記錄到3.0級以上余震28次，其中3.0～3.9級地震18次，4.0～4.9級地震8次，5.0～5.9級地震2次，分別是2022年10月22日5.0級地震和2023年1月26日5.6級地震。瀘定6.8級地震的震中位于鮮水河斷裂帶，破裂過程以走滑運動為主，兼具部分逆沖分量［3］，破裂方向為NW-SE向，與鮮水河斷裂方向基本一致，最大滑動量達1.84 m，破裂總長度約40 km。

鮮水河斷裂在青藏高原的內部形變中起到了關鍵的作用，是南北地震帶上一條非常活躍的大型左旋走滑構造帶［4］。因此，其整體滑動速率較高，特別是北西段［5］。從歷史地震活動的強度看，7級以上地震大多也發生在鮮水河斷裂帶北西段，特別是1900年以來接連發生過4次7級以上地震，其中最大地震為1973年2月6日四川爐霍7.6級地震，其指示北西段的地震活動強度明顯高于東南段。頻繁發生的強破壞性地震表明，鮮水河斷裂帶北西段及周邊區域應力積累水平較高［6］，其東南段則一直處于強震空白區。此次瀘定6.8級地震恰好發生在鮮水河斷裂帶的磨西段，即東南端，填補了東南段的強震空白區，但該區段后續潛在強震危險性依然非常高，需要持續關注。

研究表明，區域地震活動的震級、頻次服從古登堡-里克特（G-R）關系式，此式也是地震研究中重要的統計學公式。該關系式中b值具有明確的物理意義，是地震危險性分析中非常重要的地震學參數。相關研究表明［7-8］，b值在幾何意義上是地震斷裂結構方式的度量，和地下介質的性質、應力狀態及巖石破裂等有非常密切的關系。因此，可以結合地震的分布情況，對不同的區域開展b值空間掃描，將低b值作為高有效剪應力的判斷指標，指示高應力為閉鎖區，反之為非閉鎖區［9］。國內外眾多學者的研究結果也表明，區域斷裂活動性、強震危險性、地震深度、震源破裂程度等與b值的高低存在較好的相關性［10-13］。眾多的震例回溯研究指示，強震前b值會出現顯著的下降變化［14-17］。因此，b值的高低能衡量一個區域接近巖石破裂強度極限的程度［18］和承受的平均應力水平，是監視破壞性地震孕育的有效手段之一。

本文擬利用瀘定地震震源區及其周邊區域在2012年以來記錄的小震資料，通過分析2022年瀘定6.8級地震及其5.6級強余震前后b值的時空變化特征，研究瀘定地震震源區及周邊區域的應力變化過程。

1 資料選取

以瀘定6.8級地震震源區及周邊區域（29.28°～30.04°N，101.82°～102.28°E）作為研究區，選取2012年1月1日—2023年2月8日的地震資料。目錄通過全國統一編目系統（http：//10.5.160.18/console/exit.action）下載。研究區活動構造復雜，斷裂交匯，是天然地震孕育的優勢場所，沿鮮水河斷裂磨西段地震多發（圖1）。

2012年至瀘定6.8級地震前，研究區以小震活動為主［圖2（a）、（c）］，小震活動頻次不高［月頻次圖2（b）］，偶有起伏。2016年小震活動頻次明顯升高，但強度不大，最大震級為ML4.4。2022年至瀘定6.8級地震發生前，研究區小震活動頻次略有升高，強度依然不大。瀘定6.8級地震發生后，余震活動顯著增強，并先后發生了5.0級和5.6級兩次較強余震，震中分別距離主震5 km和8 km，時間分別發生在主震后47 d和143 d。

2 研究方法

2.1 方法選擇

b值計算通常有2種擬合方法，即最小二乘法和極大似然法。二者最大的區別是最小二乘法更重視含有豐富信息的大地震，給予權重較大，而數量較多的小震反而權重較小;極大似然法則相對較平均，對所有的地震震級用同樣的權重并求平均，即對數量較多的小地震進行信息加權。事實證明，最小二乘法有一定的局限性［4］，本文采用極大似然法［19］，從而得到：

b=lge-MC （1）

式中：MC為最小完備震級;為所有大于MC地震震級的平均震級。b值的誤差估計由改進公式得到［20］，即：

δb=2.3b2∑ni=1（Mi-）2n（n-1） （2）

式中：n為樣本個數，是每一個震級檔在完整時間段內所記錄的地震總數；Mi（i=1，2，3，4，…，n）則表示不同的震級。

在計算b值時，最小完備震級MC是一個非常重要的基礎參數［21］，一般指能夠被臺網完整記錄到的最小地震震級，是一個臨界震級。在b值的實際計算中，為了保證計算的準確性和科學性，一般會舍棄小于MC的地震事件。

2.2 最小完備震級分析

最小完備震級MC指能夠被完全記錄到的最小震級，其和區域地震臺網的布局、密度和人為因素等密切相關［22］。地震臺站空間布局的非均勻性、地震事件分析過程中的人為選擇等都會對地震目錄的完整性造成影響，所以MC往往也存在區域差異。

對研究區以0.02°×0.02°的網格進行掃描，計算區域內MC，結果如圖3所示。主震震中周邊MC在ML0.6左右;震中南端靠近石棉的地區MC較高，介于ML1.2～1.5;康定南偏東區域MC也略高，介于ML0.8～1.0之間。因此，研究區MC整體介于ML0.6～1.5，反映了地震監測能力存在區域差異。

2.3 背景b值計算

前人大量的研究表明，低b值閾值的高低會因研究區構造背景和震源深度等的不同而有所差異［23］。因此，首先要確定研究區域低b值的閾值，即研究區域的背景b值。本文選取了2012年以來研究區全部地震資料，利用極大似然法和標準差計算區域平均b值（公式見2.1節），并將計算結果作為研究區的背景b值。

研究表明，G-R關系式在低震級段和高震級段會出現突然的轉折，即“掉頭”和“擺尾”現象［24］，這種突然的轉折對b值的準確計算會造成一定的影響。因此，在計算背景b值時不但需要考慮最小完備震級的影響，還需要考慮高震級地震的影響，故本文只選擇了ML≤5.0的地震參與計算（圖4）。

為了進一步消除低震級段“掉頭”現象，需要分析b值和截至震級（Mcut-off）的關系，并挑選標準差最小、b值相對穩定的震級作為截至震級。結果顯示（圖4），當Mcut-off ≤ML1.0時，b值標準差隨著截止震級的增大而逐漸減小;當ML1.1≤Mcut-off≤ML2.4時，b值標準差最小，且b值結果也相對穩定;當Mcut-off≥ML2.5時，b值標準差逐漸增大，b值計算結果上下波動，變化較大。因此，當Mcut-off在ML1.1至ML2.4時計算的b值最為合理。

震級-頻度曲線的結果顯示（圖5），當最小完備震級等于ML1.3時，整個曲線擬合度最高。該結果恰好也落在上述b值計算結果穩定時的截至震級合理區間內。通過計算，研究區內b值及其標準差分別為0.95和0.05。綜上所述，b值大小可表示為0.95±0.05。有研究表明，余震個數的增加在一定程度上也會對b值的計算造成一定的影響，筆者刪除了瀘定6.8級地震后所有的余震，并重新計算了b值。結果顯示，去掉余震的b值和沒有去掉余震的b值僅差0.03。因此，在時間尺度夠長，且b值變化不大的情況下，為保證全文資料使用的一致性，本文最終選擇0.9作為低b值異常的閾值。

3 結果分析

3.1 b值時間變化特征分析

為了分析瀘定6.8級地震和5.6級余震前后b值隨時間變化的特征，本文以500個地震樣本為一個時間窗進程，以50個地震樣本的步長逐步滑動樣本窗，分別計算每一個窗口內的b值，并繪制b值時間進程曲線（圖6）。結果表明，2012年1月1日—

2022年4月20日間b值偏高（大于背景值0.9）;2022年4月20日后b值出現明顯下降，降至0.8以下，直至2022年9月5日發生瀘定6.8級地震;震

后b值快速恢復高值，持續一個月左右，b值再次下降，直至2023年1月26日發生5.6級強余震；強余震后，b值再次恢復高值。根據b值的時序變化過程，本文以2022年4月20日、2022年9月5日、2023年1月26日為時間節點，將b值隨時間的變化劃分為4個時段：2012年1月1日—2022年4月20日，2022年4月21日—2022年9月4日，2022年9月5日—2023年1月25日，2023年1月26日—2023年2月8日。

3.2 b值空間變化

根據b值的時序變化及地震發生的時間節點，筆者繼續以0.02°×0.02°的網格尺度，分別計算了4個不同時段研究區的b值（圖7）。結果顯示，2012年1月1日—2022年4月20日，主震震中附近b值偏高，僅存在局部的低b值區，主震以北靠近康定的區域存在明顯的低b值區，顯示該區域應力偏高［圖7（a）］。2012年1月1日—2022年9月4日，相對上一時段，主震東南區域存在明顯的低b值異常區，其他區域b值分布變化不大。2012年1月1日—2023年1月25日，瀘定6.8級地震發生后，震中區附近低b值區擴大，尤其是震中東南區域低b值異常仍在持續。2012年1月1日—2023年2月8日，瀘定5.6級強余震發生后，主震震中區附近低b值異常未出現明顯緩解。總體來看，瀘定6.8級地震和5.6級強余震均發生在低b值與高b值交匯過渡區域，且主震以及強余震發生后，震中區應力仍顯示偏高，可能意味著余震還將持續一段時間。

為了進一步探究以上4個時間段內研究區域的b值變化，詳細分析主震以及強余震前是否存在顯著的b值下降，筆者分別計算了以上4個時段相鄰時段的b值差值（圖8）。結果顯示，瀘定6.8級地震前5個月左右，震中區附近出現小范圍b值下降，而震中東南區域出現大范圍b值下降［圖8（a）］，主震發生后至5.6級強余震前，震中區附近b值下降區范圍擴大且分散在多個區域，震中東南區域b值下降消失［圖8（b）］；5.6級強余震后，b值下降區范圍明顯縮小，僅在該地震震中附近出現小范圍b值下降區［圖8（c）］。由此可以看出：（1）瀘定6.8級地震前5個月，不但震中附近出現應力升高，而且該異常也集中出現在震中東南區域。通過研究瀘定地震的破裂過程發現，主震引起的斷裂面呈NW-SE走向，破裂由震中主要向SE方向延伸。這一認識說明主震東南b值下降區與破裂擴展區及余震展布是吻合的。（2）主震發生后，震中區b值依舊存在下降異常，且范圍較大，顯示區域應力水平依舊偏高，余震將持續發生，并可能發生強余震，這一現象一直持續到2023年1月26日5.6級余震發生前。（3）5.6級余震發生后，b值下降異常區明顯縮小，預示著余震可能還將持續，但再次發生強余震的可能性較小。綜合來看，主震和強余震均發生在b值下降異常區或其邊緣，說明b值下降異常區的地震危險性較其他區域要強，這一結果也與巖石受力變形破壞試驗的結果一致［25-26］。

4 可靠性分析

4.1 低b值閾值的選擇

研究表明，不同區域因構造背景和介質等的差異，低b值閾值也不同［27］。一般認為，b值的低值閾值在0.7左右［28-29］，有時甚至更低［30］。但本文通過計算給出的低b值閾值為0.9，相對較大。分析原因可能與研究區地震震源深度、震源機制及構造等方面有關。由于地殼淺部巖石的壓力值較低、高低速體分布不均、介質非均勻性強，是應力發生的主要構造區，大部分地震主要發生在這一區域內，密集的淺部地震可能會引起區域b值升高［31］。研究區的重定位結果顯示［32］，瀘定地震主震及其余震區內91%以上的地震集中在深度10 km以內，其中6～8 km深度最為集中，平均深度為7.7 km，這一現象可能會導致研究區的低b值閾值偏大。另外，震源機制也是影響b值閾值高低的因素，相關研究結果證實，走滑型地震的b值高于逆斷型地震而低于正斷型地震，約為0.8。而瀘定地震及其余震的震源機制解多為走滑型地震［33-35］，并帶部分逆沖分量，這與前人關于走滑型地震低b值閾值偏高的結果相吻合。

4.2 空間掃描結果分析

b值空間掃描結果的準確性和地震分布的密度密切相關，而地震的發生又受控于活動斷裂。因此，地震的空間分布有極大不均勻性和空間差異性。在實際的b值計算中，需要根據地震的分布設定不同的b值參數。一般而言，掃描參數的設置有固定搜索半徑和固定計算樣本量兩種方式，其中固定計算樣本量根據地震的分布有變化的選擇搜索半徑，而固定搜索半徑顧名思義就是選擇一個固定值作為半徑，其缺點是邊緣或者稀疏區不能滿足樣本量的基本要求。

本文采用固定計算樣本量的方式，即根據地震事件的分布情況選擇不同的搜索半徑。在地震密度越大的區域，設置較小的搜索半徑；相反，地震密度越小，為了保證樣本量則擴大搜索半徑［15，32］，從而獲得更合理的b值分布。在利用2012年1月—2023年2月的地震目錄進行b值空間掃描時，得到了b值標準差的空間分布圖像，掃描尺度為0.02°×0.02°，每個節點參與計算b值的地震數固定為160個，并滿足大于最小完備性震級（MC）的地震數至少為30（圖9）。結果顯示，研究區內b值標準差整體較小，大部分地區都在0.1以下，特別是低b值異常區的標準差，僅在研究區邊緣地震稀疏區域標準差偏大。所以，本文計算的b值空間掃描結果是可靠的。

5 結論

本文選取2012年1月1日—2023年2月8日瀘定6.8級地震震中及其周邊區域（29.28°～30.04°N，101.82°～102.28°E）記錄到的地震資料，選取合適的最小完備震級并利用極大似然法計算了研究區的b值背景。擬合結果顯示，低b值的閾值為0.9。分析閾值偏高的原因可能和研究區的震源深度及震源機制有關。隨后對瀘定6.8級地震前后、5.6級余震前后4個不同時段的b值空間分布特征和空間變化特征進行了詳細分析，主要得到以下幾個結論：

（1） 瀘定6.8級地震前5個月左右，震中區附近出現小范圍b值下降，主要集中在震中東南區域，可能預示著震中附近出現應力顯著升高和局部巖石受力變形接近臨界狀態，與主震引起的斷層面破裂擴展區及余震分布是吻合的。

（2） 主震發生后至5.6級強余震前，震中區b值依舊存在下降異常，且范圍較大，顯示區域應力水平依舊偏高，余震將持續發生，并可能發生強余震。

（3） 5.6級余震發生后，b值下降異常區明顯縮小，預示著余震可能還將持續，但再次發生強余震的可能性較小。

（4） 瀘定6.8級地震和5.6級強余震發生前后b值的變化特征顯示，結合b值平面分布圖和變化圖可以有效判定一個區域的地震危險性。

參考文獻（References）

［1］ 梁明劍，陳立春，冉勇康，等.鮮水河斷裂帶雅拉河段晚第四紀活動性［J］.地震地質，2020，42（2）：513-525.

LIANG Mingjian，CHEN Lichun，RAN Yongkang，et al.Late-Quaternary activity of the Yalahe fault of the Xianshuihe fault zone，eastern margin of the Tibet Plateau［J］.Seismology and Geology，2020，42（2）：513-525.

［2］ 尹福光，潘桂棠，孫志明.西南三江構造體系及演化、成因［J］.沉積與特提斯地質，2021，41（2）：265-282.

YIN Fuguang，PAN Guitang，SUN Zhiming.Genesis and evolution of the structural systems during the Cenozoic in the Sanjiang orogenic belt，Southwest China［J］.Sedimentary Geology and Tethyan Geology，2021，41（2）：265-282.

［3］ 戴丹青，孫麗，楊志高.2022年9月5日四川瀘定MW6.6地震破裂過程研究［J］.地震地磁觀測與研究，2022，43（5）：186-192.

DAI Danqing，SUN Li，YANG Zhigao.Rupture process of the 5 September 2022 MW6.6 Luding earthquake in Sichuan［J］.Seismological and Geomagnetic Observation and Research，2022，43（5）：186-192.

［4］ 曾憲偉，龍鋒，任家琪，等.2019年6月17日長寧MS6.0地震前后b值時空變化分析［J］.地震，2020，40（3）：1-14.

ZENG Xianwei，LONG Feng，REN Jiaqi，et al.Spatial and temporal variation of b value before and after the Changning MS6.0 earthquake on June 17，2019［J］.Earthquake，2020，40（3）：1-14.

［5］ 熊探宇，姚鑫，張永雙.鮮水河斷裂帶全新世活動性研究進展綜述［J］.地質力學學報，2010，16（2）：176-188.

XIONG Tanyu，YAO Xin，ZHANG Yongshuang.A review on study of activity of Xianshuihe fault zone since the Holocene［J］.Journal of Geomechanics，2010，16（2）：176-188.

［6］ 張培震.青藏高原東緣川西地區的現今構造變形、應變分配與深部動力過程［J］.中國科學（D輯：地球科學），2008，38（9）：1041-1056.

ZHANG Peizhen.Current tectonic deformation，strain distribution and deep dynamic process in western Sichuan on the eastern margin of Qinghai—Tibet Plateau［J］.Science in China （Series D：Earth Sciences），2008，38（9）：1041-1056.

［7］ WIEMER S.Minimum magnitude of completeness in earthquake catalogs：examples from Alaska，the western United States，and Japan［J］.Bulletin of the Seismological Society of America，2000，90（4）：859-869.

［8］ WYSS M.Fractal dimension and b-value on creeping and locked patches of the San Andreas fault near Parkfield，California［J］.Bulletin of the Seismological Society of America，2004，94（2）：410-421.

［9］ ZUNIGA F R.Most- and least-likely locations of large to great earthquakes along the Pacific coast of Mexico estimated from local recurrence times based on b-values［J］.Bulletin of the Seismological Society of America，2001，91（6）：1717-1728.

［10］ 易桂喜，聞學澤，范軍，等.由地震活動參數分析安寧河—則木河斷裂帶的現今活動習性及地震危險性［J］.地震學報，2004，26（3）：294-303.

YI Guixi，WEN Xueze，FAN Jun，et al.Assessing current faulting behaviors and seismic risk of the Anninghe—Zemuhe fault zone from seismicity parameters［J］.Acta Seismologica Sinica，2004，26（3）：294-303.

［11］ 易桂喜，聞學澤，蘇有錦.川滇活動地塊東邊界強震危險性研究［J］.地球物理學報，2008，51（6）：1719-1725.

YI Guixi，WEN Xueze，SU Youjin.Study on the potential strong-earthquake risk for the eastern boundary of the Sichuan—Yunnan active faulted-block，China［J］.Chinese Journal of Geophysics，2008，51（6）：1719-1725.

［12］ AMITRANO D.Brittle-ductile transition and associated seismicity：experimental and numerical studies and relationship with the b value［J］.Journal of Geophysical Research （Solid Earth），2003，108（B1）：2044.

［13］ CHEN C C，WANG W C，CHANG Y F，et al.A correlation between the b-value and the fractal dimension from the aftershock sequence of the 1999 Chi-Chi，Taiwan，earthquake［J］.Geophys J Int，2006，167：1215-1219.

［14］ 張帆，韓曉明，陳立峰，等.鄂爾多斯地塊北緣b值的時空特征及其地震預測效能分析［J］.地震學報，2018，40（6）：785-796，832.

ZHANG Fan，HAN Xiaoming，CHEN Lifeng，et al.Spatio-temporal characteristics of b value in the northern margin of Ordos block and its earthquake prediction efficiency［J］.Acta Seismologica Sinica，2018，40（6）：785-796，832.

［15］ 韓佳東，楊建思，王偉平.2017年西藏米林MS6.9地震余震序列重定位和b值時空分布特征［J］.地震學報，2019，41（2）：169-180，277.

HAN Jiadong，YANG Jiansi，WANG Weiping.Relocation of the aftershock sequence of Milin MS6.9 earthquake in 2007 and spatio-temporal variation characteristics of b-value［J］.Acta Seismologica Sinica，2019，41（2）：169-180，277.

［16］ 曾憲偉，李文君，馬翀之，等.基于b值分析寧夏吳忠—靈武地區強震危險性［J］.地震研究，2021，44（1）：41-48.

ZENG Xianwei，LI Wenjun，MA Chongzhi，et al.Strong earthquake risk in Wuzhong—Lingwu region of Ningxia based on b-value［J］.Journal of Seismological Research，2021，44（1）：41-48.

［17］ 曾憲偉，趙小艷，李蒙亞，等.2021年5月21日漾濞4次MS≥5.0地震前后b值變化［J］.地震學報，2022，44（3）：401-412.

ZENG Xianwei，ZHAO Xiaoyan，LI Mengya，et al.Variation of b-value before and after Yangbi four MS≥5.0 earthquakes on May 21，2021［J］.Acta Seismologica Sinica，2022，44（3）：401-412.

［18］ 李全林，陳錦標，于淥，等.b值時空掃描：監視破壞性地震孕育過程的一種手段［J］.地球物理學報，1978，21（2）：101-125.

LI Quanlin，CHEN Jinbiao，YU Lu，et al.Time and space scanning of the b-value：a method for monitoring the development of catastrophic earthquakes［J］.Chinese Journal of Sinica，1978，21（2）：101-125.

［19］ AKI K.Maximum likelihood estimate of b in the formula logn=a-bm and its confidence limits［J］.Bull Earthquake Res Inst，Tokyo Univ.1965，43：237-239.

［20］ SHI Y L，BOLT B A.The standard error of the magnitude-frequency b value［J］.Bulletin of the Seismological Society of America，1982，72（5）：1677-1687.

［21］ WIEMER S.A software package to analyze seismicity：ZMAP［J］.Seismological Research Letters，2001，72（3）：373-382.

［22］ 易桂喜，聞學澤，辛華，等.龍門山斷裂帶南段應力狀態與強震危險性研究［J］.地球物理學報，2013，56（4）：1112-1120.

YI Guixi，WEN Xueze，XIN Hua，et al.Stress state and major-earthquake risk on the southern segment of the Longmen Shan fault zone［J］.Chinese Journal of Geophysics，2013，56（4）：1112-1120.

［23］ 陳時軍，馬麗.地震活動性統計學研究方法的若干進展與討論［J］.地震，2002，22（1）：32-42.

CHEN Shijun，MA Li.Review on the application of statistics in seismicity［J］.Earthquake，2002，22（1）：32-42.

［24］ 楊文，劉杰，解孟雨，等.2022年9月5日四川瀘定MS6.8地震序列重定位研究［J］.中國地震，2022，38（4）：622-631.

YANG Wen，LIU Jie，XIE Mengyu，et al.Study on relocation of the September 5，2022 Luding MS6.8 earthquake［J］.Earthquake Research in China，2022，38（4）：622-631.

［25］ 劉力強，馬勝利，馬瑾，等.不同結構巖石標本聲發射b值和頻譜的時間掃描及其物理意義［J］.地震地質，2001，23（4）：481-492.

LIU Liqiang，MA Shengli，MA Jin，et al.Temporal scanning of b-value and spectrum of AE activity for samples with different textures and their physical implications［J］.Seismology and Geology，2001，23（4）：481-492.

［26］ AMORSE D ，GRASSO J R，RYDELEK P A.On varying b-values with depth：results from computer-in-tensive tests for Southern California［J］.Geophys J Int，2010，180：347-360.

［27］ 朱艾斕，徐錫偉，甘衛軍，等.鮮水河—安寧河—則木河斷裂帶上可能存在的凹凸體：來自背景地震活動性的證據［J］.地學前緣，2009，16（1）：218-225.

ZHU Ailan，XU Xiwei，GAN Weijun，et al.The possible asperities on the Xianshuihe—Anninghe—Zemuhe fault zone：evidence from background seismicity［J］.Earth Science Frontiers，2009，16（1）：218-225.

［28］ 易桂喜，聞學澤，張致偉，等.川南馬邊地區強震危險性分析［J］.地震地質，2010，32（2）：282-293.

YI Guixi，WEN Xueze，ZHANG Zhiwei，et al.Study on potential strong earthquake risk in Mabian area，southern Sichuan［J］.Seismology and Geology，2010，32（2）：282-293.

［29］ 韓曉明，薛丁，趙星.內蒙古東部扎蘭屯地區的b值時空變化特征［J］.地震，2012，32（3）：142-149.

HAN Xiaoming，XUE Ding，ZHAO Xing.Spatiotemporal b-value variation characteristics in the Zhalantun area of the eastern part of Inner Mongolia［J］.Earthquake，2012，32（3）：142-149.

［30］ 謝卓娟，李山有，呂悅軍.滇西南地區主要活動斷裂的b值空間分布特征［J］.地球科學，2015，40（10）：1755-1766.

XIE Zhuojuan，LI Shanyou，L Yuejun.b values spatial distribution characteristics of the main active faults in southwestern Yunnan［J］.Earth Science，2015，40（10）：1755-1766.

［31］ MORI J，ABERCROMBIER E.Depth dependence of earthquake frequency magnitude distributions in California：implications for the rupture initiation［J］.Journal of Geophysical Research，1997，102（b7）：15801-15090.

［32］ CAI G Y，WANG W L，WU J P，et al.Three-dimensional body wave velocity structure and seismogenic structure for the 2022 MS6.8 Luding earthquake sequence in Sichuan，China［J］.Frontiers in Earth Science，2023，11：1099744.

［33］ SCHORLEMMER D，WIEMER S，WYSS M.Variations in earthquake-size distribution across different stress regimes［J］.Nature，2005，437（7058）：539-542.

［34］ GULIA L，WIEMER S.The influence of tectonic regions on the earthquake size distribution：a case study for Italy［J］.Geophysical Research Letters，2010，37（10）：L10305.

［35］ YANG Z G，DAI D Q，ZHANG Y，et al.Rupture process and aftershock focal mechanisms of the 2022 M6.8 Luding earthquake in Sichuan［J］.Earthquake Science，2022，35（6）：474-484.

（本文編輯：任 棟）