趙博，高原，馬延路

1 中國地震臺網(wǎng)中心,北京 100045 2 中國地震局地震預(yù)測研究所(地震預(yù)測重點實驗室),北京 100036

0 引言

根據(jù)中國地震臺網(wǎng)測定，北京時間2021年5月21日21時48分35秒，云南漾濞發(fā)生MS6.4地震，中國地震臺網(wǎng)中心(CENC)快報目錄中震源位置為(99.88°E，25.70°N，10 km)，主震前后數(shù)天前震和余震均較豐富.從5月18日18時49分ML3.8前震至5月24日8時止，中國地震臺網(wǎng)共監(jiān)測到地震500多次，其中，MS5.0以上地震3次，MS4.0～4.9地震12次，ML3.0～3.9地震13次.本次地震序列是前震-主震-余震型.

云南地區(qū)位于青藏高原東南緣，其西部為印-緬弧地區(qū)，其西北部是印度板塊與歐亞板塊的東部邊界，印度板塊在此處東向俯沖(Turner et al.,2013)，也被稱為印度板塊與青藏高原的側(cè)向碰撞帶.由于印度板塊與歐亞板塊的碰撞，導(dǎo)致青藏高原巖石圈物質(zhì)向東逃逸(Tapponnier et al.，2001)，在受到四川盆地阻擋后轉(zhuǎn)向SE向，進而向南，造成川滇菱形塊體發(fā)生順時針旋轉(zhuǎn)(張培震等，2003；Zhang et al.,2004)，地表變形與地殼應(yīng)力場都呈現(xiàn)了扇形散開的運動形態(tài)(Jin et al.,2019;Tian et al.,2019).地震層析成像研究結(jié)果發(fā)現(xiàn)青藏高原東南緣深部速度結(jié)構(gòu)異?，F(xiàn)象(Lei and Zhao,2016;劉偉等，2019)，區(qū)域分布的豐富金屬礦床和強烈的地震活動都表明，該區(qū)域存在劇烈的深部物質(zhì)運動(王瓊和高原，2014；王瓊等,2015；李秋生等，2018；高原等，2020)，地震各向異性揭示了復(fù)雜的區(qū)域性上地殼介質(zhì)變形圖像(Shi et al.,2012).漾濞MS6.4地震位于塊體交界區(qū)，周圍分布有怒江、瀾滄江、中甸、紅河等多條大型地震斷裂(圖1a).距離這次地震最近的斷裂是維西—喬后斷裂，位于該地震序列的NE側(cè)(龍鋒等，2021)，是蘭坪—思茅褶皺帶與昌都—云嶺褶皺帶的分界斷層(常祖峰等，2016)，其南部與紅河斷裂相連，北部與中甸斷裂相連，構(gòu)成了川滇塊體的西部邊界.地表形變觀測研究顯示，川滇地塊西南邊界斷裂滑動速率較低，但川滇地塊整體處于較強的構(gòu)造變形環(huán)境，且地震活動呈現(xiàn)彌散分布特征(Allen et al.,1984;Li et al.,2020;張克亮等，2021).1970年以來，本次漾濞MS6.4主震震中50 km以內(nèi)共發(fā)生過ML3.0以上地震140多次，圖1a給出了2009年1月至2021年4月30日研究區(qū)ML≥2.0地震震中分布.大地電磁研究顯示，漾濞震源區(qū)存在顯著電性橫向不均勻性，可能與區(qū)域構(gòu)造應(yīng)力作用下巖體相互作用有關(guān)(葉濤等，2021).

