3. 需氧量计算
通气与供氧计算的起点不是反应器,而是细胞本身。细胞在单位时间内消耗多少氧,决定了后续通风量、搅拌功率和传质面积的下限。描述这一需求最基础的参数是细胞比耗氧速率(specific oxygen uptake rate,常记作 qO₂ 或 OURspec),单位为 mol O₂·(g 细胞干重)-1·s-1,工程上更常用 mmol O₂·(g DW)-1·h-1。它乘以反应器内活性细胞浓度 X(g DW·L-1),就得到体积氧消耗速率 OUR(mol O₂·L-1·h-1):
OUR = q(O₂) × X
总需氧量速率 = OUR × V = q(O₂) × X × V
这个乘法关系看似简单,却是全部供氧设计的核心:氧需求是“单位细胞速率”与“细胞数量”的乘积。同一株菌在低密度和高密度培养中,qO₂ 可能相近,但 OUR 会相差一个数量级,这正是高密度发酵必须强化供氧的根本原因。
比耗氧速率由什么决定
qO₂ 不是常数,它随代谢状态、比生长速率 μ、碳源类型和产物合成途径而变化。工程上常用两条近似关系来估算:
- 维持-生长模型:qO₂ = mO₂ + (1/YO₂/X)·μ。其中 mO₂ 为维持系数,YO₂/X 为氧对细胞的表观得率。μ 越高,用于合成 ATP 的氧化越多,qO₂ 越大。
- 产物耦合项:当产物合成需要大量 NAD(P)H 或 ATP 时,应额外计入产物项,例如 qO₂ = mO₂ + μ/YO₂/X + qp/YO₂/P,qp 为比产物生成速率。
因此,产物生成对氧需求的影响有两条路径:一是产物合成本身消耗还原力,直接抬高 qO₂;二是产物积累常伴随细胞生长受抑、维持代谢占比上升,使单位细胞的氧耗不降反升。典型例子是青霉素发酵,产物合成期 μ 很低,但 qO₂ 仍维持在较高水平,供氧不能按生长速率简单下调。
典型体系的 OUR 估算案例
下面给出几组教学示例参数(属合理假设,非实测数据),用于说明量级差异。设反应器有效体积 V = 10 L,按 OUR = qO₂ × X 计算体积氧负荷,再乘 V 得总需氧速率。
| 体系 | q(O₂) / mmol·g⁻¹·h⁻¹ | 细胞浓度 X / g·L⁻¹ | OUR / mmol·L⁻¹·h⁻¹ | 10 L 总需氧 / mmol·h⁻¹ |
|---|---|---|---|---|
| 大肠杆菌高密度 fed-batch | 12 | 50 | 600 | 6000 |
| 酿酒酵母好氧培养 | 6 | 20 | 120 | 1200 |
| CHO 细胞培养 | 0.25 | 8×10⁶ cells·mL⁻¹(约 4 g·L⁻¹) | 约 1 | 约 10 |
| 青霉素产黄青霉 | 3 | 30 | 90 | 900 |
从表中可读出关键结论:大肠杆菌高密度培养的 OUR 比动物细胞培养高约两个数量级。这直接决定了前者必须用富氧空气、加压或大功率搅拌,而后者往往只需低强度通气甚至表面通气。估算时务必统一单位,尤其注意细胞浓度是干重还是湿重、细胞数如何折算为干重(CHO 常按 1×10⁶ cells ≈ 0.5 μg DW 量级折算,示例中取近似值)。
从 OUR 到供氧条件的衔接
得到 OUR 后,下一步是判断现有供氧能力是否匹配。稳态下氧传递速率 OTR 必须等于 OUR,而 OTR = kLa·(C* − CL)。因此可反推所需 kLa:
k(L)a(需求) = OUR / (C* − C_L)
以大肠杆菌示例为例,OUR = 600 mmol·L⁻¹·h⁻¹ = 0.6 mol·m⁻³·h⁻¹。若 37 ℃ 下 C* 取约 0.2 mol·m⁻³,允许 CL 不低于临界氧浓度 0.02 mol·m⁻³,则驱动力约 0.18 mol·m⁻³,所需 kLa ≈ 3.3 h⁻¹。这个数值对普通搅拌罐偏高,提示需要提高转速、通气量或采用富氧。反之,CHO 示例 OUR 约 1 mmol·L⁻¹·h⁻¹,所需 kLa 仅约 0.006 h⁻¹ 量级,常规低通气即可满足。
实际估算中,还可通过离线测量验证:在封闭或停气条件下记录溶氧下降斜率 dCL/dt,则 OUR = −dCL/dt(注意排除气液传质贡献,需在停止通气并假设 kLa 贡献可忽略的短时窗口内进行)。这一动态法得到的 OUR 可与 qO₂×X 的估算值互相校核,是工艺放大前值得做的一步。