页岩储层分形特征及孔隙结构研究

页岩储层分形特征及孔隙结构研究

对于独立储层，天然气主要以两种方式收集：（1）以自由态收集在微裂缝或微裂缝中；（2）它以吸附形式吸附在干酪根和粘土的表面。众多学者在研究中总结了页岩气的岩石物理特性,发现页岩储层的产能受其孔隙形态的影响。页岩气储层有多种孔隙类别,孔隙结构十分复杂,用霍多特(1966)的分类方法,可以把页岩气储层的孔隙分成:大孔(孔径>1000nm,主要的渗流空间)、中孔(100nm<孔径≤1000nm,层流渗流空间)、过渡孔(10nm<孔径≤100nm,气体扩散空间)、微孔(孔径≤10nm,是主要的吸咐空间)。当孔隙半径>100nm时,自由气能够在孔隙中渗流,其渗流会影响页岩气储层的渗流效应。而吸附气主要存在于半径<100nm的孔隙中,影响着页岩气的聚集和扩散。分形是指局部与整体在形态、功能和信息方面具有相似性的集合,这种自相似性仅在一定的尺度范围、一定的层次中才表现出来。分形维数可以用来定量地描述孔隙结构的复杂程度与非均质性,其测定的方法有很多,包括压汞法、图像分析法和气体吸附法等。前人结合测井和电镜扫描,已经将分形理论广泛应用于微孔隙结构分析的研究中,特别是将分形理论应用于常规储层的孔隙分布描述中。岩石孔隙的分形维数通常在2~3之间,分形维数越接近2,说明孔隙表面越光滑,储层的储集性能越好;反之,分形维数值越接近3,储层孔隙结构越复杂,储集性能越差。Katz和Krohn等通过扫描电镜观测孔隙结构,发现油藏中孔径在0.2~50μm范围内具有分形特征,分形维数为2.27~2.89。Mahamud等结合压汞实验、核磁共振、CT扫描和X射线衍射,利用分形理论研究了煤层气储层的孔隙特征。但是,结合压汞曲线、利用分形方法系统研究页岩气储层特征的文献还不多。1非润湿相汞压入岩心的体积页岩气、煤层气等非常规储层都具有复杂的孔隙结构,用于煤层气的毛管压力曲线的分形几何表达式同样适用于页岩气储层。非润湿相(汞)压入岩心的体积为式中:Vmp为非润湿相(汞)压入岩心的体积;Vt为连通孔隙总体积;pc为毛管压力值;pcmin为最小毛管压力值;D为分形维数。根据推导,非润湿相(汞)压入体积与毛管压力的关系也可以按下式计算根据岩心的毛细管曲线,可以采用(2)式求出分形维数D和pcmin。2压汞曲线特征研究样品取自四川盆地上三叠统页岩气储层。通过岩石样品的X射线衍射实验和扫描电镜测试,得到了样品的物性参数:孔隙度为2.34%~3.56%,平均孔隙度为2.88%;渗透率范围为(0.13~0.06)×10-3μm2,平均渗透率为0.036×10-3μm2。样品中石英的质量分数>50%,方解石的质量分数<10%。压汞实验研究采用的仪器是9500自动压汞仪,实验所得压汞曲线可以总结为2类,实验结果见图1。以2号岩心样品的压汞曲线为例,具体分析第一类压汞曲线的特征(图2)。压汞曲线可分成2段(Ⅰ段和Ⅱ段)。其中,第Ⅰ阶段又可进一步划分为a和b两部分。当压力>0.04MPa时,汞开始进入岩石,压汞曲线较平坦,这一阶段为a部分;当压力>0.15MPa,压力迅速增加,但汞的注入量却很少,这一阶段为b部分。第Ⅰ阶段的曲线形态与常规岩心的压汞曲线形态类似,表明汞注入岩石后,初期阶段压入微裂缝中,随着压力的升高,汞进入较小的孔隙中。当汞的驱替压力约>10MPa,压汞曲线进入第Ⅱ阶段后,压汞曲线呈现反常特征:压汞曲线为凸形,汞的注入量迅速增加,但压力增加不大。以4号岩心为例,分析第二类压汞曲线的特征(图1-B)。这类压汞曲线与第一类压汞曲线相比较,缺少第Ⅰ阶段的a部分。压汞曲线同样可以分为2个阶段,具体表现为:初期的进汞压力迅速增加,但汞的注入量却很少,表明岩心中缺少大孔和微裂缝;与第一类压汞曲线类似,当进汞压力约>10MPa时,压汞曲线进入第Ⅱ阶段后同样呈现反常特征:汞的注入量迅速增加,但压力增加不大。为了核实四川盆地上三叠统页岩气储层的压汞曲线特征是否具有代表性,我们与国外页岩储层的压汞曲线特征进行了对比,发现Barnett储层的Utica页岩的压汞曲线与本文提到的第二类曲线相似(图3),Elgmati研究的注入曲线与第一类曲线类似(图4)。3储层参数特征根据公式(2)分别计算出2号样品和4号样品的,研究它与lgpc之间的关系。图5是和lgpc的关系曲线。图5-A中显示了3段不同斜率的直线段,分别对应于第一类压汞曲线中第Ⅰ阶段的a部分、b部分以及第Ⅱ阶段。图5-B中的2段不同斜率的直线段对应于第Ⅱ类型压汞曲线的第Ⅰ阶段和第Ⅱ阶段。由公式(2)可知,直线段的斜率为D-4,通过斜率可以计算出分形维数。采用压汞曲线计算的分形维应该在2~3之间。越接近2,说明孔隙表面越光滑,储层的储集性能越好;反之,分形维数值越接近3,储层孔隙结构越复杂,储集性能越差。将图5-A和图2进行对照。图5-A中,只有当lgpc为-0.8~1(即pc约为0.15~10MPa)时,曲线的斜率为-1.5004,对应的分形维数约为2.5。这段曲线对应于图2中第Ⅰ段的b部分;图5-A中,曲线的斜率为-0.4909和-0.0236时,对应的分形维数>3,这2段对应于图2中的a部分和第Ⅱ阶段。图5-A证实了图2中a部分的压汞曲线不反映孔隙的结构特征,而是反映了微裂缝的存在:即使压力很小,汞也能够进入岩心中的微裂缝。第Ⅱ阶段出现的原因是由于压汞过程中,由于压力升高,岩石被压坏,则方程不再适用。所以在研究孔隙分布时,分形维数>3的压汞曲线不能反映页岩孔隙的特征。图5-B中,当lgpc<1(即pc<10MPa)时,曲线的斜率为-1.062,对应的分形维数约为2.938。与第一类压汞曲线相比较,储层孔隙结构复杂,储集性能较差。由油层物理学可知,毛管孔径与毛管压力之间存在以下关系式中:σ为液体界面张力;θ为润湿接触角;r为孔隙半径。其中,汞的σ=4.8×10.3N/cm2,θ=140°,则由公式(3)可将压汞曲线中毛管压力值换算成孔径。根据图1中第Ⅰ阶段和第Ⅱ阶段的分界点,发现当进汞压力约>10MPa后,各样品的压汞曲线出现异常特征。10MPa对应的孔隙半径为0.1407μm,说明压汞曲线能反映的孔隙半径的下限为0.1407μm。这也从一个方面证实了田华等人的观点,即:压汞方法能反映宏孔(孔径>50nm)的分布特征。根据前人文献中的压汞曲线(如图3、图4中的压汞曲线),也可得出类似的结论。4压汞曲线的发育时期a.采用压汞曲线可以判断页岩的微观孔缝结构类别。页岩气储层压汞曲线有2种类型:一种是有微裂缝存在的压汞曲线;另一种无微缝发育,只发育微孔和中孔。具有第一类压汞曲线的页岩,其孔隙和裂缝发育好,具有较好的开采潜力