地球物理学报  2021, Vol. 64 Issue (12): 4618-4628   PDF    
基于三维孔隙网络模型的纵波频散衰减特征分析
魏乐乐1, 甘利灯1, 熊繁升2, 孙卫涛2, 丁骞1, 杨昊1     
1. 中国石油勘探开发研究院, 北京 100083;
2. 清华大学周培源应用数学研究中心, 北京 100084
摘要:波场在含流体的孔隙介质中传播时会产生频散和衰减现象.波场的频散和衰减与孔隙介质的岩石物理属性有关,包括孔隙度、渗透率、流体属性等.现有的三维裂缝/软孔隙网络模型利用椭圆截面纵横比的变化模拟从扁裂缝、软孔隙到硬孔隙的多种情况,而未考虑同时包含孔隙和裂缝的全局性网络空间.为了更好地描述裂缝-孔隙空间,本文提出同时包含裂缝和孔隙的三维裂缝-孔隙网络模型,并给出渗透率的计算方法.通过体积平均法推导了三维裂缝/软孔隙网络模型和三维裂缝-孔隙网络模型的波动方程,利用平面波分析方法得到纵波频散/衰减曲线的表达式,同时应用数值模拟研究了总孔隙度、裂缝孔隙度、裂缝纵横比、裂缝数密度、孔隙流体黏度对纵波衰减和速度频散特征的影响.结果表明,在三维裂缝-孔隙网络模型下,总孔隙度、裂缝参数等对纵波频散衰减特征的影响与三维裂缝/软孔隙网络模型相似.具体表现为:纵波在高频段内出现频散和衰减现象.孔隙度的变化主要影响逆品质因子曲线峰值的大小;裂缝数密度主要控制速度显著变化的范围;裂缝纵横比对纵波速度和特征频率有显著影响.
关键词: 三维裂缝/软孔隙网络模型      三维裂缝-孔隙网络模型      波动方程      频散衰减     
Attenuation and dispersion characteristics of P-waves based on the three-dimensional pore network model
WEI LeLe1, GAN LiDeng1, XIONG FanSheng2, SUN WeiTao2, DING Qian1, YANG Hao1     
1. Geophysical Department, RIPED, PetroChina, Beijing 100083, China;
2. Zhou Pei-Yuan Center for Applied Mathematics, Tsinghua University, Beijing 100084, China
Abstract: When propagating in a porous medium containing fluid, the wave field will undergo dispersion and attenuation. Such phenomena are related to the petrophysical properties of the porous medium, including porosity, permeability, fluid properties and so forth. In previous work, based on the three-dimensional fracture/soft pore network model, the variation of the aspect ratio of the ellipse section is used to simulate a variety of conditions from flat fractures to soft pores to hard pores, but does not consider the global network space containing both pores and fractures. To better describe the fracture-pore space, we propose three-dimensional fracture-pore network and the method to calculate the permeability. Through volume averaging, we derive the wave equations of these two three-dimensional pore network models. Then, the expressions of dispersion/attenuation curves of P-waves are obtained by plane wave analysis. And by using numerical simulation, many influence factors are analyzed for the attenuation and dispersion of P-waves, such as total porosity, fracture porosity, fracture aspect ratio, fracture density and pore fluid viscosity. Results show that based on the 3D fracture-pore network model, the effects of total porosity and fracture parameters on the dispersion/attenuation of fast P-waves and the dispersion/attenuation trend of slow P-waves are similar to those of the 3D fracture/soft pore network model. Specifically, the attenuation and dispersion of fast P-waves occur in high-frequency bands. The change of porosity mainly affects the peak value of inverse quality factor curve. The fracture density mainly controls the range of significant change in velocity. And the fracture aspect ratio has a significant effect on the P-wave velocity and characteristic frequency.
Keywords: 3D fracture/soft pore network model    3D fracture-pore network model    Wave equation    Dispersion and attenuation    
0 引言

地震波在通过孔隙介质时会发生频散和衰减现象.研究地震波衰减和速度频散特征对岩性分类、油气识别,以及油气田开发过程中的油气藏管理都具有重要的价值(何润发等,2020凌云,2021赵平起等,2021).由于地下介质条件复杂,地震波在地下含流体孔隙介质中传播时的衰减和频散特征一直是研究的热点问题.

近几十年来,国内外学者开展了大量岩石物理性质的研究,特别是孔隙介质流体的流动引起波的频散衰减过程,发展出许多理论模型.

