3. 传质基础与界面理论

在生物反应器中,底物从液相主体进入细胞、氧气从气泡进入发酵液、代谢产物从细胞释放到液相,本质上都是物质在相内或相间的迁移过程。传质速率往往决定反应器的产能上限,因此理解传质的基本机制与界面理论,是分析溶氧、补料、产物提取等问题的共同基础。

分子扩散:传质的“底层通道”

分子扩散由浓度梯度驱动,遵循 Fick 第一定律:

J_A = -D_AB · (dC_A/dz)

其中 J_A 为组分 A 的摩尔通量(mol·m⁻²·s⁻¹),D_AB 为 A 在 B 中的分子扩散系数(m²·s⁻¹),dC_A/dz 为浓度梯度。负号表示扩散方向与浓度升高方向相反,即物质自发从高浓度向低浓度迁移。扩散系数是物性参数:小分子在低黏度液体中约 10⁻⁹ m²·s⁻¹ 量级,在细胞团或生物膜内因孔隙曲折可降至 10⁻¹⁰~10⁻¹¹ m²·s⁻¹。扩散是分子尺度的随机运动结果,不需要宏观流动,但速率随距离增大迅速衰减——扩散特征时间约为 L²/D,若 L 为 1 mm、D 为 10⁻⁹ m²·s⁻¹,则特征时间约 1000 s,说明单靠扩散难以支撑毫米级以上尺度的高耗氧需求。

对流传质:流动带来的“加速器”

对流传质是流体宏观运动携带物质迁移的过程,通量可写为:

N_A = k_L · (C_A,bulk - C_A,s)

式中 k_L 为液膜传质系数(m·s⁻¹),C_A,bulk 与 C_A,s 分别为液相主体和界面处的浓度。对流传质把复杂的流动—扩散耦合过程压缩为一个系数,便于工程计算。k_L 受流速、黏度、扩散系数和界面几何影响,通常通过无量纲关联式估算,例如:

Sh = k_L · L / D = f(Re, Sc)

Sh 为 Sherwood 数,Re 为 Reynolds 数,Sc 为 Schmidt 数。在搅拌槽中,提高转速或通气量会增大 Re,从而提高 k_L;但 k_L 的提升并非无限,因为界面附近仍存在一层难以完全消除的扩散阻力层。

界面理论:从双膜到表面更新

气液或液固界面附近的传质阻力如何分布,是界面理论要回答的核心问题。双膜理论假设界面两侧各存在一层静止的滞流膜,膜内为分子扩散,膜外为充分混合的主体,总传质阻力等于气膜阻力与液膜阻力之和:

1/K_L = 1/k_L + 1/(H · k_G)

其中 K_L 为以液相为基准的总传质系数,H 为 Henry 常数,k_G 为气膜传质系数。对于难溶气体如氧气,H 很小,1/(H·k_G) 项往往可忽略,总阻力主要由液膜控制,这也是好氧发酵中强调提高 k_L a 的原因。

双膜理论把界面视为静止,但实际气液界面不断被液流冲刷和更新。渗透理论认为流体微元在界面停留时间 t 内按非稳态扩散渗透,平均传质系数与 √(D/t) 成正比;表面更新理论进一步假设界面微元按一定年龄分布随机更新,引入表面更新频率 s,得到:

k_L = √(D · s)

该式表明,k_L 由扩散系数和表面更新频率共同决定。搅拌、气泡破裂、液滴聚并都会提高 s,从而强化传质。三种理论并非互相排斥,而是从不同尺度描述同一界面过程:双膜理论适合工程估算,渗透与表面更新理论更适合解释湍流界面的动态行为。

气液界面传质阻力分布与三种界面理论对比 气相主体(充分混合) C_G,bulk ≈ 均匀,氧分压 pO₂ 气膜阻力:1/(H·k_G),难溶气体可忽略 气膜(滞流层)——分子扩散 气液界面(C* 由 Henry 定律决定) 液膜(滞流层)——分子扩散,主要阻力 液相主体(湍流混合) C_L,bulk 为溶解氧浓度,受细胞耗氧影响 液膜阻力:1/k_L,好氧发酵中通常为控制步骤 O₂ 传质方向 O₂ 被细胞消耗 双膜理论 界面两侧静止膜 1/K_L = 1/k_L + 1/(H·k_G) 适合工程估算 渗透理论 微元停留时间 t 内非稳态扩散 k_L ∝ √(D/t) 解释短时接触传质 表面更新理论 界面微元按年龄分布随机更新 k_L = √(D·s) s 为表面更新频率

传质系数与传质阻力:把复杂过程压缩成可测参数

工程上很少直接计算界面附近的浓度分布,而是用总传质系数 K_L 或体积传质系数 k_L a 描述整体能力。k_L a 的单位是 s⁻¹ 或 h⁻¹,其中 a 为单位体积液体的气液界面面积(m²·m⁻³)。在好氧发酵中,氧传递速率 OTR 常写为:

OTR = k_L a · (C* - C_L)

C* 为与气相氧分压平衡的液相饱和浓度,C_L 为实际溶解氧浓度。若某发酵体系在 30 ℃、常压空气下 C* ≈ 7.5 mg·L⁻¹,实测 C_L = 2.0 mg·L⁻¹,k_L a = 120 h⁻¹,则:

OTR = 120 × (7.5 - 2.0) = 660 mg·L⁻¹·h⁻¹

若细胞耗氧速率 OUR 为 600 mg·L⁻¹·h⁻¹,则 OTR > OUR,溶解氧可维持;若提高细胞密度使 OUR 升至 750 mg·L⁻¹·h⁻¹,则 OTR 不足,C_L 会持续下降,直至新的平衡点或细胞缺氧。此时可通过提高转速、增大通气量或富氧通气提高 k_L a,也可降低 OUR 或提高 C*。

传质阻力具有串联性质:气膜阻力、界面阻力、液膜阻力、细胞团内扩散阻力依次叠加。对氧气这类难溶气体,液膜阻力通常占主导;对易溶气体或挥发性产物,气膜阻力可能不可忽略;对固定化细胞或颗粒污泥,内扩散阻力可能成为新的控制步骤。判断控制步骤的实用方法是比较各段阻力的数量级,或通过改变操作条件观察总传质系数是否随之显著变化。

在实验测定 k_L a 时,常用动态法:先通氮气将 C_L 降至接近零,再切换空气并记录 C_L 随时间上升曲线,由 ln(C* - C_L) 对时间作图的斜率得到 k_L a。该方法假设液相充分混合、气相氧分压恒定,且忽略探头响应滞后。

从理论到操作:如何判断与改善传质

面对一个溶氧不足的反应器,可按以下顺序排查:先确认 C* 是否因温度、压力或气体组成而偏低;再检查 k_L a 是否因转速、通气量、气泡尺寸或挡板配置而不足;然后判断 OUR 是否因补料过快、细胞浓度过高而超出供氧能力;最后考虑是否存在局部混合死区或细胞团内扩散限制。若提高转速后 k_L a 明显上升,说明液膜或界面更新是控制因素;若提高通气量但 k_L a 变化不大,则可能受气泡聚并或气膜阻力限制。

三种界面理论给出了不同的改善方向:双膜理论提示减小液膜厚度,渗透理论提示缩短接触时间或增大扩散系数,表面更新理论提示提高表面更新频率。实际反应器设计中,这些方向往往同时起作用——搅拌既减薄液膜,又提高表面更新频率,还打碎气泡增大 a,因此 k_L a 对操作条件通常非常敏感。理解传质系数与阻力的串联关系,才能在溶氧、补料和产物提取之间做出有依据的权衡。