川东北地区侏罗系凉高山组页岩储层孔隙结构及分形特征

李琦; 陈睿倩; 商斐; 李玲; 白昕

doi:10.11781/sysydz2025020323

川东北地区侏罗系凉高山组页岩储层孔隙结构及分形特征

doi: 10.11781/sysydz2025020323

李琦^{1, 2,},
陈睿倩^{1, 2, ,},
商斐³,
李玲^{1, 2},
白昕^{1, 2}

1.
中国石油大学(北京)油气资源与工程全国重点实验室, 北京 102249
2.
中国石油大学(北京) 地球科学学院, 北京 102249
3.
中国石油勘探开发研究院, 北京 100083

基金项目:

中国石油西南油气田分公司科技计划项目“川东北地区侏罗系下统凉高山组页岩油气地质特征研究” JS2022-90

详细信息

作者简介:
李琦(2000—)，女，硕士生，从事页岩储层评价相关研究。E-mail: 270331500@qq.com

通讯作者:
陈睿倩(1987—)，女，博士，副教授，从事石油地质学和地球化学研究。E-mail: richen@cup.edu.cn

中图分类号: TE122.2
计量
- 文章访问数: 176
- HTML全文浏览量: 64
- PDF下载量: 35
- 被引次数: 0
出版历程
- 收稿日期: 2024-07-02
- 修回日期: 2025-02-10
- 刊出日期: 2025-03-28

Pore structure and fractal characteristics of shale reservoirs in Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

LI Qi^{1, 2
,},
CHEN Ruiqian^{1, 2
, ,},
SHANG Fei³,
LI Ling^{1, 2},
BAI Xin^{1, 2}

1.
State Key Laboratory of Petroleum Resources and Engineering, China University of Petroleum (Beijing), Beijing 102249, China
2.
College of Geosciences, China University of Petroleum (Beijing), Beijing 102249, China
3.
Research Institute of Petroleum Exploration and Development, CNPC, Beijing 100083

摘要

摘要: 川东北地区侏罗系凉高山组作为四川盆地页岩油重点勘探层段，因勘探程度较低，目前页岩储层特征尚不明确。利用X射线衍射矿物分析、扫描电子显微镜分析、高压压汞、低温氮气吸附等实验，深入探究其储集空间类型与分形特征。研究区凉高山组页岩矿物组成主要为黏土矿物(平均含量51.57%)，其次为长英质矿物(平均含量47.11%)，碳酸盐矿物极少(平均含量2.69%)。储集空间类型主要包括黏土矿物层间孔、石英与长石粒间孔以及微裂缝，低温氮气吸附曲线形态接近国际理论与应用化学联合会页岩孔径划分方法中的Ⅳ类，指示狭缝型孔隙。按照压汞曲线形态及储层物性参数特征，将研究区页岩储层分为四类，从Ⅰ到Ⅳ类，排驱压力、中值压力增大，最大进汞饱和度减小，储层非均质性增强。“FHH”模型显示孔隙表面分形维数D_N1大于孔隙结构分形维数D_N2，表明孔隙表面的复杂程度高于内部；“含水饱和度法”计算结果显示大孔隙分形维数D₁平均值为2.991 2，小孔隙分形维数D₂平均值为2.679 2，大孔隙分形维数更接近于3且分布更集中，表明非均质性更强的大孔隙对储层的贡献更大。相关性分析显示，D₁与矿物组分含量(石英、黏土矿物)、孔喉结构参数存在相关性，证明大孔隙是研究区页岩储集与空间的主要贡献者。通过定性与定量分析对川东北地区凉高山组页岩进行储层评价，为后续该地区有利勘探层段评价与优选提供了思路。
- 页岩 /
- 低温氮气吸附 /
- 高压压汞 /
- 分形维数 /
- 储层非均质性 /
- 凉高山组 /
- 四川盆地
Abstract: The Jurassic Lianggaoshan Formation in the northeastern Sichuan Basin is a key exploration target for shale oil. However, due to limited exploration in this area, the shale reservoir characteristics remain unclear. Experiments such as X-ray diffraction mineral analysis, scanning electron microscopy analysis, high-pressure mercury intrusion, and low-temperature nitrogen adsorption were conducted to systematically study the storage space types and fractal features of the Lianggaoshan Formation shale reservoir. The primary mineral composition of the Lianggaoshan Formation reservoir in the northeastern Sichuan Basin is clay minerals, with an average content of 51.57%, followed by feldspar and quartz minerals at an average of 47.11%, while carbonate minerals are scarce, averaging 2.69%. The dominant storage space types mainly include interlayer pores of clay minerals, intergranular pores between quartz and feldspar, and micro-fractures. The low-temperature nitrogen adsorption curve of the Lianggaoshan Formation shale aligns with type Ⅳ in the classification system of the International Union of Pure and Applied Chemistry, indicating slit-type pores. Based on the morphology of mercury intrusion curves and reservoir physical properties, the reservoir is divided into four types. From type Ⅰ to type Ⅳ, drainage pressure and median pressure increase, whereas maximum mercury saturation decreases, leading to enhanced reservoir heterogeneity. The "FHH"model calculations show that the pore surface fractal dimension (D_N1) is greater than the pore structure fractal dimension (D_N2), indicating that the pore surface exhibits greater complexity than the internal pore structure. The average fractal dimension D₁ of large pores, calculated using the water saturation method, averages 2.991 2, while that of small pores (D₂) averages 2.679 2. The larger pores have a fractal dimension closer to 3 and exhibit a more concentrated distribution, indicating that highly heterogeneous large pores contribute more significantly to the reservoir. Correlation analysis shows that there is a correlation between D and the contents of minerals (quartz and clay minerals) as well as pore-throat struture parameters, proving that large pores are the main contributers to the shale reservoir space in the study area. Through qualitative and quantitative analyses, this paper conducts a reservoir evaluation of the Lianggaoshan Formation shale in the northeastern Sichuan Basin, offering insights for the subsequent evaluation and selection of favorable exploration intervals in this area.
- shale /
- low-temperature nitrogen adsorption /
- high-pressure mercury intrusion /
- fractal dimension /
- reservoir heterogeneity /
- Lianggaoshan Formation /
- Sichuan Basin
注释:

利益冲突声明/Conflict of Interests

作者陈睿倩是本刊青年编委会成员，未参与本文的同行评审或决策。

Author CHEN Ruiqian is a Young Editorial Board Member of this journal, and she did not take part in peer review or decision making of this article.

作者贡献/Authors’Contributions

李琦、陈睿倩参与论文写作和修改；李琦、李玲、白昕完成实验操作与数据分析；商斐提供样品与资料支持。所有作者均阅读并同意最终稿件的提交。

The manuscript was drafted and revised by LI Qi and CHEN Ruiqian. The experimental operation and data analysis were completed by LI Qi, LI Ling, and BAI Xin, while SHANG Fei provided sample and data support. All authors have read the final version of the paper and consented to its submission.

HTML全文

