王 亞 軍

(1.蘭州理工大學(xué) 土木工程學(xué)院, 甘肅 蘭州 730050; 2.東南大學(xué) 土木工程學(xué)院, 江蘇 南京 210096)

生物過濾系統(tǒng)中多級混合顆粒群滲流過程分析與完善

王 亞 軍1,2

(1.蘭州理工大學(xué) 土木工程學(xué)院, 甘肅 蘭州 730050; 2.東南大學(xué) 土木工程學(xué)院, 江蘇 南京 210096)

[目的] 準(zhǔn)確掌握過濾滲流場的變化規(guī)律，為更深刻理解生物過濾系統(tǒng)中流態(tài)的變化提供理論參考。[方法] 利用顆粒群加壓膨脹理論，基于對流固耦合效應(yīng)和滲流特性進行研究。[結(jié)果] 提出了生物過濾系統(tǒng)滲流過程的4階段假設(shè)，并通過試驗驗證了變化規(guī)律。[結(jié)論] 無論在清水和污水滲流過程中都出現(xiàn)的滲透系數(shù)由下降到上升再下降變化趨勢，是由于顆粒群加壓膨脹產(chǎn)生的原因。

生態(tài)過濾系統(tǒng); 滲流過程; 流固耦合; 飽和滲流; 滲透系數(shù)

文獻參數(shù)： 王亞軍.生物過濾系統(tǒng)中多級混合顆粒群滲流過程分析與完善[J].水土保持通報，2017,37(3)：329-333.DOI:10.13961/j.cnki.stbctb.2017.03.056； Wang Yajun. Analysis and improvement of mixed particle swarm percolation process in biofiltration system[J]. Bulletin of Soil and Water Conservation, 2017,37(3)：329-333.DOI:10.13961/j.cnki.stbctb.2017.03.056

由于實際反應(yīng)器內(nèi)部水流中顆粒物質(zhì)與基質(zhì)之間發(fā)生攔截、擴散、黏附等作用所導(dǎo)致的水流遲滯性，以及植物根系、微生物、外部降水和蒸發(fā)作用等因素對水流狀態(tài)的影響，使得實際反應(yīng)器流態(tài)偏離經(jīng)典理想反應(yīng)器的程度很大，偏離理想流型的原因主要是在反應(yīng)器中存在的死角、短路流及溝流等局部堵塞現(xiàn)象[1]。

生態(tài)過濾池的堵塞現(xiàn)象在小試、中試和現(xiàn)場使用中都頻繁出現(xiàn)[2-4]。過濾系統(tǒng)堵塞問題已成為制約生態(tài)過濾技術(shù)應(yīng)用和推廣的主要障礙，它解決的好壞直接決定著工程的使用壽命和應(yīng)用效益。

在下流式生態(tài)過濾池中，污染物由上及下依次流經(jīng)非飽和好氧過濾層、過渡層和飽和含水厭氧淹沒層，在底部區(qū)域延長水力停留時間來提高微生物濃度，達到提高反硝化處理效果的目的。在這里，出現(xiàn)了兩大問題，一是在底部反應(yīng)后端的反硝化階段，電子供體成為限速因子，蒙納仕大學(xué)研究團隊通過多年研究表明，可以通過在底部淹沒區(qū)添加碳源等方式彌補這個不足[5-7]；二是淹沒區(qū)的厭氧環(huán)境可能加速了系統(tǒng)的滲透堵塞現(xiàn)象，這在其他研究工作中有所闡述[8-9]。但是針對生態(tài)過濾系統(tǒng)的滲流過程研究還很缺乏，沒有完全掌握滲流場變化規(guī)律，同時增加了對水質(zhì)處理的不可靠性。所以，本文擬對滲濾過程理論加以分析與完善，并以試驗進行驗證，以期為相關(guān)研究提供參考。

1 試驗裝置參數(shù)

生態(tài)過濾系統(tǒng)裝置箱體材質(zhì)為無色透明有機玻璃，內(nèi)部尺寸為500 mm×300 mm×700 mm(長×寬×高)，箱體上部布置了一根豐型進水管均勻布水，裝置內(nèi)填料為河砂和礫石，按照澳大利亞FAWB標(biāo)準(zhǔn)指導(dǎo)方法要求的比例裝填[10]，填充高度為550 mm，各層參數(shù)見表1。采用恒流泵(Longer Pump BT100-2 J)保持裝置恒定進水，U型液體壓力計(BYY型)測定濾層內(nèi)部壓差，TST-70型滲透儀測定填料滲透系數(shù)。

