4. 氧传递理论与kLa
在好氧生物反应器中,氧气从气相进入液相、再被细胞消耗,整个过程由“气液传质”和“细胞呼吸”两个串联环节控制。当溶氧浓度低于临界值时,细胞代谢立即受限,因此工程上真正关心的不是单纯的氧溶解度,而是单位时间、单位体积内液相能向细胞供应多少氧。体积传质系数 kLa 正是描述这一供应能力的核心参数,它把“传质推动力”和“设备供氧能力”统一到一个可测、可比的指标上。
一、kLa 的物理意义:把“浓度差”换算成“供氧速率”
气液传质常用双膜理论描述:气相中的氧先溶解到气液界面,再通过液膜扩散进入主体液相。由于氧在水中的溶解度很低,液膜阻力通常占主导,因此氧传递速率 OTR(Oxygen Transfer Rate)可写为:
OTR = kL · a · (C* − CL)
= kLa · (C* − CL)
其中:
kL —— 液膜传质系数(m/h),反映单气泡界面的传质快慢
a —— 单位液体体积的气液界面面积(m²/m³)
kLa —— 体积传质系数(1/h),kL 与 a 的乘积
C* —— 与气相氧分压平衡的液相饱和溶氧浓度(mg/L)
CL —— 液相实际溶氧浓度(mg/L)
kLa 的单位是 h⁻¹,物理含义是单位浓度差驱动下、单位体积液体每小时能传递的氧量。它把 kL 和 a 合并成一个整体参数,因为在实际反应器中,气泡大小、数量、分布和液膜厚度很难单独测量,而 kLa 可以通过溶氧电极直接测定。对细胞而言,供氧是否足够取决于 OTR 与耗氧速率 OUR(Oxygen Uptake Rate)的平衡:
dCL/dt = OTR − OUR = kLa · (C* − CL) − qO₂ · X
qO₂ —— 比耗氧速率(mg O₂/(g 细胞·h))
X —— 细胞浓度(g/L)
稳态时 dCL/dt = 0,可得 CL = C* − OUR / kLa。这说明:kLa 越大,达到相同 OUR 时所需维持的溶氧浓度越高,或者说在相同溶氧下能支撑的细胞密度越大。放大中细胞浓度往往提高,OUR 随之上升,若 kLa 不能同步提高,CL 就会跌到临界值以下,导致供氧限制。
二、kLa 的测定:动态法与稳态法
kLa 不能靠公式直接算出,必须通过实验测定。工业上最常用的是动态法和稳态法,两者适用场景不同。
动态法(Dynamic Method)适用于间歇或补料分批培养。操作上先停止供气或通入氮气,让溶氧降到较低水平,然后恢复通气,记录 CL 随时间上升的曲线。在细胞耗氧可忽略或已知的条件下,由物料衡算可得:
dCL/dt = kLa · (C* − CL) − OUR
若先测 OUR(停止通气后 CL 线性下降的斜率),再恢复通气,
则 ln[(C* − CL₀)/(C* − CL)] = kLa · t
以 ln[(C* − CL)] 对 t 作图,斜率即为 −kLa。
动态法的关键是 OUR 必须单独测准,且 C* 要按实际温度、压力、气相组成修正。若细胞耗氧很强,恢复通气后 CL 上升缓慢,曲线弯曲,此时必须扣除 OUR 才能得到真实 kLa。
稳态法(Steady-State Method)适用于连续培养或耗氧稳定的体系。在稳态下 dCL/dt = 0,于是 kLa = OUR / (C* − CL)。只要测得 OUR、C* 和 CL,就能直接算出 kLa。稳态法不需要切换气体,对细胞干扰小,但它依赖 OUR 的准确测定,且只能给出当前操作条件下的 kLa,不能像动态法那样一次实验得到整条响应曲线。
两种方法各有取舍:动态法信息量大、适合间歇过程;稳态法操作简单、适合连续或高耗氧体系。实际工作中常用动态法做条件筛选,再用稳态法在目标工况下复核。
三、影响 kLa 的操作与物性因素
kLa 是 kL 与 a 的乘积,因此任何改变液膜厚度或气液界面面积的因素都会影响它。操作层面,通气量 Q 提高会增加气泡数量,a 上升,kLa 通常近似随 Q 的 0.4~0.7 次方增长;搅拌转速 N 通过剪切力把大气泡打碎成小气泡,同时减薄液膜,kLa 随 N 的 0.6~1.0 次方增长,但过高转速会损伤剪切敏感细胞,也会增加能耗。罐压提高会增大氧分压,使 C* 上升,但气泡体积被压缩、聚并行为改变,kLa 可能略有下降。通气分布器孔径越小,初始气泡越小,a 越大,但小孔易堵塞,需权衡。
物性层面,表面张力降低有利于形成小气泡,a 增大;黏度升高会增厚液膜、降低 kL,同时抑制气泡破碎,使 kLa 明显下降,这在高密度丝状菌或高固含体系中尤为突出。培养基中的盐和蛋白会抑制气泡聚并,使 a 增大,kLa 反而可能高于纯水。消泡剂虽然能控制泡沫,但会降低表面张力并促进聚并,常导致 kLa 下降 20%~50%,使用时应重新测定 kLa。温度升高会提高氧扩散系数,但也会降低氧溶解度,净效果需实测确认。
四、放大中的核心地位与实操判断
放大时,反应器体积从实验室的几升放大到工业的几十立方米,单位体积输入功率 P/V 通常下降,气泡停留时间分布变宽,混合时间延长。若只按几何相似或 P/V 相等放大,kLa 往往无法同步保持,导致大罐供氧不足。因此,kLa 常被选为关键放大准则之一:在小试中测定不同操作条件下的 kLa,建立 kLa 与 Q、N、P/V 的关联式,再在放大罐中通过调整通气与搅拌,使目标工况下的 kLa 与小试一致或略高。
实际操作中,可先在小试罐中用动态法测定 kLa,记录不同 N 和 Q 下的数据,拟合经验式 kLa = C · (P/V)^α · Vs^β(Vs 为表观气速)。放大时根据大罐的 P/V 和 Vs 反算 kLa,再与目标值比较。若不足,优先提高通气量或调整桨型;若细胞对剪切敏感,则考虑富氧通气或提高罐压。需要强调的是,kLa 测定必须在实际培养基和实际温度下进行,因为物性差异会显著改变结果。只有把 kLa 与 OUR 匹配起来,才能保证放大后溶氧不低于临界值,细胞代谢不受限。