Reflection and transmission of attenuated waves at the boundary between two dissimilar poroelastic solids saturated with two immiscible viscous fluids

The phenomenon of reflection and transmission of plane harmonic waves at the plane interface between two dissimilar poroelastic solids saturated with two immiscible viscous fluids is investigated. Both porous media are considered dissipative due to the presence of viscosity in pore‐fluids. Four attenuated (three dilatational and one shear) waves propagate in such a dissipative porous medium. A finite non‐dimensional parameter is used to define the effective connections between the surface‐pores of two media at their common interface. Another finite parameter represents the gas‐share in the saturation of pores. An attenuated wave in a dissipative medium is described through the specification of directions of propagation and maximum attenuation. A general representation of an attenuated wave is defined through its inhomogeneous propagation, i.e., different directions for propagation and attenuation. Incidence of an inhomogeneous wave is considered at the interface between two dissipative porous solids. This results in four reflected and four transmitted inhomogeneous waves. Expressions are derived for the partition of incident energy among the reflected and transmitted waves. Numerical examples are studied to determine the effects of saturating pore fluid, frequency, surface‐pore connections and wave inhomogeneity on the strengths of reflected and transmitted waves. Interaction energy due to the interference of different (inhomogeneous) waves is calculated in both the dissipative porous media to verify the conservation of incident energy.     

Love waves in an inhomogeneous fluid saturated porous layered half-space with linearly varying properties

The paper is concerned with the propagation of the Love waves in an inhomogeneous transversely isotropic fluid saturated porous layered half-space with linearly varying properties. The analysis is based on Biot's theory. Firstly, the dispersion equation in the complex form for the Love waves in an inhomogeneous porous layer is derived. Then the equation is solved by an iterative method. Detailed numerical calculation is presented for an inhomogeneous fluid saturated porous layer overlying a purely elastic half-space. The dispersion and attenuation of the Love waves are discussed. In addition, the upper and lower bounds of the Love wave speed are explored.     

平面S波在饱和多孔热弹性介质边界的反射问题研究

针对饱和多孔介质中热弹性波的传播特性问题,基于多孔介质理论和广义的热弹性模型,研究平面S波在饱和多孔热弹性介质边界上的反射问题.以考虑流-固耦合的饱和多孔介质波动方程和热-弹耦合的广义热弹性基本方程出发,建立饱和多孔介质的热-流-固耦合弹性波动模型.通过引入势函数并考虑自由透水和绝热的边界条件,经过理论推导最终给出在饱...     

Numerical study of seismic scattering and waveguide excitation in faulted coal seams

Finite‐difference P‐SV simulations of seismic scattering characteristics of faulted coal‐seam models have been undertaken for near‐surface P‐ and S‐wave sources in an attempt to understand the efficiency of body‐wave to channel‐wave mode conversion and how it depends on the elastic parameters of the structure. The synthetic seismograms clearly show the groups of channel waves generated at the fault: one by the downgoing P‐wave and the other by the downgoing S‐wave. These modes travel horizontally in the seam at velocities less than the S‐wavespeed of the rock. A strong Airy phase is generated for the fundamental mode. The velocity contrast between the coal and the host rock is a more important parameter than the density contrast in controlling the amplitude of the channel waves. The optimal coupling from body‐wave energy to channel‐wave energy occurs at a velocity contrast of 1.5. Strong guided waves are produced by the incident S‐sources for source angles of 75° to 90° (close to the near‐side face of the fault). As the fault throw increases, the amplitude of the channel wave also increases. The presence of a lower‐velocity clay layer within the coal‐seam sequence affects the waveguiding characteristics. The displacement amplitude distribution is shifted more towards the lower‐wavespeed layer. The presence of a ‘washout’ zone or a brecciated zone surrounding the fault also results in greater forward scattering and channel‐wave capture by the coal seam.     

两种不同势函数下半空间饱和多孔介质中Rayleigh波求解比较

分别对"考虑两种压缩波和幅值比例系数"和"考虑一种压缩波(P1或P2波)但不考虑幅值比例系数"两种不同势函数下的半空间饱和多孔介质中Rayleigh波求解进行详细推导,理论分析表明"考虑两种压缩波和幅值比例系数"下Rayleigh波求解推导更为严密,与饱和多孔介质中存在两种压缩波的事实相一致。在研究半空间饱和多孔介质中Rayleigh波时应采用"考虑两种压缩波和幅值比例系数"的势函数。     

The effect of inertial coupling on seismic reflection amplitudes

