Theoretical Developments in Electromagnetic Induction Geophysics with Selected Applications in the Near Surface

Near-surface applied electromagnetic geophysics is experiencing an explosive period of growth with many innovative techniques and applications presently emergent and others certain to be forthcoming. An attempt is made here to bring together and describe some of the most notable advances. This is a difficult task since papers describing electromagnetic induction methods are widely dispersed throughout the scientific literature. The traditional topics discussed herein include modeling, inversion, heterogeneity, anisotropy, target recognition, logging, and airborne electromagnetics (EM). Several new or emerging techniques are introduced including landmine detection, biogeophysics, interferometry, shallow-water electromagnetics, radiomagnetotellurics, and airborne unexploded ordnance (UXO) discrimination. Representative case histories that illustrate the range of exciting new geoscience that has been enabled by the developing techniques are presented from important application areas such as hydrogeology, contamination, UXO and landmines, soils and agriculture, archeology, and hazards and climat

Doctor of Philosophy

dissertationModeling of induced polarization (IP) phenomena is important for developing effective methods for remote sensing of subsurface geology and is widely used in mineral exploration. However, the quantitative interpretation of IP data in a complex 3D environment is still a challenging problem of applied geophysics. In this dissertation I use the regularized conjugate gradient method to determine the 3D distribution of the four parameters of the Cole-Cole model based on surface induced polarization (IP) data. This method takes into account the nonlinear nature of both electromagnetic induction (EMI) and IP phenomena. The solution of the 3D IP inverse problem is based on the regularized smooth inversion only. The method was tested on synthetic models with DC conductivity, intrinsic chargeability, time constant, and relaxation parameters, and it was also applied to the practical 3D IP survey data. I demonstrate that the four parameters of the Cole-Cole model, DC electrical resistivity, p0 (or electrical conductivity <r0 = 1 / p 0 ), chargeability, r ; time constant, r ; and the relaxation parameter, , can be recovered from the observed IP data simultaneously. There are four Cole-Cole parameters involved in the inversion, in other words, within each cell, there are DC conductivity ( ), chargeability ( ), time parameters ( ), and relaxation parameters ( ) compared to conductivity only, used in EM only inversion. In addition to more inversion parameters used in IP survey, dipole-dipole configuration which requires more sources and receivers. One the other hand, calculating Green tensor and Frechet matrix time consuming and storing them requires a lot of memory. So, I develop parallel computation using MATLAB parallel tool to speed up the calculation

Model design for algorithmic efficiency in electromagnetic sensing

The objective of the proposed research is to develop structural changes to the design and application of electromagnetic (EM) sensing models to more efficiently and accurately invert EM measurements to extract parameters for applications such as landmine detection. Two different acquisition modalities are addressed in this research: ground-penetrating radar (GPR) and electromagnetic induction (EMI) sensors. The models needed for practical three-dimensional (3D) spatial imaging typically become impractically large, with up to seven dimensions of parameters that need to be extracted. These parameters include, but are not limited to target type, 3D location, and 3D orientation. The new special structures for these models exploit properties such as shift invariance and tensor representation, which can be combined with strategic inversion techniques, including the Fast Fourier Transform and semidefinite programming. The structures dramatically reduce the amount of computation and can eliminate the need to store up to five dimensions of parameters while still accurately estimating them.Ph.D

Classification of Metallic Targets Using a Walk-Through Metal Detection Portal

