基于改进支持向量机的微震初至波到时自动拾取方法

基于改进支持向量机的微震初至波到时自动拾取方法
李铁牛1， 胡宾鑫1， 李化坤2， 耿文成2， 郝鹏程3， 纪旭波4， 孙增荣3， 朱峰1， 张华1， 阳铖权1（1. 齐鲁工业大学（山东省科学院） 激光研究所，山东 济南　250300；2. 赤峰吉隆矿业有限责任公司,
内蒙古 赤峰　024328；3. 山东盛隆安全技术有限公司，山东 济南　250032；
4. 山东恒邦矿业发展有限公司，山东 烟台　264199）
摘要：微震初至波到时拾取是实现微震震源高精度定位的重要前提。

传统的人工拾取方法效率低，而自动拾取方法在低信噪比条件下难以准确拾取初至波到时。

针对上述问题，提出了一种基于改进支持向量机（SVM ）的微震初至波到时自动拾取方法。

首先，对原始微震数据进行归一化处理、线性校正和适当裁剪，将微震数据的振幅、能量和相邻时刻的能量比作为特征对数据标记不同类别；然后采用粒子群优化（PSO ）算法和网格搜索法优化SVM 的惩罚参数和核函数参数，即先利用PSO 算法对参数进行大范围的快速定位，得到初步最优解，再以该解为初始位置重新构建参数搜索区间，设置小步长的网格搜索法对参数进行精细搜寻，得到最优参数，并将该最优参数代入SVM 模型进行训练，得到改进SVM 模型；最后根据改进的SVM 模型对微震数据进行分类识别，定义微震波第1个采样点对应的时刻为初至波到时。

采用某矿井下微震监测数据进行实验，结果表明：该方法对微震初至波到时的拾取准确率达96.5%，平均拾取误差为3.8 ms ，在低信噪比情况下仍可对微震初至波到时进行准确拾取，拾取精度高于自动拾取方法中常用的长短时窗能量比（STA/LTA ）法。

关键词：微震监测；微震初至波；初至波到时；支持向量机；网格搜索法；粒子群优化算法中图分类号：TD324 文献标志码：A
Automatic picking method of microseismic first arrival time based on improved support vector machine
LI Tieniu 1, HU Binxin 1, LI Huakun 2, GENG Wencheng 2, HAO Pengcheng 3, JI Xubo 4, SUN Zengrong 3, ZHU Feng 1, ZHANG Hua 1, YANG Chengquan 1
(1. Laser Research Institute, Qilu University of Technology(Shandong Academy of Sciences), Jinan 250300, China ；2. Chifeng Jilong Mining Co., Ltd., Chifeng 024328, China ； 3. Shandong Shenglong Safety Technology Co., Ltd.,
Jinan 250032, China ； 4. Shandong Hengbang Mining Development Co., Ltd., Yantai 264199, China)
Abstract : The microseismic first arrival time picking is an important prerequisite for the high-precision positioning of the microseismic source. The traditional manual picking method is inefficient. The automatic picking method is difficult to pick the arrival time of the first wave accurately under the condition of low signal-to-noise ratio. In order to solve the above problems, an automatic picking method of microseismic first arrival time based on improved support vector machine (SVM) is proposed. Firstly, the method carries out normalization processing, linear correction and proper clipping on original microseismic data. The method marks different categories of the data by taking the amplitude, the energy and the energy ratio of adjacent moments of the microseismic data as features. Secondly, the method adopts a particle swarm optimization (PSO) algorithm and a
收稿日期：2022-05-30；修回日期：2023-03-09；责任编辑：盛男。

基金项目：济宁市重点研发计划项目（2021AQGX001）；山东省自然科学基金博士基金项目（ZR2019BEE019）；山东省自然科学基金
重点项目（ZR2020KC012）；济南市“高校20条”资助项目（2020GXRC032）；济南市“新高校20条”资助项目（2021GXRC037）。

作者简介：李铁牛（1997—），男，山东济南人，硕士研究生，研究方向为信息处理技术，E -mail :Tieniul@ 。

通信作者：胡宾鑫
（1979—），男，浙江景宁人，副研究员，博士，研究方向为传感器技术与智能系统，E -mail :bxhu@qlu .edu .cn 。

引用格式：李铁牛，胡宾鑫，李化坤，等. 基于改进支持向量机的微震初至波到时自动拾取方法[J ]. 工矿自动化，2023，49（3）：63-69.
LI Tieniu, HU Binxin, LI Huakun, et al . Automatic picking method of microseismic first arrival time based on improved support vector machine [J ]. Journal of Mine Automation ，2023，49（3）：63-69
.
第 49 卷 第 3 期工　矿　自　动　化
Vol .49 No .32023 年 3 月
Journal of Mine Automation
Mar . 2023
文章编号：1671−251X （2023）03−0063−07
DOI ：10.13272/j.issn.1671-251x.2022050081
grid search method to optimize the penalty parameters and the kernel function parameters of the SVM. The method carries out large-range fast positioning on the parameters by using the PSO algorithm to obtain a preliminary optimal solution. Then the method re-constructs a parameter search interval by taking the solution as an initial position, sets a small-step grid search method to carry out fine searching on the parameters to obtain the optimal parameters. The method substitutes the optimal parameters into the SVM model to train, and obtains the improved SVM model. Finally, the microseismic data are classified and identified according to the improved SVM model. The time corresponding to the first sampling point of the microseismic wave is defined as the arrival time of the first wave. The microseismic monitoring data from a mine shaft is used for the experiment. The results show that the accuracy of the method for picking the microseismic first arrival time is 96.5%, and the average picking error is 3.8 ms. Under the condition of low signal-to-noise ratio, the microseismic first arrival time can still be picked accurately. The picking precision is higher than the short term average/long term average (STA/LTA) method commonly used in automatic picking methods.
Key words: microseismic monitoring; first arrival of microseism; the arrival time of the first wave; support vector machine; grid search method; particle swarm optimization algorithm
0　引言
微震监测是目前公认的用于矿山动力灾害、岩体破裂风险监测预警的有效手段[1-2]。

通过在监测区域布置多个微震检波器实时采集微震信号，在三维空间中确定微震事件的时空信息及能量，可实现对各区域岩体变形破坏活动范围、稳定性及发展趋势的定性、定量评价[3-5]。

震源定位是微震监测技术的核心，而准确高效的微震初至波（P波）到时拾取是实现震源高精度定位的重要前提[6-8]。

长期以来，微震初至波到时主要依靠地震专家人工拾取，处理效率低，面对海量数据时不易实现[9]。

近年来，随着计算机技术不断发展，微震初至波到时自动拾取方法成为研究热点[10-11]，主要有长短时窗能量比（Short Term Average/Long Term Average，STA/LTA）法[12]、BP神经网络法[13-14]、赤池信息准则（Akaike Information Criterion, AIC）法[15-16]、随机森林法[17]等。

STA/LTA法基于时间域信号能量变化设定阈值函数，其计算简单、运算速度快，但拾取精度易受人为经验所设时窗长度和阈值大小的影响，尤其对低信噪比信号的拾取准确率较低；BP神经网络法训练多层网络模型对微震信号进行分类拾取，但其学习训练速度慢且易陷入局部极小值收敛，导致拾取效率不高；AIC法的基本原理是通过求取微震信号AIC函数局部最小值确定到时，对于高信噪比信号可以实现精确拾取，但在低信噪比情况下容易对函数最小值误判，从而造成初至波的错误拾取；随机森林法通过训练大量决策树构建随机森林模型，通过模型对微震初至波到时进行精准预测，该方法运算速度快、拾取精度高，但受决策树数量的影响，在处理某些噪声较大的信号时，易出现过度拟合的情况。

支持向量机（Support Vector Machine，SVM）具有鲁棒性强、设置参数少、学习能力强等优点[18-19]，可以对小样本数据进行预分析，因此可通过建立SVM模型对微震数据进行初至波到时拾取。

但在实际应用中，SVM模型性能受参数选取的影响很大，如何获得SVM最优参数并提高参数优化速度成为关键。

因此，本文提出了一种基于改进SVM的微震初至波到时自动拾取方法。

首先对微震数据进行预处理；然后采用粒子群优化（Particle Swarm Optimization，PSO）算法和网格搜索法对SVM参数进行不同步长的优化训练，得到改进的SVM模型；最后利用改进SVM模型对微震信号的初至波到时进行拾取，可提高初至波拾取精度和速度。

1　方法原理
1.1　SVM
SVM是一种二分类模型[20]，其主要思想是在数据空间中找出一个最优分类超平面，使得两类数据产生最大分类间隔以完成分类，并借助核函数避免维数灾难问题。

