[en] The deep borehole disposal (DBD) concept for high-activity waste consists of drilling a large-diameter borehole into crystalline basement rock to a depth of about 5,000 m, placing waste packages in the lower portion of the borehole - the waste emplacement zone - and sealing and plugging the upper portion of the borehole with a combination of bentonite, cement plugs, and sand or crushed rock ballast. The waste emplacement zone is several times deeper than typical mined repositories and is well below the typical maximum depth of fresh groundwater resources. Recent work describes a generic reference case for deep borehole disposal of Cs and Sr capsules stored at the Hanford Site. Two scenarios are described: an undisturbed (nominal) scenario, and disturbed scenario, in which a waste package is stuck adjacent to a transmissive fracture zone in crystalline basement. Simulations of the processes affecting radionuclide release and transport in the two scenarios provide the basis for quantitative post-closure performance assessment (PA) of the reference case. This paper presents a probabilistic PA of the stuck waste package scenario, i.e. uncertainties affecting radionuclide mobilization and transport are propagated through 200 realizations of a predictive model to quantify uncertainties in predicted radionuclide concentrations within the borehole and the fracture zone. Simulations predict no radionuclide transport in the borehole to distances greater than 100 m above the stuck waste package. An assumed lateral pressure gradient results in predictions of limited radionuclide transport in the fracture zone at a lateral distance of 200 m from the stuck waste package. These results complement previous results for the nominal scenario and contribute to the development of a generic safety case for DBD. (authors)

[en] Density-dependent and topographically-driven subsurface brine flow was modeled in regional-scale and Deep Borehole Disposal-scale domains under varying conditions. These conditions included variable basement permeability, brine density, and hydraulic head. Findings revealed the possibility of hydraulic isolation of deep basement for at least 1 million years under appropriate conditions, which supports the viability of the Deep Borehole Disposal concept. (authors)

[en] The US Department of Energy (DOE) performed an assessment of disposal options for spent nuclear fuel (SNF) and high-level radioactive waste (HLW) that recommended the consideration of deep borehole disposal (DBD) of smaller DOE-managed waste forms, such as cesium (Cs) and strontium (Sr) capsules. To further assess the safety and viability of the DBD concept, post-closure performance assessment (PA) analyses were performed for DBD of Cs/Sr capsules. The post-closure PA included a reference design concept for Cs/Sr capsule disposal, feature, event, and process (FEP) analysis, scenario development, and modeling using the PFLOTRAN code. PA simulations were run for both nominal and disturbed scenarios for long-term radionuclide transport away from the deep disposal borehole. The nominal (expected) post-closure release scenario included FEPs for short-duration (a few hundred years) thermally-induced upward fluid flux through the borehole seals and/or disturbed rock zone (DRZ) followed by longer-term slow diffusive transport. Simulations included a baseline deterministic run and a set of probabilistic realizations (using Latin Hypercube Sampling (LHS) from parameter distributions) to examine the sensitivity and importance of the long-term radionuclide transport to processes and parameters such as: waste package durability; seal porosity and permeability, disturbed rock zone (DRZ) porosity and permeability; and radionuclide sorption in the seal materials, DRZ, and crystalline host rock. The disturbed scenario included a waste package 'stuck' in the crystalline basement above the emplacement zone near a hypothetical borehole-intersecting fracture. Simulations included two deterministic runs to examine the sensitivity to the regional head gradient. The nominal and disturbed scenario PA results suggest that a favorable safety case can be developed for DBD of Cs/Sr capsules and helped to identify key areas of uncertainty, which in turn can inform future DBD research and development (R and D) activities, including a field demonstration project. (authors)

