低渗透储层气体渗流数学模型及计算方法研究.docVIP

*本文为中国石油天然气集团公司 九五 重点科研项目 ! ! ! 低速非达西渗流机理研究部分内 容。 **李铁军 , 副教授 , 1964年出生 ; 1985年毕业于四川 大学数学系 , 1996年获油气田开发 硕士 , 现 从事教学 和数值模拟 方 面的科研工作。电话 :(0817 2643345。地址 :(637001 四川省南充市西 南石油学院计科系。 低渗透储层气体渗流数学模型及计算方法研究 * 李铁军 * * 李 允 (西南石油学院 李铁军等 . 低渗透储层气体渗流数学模型及 计算方法研究 . 天然气工业 , 2000; 20(5 :70~72 摘 要 低渗透气藏在我国分布极其广泛 , 据不完全 的统计 , 其储量已占所有气藏的 80%, 显而易见开发此 类 气藏具有重要的现实意义 。 然而 , 到目前为止 , 一般都采用常规气藏的理论来解释 低渗透的数 据 , 但 往往得不到 理 想的效果 。 这是因为低 渗透气藏具有常规 气藏 不具 备的特 殊性 , 岩 层致密 、 渗 透率 极低 、 气体 渗流 具有滑 脱效 应 等 ? 1、 2# 。 针对这种现状 , 文章在大量的国内外文献调研和实际计算的基础上 , 从研究低渗透储层气体渗流的特性 出 发 , 建立相应的数学模型 , 对方程中的系数 、 内外边界条件 进行特殊的处理 , 得到非 线性的差分 方程 , 采用牛顿迭 代 法求 解 。 通过 对川中 1口井数据的处理 , 文中模型及 其算法都 能较好 地拟合 气井的 实测数 据 , 这 项工作 为解决 此 类问题找到了一种行之有效的方法 , 而且 , 该成果容易推广到多维多相模型中去 。 主题词 低渗透油气藏 储集岩 气田开发 滑脱 数学模型 四川盆地 低渗透储层气体渗流数学模型 在建立模型时 , 需作以下假设 :? 无限大的气藏 中 1口气井以定产量生产 ; %气体渗流为等温渗流 且存在滑脱效应 ; 气体为真实气体 , 忽略岩石弹性 的影响。 1. 运动方程 V ? =- g (p grad (p (1 其中 : (p =1+p 。 2. 状态方程 = RTZ(p (2 其中 :R 、 T 为气体常数和温度 ; Z (p 为气体偏差系 数。 3. 质量守恒方程 ((p V ? + t =0(3 其中 : 为常数。 将 (1 、 (2 两式代入式 (3 , 化简即得考虑滑脱 效应的低渗透储层数学模型 : =A (p r r +B (p 2 (4 其中 :A (p =2Z(p -pZ (p , B (p =a 2 (p +Z (p -pZ (p , a 2=g 。 4. 定解条件 初始条件 p (r , t |t=0=p e (5 外边界条件 p (r , t |r ? ? =p e (6 内边界条件 r p r =r w = g q(t 2! K (p wf h (7 其中 :p 为气藏压力 , M Pa; b(p 为滑脱系数 , MPa; t 为生产时间 , d; q (t 为产 量 , 104m 3/d; r 为 气藏半 径 , m; r w 为气井半径 , m; p e 为气藏供给压力 , M Pa; p wf 为气井流压 , MPa; g 为天 然气粘度 , mPa (s; 为孔隙度 , 无因次 ; K 为岩石渗透率 , 10-3 m 2 。 方程 (4 与定解条件 (5 、 (6 和 (7 构成一个适 定的定解问题 , 该方程 (4 是一个强非线性的偏微分 方程 , 无法得到解析解。 数 值 求 解 1. 不均匀网格的形成 由于渗透区域内的压力分布与半径的关系不是 线性的关系 , 而是成对数的关系 , 故采用数值模拟柱 开发试采 天 然 气 工 业 2000年 9月 坐标系中常用的沿径向节点坐标公式 :? (ln r =ln r i+1-ln r i = ln r e -ln r i N -1 =Const (8 即 r i+1r i =r e r w 1 (9 式中 :r e 为区域半径。 要完全确定一个网格系统 , 除定义节点的坐标 外 , 还应该定义网格块的边界坐标 , 参考文献 ? 3#的 方法既不与节点的坐标发生矛盾 , 而且还满足物质 平衡的要求。即在计算传导系数时所用的块边界坐 标为 : r l i+1/2= r i+1-r i ln r i+1r i , r l i-1/2=r i -r i-1 ln r i r i-1而在计算网格体积时所用的块边界坐标为 :r 2i+1/2= r 21-2 i ln r i +1 r i 2 , r 2i-1/2= r 2 i r 21 ln r i r i-1 2 2. 差分方程形式 考虑式 (4 , 其中 A (p 、 B (p 是压力的函数