页岩孔隙结构及多层吸附四分形模型.docVIP

————————————————————————————————————————————————————— 页岩孔隙结构及多层吸附分形模型 分形是1975年由美国学者Mandelbrot[1]首先提出的。自然界中的物体形态各异，结构复杂，组合多样，远远超出了一般意义上研究的规则形状范畴。因此，仅仅采用理想的规则模型研究这些非均质性强、结构差异大的目标有很大的局限性，而这些复杂结构往往表现出 [2][3]分形特征中的幂律关系。Katz等把分形几何理论用来分析多孔介质内部的几何结构。他 们的研究表明；多孔介质的孔隙空间和孔隙界面都具有分形结构，有相同的分形维数，并且可以由分形维数来预测多孔介质的孔隙度。目前在多孔介质孔隙、渗流、吸附等方面已有许多基于分形几何学的研究。 在本章节中，将分别对多孔介质分形孔隙结构模型和具有分形表面的多层吸附分形模型进行研究，在已有模型的基础上进行修正，通过理论分析和实验验证将模型应用于泥页岩的孔隙结构和吸附特性研究上，分析分形维度对泥页岩多孔介质各种物性参数的影响。 多孔介质孔隙结构模型 Menger 海绵模型是应用最为广泛的多孔介质分形模型，Menger海绵模型是在Sierpinski 方毯的基础上在三维空间中的扩展[4]。Menger海绵模型能够对许多多孔介质进行有效的表 [5]征。Jin Yi改变了Menger 海绵模型的构造过程，构造出了具有连通结构的“SmVq”孔隙 模型，同时给出了模型分形维度的计算公式： 332logNlog(m+2q?3mq) =Dlogmlogm(1) 其中，D是分形维数；N代表剩余的小立方体个数；m是每边分割的分数。 采用该方法构造孔隙结构模型：1、将边长为R的正方体分成 方体边长为个小立方体，每个小立，沿贯穿每个面中心的相互垂直轴线挖去q个小立方体；2、在得到的小立方体基础上，重复步骤1 。 图1 两次迭代后的SmVq模型截面图 Hunt[6]指出，多孔介质多为固体介质和孔隙两相组成。如果多孔介质具有分形特征那么要么是孔隙分形要么是固相介质分形。在分形模型建立的过程中，一般对固相介质进行分形 1 描述，其思路为：在每次迭代过程中，模型由相同大小的颗粒组成而孔隙尺寸则不相同，此时固相介质分布呈现分形特征，尽管孔隙在几何表现上不是分形模型，但是其数量~尺寸分布却呈现出幂律指数关系并且分形维度和固相介质相同，所以用一个分形维度可以同时表示固相颗粒和孔隙分布的分形结构，尽管他们表述的途径不同。因此，式(1)中的D值可以表示孔隙分布的分形维度。 对于孔隙分布具有分形特征的多孔介质，其大于某一孔径的孔隙数量N与孔径r之间 [7]遵从以下关系： =N(≥r)∫rmax r=f(x)drar?D (2) 其中 a 是相关系数，f(x)为孔径分布密度函数，可表示为： 孔隙累计体积V(≤r)： dN(≥r)=?Dar?D?1 f(r)=dr3V(≤r)=frβrdr ()∫rminr(3) (4) 式(4)中，β 是与孔隙形状有关的因素。Kat[3]提出了基于分形维度计算孔隙度φ 的方法, Yu B[8]给出了更一般的形式： ?r?φ=?min??rmax?De?D (5) 其中 De 是几何空间的分形维数，三维空间 De=3，rmin 、 rmax 表示孔径分布区间。 通常，在孔径分布区间内rmin?rmax。某些情况下，孔径分布呈现出多分维的现象，此时孔隙度可表示为： ?r?φ=φ1+φ2=?min1??rmax1?3?D1?r?+?min2??rmax2?3?D2 (6) 式(6)假设在孔径分布范围内有两个分形维度：D1 和 D2，在每段分布区间内，孔隙度均可以通过(5)式分别算出。 根据式(4)，可以得到孔隙在分形几何分布下的孔径分布与累积孔隙体积关系式，对其求导可以得到： 对于SmVq模型，通过对比不同 度随着值的增大而减小，这是因为dV=?Daβr2?D dr(7) 值对分形维度的影响（图2A），可以看到，分形维值越大，模型越接近完全孔隙化，相应的固相介 值，随着m值的质减少，其复杂程度随之减小，所以分形维度变小；同时，对于固定的 增加（q值同样增加）分形维度增加，这是因为m值的增大相当于测量精度的增加，这与盒维数计算中度量尺寸的选取道理类似。 2 图2B显示的是对于同一m值，随着q值的增加分形维度的变化情况，可以看到，随着q值增加分形维度减小，同时m值变化引起的分形维度变化趋势与图2A显示的结果相同。