CN109001813A - 一种压制多次波的方法、装置及系统 - Google Patents
一种压制多次波的方法、装置及系统 Download PDFInfo
- Publication number
- CN109001813A CN109001813A CN201810775402.7A CN201810775402A CN109001813A CN 109001813 A CN109001813 A CN 109001813A CN 201810775402 A CN201810775402 A CN 201810775402A CN 109001813 A CN109001813 A CN 109001813A
- Authority
- CN
- China
- Prior art keywords
- speed
- stack
- reference lamina
- wave
- optimization
- 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 70
- 230000001629 suppression Effects 0.000 claims abstract description 64
- 238000012545 processing Methods 0.000 claims abstract description 62
- 238000001228 spectrum Methods 0.000 claims abstract description 53
- 229910052704 radon Inorganic materials 0.000 claims abstract description 32
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 32
- 238000005457 optimization Methods 0.000 claims description 54
- 238000007906 compression Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 16
- 230000006835 compression Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 18
- 239000003550 marker Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 16
- 238000012937 correction Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 12
- 230000008569 process Effects 0.000 description 11
- 238000004590 computer program Methods 0.000 description 7
- 230000009466 transformation Effects 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000005553 drilling Methods 0.000 description 4
- 238000011160 research 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
- 238000009826 distribution Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000000750 progressive effect Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 241000761557 Lamina Species 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 239000011425 bamboo Substances 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 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
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
技术领域
本发明涉及地震数据处理技术领域，特别地，涉及一种压制多次波的方法、装置及系统。
背景技术
多次波会对地震资料造成多方面的影响。比如多次波会干扰有效波反射信号并降低资料的信噪比；当多次波干扰与一次波重叠时，使得有效反射波振幅、频率及相位发生畸变，分辨率下降，影响地震属性提取、反演等地震解释研究的可靠性。
多次波压制一直是地震勘探领域中的研究难点，目前已经提出并形成了多种衰减多次波的方法。但是，在通常情况下，层间多次波速度与一次波速度差异较小、混叠在一起，在速度谱上很难找到明确的分界，难以保证压制多次波的应用效果。如果选取速度过小、时差大，压制多次波能力弱，多次波残留多，达不到好的多次波压制效果。如果压制处理过度的话，在压制多次波的同时，还会去除一部分一次波的能量，从而降低有效波的地震反射能量。因此，本技术领域亟需一种可以更加有效的压制多次波的方法。
发明内容
本说明书实施例的目的是提供一种压制多次波的方法、装置及系统，可以提高对地震数据进行多次波压制的效果。
本说明书提供一种压制多次波的方法、装置及系统是包括如下方式实现的：
一种压制多次波的方法，包括：
获取目标处理层段中标志层的时间层位；
将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度；
基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度；
根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
本说明书提供的所述方法的另一个实施例中，所述获取目标处理层段中标志层的时间层位，包括：
根据测井合成地震记录与地震道的相关性确定待多次波压制时窗；
根据所述待多次波压制时窗确定目标处理层段，获取目标处理层段中标志层的时间层位。
本说明书提供的所述方法的另一个实施例中，所述对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，包括：
基于所述标志层对所述第一优化叠加速度进行镶边加权处理，获得第二优化叠加速度；
根据所述第二优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波。
本说明书提供的所述方法的另一个实施例中，所述基于所述标志层对应的速度对所述第一优化叠加速度进行镶边加权处理，包括：
根据下述计算模型对所述第一优化叠加速度进行镶边加权处理：
其中，V2表示第二优化叠加速度，V1表示第一优化叠加速度，t表示时间，Δt表示镶边时窗，a表示速度加权系数，Ti表示第i个标志层的时间值，0＜i＜K。
本说明书提供的所述方法的另一个实施例中，所述获取目标处理层段中标志层的时间层位，包括：
根据地层的波阻抗、地震反射能量、横向连续性中的一个或者多个确定目标处理层段中的标志层；
对所述标志层进行层位追踪以及闭合，获得所述标志层的时间层位。
另一方面，本说明书实施例还提供一种压制多次波的装置，包括：
标志层获取模块，用于获取目标处理层段中标志层的时间层位；
标志层速度确定模块，用于将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度；
层控速度拾取模块，用于基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度；
多次波压制模块，用于根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
本说明书提供的所述装置的另一个实施例中，所述标志层获取模块包括：
压制时窗确定单元，用于根据测井合成地震记录与地震道的相关性确定待多次波压制时窗；
标志层获取单元，用于根据所述待多次波压制时窗确定目标处理层段，获取目标处理层段中标志层的时间层位。
本说明书提供的所述装置的另一个实施例中，所述多次波压制模块包括：
镶边加权处理单元，用于根据下述计算模型对所述第一优化叠加速度进行镶边加权处理：
其中，V2表示第二优化叠加速度，V1表示第一优化叠加速度，t表示时间，Δt表示镶边时窗，a表示速度加权系数，Ti表示第i个标志层的时间值，0＜i＜K。
另一方面，本说明书实施例还提供一种压制多次波的设备，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
获取目标处理层段中标志层的时间层位；
将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度；
基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度；
根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
另一方面，本说明书实施例还提供一种压制多次波的系统，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个实施例中所述方法的步骤
本说明书一个或多个实施例提供的一种压制多次波的方法、装置及系统，可以通过获取待多次波压制层段内的标志层，以标志层的速度作为标定速度点进行叠加速度的拾取。从而有效降低多次波形成的速度能量噪声对叠加速度拾取的影响，提高叠加速度拾取的准确性。然后，进一步利用拾取的叠加速度进行拉东变换压制多次波，在保持有效波能量的同时，提高多次波压制的效果。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本说明书提供的一种压制多次波的方法实施例的流程示意图；
图2为本说明书提供的一个实施例中地震叠加剖面和标志层时间层位示意图；
图3为本说明书提供的另一个实施例中钻井的井震标定示意图；
图4为本说明书提供的另一个实施例中某CMP点的叠加速度谱及其地震道集示意图；
图5为本说明书提供的另一个实施例中某CMP点压制多次波后的叠加速度谱及其地震道集示意图；
图6为本说明书提供的一种压制多次波的装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本说明书一个或多个实施例中的附图，对本说明书一个或多个实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是说明书一部分实施例，而不是全部的实施例。基于说明书一个或多个实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书实施例方案保护的范围。
Radon(拉东)变换的基本假设是多次波时差与一次波时差不同(通常多次波相比一次波具有更低的速度)，当变换域中多次波和一次波同相轴轴分离时，则能够对多次波和一次波进行分离，并保留一次波成分。再反变换回时间-空间域，达到多次波压制的目的。多次波速度和一次波速度差异越大，二者剩余时差越大，采用Radon 变换压制多次波的效果越好。但，通常层间多次波速度与一次波速度差异较小、混叠在一起，在速度谱上很难找到明确的分界，难以保证压制多次波的应用效果。
相应的，本说明书实施例提供了一种压制多次波的方法，可以通过首先确定目标处理层段的地震反射标志层；然后，将标志层的时间层位加载到叠加速度谱上，在标志层速度的约束下拾取叠加速度。从而可以提高一次波速度拾取的准确性，进而提高后续进行拉东变换压制多次波的效果，并可以有效保持有效波的能量。
图1是本说明书提供的所述一种压制多次波的方法实施例流程示意图。虽然本说明书提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。
具体的一个实施例如图1所示，本说明书提供的压制多次波的方法的一个实施例中，所述方法可以包括：
S2：获取目标处理层段中标志层的时间层位。
可以通过分析钻井的测井曲线和合成地震记录，选择地层特征与地质界面基本一致的反射界面确定为地震反射标志层。本说明书的一个或者多个实施例中，可以根据波阻抗、地震反射能量、横向连续性等地层特征中的一个或者多个确定标志层。如可以选择强波阻抗、较强的地震反射能量以及横向连续性好、可连续追踪的反射界面作为标志层，以保证层位时间构造的准确性和一次波速度拾取的有效性。
如在四川盆地某目标工区，可以根据钻井的测井曲线和合成地震记录，在目标处理层段寒武系和震旦系，具有下寒武系龙王庙底沧浪铺底筇竹寺底震旦系灯三段底Z2dn3和灯二底Z2dn2等几个强反射界面。这些强反射界面均是厚层低速泥岩和高速碳酸盐岩地层的分界面，在研究区横向稳定发育，在地震叠加剖面都具有较强的地震反射能量。因此，可以确定这5个地质强反射界面为标志层。
然后，可以对所述标志层进行构造解释，获得所述标志层的时间层位数据。如可以在地震叠加数据上，开展地震标志层追踪和闭合，得到相应时间层位Ti，其中 0<i<K，K为标志层个数。
如可以利用地震叠加数据，对前述已经确定的5个标志层进行层位构造解释。如图2中的地震剖面上，地震解释的时间层位T1、T2、T3、T4和T5，分别对应下寒武系龙王庙底沧浪铺底筇竹寺底震旦系灯三段底Z2dn3和灯二底Z2dn2 的地层界面。
可选的，还可以进一步对解释的时间层位Ti进行平滑处理。在地震层位追踪和解释过程中，往往会产生一些异常值和无效值，使构造时间突变，这些值容易会影响后继处理效果，可以通过对解释的时间层位进行平滑处理，以减小解释误差的影响，提高获得时间层位的准确性。
本说明书的一个实施例中，可以先根据测井合成地震记录与地震道的相关性确定待多次波压制时窗，根据所述待多次波压制时窗确定目标处理层段，可以进一步提高数据处理的效率。
可以利用测井数据计算获得测井合成地震记录。
如，可以首先利用下述公式(1)计算地层反射系数：
其中，r表示反射系数；Den表示测井密度数据，DT表示测井声波时差数据； i表示测井曲线采样点序号。
然后，可以根据公式(2)将反射系数和地震子波进行褶积得到测井合成地震记录；
其中，Sy表示测井合成地震记录，w表示地震子波，r表示反射系数，i表示采样点序号，N表示子波长度。
图3表示对上述四川盆地某待测地区A井计算获得的测井合成地震记录示意图，图3中第a、b列分别表示测井声波时差曲线、密度曲线，第c列表示通过公式(1) 计算获得的反射系数；第d列表示将反射系数和32Hz雷克子波褶积得到的测井合成地震记录。右侧纵向坐标表示地层深度(Depth)，左侧纵向坐标表示时间(Time)。
然后，可以进行测井和地震时深标定。
如可以对上述A井的测井合成地震记录和井旁实际地震道的波组特征进行对比，保持强反射同相轴的时间对齐进行时移，使测井合成地震记录和井旁地震道的相关系数最大，完成全井段的井震时深标定，建立相应深-时关系。
图3展示出了A井井震标定对比结果：在4500m以上的地层中，第d列测井合成地震和第e列实际地震道波组符合好，同相轴对应良好，表明井震标定正确。
具体实施时，如可以将叠加地震数据导入到地震解释或反演软件中(诸如OpenWorks软件中SynTool井震标定模块)，根据已知钻井资料如A井的声波时差和密度等)制作合成记录，将测井波阻抗层位和地震层位对应起来，完成井震标定，建立相应深-时关系。
之后，可以根据测井合成地震记录与地震道的相关性，确定多次波压制时窗。
对图3中A井和地震道进行井震标定对比，合成地震记录在1.9s以上的中浅层与实际地震道波组符合好，同相轴对应良好，二者互相关系数为在0.8以上。说明地震资料在该段时窗范围内地震反射质量好，以有效波为主，信噪比高。从1.9s开始，测井合成地震记录和实际地震道的波组特征差异明显。除了2.0s～2.1s之间二个波峰一致性较好之外，整体上寒武系和震旦系都难以有效对比，二者互相关系数仅为0.14，同相轴时间、能量和相位均难以有效对比(见图中大括号指示的范围)。
然后，还可以通过地震速度谱、地震反射时间形态、地震正演模拟等综合分析证实，研究区的寒武系和震旦系地震资料中存在多次波干扰，且具有能量强、分布广的特点。由于多次波的干扰，导致在寒武系和震旦系的有效反射波被掩盖，振幅畸变，速度谱上难以有效拾取速度。
最终综合井震标定结果和分析结果，确定寒武系和震旦系地震资料中存在强能量的干扰，为需要进行多次波压制处理的层段。
利用测井合成地震记录与地震道的相关性可以准确地确定待多次波压制时窗。通过首先确定待多次波压制时窗，只需对待多次波压制时窗内的地震数据进行多次波压制处理，可以进一步提高数据处理的效率。
S4：将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度。
可以获取目标处理层段的叠加速度谱。可以获取叠前地震道集数据，按一定间隔选取地震共中心点(CMP)道集，计算相应的CMP道集叠加速度谱。具体计算方法可以参考现有技术实施，这里不做累述。当然，一些实施方式中，也可以选取如地震共反射点道集，计算共反射点道集叠加速度谱，用于后续数据处理。
然后，可以将所述标志层的时间层位叠加到任一CMP道集的叠加速度谱上，获得所述标志层在相应CMP道集上的时间值，在相应CMP道集对应的叠加速度谱上拾取所述标志层的时间值对应的速度。从而确定标志层相对该CMP道集的速度。
对上述四川盆地某目标工区，该地区虽然构造形态相对简单，具体实施时，研究区的速度分析网格密度仍然设定为200mⅹ200m，以减小大网格插值造成横向速度不准的问题。图4为获取的某CMP点的叠加速度谱及其叠前道集，其中，左侧图表示叠加速度谱，右侧图表示叠前道集。图4中的相关系数表示叠加能量大小，相关系数越大，表示该速度点的叠加能量值越大。
将标志层的时间层位加载到叠加速度谱上，如可以将5个标志层时间层位T1、 T2、T3、T4和T5，加载到CGG处理软件系统中。图4中T1、T2、T3、T4和T5对应的时间值表示5个标志层在该CMP点上的时间值。
沿标志层的时间值拾取能量团中心点作为相应标志层的速度。图4中在5个时间点上都存在明显或较为明显的能量团，这些在标志层位置上的能量团对应低速泥岩和高速碳酸盐岩地层的分界面，理论上为有效波反射能量。分别在T1、T2、T3、T4和 T5层位的时间点位置上拾取速度，图4中叉号表示的速度拾取点即为相应标志层在该 CMP点上对应的速度。
参考上述方式，对目标处理层段中所有CMP点确定各标志层对应的速度数据。
S6：基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度。
可以将上述步骤中拾取的标志层对应的速度点作为标定点，在保持纵向速度趋势的情况下，从CMP道集的叠加速度谱上根据能量团拾取叠加速度。然后，可以根据标志层对应的速度以及拾取的叠加速度进行插值处理，获得目标处理层段对应的第一优化叠加速度。
由图4可知，该CMP点的速度谱上，时间层位2.0s以上的中浅层能量团聚焦清楚，速度纵向变化规律明显，速度点易于拾取，多解性少。然而在2.1s以下由于存在强能量的多次波干扰，CMP道集信噪比很低，速度谱上能量团聚焦差，有效波反射和干扰波速度差异小，难以有效分辨速度谱能量团是有效波还是干扰波，无法拾取有效反射波的叠加速度。
可以以T1、T2、T3、T4和T5层位对应的速度点为标定点，在保持纵向速度趋势的情况下，根据能量团拾取叠加速度，如图4中的黑色虚线所示。
图4中在这5个标志层之间的位置，还存在着的其他能量团，通常是噪声干扰波形成的能量团(图中低速区出现的多次波)。某些干扰波能量团强于一次有效波，如果没有标志层的约束，极可能会出现误判拾取错误的速度。通过将标志层对应的速度点作为标定点进行速度拾取，有效避免了错误拾取标志层之间的干扰波形成的能量团，从纵向上避免了拾取不合理速度点。
同时，时间层位包含了地层的三维空间横向变化信息，利用标志层的时间层位为约束进行各CMP点的叠加速度拾取，可以在保证纵向一致性的同时，还可以进一步保证叠加速度拾取同地层的横向变化的一致性，进一步约束速度拾取异常现象，提高速度拾取的准确性。从而，利用本说明书上述实施例提供的方案，可以从横向以及纵向上同时提高最终拾取的叠加速度的准确性，进而提高后续多次波压制处理的效果。
S8：根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
可以利用上述步骤拾取的所述第一优化叠加速度对叠前地震道集进行动校正，获得动校正后的叠前地震道集数据。然后，可以对动校正后的叠前地震道集数据进行拉东正变换，获得Radon域的一次波及多次波分布。可以将Radon域的多次波能量置零后，再将剩余的一次波数据进行拉东逆变换，从而获得压制多次波后的地震道集数据，实现对叠前地震道集数据的多次波压制处理，获得多次波压制后的地震数据。
一些实施方式中，如可以采用双曲线时差方程(公式(3))对叠前地震道集进行动校正，得到动校正后的叠前地震道集：
其中，Δt表示剩余时差；t0表示零炮检距反射波的双程旅行时间；x表示炮检距；V1表示第一优化叠加速度。
在地震叠前道集中，一次波通过动校正后，消除了炮检距对反射波旅行时的影响，地震反射波被拉平。而多次波往往小于一次波速度，随着炮检距增加，一次波和多次波时差也随之增大，经过上述双曲线时差方程动校正之后，多次波呈向下弯曲的曲线形态，其剩余时差可近似为公式(4)所示的抛物线方程。
τ＝t-qx2 (4)
其中，τ表示双程零炮检距旅行时间，t表示双程旅行时间，q表示曲率，x表示炮检距。
然后，可以对动校正后的叠前地震道集进行Radon正变换。如可以按公式(5) 对动校正后的叠前地震道集进行抛物线Radon正变换，得到Radon域(τ-p域)的一次波及多次波分布；
其中，m表示τ-p域的地震道，d表示t-x域的地震道，x表示炮检距，t表示双程旅行时间，q表示曲率，τ表示双程零炮检距旅行时间。
根据抛物线Radon变换方程，t-x域中的抛物线轨迹，变换到τ-p域中成为一个点。将叠前地震道集中一次波和多次波的同相轴看成抛物线轨迹。由于二者抛物线曲率不同，因此在τ-p域中，一次波和多次波对应不同的点，可以将二者进行分离。
在τ-p域中，将多次波区域的能量置零，将一次波区域按公式(6)进行Radon 反变换至时间-炮检距域(t-x域)，获得压制多次波之后的地震道集；
其中，d表示t-x域的地震道，m表示τ-p域的地震道，x表示炮检距，t表示双程旅行时间，q表示曲率，τ表示双程零炮检距旅行时间。
本说明书的另一个实施例中，还可以先对所述第一优化叠加速度进行镶边加权处理，获得第二优化叠加速度。然后，再利用第二优化叠加速度对叠前地震道集进行动校正，获得动校正后的叠前地震道集数据，对动校正后的叠前地震道集数据进行基于拉东变换的多次波压制处理。其中，所述镶边加权处理可以包括：保持各标志层对应的叠加速度不变，将其他第一优化叠加速度进行适应性增加一定数值，获得第二优化叠加速度。从而可以进一步提高最终用于拉东变换进行多次波压制的叠加速度的准确性，提高多次波压制处理效果。
本说明书的一个或者多个实施例中，可以根据下述公式(7)对所述第一优化叠加速度进行镶边加权处理，获得第二优化叠加速度V2：
其中，V2表示第二优化叠加速度，V1表示第一优化叠加速度，t表示时间，Δt表示镶边时窗，a表示速度加权系数，Ti表示第i个标志层的时间值，0＜i＜K。
其中，Δt、a的取值可以根据实际地质情况进行预先设定。a的取值越大，压制多次波能力越强，但a的取值过大，可能会过多压制有效波。具体实施时，如可以通过模拟分析压制多次波前后的地震数据来确定a的取值，以保证有效反射波能量的同时，提高多次波压制的效果。
利用本实施例的方案，通过标定标志层位置拾取的叠加速度不变的情况下，再对其他速度进行镶边加权处理，可以在保证有效反射波能量的同时，进一步提高多次波压制的效果。
对上述四川盆地某待测区域，可以利用第二优化叠加速度对叠前地震道集进行动校正处理，然后通过Radon变换对叠前道集2.2s以下的深层进行多次波压制处理。
图5为四川盆地的实施例Radon变换压制多次波后速度谱和道集。对比图4，通过多次波压制处理后，目的层速度谱能量团相对集中，相对低速的多次波能量团被压制。对比图4和图5多次波压制处理前、后的地震道集，道集中深层的大炮检距，图 5中的同相轴下拉现象得到大幅的改善，目的层井震一致性明显提高。以上对比说明，利用本说明书实施例的方案进行有效压制多次波的同时，还能够很好地保持一次有效反射能量。
利用本说明书上述实施例提供的方案，通过获取待多次波压制层段内的标志层，将所述标志层的时间层位叠加至叠加速度谱上，获取标志层对应的叠加速度。以标志层对应的叠加速度为标定点，在保持纵向速度趋势的条件下根据叠加速度谱拾取叠加速度。从而有效降低多次波形成的速度能量噪声对叠加速度拾取的影响，提高叠加速度拾取的准确性。然后，进一步基于拾取的叠加速度进行拉东变换来压制多次波，可以提高多次波压制的准确性，并同时有效保持有效波能量。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。具体的可以参照前述相关处理相关实施例的描述，在此不做一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书一个或多个实施例提供的一种压制多次波的方法，可以通过获取待多次波压制层段内的标志层，以标志层的速度作为标定速度点进行叠加速度的拾取。从而有效降低多次波形成的速度能量噪声对叠加速度拾取的影响，提高叠加速度拾取的准确性。然后，进一步利用拾取的叠加速度进行拉东变换压制多次波，在保持有效波能量的同时，提高多次波压制的效果。
基于上述所述的压制多次波的方法，本说明书一个或多个实施例还提供一种压制多次波的装置。所述的装置可以包括使用了本说明书实施例所述方法的系统、软件(应用)、模块、组件、服务器等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。具体的，图6 表示说明书提供的一种压制多次波的装置实施例的模块结构示意图，如图6所示，所述装置可以包括：
标志层获取模块102，可以用于获取目标处理层段中标志层的时间层位；
标志层速度确定模块104，可以用于将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度；
层控速度拾取模块106，可以用于基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度；
多次波压制模块108，可以用于根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
利用上述实施例的方案，可以提高多次波压制的准确性，并同时有效保持有效波能量。
本说明书的另一个实施例中，所述标志层获取模块102可以包括：压制时窗确定单元以及标志层获取单元，其中，
所述压制时窗确定单元。可以用于根据测井合成地震记录与地震道的相关性确定待多次波压制时窗；
所述标志层获取单元，可以用于根据所述待多次波压制时窗确定目标处理层段，获取目标处理层段中标志层的时间层位。
利用上述实施例的方案，可以提高待多次波压制处理的时间层段确定的准确性，同时，还可以进一步提高数据处理的效率。
本说明书的另一个实施例中，所述多次波压制模块108可以包括：镶边加权处理单元，其中，
所述镶边加权处理单元，可以根据下述计算模型对所述第一优化叠加速度进行镶边加权处理：
其中，V2表示第二优化叠加速度，V1表示第一优化叠加速度，t表示时间，Δt表示镶边时窗，a表示速度加权系数，Ti表示第i个标志层的时间值，0＜i＜K。
利用上述实施例的方案，在保持标志层速度不变的情况下，进一步提高用于多次波压制的叠加速度的准确性，进而在保持有效波能量的同时，提高最终多次波压制的效果。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本说明书一个或多个实施例提供的一种压制多次波的装置，可以通过获取待多次波压制层段内的标志层，以标志层的速度作为标定速度点进行叠加速度的拾取。从而有效降低多次波形成的速度能量噪声对叠加速度拾取的影响，提高叠加速度拾取的准确性。然后，进一步利用拾取的叠加速度进行拉东变换压制多次波，在保持有效波能量的同时，提高多次波压制的效果。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种压制多次波的设备，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
获取目标处理层段中标志层的时间层位；
将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度；
基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度；
根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如RAM、ROM等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘；利用光学方式存储信息的装置如，CD或DVD。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。
需要说明的，上述所述的处理设备根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的一种压制多次波的设备，可以通过获取待多次波压制层段内的标志层，以标志层的速度作为标定速度点进行叠加速度的拾取。从而有效降低多次波形成的速度能量噪声对叠加速度拾取的影响，提高叠加速度拾取的准确性。然后，进一步利用拾取的叠加速度进行拉东变换压制多次波，在保持有效波能量的同时，提高多次波压制的效果。
本说明书还提供一种压制多次波的系统，所述系统可以为单独的压制多次波的系统，也可以应用在地震数据处理中，或者应用在多种类型的地震数据采集系统中。所述的系统可以为单独的计算机，也可以包括使用了本说明书的一个或多个所述方法或一个或多个实施例装置的实际操作装置(如激发装置、接收装置、接收线路)等。所述压制多次波的系统可以包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个或者多个实施例中所述方法的步骤。
需要说明的，上述所述的系统根据方法或者装置实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的一种压制多次波的系统，可以通过获取待多次波压制层段内的标志层，以标志层的速度作为标定速度点进行叠加速度的拾取。从而有效降低多次波形成的速度能量噪声对叠加速度拾取的影响，提高叠加速度拾取的准确性。然后，进一步利用拾取的叠加速度进行拉东变换压制多次波，在保持有效波能量的同时，提高多次波压制的效果。
需要说明的是，本说明书上述所述的装置或者系统根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类、存储介质+程序实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
尽管本说明书实施例内容中提到如第一优化叠加速度的拾取、压制多次波等获取、定义、交互、计算、判断等操作和数据描述，但是，本说明书实施例并不局限于必须是符合标准数据模型/模板或本说明书实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书的可选实施方案范围之内。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
上述实施例阐明的系统、装置、模块或单元，可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、平板计算机或者这些设备中的任何设备的组合。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/ 或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法或者设备中还存在另外的相同要素。
本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书一个或多个实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本本说明书一个或多个实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述并不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书的实施例而已，并不用于限制本说明书。对于本领域技术人员来说，本说明书可以有各种更改和变化。凡在本说明书的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书的权利要求范围之内。
Claims (10)
1.一种压制多次波的方法，其特征在于，包括：
获取目标处理层段中标志层的时间层位；
将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度；
基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度；
根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
2.根据权利要求1所述的压制多次波的方法，其特征在于，所述获取目标处理层段中标志层的时间层位，包括：
根据测井合成地震记录与地震道的相关性确定待多次波压制时窗；
根据所述待多次波压制时窗确定目标处理层段，获取目标处理层段中标志层的时间层位。
3.根据权利要求1所述的压制多次波的方法，其特征在于，所述对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，包括：
基于所述标志层对所述第一优化叠加速度进行镶边加权处理，获得第二优化叠加速度；
根据所述第二优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波。
4.根据权利要求3所述的压制多次波的方法，其特征在于，所述基于所述标志层对所述第一优化叠加速度进行镶边加权处理，包括：
根据下述计算模型对所述第一优化叠加速度进行镶边加权处理：
其中，V2表示第二优化叠加速度，V1表示第一优化叠加速度，t表示时间，Δt表示镶边时窗，a表示速度加权系数，Ti表示第i个标志层的时间值，0＜i＜K。
5.根据权利要求1所述的压制多次波的方法，其特征在于，所述获取目标处理层段中标志层的时间层位，包括：
根据地层的波阻抗、地震反射能量、横向连续性中的一个或者多个确定目标处理层段中的标志层；
对所述标志层进行层位追踪以及闭合，获得所述标志层的时间层位。
6.一种压制多次波的装置，其特征在于，包括：
标志层获取模块，用于获取目标处理层段中标志层的时间层位；
标志层速度确定模块，用于将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度；
层控速度拾取模块，用于基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度；
多次波压制模块，用于根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
7.根据权利要求6所述的压制多次波的装置，其特征在于，所述标志层获取模块包括：
压制时窗确定单元，用于根据测井合成地震记录与地震道的相关性确定待多次波压制时窗；
标志层获取单元，用于根据所述待多次波压制时窗确定目标处理层段，获取目标处理层段中标志层的时间层位。
8.根据权利要求6所述的压制多次波的装置，其特征在于，所述多次波压制模块包括：
镶边加权处理单元，用于根据下述计算模型对所述第一优化叠加速度进行镶边加权处理：
其中，V2表示第二优化叠加速度，V1表示第一优化叠加速度，t表示时间，Δt表示镶边时窗，a表示速度加权系数，Ti表示第i个标志层的时间值，0＜i＜K。
9.一种压制多次波的设备，其特征在于，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
获取目标处理层段中标志层的时间层位；
将所述时间层位叠加到叠加速度谱上，在所述叠加速度谱上拾取所述标志层对应的速度；
基于所述标志层对应的速度在所述叠加速度谱上拾取叠加速度，获得第一优化叠加速度；
根据所述第一优化叠加速度对所述目标处理层段的叠前地震道集基于拉东变换压制多次波，获得多次波压制后的地震数据。
10.一种压制多次波的系统，其特征在于，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现权利要求1-5中任意一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810775402.7A CN109001813B (zh) | 2018-07-16 | 2018-07-16 | 一种压制多次波的方法、装置及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810775402.7A CN109001813B (zh) | 2018-07-16 | 2018-07-16 | 一种压制多次波的方法、装置及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109001813A true CN109001813A (zh) | 2018-12-14 |
| CN109001813B CN109001813B (zh) | 2020-01-07 |
Family
ID=64598364
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810775402.7A Active CN109001813B (zh) | 2018-07-16 | 2018-07-16 | 一种压制多次波的方法、装置及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109001813B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111239827A (zh) * | 2020-03-09 | 2020-06-05 | 吉林大学 | 基于局部相似系数的三维地震数据多次波压制方法 |
| CN111562616A (zh) * | 2019-02-14 | 2020-08-21 | 中国石油天然气股份有限公司 | 地震数据散射噪音压制方法及装置 |
| CN111624653A (zh) * | 2020-04-20 | 2020-09-04 | 中国石油天然气股份有限公司 | 地震速度自动拾取方法及装置 |
| CN111751870A (zh) * | 2019-03-26 | 2020-10-09 | 中国石油天然气集团有限公司 | 叠后层间多次波压制方法及装置 |
| CN112213774A (zh) * | 2020-06-22 | 2021-01-12 | 北京汤固能源科技有限公司 | 一种浅层q模型估计方法及装置 |
| CN114185094A (zh) * | 2020-09-14 | 2022-03-15 | 中国石油化工股份有限公司 | Rms-svd多次波压制方法、装置、电子设备及介质 |
| CN116819627A (zh) * | 2023-06-30 | 2023-09-29 | 中海石油(中国)有限公司深圳分公司 | 地震弱信号的增强方法、装置、设备和介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103777237A (zh) * | 2012-10-17 | 2014-05-07 | 中国石油化工股份有限公司 | 一种基于空变加权镶边波数域滤波的地表高程平滑方法 |
| CN104570110A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于纵横波匹配的多分量资料联合速度分析方法 |
| CN107656313A (zh) * | 2017-10-12 | 2018-02-02 | 中国石油天然气集团公司 | 一种拉东域多次波压制方法及装置 |
| CN107678062A (zh) * | 2017-09-15 | 2018-02-09 | 上海海洋大学 | 双曲Radon域综合预测反褶积和反馈循环方法压制多次波模型构建方法 |
-
2018
- 2018-07-16 CN CN201810775402.7A patent/CN109001813B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103777237A (zh) * | 2012-10-17 | 2014-05-07 | 中国石油化工股份有限公司 | 一种基于空变加权镶边波数域滤波的地表高程平滑方法 |
| CN104570110A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于纵横波匹配的多分量资料联合速度分析方法 |
| CN107678062A (zh) * | 2017-09-15 | 2018-02-09 | 上海海洋大学 | 双曲Radon域综合预测反褶积和反馈循环方法压制多次波模型构建方法 |
| CN107656313A (zh) * | 2017-10-12 | 2018-02-02 | 中国石油天然气集团公司 | 一种拉东域多次波压制方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| 褚玉环 等: "拉冬变换压制多次波技术在大庆探区的应用", 《大庆石油地质与开发》 * |
| 赵忠华: "地震层位约束下的速度分析技术", 《采油工程文集》 * |
| 黄晓容: "采用地震勘探法超前探测矿区隐伏导水构造", 《矿业安全与环保》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111562616A (zh) * | 2019-02-14 | 2020-08-21 | 中国石油天然气股份有限公司 | 地震数据散射噪音压制方法及装置 |
| CN111751870A (zh) * | 2019-03-26 | 2020-10-09 | 中国石油天然气集团有限公司 | 叠后层间多次波压制方法及装置 |
| CN111239827A (zh) * | 2020-03-09 | 2020-06-05 | 吉林大学 | 基于局部相似系数的三维地震数据多次波压制方法 |
| CN111239827B (zh) * | 2020-03-09 | 2021-07-30 | 吉林大学 | 基于局部相似系数的三维地震数据多次波压制方法 |
| CN111624653A (zh) * | 2020-04-20 | 2020-09-04 | 中国石油天然气股份有限公司 | 地震速度自动拾取方法及装置 |
| CN112213774A (zh) * | 2020-06-22 | 2021-01-12 | 北京汤固能源科技有限公司 | 一种浅层q模型估计方法及装置 |
| CN112213774B (zh) * | 2020-06-22 | 2024-06-04 | 北京汤固能源科技有限公司 | 一种浅层q模型估计方法及装置 |
| CN114185094A (zh) * | 2020-09-14 | 2022-03-15 | 中国石油化工股份有限公司 | Rms-svd多次波压制方法、装置、电子设备及介质 |
| CN116819627A (zh) * | 2023-06-30 | 2023-09-29 | 中海石油(中国)有限公司深圳分公司 | 地震弱信号的增强方法、装置、设备和介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109001813B (zh) | 2020-01-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109001813A (zh) | 一种压制多次波的方法、装置及系统 | |
| CN111239802B (zh) | 基于地震反射波形和速度谱的深度学习速度建模方法 | |
| US9759826B2 (en) | System and method for generating an implicit model of geological horizons | |
| CN110031896A (zh) | 基于多点地质统计学先验信息的地震随机反演方法及装置 | |
| CN105785440B (zh) | 一种矿井槽波双分量地震信号频散曲线提取方法 | |
| CN103758511B (zh) | 一种井下逆时偏移成像识别隐蔽储层的方法及装置 | |
| CN103424777B (zh) | 一种提高地震成像分辨率的方法 | |
| CN107831542A (zh) | Ddw高精度深度域井震匹配方法 | |
| CN109633745B (zh) | 一种三维构造图的制图方法及装置 | |
| CN102636809B (zh) | 一种传播角度域共成像点道集的生成方法 | |
| CN109870719A (zh) | 一种碳酸盐岩致密薄储层分布确定方法、装置及系统 | |
| CN109188520A (zh) | 薄储层厚度预测方法及装置 | |
| CN102937720A (zh) | 井控提高地震资料分辨率的方法 | |
| US20210190983A1 (en) | Full waveform inversion in the midpoint-offset domain | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN106338764B (zh) | 生物启发计算地层圈闭油气藏超剥线识别方法 | |
| CN117546051A (zh) | 使用s波速度模型和机器学习进行地震成像的方法和系统 | |
| CN101738636A (zh) | 一种三维vsp高斯束法多波联合偏移成像方法 | |
| CN105607119B (zh) | 近地表模型构建方法与静校正量求取方法 | |
| CN111983683A (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
| CN106199694A (zh) | 基于深变子波的合成记录制作方法 | |
| CN105093288B (zh) | 一种基于运动学波场属性的绕射波分离方法 | |
| CN108508481B (zh) | 一种纵波转换波地震数据时间匹配的方法、装置及系统 | |
| CN103675900A (zh) | 一种确定转换波叠前时间偏移最佳速度剖面的方法 | |
| CN107918153B (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 |