A problem of reflection and transmission of elastic waves at a plane interface between a uniform elastic solid half-space and a porous elastic half-space containing two immiscible fluids is investigated. The theory developed by Lo, Sposito and Majer for porous media containing two immiscible fluids is employed to find out the reflection and transmission coefficients. The incident wave is assumed to propagate through the uniform elastic half-space and two cases are considered. In the first case, a beam of plane longitudinal wave is assumed to be incident and in the second case, a beam of transverse wave is assumed to be incident at the interface. By taking granite as impervious elastic medium and columbia fine sandy loam containing air-water mixture as porous medium, reflection and transmission coefficients are obtained. By neglecting the inertial coupling coefficients, these coefficients are reduced to those obtained by Tomar and Arora using the theory of Tuncay and Corapcioglu. It is found that the inertial coupling parameters significantly affect the phase speeds and the amplitude ratios of the transmitted waves.     

Velocity–porosity relationships: Predictive velocity model for cemented sands composed of multiple mineral phases

Computer simulations are used to calculate the elastic properties of model cemented sandstones composed of two or more mineral phases. Two idealized models are considered – a grain‐overlap clay/quartz mix and a pore‐lining clay/quartz mix. Unlike experimental data, the numerical data exhibit little noise yet cover a wide range of quartz/cement ratios and porosities. The results of the computations are in good agreement with experimental data for clay‐bearing consolidated sandstones. The effective modulus of solid mineral mixtures is found to be relatively insensitive to microstructural detail. It is shown that the Hashin–Shtrikman average is a good estimate for the modulus of the solid mineral mixtures. The distribution of the cement phase is found to have little effect on the computed modulus–porosity relationships. Numerical data for dry and saturated states confirm that Gassmann's equations remain valid for porous materials composed of multiple solid constituents. As noted previously, the Krief relationship successfully describes the porosity dependence of the dry shear modulus, and a recent empirical relationship provides a good estimate for the dry‐rock Poisson's ratio. From the numerical computations, a new empirical model, which requires only a knowledge of system mineralogy, is proposed for the modulus–porosity relationship of isotropic dry or fluid‐saturated porous materials composed of multiple solid constituents. Comparisons with experimental data for clean and shaly sandstones and computations for more complex, three‐mineral (quartz/dolomite/clay) systems show good agreement with the proposed model over a very wide range of porosities.     

含两项不混合/单项黏性流体孔隙介质分界面上反透射系数特征研究

弹性孔隙介质分界面上的反透射系数特征,在岩性划分、流体识别、储层边界判识等方面有重要的应用.本文研究上层为含两项不混合黏性流体孔隙介质、下层为含单项黏性流体孔隙介质分界面上的反透射理论.首先根据两种孔隙介质分界面上的能量守恒得到边界条件,再将波函数、位移、应力与应变关系代入边界条件,推导出完全连通孔隙情况下,第一类纵波入射到孔隙介质分界面上的反透射系数方程.通过建立砂岩孔隙介质模型,分别分析不同孔隙流体类型、不同含油饱和度及不同入射角情况下,各类波的反透射系数特征.研究表明,第二、三类纵波反透射系数数值比第一类纵波小多个数量级,且两者对入射角的变化不敏感,但对孔隙流体性质、含油饱和度的变化较敏感,而横波反透射系数特征恰好与此相反;第一类纵波反透射系数特征比较复杂,入射角、孔隙流体的性质及含油饱和度的变化都对其产生影响.不同孔隙流体弹性物性的差异、孔隙介质中含油饱和度的变化及不同入射角引起垂向和切向应力分量的变化都会影响各类波的反透射系数特征,分析这些特征可以为研究储层含油气性提供理论基础.     

转换波和非转换波与射线参数p的关系

本文对Aki固体-固体分界面的反射和透射公式进行适当变换,使之成为一组适合不同坐标系下,固体、液体、气体和真空之间所有可能形成的分界面上,入射P、SV和SH波所形成的反射和透射公式,从而把不同边界条件的反射和透射方程统一起来,并且阐明了转换波和非转换波动力学特征的差异.指出:在不同弹性介质分界面上入射平面弹性波,转换波的反射系数和透射系数是射线参数p的奇函数;非转换波的反射系数和透射系数是射线参数p的偶函数.根据奇函数和偶函数的性质可见:垂直入射时不存在转换波.同时也为简化Zoeppricz方程,开展AVO分析打开新思路.     

Analytical expressions for Fresnel volumes and interface fresnel zones of seismic body waves. Part 2: Transmitted and converted waves. Head waves

