非均质低渗透油藏非线性流动数值模拟.pdfVIP

第 ??卷 第 ?期? 计? 算? 物? 理? ???．??，????? ????年 ?月? ????????????????????????????????????????? ???．，????? 文章 编 号 ：????—????????????—???????? 非 均质 低 渗 透 油 藏 非线 性 流 动 数值 模 拟? 姚? 军 ，? 黄? 涛 ，? 黄 朝 琴? ?中 国石 油 大 学 ?华 东 ?石 油 工程 学 院 ，山 东 青 岛 ???????? 摘 要 ：提 出 一 种适 用 于 非 均 质 低 渗 透 油 藏 非 线 性 流 动 数 值 模 拟 的 方 法 ．通 过 引 入 单 元 边 界 压 力 和 速 度 两 组 物 理? 量 ，在单个网格 单元 中构造数值计算格式 ，可精确地 获取 每个 网格单元 的压 力梯 度．以单 相非线性渗 流问题 为出发? 点 ，详 细 阐述 该方 法 的 基 本 原 理 ；通 过 数 值 算 例 验证 方法 的 正 确性 和适 用性 ．? 关键 词 ：非均质 ；低渗透油藏 ；非线性渗流 ；油藏数值模拟? 中 图分 类 号 ：?????；?????；????．??? 文 献标 志 码 ：?? ?? 引 言? 随着我 国油 气 田勘探 开发 的不 断深 入 ，低 渗透 油 藏 已成 为我 国增 储 上 产 的 主要 领 域 ? ．大 量低 渗 透渗? 流实 验和生 产 实践表 明? ：油气流 动存 在最 小启 动压 力 梯 度 ，其 渗 流 曲线 偏 离 了经 典 的 达西 定 律 ，具 有 显? 著 的非线性 特征 ．许 多专 家学 者对其 流 动机理进 行 了大 量 研究 ，提 出 了多 种非 线 性 渗流 模 式 ，包括 拟 启 动 压? 力梯 度模 型 、三 阶段模 型 、非线 性连 续性 模 型 等?。?，并 进 行 了 相应 的数 值 模 拟研 究? ．现 有 的低 渗 透 油? 藏非线 性流 动数 值模 拟 ，大都是 采用传 统 的有 限差 分法 或控 制体积 有 限差分 ．此类 方法 均基 于相邻 单元 中心? 点 的物理量 值来 构造 数值 计算 格式 ，求 取驱替 压力 梯度 ，进 而结合 渗流 曲线 来进行 判 断和求 解 ．然而 ，对 于非? 均质 低渗 透油藏 ，不 同 区域甚 至不 同 网格 单元 的渗 流模 式不 尽相 同 ，数 值模 拟 中必须 获取各 网格 单元 的压 力? 梯度 ，并选 择相 应 的渗流模 式 ，此时 上述 传统 的数值 计算 方 法 已无 法满 足要 求．? 对此 ，在控 制体 积有 限差 分法 的基础 上 ，通 过 引入界 面 中心 点压力 ，提 出 了一 种新 的数值 算 法．该 算法 基? 于单 个 网格单 元来构 造数 值计 算格 式 ，能够精 确 地获取 每一 网格 单元 的压力 梯 度 ，单 元 间则 通 过流 量连续 条? 件联 系起来 ．并 以此 为基 础 ，建 立 了一套 非均 质低 渗透 油藏 非线性 流 动数值 模拟 方法 ．? ? 渗 流 数 学 模 型? 由于受 启 动 压 力 梯度 的影 响 ，低 渗 透 油 藏 流体 流 动存 在