基于UNet++联合CliqueBlock的微地震P波初至拾取的方法

基于unet 联合clique block的微地震p波初至拾取的方法
技术领域
1.本发明属于地震监测技术领域，具体涉及一种基于unet 联合clique block的微地震p波初至拾取的方法。


背景技术：

2.微地震监测技术就是通过观测、分析由压裂产生的微小地震事件来监测地下状态的地球物理技术，对油田开发稳产、高产具有重要的意义。微地震有效信号能量较弱，信噪比较低，甚至完全淹没在噪音之中。因此，如何准确、快速的拾取低信噪比微地震信号p波初至对震源定位、裂缝预测、震源破裂机制分析具有重要的意义。
3.传统拾取方法(如sta/lta、aic准则、相关检测)等虽然取得了一定效果，但效率难以满足大数据时代下微地震信号初至拾取的需求。以卷积神经网络为代表的深度学习以其层数深、容量大的特点，在高精度、高效率处理海量数据上的优势日益凸显。
4.为提高微地震p波初至拾取的精度与效率，本专利提出一种unet 联合clique block的微地震p波初至拾取方法。


技术实现要素：

5.针对现有微地震初至拾取技术难以高效、准确地拾取海量微地震p波初至的不足，本发明利用unet 联合clique block构建全新网络mspicking-net，以进一步增强低信噪比微地震信号信噪可辨，提高微地震初至拾取的精度；为实现上述技术效果，本发明所采用的技术方案是：
6.一种基于unet 联合clique block的微地震p波初至拾取的方法，包括如下步骤：
7.s1：对压裂区域进行井下或者地面观测，采集水力压裂微地震监测数据；监测的微地震监测数据信噪比较低、有效信号能量微弱；
8.s2：利用有限差分波动方程正演方法生成不同主频、不同地层速度模型下的正演模拟信号；
9.其中有限差分波动方程正演方法选用论文《盛冠群.微地震信号检测与震源定位方法研究[d].中国石油大学(华东),2016.》中记载的交错网格有限差分正演方法；
[0010]
s3：对s2得到的正演模拟信号设置不同信噪比，连同实际监测剖面构建原始数据集，对所述原始数据集进行初至拾取，选取初至和非初至处的信号波形并分别进行标定；
[0011]
s4：在unet 网络中添加clique block块，构建mspicking-net网络；其中clique block块的数量以及mspicking-net的超参数设定均通过交叉验证方法进行确定；
[0012]
优选地，交叉验证方法选用loocv或者k-fold cross validation方法进行验证；
[0013]
在clique block中，除了输入节点外，同一块中的任何两层都是双向连接的；在一个block中，每一层即是其他层的输入，也是其他层的输出；clique block中，每一层都是前一层的输出，后一层的输入是任何其他层的输出；具有5层的团块的传播如下表1所示：
al.2018.unet :a nested u-net architecture for medical image segmentation》中的unet 网络技术；clique block选用论文《yang,y.,et al."convolutional neural networks with alternately updated clique."2018ieee/cvf conference on computer vision and pattern recognition(cvpr)ieee,2018.》中的技术。
[0029]
本发明的有益效果为：
[0030]
1、针对传统微地震初至拾取在海量微地震信号处理方面效率低下的劣势，本专利利用深度学习指导微地震信号初至拾取，实现了高效的海量微地震信号初至拾取；
[0031]
2、针对微地震信号的低信噪比环境及unet 层数较浅、难以捕获精细微地震信噪特征等缺陷，联合unet 网络和clique block构建全新网络mspicking-net，实现了深层精细微震信噪特征抽取，强化微地震低信噪比信号信噪分离，实现更快速、精细的低信噪比微地震初至拾取，进一步增强了低信噪比微地震信号信噪可辨，提高微地震初至拾取的精度。
[0032]
3，本发明为微地震监测技术指出了一个重要的应用方向，促进该技术的应用和推广，并且通过将人工智能技术应用于传统微地震监测领域，以奠定微地震智能处理算法的研究基础，更可以为服务于国家的能源战略安全提供重要的理论基础。
附图说明
[0033]
图1为本发明的工作流程示意图；
[0034]
图2为本发明中mspicking-net网络结构示意图；
[0035]
图3为本发明基于mspicking-net微地震信号初至拾取结果中信噪比为snr＝-10db的待检测信号图；
[0036]
图4为本发明基于mspicking-net微地震信号初至拾取结果中u-net 网络有效信号初至位置概率预测曲线；
[0037]
图5为本发明基于mspicking-net微地震信号初至拾取结果中mspicking-net网络有效信号初至位置概率预测曲线；
具体实施方式
[0038]
实施例1：
[0039]
如图1所示，一种基于unet 联合clique block的微地震p波初至拾取的方法，包括如下步骤：
[0040]
s1：对压裂区域进行井下或者地面观测，采集水力压裂微地震监测数据；监测的微地震监测数据信噪比较低、有效信号能量微弱；
[0041]
s2：利用有限差分波动方程正演方法生成不同主频、不同地层速度模型下的正演模拟信号；
[0042]
其中有限差分波动方程正演方法选用论文《盛冠群.微地震信号检测与震源定位方法研究[d].中国石油大学(华东),2016.》中记载的交错网格有限差分正演方法；
[0043]
s3：对s2得到的正演模拟信号设置不同信噪比，连同实际监测剖面构建原始数据集，对所述原始数据集进行初至拾取，选取初至和非初至处的信号波形并分别进行标定；
[0044]
s4：如图2所示，在unet 网络中添加clique block块，构建mspicking-net网络；其中clique block块的数量以及mspicking-net的超参数设定均通过交叉验证方法进行确
定；mspicking-net网络结构图如图2所示；
[0045]
优选地，交叉验证方法选用loocv方法进行验证；
[0046]
在clique block中，除了输入节点外，同一块中的任何两层都是双向连接的；在一个block中，每一层即是其他层的输入，也是其他层的输出；clique block中，每一层都是前一层的输出，后一层的输入是任何其他层的输出；具有5层的团块的传播如下表1所示：
[0047][0048]
表1：
[0049]
在第一阶段，输入层(x0)通过单向连接初始化该块中的所有层，将连接每个更新的层以更新下一层；从第二阶段开始，层开始交替更新，除了要更新的顶层之外的所有层都连接为底层，并且它们相应的参数也连接在一起；因此，第i层(i≥1)的输出在第k层(k≥2)循环中，可表示为：
[0050][0051]
其中，*表示利用超参数w进行卷积操作，g表示非线性激活函数；w
ij
在不同阶段保持重用；
[0052]
s5：训练集制作与网络训练：随机抽取s3得到的样本中等数量的初至位置和非初至位置信号得到训练集，完成mspicking-net训练；
[0053]
s6：输出初至位置：将待测信号输入网络，利用softmax预测待测信号每一个采样点属于初至类别和非初至类别的“二分类”概率，最后以初至类别中概率峰值处对应的点作为初至，并利用正确率与召回率曲线进行算法评价；softmax公式如下所示：
[0054][0055]
损失函数由真实值的概率分布pi(x)与网络预测概率分布qi(x)之间的交叉熵c(p,q)计算得到：
[0056][0057]
其中，x表示全卷积网络最后一层输出f(x)的每一个点；
[0058]
i表示类别：其中1表示初至，2表示非初至；
[0059]
qi(x)表示x分别从属不同类别的预测概率。
[0060]
优选地，上述基于unet 联合clique block的微地震p波初至拾取的方法，步骤s2中不同主频、不同地层速度模型包括mamousi模型、多层平层和斜层速度模型；不同速度模型生成的正演模拟信号各不相同，用来保证样本的足够的丰度。
[0061]
实施例2：
[0062]
如图3～图5所示，是实际应用中基于mspicking-net微地震信号初至拾取结果的示意图；
[0063]
图3为添加了-10db高斯噪音的正演模拟信号，该信号初至时间为1330ms，主频25hz；图4为u-net 网络有效信号初至位置概率预测曲线；图5为mspicking-net网络有效信号初至位置概率预测曲线；由图4和图5可知，unet 网络拾取的微地震初至时间为1480ms，而本专利所提出的方法拾取的微地震初至为1330ms，与模拟信号初至时间完全吻合；因此，相较于unet 网络，本专利对低信噪比方案的微地震信号初拾取结果更加可靠，再次证明了，本专利实际上实现了针对微地震信号的低信噪比特性和unet 网络的缺陷进行优化的技术效果，本专利提出的mspicking-net网络在实际使用中性能更加优越。
[0064]
上述的实施例仅为本发明的优选技术方案，而不应视为对于本发明的限制，本技术中的实施例及实施例中的特征在不冲突的情况下，可以相互任意组合。本发明的保护范围应以权利要求记载的技术方案，包括权利要求记载的技术方案中技术特征的等同替换方案为保护范围。即在此范围内的等同替换改进，也在本发明的保护范围之内。