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

活性污泥2号模型的应用与校正国际水质协会（IAWQ）于1995年推出了活性污泥2号模型[1]（ActivatedSludgeModdNo.2，ASM2），ASM2是活性污泥1号模型[2](ActivatedSludgeModelNo.1，ASM1）的扩展，并沿用了该模型的概念，它包含碳有机物氧化、脱氮和生物除磷处理过程，共有19种组分、19种反应、22个化学计量系数及42个动力学参数。ASM2无论是在污水厂的设计、运行管理、改造，还是在废水处理技术的研究和开发方面均有重要的使用价值，在国外已得到了成功广泛的使用[3]。我国在这方面起步很慢，本研究对ASM2在我国城市污水厂的实用性进行了验证。模型所用的污水进出水数据来自常州城北污水厂，污水厂采用A2／O工艺，设计流量为5×104m3／d，生物反应池的水力停留时间为11.5h，污泥负荷为0.15kg［BOD5]（kg［MLSS］.d），污泥的质量浓度为3500mg／L，固体平均停留时间为20d，污泥回流比100％，厌氧池体积为2900m3，缺氧池体积为5800m3，曝气池体积为14500m3。由于该污水厂具有除磷脱氮功能，而ASM1没有考虑除磷过程，所以采用ASM2进行模拟研究。污水处理流程的简化表示如图1所示：1模拟计算结果此稳态计算程序利用Matlab[4-6]软件为平台，并根据ASM2对常州污水厂运行条件进行简化。假定该厂的供气量可以满足好氧池中微生物生化反应所需要的氧气量，即假定曝气池中溶解氧组分的取值恒定在2.0mg／L，饱和溶解氧的质量浓度为8.637mg／L（20℃），系统的操作温度恒定，二沉池不考虑微生物的物质代谢活动，仅起固液分离作用且无活性污泥的积累。在这基础上，将2000年1－12月进水组分输人模型中，可以得到出水的CODcr，NH3－N，TN，TP的稳态结果，模拟结果见图2--图5。从图2一图5中可以看到，CODcr，NH3－N，TN，TP模拟值与污水厂的实测值基本相符。CODcr的最大相对误差36.5%发生在12月份，差值为16.11mg／L，其它11个月份相对误差在30％之内，其中1，2，3，5月的相对误差控制在10％之内，因而较好的模拟了CODcr值。NH3－N，TN和TP指标因为出水中浓度很低，模拟有一定的困难，但模拟结果表明模拟值与实测值在同一个数量级别．对于NH3－N而言，模拟曲线比较平缓，这是因为没有考虑温度对动力学参数的影响，因而模拟值与实测值差别较大。TP的模拟值与实测值的变化趋势相似，曲线很好的拟合了实际出水中TP的变化。偏差的存在有多方面的原因，一是模型假设的部分条件与实际污水厂运行不完全吻合，另外污水厂化验分析不可避免存在一定误差。尽管有一定的偏差，但模拟结果得到了良好的描述，因而利用活性污泥2号模型对城市污水厂进行模拟是可行的。2参数灵敏度和温度分析2．1参数灵敏度分析在输入IAWQ给出的模型动力学和化学计量学参数的典型值[7]后，模拟结果表明出水指标大部分得到了较好的模拟，但部分出水的模拟值与实际值有较大的误差，所以需要对一些参数进行校核，这里以灵敏度分析为基础，将模型中的动力学参数值都增加其初始值的10％，从而计算出水的CODcr，TN，TP，NN3－N的灵敏度，计算结果见表1～表4。由灵敏度计算结果表1——表4分析可知，大部分参数的灵敏度都比较低（小于1），因此采用I-AWQ的推荐参数值进行模拟是合理的。但仍有些参数对出水水质影响较大（灵敏度大于1）。对CODcr影响较大的动力学参数有4个：Kx，bPAO，异养菌KNH3－N，bH。表1水解过程动力学参数灵敏度分析参数CODcrNH3-NTNTP水解速率常数Kh00.0330.0330.067氧的饱和抑制系数Ko20.50.500.5颗粒性COD的饱和系数Kx1111缺氧水解速率降低修正因子η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.1250.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自氧菌和沉淀再溶解过程动力学参数灵敏度分析参数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对NH3－N影响较大的动力学参数有10个：分别是KX，ηfe，异养菌KNH3-N，Kp，μAUT，bAUT，自养菌KH3－N，Kpp，KPHA，异养菌KALK。对TN影响较大的动力学参数有6个：分别是KX，ηfe，bH，bPAO，bAUT，反硝化ηNO3。对TP影响较大的动力学参数有7个：分别是Kx，ηfe，异养菌KNH3-N，Kp，异养菌KALK，Kpp，KPHA。2．2温度影响分析一切生化反应都可以用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。方程；kT＝k20℃·exp（T－20))，不同动力学参数的温度修正系数a取值不同。根据IAWQ工作组捡拾的温度修正系数a，不同的动力学参数的温度修正系数见表5。表5温度修正系数a动力学参数温度修正系数a水解速率常数Kh0.041颗粒性COD的饱和系数Kx-0.11基于基质的最大生长速率μH0.069发酵的最大速率qfe0.069反硝化的速率降低修正固子ηNO3-0.069溶菌速率常数bH0.069PHA贮存的速率常数qPHA0.041PP贮存的速率常数qpp0.041聚磷菌最大生长速率μPAO0.041聚磷菌溶菌速率常数bPAO0.041XPP的分解速率常数bPP0.041XPHA的分解速率常数bPHA0.041自养菌最大生长速率μAUT0.105自养菌衰减速率bAUT0.1053模型的校正以1月份数据为例，考虑平均温度为10℃，根据表5建议的温度系数，修正动力学参数，得到校正模型，从而对1月进水数据进行再模拟，得到的出水模拟值为后模拟值，与未校正模型模拟值（前模拟值）进行比较，结果见图6——图7：从上图中可以看到，经过参数修正后的模拟，比未校正参数的模拟精度显著提高，从而说明有些参数在模型中的重要性，因而模型在实际污水厂应用中，有必要对参数进行修正。4结语本文根据自行开发的ASM2软件，对常州污水厂进行模拟研究，得出以下结论：①实际计算比较证明，利用活性污泥2号模型能充分考虑废水生物处理自身的特点和工艺特性，较为深人地反映活性污泥生物反应过程的主要行为，具有很好的应用性，适合我国污水处理厂的研究管理。②通过参数灵敏度分析，分别对影响CODcr，NH3-N，TN，TP4个指标的参数分别考虑，并考虑温度影响，对不同动力学参数进行温度修正。校正后的ASM2模型的预测精度得到显著提高。