Pore Pressure Changes Research Articles

A thermo-hydromechanical damage model and its application to a deep geological radioactive repository

This study conducts a numerical analysis of long-term thermo-hydromechanical (THM) processes, with a particular focus on deep geological disposal of radioactive waste. It emphasizes the modeling of damage zones resulting from excavation activities, as well as changes in pore pressure and temperature. The numerical model presented in this paper delineates the fundamental relations of thermo-poro-elasticity. It also introduces a double phase field model specifically formulated for rock materials under compressive stress, taking into account THM coupling and time-dependent behavior. The proposed model has been implemented in a finite element code and applied to numerical analysis of short and long terms responses around a horizontal borehole in the context of the French High-Level Waste (HLW) disposal concept. In particular, displacements, temperature changes, pore fluid pressure variations as well as evolution of induced damage zones are presented and analyzed for different periods including excavation, operation and post-closure. It is found that the occurrence of damage zones is clearly related to the spacing distance between two neighboring disposal boreholes. The time-dependent deformation of host rock plays an significant role in the long-term responses.

Analysis of Fluid-Injection-Induced Seismicity Using a Dynamic Sliding Model Incorporating the Rate- and State-Dependent Friction Law

Summary Earthquakes can be triggered by fluid injection into underground formations. Fluid injection can cause large changes in the underground volume that exert stresses on nearby preexisting faults, leading to seismic activity. Assuming an increase in underground development activities in the future, our understanding of the mechanism underlying induced seismicity must be improved, and methods must be developed to properly assess the risk of seismic events. The objective of this study is to develop a seismicity prediction model that calculates the magnitude and timing of triggered earthquakes or seismic events occurring during various subsurface fluid injection activities. We developed an injection-induced seismicity analysis model that predicts the dynamic earthquake nucleation caused by changes in stress and pore pressure that occur during various subsurface activities. The governing equations consisting of the dynamic motion of the poroelastic spring-slider system, rate and state friction laws and pore pressure diffusion equation were solved using the embedded semi-implicit Runge-Kutta (SIRK) method. The dynamic sliding model was also incorporated into the finite element method (FEM) model, considering the variations in the stresses and pore pressures in the formation. A field case study was also conducted to compare the model results with typical microseismicity responses observed from hydraulic fracturing treatments in shale fields. Contrary to the popular understanding derived from Amonton’s law, the dynamic friction model revealed that a large normal stress on the fault leads to rapid sliding. A larger normal stress accumulates a large amount of elastic energy until it slips owing to fluid injection, nucleating large seismic waves. The poroelastic spring-slider model estimated reasonable microseismic magnitudes for hydraulic fracturing treatment but overestimated the time required to trigger a microseismic event under field conditions. To improve the analysis results, the poroelastic spring-slider model was coupled with a linear elastic FEM that considered the complex interplay of stress changes from hydraulic fracturing and the associated pore pressure variation in the formation. Compared with the field data, the coupled simulation model estimated a reasonable timing for the induced microseismic events when the increasing pore pressure during hydraulic fracturing penetrated deep into the formation. These findings suggest the existence of permeable natural fractures in the formation, which intensify early frictional sliding during treatment. The seismicity prediction model presented in this study simulates the magnitude and timing of seismic nucleation, helping to manage and mitigate the environmental impacts of induced seismicity during various subsurface development activities, such as oil and gas extraction, hydraulic fracturing, geothermal, and carbon dioxide sequestration. Moreover, the case study results imply that the time series of seismic events predicted by the model can be used to understand the possible fracture geometry and extent of fluid invasion for field applications. Key Terms and Phrases induced seismicity, subsurface fluid injection, rate- and state-dependent friction law, embedded semi-implicit Runge-Kutta method, finite element analysis

Feasibility of reservoir monitoring in the Groningen gas field using ghost reflections from seismic interferometry

