Propagation of Rayleigh waves in fluid‐saturated non‐homogeneous soils with the graded solid skeleton distribution

The propagation characteristic of Rayleigh waves in a fluid‐saturated non‐homogeneous poroelastic half‐plane is addressed. Based on Biot's theory for fluid‐saturated media, which takes the inertia, fluid viscosity, mechanical coupling, compressibility of solid grains, and fluid into account, the dispersion equations of Rayleigh waves in fluid‐saturated non‐homogeneous soils/rocks are established. By considering the shear modulus of solid skeleton variation with depth exponentially, a small parameter, which reflects the relative change of shear modulus, is introduced. The asymptotic solution of the dispersion equation expressing the relationship between the phase velocity and wave number is obtained by using the perturbation method. In order to analyze the effects of non‐homogeneity on the propagation characteristic of Rayleigh waves, the variation of the phase velocity with the wave number is presented graphically and discussed through numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Hydro-micromechanical modeling of wave propagation in saturated granular crystals

Biot theory predicts wave velocities in a saturated granular medium using the pore geometry, viscosity, densities, and elastic moduli of the solid skeleton and pore fluid, neglecting the interaction between constituent particles and local flow, which becomes essential as the wavelength decreases. Here, a hydro-micromechanical model, for direct numerical simulations of wave propagation in saturated granular media, is implemented by two-way coupling the lattice Boltzmann method (LBM) and the discrete element method (DEM), which resolve the pore-scale hydrodynamics and intergranular behavior, respectively. The coupling scheme is benchmarked with the terminal velocity of a single sphere settling in a fluid. In order to mimic a small amplitude pressure wave entering a saturated granular medium, an oscillating pressure boundary on the fluid is implemented and benchmarked with the one-dimensional wave equation. The effects of input waveforms and frequencies on the dispersion relations in 3D saturated poroelastic media are investigated with granular face-centered-cubic crystals. Finally, the pressure and shear wave velocities predicted by the numerical model at various effective confining pressures are found to be in excellent agreement with Biot analytical solutions, including his prediction for slow compressional waves.     

Dynamic responses of a poroelastic half‐space from moving trains caused by vertical track irregularities

The vibrations of railway tracks on a poroelastic half‐space generated by moving trains are investigated through a vehicle–track–ground coupling model. The theoretical model incorporates a vehicle, a track, and a fully saturated poroelastic half‐space soil medium. The source of vibration excitation is divided into two components: the quasi‐static loads and the dynamic loads. The quasi‐static loads are related to the static component of the axle loads, whereas the dynamic loads are due to the dynamic wheel–rail interaction. A linear Hertizian contact spring is introduced between each wheelset and the rail to consider the dynamic loads. Biot's dynamic theory is used to characterize the poroelastic half‐space soil medium. Using the Fourier transform, the governing equations for the track–ground system are solved and the numerical results are presented for a single axle vehicle model. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium are investigated. In addition, the different roles of the moving axle loads and the roughness‐induced dynamic loads are identified. It is concluded that the vibration level of the free field off the track predicted by the poroelastic soil medium is smaller than that predicted by the elastic soil medium for vehicle speed below the Rayleigh wave speed of the poroelastic half‐space, whereas it is larger for vehicle speed above the Rayleigh wave speed. The dynamic loads play an important role in the dynamic responses of the track–ground system. Copyright © 2010 John Wiley & Sons, Ltd.     

横观各向同性饱和地基三维瞬态Lamb问题

研究了横观各向同性饱和土地基在地表动力荷载作用下的三维瞬态响应。基于饱和多孔介质的三维Biot波动理论,利用Laplace变换,建立圆柱坐标系下横观各向同性饱和土的波动方程;解耦波动方程后,根据算子理论,并借助Fourier展开和Hankel变换技术,得到瞬态荷载作用下,饱和土介质的土骨架位移和应力、孔隙水相对位移和孔隙水压力的一般解;利用一般解,给出横观各向同性饱和地基在地表集中荷载激励下的瞬态Lamb问题的解答。数值算例结果表明,采用各向同性饱和介质的动力学模型,不能准确描述具有明显各向异性特性的饱和土地基的瞬态动力特性。     

