汪晨晨，程方奎，呂錫武

(1.東南大學能源與環(huán)境學院，江蘇南京 210096；2.東南大學無錫太湖水環(huán)境工程研究中心，江蘇無錫 214000)

分散式污水處理系統(tǒng)是一種能耗和化學消耗均較小的廢水處理方法[1]。目前對于某些由于地形地勢難以建設統(tǒng)一的污水收集管網的地區(qū)，一體化處理技術因其安裝簡單和占地面積較小成為了較為合適的選擇。但如何控制一體化反應器的運行條件，使其既能達到理想的處理效果，又能夠盡可能地節(jié)約能耗，仍然是一體化反應器實際應用于農村生活污水處理中需要著重考慮的問題。本文從優(yōu)化所設計的側向回流一體化生物膜反應器對氨氮處理效果的運行條件出發(fā)，力求在滿足處理要求的同時做到節(jié)約能耗。

響應曲面法(RSM)是一種基于統(tǒng)計學的，將試驗設計和數(shù)學建模相結合的試驗條件尋優(yōu)方法[2-4]。該方法首先通過設計回歸方程擬合影響因素與響應輸出值之間的函數(shù)關系，然后對數(shù)學模型進行檢驗[5]，評價影響因素及其相互作用，進而得出預測的優(yōu)化響應量及優(yōu)化試驗條件[6]。與傳統(tǒng)的數(shù)理統(tǒng)計方法相比，RSM更為快速簡便，可以高效地確定顯著性影響因素進而優(yōu)化工藝參數(shù)，被廣泛地應用于污水治理領域[7]。

本文從自主設計的側向回流一體化生物膜反應器出發(fā)，優(yōu)化了一體化反應器的構型，在一個反應器內同時實現(xiàn)對污水中有機物、氨氮、總氮和總磷等污染物進行去除。設計的側向回流通道一方面能夠稀釋進水，提高反應器的抗沖擊負荷能力，另一方面也能充分利用回流液中的硝態(tài)氮，在填料區(qū)中的缺氧部分實現(xiàn)反硝化，增強對污染物的去除效果。為了探究反應器的最佳參數(shù)組合，利用RSM建立側向回流一體化生物膜反應器對氨氮的去除率與水力停留時間、填料填充率和氣水比的數(shù)學模型，進而確定各因素對氨氮去除率的影響大小順序。在此基礎上，研究多因素對氨氮去除率的交互影響，從而確定側向回流一體化生物膜反應器的最優(yōu)運行參數(shù)，最后比較最優(yōu)運行參數(shù)下氨氮去除率的模型預測值和實際值，反映模型的精確性。

1 試驗部分

1.1 試驗裝置

采用自主設計的側向回流一體化生物膜反應器，結構如圖1所示。

圖1 側向回流一體化生物膜反應器Fig.1 Lateral Reflux Integrated Biofilm Reactor

反應器采用亞克力材質，有效容積為6.5 L，填料區(qū)裝填的填料為改性聚氨酯海綿填料，相關性質如表1所示。

表1 改性聚氨酯海綿填料性質參數(shù)Tab.1 Property Parameters of Modified Polyurethane Sponge Packings

設計待處理污水從反應器底部進入，在氣流的推動下污水向上流動經過填料區(qū)[8]，污水中的污染物被填料上生長的生物膜吸附轉化降解，污水從而得到凈化[9-10]。由于側向回流通道的存在，填料區(qū)出水一部分回流至反應器底端，與進水充分混合。設計的側向回流通道能夠起到稀釋進水、提高一體化反應器的抗沖擊負荷能力的作用。填料區(qū)上部設置的擋板起到泥水分離的作用，并且使出水不易受氣流紊動。反應器下部設置的泥斗收集填料上脫落下沉的生物膜，通過排泥管定期排放。

1.2 試驗污水和接種污泥來源

試驗模擬污水水質為農村單戶或相鄰連戶排放的黑水經化糞池厭氧處理后與灰水混合后的水質。采用人工配水的方法，以葡萄糖作為碳源，氯化銨作為氮源，磷酸二氫鉀作為磷源，添加一定的碳酸氫鈉補充硝化反應所消耗的堿度以維持正常的反應速率，同時添加適量的微量元素維持微生物的正常新陳代謝[11]。污水中CODCr質量濃度為220～280 mg/L，氨氮質量濃度為30～45 mg/L，溫度為24～27 ℃，pH值為7.75～7.90。試驗接種的污泥提取自無錫某污水廠好氧池，污泥MLSS質量濃度為6 210 mg/L。

