Premium
Stochastic analysis of nonstationary subsurface solute transport: 1. Unconditional moments
Author(s) -
Graham Wendy,
McLaughlin Dennis
Publication year - 1989
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/wr025i002p00215
Subject(s) - monte carlo method , covariance , mathematics , perturbation (astronomy) , partial differential equation , standard deviation , advection , statistical physics , mathematical analysis , statistics , physics , thermodynamics , quantum mechanics
This paper applies stochastic methods to the analysis and prediction of solute transport in heterogeneous saturated porous media. Partial differential equations for three unconditional ensemble moments (the concentration mean, concentration covariance, and velocity concentration cross covariance) are derived by applying perturbation techniques to the governing transport equation for a conservative solute. Concentration uncertainty is assumed to be the result of unmodeled small‐scale fluctuations in a steady state velocity field. The moment expressions, which describe how each moment evolves over time and space, resemble the classic deterministic advection‐dispersion equation and can be solved using similar methods. A solution procedure based on a Galerkin finite element algorithm is illustrated with a hypothetical two‐dimensional example. For this example the required steady state velocity statistics are obtained from an infinite domain spectral solution of the stochastic groundwater flow equation. The perturbation solution is shown to reproduce the statistics obtained from a Monte Carlo simulation quite well for a natural log conductivity standard deviation of 0.5 and moderately well for a natural log conductivity standard deviation of 1.0. The computational effort required for a perturbation solution is significantly less than that required for a Monte Carlo solution of acceptable accuracy. Sensitivity analyses conducted with the perturbation approach provide qualitative confirmation of a number of results obtained by other investigators for more restrictive special cases.