表1 生態(tài)過濾池模擬參數(shù)

2 傳統(tǒng)過濾過程理論分析及完善

顆粒以自然狀態(tài)堆積，顆粒間會留有較多孔隙，若經(jīng)過適當(dāng)振實[11]，使得顆粒間的孔隙就會大為減少。但顆粒這種有序排列會受到顆粒尺寸、形狀及表面特性等多種因素干擾，不易達到最理想的有序排列(即最緊密的堆積)。為了探討多級混合不規(guī)則顆粒群填料堆積體的基本特性以及填充顆粒間隙填充物質(zhì)的數(shù)量和作用等基本規(guī)律，首先以球狀顆粒材料為基礎(chǔ)，作基礎(chǔ)計算式的推導(dǎo)和演算。

2.1 顆粒群加壓膨脹特性

宏觀現(xiàn)象源自內(nèi)部的力位形和幾何構(gòu)形的調(diào)整。Jacques Duran[11]通過二維顆粒球體模型闡明了顆粒間空隙部分在外界作用下發(fā)生變化的規(guī)律(如圖1所示)[11]。

圖1 顆粒的加壓膨脹特性

可以看出，當(dāng)上部擠壓顆粒時，則顆粒間空隙的面積S隨顆粒間距呈現(xiàn)拋物線趨勢，面積S存在一個最大值。當(dāng)顆粒堆積密度很大(即面積S小)時，擠壓可使面積S增大，即使得顆粒堆積密度降低。若顆粒堆積密度低到一定程度，則外加作用可造成顆粒密度增大。任何顆粒流動的發(fā)生都是以局部顆粒堆積密度的降低(體積膨脹)為條件的。

2.2 Reschke理論計算方法

對于宏觀顆粒，經(jīng)典力學(xué)可以給出單個顆粒運動狀態(tài)的精確解。但是，大量顆粒組成的體系具有特別的性質(zhì)和運動規(guī)律。顆粒群中顆粒無序分布。針對過濾系統(tǒng)內(nèi)，由不同種類的顆粒混合而成的填料，粒度分布不一定遵守特定的數(shù)學(xué)方程，由多組分組成的填料層，顆粒分布也是不連續(xù)的，不能用傳統(tǒng)的數(shù)學(xué)函數(shù)方程來表達。所以需要尋找出一種堆積密度數(shù)學(xué)模型不受某一固定顆粒分布方程的局限。Reschke理論計算方法[12]在Schwanda模型基礎(chǔ)上擴充的功能，此模型不受顆粒分布函數(shù)關(guān)系的限制，適用于粒度分布不連續(xù)的堆積密度的理論計算。這個數(shù)學(xué)模型除包括了顆粒分布函數(shù)關(guān)系外，還考慮了顆粒形狀和表面粗糙度、細顆粒間干擾作用等因素。

2.2.1 單一粒度理想球狀顆?？紫扼w積份數(shù)K0K0即為單一粒度顆粒間的孔隙體積數(shù)與固體物質(zhì)體積數(shù)的比值。

(1)

式中：K0——單一粒度顆?？紫扼w積份數(shù)；ε——孔隙體積數(shù)(%)；D——堆積密度(與固體物質(zhì)含量成比例)。

2.2.2 顆粒干擾作用寬度W在顆粒干擾作用比較明顯的混合填料體系中，兩粒級的粒度比例范圍稱為顆粒干擾作用寬度，此值取決于基礎(chǔ)顆粒粒度與極限顆粒粒度之比。極限顆粒粒度是指剛好沒有顆粒干擾作用的粒度。

(2)

式中：W——顆粒干擾作用寬度；Xs——基礎(chǔ)顆粒粒度；Xw——極限顆粒粒度。

2.2.3 多級混合顆?？紫扼w積份數(shù)h計算基本方程式以由粗顆粒和細顆粒兩組分組成的混合顆粒群為基礎(chǔ)，Reschke[12]認為，根據(jù)粒級之間粒度比例的不同，多級混合顆粒的孔隙體積份數(shù)如下：

(3)

式中：hs——基礎(chǔ)顆粒級的孔隙體積份數(shù)；K0——單一粒徑粉體孔隙體積份數(shù)；as,i——i粒級的顆粒干擾系數(shù)；Vi——i粒級所占整個固體體積份數(shù)；n——粒級組分數(shù)。

編制數(shù)學(xué)模型的基礎(chǔ)即是K，W，Xi/Xs和系數(shù)as,i之間的數(shù)學(xué)關(guān)系，可用圖2表示[13]。

