流程工业的厂级控制回路振荡源定位方法
流程工业的厂级控制回路振荡源定位方法
【技术领域】
[0001]本发明属于控制或调节系统领域，尤其是一种基于模型识别方法的流程工业厂级控制回路振荡源定位方法。
【背景技术】
[0002]由于现代过程工业日益复杂及控制器日常维护不足等方面的原因，电力、石油、化工、钢铁、造纸、水泥等工业过程普遍存在控制回路性能差的现象，国外有研宄表明，工业过程中有1/3的控制回路处于振荡状态中。回路振荡是一种明显的控制性能退化的表现形式，有多种可能的原因会引起控制回路振荡:控制阀的非线性故障，如摩擦、死区、反冲及粘滞等现象；控制参数调整不适当，使得系统趋于不稳而引起振荡；来自外部振荡干扰和回路自身的反馈机制，这种情况通常在控制回路间有较强耦合关联的系统中见到，一个振荡源会通过耦合路径传递到其它回路中去；高度的非线性过程；在用的MPC控制器与对象模型不匹配。控制回路中振荡的存在将导致产品质量下降、能耗增加、设备磨损加快等一系列问题。尤其在现代企业规模日益扩大的情况下，控制回路数目成千上万，一旦由于某种异常因素(例如阀门出现了粘滞或外部扰动等)导致了在这样大的控制系统中的某一位置上发生了振荡，由于过程系统中的各个处理装置间存在着复杂的耦合，则这种出现在局部位置上的振荡将被快速传递扩散开来，产生所谓的“厂级控制回路振荡”现象。一旦“厂级控制回路振荡”出现，不但振荡的各种危害将被成倍的放大，而且还造成了原始的那个振荡源所在的控制回路被淹没在许许多多具有非常相似振荡特征的控制回路当中。也就是说，一旦“厂级控制回路振荡”出现，
[0003]如果没有一种行之有效的定位厂级控制回路振荡源的方法，将导致“厂级控制回路振荡”现象将在一个相当长的时间内存在。显然，要消除“厂级控制回路振荡”所引起的这种严重危害，首先要找到其来源，然后对其定性，最后根据不同的原因采取相应的措施，加以解决。值得一提的是，发电机组等大型过程装置的运行过程常常伴随十分苛刻的生产条件或环境，如高温、高压，因而运行的安全性十分重要。由此可见，定位厂级控制回路振荡源方法的研宄对于电力、石油、化工、钢铁、造纸、水泥等工业过程的安全、稳定、高效地生产有十分重要的意义。因此，对于能够应用于实际过程对象的定位厂级控制回路振荡源的有效方法，其所具有的巨大的社会和经济效益将是不言而喻的。
[0004]目前，国内外有关厂级控制系统振荡源定位的方法主要有这样几种:SPCA方法、SICA方法、NMF方法、聚类方法(或称为:功率谱互相关分析法)、谱包络方法、基于转移熵的方法、基于互相关函数的方法、基于时域格兰杰因果性判断的方法等。我们通过仿真研宄发现，上面这些方法都有一个共同的问题，就是对不同的仿真对象，它们的正确率都不高。
【发明内容】
[0005]发明目的:提供一种流程工业的厂级控制回路振荡源定位方法，以解决现有技术检测准确率不高的问题，以实现对厂级控制回路振荡源进行有效的定位。
[0006]技术方案:一种流程工业的厂级控制回路振荡源定位方法，包括如下步骤:
[0007]步骤1、针对多回路控制系统，构建一包含多个回路信号振荡的原始回路变量数据集;
[0008]步骤2、采用基于频域非负矩阵分解或频域独立成分分析方法处理所述原始回路变量数据集，提取出包含相同频谱基的一个或多个待检回路集；每个待检回路变量数据集中包含若干控制回路的变量，每个变量对应一个控制回路，所述变量包括回路的偏差、被控量和控制器输出量；
[0009]步骤3、在待检回路集中针对每一个回路，进行回路的控制器输出量和被控量的二维联合概率分布曲面的计算，获得对应每一个待检回路集中所有回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值；
[0010]将每一个待检回路集中所计算出的各个回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，则每一个回路集的回路的控制器输出量和被控量的二维联合概率分布曲面上最大值的排序列表中，排在第一位的控制器输出量和被控量的二维联合概率分布曲面上的最大值所对应的回路，就是该回路集所包含的各个回路变量中所存在振荡的源。
[0011]有益效果:与现有的各种厂级控制回路振荡源定位方法相比，本发明进行厂级控制回路振荡源定位的准确性有了很大程度的提高。因此，利用本发明对厂级控制回路振荡源进行有效的定位，为消除厂级控制回路振荡所引起的严重危害提供了极大的便利。
【附图说明】
[0012]图1a至图1f是本发明实施例1的仿真结果示意图；图1a至图1f分别为控制回路的PV信号时域图、PV信号的频域图、总功率谱、两个振荡频率成分的强度因子分布图、0.0023Hz振荡频率成分的强度因子排序图、0.1Hz振荡频率成分的强度因子排序图；在实施例I中，控制回路2中的一阶惯性加延时对象中的开环放大系数K = 0.51、T = 50s。
[0013]图2a至图2f是本发明实施例2的仿真结果示意图；图2a至图2f分别为控制回路的PV信号时域图、PV信号的频域图、总功率谱、两个振荡频率成分的强度因子分布图、0.0023Hz振荡频率成分的强度因子排序图、0.1Hz振荡频率成分的强度因子排序图；在实施例2中，控制回路2中的一阶惯性加延时对象中的开环放大系数K = 1.01、T = 50s。
[0014]图3a至图3f是本发明实施例3的仿真结果示意图；图3a至图3f分别为控制回路的PV信号时域图、PV信号的频域图、总功率谱、两个振荡频率成分的强度因子分布图、
0.0023Hz振荡频率成分的强度因子排序图、0.1Hz振荡频率成分的强度因子排序图；在实施例3中，控制回路2中的一阶惯性加延时对象中的开环放大系数K = UT = 20s。
【具体实施方式】
[0015]结合附图描述本发明，本发明进行流程工业厂级控制回路振荡源定位的基本原理如下:
[0016]首先，假设一个具有N个回路的多回路控制系统中，发生了多回路的振荡。也就是说，存在M个发生了回路信号振荡情况的回路。显然，M SN。构建起一个仅仅包括了上述M个回路信号振荡的原始回路变量数据集A。数据集A中有M个由长度为S的采样数据序列所组成的元素。如果用矩阵Amxs来表示数据集A的话，则Amxs矩阵的维数为:MXS，即Amxs=[a1； a2,...，aM]T，其中aji = 1，2，...，Μ)表示第i个回路的采样数据序列。
[0017]其次，对数据集A，结合基于频域非负矩阵分解(NMF)或频域独立成分分析(SICA)等方法进行数据处理，提取出包含相同频谱基的一个或多个待检回路集:P(1)、
P(2)、......、P(i)、......、P(T)(这里，T:待检回路集的个数，其中T彡1，
T < Μ) ο
[0018]每一个待检回路集中包含有若干个控制回路的变量。进行频域非负矩阵分解(NMF)或频域独立成分分析(SICA)等方法进行数据处理的回路变量可以是回路偏差(E =SP-PV)或被控量(PV)，也可以是控制器的输出量(OP)。不管是那一种回路变量，每一个变量都对应一个控制回路。
[0019]设第i个待检回路集P(i)中回路的个数为Qi (这里，i的取值范围:1?T，Qi>l，Qi^ Μ)。
[0020]最后，在待检回路集中针对每一个回路，进行回路的控制器输出量(OP)和被控量(PV)的二维联合概率分布曲面的计算。
[0021]通过计算，获得对应每一个待检回路集中所有回路OP和PV的二维联合概率分布曲面上的最大值。将每一个待检回路集中所计算出的各个回路OP和PV的二维联合概率分布曲面上的最大值按从小到大的次序进行排序。则每一个回路集的回路OP和PV的二维联合概率分布曲面上最大值的排序列表中，排在第一位的OP和PV的二维联合概率分布曲面上的最大值所对应的回路，就是该回路集所包含的各个回路变量中所存在振荡的源。
[0022]经过上面三个步骤，就实现了厂级控制回路振荡源的定位。
[0023]下面用一个包含15个控制回路的仿真模型的仿真计算结果，对本发明中的厂级控制回路振荡源定位方法与基于频域非负矩阵分解(NMF)的厂级控制回路振荡源定位方法的定位效果进行对比。
[0024]仿真参数设置情况为:采用NMF方法，在仿真模型的回路2中 设置了控制阀门粘滞故障(S = 10，J = 3)，在仿真模型的回路7中设置了正弦干扰信号(幅度A= 10，频率f= 0.1)。对比控制回路2的一阶惯性加延时对象(first-order plus dead-time model，F(PDT)中的一阶惯性环节的增益和时间常数分别为:K = 0.51、T = 50s，K = 1.01、T =50s, K = 1、T = 20s等3种情况下的回路振荡指数(强度因子SF)，结果如图1、图2及图3所示。
[0025]从上面的图1、图2及图3可以看到，对应于控制回路2的粘滞所引起的振荡频率分别为0.0023Hz、0.0056Hz及0.0041Hz。但NMF方法所提出的控制回路振荡源定位规则，对应于K = 0.51、T = 50s，K = 1.01、T = 50s及K = 1、T = 20s三种情况的振荡源分别是:回路6、回路2及回路3。显然，在第I及第3这两种情况下的定位结果是错误的，而只有第2种情况下的定位结果是正确的。
[0026]如果用本发明的方法，计算过程及结果为:
[0027](I)对应于K = 0.51、T = 50s，振荡频率0.0023Hz所对应的待检回路集中包括7个回路:6、3、1、2、4、5、11(注意:对于回路15等，由于其强度因子SF相对于前面选中的7个回路的强度因子SF来讲很小，故没有被选进待检回路集中)。在待检回路集中针对每一个回路，进行回路的控制器输出量(OP)和被控量(PV)的二维联合概率分布曲面的计算，获得对应待检回路集中所有回路OP和PV的二维联合概率分布曲面上的最大值分别为:
0.24742(回路 I)、0.19408(回路 2)、0.35684(回路 3)、2.9778(回路 4)、4.5597(回路 5)、
1.164 (回路 6)、4.6308(回路 11)。
[0028]将该待检回路集中所计算出的各个回路OP和PV的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，得到:0.19408(回路2)、0.24742(回路I)、0.35684 ( 0路 3)、1.164 (回路 6)、2.9778(回路 4)、4.5597(回路 5)、4.6308(回路 11)。
[0029]按本发明的方法，则待检回路集中所有回路各自的OP和PV的二维联合概率分布曲面上的最大值排序列表(按从小到大的次序进行排序)中，排在第一位的OP和PV的二维联合概率分布曲面上的最大值(即:0.19408)所对应的是回路2，这就是该待检回路集所包含的各个回路变量当中所存在振荡的源头。这与实际情况一致，结果正确。
[0030](2)对应于K=L OUT = 50s，振荡频率0.0056Hz所对应的待检回路集中包括6个回路:2、4、1、3、6、5(注意:对于回路11等，由于其强度因子SF相对于前面选中的6个回路的强度因子SF来讲很小，故没有被选进待检回路集中)。在待检回路集中针对每一个回路，进行回路的控制器输出量(OP)和被控量(PV)的二维联合概率分布曲面的计算，获得对应待检回路集中所有回路最优广义被控对象模型辨识结果所对应的最小拟合误差分别为:0.20239(回路 I)、0.050485(回路 2)、0.58279(回路 3)、0.93445(回路 4)、5.6716(回路5)、2.0137(回路 6) ο
[0031]将该待检回路集中所计算出的各个回路OP和PV的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，得到:0.050485(回路2)、0.20239(回路I)、0.58279 ( 0路 3)、0.93445(回路 4)、2.0137(回路 6)、5.6716(回路 5)。
[0032]按本发明的方法，则待检回路集中所有回路各自的OP和PV的二维联合概率分布曲面上的最大值排序列表(按从小到大的次序进行排序)中，排在第一位的OP和PV的二维联合概率分布曲面上的最大值(即:0.050485)所对应的是回路2，这就是该待检回路集所包含的各个回路变量当中所存在振荡的源头。这与实际情况一致，结果正确。
[0033](3)对应于K=1、T = 20s，振荡频率0.0041Hz所对应的待检回路集中包括6个回路:4、2、1、5、3、6(注意:对于回路11等，由于其强度因子SF相对于前面选中的6个回路的强度因子SF来讲很小，故没有被选进待检回路集中)。在待检回路集中针对每一个回路，进行回路的控制器输出量(OP)和被控量(PV)的二维联合概率分布曲面的计算，获得对应待检回路集中所有回路最优广义被控对象模型辨识结果所对应的最小拟合误差分别为:0.49315 (回路 I)、0.088372 (回路 2)、0.65521 (回路 3)、1.7609 (回路 4)、6.3341 (回路 5)、
2.605 (回路 6)。
[0034]将该待检回路集中所计算出的各个回路OP和PV的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，得到:0.088372(回路2)、0.49315(回路I)、0.65521 (回路 3)、1.7609(回路 4)、2.605 (回路 6)、6.3341 (回路 5)。
[0035]按本发明的方法，则待检回路集中所有回路各自的OP和PV的二维联合概率分布曲面上的最大值排序列表(按从小到大的次序进行排序)中，排在第一位的OP和PV的二维联合概率分布曲面上的最大值(即:0.088372)所对应的是回路2，这就是该待检回路集所包含的各个回路变量当中所存在振荡的源头。这与实际情况一致，结果正确。
[0036]以上详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种等同变换，这些等同变换均属于本发明的保护范围。另外需要说明的是，在上述【具体实施方式】中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。
【主权项】
1.一种流程工业的厂级控制回路振荡源定位方法，其特征在于，包括如下步骤:步骤1、针对多回路控制系统，构建一包含多个回路信号振荡的原始回路变量数据集；步骤2、采用基于频域非负矩阵分解或频域独立成分分析方法处理所述原始回路变量数据集，提取出包含相同频谱基的一个或多个待检回路集；每个待检回路变量数据集中包含若干控制回路的变量，每个变量对应一个控制回路，所述变量包括回路的偏差、被控量和控制器输出量； 步骤3、在待检回路集中针对每一个回路，进行回路的控制器输出量和被控量的二维联合概率分布曲面的计算，获得对应每一个待检回路集中所有回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值； 将每一个待检回路集中所计算出的各个回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，则每一个回路集的回路的控制器输出量和被控量的二维联合概率分布曲面上最大值的排序列表中，排在第一位的控制器输出量和被控量的二维联合概率分布曲面上的最大值所对应的回路，就是该回路集所包含的各个回路变量中所存在振荡的源。
【专利摘要】本发明公开了一种流程工业的厂级控制回路振荡源定位方法，包括如下步骤：针对多回路控制系统，构建一包含多个回路信号振荡的原始回路变量数据集；采用基于频域非负矩阵分解或频域独立成分分析方法处理所述原始回路变量数据集，提取出包含相同频谱基的一个或多个待检回路集；在待检回路集中针对每一个回路，进行回路的控制器输出量和被控量的二维联合概率分布曲面的计算，获得对应每一个待检回路集中所有回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值。与现有技术相比，本发明具有准确程度高的优点，可对厂级控制回路振荡源进行有效的定位，为消除厂级控制回路振荡所引起的危害提供了极大的便利。
【IPC分类】G05B23/02
【公开号】CN104898649
【申请号】CN201510234239
【发明人】曹鸣
【申请人】东南大学
【公开日】2015年9月9日
【申请日】2015年5月8日
【技术领域】
[0001]本发明属于控制或调节系统领域，尤其是一种基于模型识别方法的流程工业厂级控制回路振荡源定位方法。
【背景技术】
[0002]由于现代过程工业日益复杂及控制器日常维护不足等方面的原因，电力、石油、化工、钢铁、造纸、水泥等工业过程普遍存在控制回路性能差的现象，国外有研宄表明，工业过程中有1/3的控制回路处于振荡状态中。回路振荡是一种明显的控制性能退化的表现形式，有多种可能的原因会引起控制回路振荡:控制阀的非线性故障，如摩擦、死区、反冲及粘滞等现象；控制参数调整不适当，使得系统趋于不稳而引起振荡；来自外部振荡干扰和回路自身的反馈机制，这种情况通常在控制回路间有较强耦合关联的系统中见到，一个振荡源会通过耦合路径传递到其它回路中去；高度的非线性过程；在用的MPC控制器与对象模型不匹配。控制回路中振荡的存在将导致产品质量下降、能耗增加、设备磨损加快等一系列问题。尤其在现代企业规模日益扩大的情况下，控制回路数目成千上万，一旦由于某种异常因素(例如阀门出现了粘滞或外部扰动等)导致了在这样大的控制系统中的某一位置上发生了振荡，由于过程系统中的各个处理装置间存在着复杂的耦合，则这种出现在局部位置上的振荡将被快速传递扩散开来，产生所谓的“厂级控制回路振荡”现象。一旦“厂级控制回路振荡”出现，不但振荡的各种危害将被成倍的放大，而且还造成了原始的那个振荡源所在的控制回路被淹没在许许多多具有非常相似振荡特征的控制回路当中。也就是说，一旦“厂级控制回路振荡”出现，
[0003]如果没有一种行之有效的定位厂级控制回路振荡源的方法，将导致“厂级控制回路振荡”现象将在一个相当长的时间内存在。显然，要消除“厂级控制回路振荡”所引起的这种严重危害，首先要找到其来源，然后对其定性，最后根据不同的原因采取相应的措施，加以解决。值得一提的是，发电机组等大型过程装置的运行过程常常伴随十分苛刻的生产条件或环境，如高温、高压，因而运行的安全性十分重要。由此可见，定位厂级控制回路振荡源方法的研宄对于电力、石油、化工、钢铁、造纸、水泥等工业过程的安全、稳定、高效地生产有十分重要的意义。因此，对于能够应用于实际过程对象的定位厂级控制回路振荡源的有效方法，其所具有的巨大的社会和经济效益将是不言而喻的。
[0004]目前，国内外有关厂级控制系统振荡源定位的方法主要有这样几种:SPCA方法、SICA方法、NMF方法、聚类方法(或称为:功率谱互相关分析法)、谱包络方法、基于转移熵的方法、基于互相关函数的方法、基于时域格兰杰因果性判断的方法等。我们通过仿真研宄发现，上面这些方法都有一个共同的问题，就是对不同的仿真对象，它们的正确率都不高。
【发明内容】
[0005]发明目的:提供一种流程工业的厂级控制回路振荡源定位方法，以解决现有技术检测准确率不高的问题，以实现对厂级控制回路振荡源进行有效的定位。
[0006]技术方案:一种流程工业的厂级控制回路振荡源定位方法，包括如下步骤:
[0007]步骤1、针对多回路控制系统，构建一包含多个回路信号振荡的原始回路变量数据集;
[0008]步骤2、采用基于频域非负矩阵分解或频域独立成分分析方法处理所述原始回路变量数据集，提取出包含相同频谱基的一个或多个待检回路集；每个待检回路变量数据集中包含若干控制回路的变量，每个变量对应一个控制回路，所述变量包括回路的偏差、被控量和控制器输出量；
[0009]步骤3、在待检回路集中针对每一个回路，进行回路的控制器输出量和被控量的二维联合概率分布曲面的计算，获得对应每一个待检回路集中所有回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值；
[0010]将每一个待检回路集中所计算出的各个回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，则每一个回路集的回路的控制器输出量和被控量的二维联合概率分布曲面上最大值的排序列表中，排在第一位的控制器输出量和被控量的二维联合概率分布曲面上的最大值所对应的回路，就是该回路集所包含的各个回路变量中所存在振荡的源。
[0011]有益效果:与现有的各种厂级控制回路振荡源定位方法相比，本发明进行厂级控制回路振荡源定位的准确性有了很大程度的提高。因此，利用本发明对厂级控制回路振荡源进行有效的定位，为消除厂级控制回路振荡所引起的严重危害提供了极大的便利。
【附图说明】
[0012]图1a至图1f是本发明实施例1的仿真结果示意图；图1a至图1f分别为控制回路的PV信号时域图、PV信号的频域图、总功率谱、两个振荡频率成分的强度因子分布图、0.0023Hz振荡频率成分的强度因子排序图、0.1Hz振荡频率成分的强度因子排序图；在实施例I中，控制回路2中的一阶惯性加延时对象中的开环放大系数K = 0.51、T = 50s。
[0013]图2a至图2f是本发明实施例2的仿真结果示意图；图2a至图2f分别为控制回路的PV信号时域图、PV信号的频域图、总功率谱、两个振荡频率成分的强度因子分布图、0.0023Hz振荡频率成分的强度因子排序图、0.1Hz振荡频率成分的强度因子排序图；在实施例2中，控制回路2中的一阶惯性加延时对象中的开环放大系数K = 1.01、T = 50s。
[0014]图3a至图3f是本发明实施例3的仿真结果示意图；图3a至图3f分别为控制回路的PV信号时域图、PV信号的频域图、总功率谱、两个振荡频率成分的强度因子分布图、
0.0023Hz振荡频率成分的强度因子排序图、0.1Hz振荡频率成分的强度因子排序图；在实施例3中，控制回路2中的一阶惯性加延时对象中的开环放大系数K = UT = 20s。
【具体实施方式】
[0015]结合附图描述本发明，本发明进行流程工业厂级控制回路振荡源定位的基本原理如下:
[0016]首先，假设一个具有N个回路的多回路控制系统中，发生了多回路的振荡。也就是说，存在M个发生了回路信号振荡情况的回路。显然，M SN。构建起一个仅仅包括了上述M个回路信号振荡的原始回路变量数据集A。数据集A中有M个由长度为S的采样数据序列所组成的元素。如果用矩阵Amxs来表示数据集A的话，则Amxs矩阵的维数为:MXS，即Amxs=[a1； a2,...，aM]T，其中aji = 1，2，...，Μ)表示第i个回路的采样数据序列。
[0017]其次，对数据集A，结合基于频域非负矩阵分解(NMF)或频域独立成分分析(SICA)等方法进行数据处理，提取出包含相同频谱基的一个或多个待检回路集:P(1)、
P(2)、......、P(i)、......、P(T)(这里，T:待检回路集的个数，其中T彡1，
T < Μ) ο
[0018]每一个待检回路集中包含有若干个控制回路的变量。进行频域非负矩阵分解(NMF)或频域独立成分分析(SICA)等方法进行数据处理的回路变量可以是回路偏差(E =SP-PV)或被控量(PV)，也可以是控制器的输出量(OP)。不管是那一种回路变量，每一个变量都对应一个控制回路。
[0019]设第i个待检回路集P(i)中回路的个数为Qi (这里，i的取值范围:1?T，Qi>l，Qi^ Μ)。
[0020]最后，在待检回路集中针对每一个回路，进行回路的控制器输出量(OP)和被控量(PV)的二维联合概率分布曲面的计算。
[0021]通过计算，获得对应每一个待检回路集中所有回路OP和PV的二维联合概率分布曲面上的最大值。将每一个待检回路集中所计算出的各个回路OP和PV的二维联合概率分布曲面上的最大值按从小到大的次序进行排序。则每一个回路集的回路OP和PV的二维联合概率分布曲面上最大值的排序列表中，排在第一位的OP和PV的二维联合概率分布曲面上的最大值所对应的回路，就是该回路集所包含的各个回路变量中所存在振荡的源。
[0022]经过上面三个步骤，就实现了厂级控制回路振荡源的定位。
[0023]下面用一个包含15个控制回路的仿真模型的仿真计算结果，对本发明中的厂级控制回路振荡源定位方法与基于频域非负矩阵分解(NMF)的厂级控制回路振荡源定位方法的定位效果进行对比。
[0024]仿真参数设置情况为:采用NMF方法，在仿真模型的回路2中 设置了控制阀门粘滞故障(S = 10，J = 3)，在仿真模型的回路7中设置了正弦干扰信号(幅度A= 10，频率f= 0.1)。对比控制回路2的一阶惯性加延时对象(first-order plus dead-time model，F(PDT)中的一阶惯性环节的增益和时间常数分别为:K = 0.51、T = 50s，K = 1.01、T =50s, K = 1、T = 20s等3种情况下的回路振荡指数(强度因子SF)，结果如图1、图2及图3所示。
[0025]从上面的图1、图2及图3可以看到，对应于控制回路2的粘滞所引起的振荡频率分别为0.0023Hz、0.0056Hz及0.0041Hz。但NMF方法所提出的控制回路振荡源定位规则，对应于K = 0.51、T = 50s，K = 1.01、T = 50s及K = 1、T = 20s三种情况的振荡源分别是:回路6、回路2及回路3。显然，在第I及第3这两种情况下的定位结果是错误的，而只有第2种情况下的定位结果是正确的。
[0026]如果用本发明的方法，计算过程及结果为:
[0027](I)对应于K = 0.51、T = 50s，振荡频率0.0023Hz所对应的待检回路集中包括7个回路:6、3、1、2、4、5、11(注意:对于回路15等，由于其强度因子SF相对于前面选中的7个回路的强度因子SF来讲很小，故没有被选进待检回路集中)。在待检回路集中针对每一个回路，进行回路的控制器输出量(OP)和被控量(PV)的二维联合概率分布曲面的计算，获得对应待检回路集中所有回路OP和PV的二维联合概率分布曲面上的最大值分别为:
0.24742(回路 I)、0.19408(回路 2)、0.35684(回路 3)、2.9778(回路 4)、4.5597(回路 5)、
1.164 (回路 6)、4.6308(回路 11)。
[0028]将该待检回路集中所计算出的各个回路OP和PV的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，得到:0.19408(回路2)、0.24742(回路I)、0.35684 ( 0路 3)、1.164 (回路 6)、2.9778(回路 4)、4.5597(回路 5)、4.6308(回路 11)。
[0029]按本发明的方法，则待检回路集中所有回路各自的OP和PV的二维联合概率分布曲面上的最大值排序列表(按从小到大的次序进行排序)中，排在第一位的OP和PV的二维联合概率分布曲面上的最大值(即:0.19408)所对应的是回路2，这就是该待检回路集所包含的各个回路变量当中所存在振荡的源头。这与实际情况一致，结果正确。
[0030](2)对应于K=L OUT = 50s，振荡频率0.0056Hz所对应的待检回路集中包括6个回路:2、4、1、3、6、5(注意:对于回路11等，由于其强度因子SF相对于前面选中的6个回路的强度因子SF来讲很小，故没有被选进待检回路集中)。在待检回路集中针对每一个回路，进行回路的控制器输出量(OP)和被控量(PV)的二维联合概率分布曲面的计算，获得对应待检回路集中所有回路最优广义被控对象模型辨识结果所对应的最小拟合误差分别为:0.20239(回路 I)、0.050485(回路 2)、0.58279(回路 3)、0.93445(回路 4)、5.6716(回路5)、2.0137(回路 6) ο
[0031]将该待检回路集中所计算出的各个回路OP和PV的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，得到:0.050485(回路2)、0.20239(回路I)、0.58279 ( 0路 3)、0.93445(回路 4)、2.0137(回路 6)、5.6716(回路 5)。
[0032]按本发明的方法，则待检回路集中所有回路各自的OP和PV的二维联合概率分布曲面上的最大值排序列表(按从小到大的次序进行排序)中，排在第一位的OP和PV的二维联合概率分布曲面上的最大值(即:0.050485)所对应的是回路2，这就是该待检回路集所包含的各个回路变量当中所存在振荡的源头。这与实际情况一致，结果正确。
[0033](3)对应于K=1、T = 20s，振荡频率0.0041Hz所对应的待检回路集中包括6个回路:4、2、1、5、3、6(注意:对于回路11等，由于其强度因子SF相对于前面选中的6个回路的强度因子SF来讲很小，故没有被选进待检回路集中)。在待检回路集中针对每一个回路，进行回路的控制器输出量(OP)和被控量(PV)的二维联合概率分布曲面的计算，获得对应待检回路集中所有回路最优广义被控对象模型辨识结果所对应的最小拟合误差分别为:0.49315 (回路 I)、0.088372 (回路 2)、0.65521 (回路 3)、1.7609 (回路 4)、6.3341 (回路 5)、
2.605 (回路 6)。
[0034]将该待检回路集中所计算出的各个回路OP和PV的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，得到:0.088372(回路2)、0.49315(回路I)、0.65521 (回路 3)、1.7609(回路 4)、2.605 (回路 6)、6.3341 (回路 5)。
[0035]按本发明的方法，则待检回路集中所有回路各自的OP和PV的二维联合概率分布曲面上的最大值排序列表(按从小到大的次序进行排序)中，排在第一位的OP和PV的二维联合概率分布曲面上的最大值(即:0.088372)所对应的是回路2，这就是该待检回路集所包含的各个回路变量当中所存在振荡的源头。这与实际情况一致，结果正确。
[0036]以上详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种等同变换，这些等同变换均属于本发明的保护范围。另外需要说明的是，在上述【具体实施方式】中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。
【主权项】
1.一种流程工业的厂级控制回路振荡源定位方法，其特征在于，包括如下步骤:步骤1、针对多回路控制系统，构建一包含多个回路信号振荡的原始回路变量数据集；步骤2、采用基于频域非负矩阵分解或频域独立成分分析方法处理所述原始回路变量数据集，提取出包含相同频谱基的一个或多个待检回路集；每个待检回路变量数据集中包含若干控制回路的变量，每个变量对应一个控制回路，所述变量包括回路的偏差、被控量和控制器输出量； 步骤3、在待检回路集中针对每一个回路，进行回路的控制器输出量和被控量的二维联合概率分布曲面的计算，获得对应每一个待检回路集中所有回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值； 将每一个待检回路集中所计算出的各个回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值按从小到大的次序进行排序，则每一个回路集的回路的控制器输出量和被控量的二维联合概率分布曲面上最大值的排序列表中，排在第一位的控制器输出量和被控量的二维联合概率分布曲面上的最大值所对应的回路，就是该回路集所包含的各个回路变量中所存在振荡的源。
【专利摘要】本发明公开了一种流程工业的厂级控制回路振荡源定位方法，包括如下步骤：针对多回路控制系统，构建一包含多个回路信号振荡的原始回路变量数据集；采用基于频域非负矩阵分解或频域独立成分分析方法处理所述原始回路变量数据集，提取出包含相同频谱基的一个或多个待检回路集；在待检回路集中针对每一个回路，进行回路的控制器输出量和被控量的二维联合概率分布曲面的计算，获得对应每一个待检回路集中所有回路的控制器输出量和被控量的二维联合概率分布曲面上的最大值。与现有技术相比，本发明具有准确程度高的优点，可对厂级控制回路振荡源进行有效的定位，为消除厂级控制回路振荡所引起的危害提供了极大的便利。
【IPC分类】G05B23/02
【公开号】CN104898649
【申请号】CN201510234239
【发明人】曹鸣
【申请人】东南大学
【公开日】2015年9月9日
【申请日】2015年5月8日
最新回复(0)