Rapid 2-D/3-D crosshole resistivity imaging using the analytic sensitivity function

We develop a method to obtain 2-D/3-D resistivity images for crosshole dc electric measurements. This imaging procedure consists of two steps: (1) an approximate resistivity image is obtained using a simple projection based on the analytic form of the sensitivity function for the specified electrode array and a modified version of an inversion formula; (2) a final image may be obtained by nonlinear regularized inversion using the approximate resistivity image as the starting model. The imaging procedure is sufficiently flexible to handle different crosshole resistivity applications. Numerical experiments for 2-D and 3-D synthetic models show that the first step can give a satisfactory image of the target between the boreholes or underground when the background medium is uniform. For an inhomogeneous background, the approximate resistivity image is not as good but still conveys a reasonable outline of the resistivity structure. The procedure has no limitation on choice of electrode array (pole-pole, pole-bipole, bipole-pole, bipole-bipole) or field surveying configurations (surface-to-surface, surface-to-borehole, multiborehole crosshole). It may be applied in a routine way for processing crosshole electric data in various applications.

[1]  R. Barker,et al.  Least-squares deconvolution of apparent resistivity pseudosections , 1995 .

[2]  R. Pallàs-Areny,et al.  Subsurface resistivity imaging using two single-step algorithms , 1998 .

[3]  Zhou Bing,et al.  Explicit expressions and numerical calculations for the Fréchet and second derivatives in 2.5D Helmholtz equation inversion , 1999 .

[4]  S. Friedel,et al.  Investigation of a Tertiary maar structure using three-dimensional resistivity imaging , 1999 .

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

[6]  Douglas W. Oldenburg,et al.  Applied geophysical inversion , 1994 .

[7]  Zhou Bing,et al.  A synthetic study on crosshole resistivity imaging using different electrode arrays , 1997 .

[8]  Douglas W. Oldenburg,et al.  Approximate inverse mappings in DC resistivity problems , 1992 .

[9]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[10]  Z. Bing,et al.  Finite element three-dimensional direct current resistivity modelling: accuracy and efficiency considerations , 2001 .

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

[12]  D. Patella.  Introduction to ground surface self‐potential tomography , 1997 .

[13]  Paolo Mauriello,et al.  Resistivity anomaly imaging by probability tomography , 1999 .

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

[15]  J. Daniels.  Three-Dimensional Resistivity and Induced-Polarization Modeling Using Buried Electrodes , 1977 .

[16]  J. Scales,et al.  Smoothing seismic tomograms with alpha-trimmed means , 1988 .

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

[18]  Nariida C. Smith,et al.  Two-Dimensional DC Resistivity Inversion for Dipole-Dipole Data , 1984, IEEE Transactions on Geoscience and Remote Sensing.

[19]  P. Mauriello,et al.  3D geoelectric tomography and archaeological applications 1 , 1998 .

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

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

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

[23]  Jeffrey J. Daniels,et al.  Borehole Resistivity And Electromagnetic Methods Applied to Mineral Exploration , 1984, IEEE Transactions on Geoscience and Remote Sensing.

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

[25]  R. Barker,et al.  Rapid least-squared inversion of apparent resisitivity pseudosections by a quasi-Newton method , 1996 .

[26]  D. Oldenburg,et al.  Inversion of induced polarization data , 1994 .

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

[28]  Biwen Xu,et al.  Archaeological investigation by electrical resistivity tomography: a preliminary study , 1991 .

[29]  W. Daily,et al.  Cross-borehole resistivity tomography , 1991 .