Fractal characteristics of pores in marine-continental transitional shale of Shanxi Formation in Kaiping syncline, North China Basin

  • Peng LI , 1, 2 ,
  • Haoran CHEN 3, 4 ,
  • Yang WANG , 3, 4 ,
  • Yanming ZHU 3, 4 ,
  • Zhenfei JIANG 3, 4 ,
  • Yi WANG 3, 4
Expand
  • 1. Hebei Province Coal Geological Exploration Institute,Xingtai 054000,China
  • 2. Hebei Province Unconventional Natural Gas Exploration and Research Center,Xingtai 054000,China
  • 3. Key Laboratory of Coalbed Methane Resources and Reservoir Formation Process,Ministry of Education,Xuzhou 221008,China
  • 4. School of Resource and Geosciences,China University of Mining and Technology,Xuzhou 221116,China

Received date: 2024-12-03

  Revised date: 2025-01-23

  Online published: 2025-02-24

Supported by

The“Three-Three-Three Talents Project” Funded by Hebei Province, China(C20221057)

the Geological Project of Public Finance Special Fund of Hebei Coal Geological Bureau(13000021MY3U4MEAG7BO9)

the National Natural Science Foundation of China(42172156)

Abstract

The study of pore structure and heterogeneity in shales is a key focus in unconventional reservoir research. Transitional shales are widely distributed in China and exhibit significant exploration potential. However, research on the pore structures of transitional shales significantly lags behind that of marine shales, particularly in the northern part of North China Basin, where resource potential is substantial but related studies are limited. This paper investigates the pore structures and heterogeneity of transitional facies shale samples from the Shanxi Formation, using core samples from Well KP1 in the Kaiping Syncline area of the northern North China Basin, employing various experimental methods. The results reveal that the shale is primarily composed of quartz and clay minerals. The total organic carbon (TOC) content ranges from 0.3% to 3.43%, with vitrinite reflectance(R O)values exceeding 1.6%, indicating a high level of thermal maturity. Pores in the shale are predominantly intergranular clay mineral pores, mainly in the form of flat, slit-shaped, with most pore sizes being smaller than 100 nm. Fractal dimensions for ultra-large pores, macropores, mesopores, and micropores range of 2.29-2.34, 2.25-2.69, 2.69-2.76, and 2.36-2.56, respectively, while the overall fractal dimension (D s) ranges from 2.60 to 2.63. The high clay content and low abundance of Type I kerogen in the study area result in poorly developed organic pores (micropores), which is the main reason for the absence of significant correlations between TOCR O, and fractal dimensions. The quartz-associated intergranular pores exhibit a highly complex pore structure, characterized by distinctly irregular pore and fracture shapes. This irregularity significantly increases the complexity of the overall pore system, influencing the storage and permeability characteristics of the reservoir. Consequently, higher clay mineral content contributes to a more complex overall pore structure. The fractal dimensions show a positive correlation with pore surface area and porosity, indicating that greater pore structure complexity enhances the internal pore space of the shale, improving its capacity to store both adsorbed and free gas. The impact of fractal dimensions on permeability varies with pore size: complex mesoporous systems increase flow resistance, hindering gas desorption and diffusion, while complex micropore structures enhance connectivity between larger pores, improving the flow capacity of the shale reservoir.

Cite this article

Peng LI , Haoran CHEN , Yang WANG , Yanming ZHU , Zhenfei JIANG , Yi WANG . Fractal characteristics of pores in marine-continental transitional shale of Shanxi Formation in Kaiping syncline, North China Basin[J]. Natural Gas Geoscience, 2025 , 36(7) : 1330 -1344 . DOI: 10.11764/j.issn.1672-1926.2025.01.009

0 引言

页岩气作为一种具有巨大潜力的非常规天然气资源,可以有效缓解我国油气资源短缺的压力。我国页岩气资源丰富,累计可采储量可达13.23×1012 m3,其中海陆过渡相页岩气可采储量占比可达20%以上,勘探潜力巨大1-4。页岩储层是页岩气生储的重要空间,页岩储层孔隙结构影响着页岩气的赋存与产出4-6。近年来的勘探研究表明,正确认识页岩孔隙结构,对页岩储层的精细化表征、探究页岩气勘探开发具有重要意义。
众多学者对于页岩孔隙结构表征方法已做了大量研究,目前常用扫描电镜(SEM)、场发射扫描电镜(FE-SEM)、原子力显微镜(AFM)等成像法表征孔隙的类型、形态与大小,常用高压压汞、低温液氮吸附等注入法与吸附法定量表征孔隙的体积、比表面积等7-16,单一的表征手段不能准确反映页岩孔隙结构特征,因此结合多种方法联合表征已成为目前研究孔隙结构的主要方法11-16。然而页岩孔隙结构的复杂程度与非均质性特征难以通过实验测试直接表征,BENOIT17提出分形理论,用于研究事物的自相似性,众多学者将分形理论引入页岩储层研究,基于流体注入法、低压气体吸附法的分形模型在页岩孔隙非均质性表征中得到了广泛应用18-23
前人24-26对沁水盆地山西组—太原组海陆过渡相页岩孔隙的分形特征进行研究,指出海陆过渡相页岩主要发育孔隙形态复杂的黏土矿物孔,认为分形维数与黏土矿物具有良好的正相关性,与有机质含量的关系不明确。魏建光等27对宁武盆地山西组海陆过渡相页岩做了部分研究工作,指出高成熟度、高有机质含量以及高黏土矿物有利于微孔形成,从而使得孔隙复杂程度升高。部分学者认为黔东北地区龙潭组海陆过渡相页岩分形维数与TOC呈弱正相关关系,与总孔体积和比表面积均呈良好的正相关关系28-29。谢卫东等30对鄂尔多斯盆地延长组、延安组海陆过渡相页岩做了大量研究,认为孔隙中吸收大量水分使得黏土矿物孔隙结构均质性升高是分形维数降低的主要原因。前人31-32对海陆过渡相页岩的研究主要集中在鄂尔多斯盆地、沁水盆地、四川盆地以及南华北盆地等,华北盆地北部石炭系—二叠系海陆过渡相页岩分布广泛且厚度较大,具有优秀的页岩气勘探开发潜力,但目前针对华北盆地北部页岩孔隙结构的相关研究较少。亟需对华北盆地北部页岩孔隙非均质性特征进行探究。
为此,本文以华北盆地北部开平向斜下二叠统山西组海陆过渡相页岩为研究对象,利用高压压汞、低温液氮吸附、扫描电镜、XRD全岩分析等方法,对页岩储层的孔隙结构、矿物组成等进行研究,并基于多种分形模型对不同尺寸孔隙分形特征进行定量表征,讨论了分形维数与矿物组成、TOC、镜质体反射率(R O)以及孔隙结构的关系,以期为华北盆地北部页岩储层特征研究提供科学依据。

1 地质背景

