分形低渗透缝洞型油藏非线性渗流规律-计算物理.PDF

第３５卷 第１期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３５，Ｎｏ．１ ２０１８年１月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｊａｎ． ，２０１８ 文章编号：１００１⁃２４６Ｘ（２０１８）０１⁃００５５⁃０９ 分形低渗透缝洞型油藏非线性渗流规律 刘化普，　 刘慧卿，　 王　 敬 （中国石油大学石油工程学院教育部重点实验室，北京　 １０２２４９） 摘　 要：引入分形理论，建立考虑低速非达西效应的分形三重介质缝洞型油藏数学模型，通过 Ｌａｐｌａｃｅ 变换及 Ｓｔｅｈｆｅｓｔ数值反演方法求出井底压力，借助Ｍａｔｌａｂ编程绘制压力动态曲线，划分渗流阶段，分析渗流规律，进行非线 性参数敏感性分析．最后结合实际算例，验证模型的正确性．结果表明：分形三重介质油藏渗流过程分为早期纯井 储，过渡流，缝洞窜流，拟径向流，基质与溶洞、裂缝窜流及总体径向流６个渗流阶段；低速非达西效应对渗流的影 响随时间的推移逐渐增大；启动压力梯度越大，总径向流阶段压力动态曲线上翘幅度越大；分形系数影响整个渗流 过程，随着分形系数的增大，裂缝迂曲程度随之增大，致使渗流阻力增加，引起压力动态曲线整体上移幅度增大． 关键词：缝洞型油藏；分形；三重介质；非达西；启动压力梯度 中图分类号：ＴＥ３４４ 文献标志码：Ａ ＤＯＩ：１０．１９５９６／ ｊ．ｃｎｋｉ．１００１⁃２４６ｘ．７５８９ ０　 引言 低渗透缝洞型油藏低孔、低渗、裂缝、溶洞结构发育，形成基质、溶洞、裂缝三重介质孔隙结构系统［１－４］． 由于其储层低孔、低渗的特性，地层流体在油藏中的渗流不再遵循达西定律，而是呈现出低速非达西效应，即 前人引入的启动压力梯度的概念［５－７］．另外，裂缝与溶洞结构的存在导致油藏具有高度的非均质性及复杂的 储渗空间．非达西因子的影响及孔隙结构的复杂性致使对其渗流规律研究比较困难．为了探究低渗透缝洞型 三重介质油藏的渗流机理，国内外很多学者做了大量的研究．针对缝洞型油藏孔隙结构的复杂性建立了三重 介质单渗与双渗渗流模型，对其拟稳态及非稳态渗流特征进行了分析［８－１３］．针对油藏低速非达西效应，引入 启动压力梯度的概念，建立非线性渗流模型，研究了考虑启动压力梯度下的低渗油藏渗流规律，且取得了较 大的进展［１４－１７］．以上均是在欧式几何的基础上对缝洞型油藏渗流机理进行的研究，而实际地层中裂缝分布 较为复杂，欧式几何空间难以描述裂缝实际分布情况．然而分形理论对不规则，不均匀的复杂系统恰好能比 较精确地描述． 自１９７９年Ｍａｎｄｅｌｂｒｏｔ［１８－１９］首次提出分形概念以来，分形理论便被运用到描述复杂物体当中， 且取得较大进展．后续很多学者将分形理论引入到渗流力学当中，建立分形渗流模型，分析渗流特征［２０－２５］． 然而目前对于分形渗流的研究，多是建立在均质及双重介质模型之上，对于分形三重介质模型的研究较少， 又由于考虑低速非达西效应，使得对其分形渗流规律的研究变得更加困难．本文在此基础之上，引入分形理 论，考虑裂缝系统分形特征及油藏低速非达西效应，建立基质、裂缝、溶洞分形三重介质模型，研究低渗透缝 洞型油藏的渗流规律．首先建立分形二叉树模型，对裂缝系统孔隙度及渗透率进行分形表征，获得裂缝系统 分形孔隙度及分形表观渗透率的表达式；在此基础之上，建立考虑启动压力梯度的分形三重介质缝洞型油藏 数学模型，通过Ｌａｐｌａｃｅ变换与Ｓｔｅｈｆｅｓｔ数值反演方法获得井底压力，借助Ｍａｔｌａｂ编程，绘制分形三重介质低 渗透缝洞型油藏的渗流特征曲线，划分渗流阶段，分析渗流规律． １　 分形非达西模型数学描述 １１　 裂缝系统分形描述 分形通常被定义为“一个粗糙或零碎的几何形状，可以分成数个部分，且每一部分都（至少近似地）是整 收稿日期：２０１６－１１－２１；修回日期：２０１７ －０３－１２ 基金项目：国家科技重大专项（２０１６ＺＸ０５０１４００３００４）资助项目 作者简介：刘化普（１９９３－），男，在读硕士生，研究方向为油藏工程及提高采收率技术，Ｅ⁃ｍａｉｌ：ｌｉｕｈｕａｐｕ＠ｆｏｘｍａｉｌ．ｃｏｍ ５６ 计　 　 算　 　 物　 　 理