In situ anion diffusion experiments using radiotracers

Diffusion experiments in compacted bentonite have been carried out in situ using the borehole laboratory CHEMLAB. The "ordinary" anion iodide and the redox-sensitive pertechnetate ion have been investigated. In spite of strongly reducing groundwater conditions, technetium was found to diffuse mostly unreduced as TcO4-, although in some spots in the compacted clay, the activity was significantly higher, which may be explained by reduction of some TcO4- by iron-containing minerals in the bentonite. The measured concentration profiles in the clay cannot be accommodated by assuming one single diffusion process. The experimental data are modeled assuming two diffusion paths, intralamellar diffusion and diffusion in external water. The apparent diffusivity for the intralamellar diffusion was found to be 8.6 x 10(-11) m2 s(-1) for iodide with a capacity factor of 0.1, while the apparent diffusivity for the diffusion in external water was found to be 5 x 10(-14) m2 s(-1) with alpha=2.26. The corresponding values for Tc were found to be Da= 6 x 10(-11) m2 s(-1), alpha=0.1 and Da= 1 x 10(-13) m2 s(-1), alpha=0.46, respectively. The diffusion constants and capacity factors obtained in this study are in accordance with data from laboratory experiments.     

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.     

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.     

Formation factor logging by electrical methods. Comparison of formation factor logs obtained in situ and in the laboratory

In this paper, a new in situ method for obtaining the formation factor, which is essential for the matrix diffusion, is described and tested in intrusive igneous rock. The method is based on electrical resistivity measurements in rock where the pore water and rock resistivities are essential parameters. The method is based on electromigration instead of diffusion as in traditional diffusion experiments. In previous works, quantitative formation factors of rock have been obtained by electrical methods in the laboratory. Here, a similar approach is used in situ. An in situ logging campaign was performed by SKB during 2000 in the 1700-m-deep borehole KLX02 in Laxemar, Sweden. The rock resistivity was measured with the slimhole Dual Laterolog from Antares. The groundwater resistivity was measured with the Difference Flow Meter from Posiva. A formation factor log was obtained with the maximum vertical resolution of 10 cm. In order to validate the log, 100 rock samples were taken from the bore core, and a formation factor log was obtained by using electrical methods in the laboratory. Both direct current (DC) and alternating current (AC) were used. The measurements on the core confirmed that the in situ log was quantitative, but with a possible systematic error. The in situ formation factors were on average about 1/3 to 1/5 of the laboratory formation factors, depending on depth.     

Modelling of silica diffusion experiments with 32Si in Boom Clay

A mathematical model describing the dissolution of nuclear glass directly disposed in clay combines a first-order dissolution rate law with the diffusion of dissolved silica in clay. According to this model, the main parameters describing the long-term dissolution of the glass are etaR, the product of the diffusion accessible porosity eta and the retardation factor R, and the apparent diffusion coefficient D(app) of dissolved silica in clay. For determining the migration parameters needed for long-term predictions, four Through-Diffusion (T-D) experiments and one percolation test have been performed on undisturbed clay cores. In the Through-Diffusion experiments, the concentration decrease after injection of 32Si (radioactive labelled silica) was measured in the inlet compartment. At the end of the T-D experiments, the clay cores were cut in thin slices and the activity of labelled silica in each slice was determined. The measured activity profiles for these four clay cores are well reproducible. Since no labelled silica could be detected in the outlet compartments, the Through-Diffusion experiments are fitted by two In-Diffusion models: one model assuming linear and reversible sorption equilibrium and a second model taking into account sorption kinetics. Although the kinetic model provides better fits, due to the sufficiently long duration of the experiments, both models give approximately similar values for the fit parameters. The single percolation test leads to an apparent diffusion coefficient value about two to three times lower than those of the Through-Diffusion tests. Therefore, dissolved silica appears to be strongly retarded in Boom Clay. A retardation factor R between 100 and 300 was determined. The corresponding in situ distribution coefficient K(d) is in the range 25-75 cm(3) g(-1). The apparent diffusion coefficient of dissolved silica in Boom Clay is estimated between 2 x 10(-13) and 7 x 10(-13) m(2) s(-1). The pore diffusion coefficient is in the range from 6 x 10(-11) to 1 x 10(-10) m(2) s(-1).     