淤砂层模拟方法对上层理想流体动力反应的影响——平面P波入射

采用解析法,以弹性半空间-淤砂层-理想流体层为对象,分别将淤砂层模拟为固体介质和流体介质,分析了平面P波入射时上层理想流体的动力反应特性。通过算例分析,比较了将淤砂层模拟为弹性固体、黏弹性固体、理想流体和黏性流体介质时与将淤砂层模拟为两相多孔介质时计算结果的差异。研究表明当渗透系数很小时,多孔介质可以近似简化为黏弹性固体,甚至线弹性固体。这一分析研究对于在高坝-库水-淤砂-地基系统地震反应分析的数值模型中采用合理的淤砂层模型,提高计算效率具有重要参考意义,而且还可以用于校核数值模型的合理性。     

Wave propagation in thermoelastic saturated porous medium

Biot’s theory for wave propagation in saturated porous solid is modified to study the propagation of thermoelastic waves in poroelastic medium. Propagation of plane harmonic waves is considered in isotropic poroelastic medium. Relations are derived among the wave-induced temperature in the medium and the displacements of fluid and solid particles. Christoffel equations obtained are modified with the thermal as well as thermoelastic coupling parameters. These equations explain the existence and propagation of four waves in the medium. Three of the waves are attenuating longitudinal waves and one is a non-attenuating transverse wave. Thermal properties of the medium have no effect on the transverse wave. The velocities and attenuation of the longitudinal waves are computed for a numerical model of liquid-saturated sandstone. Their variations with thermal as well as poroelastic parameters are exhibited through numerical examples.     

A new computational strategy for solving two‐phase flow in strongly heterogeneous poroelastic media of evolving scales

We develop a new computational methodology for solving two‐phase flow in highly heterogeneous porous media incorporating geomechanical coupling subject to uncertainty in the poromechanical parameters. Within the framework of a staggered‐in‐time coupling algorithm, the numerical method proposed herein relies on a Petrov–Galerkin postprocessing approach projected on the Raviart–Thomas space to compute the Darcy velocity of the mixture in conjunction with a locally conservative higher order finite volume discretization of the nonlinear transport equation for the saturation and an operator splitting procedure based on the difference in the time‐scales of transport and geomechanics to compute the effects of transient porosity upon saturation. Notable features of the numerical modeling proposed herein are the local conservation properties inherited by the discrete fluxes that are crucial to correctly capture the fingering patterns arising from the interaction between heterogeneity and nonlinear viscous coupling. Water flooding in a poroelastic formation subject to an overburden is simulated with the geology characterized by multiscale self‐similar permeability and Young modulus random fields with power‐law covariance structure. Statistical moments of the poromechanical unknowns are computed within the framework of a high‐resolution Monte Carlo method. Numerical results illustrate the necessity of adopting locally conservative schemes to obtain reliable predictions of secondary recovery and finger growth in strongly heterogeneous deformable reservoirs. Copyright © 2011 John Wiley & Sons, Ltd.     

横观各向同性饱和地基竖向振动的衰减特性

基于提出的横观各向同性饱和多孔介质Biot波动方程的一般解,研究了饱和半空间地基在竖向点源简谐激振荷载作用下地表振动的衰减特征,分析了激振频率以及横观各向同性饱和土介质的各向异性参数和孔隙渗透系数对地表振动特征的影响。计算结果表明,低频和高频激振时,地表位移衰减特性存在明显差异;在饱和土的各向异性参数中,纵向和水平方向动态渗透系数比值和刚度系数比值对地表位移衰减影响最大,这也说明采用各向同性饱和介质的动力学模型不能准确地描述具有明显各向异性特性的饱和土地基的动力特性。     

Seismic wave propagation in laterally inhomogeneous poroelastic media via BIEM

