CN103235881A - 一种基于最小割集的核反应堆故障监测系统 - Google Patents

一种基于最小割集的核反应堆故障监测系统 Download PDF

Info

Publication number
CN103235881A
CN103235881A CN2013101393272A CN201310139327A CN103235881A CN 103235881 A CN103235881 A CN 103235881A CN 2013101393272 A CN2013101393272 A CN 2013101393272A CN 201310139327 A CN201310139327 A CN 201310139327A CN 103235881 A CN103235881 A CN 103235881A
Authority
CN
China
Prior art keywords
node
fault
cut set
fault tree
nodes
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
Application number
CN2013101393272A
Other languages
English (en)
Other versions
CN103235881B (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.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
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 Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN201310139327.2A priority Critical patent/CN103235881B/zh
Publication of CN103235881A publication Critical patent/CN103235881A/zh
Application granted granted Critical
Publication of CN103235881B publication Critical patent/CN103235881B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Testing And Monitoring For Control Systems (AREA)

Abstract

一种基于最小割集的核反应堆故障监测系统,为核反应堆系统的故障诊断提供决策支持,该系统也能够应用于航空、航海、化工等相关领域。该系统主要由以下6个模块组成:模型解析器、任务分发器、并行机、结果收集器、割集解析器、故障诊断器。该系统具有存储效率高、计算速度快、容易部署与维护等特点。与现有技术相比,本发明利用并行机对故障树每层节点并行计算,较大程度的提高了计算速度和准确性,能够有效地对核反应堆故障进行诊断和监测。

Description

一种基于最小割集的核反应堆故障监测系统
技术领域
本发明涉及一种基于最小割集的核反应堆故障监测系统,属于复杂核反应堆系统可靠性与故障诊断领域。
背景技术
故障树分析(Fault Tree Analysis,FTA)方法最早是1964年由贝尔电话实验室,60年代初在航空航天工业中得到应用,此后有了迅速发展,并成为反应堆、化工等领域不可缺少的分析工具之一。故障树分析是在系统设计过程中通过对可能造成系统失效的各种因素进行分析,画出逻辑框图,进而确定系统故障原因的各种组合方式和发生概率,并采取相应的改进措施,是提高系统可靠性的一种设计分析方法。其中系统故障原因的各种组合方式就是故障树的最小割集,最小割集是用于寻找系统薄弱环节的最有效方法之一。基于最小割集的故障监测就是利用系统最小割集来辅助寻找系统的可能故障模式,因此最小割集的获取是该方法的基础。
随着故障树技术的广泛应用,故障树最小割集求法也得到了快速发展。目前故障树最小割集求法主要分为以下四类:(1)第一类是传统的上行法、下行法,这类方法由于对故障树的所有路径进行搜索,超大型复杂系统的故障树路径数高达几十万以上,计算速度基本无法忍受;(2)第二类是基于Petri网的方法、矩阵化方法,基于Petri网的求解方法实际上是下行法的变种,路径数量并没有减少,只是采用了关联矩阵的形式,计算效率依然很低;矩阵化方法的矩阵维度过大,导致一般的计算机内存无法存储,同时计算量也很大;(3)第三类是基于二元决策图、零压缩二元决策图的方法,这种方法尽管存储效率比上述方法较高,但是由于计算过程产生的中间结果较多,计算量非常大,计算速度仍然有待提高;(4)第四类是近年来提出的分支演绎法,分支演绎法基于一定规则选取事件组合,判断是否满足故障树结构函数,该方法的求解思路尽管与上述方法不同,但是求解速度依赖于事件组合选取的合理性,不合理情况下需要判断的次数大得惊人,其效率具有一定的不确定性。
针对大型复杂系统,如核反应堆系统的故障树模型,这些方法的计算速度和结果的精度仍然有待提高,尤其是针对核电厂、核反应实验装置这一类超大型复杂系统。经过对目前国际上主流的基于上述方法的代表性商用系统RiskSpectrum、CAFTA、XFTA进行测试,一个拥有12000个逻辑门和8000个底事件的超大核电厂故障树,需要计算10-20分钟不等,计算速度过慢,无法做到及时响应,影响对这些系统故障的及时甄别和诊断,同时由于计算量过大而不得不采取截断策略,导致了最终结果有一定误差,不够精确,影响对系统故障的判断。
发明内容
本发明的目的在于:克服现有技术的不足,提出一种基于最小割集的核反应堆故障监测系统,使计算速度和精度得到提高,从而提高了系统故障的判断效率。
本发明的技术方案如下:一种基于最小割集的核反应堆故障监测系统,如图1所示包括:
(1)模型解析器:对给定的核反应堆系统的常规故障树模型进行解析,构建一种分层结构故障树模型,分离出故障树每层节点列表,每层节点的具体信息包括:名称、节点类型、所属分层行号和父亲节点名称;
(2)任务分发器:将解析后的分层结构故障树模型,按照自底向上的顺序,依次将每层节点分发给不同的并行机节点进行处理,当该层节点全部处理完毕之后才能处理上一层节点,直到处理完顶层节点为止;
(3)并行机:配置一定数量的计算节点,每个计算节点上部署一个计算引擎,计算引擎针对任务分发器发送来的分层结构故障树模型某一层节点中的某一个节点进行处理,将以该节点为头节点的子故障树转换成零压缩二元决策图,并将这些零压缩二元决策图送至结果收集器;
(4)结果收集器:收集并行机得到的计算结果,即零压缩二元决策图,当并行机处理的不是顶节点时,此结果收集器将计算结果反馈给任务分发器,当并行机处理的是顶节点是,此结果收集器将计算结果送至割集解析器;
(5)割集解析器:将结果收集器发送的零压缩二元决策图进行解析,得到所分析核反应堆系统的最小割集,并将该最小割集送至故障诊断器;
(6)故障诊断器:根据割集解析器发送过来的核反应堆系统最小割集,进行系统故障排除与诊断工作。
其具体实现步骤如下:
(1)安全分析人员对核反应堆系统进行故障树建模,得到核反应堆系统的常规故障树模型;
(2)模型解析器对步骤1得到的该核反应堆系统常规故障树模型解析称分层结构故障树模型;
(3)任务分发器自底向上将步骤2得到分层结构故障树模型中每一层节点的每一个节点发送个并行机进行计算,结果收集器将每一层计算结果反馈给任务分发器,当并行机处理完顶节点后,结果收集器将该计算结果送至割集解析器;
(4)割集解析器对步骤3中结果收集器送来的结果进行解析,得到核反应堆系统故障树的最小割集;
(5)故障诊断器根据步骤4得到的最小割集辅助展开系统故障排出与诊断工作。
所述模型解析器中对给定的核反应堆系统的常规故障树模型进行解析的具体实现方式如下:
先定义一种存储故障树结构的基本单元,该单元表示如下:<TOP,K,{INPUT1,INPUT2…INPUTn}>,上述基本单元中TOP代表每个子故障树的顶节点,K代表该节点的类型,INPUT1,INPUT2…INPUTn代表顶门的所有输入节点,n表示该节点的输入节点个数,当k=1<n时表示该节点为“或”门,当1<k<n时表示该节点为n取k表决门,当k=n时表示该节点为“与”门;
从故障树顶节点开始自上往下分层遍历,初始时给故障树顶节点赋予层号为1,并将该节点置于下标为1的数组中;将该节点的所有的INPUT层号设为2,记录这些节点的名称、节点类型、所属分层行号、父亲节点名称,并将这些节点置于下标为2的数组中;对所有INPUT节点进行递归操作,这些节点的所有INPUT层号则设为3,记录这些节点的名称、节点类型、所属分层行号、父亲节点名称,并将所有这些节点置于下标为3的数组中;依此类推,直到所有叶子节点全部加入相应的哈希表位置中为止。
从所述模型解析器中构建一种分层结构故障树模型到所述结果收集器中得到所分析核反应堆系统的最小割集,具体实现过程如下:从底往上,任务分发器先将模型解析器构建的分层结构故障树模型最底层节点依次发送到并行机的不同计算节点,进行并行计算,直到该层所有节点全部转化为零压缩二元决策图,这些零压缩二元决策图将反馈给任务分发器,任务分发器再连同这些结果与分层结构故障树模型从底向上倒数第二层节点一起分发给并行机进行并行计算,依此类推,直到将顶节点处理完毕为止,此时所得到的零压缩二元决策图即为整个故障树的零压缩二元决策图,该过程中每层节点的计算是并行的,相较于传统的自底向上或自顶向下的串行方法节省了计算时间。
所述并行机中配置一定数量的计算节点,该数量处于故障树每层节点数的平均值与最大值之间,K为故障树层数,且每层节点数自顶向下依次为N1,N2,...,Nk,所需的并行机总计算核数在区间
Figure BDA00003078944800031
max{N1,N2,...,Nk}]内。
本发明与现有技术相比优点在于:
(1)由于故障树是一个树形逻辑结构,父亲节点与孩子节点之间有较明显的输入输出关系,父亲节点的转化依赖于孩子节点的转化结果,因此故障树最小割集求解并行化的最大难点在于分析过程中独立子任务的划分。本发明通过模型解析器将常规故障树解析成分层故障树,克服了传统最小割集求解方法在划分独立子任务方面的不足;并采用新的处理流程,利用并行机对故障树每层节点进行并行计算。经过验证,该方法能得到更好的计算速度,对同样拥有12000个逻辑门和8000个底事件的超大核电厂故障树进行计算,仅需1分钟左右。
(2)由于并行机的部署大大提高了模型处理规模,不需要通过截断策略也能解决大规模计算难题,克服传统方法采用截断策略从而降低结果精度的缺点,使得计算结果更精确,从而较大程度提高了核电系统故障的监测效率。
(3)本发明通过合理配置并行机的计算核心数量,在计算速度与计算资源利用率上找到一个恰当的平衡点,在提高计算速度的基础上同时做到充分利用计算资源;同时本发明系统易于部署,使得成本相对可控;由于计算节点之间的独立性,使得本系统还具有较高的可靠性和易维护性,即单个节点故障不会对系统整体功能产生较大的影响,便于较短时间内恢复性能。
附图说明
图1是本发明的系统结构图;
图2是本发明的系统工作总流程图;
图3是本发明的模型解析器工作流程图;
图4是本发明的任务分发器工作流程图;
图5a是本发明的并行机工作流程图,图5b是其工作原理示意图;
图6是本发明的结果收集器工作流程图;
图7是本发明的割集解析器工作流程图;
图8是本发明的故障诊断器工作流程图;
图9是本发明的一个具体应用实例的常规故障树模型。
具体实施方式
为了更好地理解本发明,先对一些基本概念进行一下解释说明。
故障树:一种表达系统失效模型的树形结构,通过对系统失效机理进行逐层展开建模,直至无需或无法继续展开的基本部件为止;
最小割集:系统的故障模式,表示能导致系统失效的若干基本部件失效的组合,而且其中任何一个部件不失效则系统也不失效;
中间事件:故障树中逻辑门的输出即为中间事件;
顶事件:故障树的顶节点,通常为系统失效,是一种特殊的中间事件;
底事件:故障树的叶子节点,即系统故障树中无需或无法继续展开的基本部件;
ITE结构:一种逻辑表达式,是if-then-else的缩写,即表示ITE顶节点发生则取左支,不发生则取右支,如f=ite(x,f1,f0)表示若x发生(即x=1),则f=f1,若x不发生(即x=0),则f=f0,不同的算法对ITE语义的解释是不一样的,比如二元决策图的ITE结构表达式语义为f=ite(x,f1,f0)=x·f1+ˉx·f0,而零压缩二元决策图的ITE结构表达式语义为f=ite(x,f1,f0)=x·f1+f0;
零压缩二元决策图:一种有向无环图,作为将故障树转化为割集的中间状态,具有压缩率高,转换效率好的优点,其计算公式如下:
f=ite(x,f1,f0)=x·f1+f0       (1)
ite(x,L1,R1)·ite(x,L2,R2)=ite(x,(L1·L2+L1·R2+R1·L2),R1·R2)   (2)
ite(x,L1,R1)+ite(x,L2,R2)=ite(x,(L1+L2),(R1+R2))    (3)
ite(x,L1,R1)·ite(y,L2,R2)=ite(x,L1·h,R1·h)      (4)
ite(x,L1,R1)+ite(y,L2,R2)=ite(x,L1,(R1+h))       (5)
其中公式(1)、(6)分别表示二元决策图和零压缩二元决策图的基本语义,f是零压缩二元决策图,x是该图的顶事件,f1代表x的左分支,f0代表x的右分支。公式(2)-(5)和(7)-(10)分别表示两个零压缩二元决策图之间的“与”和“或”操作。以二元决策图算法为例,(2)、(3)表示顶事件相同的两个二元决策图的“与”和“或”操作,(4)、(5)表示顶事件不同的两个二元决策图的“与”和“或”操作,其中x,y代表故障树底事件,而且x的事件编号小于y,L1、R1分别代表第一个参与操作的二元决策图的顶事件的左分支,L2、R2分别代表该二元决策图的顶事件的右分支,h代表第二个参与操作的顶事件编号偏大的二元决策图,即ite(y,L2,R2),“·”代表“与”操作,“+”代表“或”操作。
下面结合附图对本发明进行详细说明。
如图1所示,本发明面向核反应堆系统故障诊断的故障树最小割集并行求解系统包括模型解析器、任务分发器、并行机、结果收集器、割集解析器、故障诊断器6大模块。
图2是该系统的工作总流程图,整个实现过程如下:
(1)首先获取核反应堆系统常规故障树模型,经过模型解析器的处理,构建一种分层的故障树结构;
(2)基于上一步得到的分层故障树结构,自底向上逐层处理,由任务分发器将每一层节点发送并行机计算,若非最终结果,则由结果收集器将计算结果反馈给任务分发器进行下一步分发,直到处理完顶层节点为止;
(3)处理完顶层节点之后,将所得到的零压缩二元决策图发往割集解析器,获取系统的最小割集;
(4)将得到的最小割集发往故障诊断器,对系统发生的故障进行排查和甄别,为故障诊断提供决策支持。
上述各模块的具体实现过程如下:
1.模型解析器
该模块的实现流程如图3所示:
先定义一种存储故障树结构的基本单元,该单元表示如下:<TOP,K,{INPUT1,INPUT2,…,INPUTn}>,上述基本单元中TOP代表每个子故障树的顶节点,K代表该节点的类型,n表示该节点的输入节点个数,当k=1<n时表示该节点为“或”门,当1<k<n时表示该节点为n取k表决门,当k=n时表示该节点为“与”门。
根据外部程序构造的故障树模型,从故障树顶点开始进行先根遍历,获取分层信息,即每层对应一个一维节点数组;根据第一步的分层信息,从最下面一层开始构建故障树,拥有相同父亲的节点相邻的排列在一起;根据第k层的父亲信息,依次排列第k+1层的节点,并通过第k层节点的“父亲节点名称”属性与第k层关联起来;同时针对第k+1层节点按照第二步进行排列;如此不断重复第三步,直到故障树最上面一层。
具体的,初始时给故障树顶节点赋予层号为1,并将该节点置于下标为1的数组中;将该节点的所有的INPUT层号设为2,记录这些节点的名称、节点类型、所属分层行号、父亲节点名称,并将这些节点置于下标为2的数组中;对所有INPUT节点进行递归操作,这些节点的所有INPUT层号则设为3,记录这些节点的名称、节点类型、所属分层行号、父亲节点名称,并将所有这些节点置于下标为3的数组中;依此类推,直到所有叶子节点全部加入相应的哈希表位置中为止。
2.任务分发器
该模块的实现流程如图4所示:
从最下面一层开始,依次遍历一层的全部节点,如果是叶子节点,将其作为一个子任务放入任务列表;如果是非叶子节点,将其所有孩子的某种操作(操作类型与该节点的门类型一致)作为一个子任务放入任务列表;
将同一层的所有子任务,提交多核处理器机群计算,转换为零压缩二元决策图;
重复上面两个步骤,直到故障树最上面一层。
3.并行机及计算核心模块
该模块的实现流程图如图5a所示,一个转换实例如图5b所示:
从任务分发器的任务队列中接收第一个子任务,调用第一个计算核心,下一个子任务调用下一个计算核心,直到该层故障树任务队列中子任务全部处理完毕,或者计算核心使用完毕;
对于每个计算核心而言,若当前节点为叶子节点,直接将其用ITE结构表示,若x为某个叶子节点,那么对应的ITE结构为ite(x,1,0),这个结构即为该叶子节点对应的零压缩二元决策图,将该图反馈给任务分发器,作为中间结果提供给下一轮子任务;
若当前节点为逻辑门,则调用任务分发器接收到的从结果搜集器返回的上一论计算结果,根据该门的类型以及前两个输入的编号顺序,选择零压缩二元决策图方法中相对应的转换公式进行转换,得到以该门为顶节点的子故障树对应的零压缩二元决策图,若该门不是故障树顶节点则将结果反馈给任务分发器,作为中间结果提供给下一轮子任务,否则发往割集解析器。
4.结果收集器
该模块的实现流程如图6所示:
结果收集器从并行机收到处理结果后,判断该结果是否为顶层节点的处理结果,若不是就反馈给任务分发器,这些结果在下一轮任务分发时和子任务一起发送给并行机,用于处理上层节点对应的子故障树;
如果是顶层节点的处理结果,即最终结果,就直接发送给割集解析器即可。
5.割集解析器
该模块的实现流程如图7所示:
割集解析器收到结果收集器发送过来的最终处理结果,即整个故障树对应的零压缩二元决策图,然后开始对改图进行解析。解析过程是从零压缩二元决策图的顶节点开始,寻找该节点到左子树所有叶子节点的路径,每条路径对应一条割集,且仅将出现在左支末端的节点加入该割集;
然后对顶节点的右支进行上面同样的操作;
最后直到所有节点均被解析完毕,得到整个系统的割集,并将割集最小化后得到的最小割集发送到故障诊断器提供故障诊断的决策支持。
6.故障诊断器
该模块的实现流程如图8所示:
故障诊断器接收到系统的最小割集后,根据已知系统状态信息对最小割集进行排除,基于以往诊断经验按照重要基本部件库中的部件重要性逐个对故障诊断工程师进行提问;
从重要基本部件库取出第一个部件,询问故障诊断工程师该部件是否故障,工程师可以回答“是”、“否”、“不确定”三种答案,若得到的答案是“是”则将该部件从所有最小割集中去除(若某个最小割集仅包含1个底事件即该部件,则表明找到故障!),并对处理后的结果进行割集最小化;若得到的答案是“否”,则将包含该部件的最小割集删掉;若得到的答案是“不确定”则不对最小割集进行任何操作,并取下一个部件重新进行上述操作,直到所有部件全部处理完;
对于不在重要基本部件库中的部件,故障诊断工程师可以根据自己掌握的确切信息(比如该部件确定故障或无故障)自定义部件库重复上述操作,从剩余的最小割集中继续进行筛选,直到无法进一步筛选为止;
故障诊断工程师可按照割集概率进行排序,从大到小逐个检查各个割集中出现的部件是否故障,知道排查出故障为止。
下面以全超导托卡马克核聚变实验装置EAST的氦气压缩段故障作为例子进行详细说明。图9为核能系统可靠性分析人员建立的EAST氦气压缩段常规故障树模型,其中氦气压缩段故障为故障树的顶事件,一级压缩故障、二级压缩故障、过滤纯化系统故障等为中间事件,X1-X16为底事件,并假设其节点顺序为X1<X2<…<X16。
模型解析器将该模型解析成分层结构故障树模型,得到4层节点列表,分别是:
第一层节点共1个,列表如下:
{“氦气压缩段故障”,OR,1,“NULL”};
第二层节点共3个,列表如下:
{“一级压缩故障”,OR,2,“氦气压缩段故障”};
{“二级压缩故障”,OR,2,“氦气压缩段故障”};
{“过滤纯化系统故障”,OR,2,“氦气压缩段故障”};
第三层节点共11个,列表如下:
{“压缩机1故障”,OR,3,“一级压缩故障”};
{“压缩机2故障”,OR,3,“一级压缩故障”};
{“压缩机3故障”,OR,3,“一级压缩故障”};
{“压缩机4故障”,OR,3,“二级压缩故障”};
……
{“手动阀外漏X13”,EVENT,3,“过滤纯化系统故障”};
{“过滤器A1730堵塞”,EVENT,3,“过滤纯化系统故障”};
{“过滤器故障”,OR,3,“过滤纯化系统故障”};
第四层节点共20个,列表如下:
{“压缩机1运行故障”,EVENT,4,“压缩机1故障”};
{“电源故障”,EVENT,4,“压缩机1故障”};
……
{“过滤器A1731堵塞”,EVENT,4,“过滤器故障”};
{“过滤器A1732堵塞”,EVENT,4,“过滤器故障”}
部署并行机的总计算核心数处于区间[(1+3+11+20)/4,20]即[9,20]之间,取计算核心数为16。
任务分发器先取第四层节点的20个任务发往并行机,由于20>16,因此先将20个任务中的前16个发往并行机,根据公式1进行计算,得到16个零压缩二元决策图,结果分别为:
ZBDD401=ite(X1,1,0);
ZBDD402=ite(X2,1,0);
ZBDD403=ite(X3,1,0);
ZBDD404=ite(X2,1,0);
……
ZBDD415=ite(X9,1,0);
ZBDD416=ite(X10,1,0);
该16个任务计算完毕后再分发剩下的4个任务,根据公式1进行计算,得到4个零压缩二元决策图,结果分别为:
ZBDD417=ite(X11,1,0);
ZBDD418=ite(X12,1,0);
ZBDD419=ite(X15,1,0);
ZBDD420=ite(X16,1,0);
结果收集器将第4层节点的20个结果反馈给任务分发器,以便进行第3层节点的任务分发。任务分发器将上20个结果与第3层节点的11个任务发往并行机,根据公式1-5进行计算,得到11个零压缩二元决策图,结果分别为:
ZBDD301=ZBDD401+ZBDD402=ite(X1,1,ite(X2,1,0));
ZBDD302=ZBDD403+ZBDD404=ite(X3,1,ite(X4,1,0));
……
ZBDD309=ite(X13,1,0);
ZBDD310=ite(X14,1,0);
ZBDD311=ZBDD419+ZBDD420=ite(X15,1,ite(X16,1,0));
结果收集器将第3层节点的11个结果反馈给任务分发器,以便进行第2层节点的任务分发。任务分发器将上11个结果与第2层节点的3个任务发往并行机,根据公式2-5进行计算,得到3个零压缩二元决策图,结果分别为:
ZBDD201=ZBDD301+ZBDD302+ZBDD303
=ite(X1,1,ite(X2,1,ite(X3,1,ite(X4,1,0))));
ZBDD202=ZBDD304+ZBDD305+ZBDD306+ZBDD307
=ite(X2,1,ite(X5,1,ite(X6,1,ite(X7,1,ite(X8,1,0)))));
ZBDD203=ZBDD308+ZBDD309+ZBDD310+ZBDD311
=ite(X9,1,ite(X10,1,ite(X11,1,ite(X12,1,ite(X13,1,ite(X14,1,ite(X15,1,ite(X16,1,0))))))));
结果收集器将第2层节点的3个结果反馈给任务分发器,以便进行第1层即顶节点的任务分发。任务分发器将上3个结果与顶节点发往并行机,根据公式2-5进行计算,得到最终结果,如下:
ZBDD101=ZBDD201+ZBDD202+ZBDD203
=ite(X1,1,ite(X2,1,ite(X3,1,ite(X4,1,ite(X5,1,ite(X6,1,ite(X7,1,ite(X8,1,ite(X9,1,ite(X10,1,ite(X11,1,ite(X12,1,ite(X13,1,ite(X14,1,ite(X15,1,ite(X16,1,0))))))))))))))));
由于是最终结果,结果收集器将该结果发往割集解析器。割集解析器对ZBDD101进行解析,得到16个最小割集,结果如下:
{“压缩机1运行故障X1”};
{“电源故障X2”};
……
{“过滤器A1732堵塞X16”};
割集解析器将上述16个割集发往故障诊断器提供故障诊断的决策支持,故障诊断器收到这16个割集后,第一步根据已知系统状态信息对最小割集进行排除,假设氦气压缩段系统的重要基本部件库及其重要性排序为“电源”,“压缩机2”,“过滤器F1721”,“手动阀”,故障诊断器依次询问故障诊断工程师上述4个重要部件是否故障,若“电源”,“压缩机2”,“过滤器F1721”,“手动阀”均没有故障,则可以排除割集中的{X2},{X3},{X9},{X13}4条割集;第二步根据剩下的割集概率进行排序,假设X1概率最大,取出X1询问故障诊断工程师该部件是否故障,故障诊断工程师进行检查后发现X1确实发生了,从而诊断系统故障为“压缩机1运行故障”。
本发明通过前期试验,结果表明不仅计算速度快,能对系统所发生故障进行快速甄别和排除,还具有很高的精度,为系统安全员的最终判断提供决策支持等。
本发明未详细描述的部分属于本领域公知技术。