A laboratory technique for investigation of diffusion and transformation of volatile organic compounds in low permeability media

A laboratory diffusion cell technique that permits spatial and temporal estimates of porewater concentrations over short intervals suitable for estimation of effective diffusion coefficients (De) and degradation rate constants (k) of volatile organic compounds (VOCs) in saturated low permeability media is presented. The diffusion cell is a sealed cylinder containing vapour reservoirs for sampling, including a vapour reservoir source and an array of vapour-filled "mini-boreholes", which are maintained water- and sediment-free by slightly negative porewater pressures. The vapour reservoirs were sampled by Solid Phase Micro-Extraction (SPME), resulting in minimal disturbance to the experimental system. Porewater concentrations are estimated from the measured vapour concentrations. Experiments were conducted using a non-reactive medium and five VOCs with a range in partitioning properties. Calibration experiments showed that equilibrium partition coefficients could be used for calculating concentrations in the vapour reservoir source from concentrations in the SPME coating after a 1-min microextraction and that the reservoir concentration was insignificantly affected by sampling. However, equilibrium was not reached during the one-min extraction of the boreholes; the microextraction reduced the borehole vapour concentrations, leading to diffusion of VOCs from porewater into the vapour-filled borehole. Thus, empirical partitioning coefficients were required for the determination of porewater VOC concentrations. The experimental data and numerical modelling indicate masses extracted by SPME extraction are relatively small, with minimal perturbation on processes studied in diffusion experiments. This technique shows promise for laboratory investigation of diffusion and transformation processes in low permeability media.     

Diffusion of HTO, 36Cl- and 125I- in Opalinus Clay samples from Mont Terri. Effect of confining pressure

Diffusion coefficients (T=23 +/- 2 degrees C) and accessible porosities for HTO, 36Cl(-) and 125I(-) were measured on Opalinus Clay (OPA) samples from the Mont Terri Underground Rock Laboratory (URL) using the through-diffusion technique. The direction of transport (diffusion) was perpendicular to bedding. Special cells that allowed the application of confining pressure were designed and constructed. The pressures ranged from 1 to 5 MPa, the latter value simulating the overburden at the Mont Terri URL (about 200 m). The test solution used in the experiments was a synthetic version of the Opalinus Clay pore water, which has Na(+) and Cl(-) as the main components (I=0.42 M). The measured values of the effective diffusion coefficients (D(e)) and rock capacity factors (alpha) are: D(e)=1.2-1.5 x 10(-11) m(2) s(-1) and alpha=0.09-0.11 for HTO, D(e)=4.0-5.5 x 10(-12) m(2) s(-1) and alpha=0.05 for 36Cl(-) and D(e)=3.2-4.6 x 10(-12) m(2) s(-1) and alpha=0.07-0.10 for 125I(-). For non-sorbing tracers (HTO, 36Cl) the rock capacity factor alpha is equal to the diffusion-accessible porosity epsilon. The experimental results showed that pressure only had a small effect on the value of the diffusion coefficients. Increasing the pressure from 1 to 5 MPa resulted in a decrease of the diffusion coefficient of approximately 17% for HTO, approximately 28% for 36Cl(-) and approximately 30% for 125I(-). Moreover, the diffusion coefficients for 36Cl(-) and 125I(-) are smaller than for HTO, which is consistent with an effect arising from anion exclusion. The diffusion coefficients of HTO and 125I(-) measured in this study are in good agreement with recent measurements at three other laboratories performed within the framework of a laboratory comparison exercise. The values of the diffusion-accessible porosities show a larger degree of scatter.     

Characterisation of HTO diffusion properties by an in situ tracer experiment in Opalinus clay at Mont Terri

