EXTERNAL REPORT SCK•CEN-ER-198 14/Ssa/P-16
Speciation and solubility calculations for waste relevant radionuclides in Boom Clay
First Full Draft
Sonia Salah and Lian Wang
SCK•CEN Contract: CO-90-08-2214-00, RP.W&D.0064 NIRAS/ONDRAF contract: CCHO- 2009-0940000, LTBC02-GEO-01 Radionuclide migration and retention processes in Boom Clay
April, 2014
SCK•CEN RDD Boeretang 200 BE-2400 Mol Belgium
EXTERNAL REPORT OF THE BELGIAN NUCLEAR RESEARCH CENTRE SCK•CEN-ER-198 14/Ssa/P-16
Speciation and solubility calculations for waste relevant radionuclides in Boom Clay
First Full Draft
Sonia Salah and Lian Wang
SCK•CEN Contract: CO-90-08-2214-00, RP.W&D.0064 NIRAS/ONDRAF contract: CCHO- 2009-0940000, LTBC02-GEO-01 Radionuclide migration and retention processes in Boom Clay
April, 2014 Status: Unclassified ISSN 1782-2335
SCK•CEN Boeretang 200 BE-2400 Mol Belgium
© SCK•CEN Studiecentrum voor Kernenergie Centre d’étude de l’énergie Nucléaire Boeretang 200 BE-2400 Mol Belgium
Phone +32 14 33 21 11 Fax +32 14 31 50 21 http://www.sckcen.be
Contact: Knowledge Centre [email protected]
COPYRIGHT RULES
All property rights and copyright are reserved to SCK•CEN. In case of a contractual arrangement with SCK•CEN, the use of this information by a Third Party, or for any purpose other than for which it is intended on the basis of the contract, is not authorized. With respect to any unauthorized use, SCK•CEN makes no representation or warranty, expressed or implied, and assumes no liability as to the completeness, accuracy or usefulness of the information contained in this document, or that its use may not infringe privately owned rights. SCK•CEN, Studiecentrum voor Kernenergie/Centre d'Etude de l'Energie Nucléaire Stichting van Openbaar Nut – Fondation d'Utilité Publique ‐ Foundation of Public Utility Registered Office: Avenue Herrmann Debroux 40 – BE‐1160 BRUSSEL Operational Office: Boeretang 200 – BE‐2400 MOL
Abstract
The Boom Clay formation represents the reference host rock for the geological disposal of high- level and/or long-lived radioactive waste (HLW-LL) in Belgium. The Belgian authority responsible for the long-term management of radioactive waste is ONDRAF/NIRAS (Organisme National des Déchets Radioactifs et Matières Fissiles Enrichies/Nationale Instelling voor Radioactief Afval en Verrijkte Splijstoffen). The current focus of the research programme of ONDRAF/NIRAS (O/N) is the so-called Safety and Feasibility Case 1 (SFC-1), which represents a body of sound arguments and evidences describing, quantifying and substantiating that geological disposal of high-level and/or long-lived radioactive waste is a safe and feasible long- term solution.
In order to calculate the release and transport of the waste relevant radionuclides (RNs) through the near-field of the repository and the host clay layer (far-field), different parameters, such as concentration/solubility limits, retardation factors, diffusion accessible porosity values and pore diffusion coefficients are needed.
In this report, the solubility assessment methodology, i.e. results of the speciation calculations and derivation of the RN concentration/solubility limits applicable for "undisturbed" Boom Clay conditions (far-field) will be presented.
The solubility and speciation modelling was performed with the geochemical computer code The Geochemist's Workbench (versions 8.08, 8.10 and 8.12). The reference thermodynamic database that was used for the calculations and also developed at SCK•CEN is named MOLDATA thermodynamic database (2010_MOLDATA_nov_b.dat; MOLDATA TDB, version 1).
5
Table of Contents ABSTRACT ...... 5 TABLE OF CONTENTS ...... 6 LIST OF FIGURES ...... 8 LIST OF TABLES ...... 9 STRUCTURE OF THE DOCUMENT ...... 10 1 INTRODUCTION ...... 11 2 SYSTEM DEFINITION ...... 12
2.1 THE DISPOSAL SYSTEM ...... 12 2.2 BOOM CLAY REFERENCE PORE WATER ...... 12 2.3 COMPUTER CODE AND DATABASES ...... 13 2.4 MOLDATA ...... 13 2.5 GDP ...... 14 3 SPECIATION CALCULATIONS AND POURBAIX DIAGRAMS ...... 16
3.1 SPECIES ACTIVITIES AND ACTIVITY COEFFICIENTS ...... 17 4 SOLUBILITY ‐ DEFINITION AND THEORETICAL BACKGROUND ...... 19
4.1 SOLUBILITY CONSTANTS...... 19 4.2 SOLUBILITY AND SATURATION ...... 20 4.2.1 Calculation procedure ...... 21 4.2.2 Reasoning and approach of solid phase selection ...... 21 4.2.3 Solubility Source and Expert ranges ...... 23 4.2.4 Uncertainties ...... 23 4.3 PARAMETERS INFLUENCING SOLUBILITY ...... 25 4.3.1 Influence of particle size ...... 25 4.3.2 Ostwald rule and ripening ...... 25 4.3.3 Influence of temperature and pressure ...... 26 4.3.4 Influence of ionic strength ...... 26 4.3.5 Influence of inorganic complexation ...... 26 4.3.6 Dissolved Organic Carbon in BC porewater ...... 27 4.3.7 Influence of organic complexation ...... 28 4.3.8 Influence of colloids ...... 28 4.3.9 Influence of Eh, pH and pCO2 ...... 30 4.3.10 Influence of radiation ...... 30 4.4 THE ROLE OF SOLID SOLUTIONS ...... 31 4.4.1 Drawbacks for the application of solid solution models ...... 32 5 RESULTS ...... 33
5.1 ACTINIUM (AC) ...... 35 5.2 AMERICIUM (AM) ...... 36 5.3 BERYLLIUM (BE) ...... 41 5.4 CALCIUM (CA) ...... 43 5.5 CAESIUM (CS) ...... 45 5.6 CARBON (C) ...... 46 5.7 CHORINE (CL) ...... 48 5.8 CURIUM (CM) ...... 49 5.9 IODINE (I) ...... 50 5.10 MOLYBDENUM (MO) ...... 53 5.11 NEPTUNIUM (NP) ...... 57 5.12 NICKEL (NI) ...... 61 5.13 NIOBIUM (NB) ...... 66 5.14 PALLADIUM (PD) ...... 69
6
5.15 PLUTONIUM (PU) ...... 73 5.15.1 The redox behaviour of tri‐ and tetravalent (Np and) Pu ...... 78 5.16 PROTACTINIUM (PA) ...... 80 5.17 RADIUM (RA) ...... 83 5.18 RUBIDIUM (RB) ...... 87 5.19 SAMARIUM (SM) ...... 88 5.20 SELENIUM (SE) ...... 92 5.21 SILVER (AG) ...... 97 5.22 STRONTIUM (SR) ...... 102 5.23 TECHNETIUM (TC) ...... 104 5.24 THORIUM (TH) ...... 108 5.25 TIN (SN) ...... 113 5.26 URANIUM (U) ...... 117 5.27 ZIRCONIUM (ZR) ...... 122 6 ANNEX I: SUMMARY OF SOLUBILITY CALCULATIONS FOR PHASES COMPRISED IN THE SR & ER RANGES .... 126 7 ANNEX II: SUMMARY OF SOURCE AND EXPERT RANGES* ...... 128 8CALCULATION ANNEX III: OF THERMODYNAMIC UNCERTAINTIES...... 130 9 ANNEX IV: TECHNICAL NOTE ON TH‐DATA, I.E. BINARY TH‐CARBONATE AND TERNARY TH‐HYDROXO‐ CARBONATE COMPLEXES ...... 140 10 REFERENCES ...... 146
7
List of Figures Figure 1: Sulphur and carbon speciation as function of Eh and pH ...... 17 Figure 2: Eh-pH diagram of americium (Am-C-S-O-H) for the BC reference porewater system...... 37 Figure 3: Eh-pH diagram of beryllium (Be-C-S-O-H) for the BC reference porewater system...... 41 Figure 4: Eh-pH diagram of calcium (Ca-C-S-O-H) for the BC reference porewater system...... 43 Figure 5: Eh-pH diagram of caesium (Cs-C-S-O-H) for the BC reference porewater system...... 45 Figure 6: Eh-pH diagram of carbon (C-S-O-H) for the BC reference porewater system...... 46 Figure 7: Eh-pH diagram of chlorine (Cl-C-S-O-H) for the BC reference porewater system...... 48 Figure 8: Eh-pH diagram of iodine (I-C-S-O-H) for the BC reference porewater system...... 51 Figure 9: Eh-pH diagram of molybdenum (Mo-C-S-O-H) for the BC reference porewater system...... 54 Figure 10: Eh-pH diagram of neptunium (Np-C-S-O-H) for the BC reference porewater system...... 58 Figure 11: Eh-pH diagram of nickel (Ni-C-S-O-H) for the BC reference porewater system...... 62 Figure 12: Eh-pH diagram of niobium (Nb-C-S-O-H) for the BC reference porewater system...... 66 Figure 13: Eh-pH diagram of palladium (Pd-C-S-O-H) for the BC reference porewater system...... 70 Figure 14: Eh-pH diagram of plutonium (Pu-C-S-O-H) for the BC reference porewater system...... 74 Figure 15: Comparison of Np(III)/Np(IV) and Pu(III)/Pu(IV) speciation...... 78 Figure 16: Eh-pH diagram of protactinium (Pa-C-S-O-H) for the BC reference porewater system...... 80 Figure 17: Eh-pH diagram of radium (Ra-C-S-O-H) for the BC reference porewater system...... 83 Figure 18: Eh-pH diagram of rubidium (Rb-C-S-O-H) for the BC reference porewater system...... 87 Figure 19: Eh-pH diagram of samarium (Sm-C-S-O-H) for the BC reference porewater system...... 89 Figure 20: Eh-pH diagram of selenium (Se-C-S-O-H) for the BC reference porewater system...... 93 Figure 21: Eh-pH diagram of silver (Ag-C-S-O-H) for the BC reference porewater system...... 98 Figure 22: Silver speciation as function of total dissolved Ag concentration ...... 100 Figure 23: Eh-pH diagram of strontium (Sr-C-S-O-H) for the BC reference porewater system...... 102 Figure 24: Eh-pH diagram of technetium (Tc-C-S-O-H) for the BC reference porewater system...... 105 Figure 25: Eh-pH diagram of thorium (Th-C-S-O-H) for the BC reference porewater system...... 109 Figure 26: Eh-pH diagram of tin (Sn-C-S-O-H) for the BC reference porewater system...... 114 Figure 27: Eh-pH diagram of uranium (U-C-S-O-H) for the BC reference porewater system...... 119 Figure 28: Eh-pH diagram of zirconium (Zr-C-S-O-H) for the BC reference porewater system...... 123 Figure 29: Eh-pH diagram of thorium (Th-C-S-O-H) for the BC reference porewater system. Assumed activity of dissolved [Th] = 10-8. Database: MOLDATA. Code: The Geochemist's Workbench - 8.10 ...... 140
8
List of Tables Table 1: The (recalculated) Boom Clay water composition at 25°C ...... 13 Table 2: Solubility of Am in the BC reference porewater system at 25°C, pH 8.355 and Eh -281 mV...... 38 Table 3: Species distribution of Am in equilibrium with AmCO3OH(am,hyd)...... 38 Table 4: Solubility of Be in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 42 Table 5: Solubility of Ca in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 43 Table 6: Solubility of Mo in the BC reference porewater system at pH 8.355 and Eh -281 mV...... 54 Table 7: Solubility of Np in the BC reference porewater system at 25°C, pH 8.355 and Eh -281 mV...... 59 Table 8: Solubility of Ni in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 63 Table 9: Species distribution in equilibrium with Ni(OH)2(beta)...... 64 Table 10: Solubility of Nb in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 67 Table 11: Species distribution of Nb in equilibrium with Nb2O5(cr)...... 67 Table 12: Solubility of Pd in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 71 Table 13: Solubility of Pu in the BC reference porewater system at 25°C, pH 8.355 and Eh -281 mV...... 75 Table 14: Species distribution of Pu in equilibrium with PuO2(am,hyd)...... 75 Table 15: Solubility of Pa in the BC reference porewater system at 25°C, pH 8.355 and Eh -281 mV...... 80 Table 16: Solubility of Ra in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 84 Table 17: Species distribution of Ra in equilibrium with RaSO4(s)...... 84 Table 18: Solubility of Sm in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 89 Table 19: Solubility of Se in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 94 Table 20: Solubility of Ag in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 99 Table 21: Species distribution of Ag in equilibrium with AgCl(cr)...... 99 Table 22: Solubility of Sr in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 103 Table 23: Solubility of Tc in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 106 Table 24: Solubility of Th in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 110 Table 25: Solubility of Sn in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 115 Table 26: Species distribution of Sn in equilibrium with SnO2(am)...... 115 Table 27: Solubility of U in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 119 Table 28: Species distribution of U in equilibrium with UO2(am,hyd)...... 120 Table 29: Solubility of Zr in the BC reference porewater system at 25 °C, pH 8.355 and Eh -281 mV...... 123 Table 30: Mixed Th-hydroxo-carbonate complexes comprised in MOLDATA ...... 140 Table 31: Solubilty constants of Th(OH)4(am) in carbonate solution (log10K°s,1yz) and derived formation constants (logβ°1yz) ...... 142 Table 32: Recalculated formation constants of ternary Th(IV) complexes ...... 143 Table 33: Formation constants of ternary Th(IV) complexes recalculated by NEA ...... 144 Table 34: Recalculated formation constants logβ*1yz of Th(IV) complexes using the bicarbonate ...... 144
9
Structure of the document The present report comprises 5 main chapters and 4 appendices, which contents are briefly summarized in the following:
Chapter 1: Introduction
Chapter 2: In this chapter, first the Belgian disposal concept and Boom Clay porewater chemistry are briefly described. Afterwards, an overview about the geochemical code and the different thermodynamic databases that were used for the solubility and speciation calculations is given. Addionally, the in- house developed Geochemical Database Processor (GDP), that was used to compile the SCK•CEN reference TDB MOLDATA is presented.
Chapter 3: In this chapter, "speciation" is defined and the theoretical background for the calculation and interpretation of Pourbaix diagrams is given.
Chapter 4: In this chapter, first the solubility term as it is used in this report is defined, and then the thermodynamic background of solubility constants and saturation states is given. In the following, the solubility concept, the calculation procedure, the identification and selection of the solubility limiting solids are described. Afterwards, Source and Expert Ranges (SR and ER), and Best Estimate (BE) values are defined, and the uncertainty calculation approach that was used is illustrated. At the end of this chapter, the effects of different parameters (e.g. Eh, pH, pCO2, particle size) on solubility are discussed.
Chapter 5: This chapter comprises a description of the properties of each safety relevant radionuclide and the results of the speciation (Pourbaix diagrams calculated with the different TDB’s) and solubility calculations. The concentration limits for all radioelements and potential solubility limiting solids that were calculated with MOLDATA are presented. The reasoning that was taken as basis for the selection procedure is given and the determination of the Source and Expert ranges illustrated. Possible effects of changing/evolving conditions (pH, Eh, pCO2, ionic strength) on the stability of the selected solubility limiting solids are given, if known from literature and/or experimental results. In case experimental data were used to delimit, i.e. enlarge or restrict the SR, they are described. At the end of the element descriptions, general as well as specific remarks concerning the different thermodynamic source data(bases) and/or MOLDATA can be found.
Appendices: Annex I comprises a summarizing table of the concentration limits calculated with MOLDATA for each radionuclide and considered solid, in Annex II the Source and Expert Ranges for each radionuclide are tabulated, in Annex III details on the uncertainty calculations and thermodynamic data are given, and Annex IV represents a Technical Note, in which thorium data are presented and discussed that were published more recently than the ones comprised in the respective NEA review (Vol. 11; Rand et al., 2009) and in MOLDATA.
10
1 Introduction This report was prepared for Research Plan RP. W&D 0064 (WP 3: Radionuclide Speciation and Solubility) and can be considered as an extended update of a previously published report by Wang et al. (2000). It represents also a supporting document of the Level 4 (L4) integration report "Radionuclide migration in the far-field" (Bruggeman et al., in prep.), which corresponds to a deliverable of WP 1 and has to be prepared within the frame of the Safety and Feasibility Case 1 (SFC-1). In the latter report, the processes and mechanisms influencing the radionuclide migration/retention behaviour in the potential host formation and the selection strategy of migration parameters developed by SCK•CEN together with ONDRAF/NIRAS (O/N) are described in detail. Furthermore, a phenomenological model, which has been developed based on the research performed at SCK•CEN during the last 30 years, and which reflects the current state of knowledge with respect to the main BC relevant retention and migration processes, is put forward in this L4 report.
Main objective of the current report was to deliver a reference document, comprising information on the radionuclide speciation and solubility representative for "undisturbed" (far-field) BC conditions and to provide transparent and traceable input parameters, i.e. concentration limits for performance assessment (PA) calculations.
It should be explicitly mentioned that for the speciation and solubility calculations only inorganic ligands and complexation were considered. However, as revealed by long-term migration, batch and solubility experiments, transition metals, actinides and lanthanides are forming strong complexes with the humic substances of BC natural organic matter (NOM). The possible solubility enhancing effect related to the radionuclide association/interaction with the organic matter, which is thought to be of colloidal nature, is thus captured in transport/performance assessment (PA) calculations through the recently developed phenomenological model by Maes et al. (2011). According to this model, two species per radionuclide are allowed to migrate, i.e. the "free radionuclide" and the OM associated radionuclide (Rn-OM complex), which are both represented by their own parameter set.
Special attention has been paid to the description of the methodology put in place by SCK•CEN together with O/N to derive so-called solubility ranges (i.e. Source and Expert ranges) and Best Estimate (BE) values, as well as to the documentation of the selection procedure of the solubility controlling solids.
As such, the present work represents the basis for the numerical solubility values needed as input parameters in the ongoing safety and feasibility study and scoping calculations within the frame of SFC-1.
In total 25 elements are reviewed within this report, 23 safety-relevant radionuclides as well as calcium and carbon.
The radionuclide inventory for both spent fuel and vitrified high-level waste is provided by ONDRAF/NIRAS (https://www.nirond-km.be/gm/folder-1.11.315498, file "DCTs Version 2.2.xls") for a reference time of 10 years after unloading the spent fuel of the reactor. Based on these data, the respective inventory data are calculated taking into account a total cooling time of 60 years for category C waste (Weetjens et al., 2012).
11
2 System definition
2.1 The disposal system The recommended option for the long-term management of high-level and/or long-lived radioactive waste (HLW-LL)1 in Belgium is geological disposal in poorly indurated clay formations (i.e. Ypresian Clays, Boom Clay). Through its Waste Plan, ONDRAF/NIRAS (O/N) intends to satisfy its legal obligations, while providing the Government with all the elements needed to make a "decision in principle", in other words a general policy decision, about the long- term management of HLW-LL (ONDRAF/NIRAS, 2011). Within this framework, Boom Clay (BC) is investigated as one of the potential host formations. The Belgian concept (i.e. Supercontainer concept) for HLW-LL considers the installation of a multi-barrier repository system, which typically comprises the natural barrier, provided by the host rock and its surroundings (aquifers, biosphere), as well as the engineered barrier system (EBS). The former is also referenced as the geosphere or "far-field". The "engineered barrier system" in contrast, represents the man-made or artificial parts placed wihin a repository. The supercontainer design (Wickham et al., 2005) comprises different components, such as the vitrified HLW canisters or SF assemblies, surrounded by a carbon steel overpack, a concrete buffer and a stainless steel liner. The supercontainers are foreseen to be emplaced in galleries excavated in the host rock. Due to the plasticity of poorly indurated clays, the disposal galleries will be stabilized by a concrete lining. The space between the supercontainers and the concrete liner will be backfilled before the galleries will be sealed. The "near-field" includes the EBS and those parts of the host rock in contact with or near the EBS, whose properties have been affected by the presence of the repository.
2.2 Boom Clay reference pore water The solubility (and sorption) of the radionuclides depend on their speciation in the aqueous phase and thus on the chemistry of the porewater. In Table 1, the Boom Clay porewater composition for T=25°C is summarized. This composition was derived from the "Reference Boom Clay composition" given at T=16°C by De Craen et al. (2004). The major ion concentrations as well as the pH and Eh were calculated from/calibrated against an average (44 measurements) MORPHEUS water composition. This approach involved assuming equilibrium with the minerals calcite, pyrite, siderite, chalcedony and kaolinite, taking into account an ion exchange complex of 0.925 eq/5 kg of clay (i.e. 18.5 meq/100 g clay) and a pCO2 of 10-2.62. The same approach was applied to derive the Boom Clay porewater composition at 25°C. Since the Boom Clay reference porewater is highly supersaturated with respect to various clay minerals, mica, quartz, K-feldspar etc., these minerals were not allowed to precipitate during the simulations, in order to keep the water composition unchanged. This approach allowed radionuclide solubilities to be calculated under the reference Boom Clay conditions. It should be mentioned that for the derivation of the BC porewater composition at 25°C, the LLNL TDB (thermo.com.v.8.r6+.dat) was used.
1 Category A (LILW-SL): low- and intermediate-level short-lived waste. This category comprises radio-elements with half-lives of less than 30 years, emitting generally alpha radiation. Category B (LILW-LL): low- and intermediate-level long-lived waste. This category comprises radio-elements with half-lives of more than 30 years, emitting generally alpha radiation. Category C (HLW-SL and HLW-LL):high-level short lived and long-lived waste. This category comprises large amounts of beta and alpha emitting radio-elements having short or long half-lives. They are highly heat-generating. This waste arises from the reprocessing of irradiated nuclear fuel. Spent fuel that is not reprocessed also belongs to this category.
12
Table 1: The (recalculated) Boom Clay water composition at 25°C Species [mol/L] Al 4.63 ×10-8 -4 SiO2(aq) 1.95 ×10 Mg 6.56 × 10-5 Ca 5.04 × 10-5 -6 Fe 3.68 × 10 K 1.85 × 10-4 Na 1.55 × 10-2 Cl- 7.34 × 10-4 2- -5 SO4 2.41 × 10 -2.44 pCO2(g), atm 10
Eh, mV -281 pH 8.36
2.3 Computer code and databases The derivation of the pore water composition, as well as the speciation and solubility calculations presented within this report were performed using the recent releases of The Geochemist's Workbench code, versions 8.08, 8.10 and 8.12 (Bethke, 2010). Geochemist's Workbench represents a set of interactive software tools (ACT2, REACT, TACT and RXN) for solving problems in aqueous geochemistry, including those encountered in environmental protection and remediation, petroleum industry, and economic geology. Two of the four software tools, ACT2 and REACT, were used for the geochemical calculations. The former was applied for the construction of the so-called Pourbaix diagrams, while the latter program was used to derive the Boom Clay reference porewater composition, and to determine the prevalent speciation under Boom Clay conditions, as well as to calculate the radionuclide solubilities in this water.
With respect to the databases used for the different purposes, the following can be stated. The Pourbaix diagrams were calculated for each element with the four individual databases, i.e. the NAGRA/PSI, ANDRA, LLNL and NEA databases, as well with MOLDATA. Aim of this comparative study was to identify differences/discrepancies between the individual databases and MOLDATA, but also to "discover" possible inconsistencies and/or errors within MOLDATA.
The solubility calculations were also performed with all databases, however only the results of the calculations performed with MOLDATA will be presented in more detail, since the latter is considered to be the reference database for any SFC-1 related calculations.
2.4 MOLDATA The so-called MOLDATA project was launched in 2004 as part of the Research Plan "Geochemical and reactive transport modelling" (RP.WD.039) and has been followed-up within the frame of RP.WD.0064 "Radionuclide migration and retention processes in Boom Clay". The main objective has been the compilation of a high-quality and internally consistent thermochemical database that is to be used as reference database for geochemical and reactive transport calculations, which are increasingly applied to interpret phenomenological processes and safety assessment studies. In that sense, MOLDATA represents one of the building blocks on
13
which the PA of the radioactive waste disposal system is based. Therefore, reliability of the modeling results is essential and the quality of the thermodynamic data comprising the database has to be ensured. In 2006, ONDRAF/NIRAS defined quality measures (ONDRAF/NIRAS, 2006) in order to meet the high quality requirements of the Safety Case 1 (SFC-1) in general, and also to ensure the high-quality implementation and use of databases. Database implementation includes a verification process, meaning that the data need to be verified for their scientific correctness and need also to be well documented. Besides this, a high quality database is expected to be state-of- the-art, internally consistent and complete. The compilation of MOLDATA has been performed based on these quality measures. In order to strive for completeness, thermodynamic data from different so-called "source databases" were incorporated into MOLDATA. Thermodynamic data of the following databases have been evaluated within the frame of the MOLDATA compilation/selection procedure:
1) Lawrence Livermore National Laboratories (LLNL) database: thermo.com.v8.r6+.dat and thermo.com.r7beta.dat (Johnson et al., 1991; Delany and Lundeen, 1990; Helgeson et al., 1978); 2) NAGRA/PSI TDB: Chemical Thermodynamic Database 01/01, (Hummel et al., 2002); 3) ANDRA TDB: ThermoChimie v.5, June 2005; 4) NEA TDB’s: Guillaumont et al. (2003): Update on the Chemical Thermodynamics of Uranium, Neptunium, Plutonium, Americium and Technecium (Vol. 5); Gamsjäger et al. (2005): Chemical Thermodynamics of Nickel (Vol. 6); Olin et al. (2005): Chemical Thermodynamics of Selenium (Vol. 7); Brown et al. (2005): Chemical Thermodynamics of Zirconium (Vol. 8); and Rand et al. (2008): Chemical Thermodynamics of Thorium (Vol. 11).
For further details concerning the source databases and the MOLDATA compilation strategy, readers are referred to the report of Wang et al. (2011).
It should be clearly mentioned, that the MOLDATA compilation and selection procedure did not involve a full review of the thermodynamic data. The strategy of the selection procedure has been mainly based on incorporating the most accurate, state-of-the-art and accepted data by the scientific community. For example, with respect to radionuclides this approach means that priority was given to the NEA reviewed data. Other used criteria were the completeness and consistency within the database. When inconsistencies were encountered, the underlying thermodynamic data were reviewed in order to trace the problem and to make adaptations.
MOLDATA currently comprises in total 3855 species, with 84 basis species, 2140 aqueous species, 1436 minerals and 195 gases.
2.5 GDP The in-house developed Geochemical Data Processor (GDP) represents the tool, i.e. computer software (De Soete, 2010) that made the compilation of MOLDATA possible. As mentioned above, four source databases (NAPSI/NAGRA, ANDRA, NEA – Vol. 5-11, and LLNL), comprising a total of ~7500 species were used to compile MOLDATA. In order to process this huge dataset, a tool was needed enabling to upload and store them in a standardized format. The storage of the source data occurs in the so-called "Central Database" which can be seen as the
14
"core" of the GDP. The import of the data is generally done via a wizard like tool, by which the user - in an interactive manner - is guided through the "data upload process". As the source databases exist in different formats, i.e. GWB, PhreeqC or tabdelimited format, the data import into the Central Database requires a standardization process. This conversion of a datasets to a standardized format is done by so-called "import filters". In the same way, data can be exported/downloaded in different formats, for which so-called "export filters" were programmed. Each different format requires a specific import and export filter to be developed. For further details concerning the software, it is referred to the GDP manual (De Soete, 2010).
15
3 Speciation calculations and Pourbaix diagrams
Definition: Speciation is the distribution of aqueous species among free ions, ion pairs, and complexes (Nordstrom and Munoz, 1994).
Based on the water/rock equilibrium model, the speciation of elements (i.e. major elements, transition metals, actinides, lanthanides) can be evaluated by the use of so-called Pourbaix diagrams. A Pourbaix diagram generally illustrates the equilibrium aqueous species distribution of an element as a function of Eh and pH. Additionally, it is possible to map out fields of Eh and pH over which different solid phases are possibly stable. The vertical axis represents the Eh, which is the voltage potential with respect to the standard hydrogen electrode (SHE), as calculated by the Nernst equation:
0 0.059 aOx Eh E log n aRed where: E0 is the standard potential; n is the number of electrons involved in the reaction;
aOx is the activity of oxidized species; and
aRed is the activity of reduced species.
The horizontal axis represents the pH, with
pH = -log aH+
and aH+: activity of hydrogen ions.
The boundaries of fields within which a particular aqueous species dominates the aqueous speciation of an element, and stability fields of minerals, are represented/delimited by lines. The diagrams that are presented within this report have been calculated using the reference porewater composition for the Boom Clay (see paragraph 2.2) and radionuclide activities of 10-8. The boundary between a field representing the predominance of a particular aqueous species and a solid phase represents equilibrium for this radionuclide concentration.
It has to be mentioned, that the activity of the dissolved ions, the presence of binding agents (ligands), as well as temperature may modify a Pourbaix diagram and shift the lines/boundaries in accordance with the Nernst equation. The in-situ Boom Clay conditions, i.e. Eh -281 mV and pH ~8.4 are represented by a cross within the Pourbaix diagrams and enable the immediate recognition of the dominant species or stable solid phase.
In calculating Pourbaix diagrams, the effect of complexing anions on radionuclide speciation has to be treated with care, since some anions have varying speciation over the Eh-pH range considered. The speciation of dissolved sulphur and inorganic carbon is redox- and pH-dependent (see Figure 1). Consequently, erroneous results would be obtained by using constant 2- - concentrations for the complexing ligands SO4 and HCO3 (for example) over the entire Eh-pH
16
region. In the present report, the speciation of sulphur and inorganic carbon are taken into account by using so-called mosaic diagrams (Bethke, 2010). The chloride (Cl-) speciation, in contrast, is constant over the whole Eh-pH range.
1
HSO - 2- 4 SO4
.5 CO2(aq)
- HCO3
Eh (volts) Eh 0 2- H2S(aq) CO3
Methane µ
–.5 HS- 25 C 0 2 4 6 8 10 12 14 pH
Figure 1: Sulphur and carbon speciation as function of Eh and pH
Complexation of the radionuclides by organic ligands was not taken into account in the current calculations.
3.1 Species activities and activity coefficients In order to understand the qualitative meaning of activity coefficients, it is important to consider how solution concentration affects ion interaction. Generally, ions in solution interact with each other as well as with the water molecules and neutrally charged solute species. At low concentrations of a species i (ci), ionic interactions can be ignored. The properties of such a solution approach those of an ideal solution, and would become ideal in an infinitely diluted solution, so that:
[i] = ci
Determining equilibrium constants experimentally is commonly done by using media of high ionic strength (e.g. solutions of NaClO4, KNO3, Na2SO4). These electrolytes represent so-called "non- ideal solutions". In order to account for their deviation from ideality, the concept of activity coefficients, γ, was introduced into thermodynamics, according to which:
[i] = i × ci
Under dilute conditions, γ is expected to be close to 1 (i.e. approaching an ideal solution).
According to Coulomb's law, the electrostatic force between two point electric charges is proportional to the product of the charges and inversely proportional to the square of the distance between the two charges. Therefore, activity coefficients in dilute solutions decrease with increasing concentration, due to the fact that the coulombic forces become stronger as ions pack more closely together.
17
Different methods exist for the calculation of activity coefficient in electrolyte solutions. The most commonly used existing activity models are all based on the Debye-Hückel theory. This theory was proposed by Peter Debye and Erich Hückel in 1923 to calculate mean activity coefficients for ions in dilute solutions, and was further elaborated by Robinson and Stokes (1968), who derived the well-known Debye-Hückel equation:
A z2 I log i i 0 1 ai B I where 0 ai (given in Å): is the activity coefficient; zi is the electrical charge; A, B are constants (T-dependent; at 25°C A = 0.5092 and B = 0.3283); and I is the ionic strength of the solution (given in molal units).
The ionic strength I is defined as half the sum of the product of each species (molality mi) and the square of its charge: 1 I m z 2 i i
A disadvantage/limitation of the Debye-Hückel equation is that it becomes inaccurate at moderate ionic strength, i.e. above about 0.1 molal (Stumm and Morgan, 1996).
Another model, that can be carried to somewhat higher ionic strengths, i.e. 0.3-0.5, is the so-called Davies model. The Davies equation can be considered as a variant of the Debye-Hückel equation and is represented by:
2 I log i A zi 0.3 I 1 I
0 The equation was simplified from the original Debye-Hückel equation by noting that ai B is about 1 at 25°C.
An alternative correction term similar was also added to the Debye-Hückel equation by Helgeson (1969), to produce the so-called ‘Bdot equation’:
2 Azi I log i 0 B I 1 ai B I where: B is a coefficient that depends on the charges of the dissolved species and temperature.
This ‘Bdot’ equation can yield accurate estimates of i for ionic strengths of up to about 1 molal in NaCl-dominated solutions.
18
4 Solubility - Definition and theoretical background Definition for solubility: The solubility of a given element is the sum of the stoichiometric concentrations of all dissolved species containing the element (Garrels and Christ, 1965).
4.1 Solubility constants Dissolution and precipitation reactions of crystals/minerals represent heterogeneous equilibrium reactions, meaning they are involving solid and liquid/aqueous species, i.e. different phases. Conventionally, equilibrium constants involving a solid compound are typically denoted as "solubility constants", rather than as formation constants of the solid (Brown et al., 2005). Solubility constants can be defined as follows: