引言

地震目录是区域地震活动性研究、地震地质和地震危险性分析的基础资料,其数据完整性、准确性决定了地震工作者对区域地震活动的认识,在一定程度上决定了区域地震活动性研究和地震危险性分析的可靠程度(谭毅培等,2104a)。对于一次地震事件,完整地统计其余震序列数目、震级、发震时刻、震中地点,能使地震工作者更准确地分析地震活动规律,从而作出更准确全面的结论(谭毅培等,2104b)。目前,地震目录的统计主要靠人工分析,而大震后的波形干扰往往较大,且在余震较多的情况下易叠加,人工辨别受到很大影响。新疆台站密度相对稀疏,在距离震源远、震级微弱的情况下,台站记录的震相难以辨别。基于此,中国地震局监测预报司预报管理处组建的地震预报青年组地震序列小组研发了SEPD(Seismic Events and Phase Detection)软件(《SEPD用户手册-地震事件与震相检测软件v1.0》),可对地震序列中地震时空演化特征、发震构造几何形态的精细确定等科学问题,提供更完整可靠的地震序列目录和震相报告。谭毅培等(2014a, 2014b)通过该软件对2013年3月涿鹿微震群及2013年8月河北蔚县小震群地震序列目录进行分析,尹欣欣等(2018)通过该软件对青海门源Ms6.4地震进行分析研究,冯雪玲等(2019)用Morlet小波对微地震进行检测。

1 互相关扫描

2018年11月25日20时23分新疆博乐发生MS4.9地震。该地区位于新疆天山西段北麓,艾比湖西,准噶尔盆地西南部,地处阿拉套山和岗吉格山间谷地,地势西部高、东部低,高山、中山、低山丘陵和谷地平原呈阶梯状分布。其周边台站相对稀疏且分布不均匀(李艳永等, 2018, 2019),距离该地震最近的博乐台(BUL)超过50km(图 1),因受背景噪声、台基等因素的影响,ML1.0以下地震记录不清楚,造成人工识别震相困难。本文通过SEPD软件对该地震2018年11月25日-2018年11月27日连续波形进行互相关扫描和遗漏地震检测,从而完善余震地震目录,以使对该地震事件余震序列有更准确的认识。


图 1 震中位置示意图 Fig. 1 Station distribution of Bole earthquake

本文所用数据源自新疆台网监测中心记录的三分量数字化地震波形,并在Jopens Data Service网站下载震相报告数据进行检测。检测前需选择合适的模板地震,总体原则为波形记录清晰、信噪比高。为此,本文选取ML2.0以上的地震事件作为模板地震(表 1),为找到清晰的波形,每个地震事件均选取震中距最近的博乐台(BUL)、温泉台(WNQ)、阿拉山口台(ALS)3个台站记录的连续波形作为研究对象。

表 1 选取的模板地震 Table 1 Selected template events

互相关扫描参数的设置选用4阶Butterworth零相移数字滤波器1-4Hz滤波,截取地震观测报告信息中S波震相到达前2.0s至震相到达后2.0s的波形。参数设置完毕后,软件以选取的3个地震为模板,分别对每个模板3个台站进行扫描,并计算互相关系数(Cross- correlation,CC)。由于互相关扫描是通过匹配模板进行的,不能确保每个余震事件均与模板相似,所以可能存在漏检的情况,可通过增加模板数量减少漏检事件。

2 余震检测

连续波形和地震模板经过重采样,采样间隔由0.01s变为0.05s。本文通过取CC结果中三分向的平均值,并计算其绝对离差中位数MADMAD=median|Xi|-X),得到CC三分向平均值绝对离差中位数为0.09939。Peng等(2009)认为,如果CC之和大于9倍MAD,认为每天误检测的地震事件少于1次,可作为地震事件检测阈值。本文阈值取9倍MAD,即0.89。模板地震Eq1125202334阿拉山口台2018年11月25日20:20-21:10互相关扫描结果及阈值关系如图 2所示,图中有17个疑似地震的发震时刻。依据程序中的设定,所选3个台站CC之和大于0.89即认为检测到地震事件,可判定此时刻存在地震事件,所检测出的地震目录中没有的事件在本文中称为遗漏地震事件。为减少误差,在震相校正模块对疑似地震发震时刻进行人工标定目录遗漏地震的震相到时。


图 2 互相关扫描结果示意图 Fig. 2 Diagrammatic map of cross-correlation scanning result

搜索遗漏地震事件发震时刻P波、S波到时方法示意如图 3所示,图中灰色为连续波形,黑色为模板地震Eq1125202334,Cmax为最大互相关系数。最终检测出47个地震事件,较地震目录中15个地震事件多32个。模板地震Eq1125202334、Eq1125202932、Eq1125203419检测出的余震事件数量分别是16个、20个、28个,且每个模板地震均有独自检测出的余震事件,说明所选模板地震震级不能决定检测出余震事件的数量,而取决于余震事件与模板的互相关性,选取模板时,在保证波形清晰、信噪比较高的情况下,尽可能多地选取模板地震,以包含尽可能多的余震事件。本文选择的3个模板均为此次地震事件中震级较大的事件,且选择了震中距最近的台站波形,保证了较高的信噪比和清晰波形,最大程度保证了模板地震的必要条件。


图 3 地震发震时刻P波、S波到时方法示意图 Fig. 3 The process of detecting origin time, P-wave and S-wave arrival time of missing earthquakes

部分检测结果见表 2,其中三台互相关系数均值为1,代表模板地震检测到自身,这是判断互相关扫描结果是否准确的标准之一。结果中三台互相关系数均值明显偏小,多为0.4-0.6,这是因为检测出的地震事件震级较小,只有距离更近的2个台站有较好的互相关系数,而距离较远的ALS台站互相关系数较差,因此降低了3个台站的互相关系数均值,其余2个台站所测结果互相关系较理想。同一地震事件会被多个模板地震检测到,将不同模板地震检测到的同一时刻地震事件划分为一组,为一次地震事件。

