The inversion of data from very large three‐dimensional electrical resistivity tomography mobile surveys

ABSTRACT New developments in mobile resistivity meter instrumentation have made it possible to survey large areas with dense data coverage. The mobile system usually has a limited number of electrodes attached to a cable that is pulled along behind an operator so that a large area can be covered within a short time. Such surveys can produce three‐dimensional datasets with hundreds of thousands of electrodes positions and data points. Similarly, the inverse model used to interpret the data can have several hundred thousand cells. It is impractical to model such large datasets within a reasonable time on microcomputers used by many small companies employing standard inversion techniques. We describe a model segmentation technique that subdivides the finite‐element mesh used to calculate the apparent resistivity and Jacobian matrix values into a number of smaller meshes. A fast technique that optimizes the calculation of the Jacobian matrix values for multi‐channel systems was also developed. A one‐dimensional wavelet transform method was then used to compress the storage of the Jacobian matrix, in turn reducing the computer time and memory required to solve the least‐squares optimization equation to determine the inverse model resistivity values. The new techniques reduce the calculation time and memory required by more than 80% while producing models that differ by less than 1% from that obtained using the standard inversion technique with a single mesh. We present results using a synthetic model and a field dataset that illustrates the effectiveness of the proposed techniques.

[1]  L. S. Edwards,et al.  A modified pseudosection for resistivity and IP , 1977 .

[2]  S. Constable,et al.  Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data , 1990 .

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

[4]  Douglas W. Oldenburg,et al.  3-D inversion of airborne electromagnetic data parallelized and accelerated by local mesh and adaptive soundings , 2014 .

[5]  Oliver Kuras,et al.  Recent developments in the direct-current geoelectrical imaging method , 2013 .

[6]  Efficient 2D inversion of long ERT sections , 2014 .

[7]  D. Oldenburg,et al.  METHODS FOR CALCULATING FRÉCHET DERIVATIVES AND SENSITIVITIES FOR THE NON‐LINEAR INVERSE PROBLEM: A COMPARATIVE STUDY1 , 1990 .

[8]  Hansruedi Maurer,et al.  Geoelectric experimental design — Efficient acquisition and exploitation of complete pole-bipole data sets , 2011 .

[9]  M. Loke,et al.  Inversion of Data from Electrical Resistivity Imaging Surveys in Water-Covered Areas , 2004 .

[10]  N. Papadopoulos,et al.  Archaeological investigations in the shallow seawater environment with electrical resistivity tomography , 2015 .

[11]  R. Barker,et al.  Practical techniques for 3D resistivity surveys and data inversion1 , 1996 .

[12]  M. Loke,et al.  Rapid inversion of data from 2D resistivity surveys with electrode displacements , 2018 .

[13]  Apostolos Sarris,et al.  An algorithm for fast 3D inversion of surface electrical resistivity tomography data: application on imaging buried antiquities , 2011 .

[14]  Lucy S. Oxby,et al.  Computation of optimized arrays for 3-D electrical imaging surveys , 2014 .

[15]  D. Rucker,et al.  Smoothness‐constrained time‐lapse inversion of data from 3D resistivity surveys , 2014 .

[16]  D. Oldenburg,et al.  Fast inversion of large-scale magnetic data using wavelet transforms and a logarithmic barrier method , 2003 .

[17]  Paul B. Wilkinson,et al.  Optimized arrays for 2-D resistivity survey lines with a large number of electrodes , 2015 .

[18]  Tianbai Xu,et al.  Binning Method for Mapping Irregularly Distributed Continuous Resistivity Profiling Data onto a Regular Grid for 3-D Inversion , 2015 .

[19]  D. Oldenburg,et al.  NON-LINEAR INVERSION USING GENERAL MEASURES OF DATA MISFIT AND MODEL STRUCTURE , 1998 .

[20]  T. Dahlin,et al.  Properties and Effects of Measurement Errors on 2D Resistivity Imaging Surveying , 2003 .

[21]  David Beamish,et al.  Fundamentals of the capacitive resistivity technique , 2006 .

[22]  D. Oldenburg,et al.  A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems , 2004 .

[23]  Dale F. Rucker,et al.  Using marine resistivity to map geotechnical properties: a case study in support of dredging the Panama Canal , 2013 .

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

[25]  I. Daubechies Orthonormal bases of compactly supported wavelets , 1988 .

[26]  M. Loke Topographic Modelling in Electrical Imaging Inversion , 2000 .

[27]  Douglas W. Oldenburg,et al.  3-D inversion of induced polarization data3-D Inversion of IP Data , 2000 .

[28]  Kristofer Davis,et al.  Efficient 3D inversion of magnetic data via octree-mesh discretization, space-filling curves, and wavelets , 2013 .

[29]  T. Dahlin,et al.  Underwater ERT surveying in water with resistivity layering with example of application to site investigation for a rock tunnel in central Stockholm , 2018 .

[30]  A. Dey,et al.  Resistivity modeling for arbitrarily shaped three-dimensional structures , 1979 .

[31]  Dale F. Rucker,et al.  Three-dimensional electrical resistivity model of a nuclear waste disposal site , 2009 .