基于Markov模型与DS证据理论的安全仪表系统可靠性.docVIP

基于Markov模型与D-S证据理论的安全仪表系统可靠性 摘要：安全仪表系统（Safety Instrument System，SIS）的可靠性定量评估对其安全的提高有很重要的意思，基于Markov模型的安全仪表可靠性计算已经被大多数学者和研究机构使用，在Markov模型中都假设各个状态为一个固定的值，但在实际应用中每个状态均从“最佳”状态变为“次佳”状态。本文基于此问题将D-S证据理论(Dempster- Shafer theory of evidence)引入Markov模型，充分的运用到D-S证据理论中的信度函数和似真度函数，更加确切的反应各个状态的变化，并更准确的评估其安全可靠性。 1 引言 安全仪表系统是用来避免事故的发生或降低事故给设备、环境和人员造成危害的一种专用系统，主要应用于石化、核电、高铁等高新产业，在可靠性方面有较高的要求。安全仪表系统做为各种系统安全运行的最后一道屏障，其安全性由于自身系统结构、软件、硬件、维护及周围环境等原因无法完全保证，一旦其功能失效则会造成小到工厂设备停车，大到危害生命安全的灾难发生，所以对安全仪表系统的可靠性研究显得极其的重要。 安全仪表系统的可靠性主要是对安全仪表系统进行安全等级SIL(Safety Integrity Level)计算，而SIL与危险失效和安全失效有关，因此只需计算平均要求时失效概率和平均危险失效概率便可以评估SIS系统的可靠性。文献[1]、[2]用故障树计算安全仪表系统的可靠性，文献[3]中用可靠性框图计算，文献[4][5]中用简化公式法来计算，这些计算方法比较容易操作，但是存在多个性能指标需要多次建模、过多的条件假设、没法动态的反应系统各个状态间的变化。文献【6】、【7】、【8】等提出了用Markov模型来计算安全仪表系统的可靠性，该方法克服了上面几种方法存在的弊端，但是在计算时将Markov模型中系统各个状态值假设成为一个固定值，这不符合实际应用中设备的状态，在应用中一个系统状态值大多趋近“浴缸曲线”而不是一个固定值。 本文在此问题基础上将D-S证据理论引入Markov模型，用D-S证据理论中的信度函数和似真度函数计算安全仪表系统平均要求时失效概率的区间值，这样计算更加接近实际应用。 2 影响安全完整性的各种因素 失效模式 安全系统的失效模式分为2种:安全失效和危险失效。危险失效是指该失效会使安全系统处于危险状态或失去安全功能执行能力。安全失效是指该失效不会造成安全相关系统处于危险状态或失去安全功能执行能力,如造成安全系统误动作的失效。设某个设备的失效率为λ,其安全失效率和危险失效率分别为λS和λD,则有λ=λS+λD。在安全的前提下,误动作率越低越好。IEC 61508关注的重点是安全性(PFDavg)本文也主要计算安全仪表系统的PFDavg。 1.2 自诊断 安全设备往往具备自诊断功能。自诊断检测设备的运行状态,发现失效时报警,使设备能够尽快得到维修。自诊断不可能完全地检测到设备失效。因此,失效被分为检测到的和未检测到的。设总危险失效率为λD,检测到和未检测到的危险失效率分别为λDD和λDU。对安全失效做相似假设,则有: λD=λDD+λDU. 1.3 表决逻辑 安全系统用表决逻辑来决定冗余部分的输出,提高可靠性。如MooN,表示N个冗余中只要有M个或大于M个正常,那么整个冗余结构就能够正常完成其功能。 3 相关技术介绍 3.1 Markov模型 本节主要介绍Markov链的定义及Markov模型。 定义1.设是定义在概率空间上的随机序列。P表示概率测度是从集合F到实数域R的一个函数，集合F中的每个时间都被该函数赋予0到1之间的一个概率值，对任意以及状态当时，如果有 （1） 则称随机序列为Markov链，式（1）为Markov性或无后行[5]。 用Markov链建立模型可准确的计算出影响SIL的影响因素PFDavg等，通常Markov模型通过状态转移图来表示状态的变化。圆圈表示SIS的各个状态（包括正常状态、中间转换状态以及失效状态，）失效和维修的过程用一个带箭头的弧线表示。Markov模型状态转换图如图2所示： 图2 Markov模型状态转换图 Markov模型的求解方式很多，但工程上对通常使用离散时间矩阵相乘的方式来求其值，先以为基本时间单位，状态转移的概率是维修率或失效率组成的转移矩阵和的乘积（或）。根据Markov模型的状态转移图可得转移矩阵P。若SIS的初始状态为，则是经过一个后SIS各个状态，同理经过n个之后SIS各个状态为。下面以1oo1逻辑结构（如图3所示）的Markov模型为例进行计算。 图3 1oo1逻辑结构的Markov模型 由图可得 N维初始向量，n为系统状态个数，设第n-