2分形维数.doc

文章编号：CSTAM2014-B01-0168 四川龙王庙组缝洞型白云岩气藏储层孔隙结构分形研究 吴国铭1，刘华勋1 （1 中国科学院渗流流体力学研究所，河北廊坊 065007） 摘要： 白云岩油气藏储层缝洞孔隙发育，非均质性强，储层特性表征方法不完备。获取四川磨溪区块龙王庙组白云岩磨溪12井23块岩心压汞实验数据，绘制进汞饱和度与毛管力双对数曲线，利用分形几何理论计算得到不同孔隙结构的分形维数，可表征储层非均质性。缝洞系统分维数较大，分布范围较广，基质分维数与砂岩相似，白云岩储层非均质性强于砂岩。 关键词： 分形，白云岩，缝洞孔隙结构 1.引言 白云岩油气藏储层在化学成因主导下具有矿物成分简单，结构复杂的特征，主要表现在储层缝、洞充分发育，非均质性强。目前，还没有一套完整的方法来表征碳酸盐岩储层的孔、洞和缝，从而准确描述整个油气藏的储层特性[1]-[2]。分形几何学是研究自相似结构的有效方法，许多学者的研究结果表明碳酸盐岩储层裂缝及孔洞等孔隙结构具有自相似性，可用分维数辨识区分不同孔隙结构[3]-[5]。四川磨溪区块龙王庙组气藏储层为裂缝~孔、洞型碳酸盐岩储层，存在溶蚀孔洞型、溶蚀孔隙型和基质孔隙型三种储集类型，储层非均质性极强。为辨识表征不同孔隙结构，利用龙王庙组磨溪12井23组压汞数据计算不同尺度孔隙结构的分形维数，研究储层特性。 2.分形维数 压汞曲线是不同尺度的孔隙结构毛管力的直观表征，不同孔吼半径对应曲线上不同的进汞压力。岩心含汞饱和度定义为： （1） 式中为孔隙体积，为累计进汞体积。 将储层近似看作毛细管模型，则有： （2） 式中为不同毛细管半径，为毛细管长度，为不同半径毛细管个数。进汞累计体积毛细管压力的函数，即=。而毛管压力满足： （3） 式中为表面张力，为界面夹角，为毛管力。岩心孔隙体积为定值，由式（1）、（2）和（3）可知，进汞饱和度是毛管半径和毛细管个数的函数，即=。 分形理论中的计盒维数法[6]有： （4） 式中为半径大于r的累积孔隙体积，为最大孔隙半径，为孔径分布密度函数。根据欧式空间和分维数的关系，可知是的函数，其中D为分维数，即=。综合分析得到=。由式（3）可知，最后得到了进汞饱和度与毛管力的关系： = （5） 式中为常数。 对式（5）两边取对数，可知lg与lg为线性关系，斜率k=D-2，分维数D=k+2。由此，可通过压汞数据求取不同尺度孔隙结构的分维数，研究储层特性。 3. 结果分析 选取龙王庙组白云岩气藏磨溪12井23组压汞数据中的4组代表性数据按照式（5）在双对数坐标轴上绘制岩心进汞饱和度与进汞毛管压力的关系曲线（图1）。从图1可看出，缝洞较发育的白云岩储层进汞饱和度与毛管力双对数曲线呈两段式，由两条斜率不同的直线组成，在4号岩心的进汞饱和度与毛管力双对数曲线（图2）中尤其明显。文慧俭等认为斜率不同的两条直线分别与尺度不同的孔隙结构相关，斜率大的直线表征大孔隙，斜率小的直线表征小孔隙（图3）。 图1磨溪12井23块岩心进汞饱和度与毛管力双对数曲线 图2 4号岩心进汞饱和度与毛管力双对数曲线 图3 17号岩心分段拟合大、小孔隙分维数 分别求取23块岩心双对数曲线中两条直线对应的不同尺度孔隙结构的分形维数，得到表1所示的结果。由表1可看出，缝洞发育的白云岩储层小孔隙结构分维数在2.0～2.5之间，大孔隙结构分维数分布范围较大，在2.2～5.5之间，表明裂缝和孔洞较发育，与其它的岩心物性测试结果相符。相比于裂缝和孔洞的分维数，基质的分维数普遍较小，且分布较为集中。利用分维数可辨识储层不同尺度的孔隙结构。 表1 磨溪21井23块岩心分维数计算结果 岩心号 岩性 孔隙度% 渗透率mD 大孔径分维数 小孔径分维数 1 白云岩 6.8 0.00375 2.86 2.04 2 溶孔洞白云岩 4.1 0.00225 2.97 2.34 3 溶孔洞白云岩 10.4 15.20000 2.60 2.07 4 溶孔洞白云岩 4.6 0.00264 3.21 2.26 5 白云岩 3.7 0.00635 2.85 2.17