李 俊，王 蕾，張曉璐，劉清海，郭 堅(jiān)，黃冬梅

(1.中國(guó)電建集團(tuán)華東勘測(cè)設(shè)計(jì)研究院有限公司，浙江杭州310014；2.南京慧邦科技研究所，江蘇南京210037；3.中國(guó)計(jì)量大學(xué)質(zhì)量與安全工程學(xué)院，浙江杭州310018)

砂石料加工廢水沉淀池中顆粒絮凝沉降過(guò)程模擬研究

李 俊1，王 蕾1，張曉璐1，劉清海2，郭 堅(jiān)1，黃冬梅3

(1.中國(guó)電建集團(tuán)華東勘測(cè)設(shè)計(jì)研究院有限公司，浙江杭州310014；2.南京慧邦科技研究所，江蘇南京210037；3.中國(guó)計(jì)量大學(xué)質(zhì)量與安全工程學(xué)院，浙江杭州310018)

首先測(cè)定某水電站砂石料加工廢水處理系統(tǒng)中沉淀池內(nèi)不同深度處顆粒的粒徑分布，在此基礎(chǔ)上利用三維無(wú)網(wǎng)格擴(kuò)散限制集團(tuán)凝聚模型(DLCA)對(duì)沉淀池內(nèi)顆粒絮凝沉降過(guò)程進(jìn)行模擬計(jì)算，研究絮凝沉降過(guò)程中顆粒分形結(jié)構(gòu)的演化過(guò)程。研究結(jié)果表明：沉淀池內(nèi)部顆粒粒徑分布較為均勻，顆粒粒徑主要集中在30μm左右，中位粒徑約為11μm；顆粒沉降時(shí)間、沉降顆粒層高度及顆粒層分形維數(shù)隨初始濃度增大而增大。

砂石料加工廢水；沉淀池；分形結(jié)構(gòu)； DLCA；絮凝沉降

0 前 言

沉淀池是砂石料加工廢水處理過(guò)程中物理處理重要的環(huán)節(jié)之一，其效率直接影響出水質(zhì)量及廢水處理系統(tǒng)處理效率[1]。沉淀池中細(xì)小顆粒在水中時(shí)刻不停的做布朗運(yùn)動(dòng)，并且在重力作用下克服水的浮力下沉，細(xì)小顆粒首先凝聚在一起形成較大的絮團(tuán)，絮團(tuán)越大向下的重力越大，下降速度越快，最終沉降[2-4]。在實(shí)際砂石料加工廢水處理過(guò)程中，顆粒絮凝沉降時(shí)將形成復(fù)雜的絮凝體，然而卻很難對(duì)其準(zhǔn)確觀察。隨著計(jì)算機(jī)技術(shù)的發(fā)展，采用計(jì)算機(jī)對(duì)顆粒絮凝沉降過(guò)程進(jìn)行仿真，觀察顆粒絮凝沉降過(guò)程中絮凝體的生長(zhǎng)過(guò)程已成為研究顆粒沉降的主要手段之一。

砂石料加工廢水沉淀池中細(xì)小顆粒的絮凝沉降過(guò)程是一種典型的分形生長(zhǎng)現(xiàn)象[5-6]。目前采用計(jì)算機(jī)模擬對(duì)其進(jìn)行研究主要是基于擴(kuò)散限制凝聚(DLA)或擴(kuò)散限制集團(tuán)凝聚模型(DLCA)。DLA方法假設(shè)顆粒存在一個(gè)凝聚中心，隨機(jī)釋放的粒子向凝聚中心運(yùn)動(dòng)，不能模擬顆粒的隨機(jī)運(yùn)動(dòng)及沉降，其模擬結(jié)果與實(shí)際情況相距甚遠(yuǎn)。DLCA假設(shè)初始條件下，顆粒隨機(jī)分布，不設(shè)置固定中心，允許模擬空間內(nèi)顆粒和絮凝體隨機(jī)運(yùn)動(dòng)，形成顆粒更多、尺度更大、內(nèi)部結(jié)構(gòu)更復(fù)雜的絮團(tuán)。因此更接近于實(shí)際水中的絮凝沉降過(guò)程[7-11]。然而由于傳統(tǒng)三維DLCA計(jì)算需求巨大，因此目前主要集中在采用二維DLCA進(jìn)行模擬計(jì)算。

