CN107991704A - 自由表面多次波的压制方法和装置 - Google Patents
自由表面多次波的压制方法和装置 Download PDFInfo
- Publication number
- CN107991704A CN107991704A CN201710778207.5A CN201710778207A CN107991704A CN 107991704 A CN107991704 A CN 107991704A CN 201710778207 A CN201710778207 A CN 201710778207A CN 107991704 A CN107991704 A CN 107991704A
- Authority
- CN
- China
- Prior art keywords
- free surface
- wave
- mrow
- multiple wave
- reflective portion
- 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 62
- 238000000605 extraction Methods 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 4
- 241001269238 Data Species 0.000 claims description 2
- 230000007613 environmental effect Effects 0.000 claims description 2
- 239000000047 product Substances 0.000 description 19
- 238000010586 diagram Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 238000003825 pressing Methods 0.000 description 4
- 238000000926 separation method Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000000877 morphologic effect Effects 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 241000894007 species Species 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
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)
- Aerials With Secondary Devices (AREA)
Abstract
本发明提供了一种自由表面多次波的压制方法和装置，其中，该方法包括：将所获取的第一单炮地震数据划分为反射部分以及非反射部分；提取反射部分中的第一自由表面多次波；提取非反射部分中的第二自由表面多次波；根据第一单炮地震数据、第一自由表面多次波以及第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。在本发明实施例中，所得到的压制自由表面多次波之后的第二单炮地震数据精度更高。
Description
技术领域
本发明涉及地质勘探技术领域，特别涉及一种自由表面多次波的压制方法和装置。
背景技术
在地震资料处理过程中，多次波通常被视为一种相干噪声，特别是针对海上采集的地震数据，多次波发育尤为严重，处理过程中，如果不能采取有效的手段对多次波予以消除，最终的成像结果往往会造成成果数据出现地质假象，甚至会降低地震资料信噪比以及分辨率。多次波根据其产生机理，其种类非常复杂，针对海洋资料二样，能量最强的往往是自由表面多次波。
通常，通过采用多次波压制方案，可以有效的消除自由表面反射率对地震数据的影响。目前，常规多次波压制方法通常只需要提供震源处地震资料的波形特征及自由表面反射率，并不需要提取详细的地下地质信息，而随着地震勘探技术的不断进步，自由表面多次波压制方法也逐渐演变成了适应性的反演处问题，即：首先，估算震源与自由表面的反射率；其次，将不含多次的数据视为反演剩余量。但是此方法受计算成本及精度的限制，无法应用于三维地震数据。
发明内容
本发明提供了一种自由表面多次波的压制方法和装置，以解决现有技术无法满足三维地震数据压制多次波的问题。
本发明实施例提供了一种自由表面多次波的压制方法，可以包括：将所获取的第一单炮地震数据划分为反射部分以及非反射部分；提取所述反射部分中的第一自由表面多次波；提取所述非反射部分中的第二自由表面多次波；根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。
在一个实施例中，所述提取所述反射部分中的第一自由表面多次波包括：获取所述反射部分中入射波为反射波、出射波也为反射波的第三自由表面多次波；获取所述反射部分中入射波为反射波、出射波为非反射波的第四自由表面多次波；将所述第三自由表面多次波以及所述第四自由表面多次波的和作为所述反射波中的第一自由表面多次波。
在一个实施例中，按照以下公式提取所述反射波中的第一自由表面多次波：
其中，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第三自由表面多次波，表示第n次迭代后得到的第四自由表面多次波。
在一个实施例中，所述提取所述反射部分中的第一自由表面多次波包括：获取所述反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述反射部分的乘积，将所得到的乘积作为所述第三自由表面多次波；计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第四自由表面多次波；计算所述第三自由表面多次波以及所述第四自由表面多次波的和，将所得到的和作为所述反射部分中的第一自由表面多次波。
在一个实施例中，按照以下公式提取所述反射部分中的第一自由表面多次波：
其中，表示第n次迭代后得到的第一自由表面多次波，表示所述反射部分中第n次迭代去除自由表面多次波后的波形数据，Dr表示所述反射部分，Ds表示所述非反射部分。
在一个实施例中，所述提取所述非反射部分中的第二自由表面多次波包括：获取所述非反射部分中入射波为非反射波、出射波也为非反射波的第五自由表面多次波；获取所述非反射部分中入射波为非反射波、出射波为反射波的第六自由表面多次波；将所述第五自由表面多次波以及所述第六自由表面多次波的和作为所述非反射波中的第二自由表面多次波。
在一个实施例中，按照以下公式提取所述非反射部分中的第二自由表面多次波：
其中，表示第n次迭代后得到的第二自由表面多次波，表示第n次迭代后得到的第五自由表面多次波，表示第n次迭代后得到的第六自由表面多次波。
在一个实施例中，所述提取所述非反射部分中的第二自由表面多次波包括：获取所述非反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；计算所述非反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第五自由表面多次波；计算所述非反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第六自由表面多次波；计算所述第五自由表面多次波以及所述第六自由表面多次波的和，将所得到的和作为所述非反射部分中的第二自由表面多次波。
在一个实施例中，按照以下公式提取所述非反射部分中的第二自由表面多次波：
其中，表示第n次迭代后得到的第二自由表面多次波，表示所述非反射部分中第n次迭代去除自由表面多次波后的波形数据，Dr表示所述反射部分，Ds表示所述非反射部分。
在一个实施例中，按照以下公式计算得到压制自由表面多次波之后的第二单炮地震数据：
上式中，
上式中，表示第n次迭代压制自由表面多次波之后的第二单炮地震数据，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第二自由表面多次波，A表示表面算子，D表示所述第一单炮地震数据。
在一个实施例中，根据所述第一单炮地震数据的同相轴形态差异，将所获取的第一单炮地震数据划分为反射部分以及非反射部分。
在一个实施例中，在将所获取的第一单炮地震数据划分为反射部分以及非反射部分之前，所述方法还包括：对所述第一单炮地震数据进行预处理，所述预处理包括以下至少之一：环境校正、静校正、去噪、动校正拉伸。
在一个实施例中，在计算得到压制自由表面多次波之后的第二单炮地震数据之后，所述方法还包括：根据所述第一单炮地震数据的叠加速度，绘制所述第二单炮地震数据的叠加剖面。
本发明实施例还提供了一种自由表面多次波的压制装置，可以包括：第一数据划分模块，用于将所获取的第一单炮地震数据划分为反射部分以及非反射部分；第一波形提取模块，用于提取所述反射部分中的第一自由表面多次波；第二波形提取模块，用于提取所述非反射部分中的第二自由表面多次波；第二数据划分模块，用于根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。
在本发明实施例中，先将所获取的第一单炮地震数据划分为反射部分及非反射部分等两部分；并提取反射部分中所含第一自由表面多次波以及非反射部分中的第二自由表面多次波；最后，可以根据单炮地震数据、第一自由表面多次波以及第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。由于本申请所提出的方法先将第一单炮地震数据中的反射部分和非反射部分进行分离，再分别从反射部分及非反射部分中提取自由表面多次波，解决了现有技术中自由表面多次波压制法仅可有效预测二维地震数据反射部分多次波，而无法识别分辨三维地震数据中的非反射部分所含的多次波的问题。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的一种自由表面多次波的压制方法流程图；
图2是本申请提供的某地区原始单炮数据示意图；
图3是采用现有技术进行自由表面多次波压制后所得到的地震数据示意图；
图4是采用本申请所提出的自由表面多次波的压制方法所得到的地震数据示意图；
图5是采用现有技术进行自由表面多次波压制后所得到的叠加剖面示意图；
图6是采用本申请所提出的自由表面多次波的压制方法所得到的叠加剖面示意图；
图7是本申请提供的一种自由表面多次波的压制装置的一种结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
需要说明的是，在本申请的描述中，术语“第一”、“第二”等仅用于描述目的和区别类似的对象，两者之间并不存在先后顺序，也不能理解为指示或暗示相对重要性。此外，在本申请的描述中，除非另有说明，“多个”的含义是两个或两个以上。
现有技术中对自由表面多次波进行压制时，由于常规的自由表面多次波压制方法仅仅能识别二维地震数据中反射部分中的多次波；但是由非反射部分引起的多次波很难被预测，因而会限制整个方法对三维资料的的应用效果。针对上述现有技术中二维自由表面多次波压制法基本可以消除反射部分的多次波，但是由非反射部分引起的多次波很难被预测的缺陷，发明人提出了将第一单炮地震数据中体现三维效应最强的反射部分与非反射部分区分开来，并分别提取反射部分中的第一自由表面多次波以及非反射部分中的第二自由表面多次波，根据所提取的第一自由表面多次波以及非反射部分中的第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮数据。基于此，提出了一种自由表面多次波的压制方法，如图1所示，可以包括以下步骤：
S101：将所获取的第一单炮地震数据划分为反射部分以及非反射部分。
在本实施例中，可以利用野外地震队所提供的野外采集的原始数据，也就是第一单炮数据，并将所获取的第一单炮地震数据划分为反射部分以及非反射部分。
其中，上述第一单炮地震数据指的是二维或者三维的单炮地震数据，将此处的单炮地震数据描述为第一单炮地震数据是为了和下文中已经压制自由表面多次波之后所得到的第二单炮地震数据进行区分。此处的第一单炮地震数据以及下面的第二单炮地震数据均可以指的是单炮地震数据，只是描述的是不同方式获取的单炮地震数据，并不具有其它特殊含义。
进一步地，在本实施例中，可以利用第一单炮地震数据的同相轴形态差异，将反射部分以及非反射部分进行分离处理。通过分别去除反射部分以及非反射部分中的自由表面多次波，从而可以避免现有技术中二维自由表面多次波压制法无法预测三维地震数据中由于非反射部分所引起的多次波的缺陷。
S102：提取反射部分中的第一自由表面多次波。
S103：提取非反射部分中的第二自由表面多次波。
在本实施例中，由于在S101中已经将第一单炮地震数据中的反射部分以及非反射部分分离开来，因而，多次波的预测可以分别在不同的数据子集对反射部分以及非反射部分这两个部分进行预测和消除，可以分别提取反射部分中的第一自由表面多次波以及非反射部分中的第二自由表面多次波。值得的注意的是，上述第一自由表面多次波以及第二自由表面多次波以及下文的第三自由表面多次波、第四自由表面多次波、第五自由表面多次波以及第六自由表面多次波均可以指的是自由表面多次波。具体的，可以将反射部分中所含多次波命名为第一自由表面多次波；非反射部分中所含多次波命名为第二自由表面多次波；将入射部分为反射波、出射部分也为反射波(简称为R-R)命名为第三自由表面多次波；将入射部分为反射波、出射部分为非反射波(简称为R-S)命名为第四自由表面多次波；将入射部分为非反射波、出射部分也为非反射波(简称为S-S)命名为第五自由表面多次波；将入射部分为非反射波、出射部分为反射波(简称为S-R)命名为第六自由表面多次波。
此处仅是为了区分这些自由表面多次波均为不同方式获取的，这些参数之间并不存在先后顺序，也不能理解为指示或暗示相对重要性。
在本申请的一个实施例中，提取所述反射部分中的第一自由表面多次波可以包括：获取反射部分中入射波为反射波、出射波也为反射波的第三自由表面多次波；获取反射部分中入射波为反射波、出射波为非反射波的第四自由表面多次波；将第三自由表面多次波以及第四自由表面多次波的和作为反射波中的第一自由表面多次波。
其中，获取反射部分中入射波为反射波、出射波也为反射波的第三自由表面多次波，可以包括：获取反射部分中第n次迭代去除自由表面多次波后的波形数据；计算反射部分中第n次迭代去除自由表面多次波后的波形数据与反射部分的乘积，将所得到的乘积作为第三自由表面多次波。具体的，可以按照以下公式计算得到该第三自由表面多次波：
其中，表示第n次迭代后得到的第三自由表面多次波，表示反射波中第n次迭代去除自由表面多次波后的波形数据，Dr表示反射部分。
相应的，获取反射部分中入射波为反射波、出射波为非反射波的第四自由表面多次波，可以包括：计算反射部分中第n次迭代去除自由表面多次波后的波形数据与非反射部分的乘积，将所得到的乘积作为第四自由表面多次波。具体的，可以按照以下公式计算得到该第四自由表面多次波：
其中，表示第n次迭代后得到的第四自由表面多次波，表示反射波中第n次迭代去除自由表面多次波后的波形数据，Ds表示非反射部分。
上述的初始值，即可以是还没进行迭代时的值，即，分离后所得到的一次反射波及里面包含的多次波。上述迭代次数n，可以根据计算机的计算能力和多次波最终的预测程度进行一个大体的估算。或者可以提供一个预设的期望阈值，如果能够达到该期望阈值便停止迭代。实际应用过程中，我们通常会根据计算机的计算能力，人为给一个具体的迭代次数，如果计算机运算能力超强，像目前计算机技术中发展较为成熟的图形处理技术发展迅速，在此计算能力的保障之下，我们可以多迭代几次，只要每次迭代都能实现有效收敛，则压制多次波的精度，便会与迭代次数成正比；然而，目前地震数据都是海量运算，我们基本都会将有效迭代次数控制在10次以内，以便实现计算效果和计算效率的有效平衡。
在得到上述第三自由表面多次波以及第四自由表面多次波之后，可以计算第三自由表面多次波以及第四自由表面多次波的和，并将所得到的和作为反射波中的第一自由表面多次波。具体的，可以按照以下公式计算得到该第一自由表面多次波：
其中，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第三自由表面多次波，表示第n次迭代后得到的第四自由表面多次波。
在本申请的一个实施例中，提取所述非反射部分中的第二自由表面多次波可以包括：获取所述非反射部分中入射波为非反射波、出射波也为非反射波的第五自由表面多次波；获取所述非反射部分中入射波为非反射波、出射波为反射波的第六自由表面多次波；将所述第五自由表面多次波以及所述第六自由表面多次波的和作为所述非反射波中的第二自由表面多次波。
其中，获取非反射部分入射部分为非反射波、出射部分也为非反射波的第五自由表面多次波，可以包括：获取非反射部分中第n次迭代去除自由表面多次波后的波形数据；计算非反射部分中第n次迭代去除自由表面多次波后的波形数据与非反射部分的乘积，将所得到的乘积作为第五自由表面多次波。具体的，可以按照以下公式计算得到该第五自由表面多次波：
其中，表示第n次迭代后得到的第五自由表面多次波，表示非反射部分中第n次迭代去除自由表面多次波后的波形数据，Ds表示非反射部分。
相应的，获取非反射部分中入射部分为非反射波、出射波为反射波的第六自由表面多次波，可以包括：计算非反射部分中第n次迭代去除自由表面多次波后的波形数据与非反射波的乘积，将所得到的乘积作为第六自由表面多次波。具体的，可以按照以下公式计算得到该第六自由表面多次波：
其中，表示第n次迭代后得到的第六自由表面多次波，表示非反射波中第n次迭代去除自由表面多次波后的波形数据，Dr表示反射部分。
上述的初始值，即可以是还没进行迭代时的值，即，分离后所得到的一次反射波及里面包含的多次波。上述迭代次数n，可以根据计算机的计算能力和多次波最终的预测程度进行一个大体的估算。如果每次迭代都能实现有效收敛，则压制多次波的精度，便会与迭代次数成正比；然而，目前地震数据都是海量运算，我们基本都会将有效迭代次数控制在10次以内，以便实现计算效果和计算效率的有效平衡。
在得到上述第五自由表面多次波以及第六自由表面多次波之后，可以计算第五自由表面多次波以及第六自由表面多次波的和，并将所得到的和作为反射波中的第二自由表面多次波。具体的，可以按照以下公式计算得到该第二自由表面多次波：
其中，表示第n次迭代后得到的第二自由表面多次波，表示第n次迭代后得到的第五自由表面多次波，表示第n次迭代后得到的第六自由表面多次波。
S104：根据第一单炮地震数据、第一自由表面多次波以及第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。
在本实施例中，采用S102以及S103中的方法分别计算得到反射波中的第一自由表面多次波以及非反射部分中的第二自由表面多次波，并根据上述单炮地震数据、第一自由表面多次波以及第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。
具体的，在本实施例中，可以按照以下公式计算得到压制自由表面多次波之后的第二单炮地震数据：
上式中，
上式中，表示第n次迭代压制自由表面多次波之后的第二单炮地震数据，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第二自由表面多次波，A表示表面算子，D表示所述第一单炮地震数据。
其中，A为与角频率相关的一个标量，在预测到自由表面多次波之后，可以基于一次波能量最小的原理计算得到。
通过上述方式计算得到不同多次波：第三自由表面多次波、第四自由表面多次波、第五自由表面多次波、第六自由表面多次波，都可以通过与其相适应的处理方法进行滤除，利用这种新的多次波压制思路，非反射部分可以用不同于反射波的滤波器减去，这样就可以处理复杂的三维效应了。
在计算得到压制自由表面多次波之后的第二单炮地震数据之后，可以根据第一单炮地震数据的叠加速度，绘制第二单炮地震数据的叠加剖面。
因此，本申请达到了以下目的：
1)有效实现了对三维自由表面多次波的有效压制，但又避免了全三维应用。本发明尝试利用平面波记录中反射部分和非反射部分的同相轴形态差异(分别为似双曲线和似线性)，模拟平面波激发的记录，将体现三维效应最强的反射部分与非反射部分的分离。通过将地震数据分离成反射部分和非反射部分，多次波的预测就可以在输入数据的这些不同子集上进行。
2)基于波场分离理论的自由表面多次波的压制方法，可以得到四项预测多次波：R-R、R-S、S-R和S-S预测多次波。因为每一种预测多次波都能用不同的适应滤波器从输入数据中减去，故得到了更灵活的自由表面多次波的压制方法。
3)由本发明预测的不同子集的多次波，都可以通过与其相适应的处理被减去。利用这种新的多次波压制思路，非反射部分可以用不同于反射波的滤波器减去，这样就可以实现对三维自由表面多次波的有效压制，从而适应复杂的三维效应。
为了测试本发明的准确性和稳定性，同时，为了测试本发明的普适性，特别是针对野外实际采集资料的有效性，选取地球物理界内最典型、最复杂的西部某复杂山地区实际地震数据进行测试。具体按照以下实现步骤：
(1)加载西部该复杂山地区地震采集的单炮数据，如图2所示为该地区的原始单炮数据示意图；
(2)对加载单炮数据进行静校正、去噪等叠前预处理工作；
(3)提供相对准确的叠加速度；
(4)将完成步骤(2)的地震数据进行分离，分离出其中所包含的反射部分和非反射部分；
(5)将步骤(4)中分离出来的数据进行多次波的预测，分别预测出R-R、R-S、S-R和S-S等四种类型的多次波；
(6)通过自适应相减的思路，在不同的数据子集消除由步骤(5)所预测出来的四种多次波，便可以实现消除了多次波的三维单炮数据，如图3所示为采用现有技术进行自由表面多次波压制后所得到的地震数据示意图，如图4所示为采用本申请所提出的自由表面多次波的压制方法所得到的地震数据示意图，将图2、图3和图4的划框部分进行对比，可以看出图4中黑色连续双曲线型的同相轴比图3中少很多，图3中黑色连续双曲线型的同相轴比图2中也少很多，且能量也弱很多(就是颜色白一些)，说明图4相较于图3有了较大改进。
(7)根据步骤(3)中所得到的叠加速度，将步骤(6)中所得到的单炮数据转换成叠加剖面。如图5所示为采用现有技术进行自由表面多次波压制后所得到的叠加剖面示意图，如图6所示为采用本申请所提出的自由表面多次波的压制方法所得到的叠加剖面示意图。对比图6和图5中的划框部分可知：图6中划框部分噪声较少，且颜色较淡，说明图6相较于图5有了较大改进。即，采用本申请所提出的自由表面多次波的压制方法所得到的叠加剖面更加精确。
基于同一发明构思，本发明实施例中还提供了一种自由表面多次波的压制装置，如下面的实施例所述。由于自由表面多次波的压制装置解决问题的原理与自由表面多次波的压制方法相似，因此自由表面多次波的压制装置的实施可以参见自由表面多次波的压制方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图7是本发明实施例的自由表面多次波的压制装置的一种结构框图，如图7所示，可以包括：第一数据划分模块701、第一波形提取模块702、第二波形提取模块703、第二数据划分模块704，下面对该结构进行说明。
第一数据划分模块701，可以用于将所获取的第一单炮地震数据划分为反射部分以及非反射部分；
第一波形提取模块702，可以用于提取所述反射部分中的第一自由表面多次波；
第二波形提取模块703，可以用于提取所述非反射部分中的第二自由表面多次波；
第二数据划分模块704，可以用于根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。
在本发明实施例中，先将所获取的第一单炮地震数据划分为反射部分及非反射部分等两部分；并提取反射部分中所含第一自由表面多次波以及非反射部分中的第二自由表面多次波；最后，可以根据单炮地震数据、第一自由表面多次波以及第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。由于本申请所提出的方法先将第一单炮地震数据中的反射部分和非反射部分进行分离，再分别从非反射部分及反射部分中提取自由表面多次波，解决了现有技术中自由表面多次波压制法仅可有效预测二维地震数据反射部分多次波，而无法识别分辨三维地震数据中的非反射部分所含的多次波的问题。
利用上述各实施例所提供的自由表面多次波的压制装置的实施方式，可以自动实施所述自由表面多次波的压制方法，对自由表面多次波进行压制处理，可以不需要实施人员的具体参与，可以直接输出自由表面多次波的压制结果，操作简单快捷，有效提高了用户体验。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本申请并不局限于必须是符合本申请实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取/存储/计算等获取的实施例，仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的单元、装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (14)
1.一种自由表面多次波的压制方法，其特征在于，包括：
将所获取的第一单炮地震数据划分为反射部分以及非反射部分；
提取所述反射部分中的第一自由表面多次波；
提取所述非反射部分中的第二自由表面多次波；
根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。
2.如权利要求1所述的方法，其特征在于，所述提取所述反射部分中的第一自由表面多次波包括：
获取所述反射部分中入射波为反射波、出射波也为反射波的第三自由表面多次波；
获取所述反射部分中入射波为反射波、出射波为非反射波的第四自由表面多次波；
将所述第三自由表面多次波以及所述第四自由表面多次波的和作为所述反射波中的第一自由表面多次波。
3.如权利要求2所述的方法，其特征在于，按照以下公式提取所述反射波中的第一自由表面多次波：
<mrow>
<msubsup>
<mi>M</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>M</mi>
<mrow>
<mi>r</mi>
<mi>r</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mrow>
<mi>r</mi>
<mi>s</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
其中，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第三自由表面多次波，表示第n次迭代后得到的第四自由表面多次波。
4.如权利要求2所述的方法，其特征在于，所述提取所述反射部分中的第一自由表面多次波包括：
获取所述反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；
计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述反射部分的乘积，将所得到的乘积作为所述第三自由表面多次波；
计算所述反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第四自由表面多次波；
计算所述第三自由表面多次波以及所述第四自由表面多次波的和，将所得到的和作为所述反射部分中的第一自由表面多次波。
5.如权利要求4所述的方法，其特征在于，按照以下公式提取所述反射部分中的第一自由表面多次波：
<mrow>
<msubsup>
<mi>M</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>D</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>r</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>D</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>r</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mi>s</mi>
</msub>
</mrow>
其中，表示第n次迭代后得到的第一自由表面多次波，表示所述反射部分中第n次迭代去除自由表面多次波后的波形数据，Dr表示所述反射部分，Ds表示所述非反射部分。
6.如权利要求1所述的方法，其特征在于，所述提取所述非反射部分中的第二自由表面多次波包括：
获取所述非反射部分中入射波为非反射波、出射波也为非反射波的第五自由表面多次波；
获取所述非反射部分中入射波为非反射波、出射波为反射波的第六自由表面多次波；
将所述第五自由表面多次波以及所述第六自由表面多次波的和作为所述非反射波中的第二自由表面多次波。
7.如权利要求6所述的方法，其特征在于，按照以下公式提取所述非反射部分中的第二自由表面多次波：
<mrow>
<msubsup>
<mi>M</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>M</mi>
<mrow>
<mi>s</mi>
<mi>r</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>M</mi>
<mrow>
<mi>s</mi>
<mi>s</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
其中，表示第n次迭代后得到的第二自由表面多次波，表示第n次迭代后得到的第五自由表面多次波，表示第n次迭代后得到的第六自由表面多次波。
8.如权利要求6所述的方法，其特征在于，所述提取所述非反射部分中的第二自由表面多次波包括：
获取所述非反射部分中第n次迭代去除自由表面多次波后的波形数据，其中，n为正整数；
计算所述非反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第五自由表面多次波；
计算所述非反射部分中第n次迭代去除自由表面多次波后的波形数据与所述非反射部分的乘积，将所得到的乘积作为所述第六自由表面多次波；
计算所述第五自由表面多次波以及所述第六自由表面多次波的和，将所得到的和作为所述非反射部分中的第二自由表面多次波。
9.如权利要求8所述的方法，其特征在于，按照以下公式提取所述非反射部分中的第二自由表面多次波：
<mrow>
<msubsup>
<mi>M</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>D</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>s</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>D</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>s</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msub>
<mi>D</mi>
<mi>s</mi>
</msub>
</mrow>
其中，表示第n次迭代后得到的第二自由表面多次波，表示所述非反射部分中第n次迭代去除自由表面多次波后的波形数据，Dr表示所述反射部分，Ds表示所述非反射部分。
10.如权利要求1所述的方法，其特征在于，按照以下公式计算得到压制自由表面多次波之后的第二单炮地震数据：
<mrow>
<msubsup>
<mi>D</mi>
<mn>0</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mi>D</mi>
<mo>-</mo>
<msup>
<mi>AM</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
上式中，
上式中，表示第n次迭代压制自由表面多次波之后的第二单炮地震数据，表示第n次迭代后得到的第一自由表面多次波，表示第n次迭代后得到的第二自由表面多次波，A表示表面算子，D表示所述第一单炮地震数据。
11.如权利要求1所述的方法，其特征在于，根据所述第一单炮地震数据的同相轴形态差异，将所获取的第一单炮地震数据划分为反射部分以及非反射部分。
12.如权利要求1所述的方法，其特征在于，在将所获取的第一单炮地震数据划分为反射部分以及非反射部分之前，所述方法还包括：
对所述第一单炮地震数据进行预处理，所述预处理包括以下至少之一：环境校正、静校正、去噪、动校正拉伸。
13.如权利要求1所述的方法，其特征在于，在计算得到压制自由表面多次波之后的第二单炮地震数据之后，所述方法还包括：
根据所述第一单炮地震数据的叠加速度，绘制所述第二单炮地震数据的叠加剖面。
14.一种自由表面多次波的压制装置，其特征在于，包括：
第一数据划分模块，用于将所获取的第一单炮地震数据划分为反射部分以及非反射部分；
第一波形提取模块，用于提取所述反射部分中的第一自由表面多次波；
第二波形提取模块，用于提取所述非反射部分中的第二自由表面多次波；
第二数据划分模块，用于根据所述第一单炮地震数据、所述第一自由表面多次波以及所述第二自由表面多次波，计算得到压制自由表面多次波之后的第二单炮地震数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710778207.5A CN107991704B (zh) | 2017-09-01 | 2017-09-01 | 自由表面多次波的压制方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710778207.5A CN107991704B (zh) | 2017-09-01 | 2017-09-01 | 自由表面多次波的压制方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107991704A true CN107991704A (zh) | 2018-05-04 |
| CN107991704B CN107991704B (zh) | 2019-09-06 |
Family
ID=62028911
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710778207.5A Active CN107991704B (zh) | 2017-09-01 | 2017-09-01 | 自由表面多次波的压制方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107991704B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031412A (zh) * | 2018-07-24 | 2018-12-18 | 东北大学 | 一种弹性波被动源数据一次波估计方法 |
| CN109307886A (zh) * | 2018-08-24 | 2019-02-05 | 中国石油天然气股份有限公司 | 多次波自适应相减方法及装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050065758A1 (en) * | 2003-09-23 | 2005-03-24 | Ian Moore | Method for the 3-D prediction of free-surface multiples |
| CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
| CN103576198A (zh) * | 2012-08-02 | 2014-02-12 | 中国石油天然气集团公司 | 一种快速二维海上地震资料自由表面多次波预测方法 |
| WO2015144453A1 (en) * | 2014-03-24 | 2015-10-01 | Statoil Petroleum As | Removal of sea surface effects from seismic data |
| CN105259575A (zh) * | 2015-10-12 | 2016-01-20 | 中国石油大学(华东) | 快速3d自由表面多次波预测方法 |
-
2017
- 2017-09-01 CN CN201710778207.5A patent/CN107991704B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050065758A1 (en) * | 2003-09-23 | 2005-03-24 | Ian Moore | Method for the 3-D prediction of free-surface multiples |
| CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
| CN103576198A (zh) * | 2012-08-02 | 2014-02-12 | 中国石油天然气集团公司 | 一种快速二维海上地震资料自由表面多次波预测方法 |
| WO2015144453A1 (en) * | 2014-03-24 | 2015-10-01 | Statoil Petroleum As | Removal of sea surface effects from seismic data |
| CN105259575A (zh) * | 2015-10-12 | 2016-01-20 | 中国石油大学(华东) | 快速3d自由表面多次波预测方法 |
Non-Patent Citations (2)
| Title |
|---|
| 宋家文 等: "多次波压制的研究现状与进展", 《地球物理学进展》 * |
| 王通 等: "三维自由表面多次波去除方法", 《吉林大学学报（地球科学版）》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031412A (zh) * | 2018-07-24 | 2018-12-18 | 东北大学 | 一种弹性波被动源数据一次波估计方法 |
| CN109307886A (zh) * | 2018-08-24 | 2019-02-05 | 中国石油天然气股份有限公司 | 多次波自适应相减方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107991704B (zh) | 2019-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6782124B2 (en) | Three dimensional empirical mode decomposition analysis apparatus, method and article manufacture | |
| Hou et al. | Simultaneous multi-component seismic denoising and reconstruction via K-SVD | |
| CN107991704A (zh) | 自由表面多次波的压制方法和装置 | |
| CN109143343B (zh) | 绕射波成像方法、装置和电子设备 | |
| CN101739670B (zh) | 非局部均值空域时变图像滤波方法 | |
| US20100027848A1 (en) | Structure-independent analysis of 3-D seismic random noise | |
| AU2014342986B2 (en) | Method for post-stack noise mitigation | |
| CN111751870B (zh) | 叠后层间多次波压制方法及装置 | |
| US11662492B2 (en) | Seismic random noise attenuation | |
| Wu et al. | Attenuating seismic noise via incoherent dictionary learning | |
| CN117214950B (zh) | 一种多次波压制方法、装置、设备及存储介质 | |
| Oboué et al. | Erratic and random noise attenuation using adaptive local orthogonalization | |
| CN107065007B (zh) | 一种地震资料振幅调整方法和装置 | |
| CN107607995B (zh) | 虚反射的压制方法和装置 | |
| CN108693558B (zh) | 地震数据处理方法和装置 | |
| CN107703541B (zh) | 地层倾角的确定方法和装置 | |
| CN111527420A (zh) | 将正交化滤波应用于波场分离 | |
| CN114782666A (zh) | 识别模型训练方法及装置、车牌识别方法及装置 | |
| NO20210731A1 (en) | Fault Skeletonization For Fault Identification in a Subterranean Environment | |
| CN113687421B (zh) | 地震信号的数据处理方法、装置、电子设备及存储介质 | |
| CN116774280B (zh) | 一种鬼波压制定量质控方法、装置、电子设备及存储介质 | |
| CN104461500B (zh) | 一种信息处理方法及电子设备 | |
| CN113885080A (zh) | 沙漠区散射面波识别方法及装置 | |
| CN107015270B (zh) | 扫描信号计算方法及装置 | |
| Ali et al. | Performance evaluation of different inpainting algorithms for remotely sensed images |
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 |