1.3 試驗設計

本研究采用Box-Behnken RSM優(yōu)化側向回流一體化生物膜反應器去除氨氮的處理效果。該方法可以用來分析評價指標與影響因素之間的線性和非線性關系[12]，并且不需要運行連續(xù)測試[13]。因此，設計運用Box-Behnken法可以在最少試驗運行次數(shù)的情況下，對影響氨氮去除效果的運行條件進行準確模擬優(yōu)化，還可以分析多因素之間的交互作用[14]?；贐ox-Behnken響應曲面模型設計原理，考慮優(yōu)化側向回流一體化生物膜反應器的實際運行參數(shù)。結合前期試驗結果，選取水力停留時間、填料填充率、氣水比作為考察因素，以氨氮去除率作為模型的響應值。應用Design-Expert 12.0.3.0軟件設計3因素3水平的試驗，試驗設計的各因素和水平如表2所示。

表2 試驗設計各影響因子和水平Tab.2 Influencing Factors and Levels for Experimental Design

根據Box-Behnken中心組合原理，對響應曲面回歸建模，其模型如式(1)[15]。

(1)

其中：R——模型的響應值，即氨氮去除率；

β0——常數(shù)項表示中心點修正值；

βi、βii和βij——線性系數(shù)、平方系數(shù)和交互作用項系數(shù)；

Xi和Xj——不同因素水平變量；

ε——構建模型的殘差。

1.4 分析方法

采用哈希快速消解法測定CODCr(哈希消解器DRB200，哈?？梢姽夥止夤舛扔婦R3900)，納氏試劑分光光度法測定氨氮含量[《水質 氨氮的測定納氏試劑分光光度法》(HJ 535—2009)]，使用便攜式YSI-100pH 測定儀測定pH，YSI550A溶解氧測量儀測定DO含量。

2 結果與討論

2.1 模型建立與方差分析

試驗設計和操作結果如表3所示，共17組試驗，其中有5組為零點試驗。

由表3可知，在不同運行參數(shù)組合下，氨氮的去除率在48.19%～98.63%，分散度較大。利用Design-Expert 12.0.3.0軟件對試驗所得數(shù)據進行多元線性回歸分析，擬合可得氨氮去除率模型方程如式(2)。

表3 試驗設計方案與運行結果Tab.3 Experimental Design and Operation Results

(2)

氨氮去除率的模型預期值與實際值的對比如圖2所示，反映了模型預期值和實際值的契合程度。由圖2可知，氨氮去除率的實際值與模型的預期值較為接近，近似滿足Y=X的關系，說明模型的擬合度較好。

圖2 模型預期值和實際值對比Fig.2 Comparison between Expected Value and Actual Value of the Model

為了進一步判斷模型的顯著性，還需對上述建立的氨氮去除率模型進行方差分析及顯著性檢驗。顯著性判斷依據是P<0.000 1時，影響程度為極其顯著；0.000 10.050 0時，影響程度為不顯著。方差分析和顯著性檢驗結果如表4所示。

表4 方差分析結果Tab.4 Results of Variance Analysis

由表4可知，模型的F值為43.92，且P<0.000 1，表明該模型對氨氮去除率的影響極其顯著，擬合程度較好。

進一步對試驗的統(tǒng)計值進行分析以判斷模型的精密度，試驗的統(tǒng)計值如表5所示。

表5 試驗統(tǒng)計結果Tab.5 Statistics Results of the Experiment

綜上，從方差分析和顯著性檢驗的角度，該模型能夠較為精確地描述側向回流一體化生物膜反應器對氨氮的去除率與反應器所選取的填料填充率、運行時所采用的水力停留時間、氣水比之間的真實關系，可以運用該回歸模型優(yōu)化側向回流一體化生物膜反應器去除水中氨氮的運行條件以及對運行結果進行預測。

2.2 氨氮去除率的響應曲面分析

