氮气吸附滞后回环定量分析及其在孔隙结构表征中的指示意义
——以鄂尔多斯盆地上三叠统延长组7段为例
1.
2.
3.
Quantitative analysis of nitrogen adsorption hysteresis loop and its indicative significance to pore structure characterization:A case study on the Upper Triassic Chang 7 Member, Ordos Basin
1.
2.
3.
通讯作者: 姜福杰(1979—),男,博士、教授,油气成藏机理与分布规律。E‑mail:jiangfj@cup.edu.cn。
编辑: 张玉银
收稿日期: 2022-06-06 修回日期: 2022-12-15
基金项目: |
|
Received: 2022-06-06 Revised: 2022-12-15
第一作者简介 About authors
宋泽章(1988—),男,博士、副教授,非常规储层评价、油气运聚与成藏、机器学习。E‑mail:
低温氮气吸附滞后回环的形态和面积可以有效反映多孔介质的孔隙结构及其对吸附气体的滞留效应,但其在页岩孔隙结构定量表征中的作用常被忽视。优选鄂尔多斯盆地上三叠统延长组7段(长7段)页岩为研究对象,综合场发射扫描电镜孔隙结构定性观察、低温氮气吸附孔隙结构定量分析、“滞后回环”定量分析、总有机碳含量(TOC)分析、热解实验、X射线衍射实验等手段,对页岩在低温氮气吸附-脱附实验中能否形成“滞后回环”以及形成的“滞后回环”面积大小的决定因素进行了探索。研究结果表明:①页岩样品在低温氮气吸附-脱附实验中是否能形成滞后回环,与其比表面积、比孔容、黏土矿物含量及孔隙结构分形维数之间存在明显的正相关关系,与总有机碳含量呈明显负相关关系,而与平均孔径大小、孔隙表面分形维数、最高热解峰温和脆性矿物含量等无明显相关关系;②滞后回环的面积大小取决于两端开放的圆柱形孔、墨水瓶孔或平行板孔的发育程度,可借助滞后回环的面积大小来定量评价圆柱形孔、墨水瓶孔或平行板孔占孔隙空间的相对比例;③长7段泥页岩样品中开放的圆柱形孔、墨水瓶孔或平行板孔主要由黏土矿物晶间孔提供,滞后回环面积与泥页岩样品中黏土矿物含量之间存在明显的正相关关系。
关键词:
Low-temperature nitrogen adsorption can form a “hysteresis loop”, whose geometry and area can effectively reflect the pore structure of porous media and its retention effect on adsorbed gas. However, the role of the “hysteresis loop” in quantitative characterization of shale pore structure has often been ignored. The study aims to clarify whether shale can form “hysteresis loop” in low-temperature nitrogen adsorption-desorption experiment and the determinants on hysteresis loop’s area with experiments on the 7th member of the Upper Triassic Yanchang Formation shale (Chang 7 shale) in the Ordos Basin. Various measures are applied in the study, including qualitative observation of pore structure under field-emission scanning electron microscopy (FE-SEM) and quantitative characterization of pore structure by low-temperature nitrogen adsorption test, hysteresis-loop quantitative analysis, total organic carbon (TOC) content analysis, pyrolysis and X-ray diffraction (XRD) experiments. The following results are obtained. First, whether shale can form a hysteresis loop in the low-temperature nitrogen adsorption-desorption experiment has an apparently positive correlation with the specific surface area, specific pore volume, clay mineral content, and pore structure fractal dimension, and an evidently negative correlation with the TOC content, while no apparent correlation with the average pore size, pore surface fractal dimension, highest pyrolysis peak temperature, and content of brittle minerals. Second, the hysteresis loop area depends on the development degree of the cylindrical pores with both ends open, ink-bottle pores, or parallel plate pores, the proportion of which to the total pores can be quantitatively evaluated by the hysteresis loop area. Third, the open-ended cylindrical pores, ink-bottle pores, or parallel plate pores in the samples from Chang 7 shale are mainly of the intergranular pores in clay minerals. Therefore, there is an apparently positive correlation between the hysteresis loop area and the content of clay minerals.
Keywords:
本文引用格式
宋泽章, 阿比德·阿不拉, 吕明阳等.
SONG Zezhang, ABIDE Abula, LYU Mingyang, et al.
低温氮气吸附-脱附等温曲线直观反映了单位质量样品的吸附量与相对压力之间的关系,是氮气吸附技术得到的最基础数据。前人研究表明,氮气吸附-脱附曲线的形态受微观孔隙结构控制[6-9]。早在1940年,Brunauer等[10]就在微观材料学中提出了5类典型吸附-脱附等温曲线,并阐述了不同等温曲线与材料微观孔隙结构之间的关系。目前,应用最为广泛的是国际纯化学与应用化学联合会(IUPAC)对氮气吸附-脱附等温曲线形态的分类方法(6类):Ⅰ型一般指示多孔介质内部主要为微孔;Ⅱ型反映孔隙空间为大孔或不存在孔隙;Ⅲ型同样指示孔隙空间为大孔或不存在孔隙,但气-固相互作用相对Ⅱ型更弱;Ⅳ型说明中孔主导了孔隙空间;Ⅴ型多由微孔和中孔复合的孔隙空间(弱气-固相互作用)产生;Ⅵ型指示氮气在均匀、非孔表面的顺次多层吸附[9, 11]。
在吸附过程中,经常会发生毛细凝聚现象,即在毛细通道内的气体,不必达到过饱和状态即可发生凝聚,产生凝聚液。在开放的孔隙系统中,由于发生毛细凝聚和蒸发时,气、液两相的界面形态不同,导致气体在吸附和脱附过程中达到同一吸附量的相对压力不同,从而形成滞后回环[12]。滞后回环被广泛应用于多孔介质孔隙结构定性表征[5, 13-15]。早在1958年,De Boer[16]便归纳了5类典型的滞后回环,对应不同形态的孔隙:A型指示管状毛细孔;B型指示狭缝状毛细孔;C型指示锥形或双锥形管状毛细孔;D型指示楔形毛细孔;E型指示墨水瓶状孔隙[17]。IUPAC[8, 18]则将滞后回环分为4类:尺寸和排列都十分规则的孔隙结构常形成H1型回环;若微孔主导孔隙空间则会产生H4型回环;而当多孔介质为无规则孔隙结构时,则主要产生H2和H3型回环。吴俊等[19]则归纳了3类多孔介质孔隙形态对滞后回环产生与否的影响:①开放式孔隙,包括四边开口的平行板孔隙和两端开口的圆柱形孔隙,能够产生滞后回环;②一端封闭的孔隙,包括一端封闭的楔形孔、锥形孔、平行板孔和圆柱形孔,不能产生滞后回环;③“墨水瓶”孔,虽然一端是封闭的,但仍然可以产生一个滞后回环,并且脱附曲线通常有一个急剧下降的拐点[13]。
截至目前,滞后回环定量表征在页岩孔隙结构定量表征中的应用尚少见报道,亟需深化这方面的工作。为此,本文优选中国典型陆相页岩——鄂尔多斯盆地陇东地区上三叠统延长组7段(长7段)页岩为研究对象,使用低温氮气吸附技术表征长7段泥页岩的比表面积、比孔容、孔径分布、孔隙结构类型并定量表征滞后回环面积;而后,使用场发射扫描电镜技术对研究区(陇东地区)纳米级孔隙空间进行了系统的精细刻画;最后,综合总有机碳含量(TOC)分析、热解实验、XRD全岩分析等技术探讨有机质丰度、成熟度、矿物组成、非均质性、孔隙结构参数等对陆相页岩在氮气吸附实验中能否产生滞后回环的指示意义以及对滞后回环面积大小的控制作用。
1 地质背景
鄂尔多斯盆地位于华北地块西部,页岩油气资源丰富,是一个多期构造叠合形成的克拉通盆地[5, 23]。从中生代开始,华北克拉通逐步进入差异升降阶段,受太平洋板块向华北板块俯冲作用的影响,盆地东部在走滑-挤压作用下逐步抬升,形成了西倾的斜坡构造格局[12, 24-25]。构造上,鄂尔多斯盆地可以划分为6个一级构造单元,而研究区位于陕北斜坡西南部(图1),页岩油勘探已经取得了重大进展,发现了10亿吨级页岩油田(庆城大油田),建成了中国第一个百万吨级整装页岩油示范区[26-27]。中生界上三叠统长7段陆相暗色泥页岩是中国最典型的页岩油层系之一,具备埋藏浅(2 000~3 000 m)、厚度大(30~160 m)、分布稳定、有机质丰度高(TOC =1 %~15 %)、成熟度适中(镜质体反射率Ro=0.5 %~1.2 %)等特征[28-30]。本文所选的9块长7段泥页岩样品沉积环境均为半深湖-深湖相,埋藏深度介于1 563~3 058 m,具有较好的代表性。
图1
图1
鄂尔多斯盆地陇东地区区域构造背景及地层
a. 区域地质构造;b. 地层综合柱状图
Fig. 1
Regional tectonics and stratigraphic column of Longdong area, Ordos Basin
2 样品制备及实验过程
全岩矿物X射线衍射分析是在核工业北京地质研究院完成的。依照中国石油天然气行业标准SY/T5163—2010《沉积岩中黏土矿物和常见非黏土矿物X射线衍射分析方法》,采用Panalytical X’ Pert PRO X射线衍射仪进行测试。TOC和岩石热解分析则是在中国石油大学(北京)油气资源与探测国家重点实验室分别采用CS-230HC碳硫分析仪和Rock-Eval热解仪完成的。孔隙结构定性观察则是在中国科学院地质与地球物理研究所经氩离子抛光、镀碳后,借助Nova NanoSEM 450场发射扫描电子显微镜(FE-SEM)观察完成(最高分辨率可达1.0 nm)。
在进行氮气吸附实验前,必须对页岩样品进行前处理,尽量避免原油堵塞孔隙空间、导致氮气吸附实验表征精度降低:将泥页岩样品粉碎至60~80目,在油气资源与探测国家重点实验室使用自动多功能抽提仪进行索氏抽提。为了尽量抽提出孔隙空间内的滞留烃类,所有样品均经历了172 h索氏抽提。
氮气吸附实验是根据SY/T 6154—2019《岩石比表面积和孔径分布测定静态吸附容量法》标准,在北京市理化分析测试中心,采用Quadrasorb Station型比表面分析仪在-196 ℃下完成的。实验前,样品需在110 ℃的真空条件下经过8 h脱气处理。而后,在低温液氮条件下,将吸附质(N2)流经测试样品,设置相对压力平衡时间和差值,待吸附平衡后继续升高压力进行吸附实验;当相对压力接近液氮温度下的饱和蒸汽压(p0≈0.101 MPa)时,逐步降低压力,直至氮气从样品中脱附完成。根据氮气吸附的微孔充填和毛细凝聚理论[31],在较低的相对压力下,由于微孔吸附能力增强,氮气会优先被其捕捉;随着压力的增大,毛细凝聚从小孔到大孔逐渐推进,而脱附时则按照从大孔到小孔顺序进行[31-32]。因此,通过记录整个过程中氮气吸附量的变化可计算出不同孔径下的孔隙体积,饱和状态下的氮气总吸附量为样品内孔隙的总体积。
3 实验结果
3.1 XRD全岩定量分析
XRD全岩定量分析结果显示,长7段泥页岩矿物成分中,石英和黏土矿物明显占主导地位:石英含量在10.8 %~45.3 %,平均为29.9 %;黏土矿物含量在30.5 %~71.0 %,平均为53.4 %(表1[33];图2a)。此外,还含有少量长石(0~21.7 %,平均9.6 %)、黄铁矿(0~17.3 %,平均6.3 %)和微量的菱铁矿(0~6.5 %,平均0.7 %)。脆性矿物(石英+长石)含量较低,分布在24.2 %~52.2 %,平均值为39.5 %;而黏土矿物平均值高于脆性矿物,因此,可以推测长7段泥页岩的脆性和可压性不高。参考王玉满等[34]基于石英+长石、黏土矿物和碳酸盐岩三端元页岩岩相分类图版(图2b),长7段泥页岩主要以黏土质页岩为主,硅质页岩次之。
表1 鄂尔多斯盆地陇东地区长7段泥页岩样品XRD分析测试结果
Table 1
样品号 | 全岩矿物组分含量/% | |||||||
---|---|---|---|---|---|---|---|---|
石英 | 钾长石 | 斜长石 | 方解石 | 白云石 | 菱铁矿 | 黄铁矿 | 黏土矿物 | |
Z22-1584.65* | 22.2 | 0 | 6.8 | 0 | 0 | 0 | 0 | 71.0 |
Z22-1626.00 | 33.4 | 0 | 10.9 | 0 | 0 | 0 | 0 | 55.7 |
Z22-1632.45 | 10.8 | 4.7 | 8.7 | 0 | 0 | 6.5 | 0 | 69.3 |
LI231-2105.90 | 32.4 | 0 | 0 | 0 | 0 | 0 | 12.7 | 54.9 |
YAN56-3032.10* | 42.6 | 0 | 5.6 | 0 | 0 | 0 | 0 | 51.8 |
YAN56-3058.40 | 45.3 | 0 | 0.0 | 0 | 0 | 0 | 1.9 | 52.8 |
Z22-1563.40 | 21.2 | 0 | 18.3 | 0 | 0 | 0 | 13.5 | 47.0 |
Z22-1569.10 | 19.1 | 0 | 21.7 | 0 | 0 | 0 | 11.4 | 47.8 |
LI68-2078.90* | 42.3 | 0 | 9.9 | 0 | 0 | 0 | 17.3 | 30.5 |
图2
图2
鄂尔多斯盆地陇东地区长7段泥页岩矿物组成及岩相划分
a.矿物组成分布条形图;b.页岩岩相划分三端元图(据文献[34]修改)
Fig. 2
Mineral compositions and lithofacies division of the Chang 7 shale in Longdong area, Ordos Basin
3.2 热解实验
热解实验结果(图3)显示:长7段泥页岩的氢指数(HI)分布在94.34~365.07 mg/g,平均为211.96 mg/g;最高热解峰温(Tmax)分布在448~467 ℃,平均为457 ℃;以Ⅰ型(腐泥型)和Ⅱ1型干酪根为主。此外,长7段泥页岩样品均处于成熟阶段(Easy%Ro=0.7 %~1.3 %)。
图3
图3
鄂尔多斯盆地陇东地区长7段泥页岩干酪根类型及成熟度判别
Fig. 3
Kerogen type and maturity identification for the Chang 7 shale in Longdong area, Ordos Basin
长7段泥页岩样品的残留烃含量(S1)分布在1.00~3.35 mg/g,热解烃含量(S2)分布在3.30~49.58 mg/g,样品生烃潜量(S1+S2)平均可达22.77 mg/g,反映样品生烃潜量较好-极好。有机质丰度同样显示,长7段泥页岩样品TOC分布在3.50 %~16.36 %,平均高达9.00 %(所有样品TOC>2.00 %)属于好-极好的烃源岩(图4)。
图4
图4
鄂尔多斯盆地陇东地区长7段泥页岩烃源岩质量分析
Fig. 4
Quality analysis of source rocks in the Chang 7 shale of Longdong area, Ordos Basin
3.3 低温氮气吸附实验
3.3.1 低温氮气吸附-脱附曲线
整体上,长7段泥页岩样品氮气吸附-脱附曲线形态均具有反“S”形的三段式特征(图5):①低压段(相对压力p/p0<0.05),吸附量上升较快,吸附曲线向上微凸,表明微孔对氮气有较强的吸附作用;②过渡段(p/p0=0.05~0.80),吸附量增加缓慢,部分样品与低压段存在明显的拐点;③高压段(p/p0>0.80),吸附量急剧上升,且当相对压力接近饱和蒸汽压时没有出现平缓段,未达到吸附饱和状态,表明样品中存在一定量的中孔和大孔。
图5
图5
鄂尔多斯盆地陇东地区长7段泥页岩样品3类典型低温氮气吸附-脱附等温曲线
Fig. 5
Low-temperature nitrogen adsorption isotherms of three types for samples from the Chang 7 shale of Longdong area, Ordos Basin[回环面积单位:mL/(nm·g)。]
研究区9个典型页岩样品氮气吸附-脱附曲线表现出不同程度的“滞后现象”。依据IUPAC的等温曲线和“滞后回环”分类方案,将研究区长7段泥页岩样品分为3类。
3) C类(图5c):与A类和B类样品形成鲜明对比的是,C类样品氮气吸附-脱附分支几乎完全重合,没有出现明显的“滞后回环”,等温曲线类型可以归为Ⅱ型。
Barrett-Joyner-Halenda(BJH)和密度泛函理论(DFT)是用于定量分析氮气吸附实验结果、获取多孔介质孔径分布最常用的两类模型[35-36]。然而,基于热力学Kelvin方程的BJH模型并不适于微孔孔径分析。首先,微孔中相邻孔隙壁之间的作用力很强,基于热力学的Kelvin方程并不适用[37];其次,使用BJH模型分析相对比较小的介孔(孔径2~10 nm)时常常会低估孔径(对孔径<4 nm的介孔会造成高达20 %的偏差[38]);此外,当使用BJH模型分析有“陡降”脱附分支(本研究中A类样品)的数据时,其所得孔径分布图将出现假峰[38]。与之相比,DFT模型则因其比基于开尔文方程的模型更能准确地估计孔径分布而广受欢迎。DFT模型基于分子模拟技术,从分子尺度来描述孔隙中吸附质的状态[39],兼顾了气-固、气-液相互作用以及孔隙几何因素的影响,能够更真实地反映孔隙中束缚流体的热力学性质[33, 40],因而是一种更加有效的分析手段,能够提供更加精确的孔径分布信息。基于此,本研究中采用DFT模型对页岩样品的比表面积、比孔容、平均孔径等孔隙结构参数及孔径分布进行解释。而后,对吸附-脱附曲线进行插值,加密数据点,再以吸附曲线为基线,对脱附曲线进行定积分计算,求得滞后回环的绝对面积[mL/(nm·g)](表2)。
表2 鄂尔多斯盆地陇东地区长7段泥页岩样品氮气吸附实验孔隙结构参数统计
Table 2
滞后回环类型 | 样品编号 | 深度/m | 滞后回环面积/ (mL·nm-1·g-1) | DFT比表面积/ (m2·g-1) | DFT比孔容/ (mL·g-1) | DFT平均孔径/nm | 孔隙表面分形维数(D1) | 孔隙结构分形维数(D2) |
---|---|---|---|---|---|---|---|---|
A类 | Z22 | 1584.65 | 0.578 68 | 13.840 | 0.048 55 | 7.032 | 2.518 | 2.609 |
Z22 | 1626.00 | 0.565 96 | 15.080 | 0.048 05 | 7.032 | 2.533 | 2.605 | |
Z22 | 1632.45 | 0.711 39 | 20.210 | 0.060 20 | 2.583 | 2.579 | 2.620 | |
平均值 | 0.618 68 | 16.377 | 0.052 27 | 5.549 | 2.543 | 2.611 | ||
B类 | LI231 | 2105.90 | 0.561 16 | 9.872 | 0.046 56 | 7.032 | 2.451 | 2.452 |
YAN56 | 3032.10 | 0.425 26 | 8.828 | 0.043 66 | 6.079 | 2.450 | 2.487 | |
YAN56 | 3058.40 | 0.742 62 | 10.210 | 0.051 74 | 7.032 | 2.445 | 2.438 | |
平均值 | 0.576 35 | 9.637 | 0.047 32 | 6.714 | 2.449 | 2.459 | ||
C类 | Z22 | 1563.40 | 0.203 17 | 7.641 | 0.038 32 | 12.120 | 2.528 | 2.458 |
Z22 | 1569.10 | 0.143 76 | 5.675 | 0.026 76 | 7.032 | 2.558 | 2.479 | |
LI68 | 2078.90 | 0.147 53 | 8.459 | 0.027 35 | 2.583 | 2.627 | 2.545 | |
平均值 | 0.164 82 | 7.258 | 0.030 81 | 7.245 | 2.571 | 2.494 |
总体上,长7页岩样品的滞后回环面积分布在0.143 76~0.742 62 mL/(nm·g),比表面积分布在5.675~20.210 m2/g,而比孔容和平均孔径则分别分布在0.026 86~0.060 20 mL/g和2.583~12.120 nm范围内。A类样品各项参数最优:滞后回环面积平均为0.618 60 mL/(nm·g),比表面积平均为16.377 m2/g,比孔容平均为0.052 27 mL/g,而平均孔径为5.549 nm;B类样品中等:滞后回环面积平均为0.576 35 mL/(nm·g),比表面积平均为9.637 m2/g,比孔容平均为0.047 32 mL/g,平均孔径为6.714 nm; C类样品最差:滞后回环面积平均为0.164 82 mL/(nm·g),比表面积平均为7.258 m2/g,比孔容平均为0.030 81 mL/g,平均孔径为7.245 nm。
3.3.2 孔径分布
DFT模型孔径分布解释结果显示,研究区长7段泥页岩所有样品的孔径主要分布在1.27~77.70 nm的范围内(图6)。A类样品中孔径小于40 nm的孔隙占主导地位,同时,其孔径分布曲线与横轴包络面积在3类样品中最大,指示样品的总孔隙度最大;B类样品虽然仍以孔径小于40 nm的孔隙为主,但孔径小于40 nm的孔径分布曲线与横轴的包络面积明显变小;C类样品孔径分布曲线总体十分平缓,主峰明显不如前两类样品突出,孔体积最小。
图6
图6
鄂尔多斯盆地陇东地区长7段泥页岩3类样品孔径分布
Fig. 6
The pore size distributions of samples in three types from the Chang 7 shale of Longdong area, Ordos Basin
3.4 FHH分形分析
本研究采用FHH分形模型,选取低温氮气吸附实验吸附-脱附曲线的脱附分支进行孔隙空间分形分析(图7)。双对数坐标中散点线性拟合的皮尔森相关系数均大于0.98,验证了分形模型的有效性;散点的分布呈现明显的分段线性特征(分界点相对压力p/p0约为0.5),说明不同孔径大小的孔隙空间存在不同的分形特征。Hazra等[41]认为,低压段(p/p0<0.5,D1)主要受范德华力控制,而高压段(p/p0>0.5,D2)则主要是界面张力起作用,毛细管凝聚效应突出。Xiong等[42]将D1称为孔隙表面分形维数,对吸附性能影响更大,主要受微孔控制;而将D2称为孔隙结构分形维数,对气体渗流和比孔容影响更大,主要受黏土矿物含量和热成熟度控制。总体上(表2),A类样品D1和D2均较大,分别为2.543和2.611,表明A型样品孔隙空间具有较强的非均质性,且孔隙结构非均质性强于孔隙表面非均质性;B类样品D1和D2均偏低(D1平均值为2.449, D2平均值为2.459),表明B类样品孔隙空间非均质性相对较弱;C类样品孔隙表面非均质性较强(D1平均值为2.571),而孔隙结构非均质性(D2平均值为2.494)较弱。
图7
图7
鄂尔多斯盆地陇东地区长7段泥页岩3类样品低温氮气吸附实验脱附分支分形分析
Fig. 7
Fractal analysis on the desorption branch of low-temperature nitrogen adsorption tests for Chang 7 shale samples in three types in Longdong area, Ordos Basin
3.5 FE-SEM镜下孔隙空间定性评价
扫描电镜观察结果显示,长7段泥页岩储层孔隙空间主要由有机孔、无机孔和微裂缝组成。其中,无机孔常见粒间孔、晶间孔和溶蚀孔3种类型(图8)。
图8
图8
鄂尔多斯盆地陇东地区长7段页岩样品孔隙空间扫描电镜定性观察
a.Z22井,埋深1 584.65 m;黏土矿物晶间孔、粒间孔发育,溶蚀孔发育,可见黄铁矿;b.Z22井,埋深1 626.00 m,有机质孔发育;c.Z22井,埋深1 632.45 m,黏土矿物晶间孔发育,可见黄铁矿;d.LI231井,埋深2 105.90 m,溶蚀孔发育,可见黏土矿物、黄铁矿、石英、长石;e.YAN56井,埋深3 032.10 m,有机质孔发育,可见黄铁矿;f.YA56井,埋深3 058.40 m,黄铁矿晶间孔发育,可见黏土;g.Z22井,埋深1 563.40 m,黏土矿物晶间孔、石英晶内孔、微裂缝发育;h.Z22井,埋深1 569.10 m,有机质孔发育;i.LI68井,埋深2 078.90 m,溶蚀孔发育Qtz.石英;F.长石;Clay.黏土矿物;Py.黄铁矿; OM.有机质;OP.有机质孔; PIP.晶间孔;IP.粒间孔;DP.溶蚀孔; Mf.微裂缝
Fig. 8
Qualitative observation on the pore structure of Chang 7 shale samples under scanning electron microscope (SEM) in Longdong area, Ordos Basin
1) 有机孔。有机孔为泥页岩中有机质在热裂解生烃过程中形成的孔隙,主要发育在有机质和无机矿物间或有机质内部,多呈圆形、椭圆形等(图8b,e,h)。
2) 粒间孔。粒间孔为矿物颗粒间残余的孔隙,多受控于压实作用(图8a)。
4)溶蚀孔。溶蚀孔为石英、碳酸盐、长石等不稳定矿物受地下水或干酪根生烃形成的酸性环境影响,发生化学溶解形成的孔隙,镜下多见溶蚀而成的不规则圆形(图8a,d,i)。
5)微裂缝。微裂缝主要有成岩作用、构造作用和有机质热演化局部压力异常3种成因。本研究区内微裂缝多为成岩收缩缝(图8g)和有机质收缩缝,其中成岩收缩缝连通性好,开张度变化大,而有机质收缩缝多分布于有机质与矿物颗粒的接触界面。
4 讨论
为探索陆相页岩在氮气吸附-脱附过程中能否形成“滞后回环”以及形成滞后回环的控制因素,本研究在扫描电镜定性观察的基础上,综合氮气吸附实验、岩石热解实验、TOC分析测试、XRD实验和分形分析,对形成滞后回环的一系列因素进行定量分析(图9)。在优选比表面积、比孔容、平均孔径、D1和D2、滞后回环面积、石英、黏土矿物、脆性矿物、TOC和Tmax等参数的基础上,引入“滞后回环类型”“是否有滞后回环”两个参数;针对参数“是否有滞后回环”,A、B两类页岩样品为True(是,值为1), C类样品为False(否,值为0);而后,使用Label-Encoding对类别参数,如“滞后回环类型”,进行编码;最后,以热度图的方式展示各参数之间两两相关分析的结果(数值为皮尔森相关系数Pearson’s r)。热度图(图9)显示:长7页岩样品在低温氮气吸附-脱附实验中能否形成滞后回环,与泥页岩样品的比表面积、比孔容和黏土矿物含量之间存在明显的正相关关系,而与总有机碳含量呈明显负相关关系,与平均孔径大小、分形维数(孔隙空间非均质性)、最高热解峰温、脆性矿物含量等无明显相关关系。
图9
图9
鄂尔多斯盆地陇东地区长7段页岩滞后回环参数相关分析热度图
Fig. 9
Heat-mapped correlation on the hysteresis-loop-related parameters for the Chang 7 shale in Longdong area, Ordos Basin
4.1 滞后回环的形成与类型
在气体吸附过程中会发生毛细凝聚现象:毛细管内流体的饱和蒸汽压比平液面小,于是毛细管内的液面上升,从而发生凝聚。如果样品含中孔和大孔,就一定会发生毛细管凝聚现象[43]。Foster和Cohan [44]解释了吸附-脱附在两端开口的毛细管中的差异:吸附阶段N2首先在孔隙内壁上形成一层非常薄的吸附膜;随着吸附量增加,吸附膜逐渐增厚直至形成较厚的凹透镜的液相;此时,弯月形液面的半径等于孔半径。而在脱附阶段,随着凹透镜状的液相逐渐变薄会形成两个隔着薄液膜的弯月形液面,这种弯月形液面在吸附过程中不会出现。这种过程性差异是形成滞后回环的重要原因,只有当吸附和脱附过程中形成的凝聚液半月膜形状或者直径明显不同时,才会观察到滞后现象[44]。
4.2 矿物含量及有机质对滞后回环面积的影响
图10
图10
滞后回环面积与石英(a)、黏土矿物(b)和脆性矿物含量(c)的关系
(Pearson’s r为皮尔森相关系数。)
Fig. 10
Correlation of hysteresis loop area with quartz (a), clay mineral (b) and brittle mineral (c) contents
A类页岩样品中黏土矿物含量最高且黏土矿物晶间孔占主导地位,B类次之,C类页岩样品中黏土矿物孔则不太发育。与之对应的,随着黏土矿物晶间孔发育程度的逐渐下降,A,B,C类泥页岩样品的滞后回环面积也逐渐变小(表2)。黏土矿物晶间孔一般呈平行板状和楔形,因此,当黏土矿物的含量越高、黏土矿物晶间孔越发育,滞后现象越明显、滞后回环面积越大,即:黏土矿物含量以及黏土矿物晶间孔发育程度共同决定了滞后回环的面积大小。
图11
图11
滞后回环面积与TOC(a)和Tmax(b)的关系
Fig.11
Correlation of hysteresis loop area with TOC (a) and Tmax (b)
4.3 非均质性对滞后回环面积的影响
长7段泥页岩样品的孔隙表面分形维数D1、孔隙结构分形维数D2与滞后回环面积之间基本不存在相关性(Pearson’s r分别为-0.42和0.15;图12a1,a2)。相比黏土矿物,泥页岩孔隙空间的非均质性受有机质孔隙,特别是有机孔中微孔的发育程度影响更大[53]。Tmax与非均质性之间存在着良好的正相关关系(Pearson’s r =0.77;图9),指示成熟度越高,有机质微孔越发育,非均质性越强。由于长7段泥页岩样品的热成熟度(Tmax)与有机质含量之间并无相关关系(Pearson’s r =-0.04;图9),因此,有机质含量的高低并不能决定有机孔的发育程度,有机质含量高低与非均质性之间也无相关关系。
图12
图12
滞后回环面积(a1,a2)、比表面积(b1,b2)、比孔容(c1,c2)与分形维数的关系
Fig.12
Correlation of hysteresis loop area(a1,a2), specific surface area(b1,b2), specific pore volume(c1,c2) with fractal dimension
4.4 孔隙结构参数对滞后回环面积的影响
图13
图13
滞后回环面积与比表面积(a)、比孔容(b)和孔径(c)的关系
Fig. 13
Correlation of hysteresis loop area with specific surface area (a), pore volume (b) and pore diameter (c)
5 结论
1) 长7段泥页岩在低温氮气吸附-脱附实验中能否形成滞后回环,归根到底取决于孔隙空间的主要形态:两端开放的圆柱形孔、墨水瓶孔或平行板孔占孔隙空间比例越高,滞后效应越强,滞后回环的面积越大。可以借助滞后回环面积来定量评价圆柱形孔、墨水瓶孔或平行板孔占孔隙空间的相对比例。
2) 长7段泥页岩氮气吸附滞后回环面积与黏土矿物含量之间存在明显的正相关关系:孔隙空间中两端开放的圆柱形孔、墨水瓶孔或平行板孔主要由黏土矿物晶间孔提供。
3) 长7段泥页岩氮气吸附滞后回环面积与比表面积、比孔容之间呈较强的正相关关系:比表面积、比孔容均主要由黏土矿物提供。
参考文献
Comparison of isotherm models for hydrocarbon adsorption on activated carbon
[J]. ,
Modeling CO2 transport and sorption in carbon slit pores
[J]. ,
Microscopic diffusion of CO2 in clay nanopores
[J]. ,
Pore size distribution from water adsorption isotherm
[C]//
鄂尔多斯盆地渭北地区页岩纳米孔隙发育特征及其控制因素
[J]. ,
Characteristics and controlling factors of nanopores in shales in Weibei, Ordos Basin
[J]. ,
渤海湾盆地东营凹陷沙河街组页岩油储层微观孔隙特征
[J]. ,
Micro-pore characteristics of shale oil reservoirs of the Shahejie Formation in the Dongying Sag, Bohai Bay Basin
[J]. ,
快速热解褐煤焦的低温氮吸附等温线形态分析
[J]. ,
The configuration analysis of the adsorption isotherm of nitrogen in low temperature with the lignite char produced under fast pyrolysis
[J]. ,
Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution (IUPAC technical report)
[J]. ,
On a theory of the van der Waals adsorption of gases
[J]. ,
基于低温氮吸附实验的页岩储层孔隙分形特征
[J]. ,
Fractal characteristics of shale reservoir pores based on nitrogen adsorption
[J]. ,
陆相页岩气储层孔隙发育特征及其主控因素分析: 以鄂尔多斯盆地长7段为例
[J]. ,
Characteristics of pore development and its main controlling factors of continental shale gas reservoirs: A case study of Chang 7 member in Ordos Basin
[J]. ,
高压压汞法和氮气吸附法分析页岩孔隙结构
[J]. ,
Pore structure of shales from high pressure mercury injection and nitrogen adsorption method
[J]. ,
Reporting physisorption data for gas/solid systems with special reference to the determination of surface area and porosity (recommendations 1984)
[J]. ,
Investigation of pore structure and fractal characteristics of organic-rich Yanchang formation shale in central China by nitrogen adsorption/desorption analysis
[J]. ,
Evaluation of pore size spectrum of gas shale reservoirs using low pressure nitrogen adsorption, gas expansion and mercury porosimetry: A case study from the Perth and Canning basins, Western Australia
[J]. ,
Physisorption hysteresis loops and the characterization of nanoporous materials
[J]. ,
煤微孔隙特征及其与油气运移储集关系的研究
[J]. ,
Study on the characteristics of coal microporosity and its relationship with oil and gas transportation and storage
[J]. ,
川东北地区下侏罗统自流井组陆相页岩储层孔隙结构特征
[J]. ,
Pore structure of continental shale reservoirs in Lower Jurassic Ziliujing Formation, northeastern Sichuan Basin
[J]. ,
川东南龙马溪组页岩微观孔隙结构特征
[J]. ,
Micro-pore structure of Longmaxi shale from Southeast Sichuan Basin
[J]. ,
南川地区龙马溪组优质页岩段微观孔隙结构特征
[J]. ,
Micro pore structure characteristics of high quality shale section of Longmaxi formation in Nanchuan area
[J]. ,
鄂尔多斯盆地古生界气藏成藏模式及优势储层预测
[J]. ,
Accumulation pattern and favorable reservoir prediction of Paleozoic gas reservoirs in Ordos Basin
[J]. ,
鄂尔多斯盆地东部构造演化对上古生界大气田形成的控制作用
[J]. ,
Effect of structural evolution on the formation and distribution of Upper Paleozoic giant gas fields in the East Ordos Basin
[J]. ,
鄂尔多斯盆地延长组长73亚段砂-泥协同成岩作用
[J]. ,
Synergistic diagenesis of sandstone and mudstone in Chang 73 sub-member of Triassic Yanchang Formation in Ordos Basin
[J]. ,
鄂尔多斯盆地三叠系延长组长7段页岩油勘探发现与资源潜力
[J]. ,
Discovery and resource potential of shale oil of Chang 7 member, Triassic Yanchang Formation, Ordos Basin
[J]. ,
Closed-system pyrolysis-based hydrocarbon generation simulation and gas potential evaluation of the Shanxi Formation shales from the Ordos Basin, China
[J]. ,
鄂尔多斯盆地东南部延长组长7段富有机质页岩孔隙特征
[J]. ,
Micropore characteristics of the organic-rich shale in the 7th member of the Yanchang Formation in the southeast of Ordos Basin
[J]. ,
鄂尔多斯盆地东南部长7段页岩孔隙特征及其控制因素
[J]. ,
Pore characteristics and controlling factors of Chang 7 shale in southeastern Ordos Basin
[J]. ,
鄂尔多斯盆地南缘铜川地区三叠系延长组长7段剖面及其油气地质意义
[J]. ,
Description and its hydrocarbon geological implications of outcrop sections of Triassic Chang-7 Member in southern Ordos Basin
[J]. ,
海陆过渡相页岩储层孔隙表征方法对比研究
[J/OL]. :
Comparative study on pore characterization methods of shale reservoir
[J].
泸州地区龙马溪组深层页岩孔隙结构特征
[J]. ,
The pore structure characteristics of deep shale in Longmaxi Formation of Luzhou area
[J]. ,
A novel hybrid thermodynamic model for pore size distribution characterisation for shale
[J]. ,
川南下志留统龙马溪组页岩岩相表征
[J]. ,
Lithofacies characterization of Longmaxi Formation of the Lower Silurian, southern Sichuan
[J]. ,
Distribution characteristics and oil mobility thresholds in lacustrine shale reservoir: Insights from N2 adsorption experiments on samples prior to and following hydrocarbon extraction
[J]. ,
Lower Cretaceous gas shales in northeastern British Columbia, Part I: Geological controls on methane sorption capacity
[J]. ,
对BJH方法计算孔径分布过程的解读
[J]. ,
Interpretation of BJH method for calculating aperture distribution process
[J]. ,
Physical adsorption characterization of nanoporous materials
[J]. ,
Porosity controls and fractal disposition of organic-rich Permian shales using low-pressure adsorption techniques
[J]. ,
An investigation of fractal characteristics of marine shales in the southern China from nitrogen adsorption data
[J]. ,
采用机器学习分割算法和扫描电镜分析页岩微观孔隙结构
[J]. ,
Pore structure characterization of shales using SEM and machine learning-based segmentation method
[J]. ,
Spectrum of pore types and networks in mudrocks and a descriptive classification for matrix-related mudrock pores
[J]. ,
四川盆地元坝地区自流井组页岩储层孔隙结构特征
[J]. ,
Pore structure characteristics of shale reservoirs in the Ziliujing Formation in Yuanba area, Sichuan Basin
[J]. ,
陆相页岩油储层微观孔喉结构表征与含油性分级评价——以准噶尔盆地吉木萨尔凹陷二叠系芦草沟组为例
[J]. ,
Microscopic pore structure characterization and oil-bearing property evaluation of lacustrine shale reservoir: A case study of the Permian Lucaogou Formation in Jimsar Sag, Junggar Basin
[J]. ,
页岩中有机质孔隙非均质性的微观结构及电镜—拉曼联用研究
[J]. ,
Microstructure and SEM-Raman study of organic matter pore heterogeneity in shale
[J]. ,
Comparative study on micro-pore structure of marine, terrestrial, and transitional shales in key areas, China
[J]. ,
四川盆地东北部下侏罗统自流井组陆相页岩储层孔隙特征及形成控制因素
[J]. ,
Pore characteristics and controlling factors of continental shale reservoirs in the Lower Jurassic Ziliujing Formation, northeastern Sichuan Basin
[J]. ,
中国南方典型页岩孔隙特征差异及其控制因素
[J]. ,
Characteristic differences and controlling factors of pores in typical South China shale
[J]. ,
Heterogeneity of intergranular, intraparticle and organic pores in Longmaxi shale in Sichuan Basin, South China: Evidence from SEM digital images and fractal and multifractal geometries
[J]. ,
/
〈 | 〉 |