摘要：基于云南地震臺網(wǎng)和騰沖火山臺網(wǎng)的寬頻帶地震波形資料和震相報告文件以及川滇地區(qū)高分辨率公共速度模型2.0版，采用結(jié)合波形互相關(guān)技術(shù)的雙差相對定位法對2023年云南芒市ML4.0震群序列中442次地震進行重定位，利用CAP方法反演了地震序列中ML≥3.5地震震源機制解，研究了芒市ML4.0震群序列的發(fā)震機理，探討該震群是前震序列還是普通震群。雙差定位法獲得了363次地震的精定位結(jié)果，震中走向呈NE向，地震序列呈三角形形態(tài)分布。CAP波形反演ML4.0地震最佳雙力偶機制解為節(jié)面I走向45°，傾角70°，滑動角-61°；節(jié)面II走向167°，傾角35°，滑動角-143°，矩震級為MW3.91，最佳矩心深度為7 km。綜合精定位后地震序列的時空分布特征、震源機制解信息、S-P走時差、區(qū)域地質(zhì)構(gòu)造和地震活動性探討此次地震的發(fā)震構(gòu)造及孕震環(huán)境，結(jié)果表明：此次震群不是前震序列，而是一次普通的3～4級震群活動。地震序列在不同階段的時空演化過程可能是由至少2條存在一定夾角的斷層結(jié)構(gòu)在高應(yīng)力和區(qū)域的溫泉地?zé)峄顒右l(fā)。

關(guān)鍵詞：芒市震群；雙差定位；CAP反演；地震活動性；流體擴散

中圖分類號：P315.72"" 文獻標識碼：A"" 文章編號：1000-0666（2025）01-0010-12

doi：10.20015/j.cnki.ISSN1000-0666.2025.0002

0 引言

震群是指一組在空間和時間上叢集發(fā)生的地震序列，不同于常見的主-余型地震序列，震群內(nèi)主要地震事件的震級大小相當(dāng)、缺乏明顯的主震事件。根據(jù)陸遠忠等（1984）和《測震分析預(yù)測技術(shù)方法工作手冊》（中國地震局監(jiān)測預(yù)報司，2020）給出的震群定義，震群內(nèi)的地震事件應(yīng)在空間分布、頻度以及震級上滿足以下標準：①半徑≤20 km，與外圍地震的分布有較明顯的界限；②最大日頻次≥3次，總次數(shù)大于等于10次；③序列中最大地震的震級ML≤5.2，最大與次大地震震級差≤1.1；④至少有3個震級相差不大（震級差lt;0.6）的主要地震。

已有的觀測和研究結(jié)果表明，震群序列通常發(fā)生在流體侵入地區(qū)（如火山區(qū)、水庫區(qū)、高壓注水區(qū)等）、斷層深部無震滑移形變區(qū)和（或）多組斷層構(gòu)造交接復(fù)合的地殼介質(zhì)不均勻且處于高應(yīng)力狀態(tài)等構(gòu)造區(qū)域（Chen et al，2011；Shelly et al，2016；Dublanchet，De Barros，2020）。因此針對震群序列的研究可為深入認識震區(qū)斷層構(gòu)造、流體運移、應(yīng)力狀態(tài)以及地震危險性等問題提供重要科學(xué)依據(jù)。

震群的孕育發(fā)生是震源區(qū)及周邊區(qū)域地殼活動、應(yīng)力增強的客觀反映，一部分震群活動對周邊區(qū)域中強地震的發(fā)生具有一定的中短期指示意義（即“前兆震群”或“前震序列”），在此類震群發(fā)生后的一定時空范圍內(nèi)（一般3～12個月、距離震群幾十至數(shù)百千米）會發(fā)生破壞性中強地震（竇喜英等，2020）；而另一部分震群活動后周邊區(qū)域并沒有中強地震發(fā)生，即震群序列活動一段時間后就逐漸衰減、停止。然而究竟是哪些因素影響、控制震群周邊地區(qū)的中強地震活動，以及震群活動是否與周邊破壞性地震的發(fā)生有直接因果關(guān)系等問題目前還沒有明確定論。

目前針對震群型地震序列與后續(xù)大震發(fā)生的相關(guān)性研究主要有兩類：一類是根據(jù)震群地震序列從時、空、強表現(xiàn)出來的特征總結(jié)了判定震群是屬于前震序列還是一般性小震群的特征判斷指標，如U值（陸遠忠等，1984），K值（朱傳鎮(zhèn)，王林瑛，1989），ρ值（王煒，楊德志，1987），震源分布、b值和活動頻次等參數(shù)（王光明等，2024）；另一類是對一個大區(qū)域較大范圍的眾多地震，從總體上研究他們的特征（朱傳鎮(zhèn)等，1981；竇喜英等，2020）。然而，最近針對美國南加州地震目錄中的前震與震群地震活動性參數(shù)（地震活動頻次、時空分布等）的研究結(jié)果表明該地區(qū)的震群和前震在地震活動性方面可能沒有明顯區(qū)別（Zaccagnino et al，2024）。

