摘要：為了優化工業級過濾器布流系統結構參數及提高其過濾性能，運用多孔介質模型對纖維球過濾器的過濾性能進行了分析，研究了布流系統結構參數對過濾性能的影響，并進行了實驗驗證，得出：對于直徑為1.5 m，軸徑比約為1的過濾器，優化后吐水口直徑為15～20 mm，吐水口個數為10～12，過濾性能提高了約25%；當濾料層壓降達到15 000 Pa時，濾速和出水水質不滿足生產要求，需要對濾料層進行反沖洗。濾床截污量分布越均勻，濾床徑向過濾阻力差異越小，濾床“密封性”越好，過濾器運行越穩定。

關鍵詞：纖維球濾床；工業級過濾器；布流系統優化；過濾性能

中圖分類號： TU991.24文獻標識碼： ADOI：10.3969/j.issn.1007-791X.2024.03.0040引言

21世紀，在經濟高速發展和人們生活水平提高的同時，工業廢水、廢氣及城市生活污水的問題也日趨嚴峻，由此帶來的大量污水處理成為亟需解決的難題。纖維過濾器具有過濾速度快、過濾精度高、對水質寬容度高等優勢，在過濾行業占有重要的地位[1-3]。由于纖維過濾技術不僅能夠保證過濾效率，而且具有成本低、工藝方便等優點，國內外學者對纖維球過濾器做了大量研究工作。Hellen等[4]提出一種新的隨機模型產生方法，可以通過給定纖維的方向分布及彎曲度來產生緊密的非重疊纖維過濾介質。Wang [5]、Lux[6-7]等通過對纖維濾料的顯微圖像進行分析，研究了顆粒直徑對捕集效率，填充密度對氣相場壓力損失的影響，得到了孔隙率、纖維和氣孔的尺寸分布的各向異性等微觀結構信息。杜小磊[8]發現滌綸改性纖維球機械過濾器過率效果正比于濾床高度和纖維球填充密度。劉輝、權潔、孟慶利等[9-11]研究了入口傾角、濾層厚度、孔隙率、濾速、入口速度等參數對過濾器的過濾性能的影響，得出壓力損失與濾速和厚度近似呈線性關系，并針對各自研究的過濾器的結構參數和工藝參數進行了優化。

目前，學者對過濾器的研究以新型纖維過濾材料[12-13]和軸徑比較大的過濾器結構優化[8，9-11]為主，對軸徑比相當的過濾器的過濾性能及反洗耗能的研究較少。而軸徑比相當的過濾器，濾床的截污量沿徑向分布差異與污水布流系統息息相關。本文將針對軸徑比相當的過濾器，采用CFD技術和實驗相結合的方法對影響過濾的因素進行研究，為優化過濾器結構、提高過濾器過濾能力提供理論依據。

1過濾器過濾過程仿真模型

1.1仿真模型

本文以某公司的過濾器設備為模型基礎，忽略對過濾過程影響不大的結構，且由于過濾階段為軸對稱結構，為減少計算時間，取整體1/4區域建立模型，圖1是過濾器的三維模型。

1.2基本假設

基于過濾器過濾過程的復雜性，為了簡化計算，本文針對其數學模型進行如下假設：

1） 過濾器內流體流動為不可壓縮黏性流動；

2） 將懸浮物處理為廣義流體；

3） 忽略油污對流體流動的影響；

4） 過濾器內表面與流體之間無滑移。

1.3計算方法及邊界條件

k-ε雙方程模型在工程實際與科學研究中應用廣泛，其模擬結果與實驗吻合得很好，本文選擇標準k-ε模型來對過濾器內流體流動進行模擬。混合模型Euler多相流模型在工程數值模擬中應用非常廣泛，本文選用簡化的Mixture模型，將懸浮物看作廣義流體。過濾器濾床介質為纖維球，研究中多采用多孔介質模型來模擬纖維球[10]。另外，過濾器過濾過程流體的溫度變化不大，對流場的影響可忽略不計。

纖維球過濾器的出口為壓力出口，壓力等于標準大氣壓。入口為速度入口，根據入口的橫截面積和處理水量來計算入口的速度。其余壁面均為無滑移壁面。

2數值模擬結果分析

