Publications

Abstract not provided.

Numerical simulation of a repository for heatgenerating nuclear waste in fractured crystalline rock requires a method for simulating coupled heat and fluid flow and reactive radionuclide transport in both porous media (bentonite buffer, surface sediments) and fractured rock (the repository host rock). Discrete fracture networks (DFNs), networks of two-dimensional planes distributed in a three-dimensional domain, are commonly used to simulate isothermal fluid flow and particle transport in fractures, but unless coupled to a continuum, are incapable of simulating heat conduction through the rock matrix, and therefore incapable of capturing the effects of thermally driven fluid fluxes or of coupling chemical processes to thermal processes. We present a method for mapping a stochastically generated DFN to a porous medium domain that allows representation of porous and fractured media in the same domain, captures the behavior of radionuclide transport in fractured rock, and allows simulation of coupled heat and fluid flow including heat conduction through the matrix of the fractured rock. We apply the method within Sandia's Geologic Disposal Safety Assessment (GDSA) framework to conduct a post-closure performance assessment (PA) of a generic repository for commercial spent nuclear fuel in crystalline rock. The three-dimensional, kilometer-scale model domain contains approximately 4.5 million grid cells; grid refinement captures the detail of 3, 360 individual waste packages in 42 disposal drifts. Coupled heat and fluid flow and reactive transport are solved numerically with PFLOTRAN, a massively parallel multiphase flow and reactive transport code. Simulations of multiple fracture realizations were run to 1 million years, and indicate that, because of the channeled nature of fracture flow, thermally-driven fluid fluxes associated with peak repository temperatures may be a primary means of radionuclide transport out of the saturated repository. The channeled nature of fracture flow gives rise to unique challenges in uncertainty and sensitivity quantification, as radionuclide concentrations at any given location outside the repository depend heavily on the distribution of fractures in the domain.

Abstract not provided.

