E-mail Alert Rss
 

物探与化探, 2022, 46(6): 1369-1380 doi: 10.11720/wtyht.2022.1416

地质调查·资源勘查

X凹陷A构造低阻气层成因机理分析

程仁杰,1, 孙建孟1, 刘建新2, 迟蓬1, 吕馨頔1, 胡文亮2, 付焱鑫2, 赵文兵3

1.中国石油大学(华东) 地球科学与技术学院,山东 青岛 266580

2.中海石油(中国)有限公司 上海分公司,上海 200335

3.中国石油塔里木油田分公司勘探事业部,新疆 库尔勒 841000

Genetic mechanisms of low-resistivity gas zones in structure A of sag X

CHENG Ren-Jie,1, SUN Jian-Meng1, LIU Jian-Xin2, CHI Peng1, Lyu Xin-Di1, HU Wen-Liang2, FU Yan-Xin2, ZHAO Wen-BING3

1. School of Geoscience, China University of Petroleum (East China), Qingdao 266580, China

2. CNOOC China Limited, Shanghai Branch, Shanghai 200335, China

3. Research Institute of Petroleum Exploration and Development, Tarim Oilfield Company, PetroChina, Korla 841000, China

责任编辑: 王萌

收稿日期: 2021-08-2   修回日期: 2022-01-28  

基金资助: 青岛海洋科学与技术试点国家实验室山东省专项经费(2021QNLM20001)
国家自然科学基金(42174143)
国家自然科学基金(41874138)

Received: 2021-08-2   Revised: 2022-01-28  

作者简介 About authors

程仁杰(1998-),男,重庆云阳人,现为中国石油大学(华东)在读硕士生,主要研究方向为测井数据处理与综合解释。Email:582476598@qq.com

摘要

经勘探开发证实,X凹陷A构造H组地层Q3c上部层段存在大量低阻气层。针对该层段地质沉积环境认识不清、储层微观认识不深入、气层低阻成因尚未明确等问题,以研究区3口井的测井资料为基础,结合钻井、录井资料和大量岩石物理实验资料开展相关研究。对研究区开展了基于薄片鉴定资料的岩石学、物性特征分析;通过连井剖面、特殊测井资料并结合大量岩石物理实验对低阻气层成因机理开展研究;基于数字岩心技术构建多组分导电模型,从微观可视化尺度证实了低阻气层成因机理,开展有限元电性模拟定量分析各低阻成因对电阻率降低的贡献。研究结果表明,由于高阳离子交换容量粘土矿物的存在以及良好物性基础上发育复杂孔隙结构,共同导致了研究区Q3c上部气层的低阻响应。其中,粘土附加导电性对低阻响应的贡献为35.63%,良好物性条件下的复杂孔隙结构对低阻响应的贡献达到64.37%,电性模拟结果与测井电性特征吻合,证实了Q3c上部低阻气层成因机理。

关键词: 低阻气层; 粘土附加导电; 孔隙结构; 数字岩心; 导电模拟

Abstract

As confirmed by exploration and development, many low-resistivity gas zones exist in the upper portion of layer Q3c of the H Formation in structure A of Sag X. Given the problems such as unclear understanding of the geological sedimentary environment, insufficient microcosmic knowledge about the reservoirs, and unclear causes of the low resistivity of the gas zones, this study conducted extensive research based on the log data of three wells in the study area, as well as the data on drilling and many petrophysical experiments. Specifically, the petrological and physical property characteristics of the study area were studied using the thin-section identification data; the genetic mechanisms of low-resistivity gas zones were studied using the well tie sections and the special log data, as well as the data of many petrophysical experiments; the formation mechanisms of low-resistivity gas zones were confirmed from the microscopic visualization scale by constructing a multi-component conductivity model using the digital core technique, and the contributions of various low-resistivity geneses to the decrease in resistivity were quantitatively analyzed through the finite element-based electrical simulations. As indicated by the study results, the low-resistivity response of the gas zones in the study area is caused by the presence of clay minerals with high clay content and high cation exchange capacity and also results from the complex pore structure formed under the favorable physical property conditions in the anomalous high-pressure depositional setting. The contributions of the clay additional conductivity and the complex pore structure to the low resistivity are 35.63% and up to 64.37%, respectively. The electrical simulation results are consistent with the log-derived electrical characteristics, verifying the genetic mechanisms of the low-resistivity gas zones in the upper portion of the Q3c.

Keywords: low-resistivity gas zone; clay additional conductivity; pore structure; digital core; electrical simulation

