6533b873fe1ef96bd12d575f

RESEARCH PRODUCT

Higher-order Hamilton–Jacobi perturbation theory for anisotropic heterogeneous media: dynamic ray tracing in Cartesian coordinates

Einar IversenJoonas IlmavirtaTeemu SaksalaBjørn UrsinMaarten V. De Hoop

subject

Wave propagationBody wavesta1171body waveswave propagation010502 geochemistry & geophysics01 natural sciencesHamilton–Jacobi equationlaw.inventionGeochemistry and PetrologylawaaltoliikeCartesian coordinate system0101 mathematicsAnisotropy0105 earth and related environmental sciencesPhysics010102 general mathematicsMathematical analysista111computational seismologynumerical modellingseismic anisotropyseismologiaRay tracing (physics)Geophysicsnumeerinen analyysiapproksimointidifferentiaaliyhtälötComputational seismologynumerical approximations and analysis

description

With a Hamilton–Jacobi equation in Cartesian coordinates as a starting point, it is common to use a system of ordinary differential equations describing the continuation of first-order derivatives of phase-space perturbations along a reference ray. Such derivatives can be exploited for calculating geometrical spreading on the reference ray and for establishing a framework for second-order extrapolation of traveltime to points outside the reference ray. The continuation of first-order derivatives of phase-space perturbations has historically been referred to as dynamic ray tracing. The reason for this is its importance in the process of calculating amplitudes along the reference ray. We extend the standard dynamic ray-tracing scheme to include higher-order derivatives of the phase-space perturbations. The main motivation is to extrapolate and interpolate amplitude and phase properties of high-frequency Green’s functions to nearby (paraxial) source and receiver locations. Principal amplitude coefficients, geometrical spreading factors, geometrical spreading matrices, ray propagator matrices, traveltimes, slowness vectors and curvature matrices are examples of quantities for which we enhance the computation potential. This, in turn, has immediate applications in modelling, mapping and imaging. Numerical tests for 3-D isotropic and anisotropic heterogeneous models yield clearly improved extrapolation results for the traveltime and geometrical spreading. One important conclusion is that the extrapolation function for the geometrical spreading must be at least third order to be appropriate at large distances away from the reference ray. © The Author(s) 2018. Published by Oxford University Press on behalf of The Royal Astronomical Society. This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model)

http://hdl.handle.net/11250/2642091