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.     

The Peculiarities of Linear Wave Propagation in Double Porous Media

The features of propagation of longitudinal and transverse waves (LW and TW) in fractured porous medium (FPM) saturated with liquid are investigated by methods of multiphase mechanics. The mathematical model of FPM accounting for inequality of velocities and pressures of liquid in pores and fractures, liquid mass exchange and nonstationary interaction forces is developed. Processes of monochromatic wave propagation are studied. The dispersion relation is obtained and the effect of model parameters on wave propagation is analysed. It is established that one transverse and three longitudinal waves propagate in FPM saturated with liquid. The fastest LW is a deformational wave and the two others are filtrational. Filtrational waves attenuate much stronger than deformational and transverse waves. Distinction of velocities and pressures in liquid in various pore systems provides an explanation for the existence of the two filtrational waves in porous medium with two different characteristic sizes of pores.     

Wave propagation in fractured porous media

A theory of wave propagation in fractured porous media is presented based on the double-porosity concept. The macroscopic constitutive relations and mass and momentum balance equations are obtained by volume averaging the microscale balance and constitutive equations and assuming small deformations. In microscale, the grains are assumed to be linearly elastic and the fluids are Newtonian. Momentum transfer terms are expressed in terms of intrinsic and relative permeabilities assuming the validity of Darcy's law in fractured porous media. The macroscopic constitutive relations of elastic porous media saturated by one or two fluids and saturated fractured porous media can be obtained from the constitutive relations developed in the paper. In the simplest case, the final set of governing equations reduce to Biot's equations containing the same parameters as of Biot and Willis.Now at Izmir Institute of Technology, Anafartalar Cad. 904, Basmane 35230, Izmir, Turkey.     

Squirt-Flow in Fluid-Saturated Porous Media: Propagation of Rayleigh Waves

Propagation of attenuated waves is studied in a squirt-flow model of porous solid permeated by two different pore regimes saturated with same viscous fluid. Presence of soft compliant microcracks embedded in the grains of stiff porous rock defines the double-porosity formation. Microcracks and pores respond differently to the compressional effect of a propagating wave, which induces the squirt-flow from microcracks to pores. Elastodynamics of constituent particles in porous aggregate is represented through a single-porosity formulation, which involves the frequency-dependent complex moduli. This formulation is deduced as a special case of double-porosity formation allowing the wave-induced flow of pore-fluid. This squirt-flow model of porous solid supports the attenuated propagation of two compressional waves and one shear wave. Superposition of these body waves, subject to stress-free surface, defines the propagation of Rayleigh wave. This wave is governed by a complex irrational dispersion equation, which is solved numerically after rationalising into an algebraic equation. For existence of Rayleigh wave, a complex solution of the dispersion equation should represent a leaky wave, which decays for propagation along any direction in the semi-infinite medium. A numerical example is solved to analyse the effects of squirt-flow on phase velocity, attenuation and polarisation of the Rayleigh waves, for different combinations of parameters. Numerical results suggest the existence of an additional (second) Rayleigh wave in the squirt-flow model of dissipative porous solids.     

On dispersive propagation of surface waves in patchy saturated porous media

Frequency-dependent velocity and attenuation for Rayleigh-wave propagation along a vacuum/patchy saturated porous medium interface are investigated in the low frequency band (0.1–1000 Hz). Conventional patchy saturation models for compressional waves are extended to account for Rayleigh wave propagation along a free surface. The mesoscopic interaction of fluid and solid phases, as a dominant loss mechanism in patchy saturated media, significantly affects Rayleigh-wave propagation and attenuation. Researches on the dispersion characteristics at low frequencies with different gas fractions in patchy saturated media also demonstrate a strong correlation between the Rayleigh-wave mode and the fast compressional wave. Especially, the strongest attenuation with the maximum value of 1/Q$1/Q$ for Rayleigh waves are obtained in the frequency range of 1–200 Hz. Numerical results show that the significant dependence of velocity and attenuation on frequencies and gas fractions presents a distinctive dynamical response of Rayleigh waves in the time domain.     

