一种叠前地震资料线性干扰压制方法及系统

摘要

本发明公开了一种叠前地震资料线性干扰压制方法及系统，读取原始地震炮集资料后，仅仅对读取的地震炮集数据中的小部分炮集数据基于线性干扰压制方法进行滤波处理，然后基于处理得到的有效波数据及滤除部分构建网络训练集，基于读取的全部数据构建测试集，利用二维卷积神经网络对构建得到的训练集学习得到线性干扰模型，再对读取的全部炮集数据进行线性干扰的有效压制，本发明的处理速度远高于对每个炮集数据进行分别处理的传统方法。

权利要求

1.一种叠前地震资料线性干扰压制方法，其特征在于，包括以下步骤：S1、读取原始地震炮集数据；

S2、从步骤S1的N个原始地震炮集数据中选择炮集数据，采用线性干扰压制方法对选择的炮集数据逐个进行滤波处理，得到选择的炮集数据中包含的有效波数据{S1,…,Sk}以及有效波数据{N1,…,Nk}；扫描分析N个原始地震炮集数据中能量最强的线性干扰波视倾角信息，并沿对应的视倾角方向拉平全部原始炮集数据以及所选择的炮集数据对应的有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}，得到拉平结果{X′1,…,X′N}，{S′1,…,S′k}，{N′1,…,N′k}；

S3、对步骤S2得到的{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}进行归一化的自动增益控制处理得到{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}；将S″j∈{S″1,…,S″k},N″j∈(1) (1) (q) (q) (1){N″1,…,N″k}按不同比例缩放q次得到{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…(1) (q) (q) (1) (1) (q) (q) (1) (1), N″k,…, N″1,…, N″k}，将{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…, N(q) (q)″k,…, N″1,…, N″k}逐个对应相加得到 作为网络训练(1) (1) (q) (q)

数据集中的输入数据集，并以{ S″1,…, S″k,…, S″1,…, S″k}作为网络训练数据集中的标签数据集；以{X″1,…,X″N}作为网络测试数据集；

(1) (1) (q)

S4、将步骤S3得到的 与{ S″1,…, S″k,…, S(q)

″1,…, S″k}输入二维卷积去噪神经网络进行学习，将学习得到的模型压制网络测试数据集{X″1,…,X″N}中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始炮集数据中对应的位置，完成对线性干扰的压制。

2.根据权利要求1所述的方法，其特征在于，步骤S2具体为：设读取的原始地震炮集数据总个数为N，应用线性干扰压制方法对原始地震炮集数据中任意k个炮集数据{X1,…,Xk}逐个进行滤波处理，0＜k＜N，得到对应的有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}；

对原始地震炮集数据{X1,…,XN}逐个进行视倾角扫描分析，估计任一Xi∈{X1,…,XN}中能量最强的线性干扰波的视倾角方向，逐个沿Xi对应的扫描得到的视倾角方向将Xi中的线性干扰波同相轴拉平得到{X′1,…,X′N}；逐个沿Xj∈{X1,…,Xk}对应扫描得到的视倾角方向将Xj对应的Sj,Nj中的线性干扰波同相轴拉平得到拉平后数据。

3.根据权利要求2所述的方法，其特征在于，X′i∈{X′1,…,X′N}，S′j∈{S′1,…,S′k}，N′j∈{N′1,…,N′k}具体为：(i) (j) (j)

其中，m,n分别为Xi的时间、空间维度， x′m,n为X′i中第m行第n列的数据， s′m,n,N′m,n分别为S′j,N′j中第m行第n列的数据。

4.根据权利要求1所述的方法，其特征在于，步骤S3具体为：S301、对步骤S2拉平后的数据{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}中的任一数据进行归一化的自动增益控制，得到处理后数据{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}，将{X″1,…,X″N}作为网络测试数据集；

S302、分别按比例对{S″1,…,S″k},{N″1,…,N″k}的每个数据进行缩放相加，缩放次数为(i1)q，每次{S″1,…,S″k}缩放比例为 P,i1＝1,…,q，缩放比例P∈[Pmin,Pmax]，对应的(i1){N″1,…,N″k}缩放比例为2‑ P，构建得到的网络训练数据集中的输入数据集为(1)并将 的{ S″1,…

(1) (q) (q)

, S″k,…, S″1,…, S″k}按作为网络训练数据集中的标签数据集。

