Stochastic macrodispersion models for gravel aquifers

Statistical information on sedimentary structures determined from natural g ravel deposits in north-eastern Switzerland is used to characterize transpo rt phenomena of aquifers. These aquifers essentially consist of a backgroun d gravel matrix with embedded lenses of different hydraulic properties and internal structure. One of the embedded materials is highly conductive open framework gravel which leads to a;bimodal probability density function of the hydraulic conductivity of the overall gravel material. Based on an anal ysis of this information three approaches to modeling macrodispersion are s uggested. in the first one distinct unconditioned realizations of synthetic block shaped facies type aquifers are generated numerically such that they exhibit the same statistical properties with respect to facies geometry an d hydraulic properties as the deposits themselves. Numerical experiments si mulating saturated flow and transport were subsequently performed with the help of a three-dimensional finite element flow model and a corresponding r andom walk transport model. A total of 100 experiments allowed estimates of apparent, time dependent macrodispersivity values. In the second approach the statistical parameters characterizing the overall gravel deposit were d irectly applied to an analytical unimodal stochastic model of apparent macr odispersivity according to Dagan (1988). For the modeling of the bimodal na ture of the hydraulic conductivity an analytical anisotropic stochastic mod el was developed based on the isotropic model of Rubin (1995). This third a pproach was again applied using the appropriate statistical parameters. Dif ferences in the results of the three models are discussed. They can be attr ibuted mainly to uncertainty in the input parameters, and to the complex se dimentary structure predominant in the natural gravel deposits which were i nvestigated.