1. 引言
制浆造纸黑液是制浆过程中的核心副产物,其化学需氧量负荷占制浆废水总量的90%以上,同时富含木质素及多种钠盐等可回收资源。通过蒸发浓缩实现黑液化学品回收与热能利用,是制浆造纸行业实现清洁生产与资源循环的关键技术路径[1]-[3]。对于多效蒸发系统中,随着黑液浓度不断升高,碳酸钠、硫酸钠及草酸钠等无机盐易因超过溶解度极限而结晶析出,其形成的致密结垢层会显著地降低蒸发器的传热性能,使蒸发速率下降,增加能耗,成为制约静态蒸发技术高效应用的主要瓶颈[4]。因此,研究黑液静态蒸发速率的变化规律及结垢特性对于提升蒸发效率和优化运行工况具有重要意义。
对于黑液蒸发与结垢问题,国外研究主要集中于结垢机理及其抑制方法。Frederick和Euhus [4]指出黑液蒸发过程中主要结垢类型为碳酸钠和硫酸钠复盐,并提出通过调控悬浮晶体浓度和操作温差减缓结垢的方法,但未针对静态蒸发条件下蒸发速率的变化进行系统分析。Karlsson [5]-[7]系统分析了工业蒸发器运行数据,评估了钠盐结垢的形成条件及其对蒸发系统的影响,并指出体相晶体的存在可在一定程度上削弱结垢强度,但相关结论尚未在静态蒸发工况下得到验证。
Khafhafera [8] [9]明确了草酸钠及多种钠盐在黑液中的溶解度特性以及温度和总钠含量对结垢阈值的影响规律,但未将其与静态蒸发速率的演化过程相联系。Fatemi [10]基于夹带流反应器建立了黑液结垢倾向的表征方法,证实结垢行为与颗粒含碳量、液体含量及粒径分布密切相关,但该研究未将结垢表征与静态蒸发速率调控相结合。Nwaeri [11] [12]针对钠盐结垢进行研究,发现油酸、妥尔油盐水等添加剂能通过改变其无机盐成核条件和晶体形貌减缓钠盐结垢的生成,但未进一步分析其对蒸发速率的间接影响。
在工艺与工程层面,Cantrell和Straka [13]探讨了黑液蒸发器效数升级对蒸发效率的影响,Kim [14]和Variny [15]提出了机械蒸汽再压缩辅助蒸发工艺以提升系统能效,为不同黑液蒸发器工艺方面提供了新的思路。Karlsson [16]研究了钠盐结垢在降膜蒸发器中的溶解速率,为结垢清洗提供了依据。Chai [17]对黑液中非挥发性组分进行了分析,Singh [18] [19]系统研究了pH、温度及固含量对黑液流变特性的影响,并建立了红外加热条件下黑液蒸发的数学模型,但其与常规静态蒸发工况存在差异。
国内研究方面,刘贤淼[20]通过元素分析和X射线衍射表征明确了黑液中Na、K、S等元素及无机盐组成,但未对静态蒸发速率进行系统量化。张磊[21]针对麻浆黑液蒸发干燥过程优化了工艺参数,分析了温度和进料速率对干燥效果的影响,但对静态蒸发速率变化机制及结垢影响缺乏深入探讨。李雪[22]基于多效蒸发原理建立了工艺计算模型,优化了各效操作参数,但其研究对象为流动蒸发体系,未考虑静态蒸发条件下蒸发行为的影响。贾原媛[23]则针对三相循环流化床黑液蒸发器的防垢和强化传热开展了系统研究。位红卫[24]对碳酸钠和硫酸钠混合盐的结晶行为进行了实验研究,为理解黑液中无机盐析出提供了参考。
现有研究普遍存在结垢研究与静态蒸发速率分析相对分离的问题,同时对静态蒸发过程中速率变化规律及其与结垢关系的研究仍明显不足。基于此以工厂实际制浆造纸黑液为研究对象,明确黑液结垢组成的基础上,围绕静态蒸发速率开展系统研究。通过设计不同温度条件下的黑液静态蒸发实验,并结合XRD等手段对蒸发结垢的组成组分进行表征,观察黑液静态蒸发现象,分析黑液静态蒸发速率的变化规律,并建立反映静态蒸发规律的数学模型。本文旨在通过黑液静态蒸发实验建立黑液静态蒸发的数学模型,实现对黑液静态蒸发质量随时间变化的预测与分析,为后续流动过程中的蒸发、结垢规律提供理论依据及精准预测。
2. 组成分析及静态蒸发实验设计
2.1. 黑液组成分析
黑液的化学组成对其蒸发行为及结垢速率具有显著影响,为明确黑液静态蒸发过程中蒸发速率变化及结垢形成的内在机理,对静态蒸发实验后获得的结垢样品及其煅烧后的颗粒进行X射线衍射(XRD)分析,所用黑液取自制浆造纸生产过程中蒸发系统的现场运行工段,取样过程中采用耐碱密封容器,避免空气中水分蒸发及杂质影响样品,实验前将样品充分混匀后进行取用。
将一定质量的黑液置于耐高温容器中,在设定温度条件下于水浴锅中进行静态蒸发。蒸发结束后,从容器内壁及底部收集形成的结垢样品。
考虑到黑液中有机物与无机盐的影响,将黑液结垢样品以两种方式进行分析,第一种为静态蒸发后直接取用内壁及底部的原始结垢样品,第二种为取用适量结垢后经575℃高温煅烧处理后研磨的粉末状样品。煅烧过程用于去除结垢中的有机物成分,从而明确无机盐的种类及特征。XRD分析主要用于结垢中无机盐的定性分析,而非精确定量,其结果为后续针对主要结垢盐类的定量分析提供依据。
2.2. 黑液中碳酸钠和硫酸钠含量的滴定分析
基于XRD对黑液结垢晶相的分析,碳酸钠(Na2CO3)、硫酸钠(Na2SO4)以及氯化钾(KCl)是黑液静态蒸发结垢煅烧后的主要无机盐成分。虽然在结垢样品中检测到KCl,但由于其溶解度较高,并且在蒸发浓缩过程中不易形成稳定沉积层,主要选取Na2CO3和Na2SO4作为黑液结垢的主要研究对象,两种盐在黑液蒸发过程中溶解度变化显著,是形成管道中黑液结垢的主要无机盐成分。
本文采用滴定法对原始黑液样品中Na2CO3和Na2SO4这两种无机盐的含量进行分析。实验过程中,取定量体积的黑液样品,加入足量的滴定试剂进行操作,并通过酸碱检测仪确保离子完全沉淀,根据沉淀量计算样品中目标无机盐的含量,该方法可为后续静态蒸发实验及结垢分析提供定量依据。
2.3. 黑液静态蒸发实验
通常,开展静态蒸发实验主要是采用恒温加热方式,结合黑液的组成分析发现,其主要成分为水分,针对实际工程中其蒸发温度为80℃到120℃的范围,设计对应温度范围内恒温静态蒸发。实验过程中将一定初始质量的黑液置于开放容器中自然蒸发,并且不施加任何外部气流或搅拌,突出温度和溶液本身组成对蒸发的影响。通过设置不同初始质量和不同的温度条件,研究黑液在不同工况下的静态蒸发特征。
在蒸发过程中定时测量黑液质量变化,获得黑液质量随时间变化的数据。同时,对蒸发过程中容器内壁及底部的结垢形态进行观察,记录结垢出现的时间及其随蒸发过程的变化特征。
通过对黑液静态蒸发的研究,可以明确黑液中水分蒸发速率与无机盐浓缩速率之间的关系,从而判断Na2CO3和Na2SO4何时达到过饱和状态并开始析出结晶,该过程是结垢形成的前提条件,并为后续黑液静态蒸发速率模型的建立与验证提供可靠的实验依据。
3. 结果与分析
3.1. 黑液组成分析
黑液的组成对其静态蒸发速率及结垢具有显著影响,通过多次的黑液结垢成分XRD表征结果(如图1)分析发现,黑液结垢中Na2CO3、Na2SO4及KCl为其主要无机盐组成。
Figure 1. XRD patterns of black liquor scaling before and after calcination under static evaporation conditions
图1. 黑液静态蒸发后结垢煅烧前后的XRD谱图
结合黑液的组成特点,黑液锻烧前结垢中含有较多的有机物成分,这些有机物在蒸发浓缩过程中易发生黏附和聚集,并包覆无机盐晶体,从而导致XRD衍射峰弱化。表征结果显示(图1中的红线谱图),煅烧前的XRD谱图整体衍射峰强度较弱,特征晶峰不够清晰,这表明黑液蒸发后的结垢中含有大量无定形或低结晶度物质,对无机盐晶相的衍射信号产生了明显的屏蔽作用。无机盐的结晶作为黑液结垢的主要因素,而有机物在黑液蒸发中主要起到变黏以及吸附的作用,其对黑液结垢的产生有影响但无决定性作用,因此,单独对黑液蒸发后结垢进行XRD分析难以准确识别其中的无机盐种类。
为进一步明确结垢中的无机盐组成,对黑液蒸发后的结垢样品以575℃进行高温煅烧处理去除有机物,结果如图1中的黑线谱图,与煅烧前的黑液结垢对比,煅烧后结垢的衍射峰显著增强,背景信号明显降低,晶相特征更加清晰。通过对衍射峰位置的比对,煅烧后结垢样品中存在Na2CO3、Na2SO4以及KCl的特征衍射峰。由于KCl的溶解度非常高并且受温度影响较弱,KCl并不容易产生结晶沉积形成结垢,并且其在黑液中的成分较少,因此可以表明形成黑液蒸发结垢的无机盐主要由Na2CO3、Na2SO4构成。
3.2. 黑液结垢主要无机盐的静态滴定结果分析
基于XRD对黑液结垢晶相的分析,碳酸钠和硫酸钠为黑液结垢的主要无机盐组分,如3.1。通过对黑液的滴定分析,进行黑液中碳酸钠和硫酸钠组成的进一步确认,详尽结果如表1,为20 ml黑液中碳酸钠和硫酸钠的含量测定结果。
Table 1. Mass of 20 ml black liquor containing sodium sulfate and sodium carbonate
表1. 20 ml黑液硫酸钠和碳酸钠质量
质量 |
第一组 |
第二组 |
第三组 |
第四组 |
第五组 |
第六组 |
黑液质量/g |
19.165 |
19.085 |
19.238 |
19.21 |
19.151 |
19.217 |
碳酸钠质量/g |
0.326 |
0.32 |
0.331 |
0.33 |
0.319 |
0.323 |
硫酸钠质量/g |
0.101 |
0.103 |
0.11 |
0.106 |
0.1 |
0.105 |
通过六组重复实验,对表1中碳酸钠和硫酸钠的质量数据进行相对标准偏差(RSD)计算。碳酸钠的平均值为0.325 g,标准差为0.005,相对标准偏差为1.54%。硫酸钠的平均值为0.104 g,标准差为0.0037,相对标准偏差为3.55%。二者相对标准偏差均小于5%,实验误差较小,表明测定的碳酸钠和硫酸钠的含量较为精准,综合分析表明,黑液中碳酸钠的质量分数为1.69%,硫酸钠的质量分数为0.54%。
3.3. 黑液静态蒸发速率及结垢特征结果分析
3.3.1. 黑液静态蒸发过程中的形貌演化特征研究
以黑液结垢的组成以及黑液中主要无机盐的含量为基础,如3.1和3.2,研究黑液在静态蒸发过程中结垢行为的演化规律及其对蒸发速率变化的影响,本节在不同温度条件下开展黑液静态蒸发实验,重点分析蒸发过程中黑液宏观形貌变化、结垢形成特征,将黑液静态蒸发过程划分为自由蒸发阶段、浓缩结垢演化过渡阶段和残余质量受限阶段。
由于黑液中水分含量较高,黏度较小,蒸发过程主要受传热与水分扩散控制,此阶段内黑液质量随时间近似呈线性下降,蒸发速率维持在较高数值,且在短时间内变化幅度较小。不同温度下的曲线斜率差异主要来源于温度的影响,液面仅存在一层较薄的表面膜,对实际蒸发影响较小。随着温度由80℃提升至120℃,单位时间内黑液质量损失增大,说明此时温度是影响蒸发初期速率的主导因素,如图2,此时为黑液自由蒸发阶段。
随着蒸发的持续进行,黑液表面开始变粘稠,液面开始由流动状态向黏稠状态转变,此时黑液表面生成较厚的膜状结垢,其边缘开始贴附于容器壁,影响黑液的蒸发速率,如图3(a),为黑液浓缩结垢演化过渡阶段的前期。随后黑液继续蒸发,其表面膜不断加厚并趋于致密,呈现出更为致密的纹理,粘稠度更高,流动性更低,黑液中的无机盐开始结晶生成结垢,影响黑液的蒸发,逐渐形成至图3(b)形貌,此时为黑液浓缩结垢演化过渡阶段的后期。
Figure 2. Black liquor free evaporation stage
图2. 黑液自由蒸发阶段
(a) (b)
Figure 3. Transition stage of scale formation during black liquor concentration
图3. 黑液浓缩结垢演化过渡阶段
当黑液中大量无机盐已经超过溶解度极限,在器壁或表面析出并形成较为坚硬的结垢层,黑液表面形成密集褶皱,完全失去流动性,此时黑液质量变化趋于平缓,曲线逐渐接近稳定残余质量,如图4为黑液残余质量受限阶段。体系中存在不可进一步蒸发的残余物质,该残余质量在不同初始黑液量下存在明显差异,说明残余质量不仅与溶液组成有关,还受到初始质量的显著影响。
Figure 4. Stage of limited residual quality of black liquor
图4. 黑液残余质量受限阶段
3.3.2. 黑液静态蒸发速率以及质量变化规律分析
基于黑液静态蒸发三个蒸发阶段的演化过程,如3.3.1,在80℃、100℃和120℃条件下,对黑液进行静态蒸发实验,记录蒸发过程中黑液质量随时间的变化情况。
在不同温度条件下,黑液质量以及黑液蒸发速率随时间的变化均呈现出明显的非线性衰减特征,在黑液自由蒸发阶段,各温度条件下质量变化速率较快,表现出近似线性的下降趋势,此时黑液蒸发速率也呈现出较高的数值,但不同温度下蒸发速率存在显著差异,黑液蒸发温度越高,黑液的蒸发速率越大,80℃、100℃和120℃下黑液的蒸发速率分别接近于11 g/h、15 g/h、25 g/h。如图5为黑液在80℃、100℃和120℃条件下的质量以及蒸发速率随时间变化的曲线。
随着蒸发时间的延长,在黑液蒸发速率于5 g/h时进入浓缩结垢演化过渡阶段,黑液质量下降速率逐渐减缓,以近似反曲线的趋势蒸发至末期,并且如图5(c)在高温条件下,该趋势更加明显。最后蒸发进入黑液残余质量受限阶段,此时黑液质量逐渐趋于平缓,接近于残余质量,如图5(a)和图5(b),二者为不同初始质量黑液的蒸发过程,图5(a)在初始质量更高的条件下,黑液残余质量也越高,此时三种温度下黑液的蒸发速率均小于0.1 g/h。
(a) 80℃ (b) 100℃
(c) 120℃
Figure 5. Graph of black liquor quality and evaporation rate changes over time at 80˚C, 100˚C, and 120˚C
图5. 80℃、100℃、120℃下黑液质量以及蒸发速率随时间变化图
3.4. 黑液静态蒸发数学模型的建立与分析
黑液静态蒸发过程表现出显著的阶段性特征,其蒸发速率的变化与黑液形貌演化、黏度增长及结垢形成过程关联,如3.2与3.3黑液蒸发现象及规律的分析。针对其阶段性特征,建立考虑残余质量的黑液静态蒸发数学模型,以分析不同阶段黑液的蒸发过程,并基于上述实验结果,并对模型的适用性和预测能力进行验证。
3.4.1. 模型建立
黑液静态蒸发质量以指数形式的变化,并考虑到黑液蒸发的残余质量,如图5黑液质量随时间变化的趋势,作出如下合理假设:蒸发过程中水分不断减少,溶质不挥发;黑液中存在不可蒸发的残余质量,残余质量受初始黑液质量影响;温度对蒸发速率具有显著影响。
基于假设,并通过质量时间变化曲线引入与初始黑液质量相关的残余质量项,考虑温度对黑液蒸发的影响,黑液质量随时间的变化可表示为(1)指数模型:
(1)
其中:
:时间t的黑液质量(g);
:初始黑液质量(g);
:不可蒸发残余质量(g);
:蒸发系数(
);
:蒸发时间(h)。
在实验温度下不同初始黑液质量下的最终残余质量存在明显差异,考虑到初始质量越大,最终残余质量越高的数值关系,将不可蒸发残余质量表示为初始黑液质量的线性函数,如(2):
(2)
其中:
:不可蒸发质量比例系数
:与溶液组成相关的常数项。
考虑到不同温度对黑液蒸发速率的影响,引入蒸发系数,蒸发系数
采用指数型关系,如(3):
(3)
其中:
:指前因子(
);
:表观蒸发活化能(
);
:气体常数,8.314(
);
:绝对温度(K)。
对上述假设搭建的质量变化模型,根据3.3.2的实验数据进行非线性回归,得到对应的
、
、
,得出公式模型如(4):
(4)
该模型拟合采用决定系数R2进行评价,通过本研究中对80℃、100℃和120℃的实验数据与模型的拟合并计算决定系数R2,得出各温度条件下实验模型与蒸发速率模型的R2均大于0.90,表明模型具有较好的拟合度,拟合曲线如下图6:
(a) 80℃ (b) 100℃
(c) 120℃
Figure 6. Experimental data and model fitting plots at 80˚C, 100˚C, and 120˚C
图6. 80℃、100℃、120℃下实验数据与模型拟合图
3.4.2. 模型检验
基于获得的黑液质量模型,如3.4.1,对90℃、110℃下的黑液静态蒸发规律进行预测分析,计算后与实际数据进行对比,并进行相对误差分析,数据如表2。虽然一定时间内相对误差较高,但考虑到黑液蒸发模型预测整体蒸发趋势,为评估整体数据的精度,得到平均相对误差分别为1.6%、2.0%,平均相对误差均小于5%,验证了模型在80℃到120℃温度区间的准确性。
Table 2. Comparison of actual and predicted black liquor quality at 90˚C and 110˚C for different durations
表2. 90℃、110℃不同时间时黑液实际质量与预测质量对比
90℃时间/h |
实际质量/g |
预测质量/g |
相对误差/% |
110℃时间/h |
实际质量/g |
预测质量/g |
相对误差/% |
0 |
79.718 |
79.718 |
0.0 |
0 |
71.115 |
71.115 |
0.0 |
0.5 |
73.031 |
70.514 |
3.4 |
0.5 |
60.039 |
57.866 |
3.6 |
1 |
66.711 |
62.576 |
6.2 |
1 |
48.997 |
47.558 |
2.9 |
1.5 |
60.505 |
55.730 |
7.9 |
1.75 |
35.837 |
36.222 |
−1.1 |
2 |
54.440 |
49.825 |
8.5 |
2.5 |
24.267 |
28.443 |
−17.2 |
2.5 |
48.183 |
44.733 |
7.2 |
3 |
20.789 |
24.666 |
−18.7 |
3 |
42.454 |
40.342 |
5.0 |
4 |
17.547 |
19.441 |
−10.8 |
4 |
32.072 |
33.288 |
−3.8 |
5 |
16.110 |
16.279 |
−1.0 |
5 |
24.261 |
28.042 |
−15.6 |
6 |
15.189 |
14.364 |
5.4 |
6 |
21.051 |
24.139 |
−14.7 |
7 |
14.532 |
13.205 |
9.1 |
7 |
19.179 |
21.237 |
−10.7 |
8 |
14.064 |
12.504 |
11.1 |
8 |
18.164 |
19.078 |
−5.0 |
9 |
13.664 |
12.079 |
11.6 |
9 |
17.370 |
17.472 |
0.6 |
10 |
13.343 |
11.822 |
11.4 |
10 |
16.699 |
16.278 |
2.5 |
11 |
13.135 |
11.667 |
11.2 |
11 |
16.159 |
15.390 |
4.8 |
13 |
12.925 |
11.515 |
10.9 |
12 |
15.788 |
14.729 |
6.7 |
17 |
12.416 |
11.440 |
7.9 |
16 |
14.822 |
13.398 |
9.6 |
21 |
11.942 |
11.430 |
4.3 |
21 |
13.979 |
12.945 |
7.4 |
23 |
11.616 |
11.429 |
1.6 |
23 |
13.735 |
12.885 |
6.2 |
25 |
11.472 |
11.428 |
0.4 |
27 |
13.441 |
12.834 |
4.5 |
27 |
11.399 |
11.428 |
−0.3 |
30 |
13.264 |
12.820 |
3.3 |
29 |
11.283 |
11.428 |
−1.3 |
构建的黑液质量模型是在实验条件下得到的经验模型,能够较好反映80℃到120℃之间黑液静态蒸发过程中质量随时间变化的整体规律。该模型参数具有明确物理意义,其残余质量项反映了结垢与不可蒸发物质的存在,温度影响通过蒸发系数体现,模型形式简洁,参数可由实验直接获得,适合用于一定温度范围内的工程预测。但同时,当黑液在蒸发后期结垢极为严重或溶液组成发生显著变化时,模型预测精度可能受到一定影响,仍需在更宽工况范围内进一步修正和完善。
4. 结论与展望
4.1. 结论
为研究制浆造纸黑液静态蒸发规律及结垢组成,本文采用X射线衍射、滴定法以及静态蒸发实验的方法,并在此基础上建立了静态蒸发速率预测模型,得到如下结论:
XRD表征结果表明,黑液静态蒸发结垢主要由碳酸钠和硫酸钠构成,二者为黑液蒸发结垢的主要无机盐成分。化学滴定分析结果显示,黑液中碳酸钠和硫酸钠的质量分数分别约为1.69%和0.54%,测定结果的相对标准偏差分别为1.54%、3.55%。
静态蒸发实验结果表明,在80到120℃温度范围内,黑液蒸发过程具有显著的阶段性特征,可明确划分为自由蒸发阶段、浓缩结垢演化过渡阶段及残余质量受限阶段。自由蒸发阶段黑液蒸发速率维持在较高数值,80℃、100℃、120℃时黑液的蒸发速率分别为11 g/h、15 g/h、25 g/h,随着水分大量蒸发,进入浓缩结垢演化过渡阶段,蒸发速率降至5 g/h,最后由于黑液中大量结垢析出,进入残余质量受限阶段,其蒸发速率以小于0.1 g/h的速率变化,接近稳定残余质量。
基于黑液静态蒸发实验数据,构建了考虑时间、温度及初始质量的黑液静态蒸发模型,80℃、100℃和120℃条件范围内的R2分别为0.980、0.939、0.993,该模型能够较好描述黑液在静态蒸发过程中的质量变化规律,拟合效果较好。对90℃和110℃条件下的质量变化规律进行预测,将实验数据与模型数据进行平均相对误差分析,得到相对误差分别为1.6%、2.0%,精度较高,可用于黑液静态蒸发速率的预测与分析,并为预测不同蒸发时间下黑液的结垢量提供科学的方法。
4.2. 展望
尽管本文对黑液静态蒸发规律进行了较为系统的研究,但仍存在一定局限性。实验主要针对特定来源和组分范围内的黑液样品,未系统考察黑液组成差异如无机盐比例对蒸发速率及结垢行为的影响,而且黑液蒸发分析是基于恒温静态工况,与管流过程中的实际生产工况存在一定差异,流动过程更加复杂,其结垢速率及特征有待进一步深入探讨分析。
未来研究可以通过系统改变黑液中无机盐的组成比例,深入分析不同组分条件下黑液蒸发速率及结垢行为的差异规律,或者引入动态蒸发条件,研究流动状态对黑液蒸发及结垢形成的影响机制。通过上述研究,可以为制浆造纸黑液蒸发过程的运行优化提供更加完善的理论依据。
NOTES
*通讯作者。