一种基于递推非线性部分最小二乘的过程监控方法技术

本发明专利技术涉及一种基于递推非线性部分最小二乘的过程监控方法。传统方法采用固定的统计模型和控制限，不能及时进行模型更新从而跟踪过程的变化。本发明专利技术方法首先基于过程数据库建立基于非线性部分最小二乘的过程监控模型，同时计算控制限，然后将建立的模型应用到在线采集的实时工业过程数据，对过程进行监控。当获得新数据后，采用递推方法对模型和控制限进行更新，从而适应过程的变化。本发明专利技术提出的过程监控方法弥补了传统非线性多变量统计过程监控方法不能处理时变性的不足，可以利用新数据更新模型参数和控制限，从而适应过程的变化，提高了监控性能。

【技术实现步骤摘要】

本专利技术属于信息自动化
，涉及一种基于递推非线性部分最小二乘的过程 监控方法。
技术介绍
现代化的工业过程系统正朝着大规模、复杂性的方向发展，系统一旦发生故障就 可能造成人员和财产的巨大损失以及环境危害。这就要求对生产过程进行监控，从而确定 生产状态是否正常，出现异常工况时能及时采取相应的补救措施，减少生产停车时间，最大 限度地保障系统的安全性和可靠性。随着自动化、计算机网络及数据库技术的发展，工厂可 以直接从生产过程获得大量的实时运行数据。但是要从观测数据中实现对过程运行情况的 评估，已超出了工程师或操作员的能力范围。基于数据驱动的多变量统计过程监控技术不 依赖对象的解析模型，而充分利用工厂丰富的数据资源，可以有效的处理数据量大、数据维 数高、数据共线性等问题，在工业上得到广泛应用。很多实际工业过程具有较强的非线性， 并且随着时间推移，对象的特性和工作点都可能发生变化。传统的非线性多变量统计过程 监控方法都是采用固定的统计模型和控制限，当过程特征或操作条件发生变化时不能及时 进行模型更新从而跟踪过程的变化，从而产生漏报和误报的情况。因此，迫切需要一种既能 处理过程非线性又能处理时变性的过程监控方法。
技术实现思路
本专利技术的目的就是针对现有的过程监控技术的不足之处，提供一种基于递推非线 性部分最小二乘的过程监控方法。该方法弥补了传统非线性多变量统计过程监控方法不能 处理时变性的不足，可以利用新数据更新模型参数和控制限，从而适应过程的变化，提高了 监控性能。本专利技术方法采用数据采集、过程辨识、数据驱动等手段，首先基于过程数据库建立 基于非线性部分最小二乘的过程监控模型，同时计算控制限，然后将建立的模型应用到在 线采集的实时工业过程数据，对过程进行监控。获得新数据后，采用递推方法对模型和控制 限进行更新，从而适应过程的变化。本专利技术方法的具体步骤是步骤(1)基于过程数据库建立基于非线性部分最小二乘(NLPLS)过程监控模型， 同时计算控制限，具体方法是a.通过数据采集装置采集实时过程运行数据，将采集的实时过程运行数据作为数 据驱动的样本集合，表示为&(1)，7(1)}，其中乂1)表示第i组输入数据，y(i)表示第i组 输出数据；将输入数据构成输入矩阵X、将输出数据构成输出矩阵Y ；b.基于输入、输出数据建立非线性部分最小二乘模型，方法如下对矩阵X和Y进行归一化处理，使之均值为0，方差为1 ；将输入矩阵进行列扩展， 扩展项为径向基函数(RBF)神经网络的隐节点输出矩阵G和元素全为1的列向量1，其中G的每一行对应一个输入向量作用下的隐节点的输出，隐节点的偏置项系数为1 ；对如下增 广输入矩阵和输出矩阵进行部分最小二乘(PLS)回归{，Y}，得到的NLPLS过程监控模型表示为V" 式中，Xe表示增广输入矩阵，A和H分别为对应原始输入向量和对应RBF网络隐节 点输出向量的权值系数矩阵，b为输出偏置向量，T表示转置。式⑴中的未知参数为隐节点中心向量C、相应宽度向量ο、权值系数矩阵A与H、 模型偏置向量b，这些参数按如下步骤确定①用k-means聚类算法对输入数据进行聚类，得到隐节点中心c;该算法能确定最 优的聚类中心数，同时可使聚类中心合理地分布在数据空间中；②采用ρ近邻规则计算隐节点宽度 (2)其中N为隐节点中心的个数，Ci为距离第j个隐节点中心最近的ρ个隐节点中心；③采用PLS回归确定权值系数矩阵A、H和偏置向量b 根据得到的隐节点中心和宽度计算隐节点输出矩阵G，然后对输入矩阵进行扩展， 得到增广输入矩阵。对数据对{，Y}进行PLS回归，得到PLS模型参数 矩阵{T，W，P，B，Q}。为了在后面的模型更新中保留所有的信息，提取特征变量个数等于增 广输入矩阵的秩，而最终用于预测的模型所保留的特征向量个数a采用交叉校 验法确定，得到的参数矩阵记为{Ta，ffa, Pa, Ba, Qa}，由它们计算出PLS回归系数矩阵β， 从而得到Α，H和b。c.计算统计量SPE的控制限。SPE在i时刻的值是一个标量，它刻画了此时刻测 量值X(i)和y(i)对主元模型的偏离程度。当用PLS算法得到模型后，第i时刻SPE的值 为m 对0表示第i个输入样本的模型计算值，KO表示第i个输出样本的模型计算值， SPE控制限Qa按下式计算 其中 Aj = IItiII2(7) 、为保留的第i个特征向量，Ca是正态分布在检验水平为α下的临界值，Ca与 h0同号，k是主元模型中所保留的主元个数，η是全部主元个数。步骤(2)将建立的NLPLS过程监控模型应用到在线采集的实时工业过程数据，计 算新来数据的SPE值，并与控制限Qa进行比较若SPE ^ Qa，说明过程出现了异常；若SPE < Qa，说明过程正常，利用步骤(3)的方法对NLPLS过程监控模型进行更新，同时更新控制 限Qa。步骤(3)利用新数据结合原来的NLPLS过程监控模型，采用递推非线性部分最小 二乘算法对模型进行更新，同时更新控制限Qa，具体步骤如下a.记新得到的数据为Xl和Y1，采用与步骤(1)中一样的方法对新数据进行数据 预处理；b.判断是否增加新的隐节点如果新数据Xl与现有的RBF网络隐节点中心的距离大于设定值，则加入新的隐节 点；记新的隐节点中心矩阵为Cgmw，相应的宽度向量为σ gMW，对原有的隐节点中心矩阵Cg、 相应的宽度向量Qg和负荷矩阵P进行如下扩展 如果Xl与现有的RBF网络隐节点中心的距离小于等于设定值，则不需要增加隐节 点，Cg、保持不变；c.将Xl扩展为XeI = ，其中Gl为隐节点对于Xl的输出矩阵，令 ，对数据对 χΕ，γ}进行PLS回归，得到新的NLPLS过程监控模型 然后按照步骤(1)中的步骤③方法计算权值系数矩阵A、H和 偏置向量b ；d.基于新的NLPLS过程监控模型参数计算新的控制限Qa，用于新的数据，返回步 骤(2)。本专利技术提出的基于递推非线性部分最小二乘的过程监控方法，不依赖对象的解析 模型，充分利用工厂丰富的数据资源，有效的处理数据量大、数据维数高、数据共线性等问 题，弥补了传统非线性多变量统计过程监控方法不能处理时变性的不足，利用新数据更新 模型参数和控制限，从而能够适应过程的变化，提高了监控性能。具体实施例方式，具体步骤是步骤(1)基于过程数据库建立基于非线性部分最小二乘(NLPLS)过程监控模型， 同时计算控制限，具体方法是a.通过数据采集装置采集实时过程运行数据，将采集的实时过程运行数据作为数 据驱动的样本集合，表示为&(1)，7(1)}，其中乂1)表示第i组输入数据，y(i)表示第i组 输出数据；将输入数据构成输入矩阵X、将输出数据构成输出矩阵Y ；b.基于输入、输出数据建立非线性部分最小二乘模型，方法如下对矩阵X和Y进行归一化处理，使之均值为0，方差为1 ；将输入矩阵进行列扩展， 扩展项为径向基函数(RBF)神经网络的隐节点输出矩阵G和元素全为1的列向量1，其中 G的每一行对应一个输入向量作用下的隐节点的输出，隐节点的偏置项系数为1 ；对如下增 广输入矩阵和输出矩阵进行部分最小二乘(PLS)回归 得到的NLPLS过程监控模型表示为 式中，XE表示增广输入矩阵，A和H分别为对应原始输入向量和对应本文档来自技高网...

