[发明专利]粘弹性流体因子直接反演方法

权利要求书

1.粘弹性流体因子直接反演方法，其特征是：以地下介质的粘弹性性质为基础，推导出粘弹性地震纵波反射近似表达式，建立地震纵波反射系数与粘弹性流体因子、剪切模量和密度的关系，利用弹性波阻抗反演方法，通过叠前不同角度道集资料的联合，直接反演粘弹性流体因子。

2.根据权利要求1所述的粘弹性流体因子直接反演方法，其特征是，

(1)粘弹性流体因子纵波反射系数方程建立包括：

弹性介质用流体因子表示的反射系数方程

其中，

代入粘弹性流体因子表达式

f_ela＝f_ane-Δf_Q(2)

其中，

经过数学推导，得到用粘弹性流体因子表示的纵波反射系数方程：

其中

假设干岩石纵横波速度比不受衰减影响，将饱和岩石的纵横波速度比代入反射系数表达式，得

忽略反射系数中的虚部项，则上式可以进一步简化为

该表达式即建立了地震纵波反射系数与粘弹性流体因子、剪切模量和密度三个参数之间的理论关系，基于该表达式，即可利用叠前地震资料直接反演粘弹性流体因子；

(2)粘弹性波阻抗方程建立包括：

按照弹性波阻抗公式，粘弹性波阻抗QEI进一步表达为含粘弹性流体因子的粘弹性阻抗：

其中，

(3)粘弹性波阻抗反演包括：

从叠前地震资料中得到至少三个不同角度的粘弹性弹性阻抗体，即：QEI(θ₁)、QEI(θ₂)、QEI(θ₃)，利用这三个粘弹性弹性阻抗体提取粘弹性流体因子、剪切模量和密度，将方程(6)的非线性模型进行线性变换得到如下模型

ln(QEI)＝a(θ)ln(f_ane)+b(θ)ln(μ)+c(θ)ln(ρ)(7)

为了得到ln(f_ane)、ln(μ)和ln(ρ)，需要三个不同角度的粘弹性阻抗体，将三个角度值分别代入理论模型表达式(7)得到如下的方程组：

方程(8)写成矩阵的形式为：

求解方程(9)。