基于高分辨率时频分析的核爆地震信号能量统计初至检测算法.pdfVIP

第44卷第1期核电子学与探测技术Vol.44No.1

2024年1月NuclearElectronicsDetectionTechnologyJan.2024

基于高分辨率时频分析的核爆地震信号

能量统计初至检测算法

李鸿儒，李夕海，牛超

（火箭军工程大学，西安710025）

摘要：核爆地震波检测是禁核试核查的关键环节之一。常用的短时长时能量均值比、赤池信息准

则、傅里叶变换等方法虚警率都比较高，为此，提出一种基于高分辨率时频分析的核爆地震信号能量统

计初至检测算法。其基本原理是根据核爆地震事件信号与噪声的频域分布差异，采用窄带分频技术将

信号时域波形变换到时频域，在各频带计算局部能量均值与全局能量均值之比并进行二值化处理，当结

果大于设定的阈值时，认为此频带该时刻属于地震事件信号的一部分，如果在某一阈值下出现二值化结

果相同的情况，只需继续增大阈值进行检测，直到最终只有一个时间点检测到频段的数量最多即可，最

大值的跳跃点即认为是初值时刻。利用地下核爆产生的地震信号进行实验研究，结果表明本文方法能

够在提高检测准确率的同时解决高虚警问题，在低信噪比信号条件下具有较好的适用性与鲁棒性。

关键词：核爆；初至拾取；短时长时能量比；时频分析；分频处理

中图分类号：P31文献标志码：A文章编号：0258-0934(2024)1-0180-11

0引言法可分为在时域、频域以及时频域中进行这三

类,随着深度学习的蓬勃发展，也有越来越多的

地下核试验所激发的震源与天然地震均以学者将地震检测与深度学习结合，取得了不错

波的形式传播，所以对核试验的监测就成为地的效果。

震监测工作的重要组成部分，其中地震监测最其中时域中进行的方法主要有短时长时能

为重要的环节便是地震波检测。地震波检测原量比法（STA/LTA）、自回归赤池信息准则法

理是利用地震信号与背景噪声在幅值、频率、偏（AIC）、分形维数法与相关法等；在频域中处理

振和统计特性等特征上的差异将地震波从背景的方法主要有傅里叶变换法和瞬时频率法；时

噪声中识别出来[]。常用的地震信号检测方频分析方法主要有短时傅里叶变换法（STFT）、

小波变换法、S变换法等[2.3];深度学习方法则

主要利用卷积神经网络以及U形神经网络等

收稿日期：2023-11-17进行地震波初至检测以及震相识别[4.6]

基金项目：陕西省自然科学基金面上项目（2023-

JC-

虽然方法多种多样，但各种方法都有自己

YB-221）；陕西省自然科学基础研究计划项目（2023-

的优势和局限性：STA/LTA方法用信号短时窗

JC-YB-244)。

作者简介：李鸿儒（1999一），男，硕士，主要从事地震信能量平均值和长时窗能量平均值之比随时间的

号检测及分类识别、核爆地震模式识别；E-mail：动态变化来反映信号幅度（或能量）特征的变

163.com。化，并以此识别地震波初至时刻[7]。虽然计算

180

简单，但对于信噪比较低的信号则容易发生漏大小不随频率变化的缺点，对缓慢变化的低频

检，同时对于强干扰噪声的鲁棒性能较差。信号和迅速变化的高频信号具有很好的刻画能

AIC准则首先确定自回归模型（AR）的阶数，然力，能