CN107545375A - 一种石油石化高危设备故障风险定量化评价方法 - Google Patents

Abstract

本发明提供了一种石油石化高危设备故障风险定量化评价方法，包括如下步骤：S1：对石油石化高危设备中故障模式的风险因素进行分析，得到主要故障模式的层次结构；S2：建立高危设备系统的模糊故障树分析模型和模糊事件树分析模型；S3：根据模糊故障树分析模型和模糊事件树分析模型，建立模糊Bow‑tie Model；S4：确定故模糊障树和模糊事件树分析模型中的各事件发生概率；S5：基于模糊Bow‑tie Model的模糊风险和灵敏度分析。通过本发明，实现了风险分析中数据缺失或不确定情况下的风险定量化评价，计算出关键事件的风险系数和风险排序，可以指导我们采取适当的设备维修及风险管控措施，从而有助于增强高危设备的可靠性，降低重大风险事故的发生概率。

Description

一种石油石化高危设备故障风险定量化评价方法

技术领域

本发明涉及风险评估领域，具体涉及一种石油石化高危设备故障风险定量化评价方法。

背景技术

随着石油石化行业的迅速发展和国际竞争的日益激烈，实现零事故安全生产成为相关企业的宗旨。而对高危设备故障进行风险评价与风险管控是石油石化作业与生产中减少人员伤亡、环境污染、经济损失等重大风险事故的主要技术手段之一。但在高危设备风险定量评价过程中，由于原因事件和结果事件的相关数据存在丢失性、不确定性或未知性，致使有一部分事件的数据可以通过实测和历史数据定量计算得到；而另一部分只能通过定性的方法用语言变量进行评估。上述同时包含定量化和定性化数据的这一类问题可统称为混合型风险分析问题，然而目前现有的相关技术还没有针对石油石化高危设备故障风险定量化评价这类问题的具体方法。

针对上述问题，本发明提出了一种新的算法，并通过实际测试，验证了该算法的有效性和适用性。

发明内容

针对背景技术及现有设计的不足之处，本发明公开了一种石油石化高危设备故障风险定量化评价方法，所述方法包括如下步骤：

S1：对石油石化高危设备中故障模式的风险因素进行分析，得到主要故障模式的层次结构；

S2：建立高危设备系统的模糊故障树分析模型和模糊事件树分析模型；

S3：根据模糊故障树分析模型和模糊事件树分析模型，建立模糊Bow-tie Model；

S4：确定模糊故障树分析模型中的各事件发生概率，确定模糊事件树分析模型中的各事件发生概率；

S5：基于模糊Bow-tie Model的模糊风险和灵敏度分析，实现风险定量化评价。

优选的，所述步骤S1中包括，从设备的机理研究出发，对其进行潜在失效模式与后果分析；然后通过现场调研和专家经验，得到主要故障模式的层次结构。

优选的，所述步骤S4中，确定模糊故障树分析模型中的各事件发生概率的步骤为：

步骤S411：选择一个顶部事件，即一个风险因素，并建立逻辑故障树图表。

步骤S412：把故障树逻辑图的任意元素，即基本事件，划分到模糊事件的已知事件和主观语言评价的概率分析上，即可能性分析，判断基本事件是否均是已知的，即确定型数据，若是，转步骤S413；否则，转步骤S415；

步骤S413：计算基本事件发生概率；

步骤S414：计算顶事件发生概率；

步骤S415：判断基本事件是否均是未知的，若是，转步骤S417；否则，转步骤S416；

步骤S416：将确定型数据转化为模糊数；

步骤S417：对未知的基本事件，即模糊事件进行语言评估；

步骤S418：整合专家的意见，利用隶属函数把语言表达转化成模糊数；

步骤S419：通过和/或获得可能性故障概率；

步骤S4110：转换模糊数的故障率为模糊可能性评分；

步骤S4111：分析和解释结果；

其中，表示基本事件i的故障概率，表示以为基本事件的顶事件的概率，对应事件所对应的AND和OR分别表示逻辑门“与”、“或”关系，n是和逻辑门“与”有关的基本事件的数量。

优选的，所述步骤S416中，将确定型数据转化为模糊数P的具体方法为：用三角模糊数等价表示；

所述步骤S418中，利用隶属函数把语言表达转化成模糊数的具体方法包括：

依据隶属函数将语言定义转变为模糊数，其具体计算方式为，对于语言集I＝(i₀,i₁,…,i_m,…,i_n)，代表一组有序的语言评价值的集合，其中m,n∈N，且0＜m＜n，i_m为该语言集中的一个语言评价结果，则评价结果中的模糊数为((m-1)/n,m/n,(m+1)/n)；

