Saline groundwaters counteract up-flow of contaminants- implications for radionuclide repositories?

The high-level nuclear waste, HLW, from Swedish and Finnish reactors will be deposited in crystalline rock at depths around 500 m. The waste is enclosed in steel canisters protected against corrosion by a 5 cm thick copper shell, which ensures a lifetime far longer than 100 000 years. Should some canister be breached any leaking nuclides will have decayed to so low activity that even if they reached the biosphere, they would cause minimal risk to humans. The cost of the copper is significant. The dismantling of the nuclear reactors, with induced activity must also be disposed of and this waste volume is much larger than that of the HLW, which makes it impossible to protect it in the same way. This paper explores if by locating the waste at larger depth where the ground water is more saline, and where the hydraulic conductivity of the rock is lower up-flow of contaminated water can be ensured to be negligible because the denser water at larger depth counteracts up-flow due to negative buoyancy. Several processes that could cause local up-flow are addressed, such as infiltration of meteoric water, impact of surface topology, heat production of the waste, geothermal gradient, salinity gradient, hydraulic conductivity heterogeneities and salt migration between seeping water and salt in matrix pore water. Flow and transport simulations using data from extensive field investigations over more than ten years with scores of km deep boreholes suggest that a HLW repository at around one km depth may be sufficient to hinder up-flow to the biosphere.

Atomistic simulations of cation hydration in sodium and calcium montmorillonite nanopores.

During the last four decades, numerous studies have been directed to the swelling smectite-rich clays in the context of high-level radioactive waste applications and waste-liners for contaminated sites. The swelling properties of clay mineral particles arise due to hydration of the interlayer cations and the diffuse double layers formed near the negatively charged montmorillonite (MMT) surfaces. To accurately study the cation hydration in the interlayer nanopores of MMT, solvent-solute and solvent-clay surface interactions (i.e., the solvation effects and the shape effects) on the atomic level should be taken into account, in contrast to many recent electric double layer based methodologies using continuum models. Therefore, in this research we employed fully atomistic simulations using classical molecular dynamics (MD) simulations, the software package GROMACS along with the CLAYFF forcefield and the SPC/E water model. We present the ion distributions and the deformation of the hydrated coordination structures, i.e., the hydration shells of Na+ and Ca2+ in the interlayer, respectively, for MMT in the first-layer, the second-layer, the third-layer, the fourth-layer, and the fifth-layer (1W, 2W, 3W, 4W, and 5W) hydrate states. Our MD simulations show that Na+ in Na-MMT nanopores have an affinity to the ditrigonal cavities of the clay layers and form transient inner-sphere complexes at about 3.8 Å from clay midplane at water contents less than the 5W hydration state. However, these phenomena are not observed in Ca-MMT regardless of swelling states. For Na-MMT, each Na+ is coordinated to four water molecules and one oxygen atom of the clay basal-plane in the first hydration shell at the 1W hydration state, and with five to six water molecules in the first hydration shell within a radius of 3.1 Å at all higher water contents. In Ca-MMT, however each Ca2+ is coordinated to approximately seven water molecules in the first hydration shell at the 1W hydration state and about eight water molecules in the first hydration shell within a radius of 3.3 Å at all higher hydration states. Moreover, the MD results show that the complete hydration shells are nearly spherical with an orthogonal coordination sphere. They could only be formed when the basal spacing d001 ≥ 18.7 Å, i.e., approximately, the interlayer separation h ≥ 10 Å. Comparison between DFT and MD simulations shows that DFT failed to reproduce the outer-sphere complexes in the Stern-layer (within ∼5.0 Å from the clay basal-plane), observed in the MD simulations.

Solute transport in a single fracture involving an arbitrary length decay chain with rock matrix comprising different geological layers.