準(zhǔn)確的震源位置，特別是震源深度，以及震源機制和區(qū)域構(gòu)造應(yīng)力場對理解發(fā)震構(gòu)造、孕震機制有重要意義(許忠淮等，1989；高原等，1997；崔效峰等，2006；Zhao et al.,2013;易桂喜等，2015，2017;Xu et al.,2016;Luo et al.,2016;孫亞君等，2017).本研究首先采用雙差定位法(Waldhauser and Ellsworth，2000)對主震前、后數(shù)天(2021年5月18—24日)的地震序列進行重新定位，并利用面波振幅譜陷波特征對主震震源深度進行準(zhǔn)確定位.然后對ML≥3.5地震進行震源機制反演，最后利用震源機制解結(jié)果反演震源區(qū)主壓應(yīng)力方向.

1 數(shù)據(jù)

CENC產(chǎn)出的地震目錄中，2021年5月18—24日共記錄到495個地震，本文收集了19982條P波到時數(shù)據(jù)，13826條S波到時數(shù)據(jù)用于雙差定位.雙差定位所用到的65個國家地震臺網(wǎng)臺站和云南區(qū)域地震臺網(wǎng)臺站(圖1a)，定位使用的一維速度模型見圖1b(吳建平等，2004)，圖1c展示了利用面波振幅譜確定主震震源深度時使用的遠(yuǎn)場臺站.面波頻譜法測定震源深度時，臺站的選擇只與地震的震源機制解有關(guān)，根據(jù)Aki和Richards(2002)提出的面波遠(yuǎn)場輻射圖型公式以及Tsai和Aki(1979)對面波振幅譜與震源深度的位置關(guān)系表明，不同的震源機制解參數(shù)唯一確定了面波振幅譜陷波相出現(xiàn)的方位臺站，震源深度的變化不改變陷波相出現(xiàn)的臺站方位.大量的實測波形數(shù)據(jù)和理論波形數(shù)據(jù)(Fox et al.,2012;Heyburn et al.,2013;趙博等，2018,2019)也表明，當(dāng)震源機制解確定后，面波陷波相對應(yīng)頻率只隨深度而變化，方位角覆蓋不影響震源深度測定，這也是該方法不同于傳統(tǒng)的走時測定深度的特點和優(yōu)點.為了使臺站記錄的面波信號發(fā)育完整且波形質(zhì)量較高，選擇震中距大于500 km的國家臺站.對于MS6.4主震及部分前震和余震震源機制解，我們利用區(qū)域臺站記錄的波形進行反演(參見圖1a).

2 漾濞MS6.4地震序列定位

震源位置的準(zhǔn)確測定是地震學(xué)的一個基礎(chǔ)性問題.基于走時方程的蓋革方法(Geiger,1912)提出后，地震學(xué)家提出和發(fā)展了多種地震定位方法.Richards和Shearer(2000)發(fā)展了SSST方法(Source-Specific Station Terms method)，該方法適用于大區(qū)域地震定位，特別是存在橫向速度不均勻介質(zhì)時，添加了臺站-震源修正項，以減小臺站到時的固有誤差(Lin et al.,2007).許力生等(2013a,b)發(fā)展的逆時成像技術(shù)直接利用波形進行定位，避免了線性化過程的誤差，提高了定位精度.Waldhauser和Ellsworth(2000)發(fā)展了雙差定位算法，通過地震對在同一臺站的走時雙差，消除了地震波傳播路徑上的速度擾動，最大程度地減小了速度模型對定位的影響.雙差定位法是一種相對定位方法，在國內(nèi)外得到廣泛應(yīng)用.例如，國外學(xué)者利用雙差定位法對美國加州地區(qū)的地震目錄進行了重新定位(Waldhauser and Ellsworth,2000;Schaff and Waldhauser,2005;Waldhauser and Schaff,2008；Waldhauser,2009)，獲得該區(qū)域準(zhǔn)確的地震目錄.國內(nèi)學(xué)者利用該方法對川滇地區(qū)以及中國大陸其他地區(qū)地震進行了重新定位，獲得了汶川地震、玉樹地震、蘆山地震和九寨溝地震等序列更可靠的震源位置(陳九輝等，2009；吳建平等，2009；房立華等，2011；Zhao et al.,2011,2012；易桂喜，2012，2016，2017；陳晨和胥頤，2013；張廣偉和雷建設(shè)，2013;趙博等，2013).