Claims (4)

1.一种基于最小割集的核反应堆故障监测系统,其特征在于包括:
模型解析器:对给定的核反应堆系统的常规故障树模型进行解析,构建一种分层结构故障树模型,分离出故障树每层节点列表,每层节点的具体信息包括:名称、节点类型、所属分层行号和父亲节点名称;
任务分发器:将解析后的分层结构故障树模型,按照自底向上的顺序,依次将每层节点分发给不同的并行机节点进行处理,当该层节点全部处理完毕之后才能处理上一层节点,直到处理完顶层节点为止;
并行机:配置一定数量的计算节点,每个计算节点上部署一个计算引擎,计算引擎针对任务分发器发送来的分层结构故障树模型某一层节点中的某一个节点进行处理,将以该节点为头节点的子故障树转换成零压缩二元决策图,并将这些零压缩二元决策图送至结果收集器;
结果收集器:收集并行机得到的计算结果,即零压缩二元决策图,当并行机处理的不是顶节点时,此结果收集器将计算结果反馈给任务分发器,当并行机处理的是顶节点是,此结果收集器将计算结果送至割集解析器;
割集解析器:将结果收集器发送的零压缩二元决策图进行解析,得到所分析核反应堆系统的最小割集,并将该最小割集送至故障诊断器;
故障诊断器:根据割集解析器发送过来的核反应堆系统最小割集,进行系统故障排除与诊断工作。
2.根据权利要求1所述的基于最小割集的核反应堆故障监测系统,其特征在于:所述模型解析器中对给定的核反应堆系统的常规故障树模型进行解析的具体实现方式如下:
先定义一种存储故障树结构的基本单元,该单元表示如下:<TOP,K,{INPUT1,INPUT2…INPUTn}>,上述基本单元中TOP代表每个子故障树的顶节点,K代表该节点的类型,INPUT1,INPUT2…INPUTn代表顶门的所有输入节点,n表示该节点的输入节点个数,当k=1<n时表示该节点为“或”门,当1<k<n时表示该节点为n取k表决门,当k=n时表示该节点为“与”门;
从故障树顶节点开始自上往下分层遍历,初始时给故障树顶节点赋予层号为1,并将该节点置于下标为1的数组中;将该节点的所有的INPUT层号设为2,记录这些节点的名称、节点类型、所属分层行号、父亲节点名称,并将这些节点置于下标为2的数组中;对所有INPUT节点进行递归操作,这些节点的所有INPUT层号则设为3,记录这些节点的名称、节点类型、所属分层行号、父亲节点名称,并将所有这些节点置于下标为3的数组中;依此类推,直到所有叶子节点全部加入相应的哈希表位置中为止。
3.根据权利要求1所述的基于最小割集的核反应堆故障监测系统,其特征在于:从所述模型解析器中构建一种分层结构故障树模型到所述结果收集器中得到所分析核反应堆系统的最小割集具体实现过程如下:从底往上,任务分发器先将模型解析器构建的分层结构故障树模型最底层节点依次发送到并行机的不同计算节点,进行并行计算,直到该层所有节点全部转化为零压缩二元决策图,这些零压缩二元决策图将反馈给任务分发器,任务分发器再连同这些结果与分层结构故障树模型从底向上倒数第二层节点一起分发给并行机进行并行计算,依此类推,直到将顶节点处理完毕为止,此时所得到的零压缩二元决策图即为整个故障树的零压缩二元决策图,该过程中每层节点的计算是并行的,相较于传统的自底向上或自顶向下的串行方法节省了计算时间。
4.根据权利要求1所述的基于最小割集的核反应堆故障监测系统,其特征在于:所述并行机中配置一定数量的计算节点,该数量处于故障树每层节点数的平均值与最大值之间,K为故障树层数,且每层节点数自顶向下依次为N1,N2,...,Nk,所需的并行机总计算核数在区间
Figure FDA00003078944700021
max{N1,N2,...,Nk}]内。
CN201310139327.2A 2013-04-21 2013-04-21 一种基于最小割集的核反应堆故障监测系统 Expired - Fee Related CN103235881B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310139327.2A CN103235881B (zh) 2013-04-21 2013-04-21 一种基于最小割集的核反应堆故障监测系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310139327.2A CN103235881B (zh) 2013-04-21 2013-04-21 一种基于最小割集的核反应堆故障监测系统

