Crosswell electromagnetic inversion using integral and differential equations

The crosswell electromagnetic (EM) inverse problem is solved with an integral‐equation (IE) formulation using successive Born approximations in the frequency domain. Because the inverse problem is nonlinear, the predicted fields and Green’s functions are continually updated. Updating the fields and Green’s functions relates small changes in the predicted data to small changes in the model parameters through Frechet kernels. These fields and Green functions are calculated with an efficient 3-D finite‐difference solver. Since the resistivity is invariant along strike, the 3-D fields are integrated along strike so the 2-D kernels can be assembled. At the early stages of the inversion, smoothing of the electrical conductivity stabilizes the inverse solution when it is far from convergence. As the solution converges, this smoothing is relaxed and more effort is made to reduce the data misfit. Bounds on the conductivity are included in the solution to eliminate unrealistic estimates. The robustness of the inver...

[1]  T. Dobecki,et al.  Geotechnical and groundwater geophysics , 1985 .

[2]  Carlos Torres-Verdín,et al.  An Approach to Nonlinear Inversion With Applications to Cross-Well EM Tomography , 1993 .

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

[4]  Abelardo Ramirez,et al.  Evaluation of electromagnetic tomography to map in situ water in heated welded tuff , 1989 .

[5]  G. Backus,et al.  Uniqueness in the inversion of inaccurate gross Earth data , 1970, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

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

[7]  M. Heroux,et al.  A parallel preconditioned conjugate gradient package for solving sparse linear systems on a Cray Y-MP , 1991 .

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

[9]  T. Sarkar On the Application of the Generalized BiConjugate Gradient Method , 1987 .

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

[11]  D. Alumbaugh,et al.  Iterative Electromagnetic Born Inversion Applied to Earth Conductivity Imaging , 1993 .

[12]  G. Backus,et al.  The Resolving Power of Gross Earth Data , 1968 .

[13]  J. T. Smith Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator , 1996 .

[14]  M. Toksöz,et al.  Simultaneous reconstruction of permittivity and conductivity for crosshole geometries , 1990 .

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

[16]  Weng Cho Chew,et al.  Simultaneous reconstruction of permittivity and conductivity profiles in a radially inhomogeneous slab , 1986 .

[17]  G. E. Archie The electrical resistivity log as an aid in determining some reservoir characteristics , 1942 .

[18]  G. Xie,et al.  A new approach to imaging with low‐frequency electromagnetic fields , 1993 .

[19]  Stephen K. Park,et al.  Inversion of pole-pole data for 3-D resistivity structure beneath arrays of electrodes , 1991 .

[20]  G. Backus,et al.  Numerical Applications of a Formalism for Geophysical Inverse Problems , 1967 .

[21]  H. F. Morrison,et al.  Crosswell electromagnetic tomography: System design considerations and field results , 1995 .

[22]  G. W. Hohmann,et al.  Three-dimensional Electromagnetic Inversion , 1988 .

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