鄂尔多斯盆地北部上古生界致密砂岩气藏含气层位多,勘探潜力大,是油气增储上产的重要区域。目前,在鄂尔多斯盆地东北部的神府地区太原组砂岩段已钻获多口气井,并在部分井见到工业气流,但由于储层普遍致密,且非均质性强,勘探开发难度较大。如何厘清研究区储层致密化过程,是致密砂岩气藏深入勘探亟须解决的问题。通过致密砂岩成岩演化分析,建立孔隙度演化模型,对储层致密化过程进行模拟,是解决这一问题的重要技术手段[1-2]。
众多学者[3-8]对孔隙度演化模型的研究已较为深入,Athy[3]基于对泥岩的研究提出了孔隙度与埋深具有指数关系,而通过研究发现砂岩也具有类似关系式[4];刘震等[5-6]提出孔隙度是埋深和时间的双元函数,并以此建立了孔隙度演化模型,但他们是将各类胶结物作了统一处理,模糊了产生各类胶结物的阶段性,因有些胶结物可以在特定时期内大量产生,所以,孔隙度与时间可能并不简单的符合双元函数关系;孟元林等[7]引入成岩指数,用于在钻前对孔隙度进行预测,但并未分析各种成岩作用对孔隙度的影响;王瑞飞等[8]探讨了各种成岩作用对孔隙度的影响,但并未对成岩序列进行系统分析。在大量分析测试资料的基础上,根据铸体薄片、扫描电镜、阴极发光、高压压汞等分析,对研究区致密砂岩储层进行孔隙度演化分析,并建立孔隙度演化模型,以期为有利储层的预测提供依据。
1 地质概况神府地区位于鄂尔多斯盆地东北部,横跨晋西挠褶带和伊陕斜坡2个构造单元[图 1(a)],自下而上发育马家沟组(O1m)、本溪组(C2b)、太原组(P1t)、山西组(P1s)、石盒子组(P2sh)、石千峰组(P3sh)、刘家沟组(T1l)、和尚沟组(T2h)、纸坊组(T2z)、延长组(T3y)、第四系(Q)等。研究区太原组沉积时期地势平坦,水体较浅,水动力较弱,潮坪相分布广泛,发育泥坪、灰泥坪、砂坪、潮上沼泽、潮汐水道、障壁砂坝等沉积微相,在垂向上形成砂岩、泥岩、砂泥互层与煤层交互出现的岩性组合,并在砂体中发育大量流水型交错层理和羽状交错层理,整体属于低能环境。泥坪和灰泥坪以泥岩沉积为主,部分夹灰质泥岩;砂坪以粉砂岩和细砂岩沉积为主,部分夹砂泥岩薄互层;潮上沼泽以煤岩和炭质泥岩沉积为主;潮汐水道和障壁砂坝主要发育粗砂岩和中—粗砂岩,是该区较优质的储层[图 1(b)]。
统计10口井岩心和72块薄片资料发现,神府地区下二叠统太原组主要发育中—粗岩屑砂岩和长石岩屑砂岩[图 2、图 3(a)~(b)]。石英体积分数为30.3%~78.5%,平均为52.5%,长石体积分数为1.3%~18.1%,平均为10.1%;岩屑体积分数为1.1%~15.3%,平均为35.7%。长石风化程度中等,钾长石多发生高岭石化,斜长石多发生绢云母化。岩屑类型以变质岩岩屑和火成岩岩屑为主,含少量的千枚岩岩屑和泥岩岩屑。太原组砂岩填隙物体积分数约为25.86%,其中泥质杂基体积分数约为15.47%,占填隙物总量的59.9%,其胶结物主要为白云石、铁白云石、菱铁矿和高岭石等,而方解石含量较低(表 1)。颗粒分选中等—好,一般呈次圆—次棱角状,磨圆度中等—好,以颗粒支撑为主,颗粒间多为线接触,胶结类型属于孔隙式胶结。岩石结构成熟度中等,成分成熟度中等—高。
统计神府地区16口井431块钻井取心样品的物性资料可知,太原组砂岩储层孔隙度为0.2%~12.0%,平均为6.95%,渗透率为0.1~5.0 mD,平均为0.91 mD,属于较典型的致密砂岩储层。由图 4可知,当孔隙度小于6%时,部分样品的孔、渗数据比较分散,但渗透率普遍较低( < 10 mD);根据数据分布规律,拟合2条边界函数曲线,将偏离的数据分为Ⅰ,Ⅱ共2个区域。对这2个区域的岩心样品进行观察时,并未发现明显裂缝,但在进行薄片观察时发现:Ⅰ区样品含少量砾石,整体粒度偏粗,砾石边缘颗粒溶蚀强烈[图 5(a)],粒间溶孔较为发育,相邻的粒内溶孔边界往往被“溶开”而彼此连通,从而提高了岩石的渗透性;Ⅱ区样品孔隙主要为长石溶孔,孔隙中存在长石被溶蚀后残留的微晶高岭石[图 5(b)],造成喉道阻塞,使岩石渗透性降低。当孔隙度大于6%时(目前太原组实测最大孔隙度未超过16%),岩石孔隙以原生粒间孔和粒间溶孔为主,其次为少量的长石粒内溶孔,且孔隙内一般无微晶高岭石残留,孔隙度和渗透率呈较好的正相关关系,表明微裂缝对于储层渗流贡献较小,渗透率主要受基质孔喉发育程度的影响。
太原组砂岩储层具有高排驱压力、低中值半径和低退汞效率的特征(表 2)。统计研究区43块压汞资料发现:神府地区太原组砂岩储层排驱压力为0.39~18.59 MPa,平均为3.05 MPa,排驱压力较大;中值半径为0.05~0.22 μm,平均为0.11 μm,属于细喉、特细喉道型[10];退汞效率较低,为24.24%~42.64%,平均为34.75%,表明其孔隙连通性较差。毛管压力曲线能够比较直观地反映储层的孔隙结构[11-13],由图 6可见,神府地区太原组储层孔隙结构较差,属于小孔细喉和小孔微喉型孔隙结构。
神府地区太原组储层在埋藏成岩过程中经历了压实、胶结、溶蚀等成岩作用,三者对储层物性的变化均具有重要影响。
3.1.1 压实作用通过铸体薄片观察,发现太原组储层碎屑颗粒之间的接触关系多为线接触,部分样品出现线—凹凸接触[图 7(a)],此外,岩石组分中刚性颗粒石英含量较低,而塑性颗粒岩屑含量较高,如酸性喷发岩岩屑、千枚岩岩屑、片岩岩屑和云母等在岩石中的平均体积分数分别为9.75%,2.60%,0.46%和0.71%,其总和在岩屑中平均占比可达40%,这是导致太原组储层压实作用较强的重要因素。较强的压实作用使太原组储层原生孔隙大量损失,只保留少量残余原生孔,仅占现有孔隙的13.5%,且孔隙多呈三角形或不规则多边形[图 7(b)]。
神府地区太原组砂岩胶结物类型多样,发育3期碳酸盐胶结物、2期硅质胶结物和3期黏土类胶结物等。3期碳酸盐胶结物分别为:早期方解石[图 7(c)]和菱铁矿;中期白云石;晚期铁白云石和铁方解石,并可见铁白云石交代白云石[图 7(d)]。2期硅质胶结物分别为:早期石英次生加大,由压溶作用产生的硅质所形成;晚期石英次生加大,由长石溶解所产生的硅质所形成,并可见方解石包裹早期石英次生加大[图 7(e)]。3期黏土类胶结物分别为:早期伊/蒙混层(蒙皂石)黏土包壳[图 7(f)];中期长石溶蚀所产生的高岭石[图 7(g)];晚期呈发丝状分布在矿物颗粒表面的伊利石[图 7(h)],局部可见高岭石向伊利石转化以及伊利石包裹在自生石英颗粒表面。神府地区太原组储层中胶结物多呈团块状或不规则状充填于粒间孔隙,并见部分黏土矿物包裹于碎屑颗粒表面,对储层的影响以减孔作用为主。
3.1.3 溶蚀作用溶蚀作用可以使长石和酸性喷发岩岩屑等岩石组分溶解,产生次生孔隙,提高储层孔隙度,但对于储层渗透率的影响未必都是有利的。当储层矿物溶蚀作用强烈时,次生孔隙以粒间溶孔为主,粒内溶孔往往因溶蚀强烈而彼此连通,且长石溶蚀所产生的高岭石微晶也容易被流体带走[14],从而使岩石渗流能力增强;当储层矿物溶蚀程度较低时,储层次生孔隙多以孤立形式存在,长石溶蚀所产生的高岭石微晶因无法完全被孔隙流体带走而导致部分流体被残留在孔隙内,当后期储层温度、压力等环境条件改变后,高岭石微晶结晶析出,从而造成喉道阻塞,反而会降低储层的渗流能力。
溶蚀程度与储层岩性有关,同时也受到流体性质、温度和压力等条件的影响[15-16],且煤系地层含有丰富的有机质,在成岩作用早期能够产生大量有机酸,这对于储层次生孔隙的产生也具有较有利影响。神府地区太原组致密砂岩储层中各类长石和酸性喷发岩岩屑[图 7(i)]等易溶组分含量较高,在部分样品中其易溶组分体积分数可高达38%,且区内煤系地层较为发育,为该区强烈溶蚀作用的发生奠定了物质基础。
3.2 成岩序列通过对神府地区太原组泥岩镜质体反射率(Ro)、有机质最高裂解温度(Tmax)及伊/蒙混层中蒙脱石含量等的分析发现,太原组泥岩Ro为0.84%~1.21%,Tmax为450~490 ℃,通常在460 ℃左右,伊/蒙混层中蒙脱石体积分数为12%~35%。由碎屑岩成岩阶段划分标准[17]可知,神府地区太原组成岩作用阶段普遍在中成岩阶段A期,部分可达中成岩阶段B期(图 8)。根据自生矿物之间及其与碎屑颗粒之间的接触关系,确定神府地区太原组典型成岩序列为:早期压实和黏土包壳→Ⅰ期石英加大→菱铁矿和方解石胶结→长石溶蚀→自生高岭石和Ⅱ期石英加大→自生伊利石胶结→白云石、铁白云石及铁方解石胶结。碎屑颗粒首先经历早期的机械压实,由于神府地区太原组泥质杂基含量较高,碎屑颗粒常伴有黏土包壳,随着压实程度的进一步加强,碎屑颗粒发生压溶作用,其溶解所产生的二氧化硅形成了Ⅰ期的石英加大;随后发生了早期的方解石沉淀和菱铁矿的形成[18],虽然下二叠统储层早期受煤系地层影响,整体处于酸性环境[19],但由于储层非均质性强,局部仍然保留了少量的方解石胶结物;在有机酸大量充注后,酸性喷发岩岩屑和长石发生溶蚀所产生的矿物离子促进了自生高岭石和Ⅱ期石英加大的形成;随着成岩演化程度的加深,黏土矿物向伊利石转化,地层水逐渐变为碱性,含铁的碳酸盐矿物如铁方解石和铁白云石开始形成。
通过对成岩序列的研究并结合有机酸的充注过程[20-21],将胶结物划分为早期、中期和晚期3类。方解石、菱铁矿和以石英次生加大为主的硅质胶结划为早期胶结物;高岭石的生成与有机酸大量充注时间一致,归为中期胶结物;铁方解石、白云石、铁白云石和自生伊利石等的胶结主要发生在有机酸充注以后,划为晚期胶结物。
3.3 孔隙度演化结合成岩序列、薄片和物性等资料,对神府地区太原组砂岩储层孔隙度演化过程进行研究。
初始孔隙度的确定
$ {{\varphi }_{0}}=\text{ }20.91+22.90/{{S}_{0}} $ | (1) |
式中:φ0为初始孔隙度,%;S0为分选系数,由S0 = (Q1/Q2)1/2计算得出,Q1和Q2为粒度累积曲线在25%和75%处的粒径值[22-23]。神府地区太原组分选系数S0为1.57~4.95,平均为2.95,因此计算太原组初始孔隙度约为28.67%。
压实减孔量利用下式计算
$ {{\varphi }_{压}}={{\varphi }_{0}}-{{\varphi }_{粒}} $ | (2) |
$ {{\varphi }_{粒}}={{\varphi }_{胶}}+{{\varphi }_{胶后}} $ | (3) |
$ {{\varphi }_{胶后}}={{\varphi }_{残原}}+{{\varphi }_{粒间}}\cdot (100\%-\omega )+{{\varphi }_{胶溶}} $ | (4) |
式中:φ压为压实作用减小的孔隙度,%;φ粒为压实作用后粒间剩余孔隙度[24],%;φ胶为胶结物总量,%;φ胶后为胶结后的粒间孔隙度,%;φ残原为现今薄片中残余原生孔隙度,%;φ粒间为粒间溶孔所产生的孔隙度,%;φ胶溶为胶结物溶蚀孔隙度,%;ω为颗粒边缘溶蚀率,%。
由于粒间溶孔通常会包括颗粒边缘的溶蚀量和残余原生孔隙2个部分,所以引入参数ω,代表颗粒边缘溶蚀量与粒间溶孔的比值,通过对太原组大量铸体薄片观察统计,将ω的数值确定为20%。
溶蚀增孔量通过下式计算
$ {{\varphi }_{溶增}}={{\varphi }_{粒内}}+{{\varphi }_{粒间}}\cdot \omega +{{\varphi }_{胶溶}} $ | (5) |
式中:φ溶增为溶蚀作用增加的孔隙度,%;φ粒内为粒内溶孔所产生的孔隙度,%。
胶结减孔量即为胶结物总量,包括各类碳酸盐胶结物、硅质胶结物和各类自生黏土矿物。
利用实测孔隙度平均值与总面孔率平均值的比值,搭建起面孔率[式(2)~式(5)皆是通过薄片观察统计获得的面孔率]与孔隙度的桥梁[25],进而计算出以上各类孔隙度的结果(表 3)。
通过数据计算得知:太原组储层经过压实作用,损失孔隙15.50%,损失率达54.1%;压实作用后储层剩余孔隙度约为13.17%;经过早期胶结作用(硅质、方解石和菱铁矿胶结),减少孔隙2.51%,经有机酸充注,岩石中易溶组分发生溶解,溶蚀作用增加孔隙4.41%,增孔率达15.4%;经过中期胶结作用(高岭石胶结),孔隙减少1.05%,晚期胶结作用(白云石、铁白云石和铁方解石胶结)减少孔隙6.84%,胶结作用共减少孔隙10.40%,减孔率达36.3%(表 4)。经过上述多种成岩作用改造后,太原组砂岩储层平均孔隙度达到7.17%,与岩心实测平均孔隙度6.95%相比,相对误差仅为3.20%。早期压实作用与晚期碳酸盐胶结作用使储层损失了绝大部分孔隙,溶蚀作用所产生的次生孔隙是现今储层的主要储集空间类型,其中,又以具有较高渗流能力的粒间溶孔为较好的孔隙类型。
(1)神府地区太原组砂岩储层以中—粗粒岩屑砂岩和长石岩屑砂岩为主,孔隙结构主要为细孔、细喉型和特细孔、微喉型,储层物性较差,孔隙度和渗透率平均值分别为6.95%和0.91 mD,属于典型的低孔—特低孔、特低渗的致密砂岩储层。
(2)神府地区太原组砂岩储层在沉积后主要经历了压实、胶结、交代及溶解等成岩作用的改造。早期压实作用与晚期碳酸盐胶结作用是储层致密的主要原因。溶蚀作用所产生的次生孔隙是现今储层的主要储集空间类型,其中,又以粒间溶孔为较好的孔隙类型。
[1] |
姚泾利, 唐俊, 庞国印, 等. 鄂尔多斯盆地白豹-华池地区长8段孔隙度演化定量模拟. 天然气地球科学, 2013, 24(1): 38-46. YAO J L, TANG J, PANG G Y, et al. Quantitative simulation on porosity-evolution in member 8 of Yanchang Formation of Baibao-Huachi area, Ordos Basin. Natural Gas Geoscience, 2013, 24(1): 38-46. |
[2] |
廖朋, 唐俊, 王凯, 等. 砂岩成岩过程中的孔隙演化定量模拟——以鄂尔多斯盆地安塞地区长8油层组储层为例. 岩性油气藏, 2014, 26(5): 15-22. LIAO P, TANG J, WANG K, et al. Quantitative simulation on pore evolution in diagenetic process of sandstone:a case study from Chang 8 oil reservoir set in Ansai area, Ordos Basin. Lithologic Reservoirs, 2014, 26(5): 15-22. |
[3] |
ATHY L F. Density, porosity, and compaction of sedimentary rocks. AAPG Bulletin, 1930, 14(1): 1-24. |
[4] |
陈发景, 田世澄. 压实与油气运移. 武汉: 中国地质大学出版社, 1989, 113-134. CHEN F J, TIAN S C. Compaction and oil or gas migration. Wuhan: China University of Geoscience Press, 1989, 113-134. |
[5] |
刘震, 邵新军, 金博, 等. 压实过程中埋深和时间对碎屑岩孔隙度演化的共同影响. 现代地质, 2007, 21(1): 125-132. LIU Z, SHAO X J, JIN B, et al. Co-effect of depth and burial time on the evolution of porosity for clastic rocks during the stage of compaction. Geoscience, 2007, 21(1): 125-132. |
[6] |
潘高峰, 刘震, 赵舒, 等. 砂岩孔隙度演化定量模拟方法——以鄂尔多斯盆地镇泾地区延长组为例. 石油学报, 2011, 32(2): 249-256. PAN G F, LIU Z, ZHAO S, et al. Quantitative simulation of sandstone porosity evolution:a case from Yanchang Formation of the Zhenjing area, Ordos Basin. Acta Petrolei Sinica, 2011, 32(2): 249-256. DOI:10.7623/syxb201102009 |
[7] |
孟元林, 王粤川, 牛嘉玉, 等. 储层孔隙度预测与有效天然气储层确定——以渤海湾盆地鸳鸯沟地区为例. 天然气工业, 2007, 27(6): 42-44. MENG Y L, WANG Y C, NIU J Y, et al. Prediction of reservoir porosity and determination of effective gas reservoirs:taking Yuanyanggou area of Bohai Bay Basin as an example. Natural Gas Industry, 2007, 27(6): 42-44. |
[8] |
王瑞飞, 沈平平, 赵良金. 深层储集层成岩作用及孔隙度演化定量模型——以东濮凹陷文东油田沙三段储集层为例. 石油勘探与开发, 2011, 38(5): 552-559. WANG R F, SHEN P P, ZHAO L J. Diagenesis of deep sandstone reservoir and a quantitative model of porosity evolution:taking the third member of Shahejie Formation in the Wendong Oilfield, Dongpu Sag as an example. Petroleum Exploration and Development, 2011, 38(5): 552-559. |
[9] |
潘新志, 叶建平, 孙新阳, 等. 鄂尔多斯盆地神府地区中低阶煤层气勘探潜力分析. 煤炭科学技术, 2015, 43(9): 65-70. PAN X Z, YE J P, SUN X Y, et al. Analysis on exploration potential of mid-low rank coalbed methane in Shenfu area of Ordos Basin. Coal Science and Technology, 2015, 43(9): 65-70. |
[10] |
樊爱萍, 赵娟, 杨仁超, 等. 苏里格气田东二区山1段、盒8段储层孔隙结构特征. 天然气地球科学, 2011, 22(3): 482-487. FAN A P, ZHAO J, YANG R C, et al. Pore structure of reservoir rocks in Shan 1 and He 8 members, the EastⅡblock of Sulige gas field. Natural Gas Geoscience, 2011, 22(3): 482-487. |
[11] |
郝乐伟, 王琪, 唐俊. 储层岩石微观孔隙结构研究方法与理论综述. 岩性油气藏, 2013, 25(5): 123-128. HAO L W, WANG Q, TANG J. Research progress of reservoir microscopic pore structure. Lithologic Reservoirs, 2013, 25(5): 123-128. |
[12] |
李永胜, 刘学刚, 章志锋, 等. 鄂尔多斯盆地姬塬油田长8储层微观孔隙结构特征. 石油化工应用, 2013, 32(4): 63-69. LI Y S, LIU X G, ZHANG Z F, et al. Study microscopic pore structure characteristics of Chang 8 reservoir of Jiyuan oilfield by capillary curves. Petrochemical Industry Application, 2013, 32(4): 63-69. |
[13] |
张宪国, 张涛, 林承焰. 基于孔隙分形特征的低渗透储层孔隙结构评价. 岩性油气藏, 2013, 25(6): 40-45. ZHANG X G, ZHANG T, LIN C Y. Pore structure evolution of low permeability reservoir based on pore fractal features. Lithologic Reservoirs, 2013, 25(6): 40-45. |
[14] |
曹剑, 张义杰, 胡文瑄, 等. 油气储层自生高岭石发育特点及其对物性的影响. 矿物学报, 2005, 25(4): 367-373. CAO J, ZHANG Y J, HU W X, et al. Developing characteristics of kaolinite in central Junggar Basin and their effect on the reservoir quality. Acta Mineralogica Sinica, 2005, 25(4): 367-373. |
[15] |
SURDAM R C, CROSSEY L J, HAGEN E S, et al. Organicinorganic interactions and sandstone diagenesis. AAPG Bulletin, 1989, 73(1): 1-23. |
[16] |
HUNT J M. Generation and migration of petroleum from abnormally pressured fluid compartments. AAPG Bulletin, 1990, 74(1): 1-12. |
[17] |
国家经济贸易委员会. 碎屑岩成岩阶段划分标准:SY/T 5477-2003. 北京: 石油工业出版社, 2003, 1-4. State Economic and Trade Commission. The division of diagenetic stages in clastic rocks:SY/T5477-2003. Beijing: Petroleum Industry Press, 2003, 1-4. |
[18] |
张莹莹, 黄思静. 华庆地区长6油层组方解石胶结物特征. 岩性油气藏, 2012, 24(2): 48-52. ZHANG Y Y, HUANG S J. Characteristics of calcite cements of Chang 6 oil reservoir set in Huaqing area. Lithologic Reservoirs, 2012, 24(2): 48-52. |
[19] |
朱国华, 章卫平. 煤系地层砂岩成岩作用和孔隙演化研究——以长广地区龙潭组为例. 石油勘探与开发, 1993, 20(1): 42-50. ZHU G H, ZHANG W P. A study of digenesis and the evolution of porosity of the sandstones coaliferous Formations:taking Longtan group in Changguang region as an example. Petroleum Exploration and Development, 1993, 20(1): 42-50. |
[20] |
郑浚茂, 应凤祥. 煤系地层(酸性水介质)的砂岩储层特征及成岩模式. 石油学报, 1997, 18(4): 19-24. ZHENG J M, YING F X. Reservoir characteristics and diagenetic model of sandstone intercalated in coal-bearing strata. Acta Petrolei Sinica, 1997, 18(4): 19-24. DOI:10.7623/syxb199704004 |
[21] |
远光辉, 操应长, 杨田, 等. 论碎屑岩储层成岩过程中有机酸的溶蚀增孔能力. 地学前缘, 2013, 20(5): 207-219. YUAN G H, CAO Y C, YANG T, et al. Porosity enhancement potential through mineral dissolution by organic acids in the diagenetic process of clastic reservoir. Earth Science Frontiers, 2013, 20(5): 207-219. |
[22] |
BEARD D C, WEYL P K. Influence of texture on porosity and permeability of unsolidated sand. AAPG Bulletin, 1973, 57(2): 349-369. |
[23] |
SCHERER M. Parameters influencing porosity in sandstones a model for sandstone porosity prediction. AAPG Bulletin, 1987, 71(5): 485-491. |
[24] |
张兴良, 田景春, 王峰, 等. 致密砂岩储层成岩作用特征与孔隙演化定量评价——以鄂尔多斯盆地高桥地区二叠系下石盒子组盒8段为例. 石油与天然气地质, 2014, 35(2): 212-217. ZHANG X L, TIAN J C, WANG F, et al. Diagenetic characteristics and quantitative porosity estimation of tight sandstone reservoirs:a case from the 8 th member of Permian Xiashihezi Formation in the Gaoqiao region, Ordos Basin. Oil & Gas Geology, 2014, 35(2): 212-217. |
[25] |
思玉琥, 郝世彦, 张林森, 等. 延安地区上三叠统长6期储层成岩作用及孔隙演化. 特种油气藏, 2011, 18(6): 36-39. SI Y H, HAO S Y, ZHANG L S, et al. Diageneses of the Upper Triassic Chang 6 Formation in Yan' an area and pore evolution. Special Oil & Gas Reservoirs, 2011, 18(6): 36-39. |