A finite element approach to the 3D CSEM modeling problem and applications to the study of the effect of target interaction andtopography

The solution of the secondary coupled-vector potential formulation of Maxwell??s equations governing the controlled-source electromagnetic (CSEM) response of an arbitrary, threedimensionalconductivitymodelmust be calculatednumerically.The finite elementmethod is attractive, because it allows the model to be discretized into an unstructured mesh, permitting the specification of realistic irregular conductor geometries, and permitting the mesh to be refined locally, where finer resolution is needed. The calculated results for a series ofsimple test problems, ranging from one-dimensionalscalar differentialequations to three-dimensional coupled vector equations match the known analytic solutions well, with error values several orders of magnitude smaller than the calculated values. The electromagnetic fields of a fully three-dimensional CSEM model, recovered from the potentials using the moving least squares interpolation numerical differentiation algorithm, compares well with published numerical modeling results, particularly when local refinement is applied. Multiple buried conductors in a conductive host interact via mutual induction and current flow through the host due to the dissipation of charge accumulated on the conductor boundary. The effect of this interaction varies with host conductivity, transmitter frequency, and conductor geometry, orientation, and conductivity. For three test models containingtwo highly conductive plate-like targets, oriented in various geometries (parallel, perpendicular, and horizontal), mutual coupling ranges as high as twenty times the total magnetic field. The effect of varying host conductivity is significant, especially at high frequencies. Numerical modeling also shows that the vorticity of the currents density induced in a vertically oriented plate-like conductor rotates from vertical at high frequencies, to horizontal at low frequencies, a phenomenon confirmed by comparison with time domain field data collected in Brazos County, Texas. Furthermore, the effect of the presence of a simple horst on the CSEM response of a homogeneous conductive earth is significant, even when the height of the horst is only a fraction of the skin depth of the model. When the transmitter is placedon topofthe horst, the currents inducedtherein account for nearly all of the total magnetic field of the model, indicating that topography, like mutual coupling must be accounted for when interpreting CSEM data

The effects of cultural noise on controlled source electromagnetic resonses of subsurface fractures in resistive terrain

Controlled source electromagnetic (CSEM) geophysics has been used with a fair amount of success in near surface hydrogeological studies. Recently, these investigations have been conducted frequently in human impacted field sites containing cultural conductors such as metal fences and buried pipes. Cultural noise adds an element of complexity to the geological interpretation of this type of data. This research investigates the influence of mutual induction between two buried targets in a CSEM experiment. In particular, it looks at the mutual coupling between a buried cultural conductor and a geological heterogeneity. We attempt to isolate the Hz field induced by tertiary currents in targets caused by mutual coupling. This is achieved with a Texas A&M 3D CSEM finite element code, which calculates the secondary Hz fields emanating from a target buried in a halfspace. Buried geological targets and cultural conductors are modeled as volumetric slabs embedded in a halfspace. A series of models have been simulated to study the effect of varying parameters such as target conductivity, transmitter location and shape of a target on the mutual inductance. In each case, the secondary Hz field is calculated for a model with two slabs, and two models with individual slabs. The mutual coupling is calculated by removing the secondary fields from the individual slab models from the response of a two slab model. The calculations of mutual inductance from a variety of such models suggests a complicated interaction of EM fields between the two targets. However, we can explain most of these complexities by adapting a simple approach to Maxwell’s equations. Although the tertiary Hz field is complicated, it may be useful in the characterization and delineation of electrical heterogeneities in the subsurface, which can then be related to geological features such as fractures or joints. It is seen that the most important factor affecting the mutual coupling is the host conductivity. The results have also shown that mutual coupling is very sensitive to transmitter (TX) location, especially when the TX is positioned near one of the targets

Three Dimensional Controlled-source Electromagnetic Edge-based Finite Element Modeling of Conductive and Permeable Heterogeneities

Presence of cultural refuse has long posed a serious challenge to meaningful geological interpretation of near surface controlled–source electromagnetic data (CSEM). Cultural refuse, such as buried pipes, underground storage tanks, unexploded ordnance, is often highly conductive and magnetically permeable. Interpretation of the CSEM response in the presence of cultural noise requires an understanding of electromagnetic field diffusion and the effects of anomalous highly conductive and permeable structures embedded in geologic media. While many numerical techniques have been used to evaluate the response of three dimensional subsurface conductivity distributions, there is a lack of approaches for modeling the EM response incorporating variations in both subsurface conductivity σ and relative permeability μr. In this dissertation, I present a new three dimensional edge–based finite element (FE) algorithm capable of modeling the CSEM response of buried conductive and permeable targets. A coupled potential formulation for variable μ using the vector magnetic potential A and scalar electric potential V gives rise to an ungauged curl–curl equation. Using reluctivity (v=1/mu ), a new term in geophysical applications instead of traditional magnetic susceptibility, facilitates a separation of primary and secondary potentials. The resulting differential equation is solved using the finite element method (FEM) on a tetrahedral mesh with local refinement capabilities. The secondary A and V potentials are expressed in terms of the vector edge basis vectors and the scalar nodal basis functions respectively. The finite element matrix is solved using a Jacobi preconditioned QMR solver. Post processing steps to interpolate the vector potentials on the nodes of the mesh are described. The algorithm is validated against a number of analytic and multi dimensional numeric solutions. The code has been deployed to estimate the influence of magnetic permeability on the mutual coupling between multiple geological and cultural targets. Some limitations of the code with regards to speed and performance at high frequency, conductivity and permeability values have been noted. Directions for further improvement and expanding the range of applicability have been proposed

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

A Multidisciplinary Analysis of Frequency Domain Metal Detectors for Humanitarian Demining