本研究以某水電站砂石料沖洗廢水處理過(guò)程中沉淀池為研究對(duì)象，首先測(cè)定不同高度處沉淀池內(nèi)顆粒粒徑分布，基于試驗(yàn)測(cè)定結(jié)果，采用三維無(wú)網(wǎng)格DLCA模型，對(duì)沉淀池內(nèi)顆粒絮凝沉降過(guò)程進(jìn)行模擬計(jì)算，觀察其分形生長(zhǎng)過(guò)程。顆粒絮凝沉降過(guò)程一方面進(jìn)行隨機(jī)布朗運(yùn)動(dòng)，另一方面在重力作用下下沉，當(dāng)計(jì)算區(qū)域中所有顆粒均沉降到底部，計(jì)算結(jié)束。

1 絮凝沉降過(guò)程計(jì)算機(jī)模擬

采用三維無(wú)網(wǎng)格DLCA對(duì)絮凝沉降過(guò)程進(jìn)行模擬計(jì)算。模擬過(guò)程中粒子進(jìn)行無(wú)規(guī)則布朗運(yùn)動(dòng)，顆粒布朗運(yùn)動(dòng)位移采用愛(ài)因斯坦布朗運(yùn)動(dòng)位移公式計(jì)算，如下[11]：

(1)

顆粒在進(jìn)行布朗運(yùn)動(dòng)的同時(shí)，由于自身重力的作用，將克服水的浮力而不斷沉降。假設(shè)在沉淀池泥沙沉降過(guò)程中水處于靜止?fàn)顟B(tài)，則本研究采用層流區(qū)重力沉降公式計(jì)算水中顆粒沉降位移[11]：

(2)

式中，d為絮凝體直徑，m；ρs為顆粒密度，kg/m3；ρ為水的密度，kg/m3；g為重力加速度，9.81 m/s2；μ為水在溫度T時(shí)的運(yùn)動(dòng)粘度，m/s2。

計(jì)算開(kāi)始時(shí)，根據(jù)水中顆粒濃度，產(chǎn)生N個(gè)互不重疊，直徑為d的顆粒。任取某一絮凝體或顆粒進(jìn)行隨機(jī)布朗運(yùn)動(dòng)及重力沉降，當(dāng)絮凝體或顆粒與其他絮凝體或顆粒發(fā)生碰撞后，與碰撞顆粒形成新的絮凝體，不再分離，直到計(jì)算區(qū)域中所有顆粒均沉降到底部，計(jì)算結(jié)束。模擬計(jì)算過(guò)程中，顆粒直徑根據(jù)測(cè)定結(jié)果確定，計(jì)算時(shí)間步長(zhǎng)為1 s，計(jì)算區(qū)域?yàn)?00×200×200，計(jì)算區(qū)域內(nèi)粒子數(shù)采用如下公式計(jì)算[12]：

(3)

式中，c為粒子濃度；L為計(jì)算區(qū)域尺度。

對(duì)于顆粒群，有如下關(guān)系[13]：

(4)

兩邊取對(duì)數(shù)，得：

(5)

2 結(jié)果與討論

2.1 沉淀池中顆粒粒徑分布

本項(xiàng)目中設(shè)計(jì)的沉淀池容積為150 m3，廢水在沉淀池中水力停留時(shí)間約1.4 h，沉淀上清液進(jìn)入調(diào)節(jié)池。沉淀池配置有1臺(tái)鏈?zhǔn)焦文鄼C(jī)，間歇式運(yùn)行，定期將沉淀池中的污泥刮出至振動(dòng)脫水篩進(jìn)行泥水分離。沉淀池中料漿固體含量的體積分?jǐn)?shù)為15%～20%，固體量為14 t/h，料漿量為70～95 m3/h。采用激光粒徑分析儀測(cè)定0.3 m和0.7 m高度處顆粒粒徑分布。測(cè)試結(jié)果如圖1所示。由圖1可知，沉淀池兩處顆粒粒徑分布差別甚微，顆粒中位粒徑分別為10.14、11.91 μm。顆粒粒徑主要分布在0.2～100 μm之間，最大粒徑和最小粒徑分別約為95、0.2 μm。顆粒粒徑較為均勻，顆粒粒徑分布近似成偏正態(tài)分布，體積分?jǐn)?shù)為80%的粉體粒徑小于30 μm，占總體積最多的顆粒粒徑為17～29 μm。

圖1 沉淀池不同高度處粒徑分布

2.2 模擬結(jié)果分析

模擬計(jì)算中顆粒直徑取11 μm，初始狀態(tài)下顆粒在計(jì)算區(qū)域中隨機(jī)分布，計(jì)算區(qū)域內(nèi)介質(zhì)為水，初始溫度為293 K，計(jì)算體積分?jǐn)?shù)為5%、8%、10%、12%、15%、20%、25%，初始粒子數(shù)為574、919、1 149、1 378、1 723、2 297、2 871。計(jì)算結(jié)束后記錄顆粒沉降時(shí)間，沉降顆粒層厚度及分形維數(shù)。體積分?jǐn)?shù)為5%時(shí)，不同時(shí)刻顆粒絮凝體形貌如圖2所示。由圖2可知，初始時(shí)刻顆粒隨機(jī)分布于整個(gè)計(jì)算區(qū)域內(nèi)部，隨著布朗運(yùn)動(dòng)和重力沉降的進(jìn)行，顆粒碰撞，首先由多顆粒形成小絮凝體，絮凝體較之單個(gè)顆粒體積增大，質(zhì)量增大，而使其重力沉降速度加快，從而使絮凝體不斷下沉，288 s時(shí)刻顆粒主要集中在150 μm以下，588 s時(shí)刻，顆粒繼續(xù)沉降，顆粒主要集中在100 μm以下，至616 s，顆粒全部沉降至底部，泥沙高度約為80 μm。在整個(gè)沉降過(guò)程中，顆粒相互碰撞凝結(jié)，加速顆粒沉降速度。

圖2 體積分?jǐn)?shù)5%時(shí)絮凝體形貌

由對(duì)沉淀池中顆粒粒徑分布的測(cè)定可知，沉淀池中顆粒的體積分?jǐn)?shù)80%的顆粒粒徑小于30 μm，占總體積最多的顆粒粒徑為：17～29 μm。懸浮顆粒在水中做布朗運(yùn)動(dòng)和重力沉降過(guò)程中碰撞形成絮凝體，絮凝體不斷增大，做整體運(yùn)動(dòng)，由于絮凝體體積不斷增大，最終在重力的作用下，絮凝體將沉降至底部，這些絮凝體形成的結(jié)構(gòu)是典型的分形結(jié)構(gòu)，沉降結(jié)束后，絮凝體的分形維數(shù)采用式(4)計(jì)算，顆粒體積分?jǐn)?shù)15%時(shí)，分形維數(shù)計(jì)算如圖3所示。由計(jì)算可知，顆粒沉降結(jié)束后，形成的絮凝體分形維數(shù)Df為1.567 42，參數(shù)k0為-0.189 09。

圖3 體積分?jǐn)?shù)15%時(shí)分形維數(shù)

