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

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

Acoustic wave propagation in saturated porous media: reformulation of the Biot/Squirt flow theory

A model of wave propagation in fluid-saturated porous media is developed where the principal fluid/solid interaction mode affecting the propagation of the acoustic wave results from the conjunction of the Biot and the Squirt flow mechanism. The difference between the original Biot/Squirt (BISQ) flow theory and the new theory, which we call the reformulated BISQ, is that the average fluid pressure term appearing in the dynamic equation for a two component solid/fluid continuum is independent of squirt flow length. P-velocity and attenuation relate to measurable rock physical parameters: the Biot's poroelastic constants, porosity, permeability, pore fluid compressibility and viscosity. Modelling shows that velocity and attenuation dispersion obtained using the reformulated BISQ theory are of the same order of magnitude as those obtained using the original BISQ theory. Investigation on permeability effect on velocity and attenuation dispersion indicate that the transition zone in velocity and attenuation peak, occurring both at the relaxation frequency, shifts toward high frequency when permeability decreases. This behaviour agrees with Biot's theory prediction.     

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.     

Full frequency-range transient solution for compressional waves in a fluid-saturated viscoacoustic porous medium1

An analytical transient solution is obtained for propagation of compressional waves in a homogeneous porous dissipative medium. The solution, based on a generalization of Biot's poroelastic equations, holds for the low- and high-frequency ranges, and includes viscoelastic phenomena of a very general nature, besides the Biot relaxation mechanism. The viscodynamic operator is used to model the dynamic behaviour associated with the relative motion of the fluid in the pores at all frequency ranges. Viscoelasticity is introduced through the standard linear solid which allows the modelling of a general relaxation spectrum. The solution is used to study the influence of the material properties, such as bulk moduli, porosity, viscosity, permeability and intrinsic attenuation, on the kinematic and dynamic characteristics of the two compressional waves supported by the medium. We also obtain snapshots of the static mode arising from the diffusive behaviour of the slow wave at low frequencies.     

Analysis of attenuation and dispersion of propagating wave due to the coexistence of three fluid phases in the pore volume

In exploration geophysics, the efforts to extract subsurface information from wave characteristics exceedingly depend on the construction of suitable rock physics model. Analysis of different rock physics models reveals that the strength and magnitude of attenuation and dispersion of propagating wave exceedingly depend on wave-induced fluid flow at multiple scales. In current work, a comprehensive analysis of wave attenuation and velocity dispersion is carried out at broad frequency range. Our methodology is based on Biot's poroelastic relations, by which variations in wave characteristics associated with wave-induced fluid flow due to the coexistence of three fluid phases in the pore volume is estimated. In contrast to the results of previous research, our results indicate the occurrence of two-time pore pressure relaxation phenomenon at the interface between fluids of disparate nature, that is, different bulk modulus, viscosity and density. Also, the obtained results are compatible with numerical results for the same 1D model which are accounted using Biot's poroelastic and quasi-static equation in frequency domain. Moreover, the effects of change in saturation of three-phase fluids were also computed which is the key task for geophysicist. The outcomes of our research reveal that pore pressure relaxation phenomenon significantly depends on the saturation of distinct fluids and the order of saturating fluids. It is also concluded that the change in the saturation of three-phase fluid significantly influences the characteristics of the seismic wave. The analysis of obtained results indicates that our proposed approach is a useful tool for quantification, identification and discrimination of different fluid phases. Moreover, our proposed approach improves the accuracy to predict dispersive behaviour of propagating wave at sub-seismic and seismic frequencies.     

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.     

流体分布对松散介质中P波速度和衰减的影响

为了深入研究流体对岩石中弹性波速度和衰减的影响,必须考虑到流体的分布和粘性。引入气体包裹体模型来研究粘性流体的分布对松散介质中P波速度和衰减的影响,用气泡平均半径来描述流体分布的不均匀性,计算了不同气泡半径和频率下P波速度和衰减随饱和度变化的曲线,并与有效流体模型作了比较,由于流体喷流的存在会使Gassmann方程在高频下不适用,用干燥和饱和流体的P波、S波速度修正了理论曲线。测量了玻璃微珠中不同水饱和度下高频P波的速度和衰减,并尝试用峰值频率来计算衰减。此方法求出的Q和频谱比法求出的Q在干燥或饱和水时基本相同,随饱和度的变化规律也基本一致,但衰减峰的大小有差异。根据实测值得似事经修正的波速和衰减理论曲线从而估算出气泡平均半径,认为P波速度和衰减不仅与饱和度有关而且也与介质内部气体－液体压力平衡有关。     