所述步骤S418中，利用隶属函数把语言表达转化成模糊数的具体方法包括：

采用权数法进行分析，公式如下所示：

这里，A_ij是专家j为事件i的语言表达，m是事件的个数，n是专家的人数。W_j是专家j的权重因素，M_i代表事件i组合后的模糊数。

优选的，所述步骤S419中，转换模糊数的故障率为模糊可能性评分的具体方法为：将模糊数中的三个模糊参数相加，再除以3，得到最终的模糊可能性评分。

优选的，若所有的基本事件数据均是确定型，则采用计算；若基本事件数据是混合型，则采用计算。

优选的，所述步骤S4中，确定模糊事件树分析模型中的各事件发生概率的具体步骤包括：

步骤S421：对于已经识别的基本事件，找到可能导致的后果集合和没有后果的状态，从而构建事件树逻辑图；

步骤S422：用三角模糊数定义基本事件的可能性；

步骤S423：根据三角模糊数的运算规则，计算所有故障树的路径从而确定每个结果事件的概率；

步骤S424：对事件树的结果进行去模糊化，得到结果事件的概率。

步骤S425：分析和解释结果。

优选的，所述步骤S423的具体方法为：

对于已知的基本事件，用模糊数等价表示，其中P为基本事件的概率；对于未知的基本事件，依据隶属函数将语言定义转变为模糊数，其具体计算方式为，对于语言集 I＝(i₀,i₁,…,i_m,…,i_n)，代表一组有序的语言评价值的集合，其中m,n∈N，且0＜m＜n，i_m为该语言集中的一个语言评价结果，则评价结果中的模糊数为((m-1)/n,m/n,(m+1)/n)；再采用权数法进行分析，公式如下所示：

这里，A_ij是专家j为事件i的语言表达，m是事件的个数，n是专家的人数。W_j是专家j的权重因素，M_i代表事件i组合后的模糊数。

优选的，所述步骤S5中包括：将关键事件的可能性故障概率和结果事件的概率进行模糊乘积，计算基于模糊Bow-tie Model的模糊风险数再运用重心法去模糊化，即可得到最终的关键事件风险数。

优选的，所述步骤S5的灵敏度分析的具体方法为：

步骤S51：针对于某个关键事件，当所有的基本事件都存在时，可以得到模糊风险数去模糊化后，得到最终的风险数R。

步骤S52：通过消除第i个基本事件BE_i i＝(1,2,…,n)，再次运用公式进行计算，得到风险数R_i。

步骤S53：|R-R_i|表示消除第i个基本事件前后模糊数的偏差。比较偏差|R-R_i|，其值越大，则表示该基本事件的敏感度越强，即该基本事件对于关键事件越重要；相反，偏差越小，则表示该基本事件的敏感度越弱，即该基本事件对于关键事件越不重要。

步骤S54：根据偏差值的大小，对基本事件的敏感度进行排序，从而了解不同基本事件对于关键事件的重要度。

通过本发明的技术方案，实现了工业过程中数据缺失或不确定情况下的模糊风险分析，计算出关键事件的风险系数；可以通过风险系数及其风险排序，可以指导我们采取适当的预防措施，从而降低风险概率，有助于增强系统的稳定性。

附图说明

图1是本发明的具体实现流程；

图2是模糊故障树分析流程；

图3是模糊事件树分析流程；

图4是Bow-tie Model元素关系图；

图5是泥浆泵的工作原理；

图6是泥浆泵主要故障模式的层次结构图；

图7是十字头总成故障”R31为顶事件的故障树逻辑图；

图8是“连杆系统故障”R3为顶事件的故障树逻辑图；

图9是关键事件为“十字头总成故障”的事件树逻辑图；

图10是关键事件为“连杆系统故障”的事件树逻辑图。

具体实施方式

为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图说明本发明的具体实施方式。

如图1所示，本发明公开了一种石油石化高危设备故障风险定量化评价方法，所述方法包括如下步骤：

S1：对石油石化高危设备中故障模式的风险因素进行分析，得到主要故障模式的层次结构；

S2：建立高危设备的模糊故障树分析模型和模糊事件树分析模型；

S3：根据模糊故障树分析模型和模糊事件树分析模型，建立模糊Bow-tie Model；

S4：确定模糊故障树分析模型中的各事件发生概率，确定模糊事件树分析模型中的各事件发生概率；

