Solute transport properties of compacted Ca-bentonite used in FEBEX project

The present Spanish concept of a deep geological high level waste repository includes an engineered clay barrier around the canister. The clay presents a very high sorption capability for radionuclides and a very small hydraulic conductivity, so that the migration process of solutes is limited by sorption and diffusion processes. Therefore, diffusion and distribution coefficients in compacted bentonite (i.e. in "realistic" liquid to solid ratio conditions) are the main parameters that have to be obtained in order to characterise solute transport that could be produced after the canister breakdown. Through-Diffusion (TD) and In-Diffusion (ID) experiments with HTO, Sr, Cs and Se were carried out using compacted FEBEX bentonite, which is the reference material for the Spanish concept of radioactive waste disposal. Experiments were interpreted by means of available analytical solutions that allow the estimation of diffusion coefficients and, in some cases, distribution coefficients. Analytical solutions are simple to use, but rely on hypotheses that do not hold in all the experiments. These experiments were interpreted also using an automatic parameter estimation code that overcomes the limitations of analytical solutions. Numerical interpretation allows the simultaneous estimation of porosity, diffusion and distribution coefficients, accounts for the role of porous sinters and time-varying boundary concentrations, and can use different types of raw concentration data.     

Consistent interpretation of the results of through-, out-diffusion and tracer profile analysis for trace anion diffusion in compacted montmorillonite

Literature data for anion diffusion in compacted swelling clays contain systematic inconsistencies when the results of through-diffusion tests are compared with those of out-diffusion or tracer profile analysis. In the present work we investigated whether these inconsistencies can be explained by taking into account heterogeneities in the compacted samples; in particular increased porosities at the clay boundaries. Based on the combined results of out-diffusion, tracer profile analysis and the spatial distribution of the electrolyte anion in the clay, we conclude that the inconsistencies can indeed be resolved by taking into account a heterogeneous distribution of the total and the anion-accessible porosity. This, by definition, leads to a position dependence of the effective diffusion coefficient. Neglecting these effects results in a rather subordinate systematic error in the determination of effective diffusion coefficients of anions from through-diffusion tests with clay thicknesses in the centimetre range. However, stronger errors in terms of absolute values and conceptual interpretation may be introduced in out-diffusion tests and profile analyses of the diffused tracer. We recommend that anion diffusion tests should be accompanied by measurements of the total and anion-accessible porosity as a function of position in the direction of diffusion.     

Diffusion of radionuclides in concrete and concrete-bentonite systems

Various construction materials are under consideration for nuclear waste repositories. Two important materials are concrete and bentonite clay, which will act as mechanical barriers and prevent convective water flow. These barriers will also retard transport (diffusion controlled) of dissolved radionuclides by a combination of mechanical constraints and chemical interactions with the solid.An important issue is the possible change of the initial sodium bentonite into the calcium form due to interaction with calcium from the concrete. The initial leaching of concrete was studied using radioactive spiked concrete in contact with compacted bentonite.Measurement were made of the diffusion of Cs, Am and Pu into 5 different types of concrete in contact with pore water. The diffusivity measured for Cs agrees reasonably well with data found in the literature. No movement could be measured for Am and Pu (< 0.2 mm), even though the contact times were extremely long (2.5 and 5 yr, respectively). The diffusion of Na, Ca and Cs from concrete into bentonite was also measured.     

Sensitivity analysis of radionuclide migration in compacted bentonite: a mechanistic model approach

Mechanistic model calculations for the migration of Cs, Ra, Am and Pb in compacted bentonite have been carried out to evaluate sensitivities with respect to different parameter variations. A surface chemical speciation/electric double layer model is used to calculate: (i) porewater composition and radionuclide speciation in solution and at the bentonite surface, yielding the distribution of mobile and sorbed species and (ii) interaction of diffusing species with negatively charged pore walls to obtain diffusion parameters. The basic scenario considers the interaction of compacted bentonite with a fresh-type groundwater; variations include the presence of bentonite impurities and saline groundwater. It is shown that these scenarios result in significant variations of porewater composition that affect migration via three mechanisms that can partly compensate each other: (1) effects on sorption through radionuclide complexation in solution, and competition of major cations for surface sites; (2) changes in radionuclide solution speciation leading to different diffusing species under different conditions; (3) effects on diffusion through changes in the electric double layer properties of the clay pores as a function of ionic strength.     

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.     

Generation and stability of bentonite colloids at the bentonite/granite interface of a deep geological radioactive waste repository

