[发明专利]一种利用地震波阻抗定量反演孔隙度的方法

权利要求书

1.一种利用地震波阻抗定量反演孔隙度的方法，其特征在于：所述方法以地震波阻抗为输入数据，利用岩石物理方程和地震波传播速度公式构建孔隙度反演的误差函数，并在孔隙度的取值范围内构建解的可行域空间，在该可行域空间内采用非线性全局寻优算法进行迭代反演，得到最终的孔隙度。

2.根据权利要求1所述的利用地震波阻抗定量反演孔隙度的方法，其特征在于：所述方法包括：

(1)输入终止条件的值ε，在孔隙度的取值范围内随机给定孔隙度的初始值；

(2)输入固体矿物的体积模量K₀和剪切模量μ₀，输入固结系数c，计算干岩石弹性模量；

(3)计算饱和岩石弹性模量；

(4)由所述饱和岩石弹性模量计算波阻抗，所述波阻抗包括纵波阻抗I_p和横波阻抗I_s；

(5)输入比重因子λ和由叠前地震反演得到的纵波阻抗I_pi和横波阻抗I_si’建立孔隙度反演的误差函数；

(6)利用非线性全局寻优算法对所述孔隙度反演的误差函数进行迭代反演，直到所述孔隙度反演的误差函数满足终止条件；

(7)输出满足终止条件的孔隙度。

3.根据权利要求2所述的利用地震波阻抗定量反演孔隙度的方法，其特征在于：所述步骤(2)计算干岩石弹性模量是这样实现的：

所述干岩石弹性模量包括体积模量和剪切模量，分别由公式(1)和(2)计算得到：

K_dry=K₀(1-φ)／(1+cφ)  (1)

μ_dry=μ₀(1-φ)／(1+1.5cφ)  (2)

其中，K_dry和μ_dry分别为干岩石体积模量和剪切模量，φ为孔隙度。

4.根据权利要求3所述的利用地震波阻抗定量反演孔隙度的方法，其特征在于：所述步骤(3)是这样实现的：

利用岩石物理方程计算饱和岩石弹性模量，所述饱和岩石弹性模量包括体积模量和剪切模量，所述岩石物理方程为公式(3)和公式(4)：

Ksat=Kdry+(1-Kdry/K0)2/[φ/Kf+(1-φ)/K0-Kdry/K02]---(3)]]>

μ_sat=μ_dry  (4)

其中，K_f为孔隙饱含流体的体积模量，K_sat和μ_sat为饱和岩石体积模量和剪切模量。