波在双孔隙与单孔隙介质界面上的反射与透射

基于Biot理论和双重孔隙介质理论研究了弹性波在双重孔隙介质与流体饱和单一孔隙介质 界面的反射和透射问题，在界面上假定裂缝孔隙流体相对于固体骨架的位移为零，推导了反 射系数和透射系数的计算公式，数值讨论了反射系数和透射系数随入射角和频率的变化关 系. 同时，讨论了双重孔隙介质中3种压缩波(P-1, P-2和P-3波)和一种剪切波(S波) 的频散和衰减特性.     

Reflection and Transmission of Elastic Waves at the Interface Between Water and a Double Porosity Solid

The reflection and transmission of an incident plane wave at an interface between water and a fluid-saturated double porosity solid are investigated. The properties of the three compressional waves and one shear wave in double porosity media are discussed in detail. The effect of the viscosity, permeability, and porosity on the phase velocity and attenuation of four bulk waves is presented. Comparison between the double porosity theory and the Biot theory reveals that there exists obvious difference in the phase velocity, attenuation and the reflection and transmission coefficients. Two cases of boundary conditions, i.e., the sealed-pore and the open-pore boundaries, are introduced in the numerical calculation. At last, the effect of the fracture permeability and porosity on the reflection and transmission coefficients considered. It is revealed that the amplitudes of the reflected and transmitted waves strongly depend the boundary condition, fracture permeability, and porosity.     

Generalized ray expansion for pulse propagation and attenuation in fluid-saturated porous media

A theory suitable for studying pulses propagating through a layered fluid-saturated porous medium is presented. Biot's theory is used to describe the constitutive equation of a fluid-saturated porous solid. Since fast and slow compressional waves exist in a Biot solid even at normal incidence, there is mode conversion at the interface and, therefore, the transmission and reflection coefficients are 2x2 matrices. We use matrix methods in developing the solution of the wave propagation problem. A generalized ray expansion algorithm is obtained by using the Gauss-Seidel matrix iterative method. The arrivals of the fast and the slow waves are easily identified. A compact computational algorithm is developed using combinatorial analysis and the Cayley-Hamilton theorem.     

Low-Frequency Asymptotic Analysis of Seismic Reflection from a Fluid-Saturated Medium

Reflection of a seismic wave from a plane interface between two elastic media does not depend on the frequency. If one of the media is poroelastic and fluid-saturated, then the reflection becomes frequency-dependent. This paper presents a low-frequency asymptotic formula for the reflection of seismic plane p-wave from a fluid-saturated porous medium. The obtained asymptotic scaling of the frequency-dependent component of the reflection coefficient shows that it is asymptotically proportional to the square root of the product of the reservoir fluid mobility and the frequency of the signal. The dependence of this scaling on the dynamic Darcy’s law relaxation time is investigated as well. Derivation of the main equations of the theory of poroelasticity from the dynamic filtration theory reveals that this relaxation time is proportional to Biot’s tortuosity parameter.     

Elastic Waves in Swelling Porous Media

Time harmonic waves in a swelling porous elastic medium of infinite extent and consisting of solid, liquid and gas phases have been studied. Employing Eringen’s theory of swelling porous media, it has been shown that there exist three dilatational and two shear waves propagating with distinct velocities. The velocities of these waves are found to be frequency dependent and complex valued, showing that the waves are attenuating in nature. Here, the appearance of an additional shear wave is new and arises due to swelling phenomena of the medium, which disappears in the absence of swelling. The reflection phenomenon of an incident dilatational wave from a stress-free plane boundary of a porous elastic half-space has been investigated for two types of boundary surfaces: (i) surface having open pores and (ii) surface having sealed pores. Using appropriate boundary conditions for these boundary surfaces, the equations giving the reflection coefficients corresponding to various reflected waves are presented. Numerical computations are performed for a specific model consisting of sandstone, water and carbon dioxide as solid, liquid and gas phases, respectively, of the porous medium. The variations of phase speeds and their corresponding attenuation coefficients are depicted against frequency parameter for all the existing waves. The variations of reflection coefficients and corresponding energy ratios against the angle of incidence are also computed and depicted graphically. It has been shown that in a limiting case, Eringen’s theory of swelling porous media reduces to Tuncay and Corapcioglu theory of porous media containing two immiscible fluids. The various numerical results under these two theories have been compared graphically.     

