4. 单体滴加速率与残余单体峰值关系建模
为什么滴加速率决定残余单体峰值
半间歇丙烯酸聚合中,单体以滴加方式进入反应釜,体系内同时存在两个竞争过程:滴加引入单体的速率 Rin 与 聚合消耗单体的速率 Rp。游离单体浓度 M 的变化由物料衡算决定:
dM/dt = Rin − Rp − (Fout/V)·M
在正常半间歇操作(无馏出、体积近似恒定)下,可简化为 dM/dt ≈ Rin − Rp。当 Rin > Rp 时单体在釜内累积,M 上升;Rin < Rp 时 M 被消耗而下降。因此瞬时游离单体浓度本质上是滴加速率与聚合速率之差的积分,而残余单体峰值出现在累积最严重的时刻——通常是滴加中后期或滴加结束瞬间的“冲顶”。
聚合速率本身又受单体浓度与引发剂浓度控制(Rp = kp·M·(f·kd·[I]/kt)0.5),所以 Rp 会随 M 升高而增大,形成负反馈自限:M 越高,消耗越快。这一自限特性使得稳态游离单体浓度 Mss 可由 Rin = Rp(Mss) 反解得到,是建立滴加曲线模型的核心方程。
从稳态方程反解滴加速率上限
把 Rp 表达式代入稳态条件 Rin = Rp,得到给定目标游离单体浓度 Mtarget 时允许的最大滴加速率:
Rin,max = kp · Mtarget · (f·kd·[I]/kt)0.5
教学示例(假设参数,非实测):取丙烯酸水溶液聚合,kp = 1.5×104 L·mol−1·s−1,kt = 2.0×107 L·mol−1·s−1,过硫酸铵 f·kd 在 80 ℃ 取 3.0×10−4 s−1,[I] = 0.02 mol/L,Mtarget = 0.05 mol/L。则:
sqrt_term = sqrt(3.0e-4 * 0.02 / 2.0e7) = sqrt(3.0e-13) ≈ 5.48e-7
R_in,max = 1.5e4 * 0.05 * 5.48e-7 ≈ 4.1e-4 mol/(L·s)
换算为每小时:4.1e-4 * 3600 ≈ 1.48 mol/(L·h)
若釜体积 5 L、单体总投料 10 mol,则按此上限滴加需 10 / (1.48×5) ≈ 1.35 h。实际滴加时间应留 20%–30% 余量,取 1.7–1.8 h,避免因引发剂诱导期或温度波动导致 Rp 偏低时 M 冲顶。
分段滴加曲线设计
恒定滴加在反应初期(引发剂尚未充分分解、Rp 低)容易造成累积,因此采用前慢后快再收尾的分段曲线更稳妥:
| 阶段 | 时间区间 | 滴加速率(相对 Rin,max) | 目的 |
|---|---|---|---|
| 诱导段 | 0–20% | 0.5× | 等待引发剂分解、避免初期累积 |
| 主滴加段 | 20%–80% | 1.0× | 稳态运行,M 贴近 Mtarget |
| 收尾段 | 80%–100% | 0.6× | 降低 Rin,让 Rp 追平,压平峰值 |
收尾段是关键:滴加结束瞬间 Rin 骤降为 0,若此前 M 已接近限值,则残余单体峰值会出现在滴加末期而非结束之后。将收尾段速率降至 0.6×,可使 M 在滴加结束前已回落至 Mtarget 的 70% 左右,为后续保温消残留出窗口。
模型校核与参数修正
上述模型依赖 kp、kt、f·kd 的取值,实际体系中这些参数受温度、pH、中和度影响较大。工程上不必追求绝对准确,而应通过在线/离线游离单体检测校核模型:在滴加中段与收尾段各取样一次,测定游离单体浓度,与模型预测的 Mss 比较。若实测值持续高于预测值 30% 以上,说明实际 Rp 偏低,应将 Rin,max 下调 20%–30% 或提高引发剂用量;若实测值远低于预测值,则可适当加快滴加以缩短批次时间。