Necessary conditions for inverse modeling of flow through variably saturated porous media

Non-unique solutions of inverse problems arise from a lack of information that satisfies necessary conditions for the problem to be well defined. This paper investigates these conditions for inverse modeling of water flow through multi-dimensional variably saturated porous media. It shows that in order to obtain a unique estimate of hydraulic parameters, along each streamline of the flow field (1) spatial and temporal head observations must be given; (2) the number of spatial and temporal head observations required should be greater or equal to the number of unknown parameters; (3) the flux boundary condition or the pumping rate of a well must be specified for the homogeneous case and both boundary flux and pumping rate are a must for the heterogeneous case; (4) head observations must encompass both saturated and unsaturated conditions, and the functional relationships for unsaturated hydraulic conductivity/pressure head and for the moisture retention should be given, and (5) the residual water content value also need to be specified a priori or water content measurements are needed for the estimation of the saturated water content. For field problems, these necessary conditions can be collected or estimated but likely involve uncertainty. While the problems become well defined and have unique solutions, the solutions likely will be uncertain. Because of this uncertainty, stochastic approaches are deemed to be appropriate for inverse problems as they are for forward problems to address uncertainty. Nevertheless, knowledge of these necessary conditions is critical to reduce uncertainty in both characterization of the vadose zone and the aquifer, and prediction of water flow and solute migration in the subsurface.

[1]  Peter Dietrich,et al.  A travel time based hydraulic tomographic approach , 2003 .

[2]  Tian-Chyi J. Yeh,et al.  Characterization of aquifer heterogeneity using transient hydraulic tomography , 2004 .

[3]  P. Kitanidis Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[4]  R. W. Nelson,et al.  In‐place measurement of permeability in heterogeneous media: 1. Theory of a proposed method , 1960 .

[5]  C. Axness,et al.  Three‐dimensional stochastic analysis of macrodispersion in aquifers , 1983 .

[6]  T. Harter,et al.  A Numerical Model for Water Flow and Chemical Transport in Variably Saturated Porous Media , 1993 .

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

[8]  P. K. Kitanidis,et al.  Large‐scale inverse modeling with an application in hydraulic tomography , 2011 .

[9]  E. Sudicky,et al.  A view toward the future of subsurface characterization: CAT scanning groundwater basins , 2008 .

[10]  Daniel M. Tartakovsky,et al.  Conditional moment analysis of steady state unsaturated flow in bounded, randomly heterogeneous soils , 2002 .

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

[12]  James J. Butler,et al.  Steady shape analysis of tomographic pumping tests for characterization of aquifer heterogeneities , 2002 .

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

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

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

[16]  W. Li,et al.  Two‐dimensional characterization of hydraulic heterogeneity by multiple pumping tests , 2007 .

[17]  A Revil,et al.  A Potential‐Based Inversion of Unconfined Steady‐State Hydraulic Tomography , 2009, Ground water.

[18]  T. Yeh,et al.  Cokriging estimation of the conductivity field under variably saturated flow conditions , 1999 .

[19]  E. Kreyszig,et al.  Advanced Engineering Mathematics. , 1974 .

[20]  C. R. Dietrich,et al.  A stability analysis of the geostatistical approach to aquifer transmissivity identification , 1989 .

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

[22]  Jirka Simunek,et al.  Indirect estimation of soil thermal properties and water flux using heat pulse probe measurements: Geometry and dispersion effects , 2002 .

[23]  G. Newsam,et al.  Sufficient conditions for identifying transmissivity in a confined aquifer , 1990 .

[24]  Geoffrey C. Bohling,et al.  A field assessment of the value of steady shape hydraulic tomography for characterization of aquifer heterogeneities , 2007 .

[25]  S. P. Neuman,et al.  Improved forward and inverse analyses of saturated‐unsaturated flow toward a well in a compressible unconfined aquifer , 2010 .

[26]  S. Mathias,et al.  Linearized Richards' equation approach to pumping test analysis in compressible aquifers , 2006 .

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

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

[29]  J. Hadamard Sur les problemes aux derive espartielles et leur signification physique , 1902 .

[30]  Peter Dietrich,et al.  Identification of the permeability distribution in soil by hydraulic tomography , 1995 .

[31]  T.-C. Jim Yeh,et al.  A Geostatistical Inverse Method for Variably Saturated Flow in the Vadose Zone , 1995 .

[32]  William W.-G. Yeh,et al.  A stochastic inverse solution for transient groundwater flow: Parameter identification and reliability analysis , 1992 .

[33]  Steven J. Berg,et al.  Field Study of Subsurface Heterogeneity with Steady‐State Hydraulic Tomography , 2013, Ground water.

[34]  Peter Droogers,et al.  Inverse method to determine soil hydraulic functions from multistep outflow experiments. , 1994 .

[35]  T. Yeh,et al.  Stochastic inversion of pneumatic cross-hole tests and barometric pressure fluctuations in heterogeneous unsaturated formations , 2008 .

[36]  Kristopher L. Kuhlman,et al.  Basin‐Scale Transmissivity and Storativity Estimation Using Hydraulic Tomography , 2008, Ground water.

[37]  Peter Dietrich,et al.  A field assessment of high‐resolution aquifer characterization based on hydraulic travel time and hydraulic attenuation tomography , 2011 .

