[发明专利]叠前各向异性特征参数反演方法及计算机可读存储介质

权利要求书

1.一种叠前各向异性特征参数反演方法，包括：

基于叠前地震道集，获得反射系数；

针对叠前地震数据进行叠前同时反演，分别获得上层介质与下层介质的纵波速度、横波速度与密度；

基于所述反射系数、所述上层介质与下层介质的纵波速度、横波速度与密度，根据各向异性参数反演公式进行反演，获得所述上层介质与下层介质的各向异性参数。

2.根据权利要求1所述的叠前各向异性特征参数反演方法，其中，基于叠前地震道集，获得反射系数包括：

将所述叠前地震道集划分为多个方位角道集；

基于速度模型将所述多个方位角道集转换为多个入射角道集；

沿目标层针对每个入射角道集提取反射系数。

3.根据权利要求1所述的叠前各向异性特征参数反演方法，其中，所述针对叠前地震数据进行叠前同时反演，分别获得上层介质与下层介质的纵波速度、横波速度与密度包括：

对所述叠前地震数据进行叠前同时反演，获得纵波速度体、横波速度体与密度体；

将目标层划分为所述上层介质与所述下层介质，在所述上层介质与所述下层介质中分别选定一个时窗，提取所述时窗所对应的纵波速度体、横波速度体与密度体，并分别计算所述对应的纵波速度体、横波速度体与密度体的均方根，作为所述上层介质与所述下层介质的纵波速度、横波速度与密度。

4.根据权利要求1所述的叠前各向异性特征参数反演方法，其中，所述各向异性参数反演公式为：

其中，R_p表示反射振幅，i表示入射角，φ表示方位角，δ表示纵波变化系数，Δδ^(V)表示上下层介质的纵波变化系数之差，ε表示纵波各向异性参数，Δ ε^(V)表示上下层介质的纵波各向异性参数之差，γ表示横波各向异性参数，Δγ表示上下层介质的横波各向异性参数之差，表示纵波速度，Δα表示上下层介质的纵波速度之差，表示横波速度，表示纵波垂直入射时的波阻抗，ΔZ表示纵波垂直入射时的波阻抗差，表示横波的切向模量，ΔG表示横波的切向模量差。

5.根据权利要求1所述的叠前各向异性特征参数反演方法，其中，还包括：基于所述上层介质与下层介质的各向异性参数，计算所述上层介质与下层介质的各向异性参数之差。

6.一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：

基于叠前地震道集，获得反射系数；

针对叠前地震数据进行叠前同时反演，分别获得上层介质与下层介质的纵波速度、横波速度与密度；

基于所述反射系数、所述上层介质与下层介质的纵波速度、横波速度与密度，根据各向异性参数反演公式进行反演，获得所述上层介质与下层介质的各向异性参数。

7.根据权利要求6所述的计算机可读存储介质，其中，基于叠前地震道集，获得反射系数包括：

将所述叠前地震道集划分为多个方位角道集；

基于速度模型将所述多个方位角道集转换为多个入射角道集；

沿目标层针对每个入射角道集提取反射系数。

8.根据权利要求6所述的计算机可读存储介质，其中，所述针对叠前地震数据进行叠前同时反演，分别获得上层介质与下层介质的纵波速度、横波速度与密度包括：

对所述叠前地震数据进行叠前同时反演，获得纵波速度体、横波速度体与密度体；

将目标层划分为所述上层介质与所述下层介质，在所述上层介质与所述下层介质中分别选定一个时窗，提取所述时窗所对应的纵波速度体、横波速度体与密度体，并分别计算所述对应的纵波速度体、横波速度体与密度体的均方根，作为所述上层介质与所述下层介质的纵波速度、横波速度与密度。