鲜 强,冯许魁,刘永雷,吕 东,安海亭,蔡志东,王祖君,苑恒超

[1.中国石油 东方地球物理公司研究院,新疆 库尔勒 841000; 2. 中国石油 东方地球物理公司 新兴物探开发处,河北 涿州 072750]

经过十几年的油气勘探及前人研究,大家逐渐认识到塔里木盆地塔中地区鹰山组碳酸盐岩储层为潜山岩溶型油气藏。但在油气检测方面一直没能得到让人信服的结果,主要有如下几个原因：① 塔中地区地表被沙漠覆盖,地震资料信噪比低,难以满足叠前流体预测的需要;② 塔中地区储层埋藏深(>6 000 m),储集空间主要为缝、洞等溶蚀体,储层非均质性强;③ 钻遇缝洞型储层,容易发生放空,难以获得缝洞型储层的纵、横波曲线,利用纵横波曲线进行AVO正演,难以确定缝洞型储层中不同流体的AVO响应。

针对地震资料信噪比低的问题,采取小面元、高密度地震采集技术及相应的处理技术能很好的解决地震数据的信噪比低和储层成像精度差等问题。当钻遇缝洞型储层容易发生放空与漏失,孔隙度较大(>30%),横波预测困难。针对上述问题,前人做了不少研究:杨辉廷[1]、张广智[2]等人论证了修正型的Xu-white模型[3],可以用来预测缝洞型储层横波速度,Xu-white模型对于孔隙型的砂泥岩效果较好,对于碳酸盐岩缝洞型储层横波预测效果一般。杨敏[4]尝试利用DEM模型[5]进行碳酸盐岩横波预测,预测精度较高。镇晶晶对比了DEM模型与Self-consistent模型[6],认为二者假设条件和局限性相似,DEM模型与包含物加入的顺序有关,Self-consistent模型只与最后物质形式有关[7]。本文选择了Self-consistent模型进行横波预测,受孔隙大小影响较小,通过调整孔隙结构可以用于孔隙度大于临界孔隙的岩石物理建模。

AVO正演以及叠前同时反演的理论基础是Zoeppritz方程及其在不同假设条件下的简化。1985年Shuey对Zoeppritz方程进行了重新推导与简化,将地层泊松比变化与入射角的关系引入简化公式,提出了Shuey近似[8]方程。Shuey近似方程较为直观地反映了纵波反射系数与入射角及界面两侧介质的弹性参数间的关系。Connolly[9]提出了弹性波阻抗反演(EI)方法;Whitecombe等[10]提出了扩充弹性阻抗反演(EEI)方法,并应用于流体和岩性预测;Li Yongyi等人[11]利用AVO属性,对加拿大西部盆地的碳酸盐岩进行了叠前油气检测;李宗杰[12]利用正演模拟实验以及叠前弹性阻抗反演技术成功的预测碳酸盐岩储层的流体性质;唐金良[13]利用叠前反演技术,针对塔河油田碳酸盐岩进行了叠前油气检测;张远银[14]利用AVO反演预测了塔中地区碳酸盐岩含油气性。

针对上述难题,在前人研究的基础,尝试在塔中地区开展基于岩石物理建模的叠前流体检测。研究对碳酸盐岩岩石物理建模方法演进行了探讨,系统的论证了不同类型储层的AVO响应与叠前流体检测的可行性,最终开展了叠前流体检测实验,并对检测结果进行了详细分析与说明。

1 岩石物理建模

岩石物理是沟通储层与地震响应的桥梁。岩石物理建模是建立起储层弹性参数与地震响应之间的关系,通过岩性和流体的替代分析,研究岩性、物性及流体变化对地震响应的影响,来定量地描述油气与地震响应之间的对应关系[13],从而进行反演的可行性论证及流体敏感参数的优选。岩石物理建模主要包括两个方面：横波预测、流体替代。

1.1 横波预测