The possible mechanisms of colloid generation at the near field/far field interface of a radioactive repository have been investigated by means of novel column experiments simulating the granite/bentonite boundary, both in dynamic and in quasi-static water flow conditions. It has been shown that solid particles and colloids can be detached from the bulk and mobilised by the water flow. The higher the flow rate, the higher the concentration of particles found in the water, according to an erosion process. However, the gel formation and the intrinsic tactoid structure of the clay play an important role in the submicron particle generation even in the compacted clay and in a confined system. In fact, once a bentonite gel is formed, in the regions where the clay is contacted with water, clay colloids can be formed even in quasi-static flow conditions. The potential relevance of these colloids in radionuclide transport has been studied by evaluating their stability in different chemical environments. The coagulation kinetics of natural bentonite colloids was experimentally studied as a function of the ionic strength and pH, by means of time-resolved light scattering techniques. It has been shown that these colloids are very stable in low saline (approximately 1 x 10(-3) M) and alkaline (pH > or = 8) waters.     

Simultaneous estimation of effective and apparent diffusion coefficients in compacted bentonite

Effective diffusion coefficients (D(e)) are usually measured by means of "through-diffusion" experiments in which steady state is reached, and the "time-lag" methods are used to estimate the apparent diffusion coefficient (D(a)). For sorbing radionuclides (as caesium), the time needed to reach steady-state conditions is very large, and the precision in D(a) determinations is not satisfactory. In this paper, a method that allows determining simultaneously effective and apparent diffusion coefficients in compacted bentonite without reaching steady-state conditions is described. Basically, this method consists of an "in-diffusion" experiment in which the concentration profile in the bentonite sample is used to estimate D(a), and the temporal evolution of the solute concentration in the reservoir is used to estimate D(e). This method has several advantages over the typical "through-diffusion" experiments, in particular: (a) experiment duration is significantly shorter, (b) D(a) values are measured with greater precision and (c) it is not necessary to maintain a constant solute concentration in the reservoir. This new method has been used to estimate the effective and apparent diffusion coefficients for caesium in FEBEX bentonite and in order to validate it, the results have been compared with results previously obtained with standard methods. Experimental results have been satisfactorily modelled using a simple model of diffusion in porewater and the measured value of D(e)(Cs) is very similar to D(e)(HTO) in the same bentonite. There is no evidence of "surface diffusion" in FEBEX bentonite for caesium.     

Study of the contaminant transport into granite microfractures using nuclear ion beam techniques

Hydrated bentonite is a very plastic material and it is expected to enter in the rock microfractures at the granite/bentonite boundary of a deep geological high-level waste repository. This process is enhanced by the high swelling pressure of the clay. Since bentonite has a very good sorption capability for many radionuclides, the displacement of the clay might lead to a "clay-mediated" contaminant transport into the rock. The aim of this work is to study the contaminant transport into granite microfractures using nuclear ion beam techniques, and to determine to what extent the clay can favour it. To do so, bentonite previously doped with uranium, cesium and europium was put in contact with the surface of granite sheets. Granite sheets contacted with non-doped bentonite and with radionuclide solutions were also prepared as references. This allowed analysing the differences in the diffusion behaviour of the three systems: clay, radionuclides and clay plus radionuclides. A combination of Rutherford backscattering spectrometry (RBS) and other nuclear ion-beam techniques such as particle-induced X-ray emission (PIXE) and microPIXE was used to study the depth and lateral distribution of clay and contaminants inside granite. It was also tried to evaluate not only the diffusion depth and diffusion coefficients but also the different areas of the granite where the diffusants have a preferential access.     

Inverse modeling of multicomponent reactive transport through single and dual porosity media

Compacted bentonite is foreseen as buffer material for high-level radioactive waste in deep geological repositories because it provides hydraulic isolation, chemical stability, and radionuclide sorption. A wide range of laboratory tests were performed within the framework of FEBEX (Full-scale Engineered Barrier EXperiment) project to characterize buffer properties and develop numerical models for FEBEX bentonite. Here we present inverse single and dual-continuum multicomponent reactive transport models of a long-term permeation test performed on a 2.5 cm long sample of FEBEX bentonite. Initial saline bentonite porewater was flushed with 5.5 pore volumes of fresh granitic water. Water flux and chemical composition of effluent waters were monitored during almost 4 years. The model accounts for solute advection and diffusion and geochemical reactions such as aqueous complexation, acid-base, cation exchange, protonation/deprotonation by surface complexation and dissolution/precipitation of calcite, chalcedony and gypsum. All of these processes are assumed at local equilibrium. Similar to previous studies of bentonite porewater chemistry on batch systems which attest the relevance of protonation/deprotonation on buffering pH, our results confirm that protonation/deprotonation is a key process in maintaining a stable pH under dynamic transport conditions. Breakthrough curves of reactive species are more sensitive to initial porewater concentration than to effective diffusion coefficient. Optimum estimates of initial porewater chemistry of saturated compacted FEBEX bentonite are obtained by solving the inverse problem of multicomponent reactive transport. While the single-continuum model reproduces the trends of measured data for most chemical species, it fails to match properly the long tails of most breakthrough curves. Such limitation is overcome by resorting to a dual-continuum reactive transport model.     

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.     

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.     

