一种基于主成分分析的合成反应进程实时监测方法

摘要: 本发明涉及化学主成分分析,具体涉及一种基于主成分分析的合成反应进程实时监测方法。背景技术:1、主成分分析作为一种经典的多变量统计过程监测方法,因其能够将高维过程...
一种基于主成分分析的合成反应进程实时监测方法

本发明涉及化学主成分分析,具体涉及一种基于主成分分析的合成反应进程实时监测方法。

背景技术:

1、主成分分析作为一种经典的多变量统计过程监测方法,因其能够将高维过程数据投影至低维特征空间并有效提取过程主要变异方向,近年来已被广泛应用于化工合成反应的在线监测领域。

2、目前,基于主成分分析的合成反应监测技术通常采用“全局主成分分析模型+固定统计量控制限”的基本架构:首先利用历史正常批次数据建立单一主成分分析模型,然后在线计算每一采样时刻的主成分得分或平方预测误差等统计量,通过判断该统计量是否超出预设控制限来实现异常状态报警。随着工业物联网技术的发展,部分改进方案开始引入多传感器数据融合策略,将温度、压力、近红外光谱和搅拌扭矩等多源信号经时间对齐后共同输入主成分分析模型;亦有研究者提出针对反应不同阶段建立多个局部主成分分析模型的思路,以期解决全局模型无法适应过程动态特性变化的难题。

3、然而,上述现有技术仍存在以下不足,现有主成分分析模型本质上是一种纯数据驱动方法,其主成分提取方向完全取决于过程数据中的最大方差方向,缺乏对化学反应内在机理的约束。

4、当传感器因环境温度变化、探头污染或电子漂移等因素产生测量偏差时,数据方差结构会发生改变,此时主成分分析模型极易将传感器漂移所引起的数据变异误判为反应状态的异常波动,造成大量虚假报警;同时,由于主成分的物理含义不明确,操作人员难以区分统计量的超限究竟是源于真实的工艺异常还是传感器自身性能退化,严重制约了监测结果在实际工业生产中的可信度。

5、对此,本技术提出一种基于主成分分析的合成反应进程实时监测方法,用以解决上述问题。

技术实现思路

1、本发明的目的是提供一种基于主成分分析的合成反应进程实时监测方法,以解决现有技术中当传感器因环境温度变化、探头污染或电子漂移等因素产生测量偏差时,数据方差结构会发生改变,此时主成分分析模型极易将传感器漂移所引起的数据变异误判为反应状态的异常波动,造成大量虚假报警;同时,由于主成分的物理含义不明确,操作人员难以区分统计量的超限究竟是源于真实的工艺异常还是传感器自身性能退化,严重制约了监测结果在实际工业生产中的可信度的问题。

2、为了实现上述目的,本发明提供如下技术方案:

3、一种基于主成分分析的合成反应进程实时监测方法,包括:

4、获取合成反应过程中多源传感器输出的异步时序信号,对所述异步时序信号进行时间对齐与采样频率归一化,得到同步多源数据矩阵;

5、根据所述同步多源数据矩阵中的温度与压力数据计算当前反应热释放速率,并通过所述反应热释放速率划分反应诱导期、主反应期和终止期,得到当前时刻的阶段类别标签;

6、基于所述阶段类别标签,调用对应的局部主成分分析模型对所述同步多源数据矩阵进行主成分提取,得到表征反应状态的低维主成分得分向量;

7、根据所述低维主成分得分向量,并依据从所述同步多源数据矩阵中实时计算得到的反应转化率与瞬时放热速率构建机理方向向量,采用givens旋转法将所述低维主成分得分向量的载荷空间向所述机理方向向量逐步旋转,旋转步长由当前角度偏差自适应决定,得到受反应动力学方向约束的增强主成分轨迹;

8、根据所述增强主成分轨迹在空间中的移动方向和收敛速度,预测未来预设时间窗口内的轨迹落点,并依据所述轨迹落点与安全边界的距离输出反应终点预警信号。

9、进一步地,所述对所述异步时序信号进行时间对齐与采样频率归一化时,包括:

10、以近红外光谱的采集时刻为全局基准时间戳,对温度和搅拌扭矩信号进行三次样条插值,并依据各传感器预设的仪器精度系数对插值后的数据进行可信度加权,得到所述同步多源数据矩阵。

11、进一步地,所述通过反应热释放速率划分反应诱导期、主反应期和终止期时,阶段判据至少包括反应热释放速率判据,并辅以近红外光谱转化率判据进行综合判定,包括:

12、根据温度与压力数据计算当前时刻的反应热释放速率,同时根据近红外光谱特征吸收峰面积计算当前时刻的实时转化率;

13、当所述反应热释放速率对时间的一阶导数由正转负且持续下降超过第一阈值,并且所述实时转化率达到预设转化率门限时,判定进入主反应期;

14、当所述反应热释放速率下降至初始峰值的第一百分比以下且所述实时转化率的变化率小于预设变化率门限时,判定进入终止期,并根据判定结果更新所述阶段类别标签。

