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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗
NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗
DOE OSTI2022

Utah FORGE Hydrothermal Friction-Hydraulic Transmissivity Laboratory Experiments

Submission includes data from laboratory slide-hold-slide tests, combined with flow through tests, conducted on Westerly granite with 30 degree sawcut. Tests were conducted with a constant confining pressure of 30 MPa with an average pore pressure of 10 MPa at temperatures of 23 and 200 degC. Three fluid flow conditions were examined (1) no flow, (2) cycled flow, and (3) continuous flow. Data were collected to asses the effect of temperature and pore fluid on frictional healing rates in granite at geothermal conditions. Data is available in XML and JSON data types.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Cement stress and microstructure evolution during curing in semi-rigid high-pressure environments

The process of cement slurry curing into a solid has a direct effect on the initial state of stress in a composite concrete and steel structure. Failure models for critical infrastructure and wellbores could benefit from measurements of this initial stress because it greatly influences structural performance. Here, we present measurements of the initial stress state (total stress, effective stress, and pore pressure) of multiple cement formulations cured under a high-pressure of 40 MPa within a semi-rigid steel pipe. Our results show that cements cured in this restrained environment exhibit dense microstructures with low permeability (nano-Darcy), develop an anisotropic stress state greater than zero and less than the curing pressure, and undergo a staged curing process where hydration and mechanical behaviors evolve from a slurry into a pre-stressed porous medium. Both drained and undrained scenarios were investigated, revealing a feedback from the semi-rigid boundary towards improved cement performance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Utah FORGE: Fluid Injection Induced Shearing Experiments on Fractured Granitoid at Elevated Temperatures

This repository contains experimental data from a series of fluid injection-induced shearing tests conducted on Utah FORGE granitoid. The experiments were performed using an aluminum triaxial pressure vessel (TEMCO) apparatus at Pennsylvania State University. The primary aim was to investigate the impact of temperature on fault seismicity and to explore how different pre-stress ratios influence shearing behavior. The data includes results from experiments conducted under varying conditions of temperature, pre-stress ratios, and pore pressure increments. The rock samples used in these experiments were 60-grit granitoid with a single inclined fracture oriented at 60 degrees with respect to the horizontal cross-section. The dataset also includes measurements of normal stiffness at different temperatures and images of the experimental apparatus. The included README file details each experiential setup and outlines which data files represent which experimental conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗
DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2021

Predicting Fluid Flow Regime, Permeability, and Diffusivity in Mudrocks from Multiscale Pore Characterisation

In geoenergy applications, mudrocks prevent fluids to leak from temporary (H 2 , CH 4 ) or permanent (CO 2 , radioactive waste) storage/disposal sites and serve as a source and reservoir for unconventional oil and gas. Understanding transport properties integrated with dominant fluid flow mechanisms in mudrocks is essential to better predict the performance of mudrocks within these applications. In this study, small-angle neutron scattering (SANS) experiments were conducted on 71 samples from 13 different sets of mudrocks across the globe to capture the pore structure of nearly the full pore size spectrum (2 nm–5 μm). We develop fractal models to predict transport properties (permeability and diffusivity) based on the SANS-derived pore size distributions. The results indicate that transport phenomena in mudrocks are intrinsically pore size-dependent. Depending on hydrostatic pore pressures, transition flow develops in micropores, slip flow in meso- and macropores, and continuum flow in larger macropores. Fluid flow regimes progress towards larger pore sizes during reservoir depletion or smaller pore sizes during fluid storage, so when pressure is decreased or increased, respectively. Capturing the heterogeneity of mudrocks by considering fractal dimension and tortuosity fractal dimension for defined pore size ranges, fractal models integrate apparent permeability with slip flow, Darcy permeability with continuum flow, and gas diffusivity with diffusion flow in the matrix. This new model of pore size-dependent transport and integrated transport properties using fractal models yields a systematic approach that can also inform multiscale multi-physics models to better understand fluid flow and transport phenomena in mudrocks on the reservoir and basin scale.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Robust In-Situ Strain Measurements to Monitor CO 2 Storage