华北盆地位于我国华北地区,具有复杂的构造演化史,盆地如今构造格局主要是经过印支运动、燕山运动与喜马拉雅运动改造后形成33-34。本次研究区开平向斜位于华北盆地东北部唐山坳陷区,研究区内构造以开平向斜为主构造,西北翼地层急倾斜局部倒转,东南翼地层平缓而多褶曲,自北向南发育有碾子庄西背斜、高各庄北向斜、小齐各庄西背斜、杨家泊背斜。断裂构造发育走向以北东向为主,其次为北西向35图1)。石炭世—二叠世期间,研究区处于海陆交替阶段,是盆地环境由陆表海向陆内湖过渡的关键时期,山西组沉积环境为三角洲相沉积,沉积地层主要以黑色页岩与煤层为主36。根据研究区KP1钻井信息揭露(图1),山西组主要以砂岩、泥页岩以及薄层煤层为主,厚度为69.50 m。
图1 开平向斜构造位置及KP1井地层综合柱状图[图(a)据文献[37]修编);图(b),图(c)据文献[38]修编]

Fig.1 Location of Kaiping Syncline structure and strata comprehensive histogram in Well KP1[(a) Revised according to Ref.[37]; (b) (c) Revised from Ref. [38]]

2 样品采集与研究方法

本文研究采样点位于华北盆地北部开平向斜KP1井(北纬39°45′46′′,东经118°21′25′′),地理位置位于河北省唐山市丰南区,构造位置属于开平向斜东翼。本文研究共收集和分析了8个具有代表性的KP1井山西组页岩样品,采样深度为1 705.20~1 791.69 m。样品的具体信息见表1,主要岩性包括页岩和炭质页岩。为了便于对孔隙分布特征进行研究,本文根据国际理论和应用化学协会(IUPAC)的孔隙分类,并结合研究区页岩孔隙特征,将页岩孔隙划分为微孔(<2 nm)、介孔(2~50 nm)、宏孔(50~1 000 nm)与超大孔(>1 000 nm)。页岩样品实验分析包括总有机碳含量检测、孔—渗测试、XRD全岩矿物定量分析检测、扫描电镜、高压压汞测试以及低温液氮吸附测试等,以上分析测试均在山东省煤田地质规划勘察研究院实验室进行。
表1 研究区页岩样品信息

Table 1 Shale sample information in the study area

样品编号 井号 层位 深度/m
KP1-D-Y3 KP1 山西组 1 705.20
KP1-D-Y9 KP1 山西组 1 717.05
KP1-D-Y12 KP1 山西组 1 720.37
KP1-D-Y13 KP1 山西组 1 722.17
KP1-D-Y14 KP1 山西组 1 723.72
KP1-D-Y17 KP1 山西组 1 727.18
KP1-D-Y18 KP1 山西组 1 728.37
KP1-D-Y37 KP1 山西组 1 791.69
采用CS-230碳硫分析仪对岩石的总有机碳含量进行测定,测试依据国家标准《沉积岩中总有机碳测定》(GB/T 19145—2022),通过固态红外吸收法计算总有机碳含量。采用HKC-2型氦孔隙度测定仪对页岩的孔隙度进行测定,测试依据国家标准《页岩孔隙度、渗透率和饱和度测定》(GB/T 34533—2023),把一定压力的氦气充入到参比室中,然后将参比室的气体扩散到与其相连的样品室内,待压力平衡后,根据波义尔定律,计算页岩孔隙度。采用STY-2岩心渗透率测定仪对页岩渗透率进行测定,测试依据国家标准《页岩孔隙度、渗透率和饱和度测定》(GB/T 34533—2023),将岩样置于页岩岩心夹持器中,施加围压密封。在页岩岩样一端注入一定压力的气体,待压力稳定后,获取页岩岩样另一端的气体在单位时间内的流量,利用达西定律计算出页岩岩样的渗透率。采用蔡司Supra 55扫描电子显微镜对页岩微观特征进行分析,测试依据中华人民共和国石油天然气行业标准《岩石样品扫描电子显微镜分析方法》(SY/T 5162/2020),对页岩矿物组成、孔隙结构进行观察。采用Rigaku SmartLab9 X射线衍射仪对页岩矿物组成进行定量分析,整个实验过程依照中华人民共和国石油天然气行业标准《沉积岩中黏土矿物和常见非黏土矿物X射线衍射分析方法》(SY/T5163—2018)进行。页岩高压压汞测试采用Micromeritics Autopore9510孔隙仪,实验前将样品处理为边长1 cm×1 cm×1 cm立方体,实验设置进汞压力范围为0.001~140 MPa,通过实验获取样品进汞曲线,并基于进汞曲线计算孔径分布特征,具体实验操作依照中华人民共和国能源行业标准《页岩全孔径分布的测定压汞—吸附联合法》(NB/T 14008—2021)。低温液氮吸附实验采用ASAP2020比表面积与孔径分析仪,实验脱气条件250 ℃,P/P 0值介于0.01~0.995之间,在临界温度下测量氮气的吸附量和脱附量,基于BJH模型计算孔径分布特征,具体实验操作方法依照中华人民共和国石油天然气行业标准《岩石比表面积和孔径分布测定静态吸附容量法》(SY/T 6154—2019)。

3 实验结果

3.1 矿物成分和有机地球化学特征

XRD全岩矿物定量分析结果如图2所示。研究区页岩样品主要由石英与黏土矿物组成,长石含量较低,部分含白云石、菱铁矿与黄铁矿。脆性矿物含量(石英+长石+黄铁矿+白云石)介于36.40%~53.00%之间,平均值为44.79%,黏土矿物含量介于35.10%~51.60%之间,平均值为43.76%。前人研究表明,脆性矿物含量高的页岩储层有利于天然裂缝的发育以及后期压裂工程17,研究区页岩样品脆性矿物含量均大于40.00%,达到有利层段优选标准,但黏土矿物的发育会使压裂难度增大5
图2 山西组页岩样品矿物组成

Fig.2 Mineral content distribution of shale samples from Shanxi Formation

有机质含量与成熟度是评价烃源岩生烃能力的重要指标,研究区页岩TOC含量介于0.3%~3.43%之间,平均值为1.59%,根据煤系泥页岩有机质丰度评价标准可知,样品属于一般烃源岩、好烃源岩。样品R O值介于1.8%~2.23%之间,平均值为1.97%,均已达到高成熟阶段。页岩干酪根组分以腐泥组为主,不含壳质组,根据干酪根判别图(图3)显示,均属于混合型(Ⅱ型)干酪根。
图3 页岩干酪根类型判别图

Fig.3 Discriminant diagram of shale kerogen types

3.2 页岩孔隙综合表征

3.2.1 扫描电镜下页岩孔隙特征

扫描电镜结果显示,页岩样品主要发育无机孔,其中黏土矿物粒间孔最为发育,有机质孔发育较少,镜下少见。粒间孔孔径范围从几十纳米到几微米不等,形状不规则,分布不均,是页岩气储存的主要空间[图4(a),图4(b),图4(g)]。此外石英与黏土矿物以及有机质之间的粒间孔缝较为发育,主要呈尺寸较大且形状不规则的特点,细长状、椭圆状、锯齿状孔缝均有发育[图4(c),图4(h),图4(i)],在一些矿物内部(闪锌矿)发育晶内孔缝,其特点为发育规模较少,分布不均[图4(d)]。层间裂隙沿层理表现为弯曲的线状,长度一般几十到几百微米不等,宽度为几微米[图4(e)]。溶蚀孔多表现为孤立孔,主要呈雨痕状溶蚀,一般认为是溶蚀作用较弱或溶蚀时间较短的环境形成39,此类孔隙具有形状规则的特点,整体呈圆形或椭圆形,大小从0.2~2 μm不等[图4(f),图4(h),图4(i)]。部分样品发育有草莓状黄铁矿[图4(c)]。
图4 山西组页岩样品扫描电镜图像

(a),(b)大量发育的黏土矿物粒间孔隙;(c)草莓状黄铁矿与有机质共生;(d)闪锌矿晶内孔缝;(e)大尺寸的层间裂隙;(f)椭圆状、圆状溶蚀孔隙;(g)黏土矿物与有机质共生,粘土矿物内发育粒间孔;(h)石英与有机质共生,形成粒间孔缝;(i)石英与黏土矿物共生,形成粒间孔缝

Fig.4 SEM images of shale samples from Shanxi Formation

3.2.2 高压压汞实验

页岩样品的高压压汞曲线形态相似,整体呈3段形态[图5(a)]。第一阶段初始上升段(0 MPa<P<3 MPa),随着压力增大,样品的汞饱和度增长缓慢,表明样品宏孔体积较小,第二阶段中间平稳段(3 MPa<P<90 MPa),随着压力增大,样品进汞饱和度稳定增长,表明此阶段对应的孔隙体积较为发育,第三阶段末端上翘段(90 MPa<P<140 MPa),样品的进汞饱和度增加速率相对降低,表明此阶段对应孔隙体积较小。页岩样品的孔体积分布曲线整体呈双峰状,主要发育孔径小于100 nm孔隙,其中小于10 nm的孔隙最为发育,大于100 nm的孔隙发育较少[图5(b)]。
图5 基于高压压汞法毛管压力曲线(a)及孔径分布曲线(b)

Fig.5 Capillary pressure curve(a) and pore size distribution curve(b) based on high pressure mercury intrusion method

3.2.3 低温液氮吸附实验

依照IUPAC等温线分类标准(图6),页岩样品的低温液氮吸附—脱附曲线与Ⅳ型相近,Ⅳ型等温线一般认为由介孔类材料产生。由于毛细管凝聚作用,脱附曲线往往不与吸附曲线重合,产生滞后环,滞后环类型被认为可以反映孔隙结构特征40-41,样品的滞后环类型近似为H3型和H4型,H3型滞后环一般由片状颗粒材料,如黏土、裂隙孔产生,在较高相对压力区域没有表现出饱和吸附,反映的孔主要包括平板狭缝结构、裂缝等。H4型滞后环是I型和II型吸附等温线的复合,H4型多出现在微孔和介孔介质中。
图6 IUPAC等温线及滞后环分类标准

Fig.6 Classification criteria for IUPAC isotherms and hysteresis loops

页岩的液氮吸附—脱附曲线如图7所示,根据曲线形态可分为3个阶段:第一阶段(0<P/P 0<0.45),此阶段页岩样品吸附量增长缓慢,呈轻微上凸状,对应液氮在孔隙表面发生单分子层吸附;第二阶段(0.45<P/P 0<0.9),样品吸附速度加快,说明吸附进入多分子层吸附阶段;第三阶段(0.9<P/P 0<0.995),随着相对压力增加,样品吸附量迅速增加,并且随着压力增加,并未出现明显饱和吸附现象。
图7 山西组页岩样品低温液氮吸附—脱附曲线

(a) H3型回滞环样品;(b) H4型回滞环样品

Fig.7 Adsorption-desorption isotherms of shale samples from Shanxi Formation

大量研究表明,根据BJH模型计算孔体积分布特征时,利用液氮吸附脱附曲线的脱附支数据计算得到的样品孔体积分布曲线会得到假峰,峰值在4 nm左右。出现假峰现象的原因是液氮的脱附过程具有滞后性,当压力开始降低时,较大孔隙的液氮会立刻脱附,而较小的孔隙只有当压力降低到某一阈值时才会发生液氮脱附,当开始脱附之后,脱附曲线突然下降,并与吸附曲线重合,这一现象出现时对应的孔径范围在4 nm左右2542,故选取吸附支数据并基于BJH模型计算得到页岩样品的孔体积分布曲线,页岩样品的孔体积分布曲线基本一致,主要发育小于10 nm的孔隙,其中1~2 nm的孔隙最为发育(图8)。
图8 基于液氮吸附法孔径分布曲线

(a) H3型回滞环样品;(b) H4型回滞环样品

Fig.8 Pore size distribution curve based on nitrogen adsorption method

4 讨论

4.1 孔隙分形维数

分形理论可用于定量描述孔隙结构复杂程度,分形特征可以用分形维数定量表示,一般认为页岩孔隙结构的分形维数介于2~3之间,分形维数越接近2,页岩孔隙结构的非均质性越弱,孔隙结构越简单,分形维数越接近3,页岩孔隙结构的非均质性越强,孔隙结构越复杂43-44。并且随着研究的深入,众多学者发现分形维数可以较好地反映页岩孔隙结构特征,尤其是在评价页岩气储集空间方面有着重要意义45-48。较高的分形维数往往反映页岩孔隙的复杂连接结构、微裂缝的发育等,而较低分形维数则可能表明孔隙结构相对简单,流体流动阻力小等。此外,分形维数的变化还能在一定程度上反映区域沉积环境及成岩作用等的变化特征。基于此,本文研究选取不同的分形模型并结合高压压汞、低温液氮吸附实验数据,计算不同尺寸孔隙的分形维数,表征研究区页岩孔隙结构特征。

4.1.1 基于高压压汞实验的孔隙分形特征

利用高压压汞数据,结合分形模型,计算大于50 nm的孔隙的分形维数,计算公式如下48
S H g = P c - ( 2 - D )
式(1)两边同时取对数可得:
L g S H g = ( D - 2 ) P c
式中: S H g为进汞饱和度,%;P c为气水毛管压力,MPa,D为分形维数,无量纲。
根据式(2)可知,如果页岩储层的孔隙结构具有分形特征,则在双对数坐标中,毛管压力与汞饱和度成直线关系,分形维数D=k+2,k为双对数坐标下拟合直线的斜率。
由页岩样品LgP c与LgS分形特征图可知(图9),样品的拟合曲线呈现明显的两段性,可划分为宏孔段与超大孔段(本文将50~1 000 nm的孔隙定义为宏孔,1 000~10 000 nm的孔隙定义为超大孔),两段曲线的拟合程度均较高,其中超大孔段决定系数R 2介于0.983~0.998之间,宏孔段决定系数R 2介于0.877~0.999之间,表明研究区页岩样品的超大孔与宏孔具有明显的分形特征。超大孔段分形维数D 1值介于2.29~2.34之间,平均值为2.31,宏孔段分形维数D 2值介于2.25~2.69之间,平均值为2.52。
图9 页岩样品LgP c与LgS分形特征

Fig.9 Fractal characteristics of LgP c and LgS of shale samples

4.1.2 基于液氮吸附实验的孔隙分形特征

