A stochastic‐conceptual analysis of one‐dimensional groundwater flow in nonuniform homogeneous media

The most realistic representation of a naturally occurring porous medium is a stochastic set of macroscopic elements in which the values of the three basic hydrogeologic parameters (hydraulic conductivity K , compressibility α, and porosity n ) are defined by frequency distributions. A homogeneous formation under this representation is one in which the frequency distributions do not change through space. All soils and geologic formations, even the ones that are homogeneous, show random variations in the values of the hydrogeological parameters through space; that is, they are nonuniform, and a measure of the nonuniformity is provided by the standard deviation of the frequency distributions. If K and α are log normally distributed and n is normally distributed, and if we define Y = log K and C = log α, then the parameters Y, C, and n can be generated from a multivariate normal density function with means μ y , μ c , and μ n , standard deviations σ y , σ c , and σ n , and correlation coefficients ρ yc , ρ yn , and ρ cn The analysis of groundwater flow in nonuniform media requires a stochastic‐conceptual approach in which the effects of stochastic parameter distributions on predicted hydraulic heads are analyzed with the aid of a set of Monte Carlo solutions to the pertinent boundary value problems. In this study, two one‐dimensional saturated flow problems are analyzed: steady state flow between two specified heads and transient consolidation of a clay layer. The primary output is the statistical distribution of hydraulic head ϕ, through space and time, as indicated by the mean values and their standard deviations S ϕ¯ ( x , t ) Results show that the standard deviations of the input hydrogeologic parameters, particularly σ y and σ c , are important index properties; changes in their values lead to different responses for even when the means μ y , μ c , and μ n are fixed. The degree of uncertainty associated with hydraulic head predictions increases as the degree of nonuniformity of the porous medium increases. For large values of σ y and σ c it becomes virtually impossible to obtain meaningful hydraulic head predictions. For transient flow the output distribution of hydraulic head values is almost never normal; in some cases it approaches a uniform distribution. The results of this study throw into question the validity of the hidden assumption that underlies all deterministic groundwater modeling; namely, that it is possible to select a single value for each flow parameter in a homogeneous but nonuniform medium that is somehow representative and hence define an ‘equivalent’ uniform porous medium. For transient flow there may be no way to define an equivalent medium. The fact that nine index parameters rather than three are required to describe a nonuniform geologic formation, the large uncertainties in predicted hydraulic heads for relatively simple flow problems in nonuniform soils, and the contention that there may be no simple way to define an equivalent uniform porous medium all have important implications in the development of groundwater flow theory and in its most fundamental applications.