四川盆地是我国西南地区重要的含油气盆地，本研究的目的层位于川东北地区侏罗系凉高山组，也是该地区的重点勘探层位。前人对于四川盆地的勘探多集中在川中、川南地区，并取得重大勘探成果。川东北地区侏罗系由于埋深较浅、岩性致密，长期以来都被视作“过路层”^[1-2]，研究较少。近年来，川东仪陇—平昌地区平安1井、川北巴中地区巴中1HF井均获日产油量万余吨、日产气量百余立方米，稳产形势良好，展示了侏罗系凉高山组巨大的勘探前景^[3-5]。前人已对凉高山组页岩的沉积环境和层序特征等方面进行了研究，然而对于该页岩储层的空间类型和定量表征尚未明确。

传统基于欧氏几何的孔隙表征方法难以实现页岩储层非均质性的定量评价，而分形理论通过分形维数这一标度参数，可有效描述孔隙表面的粗糙度与孔隙结构的复杂程度，为页岩储集空间定量表征提供了新的理论框架^[6]。近年来，该理论已在页岩储层研究中展现出独特的应用优势^[7-8]。李小明等^[7]通过低温氮气吸附实验揭示出湖北荆门探区五峰组—龙马溪组页岩的分形维数与总有机碳(TOC)含量呈显著正相关，而与黏土矿物含量及平均孔径存在显著负相关关系；王斌等^[8]利用低温氮气吸附实验与核磁共振技术相结合，发现川东南石牛栏组页岩不同尺度孔隙分形特征存在显著差异，其中微孔(< 2 nm)与介孔(2~50 nm)分形维数与TOC含量呈正相关，与脆性矿物含量呈负相关；而宏孔(>50 nm)分形维数则主要受脆性矿物含量控制，与TOC含量无显著统计关联。

本次研究通过X射线衍射(XRD)矿物分析、扫描电子显微镜(SEM)分析、高压压汞、低温氮气吸附等实验方法，明确研究区页岩储集空间类型、计算储层分形维数，进而实现储层的定量表征，以期为后续研究区有利储层评价与优选提供思路。

1. 区域地质背景

四川盆地形成于印支期，由喜马拉雅综合褶皱运动形成，主要划分为六个构造区^[9](图 1)。研究区位于川东高陡构造区，受印支、燕山、喜马拉雅三期构造运动影响，发育一系列北东向大型隔挡式高陡构造。其中，背斜狭窄，核部主要为三叠系地层；向斜宽缓，核部出露中侏罗统地层^[10-12]。根据岩性、测井曲线和层序地层特征，将研究区凉高山组分为凉一段、凉二段及凉三段(图 2)。其中，凉一段主要发育滨浅湖亚相沉积，岩性以一套厚层深灰色泥页岩夹粉砂岩为主；凉二段主要发育滨浅湖亚相、三角洲前缘亚相以及深湖亚相沉积，表明凉二段沉积期发生了一次大规模海侵运动，岩性以泥页岩与粉砂岩互层为主^[13]；凉三段主要发育曲流河三角洲亚相沉积，岩性以黑色泥页岩与细砂岩互层为主。

图 1 四川盆地构造分区及研究区位置

据参考文献[11]修改。

Figure 1. Tectonic divisions of Sichuan Basin and location of study area

下载: 全尺寸图片幻灯片

图 2 川东北侏罗系凉高山组地层柱状图

据参考文献[11]修改。

Figure 2. Stratigraphic column of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

2. 样品及实验方法

本研究共采集样品31块，包括会浅1井19块、中兴1井1块以及五宝浅020-H₁井11块(图 1)。所有样品均开展XRD矿物分析、铸体薄片鉴定以及SEM分析，同时选择其中16块样品进行高压压汞实验以及低温氮气吸附实验。XRD矿物分析分别测定全岩矿物及黏土矿物含量，采用荷兰帕纳科X射线衍射仪X’Pert MPD PRO衍射仪，测试标准依据《矿物定量分析X射线衍射法：GB/T 30903—2014》；岩石薄片鉴定采用光学显微镜MIN-017，放大倍数在40~1000倍之间，测试标准依据《岩石薄片鉴定方法：GB/T 17412.1—1998》；SEM分析采用蔡司Merlin仪器，加速电压5~10 kV，测试标准依据《扫描电子显微镜分析方法通则：GB/T 17359— 2012》；高压压汞实验采用全自动压泵仪PoreMaster 60，测试标准依据《压汞法测定岩石孔隙结构：GB/T 21650.1—2008》；低温氮气吸附实验采用Micromeritics ASAP 2460比表面积和孔径分析仪，测试标准依据《多孔材料比表面积和孔隙结构的测定：GB/T 21650.2—2008》。以上实验均在中国石油大学(北京)油气资源与工程全国重点实验室完成。

3. 储层特征

3.1 岩石学特征

研究区凉高山组页岩矿物组成主要为黏土矿物，平均含量为51.57%；其次为石英，平均含量为40.85%；长石平均含量为6.26%；碳酸盐矿物(包括方解石、白云石及文石等)含量最低，平均为2.69%；部分样品还含有少量的黄铁矿(图 3)。研究区会浅1井、中兴1井、五宝浅020-H₁井等3口井样品的矿物组成具有相似性，黏土矿物与长英质矿物含量都主要分布在25%~75%，碳酸盐矿物含量最低，分布在0~25%(图 4)。

图 3 川东北侏罗系凉高山组矿物组成条形图

Figure 3. Bar chart of mineral composition in Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 4 川东北侏罗系凉高山组矿物组成三角图

Figure 4. Ternary diagram of mineral composition in Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

3.2 储集空间类型

铸体薄片鉴定及SEM分析结果显示，研究区凉高山组页岩储集空间类型主要包括无机孔、有机孔及微裂缝。

3.2.1 无机孔

根据油气赋存位置，研究区无机孔隙又可进一步划分为粒间孔、粒内孔以及晶间孔。粒间孔主要包括黏土矿物层间孔、石英与长石粒间孔，孔径在几百纳米到几微米之间。研究区黏土矿物层间孔多呈狭缝型(图 5a)，由于黏土矿物塑性强、抗压能力弱，网状层间孔不易保存，因此研究区网状黏土矿物网状层间孔较少(图 5b)^[14]。黏土矿物层间孔含量多、比表面积大，有利于提高储层油气的渗流能力。长石与石英粒间孔孔隙大、数量多(图 5c)，半径可达50~100 μm，为油气储集贡献了主要空间。

图 5 川东北侏罗系凉高山组页岩储集空间类型

a.狭缝状黏土矿物层间孔，会浅1井，1 764.20 m；b.网状黏土矿物层间孔，五宝浅020-H₁井，3 183.66 m；c.石英与长石粒间孔，可见缩颈型喉道，会浅1井，1 898.00 m；d.草莓状黄铁矿，会浅1井，1 765.28 m；e.长石内溶蚀孔，会浅1井，1 765.28 m；f.岩屑溶蚀孔，五宝浅020-H₁井，3 177.36 m；g.块状有机质，可见有机质收缩缝，中兴1井，1 901.35 m；h.凹坑状有机质孔，会浅1井，1 902.01 m；i.条带状有机质，可见狭缝状有机质孔，五宝浅020-H₁井，3 187.00 m；j.构造微裂缝，可见黄铁矿晶体，中兴1井，1 905.48 m；k.构造微裂缝，五宝浅020-H₁井，3 180.24 m；l.有机质收缩缝，可见层状黏土矿物，中兴1井，1 897.65 m。