為了更好地研究水力停留時間X1、填料填充率X2和氣水比X3這3個影響因素及其交互作用對氨氮去除率的影響，建立氨氮去除率模型的響應曲面圖，如圖3～圖5所示。

圖3 水力停留時間與填料填充率對氨氮去除率的影響Fig.3 Effect of Hydraulic Retention Time and Filler Packing Rate on Ammonia Nitrogen Removal Rate

圖4 水力停留時間與氣水比對氨氮去除率的影響Fig.4 Effect of Hydraulic Retention Time and Gas Water Ratio on Ammonia Nitrogen Removal Rate

圖5 填料填充率與氣水比對氨氮去除率的影響Fig.5 Effect of Filler Packing Rate and Gas Water Ratio on Ammonia Nitrogen Removal Rate

由圖3可知，隨著水力停留時間的增加，氨氮的去除率逐漸增加，增加趨勢較為穩(wěn)定。隨著反應器內填料填充率的增加，剛開始氨氮的去除率迅速增加，當填料填充率大于38%時，氨氮去除率的增加趨于平緩。且氨氮去除率隨填料填充率的變化比隨水力停留時間的變化更敏感，兩者的交互作用極其顯著，這與方差分析的結果相一致。

由圖4可知，氨氮的去除率隨著水力停留時間X1和氣水比X3的增大而增大，兩者的作用效果類似，說明水力停留時間和氣水比對側向回流一體化反應器去除氨氮的效果具有協(xié)同作用。

由圖5可知，氨氮去除率隨著氣水比的增加而增大，隨著填料填充率的增加而增大。在填料填充率一定時，氨氮去除率隨著氣水比的變化趨勢近似線性增長，增長速率較快。在氣水比較低時，隨著填料填充率的增加，氨氮去除率增長的速度較快，當氣水比較高時，隨著填料填充率的增加，氨氮去除率的增長較為平緩。當填料填充率為32%、氣水比大于7.5時即可保證氨氮的去除率大于85.00%。

綜合方差分析的結果以及對響應曲面圖的分析可知，3個試驗因素對氨氮去除率的顯著性影響順序為填料填充率>水力停留時間>氣水比，其中填料填充率對氨氮去除率的影響最顯著。當水力停留時間大于6 h、填料填充率大于32%、氣水比大于7.5時，側向回流一體化反應器對氨氮的去除率可以保持在85.00%以上。

2.3 優(yōu)化參數(shù)及模型的驗證

根據Design-Expert 12.0.3.0軟件建立的側向回流一體化生物膜反應器去除氨氮的數(shù)學模型，對模型預測結果分析表明，當氨氮的進水質量濃度為30～45 mg/L時，側向回流一體化生物膜反應器在水力停留時間為6.98 h、填料填充率為39.67%、氣水比為7.73的條件下，預測對氨氮的去除率為97.70%。在此條件下進行了3組平行試驗，實際試驗結果表明氨氮的平均去除率為98.09%，與模型預測值偏差僅為0.39%。因此，該模型對側向回流一體化生物膜反應器去除氨氮的運行條件的優(yōu)化和結果預測較為精確。

3 結論

(1)氨氮是農村分散式污水重要的排放指標之一，基于Box-Behnken RSM構建以水力停留時間、填料填充率和氣水比為影響因素，氨氮去除率為響應值的側向回流一體化生物膜反應器氨氮處理效率預測模型。模型的R2為0.982 6，說明運用模型所得的預期值與實際值擬合程度較好，預測結果準確可靠。

(2)對回歸方程和模型的方差進行分析可知，各影響因素對響應量的影響次序為：填料填充率>水力停留時間>氣水比。各影響因子之間的交互作用均為顯著，其中填料填充率和水力停留時間之間的交互作用最為顯著。

(3)綜合考慮，當水力停留時間為6.98 h、填料填充率為39.67%、氣水比為7.73時，在氨氮初始質量濃度為30～45 mg/L的進水條件下，運用模型預測側向回流一體化生物膜反應器對氨氮的去除率為97.70%。在此條件下進行了3次平行試驗，實際試驗結果表明氨氮的平均去除率為98.09%，與模型預測值偏差值僅為0.39%，說明模型的預測較為精確。