Metal detectors have been used for a long time for treasure hunting, security screening, and finding buried objects such as landmines or unexploded ordnance. Walk-through metal detection (WTMD) portals are used for making sure that forbidden or threatening metallic items, such as knives or guns, are not carried into secure areas at critical locations such as airports, court rooms, embassies, and prisons.The 9/11 terrorist act has given rise to stricter rules for aviation security worldwide, and the ensuing tighter security procedures have meant that passengers face more delays at airports. Moreover, the fear of terrorism has led to the adoption of security screening technology in a variety of places such as railway and coach stations, sports events, malls, and nightclubs.However, the current WTMD technology and scanning procedures at airports require that all metallic items be removed from clothing prior to scanning, causing inconvenience. Furthermore, alarms are triggered by innocuous items such as shoe shanks and artificial joints, along with overlooked items such as jewellery and belts. These lead to time- consuming, manual pat-down searches, which are found inconvenient, uncomfortable, and obtrusive by some.Modern WTMD portals are very sensitive devices that can detect items with only small amounts of metal, but they currently lack the ability to further classify the detected item. However, if a WTMD portal were able to classify objects reliably into, e.g., “knives”, “belts”, “keys”, the need for removing the items prior to screening would disappear, enabling a paradigm shift in the field of security screening.This thesis is based on novel research presented in five peer-reviewed publications. The scope of the problem has been narrowed down to a situation in which only one metallic item is carried through the portal at a time. However, the methods and results presented in this thesis can be generalized into a multi-object scenario. It has been shown that by using a WTMD portal and the magnetic polarisability tensor, it is possible to accurately distinguish between threatening and innocuous targets and to classify them into 10 to 13 arbitrary classes. Furthermore, a data library consisting of natural walk-throughs has been collected, and it has been demonstrated that the walk-through data collected with the above portal are subject to phenomena that might affect classification, in particular a bias and the so-called body effect. However, the publications show that, by using realistic walk-through data, high classification accuracy can be maintained regardless of the above problems. Furthermore, a self-diagnostics method for detecting unreliable samples has also been presented with potential to significantly increase classification accuracy and the reliability of decision making.The contributions presented in this thesis have a variety of implications in the field of WTMD-based security screening. The novel technology offers more information, such as an indication of the probable cause of the alarm, to support the conventional screening procedure. Moreover, eliminating the need for removing all metallic items prior to screening enables design of new products for scenarios such as sports events, where conventional screening procedures might be inconvenient, creating thus new business possibilities for WTMD manufacturing companies.The positive results give rise to a variety of future research topics such as using wideband data, enabling simultaneous classification of multiple objects, and developing the portal coil design to diminish signal nonlinearities. Furthermore, the ideas and the basic principles presented in this thesis may be applied to other metal detection applications, such as humanitarian demining

Geophysics for Mineral Exploration

This Special Issue contains ten papers which focus on emerging geophysical techniques for mineral exploration, novel modeling, and interpretation methods, including joint inversions of multi physics data, and challenging case studies. The papers cover a wide range of mineral deposits, including banded iron formations, epithermal gold–silver–copper–iron–molybdenum deposits, iron-oxide–copper–gold deposits, and prospecting forgroundwater resources

