Loading [MathJax]/jax/output/HTML-CSS/jax.js

BAF曝气管道布置方式的数值模拟

武魏魏, 叶家盛, 郑俊, 张德伟, 黄河清. BAF曝气管道布置方式的数值模拟[J]. 环境工程学报, 2019, 13(11): 2668-2674. doi: 10.12030/j.cjee.201811187
引用本文: 武魏魏, 叶家盛, 郑俊, 张德伟, 黄河清. BAF曝气管道布置方式的数值模拟[J]. 环境工程学报, 2019, 13(11): 2668-2674. doi: 10.12030/j.cjee.201811187
WU Weiwei, YE Jiasheng, ZHENG Jun, ZHANG Dewei, HUANG Heqing. Numerical simulation on the layout of BAF aeration pipeline[J]. Chinese Journal of Environmental Engineering, 2019, 13(11): 2668-2674. doi: 10.12030/j.cjee.201811187
Citation: WU Weiwei, YE Jiasheng, ZHENG Jun, ZHANG Dewei, HUANG Heqing. Numerical simulation on the layout of BAF aeration pipeline[J]. Chinese Journal of Environmental Engineering, 2019, 13(11): 2668-2674. doi: 10.12030/j.cjee.201811187

BAF曝气管道布置方式的数值模拟

    作者简介: 武魏魏(1993—),男,硕士研究生。研究方向:水污染防治。E-mail:1943960036@qq.com
    通讯作者: 黄河清(1964—),男,博士,教授。研究方向:环境流体力学。E-mail:heqing@ahut.edu.cn
  • 中图分类号: X703

