Modeling solute diffusion in the presence of pore-scale heterogeneity: method development and an application to the Culebra dolomite member of the Rustler Formation, New Mexico, USA

Previous studies have revealed the presence of pore-scale variability in diffusivity in the Culebra (dolomite) member of the Rustler Formation, NM. In this study, eight laboratory-scale diffusion experiments on five Culebra samples were analyzed using a methodology for modeling solute diffusion through porous media in the presence of multiple matrix diffusivities, Dp. A lognormal distribution of Dp is assumed within each of the lab samples. The estimated standard deviation (sigma d) of ln(Dp) within each sample ranges from 0 to 1, with most values lying between 0.5 and 1. The variability over all samples leads to a combined sigma d in the range of 1.0-1.2, which is consistent with the distribution of independently determined formation factor measurements for similar Culebra samples. A comparison of our estimation results to other rock properties suggests that, at the lab-scale, the geometric mean of Dp increases with bulk porosity and the quantity of macroscopic features such as vugs and fractures. However, sigma d appears to be determined by variability within such macroscopic features and/or by micropore-scale heterogeneity. In addition, comparison of these experiments to those at larger spatial scales suggests that increasing sample volume results in an increase in sigma d.     

Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted Äspö Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175].The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments.While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption K_ds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.     

Solute transport in crystalline rocks at Aspö--I: geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted Asp? Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours-days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption K(d)s are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.     

Radionuclide transport and retardation in rock fracture and crushed rock column experiments

Transport and retardation of non-sorbing tritiated water and chloride and slightly sorbing sodium was studied in Syyry area SY-KR7 mica gneiss, in altered porous tonalite and in fresh tonalite. Experiments were performed using dynamic fracture and crushed rock column methods. Static batch method for sodium was introduced to compare retardation values from static and dynamic experiments. The ¹⁴C-PMMA method was used to study the pore structure of matrices. The pore aperture distribution was evaluated from Hg-porosimetry determinations and the surface areas were determined using the B.E.T. method. The flow characteristics and transport behavior of tracers were interpreted using a numerical compartment model for dispersion. The effect of matrix diffusion was calculated using an analytical solution to the advection-matrix diffusion problem in which surface retardation was taken into account. Radionuclide transport behavior in rock fractures was explained on the basis of rock structure.     

Tracer diffusion from a horizontal fracture into the surrounding matrix: measurement by computed tomography

The vertical diffusion of NaI solution from a horizontal fracture into and within the surrounding matrix was tracked and quantified over time using an artificially fractured chalk core (30x5 cm) and a second-generation X-ray computed tomography (CT) scanner. The different tracer-penetration distances imaged in the matrix above and below the horizontal fracture are indicative of a greater tracer mass penetrating into the lower matrix. The enhanced transport in the matrix below the fracture was related to the Rayleigh-Darcy instability induced by the density differences between the heavier tracer solution in the fracture (1.038) and the distilled water that had initially resided in the matrix. Our observations suggest that below the fracture, the tracer is propagated by an advection-diffusion process that is characterized by both higher rates and higher concentrations relative to its propagation by diffusion above the fracture. The experimental results suggest that the prediction of contaminant migration in a rock intersected by both vertical and horizontal (e.g. along bedding planes) fractures is difficult because of density effects that result in different solute-penetration rates.     

An evaluation of contaminant migration patterns at two waste disposal sites on fractured porous media in terms of the equivalent porous medium (EPM) model

Contamination has occurred many non-indurated and bedrock systems wherein the groundwater flows almost exclusively through a network of connected, open fractures. The matrix surrounding the fractures often possesses porosity which allows contaminant diffusion into the matrix. If the diffusion rates are fast relative to the fracture groundwater velocity, transport effects may be predicted by considering the system to be an equivalent porous medium (EPM). The rapidity with which fracture/immobile-matrix equilibrium is established will be determined in part by the: fracture aperture (2b); interfracture spacing (2B); porosity in the immobile matrix (θ_im); and the matrix diffusion coefficient (D′). Two systems which are characterized by very different values of the above parameters have been studied by our laboratories. At Alkali Lake, Oregon, the EPM approach describes contaminant transport well. At Bayview Park, Ontario, the EPM approach is not appropriate. Several features of the two sites are compared to illustrate the different nature of these two sites. These features include: (1) natural characteristics of the groundwater systems; (2) contaminant distributions; (3) observed transport; and (4) computed fracture/immobile-matrix diffusion times.     

