基于相速度的横向各向同性介质三维射线追踪方法

针对传统各向异性射线追踪需反复做刚度张量旋转与Christoffel本征系统计算的问题，提出一种直接由qP波相速度及其空间与方向导数构建射线方程的三维追踪方法，避免Bond变换，并采用四阶Runge–Kutta格式积分。均质、层状与构造复杂模型的数值试验结果与解析解、快速行进法结果及模拟qP波前吻合良好。