3. 反应釜升温阶段温度程序监控与超调识别

升温程序监控对象与判据设计

环氧树脂合成在投料完成后进入升温阶段，反应釜温度通常按“阶梯升温 + 保温”程序执行。本例监控对象是釜内温度 T_reactor（由插入式 Pt100 测得）与夹套温度 T_jacket（导热油/蒸汽侧），并以程序设定值 SP 作为比较基准。核心判据有三类：

• 升温速率偏差：实际速率 v = ΔT/Δt 与程序速率 v_SP 比较，偏差超过 ±20% 且持续 ≥ 2 min 记为速率异常。
• 温度超调：进入保温段后 T_reactor 超过 SP 的幅度 ΔT_os = T_max − SP，ΔT_os > 2.0 ℃ 即判为超调事件。
• 超调持续时间：从首次越过 SP+2.0 ℃ 到回落至 SP+0.5 ℃ 以内的时间 t_os，用于评估热惯性影响。

采样周期取 5 s，程序段以“目标温度 + 升温速率 + 保温时长”描述。下面给出可直接运行的监控脚本，输入为按时间排列的温度采样序列。

监控脚本实现

脚本按程序段解析设定值，逐点计算实际升温速率、速率偏差、超调幅度与持续时间，并输出事件记录。示例数据为教学构造，用于演示判据触发逻辑。

import math

# 程序段：目标温度(℃), 升温速率(℃/min), 保温时长(min)
program = [
    {"target": 60.0, "rate": 1.0, "hold": 10},
    {"target": 90.0, "rate": 0.8, "hold": 15},
    {"target": 120.0, "rate": 0.6, "hold": 20},
]

SAMPLE_DT = 5.0 / 60.0          # 采样周期 5 s，换算为 min
RATE_TOL = 0.20                 # 速率偏差 ±20%
RATE_HOLD_MIN = 2.0             # 速率异常持续 ≥2 min
OS_LIMIT = 2.0                  # 超调阈值 ℃
OS_RECOVER = 0.5                # 回落判定带 ℃

def build_sp_curve(program, dt):
    """生成 (时间min, SP) 序列，含升温段与保温段"""
    sp, t = [], 0.0
    cur = 25.0
    for seg in program:
        dur = (seg["target"] - cur) / seg["rate"]
        n = max(1, int(round(dur / dt)))
        for i in range(1, n + 1):
            sp.append((t + i * dt, cur + (seg["target"] - cur) * i / n))
        t += n * dt
        cur = seg["target"]
        nh = max(1, int(round(seg["hold"] / dt)))
        for i in range(1, nh + 1):
            sp.append((t + i * dt, cur))
        t += nh * dt
    return sp

def monitor(sp_curve, t_reactor):
    """t_reactor: 与 sp_curve 等长的实际温度序列"""
    events = []
    rate_bad_start = None
    os_start = None
    os_peak = -1e9
    for i in range(1, len(sp_curve)):
        t, sp = sp_curve[i]
        t_prev, sp_prev = sp_curve[i - 1]
        v_sp = (sp - sp_prev) / SAMPLE_DT
        v_act = (t_reactor[i] - t_reactor[i - 1]) / SAMPLE_DT
        # 速率偏差（仅升温段）
        if v_sp > 1e-6:
            dev = (v_act - v_sp) / v_sp
            if abs(dev) > RATE_TOL:
                if rate_bad_start is None:
                    rate_bad_start = t
            else:
                if rate_bad_start is not None and t - rate_bad_start >= RATE_HOLD_MIN:
                    events.append({"type": "rate", "start": rate_bad_start,
                                   "end": t, "dur": t - rate_bad_start})
                rate_bad_start = None
        # 超调识别（保温段）
        if abs(v_sp) < 1e-6:
            if t_reactor[i] > sp + OS_LIMIT:
                if os_start is None:
                    os_start = t
                os_peak = max(os_peak, t_reactor[i])
            elif os_start is not None and t_reactor[i] <= sp + OS_RECOVER:
                events.append({"type": "overshoot", "start": os_start,
                               "end": t, "dur": t - os_start,
                               "amp": os_peak - sp})
                os_start, os_peak = None, -1e9
    return events

# 教学示例：构造带超调的实际温度（示例数据，非实测）
sp_curve = build_sp_curve(program, SAMPLE_DT)
t_reactor = []
for i, (t, sp) in enumerate(sp_curve):
    lag = 1.5 * math.exp(-0.02 * i) * math.sin(0.15 * i)   # 简化热惯性扰动
    t_reactor.append(sp + lag)

for e in monitor(sp_curve, t_reactor):
    if e["type"] == "overshoot":
        print(f"超调事件: 幅度 {e['amp']:.2f} ℃, 持续 {e['dur']:.1f} min")
    else:
        print(f"速率异常: {e['start']:.1f}~{e['end']:.1f} min, 持续 {e['dur']:.1f} min")

结果判读与工艺评估记录

脚本输出的事件记录需与工艺窗口对照：超调幅度 ΔT_os 反映热惯性大小，持续时间 t_os 反映夹套热量撤除的响应快慢。若 ΔT_os 接近或超过 2.0 ℃ 且 t_os 偏长，说明升温末段夹套温度设定偏高或保温切换滞后，应在下一批次下调夹套目标温度或提前切换。速率异常则提示导热油流量波动或测温滞后，需核查夹套回路。