吳 恒,趙婷婷,陳愛玲,張 千
(重慶理工大學(xué) 化學(xué)化工學(xué)院, 重慶 400054)
水體富營養(yǎng)化是發(fā)生范圍最廣、危害程度最大的主要污染問題之一,而地表水體中氮素積累的重要原因之一是城鎮(zhèn)污水廠尾水過量氮排放[1-2]。為降低氮素排放,城鎮(zhèn)污水廠被要求進(jìn)行提標(biāo)改造,尾水水質(zhì)必須達(dá)到一級A排放標(biāo)準(zhǔn)[3-4],但城鎮(zhèn)生活污水低碳化趨勢給污水廠的提標(biāo)改造帶來了巨大壓力[5-6]。目前,為了應(yīng)對低碳源條件下總氮去除不達(dá)標(biāo)的問題,大多數(shù)城鎮(zhèn)污水廠采取增設(shè)深度處理工藝的方式,如生物濾池工藝,但該工藝在處理低碳污水時(shí)存在反硝化過程缺少電子供體的問題,導(dǎo)致脫氮效果差[7-8]。為此,城鎮(zhèn)污水廠常采用外加碳源的方式對工藝補(bǔ)充碳源,雖然保證了脫氮效果,但易造成二次污染、運(yùn)行成本高[9-10],對適用于低碳源廢水處理的脫氮新技術(shù)的研發(fā)顯得尤為急迫。
鑒于此,研究以PBS作為BDPs,構(gòu)建了SPD脫氮工藝,通過序批式實(shí)驗(yàn)探究了該工藝單位生物膜硝化、反硝化性能,同時(shí),通過非線性擬合獲得了碳源靜態(tài)釋放動力學(xué)方程。此外,在單級經(jīng)驗(yàn)?zāi)P突A(chǔ)上推導(dǎo)了多級經(jīng)驗(yàn)?zāi)P捅磉_(dá)式,并將實(shí)測數(shù)據(jù)與預(yù)測數(shù)據(jù)進(jìn)行了相對誤差計(jì)算以判斷其準(zhǔn)確性,以期為SPD脫氮工藝的工程應(yīng)用提供實(shí)踐與理論基礎(chǔ)。
如圖1(a)所示,SPD反應(yīng)器由厚度為5 mm的有機(jī)玻璃材料構(gòu)成,高20 cm,直徑18 cm。反應(yīng)器共分2層,外層水浴保溫層,保持反應(yīng)器溫度恒定;內(nèi)層內(nèi)徑8 cm,為反應(yīng)器主體部分,PBS固體碳源填充率30%,廢水有效體積為0.58 L,底部安裝微孔曝氣盤提供氧氣。 如圖1(b)所示,浸出性能實(shí)驗(yàn)裝置采用密閉磨口錐形瓶,短管與彎長管分別用于取樣與曝氮?dú)狻?/p>
圖1 SPD工藝實(shí)驗(yàn)裝置示意圖
PBS購自深圳光華偉業(yè)公司,具體參數(shù)如表1所示。
表1 碳源填料參數(shù)
(1)
(2)
浸出性能實(shí)驗(yàn):在溫度為25~30 ℃條件下,取250 mL磨口錐形瓶進(jìn)行PBS碳源的浸出性研究,取10.00 g PBS放置于加有100 mL蒸餾水的錐形瓶中,排空后密封,置于暗處,每天換蒸餾水,第1 d測定COD(chemical oxygen demand,COD)值,此后每隔5 d測1次,連續(xù)測52 d[19]。
圖2 脫氮動力學(xué)測定曲線
探究了PBS在靜態(tài)條件下的COD的釋放規(guī)律,由圖3可以看出,PBS第1 d COD釋放量較低,第6 d COD釋放達(dá)到峰值295.83 mg/L,這可能是此類聚合材料被浸泡后逐漸釋放COD的緩釋達(dá)峰過程,可見PBS的緩釋達(dá)峰過程為6 d,第16 d后COD釋放趨勢逐漸平穩(wěn),第41~52 d,碳源釋放濃度穩(wěn)定在10~20 mg/L之間,第52 d僅為12.50 mg/L。有研究報(bào)道碳源在有微生物的條件下可以被靶向降解[16-17],而此條件下采用蒸餾水進(jìn)行實(shí)驗(yàn),無微生物干擾,因此,碳源釋放濃度逐漸降低且趨于平緩。但PBS在初期COD釋放波動較大,說明在工藝啟動初期可能會造成有機(jī)物超標(biāo)的風(fēng)險(xiǎn)。需要注意的是,初期的過量碳源釋放可以通過增加曝氣量或延長水力停留時(shí)間來提高碳源利用率和COD去除率[25]。將第6~52 d的碳源浸出數(shù)據(jù)進(jìn)行非線性擬合,可以發(fā)現(xiàn)PBS的碳源浸出規(guī)律符合對數(shù)方程y=-72.15ln(x)+273.02,說明在啟動6 d以后,PBS在不同時(shí)間釋放的COD濃度可以進(jìn)行預(yù)判,對于基于PBS的SPD工藝啟動過程具有一定的指導(dǎo)作用,以防止工藝在啟動初期存在出水有機(jī)物超標(biāo)的風(fēng)險(xiǎn)。因此,SPD工藝啟動初期需要采取措施管控出水有機(jī)物濃度,而6 d以后可以通過模型進(jìn)行COD釋放預(yù)判。
圖3 碳源靜態(tài)浸出性能曲線
米門方程常用于生物膜反應(yīng)器基質(zhì)降解過程的推導(dǎo),本文也根據(jù)該方程進(jìn)行模型推導(dǎo),在實(shí)驗(yàn)數(shù)據(jù)基礎(chǔ)上計(jì)算模型參數(shù),從而獲得經(jīng)驗(yàn)?zāi)P?。齊勇等[26]采用米門方程進(jìn)行了推導(dǎo),認(rèn)為在反應(yīng)器內(nèi)微生物系統(tǒng)對機(jī)智的降解過程符合一次反應(yīng)關(guān)系,其生物反硝化動力學(xué)表達(dá)式為:
V=Vmax·Se/(Ks+Se)
(3)
將V=-ds/dt代入式(3)得:
-ds/dt=Vmax·S/(Ks+Se)
(4)
-ds/dt=Vmax·S/Ks
(5)
令K=Vmax/Ks,
得到
-ds/dt=K·S
(6)
當(dāng)水力混合特性為推流時(shí),積分得:
Se=S0exp(-Kt)
(7)
2) 反應(yīng)器的特性接近推流。本研究采用的小試在時(shí)間上符合推流模式,符合要求;
對式(7)兩邊求對數(shù)得:
ln(Se/S0)=-Kt
(8)
根據(jù)ln(Se/S0)對t作圖,見圖4,得出K=0.08,因此可以得出:
S1=S0exp(-0.08t)
(9)
圖4 ln(Se/S0)與t之間的關(guān)系圖
為了對本單級經(jīng)驗(yàn)?zāi)P瓦M(jìn)行拓展,結(jié)合了相關(guān)的理論基礎(chǔ)進(jìn)行推導(dǎo)。徐斌等[27]對連續(xù)流的生物膜反應(yīng)器進(jìn)行了經(jīng)驗(yàn)?zāi)P屯茖?dǎo),本實(shí)驗(yàn)在其基礎(chǔ)上進(jìn)行了調(diào)整。在本實(shí)驗(yàn)中,不考慮反應(yīng)器內(nèi)的水動力學(xué)變化,僅從宏觀上以連續(xù)運(yùn)行的方式進(jìn)行考慮,可以得到方程(10),其中Q為進(jìn)水流量;A1和A2分別為第1級、第2級反應(yīng)器體內(nèi)填料總體積;S0、S1和S2分別為原水、第1級和第2級反應(yīng)器內(nèi)的硝酸鹽濃度,V1和V2分別為第1級、第2級反應(yīng)器的體積。
QS0-V1A1=QS1
(10)
QS1-V2A2=QS2
(11)
對于多級串聯(lián)反應(yīng)器,其中的K級反應(yīng)動力學(xué)方程為:
QSk-1-VkAk=QSk
(12)
其中:Vk需要進(jìn)行推導(dǎo)。
實(shí)際上,根據(jù)式(9)可以進(jìn)行以下變形:
(S1-S0)/t=[S0exp(-0.08t)-S0]/t
可以得出:
V=[S0exp(-0.08t)-S0]/t
(13)
QSk-1-[S0exp(-0.08t)-S0]/tAk=QSk
(14)
S1=(QS0-V1A1)/Q
S2=(QS0-V1A1-V2A2)/Q
S3=(QS0-V1A1-V2A2)/Q-V3A3/Q=
(QS0-V1A1-V2A2-V3A3)/Q
…
即為多級模型的具體表達(dá)式:
Sk=(QS0-V1A1-V2A2-…-VkAk)/Q
(15)
因此,根據(jù)上述可以整理得出基于米門方程的SPD工藝多級經(jīng)驗(yàn)?zāi)P蜑镼Sk-1-[S0exp(-0.08t)-S0]/tAk=QSk,具體表達(dá)式為Sk=(QS0-V1A1-V2A2-…-VkAk)/Q。
圖5 模型預(yù)測實(shí)測對比圖
2)在SPD工藝啟動初期需要采取措施以管控碳源出水濃度,如增加曝氣、延長水力停留時(shí)間等。碳源釋放峰值期后的碳源浸出規(guī)律符合對數(shù)方程y=-72.15ln(x)+273.02,可以據(jù)此預(yù)判碳源釋放情況。
3) 基于米門方程推導(dǎo)出SPD工藝的單級經(jīng)驗(yàn)?zāi)P蜑镾1=S0exp(-0.08t),而多級經(jīng)驗(yàn)?zāi)P捅磉_(dá)式為QSk-1-[S0exp(-0.08t)-S0]/tAk=QSk。
4) 該模型最大誤差為13.80%,隨著運(yùn)行時(shí)間越久會出現(xiàn)2個(gè)趨勢:誤差隨運(yùn)行時(shí)間的變化可能越大,由正誤差轉(zhuǎn)為負(fù)誤差。模型在長時(shí)間預(yù)測方面需要進(jìn)行改進(jìn),但是對于SPD工藝的應(yīng)用具有指導(dǎo)意義。