3. 生物反应器放大中的混合与传质数据解析

放大失败往往不是“菌种不行”，而是小罐里被掩盖的传递瓶颈在大罐里暴露出来。小试反应器通常混合快、供氧足、剪切均匀，而工业罐体积增大后，单位体积功率输入下降、气泡停留时间变长、局部死区与剪切峰值同时出现。因此，放大的核心工作是把混合、传质、剪切三类参数从“现象”变成“可比较的数据”，再判断哪些参数应保持恒定、哪些只能折中。

一、三类参数各自描述什么

体积传质系数 kLa 描述氧从气相进入液相的速率，单位 s⁻¹ 或 h⁻¹。它由 kL（液膜传质系数）与 a（比表面积）组成，工程上更关注 kLa，因为通气量与搅拌主要改变的是气泡比表面积。kLa 直接决定临界溶氧能否维持，是需氧发酵放大的第一约束。

混合时间 t_m 通常定义为示踪剂加入后，某监测点浓度达到最终稳定值 ±5% 所需的时间。它反映的是宏观混合能力，决定补料、酸碱、消泡剂加入后是否形成局部过浓区。混合时间过长会导致局部底物抑制、pH 波动和副产物积累。

剪切力 对丝状菌、动物细胞和剪切敏感菌株尤为关键。工程上常用叶端线速度 v_tip = πND 与平均剪切速率 γ̇ ≈ N（N 为转速）作为可比指标，桨叶区局部剪切峰值远高于平均值，因此放大时更应关注桨叶附近的能量耗散率 ε。

二、kLa 的测量与计算

实验室最常用动态法：先通氮气或停止通气，让溶氧降至接近零，再恢复通气并记录 DO 随时间上升曲线。对液相氧衡算积分可得：

ln(C* - C_L) = -kLa · t + 常数
以 ln(C* - C_L) 对 t 作直线，斜率绝对值即为 kLa

其中 C* 是与气相平衡的饱和溶氧浓度，C_L 是实测液相溶氧。若用空气，C* 取该温度、压力、盐度下的饱和值；若用富氧，需按分压重新计算 C*。示例：某 5 L 罐在 30 ℃、1 vvm、300 rpm 下，DO 从 5% 升至 80% 饱和用时 42 s，取 C* = 100%，则 ln(0.95/0.20) ≈ 1.56，kLa ≈ 1.56/42 ≈ 0.037 s⁻¹ ≈ 134 h⁻¹。该数值属于教学示例，实际应以多点拟合与重复实验为准。

kLa 与操作变量的经验关系常写成：

kLa = C · (P_g/V)^α · v_s^β
P_g/V：单位体积通气功率（W/m³）
v_s：表观气速（m/s）
α、β：由实验拟合，常在 0.4～0.8 与 0.4～0.7 之间

放大时若保持 P_g/V 与 v_s 同时恒定，kLa 通常可近似保持；但大罐中静压使底部 C* 高于顶部，实际平均推动力需按分段或对数平均处理。

三、混合时间的测量与规模相关性

混合时间常用电导法或pH 示踪法：在液面附近瞬时加入 NaCl 或酸碱，在远离加料点的位置记录电导/pH 变化，取达到最终值 ±5% 的时间。测量点应避开加料口与桨叶直接冲刷区，否则会低估混合时间。

混合时间与单位体积功率、罐径的关系可近似表示为：

t_m ∝ (V)^n ，层流/过渡区 n ≈ 0.3～0.5
t_m ∝ D^(2/3) · (P/V)^(-1/3) （湍流区常用形式）

这意味着体积放大 1000 倍时，混合时间可能增加数倍到十余倍。若小罐 t_m = 5 s，放大后可能达到 30～60 s，补料脉冲就必须改为连续或高频小量补加，否则局部浓度峰值会触发副反应。

四、剪切力的估算与判据

桨叶区最大能量耗散率可近似用 ε_max ≈ K · N³D² 估算，K 与桨型有关。叶端线速度 v_tip = πND 是更直观的判据：动物细胞通常要求 v_tip < 1.5～2.5 m/s，丝状菌可耐受更高但需关注菌丝断裂。放大时若保持 v_tip 恒定，则 N 随 D 增大而下降，但 P/V 会显著降低，kLa 与混合时间随之恶化。因此实际放大常采用“kLa 优先、剪切设上限”的折中策略。

五、规模相关性分析与放大决策

把上述参数整理成规模相关性表，是放大决策的直接依据。下表为教学示例，数值为合理假设，用于说明比较方法，不代表实测结果。

从表中可读出关键矛盾：若按恒定 v_tip 放大，生产罐转速会降到很低，P/V 与 kLa 明显不足；若按恒定 P/V 放大，转速与 v_tip 又会超出剪切上限。因此实际放大通常先确定剪切上限，再在该上限内尽量提高 P/V 与 v_s，使 kLa 满足临界溶氧要求；混合时间的恶化则通过连续补料、多点加料、提高循环次数来补偿，而不是单纯依赖提高转速。

最后需要强调，kLa、混合时间与剪切并非独立：提高转速同时改善 kLa 与 t_m，却恶化剪切；增大通气改善 kLa 与混合，却可能带来泡沫、气泡聚并和氧分压变化。放大数据解析的价值，正在于把这些耦合关系量化，找到满足供氧与混合、同时不越过剪切上限的操作窗口。