Simultaneous sorption of lead and chlorobenzene by organobentonite

Clays or organoclays have been used as a barrier to prevent the transport of hazardous contaminants in landfills. However, clays are known to effectively sorb mostly inorganic contaminants, while organoclays are mainly used for organic contaminants. Since the organoclays are basically clay particles modified with cationic surfactants, there might exist an optimal coverage of cationic surfactant on the clay particles to sorb both inorganic and organic contaminants. In order to determine the optimal mass of cationic surfactants on the bentonites, sodium bentonites were treated with various ratios of hexadecyltrimethylammonium (HDTMA) to bentonites. Chlorobenzene and lead were selected as representative contaminants. When either chlorobenzene or lead exists as a single contaminant, chlorobenzene sorption increased with increasing HDTMA to bentonite ratios, and lead sorption decreased with increasing HDTMA to bentonite ratios. Sorption of chlorobenzene was a function of HDTMA coverage on the bentonites, while lead sorption was much more influenced by the initial lead concentration rather than the mass of HDTMA added to the bentonites.     

Diffusive transport through compacted Na- and Ca-bentonite

The effect of exchangeable cation — Na⁺ and Ca ²⁺ — on the diffusive transport of I⁻, Sr ²⁺ and ³H (as HTO) in compacted bentonite was examined using a through-diffusion method. Total intrinsic diffusion coefficients, D_i, were determined from the steady-state flux of the diffusants through the clays, and apparent diffusion coefficients, D_a, were obtained from the time lag technique. The clays were compacted to a dry bulk density of 1.3 Mg/m³, and Na-bentonite was saturated with a solution of 100 mol NaCl/m3 and Ca-bentonite with one of 50 mol CaCl₂/m³. The D_i values for all diffusants are 2 to 6 times higher in the Ca- than Na-clay. We attribute this to the larger quasicrystal, or particle, size of Ca- compared to Na-bentonite. Hence, Ca-bentonite has a greater proportion of relatively large pores; this was confirmed by Hg intrusion porosimetry. This means the diffusion pathways in Ca-bentonite are less tortuous than those in Na-bentonite. Moreover, in some cases the effective porosity, or the porosity available for diffusive transport, may be greater in Ca-bentonite. The D_a values are inversely proportional to the distribution coefficients of the diffusants with the clays.     

Porewater chemistry in compacted re-saturated MX-80 bentonite

Bentonites of various types are being investigated in many countries as backfill materials in high-level radioactive waste disposal concepts. Being able to understand the chemistry of the porewater in compacted bentonite is very important since it is critical to predicting radionuclide solubilities and to the synthesis of sorption data bases, and hence to repository safety studies. In this paper, porewater compositions in compacted bentonites are calculated, taking into consideration such factors as montmorillonite swelling, semi-permeable membrane effects, very low "free water" volumes, and the highly effective buffering characteristics of the exchangeable cations and the amphoteric edge sites. The former buffer the cation concentrations and the latter fix the pH in the porewater of a re-saturated bentonite. The above considerations are used in conjunction with previously measured physico-chemical characterisation data on MX-80 powder to calculate porewater compositions in compacted bentonites. For the MX-80 material specified, the porewaters calculated for initial dry densities between 1200 and 1600 kg m(-3) had relatively high ionic strengths (I approximately 0.3 M), similar cation concentrations and a pH equal to 8.0. The porewaters changed from being Na(2)SO(4)-rich at 1200 kg m(-3) to a NaCl/Na(2)SO(4) type water at 1600 kg m(-3).     

An integrated sorption-diffusion model for the calculation of consistent distribution and diffusion coefficients in compacted bentonite