據(jù)云南地震臺網(wǎng)正式測定，2022年12月12日—2023年7月12日在云南省德宏州芒市發(fā)生了一系列ML≥3.0中小地震，并伴隨發(fā)生多次小震活動。截至2023年7月31日該區(qū)域共發(fā)生地震442次，其中MLlt;2.0地震387次，2.0≤MLlt;3.0地震45次，3.0≤MLlt;4.0地震9次，ML≥4.0地震1次。此次芒市地震序列中最大震級為ML4.0，次大震級為ML3.9，兩者震級差為0.1級，中間未出現(xiàn)超過15日無ML≥1.5地震發(fā)生的情況，這些特征符合《測震分析預(yù)測技術(shù)方法工作手冊》（中國地震局監(jiān)測預(yù)報司，2020）中關(guān)于震群的定義，為了研究和表述方便本文將其命名為2023年芒市ML4.0震群。

為了深入研究2023年芒市ML4.0震群序列的發(fā)震機理，并從震源參數(shù)角度探討該震群是前震序列還是普通震群，本文基于寬頻帶數(shù)字地震波形和區(qū)域速度模型，采用結(jié)合波形互相關(guān)技術(shù)的雙差相對定位法（Waldhauser，2000）對此次芒市ML4.0震群地震序列進行重定位，利用CAP方法（Zhao，Helmberger，1994；Zhu，Helmberger，1996）反演了地震序列ML≥3.5地震震源機制解，綜合震源位置、震源機制解等震源參數(shù)、震區(qū)地質(zhì)構(gòu)造特征和歷史地震活動，探討此次ML4.0震群的成因機理，并初步研判其與臨近區(qū)域未來潛在中強地震的活動關(guān)系。

1 數(shù)據(jù)選取

2023年芒市ML4.0震群震中附近200 km范圍內(nèi)布設(shè)有云南地震臺網(wǎng)17個固定地震臺站（其中有4個火山臺網(wǎng)子臺）、12個騰沖火山臺網(wǎng)子臺（圖1、圖2）。2023年1月9日ML4.0地震發(fā)生后，騰沖地震監(jiān)測中心站在震源區(qū)又布設(shè)了2個短周期地震臺用于加密觀測。云南地震臺網(wǎng)提供的目錄顯示，2022年12月24日—2023年3月2日的ML≥3.0地震震源深度均接近20 km（表1），與中國地震臺網(wǎng)正式測定和公布的10 km震源深度存在明顯差異。

本文利用CAP方法反演2022年12月—2023年7月研究區(qū)內(nèi)ML≥3.5地震事件（表1）的震源機制解和最佳矩心深度，選取震中距35～300 km的寬頻帶數(shù)字波形，基于給定的2種地殼速度模型進行反演（圖3a）：①根據(jù)Brian Chen和斯坦福大學(xué)地殼地球物理實驗室分享的Matlab腳本https：//igppweb.ucsd.edu/～gabi/crust1.html#smatlab.，提?。?4.5°N，98.5°E）區(qū)域內(nèi)的全球地殼模型CRUST1.0速度模型，記為模型A（Laske et al，2013）；②川滇地區(qū)高分辨率公共速度模型2.0版本（SWChinaCVM-2.0https：//github.com/liuyingustc/SWChinaCVM-V2.0.，記為模型B（Liu et al，2024）。

本文收集整理了2022年12月—2023年7月云南地震臺網(wǎng)記錄到的發(fā)生在芒市ML4.0主震震中位置（24.270°N，98.632°E）附近20 km范圍內(nèi)的442次地震（圖2）的觀測報告，并根據(jù)地震觀測報告中的發(fā)震時刻、震中、震級等信息，補充了騰沖火山臺網(wǎng)正式編目觀測報告中的震相到時資料（圖3b）。同時，從云南地震臺網(wǎng)和騰沖火山臺網(wǎng)波形資料中選取至少同時有3個臺站記錄到的地震事件波形用于波形互相關(guān)分析，以進一步提高后續(xù)的地震定位精度。

2 研究結(jié)果

2.1 CAP方法反演震源機制解及矩心深度

挑選2023年芒市ML4.0震群中4次ML≥3.5地震（表1）在震中距300 km范圍內(nèi)的寬頻帶數(shù)字波形，利用云南地震臺產(chǎn)出的觀測報告在地震波形中寫入P、S震相到時，對其做去均值、去線性趨勢及去除儀器響應(yīng)等數(shù)據(jù)預(yù)處理后，將速度記錄從NS-EW-UD分量旋轉(zhuǎn)到大圓弧路徑上，從而獲得R-T-Z分量的波形記錄；分別利用模型A和模型B，采用F-K方法（Zhu，Rivera，2002）計算2種速度模型在震中距35～300 km、震源深度1～25 km的格林函數(shù)庫用于后續(xù)的CAP反演工作（圖4）。

