空心圆柱饱和多孔介质热固结问题的解析解

基于考虑热渗效应和等温热流效应的热-水-力耦合的线性热弹性固结控制方程,建立无限长空心圆柱饱和多孔介质热固结问题的一种理论求解方法。该方法先给出Laplace变换域上的解,然后,利用Stehfest法求其数值逆变换。该理论解考虑了空心圆柱体内、外透水界面随时间变化的外力和温度荷载耦合作用过程。最后,通过一算例分析了饱和多孔介质的热固结特征,给出其温度、孔压、位移和应力的演化规律     

Analytical study of freezing behavior of a cavity in thermo-poro-elastic medium

The goal of this paper is to present an analytical solution to have a first insight of the impact of ice formation on the surrounding porous rock on underground cavities like reservoir, pipes, tunnels or wellbores. Among the other analytical solutions found in the literature on this topic, the originality of this work resides in the rigorous theoretical framework of poromechanics, which considers the coupling between liquid water and ice crystal under thermodynamic equilibrium. Liquid water transport, thermal conduction, and elastic properties of the phases are also considered. Two analytical solutions are presented, based on a linearization of the system of governing equations. The first one deals with a spherical cavity within an infinite porous medium leading to an exact analytical solution. It allows validating the Stehfest’s algorithm on the numerical inversion of Laplace Transform, used in the second analytical solution, which considers a cylindrical excavation. The validity of this solution is assessed by comparing its results to that issued from a numerical resolution of the nonlinear system of equations. The analytical solution is then ultimately used to identify the influence of key parameters like the thermal/hydraulic conductivities, the amount of ice formed and the thermal dilatation coefficients on the mechanical response of a cylindrical cavity submitted to an internal frost.     

Advective transport from a penny-shaped crack in a porous medium and an associated uniqueness theorem

This paper examines the problem of the advective transport of a contaminant from sources in the shape of either a penny-shaped crack or an elongated needle-shaped cavity located in a porous medium of infinite extent. The advective transport is induced by Darcy flow in the porous medium, where the internal boundary is maintained at a constant potential. The paper presents an approximate analytical solution to this problem, which is deduced from a formulation that models a cavity in the shape of either an oblate or a prolate spheroid. The results also represent one of the few spatially three-dimensional exact analytical solutions for the, albeit linear, hyperbolic problem governing the contaminant transport problem. The paper also presents a canonical proof of uniqueness for advective contaminant transport problems associated with media of infinite extent. Copyright © 2004 John Wiley & Sons, Ltd.     

Exact solutions for one‐dimensional consolidation of single‐layer unsaturated soil

Based on the Fredlund consolidation theory of unsaturated soil, exact solutions of the governing equations for one‐dimensional consolidation of single‐layer unsaturated soil are presented, in which the water permeability and air transmission are assumed to be constants. The general solution of two coupled homogeneous governing equations is first obtained. This general solution is expressed in terms of two functions psi₁ and ψ₂, where ψ₁ and ψ₂, respectively, satisfy two second‐order partial differential equations, which are in the same form. Using the method of separation of variables, the two partial differential equations are solved and exact solutions for three typical homogeneous boundary conditions are obtained. To obtain exact solutions of nonhomogeneous governing equations with three typical nonhomogeneous boundary conditions, the nonhomogeneous boundary conditions are first transformed into homogeneous boundary conditions. Then according to the method of undetermined coefficients and exact solutions of homogenous governing equations, the series form exact solutions are put forward. The validity of the proposed exact solutions is verified against other analytical solutions in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Thermomechanical response of a poroelastic half-space soil medium subjected to time harmonic loads

The thermomechanical responses of a porous elastic medium subjected to time harmonic loads (normal force and thermal source) are investigated analytically in the context of generalised thermoelastic theory with one relaxation time. The material of the foundation, obeying Biot’s dynamic poroelastic theory, is idealised as a uniform, fully saturated poroelastic half-space stratum. The coupled governing equations are established based on Biot’s dynamic poroelastic theory and on generalised thermoelastic theory. Assuming the disturbances to be harmonically time dependent, the general solutions of stress, displacement, temperature distribution and excess pore water pressure are deduced using the Fourier transform, and the transformed solutions are numerically inverted. The differences among the coupled thermo-hydro-mechanical dynamic model (THMD), the hydro-mechanical dynamic model (HMD) and the thermo-elastic dynamic model (TMD) are discussed. In addition, the effects of the thermal loading frequency on the displacement, stress, temperature distribution and excess pore water pressure components are analysed in the numerical results.     

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.     

Mode of a spherical cavity’s thermo-elastodynamic response in a saturated porous medium for non-torsional loads

The thermo-hydro-elastodynamic model (THED), in which the thermo-osmosis and thermal-filtration phenomenon in a two-phase porous thermoelastic medium can be considered, was previously presented by the authors [Ganbin Liu, Kanghe Xie, Rongyue Zheng. Model of nonlinear coupled thermo-hydro-elastodynamics response for a saturated poroelastic medium. Sci China (Ser E) 2009;52(8):2373–83] and is used in this paper to investigate the thermo-elastodynamic response of a spherical cavity in a saturated poroelastic medium when subjected to a time-dependent non-torsional thermal/mechanical source. The Separated Variable Method is introduced, and the non-zero displacement potentials are expanded in terms of the Legendre polynomial. Solutions for the displacement, temperature increment, pore pressure and stress are obtained in the domain of the Laplace transform. Numerical results are also performed for different modes and are compared with the results of the thermoelastic model (TED) to ascertain the validity and the difference between these two models.     

