2.5-D Resistivity Tomography using boundary integral equations

Abstract

DC Resistivity Tomography is a non-linear inversion problem. So far there are mainly two kinds of inversion methods, based on the finite-element method and alpha centers method. In this paper, the disadvantages of these two kinds of methods were analysed, and a new method of forward modeling and inversion (Tomography) based on boundary integral equations was proposed. This scheme successfuly overcomes the difficulties of the two formarly methods. It isn’t necessary to use the linearization approximation and calculate the Jacobi matrix. Numerical modeling results given in this paper showed that the computation speed of our method is fast, and there is not any special requirement for initial model, and satisfying results of tomography can be obtained in the case of great contrast of conductivity. So it has wide applications.

Keywords

2.5-D problem / boundary integral equations / resistivity tomography

References

[1] SasakiYutaka. 3-D resistivity inversion using the finite-element method. Geophysics, 1994, 59(11): 1839-1848

[2] ShimaHiromasa. Two-dimensional automatic resistivity inversion technique using alpha centers. Geophysics, 1990, 55(6): 682-694

[3] ShimaHiromasa. 2-D and 3-D resistivity image reconstruction using crosshole data. Geophysics, 1992, 57(10): 1270-1281

[4] MaoXianjin, BaoGuangshu. A new method for forward modeling in 2.5-D resistivity problem. J Cent South Univ Technol, 1997, 28(4): 307-310(in Chinese)

[5] SpiegelR J. Modeling resistivity anomalies from localized voids under irregular terrain. Geophysics, 1980, 45(7): 1164-1183