PDF (10296KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

程仁杰, 孙建孟, 刘建新, 迟蓬, 吕馨頔, 胡文亮, 付焱鑫, 赵文兵. X凹陷A构造低阻气层成因机理分析[J]. 物探与化探, 2022, 46(6): 1369-1380 doi:10.11720/wtyht.2022.1416

CHENG Ren-Jie, SUN Jian-Meng, LIU Jian-Xin, CHI Peng, Lyu Xin-Di, HU Wen-Liang, FU Yan-Xin, ZHAO Wen-BING. Genetic mechanisms of low-resistivity gas zones in structure A of sag X[J]. Geophysical and Geochemical Exploration, 2022, 46(6): 1369-1380 doi:10.11720/wtyht.2022.1416

0 引言

X凹陷位于东海陆架盆地东北部,是在前古近系基底上经历了断陷、拗陷与区域沉降3个大的演化阶段发育而来的具有东断西超特征的大型凹陷,是中国近海油气资源丰富且具有较大勘探开发潜力的新近系含油气凹陷之一[1-3]。X凹陷整体呈NNE方向展布,具有东西分带、南北分块的特点,自西向东分为斜坡带、西次凹、中央反转带、东次凹及东部断阶带5个次级构造带。

A构造位于中央反转构造带,H组为该气田低渗致密砂岩天然气勘探开发的主要目的层,H组沉积期同时存在3个方向的物源供给,为一套中—大型辫状河三角洲前缘沉积[4-5]。H组碎屑沉积物受搬运距离的影响,纵向上,自下而上沉积碎屑成熟度逐渐降低。研究区主要产气层段埋深均大于3 400 m,地层受到了强烈的压实作用,加上成岩过程中两次幕式抬升导致成岩环境的多次改变,储层具有物性差、强非均质性、压力系统异常、孔隙结构复杂等特征。在勘探开发过程中发现了诸多测井复杂响应层段,如低阻气层、层内或相邻层间相对低阻气层、高阻干层等。A构造3口井的同一地质层系Q3c上部均发现了低阻气层,针对A-2井的低阻气层,多位学者开展研究得到不同的结论[6-7],有学者认为导电矿物、粘土附加导电性和高矿化度泥浆侵入是气层低阻的成因,另外有学者认为粘土附加导电性和高不动水饱和度导致了气层的低阻响应。针对相同沉积背景下研究区内同一地质层系出现的低阻气层存在沉积背景不清楚、成因机理尚未明确、低阻气层微观成因机理以及导电机理认识不足等问题,本文以丰富的测井资料为基础,结合大量的岩石物理实验资料,对比分析研究区Q3c上部低阻气层的成因机理,基于矿物扫描、CT扫描等可视化手段建立三维数字岩心对低阻气层导电机理进行定量化研究。

1 储层特征

1.1 研究区岩石学特征

研究区A构造位于中央反转构造带中北部,其H组地层沉积期间来源于北部的母岩区碎屑物质顺着自北向南狭长河道长距离搬运,形成了H组砂岩成熟度由浅到深逐渐增大的空间展布[8]。统计了A构造3口井共计554颗岩心薄片鉴定资料(表1),研究区H组砂岩石英含量为63%~65%、长石含量为15%~18%、岩屑含量为18%~20%,具有石英含量高、杂基和胶结物含量低的特征,按照传统命名规则,研究区砂岩以长石岩屑质石英砂岩为主。纵向上Q6到Q3稳定矿物组分含量逐渐降低,碎屑粒度呈现变粗的趋势,分选和磨圆程度逐渐变差,反映随着搬运距离变短,碎屑物质成分成熟度和结构成熟度逐渐降低。

表1   A构造3口井薄片鉴定分析统计

Table 1  Analysis and statistics of thin sections of 3 wells in structure A

层位岩心数量矿物含量均值/%最大粒度
均值/mm
主要粒径
分布/mm
分选性磨圆度
石英钾长石斜长石岩屑
Q318363.06.289.3520.11.550.30~0.58中—好次棱—次圆
Q429265.16.669.6618.60.830.20~0.48好—中次棱—次圆
Q5664.26.5011.0018.30.570.15~0.36次棱—次圆
Q65864.56.169.7819.60.770.15~0.40次棱—次圆

新窗口打开| 下载CSV


1.2 研究区物性特征

研究区Q3、Q4层段为凝析气勘探开发主要目的层段,对A-2井Q3、Q4层段的岩心物性分析资料进行统计(表2),统计结果表明目的层段物性整体较差,属于低孔、低渗—特低渗储层。Q3c层间物性差异较大,Q3c上部地层孔隙度均值比下部大3.5%、渗透率均值比下部大0.9 mD。观察Q3层段岩心铸体薄片(图1),发现Q3层段碎屑颗粒之间多呈线—点接触关系,反映成岩过程中较强的压实作用,地层原生粒间孔隙在上覆地层强压实作用下大大减少,孔隙类型以溶蚀粒间孔和粒内溶孔为主,Q3c上部3 607 m外的薄片仍观察到一定量的原生粒间孔隙。有学者统计发现,H组储层次生溶孔对总孔隙度的贡献达到甚至超过50%[9]。长石及易溶岩屑等矿物组分在酸性成岩环境下易发生蚀变,溶蚀产生的溶蚀粒间孔、粒内溶孔和铸模孔对储层物性起到积极的改造作用。

表2   A-2井岩心物性分析资料统计

Table 2  Analysis of core porosity and permeability of A-2 well

层段埋深/m样品数孔隙度分布/%孔隙度均值/%渗透率均值/mD
Q3b3429~35502948.1~9.08.41.50
Q3c上部3592~360868.5~19.314.52.54
Q3c下部3611~3659148.3~12.011.01.63
Q4b3810~38682338.6~12.19.7290.911

新窗口打开| 下载CSV


图1

图1   Q3c层段铸体薄片(10倍目镜×2倍物镜,单偏光)

a—A-2井Q3c上部3 600 m;b—A-2井Q3c上部3 607 m;c—A-2井Q3c下部3 634 m

Fig.1   Casting thin sections of layer Q3c(20 times eyepiece×2 times lbjective lens,single polarized light)

a—well A-2, the upper part of the Q3c,3 600 m; b—well A-2, the upper part of the Q3c, 3 607 m; c—well A-2, the lower part of the Q3c,3 634 m


图1a、b两块铸体薄片取自Q3c上部低阻气层,对比图1c下部3 634 m高阻气层的铸体薄片,发现低阻气层薄片整体更“脏”,泥质填隙物更多,且以条带状和结核状富集于岩石颗粒之间,且通过图1b观察到低阻气层局部保留了一定量的原生粒间孔,说明Q3c上部低阻气层的岩性和孔隙结构与下部高阻气层存在较大差异。

1.3 常规测井响应特征

研究区A构造3口井的连井剖面显示Q3c层段为一大套连续砂体沉积,GR曲线整体呈箱型、近似钟型,大套砂体内夹有多个泥质薄层,底部为一套厚泥岩沉积(图2)。图中每口井第1道为岩性测井曲线道、第2道为孔隙度测井曲线道、第3道为深度道、第4道为地层层系道、第5道为气测录井道、第6道为电阻率测井曲线道,A-2井第7道为随钻电阻率曲线道。3口井的深、中、浅3条电阻率曲线在储层段均呈现明显的侵入特征,且在Q3c上部储层段出现异常的低阻响应,深侧向电阻率仅为6~8 Ω·m,而下部储层电阻率达到30~70 Ω·m,上部储层段电阻率远低于下部储层电阻率,甚至低于泥岩段的电阻率(泥岩段电阻率为20~25 Ω·m)。Q3c储层段三孔隙度曲线反映良好的物性特征,补偿中子曲线(CNCF)与补偿密度曲线(ZDEN)反向刻度显示明显的气层响应特征。气测录井显示全层段含气性好,部分深度段全烃含量达到10%,C1含量占全烃含量90%以上。A-2井Q3c上部3 606.5 m深度点进行了MDT取样,取样流体中C1含量高达93.2%,勘探及开发均证实研究区A构造Q3c上部发育有大量具有潜在商业价值的低阻气层。

图2

图2   A构造3口井连井剖面

Fig.2   The inter-well profiles of 3 wells in the A structure


2 低阻成因机理分析

2.1 低阻成因概述

低阻油气层由于低电阻率特征,在测井解释流程中很容易将其与水层混淆,解释含水饱和度会远远大于地层真实含水饱和度,导致流体性质的判断错误。低阻油气层在全球各大油田均被大量发现,随着相关学者对低阻油气层难题的攻关,目前已经建立起完整的低阻油气层成因机理理论体系。国内有关学者将低阻油气层成因机理分为内因和外因两部分[10]。内因包括:岩性细、泥质含量重、粘土附加导电性以及复杂孔隙结构等;外因则包括:泥浆侵入、砂泥岩薄互层、地层水矿化度差异、导电矿物等。经过调研发现大量低阻油气层成因并非单一因素,而是归结于多种成因的共同作用。基于已建立的低阻成因理论体系,结合测井资料和岩石物理实验数据,开展A构造Q3c上部低阻气层成因机理的研究。

2.2 低阻气层成因机理研究

2.2.1 低阻外部因素排除

A-2井采用随钻测井仪器测得Q3c上部储层电阻率为6~8 Ω·m,不同垂向探测深度的随钻测井电阻率曲线侵入特征不明显,且随钻测井电阻率测量值与电缆测井深侧向电阻率值一致,因此深探测电阻率测井曲线反映地层真实电阻率,排除泥浆侵入导致低阻这一外因。

图3为A-2井Q3c层段高、低阻两个井段的FMI成像图像,其中左道为静态成像图,右道为动态成像图。FMI成像图像显示Q3c层段地层岩性以细砂岩为主,偶见砾石和泥质薄层发育,具有强非均质性特征。Q3c上部3 596~3 602 m低阻气层的FMI动态图像显示存在多个泥质薄层,不存在明显的砂泥岩交互。由常规测井可知Q3c纯泥岩段电阻率值高达20~25 Ω·m,而Q3c上部低阻层段电阻率甚至低于纯泥岩段电阻率,故排除砂泥岩薄互层导致低阻这一外因。

图3

图3   A-2井Q3c层段FMI成像图像

a—3 596~3 602 m低阻气层段FMI成像; b—3 606~3 612 m高阻气层段FMI成像

Fig.3   FMI dynamic image of Q3c formation of Well A-2

a—FMI image of 3 596~3 602 m low-resistance gas interval;b—FMI image of 3 606~3 612 m high resistance gas interval


连井剖面显示3口井Q3c上部低阻气层和下部高阻气层之间不存在明显隔层,A-2井3 606~3 612 m的FMI动态图像显示上部低阻气层和下部高阻气层之间偶见泥质薄层发育,Q3c层段纵向连通性整体较好,上部低阻气层和下部高阻气层为同一套气水系统,故排除由于层间地层水矿化度差异导致低阻这一外因。

统计研究区Q3层段共计183块薄片鉴定资料并未发现黄铁矿等导电矿物存在,故排除导电矿物存在导致低阻这一外因。

2.2.2 高粘土含量及粘土附加导电性

研究区H组地层为一套中—大型辫状河三角洲前缘沉积,Q3段主要发育水下分流河道和分流间湾两类沉积微相。FMI动态图像观察到Q3c砂体交错层理非常发育,且在砂体底部可观察到明显的冲刷面,反映Q3c砂体为一套沉积水动力整体较强的水下分流河道沉积微相。在砂岩沉积间歇期,由于水动力环境弱,细粒沉积物及粘土物质大量沉积发育分流间湾沉积微相,如Q3c底部以及Q3b底部均发育以灰色泥岩、粉砂质泥岩为主的分流间湾沉积微相。图2这3口井的连井剖面显示Q3c砂岩段的GR曲线整体呈“箱型”或“钟型”特征,Q3c上部低阻气层的GR值有变大的趋势,反映沉积水动力环境逐渐变弱,沉积微相类型由水下分流河道向分流间湾逐渐过渡。由于沉积水动力环境较弱,Q3c上部低阻气层砂体相较于下部高阻气层段泥质薄层更发育、粘土矿物含量更高。

通过图1薄片鉴定资料发现,Q3c上部低阻气层岩心铸体薄片相较于下部高阻气层岩心整体更“脏”,部分泥质呈颗粒状富集在孔隙中。对比Q3c上部低阻气层和下部高阻气层各3颗X衍射—全岩实验资料(表3),低阻气层的粘土含量高于高阻气层近10%,石英和钾长石含量低于高阻气层近5%,其他矿物含量无明显差异。X衍射—粘土实验发现高、低阻岩心粘土矿物类型也存在明显差异(图4)。低阻气层粘土矿物类型以伊蒙混层和绿泥石为主,伊蒙混层含量占粘土矿物总量43%,绿泥石含量占粘土矿物总量31.3%,伊蒙混层中蒙脱石含量仅为16.7%;高阻气层粘土矿物以绿泥石和伊利石为主,分别占粘土矿物总量46%和27.7%。通过扫描电镜观察高、低阻气层岩心微观特征(图5),高阻气层岩心颗粒表面及粒间孔内分布大量叶片状绿泥石,低阻气层岩心粒间孔之间中发育大量的絮状、搭桥状伊利石,颗粒表面还分布一定量的叶片状绿泥石。

表3   X衍射—全岩实验结果对比

Table 3  Comparison of X-ray diffraction-whole rock experiment results

井次岩性层位深度/m矿物种类和含量/%
石英钾长石斜长石方解石白云石粘土总量
A-2细砂岩Q3c上部
低阻气层
360062.25.715.91.40.614.2
359660.24.918.31.20.315.1
A-4370959.75.714.11.00.519.0
A-2Q3c下部
高阻气层
364565.69.615.41.81.26.4
362762.410.816.70.30.89.0
A-4373269.010.613.62.01.03.8

新窗口打开| 下载CSV


图4

图4   高、低阻气层各类型粘土矿物相对含量直方图

Fig.4   Histogram of relative content of each clay mineral type


图5

图5   高、低阻气层扫描电镜微观特征对比

a—A-2井高阻气层扫描电镜;b—A-2井低阻气层扫描电镜;c—A-2井低阻气层扫描电镜

Fig.5   Comparison of microscopic features of scanning electron microscope

a—SEM of high-resistance gas layer in well A-2;b—SEM of low-resistance gas layer in well A-2;c—SEM of low-resistance gas layer in well A-2


不同类型粘土矿物具有不同阳离子交换能力,通常用阳离子交换容量CEC进行表征,各类型粘土矿物的阳离子交换容量见表4,其中蒙脱石具有最大的阳离子交换容量,绿泥石的CEC几乎接近于0。Q3c上部低阻气层粘土矿物含量高,且粘土矿物类型以伊蒙混层为主,伊蒙混层阳离子交换容量大,粘土矿物附加导电性大大降低地层电阻率。

表4   各类型粘土矿物阳离子交换容量

Table 4  Cation exchange capacity of various types of clay minerals

粘土矿物类型理论CEC范围/
(mmol·g-1)
CEC均值/
(mmol·g-1)
蒙脱石0.80~1.501.20
伊利石0.10~0.400.25
高岭石0.03~0.250.09
绿泥石接近于00

新窗口打开| 下载CSV


考虑在相同总孔隙度和骨架密度地层条件下,由粘土矿物总量以及粘土矿物类型导致的导电效率差异。根据X衍射实验得到的Q3c上部低阻气层和下部高阻气层粘土矿物总量以及各类型粘土矿物相对含量,结合式(1)计算得到Q3c上部低阻气层阳离子交换容量Qv为0.292 mmol/g,下部高阻气层Qv为0.082 mmol/g。考虑温度差异不大,地层水矿化度一致的情况下交换阳离子的当量电导率B为定值,通过式(2)对比上部低阻气层和下部高阻气层粘土矿物的导电效率,结果表明低阻气层的粘土导电效率是高阻气层的3.56倍。因此,高粘土含量以及粘土附加导电性是形成Q3c上部低阻气层的重要成因。

Qv=CEC·ρG·(1-φt)φt
Rclay=1B·Qv

式中:Rclay为粘土矿物电阻率,Ω·m;B为交换阳离子的当量电导率,S·cm3/(mmol·m);Qv为岩石的阳离子交换容量,mmol/cm3;φt为岩石总孔隙度,式中取值为0.25;CEC为阳离子交换容量,mmol/g;ρG为骨架密度,g/cm3,式中取值为2.65。

2.2.3 高压异常及复杂孔隙结构

统计A-2井共计34个深度点的MDT测压值,Q3c下部3 633~3 634 m高阻气层3个有效测试点的地层压力均值为5 157 psi,而Q3c上部3 606.5~3 607.5 m低阻气层5个测试点的地层压力均值为5 297 psi,比下部高阻气层高140 psi,压力最高达到5 460 psi,压力系数达到1.064 5。根据前人的分类标准[11],压力系数在1.06~1.38之间的地层归类为高压异常地层,压力系数大于1.38为异常高压地层,地层压力趋势图表明Q3c上部低阻气层存在高压异常(图6)。

图6

图6   地层压力趋势

Fig.6   Formation pressure trend diagram


近年来,针对不同岩性地层的超压成因研究取得了重要进展,国内有学者将超压成因分为5种成因类型[12]:不均衡压实、流体膨胀、成岩作用、构造挤压、压力传递。砂岩等渗透性岩层的超压成因在之前大多研究中仅归结于压力传递,认为其异常压力来源于邻近泥质岩层或者远距离超压源的压力传递[13]。随着泥质岩层中生烃作用以及成岩作用等超压成因的重要性和普遍性逐渐被认知和证实,渗透性砂岩层中的超压成因应当被进一步认知。

鉴于H组泥岩总体为非烃源岩,且埋深大于3 400 m,不存在机械压实不均衡的因素。通过上述沉积相带以及X衍射资料的分析证实Q3c上部低阻气层相较于下部高阻气层具有更高的粘土矿物含量,因此本文重点研究成岩作用中蒙脱石—伊利石转化对Q3c上部低阻气层高压异常现象的影响。

据相关学者研究[14-16],成岩作用中蒙脱石—伊利石转化过程主要受随埋深和温度增加导致的矿物组分变化控制。粘土转化过程中的增压效应主要包括4个方面:转化过程使得粘土颗粒定向性得到增强,压实系数增大促进超压的形成;转化过程中蒙脱石层间吸附水和矿物转换产生的水大量脱出,致使地层产生附加的孔隙流体压力;转化过程中钾长石、云母等矿物遭到破坏和溶解,发生骨架弱化,上覆地层应力转移到孔隙流体系统产生附加的孔隙流体压力;伊利石矿物以及自生矿物的沉淀增强体系封闭性,保证了地层压力系统的稳定。普遍认为,当温度略高于60 ℃时,蒙脱石开始向伊利石转化。MDT测试Q3c地层温度在132.1~135.2 ℃,且Q3c上部地层中蒙脱石含量仅占粘土总量的7%,伊利石含量占比则达到56%,证实地层成岩作用处于蒙脱石向伊利石转化末期,粘土矿物转化过程中地层处于高压异常,在高压孔隙流体的保护下地层受到上覆地层的压实减轻,地层原生孔隙得到很好的保护。

通过图1铸体薄片对比,发现Q3c下部3 634 m高阻气层孔隙类型以小孔径的溶蚀粒间孔为主;上部3 607 m低阻气层各尺寸孔隙均有发育,孔隙类型以原生粒间孔和溶蚀粒间孔为主。A-2井核磁共振测井显示3 618~3 636 m高阻气层核磁T2谱主要分布在200 ms以内,谱峰相对较低,而3 592~3 608 m低阻气层T2谱谱峰则较饱满,T2谱分布范围更广(图7)。在高压异常地层条件下,铸体薄片、物性实验以及核磁共振测井均显示Q3c上部低阻气层孔隙更加发育,相较于Q3c下部高阻气层较差的物性,上部低阻气层较好的孔隙条件为优良导电路径的形成提供了很好的基础。

图7

图7   A-2井Q3c层段核磁共振测井响应

Fig.7   NMR logging response in Q3c section of Well A-2


通过对比Q3c上部低阻气层和下部高阻气层的岩心压汞实验数据以及核磁共振孔径分布,分析高、低阻地层孔隙结构差异(图8)。高、低阻4颗岩心压汞实验对比发现:两块低阻岩心的比表面积均值达到3 500 m2/g、体面积比均值达到7 871 m2/cm3,比两块高阻岩心的比表面积和体面积比均值分别大1 200 m2/g、43 861 m2/cm3;低阻气层段两颗岩心压汞曲线呈现明显的双平台特征,说明低阻岩心内部至少存在两种不同尺寸的孔隙分布。通过对核磁T2谱划定截止值区间分别得到低阻气层和高阻气层各尺寸孔喉分布,低阻气层各尺寸孔隙均有发育,孔隙类型以小孔隙为主,还存在一定量的中—大孔隙;高阻气层的各尺寸孔隙占比直方图近似呈正态分布,孔隙类型以小孔隙为主。通过核磁共振测井和压汞实验资料综合分析,相较于下部高阻气层,上部低阻气层更加复杂的孔隙结构,大量小孔隙与一定量的大孔隙组成复杂的双孔隙系统,这类地层具有更大的比表面积和体面积比。当气体充注入地层时,大尺寸孔隙内的地层水被优先驱替,随着气驱水过程的持续,充注动力不足导致微—小孔隙内部、细小喉道内以及碎屑颗粒表面仍保留大量束缚水。

图8

图8   高、低阻气层压汞曲线及核磁测井各尺寸孔隙占比

a—高、低阻各两块岩心压汞曲线对比;b—低阻气层中各尺寸孔隙占比;c—高阻气层中各尺寸孔隙占比

Fig.8   Mercury injection curve and the proportion of pores of various sizes in NMR logging

a—comparison of mercury intrusion curves of two cores with high and low resistance; b—the proportion of pores of each size in the low-barrier gas layer; c—the proportion of pores of each size in the high-barrier gas layer


Q3c上部低阻气层在高压异常条件的保护下仍具有较好的物性条件,为上部低阻气层提供了很好的导电路径。相较于下部高阻气层单一的孔隙结构,Q3c上部低阻气层由于复杂的孔隙结构具有更大的比表面积和体面积比,从而导致更高的束缚水饱和度。在较好的物性基础上,复杂孔隙结构导致更高的束缚水饱和度为Q3c上部低阻气层提供了更好的导电路径。因此,地层具有较好的物性基础上发育复杂孔隙结构是Q3c上部低阻气层的又一重要成因。

3 基于数字岩心微观可视化的定量研究

3.1 三维数字岩心模型建立

选取高、低阻气层两块岩心开展X射线CT扫描实验,扫描分辨率为1 μm,将CT扫描获取的灰度图像进行分割处理,构建了两块三维数字岩心模型。对高、低阻两块数字岩心进行组分划分得到岩石各组分分布形态以及相对含量,模型尺寸为500 μm×500 μm×500 μm,其中白色为石英、红色为长石、蓝色为伊利石、绿色为绿泥石、黑色为孔隙(图9)。通过数值模拟得到两块数字岩心的孔、渗参数与岩石物理实验结果具有较好的一致性(表5)。

图9

图9   高、低阻两块数字岩心的矿物组分分布以及孔隙网络模型

a—低阻岩心全貌;b—低阻岩心长石组分;c—低阻岩心伊利石组分;d—低阻岩心孔隙网络模型;e—高阻岩心全貌;f—高阻岩心长石组分;g—高阻岩心伊利石组分;h—高阻岩心孔隙网络模型

Fig.9   Mineral composition distribution and pore network model

a—overall picture of low resistivity core; b—feldspar composition in low resistivity cores; c—illite composition in low resistivity cores; d—pore network model for low resistivity cores; e—overall view of high resistivity core; f—feldspar composition in high resistivity core; g—illite components in high resistivity cores; h—pore network model for high resistivity cores


表5   高、低阻两块数字岩心各矿物组分含量及物性参数

Table 5  Contents of mineral components and porosity and permeability parameters

岩心类别石英/%长石/%伊利石/%绿泥石/%孔隙度/%渗透率/mD
数字岩心物性实验数字岩心物性实验
低阻岩心54.8720.138.692.5112.5813.201.681.82
高阻岩心57.6521.292.515.238.428.540.971.03

新窗口打开| 下载CSV


3.2 导电机理及电性模拟

通过高、低阻气层两块数字岩心模型,给定固定岩石矿物组分、孔隙内地层水和气体相应的电导率张量,采用有限元方法计算多组分三维数字岩心的电阻率[17-18]。根据实际地层测试资料,设定地层水电阻率为0.1 Ω·m进行电性模拟,各粘土矿物组分的电阻率取值参考表4中的阳离子交换容量CEC均值,结合式(1)、式(2)计算得到,其余岩石矿物和流体组分的电导率设为0 S/m。模拟得到高、低阻岩心在不同含水饱和度条件下的微观多组分分布和导电效率差异如图10,图中灰色为石英、红色为长石、绿色为绿泥石、黄色为伊利石、蓝色为地层水、粉色为气体。

图10

图10   高、低阻岩心电性模拟对比

a—饱含水高阻岩心各组分分布;b—高阻岩心饱含水情况下的电流分布;c—饱含水低阻岩心的组分分布;d—低阻岩心饱含水情况下的电流分布;e—含水饱和度为40%时高阻岩心各组分分布;f—高阻岩心40%含水饱和度条件下的电流分布;g—含水饱和度为40%时低阻岩心各组分分布;h—低阻岩心40%含水饱和度条件下的电流分布

Fig.10   Simulation and comparison of electrical properties of high and low resistivity cores

a—distribution of components in water-saturated high-resistance cores; b—current distribution in the case of high-resistance core saturated with water; c—component distribution of water-saturated low-resistance cores; d—current distribution in the case of water-saturated low-resistance cores; e—distribution of components of high resistivity core when water saturation is 40%; f—the current distribution of the high-resistance core under the condition of 40% water saturation; g—distribution of components in low-resistance cores when water saturation is 40%; h—the current distribution of the low-resistance core at 40% water saturation


饱含水条件下高阻岩心的粒间溶孔和长石粒内溶孔被地层水占据,孔隙内的地层水为主要的导电路径,模拟电阻率为8.87 Ω·m;40%含水饱和度条件下,气体占据了岩心内的中—小孔隙,微孔隙、角隅及吼道内仍保留束缚地层水,碎屑颗粒表面由于润湿性保留了一层水膜。高阻岩心电流密度分布显示,低含水饱和度情况下水膜、角隅以及微孔隙内的地层水构成其主要的导电路径,导电路径连通性较差,电流密度低,模拟电阻率高达46.05 Ω·m。

低阻岩心的伊利石含量明显高于高阻岩心,并且低阻岩心的长石组分相较于高阻岩心呈现更明显的分散状,孔隙网络模型显示低阻岩心内存在大量的微小孔隙,且发育一定量的中—大孔隙,表明低阻岩心孔隙结构复杂。在饱含水条件下,低阻岩心孔隙内的地层水和伊利石粘土矿物构成了良好的导电路径,模拟电阻率仅为2.67 Ω·m;40%含水饱和度条件下,气体聚集在大尺寸孔隙内,连通性不好的中—小孔隙内仍保留大量的束缚水,加之高阳离子交换容量的伊利石矿物大量存在,在低含水饱和度时仍然构成了很好的导电路径。模拟40%含水饱和度时低阻岩心的电阻率仅为10.46 Ω·m。

高、低阻两块数字岩心电性模拟结果见图11。两块岩心电阻率随含水饱和度变化差异明显,高阻岩心含水饱和度由100%减小到40%时,电阻率由8.87 Ω·m增大到46.05 Ω·m;低阻岩心含水饱和度由100%减小到40%时,电阻率由2.67 Ω·m增大到10.46 Ω·m。高、低阻两块岩心的饱和度指数分别为1.818、1.505,低阻岩心的饱和度指数明显低于高阻岩心。在三维数字岩心可视化基础上进行多组分导电模拟研究高、低阻气层导电机理和导电效率,研究发现基于数字岩心的多组分导电模拟结果与测井电性特征一致,证实研究区Q3c上部低阻气层由粘土附加导电性和复杂孔隙结构共同作用导致。

图11

图11   高、低阻两块数字岩心电性模拟结果

a—高、低阻两块岩心不同含水饱和度下电阻率的变化;b—高、低阻两块岩心电阻增大率RI与含水饱和度的拟合

Fig.11   Electrical simulation results of two digital cores with high and low resistance

a—resistivity changes of two cores with high and low resistivity under different water saturation;b—fitting diagram of resistance increase rate RI and water saturation of two cores with high and low resistivity


3.3 基于数字岩心的低阻各因素贡献定量分析

以上述2类数字岩心为基础,采用过程法分别建立相同物性条件下,阳离子交换容量不同的2块数字岩心,采用有限元算法计算4块数字岩心在不同含水饱和度条件下的电阻率。以实际地层条件下高阻气层岩心的电阻率为标准,研究不同孔隙度、阳离子交换容量和束缚水饱和度导致的电阻率变化。

通过对比相同物性条件下,阳离子交换容量变化导致的电阻率差异来定量分析粘土附加导电性对低阻的贡献;通过对比相同阳离子交换容量条件下,不同物性基础上由束缚水饱和度差异导致的电阻率变化来定量研究物性差异下复杂孔隙结构对低阻的贡献。根据测试及水分析资料,模拟得到4块数字岩心在130 ℃地层温度、0.1 Ω·m地层水电阻率条件下电阻率随含水饱和度的变化如图12。图中4块岩心电阻率随含水饱和度的变化均呈幂函数关系。

图12

图12   4块岩心电阻率随含水饱和度的变化

Fig.12   Variation of resistivity of 4 cores with water saturation


以核磁共振测井计算得到的束缚水饱和度为参考,认为高、低阻岩心束缚水饱和度分别为0.4、0.6时的电阻率为实际地层条件下的电阻率,模拟得到实际高、低阻气层电阻率分别为39.41 Ω·m、5.82 Ω·m,电阻率差异为33.59 Ω·m。

模拟相同物性条件下由于阳离子交换容量差异导致的电阻率变化研究粘土附加导电性对低阻的贡献。由图12可知,孔隙度为8.42%岩心的阳离子交换容量由0.082 mmol/cm3增加到0.292 mmol/cm3时,模拟得到不同含水饱和度情况的电阻率均大幅降低。当含水饱和度为0.4时,电阻率由39.41 Ω·m降低为22.59 Ω·m,电阻率下降16.82 Ω·m。

图12可知,当含水饱和度为0.4时,孔隙度为8.42%的岩心阳离子交换容量由0.082 mmol/cm3增加到0.292 mmol/cm3时,电阻率由39.41 Ω·m降低为22.59 Ω·m,电阻率下降16.82 Ω·m。

考虑相同阳离子交换容量条件下物性差异和孔隙结构差异对电阻率降低的贡献。由上文可知,孔隙结构差异导致束缚水饱和度存在明显差异,通过模拟相同物性条件下,不同束缚水饱和度条件下的电阻率变化研究孔隙结构差异对低阻的影响。模拟结果显示当孔隙度为8.42%、阳离子交换容量为0.082 mmol/cm3的岩心束缚水饱和度由0.4提高到0.6时,电阻率由46.05 Ω·m降低到23.52 Ω·m,电阻率下降22.53 Ω·m;当孔隙度为12.58%,阳离子交换容量为0.082 mmol/cm3的岩心束缚水饱和度由0.4提高到0.6时,电阻率由18.24 Ω·m降低到9.02 Ω·m,电阻率下降9.22 Ω·m。相较于Q3c下部高阻气层,上部低阻气层发育更好的物性和复杂的孔隙结构,模拟孔隙度由8.41%增加到12.58%、束缚水饱和度由0.4提高到0.6时岩心电阻率的变化,研究在实际地层条件下物性和孔隙结构对电阻率的影响。模拟结果显示电阻率由39.41 Ω·m下降为9.02 Ω·m,电阻率下降30.40 Ω·m。对比前两组实验数据发现:对于物性越差的地层,孔隙结构差异对其电阻率的影响更大;对比第三组实验和前两组实验数据发现:在更好的物性条件前提下,复杂的孔隙结构对电阻率的影响远大于单一复杂孔隙结构对电阻率的影响。因此,在考虑复杂孔隙结构对低阻油气藏的影响时,需重点考虑地层物性条件是否存在差异这一前提,而非将低阻成因直接归因于单一复杂孔隙结构。

通过有限元方法计算得到高、低阻气层电阻率差异为33.59 Ω·m,粘土附加导电性导致电阻率降低16.82 Ω·m,物性差异条件下复杂孔隙结构导致电阻率下降30.40 Ω·m,上述实验各参数取值及各低阻因素导致的电阻率下降量见表6。将两个低阻成因导致电阻率下降的值进行归一化处理,得到各个成因对研究区Q3c上部气层低阻响应的贡献。研究结果表明粘土附加导电性对低阻响应的贡献为35.63%,良好物性条件下复杂孔隙结构对低阻响应的贡献高达64.37%。

表6   各低阻成因电阻率下降量分析及参数选取

Table 6  Analysis of the decrease in resistivity of each low resistance cause and parameter selection

粘土附加导电性定量条件φ=0.0842、Sw=0.4
ΔR1(Qv=0.082→0.292)/(Ω·m)16.82
物性/孔隙结构差异定量条件Qv=0.082
孔隙度变量φ=0.0842φ=0.1258φ=0.08→0.12
ΔR2(Sw=0.4→0.6)/Ω·m22.539.2230.4

新窗口打开| 下载CSV


4 结论及认识

1)通过测井资料、气测录井和生产测试资料证实A构造3口井同一地质层位Q3c上部均出现低阻气层。以丰富的测井资料和大量岩石物理实验资料为研究依据,分别研究了外因和内因对地层电阻率的影响,进而确定低阻成因。研究发现低阻气层粘土含量高,粘土矿物类型以高阳离子交换容量的伊蒙混层为主,定量分析发现低阻气层粘土矿物导电效率是高阻气层的3.7倍。高粘土含量以及粘土附加导电性是形成Q3c上部低阻气层的重要成因。

