非均匀饱含黏性流体孔隙介质中声波传播及井孔声场分析*

2023-03-17彭凡张秀梅刘琳王秀明

物理学报 2023年5期

彭凡张秀梅† 刘琳王秀明

1) (中国科学院声学研究所,声场声信息国家重点实验室,北京 100190)

2) (中国科学院大学,北京 100049)

3) (中国科学院声学研究所,北京市海洋深部钻探研究中心,北京 100190)

声波在饱含流体孔隙介质中的传播特性与流体的黏滞性及孔隙介质的非均匀性密切相关.本文在Biot 理论基础上,考虑了孔隙流体的剪切应力及孔隙结构的非均匀性,采用含黏性流体孔隙介质中的波动理论,研究了孔隙介质中四种体波的频散和衰减特性,分析了慢横波对快纵波转换散射的影响,进一步推导了孔隙地层井孔中的模式波及其声场的解析解,研究了非均匀孔隙介质中井孔模式波和波列的特征.研究结果表明,含黏性流体孔隙介质中存在慢横波,慢横波的频散很强,其传播特征受到介质孔隙度、渗透率及孔隙流体黏度的影响.在非均匀孔隙介质中,与慢横波相关的剪切应力平衡过程不仅导致快纵波的频散和衰减,还会影响井孔伪瑞利波及斯通利波的传播特征.本文的工作完善了孔隙介质中声波传播的物理机制,为孔隙地层井孔声波的解释与应用提供了理论指导.

1 引言

孔隙介质一般是由固体骨架和充填在骨架间的孔隙流体构成,在地下富含油气的储层、工程生产和日常生活等涉及的材料中广泛分布.因此,研究声波在孔隙介质中的传播规律在油气勘探、材料工程及无损检测等方面具有重要作用[1-3].Biot[4-6]早年建立的单一流体饱和的孔隙介质弹性波动力学理论,又称为Biot 理论,目前已成为研究孔隙介质中声波传播的基本理论.Rosenbaum[7]利用Biot理论研究了孔隙地层包围的井孔中声波的传播,形成了Biot-Rosenbaum 理论.应用这些理论,学者们更加深入地研究了井孔声波响应及其与孔隙地层中的岩性和物性参数之间的联系,发展了声波测井数据处理方法并建立了岩石储层物理解释模型.例如: 建立了利用声波时差评估地层孔隙度的公式;提出了利用井孔模式波反演地层渗透率、利用测井声波的衰减获取孔隙流体信息等方法[8-10].然而,理论和实际应用研究表明,基于Biot 理论得到的一些结果与实际观测结果存在很大的差异,多数情况下低估了实际孔隙介质中的频散与衰减[11-13].究其原因之一,Biot 理论假设所研究的孔隙介质具有宏观均匀特征,且孔隙中的填充流体为理想流体,这两个方面的假设往往与实际介质的特性不符.自然界孔隙介质中的流体往往具有黏性应力,并且孔隙介质中普遍存在孔隙结构非均匀的特性[14,15].Biot 理论忽视了这两方面的影响,这可能使得相关的理论和方法在描述实际问题时呈现明显的不足.随着非常规油气勘探开发的深入,对复杂孔隙地层井孔声场理论的改善越来越重要,这将为井孔声波的处理和解释开辟新的应用前景.

孔隙流体黏滞性影响方面,Liu 和 Katsube[16]提出由于流体的黏滞效应,孔隙介质中固、液两相除相对平动外,还存在相对转动,因此,固液分界面上会出现一种新横波,它的高频散、强衰减也是孔隙介质中波传能量耗散的原因之一.魏修成等[17]基于Biot 理论,提出了含黏性流体VTI孔隙介质模型,理论上推导出两种新的慢SV 波和慢SH 波的存在.Gao[18]基于Biot 理论,从液体黏弹性本构关系出发,建立了黏性多孔模型,该模型中会出现一种新的慢横波,主要在流体中传播.Sahay[19]从孔隙尺度的本构关系出发,利用体积平均法和牛顿流变理论,推导出了均匀孔隙介质中黏性修正的Biot 理论(viscosity-extended Biot theory,VEB理论),预测了慢横波的存在并描述了其传播特性.

在孔隙介质的非均匀性的影响方面,也有很多学者做了相关工作.由于孔隙内部结构的非均匀性,孔隙流体在波激励下将发生局部流体流动,造成弹性波能量的损失.在Biot 理论的框架内,这种机制可以理解为快纵波散射转化为慢纵波的过程[20,21].Mavco 和Nur[22]基于单个孔隙的非均匀结构提出了喷射流理论.Dvorkin 和 Nur[23]同时考虑了含流体岩石内部的Biot 流与喷射流,提出了BISQ 理论.崔志文等[24]基于BISQ 模型研究了非牛顿流体饱和孔隙介质中弹性波的频散与衰减.Berryman和Wang[25]采用双孔模型研究了双孔双渗介质中弹性波的频散与衰减.在真实地球材料中,孔隙结构的非均质性常以随机方式在空间上分布,随机方法具有更广泛的适用性[26,27].Müller 和Gurevich[21]采用随机分布函数研究了快纵波到慢纵波的转换散射效应,利用统计波传播理论分析了相干波场的动态特性.Müller 和Sahay[28,29]将该方法扩展到含黏性流体孔隙介质中,基于VEB 理论研究了慢横波的散射对快纵波及慢纵波的影响,该过程涉及孔隙流体黏性边界层内部的剪应力平衡,与孔隙尺度的非均质性有关.

综上所述,针对Biot 理论中未考虑流体黏滞性及非均匀性的影响,前人分别开展了一系列的工作.其中,VEB 理论从孔隙尺度出发引入孔隙流体的黏性应力,弥补了Biot 理论的本构关系中对流体黏性应力项的缺失;同时,随机转换散射理论可以分析由非均匀孔隙结构导致的散射衰减,解决了Biot 理论未考虑介质非均匀性的问题.因此,本文的工作综合这两方面工作的优势,深入地研究了流体黏滞性及孔隙结构的非均匀性对孔隙介质中声波传播的影响,并进一步研究了井孔声场中模式波及全波的传播特征.本文的基本结构如下: 首先,描述了VEB 理论的固相与流相的本构关系,结合运动方程,利用平面波分析的方法计算了含黏性流体孔隙介质中四种体波的波数,着重分析了慢横波的传播特征随孔隙介质参数的变化;接着,针对非均匀孔隙结构的影响,利用随机转换散射理论,分析了不同非均匀孔隙介质中慢横波的扩散对弹性波的影响;最后,推导了含黏性流体孔隙介质井孔声场理论,计算了非均匀孔隙地层中井孔模式波的频散特征及全波波形,分析了孔隙流体黏度对井孔声波的影响.

2 含黏性流体孔隙介质波动理论

2.1 VEB 理论

自然界孔隙介质中的流体多为牛顿流体,不仅存在流体压力,还存在黏性应力.Biot 理论中的本构方程忽略了流体黏性应力.为此,Sahay[19]从孔隙尺度出发,提出了VEB 理论,推导出宏观均匀含黏性流体孔隙介质的本构关系:

与Biot 理论对比可发现,VEB 理论中孔隙流体的本构方程多了与黏滞流体速率相关的流体剪切应力项.对于饱含流体孔隙介质的研究,体积黏度可忽略,因此,本文主要考虑流体剪切应力的影响.将本构关系((1)式—(3)式)代入流体饱和孔隙介质的运动方程[6]:

式中,τ为总应力;地层密度ρ=ρfη+ρs(1-η),ρf为孔隙流体密度,ρs为骨架密度;k0为静态渗透率(单位为D,1 D=0.986923 × 10—12m2),附加质量ρa=(S-1)ρfη;S为弯曲度;u为骨架位移矢量,w=η(U -u) 为渗流位移矢量,U为孔隙流体位移矢量,上标点号表示对时间求导.假设位移矢量依 e-iωt随时间变化,可以得到以u和w为基本量的弹性动力学方程:

在横波波数的求解方法中,利用渗流位移矢量求得的横波波数存在缺失,因此对于横波波数的求解使用的是(8)式中的流体位移矢量,得到横波波数满足的方程为:

对方程进行求解可以得到两种纵波和两种横波的波数.

2.2 体波频散与衰减分析

为了分析含黏性流体孔隙介质中体波的传播特征,根据2.1 节中求得的体波波数,进一步通过下式计算体波的相速度和衰减:

其中ki表示四种体波的波数;V为相速度;Q -1为衰减.图1 展示了VEB 理论与Biot 理论中体波的频散与衰减曲线,计算时所用的相关参数如表1 所列[30].在本文后续的计算中,若无特别说明,选取的基础参数均如表1 所列.由图1(a)—(c)可以看出,均匀孔隙介质中,VEB 理论计算的快横波、快纵波、慢纵波与Biot 理论一致.此外,由于VEB理论中考虑了流体的黏性应力,因此还多了一种慢横波(图1(d)).慢横波的频散很强,其相速度随频率的增大而增大,衰减随频率的增大先减小,当频率到达Biot 弛豫频率(10.6 kHz)后,衰减系数保持不变,此时慢横波为扩散波.Biot 弛豫频率是划分黏性主导区与惯性主导区的特征频率[5].慢横波的相速度小,波长短,它在其产生源附近迅速衰减,因此很难被直接地观测到.

表1 孔隙介质参数Table 1.Parameters of porous media.

图1 四种体波的相速度与衰减 (a) 快纵波;(b) 快横波;(c) 慢纵波;(d) 慢横波Fig.1.Phase velocity and attenuation of four kinds of body waves: (a) Fast P-wave;(b) fast S-wave;(c) slow P-wave;(d) slow S-wave.

为了进一步研究慢横波在不同孔隙介质中的传播特征,接着分析了孔隙介质物性参数(孔隙度、渗透率、孔隙流体黏度)的变化对慢横波的影响,孔隙度分别取0.1,0.2,0.3;渗透率分别取0.01,0.10 和1.00 D;流体黏度分别取0.0001,0.0010 和0.0100 Pa·s,结果如图2 所示.

从图2 中可以看出,慢横波的相速度随孔隙度的减小而增大,随渗透率减小而减小,但高达一定频率后,孔隙度和渗透率的变化对慢横波的相速度没有影响.低于一定频率时(约10 kHz),黏度越大,慢横波相速度速度越小,高于一定频率时(约80 kHz),黏度越大,相速度越大.对于慢横波的衰减,中低频率下,孔隙度越小,渗透率越大,流体黏度越小,衰减越小.在图2(c)中,频率为10 MHz时,慢横波的速度可达10 m/s.经计算,频率更高时,慢横波的速度可达到几十上百米每秒,这与卢明辉等[31]的研究结果相近.因此,在多孔材料的无损检测中,检测频率较高时,更容易检测到慢横波并进一步研究其传播特性.

图2 慢横波的相速度与衰减随孔隙介质参数的变化 (a) 孔隙度;(b) 渗透率;(c) 流体黏度Fig.2.Variation of phase velocity and attenuation of slow shear wave with porous media parameters: (a) Porosity;(b) permeability;(c) fluid viscosity.

2.3 非均匀介质中慢横波对纵波的影响

孔隙介质中,流体的速度场由两部分组成: 势流场与黏性边界层流场[32].储层岩石是常见的孔隙介质,实验室测量表明,由于不同的粒度和流体流动通道,岩石存在微观尺度上的不均匀性[33].非均匀性将导致黏性边界层中波诱导的涡度扩散,黏性流体中的剪切应力通过涡度扩散平衡,并引起的弹性波衰减和频散,该过程可视为弹性波到慢剪切波的转换散射过程[34].Müller 和Sahay[28]将非均匀孔隙介质假设为随机介质,基于VEB 理论,利用随机过程理论和统计平滑的方法,研究了快纵波到慢横波的转换散射过程,导出了快纵波有效波数表达式:

图3 均匀及不同非均匀孔隙介质中快纵波的相速度与衰减Fig.3.The phase velocity and attenuation of fast P wave in homogeneous porous medium and different inhomogeneous porous media.

除了与快纵波的相互作用,慢纵波与慢横波也将产生转换散射作用.该过程主要发生在频率高于Biot 弛豫频率时,此时慢纵波的衰减较小,为传播的波,相应有效慢纵波的表达式与(13)式类似.慢纵波主要在流体中传播,受到随频率变化的孔隙流体的运动状态的影响.Müller 和Sahay[29]将慢纵波对慢横波的散射效应转换为动态渗透率,提出随机动态渗透率模型:

图4 k(ω)的归一化幅值与相位Fig.4.Normalized amplitude and phase of k (ω) .

3 非均匀孔隙地层井孔声场

声波在井孔中的传播是声学测井的物理基础,能否获取准确的地层信息,关键在于理论分析得到孔隙介质井孔声场与实际情况是否相符.以往基于Biot 理论的井孔声场理论未考虑孔隙流体黏性应力及地层非均质性的影响,存在一定的误差[35,36].本节根据VEB 理论,推导了含黏性流体的孔隙介质地层中单极子源井孔声场的解析解,并研究了非均匀孔隙介质中井孔模式波及全波波形的特征.

3.1 井孔声场理论

在柱坐标系(r,θ,z)下,建立了含黏性流体的非均匀孔隙地层井孔模型,如图5 所示,井孔半径a为0.1 m,声场由位于井轴上的单极子声源激发,井内充满理想流体,井外是无限大黏性流体饱和的非均匀孔隙地层.

图5 非均匀孔隙介质地层井孔模型Fig.5.Borehole model of heterogeneous porous media formation.

由于井内为理想流体,因此井内声场只有压缩波,对应的势函数为φ.而井外孔隙介质地层中不仅含有压缩波,还有剪切波,因此井外孔隙地层的势函数为φ p1,φp2,ψs1,ψs2.结合轴对称井孔问题求解方法,得到井内总声场表达式如下:

其中k为轴向波数;表示井孔流体的径向波数;kf为井孔流体波数,Km(x),Im(x)为Bessel 函数.同理,井外无源区域的声场势函数分别为

流体轴向剪应力为零

m矩阵中元素的表达式见附录.解此方程,可得Af,Ap1,Ap2,As1,As2.至此,井内、井外声场就完全确定了.

3.2 井孔模式波及波列特征

为了分析非均匀孔隙介质中,慢横波的散射对井孔声场的影响,本节将基于VEB 理论的有效快纵波波数((13)式)和动态渗透率模型((14)式)代入声场的计算中,分析三种非均匀程度的孔隙地层中模式波及全波的特征.首先令m矩阵的行列式为零,可得到模式波波数满足的频散方程.单极子声源在快速地层中可以激发出伪瑞利波和斯通利波,伪瑞利波存在截止频率.图6(a)、图6(b)和图6(c)分别显示了斯通利波,一阶伪瑞利波,二阶伪瑞利波在不同非均匀程度的孔隙地层中的相速度与衰减.由图6(a)、图6(b)和图6(c)中可观测到,与均匀情况相比,模式波的相速度变化不大,且非均匀介质中的相速度略大于均匀介质的相速度.此外,衰减的变化与相速度相比更为明显,且非均匀介质中斯通利波与伪瑞利波的衰减比均匀介质中的小,但是孔隙介质非均匀性程度的大小对模式波特征的影响不大.在频率为10 kHz 下,斯通利波在均匀介质中的衰减与其在非均匀程度为0.15 的介质中的衰减相比相差17.0%,一阶伪瑞利波的衰减在两种介质下相比相差36.8%.模式波的变化与孔隙地层中的慢纵波的传播密切相关.由于孔隙地层的渗透性,井孔模式波的能量被地层慢纵波带入孔隙地层中,从而造成模式波的频散和衰减,渗透性越大,模式波的衰减越强.慢纵波与模式波的相互作用与动态渗透率有关,非均匀孔隙介质中,动态渗透率的变化受到由慢纵波到慢横波的转换散射过程控制,进一步影响到井孔模式波的特征.由图4可知,随着频率的增大,动态渗透率归一化幅值变小,表示地层慢纵波带走的能量减少,从而斯通利波与伪瑞利波的衰减减小.而均匀介质中并未考虑动态渗透率的影响,因此计算得到的模式波有较大的衰减.

声波全波列波形能直观反应地层纵波、横波、斯通利波等波形的幅度与到时.本节进一步利用实轴积分法,计算了非均匀含黏性流体孔隙介质中的井孔全波波形,并与均匀介质中的波形进行对比.声源采用余弦包络脉冲形式,对应的函数表达式为

式中,f0为声源中心频率;T为声源脉冲宽度.

图7(a)显示了常规单极子测井的井孔全波波形,其中f0为12 kHz,T为0.3 ms,源距为4 m.图7(b)为8 道阵列全波波形,第一道波形的源距为3.048 m,相邻道之间的间距为0.1524 m.从图7(a)可以看出,全波波形按先后顺序依次为快纵波、快横波、伪瑞利波和斯通利波.由于慢横波的散射作用,井孔接收到的快纵波的到达时间发生了明显的滞后,声压幅度也减小,介质非均匀程度越大,快纵波到时滞后越多,衰减越明显,这与图3的计算结果相符.同时能看出,非均匀介质中的伪瑞利波与斯通利波的幅度与均匀介质中相比,明显变大,而非均匀程度对模式波的影响不大,这与图6展示的结果相符.图7(b)中不同源距下波形的变化与图7(a)中单道波形的变化一致.

图6 模式波的相速度与衰减 (a) 斯通利波;(b) 一阶伪瑞利波;(c) 二阶伪瑞利波Fig.6.Phase velocity and attenuation of mode waves: (a) Stoneley wave;(b) the first order pseudo Rayleigh wave;(c) the second order pseudo Rayleigh wave.

图7 全波波形 (a) 源距为4 m 时的单道波形;(b) 8 道阵列全波波形Fig.7.Full wave waveform: (a) Single channel waveform at source distance of 4 m;(b) 8-channel array full wave waveform.

以上结果表明,在含黏性流体孔隙介质中,非均匀孔隙尺度上慢横波的散射作用将对快纵波和慢纵波产生影响,进一步导致井孔模式波传播特征的改变.快纵波受到孔隙介质非均匀程度的影响较大,而模式波受到非均匀程度的影响较小.均匀介质假设下的井孔声波中,快纵波的幅值偏大,而伪瑞利波及斯通利波的幅值偏小,这是该假设在实际应用与解释中存在很大误差的原因之一.对孔隙介质非均质性的考虑完善了Biot 理论的不足,相应的理论结果对实际孔隙地层中井孔波形的分析具有更好的指导意义.

为了更全面地认识非均匀饱含黏性流体孔隙介质井孔声波的特征,接着分析了全波波形随孔隙流体黏度的变化.一般而言,地层温度随着深度的增加而变高,这使孔隙地层中流体的黏度也随之变化.表2 给出了水在三种不同温度下的黏度.假设孔隙介质非均匀程度Δ=0.05 .图8(a)和图8(b)分别显示了低频斯通利波测井下(f0=2 kHz,T=1.5 ms)及常规单极子声波测井下(f0=12 kHz,T=0.3 ms)的井孔波形随流体黏度的变化.从图8(a)和图8(b)可以看出,低频下,斯通利波的幅度随着黏度的减小而减小,而高频下,伪瑞利波及斯通利波的幅度随着黏度的减小而增大,快纵波随黏度的变化不明显.这表明,孔隙介质中流体的黏度对井孔接收到的模式波有很大的影响,声源的中心频率不同,黏度对模式波的影响不同.斯通利波在孔隙介质渗透率反演中具有重要的作用,图9 进一步给出了斯通利波相速度与衰减随流体黏度的变化.其中,低频时,流体黏度越大,斯通利波的速度越小,高频时,流体黏度对速度基本没有影响;而对于衰减,低频范围流体黏度的变化对衰减的影响很大,且黏度越大,衰减越小,高频时,黏度的影响减小,但黏度越大,衰减越大.因此,在应用井孔斯通利波进行渗透率等地层信息的提取中,有必要考虑和消除孔隙流体黏度的影响,以减少参数反演的误差.