裂缝型储层五维地震有效压力参数预测

张洪学,印兴耀,李 坤,姜 曼

(1.中国石油大学(华东)地球科学与技术学院,山东青岛266580;2海洋国家实验室海洋矿产资源评价与探测技术功能实验室,山东青岛266071)

宽方位五维地震数据是指通过地震宽(全)方位观测系统采集,并经过炮检距向量片(OVT)等技术处理获得的包含五维(即空间三维+炮检距+方位角)信息的叠前地震道集,宽方位五维地震数据逐渐成为复杂油气勘探及高精度地震解释的重要基础数据[1-3]。针对储层的有效压力预测在油气勘探开发中降低钻井成本、保障钻井安全、保护油气层等方面发挥着关键作用[4-6]。国内外学者针对压力的岩石物理研究已经开展了大量基础与实践研究,主要包括纵横波速度、岩石模量或速度比与有效压力之间的经验关系式[7-11],有效压力与软孔隙纵横比变化的关系式[12]或气饱和超压页岩岩石物理模型[13]等,但在实际裂缝型储层有效压力预测方面仍面临挑战,瓶颈在于现有有效压力预测方法多以岩石物理模型为主,或通过有效压力与上覆地层压力与孔隙压力的关系间接计算,缺少针对综合五维地震资料进行裂缝型储层有效压力稳定预测的研究。

裂缝型储层五维地震有效压力预测的关键是寻找对有效压力敏感的参数及建立相应表征反射系数方程。学者们基于岩石物理实验建立了岩石模量与有效压力之间的线性关系[11]或对数关系[14],在一定程度上可以通过岩石的模量反映其受到的有效压力情况。借助统计岩石物理模型将岩石物理学关系和地质统计学相结合,可以建立起岩石弹性参数和地震响应之间的联系[15],进而指导有效压力的预测。但统计岩石物理模型是通过数据拟合和随机误差获得的[16],具有较强的随机性。SMITH[17]建立了有效压力与孔隙度间的指数关系式,为有效压力参数的预测提供了良好的途径。目前广泛使用间接方法预测有效压力,如利用纵横波速度进行孔隙压力预测,利用地层埋藏深度求取上覆地层压力,然后再进行有效压力的求解[18],该方法容易受其它干扰因素(岩性和断裂等)影响,因此该方法需要在预测结果的基础上考虑校正因子[19]。此外,从地震岩石物理理论出发,构建压力敏感参数的刚度矩阵,基于Born近似和稳相法求解,可推导与待反演参数有关的反射系数近似方程,再通过概率化方法可进行反演预测[20],目前针对有效压力的研究已从各向同性储层[21]发展到各向异性储层[22]。本文的特点是在于联合有效压力关系式[17]、Nur临界孔隙度模型[23]以及各向异性Gassman方程[24]构建新的刚度矩阵方程,优化参数组成,构建OA介质有效压力敏感参数表征的方位反射系数方程。

裂缝型储层的特殊性以及地下埋藏条件的复杂性给裂缝型储层五维地震反演带来了挑战,多裂缝参数同时反演更是导致储层参数精度下降。针对裂缝型储层,常用的直接反演方法是弹性阻抗反演方法[25-39]和基于贝叶斯框架的AVAZ反演方法[40-42],除单纯利用地震纵波信号反演外,学者们也进行了多波联合反演[42]及非线性反演方法研究[43-44]。时间域地震反演的抗噪性能相比频率域有所增强,但是时间域地震反演的分辨率不及频率域地震反演。因此在裂缝型储层压力预测中结合时频域地震反演的优势[45],充分利用五维地震资料中多波(纵、横)多特征(振幅、频率及相位)信息,采取不同的反演策略,可提高有效压力预测的可靠性。

根据各向异性Gassman方程及Schoenberg线性滑动模型,联合有效压力关系式[17]和Nur临界孔隙度模型[23]构建OA介质新的岩石物理刚度矩阵方程,推导出利用有效压力参数及流体模量和裂缝参数直接表征的新的裂缝型储层五维地震反射系数方程,并分析方程的适用条件,从理论上阐明了采用该方程开展裂缝型储层有效压力参数预测及流体识别的可行性。在此基础上,创新了有效压力参数预测方法,采用五维地震时频域分步反演的方法,基于分步反演及时频域联合反演策略,在贝叶斯理论框架下构建五维地震时频域分步反演目标泛函,基于五维地震频率和振幅信息及岩石、测井和地质资料,实现了裂缝型储层有效压力参数及流体参数和裂缝法向及切向裂缝弱度参数的直接预测。

