Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pore pressure”

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.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measuring Fault Zone and Host Rock Hydraulic Properties Using Tidal Responses

Pore pressure diffusion along faults influences induced seismicity and rupture mechanics. In situ hydraulic diffusivity measurements along faults are rare and generally lower than inferred from seismicity migration. Here we use the tidal response of deep geothermal boreholes to measure fault diffusivity and permeability. Initial interpretations of the observation with a homogeneous confined aquifer model result in diffusivities of 10 -3 –10 -1 m 2 /s. However, this model mixes signals from both the conduit and the host rock. We develop a model for tidal response with a fault passing through the aquifer based on the fault-guided fracture network and solve for hydraulic properties in both the fault and the host rock. The resulting fault permeability is 2 × 10 -14 –7 × 10 -14 m 2 (90% confidence interval[CI]) and fault diffusivity is 0.08–0.33 m 2 /s (90% CI), which is 2 orders of magnitude higher than the host rock diffusivity in some wells, thus highlighting the role of faults as fluid conduits.

58 GEOSCIENCES↗
DOE OSTI2022

Monitoring Seasonal Shear Wave Velocity Changes in the Top 6 m at Garner Valley in Southern California With Borehole Data

Abstract Subsurface structures play important roles in seismic ground motion, crustal hydrology, stability of the built environment, and more. Constraining temporal changes of subsurface shear wave velocity ( V S ) can provide useful information to all these topics and the growing field of hydrological monitoring with seismic velocity. Using borehole records at Garner Valley, CA, we estimate seasonal subsurface V S variations from impulse response functions (IRFs) of earthquake data (2005–2018) along with IRFs and cross‐correlation of cross‐hole experiment data (2015–2018). The inferred V S variations are up to ∼25% in the top 6 m and ∼10% at 2–5 m in depth. The V S variations correlate strongly with the water table depth changes, suggesting that the changes are mostly due to fluctuations of pore pressure in the shallow material. The shallow velocity changes alter the near‐surface conditions, can affect seismic hazard estimation, and may be improperly attributed to deeper processes without careful analysis.

58 GEOSCIENCES↗
DOE OSTI2023

Understanding Subsurface Fracture Evolution Dynamics Using Time‐Lapse Full Waveform Inversion of Continuous Active‐Source Seismic Monitoring Data

Abstract Predicting the behavior, geometry, and flow properties of subsurface fractures remains a challenging problem. Seismic models that can characterize fractures usually suffer from low spatiotemporal resolution. Here, we develop a correlative double‐difference time‐lapse full waveform inversion of continuous active source seismic monitoring data for determining high‐spatiotemporal‐resolution time‐lapse V p models of in‐situ fracture evolution at a shallow contamination site in Wyoming, USA. Assisted by rock physics modeling, we find that (a) rapidly increasing pore pressure initializes and grows the fracture, increasing the porosity slightly (from ∼13.7% to ∼14.6%) in the tight clay formation, thus decreasing V p (∼50 m/s); (b) the fluid injection continues decreasing V p , likely through the introduction of gas bubbles in the injectate; and (c) final V p reductions reach over ∼150 m/s due to a posited ∼4.5% gas saturation. Our results demonstrate that high‐resolution V p changes are indicative of mechanical and fluid changes within the fracture zone during hydrofracturing.

58 GEOSCIENCES↗
DOE OSTI2022

Influence of Fault Architecture on Induced Earthquake Sequence Evolution Revealed by High-Resolution Focal Mechanism Solutions

The increasing seismicity and improved seismic observation network in recent years provide an opportunity to explore factors that influence the triggering processes, spatiotemporal evolution, and maximum magnitude of induced sequences. We map the fault architecture and stress state of four induced sequences in Oklahoma to determine their influence on the seismicity. We systematically relocate the earthquakes and compute hundreds of focal mechanisms of small to medium events (1.0 < M < 5.1) using various techniques, including machine learning, for the Guthrie, Woodward, Cushing, and Fairview sequences. The detailed fault geometry and spatiotemporal evolution of seismicity and stress states reveal different dominant driving forces for each sequence. In Cushing and Fairview (largest event ≥M5.0), the main fault structures are near-vertical narrow strike-slip faults, with most of the small earthquake fault planes optimally oriented. The two sequences exhibit discontinuous temporal migration but strong earthquake self-driven rupture growth. In Guthrie and Woodward (largest event <M5.0), the two sequences show more complex diffuse fault structures with varying dipping angles along depth. The inverted focal mechanisms show a mix of strike-slip faulting and normal faulting in both sequences, and the normal faulting events are less optimally oriented than strike-slip events. The two sequences are dominated by continuous diffusive migration in time driven by pore pressure propagation. The above results suggest that fault architecture and stress state influence sequence evolution, major driving forces, and possibly maximum magnitude.

