岩性油气藏  2019, Vol. 31 Issue (5): 70-83       PDF    
×
低渗—特低渗白云岩储层成岩相分析及测井识别——以伊陕斜坡马五段为例
李百强1, 张小莉1, 王起琮2, 郭彬程3, 郭艳琴2, 尚晓庆2, 程浩2, 卢俊辉1, 赵希1    
1. 西北大学 地质学系/二氧化碳捕集与封存技术国家地方联合工程研究中心, 西安 710069;
2. 西安石油大学 地球科学与工程学院, 西安 710065;
3. 中国石油勘探开发研究院, 北京 100083
摘要: 为研究低渗—特低渗白云岩储层的成岩相测井识别方法并进行定量划分,明确白云岩成岩相的空间展布及其有利区,依据岩心观察、普通薄片和铸体薄片鉴定、以及阴极发光分析等,结合X射线衍射和物性资料,对伊陕斜坡中东部马五段白云岩储层的成岩环境及成岩作用进行剖析,通过测井曲线的交会分析与软件的综合应用,构建了低渗—特低渗白云岩储层成岩相的定量识别标准。结果表明,低渗—特低渗白云岩储层成岩相可划分为早期大气淡水溶蚀亚相、表生期大气淡水溶蚀亚相、浅埋藏活跃回流渗透云化亚相、浅埋藏隐伏回流渗透云化亚相、钙质胶结相和泥质充填相等6类;以Pe-DENRLLD-DEN交会及RLLD-AC交会为代表的测井参数交会图法及软件综合应用是进行低渗—特低渗白云岩储层不同成岩相单井纵向上定量识别与连续划分的有效方法;白云岩储层成岩相在空间的分布具有纵向分带、横向分区的典型特征。有利储层主要分布于早期大气淡水溶蚀亚相和浅埋藏活跃回流渗透云化亚相,其次为表生期大气淡水溶蚀亚相,浅埋藏隐伏回流渗透云化亚相储层分布最少。
关键词: 低渗-特低渗      白云岩储层      成岩相      测井识别      马五段      鄂尔多斯盆地     
Analysis and logging identification of diagenetic facies of dolomite reservoir with low and ultra-low permeability:a case study from Ma 5 member in Yishan slope, Ordos Basin
LI Baiqiang1, ZHANG Xiaoli1, WANG Qicong2, GUO Bincheng3, GUO Yanqin2, SHANG Xiaoqing2, CHENG Hao2, LU Junhui1, ZHAO Xi1     
1. Department of Geology, Northwest University/National & Local Joint Engineering Research Center of Carbon Capture and Storage Technology, Xi'an 710069, China;
2. School of Earth Sciences and Engineering, Xi'an Shiyou University, Xi'an 710065, China;
3. PetroChina Research Institute of Petroleum Exploration and Development, Beijing 100083, China
Abstract: In order to study the logging identification method and quantitative classification of diagenetic facies of dolostone reservoir with low and ultra-low permeability and clarify the spatial distribution characteristics of diagenetic facies and the favorable exploration area of dolostone reservoir". Based on core observation, identification of ordinary and casting thin section, and cathodoluminescence analysis, combined with X-ray diffraction and physical property data, the diagenetic environment and diagenesis of dolomite reservoir with low and ultralow permeability of Ma 5 member in Yishan slope were analyzed. Then the quantitative identification standard for diagenetic facies of dolomite reservoir with low and ultra-low permeability was established through the crossplot analysis of conventional logging curves and application of software program. The results show that the diagenetic facies of the dolomite reservoirs with low and ultra-low permeability can be divided into six types:early atmospheric freshwater dissolution, supergene atmospheric freshwater dissolution, shallow buried active reflux seepage dolomitization, shallow buried latent reflux seepage dolomitization, calcareous cementation and clay minerals filling. The crossplot method of logging parameters represented by Pe and DEN, RLLD and DEN, and RLLD and AC, and the comprehensive application of software program are effective approach for vertical quantitative identification and continuous division of different diagenetic facies of dolostone reservoir with low and ultralow permeability in single well. The spatial distribution of diagenetic facies of dolostone reservoir shows a typical characteristic of longitudinal zonation and planar partition. The favorable reservoirs are mainly developed in diagenetic facies of early atmospheric freshwater dissolution and shallow buried active reflux seepage dolomitization, followed by supergene atmospheric freshwater dissolution, and there is less reservoir deve-loped in shallow buried latent reflux seepage dolomitization diagenetic subfacies.
Key words: low and ultra-low permeability      dolomite reservoir      diagenetic facies      logging identification      Ma 5 member      Ordos Basin     
0 引言

成岩作用是低渗—特低渗碳酸盐岩储集空间最终形成及定格的主要因素,具体表现为成岩作用的类型和强度控制了储集空间的形成和展布[1-3],同时,由于碳酸盐岩成岩作用十分复杂,造成储层具有较强的非均质性[4]。因此,寻找能够有效解决控制复杂成岩作用的碳酸盐岩储层展布预测问题的新途径和新方法就显得尤为重要[5-6]。成岩相的概念最早由美国学者Railsback [7]提出,并被认为能够有效解决碳酸盐岩储层非均质性问题。随后,陈彦华等[5]引进了成岩相的概念,并将其定义为成岩环境和成岩产物的综合,是反映成岩环境的物质表现。经过众多学者对碳酸盐岩成岩相的不断探索,认为只有搞清楚成岩相的展布特征及规律,才能深刻理解和研究储集体的内部非均质性及其分布特征。因此,对于强非均质性的低渗—特低渗白云岩储层,要明确其储层的展布规律,就必须先对其成岩相进行有效识别及预测。

