Solution of the time-domain inverse resistivity problem in the model reduction framework part I. one-dimensional problem with SISO data

摘要

Many time-domain problems in engineering applications can be described by means of parameter-dependent time-invariant dynamic systems. We are interested in parameter estimation, by fitting available transient measurements using the nonlinear least square method. As the main application, we consider the control source electromagnetic (CSEM) method of geophysical exploration governed by the diffusion Maxwell system, where the unknown parameters describe the spatial distribution of electrical resistivity. We propose a novel model reduction approach for constructing an efficient approximation of the Jacobian. The reduction is based on subspace projection of the state variable, and it allows us to split the time and space dependence of the derivative. We examine several popular rational Krylov subspaces and single out the H2-optimal subspace that not only minimizes the approximation error but completely annuls its influence on the inversion result. Preliminary numerical experiments with a simplified one-dimensional, single-input/single-output CSEM setting are reported to validate our strategy.