三维层状孔隙介质中弹性波的一种积分表达式I：理论

孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen（LAC）广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段.     

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

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

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.     

Seismic reflection and transmission coefficients at an air-water interface of saturated porous soil

Based on the modified Biot＇s theory of two-phase porous media, a study was presented on seismic reflection and transmission coefficients at an air-water interface of saturated porous soil media. The major differences between air-saturated soils and water-saturated soils were theoretically discussed, and the theoretical formulas of reflection and transmission coefficients at an air-water interface were derived. The characteristics of propagation and attenuation of elastic waves in air-saturated soils were given and the relations among the frequency, the angle of incidence and the reflection, transmission coefficients were analyzed by using numerical methods. Numerical results show that the propagation characteristic of the wave in air-saturated soils is great different from that in water-saturated soils. The frequency and the angle of incidence can have great influences on the reflection and transmission coefficients at interface. Some new cognition about the wave propagation is obtained and the study suggests that we may carefully pay attention to the influence of air on the dynamic analysis of seismic wave.     

双相各向异性介质中弹性波方程的有限元解法及波场模拟

基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波.     

饱和度对波在土层交界面的反射、透射系数的影响

把均匀流体概念引入Biot两相多孔介质动力理论中，用Biot的两相多孔介质模型模拟不完全饱和土层，给出SV波，P波从不完全饱和土层入射到弹性土层时，在土层交界面上反射，透射系数的表达式，结果表明与完全饱和相比，饱和度发生很小的变化就会对交界面上反射，透射系数产生很大的影响。今后应该重视饱和度变化对地震动力响应的影响。     

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.     

流体饱和多孔隙介质弹性波方程边界元解法研究

基于流体饱和多孔隙各向同性介质模型,本文首先推导了流体饱和多孔隙介质中弹性波传播的频率域系统动力方程及边界积分方程,然后给出了流体饱和多孔隙介质弹性波方程的基本解,最后,利用本文给出的边界元方法对流体饱和多孔隙各向同性介质中的弹性波传播进行了数值模拟.结果表明：不论是从固相位移,还是液相位移的地震合成记录都能看到明显的慢速P波,本文提出的流体饱和多孔隙介质弹性波边界元法是有效可行的.     

双相介质分界面上弹性波的反射与透射

本文基于Biot理论,推导出Zoeppritz形式的双相介质分界面上弹性波的反射与透射公式,对单相Zoeppritz公式与双相反射系数公式进行了比较,对双相介质含油、水、气不同流体时的反射规律以及孔隙度、渗透率、饱和度等储层参数对纵波反射的影响进行了研究.数值模拟分析表明,双相与单相反射公式的主要差异在于双相介质反射公式中考虑了液相、固液耦合相弹性模量的影响;油砂、水砂、气砂岩的慢P波反射差异明显;快P波反射对孔隙度的变化敏感,饱和度次之,对渗透率和频率的变化不敏感.     

油水两相渗流多孔介质弹性波传播动力学模型研究

弹性波在储层渗流场中的传播与衰减规律是研究波场强化采油动力学机理的重要基础.基于等效流体理论和饱和静态流体弹性波传播Biot理论,建立油水两非混相流体渗流条件下储层多孔介质中弹性波传播的动力学模型,通过算例求解与分析,发现含油水两相渗流储层多孔介质中同时存在着3种纵波P1、P2、P3和1种横波S;受频率和含水饱和度的影响,各波波速和品质因子呈现出不同变化规律,4种体波波速与频率、饱和度正相关,P1、P2波品质因子与饱和度正相关,P3和S波品质因子与饱和度负相关;最后,通过与传统静态弹性波模型结果对比,进一步分析了宏观渗流场对弹性波传播特征的影响规律,为揭示低频人工地震波辅助强化采油技术的动力学机理和工艺参数优化提供了重要理论依据.     

基于BISQ模型的三维双相各向异性介质数值模拟

Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.     

双相介质的AVO资料中蕴含的流体信息的定量分析

