0 引言

据中国地震台网中心(CENC)正式测定，北京时间2020年5月18日21时47分云南昭通市巧家县发生M_S5.0地震，震中位于27.18°N、103.16°E，震源深度8km。此次地震震级不大，但由于震源较浅，云南昭通、四川宜宾与西昌等地震感明显，造成了一定程度的人员伤亡和财产损失。

巧家M_S5.0地震发生于川滇菱形块体东缘的南北地震带中南段的NE向早更新世活动的莲峰断裂和NE向昭通-鲁甸断裂之间(徐锡伟等，2014)(图 1(a))。这2条断裂带属于川滇块体-华南地块边界带，也是大凉山次级块体与华南地块之间的边界带。由于受印度板块NE向俯冲碰撞的影响，研究区域内构造运动较强烈，中强地震活跃。该区内曾发生1974年5月11日云南永善大关M_S7.1地震。同时，自2000以来，沿NE向昭通-鲁甸地震断裂附近地震活动极其活跃。例如，先后发生2003年鲁甸M_S5.0和M_S5.1地震、2004年鲁甸M_S5.6地震，以及距离此次震中东南侧不到20km的2014年鲁甸M_S6.5地震。

① https://www.eq-igl.ac.cn/

巧家M_S5.0地震发生后，震源区余震活动持续不断。据云南省区域地震台网测定，震后24h内，震源区及邻近区域(27°N~27.3°N，103°E~103.3°E)共记录到M_L≥1.8余震48个，其中最大震级为M_L3.8。然而，受背景噪声干扰、波形相互叠加等因素影响，传统的地震识别方法会遗漏较多的余震事件，特别是主震发生后数分钟到数小时内，常规地震目录遗漏的事件较多。

模板匹配定位方法是一种有效识别小微震事件的可靠方法(Zhang et al，2015)，基于较为完备的地震序列的时空分布形态，其可以直接反映发震断层的深部构造及主震发震机制等。因此，本文使用巧家M_S5.0地震震源区附近台站记录的连续波形数据，通过模板匹配定位方法，检测主震发震后24h内常规地震目录遗漏的地震事件，研究巧家地震震后24h序列的特征和发震构造，并进一步探讨该地震与2014年鲁甸M_S6.5地震之间的关联以及发震断层与周围断裂带之间的相互作用。

1 数据与方法

本研究使用巧家地震发生后24h内、云南和四川区域地震台网记录的震中距在1°以内的5个台站的连续波形资料，进行余震的模板匹配识别，所使用地震台站分布如图 1(b)所示。参考云南省地震局提供的快报地震目录，选取主震后M_L≥1.8的48个地震事件波形作为模板(图 1(b))。

本研究采用张淼(2015)提出的匹配定位(Match and Locate，简称M & L)方法，对常规地震目录中遗漏的地震事件进行微地震检测，该方法利用地震模板事件波形与可能的微震信号做互相关叠加来检测小微地震事件。在波形叠加之前，需要对模板地震周围的三维空间网格进行扫描搜索，计算可能的地震位置与参考位置之间在同一个台站上的相对走时差，根据走时差来对互相关波形进行校正叠加(张淼，2015)。因此，相比传统的模板检测方法，M&L方法不仅对微震检测更加有效，而且能同时获得高精度的地震位置信息，此外还可以检测到距离模板事件较远的地震事件。

检测之初，需要对连续波形进行去均值、去线性趋势等预处理。进行微震检测时，地震波的滤波频段设为2~8Hz；以参考模板事件震源位置为中心，在三维网格空间内搜索微震位置，经度和纬度方向的搜索范围为0.05°，搜索间隔为0.01°；深度方向的搜索范围为2km，搜索间距为0.1km；计算模板事件与每个可能微震记录的走时差，然后按照走时差对互相关波形进行走时校正并叠加。

模板事件波形和连续波形滑动互相关时，一般选取振幅最大的震相作为参考震相，如地方震的Sg震相。在模板匹配中，基于张广伟等(2014)对2014年鲁甸地震序列进行地震定位和震源机制所使用的同一区域速度模型，计算Sg的理论到时和慢度参数，以Sg波理论到时的前1s和后3s作为模板波形的互相关窗口。