1 方法原理

1.1 OA介质五维地震有效压力参数刚度矩阵构建及分析

各向异性Gassmann方程提供了描述裂缝型储层等效参数的理论基础,GUREVICH[46]结合线性滑动理论将各向异性Gassmann方程推导为:

(1)

(2)

其中,

(3a)

(3b)

(3c)

(3d)

结合Nur临界孔隙度模型[23],利用临界孔隙度参数φc可以将干岩石等效体积模量表示为基质矿物模量与孔隙度的形式,即

(4)

SMITH[17]提出了一种简易的指数形式的用孔隙度表示的有效压力参数函数,即

φ≈φ0e-αPeff

(5)

式中:φ0为初始孔隙度;α是有效压力系数,与岩石压实程度有关。令α=0.096MPa-1[49],Peff为有效压力,其矢量形式为有效应力,指的是地层岩石骨架和岩石基质所受的应力[50]。由公式(5)可知,储层有效压力与储层孔隙度存在指数关系,该经验公式适用于致密储层及裂缝发育的页岩储层。为便于后续反演研究,构造有效压力敏感参数Pe=e-αPeff[21],有效压力敏感参数Pe仅与有效压力系数α和有效压力Peff相关;有效压力敏感参数Pe可反映储层有效压力的大小。进而公式(5)可表示为孔隙度与有效压力的线性关系φ≈φ0Pe;有效压力敏感参数也可表示为Pe≈φ/φ0,由储层孔隙度求得。另外,Peff=Pz-Pp,其中,Pz为上覆地层压力或垂直压力。其矢量形式为垂直应力,是由上覆地层的岩石骨架,岩石基质及裂缝、孔隙中流体所受的重力对地层作用产生的;Pp为孔隙压力,是由裂缝、孔隙中的流体对地层作用产生的[49]。

考虑到流体体积模量通常远远小于基质固体颗粒的体积模量,即Kf≤Km[48],结合公式(4)和公式(5)可以得出:

(6)

(7)

将公式(3)代入公式(1),结合公式(6)和公式(7),可获得裂缝介质纵波模量、有效压力参数、两组裂缝的干裂缝弱度参数及流体体积模量表征的饱和流体OA介质的刚度矩阵,即:

(8)

(9a)

(9b)

(9c)

(9d)

(9e)

(9f)

(9g)

(9h)

(9i)

设置模型验证上述模型近似刚度参数的准确性。模型储层孔隙中饱含气和水,含水饱和度Sw分别为30%和70%,矿物成分由石英和黏土组成,比例为1∶1,使用Wood公式计算岩石的混合流体模量,使用VRH方程计算矿物平均模量[48],裂缝弱度参数分别为δN1=0.02和δN1=0.04及δN2=0.06和δN2=0.12,进而对弱各向异性近似刚度方程(9)与精确刚度表征方程(2)进行对比,如图1～图8是不同含气饱和度及裂缝弱度情况下,裂缝岩石弱各向异性近似(粉色)与精确(蓝色)刚度参数随孔隙度变化的对比,孔隙度变化范围为0～0.2,比较图1至图8可知,在含水饱和度不变的情况下,法向裂缝弱度越小,则近似刚度参数越准确,故推导的近似公式更适用于弱各向异性介质;图1至图8对比可知,在裂缝弱度不变的情况下,含水饱和度越低,则近似刚度参数越准确,故推导的近似公式更适用于低体积模量裂缝型储层。

图1 裂缝模型精确刚度矩阵参数(蓝色)与近似刚度矩阵参数(粉色)对比(Sw=0.70,δN1=0.02,δN2=0.06)

图2 裂缝模型精确刚度矩阵参数(蓝色)与近似刚度矩阵参数(粉色)对比(Sw=0.70,δN1=0.02,δN2=0.12)

图3 裂缝模型精确刚度矩阵参数(蓝色)与近似刚度矩阵参数(粉色)对比(Sw=0.30,δN1=0.02,δN2=0.06)

图4 裂缝模型精确刚度矩阵参数(蓝色)与近似刚度矩阵参数(粉色)对比(Sw=0.30,δN1=0.02,δN2=0.12)

图5 裂缝模型精确刚度矩阵参数(蓝色)与近似刚度矩阵参数(粉色)对比(Sw=0.70,δN1=0.04,δN2=0.06)

