中江气田沙溪庙组气藏致密砂岩储层测井评价

2020-02-07刘成川黎华继王志章胡华伟李定军

石油物探 2020年1期

刘成川,陈俊,黎华继,王志章,胡华伟,李定军

(1.中国石油化工集团有限公司西南油气分公司勘探开发研究院,四川成都610100;2.中国石油大学(北京)地球科学学院,北京102200)

中江气田沙溪庙组气藏(以下简称中江沙溪庙组气藏)是受断砂配置控制的以构造-岩性气藏为主的远源次生气藏。该气藏以浅水三角洲平原-前缘分流河道砂沉积为主要沉积物。储层具有岩性复杂、厚度薄,致密(孔隙度Ф约为8.66%,渗透率K约为0.21mD(1mD=0.987×10-3um2)),非均质性强的特点。储集空间类型多样、孔隙结构复杂,次生孔隙占据重要地位,孔径极差大,孔渗关系差,与大牛地气田下石盒子组(Ф约为8.6%,K约为0.81mD)和苏里格气田石盒子组(Ф约为10.3%,K约为0.57mD)气藏的情况类似,是典型的低孔、低渗、低含气丰度的“三低”气藏[1]。测井储层参数是测井综合评价的重要依据,气藏的复杂性导致了测井响应特征复杂多样,使得孔隙度、渗透率、饱和度等测井储层参数精度低,造成气、水层误判,从而影响了气田的整体评价与开发。

由于“三低”油气藏与常规储层差异明显,学者们进行了大量的测井精细解释模型研究,如建立孔隙度复相关经验关系式,在渗透率模型中考虑孔隙喉道中粘土的影响,采用不同地质模型和数学方法评价饱和度等。随着电成像、核磁共振等测井新技术的迅速发展,目前主流的研究方法利用神经网络等算法与核磁共振T2全谱数据预测致密砂岩储层渗透率[2]。利用Archie公式研究含水饱和度时,选取准确可靠的岩电参数是研究的关键。学者们对影响岩电参数的因素已经基本达成共识,由于受岩石物性、粒度、孔隙结构等众多因素的影响,选用统一的岩电参数不能满足饱和度评价的精度要求,因此各油气田针对各自的地质条件开展孔隙度、流体矿化度、地应力等参数与重要岩电参数的关系研究,并建立饱和度指数、胶结指数计算模型,力求提高含气(油)饱和度的计算精度[3-5]。

中江气田目前已进入开发阶段,因为仅采用常规测井方法(主要为自然伽马、深浅侧向电阻率、补偿声波、补偿中子以及密度测井等),故不能直接应用基于核磁共振和电成像的测井评价方法,因此需要完善基于常规测井资料的储层参数解释模型,为地质建模提供必要的基础资料。

1 储层“四性”关系

1.1 岩性和物性的关系

储层岩石粒度是决定储层物性、电性的关键因素之一。从图1可以看出,当储层岩石碎屑颗粒细砂含量>75%时,测井响应为高自然伽马、高电阻率特征,岩心分析孔隙度Por<4.0%,岩心分析渗透率Perm<0.1mD;当储层岩石碎屑颗粒粗砂含量>75%时,自然伽马值<55API,电阻率<15Ω·m,声波时差>70μs/ft(1μs/ft≈3.28μs/m),孔隙度>10.00%,渗透率为0.1～1.0mD。总体而言,岩石粒度越粗,储层物性条件越好,粗、中粒砂岩是优质储层发育的基础。

图1 储层岩石粒度-电性-物性关系(GR—自然伽马;RD—深侧向电阻率;RS—浅侧向电阻率;AC—补偿声波;CNL—补偿中子;DEN—补偿密度;M2—储层岩石碎屑颗粒粗砂含量;M3—储层岩石碎屑颗粒中砂含量;M4—储层岩石碎屑颗粒细砂含量;PORC—岩心实验分析孔隙度;PERMC—岩心实验分析渗透率)