Publications (2)

Publication Number Publication Date
CN103235881A true CN103235881A (zh) 2013-08-07
CN103235881B CN103235881B (zh) 2016-07-20

Family

ID=48883921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310139327.2A Expired - Fee Related CN103235881B (zh) 2013-04-21 2013-04-21 一种基于最小割集的核反应堆故障监测系统

Country Status (1)

Country Link
CN (1) CN103235881B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103400050A (zh) * 2013-08-21 2013-11-20 中国科学院合肥物质科学研究院 一种多人协同的核反应堆风险确定方法及系统
CN103440419A (zh) * 2013-08-31 2013-12-11 中国科学院合肥物质科学研究院 一种基于故障树和层次分析法的可靠性分配系统及分配方法
CN104268432A (zh) * 2014-10-16 2015-01-07 中国科学院合肥物质科学研究院 一种基于故障树分析的共因失效可靠性分配系统
CN104850750A (zh) * 2015-05-27 2015-08-19 东北大学 一种核电站反应堆保护系统可靠性分析方法
CN106168797A (zh) * 2016-05-25 2016-11-30 哈尔滨工程大学 一种模块化获取核电站故障树顶事项失效概率的方法
CN106548279A (zh) * 2016-11-02 2017-03-29 中广核研究院有限公司 一种用于先进反应堆总体设计阶段的安全分析方法
CN107808214A (zh) * 2017-10-22 2018-03-16 南京理工大学 车间制造系统的启发式二叉决策图变量序优化表示方法
CN109165124A (zh) * 2018-08-07 2019-01-08 南京翼辉信息技术有限公司 一种基于故障树的嵌入式系统硬件故障检测及处理方法
CN110009228A (zh) * 2019-04-04 2019-07-12 中国核动力研究设计院 概率论与确定论相结合的核电厂多样性保护系统设计方法
CN110045715A (zh) * 2019-04-15 2019-07-23 中国民航大学 基于Petri网与改进二元决策图模型的最小割集求解方法
CN110580559A (zh) * 2018-06-08 2019-12-17 大陆泰密克汽车系统(上海)有限公司 随机硬件失效指标获取方法及装置
CN115390550A (zh) * 2022-10-27 2022-11-25 潍柴雷沃智慧农业科技股份有限公司 联合收割机的故障诊断方法、系统、控制器和联合收割机

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100082318A1 (en) * 2008-09-29 2010-04-01 Korean Atomic Energy Research Institute Single quantification method of external event psa model containing multi-compartment propagation scenarios
CN102184291A (zh) * 2011-05-06 2011-09-14 北京航空航天大学 一种针对全测试无反馈系统的系统级故障诊断方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100082318A1 (en) * 2008-09-29 2010-04-01 Korean Atomic Energy Research Institute Single quantification method of external event psa model containing multi-compartment propagation scenarios
CN102184291A (zh) * 2011-05-06 2011-09-14 北京航空航天大学 一种针对全测试无反馈系统的系统级故障诊断方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
LISA M.BARTLETT等: "Choosing a Heuristic for the "Fault Tree to Binary Decision Diagram" Conversion,Using Neural Networks", 《IEEE TRANSACTIONS ON RELIABILITY》, vol. 51, no. 3, 30 September 2002 (2002-09-30), pages 344 - 349 *
刘萍: "大型故障树分析的算法研究与系统设计", 《万方数据企业知识服务平台》, 25 May 2011 (2011-05-25), pages 33 - 66 *
吴宜灿等: "大型集成概率安全分析软件系统的研究与发展", 《核科学与工程》, vol. 27, no. 3, 15 September 2007 (2007-09-15), pages 270 - 276 *
吴宜灿等: "秦山三期重水堆核电站风险监测器研发进展", 《核科学与工程》, vol. 31, no. 1, 15 March 2011 (2011-03-15) *
殷园等: "多核并行计算技术在风险监测系统中的应用", 《核科学与工程》, vol. 31, no. 2, 15 June 2011 (2011-06-15), pages 175 - 179 *
汪进等: "适于风险监测系统的零压缩二元决策图基本事件排序方法", 《核科学与工程》, vol. 30, no. 4, 15 December 2010 (2010-12-15), pages 360 - 364 *
王立柱: "《C/C++数据结构》", 31 August 2003, article "树的存储与遍历", pages: 272-278 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103400050A (zh) * 2013-08-21 2013-11-20 中国科学院合肥物质科学研究院 一种多人协同的核反应堆风险确定方法及系统
CN103400050B (zh) * 2013-08-21 2016-07-06 中国科学院合肥物质科学研究院 一种多人协同的核反应堆风险确定方法及系统
CN103440419A (zh) * 2013-08-31 2013-12-11 中国科学院合肥物质科学研究院 一种基于故障树和层次分析法的可靠性分配系统及分配方法
CN103440419B (zh) * 2013-08-31 2016-09-07 中国科学院合肥物质科学研究院 一种基于故障树和层次分析法的可靠性分配系统及分配方法
CN104268432A (zh) * 2014-10-16 2015-01-07 中国科学院合肥物质科学研究院 一种基于故障树分析的共因失效可靠性分配系统
CN104268432B (zh) * 2014-10-16 2017-02-08 中国科学院合肥物质科学研究院 一种基于故障树分析的共因失效可靠性分配系统
CN104850750A (zh) * 2015-05-27 2015-08-19 东北大学 一种核电站反应堆保护系统可靠性分析方法
CN104850750B (zh) * 2015-05-27 2017-10-24 东北大学 一种核电站反应堆保护系统可靠性分析方法
CN106168797A (zh) * 2016-05-25 2016-11-30 哈尔滨工程大学 一种模块化获取核电站故障树顶事项失效概率的方法
CN106548279A (zh) * 2016-11-02 2017-03-29 中广核研究院有限公司 一种用于先进反应堆总体设计阶段的安全分析方法
CN107808214A (zh) * 2017-10-22 2018-03-16 南京理工大学 车间制造系统的启发式二叉决策图变量序优化表示方法
CN107808214B (zh) * 2017-10-22 2022-05-27 南京理工大学 车间制造系统的启发式二叉决策图变量序优化表示方法
CN110580559A (zh) * 2018-06-08 2019-12-17 大陆泰密克汽车系统(上海)有限公司 随机硬件失效指标获取方法及装置
CN110580559B (zh) * 2018-06-08 2023-08-11 大陆泰密克汽车系统(上海)有限公司 随机硬件失效指标获取方法、装置、设备及存储介质
CN109165124A (zh) * 2018-08-07 2019-01-08 南京翼辉信息技术有限公司 一种基于故障树的嵌入式系统硬件故障检测及处理方法
CN110009228A (zh) * 2019-04-04 2019-07-12 中国核动力研究设计院 概率论与确定论相结合的核电厂多样性保护系统设计方法
CN110009228B (zh) * 2019-04-04 2021-04-02 中国核动力研究设计院 概率论与确定论相结合的核电厂多样性保护系统设计方法
CN110045715A (zh) * 2019-04-15 2019-07-23 中国民航大学 基于Petri网与改进二元决策图模型的最小割集求解方法
CN110045715B (zh) * 2019-04-15 2021-10-01 中国民航大学 基于Petri网与改进二元决策图模型的最小割集求解方法
CN115390550A (zh) * 2022-10-27 2022-11-25 潍柴雷沃智慧农业科技股份有限公司 联合收割机的故障诊断方法、系统、控制器和联合收割机

