专利摘要

本发明公开了一种基于SSNet模型的地震事件到时识别方法，包括：(1)收集大量的地震监测台站地面震动波形记录以及相应的地震目录；(2)提出SSNet模型解决地震事件到时识别问题；(3)利用地震事件样本以及随机采集的非地震事件地面震动样本对SSNet模型进行训练，得到用于地震事件检测的识别模型；(4)利用地震事件样本以及地震震相到时记录对SSNet模型进行训练，得到用于地震震相到时识别的识别模型。本发明根据地震波形数据的特点，针对性地设计了深度网络模型，综合利用卷积神经网络等技术抽取地震波形数据的特征，捕捉数据的特点，提高了对地震事件检测任务的准确率和震相到时拾取任务的精度。

权利要求

1.一种基于SSNet模型的地震事件到时识别方法，包括如下步骤：

(1)收集地震监测台站监测到的地面震动波形记录以及相应的地震目录，所述地震目录包括监测到的一次地震事件中P波和S波到达台站的时间；

(2)根据步骤(1)所收集的数据构建训练所需的样本，包括大量对应地震事件的正样本以及大量对应非地震事件的负样本；

(3)利用正样本和负样本训练SSNet模型，得到用于地震事件检测的识别模型，进而向该模型输入地面震动波形数据即可输出判断是否为地震事件；

(4)利用正样本训练SSNet模型，得到用于地震震相到时的预测模型，进而向该模型输入地面震动波形数据即可输出得到P波和S波到达台站的时间。

2.根据权利要求1所述的地震事件到时识别方法，其特征在于：所述正样本包含地面震动波形数据以及对应实际监测到的一次地震事件中P波和S波到达台站的时间，负样本只包含地面震动波形数据；正样本的地面震动波形数据长度为30s，即从P波到达台站之前的5～15s截取起点，从P波到达台站之后的15～25s截取终点；负样本的地面震动波形数据长度也为30s，即从不涉及任何地震事件的连续地面震动波形记录中随机选取确定起点，截取30s长度的地面震动波形数据。

3.根据权利要求1所述的地震事件到时识别方法，其特征在于：所述步骤(2)中保证正样本和负样本的数量相同。

4.根据权利要求1所述的地震事件到时识别方法，其特征在于：所述步骤(3)中对SSNet模型进行训练的具体过程如下：

3.1构建SSNet模型，并初始化该模型参数，包括每一层的权值矩阵、学习率、优化方法以及最大迭代次数；

3.2将所有样本中的地面震动波形数据逐一输入至SSNet模型中进行训练，计算SSNet模型输出的识别结果与对应样本的真值之间的误差函数，其中模型输出的识别结果为样本属于地震事件的概率值，正样本的真值为1，负样本的真值为0；进而通过反向传播对整个SSNet模型的参数不断进行调整更新，直至误差函数收敛或达到最大迭代次数，从而完成训练得到用于地震事件检测的识别模型。

5.根据权利要求1所述的地震事件到时识别方法，其特征在于：所述步骤(4)中对SSNet模型进行训练的具体过程如下：

4.1构建SSNet模型，并初始化该模型参数，包括每一层的权值矩阵、学习率、优化方法以及最大迭代次数；

4.2将正样本中的地面震动波形数据逐一输入至SSNet模型中进行训练，计算SSNet模型输出的预测结果与对应样本的真值之间的误差函数，其中模型输出的预测结果为P波和S波到达台站的时间，样本的真值为实际监测到的一次地震事件中P波和S波到达台站的时间；进而通过反向传播对整个SSNet模型的参数不断进行调整更新，直至误差函数收敛或达到最大迭代次数，从而完成训练得到用于地震震相到时的预测模型。

6.根据权利要求1所述的地震事件到时识别方法，其特征在于：所述SSNet模型从输入到输出由尺度选择卷积块SSBlock1、平均池化层P1、尺度选择卷积块SSBlock2、平均池化层P2、尺度选择卷积块SSBlock3、平均池化层P3、尺度选择卷积块SSBlock4、平均池化层P4以及全连接层依次连接组成。

7.根据权利要求6所述的地震事件到时识别方法，其特征在于：所述尺度选择卷积块SSBlock1～SSBlock4结构相同均由三个尺度选择卷积单元SSCell1～SSCell3级联构成，每个尺度选择卷积单元包含普通卷积层和空洞卷积层，普通卷积层采用大小为3的卷积核对卷积单元的输入以步长为1进行卷积操作并经激活函数Relu处理后输出；空洞卷积层采用大小为3的卷积核对普通卷积层的输出以步长为1进行卷积操作，其输出分成两路分别经激活函数tanh和sigmoid处理后进行点乘，得到的结果与普通卷积层的输出相加并经批归一化处理后作为卷积单元的输出；SSCell1中空洞卷积层的空洞大小为1，SSCell2中空洞卷积层的空洞大小为2，SSCell3中空洞卷积层的空洞大小为4。