Figure 5. Types of reservoir space in shale from Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

研究区页岩晶间孔主要为黄铁矿晶间孔，与有机质的发育密切相关。研究区大部分黄铁矿呈草莓状(图 5d)，少部分呈他形。

研究区粒内孔主要为矿物内部溶蚀孔，包括长石粒内溶孔及岩屑粒内溶孔(图 5e-f)，为不稳定矿物颗粒(如长石、黏土矿物等)内部遭到溶蚀而形成的孔隙，一般连通性较弱，对储层渗流能力影响不大。

3.2.2 有机孔

有机质孔分布受有机质丰度、有机质类型以及固体沥青等因素的影响。研究区有机质丰度较高，有机质孔较为发育，但由于孔隙过小，因此对储层的贡献不大。有机质孔包括因生烃作用而产生的有机质热成因孔以及生烃作用导致有机质体积收缩而产生的有机质收缩缝^[15]。由于研究区凉高山组烃源岩成熟度较高，因此有机孔多为生烃完成后有机质收缩而产生的有机质收缩缝，多在黏土矿物与有机质的交界处(图 5g)，长度可达5~10 μm。此外，由于凉高山组页岩干酪根多为Ⅲ型，以产气为主，因此研究区有机质热成因孔多为纳米级气孔，形状为狭长状、凹坑状(图 5h-i)。

3.2.3 微裂缝

微裂缝作为连接储层微观孔隙与宏观裂缝的通道，对油气存储和渗流都具有重要意义^[16]。研究区主要发育构造裂缝以及成岩微裂缝，其中构造裂缝主要受成岩作用、矿物脆性特征、地层压力及构造作用控制。由于凉高山组储层脆性矿物含量高以及川东北地区位于褶皱作用形成的高陡构造带上，因此构造裂缝在研究区内最为常见，长度达几十至数百微米(图 5j-k)。成岩微裂缝与成岩过程中的脱水作用有关，研究区凉高山组储层主要发育有机质收缩缝(图 5l)。

3.3 物性特征

由于研究区凉高山组页岩样品含有较多微裂缝，抗压实能力弱，因此本研究采用高压压汞法获取储层物性。实验结果表明(表 1)，16块样品孔隙度介于1.96%~11.45%，渗透率介于(0.015 3~3.33)× 10^-3 μm²；排驱压力小，介于0.182 2~1.838 1 MPa，表明孔隙半径越大，岩石渗透性越好；最大进汞饱和度大，介于83.02%~91.38%，表明残余油饱和度小。

表 1 川东北侏罗系凉高山组页岩基本物性及孔喉结构参数

Table 1. Basic physical properties and pore-throat structure parameters of shale from Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

样品编号	分类	井深/m	孔隙度/%	渗透率/10^-3 μm²	排驱压力/MPa	中值压力/MPa	最大进汞饱和度/%	分选系数
WB-5	Ⅰ	3 184.66	4.81	0.003 01	0.182 6	1.453 6	90.73	2.798 2
HQ1-2	Ⅰ	1 765.28	2.47	0.060 80	0.292 0	2.135 0	89.89	2.672 6
HQ1-3	Ⅰ	1 889.17	2.30	3.330 00	0.191 6	3.106 2	88.79	2.923 3
WB-4	Ⅰ	3 183.66	4.35	0.086 90	0.182 2	1.566 1	87.75	2.738 1
WB-6	Ⅱ	3 187.00	3.54	0.050 10	0.463 4	3.554 6	91.38	2.752 3
WB-1	Ⅱ	3 176.10	2.10	0.052 00	0.463 9	4.415 2	88.33	2.944 2
HQ1-1	Ⅱ	1 764.20	1.96	0.048 90	0.459 2	3.899 9	91.44	2.623 0
WB-2	Ⅱ	3 177.36	8.08	0.038 00	0.458 9	4.271 6	88.45	2.763 3
WB-3	Ⅱ	3 179.54	8.88	0.033 00	0.457 2	2.056 9	91.73	2.348 5
HQ1-8	Ⅲ	1 902.01	6.60	0.013 20	0.730 0	10.150 4	87.08	2.911 0
HQ1-7	Ⅲ	1 900.39	7.80	0.005 44	0.734 9	22.598 3	87.62	2.874 1
ZX1-1	Ⅳ	1 901.35	11.45	0.014 40	1.135 5	12.969 6	83.02	2.870 0
HQ1-4	Ⅳ	1 894.61	2.15	0.004 69	1.158 8	13.177 9	85.87	2.718 2
HQ1-5	Ⅳ	1 897.15	3.97	0.006 29	1.158 9	4.298 5	87.05	2.494 2
HQ1-9	Ⅳ	1 915.01	5.79	0.006 65	1.160 1	47.529 0	84.67	2.984 1
HQ1-6	Ⅳ	1 898.00	10.19	0.001 53	1.838 1	11.855 9	83.71	2.622 5

下载: 导出CSV

| 显示表格

根据压汞曲线形态及储层物性参数特征，将研究区凉高山组储层分为Ⅰ、Ⅱ、Ⅲ、Ⅳ四类。Ⅰ类储层压汞曲线呈现“低、长、平缓”的特征(图 6a)，排驱压力小，不超过0.3 MPa，储层物性好；Ⅱ类压汞曲线形态与Ⅰ类储层类似(图 6b)，但排驱压力稍大，介于0.4~0.5 MPa，最大进汞饱和度较Ⅰ类储层更大，储层物性较好；Ⅲ类储层压汞曲线呈现“较高、较陡”的特点(图 6c)，没有明显的平台段，排驱压力介于0.7~0.8 MPa，分选系数较大，储层非均质性较强；Ⅳ类储层压汞曲线呈现“最高、最陡、平缓段短”的特点(图 6d)，排驱压力均大于1 MPa，最大进汞饱和度是四类储层中最小的，储层物性最差。

图 6 川东北侏罗系凉高山组四类储层压汞曲线

Figure 6. High-pressure mercury intrusion curves of four types of reservoirs in Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

4. 储层孔隙结构

根据国际理论与应用化学联合会(IUPAC)页岩孔径划分方法^[17]，本文将研究区储层孔隙划分为微孔(< 2 nm)、介孔(2~50 nm)和宏孔(>50 nm)，并对样品全孔径进行定量分析。由于氮气吸附曲线在尾端压力达1 MPa时也未呈现吸附饱和现象，证明样品中存在宏孔。因此，本研究利用低温氮气吸附实验表征微孔和介孔的分形维数，利用高压压汞实验表征宏孔的分形维数。

4.1 基于“FHH”模型计算氮气吸附曲线分形维数

