1887
Volume 54, Issue 4
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

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 Fréchet 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 Fréchet 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.

Loading

Article metrics loading...

/content/journals/10.1111/j.1365-2478.2006.00547.x
2006-06-29
2024-04-25
Loading full text...

Full text loading...

References

  1. DCIP3D
    DCIP3D2004. A Program Library for Forward Modelling and Inversion of DC Resistivity and Induced Polarization Data Over 3D Structures, Version 2.0 . Developed under the consortium research project: Inversion of 3D resistivity and induced polarization data. UBC Geophysical Inversion Facility, University of British Columbia .
    [Google Scholar]
  2. DeyA. and MorrisonH.F.1979. Resistivity modeling for arbitrarily shaped three‐dimensional structures. Geophysics44, 753–780.
    [Google Scholar]
  3. GolubG.H. and Van LoanC.F.1996. Matrix Computations , 3rd edn.John Hopkins University Press, ISBN 0801854148.
    [Google Scholar]
  4. HohmannG.W. and RaicheA.P.1988. Inversion of controlled‐source electromagnetic data. In: Electromagnetic Methods in Applied Geophysics , Vol. 1 (ed. M.N.Nabighian ), pp. 469–503. Investigations in Geophysics 3. Society of Exploration Geophysicists , ISBN 1560800690.
    [Google Scholar]
  5. LiY. and OldenburgD.W.1991. Aspects of charge accumulation in d.c. resistivity experiments. Geophysical Prospecting39, 803–826.
    [Google Scholar]
  6. LiY. and OldenburgD.W.2000. 3‐D inversion of induced polarization data. Geophysics65, 1931–1945.
    [Google Scholar]
  7. LiY. and SpitzerK.2002. Three‐dimensional DC resistivity forward modelling using finite elements in comparison with finite‐difference solutions. Geophysical Journal International151, 924–934.
    [Google Scholar]
  8. LowryT., AllenM.B. and ShiveP.N.1989. Singularity removal: A refinement of resistivity modeling techniques. Geophysics54, 766–774.
    [Google Scholar]
  9. McGillivrayP.R.1992. Forward modeling and inversion of DC resistivity and MMR data . PhD thesis, University of British Columbia .
    [Google Scholar]
  10. McGillivrayP.R. and OldenburgD.W.1990. Methods for calculating Fréchet derivatives and sensitivities for the non‐linear inverse problem: A comparative study. Geophysical Prospecting38, 499–524.
    [Google Scholar]
  11. PressW.H, TeukolskyA.A., VetterlingW.T. and FlanneryB.P.1992. Numerical Recipes in Fortran 77, The Art of Scientific Computing , 2nd edn.Cambridge University Press, ISBN 052143064X.
    [Google Scholar]
  12. SnyderD.D.1976. A method for modeling the resistivity and IP response of two‐dimensional bodies. Geophysics41, 997–1015.
    [Google Scholar]
  13. SpitzerK.1995. A 3‐D finite difference algorithm for DC resistivity modelling using conjugate gradient methods. Geophysical Journal International123, 903–914.
    [Google Scholar]
  14. SpitzerK.1998. The three‐dimensional DC sensitivity for surface and subsurface sources. Geophysical Journal International134, 736–746.
    [Google Scholar]
  15. SpitzerK., ChouteauM. and BoulangerO.1999. Grid‐independent electrode positioning for 3D DC resistivity and IP forward modeling. 2nd International Symposium on Three Dimensional Electromagnetics, Salt Lake City , UT , USA , Extended Abstracts, 189–192.
  16. SpitzerK. and WurmstichB.1999. Speed and accuracy in 3‐D resistivity modeling. In: Three‐Dimensional Electromagnetics (eds  M.Oristaglio and B.Spies ), pp. 161–176. Geophysical Developments 7. Society of Exploration Geophysicists , ISBN 1560800798.
    [Google Scholar]
  17. TelfordW.M., GeldartL.P. and SheriffR.E.1990. Applied Geophysics , 2nd edn.Cambridge University Press, ISBN 0521339383.
    [Google Scholar]
  18. WuX., XiaoY., QiC. and WangT.2003. Computations of secondary potential for 3D resistivity modelling using an incomplete Choleski conjugate‐gradient method. Geophysical Prospecting51, 567–577.
    [Google Scholar]
  19. ZhaoS. and YedlinM.J.1996. Some refinements on the finite‐difference method for 3‐D dc resistivity modeling. Geophysics61, 1301–1307.
    [Google Scholar]
  20. ZhouB. and GreenhalghS.A.2001. Finite element three‐dimensional direct current resistivity modelling: Accuracy and efficiency considerations. Geophysical Journal International145, 679–688.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/j.1365-2478.2006.00547.x
Loading
/content/journals/10.1111/j.1365-2478.2006.00547.x
Loading

Data & Media loading...

  • Article Type: Research Article

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error