SUMMARY Seismic interferometry (SI) retrieves new seismic responses, for example reflections, between either receivers or sources. When SI is applied to a reflection survey with active sources and receivers at the surface, non-physical (ghost) reflections are retrieved as well. Ghost reflections, retrieved from the correlation of two primary reflections or multiples from two different depth levels, are only sensitive to the properties in the layer that cause them to appear in the result of SI, such as velocity, density and thickness. We aim to use these ghost reflections for monitoring subsurface changes, to address challenges associated with detecting and isolating changes within the target layer in monitoring. We focus on the feasibility of monitoring pore-pressure changes in the Groningen gas field in the Netherlands using ghost reflections. To achieve this, we utilize numerical modelling to simulate scalar reflection data, deploying sources and receivers at the surface. To build up subsurface models for monitoring purposes, we perform an ultrasonic transmission laboratory experiment to measure S-wave velocities at different pore pressures. Applying SI by autocorrelation to the modelled data sets, we retrieve zero-offset ghost reflections. Using a correlation operator, we determine time differences between a baseline survey and monitoring surveys. To enhance the ability to detect small changes, we propose subsampling the ghost reflections before the correlation operator and using only virtual sources with a complete illumination of receivers. We demonstrate that the retrieved time differences between the ghost reflections exhibit variations corresponding to velocity changes inside the reservoir. This highlights the potential of ghost reflections as valuable indicators for monitoring even small changes. We also investigate the effect of the sources and receivers’ geometry and spacing and the number of virtual sources and receivers in retrieving ghost reflections with high interpretability resolution.

Azimuthal pore pressure response to teleseismic waves: effects of damage and stress anisotropy

SUMMARY Pore pressure oscillations induced by stress variations, including propagating seismic waves from remote earthquakes, have been widely observed in various groundwater systems. The monitored pressure change in wells shows significant water-level oscillations to volumetric strain as well as to S and Love waves. Recent observations demonstrated azimuthal dependence of the pore pressure oscillations with respect to stress indicators and fault zone orientation. Within the fault zone, damage-induced anisotropy is the result of the alignment and orientation of cracks and other internal flaws within the rock. In this work, we provide a complete quantitative description of the pore pressure changes induced by passing seismic waves associated with different orientations and values of principal stress and damage tensor components. The model quantifies the azimuthal dependence of the pore pressure response by a non-dimensional ratio defined as the amplitude of the pressure oscillations induced by a shear strain normalized to the volumetric strain. Three angles and two values are needed to calculate the azimuthal dependence of the pore pressure response: the angle between the directions of the maximum horizontal stress and the seismic event; fault zone orientation; microcrack orientation within the fault zone; and damage and stress values. The model predicts that maximum pore pressure response occurs when microcracks and maximum horizontal stress are in the same orientation, high damage and high stress anisotropy. By adjusting these quantities, we recalculate results of recent seismological studies in the Arbuckle disposal well, Osage County, Oklahoma. The presented model successfully predicts the observed azimuthal dependence in wave-induced fluid pressure response and relates the anisotropic response to tectonic indicators such as the orientations of the maximum horizontal stress, fault zone, and microfractures.

Decomposing significant factors of Coulomb stress and its components in injection-induced seismicity

Injection-induced seismicity has been a focus of industry for decades as it poses great challenges to the associated risk mitigation and hazard assessment. The response surface methodology is integrated into the geo-mechanical model to analyze the effects of multiple factors on induced seismicity during the post shut-in period. We investigate the roles of poroelastic stress and pore pressure diffusion and examine the differences in the controlling mechanism between fault damage zones and the fault core. A sensitivity analysis is conducted to rank the selected factors, followed by a Box‒Behnken design to form response surfaces and formulate prediction models for the Coulomb stress and its components. Reservoir properties significantly affect the potentials of induced seismicity in the fault by changing pore pressure diffusion, which can be influenced by other factors to varying degrees. Coulomb stress is greater in pressurized damage zones than in fault cores, and the seismicity rate exhibits a consistent variation. Poroelastic stress plays a similar role to pore pressure diffusion in the stability of the fault within the pressurized damage zones. However, pore pressure diffusion dominates in the fault core due to the low rigidity, which limits the accumulation of elastic energy caused by poroelastic coupling. The slip along the fault core is a critical issue to consider. The potential for induced seismicity is reduced in the right damage zones as the pore pressure diffusion is blocked by the low-permeability fault core. However, poroelastic stressing still occurs, and in deep basements, the poroelastic effect is dominant even without a direct increase in pore pressure. The findings in this study reveal the fundamental mechanisms behind injection-induced seismicity and provide guidance for optimizing injection schemes in specific situations.

