Theoretical development of the differential scattering decomposition for the 3D resistivity experiment

In any numerical solution of the DC resistivity experiment, care must be taken to deal with strong heterogeneity of electrical conductivity. In order to examine the importance of conductivity contrasts, we develop a scattering decomposition of the DC resistivity equation in the sparse differential domain as opposed to the traditional dense integral formulation of scattering-type equations. We remove the singularity in the differential scattered series via separation of primary and secondary conductivity, thereby avoiding the need to address the singularity in a Green's function. The differential scattering series is observed to diverge for large conductivity contrasts and to converge for small contrasts. We derive a convergence criterion, in terms of matrix norms for the weak-form finite-volume equations, that accounts for both the magnitude and distribution of heterogeneity of electrical conductivity. We demonstrate the relationship between the differential scattering series and the Frechet derivative of the electrical potential with respect to electrical conductivity, and we show how the development may be applied to the inverse problem. For linearization associated with the Frechet derivative to be valid, the perturbation in electrical conductivity must be small as defined by the convergence of the scattered series. The differential scattering formulation also provides an efficient tool for gaining insight into charge accumulation across contrasts in electrical conductivity, and we present a derivation that equates accumulated surface charge density to the source of scattered potential.

[1]  Yuguo Li,et al.  Three‐dimensional DC resistivity forward modelling using finite elements in comparison with finite‐difference solutions , 2002 .

[2]  P. McGillivray,et al.  Forward modeling and inversion of DC resistivity and MMR data , 1992 .

[3]  Xiaoping Wu,et al.  Computations of secondary potential for 3D DC resistivity modelling using an incomplete Choleski conjugate‐gradient method , 2003 .

[4]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[5]  Matthew J. Yedlin,et al.  Some refinements on the finite-difference method for 3-D dc resistivity modeling , 1996 .

[6]  A. Dey,et al.  Resistivity modeling for arbitrarily shaped three-dimensional structures , 1979 .

[7]  Gene H. Golub,et al.  Matrix computations , 1983 .

[8]  Douglas W. Oldenburg,et al.  3-D inversion of induced polarization data3-D Inversion of IP Data , 2000 .

[9]  D. Oldenburg,et al.  3-D inversion of induced polarization data , 2001 .

[10]  K. Spitzer A 3-D FINITE-DIFFERENCE ALGORITHM FOR DC RESISTIVITY MODELLING USING CONJUGATE GRADIENT METHODS , 1995 .

[11]  Klaus Spitzer,et al.  The three‐dimensional DC sensitivity for surface and subsurface sources , 1998 .

[12]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[13]  D. Oldenburg,et al.  ASPECTS OF CHARGE ACCUMULATION IN d. c. RESISTIVITY EXPERIMENTS1 , 1991 .

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

[15]  D. D. Snyder A method for modeling the resistivity and IP response of two-dimensional bodies , 1976 .

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

[17]  Peter N. Shive,et al.  Singularity removal: A refinement of resistivity modeling techniques , 1989 .

[18]  M. Ali Ak,et al.  ‘European’ Association of Geoscientists & Engineers , 1997 .

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