CN115932959A - 一种基于UNet的面波频散质量快速评估方法 - Google Patents

Abstract

本发明属于地震波数据处理领域，公开了一种基于UNet的地震面波频散质量快速评估方法，其首先采集面波频散数据，对每个频散点人为标记合格标签，并将其分成训练集样本和测试集样本；然后构建面波频散质量评估的UNet模型，并利用训练集样本对UNet模型进行训练；最后将周期、面波频散、以及频散的变化率作为训练好的UNet模型的输入，得到质量评估后的面波频散的UNet预测结果。本发明所述的方法能够在短时间内完成批量频散数据的处理，不仅能够作为传统方法的替代，且相对于传统的人工识别方法极大地提升了数据处理的效率，为面波频散的实时反演奠定基础。

Description

一种基于UNet的面波频散质量快速评估方法

技术领域

本发明属于地震波数据处理领域，具体是涉及一种基于UNet的地震面波频散质量快速评估方法。

背景技术

面波频散方法被广泛地应用于地壳结构探测、工程与环境勘探领域。来自不同周期的面波频散数据对不同深度的速度结构具有不同的敏感性，短周期的频散数据对浅部结构敏感而长周期的频散对深部的结构敏感；因此，面波频散曲线分析是研究壳幔结构的有效方法。通过反演频散曲线，可以得到一维S波速度结构，二维和三维的速度结构可以由一维速度结构组合起来。通过反演得到可靠的速度结构的重要前提是提取可信赖的频散曲线。目前主要通过两个地震台站获得的地震事件或者噪声数据进行互相关方法来获得面波频散，然后将工区划分为多个网格，使用走时层析成像技术可以得到每个网格上的频散曲线。

面波频散的提取主要包含四个步骤(以双台法为例)：1.对两个台站的波形进行滤波，然后对滤波后的地震波行进行互相关计算，并且得到相关时差，并进一步计算出某个频率下两个台站之间的地震波速度；2.以频率为横坐标，以速度为纵坐标，以互相关值为能量值，绘制出面波的二维速度谱；3.人工挑选(或者计算机自动识别)出速度谱上能量最大区域，形成一条曲线，即频散曲线；4.人工质量筛选，由于计算误差，或者原始波形的信噪比较低等问题，通常需要对得到的频散曲线进一步筛选，得到光滑的，速度值符合基本地质规律的频散曲线。在上述的步骤3中，由于存在计算不稳定，精度较低等问题，或者噪声的存在，提取的相速度可能存在不连续现象，且频散提取的过程需要人机交互进行频散拾取，这不可避免地引入了主观因素；某些自动频散拾取技术也可能导致部分周期相速度不合理的情况。

发明内容

为解决上述技术问题，本发明提供了一种基于UNet的面波频散质量快速评估方法，通过实际数据建立起训练集和测试集，然后构建了UNet实现对面波频散质量进行自动评估，极大地提升了数据处理的效率。

本发明所述的一种基于UNet的面波频散质量快速评估方法，步骤为：

步骤1、采集面波频散数据，对每个频散点人为标记合格标签，并将其分成训练集样本和测试集样本；

步骤2、构建面波频散质量评估的UNet模型，并利用训练集样本对UNet模型进行训练；；

步骤3、将面波频散的周期、面波频散、以及频散的变化率作为训练好的UNet模型的输入，得到质量评估后的面波频散的UNet预测结果。

进一步的，步骤1中，每个训练样本需要人为指定合格的区域(标签)，将合格的频散点标记为1，不合格的标记为0。

进一步的，所述UNet网络模型左右各4层，一共包含15个卷积层，3个池化层，3个转置卷积层；选取RELU激活函数作用于卷积层的输出；卷积操作之后使用Batch_normilization对数据进行标准化。

进一步的，所述UNet模型的输入层包括三个通道：周期，面波频散，以及频散的变化率；输入层一共48个神经元，不足48个频散点的数据充零；输出层为每个周期点对应的标签，质量合格标签为1，否则标签为0；

为了能够更好的使输出层与标签匹配，对输出层作用于sigmod激活函数，使得最终输出神经元的值在0～1之间，sigmod激活函数的数学表达为：

σ(x)＝1/(1+exp(-x))