Effects of local thermal non-equilibrium on hydraulic stimulation efficiency of enhanced geothermal systems

Motivated by hydraulic stimulation of enhanced geothermal systems, the present paper investigates the coupled thermo-hydro-mechanical response of a geothermal well imbedded in a thermoporoelastic medium, subjected to a non-isothermal fluid flux and convective cooling on the borehole surface. Our focus centers on the effect of local thermal non-equilibrium (LTNE) on the temporal-spatial evolution of temperatures, pore pressure, and stresses, where the solid and fluid phases have two distinct temperatures and local heat transfer between the two phases is addressed. We employ integral transform and load decomposition techniques to derive analytical solutions in the Laplace domain. This methodology allows us to disentangle and separate the individual contributions to changes in pore pressure and stresses from fluid injection and convective heat transfer. The results reveal that compared to the classical local thermal equilibrium model, the thermally induced pore pressure is slightly lower under LTNE conditions. The LTNE has a significant influence on the temporal evolution of thermally induced stresses, especially in the vicinity of the wellbore.

Sealing capacity evaluation of underground gas storage under intricate geological conditions

Evaluating underground gas storage (UGS) sealing capacity is essential for its safe construction and operational efficiency. This involves evaluating both the static sealing capacity of traps during hydrocarbon accumulation and the dynamic sealing capacity of UGS under intensive gas injection and withdrawal, and alternating loads. This study detailed the methodology developed by Sinopec. The approach merges disciplines like geology, geomechanics, and hydrodynamics, employing both dynamic-static and qualitative-quantitative analyses. Sinopec's evaluation methods, grounded in the in situ stress analysis, include mechanistic studies, laboratory tests, geological surveys, stress analysis, and fluid-solid interactions. Through tests on the static and dynamic sealing capacity of UGS, alongside investigations into sealing mechanisms and the geological and geomechanical properties of cap rocks and faults, A geomechanics - rock damage - seepage mechanics dynamic coupling analysis method has been developed to predict in situ stress variations relative to pore pressure changes during UGS operations and evaluate fault sealing capacity and cap rock integrity, thereby setting the maximum operational pressures. Utilizing this evaluation technique, Sinopec has defined performance metrics and criteria for evaluating the sealing capacity of depleted gas reservoirs, enabling preliminary sealing capacity evaluations at UGS sites. These evaluations have significantly informed the design of UGS construction schemes and the evaluation of fault sealing capacity and cap rock integrity during UGS operations.

On the Derivation of Closed‐Form Expressions for Displacements, Strains, and Stresses Inside Poroelastic Reservoirs

AbstractWe critically review the derivation of closed‐form analytical expressions for elastic displacements, strains, and stresses inside a subsurface reservoir undergoing pore pressure changes using inclusion theory. Although developed decades ago, inclusion theory has been used recently by various authors to obtain fast estimates of depletion‐induced and injection‐induced fault stresses in relation to induced seismicity. We therefore briefly address the current geomechanical relevance of this method, and provide a numerical example to demonstrate its use to compute induced fault stresses. However, the main goal of our paper is to correct some erroneous assumptions that were made in earlier publications. While the final expressions for the poroelastic stresses in these publications were correct, their derivation contained conceptual mistakes due to the mathematical subtleties that arise because of singularities in the Green's functions. The aim of our paper is therefore to present the correct derivation of expressions for the strains and stresses inside an inclusion and to clarify some of the results of the aforementioned studies. Furthermore, we present two conditions that the strain field must satisfy, which can be used to verify the analytical expressions.

Technology Focus: Seismic (February 2024)