注：K0為單一粒徑粉體孔隙體積份數(shù)；as,i為i粒級的顆粒干擾系數(shù)；Xi為i粒級的顆粒粒度；Xs為基礎(chǔ)顆粒粒度；W0為單一粒的顆粒干擾作用寬度；W為顆粒干擾作用寬度。

圖2 顆粒干擾系數(shù)as,i與Xi/Xs的關(guān)系

圖中以Xi=Xs為分界點，左側(cè)為XiXs，當(dāng)+W0的對應(yīng)as,i值等于K0時即為第2種類型。

當(dāng)某一粒級作為基礎(chǔ)顆粒時其他各粒級均為填充顆粒，由基礎(chǔ)顆粒的孔隙體積份數(shù)K0減去其他各級填充顆粒的相應(yīng)固體體積份數(shù)Vi，即可得出基礎(chǔ)顆粒級的孔隙體積份數(shù)hs值，如式(3)所示。

2.3 過濾滲流過程機理分析與完善

2.3.1 過濾滲流過程理論 早在1958—1973年敏茨和艾弗斯位學(xué)者就濾料層的過濾過程進行了長期爭論[14]，且至今沒有統(tǒng)一的定論；在這期間，1965年Jones等[8], Kehat等[15]提出了在好氧條件下基質(zhì)堵塞過程經(jīng)歷了3個階段，第1個階段是滲透速率接近系統(tǒng)開始運行時的滲透速率，但呈現(xiàn)逐漸下降趨勢；第2階段是一個實質(zhì)性的平穩(wěn)下降階段；第3階段是直至持續(xù)堵塞發(fā)生。而在厭氧條件下，沒有以上的第2階段，第1和第3階段也不是很明顯辨別。

2.3.2 過濾滲流過程假設(shè) 經(jīng)典滲流力學(xué)一般假定流體流經(jīng)的固體介質(zhì)是完全剛性的，即在孔隙流體壓力變化過程中，固體骨架不產(chǎn)生任何彈性或者塑性變形。但是，在實際的滲流過程中，由于孔隙流體壓力的變化，一方面要引起多孔介質(zhì)骨架有效應(yīng)力變化，由此導(dǎo)致滲透率的變化，另一方面，這些變化又反過來影響孔隙流體的流動和壓力的重新分布，這就必須讓研究者考慮孔隙流體在介質(zhì)中的流動規(guī)律的同時，也要考慮滲流場對介質(zhì)本身變形造成的影響。所以，僅從污染介質(zhì)干擾濾料層角度來分析水體流經(jīng)填料的滲濾過程是不充分的，作者認為還需要考慮運行過程中孔隙的三維幾何形態(tài)變化。為此本文對生物過濾系統(tǒng)堵塞過程進行了重新詮釋。提出2種情況的假設(shè)：

第一，清水滲流過程經(jīng)歷4個階段，第1階段是滲透系數(shù)逐漸下降趨勢；第2階段是滲透系數(shù)有急劇上升勢頭；第3階段是滲透系數(shù)有稍微下降趨勢；第4階段滲透系數(shù)平穩(wěn)不再發(fā)生大的波動變化。

第二，污水滲流過程經(jīng)歷4個階段，第1階段是滲透系數(shù)逐漸下降趨勢；第2階段是滲透系數(shù)有急劇上升勢頭；第3階段是滲透系數(shù)開始遲緩下降趨勢；第4階段是滲透系數(shù)持續(xù)下降，逐漸進入堵塞階段。可能的原因：第1和第2階段，顆粒群緊密堆積下受到水壓力作用導(dǎo)致空間構(gòu)造改變，孔隙壓縮，同時污水?dāng)y帶的顆粒物通過截留、沉淀、擴散和橫向遷移吸附在濾料表層，微生物為生長掛膜階段，生長量大于脫落量；第3階段，由于顆粒加壓膨脹特性，顆?？臻g出現(xiàn)由小變大的趨勢，隨后又擠壓變小，這期間受到外加水力剪切和顆粒摩擦等的作用，濾料表層的顆粒物和生物膜發(fā)生脫落；第4階段，顆粒構(gòu)造基本穩(wěn)定，隨濾料層深度的加大，脫落下的顆粒物和生物膜持續(xù)累積，在濾料層局部不同程度地形成不透水區(qū)域。

3 結(jié)果與分析

3.1 混合填料孔隙率計算

