Fractal characteristics of pore throat structure of Chang 8 Member tight sandstone reservoir in Wuqi-Zhidan area, Ordos Basin

  • Bin LUO , 1, 2 ,
  • Rongde TAO 3 ,
  • Xiaoli CUI 1, 2 ,
  • Honggang XIN 1, 2 ,
  • Weidong DAN 1, 2
Expand
  • 1. Research Institute of Exploration and Development,PetroChina Changqing Oilfield Company,Xi'an 710018,China
  • 2. National Engineering Laboratory for Exploration and Development of Low Permeability Oil and Gas Fields,Xi'an 710018,China
  • 3. The 8th Oil Production Plant of Changqing Oilfield Company,PetroChina,Xi'an 710021,China

Received date: 2022-02-19

  Revised date: 2023-02-23

  Online published: 2023-07-28

Supported by

The Forward-looking Fundamental Major Science and Technology Objects in the Upstream Field of PetroChina during the 14th Five Year Plan Period(2021DJ0402)

Highlights

In order to accurately characterize the pore throat structure of tight reservoir in the Chang 8 Member of the Wuqi-Zhidan area of the Ordos Basin, based on the analyses of cast thin section, scanning electron microscopy, X-ray diffraction, high-pressure mercury injection, constant velocity mercury injection, and oil-water phase infiltration, the pore throat fractal characteristics and throat fractal characteristics of tight reservoirs in Chang 8 Member were studied by using fractal theory. The relationship between fractal characteristics and reservoir porosity, permeability, and microscopic pore throat characteristics was systematically analyzed. The results show that the lithology of the study area is dominated by lithic feldspar sandstone, with an average porosity of 8.2% and an average permeability of 0.16×10-3 μm2. The reservoirs are generally tight. The pore types of the reservoir are mainly intergranular pores and dissolution pores. According to the morphology of capillary pressure curve and high-pressure mercury injection parameters, the pore throat structure of the reservoir is divided into two types: type I and type II. The constant rate mercury intrusion experiment showed that the mercury saturation of the samples in the study area was mainly contributed by the throat. The fractal dimension of throat calculated by constant velocity mercury injection has two characteristics:Fine throat section D T1 and coarse throat section D T2. The full aperture fractal dimension and throat fractal dimension have good negative correlation with physical properties, mercury injection saturation and mercury removal efficiency. The throat fractal dimension has a good negative correlation with the mercury injection saturation of the throat and the radius of the mainstream throat. The more homogeneous the pore throat is, the more difficult it is for the oil phase fluid to pass through the water-bearing reservoir and the lower the relative permeability of the oil phase.

Cite this article

Bin LUO , Rongde TAO , Xiaoli CUI , Honggang XIN , Weidong DAN . Fractal characteristics of pore throat structure of Chang 8 Member tight sandstone reservoir in Wuqi-Zhidan area, Ordos Basin[J]. Natural Gas Geoscience, 2023 , 34(8) : 1305 -1315 . DOI: 10.11764/j.issn.1672-1926.2023.03.012

0 引言

