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.
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.
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.
Geophysical techniques are often implemented as quick and inexpensive ways to locate and characterize fractures in the subsurface, which is important for a number of geoscience fields. Seismic velocities are the most widely used proxies for identification of fractures, but the correlation is not always well-defined. In this study we present material property data: unconfined compressive strength (UCS), bulk density (ρ), Young's modulus (E), Poisson's ratio (ν), P wave velocity (Vp), and S wave velocity (Vs), in conjunction with microfracture densities measured on samples of granite collected before and after underground chemical explosions. Results indicate the relationship between fractures and material properties is complex, even in this single-lithology environment. We interpret that this complexity arises from varying fracture mechanisms (e.g. dilation-inducing fractures vs compression-inducing fractures) in different parts of the core, due to differences in stress conditions. Additional complexity may result from chemical interactions between the fresh fractures and the fluids in the area. Water content appears to have a significant, if not dominant, role in the unconfined compressive strength (UCS) of the samples. We suggest caution when using elastic property measurements as a proxy for fracturing in areas of explosion-induced damage, or in other areas where a variety of mechanisms induce fracturing.
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.
To increase understanding of damage associated with underground explosions, a field test program was developed jointly by Sandia and Pacific Northwest National Laboratories at the EMRTC test range in Socorro, NM. The Blue Canyon Dome test site is underlain by a rhyolite that is fractured in places. The test system included deployment of a defined array of 64 probes in eight monitoring boreholes. The monitoring boreholes radially surround a central near vertical shot hole at horizontal distances of 4.6m and 7.6m in cardinal and 45 degrees offset to cardinal directions, respectively. The probes are potted in coarse sand which touches/accesses the rhyolite and are individually accessed via nylon tubing and isolated from each other by epoxy and grout sequences. Pre and post chemical explosion air flow rate measurements, conducted for ~30-45 minutes from each probe, were observed for potential change. The gas flow measurement is a function of the rock mass permeability near a probe. Much of the flow rate change is at depth station 8 (59.4m) and is in the SE quadrant. Flow rate changes are inferred to be caused by the chemical explosion which may have opened pre-existing fractures, fractured the rock and/or caused block displacements by rotations and translations. The air flow rate data acquired here may enable a relationship and/or calibration to rock damage to be developed.
A critical component of the Underground Nuclear Explosion Signatures Experiment (UNESE) program is a realistic understanding of the post-detonation processes and changes in the environment that produce observable physical and radio-chemical signatures. Rock and fracture properties are essential parameters for modeling underground nuclear explosions. In response to the need for accurate simulations of physical and radio-chemical signatures, an experimental program to determine porosity, hydrostatic and triaxial compression, and Brazilian disc tension properties of P-Tunnel core was developed and executed. This report presents the results from the experimental program. Dry porosity for P-Tunnel core ranged from 8.7%-55%. Based on hydrostatic testing, bulk modulus was shown to increase with increasing confining pressure and ranged from 1.3GPa-42.3GPa. Compressional failure envelopes, derived from wet samples, are presented for P-Tunnel lithologies. Brazilian disc tension tests were conducted on wet samples and, along with triaxial tests, are compared with dry tests from the first UNESE test bed, Barnwell. P-Tunnel core disc tension test strength varied nearly two orders of magnitude between lithologies (0.03MPa-2.77MPa). Material tested in both tension and compression is weaker wet than dry with the exception of Strongly Welded Tuff in compression which is nearly identical in compressive strength for confining pressures of OMPa and 1 OOMPa. In addition to the inherent material properties of the rocks, fractures within the samples were quantified and characterized, in order to identify differences that might be caused by the explosion-induced damage. Finally, material property determinations are linked to optical microscopy observations. The work presented here is part of a broader material characterization effort; reports are referenced within.
Two blocks of alluvium were extensively tested at the Sandia National Laboratories Geomechanics laboratory. The alluvium blocks are intended to serve as surrogate material for mechanical property determinations to support the SPE DAG experimental series. From constant mean stress triaxial testing, strength failure envelopes were parameterized and are presented for each block. Modulus and stress relationships are given including bulk modulus versus mean stress, shear modulus versus shear stress, Young's modulus versus axial stress and Poisson's ratio versus axial stress. In addition, P-&S-wave velocities, and porosity, determined using helium porosimetry, were obtained on each block. Generally, both Young's modulus and Poisson's ratio increase with increasing axial stress, bulk modulus increases with increasing pressure, and increases more dramatically upon pore crush, shear modulus decreases with increasing shear stress and then appears to plateau. The Unconfined Compressive Strength for the BM is in the range of 0.5-0.6, and for SM in the range of 2.0-2.6 MPa. The confined compressive strength increases with increasing confining pressure, and the BM alluvium is significantly weaker compared to SM alluvium for mean stress levels above 8 MPa.
The transport properties of porous geological media are of fundamental importance when modeling the migration of chemical and radiological species in subterranean systems. Due to their relatively high mobility, short-lived noble gas species are of particular interest as detection of these species at the surface is a tell-tale indicator of recent nuclear activity. However, determining the diffusivity of these species is challenging due to their inert and transparent nature, requiring chemically insensitive techniques, such as mass spectroscopy, to quantify noble gas concentrations. The work described herein details recent advances in the methodology for determining diffusivity on porous media and results obtained on samples relevant to the UNESE project.
Understanding the dynamic behavior of geomaterials is critical for refining modeling and simulation of applications that involve impacts or explosions. Obtaining material properties of geomaterials is challenging, particularly in tension, due to the brittle and low-strength nature of such materials. Dynamic split tension technique (also called dynamic Brazilian test) has been employed in recent decades to determine the dynamic tensile strength of geomaterials. This is primarily because the split tension method is relatively straightforward to implement in a Kolsky compression bar. Typically, investigators use the peak load reached by the specimen to calculate the tensile strength of the specimen material, which is valid when the specimen is compressed at quasi-static strain rate. However, the same assumption cannot be safely made at dynamic strain rates due to wave propagation effects. In this study, the dynamic split tension (or Brazilian) test technique is revisited. High-speed cameras and digital image correlation (DIC) were used to image the failure of the Brazilian-disk specimen to discover when the first crack occurred relative to the measured peak load during the experiment. Differences of first crack location and time on either side of the sample were compared. The strain rate when the first crack is initiated was also compared to the traditional estimation method of strain rate using the specimen stress history.
Understanding the dynamic behavior of geomaterials is critical for refining modeling and simulation of applications that involve impacts or explosions. Obtaining material properties of geomaterials is challenging, particularly in tension, due to the brittle and low-strength nature of such materials. Dynamic split tension technique (also called dynamic Brazilian test) has been employed in recent decades to determine the dynamic tensile strength of geomaterials. This is primarily because the split tension method is relatively straightforward to implement in a Kolsky compression bar. Typically, investigators use the peak load reached by the specimen to calculate the tensile strength of the specimen material, which is valid when the specimen is compressed at quasi-static strain rate. However, the same assumption cannot be safely made at dynamic strain rates due to wave propagation effects. In this study, the dynamic split tension (or Brazilian) test technique is revisited. High-speed cameras and digital image correlation (DIC) were used to image the failure of the Brazilian-disk specimen to discover when the first crack occurred relative to the measured peak load during the experiment. Differences of first crack location and time on either side of the sample were compared. The strain rate when the first crack is initiated was also compared to the traditional estimation method of strain rate using the specimen stress history.
The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This report summarizes unconfined compression testing (UCS) from coreholes U-15n#12 and U-15n#13 and compares those datasets to UCS results from coreholes U-15n and U-15n#10. U-15n#12 corehole was drilled at -60° to the horizontal and U-15n#13 was drilled vertically in granitic rock (quartz monzonite) after the third SPE shot. Figure 1 illustrates at the surface, U 15n#12 and U-15n#13 coreholes were approximately 30 meters and 10 meters from the central SPE borehole (U-15n) respectively. Corehole U-15n#12 intersects the central SPE borehole (U 15n) at a core depth of 174 feet (approximately 150 feet vertical depth). The location of U 15n#12 and U-15n#13 is the site of the first, second and third SPE's, in Area 15 of the NNSS.
Dynamic Brazilian tension (DBR) tests from core hole U-15n are part of a larger material characterization effort for the Source Physics Experiment (SPE) project. This larger effort encompasses characterizing Climax Stock granite rock from the Nevada National Security Site (NNSS) both before and after each SPE shot. The current test series includes DBR tests on dry intact granite and fault material at depths of -85 and -150 ft.
Triaxial compression tests from core hole U-15n are part of a larger material characterization effort for the Source Physics Experiment (SPE) project. This larger effort encompasses characterizing Climax Stock granite rock from the Nevada National Security Site (NNSS) both before and after each SPE shot. The current test series includes triaxial compression tests on dry and saturated intact granite and fault material at 100, 200, 300, and 400 MPa confining pressure.
The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This memorandum reports on an initial phase of unconfined compression testing from corehole U-15n#10. Specimens tested came from the U-15n#10 core hole, which was drilled at -60° to the horizontal in granitic rock (quartz monzonite) after the second SPE shot (SPE-2). Figure 1 illustrates at the surface, the core hole was approximately 90 feet from the central SPE borehole. Corehole U 15n#10 intersects the central SPE borehole (U-15n) at a core depth of 170 feet (approximately 150 feet vertical depth) which is within the highly damaged zone of SPE-2. The U-15n#10 location is the site of the first, second and third SPE's, in Area 15 of the NNSS.
The Nevada National Security Site (NNSS) will serve as the geologic setting for a Source Physics Experiment (SPE) program. The SPE will provide ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central bore hole at varying depths and recording ground motions in instrumented bore holes located in two rings around the source positioned at different radii. Modeling using advanced simulation codes will be performed both a priori and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented bore holes including the geomechanical behavior of the site rock/structural features. This report presents a limited scope of work for an initial phase of primarily unconfined compression testing. Samples tested came from the U-15n core hole, which was drilled in granitic rock (quartz monzonite). The core hole was drilled at the location of the central SPE borehole, and thus represents material in which the explosive charges will be detonated. The U-15n location is the site of the first SPE, in Area 15 of the NNSS.
A critical component of the Underground Nuclear Explosion Signatures Experiment (UNESE) program is a realistic understanding of the post-detonation processes and changes in the environment that produce observable physical and radio-chemical signatures. Rock and fracture properties are essential parameters for any UNESE test bed. In response to the need for accurate modeling scenarios of these observations, an experimental program to determine the permeability and direct shear fracture properties of Barnwell core was developed. Room temperature gas permeability measurements of Barnwell core dried at 50degC yield permeability ranging from 6.24E-02 Darcys to 6.98E-08 Darcys. Friction angles from the direct shear tests vary from 28.1deg to 44.4deg for residual shear strength and average 47.9deg for peak shear strength. Cohesion averaged 3.2 psi and 13.3 psi for residual and peak shear strength values respectively. The work presented herein is the initial determination of an ongoing broader material characterization effort.
We characterize geomechanical constitutive behavior of reservoir sandstones at conditions simulating the “Cranfield” Southeast Regional Carbon Sequestration Partnership injection program. From two cores of Lower Tuscaloosa Formation, three sandstone lithofacies were identified for mechanical testing based on permeability and lithology. These include: chlorite-cemented conglomeratic sandstone (Facies A); quartz-cemented fine sandstone (Facies B); and quartz- and calcite-cemented very fine sandstone (Facies C). We performed a suite of compression tests for each lithofacies at 100 °C and pore pressure of 30 MPa, including hydrostatic compression and triaxial tests at several confining pressures. Plugs were saturated with supercritical CO2-saturated brine. Chemical environment affected the mechanical response of all three lithofacies, which experience initial plastic yielding at stresses far below estimated in situ stress. Measured elastic moduli degradation defines a secondary yield surface coinciding with in situ stress for Facies B and C. Facies A shows measurable volumetric creep strain and a failure envelope below estimates of in situ stress, linked to damage of chlorite cements by acidic pore solutions. The substantial weakening of a particular lithofacies by CO2 demonstrates a possible chemical-mechanical coupling during injection at Cranfield with implications for CO2 injection, reservoir permeability stimulation, and enhanced oil recovery.