周 振,胡大龍,吳志超,王志偉,杜興治,蔣玲燕,喬衛(wèi)敏 (.上海電力學院環(huán)境與化學工程學院,上海 00090；.同濟大學環(huán)境科學與工程學院,污染控制與資源化研究國家重點實驗室,上海 0009；.上海城投污水處理有限公司,上海 00)

基于數(shù)學模型的多模式AAO系統(tǒng)運行優(yōu)化研究

周 振1*,胡大龍1,吳志超2,王志偉2,杜興治3,蔣玲燕3,喬衛(wèi)敏1(1.上海電力學院環(huán)境與化學工程學院,上海 200090；2.同濟大學環(huán)境科學與工程學院,污染控制與資源化研究國家重點實驗室,上海 200092；3.上海城投污水處理有限公司,上海 201203)

構建了綜合考慮出水水質、污泥產量和系統(tǒng)能耗的運行成本指數(shù)CPI,并利用數(shù)學模型對多模式厭氧/缺氧/好氧(AAO)系統(tǒng)的AAO、倒置AAO(RAAO)和缺氧/好氧(AO)3種模式進行對比優(yōu)化研究.在處理成本相近的前提下,在污泥齡5～25d范圍內AAO模式的污染物去除效率和聚磷菌濃度明顯高于RAAO和AO模式.在運行模式篩選的基礎上,通過回流比優(yōu)化確定了排放標準約束下AAO工藝的達標運行區(qū)域和最佳運行工況.動態(tài)模擬結果表明,優(yōu)化工況能夠顯著改善出水水質,出水高于一級A的時間由78.4%下降至37.7%,CPI降低3.9%.

污水處理；優(yōu)化；厭氧/缺氧/好氧工藝；活性污泥；數(shù)學模型

近年來,水資源的短缺和水污染的加劇使得脫氮除磷成為污水處理廠的首要任務.作為最為廣泛應用的污水脫氮除磷技術,厭氧/缺氧/好氧(AAO)工藝通過不同的反應條件營造和復雜的操作變量控制實現(xiàn)了污水中COD和營養(yǎng)元素的同步去除[1-4].張波[5]通過 AAO 工藝厭氧和缺氧段的倒置提出了倒置 AAO(RAAO)工藝,該工藝支持者認為其具有前置反硝化強化脫氮、“群體效應”強化除磷、避免電子受體干擾厭氧釋磷、厭氧合成的聚 β羥基烷酸(XPHA)高效利用等優(yōu)點[6-8].然而,上述論述更多是基于推斷,AAO 和RAAO的運行特性尚缺乏系統(tǒng)的對比分析.

多模式AAO系統(tǒng)可根據(jù)季節(jié)變化、排放要求和進水水質在AAO、RAAO和缺氧/好氧(AO)3種模式間靈活切換[9](圖1).由于AAO工藝受污泥齡(SRT)、污泥回流(RAS)、混合液回流(MLR)和曝氣強度等條件影響很大[10-13],多模式 AAO系統(tǒng)在增加系統(tǒng)靈活性的同時也帶來了運行管理的挑戰(zhàn),通過實驗優(yōu)化存在很大難度.數(shù)學模擬則能很方便地優(yōu)化污水處理工藝,確定最佳運行模式和工況條件[10,14-16].本研究擬在綜合分析出水水質、污泥產量和系統(tǒng)能耗的基礎上構建運行成本指數(shù),并基于此對多模式AAO系統(tǒng)的3種運行模式對比分析與優(yōu)化以確定最佳運行條件.

1 材料與方法

1.1 多模式AAO系統(tǒng)

以白龍港污水處理廠多模式AAO系統(tǒng)為優(yōu)化對象,其工藝流程如圖1所示.該系統(tǒng)可通過閥門V1和V2的啟閉在3種模式間切換: V1關閉V2開啟,為AAO模式;V1開啟V2關閉,為RAAO模式; V1和V2同時關閉,為AO模式.

圖1 多模式AAO系統(tǒng)工藝流程Fig.1 Flow chart of the multi-mode anaerobic/anoxic/oxic (MAAO)system

多模式 AAO系統(tǒng)處理水量為 200000m3/d,選擇、厭氧、缺氧和好氧水力停留時間分別為0.3、1.3、2.3和 8.1h.系統(tǒng)運行工況為進水分配1:1,MLR和RAS比均為50%,SRT為15d,好氧區(qū)溶解氧(SO)為2.0mg/L,水溫20.℃進水水質如表1所示, COD組分參照Zhou等[17].進水SO、氮氣(SN2)、硝化菌(XAUT)、聚磷酸鹽(XPP)和 XPHA輸入值均為活性污泥模型[18]推薦值 0.001mg/L,金屬氫氧化物(XMeOH)和磷酸鹽(XMeP)取0.01mg/L.

表1 白龍港污水處理廠進水水質特性Table 1 Influent characteristics of the Bailonggang wastewater treatment plant

1.2 模型構建與參數(shù)校正

模型構建采用WEST 3.7.5軟件完成.多模式AAO系統(tǒng)水力模型采用串聯(lián)反應器模式.生化反應模型選擇 ASM2dTemp模型,該模型可模擬不同溫度下的碳氧化、硝化-反硝化以及包括反硝化除磷在內的生物除磷過程等.二沉池模型選用一維分層模型–雙指數(shù)模型[19].本研究所用模型及其參數(shù)校正與驗證參見文獻[9,20].

1.3 模擬結果的評價指標

本研究在COST 624/682基準[21]及Gernaey等[22]模擬經驗的基礎上建立了針對水質、能耗、污泥產量等方面的評價指標體系,該體系主要用于初始時刻 t0到終止時刻 tf的動態(tài)模擬,其中部分評價指標亦可簡化后用于穩(wěn)態(tài)模擬評估.

1.3.1 出水水質指數(shù)(EQ)出水水質指數(shù) EQ(kg污染物/d)計算如下:

式中:Qe為出水流量,m3/d;污染物負荷 PU(t)為各組分污染物之和.

各組分污染物負荷PUk(mg/L)可計算如下

式中:βk和 Ck分別為出水組分權重因子和濃度.西歐[21-22]通 常 取 βTSS=2;βCOD=1;βBOD=2;βTN=20;βTP=20.而根據(jù)我國《排污費征收標準管理辦法》[23]核算,βTSS=0.25;βCOD=1;βBOD=2;βTN=0;βTP=4,該辦法更側重于有機物去除.本文綜合考慮“十二五”氨氮減排要求和脫氮除磷,擬取值如下:βTSS=0.25;βCOD=1;βBOD=2;βNH=30;βTN=20;βTP=20.

1.3.2 出水超標(EV)出水超標評估是污水處理廠優(yōu)化的重要內容[21].本文擬按照《城鎮(zhèn)污水處理廠污染物排放標準》(GB 18918-2002)[24]一級A標準優(yōu)化評估.出水超標以超標次數(shù)(NV)和超標時間百分比(PV)定量衡量.NV指超標時間段的個數(shù), PV則是一段時間內超標時間所占百分比.比較而言,PV更能準確反映出水超標的時間長度,兩者取間隔15min的輸出數(shù)據(jù)計算.

1.3.3 運行成本 運行成本需考慮污泥產量、水泵能耗、曝氣能耗和混合能耗 4方面.污泥產量PAS(kg/d)為剩余污泥排放總量,可計算為:

式中:Qw為剩余污泥排放流量,m3/d;Xw為剩余污泥濃度,mg/L.水泵能耗PE(kW·h/d)計算為:

式中:Qa和Qr分別為MLR和RAS流量,m3/d.曝氣能耗 AE(kW·h/d)計算應基于COST基準并考慮反應器體積[25]:

式中:KLa,i和 Vi為反應器 i的傳質系數(shù)和體積,Vref為 1333m3.厭氧段和缺氧段需要攪拌混合以防污泥沉降.混合能耗 ME(kW·h/d)是反應區(qū)體積的函數(shù):

其中,f(Vi)定義如下

1.3.4 運行成本性能參數(shù) 運行成本性能參數(shù)(CPI)可作為污水處理廠運行情況的總體衡量指標,計算如下:

式中:α 值采用如下推薦值[22]:αEQ=50;αAE=αPE=αME=25;αAS=75.出水污染物超標罰款(EF)部分的α值為 αEQ的 2倍.基于這一定義,出水中每一單位污染物排放均需繳納一定數(shù)量的罰金.

2 結果與討論

2.1 不同SRT下3種模式的運行特性對比

圖2為不同SRT下多模式AAO系統(tǒng)3種運行模式出水水質模擬值的對比.當 SRT范圍為5～25d時,20℃下3種模式的出水COD模擬值均在 32～39mg/L間波動,低于一級 A標準限定值.其中,AAO和RAAO的模擬曲線幾乎重合,而AO工藝的模擬值則略高.

圖2 多模式AAO系統(tǒng)出水水質隨SRT的變化Fig.2 Variations of effluent pollutants of the MAAO system with SRT

由圖 2(b)和 2(c)知,20℃時 3種模式的出水SNH、TN與SRT曲線呈類似趨勢,這是因為硝化是脫氮的前提[1,26].兩條曲線拐點均在7.5d左右:當SRT由7.5d降低到5.0d時,出水SNH和TN會迅速上升;當SRT＞7.5d時,出水SNH和TN分別低于2.0和14.0mg/L且緩慢降低.這也說明活性污泥系統(tǒng)存在硝化所需的最小SRT[27],且SRT＞7.5d時3種模式的SNH和TN均能穩(wěn)定達到一級A的要求.值得注意的是,低泥齡時 RAAO 工藝更容易出現(xiàn)硝化和脫氮效果惡化的現(xiàn)象.

由圖2(d)知,3種模式下出水TP模擬值均隨SRT的升高而升高,這說明低SRT有利于生物除磷.AAO出水TP要明顯低于RAAO和AO.如果要達到一級A對出水氮和磷的要求,AAO和AO模式的SRT范圍分別為6～15和6～10d,而RAAO在5～25d的范圍內無法達到一級A的要求.

圖3為20℃時3種模式下好氧段混合液懸浮固體(MLSS)及微生物濃度隨SRT的變化.在3種運行模式下,MLSS均隨SRT增長近似呈線性增長,且 3者相差不大.AO 模式下異養(yǎng)菌(XH)濃度要高于RAAO和AAO,在SRT升高時該趨勢尤為明顯.3種模式下XAUT濃度均隨SRT增加而增加,其中RAAO的XAUT略低于其他兩種模式.在3種運行模式中,AAO的XPAO濃度明顯高于其它兩種模式,而AO又高于RAAO,這一模擬結果也可解釋圖2(d)中3種模式出水TP的差異.

圖4為3種模式綜合評價指標CPI隨SRT的變化趨勢對比.由圖4知,CPI-SRT曲線在SRT為 7.5d時存在極小值,在 SRT由 7.5d增加到25.0d時,CPI基本不變.CPI-SRT曲線與出水SNH和TN曲線類似,這說明3種模式下運行成本主要受出水水質影響.比較而言,RAAO模式的CPI最高,而AAO模式略高于AO模式.

圖4 多模式AAO系統(tǒng)CPI隨SRT的變化Fig.4 Variations of CPI in the MAAO system with SRT

2.2 多模式AAO系統(tǒng)的長期運行優(yōu)化

通過 3種模式對比分析知,在處理成本相近時,AAO的污染物去除效率明顯高于 RAAO和AO,特別是對生物除磷而言.因此,本部分運行優(yōu)化將針對AAO模式.由圖2和4知,為穩(wěn)定達到一級A,SRT可設定為12d.除SRT外,RAS和MLR也會顯著影響 AAO的效果[1,10,12-13].本部分將模擬分析不同回流比對系統(tǒng)性能的影響,并通過構建系統(tǒng)操作圖確定最優(yōu)工況.RAS和MLR變化范圍分別為50%～125%和0～200%,步長取25%.

圖5為20℃時MLR和RAS對AAO系統(tǒng)出水TN、TP和CPI的影響.由于TN和TP在脫氮除磷系統(tǒng)中的重要性,圖5(a)可作為AAO的達標排放操作圖.曲線TP=0.5mg/L、TN=15mg/L和操作圖邊界間的區(qū)域可作為一級A達標運行區(qū)域;而曲線TP=1.0mg/L和邊界圍成的區(qū)域則是一級B達標運行區(qū).圖5(a)中所有MLR和RAS值均可滿足二級標準的要求.在圖 5(b)的等高線圖中存在CPI為16.27k€/d的最小值點(125,75),該點在一級B達標運行區(qū)內;如果考慮一級A的達標要求, CPI在(150,50)處取最小值16.57k€/d.

圖5 AAO系統(tǒng)出水TN、TP和CPI隨MLR和RAS的變化(SRT=12d,20)℃Fig.5 Variations of effluent TN, TP and CPI with MLR and RAS in the AAO system (SRT=12d, 20)℃

圖6為AAO系統(tǒng)按照基本工況與一級A約束優(yōu)化工況運行時的出水水質對比.由圖 6的動態(tài)模擬數(shù)據(jù)計算得到的出水水質、污泥產量和系統(tǒng)能耗指標如表2所示.由圖6可知,優(yōu)化工況的出水 COD略低于基本工況.由于優(yōu)化工況的SRT(12d)低于基本工況(15d),所以優(yōu)化工況下出水SNH略高.優(yōu)化工況降低SRT實際上是在排放標準約束下脫氮和除磷效率的有效平衡.圖 6中優(yōu)化工況的出水TN和TP均低于基本工況.由表2可知,工況優(yōu)化后EQ由62.70t/d降低至54.18t/d,出水水質明顯改善,出水超標時間由 78.4%下降至37.7%,CPI指數(shù)降低了3.9%.

圖6 20℃時AAO系統(tǒng)出水水質的優(yōu)化對比Fig.6 Effluent quality of the AAO system under basic and optimal operational mode

表2 不同工況條件下AAO系統(tǒng)性能參數(shù)匯總Table 2 Performance indices of the AAO system under different operational conditions

3 結論

3.1 在處理成本相近的前提下, SRT 5～25d范圍內AAO模式的污染物去除效率和聚磷菌濃度明顯高于RAAO和AO模式.

3.2 通過內外回流比優(yōu)化確定了一級A約束條件下AAO工藝的達標運行區(qū)域,并以CPI指數(shù)為目標函數(shù)確定最佳運行工況為 SRT=12d,MLR=150%, RAS=50%.

3.3 動態(tài)模擬結果表明,優(yōu)化工況能夠顯著改善出水水質,出水超過一級A的時間由78.4%下降至37.7%,CPI指數(shù)降低了3.9%.

[1]Tchobanoglous G, Burton F L, Stensel H D. Wastewater Engineering: Treatment and Reuse (4th)[M]. New York: Metcalf and Eddy Incorporation, 2003.

[2]Pai T Y. Modeling nitrite and nitrate variations in A2O process under different return oxic mixed liquid using an extended model[J]. Process Biochem., 2007,42(6):978-987.

[3]Bekir Ersu C, Ong S K, Arslankaya E, et al. Comparison of recirculation configurations for biological nutrient removal in a membrane bioreactor [J]. Water Res., 2008,42(6):1651-1663.

[4]曾 薇,李 磊,楊瑩瑩,等.A2O 工藝處理生活污水短程硝化反硝化的研究 [J]. 中國環(huán)境科學, 2010,30(5):625-632.

[5]張 波.城市污水生物脫氮除磷技術工藝與機理研究 [D]. 上海:同濟大學, 1996.

[6]陳洪斌,唐賢春,何群彪,等.倒置AAO工藝聚磷微生物的吸磷行為 [J]. 中國環(huán)境科學, 2007,27(1):49-53.

[7]傅 鋼,董 濱,周增炎,等.倒置AAO工藝的設計特點與運行參數(shù) [J]. 中國給水排水, 2004,20(9):53-55.

[8]張 波,高廷耀.倒置A2/O工藝的原理與特點研究 [J]. 中國給水排水, 2000,16(7):11-15.

[9]周 振,吳志超,王志偉,等. 多模式厭氧/缺氧/好氧污水處理工藝的穩(wěn)態(tài)與動態(tài)模擬 [J]. 環(huán)境科學, 2013,34(4): 1442-1447.

[10]Fang F, Ni B, Li W, et al. A simulation-based integrated approach to optimize the biological nutrient removal process in a full-scale wastewater treatment plant [J]. Chem. Eng. J., 2011,174(2):635-643.

[11]Kennedy K J, McHarg A M. Optimization of municipal wastewater biological nutrient removal using ASM2d [J]. J.Environ. Eng. Sci., 2007,6(1):31-43.

[12]Makinia J, Rosenwinkel K H, Spering V. Long-term simulation of the activated sludge process at the Hanover-Gümmerwald pilot WWTP [J]. Water Res., 2005,39(8):1489-1502.

[13]Ni B J, Yu H Q. Model-based analysis on growth of activated sludge in a sequencing batch reactor [J]. Appl. Microbiol.Biotechnol., 2007,77(3):723-731.

[14]Xie W M, Zhang R, Li W W, et al. Simulation and optimization of a full-scale Carrousel oxidation ditch plant for municipal wastewater treatment [J]. Biochem. Eng. J., 2011,56(1/2):9-16.