Difference between the longitudinal Frenkel-Biot waves in water-and gas-saturated porous media

The problem of the propagation of longitudinal Biot waves in a porous medium saturated with a weakly compressible liquid (water) or a gas is considered theoretically. The frequency dependence of the phase velocities and damping coefficients is investigated numerically. It is shown that for a certain relationship between the parameters of the porous medium and the saturating fluid there is a “critical” frequency at which the properties of longitudinal waves of both kinds are identical. An analytical expression for this “critical” frequency is obtained. It is shown that for a gas-saturated porous medium, at a certain frequency, in both longitudinal waves the relative gas-matrix motion changes type. Assuming that the saturating-gas behavior corresponds to an adiabatic equation of state, an estimate is obtained for the threshold pore pressure necessary for the restructuring of the relative motion. The wave associated with matrix deformation is shown to have a high damping coefficient in a porous medium saturated with a weakly compressible liquid (water in the case considered) but to be only weakly damped in a gas-saturated porous medium.     

Propagation and Evolution of Wave Fronts in Two-Phase Porous Media

An investigation is made into the propagation and evolution of wave fronts in a porous medium which is intended to contain two phases: the porous solid, referred to as the skeleton, and the fluid within the interconnected pores formed by the skeleton. In particular, the microscopic density of each real material is assumed to be unchangeable, while the macroscopic density of each phase may change, associated with the volume fractions. A two-phase porous medium model is concisely introduced based on the work by de Boer. Propagation conditions and amplitude evolution of the discontinuity waves are presented by use of the idea of surfaces of discontinuity, where the wave front is treated as a surface of discontinuity. It is demonstrated that the saturation condition entails certain restrictions between the amplitudes of the longitudinal waves in the solid and fluid phases. Two propagation velocities are attained upon examining the existence of the discontinuity waves. It is found that a completely coupled longitudinal wave and a pure transverse wave are realizable in the two-phase porous medium. The discontinuity strength of the pore-pressure may be determined by the amplitude of the coupled longitudinal wave. In the case of homogeneous weak discontinuities, explicit evolution equations of the amplitudes for two types of discontinuity waves are derived.     

A mixture theory analysis for the surface-wave propagation in an unsaturated porous medium

The mixture theory is employed to the analysis of surface-wave propagation in a porous medium saturated by two compressible and viscous fluids (liquid and gas). A linear isothermal dynamic model is implemented which takes into account the interaction between the pore fluids and the solid phase of the porous material through viscous dissipation. In such unsaturated cases, the dispersion equations of Rayleigh and Love waves are derived respectively. Two situations for the Love waves are discussed in detail: (a) an elastic layer lying over an unsaturated porous half-space and (b) an unsaturated porous layer lying over an elastic half-space. The wave analysis indicates that, to the three compressional waves discovered in the unsaturated porous medium, there also correspond three Rayleigh wave modes (R1, R2, and R3 waves) propagating along its free surface. The numerical results demonstrate a significant dependence of wave velocities and attenuation coefficients of the Rayleigh and Love waves on the saturation degree, excitation frequency and intrinsic permeability. The cut-off frequency of the high order mode of Love waves is also found to be dependent on the saturation degree.     

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

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

Reflection and Transmission of Elastic Waves from the Interface of a Fluid-saturated Porous Solid and a Double Porosity Solid

The reflection and transmission characteristics of an incident plane P1 wave from the interface of a fluid-saturated single porous solid and a fluid-saturated double porosity solid are investigated. The fluid-saturated porous solid is modeled with the classic Biot’s theory and the double porosity medium is described by an extended Biot’s theory. In a double-porosity model with dual-permeability there exist three compressional waves and a shear wave. The effects of the incident angle and frequency on amplitude ratios of the reflected and transmitted waves to the incident wave are discussed. Two boundary conditions are discussed in detail: (a) Open-pore boundary and (b) Sealed-pore boundary. Numerical results reveal that the characteristics of the reflection and transmission coefficients to the incident angle and the frequency are quite different for the two cases of boundary conditions. Properties of the bulk waves existing in the fluid-saturated porous solid and the double porosity medium are also studied.     

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.     

