李 鴻, 藺小林*, 李建全, 裴立影

(1.陜西科技大學(xué) 數(shù)學(xué)與數(shù)據(jù)科學(xué)學(xué)院, 陜西 西安 710021; 2.陜西科技大學(xué) 環(huán)境科學(xué)與工程學(xué)院, 陜西 西安 710021)

0 引言

全球性水污染問(wèn)題已是人類(lèi)面臨的主要問(wèn)題之一,它對(duì)人類(lèi)生存和社會(huì)經(jīng)濟(jì)發(fā)展構(gòu)成越來(lái)越嚴(yán)重的威脅.防治水環(huán)境污染,保護(hù)水資源,走可持續(xù)發(fā)展的道路已經(jīng)成為人類(lèi)共同追求的目標(biāo).我國(guó)對(duì)環(huán)境治理也越來(lái)越重視,并把環(huán)境保護(hù)具體計(jì)劃提上日程.由于我國(guó)地表水體污染主要是由工業(yè)廢水和城市生活污水的排放所造成的,污水處理廠的建設(shè)成為改善我國(guó)水環(huán)境的重要措施之一.活性污泥法通過(guò)近百年的發(fā)展已成為城市污水處理廠最為廣泛使用的方法之一[1].

然而,如何合理實(shí)現(xiàn)污水處理廠工藝設(shè)計(jì),優(yōu)化運(yùn)行管理,在保證出水水質(zhì)的條件下減少運(yùn)行費(fèi)用,降低資源的浪費(fèi),提高能源利用率是順應(yīng)我國(guó)當(dāng)前保護(hù)城市環(huán)境、節(jié)能降耗發(fā)展要求的一項(xiàng)重要課題.以數(shù)學(xué)模型為基礎(chǔ)的模擬與控制是實(shí)現(xiàn)這目標(biāo)的有效手段,準(zhǔn)確、貼近實(shí)際的機(jī)理模型又是基礎(chǔ)保證.數(shù)學(xué)模型能夠較全面地描述復(fù)雜系統(tǒng)內(nèi)的多種生物反應(yīng)過(guò)程,不僅可以預(yù)測(cè)污水處理廠的出水情況,同時(shí)也能分析污水處理廠的運(yùn)行狀況,輔助工藝選擇、方案確定、管理與指導(dǎo)污水廠的操作與運(yùn)行,實(shí)現(xiàn)污水處理廠的最優(yōu)化設(shè)計(jì)與節(jié)能運(yùn)行管理.污水處理方法在早期主要是以物理方法去除污水中的雜質(zhì)和污水中主要含有的碳.隨著人們對(duì)污水回收利用要求的不斷提高,污水處理逐漸發(fā)展成為以生物處理為主的二級(jí)處理.在生物處理中,最為普遍采用的方法是活性污泥法.我國(guó)已建和在建的城市污水處理廠中,處理污水所采用的各種類(lèi)型的活性污泥法占90%以上.

活性污泥法經(jīng)過(guò)近百年的實(shí)踐和發(fā)展,已經(jīng)成為國(guó)內(nèi)外應(yīng)用最普遍的污水處理技術(shù)之一.活性污泥數(shù)學(xué)模型的發(fā)展經(jīng)歷了從傳統(tǒng)靜態(tài)模型到動(dòng)態(tài)模型的發(fā)展過(guò)程,如早期Eckenfelder模型、Mc Kinney模型和Lawrence -McCarty模型等靜態(tài)模型,以及Andrews模型、WRc模型和ASM系列模型等動(dòng)態(tài)模型.然而,污水處理系統(tǒng)是一個(gè)多變量、非穩(wěn)定、時(shí)變的復(fù)雜系統(tǒng),如水質(zhì)水量動(dòng)態(tài)變化、反應(yīng)機(jī)理復(fù)雜、存在許多干擾因素和不確定性因素等,因此,以數(shù)學(xué)模型為基礎(chǔ)、結(jié)合污水處理工藝原理的系統(tǒng)仿真方法,成為對(duì)復(fù)雜的污水處理系統(tǒng)進(jìn)行動(dòng)態(tài)分析的適用方法.