2023 wrapped up nicely with COP28 in Dubai. The key highlights were the commitments by both governments and the private sector to zero-methane production by 2030, net-zero greenhouse gas emissions from operations by 2050, and triple renewable power-generation capacity and double energy efficiency in this decade. There has been tremendous progress on our collective journey toward cleaner energy production, and now we must take it a notch higher. In the past couple of years, JPT’s Seismic feature has been marked by highly innovative seismic techniques and machine-learning work flows focused on reduction of carbon footprint; frontier exploration, expanding heartlands and brownfield development; and improved seismic survey, imaging, and attribute analysis for reservoir characterization and fracture diagnostics, which are in line with the COP28 pledges. This year’s feature is no different. Recent technical-paper titles have further shown the steady increase in the application of advanced seismic techniques and machine learning to mature “stranded” and “advantaged” hydrocarbon-bearing accumulations to production; improve carbon capture, storage, and leak detection; and analyze naturally and artificially induced fractured reservoirs and seismicity. The selected papers are reflective of these themes. In addition, some remarkable applications of seismic techniques are presented in the recommended papers. They include the following: - 2D nonlinear seismic site response analyses used to determine if the capacity of a silty sedimentation could support the load of a pipeline end termination in the event of an earthquake - The results of a pilot test of the Marine Vibrator, an industry-leading seismic-survey technique - High-resolution modeling to infer if an increase in seismicity is caused by an increase in pore pressure from wastewater injection in the Fort Worth Basin of north central Texas - 4D time/depth shifts and strain signals successfully used to estimate pressure depletion in the Gulf of Mexico Others are applying advanced measurement and modeling techniques to optimally place a long horizontal injector well in multilayered heterogeneous sandstone reservoirs and predicting reservoir presence from seismic velocity and pore-pressure prediction. These topics will supply the insights that will help us achieve the COP28 targets while we look forward to more exciting innovations in 2024. Recommended additional reading at OnePetro: www.onepetro.org. OTC 32267 Seismic Stability Assessment of a Mudmat on Liquefiable Seabed by Brian Carlton, Norwegian Geotechnical Institute, et al. OTC 32393 Predicting Reservoir Presence From Seismic Velocity Mapping and Pore-Pressure Prediction by Matthew James Legg, Shell OTC 32741 Marine Vibrator Milestone: A Pilot Seismic Survey by R. Alfaro, TotalEnergies, et al. SPE 212951 High-Resolution Modeling of Pore-Pressure Change, Fault Slip Potential, and Induced Seismicity in the Fort Worth Basin by Changqing Yao, Texas A&M University, et al. URTeC 3870451 Impact of Fracture Roughness on Fines Migration and Fracture-Aperture Growth in Calcareous Shale Rocks During Acidized Corefloods by Hasan J. Khan, King Fahd University of Petroleum and Minerals, et al.

An interpretation of data after performing Step-rate-test in multi-layered, unconsolidated and weakly consolidated reservoirs

The prerequisite for this study is the difficulty in interpreting data after stimulation of the bottomhole formation zone has been carried out. The objectives of this work are the applicability and interpretation of Step-rate-test (SRT) data to evaluate suitable injection regimes in multilayered, unconsolidated and weakly consolidated reservoirs and the demonstration of the correlation between injection rate and pressure in a graphical format. We investigated the stability of the reservoir structure to pressure gradients by increasing pressure and injecting fluid through the injection well stock. We concluded that this stability does not improve over time. This means that the practical significance of this method lies in its ability to alter the pressure gradient when injecting fluid to maintain a constant flow rate. The constancy of the pressure gradient has a major impact on the stimulation of the bottomhole formation zone and therefore it can be concluded that the reservoir properties in the bottomhole formation zone are dynamic in nature and can change over time. It should be noted that the study found that a change in pore pressure within the loose rock can also affect the phase (effective) permeability, which in turn can affect the required pressure gradients during injection.

A multi-objective optimization and multi-attribute decision-making analysis for technical-thermodynamic-economic evaluation considering the rock damage on production performance of hot dry rock geothermal resources

