多因素对地下污水厂H2S与NH3浓度的影响

2018-08-31杨鹏刘雪峰门玉葵刘承东

土木与环境工程学报 2018年5期

杨鹏，刘雪峰，门玉葵，刘承东

(1.华南理工大学电力学院，广州 510641； 2.广东省市政工程设计研究总院,广州 510060)

温度、相对湿度、排风量是脱水机房恶臭逸散的主要影响因素，温度过高、排风量不足都会造成脱水机房内恶臭气体不断聚集，恶臭气体主要成分为H2S、NH3，臭阈值低、毒性大[1-3]。若不妥善治理，将会严重影响工作人员的身体健康，危害较大[4-5]。为保证地下工作环境的安全，设计的排风量普遍过大，污水厂通风系统的投入和运行费非常高，造成能源浪费[6]。

污水处理厂的除臭系统属于末端治理：对于已产生的恶臭气体先进行密闭、收集，再通过化学、生物等方法对其进行治理。治理环节分为臭源密闭处理、恶臭气体收集输送及生化处理[7]。目前，大部分研究主要在恶臭气体的生化处理方面，通过优化通风除臭系统的研究还比较少[8]。影响内环境H2S、NH3浓度场的因素很复杂[9]，其中，温度、湿度和排风量能够通过通风除臭系统进行控制。排风系统中除臭工艺的阻力非常大，风量过大造成排风系统能耗较高[7,10-11]。H2S和NH3逸散本质是浓度差引起的质量扩散，研究表明，扩散系数受温度的影响比较大，温度越高，扩散系数越大[12]。这些因素对脱水车间H2S和NH3浓度的影响效应不同，根据权重效应合理选择通风系统控制目标能够达到节能的目的。因此，确定H2S和NH3的主要影响因素、权重及分布逸散规律，对地下污水处理厂恶臭气体的控制和节能具有重要意义。

数值仿真因其可视化、简单性和适应性受到研究者和设计者的青睐[13-14]。RNGk-ε湍流模型和组分运输模型被广泛应用[15-16]。本文采用这两种模型对地下污水处理厂污泥脱水车间不同温度、湿度和排风量下的H2S、NH3浓度分布仿真模拟，结合正交试验进行极差和方差分析，确定各因素影响权重，进一步分析主要影响因素与脱水机房H2S、NH3浓度和通风能耗的关系，通过综合评价指标对排风量进行优化，确定出污泥脱水机房合理的通风量。

1 脱水机房计算模型

1.1 几何模型

以广州某地下污水厂污泥脱水机房为原型研究对象，如图1所示，其长、宽、高分别为2 250、1 950、6 300 mm。由于整个脱水机房位于地下，温度相对比较稳定，忽略壁面热传递。脱水过程中会产生大量恶臭气体，虽脱水机进行了密封处理，但存储污泥的3个料斗仓裸露在内环境，不断地向周围逸散NH3和H2S，料斗口的长、宽、高分别为750、750、2 300 mm。为控制脱水机房内污染物浓度，排风系统以一定的排风量进行排风，防止恶臭气体的聚集，2个排风口位置在靠近料斗口的壁面侧，其长、宽分别为300、300 mm，中心高度2 000 mm。

图1 广州某地下污水厂脱水机房3维图Fig.1 3-D model of Guangzhou underground sewage

六面体网格在计算精度、计算速度上都优于非结构化网格，故采用六面体结构网格对计算域进行离散。考虑到门、料斗口和排风口处扰动较大，为了提高计算收敛性，获取准确的区域参数对其进行局部加密，最终网格数量为251万，如图2所示。

图2 脱水机房六面体结构化网格Fig.2 Dewatering room hexahedral grid

1.2 模型简化

为简化计算做如下假设：室内空气为不可压缩流体；稳态计算各边界参数不变；脱水机房温度与临室温度相同，不考虑壁面热量传递，即壁面为绝热壁面；忽略机房内设备产生的热量。

1.3 数学模型

在考虑能量守恒、动量守恒和质量守恒的基础上还需保证组分守恒，其通用控制方程可表示为

(1)

式中：φ为通用变量u、v、w、T、τp；Γφ为广义扩散系数；Sφ为广义源项。

选用RNGk-ε湍流模型以提高计算精度，其流动方程为

(2)

式中：表示湍流动能；表示湍流耗散项。

(3)

常数c1、c2、ak、aε具体数值详见参考文献[17]。

2 多因素正交仿真实验

2.1 实验方案

影响脱水车间H2S和NH3浓度分布的因素有温度、湿度和排风量，地下污水厂全年运行数据显示，脱水机房内温度范围10～33 ℃，相对湿度60%～90%，风机是变频运行，最大风量对应的排风口风压为-25 Pa，最小风量对应的风压为-15 Pa。考虑二阶交叉因素共有6个，如果每个因素取2个水平进行仿真模拟，要做26=64次试验，需要花费相当多的时间。正交试验法可通过少量实验次数对多因素进行极差分析、方差分析，并获得有价值的科学规律[18-19]。本文采用L8(27)正交试验，每个因素取最大值和最小值2水平，共需要试验8次。表1是正交试验的各因素水平值。