The Deep Borehole Field Test (DBFT) is a planned multi-year project led by the US Department of Energy's Office of Nuclear Energy to drill two boreholes to 5 km total depth into crystalline basement in the continental US. The purpose of the first characterization borehole is to demonstrate the ability to characterize in situ formation fluids through sampling and perform downhole hydraulic testing to demonstrate groundwater from 3 to 5 km depth is old and isolated from the atmosphere. The purpose of the second larger-diameter borehole is to demonstrate safe surface and downhole handling procedures. This paper details many of the drilling, testing, and characterization activities planned in the first smaller-diameter characterization borehole.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Deep geological disposal of nuclear waste in clay/shale/argillaceous rock formations has received much consideration given its desirable attributes such as isolation properties (low permeability), geochemically reduced conditions, slow diffusion, sorbtive mineralogy, and geologically widespread (Jové Colón et al., 2014). There is a wealth of gained scientific expertise on the behavior of clay/shale/ argillaceous rock given its focus in international nuclear waste repository programs that includes underground research laboratories (URLs) in Switzerland, France, Belgium, and Japan. Jové Colón et al. (2014) have described some of these investigative efforts in clay rock ranging from site characterization to research on the engineered barrier system (EBS). Evaluations of disposal options that include nuclear waste disposition in clay/shale/argillaceous rock have determined that this host media can accommodate a wide range of waste types. R&D work within the Used Fuel Disposition Campaign (UFDC) assessing thermal effects and fluid-mineral interactions for the disposition of heat-generating waste have so far demonstrated the feasibility for the EBS and clay host rock to withstand high thermal loads. This report represents the continuation of disposal R&D efforts on the advancement and refinement of coupled Thermal-Hydrological-Mechanical-Chemical (THMC), hydrothermal experiments on clay interactions, used fuel degradation (source term), and thermodynamic modeling and database development. The development and implementation of a clay/shale/argillite reference case described in Jové Colón et al. (2014) for FY15 will be documented in another report (Mariner et al. 2015) – only a brief description will be given here. This clay reference case implementation is the result of integration efforts between the GDSA PA and disposal in argillite work packages. The assessment of sacrificial zones in the EBS is being addressed through experimental work along with 1D reactive-transport and reaction path modeling. The focus of these investigations into the nature of sacrificial zones is to evaluate the chemical effects of heterogeneous chemical reactions at EBS interfaces. The difference in barrier material types and the extent of chemical reactions within these interfacial domains generates changes in mineral abundances. These mineralogical alterations also result in volume changes that, although small, could affect the interface bulk porosity. As in previous deliverables, this report is structured according to various national laboratory contributions describing R&D activities applicable to clay/shale/argillite media.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as for vapor-phase contaminant transport and/or containment of noble gases in the subsurface. In this second progress report of FY22, we present two ongoing activities related to imbibition testing on volcanic rock samples. We present the development of a new analytical solution predicting the temperature response observed during imbibition into dry samples, as discussed in our previous first progress report for FY22. We also illustrate the use of a multi-modal capillary pressure distribution to simulate both early- and late-time imbibition data collected on tuff core that can exhibit multiple pore types. These FY22 imbibition tests were conducted for an extended period (i.e., far beyond the time required for the wetting front to reach the top of the sample), which is necessary for parameter estimation and characterization of two different pore types within the samples.

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as containment transport and/or containment of gases in the subsurface. To estimate rock fluid flow properties and subsequently predict physically realistic processes such as patterns and timing of water, gas, and energy (e.g., heat) movement in the subsurface, laboratory spontaneous water imbibition with simultaneous temperature measurement and numerical modeling methods are presented in the FY22 progress report. A multiple-overlapping-continua conceptual model is used to explain and predict observed complex multi-phenomenological laboratory test behavior during spontaneous imbibition experiments. This report primarily addresses two complexities that arise during the experiments: 1) capturing the late-time behavior of spontaneous imbibition tests with dual porosity; and 2) understanding the thermal perturbation observed at or ahead of the imbibing wetting front, which are associated with adsorption of water in initially dry samples. We use numerical approaches to explore some of these issues, but also lay out a plan for further laboratory experimentation and modeling to best understand and leverage these unique observations.

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

Abstract not provided.

This report describes the development of a comprehensive catalogue of generic features, events, and processes (FEPs) that are potentially important for the post-closure performance of a repository for high-level radioactive waste (HLW) and spent nuclear fuel (SNF) in salt (halite) host rock. The FEPs and other supporting information have been entered into a “SaltFEP” Database. The generic salt repository FEPs include consideration of relevant FEPs from a number of U.S., Dutch, German, and international FEP lists and should be a suitable starting point for any repository program in salt host rock. The salt FEP catalogue and database employ a FEP classification matrix approach that is based on the concept that a FEP is typically a process or event acting upon or within a feature. The FEP matrix provides a two-dimensional structure consisting of a Features/Components axis that defines the “rows” and a Processes/Events axis that defines the “columns” of the matrix. The design of the FEP classification matrix is consistent with repository performance assessment – the Features/Components axis is organized vertically to generally correspond to the direction of potential radionuclide migration (from the waste to the biosphere) and the Processes/Events axis is designed to represent the common two-way couplings between thermal processes and other processes (such as thermal-mechanical or thermal-hydrologic processes). Related FEPs can be easily identified – related FEPs will typically be grouped in a single matrix cell or aligned along a common row (Feature/Component) or column (Process/Event). The online SaltFEP database can be downloaded from www.saltfep.org. It contains the FEP matrix, the FEPs, and the associated processes for each FEP. It provides a starting point to create and document site-specific individual FEPs. Furthermore, the FEP matrix is connected to the Salt Knowledge Archive (SKA), a database of about 20,000 references and documents representing the historical knowledge on radioactive disposal in salt. This work is the result of an ongoing collaboration between researchers in the U.S., the Netherlands, and Germany, and supports the NEA Salt Club Mandate. It builds upon prior work which is documented.

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.