S5：基于模糊Bow-tie Model的模糊风险和灵敏度分析，实现风险定量化评价。

优选的，所述步骤S1中包括，从设备的机理研究出发，对其进行潜在失效模式与后果分析；然后通过现场调研和专家经验，得到主要故障模式的层次结构。

优选的，所述步骤S4中，确定模糊故障树分析模型中的各事件发生概率的步骤为：

步骤S411：选择一个顶部事件，即一个风险因素，并建立逻辑故障树图表。

步骤S412：把故障树逻辑图的任意元素，即基本事件，划分到模糊事件的已知事件和主观语言评价的概率分析上，即可能性分析，判断基本事件是否均是已知的，即确定型数据，若是，转步骤S413；否则，转步骤S415；

步骤S413：计算基本事件发生概率；

步骤S414：计算顶事件发生概率；

步骤S415：判断基本事件是否均是未知的，若是，转步骤S417；否则，转步骤S416；

步骤S416：将确定型数据转化为模糊数；

步骤S417：对未知的基本事件，即模糊事件进行语言评估；

步骤S418：整合专家的意见，利用隶属函数把语言表达转化成模糊数；

步骤S419：通过和/或获得可能性故障概率；

步骤S4110：转换模糊数的故障率为模糊可能性评分；

步骤S4111：分析和解释结果；

其中，表示基本事件i的故障概率，表示以为基本事件的顶事件的概率，对应事件所对应的AND和OR分别表示逻辑门“与”、“或”关系，n是和逻辑门“与”有关的基本事件的数量。具体可以参见图2。

优选的，所述步骤S416中，将确定型数据转化为模糊数P的具体方法为：用三角模糊数等价表示；

所述步骤S418中，利用隶属函数把语言表达转化成模糊数的具体方法包括：

依据隶属函数将语言定义转变为模糊数，其具体计算方式为，对于语言集 I＝(i₀,i₁,…,i_m,…,i_n)，代表一组有序的语言评价值的集合，其中m,n∈N，且0＜m＜n，i_m为该语言集中的一个语言评价结果，则评价结果中的模糊数为((m-1)/n,m/n,(m+1)/n)；

所述步骤S418中，利用隶属函数把语言表达转化成模糊数的具体方法包括：

采用权数法进行分析，公式如下所示：

这里，A_ij是专家j为事件i的语言表达，m是事件的个数，n是专家的人数。W_j是专家j的权重因素，M_i代表事件i组合后的模糊数。

优选的，所述步骤S419中，转换模糊数的故障率为模糊可能性评分的具体方法为：将模糊数中的三个模糊参数相加，再除以3，得到最终的模糊可能性评分。

优选的，若所有的基本事件数据均是确定型，则采用计算；若基本事件数据是混合型，则采用计算。

优选的，所述步骤S4中，确定模糊事件树分析模型中的各事件发生概率的具体步骤包括：

步骤S421：对于已经识别的基本事件，找到可能导致的后果集合和没有后果的状态，从而构建事件树逻辑图；

步骤S422：用三角模糊数定义基本事件的可能性；

步骤S423：根据三角模糊数的运算规则，计算所有故障树的路径从而确定每个结果事件的概率；

步骤S424：对事件树的结果进行去模糊化，得到结果事件的概率。

步骤S425：分析和解释结果。如图3所示。

优选的，所述步骤S423的具体方法为：

对于已知的基本事件，用模糊数等价表示，其中P为基本事件的概率；对于未知的基本事件，依据隶属函数将语言定义转变为模糊数，其具体计算方式为，对于语言集 I＝(i₀,i₁,…,i_m,…,i_n)，代表一组有序的语言评价值的集合，其中m,n∈N，且0＜m＜n，i_m为该语言集中的一个语言评价结果，则评价结果中的模糊数为((m-1)/n,m/n,(m+1)/n)；再采用权数法进行分析，公式如下所示：

这里，A_ij是专家j为事件i的语言表达，m是事件的个数，n是专家的人数。W_j是专家j的权重因素，M_i代表事件i组合后的模糊数。

优选的，所述步骤S5中包括：将关键事件的可能性故障概率和结果事件的概率进行模糊乘积，计算基于模糊Bow-tie Model的模糊风险数再运用重心法去模糊化，即可得到最终的关键事件风险数。

优选的，所述步骤S5的灵敏度分析的具体方法为：

步骤S51：针对于某个关键事件，当所有的基本事件都存在时，可以得到模糊风险数去模糊化后，得到最终的风险数R。

步骤S52：通过消除第i个基本事件BE_i i＝(1,2,…,n)，再次运用公式进行计算，得到风险数R_i。

步骤S53：|R-R_i|表示消除第i个基本事件前后模糊数的偏差。比较偏差|R-R_i|，其值越大，则表示该基本事件的敏感度越强，即该基本事件对于关键事件越重要；相反，偏差越小，则表示该基本事件的敏感度越弱，即该基本事件对于关键事件越不重要。

步骤S54：根据偏差值的大小，对基本事件的敏感度进行排序，从而了解不同基本事件对于关键事件的重要度。

关于部分符号的计算规则，即加法乘法和减法Θ，定义为：

设两个任意的三角模糊数：和并且l₁,m₁,u₁,l₂,m₂,u₂均为实数，则

加法运算规则：

乘法运算规则：

减法运算规则：

关于语言变量的模糊化方法，由于风险分析的语言评估具有很强的主观性质以及数据信息的普遍缺乏，因此风险参数，即风险因素准确发生的可能性和后果严重性的数据测量通常是很困难的。合理以及合适的表达这些参数的方式常常是用定性的语言变量，特别是根据专家的判断来估计事件发生的可能性。例如：主观的语言变量像非常低，低，中，高和非常高等这样的变量经常被使用。此外，评估后果的严重性可能使用变量如轻微，次要，中度，关键和灾难性等。

上述主观的语言变量可以依据隶属函数进一步定义，比如三角形模糊数，有语言集 I＝(i₀,i₁,…,i_m,…,i_n)代表一组有序的语言评价值的集合，其中m,n∈N,且0＜m＜n，i_m为该语言集中的一个语言评价结果，那么在该评价结果中的三角模糊数可以表示为 ((m-1)/n,m/n,(m+1)/n)。例如，评价值“中等”是语言集(非常低，低，中等，高，非常高) 中的一个评价值，则令“非常低＝0，低＝1，中等＝2，高＝3，非常高＝4”,此时m＝2,n＝4，那么可以计算得到属于评价值“中等”的三角模糊数：(0.25,0.50,0.75)。如下所示：

本发明所述的Bow-tie Model如图4所示，模糊Bow-tie Model包括5个基本元素，图4 展示了这些元素之间的基本关系。

(1)原因：导致事件失效、故障、或者一个组件人为错误的根本原因，被称为基本事件。

(2)故障树：故障树图形表示导致最不期望发生事件或待研究事件的因果关系路径。最不期望发生事件或待研究事件是顶事件和不同原因的相互作用，通过基本事件、中间事件和逻辑描述门表示。

(3)关键事件：在蝴蝶结图表中，故障树中的顶事件是事件树中的初始事件，在蝴蝶结图表中称为关键事件。

(4)事件树：事件树序列可能产生关键事件中安全功能(如报警，自动停机)或意外因素(如起火、爆炸)的判断障碍(如成功/失败，正确/错误)。

(5)结果事件：通过模糊Bow-tie Model，从关键事件中得到的最后结果称为结果事件。

为了更好的解释本发明，下面以泥浆泵为例进行阐述。

泥浆泵是由动力段和液力端两大部分组成。动力端的作用是将动力机的回转运动转变为活塞的直线往复运动。它包括传动离合装置、变速减速装置和曲柄连杆。液力端油泵头体、缸套、活塞杆吸入阀和排出阀等组成，它的作用是通过活塞在缸套中作往复运动形成液缸容腔变化，从而完成能量转化，实现吸入和排出液体。泥浆泵的原理如图5所示，附图中的各部件分别为：1.滤水器，2.吸入阀，3.泵缸，4.活塞，5.活塞杆，6.十字头，7.连杆，8.曲柄轴，9.曲柄，10.排出阀，11.排出管道。

首先从泥浆泵系统的机理研究出发，对其进行潜在失效模式与后果分析，根据功能性质划分了6个子系统作为泥浆泵的1级故障模式R₁—R₆；再依据各子系统所属部件划分其2级故障模式共28个R₁₁—R₆₆；然后通过现场调研和专家经验，得到泥浆泵主要故障模式的层次结构如图6所示。

其中，泥浆泵系统中属于1级故障模式有提供扭矩的齿轮系统R₁，提供往复力的轴承系统R₂，提供直线运动的连杆系统R₃，提供吸排作用的活塞系统R₄，提供吸入效果的阀盖系统 R₅以及提供安全功能的压力系统R₆。而每个1级故障模式都可以继续深入分析其下属部件，再次分解为多个2级故障模式如图6所示的R₁₁-R₆₆。

