CN112213780B - 一种计及弹性阻抗二阶梯度的储层参数反演方法 - Google Patents
一种计及弹性阻抗二阶梯度的储层参数反演方法 Download PDFInfo
- Publication number
- CN112213780B CN112213780B CN202010946137.1A CN202010946137A CN112213780B CN 112213780 B CN112213780 B CN 112213780B CN 202010946137 A CN202010946137 A CN 202010946137A CN 112213780 B CN112213780 B CN 112213780B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- elastic impedance
- modulus
- rock
- porosity
- 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 44
- 239000011435 rock Substances 0.000 claims abstract description 66
- 239000012530 fluid Substances 0.000 claims abstract description 49
- 230000014509 gene expression Effects 0.000 claims abstract description 29
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 24
- 239000011148 porous material Substances 0.000 claims abstract description 9
- 239000011707 mineral Substances 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 11
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 11
- 238000004088 simulation Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000002310 reflectometry Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 101100419874 Caenorhabditis elegans snr-2 gene Proteins 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002637 fluid replacement therapy Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 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
- 230000000704 physical effect Effects 0.000 description 1
- 239000000126 substance Substances 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/282—Application of seismic models, synthetic seismograms
-
- 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
-
- 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
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
-
- 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
- G01V2210/6244—Porosity
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
技术领域
本发明涉及储层参数求解领域,尤其是涉及一种计及弹性阻抗二阶梯度的储层参数反演方法。
背景技术
流体类型识别和储层孔隙度预测是油藏描述的重要任务。流体替换模型将岩石的弹性模量与岩石的矿物类型、孔隙度、流体体积模量连接起来,为通过弹性参数预测进而计算孔隙度和流体体积模量奠定基础(Gassmann,1951)。
为了从地震数据中估测油气藏的性质,地球物理学家提出了不同类型的反射系数表达式。Aki和Richards(1980)提出了纵横波速度反射率表征的纵波反射系数近似式;同时,基于速度反射率表征的反射系数近似公式,针对不同的参数估测要求,以纵横波阻抗、模量、拉美参数、流体因子等表征的反射系数也相继被提出。类比于声波阻抗,Connolly(1999)提出了随入射角度变化的弹性阻抗。基于弹性阻抗概念,地球物理学家们开展了两步反演方法的研究:1)利用叠前地震数据预测得到弹性阻抗;2)利用估测的弹性阻抗提取储层弹性和物性参数(Zong et al.,2013;Chen et al.,2018)。
在储层弹性和物性参数求解过程中,贝叶斯框架下的叠前反演被广发应用。输入不同角度的叠前地震振幅数据,通过概率估计和约束,实现储层弹性参数和悟性参数的反演。然而反演过程中往往忽略了地震数据的二阶梯度影响,在一定程度上降低了反演的可靠性。
发明内容
本发明的目的就是为了克服上述现有技术存在忽略了地震数据的二阶梯度影响的缺陷而提供一种计及弹性阻抗二阶梯度的储层参数反演方法。
本发明的目的可以通过以下技术方案来实现:
一种计及弹性阻抗二阶梯度的储层参数反演方法,包括根据叠前地震数据预测弹性阻抗,利用所述弹性阻抗反演预测储层参数,所述利用弹性阻抗反演预测储层参数具体为,根据所述弹性阻抗的一阶梯度和二阶梯度,反演预测储层参数。
进一步地,所述储层参数的预测表达式为:
m=m0+γΔm
Δm=-H-1g
式中,m为未知数变量,即为所述储层参数的反演预测值,m0为未知数变量的初始模型,γ为步长,Δm为未知数变量的迭代值,H为弹性阻抗的二阶梯度,g为弹性阻抗的一阶梯度。
进一步地,所述弹性阻抗的二阶梯度的计算表达式为:
H≈ggT
式中,T为转置符号。
进一步地,所述弹性阻抗由密度、修正孔隙度及流体体积模量表征,所述弹性阻抗的一阶梯度的计算表达式为:
式中,EI为弹性阻抗的预测值,ρ为密度,φ为岩石的孔隙度,Kf为流体的体积模量,Δd为反演得到的EI和模拟得到的EI之间差异变量。
进一步地,所述计及弹性阻抗二阶梯度的储层参数反演方法面向内幕型孔缝储层,所述弹性阻抗的预测表达式为:
式中,EI为弹性阻抗的预测值,ρ为密度,θ为入射角,gsat为饱和岩石的剪切模量和纵波模量比值,gdry为干岩石的剪切模量和纵波模量比值,φn为修正的孔隙度,a为修正的孔隙度与岩石孔隙度之间线性关系的梯度值,Kf为流体的体积模量。
进一步地,所述弹性阻抗的预测表达式根据所述储层的反射系数近似式建立,所述储层的反射系数近似式根据所述储层的饱和岩石体积模量计算式建立。
进一步地,所述饱和岩石体积模量计算式为:
式中,Ksat为饱和岩石的体积模量,Kdry为干岩石的体积模量,Kf为流体的体积模量,K0为岩石矿物的体积模量,φ为岩石的孔隙度。
进一步地,所述计及弹性阻抗二阶梯度的储层参数反演方法面向含气储层,所述饱和岩石体积模量计算式为:
Ksat=Kdry+φnKf
式中,Ksat为饱和岩石的体积模量,Kdry为干岩石的体积模量,Kf为流体的体积模量,φn为修正的孔隙度,φ为岩石的孔隙度。
进一步地,所述储层的反射系数近似式的表达式为:
pμ(θ)=-2gsat sin2θ
式中,RPP为反射系数近似计算结果,θ为入射角,gsat为饱和岩石的剪切模量和纵波模量比值,gdry为干岩石的剪切模量和纵波模量比值,M0为岩石矿物的纵波模量,μ0为岩石矿物的剪切模量,ρ为密度,φn为修正的孔隙度,a为修正的孔隙度与岩石孔隙度之间线性关系的梯度值,Kf为流体的体积模量。
进一步地,所述计及弹性阻抗二阶梯度的储层参数反演方法面向内幕型孔缝储层,所述储层的反射系数近似式的表达式为:
与现有技术相比,本发明具有以下优点:
(1)本发明考虑了地震数据的二阶梯度影响,根据弹性阻抗的一阶梯度和二阶梯度,反演预测储层参数,具有较高的稳定性和可靠性。
(2)本发明首先采用干岩石体积模量和孔隙度之间的关系,提出修正孔隙度概念,并利用修正孔隙度和流体体积模量表征饱和岩石的体积模量,可以直观表征流体和孔隙变化对岩石弹性性质的影响;
其次,推导出一项利用岩石矿物模量、密度、修正孔隙度及流体模量表征的纵波反射系数公式,面向内幕型储层,舍掉矿物类型变化对反射系数的影响,并建立密度、修正孔隙度及流体体积模量表征的弹性阻抗,为模拟流体和孔隙对地震波反射振幅和弹性阻抗影响奠定了理论基础;
最后,建立了一套输入不同角度的纵波叠前地震道集预测弹性阻抗,进而利用弹性阻抗二阶梯度预测岩石密度、修正孔隙度和流体体积模量的方法和技术流程,具有较高的稳定性和可靠性。
(3)本发明方法可以稳定地估测储层的密度、修正孔隙度和流体模量,具有较好的抗噪性,同时实际数据预测表明,该方法能够提供比较合理的孔隙度和流体模量结果,对油气储层的储集空间及流体类型预测具有十分重要的意义。
附图说明
图1为本发明计及弹性阻抗二阶梯度的储层参数反演方法的流程示意图;
图2为将该发明应用于含噪声合成数据(信噪比2)得到的反演值与模型真实值之间的对比图;
图3为将本发明应用于实际地震数据反演得到的密度、修正孔隙度和流体体积模量示意图,其中(a)为密度示意图,(b)为修正孔隙度示意图,(c)为流体体积模量示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例1
本实施例提供一种计及弹性阻抗二阶梯度的储层参数反演方法,包括根据叠前地震数据预测弹性阻抗,利用弹性阻抗反演预测储层参数,利用弹性阻抗反演预测储层参数具体为,首先推导出一种利用修正孔隙度和流体体积模量表征的饱和流体岩石体积模量的简化近似式;其次,利用岩石体积模量简化近似式,针对内幕型孔缝储层,推导出利用密度、修正孔隙度和流体体积模量表征的反射系数近似式,并建立弹性阻抗表达式;最后,不同角度的叠前地震道集预测弹性阻抗体,利用弹性阻抗的二阶梯度,反演预测储层的密度、修正孔隙度和流体体积模量。
下面分别对本实施例方法和实验仿真结果进行具体描述。
一、本实施例方法的具体过程
本实施例计及弹性阻抗二阶梯度的储层参数反演方法具体包括以下步骤:
S1:以修正孔隙度和流体体积模量表征的饱和岩石体积模量简化近似式;
具体运算为:
根据流体替换模型,饱和岩石的体积模量为:
其中Ksat是饱和岩石的体积模量,Kdry为干岩石的体积模量,Kf是流体的体积模量,K0是岩石矿物的体积模量,φ是岩石的孔隙度。利用干岩石体积模量和孔隙度的关系,对饱和流体岩石体积模量进行改写:
在Kf<<K0假设条件下,面向含气储层,对饱和岩石体积模量进行简化:
Ksat=Kdry+φnKf (3)
S2:含修正孔隙度和流体模量的反射系数和弹性阻抗推导,过程如下:
以推导的简化饱和岩石体积模量为基础,推导获得利用岩石矿物体积模量、修正孔隙度及流体体积模量反射率表征的反射系数近似公式:
其中,
pμ(θ)=-2gsat sin2θ
式中,RPP为反射系数近似计算结果,θ为入射角,gsat为饱和岩石的剪切模量和纵波模量比值,gdry为干岩石的剪切模量和纵波模量比值,M0为岩石矿物的纵波模量,μ0为岩石矿物的剪切模量,ρ为密度,φn为修正的孔隙度,a为修正的孔隙度与岩石孔隙度之间线性关系的梯度值,a=0.1141,Kf为流体的体积模量。
针对内幕型孔缝储层,舍去岩石矿物模量M0和μ0的影响,得到简化的反射系数近似式:
由简化的反射系数近似式,推导得到弹性阻抗公式:
S3:纵波和横波反射振幅频率分量反演预测最大衰减参数值;
为了估测储层密度,修正的孔隙度和流体体积模量,建立了一种利用叠前地震道集预测弹性阻抗EI,进而利用预测的EI的一阶和二阶梯度估测未知数变量m(包含密度、修正孔隙度和流体模量)的方法。其中,EI对于未知数变量m的一阶和二阶梯度表示为:
H≈ggT (8)
其中,g是EI对于m的一阶梯度,gT是g的转置,Δd是反演得到的EI和模拟得到的EI之间差异变量,H为近似Hessian矩阵.为了得到反演值,采用全牛顿步长方法进行计算:
m=m0+γΔm (9)
其中m0为未知参数的初始模型,可以由测井和岩石物理模型估测并平滑处理得到,γ为步长,Δm是未知数参数变量m的迭代,计算如下:
Δm=-H-1g (10)
二、实验仿真结果
为了验证本发明在预测储层的密度、修正孔隙度及流体体积模量的可靠性,将本发明分别应用于合成地震数据和实际地震数据,进行处理分析。图2为将本发明应用于含随机噪声的合成地震数据(信噪比为2)时,反演得到的未知参数的真实值和预测值的对比。可以看出,在地震数据含有适量的随机噪声时,该发明得到的密度、修正孔隙度和流体模量的反演结果和真实值之间仍能较好吻合。图3为将该发明应用实际地震数据得到的密度、修正孔隙度和流体体积模量的反演结果,图中虚线框为储层位置。可以看出,利用本发明估测得到的密度和流体模量结果呈现低值,储层修正孔隙度呈现高值,合理可靠。图2和图3的数据处理效果验证了本发明中建立的反演方法预测孔隙度和流体体积模量的稳定性与可靠性,具有较大的应用潜力。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思做出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (1)
1.一种计及弹性阻抗二阶梯度的储层参数反演方法,包括根据叠前地震数据预测弹性阻抗,利用所述弹性阻抗反演预测储层参数,其特征在于,所述利用弹性阻抗反演预测储层参数具体为,根据所述弹性阻抗的一阶梯度和二阶梯度,反演预测储层参数;
所述储层参数的预测表达式为:
m=m0+γΔm
Δm=-H-1g
式中,m为未知数变量,即为所述储层参数的反演预测值,m0为未知数变量的初始模型,γ为步长,Δm为未知数变量的迭代值,H为弹性阻抗的二阶梯度,g为弹性阻抗的一阶梯度;
所述弹性阻抗的二阶梯度的计算表达式为:
H≈ggT
式中,T为转置符号;
所述弹性阻抗由密度、修正孔隙度及流体体积模量表征,所述弹性阻抗的一阶梯度的计算表达式为:
式中,EI为弹性阻抗的预测值,ρ为密度,φ为岩石的孔隙度,Kf为流体的体积模量,Δd为反演得到的EI和模拟得到的EI之间差异变量;
所述计及弹性阻抗二阶梯度的储层参数反演方法面向内幕型孔缝储层,所述弹性阻抗的预测表达式为:
式中,EI为弹性阻抗的预测值,ρ为密度,θ为入射角,gsat为饱和岩石的剪切模量和纵波模量比值,gdry为干岩石的剪切模量和纵波模量比值,φn为修正的孔隙度,a为修正的孔隙度与岩石孔隙度之间线性关系的梯度值,Kf为流体的体积模量;
所述弹性阻抗的预测表达式根据所述储层的反射系数近似式建立,所述储层的反射系数近似式根据所述储层的饱和岩石体积模量计算式建立;
所述饱和岩石体积模量计算式为:
式中,Ksat为饱和岩石的体积模量,Kdry为干岩石的体积模量,Kf为流体的体积模量,K0为岩石矿物的体积模量,φ为岩石的孔隙度;
所述计及弹性阻抗二阶梯度的储层参数反演方法面向含气储层,所述饱和岩石体积模量计算式为:
式中,Ksat为饱和岩石的体积模量,Kdry为干岩石的体积模量,Kf为流体的体积模量,φn为修正的孔隙度,φ为岩石的孔隙度;
所述储层的反射系数近似式的表达式为:
pμ(θ)=-2gsatsin2θ
式中,RPP为反射系数近似计算结果,θ为入射角,gsat为饱和岩石的剪切模量和纵波模量比值,gdry为干岩石的剪切模量和纵波模量比值,M0为岩石矿物的纵波模量,μ0为岩石矿物的剪切模量,ρ为密度,φn为修正的孔隙度,a为修正的孔隙度与岩石孔隙度之间线性关系的梯度值,Kf为流体的体积模量;
所述计及弹性阻抗二阶梯度的储层参数反演方法面向内幕型孔缝储层,所述储层的反射系数近似式的表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010946137.1A CN112213780B (zh) | 2020-09-10 | 2020-09-10 | 一种计及弹性阻抗二阶梯度的储层参数反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010946137.1A CN112213780B (zh) | 2020-09-10 | 2020-09-10 | 一种计及弹性阻抗二阶梯度的储层参数反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112213780A CN112213780A (zh) | 2021-01-12 |
CN112213780B true CN112213780B (zh) | 2021-11-09 |
Family
ID=74050653
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010946137.1A Active CN112213780B (zh) | 2020-09-10 | 2020-09-10 | 一种计及弹性阻抗二阶梯度的储层参数反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112213780B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112764106B (zh) * | 2021-01-28 | 2023-09-29 | 中海石油深海开发有限公司 | 一种流体类型的确定方法、装置、设备及存储介质 |
CN112965103B (zh) * | 2021-02-09 | 2022-08-23 | 中国石油大学(华东) | 一种多重孔隙储层叠前地震概率化多道反演方法 |
CN113552630B (zh) * | 2021-08-13 | 2022-03-04 | 广州海洋地质调查局 | 基于弹性阻抗的未固结地层渗透率预测方法及处理终端 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887132A (zh) * | 2009-05-15 | 2010-11-17 | 中国石油天然气股份有限公司 | 一种井震结合定量预测砂岩储层流体饱和度的方法 |
CN102495424B (zh) * | 2011-10-20 | 2013-06-26 | 西北大学 | 一种利用佐普立兹阻抗进行储层预测的方法 |
CN102736103A (zh) * | 2011-12-23 | 2012-10-17 | 电子科技大学 | 基于角度梯度弹性阻抗的储层预测方法 |
CN104570065B (zh) * | 2013-10-09 | 2017-05-03 | 中国石油化工股份有限公司 | 一种利用地震波阻抗定量反演孔隙度的方法 |
CN104597491B (zh) * | 2015-01-28 | 2017-07-21 | 中国石油大学(华东) | 一种基于方位弹性阻抗差异的裂缝储层参数反演方法 |
CN106405639B (zh) * | 2015-07-30 | 2018-11-23 | 中国石油化工股份有限公司 | 一种叠前地震储层岩性参数的反演方法 |
US10739480B2 (en) * | 2017-03-24 | 2020-08-11 | Exxonmobil Upstream Research Company | Full wavefield inversion with reflected seismic data starting from a poor velocity model |
CN110954948B (zh) * | 2018-09-27 | 2021-11-05 | 中国石油化工股份有限公司 | 一种岩石物理约束储层物性参数反演方法及系统 |
-
2020
- 2020-09-10 CN CN202010946137.1A patent/CN112213780B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112213780A (zh) | 2021-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112213780B (zh) | 一种计及弹性阻抗二阶梯度的储层参数反演方法 | |
Puryear et al. | Layer-thickness determination and stratigraphic interpretation using spectral inversion: Theory and application | |
RU2615591C1 (ru) | Многопараметрическая инверсия через зависящую от сдвига упругую полноволновую инверсию (fwi) | |
Chai et al. | Sparse reflectivity inversion for nonstationary seismic data | |
Carcione et al. | Computational poroelasticity—A review | |
Carcione | Theory and modeling of constant-Q P-and S-waves using fractional time derivatives | |
CN108535775B (zh) | 非平稳地震资料声波阻抗反演方法及装置 | |
Toverud et al. | Comparison of seismic attenuation models using zero-offset vertical seismic profiling (VSP) data | |
EP3259620B1 (en) | Multistage full wavefield inversion process that generates a multiple free data set | |
CN112882100B (zh) | 一种储层参数确定方法、装置、电子设备和存储介质 | |
CN112230283A (zh) | 一种基于测井曲线支持向量机建模的地震孔隙度预测方法 | |
CN108508489B (zh) | 一种基于波形微变化匹配的地震反演方法 | |
CN110618450B (zh) | 基于岩石物理建模的致密储层智能化含气性预测方法 | |
CN104483704A (zh) | 基于avo异常类型约束的剩余相位校正方法 | |
CN114488302A (zh) | 一种原位各向异性地应力场预测方法及系统 | |
US9977143B2 (en) | Method for attenuating multiples | |
CN104237936A (zh) | 一种油气检测的频变反演方法 | |
CN111427086B (zh) | 利用纵横波振幅频率分量预测地层最大衰减参数的方法 | |
Wu et al. | Iterative deblending based on the modified singular spectrum analysis | |
CN114779332A (zh) | 地震数据沉积背景去除方法、装置及电子设备 | |
CN102967884B (zh) | 波阻抗反演数据可靠性评价方法及装置 | |
CN115373022A (zh) | 一种基于振幅相位校正的弹性波场Helmholtz分解方法 | |
Amundsen et al. | Data-driven inversion/depth imaging derived from approximations to one-dimensional inverse acoustic scattering | |
Zand et al. | Integrated algorithm for high‐resolution crustal‐scale imaging using complementary OBS and streamer data | |
CN110333534B (zh) | 一种基于Biot理论的贝叶斯时移AVO反演方法及系统 |
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 |