泥水自循環(huán)反應(yīng)器的水力流態(tài)特性及生活污水處理效能

陳思遠(yuǎn)，呂錫武，戴喆秦，許 坤

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

隨著國民經(jīng)濟的發(fā)展，小城鎮(zhèn)污水排放在城鎮(zhèn)污水中比重加大，且呈分散型，并且這些地區(qū)沒有配置相關(guān)的市政管網(wǎng)，污水隨意排放，對周圍的人居環(huán)境帶來了不良影響，迫切需要對這些分散式污水進行處理達標(biāo)后排放。針對這種分散式的污水，本文主要探討了一種具有良好除磷脫氮效能、占地和投資省、運行費用低、自動化水平高的泥水自循環(huán)活性污泥工藝[1]。

該工藝通過5個矩形池子厭氧、缺氧、好氧和沉淀交替，可實現(xiàn)一體化除磷脫氮功能。其中有兩個池子交替作為出水池，依靠水流換向流動實現(xiàn)污泥和混合液的循環(huán)流動和自動回流，沒有強制性污泥回流和混合液回流措施[2-5]。泥水自循環(huán)活性污泥工藝是在傳統(tǒng)AO工藝上的改革創(chuàng)新，具有抗沖擊負(fù)荷能力強、節(jié)省能耗、脫氮除磷效率高等優(yōu)勢。

呂錫武教授科研團隊研究發(fā)現(xiàn)[1-5]，泥水自循環(huán)活性污泥工藝為一體化設(shè)置，恒水位運行，在技術(shù)經(jīng)濟方面與常規(guī)工藝有很大的優(yōu)越性：工藝占地面積為0.30～0.60 m2/m3，與氧化溝工藝的0.70～1.0 m2/m3相比，單位面積設(shè)計污水處理能力提高了25%；工藝單位運行成本為0.49元/m3左右，比AAO的0.97元/m3節(jié)約了50%；由于采用智能型的電腦專家控制系統(tǒng)，管理操作人員人數(shù)只有同等規(guī)模常規(guī)污水處理廠的1/2～1/8。此外，泥水自循環(huán)活性污泥工藝為反硝化菌和聚磷菌提供了合適的厭氧/缺氧、好氧的交替環(huán)境，讓反硝化菌和聚磷菌在系統(tǒng)中有競爭優(yōu)勢，使脫氮除磷效率高于常規(guī)工藝，尤其是TP的去除效率較一般活性污泥法提高了15%左右。

泥水自循環(huán)活性污泥工藝相較于常規(guī)工藝，污染物去除率較高、運行費用較低，具有非常好的推廣運用前景。但是，泥水自循環(huán)活性污泥工藝作為一種新工藝，反應(yīng)器內(nèi)部構(gòu)造相對復(fù)雜，采用多點進水會使水力條件更復(fù)雜，而反應(yīng)器的水力流態(tài)特性和混合程度是影響物質(zhì)傳輸、基質(zhì)和微生物的接觸程度、生化性能和反應(yīng)效率的最基本因素[6-8]。鑒于此，為進一步研究反應(yīng)器內(nèi)部水力狀況，提高泥水自循環(huán)工藝反應(yīng)器的運行效率，通過脈沖響應(yīng)法和示蹤試驗，采用停留時間分布法(RTD)[9]，對泥水自循環(huán)反應(yīng)器(AMTR)進行不同水力停留時間(HRT)下清水試驗的水流流態(tài)的研究，探討反應(yīng)器在空載條件下的水力特性，并考察AMTR穩(wěn)定運行情況下對城市生活污水處理的效能。

1 試驗部分

1.1 裝置