2)研究区Q3c上部低阻气层中的粘土矿物处于成岩作用转化末期,蒙脱石向伊利石转化过程中地层保持高压异常,有利于地层原生粒间孔隙的保存,为上部低阻气层优良导电路径的形成提供了很好的物性基础。在此基础上,低阻气层内部发育的大量中—小孔隙与一定量的原生粒间大孔隙组成的双孔隙系统大大提高了地层孔隙结构的复杂性,低阻地层具有更高的束缚水饱和度,这是形成Q3c上部低阻气层的另一重要成因。

3)建立了高、低阻气层两块数字岩心,在三维数字岩心可视化基础上进行多组分导电模拟。对低阻气层的导电机理和导电效率开展研究发现:低阻岩心在高含气饱和度时,水膜、角隅和微孔隙内残余的束缚水构成其主要导电路径,导电效率低;低阻岩心内部存在大量的伊利石粘土矿物,且由于孔隙结构复杂,连通性较差的中—小孔隙内仍保留大量束缚水,在高含气饱和度情况下,中小孔内的束缚水和高阳离子交换容量的伊利石矿物构建了优良的导电路径,模拟结果与测井电性特征吻合,证实了粘土附加导电性和复杂孔隙结构是导致研究区Q3c上部气层低阻的主要原因。

