Processing math: 100%

应用于臭氧消毒系统的新型静态混合器结构设计

胡珊, 毛澍洲, 邱光宇, 许小红, 张波. 应用于臭氧消毒系统的新型静态混合器结构设计[J]. 环境工程学报, 2020, 14(11): 3201-3207. doi: 10.12030/j.cjee.201912068
引用本文: 胡珊, 毛澍洲, 邱光宇, 许小红, 张波. 应用于臭氧消毒系统的新型静态混合器结构设计[J]. 环境工程学报, 2020, 14(11): 3201-3207. doi: 10.12030/j.cjee.201912068
HU Shan, MAO Shuzhou, QIU Guangyu, XU Xiaohong, ZHANG Bo. Structural design of new static mixer for ozone disinfection system[J]. Chinese Journal of Environmental Engineering, 2020, 14(11): 3201-3207. doi: 10.12030/j.cjee.201912068
Citation: HU Shan, MAO Shuzhou, QIU Guangyu, XU Xiaohong, ZHANG Bo. Structural design of new static mixer for ozone disinfection system[J]. Chinese Journal of Environmental Engineering, 2020, 14(11): 3201-3207. doi: 10.12030/j.cjee.201912068

应用于臭氧消毒系统的新型静态混合器结构设计

    作者简介: 胡珊(1996—),女,硕士研究生。研究方向:水处理技术。E-mail:1069828090@qq.com
    通讯作者: 张波(1975—),男,博士,副教授。研究方向:水处理技术及环保工艺设备开发。E-mail:zhangbo@ujs.edu.cn
  • 基金项目:
    江苏省研究生科研与实践创新计划项目(SJCX19_0567);常州市科技支撑项目(CE20185005)
  • 中图分类号: TQ051.71