In the long-term mining of geothermal resources in hot dry rock (HDR), the change of thermal stress and pore pressure will increase fracture conductivity evolution, further improving production performance. The optimization and decision-making of the development scheme based on the impact of damage from fractures have yet to be reported. The damage to fractures is essential in designing and adjusting geothermal resource development schemes, particularly in selecting optimal schemes. Therefore, the production performances of HDR resources under different parameters are analyzed to establish a database. Then, minimizing flow resistance, maximizing net power, and maximizing economic benefits are set as optimization goals. Various injection-mining parameters and fracture characteristics are treated as decision variables. Multi-objective optimization and multi-attribute decision analysis is conducted to obtain optimal schemes. Finally, optimal schemes are evaluated and compared, considering damage and non-damage scenarios. Results show that the NSGA-II algorithm is more suitable for optimizing geothermal development questions. Net power and economic benefits of the optimal scheme considering damage increase by 45.84 % and 21.35 % compared to the control scheme with damage. For the non-damage scenario, the above values increased by 31.55 % and 5.15 %, respectively. Compared to not considering the damage, higher mass flow and well spacing of optimal scheme can be selected for the case when damaged. Moreover, the parametric design of the optimal scheme becomes more conservative as the production cycle increases.

A novel apparent permeability model for shale considering the influence of multiple transport mechanisms

Changes in pore pressure during the extraction of shale gas lead to dynamic alterations in the pore structure and permeability, making it challenging to gain a comprehensive understanding of the flow behaviors of shale gas. The pore structure of shale is complex, with a variety of storage modes and gas transport processes constrained by a number of factors. For instance, when gas flows through a transport channel with a finite length, it is imperative to take into account the flow loss caused by the bending of inlet and outlet streamlines, prior models typically neglect the impact of end effects, resulting in an exaggerated estimation of the shale permeability. Furthermore, a decrease in pore pressure corresponds to an increase in the Knudsen number, resulting in the breakdown of the continuity assumption of the Navier–Stokes equation, this signifies the gradual shift of the transport regimes from continuum flow to other transport regimes. The gas flow process is nonlinear due to the alternating impact of multicomponent transport mechanisms and various microscale effects. In this paper, we presented a novel apparent permeability model for shale that incorporates the impact of real gas effect, end effects, transport regimes, adsorption, and effective stress. First, we assumed the channel for shale gas transport to be circular pore and calculated the viscosity under the influence of a real gas effect as well as the corresponding Knudsen number. Subsequently, building upon the foundation of the slip model, we introduce the influence of the end effects to establish a bulk phase permeability for shale, further considering the impact of surface diffusion. Then, the pore radius was quantified under the influences of adsorption and effective stress. Using the intrinsic correlation between permeability and pore radius as a bridge, a shale apparent permeability model was further derived. The model encompasses various transport regimes and microscale effects, replicating the gas flow behaviors in shale. The new model was verified through comparison with published experimental data and other theoretical models, while analyzing the evolution of apparent permeability. Additionally, this paper discusses the influence of various factors, including end effects, pore radius, internal swelling coefficient, sorption-induced strain, and model-related parameters on the shale apparent permeability.

A fully coupled model for predicting geomechanical and multiphase flow behaviour in fractured rocks

Geomechanical and multiphase flow characteristics are essential in recovering oil from naturally fractured rocks during hydrocarbon production because of changes in pore pressure and tension within the rock. It is a well-established fact that the geomechanical and multiphase flow characteristics of fractured rocks are interdependent on each other. Evaluation of these characteristics, for hydrocarbons displaced by water in fractured rocks under external stress loading, is severely lacking in published literature. This study aims to develop a novel numerical framework for a fully coupled model of fractured rocks, taking into consideration the pore pressure and porous media discontinuity at the fracture-matrix interface, along with an expanded Darcy's equation. The fully coupled Finite Element Method (FEM) and Computational Fluid Dynamics (CFD) model developed in this study is shown to accurately predict geomechanical and multiphase flow behaviour at the fracture-matrix interface. The results show that as external stress loading on the fractured rock increases, the porosity and permeability of the rock matrix decrease, capillary pressure at the fracture-matrix interface decreases, and the relative permeability curves shift to the right, indicating a water-soaked fracture-matrix interface. The findings of this study can be used to develop innovative strategies for enhanced oil recovery from fractured rocks.

AN ALTERNATIVE PROCESS-BASED APPROACH TO PREDICTING THE RESPONSE OF WATER-SATURATED POROUS MEDIA TO HARMONIC HYDRODYNAMIC LOADS

