EMDPLER: A F77 program for modeling the EM response of dipolar sources over the non-magnetic layer earth models

A program EMDPLER coded in FORTRAN 77 is presented for modeling the electromagnetic (EM) response of a wide class of dipolar sources (namely, vertical magnetic dipole, horizontal magnetic dipole and horizontal electric dipole) over the layered earth model in induction ([email protected]?50kHz) as well as higher frequency ([email protected][email protected]?1000kHz) regions. The program is based on a semi-analytical numerical approach for evaluating infinite integrals occurring in the expressions of EM field components. It is capable of computing the EM response with or without the displacement current factor for any arbitrary position of source and/or receiver either in air or on the surface of the model due to all the conventional EM dipole sources where other available algorithms fail due to one or other reasons. The program is well documented, easy to understand and amenable to modification for the user's specification. The validity and accuracy of the program is demonstrated by computing the EM response of dipolar sources over homogeneous and multi-layer earth models, and comparing the computed results with those of published results in the literature. Response curves depict their characteristic variations. The computed results match very well with the published results for the induction region ([email protected]?50kHz) and are extension of their characteristic variation in the higher frequency ([email protected][email protected]?1000kHz) region. The matching of computed results with the published results illustrates the validity of the program.

[1]  M. Poddar,et al.  A rectangular loop source of current on a two-layered earth , 1982 .

[2]  A. Sinha,et al.  Electromagnetic fields of oscillating magnetic dipoles placed over a multilayer conducting earth , 1973 .

[3]  James R. Wait,et al.  Mutual Electromagnetic Coupling of Loops Over a Homogeneous Ground , 1955 .

[4]  A. Sommerfeld.  Über die Ausbreitung der Wellen in der drahtlosen Telegraphie , 1909 .

[5]  G. Keller,et al.  Frequency and transient soundings , 1983 .

[6]  A. Sinha,et al.  Electromagnetic fields of an oscillating magnetic dipole over an anisotropic earth , 1968 .

[7]  A. A. Kaufman,et al.  The theory of inductive prospecting , 2001 .

[8]  H. Frank Morrison,et al.  ELECTROMAGNETIC FIELDS ABOUT A LOOP SOURCE OF CURRENT , 1970 .

[9]  John Knight,et al.  Transient electromagnetic calculations using the Gaver-Stehfest inverse Laplace transform method , 1982 .

[10]  M. Poddar.  A rectangular loop source of current on multilayered earth , 1983 .

[11]  Alan D. Chave,et al.  Numerical integration of related Hankel transforms by quadrature and continued fraction expansion , 1983 .

[12]  P. K. Bhattacharya,et al.  Electric dipole over an anisotropic and inhomogeneous earth , 1967 .

[13]  H. Johansen,et al.  FAST HANKEL TRANSFORMS , 1979 .

[14]  D. Oldenburg,et al.  Simultaneous 1 D inversion of loop – loop electromagnetic data for magnetic susceptibility and electrical conductivity , 2003 .

[15]  H. F. Morrison,et al.  QUANTITATIVE INTERPRETATION OF TRANSIENT ELECTROMAGNETIC FIELDS OVER A LAYERED HALF SPACE , 1969 .

[16]  P. Gupta,et al.  Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function , 2006 .

[17]  James R. Wait,et al.  The Magnetic Dipole Over the Horizontally Stratified Earth , 1951 .

[18]  Douglas W. Oldenburg,et al.  Simultaneous 1D inversion of loop loop electromagnetic data for magnetic susceptibility and electrical conductivity , 2003 .

[19]  Walter L. Anderson,et al.  Algorithm 588: Fast Hankel Transforms Using Related and Lagged Convolutions , 1982, TOMS.

[20]  K. Mallick.  Electromagnetic response of a layered transitional earth-infinite cable , 1971 .

[21]  D. Guptasarma,et al.  New digital linear filters for Hankel J0 and J1 transforms , 1997 .

[22]  James R. Wait.  Induction by an oscillating magnetic dipole over a two layer ground , 1959 .

[23]  D. Ghosh,et al.  COMPUTATION OF TYPE CURVES FOR ELECTROMAGNETIC DEPTH SOUNDING WITH A HORIZONTAL TRANSMITTING COIL BY MEANS OF A DIGITAL LINEAR FILTER , 1972 .

[24]  Niels B. Christensen,et al.  Optimized fast Hankel transform filters , 1990 .

[25]  A. Sinha,et al.  Vertical electric dipole over an inhomogeneous and anisotropic earth , 1969 .

[26]  F. Frischknecht,et al.  Fields about an oscillating magnetic dipole over a two-layer earth, and Application to ground and airborne electromagnetic surveys , 1967 .

[27]  G. W. Hohmann,et al.  4. Electromagnetic Theory for Geophysical Applications , 1987 .

[28]  T. Mogi,et al.  Effective skin depth of EM fields due to large circular loop and electric dipole sources , 2003 .

[29]  James R. Wait,et al.  HIGH‐FREQUENCY ELECTROMAGNETIC COUPLING BETWEEN SMALL COPLANAR LOOPS OVER AN INHOMOGENEOUS GROUND , 1972 .

[30]  Walter L. Anderson,et al.  Computer program; numerical integration of related Hankel transforms of orders O and 1 by adaptive digital filtering , 1979 .

[31]  R. K. Verma.  TOPOGRAPHIC EFFECTS ON ELECTROMAGNETIC DEPTH SOUNDING SYSTEMS , 1973 .

[32]  Nagendra Pratap Singh,et al.  EMLCLLER—a program for computing the EM response of a large loop source over a layered earth model☆ , 2003 .