基于卷积降噪自编码器的地震数据去噪

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

2020年12
月第55卷 第6期 
*湖北省武汉市蔡甸区蔡甸街大学路111号长江大学地球物理与石油资源学院,430100。

Email:cwjycd@163.com本文于2019年11月11日收到,最终修改稿于2020年9月8日收到。

本项研究受国家自然科学基金项目“基于经验模态分解的自由表面多次波衰减方法研究”(41804140
)和“碳酸盐岩不同孔隙结构多尺度三维数字岩心建模方法研究”(41674136
)联合资助。

·智能地球物理·
文章编号:1000-7210(2020)06-1210-
10基于卷积降噪自编码器的地震数据去噪
宋 辉①② 高 洋③ 陈 伟*①④ 张 翔①②
(①油气资源与勘探技术教育部重点实验室(
长江大学),湖北武汉430100;②长江大学地球物理与石油资源学院,湖北武汉430100;③中国石油大学(
北京)CNPC物探重点实验室,北京102249;④非常规油气湖北省协同创新中心,湖北武汉430100)宋辉,高洋,陈伟,张翔.基于卷积降噪自编码器的地震数据去噪.石油地球物理勘探,2020,55(6):1210-1219.摘要 噪声压制是地震勘探中一个长期存在的问题,虽然一些传统方法能够压制数据中的噪声,但存在有效信号丢失、噪声残留等问题。

为此,提出了一种基于卷积降噪自编码器的无监督地震数据去噪算法。

该算法首先对地震数据进行一定程度的随机损坏,然后将损坏后的地震数据输送到编、解码框架。

编码框架负责捕捉地震数据波形特征,据此消除噪声;解码框架能够对特征图进行扩大并恢复地震数据细节信息,从而得到重构的地震数据。

最后,将重构地震数据与原始地震数据之间的误差作为收敛代价进行模型训练。

考虑到地震数据的复杂性与特殊性,在编码和解码阶段使用了多尺度卷积模块提取地震数据特征。

合成数据与实际数据的验算结果表明,该方法在保护地震信号的同时能够有效压制随机噪声、提高地震信号的信噪比。

关键词 无监督学习 卷积神经网络 降噪自编码器 地震数据 去噪
中图分类号:P631 文献标识码:A doi:10.13810/j
.cnki.issn.1000-7210.2020.06.0060 引言
随着油气勘探环境变得愈加复杂,
地震数据在采集过程中受到的噪声干扰也愈加严重[
1]。

地震数据中的噪声可分为两类:
相干噪声和随机噪声。

相干噪声具有一定的波形特征,而随机噪声没有固定的波形特征,与有效信号随机混合在一起,影响高分辨率处理、属性分析、反演等。

噪声压制在地球物理学中是一个经典的问题,已提出了很多方法。

传统的噪声压制方法是基于地震数据的特性,如可预测性、稀疏性等。

基于预测性的去噪方法利用有效信号线性或拟线性生成预测滤
波器提高信噪比[2-
3],但不能处理非线性地震信号。

基于稀疏变换的去噪方法利用信号在变换域内具有
较好的稀疏性压制噪声,例如小波变换[4-
5]、曲波变换[6-
7]、Shearlet变换[8-9]等。

由于它们建立在固定的
变换基函数上,因此不能自适应地处理结构复杂的地震数据。

字典学习是一种基于数据驱动的算法,使用可学习的字典代替传统的变换基函数,因此可以自适
应地表示数据[
10-
11]。

中值滤波在图像处理中是一种经典的平滑去噪方法,利用信号与噪声之间的统计差
异抑制噪声,常用于地震数据去噪[
12
]。

以上方法噪声压制效果受限于模型假设和参数设置等因素。

近年来,深度学习在图像处理、语音处理等领域
取得了显著成果[
13-
15]。

通过构建多个处理层逐步实现抽象的特征表示,可以完成复杂的分类或预测等
任务[16-
19]。

目前,基于标签数据的深度学习方法已经成功应用于地震随机噪声压制[20]
,根据生成标签
数据的策略,主要可以分为两类:①将传统去噪方法的结果作为标签数据训练神经网络,然后将训练
好的网络模型用于实际地震数据噪声压制[
21-
22];②将创建的合成地震数据作为标签数据训练神经网
络[23
],然后将训练好的网络模型用于实际地震数据噪声压制[
24-
25]。