Low-frequency dilatational wave propagation through unsaturated porous media containing two immiscible fluids

An analytical theory is presented for the low-frequency behavior of dilatational waves propagating through a homogeneous elastic porous medium containing two immiscible fluids. The theory is based on the Berryman–Thigpen–Chin (BTC) model, in which capillary pressure effects are neglected. We show that the BTC model equations in the frequency domain can be transformed, at sufficiently low frequencies, into a dissipative wave equation (telegraph equation) and a propagating wave equation in the time domain. These partial differential equations describe two independent modes of dilatational wave motion that are analogous to the Biot fast and slow compressional waves in a single-fluid system. The equations can be solved analytically under a variety of initial and boundary conditions. The stipulation of “low frequency” underlying the derivation of our equations in the time domain is shown to require that the excitation frequency of wave motions be much smaller than a critical frequency. This frequency is shown to be the inverse of an intrinsic time scale that depends on an effective kinematic shear viscosity of the interstitial fluids and the intrinsic permeability of the porous medium. Numerical calculations indicate that the critical frequency in both unconsolidated and consolidated materials containing water and a nonaqueous phase liquid ranges typically from kHz to MHz. Thus engineering problems involving the dynamic response of an unsaturated porous medium to low excitation frequencies (e.g., seismic wave stimulation) should be accurately modeled by our equations after suitable initial and boundary conditions are imposed.     

Exact solutions to one-dimensional transient response of incompressible fluid-saturated single-layer porous media

Based on the Biot theory of porous media,the exact solutions to onedimensional transient response of incompressible saturated single-layer porous media under four types of boundary conditions are developed.In the procedure,a relation between the solid displacement u and the relative displacement w is derived,and the well-posed initial conditions and boundary conditions are proposed.The derivation of the solution for one type of boundary condition is then illustrated in detail.The exact solutions for the other three types of boundary conditions are given directly.The propagation of the compressional wave is investigated through numerical examples.It is verified that only one type of compressional wave exists in the incompressible saturated porous media.     

Compressional waves in fluid-saturated porous solid containing a penny-shaped crack

Diffraction of normal compression waves by a penny-shaped crack in a fluid-saturated porous medium is investigated. Two wave types are considered, namely, compressional wave of the first kind, and the second kind. The former, also known as fast wave, propagates primarily through the solid, whereas the latter or slow wave, propagates mainly in the fluid. Each wave propagates in the medium along with induced wave of the same type in the companion constituent of the material. Application of Biot’s theory in conjunction with integral transform technique reduces the problem to a mixed boundary-value problem whose solution is in turn governed by a Fredholm integral equation of the second kind. Near-field and far-field solutions are obtained in terms of the dynamic stress-intensity factor and the scattering cross section, respectively. They are of particular importance to the linear elastic fracture mechanics (LEFM) and in the scattering theory of elastic waves. The mode I stress-intensity factors are computed numerically for a set of selected material property values, and shown graphically for various mass density and viscosity-to-permeability ratios. The obtained results reveal significant impact of the presence of pore fluid upon the stress-intensity factors, both magnitudes and frequencies at their peak values. The influence of the fluid is also observed from the calculated scattering cross sections of the scattered far-field. Accuracy of the present solution procedure is verified by comparing the numerical results with existing results in the limiting case of dry elastic materials.     

多孔饱和半空间上刚体垂直振动的轴对称混合边值问题

研究圆柱形刚体在多孔饱和半空间上的垂直振动．首先应用Ｈａｎｋｅｌ变换求解多孔饱和固体的动力基本方程———Ｂｉｏｔ波动方程．然后按混合边值条件建立多孔饱和半空间上刚体垂直振动的对偶积分方程，用Ａｂｅｌ变换化对偶积分方程为第二类Ｆｒｅｄｈｏｌｍ积分方程．文末给出了多孔饱和半空间表面动力柔度系数的计算曲线．     