15、进一步地,所述调用对应的局部主成分分析模型对同步多源数据矩阵进行主成分提取时,当所述阶段类别标签处于诱导期与主反应期的交界区间内,通过计算所述反应热释放速率的变化率斜率确定模糊隶属度系数;

16、依据所述模糊隶属度系数同时运行诱导期模型与主反应期模型,并将两个模型输出的初始主成分向量进行加权求和,得到平滑过渡的所述低维主成分得分向量,以消除模型硬切换引起的监测数值阶跃。

17、进一步地,所述动力学特征参数包括当前时刻的反应转化率和瞬时放热速率,其中所述反应转化率通过近红外光谱特征吸收峰面积的实时消减值换算得到,所述瞬时放热速率由当前时刻的釜内温度与夹套温差结合传热系数计算得到;

18、将所述反应转化率与瞬时放热速率的乘积作为旋转基准量,以所述旋转基准量在当前时刻之前预设时间窗口内的变化梯度方向作为机理方向向量,采用givens旋转法构建正交旋转矩阵;

19、通过所述正交旋转矩阵对所述低维主成分得分向量的原始坐标轴进行旋转,使第一主成分的载荷方向与所述机理方向向量逐步对齐,每次旋转步长不超过预设角度阈值,得到所述受反应机理约束的增强主成分轨迹。

20、进一步地,所述正交旋转矩阵的构建过程中,进一步限定仅对所述低维主成分得分向量中累积方差贡献率大于预设支撑阈值的主成分维度执行所述坐标轴旋转;

21、旋转后计算当前增强主成分轨迹的数据重构残差,若所述重构残差超出动态容差带,则进一步计算所述同步多源数据矩阵中各原始变量维度的重构残差贡献度,将贡献度最大的一个或多个变量所对应的传感器通道判定为漂移源;

22、针对判定出的漂移源所对应的传感器通道,计算该通道的漂移补偿因子,并利用所述漂移补偿因子对该通道的原始数据进行反向修正。

23、进一步地,所述预测未来预设时间窗口内的轨迹落点时,利用所述增强主成分轨迹在最近历史序列中的一阶差分与二阶差分构建带衰减因子的预测模型,当所述二阶差分超过预设幅度阈值时自动减小二阶项的权重,得到预测轨迹点集;

24、以历史正常批次的增强主成分轨迹包络线作为所述安全边界,计算所述预测轨迹点集中各点至所述安全边界的最短空间距离,并记录所述最短空间距离的时间变化序列。

25、进一步地,所述反应终点预警信号包括阶段演进提示信号和终点确认信号;

26、当所述最短空间距离的时间变化序列呈递减趋势且预测轨迹的移动方向指向安全边界外侧时,输出所述阶段演进提示信号,并缩短所述预设时间窗口的时长;

27、当所述最短空间距离在连续多个预测周期内均保持在收敛静止阈值内,且所述增强主成分轨迹的模长变化率小于终止速率门限时,输出所述终点确认信号。

28、进一步地,在获取所述同步多源数据矩阵之后、调用局部主成分分析模型之前,还根据所述同步多源数据矩阵中近红外光谱数据所对应的特征吸收峰面积,计算当前时刻关键中间体的富集变化速率;

29、利用所述富集变化速率修正所述反应热释放速率的划分边界:

30、当基于反应热释放速率判定的诱导期终止时刻与基于关键中间体富集变化速率峰值判定的诱导期终止时刻之间的偏差超过预设偏差阈值时,以所述关键中间体富集变化速率达到峰值时对应的时刻作为诱导期终止的精确锚点,并依据所述精确锚点强制校正所述阶段类别标签的更新时间。

31、进一步地,在得到所述低维主成分得分向量后,根据所述同步多源数据矩阵中不同光谱波长变量的载荷系数动态变化趋势,筛选出对当前阶段类别标签敏感度最高的特征波长子集;

32、对所述特征波长子集内的变量赋予高于默认值的协方差权重,并以增权后的协方差矩阵重新执行所述主成分提取,使更新后的低维主成分得分向量强化对中间体结构转化的响应灵敏度。

33、与现有技术相比,本发明提供的一种基于主成分分析的合成反应进程实时监测方法,通过将反应动力学特征作为物理偏置项引入主成分提取过程,对数据驱动的主成分进行机理约束下的坐标旋转,使提取出的增强主成分轨迹既能反映数据统计变异,又能遵循化学反应热力学与动力学基本规律,解决了传统主成分分析因缺乏物理约束而将传感器漂移误判为工艺异常的缺陷,显著提升了监测结果的可靠性。

34、本发明采用基于模糊隶属度的阶段模型软切换策略,并结合主成分轨迹在空间中的收敛趋势进行动态终点预测,能够在不依赖固定时间阈值的前提下自适应地识别诱导期、主反应期与终止期的模糊边界,同时通过轨迹落点与安全边界的距离演化提前发出分阶段预警信号,有效避免了因模型硬切换导致的监测数值阶跃和终点误判问题,为合成反应全过程提供了连续、平稳的智能监测决策链。