1887
Volume 51, Issue 6
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

An accurate and efficient 3D finite‐difference (FD) forward algorithm for DC resistivity modelling is developed. In general, the most time‐consuming part of FD calculation is to solve large sets of linear equations: =, where is a large sparse band symmetric matrix. The direct method using complete Choleski decomposition is quite slow and requires much more computer storage. We have introduced a row‐indexed sparse storage mode to store the coefficient matrix and an incomplete Choleski conjugate‐gradient (ICCG) method to solve the large linear systems. By taking advantage of the matrix symmetry and sparsity, the ICCG method converges much more quickly and requires much less computer storage. It takes approximately 15 s on a 533 MHz Pentium computer for a grid with 46 020 nodes, which is approximately 700 times faster than the direct method and 2.5 times faster than the symmetric successive over‐relaxation (SSOR) conjugate‐gradient method. Compared with 3D finite‐element resistivity modelling with the improved ICCG solver, our algorithm is more efficient in terms of number of iterations and computer time. In addition, we solve for the secondary potential in 3D DC resistivity modelling by a simple manipulation of the FD equations. Two numerical examples of a two‐layered model and a vertical contact show that the method can achieve much higher accuracy than solving for the total potential directly with the same grid nodes. In addition, a 3D cubic body is simulated, for which the dipole–dipole apparent resistivities agree well with the results obtained with the finite‐element and integral‐equation methods. In conclusion, the combination of several techniques provides a rapid and accurate 3D FD forward modelling method which is fundamental to 3D resistivity inversion.

Loading

Article metrics loading...

/content/journals/10.1046/j.1365-2478.2003.00392.x
2003-10-23
2024-04-26
Loading full text...

Full text loading...

References

  1. AjizM.A. and JenningsA. 1984. A robust incomplete Choleski conjugate gradient algorithm. International Journal of Numerical Methods in Engineering20, 949–960.
    [Google Scholar]
  2. BentleyJ.1986. Programming Pearls . Addison‐Wesley Publishing Co., Reading , MA .
    [Google Scholar]
  3. CoggonJ.H.1971. Electromagnetic and electrical modelling by the finite element method. Geophysics36, 132–155.
    [Google Scholar]
  4. DanielsJ.J.1977. Three‐dimensional resistivity and induced polarization modelling using buried electrodes. Geophysics42, 1006–1019.
    [Google Scholar]
  5. DeyA. and MorrisonH.F.1979. Resistivity modelling for arbitrarily shaped three‐dimensional structures. Geophysics44, 753–780.
    [Google Scholar]
  6. DieterK., PatersonN.R. and GrantF.S.1969. IP and resistivity type curves for three‐dimensional bodies. Geophysics34, 615–632.
    [Google Scholar]
  7. FoxR.C., HohmannG.W., KillpackT.J. and RijoL.1980. Topographic effects in resistivity and induced‐polarization surveys. Geophysics45, 75–93.
    [Google Scholar]
  8. FuLiangkui1983. Lectures on Geoelectrical Prospecting , pp. 50–56. Geology Press, Beijing .
    [Google Scholar]
  9. HestenseM. and StiefelE.1952. Methods of conjugate gradients for solving linear systems. Journal of Research of the National Bureau of Standards49, 409–436.
    [Google Scholar]
  10. HohmannG.W.1975. Three‐dimensional induced polarization and electromagnetic modelling. Geophysics40, 309–324.
    [Google Scholar]
  11. HolcombeH.T. and JiracekG.R.1984. Three‐dimensional terrain correction in resistivity surveys. Geophysics49, 439–452.
    [Google Scholar]
  12. KershawD.S.1978. The incomplete Choleski conjugate gradient method for the iterative solution of systems of linear equations. Journal of Computational Physics26, 43–65.
    [Google Scholar]
  13. LeeT.1975. An integral equation and its solution for some two‐ and three‐dimensional problems in resistivity and induced polarization. Geophysical Journal of the Royal Astronomical Society42, 81–95.
    [Google Scholar]
  14. LokeM.H. and BarkerR.D.1996. Practical techniques for 3D resistivity surveys and data inversion. Geophysical Prospecting44, 499–523.
    [Google Scholar]
  15. ManteuffelT.A.1980. An incomplete factorization technique for positive definite linear systems. Mathematics of Computation34, 473–497.
    [Google Scholar]
  16. MeijerinkJ.A. and Van Der VorstH.A.1977. An iterative solution method for linear systems of which the coefficient matrix is a symmetric matrix. Mathematics of Computation31, 148–162.
    [Google Scholar]
  17. NewmanG.A. and AlumbaughD.L.2000. Three‐dimensional magnetotelluric inversion using non‐linear conjugate gradients. Geophysical Journal International140, 410–424.
    [Google Scholar]
  18. OkabeM.1981. Boundary element method for the arbitrary inhomogeneities problem in electrical prospecting. Geophysical Prospecting29, 39–59.
    [Google Scholar]
  19. O'NeillD.J. et al . 1984. A digital linear filter for resistivity sounding with a generalized electrode array. Geophysical Prospecting32, 105–123.
    [Google Scholar]
  20. OppligerG.L.1984. Three‐dimensional terrain corrections for mise‐a‐la‐masse and magnetometric resistivity surveys. Geophysics49, 1718–1729.
    [Google Scholar]
  21. PapadrakakisM. and DracopoulosM.C.1991. Improving the efficiency of incomplete Choleski preconditionings. Communications in Applied Numerical Methods7, 603–612.
    [Google Scholar]
  22. 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]
  23. PridmoreD.F., HohmannG.W., WardS.H. and SillW.R.1981. An investigation of finite element modelling for electrical and electromagnetic data in three dimensions. Geophysics46, 1009–1024.
    [Google Scholar]
  24. ScribeH.1981. Computations of the electrical potential in the three‐dimensional structure. Geophysical Prospecting29, 790–802.
    [Google Scholar]
  25. SmithJ.T.1996. Conservative modelling of 3D electromagnetic fields, Part II: biconjugate gradient solution as an accelerator. Geophysics61, 1319–1324.
    [Google Scholar]
  26. SpitzerK.1995. A 3D finite difference algorithm for DC resistivity modelling using conjugate gradient methods. Geophysical Journal International123, 903–914.
    [Google Scholar]
  27. VargaR.S.1960. Factorization and normalized iterative methods. In: Boundary Problems in Differential Equations (ed. R.E.Langer ), pp. 121–142. University of Wisconsin Press, Madison , WI .
    [Google Scholar]
  28. XuS.Z., GuoZ.C. and ZhaoS.K.1988. An integral formulation for three‐dimensional terrain modelling for resistivity surveys. Geophysics53, 546–552.
    [Google Scholar]
  29. ZhangJ., MackieR.L. and MaddenT.R.1995. Three‐dimensional resistivity forward modelling and inversion using conjugate gradients. Geophysics60, 1313–1325.
    [Google Scholar]
  30. ZhaoS.K. and YedlinM.J.1996. Some refinements on the finite‐difference method for 3D DC resistivity modelling. Geophysics61, 1301–1307.
    [Google Scholar]
  31. 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.1046/j.1365-2478.2003.00392.x
Loading
/content/journals/10.1046/j.1365-2478.2003.00392.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