振幅随偏移距变化(AV0)弹性参数的反演是储层预测的一个重要方法.利用正演技术分析双相介质的地震资料中含有的流体信息,有助于提高AVO弹性参数反演的可靠性.本文应用Russell推出的计算纵波反射系数的AVO线性近似公式,对双相介质的纵波资料中所含有的流体信息进行定量分析,研究从AV0资料的反演中可以确定哪些弹性参数或它们的组合关系.文中考虑的弹性参数有流体因子、剪切模量和密度三者的相对扰动量.在时间城内对反射系数的正演计算公式进行线性化,得到灵敏度矩阵,并用奇异值分解(SVD)计算其特征值和特征向量采分析弹性参数的反演.由计算结果可知,在入射角较小时由AV0资料反演只能得到流体因子和密度两者的相对扰动量;若想同时获得三个参数的相对扰动量,就必须增大入射角.     

饱和多孔微极介质的波动方程及其势函数方程

土是由一定尺寸大小颗粒所构成的多孔介质,具有明显的颗粒特性,当土颗粒间的孔隙被流体（如水或油）充满时则成为饱和土.利用微极理论和Biot波动理论的研究成果,把饱和土中多孔固体骨架部分近似地视为微极介质,孔隙中的流体部分视为质点介质,获得饱和多孔微极介质的弹性波动方程.借鉴Greetsma理论,建立了饱和多孔微极介质弹性本构方程力学参数与相应单相介质弹性参数的相互关系,使饱和多孔微极介质弹性波动方程中的物理参数具有明确的物理意义,易于在试验中确定.运用场论理论把饱和多孔微极介质的波动方程简化为势函数方程,建立了饱和多孔微极介质中五种弹性波的弥散方程,数值分析了五种简谐体波在无限饱和多孔微极介质中的传播特性. 结果表明,P1波、P2波和剪切S1波的波速弥散曲线与经典饱和多孔介质基本相同,当频率小于临界频率ω0时旋转纵波θ波和横波S2波不存在,当频率大于临界频率ω0时,θ波和S2波的传播速度随频率增加而减小.     

基于临界孔隙度模型的地震波传播

基于岩石物里学中临界孔隙度模型,建立一种简洁的均匀弹性流体饱和孔隙介质模型,进行地震波传播研究.首先定义了构建目标模型的基本力学模型:介绍了全孔隙度区间内基本力学模型和目标孔隙介质的含义,其中基本力学模型除了完全弹性固体模型S和完全弹性流体模型F还包括临界孔隙模型C.然后通过等效力学模型推出了目标力学模型介质本构关系的组分表达形式.文中分别通过直接求取弹性参数的表达形式和运用应力应变关系两种方法得到介质模型的本构关系,进而得到该模型波动方程的组分表达形式.最后对这种介质模型进行了地震波传播的数值模拟,结合模拟结果分析孔隙对地震波传播的影响.     

淤积泥砂对水平地运动作用时刚性坝面动压力的影响

文章基于固体介质，流体介质，液固两相饱和介质的时域显示有限元波动分析方法，研究了可压缩库水条件下淤积泥砂对水平地运动作用时刚性坝面上动压力的影响，研究中分析了作为两相介质处理的淤积泥砂的饱和度、厚度、渗透系数和孔隙率对坝面动压力的影响，还比较了作为流程两相介质，单相固体介质、单相重流体介质处理的不同淤积泥砂层模型间的计算结果差异。     

Dynamic stiffness matrix of a poroelastic multi-layered site and its Green＇s functions

Few studies of wave propagation in layered saturated soils have been reported in the literature. In this paper, a general solution of the equation of wave motion in saturated soils, based on one kind of practical Biot‘s equation,was deduced by introducing wave potentials. Then exact dynamic-stiffness matrices for a poroelastic soil layer and halfspace were derived, which extended Wolf‘s theory for an elastic layered site to the case of poroelasticity, thus resolving a fundamental problem in the field of wave propagation and soil-structure interaction in a poroelastic layered soil site. By using the integral transform method, Green‘s functions of horizontal and vertical uniformly distributed loads in a poroelastic layered soil site were given. Finally, the theory was verified by numerical examples and dynamic responses by comparing three different soil sites. This study has the following advantages: all parameters in the dynamic-stiffness matrices have explicitly physical meanings and the thickness of the sub-layers does not affect the precision of the calculation which is very convenient for engineering applications. The present theory can degenerate into Wolf‘s theory and yields numerical results approaching those for an ideal elastic layered site when porosity tends to zero.     

