2022年1月8日青海门源MS6.9地震序列重定位和震源机制解研究＊

许英才 郭祥云 冯丽丽

1）中国银川 750001 宁夏回族自治区地震局

2）中国北京 100081 中国地震局地球物理研究所

3）中国西宁 810001 青海省地震局

引言

据中国地震台网测定，北京时间2022年1月8日01时45分27秒青海省海北藏族自治州门源县（37.77°N，101.26°E）发生MS6.9地震，震中距离青海门源、青海西宁和甘肃兰州分别为54 km，141 km和278 km.震中附近的甘肃省河西五市震感强烈（中国新闻网，2022），而且兰州、白银、定西、平凉等市震感明显，震中附近区域的冷龙岭断裂西段产生了一条长约22 km的地表破裂带（青海省地震局，2022），在其地表及冰面可见挤压鼓包及连续张裂隙等其它破裂组合形式，鼓包高约1.5 m；震中极震区烈度为Ⅸ度（中国地震局，2022），等震线的椭圆长轴呈WNW走向且长约200 km，其短轴长约153 km.截至2022年1月12日24时，震中区共发生ML1.0以上余震606次（图1），其中ML5.0—5.9地震3次，ML4.0—4.9地震7次，ML3.0—3.9地震40次，ML2.0—2.9地震187次，ML1.0—1.9地震369次，最大余震为1月12日18时20分MS5.2地震.本次门源MS6.9地震发生在青藏高原东北缘，由于青藏地块NE向的强烈挤压作用（Tapponnieret al，2001；李祥等，2016；石富强等，2018），该区域受到阿拉善地块阻挡，从而导致地壳物质沿SE方向的挤出，形成了以左旋走滑及逆冲断裂系为主的构造体系，其最大主压应力轴呈现NE到ENE方向的变化.该地震发生在冷龙岭断裂（F1）与托莱山断裂（F2）之间，也位于对地壳速度场具有强烈控制作用的广义海原断裂带西段，震中50 km内发震时间最近的为2016年青海门源MS6.4地震（震源机制为逆冲型），震中100 km内震级最大的为1927年甘肃古浪M8.0地震.此次门源MS6.9地震发生在M6.0地震危险区（徐锡伟等，2017），该地震所处的河西走廊是“一带一路”的重要纽带，也是很多重大国家工程的必经之路，因此，研究该地震的发震构造以及如何进一步评估该区未来的强震活动趋势，成为震后亟需研判的重要科学问题.

图1 门源MS6.9地震震中和MS≥6.0历史强震分布（a）、地震序列分布（b）及其M-t图和日频次图（c）F1：冷龙岭断裂；F2：托莱山断裂；F3：昌马—俄博断裂；F4：祁连山北缘断裂；F5：皇城—双塔断裂Fig. 1 Epicenter of MS6.9 Menyuan earthquake and distribution of historical MS≥6.0 strong earthquakes near the main shock epicenter （a），earthquake sequence （b） and their M-t plot and daily frequencies （c）F1：Lenlongling fault；F2：Tuolaishan fault；F3：Changma-Ebo fault；F4：Northern Qilianshan fault；F5：Huangcheng-Shuangta fault

精确确定主震与余震序列的空间展布及震源参数对于震后趋势的判定具有重要意义.为此，本文拟基于青海、甘肃、宁夏及内蒙古区域地震台网的震相观测资料和波形数据，利用多阶段定位（multi-step locating）方法（Longet al，2015）对门源MS6.9地震序列进行重定位，通过gCAP （generalized cut and paste）方法（Zhu，Ben-zion，2013）反演主震和MS3.4以上余震序列的震源机制解和矩心深度，然后采用应力场与震源机制模拟的方法（万永革，2020）计算已有区域应力场体系在门源MS6.9地震震源机制两个节面上产生的相对剪应力和相对正应力，最后结合区域地质构造资料初步分析此次门源MS6.9地震的发震机制，以期为深入理解门源地震序列的活动特征、区域断裂的活动方式以及进一步评估该区未来强震活动趋势等提供参考.

1 门源MS6.9地震序列重定位