2.1吐水口直徑對過濾性能的影響

本節研究吐水口直徑對過濾性能的影響，以某公司的現場參數：濾速300 m3/h、濾床孔隙率93%、濾床厚度0.6 m、吐水口個數8個為初始參數。圖2～3為不同吐水口直徑對過濾器截污性能、流場和壓降的影響。

由圖2可知，吐水口直徑為15 mm和20 mm時，污水從吐水口流出后，下漩渦區域遠大于上漩渦，湍流強度小于上漩渦，這樣在濾床上表面附近的區域，污水流動較為平緩，有助于污水中的懸浮物在自重和水壓下均勻地沉降到濾床上。當吐水口直徑為25 mm和30 mm時，由于吐水口直徑較大，污水射出速度較小，流出后“傾瀉”到濾床中部，形成左右兩個漩渦，湍流強度中間高兩邊低，懸浮物在兩側堆積，濾床截污分布的均勻性較差。

由圖3可知，濾床截污量分布的越均勻，過濾壓降越大。這是因為，截污量分布不均勻時，污水在濾床截污量少的區域流動較快，過濾阻力小，因此壓力損失就小。濾床截污量分布越均勻，濾床沿徑向過濾阻力相差越小，壓力損失越均勻，越有利于濾床的穩定。

圖4、5和6分別為不同吐水口直徑下的濾床上表面湍動能徑向分布曲線、濾床中截面截污量徑向分布曲線和濾床壓力損失曲線（其他工藝參數恒定：濾速300 m3/h、濾床孔隙率93%、濾床厚度0.6 m、吐水口個數8個）。

由圖4和5可知，濾床上表面沿徑向湍動能的大小直接影響著徑向截污量分布，在湍動能大的位置，截污量少，湍動能小的位置，截污量多。

由圖5和圖6可知，濾床在截污量多的區域壓降大，截污量少的區域壓降小，濾床沿徑向截污量分布越均勻，壓降大小沿徑向分布就越均勻。這是因為截污量多的區域過流阻力大，因此壓降梯度大，同樣的厚度下壓降就大。當吐水口直徑為15 mm和20 mm時，濾床徑向方向不同位置截污量和壓降分布相對均勻，有利于濾床高度和過濾過程的的穩定。因此，在吐水口個數一定的前提下，此過濾器的污水布流系統的吐水口直徑宜在15～20 mm之間。

2.2吐水口個數對過濾性能的影響

本節研究吐水口個數對過濾性能的影響，吐水口直徑為15 mm，其他工藝參數同2.1節。圖7～9為不同吐水口個數下濾床上表面湍動能徑向分布曲線、濾床中截面截污量徑向分布曲線和吐水口個數對濾床壓力損失的影響曲線。

由圖7和8可知，隨著吐水口個數的增加，濾床上表面沿徑向湍動能的大小差異越來越小，濾床徑向截污量分布也越來越均勻，吐水口個數為10～12個時，濾床徑向截污量分布差異較小。

由圖8和圖9可知，濾床不同位置的截污量隨著吐水口個數的增加而增加，濾床壓力損失也越來越大，但濾床徑向截污量和壓力損失差異變小。這是因為，吐水口個數越多，污水布流越均勻，濾床徑向截污量分布越均勻，濾床厚度方向的壓力梯度損失沿徑向分布越均勻，濾床不容易“泄壓”，“保壓性”越好濾床的壓力損失就越大。

3過濾實驗及結果分析

本文中的過濾器以秦皇島某公司生產的過濾器為原型，實驗裝置的示意圖如圖10所示。本試驗測定的主要參數有：進、出水各取樣口濁度，濾床各測壓口壓力，出水口濾速等。濁度由在線濁度分析儀測定，壓力采用測壓管測量，濾速采用一體型電磁流量計測量。

3.1壓降實驗

圖11是過濾器壁面附近濾料層縱向不同位置不同時刻的靜壓實驗數據和仿真數據。

