A parallel finite-difference approach for 3D transient electromagnetic modeling with galvanic sources

A parallel finite-difference algorithm for the solution of diffusive, three-dimensional (3D) transient electromagnetic field simulations is presented. The purpose of the scheme is the simulation of both electric fields and the time derivative of magnetic fields generated by galvanic sources (grounded wires) over arbitrarily complicated distributions of conductivity and magnetic permeability. Using a staggered grid and a modified DuFort-Frankel method, the scheme steps Maxwell's equations in time. Electric field initialization is done by a conjugate-gradient solution of a 3D Poisson problem, as is common in 3D resistivity modeling. Instead of calculating the initial magnetic field directly, its time derivative and curl are employed in order to advance the electric field in time. A divergence-free condition is enforced for both the magnetic-field time derivative and the total conduction-current density, providing accurate results at late times. In order to simulate large realistic earth models, the algorithm has been designed to run on parallel computer platforms. The upward continuation boundary condition for a stable solution in the infinitely resistive air layer involves a two-dimensional parallel fast Fourier transform. Example simulations are compared with analytical, integral-equation and spectral Lanczos decomposition solutions and demonstrate the accuracy of the scheme.

[1]  Three-dimensional inversion of transient-electromagnetic data: A comparative study , 2003 .

[2]  G. W. Hohmann,et al.  A finite-difference, time-domain solution for three-dimensional electromagnetic modeling , 1993 .

[3]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[4]  Neubauer,et al.  A first attempt at monitoring underground gas storage by means of time‐lapse multichannel transient electromagnetics , 2000 .

[5]  K. Strack,et al.  Exploration with Deep Transient Electromagnetics , 1992 .

[6]  Hiroshi Amano,et al.  2.5‐D inversion of frequency‐domain electromagnetic data generated by a grounded‐wire source , 2002 .

[7]  C. H. Stoyer,et al.  Finite-difference calculations of the transient field of an axially symmetric Earth for vertical magnetic dipole excitation , 1983 .

[8]  Anthony Skjellum,et al.  Writing libraries in MPI , 1993, Proceedings of Scalable Parallel Libraries Conference.

[9]  E. Haber.  A mixed finite element method for the solution of the magnetostatic problem with highly discontinuous coefficients in 3D , 2000 .

[10]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[11]  H. Scriba,et al.  COMPUTATION OF THE ELECTRIC POTENTIAL IN THREE‐DIMENSIONAL STRUCTURES* , 1981 .

[12]  Gerald W. Hohmann,et al.  Diffusion of electromagnetic fields into a two-dimensional earth; a finite-difference approach , 1984 .

[13]  Bülent Tezkan,et al.  Interpretation of long‐offset transient electromagnetic data from the Odenwald area, Germany, using two‐dimensional modelling , 2000 .

[14]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[15]  L. Knizhnerman,et al.  Spectral approach to solving three-dimensional Maxwell's diffusion equations in the time and frequency domains , 1994 .

[16]  J. Shadid,et al.  Three‐dimensional wideband electromagnetic modeling on massively parallel computers , 1996 .

[17]  Gregory A. Newman,et al.  Three‐dimensional massively parallel electromagnetic inversion—I. Theory , 1997 .

[18]  Survey design for multicomponent electromagnetic systems , 1984 .

[19]  Gregory A. Newman,et al.  Transient electromagnetic response of a three-dimensional body in a layered earth , 1986 .

[20]  Michael Commer,et al.  New advances in three dimensional transient electromagnetic inversion , 2004 .

[21]  S. Frankel,et al.  Stability conditions in the numerical treatment of parabolic differential equations , 1953 .

[22]  R. N. Edwards,et al.  2. The Magnetometric Resistivity Method , 1991 .

[23]  A. Dey,et al.  Resistivity modeling for arbitrarily shaped three-dimensional structures , 1979 .

[24]  K. Spitzer.  A 3-D FINITE-DIFFERENCE ALGORITHM FOR DC RESISTIVITY MODELLING USING CONJUGATE GRADIENT METHODS , 1995 .

[25]  Frank Dale Morgan,et al.  Modeling of streaming potential responses caused by oil well pumping , 1994 .

[26]  Vladimir Druskin,et al.  INTERPRETATION OF 3-D EFFECTS IN LONG-OFFSET TRANSIENT ELECTROMAGNETIC (LOTEM) SOUNDINGS IN THE MUNSTERLAND AREA/GERMANY , 1992 .

[27]  Douglas W. Oldenburg,et al.  3D forward modelling of time domain electromagnetic data , 2002 .

[28]  G. Pinder,et al.  Numerical solution of partial differential equations in science and engineering , 1982 .

[29]  Jopie I. Adhidjaja,et al.  A finite-difference algorithm for the transient electromagnetic response of a three-dimensional body , 1989 .

[30]  Sofia Davydycheva,et al.  A Finite Difference Scheme for Elliptic Equations with Rough Coefficients Using a Cartesian Grid Nonconforming to Interfaces , 1999 .