A long-term single borehole diffusion experiment using tritiated water as tracer was carried out in Opalinus clay, an argillaceous rock formation that is accessible at the Mont Terri Underground Research Laboratory, situated in the Swiss Jura. The tracer was diluted in reconstituted formation water and introduced into a packed-off section of a borehole located in saturated rock. Pressure in this interval was maintained equal to the pore pressure of the surrounding rock in order to prevent any hydraulic gradient around the borehole and to avoid advective transport processes. The evolution of the tracer concentration in the injection system was monitored over time. After 1 year of diffusion, the claystone surrounding the interval was retrieved by overcoring the whole borehole and packer system, and by an adjacent oblique borehole. Compressed air was used as drilling fluid to reduce rock disturbances. The recovered overcore was sampled along profiles perpendicular to the borehole wall with a view to determining the tracer-concentration profiles in the rock. To avoid further evaporation of tritiated water, subsamples were immediately transferred into polyethylene bottles and disaggregated by adding a known amount of tracer-free water. Fifteen profiles were determined and showed a decreasing tracer concentration with distance into the rock. The pore-water contents were constant along those profiles, confirming that only very little water was lost during overcoring operations. The evolution of tritium-tracer concentration in the injection system over time and in situ profiles were interpreted with a 3-D numerical simulation of the experiment. That allowed for the identification of the transport parameters (orthotropic diffusion tensor and porosity) by minimising the relative quadratic error between the experimental and simulated data. The fitting is good and the results are consistent with data obtained on drill-core samples. The result of tritiated water is discussed regarding (1) the potential effect of mechanical and/or chemical disturbances around the injection borehole and (2) the specific behaviour of tritiated water.     

Prediction of some in situ tracer tests with sorbing tracers using independent data

Some recent converging tracer tests with sorbing tracers at the Asp? Hard Rock Laboratory in Sweden, the TRUE tests, have been predicted using only laboratory data and hydraulic data from borehole measurements. No model parameters were adjusted to obtain a better fit with the experiments. The independent data were fracture frequency and transmissivity data obtained in the field and laboratory data on sorption and matrix diffusion. Transmissivity measurements in five boreholes in the rock volume containing the region surrounding the injection and collection points show that there is a high frequency of water conducting fractures. Of 162 packed off sections with 0.5 m packer distances, 112 were found to have a transmissivity above the detection limit. The specific flow-wetted surface (FWS) of the rock mass could be estimated from these data. The transmissivities were found to be reasonably well described by a lognormal distribution. Laboratory data on diffusion and sorption properties together with the hydraulic data were used to "predict" the residence time distribution (RTD) of the sorbing tracers. The results were compared with the experimental breakthrough curves. In these experiments, the water residence time is very small compared to the residence time of the sorbing tracers due to their diffusion and sorption within the rock matrix. We thus could neglect the influence of the water residence time in our predictions. Therefore, no information on water residence times or on "dispersion" was needed. The dispersion of the sorbing tracers is caused by the different sorbing tracer residence times in different pathways. The sorbing tracer residence time is determined by the ratio of flowrate to the flow-wetted surface in the different pathways and not by the water residence time. Assuming a three-dimensional flow pattern and using the observed fracture frequency and flowrate distribution, breakthrough curves for three strongly sorbing tracers were predicted. Only the laboratory data, the transmissivity measurements and the pumping flowrate were used in the predictions. No information on the water residence time as obtained by the nonsorbing tracers was used. The predictions were surprisingly accurate.     

Estimation of mechanical dispersion and dispersivity in a soil-gas system by column experiments and the dusty gas model

In a previous study, column experiments were carried out with Toyoura sand (permeability 2.05×10(-11)m(2)) and Toyoura sand mixed with bentonite (permeability 9.96×10(-13)m(2)) to obtain the molecular diffusion coefficient, the Knudsen diffusion coefficient, the tortuosity for the molecular diffusion coefficient, and the mechanical dispersion coefficient of soil-gas systems. In this study, we conducted column experiments with field soil (permeability 2.0×10(-13)m(2)) and showed that the above parameters can be obtained for both less-permeable and more-permeable soils by using the proposed method for obtaining the parameters and performing column experiments. We then estimated dispersivity from the mechanical dispersion coefficients obtained by the column experiments. We found that the dispersivity depended on the mole fraction of the tracer gas and could be represented by a quadratic equation.     

