各向异性渗透率油藏数值模拟.pdfVIP

[文章编号]1001．246x(2007)03．0295．06 各向异性渗透率油藏数值模拟 刘月田1， 徐明旺2， 彭道贵3， 张宏伟4， 董 彬4 (1．中国石油大学石油工程教育部重点实验室，北京102249；2．辽河石油勘探开发研究院，辽宁盘锦124010； 3．胜利油田地质科学研究院，山东东营257060；4．冀东石油勘探开发研究院，河北唐山063004) [摘 要】 建立三维三相各向异性油藏渗流数值模拟模型并应用于油藏开发．模型包含完全各向异性渗透率张 量。采用十九点有限差分格式和全隐式迭代方法求解．可以描述不同区域不同主方向的各向异性渗透率分布， 准确计算各向异性油藏内流体运动过程、剩余油分布和压力分布． [关■词】 各向异性；油藏；数值模拟；剩余油；开发调整 [中圈分类号] 佃19 [文献标识码】 A 0 引言 目前，我国石油工业已经全面进入复杂油气藏开发阶段，其中各向异性油藏的开发问题日益突出．各向 异性油藏就是渗透率具有方向性的油藏，既包括低渗砂岩油藏、碳酸盐岩油藏等裂缝性油藏，也包括注水开 发中油水流动具有明显方向性的特殊沉积相油藏．各向异性油藏在石油储量中占有很大比重．现有的油藏数 值模拟软件无法准确模拟一般各向异性油藏内流体的渗流规律，在预测水的突破时间、压力分布及剩余油分 布时存在误差，从而在指导油藏工程措施如前期开发方案的部署、中后期开发方案的调整时有可能出现重大 失误，造成严重的经济损失．为解决这一棘手问题，首先要研究各向异性介质中的渗流规律，得到正确描述油 气水渗流运动的数学模型，然后选择合理的算法，建立数值模型及其求解过程，最后开发实用的各向异性油 藏数模软件，用于各向异性油藏开发实践． 1各向异性油藏渗流数学模型 各向异性油藏中，渗流速度的方向由位势梯度的方向确定，但两者通常不平行，当且仅当位势梯度方向 在渗透率主轴上时，渗流速度的方向才和位势梯度方向平行．渗流速度的大小由位势梯度的大小和方向共同 决定，即无论位势梯度的大小还是方向发生变化，渗流速度的大小都会改变．这是各向异性油藏渗流主要的 基本特征，它存在于各向异性油藏渗流与开发的各种复杂问题之中．解决的关键就在于正确认识、掌握和利 用这一特征…． 一般数值模拟模型和软件在处理各向异性油藏时，假设油藏内各区域的渗透率主方向都相同，即认为油 藏各区域的裂缝方向或沉积方向皆一致，所以只是“简化各向异性模型”[2“1．本文研究建立三维三相油藏渗 流”完全各向异性模型”． 1．1基本假设条件 基本假设条件如下： ①油藏渗透率为完全各向异性； ②油藏中最多只有油、气、水三相，每一相的渗流均遵循广义达西定律； ③油藏中的烃类只含有油、气两个组分； ④气体溢出或溶解于油的过程瞬时完成，即认为油藏中油、气两相始终平衡； 【牧稿日期]2006一01—26；【謦回日期】2006—0r7一17 [基盒项目]国家自然科学基金(批准资助项目 [作者简介]刘月田(1965一)，男，河北无极，教授，博士，从事油藏管理、油藏数值模拟及油田开发设计研究，北京市昌平区府学路18号 万方数据万方数据 296 计 算 物 理 第24卷 ⑨油、水之间不互浴，天然气不溶于水； ⑥油藏中的渗流是等温渗流． 1．2建立微分方程组 在直角坐标系(算，y，z)中油相、气相、水相渗流速度分别为 D]， P。一(y。+79d)V K=一字置·[V D]， (1) P。一)，。V 匕=一争K·[V P。一7。v” “ D]， y。：一孚置．[v