Determining maximum shear stress in confined substrate from elastic wave reflection coefficient

The relationship between the maximum shear stress in a substrate solid and the elastic wave reflection coefficient from the interface between the substrate solid and an overlying solid half-space is investigated. Both substrate and overlying solid media are assumed to be initially isotropic and stress-free. Then as the substrate is subjected to horizontal confined stresses it becomes anisotropic. It is shown that longitudinal and shear wave reflection coefficients are related to the degree of stress induced anisotropy in the substrate medium. From this relation the confined stress level and the maximum shear stress generated on the vertical planes of the substrate are estimated. Authors in their previous investigation computed plane wave reflection coefficient in a biaxially compressed solid substrate immersed in a fluid. This paper reports for the first time how the maximum shear stress in a biaxially compressed substrate medium can be measured from the plane wave reflection coefficients when the overlying medium is also a solid half-space.     

Reflection and refraction of plane harmonic waves at an interface between elastic solid and porous solid saturated by viscous liquid

A general solution of Biot's field equations governing small motions of a porous solid saturated by viscous liquid is employed to study the reflection and refraction at the interface between an elastic solid and a liquid-saturated porous solid. The incident wave is assumed to be plane and homogeneous, propagating through the isotropic elastic solid. The poroelastic solid is considered to be a dissipative one. Amplitude and energy ratios are computed numerically for a particular model. With first-order corrections for the porosity of solid and viscosity of liquid, the limiting cases of low and high frequencies are computed.     

Green＇s functions for uniformly distributed loads acting on an inclined line in a poroelastic layered site

Based on one type of practical Biot＇s equation and the dynamic-stiffness matrices ofa poroelastic soil layer and half-space, Green＇s functions were derived for unitformly distributed loads acting on an inclined line in a poroelastie layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green＇s functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.     

Interaction of multiple courses of wave‐induced fluid flow in layered porous media

Different theoretical and laboratory studies on the propagation of elastic waves in layered hydrocarbon reservoir have shown characteristic velocity dispersion and attenuation of seismic waves. The wave‐induced fluid flow between mesoscopic‐scale heterogeneities (larger than the pore size but smaller than the predominant wavelengths) is the most important cause of attenuation for frequencies below 1 kHz. Most studies on mesoscopic wave‐induced fluid flow in the seismic frequency band are based on the representative elementary volume, which does not consider interaction of fluid flow due to the symmetrical structure of representative elementary volume. However, in strongly heterogeneous media with unsymmetrical structures, different courses of wave‐induced fluid flow may lead to the interaction of the fluid flux in the seismic band; this has not yet been explored. This paper analyses the interaction of different courses of wave‐induced fluid flow in layered porous media. We apply a one‐dimensional finite‐element numerical creep test based on Biot's theory of consolidation to obtain the fluid flux in the frequency domain. The characteristic frequency of the fluid flux and the strain rate tensor are introduced to characterise the interaction of different courses of fluid flux. We also compare the behaviours of characteristic frequencies and the strain rate tensor on two scales: the local scale and the global scale. It is shown that, at the local scale, the interaction between different courses of fluid flux is a dynamic process, and the weak fluid flux and corresponding characteristic frequencies contain detailed information about the interaction of the fluid flux. At the global scale, the averaged strain rate tensor can facilitate the identification of the interaction degree of the fluid flux for the porous medium with a random distribution of mesoscopic heterogeneities, and the characteristic frequency of the fluid flux is potentially related to that of the peak attenuation. The results are helpful for the prediction of the distribution of oil–gas patches based on the statistical properties of phase velocities and attenuation in layered porous media with random disorder.