决定饱和岩土材料变形的有效应力及孔压系数

决定饱和岩土材料变形的有效应力 以及孔压系数B的表达式主要有2种不同的形式,为辨析其合理性,分析了不排水条件下,孔压增量和净围压增量作用下固相与骨架的体积压缩系数之间的关系,以及分别由孔压增量和净围压增量作用引起的液相、固相和骨架的体积增量,指出了 和B的两种形式的表达式在推求过程中各自存在的问题,并通过不排水和排水条件下的变形等效关系,建立了新的有效应力 和孔压系数B的表达式。     

一种适用于孔隙体积应变的有效应力方程

土力学奠基石Terzaghi有效应力原理被广泛应用于油藏孔隙和渗透率应力敏感研究中,然而其对于岩石孔隙体积应变的适用性存在争议。对颗粒不可压缩和颗粒可压缩的多孔介质分别进行了受力分析,推导了总体积、颗粒骨架、孔隙体积的有效应力表达式,与Biot、Skepmton有效应力方程对比,建立了适用于孔隙体积应变的新型有效应力方程,并进行了试验论证和应用举例。研究表明:在颗粒不可压缩多孔介质中,有效应力为超出平衡孔隙流压之外的颗粒间宏观等效应力;在颗粒可压缩变形多孔介质中,有效应力为其相同应变下的等效应力,有3种有效应力分别适用于总体积应变、颗粒体积应变、孔隙体积应变;新提出的孔隙体积有效应力方程与孔隙度、岩石总体积压缩系数、颗粒压缩系数、总应力和流压相关,4个理论计算式计算结果在3种多孔介质试验测试结果中的偏差均在5%以内;孔隙体积有效应力系数解决了如何定量增总应力来等效模拟储层降流压生产过程这一关键问题,3个压缩系数关系式理论计算准确方便。     

Experimental study of poromechanical behavior of saturated claystone under triaxial compression

This paper presents experimental results of drained and undrained triaxial compression tests of saturated Meuse–Haute/Marne claystone. The emphasis is to study the evolution of pore pressure with growth of microcracks and the effect of pore pressure on mechanical behavior. Basic mechanical responses are first investigated through drained triaxial compression tests, showing nonlinear stress strain relations, volumetric dilatancy and pressure sensitivity. In undrained triaxial compression tests, the pore pressure exhibits a transition from increase to decrease due to the transition from volumetric compressibility to dilatancy caused by the growth of microcracks. The failure surfaces, determined by total stress and Terzaghi’s effective stress under undrained condition, are compared with the one under drained condition.     

具有球形空腔饱和土分数导数 黏弹性土的动力特性

将土骨架视为具有分数阶导数本构关系的黏弹性体,基于Biot两相饱和介质模型,建立具有球形空腔饱和分数导数黏弹性土体稳态振动的控制方程。通过引入势函数,得到球对称情形下具有球形空腔饱和分数导数黏弹性土体的位移、应力和孔隙流体压力等解析表达式。考察分数导数模型参数和饱和土参数等对土体振动特性的影响,结果表明,流体压缩性对饱和土体的动力特性有显著影响,而土骨架压缩性和流-固耦合系数的影响相对较小;分数导数阶数对土体动力特性的影响与材料参数比的取值有关。同时,边界不排水条件下饱和土体的动力响应大于排水条件下饱和土体的动力响应。     

Stress dependent thermal pressurization of a fluid-saturated rock

Summary  Temperature increase in saturated porous materials under undrained conditions leads to thermal pressurization of the pore fluid due to the discrepancy between the thermal expansion coefficients of the pore fluid and of the solid matrix. This increase in the pore fluid pressure induces a reduction of the effective mean stress and can lead to shear failure or hydraulic fracturing. The equations governing the phenomenon of thermal pressurization are presented and this phenomenon is studied experimentally for a saturated granular rock in an undrained heating test under constant isotropic stress. Careful analysis of the effect of mechanical and thermal deformations of the drainage and pressure measurement system is performed and a correction of the measured pore pressure is introduced. The test results are modelled using a non-linear thermo-poro-elastic constitutive model of the granular rock with emphasis on the stress-dependent character of the rock compressibility. The effects of stress and temperature on thermal pressurization observed in the test are correctly reproduced by the model. Correspondence: Dr. S. Ghabezloo, CERMES, Ecole Nationale des Ponts et Chaussées, 6-8 avenue Blaise Pascal, Cité Descartes, 77455 Champs-sur-Marne, Marne la Vallée cedex 2, France     

