流体流经自然岩石裂缝的数值模拟.pdfVIP

维普资讯 ■ 油气国开发 ■低渗透油气田20O6年6月第11卷第1—2期 流体流经 自然岩石裂缝的数值模拟 贺成才 (西南石油学院) 摘 要 避开 自然裂缝面的复杂数学描述 ，直接利用正态分布数值生成 自然岩石裂缝 的宽度 ．然后利用Reynolds方 程对流体流经狭缝的描述 ，采用先进 的线性化方法与数值迭代技术，对变系数的 Laplace方程求解 ，最后与惯常的平 板模型计算的流量结果进行 了对比，二者相差甚大。本文的研 究为裂缝性油气藏开发 中矿场设计水力参数提供 了科 学依据。 关键词 自然裂缝 裂缝宽度 数值模拟 牛顿流体 随着常规油气藏资源的逐渐枯竭 ，非常规油气 缝宽度h的估计 ，以此作为屏蔽暂堵技术实施时的 藏的开发越来越重要 ，然而非常规油气藏大多为基 依据等。 质渗透率很低的裂缝性油气藏 。所以研究裂缝性油 本文利用经典 的润滑理论 中计算狭缝水力参数 气藏中 自然裂缝 的导流能力十分重要 。在经典的教 的雷诺方程及现代先进数值计算方法 ，研究 自然岩 科书中，描述流体在 自然岩石裂缝中流动几乎都采 石裂缝 的导流能力 ，在 500x500的大范围网格上进 用 的所谓平板模型假设 ，即把 自然裂缝假设为两个 行计算 。并与经典的平板模型计算结果进行 比较 ，文 相互平行 的有一定距离的光滑平面。如图 1所示，对 献[1]、[4】、[6】、[7】、[81同样研究了这一问题 ，由于使用 于流体沿压力差方向流经 自然岩石裂缝流动 ，利用 的传统数值计算方法 ，因此在计算机上实现有很大 的局限性 ，对计算机的内存与运算速度均有较高的 要求 ，它们仅能在 65x65的网格上进行数值模拟 ，而 在如此小的网格上进行数值计算 ，可能在边界效应 影响很大的区域内进行数值模拟 。其结果可能与真 实情况相差甚远；另外 ，最重要的缺点是使用传统的 数值计算方法可能使计算过程中途中断 ，无功而返。 1裂缝宽度数据的数值生成 众所周知 ，一个 自然岩石裂缝是由两个十分粗 糙的岩石面叠合而成。显然要形成一个 自然岩石裂 稳定层流的 Navier—Stokes方程可解得 ： 缝就必须分两步：首先数值生成两个裂缝壁面 。然后 两个裂缝壁面按一定的缝宽迭合在一起。文献[21、[51 Q：茜一 (1)均是按此思路数值生成 自然岩石裂缝的。但是 ，这样 式中：()_一 单位时间内的体积流量 ，m3／s； 做存在很大缺点。首先。无论是利用普通的随机生成 — — 流体的动力粘度，Pa·S； 还是利用分形的傅里叶滤波法生成方法 。裂缝壁面 压力 ，Pa；