Full waveform inversion guided by travel time tomography
From MaRDI portal
(Redirected from Publication:5372644)
eikonal equationfast marchingfull waveform inversionGauss-NewtonHelmholtz equationinverse problemsparallel computationsshifted Laplacian multigridtravel time tomography
PDEs in connection with geophysics (35Q86) Inverse problems for PDEs (35R30) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Multigrid methods; domain decomposition for boundary value problems involving PDEs (65N55) Seismology (including tsunami modeling), earthquakes (86A15) Inverse problems in geophysics (86A22)
Abstract: Full waveform inversion (FWI) is a process in which seismic numerical simulations are fit to observed data by changing the wave velocity model of the medium under investigation. The problem is non-linear, and therefore optimization techniques have been used to find a reasonable solution to the problem. The main problem in fitting the data is the lack of low spatial frequencies. This deficiency often leads to a local minimum and to non-plausible solutions. In this work we explore how to obtain low frequency information for FWI. Our approach involves augmenting FWI with travel time tomography, which has low-frequency features. By jointly inverting these two problems we enrich FWI with information that can replace low frequency data. In addition, we use high order regularization, in a preliminary inversion stage, to prevent high frequency features from polluting our model in the initial stages of the reconstruction. This regularization also promotes the non-dominant low-frequency modes that exist in the FWI sensitivity. By applying a joint FWI and travel time inversion we are able to obtain a smooth model than can later be used to recover a good approximation for the true model. A second contribution of this paper involves the acceleration of the main computational bottleneck in FWI--the solution of the Helmholtz equation. We show that the solution time can be reduced by solving the equation for multiple right hand sides using block multigrid preconditioned Krylov methods.
Recommendations
Cites work
- 3D Frequency-Domain Seismic Inversion with Controlled Sloppiness
- A block version of BiCGSTAB for linear systems with multiple right-hand sides
- A fast marching algorithm for the factored eikonal equation
- A fast marching level set method for monotonically advancing fronts.
- A fast method for the solution of the Helmholtz equation
- A Flexible Inner-Outer Preconditioned GMRES Algorithm
- A fully asynchronous multifrontal solver using distributed dynamic scheduling
- A Multigrid Tutorial, Second Edition
- A new level-dependent coarse grid correction scheme for indefinite Helmholtz problems.
- A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion
- A Novel Multigrid Based Preconditioner For Heterogeneous Helmholtz Problems
- A parallel sweeping preconditioner for heterogeneous 3D Helmholtz equations
- A penalty method for PDE-constrained optimization in inverse problems
- An improved two-grid preconditioner for the solution of three-dimensional Helmholtz problems in heterogeneous media.
- An Iterative Method for Nonsymmetric Systems with Multiple Right-Hand Sides
- An optimal transport approach for seismic tomography: application to 3D full waveform inversion
- Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems
- Complex additive geometric multilevel solvers for Helmholtz equations on spacetrees
- Computational Methods for Inverse Problems
- Computational methods in geophysical electromagnetics
- Factored singularities and high-order Lax-Friedrichs sweeping schemes for point-source traveltimes and amplitudes
- Fast Marching Methods
- Fast sweeping method for the factored eikonal equation
- Fast sweeping methods for factored anisotropic eikonal equations: multiplicative and additive factors
- Flexible variants of block restarted GMRES methods with application to geophysics
- GPU implementation of a Helmholtz Krylov solver preconditioned by a shifted Laplace multigrid method
- High-order factorization based high-order hybrid fast sweeping methods for point-source eikonal equations
- scientific article; zbMATH DE number 45848 (Why is no real title available?)
- scientific article; zbMATH DE number 46855 (Why is no real title available?)
- scientific article; zbMATH DE number 1561761 (Why is no real title available?)
- scientific article; zbMATH DE number 2087097 (Why is no real title available?)
- jInv -- a flexible Julia package for PDE parameter estimation
- Non-Galerkin multigrid based on sparsified smoothed aggregation
- Nonlinear total variation based noise removal algorithms
- Numerical Optimization
- On Improving Linear Solver Performance: A Block Variant of GMRES
- On optimization techniques for solving nonlinear inverse problems
- Recursive Krylov-based multigrid cycles
- Shifted-Laplacian preconditioners for heterogeneous Helmholtz problems
- Statistical and computational inverse problems.
Cited in
(26)- Preconditioner for estimation of multipole sources via full waveform inversion
- Enabling numerically exact local solver for waveform inversion -- a low-rank approach
- Two-grid full-waveform inversion in the time domain
- Accurate and efficient velocity estimation using transmission matrix formalism based on the domain decomposition method
- Full Waveform Inversion Using Extended and Simultaneous Sources
- Accounting for the definition domain of the forward map in traveltime tomography—application to the inversion of prismatic reflections
- Total variation regularization for seismic waveform inversion using an adaptive primal dual hybrid gradient method
- Total variation regularization strategies in full-waveform inversion
- Wavefield reconstruction inversion: an example
- A source independent proximal gradient method to full waveform inversion with sparsity constraints in frequency domain
- An all-at-once approach to full waveform inversion in the viscoelastic regime
- On the time-domain full waveform inversion for time-dissipative and dispersive poroelastic media
- Multimodal 3D Shape Reconstruction under Calibration Uncertainty Using Parametric Level Set Methods
- Source-IndependentAmplitude-SemblanceFull-Waveform InversionusingaHybridTime-andFrequency-Domain Approach
- Concepts in the Direct Waveform Inversion (DWI) Using Explicit Time-Space Causality
- CMFWI: coupled multiscenario full waveform inversion
- jInv -- a flexible Julia package for PDE parameter estimation
- A hybrid shifted Laplacian multigrid and domain decomposition preconditioner for the elastic Helmholtz equations
- SMIwiz: an integrated toolbox for multidimensional seismic modelling and imaging
- Efficient 1.5D full waveform inversion in the Laplace-Fourier domain
- Isogeometric multi-resolution full waveform inversion based on the finite cell method
- LFA-tuned matrix-free multigrid method for the elastic Helmholtz equation
- Optimising seismic imaging design parameters via bilevel learning
- Model order reduction for seismic applications
- Dispersion-corrected compact finite difference discretizations for the elastic Helmholtz equation in two and three dimensions
- Application of optimal basis functions in full waveform inversion
This page was built for publication: Full waveform inversion guided by travel time tomography
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5372644)
