唐 杰 孟 涛 张文征 陈学国

(①中国石油大学(华东)地球科学与技术学院，山东青岛 266580；②中国石化胜利油田分公司勘探开发研究院,山东东营 257015)

0 引言

提高地震资料的信噪比是地震数据处理的重要任务。随着勘探、开发的精度提高及面对的勘探目标越来越复杂，对地震数据去噪技术的要求也逐步提高[1]。传统地震数据去噪方法通过寻找和识别地震数据中的冗余信息，利用地震数据自身的数据结构特点去噪[2-3]。随着稀疏表示理论的发展，基于稀疏表示的去噪方法逐渐进入人们的视野。曹静杰等[4]提出一种基于曲波变换的稀疏变换自适应去噪方法，只需设置初始阈值和迭代次数，便可通过迭代过程中解的稀疏性与拟合误差之间的内在关系确定合适的阈值，进而压制随机噪声。张华等[5]将非均匀快速傅里叶变换引入多尺度、多方向二维曲波变换，应用二维非均匀曲波变换压制地震随机噪声，在将非均匀采样地震数据内插为均匀采样数据的同时，有效地压制了噪声。与传统方法相比，稀疏表示方法将去噪问题转化为优化问题，基于稀疏变换的方法使用专门设计的变换稀疏表示地震数据。为了提升字典表示的准确性，可以采用学习训练的方法构造表示字典，与其他传统的字典构造方法(如曲波变换以及Shearlet变换等)[6-8]相比，可准确地稀疏表示地震数据。李勇等[9]提出了一种数据驱动与模型驱动联合的模型约束下的在线字典学习去噪方法，先通过模型驱动方式获得一个较优质的学习样本以构建字典再进行去噪处理，能在高噪声背景下有效地提取地震数据中的弱信号。张良等[10]等使用一种基于双稀疏字典和快速迭代收缩阈值算法(FISTA)的方法对地震数据去噪，能够自适应地稀疏表示地震数据，当地震数据较复杂时的信噪比更高，不仅去噪处理速度较快，而且克服了字典学习缺少先验约束的不足。近年来基于过完备字典信号稀疏表示(K-SVD)的一系列算法被广泛使用和扩展，为了获得良好的降噪效果，这些字典学习方法需要根据经验反复测试并微调参数[11]。

随着计算机硬件技术的飞速发展，图形处理能力得到极大提升，深度学习网络自2010年以来成为热门话题，深度置信网络[12]、堆叠式自动编码器[13]和深度卷积神经网络(CNN)[14]等大量的深度学习网络应运而生。传统去噪方法通过寻找和识别地震数据中的冗余信息，利用地震数据自身的数据结构特点去噪，主要集中在数据建模和优化。深度学习网络走了一条完全不同的道路，直接针对推理阶段，学习其参数以优化端到端性能。其中CNN已在解决不同领域的高度非线性计算机视觉问题方面取得了成功[15]。Zhang等[16]使用具有17个卷积层的去噪卷积神经网络(DnCNN，Denoising Convolutional Neural Networks)将噪声作为输出(即残差学习方法)对数据降噪，降噪效果较好且提高了训练速度。与具有单层分解的字典学习相比，具有更深层次的深度学习网络能够在具有不同抽象级别的地震数据训练集中获得更好的处理效果[17-19]。由于单独使用稀疏字典和深度学习方法都存在局限性，因此人们探索将二者结合的去噪方法。将稀疏表示结构嵌入深度学习网络体系中，以使整个算法既有深度方法的灵活性，又有稀疏表示方法带来的信息。

Scetbon等[20]提出了基于深度学习的过完备字典信号稀疏表示算法(Deep-KSVD)用于图像数据去噪，取得了理想的去噪效果。以往在地震数据去噪过程中使用的仅仅是K-SVD去噪算法。本文进一步研究了基于Deep-KSVD的地震数据随机噪声压制方法，将K-SVD去噪算法与深度学习网络相结合，综合考虑深度学习网络与稀疏表示方法的优点，使Deep-KSVD算法更有效。