4.1.1 氮气吸附脱附曲线

本研究16块页岩样品的氮气吸附曲线整体上均呈现反“S”型(图 7)，与IUPAC分类中的Ⅳ类相符合，存在明显的吸附台阶，表明页岩样品主要发育介孔，并存在大量狭缝型孔隙。在相对压力(P/P₀)＜0.45时，曲线呈上凸缓慢增长，代表该阶段是由范德华力引起的单层氮气吸附；在0.45＜P/P₀＜0.9时，曲线呈下凹型迅速上升，指示该阶段是由毛管凝聚力引起的多层氮气吸附，并出现回滞环；在P/P₀＞0.9时，曲线急速上升，直到P/P₀达到1时也未出现吸附饱和现象，证明页岩样品中存在部分宏孔。从Ⅰ类到Ⅳ类储层，尾端吸附量逐渐降低，表明页岩储层物性逐渐变差(图 7)。

图 7 川东北侏罗系凉高山组页岩低温氮气吸附脱附曲线

Figure 7. Low-temperature nitrogen adsorption and desorption curves of Jurassic Lianggaoshan Formation shale, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

4.1.2 “FHH”模型分形维数

本研究采用Pfeifer P的Frenkel-Halsey-Hill(FHH)模型^[18]计算氮气吸附曲线分形维数。FHH模型中吸附于分形表面上的气体表示为：

$\ln V=K \ln \left(\ln P / P_0\right)+\mathrm{C}$

(1)

式中: V为平衡压力P对应的吸附体积，单位cm³/g；P为平衡压力，单位MPa；P₀为吸附气体的饱和蒸汽压力，单位MPa；K为线性关系系数，其值与吸附机理有关；C为常数。分形维数D=K+3，D介于2~3之间。其中，D接近2表明孔隙光滑、结构简单，D接近3表明孔隙结构复杂、非均质性强。

从样品相对压力与吸附量拟合的双对数曲线(图 8)来看，分形维数拟合效果较好，同时样品均在P/P₀=0.45左右出现滞后环，因此以P/P₀=0.45为界将分形维数分为两个不同区域。其中，P/P₀＜0.45为低比压区，是以范德华力主导的单层分子吸附，代表了孔隙表面的分形维数，记为D_N1；P/P₀＞0.45为高比压区，是以毛细管凝聚力主导的多层分子吸附，代表了孔隙结构的分形维数，记为D_N2。

图 8 川东北侏罗系凉高山组四类储层氮气吸附分形结果

Figure 8. Fractal analysis results of nitrogen adsorption in four types of reservoirs of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

结合分形维数计算结果(表 2)来看，D_N1在2.485 8~2.712 6之间，平均值为2.654 8；D_N2在2.552 9~2.789 5之间，平均值为2.615 6。其中，Ⅰ类储层D_N1均值为2.607 1，D_N2均值为2.591 0；Ⅱ类储层D_N1均值为2.649 1，D_N2均值为2.659 6；Ⅲ类储层D_N1均值为2.698 6，D_N2均值为2.596 6；Ⅳ类储层D_N1均值为2.681 2，D_N2均值为2.598 9。从Ⅰ到Ⅳ类储层，物性逐渐变差，其中Ⅲ类储层物性最差。Ⅱ类储层孔隙表面复杂程度大于孔隙内部结构(D_N1＞D_N2)，推测可能是因为Ⅱ类储层TOC含量较高，有机质孔数量多，从而导致孔隙表面复杂程度更强。

表 2 川东北侏罗系凉高山组页岩氮气吸附分形维数计算结果

Table 2. Fractal dimension calculation of nitrogen adsorption in shale from Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

样品编号	分类	D_N1	R₁²	D_N2	R₂²
WB-5	Ⅰ	2.640 7	0.997 1	2.599 9	0.999 7
HQ1-2	Ⅰ	2.649 3	0.997 9	2.601 9	0.999 0
HQ1-3	Ⅰ	2.652 7	0.986 7	2.584 0	0.996 0
WB-4	Ⅰ	2.485 8	0.991 3	2.578 1	0.992 0
WB-6	Ⅱ	2.607 6	0.994 6	2.599 9	0.952 4
WB-1	Ⅱ	2.665 6	0.994 4	2.650 3	0.996 6
HQ1-1	Ⅱ	2.667 9	0.994 7	2.658 3	0.988 0
WB-2	Ⅱ	2.610 1	0.994 8	2.789 5	0.961 5
WB-3	Ⅱ	2.694 1	0.993 2	2.599 8	0.999 3
HQ1-7	Ⅲ	2.712 6	0.994 6	2.589 3	0.998 5
HQ1-8	Ⅲ	2.684 5	0.990 7	2.603 9	0.999 4
HQ1-6	Ⅳ	2.703 2	0.990 8	2.613 6	0.999 3
HQ1-5	Ⅳ	2.693 0	0.992 2	2.614 4	0.999 0
HQ1-4	Ⅳ	2.665 4	0.998 5	2.620 6	0.999 5
ZX1-1	Ⅳ	2.682 4	0.992 4	2.552 9	0.993 8
HQ1-9	Ⅳ	2.662 1	0.996 4	2.592 8	0.999 4

下载: 导出CSV

| 显示表格

4.2 基于“含水饱和度法”计算压汞曲线分形维数

本研究采用贺承祖等^[19-20]基于高压压汞曲线创建的含水饱和度法计算分形维数，具体公式如下：

$S_{\mathrm{w}}=1-S_{\mathrm{Hg}}$

(2)

$S_{\mathrm{w}}=\left(P_{\mathrm{c}} / P_{\min }\right)^{D-3}$

(3)

两边同时取对数：

$\lg S_{\mathrm{w}}=(3-D) \lg P_{\min }+(D-3) \lg P_{\mathrm{c}}$

(4)

式中：P_c为半径为r的孔隙对应的毛细管压力，单位MPa；P_min为最大孔喉半径r_max对应的毛细管压力，单位MPa；S_Hg为压力P_c时所对应的进汞饱和度，单位%；S_w为压力P_c对应的润湿相饱和度，单位%；D为分形维数。

由式(4)可得lgS_w与lgP_c呈线性关系，斜率k为D-3，因此D=k+3。

根据上述分形维数计算公式，对本次高压压汞实验数据进行处理。结果显示，lgS_w与lgP_c拟合的曲线具有明显的分段性，即分形曲线存在明显的的转折点(图 9)。其中位于转折点压力前近乎水平段，代表大孔隙的分形维数，记为D₁；位于转折点压力后倾斜段，代表小孔隙分形维数，记为D₂。

图 9 川东北侏罗系凉高山组四类储层高压压汞分形结果

Figure 9. Fractal results based on high-pressure mercury intrusion in four types of reservoirs of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

