Show simple item record

hal.structure.identifierInstitut des Sciences de la Terre [ISTerre]
dc.contributor.authorMétivier, Ludovic
HAL ID: 21579
hal.structure.identifierInstitut des Sciences de la Terre [ISTerre]
dc.contributor.authorBrossier, Romain
HAL ID: 21608
ORCID: 0000-0002-7195-8123
hal.structure.identifierCEntre de REcherches en MAthématiques de la DEcision [CEREMADE]
dc.contributor.authorMérigot, Quentin
HAL ID: 235
hal.structure.identifierLaboratoire Jean Kuntzmann [LJK]
dc.contributor.authorOudet, Edouard
hal.structure.identifierInstitut des Sciences de la Terre [ISTerre]
dc.contributor.authorVirieux, Jean
dc.subjectoptimal transporten
dc.subjectfull waveform inversionen
dc.subjectnon-smooth convex optimizationen
dc.titleAn optimal transport approach for seismic tomography: application to 3D full waveform inversionen
dc.typeArticle accepté pour publication ou publié
dc.description.abstractenThe use of optimal transport distance has recently yielded significant progress in image processing for pattern recognition, shape identification, and histograms matching. In this study, the use of this distance is investigated for a seismic tomography problem exploiting the complete waveform; the full waveform inversion. In its conventional formulation, this high resolution seismic imaging method is based on the minimization of the L2 distance between predicted and observed data. Application of this method is generally hampered by the local minima of the associated L2 misfit function, which correspond to velocity models matching the data up to one or several phase shifts. Conversely, the optimal transport distance appears as a more suitable tool to compare the misfit between oscillatory signals, for its ability to detect shifted patterns. However, its application to the full waveform inversion is not straightforward, as the mass conservation between the compared data cannot be guaranteed, a crucial assumption for optimal transport. In this study, the use of a distance based on the Kantorovich–Rubinstein norm is introduced to overcome this difficulty. Its mathematical link with the optimal transport distance is made clear. An efficient numerical strategy for its computation, based on a proximal splitting technique, is introduced. We demonstrate that each iteration of the corresponding algorithm requires solving the Poisson equation, for which fast solvers can be used, relying either on the fast Fourier transform or on multigrid techniques. The development of this numerical method make possible applications to industrial scale data, involving tenths of millions of discrete unknowns. The results we obtain on such large scale synthetic data illustrate the potentialities of the optimal transport for seismic imaging. Starting from crude initial velocity models, optimal transport based inversion yields significantly better velocity reconstructions than those based on the L2 distance, in 2D and 3D contexts.en
dc.relation.isversionofjnlnameInverse Problems
dc.relation.isversionofjnlpublisherIOP Scienceen

Files in this item


There are no files associated with this item.

This item appears in the following Collection(s)

Show simple item record