《一种基于贝叶斯理论的石化设备失效率推断方法.pdf》由会员分享，可在线阅读，更多相关《一种基于贝叶斯理论的石化设备失效率推断方法.pdf（10页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN103336903A43申请公布日20131002CN103336903ACN103336903A21申请号201310272526022申请日20130701G06F19/0020110171申请人中国石油大学（华东）地址266580山东省青岛市青岛经济技术开发区长江西路66号72发明人王海清汪航尚胜美张玉涛刘俊芳刘芳54发明名称一种基于贝叶斯理论的石化设备失效率推断方法57摘要本发明公开了一种基于贝叶斯理论的石化设备失效率推断方法。以国外商业可靠性数据库作为先验信息，以国内的现场可靠性数据作为抽样信息，利用泊松分布模拟设备失效过程，并引入技术差距系数，通过贝叶斯方法导出动。
2、态的石化设备失效率后验分布，从而反映国内设备水平的不同失效模式下的失效率估计。企业可根据调整后的失效率信息对关键设备和易发失效模式进行分析，有针对性的进行预防维护，并对备件管理进行优化，对探索如何以较低的投入和代价建立适合国内生产操作和管理状况的可靠性数据库均具有积极意义。51INTCL权利要求书2页说明书4页附图3页19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书4页附图3页10申请公布号CN103336903ACN103336903A1/2页21一种基于贝叶斯理论的石化设备失效率推断方法，其特征在于包括如下步骤步骤一、确定样本似然函数；假定所研究的石化设备具有恒定的失效。
3、率，其寿命分布满足指数分布，则所抽取样本的似然函数为（10）样本似然函数中各参数的意义为W同一设备样本总数；D失效模式；D失效模式D的失效率；T观察间隔时间；KDW个样本失效模式D在T时间内出现次数，为非负整数；KDKD的阶乘；PXKD失效模式D出现KD次的概率；步骤二、确定先验分布密度函数；因为样本似然函数是参数为D的泊松分布，根据泊松分布均值D的共轭先验分布为伽马分布，故确定先验分布密度函数为（20）先验分布函数中各参数的意义为伽马分布中形状参数；伽马分布中尺度参数；伽马函数；DD的先验分布密度；步骤三、确定设备技术差距系数；当使用国外可靠性数据库作为先验信息时，为先验信息补充一个技术差距。
4、系数C以调整时间差对设备失效带来的影响，使用模糊综合估计法得到当前国内与国外设备技术差距的年限，即YA1Z1A2Z2ANZN（30）其中各参数的意义为N第N个专家；AN专家N的权重，由从业年数和业务能力确定，0AN1；ZN专家N估计的落后年数；将综合估计落后年数Y与设备数据获得的时间差比值作为技术差距系数C的近似估计，即（31）其中，C0，Y表示综合估计落后年数，Y1为获得先验数据的年限中值，Y2表示样本获得年限中值；步骤四、确定先验分布超参数；权利要求书CN103336903A2/2页3伽马先验分布GA（，）中有两个超参数，选择先验矩法确定超参数，可得伽马先验分布GA（，）的期望和方差，即（。
5、40）解之，可得超参数与的估计为（41）其中各参数意义为D先验分布中失效模式D的失效率；SD先验分布中失效模式D的标准差；步骤五、确定后验分布密度函数；由样本似然函数式（10）和先验分布密度函数式（20）可以写出W个样本和参数D的联合密度函数为（50）由贝叶斯公式可得D的后验分布密度函数为（51）其中，MKD是W个样本的边缘密度函数，与D无关，表示式子两边只差一个不依赖于D的常数因子，由后验分布可知失效率D的后验分布服从伽马分布，即GA（KD，WT），故此W个样本出现同一失效模式次数的后验密度为（52）由后验分布密度函数可得样本后验分布期望和标准差分别为（53）其中，WT表示W个样本的累计观察。
6、日历时间或累计观察工作时间，将公式（41）带入（53）可得经过贝叶斯方法调整后的设备失效率。权利要求书CN103336903A1/4页4一种基于贝叶斯理论的石化设备失效率推断方法发明领域0001本发明涉及一种基于贝叶斯理论的石化设备失效率推断方法，它主要应用于石化设备失效率推断领域。背景技术0002石化行业具有较高作业风险，做好安全技术和管理工作对于减少财产和人员损失具有重要意义。设备是一切石化生产的工具和基础，从石油钻采、运输、加工到化工成品的一系列过程中会用到各种各样的设备，设备本身的可靠性在安全生产过程中起着重要作用。针对设备可靠性，美国EXIDA公司发布了安全设备可靠性手册，美国CCP。
7、S发布了过程设备可靠性数据指南。目前，国际上最著名的设备可靠性数据库是挪威DNV发布的OREDAOFFSHORERELIABILITYDATA数据库。2009年，DNV发布了OREDA第5版，其中卷1为陆上设备，包含了260多种设施、16000多个设备、38000多种失效、68000多个维修记录，卷2为海底设备，该数据库包含了安全设备系统、电气设备、机电设备、机械装备、海底装备、勘探和生产设备的大量可靠性和维护数据。0003目前，我国对工业安全相关系统设备的可靠性数据信息的采集和整理工作的还没有系统展开，还没有建立可供使用的工业可靠性数据库，在可靠性数据的获取、计算与分析等方面的技术方法也缺乏。
8、。本发明以国外商业可靠性数据库为基础，利用贝叶斯方法间接调整得出反映特定生产状况的我国石化设备失效率估计，对于设备安全评估和安全策略的制定，以及对探索如何以较低的投入和代价建立适合国内生产操作和管理状况的可靠性数据库均具有积极意义。发明内容0004本发明的一个目的，是将国外已经积累的大量石化设备可靠性数据作为先验信息，以国内的现场可靠性数据作为抽样信息，通过贝叶斯方法调整得到反映国内设备水平的不同失效模式下的失效率估计。0005为了达到以上目的，本发明提供了一种基于贝叶斯理论的石化设备失效率推断方法，包括如下步骤0006步骤一、确定样本似然函数；0007假定所研究的石化设备具有恒定的失效率，其。
9、寿命分布满足指数分布，则所抽取样本的似然函数为0008（10）0009样本似然函数中各参数的意义为0010W同一设备样本总数；0011D失效模式；0012D失效模式D的失效率；说明书CN103336903A2/4页50013T观察间隔时间；0014KDW个样本失效模式D在T时间内出现次数，为非负整数；0015KDKD的阶乘；0016PXKD失效模式D出现KD次的概率；0017步骤二、确定先验分布密度函数；0018因为样本似然函数是参数为D的泊松分布，根据泊松分布均值D的共轭先验分布为伽马分布，故确定先验分布密度函数为0019（20）0020先验分布函数中各参数的意义为0021伽马分布中形状参数。
10、；0022伽马分布中尺度参数；0023伽马函数；0024DD的先验分布密度；0025步骤三、确定设备技术差距系数；0026当使用国外可靠性数据库作为先验信息时，为先验信息补充一个技术差距系数C以调整时间差对设备失效带来的影响，使用模糊综合估计法得到当前国内与国外设备技术差距的年限，即0027YA1Z1A2Z2ANZN（30）0028其中各参数的意义为0029N第N个专家；0030AN专家N的权重，由从业年数和业务能力确定，0AN1；0031ZN专家N估计的落后年数；0032将综合估计落后年数Y与设备数据获得的时间差比值作为技术差距系数C的近似估计，即00330034（31）0035其中，C0，。
11、Y表示综合估计落后年数，Y1为获得先验数据的年限中值，Y2表示样本获得年限中值；0036步骤四、确定先验分布超参数；0037伽马先验分布GA（，）中有两个超参数，选择先验矩法确定超参数，可得伽马先验分布GA（，）的期望和方差，即0038（40）0039解之，可得超参数与的估计为0040（41）说明书CN103336903A3/4页60041其中各参数意义为0042D先验分布中失效模式D的失效率；0043SD先验分布中失效模式D的标准差；0044步骤五、确定后验分布密度函数；0045由样本似然函数式（10）和先验分布密度函数式（20）可以写出W个样本和参数D的联合密度函数为0046（50）004。
12、7由贝叶斯公式可得D的后验分布密度函数为0048（51）0049其中，MKD是W个样本的边缘密度函数，与D无关，表示式子两边只差一个不依赖于D的常数因子，由后验分布可知失效率D的后验分布服从伽马分布，即GA（KD，WT），故此W个样本出现同一失效模式次数的后验密度为0050（52）0051由后验分布密度函数可得样本后验分布期望和标准差分别为0052（53）0053其中，WT表示W个样本的累计观察日历时间或累计观察工作时间，将公式（41）带入（53）可得经过贝叶斯方法调整后的设备失效率。附图说明0054图1是本发明的一个实施方法的流程图；0055图2是实施例中与国外装备技术水平差距年限的模糊综合。
13、估计；0056图3是实施例中OREDA按失效模式统计的压缩机部分失效率；0057图4是实施例中某石化企业经贝叶斯方法调整后压缩机部分失效率。具体实施方式0058以下结合附图和实施例对本发明的技术方案做进一步的说明。0059如附图3（OREDA编号为11111）所示为来自OREDA（数据采集时间19932000年，可知Y1为1996年6月）数据库的离心式电驱动压缩机（1001000KW）失效模式统计；附图4为各失效模式发生次数和观察时间为来自国内某石化企业2个站点的3个同类型压缩机整理数据，数据统计时间段为20092010年，可知Y2为2009年6月，累计日历时间为00505106小时，累计工作。
14、时间为00366106小时；根据公式（30）和附图2是由模糊综合估计法得到当前国内与国外设备技术差距的年限为Y182年，则由技术差距系数计算式（31）可得设备技术差距系数C1385。说明书CN103336903A4/4页70060以图3中的“未按照要求启动”失效模式为例，累计日历时间下观察到发生失效次数为12次，平均失效率D6968次106H，标准差SD2032次106H，C1385，带入公式（41）得到00610062由附图4中对应的“未按照要求启动”现场统计数据可知累计日历时间下观察到发生失效次数KD0，W个样本的累计观察日历时间WT00505106小时，然后由公式（53）可得后验平均失效率ED6823次106H、标准差DD1990次106H。0063同理，可得其他失效模式下的后验平均失效率核标准差，OREDA数据经贝叶斯方法调整后的失效率如附图（4）所示，通过实施例证明了基于贝叶斯理论的石化设备失效率推断方法是可行的。说明书CN103336903A1/3页8图1图2说明书附图CN103336903A2/3页9图3说明书附图CN103336903A3/3页10图4说明书附图CN103336903A10。