２０２０年１２
月第５５卷　第６期　
＊山东省青岛市长江西路６６号中国石油大学（华东）地球科学与技术学院，２６６５８０。Ｅｍａｉｌ：ｔａｎｇｊｉｅ＠ｕｐ
ｃ．ｅｄｕ．ｃｎ本文于２０２０年４月２７日收到，
最终修改稿于同年９月７日收到。本项研究受国家自然科学基金项目“基于微地震数据的致密油气储层裂纹演化分形特征研究”（４１５０４０９７
）及“深度偏移地震数据特征剖析与深度域直接反演方法研究”（４１８７４１５３
）联合资助。·智能地球物理·
文章编号：１０００－７２１０（２０２０）０６－１２０２－
０８利用基于深度学习的过完备字典信号
稀疏表示算法压制地震随机噪声
唐　杰＊①　孟　涛①　张文征①　陈学国②
（①中国石油大学（
华东）地球科学与技术学院，山东青岛２６６５８０；②中国石化胜利油田分公司勘探开发研究院，山东东营２５７０１５）唐杰，
孟涛，张文征，陈学国．利用基于深度学习的过完备字典信号稀疏表示算法压制地震随机噪声．石油地球物理勘探，２０２０，５５（６）：１２０２－
１２０９．摘要　曲波变换去噪处理使同相轴在断层等不连续区域发生畸变，对有效信号产生干扰。基于过完备字典信号稀疏表示（Ｋ－ＳＶＤ）需要人工反复调整参数才能改善去噪效果。为此，将Ｋ－ＳＶＤ去噪算法与深度学习网络相结合，综合考虑深度学习网络与稀疏表示方法的优点，研究了基于深度学习的过完备字典信号稀疏表示（Ｄｅｅｐ
－ＫＳＶＤ）的地震数据随机噪声压制方法。为了使该网络有能力学习参数，在追踪阶段用一个等价的
可学习的替代方案代替正交匹配追踪（ＯＭＰ）算法。计算过程包括将数据分解为重叠的数据块、通过适当的追踪对每个数据块去噪以及通过去噪后的数据块加权重建整个数据，去噪处理包括稀疏编码、正则化系数估计以及数据块重建三个部分。模型数据和实际数据测试结果表明：当Ｄｅｅｐ－ＫＳＶＤ网络训练完成后，给定含噪数据，能够自适应地衰减地震噪声，并保护有效不连续性信息及数据结构特点，无需再进行参数调整；与Ｋ－ＳＶＤ去噪方法相比，Ｄｅｅｐ－ＫＳＶＤ去噪方法的噪声压制效果更好，可提高全频带数据的信噪比。关键词　Ｄｅｅｐ
－ＫＳＶＤ　随机噪声　压制　深度学习　稀疏表示中图分类号：Ｐ６３１  文献标识码：Ａ  ｄｏｉ：１０．１３８１０／ｊ
．ｃｎｋｉ．ｉｓｓｎ．１０００－７２１０．２０２０．０６．００４０　引言
提高地震资料的信噪比是地震数据处理的重要任务。随着勘探、开发的精度提高及面对的勘探目标越来越复杂，对地震数据去噪技术的要求也逐步
提高［
１］
。传统地震数据去噪方法通过寻和识别地震数据中的冗余信息，利用地震数据自身的数据结
构特点去噪［２－
３］。随着稀疏表示理论的发展，基于稀
疏表示的去噪方法逐渐进入人们的视野。曹静杰
等［４］
提出一种基于曲波变换的稀疏变换自适应去噪
方法，只需设置初始阈值和迭代次数，便可通过迭代过程中解的稀疏性与拟合误差之间的内在关系确定
合适的阈值，进而压制随机噪声。张华等［５］将非均
匀快速傅里叶变换引入多尺度、多方向二维曲波变换，应用二维非均匀曲波变换压制地震随机噪声，在
将非均匀采样地震数据内插为均匀采样数据的同
时，有效地压制了噪声。与传统方法相比，稀疏表示方法将去噪问题转化为优化问题，基于稀疏变换的方法使用专门设计的变换稀疏表示地震数据。为了提升字典表示的准确性，可以采用学习训练的方法构造表示字典，与其他传统的字典构造方法（如曲波
变换以及Ｓｈｅａｒｌｅｔ变换等）［６－
８］相比，可准确地稀疏表示地震数据。李勇等［９］
提出了一种数据驱动与模
型驱动联合的模型约束下的在线字典学习去噪方法，先通过模型驱动方式获得一个较优质的学习样本以构建字典再进行去噪处理，能在高噪声背景下
有效地提取地震数据中的弱信号。张良等［１０］
等使
用一种基于双稀疏字典和快速迭代收缩阈值算法（ＦＩＳＴＡ）
的方法对地震数据去噪，能够自适应地稀疏表示地震数据，当地震数据较复杂时的信噪比更高，不仅去噪处理速度较快，而且克服了字典学习缺
　第５５卷　第６期唐杰，等：利用基于深度学习的过完备字典信号稀疏表示算法压制地震随机噪声１２０３
少先验约束的不足。近年来基于过完备字典信号稀疏表示（Ｋ－ＳＶＤ）的一系列算法被广泛使用和扩展，为了获得良好的降噪效果，这些字典学习方法需要根据经验反复测试并微调参数［１１］。
随着计算机硬件技术的飞速发展，图形处理能力得到极大提升，深度学习网络自２０１０年以来成为热门话题，深度置信网络［１２］、堆叠式自动编码器［１３］和深度卷积神经网络（ＣＮＮ）［１４］等大量的深度学习网络应运而生。传统去噪方法通过寻和识别地震数据中的冗余信息，利用地震数据自身的数据结构特点去噪，主要集中在数据建模和优化。深度学习网络走了一条完全不同的道路，直接针对推理阶段，学习其参数以优化端到端性能。其中ＣＮＮ已在解决不同领域的高度非线性计算机视觉问题方面取得了成功［１５］。Ｚｈａｎｇ等［１６］使用具有１７个卷积层的去噪卷积神经网络（ＤｎＣＮＮ，ＤｅｎｏｉｓｉｎｇＣｏｎｖｏｌｕｔｉｏｎａｌ　Ｎｅｕｒａｌ　Ｎｅｔｗｏｒｋｓ）将噪声作为输出（即残差学习方法）对数据降噪，降噪效果较好且提高了训练速度。与具有单层分解的字典学习相比，具有更深层次的深度学习网络能够在具有不同抽象级别的地震数据训练集中获得更好的处理效果［１７－１９］。由于单独使用稀疏字典和深度学习方法都存在局限性，因此人们探索将二者结合的去噪方法。将稀疏表示结构嵌入深度学习网络体系中，以使整个算法既有深度方法的灵活性，又有稀疏表示方法带来的信息。
Ｓｃｅｔｂｏｎ等［２０］提出了基于深度学习的过完备字典信号稀疏表示算法（Ｄｅｅｐ－ＫＳＶＤ）用于图像数据去噪，取得了理想的去噪效果。以往在地震数据去噪过程中使用的仅仅是Ｋ－ＳＶＤ
去噪算法。本文进一步研究了基于Ｄｅｅｐ－ＫＳＶＤ的地震数据随机噪声压制方法，将Ｋ－ＳＶＤ去噪算法与深度学习网络相结合，综合考虑深度学习网络与稀疏表示方法的优点，使Ｄｅｅｐ－ＫＳＶＤ算法更有效。
１　算法原理
１．１　Ｋ－ＳＶＤ算法
Ｋ－ＳＶＤ算法由Ａｈａｒｏｎ等［２１］、Ｅｌａｄ等［２２］根据误差最小原则提出，利用基向量的线性组合表示输入信号的数据特征，通过数据训练不断更新字典基向量及其表示系数，将逼近误差调整到最小，完成对
数据的充分表示。Ｋ－ＳＶＤ算法包括稀疏编码和字典更新两个基本步骤。
１．１．１　稀疏编码
将地震数据ｘ表示为一个大小为槡ｐ　×槡ｐ的小数据块，按字典顺序排序成长度为ｐ的列向量（也称为原子）。根据稀疏表示模型，假设ｘ是从预先指定的字典Ｄ∈Ｒｐ×ｍ（ｍ为字典中的原子数量）中获取的ｓ ｐ个列向量的线性组合，即ｘ＝Ｄα，其中α∈Ｒｍ是一个包含ｓ个非零元素的稀疏向量，用‖α‖０＝ｓ表示。将ｘ加入加性零均值高斯白噪声获得含噪结果ｙ，其噪声标准差为σ。噪声估计值的稀
疏编码＾α通过求解下式获得［２３－２４］
  ＾α＝ａｒｇ　ｍｉｎ
