X非牛顿流体的差分方法研究.pdfVIP

第 22 卷 　第 1 期 　　　　　　　　　 Vol. 22 　No. 1 西 南 石 油 学 院 学 报 　　　　　　　　　 　2000 年 　2 月 　　　　　　　Journal of Southwest Petroleum Institute 　　　　　　　　　Feb 　2000 　 ( ) 　　文章编号 : 1000 - 2634 2000 01 - 0005 - 04 非牛顿流体的差分方法研究 1 1 1 1 2 2 张茂林 ,吴清松 ,李传亮 ,梅海燕 ,孙 雷 ,孙良田 ( 1. 中国科学技术大学热科学和能源工程系 ,安徽 合肥 230027 ; 2. 西南石油学院 , 四川 南充 637001) 摘要 : 用解析方法求解非牛顿流体的流动方程是非常困难的 ,工程上通常采用有限差分方法对其进行数值求解 。为 避免用势梯度分量差分方法求解非牛顿流体渗流问题所带来的偏差 ,提出了求解非牛顿流体渗流问题的势梯度模差 分方法 ;从理论上分析了势梯度分量差分方法存在的缺陷 , 即非牛顿流体的渗流微分方程与之相对应的分量差分方 程不相容 。最后以Bingham 原油为例揭示了二种方法在数值上的差别 。 关键词 : 非牛顿流体 ;渗流 ;数值模拟 ;差分方法 中图分类号 : TE319 　O357 　　　　　　　　　文献标识码 : A 引　言 Bingham 流体的本构方程可以写成 γ τ τ μ ( ) 　　 - 0 = 2 严格说来 ,所有流体都不是牛顿流体 。在某些 当流体在多孔介质中流动时 ,渗流速度与压力 特定条件下 ,流体的粘度随剪切速率变化不大 ,近似 梯度的关系曲线如图2 所示 。 表现出牛顿流体的性质 ,此时 ,把粘度处理成常数会 给工程计算带来许多方便 ;但在某些情况下 ,把流体 粘度处理成常数会给工程计算带来较大的误差 ,而 不得不选用非牛顿流体来描述流体的流变行为 。 大多数油藏流体都不是严格的牛顿流体 。如地 层原油表现出 Bingham 流体的性质[1 ,2 ] ; 聚合物驱 替溶液的粘度随剪切速率的增大而减小 ,表现出幂 律型流体的性质[3～6 ] 。 用解析方法求解非牛顿流体的流动方程是非常 图 1 　牛顿流体和 Bingham 流体的流变曲线图 困难的 ,工程上通常采用有限差分方法对其进行数 值求解[6 ] 。如何把流体的非牛顿行为用差分方法恰 当地表示出来 ,是油藏数值模拟需要研究的课题之 一 。本文以Bingham 流体为例 ,对非牛顿流体的差 分方法进行了研究 ,从理论上分析了势梯度分量差 分方法的错误 ,并提出了正确的势梯度模差分方法 , 最后给出了显示二种方法差别的数值计算实例 。 1 　达西方程 图2 　渗流速度与压力梯度的关系曲线图 牛顿流体和非牛顿 Bingham 流体的流变曲线如 图 1 所示 。牛顿流体的本构方程可以写成 对于牛顿流体可以写出下面形式的达西方程