VES/TEM 1D joint inversion by using Controlled Random Search (CRS) algorithm

Abstract Electrical (DC) and Transient Electromagnetic (TEM) soundings are used in a great number of environmental, hydrological, and mining exploration studies. Usually, data interpretation is accomplished by individual 1D models resulting often in ambiguous models. This fact can be explained by the way as the two different methodologies sample the medium beneath surface. Vertical Electrical Sounding (VES) is good in marking resistive structures, while Transient Electromagnetic sounding (TEM) is very sensitive to conductive structures. Another difference is VES is better to detect shallow structures, while TEM soundings can reach deeper layers. A Matlab program for 1D joint inversion of VES and TEM soundings was developed aiming at exploring the best of both methods. The program uses CRS – Controlled Random Search – algorithm for both single and 1D joint inversions. Usually inversion programs use Marquadt type algorithms but for electrical and electromagnetic methods, these algorithms may find a local minimum or not converge. Initially, the algorithm was tested with synthetic data, and then it was used to invert experimental data from two places in Parana sedimentary basin (Bebedouro and Pirassununga cities), both located in Sao Paulo State, Brazil. Geoelectric model obtained from VES and TEM data 1D joint inversion is similar to the real geological condition, and ambiguities were minimized. Results with synthetic and real data show that 1D VES/TEM joint inversion better recovers simulated models and shows a great potential in geological studies, especially in hydrogeological studies.

[1]  M. Sambridge.  Geophysical inversion with a neighbourhood algorithm—II. Appraising the ensemble , 1999 .

[2]  Josef Tvrdík,et al.  The controlled random search algorithm in optimizing regression models , 1995 .

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

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

[5]  J. Porsani,et al.  TDEM survey in urban environmental for hydrogeological study at USP campus in São Paulo city, Brazil , 2012 .

[6]  Michel Menvielle,et al.  Stochastic interpretation of magnetotelluric data, comparison of methods , 2007 .

[7]  O. Koefoed,et al.  a Fast Method for Determining the Layer Distribution from the Raised Kernel Function in Geoelegtrical SOUNDING , 1970 .

[8]  D. Ghosh.  THE APPLICATION OF LINEAR FILTER THEORY TO THE ' DIRECT INTERPRETATION OF GEOELECTRICAL RESISTIVITY SOUNDING MEASUREMENTS * , 1971 .

[9]  J. Porsani,et al.  TDEM survey in an area of seismicity induced by water wells in Paraná sedimentary basin, Northern São Paulo State, Brazil , 2012 .

[10]  Fethi Tarik Bendimerad,et al.  Controlled Random Search Optimization for Linear Antenna Arrays , 2006 .

[11]  Gerald W. Hohmann,et al.  Nonlinear magnetic inversion using a random search method , 1983 .

[12]  Wyn L. Price,et al.  A Controlled Random Search Procedure for Global Optimisation , 1977, Comput. J..

[13]  S. Rhee,et al.  Determination of optimal welding conditions with a controlled random search procedure , 2005 .

[14]  Albert Tarantola,et al.  Monte Carlo sampling of solutions to inverse problems , 1995 .

[15]  Gongwen Wang,et al.  3D geological modeling based on gravitational and magnetic data inversion in the Luanchuan ore region, Henan Province, China , 2012 .

[16]  D. Ghosh,et al.  Inverse filter coefficients for the computation of apparent resistivity standard curves for a horizontally stratified earth , 1971 .

[17]  H. El-Kaliouby,et al.  Comparative study of local versus global methods for 1D joint inversion of direct current resistivity and time-domain electromagnetic data , 2010 .

[18]  Niels B. Christensen,et al.  Optimized fast Hankel transform filters , 1990 .

[19]  François Baumgartner,et al.  CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys , 2006, Comput. Geosci..

[20]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[21]  M. Goldman,et al.  On reducing ambiguity in the interpretation of transient electromagnetic sounding data1 , 1994 .

[22]  M. Schmutz,et al.  Determination of the 3D structure of an earthflow by geophysical methods: The case of Super Sauze, in the French southern Alps , 2009 .

[23]  Mrinal K. Sen,et al.  Nonlinear one-dimensional seismic waveform inversion using simulated annealing , 1991 .

[24]  O. Koefoed.  A note on the linear filter method of interpreting resistivity sounding data , 1972 .

[25]  M. Descloîtres,et al.  Mapping Coastal Aquifers by Joint Inversion of DC and TEM Soundings‐Three Case Histories , 2001 .

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

[27]  Radek Krpec,et al.  Stochastic algorithms in nonlinear regression , 2000 .

[28]  M. Sambridge,et al.  Monte Carlo analysis of inverse problems , 2002 .

[29]  J. Porsani,et al.  Geoelectrical characterization using joint inversion of VES/TEM data: a case study in Paraná Sedimentary Basin, São Paulo State, Brazil , 2014 .

[30]  H. Johansen,et al.  AN INTERACTIVE COMPUTER/GRAPHIC-DISPLAY-TERMINAL SYSTEM FOR INTERPRETATION OF RESISTIVITY SOUNDINGS* , 1975 .

[31]  M. Sambridge.  Geophysical inversion with a neighbourhood algorithm—I. Searching a parameter space , 1999 .

[32]  T. Alkhalifah,et al.  Residual dip moveout in VTI media , 2005 .

[33]  John F. Ferguson,et al.  Constrained inversion of seismic refraction data using the controlled random search , 2000 .

[34]  Max A. Meju,et al.  Simple relative space–time scaling of electrical and electromagnetic depth sounding arrays: implications for electrical static shift removal and joint DC‐TEM data inversion with the most‐squares criterion , 2005 .

[35]  Olivier Maquaire,et al.  Joint Electrical and Time Domain Electromagnetism (TDEM) Data Inversion Applied to the Super Sauze Earthflow (France) , 2000 .