第一类方法制作的标签数据并不太准确,因此会影响网络模型的去噪性能。

第二类方法制作的标签数据虽然准确,但是合成数据与实际地震数据差异大,影响网络模型的泛化能力。

基于标签数据的深度学习方法取得了良好的去噪效果,但是在制作复杂而庞大的训练集上也会浪费大量的资源。


 第55卷 第6期宋辉,等:基于卷积降噪自编码器的地震数据去噪

211 此,开发不依赖于标签数据的无监督学习方法对于地震数据噪声压制具有重要意义。

为了改善地震数据的噪声压制效果,本文提出了一种基于卷积降噪自编码器的地震噪声压制算法。

该算法基于无监督学习,能够直接从含噪的地震数据中恢复无噪的地震数据,省去了标注数据所需的大量工作。

1 理论
1.1 自编码器
自编码器属于无监督学习,它利用反向传播算法对无标签的输入数据进行重构。

典型的自编码器如图1所示,可以看出它包括两部分:编码框架和解码框架。

编码框架负责将输入数据映射为潜在空间的特征表达,解码框架则负责将特征表达解码为与原始输入数据相近的输出数据。

自编码器的目标是使输出数据尽可能等于其输入数据,但是它不能有效提取地震数据特征。

对于无标签的输入数据,自编码器按以下方式重建输入数据
H=ξθ1
(P)=σ(W1P+b1)(1
)Q=ξθ2
(H)=σ(W2H+b2)(2
)式中:P、H与Q分别为输入数据、
特征表达、输出数据;σ为sigmoid激活函数;W1与b1分别为输入层与隐藏层之间的权重矩阵与偏置;W2与b2分别
为隐藏层与输出层之间的权重矩阵与偏置;θ1=[W1,b1]和θ2=[W2,b2]
分别为编码参数和解码参数;ξθ1(P)和ξθ2
(H)分别为编码函数和解码函数。

模型训练的目的是为了优化模型参数[θ1,θ2]
,图1 自编码器模型示意图
使重建数据Q与输入数据P尽可能地接近,即
[θ1,θ2]=arg minθ1,
θ2
L P,ξθ2[ξθ1(P{})](3
)式中L(·)表示L2范数,用于衡量Q与P之间的重构误差。

1.2 降噪自编码器
降噪自编码器是自编码器的一个变体。

与自编码器不同的是,降噪自编码器通过训练损坏的输入数据进行特征学习,其核心思想是提取数据的鲁棒性特征。

典型的降噪自编码器如图2所示,可以看出,降噪自编码器与自编码器对输入层的处理不同,降噪自编码器按照一定的概率将输入节点置0,如果这种损坏的概率为0,降噪自编码器就退化为自编码器。

图2 降噪自编码器模型示意图
为了定性地比较这两种自编码器的学习效果,分别设计了对应的两种自编码器网络,网络层数设置为3,其中隐藏层单元数为64,并将这两种自编码器输入层与隐层之间的权重矩阵进行可视化,结果如图3所示。

从图3可以看出,自编码器训练后得到的权重是杂乱的,含有大量噪声,而降噪自编码器训练后得到的权重含有明显的结构特征,噪声较弱(这两种自编码器的输入数据均来源于本文实验部分的合成数据,权重尺寸为48×48,降噪自编码器的损坏程度为60%)。

可以看出,降噪自编码器能够提取、编码出具有鲁棒性的特征,具有很强的网络
表达能力。

 石油地球物理勘探2020年 1212
图3 自编码器(a)和降噪自编码器(b)的权重矩阵可视化结果
 第55卷 第6期宋辉,等:基于卷积降噪自编码器的地震数据去噪

213 1.3 卷积降噪自编码器
全连接形式的降噪自编码器的输入是一维形式,因此对图像处理存在局限性。

在图像中,局部的像素联系较为紧密,相距较远的像素关联性不大,而全连接神经网络的结构里下层神经元与所有上层神经元都能够形成连接,促使网络学习图像全局结构,忽略了图像的局部相关性,同时造成了参数糅杂。

近年来,卷积神经网络在图像处理方面应用广泛,这归因于它能够通过卷积核捕获原始图像中的局部特征,并且卷积核权重共享,大大减少了参数量。

卷积形式的降噪自编码器在图像处理方面已经取得了不错的效果,但还未被用于地震数据去噪。

本文设计的卷积降噪自编码器如图4所示。