The goal of this project was to develop and demonstrate robust instrumentation to monitor the in-situ strain tensor in order to improve the reliability and security of CO 2 storage in geologic formations. We met the original goals of the project and the major overarching accomplishment is the advancement of strain tensor monitoring from an intriguing concept to a commercially available technology with a solid foundation of novel instruments supported by theoretical analyses and validation experiments. The main accomplishments of the project are summarized below. We designed, built and evaluated nine new optical fiber strainmeters and tiltmeters using Michelson interferometers to measure deformation with ultra-high resolution at both shallow and deep point locations in the subsurface. These are the robust strainmeters that motivated the title of the project. We designed, built and evaluated a novel method of measuring distributed strain in optical fibers with nanostrain resolution, and cm-scale location, and sampling into the seismic band. The new method is called Coherence-length-gated Microwave Photonics Interfereometry (CMPI). CMPI technology has advantages over existing commercial DAS and DSS methods. We developed and demonstrated capabilities to deploy instruments in the field and used them to measure strain caused by ambient signals like barometric pressure and tides, as well as induced signals like surface loading and pore pressure changes from pumping tests. We deployed a working strainmeter at 1,700 ft depth, slightly above an active reservoir. This is to our knowledge the greatest depth a strainmeter has been deployed and the techniques we used can readily be extended to greater depths. Optical fiber borehole tensor strainmeter techology was advanced from a TRL 4 at the start, to a TRL of 7 at the conclusion of the project. The project included advances in simulations and theoretical analyses. We developed and demonstrated a computational workflow that uses machine learning to reduce the computational requirements and make it practical to use Bayesian inversion to solve large numerical poroelastic analyses needed to interpret strain tensor field data. We evaluated the strain tensor fields and time series that would be caused by leaks of CO 2 or other fluids from reservoirs. These simulations demonstrated that signals from leaks could be measured with instruments developed for the project, opening a potentially new method for ensuring storage security. We showed that strains in caprock can be used to estimate pressure in a reservoir. This avoids the need to drill monitoring wells into the reservoir, and it expands the capabilities of monitoring in the caprock. The project includes a derivation and application of a novel analytical solution to the strains in the vicinity of a pressurized poroelastic inclusion. This solution explains field data measured during injeciton tests at the North Avant Field, and it will simplify future interpretation of strain tensor data. The project included a broad range of experiments, and of the most significant is the characterization of the strain tensor at an array three strainmeters during six injection tests at the North Avant Field, Oklahoma. This demonstrated repeatability of the strain signal measured by the new instruments developed for the project, and it showed similarities between the strain signal at shallow depths and pressure in the underlying reservoir. We also demonstrated that useful strain data can be measured at reservoir depths. This confirms that strain tensor data can be measured throughout the caprock over a reservoir. The project demonstrated the feasibility of using the strain tensor and distributed strain measured in caprock during a variety of different well tests where the pumping rate was constant, sinusoidal and positive, or a periodic square wave with zero net rate. This further strengthens the validity of using strain data to characterize reservoirs and aquifers. We also demonstrated that strain caused be fluctuations of air pressure and water pressure in the vadose zone can be measured and interpreted, suggesting that high resolution distributed strain measurements hold promise for monitoring the vadose zone. The project partially supported nine graduate students in the Environmental Engineering, Hydrogeology, Electrical Engineering programs at Clemson University. The research was described in nine journal papers, 23 talks and conference abstracts. Additional journal papers are in preparation. A new company called Tensora was started to provide strainmeter technology for commercial applications.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗
NASA NTRS1995

Finite Element Analysis of Poroelastic Composites Undergoing Thermal and Gas Diffusion

A theory for time-dependent thermal and gas diffusion in mechanically time-rate-independent anisotropic poroelastic composites has been developed. This theory advances previous work by the latter two authors by providing for critical transverse shear through a three-dimensional axisymmetric formulation and using it in a new hypothesis for determining the Biot fluid pressure-solid stress coupling factor. The derived governing equations couple material deformation with temperature and internal pore pressure and more strongly couple gas diffusion and heat transfer than the previous theory. Hence the theory accounts for the interactions between conductive heat transfer in the porous body and convective heat carried by the mass flux through the pores. The Bubnov Galerkin finite element method is applied to the governing equations to transform them into a semidiscrete finite element system. A numerical procedure is developed to solve the coupled equations in the space and time domains. The method is used to simulate two high temperature tests involving thermal-chemical decomposition of carbon-phenolic composites. In comparison with measured data, the results are accurate. Moreover unlike previous work, for a single set of poroelastic parameters, they are consistent with two measurements in a restrained thermal growth test.

