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

Buy & download fulltext article:


Price: $47.00 plus tax (Refund Policy)


Fermat's variational principle states that the signal propagates from point S to R along a curve which renders Fermat's functional {\mathcal{I}}(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ℒ(xk,xk')du, taken along curve l, where ℒ(xk,xk') is the relevant Lagrangian, xk are coordinates, u is a parameter used to specify the position of points along l, and xk' = dxk÷du. If Lagrangian ℒ(xk,xk') is a homogeneous function of the first degree in xk', Fermat's principle is valid for arbitrary monotonic parameter u. We than speak of the first-degree Lagrangian ℒ(1)(xk,xk'). It is shown that the conventional Legendre transform cannot be applied to the first-degree Lagrangian ℒ(1)(xk,xk') to derive the relevant Hamiltonian ℋ(1)(xk,pk), and Hamiltonian ray equations. The reason is that the Hessian determinant of the transform vanishes identically for first-degree Lagrangians ℒ(1)(xk,xk'). 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 ℒ(2). Parameter u along the curves is taken to correspond to travel time τ, and the second-degree Lagrangian ℒ(2)(xk,{\dot x}k) is then introduced by the relation ℒ(2)(xk,{\dot x}k) = ½[ℒ(1)(xk,{\dot x}k)]2, with {\dot x}k = dxk÷dτ. The second-degree Lagrangian ℒ(2)(xk,{\dot x}k) yields the same Euler/Lagrange equations for rays as the first-degree Lagrangian ℒ(1)(xk,{\dot x}k). The relevant Hessian determinant, however, does not vanish identically. Consequently, the Legendre transform can then be used to compute Hamiltonian ℋ(2)(xk,pk) from Lagrangian ℒ(2)(xk,{\dot x}k), and vice versa, and the Hamiltonian canonical equations can be derived from the Euler-Lagrange equations. Both ℒ(2)(xk,{\dot x}k) and ℋ(2)(xk,pk) can be expressed in terms of the wave propagation metric tensor gij(xk,{\dot x}k), which depends not only on position xk, but also on the direction of vector {\dot x}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)(xk,{\dot x}k), not to the first-degree Lagrangian ℒ(1)(xk,{\dot x}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

Related content


Free Content
Free content
New Content
New content
Open Access Content
Open access content
Subscribed Content
Subscribed content
Free Trial Content
Free trial content

Text size:

A | A | A | A
Share this item with others: These icons link to social bookmarking sites where readers can share and discover new web pages. print icon Print this page