4. 化工因子的识别与水平设定
全因子设计的第一步不是排表,而是把“工艺语言”翻译成“因子语言”。在化工过程中,温度、压力、浓度、流量、催化剂用量这些量几乎总在同时变化,但并非每一个都值得进入试验设计。因子识别要回答两个问题:哪些变量对目标响应有实质影响、且能在试验中独立调节;水平设定要回答另一个问题:每个因子取哪几个数值,才能让响应变化既被激发出来、又不至于跑出安全或工艺可行边界。这两步做偏,后续的效应估计、方差分析和交互作用解释都会失去意义。
一、从工艺变量到可控因子
识别因子时,建议按“影响机理—可调性—可测性—独立性”四条线索逐一筛查。影响机理指该变量是否通过反应速率、平衡、传质或传热直接作用于响应;可调性指在现有装置上能否在试验期间稳定设定并保持;可测性指设定值和实际值能否被准确记录;独立性指它是否与其它因子强绑定,例如提高温度往往同时改变副反应选择性,若两者无法解耦,就应合并为一个复合因子或固定其一。
以某液相酯化反应为例,候选变量包括反应温度、系统压力、醇酸摩尔比(浓度类因子)、进料流量、催化剂用量和搅拌转速。经筛查:温度直接影响速率常数,可调且可测,保留;压力主要影响低沸点组分的回流与气相组成,在常压附近变化不显著,暂固定;摩尔比决定平衡转化率,保留;进料流量改变停留时间,与反应器负荷耦合,若装置能稳定控温则可作为因子;催化剂用量影响速率但不改变平衡,保留;搅拌转速在湍流区已无明显传质限制,固定。最终得到四个因子:温度、摩尔比、流量、催化剂用量。这一筛选过程本身就是工程判断,不能仅凭“变量多就都放进去”。
二、水平设定的工程依据
水平设定的核心是“范围要覆盖真实响应曲率,又不能越界”。常用做法是先确定安全操作窗口和工艺可行窗口,再在窗口内取低、高两个水平。低水平不宜贴着设备下限,高水平不宜贴着安全上限,两端都应留出约10%~20%的余量,防止试验波动导致越界。对于已知存在最优点的因子(如温度过高会引发副反应或催化剂失活),水平应围绕当前操作点对称展开,而不是单侧拉伸。
仍以上述酯化反应为例,当前操作温度为85 ℃,工艺允许60~110 ℃,催化剂用量为醇质量的0.5%,可行范围0.2%~1.2%,摩尔比为1.2,可行1.0~2.0,流量为2.0 L/h,可行1.2~2.8 L/h。若直接取可行范围端点,温度跨度达50 ℃,可能使低水平反应过慢、高水平副产物激增,响应出现强烈非线性甚至越界。更稳妥的做法是围绕当前点取±10~15 ℃,即75 ℃与95 ℃;催化剂取0.3%与0.8%;摩尔比取1.1与1.5;流量取1.6与2.4 L/h。这样每个因子的水平跨度既能激发出可观测的响应差异,又落在已验证的稳定操作区内。
三、水平范围过窄与过宽的典型失效
水平范围过窄,是化工试验设计中最常见的隐性错误。若温度只取83 ℃与87 ℃,响应差异可能被测量噪声完全淹没,效应估计的置信区间变宽,主效应显著性检验失效,交互作用更无从识别。判断依据是:水平跨度应使响应变化至少达到测量重复性标准差的3~5倍。若预试验显示该跨度下转化率变化不足1个百分点,而重复测量波动就有0.8个百分点,就应扩大水平或改用更精确的分析方法。
水平范围过宽同样危险。温度从60 ℃直接跳到110 ℃,低水平可能反应未启动、高水平可能焦化结垢,两端响应都偏离线性区,效应估计被曲率污染,甚至出现“高水平反而变差”的假象。更隐蔽的问题是越界后工艺机理改变,例如高压下气相组分冷凝、催化剂浓度超过溶解度上限,此时因子已不再是同一物理量,全因子模型的前提被破坏。因此水平设定后应做一次边界校核:两端水平是否都在已验证的稳定操作记录内、是否有安全联锁或报警会触发、取样与分析是否仍具代表性。
四、水平编码与试验前确认
确定实际水平后,应尽早转换为编码水平(低=-1,高=+1),便于后续效应计算与模型拟合。编码不改变试验内容,但能避免量纲差异导致的系数不可比。转换关系为:编码值 =(实际值 − 中心值)/(半跨度)。以上述温度为例,中心值85 ℃、半跨度10 ℃,则75 ℃对应−1,95 ℃对应+1。
因子 低水平(-1) 中心(0) 高水平(+1) 半跨度
反应温度/℃ 75 85 95 10
催化剂用量/% 0.3 0.55 0.8 0.25
醇酸摩尔比 1.1 1.3 1.5 0.2
进料流量/(L/h) 1.6 2.0 2.4 0.4
正式试验前还应完成三项确认:其一,每个水平点能否在装置上稳定保持至少一个取样周期,温度波动建议控制在±1 ℃、流量±3%以内;其二,分析方法的重复性与水平跨度匹配,必要时先做重复测定估计噪声;其三,记录实际值而非仅记录设定值,因为化工装置的设定值与实测值常有偏差,后续分析应以实测值为准。完成这些确认后,因子与水平才真正具备进入全因子表的资格,后续的主效应与交互作用分析才有可靠的比较基础。