CN111781642A - 地震数据层间多次波衰减方法及装置 - Google Patents
地震数据层间多次波衰减方法及装置 Download PDFInfo
- Publication number
- CN111781642A CN111781642A CN201910264852.4A CN201910264852A CN111781642A CN 111781642 A CN111781642 A CN 111781642A CN 201910264852 A CN201910264852 A CN 201910264852A CN 111781642 A CN111781642 A CN 111781642A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- dip
- attenuation
- filtering
- previous reference
- 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
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
- G01V1/364—Seismic filtering
-
- 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
- G01V2210/324—Filtering
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为本发明实施例中地震数据层间多次波衰减装置结构图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
为了对叠前偏移之后的地震数据进行层间多次波衰减，有效去除地震数据中各偏移距范围内的层间多次波，克服近偏移距范围内的层间多次波因与有效信号之间的时差小而无法被完全衰减的问题，从而获得更好的衰减效果，本发明实施例提供一种地震数据层间多次波衰减方法，如图1所示，该方法可以包括：
步骤101、获取叠前偏移处理后的地震数据，对所述地震数据进行倾角滤波；
步骤102、从倾角滤波后的地震数据中，选取倾角角度最小的设定空间长度的地震数据，作为前次参考地震数据；
步骤103、沿地震数据的空间方向，按如下方式对倾角滤波后的地震数据进行衰减处理：
选取与前次参考地震数据相邻的设定空间长度的地震数据，其中，选取的与前次参考地震数据相邻的地震数据中包含设定数量的前次参考地震数据；
对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减；
步骤104、将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理。
由图1所示可以得知，本发明实施例通过获取叠前偏移处理后的地震数据，对所述地震数据进行倾角滤波，从倾角滤波后的地震数据中，选取倾角角度最小的设定空间长度的地震数据，作为前次参考地震数据，沿地震数据的空间方向，按如下方式对倾角滤波后的地震数据进行衰减处理：选取与前次参考地震数据相邻的设定空间长度的地震数据，其中，选取的与前次参考地震数据相邻的地震数据中包含设定数量的前次参考地震数据，对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减，将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理。本发明实施例首先利用倾角滤波对叠前偏移处理后的地震数据进行初步的层间多次波衰减，以去除处于中远偏移距范围的层间多次波，然后通过选取倾角角度最小的地震数据，确定层间多次波衰减最完全的空间位置，以该位置作为起始位置，沿地震数据的空间方向向残留多次波较多的地震道滑动，选取相邻的地震数据，在选取的地震数据中，多次波相当于随机干扰，通过进行随机噪声衰减的方式可以将地震数据中残留在近偏移距的层间多次波去除，从而克服了近偏移距范围内的层间多次波因与有效信号之间的时差小而无法被完全衰减的问题，获得了更好的衰减效果。
具体实施时，获取叠前偏移处理后的地震数据，对所述地震数据进行倾角滤波。
实施例中，首先获取叠前偏移处理后的地震数据，对叠前偏移处理后的地震数据按共反射点和偏移距进行排序，然后对排序后的地震数据进行倾角滤波，以去除处于中远偏移距范围的层间多次波。
具体实施时，从倾角滤波后的地震数据中，选取倾角角度最小的设定空间长度的地震数据，作为前次参考地震数据。
利用现有技术对叠前偏移之后的地震数据进行衰减处理后，大部分层间多次波可以被有效的衰减，尤其是处于中远偏移距范围的层间多次波，可以被完全衰减，这是因为随着偏移距的增加，层间多次波与有效信号的时差也在不断增加，而层间多次波与有效信号的时差越大，越有利于获得更好的衰减效果，因此在中远偏移距范围内的层间多次波可以被完全衰减。而对于近偏移距范围内的地震数据，由于层间多次波与有效信号之间的时差小，层间多次波的衰减效果不好，无法被完全衰减。发明人发现，多次波虽然是相关信号，具备较强的周期性，但是仅仅对当前存在多次波的地震道成立，相对于不含多次波的信号来说，多次信号为随机信号。当中远偏移距层间多次波被衰减后，近偏移距多次波相对于中远偏移距地震信号来说为随机干扰，可以通过道集排序加随机干扰衰减的方式加以去除。因此，本发明实施例首先利用倾角滤波对叠前偏移处理后的地震数据进行初步的层间多次波衰减，以去除处于中远偏移距范围的层间多次波，然后通过选取倾角角度最小的地震数据，确定层间多次波衰减最完全的空间位置，以该位置作为起始位置，沿地震数据的空间方向向残留多次波较多的地震道滑动，选取相邻的地震数据，在选取的地震数据中，多次波相当于随机干扰，通过进行随机噪声衰减的方式可以将地震数据中残留在近偏移距的层间多次波去除，从而克服了近偏移距范围内的层间多次波因与有效信号之间的时差小而无法被完全衰减的问题，获得了更好的衰减效果。
实施例中，倾角角度越小，代表地震数据中包含的层间多次波越少，通过选取倾角角度最小的设定空间长度的地震数据，可以找到地震数据中进行倾角滤波后层间多次波衰减最完全的数据段，作为前次参考地震数据。
实施例中，按如下方法从倾角滤波后的地震数据中选取倾角角度最小的设定空间长度的地震数据：沿地震数据的空间方向，顺序选取设定空间长度的地震数据，确定每一设定空间长度的地震数据的倾角角度，选取最小倾角角度所对应的设定空间长度的地震数据。
具体实施时，沿地震数据的空间方向，按如下方式对倾角滤波后的地震数据进行衰减处理：
选取与前次参考地震数据相邻的设定空间长度的地震数据，其中，选取的与前次参考地震数据相邻的地震数据中包含设定数量的前次参考地震数据；
对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减。
实施例中，首先选取与前次参考地震数据相邻的设定空间长度的地震数据，其中，选取的与前次参考地震数据相邻的地震数据中包含设定数量的前次参考地震数据，然后对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减。举一例，设前次参考地震数据的空间长度为w，中间地震道为D(i)，其中i为前次参考地震数据中间地震道的偏移距，前次参考地震数据的地震道定义为由D(i-w/2)开始，D(i+w/2)截止，令最小偏移距为Offsetmin，最大偏移距为Offsetmax，i分别向着Offsetmin和Offsetmax变化，随之选取到的与前次参考地震数据相邻的地震数据整体随之发生变化。
实施例中，选取与前次参考地震数据相邻的设定空间长度的地震数据之后，对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减。随机噪声衰减为现有技术，本领域技术人员能够通过查阅资料了随机噪声衰减方法，本发明不再进行具体说明。
具体实施时，将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理。
实施例中，在对前次参考地震数据进行随机噪声衰减后，其中的层间多次波已被去除，然后将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理。
实施例中，完成全部倾角滤波后的地震数据的衰减处理之后，将倾角滤波后的地震数据和衰减处理后的结果作差，对作差结果进行倾角滤波处理，将衰减处理后的结果和对作差结果进行倾角滤波处理之后的结果求和。发明人发现，进行多次波衰减后，难免会将有效信号也当成层间多次波干扰衰减了，这部分被随之衰减掉的有效信号为残留有效信号，这部分残留有效信号呈现水平方向的特征，因此，本发明实施例对完成全部倾角滤波后的地震数据的衰减处理之后的数据进行倾角滤波，将这部分有效信号提取出来，再加到多次波衰减以后的结果上。
下面给出一个具体实施例，说明本发明实施例中地震数据层间多次波衰减方法的具体应用。在本具体实施例中，首先获取叠前偏移处理后的地震数据，对叠前偏移处理后的地震数据按共反射点和偏移距进行排序，如图2所示，然后对排序后的地震数据进行倾角滤波，以去除处于中远偏移距范围的层间多次波，结果如图3所示。从倾角滤波后的地震数据中，选取倾角角度最小的设定空间长度的地震数据，作为前次参考地震数据，如图3中黑色粗线框出的部分。选取与前次参考地震数据相邻的设定空间长度的地震数据，对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减。将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理，结果如图4所示。完成全部倾角滤波后的地震数据的衰减处理之后，将倾角滤波后的地震数据和衰减处理后的结果作差，对作差结果进行倾角滤波处理，将衰减处理后的结果和对作差结果进行倾角滤波处理之后的结果求和，结果如图5所示。图6为未进行层间多次波衰减的叠加剖面，从剖面中可以看到同相轴非常多，已经远远超过实际地层反射应该有的数量，图7为按本发明实施例的层间多次波衰减方法进行处理后的叠加剖面，层间多次波衰减后波组特征更加明显，地层假象得到消除。图8为叠前时间偏移后的共反射点道集，其中存在大量层间多次，尤其近偏移距受到多次干扰的影响有效信号被完全压制，图9为按本发明实施例的层间多次波衰减方法进行处理后的共反射点道集，近偏移距范围内的有效信号得到恢复。
基于同一发明构思，本发明实施例还提供了一种炼化装置尾气污染物在线监测数据预警装置，如下面的实施例所述。由于这些解决问题的原理与炼化装置尾气污染物在线监测数据预警方法相似，因此装置的实施可以参见方法的实施，重复之处不再赘述。
图10为本发明实施例中地震数据层间多次波衰减装置的结构图，如图10所示，该装置包括：
地震数据获取模块110，用于获取叠前偏移处理后的地震数据，对所述地震数据进行倾角滤波；
参考数据获取模块210，用于从倾角滤波后的地震数据中，选取倾角角度最小的设定空间长度的地震数据，作为前次参考地震数据；
衰减处理模块310，用于沿地震数据的空间方向，按如下方式对倾角滤波后的地震数据进行衰减处理：
选取与前次参考地震数据相邻的设定空间长度的地震数据，其中，选取的与前次参考地震数据相邻的地震数据中包含设定数量的前次参考地震数据；
对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减；
衰减迭代模块410，用于将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理。
一个实施例中，数据获取模块110进一步用于：对所述地震数据进行倾角滤波之前，对所述叠前偏移处理后的地震数据按共反射点和偏移距进行排序；
对所述地震数据进行倾角滤波，包括：对排序后的地震数据进行倾角滤波。
一个实施例中，参考数据获取模块210进一步用于：按如下方法从倾角滤波后的地震数据中选取倾角角度最小的设定空间长度的地震数据：沿地震数据的空间方向，顺序选取设定空间长度的地震数据，确定每一设定空间长度的地震数据的倾角角度，选取最小倾角角度所对应的设定空间长度的地震数据。
一个实施例中，地震数据层间多次波衰减装置还包括：
结果校正模块510，用于完成全部倾角滤波后的地震数据的衰减处理之后，将倾角滤波后的地震数据和衰减处理后的结果作差，对作差结果进行倾角滤波处理，将衰减处理后的结果和对作差结果进行倾角滤波处理之后的结果求和。
综上所述，本发明实施例通过获取叠前偏移处理后的地震数据，对所述地震数据进行倾角滤波，从倾角滤波后的地震数据中，选取倾角角度最小的设定空间长度的地震数据，作为前次参考地震数据，沿地震数据的空间方向，按如下方式对倾角滤波后的地震数据进行衰减处理：选取与前次参考地震数据相邻的设定空间长度的地震数据，其中，选取的与前次参考地震数据相邻的地震数据中包含设定数量的前次参考地震数据，对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减，将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理。本发明实施例首先利用倾角滤波对叠前偏移处理后的地震数据进行初步的层间多次波衰减，以去除处于中远偏移距范围的层间多次波，然后通过选取倾角角度最小的地震数据，确定层间多次波衰减最完全的空间位置，以该位置作为起始位置，沿地震数据的空间方向向残留多次波较多的地震道滑动，选取相邻的地震数据，在选取的地震数据中，多次波相当于随机干扰，通过进行随机噪声衰减的方式可以将地震数据中残留在近偏移距的层间多次波去除，从而克服了近偏移距范围内的层间多次波因与有效信号之间的时差小而无法被完全衰减的问题，获得了更好的衰减效果。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种地震数据层间多次波衰减方法，其特征在于，包括：
获取叠前偏移处理后的地震数据，对所述地震数据进行倾角滤波；
从倾角滤波后的地震数据中，选取倾角角度最小的设定空间长度的地震数据，作为前次参考地震数据；
沿地震数据的空间方向，按如下方式对倾角滤波后的地震数据进行衰减处理：
选取与前次参考地震数据相邻的设定空间长度的地震数据，其中，选取的与前次参考地震数据相邻的地震数据中包含设定数量的前次参考地震数据；
对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减；
将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理。
2.如权利要求1所述的方法，其特征在于，对所述地震数据进行倾角滤波之前，对所述叠前偏移处理后的地震数据按共反射点和偏移距进行排序；
对所述地震数据进行倾角滤波，包括：对排序后的地震数据进行倾角滤波。
3.如权利要求1所述的方法，其特征在于，按如下方法从倾角滤波后的地震数据中选取倾角角度最小的设定空间长度的地震数据：沿地震数据的空间方向，顺序选取设定空间长度的地震数据，确定每一设定空间长度的地震数据的倾角角度，选取最小倾角角度所对应的设定空间长度的地震数据。
4.如权利要求1所述的方法，其特征在于，还包括：
完成全部倾角滤波后的地震数据的衰减处理之后，将倾角滤波后的地震数据和衰减处理后的结果作差，对作差结果进行倾角滤波处理，将衰减处理后的结果和对作差结果进行倾角滤波处理之后的结果求和。
5.一种地震数据层间多次波衰减装置，其特征在于，包括：
地震数据获取模块，用于获取叠前偏移处理后的地震数据，对所述地震数据进行倾角滤波；
参考数据获取模块，用于从倾角滤波后的地震数据中，选取倾角角度最小的设定空间长度的地震数据，作为前次参考地震数据；
衰减处理模块，用于沿地震数据的空间方向，按如下方式对倾角滤波后的地震数据进行衰减处理：
选取与前次参考地震数据相邻的设定空间长度的地震数据，其中，选取的与前次参考地震数据相邻的地震数据中包含设定数量的前次参考地震数据；
对选取的与前次参考地震数据相邻的地震数据进行随机噪声衰减；
衰减迭代模块，用于将进行随机噪声衰减后的地震数据作为新的前次参考地震数据，重新沿地震数据的空间方向，按上述方式对倾角滤波后的地震数据进行衰减处理，直至完成全部倾角滤波后的地震数据的衰减处理。
6.如权利要求5所述的装置，其特征在于，所述数据获取模块进一步用于：对所述地震数据进行倾角滤波之前，对所述叠前偏移处理后的地震数据按共反射点和偏移距进行排序；
对所述地震数据进行倾角滤波，包括：对排序后的地震数据进行倾角滤波。
7.如权利要求5所述的装置，其特征在于，所述参考数据获取模块进一步用于：按如下方法从倾角滤波后的地震数据中选取倾角角度最小的设定空间长度的地震数据：沿地震数据的空间方向，顺序选取设定空间长度的地震数据，确定每一设定空间长度的地震数据的倾角角度，选取最小倾角角度所对应的设定空间长度的地震数据。
8.如权利要求5所述的装置，其特征在于，还包括：
结果校正模块，用于完成全部倾角滤波后的地震数据的衰减处理之后，将倾角滤波后的地震数据和衰减处理后的结果作差，对作差结果进行倾角滤波处理，将衰减处理后的结果和对作差结果进行倾角滤波处理之后的结果求和。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910264852.4A CN111781642B (zh) | 2019-04-03 | 2019-04-03 | 地震数据层间多次波衰减方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910264852.4A CN111781642B (zh) | 2019-04-03 | 2019-04-03 | 地震数据层间多次波衰减方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111781642A true CN111781642A (zh) | 2020-10-16 |
| CN111781642B CN111781642B (zh) | 2022-07-05 |
Family
ID=72754743
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910264852.4A Active CN111781642B (zh) | 2019-04-03 | 2019-04-03 | 地震数据层间多次波衰减方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111781642B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113514888A (zh) * | 2021-06-18 | 2021-10-19 | 中国石油化工股份有限公司 | 一种基于非剩余时差的多次波压制方法 |
| CN115993654A (zh) * | 2021-10-20 | 2023-04-21 | 中国石油化工股份有限公司 | 有效压制道集近偏移距剩余多次波的方法、装置及设备 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1997011390A2 (en) * | 1995-09-19 | 1997-03-27 | Exxon Production Research Company | Multiple suppression in geophysical data |
| CN102854533A (zh) * | 2011-07-01 | 2013-01-02 | 中国石油化工股份有限公司 | 一种基于波场分离原理提高地震资料信噪比的去噪方法 |
| CN103605157A (zh) * | 2013-10-14 | 2014-02-26 | 中国石油天然气股份有限公司 | 衰减近地表散射波的方法 |
| CN104081226A (zh) * | 2011-12-15 | 2014-10-01 | 沙特阿拉伯石油公司 | 用于地震数据处理的迭代倾角导向中值滤波器 |
| US20170068009A1 (en) * | 2015-09-04 | 2017-03-09 | Westerngeco L.L.C. | Method for attenuating multiples |
| CN108196304A (zh) * | 2018-01-09 | 2018-06-22 | 中国石油天然气股份有限公司 | 一种多次波压制方法及装置 |
-
2019
- 2019-04-03 CN CN201910264852.4A patent/CN111781642B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1997011390A2 (en) * | 1995-09-19 | 1997-03-27 | Exxon Production Research Company | Multiple suppression in geophysical data |
| CN102854533A (zh) * | 2011-07-01 | 2013-01-02 | 中国石油化工股份有限公司 | 一种基于波场分离原理提高地震资料信噪比的去噪方法 |
| CN104081226A (zh) * | 2011-12-15 | 2014-10-01 | 沙特阿拉伯石油公司 | 用于地震数据处理的迭代倾角导向中值滤波器 |
| CN103605157A (zh) * | 2013-10-14 | 2014-02-26 | 中国石油天然气股份有限公司 | 衰减近地表散射波的方法 |
| US20170068009A1 (en) * | 2015-09-04 | 2017-03-09 | Westerngeco L.L.C. | Method for attenuating multiples |
| CN108196304A (zh) * | 2018-01-09 | 2018-06-22 | 中国石油天然气股份有限公司 | 一种多次波压制方法及装置 |
Non-Patent Citations (16)
| Title |
|---|
| SHOU HAO ET AL.: "Application of CRS based migration in low SNR region", 《SPG/SEG SHENZHEN 2011 INTERNATIONAL GEOPHYSICAL CONFERENCE TECHNICAL PROGRAM EXPANDED ABSTRACTS》 * |
| SHOU HAO ET AL.: "Application of CRS based migration in low SNR region", 《SPG/SEG SHENZHEN 2011 INTERNATIONAL GEOPHYSICAL CONFERENCE TECHNICAL PROGRAM EXPANDED ABSTRACTS》, 31 December 2011 (2011-12-31), pages 1404 - 1407 * |
| WILLIAM H ET AL.: "Some remarks on surface multiple attenuation", 《GEOPHYSICS》 * |
| WILLIAM H ET AL.: "Some remarks on surface multiple attenuation", 《GEOPHYSICS》, vol. 63, no. 2, 30 April 1998 (1998-04-30), pages 772 - 789, XP002314826, DOI: 10.1190/1.1444377 * |
| 冉建斌等: "冀东滩海南堡2号构造复杂油气储层叠前地震描述技术及效果", 《石油地球物理勘探》 * |
| 冉建斌等: "冀东滩海南堡2号构造复杂油气储层叠前地震描述技术及效果", 《石油地球物理勘探》, vol. 43, no. 01, 29 February 2008 (2008-02-29), pages 59 - 68 * |
| 戴晓峰等: "Ｒ ａｄｏｎ变换压制层间多次波技术在高石梯—磨溪地区的应用", 《石油学报》 * |
| 戴晓峰等: "Ｒ ａｄｏｎ变换压制层间多次波技术在高石梯—磨溪地区的应用", 《石油学报》, vol. 39, no. 9, 30 September 2018 (2018-09-30), pages 1028 - 1036 * |
| 李丽青等: "海上单道地震资料中多次波的衰减", 《石油物探》 * |
| 李丽青等: "海上单道地震资料中多次波的衰减", 《石油物探》, vol. 46, no. 05, 30 September 2007 (2007-09-30), pages 457 - 462 * |
| 王典等: "反射法地震勘探噪声消除技术研究", 《地球物理学进展》 * |
| 王典等: "反射法地震勘探噪声消除技术研究", 《地球物理学进展》, vol. 21, no. 03, 30 September 2006 (2006-09-30), pages 957 - 970 * |
| 苑益军等: "去噪技术在地震资料处理中的应用", 《东华理工学院学报(自然科学版)》 * |
| 苑益军等: "去噪技术在地震资料处理中的应用", 《东华理工学院学报(自然科学版)》, vol. 28, no. 01, 31 March 2005 (2005-03-31), pages 12 - 16 * |
| 马继涛等: "基于一次波逆时理论的层间多次波衰减方法", 《石油地球物理勘探》 * |
| 马继涛等: "基于一次波逆时理论的层间多次波衰减方法", 《石油地球物理勘探》, vol. 48, no. 2, 30 April 2013 (2013-04-30), pages 181 - 186 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113514888A (zh) * | 2021-06-18 | 2021-10-19 | 中国石油化工股份有限公司 | 一种基于非剩余时差的多次波压制方法 |
| CN113514888B (zh) * | 2021-06-18 | 2024-04-16 | 中国石油化工股份有限公司 | 一种基于非剩余时差的多次波压制方法 |
| CN115993654A (zh) * | 2021-10-20 | 2023-04-21 | 中国石油化工股份有限公司 | 有效压制道集近偏移距剩余多次波的方法、装置及设备 |
| CN115993654B (zh) * | 2021-10-20 | 2025-10-17 | 中国石油化工股份有限公司 | 有效压制道集近偏移距剩余多次波的方法、装置及设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111781642B (zh) | 2022-07-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110658557B (zh) | 基于生成对抗网络的地震数据面波压制方法 | |
| CN112013286B (zh) | 管道泄漏点的定位方法及装置、存储介质、终端 | |
| CN105510975B (zh) | 提高地震数据信噪比的方法及装置 | |
| CN108267784A (zh) | 一种地震信号随机噪声压制处理方法 | |
| CN108828670B (zh) | 一种地震数据降噪方法 | |
| CN113885077A (zh) | 一种基于深度学习的多震源地震数据分离方法 | |
| CN105445801B (zh) | 一种消除二维地震资料随机噪音的处理方法 | |
| CN111781642A (zh) | 地震数据层间多次波衰减方法及装置 | |
| CN111337981A (zh) | 绕射波成像方法、装置和电子设备 | |
| CN112926232B (zh) | 一种基于分层融合的地震低频分量恢复方法 | |
| CN112255689B (zh) | 多次波发育区地震数据保真速度分析方法 | |
| CN106950600A (zh) | 一种近地表散射面波的去除方法 | |
| CN109212609B (zh) | 基于波动方程延拓的近地表噪音压制方法 | |
| CN110008633B (zh) | 基于人工智能深度神经网络的公路噪音压制方法及系统 | |
| CN112782766B (zh) | 去除地震资料侧源干扰的方法及装置 | |
| CN111562616B (zh) | 地震数据散射噪音压制方法及装置 | |
| CN112764109B (zh) | 偶极横波反射波的分离提取方法及装置 | |
| CN113946778A (zh) | 欠采样混合震源数据处理方法及装置 | |
| Chen et al. | Learning the blending spikes using sparse dictionaries | |
| CN117388920A (zh) | 一种基于迭代地震干涉法的近偏移距数据重建方法 | |
| CN112147694B (zh) | 一种面向主成分的avo趋势去除多次波的方法 | |
| CN108957522A (zh) | 一种基于克希霍夫积分的海底地震仪数据多次波预测方法 | |
| US9330441B2 (en) | Automated selection of filter parameters for seismic analysis | |
| CN113093282A (zh) | 一种基于几何模态特征并行网络的沙漠数据消噪方法 | |
| CN119738871B (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 |