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

被引:23
|
作者
Borcea, Liliana [1 ]
Druskin, Vladimir [2 ]
Mamonov, Alexander V. [3 ]
Zaslavsky, Mikhail [2 ]
机构
[1] Univ Michigan, Dept Math, Ann Arbor, MI 48109 USA
[2] Schlumberger Doll Res Ctr, Cambridge, MA 02139 USA
[3] Schlumberger, Houston, TX 77042 USA
关键词
inverse problem; parabolic equation; model reduction; rational Krylov subspace projection; CSEM; GAUSSIAN SPECTRAL RULES; DOMAIN; ALGORITHM; FRAMEWORK; GRIDS;
D O I
10.1088/0266-5611/30/12/125011
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
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.
引用
收藏
页数:33
相关论文
共 50 条
  • [41] Numerical solution and distinguishability in time fractional parabolic equation
    Ali Demir
    Fatma Kanca
    Ebru Ozbilge
    Boundary Value Problems, 2015
  • [42] NUMERICAL CRITICAL VALUE OF THE GENERALIZED PARABOLIC MEMS EQUATION
    Nabongo, Diabate
    JOURNAL OF MATHEMATICAL ANALYSIS, 2016, 7 (04): : 45 - 63
  • [43] On MIMO model reduction by the weighted equation-error approach
    Krajewski, Wieslaw
    Viaro, Umberto
    NUMERICAL ALGORITHMS, 2007, 44 (01) : 83 - 98
  • [44] Determining a potential of the parabolic equation from partial boundary measurements
    Fan, Jiaming
    Duan, Zhiwen
    INVERSE PROBLEMS, 2021, 37 (09)
  • [45] The partial boundary value conditions of nonlinear degenerate parabolic equation
    Yuan Zhi
    Huashui Zhan
    Boundary Value Problems, 2022
  • [46] A Benchmark for the Bayesian Inversion of Coefficients in Partial Differential Equations
    Aristoff, David
    Bangerth, Wolfgang
    SIAM REVIEW, 2023, 65 (04) : 1074 - 1105
  • [47] IDENTIFICATION OF A PARAMETER IN FOURTH-ORDER PARTIAL DIFFERENTIAL EQUATIONS BY AN EQUATION ERROR APPROACH
    Bush, Nathan
    Jadamba, Baasansuren
    Khan, Akhtar A.
    Raciti, Fabio
    MATHEMATICA SLOVACA, 2015, 65 (05) : 1209 - 1221
  • [48] Numerical Method Research of Partial Differential Equations Inverse Problem
    He, Yali
    Peng, Yamian
    Feng, Lichao
    INTELLIGENT STRUCTURE AND VIBRATION CONTROL, PTS 1 AND 2, 2011, 50-51 : 455 - 458
  • [49] The partial boundary value conditions of nonlinear degenerate parabolic equation
    Zhi, Yuan
    Zhan, Huashui
    BOUNDARY VALUE PROBLEMS, 2022, 2022 (01)
  • [50] Difference inversion model of wave equation
    Wang De-ming
    APPLIED MATHEMATICS AND MECHANICS-ENGLISH EDITION, 2008, 29 (03) : 361 - 366