強(qiáng)生銀， 王宏偉， 溫瑞智， 李春果， 任葉飛

中國(guó)地震局工程力學(xué)研究所， 中國(guó)地震局地震工程與工程震動(dòng)重點(diǎn)實(shí)驗(yàn)室， 哈爾濱 150080

0 引言

據(jù)中國(guó)地震臺(tái)網(wǎng)測(cè)定，自2021年5月18日以來(lái)云南省大理白族自治州漾濞彝族自治縣連續(xù)發(fā)生多次MS3.0～4.9地震，并在5月21日21∶21∶25發(fā)生MS5.6地震后，緊接著于21∶48∶34在其西北約6.7 km處發(fā)生MS6.4地震，此次地震是漾濞地震序列的主震，震中位于北緯25.67°、東經(jīng)99.87°，距離漾濞縣城約9.0 km，震源深度8 km，主震震中處于此次地震序列的最西北端.據(jù)Global CMT震源機(jī)制解可知(https:∥www.globalcmt.org/)，漾濞MS6.4地震為走滑型地震，矩震級(jí)MW=6.1.截至6月1日已發(fā)生MS≥3.0前震和余震共46次，其中MS5.0～5.9地震3次、MS4.0～4.9地震13次，最大震級(jí)余震為5月21日22∶31∶10發(fā)生的MS5.2地震，大約位于此次地震序列最東南端，距離主震震中約13.5 km.此次地震序列的震中呈北西-南東向分布(見(jiàn)圖1b)，傾向南西，破裂尺度約為20 km，集中分布在6～15 km深度內(nèi).漾濞地震序列是典型的前震-主震-余震型地震序列.

通過(guò)震后現(xiàn)場(chǎng)震害調(diào)查，參照地震構(gòu)造背景、余震分布、震源機(jī)制解、儀器烈度，并結(jié)合強(qiáng)震動(dòng)觀測(cè)數(shù)據(jù)，云南省地震局發(fā)布了云南漾濞MS6.4地震宏觀烈度圖(http:∥www.yndzj.gov.cn/yndzj/_300559/_300651/629959/index.html)，此次地震最高烈度為Ⅷ度，等震線長(zhǎng)軸呈北西-南東走向.截至5月22日，地震共造成3人死亡，28人受傷，192間房屋倒塌，13090間房屋一般損壞(http:∥www.dali.gov.cn/dlrmzf/c105806/202105/8d6fcefcca214dd6b3933a49ef0eb9f0.shtml).

初步研究表明漾濞MS6.4地震發(fā)生于紅河斷裂北段的維西—喬后斷裂西側(cè)的一條北西-南東向右旋走滑次生斷裂(https:∥www.eq-igl.ac.cn/zhxw/info/2021/33882.html，見(jiàn)圖1a)，紅河斷裂是青藏高原東南部的川滇菱形塊體的西南側(cè)邊界.該區(qū)域是印度板塊與歐亞板塊碰撞擠壓的邊界地帶,川滇菱形塊體是青藏高原東南緣擠出最強(qiáng)烈的活動(dòng)塊體,錯(cuò)綜復(fù)雜的地質(zhì)構(gòu)造和強(qiáng)烈的斷裂帶活動(dòng)使之成為我國(guó)地震活動(dòng)最頻繁的地區(qū)之一(闞榮舉等,1977; 汪一鵬等, 2003; 皇甫崗等, 2007).1925年大理7.0級(jí)地震和1996年麗江7.0級(jí)地震是漾濞MS6.4地震震中附近發(fā)生的有記錄的最大震級(jí)歷史地震(見(jiàn)圖1a)，近年來(lái)震中附近中強(qiáng)地震活動(dòng)十分強(qiáng)烈(例如：2013年洱源MS5.0地震、2016年云龍MS5.0地震、2017年漾濞MS5.1地震).

