2. 氧传递基础

在好氧生物反应器中,氧气从气相进入细胞内的过程并不是一步完成的,而是一条跨越多个相界面的串联路径。理解这条路径上每一段的阻力大小,是判断“供氧瓶颈在哪里”的前提。本节围绕三个核心问题展开:氧在气—液—固三相中究竟怎么走、双膜理论如何把这条路径简化为可计算的阻力模型、以及体积传质系数 kLa 为什么成为工程上最常用的供氧能力指标。

一、氧从气泡到细胞的传递路径

以典型的好氧发酵或细胞培养为例,空气经分布器形成气泡进入液相,氧分子首先从气泡内部穿过气膜到达气液界面,再穿过液膜进入主体液相;随后氧以溶解氧的形式随液体流动和湍流扩散到达菌体或细胞团表面,穿过细胞团外部的“液膜”和颗粒内部的孔隙扩散,最终穿过细胞膜进入胞内被呼吸链消耗。若存在生物膜或絮体,颗粒内部还会出现明显的扩散梯度,甚至出现氧无法穿透的厌氧核心。

这条路径可以概括为:气相主体 → 气膜 → 气液界面 → 液膜 → 液相主体 → 颗粒外液膜 → 颗粒内扩散 → 细胞膜 → 胞内反应。其中气膜阻力通常远小于液膜阻力,因为氧在水中的溶解度很低(20 ℃、1 atm 下约 8~9 mg/L),氧的传递几乎总是受液膜控制。因此工程上把注意力集中在“气液传质”这一段,而把颗粒内扩散作为第二层需要考虑的问题。

需要强调的是,氧在液相中的“传递”与“混合”是两件事:搅拌或通气可以把溶解氧输运到反应器各处,但真正决定氧能否快速溶入液体的,是气液界面附近的传质速率。这也是为什么单纯加大通气量有时收效有限——气泡停留时间短、界面更新不充分时,增加的空气可能直接逸出。

二、双膜理论:把串联阻力变成可算的式子

双膜理论(Two-Film Theory)由 Whitman 于 1923 年提出,其核心假设是:在气液界面两侧各存在一层滞流膜,膜内传质仅靠分子扩散,膜外主体则因湍流而浓度均匀;界面处气液两相达到瞬时平衡。基于这一假设,氧的传质通量可写为:

N = k_L · (C* − C_L)

其中:
N    —— 氧传递通量,mol/(m²·s) 或 kg/(m²·s)
k_L  —— 液膜传质系数,m/s
C*   —— 与气相氧分压平衡的饱和溶解氧浓度,mg/L
C_L  —— 液相主体实际溶解氧浓度,mg/L
(C* − C_L) —— 传质推动力

把单位体积反应器内的界面面积 a(m²/m³)乘进去,就得到单位体积的氧传递速率 OUR(Oxygen Uptake Rate 的供给侧,常记为 OTR):

OTR = k_L · a · (C* − C_L) = k_La · (C* − C_L)

k_La —— 体积传质系数,单位 1/s 或 1/h
a    —— 比界面面积,m²/m³

双膜理论的价值在于把复杂的界面现象压缩成两个可测参数:kLa 和 C*。C* 由气相氧分压和温度决定,可用亨利定律估算;kLa 则与设备结构、通气量、搅拌强度、介质性质密切相关,通常通过实验测定。需要注意的是,双膜理论是简化模型,对高湍流、表面更新快的体系(如鼓泡塔、气升式反应器)预测精度有限,表面更新理论和渗透理论对此有修正,但 kLa 作为工程参数的形式仍然通用。

三、kLa 的定义与物理意义

kLa 是液膜传质系数 kL 与比界面面积 a 的乘积,量纲为时间的倒数。它的物理意义可以这样理解:在单位推动力(C* − CL = 1 mg/L)下,单位体积液体单位时间内能溶入的氧量。kLa 越大,说明该反应器在相同溶解氧亏缺下供氧能力越强。

kLa 同时受“膜”和“面积”两个因素影响:kL 取决于液相湍动程度和氧的分子扩散系数,搅拌增强、黏度降低都会提高 kL;a 取决于气泡尺寸和持气量,气泡越小、数量越多,比表面积越大。因此在实际操作中,提高转速往往同时改善 kL 和 a,而增大通气量主要增加气泡数量和气含率,但过高气速会导致气泡聚并、停留时间缩短,kLa 反而可能下降。此外,发酵液中的盐、表面活性剂、菌体和高黏度多糖会显著降低 kLa,这也是“同一设备在纯水与真实发酵液中 kLa 差别很大”的原因。

工程上常用动态法测定 kLa:先通氮气将溶解氧降至接近零,再切换空气,记录溶解氧随时间上升的曲线,用下式拟合:

ln(C* − C_L) = ln(C* − C_0) − k_La · t

以 ln(C* − C_L) 对 t 作图,斜率绝对值即为 k_La。
示例(合理假设):C* = 8.0 mg/L,C_0 = 0.5 mg/L,
测得 t = 20 s 时 C_L = 4.2 mg/L,t = 40 s 时 C_L = 6.1 mg/L,
由两点估算 k_La ≈ 0.035 1/s ≈ 126 1/h。

得到 kLa 后,就可以把它与微生物的耗氧速率 OUR 联系起来判断供氧是否充足。稳态下 OTR = OUR,即 kLa·(C* − CL) = OUR。若已知 OUR 和 kLa,可反算出维持代谢所需的 CL;若算出的 CL 低于临界溶解氧浓度(许多好氧菌约 1~2 mg/L,高密度培养时更高),说明当前 kLa 不足,需要提高搅拌或通气,或降低装液量、改用富氧空气。

判断供氧瓶颈的实用思路:先测 kLa 和 C*,再估算 OUR,比较 kLa·C* 与 OUR 的大小。若 kLa·C* 已接近 OUR,说明即使 CL 降到零也无法满足需求,必须提高 kLa 或提高 C*(加压、富氧);若 kLa·C* 明显大于 OUR,则瓶颈可能在混合或局部死区,而非总传质能力。
氧传递路径与双膜理论浓度分布 气—液—固三相传递路径 气相主体 O₂ 分压 pO₂ 气膜 液膜(主阻力) 液相主体 溶解氧 C_L 湍流混合 菌体/絮体 内扩散 颗粒外液膜 细胞膜 胞内呼吸消耗 双膜理论下的氧浓度分布(示意) 距气液界面的位置 → 氧浓度 气相 气膜 液膜 液相主体(C_L 基本均匀) 界面 液膜边界 C*(饱和浓度) C_L(实际浓度) 推动力 C* − C_L OTR = k_L·a·(C* − C_L) k_La 越大,单位体积供氧能力越强

把上述关系落到实际操作上,可以形成一条清晰的判断链:先确认气相氧分压对应的 C*,再通过动态法或稳态法测得当前工况下的 kLa,结合工艺要求的 OUR 计算所需 CL,最后与临界溶解氧比较。若 CL 偏低,优先检查气泡尺寸与分布器状态、搅拌桨型与转速、装液量与挡板配置,再考虑通气量与富氧措施。需要提醒的是,kLa 的测定必须在与实际发酵液性质接近的介质中进行,纯水数据只能作为设备基准,不能直接用于工艺放大计算。

综上,氧传递的本质是一段以液膜阻力为主的串联过程,双膜理论给出了“推动力÷阻力”的定量框架,而 kLa 把界面面积与液膜传质系数合并为一个可测、可比的工程参数。掌握这三者的关系,才能在通气量、搅拌功率与溶解氧控制之间做出有依据的取舍。