Solute transport in three-dimensional heterogeneous media with a Gaussian covariance of log hydraulic conductivity

The analytical solutions for the velocity covariance, u(ij), one particle displacement covariance X-ij, and the macrodispersivity tensor alpha(ij) defined as (0.5/mu)(dX(ij)/dt), were derived in three-dimensional heterogeneous media. A Gaussian covariance function (GCF) of logarithmic hydraulic conductivity, log K, Was used, assuming uniform mean flow and first-order approximation in log-conductivity variance, where mu is the magnitude of the mean flow velocity mu. Based on these solutions, the time-dependent ensemble averages of the second spatial moments, Z(ij) = [A(ij)] - A(ij) (0) = X-ij - R-ij and the effective dispersivity tensor gamma(ij), defined as (0.5/mu)(d[A(ij)]/dt), were evaluated for a finite line source either normal or parallel to mu, where A(ij)(0) is the initial value of the second spatial moment of a plume, A(ij), and R-ij is the plume centroid covariance. The results obtained in this study were compared with previous results for an exponential covariance function (ECF). It was found that in a stationary log K field the spreading of a solute plume depends not only on the variance and integral scale of the log K field but also on the shape of its covariance function. The more strongly correlated the hydraulic conductivities at short separation distances are, the faster Z(ii) and gamma(ii) grow at early time. Also, the earlier that y(ii) approaches its asymptote or peak, and the higher the peak is, the larger the negative transverse dispersivity. More importantly, the ergodic limits for GCF are reached faster than those for ECF, as the initial size of a plume increases. The ergodic limit X-11 for GCF is slightly larger than that for ECF, but X-22 and X-33 are significantly smaller than those for ECF even though the asymptotic alpha(ii) for GCF is the same with that for ECF.