根據(jù)中國(guó)地震局工程力學(xué)研究所提供的強(qiáng)震動(dòng)觀測(cè)數(shù)據(jù)，漾濞MS6.4地震共觸發(fā)28個(gè)自由場(chǎng)強(qiáng)震動(dòng)觀測(cè)臺(tái)站，除053YPX臺(tái)站只獲取了兩個(gè)水平分量的記錄外，其余臺(tái)站均獲得了完整的三分量強(qiáng)震動(dòng)加速度記錄.此次地震觸發(fā)臺(tái)站分布如圖1a所示，距震中100 km內(nèi)僅觸發(fā)了5個(gè)臺(tái)站(053YBX、053DLY、053YPX、053BTH、053BCJ)，距離震中最近的053YBX臺(tái)站(震中距7.9 km)獲得了此次地震的最大峰值地面加速度(PGA)記錄，其東西(EW)、南北(NS)、豎向(UD)的PGA分別為379.9 cm·s-2、720.3 cm·s-2和448.4 cm·s-2(未處理原始記錄).本文依次采用基線校正、記錄波形首尾增加余弦窗并補(bǔ)零、巴特沃斯非因果帶通濾波處理強(qiáng)震動(dòng)記錄，濾波帶寬為0.08～30 Hz，處理后記錄有效頻帶為0.1～24 Hz.

圖1 漾濞地震序列震中位置、漾濞MS6.4主震觸發(fā)強(qiáng)震動(dòng)臺(tái)站分布(a)及其宏觀地震烈度(b)Fig.1 Epicenters of the Yangbi earthquake sequence, the triggered strong-motion observation stations during the MS6.4 mainshock (a), and the macroseismic intensity map (b)

漾濞MS6.4地震獲取了較多的強(qiáng)震動(dòng)觀測(cè)記錄，但多數(shù)為遠(yuǎn)場(chǎng)記錄，為預(yù)測(cè)此次地震的地震動(dòng)及其可能破壞情況，本文采用隨機(jī)有限斷層三維模擬方法，并結(jié)合強(qiáng)震動(dòng)觀測(cè)記錄，模擬了震中附近區(qū)域2823個(gè)虛擬觀測(cè)點(diǎn)的三分量地震動(dòng)加速度時(shí)程，基于模擬記錄預(yù)測(cè)了此次地震的烈度分布.

1 隨機(jī)有限斷層三維地震動(dòng)模擬方法

(1a)

(1b)

(1c)

2 漾濞MS6.4地震的地震動(dòng)模擬

2.1 隨機(jī)有限斷層模擬輸入?yún)?shù)

根據(jù)Global CMT提供的震源機(jī)制解，漾濞MS6.4地震矩震級(jí)MW=6.1，走向角、傾角和滑動(dòng)角分別為315°、86°和168°，根據(jù)地震破裂面幾何(沿走向的長(zhǎng)度L和沿傾向的寬度W)與MW、破裂類型的經(jīng)驗(yàn)關(guān)系(Wells and Coppersmith, 1994)及地殼孕震層對(duì)破裂面寬度的限制(Pacheco et al., 1992; Shaw and Wesnousky, 2008)，估計(jì)漾濞MS6.4地震震源破裂面的長(zhǎng)度L=16.5 km、寬度W=7.7 km，根據(jù)Somerville等(1999)建立的平均滑動(dòng)位錯(cuò)(D)與矩震級(jí)經(jīng)驗(yàn)關(guān)系，設(shè)定D=38.74 cm，將破裂面沿走向和沿傾向劃分為1.1 km×1.1 km共105個(gè)子斷層.根據(jù)主余震空間分布(見(jiàn)圖1b)，漾濞MS6.4地震近似為自西北向東南方向的單側(cè)破裂事件，起始破裂點(diǎn)位置沿走向坐標(biāo)為1.0L，Mai等(2005)對(duì)不同破裂類型地震的起始破裂點(diǎn)位置的統(tǒng)計(jì)研究表明，走滑地震和地殼內(nèi)傾滑地震的破裂面上起始破裂點(diǎn)沿傾向的位置近似服從Weibull分布，沿傾向?yàn)?.6W時(shí)發(fā)生破裂的概率最大，因此設(shè)定漾濞MS6.4地震起始破裂點(diǎn)在破裂面上沿走向和沿傾向的位置分別為16.5 km和4.62 km(見(jiàn)圖3)，破裂速度為0.8β(β為震源處介質(zhì)剪切波速)，破裂面埋深3.39 km，地震未破裂至地表面.

