水平层状VTI介质高效射线追踪计算方法

准确、快速计算水平层状VTI(vertical transverse isotropic)介质的走时和射线路径是实现VTI介质地震叠前时间偏移和微地震监测速度建模的基础.水平层状VTI介质射线追踪的常用算法包括打靶法和最短路径法.打靶法原理简单、实现容易,但在长偏移距时存在微小初射角度变化导致射线路径巨大扰动的问题;最短路径法原理直观,但需要沿射线群速度方向(群角)计算走时,计算特定群角对应的群速度值需要先搜索对应的相速度角度(相角),显著降低了计算效率.本文综合打靶法和最短路径法的优点,从震源同时出射一系列射线,以这些射线到达界面的交点构成稀疏的网格节点,自适应加密实现检波点周围网格节点间距高于精度要求,采用插值算法获得检波点的走时和动态网格方法获得不位于网格节点上的检波点的射线路径.本文方法严格计算已知相角对应的群角和群速度,未使用弱各向异性群速度近似公式,适用于任意强度各向异性VTI介质qP、qSV和qSH直达波和反射波射线追踪;以相角确定震源出射射线,不用遍历群角和群速度对应关系.同传统最短路径方法的数值实验对比表明本文方法具有高精度和高效性,非常适合于需要多次正演计算的地震叠前时间偏移和微地震实时监测问题.