[en] A performance assessment (PA) is an important component of a comprehensive safety analysis for a nuclear waste repository. In a PA, probabilistic simulations of the total repository system are performed, and results are evaluated against performance metrics. Uncertainty and sensitivity analyses help prioritize additional research and model development. The United States Department of Energy (DOE) has been developing a state of the art PA simulation software tool-kit, Geologic Disposal Safety Assessment (GDSA) Framework, that couples increasingly higher fidelity models of subsystem processes into total system PA simulations. Over the past several years, PAs of generic repositories in three geologic media (salt, shale, and crystalline rock) have demonstrated ongoing developments in capability. The current PA of a nuclear waste repository in a generic shale formation showcases GDSA Framework, including capabilities in domain discretization (Cubit), multi-physics simulations (PFLOTRAN), uncertainty and sensitivity analysis (Dakota), and visualization (Paraview). The generic shale reference case considers the disposal of 22,000 metric tons heavy metal of commercial spent nuclear fuel (SNF) in a generic shale formation. PA simulations account for the thermal load and radionuclide inventory of the waste form, components of the engineered barrier system including waste package and buffer, and components of the natural barrier system including the host rock shale and underlying and overlying stratigraphic units. Two repository layouts are considered, one for emplacement of waste packages containing 12 pressurized water reactor (PWR) assemblies, and one for 4-PWR waste package emplacement. Model domains are half-symmetry and contain between 7 and 22 million grid cells. Grid refinement captures the detail of individual waste packages, emplacement drifts, access drifts, and shafts. Simulations are run in a high-performance computing (HPC) environment on 512 or 2048 processes. The governing equations describing coupled heat and fluid flow and reactive transport are solved with PFLOTRAN, an open-source, massively parallel multiphase flow and reactive transport code. Additional simulated processes include waste package degradation; waste form dissolution; radioactive decay and ingrowth in aqueous, solid, and adsorbed phases; and a simple biosphere based on a pumping well. Simulations are run to 106 y; dose and radionuclide concentrations are observed within aquifers at a point approximately 5 km downgradient of the repository. Dakota is used to sample likely ranges of input parameters including waste form and waste package degradation rates and properties of engineered and natural materials to quantify uncertainty in predicted concentrations and sensitivity to input parameters. Given the assumptions of the reference case, results of 12-PWR and 4-PWR simulations are very similar. Dose at the pumping well 5-km downgradient of the repository does not exceed 10-9 Sv/y. I-129 is the largest contributor to dose. The dose due to Cl-36, the next largest contributor, is orders of magnitude smaller. At the pumping well, the concentration of I-129 5-km downgradient does not exceed 10-14 mol/L. Rank correlation coefficients indicate that the maximum concentration of I-129 within the aquifer is sensitive to shale (repository host rock) porosity and aquifer permeability. These results affirm that HPC-capable codes can be used to simulate important multi-physics couplings directly in a total system performance assessment of a geologic repository, and can be used in prioritization of future research and model development. (authors)

[en] Waste Isolation Pilot Plant (WIPP) performance assessment (PA) calculations estimate the probability of radionuclide release from the repository to the land surface and across the land withdrawal boundary for a regulatory period of 10,000 years after facility closure. Simulations of flow and transport in the repository and the surrounding Salado Formation are foundational to the PA. Because proposed additional waste emplacement panels would result in an asymmetric repository layout, the US Department of Energy (DOE) is preparing to transition to use of a three-dimensional (3-D) model domain for simulation of flow and transport instead of the two-dimensional (2-D) flared grid domain currently used. DOE has charged Sandia National Laboratories with developing the capability necessary to simulate processes affecting flow and transport in the WIPP in PFLOTRAN, an open-source massively parallel multi-phase flow and reactive transport code. The new flow and transport capabilities developed in PFLOTRAN incorporate WIPP-specific process models and will replace the 2-D simulators (BRAGFLO and NUTS) that are currently utilized for Salado flow and transport calculations in WIPP PA. The focus of this paper is on the development of a new Nuclear Waste Transport (NWT) mode in PFLOTRAN that has all of the capabilities necessary for Salado transport simulations, including the ability to handle complete dry-out (100% gas saturation) of arbitrary cells in the model domain, radionuclide mass conservation at step changes in porosity associated with borehole intrusion, and the ability to calculate fluxes on a flared grid. The new PFLOTRAN transport capability and a suite of verification tests were designed around a list of functional requirements for WIPP PA calculations. (authors)