Summary Fresnel volumes and interface Fresnel zones of transmitted and head waves are studied. The relation derived for transmitted waves may also be used for converted reflected waves. Considerable attention is devoted to the penetration of Fresnel volumes across structural interfaces, particularly for head waves.     

Reflection of elastic waves at an interfaces between two porous half-spaces filled with different fluids

The problem about the refraction of elastic waves at the interface of two half-spaces (rocks) filled with immiscible fluids is solved with the use of modified boundary conditions of the dynamics of saturated porous media [Nagy and Nayfeh, 1995]. The solution is obtained within the framework of the Frenkel’-Biot theory with allowance for the surface tension at the interface of fluids for a porous medium formed by two half-spaces, which differ only in the properties of the fluids filling them. Practically important cases of the reflection from the fluid-gas (water-air) and fluid-fluid (oil-water) interfaces are considered in detail. The calculations are performed for both harmonic waves and pulses. The possibility to determine in principle the structural factor characterizing the pore-space geometry of rocks from measured dynamic parameters of reflected waves is shown.     

Oblique wave scattering by undulating porous bottom in a two-layer fluid: Fourier transform approach

The problem involving scattering of oblique waves by small undulation on the porous ocean bed in a two-layer fluid is investigated within the framework of linearised theory of water waves where the upper layer is free to the atmosphere. In such a two-layer fluid, there exist waves with two different wave numbers (modes): wave with lower wave number propagates along the free surface whilst that with higher wave number propagates along the interface. When an oblique incident wave of a particular mode encounters the undulating bottom, it gets reflected and transmitted into waves of both modes so that some of the wave energy transferred from one mode to another mode. Perturbation analysis in conjunction with Fourier transform technique is used to derive the first-order corrections of velocity potentials, reflection and transmission coefficients at both modes due to oblique incident waves of both modes. One special type of undulating bottom topography is considered as an example to evaluate the related coefficients in detail. These coefficients are shown in graphical forms to demonstrate the transformation of water wave energy between the two modes. Comparisons between the present results with those in the literature are made for particular cases and the agreements are found to be satisfactory. In addition, energy identity, an important relation in the study of water wave theory, is derived with the help of the Green’s integral theorem.     

Reflection and refraction at an imperfectly bonded interface between poroelastic solid and cracked elastic solid

Porous solid is in contact with a cracked elastic solid at a plane interface between them. For the presence of vertically aligned microcracks, the elastic solid behaves transversely isotropic to wave propagation. The coefficients of elastic anisotropy depend on the crack density and crack porosity in the medium. A loose bonding is considered between the two solids so that a limiting case could be the welded contact. At the plane interface, the imperfection in welded bonding is represented by tangential slipping and, hence, results in the dissipation of a part of strain energy. Three types of waves propagate in an isotropic fluid-saturated porous medium, which are considered for incidence at the interface. Incidence of a wave results in three reflected waves and two refracted waves. Partition of incident energy among the reflected and refracted waves is studied for each incidence, varying from normal to grazing directions. Numerical example calculates the energy shares of reflected and refracted waves at the plane interface between water-saturated sandstone and basalt. These energy shares are computed and analyzed for different values of crack parameters as well as loose bonding parameter.     

THE SPECTRAL FUNCTION OF A VERTICALLY INHOMOGENEOUS MEDIUM1

Elastic, acoustic and electromagnetic waves in media consisting of vertically inhomogeneous layers are considered in a common formulation. The spectral function of a vertically inhomogeneous medium is the downward energy flux due to an impulsive source at the top of the first layer. A propagation-invariant form is used to derive several identities for the reflection and transmission matrices. When the top layer is bounded by a free surface, one of the expressions reduces to a formula derived by Kunetz for the one-dimensional wave equation. A source radiating upwards and downwards gives a discontinuity in the propagation-invariant form which is equal to the source energy. A new formula is derived for when the source is located just beneath the top interface of the layers.     

PP amplitude bias caused by interface scattering: are diffracted waves guilty?

This paper is concerned with the problem of interpretation of anomalous seismic amplitudes, induced by the amplitude‐scattering phenomenon. This phenomenon occurs in the vicinity of a crack distribution at the interface between elastic layers. The purpose of this work is to obtain a better understanding of the physics of this distinctive phenomenon, in order to interpret correctly the amplitudes of the reflected events. By analogy with studies in optics and in acoustics, we suggest that diffraction is widely involved in the amplitude‐scattering phenomenon. Analytical evaluation of the amount of energy carried by the reflected and the diffracted waves shows that neglecting diffraction in numerical models leads to local underestimation of the amplitude of waves reflected at interfaces with gas‐filled crack distribution.     