A thermodynamic sorption model and a diffusion model based on electric double layer (EDL) theory are integrated to yield a surface chemical model that treats porewater chemistry, surface reactions, and the influence of charged pore walls on diffusing ions in a consistent fashion. The relative contribution of Stern and diffuse layer to the compensation of the permanent surface charge represents a key parameter; it is optimized for the diffusion of Cs in Kunipia-F bentonite, at a dry density of 400 kg/m3. The model is then directly used to predict apparent diffusivities (Da) of Cs, Sr, Cl-, I- and TcO4- and corresponding distribution coefficients (Kd) of Cs and Sr in different bentonites as a function of dry density, without any further adjustment of surface chemical and EDL parameters. Effective diffusivities (De) for Cs, HTO, and TcO4- are also calculated. All calculated values (Da, De, Kd) are fully consistent with each other. A comparison with published, measured data shows that the present model allows a good prediction and consistent explanation of (i) apparent and effective diffusivities for cations, anions, and neutral species in compacted bentonite, and of (ii) Kd values in batch and compacted systems.     

Experimental and modeling studies on sorption and diffusion of radium in bentonite

The sorption and desorption behavior of radium on bentonite and purified smectite was investigated as a function of pH, ionic strength and liquid to solid ratio by batch experiments. The distribution coefficients (Kd) were in the range of 10(2) to > 10(4) ml g-1 and depended on ionic strength and pH. Most of sorbed Ra was desorbed by 1 M KCl. The results for purified smectite indicated that Ra sorption is dominated by ion exchange at layer sites of smectite, and surface complexation at edge sites may increase Ra sorption at higher pH region. Reaction parameters between Ra and smectite were determined based on an interaction model between smectite and groundwater. The reaction parameters were then used to explain the results of bentonite by considering dissolution and precipitation of minerals and soluble impurities. The dependencies of experimental Kd values on pH, ionic strength and liquid to solid ratio were qualitatively explained by the model. The modeling result for bentonite indicated that sorption of Ra on bentonite is dominated by ion exchange with smectite. The observed pH dependency was caused by changes of Ca concentration arising from dissolution and precipitation of calcite. Diffusion behavior of Ra in bentonite was also investigated as a function of dry density and ionic strength. The apparent diffusion coefficients (Da) obtained in compacted bentonite were in the range of 1.1 x 10(-11) to 2.2 x 10(-12) m2 s-1 and decreased with increasing in dry density and ionic strength. The Kd values obtained by measured effective diffusion coefficient (De) and modeled De were consistent with those by the sorption model in a deviation within one order of magnitude.     

Kinetics of metribuzin release from bentonite-polymer composites in water

A series of bentonite polymer-composites (BPCs) loaded with metribuzin were studied for their controlled release in aqueous medium. The release of active ingredient from BPCs was significantly lower as compared to commercial metribuzin formulation. The results revealed that the cumulative metribuzin release was highest (81%) from the BPCs containing 8% clay (commercial bentonite) and 2% metribuzin which correspond to the lowest (14 days) half-life values i.e., time required for 50% release of active ingredient (t_1/2). The metribuzin release from the BPCs decreased with increased concentration of clays in polymer matrix and the release was further decreased with BPCs prepared with pure nano-bentonite. BPCs containing 12% clay and 2% metribuzin showed maximum t_1/2 values i.e., 25 and 51 days for commercial bentonite and pure nano-bentonite as clay sources, respectively. The differential behaviour in the metribuzin release rates from BPCs was ascribed due to variations in crosslinking of metribuzin in the composites. As metribuzin release was found to be slower in BPCs compared to commercial formulation, it could be used for control of weeds tailored to different crops.     

Diffusion mechanism of chloride ions in sodium montmorillonite

For safety assessment of geological disposal of HLW, it is necessary to understand the diffusion mechanism of radionuclides in compacted bentonite. In this study, the diffusion behavior of chloride ions in compacted montmorillonite was studied from the viewpoints of the activation energy for apparent diffusion and the basal spacing of the compacted montmorillonite. A unique change in the activation energy as a function of the dry density of the montmorillonite was found. The activation energy decreased from 17.4 to 13.5 kJ mol-1 as the dry density increased from 0.7 to 1.0 Mg m-3 and increased to 25.1 kJ mol-1 at dry densities above 1.0 Mg m-3. The basal spacing of 1.88 nm, corresponding to the three-water layer hydrate state, was not observed by X-ray diffraction (XRD) until the dry density increased to 1.0 Mg m-3, where the minimum activation energy was obtained. On the other hand, a basal spacing of 1.56 nm, corresponding to the two-water layer hydrate state, was observed at the dry densities of 1.4 Mg m-3 and above, where the activation energies were more than 22 kJ mol-1. These experimental results suggest that there are at least two additional diffusion processes that can raise or reduce the activation energy and are affected by water in the region adjacent to the montmorillonite surfaces. If the "Grahame model" can be introduced to describe the electrical double layer, surface diffusion will be considered the possible predominant diffusion process, even for anions like chloride ions.