58 GEOSCIENCES↗
DOE OSTI2023

Joint Geophysical and Numerical Insights of the Coupled Thermal‐Hydro‐Mechanical Processes During Heating in Salt

Abstract Salt offers an optimal medium for the permanent isolation of heat‐producing radioactive waste due to its impermeability, high thermal conductivity, and ability to close fractures through creep. A thorough understanding of the thermal‐hydrological‐mechanical (THM) processes, encompassing brine migration, is fundamental for secure radioactive waste disposal within salt formations. At the Waste Isolation Pilot Plant (WIPP), we conducted joint in situ geophysical monitoring experiments during active heating to investigate brine migration near excavations. This experiment incorporated electrical resistivity tomography (ERT) alongside high‐resolution fiber‐optic‐based distributed temperature sensing within a controlled heating experiment. Additionally, discrete element model (DEM) based numerical simulations were conducted to simulate THM processes during heating, providing a more mechanistic understanding of the coupled processes leading to the observed changes in geophysical measurements. During heating, resistivity shifts near the heater were reasonably explained by temperature effects. However, in more distant, cooler regions, the resistivity decrease exceeded predictions based solely on temperature. DEM simulations highlighted brine migration, propelled by pore pressure gradients, as the likely primary factor contributing to the additional resistivity decline beyond temperature effects. The comparison between the predicted ERT responses and observations was much improved when considering the effects of brine migration based on the DEM simulations. These geophysical and simulation findings shed light on brine migration in response to salt heating, enhancing our understanding of the coupled THM processes in salt for safe radioactive waste disposal.

58 GEOSCIENCES↗
DOE OSTI2024

Relating Hydro–Mechanical and Elastodynamic Properties of Dynamically Stressed Tensile–Fractured Rock in Relation to Applied Normal Stress, Fracture Aperture, and Contact Area

We exploit nonlinear elastodynamic properties of fractured rock to probe the micro-scale mechanics of fractures and understand the relation between fluid transport and fracture aperture under dynamic stressing. Experiments were conducted on rough, tensile-fractured Westerly granite subject to triaxial stresses. We measure fracture permeability for steady-state fluid flow with deionized water. Pore pressure oscillations are applied at amplitudes ranging from 0.2 to 1 MPa at 1 Hz frequency. During dynamic stressing we transmit ultrasonic signals through the fracture using an array of piezoelectric transducers (PZTs) to monitor evolution of interface properties. We examine the influence of fracture aperture and contact area by conducting measurements at effective normal stresses of 10–20 MPa. Additionally, the evolution of contact area with stress is characterized using pressure sensitive film. These experiments are conducted separately with the same fracture and map contact area at stresses from 9 to 21 MPa. The measurements are a proxy for “true” contact area for the fracture surface and we relate them to elastic properties using the calculated PZT sensor footprints via numerical modeling of Fresnel zones. We compare the elastodynamic response of the fracture using the stress-induced changes in ultrasonic wave velocities for transmitter-receiver pairs to image spatial variations in contact properties. We show that nonlinear elasticity and permeability enhancement decrease with increasing normal stress. Additionally, post-oscillation wave velocity and permeability exhibit quick recoveries toward pre-oscillation values. Estimates of fracture contact area (global and local) demonstrate that the elastodynamic and permeability responses are dominated by fracture topology.

58 GEOSCIENCES↗
DOE OSTI2021

Photoporomechanics: An Experimental Method to Visualize the Effective Stress Field in Fluid-Filled Granular Media

Effective stress governs the mechanical behavior of porous media. In this study, we use photoelasticimetry to visualize the evolving effective stress field in fluid-filled granular media in processes that couple fluid flow and mechanical deformation. We refer to this experimental method as photoporomechanics. We develop a fabrication process to produce millimeter-scale residual-stress-free photoelastic spheres with high geometric accuracy. We use color to quantify the forces acting on the particles over a wide range of forces, while using light intensity for a small range of forces. We then provide an application of photoporomechanics to illustrate the evolution of effective stress during one-dimensional consolidation: a process by which the stresses caused by a sudden load are gradually transmitted through a fluid-filled granular pack as the fluid drains and excess pore pressures dissipate. We report our technique provides a powerful experimental model system to study the grain-scale underpinning of coupled solid-fluid processes in granular media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Environmental cell for in situ X-ray synchrotron micro-CT imaging with simultaneous acoustic measurements