A model is developed to describe solute transport and retention in fractured rocks. It accounts for advection along the fracture, molecular diffusion from the fracture to the rock matrix composed of several geological layers, adsorption on the fracture surface, adsorption in the rock matrix layers and radioactive decay-chains. The analytical solution, obtained for the Laplace-transformed concentration at the outlet of the flowing channel, can conveniently be transformed back to the time domain by the use of the de Hoog algorithm. This allows one to readily include it into a fracture network model or a channel network model to predict nuclide transport through channels in heterogeneous fractured media consisting of an arbitrary number of rock units with piecewise constant properties. More importantly, the simulations made in this study recommend that it is necessary to account for decay-chains and also rock matrix comprising at least two different geological layers, if justified, in safety and performance assessment of the repositories for spent nuclear fuel.

The weighted correlation approach for density functional theory: a study on the structure of the electric double layer.

Within the framework of density functional theory, a weighted correlation approach is developed in order to obtain the density distributions of an inhomogeneous fluid. It results in a formally exact expression, by means of the concept of a weighted pair correlation function, used to evaluate the change of the single-particle direct correlation function of the system relative to that of a reference state. When applying the approach for practical use, however, an approximation of the pair correlation function has to be made, along with an appropriate definition of a weight function. Noticeably, combining this approach with fundamental measure theory gives rise to a new method, which we call the FMT/WCA-k(2) approach, for studying the structural and thermodynamic properties of a charged hard-sphere fluid subjected to a spatially varying external potential. Application of the FMT/WCA-k(2) approach in a range of electrolyte concentrations and surface charge densities, against the Monte Carlo simulations, shows that it is superior to the typical approaches of density functional theory in predicting the ionic density profiles of both counter-ions and co-ions near a highly charged surface. It is capable of capturing the fine features of the structural properties of the electric double layers, to well reproduce the layering effect and the charge inversion phenomenon, also in strongly coupled cases where divalent counter-ions are involved. In addition, it is found that the FMT/WCA-k(2) approach even has an advantage over the anisotropic, hyper-netted chain approaches in giving better agreement with the Monte Carlo results.

A novel method to describe the interaction pressure between charged plates with application of the weighted correlation approach.

Based on the Euler-Lagrange equation for ion density distribution in an inhomogeneous, charged, and hard-sphere fluid, a novel method is proposed to determine the interaction pressure between charged plates. The resulting expression is a sum of distinct physical contributions to the pressure, which involves different contributions to the single-particle direct correlation function. It can, therefore, be conveniently used in any density functional approach to facilitate analysis of the pressure components. In this study, the so-called fundamental measure theory (FMT)∕weighted correlation approach (WCA) approach is applied to estimate both the hard-sphere and the electric residual contributions to the single-particle direct correlation function, upon the calculation of the ionic density profiles between charged plates. The results, against the Monte Carlo simulations, show that the FMT∕WCA approach is superior to the typical FMT∕mean spherical approximation approach of the density functional theory in predicting the interaction pressure between charged plates immersed in an electrolyte solution upon various conditions in the primitive model. The FMT∕WCA approach can well capture the fine features of the pressure-separation dependence, to reproduce not only the shoulder shape and the weak attractions in monovalent electrolytes but also the strongly oscillatory behavior of pressure in divalent electrolytes where pronounced attractions are observed. In addition, it is found that the FMT∕WCA approach even has an advantage over the anisotropic, hyper-netted chain approach in that it agrees with the Monte Carlo results to a very good extent with, however, much less computational effort.

A novel approach to determine the critical coagulation concentration of a colloidal dispersion with plate-like particles.

The critical coagulation concentration (ccc) of counterions is commonly described by the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory on the basis of a static force balance. It can, however, also be estimated from a kinetic point of view by studying the process of colloidal coagulation, or from a dynamic point of view by considering colloidal transport in nonequilibrium systems where other processes such as diffusion and the influence of gravity come into play. In particular, in a test tube where colloidal expansion takes place, the ccc can be interpreted as the electrolyte concentration below which expansion of colloids would always lead to full access to the entire volume of the test tube and above which a sharp boundary is established between a colloidal gel and pure water. On the basis of this perception and the dynamic force balance model that we developed to describe colloidal expansion in a test tube, accounting for the effects of particle diffusion and gravity in contrast to the DLVO theory, we propose an alternative way to assess the ccc of counterions. We also derive an approximate expression for the case of homointeraction at constant charge for montmorillonite. The estimated ccc values agree quite well with those observed experimentally for both Na(+) and Ca(2+) counterions for montmorillonite dispersions, at pH approximately 6.5. This is in contrast to the DLVO theory, which overpredicts the ccc by about 2 orders of magnitude. In addition, the detailed analyses suggest that the ccc of counterions decreases with increasing surface area and with the thickness of the particles. For montmorillonite, the ccc is nearly independent of the surface charge density of the particles for the range of typical charge densities.

