No abstract

No abstract

We introduce a method for identifying the transverse dispersion coefficient in laboratory experiments based on the analytical solution of a pulse injection of a nonreactive solute in a soil column (cylindrical geometry) packed with a homogeneous porous medium. This method takes into account the effect of boundary conditions such as no flux on the column perimeter, and it does not need a priori knowledge of the longitudinal dispersion coefficient. Numerical applications of the method show that it is stable and robust and that the results are reasonably in accordance with those found using the classical maximum likelihood method.

Despite the availability of numerical models, interest in analytical solutions of multidimensional advection-dispersion systems remains high. Such models are commonly used for performing Tier I risk analysis and are embedded in many regulatory frameworks dealing with groundwater contamination. In this work, we develop a closed-form solution of the three-dimensional advectiondispersion equation with exponential source decay, first-order reaction, and retardation, and present an approach based on some ease of use diagrams to compare it with the integral open form solution and with earlier versions of the closed-form solution. The comparison approach focuses on the relative differences associated with source decay and the effect of simulation time. The analysis of concentration contours, longitudinal sections, and transverse sections confirms that the closed-form solutions studied can be used with acceptable approximation in the central area of a plume bound transversely within the source width, both behind and beyond the advective front and for concentration values up to two orders of magnitude less than the initial source concentration. As the proposed closed-form model can be evaluated without nested numerical computations and with simple mathematical functions, it can be very useful in risk assessment procedures.

[1] We analyze the influence of the formation structure upon transport of nonreactive and kinetically sorbing solutes in heterogeneous porous formations. Formation structure is represented through a bimodal distribution of hydraulic conductivity with two scales of spatial variability, which is a simplified, yet realistic, model of spatial variability. Kinetic sorption is modeled using a linear reversible rate expression. We derive new closed-form solutions of the second ergodic spatial moments of a kinetically sorbing solute. Although these solutions are obtained through a perturbation expansion in term of logconductivity variance, they are not limited to any particular model of spatial variability. We show that the formation structure, as described by the bimodal model, has a significant impact on plume moments of a non-reactive tracer at moderate contrast between mean hydraulic conductivity of the two hydrofacies, when the two scales of variability have comparable weight. On the contrary, the impact diminishes when bimodality weakens at both low and high contrast. As expected, the large time asymptotic longitudinal macrodispersion coefficient is influenced by the structure only through the variance and integral scale of log-conductivity, and thus the result for the bimodal model is identical to that of a unimodal Gaussian field sharing the same second-order statistics. The impact of kinetic sorption depends upon the plume travel time since injection. At early times, sorption kinetics tends to accentuate the difference between the longitudinal macrodispersivity in bimodal and unimodal fields, while at later times the difference is minimal because overall macrodispersivity is dominated by the spreading caused by nonequilibrium sorption. In a bimodal formation transverse macrodispersivity peaks later and declines to zero slower than in a Gaussian unimodal field with the same second-order statistics. These results demonstrate that even a very simplified facies model, such as the bimodal model relying on a statistical distribution of high-and low-conductivity zones, has important impacts upon global parameters describing the migration of nonreactive and reactive tracers. Consequently, larger effects are expected for quantities such as early arrival times and BTC tailing which are more directly influenced by the spatial patterns of high hydraulic conductivity than global parameters, such as macrodispersivity. Furthermore, these solutions contribute to bridge the gap between stochastic modeling and applications, since bimodal fields can be obtained to honor soft data typically available in applications, for example through conditional simulations using stratigraphic information from well logs.Citation: Massabó, M., A. Bellin, and A. J. Valocchi (2008), Spatial moments analysis of kinetically sorbing solutes in aquifer with bimodal permeability distribution, Water Resour. Res., 44, W09424,