A computationally practical method for stochastic groundwater modeling

Shu-Guang Li,Dennis McLaughlin,Hua-Sheng Liao
DOI: https://doi.org/10.1016/j.advwatres.2003.08.003
IF: 5.361
2003-01-01
Advances in Water Resources
Abstract:Stochastic theories of subsurface flow and transport have changed the way we think about heterogeneity but have not had much impact on practical groundwater modeling. Most numerical models still provide no information on prediction uncertainty. This gap between theory and practice is due largely to the excessive computational demands of available numerical methods for solving stochastic problems. The two primary alternatives, Monte Carlo and classical perturbation methods, both require solving large numbers of equations on fine grid spacings that are smaller than the log conductivity correlation length, no matter how smooth the mean flow and variance dynamics of interest are. This results in computational times that are many orders-of-magnitude greater than those required for conventional deterministic simulations. In addition, classical perturbation methods need much more real memory than comparable deterministic simulations. In this paper, we present an innovative approach for stochastic groundwater modeling––one that finally allows taking advantage of the scale disparity between the mean distributions and small-scale process. The new approach, which is called the nonstationary spectral method, implicitly separates the problem into rapidly varying random (stationary) part and a slowly varying deterministic (nonstationary) part and only the latter needs to be solved numerically. A general analysis and a specific numerical example both demonstrate that the nonstationary spectral method is much more efficient than the available alternatives. The methodological innovation dramatically increases the size and range of the problems that can be solved using stochastic methods and represents a major step forward in our effort to make stochastic uncertainty analysis a routine part of groundwater modeling.
What problem does this paper attempt to address?