Salamon, N. J.↗
DOE OSTI2024

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗
DOE OSTI2022

3D Simulations of gas injection on callovo-oxfordian claystone assuming spatial heterogeneity and anisotropy

A series of gas injection tests on Callovo-Oxfordian (COx) claystone from the Bure underground research laboratory (URL) in France were carried out at the British Geological Survey (BGS). The tests were performed using a triaxial apparatus specifically designed to capture small volumetric strains induced by the injected gas flow and consequent material dilatancy. The long-duration experiments were monitored throughout. Measurements also included pressure, stresses (axial and radial stresses prescribed for each test stage), rate of gas inflow, gas outflow volume as well as pore-pressures observed at various points of the sample. A coupled hydro-gas-mechanical 3D numerical model has been developed to simulate the tests. Initial permeability is assumed heterogeneous throughout the specimen and embedded fractures are incorporated in the formulation. Gas pressure-induced deformations during the test lead to variations of permeability due to changes in matrix porosity and, especially, fracture aperture as well as fracture orientation due to material anisotropy. A programme of sensitivity analyses involving the variation of different aspects and parameters of the model contributes to a better understanding of the phenomena and highlights its complexity. The model is able to reproduce the observed behaviour of the tests.

58 GEOSCIENCES↗
DOE OSTI2025

The Effect of Hydrothermal Alteration and Microcracks on Hydraulic Properties and Poroelastic Deformation: A Case Study of the Blue Mountain Geothermal Field

Abstract Geothermal energy plays a vital role in decarbonizing electricity and heat supply. Effective utilization of geothermal resources hinges on identifying or generating permeable reservoir zones and understanding how effective pressure variations affect fluid circulation and reservoir properties by poroelastic deformation. Hydrothermal alteration can modify the petrophysical properties of geothermal reservoir rocks, which may increase or decrease its productivity. Understanding these alteration effects is essential to predict and optimize long‐term sustainable geothermal operations. Here, we investigate the impact of hydrothermal alteration on poroelastic and hydraulic properties of diverse lithologies in a series of deformation tests performed at several confining (0–80 MPa) and pore pressure (10–30 MPa) levels. Experimental results of hydrothermally altered dikes and phyllites obtained from the Blue Mountain geothermal field (Nevada, USA) are compared to thermally cracked La Peyratte granite (France) and correlated with petrophysical properties, mineral composition, and microstructures. Argillic alteration of dikes increases porosity and storage capacity but lowers thermal conductivity and increases pore compressibility. Conversely, silicate precipitation in phyllites increases stiffness and thermal conductivity but also reduces porosity and permeability. Experimentally determined effective pressure coefficients range from 0.1 to 0.9, differ for permeability and volumetric strain and decrease with increasing effective pressure. The presence of compliant microcracks and crack‐like pores significantly increases the stress sensitivity of La Peyratte granite and silicified phyllites. This study demonstrates how thermal and chemical alteration impacts poromechanical and petrophysical characteristics of geothermal targets, which ultimately govern reservoir stability and subsidence, induced seismicity as well as fluid and heat extraction efficiency during geothermal operations.

