Fermat's Variational Principle for Anisotropic Inhomogeneous Media

Fermat's variational principle states that the signal propagates from point S to R along a curve which renders Fermat's functional (l) stationary. Fermat's functional (l) depends on curves l which connect points S and R, and represents the travel times from S to R along l. In seismology, it is mostly expressed by the integral (l) = (x^k,x^k')du, taken along curve l, where (x^k,x^k') is the relevant Lagrangian, x^k are coordinates, u is a parameter used to specify the position of points along l, and x^k' = dx^k÷du. If Lagrangian (x^k,x^k') is a homogeneous function of the first degree in x^k', Fermat's principle is valid for arbitrary monotonic parameter u. We than speak of the first-degree Lagrangian ⁽¹⁾(x^k,x^k'). It is shown that the conventional Legendre transform cannot be applied to the first-degree Lagrangian ⁽¹⁾(x^k,x^k') to derive the relevant Hamiltonian ⁽¹⁾(x^k,p_k), and Hamiltonian ray equations. The reason is that the Hessian determinant of the transform vanishes identically for first-degree Lagrangians ⁽¹⁾(x^k,x^k'). The Lagrangians must be modified so that the Hessian determinant is different from zero. A modification to overcome this difficulty is proposed in this article, and is based on second-degree Lagrangians ⁽²⁾. Parameter u along the curves is taken to correspond to travel time , and the second-degree Lagrangian ⁽²⁾(x^k,^k) is then introduced by the relation ⁽²⁾(x^k,^k) = [⁽¹⁾(x^k,^k)]², with ^k = dx^k÷d. The second-degree Lagrangian ⁽²⁾(x^k,^k) yields the same Euler/Lagrange equations for rays as the first-degree Lagrangian ⁽¹⁾(x^k,^k). The relevant Hessian determinant, however, does not vanish identically. Consequently, the Legendre transform can then be used to compute Hamiltonian ⁽²⁾(x^k,p_k) from Lagrangian ⁽²⁾(x^k,^k), and vice versa, and the Hamiltonian canonical equations can be derived from the Euler-Lagrange equations. Both ⁽²⁾(x^k,^k) and ⁽²⁾(x^k,p_k) can be expressed in terms of the wave propagation metric tensor g_ij(x^k,^k), which depends not only on position x^k, but also on the direction of vector ^k. It is defined in a Finsler space, in which the distance is measured by the travel time. It is shown that the standard form of the Hamiltonian, derived from the elastodynamic equation and representing the eikonal equation, which has been broadly used in the seismic ray method, corresponds to the second-degree Lagrangian ⁽²⁾(x^k,^k), not to the first-degree Lagrangian ⁽¹⁾(x^k,^k). It is also shown that relations ⁽²⁾(x^k,^k) = ; and ⁽²⁾(x^k,p_k) = are valid at any point of the ray and that they represent the group velocity surface and the slowness surface, respectively. All procedures and derived equations are valid for general anisotropic inhomogeneous media, and for general curvilinear coordinates xⁱ. To make certain procedures and equations more transparent and objective, the simpler cases of isotropic and ellipsoidally anisotropic media are briefly discussed as special cases.