Seismic responses of a hemispherical alluvial valley to SV waves： a three-dimensional analytical approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier–Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot’s dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.The project was supported by the National Natural Science Foundation of China (50478062 and 10532070) and Open Fund at the Key Laboratory of Urban Security and Disaster Engineering (Beijing University of Technology), Chinese Ministry of Education. The English text was polished by Keren Wang.     

Transverse wave at a plane interface between isotropic elastic and unsaturated porous elastic solid half-spaces

Based on the poroelasticity theory, this article investigates the reflection and transmission characteristics of an incident plane transverse wave at a plane interface between an isotropic elastic half-space and an unsaturated poroelastic solid half-space. For this purpose, the effect of the saturation degree and frequency on the properties of the four bulk waves in unsaturated porous medium, i.e., three longitudinal waves and one transverse wave, are discussed at first. Two general cases of mode conversion are considered: (i) The initial transverse wave is incident from an unsaturated poroelastic half-space to the interface, and (ii) the initial transverse wave is incident from an elastic solid half-space to the interface. The expressions for the partition of energy at the interface during transmission and reflection process of waves are presented in explicit forms. At last, numerical computations are performed for these two cases and the results obtained are depicted, respectively. The variation of the amplitude ratios and energy ratios with the saturation degree and incident angle is illustrated in detail. It is also verified that, at the interface, the sum of energy ratios is approximately equal to unity as expected.     

On solving the problem of reflection of linear waves in a fluid from a porous half-space saturated by this fluid

Solutions of the problem of reflection of a stepwise pressure wave in a linearly compressed fluid from a flat boundary of a porous medium of infinite length saturated by the same fluid are obtained in the acoustic approximation. Based on analytical solutions, a numerical analysis is performed to reveal the specific features of the reflected and incident waves, depending on porosity and permeability of the porous half-space. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 16–26, September–October, 2006.     

Seismic Reflection from an Interface Between an Elastic Solid and a Fractured Porous Medium with Partial Saturation

The theory of Tuncay and Corapcioglu (Transp Porous Media 23:237–258, 1996a) has been employed to investigate the possibility of plane wave propagation in a fractured porous medium containing two immiscible fluids. Solid phase of the porous medium is assumed to be linearly elastic, isotropic and the fractures are assumed to be distributed isotropically throughout the medium. It has been shown that there can exist four compressional waves and one rotational wave. The phase speeds of these waves are found to be affected by the presence of fractures, in general. Of the four compressional waves, one arises due to the presence of fractures in the medium and the remaining three are those encountered by Tuncay and Corapcioglu (J Appl Mech 64:313–319, 1997). Reflection and transmission phenomena at a plane interface between a uniform elastic half-space and a fractured porous half-space containing two immiscible fluids, are analyzed due to incidence of plane longitudinal/transverse wave from uniform elastic half-space. Variation of modulus of amplitude and energy ratios with the angle of incidence are computed numerically by taking the elastic half-space as granite and the fractured porous half-space as sandstone material containing non-viscous wetting and non-wetting fluid phases. The results obtained in case of porous half-space with fractures, are compared graphically with those in case of porous half-space without fractures. It is found that the presence of fractures in the porous half-space do affect the reflection/transmission of waves, which is responsible for raising the reflection and lowering the transmission coefficients.     

Dynamic Analysis of an Infinite Cylindrical Hole in a Saturated Poroelastic Medium

In this paper, the dynamic response of an infinite cylindrical hole embedded in a porous medium and subjected to an axisymmetric ring load is investigated. Two scalar potentials and two vector potentials are introduced to decouple the governing equations of Biot’s theory. By taking a Fourier transform with respect to time and the axial coordinate, we derive general solutions for the potentials, displacements, stresses and pore pressures in the frequency-wave-number domain. Using the general solutions and a set of boundary conditions applied at the hole surface, the frequency-wave-number domain solutions for the proposed problem are determined. Numerical inversion of the Fourier transform with respect to the axial wave number yields the frequency domain solutions, while a double inverse Fourier transform with respect to frequency as well as the axial wave number generates the time-space domain solution. The numerical results of this paper indicate that the dynamic response of a porous medium surrounding an infinite hole is dependant upon many factors including the parameters of the porous media, the location of receivers, the boundary conditions along the hole surface as well as the load characteristics.     

An Asymptotic Model of Seismic Reflection from a Permeable Layer