表1 正交试验各因素水平值Table 1 Levels of factors in orthogonal test

注：排风量的驱动力是风机提供的静压，本次模拟的排风口边界条件为压力出边界。

2.2 实验结果与分析

2.2.1 极差结果分析按照L8(27)正交表进行仿真试验，取料斗口周围15个点的仿真结果对H2S、NH3逸散规律进行研究，如图3所示。

图3 仿真结果分析点分布Fig.3 The distribution of simulation analysis

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。

图4 各因素均值极差Fig.4 Mean value of extreme difference

因素试验号1234567A(温度)B(相对湿度)A×BC(排风压力)A×CB×C误差指标/(mg·m-3)H2SNH3111111114.6616.37211122224.1414.55312211223.5312.41

续表2

2.2.2 方差结果分析MSj是单因素j变化引起的差异，MSe是试验误差离差平方和，为提高检验的灵敏度，当MSj<2MSe将该因素并入误差。两水平正交试验，各因素离差平方和为

(4)

式中：n为试验总次数，n=8；a为每个水平试验次数，a=4。

表3是点8处H2S和NH3浓度的方差计算，计算结果FA×C>F0.05(1,4)，表明温度和排气压力交互效应A×C对结果有显著影响。

表3 正交试验方差分析Table 3 Variance analysis of orthogonal test

图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浓度影响效果相反。

图5 各因素显著性Fig.5 The significance of each

图6 温度和排气压力交互作用Fig.6 Interaction between temperature and

2.3 实验验证

对广州某地下污水厂污泥脱水机房进行实地测量，测量时间为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%。对比结果：在允许误差条件下，仿真实验与实际测量基本一致。

图7 实验测量与模拟值得对比图Fig.7 Comparison of experimental and simulated

3 多因素影响效应

3.1 工况设定

为了确定温度和排风量对H2S和NH3浓度和能耗的影响，对温度和排风量进一步细化分析，工况情况如表4。

表4 参数工况值Table 4 Value of various conditions

以人员呼吸z=1.7 m水平高度进行分析，图8是部分工况下H2S的浓度变化情况，从图中可以看出，排风量的变化或温度的变化都会对H2S浓度分布造成影响。

图8 不同工况下H2S浓度Fig.8 Changes of H2S concentration under

3.2 温度影响

图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 ℃是脱水机房恶臭污染控制最适宜的温度，污水厂脱水机房在满足室内环境要求的条件下可适当降低除臭风量以降低运行能耗。

图9 不同温度下H2S浓度值对应区域面积比Fig.9 Area ratio of H2S concentration at

3.3 排风量影响

图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的浓度水平，抑制浓度的聚集。

图10 25 ℃下不同H2S浓度对应区域面积Fig.10 Area of H2S concentration at 25

研究表明，污染物浓度与换气次数存在一定的关系[20]。

t→∞y平衡=τσ

由此可得：

式中：y为污染物浓度；y初为污染物初始浓度；y平衡为污染物稳定浓度；τ为换气指数，换气次数倒数；σ为污染物释放强度。

图11 不同换气次数下H2S和NH3浓度Fig.11 H2S and NH3 concentration under different

排风量增加虽然能够降低脱水机房内恶臭浓度，但是，同时风口风速随之增大，必定会造成能耗增加。能耗正比于风管阻力，单位管长的摩擦阻力损失由式(5)计算[21]。

(5)

式中：pm为单位管长摩擦阻力损失，Pa；d风管当量直径，m；v为风速，m/s；λ为摩擦阻力系数，由莫迪公式(6)求得。

(6)

脱水机房的通风能耗为

(7)

式中：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倍，随着通风量增加能耗增加显著。

图12 不同换气次数下的通风强度与通风能耗Fig.12 Ventilation strength and ventilation energy consumption under different ventilation

图13 单位能耗与H2S和NH3浓度关系Fig.13 Relationship between unit energy consumption and H2S and NH3

设计规范《城镇污水处理厂污泥处理技术规程》要求污泥脱水机房换气次数不低于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所示。通过继续增加排风量的方式改善空气品质显然是不可取的。

3.4 综合评价指标

为了反映脱水机房恶臭污染水平和能耗水平，参考《实用供热空调设计手册》中IAQ污染物评价方法，提出综合污染水平指标及综合能耗指标。

(8)

(9)

式中:C为脱水机房不同换气次数下平均浓度；S是《恶臭污染物排放标准》(GB 14554—93)中厂界一级排放标准值。Ec是实际排风能耗，EH2S、ENH3是达到厂界一级排放标准值S需要的能耗。

结合污染水平和能耗水平，引入综合评价指标作为脱水机房通风除臭性能统一量化评价指标。

K=exp(I)+exp(E)

(10)

从式(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次换气次数最为合理。