为了建立泥浆泵系统的模糊故障树分析模型，首先对图6中的所有2级故障模式R₁₁-R₆₆，分别以各自的故障模式为顶事件，下属的原因风险因素为基本事件建立其模糊故障树分析模型，得到所有2级故障模式原因风险因素的发生概率；然后，同理将所有1级故障模式R₁-R₆，分别以各自的故障模式为顶事件，下属的2级故障模式作为基本事件建立其模糊故障树分析模型，由此得到所有1级故障模式原因风险因素的发生概率；最后对泥浆泵系统的主要故障模式进行原因风险分析并对主要故障模式的原因风险因素发生概率进行排序。

下面分别以泥浆泵系统中的2级故障模式十字头总成故障R₃₁和1级故障模式连杆系统故障为例R₃建立以R₃₁和R₃为顶事件的模糊故障树分析模型。十字头总成在泥浆泵系统中是起导向作用的连杆部件，也是泥浆泵系统中主要故障模式中的重要风险因素之一。在实际生产中，十字头总成是泥浆泵的关键设备，易损件多，工况复杂，有很多因素会造成十字头总成功能的失效，例如十字头和拉杆连接失效、十字头拉杆连杆螺栓损坏、十字头和导板间间隙大、十字头不对称、导板拉断、轴承故障、拉杆变形、卡箍变形以及轴承位损坏等等。

因此，为了合理选择以“十字头总成故障”即R₃₁为顶事件的基本事件，通过现场调研以及相关专家的经验确认，最终选择了9个主要因素为基本事件，分别是十字头和拉杆连接失效BE_1、十字头拉杆连接螺栓损坏BE₂、十字头和导板间间隙大BE₃、十字头不对称BE₄、导板拉断BE₅、轴承故障BE₆、拉杆变形BE₇、卡箍变形BE₈以及轴承位损坏BE₉，从而建立了以“十字头总成故障”R₃₁为顶事件的故障树逻辑图，如图7所示。

连杆系统是泥浆泵的关键成分，也是非常重要的风险因素。它作为泥浆泵系统中的1级故障模式之一，下属它的2级故障模式即为其基本事件。因此，由图6可知，以“连杆系统故障”R₃为顶事件的基本事件依次为十字头总成故障BE₁，导板总成故障BE₂，拉杆故障BE₃以及拉缸密封件BE₄，从而建立了以“连杆系统故障”R₃为顶事件的故障树逻辑图，如图8 所示。

根据定义，模糊Bow-tie Model中故障树的顶事件和事件树的初始事件应该是一致的，即称为关键事件。因此，为了建立泥浆泵系统的模糊事件树分析模型，首先对图6中的所有2 级故障模式R_11-R₆₆，分别以各自的故障模式为关键事件建立模糊事件树分析模型，分析各自事件发生可能产生的所有后果，从而计算各自事件发生的后果严重度；然后，同理将所有1 级故障模式R₁-R₆，分别以各自的故障模式为关键事件建立模糊事件树分析模型，计算各自事件发生的后果严重度，得到所有1级故障模式后果严重度；最后对泥浆泵系统的主要故障模式进行结果风险分析并对主要故障模式的后果严重度进行排序。

下面同样分别以泥浆泵系统中的2级故障模式十字头总成故障R₃₁和1级故障模式连杆系统故障R₃为例建立以R₃₁和R₃为关键事件的模糊事件树分析模型。事件树逻辑图表示初始事件，即关键事件发生故障后所有可能导致的后果。通过现场勘查以及文献调研，我们得知在实际生产中，十字头总成发生故障所能导致的后果是不能正常提供往复力，直接导致泥浆泵停机；如果能够提供往复力，但可能造成机器运行不平稳，直接导致泥浆泵停机；当然，也有可能运行平稳进而泥浆泵正常运行这三种情况可能发生的情况。图9展现了关键事件为“十字头总成故障”的事件树逻辑图。

连杆系统R₃主要为泥浆泵提供直线运动的功能，它作为泥浆泵系统的1级故障模式之一，包含了2级故障模式十字头总成故障R₃₁。因此，R₃₁事件发生可能产生的所有结果，对于R₃事件也会同样存在。同样，通过现场勘查以及文献调研，我们得知在实际生产中，连杆系统故障所能导致的后果情况有四种。图10展现了关键事件为“连杆系统故障”的事件树逻辑图。

通过图7，我们可以得到顶事件为“十字头总成故障”的基本事件表如下：