在CAP反演過程中，將Pnl波和面波的截取窗長分別設(shè)為30 s和70 s，反演權(quán)重分別設(shè)為1和0.5。為了提高波形信噪比，并考慮到3次地震事件的震級均不超過ML4.0，本文在反演過程中采用相對高頻的濾波頻段（李姣等，2021），即對Pnl部分以0.15～0.3 Hz、面波部分以0.08～0.15 Hz進行4階Butterworth帶通濾波，最后在斷層走向（0°～360°）、傾角（0°～90°）、滑動角（-90°～90°）和震源深度（1～25 km）空間范圍內(nèi)搜索最佳的震源機制解和震源矩心深度。經(jīng)過反演，本文得到了3次ML≥ 3.5地震的震源機制解和震源深度（表2）。由于2023年4月23日ML3.9地震發(fā)生前10 s發(fā)生了一次ML3.3地震，兩次地震波形疊加，故未得到該地震準確可靠的震源機制解反演結(jié)果。

CAP方法反演表明，本文得到的3次地震的最佳矩心深度（5～12 km），比云南地震臺網(wǎng)地震目錄中的震源深度（20 km左右）要淺許多（表2），并且在同一速度模型下反演得到3次地震的震源機制解存在一定差異（最大相差約20°）。對比分析此次地震在2種速度模型下的CAP反演擬合誤差，以及對理論波形（圖5中紅色波段）和實際觀測波形（圖5中黑色波段）的互相關(guān)系數(shù)、不同波段滑移時間量等表明，基于模型B的反演結(jié)果更優(yōu)，表明該模型應(yīng)更接近研究區(qū)地下真實速度結(jié)構(gòu)，故本文采用模型B用于后續(xù)研究。

2.2 震群序列的雙差重定位

精確的震源位置對于地球內(nèi)部結(jié)構(gòu)和地震孕育過程認知、斷層精細結(jié)構(gòu)探測、地震預(yù)警研究、震后預(yù)測、救援、資源能源勘探開發(fā)等眾多科學(xué)技術(shù)問題至關(guān)重要（侯新榮等，2024）。目前云南地震臺網(wǎng)和騰沖火山臺網(wǎng)日常分析使用的地震定位方法是經(jīng)典的單純型絕對定位法，該方法的優(yōu)點是計算速度快、適用范圍廣，缺點是定位結(jié)果依賴于區(qū)域速度模型的準確性，且實際工作中震源深度反演結(jié)果不穩(wěn)定。近年來應(yīng)用廣泛的雙差定位法，以觀測報告數(shù)據(jù)和（或）波形互相關(guān)數(shù)據(jù)的P、S震相走時差作為輸入數(shù)據(jù)，對速度模型的依賴相對較小；波形互相關(guān)技術(shù)可將P、S震相走時差精確到毫秒級，在臺站密集區(qū)域地震之間的相對定位誤差可降低到幾十米甚至米級范圍（房立華等，2013；王清東，2015；李姣等，2020；王光明等，2022；Liu et al，2024）。

考慮到此次震群中有部分地震事件被騰沖火山臺網(wǎng)記錄編目，本文將云南地震臺網(wǎng)和騰沖火山臺網(wǎng)的觀測報告合并后再進行重定位研究，具體原則及流程為：針對同一地震事件，在云南地震臺網(wǎng)震相到時資料的基礎(chǔ)上加入騰沖火山臺網(wǎng)震相到時資料，將2個報告中地震事件發(fā)震時刻統(tǒng)一后獲得火山臺站針對該地震事件記錄到的Pg和Sg震相走時，經(jīng)過處理得到了10 406條P波和8 759條S波走時差數(shù)據(jù)用于雙差定位。本文還根據(jù)合并后的地震觀測報告整理事件波形用于波形互相關(guān)分析，參考姜金鐘等（2016）和Jiang 等（2019）在云南地區(qū)地震事件波形互相關(guān)分析中的參數(shù)設(shè)置，選擇P、S波到時前0.5 s、后2.5 s的窗長波形用于互相關(guān)分析，并挑選互相關(guān)系數(shù)大于0.7的走時差計算結(jié)果，經(jīng)過互相關(guān)分析，本文得到了27 003條P波和12 712條S波走時差數(shù)據(jù)用于雙差定位。

本文選用速度模型B，對2023年芒市ML4.0震群442個多臺記錄到的地震事件進行了結(jié)合波形互相關(guān)的雙差定位。考慮到此次震群周邊臺站分布相對較為稀疏，雙差定位時設(shè)置最小連接數(shù)和最小觀測數(shù)為6，震源間距小于10 km，事件對到臺站的最大距離小于200 km；設(shè)定 P 波到時的權(quán)重為 1.0，S 波震相到時的權(quán)重為 0.5，經(jīng)過2輪共20次迭代反演，本文最終得到363個地震事件精確的相對位置，重定位成功率為82.1%。圖6給出了重定位后2023年芒市ML4.0震群的震中及震源深度分布圖，并沿震群中最大的ML4.0地震的震源機制解節(jié)面Ⅰ走向（45°）繪制了EE′剖面線、以節(jié)面Ⅱ走向（135°）繪制了AA′剖面線（圖6a），從圖中可以看出震中位置以類似三角形的形態(tài)分布。根據(jù)AA′剖面線提取數(shù)據(jù)繪制了距離隨震源深度的分布關(guān)系（圖6b），震源深度主要集中于6～12 km，幾次震級較大的地震的震源深度都較深，在12～13 km，且震源位置在深度上呈一定角度的交叉線性分布。綜合分析認為，此次震群的發(fā)震構(gòu)造可能為至少兩條存在一定夾角的斷層。