x i y i
∈
y i∈{−1,1}
对于线性可分问题，给定训练样本（，），其中，x i为第i（i=1，2，…，m，m为样本数量）个样本，x i R n（n为空间维数），y i为第i个样本的正负类别，。

SVM在n维空间中构造一个目标函数，寻找一个符合分类目标的超平面：
ω
式中：为超平面的法向量；x={x1，x2，…，x m}；b为分类阈值。

• 64 •工矿自动化第 49 卷
要将所有样本正确区分，需满足以下约束条件：
同时使得最优分类超平面为最大分类间隔的超平面，如图1所示。

图 1 最优分类超平面
Fig. 1 Optimal classification hyperplane
对于线性不可分问题，可转换为极小值求解问题。

ϕξξi ≥0ξ={ξ1,ξ2,···,ξm }式中：（·）为控制分类间隔的目标函数；i 为第i 个样本的松弛变量，；；C 为惩罚参数，其作用是控制目标函数中两项（“寻找空间间隔最大的超平面”和“保证分类数据点偏差量最小”）之间的权重。

K (x i ,x )如果某个分类问题在初始定义的空间中是线性不可分的，可以通过引入核函数将该问题转换到一个更高维的空间中，此时对应的判别函数为
a i 式中为二次规划优化函数求解的Lagrange 因子。

径向基核函数（Radial Basis Function ，RBF ）具有分类准确度高、局部拟合效果好等优势，在大多数条件下具有普适性，因此本文采用该函数。

式中g 为核函数参数，该参数直接影响SVM 模型性能。

SVM 分为3层，如图2所示。

第1层为输入向量x ，第2层为基于m 个支持向量的非线性变换（内积），第3层为输出层。

对惩罚参数C 和核函数参数g 的优化是提高SVM 模型性能的关键。

1.2　PSO 算法
PSO 算法是一种基于群体智能的全局随机搜索
算法[21]，其主要思想是通过鸟类群体中个体之间的
信息共享和相互协作来寻找最优解。

设粒子数量为Q ，迭代次数为N max ，粒子当前第N 次迭代后对应的移动速度和位置分别为V N
和X N 。

每个粒子的更新公式为
式中：V N +1为更新后粒子的移动速度；X N +1为更新后粒子的当前位置；λ为惯性权值；c 1，c 2为学习因子；r 1，r 2为[0，1]区间内不同的随机数；P a ，P b 分别为粒子的个体最优解和种群最优解。

1.3　网格搜索法
网格搜索法可以演算出超参值的最优组合，是
使用最广泛的SVM 参数寻优算法之一[22]。

其基本原理：在指定的参数取值范围内，使待搜索参数按照规定的步长依次调整参数的取值，循环遍历所有可能的取值点进行排列组合形成“网格”，同时将每次的参数组合代入SVM 模型中进行交叉验证，最终将使SVM 模型分类效果最好的参数组合作为全局最优解。

1.4　微震初至波到时拾取流程
在实际应用中，PSO 算法虽然可以快速收敛，但
容易陷入局部最优问题；网格搜索法通过穷举遍历可以准确找到全局最优解，但算法初期耗时较长。

因此，本文将PSO 算法与网格搜索法相结合对SVM 进行改进，以更快速地获得最优参数，兼顾收敛性和运行效率，提高SVM 分类性能。

基于改进SVM 的微震初至波到时自动拾取流程如图3所示，具体步骤如下。

Step1：数据预处理。

先对微震数据进行归一化处理，使数据特征具有同等的度量尺度，提高训练模型的收敛速度。

再对数据进行线性校正来降低误差，保证数据的真实性。

因微震检波器采集的数据窗口过长，不利于SVM 模型训练，需对数据进行适
1
2
m
图 2 SVM 结构
Fig. 2 Support vector machine structure
2023 年第 3 期李铁牛等： 基于改进支持向量机的微震初至波到时自动拾取方法
• 65 •
当剪裁。

当微震初至波来临时，微震信号波形的振幅和能量显著增大，相邻时刻能量比（后一时刻能量与前一时刻能量的比值）也相应增大。

因此提取微震信号波形的振幅、能量和相邻时刻能量比作为特征，对微震数据进行特征标记并存储样本标签。

数据预处理流程如图4所示。

图 4 数据预处理流程
Fig. 4 Data preprocessing process
Step2：改进SVM模型建立。