Rock matrix diffusivity determinations by in-situ electrical conductivity measurements

A fast method to determine rock matrix diffusion properties directly in the bedrock would be valuable in the investigation of a possible site for disposal of nuclear waste. An "effective diffusivity borehole log" would provide important information on the variability of this entity over the area studied. As opposed to traditional matrix diffusion laboratory experiments, electrical conductivity measurements are fast, inexpensive and also easy to carry out in-situ. In this study, electrical resistivity data from borehole logging, as well as from measurements on the actual core, is evaluated with the purpose of extracting matrix diffusivity data. The influence of migration of ions in the electrical double layer, which can be of great importance in low ionic strength pore water, is also considered in evaluating the in-situ data to accurately determine the effective pore diffusivity. The in-situ data compare fairly well to those measured in the rock core.     

Apparent diffusion coefficients and chemical species of neptunium (V) in compacted Na-montmorillonite

Diffusion of neptunium (V) in compacted Na-montmorillonite was studied through the non-steady state diffusion method. In this study, two experimental attempts were carried out to understand the diffusion mechanism of neptunium. One was to establish the diffusion activation energy, which was then used to determine the diffusion process in the montmorillonite. The other was the measurement of the distribution of neptunium in the montmorillonite by a sequential batch extraction. The apparent diffusion coefficients of neptunium in the montmorillonite at a dry density of 1.0 Mg m-3 were from 3.7 x 10(-12) m2 s-1 at 288 K to 9.2 x 10(-12) m2 s-1 at 323 K. At a dry density of 1.6 Mg m-3, the apparent diffusion coefficients ranged between 1.5 x 10(-13) m2 s-1 at 288 K and 8.7 x 10(-13) m2 s-1 at 323 K. The activation energy for the diffusion of neptunium at a dry density of 1.0 Mg m-3 was 17.5 +/- 1.9 kJ mol-1. This value is similar to those reported for diffusion of other ions in free water, e.g., 18.4 and 17.4 kJ mol-1 for Na+ and Cl-, respectively. At a dry density of 1.6 Mg.m-3, the activation energy was 39.8 +/- 1.9 kJ mol-1. The change in the activation energy suggests that the diffusion process changes depending on the dry density of the compacted montmorillonite. A characteristic distribution profile was obtained by the sequential extraction procedure for neptunium diffused in compacted montmorillonite. The estimated fraction of neptunium in the pore water was between 3% and 11% at a dry density of 1.6 Mg m-3 and at a temperature of 313 K. The major fraction of the neptunium in the montmorillonite was identified as neptunyl ions sorbed on the outer surface of the montmorillonite. These findings suggested that the activation energy for diffusion and the distribution profile of the involved nuclides could become powerful parameters in understanding the diffusion mechanism.     

Direct simulation of heterogeneous diffusion and inversion procedure applied to an out-diffusion experiment. Test case of Palmottu granite

An out-diffusion laboratory experiment using a non-reactive tracer was fitted using the Time Domain Diffusion (TDD) method. This rapid particle tracking method allows simulation of the heterogeneous diffusion based on pore-scale images and local values of diffusivities. The superimposed porosity and mineral 2D maps act as computation grids to condition diffusion pathways. We focused on a Palmottu granite sample, in which the connected pore space has a composite microstructure with cracks linking microporous minerals and is above the percolation threshold. Three main results were achieved: (i) When compared to the fitting obtained with one coefficient (best mean square residual R = 1.6 x 10(-2)), diffusion is shown to be suitably characterised with two coefficients related to cracks and microporous minerals (best R = 6.5 x 10(-4)), (ii) rather than imposing a local apparent diffusion coefficient D(a) independent of the local porosity Phi, a best fit is obtained by applying Archie's relationship D(a) = D(0) x G with G = Phi(m) to each pixel of the calculation grids (G is the geometry factor, D(0) is the diffusion coefficient in free fluid, and m is Archie's exponent), and (iii) the order of magnitude of the fitted diffusion coefficient or Archie's exponents (m=0 for microcracks and m=1.82 for microporous minerals) is physically realistic.     

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.     

