[发明专利]层析偏移速度分析中的射线追踪方法及装置

权利要求书

1.一种层析偏移速度分析中的射线追踪方法，其特征在于，包括：

获取角度域共成像点道集；所述角度域共成像点道集中的横坐标具有相同的第一角度间隔；

对于所述角度域共成像点道集中深度大于深度阈值的每个反射点，根据反射点对应的深度、最小反射角、最大反射角和反射角及深度阈值处反射点对应的深度、最小反射角、最大反射角和反射角，获得反射点对应的第二角度间隔；

根据每个反射点对应的第二角度间隔，对应调整所述角度域共成像点道集中的横坐标间隔，以使调整后的角度域共成像点道集中不同深度的反射点对应相同数量的反射射线条数；

根据调整后的角度域共成像点道集进行射线追踪；其中，所述根据反射点对应的深度、最小反射角、最大反射角和反射角及深度阈值处反射点对应的深度、最小反射角、最大反射角和反射角，获得反射点对应的第二角度间隔，包括：

根据公式获取反射点对应的第二角度间隔；

其中，θ^max为大于深度阈值的某个反射点对应的最大反射角，θ^min为大于深度阈值的某个反射点对应的最小反射角，为深度阈值处反射点对应的最大反射角，为深度阈值处反射点对应的最小反射角，θ为大于深度阈值的某个反射点对应的反射角，θ_thsh为深度阈值处反射点对应的反射角，Δθ为大于深度阈值的某个反射点对应的角度间隔，Δθ_thsh深度阈值处反射点对应的角度间隔。

2.如权利要求1所述的层析偏移速度分析中的射线追踪方法，其特征在于，所述获取角度域共成像点道集，包括：

利用速度模型对叠前地震数据进行叠前深度偏移，得到角度域共成像点道集。

3.如权利要求1所述的层析偏移速度分析中的射线追踪方法，其特征在于，所述根据调整后的角度域共成像点道集进行射线追踪，包括：

确定当前速度模型中反射波射线路径的第一旅行时，并确定对应的真实速度模型中反射波射线路径的第二旅行时；

根据所述第一旅行时和所述第二旅行时的差值建立层析方程组，并求解所述层析方程组，并根据求解结果更新所述当前速度模型。

4.如权利要求3所述的层析偏移速度分析中的射线追踪方法，其特征在于，所述确定当前速度模型中反射波射线路径的第一旅行时，包括：

根据公式确定当前速度模型中反射波射线路径的第一旅行时；

其中，Γ(θ₀)为反射角为θ₀射线的从震源到反射点再到观测点的射线路径，dl为射线路径l的微分，s₀为当前速度模型中的介质慢度，t^cal为当前速度模型中反射波射线路径的第一旅行时。

5.如权利要求3所述的层析偏移速度分析中的射线追踪方法，其特征在于，所述确定对应的真实速度模型中反射波射线路径的第二旅行时，包括：

根据公式确定当前速度模型中反射波射线路径的第二旅行时；

其中，Γ(θ_t)为反射角为θ_t射线的从震源到反射点再到观测点的射线路径，dl为射线路径l的微分，s_t为真实速度模型中的介质慢度，t^obs为真实速度模型中反射波射线路径的第二旅行时。

6.如权利要求3所述的层析偏移速度分析中的射线追踪方法，其特征在于，假定同一反射角度射线在真实速度中射线路径和当前速度中射线路径相同，当前反射点建立的层析方程组包括

其中，j为当前反射点第j条射线，为当前速度模型中反射角为θ^min+jΔθ的射线从震源到当前反射点再到观测点的射线路径，Δs_i为速度模型离散化后第i个网格点处当前速度模型与真实速度模型慢度之差，l_i,j为第j条射线在当前离散化速度模型第i个网格中的射线长度，Δt_j为当前反射点第j条射线第一旅行时与第二旅行时之差。