[en] The U.S. Department of Energy Office of Nuclear Energy (DOE-NE) performed an assessment of disposal options that recommended the consideration of deep borehole disposal of smaller DOE-managed waste forms, such as cesium (Cs) and strontium (Sr) capsules. To assess the feasibility of the deep borehole disposal concept for Cs/Sr capsules, safety case considerations are identified and examined. A safety case includes quantitative (e.g., safety assessments) and qualitative information related to both pre-closure (operational) and post-closure safety. For deep borehole disposal of Cs/Sr capsules, pre-closure safety considers potential hazards associated with waste package surface handling and downhole emplacement activities; post-closure safety considers scenarios for long-term radionuclide transport to the biosphere. A preliminary list of qualitative indicators of pre-closure and post-closure safety is presented. Research, development, and demonstration (RD and D) activities that can provide a quantitative basis for these indicators are discussed. In particular, DOE-NE has initiated a Deep Borehole Field Test, using surrogate test packages without radioactive waste, to further investigate various aspects of site characterization, deep drilling, and waste package handling and emplacement that would be needed for a full-scale disposal facility for Cs/Sr capsules or other DOE-managed waste forms. Finally, preliminary results from pre-closure hazard analyses and post-closure performance assessment (PA) calculations are presented as examples of supporting information for the quantitative safety case considerations. These preliminary results suggest that a favorable safety case can be developed for deep borehole disposal of Cs/Sr capsules. (authors)

[en] Deep borehole disposal is one of the potential options for safely isolating high-level radioactive waste. It is expected that existing drilling technology can provide reliable and cost-effective construction of suitable deep boreholes. In addition, favorable disposal conditions such as low permeability host rock, high salinity, and geochemically reducing conditions, exist at depth in many locations. Coupled thermal-hydrologic processes induced by heat from the radioactive waste may impact fluid flow and the associated migration of radionuclides. This work looks at those processes as was also done with recent studies. Simulations of thermal-hydrology for the emplacement of cesium and strontium capsules in a deep borehole are presented. The simulations looked at disposal options such as different disposal configurations and aging of waste to reduce thermal effects. The simulations studied temperature and fluid flux in the vicinity of the borehole. Simulation results include temperature and vertical flux profiles around the borehole at selected depths. Of particular importance are peak temperatures, and fluxes above the disposal zone. (authors)

[en] Recently the U.S. made a policy decision to begin R and D for geologic disposal of defense-related HLW/SNF in a facility separate from commercially generated waste. A deep geologic repository for the disposal of DOE-managed HLW, and some thermally cooler DOE-managed SNF, arising from defense and DOE R and D activities, is part of a stepwise, phased approach for disposal of the Nation's nuclear waste, as recommended by the Blue Ribbon Commission on America's Nuclear Future. The work discussed here focuses on post-closure safety (or performance) assessment for such a Defense Waste Repository (DWR), which is divided into four major activities: (1) development of generic reference cases (i.e., knowledge or technical bases for 'generic' or 'non-site-specific' deep geologic repositories); (2) features, events, and processes (FEPs) analyses and screening to support the technical bases and the performance assessment (PA) model; (3) performance evaluation of alternative EBS design concepts; and (4) post-closure safety analyses of the repository system under consideration. Using the known inventory of defense-related SNF, as well as defense-related HLW stored at the Savannah River and Hanford sites, the Geologic Disposal Safety Assessment (GDSA) Framework modeling and software system has been applied to simulate the potential performance of a DWR in crystalline host rock, resulting in a suite of single-realization (i.e., deterministic) and multi-realization (i.e., probabilistic) 3-D post-closure system analyses, over a performance period of one million years. Two types of emplacement concepts are examined, including single-canister vertical-borehole emplacement for the hotter defense SNF waste (KBS-3V concept) and multi-canister horizontal emplacement for defense HLW (similar to Yucca Mountain co-disposal waste packages). Sensitivity analyses examine the effect of key uncertain parameters on repository performance, including the effects of fracture distribution, waste package degradation rate, buffer and disturbed rock zone (DRZ) properties, and sorption parameters. Initial results indicate that a crystalline host rock with a connected fracture system may require additional safety features to ensure robustness of the isolation safety function, such as a deep unsaturated zone, a sufficiently thick sedimentary overburden, and/or a disposal overpack with a very slow corrosion rate. None of these should pose an undue obstacle for successful disposal. (authors)