活性污泥系統(tǒng)的模型化對(duì)水處理設(shè)計(jì)和研究具有重要的意義.活性污泥系統(tǒng)生物去除有機(jī)碳、氮、磷及其調(diào)控機(jī)理模型研究是當(dāng)今國(guó)際水協(xié)(IWA) 污水生物處理設(shè)計(jì)與運(yùn)行數(shù)學(xué)模型課題組的前沿性基礎(chǔ)課題,該研究成果可以應(yīng)用于污水處理廠運(yùn)行的優(yōu)化和節(jié)能降耗.

國(guó)內(nèi)許多學(xué)者對(duì)污水處理機(jī)理進(jìn)行了深入研究.針對(duì)昆山錦溪污水處理二期工程的改造,周軍[2]通過(guò)控制污泥齡、污泥回流比、好氧工藝中的溶氧區(qū)溶解氧、厭氧區(qū)溶解氧等關(guān)鍵參數(shù),為實(shí)現(xiàn)水質(zhì)達(dá)標(biāo)提供了良好的操作方案;王巍巍等[3]利用GPS-X軟件,基于ASM2d模型,分析了分段進(jìn)水A2/O工藝與常規(guī)A2/O工藝的出水效果,通過(guò)數(shù)學(xué)模型的模擬,確定了最佳配水比對(duì)原有工藝進(jìn)行改造;馬昭等[4]的研究報(bào)道西安市某污水處理廠奧貝爾氧化溝污水處理工藝,利用ASM2d模型建模校正后對(duì)處理工藝進(jìn)行了優(yōu)化分析,優(yōu)化了氧化溝污泥回流比與BOD5污泥負(fù)荷率;郭彥雪等[5]對(duì)國(guó)際水協(xié)推出的活性污泥2D號(hào)模型ASM2d進(jìn)行優(yōu)化,以菌膠團(tuán)的厭氧核為著手點(diǎn),改良ASM2d模型后模擬COD、N和P在氧化溝中的降解過(guò)程,從而改良氧化溝工藝;陳文亮等[6]用ASM2d數(shù)學(xué)模型研究六箱一體化工藝在脫氮除磷中的應(yīng)用,校正后的模型對(duì)出水水質(zhì)的核心影響因素是來(lái)自于碳源不足;關(guān)夢(mèng)龍等[7]結(jié)合COD、OUR數(shù)據(jù)修正ASM3模型,該模型能準(zhǔn)確模擬異養(yǎng)微生物生長(zhǎng)與外源底物之間的關(guān)系,對(duì)于污水的過(guò)程處理分析提供了有效的借鑒方案;隋軍等[8]基于活性污泥模型ASM2d對(duì)SBR工藝的兩步硝化建模,通過(guò)模型準(zhǔn)確獲取了模擬硝化過(guò)程從啟動(dòng)到短程硝化的過(guò)程參數(shù),為模型模擬硝化作用提供了良好借鑒;張學(xué)穩(wěn)等[9]基于BSM平臺(tái),研究了污水處理過(guò)程的溫室氣體排放情況,并用軟件模擬不同環(huán)境下的溫室氣體構(gòu)成,提出了一些建議.以上文獻(xiàn)中,均沒(méi)有考慮用數(shù)學(xué)的方法分析不同條件下曝氣池的穩(wěn)定性情況.

對(duì)于部分學(xué)者提出的污水處理數(shù)學(xué)模型進(jìn)行數(shù)學(xué)方法的研究,可以更深入理解活性污泥污水處理的基本原理和處理機(jī)制,結(jié)合實(shí)際實(shí)驗(yàn)數(shù)據(jù)對(duì)參數(shù)進(jìn)行分析,可以預(yù)測(cè)未來(lái)一段時(shí)間的曝氣池中各個(gè)組分濃度的變化趨勢(shì),從而根據(jù)結(jié)果來(lái)調(diào)整反應(yīng)過(guò)程的參數(shù),更快的達(dá)到平衡態(tài)或調(diào)整平衡點(diǎn)的具體位置.本文基于經(jīng)典微生物和底物的物料守恒數(shù)學(xué)模型,分析了系統(tǒng)平衡點(diǎn)的性態(tài),用數(shù)值模擬來(lái)驗(yàn)證分析所得結(jié)果的正確性,對(duì)系統(tǒng)中參數(shù)的相互關(guān)系進(jìn)行了討論,對(duì)污水處理過(guò)程的優(yōu)化有一定指導(dǎo)意義.

1 物料守恒數(shù)學(xué)模型

圖1表示的是無(wú)回流完全混合污水處理過(guò)程的框圖.

圖1 污水處理基本流程

對(duì)應(yīng)系統(tǒng)的微分方程如下:

(1)

污泥齡(SRT)表示活性污泥在生化反應(yīng)池的平均停留時(shí)間,表達(dá)式為生化反應(yīng)池中的活性污泥總量除以每日排出的活性污泥量.水力停留時(shí)間(HRT)表示生化反應(yīng)池中的污水完全更換一輪所需要的時(shí)間,表達(dá)式為生化反應(yīng)池的容積除以向該池通水的流量.具體表達(dá)式如下:

(2)

考慮到進(jìn)水微生物濃度很低,可以忽略不計(jì)(即令X0=0),假設(shè)進(jìn)水有機(jī)質(zhì)濃度恒定.將公式(2)代入到系統(tǒng)(1)后,則上述方程組可以化簡(jiǎn)為如下系統(tǒng):

(3)

系統(tǒng)(3)中:X表示生化反應(yīng)池中的微生物濃度,S表示生化反應(yīng)池中的底物濃度.

2 平衡點(diǎn)的存在性及其穩(wěn)定性

2.1 平衡點(diǎn)的存在性

(4)

和

(5)

其中,平衡點(diǎn)(4)表示的平衡點(diǎn)為曝氣池中微生物濃度最終歸零時(shí)的污水處理的情況.

下面研究平衡點(diǎn)(5).考慮到模型表述的是無(wú)回流的完全混合情況,有如下結(jié)論[10]:

SRT=HRT

(6)

將公式(6)代入平衡點(diǎn)(5)進(jìn)一步化簡(jiǎn)得到:

(7)

由平衡點(diǎn)(7)得到模型存在正平衡點(diǎn)的條件為

(8)

該條件的實(shí)際含義是:比底物利用率和微生物產(chǎn)率的乘積減去微生物的內(nèi)源呼吸速率要比污泥齡的倒數(shù)大,且初始底物濃度要比平衡狀態(tài)的底物濃度高,此時(shí)系統(tǒng)會(huì)存在正平衡點(diǎn);反之,系統(tǒng)不存在正平衡點(diǎn).

綜合以上討論,容易得到如下結(jié)論:

定理1系統(tǒng)(3)的平衡點(diǎn) (4)一直存在,當(dāng)參數(shù)滿(mǎn)足條件(8)時(shí),系統(tǒng)(3)的平衡點(diǎn)(5)存在.

2.2 穩(wěn)定性分析

系統(tǒng)(3)考慮的是無(wú)回流的完全混合情況,在此情況下,平衡點(diǎn)(7)是平衡點(diǎn)(5)的簡(jiǎn)化形式.下面將對(duì)平衡點(diǎn)(4)、(5)進(jìn)行詳細(xì)分析.

對(duì)于系統(tǒng)(3),求得其雅可比矩陣如下

J=

(9)