本文對CENC地震目錄中2021年5月18—24日的495個(ML1.5以上)地震進行了重新定位.雙差定位時共形成2410個地震對，用到P波和S波震相到時33808個，其中P波到時19982個，S波到時13826個.由于P波在拾取時較S波精確，我們賦予P波到時在定位中的權(quán)值為1.0，S波到時權(quán)值為0.7；定位中設(shè)定待定地震對之間的最大距離為5 km，且有至少8個到時差數(shù)據(jù).采用最小二乘正交三角分解法(LSQR)進行10次迭代，一共定位出331個重新定位結(jié)果(圖2).本研究雙差定位應(yīng)用的速度模型參見圖1b(吳建平等，2004).

圖1 研究區(qū)域構(gòu)造背景、初始震中位置、地震臺站分布及速度模型(a)黑色粗線為塊體邊界(張培震等，2003)，黑色細(xì)線為地震斷層，黑色三角形為雙差定位和震源機制反演用到的地震臺站，黑色圓圈為歷史地震，彩色圓點為本次地震序列分布(CENC地震目錄).F1 怒江斷裂，F(xiàn)2 瀾滄江斷裂，F(xiàn)3 中甸斷裂，F(xiàn)4 紅河斷裂，F(xiàn)5 程海斷裂，F(xiàn)6 維西—喬后斷裂，F(xiàn)7 紅河斷裂北段.(b)雙差定位用到的一維速度模型(吳建平等，2004)，藍色為P波速度，紅色為S波速度.(c)紅色三角形為面波振幅譜定位震源深度用到的國家臺站.Fig.1 Tectonic settings of study region,initial distribution of earthquake locations and seismic stations,and the used velocity model(a)The black thick lines denote boundaries of tectonic blocks (Zhang et al.,2003)and the black thin lines denote the seismic faults.The black triangles denote the seismic stations used in hypoDD location and focal mechanisms inversion.The black circles denote the history earthquakes and color dots denote Yangbi MS6.4 earthquake sequence (earthquake locations is from CENC catalog).F1 Nujiang fault,F2 Lancangjiang fault,F3 Zhongdian fault,F4 Honghe fault,F5 Chenghai fault,F6 Weixi-Qiaohou fault,F7 the north part of Honghe fault.(b)1D velocity model used in hypoDD.The blue line indicates the velocity of P phase,and the red line indicates the velocity of S phase.(c)The red triangles denote the seismic stations used in source depth determined by surface wave amplitude spectrum location.

定位結(jié)果顯示，漾濞MS6.4地震序列呈NW-SE向展布(圖2a)，主震位置為(99.869°E，25.689°N，8.8 km)，位于余震序列條帶的NW端(圖2a五角星位置)，余震序列向SE方向延展.地震序列條帶的NE側(cè)是維西—喬后斷裂，其走向與地震序列震中分布一致，但其地表位置與震中分布相距較遠(yuǎn)，約為15 km.根據(jù)本文的定位結(jié)果，可以推斷維西—喬后斷裂并非這次漾濞MS6.4地震的發(fā)震斷裂.地震發(fā)生后，中國地震局組織了這次地震的科學(xué)考察，初步結(jié)論認(rèn)為這次漾濞MS6.4地震發(fā)生在維西—喬后斷裂SW側(cè)的一條次級斷裂上，與本文定位結(jié)果基本一致.

