基于CFD技術的臭氧給水處理工藝研究進展

陳翔宇， 姚娟娟， 王慶濤， 劉 存

(重慶大學 三峽庫區生態環境教育部重點實驗室， 重慶 400045)

現今，由于水環境的惡化，加劇了用水需求與水源污染的矛盾，使水資源問題成為制約經濟發展和人民生活水平改善的環境問題之一[1]。人工合成的有機物是水體有機污染的主要成分，在環境中不易被去除，對人體健康和水環境均構成了嚴重威脅。目前，常規的給水處理工藝難以去除水體中的痕量有機污染物，因此迫切需要對水源水進行深度處理。臭氧處理以其對水中有機物的高效去除以及良好的氣味控制而在給水處理中越發受到重視[2]。作為一種新興的水處理研究方法，計算流體動力學(Computational Fluid Dynamics，CFD)已經較多地被應用于臭氧工藝的優化與研究當中。本文對CFD技術在臭氧水處理工藝的研究與應用進行了總結。

1 臭氧給水處理工藝與CFD簡介

臭氧處理工藝在傳入我國后取得了良好的給水處理效果。臭氧工藝的主要反應器為臭氧接觸池，目前全國已建及在建具有臭氧工藝的水廠共75個，大部分集中在東南沿海地區[3]，現今國內外的常用池型為隔板式接觸池。然而它在運行時往往存在相當大體積的回流區與死水區，極大影響了運行效率。而且因池體封閉及臭氧的毒性，目前對池內的化學反應及有機物降解狀況缺少詳細的了解。因此，如何提高運行效率以及深入探索池內的化學反應機理成為目前國際上臭氧工藝研究的一個重點[4]。

CFD是利用計算機的數值計算和圖像顯示，對包含流體流動等的相關物理系統所做的分析。目前，CFD已經越來越多地被應用于給水處理的研究，許多由于條件限制而無法進行實驗的難題因而得到解決，并且節省了大量的實驗費用[5-7]。由于它與實際工程的高度貼合性及優良的經濟效益，近十幾年CFD技術也不斷地被國內外研究者應用于對臭氧接觸池的研究之中[8]。

2 CFD對臭氧接觸池內流場的表征與建模

2.1 湍流模型的選取與建立

建立精確的CFD模型是研究成功的必要前提。CFD技術發展至今，已經基本發展出了完善的湍流模型，而應用最廣泛的是k-ε模型[9]。在現有的湍流模擬中，為方便研究，包含了如下基本假設：(1)反應器內各相連續，均符合連續性方程；(2)反應器內氣液兩相體積比之和恒為1；(3)臭氧曝氣時所產生氣泡為球形且尺寸恒定，水壓對氣泡的影響不計；(4)各相在模型網格內的運動均遵守動量守恒定律等。發展較為完備的數值模擬法有Reynolds平均方法(RANS)、大渦模擬(LES)和直接數值模擬(DNS)3種。RANS法通過對Navier-Stokes(N-S)基本流動方程進行系統平均后得到湍流平均量方程，目前使用最為廣泛。k-ε模型即為RANS的子模型；LES法通過對N-S方程進行低通濾波得到描述湍流大尺度運動的方程；DNS法則是直接求解完整的三維非定常N-S方程組，它的精確度最高，但計算耗時耗力，目前很少使用[10]。

雖然目前RANS模型應用最多，但Kim等[11]用LES湍流模型來模擬池內湍流后，表明LES模型在池內呈現出了被RANS模型忽視的細微渦流，這些渦流對近壁面流態產生了影響，且LES模擬出的示蹤劑停留時間分布曲線(RTD)與中試規模得出的實驗數據更為吻合。Zhang等[12]同時運用RANS與LES模型對同一池型進行模擬對比，發現LES模型確實可以表征出RANS模型所忽略的的小渦流，但是RANS模型對池內整體渦流和死水區的表征則更為準確，這也表明RANS的整體模擬結果優于LES模型。從現有研究來看，RANS模型與LES模型各有優劣，但目前LES模型的應用仍會大量消耗計算機內存與計算時間，因此還很難在工程應用中推廣。

