低渗透储层岩石孔隙结构的分形模型.docxVIP

低渗透储层岩石孔隙结构的分形模型 1 影响收率的因素 储层岩石的空隙结构不仅影响石油的沉积，而且影响油气井的产量和最终采收率。储层岩石的孔隙结构是一种分形结构,分形几何的出现为孔隙结构的研究提供了一种新的方法,可以用分形维数来定量描述,它不仅能描述孔隙结构的孔隙大小分布和复杂程度,而且能描述储集层岩石的粒度组成等。 2 分形维数法。对于完备的、断裂 分形几何是1975年由法国著名的数学家曼德勃罗特(Mandelbrot)从拉丁文fractus(破碎的、断裂的)创造的新词,意思是破碎、细片、分数、分级等等。它研究的是自然界中不稳定的、非常不规则的、但是具有自相似性、自反演性的不规则现象,可以用分形维数来表示,即 或者写为 式中:ε为测度,N(ε)为在测度下所得到的度量值,D为研究对象的分形维数。在双对数坐标下,N(ε)-ε为一直线,其斜率就是研究对象的分形维数。 3 储层岩石中主孔隙结构模型 由于储层岩石的孔隙结构具有分形特征,根据分形几何原理,储层中孔隙半径大于r的孔隙数目N(r)与半径r之间有如下的关系: 式中:rmax为最大孔隙半径,f(r)为孔径分布密度函数,D为孔隙分维数,α为比例系数。将式(2)对r求导,可得到孔径分布的密度函数f(r)的表达式: 将式(3)代入储层中孔隙半径小于r的孔隙累积体积V(r)的表达式,并积分可得 式中:rmin为最小孔隙半径,β为与孔隙结构有关的常数(如孔隙为立方体,β=1;孔隙为球体,β=4π/3)。同理,可得到储层的总孔隙累积体积V为 根据V(r)与V,可得到孔隙半径小于r的累积孔隙体积分数的表达式为 由于rminrmax,则上式可以简化为 根据Laplace方程: 式中:Pc为孔隙半径r所对应的毛管压力;σ,θ分别为液体的表面张力和润湿接触角。将式(8)代入式(7),得 式中:Pmin为最大孔隙半径rmax所对应的毛细管压力,即入口毛细管压力;S为毛细管压力为Pc时储层岩石中润湿相的含水饱和度。所以,式(9)即为毛细管压力曲线分形几何公式。 4 孔隙概率密度分布 用式(3)可计算孔隙分布的概率密度函数,但在该表达式中有一个比例常数α很难确定,为此,采用古典型概率分布的定义方法来计算孔隙的概率密度分布。 设孔隙半径r的范围为[rmin,rmax],按古典型概率分布的定义,有 则分布函数为 将式(2)代入式(10)中并整理可以得到孔隙分布函数为 对式(11)求导数可以得到孔隙分布的概率密度函数表达式为 5 计算分形维数d和最小毛细管压力amn 分形维数的计算可以采用图版法和拟合法。图版法的基本原理是在双对数坐标系中作S和pc的关系曲线,可得到S和pc在双对数坐标系中为一条直线关系,根据直线的斜率和截距可以计算出分形维数D和最小毛细管压力Pmin。这种图版法简便快捷,但是精度较低,自动化程度低。为了提高精度,这里利用线性回归拟合法来计算分形维数D和最小毛细管压力Pmin。 对式(9)两边分别取对数,可得 设lgS=y,lgpc=x,3-D=a,lgpmin=b,则式(13)变为y=ab-ax,通过线性回归,可得到精确的D和Pmin。 6 储层岩石分形 在三维欧氏空间中孔隙结构的分形维数是介于2和3之间的分数。分形维数越接近2,即分形维数越小,说明孔喉表面越光滑,均质性较强,岩石的储集性能越好;反之,分形维数越接近3,即分形维数越大,孔喉表面越不光滑,岩石的储集性能越差,表明该分形结构的复杂程度越大,非均质性越强。根据这一性质,对不同低渗透岩石的岩样进行了研究。表1为根据不同低渗透岩样的毛细管压力资料计算出来的储层岩石的孔隙分形维数和孔隙结构参数。 从表1中可以看出,孔隙结构分形维数都接近3,均在2.7以上。说明孔隙结构的复杂程度较大,非均质性也很严重。回归方程的相关系数大部分都在-0.99以上,相关性比较好,这说明用毛细管压力资料来计算孔隙结构的分形维数是可行的。 图1,2分别为2,9号岩样毛管压力pc和含水饱和度Sw的双对数曲线。很明显,在双对数坐标中,pc和Sw呈很好的直线关系,拟合程度比较好,说明孔隙结构具有较好的分形特征。 7 储层岩石孔隙大小分布 根据式(13),通过线性回归,可以计算出孔隙结构的分形维数,再根据式(12)就可以计算出储层岩石孔隙大小的概率密度分布,由此,可以计算出孔喉概率分布。图3,4分别为2,9号岩样的孔隙大小分布曲线。 从图3,4中可以看出,岩样中孔隙主要集中在孔隙半径较小的区域,而孔隙半径较大的孔隙出现的概率比较低,呈中间大,两头小的分布。 8 孔隙大小分布的测量 (1)建立了毛细管压力曲线和孔隙大小分布的分形几何模型,在此基础上计算了孔隙结构的分形维数。 (2)利用孔隙结构的分形维数既可以定量描述孔隙结构参数,也可以定量描述孔隙结构的复杂程度。 (3)