3 討論

3.1 震群的形成與區(qū)域地質(zhì)構(gòu)造和中強地震的關(guān)系

大量觀測和研究表明，震群活動與周邊地區(qū)中強地震的發(fā)生存在一定聯(lián)系，是地殼活動增強的客觀反映（竇喜英等，2020）。芒市ML4.0震群發(fā)生于龍陵—瑞麗斷裂、畹町?dāng)嗔押团瓟嗔褔薜娜菂^(qū)域內(nèi)，該區(qū)1976年發(fā)生了龍陵MS7.3和MS7.4雙震型地震（圖1，紫色震源機制解位置），朱傳鎮(zhèn)等（1981）和竇喜英等（2020）研究發(fā)現(xiàn)中強地震發(fā)生后周圍區(qū)域因應(yīng)力調(diào)整容易誘發(fā)震群活動。但查閱云南歷史地震目錄后發(fā)現(xiàn)，截至芒市ML4.0震群發(fā)生前這個區(qū)域并沒有類似的震群形成過，僅在與怒江斷裂交匯的區(qū)域偶有ML3.0左右小地震發(fā)生。

在地震活動性研究及地震預(yù)報震情跟蹤工作中，前震活動是研究人員認為最有效的短臨預(yù)測指標之一，并且前震序列通常以震群形式出現(xiàn)，如1975年海城7.3級地震和2021年漾濞6.4級地震震前幾天至幾小時內(nèi)都有一系列3～4級震群活動（Jones，1982；Zhou et al，2022）。然而識別前震活動也是一個世界性難題，目前研究人員常用震源機制一致性（陳颙，1978）、精確的震源時空遷移特征（Zhou et al，2022）、b值（Wetzler et al，2023）等數(shù)字地震學(xué)方法研究識別前震活動。本文利用區(qū)域?qū)掝l帶數(shù)字波形和最新速度模型得到了2023年芒市ML4.0震群的精確震源位置和較大地震的震源機制解等參數(shù)，結(jié)果表明此次地震序列中3次ML≥3.5地震的震源機制解差異明顯（同一反演模型下最大相差約20°）。

本文根據(jù)震群中5次ML≥3.4地震（表1）的發(fā)震時間將整個震群活動分為6個階段（整個震群持續(xù)時間為213天，第一階段為第一個地震至表1中3號事件開始前，其他5個階段均是根據(jù) 表1中編號為3、4、5、8、10的事件作為起始時間進行分段的，圖8a～f分段同理），計算每個地震事件被震中距最近的芒市臺（MAS）記錄到的S震相和P震相的到時差（S-P分布）。結(jié)果顯示震群中地震事件在空間分布上逐漸分散、一致性逐漸降低（圖7），震源機制一致性分析結(jié)果表明此次震群是一次普通震群活動，即震源區(qū)后續(xù)一段時間內(nèi)發(fā)生震級更大的破壞性地震的可能性較低。

盡管2023年芒市ML4.0震群結(jié)束之后119天（2023年12月2日），在距ML4.0地震震中位置約55 km處發(fā)生了芒市MS5.0破壞性地震（圖1中紅色震源機制解根據(jù)顧慧冬等（2024）研究結(jié)果繪制），震源機制解結(jié)果顯示該地震主要以走滑型為主，由于該地震震中位置緊鄰龍江水庫，顧慧冬等（2024）研究認為，流體在庫區(qū)原有具備發(fā)生中強地震潛力斷層的裂隙中滲透，促使了該地震的發(fā)生；而2023年芒市ML4.0震群震源機制解以正斷為主，同時精定位結(jié)果顯示此次震群的發(fā)震構(gòu)造可能為至少2條存在一定夾角的斷層。因此，綜合對比兩次地震的震源機制解和震源位置等信息后初步認為，兩次地震序列的發(fā)生沒有直接關(guān)聯(lián)，但兩次地震序列密切的時空關(guān)系反映了芒市及其臨近區(qū)域可能確實處于應(yīng)力高度集中狀態(tài)，后續(xù)應(yīng)加強對該區(qū)域地震活動的監(jiān)測研究。

3.2 震群地震序列的遷移與流體的相關(guān)性

