Publications

Abstract not provided.

Abstract not provided.

This report summarizes the 2020 fiscal year (FY20) status of the borehole heater test in salt funded by the US Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) campaign. This report satisfies SFWST level-two milestone number M2SF-20SNO10303032. This report is an update of an August 2019 level-three milestone report to present the final as-built description of the test and the first phase of operational data (BATS la, January to March 2020) from the Brine Availability Test in Salt (BATS) field test.

In this study, I present experimental results on the equilibrium between boracite [Mg3B7O13Cl(cr)] and kurnakovite [chemical formula, Mg2B6O11.15H2O(cr); structural formula, MgB3O3(OH)5.5H2O(cr)] at 22.5 ± 0.5 °C from a long-term experiment up to 1629 days, approaching equilibrium from the direction of supersaturation, Mg3B7O13Cl(cr) + H+ + 2B(OH)4 + 18H2O(1) . 3MgB3O3(OH)5.5H2O(cr) + Cl . Based on solubility measurements, the 10-based logarithm of the equilibrium constant for the above reaction at 25 °C is determined to be 12.83 ± 0.08 (2s). Based on the equilibrium constant for dissolution of boracite, Mg3B7O13Cl(cr) + 15H2O(l) = 3Mg2+ + 7B(OH)4 + Cl + 2H+ at 25 °C measured previously (Xiong et al. 2018) and that for the reaction between boracite and kurnakovite determined here, the equilibrium constant for dissolution of kurnakovite, MgB3O3(OH)5.5H2O(cr) = Mg2+ + 3B(OH)4 + H+ + H2O(1) is derived as 14.11 ± 0.40 (2s). Using the equilibrium constant for dissolution of kurnakovite obtained in this study and the experimental enthalpy of formation for kurnakovite from the literature, a set of thermodynamic properties for kurnakovite at 25 °C and 1 bar is recommended as follows: ΔH0f = 4813.24 ± 4.92 kJ/mol, .G0f = 4232.0 ± 2.3 kJ/mol, and S0 = 414.3 ± 0.9 J/(mol.K). Among them, the Gibbs energy of formation is based on the equilibrium constant for kurnakovite determined in this study; the enthalpy of formation is from the literature (Li et al. 1997), and the standard entropy is calculated internally with the Gibbs-Helmholtz equation in this work. The thermodynamic properties of kurnakovite estimated using the group contribution method for borate minerals based on the sums of contributions from the cations, borate polyanions, and structural water to the thermodynamic properties from the literature (Li et al. 2000) are consistent, within their uncertainties, with the values listed above. Since kurnakovite usually forms in salt lakes rich in sulfate, studying the interactions of borate with sulfate is important to modeling kurnakovite in salt lakes. For this purpose, I have re-calibrated our previous model (Xiong et al. 2013) describing the interactions of borate with sulfate based on the new solubility data for borax in Na2SO4 solutions presented here.

