Adsorption of noble gas fission products onto naturally occurring minerals is of interest for its potential to retain or retard emissions from nuclear fuel reprocessing operations or underground nuclear explosions. However, experimental studies of trace noble gas adsorption in the presence of air and water have largely focused on synthetic materials, such as activated carbon or metal-organic frameworks. Here, adsorption of Kr and Xe onto the naturally occurring zeolitic mineral clinoptilolite is studied in the presence of nitrogen and water. By varying the composition of the gas phase and monitoring the change in the combined adsorbate mass, the adsorbed concentration of noble gas is calculated gravimetrically. For dry clinoptilolite, the concentration of adsorbed Kr and Xe is linearly correlated with noble gas pressure and Henry's Law appears satisfactory, despite the presence of nitrogen at atmospheric pressures. However, the presence of water significantly reduces the adsorbed concentration of both Kr and Xe, which is typical in nanoporous sorbents. Here, an empirical bivariate model is presented, combining the Henry's Law adsorption model for a dry adsorbent with the exponential reduction in the presence of water, as reported by Lungu and Underhill in 1999. This model provides a means to estimate the adsorbate concentration at the trace partial pressures and higher water contents relevant to field-scale modeling of fission gas transport through the vadose zone.
Geogenic Helium-4 (4He) in-situ increases locally in regions of large deformation generated naturally or anthropogenically. This gas release by deformation is a potential geochemical precursor signal for subsurface deformation. To evaluate the applicability of 4He degassing for correlating deformation in different lithologies, we conducted high force crush tests, up to 97,800 N axial load, to assess the total 4He released during fragmentation of the rocks. We observed that the highest 4He released occurred in the sedimentary rocks and that release correlated strongly with lithologic age and U/Th content. Microstructural changes of the pre- and post-test rocks indicate that the degree of grain size reduction relates directly to the total 4He released during crushing. The range of in-place 4He was calculated based XRF measurements of uranium and thorium in each lithology, with the results indicating that the majority of the trapped 4He was not released. However, the 4He released by deformation depended upon how the each rock deformed during deformation and the degree of grain size reduction. We postulate that 4He precursor signals can be used to understand subsurface deformation only if geomechanical and geochemical conditions for 4He enrichment in a lithology are met.
This report presents development of a new meso-scale geomechanics capability using an 8.84-m radius centrifuge for physically simulating underground chemical explosions and stress cage processes in scaled models.
A series of drained and undrained water-saturated constant mean-stress tests were performed to investigate the strength, elasticity, and poroelastic response of a water-saturated high porosity nonwelded tuff. Drained strengths are found to increase with increasing effective confining pressures. Elastic moduli increase with increasing mean stress. Undrained strengths are small due to development of high pore pressures that generate low effective confining pressures. Skempton’s values are pressure dependent and appear to reflect the onset of inelastic deformation. Permeabilities decrease after deformation from ∼ 10–14 to ∼ 10–16 m2 and are a function of the applied confining pressure. Deformation is dominated by pore collapse, compaction, and intense microfracturing, with the undrained tests favoring microfracture-dominant deformation and the drained tests favoring compaction-dominant deformation. These property determinations and observations are used to develop/parameterize physics-based models for underground explosives testing.
Underground chemical explosive experiments such as LYNM PE1 generate large multiphenomenological datasets, require complex site preparation and build out, and utilize cutting edge models and analysis techniques to analyze and simulate the explosion-induced signals. This wide range of outcomes makes it a necessity to thoroughly characterize the testbed in advance of experiments in a way that complements the wide suite of data being generated. Here, we present a broad overview of the site characterization work and data collection that was conducted before Experiment A, which is the first in a series of three PE1 experiments. This work includes, but is not limited to, geologic mapping, physical sample collection, analysis of material properties, geophysical borehole logging, and in-situ measurements. This information was collected by a large, dedicated team and was used to inform site construction, finalize instrumentation placement, generate Geologic Framework Models, feed pre-experiment predictions, and facilitate post-experiment data analysis
This data documentation report describes geologic and hydrologic laboratory analysis and data collected in support of site characterization of the Physical Experiment 1 (PE1) testbed, Aqueduct Mesa, Nevada. The documentation includes a summary of laboratory tests performed, discussion of sample selection for assessing heterogeneity of various testbed properties, methods, and results per data type.
The goal of this work is to provide a database of quality-checked seismic parameters that can be integrated with the Geologic Framework Model (GFM) for the LYNM-PE1 (Low Yield Nuclear Monitoring – Physical Experiment 1) testbed. We integrated data from geophysical borehole logs, tabletop measurements on collected core, and laboratory measurements. We reviewed for internal consistency among each measurement type, documented the caveats of measurement conditions, and integrated lithologic logs to check the validity of outlier values. The resulting consolidated parameter tables can be used as inputs for modeling and analysis codes and are designed to interface with the GFM, which is being actively developed.
Geogenic gases often reside in intergranular pore space, fluid inclusions, and within mineral grains. In particular, helium-4 (4He) is generated by alpha decay of uranium and thorium in rocks. The emitted 4He nuclei can be trapped in the rock matrix or in fluid inclusions. Recent work has shown that releases of helium occur during plastic deformation of crustal rocks above atmospheric concentrations that are detectable in the field. However, it is unclear how rock type and deformation modalities affect the cumulative gas released. This work seeks to address how different deformation modalities observed in several rock types affect release of helium. Axial compression tests with granite, rhyolite, tuff, dolostone, and sandstone - under vacuum conditions - were conducted to measure the transient release of helium from each sample during crushing. It was found that, when crushed up to 97500 N, each rock type released helium at a rate quantifiable using a helium mass spectrometer leak detector. For plutonic rock like granite, helium flow rate spikes with the application of force as the samples elastically deform until fracture, then decays slowly until grain breakdown comminution begins to occur. Both the rhyolite and tuff do not experience such large spikes in helium flow rate, with the rhyolites fracturing at much lower force and the tuffs compacting instead of fracturing due to their high porosity. Both rhyolite and tuff instead experience a lesser but steady helium release as they are crushed. The cumulative helium release for the volcanic tuffs varies as much as two orders of magnitude but is fairly consistent for the denser rhyolite and granite tested. The results indicate that there is a large degassing of helium as rocks are elastically and inelastically deformed prior to fracturing. For more porous and less brittle rocks, the cumulative release will depend more on the degree of deformation applied. These results are compared with known U/Th radioisotopes in the rocks to relate the trapped helium as either produced in the rock or from secondary migration of 4He.
Geogenic gases often reside in intergranular pore space, fluid inclusions, and within mineral grains. In particular, helium-4 (4He) is generated by alpha decay of uranium and thorium in rocks. The emitted 4He nuclei can be trapped in the rock matrix or in fluid inclusions. Recent work has shown that releases of helium occur during plastic deformation of crustal rocks above atmospheric concentrations that are detectable in the field. However, it is unclear how rock type and deformation modalities affect the cumulative gas released. This work seeks to address how different deformation modalities observed in several rock types affect release of helium. Axial compression tests with granite, rhyolite, tuff, dolostone, and sandstone - under vacuum conditions - were conducted to measure the transient release of helium from each sample during crushing. It was found that, when crushed up to 97500 N, each rock type released helium at a rate quantifiable using a helium mass spectrometer leak detector. For plutonic rock like granite, helium flow rate spikes with the application of force as the samples elastically deform until fracture, then decays slowly until grain breakdown comminution begins to occur. Both the rhyolite and tuff do not experience such large spikes in helium flow rate, with the rhyolites fracturing at much lower force and the tuffs compacting instead of fracturing due to their high porosity. Both rhyolite and tuff instead experience a lesser but steady helium release as they are crushed. The cumulative helium release for the volcanic tuffs varies as much as two orders of magnitude but is fairly consistent for the denser rhyolite and granite tested. The results indicate that there is a large degassing of helium as rocks are elastically and inelastically deformed prior to fracturing. For more porous and less brittle rocks, the cumulative release will depend more on the degree of deformation applied. These results are compared with known U/Th radioisotopes in the rocks to relate the trapped helium as either produced in the rock or from secondary migration of 4He.
A credible simulation of disposal room porosity at the Waste Isolation Pilot Plant (WIPP) requires a tenable compaction model for the 55-gallon waste containers within the room. A review of the legacy waste material model, however, revealed several out-of-date and untested assumptions that could affect the model’s compaction behavior. For example, the legacy model predicted non-physical tensile out-of-plane stresses under plane strain compression. (Plane strain compression is similar to waste compaction in the middle of a long drift.) Consequently, a suite of new compaction experiments were performed on containers filled with surrogate, non-degraded, waste. The new experiments involved uniaxial, triaxial, and hydrostatic compaction tests on quarter-scale and full-scale containers. Special effort was made to measure the volume strain during uniaxial and triaxial tests, so that the lateral strain could be inferred from the axial and volume strain. These experimental measurements were then used to calibrate a pressure dependent, viscoplastic, constitutive model for the homogenized compaction behavior of the waste containers. This new waste material model’s predictions agreed far better with the experimental measurements than the legacy model’s predictions, especially under triaxial and hydrostatic conditions. Under plane strain compression, the new model predicted reasonable compressive out-of-plane stresses, instead of tensile stresses. Moreover, the new model’s plane strain behavior was substantially weaker for the same strain, yet substantially stronger for the same porosity, than the legacy model’s behavior. Although room for improvement exists, the new model appears ready for prudent engineering use.
The Source Physics Experiment (SPE) is a long-term NNSA research and development effort designed to improve nonproliferation verification and monitoring capabilities. The overarching goals of the SPE program are to improve understanding of prompt signals and physical signatures that develop from underground chemical explosions and associated modeling capabilities. Our work focuses on a primary factor controlling chemical explosion induced signals and signatures: the material properties of the rocks in which the chemical explosion takes place. This document reports on material property determinations of legacy core USGS Test Well F and outcrop analogs for the subsurface stratigraphy for the third phase of SPE in the Rock Valley (RV) area of the Nevada National Security Site (NNSS). The objective of this work is to establish a baseline set of lithologic descriptions and material properties expected prior to observatory borehole drilling in support of the SPE-RVDC (Rock Valley Direct Comparison) experiment. We determine for each rock type the compressional failure envelope, elastic properties as a function of stress (bulk modulus versus mean stress, shear modulus versus shear stress, Young’s modulus versus axial stress and Poisson’s ratio versus axial stress), indirect tensile strength, and porosity. Geologic characterization, both at the core-scale and microscale, provides context for using the data in modeling efforts and to inform interpretations for the material properties testing.
Prediction of flow, transport, and deformation in fractured and porous media is critical to improving our scientific understanding of coupled thermal-hydrological-mechanical processes related to subsurface energy storage and recovery, nonproliferation, and nuclear waste storage. Especially, earth rock response to changes in pressure and stress has remained a critically challenging task. In this work, we advance computational capabilities for coupled processes in fractured and porous media using Sandia Sierra Multiphysics software through verification and validation problems such as poro-elasticity, elasto-plasticity and thermo-poroelasticity. We apply Sierra software for geologic carbon storage, fluid injection/extraction, and enhanced geothermal systems. We also significantly improve machine learning approaches through latent space and self-supervised learning. Additionally, we develop new experimental technique for evaluating dynamics of compacted soils at an intermediate scale. Overall, this project will enable us to systematically measure and control the earth system response to changes in stress and pressure due to subsurface energy activities.
Mechanical properties on alluvium blocks and core samples were determined to support the Source Physics Experiment Dry Alluvium Geology experimental series. Because material was not available directly from the experimental location, the alluvium blocks and core samples are intended to serve as surrogate material . P - and S - wave velocity was measured on cubes cut from the alluvium blocks and core with the intention to study variation from water content and measured direction (material anisotropy). Indirect tensile tests were conducted dry and with moisture ranging from 6 to 9.1%. For the range of water content tested, increasing moisture level resulted in slower P - and S - wave velocities. P - and S - wave variability is less influenced by material heterogeneity than moisture content. P - wave velocity ranges from 629 m/s to 2599 m/s and S - wave velocity ranges from 288 m/s to 1200 m/s. Counter to the velocity measurement findings, material variability on indirect tensile strength has a greater effect than moisture content. Compared to dry strength and at moisture levels from 6 to 9% the block's tensile strength was lowered by at least a factor of 5. Indirect tensile strength for the first block averaged 0.35 MPa and 0.25 MPa for dry and 8.9% moisture respectively. For the second block indirect tensile strength averaged 0.05 MPa for both dry and 6.4% moisture.
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.
Natural and induced fractures are potential preferential pathways for migration of radioactive gases to earths surface from underground nuclear explosions (UNEs). This report documents X-ray computed tomography (XRCT) imaging on 26 samples of rock core that was collected to support the Underground Nuclear Explosion Signatures Experiment (UNESE) program. The XRCT datasets are intended to help fill a data gap on the three-dimensional (3D) characteristics of natural and/or induced fractures at the centimeter and smaller scale, which may strongly influence multiphase fluid flow and transport properties of preferential flow paths and interaction with the matrix of the surrounding host rock. Pre- and post-UNE rock samples were carefully chosen to enable comparison of fractures as a function of lithologic and petrophysical properties, as well as distance to the past UNEs. This report serves as documentation for the data, including an introduction with the research motivation, a methods and materials section, descriptions of the XRCT datasets without post-processing, and recommendations for 3D quantification via image analysis and digital rock physics.