Methods have been developed to predict how hydrodynamic loads acting on nearly saturated porous media are transmitted to the subsoil. In line with the effective stress principle of Terzaghi, these methods apply the boundary conditions that the effective stresses at the surface of a porous medium are zero, and that the pore water pressures carry the full load. Here, a new approach is presented which is based on defining a stress and a stress gradient as boundary conditions. The stress gradient follows from the momentum balance equation, thereby assuring that the solution abides by d'Alembert's principle of minimization of virtual work. The corresponding solution is in full accordance with the volume and momentum balance equations of the linear elastic soil matrix and the volume and momentum balance equations of the pore water across the computational domain. The new method is thereby able to correctly reproduce measurements of pore pressure changes due to hydrodynamic loads under the assumption of a porous medium consisting of incompressible particles and pore water which could either be compressible or incompressible. The advantage of the proposed method is that it requires one less boundary condition at the surface of the porous medium. The method is therefore able to predict the magnitude of the effective stresses on a soil surface. Due to the ability to retain the assumption of incompressible water, the method has also become independent on a calibration parameter. The results of the method induce questions with respect to the validity of Terzaghi's principle of effective stress at the boundary when porous media are subjected to hydrodynamic loads.

Seismic Velocity Variations Observed Prior to the La Palma Volcano Eruption on 19 September 2021, in Cumbre Vieja, Canary Islands (Spain)

Abstract In the recent years, coda-wave interferometry from seismic noise correlation functions has been increasingly used for volcanic eruption forecasting through velocity changes observed in the crust. Because La Palma Island in the Canary Archipelago is very well instrumented, we studied the possible velocity variations related to the last Cumbre Vieja eruption on 19 September 2021, aiming to obtain clear variations in the seismic velocity. For this purpose, we used the moving-window cross-spectral analysis technique for seismic noise within the 0.1–1.0 Hz frequency interval for determining two- and single-station cross-component correlations. During the 2018–2022 observation period, we first detected a seasonal seismic velocity variation possibly caused by annual rainfall and the induced pore pressure change. On 12 September 2021, a dramatic decrease in the velocity of −0.15% was detected, leading to the volcanic eruption at Cumbre Vieja seven days later. The results are compatible with those of models proposed for rapid magma migration from a shallow reservoir at 11 km to the surface.

Three-dimensional hydro-mechanical coupling numerical simulation of shield-driven cross-river twin tunnels: A case study

With the rapid development of urban underground space, the construction of shield-driven cross-river twin tunnels is increasing, and the complex hydro-mechanical coupling effects and twin-tunnel interactions bring huge construction risks to such projects, which have attracted more and more attention. This study aims to understand the excavation effects induced by shield driving of cross-river twin tunnels through numerical simulation. A refined three-dimensional numerical model based on the fully coupled hydro-mechanical theory is established. The model considers the main components of the slurry pressure balance shield (SPBS) machine, including support force, jacking thrust, grouting pressure, shield-rock interaction and lining-grouting interaction, as well as the detailed construction process. The purpose is to examine the excavation effects during construction, including rock deformation around tunnels, the change in pore pressure, and the response of the lining. The results show the influence range of twin-tunnel excavation on rock deformation and pore pressure, as well as the modes of lining response. In addition, this study also systematically investigates the effects of water level fluctuation and burial depth on twin-tunnel excavation. The results indicate that the increase of water level or burial depth will enhance the excavation effects and strengthen the twin-tunnel interactions. These results provide useful insights for estimating the construction impact range and degree of twin tunnels, and serve as basic references for the design of cross-river twin tunnels.

The Porothermoelastic Stress Evolution of Inhomogeneous Solid due to Near-Surface Reservoir