Numerical simulation on the layout of BAF aeration pipeline

    Corresponding author: HUANG Heqing, heqing@ahut.edu.cn
  • 摘要: 针对曝气生物滤池(BAF)曝气不均导致的处理效率低下的问题,基于Fluent软件,采用标准k-ε模型对BAF多种曝气管道的布置方式进行了不同维度的数值模拟研究。利用二维数值模拟,通过对单边枝状、中心枝状和环状3种曝气方式的流场分析,结果表明,单边枝状和中心枝状对应的曝气孔气量分布较均匀,而环状的4个弯角处由于相对低压造成曝气量少或无曝气;进一步就原型双边枝状曝气的三维数值模拟表明,各曝气口的曝气量都接近理论均值,标准偏差为7.67×10−5 m3·s−1,可以取得较好的均匀曝气效果。本研究可以有效地助力曝气生物滤池的选型设计,以较低的成本实现高效节能的目标。
  • 厌氧氨氧化工艺氮去除负荷高、无需外加有机碳源、污泥产量低、运行成本低,但因厌氧氨氧化菌生长速度慢、倍增时间长,且对环境条件变化较为敏感,使得厌氧氨氧化工艺启动耗时较长,这极大地限制了厌氧氨氧化技术的工程应用[1]。Dokhaven污水厂的厌氧氨氧化(anaerobic ammonia oxidation,ANAMMOX)工艺耗时3.5 a成功启动[2]。因此,厌氧氨氧化菌(anaerobic ammonia oxidizing bacteria,AnAOB)的高效富集、ANAMMOX工艺的快速启动及稳定运行引起研究者的广泛关注。

    ANAMMOX工艺的快速启动与接种污泥性质、反应器类型密切相关[3]。不同类型的好氧、厌氧污泥均可用于ANAMMOX工艺启动。好氧污泥虽菌群丰富,但其所含厌氧菌属较少,启动期会相应较长;厌氧反硝化污泥中含有的反硝化菌与AnAOB同属厌氧菌,可省去好氧转向厌氧环境时的污泥适应阶段,完成AnAOB的富集,可缩短工艺启动市场,但其沉降性能较差,启动过程中易出现污泥流失[4];而经厌氧消化后的剩余污泥具有良好沉降性能,且来源广泛、易于获取,污泥碳氮比含量低、高氨氮适应性强、厌氧菌含量高且与AnAOB的代谢基质相近,较适合作为ANAMMOX工艺启动的接种污泥。反应器类型直接影响污泥固体停留时间(solid retention time,SRT)、水流上升流速等,进而影响AnAOB富集、工艺启动速度、工艺运行稳定性等;序批式反应器(sequencing batch reactor,SBR)、序批式生物膜反应器(sequencing batch biofilm reactor,SBBR)、膜生物反应器(membrane bioreactor,MBR)、生物滤池均可成功启动ANAMMOX工艺[5-7]。上流式厌氧污泥床(upflow anaerobic sludge blanket,UASB)反应器作为第2代厌氧反应器的杰出代表,具有较好的污泥持留能力和基质传质效果,可为AnAOB的生长提供良好的环境,其在ANAMMOX工艺运行中的优势已逐步显现。WANG等[8]将厌氧氨氧化颗粒污泥接种于UASB反应器(22 L),经过178 d启动及稳定运行,氮容积负荷(nitrogen loading rate,NLR)和氮去除负荷(nitrogen removal rate,NRR)可高达8.25 kg ·(m3·d)−1(以N计)和6.93 kg·(m3·d)−1(以N计)。

    因此,本研究以厌氧消化污泥为接种污泥,以UASB反应器为反应装置启动ANAMMMOX工艺,研究了ANAMMOX-UASB反应器启动过程的生物特性,考察了ANAMMOX- UASB启动运行过程中的污泥表观形态、脱氢酶及胞外聚合物的变化,分析了启动过程中胞外聚合物结构、组成的特征,解析了功能菌群动态演替规律,以期探明ANAMMOX-UASB启动过程中的生物特性,为厌氧氨氧化工艺的快速启动与工程应用提供参考。

    所采用的实验装置如图1所示,UASB反应器呈圆柱结构,由有机玻璃制成,总容积为3.5 L,有效容积为1.4 L,其中反应区内径6 cm,高度50 cm。模拟废水由蠕动泵从UASB反应器底部连续泵入,从下至上依次经反应区、沉淀区后,最后经出水口排出;产生的气体经顶部排气孔排出。反应区部分采用黑布包裹以避免光照对AnAOB产生的不利影响。反应器工作温度为(35±1) ℃。

    图 1  实验装置示意图
    Figure 1.  Schematic diagram of experimental equipment

    实验中所用接种污泥取自青岛市某污水处理厂厌氧消化罐,污泥形态呈黑褐色絮状。混合液悬浮固体质量浓度(mixed liquid suspended solids,MLSS)约为47.39 g·L−1,混合液挥发性悬浮固体质量浓度(mixed liquor volatile suspended solids,MLVSS)约为22.77 g·L−1,MLVSS/MLSS为0.48。接种污泥量为1.2 L。

    实验废水采用人工模拟废水,NH+4-N和NO2-N分别由NH4Cl和NaNO2按需配置,其他主要成分为27 mg·L−1 KH2PO4、500 mg·L−1 NaHCO3、180 mg·L−1 CaCl2·2H2O、300 mg·L−1 MgSO4·7H2O。微量元素Ⅰ、Ⅱ添加量均为1 mL·L−1。微量元素Ⅰ的组成成分及含量为5 mg·L−1 EDTA、5 mg·L−1 FeSO4·7H2O、微量元素Ⅱ为15 mg·L−1 EDTA、0.99 mg·L−1 MnCl2·4H2O、0.25 mg·L−1 CuSO4·5H2O、0.43 mg·L−1 ZnSO4·7H2O、0.014 mg·L−1 H3BO4、0.19 mg·L−1 NiCl2·6H2O、0.22 mg·L−1 Na2MoO4·2H2O。废水pH为7.5~7.8。

    NH+4-N采用纳氏试剂分光光度法测定;NO2-N采用N-(1-萘基)乙二胺分光光度法测定;NO3-N采用紫外分光光度法测定;pH采用玻璃电极法测定;SS和VSS采用重量法[9]测定。脱氢酶(dehydrogenase activity,DHA)的提取与测定采用TTC还原法[10]。胞外聚合物(extracellular polymeric substance,EPS)提取与测定参照WANG等[11]的方法分层提取松散型胞外聚合物(loosely bound EPS,LB-EPS)与紧密型胞外聚合物(tightly bound EPS,TB-EPS)。采用硫酸-蒽酮比色法[12]测定多糖(polysaccharide,PS)含量,采用改进的Folin-酚试剂法[13]测定蛋白质(protein,PN)的含量。

    三维荧光光谱(three-dimensional excitation-emission matrix fluorescence spectroscopy,3D-EEM)采用荧光光度计(F-4 600, Hitachi, Japan)测定样品中LB-EPS和TB-EPS的三维荧光光谱[14]。激发波长(Ex)和发射波长(Em)分别为200~450 nm和240~550 nm,扫描间隔为5 nm,扫描速度采用1 200 nm·min−1。用超纯水作为空白样品校正水的拉曼散射。实验数据采用Origin绘图分析。

    聚合酶链式反应(polymerase chain reaction,PCR)及高通量测序:采用Power Soil DNA Isolation Kit按照操作步骤提取样品中DNA,使用1%(质量分数)的琼脂糖凝胶电泳检测DNA浓度和质量。采用细菌16S rRNA基因V3+V4扩增区域,引物为341F(CCTACGGGNGGCWGCAG)和805R(GACTACHVGGGTATCTAATCC),对DNA进行PCR扩增,扩增体系及扩增条件按彭广生等[15]方法进行。PCR扩增产物使用OMEGA胶回收纯化试剂盒纯化后4 ℃保存。样本交由北京诺禾致源生物信息科技有限公司进行DNA提取和测序。利用Illumina HiSeq高通量测序技术在HiSeq 2500系统进行测序。

    ANAMMOX-UASB反应器启动过程共持续250 d,整个过程分为菌体水解期(阶段Ⅰ)、活性迟滞期(阶段Ⅱ)、活性提高期(阶段Ⅲ)和稳定运行期(阶段Ⅳ和阶段Ⅴ)[16]。启动过程运行情况见表1

    表 1  ANAMMOX-UASB启动过程中各阶段运行情况
    Table 1.  Performance of each stage during ANAMMOX-UASB startup process
    阶段运行时间/dHRT/h进水N质量浓度/(mg·L−1)出水N质量浓度/(mg·L−1)TN去除率/%
    NH+4-NNO2-NNH+4-NNO2-NNO3-N
    1~13245050152.36±85.151.22±1.164.07±1.76−57.65±84.85
    14~532450504.68±7.927.91±6.313.4±1.7384.19±10.6
    54~1232450~30050~3000.17±0.400.37±0.4214.71±18.296.6±2.75
    124~16024~123003001.86±2.620.83±0.7756.21±0.7790.18±2.49
    161~25012~63003008.55±7.243.22±2.0267.54±6.986.78±1.72
     | Show Table
    DownLoad: CSV

    ANAMMOX-UASB反应器接种污泥为絮状黑褐色厌氧消化污泥(图2(a))。随着启动运行,污泥内有机物不断消耗、异养菌逐渐死亡,导致反应器污泥层高度逐渐下降。第53天,絮状污泥颜色逐渐转变为黄褐色,且部分区域污泥呈现团聚体状。自第100天起,反应器内出现小颗粒状污泥,且部分已呈现浅红色(图2(b)),结合此过程反应器脱氮性能可断定厌氧氨氧化反应已明显显现。此时,反应器内的污泥与郑平等[17]报导的厌氧氨氧化颗粒污泥形状一致,说明絮状污泥逐渐颗粒化,厌氧氨氧化污泥初步形成。颗粒污泥呈现红色是由于厌氧氨氧化菌体内含有丰富的细胞色素c,污泥发红程度可以反映厌氧氨氧化菌的富集程度,可以作为肉眼判断厌氧氨氧化活性的依据[18]。随着ANAMMOX-UASB反应器继续运行,絮状污泥颗粒化进一步加强,至第250天时,反应器内以不规则状的红色颗粒污泥和褐色絮状污泥为主(图2(c)),此时NRR维持在2 kg·(m3·d)−1左右,厌氧氨氧化菌已成功富集。

    图 2  ANAMMOX-UASB启动过程的污泥表观形态
    Figure 2.  Sludge apparent morphology during ANAMMOX-UASB startup process

    脱氢酶主要参与微生物降解有机物的过程[10],也可以参与某些细胞的合成。一氧化碳脱氢酶是乙酰辅酶A途径的关键酶。尽管厌氧氨氧化菌不以有机物为碳源,但仍有脱氢酶活性。ANAMMOX-UASB反应器启动过程中脱氢酶活性变化,可以反映出厌氧氨氧化工艺启动过程中污泥中异养菌的消长情况,结果见图3。由图3可知,因接种污泥为厌氧消化污泥,其异养微生物丰度较高,脱氢酶活性为3 909.51 μg·(h·g)−1(以TF计)。在第53天,污泥的脱氢酶活性降低至2 788.809 μg·(h·g)−1,这主要是由内源有机物量降低,异养微生物逐步衰亡、自溶等造成。随着启动过程的继续(53~250 d),厌氧氨氧化菌丰度逐步增加,异养菌含量进一步降低,脱氢酶活性大幅度降低,最终降至72.13 μg·(h·g)−1。LIN等[19]利用厌氧氨氧化污泥作为接种污泥启动厌氧氨氧化工艺时,启动前脱氢酶活性为684 μg·(h·g)−1,成功启动后脱氢酶活性降为252 μg·(h·g)−1。KIM等[20]的研究表明,异养菌脱氢酶参与了有机化合物的同化和异化反应,自养菌的脱氢酶主要参与有机化合物的同化,异养菌的脱氢酶活性远高于自养菌脱氢酶。在厌氧氨氧化工艺启动中,因进水基质无有机物,污泥中异养菌大量裂解死亡,导致脱氢酶活性大幅下降;厌氧氨氧化菌虽然不利用有机物进行分解代谢,但仍以CO2作为无机碳源,通过乙酰辅酶A途径合成自身细胞物质,故工艺在成功启动后仍具有一定的脱氢酶活性。

    图 3  ANAMMOX-UASB启动过程中的脱氢酶变化
    Figure 3.  Change in DHA activity during ANAMMOX-UASB startup process

    1) EPS含量。EPS是厌氧氨氧化污泥的重要组成部分,厌氧氨氧化细菌会嵌入到由细胞和EPS组成的聚集体中[21]。EPS具有流动性的双层结构,是由松散附着的外层(LB)和紧密黏附的内层(TB)组成,EPS的主要成分是蛋白质(PN)和多糖(PS)[22]。ANAMMOX-UASB反应器启动过程中EPS含量及组成见图4。由图4可知,第0~53天,EPS含量快速降低,由初期的34.91 mg·g−1(以VSS计)(0 d)下降为14.48 mg·g−1(53 d),PN/PS由7.88(0 d)升至11.58(53 d);EPS中TB/LB由1.19变为1.07。第53~123天,EPS含量略降至14.19 mg·g−1,PN/PS由11.58升至14.25;EPS中TB/LB为6.21。第123~160天,EPS含量逐渐上升至28.02 mg·g−1,PN/PS进一步升至18.59;EPS中TB/LB为4.15。第160~250天,EPS含量继续升至39.21 mg·g−1,PN/PS升至23.20;EPS中TB/LB变为3.12。EPS中PN含量变化大,PS含量变化相对较小,且在启动过程中PN/PS持续上升。启动过程中EPS均以TB-EPS为主;LB-EPS中LB-PN/LB平均值93.8%,TB-EPS中TB-PN/TB平均值为93.1%,各分层EPS中均以PN为主,表明EPS中蛋白质在ANAMMOX-UASB启动过程中起关键作用,这与GUO等[23]报道的研究结果一致。大量的PN产生可能来源于细胞死亡后的释放、细胞的分泌以及EPS中含有大量胞外酶所致。PN是疏水性基团,PS是亲水性基团,MA等[24]的研究表明,EPS中PN的含量与污泥的絮凝和沉降能力成正相关,较高的PN含量可以增强微生物聚集体的凝聚性。DONG等[25]的研究表明,PN中带正电荷的氨基可以中和羧基的负电荷,从而加剧污泥的絮凝。CHEN等[21]发现,PN/PS越大,污泥的沉降性能越好。在ANAMMOX-UASB启动过程中,PN/PS由启动前的7.88增至23.20,与絮状污泥逐渐颗粒化及其沉降性能逐渐增强的变化趋势一致。

    图 4  ANAMMOX-UASB启动过程中的EPS含量及组分变化
    Figure 4.  Changes of EPS contents and components during startup process of ANAMMOX-UASB

    2) EPS组分的荧光特性。污泥EPS中含有大量荧光特性的物质,可利用三维荧光分光光度法分析EPS组成特征。启动过程中第0、53、123天的LB-EPS和TB-EPS三维荧光光谱如图5所示。由图5可知,ANAMMOX-UASB在启动过程中共观察到5个荧光峰(A~E)。其中,荧光峰A(220~225 nm/330~340 nm)为芳香族蛋白类物质;荧光峰B(260~275 nm/335~340 nm)为色氨酸蛋白类物质;荧光峰C(270~275 nm/300~305 nm)与酪氨酸蛋白类物质有关;荧光峰D(220 nm/310~320 nm)与简单芳环蛋白类物质有关,荧光峰E(270 nm/425 nm)为富里酸类物质[26]

    图 5  启动过程中LB-EPS和TB-EPS的三维荧光光谱
    Figure 5.  Three-dimensional fluorescence spectra of LB-EPS and TB-EPS during startup

    荧光峰A、B在污泥启动过程的LB-EPS和TB-EPS中均存在,说明在启动过程中的EPS均以蛋白质为主;虽然LB-EPS和TB-EPS所处位置不同,但其组分和结构基本相同。在LB-EPS中仅存在荧光峰A和B,无其他荧光峰,且启动过程中LB-EPS的结构与组成并未发生较大变化。在第53天时TB-EPS中,除存在荧光峰A和B,还出现荧光峰D,说明此阶段有简单芳环蛋白类物质产生;第123天时,出现荧光峰C和E,荧光峰B消失,说明部分色氨酸蛋白类物质转化为酪氨酸蛋白类物质[27]。三维荧光结果表明,在ANAMMOX-UASB启动过程中污泥EPS中的TB-EPS的结构和组成均会发生变化。

    ANAMMOX-UASB启动过程中LB-EPS和TB-EPS中各峰光谱参数见表2。由表2可知,荧光峰会向更长的波长移动(称为红移),或向更短的波长移动(称为蓝移)。在LB-EPS中,第123天荧光峰A沿Ex方向蓝移5 nm,沿Em方向红移5 nm;荧光峰B沿Ex方向蓝移15 nm,沿Em方向红移5 nm。在TB-EPS中,第53天时荧光峰A沿Em方向红移10 nm,但至第123天时沿Em方向蓝移10 nm,与启动前相同;荧光峰B沿Em方向红移5 nm。上述研究结果表明,波长移动与EPS中各成分的结构变化密切相关,其中波长红移与含羰基取代基、羟基、烷氧基、氨基和羧基含量增加相关联,而波长蓝移则与芳香酯基减少,芳香环数下降,共轭键数的降低和羰基、羟基和氨基的减少直接相关。荧光峰波长及峰强度的变化说明ANAMMOX-UASB启动过程中污泥EPS中的组分、结构均发生变化。

    表 2  ANAMMOX-UASB启动过程中LB-EPS和TB-EPS的荧光光谱参数
    Table 2.  Fluorescence spectral parameters of LB-EPS and TB-EPS during startup process of ANAMMOX-UASB
    样品时间/d峰A峰B峰C峰D峰E
    Ex/Em强度Ex/Em强度Ex/Em强度Ex/Em强度Ex/Em强度
    LB0225/330274.5275/335175.1
    53225/330193.4275/335136.2
    123220/33551.0260/34065.3
    TB0220/330176.5275/335109.0
    53220/340296.4275/340210.0220/310260.0
    123220/33047.7270/30048.0220/31047.8270/42515.3
     | Show Table
    DownLoad: CSV

    1)微生物多样性和丰富度分析。ANAMMOX-UASB启动过程中的微生物多样性和丰富度如表3所示。由表3可知,所有样品的覆盖率均大于99%。这表明样品数据涵盖的物种充足,微生物群落丰富度和多样性结果具有较高的可靠性和真实性[28]。ACE和Chao1指数可以表征物种的丰富度,其值越大,所含物种越丰富。Shannon指数常用来评价生物群落组成复杂程度,其值越大,表明群落复杂程度越高[29]。Simpson指数更倾向于反应群落的均匀性,其值越大,表示优势菌群占总体生物菌群比例越大[30]。由表3可知,接种污泥的Chao1和ACE指数分别为594.16和583.74,第53天的Chao1和ACE指数分别上升至最大值,为695.08和682.31。这表明此时的群落中所含物种最丰富,结合启动中出水指标可得,此时反硝化菌、氨化菌、厌氧氨氧化菌共存。第250天分别下降到351.00和312.34,表明工艺启动成功后,微生物以自养厌氧氨氧化菌为主,群落丰富度有所下降。接种污泥的Shannon指数为5.19,第53天上升至最大值为6.06,表明此时群落复杂度最大;第250天降低至4.43,其原因是反应器成功启动后,厌氧氨氧化菌的丰度不断升高,其他微生物的丰度逐渐降低,导致微生物群落复杂度有所下降。接种污泥的Simpson指数为0.92,第53天上升至最大值0.96,表明此时反应器内优势菌群所占的比例也最大;第250天下降为0.87,较接种污泥略有降低,表明反应器启动成功后,优势菌群所占的比例有所下降。

    表 3  ANAMMOX-UASB启动过程中微生物多样性和丰富度
    Table 3.  Microbial diversity and abundance during startup process of ANAMMOX-UASB
    时间/dChao1ACEShannonSimpson覆盖率/%
    0594.16583.745.190.9299.5
    53695.08682.316.060.9699.5
    123343.23328.194.780.9099.8
    250351.00312.344.430.8799.7
     | Show Table
    DownLoad: CSV

    2)微生物门水平物种丰度分析。将微生物在门水平上丰度前10的物种绘制成物种相对丰度柱形累加图(图6)。主要为变形菌门(Proteobacteria)、绿弯菌门(Chloroflexi)、厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)、Latescibacteria菌门、Cloacimonetes菌门、浮霉菌门(Planctomycetes)、放线菌门(Actinobacteria)和疣微菌门(Verrucomicrobia)。其中浮霉菌门和变形菌门为脱氮功能菌。变形菌门在各个样品中占最大比例,在启动过程中丰度逐渐上升,由7.22%上升到54.40%。曹雁等[31]利用上流式厌氧过滤床(upflow blanket filter,UBF)反应器启动厌氧氨氧化工艺,启动成功后变形菌门增加了12%。杨开亮等[32]利用SBR接种活性污泥启动厌氧氨氧化工艺时,变形菌门在启动过程中丰度逐渐下降。上述研究结果不同,可能与反应器类型、接种污泥类型、运行参数及环境条件等不同有关。沈耀良等[33]在厌氧折流板(anaerobic baffled reactor,ABR)反应器中发现各隔室中变形菌门为主要菌种,其在脱氮方面占有重要地位。绿弯菌门是厌氧氨氧化系统中常见的伴生菌门,接种污泥的丰度为31.23%,第53天下降为20.32%,第123天反应器启动成功时丰度上升为28.52%。CHEN等[34]研究发现,绿弯菌门是兼性厌氧菌,在厌氧氨氧化污泥中起到支撑和骨架作用,可作为形成小颗粒污泥的框架粒子,也可作为颗粒污泥的载体,从而有利于厌氧氨氧化颗粒污泥的形成。MIAO等[35]的研究表明,厚壁菌门对反硝化有重要作用,厚壁菌门丰度由28.27%下降为5.55%,说明ANAMMOX-UASB启动过程中异养反硝化菌因环境条件改变大量裂解死亡。浮霉菌门为主要的自养脱氮功能菌群,厌氧氨氧化菌归属此菌门;接种污泥浮霉菌门丰度仅为0.34%,第123天启动成功时其丰度增至1.89%,表明反应器内富集了一定量的厌氧氨氧化菌。其他研究者也发现,厌氧氨氧化工艺启动过程中变形菌门丰度始终高于浮霉菌门丰度。如朱彤等[36]发现,在厌氧氨氧化反应器启动成功后,虽然变形菌门丰度有所下降,但仍高于浮霉菌门。

    图 6  ANAMMOX-UASB启动过程中微生物门水平组成
    Figure 6.  Composition of microorganisms at the phylum level during startup process of ANAMMOX-UASB

    3)微生物属水平物种丰度分析。由于厌氧氨氧化菌属于浮霉菌门,因此,对浮霉菌门进行属水平的分析(图7)。浮霉菌门有unidentified_PhycisphaeraeCandidatus_AnammoximicrobiumCandidatus_BrocadiaThermoguttaPirellulaRhodopirellula。在已知的6种厌氧氨氧化菌属中,启动过程中存在2种厌氧氨氧化菌属,分别为Candidatus AnammoximicrobiumCandidatus Brocadia。接种污泥中Candidatus Anammoximicrobium的丰度为40.94%,第53天下降到30.77%,第123天下降为0。杨瑞丽等[29]进行厌氧氨氧化工艺启动时,观察到Candidatus Anammoximicrobium丰度随NLR增大而逐渐降低,第172天下降到0,其可能的原因是Candidatus Anammoximicrobium生长速率较慢、亚硝酸盐亲和力较低以及对水质变化的抵抗力较弱。第123天,Candidatus Brocadia是反应器内唯一的厌氧氨氧化菌,其丰度为12.15%,第250天,下降到11.63%。曹雁等[31]利用UBF反应器培养厌氧氨氧化细菌时,Candidatus Brocadia的丰度由0.01%增加到1.00%。VAN DER STAR等[37]发现Candidatus Brocadia适合存在于高含氮浓度废水中。本实验中,进水氨氮质量浓度由50 mg·L−1提高到300 mg·L−1,属于高含氮废水,更适合Candidatus Brocadia生长,故Candidatus Brocadia启动过程中可得到有效富集,成为优势菌属。有研究[38]表明,由于生存环境的不同,厌氧氨氧化菌的群落结构存在差异,在稳定的生长环境中,通常只有1个属种的厌氧氨氧化菌占优势,这与本研究结果类似。

    图 7  ANAMMOX-UASB启动过程中浮霉菌门属水平组成
    Figure 7.  Composition of microorganisms of Planctomycetes at the genus level during startup process of ANAMMOX-UASB

    1) ANAMMOX-UASB启动过程中絮状污泥逐渐颗粒化,最终以不规则状的红色颗粒污泥和褐色絮状污泥为主。微生物的脱氢酶活性由3 909.51 μg·(h·g)−1持续降至72.13 μg·(h·g)−1。EPS主要以TB-EPS为主,占比为75.7%;且LB-EPS和TB-EPS中PN占比稳步增大,分别增至99.6%和94.7%,PN在厌氧氨氧化工艺启动过程中起关键作用。

    2)在ANAMMOX-UASB启动过程中,Chao1、ACE、Shannon和Simpson指数均呈现先升后降的趋势,启动过程中污泥的优势菌门为变形菌门(Proteobacteria)、绿弯菌门(Chloroflexi)、厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)、浮霉菌门(Planctomycetes)。浮霉菌门中Candidatus Anammoximicrobium丰度逐渐降低直至消失,而Candidatus Brocadia丰度最终增至12.15%。

  • 图 1  中心枝状及单边枝状曝气方式的速度分布

    Figure 1.  Velocity distribution of center dendritic and unilaterally dendritic pipe systems

    图 2  环状曝气方式的速度分布

    Figure 2.  Velocity distribution of ring-shaped pipe system

    图 3  原型曝气管道的三维数值模拟几何模型及模拟进气后的压力分布

    Figure 3.  Geometric model and pressure distribution of the 3-D numerical simulation of prototype pipe system

    图 4  最近、最远立管与上部主干管交界处及下部曝气管交界处的流场流线及速度云图

    Figure 4.  Streamlines and velocity nephograms of the junctions between the closest andfarthest stand pipes and main pipe or lower aeration pipes

    表 1  数值模拟模型参数

    Table 1.  Model parameters of numerical simulation

    重力加速度/(m·s−2)空气系数湍流模型常数时间步长/s
    ρ/(kg·m−3)µ/(kg·(m·s)−1)CµC1εC2ε
    9.81.2251.79×10−50.091.441.920.5
    重力加速度/(m·s−2)空气系数湍流模型常数时间步长/s
    ρ/(kg·m−3)µ/(kg·(m·s)−1)CµC1εC2ε
    9.81.2251.79×10−50.091.441.920.5
    下载: 导出CSV

    表 2  立管及其对应的中间位置曝气管的3个曝气孔的曝气量

    Table 2.  Flow rate of the three aeration holes located on the standing pipes and their corresponding center aeration pipes m3·s−1

    序号立管流量曝气管曝气孔流量
    最近居中最远
    1号0.008 210.000 980.000 840.000 82
    M1号0.008 240.001 000.000 850.000 80
    M2号0.008 280.001 000.000 860.000 84
    M3号0.008 290.001 000.000 830.000 81
    M4号0.008 310.000 990.000 840.000 83
    2号0.008 340.000 980.000 850.000 83
    序号立管流量曝气管曝气孔流量
    最近居中最远
    1号0.008 210.000 980.000 840.000 82
    M1号0.008 240.001 000.000 850.000 80
    M2号0.008 280.001 000.000 860.000 84
    M3号0.008 290.001 000.000 830.000 81
    M4号0.008 310.000 990.000 840.000 83
    2号0.008 340.000 980.000 850.000 83
    下载: 导出CSV
  • [1] 张杰, 陈秀荣. 曝气生物滤池反冲洗的特性[J]. 环境科学, 2003, 24(5): 86-91. doi: 10.3321/j.issn:0250-3301.2003.05.016
    [2] 张玉平. 两段上向流曝气生物滤池(TUBAF)处理城市污水试验研究[D]. 重庆: 重庆大学, 2006.
    [3] 陈康. 曝气生物滤池在城市污水处理工程中的应用[D]. 广州: 华南理工大学, 2013.
    [4] 廖腾毅. 基于生物膜分层理论的曝气生物滤池Fluent模拟研究[D]. 长沙: 湖南大学, 2014.
    [5] 王宏, 张林军, 宫秀学, 等. 上向流与下向流BAF串联处理生活污水的特性研究[J]. 中国给水排水, 2010, 26(7): 62-64.
    [6] 郑俊, 程寒飞, 王晓焱. 上流式曝气生物滤池工艺处理生活污水[J]. 中国给水排水, 2001, 17(1): 31-33.
    [7] 黄宇萍. 曝气生物滤池设计与施工的若干问题[J]. 中国给水排水, 2005, 21(8): 62-65. doi: 10.3321/j.issn:1000-4602.2005.08.017
    [8] 张诗华, 郑俊. 曝气生物滤池(BAF)及其设计[J]. 市政技术, 2007, 25(4): 270-272. doi: 10.3969/j.issn.1009-7767.2007.04.007
    [9] 谢观体, 陆少鸣. 上向流曝气生物滤池新冲洗方式的研究[J]. 中国给水排水, 2010, 26(19): 92-94.
    [10] 邱立平, 王广伟, 张守彬, 等. 上向流曝气生物滤池反冲洗实验研究[J]. 环境工程学报, 2011, 5(7): 1522-1526.
    [11] 杨婧雯, 江霜英, 刘昭丽. 上向流曝气生物滤池的低温启动试验[J]. 水处理技术, 2011, 37(1): 112-115.
    [12] 李微, 魏炜, 张传彬. 曝气量对间歇曝气/曝气生物滤池效能影响[J]. 环境科学与技术, 2010, 33(s1): 104-106.
    [13] 杨春娣. 曝气生物滤池中曝气方式对污水处理效率的影响研究[D]. 西安: 西安理工大学, 2008.
    [14] 刘玉玲, 白戈, 邵世鹏, 等. 曝气池曝气管布置方式对流速分布影响的数值模拟研究[J]. 水力发电学报, 2016, 35(7): 84-90. doi: 10.11660/slfdxb.20160709
    [15] 黄河清. 环境流体力学[M]. 合肥: 合肥工业大学出版社, 2013.
    [16] SAJJADI B. Solid-liquid mixing analysis in stirred vessels[J]. Reviews in Chemical Engineering, 2015, 31(2): 119-147.
    [17] 段中喆. ANSYS FLUENT流体分析与工程实例[M]. 北京: 电子工业出版社, 2015.
    [18] 唐家鹏. FLUENT 14.0超级学习手册[M]. 北京: 人民邮电出版社, 2013.
    [19] 张维佳. 水力学[M]. 北京: 中国建筑工业出版社, 2010.
  • 加载中