本文利用多阶段定位方法（Longet al，2015；龙锋等，2021）对门源地震序列早期余震进行重定位.与仅使用双差定位不同，多阶段定位方法的优势在于可确保初始值和模型更准确，从而使定位结果更可靠.为确保获取足够多的样本使定位结果更为稳健，同时考虑到震中附近区域台网台站空间分布的不均匀性，本文收集了2011年1月1日至2022年1月12日门源地区ML≥2.0地震事件的震相到时资料进行定位，定位过程如下：① 初始P波速度模型（表1）参考周民都等（2006）的体波走时反演结果和王椿镛等（1995）的地震测深剖面结果，地壳厚度和平均波速比分别取56 km和1.73 km/s （杨志高，张雪梅，2016；左可桢，陈继锋，2018），通过利用Hypo2000 （Klein，2002）对地震序列中至少具有6个以上台站记录且震中距不超过240 km的ML≥2.0地震事件进行初定位，满足该条件的地震有339次（其中属于门源MS6.9地震序列的地震有173次），Pg和Sg震相分别有4425条和3383条，基于和达曲线剔除2.5倍均方差之外的离散数据（图2a）；② 筛选出具有6个以上台站记录且最大台站方位角间隙小于150°的震相信息，然后根据初始速度模型（表1）对不同深度的速度值进行适当细分（共26层），即5 km深度以内用0.5 km步长细分，5—10 km深度用1—2 km步长细分，大于10 km深度的主要采用2—3 km步长进行细分，之后再采用Velest方法（Kissling，1988；Kisslinget al，1994）进行最小一维速度模型和台站校正反演；③ 根据反演后的速度模型及台站校正，再使用Hypo2000（Klein，2002）进行二次定位；④ 采用10 km的搜索半径，基于上述所得最小一维速度模型（表2，图2b）对震相走时资料进行HypoDD定位（Waldhauser，Ellsworth，2000），采取三轮最小平方QR分解（sparse linear equations and least squares，缩写为LSQR）法（Paige，Saunders，1982）迭代得到双差定位的结果，最终获得300次地震的精确位置，其中门源地震序列事件157次.东西向、南北向和垂直向的平均定位误差分别为0.48 km，0.47 km和0.76 km，平均走时残差为0.22 s.

表1 门源地区的P波初始速度模型Table 1 Initial P-wave velocity model in Menyuan area

图2 门源震源区的震相和达曲线（a）及最小一维速度模型（b）（红虚线代表拟合直线的2.5倍均方差）Fig. 2 Wadadi diagram （a） and minimum 1-D velocity model （b） of Menyuan source region（Red dashed lines represent the limits for 2.5 RMS of the fitting line）

表2 门源地区的P波最小一维速度模型Table 2 Minimum1-D P-wave velocity model in Menyuan area

重定位后的门源MS6.9地震序列在空间上相对集中（图3，4），并大体呈现WNW-ESE向的展布，其中地震序列走向在余震区的中段略向东南方向倾斜，AA′剖面显示整个序列相对集中于长约30 km的区域，而震后24个小时内的序列长度约为25 km，与地表破裂带的长度22 km大体一致.地震序列早期地震活动主要分布在余震区的西段和中段，主震发生后24小时内余震数量显著高于后续4天内的数量，其中门源1月8日01时45分的MS6.9主震及紧跟其后的02时09分的MS5.1余震均位于西段.该地震序列整体上略呈双侧破裂的特征，其中西向破裂较东向破裂明显，余震区整体地震活动主要介于5—10 km深度之间，其中优势震源深度分布在7—8 km之间，约占整个序列的68%.而地震序列后期地震活动主要集中在余震区的东段，最大余震—1月12日18时20分MS5.2地震即位于余震区东段.BB′，CC′及DD′剖面均显示：主震的发震断层倾角较陡，近似于直立且余震区西段、中段和东段的剖面深度较为平行；结合AA′剖面来看，余震活动的空间范围自西向东逐步有所收敛，并且在约10—13 km深度之间存在小震稀疏段，小于10 km深度的余震活动不仅频次相对最多而且主要以序列的早期地震活动为主，约占整个序列的86%，而大于 13 km 深度的地震活动频次相对较少，主要以序列的中、后期活动为主，仅占整个序列的14%.Yang等（2022）基于ML1.0以上的地震观测报告，通过双差重定位得到的门源MS6.9主震及其余震震源深度分布也主要集中于5—10 km深度之间，而在11—20 km深度之间的余震十分稀少，本文ML2.0以上地震的重定位结果与其大体一致.

