非饱和土中热弹性波的传播特性分析

基于非饱和多孔介质的研究成果,考虑热效应和孔隙流体迂曲度的影响,研究了非饱和土中热弹性波的传播特性。利用非饱和土中耦合热的固-液-气三相介质的质量平衡方程、渗流连续方程、动量平衡方程和广义非Fourier热传导定律,建立了问题的热弹性波动方程。通过引入势函数,经过理论推导给出了非饱和土中热弹性波的弥散特征方程。结合数值算例,分析了几类热弹性波的波速和衰减系数随迂曲度、热膨胀系数和介质温度等热物理参数的变化规律。结果表明：孔隙水迂曲度的增大将引起P1波、P3波和S波的波速增大,而孔隙气体迂曲度的增大仅使得P2波的波速增大;热膨胀系数的增大将造成P1波波速的增大和热（T）波波速的减小;介质温度的升高将引起各类热弹性波波速的增大;频率、热膨胀系数和介质温度的变化对各类热弹性波的衰减系数均有较大影响,不可忽视。     

Wave propagation in a general anisotropic poroelastic medium: Biot’s theories and homogenisation theory

Anisotropic wave propagation is studied in a fluid-saturated porous medium, using two different approaches. One is the dynamic approach of Biot’s theories. The other approach known as homogenisation theory, is based on the averaging process to derive macroscopic equations from the microscopic equations of motion. The medium considered is a general anisotropic poroelastic (APE) solid with a viscous fluid saturating its pores of anisotropic permeability. The wave propagation phenomenon in a saturated porous medium is explained through two relations. One defines modified Christoffel equations for the propagation of plane harmonic waves in the medium. The other defines a matrix to relate the relative displacement of fluid particles to the displacement of solid particles. The modified Christoffel equations are solved further to get a quartic equation whose roots represent complex velocities of the four attenuating quasi-waves in the medium. These complex velocities define the phase velocities of propagation and quality factors for attenuation of all the quasi-waves propagating along a given phase direction in three-dimensional space. The derivations in the mathematical models from different theories are compared in order to work out the equivalence between them. The variations of phase velocities and attenuation factors with the direction of phase propagation are computed, for a realistic numerical model. Differences between the velocities and attenuations of quasi-waves from the two approaches are exhibited numerically.     

Polarisations of quasi-waves in a general anisotropic porous solid saturated with viscous fluid

Wave propagation is studied in a general anisotropic poroelastic solid saturated with a viscous fluid flowing through its pores of anisotropic permeability. The extended version of Biot’s theory is used to derive a system of modified Christoffel equations for the propagation of plane harmonic waves in such media. The non-trivial solution of this system is ensured by a biquadratic equation whose roots represent the complex velocities of four attenuating quasi-waves in the medium. These complex velocities define phase velocity and attenuation of each quasi-wave propagating along a given phase direction in three-dimensional space. The solution itself defines the polarisations of the quasi-waves along with phase shift. The variations of polarisations of quasi-waves with their phase direction, are computed for a realistic numerical model.     

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.     

饱和冻土中弹性体波传播特性影响参数研究

基于Leclaire对饱和双相孔隙弹性介质Biot模型的扩展,研究含有两种不同固相组分的三相多孔弹性介质中体波的传播特性。以饱和冻土为例,分析了各相体积分数、颗粒形状,接触参数等因素对波动方程中惯性参数、黏性参数、刚度参数的影响;对该三相介质模型进行了退化,分析了孔隙中只含液态水或固态冰时体波的特性;以饱和冻土为例,通过数值计算,探讨了饱和冻土中体波的相速度和衰减系数与胶结参数、接触参数、频率、饱和度、孔隙率等参数的关系。结果表明：与一般的饱和土不同,饱和冻土中存在5种体波,即3种纵波和2种横波;5种体波均具有弥散性和衰减性,且P1波、S1波弥散性和衰减性远小于P2、P3、S2波;胶结参数、饱和度、孔隙率对5种体波的传播特性影响显著,接触参数对传播特性影响较小。     

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.     

Propagation of Rayleigh waves in unsaturated porothermoelastic media

Based on the frame of elastic theory for unsaturated porous medium, considering the influence of thermal effect, the propagation characteristics of Rayleigh wave in unsaturated porous media are studied. Firstly, the thermoelastic wave equations for three-phase porous media are established, in which the mass balance equations, generalized Darcy law, momentum balance equations, and generalized non-Fourier heat conduction law are taken into account. Secondly, through theoretical derivation, the dispersion equation of Rayleigh wave for unsaturated porothermoelastic media is given by introducing the potential functions. Finally, the variations of the phase velocity of Rayleigh wave are analyzed with numerical examples. The results show that the thermal conductivity has little effect on the phase velocity of Rayleigh wave. The phase velocity of Rayleigh wave increase with increasing of the thermal expansion coefficient and media temperature.     

Surface load dynamic solution of saturated transversely isotropic multilayer half-space

This paper treats the dynamic response of a multilayered transversely isotropic fluid saturated poroelastic half-space under surface time-harmonic traction. The governing system of partial differential equations is uncoupled with the use of a set of physically meaningful and complete potential functions that decompose different body waves in a saturated poroelastic transversely isotropic medium. After expressing the equations in the Hankel-Fourier domain, a proper algebraic factorization is applied to generate reflection and transmission matrices for decomposed waves. All responses including displacements, stresses, and pore fluid pressure for both general patch load and point load are presented in the form of semi-infinite line integrals. The verification of the method is confirmed with the degeneration of the solutions presented here to the existing solutions for dried both homogeneous and multilayered elastic half-spaces as well as poroelastic half-space. Selected numerical results are depicted to investigate the effects of layering and pore pressure on responses of a transversely isotropic poroelastic medium. The load distribution effects are studied by comparison of the patch and point load responses. Also, resonance notion and effective parameters on this phenomenon such as layering system and anisotropy contrast are discussed. Significant influence of materials and layering configuration on number and amplitude of resonances depicted through the numerical evaluation.     

简谐线源荷载作用下热流固耦合地基的动力响应

基于Biot波动理论和广义热弹性理论,对简谐线源荷载（力荷载和热荷载）作用下的热-流-固耦合地基的动力响应问题进行了研究。将地基看成是均质各向同性、完全饱和的多孔半空间介质,利用无量纲化和Fourier变换方法对热-流-固耦合控制方程进行简化,得到了变换域内应力分量、位移分量、温度分布及超孔隙水压力的一般解,并利用Fourier逆变换得到了相应的积分形式解答。通过数值计算对按热-流-固耦合理论、饱和多孔弹性理论和热弹性理论得到的结果进行了比较,同时分析了热-流-固耦合条件下热荷载激振频率对竖向应力、竖向位移、温度分布以及超孔隙水压力的影响。     

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

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

Scattering of elastic wave by a partially sealed cylindrical shell embedded in poroelastic medium

Scattering of an elastic wave by a cylindrical shell embedded in poroelastic medium is investigated theoretically with the assumption that the shell material is also a porous elastic medium. The porous medium is modellized via Biot's theory and the scattering by cylindrical shell is expressed by the definition of scattering matrix. The normal mode expansion technique is employed for analysing the scattering field, and the asymptotic solutions of displacements, stresses and pore pressure are derived. Two limiting cases‐scattering by a poroelastic cylinder in Biot medium and a elastic cylindrical shell in elastic medium are obtained from the general solutions. The dispersion curves of displacement amplitude at the interface of shell and medium is compared with the case of elastic shell. The scattering amplitude associated with the fast, slow and transverse waves are identified by numerical simulation. Furthermore, the influence of the poroelastic property of shell material on scattering amplitude is analysed. Copyright © 2005 John Wiley & Sons, Ltd.     

Existence of longitudinal waves in pre-stressed anisotropic elastic medium

In a pre-stressed anisotropic elastic medium, three types of quasi-waves propagate along an arbitrary direction. In general, none of the waves is truly longitudinal. The present study finds the specific directions in a pre-stressed anisotropic elastic medium along which longitudinal waves may propagate. This paper demonstrates how the propagation of longitudinal waves is affected by various pre-stresses present in the medium. The study establishes the explicit expressions defining the existence and propagation of longitudinal waves in pre-stressed anisotropic elastic medium. These expressions involve not only the direction and elastic stiffness of the medium, but also the pre-stresses present in the medium. Changes in conditions for the existence of longitudinal waves in orthotropic, monoclinic and triclinic anisotropies are discussed in detail. The most important part of the paper is a practical aspect suggested to calculate the specific directions for the existence of longitudinal waves in pre-stressed anisotropic elastic medium. In this approach, only those parameters are used that can be observed by the receiver in a geophysical experiment of wave propagation. The existence of longitudinal waves has been shown graphically using a numerical example for three types of anisotropic symmetries in elastic medium.     

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.     

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.     