5.根据权利要求4所述的方法，其特征在于，步骤S301中，X″i∈{X″1,…,X″N},S″j∈{S″1,…,S″k},N″j∈{N″1,…,N″k}具体为：X″i＝AGC[X′i]

S″j＝AGC[S′j]

N″j＝AGC[N′j]

其中，AGC[]为归一化自动增益控制算子。

6.根据权利要求4所述的方法，其特征在于，步骤S302中，缩放比例最大为Pmax＞1，最小为Pmin≥0。

(i1)

7.根据权利要求4所述的方法，其特征在于，步骤S302中，(i1) (i1)

S″j∈{ S″1,…, S″k}的表达式分别如下：

(j) (j)

其中， s″m,n, N″m,n分别为S″j,N″j中第m行第n列的数据，i1＝1…,q。

8.根据权利要求1所述的方法，其特征在于，步骤S4中，构建得到的网络训练集中的输(1) (1)入数据集为 网络训练集中的标签数据集为{ S″1,…, S″k,…(q) (q)

, S″1,…, S″k}，网络测试数据集为{X″1,…,X″N}，将输入数据集(1) (1) (q) (q)

与标签数据集{ S″1,…, S″k,…, S″1,…, S″k}分块输入二维去噪卷积神经网络，二维去噪卷积神经网络通过学习获取线性干扰模型，再对网络测试数据集{X″1,…,X″N}进行线性干扰压制，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理放回原始炮集数据中的对应位置，得到完成对线性干扰的压制后的结果{X″′1…,X″′N}。

9.一种叠前地震资料线性干扰压制系统，其特征在于，包括：读取模块，读取原始地震炮集数据；

拉平模块，从读取模块的N个原始地震炮集数据中选择炮集数据，采用线性干扰压制方法对炮集数据逐个进行滤波处理，得到有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}；扫描分析N个原始地震炮集数据中能量最强的线性干扰波视倾角信息，并沿对应的视倾角方向拉平全部原始炮集数据以及选择的炮集数据对应的有效波数据{S1,…,Sk}和线性干扰波数据{N1,…,Nk}，得到拉平结果{X′1,…,X′N}，{S′1,…,S′k}，{N′1,…,N′k}；

处理模块，对拉平模块得到的{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}进行归一化的自动增益控制处理得到{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}；将S″j∈{S″1,…,S″k},(1) (1) (q) (q) (1)N″j∈{N″1,…,N″k}按不同比例缩放q次得到{ S″1,…, S″k,…, S″1,…, S″k}和{ N(1) (q) (q) (1) (1) (q) (q) (1)″1,…, N″k,…, N″1,…, N″k}，将{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…(1) (q) (q), N″k,…, N″1,…, N″k}逐个对应相加得到 作为网络(1) (1) (q) (q)

训练数据集中的输入数据集，并以{ S″1,…, S″k,…, S″1,…, S″k}作为网络训练数据集中的标签数据集；以{X″1,…,X″N}作为网络测试数据集；

(1) (1)

压制模块，将处理模块得到的 与{ S″1,…, S″k,…(q) (q)

, S″1,…, S″k}输入二维卷积去噪神经网络进行学习，将学习得到的模型压制网络测试数据集{X″1,…,X″N}中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始炮集数据中的对应位置，完成对线性干扰的压制。

说明书

一种叠前地震资料线性干扰压制方法及系统

技术领域

[0001] 本发明属于地震勘探数据处理技术领域，具体涉及一种叠前地震资料线性干扰压制方法及系统。

背景技术

[0002] 地震资料中包含的干扰波主要分为两种：不规则干扰波和规则干扰波，其中在地震剖面上呈明显线性结构的规则干扰波称为线性干扰波。线性干扰波在从我国大部分路上采集区(例如：山地、沙丘区域以及地表覆盖较薄的西部地区)采集得到的地震资料中普遍存在。它的存在会严重降低地震资料的信噪比，甚至完全淹没有效信号，使后续一系列处理以及最终的地质解释不准确。例如在滕县煤田微山湖某矿区地震勘探项目中，采集的单炮记录由于较强的多次反射折射波(一种线性干扰波)对来自主煤层的有效信号的干涉作用，导致有效信号能量变弱，单炮记录信噪比低；沙丘鸣震线性干扰波往往淹没反射波场。因此，采取一定的技术方法对线性干扰波进行有效压制是十分必要的工作。

