陳曉龍，楊海真，顧國維
（同濟大學污染控制與資源化研究國家重點實驗室，上海 200092）

　　摘　要：采用國際水質協會提出的活性污泥2號模型，以MATLAB為工具編制計算軟件，采用該軟件對常州市污水廠進行計算，在輸入典型的動力學和化學計量參數值后，計算結果表明出水CODcr指標大部分得到較好的模擬，但NH₃－N、TP和TN出水指標模擬值與實際值有差異。通過對參數靈敏度的分析，在考慮溫度影響的條件下，進行模型校正，模型校正后，預測精度得到顯著提高。
　　關鍵詞：城市污水；污水處理；活性污泥法；數學模型
　　中圖分類號：X703
　　文獻標識碼：A
　　文章編號：1009－2455（2003）01－0001－04

Application and Calibration of Activated Sludge Model No. 2
CHEN Xiao-long, YANG Hai-zhen, GU Guo-wei
(State Key fuboratory Of Pollution Control and Resource Reuse, Tongji University, Shanghai 200092, China)

　　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 the 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 temperature, 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號模型^[2](Activated Sludge Model No.1，ASM1）的擴展，并沿用了該模型的概念，它包含碳有機物氧化、脫氮和生物除磷
處理過程，共有19種組分、19種反應、22個化學計量系數及42個動力學參數。ASM2無論是在污水廠的設計、運行管理、改造，還是在廢水處理技術的研究和開發方面均有重要的使用價值，在國外已得到了成功廣泛的使用[3]。我國在這方面起步很慢，本研究對ASM2在我國城市污水廠的實用性進行了驗證。
　　模型所用的污水進出水數據來自常州城北污水廠，污水廠采用A²／O工藝，設計流量為5×10⁴m³／d，生物反應池的水力停留時間為11.5h，污泥負荷為0.15kg［BOD₅]（kg［MLSS］.d），污泥的質量濃度為3500mg／L，固體平均停留時間為20d，污泥回流比100％，厭氧池體積為2900m³，缺氧池體積為5800m³，曝氣池體積為14500m³。由于該污水廠具有除磷脫氮功能，而ASM1沒有考慮除磷過程，所以采用ASM2進行模擬研究。污水處理流程的簡化表示如圖1所示：

1 模擬計算結果

　　此穩態計算程序利用Matlab^[4-6]軟件為平臺，并根據ASM2對常州污水廠運行條件進行簡化。假定該廠的供氣量可以滿足好氧池中微生物生化反應所需要的氧氣量，即假定曝氣池中溶解氧組分的取值恒定在2.0mg／L，飽和溶解氧的質量濃度為8.637mg／L（20℃），系統的操作溫度恒定，二沉池不考慮微生物的物質代謝活動，僅起固液分離作用且無活性污泥的積累。在這基礎上，將2000年1－12月進水組分輸人模型中，可以得到出水的CODcr，NH₃－N，TN，TP的穩態結果，模擬結果見圖2--圖5。

　　

　　從圖2一圖5中可以看到，CODcr，NH₃－N，TN，TP模擬值與污水廠的實測值基本相符。CODcr的最大相對誤差36.5%發生在12月份，差值為16.11mg／L，其它11個月份相對誤差在30％之內，其中1，2，3，5月的相對誤差控制在10％之內，因而較好的模擬了CODcr值。
　　NH₃－N，TN和TP指標因為出水中濃度很低，模擬有一定的困難，但模擬結果表明模擬值與實測值在同一個數量級別．對于NH3－N而言，模擬曲線比較平緩，這是因為沒有考慮溫度對動力學參數的影響，因而模擬值與實測值差別較大。TP的模擬值與實測值的變化趨勢相似，曲線很好的擬合了實際出水中TP的變化。
　　偏差的存在有多方面的原因，一是模型假設的部分條件與實際污水廠運行不完全吻合，另外污水廠化驗分析不可避免存在一定誤差。盡管有一定的偏差，但模擬結果得到了良好的描述，因而利用活性污泥2號模型對城市污水廠進行模擬是可行的。