2.2 接觸池運行效率評價指標的建立與發展

水力效率的評估需要直觀的評價指標，目前對臭氧接觸池的水力效率評價指標較為常用的是：(1)RTD分布曲線：即累積液齡分布函數F(t)，表示流過系統的物質粒子中停留時間小于t的物質粒子所占的百分率。(2)T10/HRT:即T10是從入口投加示蹤劑后，10%示蹤劑的流出時間,HRT是示蹤劑的平均水力停留時間。如果是理想推流(活塞流)則T10/HRT=1。T10/HRT是目前最常用的評價指標；(3)CT值：即水中臭氧濃度對其作用時間的積分值。實際運行中，池中臭氧濃度分布往往不均勻，因此CT值計算還應考慮到空間的濃度變化、臭氧衰減等因素，計算非常復雜。目前普遍采用的是美國環保局(EPA)推薦的簡化計算方法即計算出水CT10值[13]。C為出水剩余臭氧濃度。

然而以上指標均存在一定不足。RTD曲線會受到模型網格密度的影響，T10/HRT則只單一考慮到了水力指標而并未考慮到池內臭氧的氧化效率；對于CT10而言，池內流態并非推流，存在回流與死水區，T10會使得實際CT值被低估[14]。且出水口處的臭氧濃度無法代表全部接觸池中的臭氧濃度，因臭氧在水中不斷衰減，因此從進水口處流到出水口處，臭氧的濃度在時時降低，因此如果只是采用出水臭氧濃度替代C值，將同樣在相當程度上低估池內實際的CT值。

因此，有研究者意識到需要建立更為全面的評價指標。Zhang等[15]對Tampa水廠的臭氧接觸池進行2D建模，比較了RTD、T10/HRT、CT10各自的優點與局限，以全面的角度，提出一種融合了環境、經濟、能耗等各種因素在內的綜合加權指標——CI(composite indicator)，用以綜合評估接觸池的運行效率及其所帶來的環境、電能消耗、經濟費用等負擔。實驗證明這一指標的運用兼顧了接觸池的高效運行與水廠的經濟效益，具有很高的推廣價值。

評價指標隨著CFD建模的不斷發展而變化。CFD建模經歷了4個階段：(1)對池內流場的還原；(2)模擬示蹤劑運輸以評價水力效率；(3)對池內臭氧濃度及衰減的表征；(4)對池內有機物降解及相關化學反應的模擬。評價指標如今日趨豐富與多元化，因此對于接觸池運行的評估已不能依靠單一指標，而應該將多種指標綜合運用。

3 應用CFD對接觸池運行參數及結構的模擬優化

3.1 對接觸池運行參數及結構的優化

隔板式臭氧接觸池在運行時往往存在大片的回流與死水區，在很大程度上影響到接觸池的運行效率，池內的氣液傳質會受到很大影響，進而影響到接觸池內臭氧的脫色除嗅及有機物的去除，因此優化接觸池的結構與外觀是十分必要的。優化結構的主要目的是使得池內的流態最大程度地接近推流，使池內氣液傳質更加充分等。在實際工程中，改動結構會耗費巨大成本，因此應用CFD進行模擬優化則具有巨大的優勢[16]。

一個普遍被采用的改進方法就是在池內增設導流板或橫擋板，目的是將流場進行分割，擴大水流的有效流動范圍，以減少池內的回流區和死水區。Wols等[17]對荷蘭典型的臭氧接觸池型進行二維建模后，在池體內隔板的上、下方設置一定角度的導流板，之后又增設橫擋板，利用示蹤劑法繪制出不同改進后的RTD分布曲線，結果發現無論設置導流板還是橫擋板都能夠極大提高池內的水力效率。但設置斜向導流板后的RTD曲線相較于橫擋板更為接近理想推流，池體水力效率相較于原始池型增加了近50%。紀家林[18]對深圳筆架山水廠的主接觸池進行了二維模擬，對在池內單獨增加導流板或橫擋板和同時增加導流、橫擋板3種方案進行了對比，以T10/HRT值作為評價指標。結果發現同時增設導流與橫擋板效果最佳，通過該方案使得該水廠的T10/HRT在原設計基礎上提高了近70%。陸麗等[19]通過不同池型內流場的模擬和HRT的計算，發現水平擋板與導流板相結合時效果卻不如單獨設置水平擋板，設置水平擋板后水力效率比兩者結合時要高出近30%。

