A model reduction approach to numerical inversion for a parabolic partial differential equation

Abstract

We propose a novel numerical inversion algorithm for the coefficients of parabolic partial differential equations, based on model reduction. The study is motivated by the application of controlled source electromagnetic exploration, where the unknown is the subsurface electrical resistivity and the data are time resolved surface measurements of the magnetic field. The algorithm presented in this paper considers inversion in one and two dimensions. The reduced model is obtained with rational interpolation in the frequency (Laplace) domain and a rational Krylov subspace projection method. It amounts to a nonlinear mapping from the function space of the unknown resistivity to the small dimensional space of the parameters of the reduced model. We use this mapping as a nonlinear preconditioner for the Gauss-Newton iterative solution of the inverse problem. The advantage of the inversion algorithm is twofold. First, the nonlinear preconditioner resolves most of the nonlinearity of the problem. Thus the iterations are less likely to get stuck in local minima and the convergence is fast. Second, the inversion is computationally efficient because it avoids repeated accurate simulations of the time-domain response. We study the stability of the inversion algorithm for various rational Krylov subspaces, and assess its performance with numerical experiments.

31 pages, 9 figures, 2 tables

References in corpus (1)

• Inverse boundary value problem for the Helmholtz equation: Multi-level approach and iterative reconstruction

Cited by in corpus (15)

• A convergent numerical method to recover the initial condition of nonlinear parabolic equations from lateral Cauchy data
• Direct, nonlinear inversion algorithm for hyperbolic problems via projection-based model reduction
• 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
• A new algorithm to determine the creation or depletion term of parabolic equations from boundary measurements
• Reduced order models for spectral domain inversion: Embedding into the continuous problem and generation of internal data
• On extension of the data driven ROM inverse scattering framework to partially nonreciprocal arrays
• Nonlinear seismic imaging via reduced order model backprojection
• A discrete Liouville identity for numerical reconstruction of Schrödinger potentials
• Reduced Order Modeling for First Order Hyperbolic Systems with Application to Multiparameter Acoustic Waveform Inversion
• Randomization for the Efficient Computation of Parametric Reduced Order Models for Inversion
• Robust Parameter Inversion using Adaptive Reduced Order Models
• Inverse scattering for Schrödinger equation in the frequency domain via data-driven reduced order modeling