使用3层卷积层和池化层作为编码框架,3层上采样层和卷积层作为解码框架。

在编码框架中,卷积层作为特征提取层,用于捕捉地震数据波形特征,而池化层作为特征压缩层,一方面能够减小特征图的尺寸,降低网络计算量;另一方面能够提取重要的地震数据特征,有效降低噪声成分。

本文将每层卷积层设置24个卷积核,卷积核的移动步长设置为1。

因此,48×48×1的地震数据经过编码框架以后就被压缩为6×6×24的压缩特征表达。

压缩特征表达保留了地震数据最重要的信息,但是丢失了大量细节信息。

因此,解码框架承担着扩大特征图与恢复地震数据细节信息的任务。

图4 卷积降噪自编码器模型
地震数据的特征复杂,
具有多尺度特征,而单一尺寸的卷积核只能学习到特定尺度的特征。

因此,本文设计了多尺度卷积模块,一个多尺度模块相当于三个卷积核,如图5所示。

多尺度卷积模块就是将不同大小的卷积核构成的网络并联起来,既有效解决了卷积核大小的选择问题,又能够有效捕捉地震数据的复杂特征。

图5 多尺度卷积模块
2 实验
本文使用数据分块的训练方式去噪。

即采用一个固定尺寸的窗口对地震数据进行滑移,窗口按固
定的滑移步长每滑移一次,产生一个样本数据。

窗口尺寸代表了样本数据的尺寸,窗口尺寸一旦确定,滑移步长就决定了样本数据的数量。

首先,分块训练的方式较好地适应了地震数据的关键信息只与局部数据有关,避免了在训练过程中计算资源的浪费;其次,该训练策略会产生大量的样本数据,有利于训练出良好的模型参数。

在本文实验中,将窗口尺寸设置为48×48,滑移步长设置为4。

由于无监督学习的性质,本文实验不需要单独的训练集训练模型参数,模型参数的更新是通过对原始含噪数据的自我学习完成,利用早停止原则控制模型训练的结束。

使用合成地震数据与实际地震数据测试本文方
法的去噪性能,并与多道奇异值分析(MSSA)
[26]
、小波变换[27]
、f-
x反褶积[3]等方法进行对比。

其中小波变换采用bior2.4小波基,f-
x反褶积预测滤波器的长度设置为6,MSSA要保留的奇异值数量设置为5。

2.1 合成地震数据
合成地震数据(图6a)共88道,每道468个样
 1214 石油地球物理勘探2020年 
点,采样间隔为1ms。

该合成数据包含线性同相轴、曲线同相轴、间断同相轴以及断层;当加入噪声后(图6b),同相轴的连续性变差。

使用信噪比(SNR)定量地评价不同方法的去噪性能
SNR=20lg
‖dcle
an‖2‖ddenoised-dclean‖
2(4
)式中dclean和ddenoise分别表示无噪地震数据与去噪后的地震数据。

加噪合成地震数据的SNR为2.04dB。

由于网络结构的复杂性,
首先测试卷积层数、卷积核数量、卷积核尺寸以及输入数据的损坏程度对去噪结果的影响。

(1
)卷积层数。

将卷积层数从2增加到8,合成数据去噪后的SNR测试结果如图7a所示。

可见随着层数的增加,信噪比先增加后降低;当卷积层数为6时,网络去噪性能最佳,因此将卷积层数设置为6。

图6 合成地震数据(a)及其加噪结果(b
) (
2)卷积核数量。

将卷积核数量从6增加到72,合成数据去噪后的SNR测试结果如图7b所示。

可见,当卷积核数量为24时网络去噪性能最佳,因此将卷积核数量设置为24。

(3
)卷积核尺寸。

将卷积核尺寸设置为3×3、5×5、7×7或三者的组合,合成数据去噪后的SNR
测试结果如图7c所示。

可见,多尺度的卷积核能够获得更优的去噪结果,因此本文采用多尺度卷积核。

(4
)输入数据的损坏程度。

本文测试了输入数据的损坏程度对去噪结果的影响,其结果如图7d所示。

由图7d可知,对输入数据进行一定程度地损坏可以实现更好的去噪性能。

根据测试结果,本文将输入数据的损坏程度设置为60%。

上述四种网络参数确定之后,运行该网络可获得最终的去噪结果,并与MSSA、小波变换和f-x反褶积等去噪方法进行比较(
图8)。