鄂尔多斯盆地吴起—志丹地区延长组8段(简称长8段)资源潜力巨大,是该地区下一步增储上产的重要接替层系,对陕北油区实现长期稳产具有重要的现实意义。研究区长8段储层普遍致密,孔隙度主要分布在6%~12%之间,渗透率分布在(0.1~0.3)×10-3 μm2之间,是典型的致密储层。致密储层具有复杂的孔喉结构特征,相较于孔喉大小,孔喉结构的复杂程度对于表征孔喉结构也至关重要。已有大量学者运用多种研究手段对孔喉结构的复杂程度进行了研究,认为孔喉结构的复杂程度受孔喉半径大小、孔喉的配置关系以及岩石矿物组成等多种因素耦合控制。自1986年MANDELBRO1提出分形理论以来,已逐渐成为系统科学研究中的前沿领域,分数维度站在数学角度上进行描述和研究,展现事物中的规律性、复杂性和自相似性。由于分形维数在研究多孔介质上具有得天独厚的优势,所以YUE等2利用分形理论提出了混凝土—岩石界面过渡区的渗透率和孔隙结构数学模型;薛帆等3以北洛河流域为研究区,基于Budyko假设和分形理论,采用弹性系数法对研究区3种不同地貌和植被类型区的水砂变化进行归因识别;ZHAO等4基于分形维数来表征添加再生废砖微粉复合土的干燥开裂特性;HUANG等5基于分形理论研究岩石内部孔喉结构,定量表征储层孔喉结构的复杂程度。砂岩储层分形维数的计算方法多样,其中,利用高压压汞、恒速压汞以及核磁共振等实验结果计算储层孔喉结构分形维数是常用的方法6-7。而对于致密储层而言,孔喉大小以及喉道分布对储层的储集性能具有显著影响8-10,因此,本文通过对研究区长8段致密砂岩储层开展铸体薄片、扫描电镜、X射线衍射、高压压汞、恒速压汞等测试,在明确储层的孔喉类型及特征基础上,结合高压压汞和恒速压汞对全孔径分形特征和喉道分形特征的表征,探讨孔喉结构分形特征与物性、压汞参数、渗流的关系,为研究区长8段油藏的效益开发提供理论依据。

1 区域地质概况

鄂尔多斯盆地位于我国中东部,是一个中、新生代稳定沉降、大型的内陆汇水盆地。上三叠统延长组是在鄂尔多斯盆地坳陷持续发展和稳定沉降过程中沉积的以河流—湖泊相为特征的陆源碎屑岩系,是盆地重要的致密砂岩储层发育层位11。吴起—志丹地区位于鄂尔多斯盆地中部的陕北斜坡(图1),长8段是该地区主要的含油气层系,根据沉积旋回可将长8段划分为长81、长82共2个亚段。研究区长8段沉积属于北东物源体系,为三角洲前缘亚相沉积,储层岩石类型以岩屑长石砂岩为主,其次为长石岩屑砂岩,还有少量的长石砂岩。砂岩碎屑成分中石英占比13%~66%,长石(主要是斜长石)占比17%~77%,岩屑占比4%~42%,岩屑类型主要以变质岩岩屑为主,沉积粒度细,总体以细砂为主,分选中等—好,磨圆主要为次棱角状—次圆状。储层平均孔隙度为8.2%,平均渗透率为0.16×10-3 μm2,属于典型的低孔超低渗透储层。
图1 研究区区域位置

Fig.1 Location of the study area

2 孔喉类型及特征

铸体薄片、扫描电镜等岩相学分析表明,研究区长8段致密砂岩样品的面孔率平均为2.08%,孔隙类型以粒间孔和溶蚀孔为主,晶间孔、微孔隙不发育。其中,粒间孔分布不均一,非均质性强,多为压实或胶结残余孔隙,呈三角形或多边形,孔隙边缘平直,内部洁净,无明显的溶蚀痕迹,填充孔隙的多为薄膜式胶结的绿泥石和铁方解石等胶结物。扫描电镜分析显示,原生孔隙周围常见绿泥石膜,有利于原生孔隙的保存,有时粒间充填黏土矿物,可见高岭石伊利石化,激光共聚焦显微镜下原生粒间孔边缘平直部分被绿泥石膜包裹,一般分布在颗粒点接触之间(图2);溶蚀孔主要是由长石、岩屑、杂基和其他胶结物等易溶矿物在酸性溶液作用下溶解形成孔隙形态复杂的次生孔隙(图3)。
图2 吴起—志丹地区原生孔隙特征

(a)残余粒间孔发育,X410井,1 994.70 m,长81亚段;(b)残余粒间孔,方解石胶结物发育,X410井,2 022.80 m,长82亚段;(c)残余粒间孔,X410井,2 023.45 m,长82亚段;(d)残余粒间孔被自生石英充填,发育绿泥石薄膜,G6井,1 970.40 m,长81亚段;(e)粒间伊利石化的高岭石,G6井,2 001.43 m,长82亚段;(f)粒间孔隙周围绿泥石膜发育,G6井,2 001.43 m,长82亚段;(g)三角形粒间孔发育,G6井,1 971.20 m,长81亚段;(h)粒间孔边缘平直部分被绿泥石膜包裹,X522井,2 031.95 m,长81亚段;(i)粒间孔隙发育,X522井,2 031.95 m,长81亚段

