3.3 生物反应器工程

3.3 生物反应器工程

本节摘要：SOURCE 3.3：反应器需匹配 kLa、混合、剪切；尺度放大时流体场剧变；P/V、OUR、kLa 是常用尺度无关参数。

动手：kLa 与 OUR 匹配

OUR ≈ OTR = kLa · (C* - CL)
# 若 CL→0，需 kLa ≥ OUR/C*
# 搅拌+通气提高 kLa；过高剪切损伤细胞
​

反应器类型谱系：不止搅拌罐

反应器是「生命环境模拟器」，不是容器。不同类型为不同生物与产物服务：

• 机械搅拌罐（STR）：最通用，混合与传质可独立调（搅拌→剪切、通气→kLa），用于细菌、酵母、CHO 的绝大多数工艺。
• 气升式（ALR）：靠通气驱动液体循环，无搅拌轴，剪切温和，适合剪切敏感的真菌与植物细胞，但粘度受限。
• 膜式反应器：细胞固定于膜表面或膜间，适合贴壁细胞与连续培养。
• 光生物反应器（PBR）：透明材质 + 光照系统，用于微藻光合培养，光在液层内的衰减是主要约束。
• 一次性（SUB）：聚合物袋代替不锈钢罐，省去清洗灭菌验证，灵活性高，在中试与 CGT（细胞基因治疗）中日益普及。

三角平衡：混合、传氧、剪切

反应器工程的核心矛盾可概括为"三角平衡"：

• 混合：搅拌提供混合，消除温度、pH、底物浓度梯度。但混合时间（θmix）随体积增大而变长——大型罐里，靠近加料点的细胞与远离加料点的细胞经历的底物浓度完全不同。
• 传氧：用 kLa（体积传质系数，h⁻¹）衡量反应器供氧能力。提高 kLa 的常用手段是增大搅拌功率、提高通气量、使用富氧气体、提高罐压。
• 剪切：搅拌与气泡破裂都会产生剪切力。过度搅拌会损伤细胞（尤其动物细胞与菌丝体），表现为活力下降、产物丢失，甚至膜损伤——2.3 节氨基酸厂 18% 碳异常逸出的根因就是剪切致膜损伤。

三者互相牵制：加搅拌提 kLa 就加大剪切；加大通气能提 kLa 但气泡破碎也伤细胞；改善混合要加搅拌又改变剪切场。反应器工程师的工作就是在这三角中为特定细胞找到一个可操作窗口。

搅拌生物反应器：混合—传氧—剪切三角

kLa 的测定与工程估算

实验上常用动态法测 kLa：停止通气让溶氧自然下降，记录 CL 随时间变化，再恢复通气，由 OTR 方程拟合出 kLa。工程上也有经验关联式（如 Van't Riet 式），把 kLa 与单位体积功率 P/V、表观气速 vs 关联起来：

kLa ≈ a · (P/V)^b · (vs)^c

其中指数 b、c 因体系而异（典型 b≈0.4–0.7，c≈0.4–0.6）。放大时若已知小罐的 P/V 与 vs，可外推预测大罐的 kLa，再验证 OUR 是否匹配。

# kla_scale.py —— 用经验关联式估算放大罐 kLa
def estimate_kla(pv, vs, a=0.026, b=0.4, c=0.5):
    """Van't Riet 类关联式（水/发酵液近似）"""
    return a * (pv ** b) * (vs ** c)

small = {"pv": 1.0, "vs": 0.1}   # 实验室罐 W/kg, m/s
large = {"pv": 0.8, "vs": 0.05}  # 生产罐（同 P/V 不同通气）

kla_small = estimate_kla(**small)
kla_large = estimate_kla(**large)
print(f"小罐 kLa ≈ {kla_small:.2f} h^-1")
print(f"大罐 kLa ≈ {kla_large:.2f} h^-1")
# 若大罐 OUR 需求更高而 kLa 不足，就需提高 P/V 或 vs
​

注意：关联式只做粗估。真实工艺放大必须用 CFD（计算流体力学）校核流场，并做混合时间实验确认无死角。