還有一種較常采用的方法就是增加池內反應室的個數。Kim等[14]采用LES模型模擬發現適當減小隔室長度，即在池型長度不變的情況下增加反應室數目能提高池內的水力效率，隨后Zhang等[15]分別對增設擋板與增加反應室數目這兩種優化方案進行了對比，發現增加隔室數目可將池內的水力效率達到相當理想的水平，但是過多增加反應室數目會引起過高的運行能耗，提高了運營成本。因此作者創立了新的指標評測后發現應當將反應室數目控制在一定范圍內，在保證水力效率提高的同時減少不必要的能量損失。

除此之外， 不少研究者將CFD模擬與水廠實測、中試實驗相結合，對影響臭氧接觸池效率的相關因素做了深入的探究，積累了豐富的成果。表2總結了眾多研究者的相關結論。

3.2 對臭氧在池內的傳質、衰減及化學反應的模擬

臭氧在池內如何衰減與傳質的問題歷來是研究的重點，因為它直接關系到接觸池的運行效率以及相關消毒副產物如溴酸鹽的生成情況。有些研究者[17,29]將CFD技術與中試實驗相結合，在一定簡化的基礎上，得出池內的臭氧衰減與溴酸鹽生成均符合一級反應動力學模型，而池內的有機物降解過程符合Chick-Watson二級反應動力學。Zhang等[22]對美國Tampa水廠臭氧接觸池進行三維建模，并通過開源數據編程向池內引入臭氧的衰減、傳質及溴酸鹽濃度變化，模擬得出的結果與水廠的定點取樣值基本吻合，驗證了臭氧與溴酸鹽反應動力學模型的準確性。雖然有少數研究認為臭氧的自衰減更符合二級反應動力學。不過，仍然有很多研究者[4,30-31]證明以一級反應動力學來預測臭氧的衰減效果更理想，前提是被氧化的有機物不會過快地自我分解。

表2 臭氧接觸池相關影響因素模擬研究總結

雖然Chick-Watson二級反應動力學以其高度簡化和普適性被許多研究者所采用，但是它沒有考慮有機物初始濃度對降解的影響，所以計算結果往往與實際有所偏差。實際工程中的反應級數很可能介于一級與二級之間。戚圣琦[32]在前人CFD研究的基礎上自制中試實驗裝置模擬反應池內有機物的降解情況，發現反應等級與水源水種類密切相關，當以地下水作為水源時，有機物降解符合Chick-Watson二級反應動力學，但當以腐殖酸配水作為水源時，反應動力學則介于一二級之間。

表3列出了目前臭氧接觸池CFD模擬研究所常用的不同物質對應的反應動力學方程。從表3中可以看出，各物質相關的反應常數對模擬預測的準確性有著重要的影響。而反應常數取決于特定的環境參數，如pH、溫度等。因此在模擬時，需要基于所研究接觸池的具體環境，仔細選用合適的反應常數。

表3 臭氧接觸池內不同物質的反應動力學模型

4 結 論

隨著數值模擬的日益進步及計算機運算能力的增強，CFD在臭氧接觸池研究中的應用已經越發廣泛。然而，CFD的應用仍然存在如下問題：

(1)目前研究偏向于模擬水流，即單一的液相，而忽視了氣液兩相的相互作用，未來應當更多地模擬氣液兩相流，將氣泡尺寸、氣液傳質等相互影響因素加入CFD模型中。

(2)池內的實際臭氧氧化過程往往由一系列元素之間的化學反應組成，反應系統復雜，因此有必要探究出較為完整的反應機制。

(3)池內有機物的降解機理至今尚未完全明確，降解動力學常數往往會隨著有機物種類、pH、溫度等因素產生變化，具有不確定性，且降解反應究竟屬于一級、二級甚至更高級反應仍有待探究。