[0003] 传统压制线性干扰的方法如：频率波数域滤波、Cadzow滤波、拉东变换方法，虽取得了一定的压制效果，但当待处理的地震资料数量庞大，地下介质以及地形变化较大时，由于线性干扰波和有效信号在空间上的较大变化，传统方法压制线性干扰波的效果有限且处理时间过长。此时，在地震资料处理方面，深度学习中的神经网络方法在处理海量地震数据方面的高效性以及不会受限于某种特定的地形地质条件较传统方法有巨大的优势。并且，有监督的卷积神经网络(包括二维以及三维去噪卷积神经网络)已经被用于规则干扰波的压制，并在部分方面取得了超越传统方法的效果。

[0004] 现有技术：

[0005] Cadzow滤波方法是通过对输入数据堆砌高维矩阵加强信号特征，再基于对高维矩阵的奇异值分解得到对应数据中包含的不同信号的奇异值及其对应的奇异向量，通过分别用干扰波和有效信号对应的特征值和特征向量对信号进行重建来分离两者的一种方法。

[0006] 现有技术的缺点：

[0007] 1、用以重构有效信号和干扰波的奇异值个数需要反复试验获取，不准确会造成压制效果差或者损伤有效信号。

[0008] 2、处理海量地震数据时，计算速度慢。

发明内容

[0009] 本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种叠前地震资料线性干扰压制方法及系统，用于压制叠前地震资料线性干扰时可以避免参数反复试验选择，通过少量原始炮集数据构建泛化性较好的网络训练测试集，更快速地压制海量炮集数据中的线性干扰。

[0010] 本发明采用以下技术方案：

[0011] 一种叠前地震资料线性干扰压制方法，包括以下步骤：

[0012] S1、读取原始地震炮集数据；

[0013] S2、从步骤S1的N个原始地震炮集数据中选择炮集数据，采用线性干扰压制方法对炮集数据逐个进行滤波处理，得到其中包含的有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}；扫描分析N个原始地震炮集数据中能量最强的线性干扰波视倾角信息，并沿对应的视倾角方向拉平全部原始炮集数据以及选择的炮集数据对应的有效波数据{S1,…,Sk}和线性干扰波数据{N1,…,Nk}，得到拉平结果{X′1,…,X′N}，{S′1,…,S′k}，{N′1,…,N′k}；

[0014] S3、对步骤S2得到的{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}进行归一化的自动增益控制处理得到{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}；将S″j∈{S″1,…,S″k},N″j∈(1) (1) (q) (q) (1){N″1,…,N″k}按不同比例缩放q次得到{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…(1) (q) (q) (1) (1) (q) (q) (1) (1)
, N″k,…, N″1,…, N″k}，将{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…, N(q) (q)
″k,…, N″1,…, N″k}逐个对应相加得到 作为网络训练
数据集中的输入数据集，并以 作为网络训练数据集中的标
签数据集；以{X″1,…,X″N}作为网络测试数据集；

[0015] S4、将步骤S3得到的 与输入二维卷积去噪神经网络进行学习，将学习得到的模型压制网络测试数据集{X″1,…,X″N}中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始炮集数据中的对应位置，完成对线性干扰的压制。

[0016] 具体的，步骤S2具体为：

[0017] 设读取的原始地震炮集数据总个数为N，应用线性干扰压制方法对原始地震炮集数据中任意k个炮集数据{X1,…,Xk}逐个进行滤波处理，0＜k＜N，得到对应的有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}；

[0018] 对原始地震炮集数据{X1,…,XN}逐个进行视倾角扫描分析，估计任一Xi∈{X1,…,XN}中能量最强的线性干扰波的视倾角方向，逐个沿Xi对应的扫描得到的视倾角方向将Xi中的线性干扰波同相轴拉平得到{X′1,…,X′N}；逐个沿Xj∈{X1,…,Xk}对应扫描得到的视倾角方向将Xj对应的Sj,Nj中的线性干扰波同相轴拉平得到拉平后数据。

[0019] 进一步的，X′i∈{X′1,…,X′N}，S′j∈{S′1,…,S′k}，N′j∈{N′1,…,N′k}具体为：

[0020]

[0021]

[0022]

[0023] 其中，m,n分别为Xi的时间、空间维度，(i)x′m,n为X′i中第m行第n列的数据，(j)s′m,n(j), N′m,n分别为S′j,N′j中第m行第n列的数据。

[0024] 具体的，步骤S3具体为：

