一种地震叠前时间偏移中的射线追踪方法

摘要

本发明提供了一种地震叠前时间偏移中的射线追踪方法，包括：步骤1，获取M个反射界面的模型参数以及N个炮检距参数；步骤2，计算来自K个界面第j个炮检距的反射波传播路径初值参数a_i⁽¹⁾及传播时间步骤3，根据设定的迭代系数，利用迭代法计算来自K个界面第j个炮检距的新的反射波传播路径初值参数a_i⁽²⁾，并计算来自K个界面第j个炮检距的新的反射波传播时间步骤4，如果则即为来自K个界面第j个炮检距的反射波传播时间，否则，返回步骤3继续迭代；步骤5，如果j<N，则转向步骤2‑4继续计算，否则，来自K个界面反射波传播时间计算完成；步骤6，如果K<M，则转向步骤2‑5继续计算，否则，最终的反射波传播时间计算完成。

主权项

一种地震叠前时间偏移中的射线追踪方法，其特征在于，所述方法包括：步骤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继续计算，否则，最终的反射波传播时间计算完成；其中：在所述步骤1中，获取的M个反射界面的模型参数包括：各层层速度v_k以及各层厚度h_k，其中，k＝1,2,...,M；在所述步骤2中，首先依据下述公式计算第k＝1个界面的各个炮检距的反射波时间：其中，x_j为第k＝1个界面上的N个由小到大排列的炮检距参数，v_k为第k＝1个界面上的层速度，h_k为第k＝1个界面的层厚度；在所述步骤2中，根据K个界面的模型参数和炮检距参数，计算来自K个界面第j个炮检距的反射波传播路径初值参数a_i⁽¹⁾及传播时间包括：如果j＝1，则${a_i}^{\left(1\right)}=\frac{v_i{h}_i}{\underset{l=1}{\overset{M}{\Sigma }}{v}_l{h}_l}{x}_1,$否则，${a_i}^{\left(1\right)}={a_i}^{\left(1\right)}+\frac{v_i{h}_i}{\underset{l=1}{\overset{M}{\Sigma }}{v}_l{h}_l}(x_j-x_{j-1}),i=1,2,...,M;$以及，$t_{k,j}^{\left(1\right)}=\underset{i=1}{\overset{M}{\Sigma }}{({\left[a_i^{\left(1\right)}\right]}^2+4h_i^2)}^{1/2}/v_i;$在所述步骤3中，如果步骤2中的迭代次数小于最大迭代次数iter，则根据设定的迭代系数，利用迭代法计算来自K个界面第j个炮检距的新的反射波传播路径初值参数a_i⁽²⁾，包括：当v_i≥v_i+1时，i＝1,...,M‑1，迭代系数$c_i={(r_i^2+r_i^2\frac{{\left[a_i^{\left(1\right)}\right]}^2}{h_i^2}-\frac{{\left[a_i^{\left(1\right)}\right]}^2}{h_{i+1}^2})}^{1/2},r_i^2=\frac{v_i^2{h}_i^2}{v_{i+1}^2{h}_{i+1}^2};$当v_i＜v_i+1时，i＝1,...,M‑1，迭代系数$c_i={(r_i^2+r_i^2\frac{{\left[a_{i+1}^{\left(1\right)}\right]}^2}{h_{i+1}^2}-\frac{{\left[a_{i+1}^{\left(1\right)}\right]}^2}{h_i^2})}^{1/2},r_i^2=\frac{v_{i+1}^2{h}_{i+1}^2}{v_i^2{h}_i^2};$迭代系数d_i＝a_i⁽¹⁾+a_i+1⁽¹⁾；则根据设定的迭代系数c_i和d_i，计算来自K个界面第j个炮检距的新的反射波传播路径初值参数：${a_i}^{\left(2\right)}=\frac{c_i}{1+c_i}{d}_i,{a_{i+1}}^{\left(2\right)}=\frac{1}{1+c_i}{d}_i,i=1,...,M-1.$