CN104238516B - 一种锅炉系统设备状态监测方法 - Google Patents

一种锅炉系统设备状态监测方法 Download PDF

Info

Publication number
CN104238516B
CN104238516B CN201410468800.6A CN201410468800A CN104238516B CN 104238516 B CN104238516 B CN 104238516B CN 201410468800 A CN201410468800 A CN 201410468800A CN 104238516 B CN104238516 B CN 104238516B
Authority
CN
China
Prior art keywords
mtd
mrow
msubsup
state
msup
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
CN201410468800.6A
Other languages
English (en)
Other versions
CN104238516A (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.)
Peisen Fujian Environmental Technology Co ltd
Original Assignee
Xiamen University
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 Xiamen University filed Critical Xiamen University
Priority to CN201410468800.6A priority Critical patent/CN104238516B/zh
Publication of CN104238516A publication Critical patent/CN104238516A/zh
Application granted granted Critical
Publication of CN104238516B publication Critical patent/CN104238516B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明涉及锅炉技术领域,具体涉及一种锅炉系统设备状态监测方法。本发明提供了基于不同时刻系统状态的变化,监控设备状态变化概率的方法,通过对锅炉设备中的系统状态进行监控,根据系统状态判断锅炉中各重要设备的工作状态,从而有效地定位故障位置,并通过对所有系统状况的变化进行分析,预测出各种设备状态出现的概率。通过这种方法,可以帮助决策者更加及时地定位故障并对可能发生的故障做好充足的准备。

Description

