Two‐dimensional joint inversion of radiomagnetotelluric and direct current resistivity data

An algorithm for the two-dimensional (2D) joint inversion of radiomagnetotelluric and direct current resistivity data was developed. This algorithm can be used for the 2D inversion of apparent resistivity data sets collected by multi-electrode direct current resistivity systems for various classical electrode arrays (Wenner, Schlumberger, dipole-diplole, pole-dipole) and radiomagnetotelluric measurements jointly. We use a finite difference technique to solve the Helmoltz and Poisson equations for radiomagnetotelluric and direct current resistivity methods respectively. A regularized inversion with a smoothness constrained stabilizer was employed to invert both data sets. The radiomagnetotelluric method is not particularly sensitive when attempting to resolve near-surface resistivity blocks because it uses a limited range of frequencies. On the other hand, the direct current resistivity method can resolve these near-surface blocks with relatively greater accuracy. Initially, individual and joint inversions of synthetic radiomagnetotelluric and direct current resistivity data were compared and we demonstrated that the joint inversion result based on this synthetic data simulates the real model more accurately than the inversion results of each individual method. The developed 2D joint inversion algorithm was also applied on a field data set observed across an active fault located close to the city of Kerpen in Germany. The location and depth of this fault were successfully determined by the 2D joint inversion of the radiomagnetotelluric and direct current resistivity data. This inversion result from the field data further validated the synthetic data inversion results.

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

[2]  B. Tezkan,et al.  Radiomagnetotelluric and Direct Current Resistivity measurements for the characterization of conducting soils , 2007 .

[3]  Per Christian Hansen,et al.  REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems , 1994, Numerical Algorithms.

[4]  M. Everett,et al.  Two-Dimensional Nonlinear Magnetotelluric Inversion Using a Genetic Algorithm , 1993 .

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

[6]  M. Bastani,et al.  Groundwater exploration using combined controlled-source and radiomagnetotelluric techniques , 2005 .

[7]  S. Friedel,et al.  Resolution, stability and efficiency of resistivity tomography estimated from a generalized inverse approach , 2003 .

[8]  Bülent Tezkan,et al.  A Review Of Environmental Applications Of Quasi-Stationary Electromagnetic Techniques , 1999 .

[9]  E. Auken,et al.  A survey of current trends in near-surface electrical and electromagnetic methods , 2006 .

[10]  K. Vozoff,et al.  The joint use of coincident loop transient electromagnetic and Schlumberger sounding to resolve layered structures , 1985 .

[11]  J. D. McNeill,et al.  7. Geological Mapping Using VLF Radio Fields , 1991 .

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

[13]  Christoph Schwarzbach,et al.  Two‐dimensional inversion of direct current resistivity data using a parallel, multi‐objective genetic algorithm , 2005 .

[14]  L. Ahorner.  Untersuchungen zur quartären Bruchtektonik der Niederrheinischen Bucht , 1962 .

[15]  M. N. Berdichevsky,et al.  Basic principles of interpretation of magnetotelluric sounding curves , 1976 .

[16]  Kalyanmoy Deb,et al.  A fast and elitist multiobjective genetic algorithm: NSGA-II , 2002, IEEE Trans. Evol. Comput..

[17]  M. Meju,et al.  Joint inversion of TEM and distorted MT soundings, Some effective practical considerations , 1996 .

[18]  M. Meju,et al.  Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints , 2004 .

[19]  T. Harinarayana.  Combination of EM and DC Measurements For Upper Crustal Studies , 1999 .

[20]  M. E. Candansayar,et al.  Deep resistivity structure of the Dikili‐Bergama region, west Anatolia, revealed by two‐ dimensional inversion of vertical electrical sounding data , 2006 .

[21]  P. Turberg,et al.  Joint application of radio-magnetotelluric and electrical imaging surveys in complex subsurface environments , 1996 .

[22]  Niklas Linde,et al.  Characterization of a fractured granite using radio magnetotelluric (RMT) data , 2004 .

[23]  A. T. Basokur,et al.  Detecting small‐scale targets by the 2D inversion of two‐sided three‐electrode data: application to an archaeological survey , 2001 .

[24]  D. Oldenburg,et al.  Estimating depth of investigation in DC resistivity and IP surveys , 1999 .

[25]  Mary M. Poulton,et al.  Computational Neural Networks for Geophysical Data Processing , 2001 .

[26]  David L.B. Jupp,et al.  Joint Inversion of Geophysical Data , 2007 .

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

[28]  Bülent Tezkan,et al.  A comparison of different radiomagnetotelluric data inversion methods for buried waste sites , 2006 .

[29]  Gregory A. Newman,et al.  3D inversion of a scalar radio magnetotelluric field data set , 2003 .

[30]  D. Eaton,et al.  Magnetotelluric and teleseismic study across the Snowbird Tectonic Zone, Canadian Shield: A Neoarchean mantle suture? , 2002 .

[31]  Douglas W. Oldenburg,et al.  Calculation of sensitivities for the frequency-domain electromagnetic problem , 1994 .

[32]  Bülent Tezkan,et al.  Two-dimensional radiomagnetotelluric investigation of industrial and domestic waste sites in Germany , 2000 .

[33]  A. Hördt,et al.  A joint application of radiomagnetotellurics and transient electromagnetics to the investigation of a waste deposit in Cologne (Germany) , 1996 .

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

[35]  Gregory A. Newman,et al.  Image appraisal for 2-D and 3-D electromagnetic inversion , 2000 .

[36]  A. Tikhonov,et al.  Numerical Methods for the Solution of Ill-Posed Problems , 1995 .

[37]  M. Emin Candansayar,et al.  Two‐dimensional inversion of magnetotelluric data with consecutive use of conjugate gradient and least‐squares solution with singular value decomposition algorithms , 2007 .

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

[39]  Jonathan Chambers,et al.  Improved strategies for the automatic selection of optimized sets of electrical resistivity tomography measurement configurations , 2006 .

[40]  T. Czerwiński,et al.  Regionalne badania magnetotelluryczne w polskich Karpatach wschodnich , 2003 .