[发明专利]基于叠前叠后联合反演的井段气层识别方法及系统

权利要求书

1.一种井段气层识别方法，其特征在于，包括如下步骤：

获取井段的叠前地震数据，所述叠前地震数据包括二维地震剖面；

获取二维地震剖面中井段的角度道集；

从所述角度道集中选择角度子波；

对各角度子波进行叠前地震反演，得到对应的横波速度；

将所述横波速度与测井数据进行吻合，根据吻合结果选择最优角度子波；

通过Zoeppritz方程计算井段的反射系数；

对所述反射系数、最优角度子波进行褶积，得到直接反演横波速度；

根据所述直接反演横波速度和约束公式计算得到纵横波速度比；

根据所述纵横波速度比和二维地震剖面中单个测井的含气饱和度计算二维地震剖面的含气饱和度；

通过所述二维地震剖面的含气饱和度对气层进行识别。

2.根据权利要求1所述的一种井段气层识别方法，其特征在于，所述Zoeppritz方程为：

其中，R_PP代表纵波入射纵波反射的反射系数，R_PS代表纵波入射横波反射的反射系数，T_PP代表纵波入射纵波反射的透射系数，T_PS代表纵波入射横波反射的透射系数，v代表速度，ρ₁代表上覆介质的密度，ρ₂代表下伏介质的密度，v_s1代表上覆介质的横波速度，v_s2代表下伏介质的横波速度，v_p1代表上覆介质的纵波速度，v_p2代表下伏介质的纵波速度，θ₁和θ₂分别代表纵波的反射角和透射角，和分别代表横波的反射角和透射角。

3.根据权利要求1所述的一种井段气层识别方法，其特征在于，所述约束公式为：

V_S＝1775.2+0.08818×V_P，

其中，V_S为横波速度，V_P为纵波速度。

4.根据权利要求3所述的一种井段气层识别方法，其特征在于，所述纵波速度的范围为3200.0-4100.0，所述横波速度的范围为1600.0-2300.0，所述横波速度偏差的范围为-200.0-200.0。

5.根据权利要求1所述的一种井段气层识别方法，其特征在于，所述角度道集的角度范围为3～45°。

6.根据权利要求5所述的一种井段含气饱和度计算方法，其特征在于，所述角度子波包括3-9°范围的角度子波、18-24°范围的角度子波和33-39°范围的角度子波。

7.一种井段气层识别系统，其特征在于，所述系统包括：

第一获取模块：用于获取井段的叠前地震数据，所述叠前地震数据包括二维地震剖面；

第二获取模块：用于获取二维地震剖面中井段的角度道集；

选择模块：用于从所述角度道集中选择角度子波；

反演模块：用于对各角度子波进行叠前地震反演，得到对应的横波速度；

吻合模块：用于将所述横波速度与测井数据进行吻合，根据吻合结果选择最优角度子波；

第一计算模块：用于通过Zoeppritz方程计算井段的反射系数；

褶积模块：用于对所述反射系数、最优角度子波进行褶积，得到直接反演横波速度；

第二计算模块：用于根据所述直接反演横波速度和约束公式计算得到纵横波速度比；

第三计算模块：用于根据所述纵横波速度比和二维地震剖面中单个测井的含气饱和度计算二维地震剖面的含气饱和度；

识别模块：用于通过所述二维地震剖面的含气饱和度对气层进行识别。

8.一种井段气层识别系统，其特征在于，所述系统包括处理器和存储介质；

所述存储介质用于存储指令；

所述处理器用于根据所述指令进行操作以执行根据权利要求1-6任一项所述方法的步骤。

9.计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1-6任一项所述方法的步骤。