储层物性参数直接反演方法、装置、电子设备及介质制造方法及图纸

本申请公开了一种储层物性参数直接反演方法、装置、电子设备及介质。该方法包括：根据工区内的地震数据、储层物性测井数据，建立储层物性参数的初始背景模型；构建储层物性参数表征的地震响应特征方程；根据初始背景模型与地震响应特征方程，计算初始背景模型对应的正演地震数据；根据工区内的地震数据、声波测井数据，建立纵、横波速度和密度的低频模型；根据初始背景模型、低频模型与正演地震数据，开展储层物性反演，得到反演结果；将反演结果与初始背景模型融合，获得最终储层物性反演结果。本发明专利技术实现物性参数的直接预测，解决目前利用两步反演方法预测物性参数带来的不确定性强、精度低的问题，为储层描述和油气开发奠定坚实基础。基础。基础。

【技术实现步骤摘要】
储层物性参数直接反演方法、装置、电子设备及介质


[0001]本专利技术涉及油气地球物理
，更具体地，涉及一种储层物性参数直接反演方法、装置、电子设备及介质。

技术介绍

[0002]储层物性参数可以表征储层的基本性质，决定了油气藏的储量、品质、开发的难易程度等，在油气藏勘探开发中起到关键作用。因此，储层物性预测一直是石油勘探领域的热点和难点。
[0003]目前常用的储层物性参数预测方法有三类。第一类是地震多属性物性反演，其思路是建立物性和地震属性之间定量的关系，进而反演物性参数；第二类是地质统计学物性反演，主要是通过地质统计学随机模拟的思路实现物性参数的反演。它们的不足之处在于：(1)岩石物理的机理不清楚，(2)需要大量的测井数据。第三类是基于岩石物理的物性反演，包括：统计岩石物理建模方法、地质统计学物性参数反演方法。但是这类方法大都是先地震反演、后物性反演的两步反演，缺乏对于误差积累的考虑，每一步反演都不可避免地产生误差，两步法严重放大了最终结果中包含的误差和不确定性，影响了物性参数预测的精度和稳定性，进而对后续的储层研究产生影响。
[0004]要克服以上几类常见方法存在的缺陷，需要发展从地震数据直接预测物性参数的方法。目前，基于地震资料针对储层物性参数的预测研究工作还处于初步探索阶段。
[0005]因此，有必要开发一种基于岩石物理模型的储层物性参数直接反演方法、装置、电子设备及介质。
[0006]公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

技术实现思路

