Relocation of the 2014 Ludian MS6.5 earthquake and its large aftershock sequence
-
摘要: 本文采用云南测震台网的观测报告数据，利用双差定位方法对2014年鲁甸MS6.5地震及其强余震序列进行了重定位，获得了3 658个地震事件的震源参数。重定位后地震序列的震中分布显示，余震分布存在两个优势方向，分别为近EW向和SES向，呈共轭型分布，近EW向条带展布长度为30 km，SES 向条带展布长度为20 km；震源深度的分布显示，地震序列总体表现为主震附近震源较深，沿近EW向和SES向逐渐变浅，地震序列的震源深度主要分布在4—20 km范围内。截至2017年2月28日，鲁甸MS6.5地震震源区共发生（同一天发生的一组地震算一次）MS≥4.5强余震4次。重定位后的鲁甸4次强余震序列震中分布存在差异：2014年9月10日和10月27日两次强余震序列的展布特征与主震相同，而2016年和2017年另外两次强余震的后续余震仅分布在强余震的周边，与主震序列明显不同。综合重定位后余震序列分布、震源区地质调查资料以及前人研究认为，鲁甸地震的4次强余震序列是区域应力场和主震引发的震源区应力场共同作用的结果，2014年9月10日和10月27日的两次强余震序列主要受主震引发的震源区应力场的影响；而2016年和2017年两次强余震序列则主要受区域应力场的影响。Abstract: We use catalog data recorded by Yunnan seismic network to relocate the 2014 Ludian MS6.5 earthquake and its aftershock sequence with double-difference location algorithm, and obtain focal parameters of 3 658 earthquakes. The results show that the relocated epicenter distribution of Ludian MS6.5 earthquake sequence exhibits two predominant directions, approxi-mately EW and SES, and conjugated to each other. The length of the EW strip is about 30 km, while the length of SES strip is about 20 km. The Ludian earthquake sequence mainly distri-buted between 4−20 km in depth. It is deeper around the mainshock and becomes shallower along the strips. There are fourMS≥4.5 large aftershocks （the earthquakes that occurred on the same day counts as one） occurred in the Ludian MS6.5 earthquake source region till February 28, 2017. The distributions of four strong aftershock sequences are different from each other after relocation. The distributions of aftershock sequences of MS4.7 occurred in September 10 and MS4.6 occurred in October 27 are similar with the mainshock sequence, while the aftershock sequences located around 2016 MS4.7 earthquake and 2017 MS4.9 earthquake are different from the mainshock sequence. Synthesizing the distribution of relocated Ludian MS6.5 earthquake sequence, geologic data in studied area and previous researches, we draw the following conclusion: the distribution of four strong aftershock sequences are resulted from combined action of regional stress field and source stress field caused by the mainshock, the strong aftershock sequences occurred in September 10 and October 27 in 2014 are mainly influenced by source stress field induced by the mainshock, while the strong aftershock sequences occurred in 2016 and 2017 are mainly influenced by regional stress field.
-
引言
据中国地震台网测定，北京时间2014年8月3日16时30分，云南省昭通市鲁甸县发生MS6.5地震，震中位置为（103.34°E，27.10°N），震源深度为6.5 km。根据云南省测震台网测定，截至2017年2月28日在鲁甸震源区共记录到ML≥1.0地震6 265次，其中ML1.0—1.9地震4 831次，ML2.0—2.9地震1 279次，ML3.0—3.9地震140次，ML4.0—4.9地震14次，ML≥6.0地震1次，其中最大余震为2017年2月8日的MS4.9地震。
鲁甸地震震源区所处的昭通—莲峰断裂带位于川滇地块与华南地块之间的边界带上，是大凉山次级地块向SE向运动的前缘，起着调节地块SE向运动的作用（闻学泽等，2013；张彦琪等，2016）。该断裂带自2003年以来强震活动明显增加，共有5次中强地震（双震或震群计为一次地震）发生，分别为2003年鲁甸MS5.0和MS5.1双震、2004年鲁甸MS5.6地震、2006年盐津MS5.1双震、2012年彝良MS5.7和MS5.6双震和2014年鲁甸MS6.5地震。其中仅鲁甸地区就发生了3次MS≥5.0强震，并且在2014年鲁甸MS6.5地震发生两年之后连续发生了多次MS≥4.0的晚期强余震。闻学泽等（2013）基于区域活动构造和动力学、重定位后的小震分布和震源机制解、历史地震资料和GPS形变场等多学科信息综合研究认为，昭通断裂带存在强震或大震的中长期危险背景。徐锡伟等（2014）关于鲁甸地震震源区库仑应力的计算结果表明，鲁甸地震可能对周边活动断层系历史地震空段产生应力加载作用，因此该地区中强地震的危险性不容忽视。
2014年鲁甸MS6.5地震发生后，王未来等（2014）、张广伟等（2014）、徐甫坤等（2014）和房立华等（2014）均使用双差定位方法对鲁甸地震的余震序列进行精定位研究，重定位结果显示鲁甸地震余震序列的震中位置呈共轭型分布，分别为近EW向和SES向两支；结合震源区断层分布、主震震源机制解和烈度分布等资料分析认为，鲁甸地震的发震断层为NNW向的包谷垴—小河断裂。张勇等（2014）使用近震资料和远震体波资料反演得到鲁甸地震的破裂过程，结果显示鲁甸地震破裂开始0—2 s内，破裂主要发生在近EW向的断层上，随后在近NS向断层开始破裂，并释放了大部分地震矩，分析认为鲁甸地震是一次在NW向主压应力与NE向主张应力的统一应力场下两条共轭断层先后发生破裂的一次复杂地震事件。何骁慧等（2015）利用主震与参考地震之间的P波到时差及CAP （cut and paste）方法反演得到的时移信息，测量得出鲁甸地震质心位置与破裂起始位置之间的差异，确定地震断层面解中走向为160°（近NS向）的节面为真实破裂面。李艳娥等（2015）使用主地震相对定位法对鲁甸地震序列进行重新定位，并借助时空图像分析方法对鲁甸地震的破裂过程进行了研究，认为本次地震的主破裂面优势方向为NW向，其发震断层为包谷垴—小河断裂的可能性很大。
此前对于鲁甸地震的研究均集中于2014年8月3日鲁甸MS6.5地震的主震破裂过程以及两个月内的余震精确定位（表1），对于2014年9月之后的余震序列研究较少。根据云南测震台网数据显示，鲁甸余震序列迄今仍在持续活动，而且发生了多次晚期强余震（表2）。鉴于此，本文拟使用双差定位方法对2014年8月3日至2017年2月28日期间的鲁甸MS6.5地震余震序列进行重定位，给出鲁甸地震余震序列的时空分布特征，并在此基础上对4次强余震（同一天的发生的一组地震算一次）序列的时空分布及其与主震序列的关系予以分析，结合地质等相关资料研究鲁甸余震序列发展、余震时空分布与震后应力场变化的关系，以期为滇东北地区强震序列发展及震后区域应力场变化研究提供震例依据。
表 2 云南测震台网测定的2014年9月至2017年2月间鲁甸地区强余震的震源参数Table 2. Focal parameters of the strong aftershocks from September 2014 to February 2017 recorded by Yunnan seismic network发震时刻 北纬/° 东经/° 深度/km MS 年−月−日 时：分：秒 2014−09−10 16：59：29 27.11 103.37 9.0 4.7 2014−10−27 00：02：28 27.11 103.33 6.0 4.6 2016−08−12 14：56：54 27.03 103.37 14.0 4.7 2016−08−12 19：25：20 27.03 103.39 7.0 4.7 2017−02−08 19：11：38 27.07 103.35 10.0 4.9 1. 数据和方法
1.1 台站分布和数据
2014年鲁甸MS6.5地震发生后，云南省地震局在震源区附近布设了3套流动台，结合原有的昭通市小孔径测震台网和云南区域地震台网台站，在震源区形成了密集且包围较好的观测台网（赵小艳等，2014）。图1为鲁甸地震震源区附近的台站分布。
![]()
图 1 2014年鲁甸MS6.5地震震源区附近的台站分布Figure 1. The distribution of stations near Ludian 2014 MS6.5 earthquake source region本文采用云南测震台网观测报告数据，选择被4个及以上台站记录到的ML≥1.0地震事件，经处理后共有4 236个地震事件成群。云南测震台网原始地震定位结果如图2所示，可以看出：鲁甸地震序列的震中位置呈共轭型分布，在主地震序列周围有很多小震级余震散乱分布；地震序列震源深度主要分布在4—20 km范围内，主震震源深度为13 km。
![]()
图 2 云南地震台网测定的鲁甸地震序列原始定位结果Figure 2. Initial location of Ludian earthquake sequence by Yunnan seismic network尽管双差定位方法可以消除由震源与台站之间速度模型的不确定性产生的误差，但却非常依赖震源所在处的速度模型，而且在地震定位中震源深度与速度模型相关性较大，因此使用双差定位方法时需要选取恰当的震源区速度模型来提高定位精度（王光明等，2015；王清东等，2015）。王未来等（2014）参考天然地震和人工地震测深的结果构建了鲁甸地区的一维速度模型（表3），本文选取该模型对鲁甸地震序列进行重新定位。
表 3 鲁甸地区速度模型Table 3. Velocity model of Ludian area地壳厚度/km vP/（km•s−1） vS/（km•s−1） 0—2 5.50 3.179 2—8 5.75 3.324 8—16 6.00 3.468 16—24 6.30 3.642 24—32 5.60 3.240 32—46 6.70 3.873 46＜ 7.70 4.451 为了检查观测报告中震相数据的可靠性，本文绘制了P波和S波的震相走时曲线，并且使用所选鲁甸地区速度结构模型绘制了理论走时曲线，如图3所示。从图中可以非常清楚地区别P波与S波震相的走时曲线，并且震相的离散度较小，表明观测报告中震相数据的可靠性较高；观测震相走时与理论走时曲线吻合较好，基本无偏差，证明定位所使用的鲁甸地区速度模型精度较高。根据震相走时曲线判断，鲁甸地区地震记录的可靠性较好，因此本文选取震中距350 km以内的震相数据进行精定位。为提高反演方程的稳定性和计算效率，我们将ph2dt软件控制文件中MAXOBS 参数修改为20，然后使用ph2dt软件对所获得的观测报告中P波和S波到时进行筛选，最终得到22万2 210条P波走时差和20万6 275条波走时差数据。
![]()
图 3 鲁甸地震序列P波和S波的震相走时曲线Figure 3. Travel time curves of P and S phases for the 2014Ludian earthquake sequence1.2 双差定位法
双差定位方法是Waldhauser和Ellsworth （2000）提出的一种基于主事件定位法的相对定位方法，该方法反演的是在一组丛集的地震中每次地震相对于该丛集矩心的相对位置，不需要利用主事件，因此避免了主事件定位法空间跨度小的问题（杨智娴等，2003）。双差定位方法是使用相邻地震事件相对于同一台站的观测走时差与理论走时差之差（双差）来确定震源的位置，可以消除由震源与台站之间速度模型不确定性所产生的误差，并有效地降低由于对地壳速度结构数据不够精确引起的误差。双差定位方法提出后在国内外获得了广泛的应用（杨智娴等，2003；Hauksson，Shearer，2005；黄媛等，2008；Waldhauser，Schaff，2008；陈晨，胥颐，2013；王未来等，2014；叶庆东等，2017）。
2. 精定位结果
2014年鲁甸MS6.5地震序列重定位后获得了3 658次地震的震源参数，地震序列的震中分布和震源深度剖面如图4和图5所示。鲁甸MS6.5地震重定位后的震中位置为（103.35°E，27.11°N），震源深度为19.93 km，发震时刻为北京时间2014年8月3日16时30分10.58秒。
![]()
图 5 重定位后鲁甸地震序列沿AA′，BB′，CC′，DD′，EE′和FF′深度剖面的分布Figure 5. Distribution of relocated Ludian earthquake sequence along the depth profiles ofAA′，BB′，CC′，DD′，EE′ and FF′从图4的震中分布可以看出，重定位后整个地震序列向内紧缩，地震序列的近EW向和SES向两个优势分布方向更加清晰，呈共轭型分布。近EW向余震展布长约30 km、宽约6 km，SES向余震展布长约20 km，宽约8 km，本文重定位后余震展布长度与其它研究结果（表4）相比更长，而且近EW向长于SES向。仔细观察震中（图4）和震源深度（图5中BB′，CC′剖面）分布可以看出，近EW向分支西端长约15 km一带地震数量较少，震中分布散乱，东西向展布形态不明显；震源深度较深，与余震序列向西逐渐变浅的主体趋势相悖。已有研究表明，早期余震主要集中分布在主震破裂面附近，而后期的余震有可能是主震触发的其它断层上的地震（Mendoza，Hartzell，1988；Freed，2005）。因此，我们猜想这一部分地震可能是余震序列向西扩展时触动次级或更小断层产生破裂的结果。
表 4 不同研究给出的余震分布特征Table 4. Aftershock distribution characteristics from different researchs来源 近EW向余震
分布长度/kmSES向余震分布
长度/km王未来等（2014） 15 15 张广伟等（2014） 17 22 徐甫坤等（2014） – 15 房立华等（2014） 16 16 本文结果 30 20 为分析地震序列在深部的展布形态，绘制了6条震源深度剖面，其中：AA′和BB′分别是与近EW向分支和SES向分支平行的两个深度剖面；CC′，DD′与BB′垂直；EE′，FF′与AA′垂直。从图5可以看出：AA′剖面中震源深度沿SES向逐渐变浅，且余震主要分布在主震的东南侧；BB′剖面中震源深度沿近EW向逐渐变浅，余震主要分布在主震的西侧。综合分析AA′和BB′两剖面震源深度分布趋势可知，地震序列总体表现为主震附近震源深度较深，并且沿近EW向和SES向震源深度逐渐变浅，整个地震序列的震源深度主要分布在4—20 km范围内。
垂直于近EW向的剖面CC′和DD′显示出鲁甸地震序列近EW向分支不同段的结果，如图5c，d所示。可以看出：余震在CC′剖面呈明显的两部分，且震源深度优势分布范围为10—20 km；余震在DD′剖面中较为集中分布，与地表垂直，震源深度分布优势范围为4—15 km，与BB′中地震序列主体部分相同。
垂直于SES向的剖面EE′和FF′显示出鲁甸地震序列SES向分支不同段的结果，如图5e，f所示，可以看出：余震序列在两个剖面中震源深度分布形态相似，近乎垂直于地表，表明控制此分支的断层为近纯走滑型，其中EE′剖面震源的优势分布深度处于4—17 km，而FF′剖面震源优势分布深度范围为3—15 km，与AA′剖面结果一致，表明余震序列沿SES向扩展过程中震源深度逐渐变浅。
剖面DD′，EE′和FF′的余震分布表明鲁甸地震的后续余震发生在两条共轭断层上，且这两条断层均为走滑断层。张勇等（2015）通过震源破裂过程反演结果分析认为，鲁甸地震是一次在统一应力场下发生的NW向主压应力和NE向主张应力的两条共轭断层先后破裂的一次复杂地震事件，重定位后余震的空间分布再次验证了该观点。2014年鲁甸MS6.5地震与2011年3月10日盈江MS5.8地震的发震构造相似，Lei等（2012）使用双差定位方法对该地震进行精确定位，由重定位后余震序列空间分布特征推测出盈江地震可能引发了共轭断层系统的破裂。
3. 强余震序列分析
截至2017年2月28日，除鲁甸MS6.5地震之外，震源区共发生4次（同一天发生的一组地震算一次）MS≥4.5强余震，4次强余震发生后30天内地震序列震中分布如图6所示。
![]()
图 6 4次强余震发生后30天内小震震中分布Figure 6. Distribution of earthquakes occurred in 30 days after four large aftershocks由图6a和图6b可以看出，2014年两次强余震的震中分布具有以下3个特征：① 两次地震发生的位置不同，2014年9月10日MS4.7地震发生在两条分支的交界处，距离主震较近，2014年10月27日MS4.6地震发生在近EW向分支的中心位置；② 两次地震发生后30天的余震展布与主震序列相同，均存在近EW向和SES向两条共轭的余震分支；③ 30天内余震数目不同，9月10日MS4.7地震较10月27日MS4.6地震的余震数目更多，且在两条分支的端点处较为集中。上述特征的形成是由于两次强余震距离主震时间较近（38天和85天），余震序列仍然在发育，因此尽管震源位置不同，但是后续30天的余震展布均与主震序列相同；另外由于9月10日MS4.7强余震在时间上距离主震更近，因此其后续余震数量明显多于10月27日MS4.6强余震。
2016年8月12日MS4.7 强余震和2017年2月28日MS4.9强余震发生后30天内，后续余震基本发生在地震周边（图6c，d），而且序列衰减较快，明显不同于2014年两次强余震。2016年8月12日两次MS4.7地震发生在主震序列的SES向分支，在近EW向分支上余震较少，且分布散乱；2017年2月8日MS4.9地震发生在近EW向分支上，而在SES向几乎无余震。
4次强余震序列的震源深度统计如图7所示，2014年两次强余震序列的震源深度均集中在4—16 km范围内，与主震序列深度分布一致；2016年和2017年两次地震的余震较2014年地震明显减少，而且震源深度无明显优势分布，与主震序列深度的分布特征不同。
![]()
图 7 4次强余震发生后30天内震源深度统计Figure 7. Focal depth histogram of earthquakes occurred in 30 days after four large aftershocks对比4次强余震序列的震中分布和深度统计发现，随着发震时刻距离主震越来越远，强余震后地震序列具有以下特征：① 地震数量越来越少；② 震源深度的优势分布越来越深；③ 地震序列越来越集中发生在强余震周围。上述特征表明鲁甸MS6.5地震引发的震源区应力场随着时间的推移逐渐变小。
4. 讨论与结论
4.1 重定位结果讨论
鲁甸地震序列精定位已有较多研究结果，本文将震中分布和误差统计结果进行对比分析，给出本文与前人研究结果的异同。
重定位后鲁甸地震序列呈共轭型分布，存在两条优势分布方向，分别为近EW向和SES向，与前人的研究相同。但是地震序列的余震展布范围与前人的研究存在差异（表4），分析认为造成差异的原因主要有两点：① 不同研究所使用的数据的时间范围不同（表1），时间范围越长，震源区产生的余震数量越多，余震沿近EW向和SES向的展布长度越长；② 重定位后地震序列空间分布（图4，5）显示近EW向分支存在一个长约15 km的余震区，相比于其它区域，该余震区地震数量较少，震中分布散乱，东西向展布形态不明显，震源深度较深，与近EW向分支震源深度向西逐渐变浅的主体趋势相悖。研究表明，早期余震主要集中分布在主震破裂面附近，而后期的余震有可能是由主震触发的其它断层上的地震（Mendoz，Hartzell，1988；Freed，2005）。因此，我们猜想这一部分地震可能是余震序列向西扩展时触动次级或更小断层产生破裂的结果。如果不考虑这一部分地震，重定位后鲁甸地震序列近EW向分支和SES向分支展布长度相同，均为20 km左右，与前人研究结果相符。
表 1 不同研究分析鲁甸MS6.5地震序列所选取的精定位方法和时间范围Table 1. Relocation method and time range of Ludian MS6.5 earthquake sequence from several researchs来源 研究方法 时间范围 王未来等（2014） 双差定位 2014−08−03—2014−08−19 张广伟等（2014） 双差定位 2014−08−03—2014−08−07 徐甫坤等（2014） 双差定位 2014−08−03—2014−08−13 房立华等（2014） 双差定位 2014−08−03—2014−10−03 李艳娥等（2015） 主事件定位 2014−08−03—2014−09−30 ![]()
图 4 鲁甸地震序列重定位后的震中分布Figure 4. The epicentral distribution of Ludian earthquake sequence after relocation不同研究人员给出的重定位误差如表5所示，与前人研究相比，本文的重定位误差较小，但是误差水平为同一量级。对比表中不同方向重定位误差发现，所有研究的水平向误差均小于垂直向误差；王未来等（2014）和房立华等（2014）的近EW向重定位误差大于NS向，而张广伟等（2014）、徐甫坤等（2014）以及本文结果的近EW向重定位误差小于NS向，本文重定位误差形成上述特征可能是由于所使用的台站分布不均匀造成的，震源区NW向无台站分布。
表 5 不同研究给出的重定位误差Table 5. Relocation error from different researchs来源 误差 东西向/km 南北向/km 深度/km 发震时刻/s 王未来等（2014） 0.800 0.712 1.520 0.220 张广伟等（2014） 0.44 0.57 1.23 – 徐甫坤等（2014） 0.80 0.90 1.50 0.33 房立华等（2014） 0.42 0.37 0.74 0.18 本文结果 0.274 0.272 0.623 0.226 4.2 鲁甸地区地质结构与强余震序列分析
鲁甸MS6.5地震震中位于青藏高原东南缘与华南地块的交界部位，属于大凉山次级地块。大凉山次级地块在巴颜喀拉地块、川滇菱形地块和华南地块之间的地壳形变中起到了重要的调节作用，因此该次级地块不仅受到川滇菱形地块在安宁河断裂带上的挤压作用，还受到巴颜喀拉地块在龙门山断裂带上的挤压作用（闻学泽等，2013；徐锡伟等，2014；程佳等，2016）。在这种挤压和NNW向的左旋走滑作用下，大凉山地块内部主要以NNW向的左旋走滑断裂为主，包括大凉山断裂、马边断裂等；这些左旋运动在断裂南端被NE向的昭通—鲁甸断裂带和莲峰断裂带吸收并转化为挤压活动特征（程佳等，2014；徐锡伟等，2014）。
昭通—鲁甸断裂带和莲峰断裂带是大凉山次级地块的南边界，起着吸收和调节凉山次级活动地块SE向运动的作用（闻学泽等，2013；常祖峰等，2014；张彦琪等，2016）。在NW-SE向主压应力与NE-SW向主张应力的统一应力场下，昭通—鲁甸断裂带和莲峰断裂带以逆冲走滑活动为主，鲁甸MS6.5地震发生在两条断裂带之间、NNW走向的包谷垴—小河断裂带上（房立华等，2014；徐锡伟等，2014；张勇等，2015）。张勇等（2015）和程佳等（2016）分别计算了鲁甸地震的震源破裂过程和库仑应力，结果显示鲁甸地震是一次在统一的区域应力场作用下发生在两条共轭断裂上的复杂地震事件。破裂最初发生在近EW向的断层上，之后触发了与其共轭的NNW向断层的破裂，本文重定位后地震序列的震中分布与上述结果相吻合。
本文重定位后的余震震中分布显示，近EW向余震分支延伸范围较大，而SES向余震分布范围基本稳定。结合重定位后的地震序列时空分布、鲁甸地区地质结构和相关研究，我们猜想近EW向分支15 km外的地震可能是余震序列向西扩展时触动次级或更小断层产生破裂的结果，其发震时刻主要集中在主震发生后100—600天之间，后续余震较少；SES向分支属于包谷垴—小河断裂南段，该分支在向南扩展的过程中可能遇到障碍体，这些地方或是较为软弱的区域，难以积累足够的应力形成地震，或者是较为坚硬的区域，目前正处于闭锁积累能量的过程（徐甫坤等，2014）。由于目前掌握的资料较少，我们认为其未来的地震危险性需要更进一步探究。
王清东（2015）对盈江、腾冲、景谷等地震序列重定位后得出，地震序列的优势分布方向与附近的断裂和区域应力场有关，但是大的地震会破坏区域应力场和附近的地质结构，使得地震序列呈现新的优势分布方向，因此，余震序列的分布同样受到由主震或大的余震引发的震源区应力场控制。刘蒲雄等（1996）认为，超晚期强余震是在区域应力增强的背景下发生的，因此超晚期强余震不宜看作是震源区剩余能量的释放，主要能量似乎是由外部补给的。
结合上述研究分析鲁甸地震的4次强余震序列，可以得出以下结论：鲁甸MS6.5主震的发生破坏了原有的区域应力场，改变了震源区附近的地质结构，形成了震源区应力场；由于2014年9月和10月两次强余震的发震的时刻距离主震较近，后续余震仍然主要受由主震引发的震源区应力场的控制，从而形成尽管震源位置不同、但后续余震分布相同的现象；随着时间的推移，主震引发的震源区应力场减弱，2016年8月和2017年2月发生的强余震可能是由于鲁甸地区的区域应力增强所引发的，后续余震分布主要受区域应力场和附近断层的控制，而这两次强余震分别发生在鲁甸震源区两条共轭的断层上，因此其后续余震分布特征不同。
4.3 结论
本文采用双差定位方法对2014年鲁甸MS6.5地震序列进行重定位，主要结论如下：
1） 鲁甸地震序列存在两个优势分布方向，分别为近EW向和SES向，呈共轭型分布；近EW向余震展布长约30 km，宽约6 km，SES向余震展布长约20 km，宽约8 km；鲁甸地震序列总体表现为主震附近震源深度较深，沿近EW向和SES向两个分支震源深度逐渐变浅，整个地震序列震源深度主要分布在4—20 km范围内。
2） 鲁甸地震序列近EW向分支西段存在一部分震中分布散乱、震源深度较深的地震，与该分支主体分布明显不同，结合前人研究猜想这部分地震可能是余震序列向西扩展时触动次级或更小断层产生破裂的结果。
3） 鲁甸地震的4次强余震序列是区域应力场和主震引发的震源区应力场共同作用的结果，其中2014年9月和10月的两次强余震序列主要受到主震引发的震源区应力场的影响，而2016年和2017年两次强余震序列主要受到区域应力场的影响。
拉蒙特-多尔蒂地球观测站的Felix Waldhuaser教授提供了双差定位程序，云南省测震台网及昭通市测震台网提供了观测报告，审稿专家提出了宝贵意见，作者在此一并表示感谢。
-
![]()
图 1 2014年鲁甸MS6.5地震震源区附近的台站分布
Figure 1. The distribution of stations near Ludian 2014 MS6.5 earthquake source region
![]()
图 2 云南地震台网测定的鲁甸地震序列原始定位结果
Figure 2. Initial location of Ludian earthquake sequence by Yunnan seismic network
![]()
图 3 鲁甸地震序列P波和S波的震相走时曲线
Figure 3. Travel time curves of P and S phases for the 2014Ludian earthquake sequence
![]()
图 5 重定位后鲁甸地震序列沿AA′，BB′，CC′，DD′，EE′和FF′深度剖面的分布
Figure 5. Distribution of relocated Ludian earthquake sequence along the depth profiles ofAA′，BB′，CC′，DD′，EE′ and FF′
![]()
图 6 4次强余震发生后30天内小震震中分布
Figure 6. Distribution of earthquakes occurred in 30 days after four large aftershocks
![]()
图 7 4次强余震发生后30天内震源深度统计
Figure 7. Focal depth histogram of earthquakes occurred in 30 days after four large aftershocks
![]()
图 4 鲁甸地震序列重定位后的震中分布
Figure 4. The epicentral distribution of Ludian earthquake sequence after relocation
表 2 云南测震台网测定的2014年9月至2017年2月间鲁甸地区强余震的震源参数
Table 2 Focal parameters of the strong aftershocks from September 2014 to February 2017 recorded by Yunnan seismic network
发震时刻 北纬/° 东经/° 深度/km MS 年−月−日 时：分：秒 2014−09−10 16：59：29 27.11 103.37 9.0 4.7 2014−10−27 00：02：28 27.11 103.33 6.0 4.6 2016−08−12 14：56：54 27.03 103.37 14.0 4.7 2016−08−12 19：25：20 27.03 103.39 7.0 4.7 2017−02−08 19：11：38 27.07 103.35 10.0 4.9 
下载: 导出CSV
表 3 鲁甸地区速度模型
Table 3 Velocity model of Ludian area
地壳厚度/km vP/（km•s−1） vS/（km•s−1） 0—2 5.50 3.179 2—8 5.75 3.324 8—16 6.00 3.468 16—24 6.30 3.642 24—32 5.60 3.240 32—46 6.70 3.873 46＜ 7.70 4.451
下载: 导出CSV
表 4 不同研究给出的余震分布特征
Table 4 Aftershock distribution characteristics from different researchs
来源 近EW向余震
分布长度/kmSES向余震分布
长度/km王未来等（2014） 15 15 张广伟等（2014） 17 22 徐甫坤等（2014） – 15 房立华等（2014） 16 16 本文结果 30 20
下载: 导出CSV
表 1 不同研究分析鲁甸MS6.5地震序列所选取的精定位方法和时间范围
Table 1 Relocation method and time range of Ludian MS6.5 earthquake sequence from several researchs
来源 研究方法 时间范围 王未来等（2014） 双差定位 2014−08−03—2014−08−19 张广伟等（2014） 双差定位 2014−08−03—2014−08−07 徐甫坤等（2014） 双差定位 2014−08−03—2014−08−13 房立华等（2014） 双差定位 2014−08−03—2014−10−03 李艳娥等（2015） 主事件定位 2014−08−03—2014−09−30
下载: 导出CSV
表 5 不同研究给出的重定位误差
Table 5 Relocation error from different researchs
来源 误差 东西向/km 南北向/km 深度/km 发震时刻/s 王未来等（2014） 0.800 0.712 1.520 0.220 张广伟等（2014） 0.44 0.57 1.23 – 徐甫坤等（2014） 0.80 0.90 1.50 0.33 房立华等（2014） 0.42 0.37 0.74 0.18 本文结果 0.274 0.272 0.623 0.226
下载: 导出CSV
-
常祖峰,周荣军,安晓文,陈宇军,周青云,李鉴林. 2014. 昭通—鲁甸断裂晚第四纪活动及其构造意义[J]. 地震地质,36(4):1260–1279 Chang Z F,Zhou R J,An X W,Chen Y J,Zhou Q Y,Li J L. 2014. Late-Quaternary activity of the Zhaotong-Ludian fault zone and its tectonic implication[J]. Seismology and Geology,36(4):1260–1279 (in Chinese)
陈晨,胥颐. 2013. 芦山MS7.0级地震余震序列重新定位及构造意义[J]. 地球物理学报,56(12):4028–4036 Chen C,Xu Y. 2013. Relocation of the Lushan MS7.0 earthquake sequence and its tectonic implication[J]. Chinese Journal of Geophysics,56(12):4028–4036 (in Chinese)
程佳,刘杰,徐锡伟,甘卫军. 2014. 大凉山次级地块内强震发生的构造特征与2014年鲁甸6.5级地震对周边断层的影响[J]. 地震地质,36(4):1228–1243 Cheng J,Liu J,Xu X W,Gan W J. 2014. Tectonic characteristics of strong earthquakes in Daliangshan sub-block and impact of the MS6.5 Ludian earthquake in 2014 on the surrounding faults[J]. Seismology and Geology,36(4):1228–1243 (in Chinese)
程佳,徐锡伟,刘杰. 2016. 2014年鲁甸6.5级地震成因、破裂特征和余震分布特征的库仑应力作用[J]. 地球物理学报,59(2):556–567 Cheng J,Xu X W,Liu J. 2016. Cause and rupture characteristics of the 2014 Ludian MS6.5 mainshock and its aftershock distribution using the Coulomb stress changes[J]. Chinese Journal of Geophysics,59(2):556–567 (in Chinese)
房立华,吴建平,王未来,吕作勇,王长在,杨婷,钟世军. 2014. 云南鲁甸MS6.5地震余震重定位及其发震构造[J]. 地震地质,36(4):1173–1185 Fang L H,Wu J P,Wang W L,Lü Z Y,Wang C Z,Yang T,Zhong S J. 2014. Relocation of the aftershock sequence of the MS6.5 Ludian earthquake and its seismogenic structure[J]. Seismology and Geology,36(4):1173–1185 (in Chinese)
何骁慧,倪四道,刘杰. 2015. 2014年8月3日云南鲁甸M6.5地震破裂方向性研究[J]. 中国科学: 地球科学,45(3):253–263 He X H,Ni S D,Liu J. 2015. Rupture directivity of the August 3rd, 2014 Ludian earthquake (Yunnan, China)[J]. Science China Earth Sciences,58(5):795–804
黄媛. 2008. 结合波形互相关技术的双差算法在地震定位中的应用探讨[J]. 国际地震动态,(4):29–34 Huang Y. 2008. Study on the application and development of the DD algorithm with cross correlation of waveform data in the earthquake location[J]. Recent Developments in World Seismology,(4):29–34 (in Chinese)
刘蒲雄,陈修启,吕晓健,韩丹. 1996. 地震序列的后续显著地震的预测研究[J]. 地震学报,18(1):27–33 Liu P X,Chen X Q,Lü X J,Han D. 1996. Prediction Research of significant aftershocks of earthquake sequence[J]. Acta Seismologica Sinica,18(1):27–33 (in Chinese)
李艳娥,陈学忠,陈丽娟,郭祥云. 2015. 2014年8月3日云南鲁甸6.5级地震序列破裂过程研究[J]. 地球物理学报,58(9):3232–3238 Li Y E,Chen X Z,Chen L J,Guo X Y. 2015. Investigation on the progress of the Ludian MS6.5 earthquake sequence on 3 August, 2014 in Yunnan Province[J]. Chinese Journal of Geophysics,58(9):3232–3238 (in Chinese)
王光明,朱良宝,苏有锦,王清东. 2015. 2012年6月24日宁蒗—盐源5.7级地震精确定位[J]. 地震学报,37(5):733–746, 885 Wang G M,Zhu L B,Su Y J,Wang Q D. 2015. Accurate locating of the June 24, 2012 Ninglang-Yanyuan MS5.7 earthquake[J]. Acta Seismologica Sinica,37(5):733–746, 885 (in Chinese)
王清东,朱良宝,苏有锦,王光明. 2015. 2012年9月7日彝良地震及余震序列双差定位研究[J]. 地球物理学报,58(9):3205–3221 Wang Q D,Zhu L B,Su Y J,Wang G M. 2015. Double-difference relocation of the 7 September 2012 Yiliang earthquake and its aftershock sequence[J]. Chinese Journal of Geophysics,58(9):3205–3221 (in Chinese)
王清东. 2015. 利用双差定位法研究云南地震分布的活动构造意义[D]. 武汉: 武汉大学: 122–123. Wang Q D. 2015. Using Double-Difference Location Method to Research the Active Tectonic Significance of Yunnan Earthquakes Distribution[D].Wuhan: Wuhan University: 122–123 （in Chinese）.
王未来,吴建平,房立华,来贵娟. 2014. 2014年云南鲁甸MS6.5地震序列的双差定位[J]. 地球物理学报,57(9):3042–3051 Wang W L,Wu J P,Fang L H,Lai G J. 2014. Double difference location of the Ludian MS6.5 earthquake sequences in Yunnan Province in 2014[J]. Chinese Journal of Geophysics,57(9):3042–3051 (in Chinese)
闻学泽,杜方,易桂喜,龙锋,范军,杨攀新,熊仁伟,刘晓霞,刘琦. 2013. 川滇交界东段昭通、莲峰断裂带的地震危险背景[J]. 地球物理学报,56(10):3361–3372 Wen X Z,Du F,Yi G X,Long F,Fan J,Yang P X,Xiong R W,Liu X X,Liu Q. 2013. Earthquake potential of the Zhaotong and Lianfeng fault zones of the eastern Sichuan-Yunnan border region[J]. Chinese Journal of Geophysics,56(10):3361–3372 (in Chinese)
徐甫坤,李静,苏有锦. 2014. 2014年云南鲁甸6.5级地震序列重定位研究[J]. 地震研究,37(4):515–522 Xu F K,Li J,Su Y J. 2014. Relocations of Yunnan Ludian MS6.5 earthquake sequence in 2014[J]. Journal of Seismological Research,37(4):515–522 (in Chinese)
徐锡伟,江国焰,于贵华,吴熙彦,张建国,李西. 2014. 鲁甸6.5级地震发震断层判定及其构造属性讨论[J]. 地球物理学报,57(9):3060–3068 Xu X W,Jiang G Y,Yu G H,Wu X Y,Zhang J G,Li X. 2014. Discussion on seismogenic fault of the Ludian MS6.5 earthquake and its tectonic attribution[J]. Chinese Journal of Geophysics,57(9):3060–3068 (in Chinese)
杨智娴,陈运泰,郑月军,于湘伟. 2003. 双差地震定位法在我国中西部地区地震精确定位中的应用[J]. 中国科学:D辑,33(增刊1):129–134 Yang Z X,Chen Y T,Zheng Y J,Yu X W. 2003. Accurate relocation of earthquakes in central-western China using the double-difference earthquake location algorithm[J]. Science in China:Series D,46(S2):181–188
叶庆东,丁志峰,王生文,余大新,郑晨. 2017. 汶川地震科学钻探3号井孔附近微震震源参数确定及其意义[J]. 地球物理学报,60(7):2716–2732 Ye Q D,Ding Z F,Wang S W,Yu D X,Zheng C. 2017. Determining the source parameters of the microearthquakes near the third borehole of the Wenchuan earthquake fault scientific drilling (WFSD-3) and its implications[J]. Chinese Journal of Geophysics,60(7):2716–2732 (in Chinese)
张广伟,雷建设,梁姗姗,孙长青. 2014. 2014年8月3日云南鲁甸MS6.5级地震序列重定位与震源机制研究[J]. 地球物理学报,57(9):3018–3027 Zhang G W,Lei J S,Liang S S,Sun C Q. 2014. Relocation and focal mechanism solutions of the 3 August 2014 Ludian, Yunnan MS6.5 earthquake sequence[J]. Chinese Journal of Geophysics,57(9):3018–3027 (in Chinese)
张彦琪,李西,谢英情,于江,陈坤华. 2016. 以鲁甸地震和彝良地震为例分析云南昭通地质构造与地震灾害特点[J]. 地震研究,39(2):270–278 Zhang Y Q,Li X,Xie Y Q,Yu J,Chen K H. 2016. Analysis on seismotectonic background and earthquake hazard characteristic in Zhaotong, Yunnan: Taking Ludian MS6.5 earthquake and Yiliang MS5.7,5.6 earthquake as examples[J]. Journal of Seismological Research,39(2):270–278 (in Chinese)
张勇,许力生,陈运泰,刘瑞丰. 2014. 2014年8月3日云南鲁甸MW6.1 (MS6.5)地震破裂过程[J]. 地球物理学报,57(9):3052–3059 Zhang Y,Xu L S,Chen Y T,Liu R F. 2014. Rupture progress of the 3 August 2014 Ludian, Yunnan, MW6.1 (MS6.5) earthquake[J]. Chinese Journal of Geophysics,57(9):3052–3059 (in Chinese)
张勇,陈运泰,许力生,魏星,金明培,张森. 2015. 2014年云南鲁甸MW6.1地震: 一次共轭破裂地震[J]. 地球物理学报,58(1):153–162 Zhang Y,Chen Y T,Xu L S,Wei X,Jin M P,Zhang S. 2015. The 2014 MW6.1 Ludian, Yunnan, earthquake: A complex conjugated ruptured earthquake[J]. Chinese Journal of Geophysics,58(1):153–162 (in Chinese)
赵小艳,韩立波,徐甫坤,刘自凤,孙楠,李永莉. 2014. 2014年云南鲁甸6.5级地震序列跟踪分析研究[J]. 地震研究,37(4):508–514 Zhao X Y,Han L B,Xu F K,Liu Z F,Sun N,Li Y L. 2014. Relocations of Yunnan Ludian MS6.5 earthquake sequences in 2014[J]. Journal of Seismological Research,37(4):508–514 (in Chinese)
Freed A M. 2005. Earthquake triggering by static, dynamic, and postseismic stress transfer[J]. Annu Rev Earth Planet Sci,33(1):335–367
Hauksson,E,Shearer,P. 2005. Southern California hypocenter relocation with waveform cross-correlation, part 1: Results using the double-difference method[J]. Bull Seismol Soc Am,95(3):896–903
Lei J S,Zhang G W,Xie F R,Li Y,Su Y J,Liu L F,Ma H H,Zhang J W. 2012. Relocation of the 10 March 2011 Yingjiang, China, earthquake sequence and it’s tectonic implications[J]. Earthquake Science,25(1):103–110
Mendoza C,Hartzell S H. 1988. Aftershock patterns and mainshock faulting[J]. Bull Seismol Soc Am,78(4):1438–1449
Waldhauser F,Ellsworth W L. 2000. A double-difference earthquake location algorithm: Method and application to the northern Hayward fault, California[J]. Bull Seismol Soc Am,99(6):1353–1368
Waldhauser F,Schaff D P. 2008. Large-scale relocation of two decades of Northern California seismicity using cross-correlation and double-difference methods[J]. J Geophys Res,113(B8):B08311
-
期刊类型引用(7)
1. 张演，姜金钟，王光明，张帅，何德强，周青云. 云南盐津地区两种地震序列的地震活动性分析. 地震研究. 2024(02): 178-190 . 
百度学术
2. 宋倩，梁姗姗，于湘伟，章文波，李春来. 云南2020年巧家M_S5.0地震序列发震构造及其与2014年鲁甸M_S6.5地震的关系. 地球物理学报. 2022(04): 1303-1324 . 百度学术
3. 曹颖，黄江培，钱佳威，付虹. 利用时移层析成像方法揭示与2014年云南鲁甸M_S6.5地震有关的P波速度变化. 地球物理学报. 2021(05): 1569-1584 . 百度学术
4. 孙楠，付虹，彭关灵，刘自凤. 云南地区强震前小震调制比特征分析. 地震地磁观测与研究. 2021(05): 82-89 . 百度学术
5. 杨帆，盛书中，胡晓辉，崔华伟. 2020年云南巧家M5.0地震震源机制解与震源区应力场研究. 地震. 2021(04): 93-105 . 百度学术
6. 李姣，姜金钟，杨晶琼. 2017年漾濞M_S4.8和M_S5.1地震序列的微震检测及重定位. 地震学报. 2020(05): 527-542+507 . 本站查看
7. 鞠慧超，张元生，刘旭宙. 2019年甘肃张掖M5.0地震余震序列重定位. 地震工程学报. 2019(06): 1396-1401 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 1557
- HTML全文浏览量: 748
- PDF下载量: 93
- 被引次数: 9
