[发明专利]一种基于自适应混合骨架参数的高精度岩石物理建模方法

权利要求书

1.一种基于自适应混合骨架参数的高精度岩石物理建模方法，其特征在于，所述高精度岩石物理建模方法包括如下步骤：

S1：利用Krief经验公式表征混合骨架基质模量K₀,μ₀与干岩石基质模量K_dry,μ_dry的关系式；

S2：将步骤S1获取的关系式和Wood公式计算流体体积模量K_f代入Gassmann方程和纵波速度计算公式推导得到以混合骨架基质模量K₀,μ₀为反演模型，实测纵波与预测纵波误差最小的目标函数；

S3：利用Levenberg-Maquardt算法迭代反演计算步骤S2得到的目标函数，得到混合骨架基质模量K₀,μ₀；

S4：将步骤S3得到的混合骨架基质模量K₀,μ₀代入步骤S1混合骨架基质模量K₀,μ₀与干岩石基质模量K_dry,μ_dry的关系式中，计算得到干岩石基质模量K_dry,μ_dry；

S5：将步骤S3得到的混合骨架基质模量K₀,μ₀和步骤S4得到的干岩石基质模量K_dry,μ_dry代入Gassmann方程计算得到饱含流体岩石模量K_sat,μ_sat，最后利用饱含流体岩石模量K_sat,μ_sat计算得到横波速度V_s；

所述目标函数的计算步骤为：

μ_sat＝μ_dry (2-3)

所述混合骨架基质模量K₀,μ₀与干岩石基质模量K_dry,μ_dry的关系的Krief经验公式为：

其中：φ是孔隙度，K₀,μ₀是混合骨架基质模量，K_dry,μ_dry是干岩石模量；

A1：将(1-1)、(2-1)、(2-2)和(2-3)代入(2-4)，得到(2-5)，其中a是一个中间变量；

A2：令求偏导得：

A3：建立目标函数：

其中：K_sat,μ_sat是饱含流体岩石模量，K_f是混合流体体积模量，S_w、S_o、S_g、K_w、K_o、K_g分别是实测含水饱和度、含油饱和度、含气饱和度、水体积模量、气体积模量、油体积模量，ρ是实测密度值，λ₁是K₀和μ₀的平滑约束因子，取默认值0.001。

2.如权利要求1所述的一种基于自适应混合骨架参数的高精度岩石物理建模方法，其特征在于，所述计算混合骨架基质模量K₀,μ₀通过利用Levenberg-Maquardt算法进行二元迭代同时反演求解得到的目标函数，具体包括如下步骤：

设反演目标m＝(K₀,μ₀)初始值；

B1：计算初始的f(K₀,μ₀)与实测的差值e，

B2：利用式(2-6)建立偏导矩阵利用B1步骤中的扰动量e计算K₀,μ₀的扰动量Δm，Δm＝(G^TG+λ₁I)^-1G^Te；

B3：更新m＝m+Δm；

B4：计算累积误差并判断累积误差是否小于预设误差门槛值，如果小于门槛值，则反演计算结束，如果大于门槛值，则重复执行B1-B4步骤。

3.如权利要求2所述的一种基于自适应混合骨架参数的高精度岩石物理建模方法，其特征在于，所述利用饱含流体岩石模量K_sat,μ_sat计算得到横波速度V_s，具体为：

将获得的饱含流体岩石模量K_sat,μ_sat代入公式：计算得到横波速度V_s。