3. 高密度发酵的化学计量与能量学
高密度发酵的核心矛盾，是细胞在极高生物量浓度下对碳源、氧和能量的需求被同步放大。若只凭经验补料，很容易出现“碳源够、氧不够”或“氧够、能量被无效消耗”的局面。要定量把握这一过程，必须建立一套从碳源分配到氧耗、再到热量释放的化学计量与能量学估算框架。本节围绕这条主线，给出可操作的估算方法。
一、碳源分配：从消耗到产物的物料平衡
碳源进入反应器后，主要流向四个去向：合成菌体（生物量）、合成产物、氧化为CO₂、以及以副产物或胞外代谢物形式残留。以葡萄糖为碳源、好氧发酵为例，可用简化的宏观化学计量式描述：
a·C6H12O6 + b·O2 + c·NH3
→ d·CH(1.8)O(0.5)N(0.2) + e·CO2 + f·H2O + g·产物其中 CH(1.8)O(0.5)N(0.2) 是典型细菌生物量的元素近似式。系数 a、b、c、d 等由实验测得的比消耗速率和比生成速率确定。工程上更常用“碳源分配系数”来量化：
- 生物量碳转化率 YX/S：每克碳源生成多少克菌体（以碳计或干重计）；
- 产物碳转化率 YP/S：每克碳源生成多少克目标产物；
- CO₂ 碳分率：被完全氧化损失的碳比例，直接反映能量浪费程度。
例如，某重组大肠杆菌高密度发酵中，实测葡萄糖比消耗速率 qS=0.35 g/(g·h)，比生长速率 μ=0.18 h⁻¹，菌体碳含量约 0.47 g C/g DCW，葡萄糖碳含量 0.40 g C/g。则用于合成的碳分率约为 (0.18×0.47)/(0.35×0.40)≈0.60，意味着约 40% 的碳进入氧化或副产物途径。这一比例是判断补料策略是否合理的第一把尺子。
二、ATP维持：高密度下的“隐形消耗”
细胞用于维持渗透压、修复大分子、维持膜电位的能量称为维持能，常用维持系数 mATP 表示，单位为 mmol ATP/(g DCW·h)。高密度发酵中，维持能占总 ATP 消耗的比例随菌体浓度升高而显著上升，因为单位体积内细胞数量多，维持负担被放大。
ATP 消耗可拆为三部分：生长相关消耗、产物合成消耗、维持消耗。以 mmol ATP/(g DCW·h) 计：
q_ATP = μ/Y_ATP^max + q_P/Y_ATP^P + m_ATP其中 YATPmax 为最大生长得率（约 10–12 g DCW/mol ATP），YATPP 为产物合成得率。若某工艺 μ=0.15 h⁻¹、mATP=2.5 mmol/(g·h)，则维持项已占 qATP 的相当比例。降低维持消耗的常用手段包括：控制比生长速率在较低水平、避免温度与 pH 剧烈波动、减少无效循环。但需注意，μ 过低会降低比生产强度，需在“维持占比”与“体积产率”之间取平衡。
三、氧耗估算：OUR 与 OTR 的匹配
氧耗速率 OUR 可由碳源消耗速率和氧的化学计量系数估算。对葡萄糖完全氧化，1 mol 葡萄糖理论耗氧 6 mol。实际发酵中，可用下式快速估算：
OUR = q_S × X × (b/a) × (1/32) [g O2/(L·h)]其中 X 为菌体浓度（g DCW/L），b/a 为每摩尔碳源对应的耗氧摩尔数。若 qS=0.35 g/(g·h)、X=80 g/L、b/a≈3.2 mol O₂/mol 葡萄糖，则 OUR≈0.35×80×3.2×(1/180)×32≈15.9 g O₂/(L·h)。这一数值必须与反应器的氧传递速率 OTR 匹配：
OTR = k_L a × (C* - C_L)高密度发酵中 CL 通常控制在临界氧浓度以上（约 20%–30% 饱和），若 kLa 不足，OUR 会被迫下降，导致碳源溢流、副产物积累。因此，氧耗估算的落脚点是：先算 OUR，再反推所需 kLa，最后通过通气量、罐压、搅拌转速实现。
四、热量释放：从氧耗到冷却负荷
好氧发酵的热量释放与氧耗几乎成正比，因为大部分自由能以氧为最终电子受体释放。工程上常用近似值：每消耗 1 mol O₂ 释放约 440–480 kJ 热量，取 460 kJ/mol O₂ 作为估算基准。将 OUR 换算为摩尔氧：
Q = OUR / 32 × 460 [kJ/(L·h)]以前述 OUR≈15.9 g O₂/(L·h) 为例，Q≈15.9/32×460≈228 kJ/(L·h)。若反应器体积 5 m³，则总放热约 1.14×10⁶ kJ/h，相当于约 317 kW。这一量级必须由冷却系统带走，否则温度会迅速失控。实际估算中还需叠加搅拌产热和通气带出的显热，通常取 10%–20% 余量。
五、参数关系与判断依据
将上述关系整理为可操作的判断表，便于在工艺开发中快速定位问题：
| 观察指标 | 估算关系 | 异常方向 | 可能原因 |
|---|---|---|---|
| 碳源分配碳分率 | YX/S、CO₂ 碳分率 | CO₂ 分率 > 50% | 供氧不足、溢流代谢 |
| ATP 维持占比 | mATP / qATP | 维持占比 > 40% | 高密度、环境胁迫 |
| OUR 与 OTR | OUR = qS·X·b/a/32 | OUR > OTR | kLa 不足、CL 过低 |
| 热量释放 Q | Q = OUR/32×460 | 冷却负荷超设计 | 冷却面积或水温不足 |
实际执行时，建议每批次记录 qS、X、OUR、Q 和冷却水进出口温差，形成趋势曲线。当 OUR 曲线出现平台而碳源补加速率未降时，往往意味着氧传递已成为限制因素；此时应优先提高 kLa，而不是继续增加补料。若维持能占比持续偏高，可尝试降低比生长速率或采用分阶段温控策略，但需同步评估体积产率是否仍满足要求。
掌握碳源分配、ATP 维持、氧耗与热量释放的定量估算，本质上是把高密度发酵从“看现象”推进到“算平衡”。只有碳、能、氧、热四条线同时闭合，才能在高菌体浓度下维持稳定、可重复的生产性能。