8.根据权利要求6所述的地震事件到时识别方法，其特征在于：所述尺度选择卷积块SSBlock1中普通卷积层和空洞卷积层所采用的卷积核数量均为32，尺度选择卷积块SSBlock2中普通卷积层和空洞卷积层所采用的卷积核数量均为64，尺度选择卷积块SSBlock3中普通卷积层和空洞卷积层所采用的卷积核数量均为128，尺度选择卷积块SSBlock4中普通卷积层和空洞卷积层所采用的卷积核数量均为256。

9.根据权利要求6所述的地震事件到时识别方法，其特征在于：所述平均池化层P1～P4结构相同均采用大小为4的池化窗口，步长为4，padding模式为same，即在窗口长度不足4时亦计算得到输出。

10.根据权利要求6所述的地震事件到时识别方法，其特征在于：所述全连接层采用三层神经网络结构，从输入至输出各层的神经元数量分别为256、64、16且均采用激活函数Relu处理。

说明书

技术领域

本发明属于地震预警技术领域，具体涉及一种基于SSNet模型的地震事件到时识别方法。

背景技术

地质工作者多年来利用地震监测台站所记录的连续波形数据在地震预警、地震速报、地震机制等领域进行了大量的工作和研究；其中地震事件检测、震相自动拾取以及震级的快速估算等任务是相关研究的重点和热点。地震事件检测是指对于地震监测台站记录的连续波形数据，检测出其中某一段数据是否属于一个地震事件，从而给出某个时刻是否发生地震的判断；震相的自动拾取是指在我们已经知道某段波形数据属于地震事件之后，给出确切的地震波传播到达台站的时间。

对于地震监测台站波形数据的研究本质上是数据驱动的，研究者可以从大量的监测台站的数据中挖掘、识别地震的特性、模式。传统上地质研究者多从参数拟合的角度出发，以期找到合适的参数来判断一段地面震动数据是否从属于一个地震事件，描述地震的震相、震级等参数与地震波形数据之间的关系。近年来，随着深度学习技术的飞速发展，越来越多的研究者开始在地震波形数据研究的各个任务上应用深度学习。

例如2018年Perol等人在science子刊上发表的文献[Perol T,Gharbi M,DenolleM.Convolutional neural network for earthquake detection and location[J].Science Advances,2018,4(2):e1700578.]提出使用一个8层的卷积神经网络来判断一段台站地面震动数据是否为地震，并且判断地震震中的大致方位。2019年Zhou等人发表的文献[Zhou Y,Yue H,Kong Q,et al.Hybrid Event Detection and Phase-PickingAlgorithm Using Convolutional and Recurrent Neural Networks[J].SeismologicalResearch Letters,2019,90(3):1079-1087.]提出使用循环神经网络等方法来拾取震相的到时等等。这些方法将深度学习技术应用在了地震波形数据上，并取得了超过传统方法的效果，使得这个领域被越来越多人注意到，具有很重要的意义。但是这些现有技术仍然停留在套用基础的神经网络模型来解决地震领域问题，缺少针对于地震数据特点的模型设计、改进与调优。

发明内容

鉴于上述，本发明提供了一种基于SSNet模型的地震事件到时识别方法，其通过大量地震监测台站的地面震动连续波形记录，训练SSNet模型，提高对地震事件的到时识别能力。

一种基于SSNet模型的地震事件到时识别方法，包括如下步骤：

(1)收集地震监测台站监测到的地面震动波形记录以及相应的地震目录，所述地震目录包括监测到的一次地震事件中P波和S波到达台站的时间；

(2)根据步骤(1)所收集的数据构建训练所需的样本，包括大量对应地震事件的正样本以及大量对应非地震事件的负样本；

(3)利用正样本和负样本训练SSNet模型，得到用于地震事件检测的识别模型，进而向该模型输入地面震动波形数据即可输出判断是否为地震事件；

(4)利用正样本训练SSNet模型，得到用于地震震相到时的预测模型，进而向该模型输入地面震动波形数据即可输出得到P波和S波到达台站的时间。

