[发明专利]一种用于求取裂缝各向异性梯度的方法

摘要

本发明公开了一种用于求取裂缝各向异性梯度的方法，包括：基于采集的地震数据获取叠前方位道集角度信息；基于角度信息对叠前方位道集进行叠加处理以获取对应的叠后方位地震数据体；对叠后方位地震数据体进行稀疏脉冲反演以获取对应的反射系数体；基于简化的Ruger HTI介质纵波近似反射系数方程对反射系数体和角度信息进行处理以获取HTI介质上下两层纵波变异系数差值和相关变量；基于HTI介质上下两层纵波变异系数差值和相关变量计算各向异性梯度。本发明以叠后方位数据为基础，选取地震数据品质高的部分作为反演输入数据，避免了叠前地震资料对计算结果的直接影响，反演方法稳定性高、信噪比高。

主权项

1.一种用于求取裂缝各向异性梯度的方法，包括以下步骤：步骤一、基于采集的地震数据获取叠前方位道集的角度信息，所述角度信息包括叠前方位道集的方位角和入射角；步骤二、基于所述角度信息对所述叠前方位道集进行叠加处理以获取对应的叠后方位地震数据体；步骤三、对所述叠后方位地震数据体进行稀疏脉冲反演以获取对应的反射系数体；步骤四、基于简化的Ruger HTI介质纵波近似反射系数方程对所述反射系数体和所述角度信息进行处理以获取各向异性参数和相关变量；步骤五、基于所述各向异性参数和所述相关变量计算各向异性梯度ΔΓ；其中，在步骤四中，进一步包括：基于所述入射角信息获取叠前方位道集的最大入射角θ2和最小入射角θ1，并设定累加的间隔，将所述简化的Ruger HTI介质纵波近似反射系数方程按入射角方向累加以得到：其中， C = Σ θ 1 θ 2 sin 2 θ , ]]> D = 1 2 ( Δ ρ ρ ‾ + Δ α α ‾ ) + 1 2 { Δ α α ‾ - ( 2 β ‾ α ‾ ) 2 ( Δ ρ ρ ‾ + 2 Δ β β ‾ ) } sin 2 θ + 1 2 Δ α α ‾ sin 2 θtan 2 θ ; ]]>基于按入射角方向累加后的简化的Ruger HTI介质纵波近似反射系数方程，任选三个不同方位角及对应的反射系数体进行处理以得到所述HTI介质上下两层纵波变异系数差值Δδ(V)和所述相关变量同时，θ为入射角，为方位角，为HTI介质上下两层纵波速度平均值，为HTI介质上下两层横波速度平均值，为HTI介质上下两层密度平均值，Δα为HTI介质上下两层纵波速度差值，Δβ为HTI介质上下两层横波速度差值，Δρ为HTI介质上下两层密度差值，Δδ(V)为HTI介质上下两层纵波变异系数差值，Δγ为HTI介质上下两层横波各向异性差值。