图3 2022年1月8日至12日门源MS6.9地震序列重定位震中分布（a） 重定位后地震随震源深度分布图；（b） 重定位后地震随距离主震的离逝时间t分布图Fig. 3 Relocated epicenters of the MS6.9 Menyuan earthquake sequence in the period of 8 to 12 January 2022（a） Epicenters distribution with focal depth after relocation；（b） Epicenters distribution of with elapsed time t from the origin time of the main shock

2 主震及MS≥3.4余震的震源机制计算

本文使用gCAP方法（Zhu，Helmberger，1996；Zhu，Ben-zion，2013）反演门源MS6.9主震及余震序列中MS≥3.4地震的震源机制和震源矩心深度，基于本文多阶段定位产出的最小一维速度模型（模型3，表2和波速比1.73），利用f-k方法（即频率-波数法）计算格林函数（Zhu，Rivera，2002），其中对于门源MS6.9主震本文还计算了模型1 （表1和波速比1.75）和模型2 （表1和波速比1.73）的格林函数，以分析这三种不同的速度模型对门源MS6.9主震震源机制及其深度结果的影响.反演过程中，数据预处理流程如下：对于主震，选取震中距处于130—360 km范围内的台站，对波形进行去均值、去趋势、去倾斜和去除仪器响应，重采样频率设为0.2 s，并将三分量的水平东西向（E）、水平南北向（N）和垂直向（Z）旋转至径向（R）、切向（T）和垂直向（Z），然后人工标注P波到时.使用gCAP反演时，设定体波窗长为35 s，面波窗长为70 s，并对体波和面波进行四阶巴特沃斯（Butterworth）带通滤波，其滤波频段分别为0.04—0.2 Hz和0.05—0.1 Hz，断层面参数走向、倾角和滑动角的搜索步长均为5°，震级搜索步长为0.1.

图4 沿地震序列长轴走向的震源深度剖面AA′以及垂直余震区长轴走向的震源深度剖面BB′，CC′和DD′ （剖面宽度为剖面线两侧各10km）Fig. 4 Source depth profiles AA′ along the major axis of the earthquake sequence and BB′，CC′ and DD′ along the strike perpendicular to the major axis of aftershock region （Projection width for each side is 10 km for section）

通过计算得到门源MS6.9主震和MS≥3.4余震的震源机制和震源矩心深度结果，为对比速度模型对震源机制解和深度的影响，本文基于模型1—3计算了主震gCAP反演结果，如图5所示.由该图可知，模型1结果中震源矩心深度为3 km，其波形拟合图整体拟合率大于0.7的占73%，模型2拟合率大于0.7的也占73%，模型3拟合率大于0.7的可达78%.这表明，相对于模型1和模型2来说，模型3不仅拟合率相对较高，各个理论波形相对实际波形的平均移动时间也较短，而且模型3的误差曲线在形态上相对更光滑，整体来看模型3与实际观测波形的拟合相对较好，更接近震源区的真实地壳速度结构，因此不同的速度模型对震源机制和深度均有一定的影响.本文采用模型3的反演结果4 km作为主震的最终震源矩心深度.

图5 基于模型3 （c）的门源MS6.9地震震源机制-深度误差图（左）和最佳深度处波形拟合图（右）拟合波形下方的两行数字分别为理论波形（红色）相对实际波形（黑色）的移动时间(单位：s)以及二者的相关系数，波形左侧第一行给出了台站名和方位角（单位：度），第二行给出了震中距（单位：km）和相对偏移时间（单位：s），台站波形按震中距排列Fig. 5 Source mechanism-depth error diagram （left） and waveform fitting diagram at optimum depth （right） of the MS6.9 Menyuan earthquake with velocity model 3 （c）The numbers of two rows beneath the traces are the time shifts （in second） of synthetics （red） relative to the observations（black） and the corresponding cross-correlation coefficients，respectively. The upper-left corner are stations and azimuths （in degree），respectively. The lower-left corner numbers represent the epicentral distance （in km）and the relative offset time （in second）. The waveforms of stations are sorted in epicentral distance

