Eigenrays in 3D heterogeneous anisotropic media,Part I: Kinematics

We present a new ray bending approach, referred to as the Eigenray method, for solving two‐point boundary‐value kinematic and dynamic ray tracing problems in 3D smooth heterogeneous general anisotropic elastic media. The proposed Eigenray method is aimed to provide reliable stationary ray path solutions and their dynamic characteristics, in cases where conventional initial‐value ray shooting methods, followed by numerical convergence techniques, become challenging. The kinematic ray bending solution corresponds to the vanishing first traveltime variation, leading to a stationary path between two fixed endpoints (Fermat's principle), and is governed by the nonlinear second‐order Euler–Lagrange equation. The solution is based on a finite‐element approach, applying the weak formulation that reduces the Euler–Lagrange second‐order ordinary differential equation to the first‐order weighted‐residual nonlinear algebraic equation set. For the kinematic finite‐element problem, the degrees of freedom are discretized nodal locations and directions along the ray trajectory, where the values between the nodes are accurately and naturally defined with the Hermite polynomial interpolation. The target function to be minimized includes two essential penalty (constraint) terms, related to the distribution of the nodes along the path and to the normalization of the ray direction. We distinguish between two target functions triggered by the two possible types of stationary rays: a minimum traveltime and a saddle‐point solution (due to caustics). The minimization process involves the computation of the global (all‐node) traveltime gradient vector and the traveltime Hessian matrix. The traveltime Hessian is used for the minimization process, analysing the type of the stationary ray, and for computing the geometric spreading of the entire resolved stationary ray path. The latter, however, is not a replacement for the dynamic ray tracing solution, since it does not deliver the geometric spreading for intermediate points along the ray, nor the analysis of caustics. Finally, we demonstrate the efficiency and accuracy of the proposed method along three canonical examples.