活性污泥2号模型的应用与校正

1、活性污泥2号模型的应用与校正    摘要：采用国际水质协会提出的活性污泥2号模型，以MATLAB为工具编制计算软件，采用该软件对常州市污水厂进行计算，在输入典型的动力学和化学计量参数值后，计算结果表明出水CODcr指标大部分得到较好的模拟，但NH3N、TP和TN出水指标模拟值与实际值有差异。通过对参数灵敏度的分析，在考虑温度影响的条件下，进行模型校正，模型校正后，预测精度得到显著提高。 关键词：城市污水 活性污泥法 数学模型  Application and Calibration of Activated Sludge Model No. 2

2、Abstract: The activated sludge model No.2 given by IAWQ and the calculating software developed by using MATLAB as the tool were used in calculating Changzhou Sewage Treatment Plant. With the input of typical ki-netic and stoichiometric parameters, the results of the calculation showed that most of t

3、he CODcr, figures of the outlet water had been well simulated while there had been differences between the simulation values and the actual values of the NH3-N, TP and TN figures of the outlet water. Having analyzed the sensitivity of the parametersand taken into consideration the effect of temperat

4、ure, the model was then calibrated. After the calibration of themodel, the accuracy of simulation was significantly improved.Key words: municipal sewage; sewage treatment; activated sludge method; mathematical model  国际水质协会（IAWQ）于1995年推出了活性污泥 2号模型1（Activated Sludge Modd No.2，ASM2），ASM2是活性污泥1号模型