将预处理后的微震数据划分为训练数据和测试数据。

在训练数据中随机抽取部分数据作为训练样本，采用PSO算法对（C，g）进行大范围的快速定位，得到（C，g）的初步最优解；之后以该解为初始位置重新构建参数搜索区间，设置小步长的网格搜索法对（C，g）进行精细搜索，获得的最终结果即最优参数。

将该最优参数代入SVM模型进行训练，得到改进SVM模型。

Step3：分类预测。

将微震测试数据代入训练好的改进SVM模型中进行数据分类预测。

数据被分为微震波和无效波2类，将微震波第1个采样点对应的时刻定义为初至波到时。

2　方法验证
2.1　实验参数
采用某矿井下微震监测系统记录的微震数据来验证基于改进SVM的微震初至波到时自动拾取方法的有效性。

将微震初至波来临时刻到波形结束时刻的数据标记为1，其余数据均标记为0，保存微震特征样本，形成数据集。

该数据集含有2 000个样本点，按7∶3的比例随机划分为训练集和测试集。

设置PSO算法的最大迭代次数为100，种群数量为30。

设置（C，g）参数的搜索区间为[2−7，27]，学习因子c1=1.8、c2=1.8，采用5折交叉验证对训练样本进行测试。

基于PSO算法得出初步最优参数（C，g）=（5.632 4，0.872 3）。

以该优化结果为中心起始点，缩小参数坐标系网格的取值范围，进行下一次小步距网格精细搜索。

设置C的取值区间为[2，24]，g的取值区间为[0.001，2]，搜索步距为0.001。

经多次搜索更新，得到全局最优参数，如图5所示。

可看出当C=5.018 7，g=0.860 8时，损失函数值最小，准确率达97.6%。

0.050
观测到的点
Best C=5.018 7 g=0.860 8 CVAccuracy=97.6%
模型均值
下一个点
模型最小可行点
0.055
0.045
0.040
0.035
损
失
函
数
值
0.030
0.025
0.9
1.0
0.80.7
g0.60.56C
7
8
5
图 5 网格搜索法参数寻优结果
Fig. 5 Parameter optimization results of grid search algorithm 2.2　实验设置
在Windows 10系统、基于Pycharm的Python3.8集成开发环境下，结合Scikit learn机器学习开发包进行实验。

用微震初至波到时拾取方法中最常用的STA/LTA法与本文方法进行对比。

设置STA/LTA 法的长时窗为180 ms，短时窗为30 ms，触发阈值为3.8。

设置改进SVM的惩罚参数C=5.018 7，核函数参数g=0.860 8。

微震信号在岩体中平均波速一般为3 500 m/s，若要求震源自动定位精度达20 m，则走时误差为5.7 ms。

以人工拾取结果为标准，适当放宽定位精度要求，视拾取误差在8 ms以内为正确拾取。

引入拾取准确率（本文方法正确拾取初至波的数量与人工拾取初至波的数量之比）指标来评估本文方法对初至波到时的拾取效果。

图 3 基于改进SVM的微震初至波到时自动拾取流程
Fig. 3 Automatic picking process of microseismic first-arrival time
based on improved support vector machine
• 66 •工矿自动化第 49 卷
2.3　实验结果与分析
为验证本文方法在低信噪比条件下的拾取效
果，对原始微震信号进行加噪处理，得到信噪比分别
为6和1.8的微震数据。

人工拾取方法、STA/LTA 法和本文方法的微震初至波到时拾取结果如图6所示。