目前FHH(Frenkel Halsey Hill)模型分形维数计算方法被广泛应用于处理液氮吸附数据,该模型计算公式如下49
L n V = D - 3 L n   L n P 0 P + C
式中:V表示的是在压力为P时的吸附体积,cm3/g;P为实验平衡时的氮气分压,MPa;P 0为是液氮温度下氮气的饱和蒸汽压,MPa;D是分形维数,无量纲;C为常数。
由页岩样品LnV和Ln(Ln(P 0 /P))分形特征图可知(图10),样品的拟合曲线呈现明显的两段性,可划分为介孔段与微孔段。两段曲线的拟合程度均较高,其中介孔段决定系数R 2值介于0.938~0.977之间,微孔段决定系数R 2值介于0.987~0.998,表明研究区页岩样品的介孔和微孔具有明显的分形特征。介孔段分形维数D 3值介于2.69~2.76之间,平均值为2.70,微孔段分形维数D 4值介于2.36~2.56之间,平均值为2.50。
图10 页岩样品LnV和Ln(Ln(P/P 0))分形特征

Fig.10 Fractal characteristics of LnV and Ln(Ln(P/P 0)) of shale samples

4.1.3 综合分形维数计算

图11(a)可知,页岩样品的孔隙体积主要由介孔(46.48%~65.46%)与微孔(20.06%~38.91%)贡献,超大孔(2.30%~3.74%)与宏孔(4.38%~21.37%)对孔隙体积的贡献极少。由图11(b)可知,样品分形维数整体呈D 3>D 2>D 4>D 1的趋势,说明本文研究区页岩样品中介孔孔隙结构较为复杂,原因在于样品的介孔数目较多,且主要为黏土矿物粒间孔,具有形状、尺寸非均质性强的特点,从而导致孔隙结构复杂程度增大。
图11 不同孔径孔隙体积占比(a)及分形维数(b)箱线图

(a) 体积占比;(b) 分形维数

Fig.11 Box plot of the proportion of pore volume with different pore sizes (a) and fractal dimension (b)

页岩超大孔的分形维数最低,表明超大孔(>1 000 nm)孔隙结构较为简单,这与前人研究具有一定差异,大部分学者认为页岩孔隙中较大孔隙结构具有多尺度特征,往往更加复杂,从而导致分形维数较大1850,而研究区超大孔分形维数低的原因可能为超大孔的数目极少,仅占总孔隙体积的2.3%~3.6%,较少的超大孔之间孔隙结构差异性较低,并且根据镜下特征可知,超大孔主要包括黏土矿物孔隙、溶蚀孔、微裂隙等,除黏土矿物孔隙以及石英粒间孔结构较为复杂外,溶蚀孔与微裂隙等孔隙结构均较为简单,形状规则,从而表现为超大孔分形维数最低。
为了更加准确地表征孔隙结构的复杂程度,依据不同孔径段的体积占比作为权重,计算样品的综合分形维数。权重计算公式如下:
D s = i = 1 n V i V · D i
式中:D s为综合分形维数,无量纲;Vi 为不同孔径段的体积,cm3/g;V为样品孔隙总体积,cm3/g;Di 为不同孔径段的分形维数,无量纲。
由计算结果(表2)可知,研究区页岩样品综合分形维数D s值介于2.60~2.63之间,平均值为2.62。
表2 页岩样品分形维数统计

Table 2 Fractal dimension statistics of shale samples

样品编号 超大孔 宏孔 介孔 微孔 D s
D 1 R 2 D 2 R 2 D 3 R 2 D 4 R 2
KP1-D-Y3 2.297 0.998 2.686 0.999 2.760 0.952 2.359 0.987 2.609
KP1-D-Y9 2.344 0.987 2.521 0.917 2.717 0.955 2.514 0.995 2.615
KP1-D-Y12 2.337 0.992 2.246 0.938 2.711 0.953 2.487 0.989 2.598
KP1-D-Y13 2.335 0.989 2.526 0.877 2.710 0.960 2.525 0.995 2.628
KP1-D-Y14 2.327 0.994 2.330 0.925 2.686 0.97 2.555 0.998 2.617
KP1-D-Y17 2.291 0.998 2.582 0.939 2.689 0.977 2.553 0.997 2.634
KP1-D-Y18 2.301 0.988 2.560 0.934 2.717 0.938 2.474 0.995 2.615
KP1-D-Y37 2.287 0.994 2.689 0.945 2.703 0.958 2.528 0.997 2.627

4.2 分形维数与TOCR O的关系

通过分形维数D s与页岩TOCR O的相关性分析可知,页岩样品TOCR OD s无明显相关关系(图12)。前人研究表明,TOC是影响海相页岩孔隙复杂程度的重要因素,在有机质热演化阶段,随着成熟度的增长,有机质内会形成大量微孔,使得页岩内部孔隙的复杂程度增加184651。虽然本研究区海陆过渡相页岩样品处于过成熟阶段,但TOCR O与页岩分形维数并无正相关关系。前人对华北盆地研究表明,华北盆地页岩孔隙具有明显的有机质孔不发育的特点,原因可能是有机质主要为Ⅱ型、III型干酪根,生烃潜量较低,不易发育有机孔隙,或因黏土矿物含量较高,抑制了有机孔隙的发育52-53,同时也有部分学者通过干酪根提取对有机质孔进行研究,认为对于海陆过渡相页岩,孔隙复杂程度与TOC含量无关,而与有机质类型有关54。根据本文研究测试结果可知,研究区页岩样品有机质孔明显不发育,黏土矿物含量高,惰质组含量较高,干酪根类型以Ⅱ型为主,富生烃潜力的Ⅰ型干酪根较少,这一特点导致了研究区有机孔隙发育较差,较少的有机孔隙对整个页岩孔隙网络的复杂程度影响有限,从而表现出TOCR O与分形维数无明显相关关系。
图12 D sTOCR O的关系

(a) D sTOC相关图解;(b) D sR O相关图解

Fig.12 The relationship between D s and TOCR O

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

分形维数可以定量表征页岩孔隙结构复杂程度,成岩矿物是页岩无机孔隙的主要载体,因此矿物组分会在一定程度上对孔隙复杂程度产生影响。通过建立矿物组分与各分形维数之间的Pearson相关系数矩阵探究各参数间相关关系。基于前人研究,认为相关系数|r|<0.4为不相关或极弱相关,|r|≥0.4为具有相关性关系55。由图13可知,长石含量与分形维数无明显的相关性关系,石英含量与D 1呈正相关关系(r=0.73),黏土矿物含量与D s呈正相关关系(r=0.70)。长石虽然易在溶蚀作用下形成较大的孔隙,但根据扫描电镜及XRD全岩分析可知,研究区页岩长石含量较低,且由长石形成的溶蚀孔数量极少,长石对孔隙复杂程度贡献低,故表现为无明显相关性;研究区与石英相关的孔隙结构主要有粒间孔隙、粒内溶蚀孔等,根据镜下特征,虽然石英溶蚀孔隙结构较为规则,但石英与黏土矿物、有机质复合体的孔隙结构较为复杂,特别是石英与黏土矿物、有机质结合部位,形成的较大尺寸的粒间孔缝形状明显不规则,细长状、椭圆状、锯齿状孔缝均有发育,这一特点增大了孔隙结构的复杂程度,且前人研究表明,石英作为刚性骨架对原生矿物孔隙具有一定的支撑作用,减小了压实作用对页岩孔隙造成的影响,增大了孔隙结构复杂程度56-57,故石英与超大孔分形维数呈一定的正相关关系;综合分形维数与黏土矿物含量呈正相关关系,在地质过程中,黏土矿物由于其塑性高,孔隙容易发生形变,导致孔隙结构不规则58-59。本研究结果与前人研究一致,研究区页岩黏土矿物含量较高,主要发育黏土矿物粒间孔,孔隙形状呈明显的不规则、大小、形态多样的特点,因此黏土矿物含量越高,整体孔隙结构越复杂,页岩分形维数越大。
图13 分形维数与矿物组成、孔隙结构参数的相关关系矩阵