首先對反應(yīng)器運行初期混合填料的孔隙率進行理論計算，過濾系統(tǒng)內(nèi)填料是由8種不同粒級組成的顆粒群混合而成，在研究運行過程中多級混合填料顆粒間孔隙體積的變化，首先需較為科學(xué)地計算出混合填料的孔隙體積，孔隙體積間接表明了水流孔道的變化情況。

3.1.1 多級混合顆?？紫扼w積份數(shù)計算 各平均粒度X分別為0.25，0.50，1.00，2.00，3.00，6.00，12.00，15.00 mm；試驗測得各粒度的孔隙率ε分別為29.5%，31.5%，32%，33.6%，35.8%，37%，43%，44.5%；各粒級孔隙體積分數(shù)K0,i根據(jù)公式(1)計算所得分別為41.84%，45.99%，47.06%，50.60%，55.76%，58.73%，75.44%，80.18%；顆粒干擾作用寬度W0為0.9。以每一粒級都分別作一次基礎(chǔ)顆粒，計算出相應(yīng)的as,i與Vi的乘積，結(jié)果如表2所示。

表2 孔隙體積份數(shù)計算模型

3.1.2 混合填料孔隙率計算 從表2中找出基礎(chǔ)顆?？紫扼w積份數(shù)最大(hs值)的粒級作為混合粉體堆積密度D的計算依據(jù)，可以看出，粒級3的hs最大，為0.161，且有：

則得到混合填料的堆積密度為0.861 3。最終得出混合填料的孔隙率為13.87%。通過表觀密度和堆積密度的測定計算出混合填料的孔隙率為14.54%，相對誤差為4.61%，因此說明按以上方法計算出混合填料孔隙率是較為精確的。

3.2 滲透特性測定

兩套生態(tài)過濾池同步進行清水和污水滲透試驗，在1.0 m3/(m2·d)進水負荷水平下連續(xù)運行。得出相同水力負荷條件下滲透系數(shù)在時間序列上的變化情況，如圖3所示。

圖3 滲透系數(shù)-時間的關(guān)系曲線

在運行初期，清水和污水通過介質(zhì)的滲透系數(shù)后表現(xiàn)出總的下降趨勢；這是由于在試驗裝置持續(xù)運行導(dǎo)致原有顆粒群間的構(gòu)造發(fā)生變化，產(chǎn)生顆粒群骨架流變和局部孔隙塌陷的現(xiàn)象，新形成的空隙比原來的空隙要小，水流在其中流動水頭損失增大，因此滲透系數(shù)減小。隨時間推移，清水和污水流經(jīng)介質(zhì)的滲透系數(shù)均表現(xiàn)出波動上升，這是因為顆粒間空隙變大的區(qū)域多于塌陷區(qū)域，水流通道變寬，滲透系數(shù)呈現(xiàn)曲折爬上趨勢。在長期連續(xù)進水后，兩者出現(xiàn)了不同的變化趨勢，清水流經(jīng)介質(zhì)時，滲透系數(shù)先略有下降后平穩(wěn)，這是由于顆粒間孔隙基本全部塌陷，滲透系數(shù)下降至穩(wěn)定，當(dāng)水頭壓力不足以使顆粒介質(zhì)之間的結(jié)構(gòu)發(fā)生新變化，滲透水流在其中流動的路徑變得較為穩(wěn)定，滲透系數(shù)較為平穩(wěn)。污水流經(jīng)介質(zhì)時，滲透系數(shù)逐漸下降，這是由于隨運行時間的延長，根據(jù)熱力學(xué)基本原理，一個封閉體系總會趨于自由能最小的狀態(tài)，最終都會通過擴散等熱運動，使系統(tǒng)的熵達到最大，形成均勻分布的平衡態(tài)，顆粒構(gòu)造最終穩(wěn)定。此時，污水?dāng)y帶的顆粒污染物和逐漸成熟的生物膜成為制約滲透的限制因子，致使產(chǎn)生堵塞現(xiàn)象發(fā)生。這也支撐了過濾堵塞過程假設(shè)的成立。

4 結(jié) 論

(1) 借助顆粒群加壓膨脹理論，基于對流固耦合效應(yīng)的思考和滲流特性的研究，提出了生態(tài)過濾系統(tǒng)滲流過程的4階段假設(shè)：清水滲流經(jīng)歷滲透系數(shù)下降、上升、下降和平穩(wěn)4階段過程；污水滲流經(jīng)歷滲透系數(shù)下降、上升、平穩(wěn)和下降4階段過程。

(2) 借助Reschke理論計算方法對多級混合填料孔隙率進行了理論計算，并與實測數(shù)據(jù)對比，相對誤差為4.61%，因此說明按此方法計算出混合填料孔隙率較為精確。

