CN111562616A - 地震数据散射噪音压制方法及装置 - Google Patents
地震数据散射噪音压制方法及装置 Download PDFInfo
- Publication number
- CN111562616A CN111562616A CN201910113832.7A CN201910113832A CN111562616A CN 111562616 A CN111562616 A CN 111562616A CN 201910113832 A CN201910113832 A CN 201910113832A CN 111562616 A CN111562616 A CN 111562616A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- noise
- dimensional
- transform
- scattering noise
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000001629 suppression Effects 0.000 title claims description 47
- 229910052704 radon Inorganic materials 0.000 claims abstract description 194
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 191
- 230000009466 transformation Effects 0.000 claims abstract description 63
- 238000012937 correction Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 13
- 230000000116 mitigating effect Effects 0.000 claims 2
- 238000005070 sampling Methods 0.000 abstract description 11
- 238000012545 processing Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 16
- 230000002238 attenuated effect Effects 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 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
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/46—Radon transform
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
技术领域
本发明涉及地震数据处理技术领域，尤其涉及一种地震数据散射噪音压制方法及装置。
背景技术
散射波广泛存在于沙漠地表等有地形横向变化和速度纵横向变化的地震数据中。在地震数据处理过程中需要将散射波作为噪声进行处理。由于散射噪音的波形为双曲线且频率和能量与反射信号差异较小，因此难以被有效衰减。未被有效去除的散射波在进行完常规处理生成叠加剖面后往往以随机干扰和近似线性干扰的形式存在于叠加剖面上。
现有技术中，通常采用线性干扰压制方法对地震数据中的散射波进行压制，但是，由于散射噪音的形态在空间上为双曲球体，通常形成一次波的近地表条件还易于产生多次波，因此散射噪音中除了包含一次波，还常会伴随多次波，而采用线性干扰压制方法无法有效压制伴随产生的多次波。
拉东变换是指沿着某一积分路径将信号叠加，然后再沿相同的积分路径将叠加后的信号进行恢复，其中沿某一积分路径将信号叠加为拉东正变换，沿相同的积分路径将叠加后的信号进行恢复为拉东反变换。如果原始信号上带有散射噪音(一次波或多次波)，叠加后散射噪音能够相互抵消，在拉东反变换后相互抵消的散射噪音不再反变换回原始信号，因此对地震数据进行拉东变换可以有效压制散射噪音中的一次波和多次波。三维拉东变换相较于二维拉东变换，精度高，压制地震数据中散射噪音的效果好。
在应用三维拉东变换时，叠加过程中拉东域的取值范围越大则对原始空间域的数据描述越准确，然而由于取值范围不能无限制的增大，总会存在一个截断值，该截断值以外的数据无法通过三维拉东变换转换到拉东域，也即三维拉东变换采样不足，引起信息泄露问题。
发明内容
本发明实施例提供一种地震数据散射噪音压制方法，用以压制地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题，该方法包括：
获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音；
对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音；
从第一地震数据中去除第一散射噪音和第二散射噪音。
本发明实施例提供一种地震数据散射噪音压制装置，用以压制地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题，该装置包括：
第一噪音确定模块，用于获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音；
第二噪音确定模块，用于对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音；
噪音压制模块，用于从第一地震数据中去除第一散射噪音和第二散射噪音。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如下方法：
获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音；
对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音；
从第一地震数据中去除第一散射噪音和第二散射噪音。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行如下方法的计算机程序：
获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音；
对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音；
从第一地震数据中去除第一散射噪音和第二散射噪音。
本发明实施例通过获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音，然后对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音，最后从第一地震数据中去除第一散射噪音和第二散射噪音。本发明实施例通过两组三维拉东正变换与三维拉东反变换，完整提取出地震数据中存在的全部散射噪声，然后从地震数据中去除全部散射噪声，获得完整有效的地震信息，有效压制了地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中地震数据散射噪音压制方法示意图；
图2为本发明实施例中对第一地震数据进行线性干扰压制处理和强能量干扰压制处理后变换至十字排列域的结果示意图；
图3为本发明实施例中对线性干扰压制处理和强能量干扰压制处理结果进行部分动校正后的结果示意图；
图4为本发明实施例中第二地震数据示意图；
图5为本发明实施例中第一散射噪音示意图；
图6为本发明实施例中对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换后的结果示意图；
图7为本发明实施例中第二散射噪音示意图；
图8为本发明实施例中从第一地震数据中去除第一散射噪音和第二散射噪音后的结果示意图；
图9为本发明实施例中对去除第一散射噪音和第二散射噪音之后的第一地震数据进行反动校正后的结果示意图；
图10为本发明实施例中所衰减的散射噪音单炮剖面示意图；
图11为本发明实施例中所衰减的散射噪音时间切片示意图；
图12为本发明实施例中散射噪音压制前的叠加剖面示意图；
图13为本发明实施例中散射噪音压制后的叠加剖面示意图；
图14为本发明实施例中散射噪音在叠加剖面上的形态示意图；
图15为本发明实施例中地震数据散射噪音压制装置结构图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
为了压制地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题，本发明实施例提供一种地震数据散射噪音压制方法，如图1所示，该方法可以包括：
步骤101、获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音；
步骤102、对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音；
步骤103、从第一地震数据中去除第一散射噪音和第二散射噪音。
由图1所示可以得知，本发明实施例通过获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音，然后对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音，最后从第一地震数据中去除第一散射噪音和第二散射噪音。本发明实施例通过两组三维拉东正变换与三维拉东反变换，完整提取出地震数据中存在的全部散射噪声，然后从地震数据中去除全部散射噪声，获得完整有效的地震信息，有效压制了地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题。
具体实施时，获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音。发明人发现，现有技术中通常采用线性干扰压制方法对地震数据中的散射波进行压制，但是，由于散射噪音的形态在空间上为双曲球体，通常形成一次波的近地表条件还易于产生多次波，因此散射噪音中除了包含一次波，还常会伴随多次波，而采用线性干扰压制方法无法有效压制伴随产生的多次波。三维拉东变换是指沿着某一积分路径将信号叠加，然后再沿相同的积分路径将叠加后的信号进行恢复，如果原始信号上带有散射噪音(一次波或多次波)，叠加后散射噪音能够相互抵消，在反变换后相互抵消的散射噪音不再反变换回原始信号。基于此，本发明实施例对地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，可以有效提取散射噪音中的一次波和多次波。
实施例中，对第一地震数据进行三维拉东正变换之前，将第一地震数据变换至十字排列域。发明人发现，线性干扰在十字排列域会变成圆锥体的形态，便于去噪。基于此，本发明实施例对第一地震数据进行三维拉东正变换之前，通过对第一地震数据按照炮线和检波线道头进行重新排序将数据由炮域转换到十字排列域，将线性干扰转变为圆锥体形态，便于去噪。
实施例中，对第一地震数据进行三维拉东正变换之前，进行第一地震数据线性干扰压制处理和强能量干扰压制处理。本发明实施例在对第一地震数据进行三维拉东正变换之前，进行第一地震数据线性干扰压制处理和强能量干扰压制处理，可以有效衰减部分线性和强能量干扰。线性干扰压制处理和强能量干扰压制处理为现有技术，本领域技术人员能够通过查阅资料了解线性干扰压制处理和强能量干扰压制处理，本发明不再进行具体说明。
实施例中，按如下方法确定第一散射噪音：对第一地震数据和第二地震数据作差，将作差的结果作为第一散射噪音。
实施例中，对第一地震数据进行三维拉东正变换之前，对第一地震数据进行部分动校正。部分动校正是对地震数据利用速度分析得到的均方根速度进行动校正，发明人发现，如果采用的均方根速度合理则动校正后的地震同相轴为水平的，如果采用的均方根速度比合理的速度小或者大，则同相轴近于水平而无法完全水平，这种动校正为部分动校正，完全动校正在将同相轴校正为水平状态的同时会失去部分浅层地震数据，而部分动校正则可以避免失去这部分地震数据。需要说明的是，速度分析指对地震数据按照不同的速度进行扫描叠加，由叠加结果形成速度谱，人工在速度谱上拾取速度。速度分析为现有技术，本领域技术人员能够通过查阅资料了解速度分析，本发明不再进行具体说明。
实施例中，三维拉东正变换为三维抛物拉东正变换，所述三维拉东反变换为三维抛物拉东反变换。发明人发现，地震资料的同相轴为双曲线的形式，因此需要采用三维双曲拉东变换对地震数据进行处理，但是三维双曲拉东正变换和三维双曲拉东反变换公式复杂，不易实际应用。基于此，本发明实施例通过对地震数据进行部分动校正，将双曲线形式的同相轴转换为抛物线形式，然后采用计算精度和效率都非常高的三维抛物拉东正变换和三维抛物拉东反变换对地震数据进行处理，有效简化了运算过程，提高了计算精度和效率。
具体实施时，对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音。发明人发现，在三维拉东变换的叠加过程中，拉东域的取值范围越大则对原始空间域的数据描述越准确，然而由于取值范围不能无限制的增大，总会存在一个截断值，该截断值以外的数据无法通过三维拉东变换转换到拉东域，也即三维拉东变换采样不足，引起信息泄露问题。基于此，本发明实施例在对地震数据进行一组三维拉东正变换与三维拉东反变换提取出一部分散射噪音之后，继续对这部分散射噪音进行三维拉东正变换与三维拉东反变换，提取地震数据中残留的散射噪音，通过两组三维拉东正变换与三维拉东反变换，完整提取出地震数据中存在的全部散射噪声，然后从地震数据中去除全部散射噪声，即可获得完整有效的地震信息，有效压制了地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题。
实施例中，按如下方法确定第二散射噪音：对第一散射噪音和三维拉东反变换结果作差，将作差的结果作为第二散射噪音。
具体实施时，从第一地震数据中去除第一散射噪音和第二散射噪音。通过两组三维拉东正变换与三维拉东反变换，完整提取出地震数据中存在的全部散射噪声，即第一散射噪音和第二散射噪音，进而从地震数据中去除全部散射噪声，即可获得完整有效的地震信息，有效压制了地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题。
实施例中，由于对第一地震数据进行三维拉东正变换之前，对第一地震数据进行了部分动校正，因此在第一地震数据去除第一散射噪音和第二散射噪音之后，对结果数据进行反动校正。
下面给出一个具体实施例，说明本发明实施例中地震数据散射噪音压制方法的具体应用。首先进行线性干扰压制处理和强能量干扰压制处理，然年将结果按照炮线和检波线进行排序得到十字排列域的道集数据，如图2所示，其中，横坐标为空间坐标，用来衡量空间相对位置关系，纵坐标为时间(秒)，灰度值为地震数据的振幅值，表示地下截至反射系数的强弱。对十字排列域道集进行速度分析得到均方根速度，对十字排列域道集按速度分析得到的均方根速度乘以在0-1之间的比例系数后进行部分动校正，比例系数的选择原则为动校正后的同相轴为近似抛物线且浅层初至波不会因为动校拉伸被破坏。图3为对线性干扰压制处理和强能量干扰压制处理结果进行部分动校正后的结果示意图。对部分动校正后的数据进行一组三维抛物拉东正变换和三维抛物拉东反变换，得到第二地震数据，如图4所示。由部分动校正后的数据与第二地震数据相减得到两者残差，也即第一散射噪音，如图5所示。对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到的结果如图6所示。将三维拉东反变换的结果与第一散射噪音相减得到两者残差，也即第二散射噪音，如图7所示。从第一地震数据中去除第一散射噪音和第二散射噪音，得到的结果如图8所示。最后，对去除第一散射噪音和第二散射噪音之后的第一地震数据进行反动校正，结果如图9所示。为验证本发明实施例的效果，下面以图4所示的道集为例，采用地震数据散射噪音压制方法，可得到如图9所示的道集。图10为图4与图9所示数据的差，为衰减掉的散射噪音，其中多为高频的双曲干扰，图11为对衰减掉的散射噪音进行的时间切片，切片时间为1700ms，从切片可以看出散射噪音在切片上呈现圆环状形态，其中，横纵坐标表示地震数据纵横向采集测线编号，实际和大地坐标存在一一对应关系，用来描述空间位置，图形灰度值表示三维地震数据在时间上的振幅切片，振幅表示地下介质反射系数的强弱。图10和图11均只见到衰减掉的散射噪音，未见有效反射信号被损伤。图12为散射噪音压制前叠加剖面，其中，横坐标为空间坐标，用来衡量空间相对位置关系，纵坐标为时间(秒)，灰度值为地震数据的振幅值，表示地下截至反射系数的强弱，剖面上信噪比低，存在大量高频线性和随机干扰，这些干扰实际是散射噪音在叠加后的变现形式，图13为对叠加前的道集数据进行散射噪音压制后再进行叠加得到的叠加剖面，与图12相比较反射信号信噪比得到明显提高，图14为衰减的散射噪音进行叠加后得到的叠加剖面，可见散射噪音再叠加剖面上的形态为双曲线，由于有过多的散射噪音叠和在一起，部分表现出随机干扰的形态。
基于同一发明构思，本发明实施例还提供了一种地震数据散射噪音压制装置，如下面的实施例所述。由于这些解决问题的原理与地震数据散射噪音压制方法相似，因此装置的实施可以参见方法的实施，重复之处不再赘述。
图15为本发明实施例中地震数据散射噪音压制装置的结构图，如图15所示，该装置包括：
第一噪音确定模块151，用于获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音；
第二噪音确定模块152，用于对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音；
噪音压制模块153，用于从第一地震数据中去除第一散射噪音和第二散射噪音。
一个实施例中，第一噪音确定模块151进一步用于，对第一地震数据进行三维拉东正变换之前，将第一地震数据变换至十字排列域。
一个实施例中，第一噪音确定模块151进一步用于，对第一地震数据进行三维拉东正变换之前，进行第一地震数据线性干扰压制处理和强能量干扰压制处理。
一个实施例中，第一噪音确定模块151进一步用于，按如下方法确定第一散射噪音：对第一地震数据和第二地震数据作差，将作差的结果作为第一散射噪音；
第二噪音确定模块152进一步用于，按如下方法确定第二散射噪音：对第一散射噪音和三维拉东反变换结果作差，将作差的结果作为第二散射噪音。
一个实施例中，第一噪音确定模块151进一步用于，对第一地震数据进行三维拉东正变换之前，对第一地震数据进行部分动校正；
噪音压制模块153进一步用于，对去除第一散射噪音和第二散射噪音之后的第一地震数据进行反动校正。
一个实施例中，所述三维拉东正变换为三维抛物拉东正变换，所述三维拉东反变换为三维抛物拉东反变换。
综上所述，本发明实施例通过获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音，然后对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音，最后从第一地震数据中去除第一散射噪音和第二散射噪音。本发明实施例通过两组三维拉东正变换与三维拉东反变换，完整提取出地震数据中存在的全部散射噪声，然后从地震数据中去除全部散射噪声，获得完整有效的地震信息，有效压制了地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题。实施例中，通过对第一地震数据按照炮线和检波线道头进行重新排序将数据由炮域转换到十字排列域，将线性干扰转变为圆锥体形态，便于去噪；在对第一地震数据进行三维拉东正变换之前，进行第一地震数据线性干扰压制处理和强能量干扰压制处理，有效衰减部分线性和强能量干扰；对三维地震数据进行处理，并采用三维拉东变换，提高地震数据中散射噪音的压制效果；对第一地震数据进行部分动校正，避免完全动校正将同相轴校正为水平状态的同时失去部分浅层地震数据；对地震数据进行部分动校正，将双曲线形式的同相轴转换为抛物线形式，然后采用计算精度和效率都非常高的三维抛物拉东正变换和三维抛物拉东反变换对地震数据进行处理，有效简化了运算过程，提高了计算精度和效率；在对地震数据进行一组三维拉东正变换与三维拉东反变换提取出一部分散射噪音之后，继续对这部分散射噪音进行三维拉东正变换与三维拉东反变换，提取地震数据中残留的散射噪音，通过两组三维拉东正变换与三维拉东反变换，完整提取出地震数据中存在的全部散射噪声，然后从地震数据中去除全部散射噪声，即可获得完整有效的地震信息，有效压制了地震数据中存在的散射噪音，在去除散射噪声的同时保证地震信息的完整性，避免因三维拉东变换采样不足引起的信息泄露问题。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (14)
1.一种地震数据散射噪音压制方法，其特征在于，包括：
获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音；
对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音；
从第一地震数据中去除第一散射噪音和第二散射噪音。
2.如权利要求1所述的方法，其特征在于，对第一地震数据进行三维拉东正变换之前，将第一地震数据变换至十字排列域。
3.如权利要求1所述的方法，其特征在于，对第一地震数据进行三维拉东正变换之前，进行第一地震数据线性干扰压制处理和强能量干扰压制处理。
4.如权利要求1所述的方法，其特征在于，按如下方法确定第一散射噪音：对第一地震数据和第二地震数据作差，将作差的结果作为第一散射噪音；
按如下方法确定第二散射噪音：对第一散射噪音和三维拉东反变换结果作差，将作差的结果作为第二散射噪音。
5.如权利要求1所述的方法，其特征在于，对第一地震数据进行三维拉东正变换之前，对第一地震数据进行部分动校正；
对去除第一散射噪音和第二散射噪音之后的第一地震数据进行反动校正。
6.如权利要求5所述的方法，其特征在于，所述三维拉东正变换为三维抛物拉东正变换，所述三维拉东反变换为三维抛物拉东反变换。
7.一种地震数据散射噪音压制装置，其特征在于，包括：
第一噪音确定模块，用于获取第一地震数据，对第一地震数据进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，得到第二地震数据，根据第一地震数据和第二地震数据，确定第一散射噪音；
第二噪音确定模块，用于对第一散射噪音进行三维拉东正变换，然后对三维拉东正变换结果进行三维拉东反变换，根据三维拉东反变换结果和第一散射噪音，确定第二散射噪音；
噪音压制模块，用于从第一地震数据中去除第一散射噪音和第二散射噪音。
8.如权利要求7所述的装置，其特征在于，第一噪音确定模块进一步用于，对第一地震数据进行三维拉东正变换之前，将第一地震数据变换至十字排列域。
9.如权利要求7所述的装置，其特征在于，第一噪音确定模块进一步用于，对第一地震数据进行三维拉东正变换之前，进行第一地震数据线性干扰压制处理和强能量干扰压制处理。
10.如权利要求7所述的装置，其特征在于，第一噪音确定模块进一步用于，按如下方法确定第一散射噪音：对第一地震数据和第二地震数据作差，将作差的结果作为第一散射噪音；
第二噪音确定模块进一步用于，按如下方法确定第二散射噪音：对第一散射噪音和三维拉东反变换结果作差，将作差的结果作为第二散射噪音。
11.如权利要求7所述的装置，其特征在于，第一噪音确定模块进一步用于，对第一地震数据进行三维拉东正变换之前，对第一地震数据进行部分动校正；
噪音压制模块进一步用于，对去除第一散射噪音和第二散射噪音之后的第一地震数据进行反动校正。
12.如权利要求11所述的装置，其特征在于，所述三维拉东正变换为三维抛物拉东正变换，所述三维拉东反变换为三维抛物拉东反变换。
13.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至6任一所述方法。
14.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至6任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910113832.7A CN111562616B (zh) | 2019-02-14 | 2019-02-14 | 地震数据散射噪音压制方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910113832.7A CN111562616B (zh) | 2019-02-14 | 2019-02-14 | 地震数据散射噪音压制方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111562616A true CN111562616A (zh) | 2020-08-21 |
| CN111562616B CN111562616B (zh) | 2023-06-30 |
Family
ID=72071304
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910113832.7A Active CN111562616B (zh) | 2019-02-14 | 2019-02-14 | 地震数据散射噪音压制方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111562616B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114740530A (zh) * | 2021-01-07 | 2022-07-12 | 中国石油天然气股份有限公司 | 基于双曲时窗约束的中高频拟线性噪声压制方法及装置 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004020972A2 (en) * | 2002-08-30 | 2004-03-11 | Robinson, John, M. | Removal of noise from seismic data using improved radon transformations |
| GB0426723D0 (en) * | 2001-01-23 | 2005-01-12 | Pgs Americas Inc | Weighted slant stack for attenuating seismic noise |
| US20050174885A1 (en) * | 2005-03-04 | 2005-08-11 | Robinson John M. | Radon transformations for removal of noise from seismic data |
| BRPI0701616A2 (pt) * | 2007-04-13 | 2008-11-25 | Petroleo Brasileiro Sa | mÉtodo geofÍsico para atenuaÇço de méltiplas |
| CN101334483A (zh) * | 2008-06-13 | 2008-12-31 | 徐基祥 | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 |
| CN102401908A (zh) * | 2010-09-07 | 2012-04-04 | 中国石油天然气集团公司 | 一种利用不同模加权稀疏的抛物拉东变换压制多次波的方法 |
| CN103869364A (zh) * | 2014-03-25 | 2014-06-18 | 中国石油大学(华东) | 一种基于双项抛物拉东变换的多次波压制方法 |
| CN106597539A (zh) * | 2016-12-28 | 2017-04-26 | 中国石油化工股份有限公司 | 针对黄土塬地区的曲波域Radon变换噪声压制方法 |
| CN107656313A (zh) * | 2017-10-12 | 2018-02-02 | 中国石油天然气集团公司 | 一种拉东域多次波压制方法及装置 |
| CN109001813A (zh) * | 2018-07-16 | 2018-12-14 | 中国石油天然气股份有限公司 | 一种压制多次波的方法、装置及系统 |
-
2019
- 2019-02-14 CN CN201910113832.7A patent/CN111562616B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0426723D0 (en) * | 2001-01-23 | 2005-01-12 | Pgs Americas Inc | Weighted slant stack for attenuating seismic noise |
| WO2004020972A2 (en) * | 2002-08-30 | 2004-03-11 | Robinson, John, M. | Removal of noise from seismic data using improved radon transformations |
| US20050174885A1 (en) * | 2005-03-04 | 2005-08-11 | Robinson John M. | Radon transformations for removal of noise from seismic data |
| BRPI0701616A2 (pt) * | 2007-04-13 | 2008-11-25 | Petroleo Brasileiro Sa | mÉtodo geofÍsico para atenuaÇço de méltiplas |
| CN101334483A (zh) * | 2008-06-13 | 2008-12-31 | 徐基祥 | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 |
| CN102401908A (zh) * | 2010-09-07 | 2012-04-04 | 中国石油天然气集团公司 | 一种利用不同模加权稀疏的抛物拉东变换压制多次波的方法 |
| CN103869364A (zh) * | 2014-03-25 | 2014-06-18 | 中国石油大学(华东) | 一种基于双项抛物拉东变换的多次波压制方法 |
| CN106597539A (zh) * | 2016-12-28 | 2017-04-26 | 中国石油化工股份有限公司 | 针对黄土塬地区的曲波域Radon变换噪声压制方法 |
| CN107656313A (zh) * | 2017-10-12 | 2018-02-02 | 中国石油天然气集团公司 | 一种拉东域多次波压制方法及装置 |
| CN109001813A (zh) * | 2018-07-16 | 2018-12-14 | 中国石油天然气股份有限公司 | 一种压制多次波的方法、装置及系统 |
Non-Patent Citations (3)
| Title |
|---|
| YARU XUE等: "Amplitude-preserving nonlinear adaptive multiple attenuation using the high-order sparse Radon transform", 《JOURNAL OF GEOPHYSICS AND ENGINEERING》 * |
| 巩向博等: "混合域高分辨率双曲Radon变换及其在多次波压制中的应用", 《石油物探》 * |
| 王立歆等: "局域双曲线Radon变换压制散射波研究与应用", 《CT理论与应用研究》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114740530A (zh) * | 2021-01-07 | 2022-07-12 | 中国石油天然气股份有限公司 | 基于双曲时窗约束的中高频拟线性噪声压制方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111562616B (zh) | 2023-06-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10677949B2 (en) | Anisotropy matching filtering for attenuation of seismic migration artifacts | |
| AU2013201072B2 (en) | Method and apparatus for automated noise removal from seismic data | |
| CN110907995B (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
| CN110261910A (zh) | 基于自适应稀疏s变换的地震数据面波去除方法 | |
| CN111736224B (zh) | 一种压制叠前地震资料线性干扰方法、存储介质及设备 | |
| CN111708087A (zh) | 一种基于DnCNN神经网络对地震数据噪声压制的方法 | |
| CN111562616B (zh) | 地震数据散射噪音压制方法及装置 | |
| CN110749925A (zh) | 一种宽频逆时偏移成像处理方法 | |
| Zhang et al. | A reverberation noise suppression method of sonar image based on shearlet transform | |
| CN110907989A (zh) | 重建拟地面地震反射波成像方法及系统 | |
| CN105676292A (zh) | 一种基于二维曲波变换的三维地震数据去噪方法 | |
| NO20130064A1 (no) | Innretning og fremgangsmate for a estimere tidsforskyvning. | |
| CN112198547A (zh) | 深层或超深层地震资料处理方法及装置 | |
| CN113866821B (zh) | 一种基于照明方向约束的被动源干涉偏移成像方法和系统 | |
| CN113156514B (zh) | 基于主频波数域均值滤波的地震数据去噪方法及系统 | |
| CN109212609A (zh) | 基于波动方程延拓的近地表噪音压制方法 | |
| WO2022250801A1 (en) | Seismic first breaks onset times determination by interface tracking using level-sets method | |
| CN111352158B (zh) | 地震信号增强方法及装置 | |
| CN111781644B (zh) | 浅中地层地震数据的线性干扰衰减方法及装置 | |
| CN110967733B (zh) | 一种改进的震源逆时定位方法及计算机可读存储介质 | |
| Jixiang et al. | A near-surface seismic scattered wave separation method | |
| CN113703039B (zh) | 逆时偏移成像方法及装置 | |
| CN113391348B (zh) | 用于叠前反演的共反射点道集构建方法及装置 | |
| CN113945979B (zh) | 一种基于数据重构的三维锥形滤波方法及装置 | |
| CN112415592B (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 |