4. 温度控制系统中的静态与动态特性辨识
在精馏塔温度控制中,塔顶与塔釜温度通道的静态与动态特性决定了控制器的整定方向。静态特性描述稳态下操纵变量(如回流量、再沸器蒸汽量)与温度之间的增益关系;动态特性则描述温度如何随时间响应操纵变量的变化。工程上最常用的动态模型是一阶加纯滞后(FOPDT)模型:
G(s) = K · e^(-τs) / (T·s + 1)
K —— 过程增益(℃ / % 或 ℃ / (m³/h))
T —— 时间常数(min 或 s)
τ —— 纯滞后时间(min 或 s)
这三个参数分别对应三个可观察的物理现象:K 决定阶跃后温度最终变化多少,T 决定温度变化快慢,τ 决定操纵变量动作后温度多久才开始响应。精馏塔温度通道的典型特点是:增益中等、时间常数较大(数分钟到数十分钟)、纯滞后明显(受塔板数、回流路径和测量点位置影响),因此辨识时必须保证阶跃幅度足够大、观察时间足够长。
阶跃响应测试的实施要点
辨识的第一步是设计一次可用的阶跃测试。以塔顶温度通道为例,假设正常工况下回流量稳定在 12 m³/h,塔顶温度为 78.5 ℃。测试时保持再沸器蒸汽量、进料量和进料组成不变,将回流量从 12 m³/h 阶跃增加到 14 m³/h(约 17% 的阶跃),并持续记录塔顶温度直至其重新稳定。关键操作要求如下:
- 基线稳定:阶跃前至少观察 10~15 分钟,确认温度波动小于 0.2 ℃,否则基线漂移会污染辨识结果。
- 单变量原则:测试期间不得同时调整蒸汽量、进料或采出量,否则无法把温度变化归因于回流量。
- 采样与记录:温度采样周期建议取预估时间常数的 1/10 左右,例如预估 T≈8 min,则每 30~60 s 记录一次,并同步记录回流量实际值以确认阶跃是否到位。
- 观察时长:至少记录到温度变化达到总变化量的 95% 以上,通常需要 3T + τ 以上的时间。
下图给出了一次塔顶温度阶跃测试的响应曲线与 FOPDT 参数在图上的对应位置,用于说明 K、T、τ 如何从数据中读出。
从数据估计模型参数
工程上最常用的是两点法(也称切线法或 63.2% 法)。以图中示例数据为例:回流量阶跃 Δu = +2 m³/h,塔顶温度由 78.5 ℃ 降至 77.3 ℃,则稳态增益 K = ΔT/Δu = (−1.2)/2 = −0.6 ℃/(m³/h)。负号表示回流量增大使塔顶温度下降,这一符号约定必须在后续控制器整定中保持一致。
时间常数 T 的读取方法是:找到温度变化达到总变化量 63.2% 的时刻。总变化 1.2 ℃,63.2% 即 0.758 ℃,对应温度 78.5 − 0.758 ≈ 77.74 ℃。若该点出现在阶跃后第 9.5 min,而温度开始响应的时刻为阶跃后第 2.5 min,则纯滞后 τ = 2.5 min,时间常数 T = 9.5 − 2.5 = 7.0 min。于是塔顶温度通道模型为:
G_top(s) = -0.6 · e^(-2.5s) / (7.0·s + 1) (时间单位:min)
塔釜温度通道的辨识方法完全相同,但对象特性差异明显。再沸器蒸汽量对塔釜温度的影响通常增益更大、纯滞后更长(蒸汽经再沸器传热到塔釜液相测温点存在热传导与混合滞后)。例如某示例中蒸汽量阶跃 +5%,塔釜温度由 102.0 ℃ 升至 103.5 ℃,K = 1.5/5 = +0.3 ℃/%,τ ≈ 4 min,T ≈ 12 min,则:
G_bottom(s) = 0.3 · e^(-4s) / (12·s + 1)
对比两个通道可见:塔顶通道增益为负、时间常数较小、滞后较短,响应相对快;塔釜通道增益为正、时间常数大、滞后长,响应迟缓。这一差异直接决定了塔顶温度回路可以整定得稍快,而塔釜温度回路必须留出更大的积分时间并适当降低比例增益,否则容易振荡。
模型建立后的校验
参数估计完成后不能直接投入使用,应做一次校验:用另一组不同幅度的阶跃数据(例如回流量 +1 m³/h)代入模型预测温度变化过程,与实测曲线对比。若预测的稳态变化量与实测偏差在 10% 以内、63.2% 点时刻偏差在 1 min 以内,可认为模型满足控制整定精度要求;若偏差过大,应检查测试期间是否存在进料扰动、采样是否过疏或阶跃是否未真正到位。校验通过的模型即可作为塔顶、塔釜温度控制器参数整定与后续前馈补偿设计的依据。