Anomalous transport of colloids and solutes in a shear zone

Transport experiments with colloids and radionuclides in a shear zone were conducted during the Colloid and Radionuclide Retardation experiment (CRR) at Nagra's Grimsel Test Site. Breakthrough curves of bentonite colloids and uranine, a non-sorbing solute, were measured in an asymmetric dipole flow field. The colloid breakthrough is earlier than that of uranine. Both breakthrough curves show anomalously long late time tails and the slope of the late time tails for the colloids is slightly higher. Anomalous late time tails are commonly associated with matrix diffusion processes; the diffusive interaction of solutes transported in open channels with the adjacent porous rock matrix or zones of stagnant water. The breakthrough curves for different colloid size classes are very similar and show no signs of fractionation due to their (size-dependent) diffusivity. It is proposed that tailing of the colloids is mainly caused by the structure of the flow field and that for the colloid transport, matrix diffusion is of minor importance. This has consequences for the interpretation of the uranine breakthrough. Comparisons of experimental results with numerical studies and with the evaluation of the colloid breakthrough with continuous time random theory imply that the tailing in the conservative solute breakthrough in this shear zone is not only caused by matrix diffusion. Part of the tailing can be attributed to advective transport in fracture networks and advection in low velocity regions. Models based on the advection-dispersion equation and matrix diffusion do not properly describe the temporal and spatial evolution of colloid and solute transport in such systems with a consistent set of parameters.     

The effect of a biofilm on solute diffusion in fractured porous media

At sites in fractured rock where contamination has been exposed to the rock matrix for extended periods of time, the amount of contaminant mass residing in the matrix can be considerable. Even though it may be possible to diminish concentrations by the advection of clean water through the fracture features, back diffusion from mass held in the matrix will lead to a continuing source of contamination. In such an event, the development of a biofilm (a thin film of microbial mass) on the wall of the fractures may act to limit or prevent the back diffusion process. The objective of this preliminary study is to explore the influence imparted by the presence of a biofilm on the process of matrix diffusion. The investigation was conducted using radial diffusion cells constructed from rock core in which biofilm growth was stimulated in a central reservoir. Once biofilms were developed, forward diffusion experiments were conducted in which a conservative solute migrated from the central reservoir into the intact rock sample. Diffusion experiments were performed in a total of 11 diffusion cell pairs where biofilm growth was stimulated in one member of the pair and inhibited in the other. The effect of the presence of a biofilm on tracer diffusion was determined by comparison of the diffusion curves produced by each cell pair. A semi-analytical model that accounts for the presence of a biofilm was used to investigate the effect of the biofilm on mass transfer due to changes in the effective porosity, effective diffusion coefficient, and the depth of penetration of the biofilm into the intact rock. The results show that the biofilm acted to plug the rock matrix, rather than forming a discrete layer on the reservoir surface. The reduction in effective porosity due to the biofilm ranged from 6% to 52% with the majority of the samples in the 30% to 50% range. Based on the present results, with more efficient biofilm stimulation, it is reasonable to assume that a more complete plugging of the microcrack porosity might be possible, leaving a much thicker and efficient barrier than could be achieved via a surface biofilm.     

Estimating the hydraulic properties of the fracture network in a sandstone aquifer

The potential hydraulic behaviour of the fracture network in a major Triassic sandstone aquifer in the UK has been evaluated. The properties of the fracture network were determined using results from detailed scan line surveys at 10 sites, television and geophysical borehole logging, and packer testing. Six sets of discontinuities common to all sites were identified and statistically characterised (dip, strike, orientation, density, size, and estimated transmissivity). A discrete fracture network model was then used stochastically to investigate the properties of the network. In general, the network is poorly connected: it is estimated that 9% of the discontinuities intersecting boreholes are transmissive. The hydraulic behaviour of the network is generally dominated by one sub-horizontal bedding plane fracture set, although when present, a relatively infrequent north-south striking, sub-vertical set modifies the bulk flow properties significantly. Ignoring this latter set, the network's minimum representative volume is about 35 x 35 x 35 m. The upscaled permeability is anisotropic, being typically 23 times greater in the horizontal than in the vertical. Tortuosity in the north-south direction is around 1.6.     

An analytical model for the interpretation of pulse injection experiments performed for testing the spatial variability of clay formations

This paper presents an analytical model to describe pulse injection experiments. This model solves the advection-diffusion equation while taking into account back diffusion from the clay core to the inlet and from the outlet to the clay core. In most analytical models, back diffusion is neglected. For sufficiently high Péclet numbers, this is a good approximation. However, in experiments where the Péclet number is low, back diffusion is important and must be taken into account. An additional advantage of the present model is that both concentration and flux are conserved at the inlet and at the outlet of the clay core. This model is used to fit pulse injection experiments with iodide and tritiated water (HTO) in clay cores. The (new) model is required for fitting the experimental results since in clay layers advection is very slow leading to a low Péclet number. The experiments are performed on clay cores taken from different depths from the Boom Clay and the Ypres Clay layer under the site of the nuclear power plant of Doel (Belgium). The quality of all fits is excellent and the obtained parameter values are coherent. For HTO, the fitted value for the diffusion accessible porosity is consistent with measurements of the water content in Ypres Clay cores. In both types of clays, the apparent diffusion coefficient at zero flow is between 10(-10) and 2 x 10(-10) m(2)/s for iodide and between 2 x 10(-10) and 3 x 10(-10) m(2)/s for HTO. The dispersion length is in the order of 10(-3) m. The average value for the diffusion accessible porosity is between 0.35 and 0.4 for HTO and between 0.2 and 0.25 for iodide.     

Gas/particle partitioning of atmospheric PCDD/Fs in a satellite town in Eastern China

Gas/particle partitioning of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) in ambient air was investigated in a satellite town in Eastern China from April 2007 to January 2008 comprehending large temperature variations (from 3 to 34 °C, daily average). Molecular weight, molecular structure and ambient temperatures are the three major factors that govern the gas/particle partitioning of atmospheric PCDD/Fs throughout the year. Generally, good agreements were obtained (except for winter) between measured particulate fractions and theoretical estimates of both the Junge–Pankow adsorption model and Harner Bidleman absorption model using different sets of subcooled liquid vapor pressure and octanol–air partition coefficient (K_oa), respectively. Models utilizing estimates, derived from gas chromatographic retention indices (GC-RIs), are more accurate than that of entropy-based. Moreover, during winter, the K_oa-based model using the GC-RIs approach performs better on lower chlorinated PCDD/Fs than that of -based. Furthermore, possible sources of mismatch between measured and predicted values in winter (3–7 °C) were discussed. Gas adsorption artifact was demonstrated to be of minor importance for the phenomena observed. On the other hand, large deviations of slopes (m_r) and intercepts (b_r) in logK_p vs. plots from theoretical values are observed in the literature data and these are found to be linearly correlated with ambient temperatures (P<0.001) in this study. This indicates that the non-equilibrium partitioning of PCDD/Fs in winter may be significantly influenced by the colder temperatures that may have slowed down the exchange between gaseous and particulate fractions.