分形维数计算结果(表 3)显示，全部样品D₁平均值为2.991 2、D₂平均值为2.679 2，D₁均大于D₂且更接近3，表明大孔隙的非均质性更强、孔喉结构更复杂。研究区页岩大孔隙类型复杂多变，包括粒间孔、粒内孔、晶间孔以及微裂缝等，同时在压实作用后期容易被黏土矿物充填导致结构更复杂、非均质性更强。此外，大孔隙分形维数分布相对更集中(表 3)，表明储层大孔隙更为发育，对页岩储集性的贡献也相对更大。因此，本研究认为大孔隙是控制凉高山组页岩储集性的主要因素。

表 3 川东北侏罗系凉高山组页岩高压压汞分形维数计算结果

Table 3. Fractal dimension calculation based on high-pressure mercury intrusion in shale from Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

样品编号	分类	D₁	R₁²	D₂	R₂²
WB-5	Ⅰ	2.995 6	0.327 9	2.694 3	0.972 6
HQ1-2	Ⅰ	2.971 0	0.292 5	2.660 8	0.971 8
HQ1-3	Ⅰ	2.976 8	0.454 9	2.677 0	0.966 5
WB-4	Ⅰ	2.983 7	0.380 8	2.719 1	0.992 0
WB-6	Ⅱ	2.996 7	0.407 4	2.685 8	0.981 8
WB-1	Ⅱ	2.984 2	0.265 5	2.700 3	0.918 3
HQ1-1	Ⅱ	2.996 9	0.392 4	2.619 8	0.972 9
WB-2	Ⅱ	2.987 9	0.270 7	2.675 3	0.971 4
WB-3	Ⅱ	2.993 2	0.298 2	2.617 3	0.992 2
HQ1-7	Ⅲ	2.989 9	0.642 9	2.722 1	0.976 3
HQ1-8	Ⅲ	2.999 8	0.230 4	2.746 1	0.965 2
HQ1-6	Ⅳ	2.996 3	0.188 7	2.631 2	0.983 7
HQ1-5	Ⅳ	2.998 7	0.284 2	2.665 9	0.985 9
HQ1-4	Ⅳ	2.999 8	0.231 9	2.664 9	0.955 1
ZX1-1	Ⅳ	2.988 0	0.279 9	2.687 5	0.926 5
HQ1-9	Ⅳ	2.998 3	0.351 2	2.699 1	0.826 3

下载: 导出CSV

| 显示表格

具体来看(表 3)，Ⅰ类储层分形曲线转折点压力为0.458 8 MPa，D₁介于2.971 0~2.995 6之间，D₂介于2.660 8~2.719 1之间(图 9a)；Ⅱ类储层分形曲线转折点压力为0.623 9 MPa，D₁在2.984 2~2.996 9之间，D₂在2.617 3~2.700 3之间(图 9b)；Ⅲ类储层分形曲线转折点压力为0.730 0 MPa，D₁大于D₂(图 9c)；Ⅳ类储层分形曲线转折点压力为1.919 6 MPa，D₁介于2.988 0~2.999 8之间，D₂介于2.631 2~2.699 1之间(图 9d)。从Ⅰ类到Ⅳ类储层，分形曲线转折点处的毛细管压力增大、排驱压力增大，表明孔喉半径减小、孔隙结构复杂性增加、储层非均质性增强(图 9)。整体而言，凉高山组页岩储层孔隙分形维数均大于2.5，表明孔隙结构复杂、非均质性强、储层物性差。

5. 分形维数的影响因素

前人研究表明页岩储层分形维数受多方面因素影响，包括TOC含量、矿物组分、储层物性、孔喉结构参数等^[21]。本研究分别讨论矿物组分、储层物性、孔喉参数以及TOC含量对研究区侏罗系凉高山组页岩储层分形维数的影响。

5.1 分形维数与矿物组分的关系

分形维数与矿物组分含量相关性分析结果表明，D₁与石英、黏土矿物含量拟合关系较好，相关系数R²分别为0.479 3与0.492 4。其中，D₁与石英含量呈负相关(图 10a)。石英抗风化能力强，表面光滑平整、不易被溶蚀，石英形成的微米级晶间孔以及不同种类型的喉道构成了凉高山组储层的主要储集空间。相反地，D₁与黏土矿物含量呈正相关(图 10b)。黏土矿物塑性强、抗压实能力弱，形成的孔隙不易保存。同时，黏土矿物形成的层间孔与网状孔连通性、储集性均较差。此外，黏土矿物中的蒙脱石遇水还会膨胀，进一步阻碍储层孔隙间的连通。D₂与石英、黏土矿物含量均无明显相关关系(图 10c-d)，证明大孔隙是研究区侏罗系凉高山组页岩储集空间的主要贡献者。

图 10 川东北侏罗系凉高山组矿物组分含量与分形维数关系

Figure 10. Correlation between mineral composition and fractal dimension of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

5.2 分形维数与储层物性的关系

分形维数与储层物性相关性分析结果表明，D₁与孔隙度几乎没有相关性(图 11a)，与渗透率呈弱的负相关性，R²为0.397 5(图 11b)。渗透率低表明储层孔隙连通性及渗流能力差，因此储层非均质性增强、分形维数增加。D₁与储层物性相关性差的原因主要包括：含水饱和度法计算分形维数的前提为R_max＞＞R_min，研究区凉高山组页岩的孔喉半径介于153.125 0~0.003 7 μm(平均17.660 0 μm)，基本满足这一前提条件，但含水饱和度法所求的S实际上为孔隙半径小于r的孔隙所占的体积比，然而这部分小孔隙却因含有束缚水或是被黏土矿物充填而失去了分形特征^[22]。尽管如此，含水饱和度法计算的分形维数对孔隙结构和储集性较好的储层依然具有参考意义^[23]。

图 11 川东北侏罗系凉高山组储层物性与分形维数关系

Figure 11. Correlation between reservoir physical properties and fractal dimension of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

5.3 分形维数与孔喉参数的关系

孔喉结构参数影响储层的结构、连通性、分选程度等，进而影响储层的分形维数^[24]。D₁与反映孔喉连通特征的排驱压力、中值压力均呈正相关，相关系数R²分别为0.494 5、0.567 0(图 12a、c)，表明随排驱压力、中值压力增大，最大进汞饱和度变小，孔喉结构、连通性变差，储层非均质性增强、分形维数增大。同时，D₁与反映孔隙分布集中程度的分选系数呈负相关，相关系数R²为0.582 4(图 12e)，表明孔隙分布越集中，分选系数越小，分形维数越大。D₂与孔喉参数相关性差(图 12b、d、f)，同样表明大孔隙是储层物性的主要贡献者。

图 12 川东北侏罗系凉高山组孔喉结构参数与分形维数关系

Figure 12. Correlation between pore-throat structure parameters and fractal dimension of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

5.4 分形维数与TOC含量的关系