This thesis details an analysis of metal detectors (low frequency electromagnetic induction devices) with emphasis on Frequency Domain (FD) systems and the operational conditions of interest to humanitarian demining. After an initial look at humanitarian demining and a review of their basic principles we turn our attention to electromagnetic induction modelling and to analytical solutions to some basic FD direct (forward) problems. The second half of the thesis focuses then on the analysis of an extensive amount of experimental data. The possibility of target classification is first discussed on a qualitative basis, then quantitatively. Finally, we discuss shape and size determination via near field imaging

Study on 3D forward modeling & inversion of surface-borehole electromagnetic data.

The purpose of this research is to develop an interpretive tool to meet the requirements of deep mineral exploration. Therefore, we carried out a series of research work as part of a doctoral training program and achieved the relevant objectives below. The core of this doctoral thesis is the development of 3D modeling tools to interpret the electromagnetic data collected in boreholes. First, a 3D model creation tool is designed, with which we can easily build a 3D geological model from sections and quickly discretize it. The sections could be true geological cross-sections or from a conceptual geological model. The utility of this tool is to facilitate the tests of the algorithms developed within the framework of this thesis, in order to model the electromagnetic responses in various geological situations and allow to easily change the parameters of the geophysical measurement system. Two parallelization algorithms, MPI-based and hybrid MPI/OpenMP-based methods, are designed for surface borehole time domain electromagnetic (BHTEM) forward modeling. The BHTEM responses are calculated from anomalous regions distributed in a 3D model (discretized into cells). The forward modeling additionally uses multiple meshes, fine meshes are used for the anomalous region in the high-frequency range and coarser meshes for geological background in the low-frequency range. Based on varying meshes for different frequency ranges, the parallel computation greatly reduces the computation time of the TEM forward modeling. An optimal survey design benefits from quick forward modeling. We found that the target BHTEM response depends upon the transmitting pulse width, target time constant, and the duration of measurement time. We proposed the formula with respect to the three variables to design optimal pulse widths in advance for different off-times in order to maximize the efficiency of TEM measurement in the field. Finally, a 3D BHTEM inversion algorithm is developed based on the Gauss-Newton method with high spatial resolution. By introducing the isosurface, neighborhood anomalies search, 3D trace envelope, and false targets elimination into the inversion process, the predicted model is improved through iterations and interactions between the computation and the user intervention

Imaging of buried utilities by ultra wideband sensory systems

Third-party damage to the buried infrastructure like natural gas pipelines, water distribution pipelines and fiber optic cables are estimated at $10 billion annually across the US. Also, the needed investment in upgrading our water and wastewater infrastructure over the next 20 years is estimated by Environmental Protection Agency (EPA) at$400 billion, however, non-destructive condition assessment technologies capable of providing quantifiable data regarding the structural integrity of our buried assets in a cost-effective manner are lacking. Both of these areas were recently identified several U.S. federal agencies as \u27critical national need\u27. In this research ultra wideband (UWB) time-domain radar technology was adopted in the development of sensory systems for the imaging of buried utilities, with focus on two key applications. The first was the development of a sensory system for damage avoidance of buried pipes and conduits during excavations. A sensory system which can be accommodated within common excavator buckets was designed, fabricated and subjected to laboratory and full-scale testing. The sensor is located at the cutting edge (teeth), detecting the presence of buried utilities ahead of the cutting teeth. That information can be used to alert the operator in real-time, thus avoiding damage to the buried utility. The second application focused on a sensory system that is capable of detecting structural defects within the wall of buried structures as well as voids in the soil-envelope encasing the structure. This ultra wideband sensory system is designed to be mounted on the robotic transporter that travels within the pipeline while collecting data around the entire circumference. The proposed approach was validated via 3-D numerical simulation as well as full-scale experimental testing

Experimental time-domain controlled source electromagnetic induction for highly conductive targets detection and discrimination

The response of geological materials at the scale of meters and the response of buried targets of different shapes and sizes using controlled-source electromagnetic induction (CSEM) is investigated. This dissertation focuses on three topics; i) frac- tal properties on electric conductivity data from near-surface geology and processing techniques for enhancing man-made target responses, ii) non-linear inversion of spa- tiotemporal data using continuation method, and iii) classification of CSEM transient and spatiotemporal data. In the first topic, apparent conductivity profiles and maps were studied to de- termine self-affine properties of the geological noise and the effects of man-made con- ductive metal targets. 2-D Fourier transform and omnidirectional variograms showed that variations in apparent conductivity exhibit self-affinity, corresponding to frac- tional Brownian motion. Self-affinity no longer holds when targets are buried in the near-surface, making feasible the use of spectral methods to determine their pres- ence. The difference between the geology and target responses can be exploited using wavelet decomposition. A series of experiments showed that wavelet filtering is able to separate target responses from the geological background. In the second topic, a continuation-based inversion method approach is adopted, based on path-tracking in model space, to solve the non-linear least squares prob- lem for unexploded ordnance (UXO) data. The model corresponds to a stretched- exponential decay of eddy currents induced in a magnetic spheroid. The fast inversion of actual field multi-receiver CSEM responses of inert, buried ordnance is also shown. Software based on the continuation method could be installed within a multi-receiver CSEM sensor and used for near-real-time UXO decision. In the third topic, unsupervised self-organizing maps (SOM) were adapted for data clustering and classification. The use of self-organizing maps (SOM) for central- loop CSEM transients shows potential capability to perform classification, discrimi- nating background and non-dangerous items (clutter) data from, for instance, unex- ploded ordnance. Implementation of a merge SOM algorithm showed that clustering and classification of spatiotemporal CSEM data is possible. The ability to extract tar- get signals from a background-contaminated pattern is desired to avoid dealing with forward models containing subsurface response or to implement processing algorithm to remove, to some degree, the effects of background response and the target-host interactions