活性污泥2號(hào)模型的應(yīng)用與校正

1、活性污泥2號(hào)模型的應(yīng)用與校正    摘要：采用國(guó)際水質(zhì)協(xié)會(huì)提出的活性污泥2號(hào)模型，以MATLAB為工具編制計(jì)算軟件，采用該軟件對(duì)常州市污水廠進(jìn)行計(jì)算，在輸入典型的動(dòng)力學(xué)和化學(xué)計(jì)量參數(shù)值后，計(jì)算結(jié)果表明出水CODcr指標(biāo)大部分得到較好的模擬，但NH3N、TP和TN出水指標(biāo)模擬值與實(shí)際值有差異。通過對(duì)參數(shù)靈敏度的分析，在考慮溫度影響的條件下，進(jìn)行模型校正，模型校正后，預(yù)測(cè)精度得到顯著提高。 關(guān)鍵詞：城市污水 活性污泥法 數(shù)學(xué)模型  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  國(guó)際水質(zhì)協(xié)會(huì)（IAWQ）于1995年推出了活性污泥 2號(hào)模型1（Activated Sludge Modd No.2，ASM2），ASM2是活性污泥1號(hào)模型

5、2(Activated Sludge Model No.1，ASM1）的擴(kuò)展，并沿用了該模型的概念，它包含碳有機(jī)物氧化、脫氮和生物除磷處理過程，共有19種組分、19種反應(yīng)、22個(gè)化學(xué)計(jì)量系數(shù)及42個(gè)動(dòng)力學(xué)參數(shù)。ASM2無論是在污水廠的設(shè)計(jì)、運(yùn)管理、改造，還是在廢水處理技術(shù)的研究和開發(fā)方面均有重要的使用價(jià)值，在國(guó)外已得到了成功廣泛的使用3。我國(guó)在這方面起步很慢，本研究對(duì)ASM2在我國(guó)城市污水廠的實(shí)用性進(jìn)行了驗(yàn)證。模型所用的污水進(jìn)出水?dāng)?shù)據(jù)來自常州城北污水廠，污水廠采用A2O工藝，設(shè)計(jì)流量為5×104m3d，生物反應(yīng)池的水力停留時(shí)間為11.5h，污泥負(fù)荷為0.15kgBOD5（kgMLSS

6、.d），污泥的質(zhì)量濃度為3500mgL，固體平均停留時(shí)間為20d，污泥回流比100，厭氧池體積為2900m3，缺氧池體積為5800m3，曝氣池體積為14500m3。由于該污水廠具有除磷脫氮功能，而ASM1沒有考慮除磷過程，所以采用ASM2進(jìn)行模擬研究。污水處理流程的簡(jiǎn)化表示如圖1所示： 1 模擬計(jì)算結(jié)果此穩(wěn)態(tài)計(jì)算程序利用Matlab4-6軟件為平臺(tái)，并根據(jù)ASM2對(duì)常州污水廠運(yùn)行條件進(jìn)行簡(jiǎn)化。假定該廠的供氣量可以滿足好氧池中微生物生化反應(yīng)所需要的氧氣量，即假定曝氣池中溶解氧組分的取值恒定在2.0mgL，飽和溶解氧的質(zhì)量濃度為8.637mgL（20），系統(tǒng)的操作溫度恒定，二沉池不考慮微生物的物質(zhì)

7、代謝活動(dòng)，僅起固液分離作用且無活性污泥的積累。在這基礎(chǔ)上，將2000年112月進(jìn)水組分輸人模型中，可以得到出水的CODcr，NH3N，TN，TP的穩(wěn)態(tài)結(jié)果，模擬結(jié)果見圖2-圖5。 從圖2一圖5中可以看到，CODcr，NH3N，TN，TP模擬值與污水廠的實(shí)測(cè)值基本相符。CODcr的最大相對(duì)誤差36.5%發(fā)生在12月份，差值為16.11mgL，其它11個(gè)月份相對(duì)誤差在30之內(nèi)，其中1，2，3，5月的相對(duì)誤差控制在10之內(nèi)，因而較好的模擬了CODcr值。NH3N，TN和TP指標(biāo)因?yàn)槌鏊袧舛群艿?，模擬有一定的困難，但模擬結(jié)果表明模擬值與實(shí)測(cè)值在同一個(gè)數(shù)量級(jí)別對(duì)于NH3N而言，模擬曲線比較平緩，這是因

8、為沒有考慮溫度對(duì)動(dòng)力學(xué)參數(shù)的影響，因而模擬值與實(shí)測(cè)值差別較大。TP的模擬值與實(shí)測(cè)值的變化趨勢(shì)相似，曲線很好的擬合了實(shí)際出水中TP的變化。偏差的存在有多方面的原因，一是模型假設(shè)的部分條件與實(shí)際污水廠運(yùn)行不完全吻合，另外污水廠化驗(yàn)分析不可避免存在一定誤差。盡管有一定的偏差，但模擬結(jié)果得到了良好的描述，因而利用活性污泥2號(hào)模型對(duì)城市污水廠進(jìn)行模擬是可行的。 2 參數(shù)靈敏度和溫度分析 21 參數(shù)靈敏度分析在輸入IAWQ給出的模型動(dòng)力學(xué)和化學(xué)計(jì)量學(xué) 參數(shù)的典型值7后，模擬結(jié)果表明出水指標(biāo)大部分得到了較好的模擬，但部分出水的模擬值與實(shí)際值有 較大的誤差，所以需要對(duì)一些參數(shù)進(jìn)行校核，這里以 靈敏度分析為基礎(chǔ)

