Direct, nonlinear inversion algorithm for hyperbolic problems via projection-based model reduction

Abstract

We estimate the wave speed in the acoustic wave equation from boundary measurements by constructing a reduced-order model (ROM) matching discrete time-domain data. The state-variable representation of the ROM can be equivalently viewed as a Galerkin projection onto the Krylov subspace spanned by the snapshots of the time-domain solution. The success of our algorithm hinges on the data-driven Gram--Schmidt orthogonalization of the snapshots that suppresses multiple reflections and can be viewed as a discrete form of the Marchenko--Gel'fand--Levitan--Krein algorithm. In particular, the orthogonalized snapshots are localized functions, the (squared) norms of which are essentially weighted averages of the wave speed. The centers of mass of the squared orthogonalized snapshots provide us with the grid on which we reconstruct the velocity. This grid is weakly dependent on the wave speed in traveltime coordinates, so the grid points may be approximated by the centers of mass of the analogous set of squared orthogonalized snapshots generated by a known reference velocity. We present results of inversion experiments for one- and two-dimensional synthetic models.

54 pages, 6 figures fixed typos and small errors expanded several sections to aid in understanding

References in corpus (1)

• Nonlinear seismic imaging via reduced order model backprojection

Cited by in corpus (11)

• Untangling the nonlinearity in inverse scattering with data-driven reduced order models
• A nonlinear method for imaging with acoustic waves via reduced order model backprojection
• Lippmann-Schwinger-Lanczos algorithm for inverse scattering problems
• Robust nonlinear processing of active array data in inverse scattering via truncated reduced order models
• Reduced order models for spectral domain inversion: Embedding into the continuous problem and generation of internal data
• Reduced order model approach for imaging with waves
• On extension of the data driven ROM inverse scattering framework to partially nonreciprocal arrays
• Multi-scale S-fraction reduced-order models for massive wavefield simulations
• Nonlinear seismic imaging via reduced order model backprojection
• Reduced Order Modeling for First Order Hyperbolic Systems with Application to Multiparameter Acoustic Waveform Inversion
• Inverse scattering for Schrödinger equation in the frequency domain via data-driven reduced order modeling