采用同樣的方法計(jì)算，得出顆粒體積分?jǐn)?shù)從5%～25%時(shí)，顆粒完全沉降時(shí)間，絮凝體高度，顆粒層分形維數(shù)。顆粒沉降時(shí)間與顆粒高度隨體積分?jǐn)?shù)變化曲線如圖4所示，由圖4可知，隨著體積分?jǐn)?shù)的增大，所需沉降時(shí)間不斷增加，沉降時(shí)間與顆粒體積分?jǐn)?shù)呈線性關(guān)系增加，由最小二乘法擬合可得本文計(jì)算條件下顆粒沉降時(shí)間(T)與體積分?jǐn)?shù)(c)的關(guān)系為T(mén)=90.74+108.17c。由圖4可知沉降結(jié)束后，顆粒高度隨體積分?jǐn)?shù)的增大而增加，當(dāng)計(jì)算區(qū)域內(nèi)顆粒體積分?jǐn)?shù)較低時(shí)，沉降結(jié)束后，出現(xiàn)明顯的顆粒區(qū)和清水區(qū)，上層較大區(qū)域內(nèi)不存在顆粒，隨著體積分?jǐn)?shù)的增大，顆粒和水的界面不斷上移，當(dāng)顆粒體積分?jǐn)?shù)大于20%時(shí)，顆粒高度大于180 μm，即當(dāng)顆粒初始體積分?jǐn)?shù)大于20%時(shí)，沉降后顆?；境錆M了整個(gè)計(jì)算區(qū)域。

圖4 顆粒沉降時(shí)間與顆粒高度隨體積分?jǐn)?shù)變化曲線

圖5為計(jì)算結(jié)束后，顆粒絮凝體分形維數(shù)隨體積分?jǐn)?shù)變化曲線。由圖5可知，顆粒絮凝體分形維數(shù)Df隨體積分?jǐn)?shù)的增大而增大。這是因?yàn)楫?dāng)計(jì)算區(qū)域內(nèi)顆粒體積分?jǐn)?shù)增大時(shí)，顆粒數(shù)增多，顆粒在布朗運(yùn)動(dòng)和沉降過(guò)程中碰撞幾率增大，使得形成的絮凝體顆粒數(shù)增多，體積增大，最終形成的絮凝體密實(shí)、結(jié)構(gòu)更加復(fù)雜。此計(jì)算結(jié)果與文獻(xiàn)[11]結(jié)果一致。

圖5 分形維數(shù)隨體積分?jǐn)?shù)變化曲線

3 結(jié) 論

本研究首先測(cè)定某砂石料加工廢水沉淀池中0.7 m和0.3 m高度處顆粒粒徑分布，基于粒徑分布的結(jié)果，采用三維無(wú)網(wǎng)格DLCA模擬沉淀池中顆粒絮凝沉降過(guò)程，顆粒在沉降過(guò)程中主要為布朗運(yùn)動(dòng)和重力沉降。模擬結(jié)果能夠形象的反應(yīng)沉淀池中顆粒懸浮體系的絮凝沉降過(guò)程，所得結(jié)論如下：

(1) 由粒徑分布測(cè)定結(jié)果可知，顆粒粒徑主要集中在30 μm左右，中位粒徑約為11 μm，沉淀池0.7 m和0.3 m高度處顆粒粒徑無(wú)明顯差異，該沉淀池中砂石料加工廢水的顆粒比較均勻；

(2) 在沉降過(guò)程中，顆粒首先形成小的絮凝體，小的絮凝體繼續(xù)運(yùn)動(dòng)，不斷碰撞其他顆粒或絮凝體，絮凝體體積不斷增大，從而使絮凝體向下沉降速度

加快，最終完全沉降到計(jì)算區(qū)域底部，所需沉降時(shí)間和沉降顆粒層高度隨顆粒初始體積分?jǐn)?shù)增大而增大；

(3) 采用DLCA計(jì)算得到的顆粒層絮凝體為分形結(jié)構(gòu)，分形維數(shù)隨顆粒初始體積分?jǐn)?shù)增大而增大，初始體積分?jǐn)?shù)越大形成的絮凝體結(jié)構(gòu)越復(fù)雜，絮凝體越致密。

[1]蔡金傍， 段祥寶， 朱亮. 沉淀池水流數(shù)值模擬[J]. 重慶建筑大學(xué)學(xué)報(bào)， 2003， 25(4)： 64-69．

[2]蔣衛(wèi)剛， 顧國(guó)維， 俞國(guó)平. 二沉池動(dòng)態(tài)仿真模型研究[J]. 環(huán)境工程學(xué)報(bào)， 2007， 1(3)： 21-24．

