面向高端燃煤发电装备的变量协同动态报警阈值优化方法技术

本发明专利技术公开了一种面向高端燃煤发电装备的变量协同动态报警阈值优化方法。本发明专利技术针对高端燃煤发电过程中的磨煤机设备，考虑了非平稳变工况运行下的变量协同动态报警阈值优化。借助时序分割有序聚类方法将非平稳过程的历史正常数据划分至若干工况(包括稳态、过渡过程)，每种工况用一个多元高斯分布模型进行数学表征。在线时，根据时序分割有序聚类方法将当前数据分类到所属工况，选取当前工况的多元高斯模型，从而根据条件概率进行阈值设计。该方法可满足全工况下的准确报警，对保证过程安全运行和人员安全有重大意义。全运行和人员安全有重大意义。全运行和人员安全有重大意义。

【技术实现步骤摘要】
面向高端燃煤发电装备的变量协同动态报警阈值优化方法


[0001]本专利技术属于非平稳高端燃煤发电过程报警管理领域，特别是涉及了一种变量协同动态报警阈值优化方法。

技术介绍

[0002]报警系统对现代工业过程的安全、高效运行起到至关重要的作用。例如分布式控制系统(DCS)和监测控制和数据采集系统(SCADA)，对一些重要的过程变量与其配置的高(低)报警阈值进行实时比较，如果变量超出阈值，则触发报警以告知操作人员。
[0003]一般来说，报警阈值是根据出厂保护定值或者在调试期间仔细确定的，如若设置不当，会导致误报和漏报的情况发生，关乎到整个报警系统的性能。然而工业过程投入生产运行后，缺乏对报警阈值的定期、大规模的维护与更改。目前，现有的报警系统普遍存在的问题是报警泛滥，产生的报警分为干扰报警和正确报警两类。首先，大量的报警属于干扰报警一类，它们不提供任何有用的信息，只会分散操作员的注意力。其次，即便所有报警都是正确的，操作员也无法在短时间内处理大量的报警信息，只能忽略许多已发生的报警。在真正发生重大灾害时，故障根源很可能埋没在大量的报警之中，无法及时被定位和应对。
[0004]针对上述报警系统的报警泛滥问题，学术界提出了很多基于数理统计的阈值优化设计方法，例如最为经典的是三倍标准差法准则。另外一种主流的阈值设计思路是对误报率和漏报率进行权衡，此时的阈值设计实际上是一个优化问题。然而，上述阈值设计方法都是从单变量角度出发，忽略了工业过程固有的耦合性。该特性在进行报警阈值设计时不容忽视，否则会容易引发误报率和漏报率。为此，有学者从变量协同的角度出发，分别计算两个或多个过程变量之间与报警数据之间的相关性，通过设置阈值使差异尽可能的小，说明报警情况比较反应过程特性，但无法根据不同工况来动态优化阈值。

技术实现思路

