[发明专利]HTI介质AVAZ定量反演方法及系统

权利要求书

1.一种HTI介质AVAZ定量反演方法，包括：

步骤1：基于各向异性项B_ani和方位角建立包括地震道N个时间采样的初始模型m；

步骤2：根据对应于不同入射角的不同方位角的地震记录差，获得观测数据；

步骤3：根据方位角反射系数差矩阵与地震子波，获得地震合成记录，进而获得地震合成数据；

步骤4：根据所述观测数据与所述地震合成数据，建立反演目标函数F(m+Δm)并进行反演，其中，Δm是模型修正量；

步骤5：当F(m+Δm)不满足设定精度要求时，将m修正为m+Δm，重复步骤3-5；满足设定精度要求时，停止迭代，将迭代模型作为反演结果输出；

其中，入射角为θ_i时方位角反射系数差矩阵为：

其中，ΔR(θ_i)为入射角为θ_i时方位角反射系数差矩阵，为入射角为θ_i时方位角为和的反射系数差，θ_i为第i个入射角，i＝1，…，I，I为入射角的数量，K为方位角的数量；

其中，所述反演目标函数为：

其中，F(m+Δm)为反演目标函数，Δm为模型修正量，μ为阻尼系数，S_obs为观测数据。

2.根据权利要求1所述的HTI介质AVAZ定量反演方法，其中，所述初始模型m为：

其中，m_n为第n个采样时间点的初始模型：

其中，B_ani为各向异性项，为第n个采样时间点的方位角，n＝1，…，N，N为采样时间点的数量。

3.根据权利要求1所述的HTI介质AVAZ定量反演方法，其中，入射角为θ_i时的方位角记录差为：

其中，ΔS(θ_i)为入射角为θ_i时的方位角记录差，为入射角为θ_i时方位角为和的地震记录差，θ_i为第i个入射角，i＝1，…，I，I为入射角的数量，为第k个方位角，k＝1，…，K，K为方位角的数量。

4.根据权利要求3所述的HTI介质AVAZ定量反演方法，其中，所述观测数据为：

S_obs＝[ΔS(θ₁) ΔS(θ₂)… ΔS(θ_i)… ΔS(θ_I)] (4)

其中，S_obs为观测数据。

5.根据权利要求1所述的HTI介质AVAZ定量反演方法，其中，方位角和的反射系数差为：

其中，为第k个方位角，k＝1，…，K，K为方位角的数量，为自然坐标系下垂直裂缝的方位角，θ为入射角。

6.根据权利要求1所述的HTI介质AVAZ定量反演方法，其中，地震合成记录为：

ΔS_model(θ_i)＝W(θ_i)*ΔR(θ_i) (7)

其中，ΔS_model(θ_i)为入射角为θ_i时的地震合成记录，W(θ_i)为入射角为θ_i时的地震子波。

7.根据权利要求6所述的HTI介质AVAZ定量反演方法，其中，所述地震合成数据为：

S(m)＝[ΔS_model(θ₁) ΔS_model(θ₂)… ΔS_model(θ_I)] (8)

其中，S(m)为地震合成数据。

8.一种HTI介质AVAZ定量反演系统，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现以下步骤：

步骤1：基于各向异性项B_ani和方位角建立包括地震道N个时间采样的初始模型m；

步骤2：根据对应于不同入射角的不同方位角的地震记录差，获得观测数据；

步骤3：根据方位角反射系数差矩阵与地震子波，获得地震合成记录，进而获得地震合成数据；

步骤4：根据所述观测数据与所述地震合成数据，建立反演目标函数F(m+Δm)并进行反演，其中，Δm是模型修正量；

步骤5：当F(m+Δm)不满足设定精度要求时，将m修正为m+Δm，重复步骤3-5；满足设定精度要求时，停止迭代，将迭代模型作为反演结果输出；

其中，入射角为θ_i时方位角反射系数差矩阵为：

其中，ΔR(θ_i)为入射角为θ_i时方位角反射系数差矩阵，为入射角为θ_i时方位角为和的反射系数差，θ_i为第i个入射角，i＝1，…，I，I为入射角的数量，K为方位角的数量；

其中，所述反演目标函数为：

其中，F(m+Δm)为反演目标函数，Δm为模型修正量，μ为阻尼系数，S_obs为观测数据。