【技术实现步骤摘要】
一种孔隙参数与储层参数的地震岩石物理解析反演方法
[0001]本专利技术属于非常规储层地震勘探
，特别涉及孔隙参数与储层参数的地震岩石物理解析反演方法。
技术介绍
[0002]地震岩石物理反演结合了地震反演技术与岩石物理模型，通过地震数据或属性定量预测储层参数，旨在为油气储层的预测、勘探及开发提供技术支撑。另一方面，随着油气勘探的不断推进，常规储层日益枯竭，致密砂岩、碳酸盐岩等非常规储层已称为增储上产的重点领域，但此类储层普遍发育复杂孔隙结构，严重制约了常规地震岩石物理反演方法的适用性。
[0003]常规地震岩石物理反演方法主要预测储层参数(孔隙度、含水饱和度等)，例如文献[1]‑
[4]公开的现有技术，但上述方法在岩石物理建模中通常将孔隙参数(孔隙纵横比)作为已知参数，并在反演过程中采用固定的孔隙纵横比。研究发现，岩石的孔隙纵横比对地震响应具有重要影响，尤其针对发育复杂孔隙结构的储层岩石[5]‑
[6]；此外，非常规储层由于地质构造及沉积环境等因素，储层岩石的孔隙结构异常复杂[7]‑
...
【技术保护点】
【技术特征摘要】
1.一种孔隙参数与储层参数的地震岩石物理解析反演方法，其特征在于，具体包括如下步骤：步骤一，采用地震道集数据，通过叠前地震反演，获取观测弹性参数模型；步骤二，设置基本岩石物理参数，基于测井数据，反演得到井旁孔隙纵横比；步骤三，计算模型参数的概率分布函数，构建模型参数的初始模型；步骤四，利用地震岩石物理正演算子模拟弹性参数，计算模拟弹性参数与观测弹性参数的误差；步骤五，基于初始模型，计算精确正演算子对模型参数的偏导矩阵，构建线性正演算子；步骤六，基于线性正演算子，解析计算模型参数的后验期望值；步骤七，更新模型参数的初始模型为后验期望值，重复步骤四至步骤六，直至模拟与观测弹性参数的误差降低到预设范围内，迭代停止，输出模型参数的后验期望值，即为孔隙参数与储层参数的反演结果。2.根据权利要求1所述的孔隙参数与储层参数的地震岩石物理解析反演方法，其特征在于，步骤一具体包括：获取工区PP波地震道集数据，通过叠前地震反演，构建目标函数其中，d为地震道集数据，m为弹性参数，为弹性参数的平均值，Σ
m
为弹性参数的协方差矩阵，G为Zoeppritz地震正演模型，σ
d
为地震数据噪声的标准差，T为转置符号；弹性参数m包括纵波速度、横波速度和密度；采用贝叶斯线性反演方法，得到满足式(1)目标函数的最优解，获取纵波速度、横波速度和密度参数模型。3.根据权利要求2所述的孔隙参数与储层参数的地震岩石物理解析反演方法，其特征在于，所述步骤二具体包括：设置基本岩石物理参数，包括岩石基质的体积模量、岩石基质的剪切模量、岩石基质的密度、烃类的体积模量、烃类的密度、卤水的体积模量和卤水的密度；基于测井观测的纵波速度与横波速度，通过匹配模型预测与井旁观测数据，反演得到井旁孔隙纵横比，构建目标函数其中，α为孔隙纵横比，σ1与σ2为误差项权系数，和分别为模型预测和井旁观测的纵波速度，和分别为模型预测和井旁观测的横波速度，T为转置符号；采用粒子群优化算法，得到满足式(2)目标函数的最优解，获取井旁孔隙纵横比参数。4.根据权利要求3所述的孔隙参数与储层参数的地震岩石物理解析反演方法，其特征在于，步骤三具体包括：
通过测井观测数据获得井旁孔隙度与井旁含水饱和度，结合反演得到的井旁孔隙纵横比，构建模型参数的混合高斯模型，其概率分布函数为；其中，z为模型参数，N
k
表示第k个高斯分布，和λ
k
分别为N
k
的期望、协方差和权系数，C为高斯分布个数；模型参数包括孔隙纵横比、孔隙度与含水饱和度；通过期望最大算法计算混合高斯模型，得到C个高斯分布的期望、协方差和权系数，再通过测井插值，得到模型参数的初始模型。5.根据权利要求4所述的孔隙参数与储层参数的地震岩石物理解析反演方法，其特征在于，步骤四具体包括：利用解耦微分等效介质模型计算岩石骨架的弹性模量K
d
(φ)＝K
m
(1
‑
φ)
p
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)μ
d
(φ)＝μ
m
(1
‑
φ)
q
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)其中，φ为孔隙度，K
d
和μ
d
分别为岩石骨架的体积模量和剪切模量，K
m
和μ
m
分别为岩石基质的体积模量和剪切模量，p和q为孔隙结构系数；基于岩石骨架的弹性模量，利用Gassmann方程计算含流体岩石的弹性模量μ
s
＝μ
d
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)其中，K
s
和μ
s
分别为含流体岩石的体积模量和剪切模量，K
f
为混合流体的体积模量，采用Wood模型计算其中，K
w
和K
h...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。