2. 稳态培养的数学基础
连续培养的核心问题可以归结为一句话:在持续进料与出料的条件下,如何让反应器内的细胞浓度和限制性底物浓度同时保持不随时间变化。要回答这个问题,必须先把几个基本量定义清楚,再用物料衡算把它们联系起来。本节从最基本的比生长速率出发,逐步推导出连续培养的稳态方程,并说明Monod方程在其中扮演的角色。
一、比生长速率与稀释率
比生长速率 μ 定义为单位时间内单位细胞量所增加的细胞量,即 μ = (1/X)·(dX/dt),单位是 h⁻¹。它描述的是“每个细胞平均每小时的增殖能力”,而不是绝对增长量。在分批培养中,μ 随底物浓度和环境条件变化;在连续培养中,稳态时 μ 被外部操作条件锁定,这正是连续培养可控性的来源。
稀释率 D 定义为进料体积流量 F 与反应器有效体积 V 之比:D = F/V,单位同样是 h⁻¹。它的物理含义是“每小时被新鲜培养基替换掉的反应器体积比例”。D 的倒数就是水力停留时间 τ = 1/D = V/F,表示一个流体微元在反应器内的平均停留时间。例如 V = 5 L、F = 1 L/h 时,D = 0.2 h⁻¹,τ = 5 h。这三个量是连续培养最基础的参数,后续所有推导都围绕它们展开。
二、连续培养的物料衡算
对恒化器做细胞和限制性底物的衡算,是推导稳态方程的关键。先约定参照系:以反应器内液体体积 V 为控制体,进料和出料体积流量相等(F_in = F_out = F),因此 V 恒定,这是连续培养能维持稳态的前提。
对细胞做衡算:积累速率 = 流入 − 流出 + 净生长。进料中通常不含细胞(无菌培养基),所以流入项为 0;流出项为 F·X;净生长项为 μ·X·V。于是
V·(dX/dt) = 0 − F·X + μ·X·V
两边除以 V,并代入 D = F/V,得到细胞衡算的微分方程:
dX/dt = (μ − D)·X
对限制性底物做衡算:积累速率 = 流入 − 流出 − 消耗。进料底物浓度为 S₀,反应器内浓度为 S,消耗速率用比消耗速率 q_S 表示,则
dS/dt = D·(S₀ − S) − q_S·X
这两个方程就是连续培养的基本方程。它们把操作参数 D、进料浓度 S₀ 与状态变量 X、S 联系起来,是后续所有稳态分析的基础。
三、稳态条件与基本方程推导
稳态意味着 dX/dt = 0 且 dS/dt = 0。先看细胞衡算:dX/dt = (μ − D)·X = 0。由于稳态下 X > 0,必然有
μ = D
这是连续培养最重要的结论:稳态时细胞的比生长速率恰好等于稀释率。操作者通过调节 F 来设定 D,也就等于设定了细胞被迫维持的 μ。若 D 大于当前条件下 μ_max,细胞被洗出,X 会持续下降直至趋近于零,这就是“洗出”现象。
再看底物衡算:dS/dt = D·(S₀ − S) − q_S·X = 0。引入细胞对底物的得率系数 Y(定义为生成单位细胞量所消耗的底物量,Y = −ΔX/ΔS),则 q_S = μ/Y。代入 μ = D,得到
X = Y·(S₀ − S)
这个式子说明:稳态细胞浓度由进料底物浓度 S₀ 和残余底物浓度 S 之差决定,再乘以得率系数 Y。S 越小,X 越大,但 S 不可能为零,因为细胞需要一定底物维持生长。
四、Monod方程与稳态解
Monod方程把比生长速率与限制性底物浓度联系起来:
μ = μ_max · S / (K_s + S)
其中 μ_max 是最大比生长速率,K_s 是半饱和常数,数值上等于 μ = μ_max/2 时的底物浓度。把稳态条件 μ = D 代入 Monod 方程,可以解出稳态残余底物浓度:
D = μ_max · S / (K_s + S)
S = K_s · D / (μ_max − D)
这个式子非常有用:只要知道 μ_max 和 K_s,就能根据设定的 D 直接算出稳态 S,再代入 X = Y·(S₀ − S) 算出稳态 X。例如取 μ_max = 0.5 h⁻¹、K_s = 0.1 g/L、Y = 0.5 g/g、S₀ = 10 g/L,当 D = 0.2 h⁻¹ 时,S = 0.1×0.2/(0.5−0.2) ≈ 0.067 g/L,X = 0.5×(10−0.067) ≈ 4.97 g/L。这组示例参数说明:D 越接近 μ_max,S 上升越快,X 下降越快,系统越接近洗出边界。
五、参数关系与操作判断
把上述关系整理成一张参数影响表,有助于快速判断操作方向:
| 操作变化 | D 变化 | 稳态 S | 稳态 X | 判断 |
|---|---|---|---|---|
| 增大 F | D 增大 | 升高 | 降低 | 接近洗出 |
| 减小 F | D 减小 | 降低 | 升高 | 底物更充分 |
| 提高 S₀ | D 不变 | 不变 | 升高 | X 随 S₀ 线性增加 |
| D ≥ μ_max | — | 趋近 S₀ | 趋近 0 | 洗出 |
实际设定 D 时,通常取 μ_max 的 40%~70% 作为安全操作区间,既保证较高细胞浓度,又留出抵抗扰动的余量。若 D 设定过高,底物来不及被消耗就随出料流失,S 升高、X 下降,最终洗出;若 D 设定过低,虽然 X 高,但设备产能利用率低,且可能因维持能耗占比增大而降低得率。理解 μ = D、X = Y·(S₀ − S) 和 Monod 方程这三条关系,就能在给定 μ_max、K_s、Y、S₀ 的条件下,快速估算稳态 X 和 S,并判断当前 D 是否处于安全操作范围。