為了更好地分析此次震群的時空演化特征，對整個震群活動6個階段重定位后的震源位置分別進行分析（圖8a～f），結(jié)果顯示第一階段的地震事件主要發(fā)生在一條走向為北東方向的發(fā)震構(gòu)造上（BB′，圖8a）；自第二階段開始，在震源區(qū)南部的另一條發(fā)震構(gòu)造（CC′）上逐漸有地震活動（圖8b）；自第四階段開始，地震在兩條發(fā)震構(gòu)造上交替活動，并由震源區(qū)北東方向逐漸向震源區(qū)西南方向遷移和集中（圖8d～f），且5次ML≥3.4地震均發(fā)生在該區(qū)域。綜合分析認為，BB′斷裂上第一階段的地震活動，使得與之共軛的CC′斷裂應(yīng)力平衡狀態(tài)被打破，誘發(fā)了CC′斷裂上的第一個ML3.4地震。結(jié)合圖8g可明顯看到ML3.4地震的震源深度在逐漸加深，隨著時間的推移，兩條斷裂應(yīng)力不斷調(diào)整，接著觸發(fā)了第三階段和第四階段地震事件的發(fā)生，此時ML≥3.0地震的震源深度在12～14 km，在ML4.0地震發(fā)生后，明顯看到余震主要在兩條斷裂之間活動。第五和第六階段之后震源深度逐漸變淺。震級較大的地震幾乎都發(fā)生在CC′斷裂上，且震源深度比BB′斷裂上的地震深，不同階段的震源時空演化特征表明了該區(qū)域復(fù)雜的發(fā)震構(gòu)造和應(yīng)力積累情況。

震群的時空分布與流體的擴散有較強的相關(guān)性（Shelly et al，2016；肖陽等，2024），鉆孔注水誘發(fā)地震和水庫誘發(fā)地震證實了震群時空分布符合流體孔隙壓力擴散特征（Shapiro et al，1997，2003；Talwan，Acree，1984）。本文基于Shapiro等（1997）提出的流體擴散公式（r=4πDt，其中r表示距第一個地震事件的距離；t表示時間；D表示擴散系數(shù)）和雙差重定位結(jié)果分析了本次震群是否存在流體驅(qū)動機制。根據(jù)流體擴散公式計算擬合得到的擴散系數(shù)D為0.578 7 m2/s，該值在地震發(fā)生后第3～70天的擬合效果最佳（圖9a），表明流體可能參與了芒市ML4.0震群序列的時空遷移，是這次震群形成的主要驅(qū)動力（Scholz，2019）；在70天以后r-t關(guān)系不再符合流體擴散公式，擬合曲線更趨于線性（圖9b），表明70天之后的地震活動可能沒有流體的參與或者流體的影響變得很微弱。

此次芒市ML4.0震群震源區(qū)及鄰近地區(qū)（“小滇西”地區(qū)）是云南地區(qū)震群活動最頻繁的地區(qū)，自有記錄以來該地區(qū)就發(fā)生了多次震群和（或）雙震型地震序列，如1928—1933年發(fā)生了騰沖7次6級以上地震和一系列5級以上地震、1976年龍陵7.3級和7.4級地震（王彬，鄧瑞生，2021）、2008、2011、2014和2021年在盈江地區(qū)陸續(xù)發(fā)生了4次5級地震和2次6級地震活動。該地區(qū)還是中國大陸地區(qū)大地?zé)崃鞲咧祬^(qū)（汪集旸，黃少鵬，1990），區(qū)域內(nèi)高溫溫泉和巖漿活動廣泛分布，是震群序列活動的典型構(gòu)造地區(qū)（Shelly et al，2016；Yoshida et al，2017；Liu et al，2024）。綜合地震震源位置、震源機制解、地震活動的時空演化等研究結(jié)果和溫泉地?zé)帷^(qū)域構(gòu)造特征等資料，本文初步分析認為此次震群發(fā)生在至少2條存在一定夾角、處于高應(yīng)力狀態(tài)的斷層上，地下流體活動可能促使了地震序列在不同階段的時空演化過程。

4 結(jié)論

本文基于云南地震臺網(wǎng)和騰沖火山臺網(wǎng)的寬頻帶地震波形資料和震相報告文件，利用雙差定位方法和基于密集流動臺陣構(gòu)建的川滇地區(qū)高分辨率公共速度模型2.0版本提供的芒市及鄰區(qū)地殼速度模型，對2023年芒市ML4.0震群地震序列進行重定位。利用CAP方法對序列中ML≥3.5地震的震源機制解及震源深度進行反演，得到了震源機制解、最優(yōu)震源深度和矩震級等參數(shù)。綜合震群序列震源分布、震源機制解和區(qū)域地質(zhì)構(gòu)造等研究結(jié)果和資料，得出以下主要結(jié)論：

（1）CAP方法反演結(jié)果表明，3次ML≥ 3.5地震的震源機制解相差較大（同一反演模型下最大相差約20°）；雙差定位結(jié)果表明震群在水平方向和垂直方向均呈一定角度的線性分布，表明此次震群的發(fā)震構(gòu)造可能為至少2條存在一定夾角的斷層。

（2）重定位結(jié)果、震源機制和S-P走時差等反映的震源機制一致性表明，此次震群不是前震序列，而是一次普通3～4級震群活動，即震源區(qū)后續(xù)一段時間內(nèi)不太可能發(fā)生震級更大的破壞性地震。