[0025] S301、对步骤S2拉平后的数据{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}中的任一数据进行归一化的自动增益控制，得到处理后数据{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}，将{X″1,…,X″N}作为网络测试数据集；

[0026] S302、分别按比例对{S″1,…,S″k},{N″1,…,N″k}的每个数据进行缩放相加，缩放(i1)次数为q，每次{S″1,…,S″k}缩放比例为 P,i1＝1,…,q，缩放比例P∈[Pmin,Pmax]，对应的(i1)
{N″1,…,N″k}缩放比例为2‑ P，构建得到的网络训练数据集中的输入数据集为
并将 的
按作为网络训练数据集中的标签数据集。

[0027] 进一步的，步骤S301中，X″i∈{X″1,…,X″N},S″j∈{S″1,…,S″k},N″j∈{N″1,…,N″k}具体为：

[0028] X″i＝AGC[X′i]

[0029] S″j＝AGC[S′j]

[0030] N″j＝AGC[N′j]

[0031] 其中，AGC[]为归一化自动增益控制算子。

[0032] 进一步的，步骤S302中，缩放比例最大为Pmax＞1，最小为Pmin≥0。

[0033] 进一步的，步骤S302中， (i1)S″j∈{(i1)S″1,…,(i1)S″k}的表达式分别如下：

[0034]

[0035]

[0036] 其中，(j)s″m,n,(j)N″m,n分别为S″j,N″j中第m行第n列的数据，i1＝1…,q。

[0037] 具体的，步骤S4中，构建得到的网络训练集中的输入数据集为网络训练集中的标签数据集为 网络测试数据集为{X″1,…,X″N}，
(1) (1) (q) (q)
将输入数据集 与标签数据集{ S″1,…, S″k,…, S″1,…, S″k}
分块输入二维去噪卷积神经网络，二维去噪卷积神经网络通过学习获取线性干扰模型，再对网络测试数据集{X″1,…,X″N}进行线性干扰压制，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始炮集数据中的对应位置得到完成对线性干扰的压制后的结果{X″′1…,X″′N}。

[0038] 本发明的另一技术方案是，一种叠前地震资料线性干扰压制系统，包括：

[0039] 读取模块，读取原始地震炮集数据；

[0040] 拉平模块，从读取模块的N个原始地震炮集数据中选择炮集数据，采用线性干扰压制方法对炮集数据逐个进行滤波处理，得到选择的炮集数据中包含的有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}；扫描分析N个原始地震炮集数据中能量最强的线性干扰波视倾角信息，并沿对应的视倾角方向拉平全部原始炮集数据以及选择的炮集数据中包含的有效波数据{S1,…,Sk}和线性干扰波数据{N1,…,Nk}，得到拉平结果{X′1,…,X′N}，{S′1,…,S′k}，{N′1,…,N′k}；

[0041] 处理模块，对拉平模块得到的{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}进行归一化的自动增益控制处理得到{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}；将S″j∈{S″1,…,(1) (1) (q) (q)S″k},N″j∈{N″1,…,N″k}按不同比例缩放q次得到{ S″1,…, S″k,…, S″1,…, S″k}和(1) (1) (q) (q) (1) (1) (q) (q) (1)
{ N″1,…, N″k,…, N″1,…, N″k}，将{ S″1,…, S″k,…, S″1,…, S″k}和{ N(1) (q) (q)
″1,…, N″k,…, N″1,…, N″k}逐个对应相加得到 作为
(1) (1) (q) (q)
网络训练数据集中的输入数据集，并以{ S″1,…, S″k,…, S″1,…, S″k}作为网络训练数据集中的标签数据集；以{X″1,…,X″N}作为网络测试数据集；

[0042] 压制模块，将处理模块得到的 与输入二维卷积去噪神经网络进行学习，将学习得到的模型压制网络测试数据集{X″1,…,X″N}中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始炮集数据中的对应位置，完成对线性干扰的压制。

[0043] 与现有技术相比，本发明至少具有以下有益效果：

[0044] 本发明一种叠前地震资料线性干扰压制方法，读取原始地震炮集资料后，仅仅对读取的地震炮集数据中的小部分炮集数据基于线性干扰压制方法进行滤波处理。基于处理得到的有效波数据和线性干扰波数据构建网络训练集，并基于读取的全部数据构建测试集。通过对一小部分数据进行滤波处理构建得到的训练集学习得到的线性干扰模型对读取的全部炮集数据进行线性干扰的有效压制，使得神经网络的处理速度远高于对每个炮集数据进行分别处理的传统方法。