通过现场调研以及专家询问，得知上述9个基本事件中BE₃、BE₄、BE₆—BE₉可以根据各自频数得到各自的发生概率。而由于实际生产中，数据的缺乏和丢失是不可避免的，上述基本事件中BE₁、BE₂和BE₅无法得到确定的发生概率值。

为了评估风险因素，即“十字头总成故障”，本发明采用专家打分法对基本事件BE₁、BE₂以及BE₅进行评价从而确定它们的发生可能性。为此邀请了3位现场经验丰富的专家进行判断，其中博士研究生1位，硕士研究生2位，这些专家都有相同的评分权重，并尊重彼此的评价结果。最后利用权数法整合专家的意见，利用隶属函数把语言表达转化成三角模糊数。

通过调研，得到6个基本事件的确定发生概率以及3个基本事件的模糊评估。基于此混合型问题，为了便于运用或进行计算，其中，将确定发生概率同样转换为三角模糊数，得到如下：

由图7可知，顶事件为“十字头总成故障”R₃₁是用逻辑门“或”连接基本事件。因此通过以及运用MATLAB软件，我们可以计算得到顶事件R₃₁的模糊发生概率记作再根据公式去模糊化，得到最终顶事件R₃₁的发生概率L_R31＝0.7487。同理，可以计算泥浆泵系统中1级故障模式“连杆系统故障”R₃下的所有2级故障模式的模糊发生概率，结果如下：

通过上表和图8以“连杆系统故障”R₃为顶事件的故障树逻辑图，运用公式及运用MATLAB软件，可以计算得到顶事件R₃的发生概率 L_R3＝0.9771。同理，可以计算泥浆泵主要故障模式中所有1级故障模式的发生概率，并进行原因风险排序。结果表明泥浆泵系统中1级故障模式原因风险排序依次为： R₃＝R₄>R5>R₁>R₂>R₆，风险概率依次为：0.9771、0.9771、0.9614、0.8296、0.6783和0.6518。因此，R₃和R₄即连杆系统和齿轮系统故障模式失效的概率最大，应加强检查维修的频率。

通过图9可知，关键事件为“十字头总成故障”的事件树逻辑图展现的所有可能发生的情况共有三种。同样通过调研，P₂和S₃可以得到确定的后果严重度数据，而P₂₁则无法获得确定型的后果严重度数据。因此，我们根据专家打分，转换为三角模糊数，从而得到事件树可能性的相关数据，如下所示：

由公式运用MATLAB软件，我们可以计算得到：

由此，我们可以得到因为“十字头总成故障”系统失效模糊概率为：

或

由于(1)式的数据均是确定型数据转换为模糊数从而进行的运算，因此若所有的基本事件数据均是确定型，则采用公式(1)的结果；而(2)式的数据是混合型，即同时包含确定型数据和模糊数，因此若基本事件数据是混合型，则采用公式(2)的结果。

针对泥浆泵系统中关键事件是“十字头总成故障”的事件树分析，其基本事件属于混合型问题，因此采用式子(2)的计算结果，从而得到关键事件为“十字头总成故障”R₃₁的模糊后果严重度再根据去模糊化，得到最终顶事件R₃₁的后果严重度S_(R31)＝0.4375。同理，可以计算泥浆泵系统中1级故障模式“连杆系统故障”R₃下的所有2级故障模式的模糊后果严重度，结果如下：

通过上表和图10以“连杆系统故障”R₃为关键事件的事件树逻辑图，运用及运用MATLAB软件，可以计算得到关键事件R₃的后果严重度S_R3＝0.5978。同理，可以计算泥浆泵主要故障模式中所有1级故障模式的后果严重度，并进行结果风险排序。结果表明泥浆泵系统中1级故障模式结果风险排序依次为： R₁>R₂>R₆>R₅>R₄>R₃，后果严重度依次为：0.9251、0.8618、0.8129、0.7596、0.6424和0.5978。因此，R₁即齿轮系统故障模式失效时造成的危害最大，应加强安全防护措施。

针对泥浆泵系统2级故障模式中关键事件为“十字头总成故障”R₃₁的模糊故障树分析结果和模糊事件树分析结果，得到以及因此，根据公式可以计算出R₃₁的模糊风险系数最后运用重心法去模糊化，得到最终的风险系数R_(R31)＝0.306。

同样，针对泥浆泵系统1级故障模式中关键事件为“连杆系统故障”R₃的模糊故障树分析结果和模糊事件树分析结果，得到L_R3＝0.9771以及S_R3＝0.5978。因此，根据公式R＝L×S (若故障树分析最终结果和事件树分析最终结果均是确定数)，可以计算出R₃的风险系数 R_(R3)＝L_(R3)×S_(R3)＝0.5841。

