3-D Finite-Difference Modeling of the Magnetic Field in Geoelectromagnetic Induction
J. T. Weaver, A. K. Agarwal, X. H. Pu, 1999. "3-D Finite-Difference Modeling of the Magnetic Field in Geoelectromagnetic Induction", Three-Dimensional Electromagnetics, Michael Oristaglio, Brian Spies
Download citation file:
A versatile finite-difference program has been developed for solving 3-D magnetotelluric problems in which the inducing field is uniform and horizontal. Finite-difference equations for the magnetic-field components are obtained by integrations over rectangular volume elements surrounding each node of the grid. This procedure avoids the ambiguity about whether to define weighted average conductivities or resistivities at the nodes of the numerical grid. The algorithm can handle an anomalous thin sheet at the surface of the general 3-D structure, thereby modeling both near-surface and deep-seated conductivity anomalies for certain ranges of period. The anomalous structure is allowed to approach different 2-D configurations at infinity in all four horizontal directions, with the result that both E-polarization and B-polarization solutions arise as boundary conditions on the grid. Integral boundary conditions on the surface of the Earth and at the top of the basement half-space (of uniform conductivity) reduce the overall size of the grid by eliminating the air and basement layers from the numerical solution. Recent developments include an automatic grid generator designed for geoelectromagnetic induction problems, the treatment of difficulties that can occur when the 2-D E-polarization problems on the side boundaries of the grid are solved in terms of the magnetic field, and a new method for calculating the electric-field components.
Figures & Tables
In 1975 Jerry Hohmann published a paper1 that described his numerical implementation of an integral-equation method for three-dimensional electromagnetic (3-D EM2) modeling. The matrix equation for the simple model that he studied—a half-space containing a rectangular body discretized into 100 cubic cells—barely fit into the computer (a UNIVAC 1108 at the University of Utah). Coaxing interesting and correct results from the model and method clearly comprised much of the art and fun of the paper. And winding through the paper’s 50 or so equations and nearly 20 figures was a clear message: 3-D EM is different!
Three-dimensional electromagnetics is qualitatively different with new phenomena3 and new challenges to our understanding of how electromagnetic fields interact with Earth and other conductive bodies (including our own). In subsequent years, Jerry with his students and colleagues pursued these challenges across many fields—mining geophysics, geothermal exploration, magnetotelluric crustal studies, environmental geophysics, oil and gas exploration—in both the time and frequency domains. Of his 51 articles4 in journals and monographs, more than half dealt with three-dimensional electromagnetics.
In 1995, 20 years after Jerry’s classic paper (and three years after his death from cancer in May, 1992), nearly 200 scientists from around the world gathered at Schlumberger–Doll Research in Ridgefield, Connecticut, for a symposium in his memory, the (first) International Symposium on Three-Dimensional Electromagnetics. More than 70 papers were presented in oral and poster sessions during three days organized around the themes: Modeling, Inversion, and Practice. The quality of the work presented, the liveliness of the discussions, and the demand for the symposium proceedings were the impetus for this new volume. We invited the authors to submit longer, more tutorial versions of their articles for a book to be published by the Society of Exploration Geophysicists (SEG) in the series Geophysical Developments.
As is evident from the size of this volume, we were overwhelmed by the response. We hope that readers will find the contents equally weighty. The 44 articles collected here are the work of 97 authors, representing 55 different institutions (universities, government or industrial research labs) from 13 countries around the world. All have been reviewed and edited according to the strict standards of SEG’s lead journal, Geophysics. They represent the state of the art in 3-D EM at the time final revisions were received (from the fall of 1997 through the spring of 1998).
The lead article addresses one of Jerry’s favorite subjects—the need for independent checks on any numerical calculation; it shows how far we have come since 1975 and how far we still are from routine, confident use of 3-D EM models. We have grouped the remaining articles into nine sections:
3-D EM and parallel computers
Magnetotellurics and global induction
Mining and exploration geophysics
Borehole geophysics and logging
This division into techniques and applications is naturally very rough; many articles could easily appear in two or three different sections. The subjects covered in this volume touch, we believe, on every major technique being used today to compute, analyze, visualize, and understand 3-DEM fields in every major application of electrical geophysics (and in two applications outside geophysics: the interaction of 3-DEM fields with the human body and the non-destructive testing of aircraft). The late 1980’s saw the rapid development of 3-D seismics, which has revolutionized exploration for oil and gas in the 1990’s. The early years of the new millenium may see another revolution brought about by the rapid advances now occurring in 3-D EM.
Ridgefield, Connecticut USA
23 May 1998