This work addresses in‐plane pressure P and vertically polarized shear SV seismic wave propagation in a finite, laterally inhomogeneous, multilayered poroelastic geological region resting on the homogeneous elastic half‐space. The particular approach followed here is based on a combination of the (i) viscoelastic approximation (isomorphism) to Biot's equations of dynamic poroelasticity and on the (ii) boundary integral equation method (BIEM) using frequency‐dependent fundamental solutions of the governing wave equations. The problem is formulated under plane strain conditions and time‐harmonic motions are assumed. Validation of the viscoelastic isomorphism and verification of the BIEM is done by solution of benchmark examples. These simulation studies reveal that the proposed methodology is able to depict a sensitivity of the seismic signals recovered to the following parameters: (i) poroelastic properties of fluid saturated layers; (ii) lateral geological inhomogeneity; (iii) surface topography and (iv) frequency content and direction of the incident wave. It is concluded that the combination of viscoelastic isomorphism with BIEM software provides an effective numerical tool for evaluating site‐effect phenomena in multilayered, fluid saturated geological regions with complex geometry. The numerical results obtained demonstrate that dynamic poroelasticity interacting with other physical peculiarities of the Earth's surface layers, such as lateral heterogeneity, material properties along the wave path, local geological profile and type of elastic wave, gives rise to complex seismic signals on the free surface at the site of interest. Copyright © 2010 John Wiley & Sons, Ltd.     

Poroelastic behaviour of saturated brittle rock with anisotropic damage

A new anisotropic poroelastic damage model is proposed for saturated brittle porous materials. The model is formulated in the framework of the continuum damage mechanics. A second‐rank symmetric tensor is used to characterize material damage due to oriented microcracks. The classic Biot poroelastic theory is then extended to include poroelastic damage coupling. Both the deterioration of elastic properties and poroelastic coefficients is taken into account. A suitable procedure for determination of model parameters from standard laboratory tests is presented. The validity of the model is tested through comparison between numerical predictions and experimental data in various loading conditions. The overall performance of the model is evaluated. The choice of relevant effective stress for the microcrack propagation criterion in saturated cohesive geomaterials is discussed. Copyright © 2000 John Wiley & Sons, Ltd.     

半空间饱和土中圆形衬砌对弹性压缩波的散射

借助于Biot 波动理论和弹性波的传播理论,采用复变函数和多级坐标法,对半空间饱和土中圆形衬砌结构对弹性稳态压缩波的散射问题进行求解和分析。利用一个半径很大的圆弧来逼近半空间直边界,将待解问题转化为稳态弹性压缩波在一个大圆孔和一个弹性衬砌结构的散射问题。通过引入势函数,将饱和土的Biot波动方程和衬砌的弹性波动方程解耦成Helmholtz 方程,借助复变函数级数展开便可以预先写出该组Helmholtz方程的通解。然后,通过引用复变量,把饱和土和衬砌结构中的应力、位移及孔压用设定的势函数表示出来,再利用半空间饱和土和衬砌结构的连续性条件和近似直边界的圆弧边界和衬砌内边界的边界条件求解出该组势函数的特解。最后,利用势函数的特解,得到饱和土中的位移,应力和孔压及衬砌结构的位移和应力;变换不同的参数求解衬砌结构内外边界的动应力和孔压的集中系数,通过对算例结果的分析得出一系列有益的结论。     

爆轰荷载作用下球空腔热流固耦合动力响应

基于饱和多孔弹性介质热流固耦合动力模型(THMD)及控制方程,研究了爆轰荷载作用下球空腔内壁面受到随时间变化热、力冲击作用下的热流固耦合动力响应。利用Laplace变换技术,对控制方程进行解耦,获得了温度、位移、孔隙水压力和应力积分形式解。采用Laplace逆变换得到数值结果,分析了热流固耦合参数的影响特性,并将热流固耦合结果与热弹性条件下的结果进行比较,以验证热流固耦合条件下结果的正确性。     

A general viscous‐spring transmitting boundary for dynamic analysis of saturated poroelastic media