A dynamic force balance model for colloidal expansion and its DLVO-based application.

A force balance model that describes the dynamic expansion of colloidal bentonite gels/sols is presented. The colloidal particles are assumed to consist of one or several thin sheets with the other dimensions much larger than their thickness. The forces considered include van der Waals force, diffuse double layer force, thermal force giving rise to Brownian motion, gravity, as well as friction force. The model results in an expression resembling the instationary diffusion equation but with an immensely variable diffusivity. This diffusivity is strongly influenced by the concentration of counterions as well as by the particle concentration in the colloid gel/sol. The properties of the model are explored and discussed, exemplified by the upward expansion of an originally highly compacted bentonite tablet in a test tube. Examples are presented for a number of cases with ionic concentrations varying between very dilute waters up to several molar of counterions. The volume fraction of particles ranges from 40% to very dilute sols.

Fast method for simulation of radionuclide chain migration in dual porosity fracture rocks.

In fractured rocks with a porous rock matrix such as granites, radionuclides will flow with the water in the fracture network. The nuclides will diffuse in and out the rock matrix where they can sorb and be considerably retarded compared to the water velocity. A water parcel entering the network will mix and split at the fracture intersections and parts of the original parcel will traverse a multitude of different fractures. The flowrates, velocities, sizes and apertures of the fractures can vary widely. Normally one must solve the transport equations for every fracture and use the effluent concentration as inlet condition to the next fracture and so on. It is shown that under some weakly simplified conditions it suffices to determine one single parameter group containing information on the flow wetted surface that a water parcel contacts along the entire path. It is also shown how this can be obtained. Then, solving the transport equations only once for time and location along the path gives the concentration and nuclide flux of every nuclide in the chain everywhere along a path. The same solution actually is valid for every path in the network. This dramatically reduces the computation effort. The same approach can be used for models based on streamtubes.

Through-electromigration: a new method of investigating pore connectivity and obtaining formation factors.

The retardation of radionuclides and other contaminants in fractured crystalline rock is strongly associated with the diffusive properties of the rock matrix. At present, the scientific community is divided concerning the question of long-range pore connectivity in intrusive igneous rock. This paper presents a fast new method, called the through-electromigration method, of obtaining formation factors and investigating pore connectivity. The method involves the migration of an ionic tracer through a rock sample with an electrical potential gradient as the main driving force. The method is analogous to the through-diffusion method but the experimental time is reduced by orders of magnitude. This enables investigations of pore connectivity, as measurements can be made on longer samples. In a preliminary investigation, the new method is compared to the traditional through-diffusion method as well as to rock resistivity methods. The diffusive properties of nine granitic rock samples from Laxemar in Sweden, ranging from 15 to 121 mm in length, have been investigated and the results are compared.

Modelling radionuclide transport for time varying flow in a channel network.

