CN108257365B - 一种基于全局不确定性证据动态融合的工业报警器设计方法 - Google Patents
一种基于全局不确定性证据动态融合的工业报警器设计方法 Download PDFInfo
- Publication number
- CN108257365B CN108257365B CN201810082465.4A CN201810082465A CN108257365B CN 108257365 B CN108257365 B CN 108257365B CN 201810082465 A CN201810082465 A CN 201810082465A CN 108257365 B CN108257365 B CN 108257365B
- Authority
- CN
- China
- Prior art keywords
- alarm
- evidence
- global
- otp
- industrial
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 230000004927 fusion Effects 0.000 title claims abstract description 15
- 238000013461 design Methods 0.000 title claims abstract description 13
- 230000002159 abnormal effect Effects 0.000 claims description 19
- 238000005070 sampling Methods 0.000 claims description 11
- 238000012544 monitoring process Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 2
- 230000001149 cognitive effect Effects 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
- G08B29/00—Checking or monitoring of signalling or alarm systems; Prevention or correction of operating errors, e.g. preventing unauthorised operation
- G08B29/18—Prevention or correction of operating errors
- G08B29/185—Signal analysis techniques for reducing or preventing false alarms or for enhancing the reliability of the system
- G08B29/188—Data fusion; cooperative systems, e.g. voting among different detectors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Security & Cryptography (AREA)
- Alarm Systems (AREA)
- Monitoring And Testing Of Exchanges (AREA)
Abstract
本发明公开了一种基于全局不确定性证据动态融合的工业报警器设计方法,属于工业报警器设计技术领域。本发明将基于连续型的Sigmoid隶属度函数和Gaussian型隶属度函数,将采集到的过程变量数据转换成每时刻的报警证据,在这个过程中将全集的不确定性信息赋予一定的数值,更能体现对不确定性证据的处理;采用条件化证据线性更新规则,将当前每时刻的报警证据与前一时刻的全局报警证据进行更新融合,最终得到当前时刻的全局报警证据,再依据相关的决策定理进行报警决策。在更新融合的过程中有效的降低了过程变量的随机不确定性和认知不确定性,提高了报警器的精度。
Description
技术领域
本发明涉及一种基于全局不确定性的证据动态融合的工业报警器的设计方法,属于工业报警器设计技术领域。
背景技术
大型工业工厂由数千个物理对象,传感器,通信网络,执行器和控制回路组成。这些***组件中的任何异常或故障都可能导致故障,非计划停车甚至会导致一些危险事故的发生。作为第一个安全屏障,报警***需要检测工厂中的异常或故障,并生成警报通知操作员。通过了解和确认报警信息,操作人员可以调查报警原因,并采取必要的纠正措施,使其恢复正常运行,避免或减少不必要的损失。
报警产生的最常见做法是将过程变量与阈值(也称为跳闸点,控制极限或报警极限)进行比较。对于操作人员来说,满意的情况是只对任何异常情况发出警报,花费大约10分钟处理警报,即在工厂正常运行期间,每小时不会收到超过6个警报。然而,由于过程变量受到的噪音干扰和不确定性因素的影响。因此,在实践中,操作员遇到的大部分警报都是虚假或误报,包括误报和漏报。在传统的报警器设计的方法中并不能有效的解决不确定性因素对报警的影响,这样报警器的误报率与漏报率并不能有效的降低,因此达不到提高报警器精度的目的。
另外已公开的《一种基于报警证据融合的工业报警器设计方法》发明专利中,只将正常状态和异常状态的报警证据纳入到证据体系中,在过程变量转换中使用的是分段的模糊隶属度函数。
发明内容
本发明针对现有技术的不足,提供了一种基于全局不确定性证据动态融合的工业报警器设计方法。
本发明将全局的不确定性证据也纳入证据体系,利用连续的Sigmoid隶属度函数和Gaussian型隶属度函数,将过程变量的信息无缺失的转换成报警证据,通过归一化处理得到每时刻的报警证据;然后利用条件化的线性证据更新将前一时刻的全局报警证据与当前时刻的报警证据进行更新融合,在更新融合的过程中计算证据的相似度使用的是S型的函数,这一点也区别于上述已公开的发明,最终得到的全局的报警证据,最后基于pignistic概率进行报警决策。
本发明包括以下各步骤:
(1)设定某一工业报警器的辨识框架为θ={NA,A},其中NA表示工业报警器处于正常状态,A表示工业报警器处于异常状态亦即报警状态,θ表示处于全局不确定状态。
(2)设x为工业报警器所检测的过程变量,令x(t),t=1,2,3,….是传感器检测设备的在线序列,t为采样的时刻,采样的个数由工业报警器的监测周期和监控计算机设备来决定,在采集到的序列中,定义xotp为最优的阈值,记max(x(t))和min(x(t))分别是序列的最大值和最小值,其中xotp∈[min(x(t)),max(x(t))]。利用Sigmoid隶属度函数和Gaussian型隶属度函数分别构建正常状态NA、异常状态A和全局状态θ下的模糊隶属度函数μNA(x(t))、μA(x(t))和μθ(x(t)),如式(1)-(3)所示:
其中a是一可调参数,当a取不同值时,隶属度函数会有不同的变化趋势,δ表示高斯型隶属度函数的方差,xotp表示最优的阈值。
(3)将t时刻的过程变量x值x(t)分别带入步骤(2)中关于正常状态NA、异常状态A和全局不确定状态下θ的模糊隶属度函数μNA(x)、μA(x)和μΘ(x)中,即可得到关于报警证据mt(NA)、mt(A)和mt(θ),分别如下(4)-(6)式所示。
将报警证据归一化,令Mt(x(t))=(mNA(x(t))+mA(x(t)+mθ(x(t)),则每一项归一化后的为,mt(NA)=(mNA(x(t))/Mt(x(t)))、mt(A)=(mA(x(t))/Mt(x(t)))、mt(θ)=(mθ(x(t))/Mt(x(t))),则称mt(NA)、mt(A)和mt(θ)是t时刻的报警证据。
(4)在步骤(3)的基础上得到了每一时刻的报警证据,在之后利用报警证据的动态更新得到全局报警证据,基本思想是利用每时刻的报警证据和前一时刻的全局报警证据进行加权融合,这样就可以得到全局报警证据,记全局报警证据为m0:t(m0:t(A),m0:t(NA),m0:t(θ)),具体步骤如下:
1)当t=1时,因在t=1之前没有相关的全局报警证据信息,所以t=1时全局报警证据是每一时刻归一化后的报警证据,即m0:1(m1(A),m1(NA),m1(θ))。
2)当t=2时全局报警证据利用t=1时刻的全局报警证据和t=2时刻的报警证据进行加权融合,得到t=2时刻的全局报警证据m0:2(m0:2(A),m0:2(NA),m0:2(θ)),其中m0:2(A)的计算如(7)式所示:
m0:2(A)=τ2m0:1(A)+ξ2m2(A|D) D∈(A,NA,θ) (7)
其中τ2=0.5表示对m0:1(A)的线性加权值,ξ2=0.5表示当前报警证据的加权值,同理可以得到m0:2(NA)和m0:2(θ)的全局报警证据如式(8)、(9)。
m0:2(NA)=τ2m0:1(NA)+ξ2m2(NA|D) D∈(A,NA,θ) (8)
m0:2(θ)=τ2m0:1(θ)+ξ2m2(θ|D) D∈(A,NA,θ) (9)
若m2(NA)>m2(A)和m2(NA)>m2(θ)则D=NA,那么m2(A|NA)=0,m2(NA|NA)=1,m2(θ|NA)=0;若m2(A)>m2(NA)和m2(A)>m2(θ)则D=A,那么m2(NA|A)=0,m2(NA|A)=1,m2(θ|A)=0;若m2(θ)>m2(NA)和m2(θ)>m2(A)则D=θ,那么m2(NA|θ)=m2(NA)/Pl(θ),m2(A|θ)=m2(A)/Pl(θ),m2(θ|θ)=m2(θ)/Pl(θ);其中Pl表示似真函数。
3)当t≥3时,m0:t=(m0:t(A),m0:t(NA),m0:t(θ)),其中m0:t(A)由(10)式确定
m0:t(A)=τtm0:t-1(A)+ξtmt(A|D) D∈(A,NA,θ) (10)
同理可以得到m0:t(NA)和m0:t(θ)的全局报警证据,其中τt表示上一时刻全局报警证据的权重因子,ξt表示当前时刻下的报警证据的权重因子。τt和ξt的计算如下面步骤所示:
a)当t≥3之后,可以得到每一时刻的报警证据mt和全局t-1以及t-2时刻的全局报警证据m0:t-1和m0:t-2,可以得到mt和m0:t-1以及m0:t-2两两全局报警证据之间的余弦值,其中其余弦值计算如下式(11)-(13)所示:
得到两两全局报警证据之间的余弦值,进而得到两两全局报警证据之间的证据距离d(mi,mj)=1-cos(mi,mj),其中i,j=t,t-1,t-2。
b)在上述步骤a)中,求得了两两全局报警证据距离,那么进一步得到两两全局报警证据之间的相似度Sim(mi,mj),所以mt、m0:t-1和m0:t-2之间的相似度如式(14)-(16)所示。
c)按照b)步骤可以得到两两全局报警证据之间的相似度,因此可以进一步的得到mt、m0:t-1和m0:t-2每个证据相对于其他两个证据的支持度,所以每个证据的支持度(17)-(19)式所示:
Sup(mt)=Sim(mt,m0:t-1)+Sim(mt,m0:t-2) (17)
Sup(mt-1)=Sim(mt,m0:t-1)+Sim(mt-1,m0:t-2) (18)
Sup(mt-2)=Sim(mt,m0:t-2)+Sim(mt-1,m0:t-2) (19)
d)在c)步骤的基础上可以得到每个证据相对于与其他两个证据的支持度,因此最后可以得到证据动态更新的权重因子τt和ξt,如下式(20)和(21)。
显然τt+ξt=1;
若mt(NA)>mt(A)和mt(NA)>mt(θ)则D=NA,那么mt(A|NA)=0,mt(NA|NA)=1,mt(θ|NA)=0;若mt(A)>mt(NA)和mt(A)>mt(θ)则D=A,那么mt(NA|A)=0,mt(NA|A)=1,mt(θ|A)=0;若mt(θ)>mt(NA)和mt(θ)>mt(A)则D=θ,那么mt(NA|θ)=mt(NA)/Pl(θ),mt(A|θ)=mt(A)/Pl(θ),mt(θ|θ)=mt(θ)/Pl(θ);其中Pl表示似真函数。
(5)根据步骤(4)得到t时刻的全局报警证据
m0:t=(m0:t(A),m0:t(NA),m0:t(θ)),将m0:t转化为相应的pignistic概率如下
(6)这步主要求得在误报漏报和漏报率最小时,找到最优阈值xotp,其基本步骤如下:
1)从在步骤(2)中所采集的过程变量x(t)中选取Num个过程变量,其中处于正常状态的过程变量有NumNA=0.5*Num,处于异常状态的过程变量有NumA=0.5*Num。
2)统计出1)中的过程变量,找到过程变量中最大值max(x(t))和最小值min(x(t)),其中要找的最后阈值xotp∈[min(x(t)),max(x(t))],搜寻最优阈值xotp按照一定的间隔,其中k=10,20,….,k的值越大搜寻的精度就会越高。
3)将1)中获得的采样序列x(t)′带入步骤(3)中得到报警证据
然后将报警证据进行归一化得到每一时刻的报警证据。
4)将上一步得到的每一时刻的报警证据带入步骤(4)中进行证据的动态更新,最终得到全局报警证据m0:t=(m0:t(A),m0:t(NA),m0:t(θ))。
5)按照2)中确定搜寻间隔后将xotp=min(x(t))+(n-1)Δ,根据4)中得到的全局报警证据进行报警决策,统计出被测试序列的误报个数M1和漏报个数M2,误报率和漏报率的计算如下:
FAR=M1/Num (26)
MAR=M2/Num (27)
逐个的按照一定的间隔进行搜寻xotp,得到在不同的xotp下的G=(FAR2+MAR2)0.5,找到使得G最小的xotp,就是所要搜寻的最优阈值。
(7)在上一步中得到最优阈值xotp后,即可按照步骤(3)-(5)获得报警结果。
本发明的有益效果:本发明将全局的不确定性纳入到证据体系当中,利用动态的证据更新融合,获得有效的报警证据,从而在报警决策上与工况状态下的描述更加接近,提高报警器的报警精度。同时本发明解决了不确定性因素对报警器的影响以及过程变量信息缺失的问题。根据本发明方法编写的程序可以在计算机上运行,并结合相关的传感器、数据采集器以及数据存储等硬件组成在线的报警***,实现对设备状况的在线监测。
附图说明
图1是本发明方法的总体流程框图。
图2是本发明方法的实施例中过程变量x的在线测量序列。
图3是本发明方法的实施例中过程变量x的测试样本序列。
图4是本发明方法的实施例中最优阈值在xotp=0.82时对应给出的模糊隶属度函数。
具体实施方式
本发明提出的一种基于全局不确定性的证据动态融合的工业报警器设计方法,其总体流程框图如图1所示,包括以下各步骤:
(1)设定某一工业报警器的辨识框架为θ={NA,A},其中NA表示工业报警器处于正常状态,A表示工业报警器处于异常状态亦即报警状态,θ表示处于全局不确定状态。
(2)设x为工业报警器所检测的过程变量,令x(t),t=1,2,3,….是传感器检测设备的在线序列,t为采样的时刻,采样的个数由工业报警器的监测周期和监控计算机设备来决定,在采集到的序列中,定义xotp为最优的阈值,记max(x(t))和min(x(t))分别是序列的最大值和最小值,其中xotp∈[min(x(t)),max(x(t))]。利用Sigmoid隶属度函数和Gaussian型隶属度函数分别构建正常状态NA、异常状态A和全局状态θ下的模糊隶属度函数μNA(x(t))、μA(x(t))和μθ(x(t)),如式(1)-(3)所示
其中a是一可调参数,当a取不同值时,隶属度函数会有不同的变化趋势,δ表示高斯型隶属度函数的方差,xotp表示最优的阈值。
需要说明的是:高斯型隶属度函数δ的大小决定全局不确定信息的宽窄程度,主要的模糊不确定信息集中在xotp的左右,在这里δ∈(0,0.5)。
(3)将t时刻的过程变量x值x(t)分别带入步骤(2)中关于正常状态NA、异常状态A和全局不确定状态下θ的模糊隶属度函数μNA(x)、μA(x)和μΘ(x)中,即可得到关于报警证据mt(NA)、mt(A)和mt(θ),分别如下(4)-(6)式所示。
将报警证据归一化,令Mt(x(t))=(mNA(x(t))+mA(x(t)+mθ(x(t)),则每一项归一化后的为,mt(NA)=(mNA(x(t))/Mt(x(t)))、mt(A)=(mA(x(t))/Mt(x(t)))、mt(θ)=(mθ(x(t))/Mt(x(t))),则称mt(NA)、mt(A)和mt(θ)是t时刻的报警证据。
需要说明的是,基于(4)-(6)式构造的模糊隶属度函数求解得到的mt(NA)、mt(A)和mt(θ)分别表示对“正常状态下”、“异常状态下”和“全局不确定状态下”的支持程度。
(3)将t时刻的过程变量x值x(t)分别带入步骤(2)中关于正常状态NA、异常状态A和全局状态下θ的模糊隶属度函数μNA(x)、μA(x)和mθ(x)中,即可得到关于报警证据mt(NA)、mt(A)和mt(θ),分别如下(4)-(6)式所示。
将报警证据归一化,令Mt(x(t))=(mNA(x(t))+mA(x(t)+mθ(x(t)),则每一项归一化后的为,mt(NA)=(mNA(x(t))/Mt(x(t)))、mt(A)=(mA(x(t))/Mt(x(t)))、mt(θ)=(mθ(x(t))/Mt(x(t))),则称mt(NA)、mt(A)和mt(θ)是t时刻的报警证据。
(4)在步骤(3)的基础上得到了每一时刻的报警证据,在之后利用报警证据的动态更新得到全局报警证据,基本思想是利用每时刻的报警证据和前一时刻的全局报警证据进行加权融合,这样就可以得到全局报警证据,记全局报警证据为m0:t(m0:t(A),m0:t(NA),m0:t(θ)),具体的的步骤如下:
1)当t=1时,因在t=1之前没有相关的全局报警证据信息,所以t=1时全局报警证据是每一时刻归一化后的报警证据,即m0:1(m1(A),m1(NA),m1(θ))。
2)当t=2时全局报警证据利用t=1时刻的全局报警证据和t=2时刻的报警证据进行加权融合,得到t=2时刻的全局报警证据m0:2(m0:2(A),m0:2(NA),m0:2(θ)),其中m0:2(A)的计算如(7)式所示:
m0:2(A)=τ2m0:1(A)+ξ2m2(A|D) D∈(A,NA,θ) (7)
其中τ2=0.5表示对m0:1(A)的线性加权值,ξ2=0.5表示当前报警证据的加权值,同理可以得到m0:2(NA)和m0:2(θ)的全局报警证据如式(8)、(9)。
m0:2(NA)=τ2m0:1(NA)+ξ2m2(NA|D) D∈(A,NA,θ) (8)
m0:2(θ)=τ2m0:1(θ)+ξ2m2(θ|D) D∈(A,NA,θ) (9)
若m2(NA)>m2(A)和m2(NA)>m2(θ)则D=NA,那么m2(A|NA)=0,m2(NA|NA)=1,m2(θ|NA)=0;若m2(A)>m2(NA)和m2(A)>m2(θ)则D=A,那么m2(NA|A)=0,m2(NA|A)=1,m2(θ|A)=0;若m2(θ)>m2(NA)和m2(θ)>m2(A)则D=θ,那么m2(NA|θ)=m2(NA)/Pl(θ),m2(A|θ)=m2(A)/Pl(θ),m2(θ|θ)=m2(θ)/Pl(θ);其中Pl表示似真函数。
3)当t≥3时,m0:t=(m0:t(A),m0:t(NA),m0:t(θ)),其中m0:t(A)由(10)式确定
m0:t(A)=τtm0:t-1(A)+ξtmt(A|D) D∈(A,NA,θ) (10)
同理可以得到m0:t(NA)和m0:t(θ)的全局报警证据,其中τt表示上一时刻全局报警证据的权重因子,ξt表示当前时刻下的报警证据的权重因子。τt和ξt的计算如下面步骤所示:
a)当t≥3之后,可以得到每一时刻的报警证据mt和全局t-1以及t-2时刻的全局报警证据m0:t-1和m0:t-2,可以得到mt和m0:t-1以及m0:t-2两两全局报警证据之间的余弦值,其中其余弦值计算如下式(11)-(13)所示:
因此得到两两全局报警证据之间的余弦值,进而得到两两全局报警证据之间的证据距离d(mi,mj)=1-cos(mi,mj),其中i,j=t,t-1,t-2。
b)在上述步骤a)中,求得了两两全局报警证据之间的证据距离,那么进一步得到两两全局报警证据之间的相似度Sim(mi,mj)=f(d(mi,mj)),所以mt、m0:t-1和m0:t-2之间的相似度如式(14)-(16)所示。
c)按照b)步骤可以得到两两全局报警证据之间的相似度,因此可以进一步的得到mt、m0:t-1和m0:t-2每个证据相对于其他两个证据的支持度,所以每个证据的支持度(17)-(19)式所示:
Sup(mt)=Sim(mt,m0:t-1)+Sim(mt,m0:t-2) (17)
Sup(mt-1)=Sim(mt,m0:t-1)+Sim(mt-1,m0:t-2) (18)
Sup(mt-2)=Sim(mt,m0:t-2)+Sim(mt-1,m0:t-2) (19)
d)在c)步骤的基础上可以得到每个证据相对于与其他两个证据的支持度,因此最后可以得到证据动态更新的权重因子τt和ξt,如下式(20)和(21)。
显然τt+ξt=1;若mt(NA)>mt(A)和mt(NA)>mt(θ)则D=NA,那么mt(A|NA)=0,mt(NA|NA)=1,mt(θ|NA)=0;若mt(A)>mt(NA)和mt(A)>mt(θ)则D=A,那么mt(NA|A)=0,mt(NA|A)=1,mt(θ|A)=0;若mt(θ)>mt(NA)和mt(θ)>mt(A)则D=θ,那么mt(NA|θ)=mt(NA)/Pl(θ),mt(A|θ)=mt(A)/Pl(θ),mt(θ|θ)=mt(θ)/Pl(θ);其中Pl表示似真函数。
为了加深对步骤(4)的理解,这里举个例子解释证据的动态更新的整个过程,假设参数a=8,t=1,2,3,4,5时刻的过程变量对应的每个时刻的报警证据如下表1所示:
表1过程变量每时刻报警证据
t | m<sub>t</sub>=(m<sub>t</sub>(A),m<sub>t</sub>(NA),m<sub>t</sub>(Θ)) | t | m<sub>t</sub>=(m<sub>t</sub>(A),m<sub>t</sub>(NA),m<sub>t</sub>(Θ)) |
1 | (0.1,0.7,0.2) | 4 | (0.32,0.58,0.1) |
2 | (0,2,0.6,0.2) | 5 | (0.7,0.16,0.14) |
3 | (0.56,0.23,0.21) | 6 |
根据步骤(4)计算5个时刻的全局报警证据
当t=1时,根据步骤1)知全局报警证据等于当前时刻的报警证据可得m0:1=m1=(0.1,0.7,0.2)。
当t=2时,根据步骤2)得到m2(NA)>m2(A),m2(NA)>m2(θ)则D=NA,m2(A|NA)=0,m2(NA|NA)=1,m2(θ|NA)=0,m0:2(A)=τ2m0:1(A)+ξ2m2(A|D)=0.05m0:2(NA)=τ2m0:1(NA)+ξ2m2(NA|D)=0.85,m0:2(θ)=τ2m0:1(θ)+ξ2m2(θ|D)=0.1,其中的ξ2=τ2=0.5,所以则有m0:2=(0.05,0.85,0.1)。
当t=3时,根据步骤3)可知,d(m0:1,m0:2)=0.6106,d(m0:2,m3)=0.811,d(m0:1,m3)=0.836,根据子步骤b)可以得到sim(m0:1,m0:2)=0.292,sim(m0:2,m3)=0.076,sim(m0:1,m3)=0.063,再根据子步骤c)可以得到Sup(m3)=0.139,Sup(m0:2)=0.368,Sup(m0:1)=0.355。再根据子步骤d)可以得到ξ3=0.162,τ3=0.838,又因为m3(NA)>m3(A),m3(NA)>m2(θ)则D=NA,m3(A|NA)=0,m3(NA|NA)=1,m3(θ|NA)=0,因此可以得到全局的报警m0:3=(0.203,0.713,0.084)
当t=4时,根据步骤3)可知,d(m0:2,m0:3)=0.611,d(m0:3,m4)=0.660,d(m0:2,m4)=0.656,根据子步骤b)可以得到sim(m0:2,m0:3)=0.292,sim(m0:3,m4)=0.217,sim(m0:2,m4)=0.223,再根据子步骤c)可以得到Sup(m4)=0.44,Sup(m0:3)=0.509,Sup(m0:2)=0.515。再根据子步骤d)可以得到ξ4=0.301,τ4=0.699,又因为若m4(A)>mt(NA)和m4(A)>mt(θ)则D=A,那么m4(NA|A)=0,m4(NA|A)=1,m4(θ|A)=0;因此可以得到全局的报警m0:4=(0.142,0.799,0.059)。
当t=5时,根据步骤3)可知,d(m0:3,m0:4)=0.613,d(m0:4,m5)=0.818,d(m0:3,m5)=0.847,根据子步骤b)可以得到sim(m0:3,m0:4)=0.288,sim(m0:4,m5)=0.073,sim(m0:3,m5)=0.058,再根据子步骤c)可以得到Sup(m5)=0.131,Sup(m0:4)=0.361,Sup(m0:3)=0.364。再根据子步骤d)可以得到ξ5=0.301,τ5=0.699,又因为m5(NA)>m5(A),m5(NA)>m5(θ)则D=NA,m5(A|NA)=0,m5(NA|NA)=1,m5(θ|NA)=0,因此可以得到全局的报警m0:5=(0.276,0.674,0.05)。
(5)根据步骤(4)可以得到全局的报警证据m0:t=(m0:t(A),m0:t(NA),m0:t(θ)),根据pignistic概率进行决策,如果若则工业报警器发出报警,则说明处于异常状态;反之,不发出报警,说明处于正常状态。
在(4)中的例子中给出5个全局的报警证据可以对其进行决策,可以给出报警结果如下表。
表2报警结果输出
实验分析:通过表中的数据可以看出,在正常状态下,利用每时刻的报警证据进行决策会有一定的误差,而基于本发明方法融合后的全局报警证据更加真实的反映了现实的工况(正常状态)。
(6)这步主要求得在误报漏报和漏报率最小时,找到最优阈值xotp,其基本步骤如下:
1)从在步骤(2)中所采集的过程变量x(t)中选取Num个过程变量,其中处于正常状态的过程变量有NumNA=0.5*Num,处于异常状态的过程变量有NumA=0.5*Num。
2)统计出1)中的过程变量,找到过程变量中最大值max(x(t))和最小值min(x(t)),其中要找的最后阈值xotp∈[min(x(t)),max(x(t))],搜寻最优阈值xotp按照一定的间隔,其中k=10,20,….,k的值越大搜寻的精度就会越高。
3)将1)中获得的采样序列x(t)′带入步骤(3)中得到报警证据
然后将报警证据进行归一化得到每一时刻的报警证据。
4)将上一步得到的没每一时刻的报警证据带入步骤(4)中进行证据的动态更新,最终得到全局报警证据m0:t=(m0:t(A),m0:t(NA),m0:t(θ))。
5)按照2)中确定搜寻间隔后将xotp=min(x(t))+Δ,根据4)中得到的全局报警证据进行报警决策,统计出被测试序列的误报个数M1和漏报个数M2,误报率和漏报率的计算如下:
FAR=M1/Num (26)
MAR=M2/Num (27)
逐个的按照一定的间隔进行搜寻xotp,得到在不同的xotp下的G=(FAR2+MAR2)0.5,找到使得G最小的xotp,就是所要搜寻的最优阈值。
需要说明,在此步骤中选取的测试样本个数Num的个数越大越好,在求取xotp时的间隔k也应尽可能的大,这样做得到的最优阈值更加精确,使得得到的报警结果误报率和漏报率更加接近真实状况。
(7)在上一步中得到最优阈值xotp后,即可按照步骤(3)-(5)获得报警结果。
以下结合附图,详细的介绍本发明方法的实施例:
本发明方法的主要流程图如图1所示,其主要内容为:通过实时监测采集相关的过程变量信息,将过程变量带入到S型和高斯型的模糊隶属度函数中去,得到每时刻的报警证据,再根据动态的证据更新得到全局的报警证据,最后利用pignistic概率进行报警决策。在将过程变量转换成每时刻的报警证据时,将全局的不确定性纳入到证据中去,有效的解决了全局的不确定性,利用连续型的模糊隶属函数转换,这样使得获得的信息更加完整,这样能够提升报警器的报警性能。在过程变量的测试样本中给出了在最低误报率和漏报率下获得最优阈值的方法。
以下结合附图2中在线监测得到的过程变量数据,给出最佳实施例,详细的介绍本发明的相关详细步骤,通过实验数据在最优阈值的条件下,得到的误报率与漏报率要比在传统的融合方法(分段式的模糊隶属度函数以及未考虑全局的不确定性)上更加低。
1、在线监测获得过程变量的信息x(t)
过程变量采集到的序列如图2所示,时间t的最大值为4000,通过统计可以得到该过程变量x(t)的最大值是max(x(t))=7.123,min(x(t))=-5.329。
2、求取在最小误报率和漏报率下的最优阈值xotp
根据步骤(6)中的子步骤1),经过分析,图2所示的在线采样序列x(t)中,x(2000)到x(2500)和x(3000)到x(3500)这1000个采样数据处于正常的工作状态,x(2500)到x(3000)和x(3500)到x(4000)这1000个采样数据处于异常的工作状态。将这两组共2000个采样数据排列成一段序列x′(t),如附图3所示,t=1,2,3….,num=2000,这2000个数据x′(t)∈[-5.329,7.123],其中前1000个是设备处于正常运行状态下的数据M1,后1000个是设备处于异常状态下的数据M2,num=M1+M2。
根据步骤(6)的子步骤3)-5)可得,当xotp=0.82,a=8,δ=0.05时,在最优阈值下(如图4所示)其对应的G=(FAR2+MAR2)0.5=0.0264,其中FAR=2.3%,MAR=1.3%。
3、在获得过程变量x(t)的最优阈值后,即可按照步骤(2)-(5)对监测的在线序列x(t),t=1,2,3…得到报警结果。在传统的融合报警器的设计方法中,在获得每时刻的报警证据是利用分段式的模糊隶属函数,而且没有考虑到全局的不确定性,如果按照所参考的发明专利,在区间[min(x(t))=-5.329,max(x(t))=7.123]中搜索使得G=(FAR2+MAR2)0.5最小,此时xotp=0.76,而在此时对应的误报率与漏报率是FAR=6.9%,MAR=5.5%。从最后的报警结果可以看出,本发明在最优阈值的下产生的误报率和漏报率都比传统的方法更加精确。
Claims (1)
1.一种基于全局不确定性证据动态融合的工业报警器设计方法,其特征在于该方法包括以下各步骤:
步骤(1)设定某一工业报警器的辨识框架为θ={NA,A},其中NA表示工业报警器处于正常状态,A表示工业报警器处于异常状态亦即报警状态,θ表示处于全局不确定状态;
步骤(2)设x为工业报警器所检测的过程变量,令x(t),t=1,2,3,….是传感器检测设备的在线序列,t为采样的时刻,采样的个数由工业报警器的监测周期和监控计算机设备来决定,在采集到的序列中,定义xotp为最优的阈值,记max(x(t))和min(x(t))分别是序列的最大值和最小值,其中xotp∈[min(x(t)),max(x(t))];利用Sigmoid隶属度函数和Gaussian型隶属度函数分别构建正常状态NA、异常状态A和全局不确定状态θ下的模糊隶属度函数μNA(x(t))、μA(x(t))和μθ(x(t)),如式(1)-(3)所示
其中a是一可调参数,当a取不同值时,隶属度函数会有不同的变化趋势,δ表示高斯型隶属度函数的方差;
步骤(3)将t时刻的过程变量x值x(t)分别带入步骤(2)中关于正常状态NA、异常状态A和全局不确定状态下θ的模糊隶属度函数μNA(x(t))、μA(x(t))和μθ(x(t))中,即可得到关于报警证据mt(NA)、mt(A)和mt(θ),分别如下(4)-(6)式所示;
将报警证据归一化,令Mt(x(t))=(mNA(x(t))+mA(x(t)+mθ(x(t)),则每一项归一化后的为,mt(NA)=(mNA(x(t))/Mt(x(t)))、mt(A)=(mA(x(t))/Mt(x(t)))、mt(θ)=(mθ(x(t))/Mt(x(t))),则称mt(NA)、mt(A)和mt(θ)是t时刻的报警证据;
步骤(4)在步骤(3)的基础上得到了每一时刻的报警证据,在之后利用报警证据的动态更新得到全局报警证据,基本思想是利用每一时刻的报警证据和前一时刻的全局报警证据进行加权融合,这样就可以得到全局报警证据,记全局报警证据为m0:t(m0:t(A),m0:t(NA),m0:t(θ)),具体步骤如下:
1)当t=1时,因在t=1之前没有相关的全局报警证据信息,所以t=1时全局报警证据是每一时刻归一化后的报警证据,即m0:1(m1(A),m1(NA),m1(θ));
2)当t=2时全局报警证据利用t=1时刻的全局报警证据和t=2时刻的报警证据进行加权融合,得到t=2时刻的全局报警证据m0:2(m0:2(A),m0:2(NA),m0:2(θ)),其中m0:2(A)的计算如(7)式所示:
m0:2(A)=τ2m0:1(A)+ξ2m2(A|D) D∈(A,NA,θ) (7)
其中τ2=0.5表示对m0:1(A)的线性加权值,ξ2=0.5表示当前报警证据的加权值,同理可以得到m0:2(NA)和m0:2(θ)的全局报警证据如式(8)、(9);
m0:2(NA)=τ2m0:1(NA)+ξ2m2(NA|D) D∈(A,NA,θ) (8)
m0:2(θ)=τ2m0:1(θ)+ξ2m2(θ|D) D∈(A,NA,θ) (9)
若m2(NA)>m2(A)和m2(NA)>m2(θ)则D=NA,那么m2(A|NA)=0,m2(NA|NA)=1,m2(θ|NA)=0;若m2(A)>m2(NA)和m2(A)>m2(θ)则D=A,那么m2(NA|A)=0,m2(NA|A)=1,m2(θ|A)=0;若m2(θ)>m2(NA)和m2(θ)>m2(A)则D=θ,那么m2(NA|θ)=m2(NA)/Pl(θ),m2(A|θ)=m2(A)/Pl(θ),m2(θ|θ)=m2(θ)/Pl(θ);其中Pl表示似真函数;
3)当t≥3时,m0:t=(m0:t(A),m0:t(NA),m0:t(θ)),其中m0:t(A)由(10)式确定
m0:t(A)=τtm0:t-1(A)+ξtmt(A|D) D∈(A,NA,θ) (10)
同理可以得到m0:t(NA)和m0:t(θ)的全局报警证据,其中τt表示上一时刻全局报警证据的权重因子,ξt表示当前时刻下的报警证据的权重因子;τt和ξt的计算如下面步骤所示:
a)当t≥3之后,可以得到每一时刻的报警证据mt和全局t-1以及t-2时刻的全局报警证据m0:t-1和m0:t-2,得到mt和m0:t-1以及m0:t-2两两证据之间的余弦值,其中其余弦值计算如下式(11)-(13)所示:
因此得到两两证据之间的余弦值后,进而得到两两证据之间的证据距离d(mi,mj)=1-cos(mi,mj),其中i,j=t,t-1,t-2;
b)在上述步骤a)中,求得两两证据之间的证据距离后,进一步得到两两证据之间的相似度Sim(mi,mj)=f(d(mi,mj)),则mt、m0:t-1和m0:t-2之间的相似度如式(14)-(16)所示;
c)按照b)步骤得到两两证据之间的相似度,进一步得到mt、m0:t-1和m0:t-2每个证据相对于其他两个证据的支持度,所以每个证据的支持度(17)-(19)式所示:
Sup(mt)=Sim(mt,m0:t-1)+Sim(mt,m0:t-2) (17)
Sup(mt-1)=Sim(mt,m0:t-1)+Sim(mt-1,m0:t-2) (18)
Sup(mt-2)=Sim(mt,m0:t-2)+Sim(mt-1,m0:t-2) (19)
d)在c)步骤的基础上得到每个证据相对于与其他两个证据的支持度,最后得到证据动态更新的权重因子τt和ξt,如下式(20)和(21);
显然τt+ξt=1;
若mt(NA)>mt(A)和mt(NA)>mt(θ)则D=NA,那么mt(A|NA)=0,mt(NA|NA)=1,mt(θ|NA)=0;若mt(A)>mt(NA)和mt(A)>mt(θ)则D=A,那么mt(NA|A)=0,mt(NA|A)=1,mt(θ|A)=0;若mt(θ)>mt(NA)和mt(θ)>mt(A)则D=θ,那么mt(NA|θ)=mt(NA)/Pl(θ),mt(A|θ)=mt(A)/Pl(θ),mt(θ|θ)=mt(θ)/Pl(θ);其中Pl表示似真函数;
步骤(5)根据步骤(4)得到t时刻的全局报警证据m0:t=(m0:t(A),m0:t(NA),m0:t(θ)),将m0:t转化为相应的pignistic概率如下
步骤(6)主要求得在误报漏报和漏报率最小时,找到最优阈值xotp,其基本步骤如下:
1)从在步骤(2)中所采集的过程变量x(t)中选取Num个过程变量,形成采样序列x(t)′,其中处于正常状态的过程变量有NumNA=0.5*Num,处于异常状态的过程变量有NumA=0.5*Num;
2)统计出1)中选取的过程变量,找到过程变量中最大值max(x(t))和最小值min(x(t)),其中要找的最优阈值xotp∈[min(x(t)),max(x(t))],搜寻最优阈值xotp按照一定的间隔,其中k=10,20,….;k的值越大搜寻的精度就会越高;
3)将1)中获得的采样序列x(t)′带入步骤(3)中得到报警证据
然后将报警证据进行归一化得到每一时刻的报警证据;
4)将上一步得到的每一时刻的报警证据带入步骤(4)中进行证据的动态更新,最终得到全局报警证据m0:t=(m0:t(A),m0:t(NA),m0:t(θ));
5)按照2)中确定搜寻间隔后将xotp=min(x(t))+(n-1)△,根据4)中得到的全局报警证据进行报警决策,统计出被测试序列的误报个数M1和漏报个数M2,误报率和漏报率的计算如下:
FAR=M1/Num (26)
MAR=M2/Num (27)
逐个的按照一定的间隔进行搜寻xotp,得到在不同的xotp下的G=(FAR2+MAR2)0.5,找到使得G最小的xotp,就是所要搜寻的最优阈值;
步骤(7)得到最优阈值xotp后,按照步骤(3)-(5)获得报警结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810082465.4A CN108257365B (zh) | 2018-01-29 | 2018-01-29 | 一种基于全局不确定性证据动态融合的工业报警器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810082465.4A CN108257365B (zh) | 2018-01-29 | 2018-01-29 | 一种基于全局不确定性证据动态融合的工业报警器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108257365A CN108257365A (zh) | 2018-07-06 |
CN108257365B true CN108257365B (zh) | 2020-04-24 |
Family
ID=62742490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810082465.4A Active CN108257365B (zh) | 2018-01-29 | 2018-01-29 | 一种基于全局不确定性证据动态融合的工业报警器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108257365B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109145972B (zh) * | 2018-08-09 | 2021-08-03 | 杭州电子科技大学 | 一种船舶电力推进***变频器报警器设计方法 |
CN109765786B (zh) * | 2019-01-25 | 2022-03-01 | 杭州电子科技大学 | 一种基于证据滤波的电推船电机转轴不平衡故障检测方法 |
CN109932908B (zh) * | 2019-03-20 | 2022-03-01 | 杭州电子科技大学 | 一种基于报警信度融合的多向主元分析过程监测方法 |
CN109946542B (zh) * | 2019-03-22 | 2021-01-01 | 杭州电子科技大学 | 基于长短时信息融合的逆变器直流回路电容故障检测方法 |
CN111443686B (zh) * | 2020-03-23 | 2021-08-10 | 杭州电子科技大学 | 基于多目标优化与证据迭代更新的工业报警器设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6415276B1 (en) * | 1998-08-14 | 2002-07-02 | University Of New Mexico | Bayesian belief networks for industrial processes |
CN102651028A (zh) * | 2012-04-09 | 2012-08-29 | 云南大学 | 基于d-s证据理论的不确定性数据世系查询处理方法 |
CN102789676A (zh) * | 2012-08-10 | 2012-11-21 | 杭州电子科技大学 | 一种基于报警证据融合的工业报警器设计方法 |
CN107065834A (zh) * | 2017-05-25 | 2017-08-18 | 东北大学 | 湿法冶金过程中浓密机的故障诊断方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130041625A1 (en) * | 2011-08-11 | 2013-02-14 | International Business Machines Corporation | Advanced Statistical Detection of Emerging Trends |
-
2018
- 2018-01-29 CN CN201810082465.4A patent/CN108257365B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6415276B1 (en) * | 1998-08-14 | 2002-07-02 | University Of New Mexico | Bayesian belief networks for industrial processes |
CN102651028A (zh) * | 2012-04-09 | 2012-08-29 | 云南大学 | 基于d-s证据理论的不确定性数据世系查询处理方法 |
CN102789676A (zh) * | 2012-08-10 | 2012-11-21 | 杭州电子科技大学 | 一种基于报警证据融合的工业报警器设计方法 |
CN107065834A (zh) * | 2017-05-25 | 2017-08-18 | 东北大学 | 湿法冶金过程中浓密机的故障诊断方法 |
Non-Patent Citations (2)
Title |
---|
一种基于证据理论的工业报警器设计方法;宋晓静等;《杭州电子科技大学学报》;20111231;第31卷(第6期);全文 * |
基于置信规则库推理的证据滤波报警器设计;徐海洋等;《山东科技大学学报》;20170831;第36卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108257365A (zh) | 2018-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108257365B (zh) | 一种基于全局不确定性证据动态融合的工业报警器设计方法 | |
US11586913B2 (en) | Power equipment fault detecting and positioning method of artificial intelligence inference fusion | |
Sun et al. | Fault detection and identification using Bayesian recurrent neural networks | |
CN111665066B (zh) | 基于卷积神经网络的设备故障自适应上下预警界生成方法 | |
BR112021013786A2 (pt) | Método e sistema para detectar anomalias ou indicações prematuras de falha de equipamento | |
CN112799898B (zh) | 基于分布式故障检测的互联***故障节点定位方法及*** | |
CN115169479A (zh) | 污水处理过程远程监控方法、***及存储介质 | |
EP1960853A1 (en) | Evaluating anomaly for one-class classifiers in machine condition monitoring | |
US12007745B2 (en) | Apparatus for predicting equipment damage | |
CN117270514B (zh) | 基于工业物联网的生产过程全流程故障检测方法 | |
CN111964909A (zh) | 滚动轴承运行状态检测方法、故障诊断方法及*** | |
CN112734977B (zh) | 一种基于物联网的设备风险预警***及算法 | |
CN113721182B (zh) | 一种电力互感器性能在线监测结果可靠度评估方法及*** | |
Nadir et al. | Utilizing principal component analysis for the identification of gas turbine defects | |
CN116907772A (zh) | 桥梁结构监测传感器的自诊断与故障源鉴别方法及*** | |
CN110057588B (zh) | 基于奇异值与图论特征融合的轴承早期故障检测与诊断方法及*** | |
CN116910677A (zh) | 一种工业仪表故障诊断方法及*** | |
US11339763B2 (en) | Method for windmill farm monitoring | |
CN116360384A (zh) | 用于诊断和监控信息物理***中的异常的***和方法 | |
CN113177484B (zh) | 基于lof自编码的机械故障智能诊断方法 | |
Dash et al. | A Comparison of Model-Based and Machine Learning Techniques for Fault Diagnosis | |
Wang et al. | Fault detection based on optimal transport theory | |
RU2784981C1 (ru) | Способ диагностики и мониторинга аномалий в кибер-физической системе | |
CN114637793B (zh) | 一种基于大数据分析的装备故障频发区域定位方法 | |
CN117723782B (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 |