破裂面上滑動(dòng)分布和破裂過(guò)程對(duì)近斷層地震動(dòng)具有顯著影響，目前尚沒(méi)有此次地震的震源破裂模型反演結(jié)果，本文采用了滑動(dòng)位錯(cuò)、滑動(dòng)速度、破裂速度以及破裂時(shí)間函數(shù)等在破裂面上隨機(jī)分布的震源運(yùn)動(dòng)學(xué)破裂模型(Graves and Pitarka 2010, 2015, 2016).王宏偉等(2021b)、Wang和Wen(2019)已采用該方法建立了2020年伽師MS6.4地震和2014年魯?shù)镸S6.6地震的震源運(yùn)動(dòng)學(xué)破裂隨機(jī)模型，并應(yīng)用于地震動(dòng)場(chǎng)的隨機(jī)有限斷層模擬.本文生成了100個(gè)漾濞MS6.4地震的震源破裂隨機(jī)模型.

根據(jù)Crust1.0全球地殼模型給出了研究區(qū)域地殼介質(zhì)水平成層模型，其中震源處介質(zhì)P波和S波波速分別為6.10 km·s-1和3.55 km·s-1，介質(zhì)密度2.74 g·cm-3(https:∥igppweb.ucsd.edu/～gabi/crust1.html).本文采用Xu等(2010)基于我國(guó)川滇地區(qū)中小震強(qiáng)震記錄建立的四段線性幾何擴(kuò)散模型(見(jiàn)表1)及S波品質(zhì)因子(QS=180f0.5)表示區(qū)域地殼介質(zhì)對(duì)地震波的衰減，其中P波品質(zhì)因子(QP)根據(jù)王勤彩等(2005)給出的云南地區(qū)QS/QP=1.5確定，即QP=120f0.5.采用王宏偉等(2021a)基于四川地區(qū)地震動(dòng)記錄建立的地震動(dòng)路徑持時(shí)經(jīng)驗(yàn)?zāi)Ｐ捅硎狙∕S6.4地震的地震動(dòng)路徑持時(shí).

場(chǎng)地放大效應(yīng)可表示為地殼介質(zhì)放大和淺層場(chǎng)地放大兩部分.地殼介質(zhì)放大采用美國(guó)國(guó)家地震減災(zāi)計(jì)劃(NEHRP)中B-C類場(chǎng)地分界對(duì)應(yīng)場(chǎng)地條件的放大效應(yīng)經(jīng)驗(yàn)?zāi)Ｐ?Atkinson and Boore, 2006)；淺層場(chǎng)地放大效應(yīng)采用Seyhan和Stewart(2014)建立的VS30(地表以下30 m深度范圍內(nèi)土層的等效剪切波速)、PGAr(VS30=760 m·s-1的參考基巖面的水平向峰值加速度中位值)相關(guān)的線性和非線性淺層場(chǎng)地效應(yīng)經(jīng)驗(yàn)?zāi)Ｐ捅硎?，本文中虛擬觀測(cè)點(diǎn)的VS30根據(jù)基于地形數(shù)據(jù)的全球VS30估計(jì)值確定(Heath et al., 2020，https:∥earthquake.usgs.gov/data/vs30/).與場(chǎng)地相關(guān)的地震動(dòng)高頻衰減參數(shù)κ=0.025 s(Xu et al., 2010).

地震應(yīng)力降表示地震前后破裂面上應(yīng)力變化，該參數(shù)對(duì)高頻(高于拐角頻率)地震動(dòng)具有控制性影響，但依靠目前技術(shù)條件尚無(wú)法直接測(cè)量，因此本文采用Yenier和Atkinson(2015)提出的方法對(duì)漾濞MS6.4地震的應(yīng)力降進(jìn)行估計(jì).我國(guó)大陸地區(qū)中小地震的地震應(yīng)力降主要分布在0.1～20.0 MPa(趙翠萍等，2011)，據(jù)此我們首先設(shè)定地震應(yīng)力降為1～20 MPa范圍內(nèi)固定間隔1 MPa的20個(gè)可能值，針對(duì)每次震源破裂隨機(jī)過(guò)程，利用隨機(jī)有限斷層方法模擬28個(gè)觸發(fā)觀測(cè)臺(tái)站的三分量地震動(dòng)加速度時(shí)程，并計(jì)算模擬記錄0.05～10 s的5%阻尼比的擬加速度反應(yīng)譜(PSA)，以100次隨機(jī)破裂過(guò)程的模擬記錄PSA的對(duì)數(shù)平均值表示模擬記錄的PSA，計(jì)算各臺(tái)站模擬記錄PSA殘差，即log10(PSAobs)-log10(PSAsim)，其中PSAobs和PSAsim分別表示觀測(cè)與模擬記錄的PSA，最后給出28個(gè)觸發(fā)臺(tái)站在0.05～10.0 s周期段模擬記錄PSA殘差均值及標(biāo)準(zhǔn)差.根據(jù)所有臺(tái)站模擬記錄PSA殘差的均值在0.05～10 s周期段的均值及標(biāo)準(zhǔn)差隨應(yīng)力降的變化(見(jiàn)圖2a、b、c)估計(jì)地震應(yīng)力降(Δσ)，本文以NS和UD分量模擬記錄PSA殘差標(biāo)準(zhǔn)差達(dá)到最小值時(shí)所對(duì)應(yīng)的應(yīng)力降和殘差均值分別作為地震應(yīng)力降和譜形一致調(diào)整系數(shù)，即Δσ=16 MPa，EW、NS、UD分量的譜形一致調(diào)整系數(shù)分別為100.00003、100.0075、100.27277.