同理，我们可以得到泥浆泵系统中所有2级故障模式的模糊风险系数和1级故障模式的风险系数，并进行最终风险系数的排序结果如下所示：

从上表可知：基于模糊Bow-tie Model的模糊风险分析，泥浆泵系统1级故障模式中R₁的风险系数最高，为0.7675，其次是R₅、R₄、R₂、R₃和R₆，风险系数依次是0.7303、0.6277、0.5846、0.5841和0.5298。因此，可以重点对R₁故障模式采取适当的预防措施，从而降低风险概率。

实际工业生产中，为了了解部件对于系统的敏感性以及重要性，常常针对某一关键事件，对其基本事件进行灵敏度分析。下面以泥浆泵系统中2级故障模式“十字头总成故障”R₃₁作为关键事件为例，分析其基本事件对于十字头总成的重要性。

首先针对关键事件R₃₁，当所有的基本事件都存在时，已经计算得到其风险数为R_R31＝0.7279；然后消除关键事件中的第i基本事件(i＝1,2,…,9)，再次运用公式三角模糊数的运算规则、进行计算，重新得到新的风险数R_BEi；接着将新的风险数与原风险数进行比较，其偏差|R_BE(R31)-R_BEi|则表示基本事件i对关键事件的重要性，即偏差指数越大，则表示该基本事件的敏感度越强，即该基本事件对于关键事件越重要；相反，偏差越小，则表示该基本事件的敏感度越弱，即该基本事件对于关键事件越不重要；最后，根据偏差值的大小，对基本事件的敏感度进行排序，从而了解不同基本事件对于关键事件R₃₁的重要度。通过MATALB进行编程运算，得到关键事件为“十字头总成故障”中基本事件的重要度，如下所示。

由上表可知，第三个基本事件BE₃的偏差最大，即它对关键事件十字头总成最重要，其次是BE₅、BE₄；而BE₇、BE₁以及BE₅的重要度相对较低。

需要说明的是，对于前述的各个方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本申请并不受所描述的动作顺序的限制，因为依据本申请，某一些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于优选实施例，所涉及的动作和单元并不一定是本申请必须的。

在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详细描述的部分，可以参见其他实施例的相关描述。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、ROM、RAM等。

以上所揭露的仅为本发明较佳实施例而已，当然不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。

Claims (10)

1.一种石油石化高危设备故障风险定量化评价方法，其特征在于，所述方法包括：

S1：对石油石化高危设备中故障模式的风险因素进行分析，得到主要故障模式的层次结构；

S2：建立高危设备的模糊故障树分析模型和模糊事件树分析模型；

S3：根据模糊故障树分析模型和模糊事件树分析模型，建立模糊Bow-tie Model；

S4：确定模糊故障树分析模型中的各事件发生概率，确定模糊事件树分析模型中的各事件发生概率；

S5：基于模糊Bow-tie Model的模糊风险和灵敏度分析，实现风险定量化评价。

2.一种如权利要求1所述的石油石化高危设备故障风险定量化评价方法，其特征在于，所述步骤S1中包括，从设备的机理研究出发，对其进行潜在失效模式与后果分析；然后通过现场调研和专家经验，得到主要故障模式的层次结构。

3.一种如权利要求1所述的石油石化高危设备故障风险定量化评价方法，其特征在于，所述步骤S4中，确定模糊故障树分析模型中的各事件发生概率的步骤为：

步骤S411：选择一个顶部事件，即一个风险因素，并建立逻辑故障树图表。

步骤S412：把故障树逻辑图的任意元素，即基本事件，划分到模糊事件的已知事件和主观语言评价的概率分析上，即可能性分析，判断基本事件是否均是已知的，即确定型数据，若是，转步骤S413；否则，转步骤S415；

步骤S413：计算基本事件发生概率；

步骤S414：计算顶事件发生概率；

步骤S415：判断基本事件是否均是未知的，若是，转步骤S417；否则，转步骤S416；

步骤S416：将确定型数据转化为模糊数；

步骤S417：对未知的基本事件，即模糊事件进行语言评估；

步骤S418：整合专家的意见，利用隶属函数把语言表达转化成模糊数；

步骤S419：通过和/或获得可能性故障概率；

步骤S4110：转换模糊可能性概率为模糊可能性评分；

步骤S4111：分析和解释结果；

