脱硫废水烟气蒸发过程中氯离子气固相分配模型及实验研究

2019-05-08马双忱万忠诚张净瑞陈嘉宁

煤炭学报 2019年4期

马双忱，柴晋，万忠诚，武凯，张净瑞，陈嘉宁

(1.华北电力大学环境科学与工程学院，河北保定 071003; 2.盛发环保有限公司，福建厦门 361000)

2015-04-16，国务院发布了《水污染防治行动计划》[1]，国家将强化对各类水污染的治理力度，脱硫废水为电厂的终端废水，水质恶劣，成分复杂，脱硫废水零排放成为电厂水污染深度治理的首要任务;2017-06-01，环保部印发《火电厂污染防治可行性技术指南》[2]，明确提出了烟气蒸发干燥、热法结晶工艺进行废水近零排放的技术路线。

在众多零排放技术路线中，高温烟气蒸发技术受到普遍关注[3]，该技术从SCR后、空预器前抽取干燥过程所需的热烟气，经雾化器雾化的脱硫废水以细小液滴的形式和热烟气在烟道内相互接触、传质传热，完成干燥过程，蒸发产物与烟气一起返回原烟道进入除尘器，蒸发产物随粉尘一起被捕集，但脱硫废水蒸发过程中会释放部分气态HCl，而除尘器对气态物质去除效率极低[4-6]，该部分气态HCl将进入脱硫系统，无疑会影响脱硫系统的氯平衡和水平衡，当氯挥发量过大时，会增加脱硫废水排放量;另外，浆液中Cl-浓度过高对石膏品质、脱硫系统腐蚀和脱硫效率等也会产生不利影响[7-9]。在废水零排放的发展趋势下，如何量化及预测不同工况条件下氯挥发量，建立氯离子气固相分配模型成为了脱硫废水高温烟气蒸发技术的核心问题。

关于氯离子的相间分配模型，目前主要有活度系数法和状态方程法2种研究手段[10-14]，国外学者主要利用活度系数法计算氯的气液相分配，定义气液相分配系数KD，即

(1)

其中，mclvap为气相Cl-浓度;mclliq为液相Cl-浓度;γclliq为氯离子液相活度系数。J M SIMOMSON[15]探究了50～350 ℃温度条件下氯离子的气液相分配规律，首次提出式(2)氯离子的气液相分配经验模型，同时提出氯离子的气液相分配受到溶液温度、pH、氯离子含量、含盐量及盐成分等多方面因素影响[16]。

(2)

其中，T为蒸发温度;ρw为水的密度;A，B，C和D为分配常数。M S GRUSZKIEWICZ[17-19]和LIU C K[20]将上述理论研究与电厂锅炉水高温高压工况下的溶解携带理论[21-22]相结合，得出式(3)关于温度的简化模型，即氯离子气液相分配系数的对数与温度的倒数成正比。

(3)

式中，Tc,HCl和ρc,HCl分别为HCl的临界温度和临界密度;ρw为水的密度;A和B为分配常数。

一般状态下氯离子通常是以HCl和NH4Cl的形式挥发到气相，在低氨浓度时，以HCl的形式为主[18]，脱硫废水中NH4+平均含量为9.39 mg/L，属于低氨水质;Otakar JONAS[23]研究了纯水的pH随温度的变化，纯水的pH随温度的升高呈现出降低的趋势，在430 ℃时达到最低pH为5.6，也就是说，水在高温下有“酸化”的倾向，当水中含有Cl-时，水解离生成的H+易与Cl-结合，以气态HCl的形式挥发溢出，另外NaCl,CaCl2,MgCl2等可能在高温条件下水解生成气态HCl。氯离子挥发机理总结如下[18,23-24]:

(1)

(2)

(3)

(4)

反应(1)为氯离子气相转化的主要形式，反应(2)～(4)可能在蒸发温度较高时发生，且反应强度不高。