在微震检测过程中，以叠加相关波形的平均互相关系数(CC，CC₁＞CC₂)和信噪比(SNR)作为地震检测的标准：当求得的平均互相关系数(CC)大于CC₁时，则仅用互相关系数作为判定标准；当求得的平均互相关系数(CC)大于CC₂而小于CC₁时，则采用互相关系数和信噪比联合约束作为判定标准。在本研究中，互相关系数CC₁、CC₂和波形信噪比经验性地设置为0.35、0.3和10倍背景相关系数(陈安国，2018；郭铁龙等，2020)。当检测到微震事件时，程序自动根据选取的台站分量中参考震相的振幅比中位数来确定微震震级(Peng et al，2009；Meng et al，2013；张淼，2015)。地震事件的检测实例如图 2所示，图中以2020年5月19日12时33分发生的2.5级地震事件为模板，检测到2020年5月18日21时58分发生的1.8级地震事件，地震模板事件波形与检测的微震信号叠加平均相关系数为0.4711，说明各个台站地震波形与2.5级模板事件的相关性较好，意味着检测到一次微地震。

最小完整震级(M_c)可以用来表征地震目录的完整程度，是评估地震台网检测能力的一项重要指标。M_c越小，地震目录越完整，地震台网监测能力越强。相关研究表明，利用模板匹配技术能有效降低目录的最小完整震级(Schaff，2008；张广伟，2016；韩佳东等，2019；徐志国等，2020)。本文基于台网常规识别地震目录和微震检测遗漏地震目录，采用最大曲率法分别计算其最小完整震级M_c(Wiemer et al，2000)，以此来评估2个地震目录的完整性。对于b值的计算，采用ZMAP程序包(Wiemer et al，2000)并基于最大似然法(Aki，1965；Shi et al，1982；Marzocchi et al，2003)得到，其中，误差估计为$\delta b = 2.3{b^2}\sqrt {\frac{{\sum_{i = 1}^n {{{\left({{M_i} - \bar M} \right)}^2}} }}{{n(n - 1)}}} $，式中，n为样本个数，M_i为第i个地震事件的震级，M为平均震级。随后，结合台网目录和微震检测目录，采用震级大于最小完整震级的目录，进一步分析巧家M_S5.0地震震后24h的地震序列活动特征。

2 结果分析 2.1 地震事件检测

利用模板匹配技术对巧家M_S5.0地震后24h内的遗漏地震进行检测和识别，共检测到327个微震事件(不包含48次自检测事件)。据云南省地震局提供的常规地震分析目录，震后24h内相同震源区范围内共发生147个地震事件，其中可定位地震事件134个，单台事件13个，微震模板匹配检测技术所确定的事件个数约为常规地震目录(可定位个数)的2.4倍。

为评估匹配检测所获得的地震目录质量，我们分析了新旧目录的完整性，检测后地震目录的最小完整震级由台网目录的M_L1.9降至检测目录的M_L1.1(图 3(a)、3(b))，表明模板匹配技术可以有效提高区域台网实际监测能力。随后，以各自的最小完整震级为起始震级分别计算b值，获得检测前后巧家地震序列的b值，检测前b值为1.04，误差为±0.21，而检测后b值为0.77，误差为±0.05。比较遗漏事件地震补充前后的b值结果，检测后的b值与误差均明显降低，b值的降低也说明检测后的地震目录小震级事件更多，计算的b值更为可靠。

2.2 地震序列活动特征

为了更好地展示巧家地震震源区的地震序列活动特征，我们给出了微震检测后震源区M-t分布、地震深度分布、深度统计以及震源深度与时间的关系(图 4)。从时间分布(图 4(a))可以看出，巧家地震序列总体呈现衰减趋势，无论是震级大小还是频度均相应降低，2020年5月19日4点35分发生了一次M_L3.8地震，使得余震水平略有增加；从深度分布(图 4(b)~4(d))可以看出，余震震源深度集中在5~15km，平均震源深度为9.3km，震后24h内余震震源深度的每小时平均值由9.8km减至5.6km，总体呈现出由深变浅的趋势。

② https://www.cenc.ac.cn

2.3 发震构造分析