由圖11可知，實驗數據和仿真數據的變化規律基本一致，證明本文建立的仿真模型是準確的。濾料層厚度方向的靜壓與距濾料層底部距離的大小近似成拋物線增長關系，隨著過濾時間的延長，濾料層厚度方向各位置靜壓逐漸變大，且中上層靜壓增長明顯高于中下層，不同位置的靜壓壓差逐漸變大。這是因為，在過濾剛開始的階段，由于濾料層的孔隙率較大，且濾料較為干凈，濾料層上附著的懸浮物累積量較少，濾料層的孔隙率變化不大，因此在剛開始的階段，濾料層不同位置的壓差較小，隨著濾料層上附著懸浮物量的增長，且由于懸浮物首先附著在濾料的中上層，再逐漸向中下層滲透，中上層的孔隙率隨著過濾時長的增加急劇變小，這樣中上層的過濾阻力明顯大于中下層，因此中上層的靜壓增長越來越快，導致濾料上下層的壓差越來越大。

為了驗證優化后的過濾器吐水口系統對過濾性能的提升效果，參考仿真分析結果：吐水口直徑在15～20 mm，吐水口個數在10～12個，選取吐水口直徑為20 mm、吐水口個數為10作為優化后的過濾器結構參數，進行過濾性能對比實驗。優化前后的吐水口系統結構參數如表1所示。

圖12是過濾時長6 h，過濾器吐水口結構參數優化前后濾床的壓力損失的實驗數據。

由圖12可知，采用優化后的吐水口進行過濾實驗，整個濾床的壓力損失沿濾料徑向方向更為均勻，優化前壓力損失方差為328 055.6，優化后為77 247.2，這樣有利于濾床徑向方向截留積泥的均勻分布，以便濾床更加穩定地工作。

3.2過濾性能實驗

圖13和圖14分別是采用優化后的吐水口布流系統進行過濾實驗，得出的壓降-時間關系曲線圖和濾速-時間關系曲線圖，圖15是優化前和優化后出口濁度值-時間關系曲線圖。

由圖13、圖14和圖15可知，采用優化后的布流系統進行過濾實驗，當濾料層的壓降損失超過8 000 Pa后，隨著過濾時間的增加，濾料層壓降損失快速增加，當壓降達到15 000 Pa附近時，濾料阻塞嚴重，濾速急劇減小，且出水質量嚴重下降，達不到工業排放標準。此時，需要停止過濾階段，對濾料層進行反沖洗過程。另外，由圖15可知，優化后過濾時長為5個多小時，優化前過濾時長大約為4小時，通過優化布流系統提高了過濾時長，過濾性能提升了約25%。

4結論

本文以直徑為1.5 m，軸徑比約為1的過濾器為研究對象，采用多孔介質和多相流模型，研究了污水布流系統和填充密度等因素對濾床過濾性能的影響，對上述系統結構參數進行了優化，并開展了過濾實驗，得出：

1） 濾床的靜壓大小與濾床截污量成正相關，截污量多的區域，靜壓和靜壓梯度變化大；濾床截污分布越均勻，濾床徑向過濾阻力差異越小，濾床“密封性”越好，越不易“泄壓”，有利于過濾的穩定運行；濾床在截污量多的區域壓降大，截污量少的區域壓降小。

2） 對纖維球過濾器的布流系統結構參數和濾床填充密度進行了優化，得出：過濾器布流系統吐水口直徑在15～20 mm之間，吐水口個數在10～12個，優化后的過濾器過濾性能提升了約25%。

3） 濾料層厚度方向的靜壓與距濾料層底部距離的大小近似成拋物線增長關系；采用優化后的吐水口，濾床壓力損失沿徑向方向更均勻，有利于濾床徑向截留積泥的均勻分布，有利于濾床的穩定；當濾料層的壓降損失超過8 000 Pa后，濾料層壓降損失快速增加，當壓降達到15 000 Pa附近時，濾速急劇減小，且出水質量嚴重下降，需停止過濾，對濾料進行反沖洗。

參考文獻

[1] 徐芳芳. 纖維材料過濾特性的數值模擬[D]. 上海： 東華大學， 2012： 3-8.

XU F F. Numerical simulation on the filtration characteristic of fiber filter[D]. Shanghai： Donghua University，2012： 3-8.

[2] 趙軍. 纖維球過濾器的結構設計與數值模擬[D]. 北京： 北京化工大學， 2012： 1-3.

ZHAO J. Structural design and numerical simulation of fiber ball filter[D]. Beijing： Beijing University of Chemical Technology， 2012： 1-3.

