Eigenrays in 3D heterogeneous anisotropic media: Part VII -- Dynamics, Finite-element implementation (2003.11418v1)

Abstract: In this part, we apply the same finite-element approach, used in Part III for the vanishing first traveltime variation (to obtain the stationary rays), for the second traveltime variation, in order to compute the dynamic characteristics along the stationary ray. The finite-element solver involves application of the weak formulation and the Galerkin method to the linear second-order Jacobi ordinary differential equation (derived in Part V), yielding an original linear algebraic equation set for dynamic ray tracing. In our formulation, the resolving matrix of the linear equation set coincides with the global traveltime Hessian computed for the kinematic ray tracing, making the solution of the dynamic problem straightforward. The proposed method is unconditionally stable (the solution does not explode when the intervals between the nodes are increased) and is more accurate than the commonly used numerical integration (e.g., Runge-Kutta) methods, in particular, for stationary rays passing through heterogeneous anisotropic models with complex wave phenomena.