[0045] 进一步的，步骤S2中仅用线性干扰压制方法对任意k个炮集数据进行滤波处理，可以快速获得少量有效波数据用以构建较为理想的标签数据。

[0046] 进一步的，在步骤S2中通过拉平处理加强线性干扰波的沿地震道空间方向的相干特征，便于后期网络对线性干扰特征进行提取，获得更为完备的线性干扰波模型。

[0047] 进一步的，设置X′i∈{X′1,…,X′N}，S′j∈{S′1,…,S′k}，N′j∈{N′1,…,N′k}用于后续计算。

[0048] 进一步的，步骤S3对现有的线性干扰压制方法估计得到的有效波数据和线性干扰波数据进行放缩相加构建多种信噪比不同的网络训练集，增强网络泛化性。

[0049] 进一步的，在步骤S301中通过对经过拉平及归一化的自动增益控制处理读取的全部地震炮集数据以及选取的部分炮集数据经过滤波处理得到的有效波数据和线性干扰波数据，放大线性干扰的细节特征，便于后期网络对线性干扰特征进行充分提取。

[0050] 进一步的，通过步骤S302利用对经过一系列处理的有效波数据和线性干扰波数据进行放缩组合，构建信噪比不同的训练集，加强网络泛化性。

[0051] 进一步的，步骤S4利用二维卷积神经网络对构建的网络训练数据集进行学习，得到较为完备的线性干扰波模型，基于此模型可以快速对基于全部原始炮集数据构建网络测试集进行处理，极大地提高压制效率。

[0052] 综上所述，本发明用于压制叠前地震资料线性干扰时可以避免参数反复试验选择，通过少量原始炮集数据构建多种信噪比的网络训练测试集，训练二维卷积神经网络得到泛化性较好的线性干扰波模型，更快速地压制海量炮集数据中的线性干扰，并且不会因地形地质条件变化导致的有效信号及线性干扰空间变化较大而降低压制效果。

[0053] 下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

[0054] 图1为某油田的一个共含有17个炮集数据的原始地震资料中前4个炮集数据剖面图；

[0055] 图2为图1所示数据经过线性干扰压制方法处理后的有效波数据剖面图；

[0056] 图3为图1所示数据与图2所示数据的差值剖面图；

[0057] 图4为本发明处理此实际炮集数据流程中部分结果展示图，其中，(a)为利用图1所示数据构建的网络测试集；(b)为利用图2所示数据构建的网络训练集中的标签数据集；(c)为利用图2所示数据以及图3所示数据构建的网络训练集中的输入数据集；

[0058] 图5为测试本发明所使用的二维卷积神经网络结构图

[0059] 图6为本发明估计得到的原炮集数据中包含的线性干扰以及通过本发明压制线性干扰后的炮集数据剖面图，其中，(a)为图1所示数据通过本发明压制线性干扰后的炮集数据剖面图，(b)为本发明估计得到的图(a)中包含的线性干扰剖面图；

[0060] 图7为本发明流程图。

具体实施方式

[0061] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0062] 在本发明的描述中，需要理解的是，术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在，但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。

[0063] 还应当理解，在本发明说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本发明。如在本发明说明书和所附权利要求书中所使用的那样，除非上下文清楚地指明其它情况，否则单数形式的“一”、“一个”及“该”意在包括复数形式。

[0064] 还应当进一步理解，在本发明说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合，并且包括这些组合。

[0065] 在附图中示出了根据本发明公开实施例的各种结构示意图。这些图并非是按比例绘制的，其中为了清楚表达的目的，放大了某些细节，并且可能省略了某些细节。图中所示出的各种区域、层的形状及它们之间的相对大小、位置关系仅是示例性的，实际中可能由于制造公差或技术限制而有所偏差，并且本领域技术人员根据实际所需可以另外设计具有不同形状、大小、相对位置的区域/层。

[0066] 请参阅图6，本发明提供了一种叠前地震资料线性干扰压制方法，构建网络训练测试集并输入二维去噪卷积神经网络中学习训练集中的线性干扰模型并基于此模型对测试集中的线性干扰进行压制；得到压制线性干扰后的炮集数据用于后期速度分析、叠加处理、成像、地质结构分析。

[0067] 请参阅图7，本发明一种叠前地震资料线性干扰压制方法，包括以下步骤：

[0068] S1、读取原始地震炮集资料；