其中，x为输入。

进一步的，UNet训练的损失函数需要极小化向量差的二范数：

其中，下标i表示第i个训练样本，N为训练样本的总数；D为由深度学习实现的预测结果；label_i为人为给定的标签，表征每个采样点数据的质量，d_i为输入频散数据，其包含三个通道；θ为网络反向传播中需要更新的超参数。

本发明所述的有益效果为：本发明提出了利用UNet来对地震面波频散数据进行快速评估，通过实际数据建立起训练集和测试集，然后构建了UNet模型实现对面波频散质量进行自动评估。本发明所述的方法简洁高效，可以在模型训练好之后进行快速预测，在短时间内完成批量频散数据的处理，实现对面波频散曲线的质量进行快速评估，剔除掉不符合要求的数据部分，保留下优质的数据。该方法不仅能够作为常规方法的替代，且相对于传统的人工识别方法极大地提升了数据处理的效率，为面波频散的实时反演奠定基础。

附图说明

图1是采集面波频散数据的台站分布示意图；

图2是本发明所述的UNet网络模型结构示意图；

图3是面波频散的UNet预测结果示意图；

图4是测试集的预测结果示意图。

具体实施方式

为了使本发明的内容更容易被清楚地理解，下面根据具体实施例并结合附图，对本发明作进一步详细的说明。

本发明所述的一种基于UNet的面波频散质量快速评估方法，具体步骤如下：

1.数据准备

为了得到面波频散数据，目前主要通过两个地震台站获得的地震事件或者噪声数据进行互相关。然后将工区划分为多个网格，使用走时层析成像技术可以得到每个网格上的频散曲线。本发明的面波频散通过噪声计算得到，一共采集了2013年12月至2015年3月(15个月)的垂直分量连续时间序列，由668个宽带和超宽带便携式台站(中国地震台阵，二期)和62个永久宽带台站记录(如图1所示)。通过互相关等方法估计经验格林函数，最后提取出4160对频散数据用于质量评估的训练(Pan et al.,2019)。每个训练样本需要人为指定合格的区域(标签)，将合格的频散点标记为1，不合格的标记为0。其中的85％用于训练集，15％用于测试集。最后可以将筛选之后的频散结果进行层析成像，获得整个区域的频散。

2.UNet构建

面波频散的挑选需要对每一个频散数据的每一个频点质量进行判断，只有当某条频散曲线所有周期的数据都较差，才全部舍去；因此，面波频散的质量评估不是一个简单的二分类问题。这里选用UNet神经网络来实现接收函数的质量评估，模型的结构如图2所示。网络左右各4层，一共包含15个卷积层，3个池化层，3个转置卷积层。由于面波频散数据中没有负数，选取RELU激活函数作用于卷积层的输出，该激活函数能够提高网络的非线性能力，同时避免了训练过程中梯度爆炸和梯度消失问题。卷积操作之后使用了Batch_normilization对数据进行标准化，以进一步防止梯度消失或梯度爆炸现象，同时可以增加正则化效果。

输入层包括三个通道：面波频散的周期，面波频散，以及频散的变化率。输入层一共48个神经元，不足48个频散点的数据充零，这意味着网络适用于具有任意周期的低于48个频散数据的质量评估。输出层为每个周期点对应的标签，质量合格标签为1，否则标签为0。为了能够更好的使输出层与标签匹配，对输出层作用于sigmod激活函数，使得最终输出神经元的值在0～1之间，sigmod激活函数的数学表达为：

σ(x)＝1/(1+exp(-x))

其中，x为输入。

UNet训练的损失函数，需要极小化向量差的二范数：

这里下标i表示第i个训练样本，N为训练样本的总数；D为由深度学习实现的预测结果；label_i为人为给定的标签，表征每个采样点数据的质量，d_i为输入频散数据，其包含三个通道；θ为网络反向传播中需要更新的超参数。

在训练过程中，使用反向传播算法Adam(Kingma and Ba,2014)来更新神经元的权值，学习率为0.001。根据目标函数随迭代次数的下降曲线，选择迭代70次以保证目标函数充分下降。

3.质量预测