1 算法原理

1.1 K-SVD算法

K-SVD算法由Aharon等[21]、Elad等[22]根据误差最小原则提出，利用基向量的线性组合表示输入信号的数据特征，通过数据训练不断更新字典基向量及其表示系数，将逼近误差调整到最小，完成对数据的充分表示。K-SVD算法包括稀疏编码和字典更新两个基本步骤。

1.1.1 稀疏编码

(1)

1.1.2 字典更新

(2)

式中：μ为正则化参数；Rk∈Rp×N表示从数据中提取的第k块算子。可将D的计算嵌入贝叶斯公式中

(3)

在这种情况下，D使用所有现有的噪声数据块从Y本身学习。采用块坐标最小化将D初始化为过完备的离散余弦变换(DCT)矩阵，并设置X=Y，然后在数据块的正交匹配追踪(OMP)计算结果和K-SVD更新的D之间迭代。在循环之后，D能够适应被处理的数据。

在实际应用时，不断循环稀疏编码和字典更新这两个步骤，直至达到设定逼近误差标准或迭代次数。基于构造学习字典的稀疏表示方法通过学习训练算法从输入数据中提取数据特征，完成对字典基函数的不断更新和完善。因此，相较于传统方法，该方法具有更好的数据适应性，可提高对不同数据结构特点的表示精度，从而改善去噪效果。

1.2 Deep-KSVD去噪原理

Deep-KSVD方法结合了K-SVD稀疏去噪算法和深度学习的思想[20]，为了使该网络有能力学习参数，在追踪阶段用一个等价的可学习的替代方案代替OMP算法。计算过程包括将数据分解为重叠的数据块、通过适当的追踪对每个数据块去噪以及通过去噪后的数据块加权重建整个数据[25-26]。图1为Deep-KSVD网络结构，其去噪处理包括稀疏编码、λ估计以及数据块重建三个部分，分述如下。

图1 Deep-KSVD网络结构

1.2.1 稀疏编码

(4)

通过迭代软阈值算法(ISTA)求解上式，能够保证收敛到全局最优

(5)

在求解过程中采用梯度下降法迭代求取最优解，这样可以使稀疏编码阶段具有可学习性[29-30]。

1.2.2λ估计

1.2.3 数据块重建

(6)

1.3 网络训练

本文基于Pytorch构建网络模型，使用一台CPU为i7-7700、内存为20G、系统为Windows10专业版64位操作系统的台式电脑训练网络模型。利用2770个大小为128×128的地震数据建立数据集，其中2500个数据组成训练集，270个数据组成测试集，训练和测试集是严格不相交的，训练阶段随机抽取训练集中的数据，添加高斯白噪声。为了更合理地选择网络参数，为网络选取L1、L2和Smooth L1等损失函数进行训练并测试(图2)。可见，当选取Smooth L1损失函数时，训练(图2a)和测试(图2b)的函数值收敛更好，即： 迭代次数为0～300时函数值逐步降至约0.005； 迭代次数为300～2500时损失函数值趋于稳定。初步说明网络模型得到有效训练。

图2 训练(a)与测试(b)损失函数曲线

在训练集和测试集中添加均值为0、噪声水平为数据幅值25%的高斯白噪声训练和测试网络。使用适应性矩估计(Adam)优化函数，将学习率设置为1×10-4，使用过完备DCT矩阵作为字典D，使用D的平方谱范数初始化稀疏编码参数c，其他参数随机初始化

2 模型数据测试

为了验证网络模型的训练效果，使用添加高斯随机噪声的不同信噪比(SNR)的含噪模型数据测试网络的去噪能力(图3)，结果表明，不同信噪比模型数据经过Deep-KSVD算法去噪后显著提高了信噪比，进一步说明网络模型的训练效果很好。

图3 SNR=3、4、5、6、7、8的模型数据去噪前、后信噪比

