在地震资料采集过程中，噪声会降低地震数据的质量，因此压制噪声很关键[1-3]。在老油田的勘探开发过程中，地震资料采集经常受到抽油机噪声的干扰。抽油机噪声的频带很窄，能量主要集中在个别频率[4]，且其频带范围完全与有效信号重叠，导致抽油机噪声的压制十分困难。

为了压制抽油机噪声，王鑫[5]提出了检波器的三维矢量组合，通过频率方向特性压制抽油机噪声。但在实际中寻找最优的检波器排列组合比较困难，并且不能保证检波器按预先设计埋置。赵彦青等[6]根据盲源分离原理和信噪比模型，将有效数据与抽油机噪声分离。李学良等[7]利用抽油机噪声与有效数据之间不同的空间特征，提出了一种基于形态成分分析的抽油机噪声迭代衰减方法。这类方法都是基于一些假设或先验信息，然后根据经验去调整参数，所以在实际数据处理时速度较慢且普适性较差。

随着深度学习技术在计算机视觉和自然语言处理等方面的快速发展，它在地震勘探领域也得到推广应用。目前，深度学习技术主要用于断层自动检测[8]、地震数据重建[9]、岩相识别[10]和噪声衰减[11]等。在地震资料噪声压制方面的应用多集中于随机噪声，而在更为复杂和具有挑战性的抽油机噪声方面的研究仍然较少。

近年来，深度学习去噪方法主要基于卷积神经网络(Convolutional Neural Networks,CNN)[12-16]和生成对抗网络(Generative Adversarial Networks,GAN)[17-21]。韩卫雪等[12]将CNN应用于随机噪声压制。王钰清等[13]在CNN的基础上加入了数据增广。Zhao等[15]从块大小方面对原始去噪卷积神经网络(Denoising Convolutional Neural Network,Dn CNN)进行了修改，使其适用于低频沙漠随机噪声抑制。Yuan等[20]将残差学习引入循环生成对抗网络(Cycle-Consistent Adversarial Networks,CycleGAN)以提高训练效率。这些方法在随机噪声和地滚噪声压制中都取得了较好的效果。

抽油机噪声是一种相干噪声，主要来源于抽油过程中电机的振动，其传播路径是以抽油机为中心沿地面向四周传播。主频为20～23 Hz，衰减速度比地震波更快。抽油机噪声在炮集上表现为强振幅特征。距离越远的检波器接收到的抽油机噪声能量越小，但在同一位置的检波器接收的抽油机噪声能量变化不大。当抽油机与检波器组合方向垂直时，抽油机噪声呈双曲线型分布规律。由于抽油机噪声频率和分布特性与随机噪声、地滚噪声不同，因此针对后两种噪声类型的深度学习方法并不能够有效地压制抽油机噪声。在针对抽油机噪声压制的方法中，Ren等[22]提出了多层生成器网络(Multi-layer Generator Network,MLGNet)，效果较好，但其网络模型计算量较大，训练时间较长，且存在有效数据失真的问题。

为此，考虑到抽油机噪声的频率和局部分布的特性，本文提出了一种基于多尺度窗口生成器网络(Multi-scale Window Generator Network,MSWG-Net)，该网络采用双层结构，每层分别有编码器和解码器，且都使用了残差连接。编码器采用的反瓶颈设计可以提取到有效数据更多的纹理特征；解码器采用更少的卷积层，在保证噪声压制效果的同时减少了模型的参数量。同时采用多尺度窗口，将不同尺度窗口提取的潜在特征进行融合，扩大模型的感受野，将抽油机噪声逐步进行压制。利用模拟数据和实际数据，并与Dn CNN、GAN、MLGNet的去噪结果进行比较，以测试该方法的噪声压制效果。

1、方法原理

1.1 生成器结构

MSWGNet关键之一是生成器的架构，每个生成器均包含一个编码器和一个解码器（图1）。其中，编码器通过层层卷积提取特征，解码器使用这些特征来重构数据。

编码器（图1a）输入的数据为256×256矩阵。采用32个3×3卷积核处理，再经过3个残差块、下采样、18个残差块、下采样、3个残差块，输出64×64矩阵的数据。