Transverse dispersion of non-reactive tracers in porous media: a new nonlinear relationship to predict dispersion coefficients

Assessing the potential of natural attenuation in groundwater relies on the ability to predict and quantify the processes that occur in contaminant plumes. Transverse dispersion is a significant mass transfer mechanism for mixing of electron acceptors and donors and thus may control the lengths of steady state plumes. Laboratory experiments were carried out using a 2-dimensional acrylic glass tank filled with glass beads, quartz sand and field site material as porous media. Flow velocities and grain sizes were varied in order to cover a large range of Peclet numbers including typical field scenarios. The laboratory study was extended by a comprehensive literature search to compare the new results with earlier work. As a result we propose a new empirical relationship for prediction of transverse dispersion coefficients (Dt) which is based on the Peclet number (Pe). This new relationship indicates a nonlinear dependency on the flow velocity (nu a) and grain size (d), namely a relative decrease of the dispersion coefficient with increasing flow velocity in relatively fast flowing water: Dt/Daq=Dp/Daq+0.28(Pe)0.72 (with Pe=nu a d/Daq; Daq and Dp denote the aqueous and pore diffusion coefficients, resp.).     

Matrix diffusion coefficients in volcanic rocks at the Nevada test site: influence of matrix porosity, matrix permeability, and fracture coating minerals

Diffusion cell experiments were conducted to measure nonsorbing solute matrix diffusion coefficients in forty-seven different volcanic rock matrix samples from eight different locations (with multiple depth intervals represented at several locations) at the Nevada Test Site. The solutes used in the experiments included bromide, iodide, pentafluorobenzoate (PFBA), and tritiated water ((3)HHO). The porosity and saturated permeability of most of the diffusion cell samples were measured to evaluate the correlation of these two variables with tracer matrix diffusion coefficients divided by the free-water diffusion coefficient (D(m)/D*). To investigate the influence of fracture coating minerals on matrix diffusion, ten of the diffusion cells represented paired samples from the same depth interval in which one sample contained a fracture surface with mineral coatings and the other sample consisted of only pure matrix. The log of (D(m)/D*) was found to be positively correlated with both the matrix porosity and the log of matrix permeability. A multiple linear regression analysis indicated that both parameters contributed significantly to the regression at the 95% confidence level. However, the log of the matrix diffusion coefficient was more highly-correlated with the log of matrix permeability than with matrix porosity, which suggests that matrix diffusion coefficients, like matrix permeabilities, have a greater dependence on the interconnectedness of matrix porosity than on the matrix porosity itself. The regression equation for the volcanic rocks was found to provide satisfactory predictions of log(D(m)/D*) for other types of rocks with similar ranges of matrix porosity and permeability as the volcanic rocks, but it did a poorer job predicting log(D(m)/D*) for rocks with lower porosities and/or permeabilities. The presence of mineral coatings on fracture walls did not appear to have a significant effect on matrix diffusion in the ten paired diffusion cell experiments.     

Influence of the mode of matrix porosity determination on matrix diffusion calculations

The theoretical basis for matrix diffusion in fractured rocks and the methodology for the determination of diffusion coefficients in the laboratory are well established. One significant problem, however, remains in that it is difficult to quantify the degree of sample disturbance affecting the geometrical, geophysical and hydraulic properties of the rock matrix. A new technique, with in situ rock impregnation with resin, for examining the diffusion-accessible rock matrix has been developed and successfully adopted to the rock matrix behind a water-conducting fracture in host crystalline rocks at Nagra's Grimsel Test Site in Switzerland and JNC's Kamaishi In Situ Test Site in Japan. In line with the results of a large number of natural analogue and laboratory studies, the existence of an in situ interconnected pore network was substantiated. Matrix porosities determined on the laboratory samples from both the sites are 1.5-3 times higher than in situ values, irrespective of the technique applied. On the Grimsel granodiorite matrix, matrix porosity existing in situ and artefacts of stress release and physical disturbance, induced by sampling and sample preparation, were clearly distinguished, allowing in situ porosity to be quantified. Laboratory work with conventional techniques tends to overestimate the porosity of the rock matrix, hence leading to an overestimation of in situ matrix diffusion. The implications of these differences to a repository performance assessment are assessed with a couple of examples from existing assessments, and recommendations for future approaches to the examination of in situ matrix porosity are made.     