Response of saturated and nearly saturated porous media: Different formulations and their applicability

The response of saturated porous medium is of significant interest in many fields ranging from geomechanics to biomechanics. Biot was the first to formulate the basic equations governing the process of coupled flow and deformation in porous media. Depending on the nature of loading vis‐à‐vis the characteristics of the media, different formulations (fully dynamic, partly dynamic, quasi‐static) are possible. In this study, analytical solutions are developed for the response of saturated and nearly saturated porous media under plane strain condition. The solutions for different formulations are developed in terms of non‐dimensional parameters. The response is studied for various conditions and the regions of validity for various formulations are identified in a parametric space. An assessment of the needed formulation for few important problems is also presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Modeling heat and mass transfer during ground freezing subjected to high seepage velocities

Natural or induced groundwater flow may negatively influence the performance of artificial ground freezing: high water flow velocities can prevent frozen conditions from developing. Reliable models that take into consideration hydraulic mechanisms are then needed to predict the ground freezing development. For forty years, numerous thermo-hydraulic coupled numerical models have been developed. Some of these models have been validated against experimental data but only one has been tested under high water flow velocity conditions. This paper describes a coupled thermo-hydraulic numerical model completely thermodynamically consistent and designed to simulate artificial ground freezing of a saturated and non-deformable porous medium under seepage flow conditions. On some points, less restrictive assumptions than the ones usually used in the literature are considered. As for the constant-porosity assumption, its validity is verified. The model appears to be well validated against analytical solutions and a three-dimensional ground freezing experiment under high seepage flow velocity conditions. It is used to highlight key thermo-hydraulic mechanisms associated with phase change in a porous medium.     

饱和多孔介质热-水-力控制方程耦合项的意义及耦合影响分析

基于饱和多孔介质完全耦合条件下热-水-力控制方程的一般形式,对各耦合项的物理意义进行了分析。以等温条件下水-力耦合的Biot渗透固结理论和等压条件下热-力耦合的热弹性理论两种退化形式为基础,给出耦合系数的物理意义及其函数表达式。研究表明,其中一些耦合系数实际上有相同的物理意义和数学形式,因而可使控制方程得到简化。针对不同的介质特性以及热、水、力条件,给出几种特定情形下控制方程的简化形式。以无限大多孔介质中球域热源的热力学响应分析为例,对一些耦合项的影响程度进行了比较研究。     

Numerical modeling of multiphase fluid flow in deforming porous media: A comparison between two- and three-phase models for seismic analysis of earth and rockfill dams

In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior.     

Poroelastic model for pile–soil interaction in a half‐space porous medium due to seismic waves

In this paper, frequency domain dynamic response of a pile embedded in a half‐space porous medium and subjected to P, SV seismic waves is investigated. According to the fictitious pile methodology, the problem is decomposed into an extended poroelastic half‐space and a fictitious pile. The extended porous half‐space is described by Biot's theory, while the fictitious pile is treated as a bar and a beam and described by the conventional 1‐D structure vibration theory. Using the Hankel transformation method, the fundamental solutions for a half‐space porous medium subjected to a vertical or a horizontal circular patch load are established. Based on the obtained fundamental solutions and free wave fields, the second kind of Fredholm integral equations describing the vertical and the horizontal interaction between the pile and the poroelastic half‐space are established. Solution of the integral equations yields the dynamic response of the pile to plane P, SV waves. Numerical results show the parameters of the porous medium, the pile and incident waves have direct influences on the dynamic response of the pile–half‐space system. Significant differences between conventional single‐phase elastic model and the poroelastic model for the surrounding medium of the pile are found. Copyright © 2007 John Wiley & Sons, Ltd.     

A semianalytical solution for one-dimensional transient wave propagation in a saturated single-layer porous medium with a fluid surface layer

One-dimensional transient wave propagation in a saturated single-layer porous medium with a fluid surface layer is studied in this paper. An analytical solution for a special case with a dynamic permeability coefficient k_f → ∞ and a semianalytical solution for a general case with an arbitrary dynamic permeability coefficient are presented. The eigenfunction expansion and precise time step integration methods are employed. The solution is presented in series form, and thus, the long-term dynamic responses of saturated porous media with small permeability coefficients can be easily computed. We first transform the nonhomogeneous boundary conditions into homogeneous boundary conditions, and then we obtain the eigenvalues and orthogonal eigenfunctions of the fluid–solid system. Finally, the solutions in the time domain are developed. As the model is one dimensional, geometric attenuation is absent, and only the attenuation in the saturated porous medium is considered. We can apply this model to analyse the influences of different seabed types on the propagation of acoustic waves in the fluid layer, which is very important in ocean acoustics and ocean seismic. This solution can also be employed to validate the accuracies of various numerical methods.     