每个残差块都设计了反瓶颈结构。通常情况下，数据维度越低，卷积层乘法计算量越小，数据特征提取效果越差。因此，提取更多的数据特征需要增大数据的维度。反瓶颈结构可以在参数量较少的情况下提取数据中更多的信息。反瓶颈结构包含了三层卷积和两层Re LU激活函数，残差块中间层卷积核的数量为两端卷积核数量的4倍。第一层和第三层均使用1×1卷积核，主要是对不同通道的像素点进行特征融合，升高或降低通道数；第二层卷积主要是进行局部特征的提取。

图1 生成器中的编码器（a）和解码器（b）结构

Conv、Deconv分别为卷积、反卷积，1×1、3×3、4×4都为卷积核尺寸；1、32、64、128、256、512都为卷积核数量；s=2表示卷积步长为2。×2、×3、×18表示残差块个数分别为2、3、18;Re LU为激活函数。

解码器（图1b）输入的数据为64×64的矩阵，经过2个残差块、上采样、2个残差块、上采样、2个残差块，最后采用32个3×3的卷积核处理，输出数据为256×256的矩阵。残差块包含了两层卷积和一层Re LU激活函数，其两层卷积的卷积核数量保持不变。解码器的残差块相较于编码器减少了卷积层数，降低了模型的参数量以进一步提高模型的性能。

编码器和解码器中的每个块都使用了残差连接。一方面，在深层网络中引入残差连接，梯度可以更轻松地传播到较浅的层，从而解决梯度消失问题；另一方面，残差连接允许网络直接学习恒等映射，可以更好地利用每个编码器提取的特征，将低层提取到的特征迁移到高层语义信息中。利用编码器和解码器之间残差块数量的不对等，可以在保持尽可能低的计算成本的同时达到较好的噪声压制效果。

1.2 多尺度窗口生成器网络结构

本文提出的多尺度窗口生成器网络结构如图2所示，网络使用不同大小的窗口对地震数据进行特征提取。构建的网络采用了双层（生成器）设计。

图2 多尺度窗口生成器网络结构

第一层在输入数据上划分为四个尺度相同的窗口，窗口内的数据分别放入生成器，使其专注于细节特征的提取，即

式中：Ei、Di分别表示第i层的编码器、解码器；N1,j表示在第1层的第j个窗口噪声数据；F1,j表示编码后的第1层的第j个窗口的特征图；F1表示第1层拼接后的特征图；Ri表示解码后的第i层的数据；©表示小窗口间的拼接。

第二层使用与原始尺度相同的窗口，使其对整个数据的特征建立关联，即

式中：N2表示第2层的噪声数据；F2表示第2层编码后的特征图；“+”表示两个数据逐点的叠加。

将第一层与第二层得到的特征数据相加并输入至第2层的解码器，得到去噪后的结果

可以看出，第二层利用了第一层提取到的潜在特征，完成了对抽油机噪声的压制。

这些从多个图像块中提取的潜在特征扩大了网络的感受野，保留了丰富的地震数据的细节。为了进一步保留细节特征，每层之间没有进行下采样操作，每层输入数据块的总大小保持一致。这样不仅保留了有效信号的细节特征，学习到有效数据的分布规律，同时提高了对抽油机噪声的压制效果。通过第一层解码得到的地震数据块与第二层待输入的地震数据进行逐点叠加，同时第一层编码器提取的特征与第二层编码器提取的特征进行逐点叠加，利用各层潜在特征的可重用性，将第一层提取到的潜在特征向第二层传递，使有效信号极大程度地得到保留,从而压制了抽油机的噪声。

1.3 评价指标

采用三个指标评价本文网络效果，即信噪比(SNR)、均方误差(MSE)和结构相似性(SSIM)。它们可分别定义为

式中：s为纯净地震数据；x为含噪声数据；μs和μx分别代表s和x的平均值；σs、σx分别代表s、x的标准差；σxs代表x和s的协方差；c1、c2为常数，以避免分母为0带来的错误。

2、数据测试