Structural design of new static mixer for ozone disinfection system

    Corresponding author: ZHANG Bo, zhangbo@ujs.edu.cn
  • 摘要: 臭氧水消毒过程中,液液混合均匀性是关键。为提高臭氧水消毒过程的混合均匀性,优化设计了一种旋流扩散静态混合器。通过数值模拟得到旋流扩散静态混合器中臭氧水浓度及压力分布,以旋流扩散静态混合器前后的不均匀系数及压损系数作为评价指标,探究混合元件间距、夹角等因素对臭氧水混合均匀性及压力损失的影响,从而优化混合器结构设计,并与传统SK 型混合器进行比较。模拟结果表明:最优结构布局参数为相邻蘑菇头间距40 mm,交错角0°,十字交叉板夹角120°。在此情况下,不均匀系数由1.92降至0.074,压力损失系数为4.96%。新型混合器相比传统混合器具有更好的混合效果,且十字交叉板夹角大小对不均匀系数、压力损失的影响较大。
  • 油浸泥土主要来源于石油开采过程产生的落地原油和含油矿渣、污泥、垃圾的堆置,以及采油废水滥排、石油泄露事件等[1-2]。石油污染使土壤的通透性降低,影响土壤中植物和微生物的生长。油浸泥土若不妥善处置,将造成污染扩散,影响周边环境和人体健康[3]。油浸泥土内部的轻质油分挥发出来,可进入到大气循环当中,并且随着大气的运动,造成污染扩散。油浸泥土中烃类可渗入地下水,导致地下水水质恶化。因此,油浸泥土已成为巨大的环境隐患,针对其处理处置已引起广泛关注[4-5]

    目前,油浸泥土的处置技术主要包括萃取技术、生物处理技术和热脱附技术。热脱附技术是一项利用燃气或者电能等能源对污染土壤加热,并对挥发后的污染物收集和处理的土壤修复技术,热脱附技术对有机污染土壤具有较好的修复效果[6-9]。根据修复方式,热脱附技术分为异位热脱附和原位热脱附技术,异位热脱附技术是将土壤开挖后进行处理,不影响土地的再利用。而原位热脱附技术是一种在污染场地上进行原地修复的技术,具有对土壤扰动小、二次污染可控、处理范围广等优点,尤其是对土壤污染较深、渗透低、黏性强、开挖难度大、不利于异位开挖修复的场地,原位热脱附技术优势更加明显[10-11]。热脱附技术作为一种重要的非燃烧技术,在有机污染土壤修复及含油污泥处置领域具有较好的应用前景。热脱附技术的影响因素较多,如反应条件、反应设备、土壤特性等[1, 12-15]。关于此技术,国外已有较多的研究和应用,但在我国起步较晚,相关研究较少,且机理研究也不够深入。张学良等[16]研究了热脱附技术对苯、氯苯等污染物的处理效果,热脱附对不同有机物的去除率均接近100%,且发现温度、停留时间、含水率、孔隙率等因素均可影响热脱附效果。王锦淮[17]开展了热脱附技术的中试实验,结果表明苯胺、氯苯、1,2-二氯苯和1,4-二氯苯等污染物接近完全去除。BAKER等[18]的研究结果表明,在300 °C的加热温度下,热脱附处理苯、乙苯、萘烯、甲苯效果较好。截至目前,国内外关于热脱附技术的研究及应用主要集中于多环芳烃类及苯系物类污染土壤的处理上,很少应用于油浸泥土的修复,因此,开展针对热脱附技术应用于油浸泥土的影响因素的研究十分必要[19]

    新疆油田是我国历史久远的大型采油区,石油开采后遗留大量油坑和油浸泥土。本研究在分析新疆某采油区油浸泥土性质的基础上,采用热解炉实验装置,对不同初始含油率(油浸泥土中的石油烃含量)、含水率的油浸泥土进行了热脱附模拟实验,通过测定不同加热温度和加热时间下的固相含油率、能耗等参数,研究了热脱附过程中加热温度、时间、初始含水率、含油率等对热修复效果的影响,为热脱附技术在油浸泥土修复领域的应用提供参考。

    油浸泥土样品采自中国西北新疆某油坑污染场地,采样深度0~30 cm。油浸泥土呈黑色、黏稠、半流质状,称取40 kg油浸泥土,将油浸泥土中的树叶、树枝和石子等异物除去,将油浸泥土充分搅拌、混合均匀,放置于密封箱中保存。原始油浸泥土样品外观见图1。理化性质如下:土壤质地为壤土,pH 7.5,含水率20%,总有机碳0.324 8 g·kg−1,阳离子交换容量23.5 cmol·kg−1

    图 1  原始油浸泥土样品
    Figure 1.  Raw soil sample

    本研究装置为自制热解炉小试装置,热解炉主要由加热系统、冷凝系统、液体收集系统、气体收集系统4个系统组成,见图2

    图 2  热解炉装置
    Figure 2.  Diagram of pyrolysis furnace

    工艺流程如图3所示。将一定量的供试泥土置于热解炉中,不断通入氮气,使热解炉内保持无氧状态。采用燃气加热的方式对炉内油浸泥土加热,通过温度自动控制系统进行温度调控,使炉内油浸泥土的温度达到实验设定温度。经过一定的停留时间后,将油浸泥土从出料口取出,并进行含油率检测,挥发出的高温混合气体进入冷凝系统,收集到的油水混合物静置分层,分别取上层油相和下层水相用于制备不同含油率和含水率的供试油浸泥土。

    图 3  热解炉装置原理图
    Figure 3.  Schematic diagram of pyrolysis furnace

    称取500 g(干质量)供试油浸泥土,从进料口加入炉筒中,分别设定不同的加热温度(200、250、300、350、400 °C)和停留时间(0.5、1、2、3、4、5、6 h),待出料冷却后,利用红外分光光度法进行固相含油率检测。每次实验设3组平行样,相对标准偏差控制在5%以内。

    按照一定比例,将热解炉小试装置回收的水、油添加于上述供试油浸泥土中,制备不同含油率(5%、10%、15%、20%、25%、30%)和含水率(5%、10%、20%、30%、40%、50%)的供试油浸泥土。在本研究中,设置加热温度300 °C、停留时间4 h的工艺条件,利用不同初始含油率和不同初始含水率的土壤进行热脱附处理,利用红外分光光度法对处理后油浸泥土进行固相含油率检测。每次实验设3组平行样,相对标准偏差控制在5%以内。

    为研究加热温度对于油浸泥土热脱附处理效果的影响,将供试油浸泥土(含油率11.3%,含水率20%)置于热解炉中,通过自动调温系统维持炉内油浸泥土温度为200、220、250、300、350、400 °C,停留4 h后卸料。待出料冷却后,测定固相含油率。

    为研究加热时间对于油浸泥土热脱附处理效果的影响,将供试油浸泥土(含油率11.3%,含水率20%)置于热解炉中,通过自动调温系统维持炉内油浸泥土温度为300 °C,分别在停留0.5、1、2、3、4、5、6 h后卸料。待出料冷却后,测定固相含油率。

    为研究油浸泥土初始含油率对于热脱附处理效果的影响,分别称取500 g(干质量)含水率为20%,含油率分别为5%、10%、15%、20%、25%、30%的油浸泥土置于热解炉内,通过自动调温系统维持炉内油浸泥土温度为300 °C,热脱附处理4 h后出料。待出料冷却后,测定固相含油率。

    为探究油浸泥土初始含水率对于油浸泥土热脱附处理效果和热脱附能耗的影响,分别称取500 g(干质量)含油率均为11.3%,含水率分别为5%、10%、20%、30%、40%、50%的油浸泥土置于热解炉内,通过自动调温系统维持炉内油浸泥土温度为300 °C,热脱附处理4 h后出料。待出料冷却后,测定固相含油率。

    油浸泥土含油率的测定采用萃取结合红外测油仪测定的方法[20]。称取10 g(干质量)样品于100 mL具塞磨口锥形瓶中,加入10 g无水硫酸钠,混合均匀,加入25 mL四氯化碳,盖好塞子轻轻摇匀。将锥形瓶置于50~60 °C水浴中,以50~60 r·min−1振荡8 h,然后静置1 h。将上清液倾倒于50 mL 容量瓶中,再向100 mL具塞磨口锥形瓶中加入20 mL四氯化碳,以50~60 r·min−1振荡8 h,然后静置1 h;再次将上清液倾倒于50 mL容量瓶中,定容,利用红外测油仪对上清液进行定量测定分析。

    采用 Excel 2010 和 Origin7.5软件对实验数据进行处理、数据分析、绘图。

    石油烃去除率根据式(1)进行计算。

    η=CtC0C0×100% (1)

    式中:η为石油烃去除率;C0为油浸泥土初始含油率;Ct为处理后油浸泥土含油率。

    在不同温度条件下,热脱附处理后样品的表观性状如图4所示。原始样品(如图1所示)由于含油率较高而呈黏结状态。油浸泥土的最大污染特征表现为土壤表里的贯通性,石油烃灌满于一定面积、深度土壤中几乎所有孔隙,堵塞绝大多数的土壤气孔;同时由于石油烃的黏稠性,将土壤中原本散状的土壤颗粒胶黏在一起,改变了土壤原有的结构特征。经过4 h热处理后,热脱附温度在250 °C以下的油浸泥土仍呈现黏结块状,从表观上来看仍还有一定量的石油烃。

    图 4  不同加热温度下处理后土壤的表观性状
    Figure 4.  Appearance properties of treated soil at different heating temperatures

    在不同加热温度下,经过热脱附处理后油浸泥土含油率的变化情况如图5所示。设置加热温度为200 °C,经过4 h处理后,油浸泥土含油率降至2.7%;在热脱附温度从200 °C升高至300 °C的过程中,含油率可持续降低至0.7%。在此温度区间内,升高加热温度,石油烃的去除率增大,处置后油浸泥土样品的含油率降低。当加热温度为300 °C时,含油率仅为0.7%,石油烃去除率可达到94%;继续升高温度(>300 °C),油浸泥土的含油率基本保持稳定(含油率保持在0.6%左右)。综合考虑热脱附效果及资源消耗2个因素,该油浸泥土的最佳热脱附处理温度为300 °C。

    图 5  不同温度条件下热脱附处理后油浸泥土样品的含油率
    Figure 5.  Petroleum content in treated soil at differentheating temperatures

    油浸泥土经过高温加热,土壤中的石油烃会热解吸挥发出来。并且当加热温度提高时,石油烃的C—C键和C—H键易断裂而发生裂解,石油烃链缩短,挥发性增强,从油浸泥土中解吸出来,实现油土分离[17]。本研究结果表明,升高加热温度可以促进油浸泥土中石油烃的去除。这是由于石油烃分子中的化学键在较高温度下更容易断裂,加热温度越高,产生裂解的石油烃种类越多。但当温度在300 °C以上时,油浸泥土含油率无显著变化,这说明在初始含油率保持不变的情况下,加热温度超过一定值时,对热脱附效果的影响有所减弱。王嫣云等[21]研究了不同热解终温对石油烃去除率的影响,结果表明,随着温度升高,热解残渣的含油率逐渐降低,但当热解温度为260 °C时,继续升温,对油浸泥土含油率则无明显影响。 LUNDIN等[22]的研究表明,高温热解有利于垃圾焚烧灰中多氯联苯和多氯联苯并呋喃的去除,热解温度从300 °C升高至500 °C时,热解后固态残渣污染物含量明显降低。热脱附温度的选择与污染物本身的特性有关。对于被苯等挥发性有机物污染的地块,温度上升至100 °C,就有较好的修复效果;而对于多环芳烃、石油烃等半挥发性有机物分子,通常需要升温至300~500 °C,才能达到相同的修复效果[16, 23]

    在不同加热时间下,油浸泥土热脱附处理后样品的表观性状如图6所示。加热时间(0~1 h)较短时,部分样品仍存在因石油烃的黏稠性而导致的黏结状态;经过4 h热脱附处理后,油浸泥土团聚现象消失,松散程度高,形成均匀细颗粒。

    图 6  不同加热时间下处理后土壤的表观性状
    Figure 6.  Appearance properties of treated soil at different heating times

    对于不同的加热时间,样品含油率检测结果如图7所示。在加热时间为0~0.5 h时,样品中石油烃去除率较明显,含油率从11.3%降至2.6%,这是因为样品中的轻质油挥发速度较快。从0.5 h继续延长加热时间至4 h的过程中,样品含油率从2.6%降低至0.7%;在4~6 h的加热过程中,样品含油率从0.66%降低至0.56%,这表明热脱附处理4 h后,继续延长加热时间,油浸泥土中的含油率基本保持不变,因此,在此时间范围内,加热时间对于样品石油烃的去除效果影响较小。张攀等[11]在进行热脱附处理硝基苯研究时也得到类似的结果,硝基苯的热脱附效率随着热脱附时间的延长而增加,但当加热时间超过30 min时,热脱附效率不再有显著增加。LAM等[24]发现,石油污染土壤含油率会随着加热时间的增加而明显降低,但加热时间超过一定值时,热脱附效果基本不再受加热时间的影响。加热时间对热脱附效果的影响也取决于热脱附温度,较低的热脱附温度需要较长的加热时间才能达到相同的处理效果[25]。综合考虑热脱附效果、资源消耗及运行成本等因素,该油浸泥土的最佳热脱附加热时间为4 h。

    图 7  不同加热时间条件下热脱附处理后油浸泥土样品的含油率
    Figure 7.  Petroleum content in treated soilat different heating times

    图8可知,在设置的6个油浸泥土初始含油率浓度梯度(5%、10%、15%、20%、25%、30%)下,随初始含油率的升高,处理后的油浸泥土固相含油率呈现升高趋势,但石油烃去除率基本保持不变。因此,在初始含油率为5%~25%时,热脱附处理后,油浸泥土固相含油率在2%以下;但当初始含油率升高至30%时,石油烃去除率为91%,处理后的样品含油率为2.75%,尚未达到《油气田含油污泥综合利用污染控制要求》(DB 65/T 3998-2017)的修复标准。

    图 8  不同初始含油率条件下热脱附处理后油浸泥土样品的含油率及去除率
    Figure 8.  Petroleum Content in treated soil and removal rate at different initial petroleum contents of soil

    王嫣云等[21]的研究表明,初始含油率对高温下石油烃去除率影响较小,且去除率于300 °C后保持稳定。这主要是由于初始含油率越高,轻质油含量越高,在较低温度下,有更多组分析出; 含油土壤其他组分在温度升至较高时才逐渐析出,此时饱和去除率已较高,因此,石油烃组分去除率不再随初始含油率变化而变化。王瑛等[26]的研究也表明,土壤中初始DDT含量对DDT的去除率没有显著影响。

    在加热温度为300 °C,加热时间为4 h的工艺条件下,初始含水率对供试油浸泥土(含油率11.3%)热脱附效果和天然气用量的影响结果如图9所示。在初始含水率为10%~50%时,随着含水率的升高,石油烃去除率降低,经处理后,油浸泥土固相含油率有所升高。在油浸泥土初始含水率超过40%时,处理后含油率不能达到修复标准。由此可见,油浸泥土初始含水率过高或过低,均会影响热脱附处理效果。高含水率不利于热脱附处理,高含水率的油浸泥土中的水分会以水膜形式存在于油浸泥土表层,在热脱附过程中水脱离会吸收部分热能,这使石油烃得到的能量降低,导致石油烃热脱附效率降低。此外,含水率过高会造成油浸泥土黏结使其通透性降低,堵塞石油烃扩散的通道,影响石油烃的热脱附效果[27]。油浸泥土含水率较低也不利于热脱附处理,油浸泥土中水分子作为极性分子相较于石油烃更容易与土壤上位点结合,占据更多的结合位点,含水率过低时,石油烃分子大量吸附于油浸泥土上,故导致热脱附处理效果较差[28]

    图 9  不同初始含水率条件下热脱附处理后油浸泥土样品的含油率及能耗
    Figure 9.  Petroleum content and energy consumption of treated soil at different moisture contents

    图9所示,天然气用量与油浸泥土初始含水率呈现反比关系,油浸泥土含水率影响热脱附能耗。初始含水率为50%时,天然气使用量为72 m3·h−1(标准状况);当初始含水率降至5%时,天然气使用量降低至20 m3·h−1 (标准状况),高含水率的油浸泥土热脱附过程的耗能是低含水率油浸泥土的3倍之多。虽然过高含水率的油浸泥土热容量高,热损失较慢,但较高的含水率不利于其中污染物的脱附,且易消耗大量的能源,工程成本较高。因此,油浸泥土初始含水率是热脱附技术应用的重点考虑因素之一[29]

    1) 热脱附技术可以用于油浸泥土中石油烃污染物的去除,且热脱附的处理效果与加热温度、加热时间、土壤的初始含油率及含水率有关。

    2) 利用本研究热解炉小试装置,对于初始含油率为11.3%、含水率为20%的油浸泥土,综合考虑样品热脱附处理效果、资源消耗及运行成本,加热温度300 °C、加热时间4 h为最佳工艺处理条件。在此最佳条件下,处理后的油浸泥土含油率仅为0.7%,符合油泥处置标准。

    3) 油浸泥土初始含水率和含油率是决定是否采用热脱附技术的关键因素。在一定温度和时间范围内,处理后油浸泥土的固相含油率随着初始含油率的升高而升高,不达标的风险增大。适中的初始含水率(10%~20%)有利于热脱附处理,当初始含水率超过40%时,热脱附效果较差且能耗较高。

  • 图 1  旋流扩散混合器与传统混合器模型图对比

    Figure 1.  Comparison of cyclone diffusion mixer and traditional mixer model

    图 2  扰流柱结构示意图

    Figure 2.  Structure schematic of turbulence column

    图 3  计算域及网格划分

    Figure 3.  Calculation domain and grid division

    图 4  120°和60°交叉板倾角下压力分布

    Figure 4.  Pressure distributions at the angels of 120° and 60° for the cross plates

    图 5  旋流扩散静态混合器与SK型混合器模拟结果

    Figure 5.  Simulation results of swirl diffusion static mixer and SK mixer

    图 6  SK型和新型静态混合器浓度分布

    Figure 6.  Concentration distributions of SK type and new type static mixers

    表 1  正交实验结果

    Table 1.  Results of orthogonal experiment

    序号(A)相邻蘑菇头间距/mm(B)相邻蘑菇头交错角/(°)(C)十字交叉板夹角/(°)不均匀系数压力损失系数/%
    1250600.05722.74
    2257.5900.09810.33
    325151200.1264.90
    4300900.10910.36
    5307.51200.1095.01
    63015600.06722.05
    74001200.0744.96
    8407.5600.02422.47
    94015900.09010.38
    序号(A)相邻蘑菇头间距/mm(B)相邻蘑菇头交错角/(°)(C)十字交叉板夹角/(°)不均匀系数压力损失系数/%
    1250600.05722.74
    2257.5900.09810.33
    325151200.1264.90
    4300900.10910.36
    5307.51200.1095.01
    63015600.06722.05
    74001200.0744.96
    8407.5600.02422.47
    94015900.09010.38
    下载: 导出CSV

    表 2  正交实验结果极差分析

    Table 2.  Range analysis of orthogonal experiment results

    因素不均匀系数压力损失系数
    k1k2k3Rk1k2k3R
    A0.0940.0950.0630.03212.65712.47312.6030.184
    B0.0800.0770.08940.01712.68712.60312.4430.244
    C0.0490.0990.1030.05422.42010.3574.95717.463
      注:不均匀系数的实验结果顺序为C>A>B;压力损失系数的实验结果顺序为C>B>A。
    因素不均匀系数压力损失系数
    k1k2k3Rk1k2k3R
    A0.0940.0950.0630.03212.65712.47312.6030.184
    B0.0800.0770.08940.01712.68712.60312.4430.244
    C0.0490.0990.1030.05422.42010.3574.95717.463
      注:不均匀系数的实验结果顺序为C>A>B;压力损失系数的实验结果顺序为C>B>A。
    下载: 导出CSV

    表 3  Topsis多指标模型评价结果

    Table 3.  Evaluation results of Topsis multi-index model

    序号不均匀系数压力损失系数Di+DiTopsis评价指数Ci
    10.7860.4790.3460.1450.296
    20.6340.7630.1850.2430.568
    30.5270.8880.2140.3390.613
    40.5900.7630.2070.2380.535
    50.5920.8850.1780.3390.656
    60.7490.4950.3380.1250.270
    70.7240.8860.1030.3550.775
    80.9090.4850.3340.2140.391
    90.6620.7620.1730.2470.587
    序号不均匀系数压力损失系数Di+DiTopsis评价指数Ci
    10.7860.4790.3460.1450.296
    20.6340.7630.1850.2430.568
    30.5270.8880.2140.3390.613
    40.5900.7630.2070.2380.535
    50.5920.8850.1780.3390.656
    60.7490.4950.3380.1250.270
    70.7240.8860.1030.3550.775
    80.9090.4850.3340.2140.391
    90.6620.7620.1730.2470.587
    下载: 导出CSV
  • [1] 宝露尔, 张海峰, 杨敏. 三卤甲烷和卤乙腈类消毒副产物在14个饮用水厂出水中的浓度水平及原水氯化/氯胺化中的生成势[J]. 环境工程学报, 2017, 11(11): 5843-5847. doi: 10.12030/j.cjee.201604003
    [2] 赵玉丽, 李杏放. 饮用水消毒副产物: 化学特征与毒性[J]. 环境化学, 2011, 30(1): 20-33.
    [3] 邱建华. 游泳池水消毒系统设计问题初探[J]. 科技经济市场, 2012(12): 70-72. doi: 10.3969/j.issn.1009-3788.2012.12.028
    [4] 方志林. SK型静态混合器的数值模拟[J]. 广东石油化工学院学报, 2013, 23(4): 84-86. doi: 10.3969/j.issn.2095-2562.2013.04.023
    [5] SOMAN S S, MADHURANTHAKAM C M R. Effects of internal geometry modifications on the dispersive and distributive mixing in static mixers[J]. Chemical Engineering and Processing, 2017, 122: 31-43. doi: 10.1016/j.cep.2017.10.001
    [6] 任新林, 梅毅, 冯梦黎, 等. SK静态混合器对工业磷酸脱砷的过程强化研究[J]. 化工学报, 2018, 69(S2): 218-225.
    [7] 吴雨, 张力钧, 宋忠俊. 静态混合器在石油化工中的应用[J]. 天然气与石油, 2014, 32(3): 23-26. doi: 10.3969/j.issn.1006-5539.2014.03.007
    [8] DIEZ A M, MOREIRA F C, MARINHO B A, et al. A step forward in heterogeneous photocatalysis: Process intensification by using a static mixer as catalyst support[J]. Chemical Engineering Journal, 2018, 343: 597-606. doi: 10.1016/j.cej.2018.03.041
    [9] 刘佳薇. 基于CFD的跑道池式微藻反应器中叶片型静态混合器的优化设计[D]. 武汉: 华中科技大学, 2019.
    [10] 陶雪峰. 低雷诺数湍流控制方法中混合阶段的试验与数值模拟研究[D]. 杭州: 杭州电子科技大学, 2019.
    [11] HIRSCHBERG S, KOUBEK R, MOSER F, et al. An improvement of the Sulzer SMXTM static mixer significantly reducing the pressure drop[J]. Chemical Engineering Research & Design, 2009, 87(4): 524-532.
    [12] HOSSEINI S M, RAZZAGHI K, SHAHRAKI F. Design and characterization of a low-pressure-drop static mixer[J]. AIChE Journal, 2019, 65(3): 1126-1133. doi: 10.1002/aic.16505
    [13] 张春成. SK型静态混合器污泥与药剂混合性能研究[D]. 长沙: 中南大学, 2014.
    [14] 胡坤, 李振北. ANSYS ICEM CFD工程实例详解[M]. 北京: 人民邮电出版社, 2014.
    [15] 王泽斌, 宋涛, 张春梅. SV型静态混合器数值模拟湍流模型选择[J]. 当代化工, 2016, 45(1): 166-168. doi: 10.3969/j.issn.1671-0460.2016.01.052
    [16] 王华然, 王尚, 李昀桥, 等. 臭氧在水中的溶解特性及其影响因素研究[J]. 中国消毒学杂志, 2009, 26(5): 481-483.
    [17] 王修纲. 静态混合环管反应器的流体力学研究及应用[D]. 天津: 天津大学, 2017.
    [18] 周锐. 基于FLUENT对管式静态混合器的优化设计[D]. 武汉: 华中科技大学, 2018.
    [19] 李婷. 基于熵权的TOPSIS模型在水污染风险评价中的应用[J]. 黑龙江水利科技, 2015, 43(7): 16-17. doi: 10.3969/j.issn.1007-7596.2015.07.006
    [20] YU W Z, TANG D S. Application of TOPSIS model based on vague set entropy in the evaluation of groundwater quality[J]. Advanced Materials Research, 2013, 712-715(1): 452-456.
  • 加载中