Fig.2 Characteristics of primary pores in Wuqi-Zhidan area

图3 吴起—志丹地区次生孔隙特征

(a)长石边缘发生溶蚀,G6井,1 967.65 m,长81亚段;(b)被溶蚀的长石颗粒边缘常呈港湾状或不规则状,X522井,2 031.95 m,长81亚段;(c)长石边缘发生溶蚀形成粒间溶扩孔隙,X410井,1 994.15 m,长81亚段;(d)孔隙中溶蚀残余发育,G6井,1 967.15 m,长81亚段;(e)长石溶孔中充填自生石英,G6井,1 970.40 m,长81亚段;(f)长石溶孔中充填高岭石,G6井,2 001.43 m,长82亚段;(g)溶孔非均质性分布,G6井,1 970.40 m,长81亚段;(h)溶蚀形成纳米级孔喉,G6井,1 967.15 m,长81亚段;(i)溶蚀孔发育,G6井.1 967.15 m,长81亚段

Fig.3 Characteristics of secondary pores in Wuqi-Zhidan area

3 储层微观孔喉结构特征

3.1 高压压汞表征孔喉结构

大量的高压压汞资料统计结果表明,排驱压力平均值为1.6 MPa,进汞饱和度平均值为85.5%,孔喉半径小,主要分布在0.01~1 μm之间,中值半径平均值为0.07 μm。分选整体一般,分选系数为1.78~3.12,平均值为2.31。连通性中等,最大进汞饱和度分布在73%~90%之间。
以毛管压力曲线特征为主,结合排驱压力、最大进汞饱和度、退汞效率、中值压力、中值半径以及分选系数等孔喉结构等参数进行综合研究,将研究区长8段储层孔喉结构可划分为I类和II类2类(图4表1)。I类孔喉结构的进汞曲线整体分布在偏左下方,排驱压力较低,平均为1.1 MPa,中值压力平均为5.71 MPa,中值半径平均为0.13 μm,相较于II类,这类样品孔喉较大,孔喉配伍性较好,分选系数平均为1.9,最大进汞饱和度表现出较高的特征,平均为88.2%,退汞效率好,连通性好,整体上属于相对较好的储层,孔隙度平均为9.94%,渗透率平均为0.28×10-3 μm2。II类孔喉结构的进汞曲线逐渐偏向右上方,中值半径、最大进汞饱和度、退汞效率较I类低,而排驱压力、中值压力较I类高,储层物性低、品质差。
图4 高压压汞分类

Fig.4 High-pressure mercury intrusion classification

表1 高压压汞分类参数

Table 1 High-pressure mercury intrusion classification parameter

