应力诱导饱和流体岩石弹性波各向异性

各向同性岩石中裂隙随机分布,施加单轴方向的压应力后,与压应力方向垂直的缝隙会优先闭合.应力作用下裂隙的闭合导致岩石呈现出横向各向异性的特征,进而影响弹性波的传播特征.前人对干燥岩石中应力引起的弹性波速度各向异性开展了大量研究,但未考虑充填流体对于应力诱导各向异性和弹性波传播的影响.弹性波通过应力作用的岩石时会诱导微孔隙结构上的"挤喷流"效应,从而造成弹性波显著的频散和衰减.在本文中,我们分析应力诱导各向异性岩石中挤喷流对弹性波传播的影响.利用硬币型裂隙模型,同时考虑挤喷流和应力诱导各向异性的影响,我们建立了应力作用下含流体岩石各向异性模型.基于本文的模型,我们研究了裂隙介质中单轴应力和充填流体对弹性波传播的影响.我们的模型成功预测了弹性波速度随应力的增加,同时表明了挤喷流造成纵波的各向异性频散更显著.     

周期性层状含孔隙、裂隙介质模型纵波衰减特征

地震波在含孔隙、裂隙斑块饱和介质传播过程中会诱发多个尺度孔隙流体流动而产生衰减和速度频散.在含有宏观尺度“Biot流”和介观尺度“局域流”衰减诱导机制的周期性层状孔隙介质模型基础上,引入了微观尺度硬币型和尖灭型裂隙“喷射流”的影响,构建了周期性层状含孔隙、裂隙介质模型.利用双解耦弹性波动方程的方法数值计算了该模型地震频带的纵波衰减和速度频散并与周期性层状孔隙介质模型做了对比研究.分析了该模型在不同裂隙参数（裂隙密度、裂隙纵横比）及裂隙体积含量下的纵波衰减和频散特征,裂隙密度越高对于纵波衰减和频散的影响越大,裂隙纵横比越小,由裂隙引起的纵波衰减部分向高频段移动,裂隙体积含量越少,纵波衰减先降低后小幅增加再降低,频散速度增加,并逐渐接近于周期性层状孔隙介质模型的纵波衰减和频散速度曲线.最后研究了周期性层状含孔隙、裂隙介质模型有效平面波模量的高低频极限以及流固相对位移在该模型中的分布特征.     

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

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

基于三维孔隙网络模型的纵波频散衰减特征分析

波场在含流体的孔隙介质中传播时会产生频散和衰减现象.波场的频散和衰减与孔隙介质的岩石物理属性有关,包括孔隙度、渗透率、流体属性等.现有的三维裂缝/软孔隙网络模型利用椭圆截面纵横比的变化模拟从扁裂缝、软孔隙到硬孔隙的多种情况,而未考虑同时包含孔隙和裂缝的全局性网络空间.为了更好地描述裂缝-孔隙空间,本文提出同时包含裂缝和孔隙的三维裂缝-孔隙网络模型,并给出渗透率的计算方法.通过体积平均法推导了三维裂缝/软孔隙网络模型和三维裂缝-孔隙网络模型的波动方程,利用平面波分析方法得到纵波频散/衰减曲线的表达式,同时应用数值模拟研究了总孔隙度、裂缝孔隙度、裂缝纵横比、裂缝数密度、孔隙流体黏度对纵波衰减和速度频散特征的影响.结果表明,在三维裂缝-孔隙网络模型下,总孔隙度、裂缝参数等对纵波频散衰减特征的影响与三维裂缝/软孔隙网络模型相似.具体表现为:纵波在高频段内出现频散和衰减现象.孔隙度的变化主要影响逆品质因子曲线峰值的大小;裂缝数密度主要控制速度显著变化的范围;裂缝纵横比对纵波速度和特征频率有显著影响.     

数字岩心宽频带动态应力应变模拟方法及其对含裂隙致密岩石频散和衰减特征的表征

利用数字岩石物理技术表征复杂非均质多孔岩石跨频段的频散和衰减特征对于综合利用多尺度的地球物理数据进行地层非均质性的刻画具有重要的意义.现有的描述波致流体流动引起的频散和衰减效应的数字岩心动态应力应变模拟方法主要为单频率模拟方法,需要在不同的频率进行多次模拟才能刻画频散和衰减特征.本文提出了宽频带动态应力应变模拟方法,通过给数字岩心加载一个快速趋于恒定的宏观应变,采用波场正演技术求解数字岩心内部流固耦合的应变场和应力场,模拟数字岩心内部的应力松弛过程,从而通过一次模拟计算目标频段范围内连续的频散和衰减曲线.该方法可以成功地用于刻画含裂隙致密岩石挤喷流效应引起的速度频散和衰减特征,并通过数值模拟较为系统地揭示致密地层中控制挤喷流效应的主控物理因素,这些认识也与现有挤喷流效应的理论模型有较好的吻合.     

流体饱和周期性层状孔隙介质模型平行层面传播的纵波频散和衰减研究

孔隙介质中的地震波传播一直是油气地震勘探领域的研究热点和难点问题.该科学难题源自不同尺度的裂隙、孔隙、溶洞与岩石骨架之间的耦合作用,导致地震波场特征复杂.目前相关的研究主要集中于探索孔隙介质中地震波的传播机制及地震响应的特征与变化规律,包括对地震波在复杂孔隙介质中传播,进行比较精确的数学物理描述以及数值实现.地球物理学家们集中于研究垂直于地层层面方向入射的地震波频散和衰减,而忽略了实际地球介质中的地震波是以任意角度(方向)入射并进行传播的普遍性情况.在前人的研究基础上,本文的创新之处在于将纵波的入射方向扩展到平行于流体饱和的周期性层状孔隙介质模型层面方向.针对流体饱和的周期性层状孔隙介质模型,提出了介观波致流(Wave-induced Fluid Flow, WIFF)对流体饱和孔隙层状介质中平行于层面方向入射的纵波频散、衰减及频变各向异性的新模型.利用准静态Biot孔弹性方程推导出了模型的孔隙压力、流体流动速度、平均应力和平均应变等物理量的解析表达式,进而得到流体饱和的周期性层状孔隙介质复纵波模量的精确解析解．然后,利用复纵波模量讨论了纵波速度频散、衰减和频变各向异性特征,讨论了背景...     

基于微观孔隙结构特征的速度频散和衰减模拟

孔隙尺度的喷射流流动是引起地震波速度频散和衰减的重要机制之一.目前,大多数喷射流模型仅考虑硬孔隙与微裂隙之间的局部流动,而忽略了具有不同孔隙纵横比微裂隙间的喷射流作用.为了研究各种类型孔隙间的流体流动效应,本文对经典喷射流模型进行了扩展,通过结合等效介质理论和孔隙结构模型,根据从干燥岩石超声速度-压力曲线中提取的微裂隙孔隙纵横比分布,求取出岩石中各种微裂隙的体积压缩系数,并在此基础上,利用孔隙空间的压力松弛效应对微裂隙间的喷射流效应进行了模拟,并运用Biot理论描述了硬孔隙间的宏观流动效应.扩展后的理论模型不仅考虑了微裂隙与硬孔隙间的局部流动、硬孔隙与硬孔隙间的Biot宏观流,还加入了微裂隙与微裂隙间的喷射流作用,且模型的高、低频极限始终与Mavko-Jizba理论和Gassmann方程保持一致.模型应用分析发现,对于砂岩和大部分致密灰岩样品,扩展模型均能给出与超声实验测量数据吻合良好的估计结果.此外,扩展模型预测的速度频散及衰减表明,喷射流机制在地震和测井频段发挥着重要作用,其速度频散曲线由低频至高频呈逐渐增大趋势,不具有明显的快速变化特征,与经典频散曲线形态存在显著差异;在低有效压力下,频散和衰减程度较大,喷射流机制发挥主要作用,而随着有效压力的增加,Biot宏观流机制开始占主导,频散和衰减程度逐渐减小.     

含流体孔隙介质中面波的传播特性及应用

基于单相介质中地震波理论的高频面波法已广泛应用于求取浅地表S波的速度.然而水文地质条件表明,普遍的浅地表地球介质富含孔隙.孔隙中充填的流体会显著地影响面波在浅地表的传播,进而造成频散和衰减的变化.本文研究了地震勘探频段内针对含流体孔隙介质边界条件的面波的传播特性.孔隙流体在自由表面存在完全疏通、完全闭合以及部分疏通的情况.孔隙单一流体饱和时,任何流体边界条件下存在R1模式波,与弹性介质中的Rayleigh波类似,相速度稍小于S波并在地震记录中显示强振幅.由于介质的内在衰减,R1在均匀半空间中也存在频散,相速度和衰减在不同流体边界下存在差异.Biot固流耦合系数(孔隙流体黏滞度与骨架渗透率之比)控制频散的特征频率,高耦合系数会在地震勘探频带内明显消除这种差异.介质的迂曲度等其他物性参数对不同流体边界下的R1波的影响也有不同的敏感度.完全闭合和部分疏通流体边界下存在R2模式波,相速度略低于慢P波.在多数条件下,如慢P波在时频响应中难以观察到.但是在耦合系数较低时会显现,一定条件下甚至会以非物理波形式接收R1波的辐射,显示强振幅.浅表风化层低速带存在,震源激发时的运动会显著影响面波的传播.对于接收点径向运动会造成面波的Doppler频移,横向运动会造成面波的时频畸变.孔隙存在多相流体时,中观尺度下不均匀斑块饱和能很好地解释体波在地震频带内的衰减.快P波受到斑块饱和显著影响,R1波与快P波有更明显关联,与完全饱和模型中不同,也更易于等效模型建立.频散特征频率受孔隙空间不同流体成分比例变化的控制,为面波方法探测浅地表流体分布与迁移提供可能性.通常情况孔隙介质频散特征频率较高,标准线性黏弹性固体可以在相对低频的地震勘探频带内等效表征孔隙介质中R1波的传播特征,特别在时域,可在面波成像反演建模中应用.     

孔隙介质弹性波频散—衰减理论模型

储层地球物理学中,孔隙介质的各类弹性波模型常用于了解地层岩石物理性质.本文介绍了含油、气、水等物质的多相孔隙介质弹性波频散和衰减研究进展,给出了流体饱和与部分饱和孔隙介质中波传播的物理模型综述.根据孔隙介质中的固、流体分布情况,从相关基础理论和实验研究工作等方面出发,在宏观、微观和介观尺度上对流体替换、Biot孔隙力学、喷射流、Biot喷射流(BISQ)、等效球体癍块饱和、双重孔隙介质局部流动等现有主要合流体孔隙介质速度频散和衰减理论进行了回顾.研究表明,应力松弛过程是弹性波频散和衰减的基本机理,该过程由平衡特征时间刻画.该特征时间与孔隙介质的渗透率、流体粘性和体积模量紧密相关.当波频较低时,特征时间小于波周期,压力平衡得以发生,可以用等效流体模型描述波速;反之,当波频较高时,局部压差始终保持较高水平。整个骨架体积模量升高,等效模型面临困难,发展出斑块饱和模型.在分析了各类模型理论框架适用性以及所面临困难后,我们对未来研究方向给出了一些有意义的探讨.     

有效应力变化对致密砂岩孔隙结构及弹性波响应的影响规律

致密砂岩气藏具有裂缝发育和有效应力高的特征,研究不同有效压力下孔、裂隙介质地震波传播特征,有利于地震解释与地下储层的识别.但是前人的研究较少考虑岩石内部微观孔隙结构特征与孔隙、裂隙间流体流动的关系.本文首先通过选取四川盆地典型致密砂岩岩样,在不同有效压力下对岩石样本进行超声波实验测量.然后基于实验测得的纵、横波速度进行裂隙参数反演,得到不同有效压力下致密砂岩样本的裂隙孔隙度.再将裂隙孔隙度和样本岩石物理参数代入双重孔隙介质模型,模拟得到不同有效压力下饱水致密砂岩样本纵横波速度频散和衰减的变化规律.结果表明模型预测的速度频散曲线与纵波速度实验测量结果能够较好的吻合.最后统计分析了致密砂岩裂隙参数,得到了致密砂岩储层裂隙参数随有效压力及孔隙度变化特征.依据实际岩石物理参数建立模型,其裂隙参数三维拟合结果能够较好描述致密砂岩裂隙结构与孔隙度、应力的关联,可为实际地震勘探中预测储层裂缝性质提供基础依据.     

The influence of mesoscopic flow on the P-wave attenuation and dispersion in a porous media permeated by aligned fractures

In sedimentary rocks attenuation/dispersion is dominated by fluid-rock interactions. Wave-induced fluid flow in the pores causes energy loss through several mechanisms, and as a result attenuation is strongly frequency dependent. However, the fluid motion process governing the frequency dependent attenuation and velocity remains unclear. We propose a new approach to obtain the analytical expressions of pore pressure, relative fluxes distribution and frame displacement within the double-layer porous media based on quasi-static poroelastic theory. The dispersion equation for a P-wave propagating in a porous medium permeated by aligned fractures is given by considering fractures as thin and highly compliant layers. The influence of mesoscopic fluid flow on phase velocity dispersion and attenuation is discussed under the condition of varying fracture weakness. In this model conversion of the compression wave energy into Biot slow wave diffusion at the facture surface can result in apparent attenuation and dispersion within the usual seismic frequency band. The magnitude of velocity dispersion and attenuation of P-wave increases with increasing fracture weakness, and the relaxation peak and maximum attenuation shift towards lower frequency. Because of its periodic structure, the fractured porous media can be considered as a phononic crystal with several pass and stop bands in the high frequency band. Therefore, the velocity and attenuation of the P-wave show an oscillatory behavior with increasing frequency when resonance occurs. The evolutions of the pore pressure and the relative fluxes as a function of frequency are presented, giving more physical insight into the behavior of P-wave velocity dispersion and the attenuation of fractured porous medium due to the wave-induced mesoscopic flow. We show that the specific behavior of attenuation as function of frequency is mainly controlled by the energy dissipated per wave cycle in the background layer.     

Velocity variations in carbonate rocks due to dual porosity and wave-induced fluid flow

Stiffness variations in carbonates may be described as resulting from different concentrations of flat compliant pores or cracks, which can have a significant effect on the effective stiffness and acoustic properties (e.g., velocities and attenuations) of dry as well as saturated carbonates, although they carry extremely little porosity. As shown in this paper, the effects of dual porosity and wave-induced fluid flow or pore pressure communication may also play a significant role. On the basis of a previously published T-matrix approach to model the effective viscoelastic properties of cracked porous media, we illustrate the (frequency-dependent) effects of wave-induced fluid flow (mainly squirt flow) or pore pressure communication for a model structure consisting of a mixture of fluid-saturated porous grains and fluid-saturated cavities (vugs, etc.) that are embedded in a solid matrix associated with carbonates. We assume that the pores within the porous grains are decoupled from the pores in the solid matrix (and possibly saturated with different fluids) but that each pore system at the micro and/or mesoscale may or may not be connected. For each of four different connectivity models, we present numerical results for four different cases of microstructure (that emphasize the importance of cracks and flat compliant pores). Our numerical results indicate that the velocity and attenuation spectra of carbonates vary significantly, even when the crack density and all other volume concentrations are constant.     

Seismic dispersion and attenuation in layered porous rocks with fractures of varying orientations

Wave‐induced fluid flow plays an important role in affecting the seismic dispersion and attenuation of fractured porous rocks. While numerous theoretical models have been proposed for the seismic dispersion and attenuation in fractured porous rocks, most of them neglect the wave‐induced fluid flow resulting from the background anisotropy (e.g. the interlayer fluid flow between different layers) that can be normal in real reservoirs. Here, according to the theories of poroelasticity, we present an approach to study the frequency‐dependent seismic properties of more realistic and complicated rocks, i.e. horizontally and periodically layered porous rock with horizontal and randomly orienting fractures, respectively, distributed in one of the two periodical layers. The approach accounts for the dual effects of the wave‐induced fluid flow between the fractures and the background pores and between different layers (the interlayer fluid flow). Because C₃₃ (i.e., the modulus of the normally incident P‐wave) is directly related to the P‐wave velocity widely measured in the seismic exploration, and its comprehensive dispersion and attenuation are found to be most significant, we study mainly the effects of fracture properties and the stiffness contrast between the different layers on the seismic dispersion and attenuation of C₃₃. The results show that the increasing stiffness contrast enhances the interlayer fluid flow of the layered porous rocks with both horizontal and randomly orienting fractures and weakens the wave‐induced fluid flow between the fractures and the background pores, especially for the layered porous rock with horizontal fractures. The modelling results also demonstrate that for the considered rock construction, the increasing fracture density reduces the interlayer fluid flow while improves the dispersion and attenuation in the fracture‐relevant frequency band. Increasing fracture aspect ratio is found to reduce the dispersion and attenuation in the fracture‐relevant frequency band only, especially for the layered porous rock with horizontal fractures.     

孔、裂隙并存地层中的声波测井理论及多极子声场特征

在油、气储层的勘探和开发中观察到的一个现象是储层岩石中普遍存在孔隙和裂隙.随着近年来孔、裂隙介质弹性波动理论的进展,我们可以将此理论应用于测井技术,以此来指导从声波测井中测量孔、裂隙地层的声学参数.本文计算了孔、裂隙地层里充流体井眼中的多极子声场,分析了声场随裂隙介质的两个主要参数(即裂隙密度和裂隙纵横比)的变化特征.井孔声场的数值计算表明裂隙密度可以大幅度地降低井中声波纵、横波的波速和振幅.随着裂隙密度的增加,在测井频段内也可以看到纵、横波速的频散现象(这种频散在孔隙地层中一般是观察不到的).本文还研究了多极子模式波 (即单极的Stoneley波、伪瑞利波以及偶极的弯曲波)随裂隙参数的变化特征.结果表明,这些模式波的振幅激发和速度频散都受裂隙密度的影响.裂隙密度越高影响越大.此外,裂隙还对模式波的传播造成较大的衰减.相对裂隙密度而言,裂隙纵横比是一个频率控制参数,它控制裂隙对声场影响的频率区间.本文的分析结果对裂缝、孔隙型地层的声波测井具有指导意义.     

Influence of frequency-dependent anisotropy on seismic amplitude-versus-offset signatures for fractured poroelastic rocks

Frequency-dependent amplitude variation with offset offers an effective method for hydrocarbon detections and analysis of fluid flow during production of oil and natural gas within a fractured reservoir. An appropriate representation for the frequency dependency of seismic amplitude variation with offset signatures should incorporate influences of dispersive and attenuating properties of a reservoir and the layered structure for either isotropic or anisotropic dispersion analysis. In this study, we use an equivalent medium permeated with aligned fractures that simulates frequency-dependent anisotropy, which is sensitive to the filled fluid of fractures. The model, where pores and fractures are filled with two different fluids, considers velocity dispersion and attenuation due to mesoscopic wave-induced fluid flow. We have introduced an improved scheme seamlessly linking rock physics modelling and calculations for frequency-dependent reflection coefficients based on the propagator matrix technique. The modelling scheme is performed in the frequency-slowness domain and can properly incorporate effects of both bedded structure of the reservoir and velocity dispersion quantified with frequency-dependent stiffness. Therefore, for a dispersive and attenuated layered model, seismic signatures represent a combined contribution of impedance contrast, layer thickness, anisotropic dispersion of the fractured media and tuning and interference of thin layers, which has been avoided by current conventional methods. Frequency-dependent amplitude variation with offset responses was studied via considering the influences of fracture fills, layer thicknesses and fracture weaknesses for three classes amplitude variation with offset reservoirs. Modelling results show the applicability of the introduced procedure for interpretations of frequency-dependent seismic anomalies associated with both layered structure and velocity dispersion of an equivalent anisotropic medium. The implications indicate that anisotropic velocity dispersion should be incorporated accurately to obtain enhanced amplitude variation with offset interpretations. The presented frequency-dependent amplitude variation with offset modelling procedure offers a useful tool for fracture fluid detections in an anisotropic dispersive reservoir with layered structures.     

Frequency-dependent PP and PS reflection coefficients in fractured media

When a porous layer is permeated by mesoscale fractures, wave-induced fluid flow between pores and fractures can cause significant attenuation and dispersion of velocities and anisotropy parameters in the seismic frequency band. This intrinsic dispersion due to fracturing can create frequency-dependent reflection coefficients in the layered medium. In this study, we derive the frequency-dependent PP and PS reflection coefficients versus incidence angle in the fractured medium. We consider a two-layer vertical transverse isotropy model constituted by an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce the intrinsic dispersion due to fracturing in the sand layer. Based on the series coefficients that control the behaviour of velocity and anisotropy parameters in the fractured medium at low frequencies, we extend the conventional amplitude-versus-offset equations into frequency domain and derive frequency-dependent amplitude-versus-offset equations at the elastic–anelastic surface. Increase in fracture length or fracture density can enlarge the frequency dependence of amplitude-versus-offset attributes of PP and PS waves. Also, the frequency dependence of magnitude and phase angle of PP and PS reflection coefficients increases as fracture length or fracture density increases. Amplitude-versus-offset type of PP and PS reflection varies with fracture parameters and frequency. What is more, fracture length shows little impact on the frequency-dependent critical phase angle, while the frequency dependence of the critical phase angle increases with fracture density.     

Comparative review of theoretical models for elastic wave attenuation and dispersion in partially saturated rocks

Saturation of porous rocks with a mixture of two fluids (known as partial saturation) has a substantial effect on the seismic waves propagating through these rocks. In particular, partial saturation causes significant attenuation and dispersion of the propagating waves, due to wave-induced fluid flow. Such flow arises when a passing wave induces different fluid pressures in regions of rock saturated by different fluids. As partial fluid saturation can occur on different length scales, attenuation due to wave-induced fluid flow is ubiquitous. In particular, mesoscopic fluid flow due to heterogeneities occurring on a scale greater than porescale, but less than wavelength scale, is responsible for significant attenuation in the frequency range from 10 to 1000 Hz.Most models of attenuation and dispersion due to mesoscopic heterogeneities imply that fluid heterogeneities are distributed in a periodic/regular way. In 1D this corresponds to periodically alternating layering, in 3D as periodically distributed inclusions of a given shape (usually spheres). All these models yield very similar estimates of attenuation and dispersion.Experimental studies show that mesoscopic heterogeneities have less idealized distributions and that the distribution itself affects attenuation and dispersion. Therefore, theoretical models are required which would simulate the effect of more general and realistic fluid distributions.We have developed two theoretical models which simulate the effect of random distributions of mesoscopic fluid heterogeneities. The first model assumes that one fluid forms a random ensemble of spherical inclusions in a porous medium saturated by the other fluid. The attenuation and dispersion predicted by this model are very similar to those predicted for 3D periodic distribution. Attenuation (inverse quality factor) is proportional to ω at low frequencies for both distributions. This is in contrast to the 1D case, where random and periodically alternating layering shows different attenuation behaviour at low frequencies. The second model, which assumes a 3D continuous distribution of fluid heterogeneities, also predicts the same low-frequency asymptote of attenuation. However, the shapes of the frequency dependencies of attenuation are different. As the 3D continuous random approach assumes that there will be a distribution of different patch sizes, it is expected to be better suited to modelling experimental results. Further research is required in order to uncover how to relate the random functions to experimentally significant parameters.