Modified total variation methods for three‐dimensional electrical resistance tomography inverse problems

Images from geophysical surveys and data inversion can be valuable information for hydrogeologic site characterization and monitoring. This paper focuses on solving nonlinear electrical resistance tomography (ERT) inverse problems; the methods are applicable to other inverse problems, such as groundwater flow. The total variation (TV) method is used to constrain ERT inverse solutions. The classic TV functions are modified to obtain differentiability with respect to parameters, which allows the use of gradient search optimization methods. The modified total variation (mTV) method constrains ERT inverse solutions effectively. Sometimes the mTV constraints cause important parameter features to be smoothed out. To eliminate excessive smoothing, two successive partial variation relaxation (SPVR) methods are devised for use with the mTV method. Computational experiments demonstrate that mTV works well in identifying very sharp discontinuities and smoothly distributed parameters and that SPVR improves mTV solutions for both problem types. The mTV and SPVR methods are applied to synthetic layered systems with and without measurement error, and a comparison of results for a large three-dimensional ERT problem with block anomaly presented by Zhang et al. [1995] demonstrates the superiority of the mTV and SPVR methods. A final example examines a smoothly varying case and explores the impact of two computational parameters introduced in the development of mTV and SPVR. The example shows also that SPVR can be used directly to invert data and not necessarily as an adjunct to mTV.

[1]  M. A. Medina,et al.  Incorporating Auxiliary Geophysical Data into Ground‐Water Flow Parameter Estimation , 1997 .

[2]  Jie Zhang,et al.  3-D resistivity forward modeling and inversion using conjugate gradients , 1995 .

[3]  P. Mäkilä,et al.  Techniques in the Maximum Likelihood Estimation of the Covariance Matrix , 1996 .

[4]  Irshad R. Mufti,et al.  FINITE‐DIFFERENCE RESISTIVITY MODELING FOR ARBITRARILY SHAPED TWO‐DIMENSIONAL STRUCTURES , 1976 .

[5]  Mark K. Segar.  A SLAP for the masses , 1989 .

[6]  W. Yeh.  Review of Parameter Identification Procedures in Groundwater Hydrology: The Inverse Problem , 1986 .

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

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

[9]  Peter K. Kitanidis,et al.  The minimum structure solution to the inverse problem , 1997 .

[10]  Kohn,et al.  Variational constraints for electrical-impedance tomography. , 1990, Physical review letters.

[11]  D. Goldfarb.  A family of variable-metric methods derived by variational means , 1970 .

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

[13]  Rexford M. Morey,et al.  TOMOGRAPHIC SITE CHARACTERIZATION USING CPT, ERT, AND GPR , 1999 .

[14]  P. M. Berg,et al.  A total variation enhanced modified gradient algorithm for profile reconstruction , 1995 .

[15]  M. Eppstein,et al.  Efficient three‐dimensional data inversion: Soil characterization and moisture Monitoring from cross‐well ground‐penetrating radar at a Vermont Test Site , 1998 .

[16]  D. LaBrecque,et al.  Monitoring Air Sparging Using Resistivity Tomography , 1996 .

[17]  M. Eppstein,et al.  Optimal 3-D traveltime tomography , 1998 .

[18]  C. G. Broyden.  The Convergence of a Class of Double-rank Minimization Algorithms 1. General Considerations , 1970 .

[19]  R. P. Alger.  Interpretation Of Electric Logs In Fresh Water Wells In Unconsolidated Formations , 1966 .

[20]  D. Urish,et al.  Electrical resistivity—hydraulic conductivity relationships in glacial outwash aquifers , 1981 .

[21]  S. P. Neuman,et al.  Estimation of aquifer parameters under transient and steady-state conditions: 2 , 1986 .

[22]  D. Dobson,et al.  An image-enhancement technique for electrical impedance tomography , 1994 .

[23]  D. McLaughlin,et al.  A Reassessment of the Groundwater Inverse Problem , 1996 .

[24]  J. Nitao,et al.  Electrical resistivity tomography of vadose water movement , 1992 .

[25]  D. Isaacson.  Distinguishability of Conductivities by Electric Current Computed Tomography , 1986, IEEE Transactions on Medical Imaging.

[26]  Youcef Saad,et al.  A Basic Tool Kit for Sparse Matrix Computations , 1990 .

[27]  Stephen K. Park,et al.  Monitoring leaks from storage ponds using resistivity methods , 1991 .

[28]  R. W. Andrews,et al.  Sensitivity Analysis for Steady State Groundwater Flow Using Adjoint Operators , 1985 .

[29]  Douglas LaBrecque,et al.  Monitoring an underground steam injection process using electrical resistance tomography , 1993 .

[30]  Mark Pilkington,et al.  Natural Smoothness Constraints in Cross-Hole Seismic TOMOGRAPHY1 , 1992 .

[31]  D. Shanno.  Conditioning of Quasi-Newton Methods for Function Minimization , 1970 .

[32]  G. E. Archie.  The electrical resistivity log as an aid in determining some reservoir characteristics , 1942 .

[33]  W. Kelly.  Geoelectric Sounding for Estimating Aquifer Hydraulic Conductivity , 1977 .

[34]  Abelardo Ramirez,et al.  Electrical resistance tomography during in-situ trichloroethylene remediation at the Savannah River Site , 1995 .

[35]  M. Eppstein,et al.  SIMULTANEOUS ESTIMATION OF TRANSMISSIVITY VALUES AND ZONATION , 1996 .

[36]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

[37]  C. Vogel,et al.  Analysis of bounded variation penalty methods for ill-posed problems , 1994 .

[38]  Mary C. Hill,et al.  Methods and Guidelines for Effective Model Calibration , 2000 .

[39]  R. Fletcher,et al.  A New Approach to Variable Metric Algorithms , 1970, Comput. J..

[40]  S. Gutman.  Identification of discontinuous parameters in flow equations , 1990 .

[41]  C. G. Broyden.  The Convergence of a Class of Double-rank Minimization Algorithms 2. The New Algorithm , 1970 .

[42]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .