CN110648480B - 一种基于变化速率的单变量报警系统及方法 - Google Patents

一种基于变化速率的单变量报警系统及方法 Download PDF

Info

Publication number
CN110648480B
CN110648480B CN201910969211.9A CN201910969211A CN110648480B CN 110648480 B CN110648480 B CN 110648480B CN 201910969211 A CN201910969211 A CN 201910969211A CN 110648480 B CN110648480 B CN 110648480B
Authority
CN
China
Prior art keywords
rate
data
alarm
window
change
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
Application number
CN201910969211.9A
Other languages
English (en)
Other versions
CN110648480A (zh
Inventor
王建东
曲俊潇
王振
白星振
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN201910969211.9A priority Critical patent/CN110648480B/zh
Publication of CN110648480A publication Critical patent/CN110648480A/zh
Application granted granted Critical
Publication of CN110648480B publication Critical patent/CN110648480B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B7/00Signalling systems according to more than one of groups G08B3/00 - G08B6/00; Personal calling systems according to more than one of groups G08B3/00 - G08B6/00
    • G08B7/06Signalling systems according to more than one of groups G08B3/00 - G08B6/00; Personal calling systems according to more than one of groups G08B3/00 - G08B6/00 using electric transmission, e.g. involving audible and visible signalling through the use of sound and light sources
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/18Status alarms

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Alarm Systems (AREA)

Abstract

本公开提供了一种基于变化速率的单变量报警系统及方法,根据分布稳定的历史变化速率样本集合,确定速率信号的报警阈值;实时采集被监测的单变量数据,利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与所述离线单元的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警指令,并发出报警信号。克服了传统的单变量报警系统大多通过监测模拟信号的幅值变化来判断系统是否发生异常,当信号幅值超出报警阈值则触发报警的问题。

Description

一种基于变化速率的单变量报警系统及方法
技术领域
本公开属于工业大数据分析技术领域,尤其涉及一种基于变化速率的单变量报警系统及方法。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息,并不必然构成在先技术。
现代工业及电力能源行业在我国的国民生产中占有十分重要的地位,而安全问题则是其中最重要的问题。现代工业设备通常由计算机化的报警系统监控,报警系统对现代工业设施的安全和有效运行至关重要,报警系统的性能越来越受到工业界和学术界的关注。报警系统是一种监视设备运行工况并提醒工作人员及时处理异常状况的系统。
发明人在研究中发现,传统的越限报警方式存在报警延迟时间较长、受噪声干扰较大、误报警和漏报警较多等缺点,会导致报警泛滥、报警系统性能低下等问题,严重时可能会造成财产损失、人员伤亡等重大事故的发生。目前,报警系统常用滤波、延时和死区三种方法来减少漏报率、误报率,以提高报警系统性能。但这些方法仍存在一定的局限性:第一,没有对模拟信号的变化速率进行监测,容易导致误报警、漏报警、报警延迟时间较长等问题,严重影响报警系统的性能;第二,无法确定选取多少工业历史数据才能保证统计分布不再发生改变,未能准确判定数据统计分布的稳定性;第三,对于变化速率较为平稳的模拟信号,传统的越限报警方式存在较长的报警延迟时间,并且在噪声的干扰下,易造成误报警和漏报警,灵敏性较差。
发明内容
本公开为了解决上述问题,提出了一种基于变化速率的单变量报警系统及方法,本公开通过判断分布稳定的历史变化速率样本集合来确定速率信号的报警阈值,利用基于时间序列的分段线性表示法计算过程信号的变化速率并判断系统是否发生异常,提高了报警系统的性能。
根据一些实施例,本公开采用如下技术方案:
一种基于变化速率的单变量报警系统,包括在线单元、离线单元和报警单元,其中:
所述离线单元,被配置为根据分布稳定的历史变化速率样本集合,确定速率信号的报警阈值;
所述在线单元,被配置为实时采集被监测的单变量数据,利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与所述离线单元的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警指令;
所述报警单元,被配置为接收所述报警指令,并发出报警信号。
本公开利用上述技术方案,一方面能够充分利用有限的历史数据,判断过程变量的历史速率集合的分布稳定性,并通过稳定的速率集合计算速率报警阈值,合理、灵活、动态设置报警阈值;另一方面利用基于时间序列的分段线性表示法计算过程信号的变化速率并判断系统是否发生异常,提高了报警系统的性能;克服了传统的单变量报警系统大多通过监测模拟信号的幅值变化来判断系统是否发生异常,当信号幅值超出报警阈值则触发报警,容易造成误判的问题。
作为可选择的实施方式,所述离线单元包括:
初步选择模块,被配置为选取一定长度历史数据作为训练样本数据;
分段模块,被配置为选定滑动窗口宽度,对每个窗口内的数据采用自底向上法进行分段,并对各段数据进行最小二乘拟合;
计算模块,被配置为采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,并计算每段数据的变化速率,组成速率集合;
处理模块,被配置为判断速率集合是否稳定,若稳定,通过分位数确定速率报警阈值,若不稳定,逐步增加训练样本的选取长度,直到稳定。
作为可选择的实施方式,所述在线单元,包括:
采集模块,被配置为实时采集被监测的单变量数据;
计算模块,被配置为利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间;
比较判断模块,被配置为将置信区间上下限与所述离线单元的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警指令。
作为可选择的实施方式,所述报警单元为声光报警器。
一种基于变化速率的单变量报警方法,包括以下步骤:
实时采集被监测的单变量数据,利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与预先设定的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警信号。
作为可选择的实施方式,预先设定的速率报警阈值通过以下过程得到:
(1)选取一定长度的历史数据作为训练样本数据;
(2)选定滑动窗口的宽度,对每个窗口内的数据采用自底向上法进行分段,并对各段数据进行最小二乘拟合;
(3)采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,并计算每段数据的变化速率,组成速率集合;
(4)根据得到的速率集合判断是否稳定,若稳定,通过分位数确定速率报警阈值,若不稳定,则增加样本长度,重复步骤(2)-(4)。
作为可选择的实施方式,所述步骤(3)中,每个窗口内的最优分段数采用均方根误差算法进行确定,表示各窗口内分段数与加权均方根误差所构成的损失函数,当所述损失函数最小时,得到对应窗口内数据的最优分段数。
作为可选择的实施方式,所述步骤(3)中,计算每段数据的变化速率的具体过程是,根据最优分段数结果,对每段数据进行一阶最小二乘拟合,每段数据拟合直线的斜率即为的该段的变化速率。
作为可选择的实施方式,所述步骤(4)中,判断速率集合是否稳定的过程包括:将速率集合的最大值和最小值等长度间隔分为多个小区间,通过贝叶斯估计法来估计每个小区间的组频率概率密度,获得组频率的置信区间,若每个小区间的组频率均在稳定区间内,则认为速率集合稳定。
一种计算机可读存储介质,其中存储有多条指令,所述指令适于由终端设备的处理器加载并执行所述的基于变化速率的单变量报警方法。
一种终端设备,包括处理器和计算机可读存储介质,处理器用于实现各指令;计算机可读存储介质用于存储多条指令,所述指令适于由处理器加载并执行所述的基于变化速率的单变量报警方法。
与现有技术相比,本公开的有益效果为:
本公开提供一种基于信号变化速率的报警装置,通过贝叶斯估计判断历史数据所获得的速率集合分布稳定。根据分布稳定的数据百分位数获得速率报警阈值的上下限,当系统在线运行时,计算窗口内信号速率的置信区间并与报警阈值相比较,判断系统是否发生异常。克服了传统的单变量报警系统大多通过监测模拟信号的幅值变化来判断系统是否发生异常,当信号幅值超出报警阈值则触发报警,容易造成误判的问题。
形成的报警系统不再局限于传统报警系统只监测模拟信号幅值变化的模式,对信号速率的变化进行监控,尤其对于速率变化较为平稳的信号,报警系统的性能要明显优于传统报警系统,并且可以得到信号变化趋势。
本公开提供的报警系统和报警方法,可以应用至灾情预警、气象预警、工业变量监测预警等诸多方面,在工业上具有很好的应用前景。
附图说明
构成本公开的一部分的说明书附图用来提供对本公开的进一步理解,本公开的示意性实施例及其说明用于解释本公开,并不构成对本公开的不当限定。
图1为本实施例一中稳定速率分布的组频率及其置信区间;
图2为本实施例一中置信区间与组频率;
图3为本实施例一中稳定速率分布;
图4(a)-(b)为本实施例一中速率报警器结果图;
图5是本实施例一种的流程示意图。
具体实施方式:
下面结合附图与实施例对本公开作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本公开提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本公开的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
实施例一
首先,进行名词解释:
所述分段线性表示,是指将一定长度的时间序列X分割成较短的段,且每段用直线表示。
所述报警信号,是指事先给定报警阈值,当模拟信号超过报警阈值时产生数字信号为1,反之为0。
在本实施例中,公开了一种基于变化速率的单变量报警方法,旨在根据有限的历史数据,利用基于贝叶斯估计的方法,判断过程变量的历史速率集合的分布稳定性,并通过稳定的速率集合计算速率报警阈值。利用等长时间窗口,对窗口内的数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围。在当前窗口内的数据完成计算后,将时间窗口滑动一个采样时间,以达到在线监测过程信号速率变化是否异常的目的。
具体包括离线部分和在线部分,如图5所示;
其中离线部分包括:
(1)选取长度为n的历史数据
Figure BDA0002231527970000081
作为训练样本数据;
(2)选定滑动窗口的宽度为L,对每个窗口内的数据采用自底向上法进行分段,并对各段数据进行最小二乘拟合;
(3)采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,并计算每段数据的速率,组成速率集合;
(4)根据得到的速率集合判断是否稳定,若稳定,通过分位数确定速率报警阈值,若不稳定增加样本,重复步骤(2)-(4);
在线部分包括:
(1)选定滑动窗口的宽度为L,对每个窗口内的实时采集的数据采用自底向上法进行分段,并对各段数据进行最小二乘拟合;
(2)采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,并计算每段线性回归的回归参数置信区间上下限,与速率报警阈值进行比较;
(3)根据每段线性回归的回归参数置信区间上下限与速率报警阈值比较的结果判断系统是否发生异常,如果发生异常则报警。
本实施例上述方法中离线部分:
1.选定滑动窗口的宽度为L,对每个窗口内的数据采用自底向上法进行分段,并对各段数据进行最小二乘拟合:
在窗口SW(i0)(1≤i0≤n-L+1)中,即选取数据段
Figure BDA0002231527970000082
记w=SW(i0),将该窗口内的长度为N、起始点为is(i0≤is≤i0+L-1)的第s段数据用线性方程表示:
Figure BDA0002231527970000091
解得系数为:
Figure BDA0002231527970000093
其中,
Figure BDA0002231527970000094
Figure BDA0002231527970000095
2.采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,并计算每段数据的速率,组成速率集合:
将窗口w内的数据段划分为K(1≤K≤Kmax())段数据段,将其表示为
Figure BDA0002231527970000096
每段数据的结果用线性方程表示:
Figure BDA0002231527970000097
每个窗口内的最优分段数采用均方根误差算法进行确定,窗口w内分段数Kw与加权均方根误差所构成的损失函数可表示为:
Figure BDA0002231527970000098
其中,
Figure BDA0002231527970000099
为窗口w内第g段数据的数据长度;
Figure BDA00022315279700000910
Figure BDA00022315279700000911
当均方根误差最小时便可得到窗口w内数据的最优分段数
Figure BDA00022315279700000912
即:
Figure BDA00022315279700000913
根据最优分段数结果,对每段数据进行一阶最小二乘拟合,每段数据拟合直线的斜率即为的该段的变化速率,所有速率的集合即为历史速率集合{ktrain};
3.判断得到的速率集合是否稳定,若稳定,通过分位数确定速率报警阈值;
通过正常状况下历史数据{xtrain}得到的历史速率集合Ktrain:={ktrain},将{ktrain}的最大值和最小值等长度间隔分为B个小区间,设第b(1≤b≤B)个小区间内包含m个数据,这里通过贝叶斯估计法来估计每个小区间的组频率概率密度,获得组频率的置信区间;
设第b个区间的组频率为:
Figure BDA0002231527970000101
且fb~U(0,1),即离散概率密度函数为:
Figure BDA0002231527970000102
其中,M为离散点的个数,则历史速率集合{ktrain}在fb的条件分布下,似然函数为:
Figure BDA0002231527970000103
故联合分布的概率密度函数为:
Figure BDA0002231527970000104
所以,在k=m时,根据贝叶斯公式,后验分布的概率密度函数为:
Figure BDA0002231527970000105
求出此时后验分布的1-α的置信区间[fb,low,fb,high],其中α为显著性水平。故置信区间的宽度为:
W(b)=fb,high-fb,low
定义组频率稳定指标:
Figure BDA0002231527970000111
设Dth为组频率稳定阈值,当D(b)>Dth时,则认为第b组的组频率稳定;反之,则第b组的组频率不稳定。
Figure BDA0002231527970000112
Figure BDA0002231527970000113
个区间均满足D(b)>Dth,且这
Figure BDA0002231527970000114
个区间的概率密度函数所围成的面积S≥β时,记为样本S0,并继续等长度增加速率样本容量nadd,即样本容量nk=nk+nadd,构成样本SN,样本SN在S0的置信区间内的组频率之和大于β,即可判断分布已稳定,此时选择S0作为稳定分布;反之,则应继续增加样本容量nk=nk+nadd,直至满足分布稳定条件,如图1-3所示。
在线部分:
1.根据离线部分第2步在每个窗口内采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,得到每段线性模型:
Figure BDA0002231527970000115
其中,用
Figure BDA0002231527970000116
Figure BDA0002231527970000117
来表示线性模型的系数,
Figure BDA0002231527970000118
为数据通过一阶最小二乘拟合所得到的估计值;
2.计算窗口w内第q段线性回归的回归参数置信区间上下限,并与速率报警阈值相比较,即置信区间上下限为:
Figure BDA0002231527970000119
Figure BDA00022315279700001110
其中,
Figure BDA00022315279700001111
其中,il表示窗口w内第l段数据的开始位置;
Figure BDA00022315279700001112
为窗口w内第q段数据的样本容量;
Figure BDA0002231527970000121
表示时间常数的平均值;
Figure BDA0002231527970000122
则表示显著性水平为α、自由度为
Figure BDA0002231527970000123
的t统计量;
3.定义报警信号为xa,当系统未发生异常时,xa=0;当系统发生异常时,xa=1,表示为:
Figure BDA0002231527970000124
并进行报警。如图4(a)、(b)所示。
实施例二
提供一种基于变化速率的单变量报警系统,包括在线单元、离线单元和报警单元,其中:
所述离线单元,被配置为根据分布稳定的历史变化速率样本集合,确定速率信号的报警阈值;
所述在线单元,被配置为实时采集被监测的单变量数据,利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与所述离线单元的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警指令;
所述报警单元,被配置为接收所述报警指令,并发出报警信号。
实施例三:
一种计算机可读存储介质,其中存储有多条指令,所述指令适于由终端设备的处理器加载并执行所述的基于变化速率的单变量报警方法。
实施例四:
一种终端设备,包括处理器和计算机可读存储介质,处理器用于实现各指令;计算机可读存储介质用于存储多条指令,所述指令适于由处理器加载并执行所述的基于变化速率的单变量报警方法。
实施例五:
一种基于变化速率的发电厂磨煤机轴承温度报警系统,包括在线单元、离线单元和报警单元,其中:
所述离线单元,被配置为根据无报警情况下分布稳定的历史磨煤机轴承温度数据的变化速率样本集合,确定速率信号的报警阈值;
所述在线单元,被配置为获取温度传感器实时采集的磨煤机轴承温度信息,利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与所述离线单元的变化速率报警阈值相比较,判断温度信号的速率变化是否超出正常的速率变化范围,超出时生成报警指令;
所述报警单元,被配置为接收所述报警指令,并发出报警信号,提醒工作人员磨煤机可能处于异常工作状态,如果不及时采取措施处理,可能会导致磨煤机轴承损坏影响正常生产。
一种基于变化速率的磨煤机轴承温度报警方法,包括以下步骤:
根据实时采集轴承温度数据,利用等长时间窗口,对窗口内的轴承温度数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与预先设定的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警信号。
预先设定的速率报警阈值通过以下过程得到:
(1)选取一定长度的无报警情况下磨煤机轴承温度历史数据作为训练样本数据;
(2)选定滑动窗口的宽度,对每个窗口内的数据采用自底向上法进行分段,并对各段数据进行最小二乘拟合;
(3)采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,并计算每段数据的变化速率,组成速率集合;
(4)根据得到的速率集合判断是否稳定,若稳定,通过分位数确定速率报警阈值,若不稳定,则增加样本长度,重复步骤(2)-(4)。
当然,在其他实施例中,发电厂磨煤机轴承温度报警系统基于的单变量参数可以替换为主蒸汽压力等其他变量。基本过程依据变量的改变可以进行适应性的调整。
当然,在其他实施例中,可以根据被监测对象的不同,对单变量参数进行适应性调整,例如,当被监测对象是基坑时,可以替换为水位参数的报警,当被监测对象是火电机组时,可以替换为功率参数的报警等等。在此不再赘述。
本领域内的技术人员应明白,本公开的实施例可提供为方法、系统、或计算机程序产品。因此,本公开可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本公开可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本公开是参照根据本公开实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本公开的优选实施例而已,并不用于限制本公开,对于本领域的技术人员来说,本公开可以有各种更改和变化。凡在本公开的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。
上述虽然结合附图对本公开的具体实施方式进行了描述,但并非对本公开保护范围的限制,所属领域技术人员应该明白,在本公开的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。

Claims (6)

1.一种基于变化速率的单变量报警系统,其特征是:包括在线单元、离线单元和报警单元,其中:
所述离线单元,被配置为根据分布稳定的历史变化速率样本集合,确定速率信号的报警阈值;
所述在线单元,被配置为实时采集被监测的单变量数据,利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与所述离线单元的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警指令;
所述报警单元,被配置为接收所述报警指令,并发出报警信号;
所述离线单元包括:
初步选择模块,被配置为选取一定长度历史数据作为训练样本数据;
分段模块,被配置为选定滑动窗口宽度,对每个窗口内的数据采用自底向上法进行分段,并对各段数据进行最小二乘拟合;
计算模块,被配置为采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,并计算每段数据的变化速率,组成速率集合;
处理模块,被配置为判断速率集合是否稳定,若稳定,通过分位数确定速率报警阈值,若不稳定,逐步增加训练样本的选取长度,直到稳定;
所述计算模块,还被配置为每个窗口内的最优分段数采用均方根误差算法进行确定,表示各窗口内分段数与加权均方根误差所构成的损失函数,当所述损失函数最小时,得到对应窗口内数据的最优分段数;计算每段数据的变化速率的具体过程是,根据最优分段数结果,对每段数据进行一阶最小二乘拟合,每段数据拟合直线的斜率即为的该段的变化速率。
2.如权利要求1所述的一种基于变化速率的单变量报警系统,其特征是:所述在线单元,包括:
采集模块,被配置为实时采集被监测的单变量数据;
计算模块,被配置为利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间;
比较判断模块,被配置为将置信区间上下限与所述离线单元的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警指令。
3.一种基于变化速率的单变量报警方法,其特征是:包括以下步骤:
实时采集被监测的单变量数据,利用等长时间窗口,对窗口内的单变量数据进行线性分段表示,计算窗口内数据的最小二乘拟合结果及斜率的置信区间,并与预先设定的速率报警阈值相比较,判断信号的速率变化是否超出正常的速率变化范围,超出时生成报警信号;
预先设定的速率报警阈值通过以下过程得到:
(1)选取一定长度的历史数据作为训练样本数据;
(2)选定滑动窗口的宽度,对每个窗口内的数据采用自底向上法进行分段,并对各段数据进行最小二乘拟合;
(3)采用均方根误差算法在每个窗口内获得分段线性表示的最佳分段数,并计算每段数据的变化速率,组成速率集合;
(4)根据得到的速率集合判断是否稳定,若稳定,通过分位数确定速率报警阈值,若不稳定,则增加样本长度,重复步骤(2)-(4);
所述步骤(3)中,每个窗口内的最优分段数采用均方根误差算法进行确定,表示各窗口内分段数与加权均方根误差所构成的损失函数,当所述损失函数最小时,得到对应窗口内数据的最优分段数;计算每段数据的变化速率的具体过程是,根据最优分段数结果,对每段数据进行一阶最小二乘拟合,每段数据拟合直线的斜率即为的该段的变化速率。
4.如权利要求3所述的一种基于变化速率的单变量报警方法,其特征是:所述步骤(4)中,判断速率集合是否稳定的过程包括:将速率集合的最大值和最小值等长度间隔分为多个小区间,通过贝叶斯估计法来估计每个小区间的组频率概率密度,获得组频率的置信区间,若每个小区间的组频率均在稳定区间内,则认为速率集合稳定。
5.一种计算机可读存储介质,其特征是:其中存储有多条指令,所述指令适于由终端设备的处理器加载并执行权利要求3-4中任一项所述的基于变化速率的单变量报警方法。
6.一种终端设备,其特征是:包括处理器和计算机可读存储介质,处理器用于实现各指令;计算机可读存储介质用于存储多条指令,所述指令适于由处理器加载并执行权利要求3-4中任一项所述的基于变化速率的单变量报警方法。
CN201910969211.9A 2019-10-12 2019-10-12 一种基于变化速率的单变量报警系统及方法 Active CN110648480B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910969211.9A CN110648480B (zh) 2019-10-12 2019-10-12 一种基于变化速率的单变量报警系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910969211.9A CN110648480B (zh) 2019-10-12 2019-10-12 一种基于变化速率的单变量报警系统及方法

Publications (2)

Publication Number Publication Date
CN110648480A CN110648480A (zh) 2020-01-03
CN110648480B true CN110648480B (zh) 2021-01-01

Family

ID=69012542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910969211.9A Active CN110648480B (zh) 2019-10-12 2019-10-12 一种基于变化速率的单变量报警系统及方法

Country Status (1)

Country Link
CN (1) CN110648480B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111259538B (zh) * 2020-01-14 2023-08-01 河北振创电子科技有限公司 一种振动报警方法、存储介质和电子设备
CN111915858B (zh) * 2020-08-04 2022-03-29 山东科技大学 一种融合模拟量与数字量相关信息的报警方法及系统
CN112114576B (zh) * 2020-08-05 2021-07-30 山东科技大学 基于输入输出变化一致性的单反馈控制回路异常监测方法
CN113020560B (zh) * 2021-03-02 2022-10-11 武汉巨沣工程技术有限公司 大包下渣检测系统
CN113486074B (zh) * 2021-06-02 2023-11-03 北方工业大学 相关过程参量延迟时间的获取方法
CN114005256B (zh) * 2021-10-25 2023-06-20 山东科技大学 基于报警持续时间概率函数的报警延迟器设计方法及系统
CN114446033B (zh) * 2022-01-21 2023-06-20 四川省地质工程勘察院集团有限公司 一种动态报警阈值确定方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8854217B2 (en) * 2013-02-13 2014-10-07 Sensus Spectrum Llc Method and apparatus for operating an electricity meter
CN104880247B (zh) * 2015-05-04 2016-01-20 华中科技大学 一种用于旋转机械在线监测系统的组合报警方法
CN105041631B (zh) * 2015-07-21 2017-06-06 莱芜钢铁集团电子有限公司 一种气体压缩机的驱动轴振动信号的检测方法和系统
US9734691B2 (en) * 2015-12-15 2017-08-15 Intel Corporation Personalized wearable gas sensor using both the average and the rate of change of the gas level
CN110057406B (zh) * 2019-05-22 2022-05-17 西安因联信息科技有限公司 一种多尺度自适应的机械设备趋势预警方法

Also Published As

Publication number Publication date
CN110648480A (zh) 2020-01-03

Similar Documents

Publication Publication Date Title
CN110648480B (zh) 一种基于变化速率的单变量报警系统及方法
CN112348237B (zh) 一种动态钻井数据异常趋势检测方法
CN108663995B (zh) 一种工业过程变量趋势异常检测方法及装置
CN110057406B (zh) 一种多尺度自适应的机械设备趋势预警方法
CN112149860A (zh) 一种自动异常检测方法和系统
CN114135477B (zh) 一种机泵设备状态监测动态阈值预警方法
CN117972614B (zh) 一种移动电源充放电过程中的温度实时监测方法
CN112288597A (zh) 一种基于层次聚类和直方图算法的能耗在线异常检测方法
CN108108253A (zh) 一种面向多数据流的异常状态检测方法
CN115220396A (zh) 一种数控机床智能监控方法及系统
CN117422310A (zh) 一种全钢内衬层联动线群的电控决策优化方法
CN117454201B (zh) 一种智慧电网异常运行状态检测方法及系统
CN109211564B (zh) 一种用于滚珠丝杠副健康评估的自适应阈值检测方法
CN111176226A (zh) 一种基于运行工况的设备特征参数报警阈值自动分析方法
CN108508860B (zh) 一种基于耦合关系的流程工业生产系统数据监测方法
CN109087490A (zh) 一种基于报警持续时间特征的干扰报警识别方法
CN110083804B (zh) 基于条件分布回归的风电场scada数据缺失的智能修复方法
CN110207827B (zh) 一种基于异常因子提取的电气设备温度实时预警方法
CN117077051A (zh) 一种大坝监测异常数据自适应识别方法
CN113052272B (zh) 一种异常检测方法、装置、电子设备及存储介质
CN118090078B (zh) 一种闭式循环水冷系统渗漏在线监测方法
CN106291344B (zh) 基于分闸压力曲线的断路器刚分点检测方法及系统
CN118464209B (zh) 一种用于远程温感监测报警的温度数据监测优化方法
CN118442125B (zh) 基于物联网的煤炭矿井安全监测与预警系统
CN118071158B (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