α
λ‖α‖０＋
１
２
‖Ｄα－ｙ‖
（）２２（１）
  随后通过＾ｘ＝Ｄ＾α得到去噪结果＾ｘ。式（１）中λ为正则化系数。
假设基础字典Ｄ已知，将上述稀疏编码操作应用于更大的数据及其含噪数据Ｙ时会有两个未知信息：每个位置的稀疏表示和输出数据Ｘ。块坐标最小化算法首先将含噪数据Ｙ作为Ｘ的初始值，然后寻所有位置ｋ的最优αｋ。这个阶段类似对尺寸为槡ｐ×槡ｐ的每个块做一个滑动窗口稀疏编码。１．１．２　
字典更新
假定所有块的稀疏表示为｛＾αｋ｝ｋ，为了更新Ｘ，需求解下式
＾Ｘ＝ａｒｇ　ｍｉｎ
Ｘ
（μ
２
‖Ｘ－Ｙ‖２２＋
１
２
∑
ｋ
‖Ｄ＾αｋ－ＲｋＸ‖２２）（２）式中：μ为正则化参数；Ｒｋ∈Ｒｐ×Ｎ表示从数据中提取的第ｋ块算子。可将Ｄ的计算嵌入贝叶斯公式中
ｍｉｎ
｛α
ｋ
｝
ｋ
，Ｘ，Ｄ［
μ
２
‖Ｘ－Ｙ‖２２＋∑
ｋ
（λｋ‖αｋ‖０＋
１
２
‖Ｄαｋ－ＲｋＸ‖２２）］（３）  在这种情况下，Ｄ使用所有现有的噪声数据块从Ｙ本身学习。采用块坐标最小化将Ｄ初始化为过完备的离散余弦变换（ＤＣＴ）矩阵，并设置Ｘ＝Ｙ，然后在数据块的正交匹配追踪（ＯＭＰ）计算结果和Ｋ－ＳＶＤ更新的Ｄ之间迭代。在循环之后，Ｄ能够适应被处理的数据。
　１２０４　石油地球物理勘探２０２０年　
在实际应用时，不断循环稀疏编码和字典更新这两个步骤，直至达到设定逼近误差标准或迭代次数。基于构造学习字典的稀疏表示方法通过学习训练算法从输入数据中提取数据特征，完成对字典基函数的不断更新和完善。因此，相较于传统方法，该方法具有更好的数据适应性，可提高对不同数据结构特点的表示精度，从而改善去噪效果。１．２　Ｄｅｅｐ
－ＫＳＶＤ去噪原理Ｄｅｅｐ
－ＫＳＶＤ方法结合了Ｋ－ＳＶＤ稀疏去噪算法和深度学习的思想［２０］
，为了使该网络有能力学习参数，在追踪阶段用一个等价的可学习的替代方案代替ＯＭＰ算法。计算过程包括将数据分解为重叠的数据块、通过适当的追踪对每个数据块去噪以及
通过去噪后的数据块加权重建整个数据［
２５－
２６］。图１为Ｄｅｅｐ－ＫＳＶＤ网络结构，其去噪处理包括稀疏编码、λ估计以及数据块重建三个部分，
分述如下
。图１　Ｄｅｅｐ
－ＫＳＶＤ网络结构１．２．
１　稀疏编码给定一个加入标准差为σ的高斯噪声的数据块
ｙ∈Ｒ槡
ｐ×槡ｐ（保持长度为ｐ的列向量），目标是根据已知字典Ｄ∈Ｒｐ×ｍ获得其稀疏编码［２７－２８］。该问题的近似解为
  ＾α＝ａｒｇ　ｍｉｎα