由图8可见:MSSA、f-
x反褶积法去噪不足,去噪剖面上含有明显的残余噪声,且MSSA法去噪剖面的间断点不清晰;小波变换法去噪剖面上无明显残留噪声,但同相轴能量不稳定;本文方法的去噪剖面上基本无噪声残留、同相轴能量更稳定、间断点更清晰。

为了比较四种去噪方法的保幅性,计算了四种去噪方法的残差剖面(图9)。

小波变换与f-x反褶
积法的残差剖面上有明显的有效信号痕迹,表明有效信号损失严重,而本文去噪方法的残差剖面上有效信号的痕迹不明显,有效信号损失最小,说明本文方法有良好的保幅效果。

从平均振幅谱(图10)上可以看出,本文方法去噪结果的平均振幅谱最接近于原始信号的平均振幅谱,所以本文方法在去除噪声的同时能够最大程度地保护有效信号不受损失。

为了定量比较四种方法的去噪性能,分别计算了不同方法去噪结果的信噪比。

MSSA、
小波变换、f-
x反褶积和本文方法的去噪结果的信噪比分别为11.46、12.26、12.41、16.32dB。

图11为在不同噪声水平下的四种方法去噪结果的信噪比,可见四种方法去噪之后信噪比都会得到一定的提升,但本文方法始终保持最高的信噪比,说明本文方法有更强的鲁棒性。

2.2 实际地震数据
选取一段明显被随机噪声污染的实际地震剖面(图12)测试各种方法的去噪能力。

该地震数据共120道地震道,单道含548个采样数,
采样间隔为
第55卷 第6期宋辉,等:基于卷积降噪自编码器的地震数据去噪1215
图7 网络参数对去噪性能的影响
(a)卷积层数;(b)卷积核数量;(c)卷积核尺寸;(d)输入数据损坏程度
图8 合成数据不同方法的去噪结果对比
(a)MSSA法;(b)小波变换法;(c)f-x反褶积法;(d)本文方法
 石油地球物理勘探2020年 1216
图9 合成数据不同方法去噪的残差剖面
(a)MSSA法;(b)小波变换法;(c)f-x反褶积法;(d)本文方法
图10 合成数据不同方法去噪结果的平均振幅谱
(a)MSSA法;(b)小波变换法;(c)f-x反褶积法;(d)本文方法
 第55卷 第6期宋辉,等:基于卷积降噪自编码器的地震数据去噪

217 
 图11 不同噪声水平下四种方法去噪结果的信噪比
1ms。

由图12可以看出,该地震数据含有曲线同相轴以及断层构造,但受噪声干扰严重,同相轴连续性较差,弱信号识别困难。

使用与合成数据去噪相同的网络模型,并重新训练用于实际地震数据去噪。

MSSA、小波变换、f-
x反褶积和本文方法去噪所耗时间分别为0.22s、0.07s、0.08s、515.86s。

图13为四种方法实际地震数据的去噪结果。

小波变换法去噪后的地震剖面分辨率低,同相轴不清晰。

与其他方法相比,f-
x反褶积法去噪后的地震剖面中残留的噪声更多。

MSSA法的去噪结果中无明显噪声残留,但损失了某些局部细节。

本文方法去噪结果含有更丰富的局部细节,无明显噪声残留,同相轴清晰连续。

图14为四种方法对应的残差剖面,MSSA、小波变换和f-x反褶积方法的残差剖面上含有明显的有效信号的痕迹,因此这三种传统去噪方法在压制随机噪声的同时损失了部分有效信号。

而本文方
 图1
2 实际地
震剖面图13 实际地震剖面四种方法的去噪结果
(a)MSSA法;(b)小波变换法;(c)f-
x反褶积法;(d)本文方法
 1218
 石油地球物理勘探2020年
图14 实际地震数据四种去噪方法的残差剖面(a)MSSA法;(b)小波变换法;(c)f-x反褶积法;(d)本文方法
法的残差剖面上无明显有效信号的痕迹,说明本文方法在处理实际资料时具有良好的保幅特性。

3 结论
为了有效地压制地震随机噪声,本文提出了基于卷积降噪自编码器的地震数据去噪方法。

该方法利用卷积降噪自编码器的特性从含噪的地震数据中自适应地学习地震信号的特征,从而滤除较强的随机噪声。

由于使用无监督学习,本文的去噪方法不需要标注无噪的地震数据,可以更加灵活应用。