横波预测的方法通常是利用经验公式进行预测。Pickett[15],Castagna[16]等许多学者给出了不同的纵波-横波速度的关系,上述学者研究基于原生孔隙储层的纵横波关系。而本区碳酸盐岩主要以溶洞、裂缝等次生孔隙为主,利用经验公式进行横波预测存在较大误差。针对本区缝洞型储层对横波预测方法进行了适当地修正。一方面是将尝试用Self-consistent模型替换Xu-white模型,Self-consistent模型可以预测孔隙度大于30%的储层的横波速度,这对本区洞穴型储层横波预测尤为有利。另一方面是将干岩石骨架由泥岩和砂岩换成泥岩和灰岩,通过Viogt-Reuss-Hill平均计算混合矿物的体积模量和剪切模量[17]。第三,利用Self-consistent模型,通过调节矿物的横纵比来实现溶洞和裂缝对碳酸盐岩弹性模量的影响,从而计算出干岩石的弹性模量。第四,利用Wood公式[18],计算在一定压力、温度、矿化度条件下,混合流体的体积模量及密度。最后,利用Gassmann方程[19]计算饱含流体岩石的体积模量和对应的纵横波速度。

利用上述方法,计算出的横波速度受流体影响小,预测结果与实测曲线吻合度较高。而纵波受流体及裂缝影响较大,预测的速度与实测的总体趋势背景一致,局部有误差。这种结果表明改进的方法在本区进行横波预测是有效的。

1.2 岩石物性分析及流体替代

流体替代就是替代孔隙中的流体,分析饱含不同流体的对岩石弹性模量的变化引起的纵、横波速度的变化,最终定量分析流体变化引起的地震响应的变化。流体替代的核心是Gassmann方程。在前面岩石物理建模的基础上,确定岩石骨架及流体成分,保持岩石骨架不变,替换孔隙中的流体。

研究区内储层埋深6 500 m左右,油气藏的类型为凝析气藏,气油比约2 000 ∶1。目的层地层压力为60 MPa,地层温度大约为140 ℃,地层水的矿化度80 000×10-6。本区储层为深埋阶段的碳酸盐岩储层,储层孔隙主要是次生溶孔和裂缝,40%的井在钻探过程中发生了放空,孔隙度变化很大(2%～30%,局部可以达到100%)。储层岩性以灰岩为主,图1中非储层密度为2.71 g/cm3,阻抗为17 300 g·cm3·m·s-1。当岩石发生溶蚀,洞穴中充注油气后密度及阻抗会迅速下降,图2中密度下降到2.5 g/cm3,阻抗下降到12 000 g·cm3·m·s-1左右。纵波速度(vp)由6 359 m/s下降到4 700 m/s,下降幅度为25%。横波速度(vs)由3 200 m/s下降到2 700 m/s,下降幅度为15%。vp/vs由1.91下降到1.55,下降幅度为21%。

图1 塔中地区密度与阻抗交汇图Fig.1 Cross-plot between density and P-impedance in Tazhong area

图2 塔中地区纵横波速度比与阻抗交汇图Fig.2 Cross-plot between vp/vs and P-impedance in Tazhong area

图3 塔中地区某气井流体替代图 Fig.3 Fluid substitution for gas well in Tazhong area(图中粉色曲线为混合流体的替代结果,黑色为实测曲线,蓝色曲线为饱含水的预测结果。)

在较大的埋深条件下,当孔隙中充满水时,图3中纵波速度变化不大;当孔隙中含气时,纵波速度明显下降。不论孔隙中含有何种流体,横波速度整体变化较小。本区流体主要为凝析气,预测流体结果与实测接近。单纯的使用纵波阻抗难以区分储层中的流体,通过前面储层与非储层的交汇以及流体替换试验可知,纵波速度与横波速度的比则可以有效区分油气与水,因此可以选取纵波与横波的速度比作为敏感参数进行流体识别。

2 AVO理论及正演

AVO正演从实测纵横波及密度曲线出发,利用Zoeppritz方程计算出一个叠前道集。分析这个叠前道集振幅随偏移距变化的关系可以确定AVO响应特征。基于流体替代的AVO正演,可以定量分析流体变化对AVO响应的影响,最终实现反演的可行性分析。