其中，表示基本事件i的故障概率，表示以为基本事件的顶事件的故障概率，对应事件所对应的AND和OR分别表示逻辑门“与”、“或”关系，n是和逻辑门“与”有关的基本事件的数量。

4.一种如权利要求3所述的石油石化高危设备故障风险定量化评价方法，其特征在于，

所述步骤S416中，将确定型数据转化为模糊数P的具体方法为：用三角模糊数等价表示；

所述步骤S418中，利用隶属函数把语言表达转化成模糊数的具体方法包括：

依据隶属函数将语言定义转变为模糊数，其具体计算方式为，对于语言集I＝(i₀,i₁,…,i_m,…,i_n)，代表一组有序的语言评价值的集合，其中m,n∈N，且0＜m＜n，i_m为该语言集中的一个语言评价结果，则评价结果中的模糊数为((m-1)/n,m/n,(m+1)/n)；

所述步骤S418中，利用隶属函数把语言表达转化成模糊数的具体方法包括：

采用权数法进行分析，公式如下所示：

<mrow> <msub> <mi>M</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>W</mi> <mi>j</mi> </msub> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>n</mi> </mrow>

这里，A_ij是专家j为事件i的语言表达，m是事件的个数，n是专家的人数。W_j是专家j的权重因素，M_i代表事件i组合后的模糊数。

5.一种如权利要求3所述的石油石化高危设备故障风险定量化评价方法，其特征在于，

所述步骤S419中，转换模糊数的故障率为模糊可能性评分的具体方法为：将模糊数中的三个模糊参数相加，再除以3，得到最终的模糊可能性评分。

6.一种如权利要求3所述的石油石化高危设备故障风险定量化评价方法，其特征在于，

若所有的基本事件数据均是确定型，则采用计算；若基本事件数据是混合型，则采用计算。

7.一种如权利要求3所述的石油石化高危设备故障风险定量化评价方法，其特征在于，

所述步骤S4中，确定模糊事件树分析模型中的各事件发生概率的具体步骤包括：

步骤S421：对于已经识别的基本事件，找到可能导致的后果集合和没有后果的状态，从而构建事件树逻辑图；

步骤S422：用三角模糊数定义基本事件的可能性；

步骤S423：根据三角模糊数的运算规则，计算所有故障树的路径从而确定每个结果事件的概率；

步骤S424：对事件树的结果进行去模糊化，得到结果事件的概率。

步骤S425：分析和解释结果。

8.一种如权利要求7所述的石油石化高危设备故障风险定量化评价方法，其特征在于，

所述步骤S423的具体方法为：

对于已知的基本事件，用模糊数等价表示，其中P为基本事件的概率；对于未知的基本事件，依据隶属函数将语言定义转变为模糊数，其具体计算方式为，对于语言集I＝(i₀,i₁,…,i_m,…,i_n)，代表一组有序的语言评价值的集合，其中m,n∈N，且0＜m＜n，i_m为该语言集中的一个语言评价结果，则评价结果中的模糊数为((m-1)/n,m/n,(m+1)/n)；再采用权数法进行分析，公式如下所示：

<mrow> <msub> <mi>M</mi> <mi>i</mi> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>W</mi> <mi>j</mi> </msub> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>n</mi> </mrow>

这里，A_ij是专家j为事件i的语言表达，m是事件的个数，n是专家的人数。W_j是专家j的权重因素，M_i代表事件i组合后的模糊数。

9.一种如权利要求7所述的石油石化高危设备故障风险定量化评价方法，其特征在于，

所述步骤S5中包括：将关键事件的可能性故障概率和结果事件的概率进行模糊乘积，计算基于模糊Bow-tie Model的模糊风险数 再运用重心法去模糊化，即可得到最终的关键事件风险数。

10.一种如权利要求9所述的石油石化高危设备故障风险定量化评价方法，其特征在于，

所述步骤S5的灵敏度分析的具体方法为：

步骤S51：针对于某个关键事件，当所有的基本事件都存在时，可以得到模糊风险数去模糊化后，得到最终的风险数R。

步骤S52：通过消除第i个基本事件BE_i i＝(1,2,…,n)，再次运用公式进行计算，得到风险数R_i。

步骤S53：|R-R_i|表示消除第i个基本事件前后模糊数的偏差。比较偏差|R-R_i|，其值越大，则表示该基本事件的敏感度越强，即该基本事件对于关键事件越重要；相反，偏差越小，则表示该基本事件的敏感度越弱，即该基本事件对于关键事件越不重要。

步骤S54：根据偏差值的大小，对基本事件的敏感度进行排序，从而了解不同基本事件对于关键事件的重要度。