图( 4) 表( 2)
计量
  • 文章访问数:  3791
  • HTML全文浏览数:  3791
  • PDF下载数:  61
  • 施引文献:  0
出版历程
  • 收稿日期:  2018-11-29
  • 录用日期:  2019-06-05
  • 刊出日期:  2019-11-15
武魏魏, 叶家盛, 郑俊, 张德伟, 黄河清. BAF曝气管道布置方式的数值模拟[J]. 环境工程学报, 2019, 13(11): 2668-2674. doi: 10.12030/j.cjee.201811187
引用本文: 武魏魏, 叶家盛, 郑俊, 张德伟, 黄河清. BAF曝气管道布置方式的数值模拟[J]. 环境工程学报, 2019, 13(11): 2668-2674. doi: 10.12030/j.cjee.201811187
WU Weiwei, YE Jiasheng, ZHENG Jun, ZHANG Dewei, HUANG Heqing. Numerical simulation on the layout of BAF aeration pipeline[J]. Chinese Journal of Environmental Engineering, 2019, 13(11): 2668-2674. doi: 10.12030/j.cjee.201811187
Citation: WU Weiwei, YE Jiasheng, ZHENG Jun, ZHANG Dewei, HUANG Heqing. Numerical simulation on the layout of BAF aeration pipeline[J]. Chinese Journal of Environmental Engineering, 2019, 13(11): 2668-2674. doi: 10.12030/j.cjee.201811187

