3. 杂质累积的物料衡算基础

母液循环是提高结晶收率的常规手段,但它同时把杂质留在了系统内部。要判断一套循环工艺能否长期稳定运行,不能只看单批结晶效果,而必须回答三个定量问题:母液被循环放大了多少倍、杂质相对产品被富集到什么程度、在给定排放策略下杂质浓度最终会稳定在哪个水平。本节围绕这三个问题建立物料衡算框架,并给出可直接套用的计算示例。

一、先约定衡算边界与参照基准

做杂质衡算前必须先固定系统边界,否则循环比与富集因子会互相矛盾。工程上通常取“结晶器+母液储槽+循环泵”为控制体,以单位时间(或单位批次)为基准,以进料为输入、产品与排放(含母液外排、洗水、滤饼夹带)为输出。全文约定:杂质浓度以质量分数计,循环比与富集因子均为无量纲比值,累积速率以“每批进入系统的杂质质量”为正方向。若改用体积浓度或摩尔分数,需在全文与图示中同步替换,不可混用。

设进料量为 F(kg/批),进料中杂质质量分数为 xF;产品量为 P,产品中杂质质量分数为 xP;母液循环量为 R,母液中杂质质量分数为 xR;外排母液量为 B,外排液中杂质质量分数同样为 xR(同一母液池)。这些符号是后续全部计算的基础。

二、母液循环比:把“回流量”折算成放大倍数

循环比最直观的定义是回流量与产品量之比:

循环比 r = R / P

它回答的是“每产出 1 kg 产品,有多少 kg 母液被送回结晶器”。例如某批产品 P = 1000 kg,母液回流 R = 2500 kg,则 r = 2.5。若以进料为基准,也可定义 rF = R / F,两者通过收率换算,工程上更常用 r = R/P,因为它直接对应“产品被母液反复洗涤”的次数概念。

需要注意,循环比本身并不等于杂质放大倍数。循环比只描述物流量,杂质是否放大还取决于杂质在产品与母液之间的分配。把两者混为一谈,是现场估算偏差的常见来源。

三、杂质富集因子:分配系数决定富集强度

定义杂质在母液与产品中的浓度比为富集因子:

富集因子 β = x_R / x_P

β 的物理含义是“杂质更愿意留在母液里还是进入晶体”。β 越大,说明杂质越倾向于留在液相,循环时越容易被“攒”起来。β 由结晶热力学与操作条件共同决定,可通过小试测定:在同一批母液与晶体上分别取样分析 x_R 与 x_P,直接求比值。若杂质几乎不进入晶格,β 可达数十;若杂质共沉淀明显,β 可能接近 1,此时循环对杂质浓度的影响反而较小,但产品纯度本身已受威胁。

把循环比与富集因子结合,可得到单批杂质进入产品与滞留母液的分配图像:每批随进料带入的杂质,一部分随产品带走,其余留在母液中参与下一轮循环。这正是累积的根源。

母液循环系统的杂质物料衡算边界 进料 F x_F(杂质) 结晶器 固液分离 母液池 x_R 晶体 x_P 产品 P x_P(杂质) 外排母液 B x_R(杂质) 母液回流 R 循环比 r = R / P 富集因子 β = x_R / x_P β 越大,杂质越易滞留母液

四、累积速率:每批净增多少杂质

在稳态建立之前,系统内杂质总量会逐批变化。以单批为时间步,设第 n 批母液中杂质质量为 Mn,则下一批的净增量为“进料带入 − 产品带走 − 外排带走”:

ΔM = F·x_F − P·x_P − B·x_R

由于 x_P = x_R / β,可把上式改写为只含 x_R 的形式:

ΔM = F·x_F − (P/β + B)·x_R

这就是累积速率的核心表达式。括号内 (P/β + B) 可理解为系统对杂质的“排出能力”:产品带走的等效项 P/β 与外排项 B 之和。排出能力越大,累积越慢、稳态越低。举例(示例假设值):F = 1200 kg/批,x_F = 0.8%,P = 1000 kg/批,β = 8,B = 150 kg/批。若当前 x_R = 3%,则 ΔM = 1200×0.008 − (1000/8 + 150)×0.03 = 9.6 − 12.0 = −2.4 kg/批,说明此时系统杂质在下降,尚未达到稳态;若 x_R 降到 2%,则 ΔM = 9.6 − 6.0 = +3.6 kg/批,杂质在上升。可见存在一个使 ΔM = 0 的临界浓度。

五、稳态浓度预测:令累积速率为零

稳态即 ΔM = 0,解得母液杂质稳态浓度:

x_R* = F·x_F / (P/β + B)

代入上例:x_R* = 9.6 / (125 + 150) = 9.6 / 275 ≈ 3.49%。也就是说,在 β = 8、外排 150 kg/批 的条件下,母液杂质会逐步趋近约 3.5% 后稳定。若把外排量降到 80 kg/批,则 x_R* = 9.6 / (125 + 80) ≈ 4.68%,稳态明显抬高;若外排提高到 250 kg/批,x_R* ≈ 2.56%。这组对比清楚说明:外排量是压低稳态浓度最直接的调节手段,而 β 由体系性质决定,靠操作调整的空间有限。

实际预测时还需注意两点。其一,上式假设 β 与各物流量在循环过程中保持不变,若杂质累积改变了晶体生长环境,β 会漂移,需分段核算或定期用小试修正。其二,达到稳态需要一定批数,可由 ΔM 的递推关系估算趋近速度:每批 x_R 向 x_R* 逼近的比例约为 (P/β + B)/(系统母液持液量),持液量越大,趋近越慢。工程上常据此判断“开车后多少批需要加强外排”。

实用判断:先算 x_R*,再与产品杂质上限反推允许的 x_P = x_R*/β。若 x_P 已超产品指标,说明单靠外排不够,必须从源头降低 x_F 或改变结晶条件提高 β,而不是继续加大循环比。

六、把三个量串起来用

循环比 r 决定母液被重复利用的强度,富集因子 β 决定杂质滞留母液的倾向,二者共同决定累积速率 ΔM 与稳态浓度 x_R*。现场可按“取样测 β → 统计 F、P、B 与 x_F → 算 x_R* → 对比产品指标 → 调整外排或进料”的顺序闭环操作。需要强调的是,上述数值均为教学示例假设,实际装置应以本厂化验数据替换后重新计算,并保留每批的 x_R 记录以校核模型偏差。只有把衡算结果与在线或离线分析对照,稳态预测才具备指导排放策略的可信度。