Rayleigh waves in anisotropic porous media and the polarization vector method

In this paper, the propagation of Rayleigh waves in orthotropic non-viscous fluid-saturated porous half-spaces with sealed surface-pores and with impervious surface is investigated. The main aim of the investigation is to derive explicit secular equations and based on them to examine the effect of the material parameters and the boundary conditions on the propagation of Rayleigh waves. By employing the method of polarization vector the explicit secular equations have been derived. These equations recover the ones corresponding to Rayleigh waves propagating in purely elastic half-spaces. It is shown from numerical examples that the Rayleigh wave velocity depends strongly on the porosity, the elastic constants, the anisotropy, the boundary conditions and it differs considerably from the one corresponding to purely elastic half-spaces. Remarkably, in the fluid saturated porous half-spaces, Rayleigh waves may travel with a larger velocity than that of the shear wave, a fact that is impossible for the purely elastic half-spaces.     

准饱和土体中圆形衬砌对弹性波的散射

采用Vardoulakis和Beskos提出的准饱和土体的波动控制方程,根据Helmholtz矢量分解定理,得到了准饱和土中P1波(快压缩波)、P2波(慢压缩波)和S波(剪切波)的波数的势函数表达式.将准饱和土体和圆形衬砌视为各向同性的均质体,运用波函数展开法将入射波、散射波和折射波的势函数展开成Fourier-Bessel函数的级数形式,根据准饱和土体与衬砌边界处应力和位移连续及衬砌内完全自由的边界条件,得到了平面P1波入射下,准饱和土体内深埋圆形衬砌的散射系数和折射系数的理论解,通过数值计算分析了饱和度对准饱和土体和衬砌的DSCF(动应力集中因子)及准饱和土体的PPCF(孔压集中因子)的影响规律,结果表明:准饱和土体的DSCF随着饱和度的增大而减小,衬砌的DSCF基本不受饱和度的影响,而准饱和土体的PPCF则随着饱和度的增大而增大.     

Reflection and transmission of elastic waves at an elastic/porous solid saturated by two immiscible fluids

Wave propagation in a porous elastic medium saturated by two immiscible fluids is investigated. It is shown that there exist three dilatational waves and one transverse wave propagating with different velocities. It is found that the velocities of all the three longitudinal waves are influenced by the capillary pressure, while the velocity of transverse wave does not at all. The problem of reflection and refraction phenomena due to longitudinal and transverse wave incident obliquely at a plane interface between uniform elastic solid half-space and porous elastic half-space saturated by two immiscible fluids has been analyzed. The amplitude ratios of various reflected and refracted waves are found to be continuous functions of the angle of incidence. Expression of energy ratios of various reflected and refracted waves are derived in closed form. The amplitude ratios and energy ratios have been computed numerically for a particular model and the results obtained are depicted graphically. It is verified that during transmission there is no dissipation of energy at the interface. Some particular cases have also been reduced from the present formulation.     

川藏公路地质环境与整治改建方案的思考

川藏公路由于地质环境复杂、建设标准低、后遗病害多,抗灾能力差,泥石流、滑坡、山崩、雪害、水毁等自然灾害频繁发生,公路阻车断道严重。国家投入巨资进行整治改建,并取得了明显的效果,但由于自然环境特殊、影响因素复杂,许多特大型、大型工程地质病害问题还没有可行、可靠的解决方案。本文通过分析川藏公路沿线的地质环境和灾害特点,总结历年整治改建和经验的教训,提出川藏公路建设的途径、可能达到的目标和应采用的原则。     

多孔材料中声波的传播与演化

采用两相多孔介质的拉格朗日模型来描述一种理论流体充填的多孔弹性固体材料,其中孔隙度的变化满足一个附加的平衡方程。