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

Abstract

ABSTRACT

A three‐dimensional (3D) electrical resistivity modelling code is developed to interpret surface and subsurface data. Based on the integral equation, it calculates the charge density caused by conductivity gradients at each interface of the mesh, allowing the estimation of the potential everywhere without the need to interpolate between nodes. Modelling generates a huge matrix, made up of Green's functions, which is stored by using the method of pyramidal compression. The potential is compared with the analytical and the numerical solutions obtained by finite‐difference codes for two models: the two‐layer case and the vertical contact case. The integral method is more accurate around the source point and at the limits of the domain for the potential calculation using a pole‐pole array. A technique is proposed to calculate the sensitivity (Jacobian) and Hessian matrices in 3D. The sensitivity is based on the derivative with respect to the block conductivity of the potential computed using the integral equation; it is only necessary to compute the electrical field at the source location. A direct extension of this technique allows the determination of the second derivatives. The technique is compared with the analytical solutions and with the calculation of the sensitivity according to the method using the inner product of the current densities calculated at the source and receiver points. Results are very accurate when the Green's function that includes the source image is used. The calculation of the three components of the electric field on the interfaces of the mesh is carried out simultaneously and quickly, using matrix compression.

Loading

Article metrics loading...

/content/journals/10.1111/j.1365-2478.2005.00484.x
2005-06-24
2024-04-26
Loading full text...

Full text loading...