Synchrotron radiation provides the necessary spatial and temporal resolution for non-invasive operando studies of dynamic processes under complex environmental conditions. Here a new environmental cell for simultaneous in situ dynamic X-ray imaging and measuring acoustic properties of geological samples is presented. The primary purpose of this cell is to study gas-hydrate formation in porous geo-materials and its influence on their acoustic properties. The cell is designed for cylindrical samples of 9 mm in diameter, confining and pore pressures up to 12 MPa, and temperatures from –20°C to room temperature. The cell is portable and can be easily assembled and operated at different X-ray sources. This cell enables a wide range of experiments studying physical/chemical processes in the Earth subsurface that change the mechanical properties of rocks (geochemical reactions, phase transitions, etc.).

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Surface Deformation and Seismicity Linked to Fluid Injection in the Raton Basin

Abstract It is suggested that in addition to seismicity deep fluid injection may cause surface uplift and subsidence in oil and gas‐producing regions. This study uses the Raton Basin as an example to investigate the hydromechanical processes of surface uplift and subsidence during wastewater injection. The Raton Basin, in southern central Colorado and northern central New Mexico, has experienced wastewater injection related to coalbed methane and gas production starting in 1994. In this study, we estimate the extent and magnitude of total vertical deformation in the Raton Basin from 1994 to 2020 and incremental deformation between the years 2017 to 2020. Results indicate a modeled uplift as much as 15 cm occurring between 1994 and 2020. Between 2017 and 2020, up to 3 cm of uplift occurred, largely near the northwestern injection wells. Most modeled uplift between 1994 and 2020 occurred near the southern wells, where the greatest cumulative volume of wastewater was injected. However, modeled subsidence occurred around the southern and eastern wells between 2017 and 2020, after the rate of injection decreased. Modeling indicates that while the magnitude of modeled uplift corresponds to cumulative injection volume and maximum rate in the long‐term, short‐term incremental deformation (uplift or subsidence) is controlled by changes in the rate of injection. The number of yearly earthquake events follows periods of rapid modeled uplifting throughout the Basin, suggesting that measurable surface deformation may be caused by the same injection‐induced pore pressure perturbations that initiate seismicity.

Geology↗
DOE OSTI2023

Utah FORGE: Fault Reactivation Through Fluid Injection Induced Seismicity Laboratory Experiments

Included are results from shear reactivation experiments on laboratory faults pre-loaded close to failure and reactivated by the injection of fluid into the fault. The sample comprises a single-inclined-fracture (SIF) transecting a cylindrical sample of Westerly granite. All experiments are conducted at ambient temperature and follow a similar protocol: (i) application of confining stresses (3MPa) on the fault fully saturated with DI water, (ii) shear-mobilization through the increase of axial loading at a constant displacement rate until a post-peak steady-state condition is reached, (iii) reduction of axial loading and related shear stress to a prescribed fraction of the peak steady-state frictional strength (typically 60% to 90%, representing intermediate to high magnitudes) and (iv), fault reactivation triggered by a stepwise increase of pore pressure on the fault in 0.1 MPa increments held constant for 1-5 minutes. Mechanical data from three ISCO pumps connected to a Temco pressure vessel measure axial, confining, and fault-related parameters, including fluid pressure (kPa), fluid flow rate (mL/min), and axial displacement (mm). See included code for initial data analysis and visualization for select experiments. Resource names represent experiment numbers found in the "Read Me" file, which describes each experimental setup and parameters.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Monitoring of ground water table depth and soil moisture at the Point Reyes field site