2 參數靈敏度和溫度分析

2．1 參數靈敏度分析
　　在輸入IAWQ給出的模型動力學和化學計量學 參數的典型值^[7]后，模擬結果表明出水指標大部分得到了較好的模擬，但部分出水的模擬值與實際值有 較大的誤差，所以需要對一些參數進行校核，這里以 靈敏度分析為基礎，將模型中的動力學參數值都增加 其初始值的10％，從而計算出水的CODcr，TN，TP，NN₃－N的靈敏度，計算結果見表1～表4。
　　 由靈敏度計算結果表1——表4分析可知，大部 分參數的靈敏度都比較低（小于1），因此采用I- AWQ的推薦參數值進行模擬是合理的。但仍有些 參數對出水水質影響較大（靈敏度大于1）。對 CODcr影響較大的動力學參數有4個：Kx，bPAO， 異養菌KNH₃－N，b_H。

表1 水解過程動力學參數靈敏度分析 參數 CODcr NH₃-N TN TP 水解速率常數Kh 0 0.033 0.033 0.067 氧的飽和抑制系數Ko₂ 0.5 0.5 0 0.5 顆粒性COD的飽和系數Kx 1 1 1 1 缺氧水解速率降低修正因子ηNO₃ 0 0.167 0.167 0.333 硝酸鹽的飽和抑制系數KNO₃^- 0.2 0.2 0.4 0.2 厭氧水解速率降低修正因子η_fe 0 1 3 2

表2　異氧菌過程動力學參數靈敏度分析 參數 CODcr NH₃-N TN TP 基于基質的最大生長速率μ_H 0.083 0.017 0 0.017 基于發酵基質的生長飽和系數K_F 0.15 0.025 0 0.05 氨氮的飽和系數KNH₃-N 2 2 0 2 磷的飽和系數K_P 0 10 0 10 堿度的飽和系數K_ALK 0 1 0 1 乙酸的飽和系數K_A 0.025 0.025 0 0.05 發酵的最大速率q_fe 0.033 0.033 0 0.033 溶菌速率常數bH 1.5 0.25 2.25 0.25 反硝化的速率降低修正固子ηNO₃^- 0.125 0.5 3.375 0.125

表3 聚磷菌過程動力學參數靈敏度分析 參數 CODcr NH₃-N TN TP PHA貯存的速率常數q_PHA 0 0.033 0.033 0.033 聚磷酸鹽的飽和系數K_pp 0 10 0 10 PP貯存的速率常數q_pp 0.067 0.133 0.067 0.4 X_pp/X_PAO的最大比率Kmax 0 0.294 0 0.588 PHA的飽和系數K_PHA 0 10 0 20 聚磷菌最大生長速率μ_PAO 0 0 0 0.5 聚磷菌溶菌速率常數b_PAO 1 0.5 1 0.5 X_PP的分解速率常數b_PP 0 0.5 0.5 0.5 X_PHA的分解速率常數b_PHA 0.5 0.5 0 0.5

表4 自氧菌和沉淀再溶解過程動力學參數靈敏度分析 參數 CODcr NH₃-N TN TP 自養菌最大生長速率 0.1 1.3 0.7 0.2 自養菌衰減速率b_AUT 0.67 9.33 4.67 0.67 自養菌氨氮的飽和系數K_NH3-N 0 1 0.5 0.1 自養菌堿度的飽和系數KALK 0 0.6 0.2 0.2 自養菌氧的飽和系數KO₂ 0.2 0.8 0.4 0.2 磷沉淀的速率系數K_PRE 0 0.1 0 0.1 再溶解速率系數K_RED 0.167 0.167 0 0.167