Pore fluid viscosity effects on P‐ and S‐wave anisotropy in synthetic silica‐cemented sandstone with aligned fractures

Ultrasonic (500 kHz) P‐ and S‐wave velocity and attenuation anisotropy were measured in the laboratory on synthetic, octagonal‐shaped, silica‐cemented sandstone samples with aligned penny‐shaped voids as a function of pore fluid viscosity. One control (blank) sample was manufactured without fractures, another sample with a known fracture density (measured from X‐ray CT images). Velocity and attenuation were measured in four directions relative to the bedding fabric (introduced during packing of successive layers of sand grains during sample construction) and the coincident penny‐shaped voids (fractures). Both samples were measured when saturated with air, water (viscosity 1 cP) and glycerin (100 cP) to reveal poro‐visco‐elastic effects on velocity and attenuation, and their anisotropy. The blank sample was used to estimate the background anisotropy of the host rock in the fractured sample; the bedding fabric was found to show transverse isotropy with shear wave splitting (SWS) of 1.45 ± 1.18% (i.e. for S‐wave propagation along the bedding planes). In the fractured rock, maximum velocity and minimum attenuation of P‐waves was seen at 90° to the fracture normal. After correction for the background anisotropy, the fractured sample velocity anisotropy was expressed in terms of Thomsen's weak anisotropy parameters ε, γ & δ. A theory of frequency‐dependent seismic anisotropy in porous, fractured, media was able to predict the observed effect of viscosity and bulk modulus on ε and δ in water‐ and glycerin‐saturated samples, and the higher ε and δ values in air‐saturated samples. Theoretical predictions of fluid independent γ are also in agreement with the laboratory observations. We also observed the predicted polarisation cross‐over in shear‐wave splitting for wave propagation at 45° to the fracture normal as fluid viscosity and bulk modulus increases.     

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

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

基于非饱和Biot-Stoll模型的海底沉积物介质声频散特性研究

利用含气非饱和Biot-Stoll模型研究了声波在海底表层沉积物介质中的传播,讨论了骨架耗散、含气饱和度对快纵波、慢纵波和横波速度和衰减的影响,并与Biot模型的结果进行了对比.研究结果表明:孔隙流体黏滞耗散与骨架耗散共同影响声波传播速度和衰减,低频情况下骨架耗散引起的衰减占主要地位,高频情况下骨架耗散引起的衰减较小;少量气体(<1%)的引入显著改变了快纵波速度,气体含量的变化对快纵波衰减影响很大,低频情况下气体对慢纵波速度的影响不大,而对横波速度的影响较大,气体含量的变化对慢纵波和横波衰减影响较小.利用超声波测量系统测量了一例杭州湾海底沉积物样品的纵波速度和衰减,当含气量趋近0%时,Biot-Stoll模型预测的纵波速度和实验测量结果较为一致.     

An equivalent medium model for wave simulation in fractured porous rocks

Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane.     

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

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

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

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

Comparison of P‐wave attenuation models of wave‐induced flow

Wave‐induced oscillatory fluid flow in the vicinity of inclusions embedded in porous rocks is one of the main causes for P‐wave dispersion and attenuation at seismic frequencies. Hence, the P‐wave velocity depends on wave frequency, porosity, saturation, and other rock parameters. Several analytical models quantify this wave‐induced flow attenuation and result in characteristic velocity–saturation relations. Here, we compare some of these models by analyzing their low‐ and high‐frequency asymptotic behaviours and by applying them to measured velocity–saturation relations. Specifically, the Biot–Rayleigh model considering spherical inclusions embedded in an isotropic rock matrix is compared with White's and Johnson's models of patchy saturation. The modeling of laboratory data for tight sandstone and limestone indicates that, by selecting appropriate inclusion size, the Biot‐Rayleigh predictions are close to the measured values, particularly for intermediate and high water saturations.     

Properties of low‐frequency trapped mode in viscous‐fluid waveguides

We derived the velocity and attenuation of a generalized Stoneley wave being a symmetric trapped mode of a layer filled with a Newtonian fluid and embedded into either a poroelastic or a purely elastic rock. The dispersion relation corresponding to a linearized Navier–Stokes equation in a fracture coupling to either Biot or elasticity equations in the rock via proper boundary conditions was rigorously derived. A cubic equation for wavenumber was found that provides a rather precise analytical approximation of the full dispersion relation, in the frequency range of 10^?3 Hz to 10³ Hz and for layer width of less than 10 cm and fluid viscosity below 0.1 Pa· s [100 cP]. We compared our results to earlier results addressing viscous fluid in either porous rocks with a rigid matrix or in a purely elastic rock, and our formulae are found to better match the numerical solution, especially regarding attenuation. The computed attenuation was used to demonstrate detectability of fracture tip reflections at wellbore, for a range of fracture lengths and apertures, pulse frequencies, and fluid viscosity.     

