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).
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.
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.
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.
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.
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.
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.