馮金鉆,陶樂仁,黃理浩,谷志攀,2,吳生禮
(1. 上海理工大學(xué) 能源與動(dòng)力工程學(xué)院/上海市動(dòng)力工程多相流動(dòng)與傳熱重點(diǎn)實(shí)驗(yàn)室,上海 200093;2. 嘉興學(xué)院 建工學(xué)院,浙江 嘉興 314001)
污水處理處置一直是關(guān)系著人們生活和工業(yè)發(fā)展的重大問題。根據(jù)2018年12月住房和城鄉(xiāng)建設(shè)部統(tǒng)計(jì)數(shù)據(jù)顯示,2017年我國(guó)污水處理量達(dá)15743萬m·日。如果按照污水處理廠平均以負(fù)荷率80%運(yùn)行,每萬t污水產(chǎn)生6 t含水率(質(zhì)量分?jǐn)?shù))為80%的污泥來算,我國(guó)每天產(chǎn)生污泥超過7.55萬t。
污泥如果得不到妥善處理,污泥中的重金屬、細(xì)菌病原體等微生物、有機(jī)污染物將會(huì)對(duì)環(huán)境造成難以預(yù)料的損害,甚至危害到人類及其他動(dòng)植物的生存。由于污泥焚燒能徹底消滅病原體等微生物,使有機(jī)物碳化,且所釋放的能量可供熱、發(fā)電,焚燒渣可制作附屬產(chǎn)品,真正使污泥實(shí)現(xiàn)減量化、穩(wěn)定化、無害化和資源化的目的,因而被認(rèn)為是最有發(fā)展前景的污泥處理方法之一。一般經(jīng)過機(jī)械脫水的污泥含水率仍在80%以上,此時(shí)污泥的熱值接近于0,無法直接焚燒,需經(jīng)干燥預(yù)處理將污泥含水率降到50%,污泥的熱值才能利用起來。而且,干燥過程使污泥含水量降低,減小了污泥的容量,便于污泥的后續(xù)處理。
很多國(guó)內(nèi)外學(xué)者對(duì)污泥干燥展開了相關(guān)研究,大多是圍繞污泥成分分布、不同形狀污泥的干燥特性、不同干燥工藝對(duì)干燥的影響等方面去研究,缺乏對(duì)污泥干燥過程中深層次動(dòng)力學(xué)機(jī)理的研究。本實(shí)驗(yàn)研究污泥厚度和溫度對(duì)污泥干燥特性的影響,通過引入薄層干燥動(dòng)力學(xué)模型,對(duì)污泥干燥進(jìn)行動(dòng)力學(xué)模擬,得出污泥干燥的最優(yōu)模型和有效擴(kuò)散系數(shù),為污泥干燥工藝提供新的技術(shù)參數(shù)和依據(jù)。
本次實(shí)驗(yàn)污泥樣品取自嘉興市某污水廠的脫水市政污泥,呈深黑色,黏稠狀,含水率達(dá)85%。污泥的主要成分C、H、O、N、Cl和苯乙烯的質(zhì)量分?jǐn)?shù)分別為18.34%、3.82%、16.27%、1.91%、32.21%和1.69%。
實(shí)驗(yàn)中采用精度為0.01 g的德國(guó)sartorius電子天平并通過VB軟件編寫程序?qū)ξ勰噘|(zhì)量進(jìn)行30″一次數(shù)據(jù)采集。將污泥平鋪成餅狀置于溫度可控的電加熱板上。首先,控制加熱溫度為110 ℃,分別對(duì)3、6、9、12 mm厚度的污泥進(jìn)行加熱,考察污泥在不同厚度下的干燥情況;然后,控制污泥厚度為9 mm,分別設(shè)定加熱溫度為70、90、110、130 ℃,通過繪制污泥的失水曲線,探究不同溫度下污泥的干燥特性。
污泥的含水率隨時(shí)間變化,即
式中:U
為污泥的含水率,為無因次量;M
、M
、M
分別為初始時(shí)刻、干燥至t
時(shí)刻、干燥至平衡狀態(tài)下的單位干污泥含水質(zhì)量,g·g。由于M
很 小,其相對(duì)于M
和M
可以忽略,因此式(1)可簡(jiǎn)化為污泥的失水效果可用干燥速率表示,即
式中:V
為污泥的干燥速率,g·(g·min);M
為t
+Δt
時(shí)刻單位干污泥含水質(zhì)量,g·g。溫度110 ℃下不同厚度污泥的含水率隨時(shí)間的變化如圖1所示。
圖1 不同厚度的污泥含水率隨時(shí)間的變化Fig. 1 Relationship between moisture content in the sludge and time under different thicknesses
從圖1中可看出,污泥越薄,干燥至平衡含水率的時(shí)間越短。這是因?yàn)椋勰嘣奖?,水分從污泥?nèi)部遷移到污泥表面的路徑越短,傳質(zhì)阻力也越小,因此所需時(shí)間越短,水分蒸發(fā)得越快。當(dāng)污泥厚度從12 mm減小至3 mm時(shí),干燥至平衡含水率所用的時(shí)間從約200 min縮短至接近30 min。因此,將污泥鋪平使其變得更薄,會(huì)大大縮短污泥干燥時(shí)間。
不同溫度下厚度為9 mm的污泥含水率隨時(shí)間的變化如圖2所示。
從圖2中可看出,在污泥厚度相同的情況下,提高干燥溫度,污泥干燥至平衡含水率的時(shí)間越短。這是因?yàn)樘岣吒稍餃囟?,即增大了污泥?nèi)部的傳熱溫差,在傳熱溫差的驅(qū)動(dòng)作用下,污泥內(nèi)的水分遷移加快。當(dāng)溫度從70 ℃上升到130 ℃時(shí),污泥干燥至平衡含水率的時(shí)間從200 min縮短至60 min。
圖2 不同溫度下污泥含水率隨時(shí)間的變化Fig. 2 Relationship between moisture content in the sludge and time under different temperatures
圖3為不同溫度下干燥速率和時(shí)間的關(guān)系。圖4為不同溫度下干燥速率和含水率的關(guān)系。結(jié)合圖3、4可看出,當(dāng)污泥剛開始被加熱時(shí),其表面含有大量的自由水,隨著時(shí)間的增長(zhǎng),污泥表面溫度不斷上升,水分不斷蒸發(fā),干燥速率明顯上升,此為升速干燥階段。當(dāng)加熱到一定時(shí)間,污泥的自由水含量減少,間隙水和剩余的自由水變成了主要干燥對(duì)象,污泥內(nèi)部的水分向表面遷移的速率接近表面水分的蒸發(fā)速率,表現(xiàn)為近似恒速干燥,而且溫度越高,恒速干燥階段范圍越小。當(dāng)溫度達(dá)到130 ℃時(shí),恒速干燥階段已不明顯。恒速干燥階段后,污泥開始出現(xiàn)干區(qū),這增加了傳熱傳質(zhì)阻力,污泥中的水分為剩下的部分間隙水和結(jié)合水;這部分水分難以去除,因此干燥速率降低,此時(shí)為降速干燥階段。最大干燥速率出現(xiàn)在恒速干燥階段,當(dāng)溫度從70 ℃上升到130 ℃時(shí),相應(yīng)的最大干燥速率從0.02587 g·(g·min)上升到0.08158 g·(g·min)。
圖3 不同溫度下干燥速率和時(shí)間的關(guān)系Fig. 3 Relationship between drying rate and time under different temperatures
圖4 不同溫度下干燥速率和含水率的關(guān)系Fig. 4 Relationship between drying rate and moisture content under different temperatures
從圖4中還可知,在對(duì)應(yīng)的含水率區(qū)間溫度越高,恒速干燥階段范圍越小,溫度升高使得污泥更快地進(jìn)入降速干燥階段,造成更多的能量用來加熱污泥的干區(qū)。這表明,若將污泥干燥至相同的含水率,從能源利用的角度分析,溫度過高會(huì)導(dǎo)致有更多的能量消耗在污泥干質(zhì)上,不利于能源的有效利用。因此,在滿足生產(chǎn)的前提下,干燥溫度應(yīng)控制在合理的范圍。
k
、k
和k
均為模型的經(jīng)驗(yàn)系數(shù);a
、b
、c
和y
均為模型的經(jīng)驗(yàn)常數(shù)。表1 薄層干燥模型
Tab. 1 Models of thin layer drying
模型來源 表達(dá)式Lewis[9]U=exp(-kt)Page[10]U=exp(-kty)Modified Page[11]U=exp[-(kt)y]Henderson and Pabis[12]U=aexp(-kt)Logarithmic[13]U=aexp(-kt)+c Midilli[14]U=aexp(-kty)+bt Two term[15]U=aexp(-k1t)+bexp(-k2t)
以9 mm薄層污泥在70、90、110和130 ℃下的干燥曲線為例,應(yīng)用Origin軟件分別利用每個(gè)薄層干燥模型進(jìn)行數(shù)據(jù)擬合,得出模型參數(shù),并用相關(guān)系數(shù)R
、標(biāo)準(zhǔn)差S
和殘差平方和R
來評(píng)價(jià)模型的擬合優(yōu)度。R
越接近1,S
和R
越小,說明模型的擬合效果越好。表2 70℃污泥薄層干燥模型擬合結(jié)果
Tab. 2 Fitting results of thin layer drying models of the sludge at 70℃
模型來源 參數(shù)值R2SDRSS Lewis[9]k = 0.01122 0.96508 0.04429 0.03923 Page[10]k = 0.00447,y = 0.88961 0.96975 0.04122 0.03229 Modified Page[11]k = 0.01131,y = 0.88944 0.96975 0.04122 0.03229 Henderson and Pabis[12]k = 0.00992,a = 0.89912 0.98883 0.02504 0.01192 Logarithmic[13]k = 0.01007,a = 0.89418,c = 0.03767 0.98823 0.02571 0.01190 Midilli[14]k = 0.00149,y = 1.45879,a = 0.84388,b = 6.64 × 10-4 0.99988 0.00262 0.00011 Two term[15]k1 = 0.01017,k2 = -0.06407,a = 0.90575,b = 1.33 × 10-7 0.99041 0.02320 0.00915
表3 90℃污泥薄層干燥模型擬合結(jié)果
Tab. 3 Fitting results of thin layer drying models of the sludge at 90℃
模型來源 參數(shù)值R2SDRSS Lewis[9]k = 0.01122 0.96508 0.04429 0.03923 Page[10]k = 0.00447,y = 0.88961 0.96975 0.04122 0.03229 Modified Page[11]k = 0.01131,y = 0.88944 0.96975 0.04122 0.03229 Henderson and Pabis[12]k = 0.00992,a = 0.89912 0.98883 0.02504 0.01192 Logarithmic[13]k = 0.01007,a = 0.89418,c = 0.03767 0.98823 0.02571 0.01190 Midilli[14]k = 0.00149,y = 1.45879,a = 0.84388,b = 6.64 × 10-4 0.99988 0.00262 0.00011 Two term[15]k1 = 0.01017,k2 = -0.06407,a = 0.90575,b = 1.33 × 10-7 0.99041 0.02320 0.00915
表4 110℃污泥薄層干燥模型擬合結(jié)果
Tab. 4 Fitting results of thin layer drying models of the sludge at 110℃
模型來源 參數(shù)值R2SDRSS Lewis[9]k = 0.02529 0.97463 0.04065 0.03305 Page[10]k = 0.00715,y = 0.99849 0.97329 0.04170 0.03304 Modified Page[11]k = 0.02531,y = 0.99569 0.97329 0.04170 0.03304 Henderson and Pabis[12]k = 0.02303,a = 0.90947 0.98407 0.03220 0.01970 Logarithmic[13]k = 0.02448,a = 0.90025,c = 0.01737 0.98472 0.03155 0.01791 Midilli[14]k = 0.00501,y = 1.39340,a = 0.85678,b = 2.79 × 10-4 0.99946 0.00595 0.00060 Two term[15]k1 = 0.02358,k2 = -0.02739,a = 0.91548,b = 2.00 × 10-4 0.98679 0.02933 0.01463
表5 130℃污泥薄層干燥模型擬合結(jié)果
Tab. 5 Fitting results of thin layer drying models of the sludge at 130℃
模型來源 參數(shù)值R2SDRSS Lewis[9]k = 0.03755 0.96481 0.04649 0.04322 Page[10]k = 0.01625,y = 1.24201 0.97054 0.04253 0.03437 Modified Page[11]k = 0.03627,y = 1.24169 0.97054 0.04253 0.03437 Henderson and Pabis[12]k = 0.03489,a = 0.92101 0.97022 0.04277 0.03475 Logarithmic[13]k = 0.03599,a = 0.91457,c = 0.00954 0.96944 0.04332 0.03378 Midilli[14]k = 0.00326,y = 1.66047,a = 0.84789,b = 1.90 × 10-4 0.99882 0.00851 0.00123 Two term[15]k1 = 0.03528,k2 = -0.02270,a = 0.92291,b = 4.12 × 10-4 0.97058 0.04250 0.03071
表2~5為每個(gè)薄層干燥模型的擬合結(jié)果。由此可看出:Midilli模型的R
均超過0.99,在所有模型中最接近1;S
和R
在所有模型中也最小。說明Midilli模型較適用于污泥薄層干燥的模擬。圖5給出了Midilli模型擬合值和實(shí)驗(yàn)值的對(duì)比結(jié)果,其數(shù)據(jù)基本上在斜率為45°的直線上,說明Midilli模型在70~130 ℃下能夠很好地對(duì)市政污泥干燥進(jìn)行模擬。由于模型能夠描述含水率隨時(shí)間的變化,因此對(duì)于估算污泥干燥程度有重要的指導(dǎo)意義。
圖5 含水率實(shí)驗(yàn)值與Midilli模型擬合值的比較Fig. 5 Comparison between experimental values and fitting ones of moisture content by Midilli
污泥中的水分遷移是一個(gè)從內(nèi)向外擴(kuò)散的過程,根據(jù)Fick擴(kuò)散定律可以推導(dǎo)出包含有效水分?jǐn)U散系數(shù)的薄層干燥理論方程,即
式中:D
為有效擴(kuò)散系數(shù),m·s;L
為污泥厚度,實(shí)驗(yàn)中取L
= 0.009 m;n
為實(shí)驗(yàn)采樣次數(shù)。由于干燥時(shí)間較長(zhǎng),取式(4)第1項(xiàng)并兩邊取對(duì)數(shù)將其簡(jiǎn)化為
通過式(5)可以發(fā)現(xiàn),由lnU
與t
的線性關(guān)系可以計(jì)算有效擴(kuò)散系數(shù)D
。經(jīng)計(jì)算,70、90、110、130 ℃下有效擴(kuò)散系數(shù)分別為5.96 ×10、1.04 × 10、1.54 × 10、3.24 × 10m·s。由此可知,有效擴(kuò)散系數(shù)隨溫度的升高而增大。通過Arrhenius方程建立水分有效擴(kuò)散系數(shù)與溫度的關(guān)系,即
式中:D
為方程指前因子,m·s;E
為水分?jǐn)U散活化能,kJ·mol;T
為干燥溫度,K;R
為氣體常數(shù),值為8.314 × 10kJ·(mol·K)。將式(6)兩邊取對(duì)數(shù),得
由式(7)可知,lnD
與1 /T
呈線性關(guān)系,作出lnD
和1 /T
的關(guān)系圖,結(jié)果如圖6所示。根據(jù)圖中直線斜率和截距即可求得E
和D
,D
=3.33×10m·s,E
=31.32kJ·mol。由于大多數(shù)生物質(zhì)廢棄物干燥熱解的活化能在20~400 k J·mol,因此,污泥干燥活化能處在一個(gè)較低水平,說明市政污泥屬于一類相對(duì)易干燥的生物質(zhì)廢棄物。lnDeff與1/T的關(guān)系Fig. 6 Relationship between lnDeff and
本文通過對(duì)市政污泥在不同厚度和溫度下進(jìn)行的干燥實(shí)驗(yàn),引入常見的薄層干燥模型,分析市政污泥在不同厚度和溫度下的干燥特性,并計(jì)算出市政污泥在干燥過程中的水分有效擴(kuò)散系數(shù)和活化能,從而得出以下結(jié)論:
(1)污泥干燥可以近似區(qū)分為3個(gè)階段:升速干燥階段、恒速干燥階段和降速干燥階段。溫度越高,恒速干燥階段范圍越小,越快進(jìn)入降速干燥階段。從能源利用的角度分析,溫度過高會(huì)導(dǎo)致有更多的能量消耗在污泥干質(zhì)上,不利于能量的有效利用。
(2)在多個(gè)工況下利用Midilli模型擬合得到的相關(guān)系數(shù)均在0.99以上,且S
和R
均表現(xiàn)出很好的擬合優(yōu)度,表明可用Midilli模型來預(yù)測(cè)市政污泥干燥過程含水率的變化和干燥時(shí)間。(3)70~130 ℃下市政污泥水分的有效擴(kuò)散系數(shù)D
為5.96 × 10~3.24 × 10m·s,有效擴(kuò)散系數(shù)隨著溫度的升高而增大;(4)市政污泥干燥過程中水分?jǐn)U散的活化能E
為31.32 kJ·mol,處在一個(gè)較低的水平,說明市政污泥是一類相對(duì)易干燥的生物質(zhì)廢棄物。