Analysis of compression wave propagation in a poroelastic medium predicts a peak of reflection from a high-permeability layer in the low-frequency end of the spectrum. An explicit formula expresses the resonant frequency through the elastic moduli of the solid skeleton, the permeability of the reservoir rock, the fluid viscosity and compressibility, and the reservoir thickness. This result is obtained through a low-frequency asymptotic analysis of Biot’s model of poroelasticity. A review of the derivation of the main equations from the Hooke’s law, momentum and mass balance equations, and Darcy’s law suggests an alternative new physical interpretation of some coefficients of the classical poroelasticity. The velocity of wave propagation, the attenuation factor, and the wave number are expressed in the form of power series with respect to a small dimensionless parameter. The absolute value of this parameter is equal to the product of the kinematic reservoir fluid mobility and the wave frequency. Retaining only the leading terms of the series leads to explicit and relatively simple expressions for the reflection and transmission coefficients for a planar wave crossing an interface between two permeable media, as well as wave reflection from a thin highly permeable layer (a lens). Practical applications of the obtained asymptotic formulae are seismic modeling, inversion, and attribute analysis.     

Wave propagation in liquid-saturated porous solid with micropolar elastic skelton at boundary surface

The present study is concerned with the reflection and transmission of plane waves at an interface between homogenous invisicid liquid half space and a micropolar liquid-saturated porous solid half space. The reflection and transmission coefficients of various reflected and transmitted waves with the angle of incident have been obtained. Numerical calculation has been performed for amplitude ratios of various reflected and transmitted waves. Micropolarity and porosity effects on the reflection and transmission coefficients have been depicted graphically. Some particular cases have been deduced from the present formulation.     

流体饱和标准线性粘弹性多孔介质中的平面波

研究了流体饱和不可压标准线性粘弹性多孔介质中平面波的传播和反射问题．在固相骨架小变形的假定下,得到了粘弹性多孔介质中波动方程的一般解,讨论了弥散关系和波的衰减特性．结果表明：在流体饱和不可压粘弹性多孔介质中,仅存在一个耦合纵波和一个耦合横波,纵波和横波的波速、衰减率等取决于孔隙流体与固相骨架间的相互作用以及固相骨架本身的粘性．同时,研究了半空间自由边界上入射波(纵波、横波)的反射问题。得到了非均匀反射波的波速、反射系数、衰减率等的表达式及其相关的数值结果．     

Waves in Nonlocal Elastic Solid with Voids

In this paper, the governing relations and equations are derived for nonlocal elastic solid with voids. The propagation of time harmonic plane waves is investigated in an infinite nonlocal elastic solid material with voids. It has been found that three basic waves consisting of two sets of coupled longitudinal waves and one independent transverse wave may travel with distinct speeds. The sets of coupled waves are found to be dispersive, attenuating and influenced by the presence of voids and nonlocality parameters in the medium. The transverse wave is dispersive but non-attenuating, influenced by the nonlocality and independent of void parameters. Furthermore, the transverse wave is found to face critical frequency, while the coupled waves may face critical frequencies conditionally. Beyond each critical frequency, the respective wave is no more a propagating wave. Reflection phenomenon of an incident coupled longitudinal waves from stress-free boundary surface of a nonlocal elastic solid half-space with voids has also been studied. Using appropriate boundary conditions, the formulae for various reflection coefficients and their respective energy ratios are presented. For a particular model, the effects of non-locality and dissipation parameter (\(\tau \)) have been depicted on phase speeds and attenuation coefficients of propagating waves. The effect of nonlocality on reflection coefficients has also been observed and shown graphically.     

Reflection and transmission of plane waves between two different fluid-saturated porous half-spaces

The present study deals with reflection and transmission of plane waves between two different fluid-saturated porous half-spaces, where longitudinal and transversal waves impinge obliquely onto the interface. Amplitude ratios of various reflected and transmitted waves are obtained. Variations of amplitude ratios with the angle of incidence are depicted graphically. A particular case of reflection at the free surface of the fluid-saturated porous half-spaces is discussed.     

Scattering of water waves by vertical porous barriers: An analytical approach

An analytical approach is proposed here to study scattering of deep water waves by a submerged or a surface piercing vertical porous barrier. It involves a connection between two wave potentials of which one is the solution of a boundary value problem associated with wave scattering by the porous barrier and the other is the solution of a complementary type problem where barrier and gap positions are interchanged and solid barrier takes the position of the porous barrier. The connection also involves an auxiliary or a connection wave potential. The potential for the solid barrier problem involves incident wave forcing while the auxiliary potential describes a solid barrier type problem that involves a non-physical forcing. The solution procedure of Ursell (Ursell, 1947) is chosen to solve these boundary value problems explicitly in the case of normal wave incidence as it also determines necessary exact behavior of the potential at the barrier edge. The reflection coefficients are also connected and the reflection amplitudes of the normally incident wave against the vertical porous barriers are obtained analytically. Numerical results for reflection and transmission coefficients are presented.     