１２‖Ｄα－ｙ‖２２＋λ‖α‖（）
１（４
）通过迭代软阈值算法（ＩＳＴＡ）求解上式，能够保证收敛到全局最优
＾αｔ＋１＝Ｓλ
／ｃ＾αｔ－１ｃＤＴ（Ｄ＾αｔ－ｙ［］
）＾α０＝烅
烄烆０
（５
）式中：ｃ为Ｄ的平方谱范数；Ｓλ
／ｃ为分量软阈值算子；＾αｔ为第
ｔ次迭代时求得的稀疏编码。在求解过程中采用梯度下降法迭代求取最优解，
这样可以使稀疏编码阶段具有可学习性［
２９－
３０］。１．２．２　λ估计
正则化系数λ不仅依赖于σ，而且还依赖于数据块ｙ本身。根据Ｋ－ＳＶＤ去噪算法，
需要为每个数据块ｙｋ
设置λｋ，
以产生具有控制误差水平的稀疏表示‖Ｄ＾αｋ－ｙｋ‖２２≤ｐσ２。可从ｙｋ出发回归相应的正则化参数λｋ，
使用多层感知器（ＭＬＰ）网络表示为λ＝ｆθ（ｙ），其中θ是ＭＬＰ参数的向量。ＭＬＰ由三个隐藏层组成，每层都由一个全连接的线性映射和一个激活函数（ＲｅＬＵ）（除了最后一层）组成。输入层由ｐ个节点组成，ｐ为矢量化数据块的维数，输出层由单个节点组成。１．２．３　数据块重建
使用Ｄ和稀疏编码＾α重建数据块ｙ对应的无
噪组合＾ｘ。学习网络中，字典代表一组参数。Ｄｅｅｐ－ＫＳＶＤ首先将输入数据分割成完全重叠的数据块，然后通过上述的数据块去噪方案处理每个含噪数据块，最后通过加权平均这些数据块的去噪结果重建
数据［
３１－
３２］。最后一个阶段允许学习加权组合的数据块，用ｗ∈Ｒ槡ｐ×槡ｐ
表示该数据块权值矩阵，
得到重建结果
＾Ｘ＝
∑ｋ
ＲＴ
ｋ（ｗ⊙＾
ｘｋ）∑ｋ
ＲＴｋ
ｗ
（６
）式中⊙表示Ｓｃｈｕｒ乘积。Ｄｅｅｐ
－ＫＳＶＤ网络Ｆ是θ、ｃ、Ｄ和ｗ的参数化函数。给定含噪数据Ｙ，
通过＾Ｘ＝Ｆ（Ｙ）返回数据去噪结果。由上述所有参数训练Ｆ，将损失函数Ｌ＝
∑ｉ
‖Ｘｉ
－Ｆ（
Ｙｉ
）‖２
２最小化，
集合｛Ｘｉ｝ｉ代表训练数据，｛Ｙｉ｝ｉ为合成噪声形式，
由Ｙｉ＝Ｘｉ＋Ｖｉ获得，其中Ｖｉ为零均值高斯白噪声。
１．３　网络训练本文基于Ｐｙ
ｔｏｒｃｈ构建网络模型，使用一台ＣＰＵ为ｉ７－
７７００、内存为２０Ｇ、系统为Ｗｉｎｄｏｗｓ１０专业版６４位操作系统的台式电脑训练网络模型。利用２７７０个大小为１２８×１２８的地震数据建立数据集，其中２５００个数据组成训练集，２７０个数据组成测试集，训练和测试集是严格不相交的，训练阶段随机抽取训练集中的数据，添加高斯白噪声。为了更合理地选择网络参数，为网络选取Ｌ１、Ｌ２和ＳｍｏｏｔｈＬ１等损失函数进行训练并测试（图２）。可见，当选取Ｓｍｏｏｔｈ　Ｌ１损失函数时，训练（图２ａ）和测试（图２ｂ）的函数值收敛更好，即：迭代次数为０～３００时函数值逐步降至约０．００５；迭代次数为３００～
　第５５卷　第６期唐杰，等：利用基于深度学习的过完备字典信号稀疏表示算法压制地震随机噪声
１
２０５  ２５００时损失函数值趋于稳定。初步说明网络模型得到有效训练
。
图２　训练（ａ）与测试（ｂ
）损失函数曲线在训练集和测试集中添加均值为０、噪声水平为数据幅值２５％的高
斯白噪声训练和测试网络。使用适应性矩估计（Ａｄａｍ）
优化函数，将学习率设置为１×１０
－４，使用过完备ＤＣＴ矩阵作为字典Ｄ，使用Ｄ的平方谱范数初始化稀疏编码参数ｃ，
其他参数随机初始化２　模型数据测试
为了验证网络模型的训练效果，使用添加高斯随机噪声的不同信噪比（ＳＮＲ）的含噪模型数据测试网络的去噪能力（图３），结果表明，不同信噪比模型数据经过Ｄｅｅｐ－ＫＳＶＤ算法去噪后显著提高了信噪比，进一步说明网络模型的训练效果很好。
为了更形象地展示对模型数据的去噪效果，选取并展示ＳＮＲ＝８的含噪模型数据及其去噪结果（图４）。可见：对原始含噪数据（图４ａ）的曲波变换去噪结果（图４ｂ）残留噪声点，在残差剖面（图４ｅ）边缘存在有效信号；Ｋ－ＳＶＤ去噪结果（图４ｃ）残留噪声，并存在轻微的伪影现象，部分地层信息较模糊，在残差剖面（图４ｆ）上存在有效信号；Ｄｅｅｐ
－ＫＳＶＤ
去图３　ＳＮＲ＝３、４、５、６、７、８的模型数据去噪前、
后信噪比图４　模型数据的去噪结果
（ａ）原始含噪数据（ＳＮＲ＝８）；（ｂ）曲波变换去噪结果；（ｃ）Ｋ－ＳＶＤ去噪结果；（ｄ）Ｄｅｅｐ－ＫＳＶＤ去噪结果；（ｅ）图ｂ与图ａ数据之差；（ｆ）图ｃ与图ａ数据之差；（ｇ
）图ｄ与图ａ数据之差
　１２０６　石油地球物理勘探２０２０年　
噪结果（图４ｄ）及其残差剖面（图４ｇ）的断点及地层尖灭点处的同相轴没有畸变，弯曲程度较大的同相轴也无变形，说明Ｄｅｅｐ－ＫＳＶＤ方法可较好地保护断层和同相轴边缘等不连续性信息。
正则化可以产生稀疏权值３　实际数据测试
３．１　测试数据１
为了进一步检验实际数据应用效果，对Ｎ区实
际地震资料采用不同方法去噪并分析去噪效果（图５）。可见：①叠后地震剖面中存在较严重的随机干扰，同相轴模糊且同相轴边界及断点不明显、地层接触关系不清，给处理、解释带来困难（图５ａ）。②曲波变换去噪结果（图５ｂ）存在严重伪影现象，且存在多个噪声点。③与Ｋ－ＳＶＤ去噪结果（图５ｃ）相比，Ｄｅｅｐ－ＫＳＶＤ去噪结果（图５ｄ
）削弱了背景斑块的影响，由两者的残差剖面（图５ｆ、图５ｇ）可知，Ｄｅｅｐ－ＫＳＶＤ方法能够更好地保留有效信号，去噪效果更好
。
图５　不同方法的去噪结果
（ａ）叠后地震剖面；（ｂ）曲波变换去噪结果；（ｃ）Ｋ－ＳＶＤ去噪结果；（ｄ）Ｄｅｅｐ
－ＫＳＶＤ去噪结果；（ｅ）图ｂ与图ａ数据之差；（ｆ）图ｃ与图ａ数据之差；（ｇ
）图ｄ与图ａ数据之差图ａ共有３３１道数据，采样间隔为０．
００２ｓ，采样时长为０．８ｓ  图６为Ｄｅｅｐ
－ＫＳＶＤ去噪前、后信噪比谱。由图可见，去噪后中、低频范围的信噪比显著提升，进一步量化说明了Ｄｅｅｐ－ＫＳＶＤ对实际含噪数据也具有较好的去噪效果。
基于Ｄｅｅｐ－ＫＳＶＤ的去噪方法可以有效地提高资料的信噪比，以稀疏表示为基础，目标明确，结
构简洁，效果良好，主要特点为：①由于空间冗余信息存在于数据块，因此Ｄｅｅｐ
－ＫＳＶＤ网络采用稀疏建模方法；②考虑了数据块的平均效应；③Ｄｅｅｐ
－ＫＳ－ＶＤ网络利用非局部自相似性作为额外的先验信
息，因此可有效捕捉相隔很远的数据块之间的相关性，从而改进去噪效果
。
图６　Ｄｅｅｐ
－ＫＳＶＤ去噪前、后信噪比谱计算时采用相关法计算信噪比，
实际计算过程中采用相邻地震记录的最大互相关值为地震信号的
功率谱，最后分频计算信噪比谱［
２］３．２　测试数据２
对叠前数据做切除直达波和面波处理，得到含

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。