[0005]本专利技术的目的在于针对现有报警阈值设计的不合理，提供一种面向高端燃煤发电装备的变量协同动态报警阈值优化方法。该方法考虑到固定的阈值无法适应不同工况，借助时序分割有序聚类方法将非平稳过程的历史正常数据划分至若干工况(包括稳态、过渡过程)，每种工况用一个多元高斯模型表征，该工况下的每个过程变量借助条件概率分布来估计出当前其他变量取值情况下的动态阈值。在线时，根据时序分割有序聚类方法将当前数据分类到所属工况，选取当前工况的多元高斯模型，从而根据条件概率进行阈值设计。本专利技术为非平稳过程下的阈值优化工作提出了一种新的分析思路，该方法可满足全工况下的准确报警，对保证过程安全运行和人员安全有重大意义。
[0006]本专利技术的目的是通过以下技术方案来实现：
[0007]面向高端燃煤发电装备的变量协同动态报警阈值优化方法，该方法包括以下步骤：
[0008](1)获得非平稳过程数据：获取高端燃煤发电装备正常运行时J个过程变量的T(t
＝1，2，...，T)个连续采样数据，其二维数据矩阵表示为X(T
×
J)；
[0009](2)利用时序分割有序聚类方法将步骤(1)获得的非平稳过程数据聚类划分工况类别，并获得每个工况类别的多元高斯分布的数学表征其中，μ
i
表示为第i个工况类别的均值，∑
i
表示为第i个工况类别的协方差。
[0010](3)实时采集高端燃煤发电装备运行时新的过程变量数据x
new
(1
×
J)，对x
new
(1
×
J)进行时序拓展，得到{x
new
‑
ω+1
，
…
，x
new
}，ω为时序拓展阶次；并通过时序分割有序聚类方法的E步将当前新的过程变量数据分配到步骤(2)聚类的某一工况类别i中；
[0011](4)基于对应工况类别的多元高斯分布，对时序拓展后的过程变量数据进行变量协同阈值设计，具体为：将Jω维的过程变量组表示为当前要阈值设计的变量x1与除x1之外的所有过程变量的x2组合其中x1维度为1；x2维度为Jω
‑
1；
[0012]均值向量对应拆分为：
[0013][0014]协方差矩阵对应拆分为：
[0015][0016]其中，∑
11
表示变量x1与变量x1的协方差矩阵，∑
12
表示变量x1与其他变量x2的协方差矩阵，∑
21
表示其他变量x2与变量x1的协方差矩阵，∑
11
表示其他变量x2与其他变量x2的协方差矩阵。
[0017]变量x1在其他变量x2取值给定时，呈现条件高斯分布：
[0018][0019]其中，条件均值为：
[0020][0021]条件方差为：
[0022][0023]根据条件高斯分布确定变量x1的报警阈值，对每个变量逐一进行变量协同阈值设计，获得所有变量的报警阈值。
[0024]进一步地，步骤1中，测量变量包括：给煤量、磨煤机电流、旋转分离器电流、出口压力、进口一次风温、进口一次风压、密封风压力、出口风粉混合物温度。
[0025]进一步地，所述步骤2由以下子步骤来实现：
[0026](2.1)参数选择：结合贝叶斯信息准则(BIC)指标与机理事先确定类别个数I、时序拓展阶次ω、一范数正则化参数λ、时序一致约束系数β；
[0027](2.2)进行聚类：根据选定的时序拓展阶次ω，非平稳过程数据中每个样本x
t
被时序扩展成{x
t
‑
ω+1
，
…
，x
t
}的形式；非平稳过程数据的二维数据矩阵表示为X(T
×
Jω)；根据时序分割有序聚类方法来进行聚类，获得I(i＝1，2，...，I)个工况类别与其对应的多元高斯分布表征
[0028]进一步地，所述步骤4中，根据条件高斯分布确定变量x1的报警阈值，具体为：
[0029]取置信度为95％，则x1的高报警阈值为μ
1|2
+2∑
1|2
，低报警阈值为μ
1|2
‑
2∑
1|2
。
[0030]与现有技术相比，本专利技术的有益效果是：本专利技术为针对高端燃煤发电非平稳运行过程数据的报警变量阈值优化技术提供了新的研究思路。通过提出一种基于工况指示变量将原本的非平稳运行数据重构成若干个平稳的工况段，每个段内各自进行概率统计，获得该工况下的阈值。在线应用时，根据工况指示变量的值来匹配相应的阈值，实现根据当前工况进行动态的阈值设置，有助于降低误报。所提方法针对实际工业数据进行了充分的实验研究，并获得了成功应用。
附图说明
[0031]图1是本专利技术所述变量协同动态报警阈值优化的流程图。
[0032]图2是本专利技术所述工况模态划分结果与相应阈值示意图。
[0033]图3是本专利技术所述有阈值优化的故障数据应用结果图。
[0034]图本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种面向高端燃煤发电装备的变量协同动态报警阈值优化方法，其特征在于，该方法包括以下步骤：(1)获得非平稳过程数据：获取高端燃煤发电装备正常运行时J个过程变量的T(t＝1，2，...，T)个连续采样数据，其二维数据矩阵表示为X(T
×
J)；(2)利用时序分割有序聚类方法将步骤(1)获得的非平稳过程数据聚类划分工况类别，并获得每个工况类别的多元高斯分布的数学表征其中，μ
i
表示为第i个工况类别的均值，∑
i
表示为第i个工况类别的协方差。(3)实时采集高端燃煤发电装备运行时新的过程变量数据x
new
(1
×
J)，对x
new
(1
×
J)进行时序拓展，得到{x
new
‑
ω+1
，
…
，x
new
}，ω为时序拓展阶次；并通过时序分割有序聚类方法的E步将当前新的过程变量数据分配到步骤(2)聚类的某一工况类别i中；(4)基于对应工况类别的多元高斯分布，对时序拓展后的过程变量数据进行变量协同阈值设计，具体为：将Jω维的过程变量组表示为当前要阈值设计的变量x1与除x1之外的所有过程变量的x2组合其中x1维度为1；x2维度为Jω
‑
1；均值向量对应拆分为：协方差矩阵对应拆分为：其中，∑
11
表示变量x1与变量x1的协方差矩阵，∑
12
表示变量x1与其他变量x2的协方差矩阵，∑
21
表示其他变量x2与变量x1的协方差矩阵，∑
11

【专利技术属性】
技术研发人员：赵春晖，赵诣，李清毅，王文海，
申请(专利权)人：浙江大学，
类型：发明
国别省市：