CN116205626A - 核电厂非规范项目的在役检查计划制定方法和系统 - Google Patents
核电厂非规范项目的在役检查计划制定方法和系统 Download PDFInfo
- Publication number
- CN116205626A CN116205626A CN202310162520.1A CN202310162520A CN116205626A CN 116205626 A CN116205626 A CN 116205626A CN 202310162520 A CN202310162520 A CN 202310162520A CN 116205626 A CN116205626 A CN 116205626A
- Authority
- CN
- China
- Prior art keywords
- failure probability
- distribution
- equipment failure
- inspection
- service inspection
- 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.)
- Pending
Links
- 238000007689 inspection Methods 0.000 title claims abstract description 146
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000009826 distribution Methods 0.000 claims abstract description 132
- 238000012423 maintenance Methods 0.000 claims abstract description 20
- 230000006870 function Effects 0.000 claims description 59
- 238000005070 sampling Methods 0.000 claims description 22
- 238000004590 computer program Methods 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- 230000007547 defect Effects 0.000 abstract description 13
- 238000002360 preparation method Methods 0.000 abstract description 2
- 230000009471 action Effects 0.000 description 26
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000003466 welding Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000032683 aging Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 229920000609 methyl cellulose Polymers 0.000 description 2
- 239000001923 methylcellulose Substances 0.000 description 2
- 235000010981 methylcellulose Nutrition 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013477 bayesian statistics method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000000246 remedial effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000009958 sewing Methods 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Primary Health Care (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Probability & Statistics with Applications (AREA)
- Public Health (AREA)
- Algebra (AREA)
- Quality & Reliability (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Entrepreneurship & Innovation (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明涉及核电厂非规范项目的在役检查计划制定方法和系统，包括：根据检查成本、维修成本以及失效损失成本，构建损失函数；根据损失函数，计算设备失效概率的临界值；构建非规范项目的设备失效概率的威布尔分布；根据设备失效概率的威布尔分布和设备失效概率的临界值进行计算，获得在役检查计划时间的临界值；根据在役检查计划时间的临界值制定在役检查计划。本发明根据在役检查计划时间的临界值制定在役检查计划，实现在役检查计划的定量制定，从而实现了在失效概率非常低的阶段，不安排在役检查；在失效概率较低的阶段，制定抽检计划，避免过度检修；在失效概率较高的阶段，安排全检，避免设备失效前未能及时检查发现出现缺陷的风险。
Description
技术领域
本发明涉及核电厂在役检查的技术领域，更具体地说，涉及一种核电厂非规范项目的在役检查计划制定方法和系统。
背景技术
在核电厂运行阶段，机械部件可能受到诸如应力、温度、辐照、氢吸附、腐蚀、振动和磨损等多种因素，从而导致机械部件的材料性能发生变化，例如老化、脆化、疲劳以及缺陷的形成和扩展。这些缺陷是指机械部件中的不连续、不规则和损伤等，例如裂纹、磨损减薄等。因此，有必要在核电厂机械部件服役期间定期或不定期地进行检查，以发现可能存在的缺陷；或者跳跃缺陷的扩展情况，以判断这些缺陷是否对核电厂继续安全运行是否可接受，或是否有必要采取维修或更换的补救措施。
通常核电厂的在役检查项目主要是按照在役检查相关标准或规范执行。核电厂参考这些在役检查标准或规范，制定在役检查计划，并开展执行在役检查工作。但是，随着核电厂的运行，会出现不在当前标准或规范指导的在役检查范围内的机械部件会出现各缺陷，并导致失效的各类内部或外部经验反馈，针对这些经验反馈如果要进行在役检查，这些项目就属于非规范在役检查项目。
目前对于非规范项目在制定在役检查计划时，由于无规范标准可以直接参考使用，通常是定性分析层面做出在役检查计划。这类定性分析方法，通常是将经验反馈事件中失效部件与本核电厂同类部件进行多维度对比，比如设计、制造、运行环境等因素进行比较。如果存在较多相似的因素，通常从设备安全保守角度出发，会对该类设备执行在役检查，但是没有在役检查计划的定量方法，因此，通常制定的在役检查计划是固定周期，或者连续跟踪几轮检查，与设备失效概率没有较好的匹配性。会存在如下不足：
如果经验反馈事件中设备失效概率非常低，仅是个案，采用固定周期检查或者连续跟踪几轮检查，虽然更加保守，但是在设备服役早期，设备的可靠性较高大部分检查结果都是未见异常，会导致过度检修。
如果经验反馈事件中设备失效概率较高，存在普遍情况，采用固定周期检查或者连续跟踪几轮检查，虽然能降低检查成本，但是在设备服役老化期，可能存在设备失效前未能及时检查发现出缺陷的风险。
发明内容
本发明要解决的技术问题在于，提供一种核电厂非规范项目的在役检查计划制定方法和系统。
本发明解决其技术问题所采用的技术方案是：构造一种核电厂非规范项目的在役检查计划制定方法，包括以下步骤：
根据检查成本、维修成本以及失效损失成本，构建损失函数；
根据所述损失函数，计算设备失效概率的临界值；
构建非规范项目的设备失效概率的威布尔分布；
根据所述设备失效概率的威布尔分布和所述设备失效概率的临界值进行计算，获得在役检查计划时间的临界值；
根据所述在役检查计划时间的临界值制定在役检查计划。
在本发明所述的核电厂非规范项目的在役检查计划制定方法中，所述根据检查成本、维修成本以及失效损失成本，构建损失函数包括：
根据所述检查成本、所述维修成本以及所述失效损失成本，建立支付函数；
根据贝叶斯决策和所述支付函数进行转化，获得所述损失函数。
在本发明所述的核电厂非规范项目的在役检查计划制定方法中，所述根据所述损失函数，计算设备失效概率的临界值包括：
对所述损失函数进行分析，确定抽检损失函数；
根据设定损失值和所述抽检损失函数进行计算，获得所述设备失效概率的临界值。
在本发明所述的核电厂非规范项目的在役检查计划制定方法中，所述构建非规范项目的设备失效概率的威布尔分布包括：
采用贝叶期统计法，确定设备失效概率的先验分布；
根据不同时刻的样本信息和所述设备失效概率的先验分布，确定设备失效概率的后验分布；所述不同时刻的样本信息为不同时刻的部件抽检数量和部件失效数量；
根据所述设备失效概率的后验分布，构建所述设备失效概率的威布尔分布。
在本发明所述的核电厂非规范项目的在役检查计划制定方法中，所述根据不同时刻的样本信息和所述设备失效概率的先验分布，确定设备失效概率的后验分布包括：
获取不同时刻的样本信息；
根据所述设备失效概率的先验分布和贝叶斯公式，获得不同时刻的后验分布。
在本发明所述的核电厂非规范项目的在役检查计划制定方法中，所述根据所述设备失效概率的后验分布，确定所述设备失效概率的威布尔分布包括：
对所述不同时刻的后验分布进行均值处理，获得不同时刻的后验分布的失效概率均值；
根据所述后验分布的失效概率均值、采用配分布曲线法和威布尔分布法进行计算，获得所述设备失效概率的威布尔分布。
在本发明所述的核电厂非规范项目的在役检查计划制定方法中，所述根据所述后验分布的失效概率均值、采用配分布曲线法和威布尔分布法进行计算，获得所述设备失效概率的威布尔分布包括：
采用配分布曲线法和威布尔分布，确定不同时刻的失效概率；
基于所述不同时刻的后验分布的失效概率均值和所述不同时刻的失效概率进行计算，获得威布尔分布的参数；
根据所述威布尔分布的参数获得所述设备失效概率的威布尔分布。
本发明还提供一种核电厂非规范项目的在役检查计划制定系统，包括：
第一构建模块，用于根据检查成本、维修成本以及失效损失成本，构建损失函数；
第一计算模块，用于根据所述损失函数，计算设备失效概率的临界值；
第二构建模块，用于构建非规范项目的设备失效概率的威布尔分布；
第二计算模块，用于根据所述设备失效概率的威布尔分布和所述设备失效概率的临界值进行计算，获得在役检查计划时间的临界值；
检查计算制定模块，用于根据所述在役检查计划时间的临界值制定在役检查计划。
本发明还提供一种存储介质，所述存储介质存储有计算机程序，所述计算机程序适于处理器进行加载，以执行如上所述的核电厂非规范项目的在役检查计划制定方法的步骤。
本发明还提供一种电子设备，包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器通过调用所述存储器中存储的所述计算机程序，执行如上所述的核电厂非规范项目的在役检查计划制定方法的步骤。
实施本发明的核电厂非规范项目的在役检查计划制定方法和系统，具有以下有益效果：包括以下步骤：根据检查成本、维修成本以及失效损失成本，构建损失函数；根据损失函数，计算设备失效概率的临界值；构建非规范项目的设备失效概率的威布尔分布；根据设备失效概率的威布尔分布和设备失效概率的临界值进行计算，获得在役检查计划时间的临界值；根据在役检查计划时间的临界值制定在役检查计划。本发明根据在役检查计划时间的临界值制定在役检查计划，实现在役检查计划的定量制定，从而实现了在失效概率非常低的阶段，不安排在役检查；在失效概率较低的阶段，制定抽检计划，避免过度检修；在失效概率较高的阶段，安排全检，避免设备失效前未能及时检查发现出现缺陷的风险。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明实施例提供的核电厂非规范项目的在役检查计划制定方法的流程示意图；
图2是本发明实施例提供的各行动项支付函数的示意图；
图3是本发明实施例提供的损失函数的示意图；
图4是本发明实施例提供的在设定损失值下的损失函数的示意图；
图5是本发明实施例提供的设备失效概率的威布尔分布的示意图；
图6是本发明提供的一个具体实施例的设备失效概率的威布尔分布的示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参考图1，本发明提供了一种核电厂非规范项目的在役检查计划制定方法，该在役检查计划制定方法可以实现在役检查计划的定量制定。
具体的，如图1所示，该核电厂非规范项目的在役检查计划制定方法包括以下步骤：
步骤S10、根据检查成本、维修成本以及失效损失成本，构建损失函数。
可选的，本发明实施例中，根据检查成本、维修成本以及失效损失成本，构建损失函数包括：根据检查成本、维修成本以及失效损失成本，建立支付函数；根据贝叶斯决策和支付函数进行转化，获得损失函数。
具体的，首先进行参数定义：
设备失效概率可记为θ，其中，0＜θ＜1；被检部件总数为M，抽检数为n，其中，0＜n＜M。
每件设备或者部件投入的检查成本为CJ，失效设备/部件的维修成本为CY，未检查的设备/部件投入运行后出现失效后造成的全部损失(包括维修成本)，即失效损失成本为CP。一般地，CY＜CP。
接着，建立行动项：
对于采取的在役检查行动按照检查数量可以分为3类行动，具体如下：
行动项a1：执行全检，检查数即为M。
行动项a2：执行抽检，检查数即为n。
行动项a3：不做检查，检查数即为0。
然后，建立支付函数：
在支付函数建立时，未检查的设备投入运行后出现失效后造成的全部损失也计入支付。其中，各行动对应的支付函数如下：
行动项a1的支付函数(L(θ,a1))：MCJ+MθCY。
行动项a2的支付函数(L(θ,a2))：nCJ+nθCY+(M-n)θCP。
行动项a3的支付函数(L(θ,a3))：MθCP。
其中，以θ作为自变量，支付函数W为因变量，通过计算可以得到：
在θ＝CJ/(CP-CY)时，三个行动项各支付数值相同。令CJ/(CP-CY)为A，A即为行动决策临界点，定义A为检损比。其中，检损比表示检修支出与失效损失支出的比值，支付函数如图2所示。
最后，建立损失函数：
具体的，按照贝叶斯决策，根据支付函数进行转化，获得损失函数。其中，各行动项的损失函数如下：
其中，各行动项的损失函数示意图如图3所示。
步骤S20、根据损失函数，计算设备失效概率的临界值。
可选的，本发明实施例中，根据损失函数，计算设备失效概率的临界值包括：对损失函数进行分析，确定抽检损失函数；根据设定损失值和抽检损失函数进行计算，获得设备失效概率的临界值。
具体的，由图3中的各行动项的损失函数的示意图可以看出，在不考虑可接受损失值的条件下，最优决策行动为：当0<θ<A时，最优行动为a3；当A<θ<1时，最优行动为a1；当θ＝A时，三个行动均可以作为最优行动。
然而，由于最优行动不考虑损失值L0，在0<θ<A时，最优行动为a3，当A<θ<1时，最优行动为a1。随着失效概率θ的增加从不检直接到全检，虽然整体期望损失最小，但是采取这种决策，会在行动决策临界点前出现完全不做检查，从设备安全角度又不够保守，因为在A值附近，执行抽样检查相比不做检查的支付损失较小。而在行动决策临界点后又执行全检，从经济角度又存在过度检修的情况，因为在A值附近，即便执行抽样检查，相比全检的损失也较小。
本发明通过考虑可接受损失值L0的最优行动决策：即允许承受一定的成本损失下，采用抽样检查，以避免完全不做检查而造成安全方面不够保守，或者执行全检而造成的过度检修的两个极端。基于此，本发明根据抽检损失函数，在设定损失值(即可接受损失值L0)的条件下，计算设备失效概率的临界值。
其中，本发明实施例中，设备失效概率的临界值包括：第一临界值(定义为θ1)和第二临界值(定义为θ2)。具体的，令：
L(θ,a2)＝0 (1)。
根据(1)式进行计算可以得到：
由0<θ1≤θ2<1，当L0＝0时，θ1＝θ2＝A，因此，可以得到：
其中，可接受损失值L0的损失函数如图4所示。
由图4中可以看出：
在(0,θ1)区间内，部件失效概率非常低，该阶段设备失效导致损失成本期望远低于可接受损失值，且投入检查(无论是抽检还是全检)将导致检修成本超过可接受损失值，可不进行在役检查。
在[θ1,θ2]区间内，部件失效概率较低，该阶段设备失效导致损失成本期望或者投入抽样检修成本均可控制在可接受损失值范围内，可选择抽样在役检查。
在(θ2,1)区间内，部件失效概率较高，该阶段不做检查或执行抽样检查，而设备失效导致损失成本期望将高于可接受损失值，需进行全检，以便在设备失效前发现缺陷以便及时处理。
步骤S30、构建非规范项目的设备失效概率的威布尔分布。
可选的，本发明实施例中，构建非规范项目的设备失效概率的威布尔分布包括：采用贝叶期统计法，确定设备失效概率的先验分布；根据不同时刻的样本信息和设备失效概率的先验分布，确定设备失效概率的后验分布；不同时刻的样本信息为不同时刻的部件抽检数量和部件失效数量；根据设备失效概率的后验分布，构建设备失效概率的威布尔分布。
具体的，根据不同时刻的样本信息和设备失效概率的先验分布，确定设备失效概率的后验分布包括：获取不同时刻的样本信息；根据设备失效概率的先验分布和贝叶斯公式，获得不同时刻的后验分布。
其中，根据设备失效概率的后验分布，确定设备失效概率的威布尔分布包括：对不同时刻的后验分布进行均值处理，获得不同时刻的后验分布的失效概率均值；根据后验分布的失效概率均值、采用配分布曲线法和威布尔分布法进行计算，获得设备失效概率的威布尔分布。
进一步地，根据后验分布的失效概率均值、采用配分布曲线法和威布尔分布法进行计算，获得设备失效概率的威布尔分布包括：采用配分布曲线法和威布尔分布，确定不同时刻的失效概率；基于不同时刻的后验分布的失效概率均值和不同时刻的失效概率进行计算，获得威布尔分布的参数；根据威布尔分布的参数获得设备失效概率的威布尔分布。
具体的，首先：
采用贝叶斯统计方法，选择贝塔分布作为设备失效概率的先验分布，即：
(4)式中，
其中，a、b为先验分布的超参数，x为自变量。本发明实施例中，a、b的值可以根据贝叶斯统计学和经验反馈数据确定。或者，当没有任何经验反馈数据时，可以采用贝叶斯假设得到。即取a＝1，b＝1，此时，设备失效概率的先验分布为(0，1)均匀分布。
其次：
在得到设备失效概率的先验分布后，根据不同时刻的样本信息，再结合贝叶斯公式，即可得到不同时刻的后验分布。即：
(5)式中，Ni为抽检数量，Xi为失效数量，即Ni和Xi为不同时刻Ti的样本信息，其中，Ti为设备的服役时刻，i表示不同时刻的序号，即i＝1,2,..,k，共计k次样本信息。
然后：
在得到不同时刻的后验分布后，对不同时刻的后验分布进行均值处理，得到不同时刻的后验分布的失效概率均值，即：
(6)式中，
为第i时刻的后验分布的失效概率均值。
接着，根据配分布曲线法和威布尔分布法进行计算，可以获得设备失效概率的威布尔分布。即根据配分布曲线法和威布尔分布可以得到：
(7)式中，θi为设备在Ti时刻的失效概率，P表示概率，T为设备的服役累积时间，m为在Ti时刻的威布尔形状参数，η为在Ti时刻的威布尔尺寸参数。
对(7)式取两次对数后线性转换，再将(6)式代入后，采用最小二乘法可以计算得到威布尔分布的参数。其中，这里的威布尔分布的参数为拟合参数，即采用最小二乘法计算得到威布尔分布的拟合形状参数(表示为
)和威布尔分布的拟合尺寸参数(表示为/>
)。
最后，将
和/>
代入威布尔分布即可得到：
步骤S40、根据设备失效概率的威布尔分布和设备失效概率的临界值进行计算，获得在役检查计划时间的临界值。
具体的，本发明实施例中，在步骤S30中得到设备失效概率的威布尔分布后(即公式(8))，将(2)式和(3)式分别代入公式(8)即可得到在可接受损失值下的在役检查计划时间的临界值。其中，在役检查计划时间的临界值包括：第一临界时间值(表示为T1)和第二临界时间值(表示为T2)。
其中，T1和T2具体为：
或者，在其他一些实施例中，在获得设备失效概率的威尔布分布后，可以通过建模的方式快速得到第一临界时间值和第二临界时间值。其中，建模函数如图5所示。
步骤S50、根据在役检查计划时间的临界值制定在役检查计划。
具体的，本发明实施例中，在获得在役检查计划时间的第一临界时间值和第二临界时间值后，即可制定相应的在役检查计划，具体如下：
在(0，T1)区间内，设备失效概率非常低，该阶段设备失效导致损失成本期望远低于可接受损失值，且投入检查将导致检修成本超过可接受损失值，可不进行在役检查。
在[T1，T2]区间内，设备失效概率较低，该阶段设备失效导致损失成本期望或者投入抽样检修成本均可控制在可接受损失值范围内，可选择抽样在役检查。
在(T2，T寿期末)区间内，设备失效概率较高，该阶段设备失效导致损失成本期望将高于可接受损失值，需进行全检，以便在设备失效前发现缺陷以便及时处理。
通过使用本发明提供的核电厂非规范项目的在役检查计划制定方法，可以实现定量确定在役检查抽检计划，可以在失效概率非常低的阶段避免过度检修，可以在失效概率较高的阶段避免设备失效前未能及时检查发现出缺陷的风险。
下面以一个具体的实施例进行说明。
基于本发明对某核电厂某系统焊缝制定在投检查计划，具体过程如下：
步骤1.初始化参数：
某系统焊缝总数为M，M＝6；抽检焊缝数为n，n＝2；可接受损失值为L0，L0＝5000；每条焊缝投入的检查成本为CJ，CJ＝10000；通过检查发现失效焊缝进行维修成本为CY，CY＝50000；未检查的焊缝投入运行后出现失效后造成的全部损失(包括维修成本)为CP，CP＝100000；该焊缝预计最大寿期为T寿期末，T寿期末＝40。
步骤2.计算检损比A，A＝CJ/(CP-CY)＝0.2。
步骤3.在给定的可接受损失值L0下，计算设备失效概率的第一临界值θ1和第二临界值θ2，即：
步骤4.考虑到该焊缝失效情况无任何先验信息，采用贝叶斯假设，即取a＝1，b＝1，此时失效概率的先验分布为(0，1)均匀分布。经过6轮次检查的样本信息，可计算出后验分布下失效概率均值如下表1所示：
表1：检查样本信息及后验分布下失效概率均值
步骤5.根据表1的信息和威布尔分布，采用最小二乘法计算可得威布尔分布的参数，即：
进行线性显著性检验，在95％置信度下，即α＝0.05，检验值t＝3.956＞tα/2＝2.776，符合线性显著性，说明该焊缝寿命服从威布尔分布，失效概率为威布尔累积分布值。
步骤6.计算在可接受损失值L0下的T1和T2两个临界时间点，即：
T1＝9.512；T2＝14.859。
由于检查轮次Ti为整数，且从设备安全保守角度，T1和T2向下取整，取整后，T1＝9，T2＝14。
步骤7.对设备失效概率的威布尔分布进行建模，如图6所示。
步骤8.通过计算出在可接受损失值L0下的T1和T2两个临界时间点，可以在定量分析后制定在役检查计划，具体如下：
在(0，9)轮次区间内，设备失效概率非常低，该阶段设备失效导致损失成本期望远低于可接受损失值，且投入检查将导致检修成本超过可接受损失值，可不进行在役检查。
在[9，14]轮次区间内，设备失效概率较低，该阶段设备失效导致损失成本期望或者投入抽样检修成本均可控制在可接受损失值范围内，可选择抽样在役检查。
在(14，40)轮次区间内，设备失效概率较高，该阶段设备失效导致损失成本期望将高于可接受损失值，需进行在役检查，以便在设备失效前发现缺陷以便及时处理。
本发明还提供一种核电厂非规范项目的在役检查计划制定系统，具体包括：
第一构建模块，用于根据检查成本、维修成本以及失效损失成本，构建损失函数。
第一计算模块，用于根据所述损失函数，计算设备失效概率的临界值。
第二构建模块，用于构建非规范项目的设备失效概率的威布尔分布。
第二计算模块，用于根据所述设备失效概率的威布尔分布和所述设备失效概率的临界值进行计算，获得在役检查计划时间的临界值。
检查计算制定模块，用于根据所述在役检查计划时间的临界值制定在役检查计划。
具体的，这里的核电厂非规范项目的在役检查计划制定系统中各单元之间具体的配合操作过程具体可以参照上述核电厂非规范项目的在役检查计划制定方法，这里不再赘述。
另，本发明的一种电子设备，包括存储器和处理器；存储器用于存储计算机程序；处理器用于执行计算机程序实现如上面任意一项的核电厂非规范项目的在役检查计划制定方法。具体的，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括承载在计算机可读介质上的计算机程序，该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过电子设备下载和安装并且执行时，执行本发明实施例的方法中限定的上述功能。本发明中的电子设备可为笔记本、台式机、平板电脑、智能手机等终端，也可为服务器。
另，本发明的一种存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上面任意一项的核电厂非规范项目的在役检查计划制定方法。具体的，需要说明的是，本发明上述的存储介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是但不限于是电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本发明中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本发明中，计算机可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读信号介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：电线、光缆、RF(射频)等等，或者上述的任意合适的组合。
上述计算机可读介质可以是上述电子设备中所包含的；也可以是单独存在，而未装配入该电子设备中。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人士能够了解本发明的内容并据此实施，并不能限制本发明的保护范围。凡跟本发明权利要求范围所做的均等变化与修饰，均应属于本发明权利要求的涵盖范围。
Claims (10)
1.一种核电厂非规范项目的在役检查计划制定方法，其特征在于，包括以下步骤：
根据检查成本、维修成本以及失效损失成本，构建损失函数；
根据所述损失函数，计算设备失效概率的临界值；
构建非规范项目的设备失效概率的威布尔分布；
根据所述设备失效概率的威布尔分布和所述设备失效概率的临界值进行计算，获得在役检查计划时间的临界值；
根据所述在役检查计划时间的临界值制定在役检查计划。
2.根据权利要求1所述的核电厂非规范项目的在役检查计划制定方法，其特征在于，所述根据检查成本、维修成本以及失效损失成本，构建损失函数包括：
根据所述检查成本、所述维修成本以及所述失效损失成本，建立支付函数；
根据贝叶斯决策和所述支付函数进行转化，获得所述损失函数。
3.根据权利要求1所述的核电厂非规范项目的在役检查计划制定方法，其特征在于，所述根据所述损失函数，计算设备失效概率的临界值包括：
对所述损失函数进行分析，确定抽检损失函数；
根据设定损失值和所述抽检损失函数进行计算，获得所述设备失效概率的临界值。
4.根据权利要求1所述的核电厂非规范项目的在役检查计划制定方法，其特征在于，所述构建非规范项目的设备失效概率的威布尔分布包括：
采用贝叶期统计法，确定设备失效概率的先验分布；
根据不同时刻的样本信息和所述设备失效概率的先验分布，确定设备失效概率的后验分布；所述不同时刻的样本信息为不同时刻的部件抽检数量和部件失效数量；
根据所述设备失效概率的后验分布，确定所述设备失效概率的威布尔分布。
5.根据权利要求4所述的核电厂非规范项目的在役检查计划制定方法，其特征在于，所述根据不同时刻的样本信息和所述设备失效概率的先验分布，确定设备失效概率的后验分布包括：
获取不同时刻的样本信息；
根据所述设备失效概率的先验分布和贝叶斯公式，获得不同时刻的后验分布。
6.根据权利要求5所述的核电厂非规范项目的在役检查计划制定方法，其特征在于，所述根据所述设备失效概率的后验分布，确定所述设备失效概率的威布尔分布包括：
对所述不同时刻的后验分布进行均值处理，获得不同时刻的后验分布的失效概率均值；
根据所述后验分布的失效概率均值、采用配分布曲线法和威布尔分布法进行计算，获得所述设备失效概率的威布尔分布。
7.根据权利要求1所述的核电厂非规范项目的在役检查计划制定方法，其特征在于，所述根据所述后验分布的失效概率均值、采用配分布曲线法和威布尔分布法进行计算，获得所述设备失效概率的威布尔分布包括：
采用配分布曲线法和威布尔分布，确定不同时刻的失效概率；
基于所述不同时刻的后验分布的失效概率均值和所述不同时刻的失效概率进行计算，获得威布尔分布的参数；
根据所述威布尔分布的参数获得所述设备失效概率的威布尔分布。
8.一种核电厂非规范项目的在役检查计划制定系统，其特征在于，包括：
第一构建模块，用于根据检查成本、维修成本以及失效损失成本，构建损失函数；
第一计算模块，用于根据所述损失函数，计算设备失效概率的临界值；
第二构建模块，用于构建非规范项目的设备失效概率的威布尔分布；
第二计算模块，用于根据所述设备失效概率的威布尔分布和所述设备失效概率的临界值进行计算，获得在役检查计划时间的临界值；
检查计算制定模块，用于根据所述在役检查计划时间的临界值制定在役检查计划。
9.一种存储介质，其特征在于，所述存储介质存储有计算机程序，所述计算机程序适于处理器进行加载，以执行如权利要求1至7任一项所述的核电厂非规范项目的在役检查计划制定方法的步骤。
10.一种电子设备，其特征在于，包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器通过调用所述存储器中存储的所述计算机程序，执行如权利要求1至7任一项所述的核电厂非规范项目的在役检查计划制定方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310162520.1A CN116205626A (zh) | 2023-02-22 | 2023-02-22 | 核电厂非规范项目的在役检查计划制定方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310162520.1A CN116205626A (zh) | 2023-02-22 | 2023-02-22 | 核电厂非规范项目的在役检查计划制定方法和系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116205626A true CN116205626A (zh) | 2023-06-02 |
Family
ID=86518810
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310162520.1A Pending CN116205626A (zh) | 2023-02-22 | 2023-02-22 | 核电厂非规范项目的在役检查计划制定方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116205626A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118013699A (zh) * | 2024-01-03 | 2024-05-10 | 阳江核电有限公司 | 蒸汽发生器在役检查抽样方案的优化验证方法和装置 |
-
2023
- 2023-02-22 CN CN202310162520.1A patent/CN116205626A/zh active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118013699A (zh) * | 2024-01-03 | 2024-05-10 | 阳江核电有限公司 | 蒸汽发生器在役检查抽样方案的优化验证方法和装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lee et al. | Reducing mold changing time by implementing Lean Six Sigma | |
| WO2014107303A1 (en) | Probabilistic modeling and sizing of embedded flaws in nondestructive inspections for fatigue damage prognostics and structural integrity assessment | |
| CN116205626A (zh) | 核电厂非规范项目的在役检查计划制定方法和系统 | |
| CN104335056A (zh) | 测试器和材料处理设备之间用于分开和控制测试单元操作中多个实体的不同请求的插入器 | |
| Deepak Prabhakar et al. | A new model for reliability centered maintenance in petroleum refineries | |
| JP2009098093A (ja) | 設備の有効メンテナンスモニタ装置 | |
| CN104364664A (zh) | 用于创建、定义和执行spc规则决策树的算法和结构 | |
| CN115147236A (zh) | 一种处理方法、处理装置和电子设备 | |
| CN115718450A (zh) | 一种设备停线监控方法、装置、电子设备及系统 | |
| CN116720278A (zh) | 一种复合加载疲劳寿命预测方法、装置及可读存储介质 | |
| CN113554324B (zh) | 一种生产过程撤回方法、系统、设备及存储介质 | |
| CN116680196A (zh) | 一种基于人工智能的测试管理方法、装置及电子设备 | |
| CN113987031A (zh) | 一种电力大数据的数据异常判别方法 | |
| CN109039714A (zh) | 云计算系统中资源的管理方法和装置 | |
| CN110993132B (zh) | 一种核电厂支持疲劳监测功能的瞬态监测方法 | |
| CN112183555A (zh) | 焊接质量的检测方法、系统、电子设备和存储介质 | |
| CN117170949A (zh) | 测试方法、装置、电子设备及可读存储介质 | |
| CN110212585A (zh) | 基于统计分析的风电机组运行可靠性预测方法、装置及风电机组 | |
| CN115841082A (zh) | 一种加气站异常诊断系统及方法 | |
| CN113341366B (zh) | 一种用户电表状态监测方法、装置设备和存储介质 | |
| CN116542036B (zh) | 核电厂在役检查实施间隔的计算方法和装置 | |
| CN113591909A (zh) | 电力系统的异常检测方法、异常检测装置以及存储介质 | |
| CN111260133A (zh) | 用于监测研发质量的质量监测模型的优化方法及装置 | |
| JP2022072653A (ja) | 保全支援方法および保全支援装置 | |
| KR101866458B1 (ko) | 공사 마무리 및 시운전 관리 방법, 및 시스템 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |