ASBR處理高濃度有機(jī)廢水?dāng)?shù)學(xué)模型

靳文東 岳秀萍

1 概述

近十幾年來(lái)國(guó)內(nèi)外學(xué)者對(duì)ASBR反應(yīng)器的啟動(dòng)[1-7]、工藝和操作特性[8-12]、廢水處理應(yīng)用研究[13-15]等方面進(jìn)行了研究，在厭氧條件下，糖類(lèi)和簡(jiǎn)單有機(jī)物等易降解有機(jī)物能很快轉(zhuǎn)化為揮發(fā)酸，然而只有少數(shù)幾種厭氧微生物能夠?qū)⒁宜猁}代謝轉(zhuǎn)化成甲烷，容易造成VFA積累，抑制甲烷菌的活性，厭氧反應(yīng)器的“酸化”會(huì)使處理過(guò)程徹底失敗。

本文根據(jù)化學(xué)反應(yīng)動(dòng)力學(xué)基本原理，借鑒MATLAB數(shù)學(xué)工具對(duì)ASBR模擬試驗(yàn)數(shù)據(jù)進(jìn)行處理，建立ASBR處理高濃度有機(jī)廢水的反應(yīng)速率數(shù)學(xué)模型。旨在及時(shí)掌握和有效控制反應(yīng)器中底物(有機(jī)物)降解速率、產(chǎn)物(甲烷)生成速率和中間產(chǎn)物(揮發(fā)性脂肪酸)的積累程度，并了解進(jìn)水中不同的堿度對(duì)這些動(dòng)態(tài)因素的影響，以保證ASBR工藝長(zhǎng)期穩(wěn)定運(yùn)行。

2 試驗(yàn)裝置與材料

2.1 試驗(yàn)裝置

ASBR反應(yīng)器用有機(jī)玻璃柱加工而成，有效容積6 L，上部留有2 L的氣室。ASBR反應(yīng)器置于礦棉絕熱保溫的恒溫箱內(nèi)，用WMXK-01型溫控儀控制環(huán)境溫度在30℃～35℃范圍內(nèi)。混合液用氣體間歇攪拌，氣體由反應(yīng)器自身生產(chǎn)的生物氣回流。攪拌時(shí)間與間隔時(shí)間為每間歇30 min攪拌3 min，由DHS-S型數(shù)顯時(shí)間繼電器控制。氣體循環(huán)攪拌流量由轉(zhuǎn)子流量計(jì)顯示和控制，約為 300 mL/min。ASBR反應(yīng)器的循環(huán)周期8 h，每個(gè)循環(huán)周期內(nèi)反應(yīng)時(shí)間6 h，每周期的進(jìn)水和排水量均為3 L。

2.2 廢水配制與測(cè)試方法

ASBR反應(yīng)器進(jìn)水以葡萄糖作為碳源，NaHCO3用于調(diào)節(jié)pH值。進(jìn)水的堿度水平由進(jìn)水COD濃度(mg/L)和所投加的NaHCO3產(chǎn)生的堿度濃度(mgCaCO3/L)的比值來(lái)表征，即:進(jìn)水 COD/Alk。本文考察堿度水平COD/Alk分別為 2.0∶1和4.8∶1條件下的運(yùn)行特征。COD濃度測(cè)定采用重鉻酸鉀法。VFA組分與濃度、沼氣組分分析采用氣相色譜法。

3 建模過(guò)程

3.1 建模原理

厭氧硝化過(guò)程可以簡(jiǎn)化為僅涉及兩大類(lèi)微生物:即產(chǎn)酸菌和產(chǎn)甲烷菌的生長(zhǎng)與衰亡過(guò)程，葡萄糖在產(chǎn)酸菌的作用下被轉(zhuǎn)化為揮發(fā)性脂肪酸(VFA)，VFA在產(chǎn)甲烷菌的作用下被轉(zhuǎn)化為CH4和CO2。假設(shè)ASBR反應(yīng)器內(nèi)部的顆粒污泥和氣相完全混合，生物氣主要含有CH4和CO2兩種氣體組分。如果把葡萄糖兩步降解過(guò)程近似看成是兩個(gè)簡(jiǎn)單反應(yīng)的連續(xù)反應(yīng)，則可以表示為:

在這一串聯(lián)反應(yīng)中速率常數(shù)k1≠k2，厭氧處理系統(tǒng)的正常運(yùn)轉(zhuǎn)取決于產(chǎn)酸反應(yīng)速率與產(chǎn)甲烷反應(yīng)速率的相對(duì)平衡。如果A的初始濃度為 CA0，則在時(shí)刻 t的濃度CA，CB及CC間存在下列關(guān)系:

上述關(guān)系可表示為圖1。由串聯(lián)反應(yīng)條件得出下列速率關(guān)系。

由式(3)積分得:

將式(6)代入式(4)積分整理后可得:

將式(7)代入式(5)得:

寫(xiě)成下列積分形式:

進(jìn)行積分:

3.2 求解反應(yīng)速率常數(shù)

堿度COD/Alk為 2.0∶1時(shí)(即中性環(huán)境下)，k1＜k2;堿度COD/Alk為4.8∶1時(shí)(即偏酸性環(huán)境下)，k1＞k2。而且，兩種堿度條件下 k1的差值為0.26，k2的差值僅0.01，說(shuō)明堿度對(duì)產(chǎn)酸速率常數(shù)k1的影響程度大于對(duì)產(chǎn)甲烷速率常數(shù)k2的影響。將k1，k2值代入式(6)、式(7)和式(9)后，得到兩種堿度條件下的反應(yīng)速率數(shù)學(xué)模型，如表1所示。

表1 不同堿度下的反應(yīng)速率方程

3.3 模型精度驗(yàn)證

圖2，圖3分別為VFA濃度(以 CODVFA計(jì))和甲烷產(chǎn)率(以CODCH4計(jì))的實(shí)測(cè)值與數(shù)學(xué)模型計(jì)算值的對(duì)比。由圖看到，所建立的模型計(jì)算值與試驗(yàn)數(shù)據(jù)之間有較好的相關(guān)性，相關(guān)系數(shù) R平均值大于0.95，表明該模型對(duì)ASBR反應(yīng)器處理以葡萄糖為底物的配制廢水中的VFA濃度和甲烷產(chǎn)率CODCH4的模擬具有良好的適用性。

4 結(jié)語(yǔ)

本文建立了不同堿度條件下ASBR反應(yīng)器處理溶解性高濃度有機(jī)廢水的反應(yīng)速率數(shù)學(xué)模型。

當(dāng)進(jìn)水COD/Alk為2.0∶1時(shí)，反應(yīng)器內(nèi)底物、中間產(chǎn)物和最終產(chǎn)物的濃度可用下式表示:

當(dāng)進(jìn)水COD/Alk為4.8∶1時(shí)，反應(yīng)器內(nèi)底物、中間產(chǎn)物和最終產(chǎn)物的濃度可用下式表示:

與其他碳水化合物相比葡萄糖的酸化速率最大，因此該數(shù)學(xué)模型可用于高濃度有機(jī)廢水ASBR處理工藝運(yùn)行的控制和預(yù)測(cè)中。

[1] Randall A W，Dague R R.Enhancement of granulation and start-up in the anaerobic sequencing batch reactor[J].Water Environmental Research，1996，68(5):883-892.

[2] Angenent L T ，Sung S，Raskin L.Methanogenic population dynamics during start-up of a full-scale ASBR treating swine waste[J].Water Research，2002(36):4648-4654.

[3] 孫劍輝，王海燕.ASBR反應(yīng)器快速啟動(dòng)的研究[J].環(huán)境工程，2001(19):10-12.

[4] Masse D I，Masse L，Verville A，et al.The start-up of anaerobic sequencing batch reactors at 20℃and 25℃for the treatment of slaughterhouse wastewater[J].Journal of Chemical Technology&Biotechnology ，2001 ，76(4):393-400.

[5] Wirtz R A，Dague R R.Effect ofinitialseedconcentration and granulation enhancement factors on start-up of the ASBR[J].Adv Filtr Sep Technol，1995(9):654-671.

[6] Ng W J，Ong S L，Tan K Y ，et al.Toxicity assays determine the start-up strategy for an ASBR reactor(an SBR)[J].Water Science Technology，2002 ，46(11):343-348.

[7] Ong S L，Hu J Y，Ng W J，et al.Granulation enhancement in ASBR operation[J].Journal of environmental Engineering，2002，28(4):387-390.

[8] Ioannis S，David M ，Bagley.Improving anaerobic sequencing batch reactor performance by modifying operational parameters[J].Water Research ，2002，36(1):363-367.

[9] 傅大放，靳 強(qiáng)，周培國(guó)，等.厭氧序批式活性污泥工藝特性研究[J].中國(guó)給水排水，2000，16(10):1-5.

[10] Bagley D M，Brodkorb T S.Modeling microbial kinetics in an anaerobic sequencing batch reactor-model development and experimental validation[J].Water Environmental Research，1999，71(7):1320-1332.

[11] Fernandes L，Kennedy K J，Nine Z.Dynamic modeling of substrate degradation in sequencing batch anaerobic reactors(SBAR)[J].Water Research，1993，27(11):1619-1628.

[12] Kennedy K J，Hamoda M F，Droste R L.A study of anaerobic sequencing batch reactors for the treatment of soluble wastes[J].Lewis Publishers，44th Purdue Industrial Waste Conference Proceedings，Printed in U.S.A ，1990(20):1475-1482.

[13] Masse D I，Masse L，CroteauF.The effect of temperature fluctuations on psychrophilic ASBRs treating swine manure[J].Bioresource Technology，2003(89):57-62.

[14] Masse L，Masse D I.Effect of organic，particulate organic，and hydraulic shock loadson anaerobicsequencing batch reactors treating slaughterhouse wastewater at 20 ℃[J].ProcessBiochemistry，2005(40):1225-1232.

[15] Ammary B Y.Treatment of olive millwastewater using anaerobic sequencing batch reactor[J].Desalination，2005(177):157-165.

[16] 任南琪，王愛(ài)杰.厭氧生物技術(shù)原理與應(yīng)用[M].北京:化學(xué)工業(yè)出版社，2004.

[17] 許保玖，龍騰銳.當(dāng)代給水與廢水處理原理[M].北京:高等教育出版社，2000.