(4)池內湍流與化學反應之間的相互作用仍缺乏深入的研究；而池內湍流很可能對化學反應有較大的影響。

(5)目前有關臭氧接觸池的CFD模擬大多數仍停留在二維模擬，而簡化后的二維模擬往往會導致模擬結果的失準，因此未來應當以三維建模為主。

參考文獻:

[1] SUI Qian, WANG Bin, ZHAO Wentao, et al. Identification of priority pharmaceuticals in the water environment of China[J]. Chemosphere, 2012, 89 (3): 280-286.

[2] WANG Qun, YANG Zhichao, MA Jun, et al. Study on the mechanism of cerium oxide catalytic ozonation for controlling the formation of bromate in drinking water[J]. Desalination and Water Treatment, 2015, 57 (33): 15533-15546.

[3] 洪覺民. 給水深度處理技術原理與工程案例[M]. 北京:中國建筑工業出版社, 2013.

[4] ZHANG Jie, TEJADA-MARTINEZ A E, ZHANG Qiong. Developments in computational fluid dynamics-based modeling for disinfection technologies over the last two decades: A review[J]. Environmental Modelling & Software, 2014, 58(4): 71-85.

[5] YAN Hexiang, LIPEME K G, GONZALEZMERCHAN C, et al. Computational fluid dynamics modelling of flow and particulate contaminants sedimentation in an urban stormwater detention and settling basin.[J]. Environmental Science & Pollution Research International, 2014, 21(8):5347-5356.

[6] ALALM M G, NASR M, OOKAWARA S. Assessment of a novel spiral hydraulic flocculation /sedimentation system by CFD simulation, fuzzy inference system, and response surface methodology[J]. Separation and Purification Technology, 2016, 169: 137-150.

[7] GHADIRI M, FAKHRI S, SHIRAZIAN S. Modeling and CFD simulation of water desalination using nanoporous membrane contactors[J]. Industrial & Engineering Chemistry Research, 2013, 52 (9): 3490-3498.

[8] CHOI J W, KIM S S, KIM J H, et al. Study on performance improved design of pressure-type ozone contactor in multistorey water treatment plant by CFD[J]. The Korean Society of Water and Wastewater, 2014, 28 (4): 431-440.

[9] SCHMALFUβ S, SOMMERFELD M. Numerical and experimental analysis of fluid phase resonance mixers[J]. Chemical Engineering Science, 2017, 173: 570-577.

[10] DEMIREL E, ARAL M. Unified analysis of multi-chamber contact tanks and mixing efficiency based on vorticity field. Part I: Hydrodynamic Analysis[J]. Water, 2016, 8 (11): 495.

[11] KIM D J, KIM D, KIM J H, et al. Large eddy simulation of flow and tracer transport in multichamber ozone contactors[J]. Journal of Environmental Engineering, 2010, 136(1):22-31.

[12] ZHANG jie, TEJADA-MARTINEZ A E, ZHANG Qiong. Reynolds-averaged navier-stokes simulation of the flow and tracer transport in a multichambered ozone contactor[J]. Journal of Environmental Engineering-Asce, 2013, 139 (3): 450-454.

[13] ANGELOUDIS A, STOESSER T, FALCONER R A. Predicting the disinfection efficiency range in chlorine contact tanks through a CFD-based approach[J]. Water Research, 2014, 60(1): 118-129.

[14] KIM D, STOESSER T, KIM J H. Modeling aspects of flow and solute transport simulations in water disinfection tanks[J]. Applied Mathematical Modelling, 2013, 37 (16-17): 8039-8050.

[15] ZHANG Jie, TEJADA-MARTINEZ A E, LEI Hongxia, et al. Indicators for technological, environmental and economic sustainability of ozone contactors[J]. Water Research, 2016, 101: 606-616.

[16] WOLS B A. CFD in drinking water treatment[M]. Delft University of Technology, 2010.

[17] WOLS B A, Hofman J A M H, Uijttewaal W S J, et al. Evaluation of different disinfection calculation methods using CFD[J]. Environmental Modelling & Software, 2010, 25 (4): 573-582.

