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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗

DOE OSTI2025

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Effects of confinement and pressure on the structure and dynamics of carbon dioxide in silica slit pores

An understanding of carbon dioxide fluid properties within geological mesopores is important in applications ranging from carbon sequestration to shale oil recovery. Here, a molecular dynamics study is presented that aims to shed light on these systems by simulation of CO 2 fluid confined in β-cristobalite silica slit pores with different pore widths and pressure conditions. The weakly associating nature of carbon dioxide leads to little difference in structural and interfacial dynamical properties for different pore sizes and pressures. Rather, the behavior is found to be dominated by the entropic effects, namely, how the CO 2 organizes next to the silica surface. The CO 2 self-diffusion coefficient shows the strongest pore size dependence. It is strongly diminished in small pores and does not reach the bulk fluid value even in 6 nm pores. It also decreases with pressure, yielding an activation volume that increases with pore size. These results provide new insight into the behavior of a compressible fluid in nanoscale confinement.

Godahewa, Sahan M. [Univ. of Kansas, Lawrence, KS ↗

NASA NTRS2007

Micromechanics of Spray-On Foam Insulation

Understanding the thermo-mechanical response of the Space Shuttle External Tank spray-on foam insulation (SOFI) material is critical, to NASA's Return to Flight effort. This closed-cell rigid polymeric foam is used to insulate the metallic Space Shuttle External Tank, which is at cryogenic temperatures immediately prior to and during lift off. The shedding of the SOFI during ascent led to the loss of the Columbia, and eliminating/minimizing foam lass from the tank has become a priority for NASA as it seeks to resume scheduled space shuttle missions. Determining the nature of the SOFI material behavior in response to both thermal and mechanical loading plays an important role as any structural modeling of the shedding phenomenon k predicated on knowledge of the constitutive behavior of the foam. In this paper, the SOFI material has been analyzed using the High-Fidelity Generalized Method of Cells (HFGMC) micromechanics model, which has recently been extended to admit a triply-periodic 3-D repeating unit cell (RUC). Additional theoretical extensions that mere made in order to enable modeling of the closed-cell-foam material include the ability to represent internal boundaries within the RUC (to simulated internal pores) and the ability to impose an internal pressure within the simulated pores. This latter extension is crucial as two sources contribute to significant internal pressure changes within the SOFI pores. First, gas trapped in the pores during the spray process will expand or contract due to temperature changes. Second, the pore pressure will increase due to outgassing of water and other species present in the foam skeleton polymer material. With HFGMC's new pore pressure modeling capabilities, a nonlinear pressure change within the simulated pore can be imposed that accounts for both of these sources, in addition to stmdar&-thermal and mechanical loading; The triply-periodic HFGMC micromechanics model described above was implemented within NASA GRC's MAC/GMC software package, giving the model access to a range of nonlinear constitutive models for the polymeric foam skeleton material. A repeating unit cell architecture was constructed that, while relatively simple, still accounts for the geometric anisotropy of the porous foam microstructure and its thin walls and thicker edges. With the lack of reliable polymeric foam skeleton materia1 properties, many simulations were executed aimed at backing out these material properties. Then, using these properties, predictions of the thermo-mechanical behavior of the foam, including calculated internal applied pressure profiles, were performed and compared with appropriate experimental data.

Bednarcyk, Brett A.↗

DOE OSTI2020

Micro- and nano-porosity of the active Alpine Fault zone, New Zealand

Porosity reduction in rocks from a fault core can cause elevated pore fluid pressures and consequently influence the recurrence time of earthquakes. We investigated the porosity distribution in the New Zealand's Alpine Fault core in samples recovered during the first phase of the Deep Fault Drilling Project (DFDP-1B) by using two-dimensional nanoscale and three-dimensional microscale imaging. Synchrotron X-ray microtomography-derived analyses of open pore spaces show total microscale porosities in the range of 0.1 %–0.24 %. These pores have mainly non-spherical, elongated, flat shapes and show subtle bipolar orientation. Scanning and transmission electron microscopy reveal the samples' microstructural organization, where nanoscale pores ornament grain boundaries of the gouge material, especially clay minerals. Our data imply that (i) the porosity of the fault core is very small and not connected; (ii) the distribution of clay minerals controls the shape and orientation of the associated pores; (iii) porosity was reduced due to pressure solution processes; and (iv) mineral precipitation in fluid-filled pores can affect the mechanical behavior of the Alpine Fault by decreasing the already critically low total porosity of the fault core, causing elevated pore fluid pressures and/or introducing weak mineral phases, and thus lowering the overall fault frictional strength. We conclude that the current state of very low porosity in the Alpine Fault core is likely to play a key role in the initiation of the next fault rupture.

58 GEOSCIENCES↗

DOE OSTI2022

Using Direct Numerical Simulation of Pore-Level Events to Improve Pore-Network Models for Prediction of Residual Trapping of CO2

Direct numerical simulation and pore-network modeling are common approaches to study the physics of two-phase flow through natural rocks. For assessment of the long-term performance of geological sequestration of CO 2 , it is important to model the full drainage-imbibition cycle to provide an accurate estimate of the trapped CO 2 . While direct numerical simulation using pore geometry from micro-CT rock images accurately models two-phase flow physics, it is computationally prohibitive for large rock volumes. On the other hand, pore-network modeling on networks extracted from micro-CT rock images is computationally efficient but utilizes simplified physics in idealized geometric pore elements. This study uses the lattice-Boltzmann method for direct numerical simulation of CO 2 -brine flow in idealized pore elements to develop a new set of pore-level flow models for the pore-body filling and snap-off events in pore-network modeling of imbibition. Lattice-Boltzmann simulations are conducted on typical idealized pore-network configurations, and the interface evolution and local capillary pressure are evaluated to develop modified equations of local threshold capillary pressure of pore elements as a function of shape factor and other geometrical parameters. The modified equations are then incorporated into a quasi-static pore-network flow solver. The modified model is applied on extracted pore-network of sandstone samples, and saturation of residual trapped CO 2 is computed for a drainage-imbibition cycle. The modified model yields different statistics of pore-level events compared with the original model; in particular, the occurrence of snap-off in pore-throats is reduced resulting in a more frontal displacement pattern along the main injection direction. Compared to the original model, the modified model is in closer agreement with the residual trapped CO 2 obtained from core flow experiments and direct numerical simulation.

Kohanpur, Amir H.↗

DOE OSTI2026

ThermoPI—an Online Tool to Calculate Heat Transfer Through Foam Insulation

Thermal insulation materials with ultra-low effective thermal conductivity are crucial for a multitude of applications. Over the years, a significant amount of experimental work has been dedicated to creating new insulation materials with low thermal conductivity. Similarly, substantial efforts have been made to enhance the theoretical understanding of thermal transport mechanisms in thermal insulation materials and to push the boundaries of lower thermal conductivity. However, ample room remains for enhancing the thermal resistivity of closed-cell foam insulations. To aid in the development of ultra-low effective thermal conductivity foam insulations, Oak Ridge National Laboratory has introduced a unique online tool, ThermoPI. This tool calculates not only the effective thermal conductivity but also the thermal conductivity components for porous materials, including gases, solids, and radiative, based on the materials’ structural information (e.g., porosity, pore size, gas species, pressure, solid species, pore geometry, and temperature).This tool collected and improved the existing theoretical models for thermal insulation materials’ gas, solid, and radiative thermal conductivity. These improved models have been validated with experimental data. Effective medium theory and sound velocity softening effects are considered for solid thermal conductivity. The built-in solid materials include polystyrene, polyurethane, polyethylene, and silica. Users are allowed to input new solid materials that are not built inside the tool. For gas thermal conductivity, the subcontinuum Knudsen effect is considered. In addition to several built-in gases, users can input new gases. For solid and radiative thermal conductivity calculations, several models are available to select, and the tool can determine the best model to choose based on the materials that users input. Users are also allowed to change selections manually. In addition, the tool can also calculate the effect of interfacial resistance on the overall thermal conductivity of layered materials. This work will elaborate on the tool and discuss how it can guide the development of new insulation products.

Shrestha, Som [ORNL] (ORCID:0000000183993797)↗

NASA NTRS2004

Lithospheric Dynamics of Mars: Water, Flow, and Failure

Some of the largest Martian erosive features were influenced by groundwater, and include valley networks, outflow channels, and possibly landslides. We argue that hydrothermal systems attending crustal formation processes were able to drive sufficient groundwater to the surface to form the Noachian southern highlands valley networks, which show a spatial correlation to crustal magnetic anomalies, also results of crustal formation. Hydrothermal activity is quantified through numerical simulations of convection in a porous medium due to the presence of a hot intruded magma chamber. The parameter space includes magma chamber depth, volume, aspect ratio, and host rock permeability and porosity. For permeabilities as low as l0(exp -17) sq m and intrusion volumes as low as 50 km , the total discharge due to intrusions building that part of the southern highlands crust associated with magnetic anomalies spans a comparable range as the inferred discharge from the overlying valley networks. The Hesperian circum-Chryse outflow channels are further manifestations of groundwater discharge and Clifford and Parker (2001) suggest that the large volumes of water required for their formation flows beneath a confining cryosphere from the South Pole where meltwater beneath an ice cap recharges a global aquifer. We argue that recharge occurs instead over the nearby Tharsis aquifer at high obliquity, assisted by cryosphere melting due to volcanic activity. Numerical simulations quantify the strength and duration of outflow discharge given either South Polar or Tharsis recharge. The contribution of South Pole recharge given Clifford and Parker aquifer properties is negligible compared to that of the initial Tharsis inventory. Tharsis recharge, despite the restrictions of improved aquifer properties, makes a significant contribution and, unlike South Pole recharge under the same conditions, fulfills discharge requirements. Groundwater may have influenced long run-out landslide formation in the Valles Marineris. We present simulations of Martian, terrestrial, and lunar landslides that gauge the role of pore fluid pressure in reproducing accurate geometries and run-out with frictional, Bingham, and fluidization rheologies. The results indicate that pore fluid is a necessary component of Martian landslide formation and we suggest scenarios that might explain its presence.

Grimm, Robert E.↗

DOE OSTI2023

Minimizing exposure to legacy wells and avoiding conflict between storage projects: Exploring area of review as a screening tool

Elevated pressure from large-volume injection is a key driver of risk and project cost. If transmissive features (e.g., non-isolating wells or fracture systems) are present, increased injection-zone pressure can drive fluids from depth toward protected freshwater resources. In US Carbon Capture and Storage (CCS) law, the area at risk is known as the Area of Review (AoR). The size and number of potentially transmissive features to be evaluated and possibly remediated or managed is a function of the size and location of the AoR. The size of the AoR depends on several variables, including properties of the injection zone, properties of protected resources, and injection rate and duration. Evaluation of the intersection of these variables across a portfolio of sites highlights the injection zone depth and boundary conditions as top-level controls. Deep injection, use of multiple stacked injection zones, reduced injection rate and choice of injection well location can all be used to minimize AoR and the number of potentially transmissive features within it. Here, we introduce the concept of pressure space (defined as connected pore volume times pressure) as the key subsurface commodity for CO 2 storage and we suggest that it forms a more robust basis for leasing and regulation than pore space alone.

58 GEOSCIENCES↗

NASA NTRS1988

Gneiss-charnockite transformation at Kottavattam, Southern Kerala (India)

At Kottavattam, leucocratic granitic garnet-biotite gneisses (age less than 2 Ga) were partially transformed to coarse-grained charnockite along a system of conjugate fractures (N70E and N20W) and the foliation planes (N60 to 80W; dip 80 to 90 SW) about 550 m.y. ago. To examine and quantify changes in fabric, mineralogy, pore fluids and chemical composition associated with this process, large rock specimens showing gneiss-charnockite transition were studied in detail. The results of the present study corroborate the concept that charnockite formation at Kottavattam is an internally-generated phenomenon and was not triggered by the influx of carbonic fluids from a deep-seated source. It is suggested that charnockitization was caused by the following mechanism: (1) near-isothermal decompression during uplift of the gneiss complex led to an increase of the pore fluid pressure (P sub fluid greater than P sub lith) which - in a regime of anisotropic stress - triggered or at least promoted the development of conjugate fractures; (2) the simultaneous release of pore fluids from bursting fluid inclusions and their escape into the developing fracture system resulted in a drop of fluid pressure; and (3) the internal generation and buffering of the fluids and their, probably, limited migration in an entirely granitic rock system explains the absence of any significant metasomatic mass transfer.

Raith, M.↗

DOE OSTI2022

Microsized Pore Structure Determination in EPDM Rubbers Using High-Pressure 129 Xe NMR Techniques

Microsized pore parameters, such as pore size and distance between pores in a series of model EPDM rubbers, were determined in situ under the pressure of 500 psi using 129 Xe nuclear magnetic resonance (NMR) techniques: spin–lattice (T 1 ) and spin–spin (T 2 ) relaxation measurements, pulsed-field gradient (PFG) NMR, and two-dimensional exchange spectroscopy (2D EXSY). The T 1 /T 2 (>>1) ratio for the xenon confined in the pores is larger than that for nonconfined free xenon. This suggests that almost the entire pore surface interacts with xenon atoms like a closed pore. While these pores still connect each other through very narrow diffusion/exchange channels, it is possible to observe the echo decay in PFG-NMR and cross-peaks in 2D EXSY. The results show that both diffusion (D pore ≈ 2.1 × 10 –10 m 2 /s) and exchange (exchange rate, τ exch = a few tens of milliseconds) of xenon between a pore within the material and outer surface are prolonged. The exchange distances (l), which correspond to the xenon gas penetration depth, were estimated to be 70–100 μm based on the measured diffusion coefficients and exchange rate (1/τ exch ). NMR diffraction analysis reveals that pore size (a) and pore distance (b) are on the order of magnitude of micrometers and tens of micrometers, while the diffusion coefficients of xenon gas in the diffusion channels (D eff ) are about 10 –8 m 2 /s. Overall, this study suggests that the pores with a few micrometers connected through very narrow flowing channels with the length of several tens of micrometers are developed 70 to 100 μm below the rubber surface. Furthermore, the overall steady-state diffusion of xenon is slower, approximately 2 orders of magnitudes, than the diffusion in the channel between the pores. Finally, the pore and exchange distances correlated with the composition of rubbers showed that the properties of EPDM rubber as a high-pressure gas barrier could be improved by reducing the size of cracks and the depth of gas penetration by the addition of both carbon black and silica fillers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

A Dynamic Pore Network Model for Imbibition Simulation Considering Corner Film Flow

Wetting films can develop in the corners of angular pores under strong wetting conditions. Modeling the dynamics of corner film remains elusive using direct numerical simulations because of the significant scale difference between main meniscus and corner film flow. In this paper, the modified interacting capillary bundle model (ICB), developed in our previous work to describe accurately corner film dynamics in a single square tube, is incorporated into a single-pressure dynamic pore network model (DPNM) to simulate imbibition in strongly wetting porous media with corner film flow. The traditional pore network is decomposed into several layers of interacting subpore networks where the 0th layer of subpore network simulates the main meniscus flow and higher layers the corner film flow. The fluid flow between different layers is captured by interlayer throats. In addition, the snap-off mechanism caused by the thickening of wetting corner film is considered. The accuracy of the developed model is validated for four cases: spontaneous imbibition in a single square tube, wetting fluid redistribution through corner films under a capillary pressure difference, snap off in a narrow throat connecting two large pores, and imbibition dynamics in a real microfluidic porous geometry. The validated model is then used to simulate both spontaneous and controlled imbibition in a pore network with random pore size distribution. Finally, the interaction between corner film and main meniscus flow in porous media is analyzed from a pore-scale perspective.

58 GEOSCIENCES↗

DOE OSTI2022

Frictional Behavior Downdip Along the Subduction Megathrust: Insights From Laboratory Experiments on Exhumed Samples at In Situ Conditions

Subduction megathrusts exhibit a spectrum of slip modes, including catastrophic earthquakes. Although the mechanical and frictional properties of materials sampled from subduction zones have been studied extensively, few datasets have been collected for compositions and at pressure and temperature conditions representative of those in situ. The Nankai subduction zone in southwest Japan is a well-studied margin, and abundant data provide an opportunity to advance our understanding of fault and earthquake physics. Here, we use samples exhumed in the Shimanto and Sanbagawa Belts on Shikoku Island of southwest Japan that represent analogs for materials along the present-day megathrust at depths of ~5–>25 km, and we shear these at their peak in situ pressure-temperature (P-T) conditions. Effective normal stresses range from 28 to 192 MPa, and temperatures from 105°C to 470°C. We used pore fluid pressures of 45–240 MPa, corresponding to fluid overpressure ratios λ of 0.65 and 0.90. Slip velocities of 0.1–100 μm/s were used, in order to focus on the nucleation of instability and earthquakes. We found predominantly velocity-strengthening (inherently stable) behavior under all conditions for λ = 0.65. For λ = 0.90, velocity-weakening behavior was observed at 350°C, with velocity-strengthening behavior at lower and higher temperatures. The rate/state frictional stability parameter (a-b) increases with slip velocity at temperatures up to ~200°C and remains constant or decreases with slip velocity at higher temperatures. Altogether, our results demonstrate the potentially important roles of both temperature and slip velocity in controlling the distribution of stress and frictional rheology along subduction thrusts.

58 GEOSCIENCES↗

DOE OSTI2021

Sub-nanometer confinement enables facile condensation of gas electrolyte for low-temperature batteries

Abstract Confining molecules in the nanoscale environment can lead to dramatic changes of their physical and chemical properties, which opens possibilities for new applications. There is a growing interest in liquefied gas electrolytes for electrochemical devices operating at low temperatures due to their low melting point. However, their high vapor pressure still poses potential safety concerns for practical usages. Herein, we report facile capillary condensation of gas electrolyte by strong confinement in sub-nanometer pores of metal-organic framework (MOF). By designing MOF-polymer membranes (MPMs) that present dense and continuous micropore (~0.8 nm) networks, we show significant uptake of hydrofluorocarbon molecules in MOF pores at pressure lower than the bulk counterpart. This unique property enables lithium/fluorinated graphite batteries with MPM-based electrolytes to deliver a significantly higher capacity than those with commercial separator membranes (~500 mAh g −1 vs. <0.03 mAh g −1 ) at −40 °C under reduced pressure of the electrolyte.

25 ENERGY STORAGE↗

NASA NTRS1986

Knudsen diffusion - The effect of small pore size and low gas pressure on gaseous transport in soil

The analytical principles and applications of the theory of Knudsen diffusion are reviewed, with emphasis on gas transport in the soils of planetary bodies. Knudsen diffusion occurs when the mean free path of diffusing gas molecules surpasses the size of the pores through which diffusion proceeds. The process is then dominated by collisions with the pore walls. Computational techniques for deriving the Knudsen coefficient for soils with a nonreentrant cross-section shape are reviewed, along with methods of deriving a coefficient for soils which permit both Knudsen and bulk diffusion. Sample calculations for three pore-size distributions are provided to illustrate the decrease in transport efficiency with increasingly smaller soil pore sizes.

Clifford, S. M.↗

DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1985

Low-angle normal faults-low differential stress at mid crustal levels

A simple model for frictional slip on pre-existing faults that considers the local stress state near the fault and the effect of non-hydrostatic fluid pressures predicts that low-angle normal faulting is restricted to areas of the crust characterized by low differential stress and nearly lithostatic fluid pressures. The model considers frictional slip on a cohesionless low-angle normal fault governed by the failure criterion tau = mu sub f (sigma (*) sub n) =mu sub f (sigma sub n - P sub f) where tau and sigma sub n are the shear and normal stresses across the fault plane, mu sub f is the static coefficient of friction, and P sub f is the pore fluid pressure. As a first approximation, the model considers a vertical greatest principal compressive stress, sigma sub 1. It is apparent that if slip on low-angle normal faults is governed by the avove frictional failure criterion, slip on the low-angle normal fault occurs only if the least effective principal stress, sigma (*) sub 3 = sigma sub 3 - P sub f, is tensile, whenever tan superscrip -1(mu sub f d, where d is the dip of the fault. If detachment faulting occurs at any significant depth in the crust, P sub f sigma sub 3 is required. In light of this conclusion I allow P sub f to vary as necessary to allow slip on the low-angle normal fault.

Power, W. L.↗

DOE OSTI2021

The evolution of seafloor venting from hydrate-sealed gas reservoirs

Here we use 3D seismic data to show that three rows of seafloor gas mounds can be traced downward to leak points that lie at the hydrate-gas contact within three individual dipping coarse-grained sand bodies in the Terrebonne Basin, Walker Ridge block 313, northern Gulf of Mexico. We predict the overpressure within the sand bodies by assuming that the gas pressure at the vent points equal the least principal stress. We interpret that free gas accumulates at the base of the hydrate stability zone, trapped by the overlying sand which has a high methane hydrate saturation. The free gas accumulates until the gas pressure at the base of the hydrate stability zone reaches the least principal stress in the overlying mudrock, whereupon hydraulic fractures form and fluids are vented to the surface. The warm rising fluids and perhaps localized exothermic formation of hydrate raise the local salinity and temperature. This process progressively shifts the base of the hydrate stability zone to shallower depths and dissociates the hydrate seal within the sand, which creates new leak points and results in the observed migration of the seafloor vents. Within the southwest corner of the Terrebonne Basin, this process has repeated multiple times within the Blue, Orange, and Green sands. Here this study shows how the hydrate stability zone can influence the location of fluid expulsion and in turn be affected by the warm, rising, saline fluids sealed by hydrate.

58 GEOSCIENCES↗