图( 6) 表( 3)
计量
  • 文章访问数:  3699
  • HTML全文浏览数:  3699
  • PDF下载数:  48
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-12-10
  • 录用日期:  2020-05-10
  • 刊出日期:  2020-11-10
胡珊, 毛澍洲, 邱光宇, 许小红, 张波. 应用于臭氧消毒系统的新型静态混合器结构设计[J]. 环境工程学报, 2020, 14(11): 3201-3207. doi: 10.12030/j.cjee.201912068
引用本文: 胡珊, 毛澍洲, 邱光宇, 许小红, 张波. 应用于臭氧消毒系统的新型静态混合器结构设计[J]. 环境工程学报, 2020, 14(11): 3201-3207. doi: 10.12030/j.cjee.201912068
HU Shan, MAO Shuzhou, QIU Guangyu, XU Xiaohong, ZHANG Bo. Structural design of new static mixer for ozone disinfection system[J]. Chinese Journal of Environmental Engineering, 2020, 14(11): 3201-3207. doi: 10.12030/j.cjee.201912068
Citation: HU Shan, MAO Shuzhou, QIU Guangyu, XU Xiaohong, ZHANG Bo. Structural design of new static mixer for ozone disinfection system[J]. Chinese Journal of Environmental Engineering, 2020, 14(11): 3201-3207. doi: 10.12030/j.cjee.201912068

