[发明专利]一种计及弹性阻抗二阶梯度的储层参数反演方法

权利要求书

1.一种计及弹性阻抗二阶梯度的储层参数反演方法，包括根据叠前地震数据预测弹性阻抗，利用所述弹性阻抗反演预测储层参数，其特征在于，所述利用弹性阻抗反演预测储层参数具体为，根据所述弹性阻抗的一阶梯度和二阶梯度，反演预测储层参数；

所述储层参数的预测表达式为：

m＝m₀+γΔm

Δm＝-H^-1g

式中，m为未知数变量，即为所述储层参数的反演预测值，m₀为未知数变量的初始模型，γ为步长，Δm为未知数变量的迭代值，H为弹性阻抗的二阶梯度，g为弹性阻抗的一阶梯度；

所述弹性阻抗的二阶梯度的计算表达式为：

H≈gg^T

式中，T为转置符号；

所述弹性阻抗由密度、修正孔隙度及流体体积模量表征，所述弹性阻抗的一阶梯度的计算表达式为：

式中，EI为弹性阻抗的预测值，ρ为密度，φ为岩石的孔隙度，K_f为流体的体积模量，Δd为反演得到的EI和模拟得到的EI之间差异变量；

所述计及弹性阻抗二阶梯度的储层参数反演方法面向内幕型孔缝储层，所述弹性阻抗的预测表达式为：

式中，EI为弹性阻抗的预测值，ρ为密度，θ为入射角，g_sat为饱和岩石的剪切模量和纵波模量比值，g_dry为干岩石的剪切模量和纵波模量比值，φ_n为修正的孔隙度，a为修正的孔隙度与岩石孔隙度之间线性关系的梯度值，K_f为流体的体积模量；

所述弹性阻抗的预测表达式根据所述储层的反射系数近似式建立，所述储层的反射系数近似式根据所述储层的饱和岩石体积模量计算式建立；

所述饱和岩石体积模量计算式为：

式中，K_sat为饱和岩石的体积模量，K_dry为干岩石的体积模量，K_f为流体的体积模量，K₀为岩石矿物的体积模量，φ为岩石的孔隙度；

所述计及弹性阻抗二阶梯度的储层参数反演方法面向含气储层，所述饱和岩石体积模量计算式为：

式中，K_sat为饱和岩石的体积模量，K_dry为干岩石的体积模量，K_f为流体的体积模量，φ_n为修正的孔隙度，φ为岩石的孔隙度；

所述储层的反射系数近似式的表达式为：

p_μ(θ)＝-2g_satsin²θ

式中，R_PP为反射系数近似计算结果，θ为入射角，g_sat为饱和岩石的剪切模量和纵波模量比值，g_dry为干岩石的剪切模量和纵波模量比值，M₀为岩石矿物的纵波模量，μ₀为岩石矿物的剪切模量，ρ为密度，φ_n为修正的孔隙度，a为修正的孔隙度与岩石孔隙度之间线性关系的梯度值，K_f为流体的体积模量；

所述计及弹性阻抗二阶梯度的储层参数反演方法面向内幕型孔缝储层，所述储层的反射系数近似式的表达式为：