本文在工区内选择了油井与水井进行AVO正演,正演结果显示图4中油井的AVO响应是振幅随偏移距增大而增大,而图5中水井的AVO响应是振幅随偏移增大基本不变或变化不大。不同类型井的截距和梯度是有差异的,可以利用AVO特征或其它叠前道集属性在本区开展定性的流体识别研究。

3 叠前同时反演

3.1 道集能量补偿基于观测系统属性

地震数据受观测系统和偏移成像方法的影响,导致CDP道集振幅在不同偏移距位置分布不均匀,中偏移距振幅能量强,而小偏移距与大偏移距处能量弱。除了常规的处理流程,在偏移之后还需要进行道集剩余能量补偿。常规的做法是：对道集内振幅进行最大与最小值统计,在道集内部进行能量归一化,这种方法一定程度上能消除不同偏移距上的能量异常。

图4 塔中地区油井AVO正演Fig.4 AVO forward modeling for oil well in Tazhong area

图5 塔中地区水井AVO正演Fig.5 AVO forward modeling for water well in Tazhong area

为了使振幅随偏移距的变化的关系符合正演模型,本文选择了基于AVO正演的道集能量补偿方法。具体做法是对非储层段如本区奥陶系灰岩顶面进行多口井AVO正演,然后对AVO趋势进行统计,拟合一个振幅随偏移距变化的关系,最终将这个关系应用到工区道集。使得不同类型的井旁道集的振幅变化与井的AVO正演道集振幅变化趋势一致。校正前,道集振幅在中偏移距能量强,在近、远偏移距振幅能量弱,见图6。校正后,道集振幅随偏移距增加而增大,与井旁(图4)正演道集的振幅变化趋势一致(图7)。初步实现了利用井曲线进行AVO正演对地震道集数据振幅变化的质控。

3.2 叠前同时反演

地震叠前反演的理论与AVO正演及属性分析理论相似,都是Zoeppritz方程及其各种类型的简化。本文在前人研究的基础上,以新的岩石物理模型以及高密度地震资料为基础,对塔中地区流体性质进行进一步的探索。

图6 塔中地区常规油井CRP道集Fig.6 CRP gather for conventional wells in Tazhong area

图7 塔中地区油井CRP道集校正后Fig.7 CRP gather energy compensation for oil wells in Tazhong area

3.2.1 道集部分叠加

AVO属性分析能够相对快速地、半定量地进行流体检测分析,而叠前同时反演则可以准确地、定量地获得能够反映流体的弹性参数。通过叠前同时反演可以得到纵波横波以及密度等弹性参数,这些弹性参数可以较好的反映由于岩性及物性变化引起的振幅随偏移距的变化,最终预测储层中的流体。

由于本区埋藏较深,AVO响应都在中远偏移距出现异常。在偏移距3 000 m左右出现振幅的变化,在远偏移距的位置容易出现动较拉伸等现象,在近偏移距的位置有多次波及其他类型波的干扰。在选择偏移距时,适当的去掉远与近偏移距数据,去除干扰的基础上保留AVO现象。

针对部分数据体的叠加方法本文主要考虑两点,单个限偏移距叠加数据体的信噪比和反演结果的稳定性。部分叠加数据体越多,反演结果越稳定。在覆盖次数一定的基础上,部分叠加数据体越多,单个数据体的覆盖次数会降低,从而导致单个数据信噪比会下降,最终会影响反演结果的稳定性。本区选取有效偏移距为1 000～6 000 m,按照偏移距增量平均将其分为5～7个部分,然后将其叠加生成5～7个部分叠加数据体。为了求解纵、横波速度以及密度等参数,至少需要3个部分叠加数据体。利用5个部分叠加数据体进行反演能够得到稳定的横波数据,利用7个部分叠加数据体进行反演能够得到稳定的横波与密度数据。将本区道集叠加成7个和5个部分叠加数据体,7个部分叠加数据体整体信噪比较低,因此使用5个部分叠加数据体进行反演。5个部分叠加剖面整体信噪比较高,对储集体反映能力较好。

3.2.2 多子波标定

利用理论子波制作合成记录,与全叠加数据体进行比较,计算其相关系数,建立初始的时深关系。然后,对每个偏移距叠加数据体分别提取井旁道子波,修正初始的时深关系,子波提取与合成记录标定迭代进行,当标定的相关系数达到最大时,确定最终的时深关系,对应于不同偏移距叠加数据体同一口井的时深关系是唯一的[20]。