漾濞地震序列重新定位后的震中位置展現(xiàn)出清晰的條帶分布特點(圖2a)，為分析其垂向深度分布，本文做了兩個深度剖面(如圖2).A1—A2剖面顯示(圖2b)，地震序列條帶沿NW-SE向分布長約20 km.NE-SW方向的剖面B1—B2顯示(圖2c)，地震序列近于垂直分布(圖2c).與重新定位前(圖2d、2e)相比，重新定位后深度剖面有明顯的改善.定位前、后深度剖面和深度統(tǒng)計直方圖(圖3)顯示，本次地震序列震源深度主要分布在5～15 km之間.初始震源深度(CENC目錄)在9 km和11 km處有兩個峰值(圖3a)，各占地震總數(shù)的20%左右.重新定位后，震源深度主要集中在8～13 km，其中10～12 km深度的地震約占地震總數(shù)的60%左右(圖3b).雙差定位結(jié)果的深度剖面顯示，主震深度為8.8 km，在地震序列的NW端，大部分余震位于主震下方，且深度≤5 km的地震較少，本研究推測，主震從起始破裂點(NW端)向SE方向破裂，且有向深部破裂的趨勢，發(fā)震斷層可能在地表沒有出露，這與初期地質(zhì)調(diào)查結(jié)果相符合.

圖2 漾濞地震序列分布(a)重新定位后震中的水平分布.A1—A2 和 B1—B2為剖面位置.(b)和(c)是重新定位后的兩個深度剖面.(d)和(e)是重新定位前的兩個深度剖面.Fig.2 The distribution of Yangbi earthquake sequence(a)The horizontal distribution of epicenters after relocated.A1—A2 and B1—B2 are the locations of the cross sections.(b)and (c)are the depth profiles after relocation.(d)and (e)are the depth profiles before relocation.

圖3 重新定位前后深度統(tǒng)計直方圖(a)重新定位前的深度統(tǒng)計直方圖;(b)重新定位后的深度統(tǒng)計直方圖.Fig.3 The depth statistical histograms before and after relocation(a)and (b)are the depth statistical histograms before and after relocation,respectively.

3 主震震源深度確定

對于一個地震序列，主震震源深度的準(zhǔn)確定位，對理解發(fā)震構(gòu)造、孕震層深度等有重要意義.深源地震遠(yuǎn)場記錄波形的深震相pP和sP對震源深度有較強的控制，但對于地殼內(nèi)發(fā)生的淺源地震，很難識別這種上行轉(zhuǎn)換震相，因此對于深度的準(zhǔn)確確定比較困難.研究發(fā)現(xiàn)(Tsai and Aki，1970)，地震面波的振幅譜特征與震源深度有密切關(guān)系，振幅譜的陷波相(notch phase)出現(xiàn)的位置(即頻譜中notch phase對應(yīng)的頻率)隨震源深度變化敏感.Tsai和Aki(1970)分析了大量震例，認(rèn)為notch phase的位置與破裂尺度和介質(zhì)屬性等關(guān)系不大，只依賴于震源機制及震源深度.當(dāng)震源機制確定后，可以利用notch phase與震源深度的唯一依賴關(guān)系準(zhǔn)確確定震源深度.研究者們利用該原理對1999年10月22日中國臺灣集集地震MW5.9余震、2002—2009年伊朗發(fā)生的4次中等震級地震、2004年9月28日美國ParkfieldMW6.0地震、2013年蘆山MS7.0地震、2017年九寨溝MS7.0地震以及四川地區(qū)近十年中強地震做了震源深度準(zhǔn)確定位(Fox et al.,2012;Heyburn et al.,2013;趙博等，2018,2019).

Aki和Richards(2002)給出的Rayleigh面波遠(yuǎn)場輻射能量公式為

(1)

其中，ωn為角頻率，φ為震中-臺站方位角，kn為波數(shù)，h為震源深度，Mij為地震矩張量.Rayleigh面波振幅A(φ,ωn)是觀測方位與頻率的函數(shù)，受震源機制和震源深度的影響.為了快速找到遠(yuǎn)場Rayleigh面波notch phase出現(xiàn)的方位角，我們首先利用主震震源機制解(震源機制解反演見第4節(jié))，假定震源深度在10 km時，計算0°～360°范圍內(nèi)理論地震圖振幅譜(圖4).理論地震圖振幅頻譜特征顯示，方位角在80°～110°和170°～290°范圍(圖4中黑色橢圓內(nèi))面波振幅譜出現(xiàn)陷波特征，即面波能量在陷波頻率處較小，而在陷波頻率兩側(cè)能量較高.由于陷波出現(xiàn)的觀測方位不隨深度變化，并考慮到面波在遠(yuǎn)場發(fā)育完全，選擇震中距大于500 km、方位角位于80°～110°和170°～290°方位、地震記錄質(zhì)量較高的國家臺站的記錄進行深度測定.本文使用的臺站位置見圖1c.

地震面波在傳播過程中會出現(xiàn)頻散現(xiàn)象，即不同頻率的面波傳播速度不同.以福建的國家臺NPDK為例，對MS6.4主震波形進行多個頻帶的帶通濾波(圖5a)，可以發(fā)現(xiàn)不同頻率的面波出現(xiàn)的時間不同.為了準(zhǔn)確測量面波能量，并且減小其他長周期能量成分(如長周期噪聲等)，對觀測波形進行相位匹配濾波從而達到去頻散效果(圖5b)，將0.01～0.1 Hz的面波能量壓縮至窗口中心位置，并在窗口中心加500 s的漢寧窗，去除非面波能量.

圖5 觀測波形面波頻散和相位匹配濾波(a)觀測波形的面波頻散現(xiàn)象，其中每條波形右側(cè)數(shù)字為濾波頻帶;(b)對原始波形進行相位匹配濾波，從上至下分別為原始波形、相位匹配濾波后的去頻散波形、加500 s漢寧窗的去頻散波形.O為發(fā)震時刻，A為P波到時，T0為S波到時，T4為時間窗中心位置.圖中波形的記錄臺站為NPDK臺(國家臺).Fig.5 The dispersion of observed wave and phase matching filter(a)The observation waves are bandpass filtered.The filter bandwidths are on the right.(b)Phase-matched filtering of the original waveform.From top to bottom,the waves are original observed wave,dedispersion waveform after phase matching filter,and dedispersion waveform with 500 s Hanning window.O marks the original time.A and T0 mark the arrival times of P and S waves,respectively.T4 marks the center of time window.The recording station of the seismic waveform in the figure is the NPDK station (national station).

對去除頻散效應(yīng)并壓縮面波能量至?xí)r間窗中心位置的波形進行快速傅里葉變換，獲得觀測波形的振幅頻譜.然后利用震源機制解計算每個測量臺站處不同深度的理論地震圖，并以與觀測波形同樣的處理方法進行去頻散和時間壓縮以及快速傅里葉變換.圖6為四個臺站的理論波形振幅譜和觀測波形振幅譜對比.紅色觀測波形振幅譜出現(xiàn)明顯能量下降notch phase，通過計算1～30 km深度的理論波形振幅譜與觀測波形振幅譜notch phase出現(xiàn)的位置，可獲得每個臺站的最優(yōu)擬合深度(表1).取所有測量值的均值作為最終的深度測量結(jié)果，漾濞MS6.4主震震源深度為7.5 km，均方差為0.7 km.

圖6 理論波形與觀測波形振幅譜(a)—(d)分別為GZH,NPDK,SHG,XNY臺站的理論波形(藍色和黑色)和觀測波形(紅色)的振幅譜.臺站名位于左上方.理論振幅譜曲線旁邊的數(shù)字為理論震源深度，藍色波形為最優(yōu)擬合深度的振幅譜.Fig.6 The amplitude spectra of synthetic surface wave and observed wave(a)—(d)are the amplitude spectra of synthetic surface wave (black and blue curves)and observed wave (red curves)at GZH,NPDK,SHG,XNY seismic stations.The station names are at the top left.The number near each synthetic curve is the synthetic focal depth.The blue curves indicate the amplitude spectra of the best-fitting.