Fig.13 The correlation matrix diagram of fractal dimension with mineral composition and pore structure parameters

4.4 分形维数与孔隙结构的关系

分形维数在一定程度上能够反映页岩孔隙结构特征,通过分形维数与孔隙比表面积、孔隙度、渗透率以及中值半径的Pearson相关系数矩阵图对分形维数与孔隙结构的关系进行研究(图13)。页岩气在页岩孔隙中主要以吸附态为主,而比表面积是控制吸附气含量的主要因素,D s与页岩比表面积呈明显的正相关关系(r=0.72),表明样品的比表面积受整体孔隙复杂程度的影响,整个孔隙空间结构越复杂或孔隙表面的不规则性越大,页岩内部的比表面积也越大。孔隙度可以反映样品整体的孔隙空间,是影响游离气含量的因素之一,D s与孔隙度呈较弱的正相关关系(r=0.47),说明样品的孔隙度受样品整体孔隙的影响,原因在于孔隙度反映样品整体的孔隙空间特征,分形维数越大,页岩孔隙结构越复杂,孔隙内部空间越大。同时孔隙度与比表面积也呈现出一定的正相关关系,表明页岩内部孔隙空间越大,比表面积也会相应增大。渗透率可以反映流体通过孔隙的能力,是储层评价的重要指标,D 3与渗透率呈较好的负相关关系(r=-0.61),说明介孔孔隙发育越复杂,样品的渗流能力越弱,原因在于黏土矿物中介孔体系结构不规则,连通孔隙的孔喉较细,使得渗流阻力加大,不利于气体的解吸与扩散,从而抑制渗流作用,同时渗透率与D 4D s呈一定的正相关关系(r分别为0.49、0.41),前人研究认为微孔孔隙结构复杂会增加孔隙迂曲度,从而抑制渗流作用60,但同时有部分学者认为复杂微孔结构可以较好地连通页岩中的较大孔隙,增强了孔隙连通性,从而提高了页岩储层渗流能力61。中值半径与分形维数相关性关系较差,这表明中值孔径只能表示孔径分布的中值,不能反映整体孔隙的大小、分布等差异,因此中值半径对孔隙复杂程度的影响较小。

5 结论

(1)华北盆地开平向斜山西组海陆过渡相页岩矿物成分主要为黏土矿物与石英,长石、菱铁矿含量较低。页岩TOC值介于0.3%~3.43 %之间,平均值为1.59 %,R O值均大于1.6 %,处于高成熟阶段。结合扫描电镜、高压压汞实验与低温液氮吸附实验得到研究区页岩孔隙以黏土矿物粒间孔最为发育,孔隙形态以平板狭缝孔为主,主要发育小于100 nm的孔隙。
(2)研究区页岩孔隙具有明显的多重分形特征,利用多种分形模型计算不同尺寸孔径段的分形维数,得到超大孔段分形维数D 1值介于2.29~2.34之间,宏孔段分形维数D 2值介于2.25~2.69之间,介孔段分形维数D 3值介于2.69~2.76之间,微孔分形维数D 4值介于2.36~2.56之间,综合分形维数D s值介于2.60~2.63之间。
(3)研究区高黏土含量、Ⅰ型干酪根少的特点导致了研究区有机孔隙发育较差,较少的有机孔隙对整个页岩孔隙网络的复杂程度影响有限,从而表现出TOCR O与分形维数无明显相关关系。
(4)与石英相关的粒间孔孔隙结构较为复杂,形成的较大尺寸的粒间孔缝形状明显不规则,细长状、椭圆状、锯齿状孔缝均有发育,增大了孔隙结构的复杂程度。黏土矿物粒间孔孔隙形状呈明显的不规则、大小、形态多样的特点,因此黏土矿物含量越高,整体孔隙结构越复杂。
(5)分形维数可以较好地反映页岩孔隙结构特征。分形维数与孔隙比表面积、孔隙度均呈一定的正相关关系,表明孔隙结构越复杂,页岩内部孔隙空间越大,页岩储存吸附气与游离气的能力越强。不同尺寸孔隙分形维数对渗透率的影响不同,复杂介孔体系会使得渗流阻力加大,不利于气体的解吸与扩散,复杂微孔结构可以较好地连通页岩中的较大孔隙,提高了页岩储层渗流能力。
[1]
LI P, ZHANG J C, TANG X, et al. Assessment of shale gas potential of the Lower Permian transitional Shanxi-Taiyuan shales in the southern North China Basin[J].Australian Journal of Earth Sciences, 2021, 68(2):262-284.

[2]
GUO W, GAO J L, LI H, et al. The geological and production characteristics of marine-continental transitional shale gas in China:Taking the example of shale gas from Shanxi Formation in Ordos Basin and Longtan Formation in Sichuan Basin[J]. Mineral Exploration, 2023, 14(3):448-458.

[3]
尹开贵, 张朝, 单长安, 等. 昭通国家级页岩气示范区海陆过渡相页岩气的地质特征与资源潜力[J]. 天然气工业, 2021, 41(S1):30-35.

YIN K G, ZHANG C, SHAN C A, et al. Geological characteristics and resource potential of the transitional shale gas in the Zhaotong National Shale Gas Demonstration Area[J]. Natural Gas Industry, 2021, 41(S1):30-35.

[4]
王建, 郭秋麟, 赵晨蕾, 等. 中国主要盆地页岩油气资源潜力及发展前景[J].石油学报, 2023, 44(12):2033-2044.

WANG J, GUO Q L, ZHAO C L, et al. Potentials and prospects of shale oil gas resources in major basins of China[J]. Acta Petrolei Sinica, 2023, 44(12):2033-2044.

[5]
邹才能, 董大忠, 王社教, 等. 中国页岩气形成机理、地质特征及资源潜力[J]. 石油勘探与开发, 2010, 37(6):641-653.

ZOU C N, DONG D Z, WANG S J, et al. Geological characteristics formation mechanism and resource potential of shale gas in China[J].Petroleum Exploration and Development,2010,37(6):641-653.

[6]
邹才能, 朱如凯, 吴松涛, 等. 常规与非常规油气聚集类型、特征、机理及展望——以中国致密油和致密气为例[J]. 石油学报, 2012, 33(2):173-187.