A time‐domain viscous‐spring transmitting boundary is presented for transient dynamic analysis of saturated poroelastic media with linear elastic and isotropic properties. The u–U formulation of Biot equation in cylindrical coordinate is adopted in the derivation. By this general viscous‐spring boundary, the effective stress and pore fluid pressure on the truncated boundary of the computational area are replaced by a set of continuously distributed spring and dashpot elements, of which the parameters are defined assuming an infinite permeability and considering the two dilatational waves. Numerical examples demonstrate good absorption of both the two cylindrical dilatational waves by the proposed ‘drained’ boundary. For general two‐dimensional wave propagation problems, acceptable accuracy can still be achieved by setting the proposed boundary relatively far away from the scatter. Numerical comparison shows that the results obtained by using this boundary are more accurate for all permeability values than those by the traditional viscous‐spring or viscous boundaries established for u–U formulation. Copyright © 2015 John Wiley & Sons, Ltd.     

Three‐dimensional time‐harmonic fundamental solutions for a fluid‐saturated poroelastic half‐space with partially permeable free surface

By virtue of a pair of scalar potentials for the displacement of the solid skeleton and the pore fluid pressure field of a saturated poroelastic medium, an alternative solution method to the Helmholtz decomposition is developed for the wave propagation problems in the framework of Biot's theory. As an application, a comprehensive solution for three‐dimensional response of an isotropic poroelastic half‐space with a partially permeable hydraulic free surface under an arbitrarily distributed time‐harmonic internal force field and fluid sources is developed. The Green's functions for the poroelastic fields, corresponding to point, ring, and disk loads, are reduced to semi‐infinite complex‐valued integrals that can be evaluated numerically by an appropriate quadrature scheme. Analytical and numerical comparisons are made with existing elastic and poroelastic solutions to illustrate the quality and features of the solution. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamic responses of unsaturated half‐space soil to a moving harmonic rectangular load

The problem of the dynamic responses of a semi‐infinite unsaturated poroelastic medium subjected to a moving rectangular load is investigated analytical/numerically. The dynamic governing equations are obtained with consideration of the compressibility of solid grain and pore fluid, inertial coupling, and viscous drag as well as capillary pressure in the unsaturated soil, and they can be easily degraded to the complete Biot's theory. Using the Fourier transform, the general solution for the equations is derived in the transformed domain, and then a corresponding boundary value problem is formulated. By introducing fast Fourier transform algorithm, the unsaturated soil vertical displacements, effective stresses, and pore pressures induced by moving load are computed, and some of the calculated results are compared with those for the degenerated solution of saturated soils and confirmed. The influences of the saturation, the load speed, and excitation frequency on the response of the unsaturated half‐space soil are investigated. The numerical results reveal that the effects of these parameters on the dynamic response of the unsaturated soil are significant.     

移动环形荷载作用下饱和土中圆形衬砌隧洞动力响应研究

基于Biot理论,研究了饱和土中带有衬砌的圆形隧洞在移动环形荷载作用下的动力响应。假定衬砌为弹性体,土体为饱和多孔介质,引入两类势函数来表示土体、孔隙水和衬砌的位移,使隧洞的控制方程解耦。结合边界条件及连续条件,通过傅立叶变换得到频率-波数域中衬砌和土体的应力、位移和孔隙水压力解答,最后用傅立叶积分逆变换得到时-空域中的数值解。计算并比较了3种隧洞模型（弹性土体隧洞、饱和土体隧洞和饱和土衬砌隧洞）的动力响应分析。数值分析结果说明：（1）移动荷载速度对3种隧洞动力响应均具有较大影响;（2）弹性土体隧洞和饱和土体隧洞的动力响应具有明显区别,所以在富水地区的隧洞动力响应中土体应该视为饱和土体;（3）衬砌对隧洞动力响应有较大影响,故隧洞的动力分析中不能忽略衬砌作用。     

Theoretical and numerical analyses of chemical‐dissolution front instability in fluid‐saturated porous rocks

The chemical‐dissolution front propagation problem exists ubiquitously in many scientific and engineering fields. To solve this problem, it is necessary to deal with a coupled system between porosity, pore‐fluid pressure and reactive chemical‐species transport in fluid‐saturated porous media. Because there was confusion between the average linear velocity and the Darcy velocity in the previous study, the governing equations and related solutions of the problem are re‐derived to correct this confusion in this paper. Owing to the morphological instability of a chemical‐dissolution front, a numerical procedure, which is a combination of the finite element and finite difference methods, is also proposed to solve this problem. In order to verify the proposed numerical procedure, a set of analytical solutions has been derived for a benchmark problem under a special condition where the ratio of the equilibrium concentration to the solid molar density of the concerned chemical species is very small. Not only can the derived analytical solutions be used to verify any numerical method before it is used to solve this kind of chemical‐dissolution front propagation problem but they can also be used to understand the fundamental mechanisms behind the morphological instability of a chemical‐dissolution front during its propagation within fluid‐saturated porous media. The related numerical examples have demonstrated the usefulness and applicability of the proposed numerical procedure for dealing with the chemical‐dissolution front instability problem within a fluid‐saturated porous medium. Copyright © 2007 John Wiley & Sons, Ltd.     

考虑热-水-力耦合效应多孔弹性地基的动力响应

通过对Biot波动方程的修正,得到考虑热-水-力学耦合效应的多孔弹性介质动力响应的控制方程,研究了简谐均布荷载作用下地基土体的热-水-力耦合动力响应问题。利用Fourier变换技术,得到地基中的应力、位移和孔隙水压力积分形式的解答。利用Fourier逆变换得到数值结果,分析了热-水-力学耦合条件下地基土体中温度增量、应力、位移和孔隙水压力响应的分布,并讨论了热源输入的影响, 结果表明：应力、位移和孔隙水压力随 的增大而有一定的减小。     

Instability of wave propagation in saturated poroelastoplastic media

In the present work, stationary discontinuities and fluttery instabilities of wave propagation in saturated poro‐elastoplastic media are analysed in the frame of Biot theory. The generalized Biot formulations are particularly employed for simulating non‐linear coupled hydro‐mechanical behaviour of the media. Inertial coupling effect between the solid and the fluid phases of the media is also taken into account. The non‐associated Drucker–Prager criterion to describe non‐linear constitutive behaviour of pressure dependent elasto‐plasticity for the solid skeleton of the media is particularly considered. With omission of compressibility of solid grains and the pore fluid, the critical conditions of stationary discontinuities and flutter instabilities occurring in wave propagation are given in explicit forms. It is shown that when the stationary discontinuity is triggered at the surface of discontinuity there still may exist real wave speeds. The wave speeds across the stationary discontinuity surface entirely cease to be real only in non‐associated plasticity, certain ranges of value of Poisson's ratio and when compression stress normal to the surface of discontinuity dominates the stress state at the surface. It is also indicated that the fluttery instabilities, under which some wave speeds cease to be real even in strain hardening stage, may occur prior to stationary discontinuities only for non‐associated plasticity under certain conditions. These conditions are: (1) both the porosity and the Poisson's ratio possess relatively low values and (2) the deviatoric part of the effective stress normal to the surface of discontinuity is compressive. A region in the porosity–Poisson's ratio plot, in which fluttery instabilities are possible to occur, is given. Copyright © 2002 John Wiley & Sons, Ltd.     

Analytical solutions of a finite two‐dimensional fluid‐saturated poroelastic medium with compressible constituents

An analytical solution is proposed for transient flow and deformation coupling of a fluid‐saturated poroelastic medium within a finite two‐dimensional (2‐D) rectangular domain. In this study, the porous medium is assumed to be isotropic, homogeneous, and compressible. In addition, the point sink can be located at an arbitrary position in the porous medium. The fluid–solid interaction in porous media is governed by the general Biot's consolidation theory. The method of integral transforms is applied in the analytical formulation of closed‐form solutions. The proposed analytical solution is then verified against both exact and numerical results. The analytical solution is first simplified and validated by comparison with an existing exact solution for the uncoupled problem. Then, a case study for pumping from a confined aquifer is performed. The consistency between the numerical solution and the analytical solution confirms the accuracy and reliability of the analytical solution presented in this paper. The proposed analytical solution can help us to obtain in‐depth insights into time‐dependent mechanical behavior due to fluid withdrawal within finite 2‐D porous media. Moreover, it can also be of great significance to calibrate numerical solutions in plane strain poroelasticity and to formulate relevant industry norms and standards. Copyright © 2014 John Wiley & Sons, Ltd.