CN112630832A - 基于衰减因子随入射角变化的含气性预测方法和装置 - Google Patents
基于衰减因子随入射角变化的含气性预测方法和装置 Download PDFInfo
- Publication number
- CN112630832A CN112630832A CN201910954988.8A CN201910954988A CN112630832A CN 112630832 A CN112630832 A CN 112630832A CN 201910954988 A CN201910954988 A CN 201910954988A CN 112630832 A CN112630832 A CN 112630832A
- Authority
- CN
- China
- Prior art keywords
- well
- stack
- angle
- attenuation factor
- longitudinal 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 34
- 230000008859 change Effects 0.000 claims abstract description 16
- 238000000605 extraction Methods 0.000 claims description 5
- 230000008707 rearrangement Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000003068 static effect Effects 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
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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
技术领域
本发明属于石油勘探领域,更具体地,涉及一种基于衰减因子随入射角变化的含气性预测方法和一种基于衰减因子随入射角变化的含气性预测装置。
背景技术
因为叠前地震数据比叠后地震数据包含了更丰富的信息,所以目前含气性检测研究多数是基于叠前地震数据。目前,在利用叠前地震数据进行含气性检测研究时主要有以下两种思路:
一、叠前反演,其主要存在的问题是低频模型的限制以及反演方法本身的多解性;
二、AVO属性,其主要存在问题是常规的AVO属性对含气性敏感性较低。
现有的上述两种含气性预测方案均难以令人满意。
发明内容
有鉴于此,本申请提出了一种在不进行地震反演的前提下,尽量多的挖掘叠前道集中的信息来进行含气性检测的技术方案。
根据本发明的一个方面,提供了一种基于衰减因子随入射角变化的含气性预测方法,所述方法包括:收集叠前角度域道集,并确定叠前角度域道集的入射角范围;按照预定的入射角范围对叠前角度域道集进行叠加,得到多个叠加后的地震数据;提取每个叠加后的地震数据的纵波衰减因子;按照入射角重新排列所提取的纵波衰减因子,得到叠前纵波衰减因子;提取目的层处叠前纵波衰减因子随入射角变化特征以预测井的含气性。
在一种可能的实施方式中,所述方法还包括:在对叠前角度域道集进行叠加前,对叠前角度域道集进行剩余时差校正,以对叠前角度域道集按照强轴拉平。
在一种可能的实施方式中,得到的叠加后的地震数据不超过5个。
在一种可能的实施方式中,所述提取目的层处叠前纵波衰减因子随入射角变化特征以预测井的含气性,包括:得到目的层处的叠前纵波衰减因子随入射角变化的曲线;根据所述曲线的梯度预测井的含气性。
在一种可能的实施方式中,所述根据所述曲线的梯度预测井的含气性,包括:如果工区内第一井的梯度绝对值大于第二井的梯度绝对值,则判断第一井的含气性高于第二井。
在一种可能的实施方式中,所述根据所述曲线的梯度预测井的含气性,包括:判断梯度绝对值大于预设阈值的井为含气井。
根据本发明的一个方面,提供了一种基于衰减因子随入射角变化的含气性预测装置,所述装置包括:叠前道集收集单元,用于收集叠前角度域道集,并确定叠前角度域道集的入射角范围;数据叠加单元,用于按照预定的入射角范围对叠前角度域道集进行叠加,得到多个叠加后的地震数据;衰减因子提取单元,用于提取每个叠加后的地震数据的纵波衰减因子;叠前重排单元,用于按照入射角重新排列所提取的纵波衰减因子,得到叠前纵波衰减因子;预测单元,用于提取目的层处叠前纵波衰减因子随入射角变化特征以预测井的含气性。
在一种可能的实施方式中,得到的叠加后的地震数据不超过5个。
在一种可能的实施方式中,所述预测单元具体包括:变化曲线获取单元,用于得到目的层处的叠前纵波衰减因子随入射角变化的曲线;含气性预测单元,用于根据所述曲线的梯度预测井的含气性。
在一种可能的实施方式中,所述含气性预测单元包括下列中的至少一者:第一判断模块,用于如果工区内第一井的梯度绝对值大于第二井的梯度绝对值,则判断第一井的含气性高于第二井;第二判断模块,用于判断梯度绝对值大于预设阈值的井为含气井。
根据本申请,通过叠前角度域道集的变换,获得具有入射角信息的叠后纵波衰减因子,并提取叠前纵波衰减因子随入射角变化特征进行含气性检测。由于本申请的地震属性是基于叠后地震数据进行计算的,稳定性高,并且不依赖于测井数据,所以确定性强。根据本申请的技术方案在含气性地区具有极大的推广应用价值。
附图说明
通过结合附图对本申请示例性实施方式进行更详细的描述,本申请的上述以及其它目的、特征和优势将变得更加明显,其中,在本申请示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出根据本申请的一个实施例的基于衰减因子随入射角变化的含气性预测方法的流程图。
图2示出根据本申请的一个实施例的基于衰减因子随入射角变化的含气性预测装置的结构框图。
图3示出根据本申请的一个示例性实施例得到的经剩余时差校正后的叠前角度域道集。
图4示出根据本申请的一个示例性实施例得到的叠前纵波衰减因子数据体。
图5示出根据本申请的一个示例性实施例得到的某高产气井的叠前纵波衰减因子随入射角变化的曲线图。
图6示出根据本申请的一个示例性实施例得到的某产水井的叠前纵波衰减因子随入射角变化的曲线图。
具体实施方式
下面将参照附图更详细地描述本申请的优选实施方式。虽然附图中显示了本申请的优选实施方式,然而应该理解,可以以各种形式实现本申请而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本申请更加透彻和完整,并且能够将本申请的范围完整地传达给本领域的技术人员。
请参见图1。图1示出根据本申请的一个实施例的基于衰减因子随入射角变化的含气性预测方法的流程图。如图所示,所述方法包括下列步骤。
步骤102,收集叠前角度域道集,并确定叠前角度域道集的入射角范围。
可收集叠前角度域道集,可确定有效的入射角范围。
还可对入射角范围内的叠前角度域道集进行剩余时差校正,以对叠前角度域道集按照强轴拉平。
步骤104,按照预定的入射角范围对叠前角度域道集进行叠加,得到多个叠加后的地震数据。
在一种可能的实施方式中,可按照得到的叠加后的地震数据不超过5个的原则进行叠加。发明人经过深入研究后发现,叠加后的地震数据在5个以内,可保证地震数据的信噪比不会太低,便于后续对叠加后的地震数据提出纵波衰减因子属性。
步骤106,提取每个叠加后的地震数据的纵波衰减因子。
可采用本领域技术人员认为适用的任意技术手段来提取叠加后的地震数据的纵波衰减因子,例如,可以采用现有提取纵波衰减因子的软件来进行提取。
步骤108,按照入射角重新排列所提取的纵波衰减因子,得到叠前纵波衰减因子。
因为是按照入射角重新排列,所以重新排列后的数据体的线方向、道方向、纵向采样深度域采样点间隔与叠前角度域道集完全一致,只是数据体中之前用于表示地震数据振幅的数据此时用于表示纵波衰减因子。
步骤110,提取目的层处叠前纵波衰减因子随入射角变化特征以预测井的含气性。
在一种可能的实施方式中,步骤110具体包括:得到目的层处的叠前纵波衰减因子随入射角变化的曲线;根据所述曲线的梯度预测井的含气性。
发明人对大量工区进行了深入研究,认为目的层处的叠前纵波衰减因子随入射角变化的曲线的梯度与井的含气性高低有很强的关联性,因此可根据该曲线的梯度来预测井的含气性。在此处,发明人创造性地借鉴了AVO(Amplitude Variation with Offset振幅随偏移距变化)的思想,又避免了现有AVO预测含气性时敏感性低的缺陷。
在一种可能的实施方式中,所述根据所述曲线的梯度预测井的含气性,包括:
如果工区内第一井的梯度绝对值大于第二井的梯度绝对值,则判断第一井的含气性高于第二井。
发明人通过对大量工区进行研究后发现,在同一工区内,如果一口井的叠前纵波衰减因子随入射角变化的曲线的梯度的绝对值较大,则其含气性通常较高。
在一种可能的实施方式中,所述根据所述曲线的梯度预测井的含气性,包括:
判断梯度绝对值大于预设阈值的井为含气井。
这个预设阈值通常对应含气井/产水井的划分节点。例如,通常认为含气饱和度在40%以下的井可能缺乏开采价值,将其称为产水井,则这个预设阈值可大约相当于含气饱和度在40%的井的叠前纵波衰减因子随入射角变化的曲线的梯度。
需要注意地是,针对不同工区,这个预设阈值往往是不同的。每个工区的预设阈值,要根据该工区的地质信息和工程先验信息等来确定。
根据本申请,通过叠前角度域道集的变换,获得具有入射角信息的叠后纵波衰减因子,并提取叠前纵波衰减因子随入射角变化特征进行含气性检测。由于本申请的地震属性是基于叠后地震数据进行计算的,稳定性高,并且不依赖于测井数据,所以确定性强。根据本申请的技术方案在含气性地区具有极大的推广应用价值。
图2示出根据本申请的一个实施例的基于衰减因子随入射角变化的含气性预测装置的结构框图。如图所示,所述装置包括叠前道集收集单元202、数据叠加单元204、衰减因子提取单元206、叠前重排单元208和预测单元210。
叠前道集收集单元202用于收集叠前角度域道集,并确定叠前角度域道集的入射角范围。
数据叠加单元204用于按照预定的入射角范围对叠前角度域道集进行叠加,得到多个叠加后的地震数据。
衰减因子提取单元206用于提取每个叠加后的地震数据的纵波衰减因子。
叠前重排单元208用于按照入射角重新排列所提取的纵波衰减因子,得到叠前纵波衰减因子。
预测单元210用于提取目的层处叠前纵波衰减因子随入射角变化特征以预测井的含气性。
在一种可能的实施方式中,得到的叠加后的地震数据不超过5个。
在一种可能的实施方式中,所述预测单元210具体包括:变化曲线获取单元,用于得到目的层处的叠前纵波衰减因子随入射角变化的曲线;含气性预测单元,用于根据所述曲线的梯度预测井的含气性。
在一种可能的实施方式中,所述含气性预测单元包括下列中的至少一者:第一判断模块,用于如果工区内第一井的梯度绝对值大于第二井的梯度绝对值,则判断第一井的含气性高于第二井;第二判断模块,用于判断梯度绝对值大于预设阈值的井为含气井。
应用示例
在某工区中,收集的叠前角度域道集的有效入射角范围为有效入射角范围为0°-36°。对该入射角范围内的叠前角度域道集进行剩余时差校正,结果如图3所示。
将有效入射角范围划分为0°-12°、12°-24°、24°-36°,分别进行叠加,获得3个角度叠加后的地震数据。按照中心角对这3个叠加后的地震数据命名,表示为seis_6、seis_18、seis_30。
采用现有软件对seis_6、seis_18、seis_30这3个角度叠加后的地震数据提取纵波衰减因子,获得对应的纵波衰减因子数据体。将这3个叠加后的地震数据对应的纵波衰减因子命名为seis_attention_6、seis_attention_18、seis_attention_30,这些数据体代表的是不同入射角范围内的地震数据纵波衰减信息。
将获得的3个纵波衰减因子数据体,重新按照入射角信息排列,获得叠前纵波衰减因子数据体,如图4所示。该数据体的线方向、道方向、纵向采样深度域采样点间隔与角度叠加数据完全一致,只是数据体中之前用于表示地震数据振幅的数据此时用于表示纵波衰减因子。叠前纵波衰减因子数据体每一个CDP(Common Depth Point共深度道集)中对应3个入射角信息;
然后,提取工区内各口井目的层处的纵波衰减因子随入射角(偏移距)变化特征。以工区内的高产气井A井和产水井B井为例,图5示出了A井的目的层处纵波衰减因子随入射角变化属性,图6示出了B井的目的层处纵波衰减因子随入射角变化属性。从图中可以看到高产气井与产水井之间纵波衰减因子随入射角变化曲线的梯度存在明显差异,其中高产气井A井的曲线梯度为-0.5752,而产水井B井的曲线梯度为-0.3455。二者截距的差异没有梯度的差异大,高产气井A井的曲线截距为-0.1817,产水井B井的曲线截距为-0.2009。本示例中,根据目的层处叠前纵波衰减因子随入射角变化特征可以敏感、准确地预测井的含气性。
本申请可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质,其上载有用于使处理器实现本申请的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身,诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如,通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
这里参照根据本申请实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本申请的各个方面。应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机可读程序指令实现。
以上已经描述了本申请的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (10)
1.一种基于衰减因子随入射角变化的含气性预测方法,其特征在于,所述方法包括:
收集叠前角度域道集,并确定叠前角度域道集的入射角范围;
按照预定的入射角范围对叠前角度域道集进行叠加,得到多个叠加后的地震数据;
提取每个叠加后的地震数据的纵波衰减因子;
按照入射角重新排列所提取的纵波衰减因子,得到叠前纵波衰减因子;
提取目的层处叠前纵波衰减因子随入射角变化特征以预测井的含气性。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
在对叠前角度域道集进行叠加前,对叠前角度域道集进行剩余时差校正,以对叠前角度域道集按照强轴拉平。
3.根据权利要求1所述的方法,其特征在于,得到的叠加后的地震数据不超过5个。
4.根据权利要求1所述的方法,其特征在于,所述提取目的层处叠前纵波衰减因子随入射角变化特征以预测井的含气性,包括:
得到目的层处的叠前纵波衰减因子随入射角变化的曲线;
根据所述曲线的梯度预测井的含气性。
5.根据权利要求4所述的方法,其特征在于,所述根据所述曲线的梯度预测井的含气性,包括:
如果工区内第一井的梯度绝对值大于第二井的梯度绝对值,则判断第一井的含气性高于第二井。
6.根据权利要求4所述的方法,其特征在于,所述根据所述曲线的梯度预测井的含气性,包括:
判断梯度绝对值大于预设阈值的井为含气井。
7.一种基于衰减因子随入射角变化的含气性预测装置,其特征在于,所述装置包括:
叠前道集收集单元,用于收集叠前角度域道集,并确定叠前角度域道集的入射角范围;
数据叠加单元,用于按照预定的入射角范围对叠前角度域道集进行叠加,得到多个叠加后的地震数据;
衰减因子提取单元,用于提取每个叠加后的地震数据的纵波衰减因子;
叠前重排单元,用于按照入射角重新排列所提取的纵波衰减因子,得到叠前纵波衰减因子;
预测单元,用于提取目的层处叠前纵波衰减因子随入射角变化特征以预测井的含气性。
8.根据权利要求7所述的装置,其特征在于,得到的叠加后的地震数据不超过5个。
9.根据权利要求7所述的装置,其特征在于,所述预测单元具体包括:
变化曲线获取单元,用于得到目的层处的叠前纵波衰减因子随入射角变化的曲线;
含气性预测单元,用于根据所述曲线的梯度预测井的含气性。
10.根据权利要求9所述的装置,其特征在于,所述含气性预测单元包括下列中的至少一者:
第一判断模块,用于如果工区内第一井的梯度绝对值大于第二井的梯度绝对值,则判断第一井的含气性高于第二井;
第二判断模块,用于判断梯度绝对值大于预设阈值的井为含气井。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910954988.8A CN112630832B (zh) | 2019-10-09 | 2019-10-09 | 基于衰减因子随入射角变化的含气性预测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910954988.8A CN112630832B (zh) | 2019-10-09 | 2019-10-09 | 基于衰减因子随入射角变化的含气性预测方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112630832A true CN112630832A (zh) | 2021-04-09 |
CN112630832B CN112630832B (zh) | 2023-09-22 |
Family
ID=75283302
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910954988.8A Active CN112630832B (zh) | 2019-10-09 | 2019-10-09 | 基于衰减因子随入射角变化的含气性预测方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112630832B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024067835A1 (zh) * | 2022-09-30 | 2024-04-04 | 中国石油化工股份有限公司 | 基于优势入射角-频率双域衰减的含气性预测方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101359057A (zh) * | 2007-07-31 | 2009-02-04 | 中国石油天然气股份有限公司 | 一种利用中心频率随入射角变化衰减信息进行气藏检测的方法 |
CN110295892A (zh) * | 2019-06-21 | 2019-10-01 | 中国石油大学(北京) | 多极子阵列声波测井中横波衰减因子的确定方法及装置 |
-
2019
- 2019-10-09 CN CN201910954988.8A patent/CN112630832B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101359057A (zh) * | 2007-07-31 | 2009-02-04 | 中国石油天然气股份有限公司 | 一种利用中心频率随入射角变化衰减信息进行气藏检测的方法 |
CN110295892A (zh) * | 2019-06-21 | 2019-10-01 | 中国石油大学(北京) | 多极子阵列声波测井中横波衰减因子的确定方法及装置 |
Non-Patent Citations (4)
Title |
---|
HONGBING LI ET AL.: "Analysis of center-frequency versus incident angle variations to detect gas reservoir", 《SEG/SAN ANTONIO 2007 ANNUAL MEETING》 * |
HONGBING LI ET AL.: "Analysis of center-frequency versus incident angle variations to detect gas reservoir", 《SEG/SAN ANTONIO 2007 ANNUAL MEETING》, 31 December 2007 (2007-12-31), pages 1490 - 1493 * |
张克非等: "基于岩石物理模型的页岩气储层横波速度估算方法", 《科学技术与工程》 * |
张克非等: "基于岩石物理模型的页岩气储层横波速度估算方法", 《科学技术与工程》, vol. 19, no. 11, 30 April 2019 (2019-04-30), pages 10 - 15 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024067835A1 (zh) * | 2022-09-30 | 2024-04-04 | 中国石油化工股份有限公司 | 基于优势入射角-频率双域衰减的含气性预测方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112630832B (zh) | 2023-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
US11016214B2 (en) | Dolomite reservoir prediction method and system based on well and seismic combination, and storage medium | |
Ross et al. | Automatic picking of direct P, S seismic phases and fault zone head waves | |
WO2017024523A1 (zh) | 一种射线弹性参数的反演方法 | |
Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
CN106501848B (zh) | 一种隧道掘进过程中隐性断层超前物探方法 | |
CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
US11221425B1 (en) | Generating a model for seismic velocities in a subsurface region using inversion with lateral variations | |
CN111897011B (zh) | 储层孔隙特征确定方法、装置及设备 | |
FR3032532A1 (zh) | ||
CN105240006A (zh) | 一种适用于火山岩储层的油水层识别方法 | |
CN112946782B (zh) | 一种致密油气储渗体地震精细刻画方法 | |
CN114236612A (zh) | 一种针对陆内水域地震背景噪声的信号处理方法 | |
CN112630832B (zh) | 基于衰减因子随入射角变化的含气性预测方法和装置 | |
WO2014108843A1 (en) | Noise attenuation of multicomponent microseismic data | |
US8441890B2 (en) | Removing ground-roll from geophysical data using source-receiver interferometry | |
CN111077574A (zh) | 一种确定地层弹性参数的方法、装置及系统 | |
CN111812716A (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
CN110727027B (zh) | 多期次河道砂的精细刻画方法 | |
CN111077577B (zh) | 井地联合储层描述方法及装置 | |
WO2023123971A1 (zh) | 基于vsp的深度域地震剖面层位标定方法及装置 | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
CN116027454A (zh) | 一种砂岩厚度等值线图的成图方法及智能绘制装置 | |
CN101285892A (zh) | 宽频地震勘探方法 | |
CN104181610A (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 |