通过合成数据和实际数据分别测试本文方法效果。采用的工作站为linux64位操作系统，版本为Ubuntu 22.4.0,GPU为RTX3060,Py Torch版本为1.8.1,Python版本为3.7。采用Adam优化器和MSE损失函数对网络进行训练。去噪效果评价指标为MSE、SNR、SSIM。抽取单道数据进行分析，进一步明确本文方法的噪声压制效果。

2.1 合成数据实验

抽油机噪声以井为中心向四周呈放射状、双曲线型分布，其干扰半径为0～130 m[4]。因此，在模拟数据时，将井泵震源设置在距检波器130 m范围内不同距离的地面上。为了使合成数据的特征更接近于实际地震数据，正演模型选择Marmousi-II纵波速度模型。模型尺寸为271×351，时间采样间隔为0.5 ms，空间步长H为5 m，总接收时间为3.6 s，主频范围为20～30 Hz。利用Ricker子波，同时采用一阶应力—速度曲线，共获得31炮类双曲线抽油机噪声的地震记录。

将31炮合成地震数据划分为256×256的块，得到的合成训练数据共包括5632对含有和不含有抽油机噪声的数据块。将数据集按照4:1的比例划分为训练集和测试集。

为了对比方法的有效性，将本文方法与用于压制随机噪声的Dn CNN[15]、用于压制地滚噪声的GAN[20]和用于压制抽油机噪声的MLGnet[22]进行对比。

从表1可以看出，Dn CNN、GAN和MLGnet与本文网络(MSWGNet)都对抽油机噪声进行了一定程度地压制，但本文方法效果明显优于前三者。

表1 不同网络的去噪效果

由图3可见，加入抽油机噪声后的合成数据(图3a下)的信噪比为-3.92 d B，抽油机噪声与地震有效数据融为了一体。图3b上显示了Dn CNN去噪后的效果，SNR为4.04 d B，仍然有较为明显的抽油机噪声残留；图3b下为图3b上与图3a上之间的残差数据，可以看出Dn CNN去噪后有效数据的丢失。图3 c上和图3 d上为分别利用GAN和MLGnet对数据进行去噪的结果，信噪比分别为10.76 d B和11.29 d B,GAN和MLGnet的去噪性能优于Dn CNN，但仍然有抽油机噪声残留，并且损失了有效信号。图3e上说明了MSWGnet应用于图3a下中的数据后的结果，抽油机噪声在样本中被明显去除，SNR可以恢复到20.27 d B。图3a上与图3e上之间的残差结果如图3e下所示，有效数据没有明显的丢失。

综上所述，与其他三种方法相比，MSWGnet有着最佳的去噪性能，有效信号损失最小。

由去噪后数据的频率—波数谱(f-k谱)（图4）可见，Dn CNN（图4c）、GAN（图4d）、MLGNet（图4e）并没有完全压制抽油机噪声，同时产生其他频率的噪声；MSWGNet（图4f）对抽油机噪声的压制效果较好，并没有产生其他频率的噪声，去噪后的f-k谱与加噪前的f-k谱最为接近。因此，同样证明了本文方法的优越性。

图3 模拟数据和加抽油机噪声数据不同方法去噪结果（上）及残差（下）

图4 图3模拟数据及加抽油机噪声数据不同方法去噪结果的频率—波数谱

图5表示单道数据去噪结果的差异。可以看出Dn CNN、GAN、MLGNet的去噪结果与加噪前的数据有较大的差异，而MSWGNet去噪后的结果与加噪前的数据基本重合。

2.2 实际数据实验

为了使训练数据与实际数据更加接近，验证该方法对实际数据的有效性，使用实际数据制作的训练集对网络进行训练。在准备实际数据训练集时，为了得到有噪声和无噪声的数据对，提取了实际数据中与有效数据不重叠的抽油机噪声，并将其添加到没有抽油机噪声的地震数据上，产生有噪声数据，没有抽油机噪声的地震数据作为无噪声数据。通过该方法得到256×256的8352对数据块，其中6681对数据块构成训练集，1671对数据块构成测试集。在测试阶段，只需将含噪地震数据输入训练好的网络，预测后即可得到去噪后的地震数据。

图5 不同方法的单道地震数据去噪结果