Ground water table (GWT) depth and soil moisture (SM) have been monitored at several locations at the Point Reyes field site (Californian coastal grassland) from 2021 to 2024. Monitoring is still on-going and data may be added to this archive at later time. The SM data have been acquired using Teros 12 Meter soil moisture sensors placed at 10, 30, 60 and 90 cm depth at 5 locations along a small hillslope. These sensors also collect soil temperature and bulk conductance. In addition, some collocated sensors provide pore pressure and Photochemical Reflectance Index (PRI). The GWT depth has been inferred from various type of Onset pressure transducers. The pressure measurements have been corrected for atmospheric pressure variations and sensor position relative to the ground surface to infer GWT depth, as well as with RTK GPS data to infer GWT elevation. The GWT data have been acquired at 5 distinct locations from 2020 to 2024 with the sensors placed at about 4 m depth. In addition, GWT data has been acquired for the 2023-2024 period with sensors located in 1 m deep shallow wells installed near each deeper well. This data is intended to evaluate possibly different dynamic in shallow (perched) and deep aquifer. The datasets are all provided in csv format. Please note that the interpretation of the GWT data needs to be done with consideration of environmental and well characteristics at the site and uncertainty in various variables. For more information on GWT and SM data, please contact the author.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hydraulic fracture characterization by integrating multidisciplinary data from the Hydraulic Fracturing Test Site 2 (HFTS-2)

Various technologies have traditionally been used to monitor and describe hydraulic fractures from different perspectives. This work demonstrates the value of data integration for hydraulic fracture characterization when multiple data resources are available. The Hydraulic Fracturing Test Site 2 (HFTS 2) is a hydraulic fracturing research project in the Delaware Basin with multiple surveillance techniques including fiber optics sensing, microseismic, pressure/temperature gauges, etc. We integrated the multidisciplinary data from the HFTS-2 to characterize hydraulic fractures. The integrated data revealed interesting fracture propagation features including layering, vertical propagation affected by pore pressure gradient, and different microseismic activities due to difference in-situ conditions. Furthermore, these findings can be insightful for understanding hydraulic fracture propagation. The comparison among multiple surveillance data also helps us to evaluate the roles of various surveillance technologies and provides us experience to make informative decisions depending on different monitoring objectives.

58 GEOSCIENCES↗
DOE OSTI2021

Wastewater Disposal Has Not Significantly Altered the Regional Stress State in Southern Kansas

Wastewater disposal is primarily responsible for the increased seismicity rate since ~2013 in southern Kansas. Previous work that used shear-wave splitting (SWS) in southern Kansas interpreted an ~90° temporal rotation in the fast polarization direction and attributed it to increased pore pressures resulting from fluid injection. However, this interpreted rotation coincided with a change in the stations used to make the SWS measurements. We investigate the temporal variability of fast azimuths in southern Kansas by making SWS measurements on earthquake families with similar source–receiver paths recorded on a stable local seismic network. We select high-quality SWS measurements by investigating the stability of results across 65 different frequency bands between 0.5 and 15 Hz. We find that the fast polarization direction in southern Kansas is relatively constant with an average east-northeast (~N79°E) orientation between 2014 and 2017. Our fast polarization measurements are primarily a reflection of the maximum principal horizontal stress direction (SHmax). We observe a slight spatial change in SHmax to the northeast (~N55°E) near the Nemaha ridge in Oklahoma. However, we do not observe any significant temporal rotation of SHmax or variation in delay time (i.e., crack density) in southern Kansas, contrary to the earlier study. The previously interpreted ~90° rotation may either be a reflection of a very local stress change or a misinterpretation of SWS results potentially due to the use of inconsistent source–receiver paths. Our SWS measurements cover the period of peak wastewater disposal and seismicity rates and suggest an absence of significant temporal rotations in the local anisotropy and stress orientations associated with wastewater disposal.

Geochemistry & Geophysics↗
DOE OSTI2023

Modeling-Based Assessment of Deep Seismic Potential Induced by Geologic Carbon Storage

Induced seismicity is an inherent risk associated with geologic carbon storage (GCS) in deep rock formations that could contain undetected faults prone to failure. Modeling-based risk assessment has been implemented to quantify the potential of injection-induced seismicity, but typically simplified multiscale geologic features or neglected multiphysics coupled mechanisms because of the uncertainty in field data and computational cost of field-scale simulations, which may limit the reliable prediction of seismic hazard caused by industrial-scale CO 2 storage. The degree of lateral continuity of the stratigraphic interbedding below the reservoir and depth-dependent fault permeability can enhance or inhibit pore-pressure diffusion and corresponding poroelastic stressing along a basement fault. In conclusion, this study presents a rigorous modeling scheme with optimal geological and operational parameters needed to be considered in seismic monitoring and mitigation strategies for safe GCS.

58 GEOSCIENCES↗
DOE OSTI2024

