Three‐dimensional massively parallel electromagnetic inversion—I. Theory

Summary An iterative solution to the non-linear 3-D electromagnetic inverse problem is obtained by successive linearized model updates using the method of conjugate gradients. Full wave equation modelling for controlled sources is employed to compute model sensitivities and predicted data in the frequency domain with an efficient 3-D finite-difference algorithm. Necessity dictates that the inverse be underdetermined, since realistic reconstructions require the solution for tens of thousands of parameters. In addition, large-scale 3-D forward modelling is required and this can easily involve the solution of over several million electric field unknowns per solve. A massively parallel computing platform has therefore been utilized to obtain reasonable execution times, and results are given for the 1840-node Intel Paragon. The solution is demonstrated with a synthetic example with added Gaussian noise, where the data were produced from an integral equation forward-modelling code, and is different from the finite difference code embedded in the inversion algorithm

[1]  T. Habashy,et al.  Rapid 2.5‐dimensional forward modeling and inversion via a new nonlinear scattering approximation , 1994 .

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

[3]  Jie Zhang,et al.  3-D resistivity forward modeling and inversion using conjugate gradients , 1995 .

[4]  G. Newman,et al.  Three-dimensional massively parallel electromagnetic inversion—II. Analysis of a crosswell electromagnetic experiment , 1997 .

[5]  Stephen K. Park,et al.  Three-dimensional magnetotelluric modelling and inversion , 1983 .

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

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

[8]  Anthony Skjellum,et al.  Using MPI - portable parallel programming with the message-parsing interface , 1994 .

[9]  W. Menke.  Geophysical data analysis : discrete inverse theory , 1984 .

[10]  D. Oldenburg,et al.  Inversion of electromagnetic data: An overview of new techniques , 1990 .

[11]  D. Oldenburg,et al.  METHODS FOR CALCULATING FRÉCHET DERIVATIVES AND SENSITIVITIES FOR THE NON‐LINEAR INVERSE PROBLEM: A COMPARATIVE STUDY1 , 1990 .

[12]  D. Oldenburg,et al.  Approximate sensitivities for the electromagnetic inverse problem , 1996 .

[13]  R. L. Mackie,et al.  Three-dimensional magnetotelluric modelling and inversion , 1989, Proc. IEEE.

[14]  R. Mackie,et al.  Three-dimensional magnetotelluric inversion using conjugate gradients , 1993 .

[15]  T. Habashy,et al.  A two-step linear inversion of two-dimensional electrical conductivity , 1995 .