Nonlocal models provide a much-needed predictive capability for important Sandia mission applications, ranging from fracture mechanics for nuclear components to subsurface flow for nuclear waste disposal, where traditional partial differential equations (PDEs) models fail to capture effects due to long-range forces at the microscale and mesoscale. However, utilization of this capability is seriously compromised by the lack of a rigorous nonlocal interface theory, required for both application and efficient solution of nonlocal models. To unlock the full potential of nonlocal modeling we developed a mathematically rigorous and physically consistent interface theory and demonstrate its scope in mission-relevant exemplar problems.
Fractured media models comprise discontinuities of multiple lengths (e.g. fracture lengths and apertures, wellbore area) that fall into the relatively insignificant length scales spanning millimeter-scale fractures to centimeter-scale wellbores in comparison to the extensions of the field of interest, and challenge the conventional discretization methods imposing highly-fine meshing and formidably large numerical cost. By utilizing the recent developments in the finite element analysis of electromagnetics that allow to represent material properties on a hierarchical geometry, this project develops computational capabilities to model fluid flow, heat conduction, transport and induced polarization in large-scale geologic environments that possess geometrically-complex fractures and man-made infrastructures without explosive computational cost. The computational efficiency and robustness of this multi-physics modeling tool are demonstrated by considering various highly-realistic complex geologic environments that are common in many energy and national security related engineering problems.
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.
This report summarizes the international collaborations conducted by Sandia funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-22SN010303063. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), numerical model comparison (DECOVALEX) and an NEA Salt Club working group on the development of scenarios as part of the performance assessment development process. Finally, we summarize events related to the US/German Workshop on Repository Research, Design and Operations. The work summarized in this annual update has occurred during the COVID-19 pandemic, and little international or domestic travel has occurred. Most of the collaborations have been conducted via email or as virtual meetings, but a slow return to travel and in-person meetings has begun.
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.
The capability of a 1-D PFLOTRAN model to simulate the S1-3 bentonite saturation experiment has been demonstrated and validated against experimental data. Work remains to be done to refine 1-D PFLOTRAN simulations of the experiment S1-4 which include evaluation of parameter sensitivities on the prediction of material saturation and relative permeabilities. This and further testing of PFLOTRAN capabilities will be done as part of DECOVALEX 2023 Task D contributions by the SNL team in the coming months.
Two-phase fluid flow properties underlie quantitative prediction of water and gas movement, but constraining these properties typically requires multiple time-consuming laboratory methods. The estimation of two-phase flow properties (van Genuchten parameters, porosity, and intrinsic permeability) is illustrated in cores of vitric nonwelded volcanic tuff using Bayesian parameter estimation that fits numerical models to observations from spontaneous imbibition experiments. The uniqueness and correlation of the estimated parameters is explored using different modeling assumptions and subsets of the observed data. The resulting estimation process is sensitive to both moisture retention and relative permeability functions, thereby offering a comprehensive method for constraining both functions. The data collected during this relatively simple laboratory experiment, used in conjunction with a numerical model and a global optimizer, result in a viable approach for augmenting more traditional capillary pressure data obtained from hanging water column, membrane plate extractor, or mercury intrusion methods. This method may be useful when imbibition rather than drainage parameters are sought, when larger samples (e.g., including heterogeneity or fractures) need to be tested that cannot be accommodated in more traditional methods, or when in educational laboratory settings.
The construction of deep geological repositories (DGR) in salt formations requires penetrating through naturally sealing geosphere layers. While the emplaced nuclear waste is primarily protected by the containment-providing rock zone (CRZ), technical barriers are required, for example during handling. For closure geotechnical barriers seal the repository along the accesses against water or solutions from outside and the possible emission paths for radionuclides contained inside. As these barriers must ensure maintenance-free function on a long-term basis, they typically comprise a set of specialized elements with diversified functions that may be used redundantly. The effects of the individual elements are coordinated so that they are collectively referred to as the Engineered Barrier System (EBS).
Elevated temperature and pressure in the earth's subsurface alters the permeability of salt formations, due to changing properties of the salt-brine interface. Molecular dynamics (MD) simulations are used to investigate the mechanisms of temperature and pressure dependence of liquid-solid interfacial tensions of NaCl, KCl, and NaCl-KCl brines in contact with (100) salt surfaces. Salt-brine dihedral angles vary between 55 and 76° across the temperature (300-450 K) and pressure range (0-150 MPa) evaluated. Temperature-dependent brine composition results in elevated dihedral angles of 65-80°, which falls above the reported salt percolation threshold of 60°. Mixed NaCl-KCl brine compositions increased this effect. Elevated temperatures excluded dissolved Na+ ions from the interface, causing the strong temperature dependence of the liquid-solid interfacial tension and the resulting dihedral angle. Therefore, at higher temperature, pressure, and brine concentrations Na-Cl systems may underpredict the dihedral angle. Higher dihedral angles in more realistic mixed brine systems maintain low permeability of salt formations due to changes in the structure and energetics of the salt-brine interface.
This report summarizes the 2021 fiscal year (FY21) status of ongoing borehole heater tests in salt funded by the disposal research and development (R&D) program of the Office of Spent Fuel & Waste Science and Technology (SFWST) of the US Department of Energy’s Office of Nuclear Energy’s (DOE-NE) Office of Spent Fuel and Waste Disposition (SFWD). This report satisfies SFWST milestone M2SF- 21SN010303052 by summarizing test activities and data collected during FY21. The Brine Availability Test in Salt (BATS) is fielded in a pair of similar arrays of horizontal boreholes in an experimental area at the Waste Isolation Pilot Plant (WIPP). One array is heated, the other unheated. Each array consists of 14 boreholes, including a central borehole with gas circulation to measure water production, a cement seal exposure test, thermocouples to measure temperature, electrodes to infer resistivity, a packer-isolated borehole to add tracers, fiber optics to measure temperature and strain, and piezoelectric transducers to measure acoustic emissions. The key new data collected during FY21 include a series of gas tracer tests (BATS phase 1b), a pair of liquid tracer tests (BATS phase 1c), and data collected under ambient conditions (including a period with limited access due to the ongoing pandemic) since BATS phase 1a in 2020. A comparison of heated and unheated gas tracer test results clearly shows a decrease in permeability of the salt upon heating (i.e., thermal expansion closes fractures, which reduces permeability).
The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heatgenerating nuclear waste disposition in deep-seated clay/shale/argillaceous rock. International collaboration activities such as heater tests, continuous field data monitoring, and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum and discrete) to tackle issues related to flow and transport at various scales of the host-rock, its interactions with barrier materials, and EBS design concept.
The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.
We present a dynamic laboratory spontaneous imbibition test and interpretation method, demonstrated on volcanic tuff samples from the Nevada National Security Site. The method includes numerical inverse modeling to quantify uncertainty of estimated two-phase fluid flow properties. As opposed to other approaches requiring multiple different laboratory instruments, the dynamic imbibition method simultaneously estimates capillary pressure and relative permeability from one test apparatus.
This report summarizes the international collaboration work conducted by Sandia and funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-20SN010303062. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), and model comparison (DECOVALEX). Lastly, the report summarizes a newly developed working group on the development of scenarios as part of the performance assessment development process, and the activities related to the Nuclear Energy Agency (NEA) Salt club and the US/German Workshop on Repository Research, Design and Operations.
Non-uniqueness in groundwater model calibration is a primary source of uncertainty in groundwater flow and transport predictions. In this study, we investigate the ability of environmental tracer information to constrain groundwater model parameters. We utilize a pilot point calibration procedure conditioned to subsets of observed data including: liquid pressures, tritium (3H), chlorofluorocarbon-12 (CFC-12), and sulfur hexafluoride (SF6) concentrations; and groundwater apparent ages inferred from these environmental tracers, to quantify uncertainties in the heterogeneous permeability fields and infiltration rates of a steady-state 2-D synthetic aquifer and a transient 3-D model of a field site located near Riverton, Wyoming (USA). To identify the relative data worth of each observation data type, the post-calibration uncertainties of the optimal parameters for a given observation subset are compared to that from the full observation data set. Our results suggest that the calibration-constrained permeability field uncertainties are largest when liquid pressures are used as the sole calibration data set. We find significant reduction in permeability uncertainty and increased predictive accuracy when the environmental tracer concentrations, rather than apparent groundwater ages, are used as calibration targets in the synthetic model. Calibration of the Riverton field site model using environmental tracer concentrations directly produces infiltration rate estimates with the lowest uncertainties, however; permeability field uncertainties remain similar between the environmental tracer concentration and apparent groundwater age calibration scenarios. This work provides insight on the data worth of environmental tracer information to calibrate groundwater models and highlights potential benefits of directly assimilating environmental tracer concentrations into model parameter estimation procedures.
We present a new pre-processor tool written in Python that creates multicontinuum meshes for PFLOTRAN to simulate two-phase flow and transport in both the fracture and matrix continua. We discuss the multicontinuum modeling approach to simulate potentially mobile water and gas in the fractured volcanic tuffs at Aqueduct Mesa, at the Nevada National Security Site.
This interim report is an update of ongoing experimental and modeling work on bentonite material described in Jové Colón et al. (2019, 2020) from past international collaboration activities. As noted in Jové Colón et al. (2020), work on international repository science activities such as FEBEX-DP and DECOVALEX19 is either no longer continuing by the international partners. Nevertheless, research activities on the collected sample materials and field data are still ongoing. Descriptions of these underground research laboratory (URL) R&D activities are described elsewhere (Birkholzer et al. 2019; Jové Colón et al. 2020) but will be explained here when needed. The current reports recent reactive-transport modeling on the leaching of sedimentary rock.
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.
Pore-scale finite-volume continuum models of electrokinetic processes are used to predict the Debye lengths, velocity, and potential profiles for two-dimensional arrays of circles, ellipses and squares with different orientations. The pore-scale continuum model solves the coupled Navier–Stokes, Poisson, and Nernst–Planck equations to characterize the electro-osmotic pressure and streaming potentials developed on the application of an external voltage and pressure difference, respectively. This model is used to predict the macroscale permeabilities of geomaterials via the widely used Carmen–Kozeny equation and through the electrokinetic coupling coefficients. The permeability results for a two-dimensional X-ray tomography-derived sand microstructure are within the same order of magnitude as the experimentally calculated values. The effect of the particle aspect ratio and orientation on the electrokinetic coupling coefficients and subsequently the electrical and hydraulic tortuosity of the porous media has been determined. These calculations suggest a highly tortuous geomaterial can be efficient for applications like decontamination and desalination.
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.
Brine availability in salt has multiple implications for the safety and design of a nuclear waste storage facility. Brine availability includes both the distribution and transport of brine through a damaged zone around boreholes or drifts excavated into the salt. Coupled thermal, hydrological, mechanical, and chemical processes taking place within heated bedded salt are complex; as part of DECOVALEX 2023 Task E this study takes a parsimonious modeling approach utilizing analytical and numerical one-dimensional simulations to match field measurements of temperature and brine inflow around a heater. The one-dimensional modeling results presented arrive at best-fit thermal conductivity of intact salt, and the permeability and porosity of damaged salt of 5.74 W/m · K, 10−17 m2, and ≈0.02, respectively.
This final report on Laboratory Directed Research and Development (LDRD) project 209234 presents background material for electrokinetics at the pore and porous media scales. We present some theoretical developments related to uncoupling electrokinetic flow solutions, from a manuscript recently accepted into Mathematical Geosciences for publication. We present a summary of two pore-scale modeling efforts undertaken as part of the academic alliance with University of Illinois, resulting in one already submitted journal publication to Transport in Porous Media and another in preparation for submission to a journal. We finally show the laboratory apparatus built in Laboratory B59 in Building 823 and discuss some of the issues that occurred with it.
This report is a summary of the international collaboration work conducted by Sandia and funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies milestone level-three milestone M3SF-205N010303062. Several stand-alone sections make up this summary report, each completed by the participants. The first two sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), and documentation of Features, Events, and Processes (FEPs).
Detection of radioxenon and radioargon produced by underground nuclear explosions is one of the primary methods by which the Comprehensive Nuclear-Test–Ban Treaty (CTBT) monitors for nuclear activities. However, transport of these noble gases to the surface via barometric pumping is a complex process relying on advective and diffusive processes in a fractured porous medium to bring detectable levels to the surface. To better understand this process, experimental measurements of noble gas and chemical surrogate diffusivity in relevant lithologies are necessary. However, measurement of noble gas diffusivity in tight or partially saturated porous media is challenging due to the transparent nature of noble gases, the lengthy diffusion times, and difficulty maintaining consistent water saturation. Here, the quasi-steady-state Ney–Armistead method is modified to accommodate continuous gas sampling via effusive flow to a mass spectrometer. An analytical solution accounting for the cumulative sampling losses and induced advective flow is then derived. Experimental results appear in good agreement with the proposed theory, suggesting the presence of retained groundwater reduces the effective diffusivity of the gas tracers by 10–1000 times. Furthermore, by using a mass spectrometer, the method described herein is applicable to a broad range of gas species and porous media.
This report is the Task E specification (Revision 0) for DECOVALEX-2023. Task E is focused on understanding thermal, hydrological, mechanical and chemical (THMC) processes, especially related to predicting brine migration in heated salt. The main test case being used is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) in Carlsbad, New Mexico. This report provides short motivational background, a summary of relevant experiments and data, and a step-by-step plan for the analysis by the teams participating in Task E (Rev. 0 includes detailed description of steps 0 and 1). This document will be revised, and more detail will be added to later steps during DECOVALEX-2023.
Permeability of salt formations is controlled by the equilibrium between the salt-brine and salt-salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. Here, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt-brine interface to provide insight into the stability of salt systems. Mixed NaCl-KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt-brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2-3 times (5.0 M systems). Size differences in Na+ and K+ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt-brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.
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.
We present an approach to uncoupling the pair of transient governing equations used in electrokinetics (i.e., streaming potential and electroosmosis). This approach allows for the solution of two uncoupled "intermediate" equations, then the physical solution is found by recombination of these intermediate potentials through a matrix multiplication. We present numerically stable expressions for the coefficients, and an example showing electrokinetics arising from pumping a fully penetrating well in a confined aquifer, surrounded by insulating aquicludes. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. (SAND2019-8712 A)
This report presents a generic (i.e., site-independent) preliminary plan for drilling, testing, sampling, and analyzing data for a deep characterization borehole drilled into crystalline basement for the purposes of assessing the suitability of a site for deep borehole disposal (DBD). This research was performed as part of the deep borehole field test (DBFT). Based on revised U.S. Department of Energy (DOE) priorities in mid-2017, the DBFT and other research related to a DBD option was discontinued; ongoing work and documentation were closed out by the end of fiscal year (FY) 2017. This report was initiated as part of the DBFT and documented as an incomplete draft at the end of FY 2017. The report was finalized by Sandia National Laboratories in FY2018 without DOE funding, subsequent to the termination of the DBFT, and published in FY2019. This report presents a possible sampling, testing, and analysis campaign that could be carried out as part of a future project to quantify geochemical, geomechanical, geothermal, and geohydrologic conditions encountered at depths up to 5 km in crystalline basement.
This report summarizes the 2018 fiscal year (FY18) field, laboratory, and modeling work funded by the US Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) campaign as part of the Sandia National Laboratories Salt Research and Development (R&D) and Salt International work packages. This report satisfies level-two milestone M2SF-18SNO10303031and comprises three related but stand-alone sections. The first section summarizes the programmatic progress made to date in the DOE-NE salt program and its goals going forward. The second section presents brine composition modeling and laboratory activities related to salt evaporation experiments, which will be used to interpret data collected during the heater test. The third section presents theoretical and numerical modeling work done to investigate the effects brine composition have on dihedral angle and the permeability of salt.
This project plan gives a high-level description of the US Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Disposition (SFWD) campaign in situ borehole heater test project being planned for the Waste Isolation Pilot Plant (WIPP) site This plan provides an overview of the schedule and responsibilities of the parties involved. This project is a collaborative effort by Sandia, Los Alamos, and Lawrence Berkeley National Laboratories to execute a series of small-diameter borehole heater tests in salt for the DOE-NE SFWD campaign. Design of a heater test in salt at WIPP has evolved over several years. The current design was completed in fiscal year 2017 (FY17), an equipment shakedown experiment is underway in April FY18, and the test implementation will begin in summer of FY18. The project comprises a suite of modular tests, which consist of a group of nearby boreholes in the wall of drifts at WIPP. Each test is centered around a packer-isolated heated borehole (5" diameter) containing equipment for water-vapor collection and brine sampling, surrounded by smaller-diameter (2" diameter) satellite observation boreholes. Observation boreholes will contain temperature sensors, tracer release points, electrical resistivity tomography (ERT) sensors, fiber optic sensing, and acoustic emission (AE) measurements, and sonic velocity sources and sensors. These satellite boreholes will also be used for plugging/sealing tests. The first two tests to be implemented will have the packer-isolated borehole heated to 120°C, with one observation borehole used to monitor changes. Follow-on tests will be designed using information gathered from the first two tests, will be conducted at other temperatures, will use multiple observation boreholes, and may include other measurement types and test designs.
This document records the Proceedings of the 2017 gathering of salt repository nations. In a spirit of mutual support, technical issues are dissected, led capably by subject matter experts. As before, it is not possible to explore all contemporary issues regarding nuclear waste disposal in salt formations. Instead, the group focused on a few selected issues to be pursued in depth, while at the same time acknowledging and recording ancillary issues.
We use helium released during mechanical deformation of shales as a signal to explore the effects of deformation and failure on material transport properties. A dynamic dual-permeability model with evolving pore and fracture networks is used to simulate gases released from shale during deformation and failure. Changes in material properties required to reproduce experimentally observed gas signals are explored. We model two different experiments of 4He flow rate measured from shale undergoing mechanical deformation, a core parallel to bedding and a core perpendicular to bedding. We find that the helium signal is sensitive to fracture development and evolution as well as changes in the matrix transport properties. We constrain the timing and effective fracture aperture, as well as the increase in matrix porosity and permeability. Increases in matrix permeability are required to explain gas flow prior to macroscopic failure, and the short-term gas flow postfailure. Increased matrix porosity is required to match the long-term, postfailure gas flow. Our model provides the first quantitative interpretation of helium release as a result of mechanical deformation. The sensitivity of this model to changes in the fracture network, as well as to matrix properties during deformation, indicates that helium release can be used as a quantitative tool to evaluate the state of stress and strain in earth materials.
This report is a summary of the international collaboration and laboratory work funded by the US Department of Energy Office of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D work package. This report satisfies milestone levelfour milestone M4SF-17SN010303014. Several stand-alone sections make up this summary report, each completed by the participants. The first two sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS) and bedded salt investigations (KOSINA), while the last three sections discuss laboratory work conducted on brucite solubility in brine, dissolution of borosilicate glass into brine, and partitioning of fission products into salt phases.
The groundwater flow system in the Culebra Dolomite Member (Culebra) of the Permian Rustler Formation is a potential radionuclide release pathway from the Waste Isolation Pilot Plant (WIPP), the only deep geological repository for transuranic waste in the United States. In early conceptual models of the Culebra, groundwater levels were not expected to fluctuate markedly, except in response to long-term climatic changes, with response times on the order of hundreds to thousands of years. Recent groundwater pressures measured in monitoring wells record more than 25 m of drawdown. The fluctuations are attributed to pumping activities at a privately owned well that may be associated with the demand of the Permian Basin hydrocarbon industry for water. The unprecedented magnitude of drawdown provides an opportunity to quantitatively assess the influence of unplanned anthropogenic forcings near the WIPP. Spatially variable realizations of Culebra saturated hydraulic conductivity and specific storage were used to develop groundwater flow models to estimate a pumping rate for the private well and investigate its effect on advective transport. Simulated drawdown shows reasonable agreement with observations (average Model Efficiency coefficient = 0.7). Steepened hydraulic gradients associated with the pumping reduce estimates of conservative particle travel times across the domain by one half and shift the intersection of the average particle track with the compliance boundary by more than 2 km. The value of the transient simulations conducted for this study lies in their ability to (a) improve understanding of the Culebra groundwater flow system and (b) challenge the notion of time-invariant land use in the vicinity of the WIPP.
This report summarizes the first stage in a collaborative effort by Sandia, Los Alamos, and Lawrence Berkeley National Laboratories to design a small-diameter borehole heater test in salt at the Waste Isolation Pilot Plant (WIPP) for the US Department of Energy Office of Nuclear Energy (DOE-NE). The intention is to complete test design during the remainder of fiscal year 2017 (FY17), and the implementation of the test will begin in FY18. This document is the result of regular meetings between the three national labs and the DOE-NE, and is intended to represent a consensus of these meetings and discussions.
ANS IHLRWM 2017 - 16th International High-Level Radioactive Waste Management Conference: Creating a Safe and Secure Energy Future for Generations to Come - Driving Toward Long-Term Storage and Disposal
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.
ANS IHLRWM 2017 - 16th International High-Level Radioactive Waste Management Conference: Creating a Safe and Secure Energy Future for Generations to Come - Driving Toward Long-Term Storage and Disposal
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.
Deep Borehole Disposal (DBD) of high-level radioactive wastes has been considered an option for geological isolation for many years (Hess et al. 1957). Recent advances in drilling technology have decreased costs and increased reliability for large-diameter (i.e., ≥50 cm [19.7”]) boreholes to depths of several kilometers (Beswick 2008; Beswick et al. 2014). These advances have therefore also increased the feasibility of the DBD concept (Brady et al. 2009; Cornwall 2015), and the current field test design will demonstrate the DBD concept and these advances. The US Department of Energy (DOE) Strategy for the Management and Disposal of Used Nuclear Fuel and High-Level Radioactive Waste (DOE 2013) specifically recommended developing a research and development plan for DBD. DOE sought input or expression of interest from States, local communities, individuals, private groups, academia, or any other stakeholders willing to host a Deep Borehole Field Test (DBFT). The DBFT includes drilling two boreholes nominally 200m [656’] apart to approximately 5 km [16,400’] total depth, in a region where crystalline basement is expected to begin at less than 2 km depth [6,560’]. The characterization borehole (CB) is the smaller-diameter borehole (i.e., 21.6 cm [8.5”] diameter at total depth), and will be drilled first. The geologic, hydrogeologic, geochemical, geomechanical and thermal testing will take place in the CB. The field test borehole (FTB) is the larger-diameter borehole (i.e., 43.2 cm [17”] diameter at total depth). Surface handling and borehole emplacement of test package will be demonstrated using the FTB to evaluate engineering feasibility and safety of disposal operations (SNL 2016).
A modified version of a published slug test model for unconfined aquifers is applied to cross-hole slug test data collected in field tests conducted at the Widen site in Switzerland. The model accounts for water-table effects using the linearized kinematic condition. The model also accounts for inertial effects in source and observation wells. The primary objective of this work is to demonstrate applicability of this semi-analytical model to multi-well and multi-level pneumatic slug tests. The pneumatic perturbation was applied at discrete intervals in a source well and monitored at discrete vertical intervals in observation wells. The source and observation well pairs were separated by distances of up to 4 m. The analysis yielded vertical profiles of hydraulic conductivity, specific storage, and specific yield at observation well locations. The hydraulic parameter estimates are compared to results from prior pumping and single-well slug tests conducted at the site, as well as to estimates from particle size analyses of sediment collected from boreholes during well installation. The results are in general agreement with results from prior tests and are indicative of a sand and gravel aquifer. Sensitivity analysis show that model identification of specific yield is strongest at late-time. However, the usefulness of late-time data is limited due to the low signal-to-noise ratios.
The R&D program from the DOE Used Fuel Disposition Campaign (UFDC) has documented key advances in coupled Thermal-Hydrological-Mechanical-Chemical (THMC) modeling of clay to simulate its complex dynamic behavior in response to thermal and hydrochemical feedbacks. These efforts have been harnessed to assess the isolation performance of heat-generating nuclear waste in a deep geological repository in clay/shale/argillaceous rock formations. This report describes the ongoing disposal R&D efforts on the advancement and refinement of coupled THMC process models, hydrothermal experiments on barrier clay interactions, used fuel and canister material degradation, thermodynamic database development, and reactive transport modeling of the near-field under non-isothermal conditions. These play an important role to the evaluation of sacrificial zones as part of the EBS exposure to thermally-driven chemical and transport processes. Thermal inducement of chemical interactions at EBS domains enhances mineral dissolution/precipitation but also generates mineralogical changes that result in mineral H2O uptake/removal (hydration/dehydration reactions). These processes can result in volume changes that can affect the interface / bulk phase porosities and the mechanical (stress) state of the bentonite barrier. Characterization studies on bentonite barrier samples from the FEBEX-DP international activity have provided important insight on clay barrier microstructures (e.g., microcracks) and interactions at EBS interfaces. Enhancements to the used fuel degradation model outlines the need to include the effects of canister corrosion due the strong influence of H2 generation on the source term.
Sandia National Laboratories has begun developing modeling and analysis tools of flow through the cemented port ion of a cemented annulus in a Strategic Petroleum Reserve (SPR) well since August of 201 5 . The goal of this work is to develop model s and testing procedures to diagnose the health of cemented annuli at SPR sites. In Fiscal Year 2016 (FY16), we have developed several tests and associated models that we believe are sufficient for this purpose. This report outlines progress made in FY16 and future work.
Salt formations hold promise for eternal removal of nuclear waste from our biosphere. Germany and the United States have ample salt formations for this purpose, ranging from flat-bedded formations to geologically mature dome structures. As both nations revisit nuclear waste disposal options, the choice between bedded, domal, or intermediate pillow formations is once again a contemporary issue. For decades, favorable attributes of salt as a disposal medium have been extoled and evaluated, carefully and thoroughly. Yet, a sense of discovery continues as science and engineering interrogate naturally heterogeneous systems. Salt formations are impermeable to fluids. Excavation-induced fractures heal as seal systems are placed or natural closure progresses toward equilibrium. Engineering required for nuclear waste disposal gains from mining and storage industries, as humans have been mining salt for millennia. This great intellectual warehouse has been honed and distilled, but not perfected, for all nuances of nuclear waste disposal. Nonetheless, nations are able and have already produced suitable license applications for radioactive waste disposal in salt. A remaining conundrum is site location. Salt formations provide isolation and geotechnical barriers reestablish impermeability after waste is placed in the geology. Between excavation and closure, physical, mechanical, thermal, chemical, and hydrological processes ensue. Positive attributes for isolation in salt have many commonalities independent of the geologic setting. In some cases, specific details of the environment will affect the disposal concept and thereby define interaction of features, events and processes, while simultaneously influencing scenario development. Here we identify and discuss high-level differences and similarities of bedded and domal salt formations. Positive geologic and engineering attributes for disposal purposes are more common among salt formations than are significant differences. Developing models, testing material, characterizing processes, and analyzing performance all have overlapping application regardless of the salt formation of interest.
We present a method of control variates for calculating improved estimates for mean performance quantities of interest, E(PQI) , computed from Monte Carlo probabilistic simulations. An example of a PQI is the concentration of a contaminant at a particular location in a problem domain computed from simulations of transport in porous media. To simplify the presentation, the method is described in the setting of a one- dimensional elliptical model problem involving a single uncertain parameter represented by a probability distribution. The approach can be easily implemented for more complex problems involving multiple uncertain parameters and in particular for application to probabilistic performance assessment of deep geologic nuclear waste repository systems. Numerical results indicate the method can produce estimates of E(PQI)having superior accuracy on coarser meshes and reduce the required number of simulations needed to achieve an acceptable estimate.
Sandia National Laboratories has begun research on the potential use of deep boreholes for the dis¬posal of radioactive waste. Characterization activities will focus on measurements and samples that are important for evaluating the long-term iso¬lation capability of the deep borehole disposal (DBD) concept. Engineering demonstration activities will focus on providing data to evaluate the concept’s operational safety and practicality. Procurement of a scientifically acceptable deep borehole field test (DBFT) site and a site management contractor is now under way.
Deep Borehole Disposal (DBD) of high-level radioactive wastes has been considered an option for geological isolation for many years (Hess et al. 1957). Recent advances in drilling technology have decreased costs and increased reliability for large-diameter (i.e., ≥50 cm [19.7”]) boreholes to depths of several kilometers (Beswick 2008; Beswick et al. 2014). These advances have therefore also increased the feasibility of the DBD concept (Brady et al. 2009; Cornwall 2015), and the current field test, introduced herein, is a demonstration of the DBD concept and these advances.
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.
This report presents efforts to develop the use of in situ naturally-occurring noble gas tracers to evaluate transport mechanisms and deformation in shale hydrocarbon reservoirs. Noble gases are promising as shale reservoir diagnostic tools due to their sensitivity of transport to: shale pore structure; phase partitioning between groundwater, liquid, and gaseous hydrocarbons; and deformation from hydraulic fracturing. Approximately 1.5-year time-series of wellhead fluid samples were collected from two hydraulically-fractured wells. The noble gas compositions and isotopes suggest a strong signature of atmospheric contribution to the noble gases that mix with deep, old reservoir fluids. Complex mixing and transport of fracturing fluid and reservoir fluids occurs during production. Real-time laboratory measurements were performed on triaxially-deforming shale samples to link deformation behavior, transport, and gas tracer signatures. Finally, we present improved methods for production forecasts that borrow statistical strength from production data of nearby wells to reduce uncertainty in the forecasts.
The amount of brine present in domal salt formation is far less than in bedded salts (e.g., 0.01 to 0.1% compared with 1 to 3%). In salt domes, shear deformation associated with diapirism has caused existing brine to coalesce, leading to flow and expulsion. Brine migration behavior was investigated in bedded salt at WIPP (Nowak and McTigue 1987, SAND87-0880), and in domal salt at Asse (Coyle et al. 1987, BMI/ONWI-624). Test methods were not standardized, and the tests involved large diameter boreholes (17 to 36 in. diameter) and large apparatus. The tested intervals were proximal to mined openings (within approximately 1 diameter) where in situ stresses are redistributed due to excavation. The tests showed that (1) brine inflow rates can range over at least 2 orders of magnitude for domal vs. bedded salt, (2) that brine inflow is strongly associated with clay and interbedded permeable layers in bedded salt, and (3) that measurement systems can readily collect very small quantities of moisture over time frames of 2 years or longer. Brine inflow rates declined slightly with time in both test series, but neither series approached a state of apparent depletion. This range of flow magnitude could be significant to repository design and performance assessment, especially if inflow rates can be predicted using stratigraphic and geomechanical inputs, and can be shown to approach zero in a predictable manner.
In situ tests implemented in a research facility mined from salt deposits, if planned appropriately, provide an opportunity to characterize the host rock before, during, and after excavation of test rooms. Characterization of the test bed is essential to interpret structural deformation, creation and evolution of the disturbed rock zone, and measurement of first-order hydromechanical properties as the salt evolves from an impermeable undisturbed state to a more-transmissive damaged state. The strategy expounded upon in this paper describes recommended geophysical measurements to characterize the initial state of a potential test bed and its evolution over the course of a field test. Discussion includes what measurements could be made, why the measurements would be made, how they are made, and how accurately they need to be made. Quantifiable parameters will establish field-scale boundary conditions and data quality objectives to characterize the test bed in an underground salt research facility.
While deep borehole disposal of nuclear waste should rely primarily on off-the-shelf technologies pioneered by the oil and gas and geothermal industries, the development of new science and technology will remain important. Key knowledge gaps have been outlined in the research roadmap for deep boreholes (B. Arnold et al, 2012, Research, Development, and Demonstration Roadmap for Deep Borehole Disposal, Sandia National Laboratories, SAND2012-8527P) and in a recent Deep Borehole Science Needs Workshop. Characterizing deep crystalline basement, understanding the nature and role of deep fractures, more precisely age-dating deep groundwaters, and demonstrating long-term performance of seals are all important topics of interest. Overlapping deep borehole and enhanced geothermal technology needs include: quantification of seal material performance/failure, stress measurement beyond the borehole, advanced drilling and completion tools, and better subsurface sensors. A deep borehole demonstration has the potential to trigger more focused study of deep hydrology, high temperature brine-rock interaction, and thermomechanical behavior.
A multiporosity extension of classical double and triple-porosity fractured rock flow models for slightly compressible fluids is presented. The multiporosity model is an adaptation of the multirate solute transport model of Haggerty and Gorelick (1995) to viscous flow in fractured rock reservoirs. It is a generalization of both pseudo steady state and transient interporosity flow double-porosity models. The model includes a fracture continuum and an overlapping distribution of multiple rock matrix continua, whose fracture-matrix exchange coefficients are specified through a discrete probability mass function. Semianalytical cylindrically symmetric solutions to the multiporosity mathematical model are developed using the Laplace transform to illustrate its behavior. The multiporosity model presented here is conceptually simple, yet flexible enough to simulate common conceptualizations of double and triple-porosity flow. This combination of generality and simplicity makes the multiporosity model a good choice for flow modelling in low-permeability fractured rocks. Key Points: New conceptual model generalizing dual porosity flow in fractured rocks Multiporosity model is generalization of existing well-known models Semianalytical radially symmetric solutions presented for new conceptual model