Stochastic inverse mapping of hydraulic conductivity and sorption partitioning coefficient fields conditioning on nonreactive and reactive tracer test data

[1] A three-dimensional geostatistically based iterative inverse method is presented for mapping spatial distributions of the hydraulic conductivity and sorption partitioning coefficient fields by sequential conditioning on both nonreactive and reactive tracer breakthrough data. A streamline-based semianalytical simulator is adopted to simulate chemical movement in a physically and chemically heterogeneous field and serves as the forward modeling. In this study, both the hydraulic conductivity and sorption coefficient are assumed to be random spatial variables. Within the framework of the streamline-based simulator an efficient semianalytical method is developed to calculate sensitivity coefficients of the reactive chemical concentration with respect to the changes of conductivity and sorption coefficient. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by a synthetic tracer test conducted within an aquifer with distinct spatial features of physical and chemical heterogeneities. The study results indicate that our iterative stochastic inverse method is able to identify and reproduce the large-scale physical and chemical heterogeneity features. The conditional study on the geostatistical distributions of conductivity and sorption coefficient are expected to significantly reduce prediction uncertainties for solute transport in the medium.

[1]  Bill X. Hu,et al.  Nonlocal reactive transport in heterogeneous dual‐porosity media with rate‐limited sorption and interregional mass diffusion , 2001 .

[2]  Akhil Datta-Gupta,et al.  On the Sensitivity and Spatial Resolution of Transient Pressure and Tracer Data For Heterogeneity Characterization , 1997 .

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

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

[5]  R. L. Naff,et al.  Nonreactive and reactive solute transport in three-dimensional heterogeneous porous media: Mean displacement, plume spreading, and uncertainty , 1994 .

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

[7]  Akhil Datta-Gupta,et al.  A semianalytic approach to tracer flow modeling in heterogeneous permeable media , 1995 .

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

[9]  S. P. Garabedian,et al.  Large-scale dispersive transport in aquifers : field experiments and reactive transport theory , 1987 .

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

[11]  G. Dagan,et al.  Conditional estimation of solute travel time in heterogeneous formations: Impact of transmissivity measurements , 1992 .

[12]  Y. Rubin.  Transport in heterogeneous porous media: Prediction and uncertainty , 1991 .

[13]  G. Dagan,et al.  Stochastic identification of transmissivity and effective recharge in steady groundwater flow: 2. Case study , 1987 .

[14]  A. Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data—I. Theory , 1997 .

[15]  G. Marsily,et al.  Cokriged estimation of aquifer transmissivity as an indirect solution of the inverse problem: A practical approach , 1993 .

[16]  Z. J. Kabala,et al.  A stochastic model of reactive solute transport with time‐varying velocity in a heterogeneous aquifer , 1991 .

[17]  Xian-Huan Wen,et al.  Significance of conditioning to piezometric head data for predictions of mass transport in groundwater modeling , 1996 .

[18]  A. S. Cullick,et al.  Construction of geostatistical aquifer models integrating dynamic flow and tracer data using inverse technique , 2002 .

[19]  D. A. Zimmerman,et al.  A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow , 1998 .

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

[21]  J. H. Cushman,et al.  Nonlocal reactive transport with physical, chemical, and biological heterogeneity , 1997 .

[22]  E. Sudicky,et al.  Spatial Variability of Strontium Distribution Coefficients and Their Correlation With Hydraulic Conductivity in the Canadian Forces Base Borden Aquifer , 1991 .

[23]  A. D. Gupta,et al.  Stochastic modeling of spatial heterogeneities conditioned to hydraulic and tracer tests , 1994 .

[24]  Yoram Rubin,et al.  Prediction of tracer plume migration in disordered porous media by the method of conditional probabilities , 1991 .

[25]  Charles F. Harvey,et al.  Mapping Hydraulic Conductivity: Sequential Conditioning with Measurements of Solute Arrival Time, Hydraulic Head, and Local Conductivity , 1995 .

[26]  Allan L. Gutjahr,et al.  Co-kriging for stochastic flow models , 1989 .

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

[28]  P. Kitanidis,et al.  Maximum likelihood parameter estimation of hydrologic spatial processes by the Gauss-Newton method , 1985 .

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

[30]  Allan L. Gutjahr,et al.  Cross‐correlated random field generation with the direct Fourier Transform Method , 1993 .

[31]  Andrea Rinaldo,et al.  Linear equilibrium adsorbing solute transport in physically and chemically heterogeneous porous formations: 2. Numerical results , 1993 .

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

[33]  M. Marietta,et al.  Pilot Point Methodology for Automated Calibration of an Ensemble of conditionally Simulated Transmissivity Fields: 1. Theory and Computational Experiments , 1995 .

[34]  L. Gelhar,et al.  Stochastic analysis of sorption macrokinetics in heterogeneous aquifers , 1996 .

[35]  Andrés Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric head data—3. Application to the Culebra formation at the Waste Isolation Pilot Plan (WIPP), New Mexico, USA , 1998 .

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

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

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

[39]  Don W. Vasco,et al.  Inversion of pressure observations: an integral formulation , 2001 .

[40]  G. Dagan,et al.  Contaminant transport in aquifers with spatially variable hydraulic and sorption properties , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[41]  J. H. Cushman,et al.  Inverse methods for subsurface flow: A critical review of stochastic techniques , 1990 .

[42]  Peter K. Kitanidis,et al.  Sensitivity of temporal moments calculated by the adjoint-state method and joint inversing of head and tracer data , 2000 .

[43]  Andrés Sahuquillo,et al.  Joint simulation of transmissivity and storativity fields conditional to steady-state and transient hydraulic head data , 1999 .

[44]  Andrés Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data 2. Demonstration on a synthetic aquifer , 1997 .

[45]  G. Pope,et al.  Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing , 1995 .

[46]  A Comparison of Several Methods for the Solution of the Inverse Problem in Two‐Dimensional Steady State Groundwater Flow Modeling , 1986 .

[47]  R. Allen‐King,et al.  Characterizing the heterogeneity and correlation of perchloroethene sorption and hydraulic conductivity using a facies‐based approach , 1998 .

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

[49]  M. G. Marietta,et al.  Pilot Point Methodology for Automated Calibration of an Ensemble of Conditionally Simulated Transmissivity Fields: 2. Application , 1995 .

[50]  A. S. Cullick,et al.  A program to create permeability fields that honor single-phase flow rate and pressure data , 1999 .

[51]  Akhil Datta-Gupta,et al.  Integrating field production history in stochastic reservoir characterization , 1997 .

[52]  Jesús Carrera,et al.  State of the Art of the Inverse Problem Applied to the Flow and Solute Transport Equations , 1988 .