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

In this paper, a new simple, efficient and semi-automatic algorithm including depth weighting constraint is introduced for ۲-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 ۲-D forward operator formula obtained based on integral equations (IE) by Perez- Flores et al. (۲۰۰۱), for the inversion algorithm. Inversion algorithm is iterative and regularization parameter and depth weighting exponent are the critical parameters that have default values of ۰.۱ and ۱, respectively. The presented technique was used only for dipole-dipole array by Perez-Flores et al. (۲۰۰۱), 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 Res۲dmod software are utilized to investigate the performance of the algorithm through comparing the results with Res۲dinv 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.

Keywords:

DC resistivity , depth weighting , integral equation , inversion , Res۲dinv