图5 基于速度模型1 （a）和2 （b）的门源MS6.9地震震源机制-深度误差图（左）和最佳深度处波形拟合图（右）拟合波形下方的两行数字分别为理论波形（红色）相对实际波形（黑色）的移动时间（单位：s）以及二者的相关系数，波形左侧第一行给出了台站名和方位角（单位：度），第二行给出了震中距（单位：km）和相对偏移时间（单位：s），台站波形按震中距排列Fig. 5 Source mechanism-depth error diagram （left） and waveform fitting diagram at optimum depth （right） of the MS6.9 Menyuan earthquake with velocity models 1 （a） and 2 （b）The numbers of two rows beneath the traces are the time shifts （in second） of synthetics （red） relative to the observations（black） and the corresponding cross-correlation coefficients，respectively. The upper-left corner are stations and azimuths （in degree），respectively. The lower-left corner numbers represent the epicentral distance （in km）and the relative offset time （in second）. The waveforms of stations are sorted in epicentral distance

由于不同模型的震源机制解结果存在一定差异，本文将这三种不同模型的震源机制结果与国内外机构或作者的结果进行对比，如图6所示，图中同时给出了采用震源机制中心解的方法（万永革，2019）所得的这些震源机制的最小空间旋转角以表示其离散程度.本文的震源机制最小空间旋转角介于16.69°—24.84°之间，国内机构测定震源机制的最小空间旋转角范围为2.77°—41.45°，国外机构的震源机制解旋转角为5.32°—18.08°，整体看来最小空间旋转角主要相对集中在16.69°—19.37°之间，其中包括本文模型2和模型3的结果，以及CENC、赵韬等、USGS和GFZ的结果，说明这些结果给出的震源机制十分接近.最小空间旋转角的差异从侧面反映了震源机制的离散程度，这主要是由于不同机构或作者使用的资料尺度、方法和速度模型等参数不同所致.根据前文速度模型和实际波形的拟合率情况，以及结合震源机制中心解的最小空间旋转角，本文采取模型3的震源机制作为门源MS6.9地震的最终震源机制结果.

图6 本文门源MS6.9地震的震源机制结果与其他作者或机构结果的对比CENC、青海地震台、王卫民等、赵翠萍等、赵韬等的结果引自Seismology小组（2022），郭祥云等和韩立波等的结果引自中国地震局地球物理研究所（2022），其它结果引自European-Mediterranean Seismological Centre （2022）Fig. 6 Comparison of focal mechanism results of the MS6.9 Menyuan earthquake in this paper with those from other authors or institutionsThe results of CENC，Qinghai seismic network，Wang Weimin et al，Zhao Cuiping et al，Zhao Tao et al are from Seismology Group （2022），those of Guo Xiangyun et al and Han Libo et al are from Institute of Geophysics，China Earthquake Administration （2022），and others are from European-Mediterranean Seismological Centre （2022）

余震序列中MS≥3.4地震的震源机制主要为走滑型地震（表3和图7），震源矩心深度介于3—7 km，较接近于主震的深度.这些MS≥3.4地震的震源机制存在较明显的分区特征，余震区西段地震的震源机制P轴走向为NE向，东段的P轴走向主要为近EW向.

图7 2022年1月8—12日门源MS6.9主震和MS≥3.4余震的震中分布及其震源机制（震中位置据表3）Fig. 7 Epicentral distribution and focal mechanisms of the MS6.9 main shock and MS≥3.4 aftershocks in Menyuan during 8-12 January 2022（Epicenters according to Table 3）

表3 2022年1月8日至12日门源地震序列MS≥3.4地震的震源机制解Table 3 Focal mechanism solutions of the MS≥3.4 events in the Menyuan earthquake sequence from January 8 to January 12，2022

3 门源MS6.9地震的震源机制与应力场关系

本文基于李祥等（2016）的门源地区的构造应力场结果（主压应力轴的走向为37.7°，倾伏角为8.7°，主张应力轴的走向为174.3°，倾伏角为78.1°，以及应力形因子为0.69），使用震源机制和应力场模拟的方法（Wan，2010；万永革，2020），分别计算了2022年门源MS6.9和2016年门源MS6.4两次地震震源机制的两个节面所产生的相对剪应力和正应力，其中2016年门源MS6.4地震震源机制解采用尹欣欣等（2018）给出的断层面解.