合成数据和实际数据的应用结果表明,与MSSA、小波变换、f-x反褶积方法相比,本文方法去噪效果最好。

然而,本文方法虽保证了去噪效果,但效率较低,如何提高计算效率需进一步研究。

参考文献
[1] 甘利灯,张昕,王峣钧,等.从勘探领域变化看地震储层预测技术现状和发展趋势[J].石油地球物理勘探,
2018,53(1):214-225.
GAN Lideng,ZHANG Xin,WANG Yaojun,et al.
Current status and development trends of seismic re-
servoir prediction viewed from the exploration industry
[J].Oil Geophysical Prospecting,2018,53(1):214-225.[2] Canales L L.Random noise reduction[J].SEG Tech-nical Program Expanded Abstracts,1984,3:525-527.[3] Gulunay N.FXDECON and complex Wiener predictionfilter[C].SEG Technical Program Expanded Ab-
stracts,1986,5:279-281.
[4] 薛雅娟,曹俊兴.聚合经验模态分解和小波变换相结合的地震信号衰减分析[J].石油地球物理勘探,
2016,51(6):1148-1155.
XUE Yajuan,CAO Junxing.Seismic attenuation ana-
lysis using ensemble empirical mode decomposition
and wavelet transform[J].Oil Geophysical Prospec-
ting,2016,51(6):1148-1155.
[5] Mousavi S M,Langston C A.Hybrid seismic denoi-sing using higher-order statistics and improved wavelet
block thresholding[J].Bulletin of the Seismological
Society of America,2016,106(4):1380-1393.
[6] 张华,陈小宏,李红星,等.曲波变换三维地震数据去噪技术[J].石油地球物理勘探,2017,52(2):226-232.
ZHANG Hua,Chen Xiaohong,Li Hongxing,et al.3D
seismic data de-noising approach based on Curvelet
transform[J].Oil Geophysical Prospecting,2017,
52(2):226-232.
[7] Neelamani R,Baumstein A I,Gillard D G,et al.Co-
 第55卷 第6期宋辉,等:基于卷积降噪自编码器的地震数据去噪1219
