Diffusion of electromagnetic fields into a two-dimensional earth; a finite-difference approach

We describe a numerical method for time‐stepping Maxwell’s equations in the two‐dimensional (2-D) TE‐mode, which in a conductive earth reduces to the diffusion equation. The method is based on the classical DuFort‐Frankel finite‐difference scheme, which is both explicit and stable for any size of the time step. With this method, small time steps can be used at early times to track the rapid variations of the field, and large steps can be used at late times, when the field becomes smooth and its rates of diffusion and decay slow down. The boundary condition at the earth‐air interface is handled explicitly by calculating the field in the air from its values at the earth’s surface with an upward continuation based on Laplace’s equation. Boundary conditions in the earth are imposed by using a large, graded grid and setting the values at the sides and bottom to those for a haft‐space. We use the 2-D model to simulate transient electromagnetic (TE) surveys over a thin vertical conductor embedded in a half‐space...

[1]  C. Morawetz,et al.  Solving the Helmholtz Equation for Exterior Problems with Variable Index of Refraction , 1980 .

[2]  A. Bayliss,et al.  Radiation boundary conditions for wave-like equations , 1980 .

[3]  James R. Wait,et al.  Quasi-Static Transient Response of a Conducting Permeable Sphere , 1969 .

[4]  C. R. Deboor,et al.  A practical guide to splines , 1978 .

[5]  T. Lee,et al.  The effect of host rock on transient electromagnetic fields , 1979 .

[6]  Diffusion of electromagnetic fields into the Earth from a line source of current , 1982 .

[7]  Refined finite-difference simulations using local integral forms: Application to telluric fields in two dimensions , 1982 .

[8]  V. Vemuri,et al.  Digital Computer Treatment of Partial Differential Equations , 1981 .

[9]  P. B. Johns,et al.  A simple explicit and unconditionally stable numerical routine for the solution of the diffusion equation , 1977 .

[10]  G. M. Birtwistle.  The explicit solution of the equation of heat conduction , 1968, Comput. J..

[11]  The effect of a conductive overburden on the transient electromagnetic response of a sphere , 1977 .

[12]  James R. Wait.  Transient excitation of the earth by a line source of current , 1971, Proceedings of the IEEE.

[13]  G. W. Hohmann,et al.  An investigation of finite-element modeling for electrical and electromagnetic data in three dimensions , 1981 .

[14]  P. Hood,et al.  Geophysics and geochemistry in the search for metallic ores , 1979 .

[15]  A. Kaufman.  The influence of currents induced in the host rock on electromagnetic response of a spheroid directly beneath a loop , 1981 .

[16]  Y. Lamontagne,et al.  Applications of wide Band Time Domain EM Measurements in Mineral Exploration , 1975 .

[17]  Gerald W. Hohmann,et al.  ELECTROMAGNETIC SCATTERING BY CONDUCTORS IN THE EARTH NEAR A LINE SOURCE OF CURRENT , 1971 .

[18]  P. Weidelt,et al.  Response characteristics of coincident loop transient electromagnetic systems , 1982 .

[19]  S. Orszag,et al.  Approximation of radiation boundary conditions , 1981 .

[20]  R. A. Silverman,et al.  Special functions and their applications , 1966 .

[21]  Ralph A. Stephen,et al.  An implicit finite-difference formulation of the elastic wave equation , 1982 .

[22]  M. L. Oristaglio,et al.  On the approximation of finite loop sources by two‐dimensional line sources , 1984 .

[23]  O. Ladyzhenskaya.  The Boundary Value Problems of Mathematical Physics , 1985 .

[24]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

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

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

[27]  T. Lee.  Transient electromagnetic response of a sphere in a layered medium , 1975 .

[28]  The construction of effective methods for electromagnetic modelling , 1982 .

[29]  B. Spies.  Interpretation and design of time domain EM surveys in areas of conductive overburden , 1980 .