进一步地，所述正样本包含地面震动波形数据以及对应实际监测到的一次地震事件中P波和S波到达台站的时间，负样本只包含地面震动波形数据；正样本的地面震动波形数据长度为30s，即从P波到达台站之前的5～15s截取起点，从P波到达台站之后的15～25s截取终点；负样本的地面震动波形数据长度也为30s，即从不涉及任何地震事件的连续地面震动波形记录中随机选取确定起点，截取30s长度的地面震动波形数据。

进一步地，所述步骤(2)中保证正样本和负样本的数量相同。

进一步地，所述步骤(3)中对SSNet模型进行训练的具体过程如下：

3.1构建SSNet模型，并初始化该模型参数，包括每一层的权值矩阵、学习率、优化方法以及最大迭代次数；

3.2将所有样本中的地面震动波形数据逐一输入至SSNet模型中进行训练，计算SSNet模型输出的识别结果与对应样本的真值之间的误差函数，其中模型输出的识别结果为样本属于地震事件的概率值，正样本的真值为1，负样本的真值为0；进而通过反向传播对整个SSNet模型的参数不断进行调整更新，直至误差函数收敛或达到最大迭代次数，从而完成训练得到用于地震事件检测的识别模型。

进一步地，所述步骤(4)中对SSNet模型进行训练的具体过程如下：

4.1构建SSNet模型，并初始化该模型参数，包括每一层的权值矩阵、学习率、优化方法以及最大迭代次数；

4.2将正样本中的地面震动波形数据逐一输入至SSNet模型中进行训练，计算SSNet模型输出的预测结果与对应样本的真值之间的误差函数，其中模型输出的预测结果为P波和S波到达台站的时间，样本的真值为实际监测到的一次地震事件中P波和S波到达台站的时间；进而通过反向传播对整个SSNet模型的参数不断进行调整更新，直至误差函数收敛或达到最大迭代次数，从而完成训练得到用于地震震相到时的预测模型。

进一步地，所述SSNet模型从输入到输出由尺度选择卷积块SSBlock1、平均池化层P1、尺度选择卷积块SSBlock2、平均池化层P2、尺度选择卷积块SSBlock3、平均池化层P3、尺度选择卷积块SSBlock4、平均池化层P4以及全连接层依次连接组成。

进一步地，所述尺度选择卷积块SSBlock1～SSBlock4结构相同均由三个尺度选择卷积单元SSCell1～SSCell3级联构成，每个尺度选择卷积单元包含普通卷积层和空洞卷积层，普通卷积层采用大小为3的卷积核对卷积单元的输入以步长为1进行卷积操作并经激活函数Relu处理后输出；空洞卷积层采用大小为3的卷积核对普通卷积层的输出以步长为1进行卷积操作，其输出分成两路分别经激活函数tanh和sigmoid处理后进行点乘，得到的结果与普通卷积层的输出相加并经批归一化处理后作为卷积单元的输出；SSCell1中空洞卷积层的空洞大小为1，SSCell2中空洞卷积层的空洞大小为2，SSCell3中空洞卷积层的空洞大小为4。

进一步地，所述尺度选择卷积块SSBlock1中普通卷积层和空洞卷积层所采用的卷积核数量均为32，尺度选择卷积块SSBlock2中普通卷积层和空洞卷积层所采用的卷积核数量均为64，尺度选择卷积块SSBlock3中普通卷积层和空洞卷积层所采用的卷积核数量均为128，尺度选择卷积块SSBlock4中普通卷积层和空洞卷积层所采用的卷积核数量均为256。

进一步地，所述平均池化层P1～P4结构相同均采用大小为4的池化窗口，步长为4，padding模式为same，即在窗口长度不足4时亦计算得到输出。

进一步地，所述全连接层采用三层神经网络结构，从输入至输出各层的神经元数量分别为256、64、16且均采用激活函数Relu处理。

基于上述技术方案，本发明根据地震波形数据的特点，针对性地设计了深度网络模型，综合利用卷积神经网络等技术抽取地震波形数据的特征，捕捉数据的特点。同时，本发明利用大量地震数据学习了一个地震事件到时识别模型，提高了对地震事件检测任务的准确率和震相到时拾取任务的精度。

附图说明

图1为本发明地震事件到时识别方法的流程示意图。

图2为本发明SSNet模型架构示意图。

图3为SSNet模型中SSBlock以及SSCell的结构示意图。

具体实施方式

为了更为具体地描述本发明，下面结合附图及具体实施方式对本发明的技术方案进行详细说明。

如图1所示，本发明基于SSNet模型的地震事件到时识别方法，具体步骤如下：

