锥度对天然气净化用旋风分离器流场影响.pdf
2 0 1 4年 第4 3卷 第1 2期第8页 石 油 矿 场 机 械 犗 犐 犔 犉 犐 犈 犔 犇 犈 犙 犝 犐 犘犕犈 犖 犜 2 0 1 4,4 3(1 2) 8 1 2 文章编号 1 0 0 1 3 4 8 2(2 0 1 4)1 2 0 0 0 8 0 5 锥度对天然气净化用旋风分离器流场影响 李 琦1, 2, 罗 敏1, 韩传军2, 杨 雪2, 管西旗2 (1.西南油气田分公司 输气管理处, 成都6 1 0 2 1 3; 2.西南石油大学, 成都6 1 0 5 0 0) 摘要 基于计算流体力学, 采用R NG犽 犲湍流模型和离散相模型, 研究了锥度在1 6 3~1 7 5 的天然 气净化用旋风分离器的压力场、 速度场分布以及分离效率。结果发现, 升气管入口附近的部分区域 速度和压力变化最大; 锥度对筒体及小锥体区域几乎没有影响, 而对升气管和大锥体的影响很大, 但是并不改变流场的整体分布规律。随着锥度的增加, 压降呈递减趋势。气体总速度与切向速度 的变化趋势相同, 均随锥度的增大而减小。在最小粒径为5μm时, 不同锥度下旋风分离器的分离 效率均为1 0 0%, 但是在顶板附近有不同程度的颗粒堆积现象。锥度为1 6 3 和1 6 6 时, 颗粒返混现 象很严重, 极易造成旋风分离器的顶板腐蚀穿孔, 因此不适于工程实际。在剩余的3种结构中, 速度 变化相差小, 从压降和体积方面考虑, 1 7 2 的锥角结构最优, 可较大程度减少材料耗损, 节约制造成本。 关键词 旋风分离器; 锥度; 流场 中图分类号T E 9 3 1. 1 0 1 文献标识码A 犱 狅 犻 1 0. 3 9 6 9/j 檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪 .i s s n . 1 0 0 1 3 4 8 2. 2 0 1 4. 1 2. 0 0 3 图4 套管偏梯形螺纹黑皮现象 2 结论 1) 工厂自身质量意识的提高, 采购方明确的 订货技术要求, 相关机构对油管、 套管标准化的进一 步完善, 将整体提升国内油管、 套管的质量水平, 降 低用户的使用风险。 2) 对于采购方和制造方都应该在产品订货制 造前做好相应沟通, 明确相关技术参数要求, 才能对 产品质量有更加准确的把握, 不会造成产品质量信 息的不对称。 3) 制造企业、 监督部门应针对标准中存在疑 问的产品技术要求, 验证、 收集相关产品的质量数 据, 通过有效分析和证明, 得出合理的油管、 套管产 品控制方法。 4) 相关工厂、 标准研究者应研究或通过验证 的大量数据, 向A P I标准化协会提出标准修订要 求。通过专业化协会组织广泛征集意见, 将标准中 含糊或未明确的技术条款得到有效确认, 让产品制 造质量控制有据可依, 进一步提升产品的质量。 参考文献 [1] A P IS P E C5 C T,S p e c i f i c a t i o nf o rC a s i n ga n dT u b i n g [S]. 9 t he d . 2 0 1 1. [2] G B/T 2 3 0 . 32 0 1 2, 金属洛氏硬度试验.第3部分 标准硬 度块(A、B、C、D、E、F、G、H、K、N、T标尺) 的标定[S]. [3] 李昱坤, 易晓明, 李京川, 等.国内非A P I油套管质量 现状浅谈[J].石油工业技术监督, 2 0 1 2,2 8(6) 1 3 1 5. [4] S Y/T 6 8 5 7. 12 0 1 2, 含H2S油气田环境下碳钢和低 合金钢油管和套管选用推荐做法[S]. 2 0 1 2 1 2 1 3. [5] A P IT R5 C 3, 石油和天然气工业套管、 油管、 钻杆和作 为套管或油管使用的管线管的性能公式及计算[S]. 2 0 0 82 2 2 7. [6] A P IS P E C5 B,S p e c i f i c a t i o nf o rT h r e a d i n g,G a u g i n g a n d T h r e a dI n s p e c t i o no fC a s i n g,T u b i n g,a n d L i n e P i p eT h r e a d s[S]. 1 5 t he d . 2 0 0 8. [7] 吕拴录, 姬丙寅, 杨成新, 等. 2 4 4. 5mm套管偏梯形螺 纹接头犔4长度公差分析及控制[J].石油矿场机械, 2 0 1 2,4 1(6) 6 3 6 6. 收稿日期2 0 1 4 0 7 1 0 作者简介 李 琦(1 9 8 2 ) , 女, 四川邛崃人, 工程师, 硕士研究生, 主要从事油气储运技术和装备方面的研究和技术管理工 作,E m a i l l_q i @p e t r o c h i n a . c o m. c n。 犈 犳 犳 犲 犮 狋 狅 犳犜 犪 狆 犲 狉犃 狀 犵 犾 犲狅 狀犉 犾 狌 犻 犱犉 犻 犲 犾 犱 犻 狀犆 狔 犮 犾 狅 狀 犲犛 犲 狆 犪 狉 犪 狋 狅 狉狅 犳犖 犪 狋 狌 狉 犪 犾犌 犪 狊犘 狌 狉 犻 犳 犻 犮 犪 狋 犻 狅 狀 L IQ i 1,2, L UO M i n 1, HANC h u a n j u n 2, YANGX u e 2, GUANX i q i 2 (1.犗 犳 犳 犻 犮 犲 狅 犳犌 犪 狊犜 狉 犪 狀 狊 犳 犲 狉,犜 犺 犲犛 狅 狌 狋 犺 狑 犲 狊 狋犗 犻 犾犪 狀 犱犌 犪 狊犉 犻 犲 犾 犱,犆 犺 犲 狀 犵 犱 狌6 1 0 2 1 3,犆 犺 犻 狀 犪; 2.犛 狅 狌 狋 犺 狑 犲 狊 狋犘 犲 狋 狉 狅 犾 犲 狌 犿犝 狀 犻 狏 犲 狉 狊 犻 狋 狔,犆 犺 犲 狀 犵 犱 狌6 1 0 5 0 0,犆 犺 犻 狀 犪) 犃 犫 狊 狋 狉 犪 犮 狋T h ec o m p u t a t i o n a l f l u i dd y n a m i c(C F D)w a s c o n d u c t e d t oa n a l y z e t h e e f f e c t o f t a p e r a n g l eb e t w e e n1 6 3 t o1 7 5 o n f l u i d f i e l d i nc y c l o n e s e p a r a t o r,w h i c hw a su t i l i z e d f o rn a t u r a l g a sp u r i f i c a t i o nb yu s i n gR NG犽 e p s i l o nt u r b u l e n tm o d e l a n dd i s c r e t ep h a s em o d e l i nF L U E NTf l u i da n a l y s i ss o f t w a r ep l a t f o r m. S o m ec o n c l u s i o n s c o u l db e r e c e i v e d f r o mt h en u m e r i c a l s i m u l a t i o n . T h e g r e a t e s t c h a n g e so fv e l o c i t ya n dp r e s s u r ew e r eo c c u r r e dn e a r t h ee n t r a n c eo f r i s e r;a n dt h e r ew a s a l m o s tn o i n f l u e n c eo nc y l i n d e ra n ds o m ec o n ea r e aw i t hd i f f e r e n t t a p e r s;h o w e v e r,t h e i m p a c to f t a p e ro nt h er i s e ra n dm o s t c o n ea r e aw a ss i g n i f i c a n t . O nt h eo t h e rh a n d,t h ed i s t r i b u t i o no f f l o w f i e l dw o u l dn o tb e a f f e c t e db y t a p e r a s aw h o l e . T h e s e p a r a t i o ne f f i c i e n c yo f g a s s o l i d i s 1 0 0%d u e t ot h e5μm m i n i m u mp a r t i c l es i z e,b u t t h e r e i ss t a c k e dp h e n o m e n o ni n s i d et h ec y c l o n es e p a r a t o r w h i c h i s s e r i o u sn e a r t h e r o o f . F i n a l l y,t h eb e s t s t r u c t u r eo f c y c l o n e s e p a r a t o rw a s c h o s ew i t h1 7 2 t a p e rw h a tm e a n t a 0. 9 2 5mc o n e l e n g t h. T h i s s t r u c t u r eh a d t h e l o w e s t p r e s s u r ed r o pa n d l e s sm a t e r i a l c o n s u m i n g . 犓 犲 狔狑 狅 狉 犱 狊c y c l o n es e p a r a t o r;t a p e ra n g l e;f l o wf i e l d 旋风分离器是一种利用离心力分离气流中固体 颗粒或液滴的设备。由于其结构简单, 造价低廉, 无 运动部件, 操作范围广, 受温度、 压力限制小, 分离效 率高等特点, 已经成为输气站场上应用广泛的气固 分离设备。然而, 从近年清管的污物、 下游设备运行 情况来看, 输气站场的旋风分离器除尘效果并不是 十分理想。例如, 在智能检测前期的清管通球过程 中, 排出了大量的污物; 同时, 由微粒粉尘冲刷作用 造成的下游调压、 计量设备失效, 运行寿命减短等问 题普遍存在。因此, 研究旋风分离器的流场, 提高其 分离效率显得尤为重要。 影响旋风分离器流场的因素包括内部的结构参 数和外部的环境条件, 近年来许多学者对其中的一 些方面进行了研究。例如, 曹晴云等[ 1]系统研究了 不同的入口面积、 排气管直径比、 排尘管直径比对旋 风分离器的速度和压降的影响; 赵宏强等[ 2]利用大 涡模拟, 详细介绍了分离器的三维速度场; 万谷军 等[ 3]对不同温度和压力下的旋风分离器速度场进行 了分析; 谷瑞青等[ 4]等研究了升气管插入深度对分 离器流场的影响; 葛坡[ 5]等提出了一种对称多入口 结构的旋风分离器, 并对其流场进行了研究; 朱小 兵[ 6]等提出了一种用于钻井柴油机尾气净化的旋风 分离器, 完成了现场试验。 目前, 涉及到旋风分离器锥体结构的研究也有 很多。X i a n gR o n gB i a o等[ 7]研究了锥体尺寸对用 于大气采样的小型旋风分离器的影响情况; 王清华 等[ 8]总结了不同的锥体结构及其影响; 查文炜、 李昌 剑等[ 9 1 0]的结构优化设计中涉及到了锥体的相对尺 寸。但是并没有涉及到锥度对其流场的具体影响。 锥度是旋风分离器的一个重要结构参数, 它连接着 圆筒体和锥体, 是结构突变的转折点, 因此对流场的 影响不可忽视。 本文以计算流体力学为基础, 利用F L U E NT 的R NG犽 e p s i l o n湍流模型和D PM离散相模型, 研 究不同锥度下旋风分离器的流场分布情况及分离效 率。旨在找出最佳锥度, 为旋风分离器的结构优化 提供参考。 1 计算模型 1. 1 计算方法 旋风分离器的流场计算模型一般由气相流场和 两相流场[ 1 1]组成。由于旋风分离器内流体所做的 是强旋流的湍流运动, 其雷诺数很高, 因此在计算过 程中, 连续相选择R NG犽 e p s i l o n湍流模型, 离散相 采用D PM模型。 1. 2 几何模型及网格划分 采用文献[ 1 2] 中的几何模型, 是1个直切式入 口结构的旋风分离器。其筒体直径为0. 3 6m, 入口 尺寸为0. 1 8m0. 0 9m, 排尘口直径为0. 1 0m, 排 气管直径为0. 1 4 4m, 锥度为1 7 2 , 升气管底面中心 为坐标原点。采用非结构网格对计算区域进行划 分, 减少网格数量, 提高质量。图1为旋风分离器的 几何结构和计算模型。 9 第4 3卷 第1 2期 李 琦, 等 锥度对天然气净化用旋风分离器流场影响 a 几何模型 b 网格划分 图1 旋风分离器的分析模型 1. 3 边界条件设置 由于是模拟输气站场上用旋风分离器的流场, 因此气相设置为甲烷, 粉尘用与天然气输气管道所 含粉尘粒径分布接近的8 0 0目的滑石粉, 其密度为 27 0 0k g /m 3[1 3]。流场计算时采用 S I MP L E算法、 P R E S T O压力插补格式和QU I C K差分格式。气 相入口 边 界 条 件 选 择V e l o c i t y_ i n l e t, 速 度 为2 0 m/s。 气相出口设置为o u t f l o w, 直接与大气相连。 颗粒的进入设置为面射流源, 其入口速度与气相相 同, 可以避免产生相对速度滑移。颗粒加入量0. 0 3 k g /s, 粒径分布满足R o s s i n R a mm l e r, 中位粒径取 1 0μm。升气管出口设置为逃逸(e s c a p e d) , 锥底灰 斗出口设置为捕捉( t r a p p e d) [1 4]。 理论上, 不同的锥度会对旋风分离器的压力场 和速度场造成影响。保持其他条件不变, 选择不同 的锥度进行数值模拟, 从而得到流场随锥度的变化 规律, 找出最佳锥度。 文献[ 1 1] 中给出了1 6种旋风分离器的结构, 通 过测量可以得到其锥度取值为1 6 3~1 7 5 。对文献 [ 1]~[1 5] 所取的旋风分离器结构进行验证, 可以知 道锥度均在这一范围内。表1列出了本文所取的锥 度, 并通过调整旋风分离器的锥体总长度来改变其 大小。 表1 旋风分离器的锥度及对应长度 锥度/ ( )1 6 31 6 61 6 91 7 21 7 5 锥长/m 0. 4 2 50. 5 2 10. 6 6 90. 9 0 01. 4 8 6 2 计算结果分析 2. 1 计算模型验证 数值求解和分析的结果需要与实际的现象、 数 据比较, 检验模型的合理性和适用性, 这对模型是否 有用非常关键。为了验证本文的模型, 将气相流场 的数值计算结果与试验结果进行对比, 用不同气相 入口流量( 控制在9 5 0~13 5 0m 3/ h) 下总压降的变 化与文献[ 1 4] 中的试验数据进行比较, 如图2所示。 从图2中可以看出, 模拟压降与试验压降的值基本 吻合。 图2 模拟压降与试验压降对比 2. 2 锥度对压力场的影响 压降是衡量旋风分离器性能的一个重要指标, 其值越小, 则能量损失越小。图3给出了不同锥度 对旋风分离器压力场的影响。 在图3 a中可以看出动压是通过旋流运动, 由静 压转化而成的。在轴向位置狕=0m( 轴向向下为 正) 即升气管入口附近, 动压最大, 静压最小, 说明此 处旋流运动最强烈。 a 锥度为1 7 2 下的流场压力分布 01 石 油 矿 场 机 械 2 0 1 4年1 2月 b 不同锥度下的动压分布 c 不同锥度下的静压分布 d 不同锥度下的总压分布 图3 锥度对压力场的影响 不同的锥度下, 压力分布走向一致, 但是数值上 有较大差异, 尤其是在升气管和锥体部分。狕>0. 5 m以后的分离区域, 由锥度引起的动压大小交替变 化, 规律不明显, 如图3 b所示, 说明锥体分离空间长 度的改变使流体旋流中心发生不规则变化。从图 3 c和图3 d可以看出, 静压和总压的分布比较规律, 其值总体上随锥度的增加而降低。在分离器升气管 出口处, 1 7 5 锥度下的静压最小, 表2给出了5个锥 度下旋风分离器的压降。可见, 随着锥度的增加, 压 降逐渐减小。 表2 不同锥度下旋风分离器的压降值 锥度/ ( )1 6 31 6 61 6 91 7 21 7 5 压降/P a 13 7 712 4 411 4 410 7 510 0 8 2. 3 锥度对速度场的影响 流体在旋风分离器内部做的是三维强旋流的湍 流运动, 流动状态相当复杂。旋转速度越高, 颗粒被 离心力送到壁面的概率越高, 分离效率也越好。图 4给出 了 不 同 锥 度 下 旋 风 分 离 器 内 部 速 度 场 的 分布。 从图4可以看出, 总速度和动压的分布相似, 进 一步说明了动压是由旋流运动产生。切向速度最 大, 占据主导地位, 其走向与总速度近似一致; 径向 速度在零点上下波动, 正负交替变化, 表现出较强的 不稳定性。 随着锥度的变化, 旋风分离器的三维速度场发 生了较大的变化。在狕=0m( 即升气管入口处) 的 附近区域, 总速度和切向速度有最大值。锥度增加, 在狕>0. 5m( 轴向向下为正) 的区域, 速度变化规律 性不强, 总速度出现多个波峰、 波谷。1 6 3 时升气管 出口的总速度最大, 为3 1. 6 8m/s, 而1 7 5 时最小, 为2 8. 6 1m/s。切向速度与总速度的分布规律相 似, 随着锥度的增加而减小, 从1 6 3 的2 3. 9 4m/s降 低到1 7 5 的1 9. 6 3m/s。从图4 d可以看出, 在前面 4个锥度下, 轴向速度的变化很接近。 a 锥度为1 7 2 的流场各向速度 b 不同锥度下的总速度 c 不同锥度下的切向速度 11 第4 3卷 第1 2期 李 琦, 等 锥度对天然气净化用旋风分离器流场影响 d 不同锥度下的轴向速度 图4 锥度对速度场的影响 3 锥度对分离效率的影响 分离效率是衡量旋风分离器性能的最重要的指 标。通过设置最小粒径为5μm, 最大 粒 径为1 0 μm, 得到不同锥度下旋风分离器的分离效率均为 1 0 0%。同时追踪了不同锥度下颗粒的流动轨迹, 如 图5所示。可以看出, 在θ=1 6 3 和θ=1 6 6时, 颗粒 出现了返混现象, 在顶板附近和锥底一直做旋流运 动, 而其余3种情况颗粒轨迹相差不大。 虽然没有颗粒随气体流出升气管, 但是分离器 内部有不同程度的颗粒堆积, 这种情况在旋风分离 器的顶板附近最为严重, 即大量的颗粒在此处做旋 流运动, 并未向下排出, 会对旋风分离器的使用性能 造成影响, 增加顶板的腐蚀率, 减短分离器的使用寿 命。可以看出, 锥度越小, 顶板附近的返混现象越 严重。 图5 锥度对颗粒轨迹的影响 4 结论 1) 锥度的变化对升气管和锥体部分区域流场 的影响很大, 而分离器筒体部分受锥度的影响相对 较小。 2) 随着锥度的增加, 压降呈递减趋势。 3) 旋风分离器流场的总速度与切向速度的变 化趋势相同, 均随锥度的增大而减小。 4) 粉尘粒径为5~1 0μm的情况下, 不同锥 度下旋风分离器的分离效率均为1 0 0%, 但是在顶 板附近有不同程度的颗粒堆积现象, 锥度为1 6 3 和 1 6 6 时最严重, 极易造成旋风分离器的顶板腐蚀穿 孔, 因此该锥度结构的旋风分离器不适用于工程实 际应用。 5) 锥度为1 6 9 、1 7 2 、1 7 5 结构的3种旋风分 离器, 速度变化仅相差1~2m/s, 故主要从压降方 面考虑。1 7 2 的锥角结构最优, 能较大程度减少材 料耗损, 节约成本。 参考文献 [1] 曹晴云, 姬广勤, 金有海, 等.不同结构参数下旋风分离 器气相流场的数值研究[J].流体机械,2 0 0 8,3 6(6) 3 4 3 8. [2] 赵宏强, 郭艳, 蒋海华.旋风分离器内速度场的数值模 拟研究[J].现代制造工程, 2 0 1 0(4) 6 2 6 6. [3] 万古军, 孙国刚, 魏耀东, 等.温度和压力对旋风分离器 内气相流场的综合影响[J].动力工程,2 0 0 8,2 8(4) 5 7 9 5 8 4. [4] 谷瑞青, 陶华东.升气管插入深度对旋风分离器流场影 响的数值模拟[J].化工时刊, 2 0 1 3,2 7(4) 6 8. [5] 葛坡, 袁惠新, 付双成.对称多入口型旋风分离器的数 值模拟[J].化工进展, 2 0 1 2,3 1(2) 2 9 6 2 9 9. [6] 朱小兵, 吕瑞典, 何娟, 等.钻井柴油机尾气净化装置研 制[J].石油矿场机械, 2 0 1 2,4 1(6) 7 5 7 7. [7] X i a n gR o n g b i a o,P a r kSH,L e eK W. E f f e c t so f c o n e d i m e n s i o no nc y c l o n ep e r f o r m a n c e[J]. J . A e r o s o lS c i, 2 0 0 1(3 2)5 4 9 5 6 1. [8] 王清华.旋风分离器结构改进的研究现状和发展趋势 [J].锅炉技术, 2 0 0 7,3 8(2) 5 9,2 7. [9] 查文炜, 葛友华, 倪文龙, 等.基于C F D的旋风分离器 的结构优化设计[J].机械工程师, 2 0 1 0(1 2) 4 4 4 6. [1 0] 李昌剑, 陈雪莉, 于广锁, 等.基于响应曲面法径向入 口旋风分离器的结构优化[J].高校化学工程学报, 2 0 1 3,2 7(1) 2 4 3 1. [1 1] 霍 夫 曼A C.旋 风 分 离 器 原 理、 设 计 和 工 程 应 用 [M].北京 化学工业出版社, 2 0 0 4. [1 2] 赵新学, 金有海.排尘口直径对旋风分离器壁面磨损 影响的 数值 模拟 [J].机 械 工 程学 报,2 0 1 2,4 8(6) 1 4 2 1 4 8. [1 3] 吴小林, 熊至宜, 姬忠礼.天然气净化用旋风分离器气 液分离性能[J].化工学报, 2 0 1 0,6 1(9) 2 4 3 0 2 4 3 6. [1 4] 赵新学, 金有海.基于C F D的旋风分离器壁面磨损数 值预测[J].石油机械, 2 0 1 0,3 8(1 2) 4 2 4 5. [1 5] 惠胜利, 姬广勤, 金有海, 等.循环流化床锅炉用旋风 分离器性能的试验研究[J].粉体加工与处理,2 0 0 8, 1 4(2) 4 2 4 4. 21 石 油 矿 场 机 械 2 0 1 4年1 2月