Analytical solution to transport in three-dimensional heterogeneous well capture zones.

Solute transport is investigated in a heterogeneous aquifer for combined natural-gradient and well flows. The heterogeneity is associated with the spatially varying hydraulic conductivity K(x, y, z), which is modelled as a log-normal stationary-random function. As such, the conductivity distribution is characterized by four parameters: the arithmetic mean K(A), the variance sigma(Y)(2) (Y=lnK), the horizontal integral scale I of the axisymmetric log-conductivity autocorrelation and the anisotropy ratio e=I(v)/I (I(v) is the vertical integral scale). The well fully penetrates an aquifer of constant thickness B and has given constant discharge QB, while the background aquifer flow is driven by an uniform mean head-gradient, - J. Therefore, for a medium of homogeneous conductivity K(A), the steady-state capture zone has a width 2L=Q/(K(A)|J|) far from the well (herein the term capture zone is used to refer both to the zone from which water is captured by a pumping well and the zone that captures fluid from an injecting well). The main aim is to determine the mean concentration as a function of time in fluid recovered by a pumping well or in a control volume of the aquifer that captures fluid from an injecting well. Relatively simple solutions to these complex problems are achieved by adopting a few assumptions: a thick aquifer B>>I(v) of large horizontal extent (so that boundary effects may be neglected), weak heterogeneity sigma(Y)(2)<1, a highly anisotropic formation e<0.2 and neglect of pore-scale dispersion. Transport is analyzed to the first-order in sigma(Y)(2) in terms of the travel time of particles moving from or towards the well along the steady streamlines within the capture zone. Travel-time mean and variance to any point are computed by two quadratures for an exponential log-conductivity two-point covariance. Spreading is reflected by the variance value, which increases with sigma(Y)(2) and I/L. For illustration, the procedure is applied to two particular cases. In the first one, a well continuously injects water at constant concentration. The mean concentration as function of time for different values of the controlling parameters sigma(Y)(2) and I/L is determined within control volumes surrounding the well or in piezometers. In the second case, a solute plume, initially occupying a finite volume Omega(0), is drawn towards a pumping well. The expected solute recovery by the well as a function of time is determined in terms of the previous controlling parameters as well as the location and extent of Omega(0). The methodology is tested against a full three-dimensional simulation of a multi-well forced-gradient flow field test ([Lemke, L., W.B. II, Abriola, L., Goovaerts, P., 2004. Matching solute breakthrough with deterministic and stochastic aquifer models. Ground Water 42], SGS simulations). Although the flow and transport conditions are more complex than the ones pertinent to capture zones in uniform background flow, it was found that after proper adaptation the methodology led to results for the breakthrough curve in good agreement with a full three-dimensional simulation of flow and transport.

Analytical model of solute transport by unsteady unsaturated gravitational infiltration.

Penetration of reactive solute into a soil during a cycle of water infiltration and redistribution is investigated by deriving analytical closed form solutions for fluid flux, moisture content and contaminant concentration. The solution is developed for gravitational flow and advective transport and is applied to two scenarios of solute applications encountered in the applications: a finite pulse of solute dissolved in irrigation water and an instantaneous pulse broadcasted onto the soil surface. Through comparison to simulations of Richards' flow, capillary suction is shown to have contrasting effects on the upper and lower boundaries of the fluid pulse, speeding penetration of the wetting front and reducing the rate of drying. This leads to agreement between the analytical and numerical solutions for typical field and experimental conditions. The analytical solution is further incorporated into a stochastic column model of flow and transport to compute mean solute concentration in a heterogeneous field. An unusual phenomenon of plume contraction is observed at long times of solute propagation during the drying stage. The mean concentration profiles match those of the Monte-Carlo simulations for capillary length scales typical of sandy soils.