A technique for estimating one-dimensional diffusion coefficients in low-permeability sedimentary rock using X-ray radiography: comparison with through-diffusion measurements

The measurement of diffusive properties of low-permeability rocks is of interest to the nuclear power industry, which is considering the option of deep geologic repositories for management of radioactive waste. We present a simple, non-destructive, constant source in-diffusion method for estimating one-dimensional pore diffusion coefficients (D(p)) in geologic materials based on X-ray radiography. Changes in X-ray absorption coefficient (Deltamicro) are used to quantify changes in relative concentration (C/C(0)) of an X-ray attenuating iodide tracer as the tracer solution diffuses through the rock pores. Estimated values of D(p) are then obtained by fitting an analytical solution to the measured concentration profiles over time. Measurements on samples before and after saturation with iodide can also be used to determine iodide-accessible porosity (phi(I)). To evaluate the radiography method, results were compared with traditional steady-state through-diffusion measurements on two rock types: shale and limestone. Values of D(p) of (4.8+/-2.5)x10(-11) m(2).s(-1) (mean+/-standard deviation) were measured for samples of Queenston Formation shale and (2.6+/-1.0)x10(-11) m(2).s(-1) for samples of Cobourg Formation limestone using the radiography method. The range of results for each rock type agree well with D(p) values of (4.6+/-2.0)x10(-11) m(2).s(-1) for shale and (3.5+/-1.8)x10(-11) m(2).s(-1) for limestone, calculated from through-diffusion experiments on adjacent rock samples. Low porosity (0.01 to 0.03) and heterogeneous distribution of porosity in the Cobourg Formation may be responsible for the slightly poorer agreement between radiography and through-diffusion results for limestones. Mean values of phi(I) for shales (0.060) and limestones (0.028) were close to mean porosity measurements made on bulk samples by the independent water loss technique (0.062 and 0.020 for shales and limestones, respectively). Radiography measurements offer the advantage of time-saving for diffusion experiments because the experiment does not require steady-state conditions and also allows for visualization of the small-scale heterogeneities in diffusive properties within rocks at the mm to cm scale.     

Analysis of the migration of instantaneously injected cesium in artificial fractures of Lac du Bonnet granite, Manitoba, Canada

Breakthrough curves of ¹³⁷Cs and tritiated water injected instantaneously into artificial fractures in Lac du Bonnet granite were analyzed using the analytical solution for a single rock-fracture system and assuming the linear sorption isotherm of the solute. Parameters of nuclide diffusion and sorption in rock matrices, obtained by fitting, varied depending on the flow velocity in the fractures. According to theoretical calculations, different fracture flow velocities lead to different diffusion distances of nuclides in matrices at the same injection volume. As microscopic inhomogeneity is considered to exist in the rock matrix, the average diffusion-sorption characteristics of the matrix within the diffusion distance may have varied with the fracture flow velocity. Surface sorption was marked in fractures that had relatively high matrix sorption-diffusion capacities. The phenomenon was interpreted using the theoretical relationships developed between the surface sorption, matrix sorption and pore diffusion coefficient, and the porosity of matrices.The effect of the nonlinear sorption of solute was examined by numerically solving model equations that incorporate the nonlinear isotherm. This incorporation may contribute to the reduction of deviations between theoretical and experimental BTC's.     

Three-dimensional diffusion of non-sorbing species in porous sandstone: computer simulation based on X-ray microtomography using synchrotron radiation

