Data-driven selection of the minimum-gradient support parameter in time-lapse focused electric imaging

We have considered the problem of the choice of the minimum-gradient support (MGS) parameter in focused inversion for time-lapse (TL) electric resistivity tomography. Most existing approaches have relied either on an arbitrary choice of this parameter or one based on the prior information, such as the expected contrast in the TL image. We have decided to select the MGS parameter using a line search based on the value of the TL data root-mean-square misfit at the first iteration of the nonlinear inversion procedure. The latter was based on a Gauss-Newton scheme minimizing a regularized objective function in which the regularization functional was defined by the MGS functional. The regularization parameter was optimized to achieve a certain target level, following the Occam principles. We have validated our approach on a synthetic benchmark using a complex and heterogeneous model and determined its effectiveness on electric tomography TL data collected during a salt tracer experiment in fractured limestone. Our results have determined that the approach was successful in retrieving the focused anomaly and did not rely on prior information.

[1]  R. Lytle,et al.  Investigation of fracture flow paths using alterant geophysical tomography , 1986 .

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

[3]  W. Daily,et al.  The effects of noise on Occam's inversion of resistivity tomography data , 1996 .

[4]  Michael S. Zhdanov,et al.  Focusing geophysical inversion images , 1999 .

[5]  Rob S. MacLeod,et al.  Inverse electrocardiography by simultaneous imposition of multiple constraints , 1999, IEEE Transactions on Biomedical Engineering.

[6]  Douglas LaBrecque,et al.  Difference Inversion of ERT Data: a Fast Inversion Method for 3-D In Situ Monitoring , 2001 .

[7]  E. Miller,et al.  Efficient determination of multiple regularization parameters in a generalized L-curve framework , 2002 .

[8]  H. Vereecken,et al.  Imaging and characterisation of subsurface solute transport using electrical resistivity tomography (ERT) and equivalent transport models , 2002 .

[9]  Frederick D. Day-Lewis,et al.  Time‐lapse inversion of crosswell radar data , 2002 .

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

[11]  Michael S. Zhdanov,et al.  Minimum support nonlinear parametrization in the solution of a 3D magnetotelluric inverse problem , 2004 .

[12]  John P. Castagna,et al.  Inverse Spectral Decomposition , 2004 .

[13]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[14]  Jonathan B. Ajo-Franklin,et al.  Applying Compactness Constraints to Differential Traveltime Tomography , 2007 .

[15]  A. Binley,et al.  Quantitative imaging of solute transport in an unsaturated and undisturbed soil monolith with 3‐D ERT and TDR , 2008 .

[16]  James P. McNamara,et al.  Application of time-lapse ERT imaging to watershed characterization , 2008 .

[17]  A. Hördt,et al.  A new sensitivity-controlled focusing regularization scheme for the inversion of induced polarization data based on the minimum gradient support , 2008 .

[18]  J. G. Kim,et al.  4-D inversion of DC resistivity monitoring data acquired over a dynamically changing earth model , 2009 .

[19]  J. Gisbert,et al.  Automated monitoring of coastal aquifers with electrical resistivity tomography , 2009 .

[20]  H. Vereecken,et al.  Investigating Preferential Flow Processes in a Forest Soil Using Time Domain Reflectometry and Electrical Resistivity Tomography , 2010 .

[21]  A. Flores-Orozco,et al.  Temperature‐calibrated imaging of seasonal changes in permafrost rock walls by quantitative electrical resistivity tomography (Zugspitze, German/Austrian Alps) , 2010 .

[22]  A. Ward,et al.  Use of electrical imaging and distributed temperature sensing methods to characterize surface water–groundwater exchange regulating uranium transport at the Hanford 300 Area, Washington , 2010 .

[23]  I. Cho,et al.  Time-lapse inversion of 2D resistivity monitoring data with a spatially varying cross-model constraint , 2011 .

[24]  Frédéric Nguyen,et al.  Assessing the contribution of electrical resistivity tomography (ERT) and self-potential (SP) methods for a water well drilling program in fractured/karstified limestones , 2011 .

[25]  Panagiotis Tsourlos,et al.  4D active time constrained resistivity inversion , 2011 .

[26]  M. Rosas Carbajal,et al.  Focused time-lapse inversion of radio and audio magnetotelluric data , 2012, 1701.02529.

[27]  Niklas Linde,et al.  Imaging and quantifying salt-tracer transport in a riparian groundwater system by means of 3D ERT monitoring , 2012 .

[28]  Hisashi Utada,et al.  A regularized three-dimensional magnetotelluric inversion with a minimum gradient support constraint , 2012 .

[29]  Frédéric Nguyen,et al.  A salt tracer test monitored with surface ERT to detect preferential flow and transport paths in fractured/karstified limestones , 2012 .

[30]  Lee Slater,et al.  Challenges and opportunities for fractured rock imaging using 3D cross-borehole electrical resistivity , 2015 .

[31]  K. Singha,et al.  Advances in interpretation of subsurface processes with time‐lapse electrical imaging , 2015 .

[32]  Frédéric Nguyen,et al.  Quantitative temperature monitoring of a heat tracing experiment using cross-borehole ERT , 2015 .

[33]  Esben Auken,et al.  Sharp spatially constrained inversion with applications to transient electromagnetic data , 2015 .

[34]  Esben Auken,et al.  Generalized focusing of time-lapse changes with applications to direct current and time-domain induced polarization inversions , 2015 .