The US Department of Energy Office of Nuclear Energy is conducting a brine availability heater test to characterize the thermal, mechanical, hydrological and chemical response of salt at elevated temperatures. In the heater test, brines will be collected and analyzed for chemical compositions. In order to support the geochemical modeling of chemical evolutions of the brines during the heater test, we are recalibrating and validating the solubility models for the mineral constituents in salt formations up to 100°C, based on the solubility data in multiple component systems as well as simple systems from literature. In this work, we systematically compare the model-predicted values based on the various solubility models related to the constituents of salt formations, with the experimental data. As halite is the dominant constituent in salt formations, we first test the halite solubility model in the Na-Mg-Cl dominated brines. We find the existing halite solubility model systematically over-predict the solubility of halite. We recalibrate the halite model, which can reproduce halite solubilities in Na-Mg-Cl dominated brines well. As gypsum/anhydrite in salt formations controls the sulfate concentrations in associated brines, we test the gypsum solubility model in NaCl solutions up to 5.87 mol•kg-1 from 25°C to 50°C. The testing shows that the current gypsum solubility model reproduces the experimental data well when NaCl concentrations are less than 1 mol•kg-1. However, at NaCl concentrations higher than 1, the model systematically overpredicts the solubility of gypsum. In the Na - Cl - SO4 - CO3 system, the validation tests up to 100°C demonstrate that the model excellently reproduces the experimental data for the solution compositions equilibrated with one single phase such as halite (NaCl) or thenardite (Na2SO4), with deviations equal to, or less than, 1.5 %. The model is much less ideal in reproducing the compositions in equilibrium with the assemblages of halite + thenardite, and of halite + thermonatrite (Na2CO3•H2O), with deviations up to 31 %. The high deviations from the experimental data for the multiple assemblages in this system at elevated temperatures may be attributed to the facts that the database has the Pitzer interaction parameters for Cl - CO3 and SO4 - CO3 only at 25°C. In the Na - Ca - SO4 - HCO3 system, the validation tests also demonstrate that the model reproduces the equilibrium compositions for one single phase such as gypsum better than the assemblages of more than one phase.

Abstract not provided.

Abstract not provided.

This report summarizes the 2019 fiscal year (FY19) status of the borehole heater test in salt funded by the US Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) campaign. This report satisfies SFWST level-three milestone report M3SF-19SN010303033. This report is an update of the April 2019 level-two milestone report M2SF-19SNO10303031 to reflect the nearly complete as-built status of the borehole heater test. This report discusses the fiscal year 2019 (FY19) design, implementation, and preliminary data interpretation plan for a set of borehole heater tests call the brine availability tests in salt (BATS), which is funded by the DOE Office of Nuclear Energy (DOE-NE) at the Waste Isolation Pilot Plant (WIPP), a DOE Office of Environmental Management (DOE-EM) site. The organization of BATS is outlined in Project Plan: Salt In-Situ Heater Test (SNL, 2018). An early design of the field test is laid out in Kuhlman et al. (2017), including extensive references to previous field tests, which illustrates aspects of the present test. The previous test plan by Stauffer et al. (2015) places BATS in the context of a multi-year testing strategy, which involves tests of multiple scales and processes, eventually culminating in a drift-scale disposal demonstration. This level-3 milestone report is an update of a level-2 milestone report from April 2019 by the same name. The update adds as-built details of the heater test, which at the time of writing (August 2019) is near complete implementation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report discusses the fiscal year 2019 (FY19) design, implementation, and preliminary data interpretation plan for a set of borehole heater tests call the brine availability tests in salt (BATS), which is funded by the DOE Office of Nuclear Energy (DOE-NE) at the Waste Isolation Pilot Plant (WIPP). The organization of BATS is outlined in Project Plan: Salt In-Situ Heater Test. An early design of the field test is laid out in Kuhlman et al., including extensive references to previous field tests, which illustrates aspects of the present test. The previous test plan by Stauffer et al., places BATS in the context of a multi-year testing strategy, which involves tests of multiple scales and processes, possibly culminating in a drift-scale disposal demonstration.

Abstract not provided.

Uranyl ion, UO22+, and its aqueous complexes with organic and inorganic ligands, are the dominant species for transport of natural occurring uranium at the Earth surface environments. In the nuclear waste management, uranyl ion and its aqueous complexes are expected to be responsible for uranium mobilization in the disposal concepts where spent fuel is disposed in oxidized environments such as unsaturated zones relative to the underground water table. In the natural environments, oxalate, in fully deprotonated form, C2O42-, is ubiquitous, as oxalate is one of the most important degradation products of humic and fulvic acids. Oxalate is known to form aqueous complexes with uranyl ion to facilitate the transport of uranium. However, oxalate also forms solid phases with uranyl ion in certain environments, limiting the movement of uranium. Therefore, the knowledge of the stability constants of aqueous and solid uranyl oxalate complexes is important not only to the understanding of the mobility of uranium in natural environments, but also to the performance assessment of radionuclides in geological repositories for spent nuclear fuel. In this work, we present the stability constants for UO2C2O4(aq) and UO2(C2O4)22- at infinite dilution based on our evaluation of the literature data over a wide range of ionic strengths up to 9.5 mol•kg-1. We also obtain the solubility constants at infinite dilution for the following solid uranyl oxalates, UO2C2O4•3H2O and UO2C2O4•H2O, based on the solubility data in a wide range of ionic strengths up to 11 mol•kg-1. In our evaluation, we use the computer code EQ3/6 Version 8.0a. The model developed by us is expected to enable researchers to accurately assess the role of oxalate in mobilization/immobilization of uranium under various conditions including those in geological repositories.

