基于PEBI网格加密的水平井分段压裂模拟.pdfVIP

大 庆 石 油 学 院 学 报 第 34卷 第 6期 2010年 12月 JOuRNAIOFDAQINGPETR0IEUM INSTITUTE Vo1．34 No．6 Dec． 2010 基于 PEBI网格加密的水平井分段压裂模拟 吴军来 ，刘月 田 (1．中国石油大学 (北京)石油工程教育部重点实验室 ，北京 昌平 102249； 2．中 国石油 大学 (北京)石 油工程学 院 北京 昌平 102249) 摘 要 ：针对有些模拟水力压裂井方法过于简单 ，不能表征压裂井 的流动特性问题 ，通过分析模拟水力压裂井的表 皮因子法 、直角网格加密法和 PEBI网格加密法 ，提 出利用 PEBI网格加密技术建立直井 水平井压裂模型 ，同时结合试 井解释理论 ，推导 出水平井压裂多裂缝 时会出现双径向流 ，验证直井压 裂时存在 明显 的线性流和双线性流动特性．结果 表明：水平井 压裂成单缝 时 ，中部压力下降 比端部快 ，容易水淹 ；水平井压裂形成多缝 时，随着压裂缝 数量 的增 加，缝间距 离减小 ，裂缝 问的相互影响作用明显 ；早期压力下降快 ，由径向流转为双径 向流的时 间变 早 ，试井 曲线上 出现早 期曲线 向 下移动． 关 键 词 ：分段压 裂 ；水平井 ；PEBI网格；双径 向流 ；数值模拟 中图分类 号 ：TE357 文献标识码：A 文章编号 ：1000—1891(201O)06—0053—05 水力压裂作为主要的增产措施在油气 田开发过程 中大量使用 ，压裂可能产生垂直或水平 2类裂缝． 对于一般 的数值模拟软件 ，由于把裂缝和地层看作 同一渗流体系，在 网格划分 中把裂缝单独作为一排 网 格 ，而实际井筒附近形成的裂缝宽度只有 2～5mm，在实际计算中裂缝 内的渗透率和流量会很大，从而会 带来计算收敛性和稳定性问题 ，因此在数值模拟 中如何表征裂缝特性 ，达到模拟效果一直是数值模拟工作 者十分关注的问题． 随着网格技术 的发展 ，局部 网格加密、九点法、十三点法等技术先后被应用于直角网格系统中，而其他 网格技术的发展由三角形 、六边形逐步发展到现在 的Voronoi网格或 PEBI网格．如今 PEBI网格 已有 5 种类型，即六边形 、椭圆形 、长方形 、圆形和其他可变 网格形状 ，这些网格类型先后应用于油藏数值模拟中． 传统的水力压裂模拟方法只是简单的利用等效压裂效果 的思想 ，如等价有效井筒半径 、传导率修正法、井 连接因子 、表皮 因子等方法．但是有的方法 由于其复杂性很难达到模拟大模型水力压裂形成的多裂缝问 题 ，也就不能表征裂缝井的流动特征．为此 ，笔者通过分析几种水力压裂数值模拟方法l1]，提出应用大型 商业数值模拟软件 Eclipse提供的PEBI网格加密功能_2]，对水力压裂缝 ，尤其是水平井多裂缝进行模拟 研究． 1 数值模拟方法 1．1 表皮 因子法 表皮因子法_3是一种最简单的模拟油藏增产的方法，早期 的数值模拟方法模拟油气藏增产都是基于 井筒附近通过增加一个负表皮 因子．该方法的主要缺点是没有考虑裂缝方位的影 响，不能表征裂缝井 的 流动特性． 1．2 直角网格加密法 为了刻画流体在细小 的高渗裂缝 中的流动，可通过在粗网格系统中进行近井附近网格加密后修改其 渗透率来实现模拟水力裂缝的效果．直角网格加密法的优点是可模拟 3个方 向(网格坐标轴方向)，方法 收稿 日期 ：2010—05—25；审稿人 ：张继成 ；编辑 ：关开澄 基金项 目：国家科技重大专项(2oo9zxo5oo9—004—03) 作者简介 ：吴军来(1983一)，男 ，博士生 ，主要从事油气藏渗流力学及数值模拟方面 的研究 大 庆 石 油 学 院 学 报 第 34卷 2010年 相对简单 ；缺点是只能模拟 2个方位的裂缝、网格等宽度 ，不能表征裂缝从射孔到端部的差异效果． 1．3 PEBI网格加密法 PEBI网格 (非结构网格)加密法 是局部网格加密 (LGR)方法的一种 ，与直角 网格加密方法相 比， 该方法可加密成