1.234
−202 1.191
1.195
1.264
1.215
1.259
1.219
1.274
1.2271.262
1.231
1.189
1.264
−1.5
−0.50.51.5 1.213
−0.15−0.050.050.15 1.251
−0.20.2 1.235−2−1012 1.193
−0.2
00.2 1.268
−1.5
−0.50.51.5 1.218−0.15−0.050.050.15 1.246
−0.15−0.050.050.15加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)加速度/(m ·s −2)1.253
−1.5
−0.50.51.5 1.211−0.2
00.2 1.259
−2−1012 1.185
−0.6
−0.20.20.6 1.229时间/s （a ） 人工拾取方法
（b ） STA/LTA 法 (信噪比为 6)（d ） STA/LTA 法 (信噪比为 1.8)（c ） 本文方法 (信噪比为 6)（e ） 本文方法 (信噪比为 1.8)
0.4−0.6
0.60.44−2024时间/s −0.20.2−0.6
0.6−0.20.2−0.6
0.6−0.20.2−0.6
0.6−2−1012−0.2
00.20.4−1.5
−0.50.51.5−0.15−0.050.050.15−0.2
00.20.4−1.5
−0.50.51.5−0.15−0.050.050.15−0.2
00.20.4−0.4
−0.200.2
0.4 1.2加速度/(m ·s −2) 1.01.21.41.61.82.0
时间/s 1.198加速度/(m ·s −2)−0.4
−0.200.20.41.01.21.41.61.82.0
1.204加速度/(m ·s −2)−0.4
−0.200.20.41.01.21.41.61.82.0
时间/s 1.189加速度/(m ·s −2)−0.4
−0.200.20.41.01.21.41.61.82.0
时间/s
加速度/(m ·s −2) 1.01.21.41.61.82.0
−0.4
−0.200.2 1.1950.4图 6 微震信号加噪后初至波到时拾取结果
Fig. 6 First-arrival time picking results after microseismic signal plus noise
由图6可知：当信噪比为6时，本文方法和STA/LTA 法的平均拾取误差分别为3 ms 和7 ms ，均能做到有效拾取；当信噪比为1.8时，本文方法和STA/LTA 法的平均拾取误差分别为5.3 ms 和11.5 ms ，本文方法仍能做到有效拾取，而STA/LTA 法的拾取效果不太理想。

可以看出，在低信噪比条件下，本文方法的拾取误差虽然有一定增加，但仍可对微震初至波进行准确拾取。

为进一步验证本文方法的有效性，选取某矿井下12个微震检波器2022年1−6月采集的4 320组微震数据作为测试集，STA/LTA 法和本文方法的实验结果见表1。

从表1可看出，本文方法对微震初至波到时的拾取准确率高于STA/LTA 法，平均拾取准确率达96.5%，平均拾取误差仅为3.8 ms 。

随机抽取同一微震检波器（共8个通道）在不同时刻采集的微震数据，人工拾取方法、STA/LTA 法和本文方法的微震初至波到时拾取结果如图7所示。

由图7可知：当微震信号信噪比较高时，本文方法和STA/LTA 法对初至波到时的拾取误差均较低；
当微震信号信噪比较低时，本文方法的初至波到时拾取误差明显低于STA/LTA 法；本文方法和STA/
表 1 微震初至波到时统计结果
Table 1 Statistical results of microseismic first-arrival time
微震检波器序号微震信号数量/组人工
拾取数量/组自动拾取数量/组拾取准确率/%拾取误差/ms STA/LTA 法本文方法STA/LTA 法本文方法STA/LTA 法本文方法136035732335190.598.3 6.5 3.5236035531534388.796.67.3 4.2336035631834689.397.2 5.9 2.9436035230733587.295.27.4 4.5536035431633989.295.8 6.7 3.4636035231133488.394.97.5 5.1736035733134892.797.5 6.9 3.7836035531234187.996.17.1 4.3936035832935391.998.6 6.3 2.71036035431633689.394.9 6.9 3.91136035330933587.594.97.7 4.21236035632234890.497.8 6.6 3.6平均值
360
354.9
317.4
342.4
89.4
96.5
6.9
3.8
2023 年第 3 期
李铁牛等： 基于改进支持向量机的微震初至波到时自动拾取方法• 67 •
LTA 法对微震初至波到时自动拾取的平均误差分别为4 ms 和9.3 ms ，本文方法对微震初至波到时的拾取精度高于STA/LTA 法。

3　结语
结合PSO 算法和网格搜索法的优点对SVM 进
行改进，优化SVM 的惩罚因子和核参数，提高了SVM 模型性能，并将改进SVM 应用于微震初至波到时自动拾取。

与STA/LTA 法相比，该方法无需设置时窗长度和触发阈值，消除了不同时窗和阈值对拾取精度的影响，微震初至波到时拾取精度明显提升；在低信噪比条件下，该方法仍可对微震信号进行精确的初至波到时拾取。