ZOU C N, ZHU R K, WU S T, et al. Types, characteristics, genesis and prospects of conventional and unconventional hydrocarbon accumulations:Taking tight oil and tight gas in China as an instance[J].Acta Petrolei Sinica,2012,33(2):173-187.

[7]
邓泽,王红岩,姜振学,等.页岩和煤岩的孔隙结构差异及其天然气运移机理[J].天然气工业,2022,42(11):37-49.

DENG Z, WANG H, JIANG Z X, et al. Pore structure differences between shale and coal and their gas migration mechanisms[J]. Natural Gas Industry, 2022, 42(11): 37-49.

[8]
焦堃, 姚素平, 吴浩, 等. 页岩气储层孔隙系统表征方法研究进展[J]. 高校地质学报, 2014, 20(1):151-161.

JIAO K, YAO S P, WU H, et al. Advances in characterization of pore system of gas shales[J]. Geological Journal of China Universities, 2014, 20(1):151-161.

[9]
孙寅森, 郭少斌. 基于图像分析技术的页岩微观孔隙特征定性及定量表征[J]. 地球科学进展, 2016, 31(7):751-763.

SUN Y S, GUO S B. Qualitative and quantitative characterization of shale microscopic pore characteristics based on image analysis technology[J]. Advances in Earth Science, 2016, 31(7):751-763.

[10]
曹涛涛, 宋之光, 刘光祥, 等. 氮气吸附法—压汞法分析页岩孔隙、分形特征及其影响因素[J]. 油气地质与采收率, 2016, 23(2):1-8.

CAO T T,SONG Z G,LIU G X,et al.Characteristics of shale pores,fractal dimension and their controlling factors determined by nitrogen adsorption and mercury injection methods[J].Petroleum Geology and Recovery Efficiency,2016,23(2):1-8.

[11]
刘苗苗, 付小平, 谢佳彤. 四川盆地中侏罗统页岩储层特征及其影响因素[J].断块油气田, 2023, 30(6):905-913.

LIU M M, FU X P, XIE J T. The reservoirs characteristics and its influencing factors of Middle Jurassic shale in Sichuan Basin[J].Fault-Block Oil & Gas Field,2023,30(6):905-913.

[12]
陈尚斌, 朱炎铭, 王红岩, 等. 川南龙马溪组页岩气储层纳米孔隙结构特征及其成藏意义[J]. 煤炭学报, 2012, 37(3):438-444.

CHEN S B, ZHU Y M, WANG H Y, et al. Structure characteristics and accumulation significance of nanopores in Longmaxi shale gas reservoir in the southern Sichuan Basin[J].Journal of China Coal Society, 2012, 37(3):438-444.

[13]
彭艳霞, 杜玉山, 蒋龙, 等.济阳坳陷缓坡带页岩油储层微观孔隙结构及分形特征[J]. 断块油气田,2023,30(4):535-544.

PENG Y X, DU Y S, JIANG L, et al. Micropore structure and fractal characteristics of shale oil reservoir in gentle slope zone of Jiyang Depression[J]. Fault-Block Oil & Gas Field, 2023, 30(4):535-544.

[14]
付常青, 朱炎铭, 陈尚斌. 浙西荷塘组页岩孔隙结构及分形特征研究[J]. 中国矿业大学学报, 2016, 45(1):77-86.

FU C Q, ZHU Y M, CHEN S B. Pore structure and fractal features of Hetang Formation shale in western Zhejiang[J].Journal of China University of Mining & Technology,2010,45(1):77-86.

[15]
张喜淳,胡晨林,田继军,等. 川南ZX向斜五峰组—龙马溪组页岩孔隙特征及差异性[J].断块油气田,2022,29(4):469-474.

ZHANG X C,HU C L,TIAN J J,et al.Characteristics and differences of shale pores in Wufeng-Longmaxi Formation of ZX syncline,southern Sichuan Basin[J].Fault-Block Oil & Gas Field,2022,29(4):469-474.

[16]
朱文博, 张训华, 周道容, 等. 下扬子地区二叠系海相页岩孔隙特征新认识及页岩气勘探启示[J]. 天然气工业, 2021, 41(7):41-55.

ZHU W B, ZHANG X H, ZHOU D R,et al.New cognition on pore structure characteristics of Permian marine shale in the Lower Yangtze Region and its implications for shale gas exploration[J]. Natural Gas Industry, 2021, 41(7):41-55.

[17]
BENOIT M.How long is the coast of britain?statistical selfsi-milarity and fractional dimension[J].Science,1967,156:636-638.

[18]
陈居凯, 朱炎铭, 崔兆帮, 等. 川南龙马溪组页岩孔隙结构综合表征及其分形特征[J]. 岩性油气藏, 2018, 30(1):55-62.

CHEN J K, ZHU Y M, CUI Z B, et al. Pore structure and fractal characteristics of Longmaxi shale in southern Sichuan Basin[J].Lithologic Reservoirs,2018,30(1):55-62.

[19]
黄金亮, 董大忠, 李建忠, 等. 陆相页岩储层孔隙分形特征——以四川盆地三叠系须家河组为例[J]. 天然气地球科学, 2016, 27(9):1611-1618,1708.

HUANG J L, DONG D Z, LI J Z, et al. Reservoir fractal characteristics of continental shale:An example from Triassic Xujiahe Formation shale, Sichuan Basin,China[J].Natural Gas Geoscience, 2016, 27(9):1611-1618,1708.

[20]
YANG F, NING Z F, LIU H Q. Fractal characteristics of shales from a shale gas reservoir in the Sichuan Basin,China[J]. Fuel, 2014, 115:378-384.

[21]
LIU X J. XIONG J, LIANG L X. Investigation of pore structure and fractal characteristics of organic-rich Yanchang Formation shale in Central China by nitrogen adsorption/desorption analysis[J]. Plateau Meteorology, 2015, 22:62-72.

[22]
FU H J, YAN D T, YAO C P, et al. Pore structure and multi-scale fractal characteristics of adsorbed pores in marine shale:A case study of the Lower Silurian Longmaxi Shale in the Sichuan Basin,China[J].Journal of Earth Science, 2022, 33(5):1278-1290.

[23]
XI Z D, TANG S H, ZHANG S H, et al. Pore structure characteristics of marine-continental transitional shale:A case study in the Qinshui Basin,China[J].Energy &Fuels,2017,31(8):7701-8822.

[24]
闫高原, 张军建, 路冠文, 等. 沁水盆地太原组—山西组页岩孔隙分形特征[J]. 新疆石油地质, 2021, 42(5):548-553.

YAN G Y, ZHANG J J, LU G W, et al. Fractal characteristics of shale pores from Taiyuan Formation to Shanxi Formation in Qinshui Basin[J]. Xingjiang Petroleum Geology, 2021, 42(5):548-553.

[25]
郗兆栋, 唐书恒, 李俊, 等. 沁水盆地中东部海陆过渡相页岩孔隙结构及分形特征[J]. 天然气地球科学, 2017, 28(3):366-376.

XI Z D, TANG S H, LI J, et al. Investigation of pore structure and fractal characteristics of marine continental transitional shale in the east-central of Qinshui Basin[J]. Natural Gas Geoscience, 2017, 28(3):366-376.