Anisotropic and heterogeneous solids, comprising a minimum of two or more elements with different properties, appear pervasively in rock materials including the pore structure and mineral composition of quartz and granite and usually have extensive applications in construction aggregates, dimension stone, geotechnical engineering, and petrology/geology. The elastic stress fields of inhomogeneous materials or composites inevitably change due to the presence of heterogeneous microstructure under applied external conditions, and hence the total mechanical and physical properties are affected by the temperature and pore pressure fluctuations beneath the surface. In this paper, the thermoporoelastic (TPE) approach on the basis of eigenstrain concept is introduced to predict the stress fields in fluid-saturated porous geologic materials like hydrocarbon reservoirs or aquifers by accounting for the coupling between thermal, poroelastic, and mechanical effects. It is an extension of the micromechanical theory that also incorporates thermal effects like pore pressure changes and temperature alternations. In addition, the TPE approach provides an important multiphysical modeling tool for understanding subsurface fluid-rock interactions and stress states in applications like unconventional oil/gas and geothermal energy.

Numerical analysis of the post-installation consolidated response of skirted foundations

The effect of the installation process on the subsequent foundation consolidation response requires assessment when the preload is initially applied. In the framework of effective stress, the influence of large soil deformations during penetration is analysed in this study using a comparison of results from large deformation finite element (LDFE) analysis based on periodic mesh regeneration (Ritss technique) and a 'wished-in-place' small strain finite element (SSFE) model. The LDFE results demonstrate that considerations of the skirt penetration and soil plug bulge lead to different soil excess pore pressure and void ratio change distributions from that of the SSFE analysis, thereby affecting the foundation settlement and consolidated bearing capacity. The normalised fully consolidated undrained bearing capacity grows linearly with the preload ratio for preloads higher than the penetration resistance, and the capacity obtained from the LDFE is higher. The settlement ratio and the proportion of partially consolidated to fully consolidated capacity gain can be represented as a hyperbolic function of a revised dimensionless time factor. The methods integrated in this study can be used to predict the increase in the vertical undrained capacity of skirted foundations after consolidation, accounting for installation effects in engineering practice.

Stress‐Induced Anisotropic Poroelasticity in Westerly Granite

AbstractWe measured poroelastic properties of cracked granite under triaxial conditions, at elevated confining pressure and a range of differential stresses. Skempton's coefficients and undrained Young's modulus and Poisson's ratio were determined directly by recording in situ fluid pressure during rapid cycles of axial and radial stress. Drained properties were measured both statically and dynamically at ultrasonic frequencies. At a given confining pressure, increasing differential stress leads to the development of elastically transverse isotropy, with symmetry axis aligned with the compression axis. Skempton's coefficients are also anisotropic, with larger changes in pore pressure in response to radial stress (coefficient Bx) than to axial stress steps (Bz). The anisotropy in the Skempton coefficients increases with increasing differential stress, with Bz decreasing and Bx slightly increasing. The evolution of static moduli and the Skempton coefficients is well approximated by Gassmann's equation using dry moduli obtained from ultrasonic measurements. Simplified predictions of the Skempton coefficients based on crack density tensors inverted from dynamic data also shows acceptable agreement with direct observations. Perfect quantitative agreement is not reached, due to the imprecision of our dynamic measurements, model simplifications, and inherent differences between static moduli obtained using stress steps of several MPa stress and dynamic ultrasonic stress oscillations.

FEM-BEM modeling dynamic responses of saturated soil around a fluid-filled lined tunnel subjected to water hammer

Dynamic stress responses of saturated soil around a fluid-filled lined tunnel caused by a water hammer are investigated by a frequency-domain FEM-BEM coupled model. The fluid is modeled as an inviscid and compressible fluid, the lining is modeled by elastic medium and conceptualized as a hollow cylinder of finite length, and the saturated poroelastic medium is adopted to model the soil. Initially, governing equations of fluid and those of lining are solved by FEM in the frequency domain, while those of soil are solved by BEM in the same domain. In the following, fluid, lining, and soil are coupled based on the conditions of deformation compatibility and force balance on their interfaces. Water pressure (inside the tunnel), the distribution of lining displacement and dynamic stress responses of saturated soil generated by the water hammer are presented. It is concluded that the dynamic stresses and the pore pressure change periodically in saturated soil under a water hammer. Modeling soil as an elastic medium inaccurately evaluates the distribution of lining displacement. The soil permeability has a significant influence on the normal stresses of soil and pore pressure but has a slight effect on the shear stresses of soil.