A geostatistically based inverse model for electrical resistivity surveys and its applications to vadose zone hydrology

[1] A sequential, geostatistical inverse approach was developed for electrical resistivity tomography (ERT). Unlike most ERT inverse approaches, this new approach allows inclusion of our prior knowledge of general geological structures of an area and point electrical resistivity measurements to constrain the estimate of the electrical resistivity field. This approach also permits sequential inclusion of different data sets, mimicking the ERT data collection scheme commonly employed in the field survey. Furthermore, using the conditional variance concept, the inverse model quantifies uncertainty of the estimate caused by spatial variability and measurement errors. Using this approach, numerical experiments were conducted to demonstrate the effects of bedding orientation on ERT surveys and to show both the usefulness and uncertainty associated with the inverse approach for delineating the electrical resistivity distribution using down-hole ERT arrays. A statistical analysis was subsequently undertaken to explore the effects of spatial variability of the electrical resistivity-moisture relation on the interpretation of the change in water content in the vadose zone, using the change in electrical resistivity. Core samples were collected from a field site to investigate the spatial variability of the electrical resistivity-moisture relation. Numerical experiments were subsequently conducted to illustrate how the spatially varying relations affect the level of uncertainty in the interpretation of change of moisture content based on the estimated change in electrical resistivity. Other possible complications are also discussed. INDEX TERMS: 0903 Exploration Geophysics: Computational methods, potential fields; 1869 Hydrology: Stochastic processes; 1875 Hydrology: Unsaturated zone; 1866 Hydrology: Soil moisture; 3260 Mathematical Geophysics: Inverse theory; KEYWORDS: geostatistical inverse model, electrical resistivity tomography, vadose zone, resistivitymoisture relation, spatial variability, sequential/successive linear estimator

[1]  D. Oldenburg,et al.  Estimating depth of investigation in DC resistivity and IP surveys , 1999 .

[2]  R. Srivastava,et al.  A three-dimensional numerical model for water flow and transport of chemically reactive solute through porous media under variably saturated conditions , 1992 .

[3]  P. Kitanidis,et al.  An Application of the Geostatistical Approach to the Inverse Problem in Two-Dimensional Groundwater Modeling , 1984 .

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

[5]  A. Sato,et al.  Three‐dimensional spatial and temporal monitoring of soil water content using electrical resistivity tomography , 2001 .

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

[7]  Sean Hatch,et al.  Joint conditional simulations and the spectral approach for flow modeling , 1994 .

[8]  T.-C. Jim Yeh,et al.  Estimation of co-conditional moments of transmissivity, hydraulic head, and velocity fields , 1998 .

[9]  P. C. Shah,et al.  Reservoir History Matching by Bayesian Estimation , 1976 .

[10]  T. Yeh,et al.  Hydraulic tomography: Development of a new aquifer test method , 2000 .

[11]  T.-C. Jim Yeh,et al.  An inverse model for three‐dimensional flow in variably saturated porous media , 2000 .

[12]  Prem V. Sharma.  Environmental and Engineering Geophysics by Prem V. Sharma , 1997 .

[13]  E. G. Vomvoris,et al.  A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one‐dimensional simulations , 1983 .

[14]  J. A. Vargas-Guzmán,et al.  Sequential kriging and cokriging: Two powerful geostatistical approaches , 1999 .

[15]  A. Gutjahr,et al.  Stochastic Analysis of Unsaturated Flow in Heterogeneous Soils: 2. Statistically Anisotropic Media With Variable α , 1985 .

[16]  N. Sun.  Inverse problems in groundwater modeling , 1994 .

[17]  Rosemary Knight,et al.  Hysteresis in the electrical resistivity of partially saturated sandstones , 1991 .

[18]  T. Yeh,et al.  A geostatistically based inverse model for three-dimensional variably saturated flow , 1998 .

[19]  Allan L. Gutjahr,et al.  Stochastic Analysis of Unsaturated Flow in Heterogeneous Soils: 3. Observations and Applications , 1985 .

[20]  Douglas W. Oldenburg,et al.  Inversion of 3-D DC resistivity data using an approximate inverse mapping , 1994 .

[21]  Allan L. Gutjahr,et al.  Stochastic Analysis of Unsaturated Flow in Heterogeneous Soils: 1. Statistically Isotropic Media , 1985 .

[22]  Steven F. Carle,et al.  High-Performance Computational And Geostatistical Experiments For Testing The Capabilities Of 3-D Electrical Resis’Iiance Tomography , 1999 .

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

[24]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[25]  Minghui Jin,et al.  AN ITERATIVE STOCHASTIC INVERSE METHOD: CONDITIONAL EFFECTIVE TRANSMISSIVITY AND HYDRAULIC HEAD FIELDS , 1995 .

[26]  D. Oldenburg,et al.  Incorporating geological dip information into geophysical inversions , 2000 .

[27]  Allan L. Gutjahr,et al.  An Iterative Cokriging‐Like Technique for Ground‐Water Flow Modeling , 1995 .

[28]  P. Kitanidis.  Comment on “A reassessment of the groundwater inverse problem” by D. McLaughlin and L. R. Townley , 1997 .

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

[30]  Robert G. Ellis,et al.  The pole-pole 3-D Dc-resistivity inverse problem: a conjugategradient approach , 1994 .

[31]  Shuyun Liu,et al.  Effectiveness of hydraulic tomography: Sandbox experiments , 2001 .

[32]  T.-C. Jim Yeh,et al.  An iterative geostatistical inverse method for steady flow in the vadose zone , 1996 .

[33]  L. Gelhar.  Stochastic Subsurface Hydrology , 1992 .