4 震源機制與應(yīng)力場

本文利用CPS程序(Herrmann,2013)和近場波形數(shù)據(jù)進行波形擬合，反演了主震和較大余震(ML>3.5)的震源機制解.該方法首先利用區(qū)域一維速度模型(圖1b)計算格林函數(shù)，然后進行體波波形擬合，從而反演震源機制解.我們用到的觀測波形為圖1a中顯示的云南地區(qū)的國家臺和區(qū)域臺站數(shù)據(jù).圖7展示了主震最優(yōu)擬合情形下的理論波形與觀測波形對比.漾濞MS6.4地震的震源機制解為：節(jié)面I走向143°/傾角80°/滑動角-165°，節(jié)面II走向50°/傾角75°/滑動角-10°，矩震級MW6.0.根據(jù)地震序列NW-SE向分布特征，可以推斷節(jié)面I為發(fā)震斷層面，斷層面走向與維西—喬后斷裂平行，發(fā)震斷裂可能是維西—喬后斷裂的次生斷裂.主震震源機制節(jié)面傾角為80°，結(jié)合地震序列震源深度的近垂直分布特征，推斷該次生斷層為高傾角，SW傾向的斷裂.震源機制P軸方位角為近NS向(7°)，傾角近水平(18°).震源區(qū)近NS向的主壓應(yīng)力使發(fā)震斷層發(fā)生了右旋走滑錯動.圖8顯示了主震震源機制參數(shù)與震源深度之間的關(guān)系.主震震源機制反演得到的深度為7 km，該深度為矩心深度，反映了斷層釋放最大能量的位置.前文雙差定位是利用初至震相到時確定的震源深度為起始破裂點深度(8.8 km).二者相差約1.8 km，主震矩心深度略淺于起始破裂深度.

圖7 漾濞MS6.4主震觀測波形和理論波形對比圖圖中紅色為觀測波形，藍色為最優(yōu)擬合的理論波形，臺站名、方位角和震中距在波形右側(cè)；上方Z,R,T分別表示垂直向、徑向和切向分量.Fig.7 Comparison of observed and synthetic waveforms of Yangbi MS6.4 main shockThe red are the observed waveforms,and blue are the best fit synthetic waveforms.The names of seismic stations,azimuth and epicentral distance are on the right of the waves.On the top,Z,R and T indicate vertical,radial and tangential direction components,respectively.

圖8 漾濞MS6.4主震擬合度-深度分布圖Fig.8 The distribution of fitness-depth of Yangbi MS6.4 main shock

表1 各臺站深度測量結(jié)果Table 1 The result of the depth location at each station

本研究獲得主震及29個余震的震源機制解(圖9a和表2).Frohlich(1992)利用震源機制解應(yīng)力P軸、T軸、B軸的傾角大小進行地震類型分類，并用三角形圖解法來描述不同類型.本研究按照三角形圖解法對漾濞地震序列震源機制進行分類(圖9c).結(jié)果顯示，30個震源機制解中，大部分為走滑型兼有少量正斷分量.其中，8號、14號、15號和26號地震P軸傾角均大于40°，正斷分量較大.

表2 主震、前震和余震震源機制解Table 2 Focal mechanisms of the mainshock,foreshocks and aftershocks

續(xù)表2

