考虑动半径和粘度变化的稠油非牛顿不稳定渗流数学模型

计算物理 ›› 2011, Vol. 28 ›› Issue (6): 869-874.

考虑动半径和粘度变化的稠油非牛顿不稳定渗流数学模型

罗艳艳, 程林松, 黄世军

中国石油大学石油工程学院石油工程教育部重点实验室, 北京 102249

收稿日期:2010-12-03 修回日期:2011-03-28 出版日期:2011-11-25 发布日期:2011-11-25
作者简介:罗艳艳(1984-),女,博士生,主要从事油气田开发研究,E-mail:2002hoyy@163.com
基金资助:
国家973项目(2011CB707305);国家科技重大专项(2011ZX05024-005-006);国家重大专项(2011ZX05012-004)资助项目

Mathematical Model for Transient Flow of Non-Newtonian Heavy Oil with Dynamic Radius and Variable Viscosity

LUO Yanyan, CHENG Linsong, HUANG Shijun

MOE Key Laboratory of Petroleum Engineering, College of Petroleum Engineering, China University of Petroleum, Beijing 102249, China

Received:2010-12-03 Revised:2011-03-28 Online:2011-11-25 Published:2011-11-25

摘要/Abstract

摘要： 针对稠油非牛顿特征,在Bingham流体渗流方程基础上,通过对动半径和粘度进行表征,建立同时考虑启动压力梯度、动半径变化和粘度变化的非牛顿稠油不稳定渗流数学模型,完善Bingham型稠油渗流数学模型.通过空间、时间离散差分及Matlab数值计算,得到非牛顿稠油非稳态渗流地层压力分布.结果表明,相同产量下,随启动压力梯度增大,动半径向井方向移动;启动压力梯度越大,压降曲线越陡,相应近井压降越大;相同启动压力梯度下,产量越大,不同吞吐周期压力差距越大.将半径和粘度动态变化相结合,弥补了现行非牛顿稠油渗流数学模型的-个缺陷.

关键词: 动边界, 非牛顿稠油, 不稳定渗流, 压力分布, 启动压力梯度

Abstract: A mathematical model for transient flow of non-Newtonian heavy oil with dynamic radius and variable viscosity is established.The non-Newtonian transient pressure distribution of heavy oil is obtained with time and space.It shows that under same production,dynamic radius is closer to the well as threshold pressure gradient increases.The greater the threshold pressure gradient, the greater the pressure drop neflr the well.We combine dynamic radius and viscosity to avoid shortcomings of existing non-Newtonian heavy oil flow model.

Key words: dynamic radius, non-Newtonian, transient flow, pressure distribution, threshold pressure gradient

中图分类号:

TE345

罗艳艳, 程林松, 黄世军. 考虑动半径和粘度变化的稠油非牛顿不稳定渗流数学模型[J]. 计算物理, 2011, 28(6): 869-874.

LUO Yanyan, CHENG Linsong, HUANG Shijun. Mathematical Model for Transient Flow of Non-Newtonian Heavy Oil with Dynamic Radius and Variable Viscosity[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 28(6): 869-874.

[1]	王兆清, 钱航, 李金. 移动边界问题的时空域正则区域迭代配点法[J]. 计算物理, 2021, 38(1): 16-24.
[2]	刘化普, 刘慧卿, 王敬. 分形低渗透缝洞型油藏非线性渗流规律[J]. 计算物理, 2018, 35(1): 55-63.
[3]	曾现洋, 倪国喜. 振动NACA0012翼型的移动网格数值模拟[J]. 计算物理, 2016, 33(3): 266-272.
[4]	何聪鸽, 范子菲, 许安著, 方思冬, 李博凯. 特低渗透油藏面积井网见水时间计算[J]. 计算物理, 2016, 33(2): 190-196.
[5]	周春华. 一种基于预估-修正策略的非定常流网格自适应加密/稀疏方法[J]. 计算物理, 2013, 30(5): 633-641.
[6]	刘文超, 姚军, 王建忠. 低渗透多孔介质非达西渗流动边界界面追踪[J]. 计算物理, 2012, 29(6): 823-827.
[7]	王强, 胡湘渝, 姜宗林. 一种含运动固壁超声速流动的Descartes网格算法[J]. 计算物理, 2009, 26(4): 517-526.
[8]	邓树升, 耿继辉, 叶经方, 谭俊杰. 含运动变形物体流场的数值模拟[J]. 计算物理, 2009, 26(2): 200-210.
[9]	刘演华, 林建忠, 库晓珂. 格子Boltzmann方法中三维运动边界的统一模型[J]. 计算物理, 2008, 25(5): 535-542.
[10]	王兵, 许厚谦. 含有大位移动边界的复杂流场的数值模拟[J]. 计算物理, 2008, 25(4): 396-400.
[11]	张烈辉, 张锦良, 徐冰青. 变形双重介质分形气藏非线性渗流理论模型及数值研究[J]. 计算物理, 2007, 24(1): 90-94.
[12]	张玉东, 纪楚群. 包含动边界的非定常流场动网格数值模拟[J]. 计算物理, 2006, 23(2): 165-170.
[13]	倪学锋, 程林松, 李春兰, 安九泉. 注蒸汽井井筒内参数计算新模型[J]. 计算物理, 2005, 22(3): 251-255.
[14]	吴兆春. 定时间步长变坐标步长差分求解单相Stefan问题[J]. 计算物理, 2003, 20(6): 521-524.
[15]	张攀峰, 王晋军. 开槽扰流片流动特性研究[J]. 计算物理, 2003, 20(3): 226-232.