[26]
李小明, 柳吉荣, 吝文, 等. 荆门探区五峰组—龙马溪组页岩储层特征及甜点层段评价[J]. 煤田地质与勘探, 2021, 49(6):1-11.

LI X M, LIU J R, LIN W, et al. Characteristics of the shale gas reservoirs and evaluation of sweet spots in Wufeng Formation and Longmaxi Formation in Jingmen exploration area[J]. Coal Geology & Exploration, 2021, 39(1):29-45.

[27]
魏建光,唐书恒,张松航,等. 宁武盆地山西组过渡相页岩孔隙特征及影响因素[J].煤田地质与勘探,2018,46(1):78-85.

WEI J G,TANG S H,ZHANG S H,et al.Analysis on characteristics and influence factors of transitional facies shale pore in Ningwu Basin[J].Coal Geology & Exploration,2018,46(1):78-85.

[28]
邓恩德,姜秉仁,高为,等.黔西地区龙潭组煤系泥页岩孔隙结构及分形特征研究[J].煤炭科学技术,2020,48(8):184-190.

DENG E D,JIANG B R,GAO W,et al. Study on pore structure and fractal characteristics of shale from coal measures of Longtan Formation in western Guizhou[J].Coal Science and Technology, 2020, 48(8):184-190.

[29]
MA X, GUO S B, SHI D S, et al. Investigation of pore structure and fractal characteristics of marine-continental transitional shales from Longtan Formation using MICP, gas adsorption, and NMR (Guizhou, China)[J]. Marine and Petroleum Geology, 2019, 107:555-571.

[30]
谢卫东, 王猛, 王华, 等. 海陆过渡相页岩气储层孔隙多尺度分形特征[J]. 天然气地球科学, 2022, 33(3):451-460.

XIE W D, WANG M, WANG H, et al. Multi-scale fractal characteristics of pores in transitional shale gas reservoir[J]. Natural Gas Geoscience, 2022, 33(3):451-460.

[31]
安成, 柳广弟, 孙明亮, 等. 基于氮气吸附实验与分形FHH模型分析页岩孔隙结构特征——以鄂尔多斯盆地华池地区长7段为例[J]. 石油实验地质, 2023, 45(3):576-586.

AN C, LIU G D, SUN M L, et al. Analysis of shale pore structure characteristics based on nitrogen adsorption experiment and fractal FHH model:A case study of 7th member of Triassic Yanchang Formation in Huachi area, Ordos Basin[J]. Petroleum Geology & Experiment, 2023, 45(3):576-586.

[32]
张岩, 刘金城, 徐浩, 等. 陆相与过渡相煤系页岩孔隙结构及分形特征对比——以鄂尔多斯盆地东北缘延安组与太原组为例[J]. 石油学报, 2017, 38(9):1036-1046.

ZHANG Y, LIU J C, XU H, et al. Comparison between pore structure and fractal characteristics of continental and transitional coal measures shale:A case study of Yan’an and Taiyuan formations at the northeastern margin of Ordos Basin[J]. Acta Petrolei Sinica, 2017, 38(9):1036-1046.

[33]
彭艳霞. 南华北盆地太原组页岩气形成条件研究[D]. 北京:中国地质大学(北京), 2020.

PENG Y X. The Shale Gas Accumulation Conditions of Taiyuan Formation in the Southern North China Basin[D].Beijing:China University of Geosciences (Beijing), 2020.

[34]
何明喜. 南华北多期复合叠合盆地与油气[M]. 北京:地质出版社, 2012.

HE M X. Multi-Phase Composite Stacked Basins and Hydrocarbons in North and South China[M].Beijing:Geological Publishing House,2012.

[35]
ZHANG S W, SUI F G, LIN H X, et al. Petroleum Geology and Perspective Evaluation on Pre-Paleogene System in Bohai Bay Basin[M].Beijing:Geological Publishing House,2009:124-128.

[36]
董大忠, 邱振, 张磊夫, 等. 海陆过渡相页岩气层系沉积研究进展与页岩气新发现[J]. 沉积学报, 2021, 49(6):1-11.

DONG D Z, QIU Z, ZHANG L F, et al. Progress on sedimentology of transitional facies shales and new discoveries of shale gas[J].Acta Sedimentologica Sinica,2021,39(1):29-45.

[37]
JIA S X, ZHANG C K, ZHAO J R, et al. Crustal structure of the rift-depression basin and Yanshan Uplift in the Northeast Part of North China[J].Chinese Journal of Geophysics,2009, 52(1):51-63.

[38]
李澎,李正,张爱印,等.开平向斜海陆过渡相煤系气成藏特征及勘探方向[J/OL].特种油气藏,2025,(2):1-12[2025-06-19].http://kns.cnki.net/kcms/detail/21.1357.TE.20241104.1405.005.html.

LI P, LI Z, ZHANG A Y, et al. Accumulation characteristics and exploration direction of marine-continental transitional facies coal measure gas in Kaiping syncline[J/OL]. Special Oil & Gas Reservoirs,2025,(2):1-12[2025-06-19]. http://kns.cnki.net/kcms/detail/21.1357.TE.20241104.1405.005.html.

[39]
郑宇超, 林小兵, 刘莉萍, 等. 石英颗粒边界溶蚀类型特征及成因探讨[J/OL]. 沉积学报,1-18[2025-06-19].https://doi.org/10.14027/j.issn.1000-0550.2023.103.

ZHENG Y C,LIN X B,LIU L P,et al.Discussion of the charac-teristics and causes of different types of quartz grain boundary dissolution[J].Acta Sedimentologica Sinica,2023,1-18[2025-06-19].https://doi.org/10.14027/j.issn.1000-0550.2023.103.

[40]
张吉振,韩登林,林伟,等.松辽盆地白垩系嫩江组煤系页岩孔隙结构及分形特征[J].天然气地球科学,2024,35(1):119-132.

ZHANG J Z,HAN D L,LIN W, et al. Pore structure and fractal characteristics of coal-bearing Cretaceous Nenjiang shales from Songliao Basin, Northeast China[J]. Natural Gas Geoscience, 2024, 35(1):119-132.

[41]
宋泽章, 阿比德·阿不拉, 吕明阳, 等. 氮气吸附滞后回环定量分析及其在孔隙结构表征中的指示意义——以鄂尔多斯盆地上三叠统延长组7段为例[J]. 石油与天然气地质, 2023, 44(2):495-509.

SONG Z Z, ABIDE·A, LYU M Y, et al. 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[J]. Oil & Gas Geology, 2023, 44(2):495-509.

[42]
GROEN J C, PEFFER L A, PEREZ-RAMIREZ J. Pore size determination in modified micro and mesoporous materials. Pitfalls and limitations in gas adsorption data analysis[J]. Microporous Mesoporous Mater, 2003, 60 (1-3);1-17.

[43]
徐勇, 吕成福, 陈国俊, 等. 川东南龙马溪组页岩孔隙分形特征[J]. 岩性油气藏, 2015, 27(4):32-39.

XU Y, LYU C F, CHEN G J, et al. Fractal characteristics of shale pores of Longmaxi Formation in southeast Sichuan Basin[J]. Lithologic Reservoirs, 2015, 27(4):32-39.