表 2 通过模板所检测出的部分地震事件 Table 2 Partial seismic events detected by template events
3 震级的确定

震级测定使用SEPD软件震相校正模块实现。通过震相校正,检测出遗漏地震的发震时刻和震级见表 3。通过遗漏地震事件和目录地震事件对比,目录原有地震事件15个,共检测出32个遗漏地震事件,由于负震级震相特征微弱,更不易检测,为保证检测结果的可靠性,本文研究以ML0.0以上为主。

表 3 遗漏地震事件发震时刻与震级 Table 3 The origin time and magnitude of missing earthquakes events
4 地震目录完整性

相比原有地震目录,检测结果多出213%,检测出地震事件的震级为ML-0.3-1.1,其中ML0.0-0.5地震最多。地震频次与震级之间存在的关系是当前地震事件目录完整性分析的重要依据,通过遗漏地震检测,ML0.0-1.0地震序列得到补充,填补了检测前ML0.0-1.0地震缺失的情况。该次地震序列事件中ML2.0-2.9、ML3.0-3.9、ML5.0-5.9地震事件各出现一次,为去除其偶然性,本文以ML0.0-2.9地震为研究对象,在该范围下,检测后的结果使震级与频次之间具有更好的线性关系(图 4)。


图 4 遗漏地震检测前后地震频次对比图 Fig. 4 Area chart before and after missing detection

地震目录最小完整性震级表征台网监测能力水平,能填补地震目录数据空白,使地震工作者根据研究方向选择所需地震事件,以达到对地震活动性研究更详细、可靠的目的(苏有锦等,2003冯建刚等,2012)。由于数据量有限,本研究分别以原有地震目录和检测出的地震事件为基础,采用最大曲率法分析二者最小完整性震级。通过分析研究事件震后3d的检测目录,得出台网目录与检测目录最小完整性震级分别为ML1.6、ML0.8(图 5),检测出的地震目录最小完整性震级比台网目录最小完整性震级明显降低,从而提高了台网监测能力。


图 5 检测前后最小完整性震级对比 Fig. 5 Contrast of minimum integrity magnitude before and after detection
5 结论

本文通过匹配滤波技术对2018年11月25日新疆博乐MS4.9地震及余震序列中遗漏地震进行检测,并对发震时刻和震级进行确定,检测出遗漏地震32条,其中84.4%为ML0.0-1.0,9.4%小于ML0.0,较地震目录中原有15条地震多213%,使ML0.0-1.0地震目录得以完善,完善后的震级和频次呈现出更好的线性关系。最小完整性震级由ML1.6减至ML0.8,从而提高了台网监测能力。通过对连续波形的对比分析,认为本次拾取遗漏地震结果可靠。

相比人工分析地震,该技术能较大程度上排除背景噪声及大震面波干扰。在与主震相邻的时间段内,地震波形产生较大面波,并受地脉动、地震事件叠加及台站距离较远、台基不理想等因素的影响,较小地震震相难以辨别,无法分析,该软件通过连续波形与地震模板波形的对比,更易发现小震并进行分析,以减少人工识别误差,获得更完善的地震目录,使地震工作者对某一地区地震活动性的分析研究更全面准确,并使地震危险性分析更可靠。但该软件也有不足,首先,该软件依赖于模板地震,当模板地震选取不典型或数量不够时,易遗漏所要检测的地震,当余震波形特征与主震或模板不同时,易造成检测结果不准确;其次,当信噪比不高或波形记录信号与噪声频带重叠时,检测效果不理想;最后,在地震校正计算中,需凭借地震工作者的经验完成,且互相关扫描所需时间较长(如检测1d的波形大概需要1d),降低了工作效率。

致谢: 感谢谭毅培副研究员提供的SEPD软件及新疆维吾尔自治区地震局监测中心提供的连续波形资料。
参考文献
冯建刚, 蒋长胜, 韩立波, 等, 2012. 甘肃测震台网监测能力及地震目录完整性分析[J]. 地震学报, 34(5): 646-658. DOI:10.3969/j.issn.0253-3782.2012.05.006
冯雪玲, 刘东亚, 李志海, 2019. 基于小波互相关的微地震信号检测方法研究[J]. 大地测量与地球动力学, 39(9): 982-986.
李艳永, 王成虎, 杨佳佳, 2018. 呼图壁地区震源机制解及构造应力场特征分析[J]. 大地测量与地球动力学, 38(12): 1246-1250.
李艳永, 唐明帅, 乌尼尔, 2019. 新疆北天山中东段呼图壁地区震源深度的重新测定[J]. 震灾防御技术, 14(2): 352-362.
苏有锦, 李永莉, 李忠华, 等, 2003. 川滇地区区域地震目录完整性最小震级分析[J]. 地震研究, 26(S1): 10-16.
谭毅培, 曹井泉, 卞真付, 等, 2014a. 2013年8月河北蔚县小震群遗漏地震检测与发震构造分析[J]. 地震学报, 36(6): 1022-1031.
谭毅培, 曹井泉, 刘文兵, 等, 2014b. 2013年3月涿鹿微震群遗漏地震事件检测和发震构造分析[J]. 地球物理学报, 57(6): 1847-1856.
尹欣欣, 杨立明, 赵林林, 等, 2018. 青海门源Ms6.4地震微震检测与目录完备性研究[J]. 大地测量与地球动力学, 38(7): 760-764, 770.
Peng Z. G., Zhao P., 2009. Migration of early aftershocks following the 2004 Parkfield earthquake[J]. Nature Geoscience, 2(12): 877-881.