基于横向各向同性BISQ方程的弹性波传播数值模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制. 近年来，利用同时处理这两种力学机制的BISQ(Biot-Squirt)模型，弹性波衰减和频散的问题已被广泛研究；然而基于BISQ方程的波场数值模拟尚未见到公开的报道．本文从BISQ方程出发，利用交错网格方法对横向各向同性孔隙介质中不同频率和相界情况，以及双层介质中的弹性波传播进行数值模拟，研究了在同时考虑两种流动机制作用情况下地震波和声波的传播特性及传播过程中出现的各种波动现象．      

基于Biot-Squirt方程的波场模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制,对地震波和声波的传播均产生重要影响. Dvorkin和Nur提出了同时包含Biot流动和喷射流动力学机制的统一的BISQ(Biot-Squirt)模型,基于这一模型,尽管有关弹性波在多孔隙介质中的衰减和频散问题已被广泛研究,然而,基于BISQ波传播方程的波场数值模拟至今仍未见报道. 本文从同时包含两种力学机制的孔隙弹性波方程出发,利用FCT有限差分法对含流体孔隙各向同性介质中的地震波和声波进行了数值模拟,并与基于Biot流动的Biot理论之模拟结果进行比较. 数值模拟结果表明:同时包含Biot流动和喷射流动影响的地震波和声波速度比仅包含Biot流动作用的地震波和声波速度慢,慢P波的衰减比根据Biot理论模拟的慢P波衰减更强.     

基于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.     

含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究

本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.     

基于非饱和多孔隙介质BISQ模型的储层参数反演

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制．文中基于同时包含这两种力学机制的非饱和多孔隙BISQ模型,利用小生境遗传算法实现了储层参数(孔隙度、渗透率、含流体饱和度等)的反演．结果表明,本文方法在储层参数的反演过程中目标函数收敛快、且具有较强的抗“噪声”干扰性能,当观测噪声ε≤5%时,储层参数的反演精度很高．最后,通过实测数据的反演应用验证了该方法处理储层参数反演问题的有效性和非饱和BISQ模型的准确性．     

多种物理机制耦合作用下的储层介质参数反演研究

孔隙介质的黏弹性、孔隙流体的Biot流动和喷射流动是影响波传播的重要物理机制.本文分别基于弹性和黏弹性BISQ模型,利用自适应杂交遗传算法研究了多种物理机制耦合作用条件下储层介质参数反演.为了测试自适应杂交遗传算法的有效性,本文分别利用自适应杂交遗传算法和传统实数编码遗传算法对含有不同噪声的理论合成数据进行了反演试算.对比理论合成数据反演结果可知,自适应杂交遗传算法具有抗干扰能力强且收敛速度快的特点,是一种有效的储层介质参数反演方法.同时本文也利用不同频率尺度和不同温度条件下的P波和S波实测数据进行了联合反演.对比研究表明,黏弹性BISQ模型能够很好地解释不同频率尺度的波频散特征,不仅能够很好地预测P波速度,而且也能够很好地预测S波速度,从而证明了黏弹性BISQ模型能够准确地描述低频条件下的波频散.     

Propagation and attenuation of Rayleigh waves in a semi-infinite unsaturated poroelastic medium

An analytical model for describing the propagation and attenuation of Rayleigh waves along the free surface of an elastic porous medium containing two immiscible, viscous, compressible fluids is developed in the present study based on the poroelastic equations formulated by Lo et al. [Lo WC, Sposito G, Majer E. Wave propagation through elastic porous media containing two immiscible fluids. Water Resour Res 2005;41:W02025]. The dispersion equation obtained is complex-valued due to viscous dissipation resulting from the relative motion of the solid to the pore fluids. As an excitation frequency is stipulated, the dispersion equation that is a cubic polynomial is numerically solved to determine the phase speed and attenuation coefficient of Rayleigh waves in Columbia fine sandy loam permeated by an air–water mixture. Our numerical results show that, corresponding to three dilatational waves, there is also the existence of three different modes of Rayleigh wave in an unsaturated porous medium, which are designated as the R1, R2, and R3 waves in descending order of phase speed, respectively. The phase speed of the R1 wave is non-dispersive (frequency-independent) in the frequency range we examined (10 Hz–10 kHz) and decreases as water saturation increases, whose magnitude ranges from 20% to 49% of that of the first dilatational wave with respect to water content. However, it is revealed numerically that the R2 and R3 waves are functions of excitation frequency. Given the same water saturation and excitation frequency, the phase speeds of the R2 and R3 waves are found to be approximately 90% of those of the second and third dilatational waves, respectively. The R1 wave has the lowest attenuation coefficient whereas the R3 wave attenuates highest.     