　　對NH₃－N影響較大的動力學參數有 10個：分別是KX，ηfe，異養菌KNH₃-N，Kp，μAUT，bAUT，自養菌KH₃－N，Kpp，KPHA，異養菌KALK。
　　對TN影響較大的動力學參數有6個：分別是KX，ηfe，bH，bPAO，bAUT，反硝化ηNO₃。 對TP影響較大的動力學參數有7個：分別是Kx，ηfe，異養菌KNH3-N，Kp，異養菌KALK，Kpp，KPHA。
2．2 溫度影響分析
　　一切生化反應都可以用Arrehenius方程來描述溫度的影響，只是因反應不同，參數有所不同。在ASM2中，受溫度影響的動力學參數有14個。分別是：
　?、?個水解過程動力學參數：水解速率常數Kh，顆粒性COD的飽和系數K。
　?、?個異氧菌過程動力學參數：最大生長速率μH，發酵最大速率qfe，反硝化校正因子ηNO₃，溶菌速率常數bH。
　?、?個聚磷菌過程動力學參數：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 水解速率常數Kh 0.041 顆粒性COD的飽和系數Kx -0.11 基于基質的最大生長速率μ_H 0.069 發酵的最大速率q_fe 0.069 反硝化的速率降低修正固子ηNO₃^- 0.069 溶菌速率常數bH 0.069 PHA貯存的速率常數q_PHA 0.041 PP貯存的速率常數q_pp 0.041 聚磷菌最大生長速率μ_PAO 0.041 聚磷菌溶菌速率常數b_PAO 0.041 XPP的分解速率常數b_PP 0.041 XPHA的分解速率常數b_PHA 0.041 自養菌最大生長速率μ_AUT 0.105 自養菌衰減速率b_AUT 0.105

3 模型的校正

　　 以 1月份數據為例，考慮平均溫度為10℃，根據表5建議的溫度系數，修正動力學參數，得到校正模型，從而對1月進水數據進行再模擬，得到的出水模擬值為后模擬值，與未校正模型模擬值（前模擬值）進行比較，結果見圖6——圖7：

　　 從上圖中可以看到，經過參數修正后的模擬，比未校正參數的模擬精度顯著提高，從而說明有些參數在模型中的重要性，因而模型在實際污水廠應用中，有必要對參數進行修正。

4 結語

　　 本文根據自行開發的ASM2軟件，對常州污水廠進行模擬研究，得出以下結論：
　　 ①實際計算比較證明，利用活性污泥2號模型能充分考慮廢水生物處理自身的特點和工藝特性，較為深人地反映活性污泥生物反應過程的主要行為，具有很好的應用性，適合我國污水處理廠的研究管理。
　　 ②通過參數靈敏度分析，分別對影響CODcr，NH3-N，TN，TP4個指標的參數分別考慮，并考慮溫度影響，對不同動力學參數進行溫度修正。校正后的ASM2模型的預測精度得到顯著提高。

參考文獻：

　　[1]Henzn M, Gujer W, Mino T, et al. Activated Sludge Model(2):I-AWC Scientific and Technical Report No.3[R]. London: IAWQ,1995.
　　[2] Henze M, Grady C P L Jr, Gujer W, el al. Activated Sludge Model No.1, IAWC Scientific and Technical Report No.1[R].London: I-AWQ,1987.
　　[3] Ozer Cinar, Glen T Daigger, Stephen P Graef, et al . Evaluation of the IAWQ Activated Sludge Model No.2 Using Steady-State Data from Four Full-Scale Wastewater Treatment Plants[J]. Water En-vironment Research, 1978,70(6):1216-1223.
　　[4] 王沫然. MATLABSX與科學計算[M］．北京：清華大學出版 社，2000．
　　[5] 許波，劉征,MATLAB程數學應用［M］．北京：清華大學出版 社，2000.
　　[6] 潘曉輝MATLABS全攻略寶典[M].北京：中國水利水電出 版社，1999
　　[7] 張亞雷，李詠梅 活性污泥數學模型[M]．上海：同濟大學出 版社，2002

---

　　作者簡介：陳曉龍（1977－），男，浙江紹興人，現攻讀同濟大學環境工程碩士學位，主要從事污水廠自動化控制研究，電話（021）65984254。