[0069] 请参阅图1，其为读取的实际原始地震炮集数据{X1,…,XN},N＝17中的前四个炮集数据(数据量大仅取前四个作为展示)，单个炮集数据时间空间维度均分别为3501ms、256道。

[0070] S2、设读取的原始地震炮集数据总个数为N，任意选择{X1,…,XN}中部分空间维度及时间维度一致的炮集数据{X1,…,Xk},k＝5，应用线性干扰压制方法对{X1,…,Xk}逐个进行滤波处理，得到如图2所示的对应的有效波数据{S1,…,Sk}中前四个结果以及如图3所示的线性干扰波数据{N1,…,Nk}中的前四个结果，再对{X1,…,XN}逐个进行视倾角扫描分析，估计任一Xi∈{X1,…,XN}中能量最强的线性干扰波的视倾角方向以及污染区域，逐个沿Xi对应的扫描得到的视倾角方向将Xi中的污染区域线性干扰波同相轴拉平得到{X′1,…,X′N}；逐个沿Xj∈{X1,…,Xk}对应的扫描得到的视倾角方向将Xj对应的Sj,Nj中相同区域内的线性干扰波同相轴拉平得到拉平后数据{S′1,…,S′k}，{N′1,…,N′k}；

[0071] 数据X′i∈{X′1,…,X′N},S′j∈{S′1,…,S′k},N′j∈{N′1,…,N′k}表达式如下：

[0072]

[0073]

[0074]

[0075] 其中，m,n分别为Xi的时间、空间维度，(i)x′m,n为X′i中第m行第n列的数据，(j)s′m,n(j), N′m,n分别为S′j,N′j中第m行第n列的数据。

[0076] S3、对步骤S2得到的{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}进行归一化的自动增益控制处理得到{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}；将S″j∈{S″1,…,S″k},N″j∈(1) (1) (q) (q) (1){N″1,…,N″k}按不同比例缩放q次得到{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…(1) (q) (q) (1) (1) (q) (q) (1) (1)
, N″k,…, N″1,…, N″k}，将{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…, N(q) (q)
″k,…, N″1,…, N″k}逐个对应相加得到 作为网络训练
数据集中的输入数据集，并以 作为网络训练数据集中的标
签数据集；以{X″1,…,X″N}作为网络测试数据集；

[0077] S301、对{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}中的任一数据进行归一化的自动增益控制，得到处理后数据{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}，{X″1,…,X″N}为网络测试数据集(其中前四个结果如图4(a)所示)；

[0078] X″i∈{X″1,…,X″N},S″j∈{S″1,…,S″k},N″j∈{N″1,…,N″k}表达式分别如下：

[0079] X″i＝AGC[X′i]

[0080] S″j＝AGC[S′j]

[0081] N″j＝AGC[N′j]

[0082] 其中，AGC[]为归一化自动增益控制算子。

[0083] S302、分别按比例对{S″1,…,S″k},{N″1,…,N″k}每个数据进行缩放相加，缩放比(1) (2) (3) (4)例最大为2，缩放次数为q＝5，每次{S″1,…,S″k}缩放比例为 P＝1.8， P＝0.9， P＝1，(5) (1) (2) (3) (4)
P＝0.2， P＝1.1，则对应的{N″1,…,N″k}缩放比例为 P＝0.2， P＝1.1， P＝1， P＝(5)
1.8， P＝0.9，构建得到的网络训练数据的输入数据集为
(1) (1) (q)
(其中前四个结果如图4(c)所示)，并将与其对应的按比例缩放后的{ S″1,…, S″k,…,(q)
S″1,…, S″k}(其中前四个结果如图4(b)所示)作为网络训练数据集中的含标签数据集。

[0084] (i1)S″j∈{(i1)S″1,…,(i1)S″k},i1＝1…,q的表达式分别如下：

[0085]

[0086]

[0087] 其中，(j)s″m,n,(j)N″m,n分别为S″j,N″j中第m行第n列的数据。

[0088] S4、将步骤S3得到的 与{(1)S″1,…,(1)S″k,…,(q)(q)
S″1,…, S″k}输入二维卷积去噪神经网络进行学习，将学习得到的模型压制网络测试数据集{X″1,…,X″N}中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始炮集数据中对应的位置，完成对线性干扰的压制。

