三角网格差分方法在油藏数值模拟中的应用.pdfVIP

维普资讯 破， ； 南叼鹤荔域 54 1999年 l2月 三角网格差分方法在油藏数值模拟 中的应用 @ 徐怀大 孔祥言 了 ?7 一 (1)中国科学技术大学 (2)中国地质大学(北京) 丁E／7 ‘ 量 ，有 引 言 N Q (f)=∑ (f) (=o，g，w) (1) = 1 随着地震勘探分辨率的不断提高，油藏数值模拟 中建立的储集层地质模型更精细、结构更复杂，利用传 统的矩形网格剖分渗流区域 已不够灵活 (特别在边界 和断层处理上)。矩形网格可 以采用不等长间距来加 密网格，但要增加很多结点，使计算工作量增大，而且 “ 、 (a) 、 只能用折线边界来近似弯曲的渗流区域边界，对断层 圈 1 均衡单元(a)和侧 向谵■横断面(b】示意圈 一 般处理为一个特殊的渗流带。三角网格的边界本身 由图1b可看出，当面元很小时，可以认为通过断 就比较接近弯曲边界，还能方便地处理渗流区内的断 面 的侧向流量是相同的。按图lb所取坐标 层 。三角网格剖分渗流区域应遵守的原则只有两项： 系及均衡单元与流量断面的位置，根据达西定律有 ①某些面元(即三角形单元)的顶点(称为点元或结点) 不允许落在线元 (即三角形的3条边)的非端点位置 )= ㈤ 一 ddS(t) ^ 上，②面元的任一内角不应大于90。或小于30。，3条边 (2) 的长度应尽可能接近_1J。此外，为了更好地验证模型， 尽量使结点与采液井和注液井重合。因此，在不违反 为了求压力梯度J鲁(f)和蹋 (t)，设在面元卢上 以上原则的前提下，三角网格可 以方便灵活地加密重 流体压力 (，Y，t)为坐标的线性函数，即 点研究区结点。本文提出的三角网格差分方法在油藏 数值模拟中应用效果很好 ，模型计算也很方便。 (，Y，t)=。+ +cy (3) 对于渗流区域 (储集层、油藏)来讲，P( ，Y，t)是 三角网格差分方程的形成 一 个非常复杂的连续函数，几乎不可能用解析式来表 示。但当渗流区域被三角网格剖分后，(3)式的假设是 描述油藏流体运动规律的偏微分方程主要 由运动 可以接受的，只要面元取得足够小，而且满足在点元上 方程、连续方程和状态方程组成。从均衡原理来讲 ，一 等于点元的真实流体压力的条件，即 个典型单元的微分方程包括单元 的侧向流量、源汇项 及其内部储存量的变化量。本文以二维三相黑油模型 (z，Y，t)=Pf(f) (1= i，j，k) (4)