[发明专利]一种基于炮检距矢量片的方位角道集提取方法和系统

权利要求书

1.一种基于炮检距矢量片的方位角道集提取方法，其特征在于，包括以下步骤：

S1将地震数据分为若干炮检距矢量片；

S2计算每个所述炮检距矢量片的地面方位角；

S3对每个所述炮检距矢量片中地震数据进行三维Kirchhoff积分叠前时间偏移，并计算出相应的地下散射方位角和反射角；

S4保留与所述地面方位角角度一致的所述地下散射方位角；

S5根据保留的所述地下散射方位角在一平面内选取若干x轴，所述x轴用于表征发射角的角度变化，所述平面的法线方向为z轴，所述z轴用于表征地层深度变化，将步骤S3中获得的所述反射角标记在x-z坐标系中得到反射角分布图谱；

S6根据所述反射角分布图谱获得方位角道集；

所述反射角分布图谱通过以下方式获得：首先将所述平面内的所述反射角在x轴上进行标记，再根据预设的地面炮检距和初始速度模型，计算不同深度的x轴方向上的反射角，并将其标记在所述x-z坐标系中，从而获得反射角分布图谱。

2.如权利要求1所述的基于炮检距矢量片的方位角道集提取方法，其特征在于，所述散射方位角的计算公式为：

其中，P_M是入射慢度矢量与散射慢度矢量之和；P_R是散射慢度矢量；P_S是入射慢度矢量，x，y分别为以散射点为原点建立的笛卡尔坐标的x轴和y轴。

3.如权利要求2所述的基于炮检距矢量片的方位角道集提取方法，其特征在于，所述反射角的计算公式为：

其中，P_R是散射慢度矢量；P_S是入射慢度矢量；γ是散射张角。

4.一种基于炮检距矢量片的方位角道集提取系统，其特征在于，包括：

分片模块，用于将地震数据分为若干炮检距矢量片；

地面方位角计算模块，用于计算每个所述炮检距矢量片的地面方位角；

散射方位角和反射角计算模块，用于对每一个所述炮检距矢量片中地震数据进行三维Kirchhoff积分叠前时间偏移，并计算出相应的散射方位角和反射角；

地下散射方位角选取模块，保留与所述地面方位角角度一致的所述地下散射方位角；

反射角分布图谱模块，用于根据保留的所述地下散射方位角在一平面内选取若干x轴，所述x轴用于表征发射角的角度变化，所述平面的法线方向为z轴，所述z轴用于表征地层深度变化，将步骤S3中获得的所述反射角标记在x-z坐标系中得到反射角分布图谱；

方位角道集获取模块，用于根据所述反射角分布图谱获得方位角道集；

所述反射角分布图谱通过以下方式获得：首先将所述平面内的所述反射角在x轴上进行标记，再根据预设的地面炮检距和初始速度模型，计算不同深度的x轴方向上的反射角，并将其标记在所述x-z坐标系中，从而获得反射角分布图谱。

5.如权利要求4所述的基于炮检距矢量片的方位角道集提取系统，其特征在于，所述散射方位角的计算公式为：

其中，P_M是入射慢度矢量与散射慢度矢量之和；P_R是散射慢度矢量；P_S是入射慢度矢量，x，y分别为以散射点为原点建立的笛卡尔坐标的x轴和y轴。

6.如权利要求5所述的基于炮检距矢量片的方位角道集提取系统，其特征在于，所述反射角的计算公式为：

其中，P_R是散射慢度矢量；P_S是入射慢度矢量，γ是散射张角。