A dynamic two‐phase flow model for air sparging

Air sparging (AS) is an in situ soil/groundwater remediation technology, which involves the injection of pressurized air/oxygen through an air sparging well below the zone of contamination. Characterizing the mechanisms governing movement of air through saturated porous media is critical for the design of an effective cleanup treatment system. In this research, micromechanical investigation was performed to understand the physics of air migration and subsequent spatial distribution of air at pore scale during air sparging. The void space in the porous medium was first characterized by pore network consisting of connected pore bodies and bonds. The biconical abscissa asymmetric concentric bond was used to describe the connection between two adjacent pore bodies. Then a rule‐based dynamic two‐phase flow model was developed and applied to the pore network model. A forward integration of time was performed using the Euler scheme. For each time step, the effective viscosity of the fluid was calculated based on fractions of two phases in each bond, and capillary pressures across the menisci was considered to compute the pressure field. The developed dynamic model was used to study the rate‐dependent drainage during air sparging. The effect of the capillary number and geometrical properties of the network on the dynamic flow properties of two‐phase flow including residual saturation, spatial distribution of air and water, dynamic phase transitions, and relative permeability‐capillary pressure curves were systematically investigated. Results showed that all the above information for describing the air water two‐phase flow are not intrinsic properties of the porous medium but are affected by the two‐phase flow dynamics and spatial distribution of each phase, providing new insight to air sparging. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of fully saturated porous materials

Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd.     

特察模型在孔隙性岩石中应用的局限性浅析

地层压实过程中的静力平衡涉及到异常地层压力的形成和地层压力的计算方法。本文认为经典的特察模型不利于表达孔隙性岩石中微观的静力平衡;有效应力并非与孔隙度有直接对应关系;颗粒应力是岩石骨架中传递的真正应力,而且在纵向上颗粒应力变化比较复杂。这些结论影响到地层压力计算模型中的一些基本假设条件。     

Instability criteria for quasi‐saturated viscous soils

This paper presents a theoretical framework to interpret the inception of unstable undrained creep in quasi‐saturated soils. For this purpose, the effect of gas bubbles occluded in the fluid phase is embedded into an augmented compressibility of the fluid mixture, while the mechanical characteristics of the solid skeleton have been simulated through a viscoplastic strain‐hardening model. This constitutive framework has been been used to formulate a theoretical platform able to detect runaway failures resulting from extended stages of undrained creep. It is shown that the conditions identifying the onset of spontaneous accelerations are governed by the same stability index associated with the initiation of static liquefaction. At variance with soils saturated by incompressible fluids, the conditions for undrained instability are altered by the appearance of the Skempton coefficient B, thus reflecting the beneficial effect of the fluid compressibility and its ability to decrease the liquefaction potential. The capabilities of the theory are verified through a sequence of undrained creep simulations showing the transition from stable to unstable behavior resulting from an increase of the degree of saturation. The proposed findings provide a conceptual framework to interpret the effects of gas bubbles in loose soils, as well as to assess effectiveness and longevity of liquefaction mitigation strategies based on desaturation technologies.     

Coupling of solid deformation and pore pressure for undrained deformation—a discrete element method approach

This paper presents a numerical scheme for fluid‐particle coupling that uses the discrete element method by taking into consideration solid deformation and pore pressure generation. A new water particle element is introduced to calculate pore water pressure due to porosity changes. The water particle element has the same size and shape as the solid element and experiences the same amount of deformation. On the basis of the effective stress principle at the element contact, the total force is equal to the sum of the force transmitted through the solid element contact and the water particle force due to pore water pressure. Analytical solutions of traditional soil mechanics problems, such as isotropic compression and consolidated triaxial undrained test, are used to quantitatively validate the proposed model. The numerical results show good agreement between the model and the analytical solutions. The model therefore provides an effective method to calculate pore pressure in a porous medium in discrete modeling.     

Numerical study of the deformation of saturated soil in the vicinity of a vibrating pile

The paper presents a dynamic finite-element analysis of the deformation of saturated cohesionless soil in the vicinity of the toe of a vibrating cylindrical pile. The soil behaviour is described by a hypoplastic constitutive model with intergranular strain assuming locally undrained conditions for the pore fluid. Emphasis is placed on the detailed analysis of stresses and displacements in the soil. It is shown that, starting from a homogeneous stress state, the first several cycles of vibration lead to the formation of a permanent liquefaction zone with vanishing effective stresses at a certain distance from the pile. The displacement field reveals an accumulation of residual displacements in the soil in the form of rotation. The influence of the initial stress state, the pore fluid compressibility, the pile displacement amplitude and the soil density on the formation of a liquefaction zone, on the stress amplitude in the soil and on the intensity of the permanent rotation is investigated.     

