CN111948714B - 海底结点地震数据中上行波的降噪方法及降噪装置 - Google Patents
海底结点地震数据中上行波的降噪方法及降噪装置 Download PDFInfo
- Publication number
- CN111948714B CN111948714B CN201910405250.6A CN201910405250A CN111948714B CN 111948714 B CN111948714 B CN 111948714B CN 201910405250 A CN201910405250 A CN 201910405250A CN 111948714 B CN111948714 B CN 111948714B
- Authority
- CN
- China
- Prior art keywords
- wave
- frequency
- dimensional fourier
- fourier transform
- 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 63
- 230000009467 reduction Effects 0.000 title claims abstract description 48
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 39
- 230000009466 transformation Effects 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 26
- 238000000926 separation method Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 8
- 230000021615 conjugation Effects 0.000 claims 2
- 238000003384 imaging method Methods 0.000 abstract description 18
- 230000008569 process Effects 0.000 abstract description 11
- 238000004891 communication Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 230000001174 ascending effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000001131 transforming effect Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/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
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
技术领域
本发明涉及地震勘探技术领域，具体涉及一种海底结点地震数据中上行波的降噪方法及降噪装置。
背景技术
海底结点(Ocean bottom node，OBN)地震数据采集是近年出现的一种海洋地震勘探技术。海底结点地震数据采集与常规拖缆采集相比，具有以下优点：采集的数据是多分量数据，不会受到水面钻井平台、船只的影响，基于海底结点地震数据的定位也更加准确；采集的数据具有更大的炮检距和更宽的方位角，以及海底结点地震数据信噪比更高。因此，海底结点地震数据采集越来越受到石油公司的重视，对海底结点地震数据的处理技术的要求也越来越高。海底结点地震数据的处理主要包括噪音压制、上下行波分离、速度分析和上下行波成像。而在上下行波分离的过程中，可以去除海底结点地震数据上行波中检波点的多次波(多次波是地震勘探数据中出现的一种干扰波的种类)，但是上行波中炮点端的多次波仍然存在。
为了压制上行波中炮点端的多次波，目前采用预测反褶积的方式，但是上行波中炮点端的多次波的周期并不是固定的，这就造成了采用预测反褶积的方式压制上行波中炮点端的多次波效果不理想。
随着海底结点地震数据采集技术的发展，对海底结点地震数据的精度要求越来越高，因此，迫切需要一种精确更高的处理方式来压制上行波中炮点端的多次波。
发明内容
针对现有技术中的问题，本发明提供一种海底结点地震数据中上行波的降噪方法及降噪装置，能够压制上行波中炮点端的多次波，提高海底结点地震数据的信噪比，进而提高地震勘探资料成像的精度。
为实现上述目的，本发明提供以下技术方案：
第一方面，本发明提供了一种海底结点地震数据中上行波的降噪方法，包括：
确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；
根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；
基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；
对所述地震数据进行二维傅里叶反变换得到降噪的上行波。
进一步地，所述确定海底结点地震数据的下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值之前，还包括：
对海底结点地震数据进行上下行波分离处理得到上行波和下行波。
进一步地，确定下行波经二维傅里叶变换后的频率和最大幅值，包括：
确定下行波经二维傅里叶变换后的频率和波数；
根据所述经二维傅里叶变换后的频率和所述波数计算振幅值并确定所述振幅值中的最大幅值。
进一步地，所述根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子，包括：
根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率和所述最大幅值计算白噪系数；
基于所述白噪系数和经二维傅里叶变换的下行波确定反褶积算子。
第二方面，本发明还提供了一种海底结点地震数据中上行波的降噪装置，包括：
变换单元，用于确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；
运算单元，用于根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；
反褶积单元，用于基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；
反变换单元，用于对所述地震数据进行二维傅里叶反变换得到降噪的上行波。
进一步地，还包括：
分离单元，用于对海底结点地震数据进行上下行波分离处理得到上行波和下行波。
进一步地，所述变换单元包括：
变换子单元，用于确定下行波经二维傅里叶变换后的频率和波数；
振幅子单元，用于根据所述经二维傅里叶变换后的频率和所述波数计算振幅值并确定所述振幅值中的最大幅值。
进一步地，所述运算单元包括：
白噪子单元，用于根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率和所述最大幅值计算白噪系数；
算子子单元，用于基于所述白噪系数和经二维傅里叶变换的下行波确定反褶积算子。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述海底结点地震数据中上行波的降噪方法的步骤。
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述海底结点地震数据中上行波的降噪方法的步骤。
由上述技术方案可知，本发明所述的一种海底结点地震数据中上行波的降噪方法及降噪装置，通过确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；对所述地震数据进行二维傅里叶反变换得到降噪的上行波，能够压制上行波中存在的炮点端的多次波，提高海底结点地震数据的信噪比，进而提高地震勘探资料成像的精度，以根据精度更高的地震勘探资料成像结果实现更为精确且可靠的地震勘探过程。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种海底结点地震数据中上行波的降噪方法的流程示意图；
图2是本发明实施例提供的一种海底结点地震数据中上行波的降噪方法中降噪前的上行波炮集记录图；
图3是本发明实施例提供的一种海底结点地震数据中上行波的降噪方法中降噪前的上行波共反射点叠加剖面图；
图4是本发明实施例提供的一种海底结点地震数据中上行波的降噪方法中降噪后，再延拓回海面的上行波炮集记录图；
图5是本发明实施例提供的一种海底结点地震数据中上行波的降噪方法中降噪后，再延拓回海面的上行波共反射点叠加剖面图；
图6是本发明实施例提供的另一种海底结点地震数据中上行波的降噪方法的流程示意图；
图7是本发明实施例提供的一种海底结点地震数据中上行波的降噪装置的结构示意图；
图8是本发明实施例提供的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
考虑到现有的海底结点地震数据的上下行波分离过程中，仅能去除海底结点地震数据上行波中检波点的多次波，无法去除上行波中炮点端的多次波，导致海底结点地震数据的信噪比较低，进而影响地震勘探资料成像的精度。因此，本发明提供一种海底结点地震数据中上行波的降噪方法及降噪装置，通过确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；对所述地震数据进行二维傅里叶反变换得到降噪的上行波，能够压制上行波中存在的炮点端的多次波，提高海底结点地震数据的信噪比，进而提高地震勘探资料成像的精度，以根据精度更高的地震勘探资料成像结果实现更为精确且可靠的地震勘探过程。
基于上述内容，本发明实施例提供一种海底结点地震数据中上行波的降噪方法，参见图1，该降噪方法具体内容包括：
S101：确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；
在本步骤中，对下行波d(x,ti)进行傅里叶变换，实现将下行波d(x,ti)由时间域转换至频率域。其中，下行波d(x,ti)中x是叠前数据对应的位置，ti是叠前数据对应的采样时间；i是采样的序列，其满足i＝1,2,…,n，n为单道数据总的采样数目。
在对下行波d(x,ti)进行一维傅里叶变换后确定变换后的下行波的主频fc和最小频率fmin。其中，主频fc是频谱极大值在时间域所对应的频率，最小频率fmin是各个频谱值在时间域所对应的频率中的最小的频率。
在对下行波d(x,ti)进行二维傅里叶变换后得到下行波D(Fj,Kl)，其中，Fj是频率轴上第j个频率，j是频率的序号；Kl是波数轴上第l个波数，l是波数的序号。根据频率Fj和波数Kl计算下行波D(Fj,Kl)的振幅值。在计算振幅值的过程中，根据N个频率中每个频率以及每个频率对应的N个波数计算得到N2个振幅值，并从中确定最大幅值A。
S102：根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；
在本步骤中，需要先根据主频、最小频率、经二维傅里叶变换后的频率和最大幅值计算白噪系数；然后基于白噪系数和经二维傅里叶变换的下行波确定反褶积算子。
S103：基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；
S104：对所述地震数据进行二维傅里叶反变换得到降噪的上行波。
从上述描述可知，本发明实施例提供的海底结点地震数据中上行波的降噪方法，能够压制上行波中存在的炮点端的多次波，提高海底结点地震数据的信噪比，进而提高地震勘探资料成像的精度，以根据精度更高的地震勘探资料成像结果实现更为精确且可靠的地震勘探过程。
在一种可选的实施方式中，本发明实施例提供了上述海底结点地震数据中上行波的降噪方法实施例中反褶积算子的具体计算方法，包括：
(1)采用下式计算白噪系数white：
其中，A为最大幅值，Fj是频率轴上第j个频率，fc为主频和fmin为最小频率。
(2)采用下式计算反褶积算子G(Fj,Kl)：
其中，D(Fj,Kl)为二维傅里叶变换后的下行波，是D(Fj,Kl)的共轭
相对应的，上述海底结点地震数据中上行波的降噪方法实施例中地震数据的具体计算方法。
采用下式计算地震数据P(Fj,Kl)：
P(Fj,Kl)＝U(Fj,Kl)G(Fj,Kl)；
其中，U(Fj,Kl)为二维傅里叶变换后的上行波。
从上述描述可知，本发明实施例提供地震数据的计算方法，能够压制上行波中存在的炮点端的多次波，提高海底结点地震数据的信噪比，进而提高地震勘探资料成像的精度，以根据精度更高的地震勘探资料成像结果实现更为精确且可靠的地震勘探过程。
基于上述实施例并结合具体应用例，对上述实施例中的方法进行说明。
参见图2所示的上行波的炮集记录图和图3所示的上行波中共反射点叠加剖面图，从图2和图3中可以看出，反褶积前的上行波中存在明显的炮点端多次波，在1400ms以下存在多组同相轴，除了1700ms这组同相轴外，其余都是多次波，在图2和图3之间，图3所示的上行波的共反射点叠加剖面图上炮点端的多次波更加明显。
对图2和图3所对应的上行波采用上述实施例中的方法进行降噪处理，得到图4和图5所对应的上行波，其中，图4是图2对应的上行波进行降噪处理后，再延拓回海面得到的上行波中炮集记录图，图5是图3对应的上行波进行降噪处理后，再延拓回海面得到的上行波的共反射点叠加剖面图。
从图4和图5中可以看出，上行波中的炮点端多次波能量明显削弱，海底结点地震数据的信噪比大大提高，为提高后续地震数据处理和偏移成像打下了良好基础。
在上述实施例的基础上，本发明实施例提供另一种海底结点地震数据中上行波的降噪方法，参见图6，在上述实施例中步骤S100之前，还包括：
S100：对海底结点地震数据进行上下行波分离处理得到上行波和下行波。
在本步骤中，利用海底结点采集海上地震勘探数据，取得叠前海底结点地震数据并对其进行上下行波分离，得到分离后的上行波u(x,ti)和下行波d(x,ti)；
相对应的，对分离后的上行波u(x,ti)进行二维傅里叶变换得到上行波U(Fj,Kl)。
在获取上行波u(x,ti)、下行波d(x,ti)和二维傅里叶变换得到上行波U(Fj,Kl)后，请参考上述海底结点地震数据中上行波的降噪方法实施例中的步骤S101至步骤S104，在此不在赘述。
从上述描述可知，本发明实施例提供的海底结点地震数据中上行波的降噪方法，首先根据上下行波分离后得到的下行波，确定一维傅里叶变换后的下行波的主频和最小频率以及二维傅里叶变换后的下行波的频率和振幅值的最大值，计算白噪系数以及反褶积算子，利用反褶积算子来完成上下行波反褶积，进而压制上行波中存在的炮点端多次波，从而提高海底结点地震数据的信噪比，提高海底结点地震数据的成像质量，以根据精度更高的地震勘探资料成像结果实现更为精确且可靠的地震勘探过程。
本发明实施例提供一种海底结点地震数据中上行波的降噪装置，参见图7，包括：
变换单元10，用于确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；
运算单元20，用于根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；
反褶积单元30，用于基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；
反变换单元40，用于对所述地震数据进行二维傅里叶反变换得到降噪的上行波。
进一步地，还包括：
分离单元50，用于对海底结点地震数据进行上下行波分离处理得到上行波和下行波。
进一步地，所述变换单元10包括：
变换子单元101，用于确定下行波经二维傅里叶变换后的频率和波数；
振幅子单元102，用于根据所述经二维傅里叶变换后的频率和所述波数计算振幅值并确定所述振幅值中的最大幅值。
进一步地，所述运算单元20包括：
白噪子单元201，用于根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率和所述最大幅值计算白噪系数；
算子子单元202，用于基于所述白噪系数和经二维傅里叶变换的下行波确定反褶积算子。
需要说明的是，本申请提供的海底结点地震数据中上行波的降噪装置的实施例具体可以用于执行上述实施例中的海底结点地震数据中上行波的降噪方法的实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
从上述描述可知，本发明实施例提供的一种海底结点地震数据中上行波的降噪装置，通过确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；对所述地震数据进行二维傅里叶反变换得到降噪的上行波，能够压制上行波中存在的炮点端的多次波，提高海底结点地震数据的信噪比，进而提高地震勘探资料成像的精度，以根据精度更高的地震勘探资料成像结果实现更为精确且可靠的地震勘探过程。
本发明实施例提供了一种电子设备，参见图8，该电子设备可以包括：处理器(processor)810、通信接口(Communications Interface)820、存储器(memory)830和通信总线840，其中，处理器810，通信接口820，存储器830通过通信总线840完成相互间的通信。处理器810可以调用存储器830中的逻辑指令，以执行如下方法：通过确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；对所述地震数据进行二维傅里叶反变换得到降噪的上行波。
此外，上述的存储器830中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本发明实施例提供了一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述各方法实施例所提供的方法，例如包括：通过确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；对所述地震数据进行二维傅里叶反变换得到降噪的上行波。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (8)
1.一种海底结点地震数据中上行波的降噪方法，其特征在于，包括：
确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；
根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；
基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；
对所述地震数据进行二维傅里叶反变换得到降噪的上行波；
具体地，根据主频、最小频率、经二维傅里叶变换后的频率和最大幅值计算白噪系数；基于白噪系数和经二维傅里叶变换的下行波确定反褶积算子；
其中，采用下式计算白噪系数white：
；
其中，A为最大幅值，是频率轴上第/>个频率，/>为主频，/>为最小频率；
其中，采用下式计算反褶积算子：：
；
其中，为二维傅里叶变换后的下行波，/>是的/>共轭；/>是波数轴上第/>个波数，/>是波数的序号。
2.根据权利要求1所述的海底结点地震数据中上行波的降噪方法，其特征在于，所述确定海底结点地震数据的下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值之前，还包括：
对海底结点地震数据进行上下行波分离处理得到上行波和下行波。
3.根据权利要求1或2所述的海底结点地震数据中上行波的降噪方法，其特征在于，确定下行波经二维傅里叶变换后的频率和最大幅值，包括：
确定下行波经二维傅里叶变换后的频率和波数；
根据所述经二维傅里叶变换后的频率和所述波数计算振幅值并确定所述振幅值中的最大幅值。
4.一种海底结点地震数据中上行波的降噪装置，其特征在于，包括：
变换单元，用于确定下行波经一维傅里叶变换后的主频和最小频率，以及经二维傅里叶变换后的频率和最大幅值；
运算单元，用于根据所述主频、所述最小频率、所述经二维傅里叶变换后的频率、所述最大幅值以及经二维傅里叶变换的下行波来获取反褶积算子；
反褶积单元，用于基于所述反褶积算子和经二维傅里叶变换的上行波进行反褶积处理得到地震数据；
反变换单元，用于对所述地震数据进行二维傅里叶反变换得到降噪的上行波；
具体地，根据主频、最小频率、经二维傅里叶变换后的频率和最大幅值计算白噪系数；基于白噪系数和经二维傅里叶变换的下行波确定反褶积算子；
其中，采用下式计算白噪系数white：
;
其中，A为最大幅值，是频率轴上第/>个频率，/>为主频，/>为最小频率；
其中，采用下式计算反褶积算子:
;
其中，为二维傅里叶变换后的下行波，/>是/>的共轭；/>是波数轴上第/>个波数，/>是波数的序号。
5.根据权利要求4所述的海底结点地震数据中上行波的降噪装置，其特征在于，还包括：
分离单元，用于对海底结点地震数据进行上下行波分离处理得到上行波和下行波。
6.根据权利要求4或5所述的海底结点地震数据中上行波的降噪装置，其特征在于，所述变换单元包括：
变换子单元，用于确定下行波经二维傅里叶变换后的频率和波数；
振幅子单元，用于根据所述经二维傅里叶变换后的频率和所述波数计算振幅值并确定所述振幅值中的最大幅值。
7.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至3任一项所述的海底结点地震数据中上行波的降噪方法。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至3任一项所述的海底结点地震数据中上行波的降噪方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910405250.6A CN111948714B (zh) | 2019-05-16 | 2019-05-16 | 海底结点地震数据中上行波的降噪方法及降噪装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910405250.6A CN111948714B (zh) | 2019-05-16 | 2019-05-16 | 海底结点地震数据中上行波的降噪方法及降噪装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111948714A CN111948714A (zh) | 2020-11-17 |
| CN111948714B true CN111948714B (zh) | 2023-10-27 |
Family
ID=73336433
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910405250.6A Active CN111948714B (zh) | 2019-05-16 | 2019-05-16 | 海底结点地震数据中上行波的降噪方法及降噪装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111948714B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009066057A2 (en) * | 2007-11-19 | 2009-05-28 | Geco Technology B.V. | Seismic data processing method for surface related multiple attenuation |
| CN107884829A (zh) * | 2017-12-06 | 2018-04-06 | 东华理工大学 | 一种联合压制浅海obc地震资料多次波的方法 |
| CN108363096A (zh) * | 2018-02-01 | 2018-08-03 | 中国石油天然气集团有限公司 | 一种水陆检波器地震数据上下行波场分离方法及装置 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020118602A1 (en) * | 2001-02-27 | 2002-08-29 | Sen Mrinal K. | Angle dependent surface multiple attenuation for two-component marine bottom sensor data |
| US7336561B2 (en) * | 2004-09-07 | 2008-02-26 | Pgs Americas, Inc. | System for attenuation of water bottom multiples in seismic data recorded by pressure sensors and particle motion sensors |
| US8971149B2 (en) * | 2012-05-10 | 2015-03-03 | Cggveritas Services Sa | Method, device and processing algorithm for multiple and noise elimination from marine seismic data |
| US9453930B2 (en) * | 2012-11-01 | 2016-09-27 | Pgs Geophysical As | Systems and methods for high-resolution imaging using separated wavefields |
| US9784870B2 (en) * | 2013-10-22 | 2017-10-10 | Cgg Services Sas | Demultiple using up/down separation of towed variable-depth streamer data |
-
2019
- 2019-05-16 CN CN201910405250.6A patent/CN111948714B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009066057A2 (en) * | 2007-11-19 | 2009-05-28 | Geco Technology B.V. | Seismic data processing method for surface related multiple attenuation |
| CN107884829A (zh) * | 2017-12-06 | 2018-04-06 | 东华理工大学 | 一种联合压制浅海obc地震资料多次波的方法 |
| CN108363096A (zh) * | 2018-02-01 | 2018-08-03 | 中国石油天然气集团有限公司 | 一种水陆检波器地震数据上下行波场分离方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| Attenuation of free-surface multiples by up/down deconvolution for marine towed-streamer data;Mariusz Majdański,et al.;《GEOPHYSICS》;20111231;第76卷(第6期);全文 * |
| Multidimensional signature deconvolution and free-surface multiple elimination of marine multicomponent ocean-bottom seismic data;Lasse Amundsen,et al.;《GEOPHYSICS》;20011031;第66卷(第5期);全文 * |
| What Comes Up Must Have Gone Down the Principle and Application of Up-Down Deconvolution for Multiple Attenuation of Ocean Bottom Data;Yi Wang,et al.;《CSEG RECORDER》;20150918;第16-20页及附图1-11 * |
| 上下行波场反褶积衰减OBC表面多次波方法研究;王兆旗等;《地球物理学进展》;20161231(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111948714A (zh) | 2020-11-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107193045B (zh) | 一种地震数据处理方法及装置 | |
| US11493595B2 (en) | Method and apparatus for evaluating accuracy in positioning a receiver point | |
| CN111781635B (zh) | 海底四分量弹性波高斯束深度偏移方法和装置 | |
| CN111948714B (zh) | 海底结点地震数据中上行波的降噪方法及降噪装置 | |
| CN111239832B (zh) | 3d3c-vsp成像处理方法、装置及设备 | |
| CN110703332B (zh) | 一种鬼波压制方法 | |
| CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
| Warner et al. | Bowhead whale localization using asynchronous hydrophones in the Chukchi Sea | |
| CN115292886A (zh) | 一种基于声场干涉现象的弹性海底声源运动参数估计系统及估计方法 | |
| CN108363103A (zh) | 规则干扰压制方法及装置 | |
| CN111352158B (zh) | 地震信号增强方法及装置 | |
| Thomas et al. | Deep learning for denoising high-rate global navigation satellite system data | |
| CN113031071B (zh) | 地震波长波长静校正方法及装置 | |
| CN115963565B (zh) | 基于双曲矢量中值滤波器的海底多次波衰减方法及装置 | |
| CN112241022B (zh) | 基于射线密度生成层析反演模型速度界面的方法和装置 | |
| CN110488348B (zh) | 深度域地震数据的拓频处理方法、装置及存储介质 | |
| WO2024016572A1 (zh) | 噪声识别方法、装置、设备及存储介质 | |
| CN113009569B (zh) | 一种地震偏移成像方法及装置 | |
| CN112180448A (zh) | 剩余静校正方法、装置和电子设备 | |
| CN117309126B (zh) | 一种消除浅海环境水面反射的近场辐射噪声测量方法 | |
| CN112230285B (zh) | 地震数据激发延迟时间的校正方法及校正装置 | |
| CN115598696A (zh) | 水体速度模型质控方法、装置、电子设备及可读存储介质 | |
| CN116027413A (zh) | 基于反射波波形反演梯度预处理的背景速度场重构方法 | |
| CN117520736A (zh) | 海洋obn资料假频直达波干扰压制方法、系统、存储介质 | |
| CN109242781B (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 |