SEARCH · Search NASA
Results for “Subsurface Stress”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Initial Field Testing in the Deep Subsurface for the Thermal Breakout Project for Measuring In Situ Stress
The current state-of-the-art technology for in situ stress measurements involves an integrated approach that combines borehole breakout observations, drilling-induced tensile fractures, and hydraulic fracturing tests (i.e., mini-fracs); however, this methodology has several limitations that often prevent successful in situ stress measurements. One major limitation is that breakouts do not appear in every borehole and are generally a natural occurrence that cannot easily be controlled. Because borehole breakouts are used to measure the maximum horizontal in situ stress magnitude, the absence of borehole breakouts presents a major data gap for in situ stress measurements. As a response to this data gap a new thermal breakout technology is being developed that will provide a method for thermally inducing borehole breakouts and obtaining consistent measurements of the maximum horizontal stress magnitude. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which is directly correlated to the maximum horizontal stress magnitude. An initial component of the thermal breakout project involves developing a prototype downhole tool and performing field testing in the deep subsurface. The primary objective of the initial tool development and field testing is to provide a physical proof of concept for the thermal breakout technology. Three thermal breakout field tests were performed at approximately 4850 feet deep in the Sanford Underground Research Facility. Each test was performed in a borehole that had been previously used for overcoring stress measurements. The third and final test successfully created two diametrically opposed breakouts after deploying the thermal breakout tool for approximately 3 hours of continuous heating. Post-analysis of the thermally induced breakouts confirmed that the thermal breakout orientation directly corresponds to the known in situ stress orientation and magnitude.
Challenges in simulating ground interacting nuclear explosions
This paper summarizes recent above-ground nuclear explosion simulations as part of a broader effort to better characterize conditions within a fireball that may influence the chemical evolution of bomb materials and other materials entrained from the local explosion environment. A critical component of this work is validation against historic footage of atmospheric testing, requiring that we understand how the frequency-dependent sensitivity of the utilized film footage influences data captured in such images. We focus first on the early physics of a nuclear explosion in the atmosphere before discussing some of the technical challenges we seek to capture in late-time models that include more complex emplacement conditions and subsurface features. We discuss required physics packages (compressible hydrodynamics, radiation transport, as well as necessary ancillary tables such as equations of state (EOS) and opacities). Additionally, we note reasonable “shortcuts” one may make and their limitations, e.g., using ideal gas EOS, replacing spectrally resolved radiation with spectrally averaged radiation, and exchanging deterministic transport with diffusion. We then discuss an approach to achieving an equilibrated initial stress state for problems where buoyancy and subsurface lithostatic stress are important. Our methodology is presented in the context of LLNL’s ALE3D multiphysics code but may readily be implemented in other codes. In this paper, we start with a description of the challenges of NUDET simulations, followed by a presentation of the simulated intensity (flux) as it would appear on an analysis of the Dixie test. We then progressively introduce additional complexity in subsequent sections (near-surface burst and gravity initialization) before concluding.
Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir
Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.
Computational Analysis of Coupled Geoscience Processes in Fractured and Deformable Media
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.
Geologic stress modulates fluid mixing at fracture intersections
Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.
Integration of the Biot–Gassmann Fluid Substitution Method and Machine Learning-Based Velocity–Stress Relationship for Estimating In Situ Stresses
Recent advancements have shown that in situ stresses can be reliably estimated through an integrated machine/deep learning (ML/DL)-based framework, which relies on models trained and validated using true triaxial ultrasonic velocity (TUV) experimental data that involve measurements of ultrasonic velocity in saturated rocks under varying stress configurations. However, when the goal is to interpret lower frequency measurements, it may be more appropriate to run experiments on dry rocks and then obtain Biot–Gassmann-derived equivalent saturated velocities (low-frequency approximation) and employ these quantities for training ML/DL models to predict in situ stress. Whether the dispersion effect of frequency on the velocity–stress relationship substantially impacts in situ stress prediction is an important and unresolved question. This work presents an enhancement of ML/DL-based workflow by training and implementing ML/DL models using equivalent saturated acoustic velocities (low-frequency) obtained by applying Biot–Gassmann fluid substitution on the ultrasonic velocities of dry cores. The models were trained on TUV data sets derived from three subsurface cores extracted from the geothermal well 16B(78)-32 at the Utah FORGE site. Each core was subjected to 75 unique stress configurations for velocity measurement in the dry state. The ML/DL trained on the TUV data set with equivalent saturated velocities demonstrated promising performance to predict in situ stress in subsurface geological rocks using velocity–stress relationships with R 2 of 0.86, 0.971, and 0.975 and root mean squared error (RMSE) of 2.59, 1.92, and 1.80 for validation/testing phases of vertical, minimum horizontal, and maximum horizontal stress models, respectively. Additionally, interpretation and explanation by Shapley additive explanations (SHAP) analysis further improved scientific validation and model reliability for estimating in situ stresses.
Helium tribology of Inconel 617 subjected to laser peening for high temperature nuclear reactor applications
Inconel 617 is among the best candidates for utilization in high temperature gas cooled reactor tribo-components. However, the combined effects of sliding contact, along with intermittent idle times and very high temperature material degradation, deteriorates the alloy tribological performance, especially under a helium atmosphere. Laser peening is a surface treatment technique which can enhance the properties at the surface and subsurface by generating deep residual stresses and enhanced microstructure. In this work, we report the tribological behavior of regular laser peened as well as thermally-engineered laser peened Inconel 617 under helium and air atmospheres at 800 °C. In addition to friction and wear studies, the specimens are characterized by different analytical techniques to further understand the mechanisms involved in the peening process and sliding contact. Regardless of the peening process and post-process treatment types, it is observed that laser peening improves the tribological characteristics of Inconel 617. Interestingly, laser peening followed by helium thermal aging shows highly enhanced tribological behavior. This is attributed to the strengthening effect of the laser peening on the surface oxides providing an excellent and lasting protective and lubricating film under helium exposure.
Mechanical Response of Castlegate Sandstone under Hydrostatic Cyclic Loading
The stress history of rocks in the subsurface affects their mechanical and petrophysical properties. Rocks can often experience repeated cycles of loading and unloading due to fluid pressure fluctuations, which will lead to different mechanical behavior from static conditions. This is of importance for several geophysical and industrial applications, for example, wastewater injection and reservoir storage wells, which generate repeated stress perturbations. Laboratory experiments were conducted with Castlegate sandstone to observe the effects of different cyclic pressure loading conditions on a common reservoir analogue. Each sample was hydrostatically loaded in a triaxial cell to a low effective confining pressure, and either pore pressure or confining pressure was cycled at different rates over the course of a few weeks. Fluid permeability was measured during initial loading and periodically between stress cycles. Samples that undergo cyclic loading experience significantly more inelastic (nonrecoverable) strain compared to samples tested without cyclic hydrostatic loading. Permeability decreases rapidly for all tests during the first few days of testing, but the decrease and variability of permeability after this depend upon the loading conditions of each test. Cycling conditions do affect the mechanical behavior; the elastic moduli decrease with the increasing loading rate and stress cycling. The degree of volumetric strain induced by stress cycles is the major control on permeability change in the sandstones, with less compaction leading to more variation from measurement to measurement. The data indicate that cyclic loading degrades permeability and porosity more than static conditions over a similar period, but the petrophysical properties are dictated more by the hydrostatic loading rate rather than the total length of time stress cycling is imposed.
SOSAT: Geohazard Risk Assessment in Carbon Sequestration Operations
Like all operations in which fluids are injected into the subsurface, geologic carbon storage (GCS) presents inherit risks that must be assessed and mitigated to successfully deploy such technology in the field. For example, alterations to pore pressure and the stress state due to fluid injection may trigger certain geohazards, such as potential fault activation and induced seismicity or unintentional hydraulic fracturing. Here we present capabilities for assessing these potential risks using the State of Stress Analysis Tool (SOSAT): a Python library and web-based application capable of probabilistically estimating the subsurface state of stress informed by various field observations. Using posterior distributions of principal stress components, SOSAT can assess the risk of fault activation within the injection formation for either a critically oriented fault or a fault with a user-prescribed orientation. Additionally, SOSAT can estimate the probability of unintentional hydraulic fracturing of the intact reservoir rock. Here, we demonstrate these SOSAT capabilities using a hypothetical GCS site. This type of risk assessment can assist with the effective deployment of GCS technology in the field by informing safe design and management practices.
Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress
Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.
Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions
Understanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo–Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High–resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first–order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. Here, the results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.
Deep learning forecasts the spatiotemporal evolution of fluid-induced microearthquakes
Microearthquakes generated by subsurface fluid injection record the evolving stress state and permeability of reservoirs. Forecasting their spatiotemporal evolution is therefore critical for applications such as enhanced geothermal systems, carbon dioxide sequestration and other geoengineering applications. Here we propose a transformer neural network model that ingests hydraulic stimulation history and prior microearthquake observations to forecast four key quantities: cumulative microearthquake count, cumulative logarithmic seismic moment, and the 50th- and 95th-percentile extents of the microearthquake cloud. Applied to the EGS Collab Experiment 1 dataset, the model achieves R2 > 0.98 for the 1-s forecast horizon and R2 > 0.88 for the 15-s forecast horizon across all targets, and supplies uncertainty estimates through a learned standard deviation term. These accurate, uncertainty-quantified forecasts enable real-time inference of fracture propagation and permeability evolution, demonstrating the strong potential of deep-learning approaches to improve seismic-risk assessment and guide mitigation strategies in future fluid-injection operations.
Cooperative Research and Development Agreement between National Energy Technology Laboratory and VariGrid Explorations Inc. [Final Report]
The scope of this project is to advance the National Energy Technology Laboratory’s (NETL) Variable Grid Method (VGM) technology to commercial status thought enhancements upon prior work to produce a VGM application that can rapidly quantify multiple types of data and uncertainty commonly encountered by the subsurface energy industry. Along with development, the project will also include technology demonstrations, through applications using common industry data types and uncertainty metrics, to stress the value the VGM can offer the subsurface energy industry by providing an improved understanding of the strengths and limitations of their data, models, and analytical results. The resulting commercially ready VGM will offer interoperability with common industry software applications to ensure broad ‘plug and play’ functionality that will help keep costs relatively low and foster more rapid industry adoption.
3D hydrogeophysical characterization of managed aquifer recharge basins
Aquifers are increasingly stressed. Managed aquifer recharge provides a potential solution to mitigate this stress and provide sustainable groundwater resources. Subsurface properties are known to have a strong control on the infiltration rates that can be achieved. However, these properties are often highly heterogeneous and difficult to assess with conventional probing techniques. Here, we show the application of 3D geophysical imaging to assess the recharge potential and its variation across several basins used for managed aquifer recharge. We link in-situ measurements of saturated hydraulic conductivity with the electrical resistivity of the subsurface to establish petrophysical relationships and use those relationships to estimate the distribution of hydraulic conductivity throughout the five recharge basins. Our results show a considerable variability in the hydraulic properties, i.e., soil texture and saturated hydraulic conductivity, that have a direct impact on potential infiltration rates. We use the 3D hydraulic property distributions to model groundwater recharge and provide estimates for infiltration rates and volumes, and use this approach to assess the impact of management activities on groundwater recharge performance. Having such data not only enables us to predict infiltration rates, but also provides means for optimizing such water infrastructure.
Strengthening of Calcite Assemblages Through Chemical Complexation Reactions
Abstract Utilization of subsurface reservoirs, including fluid extraction/injection, can induce stress changes and modify in‐situ chemical equilibrium causing subcritical fracturing and deformation. Here, we show that chemical complexation reactions at the newly created fracture surfaces within the crack tip limit fracturing during consolidation of granular calcite assemblages. Previously, we showed that analogous chemical complexation reactions can increase in‐situ fracture toughness in pre‐fractured single crystal calcite. To test the chemical complexation effect under subsurface conditions, we consolidated samples under increasing hydrostatic pressure to induce widespread intergranular fracturing in the presence of fluid containing common anions. Measured consolidation, acoustic emissions, and microfracturing demonstrate that deformation correlates with complexation affinity for calcium in calcite and aqueous anions. We observed that deformation is lowest in the presence of strongly complexing anions (e.g., sulphate), provided sufficient anion diffusion rates. This understanding is important for predicting reservoir effective strength in the presence of complex aqueous fluids.
A Novel Rotating Flexure-Test Technique for Brittle Materials with Circular Geometries
The four-point bend test is one of the simplest and often the preferred flexural strength evaluation method for brittle materials. In this loading mode, fracture often initiates from a critical surface (or subsurface) flaw when subjected to a tensile stress state. However, if the critical flaw exists on the compression side of the test specimen, it may not activate to grow a crack and hence the resulting flexural strength will be higher than the true value. The goal of this study is to measure the true flexural strength of a solid or hollow cylindrical brittle specimen by ensuring that failure occurs at its weakest point by rotating along its longitudinal axis, thereby exposing and activating its critical surface flaw during a four-point bend test. A novel test fixture has been designed and fabricated, and the true flexural strengths of cylindrical brittle and quasi-brittle tubular specimens have been measured and compared to existing experimental data obtained through traditional four-point bend tests. In this work, experimental results showcase the orientation dependance on flexural strength for various materials. Additionally, similarities between experimental findings and those available in literature, including observations of fracture surfaces and relationships between surface roughness and material strength, are discussed.