Reflection and refraction of attenuated waves at boundary of elastic solid and porous solid saturated with two immiscible viscous fluids

The propagation of elastic waves is studied in a porous solid saturated with two immiscible viscous fluids.The propagation of three longitudinal waves is represented through three scalar potential functions.The lone transverse wave is presented by a vector potential function.The displacements of particles in different phases of the aggregate are defined in terms of these potential functions.It is shown that there exist three longitudinal waves and one transverse wave.The phenomena of reflection and refraction due to longitudinal and transverse waves at a plane interface between an elastic solid half-space and a porous solid half-space saturated with two immiscible viscous fluids are investigated.For the presence of viscosity in pore-fluids,the waves refracted to the porous medium attenuate in the direction normal to the interface.The ratios of the amplitudes of the reflected and refracted waves to that of the incident wave are calculated as a nonsingular system of linear algebraic equations.These amplitude ratios are used to further calculate the shares of different scattered waves in the energy of the incident wave.The modulus of the amplitude and the energy ratios with the angle of incidence are computed for a particular numerical model.The conservation of the energy across the interface is verified.The effects of variations in non-wet saturation of pores and frequencies on the energy partition are depicted graphically and discussed.     

Numerical analysis of the isolation of the vibration due to Rayleigh waves by using pile rows in the poroelastic medium

The isolation of the vibration due to harmonic Rayleigh waves using pile rows embedded in a saturated poroelastic half-space is investigated in this study. Based on Biot’s theory and the potential function method, the free field solution for Rayleigh waves along the surface of the poroelastic half-space is derived first. The fundamental solution for a harmonic circular patch load applied in the poroelastic half-space are obtained in terms of Biot’s theory and the integral transform method. Using Muki’s method and the fundamental solution for the circular patch load as well as the Rayleigh waves solution for the poroelastic half-space, the second kind of Fredholm integral equations in the frequency domain for pile rows are derived. Numerical solution of the integral equations yields the dynamic response of the pile–soil system to incident Rayleigh waves. Influences of various parameters on the vibration isolation effect of piles rows are investigated numerically. Numerical results suggest that for the same vibration source, the same pile rows will produce a better vibration isolation effect for the poroelastic medium than for a single phase elastic medium. Also, stiffer piles tend to have better vibration isolation effect than flexible piles. Moreover, the pile length and the spacing between neighboring piles in each pile row have significant influence on the vibration isolation effect of pile rows.     

Reflection and transmission of plane waves at the interface of an elastic half-space and a fractional order thermoelastic half-space

The problem of reflection and transmission due to longitudinal and transverse waves incident obliquely at a plane interface between uniform elastic solid half-space and fractional order thermoelastic solid half-space has been studied. It is found that the amplitude ratios of various reflected and refracted waves are functions of angle of incidence and frequency of incident wave and are influenced by the fractional order thermoelastic properties of media. The expressions of amplitude ratios and energy ratios have been computed numerically for a particular model. The variation of amplitude and energy ratios with angle of incidence is shown graphically. The conservation of energy at the interface is verified.     

Reflection and Refraction of Anti-Plane Shear Waves from a Moving Phase Boundary

The reflection and refraction of anti-plane shear waves from an interface separating half-spaces with different moduli is well understood in the linear theory of elasticity. Namely, an oblique incident wave gives rise to a reflected wave that departs at the same angle and to a refracted wave that, after transmission through the interface, departs at a possibly different angle. Here we study similar issues for a material that admits mobile elastic phase boundaries in anti-plane shear. We consider an energy minimal equilibrium state in anti-plane shear involving a planar phase boundary that is perturbed due to an incident wave of small magnitude. The phase boundary is allowed to move under this perturbation. As in the linear theory, the perturbation gives rise to a reflected and a refracted wave. The orientation of these waves is independent of the phase boundary motion and determined as in the linear theory. However, the phase boundary motion affects the amplitudes of the departing waves. Perturbation analysis gives these amplitudes for general small phase boundary motion, and also permits the specification of the phase boundary motion on the basis of additional criteria such as a kinetic relation. A standard kinetic relation is studied to quantify the subsequent energy partitioning and dissipation on the basis of the properties of the incident wave.