[发明专利]一种多重孔隙储层叠前地震概率化多道反演方法

权利要求书

1.一种多重孔隙储层叠前地震概率化多道反演方法，其特征在于，其包括如下步骤：

步骤1，推导含多重孔隙空间的岩石弹性模量表达式；

步骤2，推导多重孔隙储层物性参数表征的地震反射系数方程；

步骤3，验证地震反射系数方程的近似精度和反演可行性；

步骤4，构建待反演模型参数的后验概率密度分布及目标泛函；

步骤5，研发多马尔科夫链随机采样的叠前地震多道逐级反演算法；

步骤6，研发基于逐步模拟策略的岩石物理参数反演方法。

2.如权利要求1所述多重孔隙储层叠前地震概率化多道反演方法，其特征在于，

所述步骤1包括：

根据比奥特-伽斯曼孔隙弹性理论，饱和岩石的体积模量和剪切模量的表达式为：

K_s≈K_d+[1-K_d/K_m]²φ^-1K_f (1)，

μ_s＝μ_d (2)，

式(1)、式(2)中，K_s和μ_s为饱和岩石的体积模量和剪切模量，K_d和μ_d为干燥岩石骨架的体积模量和剪切模量，K_m为岩石基质的体积模量，φ为孔隙度，K_f为孔隙流体的体积模量；

将科斯-徐模型用于计算干燥岩石骨架的体积模量和剪切模量：

K_d/K_m≈(1-φ)^p (3),

μ_d/μ_m≈(1-φ)^q (4),

式(3)、式(4)中，p和q表示与孔隙纵横比相关的等效极化因子，μ_m为岩石基质的剪切模量；

假设不同纵横比的孔隙空间划分为小纵横比的软孔隙和大纵横比的硬孔隙，等效极化因子p和q写成，

p＝α_sν_s+(1-α_s)ν_n (5),

式(5)、式(6)中，1-α_s和α_s分别是软孔隙体积分数和硬孔隙体积分数,软孔隙体积分数是软孔隙度与总孔隙度的比值，硬孔隙体积分数是硬孔隙度与总孔隙度的比值，系数ν_s和描述了硬孔隙对干燥岩石骨架的影响，系数ν_n和描述了软孔隙对干燥岩石骨架的影响；

令，G(φ)＝(1-φ)^p，L(φ)＝(1-φ)^q，将式(3)和式(4)代入式(1)、式(2)中，得到饱和岩石的体积模量和剪切模量的表达式，

K_s≈K_mG(φ)+[1-G(φ)]²φ^-1K_f (7)，

μ_s≈μ_mL(φ) (8)，

f≈[1-G(φ)]²φ^-1K_f (9)，

式(9)中，f为比奥特-伽斯曼理论中的伽斯曼流体项。