应用于臭氧消毒系统的新型静态混合器结构设计

    通讯作者: 张波(1975—),男,博士,副教授。研究方向:水处理技术及环保工艺设备开发。E-mail:zhangbo@ujs.edu.cn
    作者简介: 胡珊(1996—),女,硕士研究生。研究方向:水处理技术。E-mail:1069828090@qq.com
  • 1. 江苏大学环境与安全工程学院,镇江 212013
  • 2. 江苏恒泰泳池科技股份有限公司,常州 213179
  • 3. 江苏大学常州工程研究院,常州 213164
基金项目:
江苏省研究生科研与实践创新计划项目(SJCX19_0567);常州市科技支撑项目(CE20185005)

摘要: 臭氧水消毒过程中,液液混合均匀性是关键。为提高臭氧水消毒过程的混合均匀性,优化设计了一种旋流扩散静态混合器。通过数值模拟得到旋流扩散静态混合器中臭氧水浓度及压力分布,以旋流扩散静态混合器前后的不均匀系数及压损系数作为评价指标,探究混合元件间距、夹角等因素对臭氧水混合均匀性及压力损失的影响,从而优化混合器结构设计,并与传统SK 型混合器进行比较。模拟结果表明:最优结构布局参数为相邻蘑菇头间距40 mm,交错角0°,十字交叉板夹角120°。在此情况下,不均匀系数由1.92降至0.074,压力损失系数为4.96%。新型混合器相比传统混合器具有更好的混合效果,且十字交叉板夹角大小对不均匀系数、压力损失的影响较大。