（3）綜合地震震源位置、震源機制解、地震活動的時空演化等研究結(jié)果和區(qū)域溫泉地?zé)帷?gòu)造特征等資料，本文初步分析認為此次震群發(fā)生在至少2條存在一定夾角、處于高應(yīng)力狀態(tài)的斷層上，地下流體活動可能促使了地震序列在不同階段的時空演化過程。

參考文獻：

安曉文，常祖峰，毛先進.2018.云南第四紀活動斷裂暨《云南第四紀活動斷裂分布圖》［M］.北京：地震出版社.An X W，Chang Z F，Mao X J，et al.2018.Quaternary active faults in Yunnan and distribution map of Quaternary active faults in Yunnan［M］.Beijing：Seismological Press.（in Chinese）

陳颙.1978.用震源機制一致性作為描述地震活動性的新參數(shù)［J］.地球物理學(xué)報，21（2）：146-159.Chen Y.1978.Consistency of focal mechanism as a new parameter in describing seismic activity［J］.Chinese Journal of Geophysics，21（2）：142-159.（in Chinese）

竇喜英，王恩利，張輝，等.2020.甘肅地區(qū)前兆震群對省內(nèi)及邊鄰中強地震的預(yù)測意義研究［J］.地球物理學(xué)進展，35（5）：1651-1658.Dou X Y，Wang E L，Zhang H，et al.2020.Study on the prediction significance of the precursory earthquake swarms to mid-strong earthquakes in Gansu and its adjacent area［J］.Progress in Geophysics，35（5）：1651-1658.（in Chinese）

房立華，吳建平，王未來.2013.四川蘆山 MS7.0 級地震及其余震序列重定位［J］.科學(xué)通報，58：1901-1909.Fang L H，Wu J P，Wang W L.2013.Relocation of the Lushan MS7.0 earthquake and its aftershock sequence in Sichuan［J］.Chinese Scientific Bulletin，58：1901-1909.（in Chinese）

顧慧冬，姜金鐘，李姣，等.2024.2023年云南芒市MS5.0地震震源深度測定及其發(fā)震機理探討［J］.地震學(xué)報，46，doi：10.11939/jass.20240026.Gu H D，Jiang J Z，Li J，et al.2024.Focal depth determination of the 5.0 magnitude earthquake in Mangshi，Yunnan，2023，and discussion of its seismicity mechanism［J］.Acta Seismologica Sinica，46，doi：10.11939/jass.20240026.（in Chinese）

侯新榮，郭振威，郭大雄，等.2024.地震定位方法最新進展綜述［J］.地球物理學(xué)進展，39（3）：959-974.Hou X R，Guo Z W，Guo D X，et al.2024.Review of recent advances in seismic location methods［J］.Progress in Geophysics，39（3）：959-974.（in Chinese）

姜金鐘，付虹，陳棋福.2016.位于構(gòu)造活躍區(qū)的小灣水庫地震活動特征——基于地震精定位的分析［J］.地球物理學(xué)報，59（7）：2468-2485.Jiang J Z，F(xiàn)u H，Chen Q F.2016.Characteristics of seismicity of the Xiaowan reservoir in an area of active tectonics from double-difference relocation analysis［J］.Chinese Journal of Geophysics，59（7）：2468-2485.（in Chinese）

李姣，姜金鐘，楊晶瓊.2020.2017 年漾濞 MS4.8 和 MS5.1 地震序列的微震檢測及重定位［J］.地震學(xué)報，42（5）：527-542.Li J，Jiang J Z，Yang J Q.2020.Microseismic detection and relocation of the 2017 MS4.8 and MS5.1 Yangbi earthquake sequence，Yunnan［J］.Acta Seismologica Sinica，42（5）：527-542.（in Chinese）

李姣，姜金鐘，王光明，等.2021.2018年云南通海2次MS5.0地震震源機制解及深度測定［J］.地震研究，44（2）：133-144.Li J，Jiang J Z，Wang G M，et al.2021.Focal mechanism and focal depth determination of the two Tonghai，Yunnan，MS5.0 Earthquakes in 2018［J］.Journal of Seismological Research，44（2）：133-144.（in Chinese）

陸遠忠，宋俊高，戴維樂.1984.一個判斷震情的指標——震群的U值［J］.地震學(xué)報，6（S1）：495-508.Lu Y Z，Song J G，Dai W L.1984.U-value of earthquake swarms as an indicator for judging earthquake conditions［J］.Acta Seismologica Sinica，6（S1）：495-508.（in Chinese）

汪集旸，黃少鵬.1990.中國大陸地區(qū)大地?zé)崃鲾?shù)據(jù)匯編（第二版）［J］.地震地質(zhì)，12（4）：351-366.Wang J Y，Huang S P.1990.Compilation of geodetic heat flow data in Chinese mainland（second edition）［J］.Seismology and Ggeology，12（4）：351-366.（in Chinese）

王彬，鄧瑞生.2021.云南兩千年地震全記錄［M］.昆明：云南科技出版社.Wang B，Deng R S.2021.Complete record of earthquakes in Yunnan province for 2000 years［M］.Kunming：Yunnan Science and Technology Press.（in Chinese）