4)模拟物性、阳离子交换容量存在差异的4块数字岩心在不同含水饱和度条件下的电阻率变化,定量分析粘土附加导电性和物性差异下复杂孔隙结构对低阻的贡献。模拟发现,对于研究区这类低孔、低渗—特低渗地层,物性条件越差,单一孔隙结构差异对其电阻率的影响越明显;而对于物性基础更好的地层,复杂孔隙结构对其电阻率的影响远大于单一复杂孔隙结构对其电阻率的影响,因此,在考虑复杂孔隙结构对低阻油气藏的影响时,需重点考虑地层物性条件是否存在差异这一前提,而非将低阻成因直接归因于单一复杂孔隙结构。模拟结果显示粘土附加导电性对研究区Q3c上部气层低阻响应的贡献为35.63%,良好物性条件下复杂孔隙结构对低阻响应的贡献高达64.37%。

参考文献

陶士振, 邹才能.

东海盆地西湖凹陷天然气成藏及分布规律

[J]. 石油勘探与开发, 2005, 32(4): 103-110.

[本文引用: 1]

Tao S Z, Zou C N.

Accumulation and disrtibution of natural gases in Xihu Sag, East China Sea Basin

[J]. Petroleum Exploration & Development, 2005, 32(4): 103-110.

[本文引用: 1]

