A numerical study of multi-parameter full waveform inversion with iterative regularization using multi-frequency vibroseis data

Jia Shi, Elena Beretta, Maarten V. de Hoop, Elisa Francini, Sergio Vessella

Research output: Contribution to journalArticle

Abstract

We study the inverse boundary value problem for time-harmonic elastic waves, for the recovery of P- and S-wave speeds from vibroseis data or the Neumann-to-Dirichlet map. Our study is based on our recent result pertaining to the uniqueness and a conditional Lipschitz stability estimate for parametrizations on unstructured tetrahedral meshes of this inverse boundary value problem. With the conditional Lipschitz stability estimate, we design a procedure for full waveform inversion (FWI) with iterative regularization. The iterative regularization is implemented by projecting gradients, after scaling, onto subspaces associated with the mentioned parametrizations yielding Lipschitz stability. The procedure is illustrated in computational experiments using the continuous Galerkin finite element method of recovering the rough shapes and wave speeds of geological bodies from simple starting models, near and far from the boundary, that is, the free surface.

Original languageEnglish (US)
Pages (from-to)89-107
Number of pages19
JournalComputational Geosciences
Volume24
Issue number1
DOIs
StatePublished - Feb 1 2020

Keywords

  • Finite element method
  • Full waveform inversion
  • Stability and convergence

ASJC Scopus subject areas

  • Computer Science Applications
  • Computers in Earth Sciences
  • Computational Theory and Mathematics
  • Computational Mathematics

Fingerprint Dive into the research topics of 'A numerical study of multi-parameter full waveform inversion with iterative regularization using multi-frequency vibroseis data'. Together they form a unique fingerprint.

  • Cite this