BAF曝气管道布置方式的数值模拟

    通讯作者: 黄河清(1964—),男,博士,教授。研究方向:环境流体力学。E-mail:heqing@ahut.edu.cn
    作者简介: 武魏魏(1993—),男,硕士研究生。研究方向:水污染防治。E-mail:1943960036@qq.com
  • 1. 安徽工业大学能源与环境学院,马鞍山 243032
  • 2. 安徽华骐环保科技股份有限公司,马鞍山 243032

摘要: 针对曝气生物滤池(BAF)曝气不均导致的处理效率低下的问题,基于Fluent软件,采用标准k-ε模型对BAF多种曝气管道的布置方式进行了不同维度的数值模拟研究。利用二维数值模拟,通过对单边枝状、中心枝状和环状3种曝气方式的流场分析,结果表明,单边枝状和中心枝状对应的曝气孔气量分布较均匀,而环状的4个弯角处由于相对低压造成曝气量少或无曝气;进一步就原型双边枝状曝气的三维数值模拟表明,各曝气口的曝气量都接近理论均值,标准偏差为7.67×10−5 m3·s−1,可以取得较好的均匀曝气效果。本研究可以有效地助力曝气生物滤池的选型设计,以较低的成本实现高效节能的目标。

English Abstract

  • 曝气生物滤池(biological aerated filter, BAF)是将生物氧化和截留吸附结合的新型好氧污水处理工艺,具有流程简单、处理效率高等优点[1]。按照其进水方式可分为上向流和下向流[2],其中,上向流滤池不仅布水和布气均匀、不易堵塞,而且防止了滤料的流失[3],因而应用较广泛[4-6]。此外,BAF池型一般为长方体,这便于施工和曝气管道布置,故很多设计研究中都选用该池型[7-8]。当污水流经滤料表面时,水中的溶解氧也向生物膜内扩散,污水中的有机物在滤料表面微生物的生化作用下得到有效降解,最终转化成CO2和H2O等代谢产物[9-11]。因此,良好的曝气可有效地提高对污染物的去除效果。本研究旨在通过数值模拟探究出较好的曝气方式及曝气效果。

    近年来,BAF曝气方面的相关研究多以实验和数值模拟为主。在BAF反冲洗的研究方面,张杰等[1]通过BAF反冲洗特性实验研究出脉冲气冲和连续水冲的组合比气水同时反洗要好;李微等[12]采用间歇曝气/曝气生物滤池工艺探究脱氮除磷的效果,通过定量控制一级及二级曝气量达到了较好的脱氮除磷效果;杨春娣[13]采用模型实验和数值模拟结合的方法,通过控制水深纵横比及曝气孔径,使气液两相流场分布均匀、混合充分;刘玉玲等[14]对推流式曝气池进行了数值模拟研究,通过对流场结构的分析探索出提高氧气利用率、减少能耗的方法。

    目前,有关BAF曝气管道布置方式的直接研究较少,采用数值模拟研究方法也很少。本研究以某污水处理厂采用前置反硝化曝气生物滤池的曝气管道系统为对象,首先进行二维数值模拟,探索中心枝状、单边枝状和环状曝气方式的优劣,然后通过三维数值模拟对拟采用的基于二维数值模拟所确定的较好曝气方式的曝气均匀性及效率进行定量确认。本研究运用数值模拟方法,将定性和定量地表达曝气系统的流场分布情况,可快速、准确地找出较好的曝气方式,从而为曝气生物滤池的选型设计提供参考。

  • 守恒方程包括质量守恒方程和动量守恒方程[15],如式(1)和式(2)所示。

    式中:ρ为流体密度,kg·m−3t为流动时间,s;xyz分别为直角坐标系的3个互相垂直方向的距离,m;uvw分别为xyz方向的速度,m·s−1

    式中:ρ为流体密度,kg·m−3µµt分别为流体分子动力黏度和涡动力黏度,kg·(m·s)−1g为重力加速度,m·s−2p为雷诺平均压力,Pa。

  • 本研究采用适用范围广、计算经济且被广泛应用于工业流场模拟的湍流标准k-ε模型,通过求解湍流动能k[16]和耗散率ε[16]方程,从而对涡黏度µt[16]进行求解。湍流动能k对应的输运方程见式(3)。耗散率ε的输运方程见式(4)。湍动黏度µt的方程见式(5)。

    式中:Gk为平均速度梯度引起的湍动能k的产生项;YM为可压湍流中脉动扩张的贡献;σkσε分别为与湍动能k和耗散率ε对应的普朗特数,取值为1.0和1.3;SkSε分别为定义的湍动能源项和湍流耗散源项;Cµ为经验常数,取值为0.09;C1ε为经验常数,取值为1.44;C2ε为经验常数,取值为1.92。

  • 在对连续性方程、动量方程求解时,都必须先确定求解初始条件和边界条件,而边界条件的设定对求解结果的影响十分关键[17-18]。本研究根据入口流量定义了曝气干管的速度入口,出口则根据曝气孔的设置深度定义为压力出口,管壁设定为给定粗糙度的固体界面。本研究采用专业模拟软件Fluent对曝气管道进行数值模拟求解。采用的数值求解方法为有限体积法,其具有很好的守恒性,可以克服泰勒公式展开出现离散的缺点且对网格的适应性很好。其中:对时变项,采用二阶精度的隐式格式的三时间段法;对对流项,采用二阶精度的迎风格式;对扩散项,采用二阶精度的中心差分格式。

    研究对应的数值模拟过程中用到的基本参数如表1所示。

  • 二维数值模拟研究了中心枝状、单边枝状和环状的3种曝气方式。中心枝状是指气体从中间干管位置进入,干管两侧支管对称分布(见图1(a));单边枝状是指气体从干管进入,支管布置在干管一侧(见图1(b));环状是指气体从干管进入,干管为环状,支管均匀分布在干管上(见图2(a))。这3种不同的曝气管道系统的干管、支管管径相同,分别为100 mm及20 mm;干管的入流速度根据所需曝气量取为8 m·s−1,出口根据水厂的水下管道设置深度定为50 000 Pa的压力出口;曝气管道上孔口按照等间距设置成圆形,孔口直径为2 mm。

    模拟所得中心及单边枝状曝气方式的二维速度场分别如图1(a)图1(b)所示。可见气体自干管的中间位置进入管道后,经支管到达曝气孔,几乎均匀地从各个曝气孔流出。

    环状曝气方式的管道及局部流场模拟结果如图2所示。图2(a)为以水厂原型的管道系统气体流动的二维速度分布。气体由右侧中间的干管入口进入,图2(a)下面b、c处分别为环状干管的弯角处,图2(b)上部d为位于环状系统中间对称面的代表部分。图2(b)弯管处的曝气量和干管直管处的图2(d)相比较少;图2(c)的干管弯管处附近的曝气孔几乎无曝气。进一步观察干管弯管处的气体流动可见,其靠近连接支管处的气体流动较快,由伯努利方程[19]可知,该处的压力小,产生了相对的负压,致使环状管网角部出气量少或者没有检测到出气量。模拟结果表明,应该避免环状曝气管网的设计,以提高滤池角部的曝气均匀性。

  • 根据二维模拟的结果,水厂拟采用如图3(a)所示的类似中心枝状的曝气管道系统。与图1(a)模型中的中心支状相比:二者相同点在于均以进气干管为对称轴,以保证干管左右两侧进气量均分;不同点在于中心支状的结构分布更均匀些,但实际安装、运行费用等方面存在不足。三维数值模拟对其曝气效果进行了模拟确认。根据其对称性,为提高模拟效益,取如图3(b)所示的一侧模拟即可。模拟的主要关注点是曝气量是否达到各个角落,其次是各个支管的曝气是否均匀。研究设计基本参数如下:水平干管、支管(立管)及下部水平曝气管的管径分别为150、50和20 mm;立管长为5.3 m,曝气管长为3 m。总曝气量为6 m3·min−1,根据对称性,模拟时半边取3 m3·min−1;曝气孔在水下5 m,取50 000 Pa为模拟的气体出口压力。由图3(b)所示的整体压力分布可以看出,管道压力自左侧气体入口至曝气管压力逐级减小,这说明曝气过程是基本均匀且为压力驱动的。

    研究分别选择离进气口最近的1号立管及最远的2号立管进行对比分析,若二者差别不大,那么基本可判定整体供气是均匀的。依次观察分析干管与立管交界处、立管与支管交界处以及曝气管的流场分布情况。

  • 与干管进气口最近的1号立管、最远的2号立管与上部主供气干管(图4(a)图4(b))及下部曝气管(图4(c)图4(d))交界处的流场如图4所示。由图4可见,在接近主进气口的1号立管(图4(a))和最远的2号立管(图4(b))的流场是有很大不同的。1号立管的进气口在主干管的下游侧进气较多,而上游侧较少;2号立管的进气口则相对均匀对称,中间多,两侧较少。但在立管的下部与曝气管交界处的流场(图4(c)图4(d)),二者却非常近似,接近于稳定管道流中的流场分布,即中间大,四周小。这可能是因为立管足够长,在流至曝气管前,使管道内的气体的流动达到了理想的稳定分布。从速度云图上看,速度的大小也是基本相同的,流量应比较接近,研究对曝气管的流量作进一步的定量检验。

  • 通过Fluent软件处理,可以检查数值模拟的各管道的流量及曝气孔的曝气量,定量地分析曝气均匀性。所有6根立管的进气量以及和它们相连接的中间曝气管的3个曝气孔的曝气量如表2所示。

    表2可以看出,各立管及曝气孔的曝气量分布基本均匀,接近理论平均值。其中,各立管流量随着离进气口越远而渐渐增加,理论均值为0.008 33 m3·s−1,但标准偏差仅为4.72×10−5 m3·s−1;而曝气孔的曝气量则是离干管越远越小的趋势,但总体差别不大,理论均值为0.009 3 m3·s−1,标准偏差为7.67×10−5 m3·s−1。各立管、曝气孔间的气量相对误差较小,可认为曝气基本达到了均衡,存在的较小差别可能为实际流动的沿途压损造成的,也有可能与湍流模型、网格质量有关,下一步可通过实测来检验各立管、曝气孔间的气量相对误差。

    综合以上流场及流量结果分析,可以认为该模型下的曝气达到了在各曝气管及曝气孔的流量的基本均衡。

  • 1)数值模拟研究为某曝气生物滤池的曝气管道系统进行了科学的选型优化。通过二维数值模拟得出:环状曝气的4个弯道处的曝气极少;而单边枝状和中心枝状的曝气方式均匀性相对较好。

    2)对原型双边枝状曝气方式三维的数值模拟显示,各曝气口的曝气量都接近均值,标准偏差为7.67×10−5 m3·s−1,故可以认为基本达到了均匀曝气。本研究可有效地助力曝气生物滤池的选型设计,以较低的成本实现高效节能的目标。

参考文献 (19)

返回顶部

目录

/

返回文章
返回