笔者以脱硫废水蒸发过程为研究对象，将传质分离过程中涉及的气液相平衡理论和状态方程法应用于复杂的混盐体系，并根据实验数据加以修正，得出适用于180～380 ℃温度条件且包含关键影响因素的氯离子气固相分配系数计算模型，实现了对脱硫废水蒸发过程中氯挥发量的量化及预测，为脱硫废水烟气蒸发技术的应用提供了重要的数据支持，有重要的理论和实用价值。

1 氯离子气固相分配模型

1.1 模型简化及假设

图1为不同体系氯离子的气固相分配机理示意图，模型简化及假设如下:

(1)由上述研究进展可知，在低氨浓度时，氯离子主要以HCl的形式挥发到气相，且纯水pH随温度升高而降低，水解离的H+易与Cl-结合，以气态HCl的形式溢出，脱硫废水属于低氨酸性水质，所以，理论模型计算时定义理想状态，假设理想状态下氯离子全部以HCl的形式挥发到气相;

图1 不同体系氯离子气固相分配机理Fig.1 Gas-solid phase partition mechanism of chloride ions in different systems

(2)将氯离子在理想状态下的挥发简化为HCl-纯水体系(以下简称HCl-H2O体系)，如图1(a)所示，在蒸发过程中HCl和H2O的传质传热分离过程，在探究脱硫废水等复杂体系时，将理想状态替换为实际的Cl--混盐体系(以下简称Cl--Salt体系)，如图1(b)所示，并对HCl-H2O体系得出的理想模型加以修正;

(3)蒸发过程中传质传热分离只在气液相接触面处进行，忽略其他散热损失。

1.2 模型计算方法

气液相平衡时，组分i在气液两相的逸度相等[25-27]，即

(4)

(5)

式中，p为总压；yi为气相组分i的摩尔分率，同理可得液相中组分i的逸度系数:

(6)

其中,xi为液相中组分i的摩尔分率。将式(5)，(6)代入式(4)中得到

(7)

则气液相分配系数可表示为

(8)

式中，逸度系数可通过状态方程来计算，常用的以V和T作为独立变量，导出计算逸度系数的普遍化方程:

其中，Vm,t为气(液)相混合物的总体积;Zm为气(液)相混合物的压缩因子;ni为组分i的物质的量。但通过式(9)并不能直接得到数值解，本文采用Van Der Waals提出的实际气体状态方程[28-29]:

(10)

Van Der Waals方程，引入了a,b两项气相混合物范德华常数，a是分子间引力对p的校正，b为分子本身的体积校正，结合经典的VDW混合规则，可以推导出可迭代计算的气、液两相逸度系数的表达式:

(11)

式中，ai，bi为气相组分i的范德华常数，a,b由下式计算:

(12)

(13)

用式(14)～(15)计算总体积Vm,t和压缩因子Zm:

(14)

(15)

其中，Tc,i为组分i的临界温度;pc,i为组分i的临界压力。将上述结果代入逸度系数表达式(11)即可求得气相、液相逸度系数。

以温度T为变量，计算180～380 ℃温度范围的KHCl和KH2O，定义氯离子气固相分配系数为KHCl和KH2O的比值，得到HCl-H2O体系条件下氯离子气固相分配系数:

(16)

Cl--Salt体系条件下氯离子气固相分配系数:

(17)

式中，ε为修正系数;φ为修正常数。

以上模型通过Matlab计算得到拟合度最高的HCl-H2O体系条件下氯离子气固相分配系数模型如图2所示，模型表达式为

lgKHCl-H2O=0.716 3e0.002 3T

(18)

2 氯离子气固相分配实验

上述模型(18)为理想状态下得出的，当用于脱硫废水等复杂的Cl--Salt体系时需要通过实验验证并予以修正。引入修正系数和修正常数对HCl-H2O

体系模型进行线性修正，形式如下:

图2 HCl-H2O体系条件下氯离子气固相分配系数模型Fig.2 Model of chloride ions gas-solid phase partition in HCl-H2O system

(19)

式中，CCl-vap为气相挥发的氯浓度；CCl-Total为液相氯的总浓度。

若确定了ε与φ的值或表达式，则实现了对于蒸发过程中氯离子气固相分配系数的量化及预测。

2.1 模型验证及修正参数的确定