Thermal effects on fluid flow and hydraulic fracturing from wellbores and cavities in low-permeability formations

The coupled heat-fluid-stress problem of circular wellbore or spherical cavity subjected to a constant temperature change and a constant fluid flow rate is considered. Transient analytical solutions for temperature, pore pressure and stress are developed by coupling conductive heat transfer with Darcy fluid flow in a poroelastic medium. They are applicable to low permeability porous media suitable for liquid-waste disposal and also simulating reservoir for enhanced oil recovery, where conduction dominates the heat transfer process. A full range of solutions is presented showing separately the effects of temperature and fluid flow on pore pressure and stress development. It is shown that injection of warm fluid can be used to restrict fracture development around wellbores and cavities and generally to optimize a fluid injection operation. Both the limitations of the solutions and the convective flow effect are addressed. Copyright © 1999 John Wiley & Sons, Ltd.     

Strain‐softening analysis of a spherical cavity

This paper studies the excavation of a spherical cavity subjected to hydrostatic initial stresses in the infinite homogeneous and isotropic rock mass with strain‐softening Mohr–Coulomb (M‐C) and Hoek–Brown (H‐B) behaviors. Numerical solutions of the spherical cavity are obtained and the application to determining stress–strain curve of strain‐softening M‐C and H‐B rock mass is studied. A closed‐form solution for the elastic–brittle–plastic medium is introduced first, and then a numerical procedure that simplifies the strain‐softening process into a series of brittle–plastic ones is presented. The approach is validated against the facts that the strain‐softening process evolves into a brittle–plastic one when the softening slope is very steep, whereas it evolves into an elasto‐plastic one when the softening slope approaches zero. Numerical solutions for the prediction of displacements and stresses around the spherical cavity in the strain‐softening M‐C and H‐B rock mass are presented. On the basis of the analysis of the spherical cavity in strain‐softening rock mass, the stress–strain relationship at an infinitesimal cube around the cavity is obtained and discussed with different evolution laws for the strength parameters considered. Copyright © 2011 John Wiley & Sons, Ltd.     

A solution around a backfilled cavity in a low‐permeability poroelastic medium with application in in situ heating tests

Thermo‐hydro‐mechanical responses around a cylindrical cavity drilled or excavated in a low‐permeability formation are studied when the cavity is subjected to a time‐dependent thermal loading. The cavity is considered backfilled after it is supported by casing or lining. Solutions of temperature, pore water pressure, stress, and displacement responses are analytically formulated based on Biot's consolidation theory with the assumption that the backfilling material, supporting material, and surrounding low‐permeability formation are poroelastic media. The solution is expressed in Laplace space, and numerical inversion techniques are used to find field variables in the real‐time domain. After the solution is verified with the numerical results, it is applied in a large‐scale in situ heating test – PRACLAY heating test – for a predictive reference calculation and an extensive parametric study. Another medium‐scale in situ heating test – ATLAS III heating test – is also analyzed using the solution, which provides reasonable agreement with measurements. The new analytical solution proves to be a convenient tool for a good understanding of the resulting coupled thermo‐hydro‐mechanical behavior and is therefore valuable for the interpretation of measured data in engineering practices and for a rational design of potential radioactive waste repositories. Copyright © 2016 John Wiley & Sons, Ltd.     

The behaviour of an elastic non-homogeneous half-space. Part I-line and point loads

In the first part of this paper solutions are developed for the response of a non-homogeneous half-space subjected to either a surface point load or a surface line load. The non-homogeneity considered is a variation in Young's modulus (E) with depth (z) which takes the form E=m_EZ^α where m_E is a constant and α is referred to as the non-homogeneity parameter. The variation of these solutions as the non-homogeneity parameter α varies between the limits of zero (homogeneous soil) to unity (Gibson soil) gives some fresh insight into both these limiting cases.     

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

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

A finite element method for localized erosion in porous media with applications to backward piping in levees

This paper presents a computational method able to effectively model both the simultaneous processes typically observed in backward erosion piping, ie, the pipe tip propagation and the conduit cross-section enlargement. The numerical method is based on the novel formulation of a problem of localized erosion along a line propagating in a multidimensional porous medium. In this line, a conduit with evolving transverse size is embedded, which conveys a multiphase flow. The two systems, porous medium and pipe, are bridged by exchange terms of multiphase fluid mass and by a shared fluid pressure field. On the contrary, different fields are considered to describe flows, which are assumed as Darcian in the porous medium and turbulent in the conduit. These two flows drive pipe propagation and enlargement, respectively, as modeled by means of proper erosion kinetic laws. The corresponding numerical formulation is based on the combination between one- and multidimensional finite elements, to model the erosion conduit and the porous medium, respectively. Several simulations are proposed to demonstrate the ability of the proposed approach in reproducing available experimental data of real-scale tests on levees. Our results point out the crucial role played by the combined influence of pipe propagation and enlargement, as well as of three-dimensional (3D) effects. We also assess the mesh independence of the proposed numerical solution, particularly as concerns the calculated pipe propagation history.