Biot流动和喷射流动耦合作用下波传播的FCT紧致差分模拟

本文从含流体多孔隙介质中同时包含Biot流动和喷射流动两种力学机制的BISQ(Biot-Squirt)方程出发,利用FCT(Flux-Corrected Transport)紧致差分方法数值模拟了在Biot流和喷射流共同作用下的波在含流体多孔隙各向同性介质中的传播.通过与仅受Biot流动作用下的波场结果对比,我们研究了...     

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

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

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.     

基于BISQ模型的双相介质位移场Green函数

以BISQ模型的波传播方程为基础,研究了弹性波在孔隙各向 同性介质中的传播速度及与流体的Biot流动和喷射流动力学机制的关系,进一步利用场势分 解和δ函数的性质,给出了同时受Biot流动和喷射流动两种动力学机制作用下,两相介 质波场位移在集中力作用下的Green函数,并讨论了喷射流动在Green函数中的表现形式和作用     

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.     

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

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

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

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

Propagation of seismic waves in patchy-saturated porous media: double-porosity representation

Propagation of harmonic plane waves is studied in a patchy-saturated porous medium. Patchy distribution of the two immiscible fluids is considered in a porous frame with uniform skeletal properties. A composition of two types of patches, connected through continuous paths, constitutes a double-porosity medium. Different compressibilities of pore-fluids in two porous phases facilitate the wave-induced fluid-flow in this composite material. Constitutive relations are considered with frequency-dependent complex elastic coefficients, which define the dissipative behaviour of porous aggregate due to the flow of viscous fluid in connected patches. Relevant equations of motion are solved to explain the propagation of three compressional waves and one shear wave in patchy-saturated porous solids. A numerical example is solved to illustrate dispersion in phase velocity and quality factor of attenuated waves in patchy-saturated porous materials. Role of fluid–solid inertial coupling in Darcy's law is emphasized to keep a check on the dispersion of wave velocities in the porous composite. Effects of patchy saturation on phase velocities and quality factors of attenuation are analysed using the double-porosity formulation as well as the reduced single-porosity equivalents.     

基于改进BISQ机制的双相HTI介质波传播伪谱法模拟与特征分析

改进BISQ(Biot-Squirt)机制在不引入特征喷流长度的情况下,将含流体孔隙介质中Biot流动和喷射流动两种重要的力学机制有机地结合起来,且各相关参数具有明确物理意义和可实现性.本文将改进BISQ机制一维孔隙流体压力公式推广到三维具有水平对称轴横向各向同性介质(HTI介质)情况,结合裂缝各向异性理论,给出了基于改进BISQ机制的双相HTI介质模型及其二维三分量波传播方程,采用伪谱法求解该方程,进行了不同相界、不同频率以及双层地质结构情况下该类介质中波场的数值模拟与特征分析.数值模拟结果表明:伪谱法模拟精度高,压制网格频散效果好,可以得到高精度的波场快照和合成记录;基于改进BISQ机制的双相HTI介质模型兼具裂缝各向异性特征和孔隙弹性特征,其为从双相各向异性理论角度深入研究裂缝性储层的地震响应奠定了理论基础.     

Pseudo Rayleigh wave in a partially saturated non-dissipative porous solid