为了详细地展示巧家M_S5.0地震发震断层面的结构特征，我们给出了2条剖面，分别为沿余震展布的SSE方向剖面(170°，纵剖面AA′)和垂直方向剖面(纵剖面BB′)，并以本次主震地震发震时间为参考零点，用非等间距时间段来显示地震序列空间分布，如图 5所示。由图 5可以看出，在时间尺度上，余震分布未看出明显的扩展过程，这可能与本次地震震级规模较小有关；在空间分布上，主震NNW向和SSE向余震分布较多，且分布较广，沿余震展布方向，AA′扩展方向地震延展范围约14km(图 5 (a)、5(b))，震源深度集中在3~15km，属于中上地壳(图 5 (c))；穿过主震且近似垂直于余震扩展方向的纵剖面BB′显示出主震由深至浅的余震活动较为连续，发震断层面较陡，呈略向西倾的趋势(图 5 (c)。为了定量研究发震断层的倾角，将模拟退火全局搜索算法和高斯牛顿局部搜索算法相结合，拟合发震断层面倾角(万永革等，2008)，结果表明巧家地震发震断层倾角约81°，与震源机制反演所得到的2个节面倾角75°和73°均较为接近。结合主震震源机制解结果，推断此次地震的发震断层为NNW-SSE向走滑性质的断裂。本研究给出的推断与云南省地震局根据现场调查^③得到的震源区地震烈度结果有较好的一致性。

③ https://www.yndzj.gov.cn

为进一步了解余震的扩展模式，将震后2h内的地震在平面及深度上的分布展示出来，如图 6所示，由图可见主震发生后0~1h，余震多集中在主震周围；主震发生后1~2h，余震扩展范围增大，主震NNW向和SSE向两侧均有余震分布，且NNW侧余震明显多于SSE侧。

对比巧家M_S5.0地震和2014年8月3日鲁甸M_S6.5地震的余震序列分布(图 7)可以看出，2次地震序列分布均展示出发震断层具有高倾角的特征，有所不同的是巧家余震优势展布方向为NW向，呈现出单侧扩散模式，而在主震区东南方向十余千米处的2014年鲁甸地震主震区的余震优势展布方向为NW向和NE向，呈明显的共轭性破裂模式。前人结合断裂构造性质、主震震源机制和余震时空分布特征等综合认为，2014年鲁甸M_S6.5地震的发震断层为NW向的包谷垴-小河断裂，该断裂为NE向昭通-鲁甸地震的次级断裂，同属小江断裂系(房立华等，2014；张广伟等，2014；张勇等，2015；徐锡伟等，2014)。

根据震源区余震展布特征和震源机制解结果，认为巧家地震的发震断层为NNW向的走滑断裂，这表明鲁甸M_S6.5地震和巧家M_S5.0地震的发震断层均为NW向走滑断裂，且此次巧家地震和鲁甸地震均发生在小江断裂、昭通-鲁甸断裂和莲峰断裂所围限的块体内。由于2次地震的震源位置较为接近，故认为这2次地震有着紧密的联系。从前人给出的地震波速结构来看，巧家地震和鲁甸地震均发生在地壳厚度和泊松比变化较为剧烈的区域(王兴臣等，2015)，且震中区域地壳内存在高、低速相间的多层层状异常结构(李永华等，2014；李大虎等，2019)，而巧家地震的震源体低速异常显著，深度可达15km(张娜等，2019)。特有的深部结构可能为此次巧家地震的发生提供了孕震环境。因此，我们认为巧家地震的发生机制与2014年8月3日云南鲁甸M_S6.5和2011年云南盈江M_S5.8地震相同(杨婷等，2016)，其与印度板块和欧亚板块相互作用形成的高温、高压深部流体作用于断裂带密切相关。值得注意的是，巧家地震序列并未显现出像鲁甸地震序列和盈江地震序列触发共轭断层的特征，可能与本研究所用资料的时间较短或者震级规模较小有关，需要进一步增加震后余震序列资料来证实。

3 结论

本文使用匹配定位方法对2020年5月18日云南巧家M_S5.0主震后24h的波形进行遗漏事件检测，总共检测出327个新的地震事件，约为台网目录的2.4倍。检测后，地震目录的完整性有了明显改善，最小完整震级由最初的M_L1.9降至M_L1.1，从而完善了巧家地震序列目录，为该地区的地震活动性和构造分析提供了丰富的微震数据。通过对检测后的地震事件进行活动性分析，初步分析了此次地震的发震构造，获得如下认识：

(1) 巧家地震序列总体呈现衰减趋势，无论是地震大小还是频度均相应降低；从深度分布可以看出，主震后24h内余震震源深度总体呈现出由深变浅的趋势。

(2) 余震在平面上显示出NNW-SSE向优势展布方向，长度约14km，震源优势深度集中在3~15km之间；深度剖面展示出主震的发震断层面较陡，并且具有向西倾的趋势。由于震源区内台站分布比较稀疏，关于此次地震的可能发震断层面，仍需要结合其他资料进一步分析。

(3) 此次巧家地震发生在小江断裂、昭通-鲁甸断裂和莲峰断裂所围限的块体内，是自2014年鲁甸6.5级地震以来在该区域发生的最大一次地震，2次地震的震中相距仅19.9km，从震源位置、震源机制类型、余震展布形态和速度结构等方面，我们均认为这2次地震关系密切，并认为本次地震的发生与印度板块和欧亚板块相互作用形成的高温、高压深部流体作用于断裂带密切相关。关于此次巧家地震和2014年鲁甸地震之间的关系还有待进一步深入研究。