何将启, 梁世友, 陈拥锋, .

东海盆地西湖凹陷新生代构造演化对油气的控制作用——以平湖组油气响应为例

[J]. 石油实验地质, 2008, 30(3): 221-226.

[本文引用: 1]

He J Q, Liang S Y, Chen Y F, et al.

Control on petroleum by Cenozoic tectonic evolution in the Xihu Sag, the East China Sea Basin: Taking petroleum response of the Pinghu Foramation as an example

[J]. Petroleum Geology & Experiment, 2008, 30(3):221-226.

[本文引用: 1]

张银国.

东海西湖凹陷花港组油气地质条件与油气分布规律

[J]. 石油实验地质, 2010, 32(3): 223-241.

[本文引用: 1]

Zhang Y G.

Petroleum Geology and hydrocarbom distribution pattern of Huagang Formation in the Xihu sag of the East China Sea

[J]. Petroleum Geology & Experiment, 2010, 32(3): 223-241.

[本文引用: 1]

董春梅, 赵仲祥, 张宪国, .

西湖凹陷中北部花港组物源及沉积相分析

[J]. 东北石油大学学报, 2018, 42(5): 25-34.

[本文引用: 1]

Dong C M, Zhao Z X, Zhang X G, et al.

Analysis of provenance and sedimentary facies of Huagang formation in the north central of Xihu Sag

