CN116258221A - Nuclear power plant spare part demand prediction method based on normal distribution - Google Patents

Nuclear power plant spare part demand prediction method based on normal distribution Download PDF

Info

Publication number
CN116258221A
CN116258221A CN202111463258.1A CN202111463258A CN116258221A CN 116258221 A CN116258221 A CN 116258221A CN 202111463258 A CN202111463258 A CN 202111463258A CN 116258221 A CN116258221 A CN 116258221A
Authority
CN
China
Prior art keywords
normal distribution
spare part
power plant
nuclear power
method based
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
Application number
CN202111463258.1A
Other languages
Chinese (zh)
Inventor
姚昊
葛旭阳
吴宝华
蔡胜武
皮敏
胡文勇
熊伟
李志鹏
杜君尧
刘忠良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CNNC Nuclear Power Operation Management Co Ltd
Original Assignee
CNNC Nuclear Power Operation Management Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CNNC Nuclear Power Operation Management Co Ltd filed Critical CNNC Nuclear Power Operation Management Co Ltd
Priority to CN202111463258.1A priority Critical patent/CN116258221A/en
Publication of CN116258221A publication Critical patent/CN116258221A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/08Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/087Inventory or stock management, e.g. order filling, procurement or balancing against orders
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Software Systems (AREA)
  • Educational Administration (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Accounting & Taxation (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Algebra (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Finance (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention belongs to the technical field of spare part management, and particularly relates to a nuclear power plant spare part demand prediction method based on normal distribution. The method comprises the following steps: step 1: acquiring a parameter mean mu and a standard deviation sigma of normal distribution according to spare part life data; step 2: acquiring expected values of failure times of spare parts in a given time interval according to normal distribution; step 3: inventory quota for the spare parts is determined based on the service level of the spare parts. The invention has the beneficial effects that: at present, the inventory quota of spare parts is determined by the nuclear power plant manually according to experience, the subjectivity is strong, the quota is conservative, and the quantitative calculation of the demand and the probability of the spare parts with the life distribution obeying normal distribution in a given time interval in the future can be realized by the method provided by the invention, so that the manual subjective judgment is reduced, and the inventory of the spare parts is reduced.

Description

Nuclear power plant spare part demand prediction method based on normal distribution
Technical Field
The invention belongs to the technical field of spare part management, and particularly relates to a nuclear power plant spare part demand prediction method based on normal distribution.
Background
In general, due to technical deficiency, economic limitations, etc., it is impossible to design a product to fully fulfill its intended function throughout its life cycle, which may lead to downtime for commercial equipment (e.g., nuclear power plants, airplanes, high-speed rails, etc.), at which point the assurance of spare parts is important. When the components are expensive, the inventory of spare parts must be properly managed, as a low inventory means an increased likelihood of waiting for spare parts, and a high inventory means too much money is spent. To ensure a certain safety stock to meet the demand for unplanned replacement of spare parts in field service work, nuclear power plants implement spare part quota management.
Spare part demand is an important input to spare part quota management, and its prediction accuracy is of great importance to reduce inventory and ensure on-site operation. The spare part demand prediction methods generally adopted mainly have two types: the first is a reliability-based method, and the second is a black box method based on spare part consumption history data. In some cases, spare part requirements present a pattern that is not well predicted by conventional methods.
Disclosure of Invention
The invention aims to provide a nuclear power plant spare part demand prediction method based on normal distribution, which can ensure spare part consumption requirements of a nuclear power plant in a certain time, rationalize spare part inventory and provide support for better developing the rated management work of the spare parts of the nuclear power plant.
The technical scheme of the invention is as follows: a nuclear power plant spare part demand prediction method based on normal distribution comprises the following steps:
step 1: acquiring a parameter mean mu and a standard deviation sigma of normal distribution according to spare part life data;
step 2: acquiring expected values of failure times of spare parts in a given time interval according to normal distribution;
step 3: inventory quota for the spare parts is determined based on the service level of the spare parts.
Step 1 fits the life data of spare parts with life obeying normal distribution to normal distribution according to reliability theory, and the specific process is as follows:
step 11: for all complete data t i Using functions
Figure BDA0003390221940000021
Calculate, recorded as LK i The method comprises the steps of carrying out a first treatment on the surface of the For the truncated data t j Use +.>
Figure BDA0003390221940000022
Calculate, recorded as LK j Wherein->
Figure BDA0003390221940000023
Step 12: all LK is taken i And LK (sum of LK) j Summing to obtain LK;
step 13: solving the mean value estimated value when LK takes the maximum value by using Excel programming solving function, matlab fsolve function and other tools
Figure BDA0003390221940000024
And standard deviation estimate +.>
Figure BDA0003390221940000025
Step 2 calculates expected failure times in a given interval (0, t) according to the normal distribution obtained in step 1, wherein the calculation formula is as follows:
Figure BDA0003390221940000026
and the step 2 is used for calculating M (t), and comprises the following steps:
step 21: dividing the interval (0, t) into N equal parts, wherein the greater the interval length deltat, namely t=N×deltat, the greater the N, the higher the calculation accuracy of M (t);
step 22: calculating an expected value of the average number of failures
Figure BDA0003390221940000027
Wherein F (t) is a normal distributed cumulative probability density function; t is t i For the position of the ith part Deltat in the interval (0, t), t i =i×Δt;
Step 23: calculating variance
Figure BDA0003390221940000031
Wherein: var [ N [ t ] ] is the variance of the number of failures that occur in the spare part during the time interval (0, t).
The step 3 includes the steps of,
step 31: assuming that there are S positions requiring the use of a spare part, the life of each spare part is L at the time of prediction i Then the average demand for all the positional spare parts after the lapse of time L is
Figure BDA0003390221940000032
Variance is
Figure BDA0003390221940000033
Step 32: calculating inventory quota D using poisson distribution p =P -1 (k%,M s ) Wherein P is -1 () An inverse function representing the Poisson's distribution cumulative density function, k being the service level to be achieved by the spare part, M s Calculating inventory quota D for parameters of poisson distribution using normal distribution N =N -1 (k%,M s ,var[N s (t)]) Wherein N is -1 () An inverse function representing a normal distribution cumulative density function, k being the service level to be achieved by the spare part, M s Is the mean value of normal distribution, var [ N ] s (t)]Is the variance of the normal distribution.
The invention has the beneficial effects that: at present, the inventory quota of spare parts is determined by the nuclear power plant manually according to experience, the subjectivity is strong, the quota is conservative, and the quantitative calculation of the demand and the probability of the spare parts with the life distribution obeying normal distribution in a given time interval in the future can be realized by the method provided by the invention, so that the manual subjective judgment is reduced, and the inventory of the spare parts is reduced.
Drawings
Fig. 1 is a schematic flow chart of a method for predicting the demand of spare parts of a nuclear power plant based on normal distribution.
Detailed Description
The invention will be described in further detail with reference to the accompanying drawings and specific examples.
The invention is suitable for demand prediction of spare parts of the nuclear power plant, such as vehicle tyre wear, transformers, alarm bulbs and the like, with the service life obeying normal distribution.
As shown in fig. 1, a nuclear power plant spare part demand prediction method based on normal distribution includes the following steps:
step 1: acquiring a parameter mean mu and a standard deviation sigma of normal distribution according to spare part life data;
fitting life data of spare parts with life obeying normal distribution to normal distribution according to reliability theory, wherein the specific process is as follows:
step 11: for all complete data t i Using functions
Figure BDA0003390221940000041
Calculate, recorded as LK i The method comprises the steps of carrying out a first treatment on the surface of the For the truncated data t j Use +.>
Figure BDA0003390221940000042
Calculate, recorded as LK j Wherein->
Figure BDA0003390221940000043
Step 12: all LK is taken i And LK (sum of LK) j Summing to obtain LK
Step 13: solving the mean value estimated value when LK takes the maximum value by using Excel programming solving function, matlab fsolve function and other tools
Figure BDA0003390221940000044
And standard deviation estimate +.>
Figure BDA0003390221940000045
Figure BDA0003390221940000046
And->
Figure BDA0003390221940000047
I.e. the parameters that need to be fitted.
Step 2: and obtaining the expected value of the failure times of the spare part in a given time interval according to the normal distribution.
And (2) calculating expected values of failure times in a given interval (0, t) according to the normal distribution obtained in the step (1), wherein a calculation formula is as follows:
Figure BDA0003390221940000048
in this embodiment, the design numerical calculation method calculates M (t), and includes the following steps:
step 21: the interval (0, t) is divided into N equal parts, and the greater the interval length Δt, i.e., t=n×Δt, the greater the calculation accuracy of M (t).
Step 22: calculating an expected value of the average number of failures
Figure BDA0003390221940000049
Wherein F (t) is a normal distributed cumulative probability density function; t is t i For the position of the ith part Deltat in the interval (0, t), t i =i×Δt。
Step 23: calculating variance
Figure BDA0003390221940000051
Wherein: var [ N [ t ] ] is the variance of the number of failures that occur in the spare part during the time interval (0, t).
Step 3: determining inventory quota for spare parts based on service level of spare parts
Step 31: assuming that there are S positions requiring the use of a spare part, the life of each spare part is L at the time of prediction i Then the average demand for all the positional spare parts after the lapse of time L is
Figure BDA0003390221940000052
Variance is
Figure BDA0003390221940000053
Step 32: calculating inventory quota D using poisson distribution p =P -1 (k%,M s ) Wherein P is -1 () An inverse function representing the Poisson's distribution cumulative density function, k being the service level to be achieved by the spare part, M s Calculating inventory quota D for parameters of poisson distribution using normal distribution N =N -1 (k%,M s ,var[N s (t)]) Wherein N is -1 () An inverse function representing a normal distribution cumulative density function, k being the service level to be achieved by the spare part, M s Is the mean value of normal distribution, var [ N ] s (t)]Is the variance of the normal distribution.

Claims (10)

1. The nuclear power plant spare part demand prediction method based on normal distribution is characterized by comprising the following steps of:
step 1: acquiring a parameter mean mu and a standard deviation sigma of normal distribution according to spare part life data;
step 2: acquiring expected values of failure times of spare parts in a given time interval according to normal distribution;
step 3: inventory quota for the spare parts is determined based on the service level of the spare parts.
2. The nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 1, wherein: step 1 fits the life data of spare parts with life obeying normal distribution to normal distribution according to reliability theory, and the specific process is as follows:
step 11: for all complete data t i Using functions
Figure FDA0003390221930000011
Calculate, recorded as LK i The method comprises the steps of carrying out a first treatment on the surface of the For the truncated data t j Use +.>
Figure FDA0003390221930000012
Calculate, recorded as LK j Wherein->
Figure FDA0003390221930000013
3. The nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 2, wherein: step 1 fits the life data of spare parts with life obeying normal distribution to normal distribution according to reliability theory, and the specific process is as follows:
step 12: all LK is taken i And LK (sum of LK) j Summing to obtain LK.
4. The nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 3, wherein: step 1 fits the life data of spare parts with life obeying normal distribution to normal distribution according to reliability theory, and the specific process is as follows:
step 13: solving the mean value estimated value when LK takes the maximum value by using Excel programming solving function, matlab fsolve function and other tools
Figure FDA0003390221930000014
And standard deviation estimate +.>
Figure FDA0003390221930000015
5. The nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 1, wherein: step 2 calculates expected failure times in a given interval (0, t) according to the normal distribution obtained in step 1, wherein the calculation formula is as follows:
Figure FDA0003390221930000021
6. the nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 5, wherein: and the step 2 is used for calculating M (t), and comprises the following steps:
step 21: the interval (0, t) is divided into N equal parts, and the greater the interval length Δt, i.e., t=n×Δt, the greater the calculation accuracy of M (t).
7. The nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 6, wherein: and the step 2 is used for calculating M (t), and comprises the following steps:
step 22: calculating an expected value of the average number of failures
Figure FDA0003390221930000022
Wherein F (t) is a normal distributed cumulative probability density function; t is t i For the position of the ith part Deltat in the interval (0, t), t i =i×Δt。
8. The nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 7, wherein: and the step 2 is used for calculating M (t), and comprises the following steps:
step 23: calculating variance
Figure FDA0003390221930000023
Wherein: var [ N [ t ] ] is the variance of the number of failures that occur in the spare part during the time interval (0, t).
9. The nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 1, wherein: the step 3 includes the steps of,
step 31: assuming that there are S positions requiring the use of a spare part, the life of each spare part is L at the time of prediction i Then the average demand for all the positional spare parts after the lapse of time L is
Figure FDA0003390221930000024
Variance is
Figure FDA0003390221930000031
10. The nuclear power plant spare part demand prediction method based on normal distribution as claimed in claim 9, wherein: the step 3 includes the steps of,
step 32: using poissonDistribution computing inventory quota D p =P -1 (k%,M s ) Wherein P is -1 () An inverse function representing the Poisson's distribution cumulative density function, k being the service level to be achieved by the spare part, M s Calculating inventory quota D for parameters of poisson distribution using normal distribution N =N -1 (k%,M s ,var[N s (t)]) Wherein N is -1 () An inverse function representing a normal distribution cumulative density function, k being the service level to be achieved by the spare part, M s Is the mean value of normal distribution, var [ N ] s (t)]Is the variance of the normal distribution.
CN202111463258.1A 2021-12-03 2021-12-03 Nuclear power plant spare part demand prediction method based on normal distribution Pending CN116258221A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111463258.1A CN116258221A (en) 2021-12-03 2021-12-03 Nuclear power plant spare part demand prediction method based on normal distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111463258.1A CN116258221A (en) 2021-12-03 2021-12-03 Nuclear power plant spare part demand prediction method based on normal distribution

Publications (1)

Publication Number Publication Date
CN116258221A true CN116258221A (en) 2023-06-13

Family

ID=86682985

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111463258.1A Pending CN116258221A (en) 2021-12-03 2021-12-03 Nuclear power plant spare part demand prediction method based on normal distribution

Country Status (1)

Country Link
CN (1) CN116258221A (en)

Similar Documents

Publication Publication Date Title
US7478347B2 (en) Semiconductor manufacturing apparatus, management apparatus therefor, component management apparatus therefor, and semiconductor wafer storage vessel transport apparatus
CN106777603B (en) intelligent substation relay protection device service life prediction method
CN101320455A (en) Spare part demand forecast method based on in-service lift estimation
CN110320893B (en) Thermal power generating unit equipment maintenance method based on RCM
CN112396215B (en) Intelligent prediction method for self-adaptive interval of residual life of equipment
CN111160741B (en) Disturbance response method for multi-level task configuration of remanufacturing workshop under complex working condition
CN114862296A (en) Spare part inventory management method and device, electronic equipment and computer storage medium
CN105373883A (en) Early-warning method for familial defects of transformer based on data mining
CN106529767A (en) Equipment health status-based secondary equipment risk assessment method
Li et al. Integrated predictive maintenance approach for multistate manufacturing system considering geometric and non-geometric defects of products
CN117094701A (en) Equipment state management system method and system and electronic equipment
CN113723827B (en) Operation risk diagnosis and operation management and control method and system for subway electromechanical equipment
Al-Duais et al. Optimal periods of conducting preventive maintenance to reduce expected downtime and its impact on improving reliability
CN116258221A (en) Nuclear power plant spare part demand prediction method based on normal distribution
CN111799798B (en) Method and system for improving accuracy of future state load flow calculation result
CN106557839B (en) Equipment maintenance strategy optimization method and system based on big data
CN116720686B (en) Production method and system for automatic product scheduling
CN116258218A (en) Nuclear power plant spare part demand prediction method based on index distribution
CN116258223A (en) Nuclear power plant spare part demand prediction method based on gamma distribution
CN116258222A (en) Nuclear power plant spare part demand prediction method based on lognormal distribution
CN112883549A (en) Method for establishing visual condition maintenance model considering random impact influence
CN116777349A (en) Spare part inventory risk management method for information system
CN116258220A (en) Nuclear power plant spare part demand prediction method
CN116258219A (en) Nuclear power plant spare part demand prediction method based on Weibull distribution
CN116128145A (en) Power equipment state maintenance strategy optimization method

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