CN107656313B - 一种拉东域多次波压制方法及装置 - Google Patents
一种拉东域多次波压制方法及装置 Download PDFInfo
- Publication number
- CN107656313B CN107656313B CN201710945351.3A CN201710945351A CN107656313B CN 107656313 B CN107656313 B CN 107656313B CN 201710945351 A CN201710945351 A CN 201710945351A CN 107656313 B CN107656313 B CN 107656313B
- Authority
- CN
- China
- Prior art keywords
- trace gather
- wave
- multiple wave
- domain
- original
- 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 61
- 238000012937 correction Methods 0.000 claims abstract description 64
- 238000004458 analytical method Methods 0.000 claims abstract description 20
- 238000006243 chemical reaction Methods 0.000 claims abstract description 20
- 230000001629 suppression Effects 0.000 claims abstract description 18
- 229910052704 radon Inorganic materials 0.000 claims abstract description 16
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 16
- 230000009466 transformation Effects 0.000 claims abstract description 12
- 238000004590 computer program Methods 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 16
- 238000012545 processing Methods 0.000 description 11
- 238000003860 storage Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000000746 purification Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000007630 basic procedure Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001131 transforming 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)
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供了一种拉东域多次波压制方法及装置，该方法包括：拾取去噪后的原始道集的第一一次波速度，并用其对原始道集作动校正；将动校正后道集作拉东变换，获得拉东域的一次波及多次波；将多次波能量置零，并将一次波及能量置零后的多次波反变换至TX域，获得新道集；将原始道集与新道集相减获得差道集；对新道集作速度分析获得第二一次波速度；用第一一次波速度对差道集作反动校正，并用第二一次波速度对反动校正后的差道集作动校正；对动校正后的差道集作FK变换以获取多次波能量区域，并将其反变换至TX域获得多次波道集；将原始道集与多次波道集相减，获得压制多次波后的道集。本申请可在准确压制多次波的同时，避免反射信息的泄漏。
Description
技术领域
本申请涉及地震数据处理技术领域，尤其是涉及一种拉东域多次波压制方法及装置。
背景技术
多次波广泛存在于地震资料当中。在一般的反射波地震勘探中，多次波属于干扰波，多次波的存在掩盖有效反射信息，或者产生虚假地层，因此多次波的压制一直是地震资料处理中的重要课题。
目前多次波压制方法主要有两大类：以自由界面多次波预测与压制(SRME)为代表的预测模型相减法以及数学变换滤波法(包括Radon，FK，KL等数学变换)。其中，拉东域(Radon域)压制多次波方法是目前最为常用的滤波法之一，已经属于地震资料处理中的常规和基本流程。该方法主要通过一次波和多次波在相同位置的速度差异进行多次波压制，多次波和有效波区域界限能否分明，依赖于获得的一次波速度是否准确。
然而，在没有进行多次波压制之前，由于多次波掩盖的因素，很难得到相对准确的一次波速度。总而言之，速度的准确与否直接影响了拉东域多次波压制的效果，而同时多次波的存在又直接影响了速度的准确性，因此要想较好的压制多次波，就很难避免有效反射信息的泄漏。由此可见，常规的拉东域多次波压制方案难以在准确压制多次波的同时，避免有效反射信息的泄漏。
发明内容
本申请实施例的目的在于提供一种拉东域多次波压制方法及装置，以实现在准确压制多次波的同时，避免有效反射信息的泄漏。
为达到上述目的，一方面，本申请实施例提供了一种拉东域多次波压制方法，包括：
对去噪后的原始道集进行速度分析，拾取第一一次波速度；
用所述第一一次波速度对所述原始道集进行动校正；
将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域；
将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集；
将所述原始道集与所述新道集进行相减，获得差道集；
对所述新道集进行速度分析，获得第二一次波速度；
用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正；
对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集；
将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
优选的，所述原始道集包括CMP道集。
优选的，所述将所述多次波区域的能量置零，包括：
将所述多次波区域中其能量值大于零的部分置零。
优选的，所述将所述原始道集与所述多次波道集相减，包括：
对所述原始道集与所述多次波道集进行相位和振幅的匹配；
将匹配后获得的原始道集与多次波道集进行相减。
优选的，所述对所述原始道集与所述多次波道集进行相位和振幅的匹配，包括：
利用最小平方滤波法对所述原始道集与所述多次波道集进行相位和振幅的匹配。
另一方面，本申请实施例还提供了一种拉东域多次波压制装置，包括：
第一波速拾取模块，用于对去噪后的原始道集进行速度分析，拾取第一一次波速度；
第一校正模块，用于用所述第一一次波速度对所述原始道集进行动校正；
拉东变换模块，用于将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域；
新道集获取模块，用于将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集；
差道集获取模块，用于将所述原始道集与所述新道集进行相减，获得差道集；
第二波速拾取模块，用于对所述新道集进行速度分析，获得第二一次波速度；
第二校正模块，用于用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正；
多次波道集获取模块，用于对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集；
多次波道集压制模块，用于将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
优选的，所述原始道集包括CMP道集。
优选的，所述将所述多次波区域的能量置零，包括：
将所述多次波区域中其能量值大于零的部分置零。
优选的，所述将所述原始道集与所述多次波道集相减，包括：
对所述原始道集与所述多次波道集进行相位和振幅的匹配；
将匹配后获得的原始道集与多次波道集进行相减。
优选的，所述对所述原始道集与所述多次波道集进行相位和振幅的匹配，包括：
利用最小平方滤波法对所述原始道集与所述多次波道集进行相位和振幅的匹配。
再一方面，本申请实施例还提供了另一种拉东域多次波压制装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：
对去噪后的原始道集进行速度分析，拾取第一一次波速度；
用所述第一一次波速度对所述原始道集进行动校正；
将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域；
将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集；
将所述原始道集与所述新道集进行相减，获得差道集；
对所述新道集进行速度分析，获得第二一次波速度；
用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正；
对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集；
将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
由以上本申请实施例提供的技术方案可见，本申请实施例在拉东域可以最大限度地切除多次波能量区域，以彻底压制多次波；而不用担心会损失掉一部分一次波，因为在后续FK域提纯过程中，通过将原始道集与新多次波道集进行相减，可以找回之前损失的一次有效波，从而即达到了在精准压制多次波的同时，又不损失有效信号的目的，因此，本申请实施例的拉东域多次波压制方法相对常规的拉东域去多次波更有优势。此外，由于后续会找回之前损失的一次有效波，使得本申请实施例可以容忍更大的一次波速度误差。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本申请一实施例的拉东域多次波压制方法的流程图；
图2为本申请一实施例中拉东域多次波求取示意图；
图3为本申请一实施例中FK域多次波提纯示意图；
图4为本申请一实施例中原始地震资料道集示意图；
图5为采用现有多次波压制技术处理图4所示原始地震资料道集而获得的道集示意图；
图6为采用本申请实施例的拉东域多次波压制方法处理图4所示原始地震资料道而获得的道集示意图；
图7为采用现有多次波压制技术获得的多次波压制后数据及多次波模型间的互相关结果示意图；
图8为采用本申请实施例的拉东域多次波压制方法获得的多次波压制后数据及多次波模型间的互相关结果示意图；
图9为本申请一实施例的拉东域多次波压制装置的结构框图；
图10为本申请另一实施例的拉东域多次波压制装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
参考图1所示，本申请实施例的拉东域多次波压制方法可以包括以下步骤：
S101、对去噪后的原始道集进行速度分析，拾取第一一次波速度。
在本申请一些实施方式中，由于多次波的干扰，拾取第一一次波速度(即一次波速度V0)可能并不准确。在本申请一示例性实施方式中，所述原始道集例如可以为共中心点(Common Middle Point，CMP)道集等。
S102、用所述第一一次波速度对所述原始道集进行动校正。
在本申请一些实施方式中，利用一次波速度V0对原始道集进行动校正时可以容忍较高的误差，即速度精度要求不高。
S103、将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域。
在本申请一些实施方式中，通过拉东变换可以将动校正后获得的道集中的一次波区域及多次波区域区分开来，以便于后续处理。
S104、将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集。
在本申请一些实施方式中，为了多次波能最少残留，可尽最大限度去除多次波区域，而不用担心有效信息的泄露，因为在后续的FK变换(即频率波数变换)提纯过程中，还可以将泄露的有效信息回填。而在拉东域切除掉多次波能量后，就可以完全压制多次波了。因此，获得的新道集中并不含多次波能量。在本申请一示例性实施方式中，所述将多次波区域置零可以是仅将所述多次波区域中其能量值大于零的部分置零(例如图2所示)。在图2中，黑竖线为Q值等于0，多次波的能量集中在Q大于0的区域，切除多次波时有两种选择：可以将大于0的区域能量全部充零；或者与黑竖线给一个小的斜坡缓冲即白线右边的能量全部充0。
S105、将所述原始道集与所述新道集进行相减，获得差道集。
在本申请一些实施方式中，如上所述，虽然上述新道集中并不含多次波能量，但由于多次波的存在导致获得的一次波速度V0有较大误差，即由于多次波的存在导致了一次波损失。因此，为了恢复一次波的损失，需要将所述原始道集与所述新道集进行相减，如此获得的差道集中就含有完整的一次波信息了。
S106、对所述新道集进行速度分析，获得第二一次波速度。
在本申请一些实施方式中，虽然多次波的存在导致了一次波损失，但得益于新道集不再有多次波的干扰，后续对这个新道集上进行速度分析，可就以得到更为准确的一次波速度V1。
S107、用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正。
S108、对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集。
如图3所示，在本申请一示例性实施方式中，一次波区域都集中在K＝0区域，提取多次波时，根据FK谱里的能量分布，将不集中在0的能量提取出来(如图中白色三角形区域)。
在本申请一些实施方式中，基于上述步骤S107和S108，即通过一次波速度V1，对差道集重新动校正，再通过FK变换将差道集变换到FK域，在FK域利用一次波和多次波能量分布在不同的区域的特点，可再次提取多次波模型，这个多次波模型相对初始多次波模型更纯，这种方法充分利用了拉东域多次波压制和FK域多次波压制技术两者的优势，在拉东域可以彻底切除多次波能量区域，不用担心损失有效波，因为在FK域提纯过程中，可以“回加”之前损失的一次有效波，即达到了在精准压制多次波的同时，又不损失有效信号的目的，因此，本申请实施例的拉东域多次波压制方法相对常规的拉东域去多次波更有优势。
S109、将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
本申请一些实施方式中，将所述原始道集与所述多次波道集相减，就可以得到压制多次波后的道集数据。在本申请一示例性实施方式中，所述将所述原始道集与所述多次波道集相减例如可以是将所述原始道集与所述多次波道集自适应相减，具体的：可先利用最小平方滤波法等对所述原始道集与所述多次波道集进行相位和振幅的匹配；然后将匹配后获得的原始道集与多次波道集进行相减。
为验证本申请实施例的效果，下面以图4所示的某工区的原始地震资料道集为例。采用现有多次波压制技术(即常规的拉东变换)处理图4所示原始地震资料道集，可得到如图5所示的道集。对比图4和图5可以看出：图5所示的道集在远偏(图5中上行箭头所指部分)依然有多次波残留，且近偏(图5中下行箭头所指部分)有同相轴泄露。而采用本申请实施例的拉东域多次波压制方法处理图4所示原始地震资料道，可得到如图6所示的道集。对比图5和图6可以看出：图6所示的道集的远偏并没有多次波残留，而且不存在同相轴泄露。
为了进一步验证本申请实施例的效果。可将采用现有多次波压制技术获得的多次波压制后数据及多次波模型进行互相关，得到如图7所示的互相关结果；并将采用本申请实施例的拉东域多次波压制方法获得的多次波压制后数据及多次波模型间进行互相关，得到如图8所示的互相关结果。从图7中可以看出，采用现有多次波压制技术获得的多次波模型及压制后数据间仍有较高的相关性，这就意味其存在多次波残留、有效信号泄露和/或有效反射损失等问题。而从图8中可以看出，采用本申请实施例的拉东域多次波压制方法获得的多次波压制后数据及多次波模型间已基本无相关性，表明多次波被较好的压制，同时有效信息也得到有效保护。
虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
参考图9所示，本申请实施例的拉东域多次波压制装置可以包括：
第一波速拾取模块91，可以用于对去噪后的原始道集进行速度分析，拾取第一一次波速度；
第一校正模块92，可以用于用所述第一一次波速度对所述原始道集进行动校正；
拉东变换模块93，可以用于将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域；
新道集获取模块94，可以用于将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集；
差道集获取模块95，可以用于将所述原始道集与所述新道集进行相减，获得差道集；
第二波速拾取模块96，可以用于对所述新道集进行速度分析，获得第二一次波速度；
第二校正模块97，可以用于用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正；
多次波道集获取模块98，可以用于对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集；
多次波道集压制模块99，可以用于将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
本申请实施例的装置与上述实施例的方法对应，因此，有关于本申请的装置细节，请参见上述实施例的方法，在此不再赘述。
参考图10所示，本申请实施例的拉东域多次波压制装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：
对去噪后的原始道集进行速度分析，拾取第一一次波速度；
用所述第一一次波速度对所述原始道集进行动校正；
将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域；
将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集；
将所述原始道集与所述新道集进行相减，获得差道集；
对所述新道集进行速度分析，获得第二一次波速度；
用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正；
对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集；
将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
本申请实施例的装置与上述实施例的方法对应，因此，有关于本申请的装置细节，请参见上述实施例的方法，在此不再赘述。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (11)
1.一种拉东域多次波压制方法，其特征在于，包括：
对去噪后的原始道集进行速度分析，拾取第一一次波速度；
用所述第一一次波速度对所述原始道集进行动校正；
将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域；
将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集；
将所述原始道集与所述新道集进行相减，获得差道集；
对所述新道集进行速度分析，获得第二一次波速度；
用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正；
对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集；
将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
2.如权利要求1所述的拉东域多次波压制方法，其特征在于，所述原始道集包括CMP道集。
3.如权利要求1所述的拉东域多次波压制方法，其特征在于，所述将所述多次波区域的能量置零，包括：
将所述多次波区域中其能量值大于零的部分置零。
4.如权利要求1所述的拉东域多次波压制方法，其特征在于，所述将所述原始道集与所述多次波道集相减，包括：
对所述原始道集与所述多次波道集进行相位和振幅的匹配；
将匹配后获得的原始道集与多次波道集进行相减。
5.如权利要求4所述的拉东域多次波压制方法，其特征在于，所述对所述原始道集与所述多次波道集进行相位和振幅的匹配，包括：
利用最小平方滤波法对所述原始道集与所述多次波道集进行相位和振幅的匹配。
6.一种拉东域多次波压制装置，其特征在于，包括：
第一波速拾取模块，用于对去噪后的原始道集进行速度分析，拾取第一一次波速度；
第一校正模块，用于用所述第一一次波速度对所述原始道集进行动校正；
拉东变换模块，用于将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域；
新道集获取模块，用于将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集；
差道集获取模块，用于将所述原始道集与所述新道集进行相减，获得差道集；
第二波速拾取模块，用于对所述新道集进行速度分析，获得第二一次波速度；
第二校正模块，用于用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正；
多次波道集获取模块，用于对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集；
多次波道集压制模块，用于将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
7.如权利要求6所述的拉东域多次波压制装置，其特征在于，所述原始道集包括CMP道集。
8.如权利要求6所述的拉东域多次波压制装置，其特征在于，所述将所述多次波区域的能量置零，包括：
将所述多次波区域中其能量值大于零的部分置零。
9.如权利要求6所述的拉东域多次波压制装置，其特征在于，所述将所述原始道集与所述多次波道集相减，包括：
对所述原始道集与所述多次波道集进行相位和振幅的匹配；
将匹配后获得的原始道集与多次波道集进行相减。
10.如权利要求9所述的拉东域多次波压制装置，其特征在于，所述对所述原始道集与所述多次波道集进行相位和振幅的匹配，包括：
利用最小平方滤波法对所述原始道集与所述多次波道集进行相位和振幅的匹配。
11.一种拉东域多次波压制装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，其特征在于，所述计算机程序被所述处理器运行时执行如下步骤：
对去噪后的原始道集进行速度分析，拾取第一一次波速度；
用所述第一一次波速度对所述原始道集进行动校正；
将动校正后获得的道集进行拉东变换，获得拉东域的一次波及多次波区域；
将所述多次波区域的能量置零，并将所述一次波区域及能量置零后的多次波区域反变换至TX域，获得新道集；
将所述原始道集与所述新道集进行相减，获得差道集；
对所述新道集进行速度分析，获得第二一次波速度；
用所述第一一次波速度对所述差道集进行反动校正，并用所述第二一次波速度对反动校正后的差道集进行动校正；
对动校正后的差道集进行FK变换，以获取多次波能量区域，并将该多次波能量区域反变换至TX域，获得多次波道集；
将所述原始道集与所述多次波道集相减，获得压制多次波后的道集。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710945351.3A CN107656313B (zh) | 2017-10-12 | 2017-10-12 | 一种拉东域多次波压制方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710945351.3A CN107656313B (zh) | 2017-10-12 | 2017-10-12 | 一种拉东域多次波压制方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107656313A CN107656313A (zh) | 2018-02-02 |
| CN107656313B true CN107656313B (zh) | 2019-08-06 |
Family
ID=61117181
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710945351.3A Active CN107656313B (zh) | 2017-10-12 | 2017-10-12 | 一种拉东域多次波压制方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107656313B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109001813B (zh) * | 2018-07-16 | 2020-01-07 | 中国石油天然气股份有限公司 | 一种压制多次波的方法、装置及系统 |
| CN111562616B (zh) * | 2019-02-14 | 2023-06-30 | 中国石油天然气股份有限公司 | 地震数据散射噪音压制方法及装置 |
| CN111596358B (zh) * | 2019-02-21 | 2023-04-07 | 中国石油天然气股份有限公司 | 一种多次波的压制方法、设备以及系统 |
| CN114428337A (zh) * | 2020-09-28 | 2022-05-03 | 中国石油化工股份有限公司 | 一种拉冬变换多次波压制方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6094620A (en) * | 1999-01-08 | 2000-07-25 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
| WO2005019869A1 (en) * | 2003-08-15 | 2005-03-03 | Westerngeco, L.L.C. | Method for attenuating water layer multiples |
| CN102401908A (zh) * | 2010-09-07 | 2012-04-04 | 中国石油天然气集团公司 | 一种利用不同模加权稀疏的抛物拉东变换压制多次波的方法 |
| CN103869364A (zh) * | 2014-03-25 | 2014-06-18 | 中国石油大学(华东) | 一种基于双项抛物拉东变换的多次波压制方法 |
-
2017
- 2017-10-12 CN CN201710945351.3A patent/CN107656313B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6094620A (en) * | 1999-01-08 | 2000-07-25 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
| WO2005019869A1 (en) * | 2003-08-15 | 2005-03-03 | Westerngeco, L.L.C. | Method for attenuating water layer multiples |
| CN102401908A (zh) * | 2010-09-07 | 2012-04-04 | 中国石油天然气集团公司 | 一种利用不同模加权稀疏的抛物拉东变换压制多次波的方法 |
| CN103869364A (zh) * | 2014-03-25 | 2014-06-18 | 中国石油大学(华东) | 一种基于双项抛物拉东变换的多次波压制方法 |
Non-Patent Citations (5)
| Title |
|---|
| multiple attenuation:coping with the spatial truncation effect in the radon transform domain;yanghua wang;《geophysical prospecting》;20031231;第75-87页 |
| Suppression of mUltiple reflections using the Radon transform;Douglas J. Foster et al.;《GEOPHYSICS》;19920331;第386-395页 |
| 基于F-K和Radon变换的多次波衰减方法;张亚伟，等;《中国地质调查》;20150331;第22-27页 |
| 线性拉东域预测反褶积在海洋多次波去除中的应用;赵昌垒，等;《地球物理学进展》;20130430;第1026-1032页 |
| 迭代抛物Radon变换法分离一次波与多次波;谢俊法，等;《石油地球物理勘探》;20140228;第76-82页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107656313A (zh) | 2018-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107656313B (zh) | 一种拉东域多次波压制方法及装置 | |
| CN108549106B (zh) | 混叠噪声压制方法及装置 | |
| CN106855640B (zh) | 一种海底电缆地震数据初至计算方法及装置 | |
| CN101877237A (zh) | 降噪装置和噪声确定方法 | |
| CN109143343B (zh) | 绕射波成像方法、装置和电子设备 | |
| CN107894613A (zh) | 弹性波矢量成像方法、装置、存储介质及设备 | |
| CN109101615A (zh) | 地震勘探数据处理方法及装置 | |
| CN110380985A (zh) | 基于交易链路的流量控制方法、装置、设备及存储介质 | |
| CN115267912A (zh) | 基于压缩感知的avo特征保持去噪方法及系统 | |
| CN109541682A (zh) | 各向同性弹性参数保幅反演方法及装置 | |
| CN105068119B (zh) | 低频地震数据中的面波的衰减方法及其装置 | |
| CN117214950B (zh) | 一种多次波压制方法、装置、设备及存储介质 | |
| CN105093328A (zh) | 一种滑动扫描谐波压制方法及装置 | |
| CN107255833B (zh) | 地震叠加剖面的确定方法和装置 | |
| CN108363103B (zh) | 规则干扰压制方法及装置 | |
| CN109856672B (zh) | 基于深度波数谱的瞬变波包提取方法、存储介质与终端 | |
| CN106094020A (zh) | 一种地震反演方法及装置 | |
| CN110163407A (zh) | 量化策略的优化方法及装置 | |
| CN107991704B (zh) | 自由表面多次波的压制方法和装置 | |
| CN107229073A (zh) | 地震数据处理方法和装置 | |
| CN112014884B (zh) | 压制近炮点强能量噪声的方法及装置 | |
| CN108845356B (zh) | 绕射波成像方法、装置和电子设备 | |
| CN111856557A (zh) | 深度域合成地震记录制作方法及装置 | |
| CN108900003A (zh) | 一种配电台区监控系统 | |
| CN111781644A (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 |