【技术实现步骤摘要】
本专利技术属于核电，具体涉及一种核电设备预防性维修模板时间代码的确定方法及装置。
技术介绍
1、预防性维修模板是对某一类型设备按照关键度分级和工作环境、工作频度，分别提供预防性维修工作项目(包括监督监测、定期任务、故障排查)和周期等建议的基础文件，为核电厂设备预防性维修大纲开发和优化提供指导。
2、预防性维修模板编制的一项重要的工作内容是确定时间代码。时间代码可以表示为设备部件特定故障模式预期最早故障时间。时间代码主要有3种类型，分别为：r，代表着随机，表示该故障模式在设备的运行周期中随机发生；uw，代表着无条件磨损，表示在设备的运行周期中该故障模式的发生时机与部件受到的特定压力(如温度、振动或污染)大小无关。如uw5，表示该故障模式预期的发生时间为5年；w，代表着磨损，表示在设备的运行周期中该设备故障模式的发生时间受到特定压力的影响。如，w7_10，表示如果存在加重压力，则预期的早期故障时间为7到10年之间。
3、目前国内核电厂预防性维修模板的时间代码主要参考美国电力科学研究院(electric power research institute，简称：epri)相关成果和相关人员经验来确定，人因依赖较大，因此，需要提高时间代码与国内核电厂预防性维修工作实践的匹配程度。
技术实现思路
1、为克服相关技术中存在的问题，提供了一种核电设备预防性维修模板时间代码的确定方法及装置。
2、根据本公开实施例的一方面，提供一种核电设备预防性维修模板时间代码的确
3、步骤1，针对待确定时间代码的部件类型，获取所述部件类型对应的故障数据和预防性更换数据，每个故障数据包括一个属于所述部件类型的故障部件的运行时间、运行环境、故障部位和故障模式，每个预防性更换数据包括一个属于所述部件类型的被预防性更换的部件的运行时间、运行环境和更换部位；
4、步骤2，根据所述故障数据和所述预防性更换数据确定所述备件类型的备件寿命服从的寿命分布函数；
5、步骤3，根据所述寿命分布函数的函数类型，以及各函数类型与各时间代码类型之间的对应关系，确定所述部件类型的时间代码类型，所述部件类型的时间代码类型为r、uw或w中的任意一种；
6、步骤4，在所述部件类型的时间代码类型为uw或w的情况下，根据所述寿命分布函数的参数确定变异系数；
7、步骤5，根据变异系数与变异阈值的关系确定时间代码区间。
8、在一种可能的实现方式中，步骤2包括：
9、步骤20，在故障部件的数量小于预设数值的情况下，所述备件类型的寿命分布函数为指数分布函数，并根据所述故障数据和所述预防性更换数据确定指数分布的参数；
10、步骤21，在故障部件的数量大于或等于预设数值的情况下，根据所述故障数据和所述预防性更换数据拟合得到所述备件类型的寿命分布函数。
11、在一种可能的实现方式中，步骤20包括：
12、步骤200，在故障部件的数量为0的情况下，采用下式确定指数分布函数的参数λ：
13、步骤201，在故障部件的数量为1的情况下，采用下式确定指数分布函数的参数λ：
14、其中，m为故障部件的数量，n为预防性更换部件的数量。
15、在一种可能的实现方式中，步骤21包括：
16、步骤210，根据下式确定序号为l的故障部件运行时间的多个可靠度非参估计值r(tl)：
17、rl＝rl-1+(m+n+1-rl-1)/(m+n+2-k)
18、r(tl)＝1-(rl-0.3)/(m+n+0.4)
19、其中，将全部故障部件和预防性更换部件按照运行时间由短至长的顺序排列，得到各故障部件和预防性更换部件的序号k，k∈[1,m+n]，k为整数，将全部故障部件按照运行时间由短至长的顺序排列，得到每个故障部件的序号l，l∈[1,m]，l为整数，r0＝0；
20、步骤211，确定多种分布函数类型，并针对每种分布函数类型，根据所述多个可靠度非参估计值和该分布函数类型对应的线性映射关系，确定该分布函数类型对应的数据组；
21、步骤212，对每种分布函数类型的数据组进行线性回归拟合，得到线性拟合参数和误差平方和，并根据该线性拟合参数确定该分布函数类型对应分布函数的参数；
22、步骤213，将误差平方和最小的分布函数作为寿命分布函数。
23、在一种可能的实现方式中，步骤211中：
24、在分布函数类型为指数分布时数据组(xl,yl)为xl＝tl,
25、在分布函数类型为威布尔分布时数据组(xl,yl)为xl＝ln(tl),
26、在分布函数类型为正态分布时数据组(xl,yl)为xl＝tl,yl＝φ-1(1-r(tl))；
27、在分布函数类型为对数正态分布时数据组(xl,yl)为：
28、xl＝ln(tl),yl＝φ-1(1-r(tl))；
29、其中，φ-1(x)为标准正态分布的反函数，标准正态分布函数为：
30、
31、在一种可能的实现方式中，步骤212包括：采用线性方程y＝ax+b归一拟合每种分布函数类型对应的数据组(xl,yl)，采用下式确定线性拟合参数a和b：
32、
33、
34、误差平方和
35、在分布函数类型为指数分布时参数λ＝b；
36、在分布函数类型为威布尔分布时参数β＝b，
37、在分布函数类型为正态分布时参数
38、在分布函数类型为对数正态分布时参数
39、在一种可能的实现方式中，步骤21还包括：
40、步骤214，根据下式求解似然函数l(θ)：
41、
42、在分布函数类型为指数分布时f(t)＝λe-λt，f(t)＝1-e-λt；
43、在分布函数类型为威布尔分布时
44、在分布函数类型为正态分布时
45、在分布函数类型为对数正态分布时
46、步骤215，采用牛顿-拉夫森法确定似然函数l(θ)的参数估计结果
47、步骤216，根据似然函数l(θ)的参数估计结果进行以下处理：
48、采用赤池信息准则确定不同分布函数对应的赤池信息值，并将赤池信息值最小的分布函数作为寿命分布函数；或者
49、采用贝叶斯信息准则确定不同分布函数对应的贝叶斯信息准则值，并将贝叶斯信息准则值最小的分布函数作为寿命分布函数。
50、在一种可能的实现方式中，步骤3包括：
51、步骤30，在故障数据的数量小于或等于1的情况下，所述部件类型对应的寿命分布函数为指数分布函数，所述部件类型的时间代码类型为r；
52、步骤31，所述故障数据的数量大于2的情况下，若所述部件类型对应的寿命分布函数为指数分布本文档来自技高网...
【技术保护点】
1.一种核电设备预防性维修模板时间代码的确定方法，其特征在于，所述方法包括：
2.根据权利要求1所述的方法，其特征在于，步骤2包括：
3.根据权利要求2所述的方法，其特征在于，步骤20包括：
4.根据权利要求2所述的方法，其特征在于，步骤21包括：
5.根据权利要求4所述的方法，其特征在于，步骤211中：
6.根据权利要求4所述的方法，其特征在于，步骤212包括：采用线性方程y＝ax+b归一拟合每种分布函数类型对应的数据组(xl,yl)，采用下式确定线性拟合参数a和b：
7.根据权利要求6所述的方法，其特征在于，步骤21还包括：
8.根据权利要求1所述的方法，其特征在于，步骤3包括：
9.根据权利要求8所述的方法，其特征在于，步骤3还包括：
10.根据权利要求1所述的方法，其特征在于，步骤4还包括：
11.一种核电设备预防性维修模板时间代码的确定装置，其特征在于，所述装置包括：
12.一种核电设备预防性维修模板时间代码的确定装置，其特征在于，所述装置
13.一种非易失性计算机可读存储介质，其上存储有计算机程序指令，其特征在于，所述计算机程序指令被处理器执行时实现权利要求1至10中任意一项所述的方法。
...【技术特征摘要】
1.一种核电设备预防性维修模板时间代码的确定方法，其特征在于，所述方法包括：
2.根据权利要求1所述的方法，其特征在于，步骤2包括：
3.根据权利要求2所述的方法，其特征在于，步骤20包括：
4.根据权利要求2所述的方法，其特征在于，步骤21包括：
5.根据权利要求4所述的方法，其特征在于，步骤211中：
6.根据权利要求4所述的方法，其特征在于，步骤212包括：采用线性方程y＝ax+b归一拟合每种分布函数类型对应的数据组(xl,yl)，采用下式确定线性拟合参数a和b：
7.根据权利要求6所述的方法，其特...
【专利技术属性】
技术研发人员：麻浩军，徐中，刘东林，李志涛，蔡胜武，刘伟，朱昌荣，聂成平，
申请(专利权)人：中核核电运行管理有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。