herent and random noise attenuation using the curve-
let transform[J].The Leading Edge,2008,27(2):
240-248.
[8] 程浩,王德利,王恩德,等.尺度自适应三维Shearlet变换地震随机噪声压制[J].石油地球物理勘探,
2019,54(5):970-978.
CHENG Hao,WANG Deli,WANG Ende,et al.Seis-
mic random noise suppression based on scale-adaptive
3D-Shearlet transform[J].Oil Geophysical Prospec-
ting,2019,54(5):970-978.
[9] 童思友,高航,刘锐,等.基于Shearlet变换的自适应地震资料随机噪声压制[J].石油地球物理勘探,
2019,54(4):744-750.
TONG Siyou,GAO Hang,LIU Rui,et al.Seismic
random noise adaptive suppression based on the
Shearlet transform[J].Oil Geophysical Prospecting,
2019,54(4):744-750.
[10] 乐友喜,杨涛,曾贤德.CEEMD与KSVD字典训练相结合的去噪方法[J].石油地球物理勘探,2019,
54(4):729-736.
YUE Youxi,YANG Tao,ZENG Xiande.Seismic de-
noising with CEEMD and KSVD dictionary combined
training[J].Oil Geophysical Prospecting,2019,
54(4):729-736.
[11] 邵婕,孙成禹,唐杰,等.基于字典训练的小波域稀疏表示微地震去噪方法[J].石油地球物理勘探,2016,
51(2):254-260.
SHAO Jie,SUN Chengyu,TANG Jie,et al.Micro-
seismic data denoising based on sparse representations
over learned dictionary in the wavelet domain[J].Oil
Geophysical Prospecting,2016,51(2):254-260.
[12] Mi Y.Median filtering in Kirchhoff migration for noisydata[C].SEG Technical Program Expanded Ab-
stracts,2000,19:822-825.
[13] Kavukcuoglu K,Sermanet P,Boureau Y L,et al.Lear-ning convolutional feature hierarchies for visual re-
cognition[C].Advances in Neural Information Pro-
cessing Systems,2010,1090-1098.
[14] Noda K,Yamaguchi Y,Nakadai K,et al.Audio-visualspeech recognition using deep learning[J].Applied
Intelligence,2015,42(4):722-737.
[15] Maggiori E,Tarabalka Y,Charpiat G,et al.Convolu-tional neural networks for large-scale remote-sensing
image classification[J].IEEE Transactions on Geo-
science and Remote Sensing,2016,55(2):645-657.[16] 宋辉,陈伟,李谋杰,等.基于卷积门控循环单元网络的储层参数预测方法[J].油气地质与采收率,2019,
26(5):73-78.
SONG Hui,CHEN Wei,LI Moujie,et al.A method to
predict reservoir parameters based on convolutional
neural network-gated recurrent unit(CNN-GRU)[J].
Petroleum Geology and Recovery Efficieney,2019,
26(5):73-78.
[17] 杨柳青,陈伟,查蓓.利用卷积神经网络对储层孔隙度的预测研究与应用[J].地球物理学进展,2019,
34(4):1548-1555.
YANG Liuqing,CHEN Wei,ZHA Bei.Prediction and
application of reservoir porosity by convolutional neu-
ral network[J].Progress in Geophysics,2019,34(4):
1548-1555.[18] 林年添,付超,张栋,等.无监督与监督学习下的含油气储层预测[J].石油物探,2018,57(4):601-610.
LIN Niantian,FU Chao,ZHANG Dong,et al.Super-
vised learning and unsupervised learning for hydrocar-
bon prediction using multiwave seismic data[J].Geo-
physical Prospecting for Petroleum,2018,57(4):601-
610.
[19] 韩卫雪,周亚同,池越.基于深度学习卷积神经网络的地震数据随机噪声去除[J].石油物探,2018,57(6):
862-869.
HAN Weixue,ZHOU Yatong,CHI Yue.Deep lear-
ning convolutional neural networks for random noise
attenuation in seismic data[J].Geophysical Prospec-
ting for Petroleum,2018,57(6):862-869.
[20] Wang E,Nealon J.Applying machine learning to 3Dseismic image denoising and enhancement[J].Inter-
pretation,2019,7(3):SE131-SE139.
[21] Liu D,Wang W,Chen W,et al.Random noise sup-pression in seismic data:What can deep learning do?
[C].SEG Technical Program Expanded Abstracts,
2018,37:2016-2020.
[22] Liu D,Wang W,Wang X,et al.Poststack seismic datadenoising based on 3-D convolutional neural network
[J].IEEE Transactions on Geoscience and Remote
Sensing,2019,58(3):1598-1629.
[23] Kim Y,Hardisty R,Marfurt K J.Seismic random noiseattenuation in f-x domain using complex-valued resi-
dual convolutional neural network[C].SEG Technical
Program Expanded Abstracts,2019,38:2579-2583.[24] Zhu W,Mousavi S M,Beroza G C.Seismic signal de-noising and decomposition using deep neural networks
[J].IEEE Transactions on Geoscience and Remote
Sensing,2019,57(11):9476-9488.
[25] Zhao Y,Li Y,Dong X,et al.Low-frequency noise sup-pression method based on improved DnCNN in desert
seismic data[J].IEEE Geoscience and Remote Sen-
sing Letters,2018,16(5):811-815.
[26] Oropeza V,Sacchi M.Simultaneous seismic data de-noising and reconstruction via multichannel singular
spectrum analysis[J].Geophysics,2011,76(3):V25-
V32.
[27] 张华,陈小宏,杨海燕.地震信号去噪的最优小波基选取方法[J].石油地球物理勘探,2011,46(1):70-75.
ZHANG Hua,CHEN Xiaohong,YANG Haiyan.Op-
timistic wavelet basis selection in seismic signal noise
elimination[J].Oil Geophysical Prospecting,2011,
46(1):70-75.
(本文编辑:宜明理)
作者简

宋辉 硕士研究生,1996年生;
2018年获长江大学地球科学学院石油
工程专业学士学位;现在该校地球物
理与石油资源学院攻读地质工程专业
硕士学位,主要研究方向为人工智能地
震数据处理。

 Ⅱ Oil Geophysical Prospecting2020 