(3) 通過連續(xù)運行下生態(tài)過濾池滲透系數(shù)的測定，其變化規(guī)律與假設(shè)滲流過程相同，并通過試驗數(shù)據(jù)支撐了過濾堵塞過程假設(shè)的成立。

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

[2] Coustumer S L, Fletcher T D, Deletic A, et al. The influence of design parameters on clogging of stormwater biofilters: A large-scale column study [J]. Water Research, 2012,46(20):6743-6752.

[3] Kandra H, Mccarthy D, Fletcher T D, et al. Assessment of clogging phenomena in granular filter media used for stormwater treatment [J]. Journal of Hydrology, 2014,512(4):518-527.

[4] Coustumer S L, Fletcher T D, Deletic A, et al. Hydraulic performance of biofilter systems for stormwater management: Influences of design and operation[J]. Journal of Hydrology, 2009,376(1/2):16-23.

[5] Zinger Y, Blecken G T, Fletcher T D, et al. Optimising nitrogen removal in existing stormwater biofilters: Benefits and tradeoffs of a retrofitted saturated zone[J]. Ecological Engineering, 2013,51(2):75-82.

[6] Glaister B J, Fletcher T D, Cook P L, et al. Co-optimisation of phosphorus and nitrogen removal in stormwater biofilters: The role of filter media, vegetation and saturated zone [J]. Water Science & Technology, 2014,69(9):1961-1969.

[7] Payne E G, Fletcher T D, Russell D G, et al. Temporary storage or permanent removal? The division of nitrogen between biotic assimilation and denitrification in stormwater biofiltration systems[J]. PloS one, 2014,9(3):1-12.

[8] Jones J. H. and Taylor G. S. Septic tank effluent percolation through sands under laboratory conditions[J].Soil Science,1965,99(5):310-309.

[9] Blazejewski R, Muratalazejewska S. Soil clogging phenomena in constructed wetlands with surface flow[J]. Water Science and Technology, 1997,35(5):183-188.

[10] Payne E G I, Hatt B E, Deletic A, et al. Adoption Guidelines for Stormwater Biofiltration Systems[M]. Melbourne, Australia: Cooperative Research Centre for Water Sensitive Cities, 2015.

[11] Duran J. Sands, Powders, and Grains: An Introduction to the Physics of Granular Materials[M]. New York: Springer, 2000.

[12] Reschke T. Schriftenreihe der zementindustrie62/2006[M]. Dusseldorf: Verlag Bau Technik GmbH, 2000.

[13] 喬齡山.水泥堆積密度理論計算方法介紹[J].水泥,2007(7):1-7.

[14] 李圭白,張杰.水質(zhì)工程學(xué)[M].北京:中國建筑工業(yè)出版社,2005:148-149.

[15] Kehat E, Lin A, Kaplan A. Clogging of filter media[J]. Industrial & Engineering Chemistry Process Design and Development, 1967,6(1):48-55.

Analysis and Improvement of Mixed Particle Swarm Percolation Process in Biofiltration System

WANG Yajun1,2

(1.SchoolofCivilEngineering,LanzhouUniversityofTechnology,Lanzhou,Gansu730050,China; 2.SchoolofCivilEngineering,SoutheastUniversity,Nanjing,Jiangsu210096,China)

[Objective] To accurately master the variation law of the filtration flow field, and to provide a theoretical reference for the deeper understanding of the change in the flow in the biofiltration system. [Methods] Based on the particle swarm pressure expansion theory, characteristics of the infiltration were studied mainly in consideration of the fluid-solid coupling effect. [Results] A four stage hypothesis in the percolation process of the biofiltration system was put forward and verified by the experiment. [Conclusion] Either for the process of water seepage or sewage seepage, the permeability coefficients changed from falling to rising and then decreasing, which is caused by the particle swarm pressure expansion.

biofiltration; infiltration process; fluid-structure interaction; saturated flow; permeability coefficient

2016-10-30

2016-11-12

“十二五”國家科技支撐計劃重大項目“長三角快速城鎮(zhèn)化地區(qū)美麗鄉(xiāng)村建設(shè)關(guān)鍵技術(shù)綜合示范”(2015BAL02B05)； 江蘇省重大科技專項課題(BE2015356-2); 江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項目(CE02-1-1); 甘肅省高等學(xué)?？蒲许椖?2016B-026)

王亞軍(1979—),男(漢族),山西省介休市人，副教授，主要從事污水資源化利用研究。E-mail:wyj79626@163.com。

B

1000-288X(2017)03-0329-05

X502