An official website of the United States government
Temperature variations in low permeable soil (e.g. clay) induce pore pressure, which is known as thermal pressurization. Previous research showed that thermal pressurization highly depends on thermal pressurization coefficient. This coefficient depends on the soil type and changes with temperature due to temperature dependency of thermal expansion coefficient of water. Thermal pressurization is often investigated through thermo-hydro-mechanical (THM) numerical modeling. THM process, with respect to thermal loading, has been examined in the literature to justify the field observations by incorporating advanced thermo-mechanical constitutive models. However, result of numerical simulations using advanced thermo-elastoplastic models still show some discrepancies with experimental and field observations. In this study, the assessment of thermal pressurization in Boom clay is scrutinized through employing a relatively simple while practical thermo-poroelastic finite element model with careful consideration of the temperature-dependent thermal, hydraulic, and mechanical properties of the medium and saturating fluid (i.e. water). The numerical model is carried out using COMSOL Multiphysics and the results of the numerical simulations are compared and validated with the ATLAS project, a large-scale experimental facility in Belgium. The results confirm that thermal and hydraulic coupling parameters are the key factors to change thermal pressurization.
more »
« less
The Pressure‐ and Temperature‐Dependence of Thermal Pressurization in Localized Faults
Key Points Changes in hydraulic diffusivity and pressurization factor during thermal pressurization (TP) balance each other in low permeability and low porosity fault rocks Hydraulic diffusional length scales as time 0.7 when considering TP parameters that depend on temperature and pressure The constant case model should be considered with ambient initial conditions and not time‐averaged ones
more »
« less
- Award ID(s):
- 2024660
- PAR ID:
- 10451402
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 128
- Issue:
- 8
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Frictional heating during earthquake rupture raises the fault-zone fluid pressure, which affects dynamic rupture and seismic radiation. Here, we investigate two key parameters governing thermal pressurization of pore fluids – hydraulic diffusivity and shear-zone half-width – and their effects on earthquake rupture dynamics, kinematic source properties, and ground motions. We conduct 3D strike-slip dynamic rupture simulations assuming a rate-and-state dependent friction law with strong velocity weakening coupled to thermal-pressurization of pore fluids. Dynamic rupture evolution and ground shaking are densely evaluated across the fault and Earth’s surface to analyze the variations of rupture parameters (slip, peak slip rate, rupture speed, and rise time), correlations among rupture parameters, and variability of peak ground velocity. Our simulations reveal how variations in thermal-pressurization affect earthquake rupture properties. We find that the mean slip and rise time decrease with increasing hydraulic diffusivity, whereas mean rupture speed and peak slip-rate remain almost constant. Mean slip, peak slip-rate, and rupture speed decrease with increasing shear-zone half-width, whereas mean rise time increases. Shear-zone half-width distinctly affects the correlation between rupture parameters, especially for parameter pairs (slip, rupture speed), (peak slip-rate, rupture speed), and (rupture speed, rise time). Hydraulic diffusivity has negligible effects on these correlations. Variations in shear-zone half-width primarily impact rupture speed, which then may affect other rupture parameters. We find a negative correlation between slip and peak slip-rate, unlike simpler dynamic rupture models. Mean peak ground velocities decrease faster with increasing shear-zone half-width than with increasing hydraulic diffusivity, whereas ground-motion variability is similarly affected by both the parameters. Our results show that shear-zone half-width affects rupture dynamics, kinematic rupture properties, and ground shaking more strongly than hydraulic diffusivity. We interpret the importance of shear-zone half-width based on the characteristic time of diffusion. Our findings may inform pseudodynamic rupture generators and guide future studies on how to account for thermal-pressurization effects.more » « less
-
The EPA’s StormWater Management Model (SWMM) has been applied across the globe for citywide stormwater modeling due to its robustness and versatility. Recent research indicated that SWMM, with proper setup, can be applied in the description of more dynamic flow conditions, such as rapid inflow conditions. However, stormwater systems often have geometric discontinuities that can pose challenges to SWMM model accuracy, and this issue is poorly explored in the current literature. The present work evaluates the performance of SWMM 5 in the context of a real-world stormwater tunnel with a geometric discontinuity. Various combinations of spatiotemporal discretization are systematically evaluated along with four pressurization algorithms, and results are benchmarked with another hydraulic model using tunnel inflow simulations. Results indicated that the pressurization algorithm has an important effect on SWMM’s accuracy in conditions of sudden diameter changes. From the tested pressurization algorithms, the original Preissmann slot algorithm was the option that yielded more representative results for a wider range of spatiotemporal discretizations. Regarding spatiotemporal discretization options, intermediate discretization, and time steps that lead to Courant numbers equal to one performed best. Interestingly, the traditional SWMM’s link-node approach also presented numerical instabilities despite having low continuity errors. Results indicated that although SWMM can be effective in simulating rapid inflow conditions in tunnels, situations with drastic geometric changes need to be carefully evaluated so that modeling results are representative.more » « less
-
This work describes an efficient means to adjust the power level of an axial piston hydraulic pump/motor. Conventionally, the displacement of a piston pump is varied by changing the stroke length of each piston. Since the losses do not decrease proportionally to the displacement, the efficiency is low at low displacements. Here, with partial-stroke piston pressurization (PSPP), displacement is varied by changing the portion of the piston stroke over which the piston is subjected to high pressure. Since leakage and friction losses drop as the displacement is decreased, higher efficiency is achieved at low displacements with PSPP. While other systems have implemented PSPP with electric or cam-actuated valves, the pump described in this paper is unique in implementing PSPP by way of a simple, robust hydro-mechanical valve system. Experimental testing of a prototype PSPP pump/motor shows that the full load efficiency is maintained even at low displacements.more » « less
-
Injection-induced seismicity and aseismic slip often involve the reactivation of long-dormant faults, which may have extremely low permeability prior to slip. In contrast, most previous models of fluid-driven aseismic slip have assumed linear pressure diffusion in a fault zone of constant permeability and porosity. Slip occurs within a frictional shear crack whose edge can either lag or lead pressure diffusion, depending on the dimensionless stress-injection parameter that quantifies the prestress and injection conditions. Here, we extend this foundational work by accounting for permeability enhancement and dilatancy, assumed to occur instantaneously upon the onset of slip. The fault zone ahead of the crack is assumed to be impermeable, so fluid flow and pressure diffusion are confined to the interior, slipped part of the crack. The confinement of flow increases the pressurization rate and reduction of fault strength, facilitating crack growth even for severely understressed faults. Suctions from dilatancy slow crack growth, preventing propagation beyond the hydraulic diffusion length. Our new two-dimensional and three-dimensional solutions can facilitate the interpretation of induced seismicity data sets. They are especially relevant for faults in initially low permeability formations, such as shale layers serving as caprock seals for geologic carbon storage, or for hydraulic stimulation of geothermal reservoirs. This article is part of the theme issue ‘Induced seismicity in coupled subsurface systems’.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2023JB026558