Classification, identification, and modeling of unexploded ordnance in realistic environments

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2008.Includes bibliographical references (p. 205-218).Recovery of buried unexploded ordnance (UXO) is very slow and expensive due to the high false alarm rate created by clutter. Electromagnetic induction (EMI) has been shown to be a promising technique for UXO detection and discrimination. This thesis uses the EMI response of buried targets to identify or classify them. To perform such discrimination, accurate forward models of buried UXO are needed. This thesis provides a survey of existing target models: the dipole model, the spheroid model, and the fundamental mode model. Then the implementation of a new model, the spheroidal mode model, is described and validated against measurements of a UXO. Furthermore, an in-depth study of the effects of permeable soil, modeled as a permeable half space, is presented. This study concludes that the discontinuity created by the air to permeable soil interface produces minimal effect in the response of a buried object. The change is limited to a magnitude shift of the real portion of the EMI response and can be reproduced by superposition of a permeable half space response on the response of the same object in frees pace. Accurate soil modeling also allows one to invert for soil permeability values from measured data if such data are in known units. However, the EMI sensor used in this study provides measurements in consistent but unknown units. Furthermore, the instrument is from a third party and is proprietary. Therefore, this thesis describes the development of a non-invasive method to model and calibrate non-adaptive instruments so that all measurements can be converted into units consistent with modeled data. This conversion factor is shown to be a constant value across various conditions, thus demonstrating its validity.(cont.) Given that now a more complete model of the measurable response of a buried UXO is implemented, this study proceeds to demonstrate that EMI responses from UXO and clutter objects can be used to identify the objects through the application of Differential Evolution (DE), a type of Genetic Algorithm. DE is used to optimize the parameters of the UXO fundamental mode model to produce a match between the modeled response and the measured response of an unknown object. When this optimization procedure is applied across a library of models for possible UXO, the correct identity of the unknown object can be ascertained because the corresponding library member will produce the closest match. Furthermore, responses from clutter objects are shown to produce very poor matches to library objects, thus providing a method to discriminate UXO from clutter. These optimization experiments are conducted on measurements of UXO in air, UXO in air but obscured by clutter fragments, buried UXO, and buried UXO obscured by clutter fragments. It is shown that the optimization procedure is successful for shallow buried objects obscured by light clutter contributing to roughly 20 dB SNR, but is limited in applicability towards very deeply buried UXO or those in dense clutter environments. The DE algorithm implemented in this study is parallelized and the optimization results are computed with a multi-processor supercomputer. Thus, the computational requirement of DE is a considerable drawback, and the method cannot be used for real time, on-site inversion of measured UXO data. To address this concern, a different approach to inversion is also implemented in this study. Rather than identifying particular UXO, one may do a discrimination between general UXO and general clutter items. Previous work has shown that the expansion coefficients of EMI responses in the spheroidal coordinate system can uniquely characterize the corresponding targets.(cont.) Therefore, these coefficients readily lend themselves for use as features by which objects can be classified as likely to be UXO or unlikely to be UXO. To do such classification, the relationship between these coefficients and the physical properties of UXO and clutter, such as differences in size or body-of-revolution properties or material heterogeneity properties, must be found. This thesis shows that such relationships are complex and require the use of the automated pattern recognition capability of machine learning. Two machine learning algorithms, Support Vector Machines and Neural Networks, are used to identify whether objects are likely to be UXO. Furthermore, the effects of small diffuse clutter fragments and uncertainty about the target position are investigated. This discrimination procedure is applied on both synthetic data from models and measurements of UXO and clutter. It is found that good discrimination is possible for up to 20 dB SNR. But the discrimination is sensitive to inaccurate estimations of a target's depth. It is found that the accuracy must be within a 10 cm deviation of an object's true depth. The general conclusion forwarded by this work is that while increasingly accurate discrimination capabilities can be produced through more detailed forward modeling and application of robust optimization and learning algorithms, the presence of noise and clutter is still of great concern. Minimization or filtering of such noise is necessary before field deployable discrimination techniques can be realized.by Beijia Zhang.Ph.D

Parallel 3-D marine controlled-source electromagnetic modelling using high-order tetrahedral Nédélec elements

We present a parallel and high-order Nédélec finite element solution for the marine controlled-source electromagnetic (CSEM) forward problem in 3-D media with isotropic conductivity. Our parallel Python code is implemented on unstructured tetrahedral meshes, which support multiple-scale structures and bathymetry for general marine 3-D CSEM modelling applications. Based on a primary/secondary field approach, we solve the diffusive form of Maxwell’s equations in the low-frequency domain. We investigate the accuracy and performance advantages of our new high-order algorithm against a low-order implementation proposed in our previous work. The numerical precision of our high-order method has been successfully verified by comparisons against previously published results that are relevant in terms of scale and geological properties. A convergence study confirms that high-order polynomials offer a better trade-off between accuracy and computation time. However, the optimum choice of the polynomial order depends on both the input model and the required accuracy as revealed by our tests. Also, we extend our adaptive-meshing strategy to high-order tetrahedral elements. Using adapted meshes to both physical parameters and high-order schemes, we are able to achieve a significant reduction in computational cost without sacrificing accuracy in the modelling. Furthermore, we demonstrate the excellent performance and quasi-linear scaling of our implementation in a state-of-the-art high-performance computing architecture.This project has received funding from the European Union's Horizon 2020 programme under the Marie Sklodowska-Curie grant agreement No. 777778. Furthermore, the research leading to these results has received funding from the European Union's Horizon 2020 programme under the ChEESE Project (https://cheese-coe.eu/ ), grant agreement No. 823844. In addition, the authors would also like to thank the support of the Ministerio de Educación y Ciencia (Spain) under Projects TEC2016-80386-P and TIN2016-80957-P. The authors would like to thank the Editors-in-Chief and to both reviewers, Dr. Martin Cuma and Dr. Raphael Rochlitz, for their valuable comments and suggestions which helped to improve the quality of the manuscript. This work benefited from the valuable suggestions, comments, and proofreading of Dr. Otilio Rojas (BSC). Last but not least, Octavio Castillo-Reyes thanks Natalia Gutierrez (BSC) for her support in CSEM modeling with BSIT.Peer ReviewedPostprint (author's final draft