Skip to Main content Skip to Navigation
Preprints, Working Papers, ...

Single pass eikonal solver in tilted transversely anisotropic media

Abstract : We present a numerical scheme to solve the eikonal equation in a Tilted Transversely Isotropic (TTI) medium. The solution to this equation corresponds to the first arrival time of seismic pressure waves in the high frequency asymptotic regime, whose propagation speed is neither isotropic nor elliptic. Instead, the speed profile is characterized by a fourth degree polynomial equation in a rotated frame, defined in terms of the Thomsen or Hooke elasticity coefficients of the geophysical medium. We show that TTI eikonal equations can be expressed as the maximum or minimum of a family of Riemannian eikonal equations, for which efficient discretizations are known. Based on this observation, we propose an original scheme that is causal, thus solvable in a single pass over the domain, and Eulerian, hence also mapping well to massively parallel architectures. Numerical experiments illustrate the method's accuracy, speed and robustness, on both a problem with analytical solution and a realistic synthetic instance, and compare a CPU with a GPU implementation, with the GPU being fifty times faster than the CPU implementation.
Document type :
Preprints, Working Papers, ...
Complete list of metadata

https://hal.archives-ouvertes.fr/hal-03655290
Contributor : François Desquilbet Connect in order to contact the contributor
Submitted on : Tuesday, May 3, 2022 - 11:06:41 AM
Last modification on : Thursday, May 5, 2022 - 3:51:41 PM

File

Manuscript.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : hal-03655290, version 2

Citation

François Desquilbet, Jean-Marie Mirebeau, Ludovic Métivier. Single pass eikonal solver in tilted transversely anisotropic media. 2022. ⟨hal-03655290v2⟩

Share

Metrics

Record views

59

Files downloads

15