4. 单级萃取分配系数敏感性分析

案例任务与衡算模型

本案例固定单级萃取的操作条件，只改变分配系数 K（萃取相中溶质浓度与萃余相中溶质浓度之比，K = y/x，稀溶液近似为常数），观察两相组成与萃取率如何随之变化。采用的教学示例条件如下：

• 原料液（稀释剂相）流量 F = 100 kg/h，溶质质量分数 x_F = 0.10（即 10 kg/h 溶质）；
• 新鲜溶剂为纯溶剂，流量 S = 60 kg/h，溶质质量分数 y_S = 0；
• 两相不互溶假设成立，F、S 在萃取前后保持不变；
• 分配系数取 K = 0.5、1.0、2.0、3.0、5.0 五个水平。

单级平衡级的总物料衡算与相平衡联立后，可直接解出萃余相组成：

溶质衡算：F·x_F + S·y_S = F·x + S·y
相平衡：  y = K·x
代入：    F·x_F = F·x + S·K·x
解得：    x = F·x_F / (F + S·K)
          y = K·x
萃取率：  E = S·y / (F·x_F) = S·K / (F + S·K)

把 F = 100、S = 60、x_F = 0.10 代入，得到 x = 10/(100 + 60K)，E = 60K/(100 + 60K)。

逐点计算与结果

按上式对五个 K 值逐点计算，结果如下表（数值为按给定输入直接算出的示例结果）：

可以看到：K 从 0.5 增至 5.0（10 倍），萃取率从 23.1% 升到 75.0%，但增幅并非线性——K 由 0.5 到 1.0 时 E 增加 14.4 个百分点，而由 3.0 到 5.0 时只增加 10.7 个百分点，说明高 K 区段继续提高分配系数的收益在递减。

敏感性度量：局部灵敏度与半值点

用差分近似求 dE/dK，可量化“每提高 1 个单位 K 带来多少萃取率增量”。以相邻两点的中心差分估算：

K=0.5→1.0 区间：ΔE/ΔK = (0.375-0.231)/(1.0-0.5) ≈ 0.288 /单位K
K=1.0→2.0 区间：ΔE/ΔK = (0.545-0.375)/(2.0-1.0) ≈ 0.170 /单位K
K=2.0→3.0 区间：ΔE/ΔK = (0.643-0.545)/(3.0-2.0) ≈ 0.098 /单位K
K=3.0→5.0 区间：ΔE/ΔK = (0.750-0.643)/(5.0-3.0) ≈ 0.054 /单位K

灵敏度随 K 增大而单调下降。另一个有用的判据是“半值点”：令 E = 0.5，解得 K_1/2 = F/S = 100/60 ≈ 1.67。即当 K 低于约 1.67 时，单级萃取连一半溶质都取不出来；K 远大于 1.67 后，E 才逐渐逼近上限 S/(S+F/K) → 1，但受溶剂比限制，实际趋近速度很慢。

结果判读与工程结论

把上述计算整理成可复用的判读要点：

1. 方向性：K 增大，萃余相 x 单调下降、萃取相 y 单调上升、萃取率 E 单调上升，三者方向一致，无交叉。
2. 敏感区划分：以半值点 K_1/2 = F/S 为界，低于该值时 dE/dK 大（本例约 0.29/单位K），高于 3 倍 F/S 后 dE/dK 降到 0.05 量级。
3. 上限约束：即使 K → ∞，单级萃取率也只能趋近 100%，但实际受溶剂比限制，K = 5 时 E 仅 75%，说明单靠分配系数无法突破溶剂比瓶颈。
4. 校核方法：任取一点回代总溶质衡算，如 K = 2 时 100×0.10 = 100×0.0455 + 60×0.0909 = 4.55 + 5.45 = 10.00 kg/h，闭合，说明计算无误。