爆破荷载应力波在无限弹性介质中传播的特征线法解

双指数爆破荷载的应力时程关系是爆破工程的数值分析常用的动力输入。其在弹性岩石中传播衰减的特征线法解是工程爆破数值分析中弹性以及弹塑性分析的重要参照,也是动力学数值分析算法开发的重要高阶校准工具。应用统一弹性波的特征线法解求解源自圆形孔洞的双指数爆破荷载在理想介质中的传播。介绍了源自圆形孔洞的双指数爆破荷载在周围弹性介质传播过程中的径向应力、周向应力、位移以及速度波的特征线法解的具体过程,在MATLAB中进行了4种波的特征线法解的数值实施。对径向应力、周向应力、位移以及速度波的特征线法解的结果进行了多方面讨论。结果表明,特征线法是求解双曲型偏微分方程的一种有效方法,双指数爆破荷载在弹性岩石中传播衰减的特征线法解可为爆破工程数值分析以及动力学数值分析算法开发提供重要的支撑。     

Elastic waves along a cylindrical borehole in a poroelastic medium saturated by two immiscible fluids

The propagation of elastic waves along a cylindrical borehole filled with/without liquid and embedded in an infinite porous medium saturated by two immiscible fluids has been studied. The theory of porous media saturated by two immiscible fluids developed by Tuncay and Corapcioglu (1997) is employed. Frequency equations determining the phase velocity of axial symmetric waves are obtained. It is found that the surface waves along cylindrical borehole are dispersive. The dispersion equation of Rayleigh-type surface waves along the boundary of a poroelastic solid half-space saturated by two immiscible fluids is also obtained. Some special cases have been deduced and the dispersion curves are obtained numerically for a peculiar model. It is found that the density of fluids affects the Rayleigh mode.     

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.     

饱和多孔介质中的非均匀平面波

讨论了两相饱和多孔介质中的非均匀平面波,在所采用的多孔介质模型中,利用了每种组分的微观密度保持不变而与体积分数相关的宏观密度是变化的假设。针对理想无粘流体充填的弹性多孔固体,深入讨论了两种非均匀平面波的传播特性以及质点运动的轨迹线。     

Numerical analysis of 1‐D compression wave propagation in saturated poroelastic media

A numerical procedure is presented for the simulation of 1‐D compression wave propagation in saturated poroelastic media. The media are modelled as a two‐phase system consisting of compressible fluid and solids. Viscous coupling forces resulting from the relative motion between phases are characterized as Darcy type. The numerical procedure can account for effects of axial strain, nonlinear material behaviour, and various drained and undrained boundary conditions. Time integration is carried out explicitly and isothermal conditions are assumed. The method is capable of modelling shock wave fronts without introducing artificial viscosity. Numerical results are in close agreement with analytical solutions for several simplified cases and indicate that mass coupling may have important effects on fluid velocity and wave speed. Corresponding effects on solid velocity and wave speed are much smaller. Numerical results also indicate that damping occurs in a saturated poroelastic column and is dependent on the value of hydraulic conductivity. Copyright © 2007 John Wiley & Sons, Ltd.     

Consolidation of a poroelastic half-space with anisotropic permeability and compressible constituents by axisymmetric surface loading

The fully coupled Biot quasi-static theory of linear poroelasticity is used to study the consolidation of a poroelastic half-space caused by axisymmetric surface loads. The fluid and solid constituents of the poroelastic medium are compressible and its permeability in the vertical direction is different from its permeability in the horizontal direction. An analytical solution of the governing equations is obtained by taking the displacements and the pore pressure as the basic state variables and using a combination of the Laplace and Hankel transforms. The problem of an axisymmetric normal load is discussed in detail. An explicit analytical solution is obtained for normal disc loading. Detailed numerical computations reveal that the anisotropy in permeability as well as the compressibilities of the fluid and solid constituents of the poroelastic medium have significant effects on the consolidation of the half-space. The anisotropy in permeability may accelerate the consolidation process and may lead to a dilution in the theoretical prediction of the Mandel-Cryer effect. The compressibility of the solid constituents may also accelerate the consolidation process.