The diffusion pathways of porous sandstone were examined by a three-dimensional (3-D) imaging technique based on X-ray computed tomography (CT) using the SPring-8 (Super Photon ring-8 GeV, Hyogo, Japan) synchrotron radiation facility. The analysis was undertaken to develop better understanding of the diffusion pathways in natural rock as a key factor in clarifying the detailed mechanism of the diffusion of radionuclides and water molecules through the pore spaces of natural barriers in underground nuclear waste disposal facilities. A cylindrical sample (diameter 4 mm, length 6 mm) of sandstone (porosity 0.14) was imaged to obtain a 3-D image set of 450(3) voxels=2.62(3) mm(3). Through cluster-labeling analysis of the 3-D image set, it was revealed that 89% of the pore space forms a single large pore-cluster responsible for macroscopic diffusive transport, while only 11% of the pore space is made up of isolated pores that are not involved in long-range diffusive transport. Computer simulations of the 3-D diffusion of non-sorbing random walkers in the largest pore cluster were performed to calculate the surface-to-volume ratio of the pore, tortuosity (diffusion coefficient in free space divided by that in porous rock). The results showed that (i) the simulated surface-to-volume ratio is about 60% of the results obtained by conventional pulsed-field-gradient proton nuclear magnetic resonance (NMR) laboratory experiments and (ii) the simulated tortuosity is five to seven times larger than the results of laboratory diffusion experiments using non-sorbing I(-) and Br(-). These discrepancies are probably attributed to the intrinsic sample heterogeneity and limited spatial resolution of the CT system. The permeability was also estimated based on the NMR diffusometry theory using the results of the random walk simulations via the Kozeny-Carman equation. The estimated permeability involved an error of about 20% compared with the permeability measured by the conventional method, suggesting that the diffusometry-based NMR well logging with gradient coils is applicable to the in-situ permeability measurement of strata. The present study demonstrated that X-ray CT using synchrotron radiation is a powerful tool for obtaining 3-D pore structure images without the beam-hardening artifacts inevitable in conventional CT using X-ray tubes.     

Changes in the pore network structure of Hanford sediment after reaction with caustic tank wastes

At the former nuclear weapon production site in Hanford, WA, caustic radioactive tank waste leaks into subsurface sediments and causes dissolution of quartz and aluminosilicate minerals, and precipitation of sodalite and cancrinite. This work examines changes in pore structure due to these reactions in a previously-conducted column experiment. The column was sectioned and 2D images of the pore space were generated using backscattered electron microscopy and energy dispersive X-ray spectroscopy. A pre-precipitation scenario was created by digitally removing mineral matter identified as secondary precipitates. Porosity, determined by segmenting the images to distinguish pore space from mineral matter, was up to 0.11 less after reaction. Erosion-dilation analysis was used to compute pore and throat size distributions. Images with precipitation had more small and fewer large pores. Precipitation decreased throat sizes and the abundance of large throats. These findings agree with previous findings based on 3D X-ray CMT imaging, observing decreased porosity, clogging of small throats, and little change in large throats. However, 2D imaging found an increase in small pores, mainly in intragranular regions or below the resolution of the 3D images. Also, an increase in large pores observed via 3D imaging was not observed in the 2D analysis. Changes in flow conducting throats that are the key permeability-controlling features were observed in both methods.     

Diffusive flux and pore anisotropy in sedimentary rocks

Diffusion of dissolved contaminants into or from bedrock matrices can have a substantial impact on both the extent and longevity of dissolved contaminant plumes. For layered rocks, bedding orientation can have a significant impact on diffusion. A series of laboratory experiments was performed on minimally disturbed bedrock cores to measure the diffusive flux both parallel and normal to mineral bedding of four different anisotropic sedimentary rocks. Measured effective diffusion coefficients ranged from 4.9×10(-8) to 6.5×10(-7)cm(2)/s. Effective diffusion coefficients differed by as great as 10-folds when comparing diffusion normal versus parallel to bedding. Differences in the effective diffusion coefficients corresponded to differences in the "apparent" porosity in the orientation of diffusion (determined by determining the fraction of pore cross-sectional area measured using scanning electron microscopy), with the difference in apparent porosity between normal and parallel bedding orientations differing by greater than 2-folds for two of the rocks studied. Existing empirical models failed to provide accurate predictions of the effective diffusion coefficient in either bedding orientation for all four rock types studied, indicating that substantial uncertainty exists when attempting to predict diffusive flux through sedimentary rocks containing mineral bedding. A modified model based on the apparent porosity of the rocks provided a reasonable prediction of the experimental diffusion data.     

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.     