[en] This work examines the possible migration of radionuclides from a deep geologic nuclear waste repository sited in fractured crystalline host rock. The key safety concern is the potential effect on waste isolation of the inter-connectivity of the fracture network, which is primarily established by the temporal evolution of the temperature and stress fields at the time of original rock intrusion and cooling. Two end members are considered in these simulations, one with a high degree of connectivity to the biosphere, such that advective transport through the fracture and fault network controls radionuclide migration, and the other with a low degree of connectivity, such that slow diffusive transport through the crystalline rock matrix is the controlling process for radionuclide migration to the biosphere. Both end members have several areally extensive, high-transmissivity deformation zones, whose distance to the repository is an important factor for waste isolation capability. Uncertainties in fracture properties (transmissivity, orientation, radius) and fracture distribution also give rise to uncertainty in waste isolation capability (i.e., the overall connectivity of the repository with the surface biosphere). These types of natural system uncertainties are generally present in any geologic site-characterization program and thus represent an important factor in assessing repository performance in hard rock environments. The waste considered in these simulations are spent fuel rod assemblies from the U.S. commercial nuclear reactor fleet. Thermal output of the waste must be considered in the simulations and can influence the rate and timing of waste package failure, waste form degradation, and fluid flux due to thermal expansion around the repository horizon. To achieve a high-fidelity representation of radionuclide transport in fractures and rock matrix, combined with thermal energy transport and fluid flow in fractures and matrix, the mathematical model is solved numerically in a parallel high-performance computing (HPC) environment on a finite volume unstructured grid consisting of approximately 5 million cells, using Geologic Disposal Safety Assessment (GDSA) Framework (https://pa.sandia.gov), an open-source performance assessment tool for deep underground disposal of nuclear waste. GDSA Framework uses PFLOTRAN to solve the balance equations on a three-dimensional grid with heterogeneous properties, using multiple processors in a parallel configuration based on domain decomposition. The fracture networks in these simulations are originally generated as discrete fracture networks (DFNs), which are sets of two-dimensional planes distributed in a three-dimensional domain. The method used in GDSA Framework maps the stochastically generated DFN to an equivalent continuous porous medium (ECPM) domain that allows for the simulation of coupled heat flow, fluid flow, and radionuclide transport, including heat conduction through the matrix of the fractured rock, which is a process not easily modeled in a DFN representation. Computational efficiency is also greatly enhanced using the ECPM method, allowing for a realistic representation and analysis of uncertainties in a multi-realization performance assessment of a deep geologic repository. The effect of fracture connectivity on the waste isolation safety function, as brought to light by these GDSA Framework simulations, points to the importance of including a realistic representation of uncertainties in fracture properties and distribution (effectively, uncertainty in spatial heterogeneity) in repository safety assessment simulations. (authors)

[en] Development of an enhanced performance assessment (PA) capability for geologic disposal of spent nuclear fuel and high-level waste has been ongoing for several years in the U.S. repository program. The new Generic Disposal System Analysis (GDSA) modeling and software framework is intended to be flexible enough to evolve through the various phases of repository activities, beginning with generic PA activities in the current Concept Evaluation phase to site-specific PA modeling in the Repository Development phase. The GDSA Framework utilizes modern software and hardware capabilities by being based on open-source software architecture and being configured to run in a massively parallel, high-performance computing (HPC) environment. It consists of two main components, the open-source Dakota uncertainty sampling and analysis software and the PFLOTRAN reactive multi-phase flow and transport simulator. Reference cases or 'generic repositories' have been, and are being developed, based on typical properties for potential salt, clay, and granite host-rock formations and corresponding engineered design concepts for each medium. Past simulations have focused on a generic repository in bedded-salt host rock, while the most recent research has focused on a reference case for a typical clay/shale host rock. A variety of single-realization (i.e., deterministic) and multi-realization (probabilistic) results for the new clay reference case are presented, including an analysis of the effects of heat generation on repository performance, assuming a 100-year out-of-reactor commercial SNF waste form. Order-of-magnitude differences between predicted radionuclide concentrations in thermal versus isothermal simulations imply that mechanistic, coupled-process modeling in three-dimensional (3-D) domains can be important for building confidence in post-closure performance assessments. (authors)