CN112379426B - 地震资料的水平线性干扰去除方法及装置 - Google Patents
地震资料的水平线性干扰去除方法及装置 Download PDFInfo
- Publication number
- CN112379426B CN112379426B CN202011192300.6A CN202011192300A CN112379426B CN 112379426 B CN112379426 B CN 112379426B CN 202011192300 A CN202011192300 A CN 202011192300A CN 112379426 B CN112379426 B CN 112379426B
- Authority
- CN
- China
- Prior art keywords
- seismic
- data
- processing
- difference value
- seismic data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000006243 chemical reaction Methods 0.000 claims abstract description 76
- 230000009466 transformation Effects 0.000 claims abstract description 50
- 238000011282 treatment Methods 0.000 claims abstract description 46
- 230000009467 reduction Effects 0.000 claims abstract description 19
- 230000008569 process Effects 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 16
- 230000008030 elimination Effects 0.000 claims description 15
- 238000003379 elimination reaction Methods 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000011946 reduction process Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 9
- 238000010586 diagram Methods 0.000 description 11
- 230000015572 biosynthetic process Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000009499 grossing Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种地震资料的水平线性干扰去除方法及装置,该方法包括:获取预设的水平线性干扰待处理区域的时窗;计算水平线性干扰待处理区域的时窗的应用时差值;对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理;根据沿预设地层同相轴方向拾取的平滑曲线,对得到的处理后的地震资料进行正向校正处理;之后对得到的地震正变换资料进行去线性干扰处理;对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料。本发明可改善地震资料去除水平线性干扰的处理效果。
Description
技术领域
本发明涉及地震勘探数据处理技术领域,尤其涉及地震资料的水平线性干扰去除方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
地震资料在处理过程中,由于种种原因会在叠加剖面上存在一些水平线性干扰,如水平线性噪声。这些水平线性干扰的存在严重影响了地震资料的最终成像效果,进而影响最终的地层构造解释。而现阶段,消除这些水平线性干扰通常是基于速度或者基于时差进行识别进而加以去除。但地震资料的水平同相轴的速度是接近无穷大的,对应的时差接近于零。此外,地震资料中的水平线性干扰很难通过去线性干扰处理加以去除,如频率波数域变换(FK变换)等方法。
总之,现有技术下地震资料去除水平线性干扰的处理效果不佳,地震资料中的水平线性干扰难以消除。
发明内容
本发明实施例提供一种地震资料的水平线性干扰去除方法,用以改善地震资料去除水平线性干扰的处理效果,该方法包括:
从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;
计算水平线性干扰待处理区域的时窗的应用时差值;所述应用时差值用于表征时窗的最大值和最小值的差值;
对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;
根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;所述正向校正处理用于对地震资料进行基于地层倾角的正向变换;
对地震正变换资料进行去线性干扰处理;
对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;所述反向校正处理用于对地震资料进行基于地层倾角的反向变换;
对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料。
本发明实施例还提供一种地震资料的水平线性干扰去除装置,用以改善地震资料去除水平线性干扰的处理效果,该装置包括:
时窗获取模块,用于从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;
应用时差值计算模块,用于计算水平线性干扰待处理区域的时窗的应用时差值;所述应用时差值用于表征时窗的最大值和最小值的差值;
时移处理和地震道长度延长处理模块,用于对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;
正向校正处理模块,用于根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;所述正向校正处理用于对地震资料进行基于地层倾角的正向变换;
去线性干扰处理模块,用于对地震正变换资料进行去线性干扰处理;
反向校正处理模块,用于对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;所述反向校正处理用于对地震资料进行基于地层倾角的反向变换;
反向时移处理和地震道长度缩减处理模块,用于对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述地震资料的水平线性干扰去除方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述地震资料的水平线性干扰去除方法的计算机程序。
本发明实施例中,从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;计算水平线性干扰待处理区域的时窗的应用时差值;所述应用时差值用于表征时窗的最大值和最小值的差值;对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;所述正向校正处理用于对地震资料进行基于地层倾角的正向变换;对地震正变换资料进行去线性干扰处理;对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;所述反向校正处理用于对地震资料进行基于地层倾角的反向变换;对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料,与现有技术中地震资料的水平线性干扰很难通过去线性干扰处理加以去除的技术方案相比,通过时移处理、地震道长度延长处理和正向校正处理,可将地震资料的水平线性干扰改变成其他角度,可使去线性干扰处理间接性地发挥作用,极大地提高了去线性干扰处理的可行性,从而改善了地震资料去除水平线性干扰的处理效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中一种地震资料的水平线性干扰去除方法的流程示意图;
图2为本发明实施例中一种地震资料的水平线性干扰去除方法的具体示例图;
图3为本发明实施例中一种地震资料的水平线性干扰去除方法的具体示例图;
图4为本发明实施例中一种地震资料的水平线性干扰去除方法的具体示例图;
图5为本发明实施例中一种地震资料的水平线性干扰去除装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
图1为本发明实施例中一种地震资料的水平线性干扰去除方法的流程示意图,如图1所示,本发明实施例提供的地震资料的水平线性干扰去除方法,可以包括:
步骤101:从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;
步骤102:计算水平线性干扰待处理区域的时窗的应用时差值;应用时差值用于表征时窗的最大值和最小值的差值;
步骤103:对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;
步骤104:根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;正向校正处理用于对地震资料进行基于地层倾角的正向变换;
步骤105:对地震正变换资料进行去线性干扰处理;
步骤106:对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;反向校正处理用于对地震资料进行基于地层倾角的反向变换;
步骤107:对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料。
具体实施时,首先从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗。
实施例中,可通过直接在叠加剖面上选取要进行水平线性干扰待处理的区域来获取水平线性干扰待处理区域的时窗。
在上述实施例中,可通过从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗,有助于在后续步骤中对水平线性干扰待处理区域进行处理。
具体实施时,在从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗后,计算水平线性干扰待处理区域的时窗的应用时差值;应用时差值用于表征时窗的最大值和最小值的差值。
实施例中,计算水平线性干扰待处理区域的时窗的应用时差值,可以包括:获取水平线性干扰待处理区域的时窗的最大值和最小值;根据获取的时窗的最大值和最小值,计算时窗的应用时差值。
在上述实施例中,通过确定水平线性干扰待处理区域的所涉及到时窗的最小时间和最大时间,计算时窗的应用时差值,,有助于在后续步骤中对水平线性干扰待处理区域进行处理。
具体实施时,在计算水平线性干扰待处理区域的时窗的应用时差值;应用时差值用于表征时窗的最大值和最小值的差值后,对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料。
实施例中,对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料,可以包括:对地震资料中地震道以应用时差值的第一预设倍数,进行向下时移处理;对向下时移处理后的地震道以应用时差值的第二预设倍数,进行地震道长度延长处理,得到处理后的地震资料。
在上述实施例中,第一预设倍数可为一倍;第二预设倍数可为两倍。地震资料中地震道以应用时差值的第一预设倍数,进行向下时移处理;对向下时移处理后的地震道以应用时差值的第二预设倍数,进行地震道长度延长处理,可以包括:将地震道整体向下时移一倍应用时差值长度,并且将道长度延长两倍应用时差值长度。具体实施时,在对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料后,根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;正向校正处理用于对地震资料进行基于地层倾角的正向变换。
实施例中,根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料,可以包括:对处理后的地震资料,沿预设地层同相轴方向拾取平滑曲线,计算该平滑曲线经每一地震道的拾取时长与该平滑曲线平均值的差值;计算该差值的平均值;根据该差值的平均值,对处理后的地震资料进行正向校正处理,得到地震正变换资料。
在上述实施例中,沿预设地层同相轴方向拾取平滑曲线,可以包括:沿某一地层倾角方向拾取平滑的曲线,其他不需要进行水平线性干扰去除处理的部分进行适当平滑外推,直至全部拾取完成;求取所拾取平滑曲线的平均值;求取平滑曲线上每一个值和该平滑曲线平均值的平均时差;使用该平均时差对处理后的地震资料进行正向校正处理。由此实现了对地震数据的地层倾角变换的正向变换。在地层倾角正向变换之后的地震数据上做相应的地震资料处理,可将地震资料的水平线性干扰改变成其他角度,可使去线性干扰处理间接性地发挥作用,极大地提高了去线性干扰处理的可行性,为FK变换等去线性处理奠定基础。
具体实施时,在根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料后,对地震正变换资料进行去线性干扰处理。
在上述实施例中,可通过去线性干扰处理对地震正变换资料中水平线性干扰加以去除,如通过频率波数域变换(FK变换)等方法。在地震正变换资料上可进行去线性等处理去除原先水平,经正向校正处理后倾斜的水平线性干扰。
具体实施时,在对地震正变换资料进行去线性干扰处理后,对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;反向校正处理用于对地震资料进行基于地层倾角的反向变换。
实施例中,对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料,可以包括:根据该差值的平均值,对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料。
在上述实施例中,使用该差值的平均值对去线性干扰处理后的地震正变换资料进行时差的反向应用。
具体实施时,在对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料后,对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料。
实施例中,对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,可以包括:对地震反变换资料中地震道以应用时差值的第一预设倍数,进行向上时移处理;对向上时移处理后的地震道以应用时差值的第二预设倍数,进行地震道长度缩短处理,得到处理后的地震资料。
在上述实施例中,将地震反变换资料整体向上时移一倍应用时差值长度,并且将道长度缩短至原先道长度,由此实现了地震数据的地层倾角变换的反向变换。通过进行地层倾角变换的反变换,可使地震反变换资料回到原始的地震数据的构造状态,经过FK变换等处理之后,水平的线性干扰,如噪声,可得到很好的消除。
下面给出一个具体实施例,来说明本发明的方法的具体应用。该具体实施例,可以包括如下步骤:
1、如图2所示,图2显示了地震资料的叠加剖面;首先在地震资料的叠加剖面上确定水平线性干扰待处理区域的时窗;图2中纵坐标是时间轴,英文是time;横坐标是地震资料叠加剖面的“共中心点”顺序号,英文简写是CMP或者CDP,全拼为common midpoint;
2、在地震资料的叠加剖面中水平线性干扰待处理区域中,查看水平线性干扰待处理区域的时窗所涉及到的最小时间Tmin和最大时间Tmax;
3、确定最小时间和最大时间的时差值为应用时差值T,其中,T=Tmax-Tmin;
4、将地震资料中地震道的叠加剖面数据整体向下时移一倍应用时差值T,并且将地震道长度延长两倍应用时差值形成处理后的地震资料;在该处理后的地震资料中道长度Ttrace=Ttrace0+2T,其中,Ttrace0为地震资料的原始地震道长度;Ttrace为处理后的地震资料的地震道长度;
5、沿处理后的地震资料的水平线性干扰待处理区域的某一倾斜同相轴拾取平滑的曲线,其他不需要处理部分进行适当平滑外推,直至全部拾取完成;
6、求取所拾取平滑曲线的平均值Tmean,其中Ti为曲线上每一道的拾取时间;n为曲线上的地震道的总道数,
7、求取平滑曲线上每一地震道的时间数值Ti和平均值Tmean的平均差值ΔTi;
8、使用该差值ΔTi对地震数据每一地震道进行校正,这样处理后的地震资料就会沿着所拾取的曲线进行拉平,进而形成一个地震正变换资料,如图3所示;
9、经过以上各步处理之后就实现了地震资料地层倾角变换的正变换;
10、在地震正变换资料上进行FK等处理去除原先为水平状态,现在为倾斜状态的水平线性干扰,形成去线性干扰处理后的地震正变换资料;
11、将去线性干扰处理后的地震正变换资料使用之前的平均时差ΔTi进行反应用,即反向校正处理,形成地震反变换资料;
12、将地震反变换资料分别进行反向时移处理和地震道长度缩减处理,如整体向上时移一倍应用时差值T,并且将道长度缩短至原始地震道长度;
13、经过以上两步就实现了地层倾角变换的反变换,得到消除水平线性干扰后的地震资料,如图4所示。图4经过地层倾角变换正变换之后的FK去除线性,有效消除了地震资料的水平噪声。
当然,可以理解的是,上述详细流程还可以有其他变化例,相关变化例均应落入本发明的保护范围。
本发明实施例中,从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;计算水平线性干扰待处理区域的时窗的应用时差值;应用时差值用于表征时窗的最大值和最小值的差值;对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;正向校正处理用于对地震资料进行基于地层倾角的正向变换;对地震正变换资料进行去线性干扰处理;对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;反向校正处理用于对地震资料进行基于地层倾角的反向变换;对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料,与现有技术中地震资料的水平线性干扰很难通过去线性干扰处理加以去除的技术方案相比,通过时移处理、地震道长度延长处理和正向校正处理,可将地震资料的水平线性干扰改变成其他角度,可使去线性干扰处理间接性地发挥作用,极大地提高了去线性干扰处理的可行性,从而改善了地震资料去除水平线性干扰的处理效果。
如上所示,本发明实施例本方法能够改变地震资料的地层倾角分布特征,从而改变地层倾角的角度,从而使FK去除线性得到好的处理效果。
本发明实施例中还提供了一种的水平线性干扰去除装置,如下面的实施例。由于该装置解决问题的原理与的水平线性干扰去除方法相似,因此该装置的实施可以参见的水平线性干扰去除方法的实施,重复之处不再赘述。
图5为本发明实施例中一种地震资料的水平线性干扰去除装置的结构示意图,如图5所示,本发明实施例提供的地震资料的水平线性干扰去除装置,可以包括:
时窗获取模块01,用于从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;
应用时差值计算模块02,用于计算水平线性干扰待处理区域的时窗的应用时差值;应用时差值用于表征时窗的最大值和最小值的差值;
时移处理和地震道长度延长处理模块03,用于对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;
正向校正处理模块04,用于根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;正向校正处理用于对地震资料进行基于地层倾角的正向变换;
去线性干扰处理模块05,用于对地震正变换资料进行去线性干扰处理;
反向校正处理模块06,用于对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;反向校正处理用于对地震资料进行基于地层倾角的反向变换;
反向时移处理和地震道长度缩减处理模块07,用于对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料。
在一个实施例中,应用时差值计算模块,具体用于:获取水平线性干扰待处理区域的时窗的最大值和最小值;
根据获取的时窗的最大值和最小值,计算时窗的应用时差值。
在一个实施例中,时移处理和地震道长度延长处理模块,具体用于:
对地震资料中地震道以应用时差值的第一预设倍数,进行向下时移处理;
对向下时移处理后的地震道以应用时差值的第二预设倍数,进行地震道长度延长处理,得到处理后的地震资料。
在一个实施例中,反向时移处理和地震道长度缩减处理模块,具体用于:
对地震反变换资料以应用时差值的第一预设倍数,进行向上时移处理;
对向上时移处理后的地震反变换资料以应用时差值的第二预设倍数,进行地震道长度缩减处理,得到消除水平线性干扰后的地震资料。
在一个实施例中,正向校正处理模块,具体用于:
对处理后的地震资料,沿预设地层同相轴方向拾取平滑曲线,计算该平滑曲线经每一地震道的拾取时长与该平滑曲线平均值的差值;
计算该差值的平均值;
根据该差值的平均值,对处理后的地震资料进行正向校正处理,得到地震正变换资料。
在一个实施例中,反向校正处理模块,具体用于:
根据该该差值的平均值,对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述地震资料的水平线性干扰去除方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述地震资料的水平线性干扰去除方法的计算机程序。
本发明实施例中,从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;计算水平线性干扰待处理区域的时窗的应用时差值;应用时差值用于表征时窗的最大值和最小值的差值;对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;正向校正处理用于对地震资料进行基于地层倾角的正向变换;对地震正变换资料进行去线性干扰处理;对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;反向校正处理用于对地震资料进行基于地层倾角的反向变换;对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料,与现有技术中地震资料的水平线性干扰很难通过去线性干扰处理加以去除的技术方案相比,通过时移处理、地震道长度延长处理和正向校正处理,可将地震资料的水平线性干扰改变成其他角度,可使去线性干扰处理间接性地发挥作用,极大地提高了去线性干扰处理的可行性,从而改善了地震资料去除水平线性干扰的处理效果。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种地震资料的水平线性干扰去除方法,其特征在于,包括:
从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;
计算水平线性干扰待处理区域的时窗的应用时差值;所述应用时差值用于表征时窗的最大值和最小值的差值;
对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;
根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;所述正向校正处理用于对地震资料进行基于地层倾角的正向变换;
对地震正变换资料进行去线性干扰处理;
对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;所述反向校正处理用于对地震资料进行基于地层倾角的反向变换;
对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料;
根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料,包括:对处理后的地震资料,沿预设地层同相轴方向拾取平滑曲线,计算该平滑曲线经每一地震道的拾取时长与该平滑曲线平均值的差值;计算该差值的平均值;根据该差值的平均值,对处理后的地震资料进行正向校正处理,得到地震正变换资料;
对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料,包括:根据该差值的平均值,对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料。
2.如权利要求1所述的方法,其特征在于,对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料,包括:
对地震资料中地震道以应用时差值的第一预设倍数,进行向下时移处理;
对向下时移处理后的地震道以应用时差值的第二预设倍数,进行地震道长度延长处理,得到处理后的地震资料。
3.如权利要求2所述的方法,其特征在于,对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料,包括:
对地震反变换资料以应用时差值的第一预设倍数,进行向上时移处理;
对向上时移处理后的地震反变换资料以应用时差值的第二预设倍数,进行地震道长度缩减处理,得到消除水平线性干扰后的地震资料。
4.一种地震资料的水平线性干扰去除装置,其特征在于,包括:
时窗获取模块,用于从地震资料的叠加剖面中,获取预设的水平线性干扰待处理区域的时窗;
应用时差值计算模块,用于计算水平线性干扰待处理区域的时窗的应用时差值;所述应用时差值用于表征时窗的最大值和最小值的差值;
时移处理和地震道长度延长处理模块,用于对地震资料中地震道以应用时差值的预设倍数分别进行时移处理和地震道长度延长处理,得到处理后的地震资料;
正向校正处理模块,用于根据沿预设地层同相轴方向拾取的平滑曲线,对处理后的地震资料进行正向校正处理,得到地震正变换资料;所述正向校正处理用于对地震资料进行基于地层倾角的正向变换;
去线性干扰处理模块,用于对地震正变换资料进行去线性干扰处理;
反向校正处理模块,用于对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料;所述反向校正处理用于对地震资料进行基于地层倾角的反向变换;
反向时移处理和地震道长度缩减处理模块,用于对地震反变换资料以应用时差值的预设倍数分别进行反向时移处理和地震道长度缩减处理,得到消除水平线性干扰后的地震资料;
正向校正处理模块,具体用于:
对处理后的地震资料,沿预设地层同相轴方向拾取平滑曲线,计算该平滑曲线经每一地震道的拾取时长与该平滑曲线平均值的差值;计算该差值的平均值;根据该差值的平均值,对处理后的地震资料进行正向校正处理,得到地震正变换资料;
反向校正处理模块,具体用于:根据该差值的平均值,对去线性干扰处理后的地震正变换资料进行反向校正处理,得到地震反变换资料。
5.如权利要求4所述的装置,其特征在于,时移处理和地震道长度延长处理模块,具体用于:
对地震资料中地震道以应用时差值的第一预设倍数,进行向下时移处理;
对向下时移处理后的地震道以应用时差值的第二预设倍数,进行地震道长度延长处理,得到处理后的地震资料。
6.如权利要求5所述的装置,其特征在于,反向时移处理和地震道长度缩减处理模块,具体用于:
对地震反变换资料以应用时差值的第一预设倍数,进行向上时移处理;
对向上时移处理后的地震反变换资料以应用时差值的第二预设倍数,进行地震道长度缩减处理,得到消除水平线性干扰后的地震资料。
7.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至3任一所述方法。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至3任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011192300.6A CN112379426B (zh) | 2020-10-30 | 2020-10-30 | 地震资料的水平线性干扰去除方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011192300.6A CN112379426B (zh) | 2020-10-30 | 2020-10-30 | 地震资料的水平线性干扰去除方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112379426A CN112379426A (zh) | 2021-02-19 |
CN112379426B true CN112379426B (zh) | 2023-09-26 |
Family
ID=74576912
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011192300.6A Active CN112379426B (zh) | 2020-10-30 | 2020-10-30 | 地震资料的水平线性干扰去除方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112379426B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
CN104977615A (zh) * | 2015-07-01 | 2015-10-14 | 中国海洋石油总公司 | 一种基于模型统计拾取的深水obc资料多次波压制方法 |
CN105911585A (zh) * | 2016-07-05 | 2016-08-31 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震记录规则干扰波的提取方法及装置 |
CN107193045A (zh) * | 2017-07-07 | 2017-09-22 | 中国石油天然气集团公司 | 一种地震数据处理方法及装置 |
CN111781644A (zh) * | 2019-04-03 | 2020-10-16 | 中国石油天然气股份有限公司 | 浅中地层地震数据的线性干扰衰减方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140078860A1 (en) * | 2012-09-19 | 2014-03-20 | Cgg Services Sa | Interference noise attenuation method and apparatus |
US10444393B2 (en) * | 2015-03-27 | 2019-10-15 | Halliburton Energy Services, Inc. | Coherent noise attenuation using statistical methods |
-
2020
- 2020-10-30 CN CN202011192300.6A patent/CN112379426B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
CN104977615A (zh) * | 2015-07-01 | 2015-10-14 | 中国海洋石油总公司 | 一种基于模型统计拾取的深水obc资料多次波压制方法 |
CN105911585A (zh) * | 2016-07-05 | 2016-08-31 | 中国石油集团东方地球物理勘探有限责任公司 | 一种地震记录规则干扰波的提取方法及装置 |
CN107193045A (zh) * | 2017-07-07 | 2017-09-22 | 中国石油天然气集团公司 | 一种地震数据处理方法及装置 |
CN111781644A (zh) * | 2019-04-03 | 2020-10-16 | 中国石油天然气股份有限公司 | 浅中地层地震数据的线性干扰衰减方法及装置 |
Non-Patent Citations (6)
Title |
---|
Nawal A. Al-Ridha等.Seismic Data Processing of Subba Oil Field in South Iraq . Iraqi Journal of Science.2018,全文. * |
付燕,赵荣椿.应用时空变-倾角扫描叠加KL变换提高地震资料信噪比.物探化探计算技术.2001,(第04期),全文. * |
崔三元等.起伏地表条件下近地表低降速带对地震资料成像影响研究与应用.SPG/SEG南京 2020 年国际地球物理会议 .2020,全文. * |
崔辉霞等.GeoEast系统海底电缆双检合成处理技术.石油地球物理勘探.2014,全文. * |
李振春等.地震数据规则化重构方法策略.中国石油大学学报( 自然科学版).2018,全文. * |
苑益军 ; 徐林 ; .叠前线性干扰自动追踪SVD压制方法.现代地质.2007,(第04期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN112379426A (zh) | 2021-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109785854B (zh) | 一种经验模态分解和小波阈值去噪相结合的语音增强方法 | |
CN110598166B (zh) | 一种自适应确定小波分层级数的小波去噪方法 | |
CN111311510B (zh) | 图像中条纹噪声的去除方法及装置 | |
CN102386889B (zh) | 基线漂移去除方法、装置及中值滤波器 | |
CN110619265A (zh) | 球磨机筒体振动信号联合去噪方法、装置及存储介质 | |
CN112379426B (zh) | 地震资料的水平线性干扰去除方法及装置 | |
CN109684937B (zh) | 一种基于fft及数学形态法的信号去噪方法及装置 | |
CN117318671B (zh) | 一种基于快速傅里叶变换的自适应滤波方法 | |
JP3108317B2 (ja) | 信号処理装置および信号処理方法 | |
CN102074241A (zh) | 一种通过快速声音波形修复实现声音还原的方法 | |
CN110161565B (zh) | 一种地震数据重建方法 | |
CN110413945B (zh) | 基于分数阶短时傅里叶变换的线性调频信号相位恢复方法 | |
CN113297195B (zh) | 一种时间序列异常检测方法、装置及设备 | |
CN111781644B (zh) | 浅中地层地震数据的线性干扰衰减方法及装置 | |
CN111064687B (zh) | 一种频率和相位补偿的方法及装置 | |
CN111781642A (zh) | 地震数据层间多次波衰减方法及装置 | |
CN111766444A (zh) | 基于综合算法的多分量线性调频信号参数估计方法及系统 | |
CN103971319A (zh) | 一种由计算机执行的图像处理方法 | |
CN114118155B (zh) | 一种火箭着陆冲击信号的处理方法、装置及计算机设备 | |
CN113138419B (zh) | 提取下行子波和衰减参数的方法、装置 | |
CN116595324B (zh) | 信号瞬态冲击始点的提取方法 | |
CN110737019B (zh) | 提高速度谱精度的处理方法及系统 | |
CN111540372B (zh) | 一种多麦克阵列降噪处理的方法及装置 | |
CN112526604B (zh) | 基于目的层谱分析的自适应低频补偿方法及系统 | |
CN113156495B (zh) | 网格层析反演反射点确定方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |