温度、相对湿度、排风量是脱水机房恶臭逸散的主要影响因素,温度过高、排风量不足都会造成脱水机房内恶臭气体不断聚集,恶臭气体主要成分为H2S、NH3,臭阈值低、毒性大[1-3]。若不妥善治理,将会严重影响工作人员的身体健康,危害较大[4-5]。为保证地下工作环境的安全,设计的排风量普遍过大,污水厂通风系统的投入和运行费非常高,造成能源浪费[6]。
污水处理厂的除臭系统属于末端治理:对于已产生的恶臭气体先进行密闭、收集,再通过化学、生物等方法对其进行治理。治理环节分为臭源密闭处理、恶臭气体收集输送及生化处理[7]。目前,大部分研究主要在恶臭气体的生化处理方面,通过优化通风除臭系统的研究还比较少[8]。影响内环境H2S、NH3浓度场的因素很复杂[9],其中,温度、湿度和排风量能够通过通风除臭系统进行控制。排风系统中除臭工艺的阻力非常大,风量过大造成排风系统能耗较高[7, 10-11]。H2S和NH3逸散本质是浓度差引起的质量扩散,研究表明,扩散系数受温度的影响比较大,温度越高,扩散系数越大[12]。这些因素对脱水车间H2S和NH3浓度的影响效应不同,根据权重效应合理选择通风系统控制目标能够达到节能的目的。因此,确定H2S和NH3的主要影响因素、权重及分布逸散规律,对地下污水处理厂恶臭气体的控制和节能具有重要意义。
数值仿真因其可视化、简单性和适应性受到研究者和设计者的青睐[13-14]。RNG k-ε湍流模型和组分运输模型被广泛应用[15-16]。本文采用这两种模型对地下污水处理厂污泥脱水车间不同温度、湿度和排风量下的H2S、NH3浓度分布仿真模拟,结合正交试验进行极差和方差分析,确定各因素影响权重,进一步分析主要影响因素与脱水机房H2S、NH3浓度和通风能耗的关系,通过综合评价指标对排风量进行优化,确定出污泥脱水机房合理的通风量。
以广州某地下污水厂污泥脱水机房为原型研究对象,如图 1所示,其长、宽、高分别为2 250、1 950、6 300 mm。由于整个脱水机房位于地下,温度相对比较稳定,忽略壁面热传递。脱水过程中会产生大量恶臭气体,虽脱水机进行了密封处理,但存储污泥的3个料斗仓裸露在内环境,不断地向周围逸散NH3和H2S,料斗口的长、宽、高分别为750、750、2 300 mm。为控制脱水机房内污染物浓度,排风系统以一定的排风量进行排风,防止恶臭气体的聚集,2个排风口位置在靠近料斗口的壁面侧,其长、宽分别为300、300 mm,中心高度2 000 mm。
六面体网格在计算精度、计算速度上都优于非结构化网格,故采用六面体结构网格对计算域进行离散。考虑到门、料斗口和排风口处扰动较大,为了提高计算收敛性,获取准确的区域参数对其进行局部加密,最终网格数量为251万,如图 2所示。
为简化计算做如下假设:室内空气为不可压缩流体;稳态计算各边界参数不变;脱水机房温度与临室温度相同,不考虑壁面热量传递,即壁面为绝热壁面;忽略机房内设备产生的热量。
在考虑能量守恒、动量守恒和质量守恒的基础上还需保证组分守恒,其通用控制方程可表示为
式中:φ为通用变量u、v、w、T、τp;Γφ为广义扩散系数;Sφ为广义源项。
选用RNG k-ε湍流模型以提高计算精度,其流动方程为
式中:表示湍流动能;表示湍流耗散项。
常数c1、c2、ak、aε具体数值详见参考文献[17]。
影响脱水车间H2S和NH3浓度分布的因素有温度、湿度和排风量,地下污水厂全年运行数据显示,脱水机房内温度范围10~33 ℃,相对湿度60%~90%,风机是变频运行,最大风量对应的排风口风压为-25 Pa,最小风量对应的风压为-15 Pa。考虑二阶交叉因素共有6个,如果每个因素取2个水平进行仿真模拟,要做26=64次试验,需要花费相当多的时间。正交试验法可通过少量实验次数对多因素进行极差分析、方差分析,并获得有价值的科学规律[18-19]。本文采用L8(27)正交试验,每个因素取最大值和最小值2水平,共需要试验8次。表 1是正交试验的各因素水平值。
按照L8(27)正交表进行仿真试验,取料斗口周围15个点的仿真结果对H2S、NH3逸散规律进行研究,如图 3所示。
8次试验结果记为y1,y2,…,y8,Kij为第j列因素第i水平所对应的H2S、NH3浓度值之和,κij=Kij,R=κijmax-κijmin,R值越大,说明该因素对试验指标的影响越大。表 2是点8处的H2S和NH3浓度的极差计算结果,从结果看出,A×C列对应R值最大,RH2S=9.4,RNH3=6.18。图 4是15个点各因素极差均值,从图中可以看出,温度和排风压力交互效应对应的极差最大,其值R=0.93,反映出温度和排风压力交互效应对H2S和NH3的浓度影响最大,其次为单因素排风压力C、温度A、相对湿度B,极差值相差不大,分别为0.89、0.71、0.65,交互因素A×B、B×C最小。因素影响权重从大到小依次为:A×C>C>A>B>A×B>B×C。
MSj是单因素j变化引起的差异,MSe是试验误差离差平方和,为提高检验的灵敏度,当MSj<2MSe将该因素并入误差。两水平正交试验,各因素离差平方和为
式中:n为试验总次数,n=8;a为每个水平试验次数,a=4。
表 3是点8处H2S和NH3浓度的方差计算,计算结果FA×C>F0.05(1, 4),表明温度和排气压力交互效应A×C对结果有显著影响。
图 5为各因素的均值极差和显著性大小,从图中可以看出,15个点的因素显著性分析中,因素A×C显著性p<0.05,占57%,其他因素小于13%。结果反映出,在这些因素中,温度和排气压力交互效应A×C对H2S和NH3浓度影响显著(p<0.05)。原因在于H2S和NH3在逸散和抽排过程中,温度和排风量存在交互效应,如图 6所示。在15 Pa排风压力作用下,温度由10 ℃升高到33 ℃,H2S和NH3浓度分别升高0.04、0.16 mg/m3,排风压力升高到25 Pa,温度由10 ℃升高到33 ℃,H2S和NH3浓度不同程度降低。不同排风压力下,温度对H2S和NH3浓度影响效果相反。
对广州某地下污水厂污泥脱水机房进行实地测量,测量时间为9:30到16:30。室外晴天,温度15.9 ℃,相对湿度59.7%,排风口静压-10 Pa。3台污泥脱水机有2台运行。在污泥脱水车间取3个高度分别为0.9、2.3、3.9 m,每个高度取9个测点坐标(x, y)依次为:
Y1(-6.6, -6.6), Y2(0, -6.6), Y3(6.6, -6.6)
Y4(-6.6, 0), Y5(0, 0), Y6(6.6, 0),
Y7(-6.6, 6.6), Y8(0, 6.6), Y9(6.6, 6.6)。
对各点的NH3和H2S的浓度进行测量。图 7为测量数据与模拟数据的对比图,由于实际情况的复杂性,实测数据与模拟数据存在一定的误差,H2S浓度平均相对误差为9.70%,NH3浓度平均相对误差为12.78%。对比结果:在允许误差条件下,仿真实验与实际测量基本一致。
为了确定温度和排风量对H2S和NH3浓度和能耗的影响,对温度和排风量进一步细化分析,工况情况如表 4。
以人员呼吸z=1.7 m水平高度进行分析,图 8是部分工况下H2S的浓度变化情况,从图中可以看出,排风量的变化或温度的变化都会对H2S浓度分布造成影响。
图 9反映出换气次数2.4次条件下,温度从15~35 ℃变化过程中,不同H2S浓度值对应区域的变化规律。15 ℃时面积比最大区域对应H2S浓度范围为0.019~0.033 mg/m3,同理, 20、25、30、35 ℃时对应H2S浓度值分别是0.013、0.026、0.039~0.046、0.039~0.052 mg/m3。从结果能够看出,温度由15 ℃升高到35 ℃,污染水平出现先降低后升高的趋势,在20 ℃时最低,与35 ℃相比降低了2.5倍。H2S浓度高于0.052 mg/m3的区域随着温度的升高不断扩大,面积比由20 ℃时的5%增大到35 ℃时的34%,高污染区域面积提高了4.8倍。可见,室外温度20 ℃是脱水机房恶臭污染控制最适宜的温度,污水厂脱水机房在满足室内环境要求的条件下可适当降低除臭风量以降低运行能耗。
图 10是25 ℃时不同排风量下H2S浓度与对应区域面积的关系图,最大面积对应的H2S浓度代表着该排风量下恶臭的水平。图中数据显示,换气次数由1.7次增加到4.7次,H2S最高浓度值出现不断降低的趋势,由0.072 mg/m3逐渐降到0.026 mg/m3。用标准方差来反映H2S浓度分布的均匀程度,换气次数由1.7次增加到4.7次,浓度方差值也不断下降,依次为70.9、69.5、63.2、40.8、39.4、30.6和30.6。从上述结果看出,排风量对地下污水厂脱水机房恶臭浓度有影响,排风量增大恶臭的排除量也增加,高浓度区域逐渐降低,并且浓度分布趋于均匀。因此,增大排风量能够降低H2S的浓度水平,抑制浓度的聚集。
研究表明,污染物浓度与换气次数存在一定的关系[20]。
由此可得:
式中:y为污染物浓度;y初为污染物初始浓度;y平衡为污染物稳定浓度;τ为换气指数,换气次数倒数;σ为污染物释放强度。
从公式可以看出,换气次数越大,脱水机房内NH3和H2S浓度y(t)越低,但恶臭浓度变化率|y′平衡|在不断降低。图 11为排风口、料斗口和内部点(-5、-2、1.7)处NH3随换气次数的变化规律,从图中曲线斜率可以看出,随着换气次数增大,NH3浓度变化率不断减小,换气次数超过2.7次,换气次数的增加对H2S和NH3浓度变化作用不大,与上述结论一致。
排风量增加虽然能够降低脱水机房内恶臭浓度,但是,同时风口风速随之增大,必定会造成能耗增加。能耗正比于风管阻力,单位管长的摩擦阻力损失由式(5)计算[21]。
式中:pm为单位管长摩擦阻力损失,Pa;d风管当量直径,m;v为风速,m/s;λ为摩擦阻力系数,由莫迪公式(6)求得。
脱水机房的通风能耗为
式中:p为全压,p=pm+1/2ρv2;G为通风量,m3/h;η1、η2为风机效率和机械效率,分别取值0.75和0.98。根据广州某地下污水厂的实际情况,L取50 m,K取0.18 mm。图 12中,换气次数与风口风速、能耗的关系显示,换气次数由1.7增加的4.7,风口风速由4.1 m/s增大到11.2 m/s,使得能耗由0.29 kW增加到8.16 kW,增加了27倍,单位风量能耗增加了6.1倍,随着通风量增加能耗增加显著。
设计规范《城镇污水处理厂污泥处理技术规程》要求污泥脱水机房换气次数不低于6次/h[22]。但研究表明,地下污水厂脱水机房换气次数小于6次条件下,NH3浓度降到0.09 mg/m3,H2S浓度降到0.012 mg/m3,已经低于《恶臭污染物排放标准》(GB 14554—93)中厂界一级排放标准。单位能耗对H2S和NH3作用分别由0.077、0.39 mg/(m3·kW)降为0.001、0.007 mg/(m3·kW),如图 13所示。通过继续增加排风量的方式改善空气品质显然是不可取的。
为了反映脱水机房恶臭污染水平和能耗水平,参考《实用供热空调设计手册》中IAQ污染物评价方法,提出综合污染水平指标及综合能耗指标。
式中:C为脱水机房不同换气次数下平均浓度;S是《恶臭污染物排放标准》(GB 14554—93)中厂界一级排放标准值。Ec是实际排风能耗,EH2S、ENH3是达到厂界一级排放标准值S需要的能耗。
结合污染水平和能耗水平,引入综合评价指标作为脱水机房通风除臭性能统一量化评价指标。
从式(10)可以看出,I值越小,污染程度越低,E越小,能耗水平越低,K值越小,脱水机房通风除臭性能越好。
根据式(8)~式(10)可得不同换气次数下通风量与K的关系,如图 14所示。从图中可以看出,在2.7换气次数,K最小为4.5,对应的H2S和NH3浓度约为0.022、0.16 mg/m3,排风口风速6.8 m/s,能耗1.5 kW。换气次数低于2.7次,脱水机房的浓度偏高,高于2.7次,风速过大,造成能耗过高,对于污水处理量10万t/d的脱水机房,排风量设置为2.7次换气次数最为合理。
评价指标方程化能够更加方便直观地判断脱水机房排风量是否合理。根据上述结果,对综合评价指标K进行非线性拟合,拟合模型选择指数函数,拟合曲线为,拟合结果为
式中:y表示综合评价指标K;x表示换气次数。
1) 通过正交仿真试验确定了温度和排风量的交互作用对污泥脱水车间H2S和NH3浓度的影响最大,且影响作用显著(p<0.05)。
2) 室外温度20 ℃时,脱水机房NH3和H2S浓度最低,有利于通风节能;换气次数超过2.7次后,通风量增加对H2S和NH3浓度变化不明显;综合评价指标K计算结果显示,对于污水处理量10万t/d的脱水机房,排风量为2.7次换气次数综合性能最好。
3) 过大的排风量对进一步改善空气品质的作用效果并不明显,反而使能耗急剧增加。因此,要降低综合评价指标K,单纯依靠机械排风是远远不够,主动控制脱水机房的恶臭气体还可以从气流组织的控制入手。