Biot(1955)建立了完全饱和孔隙介质中波传播的动力学方程.在Biot理论中,黏性的孔隙流体和固体骨架之间的相互作用导致波的速度频散和衰减.Geertsma和Smit(1961)应用Biot理论推导出纵波的速度频散和衰减的近似解.然而,越来越多的研究发现,Biot理论的这种宏观的流体流动机制很难解释实际生产中地震波的高频散和高衰减现象(Carcione et al., 2010).

为了突破Biot理论的局限性,学者们着眼于研究微观尺度的流体流动机制.Mavko和Nur(1979)建立了喷射流模型来描述弹性波在流体部分饱和的裂缝/孔隙介质中的传播,而且认为喷射流动会造成更强的衰减.欧阳芳等(2021)对经典喷射流模型进行了扩展,结合等效介质理论和孔隙结构模型,通过数值模拟研究了微观孔隙结构下的速度频散和衰减特性.

Dvorkin和Nur(1993)认为,在含有流体的岩石中,波的速度频散和衰减除了受Biot机制的影响,还受到喷射流动机制的影响.他们将上述两种机制统一起来,提出了BISQ模型,并推导了纵波相速度和品质因子的表达式,研究了渗透率等对纵波速度频散和衰减的影响.

在Biot理论的框架下,White(1975a, b)首先引入中观尺度的耗散机制,提出非均匀部分饱和的孔隙介质模型,也称为斑状饱和模型.Dutta和Odé(1979ab)通过严格求解Biot方程组,研究了在盐水充填的岩石中,衰减系数随频率、气体饱和度等的变化.Carcione等(2003)深入研究了White模型的物性及流体变化对衰减和频散的影响.王峣钧等(2014)应用斑块饱和岩石物理模型,从理论上分析了不同固结程度岩石中含气饱和度对速度和衰减的影响,结果表明局部含气储层在地震频带内会发生频散和振幅频变效应.在此基础上,他们建立了岩石物理量板,用于估计储层的含气饱和度.

Sun(2021)提出了一种BIPS(Biot-Patchy-Squirt) 模型,用于描述非混相流体饱和裂缝孔隙弹性介质中的波的频散衰减特性.研究结果表明,BIPS模型与实验室数据吻合较好,这一发现将推动上述三种频散衰减机制在波速预测中的潜在应用.

随着研究的深入,学者们开始建立和研究多重孔隙介质模型.巴晶等(2012)建立了双重孔隙结构模型,利用描述非饱和双重孔隙介质地震波传播的Biot-Rayleigh方程求解平面波解,得到纵波、横波的相速度及逆品质因子的表达式,研究了非饱和岩石中的纵波频散与衰减特征.通过对三个地区砂岩储层的分析,发现在地震频带内纵波对储层中含气较为敏感,但对含气饱和度的定量表征效率不高,低孔隙度砂岩的纵波频散和衰减在地震频带更为显著.此后,郭梦秋等(2018)进一步研究了基于此模型的含流体致密砂岩的纵波频散和衰减特征.Ba等(2017)通过分析流体非均质性与岩石组构非均质性的叠加效应,提出了一种双重双孔隙理论.该模型描述了波在不同尺度下的具有组构非均质性的斑状饱和岩石中的传播.建模和实验数据表明,岩石组构的非均质性尺度会造成波的多频频散和衰减效应,特征频率的大小也与岩石组构的非均质性尺度有关.Zhang等(2020, 2021)考虑到岩石组构具有多尺度的非均质性和自相似的分形特征.在此基础上,他们建立了多孔介质模型,通过数值模拟研究了波的频散和衰减特性与分形维数的关系.

在含油气复杂储层勘探时,不可避免的要遇到含裂缝和孔隙的储层岩石,这些岩石孔隙结构复杂、孔隙流体性质多变,导致波场信号出现频散衰减等特征.建立合理的裂缝-孔隙介质模型,是地震正演等技术的重要基础.

Chapman等(2002)提出了在微观尺度上建立孔隙弹性介质模型的方法,将模型设置成规则的立方体网格,每个网络节点代表孔隙/裂缝空间.但是此模型并未建立孔隙介质渗透率与裂缝参数等之间的几何关系.Tang(2011)Tang等(2012)提出了含裂缝孔隙介质的波动方程模型,在推导过程中直接使用了Johnson等(1987)提出的动态渗透率模型,并没有在渗透率与裂缝参数等影响因素之间建立关系,也没有给出波频散衰减与渗透率之间的关系.Song等(2020)假定裂缝的形状可以是矩形的,建立了一个介质模型来估计饱和多孔岩石中包含一个随机分布的、无限小厚度的矩形裂缝时的纵波衰减和速度频散特征.

在天然岩石中,流体流动发生的空间是全局性的裂缝-孔隙网络,而不是局部的单一孔隙或孔管束.Xiong等(2020)熊繁升等(2021)基于岩石内部的裂缝连通网络,用椭圆截面纵横比的变化模拟从扁裂缝、软孔隙到硬孔隙的多种情况,提出了具有椭圆形截面的三维裂缝/软孔隙网络模型,建立了宏观可测观量与裂缝参数之间的关系,并给出了渗透率的计算方法.该模型可以实现孔隙介质中跨尺度的整体建模,能更好地反映孔隙系统的整体效应.

考虑到熊繁升等(2021)建立的流体全饱和的三维裂缝/软孔隙网络模型没有考虑孔隙的存在,因此本文在此模型的基础上建立同时包含裂缝和孔隙的三维裂缝-孔隙网络模型,并推导出渗透率的计算方法.基于三维裂缝/软孔隙网络模型和三维裂缝-孔隙网络模型,运用体积平均的方法推导出基于渗透率的波动方程,利用平面波分析方法得到纵波频散和衰减曲线的表达式,详细研究了总孔隙度、裂缝孔隙度、裂缝纵横比、裂缝数密度以及孔隙流体黏度对快纵波频散衰减特征的影响,以及分析了慢纵波的频散衰减特征.

1 三维孔隙网络模型 1.1 三维裂缝/软孔隙网络模型

熊繁升等(2021)提出三维裂缝/软孔隙网络模型(图 1).模型中单个管道横截面是纵横比可变的椭圆形,可以模拟从圆形孔隙到狭窄裂缝的不同类型的裂缝/软孔隙网络空间.

图 1 三维裂缝/软孔隙网络模型 Fig. 1 3D fracture/soft pore network model

设单元立方体边长为l,沿XYZ方向上分别排布有MNL个椭圆截面管道,用R1R2表示裂缝椭圆截面的长短轴半径.定义如下参数:

(1) 微裂缝的纵横比:a=R2/R1.

(2) 裂缝的孔隙度:

(1)

(3) 单元体内的裂缝数密度:

(2)

熊繁升等(2021)计算三维裂缝/软孔隙介质渗透率的步骤如下:

第一步:采用椭圆柱坐标系,基于质量守恒和NS方程推导出含不可压缩牛顿流体在椭圆截面微管中的流量表达式.

第二步:利用微管两端每个节点的流量守恒条件,得到全部网络节点满足的流量-压力线性方程组.该方程组中未知量为各节点的压力,入口和出口端的压力边界条件是方程的非齐次项.

第三步:求解上述流量-压力线性方程组,得到各节点压力,进而计算得到整个三维裂缝/软孔隙网络的流量.

第四步:结合达西定律计算得到岩石样本的渗透率.

1.2 三维裂缝-孔隙网络模型

在三维裂缝/软孔隙网络模型的基础上,本文建立了同时包含孔隙和裂缝的三维裂缝-孔隙网络模型(图 2).

图 2 三维裂缝-孔隙网络模型 Fig. 2 3D fracture-pore network model

同样地,设单元立方体边长为l,沿XYZ方向上分别排布有MNL个椭圆截面管道.节点上的球形表示孔隙,连接节点的椭圆柱体通道表示裂缝.用rpore表示孔隙半径,Rasp表示平均孔隙半径与平均裂缝截面长轴半径之比.定义连接孔隙ij的裂缝截面长轴半径为

(3)

定义如下参数:

(1) 微裂缝的纵横比a:裂缝椭圆截面的长短轴半径之比的倒数;

(2) 总孔隙度:ϕ=ϕf+ϕpore

其中,裂缝所占孔隙度为

(4)

孔隙空间所占孔隙度为

(5)

式中,NfNpore分别表示裂缝和孔隙总数;AfApore分别表示裂缝和孔隙的截面积,AfaRij2;〈Lij〉表示孔隙ij间距的平均值.

(3) 单元体内的裂缝数密度:表达式同公式(2).

对于三维裂缝-孔隙网络模型的渗透率计算,仍采用1.1节的推导方法.

RiRj表示正圆球形孔隙ij的半径,η为流体黏度,则流体传导系数为,这里r取(Ri+Rj)/2.

对于同时含有裂缝和孔隙网络的介质,考虑任意一段微管,微管两端的网络节点为ij,节点处压力分别是pipj,微管中压力与流量的关系式为

(6)

其中:

(7)

(8)

这里:

(9)

(10)

(11)

式中,aij是截面纵横比;LijRij表示节点ij所连接微管的长度与截面长轴半径;ρf为微管中流体的密度;J是贝塞尔函数;c0是声波在流体中的传播速度,c0=1/(βρf),β为流体压缩系数.

接下来,对各节点列流量平衡方程式,得到全部网络节点满足的流量-压力方程组.求解该方程组得到各节点压力,进而计算得到整个三维裂缝-孔隙网络的流量.最后,结合达西定律计算得到岩石样本的渗透率.

2 基于体积平均法的波动方程推导 2.1 体积平均法

体积平均法是推导波动方程的重要方法.Whitaker(1999)详细介绍了采用体积平均法推导波动方程的原理和步骤.由于该方法通过对微观物理量在单元体上进行积分平均得到宏观物理量,因此可以考虑不同流体类型和孔隙空间的具体几何特征,可以认为是一种更一般化的建模方法.

记单元体Ω的体积为V,孔隙空间由固体骨架的体积和孔隙介质中流体的体积构成,即:

(12)

ψf为孔隙介质中与流体有关的物理量,规定在流体的外部,ψf的值为0.定义在整个区域Ω上对ψf作体积平均的方法为

(13)

定义本征体积平均量:

(14)

根据以上定义式,可知:

(15)

2.2 波动方程推导

在下述的推导过程中:假设流体为牛顿流体,不考虑流体发生正应变过程中黏性的作用,仅考虑剪切黏性;考虑孔隙介质含一种固体骨架、一种流体(或两种流体等效成一种).

uU分别表示固体、流体质点的位移向量;各物理量上方的一点表示对时间求导;以上下标中的sf分别表示固体骨架、孔隙介质中所含的流体.

下面运用体积平均法来推导基于渗透率的三维孔隙网络模型(三维裂缝/软孔隙网络模型与三维裂缝-孔隙网络模型)的波动方程.

第一步,建立固体和流体的微观运动方程(动量守恒方程).

固体骨架部分的微观运动方程为

(16)

流体的微观运动方程为

(17)

式中,ρs为固体骨架的质量密度;σ为固体应力;ρfpf为等效后流体的密度、压力;⊗表示张量积;s为流体应力.

第二步,给出固体部分和流体的本构方程.

固体部分的微观本构方程为

(18)

式中,Ks为固体颗粒的体积模量;G为固体骨架的剪切模量;e为固体的体应变.

流体本构方程为:

(19)

式中,η为等效后流体的黏度.

第三步,对微观方程做体积平均得到宏观大尺度方程.

对流体微观运动方程取体积平均后得到:

(20)

对固体微观运动方程取体积平均后得到:

(21)

式中,κ为孔隙介质的渗透率,可分别根据1.1节和1.2节的计算方法确定.

第四步,推导出流固相互作用的流体和固体的本构关系.

固体和流体的本构关系为

(22)

式中,弹性常数aij可根据固体和流体的孔隙度、体积模量来进行计算.具体表示为

(23)

这里:

(24)

等效体积模量Kb的表达式为(Mavko and Nur, 1978):

(25)

式中,v是固体骨架矿物的泊松比;Nc为单位体内的裂缝数;R1i是第i条裂缝的椭圆截面长轴半径长度;di是裂缝延伸长度;V是单元体的体积.

最后,将本构方程和体积平均后的动量守恒方程结合,得到最终形式的波动方程:

(26)

式中,ξ为流体的体应变.其余参数均沿用上文的定义.

2.3 频散/衰减表征

对波动方程两端取散度,可得到纵波方程.对于各向同性孔隙介质,假设弹性波沿各方向传播一致,利用平面波分析方法可得到纵波频率-波数方程.将纵波频率-波数方程的解记为v,根据如下定义,即可计算纵波频散/衰减曲线:

(27)

3 频散/衰减特征分析

波场在含流体的孔隙介质中传播时会产生频散和衰减现象,即波速随着频率发生改变,同时波的振幅也会发生变化.

下面基于三维裂缝/软孔隙网络模型和三维裂缝-孔隙网络模型分别研究纵波频散与衰减特征.

3.1 基于三维裂缝/软孔隙网络模型的频散衰减特征

利用前面所述频散衰减曲线的表达式,根据岩石物理参数可计算得到基于三维裂缝/软孔隙网络模型的频率依赖的纵波速度和逆品质因子曲线.数值模拟所选取的模型、矿物和流体参数具体见表 1.

表 1 模型、矿物及流体参数表 Table 1 Parameters of model, mineral and pore fluid
3.1.1 裂缝孔隙度的影响

图 3图 4给出了不同裂缝孔隙度条件下,快纵波的频散和衰减曲线.本算例中除裂缝孔隙度外,其余参数见表 1.

图 3 快纵波频散曲线随裂缝孔隙度的变化 Fig. 3 Velocity dispersion changing with fracture porosity
图 4 快纵波衰减曲线随裂缝孔隙度的变化 Fig. 4 Inverse quality factor versus fracture porosity

图 3图 4中四条曲线的裂缝孔隙度分别为0.001、0.005、0.02和0.06,对应的渗透率分别为7.00×10-17m2、1.57×10-16m2、3.14×10-16m2和5.44×10-16m2.从图 3可以看出,随着裂缝孔隙度的减小,快纵波速度增大,起始速度由3116 m·s-1增大到4863 m·s-1.而且随着裂缝孔隙度的增大,快纵波速度表现出更明显的频散现象.从图 4可以看出,随着裂缝孔隙度的减小,逆品质因子曲线的峰值下降,特征频率(逆品质因子曲线的峰值对应的频率)向低频方向移动.

3.1.2 裂缝纵横比的影响

图 5图 6给出了不同裂缝纵横比条件下,快纵波的频散和衰减曲线.本算例中除裂缝纵横比外,其余参数见表 1.

图 5 快纵波频散曲线随裂缝纵横比的变化 Fig. 5 Velocity dispersion versus aspect ratio of fracture
图 6 快纵波衰减曲线随裂缝纵横比的变化 Fig. 6 Attenuation changing with aspect ratio of fracture

图 5图 6中四条曲线的裂缝纵横比分别为0.02、0.04、0.08和0.16,对应的渗透率分别为5.44×10-16m2、3.07×10-15m2、1.73×10-14m2和9.60×10-14m2.通过计算可知,随着裂缝纵横比的下降,渗透率下降.这是一个合理的事实,因为随着裂缝纵横比的减小,裂缝趋于闭合.图 5显示,随着裂缝纵横比的增大,快纵波速度增大.这主要是因为随着纵横比的增大,干岩石的可压缩性变差,干岩石的体积模量Kb增大,从而引起快纵波速度的增大.图 6显示,随着裂缝纵横比的增大,特征频率向低频方向移动.

3.1.3 裂缝数密度的影响

图 7图 8给出了不同裂缝数密度条件下,快纵波的频散和衰减曲线.本算例中除网络节点数外,其余参数见表 1.

图 7 快纵波频散曲线随裂缝数密度的变化 Fig. 7 Velocity dispersion varying with fracture density
图 8 快纵波衰减曲线随裂缝数密度的变化 Fig. 8 Attenuation varying with fracture density

沿XYZ方向上分别排布有3×3×3、4×4×4、5×5×5和6×6×6个孔隙,即裂缝数密度分别为54、144、300和540.计算得到对应的渗透率分别为1.40×10-16m2、3.11×10-16m2、5.44×10-16m2和8.39×10-16m2.通过计算可知,渗透率随裂缝数密度的增加而增大.这表明随着裂缝之间的连通程度升高,岩石渗透率逐步增大.图 7显示,频散幅值(速度极大值与极小值的差)对裂缝数密度的变化不敏感,但频散曲线随裂缝数密度的变化发生移动.在相同的频率上,快纵波速度最大相差仅10 m·s-1.图 8显示,随着裂缝数密度的增大,特征频率向低频方向移动,而几乎不影响逆品质因子曲线的峰值.

3.1.4 孔隙流体黏度的影响

下面研究孔隙流体的黏度对快纵波频散衰减特征的影响.孔隙中的流体参数见表 2,其余参数同表 1.

表 2 孔隙流体参数表 Table 2 Parameters of pore fluid

图 9图 10可以看出,在此模型下的快纵波速度:油饱和砂岩明显大于水饱和砂岩,水饱和砂岩明显大于气饱和砂岩.这主要是因为孔隙流体黏度增大时,孔隙流体与岩石之间的耦合作用加强,相当于提高了岩石的总体刚度,从而提高了快纵波的速度.但油饱和时的速度变化率相对于水饱和时变小.相较于孔隙充填其他流体,气饱和时的快纵波频散现象最不明显,且逆品质因子曲线的峰值最小.

图 9 不同孔隙流体黏度下的快纵波频散曲线 Fig. 9 Velocity dispersion curve as function of fluid viscosity
图 10 不同孔隙流体黏度下的快纵波衰减曲线 Fig. 10 Attenuation curve under different fluid viscosity
3.1.5 慢纵波频散衰减曲线

模型内不仅存在快纵波,还存在慢纵波.将表 1中的岩石物理参数代入公式(27),得到慢纵波的频散衰减曲线,如图 11图 12所示.

图 11 慢纵波频散曲线 Fig. 11 Velocity dispersion curve of slow P-waves
图 12 慢纵波衰减曲线 Fig. 12 Attenuation curve of slow P-waves

图 11可以看出,在低频范围内(0~102 Hz)慢纵波速度基本为0,且在0~108 Hz频率范围内发生频散,速度随频率的增加而增大,在高频范围内稳定在大约740 m·s-1.从图 12可以看出,慢纵波逆品质因子曲线在0~106 Hz的频率范围内为水平线,后随着频率的增加而减小,呈直线下降趋势.

3.2 基于三维裂缝-孔隙网络模型的频散衰减特征

根据式(3)知,裂缝椭圆截面长轴半径与裂缝两端的孔隙半径有关.因此,裂缝的孔隙度ϕf和孔隙空间的孔隙度ϕpore不是互相独立的.

在实际的计算中,给定ϕpore,首先根据式(5)计算得到孔隙半径,后根据式(3)计算得到裂缝截面的长轴半径,再根据裂缝纵横比得到裂缝的短轴半径.将上述计算得到的裂缝截面的长轴半径代入式(4),计算得到ϕf.

为了满足预设的总孔隙度ϕ,可能需要对孔隙半径和裂缝半径进行重新处理.因此,在三维裂缝-孔隙网络模型中,总孔隙度ϕ与预设结果相同,但是ϕporeϕf可能与预设结果不同.

3.2.1 总孔隙度的影响

除裂缝孔隙度外,仍采用表 1中的岩石物理参数.图 13图 14给出了总孔隙度分别为0.001、0.005、0.02和0.06时的快纵波频散与衰减曲线.通过计算得到对应的渗透率为3.12×10-17m2、7.00× 10-17m2,1.40×10-16m2和2.22×10-16m2.从图中可以看出,随着总孔隙度的减小,快纵波速度增大,起始速度由3559 m·s-1增大到4367 m·s-1,逆品质因子曲线的峰值下降,特征频率向低频方向移动.

图 13 快纵波频散曲线随总孔隙度的变化 Fig. 13 Velocity dispersion curve as function of total porosity
图 14 快纵波衰减曲线随总孔隙度的变化 Fig. 14 Attenuation curve as function of total porosity
3.2.2 裂缝参数的影响

从建模过程中可以发现,三维裂缝/软孔隙网络模型可以认为是三维裂缝-孔隙网络模型的特例,即当孔隙空间所占孔隙度为0时,三维裂缝-孔隙网络模型退化为三维裂缝/软孔隙网络模型.若固定孔隙空间所占孔隙度为0,分析裂缝参数对快纵波频散衰减特征的影响,即为文章3.1节的内容.

3.2.3 慢纵波频散衰减曲线

表 1中的岩石物理参数代入公式(27),得到慢纵波的频散衰减曲线.

图 15图 16可以看出,在该模型下慢纵波频散与衰减的趋势与三维裂缝/软孔隙网络模型相似,但频散的幅值更大,发生频散的频率范围更宽,在高频范围内慢纵波的速度稳定在大约1214 m·s-1.

图 15 慢纵波频散曲线 Fig. 15 Velocity dispersion curve of slow P-waves
图 16 慢纵波衰减曲线 Fig. 16 Attenuation curve of slow P-waves
4 结论

本文基于三维裂缝/软孔隙网络模型和三维裂缝-孔隙网络模型,详细研究了总孔隙度、裂缝孔隙度、裂缝纵横比、裂缝数密度、流体黏度对快纵波频散衰减特征的影响,以及分析了慢纵波的频散衰减特征.

基于三维裂缝/软孔隙网络模型,可以发现:

(1) 随着裂缝孔隙度的增大,快纵波速度表现出更明显的频散现象.随着裂缝孔隙度的减小,快纵波速度增大,逆品质因子曲线的峰值下降,特征频率(逆品质因子曲线的峰值对应的频率)向低频方向移动.

(2) 随着裂缝纵横比的增大,快纵波速度增大,特征频率向低频方向移动,而频散幅值(速度极大值与极小值的差)和逆品质因子曲线的峰值的大小未表现出规律性的变化.

(3) 裂缝数密度的变化几乎不会影响快纵波速度、频散幅值和逆品质因子曲线的峰值,但对特征频率影响显著.随裂缝数密度的增大,特征频率向低频方向移动.

(4) 当孔隙充填油、水、气时,油饱和砂岩的快纵波速度明显大于水饱和砂岩,水饱和砂岩明显大于气饱和砂岩.气饱和时的快纵波频散现象最不明显,且逆品质因子曲线的峰值最小.

基于三维裂缝-孔隙网络模型,可以发现:

在此模型下,总孔隙度、裂缝参数等对快纵波频散衰减特征的影响以及慢纵波的频散衰减趋势与三维裂缝/软孔隙网络模型相似.

总的来说,孔隙度的变化主要影响逆品质因子曲线的峰值的大小;裂缝数密度主要控制速度显著变化的范围;裂缝纵横比的变化对纵波速度和特征频率都影响显著.

References
Ba J, Carcione J M, Cao H, et al. 2012. Velocity dispersion and attenuation of P waves in partially-saturated rocks: wave propagation equations in double-porosity medium. Chinese Journal of Geophysics (in Chinese), 55(1): 219-231. DOI:10.6038/j.issn.0001-5733.2012.01.021
Ba J, Xu W H, Fu L Y, et al. 2017. Rock anelasticity due to patchy saturation and fabric heterogeneity: a double double-porosity model of wave propagation. Journal of Geophysical Research: Solid Earth, 122(3): 1949-1976.
Biot M A. 1955. Theory of propagation of elastic waves in a fluid-saturated porous solid. Ⅰ. low-frequency range. The Journal of the Acoustical Society of America, 28(182): 168-191.
Carcione J M, Helle H B, Pham N H. 2003. White's model for wave propagation in partially saturated rocks: comparison with poroelastic numerical experiments. Geophysics, 68(4): 551-566.
Carcione J M, Morency C, Santos J E. 2010. Computational poroelasticity-a review. Geophysics, 75(5): 75A229-75A243. DOI:10.1190/1.3474602
Chapman M, Zatsepin S V, Crampin S. 2002. Derivation of a microstructural poroelastic model. Geophysical Journal International, 151(2): 427-451. DOI:10.1046/j.1365-246X.2002.01769.x
Dutta N C, Odé H. 1979a. Attenuation and dispersion of compressional waves in fluid-filled porous rocks with partial gas saturation (White model)-Part Ⅰ: Biot theory. Geophysics, 44(11): 1777-1788. DOI:10.1190/1.1440938
Dutta N C, Odé H. 1979b. Attenuation and dispersion of compressional waves in fluid-filled porous rocks with partial gas saturation (White model)-Part Ⅱ: Results. Geophysics, 44(11): 1789-1805. DOI:10.1190/1.1440939
Dvorkin J, Nur A. 1993. Dynamic poroelasticity: a unified model with the squirt and the Biot mechanisms. Geophysics, 58(4): 524-533. DOI:10.1190/1.1443435
Geertsma J, Smit D C. 1961. Some aspects of elastic wave propagation in fluid-saturated porous solids. Geophysics, 26(2): 169-181. DOI:10.1190/1.1438855
Guo M Q, Ba J, Ma R P, et al. 2018. P-wave velocity dispersion and attenuation in fluid-saturated tight sandstones: characteristics analysis based on a double double-porosity structure model description. Chinese Journal of Geophysics (in Chinese), 61(3): 1053-1068. DOI:10.6038/cjg2018L0678
He R F, Ba J, Du Q Z, et al. 2020. Seismic wave propagation theory and fluid prediction methods in two-phase media of reservoirs. Progress in Geophysics (in Chinese), 35(4): 1379-1390. DOI:10.6038/pg2020DD0177
Johnson D L, Koplik J, Dashen R. 1987. Theory of dynamic permeability and tortuosity in fluid-saturated porous media. Journal of Fluid Mechanics, 176: 379-402. DOI:10.1017/S0022112087000727
Ling Y. 2021. Well-seismic dispersion correction and its application: a case study of deep-water sandstone in West Africa. Progress in Geophysics (in Chinese), 36(4): 1554-1559. DOI:10.6038/pg2021EE0312
Mavko G M, Nur A. 1978. The effect of nonelliptical cracks on the compressibility of rocks. Journal of Geophysical Research: Solid Earth, 83(B9): 4459-4468. DOI:10.1029/JB083iB09p04459
Mavko G M, Nur A. 1979. Wave attenuation in partially saturated rocks. Geophysics, 44(2): 161-178. DOI:10.1190/1.1440958
Ouyang F, Zhao J G, Li Z, et al. 2021. Modeling velocity dispersion and attenuation using pore structure characteristics of rock. Chinese Journal of Geophysics (in Chinese), 64(3): 1034-1047. DOI:10.6038/cjg2021O0355
Song Y J, Hu H S, Han B. 2020. P-wave attenuation and dispersion in a fluid-saturated rock with aligned rectangular cracks. Mechanics of Materials, 147: 103409. DOI:10.1016/j.mechmat.2020.103409
Sun W T. 2021. On the theory of Biot-patchy-squirt mechanism for wave propagation in partially saturated double-porosity medium. Physics of Fluids, 33(7): 076603. DOI:10.1063/5.0057354
Tang X M. 2011. A unified theory for elastic wave propagation through porous media containing cracks-an extension of Biot's poroelastic wave theory. Science China Earth Sciences, 54(9): 1441-1452. DOI:10.1007/s11430-011-4245-7
Tang X M, Chen X L, Xu X K. 2012. A cracked porous medium elastic wave theory and its application to interpreting acoustic data from tight formations. Geophysics, 77(6): D245-D252. DOI:10.1190/geo2012-0091.1
Wang Y J, Chen S Q, Wang L, et al. 2014. Gas saturation analysis with seismic dispersion attribute based on patchy-saturation model. Oil Geophysical Prospecting (in Chinese), 49(4): 715-722.
Whitaker S. 1999. Dispersion in porous media. //The Method of Volume Averaging. Dordrecht: Springer: 124-160.
White J E. 1975a. Computed seismic speeds and attenuation in rocks with partial gas saturation. Geophysics, 40(2): 224-232. DOI:10.1190/1.1440520
White J E. 1975b. Low-frequency seismic waves in fluid-saturated layered rocks. The Journal of the Acoustical Society of America, 57(S1): S30. DOI:10.1121/1.1995164
Xiong F S, Sun W T, Ba J, et al. 2020. Effects of fluid rheology and pore connectivity on rock permeability based on a network model. Journal of Geophysical Research: Solid Earth, 125(3): e2019JB018857. DOI:10.1029/2019JB018857
Xiong F S, Gan L D, Sun W T, et al. 2021. Characterization of reservoir permeability and analysis of influencing factors in fracture-pore media. Chinese Journal of Geophysics (in Chinese), 64(1): 279-288. DOI:10.6038/cjg2021N0175
Zhang L, Ba J, Carcione J M, et al. 2020. Differential poroelasticity model for wave dissipation in self-similar rocks. International Journal of Rock Mechanics and Mining Sciences, 128: 104281. DOI:10.1016/J.IJRMMS.2020.104281
Zhang L, Ba J, Carcione J M. 2021. Wave propagation in infinituple-porosity media. Journal of Geophysical Research: Solid Earth, 126(4): e2020JB021266. DOI:10.1029/2020JB021266
Zhao P Q, Ni T L, Zhang J L, et al. 2021. Review of the correlation between seepage field parameters and seismic wave field parameters. Progress in Geophysics (in Chinese), 36(4): 1661-1668. DOI:10.6038/pg2021EE0478
巴晶, Carcione J M, 曹宏, 等. 2012. 非饱和岩石中的纵波频散与衰减: 双重孔隙介质波传播方程. 地球物理学报, 55(1): 219-231. DOI:10.6038/j.issn.0001-5733.2012.01.021
郭梦秋, 巴晶, 马汝鹏, 等. 2018. 含流体致密砂岩的纵波频散及衰减: 基于双重双重孔隙结构模型描述的特征分析. 地球物理学报, 61(3): 1053-1068. DOI:10.6038/cjg2018L0678
何润发, 巴晶, 杜启振, 等. 2020. 储层双相介质地震波传播理论及流体预测方法. 地球物理学进展, 35(4): 1379-1390. DOI:10.6038/pg2020DD0177
凌云. 2021. 井震频散校正及其应用: 以西非深水砂岩为例. 地球物理学进展, 36(4): 1554-1559. DOI:10.6038/pg2021EE0312
欧阳芳, 赵建国, 李智, 等. 2021. 基于微观孔隙结构特征的速度频散和衰减模拟. 地球物理学报, 64(3): 1034-1047. DOI:10.6038/cjg2021O0355
王峣钧, 陈双全, 王磊, 等. 2014. 基于斑块饱和模型利用地震波频散特征分析含气饱和度. 石油地球物理勘探, 49(4): 715-722.
熊繁升, 甘利灯, 孙卫涛, 等. 2021. 裂缝-孔隙介质储层渗透率表征及其影响因素分析. 地球物理学报, 64(1): 279-288. DOI:10.6038/cjg2021N0175
赵平起, 倪天禄, 张家良, 等. 2021. 渗流场参数与地震波场参数的关联关系研究综述. 地球物理学进展, 36(4): 1661-1668. DOI:10.6038/pg2021EE0478