[J]. Journal of Noryheast Petroleum University, 2018, 42(5): 25-34.

[本文引用: 1]

胡明毅, 柯岭, 梁建设, .

西湖凹陷花港组沉积相特征及相模式

[J]. 石油天然气学报, 2010, 32(5): 1-5.

[本文引用: 1]

Hu M Y, Ke L, Liang J S, et al.

The characteristics and pattern of sedimentary facies of Huagang Formation in Xihu Depression

[J]. Journal of Oil and Gas Technology, 2010, 32(5): 1-5.

[本文引用: 1]

姜艳娇. A地区低孔渗复杂储层导电机理及产能预测方法研究[D]. 东营: 中国石油大学(华东), 2017.

[本文引用: 1]

Jiang Y J. The research of conductive mechanism and productivity prediction of low porosity and low permeability reservoir in A area[D]. Dongying: China University of Petroleum(East China), 2017.

[本文引用: 1]

王迪, 戚家振, 陈现, .

东海N气田低阻气层成因分析及饱和度定量评价

[J]. 复杂油气藏, 2017, 10(4):7-13.

[本文引用: 1]

Wang D, Qi J Z, Chen X, et al.

Forming reason analysis and saturation quantitative evaluation of low-resistivity gas layer in N Gas field of Donghai Sea

