Three-dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example

We have developed a robust and efficient finite difference algorithm for computing the magnetotelluric response of general three-dimensional (3-D) models using the minimum residual relaxation method. The difference equations that we solve are second order in H and are derived from the integral forms of Maxwell's equations on a staggered grid. The boundary H field values are obtained from two-dimensional transverse magnetic mode calculations for the vertical planes in the 3-D model. An incomplete Cholesky decomposition of the diagonal subblocks of the coefficient matrix is used as a preconditioner, and corrections are made to the H fields every few iterations to ensure there are no H divergences in the solution. For a plane wave source field, this algorithm reduces the errors in the H field for simple 3-D models to around the 0.01% level compared to their fully converged values in a modest number of iterations, taking only a few minutes of computation time on our desktop workstation. The E fields can then be determined from discretized versions of the curl of H equations.

[1]  A. Jennings Influence of the Eigenvalue Spectrum on the Convergence Rate of the Conjugate Gradient Method , 1977 .

[2]  Roland W. Freund,et al.  Conjugate Gradient-Type Methods for Linear Systems with Complex Symmetric Coefficient Matrices , 1992, SIAM J. Sci. Comput..

[3]  P. Wannamaker,et al.  Three-dimensional magnetotelluric modeling using difference equations­ Theory and comparisons to integral equation solutions , 1993 .

[4]  D. A. H. Jacobs,et al.  A Generalization of the Conjugate-Gradient Method to Solve Complex Systems , 1986 .

[5]  A. Jennings,et al.  The solution of sparse linear equations by the conjugate gradient method , 1978 .

[6]  Tapan K. Sarkar,et al.  A limited survey of various conjugate gradient methods for solving complex matrix equations arising in electromagnetic wave interactions , 1988 .

[7]  Randall L. Mackie,et al.  Conjugate direction relaxation solutions for 3-d magnetotelluric modeling , 1993 .

[8]  G. W. Hohmann Three-Dimensional Induced Polarization and Electromagnetic Modeling , 1975 .

[9]  O. Axelsson Conjugate gradient type methods for unsymmetric and inconsistent systems of linear equations , 1980 .

[10]  Gerald W. Hohmann,et al.  Electromagnetic modeling of three-dimensional bodies in layered earths using integral equations , 1983 .

[11]  Allen Taflove,et al.  The Finite-Difference Time-Domain Method for Numerical Modeling of Electromagnetic Wave Interactions , 1990 .

[12]  L. Cagniard Basic theory of the magneto-telluric method of geophysical prospecting , 1953 .

[13]  R. Fletcher Conjugate gradient methods for indefinite systems , 1976 .

[14]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[15]  Gerald W. Hohmann,et al.  Integral equation solution for the transient electromagnetic response of a three-dimensional body in a conductive half-space , 1985 .

[16]  D. Kershaw The incomplete Cholesky—conjugate gradient method for the iterative solution of systems of linear equations , 1978 .

[17]  Philip E. Wannamaker,et al.  Advances in three-dimensional magnetotelluric modeling using integral equations , 1991 .

[18]  J. Dennis,et al.  Generalized conjugate directions , 1987 .

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

[20]  Tapan K. Sarkar,et al.  Application of Conjugate Gradient Method to Electromagnetics and Signal Analysis , 1991 .

[21]  F. W. Jones,et al.  The Perturbation of Alternating Geomagnetic Fields by Three-Dimensional Conductivity Inhomogeneities , 1972 .

[22]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[23]  C. M. Swift,et al.  On determining electrical characteristics of the deep layers of the Earth's crust , 1986 .

[24]  D. Rankin,et al.  Three-dimensional modelling in magnetotelluric and magnetic variational sounding , 1977 .