门源地区应力体系下的震源机制模拟及相对剪应力和相对正应力的结果如图8所示，可见：应力张量在2022年门源MS6.9地震震源机制节面Ⅰ（走向290°，倾角81°）上产生的相对剪应力和相对正应力分别为0.638和-0.606，剪滑角（即剪应力的滑动角）为70.7°，在节面Ⅱ（走向197°，倾角74°）上的相对剪应力和相对正应力分别为0.522和-0.175，剪滑角为125.5°；而应力张量在2016年MS6.4地震震源机制节面Ⅰ（走向143°，倾角40°）上的相对剪应力和相对正应力分别为0.995和0.057，剪滑角为106.5°，节面Ⅱ（走向347°，倾角53°）上的相对剪应力和相对正应力分别为0.833和0.105，剪滑角为91.8°.2022年门源MS6.9地震在WNW向冷龙岭断裂上产生的相对剪应力为0.638，未达到最大剪应力（最大相对剪应力为1），说明此次地震的破裂面走向（即冷龙岭断裂走向）与最大释放应力的破裂面走向有一定的偏差，构造区积累的应力未得到完全的释放.这与2016年门源MS6.4地震不同，该地震的两个震源机制节面的剪应力几乎都达到最大值，应力可能得到了充分的释放.

图8 门源地区应力体系下的震源机制模拟及相对剪应力（a）和相对正应力（b）Fig. 8 Simulated focal mechanisms and relative shear stress （a） and relative normal stress （b） under stress system in Menyuan region

4 讨论与结论

根据已有地质构造（Lasserreet al，2002；袁道阳等，2004；郭鹏等，2017），门源MS6.9地震序列所在的冷龙岭断裂是一条走向为NW60°到NW70°且长约127 km的全新世活动断裂，该断裂属于北祁连山活动断裂带的分支，主要呈以左旋走滑为主且兼具逆冲特征.重定位结果显示本次MS6.9地震序列发生在冷龙岭断裂西段，呈WNW向的优势分布，也与冷龙岭断裂的走向基本一致，而且MS6.9主震震源机制解的节面 Ⅰ 走向与冷龙岭断裂的走向基本一致，已有应力场体系在该地震震源机制WNW向的节面 Ⅰ 所产生的相对剪应力略大于NNE向的节面 Ⅱ ，为此WNW向的节面 Ⅰ 可能是门源MS6.9地震的发震构造，即该地震是在青藏高原NE向的挤压作用下，发生在冷龙岭断裂西段上的一次左旋走滑型地震.MS≥3.4余震序列的震源机制存在较明显的分区特征，即余震区东段的震源机制解P轴方向大体呈EW向，西段的震源机制解P轴方向主要呈NE向，由于P轴走向取决于破裂面的走向，因此余震震源机制的分区特征从侧面反映了东段与西段余震破裂面形状可能有所不同，东段与西段的发震构造也可能存在一定的差异，中、西段区域可能是构造变形模式的转换区域，且转换区域东、西侧的构造变形也存在一定的差异.

本文所获门源MS6.9地震的震源机制解与多个机构或作者的震源机制中心解基本一致，最小空间旋转角与其他作者或机构的大部分结果也大体一致，震源矩心深度为4 km，这与韩立波等测定的矩心深度结果4 km （中国地震局地球物理研究所，2022）相一致.本文重定位的主震震源深度为7.8 km，考虑到重定位和gCAP方法的深度结果分别代表初始破裂深度和震源矩心深度，而且该地震破裂尺度较大，破裂时间持续10 s左右（中国地震局地球物理研究所，2022），一般初始破裂深度要大于震源矩心深度，矩心深度的物理意义为地震破裂过程中的平均能量深度.王月等（2021）给出了云南漾濞MS6.4地震序列中较大地震的震源机制，其结果显示初始破裂深度与矩心深度平均差距在2—5 km之间，其中漾濞前震MS5.2地震的重定位与CAP结果深度相差9 km左右，说明初始破裂深度普遍大于矩心深度，而且重定位和波形拟合反演的不同深度可能反映了震源破裂过程比较复杂，也可能是资料限制和方法不同所导致.另外本文关于门源MS6.9地震早期序列的重定位结果显示，在震源深度处存在长约30—40 km且宽约5 km的地震稀疏段，认为该现象可能与门源地区下方存在高速异常体有关.2022年门源MS6.9地震序列所在区域处于2016年门源MS6.4地震西侧，已有研究（左可桢，陈继锋，2018）表明，2016年门源MS6.4地震西侧10—15 km地壳深度处正好存在显著的P波高速异常区且其P波速度结构异常分布特征在10 km与15 km处十分相似，而高速异常区表明该区域内的介质强度相对较高，中小地震相对偏少，而且介质不容易破裂（吴建平等，2009），这有利于更高强度的应力积累及集中释放，为强震的孕育和破裂提供条件，因此该区域存在的高速异常体所代表的障碍体也可能是2022年门源MS6.9地震发震机制的一种解释.

