Analysis of optimal transport and related misfit functions in full-waveform inversion

ABSTRACTFull-waveform inversion has evolved into a powerful computational tool in seismic imaging. New misfit functions for matching simulated and measured data have recently been introduced to avoid the traditional lack of convergence due to cycle skipping. We have introduced the Wasserstein distance from optimal transport for computing the misfit, and several groups are currently further developing this technique. We evaluate three essential observations of this new metric with implication for future development. One is the discovery that trace-by-trace comparison with the quadratic Wasserstein metric works remarkably well together with the adjoint-state method. Another is the close connection between optimal transport-based misfits and integrated techniques with normalization as, for example, the normalized integration method. Finally, we study the convexity with respect to selected model parameters for different normalizations and remark on the effect of normalization on the convergence of the adjoint...

[1]  Ru-Shan Wu,et al.  Seismic envelope inversion: reduction of local minima and noise resistance , 2015 .

[2]  J. Virieux,et al.  An optimal transport approach for seismic tomography: application to 3D full waveform inversion , 2016 .

[3]  Biondo Biondi,et al.  Tomographic full-waveform inversion (TFWI) by combining FWI and wave-equation migration velocity analysis , 2013 .

[4]  R. Pratt,et al.  INVERSE THEORY APPLIED TO MULTI‐SOURCE CROSS‐HOLE TOMOGRAPHY.: PART 1: ACOUSTIC WAVE‐EQUATION METHOD1 , 1990 .

[5]  William W. Symes,et al.  Velocity inversion by differential semblance optimization , 1991 .

[6]  Yann Brenier,et al.  A computational fluid mechanics solution to the Monge-Kantorovich mass transfer problem , 2000, Numerische Mathematik.

[7]  R. Gerhard Pratt,et al.  Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies , 2004 .

[8]  J. Virieux,et al.  Measuring the misfit between seismograms using an optimal transport distance: application to full waveform inversion , 2016 .

[9]  Brittany D. Froese,et al.  Application of the Wasserstein metric to seismic signals , 2013, 1311.4581.

[10]  A. Tarantola,et al.  Generalized Nonlinear Inverse Problems Solved Using the Least Squares Criterion (Paper 1R1855) , 1982 .

[11]  S. Operto,et al.  Which data residual norm for robust elastic frequency-domain full waveform inversion? , 2010 .

[12]  Changsoo Shin,et al.  Waveform inversion using a back-propagation algorithm and a Huber function norm , 2009 .

[13]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[14]  C. Bunks,et al.  Multiscale seismic waveform inversion , 1995 .

[15]  Brittany D. Froese,et al.  Optimal Transport for Seismic Full Waveform Inversion , 2016, 1602.01540.

[16]  Arthur B. Weglein,et al.  Inverse scattering series and seismic exploration , 2003 .

[17]  Bjorn Engquist,et al.  Application of optimal transport and the quadratic Wasserstein metric to full-waveform inversion , 2016, 1612.05075.

[18]  Giuseppe Savaré,et al.  A new class of transport distances between measures , 2008, 0803.1235.

[19]  Y. Brenier Polar Factorization and Monotone Rearrangement of Vector-Valued Functions , 1991 .

[20]  Sergey Fomel,et al.  Building good starting models for full-waveform inversion using adaptive matching filtering misfit , 2016 .

[21]  William W. Symes,et al.  Migration velocity analysis and waveform inversion , 2008 .

[22]  P. Lailly,et al.  Pre-stack inversion of a 1-D medium , 1986, Proceedings of the IEEE.

[23]  Mike Warner,et al.  Adaptive waveform inversion: Theory , 2014 .

[24]  J. Tromp,et al.  Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements , 2011 .