試驗裝置如圖1和圖2所示。反應(yīng)器由5 mm厚的有機玻璃制成，以便于觀察反應(yīng)器內(nèi)的狀態(tài)。其外觀尺寸：長×寬×高=490 mm×330 mm×440 mm，有效水深為368 mm，總?cè)莘e為65 L，有效容積為54.4 L。整個反應(yīng)器由5個格室組成，如圖1所示，從水流方向來說，1#、5#池為邊池，2#、3#、4#三池為中間池，兩個邊池和3個中間池體積大小各相同，邊池體積∶中間池體積=3∶2。五池之間水力連通，每池都設(shè)曝氣和攪拌系統(tǒng)，1#、5#兩池設(shè)有出水口和污泥排放口，交替作為攪拌/曝氣、沉淀出水池[1]。污水可進入除 3#池中的任一個，采用連續(xù)進水，周期交替運行。

圖1 泥水自循環(huán)反應(yīng)器結(jié)構(gòu)示意圖Fig.1 Sketch of AMTR Structure

圖2 水力特性試驗裝置圖Fig.2 Equipments of Experiment

1.2 試驗方法

水樣由蠕動泵定量進入AMTR進口處，整套裝置由可編程序控制器(PLC)控制，按照預(yù)先設(shè)定的運行方式，將設(shè)備的運行狀態(tài)轉(zhuǎn)換為相應(yīng)的電流信號來控制進出水電磁閥、排泥電磁閥、曝氣電磁閥、攪拌機、水泵和鼓風(fēng)機的開關(guān)及運轉(zhuǎn)。

AMTR的水力特性試驗通過采用脈沖信號的示蹤劑試驗法進行RTD的測定。以分析純NaCl為示蹤劑，待反應(yīng)器內(nèi)流速穩(wěn)定后，用蠕動泵按當(dāng)前清水流量在進水口處加入50 mL NaCl溶液，此時Cl-注入量約為30 g，再過一定時間間隔，在檢測口測定示蹤劑濃度[10]。

具體試驗步驟：(1)配置一定濃度的示蹤劑溶液(600 g/L)，密封于試劑瓶中；(2)在取樣口(第一、第三、第五隔室)取適量出水，測定該組試驗條件下反應(yīng)器內(nèi)示蹤劑的本底濃度；(3)將示蹤劑溶液用蠕動泵按當(dāng)前清水流量在進水口處加入，同時開始計時；(4)每隔10 min從取樣口取出適量水樣，測定離子濃度(扣除了水相的豐底值)并記下取樣時間，試驗?zāi)┪舶措妼?dǎo)率變化量，延長測樣時間至20 min，整個取樣歷時為反應(yīng)器理論水力停留時間的3倍[11]。

2 RTD標(biāo)準(zhǔn)化曲線的理論分析(流態(tài)試驗參數(shù)計算公式)

2.1 標(biāo)準(zhǔn)化質(zhì)量濃度

標(biāo)準(zhǔn)化質(zhì)量濃度計算如式(1)。

(1)

Ci—t時刻反應(yīng)器出水中示蹤劑的濃度。

2.2 標(biāo)準(zhǔn)化時間

標(biāo)準(zhǔn)化時間計算如式(2)。

θ=t/HRT

(2)

t—示蹤劑離開反應(yīng)器的時刻。

RTD曲線是以標(biāo)準(zhǔn)化質(zhì)量濃度C(θ)為縱軸，以標(biāo)準(zhǔn)化時間θ為橫軸，作出的C(θ)-θ曲線稱為RTD曲線[6]，此曲線所圍面積為1。這樣，對于某種運行條件下的RTD研究，以θ為橫軸、C(θ)為縱軸可繪得一系列的離散點，由這些離散點擬合所得的曲線即為該種條件下的RTD標(biāo)準(zhǔn)化曲線[9-10]。

2.3 平均停留時間

平均停留時間計算如式(3)。

(3)

其中：△t—示蹤劑測定的取樣時間間隔。

2.4 方差

方差計算如式(4)。

(4)

無因次方差計算如式(5)。

(5)

2.5 多級全混流模型

