2. 嘉兴学院 建工学院,浙江 嘉兴 314001
2. College of Civil Engineering and Architecture, Jiaxing University, Jiaxing 314001, China
污水处理处置一直是关系着人们生活和工业发展的重大问题。根据2018年12月住房和城乡建设部统计数据[1]显示,2017年我国污水处理量达15 743万m3·日−1。如果按照污水处理厂平均以负荷率80%运行,每万t污水产生6 t含水率(质量分数)为80%的污泥来算,我国每天产生污泥超过7.55万t。
污泥如果得不到妥善处理,污泥中的重金属、细菌病原体等微生物、有机污染物将会对环境造成难以预料的损害,甚至危害到人类及其他动植物的生存。由于污泥焚烧能彻底消灭病原体等微生物,使有机物碳化,且所释放的能量可供热、发电,焚烧渣可制作附属产品,真正使污泥实现减量化、稳定化、无害化和资源化的目的,因而被认为是最有发展前景的污泥处理方法之一[2]。一般经过机械脱水的污泥含水率仍在80%以上,此时污泥的热值接近于0,无法直接焚烧,需经干燥预处理将污泥含水率降到50%,污泥的热值才能利用起来[3]。而且,干燥过程使污泥含水量降低,减小了污泥的容量,便于污泥的后续处理。
很多国内外学者对污泥干燥展开了相关研究,大多是围绕污泥成分分布[4]、不同形状污泥的干燥特性[5]、不同干燥工艺对干燥的影响[6]等方面去研究,缺乏对污泥干燥过程中深层次动力学机理的研究。本实验研究污泥厚度和温度对污泥干燥特性的影响,通过引入薄层干燥动力学模型,对污泥干燥进行动力学模拟,得出污泥干燥的最优模型和有效扩散系数,为污泥干燥工艺提供新的技术参数和依据。
1 实验材料与方法 1.1 实验材料本次实验污泥样品取自嘉兴市某污水厂的脱水市政污泥,呈深黑色,黏稠状,含水率达85%。污泥的主要成分C、H、O、N、Cl和苯乙烯的质量分数分别为18.34%、3.82%、16.27%、1.91%、32.21%和1.69%。
1.2 实验装置和方法实验中采用精度为0.01 g的德国sartorius电子天平并通过VB软件编写程序对污泥质量进行30″一次数据采集。将污泥平铺成饼状置于温度可控的电加热板上。首先,控制加热温度为110 ℃,分别对3、6、9、12 mm厚度的污泥进行加热,考察污泥在不同厚度下的干燥情况;然后,控制污泥厚度为9 mm,分别设定加热温度为70、90、110、130 ℃,通过绘制污泥的失水曲线,探究不同温度下污泥的干燥特性。
1.3 数据处理方法污泥的含水率随时间变化[7],即
| $\qquad U = \frac{{({M_t} - {M_{\rm{e}}})}}{{({M_0} - {M_{\rm{e}}})}}$ | (1) |
式中:U为污泥的含水率,为无因次量;
由于
| $\qquad U = \frac{{{M_t}}}{{{M_0}}}$ | (2) |
污泥的失水效果可用干燥速率表示,即
| $\qquad V = - \frac{{({M_{t + \Delta t}} - {M_t})}}{{\Delta t}}$ | (3) |
式中:V为污泥的干燥速率,g·(g·min)−1;
温度110 ℃下不同厚度污泥的含水率随时间的变化如图1所示。
|
图 1 不同厚度的污泥含水率随时间的变化 Fig.1 Relationship between moisture content in the sludge and time under different thicknesses |
从图1中可看出,污泥越薄,干燥至平衡含水率的时间越短。这是因为,污泥越薄,水分从污泥内部迁移到污泥表面的路径越短,传质阻力也越小,因此所需时间越短,水分蒸发得越快。当污泥厚度从12 mm减小至3 mm时,干燥至平衡含水率所用的时间从约200 min缩短至接近30 min。因此,将污泥铺平使其变得更薄,会大大缩短污泥干燥时间。
2.2 不同温度下污泥的干燥特性不同温度下厚度为9 mm的污泥含水率随时间的变化如图2所示。
|
图 2 不同温度下污泥含水率随时间的变化 Fig.2 Relationship between moisture content in the sludge and time under different temperatures |
从图2中可看出,在污泥厚度相同的情况下,提高干燥温度,污泥干燥至平衡含水率的时间越短。这是因为提高干燥温度,即增大了污泥内部的传热温差,在传热温差的驱动作用下,污泥内的水分迁移加快。当温度从70 ℃上升到130 ℃时,污泥干燥至平衡含水率的时间从200 min缩短至60 min。
图3为不同温度下干燥速率和时间的关系。图4为不同温度下干燥速率和含水率的关系。结合图3、4可看出,当污泥刚开始被加热时,其表面含有大量的自由水,随着时间的增长,污泥表面温度不断上升,水分不断蒸发,干燥速率明显上升,此为升速干燥阶段。当加热到一定时间,污泥的自由水含量减少,间隙水和剩余的自由水变成了主要干燥对象,污泥内部的水分向表面迁移的速率接近表面水分的蒸发速率,表现为近似恒速干燥,而且温度越高,恒速干燥阶段范围越小。当温度达到130 ℃时,恒速干燥阶段已不明显。恒速干燥阶段后,污泥开始出现干区,这增加了传热传质阻力,污泥中的水分为剩下的部分间隙水和结合水;这部分水分难以去除,因此干燥速率降低,此时为降速干燥阶段。最大干燥速率出现在恒速干燥阶段,当温度从70 ℃上升到130 ℃时,相应的最大干燥速率从0.025 87 g·(g·min)−1上升到0.081 58 g·(g·min)−1。
|
图 3 不同温度下干燥速率和时间的关系 Fig.3 Relationship between drying rate and time under different temperatures |
|
图 4 不同温度下干燥速率和含水率的关系 Fig.4 Relationship between drying rate and moisture content under different temperatures |
从图4中还可知,在对应的含水率区间温度越高,恒速干燥阶段范围越小,温度升高使得污泥更快地进入降速干燥阶段,造成更多的能量用来加热污泥的干区。这表明,若将污泥干燥至相同的含水率,从能源利用的角度分析,温度过高会导致有更多的能量消耗在污泥干质上,不利于能源的有效利用。因此,在满足生产的前提下,干燥温度应控制在合理的范围。
2.3 干燥模型分析引入薄层干燥模型对污泥干燥曲线进行模拟分析,模型表达式如表1所示,其中:k、k1和k2均为模型的经验系数;a、b、c和y均为模型的经验常数。
|
|
表 1 薄层干燥模型 Table 1 Models of thin layer drying |
以9 mm薄层污泥在70、90、110和130 ℃下的干燥曲线为例,应用Origin软件分别利用每个薄层干燥模型进行数据拟合,得出模型参数,并用相关系数R2、标准差SD和残差平方和RSS来评价模型的拟合优度。R2越接近1,SD和RSS越小,说明模型的拟合效果越好。
表2~5为每个薄层干燥模型的拟合结果。由此可看出:Midilli模型的R2均超过0.99,在所有模型中最接近1;SD和RSS在所有模型中也最小。说明Midilli模型较适用于污泥薄层干燥的模拟。
|
|
表 2 70 ℃污泥薄层干燥模型拟合结果 Table 2 Fitting results of thin layer drying models of the sludge at 70 ℃ |
|
|
表 3 90 ℃污泥薄层干燥模型拟合结果 Table 3 Fitting results of thin layer drying models of the sludge at 90 ℃ |
|
|
表 4 110 ℃污泥薄层干燥模型拟合结果 Table 4 Fitting results of thin layer drying models of the sludge at 110 ℃ |
|
|
表 5 130 ℃污泥薄层干燥模型拟合结果 Table 5 Fitting results of thin layer drying models of the sludge at 130 ℃ |
图5给出了Midilli模型拟合值和实验值的对比结果,其数据基本上在斜率为45°的直线上,说明Midilli模型在70~130 ℃下能够很好地对市政污泥干燥进行模拟。由于模型能够描述含水率随时间的变化,因此对于估算污泥干燥程度有重要的指导意义。
|
图 5 含水率实验值与Midilli模型拟合值的比较 Fig.5 Comparison between experimental values and fitting ones of moisture content by Midilli |
污泥中的水分迁移是一个从内向外扩散的过程,根据Fick扩散定律可以推导出包含有效水分扩散系数的薄层干燥理论方程[10],即
| $ U = \frac{8}{{{ {\text{π}} ^2}}}\sum\limits_{n = 0}^\infty \left \{{\frac{1}{{{{(2n + 1)}^2}}}\exp \left[ - \frac{{{(2n + 1)^2}{{\text{π}} ^2}{D_{\rm{eff}}}t}}{{4{L^2}}}\right]} \right \} $ | (4) |
式中:Deff为有效扩散系数,m2·s−1;L为污泥厚度,实验中取L = 0.009 m;n为实验采样次数。
由于干燥时间较长,取式(4)第1项并两边取对数将其简化为[16]
| $\qquad\ln U = \ln \frac{8}{{{{\text{π}} ^2}}} - \frac{{{{\text{π}} ^2}{D_{\rm{eff}}}t}}{{4{L^2}}}$ | (5) |
通过式(5)可以发现,由lnU与t的线性关系可以计算有效扩散系数Deff。经计算,70、90、110、130 ℃下有效扩散系数分别为5.96 × 10−9、1.04 × 10−8 、1.54 × 10−8 、3.24 × 10−8 m2·s−1。由此可知,有效扩散系数随温度的升高而增大。
通过Arrhenius方程建立水分有效扩散系数与温度的关系[17],即
| $\qquad{D_{{\rm{eff}}}} = {D_0}\exp \left( - \frac{{{E_{\rm{a}}}}}{{RT}} \right)$ | (6) |
式中:D0为方程指前因子,m2·s−1;Ea为水分扩散活化能,kJ·mol−1;T为干燥温度,K;R为气体常数,值为8.314 × 10−3 kJ·(mol·K)−1。
将式(6)两边取对数,得
| $\qquad\ln {D_{{\rm{eff}}}} = \ln {D_0} - \frac{{{E_{\rm{a}}}}}{{RT}}$ | (7) |
由式(7)可知,lnDeff与
|
图 6 lnDeff与
|
本文通过对市政污泥在不同厚度和温度下进行的干燥实验,引入常见的薄层干燥模型,分析市政污泥在不同厚度和温度下的干燥特性,并计算出市政污泥在干燥过程中的水分有效扩散系数和活化能,从而得出以下结论:
(1)污泥干燥可以近似区分为3个阶段:升速干燥阶段、恒速干燥阶段和降速干燥阶段。温度越高,恒速干燥阶段范围越小,越快进入降速干燥阶段。从能源利用的角度分析,温度过高会导致有更多的能量消耗在污泥干质上,不利于能量的有效利用。
(2)在多个工况下利用Midilli模型拟合得到的相关系数均在0.99以上,且SD和RSS均表现出很好的拟合优度,表明可用Midilli模型来预测市政污泥干燥过程含水率的变化和干燥时间。
(3)70~130 ℃下市政污泥水分的有效扩散系数
(4)市政污泥干燥过程中水分扩散的活化能
| [1] |
中华人民共和国城乡和住房建设部. 中国城乡建设统计年鉴2017 [M]. 北京: 中国统计出版社, 2018 .
|
| [2] |
郭志光, 马明磊, 刘斌, 等. 污泥处理处置技术研究进展[J]. 河北地质大学学报, 2019, 42(1): 59-63. |
| [3] |
麻红磊. 城市污水污泥热水解特性及污泥高效脱水技术研究[D]. 杭州: 浙江大学, 2012.
|
| [4] |
张卫军. 城市污水污泥热干化特性及能耗研究[D]. 重庆: 重庆大学, 2009: 66 - 70.
|
| [5] |
马学文, 翁焕新. 温度与颗粒大小对污泥干燥特性的影响[J]. 浙江大学学报: 工学版, 2009, 43(9): 1661-1667. DOI:10.3785/j.issn.1008-973X.2009.09.020 |
| [6] |
DENG W Y, SU Y X. Experimental study on agitated drying characteristics of sewage sludge under the effects of different additive agents[J]. Journal of Environmental Sciences, 2014, 26(7): 1523-1529. DOI:10.1016/j.jes.2014.05.019 |
| [7] |
杨小文, 杜英豪. 国外污泥干化技术进展[J]. 给水排水, 2002, 28(2): 35-36. DOI:10.3969/j.issn.1002-8471.2002.02.011 |
| [8] |
THAKOR N J, SOKHANSANJ S, SOSULSKI F W, et al. Mass and dimensional changes of single canola kernels during drying[J]. Journal of Food Engineering, 1999, 40(3): 153-160. DOI:10.1016/S0260-8774(99)00042-4 |
| [9] |
AYENSU A. Dehydration of food crops using a solar dryer with convective heat flow[J]. Solar Energy, 1997, 59(4−6): 121-126. DOI:10.1016/S0038-092X(96)00130-2 |
| [10] |
AKGUN N A, DOYMAZ I. Modelling of olive cake thin-layer drying process[J]. Journal of Food Engineering, 2005, 68(4): 455-461. DOI:10.1016/j.jfoodeng.2004.06.023 |
| [11] |
AKPINAR E K, BICER Y, YILDIZ C. Thin layer drying of red pepper[J]. Journal of Food Engineering, 2003, 59(1): 99-104. DOI:10.1016/S0260-8774(02)00425-9 |
| [12] |
HENDERSON S M, PABIS S. Grain drying theory II: temperature effects on drying coefficients[J]. Journal of Agricultural Engineering Research, 1961(6): 169-174. |
| [13] |
TOĞRUL İ T, PEHLIVAN D. Mathematical modelling of solar drying of apricots in thin layers[J]. Journal of Food Engineering, 2002, 55(3): 209-216. DOI:10.1016/S0260-8774(02)00065-1 |
| [14] |
CELMA A R, ROJAS S, LÓPEZ F, et al. Thin-layer drying behaviour of sludge of olive oil extraction[J]. Journal of Food Engineering, 2007, 80(4): 1261-1271. DOI:10.1016/j.jfoodeng.2006.09.020 |
| [15] |
SHARAF-ELDEEN Y I, BLAISDELL J L, HAMDY M Y. A model for ear corn drying[J]. Transactions of the ASAE, 1980, 23(5): 1261-1265. DOI:10.13031/2013.34757 |
| [16] |
LERTWORASIRIKUL S. Drying kinetics of semi-finished cassava crackers: a comparative study[J]. LWT-Food Science and Technology, 2008, 41(8): 1360-1371. DOI:10.1016/j.lwt.2007.09.009 |
| [17] |
MADAMBA P S, DRISCOLL R H, BUCKLE K A. The thin-layer drying characteristics of garlic slices[J]. Journal of Food Engineering, 1996, 29(1): 75-97. DOI:10.1016/0260-8774(95)00062-3 |
| [18] |
陈允轩. 焚烧炉条件下典型城市生活垃圾干燥过程的试验研究[D]. 北京: 北京交通大学, 2007.
|
2021, Vol. 37

