A Semi-Automatic 2-D Linear Inversion Algorithm Including Depth Weighting Function for DC Resistivity Data: A Case Study on Archeological Data Sets of Pompeii

Document Type : Research Article


1 Ph.D. Student, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran

2 Associate Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran


In this paper, a new simple, efficient and semi-automatic algorithm including depth weighting constraint is introduced for 2-D DC resistivity data inversion. Inversion procedure is linear; however, DC resistivity data inversion is generally nonlinear due to the nonlinearity of Maxwell’s equation relative to resistivity (conductivity). We took the advantage of the 2-D forward operator formula obtained based on integral equations (IE) by Perez- Flores et al. (2001), for the inversion algorithm. Inversion algorithm is iterative and regularization parameter and depth weighting exponent are the critical parameters that have default values of 0.1 and 1, respectively. The presented technique was used only for dipole-dipole array by Perez-Flores et al. (2001), but here in addition to improving results for dipole-dipole array, its productivity is demonstrated for other geo-electrical arrays such as Wenner alfa, Wenner Schlumberger. Three synthetic data sets computed by Res2dmod software are utilized to investigate the performance of the algorithm through comparing the results with Res2dinv software output sections. Finally, the algorithm is applied on an archeological data set of Pompeii, which was collected by dipole-dipole array. IE inversion algorithm lead to satisfactory inversion models for both synthetic and real cases which reconstruct the subsurface better than or as well as that of the software.


Main Subjects

Auken, E. and Christiansen, A. V., 2004, Layered and laterally constrained 2D inversion of resistivity data, Geophysics, 69(3), 752–761.
Boonchaisuk, S., Vachiratienchai, C. and Siripunvaraporn, W., 2008, Two-dimensional direct current (DC) resistivity inversion: data space Occam's approach, Phys. Earth Planet. Inter, 168, 204–211.
Cockett, R., 1971, Direct Current Resistivity Inversion using Various Objective Functions.
Coggon, J.H., Electromagnetic and electrical modelling by the finite element method, Geophysics,36, 132–155.
Dieter, K., Paterson, N.R. and Grant, F.S., 1969, IP and resistivity type curves for three-dimensional bodies, Geophysics,34, 615–632.
Günther, T., Rücker, C. and Spitzer, K., 2006, Three-dimensional modelling and inversion of DC resistivity data incorporating topography—II. Inversion, Geophysical Journal International, 166, 506-17.
Jackson, P.D., Earl, S.J. and Reece, G.J., 2001, 3D resistivity inversion using 2D measurements of the electric field, Geophysical Prospecting, 49, 26-39.
Li, C.W., Bin X., Qiang, J.K. and Lü, Y.Z., 2012, Multiple linear system techniques for 3D finite element method modeling of direct current resistivity, Journal of Central South University, 19, 424-32.
Li, Y. and Oldenburg, D. W, 1996, 3-D inversion of magnetic data, Geophysics, 61, 394-408.
Li, Y. and Oldenburg, D. W., 1998, 3D inversion of gravity data, Geophysics, 63, 109-19.
Loke, M.H., Acworth, R. and Dahlin, T., 2003, A comparison of smooth and blocky inversion methods in 2D electrical imaging surveys, Exploration Geophysics, 34(1), 182-87.
Loke, M.H. and Dahlin, T., 1997, A combined Gauss-Newton and Quasi-Newton inversion method for the interpretation of apparent resistivity pseudo-sections, In 3rd EEGS Meeting,
Mufti, I.R., 1976, Finite-difference resistivity modeling for arbitrarily shaped two-dimensional structures, Geophysics,41, 62–78.
Narayan, S., Dusseault, M.B. and Nobes, D.C., 1994, Inversion techniques applied to resistivity inverse problems, Inverse Problems, 10(3), 669-686."
Perez-Flores, M.A., Méndez-Delgado S. and Gomez-Treviño, E., 2001, Imaging low frequency and dc electromagnetic fields using a simple linear approximation. Geophysics, 66, 1067–1081.
Sasaki, Y., 1994, 3D resistivity inversion using the finite-element method,Geophysics, 59, 1839–1848.
Smith, N.C. and Vozoff, K., 1984, Two dimensional DC resistivity inversion for dipole-dipole data, IEEE Trans. Gemeience Remote Sensing, 22, 21-28.
Timothy, W., Greenhalgh, S., Zhou, B., Greenhalgh, M. and Marescot, L., 2015, Resistivity inversion in 2-D anisotropic media: numerical experiments, GeophysicalJournalInternational, 201, 247-66.
Varfinezhad, R. and Oskooi, B., 2019, 2D DC resistivity forward modeling based on the integral equation method and a comparison with the RES2DMOD results, Journal of Earth and Space Physics. doi:10.22059/jesphys. 2019.260824. 1007020.
Wei, W., Wu, X., and Spitzer, K., 2013, Three-dimensional DC anisotropic resistivity modelling using finite elements on unstructured grids, GeophysicalJournalInternational, 193, 734-746.
Zhdanov, M.S., 2009, Geophysical electromagnetic theory and methods, Elsevier.