[15]Souza S, Araújo O, Coelho M. Model-based optimization of a sequencing batch reactor for biological nitrogen removal [J].Bioresour. Technol., 2008,99(8):3213-3223.

[16]Ostace G S, Baeza J A, Guerrero J, et al. Development and economic assessment of different WWTP control strategies for optimal simultaneous removal of carbon, nitrogen and phosphorus [J]. Comput. Chem. Eng., 2013,53:164-177.

[17]Zhou Z, Wu Z, Wang Z, et al. COD fractionation and parameter estimation for combined sewers by respirometric tests [J]. J.Chem. Technol. Biotechnol., 2008,83(12):1596-1601.

[18]Henze M, Gujer W, Mino M, et al. Activated sludge models ASM1, ASM2, ASM2d and ASM3 [M]. London: IWA Publishing,2000.

[19]周 振,吳志超,顧國維,等.沉淀池污泥層高度的模擬與控制[J]. 中國環(huán)境科學, 2008,28(3):274-278.

[20]Zhou Z, Wu Z, Wang Z, et al. Simulation and performance evaluation of the anoxic/anaerobic/aerobic process for biological nutrient removal [J]. Korean J. Chem. Eng., 2011,28(5):1233-1240.

[21]Copp J B. The COST simulation benchmark: description and simulator manual [M]. Luxembourg: Office for Official Publication of the European Community, 2002.

[22]Gernaey K V, van Loosdrecht M C M, Henze M, et al. Activated sludge wastewater treatment plant modelling and simulation:State of the art [J]. Environ. Modell. Softw., 2004,19(9):763-783.

[23]中華人民共和國國務院.排污費征收標準管理辦法(國務院令字第369號)[Z]. 北京:中華人民共和國國務院, 2003.

[24]GB 18918-2002 城鎮(zhèn)污水處理廠污染物排放標準 [S].

[25]Stare A, Vre?ko D, Hvala N, et al. Comparison of control strategies for nitrogen removal in an activated sludge process in terms of operating costs: a simulation study [J]. Water Res.,2007,41(9):2004-2014.

[26]Head M A, Oleszkiewicz J A. Bioaugmentation for nitrification at cold temperatures [J]. Water Res., 2004,38(3):523-530.

[27]Hao X, Van Loosdrecht M C M, Meijer S C, et al. Model-based evaluation of two BNR processes - UCT and A2N [J]. Water Res.,2001,35(12):2851-2860.

Model-based optimization of the multimode anaerobic/anoxic/aerobic system.

ZHOU Zhen1*, HU Da-long1, WU Zhichao2, WANG Zhi-wei2, DU Xing-zhi3, JIANG Ling-yan3, QIAO Wei-min1(1.College of Environmental and Chemical Engineering, Shanghai University of Electric Power, Shanghai 200090, China；2.College of Environmental Science and Engineering, State Key Laboratory of Pollution Control and Resource Reuse, Tongji University, Shanghai 200092, China；3.Shanghai Chengtou Wastewater Treatment Limited Company, Shanghai 201203, China). China Environmental Science,2014,34(7)：1734～1739

The cost performance index (CPI)of a wastewater treatment plant operating with multimode anaerobic/anoxic/oxic (AAO)system was calculated based on effluent quality, sludge production and energy consumption using a mathematical model. The multimode AAO system could be operated under three modes: AAO, reversed AAO(RAAO)and anoxic/oxic (AO), and CPI for each mode was studied and compared. On the promise of similar operation costs, the AAO mode showed higher pollutants removal efficiency and concentration of phosphorus-accumulating organisms than RAAO and AO modes at sludge retention time of 5～25d. The optimal operation parameters and qualified operational region under discharge standard restraint were established for the AAO process by regulating recirculation ratios. Dynamic simulation results showed that under optimized operating condition, effluent quality was significantly improved with percentage of violation decreased from 78.4% to 37.7%, and the CPI index was reduced by 3.9%.

wastewater treatment；optimization；anaerobic/anoxic/aerobic process；activated sludge；mathematical model

X703

A

1000-6923(2014)07-1734-06

2013-10-25

國家高技術研究發(fā)展計劃(863)項目(2012AA063403);上海市科委能力建設計劃項目(12250500900)

* 責任作者, 副教授, zhouzhen@shiep.edu.cn

周 振(1981-),男,山東臨沂人,副教授,博士,主要從事水污染控制方面的研究.發(fā)表論文60余篇.