Abstract This paper presents a method for inverting ground penetrating radargrams in terms of one-dimensional profiles. We resort to a special type of linearization of the damped E-field wave equation to solve the inverse problem. The numerical algorithm for the inversion is iterative and requires the solution of several forward problems, which we evaluate using the matrix propagation approach. Analytical expressions for the derivatives with respect to physical properties are obtained using the self-adjoint Green's function method. We consider three physical properties of materials; namely dielectrical permittivity, magnetic permeability and electrical conductivity. The inverse problem is solved minimizing the quadratic norm of the residuals using quadratic programming optimization. In the iterative process to speed up convergence we use the Levenberg–Mardquardt method. The special type of linearization is based on an integral equation that involves derivatives of the electric field with respect to magnetic permeability, electrical conductivity and dielectric permittivity; this equation is the result of analyzing the implication of the scaling properties of the electromagnetic field. The ground is modeled using thin horizontal layers to approximate general variations of the physical properties. We show that standard synthetic radargrams due to dielectric permittivity contrasts can be matched using electrical conductivity or magnetic permeability variations. The results indicate that it is impossible to differentiate one property from the other using GPR data.
[1]
Enrique Gómez-Treviño,et al.
1-D inversion of resistivity and induced polarization data for the least number of layers
,
1997
.
[2]
D. Oldenburg,et al.
METHODS FOR CALCULATING FRÉCHET DERIVATIVES AND SENSITIVITIES FOR THE NON‐LINEAR INVERSE PROBLEM: A COMPARATIVE STUDY1
,
1990
.
[3]
D. Marquardt.
An Algorithm for Least-Squares Estimation of Nonlinear Parameters
,
1963
.
[4]
Enrique Gómez-Treviño,et al.
Inversion of magnetotelluric soundings using a new integral form of the induction equation
,
1996
.
[5]
P. Gill,et al.
Fortran package for constrained linear least-squares and convex quadratic programming. User's Guide for LSSOL (Version 1. 0)
,
1986
.
[6]
E. Gómez-Treviño,et al.
Synthetic radargrams from electrical conductivity and magnetic permeability variations
,
1996
.
[7]
Kenneth Levenberg.
A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES
,
1944
.
[8]
Enrique Gómez-Treviño.
Nonlinear integral equations for electromagnetic inverse problems
,
1987
.
[9]
V. V. Rzhevskiĭ,et al.
The physics of rocks
,
1971
.