2.3 热力学与物料衡算

2.3 热力学与物料衡算

本节摘要：SOURCE 2.3：建立 dX/dt、dS/dt、dP/dt 衡算方程追踪碳流；某氨基酸厂发现 18% 碳异常逸出 CO₂，根因是剪切致膜损伤。

动手：碳平衡核算

反应器内生物量、底物、产物的变化速率，可由三个状态方程描述：

# 状态方程：批式反应器的简化衡算
# dX/dt = mu*X - kd*X          生物量：生长 - 死亡
# dS/dt = -(1/Yx_s)*mu*X - ms*X  底物：用于生长与维持
# dP/dt = qP*X                 产物：与生物量成正比
# 式中 mu 为比生长速率, kd 为比死亡速率,
# Yx_s 为生物量得率, ms 为维持系数, qP 为比生产速率
def batch_step(x, s, p, dt, mu, kd, yx_s, ms, qp):
    dx = (mu * x - kd * x) * dt
    ds = (-(1.0 / yx_s) * mu * x - ms * x) * dt
    dp = (qp * x) * dt
    return x + dx, s + ds, p + dp

x, s, p = 0.5, 20.0, 0.0      # 初始：生物量g/L, 葡萄糖g/L, 产物g/L
for _ in range(100):
    x, s, p = batch_step(x, s, p, dt=0.1,
                         mu=0.3, kd=0.01, yx_s=0.4, ms=0.02, qp=0.05)
    if s <= 0:
        break
print(f"终点: X={x:.2f} g/L, S={s:.2f} g/L, P={p:.2f} g/L")
​

热力学：可行性判官与效率天花板

物料衡算回答"物质去向"，热力学回答"反应能否发生、上限在哪"。两个基本判据：

• 吉布斯自由能 ΔG：ΔG < 0 的反应热力学可行。微生物几乎所有代谢反应都靠 ATP 水解驱动，而 ATP 的再生依赖底物氧化。工程上常据此判断一条代谢通路是否"合理"——如果一条设计通路的中间步骤 ΔG 大正，大概率行不通，需要换酶或换路线。
• 第二定律：任何能量转化都有损耗。葡萄糖完全氧化释放约 2870 kJ/mol，但微生物只能把其中一部分捕获为 ATP 化学能。葡萄糖理论 ATP 上限约 30–32 mol/mol，工业菌株常仅利用 60–70%——其余以热和有机酸耗散。

发酵热的工程后果是直接可感的：高密度培养中代谢热可达每立方米数十千瓦，生产罐需要冷却夹套或换热器及时移热，否则罐温失控。冷却能力有时反而成为放大后的实际限制因素——这一点在 3.3 节反应器工程中会再展开。

逐项碳流：衡算怎么闭合

以葡萄糖为唯一碳源，进入系统的碳应等于出系统的碳：

• 进碳 = 消耗的葡萄糖 × 每摩尔含碳摩尔数（葡萄糖 6 mol C/mol）
• 出碳 = 生物量碳 + 产物碳 + CO₂ 碳 + 溶解态副产物碳（乙酸、乳酸等）

闭合率 = 出碳/进碳 × 100%。正常情况下应达到 90% 以上；系统性偏低说明存在未识别的碳去向，可能来自：泄漏挥发、测量误差系统偏差、或未被测出的副产物。某氨基酸厂发现 18% 碳异常逸出 CO₂，追查根因是剪切致膜损伤——剪切损伤导致细胞膜通透性改变，糖代谢绕过正常通路直接氧化成 CO₂。这个案例说明，衡算方程是发现「碳泄漏」的镜子：当闭合率显著偏离 1，工艺内部必然发生了未预期的事件，值得停下追问。

补料分批培养的碳流闭合图

衡算在工艺设计中的三个用途

衡算方程组不是"交作业用的理论"，而是在工艺设计全周期反复调用的工具，至少有三个实际用途。

用途一：确定配比与投料量。知道理论得率与维持消耗，就能在项目初期估算"要产 100 kg 产品，需要多少底物、多少水、多大的罐"。这在工艺经济性评估（COGS 估算）中是第一笔账，误差往往来自对 ms 与 Yx/s 的低估。

用途二：识别瓶颈与异常。当 dP/dt 明显低于按 qP 计算的预期，或者碳闭合率持续偏低，说明系统里有未建模的环节。可能是副产物通路被激活、细胞死亡加剧、或者设备泄漏。衡算的闭合检验相当于"过程审计"，能抓住实验记录里看不出来的问题。

用途三：指导控制策略设计。由衡算方程可以反推补料速率的上限：底物消耗速率 = 生长消耗 + 维持消耗 + 产物消耗，三者都由比速率决定。知道了每个速率的量级，补料曲线就不是拍脑袋，而是由"速率预算"推导出来的。

用呼吸商做在线诊断

呼吸商 RQ 不需要离线取样，仅靠尾气分析（O₂、CO₂）即可在线计算，是成本最低的过程诊断手段：

# rq_monitor.py —— 由尾气数据计算 RQ 并预警
# CER = CO2 释放速率, OUR = O2 消耗速率（mol/m³/h）
def rq_from_gas(cer, our):
    return cer / our if our > 0 else float("nan")

# 某好氧批次的尾气数据（每小时一行）
gas_log = [
    (1.20, 1.18), (1.25, 1.24), (1.30, 1.28),
    (1.05, 1.02), (0.80, 0.78), (0.60, 0.58),
]
for h, (cer, our) in enumerate(gas_log):
    rq = rq_from_gas(cer, our)
    flag = "警告" if rq < 0.85 or rq > 1.15 else ""
    print(f"h={h}: CER={cer:.2f} OUR={our:.2f} RQ={rq:.2f} {flag}")
​

以葡萄糖为底物时 RQ 正常约 1.0。RQ 持续显著偏离 1.0，提示底物类型改变、产物类型改变或代谢异常——值得在第一时间排查。

能量衡算与呼吸商

除碳衡算外，能量衡算同样实用。呼吸商 RQ = CER/OUR（CO₂ 释放率/耗氧率）能在线反映代谢类型：以葡萄糖为底物时，好氧生长 RQ 约 1.0；转向合成富氧化合物时 RQ 下降；厌氧产乙醇时 RQ 可升至 1.2 以上。RQ 漂移可作为工艺异常的早期信号，属于"不看不知道，一看吓一跳"的廉价诊断手段。

系统边界的选择决定衡算结果的解释力。LCA 要求把上游农业、中游制造、下游处置纳入同一框架——动手做碳足迹需划清系统边界。只算发酵罐内平衡，可能得出"绿色"结论；把上游玉米淀粉水解的能耗算进去，"绿色"标签可能立不住。这正是 5.2 节会详述的全生命周期视角。

⚠️ 常见坑：只看发酵罐内平衡——上游玉米淀粉水解能耗可抵消「绿色」标签。

💡 关键直觉：衡算方程是发现「碳泄漏」的镜子。

要点串联

• 三方程是碳流追踪基础
• 热力学判据决定通路的可行性与上限
• 闭合率是工艺健康的体检指标
• 异常呼吸商提示物理损伤
• 系统边界决定 LCA 结论
• 热力学第二定律设定效率上限