CN113885075A - 一种叠后地震波相对波阻抗反演方法及系统 - Google Patents
一种叠后地震波相对波阻抗反演方法及系统 Download PDFInfo
- Publication number
- CN113885075A CN113885075A CN202111096121.7A CN202111096121A CN113885075A CN 113885075 A CN113885075 A CN 113885075A CN 202111096121 A CN202111096121 A CN 202111096121A CN 113885075 A CN113885075 A CN 113885075A
- Authority
- CN
- China
- Prior art keywords
- post
- stack seismic
- data
- wavelet packet
- wave
- 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 59
- 230000001360 synchronised effect Effects 0.000 claims abstract description 44
- 230000009466 transformation Effects 0.000 claims abstract description 42
- 238000001125 extrusion Methods 0.000 claims abstract description 38
- 230000010354 integration Effects 0.000 claims abstract description 10
- 238000009826 distribution Methods 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 230000008602 contraction Effects 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 5
- 230000006870 function Effects 0.000 claims description 4
- 230000002146 bilateral effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 7
- 230000007547 defect Effects 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/51—Migration
- G01V2210/514—Post-stack
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
技术领域
本发明涉及油气勘探地球物理处理方法技术领域，特别是涉及一种叠后地震波相对波阻抗反演方法及系统。
背景技术
本部分的陈述仅仅是提供了与本发明相关的背景技术信息，不必然构成在先技术。
叠后地震波阻抗反演是目前油气储层预测的一种重要技术，常见的叠后波阻抗反演技术如带限反演，主要是根据地层反褶积递归求出反射系数，或者是利用地震波形近似为反射系数，或者是通过建立波阻抗与反射系数的递推关系式逐层求解发射系数。该类方法的缺点是要求地震资料具有比较高的信噪比，有较宽的频带，同时保持相对振幅。
另一类常见的叠后波阻抗反演技术如测井约束下的宽带反演，主要通过各种反褶积计算得到的反射系数进行合成记录制作，通过目标函数对由井建立的波阻抗模型进行迭代修正来获得波阻抗。该类方法效果依赖于地震资料本身的品质，同样要求地震资料具有高信噪比，宽频带，保持相对振幅。
以随机反演为代表的地震约束下的测井内插外推技术，是从随机建模产生的一系列储层模型中优选出与地震数据最匹配的模型，该类方法突破了对地震宽频带的限制，但是反演效果依赖于各种统计特征，且统计特征要服从正态分布、对数正态分布或者能通过转换形成上述分布，要求工区内井点较多并且分布均匀，影响因素多，不利于广泛应用。
测井-地震联合反演技术，如特征反演，利用已知的储层参数信息建立三维属性模型，通过模型估算技术修正初始模型以匹配地震数据，得到反演结果。这类方法的主要缺点是要求井资料等已知信息较多，统计测井特征和地震特征的正确性难度较大。
针对常规叠后波阻抗反演技术如反褶积、波动方程、随机过程、特征分析以及动力学特征等受限于地震波频带宽度，很难获得高分辨率的波阻抗信息，或对井信息要求较多，在一些区域无法实现。
发明内容
为了解决上述问题，本发明提出了一种叠后地震波相对波阻抗反演方法及系统，基于同步挤压小波包变换进行叠后地震波相对波阻抗的反演，针对目标层位地震数据进行高分辨率的同步挤压小波包变换，在时频域内提取预设阈值范围内的模极大值曲线，通过选取适当的收缩母子波参数，利用同步挤压小波包反变换重构地震数据，实现重构高精度地震数据，最后结合带通积分获得相对声波阻抗。
为了实现上述目的，本发明采用如下技术方案：
第一方面，本发明提供一种叠后地震波相对波阻抗反演方法，包括：
获取目标层位的叠后地震波数据；
对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；
根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；
对重构的叠后地震波数据进行带通积分得到相对波阻抗。
作为可选择的实施方式，对叠后地震波数据逐道进行同步挤压小波包变换得到叠后地震波数据的时频能量分布，在时频能量分布中提取预设阈值范围内的模极大值曲线。
作为可选择的实施方式，对叠后地震波数据x(t)进行同步挤压小波包变换的过程为：
Tx(v,b)＝∫|Wx(p,b)|2δ(vx(p,b)-v)dp；
其中，Wx(p,b)是叠后地震波数据x(t)在点(p,b)处的小波包变换，vx(p,b)为叠后地震波数据x(t)在点(p,b)处的局域小波向量估计。
作为可选择的实施方式，叠后地震波数据x(t)的小波包变换Wx(p,b)为：
其中，
为母小波的共轭函数，ωpb(t)＝|p|sω(|p|s(x-b))e2πi(x-b)p，参数s∈(1/2,1)，|p|≥1；p,b∈R2；
叠后地震波数据x(t)的局域小波向量估计vx(p,b)为：
其中，Wx(p,b)≠0。
作为可选择的实施方式，在时频域内提取预设阈值范围内的模极大值曲线的过程为：对经同步挤压小波包变换得到的时频能量分布按列优先原则进行矩阵元素的扫描，得到尺度因子fj下，区间t∈(t0-δ,t0+δ)内，满足|W(fj,t)|≤|W(fj,t0)|的点(fj,t0)，则点(fj,t0)为同步挤压小波包变换的局部模极大值点，其中，δ为常数；在预设阈值范围内，将局部模极大值点连接后，得到模极大值曲线。
作为可选择的实施方式，重构的叠后地震波数据为：
其中，
为ωpb(t)的双边框架；Uk为第k个聚类子集。
作为可选择的实施方式，获取目标层位的叠后地震波数据的过程为：综合利用地质信息、测井以及合成地震记录，标定叠后地震波数据的目标层位，确定叠后地震波数据的分析时间范围。
第二方面，本发明提供一种叠后地震波相对波阻抗反演系统，包括：
数据获取模块，被配置为获取目标层位的叠后地震波数据；
同步挤压小波包变换模块，被配置为对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；
数据重构模块，被配置为根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；
反演模块，被配置为对重构的叠后地震波数据进行带通积分得到相对波阻抗。
第三方面，本发明提供一种电子设备，包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令，所述计算机指令被处理器运行时，完成第一方面所述的方法。
第四方面，本发明提供一种计算机可读存储介质，用于存储计算机指令，所述计算机指令被处理器执行时，完成第一方面所述的方法。
与现有技术相比，本发明的有益效果为：
本发明提供一种基于同步挤压小波包变换的叠后地震波相对波阻抗反演方法，核心是将叠后地震波数据进行同步挤压小波包变换，获得阈值范围内的模极大值曲线，结合同步挤压小波包反变换重构数据，再对重构数据进行带通积分获得相对波阻抗。无需井位信息，同时避免了传统波阻抗反演技术受限于地震波频带宽度及需要大量井位信息的缺陷，可适用于低信噪比地震数据。
本发明的叠后地震波相对波阻抗反演方法使用了具有极高时频分辨率的同步挤压小波包变换，更适合非线性非平稳地震信号的处理，抗噪声性能更强，保证计算结果的精度。
本发明的叠后地震波相对波阻抗反演方法适合地震信号大批量处理，适用于深度域地震数据处理。
本发明附加方面的优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1为本发明实施例1提供的叠后地震波相对波阻抗反演方法流程图；
图2为本发明实施例1提供的过井二维叠后偏移地震剖面；
图3为本发明实施例1提供的过井地震剖面经过同步挤压小波包变换处理后的重构剖面；
图4为本发明实施例1提供的利用叠后地震波相对波阻抗反演方法估计的波阻抗剖面；
图5为本发明实施例1提供的利用常规波阻抗技术估计的波阻抗剖面。
具体实施方式
下面结合附图与实施例对本发明做进一步说明。
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
实施例1
如图1所示，本实施例提供一种叠后地震波相对波阻抗反演方法，具体是一种基于同步挤压小波包变换的叠后地震波相对声波阻抗反演方法，包括：
S1：获取目标层位的叠后地震波数据；
S2：对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；
S3：根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；
S4：对重构的叠后地震波数据进行带通积分得到相对波阻抗。
所述步骤S1，标定叠后地震波数据的目标层位，具体包括：综合利用地质信息、测井以及合成地震记录等信息，准确标定叠后地震波数据的目标层位，同时确定叠后地震波数据的分析时间范围。
所述步骤S2中，在叠后地震波数据的分析时间范围内，利用同步挤压小波包变换得到目标层位的叠后地震波数据的时频能量分布，并在得到的时频能量分布中提取预设阈值范围内的模极大值曲线；具体包括如下步骤：
S2-1：对叠后地震波数据进行同步挤压小波包变换；同步挤压小波包变换是一种基于小波包变换和能量重排的自适应高分辨率时频分析方法；利用下式对叠后地震波数据x(t)进行同步挤压小波包变换：
Tx(v,b)＝∫|Wx(p,b)|2δ(vx(p,b)-v)dp， (1)
其中，Wx(p,b)是叠后地震波数据x(t)的小波包变换，即：
其中，
为母小波的共轭函数，ωpb(t)＝|p|sω(|p|s(x-b))e2πi(x-b)p，参数s∈(1/2,1)，|p|≥1；p,b∈R2；i为虚数单位。
vx(p,b)为叠后地震波数据x(t)在点(p,b)处的局域小波向量估计：
其中，Wx(p,b)≠0。
S2-2：从同步挤压小波包变换的时频能量分布中提取预设阈值范围内的模极大值曲线；
在时频能量分布中按列优先原则进行矩阵元素扫描，计算尺度因子fj下，区间t∈(t0-δ,t0+δ)内，满足|W(fj,t)|≤|W(fj,t0)|的点(fj,t0)，则点(fj,t0)为同步挤压小波包变换的局部模极大值点，其中，δ为常数；选取合适的阈值C，将能量占优的模极大值点连接后，得到模极大值曲线。
所述步骤S3中，针对目标层位的叠后地震波数据，选取适当的收缩母子波参数，利用同步挤压小波包反变换从模极大值曲线中计算重构数据；
可以理解的，适当的收缩母子波参数可以结合实际地震资料测试进行选择。
重构数据y(t)表示为：
其中，
为ωpb(t)的双边框架；Uk为第k个聚类子集。
所述步骤S4中，针对目标层位的叠后地震波数据，对重构数据利用带通积分器计算相对声波阻抗，完成对相对声波阻抗的反演。
本实施例基于同步挤压小波包变换算法进行叠后地震波相对声波阻抗反演，核心是通过时频域地震信号的同步挤压小波包分解，利用合理的模极大值曲线重构高分辨率地震数据，结合带通积分获得相对声波阻抗。
对上述方法进行验证，如图2所示是一条过井二维叠后偏移地震剖面，其中，黑色线为目标层上下层位，井轨迹如黑色曲线所示。如图3所示为过井地震剖面经过同步挤压小波包变换处理后的重构剖面，从图3中可以看到，重构剖面具有更高的分辨率，显示了目标层段更多细节信息。
图4为利用本实施例方法估计的波阻抗剖面(目标层段)，图4中椭圆示意范围内为测井解释含气储层所在区域。图5为利用常规地质统计波阻抗技术估计的波阻抗剖面(目标层段)，与图4相比，常规方法分辨率较低，椭圆示意区内，本实施例方法能更好的识别出储层的三个主要含气区域，并且对中间含气区域能显示出更多细节，细节清晰，常规方法分辨率较低。
实施例2
本实施例提供一种叠后地震波相对波阻抗反演系统，包括：
数据获取模块，被配置为获取目标层位的叠后地震波数据；
同步挤压小波包变换模块，被配置为对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；
数据重构模块，被配置为根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；
反演模块，被配置为对重构的叠后地震波数据进行带通积分得到相对波阻抗。
此处需要说明的是，上述模块对应于实施例1中所述的步骤，上述模块与对应的步骤所实现的示例和应用场景相同，但不限于上述实施例1所公开的内容。需要说明的是，上述模块作为系统的一部分可以在诸如一组计算机可执行指令的计算机系统中执行。
在更多实施例中，还提供：
一种电子设备，包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令，所述计算机指令被处理器运行时，完成实施例1中所述的方法。为了简洁，在此不再赘述。
应理解，本实施例中，处理器可以是中央处理单元CPU，处理器还可以是其他通用处理器、数字信号处理器DSP、专用集成电路ASIC，现成可编程门阵列FPGA或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
存储器可以包括只读存储器和随机存取存储器，并向处理器提供指令和数据、存储器的一部分还可以包括非易失性随机存储器。例如，存储器还可以存储设备类型的信息。
一种计算机可读存储介质，用于存储计算机指令，所述计算机指令被处理器执行时，完成实施例1中所述的方法。
实施例1中的方法可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器、闪存、只读存储器、可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。为避免重复，这里不再详细描述。
本领域普通技术人员可以意识到，结合本实施例描述的各示例的单元即算法步骤，能够以电子硬件或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.一种叠后地震波相对波阻抗反演方法，其特征在于，包括：
获取目标层位的叠后地震波数据；
对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；
根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；
对重构的叠后地震波数据进行带通积分得到相对波阻抗。
2.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，对叠后地震波数据逐道进行同步挤压小波包变换得到叠后地震波数据的时频能量分布，在时频能量分布中提取预设阈值范围内的模极大值曲线。
3.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，对叠后地震波数据x(t)进行同步挤压小波包变换的过程为：
Tx(v,b)＝∫|Wx(p,b)|2δ(vx(p,b)-v)dp；
其中，Wx(p,b)是叠后地震波数据x(t)在点(p,b)处的小波包变换，vx(p,b)为叠后地震波数据x(t)在点(p,b)处的局域小波向量估计。
4.如权利要求3所述的一种叠后地震波相对波阻抗反演方法，其特征在于，叠后地震波数据x(t)的小波包变换Wx(p,b)为：
其中，
为母小波的共轭函数，ωpb(t)＝|p|sω(|p|s(x-b))e2πi(x-b)p，参数s∈(1/2,1)，|p|≥1；p,b∈R2；
叠后地震波数据x(t)的局域小波向量估计vx(p,b)为：
其中，Wx(p,b)≠0。
5.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，在时频域内提取预设阈值范围内的模极大值曲线的过程为：对经同步挤压小波包变换得到的时频能量分布按列优先原则进行矩阵元素的扫描，得到尺度因子fj下，区间t∈(t0-δ,t0+δ)内，满足|W(fj,t)|≤|W(fj,t0)|的点(fj,t0)，则点(fj,t0)为同步挤压小波包变换的局部模极大值点，其中，δ为常数；在预设阈值范围内，将局部模极大值点连接后，得到模极大值曲线。
6.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，重构的叠后地震波数据为：
其中，
为ωpb(t)的双边框架；Uk为第k个聚类子集。
7.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，获取目标层位的叠后地震波数据的过程为：综合利用地质信息、测井以及合成地震记录，标定叠后地震波数据的目标层位，确定叠后地震波数据的分析时间范围。
8.一种叠后地震波相对波阻抗反演系统，其特征在于，包括：
数据获取模块，被配置为获取目标层位的叠后地震波数据；
同步挤压小波包变换模块，被配置为对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；
数据重构模块，被配置为根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；
反演模块，被配置为对重构的叠后地震波数据进行带通积分得到相对波阻抗。
9.一种电子设备，其特征在于，包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令，所述计算机指令被处理器运行时，完成权利要求1-7任一项所述的方法。
10.一种计算机可读存储介质，其特征在于，用于存储计算机指令，所述计算机指令被处理器执行时，完成权利要求1-7任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111096121.7A CN113885075B (zh) | 2021-09-15 | 2021-09-15 | 一种叠后地震波相对波阻抗反演方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111096121.7A CN113885075B (zh) | 2021-09-15 | 2021-09-15 | 一种叠后地震波相对波阻抗反演方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113885075A true CN113885075A (zh) | 2022-01-04 |
| CN113885075B CN113885075B (zh) | 2022-06-28 |
Family
ID=79009897
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111096121.7A Active CN113885075B (zh) | 2021-09-15 | 2021-09-15 | 一种叠后地震波相对波阻抗反演方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113885075B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110471104A (zh) * | 2019-08-26 | 2019-11-19 | 电子科技大学 | 基于智能特征学习的叠后地震反射模式识别方法 |
| CN110542923A (zh) * | 2019-09-02 | 2019-12-06 | 成都理工大学 | 一种快速高精度叠后地震阻抗反演方法 |
| WO2020160351A1 (en) * | 2019-02-01 | 2020-08-06 | University Of Georgia Research Foundation, Inc. | Contactless monitoring of sleep activities and body vital signs via seismic sensing |
| CN112083495A (zh) * | 2020-10-15 | 2020-12-15 | 中国石油化工股份有限公司 | 基于变分模态分解的同步压缩小波变换提高分辨率方法 |
-
2021
- 2021-09-15 CN CN202111096121.7A patent/CN113885075B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020160351A1 (en) * | 2019-02-01 | 2020-08-06 | University Of Georgia Research Foundation, Inc. | Contactless monitoring of sleep activities and body vital signs via seismic sensing |
| CN110471104A (zh) * | 2019-08-26 | 2019-11-19 | 电子科技大学 | 基于智能特征学习的叠后地震反射模式识别方法 |
| CN110542923A (zh) * | 2019-09-02 | 2019-12-06 | 成都理工大学 | 一种快速高精度叠后地震阻抗反演方法 |
| CN112083495A (zh) * | 2020-10-15 | 2020-12-15 | 中国石油化工股份有限公司 | 基于变分模态分解的同步压缩小波变换提高分辨率方法 |
Non-Patent Citations (4)
| Title |
|---|
| 唐贵基 等: "基于SSWPT的转子油膜失稳故障分析", 《动力工程学报》 * |
| 张雁: "同步压缩小波变换及其在地震资料处理中的应用研究", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 * |
| 曹军 等: "分频道积分技术在苏59气田储层预测中的应用", 《工程地球物理学报》 * |
| 裴正林: "小波理论及其在地震数据处理中的应用", 《地球物理学进展》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113885075B (zh) | 2022-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106597532B (zh) | 一种结合井资料与层位资料的叠前地震数据频带拓展方法 | |
| CN104360393A (zh) | 一种地震数据重建方法 | |
| CN109001800B (zh) | 一种基于地震数据的时频分解与气藏检测方法及系统 | |
| CN108845357B (zh) | 一种基于同步挤压小波变换估计地层等效品质因子的方法 | |
| Chen et al. | Iterative deblending using shaping regularization with a combined PNMO-MF-FK coherency filter | |
| CN104217448A (zh) | 基于迭代特征修正的磁共振快速成像方法及系统 | |
| Chen et al. | Seismic resolution enhancement by frequency-dependent wavelet scaling | |
| CN107179550A (zh) | 一种数据驱动的地震信号零相位反褶积方法 | |
| CN110687597B (zh) | 一种基于联合字典的波阻抗反演方法 | |
| Zhou et al. | Simultaneous deblending and interpolation using structure-oriented filters | |
| CN113885075B (zh) | 一种叠后地震波相对波阻抗反演方法及系统 | |
| CN115205136A (zh) | 一种基于傅里叶先验的图像去雨方法 | |
| CN106842297A (zh) | 井约束非稳态相位校正方法 | |
| Pan et al. | SVD-constrained MWNI with shaping theory | |
| CN112444867A (zh) | 一种基于压缩感知的地震资料频带拓展处理方法 | |
| CN110568489A (zh) | 一种块状介质的宽频反演方法 | |
| CN102353991A (zh) | 基于匹配地震子波的物理小波的地震瞬时频率分析方法 | |
| CN114114406B (zh) | 一种储层渗透率估计方法及装置 | |
| CN114035225B (zh) | 一种地震资料高分辨率处理bpfe方法及装置 | |
| CN113009579B (zh) | 地震数据反演方法及装置 | |
| CN111427091B (zh) | 挤压短时傅里叶变换的地震勘探信号随机噪声压制方法 | |
| CN104865603A (zh) | 一种针对倾斜层的svd滤波方法及装置 | |
| CN111665542B (zh) | 地震数据拓频方法及系统 | |
| CN114624765B (zh) | 一种相位域地震数据处理与重构方法、装置及可存储介质 | |
| CN117434592B (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 |