由储层岩心分析孔隙度-渗透率交会统计结果(图2)可知,细粒砂岩孔隙度主要为3.00%～9.00%,渗透率主要为0.01～0.20mD;中粒砂岩孔隙度主要为6.00%～15.00%,渗透率主要为0.04～1.00mD,岩石粒度对储层的储渗能力具有显著的影响。

图2 储层岩心分析孔隙度-渗透率交会统计结果

1.2 物性和含气性的关系

为进一步探索储层物性与含气饱和度之间的关系,我们对不同含气饱和度下的纵、横波速度进行了测量,结果如图3、图4、图5所示。在不同孔隙度条件下,纵波时差随着含气饱和度的增大而增大,横波速度则随着含气饱和度的增大而降低,纵横波速度比(vP/vS)随着含气饱和度的增大而降低,说明根据纵横波速度比判定储层含气性的方法适用于中江沙溪庙组气藏[6]。

值得注意的是,当储层物性较低(Φ<7%)时,纵波时差虽然随着含气饱和度的增大而增大,但增大幅度小,纵横波速度比对含气饱和度的变化敏感性较差;而当储层物性较高(Φ>7%)时,纵波时差和纵横波速度比对含气饱和度的变化反映则较为敏感。总体来说,在构造、岩性相同的条件下,中江气田沙溪庙组气藏具有物性越好,含气饱和度越高的特点。

图3 纵波时差与含气饱和度关系

图4 横波速度与含气饱和度关系

图5 纵横波速度比与含气饱和度关系

1.3 物性和含气性对电性的影响

中江沙溪庙组气藏不同物性和含气饱和度的储层电阻率差别大,范围约为15～80Ω·m。储层电阻率受岩石物性和含气饱和度的双重影响,储层电阻率随含气饱和度的增大而增大,气层和干层电阻率均随孔隙度的减小而增大。

孔隙度、胶结指数与电阻率数值模型正演模拟结果如图6所示,由图可见,孔隙度相同时,储层电阻率(Rt)值随胶结指数的增大而增大,孔隙度越低,电阻率对孔隙度、胶结指数越敏感,但当Φ>0.12时(图中红线处),电阻率对孔隙度、胶结指数的敏感性明显降低(图6中颜色越浅代表增大幅度越小,颜色越深代表增大幅度越大)。含水饱合度(Sw)分别为25%,50%,75%,100%的条件下,孔隙度、胶结指数和电阻率的数值模型正演模拟结果如图7所示,可以看出,储层电阻率受孔隙度、胶结指数、含水饱和度共同影响,三者均为储层电阻率的重要影响因素(图7中深绿色表明增大幅度大,淡紫色表明增大幅度小)。相比之下,岩石物性对电阻率的影响更为明显。

图6 孔隙度、胶结指数和电阻率的数值模型正演模拟结果

1.4 孔隙结构特性对含气性的影响

早期研究主要是将储集层的孔隙结构研究结果应用于地层流体运移及成藏规律研究,研究成果未广泛应用于油、气、水测井评价。近些年的研究发现,在低孔低渗致密储层中,孔隙结构直接影响储层产能评价和油、气、水层测井评价的准确性[7-10]。因中江沙溪庙组气藏地质特征复杂,所以孔隙结构对该气藏产能的影响尤为明显。

图7 不同含水饱和度(从左到右依次为25%,50%,75%,100%)条件下孔隙度、胶结指数和电阻率数值模型正演模拟结果

图8为同一层位、同一河道的XX306和XX302井的测井曲线。两口井的测井响应特征差异不大,XX306井岩心实验分析孔隙度约为8.58%,XX302井岩心实验分析孔隙度约为9.1%;XX306井的孔隙结构指数(储层品质因子)约为0.18,XX302井的孔隙结构指数约为0.28。XX306井测试获得天然气产量约为900m3/d,XX302井测试获得天然气产量约为4.30×104m3/d,表明低孔低渗储层的含气性及产能与孔隙结构密切相关。