门源地震震源深度处的地震稀疏段长约30—40 km，十分接近于余震序列的优势展布长度30 km，也大于地表破裂长度22 km，为此本文将长度30 km≤L≤40 km代入不同研究者统计的震源破裂长度L与震级的关系式：MS＝3.3＋2.1lgL（郭增建，秦保燕，1979）、MW＝4.33＋1.49lgL（Wells，Coppersmith，1994），MS＝1.860lgL＋3.821 （龙锋等，2006），相应的估算震级分别为MS6.4—6.7、MW6.5—6.7和MS6.6—6.8，震级差仅为0.2左右，而且十分接近2022年门源地震的实际地震震级MS6.9；再根据震源破裂面积A与震级MS的关系式MS＝0.954lgA＋4.134 （龙锋等，2006），将门源MS6.9地震早期地震序列稀疏段面积A＝30 km×5 km以及A＝40 km×5 km代入，所得震级分别为MS6.2和MS6.3，大小基本一致，这与2016年门源地震的实际地震震级MS6.4基本接近.加之现今应力场体系在门源MS6.9地震震源机制两个节面所产生的相对剪应力均未达到最大，即应力未得到完全的释放，该地区仍存在一定的应力积累，本文认为门源地区仍然存在发生强震的危险性，孕震强度水平可能为MS6.2—6.8.

本文根据2022年门源MS6.9地震序列的重定位结果、MS≥3.4地震的震源机制解以及主震震源机制两个节面的相对剪应力等，结合地质构造资料以及震级与破裂关系初步分析和讨论了2022年门源MS6.9地震的发震构造和其强震危险性，主要获得了以下结论：

1） 重定位后的门源早期地震序列157次ML≥2.0地震震中主要沿冷龙岭断裂西段呈WNW-ESE向展布，优势长轴长约30 km，宽约5 km，初始破裂深度主要集中于7—8 km之间，占序列的68%，而且震源深度由西向东整体呈逐步收敛的特征.在2022年1月8日至12日所发生的三次MS≥5.0地震中，1月8日门源MS6.9主震和MS5.1地震位于余震区的西段，1月12日MS5.2地震位于余震区东段，大体呈双侧破裂的特征.

2） 主震震源机制解的节面 Ⅰ 为走向290°、倾角81°、滑动角16°，节面 Ⅱ 为走向197°、倾角74°、滑动角171°；10次MS≥3.4余震的震源机制类型以走滑型为主，震源矩心深度为3—7 km，P轴走向从余震区西段到东段之间大体呈NE向到EW向的变化，并在东、西段呈现明显的分区特征，表明余震区东、西段的构造变形模式存在着较为显著的差异.

3） 应力场和震源机制模拟结果显示，现今已有应力场体系在2022年门源MS6.9地震震源机制解节面 Ⅰ 上产生的相对剪应力为0.638，而在节面 Ⅱ 上产生的相对剪应力为0.522，表明2022年门源MS6.9地震不是发生在构造应力场最大剪应力的释放节面上，地震破裂面的形状不是最大剪应力的断层面形状，而2016年门源MS6.4地震发生在构造应力场剪应力最大的区域，震源机制节面十分接近剪应力的最大释放节面.

4） 结合地质构造、震源机制和余震展布，本文初步认为2022年1月8日门源MS6.9地震震源机制的WNW向节面为该地震的破裂面，其破裂方式为左旋走滑，且其发震构造可能与冷龙岭断裂西段有关.

5） 根据重定位结果和震级-破裂关系以及剪应力结果可知，门源MS6.9震源区下方可能存在高速异常体，而且该地区的应力未得到充分释放，存在一定的应力积累，该地区仍然具有一定的强震危险性.

本文的图件利用GMT软件绘制.防灾科技学院万永革教授提供了震源机制中心解等程序，四川省地震局龙锋高级工程师对多阶段定位方法予以指导，中国地震局地球物理研究所“国家数字测震台网数据备份中心”（doi：10.7914/SN/CB）为本研究提供了地震波形数据，两位审稿专家为本文提供了建设性的修改意见，对本文质量提升给予了很大的帮助.作者在此一并表示感谢！