Discontinuous Galerkin frequency domain forward modelling for the inversion of electric permittivity in the 2D case

Abstract

A B S T R A C TWe have recently developed a discontinuous Galerkin frequency domain modelling algorithm for the solution of the 2D transverse magnetic Maxwell equations. This method is formulated on an unstructured triangular discretization of the computational domain and makes use of a high order polynomial interpolation of the electromagnetic field components within each triangular element. The discontinuous nature of the approximation naturally allows for a local definition of the interpolation order that is, in combination with a possibly non-conforming local refinement of the mesh, a key ingredient for obtaining a flexible and accurate discretization method. Moreover, heterogeneity of the propagation media is easily dealt with by assuming element-wise values of the electromagnetic parameters. In this paper, we propose the use of this discontinuous Galerkin frequency domain method as the forward modelling algorithm for solving the inverse problem for the electric permittivity in the 2D case. The inversion process is based on a gradient minimization technique developed by Pratt for seismological applications. Preliminary numerical results are presented for the imaging of a simplified subsurface model with the aim of assessing the performances of the proposed inversion methodology with regards to the number of frequencies, the number of recorded data and the number of sources.

Dual-Parameter Simultaneous Full Waveform Inversion of Ground-Penetrating Radar for Arctic Sea Ice

“…where λ is the regularization factor, and TV is the TV regularization operator. Since the derivative of TV is not continuous, its differentiability can be guaranteed by the following formula approximation: (18) where Ω is the imaging region. The gradient of TV regularization is shown as follows:…”

“…When FWI was initially applied to GPR, the single-parameter inversion of the permittivity was realized [18]. At present, the requirement of radar FWI has been improved from single-parameter inversion to multiparameter inversion, but this will lead to a series inversion phenomenon, and the target parameters and morphology cannot be reconstructed.…”

Optimized Schwarz Methods for the Time-Harmonic Maxwell Equations with Damping

“…We present now a second more realistic test which consists in the simulation of electromagnetic wave propagation in a heterogenous subsurface medium. This kind of simulation is very important in imaging, see [13] for details. The configuration of the subsurface is shown in Figure 6.2 and is constituted of media which are characterized by various parameter values ε and σ, ε = 1.5, 2.25, 3.25, 4.25 and σ = 10 −6 , 10 −5 , 10 −4 10 −3 .…”

Multi-parameter imaging by finite difference frequency domain full waveform inversion of GPR data: A guide for sedimentary architecture modeling

“…Along with this advancement, the development and implementation of FDFD gradient method were used in GPR data inversion (Yang et al, 2012). There are a few numbers of published articles on FWI for surface-GPR data, which considers only the permittivity modeling with fixed conductivity value (El Bouajaji et al, 2011). The great advancement in this method was done by (Lavoué et al, 2014).…”