CN111915120B - 一种基于动态正交近邻保持嵌入模型的过程监测方法 - Google Patents
一种基于动态正交近邻保持嵌入模型的过程监测方法 Download PDFInfo
- Publication number
- CN111915120B CN111915120B CN201910474174.4A CN201910474174A CN111915120B CN 111915120 B CN111915120 B CN 111915120B CN 201910474174 A CN201910474174 A CN 201910474174A CN 111915120 B CN111915120 B CN 111915120B
- Authority
- CN
- China
- Prior art keywords
- vector
- matrix
- formula
- new
- neighbor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 91
- 230000008569 process Effects 0.000 title claims abstract description 52
- 238000012544 monitoring process Methods 0.000 title claims abstract description 42
- 239000013598 vector Substances 0.000 claims abstract description 102
- 230000009466 transformation Effects 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims description 59
- 238000005070 sampling Methods 0.000 claims description 17
- 238000010606 normalization Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 238000004519 manufacturing process Methods 0.000 claims description 8
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000012549 training Methods 0.000 abstract description 6
- 238000000605 extraction Methods 0.000 abstract description 2
- 238000005065 mining Methods 0.000 abstract 1
- 239000000463 material Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 239000000498 cooling water Substances 0.000 description 4
- 239000007788 liquid Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 239000007791 liquid phase Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 241000540325 Prays epsilon Species 0.000 description 1
- 230000003416 augmentation Effects 0.000 description 1
- 238000012824 chemical production Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011112 process operation Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于动态正交近邻保持嵌入模型的过程监测方法,该方法在挖掘训练数据中潜藏的有用特征时,能同时考虑自相关特征与局部近邻结构特征并保证投影变换向量的正交特性。本发明方法首先将自相关性与局部近邻结构嵌入量化成一个目标函数;其次,在保证各投影变换向量相互正交的约束下求解该目标函数,从而利用得到的投影变换向量实施特征提取;最后,利用提取出的潜在特征以及模型误差实施过程监测。与传统方法相比:首先本发明方法中涉及的动态正交近邻保持嵌入模型的建立方法是一种全新的建模策略,它同时考虑了自相关性特征与局部近邻特征,并且保证了投影变换向量的正交特性,能够更全面地挖掘训练数据中潜藏的有用信息。
Description
技术领域
本发明涉及一种数据驱动的过程监测方法,尤其涉及一种基于动态正交近邻保持嵌入模型的过程监测方法。
背景技术
为保证生产过程的安全运行与维持产品质量的稳定性,学术界与工业界都投入了大量的人力与物力研究以故障监测为核心任务的过程监测方法。而今,工业发展已进入信息化建设,利用采样数据实施数据驱动的过程运行状态监测,已经成为工业自动化领域的研究热门之一。通常来讲,数据驱动的过程监测方法的核心思想在于:如何对过程正常数据进行有效地挖掘以提取能反应过程运行状态的潜在有用信息。然而,考虑到现代工业过程规模的复杂化趋势,采集到的工业数据所呈现出的特征往往也是非常复杂的。可以说,如何更有效地挖掘出过程数据中潜藏的有用信息,并建立更适于监测现代工业过程对象的过程监测模型,一直以来都是该研究领域所面临的主要问题。
在现有的科研文献与专利材料中,主元分析(Principal Component Analysis,缩写:PCA)与近邻保持嵌入(Neighborhood Preserving Embedding,缩写:NPE)算法等都被用于过程监测。与PCA算法提取数据方差信息不同的是,NPE方法则考虑的是原始数据点的局部近邻特征。NPE在提取原始数据中的潜在信息时,尽量保留数据点在空间距离上的分布特征。此外,由于现代工业过程的采样时间都较短,采样数据之间的时序自相关性是数据驱动的过程监测方法在建模时另一个必须考虑的一个问题。最近,有研究学者提出一类能挖掘自相关潜在特征的数据建模算法,并较好地解决了动态过程监测问题。此外,通过增广矩阵引入延时测量值可将数据的时序自相关性考虑进来,是一种最为经典的动态过程监测思路,PCA或NPE等都可以直接拓展成动态PCA与动态NPE方法。
由于局部近邻特征与自相关特征都是数据中潜藏的有用特征,因此两者都应予以考虑。然而,在现有的科研文献与专利材料中,却没有同时提取局部近邻特征与自相关特征的算法。此外,提取潜在特征的投影变换向量的正交性在建立模型时也是非常重要的一个考量,因为这直接关系着模型误差的计算。因此,有效地考虑局部近邻特征与自相关特征是动态过程监测领域继续解决的一个问题。
发明内容
本发明所要解决的主要技术问题是:如何在挖掘训练数据中潜藏的有用特征时,同时考虑自相关特征与局部近邻结构特征并保证投影变换向量的正交特性,并基于此实施对对生产过程运行状态的监测。具体来讲,本发明方法首先将自相关性与局部近邻结构嵌入量化成一个目标函数;其次,在保证各投影变换向量相互正交的约束下求解该目标函数,从而利用得到的投影变换向量实施特征提取;最后,提取出的潜在特征以及模型误差实施过程监测模。
本发明解决上述技术问题所采用的技术方案为:一种基于动态正交近邻保持嵌入模型的过程监测方法,包括以下步骤:
步骤(1):在生产过程正常运行状态下,按照采样时间先后依次采集n个样本数据组成矩阵X∈Rn×m,并计算矩阵X中各行向量的均值向量μ与标准差向量δ,其中m为测量变量的总个数、R为实数集,Rn×m表示n×m维的实数矩阵。
步骤(2):利用均值向量μ与标准差向量δ对矩阵X中各行向量实施标准化处理,从而得到矩阵其中/>为标准化处理后的数据向量,i=1,2,…,n,上标号T表示矩阵或向量的转置。
步骤(3):设置自相关阶数为d,将矩阵中第d+1行至第n行向量组成矩阵Y∈R(n -d)×m,并根据公式①构造矩阵Z:
一方面,为提取训练数据中自相关特征,可根据最大化潜在特征的时序自相关性来设计如下所示的目标函数:
上式中,表示Kronecker乘法,p∈Rm×1表示投影变换向量,β∈Rd×1为自相关系数向量,d为自相关阶数(一般可取d=2)。由此可见,上式②要求经投影变换向量p∈Rm×1转换后的潜在特征满足时序自相关性最大的要求。
另一方面,为保留训练数据中局部近邻结构特征,可借用如下所示的NPE算法的目标函数,从而在投影变换的过程中嵌入局部近邻重构信息:
上式中,矩阵M=(In-W)(In-W)T,In表示n×n维的单位矩阵,M∈Rn×n为近邻保持嵌入矩阵。根据NPE算法推理过程,现将近邻保持嵌入矩阵M∈Rn×n的获取过程介绍如下:
步骤(A):针对矩阵中的各个数据向量/>计算各数据向量与其他n-1个数据向量之间的距离/>其中j=1,2,…,n且j≠i。
步骤(B):根据这些计算出来的距离,为各数据向量找出距离其最近的k个数据向量,从而组成近邻矩阵/>其中k为近邻个数,一般可取k=5。
步骤(C):先根据公式计算向量wi∈Rk×1,再根据wi=wi/||wi||对向量wi进行归一化处理。
步骤(D):根据近邻矩阵中的k个数据向量对应于矩阵/>中的位置,相应地将向量wi中的k个元素赋值给全零矩阵W∈Rn×n中第i行的相应元素。
步骤(E):根据公式M=(In-W)(In-W)T计算近邻保持嵌入矩阵M∈Rn×n。
若想同时考虑自相关特征与局部近邻结构特征,需在求解投影变换向量时,同时满足公式②与公式③中的目标函数。由于J1是最大化问题而J2是最小化问题,因此可设计如下所示的目标函数:
为不失一般性,可设此外,为保证投影变换向量的正交化,还需增加正交约束条件:pa Tp1=pa Tp2=…=pa Tpa-1=0,其中a表示求解的第a个投影变换向量。上式④最终转变成如下所示的优化问题:
再通过拉格朗日乘子法即可求解上式⑤:构造如下所示的拉格朗日函数L:
计算函数L相对于pa与βa的偏微分:
上式中,Im与Id分别表示m×m维与d×d维的单位矩阵。根据极值求解思路,上式⑦与⑧中的偏微分都等于零。因此,可推理出:
其中,显而易见,公式⑩则定义了βa与pa之间的等式关系。在公式⑨中等号两边同时左乘pa T,即可得到:/>所以,拉格朗日乘子λ即等价于公式⑤中的目标函数值。
在公式⑨中等号两边同时依次左乘则可得到a-1个等式关系:
上式则可以等价的写成如下所示的矩阵形式:
上式中,γ(a-1)=[γ1,γ2,…,γa-1]T,A(a-1)=[p1,p2,…,pa-1],/>将公式⑨等号两边同时左乘/>后再减去公式/>即可得到:
因此,公式定义了一个常规的特征值问题,投影变换向量pa即是最大特征值所对应的特征同量。
由于矩阵G的计算涉及到自回归系数向量βa,而βa的求解需要在pa已知的前提下根据公式⑩计算得到。所以,pa与βa的求解相互耦合,可通过迭代收敛的方式予以同时求解。
综上所述,本发明方法的步骤(4)即是在求解得到A个投影变换向量p1,p2,…,pA后,建立动态正交近邻保持嵌入模型。
步骤(4):设置投影变换向量的个数为A,根据如下所示步骤(4.1)至步骤(4.6)求解得到A个投影变换向量p1,p2,…,pA后,建立动态正交近邻保持嵌入模型:其中A<m。
步骤(4.1):初始化a=1与初始化pa为任意m×1维的非零实数向量。
步骤(4.2):根据公式计算自相关系数向量βa,并根据公式βa=βa/||βa||对βa实施归一化处理,其中|| ||表示计算向量的长度。
步骤(4.3):根据公式计算矩阵G,判断是否满足条件:a<2?若是,则计算特征值问题:/>最大特征值λ所对应的特征向量pa;若否,则求解如公式/>所示特征值问题中最大特征值所对应的特征向量pa:
步骤(4.4):根据公式pa=pa/||pa||对pa实施归一化处理后,判断pa是否收敛?收敛的标准为:pa中各元素不再发生变化,若是,则得到第a个投影变换向量pa并执行步骤(4.5);若否,则返回步骤(4.2)。
步骤(4.5):判断是否满足条件:a<A?若是,则置a=a+1后初始化pa为任意m×1维的非零实数向量,并返回步骤(4.2);若否,则得到A个投影变换向量p1,p2,…pA。
步骤(4.6):根据公式分别计算得分矩阵S与误差矩阵E,那么动态正交近邻保持嵌入模型即为:/>其中P=[p1,p2,…pA]T。
步骤(5):根据公式Λ=STS/(n-1)计算协方差矩阵Λ,并根据公式ψ=diag{SΛ- 1ST}与分别计算监测指标向量ψ与Q,其中diag{ }表示将矩阵对角线上的元素组成向量的操作。
步骤(6):分别将监测指标向量ψ与Q中第n/100个最大元素记做监测指标上限ψlim与Qlim。
以上步骤(1)至步骤(6)为本发明方法的离线建模阶段,离线建模完成后,即可实施对工业过程过程对象的在线监测。
步骤(7):采集最新采样时刻的样本数据xnew∈Rm×1,并利用均值向量μ与标准差向量δ对xnew实施标准化处理得到向量
步骤(8):根据公式计算得分向量snew后,再根据公式ψnew=snewΛ-1snew T与/>计算出监测指标ψnew与Qnew的具体数值。
步骤(9):判断是否满足条件:ψnew≤ψlim且Qnew≤Qlim?若是,则当前采样时刻过程正常运行,返回步骤(7)继续实施对下一采样时刻数据的监测;若否,则当前采样时刻过程进入异常工况状态。
与传统方法相比,本发明方法的优势在于:
首先,本发明方法中涉及的动态正交近邻保持嵌入模型的建立方法是一种全新的建模策略,它同时考虑了自相关性特征与局部近邻特征,并且保证了投影变换向量的正交特性,能够更全面地挖掘训练数据中潜藏的有用信息。其次,在具体实施案例中,通过TE过程上的监测详情对比,验证了本发明方法相比于传统动态过程监测方法在故障监测上的优越性。因此,本发明方法是一种更为优选的动态过程监测方法。
附图说明
图1为本发明方法的实施流程图。
图2为建立动态正交近邻保持嵌入模型的实施流程图。
图3本发明方法与传统方法的故障监测效果对比图。
具体实施方式
下面结合附图与具体的实施案例对本发明方法进行详细的说明。
如图1所示,本发明公开一种基于动态正交近邻保持嵌入模型的过程监测方法。现结合一个具体的实施案例来阐述本发明方法的具体实施方式。
所测试的过程对象为TE过程,该过程原型是伊斯曼化工生产车间的一个实际工艺流程。目前,TE过程因其流程的复杂性,已作为一个标准实验平台被广泛用于故障检测研究。整个TE过程包括22个测量变量、12个操作变量、和19个成分测量变量。所采集的数据分为22组,其中包括1组正常工况下的数据集与21组故障数据。而在这些故障数据中,有16个是已知故障类型,如冷却水入口温度或进料成分的变化、阀门粘滞、反应动力学漂移等,还有5个故障类型是未知的。为了对该过程进行监测,选取如表1所示的33个过程变量,接下来结合该TE过程对本发明具体实施步骤进行详细的阐述。
表1:TE过程监测变量。
序号 | 变量描述 | 序号 | 变量描述 | 序号 | 变量描述 |
1 | 物料A流量 | 12 | 分离器液位 | 23 | D进料阀门位置 |
2 | 物料D流量 | 13 | 分离器压力 | 24 | E进料阀门位置 |
3 | 物料E流量 | 14 | 分离器塔底流量 | 25 | A进料阀门位置 |
4 | 总进料流量 | 15 | 汽提塔液位 | 26 | A和C进料阀门位置 |
5 | 循环流量 | 16 | 汽提塔压力 | 27 | 压缩机循环阀门位置 |
6 | 反应器进料 | 17 | 汽提塔底部流量 | 28 | 排空阀门位置 |
7 | 反应器压力 | 18 | 汽提塔温度 | 29 | 分离器液相阀门位置 |
8 | 反应器液位 | 19 | 汽提塔上部蒸汽 | 30 | 汽提塔液相阀门位置 |
9 | 反应器温度 | 20 | 压缩机功率 | 31 | 汽提塔蒸汽阀门位置 |
10 | 排空速率 | 21 | 反应器冷却水出口温度 | 32 | 反应器冷凝水流量 |
11 | 分离器温度 | 22 | 分离器冷却水出口温度 | 33 | 冷凝器冷却水流量 |
步骤(1):在TE过程正常运行状态下,按照采样时间先后依次采集n=960个样本数据x1,x2,…,xn组成矩阵X=[x1,x2,…,xn]T∈R960×33,并计算矩阵X中各行向量的均值向量μ与标准差向量δ。
步骤(2):利用均值向量μ与标准差向量δ对矩阵X中各行向量实施标准化处理,从而得到矩阵其中/>为标准化处理后的数据向量。
步骤(3):设置自相关阶数为d,将矩阵中第d+1行至第n行向量组成矩阵Y∈R(n -d)×m,并根据公式①构造矩阵Z。
步骤(4):根据上述步骤(A)至步骤(E)计算得到近邻保持嵌入矩阵M∈Rn×n。
步骤(5):设置投影变换向量的个数为A=12,根据如图2所示的实施流程求解得到A个投影变换向量p1,p2,…,pA后,建立动态正交近邻保持嵌入模型:
步骤(6):根据公式Λ=STS/(n-1)计算协方差矩阵Λ,并根据公式ψ=diag{SΛ- 1ST}与分别计算监测指标向量ψ与Q。
步骤(7):分别将监测指标向量ψ与Q中第n/100≈10个最大元素记做监测指标上限ψlim与Qlim。
离线建模阶段完成后,即可实施在线过程监测。采集TE过程在第19种故障工况条件下的测试数据,其中前160个样本数据采集自正常运行状态,故障样本数据从第161个采样时刻引入。
步骤(8):采集最新采样时刻的样本数据xnew∈R33×1,并利用均值向量μ与标准差向量δ对xnew实施标准化处理得到向量
步骤(9):根据公式计算得分向量snew后,再根据公式ψnew=snewΛ-1snew T与/>计算出监测指标ψnew与Qnew的具体数值。
步骤(10):判断是否满足条件:ψnew≤ψlim且Qnew≤Qlim?若是,则当前采样时刻过程正常运行,返回步骤(7)继续实施对下一采样时刻数据的监测;若否,则当前采样时刻过程进入异常工况状态。
将本发明方法与其他传统过程监测方法(包括NPE与动态PCA)对TE过程该故障的监测详情对比于图3中。从图3中可以明显地发现,本发明方法的故障漏报情况更少,监测效果得到了有效的改善。
上述实施案例只用来解释说明本发明的具体实施,而不是对本发明进行限制。在本发明的精神和权利要求的保护范围内,对本发明做出的任何修改,都落入本发明的保护范围。
Claims (1)
1.一种基于动态正交近邻保持嵌入模型的过程监测方法,其特征在于,包括以下步骤:
首先,离线建模阶段包括如下所示步骤(1)至步骤(5);
步骤(1):在生产过程正常运行状态下,按照采样时间先后依次采集n个样本数据组成矩阵X∈Rn×m,并计算矩阵X中各行向量的均值向量μ与标准差向量δ,其中m为测量变量的总个数、R为实数集,Rn×m表示n×m维的实数矩阵;
步骤(2):利用均值向量μ与标准差向量δ对矩阵X中各行向量实施标准化处理,从而得到矩阵其中/>为标准化处理后的数据向量,i=1,2,…,n,上标号T表示矩阵或向量的转置;
步骤(3):设置自相关阶数为d,将矩阵中第d+1行至第n行向量组成矩阵Y∈R(n-d)×m,并根据公式①构造矩阵Z:
步骤(4):根据如下所示步骤(4.1)至步骤(4.5)计算得到近邻保持嵌入矩阵M∈Rn×n:
步骤(4.1):针对矩阵中的各个数据向量/>计算各数据向量与其他n-1个数据向量之间的距离/>其中j=1,2,…,n且j≠i;
步骤(4.2):根据这些计算出来的距离,为各数据向量找出距离其最近的k个数据向量,从而组成近邻矩阵/>
步骤(4.3):先根据公式计算向量wi∈Rk×1,再根据wi=wi/||wi||对向量wi进行归一化处理;
步骤(4.4):根据近邻矩阵的k个数据向量对应于矩阵/>中的位置,相应地将向量wi中的k个元素赋值给全零矩阵W∈Rn×n中第i行的相应元素;
步骤(4.5):根据公式M=(In-W)(In-W)T计算近邻保持嵌入矩阵M∈Rn×n;
步骤(5):设置投影变换向量的个数为A,根据如下所示步骤(5.1)至步骤(5.6)求解得到A个投影变换向量p1,p2,…,pA后,建立动态正交近邻保持嵌入模型:
步骤(5.1):初始化a=1与初始化pa为任意m×1维的非零实数向量;
步骤(5.2):根据公式计算自相关系数向量βa,并根据公式βa=βa/||βa||对βa实施归一化处理,其中|| ||表示计算向量的长度,/>表示Kronecker乘法;
步骤(5.3):根据公式计算矩阵G,判断是否满足条件:a<2;若是,则计算特征值问题:/>最大特征值λ所对应的特征向量pa;若否,则求解如公式②所示特征值问题中最大特征值所对应的特征向量pa:
其中,A(a-1)=[p1,p2,…,pa-1],Im与Id分别表示m×m维与d×d维的单位矩阵;
步骤(5.4):根据公式pa=pa/||pa||对pa实施归一化处理后,判断pa是否收敛;若是,则得到第a个投影变换向量pa并执行步骤(5.5);若否,则返回步骤(5.2);
步骤(5.5):判断是否满足条件:a<A;若是,则置a=a+1后初始化pa为任意m×1维的非零实数向量,并返回步骤(5.2);若否,则得到A个投影变换向量p1,p2,…pA;
步骤(5.6):根据公式与/>分别计算得分矩阵S与误差矩阵E,那么动态正交近邻保持嵌入模型即为:/>其中P=[p1,p2,…pA]T;
步骤(6):根据公式Λ=STS/(n-1)计算协方差矩阵Λ,并根据公式ψ=diag{SΛ-1ST}与分别计算监测指标向量ψ与Q,其中diag{ }表示将矩阵对角线上的元素组成向量的操作;
步骤(7):分别将监测指标向量ψ与Q中第n/100个最大元素记做监测指标上限ψlim与Qlim;
其次,实施在线过程监测的步骤如下所示:
步骤(8):采集最新采样时刻的样本数据xnew∈Rm×1,并利用均值向量μ与标准差向量δ对xnew实施标准化处理得到向量
步骤(9):根据公式计算得分向量snew后,再根据公式ψnew=snewΔ-1snew T与计算出监测指标ψnew与Qnew的具体数值;
步骤(10):判断是否满足条件:ψnew≤ψlim且Qnew≤Qlim;若是,则当前采样时刻过程正常运行,返回步骤(8)继续实施对下一采样时刻数据的监测;若否,则当前采样时刻过程进入异常工况状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910474174.4A CN111915120B (zh) | 2019-05-24 | 2019-05-24 | 一种基于动态正交近邻保持嵌入模型的过程监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910474174.4A CN111915120B (zh) | 2019-05-24 | 2019-05-24 | 一种基于动态正交近邻保持嵌入模型的过程监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111915120A CN111915120A (zh) | 2020-11-10 |
CN111915120B true CN111915120B (zh) | 2023-09-26 |
Family
ID=73242711
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910474174.4A Active CN111915120B (zh) | 2019-05-24 | 2019-05-24 | 一种基于动态正交近邻保持嵌入模型的过程监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111915120B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106384130A (zh) * | 2016-09-22 | 2017-02-08 | 宁波大学 | 基于数据多近邻局部特征嵌入的故障检测方法 |
CN106404442A (zh) * | 2016-09-22 | 2017-02-15 | 宁波大学 | 基于数据邻域特征与非邻域特征保持的工业过程故障检测方法 |
-
2019
- 2019-05-24 CN CN201910474174.4A patent/CN111915120B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106384130A (zh) * | 2016-09-22 | 2017-02-08 | 宁波大学 | 基于数据多近邻局部特征嵌入的故障检测方法 |
CN106404442A (zh) * | 2016-09-22 | 2017-02-15 | 宁波大学 | 基于数据邻域特征与非邻域特征保持的工业过程故障检测方法 |
Non-Patent Citations (4)
Title |
---|
Statistical process monitoring based on nonlocal and multiple neighborhoods preserving embedding model;Chudong Tong 等;Journal of Process Control;第65卷;第34–40页 * |
变量加权型主元分析算法及其在故障检测中的应用;蓝艇 等;化工学报;第68卷(第8期);第3177-3182页 * |
基于潜变量自回归算法的化工过程动态监测方法;唐俊苗 等;化工学报;第70卷(第3期);第987-994页 * |
基于聚类选择k近邻的LLE 算法及故障检测;薄翠梅 等;化工学报;第67卷(第3期);第925-930页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111915120A (zh) | 2020-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108803520B (zh) | 一种基于变量非线性自相关性剔除的动态过程监测方法 | |
CN111913460B (zh) | 一种基于序列相关局部保持投影算法的故障监测方法 | |
CN109522972B (zh) | 一种基于潜变量自回归模型的动态过程监测方法 | |
CN108897286B (zh) | 一种基于分散式非线性动态关系模型的故障检测方法 | |
CN108508865B (zh) | 一种基于分散式osc-pls回归模型的故障检测方法 | |
CN108469805B (zh) | 一种基于动态性最优选择的分散式动态过程监测方法 | |
CN106940808A (zh) | 一种基于改进型主元分析模型的故障检测方法 | |
CN109669415B (zh) | 一种基于结构化典型变量分析的动态过程监测方法 | |
CN108445867B (zh) | 一种基于分散式icr模型的非高斯过程监测方法 | |
CN114757269A (zh) | 一种基于局部子空间-邻域保持嵌入的复杂过程精细化故障检测方法 | |
CN108491878B (zh) | 一种基于多个误差生成模型的故障分类诊断方法 | |
CN108572639B (zh) | 一种基于主成分自相关性剔除的动态过程监测方法 | |
CN108427398B (zh) | 一种基于分散式ar-pls模型的动态过程监测方法 | |
CN111915121B (zh) | 一种基于广义典型变量分析的化工过程故障检测方法 | |
CN111913415B (zh) | 一种基于时序数据分析的连续搅拌反应釜运行状态监测方法 | |
CN111914886B (zh) | 一种基于在线简略核学习的非线性化工过程监测方法 | |
CN111915120B (zh) | 一种基于动态正交近邻保持嵌入模型的过程监测方法 | |
CN111914206B (zh) | 一种基于动态近邻保持嵌入算法的过程监测方法 | |
CN111914384B (zh) | 一种基于自相关潜变量模型的化工过程状态监测方法 | |
CN108536943B (zh) | 一种基于多生产单元变量交叉相关解耦策略的故障监测方法 | |
CN109669412B (zh) | 一种基于新型动态独立成分分析的非高斯过程监测方法 | |
CN111913461B (zh) | 一种基于正则化gcca模型的分布式化工过程监测方法 | |
CN111694328B (zh) | 一种基于多块独立成分分析算法的分布式过程监测方法 | |
CN114707424A (zh) | 基于质量相关慢特征分析算法的化工过程软测量方法 | |
CN111913462B (zh) | 一种基于广义多块独立元分析模型的化工故障监测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |