TY - JOUR

T1 - Application of optimal transport and the quadratic Wasserstein metric to full-waveform inversion

AU - Yang, Yunan

AU - Engquist, Björn

AU - Sun, Junzhe

AU - Hamfeldt, Brittany F.

N1 - Funding Information:
We thank S. Fomel, Z. Xue, and L. Qiu for very helpful discussions, and we thank the sponsors of the Texas Consortium for Computational Seismology for financial support. J. Sun was additionally supported by the Statoil Fellows Program at the University of Texas at Austin. B. Engquist was partially supported by NSF DMS-1620396. B. Froese was partially supported by NSF DMS-1619807. Y. Yang was partially supported by a grant from the Simons Foundation (#419126, B. Froese).
Publisher Copyright:
© 2018 Society of Exploration Geophysicists.

PY - 2018/1/1

Y1 - 2018/1/1

N2 - Conventional full-waveform inversion (FWI) using the least-squares norm as a misfit function is known to suffer from cycle-skipping issues that increase the risk of computing a local rather than the global minimum of the misfit. The quadratic Wasserstein metric has proven to have many ideal properties with regard to convexity and insensitivity to noise. When the observed and predicted seismic data are considered to be two density functions, the quadraticWasserstein metric corresponds to the optimal cost of rearranging one density into the other, in which the transportation cost is quadratic in distance. Unlike the least-squares norm, the quadratic Wasserstein metric measures not only amplitude differences but also global phase shifts, which helps to avoid cycle-skipping issues. We have developed a new way of using the quadratic Wasserstein metric trace by trace in FWI and compare it with the global quadratic Wasserstein metric via the solution of the Monge- Ampère equation. We incorporate the quadratic Wasserstein metric technique into the framework of the adjoint-state method and apply it to several 2D examples. With the corresponding adjoint source, the velocity model can be updated using a quasi- Newton method. Numerical results indicate the effectiveness of the quadratic Wasserstein metric in alleviating cycle-skipping issues and sensitivity to noise. The mathematical theory and numerical examples demonstrate that the quadratic Wasserstein metric is a good candidate for a misfit function in seismic inversion.

AB - Conventional full-waveform inversion (FWI) using the least-squares norm as a misfit function is known to suffer from cycle-skipping issues that increase the risk of computing a local rather than the global minimum of the misfit. The quadratic Wasserstein metric has proven to have many ideal properties with regard to convexity and insensitivity to noise. When the observed and predicted seismic data are considered to be two density functions, the quadraticWasserstein metric corresponds to the optimal cost of rearranging one density into the other, in which the transportation cost is quadratic in distance. Unlike the least-squares norm, the quadratic Wasserstein metric measures not only amplitude differences but also global phase shifts, which helps to avoid cycle-skipping issues. We have developed a new way of using the quadratic Wasserstein metric trace by trace in FWI and compare it with the global quadratic Wasserstein metric via the solution of the Monge- Ampère equation. We incorporate the quadratic Wasserstein metric technique into the framework of the adjoint-state method and apply it to several 2D examples. With the corresponding adjoint source, the velocity model can be updated using a quasi- Newton method. Numerical results indicate the effectiveness of the quadratic Wasserstein metric in alleviating cycle-skipping issues and sensitivity to noise. The mathematical theory and numerical examples demonstrate that the quadratic Wasserstein metric is a good candidate for a misfit function in seismic inversion.

UR - http://www.scopus.com/inward/record.url?scp=85039984762&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85039984762&partnerID=8YFLogxK

U2 - 10.1190/GEO2016-0663.1

DO - 10.1190/GEO2016-0663.1

M3 - Article

AN - SCOPUS:85039984762

VL - 83

SP - R43-R62

JO - Geophysics

JF - Geophysics

SN - 0016-8033

IS - 1

ER -