Cyclic Injection Leads to Larger and More Frequent Induced Earthquakes under Volume-Controlled Conditions

As carbon storage technologies advance globally, methods to understand and mitigate induced earthquakes become increasingly important. Although the physical processes that relate increased subsurface pore pressure changes to induced earthquakes have long been known, reliable methods to forecast and control induced seismic sequences remain elusive. Suggested reservoir engineering scenarios for mitigating induced earthquakes typically involve modulation of the injection rate. Some operators have implemented periodic shutdowns (i.e., effective cycling of injection rates) to allow reservoir pressures to equilibrate (e.g., Paradox Valley) or shut-in wells after the occurrence of an event of concern (e.g., Basel, Switzerland). Other proposed scenarios include altering injection rates, actively managing pressures through coproduction of fluids, and preinjection brine extraction. Here, in this work, we use 3D physics-based earthquake simulations to understand the effects of different injection scenarios on induced earthquake rates, maximum event magnitudes, and postinjection seismicity. For comparability, the modeled injection considers the same cumulative volume over the project’s operational life but varies the schedule and rates of fluid injected. Simulation results show that cyclic injection leads to more frequent and larger events than constant injection. Furthermore, with intermittent injection scenario, a significant number of events are shown to occur during pauses in injection, and the seismicity rate remains elevated for longer into the postinjection phase compared to the constant injection scenario.

58 GEOSCIENCES↗
DOE OSTI2023

Discrete Measurements of the Least Horizontal Principal Stress from Core Data: An Application of Viscoelastic Stress Relaxation

Summary The emerging Paradox Oil Play in southeastern Utah is among the most significant unconventional plays in the western USA. The mean total undiscovered oil resources within just the Pennsylvanian Cane Creek interval of the Paradox Basin are believed to exceed 215 million barrels. However, to date, less than 5% (~9 million barrels) of the total Cane Creek resource has been produced from fewer than 40 wells, and only approximately one-half of those are horizontal wells. More than 95% of production is from the central Cane Creek Unit (CCU). Natural fractures are a key feature of many production wells, but stimulation by induced hydraulic fractures is not consistently successful. We hypothesize that more effective production in this play will rely on better fundamental characterization, especially on better quantification of the state of stress. Approximately 110 ft of core, well logs, and a diagnostic fracture injection test (DFIT) were acquired from the State 16-2 well within the CCU. With these data, we applied two methods to constrain and clarify the state of stress. The first technique, the Simpson’s coefficient method, provides lower bounds on the two horizontal principal stresses and relies on only limited data. Alternatively, the viscoelastic stress relaxation (VSR) method is used to estimate the least horizontal principal stress, building on observations that principal stresses become more isotropic as the viscous behavior of a rock is more pronounced. Results of these two methods support the hypothesis that the state of stress in the CCU of the Paradox Basin is nearly lithostatic and isotropic. Other factors consistent with this hypothesis include high formation pore pressure, which tends to reduce the possible stress states by changing the frictional failure equilibrium; lack of induced fractures in the core, which should be present in the case of stress anisotropy; and interbedded halite layers, which given their high degree of ductility, probably lead to greater VSR for the entire sedimentary package.

Energy & Fuels↗
DOE OSTI2023

M4 Summary of EBS International Activity

Thermal-Hydrologic (TH) modeling of DECOVALEX 2023, Task C has continued in FY23. This report summarizes progress in TH modeling of Step 1c, with calibration modeling and the addition of shotcrete. The work involves 3-D modeling of the full-scale emplacement experiment at the Mont Terri Underground Rock Laboratory (Nagra, 2019). While Step 1 is focused on modeling the heating phase of the FE experiment with changes in pore pressure in the Opalinus clay resulting from heating, Step 1c is focused on calibration of models using available data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

MRCI Subtask 2.2: Precambrian Basement Structure, Faulting, and Stress Final Technical Summary Report

Subtask 2.2 was largely focused on compiling information about the Precambrian basement rock that underlies the MRCI region. Understanding the nature and extent of the basement rock is important in order to mitigate certain risks associated with CO2 injection. For example, this information will provide a resource for defining areas where a storage reservoir is in direct contact with the Precambrian or indirectly via a fault. Also, an understanding of the stress field in addition to fault orientations will help determine which faults may be subject to reactivation under increased pore pressures.

Battelle,CCS,CCUS,Carbon storage,MRCI,Midwest Regi↗