Modeling the Reduction of Vapor Phase Emissions from Surface Soils Due to Soil Matrix Effects: Porosity/Tortuosity Concepts

ABSTRACT

A major route for transport of volatile organic compounds within porous media is vapor phase diffusion. The diffusion rate through a porous medium is less than that through free-air due to the decreased cross-sectional area available for gas movement and the increased path length due to pore tortuosity. Numerous empirical expressions have been published that relate the diffusion coefficient in porous media to the diffusion coefficient in free-air (unobstructed gas phase). Published measurements of relative diffusivity and air-filled porosity were combined into a database. Empirical expressions available in the literature, including the popular Millington-Quirk equation, were evaluated along with a fourth-degree polynomial expression developed by the authors to determine the best type of equation to predict relative diffusivity as a function of air-filled porosity over the domain of values for porosity ranging from 0.071 to 1 for different types of materials. Mean square deviations were used as the statistical test to compare equations. The polynomial expression developed in this project produced a significantly different effective diffusion coefficient (1.3 x 10^-6 m²/sec) compared to values of 9.2 x 10^-6 m²/sec and 3.1 x 10^-6 m²/ sec predicted by forms of the Millington-Quirk equation for a specific case.     

Field investigation into unsaturated flow and transport in a fault: model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture-matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault-matrix and fracture-matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture-matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model.     

A multi-directional tracer test in the fractured Chalk aquifer of E. Yorkshire, UK

A multi-borehole radial tracer test has been conducted in the confined Chalk aquifer of E. Yorkshire, UK. Three different tracer dyes were injected into three injection boreholes and a central borehole, 25 m from the injection boreholes, was pumped at 330 m(3)/d for 8 days. The breakthrough curves show that initial breakthrough and peak times were fairly similar for all dyes but that recoveries varied markedly from 9 to 57%. The breakthrough curves show a steep rise to a peak and long tail, typical of dual porosity aquifers. The breakthrough curves were simulated using a 1D dual porosity model. Model input parameters were constrained to acceptable ranges determined from estimations of matrix porosity and diffusion coefficient, fracture spacing, initial breakthrough times and bulk transmissivity of the aquifer. The model gave equivalent hydraulic apertures for fractures in the range 363-384 microm, dispersivities of 1 to 5 m and matrix block sizes of 6 to 9 cm. Modelling suggests that matrix block size is the primary controlling parameter for solute transport in the aquifer, particularly for recovery. The observed breakthrough curves suggest results from single injection-borehole tracer tests in the Chalk may give initial breakthrough and peak times reasonably representative of the aquifer but that recovery is highly variable and sensitive to injection and abstraction borehole location. Consideration of aquifer heterogeneity suggests that high recoveries may be indicative of a high flow pathway adjacent, but not necessarily connected, to the injection and abstraction boreholes whereas low recoveries may indicate more distributed flow through many fractures of similar aperture.     

Physicochemical controls on initiation and evolution of desiccation cracks in sand-bentonite mixtures: X-ray CT imaging and stochastic modeling

The shrink-swell behavior of active clays in response to changes in physicochemical conditions creates great challenges for construction of geotechnical barriers for hazardous waste isolation, and is of significant importance for management of agricultural and natural resources. Initiation and evolution of desiccation cracks in active clays are strongly dependent on physicochemical initial and boundary conditions. To investigate effects of bentonite content (20, 40, 60%), pore fluid chemistry (0.05 and 0.5M NaCl) and drying rates (40 and 60°C) on cracking behavior, well-controlled dehydration experiments were conducted and X-ray Computed Tomography (CT) was applied to visualize and quantify geometrical features of evolving crack networks. A stochastic model based on the Fokker-Plank equation was adopted to describe the evolution of crack aperture distributions (CAD) and to assess the impact of physicochemical factors on cracking behavior. Analyses of crack porosity and crack specific surface area showed that both clay content and temperature had larger impact on cracking than pore fluid concentration. More cracks formed at high bentonite contents (40 and 60%) and at high drying rate (60°C). The drift, diffusion and source terms derived from stochastic analysis indicated that evaporative demand had greater influence on the dynamics of the CAD than solution chemistry.