参考文献（References ）：
孙继平，程继杰. 煤矿冲击地压和煤与瓦斯突出感知
报警方法研究[J ]. 工矿自动化，2022，48（1）：1-6.SUN Jiping ，CHENG Jijie. Study on the perception and alarm method of coal mine rock burst and coal and gas outburst [J ]. Industry and Mine Automation ，2022，48（1）：1-6.[ 1 ]
马天辉，唐春安，唐烈先，等. 基于微震监测技术的岩
爆预测机制研究[J ]. 岩石力学与工程学报，2016，35（3）：470-483.
MA Tianhui ，TANG Chun'an ，TANG Liexian ，et al.Mechanism of rock burst forecasting based on micro-seismic monitoring technology [J ]. Chinese Journal of Rock Mechanics and Engineering ，2016，35（3）：470-483.
[ 2 ]
−0.2−0.400.20.4 1.191
−0.2−0.100.10.2 1.204−0.1−0.200.10.2−0.1−0.200.10.2 1.2031.106
1.155−0.10−0.0500.050.10 1.202 51.178（a ） 人工拾取方法加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)
加速度/(m ·s −2)−0.1−0.200.10.2加速度/(m ·s −2)
−0.1−0.2
00.10.2加速度/(m ·s −2)
−0.2−0.400.20.4 1.184−0.2−0.100.10.2 1.198−0.1−0.200.10.2−0.1−0.200.10.2 1.2041.11
1.148−0.10−0.0500.050.10 1.207 51.176（c ） 本文方法加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)
加速度/(m ·s −2)−0.1−0.200.10.2加速度/(m ·s −2)−0.1−0.2
00.10.2加速度/(m ·s −2)−0.2−0.400.20.4 1.211
−0.2−0.100.10.2 1.213−0.1−0.200.10.2−0.1−0.200.10.2 1.21.118
1.163−0.10−0.0500.050.10 1.217 51.185（b ） STA/LTA 法加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)
加速度/(m ·s −2)
−0.1−0.200.10.2加速度/(m ·s −2)
−0.1−0.2
00.10.2加速度/(m ·s −2)
1.01.11.21.31.41.51.61.71.81.9
2.0
−0.4−0.6
−0.200.20.40.6时间/s
1.247加速度/(m ·s −2)
1.01.11.21.31.41.51.61.71.81.9
2.0
−0.4−0.6
−0.200.20.40.6时间/s
1.247加速度/(m ·s −2)1.01.11.21.31.41.51.61.71.81.9
2.0
−0.4−0.6
−0.200.20.40.6时间/s
1.247 5加速度/(m ·s −2)
加速度/(m ·s −2)
加速度/(m ·s −2)加速度/(m ·s −2)图 7 8通道微震初至波到时拾取结果
Fig. 7 8-channel microseismic first-arrival time picking results
• 68 •
工矿自动化
第 49 卷
刘强，蔡永顺，石峰，等. 微震监测技术在某金矿的应
用与研究[J ]. 中国矿业，2020，29（增刊1）：145-148.LIU Qiang ，CAI Yongshun ，SHI Feng ，et al.Application and study of microseismic monitoring technology in a gold mine [J ]. China Mining Magazine ，2020，29（S1）：145-148.[ 3 ]
WAMRIEW D ，DORHJIE D B ，BOGOEDOV D ，et al.
Microseismic monitoring and analysis using cutting-edge technology ：a
key
enabler
for
reservoir
characterization [J ]. Remote Sensing ，2022，14（14）：3417-3438.[ 4 ]
LIU Fei ，MA Tianhui ，TANG Chun'an ，et al. Prediction
of rockburst in tunnels at the Jinping II hydropower station using microseismic monitoring technique [J ].Tunnelling and Underground Space Technology ，2018，81：480-493.[ 5 ]
李楠，王恩元，GE Maochen. 微震监测技术及其在煤
矿的应用现状与展望[J ]. 煤炭学报，2017，42（增刊1）：83-96.
LI Nan ，WANG Enyuan ，GE Maochen. Microseismic monitoring technique and its applications at coal mines ：present status and future prospects [J ]. Journal of China Coal Society ，2017，42（S1）：83-96.[ 6 ]
GESRET A ，DESASSIS N ，NOBLE M ，et al.
Propagation of the velocity model uncertainties to the seismic event
location [J ].
Geophysical
Journal
International ，2014（1）：52-66.
[ 7 ]
宋广东. 基于光纤传感的微震信号采集和识别与震源定位研究[D ]. 北京：中国矿业大学（北京），2019.SONG Guangdong. Research on acquisition and recognition of microseismic signal and source location based on optical fiber sensing [D ]. Beijing ：China University of Mining and Technology-Beijing ，2019.[ 8 ]
STEFANO R D ，ALDERSONS F ，KISSLING E ，et al.
Automatic seismic phase picking and consistent observation error assessment ：application to the Italian seismicity [J ]. Geophysical Journal International ，2006，165（1）：121-134.[ 9 ]
贾瑞生，谭云亮，孙红梅，等. 低信噪比微震P 波震相初
至自动拾取方法[J ]. 煤炭学报，2015，40（8）：1845-1852.
JIA Ruisheng ，TAN Yunliang ，SUN Hongmei ，et al.Method of automatic detection on micro-seismic P-arrival time under low signal-to-noise ratio [J ]. Journal of China Coal Society ，2015，40（8）：1845-1852.[10]
高新成，杜功鑫，王莉利，等. 深度学习在地震初至拾
取中的应用综述[J ]. 计算机技术与发展，2022，32（8）：1-6.
GAO Xincheng ，DU Gongxin ，WANG Lili ，et al.Application of deep learning in earthquake first break picking [J ]. Computer Technology and Development ，2022，32（8）：1-6.
[11]
LEE M ，BYUN J ，KIM D ，et al. Improved modified
energy ratio method using a multi-window approach for accurate arrival picking [J ]. Journal of Applied Geophysics ，2017，139：117-130.[12]
MARTON S ，LUDTKE S ，BARTELT C. Explanations
for neural networks by neural networks [J ]. Applied Sciences ，2022，12（3）：980-993.[13]
刘佳楠，武杰. 基于全卷积神经网络的地震初至波拾
取[J ]. 信息技术与网络安全，2018，37（11）：58-63.LIU Jianan ，WU Jie. Seismic first arrival picking based on
full
convolutional
networks [J ].
Information
Technology and Network Security ，2018，37（11）：58-63.[14]
李洪丽，张晗. 基于长短时均值比法、偏振法和AIC 法
三种微地震初至拾取方法的对比[J ]. 世界地质，2020，39（3）：649-655，663.
LI Hongli ，ZHANG Han. Comparison of three microseismic first arrival picking methods based on ratio of short-term average and long-term average ，polarization method and AIC method [J ]. Global Geology ，2020，39（3）：649-655，663.[15]
王洪超. 基于Fast−AIC 算法的微地震事件初至拾取及
自动识别技术研究[D ]. 长春：吉林大学，2017.WANG Hongchao. Research on automatic microseismic event detection and arrival picking based on Fast-AIC algorithm [D ]. Changchun ：Jilin University ，2017.[16]
高煜，胡宾鑫，朱峰，等. 微震初至波到时自动拾取研
究[J ]. 工矿自动化，2020，46（12）：106-110.
GAO Yu ，HU Binxin ，ZHU Feng ，et al. Research on automatic picking of microseismic first arrival [J ].Industry and Mine Automation ，2020，46（12）：106-110.[17]
李亦滔. 基于支持向量机的改进分类算法[J ]. 计算机
系统应用，2019，28（10）：145-151.
LI Yitao. Improved classification algorithm based on support vector machine [J ]. Computer System &Applications ，2019，28（10）：145-151.[18]
丁世飞，齐丙娟，谭红艳. 支持向量机理论与算法研究
综述[J ]. 电子科技大学学报，2011，40（1）：2-10.DING Shifei ，QI Bingjuan ，TAN Hongyan. An overview on theory and algorithm of support vector machines [J ].Journal of University of Electronic Science and Technology of China ，2011，40（1）：2-10.[19]
MA Qian ，SUN Cong ，CUI Baojiang ，et al. A novel
model for anomaly detection in network traffic based on kernel support vector machine [J ]. Computers &Security ，2021，104. DOI ：10.1016/J.COSE.2021.102215.[20]
LYU Zhiming ，ZHAO Jun ，WANG Wei ，et al. A
multiple surrogates based PSO algorithm [J ]. Artificial Intelligence Review ，2019，52（4）：2169-2190.[21]
SUN Yuting ，DING Shifei ，ZHANG Zichen ，et al. An
improved grid search algorithm to optimize SVR for prediction [J ]. Soft Computing ，2021（10）：5633-5644.
[22]
2023 年第 3 期
李铁牛等： 基于改进支持向量机的微震初至波到时自动拾取方法
• 69 •。