研究区凉高山组TOC含量与分形维数相关性分析表明，不同孔径的分形维数与TOC含量存在不同的相关性。其中，D_N1和TOC含量相关性较弱(图 13a)，但D_N2与TOC含量呈正相关关系(图 13b)，表明凉高山组页岩TOC含量的高低会对储层孔隙内部结构造成影响，从而进一步影响储层孔隙的分形维数。研究区凉高山组页岩成熟度较高，以Ⅲ型干酪根为主，存在较多气体排出后形成的热成因孔，孔隙形状多不规则，大小一般为纳米级，连通性较差。TOC含量高表明有机质孔增多，且有机质孔以微孔和介孔为主，因此储层非均质性增强、分形维数增大。与此同时，D₁、D₂与TOC含量间均不存在相关性(图 13c-d)，表明TOC含量主要对有机质孔隙(微孔与介孔为主)存在影响，因此与宏孔分形维数的相关性差。

图 13 川东北侏罗系凉高山组TOC含量与分形维数关系

Figure 13. Correlation between TOC content and fractal dimension of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

6. 结论

(1) 研究区凉高山组页岩矿物组成以黏土矿物(平均42.3%)与长英质矿物(平均48.1%)为主导，碳酸盐矿物仅占2.69%。页岩储集空间以黏土矿物层间孔、石英与长石粒间孔以及微裂缝为主。从Ⅰ到Ⅳ类储层，排驱压力增大，最大进汞饱和度减小，储层物性、孔喉参数变差，非均质性增强。

(2) 低温氮气吸附曲线形态符合IUPAC分类中的Ⅳ类曲线，指示狭缝型孔隙。基于“FHH”模型计算的分形维数显示，微孔与介孔分形维数间存在显著线性相关(R²>0.900 0)；孔隙表面分形维数D_N1(平均值2.654 8)大于孔隙结构分形维数D_N2(平均值2.615 6)，表明孔隙表面非均质性更强。

(3) 基于“含水饱和度法”获取的高压压汞分形曲线存在明显的转折点，从Ⅰ类到Ⅳ类储层转折点压力增大，物性变差。D₁均值为2.991 2，D₂均值为2.679 2，表明大孔隙非均质性强于小孔隙，但大孔隙分形维数分布更集中、对储层的贡献更大。

(4) 矿物组分、孔喉结构参数与D₁呈现相关性，与D_N1、D_N2相关性不明显；储层物性方面，孔渗与D₁相关性不大，排驱压力、中值压力与D₁呈正相关，分选系数与D₁呈负相关。TOC含量主要影响微孔与介孔分形维数，与D_N2呈现正相关关系。综上所述，Ⅰ类储层物性最好，勘探前景最大，可以作为下一步有利勘探层段进行优选评价。

致谢: 衷心感谢中国石油勘探开发研究院提供的样品与资料支持！

利益冲突声明/Conflict of Interests

作者陈睿倩是本刊青年编委会成员，未参与本文的同行评审或决策。

Author CHEN Ruiqian is a Young Editorial Board Member of this journal, and she did not take part in peer review or decision making of this article.

作者贡献/Authors’Contributions

李琦、陈睿倩参与论文写作和修改；李琦、李玲、白昕完成实验操作与数据分析；商斐提供样品与资料支持。所有作者均阅读并同意最终稿件的提交。

The manuscript was drafted and revised by LI Qi and CHEN Ruiqian. The experimental operation and data analysis were completed by LI Qi, LI Ling, and BAI Xin, while SHANG Fei provided sample and data support. All authors have read the final version of the paper and consented to its submission.

图 1 四川盆地构造分区及研究区位置

据参考文献[11]修改。

Figure 1. Tectonic divisions of Sichuan Basin and location of study area

下载: 全尺寸图片幻灯片

图 2 川东北侏罗系凉高山组地层柱状图

据参考文献[11]修改。

Figure 2. Stratigraphic column of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 3 川东北侏罗系凉高山组矿物组成条形图

Figure 3. Bar chart of mineral composition in Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 4 川东北侏罗系凉高山组矿物组成三角图

Figure 4. Ternary diagram of mineral composition in Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 5 川东北侏罗系凉高山组页岩储集空间类型

Figure 5. Types of reservoir space in shale from Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 6 川东北侏罗系凉高山组四类储层压汞曲线

Figure 6. High-pressure mercury intrusion curves of four types of reservoirs in Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 7 川东北侏罗系凉高山组页岩低温氮气吸附脱附曲线

Figure 7. Low-temperature nitrogen adsorption and desorption curves of Jurassic Lianggaoshan Formation shale, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 8 川东北侏罗系凉高山组四类储层氮气吸附分形结果

Figure 8. Fractal analysis results of nitrogen adsorption in four types of reservoirs of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 9 川东北侏罗系凉高山组四类储层高压压汞分形结果

Figure 9. Fractal results based on high-pressure mercury intrusion in four types of reservoirs of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 10 川东北侏罗系凉高山组矿物组分含量与分形维数关系

Figure 10. Correlation between mineral composition and fractal dimension of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 11 川东北侏罗系凉高山组储层物性与分形维数关系

Figure 11. Correlation between reservoir physical properties and fractal dimension of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 12 川东北侏罗系凉高山组孔喉结构参数与分形维数关系

Figure 12. Correlation between pore-throat structure parameters and fractal dimension of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

图 13 川东北侏罗系凉高山组TOC含量与分形维数关系

Figure 13. Correlation between TOC content and fractal dimension of Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

下载: 全尺寸图片幻灯片

表 1 川东北侏罗系凉高山组页岩基本物性及孔喉结构参数

Table 1. Basic physical properties and pore-throat structure parameters of shale from Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

样品编号	分类	井深/m	孔隙度/%	渗透率/10^-3 μm²	排驱压力/MPa	中值压力/MPa	最大进汞饱和度/%	分选系数
WB-5	Ⅰ	3 184.66	4.81	0.003 01	0.182 6	1.453 6	90.73	2.798 2
HQ1-2	Ⅰ	1 765.28	2.47	0.060 80	0.292 0	2.135 0	89.89	2.672 6
HQ1-3	Ⅰ	1 889.17	2.30	3.330 00	0.191 6	3.106 2	88.79	2.923 3
WB-4	Ⅰ	3 183.66	4.35	0.086 90	0.182 2	1.566 1	87.75	2.738 1
WB-6	Ⅱ	3 187.00	3.54	0.050 10	0.463 4	3.554 6	91.38	2.752 3
WB-1	Ⅱ	3 176.10	2.10	0.052 00	0.463 9	4.415 2	88.33	2.944 2
HQ1-1	Ⅱ	1 764.20	1.96	0.048 90	0.459 2	3.899 9	91.44	2.623 0
WB-2	Ⅱ	3 177.36	8.08	0.038 00	0.458 9	4.271 6	88.45	2.763 3
WB-3	Ⅱ	3 179.54	8.88	0.033 00	0.457 2	2.056 9	91.73	2.348 5
HQ1-8	Ⅲ	1 902.01	6.60	0.013 20	0.730 0	10.150 4	87.08	2.911 0
HQ1-7	Ⅲ	1 900.39	7.80	0.005 44	0.734 9	22.598 3	87.62	2.874 1
ZX1-1	Ⅳ	1 901.35	11.45	0.014 40	1.135 5	12.969 6	83.02	2.870 0
HQ1-4	Ⅳ	1 894.61	2.15	0.004 69	1.158 8	13.177 9	85.87	2.718 2
HQ1-5	Ⅳ	1 897.15	3.97	0.006 29	1.158 9	4.298 5	87.05	2.494 2
HQ1-9	Ⅳ	1 915.01	5.79	0.006 65	1.160 1	47.529 0	84.67	2.984 1
HQ1-6	Ⅳ	1 898.00	10.19	0.001 53	1.838 1	11.855 9	83.71	2.622 5