In situ tracer tests to determine retention properties of a block scale fracture network in granitic rock at the Aspö Hard Rock Laboratory, Sweden

Experiments were conducted at the Asp? Hard Rock Laboratory in order to improve the understanding of radionuclide retention properties of fractured crystalline bedrock in the 10-100 m scale (TRUE Block Scale Project, jointly funded by ANDRA, ENRESA, Nirex, JNC, Posiva and SKB). A series of tracer experiments were performed using sorbing tracers in three different flow paths. The different flow paths had Euclidian lengths of 14, 17 and 33 m, respectively, and one to three water conducting structures. Four tests were performed using different cocktails made up of radioactive sorbing tracers (22,24Na+, 42K+, 47Ca2+, 85Sr2+, 83,86Rb+, 131,133Ba2+ and 134,137Cs+). For each tracer injection, the breakthrough of sorbing tracers was compared to the breakthrough of a conservative tracer, 82Br-, 131I-, HTO and 186ReO4-, respectively. In the two longer flow paths, no breakthrough of 83Rb+ and 137Cs+ was observed after 8 months of pumping. Selected tracer tests were subject to basic modelling in which a one-dimensional (1D) advection-dispersion model, including surface sorption, and an unlimited matrix diffusion were used for the interpretation of the results. The results of the modelling indicated that there is a slightly higher mass transfer into a highly porous material in the block-scale experiment compared with in situ experiments performed over shorter distances and significantly higher than what would have been expected from laboratory data obtained from studies of the interactions in nonaltered intact rock.     

Solute transport in crystalline rocks at Aspö--II: blind predictions,inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish Asp? Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na(+), 85Sr(2+), 47Ca(2+)and more strongly sorbing 86Rb(+), 133Ba(2+), 137Cs(+). Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2-15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted K(d) values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25-60, but agree within a factor of 3-5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning K(d)s is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.     

The effects of matrix diffusion on radionuclide migration in rock column experiments

Rock column experiments were performed to examine the effects of matrix diffusion and hydrodynamic dispersion on the migration of radionuclides at the laboratory scale. Tritiated water and chloride transportation was studied in intact mica gneiss and in altered more porous tonalite columns with narrow flow channels. The column diffusion properties were estimated prior to water flow experiments using the gas diffusion method with helium as the tracer gas. The numerical compartment model for advection and dispersion, with and without matrix diffusion, was used to interpret the tracer transport in the columns. Matrix diffusion behavior was also distinguished from dominating hydrodynamic dispersion in rock column experiments at the slowest water flow rates.     

Tracer test at El Berrocal site

Tracer tests provide highly valuable information about the transport properties of saturated rocks which is essential to the characterization of a potential radioactive waste repository site. In the frame of El Berrocal project, a set of tracer tests was performed in a complex geometry of inclined boreholes, combined with highly fractured transmissive zones. The aims of the tracer test programme were to gain experience, knowledge and insight into field transport experiments. To achieve this a detailed programme of new instrumentation design, site characterization and laboratory tasks was developed. For field monitoring a new electronic system was developed. The system is able to measure up to 256 parameters per borehole, with surface equipment to control pumping rates and physical and chemical parameters at both injection and extraction boreholes. The experiments progressed from single borehole dilution tests under both natural flow and forced gradient conditions to convergent flow tracer tests. Dilution tests helped to discriminate the most suitable borehole sections at which to inject tracers. The tracers were selected by the results of the laboratory programme. Uranine (fluorescein), DTPA-gadolinium (diethylenetriaminopentacetic acid, gadolinium (III)), and deuterium were injected simultaneously in one borehole section and recovered at another borehole 20 m away, pumping at a flowrate of 0.1 1 min⁻¹. First results showed a thickness porosity of 1.2 × 10⁻³ m and a longitudinal dispersivity of 17.0 m using uranine data acquired over a period of 4 d, at which point the recovery concentration had reached a maximum. However, gadolinium and deuterium appeared to travel faster, arriving at peak values after only 2 d of injection.     

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.     

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.