[0089] 构建得到的网络训练集中的输入数据集为 标签(1) (1) (q) (q)
数据集为{ S″1,…, S″k,…, S″1,…, S″k}，网络测试数据集为{X″1,…,X″N}，将(1) (1) (q) (q)
与{ S″1,…, S″k,…, S″1,…, S″k}分块(自定义分
块的时间空间维度)输入二维去噪卷积神经网络，网络通过对两者间的残差学习获取线性干扰模型，再对{X″1,…,X″N}进行线性干扰压制，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始炮集数据中对应的位置，得到完成对线性干扰的压制后的结果{X″′1…,X″′N}。

[0090] 本发明再一个实施例中，提供一种叠前地震资料线性干扰压制系统，该系统能够用于实现上述叠前地震资料线性干扰压制方法，具体的，该叠前地震资料线性干扰压制系统包括读取模块、拉平模块、处理模块以及压制模块。

[0091] 其中，读取模块，读取原始地震炮集数据；

[0092] 拉平模块，从读取模块的N个原始地震炮集数据中选择炮集数据，采用线性干扰压制方法对炮集数据逐个进行滤波处理，得到有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}；扫描分析N个原始地震炮集数据中能量最强的线性干扰波视倾角信息，并沿对应的视倾角方向拉平全部原始炮集数据，得到拉平结果{X′1,…,X′N}，{S′1,…,S′k}，{N′1,…,N′k}；

[0093] 处理模块，对拉平模块得到的{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}进行归一化的自动增益控制处理得到{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}；将S″j∈{S″1,…,(1) (1) (q) (q)S″k},N″j∈{N″1,…,N″k}按不同比例缩放q次得到{ S″1,…, S″k,…, S″1,…, S″k}和(1) (1) (q) (q) (1) (1) (q) (q) (1)
{ N″1,…, N″k,…, N″1,…, N″k}，将{ S″1,…, S″k,…, S″1,…, S″k}和{ N(1) (q) (q)
″1,…, N″k,…, N″1,…, N″k}逐个对应相加得到 作为
(1) (1) (q) (q)
网络训练数据集中的输入数据集，并以{ S″1,…, S″k,…, S″1,…, S″k}作为网络训练数据集中的标签数据集；以{X″1,…,X″N}作为网络测试数据集；

[0094] 压制模块，将处理模块得到的 与输入二维卷积去噪神经网络进行学习，将学习得到的模型压制网络测试数据集{X″1,…,X″N}中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理放回原始炮集数据中的对应位置，完成对线性干扰的压制。

[0095] 本发明再一个实施例中，提供了一种终端设备，该终端设备包括处理器以及存储器，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器用于执行所述计算机存储介质存储的程序指令。处理器可能是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor、DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field‑Programmable GateArray，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等，其是终端的计算核心以及控制核心，其适于实现一条或一条以上指令，具体适于加载并执行一条或一条以上指令从而实现相应方法流程或相应功能；本发明实施例所述的处理器可以用于叠前地震资料线性干扰压制方法的操作，包括：

[0096] 读取原始地震炮集数据；从N个原始地震炮集数据中选择炮集数据，采用线性干扰压制方法对炮集数据逐个进行滤波处理，得到有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}；扫描分析N个原始地震炮集数据中能量最强的线性干扰波视倾角信息，并沿对应的视倾角方向拉平全部原始炮集数据，得到拉平结果{X′1,…,X′N}，{S′1,…,S′k}，{N′1,…,N′k}；对{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}进行归一化的自动增益控制处理得到{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}；将S″j∈{S″1,…,S″k},N″j∈{N″1,…,N(1) (1) (q) (q) (1) (1)″k}按不同比例缩放q次得到{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…, N″k,…(q) (q) (1) (1) (q) (q) (1) (1) (q)
, N″1,…, N″k}，将{ S″1,…, S″k,…, S″1,…, S″k}和{ N″1,…, N″k,…, N(q)
″1,…, N″k}逐个对应相加得到 作为网络训练数据集中的
(1) (1) (q) (q)
输入数据集，并以{ S″1,…, S″k,…, S″1,…, S″k}作为网络训练数据集中的标签数(1)
据集；以{X″1,…,X″N}作为网络测试数据集；将 与{ S
(1) (q) (q)
″1,…, S″k,…, S″1,…, S″k}输入二维卷积去噪神经网络进行学习，将学习得到的模型压制网络测试数据集{X″1,…,X″N}中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始数据中对应的位置，完成对线性干扰的压制。

