×

Seismic imaging and optimal transport. (English) Zbl 1455.86009

Summary: Seismology has changed character since 50 years ago when the full wavefield could be determined. Partial differential equations (PDE) started to be used in the inverse process of finding properties of the interior of the earth. In this paper, we will review earlier techniques focusing on Full Waveform Inversion (FWI), which is a large-scale non-convex PDE constrained optimization problem. The minimization of the objective function is usually coupled with the adjoint state method, which also includes the solution to an adjoint wave equation. The least-squares \((L^2)\) norm is the conventional objective function measuring the difference between simulated and measured data, but it often results in the minimization trapped in local minima. One way to mitigate this is by selecting another misfit function with better convexity properties. Here we propose using the quadratic Wasserstein metric \((W_2)\) as a new misfit function in FWI. The optimal map defining \(W_2\) can be computed by solving a Monge-Ampère equation. Theorems pointing to the advantages of using optimal transport over \(L^2\) norm will be discussed, and several large-scale computational examples will be presented.

MSC:

86A15 Seismology (including tsunami modeling), earthquakes
86A22 Inverse problems in geophysics
94A08 Image processing (compression, reconstruction, etc.) in information and communication theory
65K10 Numerical optimization and variational techniques
49Q20 Variational problems in a geometric measure-theoretic setting
PDFBibTeX XMLCite
Full Text: DOI arXiv