下载: 导出CSV

表 2 川东北侏罗系凉高山组页岩氮气吸附分形维数计算结果

Table 2. Fractal dimension calculation of nitrogen adsorption in shale from Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

样品编号	分类	D_N1	R₁²	D_N2	R₂²
WB-5	Ⅰ	2.640 7	0.997 1	2.599 9	0.999 7
HQ1-2	Ⅰ	2.649 3	0.997 9	2.601 9	0.999 0
HQ1-3	Ⅰ	2.652 7	0.986 7	2.584 0	0.996 0
WB-4	Ⅰ	2.485 8	0.991 3	2.578 1	0.992 0
WB-6	Ⅱ	2.607 6	0.994 6	2.599 9	0.952 4
WB-1	Ⅱ	2.665 6	0.994 4	2.650 3	0.996 6
HQ1-1	Ⅱ	2.667 9	0.994 7	2.658 3	0.988 0
WB-2	Ⅱ	2.610 1	0.994 8	2.789 5	0.961 5
WB-3	Ⅱ	2.694 1	0.993 2	2.599 8	0.999 3
HQ1-7	Ⅲ	2.712 6	0.994 6	2.589 3	0.998 5
HQ1-8	Ⅲ	2.684 5	0.990 7	2.603 9	0.999 4
HQ1-6	Ⅳ	2.703 2	0.990 8	2.613 6	0.999 3
HQ1-5	Ⅳ	2.693 0	0.992 2	2.614 4	0.999 0
HQ1-4	Ⅳ	2.665 4	0.998 5	2.620 6	0.999 5
ZX1-1	Ⅳ	2.682 4	0.992 4	2.552 9	0.993 8
HQ1-9	Ⅳ	2.662 1	0.996 4	2.592 8	0.999 4

下载: 导出CSV

表 3 川东北侏罗系凉高山组页岩高压压汞分形维数计算结果

Table 3. Fractal dimension calculation based on high-pressure mercury intrusion in shale from Jurassic Lianggaoshan Formation, northeastern Sichuan Basin

样品编号	分类	D₁	R₁²	D₂	R₂²
WB-5	Ⅰ	2.995 6	0.327 9	2.694 3	0.972 6
HQ1-2	Ⅰ	2.971 0	0.292 5	2.660 8	0.971 8
HQ1-3	Ⅰ	2.976 8	0.454 9	2.677 0	0.966 5
WB-4	Ⅰ	2.983 7	0.380 8	2.719 1	0.992 0
WB-6	Ⅱ	2.996 7	0.407 4	2.685 8	0.981 8
WB-1	Ⅱ	2.984 2	0.265 5	2.700 3	0.918 3
HQ1-1	Ⅱ	2.996 9	0.392 4	2.619 8	0.972 9
WB-2	Ⅱ	2.987 9	0.270 7	2.675 3	0.971 4
WB-3	Ⅱ	2.993 2	0.298 2	2.617 3	0.992 2
HQ1-7	Ⅲ	2.989 9	0.642 9	2.722 1	0.976 3
HQ1-8	Ⅲ	2.999 8	0.230 4	2.746 1	0.965 2
HQ1-6	Ⅳ	2.996 3	0.188 7	2.631 2	0.983 7
HQ1-5	Ⅳ	2.998 7	0.284 2	2.665 9	0.985 9
HQ1-4	Ⅳ	2.999 8	0.231 9	2.664 9	0.955 1
ZX1-1	Ⅳ	2.988 0	0.279 9	2.687 5	0.926 5
HQ1-9	Ⅳ	2.998 3	0.351 2	2.699 1	0.826 3

下载: 导出CSV

参考文献(24)