假設(shè)以N個等體積的全混流反應(yīng)器(CSTR)來模擬實際的反應(yīng)器，以反應(yīng)器串聯(lián)級數(shù)N作為模型參數(shù)，如式(6)。

(6)

N越大，串聯(lián)個數(shù)越多，越趨于理想推流流態(tài)；反之，越趨于完全混合流態(tài)。即N=1時為完全混合流態(tài)，N=∞時為理想推流流態(tài)[12]。

2.6 軸向擴散模型

軸向擴散模型在反應(yīng)器存在某種程度的返混時使用，該模型的唯一參數(shù)是D/UL(其中D為軸向分散系數(shù)，U為流體速度，L為特征長度)，其意義為軸向擴散流動與軸向?qū)α髁鲃拥南鄬Υ笮?，計算如?7)。

(7)

軸向擴散模型參數(shù)表示偏移推流的程度，D/UL越大，偏移推流的程度越大，越趨于完全混合，反之亦然[10-14](理想完全混合反應(yīng)器D/UL=∞，而理想平推流D/UL=0)。

2.7 死區(qū)率

死區(qū)率計算如式(8)～式(9)。

Vd/V=1-υaμa

(8)

(9)

其中：Vd—反應(yīng)器死區(qū)，L；

V—反應(yīng)器的有效容積，L；

υa—當(dāng)θ=2時，示蹤劑流出量占注入量的比值；

μa—C-θ曲線的平均值[8]。

3 結(jié)果與討論

3.1 同一HRT對不同工藝運行方式的影響

清水試驗是在反應(yīng)器內(nèi)不加污泥的條件下，采用自來水作為進水，加入示蹤劑進行試驗。首先，考察在同一HRT下，HRT對AMT和AAO工藝的影響。五箱正常運行，進行階段轉(zhuǎn)換；此后，通過改變PLC程控柜設(shè)置，讓反應(yīng)器一直處在階段一狀態(tài)，無階段轉(zhuǎn)化，即為簡易的AAO工藝。試驗數(shù)據(jù)標(biāo)準(zhǔn)化處理后結(jié)果如圖3所示，對應(yīng)的RTD 曲線的特征參數(shù)的計算結(jié)果如表1所示。

圖3 HRT=9.7 h時不同工藝運行方式的RTD曲線Fig.3 RTD Curves of Different Processes under HRT=9.7 h

表1 HRT=9.7 h時不同工藝運行方式RTD曲線數(shù)值分析結(jié)果Tab.1 Effect of Different Processes on RTD Curve under HRT=9.7 h

由圖3可知，AMT時的C(θ)-θ曲線峰值高于AAO，且前者峰值出現(xiàn)的時間較早。AMT時C(θ)-θ曲線不是連續(xù)曲線，主要是因為AMT運行時存在著進出水的交替，在階段交替時測量的數(shù)據(jù)存在一定的間斷。AMT時峰值出現(xiàn)的時間與HRT相差較大，因此AMT時裝置內(nèi)水流并非趨于平推流，且混合程度較大[8,15-16]。

由表1可知，代表反應(yīng)器串聯(lián)級數(shù)的N由AAO運行時的3.49降低到了AMT時的2.87，與之相反，反應(yīng)器混合程度參數(shù)D/UL由AAO時的0.173增加到了AMT時的0.223，死區(qū)率由AAO時的17.4%下降到了AMT時的14.4%。主要是因為,當(dāng)AMT運行時，進出水交替，池內(nèi)水的返混程度較高，水流在各個格室之間縮放流動且連續(xù)不斷，以至形成眾多的小旋渦，死區(qū)率及推流效果相比于呈推流狀態(tài)的AAO工藝均有下降。

3.2 不同HRT對AMTR的影響

清水條件下，加入NaCl示蹤劑，考察HRT為9.7、11.7、13.75 h時對泥水自循環(huán)反應(yīng)器流態(tài)的影響，并進行對比分析。試驗數(shù)據(jù)標(biāo)準(zhǔn)化處理后結(jié)果如圖3所示，對應(yīng)的RTD曲線的特征參數(shù)的計算結(jié)果如表2所示。