Schuster, Valerian [Helmholtz Centre Potsdam GFZ G↗
NASA NTRS2021

Hot Isostatic Pressing for Reduction of Defects in Additively Manufactured Inconel-718 Preliminary Simulation Approach and Results

Selective laser melting (SLM) of Inconel-718 (IN718) is a promising advanced manufacturing technology; however, challenges exist for fracture-critical parts where pores induced during the manufacturing process reduce fatigue life of the material. Hot isostatic pressing (HIP) is a common post processing technique used to reduce the part porosity by exposing parts to high temperature and pressures. Pore closure depends on the pore shape, trapped gas, and HIP parameters, and the process can only be optimized by understanding the effect of these variables. An analysis capability is described to better understand the conditions necessary for pore closure, and the effect of several variables on porosity are explored and compared to a baseline case. The importance of valid material response models is identified as a significant factor in model accuracy; especially the creep material response model which, in the baseline case, predicts that creep strain accounts for 87.7% of the total strain at the conclusion of the HIP cycle. Lessons learned are described, in addition to future work including model validation.

Joshua M Fody↗
DOE OSTI2026

Thermal loading effects on chalk hydromechanical behavior for nuclear waste disposal

Safe disposal of heat-generating nuclear waste depends on host rock stability under thermal, hydrological, and mechanical stresses. This study investigates the effect of thermal loading on mechanical behavior of the shallowly buried Ghareb formation chalk through triaxial and hydrostatic constant strain rate and creep tests at temperatures up to 100 ˚C and effective pressures up to 20.7 MPa. Experimental results show that thermal loading reduces the elastic moduli of chalk by 50–75%, and a transition occurs above 60 ˚C where creep rates increase rapidly. Water saturation nearly doubles the thermally induced strain compared to dry conditions and strongly decreases material rigidity. Thermal loading also leads to significant pore pressure increases under undrained conditions and leads to reductions in the apparent permeability during drained conditions. Laboratory experimental data were used to parameterize and develop a preliminary constitutive model for predicting future deformation during repository operations in the Ghareb. The strongly coupled effects – mechanical weakening, fluid pressure fluctuations, and permeability modification – demonstrate that elevated repository temperatures will have a pronounced effect on the near field Ghareb behavior during waste disposal operations. The findings indicate that the coupled interactions must be considered in predictive models and repository design to ensure long-term nuclear waste isolation and safety.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1993

The seismic response of an aquifer to the propagation of an impact generated shock wave: A possible trigger of the Martian outflow channels?

Aquifer dilation from shock waves produced by the 8.4 magnitude Alaskan earthquake of 1964 led to water and sediment ejection from the ground up to 400 km away from the earthquake's epicenter. Groundwater disturbances were observed as far away as Perry, Florida (approximately 5500 km), where well water fluctuations with an amplitude of as much as 2.3 m were reported. The martian cratering record provides evidence that the planet has experienced numerous seismic events of a similar, and often much greater, magnitude. Given this fact, and the photogeologic evidence for abundant water in the early crust, the response of a basalt aquifer to the propagation of compressional waves (P-waves) produced by impacts in the 33-1000 km diameter size range were investigated. The resulting one-dimensional changes in effective stress and pore pressure were calculated - as a function of both distance and time - based on the following assumptions: (1) that all of the seismic energy radiated by an impact is transmitted as a single compressional wave; (2) that both the host rock and groundwater are compressible; and (3) that there is no net flow between the water-filled pores.

Leyva, Ivett A.↗
DOE OSTI2024

Modeling injection-induced fault slip using long short-term memory networks

Stress changes due to changes in fluid pressure and temperature in a faulted formation may lead to the opening/shearing of the fault. This can be due to subsurface (geo)engineering activities such as fluid injections and geologic disposal of nuclear waste. Such activities are expected to rise in the future making it necessary to assess their short- and long-term safety. Here, a new machine learning (ML) approach to model pore pressure and fault displacements in response to high-pressure fluid injection cycles is developed. The focus is on fault behavior near the injection borehole. To capture the temporal dependencies in the data, long short-term memory (LSTM) networks are utilized. To prevent error accumulation within the forecast window, four critical measures to train a robust LSTM model for predicting fault response are highlighted: (i) setting an appropriate value of LSTM lag, (ii) calibrating the LSTM cell dimension, (iii) learning rate reduction during weight optimization, and (iv) not adopting an independent injection cycle as a validation set. Several numerical experiments were conducted, which demonstrated that the ML model can capture peaks in pressure and associated fault displacement that accompany an increase in fluid injection. The model also captured the decay in pressure and displacement during the injection shut-in period. Further, the ability of an ML model to highlight key changes in fault hydromechanical activation processes was investigated, which shows that ML can be used to monitor risk of fault activation and leakage during high pressure fluid injections.

58 GEOSCIENCES↗