Piecewise 1D laterally constrained inversion of resistivity data

In a sedimentary environment, layered models are often capable of representing the actual geology more accurately than smooth minimum structure models. Furthermore, interval thicknesses and resistivities are often the parameters to which non-geophysicist experts can relate and base decisions on when using them in waste site remediation, groundwater modelling and physical planning. We present a laterally constrained inversion scheme for continuous resistivity data based on a layered earth model (1D). All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. The models are regularized through laterally equal constraints that tie interface depths and resistivities of adjacent layers. Prior information, e.g. originating from electric logs, migrates through the lateral constraints to the adjacent models, making resolution of equivalences possible to some extent. Information from areas with well-resolved parameters will migrate through the constraints in a similar way to help resolve the poorly constrained parameters. The estimated model is complemented by a full sensitivity analysis of the model parameters, supporting quantitative evaluation of the inversion result. Examples from synthetic 2D models show that the model recognition of a sublayered 2D wedge model is improved using the laterally constrained inversion approach when compared with a section of combined 1D models and when compared with a 2D minimum structure inversion. Case histories with data from two different continuous DC systems support the conclusions drawn from the synthetic example.

[1]  R. A. Overmeeren,et al.  Continuous vertical electrical sounding , 1988 .

[2]  E. Auken,et al.  MODEL RECOGNITION BY 1-D LATERALLY CONSTRAINED INVERSION OF RESISTIVITY DATA , 2002 .

[3]  R. Dmowska,et al.  International Geophysics Series , 1992 .

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

[5]  M. Descloîtres,et al.  Electromagnetic Investigations of Aquifers in the Grand Brulé Coastal Area of Piton de la Fournaise Volcano, Reunion Island , 1997 .

[6]  D. Oldenburg,et al.  Inversion of Induced-Polarization Data , 1993 .

[7]  J. T. Smith,et al.  Rapid inversion of two‐ and three‐dimensional magnetotelluric data , 1991 .

[8]  T. Dahlin 2D resistivity surveying for environmental and engineering applications , 1996 .

[9]  Fernando A. Monteiro Santos,et al.  1-D laterally constrained inversion of EM34 profiling data , 2004 .

[10]  K. Taylor,et al.  Use of transient electromagnetics to define local hydrogeology in an arid alluvial environment , 1992 .

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

[12]  Esben Auken,et al.  Layered and laterally constrained 2D inversion of resistivity data , 2004 .

[13]  K. Sørensen,et al.  Pulled Array Continuous Electrical Profiling , 1996 .

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

[15]  T. Dahlin,et al.  A numerical comparison of 2D resistivity imaging with 10 electrode arrays , 2004 .

[16]  H. Johansen,et al.  FAST HANKEL TRANSFORMS , 1979 .

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

[18]  H. Johansen A MAN/COMPUTER INTERPRETATION SYSTEM FOR RESISTIVITY SOUNDINGS OVER A HORIZONTALLY STRAFIFIED EARTH* , 1977 .

[19]  T. Ormos,et al.  A new procedure for the interpretation of VES data: 1.5-D simultaneous inversion method , 1999 .

[20]  M. Ali Ak,et al.  ‘European’ Association of Geoscientists & Engineers , 1997 .

[21]  T. Dahlin The development of DC resistivity imaging techniques , 2001 .

[22]  Examples of Transient Sounding for Ground‐Water Exploration in Sedimentary Aquifers , 1987 .

[23]  D. Oldenburg,et al.  Inversion of geophysical data using an approximate inverse mapping , 1991 .

[24]  Alain Tabbagh,et al.  Potential focusing: a new multielectrode array concept, simulation study, and field tests in archaeological prospecting , 1997 .

[25]  R. Barker,et al.  Rapid least-squares inversion of apparent resistivity pseudosections , 1994 .

[26]  P. McGillivray,et al.  Forward modeling and inversion of DC resistivity and MMR data , 1992 .

[27]  Michel Chouteau,et al.  Utility of cokriging GPR and hydrogeological data for aquifer evaluation , 2000 .

[28]  T. Dahlin,et al.  A comparison of smooth and blocky inversion methods in 2-D electrical imaging surveys , 2001 .

[29]  A. Tarantola,et al.  Generalized Nonlinear Inverse Problems Solved Using the Least Squares Criterion (Paper 1R1855) , 1982 .

[30]  D. W. Hall,et al.  Geophysical Investigation of an Unconsolidated Coastal Plain Aquifer System and the Underlying Bedrock Geology in Central New Jersey , 1990 .

[31]  T. Dahlin,et al.  Assessment of Two Automated Electrical Resistivity Data Acquisition Systems for Landfill Location Surveys: Two Case Studies , 1999 .

[32]  K. Sørensen Pulled Array Continuous Electrical Profiling , 1996 .

[33]  W. Menke Geophysical data analysis : discrete inverse theory , 1984 .