[J]. Complex Hydrocarbon Reservoirs, 2017, 10(4): 7-13.

[本文引用: 1]

徐昉昊, 徐国盛, 刘勇, .

东海西湖凹陷中央反转构造带古近系花港组致密砂岩储集层控制因素

[J]. 石油勘探与开发, 2020, 47(1): 98-109.

[本文引用: 1]

Xu F H, Xu G S, Liu Y, et al.

Factors controlling the development of tight sandstone reservoirs in the Huagang Formation of the central inverted structual belt in Xihu sag, East China Sea Basin

[J]. Petroleum Exploration and Development, 2020, 47(1): 98-109.

[本文引用: 1]

刘金水, 曹冰, 徐志星, .

西湖凹陷某构造和花港组沉积相及致密砂岩储层特征

[J]. 成都理工大学学报:自然科学版, 2012, 39(2): 130-136.

[本文引用: 1]

Liu J S, Cao B, Xu Z X, et al.

Sedimentary facies and the characteristics of tight sandstone reservoirs of Huagang Formation in Xihu Depression, East China Sea Basin

[J]. Journal of Chengdu Unibersity of Technology:Natural Science Edition, 2012, 39(2): 130-136.

[本文引用: 1]

程相志. 低阻油气层识别评价技术及分布规律研究[D]. 东营: 中国石油大学(华东), 2008.