曲线类型 孔隙度/% 渗透率 /(10-3 μm2 排驱压力/MPa 中值压力/MPa 中值半径/μm 分选系数 最大进汞饱和度 /%

退汞效率

/%

I 9.94 0.28 1.1 5.71 0.13 1.9 88.2 27.68
II 8.2 0.16 3.5 16.52 0.04 2.0 81.4 24.11

3.2 恒速压汞表征孔喉结构

研究区长8段不同样品的孔隙半径分布区间较为相似,均分布于70~280 μm之间,孔隙半径平均值集中于131.87~148.32 μm之间;喉道半径分布于0.10~1.90 μm之间,喉道半径平均值集中于0.17~0.79 μm之间(表2)。
表2 恒速压汞参数

Table 2 Constant speed mercury injection parameter

井号 孔隙度/%

渗透率

/(10-3 μm2

排驱压力/MPa 最终进汞饱和度/% 总孔隙进汞饱和度/% 总喉道进汞饱和度/% 喉道半径平均值/μm 孔隙半径平均值/μm 孔喉半径比平均值
W484 14.40 0.40 0.40 59.00 23.71 35.29 0.79 140.96 215.16
X287 10.10 0.28 0.88 58.31 28.48 29.84 0.60 148.32 282.93
X420 9.00 0.17 0.94 60.33 28.68 31.64 0.47 143.84 367.54
G124 11.30 0.20 1.18 42.94 17.34 25.60 0.51 143.40 327.19
X268 8.30 0.22 1.43 48.52 12.73 35.79 0.47 141.12 374.56
X297 7.40 0.14 1.64 28.90 3.06 25.84 0.44 135.70 321.23
X71 9.0 0.10 2.13 16.79 1.42 15.37 0.30 132.77 438.19
X372 9.6 0.13 1.71 21.37 3.86 17.51 0.37 139.83 395.59
X248 8.9 0.16 2.55 12.75 0.34 12.41 0.27 131.87 486.67
X414 11.4 0.07 4.07 15.03 1.22 13.81 0.17 132.64 780.53
通过观察总进汞饱和度、喉道进汞和孔隙进汞的恒速压汞曲线,将毛管压力曲线形态可以分成孔隙区、孔喉区和喉道区12-13图5)。从图5可以看出,随着毛管压力的增加,进汞饱和度上升,在孔隙区,喉道进汞为零,总进汞饱和度只由孔隙进汞饱和度贡献;在孔喉区时,喉道开始进汞,总进汞饱和度由孔隙和喉道同时贡献;而在喉道区,毛管压力不断上升,此时孔隙进汞饱和度进入平台期,总进汞饱和度的增量为喉道进汞增量。研究区整体物性差,孔喉结构复杂,没有出现以孔隙型为主的恒速压汞曲线。
图5 恒速压汞毛管压力曲线示意

Fig.5 Schematic diagram of capillary pressure curve of constant speed mercury injection

研究区长8段储层大多为过渡型恒速压汞曲线[图6(a)],曲线上基本没有孔隙区,随着样品进汞压力逐渐增大,孔隙和喉道的进汞量同时上升,但孔隙进汞量高于喉道的进汞量,此时为孔喉区,当孔隙进汞量不再增加时进入喉道区,孔隙进汞量的曲线趋于平缓为平台期,而喉道进汞量随着毛管压力的增大持续增加,最终进汞量超过孔隙进汞量。而喉道型压汞曲线[图6(b)]与过渡型恒速压汞曲线相同,没有孔隙区,在孔喉区,随着进汞压力的增大,孔隙的进汞量微弱,远小于喉道的进汞量,到了孔隙进汞量的平台期时,此时喉道进汞量快速增加,最终远超于孔隙进汞量,喉道所占的体积要大于孔隙体积。
图6 研究区恒速压汞毛管压力曲线

(a)W484井;(b)X268井

Fig.6 Capillary pressure curve of constant velocity mercury injection in the study area

从过渡型恒速压汞曲线到喉道型压汞曲线,岩样的物性逐渐降低,恒速压汞参数除排驱压力以外,其余参数都表现出下降趋势,过渡型喉道孔隙分布更加均匀,整体孔喉配置更好,喉道型压汞曲线的样品整体由喉道主导,孔喉分布复杂,喉道在很大程度上影响了渗流能力。

3.3 分形维数表征孔喉结构

3.3.1 分形理论

岩石内部结构作为一个立体的空间,拥有复杂的孔隙、喉道大小及形态,分形理论可以有效表征储层孔喉结构的复杂程度,是连接储层孔喉大小、形状以及孔喉连通性和宏观物性特征之间的桥梁14-15
根据分形理论可推导出计算致密砂岩储层微观孔喉结构的分形维数数学模型16-18
L g 1 - S H g r = 3 - D L g r - 3 - D L g r m a x
式(1)可以看出,3-D为Lg[1-S Hg(≥r)]与Lgr拟合的曲线的斜率。即:
D=3-k
式中:r为孔喉半径,μm;S Hg为当孔喉半径为r时的进汞饱和度,%;k为Lg[1-S Hg(≥r)]与Lgr的相关关系的斜率;D为对应的分形维数。
分形维数值越接近3,说明孔喉结构越复杂,非均值性越强,表面越粗糙,越靠近2则说明孔喉表面较光滑,孔喉结构非均值性弱10

3.3.2 全孔径分形特征

通过高压压汞实验数据,绘制各岩样的Lg(1-S Hg)和Lgr的散点图,得到高压压汞实验的孔喉分形维数计算图(图7),由于高压压汞提供的毛管压力范围是从0到200.25 MPa,可表征孔喉半径范围是最大连通孔喉至3.67 nm,所以可用来表征尺度相对较小的孔喉分布特征。
图7 X297井相对较小孔喉分形特征

Fig.7 Fractal characteristics of relatively small pore throat in Well X297

恒速压汞可以表征大于63 μm的孔喉分布特征,可用来表征相对较大的孔喉的复杂程度,与高压压汞结合可以近似表征全孔径的分形特征13图8)。
图8 W484井相对较大孔喉分形特征

