地震资料的水平线性干扰去除方法及装置与流程

[0001]
本发明涉及地震勘探数据处理技术领域，尤其涉及地震资料的水平线性干扰去除方法及装置。


背景技术：

[0002]
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
[0003]
地震资料在处理过程中，由于种种原因会在叠加剖面上存在一些水平线性干扰，如水平线性噪声。这些水平线性干扰的存在严重影响了地震资料的最终成像效果，进而影响最终的地层构造解释。而现阶段，消除这些水平线性干扰通常是基于速度或者基于时差进行识别进而加以去除。但地震资料的水平同相轴的速度是接近无穷大的，对应的时差接近于零。此外，地震资料中的水平线性干扰很难通过去线性干扰处理加以去除，如频率波数域变换(fk变换)等方法。
[0004]
总之，现有技术下地震资料去除水平线性干扰的处理效果不佳，地震资料中的水平线性干扰难以消除。


技术实现要素：

[0005]
本发明实施例提供一种地震资料的水平线性干扰去除方法，用以改善地震资料去除水平线性干扰的处理效果，该方法包括：
[0006]
从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗；
[0007]
计算水平线性干扰待处理区域的时窗的应用时差值；所述应用时差值用于表征时窗的最大值和最小值的差值；
[0008]
对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料；
[0009]
根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料；所述正向校正处理用于对地震资料进行基于地层倾角的正向变换；
[0010]
对地震正变换资料进行去线性干扰处理；
[0011]
对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料；所述反向校正处理用于对地震资料进行基于地层倾角的反向变换；
[0012]
对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，得到消除水平线性干扰后的地震资料。
[0013]
本发明实施例还提供一种地震资料的水平线性干扰去除装置，用以改善地震资料去除水平线性干扰的处理效果，该装置包括：
[0014]
时窗获取模块，用于从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗；
[0015]
应用时差值计算模块，用于计算水平线性干扰待处理区域的时窗的应用时差值；所述应用时差值用于表征时窗的最大值和最小值的差值；
[0016]
时移处理和地震道长度延长处理模块，用于对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料；
[0017]
正向校正处理模块，用于根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料；所述正向校正处理用于对地震资料进行基于地层倾角的正向变换；
[0018]
去线性干扰处理模块，用于对地震正变换资料进行去线性干扰处理；
[0019]
反向校正处理模块，用于对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料；所述反向校正处理用于对地震资料进行基于地层倾角的反向变换；
[0020]
反向时移处理和地震道长度缩减处理模块，用于对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，得到消除水平线性干扰后的地震资料。
[0021]
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地震资料的水平线性干扰去除方法。
[0022]
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述地震资料的水平线性干扰去除方法的计算机程序。
[0023]
本发明实施例中，从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗；计算水平线性干扰待处理区域的时窗的应用时差值；所述应用时差值用于表征时窗的最大值和最小值的差值；对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料；根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料；所述正向校正处理用于对地震资料进行基于地层倾角的正向变换；对地震正变换资料进行去线性干扰处理；对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料；所述反向校正处理用于对地震资料进行基于地层倾角的反向变换；对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，得到消除水平线性干扰后的地震资料，与现有技术中地震资料的水平线性干扰很难通过去线性干扰处理加以去除的技术方案相比，通过时移处理、地震道长度延长处理和正向校正处理，可将地震资料的水平线性干扰改变成其他角度，可使去线性干扰处理间接性地发挥作用，极大地提高了去线性干扰处理的可行性，从而改善了地震资料去除水平线性干扰的处理效果。
附图说明
[0024]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
[0025]
图1为本发明实施例中一种地震资料的水平线性干扰去除方法的流程示意图；
[0026]
图2为本发明实施例中一种地震资料的水平线性干扰去除方法的具体示例图；
[0027]
图3为本发明实施例中一种地震资料的水平线性干扰去除方法的具体示例图；
[0028]
图4为本发明实施例中一种地震资料的水平线性干扰去除方法的具体示例图；
[0029]
图5为本发明实施例中一种地震资料的水平线性干扰去除装置的结构示意图。
具体实施方式
[0030]
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
[0031]
图1为本发明实施例中一种地震资料的水平线性干扰去除方法的流程示意图，如图1所示，本发明实施例提供的地震资料的水平线性干扰去除方法，可以包括：
[0032]
步骤101：从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗；
[0033]
步骤102：计算水平线性干扰待处理区域的时窗的应用时差值；应用时差值用于表征时窗的最大值和最小值的差值；
[0034]
步骤103：对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料；
[0035]
步骤104：根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料；正向校正处理用于对地震资料进行基于地层倾角的正向变换；
[0036]
步骤105：对地震正变换资料进行去线性干扰处理；
[0037]
步骤106：对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料；反向校正处理用于对地震资料进行基于地层倾角的反向变换；
[0038]
步骤107：对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，得到消除水平线性干扰后的地震资料。
[0039]
具体实施时，首先从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗。
[0040]
实施例中，可通过直接在叠加剖面上选取要进行水平线性干扰待处理的区域来获取水平线性干扰待处理区域的时窗。
[0041]
在上述实施例中，可通过从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗，有助于在后续步骤中对水平线性干扰待处理区域进行处理。
[0042]
具体实施时，在从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗后，计算水平线性干扰待处理区域的时窗的应用时差值；应用时差值用于表征时窗的最大值和最小值的差值。
[0043]
实施例中，计算水平线性干扰待处理区域的时窗的应用时差值，可以包括：获取水平线性干扰待处理区域的时窗的最大值和最小值；根据获取的时窗的最大值和最小值，计算时窗的应用时差值。
[0044]
在上述实施例中，通过确定水平线性干扰待处理区域的所涉及到时窗的最小时间和最大时间，计算时窗的应用时差值，，有助于在后续步骤中对水平线性干扰待处理区域进
行处理。
[0045]
具体实施时，在计算水平线性干扰待处理区域的时窗的应用时差值；应用时差值用于表征时窗的最大值和最小值的差值后，对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料。
[0046]
实施例中，对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料，可以包括：对地震资料中地震道以应用时差值的第一预设倍数，进行向下时移处理；对向下时移处理后的地震道以应用时差值的第二预设倍数，进行地震道长度延长处理，得到处理后的地震资料。
[0047]
在上述实施例中，第一预设倍数可为一倍；第二预设倍数可为两倍。地震资料中地震道以应用时差值的第一预设倍数，进行向下时移处理；对向下时移处理后的地震道以应用时差值的第二预设倍数，进行地震道长度延长处理，可以包括：将地震道整体向下时移一倍应用时差值长度，并且将道长度延长两倍应用时差值长度。具体实施时，在对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料后，根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料；正向校正处理用于对地震资料进行基于地层倾角的正向变换。
[0048]
实施例中，根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料，可以包括：对处理后的地震资料，沿预设地层同相轴方向拾取平滑曲线，计算该平滑曲线经每一地震道的拾取时长与该平滑曲线平均值的差值；计算该差值的平均值；根据该差值的平均值，对处理后的地震资料进行正向校正处理，得到地震正变换资料。
[0049]
在上述实施例中，沿预设地层同相轴方向拾取平滑曲线，可以包括：沿某一地层倾角方向拾取平滑的曲线，其他不需要进行水平线性干扰去除处理的部分进行适当平滑外推，直至全部拾取完成；求取所拾取平滑曲线的平均值；求取平滑曲线上每一个值和该平滑曲线平均值的平均时差；使用该平均时差对处理后的地震资料进行正向校正处理。由此实现了对地震数据的地层倾角变换的正向变换。在地层倾角正向变换之后的地震数据上做相应的地震资料处理，可将地震资料的水平线性干扰改变成其他角度，可使去线性干扰处理间接性地发挥作用，极大地提高了去线性干扰处理的可行性，为fk变换等去线性处理奠定基础。
[0050]
具体实施时，在根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料后，对地震正变换资料进行去线性干扰处理。
[0051]
在上述实施例中，可通过去线性干扰处理对地震正变换资料中水平线性干扰加以去除，如通过频率波数域变换(fk变换)等方法。在地震正变换资料上可进行去线性等处理去除原先水平，经正向校正处理后倾斜的水平线性干扰。
[0052]
具体实施时，在对地震正变换资料进行去线性干扰处理后，对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料；反向校正处理用于对地震资料进行基于地层倾角的反向变换。
[0053]
实施例中，对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料，可以包括：根据该差值的平均值，对去线性干扰处理后的地震正变换资料进行
反向校正处理，得到地震反变换资料。
[0054]
在上述实施例中，使用该差值的平均值对去线性干扰处理后的地震正变换资料进行时差的反向应用。
[0055]
具体实施时，在对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料后，对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，得到消除水平线性干扰后的地震资料。
[0056]
实施例中，对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，可以包括：对地震反变换资料中地震道以应用时差值的第一预设倍数，进行向上时移处理；对向上时移处理后的地震道以应用时差值的第二预设倍数，进行地震道长度缩短处理，得到处理后的地震资料。
[0057]
在上述实施例中，将地震反变换资料整体向上时移一倍应用时差值长度，并且将道长度缩短至原先道长度，由此实现了地震数据的地层倾角变换的反向变换。通过进行地层倾角变换的反变换，可使地震反变换资料回到原始的地震数据的构造状态，经过fk变换等处理之后，水平的线性干扰，如噪声，可得到很好的消除。
[0058]
下面给出一个具体实施例，来说明本发明的方法的具体应用。该具体实施例，可以包括如下步骤：
[0059]
1、如图2所示，图2显示了地震资料的叠加剖面；首先在地震资料的叠加剖面上确定水平线性干扰待处理区域的时窗；图2中纵坐标是时间轴，英文是time；横坐标是地震资料叠加剖面的“共中心点”顺序号，英文简写是cmp或者cdp，全拼为common midpoint；
[0060]
2、在地震资料的叠加剖面中水平线性干扰待处理区域中，查看水平线性干扰待处理区域的时窗所涉及到的最小时间t
min
和最大时间t
max
；
[0061]
3、确定最小时间和最大时间的时差值为应用时差值t，其中，t＝t
max-t
min
；
[0062]
4、将地震资料中地震道的叠加剖面数据整体向下时移一倍应用时差值t，并且将地震道长度延长两倍应用时差值形成处理后的地震资料；在该处理后的地震资料中道长度t
trace
＝t
trace0
+2t，其中，t
trace0
为地震资料的原始地震道长度；t
trace
为处理后的地震资料的地震道长度；
[0063]
5、沿处理后的地震资料的水平线性干扰待处理区域的某一倾斜同相轴拾取平滑的曲线，其他不需要处理部分进行适当平滑外推，直至全部拾取完成；
[0064]
6、求取所拾取平滑曲线的平均值t
mean
，其中ti为曲线上每一道的拾取时间；n为曲线上的地震道的总道数，
[0065]
7、求取平滑曲线上每一地震道的时间数值t
i
和平均值t
mean
的平均差值δt
i
；
[0066]
8、使用该差值δt
i
对地震数据每一地震道进行校正，这样处理后的地震资料就会沿着所拾取的曲线进行拉平，进而形成一个地震正变换资料，如图3所示；
[0067]
9、经过以上各步处理之后就实现了地震资料地层倾角变换的正变换；
[0068]
10、在地震正变换资料上进行fk等处理去除原先为水平状态，现在为倾斜状态的水平线性干扰，形成去线性干扰处理后的地震正变换资料；
[0069]
11、将去线性干扰处理后的地震正变换资料使用之前的平均时差δt
i
进行反应用，即反向校正处理，形成地震反变换资料；
[0070]
12、将地震反变换资料分别进行反向时移处理和地震道长度缩减处理，如整体向上时移一倍应用时差值t，并且将道长度缩短至原始地震道长度；
[0071]
13、经过以上两步就实现了地层倾角变换的反变换，得到消除水平线性干扰后的地震资料，如图4所示。图4经过地层倾角变换正变换之后的fk去除线性，有效消除了地震资料的水平噪声。
[0072]
当然，可以理解的是，上述详细流程还可以有其他变化例，相关变化例均应落入本发明的保护范围。
[0073]
本发明实施例中，从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗；计算水平线性干扰待处理区域的时窗的应用时差值；应用时差值用于表征时窗的最大值和最小值的差值；对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料；根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料；正向校正处理用于对地震资料进行基于地层倾角的正向变换；对地震正变换资料进行去线性干扰处理；对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料；反向校正处理用于对地震资料进行基于地层倾角的反向变换；对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，得到消除水平线性干扰后的地震资料，与现有技术中地震资料的水平线性干扰很难通过去线性干扰处理加以去除的技术方案相比，通过时移处理、地震道长度延长处理和正向校正处理，可将地震资料的水平线性干扰改变成其他角度，可使去线性干扰处理间接性地发挥作用，极大地提高了去线性干扰处理的可行性，从而改善了地震资料去除水平线性干扰的处理效果。
[0074]
如上所示，本发明实施例本方法能够改变地震资料的地层倾角分布特征，从而改变地层倾角的角度，从而使fk去除线性得到好的处理效果。
[0075]
本发明实施例中还提供了一种的水平线性干扰去除装置，如下面的实施例。由于该装置解决问题的原理与的水平线性干扰去除方法相似，因此该装置的实施可以参见的水平线性干扰去除方法的实施，重复之处不再赘述。
[0076]
图5为本发明实施例中一种地震资料的水平线性干扰去除装置的结构示意图，如图5所示，本发明实施例提供的地震资料的水平线性干扰去除装置，可以包括：
[0077]
时窗获取模块01，用于从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗；
[0078]
应用时差值计算模块02，用于计算水平线性干扰待处理区域的时窗的应用时差值；应用时差值用于表征时窗的最大值和最小值的差值；
[0079]
时移处理和地震道长度延长处理模块03，用于对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料；
[0080]
正向校正处理模块04，用于根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料；正向校正处理用于对地震资料进行基于地层倾角的正向变换；
[0081]
去线性干扰处理模块05，用于对地震正变换资料进行去线性干扰处理；
[0082]
反向校正处理模块06，用于对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料；反向校正处理用于对地震资料进行基于地层倾角的反向变换；
[0083]
反向时移处理和地震道长度缩减处理模块07，用于对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，得到消除水平线性干扰后的地震资料。
[0084]
在一个实施例中，应用时差值计算模块，具体用于：获取水平线性干扰待处理区域的时窗的最大值和最小值；
[0085]
根据获取的时窗的最大值和最小值，计算时窗的应用时差值。
[0086]
在一个实施例中，时移处理和地震道长度延长处理模块，具体用于：
[0087]
对地震资料中地震道以应用时差值的第一预设倍数，进行向下时移处理；
[0088]
对向下时移处理后的地震道以应用时差值的第二预设倍数，进行地震道长度延长处理，得到处理后的地震资料。
[0089]
在一个实施例中，反向时移处理和地震道长度缩减处理模块，具体用于：
[0090]
对地震反变换资料以应用时差值的第一预设倍数，进行向上时移处理；
[0091]
对向上时移处理后的地震反变换资料以应用时差值的第二预设倍数，进行地震道长度缩减处理，得到消除水平线性干扰后的地震资料。
[0092]
在一个实施例中，正向校正处理模块，具体用于：
[0093]
对处理后的地震资料，沿预设地层同相轴方向拾取平滑曲线，计算该平滑曲线经每一地震道的拾取时长与该平滑曲线平均值的差值；
[0094]
计算该差值的平均值；
[0095]
根据该差值的平均值，对处理后的地震资料进行正向校正处理，得到地震正变换资料。
[0096]
在一个实施例中，反向校正处理模块，具体用于：
[0097]
根据该该差值的平均值，对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料。
[0098]
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地震资料的水平线性干扰去除方法。
[0099]
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述地震资料的水平线性干扰去除方法的计算机程序。
[0100]
本发明实施例中，从地震资料的叠加剖面中，获取预设的水平线性干扰待处理区域的时窗；计算水平线性干扰待处理区域的时窗的应用时差值；应用时差值用于表征时窗的最大值和最小值的差值；对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理，得到处理后的地震资料；根据沿预设地层同相轴方向拾取的平滑曲线，对处理后的地震资料进行正向校正处理，得到地震正变换资料；正向校正处理用于对地震资料进行基于地层倾角的正向变换；对地震正变换资料进行去线性干扰处理；对去线性干扰处理后的地震正变换资料进行反向校正处理，得到地震反变换资料；反向校正处理用于对地震资料进行基于地层倾角的反向变换；对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理，得到消除水平线性干扰后的地震资料，与现有技术中地震资料的水平线性干扰很难通过去线性干扰处理加以去除的技术方案相比，通过时移处理、地震道长度延长处理和正向校正处理，可将地震资料的水平线性干扰
改变成其他角度，可使去线性干扰处理间接性地发挥作用，极大地提高了去线性干扰处理的可行性，从而改善了地震资料去除水平线性干扰的处理效果。
[0101]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0102]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0103]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0104]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0105]
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。