[0097] 本发明再一个实施例中，本发明还提供了一种存储介质，具体为计算机可读存储介质(Memory)，所述计算机可读存储介质是终端设备中的记忆设备，用于存放程序和数据。可以理解的是，此处的计算机可读存储介质既可以包括终端设备中的内置存储介质，当然也可以包括终端设备所支持的扩展存储介质。计算机可读存储介质提供存储空间，该存储空间存储了终端的操作系统。并且，在该存储空间中还存放了适于被处理器加载并执行的一条或一条以上的指令，这些指令可以是一个或一个以上的计算机程序(包括程序代码)。
需要说明的是，此处的计算机可读存储介质可以是高速RAM存储器，也可以是非不稳定的存储器(non‑volatile memory)，例如至少一个磁盘存储器。

[0098] 可由处理器加载并执行计算机可读存储介质中存放的一条或一条以上指令，以实现上述实施例中有关叠前地震资料线性干扰压制方法的相应步骤；计算机可读存储介质中的一条或一条以上指令由处理器加载并执行如下步骤：

[0099] 读取原始地震炮集数据；从N个原始地震炮集数据中选择炮集数据，采用线性干扰压制方法对炮集数据逐个进行滤波处理，得到选择的炮集数据中包含的有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}；扫描分析N个原始地震炮集数据中能量最强的线性干扰波视倾角信息，并沿对应的视倾角方向拉平全部原始炮集数据以及选择的炮集数据中包含的有效波数据{S1,…,Sk}以及线性干扰波数据{N1,…,Nk}，得到拉平结果{X′1,…,X′N}，{S′1,…,S′k}，{N′1,…,N′k}；对{X′1,…,X′N},{S′1,…,S′k},{N′1,…,N′k}进行归一化的自动增益控制处理得到{X″1,…,X″N},{S″1,…,S″k},{N″1,…,N″k}；将S″j∈{S″1,…,S(1) (1) (q) (q)
″k},N″j∈{N″1,…,N″k}按不同比例缩放q次得到{ S″1,…, S″k,…, S″1,…, S″k}和(1) (1) (q) (q) (1) (1) (q) (q) (1)
{ N″1,…, N″k,…, N″1,…, N″k}，将{ S″1,…, S″k,…, S″1,…, S″k}和{ N(1) (q) (q)
″1,…, N″k,…, N″1,…, N″k}逐个对应相加得到 作为
网络训练数据集中的输入数据集，并以 作为网络训练数据
集中的标签数据集；以{X″1,…,X″N}作为网络测试数据集；将
(1) (1) (q) (q)
与{ S″1,…, S″k,…, S″1,…, S″k}输入二维卷积去噪神经网络进行学习，将学习得到的模型压制网络测试数据集{X″1,…,X″N}中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理后放回原始数据中对应的位置，完成对线性干扰的压制。

[0100] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中的描述和所示的本发明实施例的组件可以通过各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0101] 将构建的网络训练集分为若干个空间时间维度相等的分块(分块时间维度为20ms，空间维度127道)输入二维卷积去噪神经网络进行学习，此次测试本发明的二维卷积去噪神经网络结构如图5所示，将学习得到的模型用于压制构建的网络测试数据集中的线性干扰，对压制线性干扰后的网络测试数据集进行逆归一化自动增益控制以及反拉平处理完成对线性干扰的压制，得到的本发明压制线性干扰后的炮集数据中的前四个结果剖面图如图6(a)所示，本发明估计得到的原炮集数据中前四个数据包含的线性干扰剖面图如图6(b)所示。

[0102] 利用本发明一种压制叠前地震资料线性干扰的神经网络方法，可以避免参数反复试验选择，并且通过少量原始炮集数据构建泛化性较好的网络训练测试集，快速压制海量炮集数据中的线性干扰，并且不会因地形地质条件变化导致的有效信号及线性干扰空间变化较大而降低压制效果。

[0103] 综上所述，本发明一种叠前地震资料线性干扰压制方法及系统具有以下优点：

[0104] 1，可以避免参数反复试验选择，并且通过少量原始炮集数据构建泛化性较好的网络训练测试集，快速压制海量炮集数据中的线性干扰；

[0105] 2，不会因地形地质条件变化导致的有效信号及线性干扰空间变化较大而降低压制效果线性。

[0106] 本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD‑ROM、光学存储器等)上实施的计算机程序产品的形式。

[0107] 本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

[0108] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

[0109] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

[0110] 以上内容仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明权利要求书的保护范围之内。