王光明，吳中海，劉昌偉，等.2022.2022 年 1 月 2 日寧蒗 MS5.5 地震序列重定位與發(fā)震構(gòu)造分析［J］.地震學(xué)報，44（4）：581-593.Wang G M，Wu Z H，Liu C W，et al.2022.Relocation and seismogenic structure analysis of the MS5.5 Ninglang earthquake sequence on January 2，2022［J］.Acta Seismologica Sinica，44（4）：581-593.（in Chinese）

王光明，歐陽鑫，劉自鳳，等.2024.維西—喬后斷裂帶及其周邊地區(qū)震群和前震序列活動特征［J］.地震學(xué)報，46（2）：192-207.Wang G M，Ouyang X，Liu Z F，et al.2024.Characteristics of earthquake swarms and foreshock sequences in Weixi-Qiaohou fault zone and its surrounding areas［J］.Acta Seismologica Sinica，46（2）：192-207.（in Chinese）

王清東.2015.利用雙差定位法研究云南地震分布的活動構(gòu)造意義［D］.武漢：武漢大學(xué).Wang Q D.2015.Using double-difference location method to research the active tectonic significance of Yunnan earthquakes distribution［D］.Wuhang：Wuhan University.（in Chinese）

王煒，楊德志.1987.利用Weibull分布研究華北地區(qū)前兆震群特征［J］.中國地震，3（4）：13-21.Wang W，Yang D Z.1987.A study on the features of precursory seismic swarms in North China by making use of Weibull distribution［J］.Earthquake Research in China，3（4）：13-21.（in Chinese）

肖陽，單斌，熊熊，等.2024.以2022年馬爾康震群為例研究震群發(fā)震機制［J］.地球物理學(xué)報，67（6）：2159-2174.Xiao Y，Shan B，Xiong X.et al.2024 Probing into the mechanism of earthquake swars：a case study of the 2022 Barkam earthquake swam［J］.Chinese Journal of Geophysics，67（6）：2159-2174.（in Chinese）

中國地震局監(jiān)測預(yù)報司.2020.地震危險性判定技術(shù)方法系列叢書—測震分析預(yù)測技術(shù)方法操作手冊［M］.北京：地震出版社.Department of Monitoring and Forecasting，China Earthquake Administration.2020.Earthquake hazard determination technical methods series—operation manual of seismic analysis and prediction technical methods［M］.Beijing：Seismological Press.（in Chinese）

朱傳鎮(zhèn)，傅昌洪，羅勝利.1981.震群與大地震關(guān)系的研究（華北地區(qū)）［J］.地震學(xué)報，3（2）：105-116.Zhu C Z，F(xiàn)u C H，Luo S L.1981.Study of earthquake swarms in relation to large earthquakes（North China area）［J］.Acta Seismologica Sinica，3（2）：105-116.（in Chinese）

朱傳鎮(zhèn)，王林瑛.1989.震群信息熵異常與地震預(yù)測［M］//許紹燮，陸遠忠，朱傳鎮(zhèn)，等.地震預(yù)報方法實用化研究文集（地震學(xué)專輯）.北京：學(xué)術(shù)書刊出版社，229-242.Zhu C Z，Wang L Y.1989.Abnormal entropy of earthquake swarm information and earthquake prediction［M］//Xu S X，Lu Y H，Zhu C Z，et al.Collection of practical research on earthquake prediction methods（Seismological Album）.Beijing：Academic Publishing House，229-242.（in Chinese）

Chen X，Shearer P M，Abercrombie R E.2012.Spatial migration of earthquakes within seismic clusters in Southern California：Evidence for fluid diffusion［J］.Journal of Geophysical Research：Solid Earth，117（B4）：1-30.

Dublanchet P，De Barros L.2021.Dual seismic migration velocities in seismic swarms［J］.Geophysical Research Letters，48（1）：e2020GL090025.

Jiang J Z，Li J，F(xiàn)u H.2019.Seismicity analysis of the 2016 MS5.0 Yunlong earthquake，Yunnan，China and its tectonic implications［J］.Pure and Applied Geophysics，176（3）：1225-1241.

Jones M，Wang B Q，Xu S X.1982.The foreshock sequence of the February 4，1975，Haicheng Earthquake（M=7.3）［J］.Journal of Geophysical Research，84（B6）：4575-4584.

Laske G，Masters G，Ma Z T，et al.2013.Update on CRUST1.0-A 1-degree global model of earth s crust［C］//Proceedings of EGU General Assembly Conference.Vienna：EGU.

Liu M，Tan Y J，Lei X L，et al.2024.Intersection between tectonic faults and magmatic systems promotes swarms with large-magnitude earthquakes around the Tengchong volcanic field，southeastern Tibetan Plateau［J］.Geology，52（4）：302-307.

Shapiro S A，Huenges E，Borm G.1997.Estimating the crust permeability from fluid-injection-induced seismic emission at the KTB site［J］.Geophysical Journal International，131（2）：15-18.

Shapiro S A，Patzig R，Rothert E，et al.2003.Triggering of seismicity by pore-pressure perturbations：permeability-related signatures of the phenomenon［J］.Pure and Applied Geophysics，160（5）：1051-1066.

Scholz C H.2019.The mechanics of earthquakes and faulting（3rd edition）［M］.New York：Cambridge University Press，512.

Shelly D R，Ellsworth W L，Hill D P.2016.Fluid-faulting evolution in high definition：Connecting fault structure and frequency-magnitude variations during the 2014 Long Valley Caldera，California，earthquake swarm［J］.Journal of Geophysical Research：Solid Earth，121（3）：1776-1795.

Shelly D R，Hill D P，Massin F，et al.2013.A fluid-driven earthquake swarm on the margin of the Yellowstone caldera［J］.Journal of Geophysical Research：Solid Earth，118（9）：4872-4886.

Talwani P，Acree S.1984.Pore pressure diffusion and the mechanism of reservoir-induced seismicity［J］.Pageoph，122（6）：947-965.

Waldhauser F.2000.A double-difference earthquake location algorithm：method and application to the Northern Hayward Fault，California［J］.Bulletin of the Seismological Society of America，90（6）：1353-1368.

Wetzler N，Brodsky E E，Chaves E J，et al.2023.Regional characteristics of observable foreshocks［J］.Bulletin of the Seismological Society of America，94（1）：428-442.

Yoshida K，Saito T，Urata Y，et al.2017.Temporal changes in stress drop，frictional strength，and earthquake size distribution in the 2011 Yamagata-Fukushima，NE Japan，earthquake swarm，caused by fluid migration［J］.Journal of Geophysical Research：Solid Earth，122（12）：10-379.

Zaccagnino D，Vallianatos F，Michas G，et al.2024.Are foreshocks fore-shocks？［J］.Journal of Geophysical Research：Solid Earth，129（2）：e2023JB027337.

Zhao L S，Helmberger D V.1994.Source estimation from broadband regional seismograms［J］.Bulletin of the Seismological Society of America，84（1）：91-104.

Zhou Y，Ren C，Ghosh A，et al.2022.Seismological characterization of the 2021 Yangbi foreshock-mainshock sequence，Yunnan，China：More than a triggered cascade［J］.Journal of Geophysical Research：Solid Earth，127（8）：e2022JB024534.

Zhu L P，Helmberger D V.1996.Advancement in source estimation techniques using broadband regional seismograms［J］.Bulletin of the Seismological Society of America 86（5）：1634-1641.

Zhu L P，Rivera L A.2002.A note on the dynamic and static displacements from a point source in multilayered media［J］.Geophys J Int，148（3）：619-627.

Research on the Characteristics of the MangshiML4.0 Earthquake Sequence in 2023

YANG Yuncun1，JIANG Jinzhong1，2，LIU Kexiang1，WANG Hualiu1，SUN Zigang1，YANG Mingchang1

（1.Yunnan Earthquake Agency，Kunming 650224，Yunnan，China）

（2.Kunming Institute of Prediction，China Earthquake Administration，Kunming 650224，Yunnan，China）

Abstract

Based on the broadband seismic waveform data and seismic phase reports produced by the Yunnan Seismic Network and the Tengchong Volcano Network，as well as the latest developed high-resolution common velocity model（Version 2.0）of the Crust in Sichuan-Yunnan region，the double-difference relative positioning method combined with waveform cross-correlation techniques is used to relocate 442 earthquakes in the Mangshi，Yunnan ML4.0 earthquake sequence in 2023.The CAP method is used to invert the focal mechanism solution of the ML≥3.5 earthquakes in the sequence，and the seismogenic mechanism of the sequence is studied in depth to explore whether this sequence is a foreshock sequence or just an ordinary earthquake swarm.Through the double-difference location method，precise locations of 363 earthquakes are obtained.Their epicenters are distributed in NE direction，and these earthquakes are distributed in a triangular shape.The optimal double-couple mechanism solution of the main shock，the ML4.0 earthquake，is obtained through CAP waveform inversion：the nodal Plane I strikes 45 °，the dip angle is 70 °，and the sliding angle is -61 °.The nodal Plane II has a strike of 167 °，an inclination of 35 °，a slip angle of -143 °，a moment magnitude（MW）of 3.91，and an optimal depth of 7 km.The temporal and spatial distribution characteristics of this earthquake sequence，the focal mechanism solution information，the travel time difference（S-P），the regional geological structure and the earthquake activity through the comprehensive precise positioning are used to explore the seismogenic structure and the seismogenic environment of this earthquake swarm.The results show that this earthquake swarm is not a foreshock sequence，but an ordinary earthquake swarm（magnitudes between 3 to 4）.This earthquake sequence occurred on at least two faults.The faults are in high stress state and form an angle.The underground fluid affected this earthquake sequence s evolution at different stages.

Keywords：the Mangshi earthquake swarm；double difference positioning；CAP inversion；earthquake activity；fluid diffusion