李小川 ,胡亚非,张巍,陈明军,李强
(1.中国矿业大学 化工学院,江苏 徐州,221008;2.中国矿业大学 煤炭加工与高效洁净利用教育部重点实验室,江苏 徐州,221008)
由于除尘效果好,湿式除尘器在矿物开采和加工利用等高产尘环节应用广泛,对湿式除尘原理、除尘效率、除尘器阻力特性、气液两相流等理论的研究也较为丰富[1-6]。随着能源危机的凸显,流体流动的能耗与除尘效率的综合评价显得尤为重要,在较高的除尘效率、较低的阻力及设备安全工况区运行是湿式除尘器运行参数设定的更佳选择,目前这方面的综合研究与评述也备受关注[7-9]。目前对湿式除尘器内部复杂多相流动特性的研究方法或使用数值模拟或基于测试并不精确的气液比,更为有效的实验手段较少,很难对除尘器效率、阻力和安全工况区进行综合评价[10-13]。本文作者应用计算机多点实时测试技术和变频调速技术在较为宽泛的试验条件下对除尘器性能进行测试,已完成湿式除尘器内部压力分布和气液耦合特性等研究[14-15],在此基础上继续从除尘器运行时液相节流液位差、气相压损、除尘效率和除尘器共振等的变化规律出发进行研究,为探讨除尘器高效、低耗、安全的综合运行参数提供参考。
节流型自激式水幕除尘器结构原理如图1所示,含尘气流在通风机抽吸作用下流入除尘器,气流经导流器引导至节流器处加速后在节流口对除尘液面形成冲击,节流口前后产生节流液位差 Δh并同时产生大量高速运动的液滴,固相粉尘在冲击液面和与液滴接触时被捕捉,清洁气体经排气口排出。通过液床高度调节节流口初始过流面积,脱水器降低排气口液滴夹带量。
节流液位差 Δh是除尘器节流口液面与接触腔液面在节流作用下产生的液位差,如图1所示,Δh表征节流(包括气液固相的传质传动过程)产生的气相总压头损失,即
式中:p1和p2分别为节流口和除尘腔液面上方的静压力,Pa;x1和x2分别为节流口和除尘腔液面在x方向的最大刻度,m;dx为微分量;ρ为气体密度,kg/m3;g为重力加速度,kg/s2。
图1 节流型自激式水幕除尘器原理图Fig.1 Structure chart of wet dust collector
Δh与气体流速v(m/s)和初始过流高度b0(mm)有关,b0表示除尘器静止时节流器末端距液面的距离,当初始液位超过节流器最低点时b0为负值。
液相捕捉粉尘的主要方式有惯性碰撞、截留和布朗扩散等,由于气液接触方式和流速等不同这三种方式对除尘效率的贡献也有差异。对节流型除尘器而言,由于 Δh的影响,气流方向在节流口发生较大角度转变(图1中虚线表示气相流线),随着气流流速和b0的变化,粉尘捕捉过程也发生着微妙的变化,当气流流速较小b0较大时,气流对液面的冲击微弱,大颗粒粉尘由于惯性碰撞进入液面被捕捉,部分小颗粒粉尘由于布朗扩散进入液面被捕捉,而绝大多数小颗粒粉尘随气流流出除尘器而未被捕捉;当气流速度增大,Δh随之不断提高时,气流对液面产生较大冲击,激发大量液滴,此时气流方向发生较大角度的改变,几乎所有大颗粒粉尘由于惯性碰撞而被捕捉,由于气流方向转变较大,更多较小粒径粉尘处于截留作用范围内被捕捉,而由于大量液滴在除尘腔中形成水幕,增加了微细粉尘布朗扩散的捕捉几率,此时即使较大的b0也不会太多的降低粉尘捕捉效率;而当气流速度进一步增加时,惯性碰撞与截留作用并不会随气流速度的增大而继续增大,反而会因气流速度过快而缩短接触时间,降低了布朗扩散对微细粉尘的捕捉效率。
除尘器整体压力损失R(Pa)主要包括:含尘气流节流过程产生的节流损失RP,脱水器阻力RT、固相粉尘携带阻力RG和其他局部阻力Rf等。
除尘器内为气-液-固多相流动,伴随有复杂的传质传动过程,分别测量RP,RT,RG和Rf较为困难,在考察除尘器阻力特性时可将整个除尘器选为控制体,气相视为不可压缩流体,用进出控制体的气相全压差进行计算:
式中:R为所选控制体全压差;p1-1和p2-2分别为控制体进、出截面的静压力;v1-1和v2-2分别为控制体进、出截面的气流平均速度。
除尘器的长×入口宽×入口高为 2 m×0.4 m×0.2 m,除尘风机最大流量3 026 m3/h,节流器末端距除尘腔底高130 mm,入口设置直径为200 mm的均速管,安装TES-1341型热线风速仪测量入口流速,风速仪后方适当位置设置粉尘投料口,JYB-DW 型微差压传感器测量除尘器1-1和2-2截面处静压力,固定在节流口前后两侧的刻度尺测量节流液位差Δh,变频器调节除尘器流速,实验系统及测点布置如图2所示。
图2 实验系统及测点布置图Fig.2 Experimental system and measuring point arrangement
除尘效率实验所选粉尘为某烧结厂皮带输送转运点布袋收尘除尘灰中粒径小于106 μm的粉尘,其中,D16<11.67 μm,D75<38.97 μm,D90<53.03 μm,D98<71.11 μm,粉尘堆积密度1.13 g/cm3,FC-1A型粉尘采样仪采集除尘器出口粉尘浓度。
实验保持节流器高度 130 mm,粉尘投放量 0.2 kg/min,测量不同b0时除尘器节流液位差-气流速度(Δh-v)关系如图 3所示,全压损失-气体流速(R-v)关系和除尘器共振点如图4所示;测量除尘器出口粉尘浓度,采样时间 5 min,计算每次测试的除尘效率如图5所示。
由图3可知:同一初始过流高度b0条件的节流液位差Δh随入口流速v增加而增大,相同v时b0越大,Δh越低。另外,当b0为-15,0和14 mm时,Δh-v曲线呈显著的线性关系;当b0为34和45 mm时,Δh-v曲线呈二次曲线关系;在v<16.8 m/s,b0≤14 mm时,Δh-v曲线可近似表示为:
式中:k为与除尘器结构特点有关的系数,s/mm;Δh0为与初始液位b0有关的参数,mm。当b0>14 mm时Δh-v曲线用二次关系表示更为贴切。
由Δh-v曲线的变化规律可知:当b0较小时,除尘器在较低流速时就使气液两相较为充分的接触,液相的节流液位差形成了气相动能与压力能转换的抑制体,当v增大时,Δh与气相静压力形成线性平衡体系,表现为Δh随v的线性增长;而当b0较大时,除尘器在较低流速时气相对液相的作用较弱,液相在一定程度上以壁面摩擦阻力抑制气相动能与压力能的转换,因此Δh与v在一定程度上表现出了二次曲线关系的增长趋势。
图3 不同b0时Δh-v关系Fig.3 Δh-v curves of dust collector on different b0
在大节流条件下气流通过液相底层的阻力特性与气流通过固体边壁的传统模型(R-v曲线呈二次关系)和文献[11]中R-v曲线呈直线关系有较大差别。由图4可知:当b0较大时,R-v曲线体现出二次曲线关系,随着b0减小,线性关系逐渐凸现。由曲线拟合数据可知:当b0=45 mm时,拟合曲线二次项系数为1.290,而当b0为0和-15 mm时,拟合曲线线性相关度都高达0.998。
图4 不同b0的除尘器R-v曲线Fig.4 R-v curves of dust collector on different b0
除尘器阻力的二次曲线与直线关系主要是节流液位差 Δh对气相动能与压力能转换抑制的结果。除尘器运行时,在节流口前后液面分别加载了不同大小的静压力,产生了液位差值,随着v增大,节流口前后静压差增大,Δh随之增大,此时气体过流面积随Δh增大而增大,过流面积的增大抑制了流动阻力系数的进一步增大,最终全压损失与v达到一种线性平衡,Δh趋于稳定。具体而言,当b0较大v较小时,气体过流面积较大,液相以近壁面摩擦作用抑制气相的流动,Δh的抑制作用表现微弱,气相阻力体现出气体流过固体边壁的二次曲线特性;当v增大到一定值后,液相近壁面摩擦作用对气相阻力的影响远弱于 Δh的变化,除尘器R-v曲线体现出了直线特性,在图4中,当速度v分别大于7.33,10.64和11.95 m/s时,b0为45,34和14 mm曲线曲率都不断变小,体现出直线特性。当b0较小时,气体过流面积较小,在较低v时液相表面摩擦作用对气相阻力的影响就远小于 Δh的变化,当v增大时,Δh随之增大来降低流动阻力系数,使全压损失与v最终稳定在线性平衡状态,这从图4中b0为0和-15 mm两曲线得到证实。
湿式除尘器的除尘效率η与气液两相接触充分程度有关。由图5可知:随着流速增加,不同b0时的η都不断升高,b0较小时的除尘效率较高且有效流速范围宽。当b0为-15,0和14 mm时,在较小流速时除尘效率就能达到 90%以上,而当流速分别大于 7.5,8.8和9.2 m/s时,η就能达到98%,随着流速进一步增加,除尘效率能达到99.5%以上;而当b0为34和45 mm时,在较小流速时除尘效率不足90%,随着流速增加虽也能达到较高的除尘效率,但有效流速范围较窄,当b0=34 mm,v>12.3 m/s时,η能达98.6%,随着流速进一步增加η能达到99.4%,而当b0=45 mm,v=15.8 m/s时,η仅为98.2%,此时的除尘器阻力处于较高值。
图5 不同b0的除尘器η-v曲线Fig.5 η-v curve of dust collector on different b0
分析可知:当b0较小时,初始过流面积较小,在较小流速情况下就产生了较大的Δh,气流在节流口处受节流液位差影响,产生较大程度的方向改变,很大部分大颗粒粉尘都通过惯性碰撞被捕捉,得到了较高的除尘效率,随着流速增加,气流对液面的冲击增大,Δh进一步增大,且气流冲击的液滴量增加,在除尘腔中微细粉尘颗粒受布朗扩散作用影响被液滴捕捉,进一步提高了除尘效率。而当b0较大时,初始过流面积较大,需要有足够大的速度才能产生较大的Δh,因此惯性碰撞作用要在较大速度时才能起到较好作用,降低了整体除尘效率,而当流速进一步增加时,由于b0较大,同时 Δh随流速增加而增加,二者形成了很大的过流面积,较多小颗粒粉尘来不及发生方向转变就随气流流过了节流口,逃脱了液相的捕捉作用,除尘效率始终处于较低水平。
除尘器固有频率与除尘器风机转动频率会在一定转速下耦合发生共振,威胁到除尘器运行安全,实验知,除尘风机运行频率在44.8~45.6 Hz时,除尘器机械振动振幅明显加大,并发生异常声响,图4中“振动点”位置为各b0时除尘器开始发生共振现象的速度位置(速度分别为15.28,13.68,12.90,12.21和10.23 m/s),进一步分析可知:除尘风机叶轮旋转的不平衡惯性力是产生共振的激励作用力。因此,除尘器应避免运行在共振频率区间内。
在兼顾安全性的同时,较低的阻力损耗和较高的除尘效率是除尘器综合运行参数的重要评价指标,由图4和图5可知:当b0较小时除尘效率η很高,但除尘器全压损失R较大;当b0适中时除尘器全压损失R整体降低,但除尘效率η在速度达到一定值时都能达到较高值;而当流速增加到较高值时,所有b0的除尘效率增加缓慢,而阻力增加较快。综合考虑共振频率、阻力和除尘效率3个因素可知:当b0为0,14和34 mm时,入口速度分别为9.31~12.21 m/s,9.22~12.90 m/s和12.30~13.68 m/s是较为安全、高效和低能耗的运行工况区,相应的除尘器全压损失分别为 554~705,430~618和442~535 Pa,除尘效率能达到98.6%以上。
(1) 节流液位差对气相阻力的增加有抑制作用,除尘器Δh-v和R-v关系由于这种抑制作用,都表现出直线和二次曲线的组合,当初始液位b0较大时呈现出二次线关系,当b0较小时呈现出直线关系。
(2) 除尘效率η随b0减小而增大,相同b0时,η随气流流速增大而增大,液相通过调节 Δh更大程度地改变了含尘气流的流动方向,增大了粉尘惯性碰撞捕捉几率,Δh与流速v在提高除尘效率方面有着同样重要的作用。
(3) 除尘器在风机叶轮旋转的不平衡惯性力激励下会发生共振,综合共振频率、阻力和除尘效率3个因素,除尘器适宜的运行工况为b0在0~34 mm、入口流速在9.22~13.68 m/s,此时除尘效率高达98.6%,全压损失为430~705 Pa。
[1]Petavratzi E, Kingman S, Lowndes I.Particulates from mining operations∶ A review of sources, effects and regulations[J].Minerals Engineering, 2005, 18(2)∶ 1183-1199.
[2]于庆波, 王新华, 王庆林, 等.自激式除尘器的除尘性能[J].东北大学学报∶ 自然科学版, 2002, 23(8)∶ 799-801.YU Qingbo, WANG Xinhua, WANG Qingli, et al.Experimental study on the performance of self-swash dust catcher[J].Journal of Northeastern University∶ Natural Science, 2002, 23(8)∶799-801.
[3]刘建, 姚海飞, 魏传光, 等.掘进工作面湿式离心除尘器的结构优化及数值模拟[J].煤炭学报, 2010, 35(3)∶ 424-428.LIU Jian, YAO Haifei, WEI Chuanguang, et al.The structural optimization and numerical simulation of wet centrifugal precipitator at working face of diving airway[J].Journal of China Coal Society, 2010, 35(3)∶ 424-428.
[4]SHEN Xuefeng, WAN Jingsong, XUE Qingguo.Basic properties of steel plant dust and technological properties of direct reduction[J].International Journal of Minerals, Metallurgy and Materials, 2011, 18(3)∶ 227-284.
[5]马中飞, 贾勇.聚丙烯酰胺复配物在自激式除尘器中的应用[J].江苏大学学报∶ 自然科学版, 2007, 28(2)∶ 147-150.MA Zhongfei, JIA Yong.Application of polyacrylamide complexation on self-swash dust catcher[J].Journal of Jiangsu University∶ Natural Science Edition, 2007, 28(2)∶ 147-150.
[6]蒋仲安, 金龙哲, 杜翠风.湿式除尘器中脱水装置性能的理论和实验[J].北京科技大学学报, 2001, 23(3)∶ 196-198.JIANG Zhongan, JIN Longzhe, DU Cuifeng.Theoretical and experimental investigation of the characteristic of dewater equipment on wet scrubbers[J].Journal of University of Science and Technology Beijing, 2001, 23(3)∶ 196-198.
[7]Dubinskaya F E.A priori estimation for coefficients in the energy theory of wet dust trapping[J].Chemical and Petroleum Engineering, 2001, 37(9/10)∶ 525-529.
[8]赵建涛, 黄戒介, 吴晋沪, 等.气固错流移动颗粒床过滤器除尘效率[J].化工学报, 2004, 55(5)∶ 721-726.ZHAO Jiantao, HUANG Jiejie, WU Jinhu, et al.Collection efficiency of cross-flow moving granular bed filter[J].Journal of Chemical Industry and Engineering (China), 2004, 55(5)∶721-726.
[9]YANG Kejun, CAO Shuyou, LIU Xingnian.Flow resistance and its prediction methods in compound channels[J].Acta Mech Sin,2007, 23(1)∶ 23-31.
[10]袁文麒, 刘遂庆.管道充水工况下气液两相流瞬态数值模拟[J].同济大学学报∶ 自然科学版, 2010, 38(5)∶ 709-715.YUAN Wenqi, LIU Suiqing.Transient numerical modeling of gas-liquid two-phase flow in water-filling pipelines[J].Journal of Tongji University∶ Natural Science, 2010, 38(5)∶ 709-715.
[11]蒋仲安.湿式除尘机理的研究与应用[D].北京∶ 中国矿业学院北京研究生部, 1994∶ 84-88.JIANG Zhongan.The study of the mechanism of wet dust suppression and its application[D].Beijing∶ China University of Mining & Technology.Beijing Graduate Program, 1994∶ 84-88.
[12]吴志均, 何顺利.低气液比携液临界流量的确定方法[J].石油勘探与开发, 2004, 31(4)∶ 108-111.WU Zhijun, HE Shunli.Determination of the critical liquid carrying flow rate at low gas liquid ratio[J].Petroleum Exploration and Development, 2004, 31(4)∶ 108-111.
[13]Minami N, Murase M, Utanohara Y, et al.Numerical simulations on countercurrent gas-liquid flow in a PWR hot leg with air-water flow in a 1/15th scale-model[J].Journal of Energy and Power Engineering, 2011(5)∶ 495-503.
[14]李小川, 胡亚非, 蒋炎炎, 等.基于脉动压力时均值的除尘器气液耦合研究[J].煤炭学报, 2011, 36(8)∶ 1412-1416.LI Xiaochuan, HU Yafei, JIANG Yanyan, et al.Study of gas-liquid coupling in dust separator based on time-averaged value of pustule pressure[J].Journal of China Coal Society, 2011,36(8)∶ 1412-1416.
[15]熊建军, 骆振福, 卓卫民, 等.两相对流自激式水幕除尘器内压力分布研究[J].中国矿业大学学报, 2010, 39(6)∶ 902-906.XIONG Jianjun, LUO Zhenfu, ZHUO Weimin, et al.Research on pressure distribution in dust chamber of two phase convective flow self-oscillating water curtain dust removal system[J].Journal of China University of Mining & Technology, 2010,39(6)∶ 902-906.