非饱和土一维固结混合流体方法的解析分析

对采用混合可压缩流体方法分析非饱和土一维固结问题的固结方程进行了求解,在得到的解析解的基础上,对影响非饱和土一维固结的因素进行了分析。分析结果表明,在采用混合流体方法计算非饱和土一维固结的孔隙水压力时,所用公式与计算饱和土一维固结的太沙基理论公式基本相同,不同之处在于引入Bishop有效应力系数来体现孔隙气对孔隙水的影响。而在非饱和土孔隙气压的计算公式中除了体现孔隙水对孔隙气的影响参数以外,还有体现孔隙气体的可压缩性对固结影响的参数。在所有影响因素中,影响非饱和土一维固结最重要的因素是孔隙流体的渗流路径。     

Numerical stability for modelling of dynamic two‐phase interaction

Dynamic two‐phase interaction of soil can be modelled by a displacement‐based, two‐phase formulation. The finite element method together with a semi‐implicit Euler–Cromer time‐stepping scheme renders a discrete equation that can be solved by recursion. By experience, it is found that the CFL stability condition for undrained wave propagation is not sufficient for the considered two‐phase formulation to be numerically stable at low values of permeability. Because the stability analysis of the two‐phase formulation is onerous, an analysis is performed on a simplified two‐phase formulation that is derived by assuming an incompressible pore fluid. The deformation of saturated porous media is now captured in a single, second‐order partial differential equation, where the energy dissipation associated with the flow of the fluid relative to the soil skeleton is represented by a damping term. The paper focuses on the different options to discretize the damping term and its effect on the stability criterion. Based on the eigenvalue analyses of a single element, it is observed that in addition to the CFL stability condition, the influence of the permeability must be included. This paper introduces a permeability‐dependent stability criterion. The findings are illustrated and validated with an example for the dynamic response of a sand deposit. Copyright © 2015 John Wiley & Sons, Ltd.     

Explicit solutions for the instantaneous undrained contraction of hollow cylinders and spheres in porous elastoplastic medium

In this article we present closed‐form solutions for the undrained variations in stress, pore pressure, deformation and displacement inside hollow cylinders and hollow spheres subjected to uniform mechanical pressure instantaneously applied to their external and internal boundary surfaces. The material is assumed to be a saturated porous medium obeying a Mohr–Coulomb model failure criterion, exhibiting dilatant plastic deformation according to a non‐associated flow rule which accounts for isotropically strain hardening or softening. The instantaneous response of a porous medium submitted to an instantaneous loading is undrained, i.e. without any fluid mass exchange. The short‐term equilibrium problem to be solved is now formally identical to a problem of elastoplasticity where the constitutive equations involve the undrained elastic moduli and particular equivalent plastic parameters. The response of the model is presented (i) for extension and compression undrained triaxial tests, and (ii) for unloading problems of hollow cylinders and spheres through the use of appropriately developed closed‐form solutions. Numerical results are presented for a plastic clay stone with strain hardening and an argilite with strain softening. The effects of plastic dilation, of the strain softening law and also of geometry of the cavity on the behaviour of the porous medium have been underlined. Analytical solutions provide valuable benchmarks enabling various numerical methods in undrained conditions with a finite boundary to be verified. Copyright © 2002 John Wiley & Sons, Ltd.     

An elasto‐plastic three phase model for partially saturated soil for the finite element simulation of compressed air support in tunnelling

This paper presents a fully coupled finite element formulation for partially saturated soil as a triphasic porous material, which has been developed for the simulation of shield tunnelling with heading face support using compressed air. While for many numerical simulations in geotechnics use of a two‐phase soil model is sufficient, the simulation of compressed air support demands the use of a three‐phase model with the consideration of air as a separate phase. A multiphase model for soft soils is developed, in which the individual constituents of the soil—the soil skeleton, the fluid and the gaseous phase—and their interactions are considered. The triphasic model is formulated within the framework of the theory of porous media, based upon balance equations and constitutive relations for the soil constituents and their mixture. An elasto‐plastic, cam–clay type model is extended to partially saturated soil conditions by incorporating capillary pressure according to the Barcelona basic model. The hydraulic properties of the soil are described via DARCY 's law and the soil–water characteristic curve after VAN GENUCHTEN . Water is modelled as an incompressible and air as a compressible phase. The model is validated by means of selected benchmark problems. The applicability of the model to geotechnical problems is demonstrated by results from the simulation of a compressed air intervention in shield tunnelling. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of a high-cycle accumulation model to the analysis of soil liquefaction around a vibrating pile toe

High residual pore pressure observed in the vicinity of piles driven in saturated soil indicates that the soil around the pile may be liquefied. In the present paper, the problem of deformation of saturated sand around a vibrating pile is formulated with the use of a high-cycle accumulation model capable of describing a large number of cycles. The problem is solved numerically for locally undrained conditions in spherically symmetric formulation suitable for the lower part of a cylindrical closed-ended pile near the toe. The aim of the study is to calculate the evolution of the liquefaction zone around the pile for a large number of cycles. A parametric study is carried out to show how the growth of the liquefaction zone depends on the pile displacement amplitude, the relative soil density, the effective stress in the far field and the pore fluid compressibility.     

Two-phase Material Point Method applied to the study of cone penetration

This paper presents numerical simulations of Cone Penetration Test (CPT) in water-saturated soft soils taking into account pore pressure dissipation during installation. Besides modelling interaction between soil skeleton and pore fluid, the problem involves large soil deformations in the vicinity of the penetrometer, soil–structure interaction, and complex non-linear response of soil. This makes such simulations challenging. Depending on the soil’s permeability and compressibility, undrained, partially drained or drained conditions might occur. Partially drained conditions are commonly encountered in soils such as silts and sand–clay mixtures. However, this is often neglected in CPT interpretation, which may lead to inaccurate estimates of soil properties. This paper aims at improving the understanding of the penetration process in different drainage conditions through advanced numerical analyses. A two-phase Material Point Method is applied to simulate large soil deformations and generation and dissipation of excess pore pressures during penetration. The constitutive behaviour of soil is modelled with the Modified Cam Clay model. Numerical results are compared with experimental data showing good agreement.     

Hydro‐mechanical modeling of cohesive crack propagation in multiphase porous media using the extended finite element method

In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd.     

Two formulations for dynamic response of a cylindrical cavity in cross‐anisotropic porous media

Two formulations for calculating dynamic response of a cylindrical cavity in cross‐anisotropic porous media based on complex functions theory are presented. The basis of the method is the solution of Biot's consolidation equations in the complex plane. Employing two groups of potential functions for solid skeleton and pore fluid (each group includes three functions), the u–w formulation of Biot's equations are solved. Difference of these two solutions refers to use of two various potential functions. Equations for calculating stress, displacement and pore pressure fields of the medium are mentioned based on each two formulations. Copyright © 2009 John Wiley & Sons, Ltd.     

An analytical study of the effect of penetration rate on piezocone tests in clay

Finite element cavity expansion analysis investigating the effect of penetration rate on piezocone tests in clay is presented. A coupled analysis was performed, in which the rate of cavity expansion was linked to the penetration rate of the cone and the cone angle, using the assumption that the deformation was wholly radial, and took place only between the cone tip and the cone shoulder. The soil was modelled using modified cam clay with two sets of parameters and varying values of overconsolidation ratio (OCR). The influence of penetration rate on the stress and pore pressure distributions was examined. For slower penetration rates, the excess pore pressure at the cone shoulder is lower since consolidation is permitted coincident with penetration. The radial profiles of post‐penetration voids ratio demonstrate that partially drained penetration is permitted by volume change in the near field, in addition to radial movement in the far field. The radial distribution of excess pore pressure after slow penetration differs from the undrained case, with a relatively low radial gradient existing at the cone face. As a result, the dissipation curves after slow penetration lag behind those following fast penetration. The cone velocity is made dimensionless by normalizing with the coefficient of consolidation and the cone diameter. ‘Backbone’ curves of normalized velocity against normalized tip resistance and excess pore pressure capturing the transition from undrained to drained penetration are derived. The normalized pore pressure backbone curve is unique, whilst the normalized tip resistance shows a small dependency on OCR. These backbone penetration curves are compared with centrifuge model piezocone tests conducted at varying rates, and subsequent dissipation tests. The numerical and experimental results suggest that the value of consolidation coefficient operative during the dissipation phase is 2–4 times higher than the virgin compression value due to changes in the operative soil stiffness, as demonstrated from the stress paths of individual soil elements. The use of multi‐rate penetration tests to deduce values of consolidation coefficient is discussed, in light of these differences. Copyright © 2005 John Wiley & Sons, Ltd.