The quadratic Wasserstein metric for earthquake location
From MaRDI portal
Publication:2311623
Abstract: In [Engquist et al., Commun. Math. Sci., 14(2016)], the Wasserstein metric was successfully introduced to the full waveform inversion. We apply this method to the earthquake location problem. For this problem, the seismic stations are far from each other. Thus, the trace by trace comparison [Yang et al., arXiv(2016)] is a natural way to compare the earthquake signals. Under this framework, we have derived a concise analytic expression of the Fr`echet gradient of the Wasserstein metric, which leads to a simple and efficient implementation for the adjoint method. We square and normalize the earthquake signals for comparison so that the convexity of the misfit function with respect to earthquake hypocenter and origin time can be observed numerically. To reduce the impact of noise, which can not offset each other after squaring the signals, a new control parameter is introduced. Finally, the LMF (Levenberg-Marquardt-Fletcher) method is applied to solve the resulted optimization problem. According to the numerical experiments, only a few iterations are required to converge to the real earthquake hypocenter and origin time. Even for data with noise, we can obtain reasonable and convergent numerical results.
Recommendations
- The Wasserstein-Fisher-Rao Metric for Waveform Based Earthquake Location
- Optimal transport for seismic full waveform inversion
- Misfit function for full waveform inversion based on the Wasserstein metric with dynamic formulation
- Application of the Wasserstein metric to seismic signals
- A new earthquake location method based on the waveform inversion
Cites work
- scientific article; zbMATH DE number 3134565 (Why is no real title available?)
- scientific article; zbMATH DE number 1243473 (Why is no real title available?)
- scientific article; zbMATH DE number 1909499 (Why is no real title available?)
- A method for the solution of certain non-linear problems in least squares
- A new earthquake location method based on the waveform inversion
- An Algorithm for Least-Squares Estimation of Nonlinear Parameters
- An optimal transport approach for seismic tomography: application to 3D full waveform inversion
- Application of the Wasserstein metric to seismic signals
- Computational high frequency wave propagation
- Gaussian beam methods for the Schrödinger equation in the semi-classical regime: Lagrangian and Eulerian formulations
- High Order Numerical Quadratures to One Dimensional Delta Function Integrals
- Optimal transport for applied mathematicians. Calculus of variations, PDEs, and modeling
- Optimal transport for seismic full waveform inversion
- Scaling algorithms for unbalanced optimal transport problems
- Stability of a 4th-order curvature condition arising in optimal transport theory
- The Eulerian Gaussian beam method for high frequency wave propagation in the reduced momentum space
Cited in
(13)- The Wasserstein-Fisher-Rao Metric for Waveform Based Earthquake Location
- Fast sinkhorn. II: Collinear triangular matrix and linear time accurate computation of optimal transport
- The Wasserstein metric matrix and its computational property
- Optimal transport natural gradient for statistical manifolds with continuous sample space
- The auxiliary function method for waveform based earthquake location
- Earthquake location methods
- The quadratic Wasserstein metric for inverse data matching
- Misfit function for full waveform inversion based on the Wasserstein metric with dynamic formulation
- Optimal transportation for electrical impedance tomography
- Greedy selection of optimal location of sensors for uncertainty reduction in seismic moment tensor inversion
- Application of the Wasserstein metric to seismic signals
- A new earthquake location method based on the waveform inversion
- On Quadratic Wasserstein Metric with Squaring Scaling for Seismic Velocity Inversion
This page was built for publication: The quadratic Wasserstein metric for earthquake location
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2311623)