Direct, Nonlinear Inversion Algorithm for Hyperbolic Problems via Projection-Based Model Reduction

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.

[1]  Liliana Borcea,et al.  Pyramidal resistor networks for electrical impedance tomography with partial boundary measurements , 2010 .

[2]  V. Druskin,et al.  A model reduction approach to numerical inversion for a parabolic partial differential equation , 2012, 1210.1257.

[3]  Liliana Borcea,et al.  Electrical impedance tomography with resistor networks , 2008 .

[4]  Evert Slob,et al.  Three-dimensional single-sided Marchenko inverse scattering, data-driven focusing, Green's function retrieval, and their mutual relations. , 2013, Physical review letters.

[5]  Vladimir Druskin,et al.  Gaussian Spectral Rules for the Three-Point Second Differences: I. A Two-Point Positive Definite Problem in a Semi-Infinite Domain , 1999, SIAM J. Numer. Anal..

[6]  M. Sondhi,et al.  Determination of vocal-tract shape from impulse response at the lips. , 1971, The Journal of the Acoustical Society of America.

[7]  W. Gautschi.  Orthogonal Polynomials: Computation and Approximation , 2004 .

[8]  B. Parlett.  The Symmetric Eigenvalue Problem , 1981 .

[9]  F. Santosa.  Numerical scheme for the inversion of acoustical impedance profile based on the Gelfand-Levitan method , 1982 .

[10]  Eric James Grimme,et al.  Krylov Projection Methods for Model Reduction , 1997 .

[11]  Liliana Borcea,et al.  INVERSE PROBLEMS PII: S0266-5611(02)33630-X Optimal finite difference grids for direct and inverse Sturm–Liouville problems , 2002 .

[12]  W. Symes,et al.  Inverse boundary value problems and a theorem of Gel'fand and Levitan , 1979 .

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

[14]  B. M. Levitan,et al.  Determination of a Differential Equation by Two of its Spectra , 1964 .

[15]  Liliana Borcea,et al.  On the continuum limit of a discrete inverse spectral problem on optimal finite difference grids , 2005 .

[16]  Robert Burridge,et al.  The Gelfand-Levitan, the Marchenko, and the Gopinath-Sondhi integral equations of inverse scattering theory, regarded in the context of inverse impulse-response problems , 1980 .

[17]  Vladimir Druskin,et al.  On combining model reduction and Gauss–Newton algorithms for inverse partial differential equation problems , 2007 .

[18]  Gene H. Golub,et al.  Calculation of Gauss quadrature rules , 1967, Milestones in Matrix Computation.

[19]  M. M. Sondhi,et al.  Inversion of the telegraph equation and the synthesis of nonuniform lines , 1971 .

[20]  I. S. Kats.  SPECTRAL FUNCTIONS OF A STRING , 1983 .

[21]  Abdigany D. Satybaev,et al.  Direct Methods of Solving Multidimensional Inverse Hyperbolic Problems , 2004 .

[22]  Vladimir Druskin,et al.  Optimal finite difference grids and rational approximations of the square root I. Elliptic problems , 2000 .

[23]  Paul Van Dooren,et al.  A rational Lanczos algorithm for model reduction , 1996, Numerical Algorithms.

[24]  T. Habashy,et al.  A generalized Gel'fand-Levitan-Marchenko integral equation , 1991 .

[25]  ERIC DE STURLER,et al.  Nonlinear Parametric Inversion Using Interpolatory Model Reduction , 2013, SIAM J. Sci. Comput..

[26]  T. Stieltjes.  Recherches sur les fractions continues , 1995 .

[27]  G. Golub,et al.  A survey of matrix inverse eigenvalue problems , 1986 .

[28]  Athanasios C. Antoulas,et al.  Approximation of Large-Scale Dynamical Systems , 2005, Advances in Design and Control.

[29]  Z. Strakos,et al.  Krylov Subspace Methods: Principles and Analysis , 2012 .

[30]  Liliana Borcea,et al.  Resistor network approaches to electrical impedance tomography , 2011 .

[31]  A. Fichtner.  Full Seismic Waveform Modelling and Inversion , 2011 .

[32]  Vladimir Druskin,et al.  Application of the Difference Gaussian Rules to Solution of Hyperbolic Problems , 2000 .

[33]  Paul Sava,et al.  Riemannian wavefield extrapolation , 2005 .

[34]  I. Gel'fand,et al.  On the determination of a differential equation from its spectral function , 1955 .

[35]  Robert Skelton,et al.  Model reductions using a projection formulation , 1987, 26th IEEE Conference on Decision and Control.

[36]  D. Sorensen,et al.  Approximation of large-scale dynamical systems: an overview , 2004 .

[37]  The Gel’fand-Levitan Method in the Inverse Scattering Problem , 1974 .

[38]  Kenneth P. Bube,et al.  The One-Dimensional Inverse Problem of Reflection Seismology , 1983 .