The evolution of clay rock/cement interfaces in a cementitious repository for low- and intermediate level radioactive waste
Highlights
► We compared geochemical evolution at host rock/concrete interfaces. ► We developed a set of thermodynamic data, valid for different materials. ► The thermodynamic model contains intrinsically varying cation exchange properties. ► The mineralogical changes at material interfaces are restricted to narrow zones. ► The calculations show massive porosity changes near the interface.
Introduction
In Switzerland a site selection process for geological repositories for spent fuel (SF), vitrified high-level waste (HLW), long-lived intermediate-level waste (ILW) and for low- and intermediate-level radioactive waste (L/ILW) (“Sectoral Plan for Geological Repositories”) is under way. The safety of the SF/HLW and the L/ILW repositories shall be ensured by a multi-barrier concept. Barriers include waste matrices, steel containments (canisters, drums), emplacement containers, cavern backfills and the host rock (Nagra, 2002a).
A L/ILW repository is characterized by rather high waste volumes on one hand and by low radiotoxicities on the other hand. The stored wastes are characterized by a large variability and heterogeneity. Apart from the embedded radioactive wastes, the principal materials in the L/ILW mainly are hydrated cements, aggregate materials and steel. Usually, the wastes are solidified/embedded in a cement/concrete matrix, which itself is included in steel drums. These drums are mainly needed for interim storage, handling and transport prior to emplacement. Parts of the wastes might also be solidified in bituminized matrices. For the construction of the caverns a reinforced concrete or shotcrete wall support is used. Specifically designed concretes (e.g. “monocorn” concrete) may be used for backfilling the remaining cavities. The concrete and mortar aggregates are considered to consist of siliceous materials (primarily quartz), but calcareous materials (lime) are an option. In summary, the backfilled caverns contain, in addition to the wastes, chemically reactive hydrated cements, significant amounts of steel (drums, tunnel support, reinforcement bars and construction materials) and concrete aggregates (Nagra, 2008a).
In Switzerland, the proposed surrounding host rocks are clay-rich sediments (Opalinus Clay, “Brown Dogger”, Effingen Member, “Helvetic Marls”) (Nagra, 2008b). At first glance, all these rocks are characterized by similar solid phase assemblages, by similar pore solution compositions and by the presence of significant fractions of strongly sorbing clay materials.
The knowledge on the temporal and spatial evolution of alterations near interfaces between different materials is important for the performance assessment of deep geological repositories for radioactive waste. The chemical status of the repository, as well as its temporal changes, is driven by the interactions between hydrated cement, steel, aggregates and clays. Considering the chemical nature of the involved major materials, two basic processes take place: an acid–base reaction (aggregate/clays–alkaline cement) on one hand, and the corrosion of iron on the other hand. All chemical reactions solely proceed if a stagnant or mobile aqueous phase is present in sufficient quantity. Chemical gradients, here in particular emanating from the strongly alkaline cementitious materials, dictate the reaction courses, and the principles of chemical thermodynamics determine the fundamental reaction products of the interactions. In terms of transport a stagnant aqueous phase is dominated by diffusion of solutes, whereas for a mobile aqueous phase advective transport may dominate.
Information from literature reviews and mass balance calculations give indications on maximum possible changes of barrier materials. To investigate the spatial and temporal evolution of the overall system, the specific influence of single processes (e.g. pore space changes) and the potential feedback of different processes, numerical models are established research tools.
Research related to cement/clay interactions is, for more than 20 years, a focal point in the field of radioactive waste management.
Gaucher and Blanc (2006) give an overview on publications on experiments, natural analogs and modeling related to cement/clay interactions. Based on their literature review they conclude that it appears premature to conclude definitely on the consequences of an alkaline disturbance in a clay medium. They state that the mineralogical consequences of an alkaline disturbance are fairly well-known, although the knowledge of the thermodynamics of clays, zeolites and cement phases needs to be improved further (see also Savage et al. (2007) on this topic). Large uncertainties remain in the field of dissolution kinetics and mineral precipitation, and especially in understanding the dissolution kinetics of montmorillonites.
Very recently Wang et al. (2010) presented the results of a feasibility study for disposing radioactive waste in a cementitious repository in Boom Clay. They give an extensive overview on work done in Belgium, Switzerland and France on cement–clay host rock interactions, a summary of studies on bentonite–cement interactions and further present results of the ECOCLAY II project. They conclude that the nature of the alkaline plume disturbance in clay materials is relatively well understood.
1D reactive transport simulations typically show that the cement–clay interface most likely clogs and that mass transport across the interface is significantly slowed down. Reactive transport calculations for various clay–cement interfaces show porosity clogging at the interface (Fernández et al., 2009, Marty et al., 2009, Smellie, 1998, Trotignon et al., 2006, Trotignon et al., 2007, De Windt et al., 2004). Trotignon et al., 2006, Trotignon et al., 2007 did numerical modeling on the durability of concrete engineered barriers in contact with mudrock (clay host rock). They think that one key process is the progressive localized cementation of the altered mudrock. Their results suggest that a sharpening of the cementation front will occur and lead to low cementation lengths in the mudrock (∼0.2m). They showed that adding kinetic constraints to mineral dissolution/precipitation does not change the alteration distances significantly, but may influence the mineral assembly to a certain degree. Traber and Mäder (2012) investigated the evolution of an Opalinus Clay–concrete interface with different mineralogical and transport scenarios. They found a porosity reduction at the interface and a limited extent of mineralogical changes (less than 0.2 m).
However, several authors recognize a strong sensitivity of the time behavior of clogging processes on the numerical mesh size and kinetic parameters (Kosakowski et al., 2009, Marty et al., 2009, Traber and Mäder, 2012, Trotignon et al., 2006).
A natural analog for the alteration of clays by high-pH solutions was found in Maqarin (northern Jordan). A very detailed description of the site can be found in Smellie (1998). At the site, hyperalkaline waters, similar to cement waters with a high pH, circulated through fractures in a clayey biomicrite. The fracture edges show dissolution of different minerals, whereas within the fractures different opening/closing episodes may have taken place with complex mineralogical sequences. The extended summary of Smellie (1998) states, that evidence from Maqarin shows
- •
that sequences of minerals predicted by thermodynamic and coupled modeling are similar to those observed in hyperalkaline alteration zones,
- •
that the rock matrix may be accessible to diffusion of aqueous species even during the phase of on-going wallrock alteration,
- •
and that narrow aperture fractures will probably be self-healing.
- •
Although the investigated system is not purely diffusion controlled, the above findings are indications that clogging at clay–cement interfaces is possible and likely over long times.
Gaboreau et al. (2011) characterized the porosity at a concrete–clay–rock interface at the Underground Research Laboratory at Tournemire after 15 years of interaction and could for the first time quantitatively verify porosity clogging. The authors could show that the perturbations (in terms of mineralogical and porosity changes) are limited to 3.5 cm and 1.5–2 cm from the interface in cement and clay rock, respectively. Both, concrete and clay rock are not homogeneous and contain a fissure network. Alterations are observable near the material interface and near the fissures. The porosity in the concrete is increased due to portlandite dissolution, whereas the porosity in the clay rock is strongly reduced, which the authors primarily attribute to the precipitation of C–(A)–S–H phases. The corresponding mineralogical data are published in Techer et al. (2012).
The above mentioned studies show, that the evolution of cement/clay interfaces seem to follow a general pattern. Nevertheless, differences in host rock and concrete transport properties, mineralogy and pore water composition might significantly change the evolution of the interfaces or the extension of mineralogical and pH changes.
We present reactive transport calculations which are specifically tailored to the host rocks and concrete materials foreseen in the Swiss disposal design. Such reactive transport calculations do not only allow refined mass balance calculations, but also introduce a qualitative view onto the system evolution as several processes and their mutual impacts can be considered simultaneously. Some of the presented models include a full coupling between porosity changes due to precipitation/dissolution of minerals and the change of porosity dependent transport parameters (effective diffusion coefficients).
Setting up these coupled models in turn required the refinement of geochemical models for cementitious (concrete) and clay materials and, particularly, the development of a consistent set of thermodynamic data valid for both types of materials. This improved geochemical setup is also used to assess the solubility of radionuclides in the cementitious near field.
Section snippets
Model setup
The model results presented in this work are only a small part of calculations that were done to evaluate the possible geochemical evolution of the L/ILW near-field of a deep geological repository in Switzerland. One task of this study, which also included literature studies and mass balance calculations, was to compare the impact of different host rocks (reflecting different siting regions) on the repository evolution.
To do so, we needed tools and models allowing for a quick and simple
Applications
Our reactive transport calculations concentrated on two scenarios: (1) transport in the repository is dominated by diffusion and (2) transport is dominated by advection.
For each transport scenario three host rocks were considered: Opalinus Clay, Effingen Member and Palfris formation. It is then possible to add up the effects of both transport types to appraise mixed diffusive–advective scenarios. Complex geochemical models in combination with several scenarios provide an overwhelming and hardly
Summary and conclusions
Present reactive transport simulations give insight into the behavior of alteration zones and the progress of reaction fronts. The results are particularly helpful when comparing the impact of different transport conditions (e.g. constant advective or diffusive transport parameters). The predicted maximal spatial extents of mineral alterations and pH increases in the host rocks are remarkably close between different studies, even if they are based on different discretization, chemical setups
Acknowledgments
Partial financial support by the Swiss National Cooperative for the Disposal of Radioactive Waste (NAGRA) is gratefully acknowledged.
References (63)
- et al.
Multicomponent diffusion of a suite of tracers (HTO, Cl, Br, I, Na, Sr, Cs) in a single sample of Opalinus Clay
Geochim. Cosmochim. Acta
(2010) Evolution of pore water chemistry during degradation of cement in a radioactive waste repository environment
Waste Manage. (Oxford)
(1992)- et al.
Ion equilibrium between montmorillonite interlayer space and an external solution—Consequences for diffusional transport
Geochim. Cosmochim. Acta
(2009) - et al.
Long-term corrosion of two nuclear waste reference glasses (MW and SON68): a kinetic and mineral alteration study
Appl. Geochem.
(2006) - et al.
15 Years of in situ cement–argillite interaction from Tournemire URL: characterisation of the multi-scale spatial heterogeneities of pore space evolution
Appl. Geochem.
(2011) - et al.
Cement/clay interactions – a review: experiments, natural analogues, and modeling
Waste Manage.
(2006) - et al.
Consistent interpretation of the results of through-, out-diffusion and tracer profile analysis for trace anion diffusion in compacted montmorillonite
J. Contam. Hydrol.
(2011) - et al.
Self-diffusion of water and its dependence on temperature and ionic strength in highly compacted montmorillonite, illite and kaolinite
Appl. Geochem.
(2008) - et al.
Diffusion of ionic tracers in the Callovo–Oxfordian clay–rock using the Donnan equilibrium model and the formation factor
Geochim. Cosmochim. Acta
(2009) - et al.
Experimentally determined swelling pressures and geochemical interactions of compacted Wyoming bentonite with highly alkaline solutions
Phys. Chem. Earth, Parts A/B/C
(2007)
A thermodynamic approach to the hydration of sulphate-resisting Portland cement
Waste Manage. (Oxford)
Thermodynamic modelling of the hydration of Portland cement
Cem. Concr. Res.
Thermodynamic modelling of the effect of temperature on the hydration and porosity of Portland cement
Cem. Concr. Res.
Influence of reaction kinetics and mesh refinement on the numerical modelling of concrete/clay interactions
J. Hydrol.
Thermodynamic properties of Portland cement hydrates in the system CaO–Al2O3–SiO2–CaSO4–CaCO3–H2O
Cem. Concr. Res.
Solubility of Fe–ettringite (Ca6[Fe(OH)6]2(SO4)3·26H2O)
Geochim. Cosmochim. Acta
Biogeochemical processes in a clay formation in situ experiment: Part E – equilibrium controls on chemistry of pore water from the Opalinus Clay, Mont Terri Underground Research Laboratory, Switzerland
Appl. Geochem.
Water and ions in clays: unraveling the interlayer/micropore exchange using molecular dynamics
Geochim. Cosmochim. Acta
Alteration of bentonite by hyperalkaline fluids: a review of the role of secondary minerals
Phys. Chem. Earth Parts A/B/C
Modeling reactive transport in non-ideal aqueous–solid solution system
Appl. Geochem.
Multicomponent reactive transport in discrete fractures – II: infiltration of hyperalkaline groundwater at Maqarin, Jordan, a natural analogue site
J. Hydrol.
Tracing interactions between natural argillites and hyper-alkaline fluids from engineered cement paste and concrete: chemical and isotopic monitoring of a 15-years old deep-disposal analogue
Appl. Geochem.
Comparison of performance of concrete barriers in a clayey geological medium
Phys. Chem. Earth Parts A/B/C
Predicting the long term durability of concrete engineered barriers in a geological repository for radioactive waste
Phys. Chem. Earth Parts A/B/C
On the relationship between porosity and electrical resistivity in cementitious systems
Cem. Concr. Res.
Thermodynamics of silicates
Physico-Chemical Characterisation and Calculated in Situ Porewater Chemistries for a Low Permeability Palfris Marl Sample from Wellenberg, PSI Report 94-19
Porewater Chemistry in Compacted Re-saturated MX-80 Bentonite: Physico-chemical Characterisation and Geochemical Modelling, PSI Report 02–10
Cited by (69)
Phase changes in cementitious materials exposed to saline solutions
2023, Cement and Concrete ResearchA systematic review of bentonite/concrete interaction system in HLW disposal repositories: Theoretical, experimental and geochemical modelling analysis
2022, Construction and Building MaterialsCitation Excerpt :The effective diffusion coefficient in De Windt et al. [138] is related to the porosity evolution according to the Archie’s law, which is similar with that in Trotignon et al. [143] and Marty et al. [159]. In Kosakowski and Berner [153], the porosity was not updated associated with the mineral transformation to give a constant water flow. That is the advective cases did not adopt a full coupling between mass transport, mineral alteration and porosity change.
Cement‑carbonate rock interaction under saturated conditions: From laboratory to modeling
2022, Cement and Concrete ResearchA new operator-splitting finite element scheme for reactive transport modeling in saturated porous media
2022, Computers and GeosciencesInverse chemical equilibrium problems: General formulation and algorithm
2022, Chemical Engineering Science