[发明专利]一种射线参数域无约束弹性参数反演方法

权利要求书

1.一种射线参数域无约束弹性参数反演方法，其特征在于：所述方法将偏移距域下的地震数据转换到射线域中，通过构建射线域下振幅变化的关系并进行模型迭代反演，得到纵波速度、横波速度以及密度。

2.根据权利要求1所述的射线参数域无约束弹性参数反演方法，其特征在于：所述方法包括以下步骤：

(1)创建射线域下振幅变化的关系，即Zoeppritz方程；

(2)将偏移距下的地震数据转换到射线域形成射线域地震数据；

(3)对于每个射线域下的地震数据进行子波估算；

(4)根据测井资料建立初始模型；

(5)对步骤(2)得到的射线域地震数据在所述初始模型的基础上进行模型迭代反演，当前后两次结果收敛时，或者迭代次数超过预先设定次数后停止迭代反演，最终得到的模型值即为反演结果，包括纵波速度、横波速度以及密度；

(6)输出纵波速度、横波速度以及密度。

3.根据权利要求2所述的射线参数域无约束弹性参数反演方法，其特征在于：所述步骤(1)中的Zoeppritz方程如下：

-pα1-1-p2γ32α12pγ2α1-1-p2γ42α121-pγ12α12-pγ4α11-p2γ22α12pγ4α12pα1·1-pγ12α121-2(pγ3α1)2γ32pγ1α1γ42γ321-p2γ22α12-γ1·γ4γ32(1-2p2γ42α12)2(pγ3α1)2-12pγ32α11-p2γ32α12γ1γ2(1-2p2γ42α12)2γ12γ42pγ4α11-p2γ42α12RppRpsTppTps=pα11-pγ12α122pα11-pγ12α121-2(pγ3α1)2---(1)]]>

其中p代表射线参数，γ₁，γ₂，γ₃和γ₄的定义如下：

γ₁＝ρ₂/ρ₁，γ2=Vp2/Vp1,]]>γ3=Vs1/Vp1,]]>γ4=Vs2/Vp1---(2)]]>

其中V_p1、V_s1、ρ₁、V_p2、V_s2、ρ₁和ρ₂分别是速度和密度，下标1代表上层，2代表下层。