9、，將模型中的動(dòng)力學(xué)參數(shù)值都增加 其初始值的10，從而計(jì)算出水的CODcr，TN，TP，NN3N的靈敏度，計(jì)算結(jié)果見表1表4。 由靈敏度計(jì)算結(jié)果表1表4分析可知，大部 分參數(shù)的靈敏度都比較低（小于1），因此采用I- AWQ的推薦參數(shù)值進(jìn)行模擬是合理的。但仍有些 參數(shù)對(duì)出水水質(zhì)影響較大（靈敏度大于1）。對(duì) CODcr影響較大的動(dòng)力學(xué)參數(shù)有4個(gè)：Kx，bPAO， 異養(yǎng)菌KNH3N，bH。表1  水解過程動(dòng)力學(xué)參數(shù)靈敏度分析 參數(shù)CODcrNH3-NTNTP水解速率常數(shù)Kh00.0330.0330.067氧的飽和抑制系數(shù)Ko20.50.500.5顆粒性COD的飽和系數(shù)Kx1111缺氧水解速率

10、降低修正因子NO300.1670.1670.333硝酸鹽的飽和抑制系數(shù)KNO3-0.20.20.40.2厭氧水解速率降低修正因子fe0132表2異氧菌過程動(dòng)力學(xué)參數(shù)靈敏度分析 參數(shù)CODcrNH3-NTNTP基于基質(zhì)的最大生長(zhǎng)速率H0.0830.01700.017基于發(fā)酵基質(zhì)的生長(zhǎng)飽和系數(shù)KF0.150.02500.05氨氮的飽和系數(shù)KNH3-N2202磷的飽和系數(shù)KP010010堿度的飽和系數(shù)KALK0101乙酸的飽和系數(shù)KA0.0250.02500.05發(fā)酵的最大速率qfe0.0330.03300.033溶菌速率常數(shù)bH1.50.252.250.25反硝化的速率降低修正固子NO3-0.12

11、50.53.3750.125表3  聚磷菌過程動(dòng)力學(xué)參數(shù)靈敏度分析 參數(shù)CODcrNH3-NTNTPPHA貯存的速率常數(shù)qPHA00.0330.0330.033聚磷酸鹽的飽和系數(shù)Kpp010010PP貯存的速率常數(shù)qpp0.0670.1330.0670.4Xpp/XPAO的最大比率Kmax00.29400.588PHA的飽和系數(shù)KPHA010020聚磷菌最大生長(zhǎng)速率PAO0000.5聚磷菌溶菌速率常數(shù)bPAO10.510.5XPP的分解速率常數(shù)bPP00.50.50.5XPHA的分解速率常數(shù)bPHA0.50.500.5表4  自氧菌和沉淀再溶解過程動(dòng)力學(xué)參數(shù)靈敏度分析 參數(shù)

12、CODcrNH3-NTNTP自養(yǎng)菌最大生長(zhǎng)速率0.11.30.70.2自養(yǎng)菌衰減速率bAUT0.679.334.670.67自養(yǎng)菌氨氮的飽和系數(shù)KNH3-N010.50.1自養(yǎng)菌堿度的飽和系數(shù)KALK00.60.20.2自養(yǎng)菌氧的飽和系數(shù)KO20.20.80.40.2磷沉淀的速率系數(shù)KPRE00.100.1再溶解速率系數(shù)KRED0.1670.16700.167對(duì)NH3N影響較大的動(dòng)力學(xué)參數(shù)有 10個(gè)：分別是KX，fe，異養(yǎng)菌KNH3-N，Kp，AUT，bAUT，自養(yǎng)菌KH3N，Kpp，KPHA，異養(yǎng)菌KALK。對(duì)TN影響較大的動(dòng)力學(xué)參數(shù)有6個(gè)：分別是KX，fe，bH，bPAO，bAUT，反硝化

13、NO3。 對(duì)TP影響較大的動(dòng)力學(xué)參數(shù)有7個(gè)：分別是Kx，fe，異養(yǎng)菌KNH3-N，Kp，異養(yǎng)菌KALK，Kpp，KPHA。22 溫度影響分析一切生化反應(yīng)都可以用Arrehenius方程來描述溫度的影響，只是因反應(yīng)不同，參數(shù)有所不同。在ASM2中，受溫度影響的動(dòng)力學(xué)參數(shù)有14個(gè)。分別是：2個(gè)水解過程動(dòng)力學(xué)參數(shù)：水解速率常數(shù)Kh，顆粒性COD的飽和系數(shù)K。4個(gè)異氧菌過程動(dòng)力學(xué)參數(shù)：最大生長(zhǎng)速率H，發(fā)酵最大速率qfe，反硝化校正因子NO3，溶菌速率常數(shù)bH。6個(gè)聚磷菌過程動(dòng)力學(xué)參數(shù)：PHA的貯存的速率常數(shù)qPHA，PP的貯存的速率常數(shù)qpp，聚磷菌最大生長(zhǎng)速率PAO，聚磷菌的溶菌速率常數(shù)bPAO，Xpp的分解速率常數(shù)bpp，聚磷菌的分解速率常數(shù)bPHA。2個(gè)自養(yǎng)菌過程動(dòng)力學(xué)參數(shù)：自養(yǎng)菌最大生長(zhǎng)速率AUT，自養(yǎng)曲衰減速率bAUT。根據(jù)Arreheniu。方程；kTk20·exp（T20)，不同動(dòng)力學(xué)參數(shù)的溫度修正系數(shù)a取值不同。根據(jù)IAWQ工作組撿拾的溫度修正系數(shù)a，不同的動(dòng)力學(xué)參數(shù)的溫度修正系數(shù)見表5。 表5 溫度修正系數(shù)a動(dòng)力學(xué)參數(shù)溫度修正系數(shù)a水解速率常數(shù)Kh0.041顆粒性COD的飽和系