cannot take into account picking effect and efficien-cy,the algorithm stability is poor,and the industri-al application has not been very mature.The firstarrival picking method based on deep learning istime-consuming and labor-intensive,the process ofdata preprocessing is cumbersome,and the networkstructure is too complex,resulting in low trainingand test efficiency.Combining the advantages of U-Net with those of SegNet,a new hybrid networkU-SegNet is constructed,and based on which firstarrivals can be picked automatically.Based on theSegNet structure,U-SegNet provides multi-scaleinformation of the encoder network by fusing jumpconnections information before the deconvolutionlayer of the decoder network to obtain better per-formance,and its upsampling operation changes thedeconvolution in U-Net to unpooling.Because thepooling index is passed to the upsampling layer,thenetwork model converges faster.Therefore,the U-SegNet network structure is more conducive tosegmenting the background noise area and the areawhere background noise and valid signal overlap,thereby improving the accuracy of first arrivalpicking.The first arrival automatic picking processbased on U-SegNet includes making a training dataset,designing a network model,training the net-work model,testing the network model and apply-ing it to real seismic data.Tests and applications ofthe U-SegNet model show that the picking efficien-cy of the proposed method is about 2.2times thatof a commercial software.It is easy to be industri-alized and has a good future in large-scale application.Keywords:seismic first arrival,pick up,deep learn-ing,U-Net,SegNet,U-SegNet
1.Northwest Branch,Research Institute of Petrole-um Exploration &Development,PetroChina,Lanzhou,Gansu 730020,China
Suppressing seismic random noise based on Deep-KSVD.TANG Jie1,MENG Tao1,ZHANG Wenzheng1,and CHEN Xueguo2.Oil Geophysical Prospecting,2020,55(6):1202-1209.
Curvelet transform denoising causes the e-vents to be distorted and interferes with the effec-tive signals in discontinuous areas such as faultzones.The algorithm of overcomplete dictionariesfor sparse representation(K-SVD)requires manu-al and repeated adjustment of parameters to im-prove the denoising effect.After comprehensivelyconsidering the advantages of deep learning net-work and sparse representation,we combined theK-SVD denoising algorithm with the deep learningnetwork,and proposed the random noise suppres-sion method based on Deep-KSVD.In order to
make the network have the ability to learn parame-ters,the OMP algorithm is replaced by an equiva-lent learnable alternative in the tracking phase.Thecalculation process includes decomposing seismicdata into overlapped data blocks,de-noising eachdata block by proper tracking,and reconstructingthe whole data by weighting the denoised datablocks.The denoising process includes three parts:sparse coding,estimation of regularized coefficientand reconstruction of data block.The test resultson model data and actual data show that aftertraining a Deep-KSVD network,for given noisy da-ta,it can adaptively attenuate the seismic noiseswithout further adjusting parameters while protec-ting the effective information of discontinuity andthe characteristics of data structure.Comparedwith the K-SVD denoising method,the Deep-KS-VD denoising method provides better effect ofnoise suppression and can improve the signal-to-noise ratio of full-band data.
Keywords:Deep-KSVD,random noise,suppression,deep learning,sparse signal
1.School of Geosciences,China University of Pe-troleum(East China),Qingdao,Shandong 266580,China
2.Research Institute of Exploration and Develop-ment,SINOPEC Shengli Oilfied Company,Dongy-ing,Shangdong 257105,China
Seismic noise suppression based on convolutional de-noising autoencoders.SONG Hui 1,2,GAO Yang3,CHEN Wei 1,4,ZHANG Xiang1,2.Oil GeophysicalProspecting,2020,55(6):1210-1219.
Noise attenuation is a long-standing problemin seismic exploration.Traditional denoising me-thods can suppress seismic noise,but they may causelost effective signals,residual noises and otherproblems.An unsupervised denoising algorithmbased on a convolutional denoising autoencoder isproposed,which can significantly improve the sig-nal-to-noise ratio of seismic data.The algorithm lo-cally and randomly damages seismic data,and thentransmits the seismic data damaged to coding anddecoding frameworks.The coding framework cap-tures the waveform of the seismic data and elimi-nates the noise.The decoding framework expandsthe feature map,recovers the details of the seismicdata.Finally,after reconstructing the seismic data,the algorithm trains a model with the error be-tween the reconstructed seismic data and the origi-nal seismic data.Considering the complexity ofseismic data,a multi-scale convolution module isrequired to extract the characteristics of seismic da-ta during coding and decoding.Applications to syn-
 Vol.55 No.6AbstractsⅢ
