《气体管流及嘴流动态优秀课件.ppt》由会员分享,可在线阅读,更多相关《气体管流及嘴流动态优秀课件.ppt(37页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、气体管流及嘴流动态第1 页,本讲稿共37 页气体管流及嘴流动态气体通过油管和油嘴的流动是气井生产系统中的重要流动过程。掌握可压缩介质天然气的稳定一维管流的基本方程,是分析气井系统中单相气体管流和气液两相管流规律的基础。重点介绍干气井井底静压、流压和输气管线流量的计算方法。第2 页,本讲稿共37 页第一节第一节气相管流基本方程气相管流基本方程将 气 相 管 流 考 虑 为 稳 定 的 一 维问 题。在 管 流 中 取 一 控 制 体(如图 示),以 管 子 轴 线 为 坐 标 轴z,规 定 坐 标 轴 正 向 与 流 向 一 致。定义 管 斜 角 为 坐 标 轴z 与 水 平 方 向的夹角。v+
2、dvz+dzdzzp+dpwpvzgAdz稳定一维气相流动第3 页,本讲稿共37 页气相管流基本方程气相管流基本方程1.连续方程假设无流体通过管壁流出和流入,由质量守恒得连续性方程即G=vA=常数上式表示任意管子截面z 上气体质量流量均保持不变。式中 气体密度,kg/m3;v 气体流速,m/s;A管子流通截面积,m2;D管子内径,m;G 气体质量流量,kg/s;v流过单位截面积的气体质量流量,kg/(m2.s)。对于等径油管,v 为常数。(1)第4 页,本讲稿共37 页2.动量方程作用于控制体的外力应等于流体的动量变化:作用于控制体的外力Fz包括:质量力(重力)沿z 轴的分力gAdzsin压力
3、pA(p+dp)A管壁摩擦阻力(与气体流向相反)wDdzw流体与管壁的摩擦应力(单位面积上的摩擦力),Pa;D控制体的周界长,m;p压力,Pa;g 重力加速度=9.81m/s2;管斜角,度。(2)第5 页,本讲稿共37 页将上述三项外力代入式(1)得实验表明,管壁摩擦应力与单位体积流体所具有的动能成正比。引入摩阻系数f,即摩阻项可表示为(3)第6 页,本讲稿共37 页动量方程即为压力梯度方程上式总压降梯度可用下式表示为三个分量之和,即重力、摩阻、动能压降梯度(分别用下标g、f 和a 表示)。其中动能项较前两项甚小,在工程计算中往往可忽略不计。第7 页,本讲稿共37 页应该强调,方程(3)的坐标
4、轴z 的正向与流体流动方向一致。管子的倾角 规定为与水平方向的夹角,对于垂直气井=90,sin=1。在气井管流计算时往往是已知地面参数,计算井底静压和流压,习惯上是以井口作为计算起点(z=0),沿井身向下为z 的正向,即与气井流动方向相反。此时,压力梯度取“+”号。(4)第8 页,本讲稿共37 页单相管流摩阻系数流 体 在 管 内 流 动,往 往 使 部 分 机 械 能 转 换 为 热 能 而 造 成 不可 逆 的 能 量 损 失。在 单 相 流 动 的 情 况 下,不 可 逆 损 失 主 要 是 摩擦 损 失,此 项 损 失 包 括 由 于 流 体 粘 滞 性 产 生 的 内 部 损 失 和
5、 管 壁形 成 的 外 部 损 失。除 层 流 外,实 际 能 量 损 失 无 法 由 理 论 计 算 确定,而 是 采 用 实 验 方 法 和 相 关 分 析 确 定 摩 阻 系 数f。摩 阻 系 数是 一 个 无 因 次 量,它 反 映 了 管 壁 剪 切 应 力 对 摩 阻 压 降 的 影 响 程度。摩 阻 系 数 是 雷 诺 数NRe和 相 对 粗 糙 度e/D 的 函 数,摩 阻 系 数f 与NRe为 双 对 数 关 系。可 分 为 四 种 流 态:层 流、临 界 流、过 渡流和紊流。第9 页,本讲稿共37 页图2 Moody 摩阻系数第10 页,本讲稿共37 页雷诺数是一无因次量,
6、定义为式中 为流体粘度,Pa.s。其余符号意义和单位同前。雷诺数表示流体惯性力与粘滞剪切力之比值,它是判别层流与紊流的重要参数。通常认为,层流与紊流的分界雷诺数为21002300。第11 页,本讲稿共37 页实际管子内壁的粗糙度e 与管子材质、制造方法以及腐蚀和结垢等情况有关。管壁粗糙度的取值往往比较困难,因为其值不是可直接测量的参数。而是根据测试的压力梯度计算其摩阻系数,由Moody 图反求有效的e/D 值。对于新油管推荐e=0.016 mm(0.0006in)。Jain(1976)摩阻系数公式第12 页,本讲稿共37 页第二节气相管流压降计算应用气体稳定管流压降梯度方程可以解决井筒或输气管
7、线的实际工作问题。根据井口参数计算井底静压或流压 已知输气管线两端压力计算其输气能力 一般的问题是根据管子两端压力或流量中的两个参数计算另一个参数由于气相管流压力梯度方程目前还不能直接用解析显式表示压力,所以发展了多种不同程度简化和近似的方法。第13 页,本讲稿共37 页一、井底静压计算根据井口参数计算井底压力,取坐标z 沿井轴向下为正,井口z=0。垂直井:测深L 等于垂深H,90,sin=1;斜直井:sin=H/L 关井静气柱:v=0(气体不流动),压降梯度方程(4)中摩阻项和动能项均为零。垂直井静气柱总压降梯度即为重位压降梯度第14 页,本讲稿共37 页式中气体密度,kg/m3;Mg天然气
8、视分子量=28.97 g,kg/kmol;g天然气相对密度;R 通用气体常数=8315Pa.m3/(kmol.K);p压力,Pa;T 温度,K;Z 气体偏差系数。垂直井静气柱压降梯度第15 页,本讲稿共37 页分离变量并积分由于T,p,Z 是沿井深变化的,为了便于直接积分,采用井筒平均温度和平均压力计算平均Z 值,积分得即第16 页,本讲稿共37 页式中pws气井井底静压,Pa;pwh气井井口静压,Pa;H井口到气层中部深度,m。井筒气柱平均温度=(Twh+Tws)/2,K;井筒气柱平均偏差系数;s指数,无因次。井底静压为第17 页,本讲稿共37 页二、井底流压计算仍以井口为计算起点,沿井深向
9、下为z的正向,与气体流动方向相反。忽略动能压降梯度,垂直气井的压力梯度方程为任意流动状态(p,T)下的气体流速可表示为式中vsc标准状态下气井流速,m/s;v任意位置处流动状态下的气体流速,m/s;psc,Tsc标准状况的压力,温度,psc=0.101MPa;Tsc=293K;qsc气井日产气量(标准状态),m3/d。第18 页,本讲稿共37 页分离变量积分第19 页,本讲稿共37 页式中 pwf气井井底流压,MPa;pwh气井井口流压,MPa;f T,p下的摩擦系数;qsc标准状态下天然气体积流量,m3/d;D 油管内径,m。其它符号及其单位与静压计算公式相同。井底流压为 第20 页,本讲稿
10、共37 页算例某垂直气井井深3000m,井口关井压力为2MPa,井筒平均温度为50,ppc=4.6MPa,Tpc=205k,g=0.65。1)计算井底静压;2)气 井 开 井 生 产,油 管 内 径 62mm,产 气 量 范 围 为1050104m3/d,计算井底流压的变化。第21 页,本讲稿共37 页井深(m)静气柱(qSC=0)流动气柱(qsc,104m3/d)10 20 30 40 500 2 2 2 2 2 2300 2.043 2.256 2.785 3.482 4.261 5.081600 2.087 2.495 3.411 4.522 5.704 6.912900 2.132 2
11、.722 3.958 5.390 6.880 8.3851200 2.179 2.941 4.455 6.161 7.914 9.6761500 2.226 3.152 4.920 6.871 8.862 10.8561800 2.275 3.358 5.361 7.541 9.753 11.9642100 2.324 3.561 5.786 8.181 10.604 13.0242400 2.375 3.761 6.198 8.801 11.427 14.0502700 2.427 3.960 6.602 9.406 12.23 15.0533000 2.4804.157 6.999 10.
12、001 13.021 16.040静、动气柱压力(MPa)第22 页,本讲稿共37 页静、流压分布曲线第23 页,本讲稿共37 页三、水平输气管线压降设 水 平 输 气 管 线 其 流 动 方 向 与 水 平 方 向 一 致,无 高 程 变 化 故不 存 在 重 位 压 降。忽 略 动 能 压 降,故 总 压 降 梯 度 为 摩 擦 压 降 梯度分离变量积分式中qsc气流量,m3/d;p1,p2分别为管线起点和终点压力,MPa;L水平管线长度,m。第24 页,本讲稿共37 页算例计算水平集气管线压降。假设产气量变化范围为10100104m3/d,计算对比50.3、62、73、106mm 管径,
13、1000m 管线的压降。已知井口起点压力为15MPa,平均温度为20,其它数据同上例。第25 页,本讲稿共37 页表3-2 1000m水平管线的摩阻压降(MPa)气流量(104m3/d)管径(mm)50.3 62 73 10610 0.205 0.070 0.030 0.00420 0.824 0.274 0.118 0.01730 1.912 0.619 0.263 0.03840 3.603 1.114 0.468 0.06750 6.239 1.775 0.735 0.105第26 页,本讲稿共37 页1000m 水平管线摩阻压降第27 页,本讲稿共37 页四、凝析气的修正对于含有重烃的
14、凝析气井,在油管的举升过程中会部分冷凝成液相,形成气液两相流动。由于这类气井的气液比很高且紊流程度严重,气液两相混合较为均匀,可视为均匀的单相流(称为拟单相流)。凝析气井的产出物包括三部分:即经地面分离器分离出的干气、凝析油罐逸出的凝析气和凝析油。当井底流压接近凝析气的上露点压力,油管内可能存在液烃。若气液比大于1780m3(标)/m3,可近似考虑为单相气体流动。计算井底流动压力须对气体流量 气体流量和相对密度 相对密度作相应的修正。第28 页,本讲稿共37 页凝析气总气量的修正式中qT修正后的总气量,m3/d;qSG 经分离器得到的干气产量,m3/d;qO 凝析油产量,m3/d;qEG 凝析
15、油的相当气相体积,m3/m3;qTG 凝析油罐逸出的凝析气量,m3/d;o凝析油罐内凝析油的相对密度;Mo凝析油罐内凝析油的平均分子量。第29 页,本讲稿共37 页凝析气复合相对密度考虑了凝析油物性和含量的综合影响,修正公式为 w复合气体相对密度;Rg地面总生产气油比,m3/m3;g地面分离器和凝析油罐气的平均相对密度;SG分离器的干气相对密度;TG凝析油罐逸出气相对密度。式中第30 页,本讲稿共37 页第四节 气体嘴流动态许多自喷井在井口都要安装节流装置油嘴,用于控制气井的产量。有多种情况要求限制气井的产量,包括防止底水锥进和地层出砂,通过调节油嘴的大小控制井口压力以满足地面设备的耐压要求或
16、防止生成水化物。节流部件种类很多,包括井口油嘴或针形阀,安装在油管鞋附近的井下油嘴,油管上部的井下安全阀。(SSSV),气举阀的气孔等。当气流通过这些流通截面突缩部件时,其流动规律基本一致,可概括为嘴流。第31 页,本讲稿共37 页左图示意圆形孔眼的油嘴,若上游压力p1保持不变,气体流量(标准状态下)将随下游压力p2的降低而增大。但当p2达到某值pc时,流量将达到最大值即临界流量。若p2再进一步降低时,流量也不再增加。流量与油嘴上下游压力比的关系如右图所示。p2p1dq0pc/p1 1p2/p10qsc嘴流示意图嘴流动态关系第32 页,本讲稿共37 页所谓“临界流”是流体在油嘴吼道里被加速到声
17、速时的流动状态。在临界流状态下,油嘴下游压力变化对气井产量没有影响,因为压力干扰向上游的传播不会快于声速。因此,为了预测嘴流动态即产量与节流压降的关系,必须确定是否为临界流状态。根据热力学原理,临界压力比为式中k为气体绝热指数当时,为临界流;否则为亚临界流第33 页,本讲稿共37 页根据气体嘴流的等熵原理,流量与压力比的关系为下式,对于亚临界流状态式中qsc通过油嘴的体积流量(标),104m3/d;p压力,MPa;d嘴眼直径,mm;T温度,K;下标1、2分别表示嘴前、嘴后位置;p2/p1压力比。第34 页,本讲稿共37 页对于临界流,嘴流最大产气量为第35 页,本讲稿共37 页算例油嘴直径分别用4、5、6、7、8mm。气体相对密度为0.7,k 为1.25,气体偏差因子为0.93,井口温度和压力分别为38C和4MPa。绘制不同嘴径的气产量与压力比的特性曲线。解:天然气的临界压力比为0.555。对每一油嘴尺寸,从0.555到1改变p2/p1值,可绘制嘴流曲线(见图)。不同油嘴尺寸的临界流量d,mm 4 5 6 7 8qmax,104m3/d0.85 1.33 1.92 2.61 3.41第36 页,本讲稿共37 页不同油嘴尺寸的嘴流动态第37 页,本讲稿共37 页
限制150内