地震目录是区域地震活动性分析最重要的基础资料，目录完整性如何直接关系到分析结果的可靠性和准确性。由于受背景噪声水平、记录仪器分布等多种因素的影响，许多微弱的地震信号被淹没在噪声中，无法直接观测。另一方面，大地震发生后的短时间内，由于受强震尾波的影响以及大量余震事件的发生，波形相互重叠，难以区分。如果能够自动识别出这些微弱的地震信号和遗漏的早期余震，对于完善地震目录、全面分析地震活动性、监测地震破裂区域的震后形变等具有重要意义^[1]。目前主要采用两种方法自动识别地震信号，一种是基于震相拾取的方法，其中最常用的为长短时间平均法（STA/LTA）^[2-3]，此类方法在数据信噪比较高时能够可靠准确地识别震相^[4]；另一种方法是基于互相关的模板匹配滤波技术^[5-7]，该方法基于滑动窗互相关检测技术，能有效地压制噪声干扰，在数据信噪比低时仍然能够提取微弱信号^[8]。目前模板匹配技术在各个领域都有应用，如火山地震的检测、水压致裂过程中诱发的地震、地震活动性研究等。在地震活动性研究方面，大多数研究者针对震群开展遗漏地震检测，完善地震目录^[8-17]，针对区域地震进行遗漏地震检测的研究不多，张雅楠使用模板匹配定位技术对上海及邻区10年的连续资料进行遗漏地震事件的扫描和检测，共识别出824个地震事件^[18]。本文将利用模板匹配滤波技术对寻乌及邻区2009—2020年的遗漏地震事件进行检测。

赣南地区是一个中强地震多发区，历史上破坏性地震频发。寻乌及邻区是赣南地区地震活动最为强烈的地区，有史以来共发生4.7级以上地震5次，最大地震为1806年会昌南6.0级地震，该地震也是江西境内有史以来记载的最大地震。寻乌及邻区小震一直较为活跃，2009—2020年共记录到地震760次。期间发生了3次小震群活动，分别为2016年7月平远M_L2.7震群、2017年11月寻乌M_L3.8震群和2018年3月安远M_L2.1震群。本文利用模板匹配方法检测寻乌及邻区（包括区域小震群）的遗漏地震，分析地震活动性参数变化，可为区域地震危险性分析和地震趋势判定提供有效的依据。

寻乌及邻区主要的控震断裂带为石城—寻乌断裂带、全南—寻乌断裂带。该地区2009年1月—2020年12月发生M_L≥1.0地震383次，其中M_L2.0～2.9地震58次，M_L3.0以上地震10次。选取距离研究区域最近的寻乌台（XUW）、安远台（ANY）和会昌台（HUC）在2009年1月—2020年12月记录到的连续波形资料，利用M_L2.0以上地震作为模板进行遗漏地震检测。资料来源于江西省测震台网。研究区地震和台站分布见图1。

采用模板匹配识别方法^[7]进行遗漏地震检测。在模板选取过程中，首先选M_L≥1.5地震进行试算，发现研究区域M_L<2.0的模板信噪比不高，M_L1.5～1.9模板很难检测到可靠的遗漏地震，最终选取68次M_L2.0以上地震作为模板，截取寻乌台（XUW）、安远台（ANY）和会昌台（HUC）三分量直达S波到时前2 s到后2 s的波形，经过零相移4阶Butterworth滤波器2～8 Hz滤波后，在连续波形上进行互相关扫描。遗漏地震检测过程分为4步：①将3个台站扫描得到的互相关系数序列按模板地震在此台站的走时平移；②计算3个台站波形互相关系数的平均数，取9倍绝对离差中位数（median absolute devi⁃ation，MAD）为阈值，挑选大于阈值或单台互相关系数大于0.6，而目录中未记录的作为疑似遗漏地震事件；③基于波形互相关震相检测技术检测疑似遗漏地震事件的震相，在距离震中最近的3个台至少检测到一个P波和S波震相的则确认为遗漏地震，扫描得到的互相关序列中互相关系数最大值位置即为遗漏事件P波、S波到时（图2）；④利用遗漏事件水平向波形S波段最大振幅与模板地震的振幅比估计遗漏事件的震级^[8，14]。

本文利用匹配滤波技术对寻乌及邻区2009年1月—2020年12月的遗漏地震事件进行检测，共检测到306次遗漏地震事件，约为台网目录的40％，与胡涛等人^[17]对青海茫崖5.0级地震震源区微震检测得到的遗漏地震比例43%相近，也与尹欣欣等人^[16]对青海门源6.4级地震微震检测得到的40%相似。区域内的寻乌震群、平远震群、安远震群分别检测到遗漏地震36个、15个、8个，为台网目录个数的¹/₂、¹/₃和¹/₅。遗漏地震比例不高的原因一是区域监测能力不高，太小的地震无法记录到信噪比高的波形，二是研究区域虽然小震频繁，但强度不高，可用作模板的地震数量有限。从3个震群的遗漏地震比例可以看出，地震数量越多且震级越大的震群检测出来的遗漏地震比例越高。

表1 寻乌震群增加遗漏地震前后的前兆参数指标Table 1 Precursor parameter index of Xunwu earthquake swarm before and after adding missing earthquakes

采用模板匹配识别方法得到的地震事件震级与台网目录的震级相差不大，86％的震级差在0.3内，表明遗漏地震的震级估计较为可靠。震级—频次统计分析表明加入遗漏事件后，M_L1.5以下地震的目录完整性有较明显的改善。增加遗漏地震后研究区域的b值从0.885变为0.917，该参数的变化对区域地震趋势分析有着重要作用。

对寻乌震群增加遗漏地震前后的前兆震群判定指标进行对比分析，发现U、F、ρ、b值都发生了变化，其中ρ值判定结果由非前兆震群变为前兆震群，说明遗漏地震的补充对前兆震群的判定结果产生了一定影响，对区域趋势判定具有重要意义。

本文对寻乌及邻区遗漏地震检测得到了大量的遗漏地震，丰富了区域地震目录资料，但是补充遗漏地震后M_L0.8以下的地震数量仍与震级—频度拟合直线相差很远，说明检测出来的遗漏地震数量有限，由于研究区域监测能力不高，太小的地震波形信噪比较低，模板匹配滤波技术无法识别出太小的地震。寻乌震群遗漏地震比例比整个研究区域的比例高，主要原因是震群地震之间的互相关程度更高。总的来说，检测出的遗漏地震数量与监测能力、地震波形互相关程度和模板地震数量关系密切。

从本研究结果来看，检测到的遗漏地震有三类：一是震级较小，只有单台记录波形清晰；二是地震信号被淹没在噪声中，无法直接观测；三是受前一个地震尾波的影响，波形相互重叠，难以区分。本研究占比最高的遗漏地震是第一类。总的来说，由于本研究区域监测能力不高，利用遗漏地震事件波形进行后续研究的作用有限，难以对遗漏地震做出精确定位，常用的做法是将遗漏地震事件震中置于互相关系数最大的模板地震震中位置^[7，20]，这样在空间分布上起不到补充作用，对发震构造的研究意义不大。本研究对遗漏地震检测的主要作用是丰富地震目录资料，使地震活动性参数计算更加准确，有利于地震趋势研判。

致谢：本文所用程序由天津市地震局谭毅培高级工程师提供，河北省地震局单连君高级工程师、安徽省地震局倪红玉高级工程师为本研究提供宝贵意见，在此表示感谢。