图6 裂缝模型精确刚度矩阵参数(蓝色)与近似刚度矩阵参数(粉色)对比(Sw=0.70,δN1=0.04,δN2=0.12)

图7 裂缝模型精确刚度矩阵参数(蓝色)与近似刚度矩阵参数(粉色)对比(Sw=0.30,δN1=0.04,δN2=0.06)

图8 裂缝模型精确刚度矩阵参数(蓝色)与近似刚度矩阵参数(粉色)对比(Sw=0.30,δN1=0.04,δN2=0.12)

1.2 裂缝型储层五维地震有效压力参数地震反射方程推导

反射系数与散射函数之间的关系可简单地表示为[51-52]:

(10)

其中,

(11)

(12a)

(12b)

式中下标符号m和n与i,j,k及l的关系可表示为:

m=iδij+(9-i-j)(1-δij)

(13a)

n=kδkl+(9-k-l)(1-δkl)

(13b)

(i,j,k,l=1,2,3)

式中:δij和δkl均为Kronecker函数。

(14)

式中:

式中:符号Δ表示裂缝型储层参数的扰动量。

将方程(14)扰度刚度矩阵代入公式(10),可以建立有效压力敏感参数、流体体积模量、介质剪切模量、密度、水平裂缝与垂直裂缝法向弱度、切向弱度等参数直接表征的反射系数方程,即:

(15)

式中:

为研究有效压力敏感参数、流体体积模量、介质剪切模量、密度、水平裂缝与垂直裂缝法向弱度、切向弱度参数变化对公式(15)的贡献度差异,逐一分析每个参数对反射系数产生的影响。图9分别展示了在参数统一变化范围(-0.30∶0.15∶0.30)内引起的地震反射系数变化特征,整体分析可以看出,各向同性参数的贡献度要远大于各向异性部分,重点分析有效压力敏感参数项,即图9c所示内容,随着入射角及参数扰动范围的变化,有效压力敏感参数对OA介质反射系数的贡献较明显,相较裂缝参数更容易实现稳定反演,其它各向同性介质参数如流体体积模量和介质剪切模量对反射系数影响同样较明显,图9d中密度参数对反射系数引起的变化较小,但仍比裂缝参数的明显。由图9e至图9h分析可知,裂缝参数对反射系数整体贡献度较低,裂缝参数之间贡献度差异不明显,随着地震波入射角的增大,裂缝参数对反射系数的贡献度逐渐增加,因此,相较小角度地震数据,裂缝参数对大角度方位地震数据更敏感,利用大入射角的宽方位五维地震资料更容易实现上述参数的稳定反演。

图9 参数统一变化范围(-0.30∶0.15∶0.30)内引起的地震反射系数变化特征a 流体体积模量;b 剪切模量;c 有效压力敏感参数;d 密度;e 水平裂缝法向弱度;f 水平裂缝切向弱度;g 垂直裂缝法向弱度;h 垂直裂缝切向弱度

1.3 裂缝型储层五维地震有效压力参数反演方法

对于有效压力敏感参数地震反射方程(15),将其中的非裂缝待反演参数的相对差值用自然函数的差值Δln(·)近似代替,即ΔKf/Kf0≈Δ(lnKf),Δμ/μ0≈Δ(lnμ),ΔPe/Pe0≈Δ(lnPe)和Δρ/ρ0≈Δ(lnρ),进而消除了反射系数方程中的分式项。针对弱各向异性裂缝型储层,各向异性特征及非均质特征较弱,储层地球物理参数变化连续,因此可将待反演参数的差值形式表示为:Δ(lnKf)≈d(lnKf),Δ(lnμ)≈d(lnμ),Δ(lnPe)≈d(lnPe)和Δ(lnρ)≈d(lnρ),其中裂缝弱度参数同样表示为ΔδN1≈dδN1,ΔδT1≈dδT1,ΔδN2≈dδN2和ΔδT2≈dδT2,上述参数表征为利用五维地震资料反演奠定了基础。

将方程(15)与时间域和频率域地震子波分别褶积,得到时间域和频率域的正演方程,可表示为:

(16)

(17)

式中:s(t)和S′(ω)分别是时间域和频率域地震记录;W(t)和W′(ω)分别为时间域和频率域子波;E=e-iωτk代表傅里叶算子,令l=1∶T,ωl表示第l个有效频率成分,D是差分算子。

将公式(16)和公式(17)写成矩阵的形式:

(18)

[Giso(t)]MNK×4K=

(19)

其中,

(20)

其中,

(21)

其中,

(22)

其中,

借助待反演参数模型协方差矩阵Ccov的特征值和特征向量,对正演算子矩阵及模型参数矩阵进行处理,以避免因参数尺度或者相关性的原因干扰反演结果,首先奇异值分解可得矩阵的特征向量U和特征值正定矩阵Σ,即Ccov=UΣUT,相关结果如下:

(23)

(24)

因此,地震记录可以写为:

(25)

式中:kron(·)表示Kronecker乘积算子,IL为L阶单位矩阵。构建协方差矩阵的方式有多种,本文是基于测井数据统计所得,L为井曲线采样点数。

(26)

(27)

2 模型试算及实际资料处理

2.1 模型试算

利用实际裂缝工区中的钻井资料合成方位地震数据,验证五维地震时频域分步反演有效压力敏感参数、流体体积模量、介质剪切模量、密度、水平裂缝与垂直裂缝法向弱度、切向弱度参数的可行性,裂缝弱度参数通过岩石物理建模利用井上钻井资料估算得到。首先,基于五维地震褶积模型合成8个方位:22.5°,45°,67.5°,90°,112.5°,135°,157.5°和180°的地震道集。在合成地震道集中添加部分随机噪声,生成不含噪声及信噪比分别为5∶1和2∶1的合成方位道集,分别如图10至图12 所示。基于上述合成的8个方位地震道集,利用五维地震时频域分步反演方法进行方位反演测试。反演结果见图13至图18,图中绿色、蓝色和红色曲线分别表示待反演参数的初始模型、模型参数真实值及模型参数的反演结果。分析不含噪声或含有较强噪声地震记录反演结果可知,有效压力敏感参数、流体体积模量、介质剪切模量、密度、水平裂缝与垂直裂缝法向弱度、切向弱度参数反演结果与实际测井数据趋势吻合较好,且数值之间误差较小,即使在噪声较高的情况下,同样可实现多参数的直接反演。因此该方法具有较高的准确度和抗噪性,验证了五维地震时频域分步反演方法直接预测上述8个参数的可行性。

图10 无噪声的不同方位地震记录a 22.5°;b 45°;c 67.5°;d 90°;e 112.5°;f 135°;g 157.5°;h 180°

图11 信噪比为5∶1的不同方位地震记录a 22.5°;b 45°;c 67.5°;d 90°;e 112.5°;f 135°;g 157.5°;h 180°

图12 信噪比为2∶1的不同方位地震记录a 22.5°;b 45°;c 67.5°;d 90°;e 112.5°;f 135°;g 157.5°;h 180°

图13 无噪声地震记录流体体积模量、介质剪切模量、有效压力参数及密度的反演结果

图14 无噪声地震记录水平裂缝法向弱度(δN1)、切向弱度(δN2)与垂直裂缝法向弱度(δT1)、切向弱度(δT2)的反演结果

图15 信噪比为5∶1的流体体积模量、介质剪切模量、有效压力参数及密度的反演结果

图16 信噪比为5∶1的水平裂缝法向弱度(δN1)、切向弱度(δN2)与垂直裂缝法向弱度(δT1)、切向弱度(δT2)的反演结果

图17 信噪比为2∶1的流体体积模量、介质剪切模量、有效压力参数及密度的反演结果

图18 信噪比为2∶1的水平裂缝法向弱度(δN1)、切向弱度(δN2)与垂直裂缝法向弱度(δT1)、切向弱度(δT2)的反演结果

2.2 实际资料处理