[18] 紀家林. 基于CFD的臭氧接觸池優化[D]. 哈爾濱：哈爾濱工業大學, 2010.

[19] 陸 麗, 鄧保慶, 李勛棟, 等. 不同臭氧接觸池性能的CFD模擬[J]. 水資源與水工程學報, 2012, 23 (3): 81-83.

[20] TALVY S, DEBASTE F, MARTINELLI L, et al. Development of a tool, using CFD, for the assessment of the disinfection process by ozonation in industrial scale drinking water treatment plants[J]. Chemical Engineering Science, 2011, 66 (14): 3185-3194.

[21] ZHANG Jianping, HUCK P M, STUBLEY G D, et al. Application of a multiphase CFD modelling approach to improve ozone residual monitoring and tracer testing strategies for full-scale drinking water ozone disinfection processes[J]. Journal of Water Supply Research and Technology-Aqua, 2008, 57 (2): 79-92.

[22] ZHANG Jie, TEJADA-MARTINEZ A E, ZHANG Qiong, et al. Evaluating hydraulic and disinfection efficiencies of a full-scale ozone contactor using a RANS-based modeling framework[J]. Water Research, 2014, 52(4): 155-167.

[23] YANG Jingxin, LI Ji, ZHU Jia, et al. A novel design for an ozone contact reactor and its performance on hydrodynamics, disinfection, bromate formation and oxidation[J]. Chemical Engineering Journal, 2017, 328: 207-214.

[24] GREENE D J, HAAS C N, FAROUK B. Computational fluid dynamics analysis of the effects of reactor configuration on disinfection efficiency[J]. Water Environment Research, 2006, 78 (9): 909-919.

[25] LI J, ZHANG J, MIAO J, et al. Application of computational fluid dynamics (CFD) to ozone contactor optimization[J]. Water Science & Technology, 2006, 6 (4): 505-512.

[26] KIM D, STOESSER T, KIM J H. The effect of baffle spacing on hydrodynamics and solute transport in serpentine contact tanks[J]. Journal of Hydraulic Research, 2013, 51 (5): 558-568.

[27] PEPLINSKI D K, DUCOSTE J J. Lessons for applying computational fluid dynamics modeling to disinfection clearwells[C]// World Water and Environmental Resources Congress，2001:1-10.

[28] DEMIREL E, ARAL M. Unified analysis of multi-chamber contact tanks and mixing efficiency evaluation based on vorticity field. Part II: Transport analysis[J]. Water, 2016, 8 (11): 495.

[29] VAN DER HELM A W C, SMEETS P W M H, BAARS E T, et al. Modeling of ozonation for dissolved ozone dosing[J]. Ozone Science & Engineering, 2007, 29 (5): 379-389.

[30] BARTRAND T A, FAROUK B, HAAS C N. Countercurrent gas/liquid flow and mixing: Implications for water disinfection [J]. International Journal of Multiphase Flow, 2009, 35(2):171-184.

[31] KIM D, CHO J. Obtaining real-time kinetic parameters on ozone decay from a full-scale ozone contactor using the axial dispersion reactor model[J]. Ozone Science & Engineering, 2014, 36 (1): 100-109.

[32] 戚圣琦. 臭氧接觸氧化模型及臭氧接觸池優化研究[D]. 北京：清華大學, 2015.

[33] HUANG T H, BROUCKAERT C J, PRYOR M, et al. Application of computational fluid dynamics modelling to an ozone contactor[J]. Water S A, 2004, 30 (1): 51-56.

[34] TANG G， ADU-SARKODIE K， KIM D， et al. Modeling cryptosporidium parvum oocyst inactivation and bromate formation in a full-scale ozone contactor[J]. Environmental Science and Technology,2005,39(23):9343-9350.

[35] ZHANG Jianping. An integrated design approach for improving drinking water ozone disinfection treatment based on computational fluid dynamics[D]. Canada:University of Waterloo, 2006.

[36] ZHANG Jian, HUCK P M, ANDERSON W B, et al. A computational fluid dynamics based integrated disinfection design approach for improvement of full-scale ozone contactor performance[J]. Ozone Science & Engineering, 2007, 29 (6): 451-460.