图8 XX306井(a)和XX302井(b)的测井曲线(SP—自然电位;SH—泥质;C1—碳酸盐含量;POR—孔隙度)

2 储层测井评价方法

2.1 孔隙度及渗透率评价

孔隙度评价是在岩心归位的基础上,将声波时差曲线与岩心孔隙度进行拟合,建立关系式并进行总孔隙度解释,再根据威利公式[11-12]计算泥质孔隙度,最后得到泥质含量校正后的有效孔隙度。本文研究区有效孔隙度模型公式为:

(1)

式中:φe为有效孔隙度;tAC为声波时差;Vsh为泥质含量。

渗透率评价即建立岩心分析的渗透率与孔隙度简单经验关系,如指数或幂函数模型,但中江沙溪庙组气藏的储层孔隙和岩石结构复杂,非均质性强,岩心分析孔隙度和渗透率在交会图上分布分散,若依据传统的单一关系建立渗透率模型,则计算的渗透率不能满足致密储层精细评价的要求。因此按照同一流动单元范围内孔隙结构特征大致相同的特点,将岩心样品分为几类流动单元并分类建立渗透率模型,从而达到提高非均质致密储层渗透率评价精度的目的,该方法称为流动单元法。

流动单元指数(FZI)定义如下:

(2)

式中:Fs为孔隙几何形状指数;Sgv为矿物颗粒比表面积。

流动单元指数将岩石矿物特征、孔喉特征相结合用于判定孔隙结构,可准确描述气藏的非均质特征。将(2)式进行反变换即可得到基于流动单元的渗透率为:

(3)

在双对数坐标系中这种孔渗模型是近似的直线(斜率相同)。因此可近似认为同一直线上的样品具有相似的孔喉特征,为同一类流动单元。根据中江沙溪庙组气藏地质及开发情况,利用该区实验分析资料计算标准化孔隙度指标、流动单元指数、储层品质因子等参数,采用K-Means聚类[13-15]分析方法,将研究区样品分为如图9所示的3类流动单元。Ⅰ类流动单元主要包括中、粗粒砂岩以及少量泥岩,门槛压力最低,中值半径最大,最大进汞饱和度最高,平均渗透率约为2.531mD;Ⅱ类流动单元主要包括中、细砂岩以及部分泥岩和泥质粉砂岩,平均渗透率约为0.353mD;Ⅲ类流动单元渗透性最差,主要为细砂岩,泥质含量较大,平均渗透率约为0.075mD。

图9 研究区样品3类流动单元划分结果

分析各流动单元的电性特征可知,AC、RD和RS的测井响应组合对流动单元敏感性高,因此,采用支持向量机(SVM)对井中取心段地层进行流动单元特征识别,并对未取心段地层进行分类预测,可实现不同流动单元的整体判断,正确率达70%。根据不同流动单元渗透率与孔隙度的关系,建立如下的不同流动单元渗透率模型公式及参数(表1)。

将本文提出的流动单元法应用于中江沙溪庙组气藏渗透率预测。如图10所示中蓝色曲线为采用传统方法预测的渗透率曲线,紫色曲线为采用本文方法预测的渗透率曲线,红色曲线为岩心分析的孔隙度及渗透率,可以看出,采用本文方法预测的渗透率更接近岩心分析渗透率。

2.2 饱和度评价方法

如图11和图12所示,中江沙溪庙组气藏岩心的地层因素与孔隙度能够较好地满足双对数关系,电阻率增大率与饱和度的双对数关系也较为明显,表明Archie公式适用于本研究区。不同河道拟合关系存在差异,说明利用包含固定胶结指数(m)和饱和度指数(n)值的Archie公式计算含气饱和度会导致结果出现偏差,因此,我们尝试引入其它测井参数修正模型表征参数,从而完善含水饱和度模型。