圖4 不同HRT時的RTD曲線Fig.4 Curves of Effect of RTD on Different HRT

表2 不同HRT時泥水自循環(huán)反應(yīng)器RTD曲線數(shù)值結(jié)果分析Tab.2 Effect of HRT on RTD Curve in AMTR

由圖4可知，三個HRT條件下，AMTR的RTD曲線均較相似，說明AMTR的結(jié)構(gòu)對標(biāo)準(zhǔn)化后的曲線影響程度一致[10]。各HRT下的C-θ曲線均為不對稱單峰曲線，反映的變化規(guī)律是在峰值之前，示蹤劑濃度C值迅速增大；峰值過后，C值逐漸減小[10]。各HRT條件下，反應(yīng)器出口處示蹤劑的濃度出現(xiàn)峰值均是在t=0.5 HRT左右，且峰高變化不大。根據(jù)眾數(shù)停留時間指數(shù)的定義，即對于平推流，該比值趨近于1；對于完全混合流態(tài)，該比值趨近于0[15]。此外，同一HRT下，單獨觀察5個箱子，發(fā)現(xiàn)每個箱子的D/UL均較大，混合程度較好，單個箱體的空間可視為CSTR，五個箱子的CSTR型單元反應(yīng)器串聯(lián)起來就接近于推流型反應(yīng)器[11]，與眾數(shù)停留時間指數(shù)相符合。因此，AMTR的流態(tài)處在完全混合與平推流之間，具有較強的抗沖擊負(fù)荷能力。

通常，串聯(lián)級數(shù)N隨著HRT的增大而增大，即反應(yīng)器的流態(tài)隨著HRT的增大由完全混合態(tài)向推流態(tài)轉(zhuǎn)變[10,16-17]。由表2可知，隨著HRT的增大，AMTR的串聯(lián)級數(shù)N由2.87降低到了2.27，混合程度指數(shù)D/UL由0.223增加到了0.317，這與常規(guī)反應(yīng)器的變化規(guī)律完全相反[17-19]。主要原因是AMTR階段轉(zhuǎn)換進出水方向，由于進出水方向的不斷交替，反應(yīng)器內(nèi)的水流返混程度劇烈，相同的運行周期條件下，HRT越大，水流的返混次數(shù)越多。因此，隨著HRT的增大，反應(yīng)器完全混合作用增強，推流作用削弱。

3.3 死區(qū)分析

示蹤試驗為清水試驗，反應(yīng)器死區(qū)指水力死區(qū)，是由于反應(yīng)器內(nèi)部結(jié)構(gòu)造成的。由圖5可知，當(dāng)HRT由9.7 h變化到10.5 h再增加到11.7 h時，死區(qū)率由14.4%減少到13.26%再減少到12.4%，呈遞減趨勢。隨后繼續(xù)增大HRT，死區(qū)率反而增大，從12.7 h時的15.06%增加到13.75 h時的16.89%。分析原因主要是，隨著HRT的增大，反應(yīng)器內(nèi)部水流的流速變慢，返混現(xiàn)象減弱，死區(qū)率減小[10,16,19]，此時HRT的變化占主導(dǎo)因素；隨著HRT進一步增大，反應(yīng)器內(nèi)部水流流速的變慢速率降低，而反應(yīng)器交替進水所導(dǎo)致的返混現(xiàn)象程度加深，使得死區(qū)率增高[20]。另外，由于1#、5#池在階段轉(zhuǎn)換中間歇性作為沉淀池，導(dǎo)致整個泥水自循環(huán)工藝死區(qū)率略高于一般反應(yīng)器。

圖5 HRT對Vd/V的影響Fig.5 Impact of HRT on Vd/V

3.4 污泥接種培養(yǎng)(效能驗證)