放大：从微升到万升的六个数量级

微升级到万升级至少 6 个数量级 尺度跨越——许多摇瓶优异性能放大即骤降。放大的本质是保持"细胞生理环境"一致，但物理环境无法等比缩放：

• 表面积/体积比下降 → 传热能力相对变弱，代谢热难以及时移除；
• 搅拌功率密度若保持恒定，叶尖速度上升 → 剪切增大；
• 混合时间延长 → 宏观梯度出现。

因此放大策略从不依赖单一参数，而是用尺度无关参数组合：P/V、kLa、叶尖速度、混合时间、OTR/OUR 余量等，按产品特性选取主控准则。例如动物细胞常用 P/V 恒定并严格控制剪切；好氧细菌发酵常用 kLa 恒定或 OTR 余量设计。

⚠️ 常见坑：只复制 rpm 和 vvm——P/V 与单位体积功率才更可比。

💡 关键直觉：反应器是「生命环境模拟器」，不是容器。

混合时间的意义与测定

混合时间（θmix）定义：在罐内某点加入示踪剂后，浓度波动衰减到设定偏差（通常 ±5%）所需时间。它在工程上的价值在于回答"加进去的东西多久才能均匀"：

• 补料、调 pH、调温，本质上都是"先混合后生效"；θmix 过长意味着加碱点附近的细胞在很长一段时间里经历着与远处完全不同的 pH。
• 放大时 θmix 随体积增大而显著变长（典型的混合时间随罐容的 1/3 到 1/2 次方增长）。生产罐里 θmix 可达数十秒到数分钟，而细胞代谢的响应时间在秒量级——宏观梯度因此不可避免。

实验测定常用 pH 脉冲法：快速注入少量碱或酸，用多点 pH 电极记录浓度均衡过程。工程上也有半经验关联式与 CFD 两种预测路径：关联式快但精度有限，CFD 能给出完整的流场与死区分布，是大型罐设计的标配手段。无论哪种方式，结论要落到"梯度是否超过细胞耐受窗口"这个工程判据上。

剪切环境的量化与缓解

剪切力的描述量级包括叶尖速度（tip speed）、能量耗散率（ε，W/kg）与剪切率（γ，s⁻¹）。不同生物敏感性差异巨大：

• 细菌、酵母：可耐受高剪切，一般不是限制因素；
• 哺乳动物细胞：无细胞壁，对剪切敏感，叶尖速度通常限制在 1–2 m/s 级别；
• 真菌菌丝体：菌丝网络易被剪切打断，但适度剪切反而有利于分散与传质——存在"最优剪切窗口"。

缓解剪切损伤的手段包括：增大叶轮直径以降低转速（保持 P/V）、采用气体分布替代高搅拌、使用一次性袋式反应器的柔性壁面、改用灌流或固定化培养避免细胞被反复夹带。

# mixing_shear.py —— 混合时间与剪切量级估算（工程粗算）
def mixing_time(v_l, pv, tip_speed):
    """经验关联：混合时间随体积与搅拌功率变化（s）"""
    theta = 8.0 * (v_l ** (1.0/3.0)) / (pv ** 0.35)
    shear_risk = "高" if tip_speed > 2.0 else "中" if tip_speed > 1.0 else "低"
    return theta, shear_risk

for vol in [5, 250, 2000, 10000]:   # 实验室→中试→生产
    theta, risk = mixing_time(vol, pv=1.0, tip_speed=1.5)
    print(f"{vol:>6} L: θmix ≈ {theta:.1f} s, 剪切风险 {risk}")
# 输出趋势：体积越大混合时间越长，剪切风险受叶尖速度主导
​

这样一组粗算就能在选罐阶段揭示"这台罐的混合能力够不够"——比盲目按体积比放大可靠得多。

核心回顾

• 反应器类型按细胞与产物特性选择
• kLa-OUR 平衡是供氧核心
• 混合/传氧/剪切三角互为约束
• 尺度无关参数指导放大
• 剪切与传质常相互制约
• 放大必须验证，不能只靠关联式外推