English Abstract

  • 臭氧消毒技术具有氧化性强、杀菌效率高、安全无卤代副产物的优势,可克服氯系消毒剂易产生余氯气味及三卤甲烷、卤乙酸等致癌副产物的缺点[1-2]。在管道分质供水系统、游泳池分流量臭氧消毒系统中,射流器或气液混合泵将臭氧溶解于旁路管道水中(循环水量的15%~25%),臭氧水进入主管道经静态混合器与未经臭氧处理的水混合,臭氧与水充分接触发挥消毒作用[3]。因此,流体混合的均匀程度对消毒效率提升起着重要作用。静态混合器与其他设备相比,具有流程简单、体积小、投资少、操作方便等优点[4-5],被广泛应用于石油、化工、环保过程不同流体的连续分散混合[6-8]

    国内外学者通过数值模拟方法对静态混合器的结构设计及流动特性进行了大量研究。刘佳薇[9]采用数值模拟优化设计光生物反应器中静态混合器的叶片结构,在不大幅提升系统功耗的前提下,改善了藻液混合效果。陶雪峰[10]研究了絮凝系统中SK型静态混合器的混合过程,对内部流场进行数值模拟,并与实验结果进行了对比分析。HIRSCHBERG等[11]研究发现在SMX元件之间增加适当间隙,可明显降低混合压降。HOSSEINI等[12]设计了低压降静态混合器,其摩擦系数小于现有SMV、KMX和挡板式静态混合器。然而,现有静态混合器仍存在以下缺点:混合单元结构复杂,流体通过时受到的剪切机械力大、能耗大[13];标准化的静态混合器内部元件规格确定,不易根据实际工况调整组合,混合均匀度有限。因此,新型静态混合器具有良好的设计开发及应用价值。

    针对现有静态混合器结构上的缺陷,本课题组设计了一种旋流扩散静态混合器。为优化其结构布局,运用3因素3水平正交实验,考察相邻蘑菇头间距、相邻蘑菇头交错角、交叉板夹角对不同液体混合效果及压力损失的综合影响。采用三维数值仿真研究混合元件前后的流场特性,优化结构参数,得到最优布局,并将其与传统SK型混合器作对比,以期为旋流扩散静态混合器设计提供参考。

  • 旋流扩散静态混合器的工作原理与传统混合器类似。不同的是,静态混合器只保留了交叉板与中心导流柱,板后区域通过增设蘑菇头状扰流柱对流动进行干扰,以提高液液混合效果。如图1所示,混合器前段的交叉板将流体分割,从两侧流入,达到分流目的。次流液体沿中心轴线进入流域,在导流柱前缘处被分散经过交叉板,加强与主流的混合作用。相比于传统混合器,旋流扩散静态混合器具有结构紧凑、混合效率高等优点。

  • 静态混合器内部多采用交叉板或扭曲叶片的形式。流体沿着扭曲叶片绕管中心轴线旋转前进,存在中心与周边流体掺混均匀度不高的情况,且流体易在壁面处形成附面层。

    旋流扩散静态混合器后段为数个蘑菇头元件均匀交错排列在管道内壁上。图2为一组混合元件,随着流速增大,混合器内部会产生激烈的湍流和剪切力,蘑菇头扰流柱打破了流体在壁面处形成的附面层。这种结构设计促进了流体被分散切割、反向旋转、径向混合的效果。

  • 本研究的重点为混合器内部流体掺混的稳态流动,故控制方程中不考虑非稳态项。数值计算中使用的控制方程是湍流状态下时均形式的Navier-Stokes方程,主要包括连续方程(式(1))、动量方程(式(2)),及Reynolds应力与湍流黏度μt和平均速度梯度关联的Boussinesq表达式(式(3))。

    式中:ρ¯uiuj为Reynolds应力或湍流应力。

    k-ε模型作为经典湍流模型之一,可提高计算速度及求解稳定性[14]。王泽斌等[15]在对SV型静态混合器进行模拟时,就湍流模型的选取指出,k-ε模型计算结果与实验吻合较好,是适合混合器流场仿真的湍流模型之一。因此,本研究选用标准k-ε模型进行数值模拟。该模型下湍动能及湍流耗散率的输运方程分别为式(4)和式(5)。

    式中:μt=ρCμk2εC1ε=1.44C2ε=1.92Cμ=0.09σk=1.0σε=1.3

  • 由于固壁及流固之间的换热对流动几乎没有影响,故计算只涉及流动区域。计算域及网格划分如图3所示。采用非结构网格进行划分,且对近壁面处网格进行局部加密,以保证壁面处y+在30~300范围内,以满足壁面函数处理的要求。交叉板及蘑菇头扰流柱区域为流动主要研究区域,因此,在网格划分过程中对该部分进行加密,以保证相对准确的流动计算。各模型平均网格量为2.87×106个,全局壁面网格符合y+在30~300范围内。

  • 数值计算在ANSYS Fluent上进行,选用稳态计算方式,求解器的相关设置分4个部分。

    1)边界条件设定。主流和次流入口均采用速度入口边界条件,分别设置为1和 3.5 m·s−1;主流进口为液态水,次流进口为臭氧水;臭氧在水中溶解度低[16],对液态水物性影响小,故臭氧水采用与液态水相同的物性参数,即密度998.2 kg·m−3、黏度0.001 Pa·s;出口采用自由出流边界(outflow),其余壁面均设为无滑移壁面。

    2)计算模型选择。湍流模型选择Standard 模型,壁面处采用壁面函数处理方式;开启组分运输模型表征2种液体的混合,组分数为2,混合物密度计算方式为体积加权混合。

    3)收敛判断准则。迭代过程中,监视模型中截面处的面积加权平均速度和平均质量分数,当其值保持不变且连续性残差、动量残差等各残差均小于10−3 时即认为收敛。

    4)数值计算方法。选择Fluent 中的压力求解器对数值模拟中涉及的各个方程进行离散,并对代数方程进行求解,以获得数值解;计算过程中的流体压力与流体速度耦合方法选择SIMPLE 算法进行计算,动量方程求解采用二阶迎风格式以保证精度。

    液液分布混合程度采用不均匀系数ψ定量计算,其定义[17-18]如式(6)~式(8)所示。

    式中:σ为截面质量分数分布方差;¯ϕ为截面质量分数算术平方根;N为截面上的节点数;ϕi为每个节点上臭氧水相的质量分数。

    混合段的压损用压力损失系数δ表示,其定义如式(9)所示。

    式中:P0为混合段进口静压,Pa;P1为混合段出口静压,Pa。

    在仿真后处理中,分别截取混合段前后2处横截面,并计算得到截面上的不均匀系数和平均压力,进而分析混合器的混合效果及压损。

  • 将相邻蘑菇头间距(A)、相邻蘑菇头交错角(B)、交叉板夹角(C)作为可变参数,采用相同的边界条件对正交实验表中9组模型开展数值模拟计算。正交实验结果如表1所示,极差分析如表2所示。

    用直观分析法对表1设计的9组正交实验结果进行分析。不均匀系数反映了混合器出口处不同因素及水平组合对水和臭氧水的混匀效果;压力损失系数反映了臭氧水经过十字交叉板、蘑菇头混合元件至混合器出口这一段的压力降低百分比。表2均值3行中k1k2k3的数据分别代表某一影响因素在某一水平下的系数均值。极差R越大说明该因素对该系数的影响越大。从表2极差结果可知,影响不均匀系数大小的各因素重要性顺序为十字交叉板夹角>相邻蘑菇头间距>相邻蘑菇头交错角;影响压力损失系数大小的各因素重要性顺序为十字交叉板夹角>相邻蘑菇头交错角>相邻蘑菇头间距。因此,交叉板的夹角对流动有较大影响。

    表1中压损系数可知,第1、6、8组压损较大,对应十字交叉板夹角均为60°;而3、5、7组压损较小,对应十字交叉板夹角均为120°。如图4所示,流体通过交叉板后流速增大,压力降低。流体在通过60°交叉板时,由于板夹角较小,板面的倾斜程度较120°低,对流动产生较大阻碍,部分流体在交叉板表面滞止,造成较大压损,且根据伯努利方程(式(10))得出,流速降低会引起静压的升高,故图4(b)整体静压高于图4(a)

    式中:Pstatic为流体静压;ρ为密度;v为流速;PTotal为流体总压;C为恒定常数。

    选用准确、客观的评价方法能为新型混合器的结构设计提供重要依据。本研究使用基于熵权的Topsis模型评价方法[19-20],即利用熵权法确定不均匀系数及压损系数2个指标的权重,并应用Topsis模型得到评价指数,对9组正交实验结果进行综合评价。

    通过Topsis多指标评价模型对各组结果进行分析,结果如表3所示。根据结果,确定第7组为最优组,其混合效果和压损都较为合理,其结构布局参数为相邻蘑菇头间距40 mm,交错角0°,十字交叉板夹角120°。

  • 图5表明各组模型计算结果之间的差异,其中第10组为传统SK型静态混合器的模拟结果。由图4可知,改进后的9组混合器,其充分发展后的不均匀系数均低于第10组,表明其自身具有很好的分散混合性能。混匀效果最好的一组相较于传统混合器其均匀程度提高了83.8%。但随着混合效果的提高,沿程压损出现一定程度的上升,增加水泵耗能。

    SK型静态混合器是使用较为广泛的一类静态混合器,由一系列相对于轴线扭转180°的螺旋板组成,不容易堵塞。与传统SK型静态混合器相比,旋流扩散静态混合器是在十字交叉切割片后,圆柱体内部一圈设置了类似于蘑菇头形状的混合单元。这种结构使液体在前进的过程中受到蘑菇头状扰流体的阻碍。由于液体流速不高,易在壁面处形成附面层,扰流柱的存在可以打破附面层,完善次流液体与主流液体的混合均匀程度。图6给出了2种混合器截面的臭氧水浓度分布,并截取了混合段前后两处横截面,分别对应图中黑色虚线位置。

    混合器(图6(a))出口横截面周边一圈浓度较中心高,是因为次流流体沿螺旋片绕流旋转,在离心力作用下中心与周边液体较难掺混。然而,由图6(b)可知,旋流扩散静态混合器出口截面浓度大小分布稳定且相同,几乎完全混合。由于液体在前进的过程中受到蘑菇头状扰流体的阻碍,水流不断地被分散、打乱、改变流向,从而获得良好的径向混合效果。

  • 1)旋流扩散静态混合器中交叉板夹角的角度对不均匀系数及压损的影响最为显著。夹角大,不均匀系数大,压损小;夹角小,不均匀系数小,压损大。通过Topsis建立的多指标评价模型,计算得到第7组参数最优,其结构布局为:相邻蘑菇头间距40 mm,交错角0°,交叉板夹角120°。在此情况下,不均匀系数为0.074,混合均匀程度提高96.15%,压力损失系数为4.96%。

    2)在相同且较宽的操作条件下,流体充分发展后的不均匀系数均达0.15以下,新型静态混合器的不均匀系数均低于传统SK型。这表明静态混合器自身具有较好的分布混合性能,且改进后的混合效果更好。

参考文献 (20)

返回顶部

目录

/

返回文章
返回