[3] 楊林. 纖維球過濾器反洗時間過長故障處理[J]. 通用機械， 2018（4）：60-66.

YANG L. Long backwash time of fiber ball filter troubleshooting[J]. General Purpose Machinery， 2018（4）： 60-66.

[4] ALTENDORF H， JEULIN D. Random-walk-based stochastic modeling of three-dimensional fiber systems[J]. Physical Review E Statistical Nonlinear amp; Soft Matter Physics， 2011， 83（4）：22-28.

[5] WANG Q， MAZEA B， TAFRESHI H V， et al. A case study of simulating submicron aerosol filtration via lightweight spun-bonded filter media[J]. Chemical Engineering Science， 2006， 61（15）： 4871-4883.

[6] LUX J， DELISE C， THIBAULT X. 3D characterization of wood based fibrous materials： an application[J]. Image Analysis amp; Stereology， 2006， 25（1）： 25-35.

[7] DELISEE C， LUX J， MALVESTIO J. 3D morphology and permeability of highly porous cellulosic fibrous material [J]. Transport in Porous Media， 2010， 8（3）： 623-636.

[8] 杜小磊. 纖維球機械過濾器處理含油廢水過程優化研究[D]. 天津： 天津工業大學， 2018： 43-58.

DU X L. Fiber ball mechanical filter oil-bearing wastewater treatment process optimization research[D]. Tianjin： Tianjin Polytechnic University， 2018： 43-58.

[9] 劉輝， 龍儼麗， 劉萍， 等. 基于CFD的排污管道過濾器數值模擬及結構優化[J]. 廣州化工， 2021， 49（1）： 97-103.

LIU H， LONG Y L， LIU P， et al. Numerical simulation and structure optimization of sewage pipeline filter based on CFD[J]. Guangzhou Chemical Industry， 2021， 49（1）： 97-103.

[10] 權潔， 陳浩， 孟慶睿. 基于Fluent的纖維過濾器內部流場數值模擬[J]. 液壓與氣動， 2019（12）： 5.

QUAN J， CHEN H， MENG Q R. Numerical simulation of flow field inside fiber filter based on Fluent[J]. Hydraulic and pneumatic， 2019（12）： 5.

[11] 孟慶利. 基于Fluent和Workbench彗星式纖維模型流場數值模擬與分析[J]. 煤炭技術， 2020， 39（12）： 152-155.

MENG Q L. Flow field numerical simulation and analysis of comet fiber based on Fluent and Workbench[J]. Coal Technology， 2020， 39（12）： 152-155.

[12] 諸文旎. 隨機纖維過濾介質模型建立及其過濾特性數值模擬研究[D]. 杭州： 浙江理工大學， 2022： 25-40.

CHU W N. Geometric modeling of random fiber filter and numerical simulation of its filtration performance[D]. Hangzhou： Zhejiang Sci-Tech University， 2022： 25-40.

[13] 付澤明. 纖維過濾材料水清洗再生特性的數值研究[D]. 西安： 西安建筑科技大學， 2022： 36-51.

FU Z M. Numerical study on regeneration characteristics of fibrous filters by water cleaning[D]. Xi’an： Xi’an University of Architecture and Technology， 2022： 36-51.

Study on distribution system optimization and

filtration performance of fiber ball filter

Abstract： In order to optimize the structural parameters of the industrial filter flow system and improve its filtration performance， the porous media model was used to analyze the filtration performance of the fiber ball filter， and the influence of the structural parameters of the flow system on the filtration behavior was analyzed， and experimental verification was carried out. The results are as follows： for the filter with diameter of 1.5 m and axial diameter ratio of about 1， the optimized outlet diameter is 15～20 mm and the number of outlet is 10～12，the filtration performance is improved about 25%. When the pressure drop of filter material reaches 15 000 Pa， the filtration rate and effluent quality do not meet the production requirements， and the filter material layer needs to be backwashed. The more uniform the distribution of filter bed pollution， the smaller the difference of filter bed radial filtration resistance， the better the \"sealing\" of filter bed， and the more stable the filter operation.

Keywords： fiber ball filter bed; industrial filter; distribution system optimization; filtration performance