3.2.3 低频模型的建立

由于地震数据的频带是有限的,缺失高频和低频信息。基于约束稀疏脉冲算法的叠前反演能够提高地震数据的带宽,需要利用工区内测井低频信息建立趋势模型进行阻抗和密度的趋势约束,保证反演结果的稳定性。同时测井资料提供了较高的垂向分辨率,较多的保留了测井的高频信息。

3.2.4 叠前同时反演及结果分析

结合偏移速度体,对不同偏移距数据体进行叠前同时反演,通过反演可以得到反映流体的弹性参数。本区储层在充填泥质、油气及水,其地震响应都是串珠状响应(图8a),在叠后反演剖面上都为低阻抗的响应(图8b)。常规的叠后反演预测储层能力较好,而反映油气的能力有限。ZG8-H1井为水井,在叠后反演上为储层响应,在vp/vs剖面上数值下降较小,由1.8下降到1.7(图8c)。ZG8-4井为低效井(油气累计产量低),在叠后反演剖面上同为储层响应,在vp/vs剖面上数值下降不明显,由1.8下降到1.65。而ZG8井为高效井(油气累计产量高),含油气后储层的vp/vs值有明显下降,由1.8下降到1.6以下。叠前反演得到的纵横波速度比能够有效的识别储层中的流体。

图8 塔中地区叠前油气预测剖面(剖面位置见图10)Fig.8 Pre-stack hydrocarbon prediction in Tazhong area(seeFig.10 for profile Cocation)a.叠后地震剖面;b.波阻抗剖面;c.纵、横波速度比剖面

图9 塔中地区ZG8井区鹰山组储层平面预测Fig.9 Reservoir prediction of the Yingshan Fm in ZG8 well-block of Tazhong area

图10 塔中地区ZG8井区鹰山组叠前油气检测平面图Fig.10 Hydrocarbon prediction plan for Yingshan Fm in ZG8 well-block of Tazhong area

塔中地区鹰山组储层为早期潜山岩溶和后期断控岩溶共同作用下的产物。塔中地区在鹰山组沉积后,经历了短暂的抬升,在混合水差异溶蚀下[21],在工区西南部形成了大面积不规则的储层,红色为洞穴型储层(图9)黄色为孔隙型储层。在工区的北部储集体以洞穴型为主,为后期断裂改造后形成断溶体[22-23],在平面上呈点状,局部串珠呈线状间接反映断裂发育方向。通过叠前流体检测后,结果显示工区南部潜山型储层区含油较差(图10)。钻井资料揭示潜山型储层在暴露期泥质充填,储集空间难以保存。工区北部为断控岩溶发育区。由于没有暴露剥蚀,同时储集空间得到了较好的保存,在多期断裂活动的改造下,储集空间得到进一步的扩大,为油气充注提供条件,为有利储层发育区。对本工区的18口井进行了预测,15口井的预测结果与钻探结果一致,吻合率达到83.3%,利用此结果提供了多个钻探目标。

4 结论

针对缝洞型碳酸盐岩储层的复杂孔隙结构,本文使用Self-consistent模型能够较为准确的预测缝洞型储层的横波速度。通过岩石物理分析及流体替代研究可以优选出反映流体变化的敏感参数,利用纵横波速度比可以有效的进行流体预测。通过AVO正演可知,碳酸盐岩缝洞型储层饱含油气时,AVO的响应为：振幅随偏移距增大而增大;储层含水时,AVO响应为：振幅随偏移距增大保持不变。对常规CRP道集进行振幅能量补偿,消除了观测系统造成的道集能量不均。在AVO正演的基础上,对实际的地震道集进行基于AVO正演的振幅能量补偿。利用多个分偏移距叠加的地震数据进行叠前同时反演,反演得到的弹性参数vp/vs,能够有效的反映研究区内流体。

致谢：感谢领导、专家以及同事的建议、指导与帮助。同时感谢中国石油塔里木油田公司对研究提供的资料支持