
Related resources
Search for item elsewhere
University researcher(s)
Academic department(s)
Better Imaging for Landmine Detection: An exploration of 3D full-wave inversion for ground-penetrating radar
[Thesis]. Manchester, UK: The University of Manchester; 2016.
Access to files
-
FULL-TEXT.PDF (pdf)
Abstract
Humanitarian clearance of minefields is most often carried out by hand, conventionally using a a metal detector and a probe. Detection is a very slow process, as every piece of detected metal must treated as if it were a landmine and carefully probed and excavated, while many of them are not. The process can be safely sped up by use of Ground-Penetrating Radar (GPR) to image the subsurface, to verify metal detection results and safely ignore any objects which could not possibly be a landmine. In this thesis, we explore the possibility of using Full Wave Inversion (FWI) to improve GPR imaging for landmine detection. Posing the imaging task as FWI means solving the large-scale, non-linear and ill-posed optimisation problem of determining the physical parameters of the subsurface (such as electrical permittivity) which would best reproduce the data. This thesis begins by giving an overview of all the mathematical and implementational aspects of FWI, so as to provide an informative text for both mathematicians (perhaps already familiar with other inverse problems) wanting to contribute to the mine detection problem, as well as a wider engineering audience (perhaps already working on GPR or mine detection) interested in the mathematical study of inverse problems and FWI.We present the first numerical 3D FWI results for GPR, and consider only surface measurements from small-scale arrays as these are suitable for our application. The FWI problem requires an accurate forward model to simulate GPR data, for which we use a hybrid finite-element boundary-integral solver utilising first order curl-conforming N\'{e}d\'{e}lec (edge) elements. We present a novel `line search' type algorithm which prioritises inversion of some target parameters in a region of interest (ROI), with the update outside of the area defined implicitly as a function of the target parameters. This is particularly applicable to the mine detection problem, in which we wish to know more about some detected metallic objects, but are not interested in the surrounding medium. We may need to resolve the surrounding area though, in order to account for the target being obscured and multiple scattering in a highly cluttered subsurface.We focus particularly on spatial sensitivity of the inverse problem, using both a singular value decomposition to analyse the Jacobian matrix, as well as an asymptotic expansion involving polarization tensors describing the perturbation of electric field due to small objects. The latter allows us to extend the current theory of sensitivity in for acoustic FWI, based on the Born approximation, to better understand how polarization plays a role in the 3D electromagnetic inverse problem. Based on this asymptotic approximation, we derive a novel approximation to the diagonals of the Hessian matrix which can be used to pre-condition the GPR FWI problem.
Layman's Abstract
Humanitarian clearance of minefields is most often carried out by hand, conventionally using a a metal detector and a probe. Detection is a very slow process, as every piece of detected metal must treated as if it were a landmine and carefully probed and excavated, while many of them are not. The process can be safely sped up by use of Ground-Penetrating Radar (GPR) to image the subsurface, to verify metal detection results and safely ignore any objects which could not possibly be a landmine. In this thesis, we explore the possibility of using Full Wave Inversion (FWI) to improve GPR imaging for landmine detection. Posing the imaging task as FWI means solving the large-scale, non-linear and ill-posed optimisation problem of determining the physical parameters of the subsurface (such as electrical permittivity) which would best reproduce the data. We consider the mathematical and implementational aspects of such an approach, presenting the first FWI results for GPR in 3D.
Additional content not available electronically
MATLAB code for simulating GPR data, solving the inverse (FWI) problem, and approximating scattered fields with polarization tensors. DOI 10.5281/zenodo.31437
Keyword(s)
3D; BEM; Boundary integral; Computational electromagnetics; Electromagnetics; FE-BI; FEBI; FEM; FWI; Finite element; Full-wave inversion; Full-waveform inversion; GPR; Ground penetrating radar; Helmholtz equation; Hessian; Ill-posed; Imaging; Inverse problems; Inverse scattering; Maxwell's equations; Nedelec elements; Non-linear; Nuisance parameters; Optimisation; Preconditioning; ROI; RWG; Region of interest; Regularization; SVD; Scattering; Sensitivity; Shape sensitivity; Spatial sensitivity; Surface equivalence; TV; Total variation; Vector wave equation; asymptotic; edge elements; l-BFGS; landmine detection; polarizability tensor; polarization tensor
Bibliographic metadata
- Inverse problems
- Ill-posed
- Non-linear
- Imaging
- Full-wave inversion
- FWI
- Full-waveform inversion
- SVD
- polarization tensor
- polarizability tensor
- asymptotic
- landmine detection
- Ground penetrating radar
- GPR
- 3D
- Helmholtz equation
- Electromagnetics
- Maxwell's equations
- Computational electromagnetics
- Vector wave equation
- Surface equivalence
- Finite element
- Nedelec elements
- RWG
- edge elements
- Boundary integral
- BEM
- FEM
- FE-BI
- FEBI
- Scattering
- Inverse scattering
- Regularization
- Total variation
- TV
- Optimisation
- Sensitivity
- Preconditioning
- Spatial sensitivity
- Shape sensitivity
- Hessian
- Nuisance parameters
- l-BFGS
- ROI
- Region of interest