2.2.1 平衡點(diǎn)(4)的穩(wěn)定性分析

下面先討論平衡點(diǎn)(4)的穩(wěn)定性.對(duì)于平衡點(diǎn)(4),將該點(diǎn)坐標(biāo)代入雅可比矩陣(9)后得到:

(10)

設(shè)λ1、λ2是雅可比矩陣(10)對(duì)應(yīng)特征方程f(λ)=|λE-J|=0的根,則有:

(11)

當(dāng)正平衡點(diǎn)(5)存在時(shí),將條件(8)與公式(11)聯(lián)立后得到λ1λ2<0,可知特征方程的兩個(gè)特征值符號(hào)情況是一正一負(fù),根據(jù)平衡點(diǎn)穩(wěn)定性定理,此時(shí)平衡點(diǎn)(4)是鞍點(diǎn).當(dāng)初值從X=0開(kāi)始運(yùn)行時(shí),系統(tǒng)會(huì)最終趨于平衡點(diǎn)(4);而當(dāng)初值從X>0開(kāi)始運(yùn)行時(shí),系統(tǒng)最終會(huì)遠(yuǎn)離平衡點(diǎn)(4).

當(dāng)正平衡點(diǎn)(5)不存在時(shí),根據(jù)條件(8)有:

(12)

將公式(12)代入到公式(11)可以得到:

(13)

可知特征方程的兩個(gè)特征值符號(hào)情況是同為負(fù),根據(jù)平衡點(diǎn)穩(wěn)定性定理,此時(shí)平衡點(diǎn)(4)穩(wěn)定,系統(tǒng)最終會(huì)趨于該平衡點(diǎn).

2.2.2 平衡點(diǎn)(5)的穩(wěn)定性分析

要對(duì)正平衡點(diǎn)(5)進(jìn)行分析,首先要保證正平衡點(diǎn)的存在,下面的分析將在滿(mǎn)足條件(8)下進(jìn)行.

設(shè)λ3、λ4是平衡點(diǎn)(5)代入雅可比矩陣(9)后對(duì)應(yīng)特征方程f(λ)=|λE-J|=0的根,則有:

(14)

下面對(duì)λ3、λ4的正負(fù)性進(jìn)行討論.

(15)

當(dāng)平衡點(diǎn)是(5)時(shí),又有:

(16)

將公式(15)、(16)代入公式(14)中,化簡(jiǎn)后有:

(17)

(18)

由公式(17)、(18)得到,此時(shí)系統(tǒng)(3)在正平衡點(diǎn)(5)處的雅可比矩陣的特征值滿(mǎn)足λ3λ4>0,且λ3+λ4<0,因此,兩個(gè)特征值都具有負(fù)實(shí)部,根據(jù)平衡點(diǎn)穩(wěn)定性理論可知,此時(shí)正平衡點(diǎn)(5)是穩(wěn)定的,從而在第一象限內(nèi),系統(tǒng)的解最終會(huì)趨于平衡點(diǎn)(5).

結(jié)合對(duì)平衡點(diǎn)(4)、(5)的分析,有如下結(jié)論:

下面來(lái)模擬兩種情形對(duì)應(yīng)的系統(tǒng)相圖.

可見(jiàn)相圖表示的規(guī)律與理論分析的結(jié)果是一致的,即不存在正平衡點(diǎn),系統(tǒng)從第一象限出發(fā)最終趨于平衡點(diǎn)(4);存在正平衡點(diǎn),系統(tǒng)從第一象限出發(fā)最終趨于平衡點(diǎn)(5).

圖 2 不同情況下系統(tǒng)(3)的相圖

3 數(shù)值模擬

圖3 不同情況下系統(tǒng)(3)的時(shí)序圖

從兩種不同情況下的時(shí)序圖可以得到,若參數(shù)不滿(mǎn)足正平衡點(diǎn)存在的條件,即不存在正平衡點(diǎn),系統(tǒng)(3)最終會(huì)趨于平衡點(diǎn)(4);若參數(shù)滿(mǎn)足正平衡點(diǎn)存在的條件,系統(tǒng)(3)的解最終會(huì)趨于正平衡點(diǎn)(5),即正平衡點(diǎn)(5)是穩(wěn)定的.實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了理論分析所得到的結(jié)論是正確的.

對(duì)實(shí)驗(yàn)數(shù)據(jù)擬合[11],文獻(xiàn)中可以得到如下數(shù)據(jù)SRT=HRT=29 h,k=0.85 d-1,kd=0.3 d-1,S0=150 mg/L,X0=1 500 mg/L,經(jīng)計(jì)算,KS=1.8 gBOD/m3,YT=0.535 gBOD/gBOD,代入到系統(tǒng)(3)得到如下模擬圖,如圖4所示.

圖4 代入實(shí)際數(shù)據(jù)的系統(tǒng)(3)時(shí)序圖

從圖4所示的模擬圖像可知,到達(dá)平衡態(tài)的出水底物濃度在實(shí)驗(yàn)數(shù)據(jù)的3.9～7.7 mg/L范圍內(nèi),平衡時(shí)的出水微生物濃度在實(shí)驗(yàn)數(shù)據(jù)的2.5～8.9 mg/L范圍內(nèi),由此可知系統(tǒng)(3)的描述是可行的.

4 部分參數(shù)對(duì)正平衡點(diǎn)的影響

由于實(shí)際工廠中更關(guān)心正平衡點(diǎn)的穩(wěn)定性和移動(dòng)趨勢(shì),所以接下來(lái)討論各個(gè)參數(shù)對(duì)正平衡點(diǎn)的移動(dòng)影響.

4.1 污泥齡和水力停留時(shí)間

首先考慮污泥齡和水力停留時(shí)間對(duì)平衡點(diǎn)的影響.固定S0=10,KS=5,k=4,YT=0.5,kd=0.1改變污泥齡和水力停留時(shí)間,觀察到系統(tǒng)(3)平衡點(diǎn)的變化過(guò)程[12,13],如圖5所示.

圖5 改變污泥齡系統(tǒng)(3)的平衡點(diǎn)移動(dòng)規(guī)律

圖5表示系統(tǒng)(3)存在正平衡點(diǎn)情況下,改變污泥齡(SRT)對(duì)應(yīng)的平衡點(diǎn)移動(dòng)規(guī)律.對(duì)于微生物,隨著污泥齡的增大,平衡狀態(tài)的微生物濃度會(huì)先變大后逐漸變小;對(duì)于底物,隨著污泥齡的增大,平衡狀態(tài)的底物濃度在逐漸減小靠近零.

4.2 生物反應(yīng)池的溫度

再考慮溫度對(duì)平衡點(diǎn)的影響.污水的溫度會(huì)影響微生物的相關(guān)動(dòng)力學(xué)參數(shù),具體公式可以根據(jù)Arrehenius方程[14],如下:

k(T)=k(20℃)exp[α(T-20℃)]

(19)

式(19)中:k(T)是T溫度下的參數(shù)值;α是參數(shù)的修正系數(shù).可以觀察到,溫度對(duì)微生物的動(dòng)力學(xué)參數(shù)影響,相當(dāng)于給微生物的動(dòng)力學(xué)參數(shù)乘以一個(gè)系數(shù),固定SRT=HRT=4,S0=10,KS=5,k=4,YT=0.5,kd=0.1改變微生物的半飽和常數(shù)、比增長(zhǎng)速率、產(chǎn)率和內(nèi)源呼吸速率的系數(shù),觀察平衡點(diǎn)的移動(dòng)情況,如圖6所示.

圖6 改變溫度的系統(tǒng)(3)平衡點(diǎn)移動(dòng)規(guī)律