【技术保护点】
一种基于递推非线性部分最小二乘的过程监控方法，其特征在于该方法的具体步骤是：步骤（１）基于过程数据库建立基于非线性部分最小二乘过程监控模型，同时计算控制限，具体方法是：ａ．通过数据采集装置采集实时过程运行数据，将采集的实时过程运行数据作为数据驱动的样本集合，表示为｛ｘ（ｉ），ｙ（ｉ）｝，其中ｘ（ｉ）表示第ｉ组输入数据，ｙ（ｉ）表示第ｉ组输出数据；将输入数据构成输入矩阵Ｘ、将输出数据构成输出矩阵Ｙ；ｂ．基于输入、输出数据建立基于非线性部分最小二乘过程监控模型，方法如下：对矩阵Ｘｈ↓［０］（ｈ↓［０］－１）／θ↓［１］↑［２］］↑［１／ｈ↓［０］］其中：θ↓［ｉ］＝＊λ↓［ｊ］↑［ｉ］（ｉ＝１，２，３）λ↓［ｊ］＝‖ｔ↓［ｉ］‖↑［２］ｈ↓［０］＝１－２θ↓［１］θ↓［３］／３θ↓［２］↑［２］ｔ↓［ｉ］为保留的第ｉ个特征向量，Ｃ↓［α］是正态分布在检验水平为α下的临界值，Ｃ↓［α］与ｈ↓［０］同号，ｋ是主元模型中所保留的主元个数，ｎ是全部主元个数；步骤（２）将建立的基于非线性部分最小二乘过程监控模型应用到在线采集的实时工业过程数据，计算新来数据的ＳＰＥ值，并与控制限Ｑ↓［α］进行比较：若ＳＰＥ≥Ｑ↓［α］，说明过程出现了异常；若ＳＰＥ＜Ｑ↓［α］，说明过程正常，利用步骤（３）的方法对基于非线性部分最小二乘过程监控模型进行更新，同时更新控制限Ｑ↓［α］；步骤（３）利用新数据结合原来的基于非线性部分最小二乘过程监控模型，采用递推非线性部分最小二乘算法对模型进行更新，同时更新控制限Ｑ↓［α］，具体步骤如下：ａ．记新得到的数据为Ｘ１和Ｙ１，采用与步骤（１）中一样的方法对新数据进行数据预处理；ｂ．判断是否增加新的隐节点：如果新数据Ｘ１与现有的径向基函数神经网络隐节点中心的距离大于设定值，则加入新的隐节点；记新的隐节点中心矩阵为Ｃ↓［ｇｎｅｗ］，相应的宽度向量为σ↓［ｇｎｅｗ］，对原有的隐节点中心矩阵Ｃ↓［ｇ］、相应的宽度向量σ↓［ｇ］和负荷矩阵Ｐ进行如下扩展：Ｃ↓［ｇ］＝［＊＊＊］，σ↓［ｇ］＝［＊＊＊］，Ｐ＝［＊＊＊］如果Ｘ１与现有的径向基函数神经网络隐节点中心的距离小于等于设定值，Ｃ↓［ｇ］、σ↓［ｇ］和Ｐ保持不变；ｃ．将Ｘ１扩展为Ｘ↓［Ｅ］１＝［１Ｘ１Ｇ１］，其中Ｇ１为隐节点对于Ｘ１的输出矩阵，令Ｘ↓［Ｅ］＝［＊＊＊］、Ｙ＝［＊＊＊］，对数据对｛Ｘ↓［Ｅ］，Ｙ｝进行部分最小二乘回归，得到新的基于非线性部分最小二乘过程...

【技术特征摘要】

【专利技术属性】
技术研发人员：李春富，郑松，魏江，郑小青，葛铭，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：86[中国|杭州]