[3]聶德明. 顆粒沉降及其在流場(chǎng)中做布朗運(yùn)動(dòng)的研究[D]. 杭州： 浙江大學(xué)， 2011．

[4]曾光明， 葛衛(wèi)華. 污水廠二維沉淀池水流和懸浮物運(yùn)動(dòng)數(shù)值模擬[J]. 中國(guó)環(huán)境科學(xué)， 2002， 22(4)： 338-341．

[5]郭航. 水中絮凝體分形成長(zhǎng)動(dòng)態(tài)仿真技術(shù)研究[D]. 哈爾濱： 哈爾濱工業(yè)大學(xué)， 2012．

[6]李冬梅， 梅勝， 譚萬(wàn)春， 等. 黃河泥沙架橋絮凝體分形結(jié)構(gòu)的動(dòng)態(tài)演變研究[J]. 給水排水， 2005， 30(11)： 1-5．

[7]蔣新. MATLAB 平臺(tái)上的 DLCA 過(guò)程模擬[J]. 計(jì)算機(jī)仿真， 2003， 20(6)： 68-70．

[8]劉東京， 邱祖民. 反應(yīng)概率置限二維 DLCA 模型分形模擬[J]. 計(jì)算機(jī)與應(yīng)用化學(xué)， 2012， 29(4)： 451-453．

[9]齊燦， 吳世先. 應(yīng)用 DLCA 模型模擬超細(xì)氣溶膠粒子的凝并行為[J]. 桂林航天工業(yè)高等專(zhuān)科學(xué)校學(xué)報(bào)， 2008， 13(2)： 31-34．

[10]溫娟， 李宏， 王萬(wàn)軍， 等. 應(yīng)用 DLCA 模型對(duì)鋼水中夾雜物粒子集團(tuán)凝聚過(guò)程動(dòng)力學(xué)行為的模擬研究[C]∥冶金研究中心 2005 年“冶金工程科學(xué)論壇” 論文集， 2005．

[11]楊鐵笙， 李富根， 梁朝皇. 粘性細(xì)顆粒泥沙靜水絮凝沉降生長(zhǎng)的計(jì)算機(jī)模擬[J]. 泥沙研究， 2005(4)： 14-20．

[12]JULLIEN R， OLIVI-TRAIN N， HASMY A， et al. Scaling theory and numerical simulations of aerogel sintering[J]. Journal of non-crystalline solids， 1995， 188(1)： 1-10．

[13]MEAKIN P. Formation of fractal clusters and networks by irreversible diffusion-limited aggregation[J]. Physical Review Letters， 1983， 51(13)： 1119．

(責(zé)任編輯 王 琪)

Simulation of Flocculation-Sedimentation Process in Clarifier of Wastewater from Sand and Gravel Processing

LI Jun1, WANG Lei1, ZHANG Xiaolu1, LIU Qinghai2, GUO Jian1, HUANG Dongmei3

(1. PowerChina Huadong Engineering Corporation Limited, Hangzhou 310014, Zhejiang, China;2. Nanjing Huibang Institute of Technology, Nanjing 210037, Jiangsu, China;3. College of Quality and Safety Engineering, China Jiliang University, Hangzhou 310018, Zhejiang, China)

The particle size distribution at different height of clarifier is firstly measured by using laser particle size analyzer, and then a three dimensional mesh free diffusion limited cluster aggregation model (DLCA) is developed to simulate the flocculation settling process of particle in clarifier. The evolution of floc fractal is analyzed. The results indicate that, (a) the particle size uniformly distributes in clarifier and mainly concentrates in about 30 μm, and the median value is about 11 μm; and (b) with the increase of initial particle concentration, the sedimentation time, particle bed height and fractal dimension are increased．

Waste water from sand and gravel processing; clarifier; fractal structural; DLCA; flocculation-sedimentation

2016-01-06

李俊(1983—)，男，四川廣元人，工程師，碩士，主要從事水利水電工程環(huán)保設(shè)計(jì)、環(huán)境影響評(píng)價(jià)工作．

X703

A

0559-9342(2016)11-0106-04