[发明专利]真地表射线追踪自适应插值方法

权利要求书

1.一种真地表射线追踪自适应插值方法，包括以下步骤：

(a)将需要进行旅行时射线追踪的空间划分为由多个空间网格构成的深度域的空间网格体；

(b)从真地表出发进行射线追踪计算，得到地震波传播中的多条旅行时射线；

(c)确定每个空间网格中穿过的旅行时射线的数量；

(d)针对穿过的旅行时射线的数量少于预定数量的空间网格，选择从与该空间网格相邻的一个或多个空间网格中穿过的三条旅行时射线，并基于预先设置的旅行时射线步长构造三维射线五面体；

(e)通过位于该射线五面体的六个顶点处的控制点的旅行时来进行自适应内插，从而获得穿过该空间网格的一条旅行时射线；

(f)反复执行操作(d)和(e)，直到穿过该空间网格的旅行时射线的数量达到所述预定数量为止。

2.根据权利要求1所述的真地表射线追踪自适应插值方法，还包括以下步骤：

(g)针对穿过的旅行时射线的数量超过所述预定数量的空间网格，舍去部分旅行时射线，以使穿过相应空间网格的旅行时射线的数量等于所述预定数量。

3.根据权利要求2所述的真地表射线追踪自适应插值方法，还包括以下步骤：

(h)确定是否对所有空间网格执行了舍去步骤或插值步骤；

(i)如果没有对所有空间网格执行了舍去或插值操作，则重复执行步骤(d)、(e)、(f)、(g)，直到对所有空间网格执行了舍去或插值操作。

4.根据权利要求1所述的真地表射线追踪自适应插值方法，其中，所述预定数量为3。

5.根据权利要求1所述的真地表射线追踪自适应插值方法，其中，所述射线步长是10毫秒至40毫秒中的任意值。

6.根据权利要求1所述的真地表射线追踪自适应插值方法，其中，在步骤(e)中，基于控制点的旅行时通过克里金插值方法获得该空间网格内的一个点的旅行时，从而获得穿过该空间网格的一条旅行时射线。

7.根据权利要求1所述的真地表射线追踪自适应插值方法，其中，空间网格的边长是5米。