Poro‐acoustoelasticity of fluid‐saturated rocks

We generalize the classical theory of acoustoelasticity to the porous case (one fluid and a solid frame) and finite deformations. A unified treatment of non‐linear acoustoelasticity of finite strains in fluid‐saturated porous rocks is developed on the basis of Biot’s theory. A strain‐energy function, formed with eleven terms, combined with Biot’s kinetic and dissipation energies, yields Lagrange’s equations and consequently the wave equation of the medium. The velocities and dissipation factors of the P‐ and S‐waves are obtained as a function of the 2nd‐ and 3rd‐order elastic constants for hydrostatic and uniaxial loading. The theory yields the limit to the classical theory if the fluid is replaced with a solid with the same properties of the frame. We consider sandstone and obtain results for open‐pore jacketed and closed‐pore jacketed ‘gedanken’ experiments. Finally, we compare the theoretical results with experimental data.     

饱和土中深埋圆柱形衬砌洞室对瞬态平面波的散射

基于Biot饱和多孔介质动力学理论,运用Laplace变换和波函数展开法,根据饱和土体与衬砌结构交界面的连续条件和衬砌结构内边界上的应力自由条件,得到饱和土中深埋圆柱形衬砌洞室对瞬态平面P波和SV波散射问题的解答,该解答可以退化成为饱和土中深埋圆柱形空穴或弹性夹塞物的情形,并很容易转换成为对稳态波散射的解.通过与已有的相关问题的解析解答进行对比,验证了该解答的正确性.同时利用Laplace逆变换的数值方法,给出了饱和土和衬砌中应力和位移场在时域内的数值解,通过算例,分析了衬砌厚度、刚度对衬砌内边界处应力集中因子的影响.     

含黏滞流体各向异性孔隙介质中弹性波的频散和衰减

基于Biot理论,考虑液相的黏弹性变形和固液相接触面上的相对扭转,提出了含黏滞流体VTI孔隙介质模型.从理论上推导出,在该模型中除存在快P波、慢P波、SV波、SH波以外,还将存在两种新横波－慢SV波和慢SH波.数值模拟分析了6种弹性波的相速度、衰减、液固相振幅比随孔隙度、频率的变化规律以及快P波、快SV波的衰减随流体性质、渗透率、入射角的变化规律.结果表明慢SV波和慢SH波主要在液相中传播,高频高孔隙度时,速度较高;大角度入射时,快P波衰减表现出明显的各向异性,而快SV波的衰减则基本不变;储层纵向和横向渗透率存在差异时,快SV波衰减大的方向渗透率高.     

Characterization of seismic hazard and structural response by energy flux

Seismic safety of structures depends on the structure's ability to absorb the seismic energy that is transmitted from ground to structure. One parameter that can be used to characterize seismic energy is the energy flux. Energy flux is defined as the amount of energy transmitted per unit time through a cross-section of a medium, and is equal to kinetic energy multiplied by the propagation velocity of seismic waves. The peak or the integral of energy flux can be used to characterize ground motions. By definition, energy flux automatically accounts for site amplification. Energy flux in a structure can be studied by formulating the problem as a wave propagation problem. For buildings founded on layered soil media and subjected to vertically incident plane shear waves, energy flux equations are derived by modeling the building as an extension of the layered soil medium, and considering each story as another layer. The propagation of energy flux in the layers is described in terms of the upgoing and downgoing energy flux in each layer, and the energy reflection and transmission coefficients at each interface. The formulation results in a pair of simple finite-difference equations for each layer, which can be solved recursively starting from the bedrock. The upgoing and downgoing energy flux in the layers allows calculation of the energy demand and energy dissipation in each layer. The methodology is applicable to linear, as well as nonlinear structures.     

SV波全反射作用下地震动的合成及振动特性

越来越多的研究表明来自基岩的地震波并不是垂直地面向上传播的。地震波在斜入射与垂直入射时所产生的场地效应有很大不同,由于存在全反射现象,SV波在斜入射时产生的场地效应更为复杂。文章基于均匀弹性半空间地震波传播理论,分别推导得到SV波入射角在小于、大于等于临界角时地震动的计算表达式,通过模型算例研究SV波全反射作用下的地震动特性。研究发现：由SV波产生的地震动主要由入射波和反射波构成,滑行波作用可以忽略;地面运动轨迹具有面波旋转振动特点;随着入射角的增大,地面震动从以水平方向振动为主逐渐过渡到以垂直方向振动为主,两个方向的振动均小于入射波峰值的2倍。