Propagation of surface waves is studied at the pervious boundary of a porous solid saturated with a mixture of two immiscible fluids. An approach, based on continuum mixture theory, is used to derive a secular equation for the propagation of harmonic waves at the stress-free plane surface of this non-dissipative medium. Numerical analysis shows that this secular equation may not represent the propagation of true surface wave in the porous aggregate. Then, this equation is solved numerically for the propagation of pseudo Rayleigh wave or the leaky surface waves. To ensure the existence of pseudo Rayleigh wave, capillary effect between two (wetting and non-wetting) pore-fluids is related to the partial saturation. Effects of porosity and partial saturation coupled with capillary effect are observed on the phase velocity of pseudo Rayleigh waves in sandstone saturated with water-CO₂ mixture.     

近地表强衰减介质中的地震波传播模型

建立符合油气储层近地表复杂介质强衰减性质的数学-物理模型是油气勘探和开发的重要课题.本文针对D'Alembert黏弹性介质模型存在的量纲不统一和不能充分刻画孔隙介质结构的不足,通过引入孔隙度和渗透率等参数修改耗散项,改进了原D'Alembert模型,获得了能精细刻画具有强吸收衰减特征的近地表复杂介质模型,即改进的近地表黏弹性模型.基于这种改进的模型,推导了波频散和衰减的表达式,并研究了孔隙度、流体黏度等物理参数对波频散和衰减的影响,获得了相应的规律性认识.为了验证新模型预测近地表介质中波衰减的有效性,本文将新模型应用于胜利油田YX工区的近地表实测数据,同时与弹性Biot模型和BISQ模型以及黏弹性BISQ模型进行了比较.结果表明,与其它三个模型比较,改进的近地表黏弹性模型能够很好地刻画近地表介质的强衰减性,而且新模型所涉及到的物理参数明显少于其它几种模型,有利于新模型在油气储层近地表复杂介质地震勘探的实际应用.     

含混合裂隙、孔隙介质的纵波衰减规律研究

地下多孔介质中的孔隙类型复杂多样,既有硬孔又有扁平的软孔.针对复杂孔隙介质,假设多孔介质中同时含有球型硬孔和两种不同产状的裂隙(硬币型、尖灭型裂隙),当孔隙介质承载载荷时,考虑两种不同类型的裂隙对于孔隙流体压力的影响,建立起Biot理论框架下饱和流体情况含混合裂隙、孔隙介质的弹性波动方程,并进一步求取了饱和流体情况下仅由裂隙引起流体流动时的含混合裂隙、孔隙介质的体积模量和剪切模量,随后,在此基础上讨论了含混合裂隙、孔隙介质在封闭条件下地震波衰减和频散的高低频极限表达式.最后计算了给定模型的地震波衰减和频散,发现地震波衰减曲线呈现"多峰"现象,速度曲线为"多频段"频散.针对该模型分析讨论了渗透率参数、裂隙纵横比参数以及流体黏滞性参数对于地震波衰减和频散的影响,表明三个参数均为频率控制参数.     

Reflection and transmission of inhomogeneous waves in a composite porous solid saturated by two immiscible fluids

This paper is concerned with reflection and transmission of a plane, elastic, and inhomogeneous wave striking obliquely at some discontinuity inside a porous medium composed of two distinct solids and saturated by two immiscible fluids. It is found that four P‐ and two SV‐waves are reflected, whereas four P‐ and two SV‐waves are transmitted at the interface. All reflected and transmitted waves are inhomogeneous in nature and specified with different directions of propagation and attenuation vectors. An expression for the Umov–Poynting energy flux vector is derived for the system. Continuity of energy flux along normal to the interface gives 12 required boundary conditions. Expressions of amplitude and energy ratios of various reflected and transmitted waves are derived. Variations in amplitude and energy coefficients of reflected and transmitted waves with angle of incidence are numerically studied for a porous matrix composed of shaley sandstone and clay, saturated with water and oil. The effects of change in oil saturation and volume fraction of clay are also observed on amplitude ratios. Numerical simulation reveals that the change in sign in the difference of capillary pressure across the interface causes jump in the values of amplitude ratios of all waves.