(1)收集地震监测台站记录的地面震动连续波形记录，用于地震事件到时识别。样本的特征维度为3×3000，其中3表示波形记录具有东西、南北以及竖直三个方向分量，3000的维度则是截取数据的时间长度为30s，数据的采样频率为100Hz。

(2)本发明采用如图2所示的SSNet模型进行地震事件到时识别，以卷积神经网络为基本方法，在此基础上加入空洞卷积、门控单元和残差结构等进行改进。SSNet模型由4个结构类似的SSBlock以及平均池化层堆叠组成，而每个SSBlock则是由三个相同的SSCell单元构成。其中，SSCell使用图3所示的卷积单元结构来替代简单卷积网络层的方案，SSCell除了常规的卷积网络层之外，还加入了空洞卷积层和门控单元和残差结构等改进，能够在计算中更好地抽取地震连续波形记录的特征。

由于使用了空洞卷积网络层，使得网络对于相对较宽的数据具有更大的视野；该层中采用了3×3的卷积，一个SSBlock里面的卷积的空洞大小依次增大，其大小为1、2、4，故其对于数据的视野大小依次为7、11、19。

门控单元层将空洞卷积的结果通过tanh和sigmoid两个不同的激活函数，之后对其结果进行点乘。由于地震数据样本的取值在数轴的两侧基本上是均匀分布的，如果使用常用的relu激活函数或其变种，那么将会损失掉一些数据中的负值信息，所以本发明采用了tanh(双曲正切函数)这一取值范围在[-1,1]的函数来作为主激活函数；而通过sigmoid函数的一路数据，起到的就是门控单元中的“开关信号”的作用，之后将这两路信号的对应位置相乘(Hadamard Product)。由于sigmoid函数的取值范围为[0,1]，所以相当于数据中的每一个值都经过了一定比例的衰减控制，最后再将门控单元的计算所得与捷径连接的数据进行拼接。

对于一个SSCell而言，假设其输入为矩阵为X，输出矩阵为Y，X与Y之间的计算关系为：

Conv_1＝conv1d(W1，X)

Conv_2＝dilated_conv1d(W2,Conv_1)

Xtanh＝tanh(Conv_2)

Xσ＝σ(Conv_2)

Y＝BN([Conv_1,Xtanh·Xσ])

上述公式中conv1d()表示一维卷积，dilated_conv1d表示带空洞的一维卷积，W1、W2为两次卷积的参数，σ()表示sigmoid函数，BN()表示归一化操作。

(3)设定模型参数，通过大量真实的地震事件到时记录，反向传播训练网络参数，得到最后的网络模型。

实验数据分为训练集、验证集与测试集，测试集是在数据预处理之前就预先划分出来的大小为总数据集20％的数据集，之后也经过了与训练集相同的数据预处理过程。在训练过程中，将数据分为训练集与验证集，其比例为4:1，实验使用的优化算法为Adam，学习率为0.001，批(batch)大小设置为128；实验不设置最大训练轮数或是迭代次数，采用了早停(EarlyStopping)的训练策略，当验证集上的损失在两个训练轮次(Epoch)中没有下降时，模型便会停止继续训练；在参数初始化上，使用Xavier初始化。

对于地震事件检测问题和震相到时拾取问题使用不同的损失函数。地震事件检测使用的是交叉熵，震相到时拾取使用的是平均平方损失(MSE)。

(4)对于地震事件检测问题，通过数据样本输入模型，通过预测层的Softmax激活函数，输出预测样本为地震事件的概率值，越大即表示该受测样本为地震事件的概率越大。

对于震相到时拾取问题，通过数据样本输入模型，最后输出P、S波到时数值即为P、S波到达时间与受测样本截取起点时间的差值。

以下通过AUC以及准确率衡量地震事件检测问题中的效果，表1为本发明方法与一种经典的地震事件检测的卷积神经网络方法的效果比较，从中可以看出本发明方法相对其在检测错误率方面有明显的降低。

表1

模型AUCAccErrorRateConvNetQuake0.999550.992920.00708SSNet0.999660.995140.00486

以下通过AUC以及准确率衡量震相到时识别问题中的效果，表2为本发明方法与现有的经典方法(包括传统方法和神经网络方法)的效果比较，从中可以看出本发明方法相对其在识别方面有明显的降低。

表2

上述的对实施例的描述是为便于本技术领域的普通技术人员能理解和应用本发明。熟悉本领域技术的人员显然可以容易地对上述实施例做出各种修改，并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此，本发明不限于上述实施例，本领域技术人员根据本发明的揭示，对于本发明做出的改进和修改都应该在本发明的保护范围之内。

一种基于SSNet模型的地震事件到时识别方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。