为了更形象地展示对模型数据的去噪效果，选取并展示SNR=8的含噪模型数据及其去噪结果(图4)。可见：对原始含噪数据(图4a)的曲波变换去噪结果(图4b)残留噪声点，在残差剖面(图4e)边缘存在有效信号； K-SVD去噪结果(图4c)残留噪声，并存在轻微的伪影现象，部分地层信息较模糊，在残差剖面(图4f)上存在有效信号；Deep-KSVD去噪结果(图4d)及其残差剖面(图4g)的断点及地层尖灭点处的同相轴没有畸变，弯曲程度较大的同相轴也无变形，说明Deep-KSVD方法可较好地保护断层和同相轴边缘等不连续性信息。

图4 模型数据的去噪结果

3 实际数据测试

3.1 测试数据1

为了进一步检验实际数据应用效果，对N区实际地震资料采用不同方法去噪并分析去噪效果(图5)。可见：①叠后地震剖面中存在较严重的随机干扰，同相轴模糊且同相轴边界及断点不明显、地层接触关系不清，给处理、解释带来困难(图5a)。②曲波变换去噪结果(图5b)存在严重伪影现象，且存在多个噪声点。③与K-SVD去噪结果(图5c)相比，Deep-KSVD去噪结果(图5d)削弱了背景斑块的影响，由两者的残差剖面(图5f、图5g)可知，Deep-KSVD方法能够更好地保留有效信号，去噪效果更好。

图5 不同方法的去噪结果

图6为Deep-KSVD去噪前、后信噪比谱。由图可见，去噪后中、低频范围的信噪比显著提升，进一步量化说明了Deep-KSVD对实际含噪数据也具有较好的去噪效果。

基于Deep-KSVD的去噪方法可以有效地提高资料的信噪比，以稀疏表示为基础，目标明确，结构简洁，效果良好，主要特点为： ①由于空间冗余信息存在于数据块，因此Deep-KSVD网络采用稀疏建模方法； ②考虑了数据块的平均效应； ③Deep-KSVD网络利用非局部自相似性作为额外的先验信息，因此可有效捕捉相隔很远的数据块之间的相关性，从而改进去噪效果。

图6 Deep-KSVD去噪前、后信噪比谱

计算时采用相关法计算信噪比，实际计算过程中采用相邻地震记录的最大互相关值为地震信号的功率谱，最后分频计算信噪比谱[2]

3.2 测试数据2

对叠前数据做切除直达波和面波处理，得到含噪数据(图7a)，其中包含较强的随机噪声。采用Deep-KSVD方法压制随机噪声，有效衰减了随机噪声且有效信号未遭到过度破坏(图7b)。因此Deep-KSVD方法可为解释、处理提供信噪比较高的地震数据。

图8为图7数据第200道的F-K谱与时频谱。

图7 叠前实际数据去噪结果

由图可见： ①原始数据中含有较强的随机噪声，导致F-K谱较杂乱(图8a)； 去噪后F-K谱得到改善(图8b)。②由去噪前(图8c)、后(图8d)时频谱可见，后者的高频段能量减少，且充分保留了低频段能量，说明有效压制了随机噪声且没有破坏低频段有效信号。

图8 图7数据第200道的F-K谱与时频谱

4 结束语

曲波变换去噪处理使同相轴在断层等不连续区域发生畸变，对有效信号产生干扰。K-SVD算法需要人工反复调整参数才能改善去噪效果。本文将K-SVD去噪算法与深度学习网络相结合，综合考虑深度学习网络与稀疏表示方法的优点，研究了基于Deep-KSVD的地震数据随机噪声压制方法。模型数据和实际数据测试结果表明，当Deep-KSVD网络训练完成后，给定含噪数据，能够自适应地衰减地震噪声，并保护有效不连续性信息及数据结构特点，无需再进行参数调整。与K-SVD去噪方法相比，Deep-KSVD去噪方法的噪声压制效果更好，可提高全频带数据的信噪比。