從圖6可以得到,溫度的適當(dāng)提高與降低,對(duì)平衡狀態(tài)的微生物濃度影響要比對(duì)平衡狀態(tài)的底物濃度影響要大,總體趨勢(shì)是在一定范圍內(nèi),溫度提高可以提升平衡狀態(tài)的底物和微生物濃度,可以根據(jù)實(shí)際情況選擇調(diào)節(jié)污水處理的溫度.

4.3 污水的溶解氧

好氧過(guò)程的污水中溶解氧的影響可以用Monod方程表示[1]

(20)

公式(20)中:μobs表示好氧微生物的比增長(zhǎng)速率,SO2是反應(yīng)器的溶解氧濃度,KS,O2是氧的飽和常數(shù).

公式(20)表明,溶解氧濃度會(huì)通過(guò)對(duì)微生物的比增長(zhǎng)速率產(chǎn)生影響,進(jìn)而使平衡點(diǎn)移動(dòng).下面固定SRT=HRT=10,S0=10,KS=5,YT=0.5,kd=0.1,KS,O2=1,用公式(20)替換掉系統(tǒng)(3)中的比增長(zhǎng)速率,改變氧氣濃度,觀察其對(duì)平衡點(diǎn)的影響.

經(jīng)計(jì)算,要使得正平衡點(diǎn)存在,需要滿(mǎn)足條件(8),此時(shí)溶解氧濃度需要滿(mǎn)足SO2>1.5,讓溶解氧濃度在該條件下逐漸增大,有如下模擬圖,如圖7所示.

圖7 改變?nèi)芙庋醯南到y(tǒng)(3)平衡點(diǎn)移動(dòng)規(guī)律

從圖7可以觀察到,溶解氧濃度在較低水平時(shí),提高溶解氧濃度可以較大幅度地減少平衡狀態(tài)的底物濃度,同時(shí)較大幅度提升平衡狀態(tài)的微生物濃度;當(dāng)溶解氧濃度較高時(shí),提升溶解氧濃度獲得的收益會(huì)逐漸低于投入的成本,需要根據(jù)實(shí)際情況考慮是否加強(qiáng)曝氣.

5 結(jié)論

在保證正平衡點(diǎn)存在的條件下,要想減少平衡狀態(tài)的底物濃度,可以適當(dāng)提高污泥齡,但需要注意不要使得污泥齡過(guò)高,否則會(huì)導(dǎo)致污泥中的絲狀菌大量繁殖,導(dǎo)致污泥的沉降性降低,對(duì)后續(xù)泥水分離產(chǎn)生影響[15];要想提高平衡狀態(tài)的微生物濃度,可以適當(dāng)提高污泥齡;同時(shí),根據(jù)平衡點(diǎn)(7)的表達(dá)式也可以得到,初始底物濃度的提高也會(huì)使得平衡狀態(tài)的微生物濃度提高.

對(duì)于反應(yīng)池的污水處理溫度,適當(dāng)提高溫度,有利于保持活性污泥的濃度和提高反應(yīng)效率;當(dāng)然,溫度對(duì)達(dá)到平衡狀態(tài)的底物濃度影響較小,考慮到成本等問(wèn)題,也可以視情況不對(duì)生化池的溫度就行調(diào)整.

對(duì)于污水中的溶解氧,由于其可以影響微生物的比增長(zhǎng)速率,使得平衡狀態(tài)的底物濃度和微生物濃度均會(huì)被影響.在低溶解氧情況下,活性污泥反應(yīng)速率受到影響,可以適當(dāng)提高污水的溶解氧,以較大幅度降低平衡時(shí)的底物濃度,并較大幅度提高微生物濃度水平;通過(guò)這種方式對(duì)污水處理過(guò)程的優(yōu)化,會(huì)在溶解氧濃度較高時(shí)收益大幅減少,實(shí)際操作中也是需要結(jié)合成本問(wèn)題進(jìn)行取舍.