[本文引用: 1]

Cheng X Z. Study of recognition technology and distribution law on low-resistivity oil reservoir[D]. Dongying: China University of Petroleum(East China), 2008.

[本文引用: 1]

杜栩, 郑洪印, 焦秀琼.

异常压力与油气分布

[J]. 地学前缘, 1995, 2(3/4):137-148.

[本文引用: 1]

Du X, Zheng H Y, Jiao X Q.

Abnormal pressure and hydrocarbon accumulation

[J]. Earth Science Frontiers, 1995, 2(3/4): 137-148.

[本文引用: 1]

赵靖舟, 李军, 徐泽阳.

沉积盆地超压成因研究进展

[J]. 石油学报. 2017, 38(9): 973-998.

DOI:10.7623/syxb201709001      [本文引用: 1]

近十多年来,超压成因研究在国际上取得了重要进展,主要表现在5个方面。1将超压按照成因分为5种类型:不均衡压实、流体膨胀、成岩作用、构造挤压、压力传递;其中流体膨胀超压的成因包括生烃作用、油裂解气、水热膨胀等,成岩作用则包括蒙脱石—伊利石转化作用等。2总结提出了6种超压成因判识方法,包括测井曲线组合分析法、鲍尔斯法(加载-卸载曲线法)、声波速度-密度交会图法、孔隙度对比法、压力计算反推法以及综合分析法。3随着实证超压成因研究方法的逐渐广泛应用,许多被普遍认为属于不均衡压实成因的超压已被完全或部分否定,生烃作用作为超压成因的重要性和普遍性正在得到愈来愈多证实,黏土矿物成岩作用特别是蒙脱石—伊利石转化作用对超压形成的重要性也已受到重视,构造挤压和压力传递成因的超压得到进一步确认,而且现已发现许多盆地或地区的超压可能为复合成因。4不同的岩性其超压成因往往有别:就泥质岩而言,烃源岩与非烃源岩的超压成因常常不同,通常烃源岩内发育的超压大多与生烃作用密切相关,有时成岩作用也具有重要贡献;非烃源岩内所发育的超压多与不均衡压实、成岩作用或压力传递等有关;砂岩等渗透性岩层的超压则多为压力传递成因。5对于烃源岩层段的超压成因分析而言,由于有机质含量对密度、声波等测井参数具有明显影响,因此在运用这些测井资料分析烃源岩层段的超压成因时,需要进行有机质含量校正,校正前后超压成因的分析结果往往不同。

Zhao J Z, Li J, Xu Z Y.

Advances in the origin of overpressure in sedimentary basins

[J]. Acta Petrolei Sinica, 2017, 38(9): 973-998.

[本文引用: 1]

Osborne M J, Swarbrick R E.

Mechanisms for generating overpressure in sedimentary basins: A reevaluation: Reply

[J]. AAPG Bulletin, 2001, 85(12): 2119-2119.

[本文引用: 1]

李超, 罗晓容, 张立宽.

泥岩化学压实作用的超压响应与孔隙压力预测

[J]. 中国矿业大学学报, 2020, 49(5): 951-973.

[本文引用: 1]

Li C, Luo X R, Zhang L K.

Overpressure responses for chemical compaction of mudstones and the pore pressure prediction. Jouranl of China University of Mining and Technology

[J]. Journal of China University of Ming and Technology, 2020, 49(5): 951-973.

[本文引用: 1]

褚庆忠.

异常压力形成机制研究综述

[J]. 天然气勘探与开发, 2001, 24(4): 38-46.

[本文引用: 1]

Chu Q Z.

Review of abnormal pressure formation mechanism

[J]. Natural Gas Exploration and Development, 2001, 24(4): 38-46.

[本文引用: 1]

查明, 曲江秀, 张卫海.

异常高压与油气成藏机理

[J]. 石油勘探与开发, 2002, 29(1): 19-23.

[本文引用: 1]

Zha M, Qu J X, Zhang W H.

The relationship between overpressure and reservoir forming mechanism

[J]. Petroleum Exploration and Development, 2002, 29(1): 19-23.

[本文引用: 1]

Yan W C, Sun J M, Zhang J Y, et al.

Studies of electrical properties of low-resistivity sandstones based on digital rock technology

[J]. Journal of Geophysics and Engineering, 2018, 15(1):153-163.

DOI:10.1088/1742-2140/aa8715      URL     [本文引用: 1]

李霞, 李潮流, 李波.

致密砂岩岩电响应规律与饱和度评价方法

[J]. 石油勘探与开发, 2020, 47(1): 202-212.

[本文引用: 1]

Li X, Li C L, Li B.

Response laws of rock electrical property and saturation evaluation method of tight sandstone

[J]. Petroleum Exploration and Development, 2020, 47(1): 202-212.

[本文引用: 1]

/

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com