A fast marching algorithm for the factored eikonal equation
From MaRDI portal
Publication:525940
DOI10.1016/J.JCP.2016.08.012zbMATH Open1360.65266arXiv1607.00973OpenAlexW2474675250MaRDI QIDQ525940FDOQ525940
Publication date: 8 May 2017
Published in: Journal of Computational Physics (Search for Journal in Brave)
Abstract: The eikonal equation is instrumental in many applications in several fields ranging from computer vision to geoscience. This equation can be efficiently solved using the iterative Fast Sweeping (FS) methods and the direct Fast Marching (FM) methods. However, when used for a point source, the original eikonal equation is known to yield inaccurate numerical solutions, because of a singularity at the source. In this case, the factored eikonal equation is often preferred, and is known to yield a more accurate numerical solution. One application that requires the solution of the eikonal equation for point sources is travel time tomography. This inverse problem may be formulated using the eikonal equation as a forward problem. While this problem has been solved using FS in the past, the more recent choice for applying it involves FM methods because of the efficiency in which sensitivities can be obtained using them. However, while several FS methods are available for solving the factored equation, the FM method is available only for the original eikonal equation. In this paper we develop a Fast Marching algorithm for the factored eikonal equation, using both first and second order finite-difference schemes. Our algorithm follows the same lines as the original FM algorithm and requires the same computational effort. In addition, we show how to obtain sensitivities using this FM method and apply travel time tomography, formulated as an inverse factored eikonal equation. Numerical results in two and three dimensions show that our algorithm solves the factored eikonal equation efficiently, and demonstrate the achieved accuracy for computing the travel time. We also demonstrate a recovery of a 2D and 3D heterogeneous medium by travel time tomography using the eikonal equation for forward modelling and inversion by Gauss-Newton.
Full work available at URL: https://arxiv.org/abs/1607.00973
seismic imagingeikonal equationGauss-Newtonfast marchingfactored eikonal equationfirst arrivaltravel time tomography
Numerical methods for inverse problems for boundary value problems involving PDEs (65N21) Finite difference methods for boundary value problems involving PDEs (65N06) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12)
Cites Work
- jInv--a Flexible Julia Package for PDE Parameter Estimation
- Nonlinear total variation based noise removal algorithms
- Viscosity Solutions of Hamilton-Jacobi Equations
- Nonstationary inverse problems and state estimation
- Computational Methods for Inverse Problems
- Fast Sweeping Algorithms for a Class of Hamilton--Jacobi Equations
- A fast sweeping method for Eikonal equations
- Efficient algorithms for globally optimal trajectories
- A fast marching level set method for monotonically advancing fronts.
- Lax-Friedrichs sweeping scheme for static Hamilton-Jacobi equations
- A Viscosity Solutions Approach to Shape-From-Shading
- Fast Marching Methods
- Computing geodesic paths on manifolds
- Fast sweeping methods for factored anisotropic eikonal equations: multiplicative and additive factors
- A fast sweeping method for static convex Hamilton-Jacobi equations
- Fast Huygens sweeping methods for Helmholtz equations in inhomogeneous media in the high frequency regime
- High-Order Factorization Based High-Order Hybrid Fast Sweeping Methods for Point-Source Eikonal Equations
- Computational Study of Fast Methods for the Eikonal Equation
- Factored singularities and high-order Lax-Friedrichs sweeping schemes for point-source traveltimes and amplitudes
- Fast sweeping method for the factored eikonal equation
- An efficient solution to the eikonal equation on parametric manifolds
- An adjoint state method for three-dimensional transmission traveltime tomography using first-arrivals
- A fast method for the solution of the Helmholtz equation
Cited In (11)
- Single pass computation of first seismic wave travel time in three dimensional heterogeneous media with general anisotropy
- Corner cases, singularities, and dynamic factoring
- Full Waveform Inversion Guided by Travel Time Tomography
- Ordered line integral methods for solving the eikonal equation
- Hybrid fast sweeping methods for anisotropic eikonal equation in two-dimensional tilted transversely isotropic media
- jInv--a Flexible Julia Package for PDE Parameter Estimation
- Multi-stencils fast marching method for factored eikonal equations with quadratic anisotropy
- Approximate Inversion of a Class of Generalized Radon Transforms
- Neural eikonal solver: improving accuracy of physics-informed neural networks for solving eikonal equation in case of caustics
- Development of efficient and robust eikonal solver variants for first-arrival seismic modeling
- An Uncertainty-Weighted Asynchronous ADMM Method for Parallel PDE Parameter Estimation
Uses Software
This page was built for publication: A fast marching algorithm for the factored eikonal equation
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q525940)
