[发明专利]一种地震叠前时间偏移中的射线追踪方法

权利要求书

1.一种地震叠前时间偏移中的射线追踪方法，其特征在于，所述方法包括：

步骤1，获取M个反射界面的模型参数以及N个炮检距参数；

步骤2，根据K个界面的模型参数和炮检距参数，计算来自K个界面第j个炮检距的反射波传播路径初值参数a_i⁽¹⁾及传播时间i=1,2,...,M，k=1,2,...,M，j=1,2,...,N；

步骤3，根据设定的迭代系数，利用迭代法计算来自K个界面第j个炮检距的新的反射波传播路径初值参数a_i⁽²⁾，并计算来自K个界面第j个炮检距的新的反射波传播时间

步骤4，如果ε为设定的误差，则即为来自K个界面第j个炮检距的反射波传播时间，否则，返回所述步骤3继续迭代；

步骤5，如果j<N，则转向步骤2-4继续计算，否则，来自K个界面反射波传播时间计算完成；

步骤6，如果K<M，则转向步骤2-5继续计算，否则，最终的反射波传播时间计算完成。

2.根据权利要求1所述的地震叠前时间偏移中的射线追踪方法，其特征在于，在所述步骤1中，获取的M个反射界面的模型参数包括：各层层速度v_k以及各层厚度h_k，其中，k=1,2,...,M。

3.根据权利要求2所述的地震叠前时间偏移中的射线追踪方法，其特征在于，在所述步骤2中，首先依据下述公式计算第k=1个界面的各个炮检距的反射波时间：

t1,j(2)=(xj2+4h12)1/2/v1,j=1,2,...,N;]]>

其中，x_j为第k=1个界面上的N个由小到大排列的炮检距参数，v_k为第k=1个界面上的层速度，h_k为第k=1个界面的层厚度。