采用估計(jì)的地震應(yīng)力降和譜形一致調(diào)整系數(shù)，針對(duì)每次震源破裂隨機(jī)過(guò)程，采用隨機(jī)有限斷層方法給出了近斷層053YBX臺(tái)站的三分量地震動(dòng)加速度時(shí)程模擬結(jié)果，并給出模擬記錄PSA的殘差絕對(duì)值及其在0.05～10 s周期段的平均值(見(jiàn)圖2d)，其中第34次震源破裂隨機(jī)過(guò)程的殘差絕對(duì)值均值最小，因此近似以此震源破裂隨機(jī)過(guò)程作為漾濞MS6.4地震震源破裂過(guò)程(見(jiàn)圖3).

圖2 所有臺(tái)站模擬記錄PSA殘差的均值在0.05～10 s周期段的均值及標(biāo)準(zhǔn)差隨應(yīng)力降的變化(a)、(b)、(c)分別為EW、NS和UD分量；(d)模擬記錄PSA的殘差絕對(duì)值在0.05～10 s周期段的均值隨隨機(jī)破裂過(guò)程的變化.Fig.2 Effect of stress drop on average and standard deviation over a period range from 0.05 to 10 s for PSA residual means for the simulated recordings at all stations(a), (b), and (c) for EW, NS, and UD components respectively. (d) Effect of stochastic rupture model on the mean of absolute residuals of PSAs of simulated recordings over a period range from 0.05 to 10 s.

圖3 漾濞MS6.4地震的震源破裂過(guò)程其中紅色五角星表示破裂起始點(diǎn)位置，黑色實(shí)線表示破裂傳播時(shí)間.Fig.3 Source rupture model of the Yangbi MS6.4 earthquakeBlack solid lines represent the rupture propagation time from the initial rupture point (highlighted by red star).

2.2 模擬結(jié)果

本文利用隨機(jī)有限斷層方法模擬了漾濞MS6.4地震中28個(gè)觸發(fā)強(qiáng)震動(dòng)臺(tái)站及震中附近區(qū)域內(nèi)(北緯23.67°—27.67°、東經(jīng)97.87°—101.87°)2823個(gè)虛擬觀測(cè)點(diǎn)的三維地震動(dòng)加速度時(shí)程，其中距離震中50、100和500 km范圍內(nèi)網(wǎng)格精度分別為0.025°、0.05°和0.1°，模擬輸入?yún)?shù)詳見(jiàn)表1.

表1 漾濞MS6.4地震隨機(jī)有限斷層模擬輸入?yún)?shù)Table 1 Input parameters for the stochastic finite-fault simulation of the Yangbi MS6.4 earthquake