5、2(Activated Sludge Model No.1，ASM1）的扩展，并沿用了该模型的概念，它包含碳有机物氧化、脱氮和生物除磷处理过程，共有19种组分、19种反应、22个化学计量系数及42个动力学参数。ASM2无论是在污水厂的设计、运管理、改造，还是在废水处理技术的研究和开发方面均有重要的使用价值，在国外已得到了成功广泛的使用3。我国在这方面起步很慢，本研究对ASM2在我国城市污水厂的实用性进行了验证。模型所用的污水进出水数据来自常州城北污水厂，污水厂采用A2O工艺，设计流量为5×104m3d，生物反应池的水力停留时间为11.5h，污泥负荷为0.15kgBOD5（kgMLSS

6、.d），污泥的质量浓度为3500mgL，固体平均停留时间为20d，污泥回流比100，厌氧池体积为2900m3，缺氧池体积为5800m3，曝气池体积为14500m3。由于该污水厂具有除磷脱氮功能，而ASM1没有考虑除磷过程，所以采用ASM2进行模拟研究。污水处理流程的简化表示如图1所示： 1 模拟计算结果此稳态计算程序利用Matlab4-6软件为平台，并根据ASM2对常州污水厂运行条件进行简化。假定该厂的供气量可以满足好氧池中微生物生化反应所需要的氧气量，即假定曝气池中溶解氧组分的取值恒定在2.0mgL，饱和溶解氧的质量浓度为8.637mgL（20），系统的操作温度恒定，二沉池不考虑微生物的物质

7、代谢活动，仅起固液分离作用且无活性污泥的积累。在这基础上，将2000年112月进水组分输人模型中，可以得到出水的CODcr，NH3N，TN，TP的稳态结果，模拟结果见图2-图5。 从图2一图5中可以看到，CODcr，NH3N，TN，TP模拟值与污水厂的实测值基本相符。CODcr的最大相对误差36.5%发生在12月份，差值为16.11mgL，其它11个月份相对误差在30之内，其中1，2，3，5月的相对误差控制在10之内，因而较好的模拟了CODcr值。NH3N，TN和TP指标因为出水中浓度很低，模拟有一定的困难，但模拟结果表明模拟值与实测值在同一个数量级别对于NH3N而言，模拟曲线比较平缓，这是因

8、为没有考虑温度对动力学参数的影响，因而模拟值与实测值差别较大。TP的模拟值与实测值的变化趋势相似，曲线很好的拟合了实际出水中TP的变化。偏差的存在有多方面的原因，一是模型假设的部分条件与实际污水厂运行不完全吻合，另外污水厂化验分析不可避免存在一定误差。尽管有一定的偏差，但模拟结果得到了良好的描述，因而利用活性污泥2号模型对城市污水厂进行模拟是可行的。 2 参数灵敏度和温度分析 21 参数灵敏度分析在输入IAWQ给出的模型动力学和化学计量学 参数的典型值7后，模拟结果表明出水指标大部分得到了较好的模拟，但部分出水的模拟值与实际值有 较大的误差，所以需要对一些参数进行校核，这里以 灵敏度分析为基础

9、，将模型中的动力学参数值都增加 其初始值的10，从而计算出水的CODcr，TN，TP，NN3N的灵敏度，计算结果见表1表4。 由灵敏度计算结果表1表4分析可知，大部 分参数的灵敏度都比较低（小于1），因此采用I- AWQ的推荐参数值进行模拟是合理的。但仍有些 参数对出水水质影响较大（灵敏度大于1）。对 CODcr影响较大的动力学参数有4个：Kx，bPAO， 异养菌KNH3N，bH。表1  水解过程动力学参数灵敏度分析 参数CODcrNH3-NTNTP水解速率常数Kh00.0330.0330.067氧的饱和抑制系数Ko20.50.500.5颗粒性COD的饱和系数Kx1111缺氧水解速率

10、降低修正因子NO300.1670.1670.333硝酸盐的饱和抑制系数KNO3-0.20.20.40.2厌氧水解速率降低修正因子fe0132表2异氧菌过程动力学参数灵敏度分析 参数CODcrNH3-NTNTP基于基质的最大生长速率H0.0830.01700.017基于发酵基质的生长饱和系数KF0.150.02500.05氨氮的饱和系数KNH3-N2202磷的饱和系数KP010010碱度的饱和系数KALK0101乙酸的饱和系数KA0.0250.02500.05发酵的最大速率qfe0.0330.03300.033溶菌速率常数bH1.50.252.250.25反硝化的速率降低修正固子NO3-0.12

11、50.53.3750.125表3  聚磷菌过程动力学参数灵敏度分析 参数CODcrNH3-NTNTPPHA贮存的速率常数qPHA00.0330.0330.033聚磷酸盐的饱和系数Kpp010010PP贮存的速率常数qpp0.0670.1330.0670.4Xpp/XPAO的最大比率Kmax00.29400.588PHA的饱和系数KPHA010020聚磷菌最大生长速率PAO0000.5聚磷菌溶菌速率常数bPAO10.510.5XPP的分解速率常数bPP00.50.50.5XPHA的分解速率常数bPHA0.50.500.5表4  自氧菌和沉淀再溶解过程动力学参数灵敏度分析 参数

12、CODcrNH3-NTNTP自养菌最大生长速率0.11.30.70.2自养菌衰减速率bAUT0.679.334.670.67自养菌氨氮的饱和系数KNH3-N010.50.1自养菌碱度的饱和系数KALK00.60.20.2自养菌氧的饱和系数KO20.20.80.40.2磷沉淀的速率系数KPRE00.100.1再溶解速率系数KRED0.1670.16700.167对NH3N影响较大的动力学参数有 10个：分别是KX，fe，异养菌KNH3-N，Kp，AUT，bAUT，自养菌KH3N，Kpp，KPHA，异养菌KALK。对TN影响较大的动力学参数有6个：分别是KX，fe，bH，bPAO，bAUT，反硝化

13、NO3。 对TP影响较大的动力学参数有7个：分别是Kx，fe，异养菌KNH3-N，Kp，异养菌KALK，Kpp，KPHA。22 温度影响分析一切生化反应都可以用Arrehenius方程来描述温度的影响，只是因反应不同，参数有所不同。在ASM2中，受温度影响的动力学参数有14个。分别是：2个水解过程动力学参数：水解速率常数Kh，颗粒性COD的饱和系数K。4个异氧菌过程动力学参数：最大生长速率H，发酵最大速率qfe，反硝化校正因子NO3，溶菌速率常数bH。6个聚磷菌过程动力学参数：PHA的贮存的速率常数qPHA，PP的贮存的速率常数qpp，聚磷菌最大生长速率PAO，聚磷菌的溶菌速率常数bPAO，Xpp的分解速率常数bpp，聚磷菌的分解速率常数bPHA。2个自养菌过程动力学参数：自养菌最大生长速率AUT，自养曲衰减速率bAUT。根据Arreheniu。方程；kTk20·exp（T20)，不同动力学参数的温度修正系数a取值不同。根据IAWQ工作组捡拾的温度修正系数a，不同的动力学参数的温度修正系数见表5。 表5 温度修正系数a动力学参数温度修正系数a水解速率常数Kh0.041颗粒性COD的饱和系