CN106368813A - 一种基于多元时间序列的异常报警数据检测方法 - Google Patents
一种基于多元时间序列的异常报警数据检测方法 Download PDFInfo
- Publication number
- CN106368813A CN106368813A CN201610772979.3A CN201610772979A CN106368813A CN 106368813 A CN106368813 A CN 106368813A CN 201610772979 A CN201610772979 A CN 201610772979A CN 106368813 A CN106368813 A CN 106368813A
- Authority
- CN
- China
- Prior art keywords
- time series
- point
- variable
- segmentation
- multivariate time
- 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.)
- Granted
Links
- 230000002159 abnormal effect Effects 0.000 title claims abstract description 35
- 238000001514 detection method Methods 0.000 title claims abstract description 28
- 230000011218 segmentation Effects 0.000 claims abstract description 52
- 230000002596 correlated effect Effects 0.000 claims abstract description 11
- 238000010219 correlation analysis Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 23
- 230000000875 corresponding effect Effects 0.000 claims description 21
- 230000007306 turnover Effects 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000007689 inspection Methods 0.000 claims description 4
- 238000002372 labelling Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000003556 assay Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims description 2
- 239000012634 fragment Substances 0.000 claims description 2
- 230000002349 favourable effect Effects 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 241001212149 Cathetus Species 0.000 description 1
- 238000003646 Spearman's rank correlation coefficient Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02B—INTERNAL-COMBUSTION PISTON ENGINES; COMBUSTION ENGINES IN GENERAL
- F02B77/00—Component parts, details or accessories, not otherwise provided for
- F02B77/08—Safety, indicating, or supervising devices
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于多元时间序列的异常报警数据检测方法,从历史数据中提取多个相关变量的数据,建立多元时间序列,将其标准化,计算正常状态下每个变量之间的符号方向;确定基于关键转折点的时间序列分段描述,设置最小时间间隔,进行关键转折点搜索;表示多元时间序列的分段线性,根据数据点到各个分段的正交距离确定拟合误差,设置损失函数阈值,优化分段数量,得到优化后的分段结果;基于优化之后分段结果,对多元时间序列的各个分段进行相关性分析,提取各个分段变量之间的符号方向,检测与正常状态下的符号方向不一致的异常数据。本发明为实现多变量报警系统的动态报警阈值设计提供有利的条件,从而减少干扰报警。
Description
技术领域
本发明涉及一种基于多元时间序列的异常报警数据检测方法。
背景技术
报警系统对保障燃煤发电机组的安全生产与高效运行发挥着至关重要的作用,由于实际工业过程中关联变量之间的相互影响,传统的单变量报警阈值设计方法可能产生大量干扰报警(漏报警和误报警)并导致“报警过多”的发生,使得现场操作人员的注意力受到影响,增大了在异常生产状况发生时做出正确处置的难度。为了实现多变量报警系统的动态报警阈值设计,寻找一种从历史数据中自动筛选出处于正常状况和异常状况数据段的检测方法是十分必要的。
发明内容
本发明为了解决上述问题,提出了一种基于多元时间序列的异常报警数据检测方法,本方法通过建立多元时间序列,从模式异常的角度提出了一种结合时间序列分段线性表示方法和定性趋势分析方法的异常数据检测算法,自动对历史数据进行异常数据检测,为实现多变量报警系统的动态报警阈值设计提供有利的条件,从而减少干扰报警,提高现场操作人员处理报警的效率,保障了生产安全性。
为了实现上述目的,本发明采用如下技术方案:
一种基于多元时间序列的异常报警数据检测方法,包括以下步骤:
(1)从历史数据中提取当前工作点之前的一定时间内的多个相关变量的数据,建立多元时间序列,将其标准化,计算正常状态下每个变量之间的符号方向;
(2)确定基于关键转折点的时间序列分段描述,设置最小时间间隔,进行关键转折点搜索;
(3)基于线性分段的多元时间序列,根据数据点到各个分段的正交距离确定拟合误差,设置损失函数阈值,优化分段数量,得到优化后的分段结果;
(4)基于优化之后分段结果,对多元时间序列的各个分段进行相关性分析,提取各个分段变量之间的符号方向,检测与正常状态下的符号方向不一致的异常数据。
所述步骤(1)中,具体包括:
(1-1)提取当前工作点之前时间长度为n的多个相关变量的原始数据,建立多元时间序列;
(1-2)求取原始数据的样本均值和标准差,将多元时间序列标准化;
(1-3)根据每两个变量的相关系数确定符号方向,构建符号方向矩阵。
所述步骤(1-3)中,利用ρΔT[Xi,Xj]代表两时刻之间的子时间序列内变量Xi和Xj的相关系数,任意两个变量Xi和Xj在同一段子时间序列内的符号方向signΔT(Xi,Xj)取值1,-1,0分别表示变量之间关系为正相关、负相关、无显著相关。
所述步骤(2)中,具体包括:
(2-1)给出基于关键转折点的时间序列分段的数学描述,即多元时间序列被K+1个关键点分为K个不重叠的时间片段;
(2-2)给出由m个变量和时间t构成的m+1维线性空间中点到直线的正交距离的数学描述;
(2-3)设置最小时间间隔,以其作为搜索关键转折点过程的停止条件。
所述步骤(3)中,利用在最小时间间隔δ的约束下得到的关键转折点集合对原始时间序列进行分段线性表示,选择恰当的分段数量以避免过拟合,实现在一定拟合误差的约束下使用较少的关键转折点作为最终的分段点。
所述步骤(3)中,具体包括:
(3-1)利用线性插值的方法,将多维空间中的数据点在所属分段首尾数据点连线上的投影作为拟合点,并将多元时间序列进行分段线性表示;
(3-2)利用正交距离表示拟合误差;
(3-3)将拟合误差作为损失函数,将不同分段数K对应的损失函数值E(K)绘制到平面直角坐标系中,选取损失函数值小于阈值对应的点所对应的分段数为优化结果。
所述步骤(3-3)中,观察损失函数值随分段数K的收敛情况,在损失函数值减小随K值增加趋于平稳的区域设置合理的损失函数阈值。
所述步骤(3-3)中,设第一个满足损失函数值小于损失函数阈值对应的点坐标为(c,E(c)),则选择该点对应的分段数K=c作为优化结果和对应的关键转折点集合Qc=[q1,…,qc+1]作为用于进行相关性趋势提取的分段点。
所述步骤(4)中,具体包括:
(4-1)计算每个分段中任意两个变量之间的相关系数;
(4-2)对变量相关性进行单边假设检验,设置显著性水平,根据单边假设检验结果和显著性水平确认变量间的相关性,确定变量符号;
(4-3)根据变量符号,构建变量符号矩阵,将变量符号矩阵和符号方向矩阵中对应位置的元素进行比较,如两者不同,则其对应的分段为异常数据。
所述步骤(4-1)中,具体方法为:
假设将时间序列T的第s个分段内变量Xi和Xj的观测值分别从小到大排列,并依次用正整数k=1,…,zs标记,记为Rk和Sk,则时间序列T中第s个分段内任意两个变量之间的Spearman样本相关系数为:
其中dk=Rk-Sk。
所述步骤(4-2)中,单边假设检验:H0:ρs[Xi,Xj]=0vs H1:ρs[Xi,Xj]>0,H0:ρs[Xi,Xj]=0vs H2:ρs[Xi,Xj]<0;
当样本个数n>10时,随机变量Us被定义为:其中,zs为第s个分段内的样本个数;给定显著性水平α,如果Us>tα(zs-2),则与H1相对的H0被拒绝,如果Us<-tα(zs-2),则与H2相对的H0被拒绝,其中tα(zs-2)表示统计量Us的分位数,此时,第s个分段内Xi和Xj的相关性被认为是显著的,符号方向signs(Xi,Xj)分别取值为1或-1,如果|Us|<tα(zs-2),无论对于H1或者H2,H0都不能被拒绝,此时变量间无显著相关性,符号方向signs(Xi,Xj)取值为0。
所述步骤(4-2)中,当样本个数n<10时,查询用于小样本容量假设检验的Spearman秩相关系数的临界值,将对应于给定zs和α的相关系数临界值表示为ρα(zs),如果|ρs[Xi,Xj]|>ρα(zs),H0被拒绝,signs(Xi,Xj)分别取值为1或-1,反之H0不能被拒绝,符号方向signs(Xi,Xj)取值为0。
本发明的有益效果为:本发明选取工业变量之间的相关性作为判断工作点状态是否异常的特征,通过建立多元时间序列,从模式异常的角度提出了一种结合时间序列分段线性表示 方法和定性趋势分析方法的异常数据检测算法,自动对历史数据进行异常数据检测,为实现多变量报警系统的动态报警阈值设计提供有利的条件,从而减少干扰报警,提高现场操作人员处理报警的效率,保障了生产安全性。
附图说明
图1为本发明所述基于工业历史数据的报警系统异常数据检测方法流程图;
图2为本发明具体实施例中变量时间序列和分段结果图;
图3为本发明具体实施例中分段数K的决策图;
图4(a)为变量在每个分段中相关性分析结果;
图4(b)为变量之间符号方向计算结果;
图5(a)为用不同线段表示的异常数据检测结果;
图5(b)为用不同数值表示的异常数据检测结果。
具体实施方式:
下面结合附图与实施例对本发明作进一步说明。
图1为本发明所述基于工业历史数据的报警系统异常数据检测方法流程图。
如图1所示,一种基于工业历史数据的报警系统异常数据检测方法,包括如下步骤:
步骤S1,从历史数据中提取当前工作点之前的时间t内的多个相关变量的数据,建立多元时间序列T’,并将其标准化为时间序列T,同时计算正常状态下各个变量之间的符号方向;
步骤S2,设置最小时间间隔δ,并基于多元时间序列T进行关键转折点搜索;
步骤S3,基于多元时间序列T线性分段的拟合误差进行分段数K的优化;
步骤S4,基于优化之后分段结果,对多元时间序列T的各个分段进行相关性分析,提取各个分段变量之间的符号方向,根据其与正常状态下的符号方向是否一致来检测异常数据。
在本发明的具体实施例中,步骤S1的具体实现为:
步骤S11,提取当前工作点之前时间长度为n的多个相关变量的原始数据,用X(t)来表示变量i在t时刻的数值,建立多元时间序列T′={X′i(t)},其中i=1,…,m,t=1,…,n,m表示变量个数,n表示时间长度。
步骤S12,将多元时间序列T’标准化为T,其中代表原始数据X′i(t)的样本均值,代表原始数据X′i(t)的样本标准差。
步骤S13,计算任意两个变量Xi和Xj在同一段子时间序列内的符号方向: 其中ΔT代表时刻t1到时刻t2的子时间序列,1≤t1≤t2≤n,ρΔT[Xi,Xj]代表ΔT内变量Xi和Xj的相关系数,signΔT(Xi,Xj)取值1,-1,0分别表示变量之间关系为正相关、负相关、无显著相关。在正常情况下,变量之间的符号方向保持不变,符号方向矩阵可被定义为:
在本发明的具体实施例中,步骤S2的具体实现为:
步骤S21,给出基于关键转折点的时间序列分段的数学描述。给定整数K,时间序列T可以被K+1个关键点分为K个不重叠的的时间片段,用S={p1,…,pK+1}表示,其中pi,i=1,…,k+1,代表第i个关键转折点的时间标签,并且有p1=1,pK+1=n。将S中第j个分段表示为sj={{Xi(t)},pj<t≤pj+1},其中j=1,…,K。定义zj=pj+1-pj为第j个分段中包含的数据点的个数。
步骤S22,给出由m个变量和时间t构成的m+1维线性空间中点到直线的正交距离的数学描述。空间中直线AB的参数方程可表示为:其中i=1,…,m。则直线AB上任意一点P0的坐标可表示为[(XiB-XiA)β+XiA,(tB-tA)β+tA]。因此,点P到直线AB的距离可被定义为: 其中, 当取极小值时对应的参数则点P到直线AB的最小距离即正交距离为
步骤S23,设置最小时间间隔δ,0<δ<n,用于减少噪音对分段结果的影响,并作为为搜索关键转折点过程的停止条件。在处理实际工业过程数据时,噪声的干扰会导致关键点之间的时间间隔过短,因此的当满足条件min(zj)<δ,j=1,…,K时停止搜索关键转折点。
在本发明的具体实施例中,步骤S3的具体实现为:
利用在最小时间间隔δ的约束下得到的关键转折点集合Q对原始时间序列进行分段线性表示,通常会导致过拟合,为了避免过拟合,有必要选择合适的分段数K,以实现在一定拟合误差的约束下使用较少的关键转折点作为最终的分段点。
步骤S31,定义多元时间序列的分段线性表示。当多元时间序列T被K+1个关键转折点p1,…,pK+1分为K段,则此时多元时间序列T的分段线性表示为:TPLR=<f1[(Xi(p1),p1),(Xi(p2),p2)],…,fK[(Xi(pK),pK),(Xi(pK+1),pK+1)]>。其中f1[(Xi(p1),p1),(Xi(p2),p2)]表示在分段[pj,pj+1]内的线性拟合函数。本发明利用线性插值的方法,将m+1维空间中的数据点(Xi(t),t),其中i=1,…,m,t=1,…,n在所属分段首尾数据点连线上的投影作为拟合点,从而得到拟合点
步骤S32,定义时间序列分段线性表示的拟合误差。对时间序列T进行分段线性表示,采用线性插值得到原始数据的拟合点,则拟合误差为:其中i=1,…,m。根据线性空间中点到直线正交距离的定义,拟合误差还可以用正交距离表示为: 其中D(j)表示数据点(Xi(j),j)到所属分段的正交距离。
步骤S33,设置损失函数阈值η,对分段数K进行优化并得到优化的分段结果。将拟合误差E作为损失函数,假定关键转折点集合Q=[q1,q2,…,ql,ql+1],计算当2≤K≤l时,不同分段数K对应的损失函数值E(K),并将所有点(K,E(K))绘制到平面直角坐标系中。
观察损失函数E随分段数K的收敛情况,在E值减小随K值增加趋于平稳的区域设置合理的损失函数阈值η。假设第一个满足E(K)<η对应的点坐标为(c,E(c)),则选择该点对应的分段数K=c作为优化结果和对应Q中的关键转折点集合Qc=[q1,…,qc+1]作为用于进行相关性趋势提取的分段点。
在本发明的具体实施例中,步骤S4的具体实现为:
步骤S41,获得任意两个变量之间的相关系数。假设将时间序列T的第s个分段内变量Xi和Xj的观测值分别从小到大排列,并依次用正整数k=1,…,zs标记,记为Rk和Sk。则时间序列T中第s个分段内任意两个变量之间的Spearman样本相关系数为: 其中dk=Rk-Sk。
步骤S42,对变量相关性进行单边假设检验。单边假设检验:H0:ρs[Xi,Xj]=0vsH1:ρs[Xi,Xj]>0,H0:ρs[Xi,Xj]=0vs H2:ρs[Xi,Xj]<0。当样本个数n>10时,随机变量Us可被定义为:其中,zs为第s个分段内的样本个数。给定显著性水平α,如果Us>tα(zs-2),则与H1相对的H0被拒绝,如果Us<-tα(zs-2),则与H2相对的H0被拒绝,其中tα(zs-2)表示统计量Us的分位数。此时,第s个分段内Xi和Xj的相关性被认为是显著的,符号方向signs(Xi,Xj)分别取值为1或-1。如果|Us|<tα(zs-2),无论对于H1或者H2,H0都不能被拒绝,此时变量间无显著相关性,符号方向signs(Xi,Xj)取值为0。
当样本个数n<10时,查询用于小样本容量假设检验的Spearman秩相关系数的临界值,将对应于给定zs和α的相关系数临界值表示为ρα(zs),如果|ρs[Xi,Xj]|>ρα(zs),H0被拒绝,signs(Xi,Xj)分别取值为1或-1,反之H0不能被拒绝,符号方向signs(Xi,Xj)取值为0。
步骤S43,根据异常数据检测规则判定异常数据。如果对于所有的i,j∈[1,2,…,m],均满足signs(Xi,Xj)=signT(Xi,Xj),则第s分段被划分为正常数据。如果存在任一i,j∈[1,2,…,m],使得signs(Xi,Xj)≠signT(Xi,Xj),则第s分段被划分为异常数据。
以下是本发明所述方法在具体示例中的应用,具体应用场景为电厂。
选定电厂中的给水泵的进口流量和给水泵汽轮机转速作为相关变量,在电厂中的一次停机事故中,选取停机前采样周期为1秒,样本容量为n=2239的流量和转速二元时间序列作为历史数据,并将其标准化,标准化之后的二元时间序列记为T=[Q(t),V(t)],其中t=1,…,n。
根据其工作原理,正常工况下给水泵进口流量(简称流量,用Q表示)与给水泵汽轮机转速(简称转速,用V表示)保持高度正相关性,因此符号方向signT(Q,V)=1,符号方向矩阵
给定最小时间间隔δ=15进行关键转折点搜索,当搜索停止时,K=116,损失函数阈值与K的分布如图2所示,根据观察,选择η=0.3作为损失函数阈值,得到优化后的分段数K=28。
根据优化后获得的分段数,绘制变量的分段时间序列图,如图3所示。其中实线代表的是各个关键转折点的位置,虚线代表相关性发生显著变化的时间段。
给定α=0.05,计算每一个分段中Q和V的样本相关系数并进行相关性检验,流量Q与转速V在每个分段中相关性分析结果如图4(a)所示,流量Q与转速V的符号方向计算结果如图4(b)所示,从图中可以看出第24,25和28分段中变量相关性发生了显著变化,即此时signs(Q,V)≠signT(Q,V)。
根据相关性趋势分析结果可知,t=1911-2126,以及t=2195-2239之间的数据被检测为异常数据,图5(a)用不同的线段标记了异常数据的检测结果,这些数据所在分段与相关性发生显著变化的第24,25和28分段是一致的,图5(b)用不同数值标记异常数据检测结果。分析可知,第一部分数据出现异常是由于位于给水泵下游的汽包压力异常升高,导致给水泵进出口压差减小,阻力增大,出现转速升高但流量却下降的现象;第二部分数据出现异常是由于前期异常没有得到及时正确处置触发了机组应急停车,这部分数据对于报警系统设计以及故障原因分析的作用十分有限。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.一种基于多元时间序列的异常报警数据检测方法,其特征是:包括以下步骤:
(1)从历史数据中提取当前工作点之前的一定时间内的多个相关变量的数据,建立多元时间序列,将其标准化,计算正常状态下每个变量之间的符号方向;
(2)确定基于关键转折点的时间序列分段描述,设置最小时间间隔,进行关键转折点搜索;
(3)基于线性分段的多元时间序列,根据数据点到各个分段的正交距离确定拟合误差,设置损失函数阈值,优化分段数量,得到优化后的分段结果;
(4)基于优化之后分段结果,对多元时间序列的各个分段进行相关性分析,提取各个分段变量之间的符号方向,检测与正常状态下的符号方向不一致的异常数据。
2.如权利要求1所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(1)中,具体包括:
(1-1)提取当前工作点之前时间长度为n的多个相关变量的原始数据,建立多元时间序列;
(1-2)求取原始数据的样本均值和标准差,将多元时间序列标准化;
(1-3)根据每两个变量的相关系数确定符号方向,构建符号方向矩阵。
3.如权利要求2所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(1-3)中,利用ρΔT[Xi,Xj]代表两时刻之间的子时间序列内变量Xi和Xj的相关系数,任意两个变量Xi和Xj在同一段子时间序列内的符号方向signΔT(Xi,Xj)取值1,-1,0分别表示变量之间关系为正相关、负相关、无显著相关。
4.如权利要求1所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(2)中,具体包括:
(2-1)给出基于关键转折点的时间序列分段的数学描述,即多元时间序列被K+1个关键点分为K个不重叠的时间片段;
(2-2)给出由m个变量和时间t构成的m+1维线性空间中点到直线的正交距离的数学描述;
(2-3)设置最小时间间隔,以其作为搜索关键转折点过程的停止条件。
5.如权利要求1所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(3)中,利用在最小时间间隔δ的约束下得到的关键转折点集合对原始时间序列进行分段线性表示,选择恰当的分段数量以避免过拟合,实现在一定拟合误差的约束下使用较少的关键转折点作为最终的分段点。
6.如权利要求1所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(3)中,具体包括:
(3-1)利用线性插值的方法,将多维空间中的数据点在所属分段首尾数据点连线上的投影作为拟合点,并将多元时间序列进行分段线性表示;
(3-2)利用正交距离表示拟合误差;
(3-3)将拟合误差作为损失函数,将不同分段数K对应的损失函数值E(K)绘制到平面直角坐标系中,选取损失函数值小于阈值对应的点所对应的分段数为优化结果;
所述步骤(3-3)中,观察损失函数值随分段数K的收敛情况,在损失函数值减小随K值增加趋于平稳的区域设置合理的损失函数阈值。
7.如权利要求6所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(3-3)中,设第一个满足损失函数值小于损失函数阈值对应的点坐标为(c,E(c)),则选择该点对应的分段数K=c作为优化结果和对应的关键转折点集合Qc=[q1,…,qc+1]作为用于进行相关性趋势提取的分段点。
8.如权利要求1所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(4)中,具体包括:
(4-1)计算每个分段中任意两个变量之间的相关系数;
(4-2)对变量相关性进行单边假设检验,设置显著性水平,根据单边假设检验结果和显著性水平确认变量间的相关性,确定变量符号;
(4-3)根据变量符号,构建变量符号矩阵,将变量符号矩阵和符号方向矩阵中对应位置的元素进行比较,如两者不同,则其对应的分段为异常数据。
9.如权利要求8所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(4-1)中,具体方法为:
假设将时间序列T的第s个分段内变量Xi和Xj的观测值分别从小到大排列,并依次用正整数k=1,…,zs标记,记为Rk和Sk,则时间序列T中第s个分段内任意两个变量之间的Spearman样本相关系数为:
其中dk=Rk-Sk。
10.如权利要求8所述的一种基于多元时间序列的异常报警数据检测方法,其特征是:所述步骤(4-2)中,单边假设检验:H0:ρs[Xi,Xj]=0 vs H1:ρs[Xi,Xj]>0,H0:ρs[Xi,Xj]=0vs H2:ρs[Xi,Xj]<0;
当样本个数n>10时,随机变量Us被定义为:其中,zs为第s个分段内的样本个数;给定显著性水平α,如果Us>tα(zs-2),则与H1相对的H0被拒绝,如果Us<-tα(zs-2),则与H2相对的H0被拒绝,其中tα(zs-2)表示统计量Us的分位数,此时,第s个分段内Xi和Xj的相关性被认为是显著的,符号方向signs(Xi,Xj)分别取值为1或-1,如果|Us|<tα(zs-2),无论对于H1或者H2,H0都不能被拒绝,此时变量间无显著相关性,符号方向signs(Xi,Xj)取值为0。
所述步骤(4-2)中,当样本个数n<10时,查询用于小样本容量假设检验的Spearman秩相关系数的临界值,将对应于给定zs和α的相关系数临界值表示为ρα(zs),如果|ρs[Xi,Xj]|>ρα(zs),H0被拒绝,signs(Xi,Xj)分别取值为1或-1,反之H0不能被拒绝,符号方向signs(Xi,Xj)取值为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610772979.3A CN106368813B (zh) | 2016-08-30 | 2016-08-30 | 一种基于多元时间序列的异常报警数据检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610772979.3A CN106368813B (zh) | 2016-08-30 | 2016-08-30 | 一种基于多元时间序列的异常报警数据检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106368813A true CN106368813A (zh) | 2017-02-01 |
CN106368813B CN106368813B (zh) | 2018-09-25 |
Family
ID=57902386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610772979.3A Active CN106368813B (zh) | 2016-08-30 | 2016-08-30 | 一种基于多元时间序列的异常报警数据检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106368813B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106778053A (zh) * | 2017-03-31 | 2017-05-31 | 国网山东省电力公司电力科学研究院 | 一种基于相关性的报警关联变量检测方法及系统 |
CN107092654A (zh) * | 2017-03-20 | 2017-08-25 | 山东科技大学 | 基于均值变化检测的报警正常与异常数据检测方法和装置 |
CN108320810A (zh) * | 2018-04-11 | 2018-07-24 | 平安科技(深圳)有限公司 | 疾病异常数据检测方法及装置、计算机装置及存储介质 |
CN108615017A (zh) * | 2018-04-28 | 2018-10-02 | 山东科技大学 | 滑动窗分段直线表示法的总拟合误差阈值估计方法与系统 |
CN109444232A (zh) * | 2018-12-26 | 2019-03-08 | 苏州同阳科技发展有限公司 | 一种多通道智能化污染气体监测装置与扩散溯源方法 |
CN109543083A (zh) * | 2018-11-19 | 2019-03-29 | 国网陕西省电力公司电力科学研究院 | 一种多元电网实时数据中异常数据的检测方法 |
CN110298767A (zh) * | 2019-05-29 | 2019-10-01 | 山东科技大学 | 一种火力发电厂时间序列变量异常监测方法及系统 |
WO2019214367A1 (zh) * | 2018-05-07 | 2019-11-14 | 山东科技大学 | 一种动态报警阈值设计和报警消除的方法与系统 |
CN110955648A (zh) * | 2019-12-18 | 2020-04-03 | 重庆大学 | 一种非等距时间序列监测数据的归一化映射处理方法 |
CN111712771A (zh) * | 2017-11-10 | 2020-09-25 | Abb瑞士股份有限公司 | 能够执行具有多个机器人的生产系统中的问题诊断的数据处理装置以及方法 |
CN112348644A (zh) * | 2020-11-16 | 2021-02-09 | 上海品见智能科技有限公司 | 一种通过建立单调正相关过滤网的异常物流订单检测方法 |
CN113449626A (zh) * | 2021-06-23 | 2021-09-28 | 中国科学院上海高等研究院 | 隐马尔科夫模型振动信号分析方法装置、存储介质和终端 |
CN114236272A (zh) * | 2021-12-02 | 2022-03-25 | 深圳市环球众一科技有限公司 | 一种电子产品的智能检测系统 |
CN115858894A (zh) * | 2023-02-14 | 2023-03-28 | 温州众成科技有限公司 | 一种可视化的大数据分析方法 |
CN117370329A (zh) * | 2023-12-07 | 2024-01-09 | 湖南易比特大数据有限公司 | 基于工业物联网的设备数据智能化管理方法及系统 |
CN117814805A (zh) * | 2024-03-05 | 2024-04-05 | 自贡市第一人民医院 | 一种临床护理设备数据智能处理方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101713395A (zh) * | 2009-11-20 | 2010-05-26 | 西安交通大学 | 基于增量式非线性流形学习的喘振监测方法 |
JP5186322B2 (ja) * | 2008-09-26 | 2013-04-17 | インターナショナル・ビジネス・マシーンズ・コーポレーション | 時系列データ解析システム、方法及びプログラム |
CN104019000A (zh) * | 2014-06-23 | 2014-09-03 | 宁夏银星能源股份有限公司 | 风力发电机组的载荷谱测定与前瞻性维护系统 |
US20150235139A1 (en) * | 2014-02-20 | 2015-08-20 | Nec Laboratories America, Inc. | Extracting interpretable features for classification of multivariate time series from physical systems |
CN105241669A (zh) * | 2015-09-09 | 2016-01-13 | 哈尔滨工业大学 | 基于比较编码的燃气轮机燃烧系统在线监测方法 |
CN105761407A (zh) * | 2016-01-06 | 2016-07-13 | 深圳大学 | 基于无线网络信号传输的室内探测火情及报警方法及系统 |
-
2016
- 2016-08-30 CN CN201610772979.3A patent/CN106368813B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5186322B2 (ja) * | 2008-09-26 | 2013-04-17 | インターナショナル・ビジネス・マシーンズ・コーポレーション | 時系列データ解析システム、方法及びプログラム |
CN101713395A (zh) * | 2009-11-20 | 2010-05-26 | 西安交通大学 | 基于增量式非线性流形学习的喘振监测方法 |
US20150235139A1 (en) * | 2014-02-20 | 2015-08-20 | Nec Laboratories America, Inc. | Extracting interpretable features for classification of multivariate time series from physical systems |
CN104019000A (zh) * | 2014-06-23 | 2014-09-03 | 宁夏银星能源股份有限公司 | 风力发电机组的载荷谱测定与前瞻性维护系统 |
CN105241669A (zh) * | 2015-09-09 | 2016-01-13 | 哈尔滨工业大学 | 基于比较编码的燃气轮机燃烧系统在线监测方法 |
CN105761407A (zh) * | 2016-01-06 | 2016-07-13 | 深圳大学 | 基于无线网络信号传输的室内探测火情及报警方法及系统 |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107092654A (zh) * | 2017-03-20 | 2017-08-25 | 山东科技大学 | 基于均值变化检测的报警正常与异常数据检测方法和装置 |
CN107092654B (zh) * | 2017-03-20 | 2019-02-15 | 山东科技大学 | 基于均值变化检测的报警正常与异常数据检测方法和装置 |
CN106778053A (zh) * | 2017-03-31 | 2017-05-31 | 国网山东省电力公司电力科学研究院 | 一种基于相关性的报警关联变量检测方法及系统 |
CN106778053B (zh) * | 2017-03-31 | 2019-04-09 | 国网山东省电力公司电力科学研究院 | 一种基于相关性的报警关联变量检测方法及系统 |
CN111712771A (zh) * | 2017-11-10 | 2020-09-25 | Abb瑞士股份有限公司 | 能够执行具有多个机器人的生产系统中的问题诊断的数据处理装置以及方法 |
CN111712771B (zh) * | 2017-11-10 | 2024-03-12 | Abb瑞士股份有限公司 | 能够执行问题诊断的数据处理装置以及方法 |
CN108320810A (zh) * | 2018-04-11 | 2018-07-24 | 平安科技(深圳)有限公司 | 疾病异常数据检测方法及装置、计算机装置及存储介质 |
CN108615017A (zh) * | 2018-04-28 | 2018-10-02 | 山东科技大学 | 滑动窗分段直线表示法的总拟合误差阈值估计方法与系统 |
WO2019214367A1 (zh) * | 2018-05-07 | 2019-11-14 | 山东科技大学 | 一种动态报警阈值设计和报警消除的方法与系统 |
CN109543083A (zh) * | 2018-11-19 | 2019-03-29 | 国网陕西省电力公司电力科学研究院 | 一种多元电网实时数据中异常数据的检测方法 |
CN109444232B (zh) * | 2018-12-26 | 2024-03-12 | 苏州同阳科技发展有限公司 | 一种多通道智能化污染气体监测装置与扩散溯源方法 |
CN109444232A (zh) * | 2018-12-26 | 2019-03-08 | 苏州同阳科技发展有限公司 | 一种多通道智能化污染气体监测装置与扩散溯源方法 |
CN110298767A (zh) * | 2019-05-29 | 2019-10-01 | 山东科技大学 | 一种火力发电厂时间序列变量异常监测方法及系统 |
CN110955648A (zh) * | 2019-12-18 | 2020-04-03 | 重庆大学 | 一种非等距时间序列监测数据的归一化映射处理方法 |
CN112348644A (zh) * | 2020-11-16 | 2021-02-09 | 上海品见智能科技有限公司 | 一种通过建立单调正相关过滤网的异常物流订单检测方法 |
CN112348644B (zh) * | 2020-11-16 | 2024-04-02 | 上海品见智能科技有限公司 | 一种通过建立单调正相关过滤网的异常物流订单检测方法 |
CN113449626B (zh) * | 2021-06-23 | 2023-11-07 | 中国科学院上海高等研究院 | 隐马尔科夫模型振动信号分析方法装置、存储介质和终端 |
CN113449626A (zh) * | 2021-06-23 | 2021-09-28 | 中国科学院上海高等研究院 | 隐马尔科夫模型振动信号分析方法装置、存储介质和终端 |
CN114236272B (zh) * | 2021-12-02 | 2022-11-08 | 深圳市环球众一科技有限公司 | 一种电子产品的智能检测系统 |
CN114236272A (zh) * | 2021-12-02 | 2022-03-25 | 深圳市环球众一科技有限公司 | 一种电子产品的智能检测系统 |
CN115858894A (zh) * | 2023-02-14 | 2023-03-28 | 温州众成科技有限公司 | 一种可视化的大数据分析方法 |
CN117370329A (zh) * | 2023-12-07 | 2024-01-09 | 湖南易比特大数据有限公司 | 基于工业物联网的设备数据智能化管理方法及系统 |
CN117370329B (zh) * | 2023-12-07 | 2024-02-27 | 湖南易比特大数据有限公司 | 基于工业物联网的设备数据智能化管理方法及系统 |
CN117814805A (zh) * | 2024-03-05 | 2024-04-05 | 自贡市第一人民医院 | 一种临床护理设备数据智能处理方法 |
CN117814805B (zh) * | 2024-03-05 | 2024-06-11 | 自贡市第一人民医院 | 一种临床护理设备数据智能处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106368813B (zh) | 2018-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106368813A (zh) | 一种基于多元时间序列的异常报警数据检测方法 | |
CN101446831B (zh) | 一种分散的过程监测方法 | |
Di Maio et al. | Fault detection in nuclear power plants components by a combination of statistical methods | |
CN110738274A (zh) | 一种基于数据驱动的核动力装置故障诊断方法 | |
Jia et al. | A deviation based assessment methodology for multiple machine health patterns classification and fault detection | |
de la Hermosa González | Wind farm monitoring using Mahalanobis distance and fuzzy clustering | |
CN104699077B (zh) | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 | |
EP2442288A1 (en) | Device abnormality monitoring method and system | |
CN106778053A (zh) | 一种基于相关性的报警关联变量检测方法及系统 | |
Cheng et al. | Monitoring influent conditions of wastewater treatment plants by nonlinear data-based techniques | |
CN105955219A (zh) | 基于互信息的分布式动态过程故障检测方法 | |
CN103776480A (zh) | 基于多次移动平均的微小故障检测方法和装置 | |
CN105404280A (zh) | 基于自回归动态隐变量模型的工业过程故障检测方法 | |
Yu et al. | Abnormal data detection for multivariate alarm systems based on correlation directions | |
Li et al. | Canonical variate residuals-based contribution map for slowly evolving faults | |
Maestri et al. | A robust clustering method for detection of abnormal situations in a process with multiple steady-state operation modes | |
WO2019043600A1 (en) | ESTIMATOR OF REMAINING USEFUL LIFE ESTIMATOR | |
CN110009126B (zh) | 基于pls模型与pca贡献度融合的在线报警分析方法 | |
Zhang et al. | Selecting optimal features for cross-fleet analysis and fault diagnosis of industrial gas turbines | |
CN104615123A (zh) | 基于k近邻的传感器故障隔离方法 | |
Kusiak et al. | Data-mining-based system for prediction of water chemistry faults | |
CN103207567A (zh) | 一种低误报率的改进主元分析过程监测方法及其监测系统 | |
CN117171655A (zh) | 一种基于频谱分析的煤矿生产状态智能判别方法及装置 | |
Martí et al. | YASA: yet another time series segmentation algorithm for anomaly detection in big data problems | |
CN103870711A (zh) | 一种油品调合压力状态识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201010 Address after: 266590 No. 579, Bay Road, Huangdao District, Shandong, Qingdao Patentee after: Shandong University of Science and Technology Address before: 100000, room 3, floor 13, building 3, 304 Feng Xiu Road, Beijing, Haidian District Patentee before: BEIJING XIETONG INNOVATION |