[38]  J. M. Shafer,et al.  Assessment of Uncertainty in Time‐Related Capture Zones Using Conditional Simulation of Hydraulic Conductivity , 1991 .

[39]  Dennis McLaughlin,et al.  Stochastic analysis of nonstationary subsurface solute transport: 1. Unconditional moments , 1989 .

[40]  Junfeng Zhu,et al.  Sequential aquifer tests at a well field, Montalto Uffugo Scalo, Italy , 2007 .

[41]  W. R. Gardner SOME STEADY‐STATE SOLUTIONS OF THE UNSATURATED MOISTURE FLOW EQUATION WITH APPLICATION TO EVAPORATION FROM A WATER TABLE , 1958 .

[42]  Jan W. Hopmans,et al.  Optimization of Hydraulic Functions from Transient Outflow and Soil Water Pressure Data , 1993 .

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

[44]  Franklin W. Schwartz,et al.  mass transport: 2. Analysis of uncertainty in prediction , 1981 .

[45]  S. P. Neuman,et al.  Conditional stochastic averaging of steady state unsaturated flow by means of Kirchhoff Transformation , 1999 .

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

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

[48]  T.-C. Jim Yeh,et al.  Stochastic modelling of groundwater flow and solute transport in aquifers , 1992 .

[49]  J. Berryman Stable iterative reconstruction algorithm for nonlinear traveltime tomography , 1990 .

[50]  T. Yeh,et al.  Effective unsaturated hydraulic conductivity of layered sands , 1990 .

[51]  Jet-Chau Wen,et al.  A simultaneous successive linear estimator and a guide for hydraulic tomography analysis , 2009 .

[52]  G. Dagan Stochastic Modeling of Groundwater Flow by Unconditional and Conditional Probabilities: The Inverse Problem , 1985 .

[53]  Tian-Chyi J. Yeh,et al.  Robustness of joint interpretation of sequential pumping tests: Numerical and field experiments , 2011 .

[54]  S. P. Neuman,et al.  Inverse stochastic moment analysis of steady state flow in randomly heterogeneous media , 2006 .

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

[56]  Peter Droogers,et al.  Inverse Method for Determining Soil Hydraulic Functions from One‐Step Outflow Experiments , 1992 .

[57]  D. McLaughlin,et al.  Stochastic analysis of nonstationary subsurface solute transport: 2. Conditional moments , 1989 .

[58]  W. Nowak,et al.  Geostatistical inference of hydraulic conductivity and dispersivities from hydraulic heads and tracer data , 2006 .

[59]  T. Yeh,et al.  Analysis of hydraulic tomography using temporal moments of drawdown recovery data , 2006 .

[60]  A. Bellin,et al.  A Bayesian approach for inversion of hydraulic tomographic data , 2009 .

[61]  Junfeng Zhu,et al.  Traditional analysis of aquifer tests: Comparing apples to oranges? , 2005 .

[62]  Yefang Jiang,et al.  Transmissivity estimation for highly heterogeneous aquifers: comparison of three methods applied to the Edwards Aquifer, Texas, USA , 2007 .

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

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

[65]  S. Yakowitz,et al.  Instability in aquifer identification: Theory and case studies , 1980 .

[66]  S. P. Neuman,et al.  Three‐dimensional saturated‐unsaturated flow with axial symmetry to a partially penetrating well in a compressible unconfined aquifer , 2007 .

[67]  Willem Bouten,et al.  Toward Improved Identifiability of Soil Hydraulic Parameters: On the Selection of a Suitable Parametric Model , 2003 .

[68]  Peter K. Kitanidis,et al.  An interactive Bayesian geostatistical inverse protocol for hydraulic tomography , 2006 .

[69]  Franklin W. Schwartz,et al.  Mass transport: 3. Role of hydraulic conductivity data in prediction , 1981 .

[70]  A. Moench,et al.  Analytical and numerical analyses of an unconfined aquifer test considering unsaturated zone characteristics , 2008 .

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

[72]  S. Nakagiri,et al.  Identifiability of Spatially-Varying and Constant Parameters in Distributed Systems of Parabolic Type , 1977 .

[73]  Yih-Chi Tan,et al.  Fusion of Active and Passive Hydrologic and Geophysical Tomographic Surveys: The Future of Subsurface Characterization , 2013 .

[74]  W. Press Numerical recipes in Fortran 77 : the art of scientific computing : volume 1 of fortran numerical recipes , 1996 .

[75]  K. Singha,et al.  Subsurface Hydrology: Data Integration for Properties and Processes , 2007 .

[76]  Hiromitsu Saegusa,et al.  Hydraulic tomography in fractured granite: Mizunami Underground Research site, Japan , 2009 .

[77]  S. P. Neuman,et al.  Role of model selection criteria in geostatistical inverse estimation of statistical data‐ and model‐parameters , 2011 .

[78]  T. Harter,et al.  Conditional stochastic analysis of solute transport in heterogeneous, variably saturated soils , 1996 .

[79]  J. P. Delhomme,et al.  Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach , 1979 .

[80]  Walter A. Illman,et al.  Three‐dimensional transient hydraulic tomography in a highly heterogeneous glaciofluvial aquifer‐aquitard system , 2011 .