震源機制解資料可以用來反演區(qū)域構(gòu)造應(yīng)力場.雖然一次地震的P軸只反映了本次破裂釋放壓應(yīng)力的方向，它與區(qū)域構(gòu)造壓應(yīng)力方向存在一定差異，但對研究地區(qū)內(nèi)較多歷史地震的P軸和T軸方位做統(tǒng)計平均，可以得出構(gòu)造應(yīng)力的方向.本文采用力軸張量計算法(鐘繼茂和程萬正，2006；Zhao et al.,2012;趙博等，2013,2019)，對漾濞地震震源區(qū)構(gòu)造應(yīng)力方向進行反演.圖9b為本次地震序列主震和余震震源機制解P軸(黑色玫瑰圖)和T軸(紅色玫瑰圖)方位角，P軸為近NS向，T軸為近EW向.由P軸、B軸和T軸的方位角和傾角很容易得到每個震源機制解在力軸坐標(biāo)系下的力軸張量，經(jīng)過坐標(biāo)轉(zhuǎn)換可以得到應(yīng)力張量在地理坐標(biāo)系下的值并可求出平均力軸張量，然后通過求解平均力軸張量的主軸方位即可獲得區(qū)域構(gòu)造應(yīng)力方向.經(jīng)過計算，本研究得到漾濞地震震源區(qū)主壓應(yīng)力方位為7°，近NS向，傾角為3°(圖9a).

近場快剪切波得到的上地殼各向異性研究結(jié)果表明，快剪切波偏振優(yōu)勢方向平行于主壓應(yīng)力方向(Crampin and Peacock,2005;Gao et al.,2011;Li et al.,2021).震源機制反演和地殼各向異性研究結(jié)果表明，應(yīng)力方向與上地殼快剪切波偏振方向有很好的一致性(Zhao et al.,2012).川滇地區(qū)的快剪切波偏振優(yōu)勢方向研究(高原等，2020)，揭示出區(qū)域主壓應(yīng)力方向的分布.張藝和高原(2017)的剪切波分裂研究結(jié)果則顯示，在這次漾濞地震震源區(qū)所在的局部區(qū)域，快剪切波偏振優(yōu)勢方向為NNE近NS方向，與本文得到的漾濞MS6.4地震震源區(qū)主壓應(yīng)力方向為7°的結(jié)果高度一致，表明區(qū)域上地殼剪切波分裂快波偏振方向可用于區(qū)域構(gòu)造應(yīng)力場的分析.Gao 等(2019)的結(jié)果揭示，這次漾濞MS6.4地震震源區(qū)位于上地殼各向異性快波方向發(fā)生變化的區(qū)域附近.

5 討論與結(jié)論

通過漾濞MS6.4地震序列的重新定位結(jié)果，揭示了本次地震的發(fā)震斷層產(chǎn)狀結(jié)構(gòu)為傾向SW的高傾角右旋走滑斷層.主震發(fā)生在淺部地殼，地震序列深度在15 km以上.雙差走時定位確定的深度為8.8 km，利用面波振幅譜陷波相定位的深度為7.5 km.本次地震震源區(qū)周邊臺站較密集，臺網(wǎng)空隙角較小，走時定位結(jié)果與振幅譜深度精細(xì)分析結(jié)果相差不大.區(qū)域應(yīng)力場表征這局部地震的斷層模式，因此震源機制解是研究區(qū)域應(yīng)力場的有效手段，結(jié)合地殼剪切波分裂數(shù)據(jù)，可以真實的反映區(qū)域的構(gòu)造應(yīng)力場.

(1)經(jīng)過雙差重新定位獲得了主震震源位置為(99.869°E，25.689°N，8.8 km).主震位于余震帶NW端，地震序列呈NW-SE向展布，長約20 km，根據(jù)主震與余震序列的位置關(guān)系，推測本次地震為單側(cè)破裂.震源深度分布在5～15 km，對重新定位后的震源深度進行統(tǒng)計(圖3b)，深度在8～13 km范圍內(nèi)的地震約占地震總數(shù)的80%左右.垂向的深度剖面顯示，地震在深度方向上近垂直分布，與高傾角的震源機制結(jié)果相一致.定位結(jié)果顯示，余震分布與其NE側(cè)近平行的維西—喬后斷裂之間有約15 km的距離，因此推斷維西—喬后斷裂不是本次地震的發(fā)震斷層，其發(fā)震斷層可能是維西—喬后斷裂SW側(cè)的一條次級斷裂.

(2)漾濞MS6.4主震震源深度較淺，屬于發(fā)生于上地殼的淺源地震.雖然雙差定位算法在一定程度上減小了速度擾動引起的定位誤差，但是所選用模型的分層結(jié)構(gòu)以及地殼淺部的低速層厚度都會在深度方向引起較大誤差.我們利用遠(yuǎn)震面波振幅譜陷波特征進行震源深度精確定位，避免了區(qū)域速度結(jié)構(gòu)不準(zhǔn)確性對震源深度定位的影響.經(jīng)過測量十余個波形質(zhì)量較高的國家臺站的面波振幅譜notch phase的位置，較精確的獲得了主震震源深度為7.5 km.該方法由于測量的是遠(yuǎn)場面波能量，能可靠地測量中強地震震源深度.

(3)漾濞MS6.4主震震源機制解顯示為一次右旋走滑型地震，矩震級為MW6.0.29個較強前震和余震類型為走滑型地震，兼少量正斷分量；其中有4個地震正斷分量較大，P軸傾角大于40°.結(jié)合定位結(jié)果的分布情況，判定主震破裂面為：走向143°、傾角80°、滑動角-165°.通過主震和較大前震、余震震源機制解的力軸信息，反演得到震源區(qū)的主壓應(yīng)力方向為近NS向(7°).此外，地殼剪切波分裂也可以反映區(qū)域主壓應(yīng)力情況(Crampin and Peacock,2005；Gao et al.,2011,2019).吳鵬等(2020)對三江地區(qū)進行了各向異性研究，得到了該區(qū)域內(nèi)的TUS臺剪切波分裂參數(shù)(圖9a)，其快剪切波偏振方向為148°，與主震震源機制解走向一致，也與該臺附近的維西—喬后斷裂、紅河斷裂北段走向一致，表明TUS臺的快剪切波偏振方向受到附近斷層影響.但張藝和高原(2017)以及高原等(2020)等利用附近區(qū)域內(nèi)更多的臺站剪切波分裂資料進行綜合分析，得到區(qū)域的快剪切波偏振方向為近NS向(圖9a中的藍色玫瑰圖)，該結(jié)果反映了區(qū)域主壓應(yīng)力方向，與本文反演的震源區(qū)主壓應(yīng)力方向高度一致.雖然剪切波分裂與震源機制反演應(yīng)力場的物理機制不同，但兩種結(jié)果都較為準(zhǔn)確的反映了區(qū)域主壓應(yīng)力的方向.

圖9 主震及余震震源機制解及震源區(qū)主壓應(yīng)力方向(a)中紅色沙灘球為主震震源機制，藍色沙灘球為兩次5級地震震源機制，綠色沙灘球為MS3.5～4.9地震震源機制.紅色短線為震源機制解P軸方位，藍色箭頭為本研究反演的主壓應(yīng)力方向，藍色玫瑰圖為剪切波分裂快波偏振方向(張藝和高原，2017；高原等，2020)，TUS臺站藍色短線為該臺剪切波分裂偏振方向(吳鵬等，2020).(b)為所有震源機制解P軸和T軸方向，黑色為P軸，紅色為T軸.(c)震源機制解分類.Fig.9 Focal mechanisms of main shock and aftershocks,and the direction of the main compressive stress in focal area(a)Red beach ball is the focal mechanism of main shock,the blue beach balls are the focal mechanism of two aftershocks with 5.0≤MS≤5.9,and green beach balls are the focal mechanism of aftershocks with 3.5≤MS≤4.9.The red bars indicate the direction of P axis of focal mechanism solutions.The blue arrow denotes the direction of the main compressive stress in this study.The blue rose diagram indicates the fast shear-wave polarization in this region (Zhang and Gao,2017;Gao et al.,2020).The blue bar indicates the fast shear-wave polarization at TUS station (Wu et al.,2020).(b)The directions of P (black)and T (red)axis of all focal mechanisms.(c)The classification of focal mechanisms.