砾岩储层地震波传播方程:三重孔隙结构模型

针对砾岩储层的砂、砾、泥三重孔隙结构特征,本文分析砾岩孔隙区域、砂岩孔隙区域以及泥岩孔隙区域相互之间的孔隙流体流动机制,将静态的砾岩骨架本构方程与动态的孔隙流体运动方程联立,提出了复杂砾岩储层的弹性波传播理论方程.采用实测砾岩储层参数,在算例中与双重孔隙介质理论进行对比分析,验证了本文理论方程的合理性;基于三重孔隙介质模型,分析不同储层环境下纵波的传播特征,结果显示:随流体黏滞系数增大,在衰减-频率轴坐标系中,砾与砂、砂与泥孔隙区域间局域流导致的两个衰减峰向低频端移动,而Biot全局流导致的衰减峰向高频端移动;嵌入体尺寸及背景相介质渗透率的变化,主要影响纵波速度频散曲线沿频率轴左、右平移,不影响波速低频、高频极限幅值;嵌入体含量及孔隙度的变化改变了岩石干骨架的弹性、密度参数,不仅影响速度频散曲线沿频率轴平移,而且影响其上、下限幅值;砾包砂包泥三重孔隙介质模型所预测的衰减曲线中,低频段"第一个衰减峰"主要由砾岩孔隙区域与砂岩孔隙区域之间的局域流导致,中间频段"第二个衰减峰"主要由砂岩孔隙区域与泥岩孔隙区域之间的局域流导致,超声频段"第三个衰减峰"由Biot全局流导致.对慢纵波传播特征的分析显示,砂岩骨架(局部孔隙度较大)内部的宏观孔隙流体流动造成的耗散明显强于砾岩与泥岩骨架.     

Seismic wave attenuation in fluid-saturated porous media

Contrary to the traditional view, seismic attenuation in Biot's theory of fluid-saturated porous media is due to viscous damping of local (not global) pore-fluid motion. Since substantial inhomogeneities in fluid permeability of porous geological materials are to be expected, the regions of highest local permeability contribute most to the wave energy dissipation while those of lowest permeability dominate the fluid flow rate if they are uniformly distributed. This dichotomy can explain some of the observed discrepancies between computed and measured attenuation of compressional and shear waves in porous earth. One unfortunate consequence of this result is the fact that measured seismic wave attenuation in fluid-filled geological materials cannot be used directly as a diagnostic of the global fluid-flow permeability.     

横向各向同性多孔介质中的地震波传播

基于各向异性多孔介质中的广义Ｂｉｏｔ理论，导出了横向各向同性多孔介质中波传播的特征方程．指出在多孔介质中有４种类型的频散和耗散波传播：准纵波ＱＰ１（快纵波）、准纵波ＱＰ２（慢纵波）、准横波ＱＳＶ和横渡ＳＨ．文中给出了４种波速度的解析表达式．数值计算频率曲线和衰减曲线与Ｓｃｈｍｉｔｔ（１９８９）用均值处理得到的结果类似．还给出了波传播过程中３种类型准体波之间的耦合系数（或称转换系数）．     

Simulation of full-waveform log in saturated cracked formations using Hudson's approach

We study the propagation of elastic waves that are generated in a fluid‐filled borehole surrounded by a cracked transversely isotropic medium. In the model studied the anisotropy and borehole axes coincide. To obtain the effective elastic moduli of a cracked medium we have applied Hudson's theory that enables the determination of the overall properties as a function of the crack orientation in relation to the symmetry axis of the anisotropic medium. This theory takes into account the hydrodynamic mechanism of the elastic‐wave attenuation caused by fluid filtration from the cracks into a porous matrix. We have simulated the full waveforms generated by an impulse source of finite length placed on the borehole axis. The kinematic and dynamic parameters of the compressional, shear and Stoneley waves as functions of the matrix permeability, crack orientation and porosity were studied. The modelling results demonstrated the influence of the crack‐system parameters (orientation and porosity) on the velocities and amplitudes of all wave types. The horizontally orientated cracks result in maximal decrease of the elastic‐wave parameters (velocities and amplitudes). Based on the fact that the shear‐ and Stoneley‐wave velocities in a transversely isotropic medium are determined by different shear moduli, we demonstrate the feasibility of the acoustic log to identify formations with close to horizontal crack orientations.