Montmorillonite with an empirical formula of Na0.2Ca0.1Al2Si4O10(OH)2(H2O)10 is a di-octahedral smectite. Montmorillonite-rich bentonite is a primary buffer candidate for high level nuclear waste (HLW) and used nuclear fuel to be disposed in mild environments. In such environments, temperatures are expected to be ≤ 90oC, the solutions are of low ionic strengths, and pH is close to neutral. Under the conditions outside the above parameters, the performance of montmorillonite-rich bentonite is deteriorated because of collapse of swelling particles as a result of illitization, and dissolution of the swelling clay minerals followed by precipitation of non-swelling minerals. It has been well known that tri-octahedral smectites such as saponite, with an ideal formula of Mg3(Si, Al)4O10(OH)2•4H2O for an Mg-end member (saponite-15A), are less susceptible to alteration under harsh conditions. Recently, Mg-bearing saponite has been favorably considered as a preferable engineered buffer material for the Swedish very deep holes (VDH) disposal concept in crystalline rock formations. In the VDH, HLW is disposed in deep holes at depth between 2,000 m and 4,000 m. At such deployment depths, the temperatures are expected to be between 100oC and 150oC, and the groundwater is of high ionic strength. The harsh chemical conditions of high pH are also introduced by the repository designs in which concretes and cements are used as plugs and buffers. In addition, harsh chemical conditions introduced by high ionic strength solutions are also present in repository designs in salt formations and sedimentary basins. For instance, the two brines associated with the salt formations for the Waste Isolation Pilot Plant (WIPP) in USA have ionic strengths of 5.82 mol•kg-1 (ERDA-6) and 8.26 mol•kg-1 (GWB). In the Asse site proposed for a geological repository in salt formations in Germany, the Q-brine has an ionic strength of ~13 mol•kg-1. In this work, we present our investigations regarding the stability of saponite under hydrothermal conditions in harsh environments.

Abstract not provided.

Abstract not provided.

In the published article (Xiong et al., 2017), there was an error for the reaction coefficient for the dissolution reaction of Ca2C10H12N2O8•7H2O(s) in the database used for the modeling. In the database for the modeling, the coefficient for water (i.e., 7H2O) was inadvertently omitted. Because of this omission, the results in Table 3 were affected. The authors wish to make the corrections to Table 3. The corrected values are tabulated in the revised Table 3. The corrected values reproduce the experimental data in MgCl2 solutions much better (see revised Fig. 5).

Abstract not provided.

Methane (CH4) and carbon dioxide (CO2), the two major components generated from kerogen maturation, are stored dominantly in nanometer-sized pores in shale matrix as (1) a compressed gas, (2) an adsorbed surface species and/or (3) a species dissolved in pore water (H2O). In addition, supercritical CO2 has been proposed as a fracturing fluid for simultaneous enhanced oil/gas recovery (EOR) and carbon sequestration. A mechanistic understanding of CH4-CO2-H2O interactions in shale nanopores is critical for designing effective operational processes. Using molecular simulations, we show that kerogen preferentially retains CO2 over CH4 and that the majority of CO2 either generated during kerogen maturation or injected in EOR will remain trapped in the kerogen matrix. The trapped CO2 may be released only if the reservoir pressure drops below the supercritical CO2 pressure. When water is present in the kerogen matrix, it may block CH4 release. However, the addition of CO2 may enhance CH4 release because CO2 can diffuse through water and exchange for adsorbed methane in the kerogen nanopores.