所有虛擬觀測(cè)點(diǎn)模擬記錄三分量的PGA和PGV如圖4所示，模擬與觀測(cè)記錄的地震動(dòng)峰值(PGA、PGV)整體上均具有較好的一致性，模擬記錄可以很好的反映地震動(dòng)峰值隨距離的衰減規(guī)律及其近場(chǎng)飽和效應(yīng).同時(shí)對(duì)比了俞言祥等(2013)為我國(guó)第五代區(qū)劃圖編制建立的青藏地震區(qū)基巖場(chǎng)地(VS30≥500 m·s-1)水平向地震動(dòng)峰值預(yù)測(cè)方程給出的2823個(gè)虛擬觀測(cè)點(diǎn)的預(yù)測(cè)中位值.設(shè)定預(yù)測(cè)方程中長(zhǎng)軸方向與破裂面走向一致，虛擬觀測(cè)點(diǎn)VS30主要在500～900 m·s-1范圍內(nèi)，只有不超過(guò)5%的虛擬觀測(cè)點(diǎn)VS30小于500 m·s-1，因此虛擬觀測(cè)點(diǎn)可近似看作位于基巖場(chǎng)地.對(duì)比發(fā)現(xiàn)，模擬記錄峰值與預(yù)測(cè)中位值的大小、距離衰減特性均具有較好的一致性，盡管預(yù)測(cè)中位值在近場(chǎng)區(qū)域(震源距<10 km)整體上高于模擬值.

圖4 模擬與觀測(cè)記錄的PGA、PGV對(duì)比Fig.4 Comparisons between PGAs, and PGVs from the simulated and observed ground motions

為反映模擬記錄對(duì)不同周期地震動(dòng)的預(yù)測(cè)水平，圖5給出了漾濞MS6.4地震中28個(gè)觸發(fā)強(qiáng)震動(dòng)臺(tái)站的模擬記錄PSA的殘差，所有臺(tái)站模擬記錄PSA殘差主要在±0.5范圍內(nèi)，所有臺(tái)站模擬記錄PSA殘差的均值在零附近波動(dòng)，主要在±0.2的較小范圍內(nèi)，且并未發(fā)現(xiàn)明顯的周期相關(guān)性，這表明模擬記錄在整個(gè)周期段(0.05～10 s)均能較好的預(yù)測(cè)地震動(dòng)，模擬與觀測(cè)記錄PSA的幅值接近且譜形相似.

圖5 模擬和觀測(cè)記錄的PSA在0.05～10 s周期段的殘差Fig.5 Residuals of the PSAs at periods of 0.05～10 s between the simulated and observed ground motions