thetic and real seismic data have proved that thenew method is more effective in preserving signalswhile suppressing noises.Its denoising result isbetter than a traditional algorithm.
Keywords:unsupervised learning,convolutionalneural network,denoising autoencoders,seismicdata,denoising
1.Key Laboratory of Exploration Technologies forOil and Gas Resources(Yangtze University),Min-istry of Education,Wuhan,Hubei 430100,China;2.College of Geophysics and Petroleum Resources,Yangtze University,Wuhan,Hubei 430100,China;3.CNPC Key Laboratory of Geophysical Prospec-ting,China University of Petroleum(Beijing),Bei-jing 102249,China;
4.Hubei Cooperative Innovation Center of Uncon-ventional Oil and Gas,Wuhan,Hubei 430100,ChinaThe shot infilling method based on seismic legacy da-ta.XU Yinpo1,2,SONG Qianggong3,PAN Yingjie1,NI Yudong1,ZOU Xuefeng4,and YU Jianpeng1.OilGeophysical Prospecting,2020,55(6):1220-1230.Based on ray tracing and wave equation illumi-nation result,the conventional method for impro-ving the imaging quality of the shadow zone of atarget layer with complex structures is to analyzethe distribution of irradiating energy in the shadowzone and select the infilled shots according to thecontribution from illuminated energy to the shadowzone.However,it is impossible to accurately createageological model due to the complexity of geolog-ical structures,which will affect the illumination a-nalysis to complex target layer and lead to it is dif-ficult to improve the imaging quality of the targetarea with complex structure by infilling shots.Based on this situation,this paper proposes themethod for improving the imaging quality of shad-ow zone of the target layer on the basis of the lega-cy data according to the local similarity theory:Firstly,the geometry should be drawn up accordingto the geological task of exploration area,and theshot interval(or shot line interval)should be assmall as possible;Secondly,calculate the local simi-larity of each trace of the common imaging point(CIP)gathers of the legacy seismic data and picktarget layer in stacked data;so the local similarityof target layer for all traces of all shots is ob-tained;Thirdly,use spatial interpolation to calcu-late the local similarity of all shots and the contri-bution of the normal shots to the target layer withthe designed geometry,and establish the relationcurve(curved face)between the local similarity ofeach CIP gather and spatial position to determinethe zones with low local similarity value;Finally,use the uniformity of local similarity of target layeras index to determine the positions and density ofinfilled shots.The verifications by theoretical dataand practical data indicate that this method can e-conomically infill suitable shots and improve theimaging effect of shadow zones of the target layerin areas with high steep dip angle and overthruststructures and whose lateral velocities change se-verely.
Keywords:geometry,infilled shots,common ima-ging point(CIP)gathers,local similarity,stackeddata,uniformity
1.Acquisition Technique Center,BGP,CNPC,Zhuozhou,Hebei 072751,China
2.Southwest Petroleum University,Sichuan Prov-ince Natural Gas Geology Key Laboratory,Cheng-du,Sichuan 610500,China
3.Research &Development Center,BGP,CNPC,Zhuozhou,Hebei 072751,China
4.Changqing Division,BGP,CNPC,Xi’an,Shaanxi710021,China
Depth-domain wavelet estimation using the subspace-constrained Huber norm.ZHANG Jie1,2,CHEN Xue-hua1,2,JIANG Wei 1,2,DAN Zhiwei 3,and XIAOWei 3.Oil Geophysical Prospecting,2020,55(6):1231-1236,1244.
Generally,in the constant-velocity depth do-main,only a few hundred meters of logging infor-mation is available,and is equivalent to 2~5timesthe length of a constant depth-domain seismicwavelet.It is difficult to estimate a reliable seismicwavelet from such a short data segment.To ad-dress this issue,we proposed a method for estima-ting the constant-velocity depth-domain wavelet u-sing the subspace-constrained Huber norm.Thecorresponding procedure can be divided into threesteps.First,the reflectivity and the seismic tracenear the well location are transformed from the realdepth domain to the constant-velocity depth do-main.Next,set a threshold and generate the initialsynthetic seismogram by convoluting the constant-velocity depth-domain reflectivity and the initialseismic wavelet.Finally,the seismic wavelet is up-dated by iterative least square method according tothe residual error of the synthetic seismogram andthe seismic trace until the iteration terminationcondition is reached.We compare the proposedmethod and the conventional least square-basedmethods through the synthetic and field seismic da-ta.The results show that the high performance ofour method for estimating the reliable wavelet fromlimited seismic and well-log data segment.The pro-posed method can be further extended to estiamte。

相关文档
最新文档