[44]
张闯辉, 朱炎铭, 刘宇, 等. 不同成熟度页岩孔隙及其分形特征[J]. 断块油气田, 2016, 23(5):583-588.

ZHANG C H, ZHU Y M, LIU Y, et al. Pore and fractal characteristics of shale in different maturity[J]. Fault-block Oil & Gas Field, 2016, 23(5):583-588.

[45]
赵越, 李磊, 司运航, 等. 浅层页岩气储层孔隙分形特征及控制因素——以云南昭通地区龙马溪组为例[J]. 吉林大学学报(地球科学版), 2022, 52(6):1813-1829.

ZHAO Y,LI L,SI Y H,et al.Fractal characteristics and contro-lling factors of pores in shallow shale gas Reservoirs:A case stu-dy of Longmaxi Formation in Zhaotong area,Yunnan Province[J].Journal of Jilin University(Earth Science Edition),2022, 52(6):1813-1829.

[46]
何陈诚, 赵永强, 俞凌杰, 等. 川东北地区二叠系大隆组深层页岩气储层孔隙结构及其分形特征[J]. 石油实验地质, 2024, 46(2):263-277.

HE C C,ZHAO Y Q,YU L J,et al. Pore structure and fractal characteristics of deep shale gas reservoirs in the Permian Dalong Formation, northeastern Sichuan Basin[J].Petroleum Geology & Experiment, 2024, 46(2):263-277.

[47]
管全中, 董大忠, 孙莎莎, 等. 深层富有机质页岩孔隙结构分形特征及其地质意义——以四川盆地威远地区下志留统龙马溪组为例[J]. 天然气工业, 2024, 44(3):108-118.

GUAN Q Z, DONG D Z, SUN S S, et al. Fractal characteristics of organic-rich shale pore structure and its geological implications:A case study of the Lower Silurian Longmaxi Formation in the Weiyuan block, Sichuan Basin[J]. Natural Gas Industry, 2024, 44(3):108-118.

[48]
KEWEN L. Analytical derivation of Brooks⁃Corey type capillary pressure models using fractal geometry and evaluation of rock heterogeneity[J]. Journal of Petroleum Science and Engineering, 2010, 73:20-26.

[49]
PFEIFERPER P, AVNIR D. Chemistry in nonintegral dimensions between two and three[J].The Journal of Chemical Physics, 1983, 79(7):3369-3558.

[50]
腾格尔, 卢龙飞, 俞凌杰, 等. 页岩有机质孔隙形成、保持及其连通性的控制作用[J]. 石油勘探与开发, 2021,48(4):687-699.

TENG G E, LU L F, YU L J, et al. Formation, preservation and connectivity control of organic pores in shale[J]. Petroleum Exploration & Development,2021,48(4):687-699.

[51]
吴伟, 梁志凯, 郑马嘉, 等.页岩储层孔隙结构与分形特征演化规律[J]. 油气地质与采收率, 2022, 29(4):35-45.

WU W,LIANG Z K,ZHENG M J,et al.Pore structures in shale reservoirs and evolution laws of fractal characteristics[J].Petroleum Geology and Recovery Efficiency,2022,29(4):35-45.

[52]
谷渊涛, 李晓霞, 万泉, 等.泥页岩有机质孔隙差异特征及影响因素分析——以我国典型海相、陆相、过渡相储层为例[J]. 沉积学报, 2021, 39(4):794-810.

GU Y T,LI X X,WANG Q,et al. On the different characteristics of organic pores in shale and their influencing factors:Taking typical marine,continental,and transitional facies reservoirs in China as examples[J].Acta Sedimentologica Sinica,2021,39(4):794-810.

[53]
郭少斌, 王子龙, 马啸. 中国重点地区二叠系海陆过渡相页岩气勘探前景[J]. 石油实验地质,2021, 43(3):377-385,414.

GUO S B,WANG Z L,MA X.Exploration prospect of shale gas with Permian transitional facies of some key areas in China[J].Petroleum Geology & Experiment,2021,43(3):377-385,414.

[54]
NING S,XIA P,HAO F,et al.Pore fractal characteristics between marine and marine-continental transitional black shales:A case study of Niutitang Formation and Longtan Formation[J].Fractal and Fractional, 2024, 8(5):288.

[55]
董立朋, 聂清浩, 孙晓坤, 等. 基于皮尔逊相关系数法的盾构掘进参数对地表沉降影响分析[J]. 施工技术,2024, 53(1):116-123.

DONG L P, NIE Q H, SUN X K, et al. Analysis of impact of shield tunneling parameters on ground settlement based on pearson correlation coefficient method[J]. Construction Technology, 2024,53(1):116-123.

[56]
王大兴, 曾治平, 胡海燕, 等. 准噶尔盆地中部下乌尔禾组深层陆相页岩孔隙结构分形特征及其地质意义[J]. 油气地质与采收率, 2024, 31(1):23-35.

WANG D X, ZENG Z P, HU H Y, et al. Fractal characteristics of pore structure of deep continental shale of Lower Wuerhe Formation in central Junggar Basin and its geological significance[J]. Petroleum Geology and Recovery Efficiency, 2024, 31(1):23-35.

[57]
王伟庆, 王学军, 李政, 等. 东营凹陷低成熟页岩储集空间与分形特征——以民丰洼陷沙四段上亚段页岩为例[J].油气地质与采收率,2025,32(2):14-24.

WANG W Q, WANG X J, LI Z, et al. Reservoir space and fractal characteristics of low-maturity shale in Dongying Depression:A case study of upper submember of the fourth member in Shahejie Formation in the Minfeng Area[J]. Petroleum Geology and Recovery Efficiency,2025,32(2):14-24.

[58]
吉利明, 邱军利, 夏燕青, 等. 常见黏土矿物电镜扫描微孔隙特征与甲烷吸附性[J]. 石油学报, 2012, 33(2):249-256.

JI L M, QIU J L, XIA Y Q, et al. Micro-pore characteristics and methane adsorption properties of common clay minerals by electron microscope scanning[J]. Acta Petrolei Sinica, 2012, 33(2):249-256.

[59]
曹江骏, 陈朝兵, 罗静兰, 等. 自生黏土矿物对深水致密砂岩储层微观非均质性的影响——以鄂尔多斯盆地西南部合水地区长6油层组为例[J]. 岩性油气藏, 2020, 32(6):36-49.

CAO J J, CHEN C B, LUO J L, et al. Impact of authigenic clay minerals on micro-heterogeneity of deep water tight sandstone reservoirs:A case study of Triassic Chang 6 oil reservoir in Heshui area, southwestern Ordos Basin[J]. Lithologic Reservoirs, 2020, 32(6):36-49.

[60]
YANG B, QU H, PU R, et al. Controlling effects of tight reservoir micropore structures on seepage ability:A case study of the Upper Paleozoic of the eastern Ordos Basin, China[J]. Acta Geologica Sinica-English Edition, 2020, 94(2):322-336.

[61]
CHENG L,GUAN F,LIU D,et al. Micropore characteristics and geological significance of shale reservoirs:Case study of Fuling Shale Gas in Sichuan Basin,China[J].Geofluids,2021,2021(1): 6636156.

Outlines

/