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

In this paper, a new simple, efficient and semiautomatic algorithm including depth weighting constraint is introduced for 2d 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 2d 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 perezflores et al. (2001), but here in addition to improving results for dipole-dipole array, its productivity is demonstrated for other geoelectrical 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.

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

In this paper, a new simple, efficient and semiautomatic algorithm including depth weighting constraint is introduced for 2D 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 2D 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 PerezFlores et al. (2001), but here in addition to improving results for dipole-dipole array, its productivity is demonstrated for other geoelectrical 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.