1. 生物反应器控制系统工程导论

生物反应器控制系统工程研究的不是“把传感器接上、把泵开起来”这么简单,而是围绕活体细胞或酶这类生物催化剂,在可控的物理化学环境中稳定地实现目标产物积累。它横跨生物工艺、检测仪表、自动控制与数据管理,核心任务可以概括为一句话:让反应器内的关键状态变量在扰动下仍保持在工艺允许的窗口内,并让整个过程可重复、可追溯、可放大。

课程目标与能力定位

本课程希望学习者最终具备三类能力。第一是建模与辨识能力:能把温度、pH、溶氧、补料等回路抽象成被控对象、执行器、传感器和控制器组成的反馈结构,理解增益、时间常数、纯滞后对控制品质的影响。第二是回路实现能力:能针对具体反应器选择检测原理、执行机构和控制策略,完成参数整定,并判断振荡、超调、积分饱和等异常。第三是系统集成与放大能力:能把单元控制回路组织成符合 GMP 或工业数据规范的监控系统,理解规模放大时混合、传质、传热时间常数的变化如何改变控制难度。

生物过程控制的特殊性

与化工过程相比,生物过程的控制难点来自“活体”本身。首先是时变与非线性:细胞在延滞期、对数生长期、稳定期和衰亡期的比生长速率、耗氧速率、产物合成速率都不同,同一套 PID 参数在批次后期可能不再适用。其次是可测性差:葡萄糖、产物浓度、细胞活性等关键变量往往无法在线直接测量,只能通过溶氧、pH、尾气 CO₂、OUR/CER 等间接量推断,形成“软测量”问题。第三是无菌约束:传感器必须可灭菌、可校准,取样会带来污染风险,因此在线检测点位有限。第四是多变量耦合:通气量同时影响溶氧和 pH(因 CO₂ 剥离),补料同时影响底物浓度和体积,温度又影响溶氧饱和度,回路之间不能孤立整定。

理解这些特殊性,是后续选择控制策略的前提:不是所有变量都适合闭环,也不是所有闭环都用 PID。对滞后大、可测性差的变量,往往采用“前馈+反馈”或“模型预测”思路;对无菌要求高的场合,优先选择非接触式或原位灭菌型传感器。

系统边界定义

做控制系统设计,第一步是明确边界。通常把生物反应器控制系统划分为三层:现场层包括反应器本体、搅拌、夹套/盘管、通气、补料泵、酸碱泵、消泡剂泵以及温度、pH、溶氧、压力、液位等传感器;控制层包括 PLC 或嵌入式控制器、PID 回路、联锁与报警逻辑;监控层包括 SCADA/上位机、批次记录、趋势曲线、电子批记录与数据归档。边界之外是公用工程(蒸汽、冷却水、压缩空气、纯水)和下游处理,它们通过接口影响反应器,但不在本课程控制回路内部。

生物反应器控制系统边界与信号流 现场层(反应器与执行器) 搅拌电机 / 转速 影响 kLa 与混合时间 夹套 / 盘管 温度 冷却水阀 / 电加热 通气 / 尾气 空气、O₂、N₂、CO₂ 补料 / 酸碱 / 消泡 蠕动泵、计量泵 传感器:T、pH、DO、P、液位 控制层(PLC / 嵌入式) PID 回路:T / pH / DO 比例-积分-微分运算 前馈 / 串级 / 比值 补料与 DO 联动 联锁与报警逻辑 超温、超压、低液位 批次阶段管理 灭菌、接种、发酵、收获 信号调理与 A/D 转换 4–20 mA / 0–10 V / 数字量 采样周期 0.5–5 s 监控层(SCADA / 上位机) 趋势曲线与实时画面 T、pH、DO、OUR、CER 批次记录与电子批记录 审计追踪、数据完整性 报警管理与报表 偏差、放行判断 数据归档与接口 MES / LIMS / 历史库 工艺配方与参数下发 设定值、报警限、阶段切换 权限与电子签名 测量 数据 控制 指令

典型反应器类型与控制需求概览

不同反应器的流体力学和传质特性不同,控制需求也不同。下表给出常见类型的对比,参数为教学示例,实际设计需以设备手册和工艺验证为准。

反应器类型 典型规模 关键控制变量 主要难点 常用策略
机械搅拌釜(STR) 1–200 m³ T、pH、DO、补料 放大后混合时间变长,DO 梯度大 串级 DO、前馈补料
气升式反应器 0.5–50 m³ 通气量、DO、pH 剪切低但传质受气速限制 通气-搅拌比值控制
鼓泡塔 1–100 m³ 气速、液位、DO 返混强,pH 分布不均 分段 pH 检测、多点 DO
一次性生物反应器 0.05–2 m³ T、pH、DO、袋压 传感器一次性、校准窗口短 预校准、单点校正
微反应器 / 微流控 µL–mL 流速、温度、停留时间 检测点少、响应快 流量闭环、快速 PID

从“是什么”到“怎么做”:一个温度回路的示例

以 5 L 机械搅拌釜温度控制为例,说明控制工程的实际落地。被控对象是反应液温度 T,执行器是夹套冷却水比例阀,传感器是 Pt100。假设工艺要求 37.0 ± 0.2 ℃,冷却水入口 12 ℃,反应放热使温度有上升趋势。

第一步是确定采样周期。温度对象时间常数约 8–15 min,纯滞后约 20–40 s,采样周期取 2 s 足够,既能捕捉变化,又不会引入过多噪声。第二步是选择控制律。由于滞后不大,先用 PI 控制,微分项容易放大传感器噪声,暂不启用。第三步是整定参数。可用阶跃响应法:手动将阀门从 30% 阶跃到 40%,记录温度响应,得到增益 K≈0.8 ℃/%、时间常数 τ≈600 s、滞后 θ≈30 s。按 Ziegler-Nichols 反应曲线法,PI 参数可初设为 Kp≈1.2τ/(Kθ)=1.2×600/(0.8×30)=30,Ti≈2.2θ≈66 s。实际投运时先降 Kp 到 15–20,观察超调和振荡,再逐步收紧。

第四步是观察与判断。投运后看三条曲线:设定值、实际温度、阀门开度。若温度在设定值附近等幅振荡,说明 Kp 偏大或 Ti 偏小;若温度缓慢爬升且阀门已全开仍不降,说明冷却能力不足或存在积分饱和;若温度稳定但阀门频繁小幅动作,说明 Kp 偏大或采样噪声大,可加滤波或死区。最终目标是温度波动控制在 ±0.2 ℃ 内,阀门动作平稳,无持续振荡。

教学示例中的 K、τ、θ 和 PID 参数均为合理假设,用于说明整定思路。真实反应器的参数必须通过现场阶跃试验或系统辨识获得,不能直接套用。

本课程的学习路径