距離震中最近的三個(gè)觀測(cè)臺(tái)站(053YBX、053DLY、053YPX)的模擬記錄三分量加速度時(shí)程、速度時(shí)程和PSA如圖6所示.三個(gè)臺(tái)站模擬與觀測(cè)記錄的加速度和速度時(shí)程的波形、幅值以及PSA的幅值、譜形均具有較好的一致性.圖6同時(shí)給出了最接近大理州市區(qū)的虛擬觀測(cè)點(diǎn)(#2106，震源距約為35.5 km)的模擬記錄時(shí)程及PSA，其EW、NS和UD分量的PGA分別為122.9、146.7和73.6 cm·s-2，PGV分別為7.8、10.5和3.8 cm·s-1.

圖6 模擬記錄加速度時(shí)程、速度時(shí)程及PSAFig.6 Acceleration and velocity time histories and PSAs for the simulated ground motions

圖7對(duì)比了漾濞MS6.4地震破裂前方(SE)和破裂后方(NW)震源距相同的虛擬觀測(cè)點(diǎn)的模擬記錄，#2504和#2740震源距均為20 km，#1726和#2101震源距均為50 km，#2740和#2101位于破裂前方，而#2504和#1726則位于破裂后方，破裂前方臺(tái)站記錄的加速度和速度時(shí)程的幅值更大、持時(shí)更短，同時(shí)其PSA幅值更高，表現(xiàn)出明顯的破裂方向性效應(yīng)，說(shuō)明模擬記錄體現(xiàn)了震源破裂方向性對(duì)地震動(dòng)的影響.

圖7 震源距相同的地震破裂前后方虛擬觀測(cè)點(diǎn)模擬記錄加速度時(shí)程、速度時(shí)程和PSAFig.7 Acceleration and velocity time histories, as well as PSAs simulated for both stations in each pair, which share the same hypocentral distances and are located in and opposite the direction of source rupture, respectively

根據(jù)《中國(guó)地震烈度表 GB/T 17742-2020》，利用虛擬觀測(cè)點(diǎn)的三分量模擬記錄得到了合成加速度和速度時(shí)程，給出了三分量合成PGA和PGV，最終計(jì)算了儀器測(cè)定地震烈度.基于虛擬觀測(cè)點(diǎn)模擬記錄的儀器測(cè)定地震烈度空間分布如圖8所示，基于模擬記錄的地震烈度空間分布的大致走向、范圍與云南省地震局正式發(fā)布結(jié)果(見(jiàn)圖8a)基本一致，斷層附近區(qū)域極震區(qū)的最大烈度可達(dá)到Ⅷ度，與宏觀烈度圖的最大烈度Ⅷ度相一致，其中漾濞縣城及鄰近鄉(xiāng)鎮(zhèn)位于Ⅷ度區(qū)，大理市區(qū)位于Ⅵ度區(qū)，均與正式發(fā)布結(jié)果相同.極震區(qū)近似為橢球狀，受震源破裂方向性影響，相同距離上震中SE方向的烈度普遍高于NW方向.此外，沿洱海西側(cè)的烈度普遍偏高(Ⅶ度)，出現(xiàn)高烈度異常區(qū)，該地區(qū)的VS30普遍較小，局部場(chǎng)地放大效應(yīng)可能是高烈度異常的主要原因.

圖8 基于模擬記錄的漾濞MS6.4地震的地震烈度與(a)云南省地震局正式發(fā)布的地震烈度圖(藍(lán)實(shí)線)、(b)基于地震烈度臺(tái)站觀測(cè)記錄的地震烈度(藍(lán)實(shí)線)對(duì)比Fig.8 Seismic intensity of the Yangbi MS6.4 earthquake derived from the simulated ground motions comparison with (a) seismic intensity officially released by the Earthquake Administration of Yunnan Province (blue lines) and (b) seismic intensity based on recordings obtained by seismic intensity stations (blue lines)

自2018年以來(lái)，國(guó)家地震烈度速報(bào)與預(yù)警工程建設(shè)初見(jiàn)成效，此次漾濞MS6.4地震中震中距100 km以內(nèi)共獲得地震烈度臺(tái)記錄86組，基于烈度臺(tái)觀測(cè)記錄的儀器測(cè)定地震烈度的等震線如圖8b所示，其等震線長(zhǎng)軸近東西向，Ⅵ、Ⅶ和Ⅷ度區(qū)范圍均偏小，與基于模擬記錄的儀器測(cè)定地震烈度和云南省地震局發(fā)布的地震烈度的結(jié)果有明顯差異，這可能與地震烈度臺(tái)站密度較小有關(guān).

3 結(jié)論

本文利用漾濞MS6.4地震的震源破裂隨機(jī)過(guò)程在28個(gè)觸發(fā)強(qiáng)震動(dòng)觀測(cè)臺(tái)站的隨機(jī)有限斷層三維地震動(dòng)模擬記錄，結(jié)合強(qiáng)震動(dòng)觀測(cè)記錄，估計(jì)了此次地震的地震應(yīng)力降及震源破裂過(guò)程，給出了漾濞MS6.4地震在震中附近范圍2823個(gè)虛擬觀測(cè)點(diǎn)的模擬記錄三分量加速度時(shí)程，結(jié)果表明：

(1) 模擬記錄PGA、PGV與觀測(cè)記錄、地震動(dòng)預(yù)測(cè)方程預(yù)測(cè)中位值均具有較好的一致性，反映了地震動(dòng)峰值隨距離的衰減規(guī)律、近場(chǎng)飽和效應(yīng)及破裂方向性效應(yīng)；

(2) 模擬記錄與觀測(cè)記錄PSA的幅值接近、譜形相似，模擬記錄在0.05～10 s周期段均很好地預(yù)測(cè)了地震動(dòng)；

(3) 基于三分量模擬記錄給出了此次地震的儀器測(cè)定地震烈度，與云南省地震局發(fā)布烈度圖的結(jié)果基本一致，極震區(qū)最大烈度可達(dá)到Ⅷ度，震源破裂方向性導(dǎo)致震中SE方向的烈度普遍高于NW方向，受局部場(chǎng)地條件影響沿洱海西側(cè)出現(xiàn)高烈度異常區(qū).

致謝感謝中國(guó)地震局工程力學(xué)研究所為本研究提供的強(qiáng)震動(dòng)觀測(cè)數(shù)據(jù)及地震烈度臺(tái)站的儀器烈度數(shù)據(jù).感謝中國(guó)地震局地球物理研究所房立華研究員提供的漾濞地震序列震源精定位結(jié)果.感謝兩位匿名評(píng)審專家對(duì)本文修改和完善提出的寶貴意見(jiàn).