目前,预测成岩相的方法为实验室岩心测试,包括薄片和铸体薄片鉴定、扫描电镜和阴极发光分析等技术[8]。其优点是直观、具体、准确度高,但因取样成本高、消耗时间长、取心资料较少等问题,无法实现成岩相在空间上的连续评价,从而给储层综合评价带来了一定的困难[9]。已有学者采用测井技术对致密砂岩储层成岩相进行了识别和评价,并初步取得了较为满意的效果[8-10],而涉及低渗透、特低渗透碳酸盐岩储层相对较少,因此,低渗—特低渗白云岩储层成岩相的测井识别可行性及可靠性也急须得到实例验证。以伊陕斜坡马五5—马五1亚段为例,通过开展低渗—特低渗白云岩储层成岩作用的研究,并进行成岩相的分类评价以及测井识别,明确储层成岩相的空间展布特征,以期能够为白云岩储层的有利区预测提供更可靠的依据以及为低渗—特低渗碳酸盐岩储层成岩相的研究提供一定的借鉴。

1 区域地质概况

研究区构造上位于鄂尔多斯盆地伊陕斜坡中东部,北起大牛地、南至富县、西起靖边、东抵延川。在岩溶古地貌单元上位于岩溶斜坡带内(图 1[11]。马家沟组沉积末期,加里东运动导致地层抬升至地表并经历了长期的暴露、溶蚀及淋滤,形成了一套储、渗能力较强、天然气相对富集的风化壳型碳酸盐岩储层,马五段位于该风化壳的顶部,因此也是马家沟组重要的产气层位[12]

下载eps/tif图 图 1 研究区位置 Fig. 1 Location of the study area

从古地理环境来看,研究区马五5—马五1沉积环境经历了开阔海台地—蒸发台地—局限台地的转变[13]。储层也经历了早期淡水溶蚀、浅埋藏期白云石化以及表生期淡水溶蚀等复杂成岩作用的改造,形成的储层岩石类型包括膏盐溶蚀泥—粉晶白云岩、残余砂屑粉—细晶白云岩、豹斑白云岩以及岩溶角砾泥—粉晶白云岩,储集空间类型为各类铸模孔、次生溶蚀孔、晶间孔、砾间孔以及微裂缝等。孔隙度整体为0.3%~9.0%,平均为2.5%,渗透率为0.003~10.000 mD,平均为0.1 mD,属于典型的低孔低渗—特低渗储层[2]

2 白云岩成岩相分类

关于碳酸盐岩成岩相的分类目前主要有2种观点:①成岩相主要指成岩作用的组合,从而提出了胶结-溶蚀、压溶-白云化等成岩相类型;②成岩相即是成岩环境的反映,因此总结出了准同生相、混合水相以及表生相等[14-16]类型。杨威等[1]认为,这2种观点均有缺陷,主要是由于同一成岩环境下可以发生不同的成岩作用,所以成岩环境不能完全反映成岩作用特征;在不同成岩环境下也可能出现相同的成岩作用,但作用的机理及对储层的意义差别很大,并提出了碳酸盐岩成岩相应与沉积相相区分,即成岩相是岩石在一定的成岩环境下各类成岩作用的综合表现。因此,划分碳酸盐岩成岩相,既要综合考虑成岩环境和成岩作用的双重效应,也要有利于进行成岩相的定性和定量识别。依据不同的成岩作用类型,可以将碳酸盐岩成岩相划分为压实相、压溶相、溶蚀相、胶结相、云化相、去云化相、膏化相、去膏化相、重结晶相、破裂相以及充注相[1]。再结合成岩环境的不同,可以将溶蚀相进一步划分为大气淡水溶蚀亚相和埋藏溶蚀亚相,将胶结相进一步划分为海底胶结亚相、淡水胶结亚相以及埋藏胶结亚相,将云化相进一步划分为准同生云化亚相、埋藏云化亚相以及混合水云化亚相[1]

采用岩心观察、铸体薄片鉴定、扫描电镜及阴极发光分析等方法,考虑岩石的结构及构造、储集空间类型等,研究马五5—马五1亚段低渗—特低渗白云岩的形成环境及其所经历的成岩作用,进而进行储层成岩相的分类。首先识别出白云岩储层主要发育溶蚀相、云化相、胶结相和充填相等四大类成岩相以及大气淡水溶蚀亚相和埋藏云化亚相2类成岩亚相。其次根据岩石的结构及构造、储集空间类型将大气淡水溶蚀亚相细分为早期大气淡水溶蚀亚相和表生期大气淡水溶蚀亚相2类,将埋藏云化亚相划分为浅埋藏活跃回流渗透云化亚相和浅埋藏隐伏回流渗透云化亚相2种,胶结相主要为钙质胶结相,充填相则以泥质充填相为主。

2.1 溶蚀相

(1)早期大气淡水溶蚀亚相。这是指形成于早期近地表淡水环境中,泥—微晶白云岩中的一部分孔隙被硬石膏及石盐组分充填或胶结,当在表生期时暴露于地表,岩溶作用使得部分石膏或石盐组分再次发生组构选择性溶解,从而形成膏模孔或盐模孔,而未被溶解的组分则仍以胶结物形式存在于白云岩中,从而形成膏盐溶蚀泥—微晶/粉晶白云岩等一系列过程的成岩环境和成岩作用及其产物的组合。这类成岩相中,代表性的岩石类型为膏盐溶蚀泥—微晶/粉晶白云岩[图 2(a)],X射线衍射分析结果表明,矿物组分以泥—微晶白云石和泥—粉晶白云石为主,其中泥—微晶白云石主要为他形晶,白云石有序度为0.72~0.73,而泥—粉晶白云石以自形和半自形晶为主,白云石有序度一般为0.75~1.00。储集空间以膏模孔和盐模孔等铸模孔[图 2(b)]为主,同时还有一些未被充填的示底构造所形成的示底孔隙。由于受到早期淡水的溶蚀作用,在阴极发光镜下,孔隙一般不发光或弱发光,而岩石基质一般发光较明亮[图 2(c)]。示底构造也是早期大气淡水溶蚀亚相的典型岩石构造证据,其下部一般为成岩阶段的物理成因的沉积物,反映在溶解过程中的重力垮塌物,上部可以为全开放或部分开放的孔隙,也可被胶结物所填充。这类成岩亚相是研究区内近地表环境的一类重要的建设性成岩相。

下载eps/tif图 图 2 伊陕斜坡马五段白云岩储层成岩相的岩石学识别标志 (a)泥晶白云岩,FG4井,3 082.03 m,马五3,单偏光;(b)泥晶白云岩中盐模孔隙,FG7井,2 960.55 m,马五1,铸体薄片;(c)受早期大气淡水影响,泥晶白云岩基质呈明亮的橙红色,盐模孔隙为暗色不发光,S157井,3 311.81 m,马五1,阴极发光;(d)岩溶角砾泥—粉晶白云岩,FG7井,2 976.2 m,马五2,单偏光;(e)岩溶角砾间孔隙,S115井,3 342.20 m,马五2,铸体薄片;(f)岩溶角砾泥晶白云岩,角砾不发光或很暗,角砾间填隙物受淡水影响发亮橙色光,S115井,3 342.20 m,马五2,阴极发光;(g)粉—细晶白云岩中的白云石晶间孔,D48井,3 014.0 m,马五5,铸体薄片;(h)粉—细晶白云岩中白云石的亮边结构,D48井,3 016.17 m,马五5,阴极发光;(i)豹斑白云岩中的豹斑构造,其中图片左半部分为豹斑(成分为白云石晶体),右侧为基质(成分为方解石),D78井,2 887.90 m,马五5,单偏光 Fig. 2 Petrological identification marks of diagenetic facies of dolomite reservoir of Ma 5 member in Yishan slope

(2)表生期大气淡水溶蚀亚相。这是在近地表淡水溶蚀作用下,与易溶类的膏盐层或盐岩层相伴生的泥—粉晶白云岩被长期、大范围地非组构选择性地溶解形成大量溶蚀孔、洞,并伴随着岩石的坍塌破裂形成角砾结构及其形成岩溶角砾泥—粉晶白云岩[图 2(d)]的过程中,成岩环境与成岩作用及其产物的组合,以岩溶角砾泥—粉晶白云岩为代表。X射线衍射分析结果表明,泥—粉晶白云石以他形晶为主,白云石有序度变化范围较大,一般为0.54~1.00。主要的储集空间类型为溶解作用形成的溶蚀孔、洞以及岩石坍塌破裂及角砾化过程中形成的角砾间孔[图 2(e)]。由于角砾很少遭受淡水影响,而角砾间孔隙中的填隙物或胶结物受淡水影响较严重,因此在阴极发光镜下,角砾不发光或发光很弱,而角砾间孔填隙物或胶结物发明亮的橙红色光[图 2(f)],这也是判断表生期岩溶作用的一个重要标志。这类成岩相为建设性成岩亚相,而其所形成的储集层物性相比早期大气淡水成岩相要差一些。

2.2 云化相

(1)浅埋藏活跃回流渗透云化亚相。这是指形成于近地表—浅埋藏超咸水环境中,蒸发台地内的浓卤水在重力和密度差的驱使下向下部渗透性较好的颗粒灰岩回流渗透并交代其中的灰质组分而发生白云化作用,最终形成残余砂屑粉—细晶白云岩,在此过程中的成岩环境、成岩作用及其产物的组合。这类成岩相中,代表性的岩石类型为残余砂屑粉—细晶白云岩。此类白云岩中白云石的盐度为2.46%~2.69%,平均为2.59%。X射线衍射分析结果表明,粉—细晶白云石自形程度相对较高(半自形—自形),白云石有序度为0.75~0.85。储集空间以白云石晶间孔[图 2(g)]为主,此类孔隙连通性较好,具有良好的储集性能。虽然活跃回流渗透白云化作用相对彻底,但是其属于组构选择性白云岩化作用,导致经过白云化流体交代后的灰岩的颗粒结构仍清晰可见,因此交代残余结构是识别这类成岩相的重要标志。此外,在阴极发光镜下,也可见到由于白云化流体交代过程中残留的灰泥质赋存于白云石表面,从而形成典型的雾心亮边结构[图 2(h)],这也是识别这类成岩亚相的重要依据。该成岩亚相是研究区进入浅埋藏环境后的一类重要的建设性成岩相。

(2)浅埋藏隐伏回流渗透云化亚相。这是指形成于近地表—浅埋藏变盐度海水环境中,蒸发台地内的变盐度海水(台地内的浓卤水遭受再一次的海侵所形成的盐度和密度均低于浓卤水但仍旧高于海水的白云化流体)在重力和密度差的驱使下向下部生物钻孔发育的泥晶灰岩进行回流渗透白云化,最终形成豹斑云岩的过程中的成岩环境、成岩作用及其产物的组合。这类成岩亚相中,豹斑白云岩具有代表性,其中白云石的盐度为2.47%~2.52%,平均为2.49%。X射线衍射分析结果表明,豹斑中白云石为他形—半自形晶,其有序度变化范围较大,主要为0.56~1.00。储集空间以隐伏回流渗透白云化作用形成的白云石晶间孔以及未被填充的生物钻孔为主,其孔隙连通性相比活跃回流渗透白云化作用所形成的晶间孔较差一些,储集性能也相对较弱。这类成岩亚相最典型的岩石结构及构造为豹斑构造(由豹斑和基质构成,其中豹斑成分为白云石,基质成分为方解石)[图 2(i)]和一些生物活动所留下的扰动构造,这也是判别这类成岩亚相的重要岩石学标志。这类成岩亚相纵向上通常分布在比浅埋藏活跃回流渗透云化亚相埋藏更深的部位,且形成的白云岩的厚度相对更薄,是研究区内进入浅埋藏阶段后的另一类重要的储层成岩相。

2.3 胶结相

此类成岩相以钙质胶结相为代表,其形成主要包括2种情况:①早期近地表环境下,泥—微晶白云岩中的孔隙被钙质胶结,在后期淡水环境下孔隙中的钙质未被溶解,从而形成白云岩孔隙中的钙质胶结现象;②白云岩在早期所形成的微裂缝中,后期被钙质所胶结,表生期也未被淡水溶解,形成白云岩微裂缝中的钙质胶结。钙质胶结相主要发育在马五1—马五4亚段,非研究区的优势成岩相。

2.4 充填相

充填相主要为泥质,形成于表生期岩溶环境。表生期岩溶作用会导致大量的岩石破裂、坍塌、溶蚀,与此同时,大气淡水会携带大量的陆源泥质灌入到所形成的溶蚀孔洞及破裂的裂缝中,从而形成白云岩溶孔、溶洞及裂缝中的泥质充填,该成岩相同样主要发育在马五1—马五4亚段。

3 成岩相测井识别 3.1 成岩相测井响应特征

受限于取心成本以及岩心资料未必连续,因此要实现成岩相的单井连续划分以及平面拓展,就必须依托测井资料建立岩心与测井资料之间的对应关系[17-18]。在岩心归位的基础上,把岩心及薄片分析所划分出的成岩相与测井资料进行对比,研究成岩相在测井曲线上的响应特征,并建立了各类白云岩成岩相的测井识别模式。已有研究表明,常规测井曲线中能够灵敏指示成岩相特征变化的曲线包括自然伽马(GR)、声波时差(AC)、补偿中子(CNL)、补偿密度(DEN)、电阻率(RLLD,RLLS)以及光电吸收截面指数(Pe[17-22]等7条测井曲线。

(1)早期大气淡水溶蚀亚相。其岩性为膏、盐溶蚀泥—微晶白云岩及泥—粉晶白云岩,由于其泥质含量及放射性矿物含量较低,在测井曲线上GR表现为相对低值;由于膏、盐溶解作用所形成的铸模孔隙会降低岩石的相对密度、增加岩石的孔隙度,在对应层段的测井曲线上DEN表现为低值,CNLAC均为相对高值(尤其是含气层段);由于岩石中含有灰质组分,这类成岩亚相对应层段的Pe值相对较高(图 3)。

下载eps/tif图 图 3 溶蚀相测井响应(FG7井) Fig. 3 Logging response of dissolution diagenetic facies

(2)表生期大气淡水溶蚀亚相。其岩性为岩溶角砾泥—粉晶白云岩,角砾之间常含岩溶坍塌所携带的陆源泥质,因此,在其对应层段测井曲线上GR表现为高值;由于岩溶作用会形成一定的溶蚀孔隙,从而造成岩石密度的降低以及孔隙度的相对增加,表现在对应层段的测井曲线上为DEN相对减小,CNLAC有所增加,但增加幅度相对早期大气淡水溶蚀亚相相对较小。由于这类成岩亚相内的岩石较少或不含灰质、石盐及硬石膏组分,因此其Pe值相对较低(图 3)。

(3)浅埋藏活跃回流渗透云化亚相。该亚相以残余砂屑粉—细晶白云岩为主,富含晶间孔,因此在测井曲线上往往表现为低GR值、中等DEN值和低—中等CNL值,电阻率值一般根据含气情况而变化,通常表现为相对较高值,且RLLDRLLS具有一定的幅度差,Pe值相对较低(图 4)。

下载eps/tif图 图 4 浅埋藏活跃回流渗透云化亚相测井响应(D48井) Fig. 4 Logging response of shallow buried active reflux seepage dolomitization diagenetic subfacies

(4)浅埋藏隐伏回流渗透云化亚相。该亚相以豹斑白云岩为主,因此,在测井曲线上GR也表现为低值。由于其白云石化程度没有活跃回流渗透白云石化作用彻底,所以DEN为中等—高值,CNL值为低—中等,电阻率值比前2类溶蚀亚相要高,但相比活跃回流渗透云化亚相要低,且RLLDRLLS幅度差不明显,Pe值相对较低,但高于浅埋藏活跃回流渗透云化亚相(图 5)。

下载eps/tif图 图 5 浅埋藏隐伏回流渗透云化亚相测井响应(D78井) Fig. 5 Logging response of shallow buried latent reflux seepage dolomitization diagenetic subfacies

(5)钙质胶结相。该相以泥—微晶白云岩孔隙、微裂缝中的钙质胶结物为代表,因此其测井响应特征为“三高三低”,即PeDENRLLD值较高,而GRCNLAC值均相对较低,此外,RLLDRLLS之间通常具有一定的幅差。

(6)泥质充填相。泥质充填相在测井上的响应特征为“两高四低”,由于此类成岩相泥质组分含量最高,因此具有最高的GR值;泥质组分一般较致密的白云岩相对疏松一些,因此其DEN值相对较低,AC值相对较高;该成岩相一般不含天然气组分,因此其CNLRLLD值也较低。此外,泥质充填相也具有相对较低的Pe值。

3.2 成岩相定量表征与识别 3.2.1 测井曲线交会图分析及定量识别

选取研究区近90个样品,采用对成岩相敏感度相对较高的GRACDENCNLRLLDRLLS以及Pe等测井曲线的交会分析,可以较为直观、清晰地区分6类成岩相之间的分界(图 6)。虽然不同成岩相分布区域可能会有一定程度的重叠,但以DENPeRLLD-PeRLLD-AC交会图为例,能够较为清晰地区分6类成岩相。由于DEN能够反映储层的总孔隙度,而Pe能够有效反映碳酸盐岩的岩性。一般情况下,早期大气淡水溶蚀亚相主要发育在膏盐溶蚀泥—粉晶白云岩,由于早期大气淡水溶蚀石膏及石盐组分形成铸模孔后,一部分可能会被后期的钙质(方解石)胶结填充,因此其Pe相对较高,而密度值中等;表生期大气淡水溶蚀亚相主要发育在岩溶角砾泥—粉晶白云岩中,岩溶作用能够形成大量的溶蚀孔洞,增加岩石孔隙度,但随之而来的泥质灌入也会填充孔隙,因此其Pe相对早期大气淡水溶蚀亚相要低,其密度变化范围相对较大;浅埋藏活跃回流渗透云化亚相主要发育在残余砂屑粉—细晶白云岩中,富含大量晶间孔,白云化程度相对彻底,因此其DENPe均相对较低,而浅埋藏隐伏回流渗透白云化亚相主要发育在豹斑白云岩中,其相比活跃回流渗透白云化亚相的云化程度弱一些,灰质组分相对更多,孔隙发育相对较差,因此具有比活跃回流渗透云化亚相更高的DEN以及相对更高的Pe;钙质胶结相在几种成岩相中Pe值最大,由于孔隙度很小,因此DEN值也相对较高;泥质充填相测井响应特征也比较明显,具有最低的DEN值以及相对较低的Pe

下载eps/tif图 图 6 常规测井曲线交会图 Fig. 6 Crossplot of conventional logging curves

在其他油田中可能含有较多的开发井,因而缺少Pe曲线等测井资料,加之DENCNLPe等测井曲线的探测范围相对较小,易受井径的影响,因此就需要借助更常用且探测范围相对较大的测井曲线资料来区分不同的白云岩成岩相类型。RLLD-ACRLLD-GR交会图也可以较为清晰地区别几类成岩相(图 6)。以RLLD-AC交会图为例,不同矿物或岩石的电阻率大小均不同,通常黏土矿物或泥质的电阻率最小,白云岩的电阻率小于灰岩的电阻率,方解石和硬石膏的电阻率相对较大。因此,在6类成岩相中,钙质胶结相的RLLD值最大,泥质充填相最小,2种云化亚相中由于岩石含有一定的灰质组分,从而其RLLD值要大于2类溶蚀亚相,表生期大气淡水溶蚀亚相中可能由于岩石所含泥质含量相对高于早期大气淡水溶蚀亚相,因此其RLLD值相对更低。

2类溶蚀亚相的白云岩中的溶蚀孔洞相对发育,导致孔隙度增加,AC相对较大,而2类白云化亚相的白云岩中晶间孔隙可能遭受灰质组分的填充,因此岩石相对致密,AC变小,钙质胶结相由于孔隙度也相对较低,因此具有相对较小的AC值,而泥质充填相可能由于泥质组分相对于白云岩较疏松,因此AC值在几类成岩相中处于最大。此外,统计各类成岩相的测井响应特征值见表 1

下载CSV 表 1 伊陕斜坡马五段白云岩储层成岩相的测井响应 Table 1 Logging response characteristics of diagenetic facies of dolomite reservoir of Ma 5 member in Yishan slope
3.2.2 测井定量识别成岩相的方法

根据上述交会图特征,采用Forward软件编写了成岩相定量划分程序,利用方程对交会图中白云岩成岩相区域进行划分,进而建立成岩相测井交会识别图版(图 7)。以DEN-PeRLLD-PeRLLD-AC交会图为例,3种识别方法的具体步骤分别如[图 8(a)(c)]所示。

下载eps/tif图 图 7 伊陕斜坡马五段白云岩储层成岩相测井识别图版 Fig. 7 Logging identification chart of diagenetic facies of dolomite reservoir of Ma 5 member in Yishan slope
下载eps/tif图 图 8 伊陕斜坡马五段白云岩储层成岩相划分流程 Fig. 8 Flow chart of quantitative division of diagenetic facies of dolomite reservoir of Ma 5 member in Yishan slope
3.2.3 成岩相单井识别及结果验证

用上述方法对Y101井进行了单井白云岩成岩相纵向划分。运用PeDEN交会图法识别出的成岩相结果与其他4条测井曲线(GRACCNLRLLD/RLLS)的形态变化均具有较好的匹配关系(图 9),也印证了Pe-DEN交会图法的可行性。早期大气淡水溶蚀亚相和表生期大气淡水溶蚀亚相分布于马五1—马五4亚段,浅埋藏活跃回流渗透云化亚相和浅埋藏隐伏回流渗透云化亚相则分布于马五5亚段,钙质胶结相主要分布于马五1—马五4亚段,泥质充填相主要分布于马五1—马五3亚段。为了进一步验证该方法的有效性,将测井解释与薄片分析的结果进行对比,可以看出,马五1亚段底部(3 249.6~3251.0 m)层段测井解释成岩相结果与薄片鉴定识别结果不符,原因是该段发生了扩径,导致Pe和DEN测井响应失真,此时就可以利用受井径影响小、探测范围相对较大的RLLDAC曲线交会图进行识别。总体上,利用薄片资料所确定的白云岩成岩相与测井解释结果的符合率可达92.85%,而测井解释的成岩相在纵向上与实测的岩心孔隙度也具有良好的对应关系,其中早期大气淡水溶蚀亚相储层物性最好,孔隙度为1.27%~8.40%,平均为4.05%;浅埋藏活跃回流渗透云化亚相次之,孔隙度为1.14%~5.83%,平均为3.14%;表生期大气淡水溶蚀亚相相对较差,孔隙度为0.36%~5.55%,平均为2.59%, 而浅埋藏隐伏回流渗透云化亚相更差,钙质胶结相和泥质充填相物性均差。此外,利用测井解释的成岩相与利用CNL-DEN曲线综合识别的岩性剖面在纵向上也具有较好的匹配性,即早期大气淡水溶蚀亚相发育在马五1—马五4亚段的泥晶白云岩中,表生期大气淡水溶蚀亚相发育在含泥白云岩或泥云岩中,浅埋藏活跃回流渗透云化亚相多集中发育在马五5亚段的粉—细晶白云岩中,而浅埋藏隐伏回流渗透云化亚相发育在豹斑云岩中,钙质胶结相以孔隙、微裂缝中的钙质充填物为主,泥质充填相以溶蚀孔、洞中的泥质最为典型。上述对比结果均表明,利用Pe-DEN交会图定量识别低渗—特低渗白云岩储层成岩相具有可靠性。在缺乏PeDEN测井资料或利用Pe-DEN交会图定量识别白云岩储层成岩相存在误差的情况下,可以用RLLD-AC交会替代以进一步提高测井解释成岩相的精度。

下载eps/tif图 图 9 Y101井白云岩储层成岩相测井解释及结果验证 Fig. 9 Logging interpretation and result verification of diagenetic facies of dolomite reservoir in well Y101
4 成岩相空间展布 4.1 成岩相剖面展布

利用单井识别白云岩成岩相的方法对伊陕斜坡近百口井的马五段进行了白云岩成岩相的纵向定量划分,并建立连井剖面(图 10)。可以看出:早期大气淡水溶蚀亚相和表生期大气淡水溶蚀亚相分布于马五1—马五4亚段,其中FG7井的马五1亚段下部、马五2亚段以及马五3亚段上部以早期大气淡水溶蚀亚相为主,而马五3亚段中下部及马五4亚段主要发育表生期大气淡水溶蚀亚相,马五5亚段则以浅埋藏活跃回流渗透白云化亚相为主,仅在顶部发育浅埋藏隐伏回流渗透白云化亚相;Y844井马五1—马五4亚段整体以早期大气淡水溶蚀亚相为主,而表生期大气淡水溶蚀亚相仅分布于马五1亚段中部、马五3亚段上部以及马五4亚段中部,马五5亚段演变为以浅埋藏隐伏回流渗透白云化亚相为主;向北至Y101井,马五1亚段顶部、马五3亚段以及马五4亚段中部演变为以表生期大气淡水溶蚀亚相为主,马五5亚段则发育浅埋藏活跃回流渗透白云化亚相;S153井的马五2、马五3亚段以及D48井的马五2亚段顶部和底部、马五3亚段以及马五4亚段表生期大气淡水溶蚀亚相发育,但D48井的单层厚度相对较小,而马五5亚段2口井均以浅埋藏活跃回流渗透云化亚相为主。整体上,马五1亚段主要发育早期大气淡水溶蚀亚相,其底部发育泥质充填相;马五2、马五3亚段在中部、北部及马五3亚段在南部较发育表生期大气淡水溶蚀亚相,且在Y101井区马五2亚段顶部发育钙质胶结相;马五4亚段在南部、北部发育表生期大气淡水溶蚀亚相,在中部主要发育早期大气淡水溶蚀亚相及钙质胶结相;马五5亚段在中南部发育浅埋藏隐伏回流渗透云化亚相,在全区主要发育浅埋藏活跃回流渗透云化亚相。

下载eps/tif图 图 10 伊陕斜坡马五段白云岩储层成岩相纵向分布 Fig. 10 Vertical distribution characteristics of diagenetic facies of dolomite reservoir of Ma 5 member in Yishan slope
4.2 成岩相平面展布及有利储层分布预测

(1)成岩相平面展布特征。在白云岩成岩相单井分析、连井剖面对比的基础上,绘制了区内各层段的白云岩成岩相平面分布图。分别以马五1亚段和马五5亚段为例,可以看出,马五1亚段表生期大气淡水溶蚀亚相的岩石为泥云岩、含泥含灰云岩以及含灰泥云岩,分布于区内西北部靠近剥蚀区的D48井西部—S153井—S257井区、西部S123井— S124井—S245井区、东南部Y61井—延川—延长— Y104井一带以及南部Y558井—富县—Yt1井一带。早期大气淡水溶蚀亚相岩石以泥—粉晶白云岩为主,分布范围相对较广,研究区西北部、东北部、中部以及西南部,连片性较好[图 11(a)]。

下载eps/tif图 图 11 伊陕斜坡中东部马五段白云岩储层成岩相图 Fig. 11 Diagenetic facies of dolomite reservoir of Ma 5 member in central-eastern Yishan slope

马五5亚段主要发育浅埋藏活跃回流渗透云化亚相和浅埋藏隐伏回流渗透云化亚相。其中,研究区北部和中部广大地区以及南部Y112井—Y135井—Y106井一线以南主要为浅埋藏活跃回流渗透云化亚相,岩石类型以含云灰岩、云灰岩以及含灰白云岩为主;浅埋藏隐伏回流渗透云化亚相分布面积相对较小,发育在区内南部Y558井—延安— Y127井—延长—Y130井一带[图 11(b)]。

(2)有利储层分布预测。利用成岩相平面展布特征与测井解释的孔隙度等值线平面叠合图对区内有利储层的分布区进行了初步预测。可以看出,马五1亚段优势成岩相类型主要为早期大气淡水溶蚀亚相和表生期大气淡水溶蚀亚相,其中早期大气淡水溶蚀亚相中储层分布面积相对较广,多呈连片式发育,该亚相储层物性较好,孔隙度通常为2.0%~4.0%;表生期大气淡水溶蚀亚相中储层分布面积相对较小,主要发育在横山地区、S124井区以及研究区的东南部,储层物性相对较差,孔隙度一般为1.5%~3.0%。马五5亚段储层的优势成岩相为浅埋藏活跃回流渗透云化亚相和浅埋藏隐伏回流渗透云化亚相,储层主要发育在浅埋藏活跃回流渗透云化亚相,呈孤立状分布于横山、S245井区、安塞西部、延川西部以及研究区南部地区,该亚相内储层物性也相对较好,孔隙度为1.0%~3.5%;浅埋藏隐伏回流渗透云化亚相储层分布面积相对较小,物性较差,孔隙度为1.0%~2.0%。

5 结论

(1)岩心观察及镜下薄片鉴定认为,伊陕斜坡中东部马五5—马五1亚段低渗—特低渗白云岩储层成岩相可划分为早期大气淡水溶蚀亚相、表生期大气淡水溶蚀亚相、浅埋藏活跃回流渗透云化亚相、浅埋藏隐伏回流渗透云化亚相、钙质胶结相和泥质充填相等6类。

(2)以Pe-DENRLLD-Pe交会为代表,配合RLLD-AC交会的测井参数交会图法及Forward软件处理划分程序综合应用是进行伊陕斜坡马五5—马五1亚段低渗—特低渗白云岩储层不同成岩相单井纵向上定量识别与连续划分的有效方法。

(3)早期大气淡水溶蚀亚相分布于伊陕斜坡中部及南部的马五1—马五2亚段以及中部的马五4亚段,而表生期大气淡水溶蚀亚相分布于南部和北部的马五3—马五4亚段以及中部的马五3亚段;浅埋藏活跃回流渗透云化亚相分布于研究区最南部、中部及北部地区的马五5亚段,而浅埋藏隐伏回流渗透云化亚相多分布在区内南部的马五5亚段;钙质胶结相主要分布于马五2亚段中、上部及马五4亚段中部;泥质充填相主要分布于马五1亚段底部。

(4)伊陕斜坡马五1亚段早期大气淡水溶蚀亚相相对更发育,分布于区内西北部、东北部、中部及西南部地区,而表生期大气淡水溶蚀亚相分布在区内靠近剥蚀区的北部D48井西部—横山地区、西部S123—S245井一带、南部Y558井—富县地区以及清涧—延川—延长—Y104井一带;马五5亚段区内除南部Y558井—延安—延长一带为浅埋藏隐伏回流渗透云化亚相发育区,其余广大地区则以浅埋藏活跃回流渗透云化亚相为主,有利储层主要分布于早期大气淡水溶蚀亚相和浅埋藏活跃回流渗透云化亚相,其次为表生期大气淡水溶蚀亚相,浅埋藏隐伏回流渗透云化亚相储层分布最少。

参考文献
[1]
杨威, 魏国齐, 金惠, 等. 碳酸盐岩成岩相研究方法及其应用:以扬子地块北缘飞仙关组鲕滩储层为例. 岩石学报, 2011, 27(3): 749-756.
YANG W, WEI G Q, JIN H, et al. Methods for diagenetic facies research on carbonate rocks and its application:Example from oolitic beach reservoir of the Feixianguan Formation in the northern margin of Yangtze block. Acta Petrologica Sinica, 2011, 27(3): 749-756.
[2]
王起琮, 赵淑萍, 魏钦廉, 等. 鄂尔多斯盆地中奥陶统马家沟组海相碳酸盐岩储集层特征. 古地理学报, 2012, 14(2): 229-242.
WANG Q C, ZHAO S P, WEI Q L, et al. Marine carbonate reservoir characteristics of the Middle Ordovician Majiagou Formation in Ordos Basin. Journal of Paleogeography, 2012, 14(2): 229-242.
[3]
李百强, 王起琮, 张小莉, 等. 鄂尔多斯盆地中南部马家沟组马五5-马五1亚段白云岩成岩相特征及其对储层的影响. 西北大学学报(自然科学版), 2018, 48(2): 255-260.
LI B Q, WANG Q C, ZHANG X L, et al. Characteristics of diagenetic facies and its effect on the reservoir of M55-M51 submembers of Majiagou Formation in central-southern Ordos Basin. Journal of Northwest University(Natural Science Edition), 2018, 48(2): 255-260.
[4]
王小芬, 杨欣. 鄂尔多斯盆地富县地区马五段碳酸盐岩成岩作用研究. 岩性油气藏, 2011, 23(3): 75-79.
WANG X F, YANG X. Carbonate diagenesis of the fifth member of Majiagou Formation in Fuxian area, Ordos Basin. Lithologic Reservoirs, 2011, 23(3): 75-79. DOI:10.3969/j.issn.1673-8926.2011.03.014
[5]
陈彦华, 刘莺. 成岩相:储集体预测的新途径. 石油实验地质, 1994, 16(3): 274-281.
CHEN Y H, LIU Y. Diagenetic facies:a new approach to the predication of reservoir rocks. Petroleum Geology & Experiment, 1994, 16(3): 274-281.
[6]
张秀莲. 碳酸盐岩成岩研究的新方法:成岩综合柱状图. 地质科学, 1995(4): 401-413.
ZHANG X L. A new method in the study of carbonatite diagenesis:Composite column of diagenesis. Chinese Journal of Geology(Scientia Geologica Sinica), 1995(4): 401-413.
[7]
RAILSBACK L B. Carbonate diagenetic facies in the Upper Pennsylvanian Dennis Formation in Iowa, Missouri, and Kansas. Journal of Sedimentary Petrology, 1984, 54(3): 986-999. DOI:10.1306/212F8549-2B24-11D7-8648000102C1865D
[8]
范宜仁, 李菲, 邓少贵, 等. 致密砂岩储层成岩相特征分析及测井识别. 测井技术, 2018, 42(3): 63-70.
FAN Y R, LI F, DENG S G, et al. Characteristics analysis of diagenetic facies in tight sandstone reservoir and its logging identification. Well Logging Technology, 2018, 42(3): 63-70.
[9]
冉冶, 王贵文, 赖锦, 等. 利用测井交会图法定量表征致密油储层成岩相:以鄂尔多斯盆地华池地区长7致密油储层为例. 沉积学报, 2016, 34(4): 694-706.
RAN Y, WANG G W, LAI J, et al. Quantitative characterization of diagenetic facies of tight sandstone oil reservoir by using logging crossplot:a case study on Chang 7 tight sandstone oil reservoir in Huachi area, Ordos Basin. Acta Sedimentologica Sinica, 2016, 34(4): 694-706.
[10]
赖锦, 王贵文, 王书南, 等. 碎屑岩储层成岩相测井识别方法综述及研究进展. 中南大学学报(自然科学版), 2013, 44(12): 4942-4953.
LAI J, WANG G W, WANG S N, et al. Overview and research progress in logging recognition method of clastic reservoir diagenetic facies. Journal of Central South University(Science and Technology), 2013, 44(12): 4942-4953.
[11]
李百强.延长探区马家沟组马五段碳酸盐岩储层分布规律及主控因素.西安: 西安石油大学, 2016.
LI B Q. Distribution law and main controlling factors of carbonate reservoir in member Ma 5 of Majiagou Formation in Yanchang exploration area. Xi'an: Xi'an Shiyou University, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10705-1016097118.htm
[12]
李维, 张军涛, 朱筱敏, 等. 鄂尔多斯盆地富县地区马五段白云岩储层特征与演化. 岩性油气藏, 2016, 28(3): 58-67.
LI W, ZHANG J T, ZHU X M, et al. Characteristics and evolution of the Ordovician dolomite reservoir of Ma 5 member in Fuxian area, Ordos Basin. Lithologic Reservoirs, 2016, 28(3): 58-67. DOI:10.3969/j.issn.1673-8926.2016.03.009
[13]
李百强, 王起琮, 张小莉, 等. 鄂尔多斯盆地中南部马五1-2亚段岩相古地理特征及其对储层的影响. 沉积学报, 2019, 37(3): 589-600.
LI B Q, WANG Q C, ZHANG X L, et al. Characteristics of lithofacies paleography and its effect on the Majiagou submember 51-2 reservoir in the central-southern Ordos Basin. Acta Sedimentologica Sinica, 2019, 37(3): 589-600.
[14]
MOU D C, BRENNER R L. Control of reservoir properties of Tensleep sandstone by depositional and diagenetic facies:Lost Soldire field, Wyoming. Journal of sedimentology, 1982, 39: 1119-1136.
[15]
SPTL C, WRIGHT V P. Groundwater dolocretes from the Upper Triassic of the Paris Basin, France:a case study of an arid, continental diagenetic facies. Sedimentology, 1992, 39: 1119-1136. DOI:10.1111/j.1365-3091.1992.tb02000.x
[16]
ANDREWS J E. Geochemical indicators of depositional and early diagenetic facies in Holocene carbonate muds and their preservation potential during stabilization. Chemical Geology, 1991, 93: 267-289. DOI:10.1016/0009-2541(91)90118-B
[17]
邹才能, 陶士振, 周慧, 等. 成岩相的形成、分类与定量评价方法. 石油勘探与开发, 2008, 35(5): 526-540.
ZOU C N, TAO S Z, ZHOU H, et al. Genesis, classification and evaluation method of diagenetic facies. Petroleum Exploration and Development, 2008, 35(5): 526-540. DOI:10.3321/j.issn:1000-0747.2008.05.002
[18]
李国欣, 赵太平, 石玉江, 等. 鄂尔多斯盆地马家沟组碳酸盐岩储层成岩相测井识别评价. 石油学报, 2018, 39(10): 1141-1154.
LI G X, ZHAO T P, SHI Y J, et al. Diagenetic facies logging recognition and evaluation of carbonate reservoirs in Majiagou Formation, Ordos Basin. Acta Petrolei Sinica, 2018, 39(10): 1141-1154. DOI:10.7623/syxb201810006
[19]
王正允, 张紫光, 陈梦蛟, 等. 冀北坳陷雾迷山组的成岩作用及成岩相特征:以辽宁省凌源县魏杖子剖面为例. 石油天然气学报, 2009, 31(4): 1-8.
WANG Z Y, ZHANG Z G, CHEN M J, et al. Diagenesis and features of diagenetic facies of Wumishan Formation of Jibei Depression:by taking the Weizhangzi section of Lingyuan county in Liaoning province for example. Journal of Oil and Gas Technology, 2009, 31(4): 1-8. DOI:10.3969/j.issn.1000-9752.2009.04.001
[20]
石玉江, 肖亮, 毛志强, 等. 低渗透砂岩储层成岩相测井识别方法及其地质意义:以鄂尔多斯盆地姬塬地区长8段储层为例. 石油学报, 2011, 32(5): 820-828.
SHI Y J, XIAO L, MAO Z Q, et al. An identification method for diagenetic facies with well logs and its geological significance in low-permeability sandstones:a case study on Chang 8 reservoirs in the Jiyuan region, Ordos Basin. Acta Petrolei Sinica, 2011, 32(5): 820-828.
[21]
张海涛, 时卓, 石玉江, 等. 低渗透致密砂岩储层成岩相类型及测井识别方法:以鄂尔多斯盆地苏里格气田下石盒子组8段为例. 石油与天然气地质, 2012, 33(2): 256-264.
ZHANG H T, SHI Z, SHI Y J, et al. Diagenetic facies types and logging identification methods for low-permeability tight sandstone reservoirs:a case study on the 8 th member of Xiashihezi Formation in Sulige Gas Field, Ordos Basin. Oil & Gas Geology, 2012, 33(2): 256-264.
[22]
赖锦, 王贵文, 信毅, 等. 库车坳陷巴什基奇克组致密砂岩气储层成岩相分析. 天然气地球科学, 2014, 25(7): 1019-1032.
LAI J, WANG G W, XIN Y, et al. Diagenetic facies analysis of tight sandstone gas reservoir of Bashijiqike Formation in Kuqa Depression. Natural Gas Geoscience, 2014, 25(7): 1019-1032.