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-batch12506006000
酿酒酵母好氧培养6201201200
CHO 细胞培养0.258×10⁶ cells·mL⁻¹(约 4 g·L⁻¹)约 1约 10
青霉素产黄青霉33090900

从表中可读出关键结论:大肠杆菌高密度培养的 OUR 比动物细胞培养高约两个数量级。这直接决定了前者必须用富氧空气、加压或大功率搅拌,而后者往往只需低强度通气甚至表面通气。估算时务必统一单位,尤其注意细胞浓度是干重还是湿重、细胞数如何折算为干重(CHO 常按 1×10⁶ cells ≈ 0.5 μg DW 量级折算,示例中取近似值)。

需氧量计算链:比耗氧速率 → 细胞密度 → 体积氧负荷 比耗氧速率 q(O₂) 受 μ、维持代谢、 产物合成还原力影响 细胞密度 X g DW·L⁻¹ 或 cells·mL⁻¹ 折算干重 反应器体积 V 有效装液量, 非罐体公称容积 × OUR = q(O₂) × X 体积氧消耗速率 mmol·L⁻¹·h⁻¹ 总需氧速率 = OUR × V mmol·h⁻¹,供氧设计的目标负荷 OUR 量级对比(示例) 大肠杆菌 600 酵母 120 青霉 90 CHO 约 1(mmol·L⁻¹·h⁻¹)

从 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⁻¹ 量级,常规低通气即可满足。

估算时先确认三件事:qO₂ 的单位基准是细胞干重还是细胞数;X 是活细胞浓度还是总细胞浓度;V 是装液体积。三者任一错位,OUR 就可能偏差数倍,后续 kLa 与通风量计算全部失真。
高密度培养后期,qO₂ 可能因产物合成或维持代谢上升,而 X 仍在增长,OUR 呈“双因素放大”。若按前期数据线性外推供氧能力,常在诱导后或补料后期出现溶氧骤降,应预留 20%~30% 的供氧余量并设置溶氧联锁。

实际估算中,还可通过离线测量验证:在封闭或停气条件下记录溶氧下降斜率 dCL/dt,则 OUR = −dCL/dt(注意排除气液传质贡献,需在停止通气并假设 kLa 贡献可忽略的短时窗口内进行)。这一动态法得到的 OUR 可与 qO₂×X 的估算值互相校核,是工艺放大前值得做的一步。