一种锅炉系统设备状态监测方法
技术领域
本发明涉及锅炉技术领域,具体涉及一种锅炉系统设备状态监测方法。
背景技术
火力发电是我国主要的发电方式,电站锅炉作为火力电站的三大主机设备之一,伴随着我国火电行业的发展而发展。《2013-2017年中国锅炉制造行业产销需求预测与转型升级分析报告》显示,至2010年底,单机容量30万千瓦及以上火电机组占全部火电机组容量的60%以上。
锅炉设备是一种受压设备,经常处于高压下运行,并受到烟气和炉水中有害杂质侵蚀和飞灰磨损,因此需要定期的对锅炉系统进行维护并在运行过程中实时监测可能发生的设备故障。准确的查找锅炉设备故障正是制约锅炉正常使用和发展的重要因素之一,一旦管理不严,使用不当就会发生事故,甚至发生重大破坏性事故。在出现事故后及时的查找事故根源的速度越快越准确,事故造成的损失就会越小。但故障后及时地查找故障根源保证其正常使用是困扰锅炉设备发展的一大难题,一般情况下都是凭借人工经验去推断,往往准确率极低。
锅炉系统状态由该系统的多个变量决定。这些变量的不同取值构成了系统的不同状态。另外,系统中的每个设备元件也有多种不同的工作状态。设备状态与系统状态之间相互影响。但目前的设备控制程序只是单一的将系统元件与系统状态建立一一映射的关系,但在实际情况中,许多系统中的重要设备与系统状态的映射关系并不是一一对应的,系统中一个设备的状态发生变化可能会引起多个系统状态同时发生变化。例如在发电厂的锅炉系统中,增加通风阀门的开度,不仅可以使锅炉中含氧量提高,同时会对锅炉中的压力产生影响。而且,在即便在与系统状态存在一一映射的系统元件中,有的设备会对系统状态产生抑制总用,有的会对系统状态产生增益作用。因此,设备状态的转移概率在监控过程中,需要充分考虑设备的种类和作用对象。
发明内容
解决上述技术问题,本发明提供了基于不同时刻系统状态的变化,监控设备状态变化概率的方法,通过对锅炉设备中的系统状态进行监控,根据系统状态判断锅炉中各重要设备的工作状态,从而有效地定位故障位置,并通过对所有系统状况的变化进行分析,预测出各种设备状态出现的概率。通过这种方法,可以帮助决策者更加及时地定位故障并对可能发生的故障做好充足的准备。
为了达到上述目的,本发明所采用的技术方案是,一种锅炉系统设备状态监测方法,方法包括以下步骤:
步骤一:将每套锅炉系统由n个系统状态变量表示,并根据这n个系统状态变量的不同取值组合来表示每套锅炉系统不同的状态,每一锅炉系统包括m个设备,Ni表示第i个设备的工作状态类型个数,而则表示第t时刻的第i个设备的工作状态,其中,i=1,2,3.......m,确定锅炉系统状态变量集合为H,其中第i个系统状态变量用Hi表示;
步骤二:将第k个锅炉系统中的系统状态进行细分,设细分后锅炉系统状态组合总数为K,第k个锅炉系统状态组合区间为Mk
步骤三:设定设备状态间的转移概率λ;
步骤四:设定设备工作状态集合Y与系统状态变量集合H的对应关系W;
步骤五:计算第i个设备在第t'时刻的设备状态到第t时刻的设备状态的概率其中t表示当前时刻,t'表示前一时刻,当前时刻t与前一时刻t'的时间间隔为单位时间△t,第i个设备在单位时间△t内的状态变化的概率相同,第i个设备在t'、t时刻的设备状态为
进一步的,所述步骤三具体为,设定第i个设备状态间的转移概率λ即第i个设备状态转移矩阵λi(i=1,2,3...,m),其计算公式如下:
其中,表示第i个设备从状态1变化到状态1的概率,表示第i个设备从状态1变化到状态2的概率,以此类推,表示第i个设备从态yi变化到状态yi的概率。
进一步的,所述步骤四具体为:该对应关系W为第i个设备在开启状态下对锅炉系统的n个系统状态变量的影响趋势,设定设备工作状态集合Y与系统状态变量集合H的对应关系W;
其公式为:
其中wij表示第i个设备对第j个系统状态变量的影响值,所述j=1,2,3......n,wij的取值范围为{0,1,-1},分别表示第i个设备设备对第j个系统状态变量没有影响、有增益影响、有抑制影响。
进一步的,所述步骤具体包括以下步骤:
步骤51:判断第i个设备在前一时刻t'是否处于正常工作状态,若是,则执行步骤53,否则执行步骤52;判断方法:设定第i个设备工作状态集合维数为Li,而设定第i个设备工作状态集合为Yi={y1,y2,y3,...,yn},其中n为自然数,则表示该第i个设备具有n个不同的工作状态,因此Li=n,若第i个设备在前一时刻t'的工作状态不属于集合Yi={y1,y2,y3,...,yn}中的正常工作状态,则判定该第i个设备不处于正常工作状态,异则反之。
步骤52:分别计算第i个设备在当前时刻t的状态的转移概率Pi,具体公式如下:
a)若当前时刻t第i个设备非正常工作:
b)若当前时刻t第i个设备正常工作:
步骤53:判断第i个设备在当前时刻t是否处于正常工作状态,若是,则计算第i个设备在当前时刻t的工作状态的转移概率Pi的公式如下:
其中当前时刻t的锅炉系统状态为Mt,前一时刻t'的锅炉系统状态为Mt',R(Mt|Mt')是指设备在系统状态Mt'→Mt的允许变化状态集合,系统状态Mt'→Mt允许变化的状态集合与W的关系为:R(Mt|Mt')∈W,表示的是第t时刻第i个设备的工作状态,这里的允许变化状态集合R(Mt|Mt')是按照且R(Mk|Mk')∈W的方式进行的。
否则,判定第i个设备在当前时刻t的工作状态是否满足设备状态变化公式:
若是,则执行步骤54,否则,执行步骤55。
步骤54:计算第i个设备在当前时刻t的工作状态的转移概率Pi的公式如下:
步骤55:计算第i个设备在当前时刻t的状态的转移概率Pi的公式如下:
本发明通过采用上述技术方案,与现有技术相比,具有如下优点:
本发明通过遍历系统中所有可能的状态变化方式,结合系统种设备对于系统中状态变量的影响方式计算设备状态转移概率,计算结果可以直接为使用者提供设备处于某种状态时系统中的状态变量最可能发生的变化,同时提供给使用者发生这种变化的概率,从而帮助使用者做好相应的准备。同时,本发明在对设备的类型进行分类,使其可以适用于系统中不同的设备控制方式,可以更有效的进行设备转移概率的计算分析。
附图说明
图1是本发明的实施例的流程图;
图2锅炉系统每种工作状态范围。
具体实施方式
现结合附图和具体实施方式对本发明进一步说明。
一种锅炉系统设备状态监测方法,参考图1所示,包括以下步骤:
步骤一:将每套锅炉系统由n个系统状态变量表示,并根据这n个系统状态变量的不同取值组合来表示每套锅炉系统不同的状态,本文中K=k=1,因此锅炉系统包括m个设备,Ni表示第i个设备的工作状态类型个数,而则表示第t时刻的第i个设备的工作状态,其中,i=1,2,3.......m,确定锅炉系统状态变量集合为H,其中第i个系统状态变量用Hi表示。作为一个具体的实施例,本发明采用为了简化计算过程,在锅炉系统中取锅炉中的蒸汽温度和蒸汽压力两个较为关键的系统状态变量,进行分析。因此变量个数n为2,n=2。则锅炉系统状态变量集合为H=[H1,H2]=[蒸汽温度,蒸汽压力]。
目前大型的电厂,多数采用亚临界、超临界机组。本例中取亚临界机组作为研究对象。亚临界中各主要系统变量的工作范围一般为:主蒸汽压力16-18MPa,蒸汽温度535-540℃。因此设定锅炉蒸汽压力正常运行范围为[16,18],蒸汽温度正常运行范围为[535,540]。
步骤二:将第k个锅炉系统中的系统状态进行细分,设细分后锅炉系统状态组合总数为K,第k个锅炉系统状态组合区间为Mk。本发明中的k=1,设定细分后锅炉系统状态组合总数K=30种工作状态,每种工作状态范围如图2所示,假定系统状态为蒸汽温度和蒸汽压力两种,将蒸汽温度的正常工作范围平均分成4份,将蒸汽压力的正常工作范围平均分成3份,因此锅炉系统整体运行状况范围(包括正常运行情况、低于正常范围、高于正常范围)可分成(4+1+1)*(3+1+1)=30个细分区间,即图2中的表格示意出的即为本发明的锅炉系统状态组合区间Mk=M1
定义锅炉系统中分别由燃料输送装置和通气装置两个设备控制这两个变量,因此该锅炉系统包括m=2个设备,Ni表示第i个设备的工作状态类型个数,其中,i=1,2,每个设备均有正常运行、故障-开、故障-关三种状态,即Ni=3,分别用1、2、3表示。而则表示第t时刻的第i个设备的工作状态,其中,i=1,2。
步骤三:设定设备状态间的转移概率λ。设定第i个设备状态间的转移概率λ即第i个设备状态转移矩阵λi(i=1,2,3...,m),其计算公式如下:
其中,表示第i个设备从状态1变化到状态1的概率,表示第i个设备从状态1变化到状态2的概率,以此类推,表示第i个设备从态yi变化到状态yi的概率。本实施例中,根据实际运行的历史数据,设定每一设备的每种状态变化的概率,设备在各种运行状态下相互转换的概率(例如,设备从正常工作状态变为故障-关工作状态下的概率,本发明中设燃料输送装置和通气装置中的工作状态相互转移概率为λ,
λ=[λ12]
其中
其中假设两种控制设备在固定时间步长△t内均为不可修复系统。
步骤四:设定设备工作状态集合Y与系统状态变量集合H的对应关系W。该对应关系W为第i个设备在开启状态下对锅炉系统的n个系统状态变量的影响趋势,设定设备工作状态集合Y与系统状态变量集合H的对应关系W;其公式为:
其中wij表示第i个设备对第j个系统状态变量的影响值,所述j=1,2,3......n,wij的取值范围为0,1,-1,分别表示第i个设备设备对第j个系统状态变量没有影响、有增益影响、有抑制影响。本实施例中,设定燃料输送装置会应影响锅炉中的温度,并且输送装置的正常运行会使锅炉中的蒸汽温度升高;而通气装置的正常运行会时锅炉中的蒸汽压力降低,因此有,
步骤五:计算第i个设备在第t'时刻的设备状态到第t时刻的设备状态的概率P,其中t表示当前时刻,t'表示前一时刻,当前时刻t与前一时刻t'的时间间隔为单位时间△t,并假设第i个设备在单位时间△t内的状态变化的概率相同。
本实施例以给定的单位时间△t=3min为观测系统状态变化的时间间隔,已知t'、t时刻的设备工作状态将设备工作状态分别遍历所有锅炉系统状态组合区间Mk,详细见图2中的蒸汽温度-蒸汽压力表,其中将系统状态细分为30个子区间。根据判断系统所处的状态(Mt,Mt'),进一步判断系统状态(Mt,Mt')是否符合燃料输送装置和通气装置的设定的运行规则W,根据设备状态转移概率计算原理计算每个设备在系统状态处于Mt'→Mt的过程的设备的变化概率部分系统状态转移概率如下表1所示,
表1
具体包括以下步骤:
步骤51:判断第i个设备在前一时刻t'是否处于正常工作状态,若是,则执行步骤53,否则执行步骤52;判断方法:设定第i个设备工作状态集合维数为Li,而设定第i个设备工作状态集合为Yi={y1,y2,y3,...,yn},其中n为自然数,则表示该第i个设备具有n个不同的工作状态,因此Li=n,若第i个设备在前一时刻t'的工作状态不属于集合Yi={y1,y2,y3,...,yn}中的正常工作状态,则判定该第i个设备不处于正常工作状态,异则反之。
本实施例中,参考图2的细分区间中,7、8、9、12、13、14、17、18、19、22、23、24为系统正常工作区间,其余均为系统不正常工作区间。
步骤52:分别计算第i个设备在当前时刻t的状态的转移概率Pi,具体公式如下:
c)若当前时刻t第i个设备非正常工作:
d)若当前时刻t第i个设备正常工作:
例如第2个设备在t'时刻的处于不正常工作状态(设为故障-关:设定设备状态为2),当第2个设备在t时刻处于正常工作状态(设定设备状态为1),则
当第2个设备在t时刻处于不正常工作状态(即故障-关:设定设备状态为2)时,有
因为根据λ2可知,设备2为不可修复系统,因此的概率为0,的概率为1。
步骤53:判断第i个设备在当前时刻t是否处于正常工作状态,若是,则计算第i个设备在当前时刻t的工作状态的转移概率Pi的公式如下:
其中当前时刻t的锅炉系统状态为Mt,前一时刻t'的锅炉系统状态为Mt',R(Mt|Mt')是指设备在系统状态Mt'→Mt的允许变化状态集合,系统状态Mt'→Mt允许变化的状态集合与W的关系为:R(Mt|Mt')∈W,表示的是第t时刻第i个设备的工作状态,这里的允许变化状态集合R(Mt|Mt')是按照且R(Mk|Mk')∈W的方式进行的。
否则,判定第i个设备在当前时刻t的工作状态是否满足如下公式:
若是,则执行步骤54,否则,执行步骤55。
本实施例中,前一时刻t'时第2个设备处于正常状态(设定设备状态为1),若在时刻t时第2个设备也处于正常工作状态(设定设备状态为1),则
满足执行步骤55,
步骤54:计算第i个设备在当前时刻t的工作状态的转移概率Pi的公式如下:
本实施例中,在时刻t时第2个设备处于非正常工作状态(设为故障-关:设定设备状态为2),根据W的设定可知设备2对系统状态变量中的蒸汽压力是抑制作用,因此,当(Mt|Mt')=(12|14)时,
步骤55:计算第i个设备在当前时刻t的状态的转移概率Pi的公式如下:
当(Mt|Mt')=(14|12)时,因此
通过遍历系统中所有可能的状态变化方式,结合系统种设备对于系统中状态变量的影响方式计算设备状态转移概率,不仅可以告诉决策者直观的知道设备处于某种状态时系统中的状态变量会发生怎么样的变化,同时告诉决策者这种变化的概率是多少,从而帮助决策者做好相应的准备。同时,在对设备的类型进行分类,使其可以适用于系统中不同的设备控制方式,可以更有效的进行设备转移概率的计算分析。
尽管结合优选实施方案具体展示和介绍了本发明,但所属领域的技术人员应该明白,在不脱离所附权利要求书所限定的本发明的精神和范围内,在形式上和细节上可以对本发明做出各种变化,均为本发明的保护范围。

Claims (4)

1.一种锅炉系统设备状态监测方法,其特征在于:包括以下步骤:
步骤一:将每一个锅炉系统由n个系统状态变量表示,并根据这n个系统状态变量的不同取值组合来表示每一个锅炉系统不同的状态,每一个锅炉系统包括m个设备,Ni表示第i个设备的工作状态类型个数,而则表示第t时刻的第i个设备的工作状态,其中,i=1,2,3.......m,确定锅炉系统状态变量集合为H,其中第i个系统状态变量用Hi表示;
步骤二:将第k个锅炉系统中的系统状态进行细分,设细分后锅炉系统状态组合总数为K,第k个锅炉系统状态组合区间为Mk
步骤三:设定设备状态间的转移概率λ;
步骤四:设定设备工作状态集合Y与系统状态变量集合H的对应关系W;
步骤五:计算第i个设备在第t'时刻的设备状态到第t时刻的设备状态的概率P,其中t表示当前时刻,t'表示前一时刻,当前时刻t与前一时刻t'的时间间隔为单位时间Δt,第i个设备在单位时间Δt内的状态变化的概率相同,第i个设备在t'、t时刻的设备状态为
2.根据权利要求1所述的一种锅炉系统设备状态监测方法,其特征在于:所述步骤三具体包括以下步骤:设定第i个设备状态间的转移概率λ即第i个设备状态转移矩阵λi(i=1,2,3...,m),其计算公式如下:
<mrow> <msup> <mi>&amp;lambda;</mi> <mi>i</mi> </msup> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mn>11</mn> <mi>i</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mn>12</mn> <mi>i</mi> </msubsup> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mrow> <mn>1</mn> <mo>,</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> </mrow> <mi>i</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mn>21</mn> <mi>i</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mn>22</mn> <mi>i</mi> </msubsup> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mrow> <mn>2</mn> <mo>,</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> </mrow> <mi>i</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mrow> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>,</mo> <mn>1</mn> </mrow> <mi>i</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mrow> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>,</mo> <mn>2</mn> </mrow> <mi>i</mi> </msubsup> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msubsup> <mi>&amp;lambda;</mi> <mrow> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> </mrow> <mi>i</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,表示第i个设备从状态1变化到状态1的概率,表示第i个设备从状态1变化到状态2的概率,以此类推,表示第i个设备从态yi变化到状态yi的概率。
3.根据权利要求1所述的一种锅炉系统设备状态监测方法,其特征在于:所述步骤四具体包括以下步骤:该对应关系W为第i个设备在开启状态下对锅炉系统的n个系统状态变量的影响趋势,设定设备工作状态集合Y与系统状态变量集合H的对应关系W,其公式为:
<mrow> <mi>W</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>w</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>w</mi> <mn>12</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>w</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>w</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>w</mi> <mn>22</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>w</mi> <mrow> <mn>2</mn> <mi>j</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>w</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>w</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>w</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中wij表示第i个设备对第j个系统状态变量的影响值,所述j=1,2,3......n,wij的取值范围为{0,1,-1}。
4.根据权利要求1所述的一种锅炉系统设备状态监测方法,其特征在于:所述步骤五具体包括以下步骤:
步骤51:判断第i个设备在前一时刻t'是否处于正常工作状态,若是,则执行步骤53,否则执行步骤52;判断方法:设定第i个设备工作状态集合维数为Li,而设定第i个设备工作状态集合为Yi={y1,y2,y3,...,yn},其中n为自然数,则表示该第i个设备具有n个不同的工作状态,因此Li=n,若第i个设备在前一时刻t'的工作状态不属于集合Yi={y1,y2,y3,...,yn}中的正常工作状态,则判定该第i个设备不处于正常工作状态,异则反之,
步骤52:分别计算第i个设备在当前时刻t的状态的转移概率Pi,具体公式如下:
a)若当前时刻t第i个设备非正常工作:
<mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>|</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>&amp;lambda;</mi> <mrow> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> <mi>i</mi> </msubsup> <mo>&amp;times;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>;</mo> </mrow>
b)若当前时刻t第i个设备正常工作:
<mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>|</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>&amp;lambda;</mi> <mrow> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> <mi>i</mi> </msubsup> <mo>&amp;times;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>;</mo> </mrow>
步骤53:判断第i个设备在当前时刻t是否处于正常工作状态,若是,则计算第i个设备在当前时刻t的工作状态的转移概率Pi的公式如下:
<mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>|</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>&amp;Element;</mo> <mi>R</mi> <mrow> <mo>(</mo> <msup> <mi>M</mi> <mi>t</mi> </msup> <mo>|</mo> <msup> <mi>M</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msup> <mo>)</mo> </mrow> </mrow> </munder> <msubsup> <mi>&amp;lambda;</mi> <mrow> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> <mi>i</mi> </msubsup> <mo>&amp;times;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>;</mo> </mrow>
其中当前时刻t的锅炉系统状态为Mt,前一时刻t'的锅炉系统状态为Mt',R(Mt|Mt')是指设备在系统状态Mt'→Mt的允许变化状态集合,系统状态Mt'→Mt允许变化的状态集合与W的关系为:R(Mt|Mt')∈W,表示的是第t时刻第i个设备的工作状态,这里的允许变化状态集合R(Mt|Mt')是按照且R(Mx|Mx')∈W的方式进行的,
否则,判定第i个设备在当前时刻t的工作状态是否满足设备状态变化公式:
<mrow> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>&amp;Element;</mo> <mi>R</mi> <mrow> <mo>(</mo> <msup> <mi>M</mi> <mi>t</mi> </msup> <mo>|</mo> <msup> <mi>M</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msup> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
若是,则执行步骤54,否则,执行步骤55,
步骤54:计算第i个设备在当前时刻t的工作状态的转移概率Pi的公式如下:
<mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>|</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>&amp;lambda;</mi> <mrow> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> <mi>i</mi> </msubsup> <mo>&amp;times;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>,</mo> </mrow>
步骤55:计算第i个设备在当前时刻t的状态的转移概率Pi的公式如下:
<mrow> <msub> <mi>P</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>|</mo> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>|</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>t</mi> <mo>&amp;prime;</mo> </msup> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <mn>0.</mn> </mrow> 2
CN201410468800.6A 2014-09-15 2014-09-15 一种锅炉系统设备状态监测方法 Active CN104238516B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410468800.6A CN104238516B (zh) 2014-09-15 2014-09-15 一种锅炉系统设备状态监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410468800.6A CN104238516B (zh) 2014-09-15 2014-09-15 一种锅炉系统设备状态监测方法

Publications (2)

Publication Number Publication Date
CN104238516A CN104238516A (zh) 2014-12-24
CN104238516B true CN104238516B (zh) 2017-10-13

Family

ID=52226812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410468800.6A Active CN104238516B (zh) 2014-09-15 2014-09-15 一种锅炉系统设备状态监测方法

Country Status (1)

Country Link
CN (1) CN104238516B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109523027B (zh) * 2018-10-22 2021-01-05 新智数字科技有限公司 一种基于贝叶斯网络的锅炉运行数据监测方法和装置
CN114397867B (zh) * 2022-03-18 2022-06-10 山西正合天科技股份有限公司 一种基于物联网的工控机控制方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE522545C2 (sv) * 2001-03-06 2004-02-17 Goalart Ab System, anordning och förfarande för diagnostisering av flödesprocesser
US7016742B2 (en) * 2002-11-27 2006-03-21 Bahelle Memorial Institute Decision support for operations and maintenance (DSOM) system
CN102195362A (zh) * 2011-05-26 2011-09-21 中国电力科学研究院 一种计及系统运行可靠性的概率动态调度方法
CN102306178A (zh) * 2011-08-25 2012-01-04 盛乐信息技术(上海)有限公司 视频推荐方法及装置
CN103336428B (zh) * 2013-06-24 2016-06-08 广东电网有限责任公司电力科学研究院 一种电站锅炉优化方法和系统

Also Published As

Publication number Publication date
CN104238516A (zh) 2014-12-24

Similar Documents

Publication Publication Date Title
CN104267668B (zh) 基于贝叶斯方法的航天阀门零件加工过程故障诊断方法
Wang et al. An online fault diagnosis tool of VAV terminals for building management and control systems
EP3203338B1 (en) Plant performance evaluation apparatus, plant performance evaluation system, and plant performance evaluation method
WO2005111806A3 (en) Sensor fault diagnostics and prognostics using component model and time scale orthogonal expansions
Liu et al. A fault diagnosis method based on signed directed graph and matrix for nuclear power plants
CN104240779B (zh) 核电站堆芯功率象限倾斜的监督方法和装置
CN102460529A (zh) 装置异常监视方法以及系统
US20130318018A1 (en) Neural network-based turbine monitoring system
Jeong et al. A heuristic method of variable selection based on principal component analysis and factor analysis for monitoring in a 300 kw mcfc power plant
CN110112442B (zh) 一种燃料电池系统控制方法及装置
CN106200624B (zh) 基于交叉分段pca的工业锅炉故障诊断方法
CN104238516B (zh) 一种锅炉系统设备状态监测方法
CN102110187A (zh) 一种基于pca和人工免疫系统的流程工业混合故障诊断方法和系统
CN111679648B (zh) 一种基于高斯过程回归的多变量闭环控制回路性能评估方法
CN106289754B (zh) 一种基于统计控制技术的火电机组部件性能变化判定方法
CN111476386A (zh) 能源系统维护方法、装置、系统以及存储介质
CN105387449B (zh) 一种在锅炉蒸汽温度控制中使用二阶微分的控制方法
CN106050435B (zh) 用于燃气涡轮的调节和控制系统、计算机程序产品及方法
Li et al. Model-based monitoring and fault diagnosis of fossil power plant process units using group method of data handling
Jain et al. Prediction of remaining useful life of an aircraft engine under unknown initial wear
CN111666673A (zh) 一种锅炉过热器寿命的监控方法、装置、存储介质及设备
Dalapatu et al. Alarm allocation for event-based process alarm systems
CN106124933A (zh) 一种基于等价输入干扰法的电力系统故障诊断方法
CN103364199A (zh) 一种燃气轮机图形化故障诊断方法
CN104165352B (zh) 一种核电站蒸发器水位监控系统及其监控方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Miao Huifang

Inventor after: Wu Yichun

Inventor after: Zheng Jianxiang

Inventor after: Xie Shan

Inventor after: Zhao Xu

Inventor before: Zhao Xu

Inventor before: Miao Huifang

Inventor before: Wu Yichun

Inventor before: Xie Shan

Inventor before: Zheng Jianxiang

Inventor before: Chen Shaomin

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230420

Address after: Unit 02, 2nd Floor, Zone 2, Building 3, No. 2311-2 Heshan Road, Huli District, Xiamen City, Fujian Province, 361000

Patentee after: Peisen (Fujian) Environmental Technology Co.,Ltd.

Address before: 361000 Siming South Road, Xiamen, Fujian Province, No. 422

Patentee before: XIAMEN University