[0007]本专利技术提出了一种储层物性参数直接反演方法、装置、电子设备及介质，其针对实际工区的储层特点建立岩石物理模型，建立地震响应与储层物性参数的直接联系，利用模型约束的贝叶斯反演理论，实现物性参数的直接预测，从而解决目前利用两步反演方法预测物性参数带来的不确定性强、精度低的问题，为储层描述和油气开发奠定坚实基础。
[0008]第一方面，本公开实施例提供了一种储层物性参数直接反演方法，包括：
[0009]根据工区内的地震数据、储层物性测井数据，建立储层物性参数的初始背景模型；
[0010]构建储层物性参数表征的地震响应特征方程；
[0011]根据所述初始背景模型与所述地震响应特征方程，计算所述初始背景模型对应的正演地震数据；
[0012]根据工区内的地震数据、声波测井数据，建立纵、横波速度和密度的低频模型；
[0013]根据所述初始背景模型、所述低频模型与所述正演地震数据，开展储层物性反演，
得到反演结果；
[0014]将所述反演结果与所述初始背景模型融合，获得最终储层物性反演结果。
[0015]优选地，构建储层物性参数表征的地震响应特征方程包括：
[0016]根据泰勒级数理论和广义线性反演理论，构建储层物性参数表征的地震响应特征广义方程；
[0017]根据地震岩石物理建模理论，构建储层物性参数表征的地震响应特征方程。
[0018]优选地，储层物性参数表征的地震响应特征广义方程为：
[0019][0020]其中，d为地震数据，W为子波，R为地震反射系数，e是地震数据中包含的噪音，V
m
、S
w
分别为储层的孔隙度、矿物含量和含水饱和度，ΔV
m
、ΔS
w
表示孔隙度、矿物含量和含水饱和度的扰动量，V
m0
、S
w0
表示孔隙度、矿物含量和含水饱和度的背景值，为背景物性参数V
m0
、S
w0
对应的地震反射系数，为反射系数对弹性参数的偏导数，V1、V2、V3表示三个不同的弹性参数，为反射系数对储层物性参数的偏导数，为弹性参数对储层物性参数的偏导数。
[0021]优选地，储层物性参数表征的地震响应特征方程为：
[0022][0023]其中，为反射系数对弹性参数的偏导数，为反射系数对弹性参数的偏导数，为反射系数对弹性参数的偏导数，V
Pi
,V
Pi+1
,V
Si
,V
Si+1
,ρ
i
,ρ
i+1
分别是不同层反射界面两侧的纵波速度、横波速度和密度，γ是纵横波速度比，为纵波速度对孔隙度、白云岩含量和含水饱和度的偏导数，
为横波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为横波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为横波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为横波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为密度对孔隙度、白云岩含量和含水饱和度的偏导数，饱和度的偏导数，饱和度的偏导数，饱和度的偏导数，饱和度的偏导数，
K
sat
和μ
sat
分别为饱和岩石的体积模量和剪切模量；K
m
和μ
m
分别为岩石基质的体积模量和剪切模量，K
cal
和μ
cal
为方解石的体积模量和剪切模量；K
dol
和μ
dol
为白云石的体积模量和剪切模量；K
f
为孔隙流体的体积模量；M
sat
为饱和岩石的纵波模量，V
dol
是白云岩含量，ρ
m
、ρ
f
、ρ
dol
、ρ
cal
、ρ
hc
、ρ
w
分别为基质矿物、孔隙流体、白云石、方解石、油气和水的密度。
[0024]优选地，根据所述初始背景模型与所述地震响应特征方程，计算所述初始背景模型对应的正演地震数据包括：
[0025]将所述初始背景模型代入岩石物理模型和所述地震响应特征方程中，通过褶积计算初始背景模型对应的正演地震数据。
[0026]优选地，根据所述初始背景模型、所述低频模型与所述地震数据，开展储层物性反演，得到反演结果包括：
[0027]建立目标函数；
[0028]将叠前地震数据CRP道集、所述正演地震数据、所述初始背景模型与所述低频模型代入所述目标函数中，通过迭代重加权最小二乘方法，求解所述目标函数，获得所述反演结果。
[0029]优选地，目标函数为：
[0030][0031]其中，α是权重因子；C和ξ为正则化约束项，W为子波矩阵，d
i
为第i个角度道且具有k个元素的地震列向量，A
i
、B
i
、C
i
为系数a(θ)、b(θ)、c(θ)组成的k
×
k阶矩阵，k阶矩阵，k阶矩阵，
[0032]作为本公开实施例的一种具体实现方式，
[0033]第二方面，本公开实施例还提供了一种储本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种储层物性参数直接反演方法，其特征在于，包括：根据工区内的地震数据、储层物性测井数据，建立储层物性参数的初始背景模型；构建储层物性参数表征的地震响应特征方程；根据所述初始背景模型与所述地震响应特征方程，计算所述初始背景模型对应的正演地震数据；根据工区内的地震数据、声波测井数据，建立纵、横波速度和密度的低频模型；根据所述初始背景模型、所述低频模型与所述正演地震数据，开展储层物性反演，得到反演结果；将所述反演结果与所述初始背景模型融合，获得最终储层物性反演结果。2.根据权利要求1所述的储层物性参数直接反演方法，其中，构建储层物性参数表征的地震响应特征方程包括：根据泰勒级数理论和广义线性反演理论，构建储层物性参数表征的地震响应特征广义方程；根据地震岩石物理建模理论，构建储层物性参数表征的地震响应特征方程。3.根据权利要求2所述的储层物性参数直接反演方法，其中，储层物性参数表征的地震响应特征广义方程为：其中，d为地震数据，W为子波，R为地震反射系数，e是地震数据中包含的噪音，V
m
、S
w
分别为储层的孔隙度、矿物含量和含水饱和度，ΔV
m
、ΔS
w
表示孔隙度、矿物含量和含水饱和度的扰动量，V
m0
、S
w0
表示孔隙度、矿物含量和含水饱和度的背景值，为背景物性参数V
m0
、S
w0
对应的地震反射系数，为反射系数对弹性参数的偏导数，V1、V2、V3表示三个不同的弹性参数，为反射系数对储层物性参数的偏导数，为弹性参数对储层物性参数的偏导数。4.根据权利要求2所述的储层物性参数直接反演方法，其中，储层物性参数表征的地震响应特征方程为：
其中，为反射系数对弹性参数的偏导数，V
Pi
,V
Pi+1
,V
Si
,V
Si+1
,ρ
i
,ρ
i+1
分别是不同层反射界面两侧的纵波速度、横波速度和密度，γ是纵横波速度比，为纵波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为纵波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为纵波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为纵波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为横波速度对孔隙度、白云岩含量和含水饱和度的偏导数，为横波速度对孔隙度、白云岩含量和含水饱和度的偏导数，
为密度对孔隙度、白云岩含量和含水饱和度的偏导数，饱和度的偏导数，K
sat
和μ
sat
分别为饱和岩石的体积模量和剪切模量；K
m
和μ
m
分别为岩石基质的体积模量和剪切模量，K
cal
和μ...

【专利技术属性】
技术研发人员：刘倩，陈天胜，
申请(专利权)人：中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：