This study is a review of recent applications to seismic imaging of optimal transport based numerical tools. Modern seismic imaging methods used in the industry rely on the interpretation of the full signal. The characterization of the subsurface mechanical properties is formulated as a PDE-constrained optimization problem, solved through local optimization strategies. The choice of the misfit function used to measure the distance between actual seismic recordings and those synthetized by the solution of wave propagation PDE is crucial. Indeed, the conventional least-squares distance function leads to a non-convex optimization problem whose solution through local optimization then strongly depends on the initial guess. Using an optimal transport distance is an interesting alternative from its convexity properties with respect to translation and dilation. Specific strategies need however to be implemented as seismic data are oscillatory, while the optimal transport theory has been developed for the comparison of positive measures. In this study we review two optimal transport based misfit functions, from their mathematical formulation to their application to field data through their numerical implementation. Advantages and drawbacks of both strategies are discussed. Numerical experiments show that they represent two interesting and complementary alternative to the classical least-squares misfit function, mitigating the dependency to the choice of the initial guess.

Keywords: Optimal transport, convexity, optimization, seismic imaging

^{1}; Romain Brossier

^{2}; Félix Kpadonou

^{3}; Jérémie Messud

^{3}; Arnaud Pladys

^{2}

@article{MSIA_2022__11_1_3_0, author = {Ludovic M\'etivier and Romain Brossier and F\'elix Kpadonou and J\'er\'emie Messud and Arnaud Pladys}, title = {A review of the use of optimal transport distances for high resolution seismic imaging based on the full waveform}, journal = {MathematicS In Action}, pages = {3--42}, publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles}, volume = {11}, number = {1}, year = {2022}, doi = {10.5802/msia.15}, language = {en}, url = {https://msia.centre-mersenne.org/articles/10.5802/msia.15/} }

TY - JOUR AU - Ludovic Métivier AU - Romain Brossier AU - Félix Kpadonou AU - Jérémie Messud AU - Arnaud Pladys TI - A review of the use of optimal transport distances for high resolution seismic imaging based on the full waveform JO - MathematicS In Action PY - 2022 SP - 3 EP - 42 VL - 11 IS - 1 PB - Société de Mathématiques Appliquées et Industrielles UR - https://msia.centre-mersenne.org/articles/10.5802/msia.15/ DO - 10.5802/msia.15 LA - en ID - MSIA_2022__11_1_3_0 ER -

%0 Journal Article %A Ludovic Métivier %A Romain Brossier %A Félix Kpadonou %A Jérémie Messud %A Arnaud Pladys %T A review of the use of optimal transport distances for high resolution seismic imaging based on the full waveform %J MathematicS In Action %D 2022 %P 3-42 %V 11 %N 1 %I Société de Mathématiques Appliquées et Industrielles %U https://msia.centre-mersenne.org/articles/10.5802/msia.15/ %R 10.5802/msia.15 %G en %F MSIA_2022__11_1_3_0

Ludovic Métivier; Romain Brossier; Félix Kpadonou; Jérémie Messud; Arnaud Pladys. A review of the use of optimal transport distances for high resolution seismic imaging based on the full waveform. MathematicS In Action, Volume 11 (2022) no. 1, pp. 3-42. doi : 10.5802/msia.15. https://msia.centre-mersenne.org/articles/10.5802/msia.15/

[1] Improving full-waveform inversion by wavefield reconstruction with alternating direction method of multipliers, Geophysics, Volume 84 (2019) no. 1, p. R139-R162 | DOI

[2] A genuinely polynomial primal simplex algorithm for the assignment problem, Discrete Appl. Math., Volume 45 (1993) no. 2, pp. 93-115 | DOI | MR | Zbl

[3] Quantitative Seismology: Theory and Methods, W. H. Freeman and Co, 1980

[4] Lecture Notes on Optimal Transport Problems, Mathematical Aspects of Evolving Interfaces (Lecture Notes in Mathematics), Volume 1812, Springer, 2003, pp. 1-52 | DOI | MR | Zbl

[5] Gradient flow of the Chapman Rubinstein Schatzman model for signed vortices, Ann. Inst. Henri Poincaré, Anal. Non Linéaire, Volume 28 (2011) no. 2, pp. 217-246 | Numdam | MR | Zbl

[6] Signature Methods for the Assignment Problem, Oper. Res., Volume 33 (1985) no. 3, pp. 527-536 | DOI | MR | Zbl

[7] Life of Field Seismic Implementation - Another “First at Valhal”, 65th EAGE Workshop, Stavanger (2003)

[8] Network Optimization: Continuous and Discrete Models, Athena Scientific, 1998

[9] The auction algorithm for the transportation problem, Ann. Oper. Res., Volume 20 (1989) no. 1, pp. 67-96 | DOI | MR | Zbl

[10] On the imaging of reflectors in the Earth, Geophysics, Volume 52 (1987) no. 7, pp. 931-942 | DOI

[11] Measure Theory. Vol. I and II, Springer, 2007 | DOI

[12] Global adjoint tomography: first-generation model, Geophys. J. Int., Volume 207 (2016) no. 3, pp. 1739-1766 | DOI

[13] Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements, Geophys. J. Int., Volume 185 (2011) no. 2, pp. 845-870 | DOI

[14] Multi-level adaptive solutions to boundary-value problems, Math. Comput., Volume 31 (1977), pp. 333-390 | DOI | MR | Zbl

[15] 2D elastic full waveform imaging of the near surface: Application to synthetic and a physical modelling data sets, Near Surface Geophysics, Volume 11 (2013), pp. 307-316 | DOI

[16] Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion, Geophysics, Volume 74 (2009) no. 6, p. WCC105-WCC118 | DOI

[17] Velocity model building from seismic reflection data by full waveform inversion, Geophysical Prospecting, Volume 63 (2015), pp. 354-367 | DOI

[18] Multiscale seismic waveform inversion, Geophysics, Volume 60 (1995) no. 5, pp. 1457-1473 | DOI

[19] Assignment Problems, Society for Industrial and Applied Mathematics, 2012 | DOI

[20] Optimal Transport Full-Waveform Inversion - Applications, 82nd EAGE Conference and Exhibition, Expanded Abstracts, European Association of Geoscientists and Engineers, 2020, p. Th Dome1 17

[21] Imaging the Earth’s interior, Blackwell Scientific Publication, 1985

[22] Proximal Splitting Methods in Signal Processing, Fixed-Point Algorithms for Inverse Problems in Science and Engineering (H. H. Bauschke; R. S. Burachik; L. ttes P. Combe; V. Elser; D. R. Luke; H. Wolkowicz, eds.) (Springer Optimization and Its Applications), Volume 49, Springer, 2011, pp. 185-212 | DOI | MR | Zbl

[23] Movie and video scale-time equalization application to flicker reduction, IEEE Trans. Image Process., Volume 15 (2006) no. 1, pp. 241-248 | DOI

[24] Texture Mapping via Optimal Mass Transport, IEEE Trans. Visual. Comput. Graph., Volume 16 (2010) no. 3, pp. 419-433 | DOI

[25] Application of the Wasserstein metric to seismic signals, Commun. Math. Sci., Volume 12 (2014) no. 5, pp. 979-988 | DOI | MR | Zbl

[26] Optimal transport for seismic full waveform inversion, Commun. Math. Sci., Volume 14 (2016) no. 8, pp. 2309-2330 | DOI | MR | Zbl

[27] Theoretical background for continental- and global-scale full-waveform inversion in the time-frequency domain, Geophys. J. Int., Volume 175 (2008), pp. 665-685 | DOI

[28] Full waveform tomography for radially anisotropic structure: New insights into present and past states of the Australasian upper mantle, Earth and Planetary Science Letters, Volume 290 (2010) no. 3-4, pp. 270-280 | DOI

[29] Two-dimensional nonlinear inversion of seismic waveforms: numerical results, Geophysics, Volume 51 (1986) no. 7, pp. 1387-1403 | DOI

[30] The role of attenuation in 2D full-waveform inversion of shallow-seismic body and Rayleigh waves, Geophysics, Volume 79 (2014) no. 6, p. R247-R261 | DOI

[31] Imaging Complex Fault Structures On-shore Oman Using Optimal Transport Full-Waveform Inversion, 82nd EAGE Conference and Exhibition, Expanded Abstracts, European Association of Geoscientists and Engineers, 2020, p. We Dome1 19

[32] Broadband FAZ land data: an opportunity for FWI, 81st EAGE Conference and Exhibition, Workshop Programme, European Association of Geoscientists and Engineers, 2019, p. WS08 11

[33] Waveform inversion via source extension (2019), pp. 4761-4766 | DOI

[34] Towards 3D 9C Elastic Full Waveform Inversion of Shallow Seismic Wavefields - Case Study Ettlingen Line, 81st EAGE Conference and Exhibition, Expanded Abstracts, European Association of Geoscientists and Engineers, 2019, p. We P01 04 | DOI

[35] Wavelengths of Earth structures that can be resolved from seismic reflection data, Geophysics, Volume 54 (1989) no. 7, pp. 906-910 | DOI

[36] Multiparameter full-waveform inversion of 3D ocean-bottom cable data from the Valhall field, Geophysics, Volume 86 (2021) no. 1, p. B15-B35 | DOI

[37] On the transfer of masses, Dokl. Akad. Nauk SSSR, Volume 37 (1942), pp. 7-8

[38] Optimal transport FWI with graph transform: Analysis and proposal of a partial shift strategy, 82nd EAGE Annual Conference and Exhibition, European Association of Geoscientists & Engineers, 2021 | DOI

[39] The Hungarian method for the assignment problem, Nav. Res. Logist. Q., Volume 2 (1955) no. 1-2, pp. 83-97 | DOI | MR | Zbl

[40] The seismic inverse problem as a sequence of before stack migrations, Conference on Inverse Scattering, Theory and application, Society for Industrial and Applied Mathematics, Philadelphia (1983), pp. 206-220 | MR

[41] Stereotomography, Geophysics, Volume 73 (2008) no. 5, p. VE25-VE34 | DOI

[42] Imaging with Kantorovich–Rubinstein Discrepancy, SIAM J. Imaging Sci., Volume 7 (2014) no. 4, pp. 2833-2859 | DOI | MR | Zbl

[43] A deconvolution-based objective function for wave-equation inversion, SEG Technical Program Expanded Abstracts 2011, Society of Exploration Geophysicists, 2011 | DOI

[44] Wave-equation traveltime inversion, Geophysics, Volume 56 (1991) no. 5, pp. 645-653 | DOI

[45] A description of transport cost for signed measures, J. Math. Sci., New York, Volume 181 (2012) no. 6, pp. 837-855 | DOI | MR | Zbl

[46] Marmousi2: An elastic upgrade for Marmousi, The Leading Edge, Volume 25 (2006) no. 2, pp. 156-166 | DOI

[47] Optimal transport in full-waveform inversion: Analysis and practice of the multidimensional Kantorovich-Rubinstein norm, Inverse Probl., Volume 37 (2021) no. 065012, pp. 1-42 | MR | Zbl

[48] Multidimensional Optimal Transport for 3D FWI: Demonstration on Field Data, 81st EAGE Conference and Exhibition, Expanded Abstracts, European Association of Geoscientists and Engineers, 2019, p. Tu R08 02

[49] On the Use of Optimal Transport Distances for a PDE-Constrained Optimization Problem in Seismic Imaging, Frontiers in PDE-Constrained Optimization (Harbir Antil; Drew P. Kouri; Martin-D. Lacasse; Denis Ridzal, eds.), Springer, 2018, pp. 377-397 | DOI | Zbl

[50] Optimal transport for mitigating cycle skipping in full waveform inversion: a graph space transform approach, Geophysics, Volume 83 (2018) no. 5, p. R515-R540 | DOI

[51] New insights on the graph space optimal transport distance for full waveform inversion, SEG Technical Program Expanded Abstracts 2021 (2021)

[52] Receiver-extension strategy for time-domain full-waveform inversion using a relocalization approach, Geophysics, Volume 87 (2022) no. 1, p. R13-R33 | DOI

[53] A graph space optimal transport distance as a generalization of ${L}^{p}$ distances: application to a seismic imaging inverse problem, Inverse Probl., Volume 35 (2019) no. 8, 085001, 49 pages | DOI | MR | Zbl

[54] Increasing the robustness and applicability of full waveform inversion: an optimal transport distance strategy, The Leading Edge, Volume 35 (2016) no. 12, pp. 1060-1067 | DOI | Zbl

[55] Measuring the misfit between seismograms using an optimal transport distance: Application to full waveform inversion, Geophys. J. Int., Volume 205 (2016), pp. 345-377 | DOI

[56] An optimal transport approach for seismic tomography: Application to 3D full waveform inversion, Inverse Probl., Volume 32 (2016) no. 11, p. 115008 | DOI | MR | Zbl

[57] Mémoire sur la théorie des déblais et des remblais, Histoire de l’Académie Royale des Sciences de Paris (1781)

[58] Updating Quasi-Newton Matrices With Limited Storage, Math. Comput., Volume 35 (1980) no. 151, pp. 773-782 | DOI | MR | Zbl

[59] Numerical Optimization, Springer, 2006

[60] A Breviary of Seismic Tomography, Cambridge University Press, 2008 | DOI

[61] A guided tour of multiparameter full waveform inversion for multicomponent data: from theory to practice, The Leading Edge, Volume 32 (2013) no. 9, pp. 1040-1054 (Special section Full Waveform Inversion) | DOI

[62] Efficient 3-D frequency-domain mono-parameter full-waveform inversion of ocean-bottom cable data: application to Valhall in the visco-acoustic vertical transverse isotropic approximation, Geophys. J. Int., Volume 202 (2015) no. 2, pp. 1362-1391 | DOI

[63] Automated colour grading using colour distribution transfer, Computer Vision and Image Understanding, Volume 107 (2007) no. 1, pp. 123-137 (Special issue on color image processing) | DOI

[64] Robust FWI with graph space optimal transport: application to 3D OBC Valhall data, Geophysics, Volume 87 (2022) no. 3, pp. 1-76

[65] On cycle-skipping and misfit function modification for full-wave inversion: Comparison of five recent approaches, Geophysics, Volume 86 (2021) no. 4, p. R563-R587 | DOI

[66] A review of the adjoint-state method for computing the gradient of a functional with geophysical applications, Geophys. J. Int., Volume 167 (2006) no. 2, pp. 495-503 | DOI

[67] Full waveform inversion of a deep water ocean bottom seismometer dataset, First Break, Volume 28 (2010), pp. 71-78

[68] FWI with Optimal Transport: A 3D Implementation and an Application on a Field Dataset, 80th EAGE Conference and Exhibition, European Association of Geoscientists and Engineers, 2018, p. We A12 02

[69] On the equality between Monge’s infimum and Kantorovich’s minimum in optimal mass transportation, Ann. Inst. Henri Poincaré, Probab. Stat., Volume 43 (2007) no. 1, pp. 1-13 | DOI | Numdam | MR | Zbl

[70] On the footprint of anisotropy on isotropic full waveform inversion: the Valhall case study, Geophys. J. Int., Volume 187 (2011), pp. 1495-1515 | DOI

[71] Multiparameter full waveform inversion of multicomponent OBC data from Valhall. Part 1: imaging compressional wavespeed, density and attenuation, Geophys. J. Int., Volume 194 (2013) no. 3, pp. 1640-1664 | DOI

[72] Joint FWI of diving and reflected waves using a graph space optimal transport distance: Synthetic tests on limited-offset surface seismic data (2020), pp. 780-784 | DOI

[73] Full-waveform inversion with an exponentially encoded optimal-transport norm, SEG Technical Program Expanded Abstracts 2017 (2017), pp. 1286-1290 | DOI

[74] Geodesic Shape Retrieval via Optimal Mass Transport, Computer Vision – ECCV 2010 (2010), pp. 771-784 | DOI

[75] Wasserstein Barycenter and Its Application to Texture Mixing, Scale Space and Variational Methods in Computer Vision (2012), pp. 435-446 | DOI

[76] The Earth Mover’s Distance as a Metric for Image Retrieval, Int. J. Comput. Vision, Volume 40 (2000) no. 2, pp. 99-121 | DOI | Zbl

[77] Optimal Transport for Applied Mathematicians: Calculus of Variations, PDEs, and Modeling, Progress in Nonlinear Differential Equations and their Applications, Springer, 2015 | DOI

[78] 2D Full Waveform Inversion of Recorded Shallow Seismic Rayleigh Waves on a Significantly 2D Structure, Proceedings of 19th European Meeting of Environmental and Engineering Geophysics, Expanded Abstracts, Bochum, Germany (2013)

[79] Acoustic Land Full-Waveform Inversion on a Broadband Land Dataset: The Impact of Optimal Transport, 81st EAGE Conference and Exhibition, Expanded Abstracts, European Association of Geoscientists and Engineers, 2019, p. Th R08 07

[80] Two-dimensional full wavefield inversion of wide-aperture marine seismic streamer data, Geophys. J. Int., Volume 151 (2002), pp. 325-344 | DOI

[81] Full waveform inversion: the next leap forward in imaging at Valhall, First Break, Volume 28 (2010), pp. 65-70 | DOI

[82] Multiparameter waveform inversion of a large wide-azimuth low-frequency land data set in Oman, Geophysics, Volume 79 (2014) no. 3, p. WA69-WA77 | DOI

[83] A Direct Method for the Discrete Solution of Separable Elliptic Equations, SIAM J. Numer. Anal., Volume 11 (1974) no. 6, pp. 1136-1150 | DOI | MR | Zbl

[84] Migration velocity analysis and waveform inversion, Geophysical Prospecting, Volume 56 (2008), pp. 765-790 | DOI

[85] Seismic tomography of the southern California crust based on spectral-element and adjoint methods, Geophys. J. Int., Volume 180 (2010), pp. 433-462 | DOI

[86] Inversion of seismic reflection data in the acoustic approximation, Geophysics, Volume 49 (1984) no. 8, pp. 1259-1266 | DOI

[87] Mitigating local minima in full-waveform inversion by expanding the search space, Geophys. J. Int., Volume 195 (2013) no. 1, pp. 661-667 | DOI

[88] A penalty method for PDE-constrained optimization in inverse problems, Inverse Probl., Volume 32 (2016) no. 1, 015007, 26 pages | MR | Zbl

[89] A correlation-based misfit criterion for wave-equation traveltime tomography, Geophys. J. Int., Volume 182 (2010) no. 3, pp. 1383-1394 | DOI

[90] Topics in optimal transportation, Graduate Studies in Mathematics, 50, American Mathematical Society, 2003

[91] Optimal transport: old and new, Grundlehren der Mathematischen Wissenschaften, Springer, 2008

[92] An introduction to Full Waveform Inversion, Encyclopedia of Exploration Geophysics, Society of Exploration Geophysicists, 2017, p. R1-1–R1-40 | DOI

[93] Reflection seismic waveform tomography, J. Geophys. Res., Volume 114 (2009) no. B3, pp. 1978-2012

[94] Adaptive waveform inversion: Theory, Geophysics, Volume 81 (2016) no. 6, p. R429-R445 | DOI

[95] Seismic envelope inversion and modulation signal model, Geophysics, Volume 79 (2014) no. 3, p. WA13-WA24

[96] Inversion on Reflected Seismic Wave, 2012 | DOI

[97] A Time-Domain Preconditioned Truncated Newton Approach to Multiparameter Visco-acoustic Full Waveform Inversion, SIAM J. Sci. Comput., Volume 40 (2018) no. 4, p. B1101-B1130 | DOI | Zbl

[98] Analysis of optimal transport and related misfit functions in full-waveform inversion, Geophysics, Volume 83 (2018) no. 1, p. A7-A12 | DOI

[99] Model recovery below reflectors by optimal-transport FWI (2018), pp. 1178-1182 | DOI

[100] Improving optimal transport based FWI through data normalization (2019), pp. 1315-1319 | DOI

[101] Application of optimal transport and the quadratic Wasserstein metric to full-waveform inversion, Geophysics, Volume 83 (2018) no. 1, p. R43-R62 | DOI

[102] Full Waveform Inversion of Diving and Reflected Waves for Velocity Model Building with Impedance Inversion Based on Scale Separation, Geophys. J. Int., Volume 202 (2015) no. 3, pp. 1535-1554 | DOI

*Cited by Sources: *