In this study, solubility measurements regarding boracite [Mg3B7O13Cl(cr)] and aksaite [MgB6O7(OH)6·2H2O(cr)] from the direction of supersaturation were conducted at 22.5 ± 0.5 °C. The equilibrium constant (log10K0) for boracite in terms of the following reaction, Mg3B7O13Clcr+15H2Ol⇌3Mg2++7BOH4 −+Cl−+2H+ is determined as −29.49 ± 0.39 (2σ) in this study. The equilibrium constant for aksaite according to the following reaction, MgB6O7OH6·2H2Ocr+9H2Ol⇌Mg2++6BOH4 −+4H+ is determined as −44.41 ± 0.41 (2σ) in this work. This work recommends a set of thermodynamic properties for aksaite at 25 °C and 1 bar as follows: ΔHf 0 = −6063.70 ± 4.85 kJ·mol−1, ΔGf 0 = −5492.55 ± 2.32 kJ·mol−1, and S0 = 344.62 ± 1.85 J·mol−1·K−1. Among them, ΔGf 0 is derived from the equilibrium constant for aksaite determined by this study; ΔHf 0 is from the literature, determined by calorimetry; and S0 is computed in the present work from ΔGf 0 and ΔHf 0. This investigation also recommends a set of thermodynamic properties for boracite at 25 °C and 1 bar as follows: ΔHf 0 = −6575.02 ± 2.25 kJ·mol−1, ΔGf 0 = −6178.35 ± 2.25 kJ·mol−1, and S0 = 253.6 ± 0.5 J·mol−1·K−1. Among them, ΔGf 0 is derived from the equilibrium constant for boracite determined by this study; S0 is from the literature, determined by calorimetry; and ΔHf 0 is computed in this work from ΔGf 0 and S0. The thermodynamic properties determined in this study can find applications in many fields. For instance, in the field of material science, boracite has many useful properties including ferroelectric and ferroelastic properties. The equilibrium constant of boracite determined in this work will provide guidance for economic synthesis of boracite in an aqueous medium. Similarly, in the field of nuclear waste management, iodide boracite [Mg3B7O13I(cr)] is proposed as a waste form for radioactive 129I. Therefore, the solubility constant for chloride boracite [Mg3B7O13Cl(cr)] will provide the guidance for the performance of iodide boracite in geological repositories. Boracite/aksaite themselves in geological repositories in salt formations may be solubility-controlling phase(s) for borate. Consequently, solubility constants of boracite and aksaite will enable researchers to predict borate concentrations in equilibrium with boracite/aksaite in salt formations.

In this work, a solubility study on brucite [Mg(OH)2(cr)] in Na2SO4 solutions ranging from 0.01 to 1.8 mol·kg−1, with 0.001 mol·kg−1 borate, has been conducted at 22.5 °C. Based on the solubility data, the Pitzer interaction parameters for MgB(OH)4 + − SO4 2− and MgB(OH)4 + − Na+ along with the formation constant for MgSO4(aq) are evaluated using the Pitzer model. The formation constant (log10β10 = 2.38 ± 0.08) for MgSO4(aq) at 25 °C and infinite dilution obtained in this study is in excellent agreement with the literature values. The experimental data on the solubility of gypsum (CaSO4·2H2O), at 25 °C, in aqueous solutions of MgSO4 with ionic strengths up to ~ 11 mol·kg−1 were analyzed using models with and without considering the MgSO4(aq) species. The model incorporating MgSO4(aq) fits better to the experimental data than the model without MgSO4(aq), especially in the ionic strength range beyond ~ 4 mol·kg−1, demonstrating the need for incorporation of MgSO4(aq) into the model to improve the accuracy.

Abstract not provided.

Abstract not provided.

Abstract not provided.