利用我国某地区裂缝型储层的五维地震资料进行处理,地质资料显示该地区储层横向连续性较好,且发育大量高角度倾角裂缝,将该储层裂缝介质等效为正交各向异性介质。选取方位角分别为φ1=45°,φ2=75°,φ3=135°和φ4=165°的4个方位地震道集,每个方位地震道集均经过了井控道集优化处理。图19 和图20是这4个方位的近、中、远部分叠加道集地震二维剖面,其部分叠加的入射角分别为10°(叠前角度为4°～16°)、22°(叠前角度为16°～28°)及34°(叠前角度为28°～40°)。图21至图24是利用井上流体体积模量、介质剪切模量、有效压力敏感参数、密度、水平裂缝与垂直裂缝法向弱度、切向弱度参数插值外推建立的初始模型剖面。图25至图28是利用五维地震时频域分步反演方法得到的反演结果剖面。上述参数反演结果整体分辨率较高,细节特征明显,趋势变化符合地震资料特征,与井上数据对比吻合程度较高。由图26b可以看出,有效压力敏感参数反演的连井剖面与测井解释保持一致。因此可以在实际资料中实现有效压力敏感参数的直接反演。为更加明确显示流体体积模量反演预测的准确性,图29为流体体积模量直接反演结果。图中增加了A井和B井的信息对比,其中B井为测试井。需要说明的是,图29反演结果是利用A井数据建立初始模型反演得到的,反演过程中没有加入测试井B的任何信息。分析图中红线标识区域可知,流体体积模量反演剖面中的低值对应该工区流体发育区域,而且反演结果与测试井上流体解释结果吻合,验证了反演结果的准确性。为进一步验证该方法的有效性与可靠性,将井旁道流体体积模量、介质剪切模量、有效压力敏感参数、密度、水平裂缝与垂直裂缝法向弱度、切向弱度反演结果与井上数据进行对比,如图30所示,各参数的反演结果与测井参数吻合程度高,两者之间的误差较小,验证了方法的可靠性。图31是根据有效压力敏感参数计算得到的储层有效压力值,图中色标红色为低值,由此可以分析反演中压力低值区为实际生产提供依据。

图19 方位角为45°和135°时的部分角度叠加道集地震数据剖面a φ1=45°时,小角度地震数据剖面;b φ1=45°时,中角度地震数据剖面;c φ1=45°时,大角度地震数据剖面;d φ2=135°时,小角度地震数据剖面;e φ2=135°时,中角度地震数据剖面;f φ2=135°时,大角度地震数据剖面

图20 方位角为75°和165°时的部分角度叠加地震数据剖面a φ2=75°时,小角度地震数据剖面;b φ2=75°时,中角度地震数据剖面;c φ2=75°时,大角度地震数据剖面;d φ4=165°时,小角度地震数据剖面;e φ4=165°时,中角度地震数据剖面;f φ4=165°时,大角度地震数据剖面

图21 流体体积模量(a)和介质剪切模量(b)的初始模型剖面

图22 密度(a)和有效压力敏感参数(b)的初始模型剖面

图23 水平裂缝法向弱度(a)和垂直裂缝法向弱度(b)的初始模型剖面

图24 水平裂缝切向弱度(a)和垂直裂缝切向弱度(b)的初始模型剖面

图25 五维地震时频域分步反演方法得到的流体体积模量(a)和介质剪切模量(b)反演结果剖面

图26 五维地震时频域分步反演方法得到的密度(a)和有效压力敏感参数(b)反演结果剖面

图27 五维地震时频域分步反演方法得到的水平裂缝法向弱度(a)和垂直裂缝法向弱度(b)反演结果剖面

图28 五维地震时频域分步反演方法得到的水平裂缝切向弱度(a)和垂直裂缝切向弱度(b)反演结果剖面

图29 流体体积模量直接反演结果

图30 原始测井数据(蓝色)与井旁模型参数(红色)反演结果对比a 流体体积模量、剪切模量、有效压力参数及密度反演结果对比;b 水平裂缝法向弱度(δN1)、切向弱度(δN2)与垂直裂缝法向弱度(δT1)、切向弱度(δT2)反演结果对比

图31 储层有效压力计算结果

3 结论

五维地震有效压力参数预测方法旨在从五维地震数据中直接获取裂缝型储层有效压力参数,以正交各向异性介质为研究对象,推导有效压力敏感参数直接表征的新的纵波反射系数弱各向异性近似方程,在此基础上,创新了贝叶斯框架下五维地震时频域分步反演方法。模型测试结果表明,即使是在高噪声(信噪比为2∶1)的情况下,该方法参数反演结果仍与真实模型参数吻合,针对小贡献度参数(裂缝参数)的反演结果平均误差也较小,验证了五维地震时频域分步反演方法的可行性与稳定性。利用经过井控道集优化处理的五维地震资料,实现了裂缝型储层有效压力参数及流体参数和裂缝法向和切向裂缝弱度参数等多参数反演预测。分析可知,反演结果与实际测井数据吻合良好,进一步说明了该方法在裂缝型储层有效压力预测中具有广泛的应用前景。需要说明的是,五维地震有效压力预测应重视地震道集的优化处理,如果面向品质较差的地震资料,该方法反演效果会降低,后续研究可以在该方面做出优化改进。