INVERSION OF ELECTRICAL RESISTIVITY TOMOGRAPHY DATA DERIVING FROM 3D STRUCTURES


P. Tsourlos
Abstract

In this work the effectiveness of 2D and 3D algorithms for inverting Electrical Resistivity Tomography (ERT) data deriving from 3D structures is studied. Further, an analysis of data-collection strategies in the case of 3D structures is being carried out. Dense 2D measurements are nconsidered a practical tool for mapping 3D structures given the current limitations in ERT hardware. To perform the tests 2D and a 3D inversion programs are used. Both schemes use a forward model based on a 2.5D and 3D finite element scheme respectively. For both the 2D and 3D cases a fully non-linear inversion scheme based on a smoothness constrained algorithm is used. The Jacobian matrix is calculated using the adjoined equation technique. Comparisons are being carried out by means of synthetic examples for 3D models and dense 2D measurements with their axis parallel to the X (X-lines) and/or Y (Y lines) directions. For the case of 3D structures and 2D inversion tests results illustrate that both X-line, Y-line measurements are required to delineate the modeling body. However, when 3D inversion is considered either Xline or Y-line measurements are adequate to produce good quality reconstructions of the subsurface. Overall, results clearly illustrate the superiority of 3D over 2D inversion schemes in the case of 3D structures both in view of quality and logistics. Despite the increased computational time required by 3D inversion schemes, good quality results can be produced. Further, 2D inversion techniques require effectively a double amount of measurements to produce acceptable results. The ongoing advancement of fast computers renders the described approach of combining dense 2-D measurement with 3D inversion practical for routine data treatment.

Article Details
  • Section
  • Geophysics
Downloads
Download data is not yet available.
References
Chambers J.E., 2001. The application of 3D electrical tomography to the investigations of brownfield sites. PhD Thesis, University of Sheffield.
Constable, S. Parker, R., & Constable C, 1987. Occam's inversion: A practical algorithm for generating smooth models from electromagnetic sounding data. Geophysics, 52:289-300.
deGroot-Hedlin, C, and Constable, S., 1990.Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data. Geophysics, 55:1613-1624.
Ellis, R., & Oldenburg D.W., 1994. Applied geophysical inversion.Geophys.J.Int., 116, 5-11.
Loke M.H. and Barker R.D.,1996a. Rapid least-squares inversion of apparent resistivity pseudosections using a quasi-Newton method. Geophysical Prospecting, 44, 131-152.
Loke M.H. and Barker R.D.,1996b. Practical techniques for 3D resistivity surveys and data inversion. Geophysical Prospecting, 44, 499-523.
McGillivray, P., & Oldenburg, D., 1990, Methods for calculating Frechet derivatives and sensitivities for the nonlinear inverse problem: A comparative study: Geophysical Prospecting, 38, 499-524.
Park, S.K. & Van, G.P. 1991. Inversion of pole-pole data for 3-D resistivity structures beneath arrays of electrodes. Geophysics, 56, 951-960.
Sasaki, Y., 1989. 2-D joint inversion of magnetotelluric and dipole-dipole resistivity data. Geophysics, 54:254-262.
Sasaki, Y., 1992. Resolution of resistivity tomography inferred from numerical simulation. Geophysical Prospecting, 40:453-464.
Tsourlos P., 1995. Modelling interpretation and inversion of multielectrode resistivity survey data. Ph.D. Thesis, University of York.
Tsourlos, P., Szymanski, J. & Tsokas G., 1999, The effect of terrain topography on commonly used resistivity arrays: Geophysics, 64, 1357-1363.
Tsourlos P. & Ogilvy R. 1999. An algorithm for the 3-D Inversion of Tomographic Resistivity and Induced Polarisation data: Preliminary Results. Journal of the Balkan Geophysical Society, 2, 2, 30-45.
Xu, B., 1993. Development of electrical resistivity imaging methods for geological and archaeological prospecting. Ph.D. Thesis, University of Durham.
Most read articles by the same author(s)