2. 传质基础与kLa定义
在发酵工程中,好氧发酵的成败往往不取决于菌种本身,而取决于氧气能否及时、足量地从气相进入液相,再被细胞利用。氧气在水中的溶解度极低,20℃、1 atm下纯水平衡溶解氧浓度仅约9 mg/L,而微生物的耗氧速率常常远高于自然溶解速率。因此,发酵罐必须依靠强制通风与搅拌,持续把气相中的氧传递到液相。描述这一传递过程快慢的核心参数,就是体积传质系数 kLa。
双膜理论:把复杂界面拆成两层膜
气液传质最经典的解释模型是双膜理论(Two-Film Theory)。它假设气液界面两侧各存在一层滞流膜:气膜和液膜,界面本身处于平衡状态,传质阻力全部集中在这两层膜上。氧气从气相主体出发,先穿过气膜到达界面,再穿过液膜进入液相主体。由于氧气难溶于水,液膜阻力通常远大于气膜阻力,因此好氧发酵中传质速率主要受液膜控制。
按照双膜理论,传质通量可写成:
N = k_L · (C* - C_L)
N :单位界面面积的传质通量,mol/(m²·s)
k_L :液膜传质系数,m/s,反映单层液膜的传递能力
C* :与气相氧分压平衡的液相饱和溶解氧浓度,mol/m³
C_L :液相主体实际溶解氧浓度,mol/m³
(C* - C_L):传质推动力,浓度差越大传质越快
这里 kL 只描述“单位面积液膜传得多快”,单位是 m/s。但发酵罐里真正决定总传氧量的,不只是单位面积传得多快,还有“有多少面积可供传递”。气泡越小、气含率越高、搅拌越强,单位液体体积所对应的气液界面面积 a 就越大,a 的单位是 m²/m³(比表面积)。把两者相乘,就得到体积传质系数 kLa,单位是 1/s(或 s⁻¹)。
kLa 的物理意义与单位
kLa 的物理意义是:单位体积发酵液、单位浓度差推动力下,每秒钟传递的氧气量。它把“膜传递能力”和“界面面积供给能力”合并成一个可测量的宏观参数。工程上常用的传质速率方程为:
OTR = k_L a · (C* - C_L)
OTR :体积传氧速率,mol/(m³·s)
k_L a:体积传质系数,s⁻¹
C* :饱和溶解氧浓度,mol/m³
C_L :实际溶解氧浓度,mol/m³
从单位看,kLa 是 s⁻¹,乘以浓度差(mol/m³)后得到 mol/(m³·s),正好是单位体积的传氧速率。这也解释了为什么 kLa 常被直接当作发酵罐“供氧能力”的标尺:在同样的 C* 和 CL 下,kLa 越大,OTR 越大,罐子能支撑的菌体浓度和耗氧速率就越高。
kL 与 a 的独立贡献
理解 kL 与 a 的独立性,是分析传质瓶颈的关键。kL 主要受液相物性(黏度、扩散系数、表面张力)和气泡周围液膜厚度影响;a 则主要受气泡尺寸分布、气含率、通气量和搅拌功率影响。两者可以独立变化:
- 只提高 kL:例如降低发酵液黏度、改善液相湍动,使液膜变薄,但气泡数量和比表面积未必增加。
- 只提高 a:例如增大通气量或提高搅拌转速,把大气泡打碎成小气泡、提高气含率,界面面积显著增加,但单层液膜传递能力未必同步提升。
- 两者同时提高:高黏度体系加消泡剂、优化搅拌桨组合,既细化气泡又强化液膜更新,kLa 提升最明显。
实际发酵中,高黏度、高菌浓体系往往出现“a 上不去、kL 也下降”的双重困境:气泡易聚并变大,气含率低,同时液膜变厚。此时单纯加大通气量可能只增加大气泡和泡沫,kLa 提升有限,甚至因泡沫带走菌体而恶化。因此辨析 kL 与 a 的贡献,能帮助判断该调搅拌、调通气、改配方还是改设备。
影响 kLa 的主要因素
从工程操作角度看,kLa 受设备、操作和物性三类因素共同影响。设备方面,搅拌桨型式与直径、挡板、分布器孔径决定气泡破碎与分散;操作方面,通气量、搅拌转速、罐压通过改变气含率和气泡尺寸影响 a,也通过湍动强度影响 kL;物性方面,发酵液黏度、表面张力、离子强度和消泡剂会改变气泡聚并行为和液膜厚度。经验上,kLa 常随搅拌功率和表观气速增大而上升,但存在边际递减,过高转速还会剪切损伤菌体、增加能耗。
需要强调的是,kLa 是“设备+操作+物料”共同决定的宏观参数,不能脱离具体发酵体系单独比较。同一台罐在清水和在高黏度发酵液中的 kLa 可能相差数倍,因此评估时应在实际发酵液或接近实际物性的条件下测定。
kLa 的测定思路(示例)
工业上最常用的是动态法(动态排气法)。以某好氧发酵罐为例(以下为教学示例参数,非实测):先通空气使溶解氧接近饱和 C*,然后停止通气并保持搅拌,菌体耗氧使 CL 按耗氧速率线性下降;随后恢复通气,记录 CL 随时间回升的曲线。在回升段,若忽略耗氧,传质方程为 dCL/dt = kLa(C*−CL),积分得 ln[(C*−CL0)/(C*−CL)] = kLa·t。以 ln(C*−CL) 对 t 作图,斜率绝对值即为 kLa。
示例:C* = 7.5 mg/L,恢复通气后测得
t=0 s, C_L=2.0 mg/L → C*−C_L=5.5
t=10 s, C_L=3.6 mg/L → C*−C_L=3.9
t=20 s, C_L=4.8 mg/L → C*−C_L=2.7
t=30 s, C_L=5.7 mg/L → C*−C_L=1.8
ln(5.5)=1.705, ln(3.9)=1.361, ln(2.7)=0.993, ln(1.8)=0.588
斜率 ≈ (0.588−1.705)/30 ≈ −0.0372 s⁻¹
故 k_L a ≈ 0.037 s⁻¹(约 2.2 min⁻¹)
实际测定时要注意:菌体耗氧会压低回升曲线,需用亚硫酸盐法或灭活菌体做对照校正;探头响应时间应远小于 1/kLa;C* 应按实际温度、罐压和培养基盐析效应修正。测得 kLa 后,结合菌体的临界耗氧速率,就能判断当前操作条件下供氧是否足够,并为放大和工艺调整提供依据。
小结:双膜理论给出了“气膜—界面—液膜”的传质图像,kL 描述液膜单位面积的传递能力,a 描述单位体积的界面供给能力,二者相乘得到 kLa(s⁻¹),直接决定体积传氧速率 OTR。掌握三者的量纲关系与独立影响因素,才能在发酵供氧不足时准确判断瓶颈在“膜”还是在“面积”,从而选择正确的调整手段。