Two-dimensional inversion of direct current resistivity data incorporating topography by using finite difference techniques with triangle cells: Investigation of Kera fault zone in western Crete

In this study, we suggest the use of a finite difference (FD) forward solution with triangular grid to incorporate topography into the inverse solution of direct current resistivity data. A new inversion algorithm was developed that takes topography into account with finite difference and finite element forward solution by using triangular grids. Using the developed algorithm, surface topography could also be incorporated by using triangular cells in a finite difference forward solution. Initially, the inversion algorithm was tested for two synthetic data sets. Inversion of synthetic data with the finite difference forward solution gives accurate results as well as inversion with finite element forward solution and requires less CPU time. The algorithm was also tested with a field data set acquired across the Kera fault located in western Crete, Greece. The fault location and basement depth of sedimentary units were resolved by the developed algorithm. These inversion results showed that if underground st...

[1]  Myeong-Jong Yi,et al.  Geophysical investigation of tumuli by means of surface 3D Electrical Resistivity Tomography , 2010 .

[2]  Bülent Tezkan,et al.  Two‐dimensional joint inversion of radiomagnetotelluric and direct current resistivity data , 2008 .

[3]  E. Haber,et al.  A GCV based method for nonlinear ill-posed problems , 2000 .

[4]  D. LaBrecque,et al.  Occam's inversion of 3-D electrical resistivity tomography. Three-dimensional Electromagnetics , 1999 .

[5]  Yutaka Sasaki,et al.  3-D resistivity inversion using the finite-element method , 1994 .

[6]  J. Coggon Electromagnetic and electrical modeling by the finite element method , 1971 .

[7]  J. Torquil Smith,et al.  The forward problem of electromagnetic induction: accurate finite-difference approximations for two-dimensional discrete boundaries with arbitrary geometry , 1997 .

[8]  I. Fountoulis,et al.  Paleoseismological investigations along the Kera fault zone, Western Crete: implications for seismic hazard assessment , 2001 .

[9]  Klaus Spitzer,et al.  The three‐dimensional DC sensitivity for surface and subsurface sources , 1998 .

[10]  Seong-Jun Cho,et al.  Three‐dimensional imaging of subsurface structures using resistivity data , 2001 .

[11]  Nariida C. Smith,et al.  Two-Dimensional DC Resistivity Inversion for Dipole-Dipole Data , 1984, IEEE Transactions on Geoscience and Remote Sensing.

[12]  J. E. Szymanski,et al.  The effect of terrain topography on commonly used resistivity arrays , 1999 .

[13]  Gerald W. Hohmann,et al.  Topographic effects in resistivity and induced-polarization surveys , 1980 .

[14]  Alan C. Tripp,et al.  Two-dimensional resistivity inversion , 1984 .

[15]  R. Barker,et al.  Rapid least-squared inversion of apparent resisitivity pseudosections by a quasi-Newton method , 1996 .

[16]  A. Dey,et al.  Resistivity modelling for arbitrarily shaped two-dimensional structures , 1979 .

[17]  M. E. Candansayar Two-dimensional individual and joint inversion of three- and four-electrode array dc resistivity data , 2008 .

[18]  Yutaka Sasaki,et al.  Two‐dimensional joint inversion of magnetotelluric and dipole‐dipole resistivity data , 1989 .

[19]  J. T. Weaver Mathematical methods for geo-electromagnetic induction , 1994 .

[20]  Lun-Tao Tong,et al.  Incorporation of topography into two-dimensional resistivity inversion , 1990 .

[21]  Alan G. Green,et al.  Experimental design: Electrical resistivity data sets that provide optimum subsurface information , 2004 .

[22]  P. Soupios,et al.  Geological and Electrical Resistivity Tomography surveys applied to model the tectonic environment of the Kastelli-Kissamou basin, northwestern Crete, Greece , 2010 .

[23]  M. E. Candansayar,et al.  Incorporating topography into 2D resistivity modeling using finite-element and finite-difference approaches , 2008 .

[24]  H. Maurer,et al.  Advances in three-dimensional geoelectric forward solver techniques , 2009 .

[25]  H. T. Holcombe,et al.  Three-dimensional terrain corrections in resistivity surveys , 1984 .