Water flowrates and flow directions may change over time in the subsurface for a number of reasons. In fractured rocks flow takes place in channels within fractures. Solutes are carried by the advective flow. In addition, solutes may diffuse in and out of stagnant waters in the rock matrix and other stagnant water regions. Sorbing species may sorb on fracture surfaces and on the micropore surfaces in the rock matrix. We present a method by which solute particles can be traced in flowing water undergoing changes in flowrate and direction in a complex channel network where the solutes can also interact with the rock by diffusion in the rock matrix. The novelty of this paper is handling of diffusion in the rock matrix under transient flow conditions. The diffusive processes are stochastic and it is not possible to follow a particle deterministically. The method therefore utilises the properties of a probability distribution function for a tracer moving in a fracture where matrix diffusion is active. The method is incorporated in a three dimensional channel network model. Particle tracking is used to trace out a multitude of flowpaths, each of which consists of a large number of channels within fractures. Along each channel the aperture and velocity as well as the matrix sorption properties can vary. An efficient method is presented whereby a particle can be followed along the variable property flowpath. For stationary flow conditions and a network of channels with advective flow and matrix diffusion, a simple analytical solution for the residence time distribution along each pathway can be used. Only two parameter groups need to be integrated along each path. For transient flow conditions, a time stepping procedure that incorporates a stochastic Monte-Carlo like method to follow the particles along the paths when flow conditions change is used. The method is fast and an example is used for illustrative purposes. It is exemplified by a case where land rises due to glacial rebound. It is shown that the effects of changing flowrates and directions can be considerable and that the diffusive migration in the matrix can have a dominating effect on the results.

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.

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.

Determination of the flow-wetted surface in fractured media.

Diffusion and sorption in the rock matrix are important retardation mechanisms for radionuclide transport in fractured media. For the conditions existing in a deep repository in crystalline rock, interaction with the rock matrix is controlled by the water flowrate in the fractures and the surface area in contact with the flowing water (the so-called "flow-wetted surface" (FWS)). The flow-wetted surface may be determined from the frequency of open fractures intersecting a borehole. The choice of packer distance used in these hydraulic measurements is crucial, however, since several open fractures may be found in one packer interval. The use of a packer distance that is too large may result in a considerable underestimation of the flow-wetted surface. This is especially important in zones with a high frequency of open fractures (fracture zones) where a small packer distance is a fundamental requirement. A large volume of hydraulic data has been compiled in Sweden from measurements using quite small packer distances. Over the last decade, the most common packer distance used for the hydraulic tests has been 3 m, although some new measurements using a shorter packer distance have also been performed. In several cases, the resolution of these measurements has been less than 0.5 m. All these data have been analysed in detail. From these data, the flow-wetted surface has been calculated and compared with the flow-wetted surface estimated in earlier studies. The results show the importance of using a small packer distance for carrying out borehole transmissivity measurements.

A stochastic multi-channel model for solute transport--analysis of tracer tests in fractured rock.

Some of the basic assumptions of the advection-dispersion model (AD-model) are revisited. This model assumes a continuous mixing along the flowpath similar to Fickian diffusion. This implies that there is a constant dispersion length irrespective of observation distance. This is contrary to most field observations. The properties of an alternative model based on the assumption that individual water packages can retain their identity over long distances are investigated. The latter model is called the multi-channel model (MCh-model). Inherent in the latter model is that if the waters in the different pathways are collected and mixed, the "dispersion length" is proportional to distance. The conditions for when non-mixing between adjacent streams can be assumed are explored. The MCh- and AD-models are found to have very similar residence time distributions (RTD) for Peclet numbers larger than 3. A generalized relation between flowrate and residence time is developed, including the so-called cubic law and constant aperture assumptions. The two models extrapolate very differently when there is strong matrix interaction. The AD-model could severely underestimate the effluent concentration of a tracer pulse and overestimate the residence time. The conditions are explored for when in-filling particles in the fracture will not be equilibrated but will act as if there was seemingly a much larger flow wetted surface (FWS). It is found that for strongly sorbing tracers, relatively small particles can act in this way for systems and conditions that are typical of many tracer tests. The assumption that the tracer residence time found by cautiously injecting a small stream of traced water represents the residence time in the whole fracture is explored. It is found that the traced stream can potentially sample a much larger fraction of the fracture than the ratio between the traced flowrate and the total pumped flowrate. The MCh-model was used to simulate some recent tracer tests in what is assumed to be a single fracture at the Aspö Hard rock laboratory in Sweden. Non-sorbing tracers, HTO and Uranin were used to determine the mean residence time and its variance. Laboratory data on diffusion and sorption properties were used to "predict" the RTD of the sorbing tracers. At least 30 times larger FWS or 1000 times larger diffusion or sorption coefficients would be needed to explain the observed BTCs. Some possible reasons for such behavior are also explored.