为了进一步验证本文网络对野外实际数据的有效性，使用带有抽油机噪声的野外实际地震数据输入训练好的网络进行处理，结果如图6所示。由图可见，Dn CNN和GAN对抽油机噪声的压制能力较弱；MLGnet对抽油机噪声进行了一定程度地压制，但噪声仍有残留，且有时会产生一些新的噪声；MSWGnet方法有效地去除了抽油机噪声，最大程度地保留了有效信号且没有产生新的噪声。

为了更清晰地展示本文方法的去噪结果，对图6中红色方框内的数据进行放大（图7）。由图可见，去噪后的结果具有更清晰的地震反射特征，且有效信号保留较好，抽油机噪声压制较为彻底。

图8给出了图6数据的f-k谱。图8e中没有发现相干噪声，这表明了本文方法的有效性。

图6 实际数据及其不同方法去噪结果

图7 图6红框部分放大显示

图8 图6实际数据及其不同方法结果的频率—波数谱

3、结束语

本文针对抽油机噪声的压制提出了一种基于多尺度窗口生成器网络，即通过使用尺度不同的窗口对含有抽油机噪声的地震数据进行特征提取，通过双层生成器与不同尺度的窗口结合可以提取到地震数据中更多信息以进行抽油机噪声压制。模拟数据和野外实际数据测试均表明，本文提出的MSWG-net与Dn CNN、GAN、MLGnet相比可以有效地对抽油机噪声进行压制，并能完整地保留有效信号。同时模型的泛化能力较强，在使用有限数据集训练的情况下，对实际数据中不同强度的噪声均可达到较好的压制效果。

参考文献:

[4]樊德华,段云卿,张婉,等.三维地震勘探抽油机噪声特征分析[J].石油地球物理勘探, 2008, 43(6):662-667.

[5]王鑫.基于检波器三维矢量组合的抽油机噪声压制方法研究[J].石油物探, 2011, 50(3):295-300, 306.

[6]赵彦青,萧蕴诗.基于盲源分离的抽油机噪声压制方法[J].计算机应用, 2014, 34(增刊1):349-351.

[7]李学良,王晓凯,赵昊,等.地震勘探中基于形态成分分析的抽油机噪声衰减方法[J].地球物理学进展,2017, 32(2):657-663.

[8]杨午阳,杨佳润,陈双全,等.基于U-Net深度学习网络的地震数据断层检测[J].石油地球物理勘探,2021, 56(4):688-697.

[9]蒋润,李振春,孙小东.人工智能地震数据重建方法现状分析[J].地球物理学进展, 2023, 38(5):2047-2062.

[10]曾丽丽,汤华贝,牛艺晓,等.基于LSTM堆叠残差网络的岩相识别方法[J].计算机与现代化, 2023(8):38-43.

[11]邱怡,包乾宗,马铭,等.基于U-Net网络的二维小波域地震数据随机噪声衰减[J].石油物探, 2023,62(5):878-890.

[12]韩卫雪,周亚同,池越.基于深度学习卷积神经网络的地震数据随机噪声去除[J].石油物探, 2018, 57(6):862-869, 877.

[13]王钰清,陆文凯,刘金林,等.基于数据增广和CNN的地震随机噪声压制[J].地球物理学报, 2019, 62(1):421-433.

[17]刘玉敏,魏海军,袁硕,等.基于卷积神经网络的地震数据去噪方法[J].吉林大学学报(信息科学版),2022, 40(2):231-239.

[18]韩建光,王卿,许媛,等.基于A-CGAN的深反射地震数据随机噪声压制方法研究[J].地质论评, 2024,70(1):228-238.

基金资助:四川省中央引导地方科技发展专项“岩石物理驱动的非常规油气地震多参数预测理论与方法研究”(2023ZYD0158); 四川省自然科学基金项目“四川盆地碳酸盐岩储层的脉冲神经网络识别理论及方法研究”(2023NSFSC0258)联合资助;

文章来源:马一凡,文武,薛雅娟,等.基于多尺度窗口生成器网络的抽油机噪声压制[J].石油地球物理勘探,2024,59(04):684-691.