Also Published As

Publication number Publication date
CN103235881B (zh) 2016-07-20

Similar Documents

Publication Publication Date Title
CN103235881A (zh) 一种基于最小割集的核反应堆故障监测系统
CN107272646B (zh) 基于专家系统的压力机故障诊断系统
CN101976221B (zh) 一种基于粒子群禁忌组合的并行测试任务调度方法及平台
CN103399820A (zh) 一种基于事件树的序列和后果分析的故障判断系统及方法
CN107679146A (zh) 电网数据质量的校验方法和系统
CN110851499A (zh) 基于知识的装配工艺设计方法、系统及介质
CN105070097A (zh) 一种营运车辆道路交通事故风险识别方法
CN105095602B (zh) 一种基于Petri网的Web服务组合演化方法
CN108920609A (zh) 基于多维度分析的电力实验数据挖掘方法
CN108009270A (zh) 一种基于分布式内存计算的文本检索方法
CN103440419A (zh) 一种基于故障树和层次分析法的可靠性分配系统及分配方法
CN108681598A (zh) 任务自动重跑方法、系统、计算机设备和存储介质
Liu et al. A concurrent fault diagnosis method of transformer based on graph convolutional network and knowledge graph
CN104317244A (zh) 一种可重构制造系统零件族构建方法
CN102799960A (zh) 面向数据模型的并行业务流程异常检测方法
CN104765688B (zh) 一种可动态部署规则的Java程序运行时验证系统
CN106056299A (zh) 一种基于流程挖掘的决策支持方法及装置
CN103400050A (zh) 一种多人协同的核反应堆风险确定方法及系统
CN107918358A (zh) 数控装备故障分析方法和装置
CN106649124A (zh) 一种基于Actor模型的并行动态符号执行方法和系统
CN110007662A (zh) 自主保障技术验证系统全飞行器健康评估与管理地面推理平台
Wang et al. Feature ordering for neural incremental attribute learning based on Fisher's linear discriminant
CN102270190B (zh) 一种复杂决策问题的计算机建模与求解处理方法
Chen et al. Research on equipment situation display based on multi-source data fusion
CN103970854B (zh) Sap ecc端与sap bw端之间增量凭证信息同步方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160720

Termination date: 20210421

CF01 Termination of patent right due to non-payment of annual fee