一种多变量化工过程的故障源和故障传播路径定位方法
【专利摘要】本发明公开了一种多变量化工过程的故障源和故障传播路径定位方法，包括：确定待监测化工过程的n个变量，并计算变量间的相关性系数；构建邻接矩阵和可达矩阵；识别源节点变量和末节点变量；删除可达矩阵的相关行和列：输出源节点变量和源节点变量，确定故障源和故障传播路径。本发明包括利用所观察变量的相关性特征，构造元素为0和1的邻接矩阵，然后通过邻接矩阵的乘方获得变量间的相互关系信息，进而构造出可达矩阵，利用可达矩阵的行列信息能够快速清晰的定位出故障信号的源头。本发明可以滑动窗口的办法，对复杂化工过程进行在线分析，对过程无干扰，及时定位并预警，能显著降低大规模化工过程的维护和故障诊断成本，提升安全性和经济效益。
【专利说明】
一种多变量化工过程的故障源和故障传播路径定位方法
技术领域
[0001] 本发明涉及石油、化工等流程工业控制系统故障诊断领域，具体是涉及一种多变 量化工过程的故障源和故障传播路径定位方法。
【背景技术】
[0002] 石油、化工等流程工业常常以多变量的单元作为被控对象。这些单元通常由几个 至几十个关键变量组成，而且这些变量间由于对象的物理特性，具备强耦合，非线性等特 点。这些特性给控制系统的设计带来困难的同时，给后期的维护也构成了挑战。
[0003] 通常的，复杂化工过程容易受到外部环境的变化，其它单元波动的影响以及本单 元的故障等因素的影响，造成变量的波动，由于耦合关系进行传播，而且由于呈现非线性、 时滞等特点使得对波动信号的分析以及对故障源的定位造难题。
[0004] 针对上述难题，学术界和工业界出现了一些解决办法，例如通过输入激励信号，辨 识被监测对象的模型变化。这一方法具有较高的可靠性，但是在工业实际应用中，激励信号 的设计是一个难题，同时人为的输入激励信号对过程本身造成干扰，这在安全性要求极高 的石油、化工过程中是很难实现的。
[0005] 另外的一些方法通过分析变量间的因果关系，从而找到故障源，这一方法的缺点 是采用了复杂的非线性时间序列分析方法，该方法需要严格的参数设定才能得到较好的效 果，从而导致该方法在实际应用中受到限制，而且不利于采用自动化的分析工具来实现。
[0006] 多变量化工过程通常伴随着强耦合性，非线性等特征，当由于传感器故障或者执 行器故障发生时，故障信息会由于系统的耦合特性而在整个系统中进行传播，危及生产质 量和生产安全。
【发明内容】
[0007] 针对现有技术存在的问题，本发明提供了一种对过程无附加干扰，通过简洁有效 的分析定位过程寻找多变量化工过程波动源头的方法，该方法无需复杂的参数设置，容易 通过计算机程序自动化实现。通过该方法的使用，能够在线分析故障源所在变量，有助于提 高系统性能维护效率，对于延长预测控制系统寿命，提高经济效益具有重要的实用价值。
[0008] -种多变量化工过程的故障源和故障传播路径定位方法，包括如下步骤：
[0009] (1)确定待监测化工过程的η个变量，η为大于1的自然数，组成变量矩阵X= [X1， X2, · · ·，Xn]，并计算Xi和Xj的相关性系数Cij，其中i, j = l，2, · · ·，n;
[0010] (2)根据置信大小，利用步骤(1)得到的相关性系数(?构建邻接矩阵A:
[0012 ] (3)利用步骤(2)中的邻接矩阵A构建可达矩阵IT;
[0013] 令p = l，进入步骤⑷；
[0014] (4)识别第p组源节点变量sP，源节点变量需要依次满足：（i)可达矩阵IT的第邱行 含有值为1的元素；（ii)同时可达矩阵的第的元素全部为0;&为1~η中的一个或多个 值；
[0015] (5)识别第ρ组末节点变量kP，末节点变量需要依次满足：（i)可达矩阵ΙΓ的第匕列 含有值为1的元素；（ii)同时可达矩阵的第匕行的元素全部为0;匕为1~η中的一个或多个 值；
[0016] (6)识别得到第ρ组源节点变量sP和末节点变量kP，同时删除可达矩阵ΙΓ的第s P行、 第sdU、第匕行、第1^列：
[0017] 如果可达矩阵ΙΓ中剩余元素不全为零，则令P = P+1，返回步骤(4);
[0018] 如果可达矩阵ΙΓ中剩余元素为零，则输出第1~p组源节点变量和末节点变量，确 定故障源和故障传播路径。
[0019] 作为优选，步骤(2)中，根据置信大小，利用步骤(1)得到的相关性系数Clj构建邻接 矩阵A的具体方法为：
[0020] 对步骤（1)所得到的相关性系数做显著性分析，令置信度高于95%的相关系数ClJ 为1，否则相关系数(?为0，构建邻接矩阵A。
[0021 ]作为优选，步骤(3)中，构建可达矩阵的方法如下：
[0022] (3-1)计算步骤(2)所得邻接矩阵的高次幂，从2次幂开始依次增加，直到第N+1次 幂，使得4^^0且4& 1 = 0,计算邻接矩阵及其所有N+1次幂的和，得到矩阵R，其表达式如下：
[0023] R=A+A2+. . ,An+1
[0024] (3-2)获得可达矩阵IT，其表达式如下：
[0025] R-= R-1。
[0026] 作为优选，步骤(6)中，故障传播路径的确定方法为：
[0027] (6-1)故障传播路径包括的总节点数为2p个，首先确定每个节点：
[0028]对于第1~p个节点，第j个节点为第j组源节点变量Sj，j为1~p的其中一个数值； [0029]对于第p+1~2p个节点，第j '个节点为第2p_( j '-1)组末节点变量k2p-(j，-υ，j '为p+ 1~2p的其中一个数值；
[0030] (6-2)以第1个节点为起始源头，将第1~2p个节点依次连接，得到故障传播路径。
[0031 ]作为优选，步骤(6)中，故障源为第1组源节点变量对应的变量:XS1。
[0032] 作为优选，所述多变量化工过程为聚丙烯生产过程。
[0033] 作为优选，所述聚丙烯生产过程采用的装置包括：
[0034] 按照物料流量依次连接的预聚单元、第一环管反应器、第二环管反应器和流化床 反应器;所述预聚单元、第一环管反应器、第二环管反应器和流化床反应器上分别设有丙烯 氢气和乙烯加料口或补加口；所述预聚单元设有催化剂加料口；所述流化床反应器上设有 聚丙烯出料口和返料口。
[0035] 本发明包括利用所观察变量的相关性特征，构造元素为0和1的邻接矩阵，然后通 过邻接矩阵的乘方获得变量间的相互关系信息，进而构造出可达矩阵，利用可达矩阵的行 列信息能够快速清晰的定位出故障信号的源头。
[0036] 本发明与现有技术相比具有的有益效果：
[0037] 1、无需外部信号激励，故对系统的没有附加扰动，能够实现完全的非侵入式检测 与诊断。
[0038] 2、计算简单，便于操作，无需复杂的算法，易于在现有的DCS工作站或者控制系统 上位机上实施。
[0039] 3、能够在故障发生早期定位到故障源变量，从而能够起到预防作用，有效增强复 杂系统的安全性和无故障运行时间。
[0040] 4、完全采用数据驱动的方法，无需具备过程耦合机理及动态特性知识。
[0041] 本发明可以滑动窗口的办法，对复杂化工过程进行在线分析，对过程无干扰，及时 定位并预警，能显著降低大规模化工过程的维护和故障诊断成本，提升安全性和经济效益。
【附图说明】
[0042] 图1为本发明双环管工艺聚丙烯(PP)生产过程示意图。
[0043]图2为实施例中得出的故障传播路径图。
【具体实施方式】
[0044]结合实施例对本发明作进一步说明：
[0045] -种多变量化工过程的故障源和故障传播路径定位方法，它的步骤如下：
[0046] 1)记所要分析的多变量化工过程具有η个待分析变量，变量名记为xi，其中i = 1， 2, . . .，n。变量数据用矩阵方式表示为变量矩阵X= [X1，X2, . . .，χη]，并计算Xi和Xj的相关性 系数cij，其中i，j = l，2,. . .，n;
[0047] 2)根据步骤1)的计算结果构造邻接矩阵，方法如下：
[0048]对步骤1)所得到的相关性系数做显著性分析，令置信度高于95%的相关系数ClJ为 1，否则为〇,并按照如下方式构造邻接矩阵A，
[0050] 3)构造可达矩阵，其方法如下：
[0051] 计算步骤2)所得邻接矩阵的高次幂，从2次幂开始依次增加，直到第N+1次幂，使得 AV(mAN+1 = 0,计算邻接矩阵及其所有N+1次幂的和，得到矩阵R，其表达式如下：
[0052] R=A+A2+---An+1
[0053] 进一步的，通过矩阵IT获得可达矩阵，其表达式如下：
[0054] IT = R-1(即矩阵R中所有元素减去1)。
[0055] (4)识别第p组源节点变量sP，源节点变量需要依次满足：（i)可达矩阵ΙΓ的第&行 含有值为1的元素；（ii)同时可达矩阵的第的元素全部为0;&为1~η中的一个或多个 值；
[0056] (5)识别第ρ组末节点变量kP，末节点变量需要依次满足：（i)可达矩阵ΙΓ的第匕列 含有值为1的元素；（ii)同时可达矩阵的第匕行的元素全部为0;匕为1~η中的一个或多个 值；
[0057] (6)识别得到第p组源节点变量sP和末节点变量kP，同时删除可达矩阵IT的第sP行、 第sdU、第匕行、第1^列：
[0058] 如果可达矩阵ΙΓ中剩余元素不全为零，则令p = p+l，返回步骤(4);
[0059] 如果可达矩阵ΙΓ中剩余元素为零，则输出第1~p组源节点变量和末节点变量，确 定故障源和故障传播路径。
[0000]其中，上述步骤(6)中，故障传播路径的确定方法为：
[0061 ] (6-1)故障传播路径包括的总节点数为2p个，首先确定每个节点：
[0062 ]对于第1~p个节点，第j个节点为第j组源节点变量Sj，j为1~p的其中一个数值； [0063]对于第p+1~2p个节点，第j '个节点为第2p_( j '-1)组末节点变量，j '为p+ 1~2p的其中一个数值；
[0064] (6-2)以第1个节点为起始源头，将第1~2p个节点依次连接，得到故障传播路径。
[0065] 步骤(6)中，故障源为第1组源节点变量对应的变量:XS1。
[0066] 实施例
[0067] 某石化公司的聚丙稀化〇17口1'(^716116，？？)生产过程工艺流程如图1所示，它采用 Spheripol工艺，使用Ziegler-Natta系列催化剂。整个生产装置主要由两个环管反应器和 一个流化床反应器组成，其中均聚反应在环管内进行而共聚反应在流化床反应器中进行。 催化剂、丙烯、乙烯和氢气经预聚单元R200处理后进入第一环管反应器R201，在环管中循环 并发生均聚反应。第一环管反应器R201的部分混合物(包括反应物和生产物)和新鲜进料气 体一起进入第二环管反应器R202继续反应。均聚反应的产物进入后续的流化床反应器，与 乙烯、丙烯等发生共聚反应。在这套装置中通过改变操作条件能够生产不同牌号的聚丙烯 产品。
[0068]熔融指数MI是表征聚合物产品质量特性一个重要指标，通常用MI来区分不同的聚 丙烯产品牌号。对于丙烯聚合这一复杂的化工生产过程，选取10个主要的过程变量进行分 析，如表1所示。
[0069] 表1过程变量描述
[0071 ]根据步骤2)提供的方法，得到邻接矩阵如下：
[0074]根据步骤3)构造可达矩阵，通过计算，得到邻接矩阵的5次方为全零矩阵，则按照 步骤3)的方法得到可达矩阵如下：
[0076] 根据步骤4)可得到故障传动的源变量为变量1，即R201的密度；
[0077] 根据步骤5)可得到故障传动的末变量为变量3,4，10;
[0078] 根据步骤6)，删除变量1，3,4，10所对应的行和列后得到矩阵如下：
[0081 ]重复步骤4)，步骤5)进一步得到故障的源变量为变量5,末变量为6);
[0082]根据步骤6,删除变量5)和变量6)对应的行和列后，得到如下矩阵：
[0084] 矩阵的元素全部为0,则分析步骤结束，通过以上的分析结果，可以得到该系统的 故障传播路径如图2所示。由图2可以明确看出该聚丙稀(pol ypropy 1 ene，PP)生产过程的故 障传播途径，其中故障源为幻，即R201的密度可能是导致故障的最初源头。
[0085] 在化工生产过程中，可采用本发明的方法进行实时监测，并通过监测结果利用现 有方法或经验进行反馈控制。或者，也可将监测的故障结果储存或者直接显示在现有的显 示设备上，供现场人员参考和进行过程控制。
[0086] 上述实施例中，得到的第1组源节点为一个，故障源路径相对较为简单，针对源节 点和末节点为多个的情况，可采用现有的统计分析和实验方法，针对节点与节点之间的影 响关系进行分析，最终得到具体的故障传播途径。利用本发明的方法，可以从整体上给出故 障传播路径，能够大幅度降低故障定位的错误率和提高故障定位效率。
【主权项】
1. 一种多变量化工过程的故障源和故障传播路径定位方法，其特征在于，包括如下步 骤： (1)确定待监测化工过程的η个变量，η为大于1的自然数，组成变量矩阵X=[X1，X2，...， Xn]，并计算^和义』的相关性系数Cij，其中i，j = l，2,…，n; ⑵根据置信大小，利用步骤⑴得到的相关性系数(?构建邻接矩阵A:1其中元素 aij为0或1; (3) 利用步骤(2)中的邻接矩阵A构建可达矩阵IT; 令P=l，进入步骤(4); (4) 识别第p组源节点变量sP，源节点变量需要依次满足：（i)可达矩阵IT的第sP行含有 值为1的元素；（ii)同时可达矩阵的第sdU的元素全部为0;&为1~η中的一个或多个值； (5) 识别第ρ组末节点变量kP，末节点变量需要依次满足：（i)可达矩阵ΙΓ的第kP列含有 值为1的元素；（ii)同时可达矩阵的第1^行的元素全部为0;匕为1~η中的一个或多个值； (6) 识别得到第ρ组源节点变量sP和末节点变量kP，同时删除可达矩阵ΙΓ的第&行、第s P 列、第Μ?、第匕列： 如果可达矩阵R_中剩余元素不全为零，则令Ρ=Ρ+1，返回步骤(4); 如果可达矩阵ΙΓ中剩余元素为零，则输出第1~Ρ组源节点变量和末节点变量，确定故障 源和故障传播路径。2. 根据权利要求1所述的多变量化工过程的故障源和故障传播路径定位方法，其特征 在于，步骤(2)中，根据置信大小，利用步骤（1)得到的相关性系数 Clj构建邻接矩阵Α的具体 方法为： 对步骤（1)所得到的相关性系数做显著性分析，令置信度高于95%的相关系数Clj为1， 否则相关系数(?为0，构建邻接矩阵A。3. 根据权利要求1所述的多变量化工过程的故障源和故障传播路径定位方法，其特征 在于，步骤(3)中，构建可达矩阵的方法如下： (3-1)计算步骤(2)所得邻接矩阵的高次幂，从2次幂开始依次增加，直到第N+1次幂，使 得AV(mAN+1 = 0,计算邻接矩阵及其所有N+1次幂的和，得到矩阵R，其表达式如下： R=A+A2+---AN1 (3-2)获得可达矩阵IT，其表达式如下： R-=R-l〇4. 根据权利要求1所述的多变量化工过程的故障源和故障传播路径定位方法，其特征 在于，步骤(6)中，故障传播路径的确定方法为： (6-1)故障传播路径包括的总节点数为2p个，首先确定每个节点： 对于第1~P个节点，第j个节点为第j组源节点变量sj，j为1~p的其中一个数值； 对于第P+1~2p个节点，第j'个节点为第2p-( j'-l)组末节点变量k2P-(j，-i)，j '为p+1~ 2p中的一个数值； (6-2)以第1个节点为起始源头，将第1~2p个节点依次连接，得到故障传播路径。5. 根据权利要求1所述的多变量化工过程的故障源和故障传播路径定位方法，其特征 在于，步骤(6)中，故障源为第1组源节点变量对应的变量:xsj?6. 根据权利要求5所述的多变量化工过程的故障源和故障传播路径定位方法，其特征 在于，所述多变量化工过程为聚丙烯生产过程。
【文档编号】G05B23/02GK105867345SQ201610171423
【公开日】2016年8月17日
【申请日】2016年3月24日
【发明人】费正顺, 陈贵, 刘康玲
【申请人】浙江科技学院