2-D and 3-D interpretation of electrical tomography measurements, Part 2: The inverse problem

The interpretation of borehole‐to‐borehole electrical measurements requires solving an inverse problem for a given class of model geometries. The usual approach to an inverse problem includes a model dependent task (i.e., forward modeling) and a generic task (i.e., an optimization scheme). We have developed an optimization algorithm using a nonlinear inversion technique. This algorithm allows recovery of a possible resistivity distribution in an investigated zone between two boreholes or in the vicinity of them. This resistivity distribution is defined as a set of 2-D or 3-D volumes of constant resistivity. The inversion procedure minimizes a least‐squares term plus a damping term. This latter term seeks to minimize the roughness of the solution. An improved form of this smoothness term may enhance the spatial resolution of the resistivity image, assuming that the resistivity contrast is known a priori. This reconstruction algorithm has been tested for both 2-D and 3-D geometries. These inversion tests we...

[1]  Kenneth Levenberg.  A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[2]  Robert G. Ellis,et al.  The pole-pole 3-D Dc-resistivity inverse problem: a conjugategradient approach , 1994 .

[3]  G. Stewart.  The economical storage of plane rotations , 1976 .

[4]  Angular singularities of elliptic problems , 1972 .

[5]  R. Kohn,et al.  Relaxation of a variational method for impedance computed tomography , 1987 .

[6]  A. George.  Nested Dissection of a Regular Finite Element Mesh , 1973 .

[7]  Douglas W. Oldenburg,et al.  Inversion of 3-D DC resistivity data using an approximate inverse mapping , 1994 .

[8]  A. Nachman,et al.  Global uniqueness for a two-dimensional inverse boundary value problem , 1996 .

[9]  Ziqi Sun.  The inverse conductivity problem in two dimensions , 1990 .

[10]  David Isaacson,et al.  Layer stripping: a direct numerical method for impedance imaging , 1991 .

[11]  W. Daily,et al.  The effects of noise on Occam's inversion of resistivity tomography data , 1996 .

[12]  D. Oldenburg,et al.  Generalized subspace methods for large-scale inverse problems , 1993 .

[13]  Vladimir Druskin,et al.  The Lanczos optimization of a splitting-up method to solve homogeneous evolutionary equations , 1992 .

[14]  Ivo Babuska,et al.  The p and h-p Versions of the Finite Element Method, Basic Principles and Properties , 1994, SIAM Rev..

[15]  J. Sylvester,et al.  A global uniqueness theorem for an inverse boundary value problem , 1987 .

[16]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[17]  Hiromasa Shima,et al.  Two-dimensional automatic resistivity inversion technique using alpha centers , 1990 .

[18]  Andrew Binley,et al.  ERT monitoring of environmental remediation processes , 1996 .

[19]  Claude Lemaréchal,et al.  Some numerical experiments with variable-storage quasi-Newton algorithms , 1989, Math. Program..

[20]  Zhiqiang Cai,et al.  The finite volume element method for diffusion equations on general triangulations , 1991 .

[21]  Douglas LaBrecque,et al.  Monitoring an underground steam injection process using electrical resistance tomography , 1993 .

[22]  Yutaka Sasaki,et al.  3-D resistivity inversion using the finite-element method , 1994 .

[23]  G. Vasseur,et al.  Three‐dimensional modeling of a hole‐to‐hole electrical method: Application to the interpretation of a field survey , 1988 .

[24]  D. Dobson,et al.  An image-enhancement technique for electrical impedance tomography , 1994 .

[25]  3-D Inversion In Subsurface Electrical Surveying—I. Theory , 1994 .

[26]  Robert V. Kohn,et al.  Numerical implementation of a variational method for electrical impedance tomography , 1990 .

[27]  Albert G. Buckley,et al.  Remark on algorithm 630 , 1989, TOMS.