氯离子气固相分配系数测定系统如图3所示，每次量取等量的样品溶液置于坩埚中并在恒温干燥箱中加热;空压机产生的空气经加热管加热后引入恒温箱中的密封玻璃罩，恒温箱进出气用玻璃棉保温;实验系统尾部的真空泵将蒸干后的气体抽送至孟氏洗瓶，挥发出的气态氯化物被洗瓶中的NaOH溶液吸收，二级洗瓶内为AgNO3，用以监测一级洗瓶吸收气相Cl-是否完全，将洗瓶置于恒温水浴维持吸收温度恒定，剩余气体排空，整个实验过程通过调节流量计维持空压机进气与真空泵排气处于动态平衡;利用皖仪IC6000离子色谱对洗气瓶内吸收液中溶解的Cl-进行定量分析，实验控制条件见表1。

表1 实验条件参数Table 1 Parameters of experimental condition

用于计算并Cl--Salt体系内氯离子气固相分配系数模型中修正系数，验证模型适用性的样品溶液分别为两个电厂的脱硫废水，部分水质参数见表2。

水样1和水样2中所含Cl-浓度用CTotal表示，按上述实验方法，分别在180，200，230，250，280，300，330，350，380 ℃温度条件下蒸干水样1和水样2。

表2 样品溶液水质参数Table 2 Water quality parameters of sample

图4 lg KCl--Salt与lg KHCl-H2O的线性修正关系Fig.4 Linear correction of lg KCl--Salt and lg KHCl-H2O

2.2 影响因素探究及修正常数的确定

针对影响氯离子气相转化的关键环境参数，如温度，水质参数pH,Cl-浓度,Na+,Ca2+和Mg2+六个因素，设计“六因素四水平”正交实验，分析各个影响因素对氯离子气固相分配的显著性。利用上述图3实验系统得出正交实验结果见表3。

K为均值，R为极差，从R值的大小可以看出各因素对氯离子气固相分配系数影响的显著性顺序:R1(0.781)>R2(0.611)>R3(0.366)>R6(0.211)>R5(0.079)>R4(0.03)，则对其影响最大的因素为温度，其次是pH，Cl-浓度，水中所含阳离子影响最大的是Mg2+，其次是Ca2+，Na+最小，相对温度和pH这两个因素而言，其余均为次要因素。

以正交实验的32组实验数据作为源数据，在Matlab中利用nlinfit多元非线性回归函数和regress多元线性回归函数对3种形式的氯离子气固相分配系数经验模型进行拟合计算。图5为3种不同形式的模型结果，结果表明，拟合度R2(1)=0.815 8，R2(2)=0.808 3，模型3的拟合度R2(3)最高，达到0.925 8，说明模型可信，满足计算要求。

利用regress函数确定模型内参数，得到氯离子气固相分配系数模型:

0.110 1pH-11.35[Cl-]-7.7[Na+]+4.39[Ca2+]+1.94[Mg2+]

(20)

用上述模型(20)的第1项和第2项减去模型(19)的第1项，并分别代入T=200，250，300和350 ℃，得到以温度T为自变量的一次函数，将该项定义为温度校正项，拟合其形式如图6所示。最终的氯离子气固相分配系数模型为

(21)

表3 正交实验结果Table 3 Results of orthogonal experimental

图5 不同形式公式的拟合情况Fig.5 Fitting results of different forms of model

图6 温度校正项拟合Fig.6 Fitting of temperature correction term

[Cl-]-7.7[Na+]+4.39[Ca2+]+1.94[Mg2+]

lgKHCl-H2O=0.716 3e0.002 3T

(22)

表4为计算得出氯离子异相分配系数与其他类似研究测定的结果比较。通过表4的比较可以看出，本实验针对的是氯离子的气固相分配，而其他研究对象为气液相;相比之下，模型适用的最低温度适中，最高温度有所提高，另外，更重要的一点是，其他研究除温度变量外，修正项使用的是液相主体的密度，而本研究的模型中，温度变量是融合在气液相平衡理论之内的，且修正常数包含了影响氯离子气固相分配的关键环境参数及水质参数，更适用于实际工况。