Fig.8 Fractal characteristics of relatively large pore throat in Well W484

通过计算,从表3可以看出,研究区长8段储层相对较小的孔喉(d)分形维数为2.47~2.71,且相关关系较好,R 2值平均为0.946 8。相对较大的孔喉分形维数(D)为2.53~2.76,且相关关系较好,R 2值平均为0.969 7。其中W484井较小孔喉的分形维数为2.47,较大孔喉的分形维数为2.61,平均为2.54,孔喉结构复杂程度最低;X414井较小孔喉的分形维数为2.71,较大孔喉的分形维数为2.76,平均为2.74,孔喉结构最复杂。
表3 全孔径分形维数参数

Table 3 Full-aperture fractal dimension parameter

井号 孔隙度/%

渗透率

/(10-3 μm2

d R 2 D R 2
W484 14.40 0.40 2.47 0.894 3 2.61 0.976 9
X287 10.10 0.28 2.57 0.947 0 2.58 0.975 9
G124 11.30 0.20 2.58 0.955 2 2.62 0.956 7
X420 9.00 0.17 2.65 0.955 7 2.53 0.977 2
X268 8.30 0.22 2.61 0.950 1 2.69 0.958 8
X297 7.40 0.14 2.61 0.947 5 2.54 0.969 1
X71 9.0 0.10 2.66 0.963 4 2.69 0.987 4
X372 9.6 0.13 2.61 0.937 8 2.71 0.973 4
X248 8.9 0.16 2.63 0.954 4 2.67 0.965 4
X414 8.4 0.07 2.71 0.962 3 2.76 0.956 7

3.3.3 喉道分形特征

高压压汞和恒速压汞联合反映了孔喉的全孔径特征,而恒速压汞可以区分喉道和孔隙,利用分形维数可以在一定程度上反映孔隙和喉道的复杂程度。恒速压汞可以很好地区分喉道和孔隙,通过对致密砂岩的孔隙数据进行分形特征研究时发现,孔隙半径较大段曲线相关系数很低,致密砂岩的孔隙分形特征多为三段式(图9),而且相较于通过压力计算的喉道半径,孔隙大小是依据等效球模型得到,在恒速进汞的过程中一次压力的下降视为一个孔隙空间,但实际可能是由一个喉道连通的多个小孔隙空间共同组成,因此,孔隙数量不能单纯地认为是进汞过程中的压力下降次数,只能反映出不同大小的喉道控制的孔隙体积,这样计算的结果用来表征孔隙的复杂程度可能不太真实,而喉道大小的计算与压力有直接的对应关系,所以此次以喉道分形维数为主。
图9 恒速压汞孔隙分形维数

Fig.9 Pore fractal dimension of constant velocity mercury injection

基于恒速压汞得到的喉道数据,根据分形维数特征表明研究区长8段致密储层样品的喉道具有两段式分形结构,分段点所对应的喉道半径平均为0.45 μm,两段曲线有良好的相关性(表4图10)。曲线的左段代表小喉道段,分形维数D T1分布在2.74~2.88之间;右段为大喉道段,分形维数D T2分布在2.93~2.99之间。这可能是由于大喉道中存在附着在表面的黏土矿物,增加了表面的粗糙程度,导致大喉道结构复杂,而小喉道主要呈片状—弯片状,黏土矿物较少,喉道表面较为光滑。
表4 恒速压汞喉道分形维数统计

Table 4 Statistical of fractal dimension of constant pressure mercury throat

井号

孔隙度

/%

渗透率

/(10-3 μm2

左段 右段
D T1 R 2 D T2 R 2
W484 14.40 0.40 2.74 0.998 6 2.93 0.951 2
X287 10.10 0.28 2.78 0.991 2 2.98 0.933 0
G124 11.30 0.20 2.74 0.985 3 2.99 0.968 6
X420 9.00 0.17 2.81 0.989 7 2.94 0.908 3
X268 8.30 0.22 2.68 0.967 9 2.98 0.991 5
X297 7.40 0.14 2.80 0.989 3 2.99 0.943 8
X71 9.0 0.10 2.84 0.989 3 2.94 0.989 3
X372 9.6 0.13 2.81 0.989 3 2.94 0.989 3
X248 8.9 0.16 2.78 0.989 3 2.99 0.989 3
X414 8.4 0.07 2.88 0.989 3 2.98 0.989 3
图10 恒速压汞喉道分形维数计算曲线

(a)W484井;(b)X297井

Fig.10 Calculation curve of fractal dimension of constant velocity mercury injection throat

4 讨论

4.1 分形维数与物性之间的关系

将测试样品的物性参数和分形维数进行对比可以看出(图11),孔喉分形维数与孔隙度、渗透率之间有较好的负相关关系,即孔喉结构越复杂,造成储层的物性越差,且对渗透率的影响更大。对恒速压汞喉道分形维数和致密砂岩的物性参数进行对比(图12),孔隙度的大小并不受喉道的复杂程度的影响,而渗透率与喉道的分形维数有一定的负相关性,同时大喉道分形维数和渗透率的相关性好于小喉道分形维数,这可能是由于小喉道相对更加均质,形态多为片状、管状,而大喉道胶结物较多,喉道更加弯曲复杂,不利于流体渗流,对比高压压汞的数据计算所得的全孔径分形维数,由恒速压汞数据所计算的喉道分形维数对渗透率的影响相对较小,说明对于整个孔喉整体来说,单一喉道本身的复杂程度对渗透率的影响程度有限,而渗透率的好坏与孔喉配置也密切相关,因此孔喉整体结构的非均质程度对渗流能力的影响较大。
图11 全孔径分形维数与物性关系

Fig.11 Relationship between full-aperture fractal dimension and physical properties

图12 喉道分形维数与物性关系

Fig.12 Relationship between throat fractal dimension and physical properties

4.2 分形维数与孔喉结构参数的关系

将高压压汞得到的最大进汞饱和度和退汞效率与对应的分形维数进行对比(图13),发现孔喉结构参数均与全孔径分形维数具有较好的相关性。由于最大进汞饱和度是一个反映孔喉储集空间大小的参数,能进入多少的汞受到孔喉体积的大小及孔喉复杂程度影响,因此其大小与分形维数呈现一个负相关趋势,孔喉越粗糙,非润湿相在进汞过程中越难形成,进入孔喉的体积量也越少。而退汞效率在一定程度上可以反映孔喉的连通程度,即孔喉结构越复杂,非均质性越强,连通性越差。而喉道分形维数和喉道进汞饱和度相关性较高(图14),表面粗糙程度不同造成的毛细管力不同,喉道表面越光滑,进汞过程中受到阻力小,造成进汞饱和度越高。主流喉道半径在一定程度上可以用来表征整个岩心渗透率起主要贡献的渗流空间喉道,通过研究分形维数和主流喉道半径关系可以看出,随着喉道分形维数的增大即喉道越复杂,主流喉道半径越小,主要贡献渗透率的渗流空间也越小,喉道样品整体的孔喉非均质性越强。
图13 全孔径分形维数与孔喉结构参数的关系

Fig.13 The relationship between full aperture fractal dimension and pore throat structure parameters

图14 喉道分形维数与孔喉结构参数的关系

Fig.14 Relationship between throat fractal dimension and pore throat structure parameters

4.3 分形维数与渗流特征的关系

通过对比油水相渗实验结果可以直观判断油水两相的渗流规律,对生产过程中的见水特征有一定指导作用。对比W484、X268、X297、X414井的相渗曲线(图15),分形维数最小的W484井相渗曲线为双曲下凹型,油相下降较为缓慢,油水两相共渗区范围较大,随着分形维数逐渐增大,孔喉结构趋于复杂,相渗曲线中的等渗点逐渐降低,束缚水饱和度逐渐升高,油水两相共渗区范围逐渐缩小,残余油时对应的含水饱和度降低,油有效渗透率逐渐降低。储层岩石组分的细微差异、矿物颗粒粒度减小,压实作用的增强,胶结物的增多、溶蚀作用的发生,使储层孔喉表面粗糙程度逐渐增大,孔喉的非均质性增强,从而导致油相流体在含水储层中通过时受到的阻力较大,造成油相相对渗透率变小,从而直接影响储层的渗流能力。
图15 研究区样品油水相渗曲线

Fig.15 Oil-water phase permeability curve of samples in the study area

5 结论

(1)鄂尔多斯盆地吴起—志丹地区长8段储层孔隙结构以小孔、微细孔喉为主,孔喉中值半径平均为0.14 μm。排驱压力平均为1.41 MPa。中值压力平均为1.85 MPa。根据高压压汞曲线形态和压汞参数将储层孔喉结构分为I类、II类2种类型,其储集性能和渗流能力依次变差,孔喉结构非均质性逐渐增强。研究区样品大多为过渡型恒速压汞曲线,部分井是为喉道型恒速压汞曲线,进汞饱和度以喉道贡献为主。
(2)高压压汞与恒速压汞计算的全孔径分形维数分布在2.47~2.76之间,恒速压汞计算的喉道分形维数具有两段特征,小喉道段分形维数D T1分布在2.74~2.88之间,大喉道段分形维数D T2分布在2.93~2.99之间。全孔径分形维数和喉道分形维数均与物性参数具有良好的负相关性,全孔径分形维数与进汞饱和度和退汞效率呈良好负相关性,喉道分形维数与喉道进汞饱和度和主流喉道半径呈良好负相关性,即孔喉整体和喉道越复杂,孔喉空间越小,连通性越差,造成储层的物性越差。储层孔喉表面粗糙程度增加、孔喉连通性减弱、非均质性增强,致使油相流体在含水储层中更加难以通过,油相相对渗透率变小。
1
MANDELBROT B B. Self-affine Fractals Sets[M]//Fractals in Physics, PIETRONERO L,TOSATTI E, North-Holland, Amsterdam: Elsevier,1990:3-28.

2
YUE J, SHENG J C,WANG H M, et al. Investigation on pore structure and permeability of concrete-rock interfacial transition zones based on fractal theory[J].Fractal and Fractional,2022,6(6):329.

3
薛帆,张晓萍,张橹,等.基于Budyko假设和分形理论的水沙变化归因识别——以北洛河流域为例[J].地理学报,2022,77(1):79-92.

XUE F, ZHANG X P,ZHANG L, et al. Attribution recognition of streamflow and sediment changes based on the Budyko hypothesis and fractal theory: A case study in the Beiluo River Basin[J]. Acta Geographica Sinica,2022,77(1):79-92.

4
ZHAO X Z,LI L C,YANG B B, et al. Application of fractal to evaluate the drying shrinkage behavior of soil composites from recycled waste clay Brick[J].Fractal and Fractional,2022,7(1):25.

5
HUANG H X, CHEN L, SUN W, et al. Pore-throat structure and fractal characteristics of Shihezi Formation tight gas sandstone in the Ordos Basin,China[J]. Fractals,2018,26(2):1840005.

6
李磊,鲍志东,李忠诚,等.致密砂岩气储层微观孔隙结构与分形特征:以松辽盆地长岭气田登娄库组为例[J/OL].天然气地球科学:1-18[2023-02-16]. http://kns.cnki.net/kcms/detail/62.1177.TE.20230210.1707.002.html.

LI L, BAO Z D,LI Z C,et al Micro pore structure and fractal characteristics of tight sandstone gas reservoir:Taking Denglouku Formation of Changling Gas Field in Songliao Basin as an examples[J/OL].Natural Gas Geoscience1-18[2023-02-16].http://kns.cnki.net/kcms/detail/62.1177.TE.20230210.1707.002.html.

7
马剑,黄志龙,吴红烛,等.莺歌海盆地东方区黄流组储层微观孔喉特征及对物性的影响[J].沉积学报,2015,33(5):983-990.

MA J, HUANG Z L, WU H Z, et al. Characteristics of reservoir microscopic pores and throats and their influence on reservoir physical properties in Huangliu Formation of DF area, Yinggehai Basin[J].Acta Sedimentologica Sinica,2015,33(5):983-990.

8
任大忠.鄂尔多斯盆地延长组致密砂岩储层微观特征[D].西安:西北大学,2015.

REN D Z. Micro-characteristics of Tight Sandstone Reservoir of Yanchang Fromation in Ordos Basin[D]. Xi’an: Northwest University,2015.

9
LI P, ZHENG M, BI H, et al. Pore throat structure and fractal characteristics of tight oil sandstone: A case study in the Ordos Basin, China[J]. Journal of Petroleum Science and Engineering,2017,149:665-674.

10
WU H, ZHANG C L, JI Y L, et al. An improved method of characterizing the pore structure in tight oil reservoirs: Integrated NMR and constant-rate-controlled porosimetry data[J].Jour-nal of Petroleum Science and Engineering,2018,166:778-796.

11
杨俊杰.鄂尔多斯盆地构造演化与油气分布规律[M].北京:石油工业出版社,2002.

YANG J J.Tectonic Evolution and Hydrocarbon Distribution in the Ordos Basin[M].Beijing: Petroleum Industry Press,2002.

12
WANG Y, ZHU Y M, LIU S M, et al. Pore characterization and its impact on methane adsorption capacity for organic-rich marine shales[J]. Fuel, 2016,181:227-237.

13
欧阳思琪,孙卫,黄何鑫.多方法协同表征特低渗砂岩储层全孔径孔隙结构:以鄂尔多斯盆地合水地区砂岩储层为例[J].石油实验地质,2018,40(4):595-604.

OUYANG S Q, SUN W, HUANG H X. Multi-method synergistic characterization of total pore structure of extra-low permeability sandstone reservoirs: Case study of the Heshui area of Ordos Basin[J]. Petroleum Geology & Experiment,2018,40(4):595-604.

14
SU P H, XIA Z H, QU L C, et al. Fractal characteristics of low-permeability gas sandstones based on a new model for mercury intrusion porosimetry[J]. Journal of Natural Gas Science & Engineering,2018,60:246-255.

15
刘立航,胡海燕,王士跞,等.吐哈盆地沙尔湖地区煤岩孔隙分形特征及其对含气性控制[J].中国科技论文,2022,17(5):495-501.

LIU L H,HU H Y,WANG S L,et al. Fractal characteristics of coal rock pores and their control on gas-bearing propeterties in Shaerhu area, Tuha Basin[J].China Sciencepaper,2022,17(5):495-501.

16
吴浩,刘锐娥,纪友亮,等.致密气储层孔喉分形特征及其与渗流的关系——以鄂尔多斯盆地下石盒子组盒8段为例[J].沉积学报,2017,35(1):151-162.

WU H, LIU R E, JI Y L, et al. Classification of pore structures in typical tight gas reservoir and its significance: A case study of the He 8 member of Upper Palaeozoic Shihezi Formation in Ordos Basin, China[J]. Acta Sedimentologica Sinica,2017,35(1):151-162.

17
张全培,王海红,刘美荣,等.超低渗透储层全孔径分布及其分形特征研究[J].中国矿业大学学报,2020,49(6):1137-1149.

ZHANG Q P,WANG H H, LIU M R, et al. Study of the full pore size distribution and fractal characteristics of ultra-low permeability reservoir[J].Journal of China University of Mining & Technology,2020,49(6):1137-1149.

18
ZHANG L H, LI J C, TANG H M, et al. Fractal pore structure model and multilayer fractal adsorption in shale[J]. Fractals,2014,22(3):90-102.

Outlines

/