Fermat's Variational Principle for Anisotropic Inhomogeneous Media

**Author: **Červený, V.

**Source:** Studia Geophysica et Geodaetica,
Volume 46, Number 3, July 2002
, pp. 567-588(22)

**Publisher: **Springer

**Abstract:**

*l*) stationary. Fermat's functional {\mathcal{I}}(

*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 {\mathcal{I}}(

*l*) = \int\limits_S^Rℒ(

*x*,

^{k}*x*)d

^{k}'*u*, taken along curve

*l*, where ℒ(

*x*,

^{k}*x*) is the relevant Lagrangian,

^{k}'*x*are coordinates,

^{k}*u*is a parameter used to specify the position of points along

*l*, and

*x*= d

^{k}'*x*÷d

^{k}*u*. If Lagrangian ℒ(

*x*,

^{k}*x*) is a homogeneous function of the first degree in

^{k}'*x*, Fermat's principle is valid for arbitrary monotonic parameter

^{k}'*u*. We than speak of the first-degree Lagrangian ℒ

^{(1)}(

*x*,

^{k}*x*). It is shown that the conventional Legendre transform cannot be applied to the first-degree Lagrangian ℒ

^{k}'^{(1)}(

*x*,

^{k}*x*) to derive the relevant Hamiltonian ℋ

^{k}'^{(1)}(

*x*,

^{k}*p*), and Hamiltonian ray equations. The reason is that the Hessian determinant of the transform vanishes identically for first-degree Lagrangians ℒ

_{k}^{(1)}(

*x*,

^{k}*x*). 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 ℒ

^{k}'^{(2)}. Parameter

*u*along the curves is taken to correspond to travel time

*τ*, and the second-degree Lagrangian ℒ

^{(2)}(

*x*,{\dot x}

^{k}^{k}) is then introduced by the relation ℒ

^{(2)}(

*x*,{\dot x}

^{k}^{k}) = ½[ℒ

^{(1)}(

*x*,{\dot x}

^{k}^{k})]

^{2}, with {\dot x}

^{k}= d

*x*÷d

^{k}*τ*. The second-degree Lagrangian ℒ

^{(2)}(

*x*,{\dot x}

^{k}^{k}) yields the same Euler/Lagrange equations for rays as the first-degree Lagrangian ℒ

^{(1)}(

*x*,{\dot x}

^{k}^{k}). The relevant Hessian determinant, however, does not vanish identically. Consequently, the Legendre transform can then be used to compute Hamiltonian ℋ

^{(2)}(

*x*,

^{k}*p*) from Lagrangian ℒ

_{k}^{(2)}(

*x*,{\dot x}

^{k}^{k}), and vice versa, and the Hamiltonian canonical equations can be derived from the Euler-Lagrange equations. Both ℒ

^{(2)}(

*x*,{\dot x}

^{k}^{k}) and ℋ

^{(2)}(

*x*,

^{k}*p*) can be expressed in terms of the wave propagation metric tensor

_{k}*g*(

_{ij}*x*,{\dot x}

^{k}^{k}), which depends not only on position

*x*, but also on the direction of vector {\dot x}

^{k}^{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 ℒ

^{(2)}(

*x*,{\dot x}

^{k}^{k}), not to the first-degree Lagrangian ℒ

^{(1)}(

*x*,{\dot x}

^{k}^{k}). It is also shown that relations ℒ

**Keywords:**
Fermat's principle;
Finsler space;
Hamiltonian;
Lagrangian;
anisotropic media;
wave propagation metric tensor

**Document Type:** Research Article

**Affiliations:**
Department of Geophysics, Charles University, Ke Karlovu 3, 121 16 Praha 2, Czech Republic vcerveny@seis.karlov.mff.cuni.cz

Publication date: July 1, 2002

**In this:**publication**By this:**publisher**By this author:**Červený, V.