模型训练好之后，将需要预测的数据按照训练集输入的格式准备好，即包括周期，频散数据以及频散的时间导数；然后导入UNet的输入层，可以快速得到质量评估结果，每个周期点的质量评估值位于0～1区间内。此时可以设定一个阈值(例如0.8)，高于该阈值保留频散值，反之删去频散值，值越高表明网络对质量的要求更高。最后可以得到质量较好的频散曲线，删去了质量较差的部分。

图3展示了三个面波频散数据的质量评估结果，其中a～d图分别展示了整体曲线质量较好，中间部分质量较好、后段质量较好、整体曲线质量较差的评估结果。每个子图的上图表示预测质量较好的频散数据以及原始频散数据对比结果，下图表示预测结果，平直折线为人为给定标签，带竖棒的虚线表示UNet预测结果，越接近于1表示质量越好，竖棒长度表示不确定性。可以看出，UNet预测值较大的区域对应着相对连续，平滑的频散曲线部分，是合格的标准频散曲线形态。根据测试结果，设定阈值0.8可以较好地筛选出合格的频散曲线，即预测结果大于0.8被认为质量合格。总体上来看，UNet预测结果与人为给定的数据吻合良好，说明了UNet筛选面波频散曲线的有效性。但是在图3中，面波频散数据的变化趋势较大，存在较大的不确定性，UNet预测整段数据不符合质量要求，虽然与人为给定标签存在差异，但是也可以认为UNet的预测结果是合理的。图4显示了全部测试集的预测结果，其中，图4(a)为原始数据，图4(b)为UNet的预测结果，可以看出，UNet筛选后的频散曲线去掉了大部分高频的不连续曲线，整体上表现得更加平滑，渐变，更加符合实际情况。

这说明了本发明所提出的利用CNN评估面波频散质量是可行的，该方法可以在短时间内对面波频散曲线的质量进行快速评估，剔除掉不符合要求的数据部分，保留下优质的数据。该方法不仅可以替代人工挑选频散曲线，且极大的提高了面波频散数据的处理效率。

以上所述仅为本发明的优选方案，并非作为对本发明的进一步限定，凡是利用本发明说明书及附图内容所作的各种等效变化均在本发明的保护范围之内。

Claims (5)

1.一种基于UNet的面波频散质量快速评估方法，其特征在于，所述方法步骤为：

步骤1、采集地震面波频散数据，对每个频散点人为标记合格标签，并将其分成训练集样本和测试集样本；

步骤2、构建面波频散质量评估的UNet模型，并利用训练集样本对UNet模型进行训练；；

步骤3、将面波频散的周期、面波频散、以及频散的变化率作为训练好的UNet模型的输入，得到质量评估后的面波频散的UNet预测结果。

2.根据权利要求1所述的一种基于UNet的面波频散质量快速评估方法，其特征在于，步骤1中，每个训练样本需要人为指定合格的标签，将合格的频散点标记为1，不合格的标记为0。

3.根据权利要求1所述的一种基于UNet的面波频散质量快速评估方法，其特征在于，所述UNet模型左右各4层，一共包含15个卷积层，3个池化层，3个转置卷积层；选取RELU激活函数作用于卷积层的输出；卷积操作之后使用Batch_normilization对数据进行标准化。

4.根据权利要求3所述的一种基于UNet的面波频散质量快速评估方法，其特征在于，所述UNet模型的输入层包括三个通道：周期，面波频散，以及频散的变化率；输入层一共48个神经元，不足48个频散点的数据充零；输出层为每个周期点对应的标签，质量合格标签为1，否则标签为0；

为了能够更好的使输出层与标签匹配，对输出层作用于sigmod激活函数，使得最终输出神经元的值在0～1之间，sigmod激活函数的数学表达为：

σ(x)＝1/(1+exp(-x))

其中，x为输入。

5.根据权利要求4所述的一种基于UNet的面波频散质量快速评估方法，其特征在于，UNet训练的损失函数需要极小化向量差的二范数：

其中，下标i表示第i个训练样本，N为训练样本的总数；D为由深度学习实现的预测结果；label_i为人为给定的标签，表征每个采样点数据的质量，d_i为输入频散数据，其包含三个通道；θ为网络反向传播中需要更新的超参数。

Images