低渗透油藏流体渗流的数学模型及算法研究.pdfVIP

石油天然气学报 (江汉石油学院学报)2013年 12月 第35卷 第 12期 JournalofOilandGasTechnology (J．JPI)Dec．2013Vo1．35No．12 低渗透油藏流体渗流的数学模型及算法研究 陈忠 (长江大学一年级教学工作部，湖北荆州434025) [摘要]在低渗透油藏开采过程 中，流体在油藏 中的渗流特征常表现为非线性渗流特性，渗流方程 已不符 合达西定律。通过分析油水两相 的渗流特 征规律和 影响流体 非线性渗 流的 因素 ，建立 了低渗透 油藏 油水 两相非线性渗流模型；然后用全 隐式差分格式将渗流方程离散化 ，构造 了非线性渗流 的数值模 型；进而 利用拟牛顿法讨论 了数值模型的计算方法，并给 出了算法 的收敛性定理。 ， [关键词]低渗透油藏；油水两相非线性渗流模型；全隐式差分方法；拟牛顿法 [中图分类号]TE3l9 [文献标志码]A [文章编号]1000—9752 (2ol3)12—0130—04 目前 ，国内外通用的商业化油藏数值模拟软件在建立数学模型时，流体在地层的渗流均考虑为达西 渗流，渗流速度与压力梯度的关系为线性关系，或者只考虑了注入单相流体的非线性渗流 ；处理油水两 相相对渗透率 曲线时，均用线性渗流方程 ，未考虑非线性渗流对油水两相相对渗透率的影响；因此常规 油藏的数值模拟方法都不适合于解决非线性渗流问题 。为此，笔者通过研究低渗透油藏的非线性渗 流规律 ，建立油水两相渗流方程；再根据建立的渗流方程 ，结合状态方程 、连续性方程、辅助方程 以及 定解条件，建立油水两相非线性渗流数学模型 ；然后将数学模型用全隐式差分方法进行离散 ，建立相应 的差分方程组 ；最后利用拟牛顿算法建立中小型数值模型的计算方法 ，并对算法的收敛性进行分析。 l 低渗透油藏流体渗流的数学模型 由于低渗透油藏具有启动压力梯度 ]，通常不符合常规油藏的达西渗流规律，若采用常规油藏的处 理方法来处理低渗透油藏，得到的数学模型将会有较大误差 ]。下面将根据低渗透油藏的非线性渗流规 律 ，即油相的渗流速度与压力梯度的非线性关系，建立相应的数学模型。 1．1 模型假设 1)油藏中有油水两相流动 ，油相的流动符合非线性流动规律 ，水相流动符合达西定律。 2)油藏的渗透率变化与压力有关 ，是油藏压力的函数 。 3)考虑油水两相毛细管压力的作用和重力的影响。 4)油水和岩石都是微可压缩的，且具备向异性。 5)油藏是等温渗流。 1．2 数学模型 油相的流动方程： 一 a ( 。 f。。g D) (1) 0 水相的流动方程 ： KK rw( 2)w 一 一 p 一p g D) (2) w 式 中：一U。、 分别为油 、水相的流速 ，m／s；n为非线性因子；K为油藏多孔介质的绝对渗透率 ，mD；K…K 分别为油、水相的相对渗透率 ，1； 为 Hamilton算子；P。、p 分别为油 、水相的压力 ，Pa；／l。、 分别为油、 [收稿 日期]2013 09—05 [基金项目]国家自然科学基金项目 61273179)；湖北省教育厅重点项目 (。 [作者简介]陈忠 (1965一)，男，1986年大学毕业 ，博士 (后)，教授 ，现主要从事最优化理论与算法方面的研究和教学工作。 第 35卷第 12期 陈忠：低渗透油藏流体渗流的数学模型及算法研究 水相的黏度 ，mPa·s；l0。、D』 分别为油、水相的密度 ，kg／m。；g为重力加速度 ，m／s；D为标高，m。 油相连续性方程 ： (a qSo) 。一 ( ㈦ 水相连续性方程 ： ( )