高速非达西气井渗流模型及压力曲线特征.pdfVIP

第 37卷第 4期 石 油 钻 探 技 术 Vo1．37．NO．4 PETR()IEUM DRIIIING TECHN1QUES Ju1．，2009 ● “863”计划专栏 高速非达西气井渗流模型及压力曲线特征 马 婧 程时清 (中国石油大学 石油工程教育部重点实验室 ，北京 昌平 102249) 摘 要 ：大量研究表明，气井在井筒周围的高速流动规律不符合达西定律 。引入非达西修正系数 ，建立了考虑 井筒储集和表皮效应的高速非达西气井试井模型，给 出了相应 的有限差分解法 ，绘制 了合非达西流动系数 的试井 典型曲线图版 。与达西流动的 曲线相比，存在高速非达西流动 效应 的试 井导数 曲线驼峰后 移，水平段 出现 时间延 后 ，非达西流动系数越大，其导数 曲线峰值越大，出现径 向流 的时间越晚，表明惯性效应本身是一种集 中于井筒的 附加流动阻力，引起 附加压降 ，气井 出现径 向流的时间也 变长 。 关键词 ：气井；试井；非达西流；数学模 型 中图分类号 ：TE319 文献标识码：A 文章编号 ：1001—0890(2009)04—0028-04 大量研究表明，由于气体的黏度低，所以实际上 温渗流且存在高速非达西效应；3)气体为真实气体 ， 气体的流动速度非常高，特别是在压力梯度最高的 忽略岩石弹性的影响；4)考虑井筒存储效应和表皮 井筒附近[1]。在这些条件下，惯性力将导致偏离线 效应 。 性的达西定律 ]。在其他因素不变的条件下，气井 1)运动方程 产量越高 ，非达西效应越强；渗透性越差，非达西效 Forcheimer建议用二次流动项来扩展达西定 应越强_l3]。非达西流 已经成为影响气藏开发的主要 律 ，以便将惯性力的影响包括进去 ]： 因素之一 ，特别是在高产气藏渗流 中尤为明显。对 ll+ “ -- gradp (1) 非达西行为描述的出发点是基于 Forcheimer方程 ， 式 中， 为气体黏度 ，mPa ·S；K 为地层渗透率 ， 利用该方程来弥补流体高速流动时达西定律的不 10一 m ；“为渗流速度，m／s；为非达西流动系数 ， 足 ，与达西方程相比该方程增加了一个非达西修正 无因次 ；D．为气体密度，kg／m。；P为气藏压力 ，MPa。 系数 。 按 Wattenbarger等人 的观点 ，将 Forcheimer 长期以来，气井试井和渗流力学计算均是通过 定律方程作为达西定律的修正公式表示出[引： 引入拟压力 ，将气藏渗流中含压力、黏度和压缩因子 K 的项线性化。由于气井渗流含高速非达西项 ，渗流 U～ grad． (2) 方程为含二次项的非线性方程 ，不能应用线性叠加 式中，是与流动速度大小有关的非达西流动效应 原理。为此，笔者通过引入非达西流动系数 ，建立了 繇 1 删 高速非达西气井试井数学模型，给出了相应的有 限 1+ 肋 “ 差分解法，并编制了计算机程序 ，绘制了含非达西流