[1]	王世谦, 陈更生, 黄先平. 四川盆地油气资源潜力及重点勘探领域[J]. 石油学报, 2005, 26(S1): 97-101. WANG Shiqian, CHEN Gengsheng, HUANG Xianping. Petroleum resource potential and major exploration areas in Sichuan Basin[J]. Acta Petrolei Sinica, 2005, 26(S1): 97-101.
[2]	王世谦, 胡素云, 董大忠. 川东侏罗系: 四川盆地亟待重视的一个致密油气新领域[J]. 天然气工业, 2012, 32(12): 22-29. WANG Shiqian, HU Suyun, DONG Dazhong. Jurassic tight oil & gas resources in East Sichuan Basin: a new exploration target[J]. Natural Gas Industry, 2012, 32(12): 22-29.
[3]	何文渊, 何海清, 王玉华, 等. 川东北地区平安1井侏罗系凉高山组页岩油重大突破及意义[J]. 中国石油勘探, 2022, 27(1): 40-49. HE Wenyuan, HE Haiqing, WANG Yuhua, et al. Major breakthrough and significance of shale oil of the Jurassic Lianggaoshan Formation in well Ping'an 1 in northeastern Sichuan Basin[J]. China Petroleum Exploration, 2022, 27(1): 40-49.
[4]	胡东风, 李真祥, 魏志红, 等. 四川盆地北部地区巴中1HF井侏罗系河道砂岩油气勘探突破及意义[J]. 天然气工业, 2023, 43(3): 1-11. HU Dongfeng, LI Zhenxiang, WEI Zhihong, et al. Breakthrough in oil and gas exploration of Jurassic channel sandstone in well Bazhong 1HF in northern Sichuan Basin and its significance[J]. Natural Gas Industry, 2023, 43(3): 1-11.
[5]	邹娟, 乐园, 金涛, 等. 川中地区侏罗系致密砂岩储集层成因及"甜点"预测[J]. 新疆石油地质, 2018, 39(5): 555-560. ZOU Juan, LE Yuan, JIN Tao, et al. Genesis and sweet spot prediction of Jurassic tight sandstone reservoir in central Sichuan Basin[J]. Xinjiang Petroleum Geology, 2018, 39(5): 555-560.
[6]	MANDELBROT B B. Fractals: form, chance and dimension[M]. New York: Freeman W H and Co., 1977.
[7]	李小明, 王亚蓉, 吝文, 等. 湖北荆门探区五峰组—龙马溪组深层页岩微观孔隙结构与分形特征[J]. 天然气地球科学, 2022, 33(4): 629-641. LI Xiaoming, WANG Yarong, LIN Wen, et al. Micro-pore structure and fractal characteristics of deep shale from Wufeng Formation to Longmaxi Formation in Jingmen exploration area, Hubei Province[J]. Natural Gas Geoscience, 2022, 33(4): 629-641.
[8]	王斌, 过敏, 施泽进, 等. 川东南地区志留系石牛栏组海相页岩储层孔隙结构特征及油气地质意义[J/OL]. 成都理工大学学报(自然科学版), (2024-08-21)[2024-10-27]. http://kns.cnki.net/kcms/detail/51.1634.N.20240821.1330.002.html. WANG Bin, GUO Min, SHI Zejin, et al. Characteristics of shale pore structure and its oil and gas geological significance of the Shiniulan Formation, southeast Sichuan Basin[J/OL]. Journal of Chengdu University of Technology (Science & Technology Edition), (2024-08-21)[2024-10-27]. http://kns.cnki.net/kcms/detail/51.1634.N.20240821.1330.002.html.
[9]	杨鑫翊, 王民, 白雪峰, 等. 川东北地区侏罗系凉高山组页岩储集空间特征及页岩油可动性探索[J]. 石油科学通报, 2024, 9(2): 196-212. YANG Xinyi, WANG Min, BAI Xuefeng, et al. Reservoir space characteristics and exploration of shale oil mobility of the Jurassic Lianggaoshan Formation shale in the northeastern Sichuan Basin[J]. Petroleum Science Bulletin, 2024, 9(2): 196-212.
[10]	王道军, 陈超, 刘珠江, 等. 四川盆地复兴地区侏罗系纹层型页岩油气富集主控因素[J]. 石油实验地质, 2024, 46(2): 319-332. doi: 10.11781/sysydz202402319 WANG Daojun, CHEN Chao, LlU Zhujiang, et al. Main controlling factors for oil and gas enrichment in Jurassic laminated shale in Fuxing area of Sichuan Basin[J]. Petroleum Geology & Experiment, 2024, 46(2): 319-332. doi: 10.11781/sysydz202402319
[11]	李楠, 洪海涛, 李国辉, 等. 四川盆地高陡构造区凉高山组页岩油气地质特征[J]. 天然气勘探与开发, 2022, 45(4): 86-95. LI Nan, HONG Haitao, LI Guohui, et al. Geological characteristics of shale oil and gas of Lianggaoshan Formation in high-steep structural zone, Sichuan Basin[J]. Natural Gas Exploration and Development, 2022, 45(4): 86-95.
[12]	郭彤楼. 多旋回盆地叠合—复合控藏在常规非常规天然气勘探中的实践[J]. 地学前缘, 2022, 29(6): 109-119. GU0 Tonglou. Superimposition and composition control on hydrocarbon reservoirs in multi-cycle basins: concept and practice in conventional and unconventional natural gas explorations[J]. Earth Science Frontiers, 2022, 29(6): 109-119.
[13]	成大伟, 张志杰, 洪海涛, 等. 四川盆地东部侏罗系凉高山组沉积、物源及盆山关系[J]. 石油勘探与开发, 2023, 50(2): 1-11. CHENG Dawei, ZHANG Zhijie, HONG Haitao, et al. Sedimentary and provenance characteristics and the basin-mountain relationship of the Jurassic Lianggaoshan Formation in eastern Sichuan Basin, SW China[J]. Petroleum Exploration and Development, 2023, 50(2): 1-11.
[14]	王羽, 金婵, 汪丽华, 等. 应用氩离子抛光—扫描电镜方法研究四川九老洞组页岩微观孔隙特征[J]. 岩矿测试, 2015, 34(3): 278-285. WANG Yu, JIN Chan, WANG Lihua, et al. Characterization of pore structures of Jiulaodong Formation shale in the Sichuan Basin by SEM with Ar-ion milling[J]. Rock And Mineral Analysis, 2015, 34(3): 278-285.
[15]	赵毅楠. 川南地区龙马溪组海相富有机质页岩储层岩石学特征及评价[D]. 成都: 成都理工大学, 2017. ZHAO Yinan. Petrological characteristics and evaluation of marine organic shale reservoirs in Longmaxi Formation, southern Sichuan[D]. Chengdu: Chengdu University of Technology, 2017.
[16]	SCHIEBER J. Common themes in the formation and preservation of intrinsic porosity in shales and mudstones: illustrated with examples across the Phanerozoic[C]//SPE Unconventional Gas Conference. Pittsburgh: SPE, 2010.
[17]	THOMMES M, KANEKO K, NEIMARK A V, et al. Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution(IUPAC Technical Report)[J]. Pure and Applied Chemistry, 2015, 87(9/10): 1051-1069.
[18]	PFEIFER P, AVNIR D. Chemistry in noninteger dimensions between two and three.I.Fractal theory of heterogeneous surfaces[J]. The Journal of Chemical Physics, 1983, 79(7): 3358-3565.
[19]	贺承祖, 华明琪. 储层孔隙结构的分形几何描述[J]. 石油与天然气地质, 1998, 19(1): 17-25. HE Chengzu, HUA Mingqi. Fractal geometry description of reservoir pore structure[J]. Oil & Gas Geology, 1998, 19(1): 17-25.
[20]	彭军, 韩浩东, 夏青松, 等. 深埋藏致密砂岩储层微观孔隙结构的分形表征及成因机理: 以塔里木盆地顺托果勒地区柯坪塔格组为例[J]. 石油学报, 2018, 39(7): 775-791. PENG Jun, HAN Haodong, XIA Qingsong, et al. Fractal characterization and genetic mechanism of micro-pore structure in deeply buried tight sandstone reservoirs: a case study of Kalpintag Formation in Shuntuoguole area, Tarim Basin[J]. Acta Petrolei Sinica, 2018, 39(7): 775-791.
[21]	赵会涛, 郭英海, 杜小伟, 等. 鄂尔多斯盆地高桥地区本溪组砂岩储层微观孔隙多重分形特征[J]. 地质科技通报, 2020, 39(6): 175-184. ZHAO Huitao, GUO Yinghai, DU Xiaowei, et al. Micro-pore multifractal characteristics of Benxi Formation sandstone reservoir in Gaoqiao area, Ordos Basin[J]. Bulletin of Geological Science and Technology, 2020, 39(6): 175-184.
[22]	DANIEL J K, ROSS R, BUSTIN M. Shale gas potential of the Lower Jurassic Gordondale Member, northeastern British Columbia, Canada[J]. Bulletin of Canadian Petroleum Geology, 2007, 55(1): 51-75.
[23]	赖锦, 王贵文, 郑懿琼, 等. 低渗透碎屑岩储层孔隙结构分形维数计算方法: 以川中地区须家河组储层41块岩样为例[J]. 东北石油大学学报, 2013, 37(1): 1-7. LAI Jin, WANG Guiwen, ZHENG Yiqiong, et al. Method For calculating the fractal dimension of the pore structure of low permeability reservoirs: a case study on the Xujiahe Formation reservoir in central Sichuan Basin[J]. Journal of Northeast Petroleum University, 2013, 37(1): 1-7.
[24]	张婷, 徐守余, 杨珂. 储层微观孔隙结构分形维数应用[J]. 大庆石油学院学报, 2010, 34(3): 44-47. ZHANG Ting, XU Shouyu, YANG Ke. Application of fractal dimension of micro-pore structure[J]. Journal of Daqing Petroleum Institute, 2010, 34(3): 44-47.