References

  1. AlfanoL.1959. Introduction to the interpretation of resistivity measurements for complicated structural conditions. Geophysical Prospecting7, 311–366.
    [Google Scholar]
  2. AlfanoL.1960. The influence of surface formations on the apparent resistivity values in electrical prospecting, part 1. Geophysical Prospecting8, 576–606.
    [Google Scholar]
  3. AlfanoL.1961. The influence of surface formations on the apparent resistivity values in electrical prospecting, part 2. Geophysical Prospecting9, 213–241.
    [Google Scholar]
  4. BeasleyC.W. and WardS.H.1986. Three‐dimensional mise‐à‐la‐masse modeling applied to mapping fracture zones. Geophysics51, 98–113.
    [Google Scholar]
  5. BhattacharyaP.K. and PatraH.P.1968. Direct Current Geoelectric Sounding: Principles and Interpretations . Volume 9 of Methods in Geochemistry and Geophysics . Elsevier Science Publishing Co.
    [Google Scholar]
  6. BingZ. and GreenhalghS.A.1999. Explicit expressions and numerical calculations for Frechet and second derivatives in 2.5D Helmholtz equation inversion. Geophysical Prospecting47, 443–468.
    [Google Scholar]
  7. BingZ. and GreenhalghS.A.2001. Finite element three‐dimensional direct current resistivity modelling: accuracy and efficiency considerations. Geophysical Journal International145, 679–688.
    [Google Scholar]
  8. BjörckA.1996. Numerical Methods for Least Squares Problems . SIAM Publications.
    [Google Scholar]
  9. BoernerD.E. and WestG.F.1989. Frechet derivatives and single scattering theory. Geophysical Journal International98, 385–390.
    [Google Scholar]
  10. Brewitt‐TaylorC.R. and WeaverJ.T.1976. On the finite difference solution of two‐dimensional induction problems. Geophysical Journal of the Royal Astronomical Society47, 375–396.
    [Google Scholar]
  11. BurtP.J. and AdelsonE.H.1983. The Laplacian pyramid as a compact image code. IEEE Transactions on CommunicationsCOM‐31, 337–345.
    [Google Scholar]
  12. CoggonJ.H.1971. Electromagnetic and electrical modeling by finite element method. Geophysics36, 132–155.
    [Google Scholar]
  13. DeyA. and MorrisonH.F.1979a. Resistivity modeling for arbitrarily shaped three‐dimensional structures. Geophysics44, 753–780.
    [Google Scholar]
  14. DeyA. and MorrisonH.F.1979b. Resistivity modelling for arbitrarily shaped two‐dimensional structures. Geophysical Prospecting27, 106–136.
    [Google Scholar]
  15. FreundR. and NachtigalN.1991a. An implementation of the QMR method based on coupled two‐term recurrences. SIAM Journal on Scientific and Statistical Computation15, 313–337.
    [Google Scholar]
  16. FreundR. and NachtigalN.1991b. QMR: A quasi‐minimal residual method for non‐Hermitian linear systems. Numerical Mathematics60, 315–339.
    [Google Scholar]
  17. GeselowitzD.B.1971. An application of electrocardiographic lead theory to impedance plethysmography. IEEE Transactions on Bio-Medical EngineeringBME‐18, 38–41.
    [Google Scholar]
  18. GolubG.H. and Van LoanC.F.1996. Matrix Computations , 3rd edn.Johns Hopkins University Press.
    [Google Scholar]
  19. HestenesM. and StiefelE.1952. Methods of conjugate gradients for solving linear systems. Journal of Research of the National Bureau of Standards49, 409–436.
    [Google Scholar]
  20. HohmannG.W.1975. Three‐dimensional induced polarization and electromagnetic modeling. Geophysics40, 309–324.
    [Google Scholar]
  21. HvozdaraM. and KaikkonenP.1998. An integral equations solution of the forward D.C. geoelectric problem for a 3‐D body of inhomogeneous conductivity buried in a halfspace. Journal of Applied Geophysics39, 95–107.
    [Google Scholar]
  22. LaBrecqueD.J., MilettoM., DailyW., RamirezA. and OwenE.1996. The effects of noise on Occam's inversion of resistivity tomography data. Geophysics61, 538–548.
    [Google Scholar]
  23. LaBrecqueD.J., MorelliG., DailyW., RamirezA. and LundegardP.1999. Occam's inversion of 3‐D electrical resistivity tomography. In: Three‐Dimensional Electromagnetics (eds M.Oristaglio and B.Spies ), pp. 575–590. Society of Exploration Geophysicists , Tulsa , OK .
    [Google Scholar]
  24. LanczosC.1961. Linear Differential Operators . Van Nostrand ( Dover , 1997).
    [Google Scholar]
  25. LesurV., CuerM. and StraubA.1999. 2‐D and 3‐D interpretation of electrical tomography measurements, Part 1: The forward problem. Geophysics64, 386–395.
    [Google Scholar]
  26. LiY. and OldenburgD.W.1991. Aspects of charge accumulation in DC resistivity experiments. Geophysical Prospecting39, 803–826.
    [Google Scholar]
  27. LiY. and OldenburgD.1997. Fast inversion of large scale magnetic data using wavelets. 67th SEG meeting, Dallas , USA , Expanded Abstracts, 490–493.
  28. LosanoV. and LagetB.1996. Fractional pyramids for color image segmentation. In: Proceedings of the IEEE Southwest Symposium on Image Analysis and Interpretation , San Antonio , Texas , USA , pp. 13–17.
  29. LowryT., AllenM.B. and ShiveP.N.1989. Singularity removal: A refinement of resistivity modeling techniques. Geophysics54, 766–774.
    [Google Scholar]
  30. McGillivrayP. and OldenburgD.W.1990. Methods for calculating Frechet derivates and sensitivities for the non‐linear inverse problem: a comparative study. Geophysical Prospecting38, 499–524.
    [Google Scholar]
  31. MuftiI.1978. A practical approach to finite‐difference resistivity modeling. Geophysics43, 930–942.
    [Google Scholar]
  32. MundryE.1984. Geoelectrical model calculations for two‐dimensional resistivity distributions. Geophysical Prospecting32, 124–131.
    [Google Scholar]
  33. MuraiT. and KagawaY.1985. Electrical impedance computed tomography based on finite element model. IEEE Transactions on Biomedical EngineeringBME‐32, 177–184.
    [Google Scholar]
  34. NocedalJ. and WrightS.J.2000. Numerical Optimization . Springer‐Verlag, Inc.
    [Google Scholar]
  35. ParkS.K. and VanG.P.1991. Inversion of pole‐pole data for 3‐D resistivity structure beneath arrays of electrodes. Geophysics56, 951–960.
    [Google Scholar]
  36. PissanetzkyS.1984. Sparse Matrix Technology . Academic Press, Inc.
    [Google Scholar]
  37. PoirmeurC. and VasseurG.1988. Three‐dimensional modeling of a hole‐to‐hole electrical method: Application to the interpretation of a field survey. Geophysics53, 402–414.
    [Google Scholar]
  38. PortniaguineO.1999. Image focusing and data compression in the solution of geophysical inverse problems . PhD thesis, University of Utah .
  39. PrattD.A.1972. The surface integral approach to the solution of the 3D resistivity problem. Bulletin of the Australian Society for Exploration Geophysics3, 33–50.
    [Google Scholar]
  40. PressW.H., TeukolskyS.A., VetterlingW.T. and FlanneryB.P.1992. Numerical Recipes in C: The Art of Scientific Computing , 2nd edn.Cambridge University Press.
    [Google Scholar]
  41. PridmoreD.F., HohmannG.W., WardS.H. and SillW.R.1981. An investigation of finite‐element modeling for electrical and electromagnetic data in three dimensions. Geophysics46, 1009–1024.
    [Google Scholar]
  42. RaicheA.P.1974. An integral equation approach to three‐dimensional modelling. Geophysical Journal of the Royal Astronomical Society36, 363–376.
    [Google Scholar]
  43. RoyA. and ApparaoA.1971. Depth of investigation in direct current methods. Geophysics36, 943–959.
    [Google Scholar]
  44. SasakiY.1992. Resolution of resistivity tomography inferred from numerical simulation. Geophysical Prospecting40, 453–463.
    [Google Scholar]
  45. SasakiY.1994. 3‐D resistivity inversion using the finite‐element method. Geophysics59, 1839–1848.
    [Google Scholar]
  46. SchulzR.1985. The method of integral equation in direct current resistivity method and its accuracy. Journal of Geophysics56, 192–200.
    [Google Scholar]
  47. ScribaH.1981. Computation of the electric potential in three‐dimensional structures. Geophysical Prospecting29, 790–802.
    [Google Scholar]
  48. SeigelH.O.1959. Mathematical formulation and type curves for induced polarization. Geophysics24, 547–565.
    [Google Scholar]
  49. SleijpenG. and FokkemaD.1993. BICGSTAB(l) for linear equations involving unsymmetric matrices with complex spectrum. Electronic Transactions on Numerical Analysis1, 11–32.
    [Google Scholar]
  50. SnyderD.D.1976. A method for modeling the resistivity and IP response of two‐dimensional bodies. Geophysics41, 997–1015.
    [Google Scholar]
  51. SonneveldP.1989. CGS, a fast Lanczos‐type solver for nonsymmetric linear systems. SIAM Journal on Scientific and Statistical Computation10, 36–52.
    [Google Scholar]
  52. SpiegelR.J., SturdivantV.R. and OwenT.E.1980. Modeling resistivity anomalies from localized voids under irregular terrain. Geophysics45, 1164–1183.
    [Google Scholar]
  53. SpitzerK.1995. A 3‐D finite‐difference algorithm for dc resistivity modelling using conjugate gradient methods. Geophysical Journal International123, 903–914.
    [Google Scholar]
  54. SpitzerK.1998. The three‐dimensional DC sensitivity for surface and subsurface sources. Geophysical Journal International134, 736–746.
    [Google Scholar]
  55. SpitzerK., ChouteauM. and BoulangerO.1999. Grid‐independent electrode positioning for 3D DC resistivity and IP forward modelling. In: Symposium 3DEM2 , Salt Lake City , Utah , pp. 189–192.
    [Google Scholar]
  56. 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. Society of Exploration Geophysicists , Tulsa , OK .
    [Google Scholar]
  57. TelfordW.M., GeldartL.P. and SheriffR.E.1990. Applied Geophysics , 2nd edn.Cambridge University Press.
    [Google Scholar]
  58. TingS. and HohmannG.1981. Integral equation modeling of three‐dimensional magnetotelluric response. Geophysics46, 182–197.
    [Google Scholar]
  59. Van Der VorstH.1992. Bi‐CGSTAB: A fast and smoothly converging variant of Bi‐CG for the solution of nonsymmetric linear systems. SIAM Journal on Scientific and Statistical Computation13, 631–644.
    [Google Scholar]
  60. WangT. and MezzatestaA.2001. 3D DC resistivity modeling using the sweeping‐seed conjugate gradient method. Geophysics66, 441–447.
    [Google Scholar]
  61. WannamakerP., HohmannG. and SanFilipoW.1984. Electromagnetic modeling of three‐dimensional bodies in layered earths using integral equations. Geophysics49, 60–74.
    [Google Scholar]
  62. WurmstichB. and MorganF.D.1994. Modeling of streaming potential responses caused by oil well pumping. Geophysics59, 46–56.
    [Google Scholar]
  63. XiongZ.1989. Short note: Electromagnetic fields of electric dipoles embedded in a stratified anisotropic earth. Geophysics54, 1643–1646.
    [Google Scholar]
  64. XiongZ.1992a. Electromagnetic modeling of 3‐D structures by the method of system iteration using integral equations. Geophysics57, 1556–1561.
    [Google Scholar]
  65. XiongZ.1992b. Short note: Symmetry properties of the scattering matrix in 3‐D electromagnetic modeling using the integral equation method. Geophysics57, 1199–1202.
    [Google Scholar]
  66. ZhangJ., MackieR.L. and MaddenT.R.1995. Three‐dimensional resistivity forward modeling and inversion using conjugate gradients. Geophysics60, 1313–1325.
    [Google Scholar]
  67. ZhaoS. and YedlinM.J.1996. Some refinements on the finite‐difference method for 3‐D resistivity modeling. Geophysics61, 1301–1307.
    [Google Scholar]
  68. ZhdanovM.S. and KellerG.V.1994. The Geoelectrical Methods in Geophysical Exploration . Volume 31 of Methods in Geochemistry and Geophysics . Elsevier Science Publishing Co.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/j.1365-2478.2005.00484.x
Loading
/content/journals/10.1111/j.1365-2478.2005.00484.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