為了驗證泥水自循環(huán)反應(yīng)器對城市生活污水的實際處理能力，對其進行除磷脫氮工藝的試驗研究，試驗地點為南京市江寧區(qū)科學(xué)園污水處理廠。

由圖6可知，污泥培養(yǎng)時間由2017年9月6日～2017年10月5日，試驗的第1 d就直接將污水處理廠 Orbal氧化溝中的混合液(MLSS=3 310 mg/L，SV=24%)加入到試驗裝置中，采取直接進水培養(yǎng)方式，從沉砂池出口取水，開始3 d間歇進水(進水3 h，停止進水1 h，進水量緩慢增加[9])。之后采取連續(xù)進水，進水流量緩慢增加至第13 d時，流量穩(wěn)定在5 L/h。2017年4月22號開始切換到自動運行狀態(tài)，設(shè)置為主體段2 h、過渡段1 h、沉淀段1 h[1-4]， HRT=11.7 h。

圖6 污泥培養(yǎng)期間污泥沉降比和污泥濃度Fig.6 SV and MLSS during Sludge Cultivation

污泥培養(yǎng)馴化過程中，反應(yīng)器內(nèi)污泥量逐漸增加，在污泥濃度提高的同時，污泥活性增強，絮凝效果改善，污泥沉降比相應(yīng)提高。培養(yǎng)開始階段不排泥，待5個反應(yīng)池平均污泥濃度達到3 000 mg/L時，通過PLC自動排泥，控制污泥齡在12 d左右，維持污泥量穩(wěn)定在3 500 mg/L。

污泥濃度由第1 d的2 101 mg/L，至后期穩(wěn)定為3 500 mg/L；污泥沉降比由第1 d的15%緩慢上升到最后的25%。

隨著污泥量的緩慢增加和混合液中微生物的多樣化，泥水自循環(huán)活性污泥系統(tǒng)的處理效果也不斷提高，COD、TN、TP的去除率均得到明顯改善。對AMTR系統(tǒng)進行了3次污染物去除率的測定，在污泥培養(yǎng)馴化結(jié)束后(即系統(tǒng)后期穩(wěn)定時)，去除率分別在85%、70%和85%左右。如圖7所示，出水COD、TN、TP基本滿足《城鎮(zhèn)污水處理廠污染物排放標(biāo)準(zhǔn)》(GB 18918—2002)一級 A 標(biāo)準(zhǔn)要求，出水效果趨于穩(wěn)定。試驗表明，泥水自循環(huán)反應(yīng)器對城市生活污水的處理效能較高。

圖7 污泥培養(yǎng)期間污染物去除率變化圖Fig.7 Change of Pollutant Removal Rate during Sludge Cultivation

4 結(jié)論

(1)同一HRT，泥水自循環(huán)工藝與AAO工藝相比，具有較低的死區(qū)率以及較高的混合程度，流態(tài)處于完全混合與平推流之間，且具有較強的抗沖擊負(fù)荷能力。

(2)HRT是影響AMTR水力特性的主要因素。隨著水力停留時間的增大，與一般的反應(yīng)器變化規(guī)律相反，AMTR的N值減小，D/UL增大，水流推流程度減弱，混合程度加強。

(3)HRT對泥水自循環(huán)工藝死區(qū)率影響顯著。隨著HRT的增大，死區(qū)率先減小后升高，在HRT=11.7 h時，取得最小死區(qū)率為12.4%，對反應(yīng)器的設(shè)計確定具有指導(dǎo)意義。

(4)泥水自循環(huán)工藝處理城市生活污水效能較高。污泥培養(yǎng)馴化結(jié)束后即系統(tǒng)后期穩(wěn)定運行時，COD、TN、TP的去除率分別為85%、70%和85%，出水穩(wěn)定，且基本滿足《城鎮(zhèn)污水處理廠污染物排放標(biāo)準(zhǔn)》(GB 18918—2002)一級 A 標(biāo)準(zhǔn)要求。