JP5345184B2 - Maintenance planning method - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a maintenance planning method capable of acquiring the probability of the vicinity of a significant state on a maintenance plan, such as whether or not an inspection is necessary, more precisely than a conventional manner without performing a large amount of arithmetic operations. <P>SOLUTION: Components or units configuring a plant or a product are classified and grouped for each of components or units whose structures are similar, and a part of the components or units is extracted for each group so that the inspection of the distribution L* of damaged states can be performed, and the distribution L'<SB POS="POST">n</SB>of the damaged states created based on the demonstration data of the components or units configuring the plant or product is stored in storage means, and the degree of matching between the distribution L* and the distribution L'<SB POS="POST">n</SB>is determined so that the L'<SB POS="POST">n</SB>whose degree of matching is the highest can be selected, and the damaged states are calculated based on the selected L'<SB POS="POST">n</SB>, and the components or units to be inspected or repaired are selected from among the non-extracted components or units based on the result. <P>COPYRIGHT: (C)2012,JPO&amp;INPIT

Description

  The present invention relates to a maintenance planning method for various plants and products such as nuclear plants, petroleum plants, and chemical plants that require maintenance over a long period of time. Especially, whether or not an accident occurs or whether inspection is necessary, etc. The present invention relates to a maintenance planning method in which the probability in the vicinity of an important state in the maintenance plan can be obtained in more detail than before and can be obtained without performing a large amount of computation.

  Various plants such as nuclear power plants, petroleum plants, chemical plants, and some large-scale products such as aircraft that have a long depreciation period have a large influence on the surroundings at the time of failure, and thus appropriate maintenance is essential. For this reason, it is common to perform periodic maintenance and inspections during the depreciation period, and to perform inspections, repairs, replacements, etc. while considering the life and fatigue level of the parts and units that make up the plant and product. Has been done. However, since such conventional periodic maintenance inspections were performed by simply deciding that they would be performed every fixed year, such as one year, three years, or five years, parts that must be replaced soon. Even if there are units or units, the replacement was carried over until the next periodic inspection, and during that time a failure occurred, causing a large amount of money to be repaired or restored, and lost credit. Conversely, even if there are parts or units that have reached the end of their life and have been replaced, inspections may be performed, or parts or units that do not need to be replaced may be replaced, resulting in unnecessary costs. .

  For this reason, conventionally, it has been studied how it is effective to perform maintenance during a long depreciation period, which is disclosed in Patent Documents 1 and 2, for example. Patent Document 1 relates to a maintenance plan in which a maintenance plan over the entire depreciation period of a plant or product can be planned quantitatively from the viewpoint of cost. The depreciation period of a plant or product is divided into short-term maintenance periods, For each short-term maintenance period, calculate the state of the units and parts that make up the plant and product based on the damage probability at the end of the period, and perform inspection and repair at the start of the short-term maintenance period based on the calculation result A maintenance plan that selects the units and parts to be performed and calculates the inspection costs and repair costs by the calculation means, and repeats them until the depreciation year, thereby selecting the maintenance units and parts and calculating the maintenance costs. Is.

  Patent Document 2 relates to a maintenance plan that predicts the risk cost due to damage, performance degradation, and function stop of equipment constituting the plant, and selects and supports the optimal operation method while comparing with the gain from operation. An event simulation means for calculating an unreliability corresponding to an operation parameter according to an event tree of a plant device based on a failure statistics database stored in advance by associating an event tree of the device and an unreliability with respect to a failure event, and an event tree Therefore, the risk cost calculation means for calculating the risk cost by accumulating the product of the unreliability and the recovery cost, and the operation method determination means for determining the suitability of the operation condition from the comparison between the risk cost and the gain planned to continue operation And the operation method is a plant according to the operation method determined by the fixed means. It relates maintenance plan which includes the operation method instruction means for instructing a specified operating conditions vessel.

JP 2003-99119 A JP 200494663 A

  However, although both the technique disclosed in Patent Document 1 and the technique disclosed in Patent Document 2 can perform a maintenance plan for a plant or equipment, for example, does it lead to an accident or is inspection necessary? Even if you want to find out only the probabilities of states that are important in the maintenance plan, such as whether or not, in detail, random numbers must be generated uniformly. Because it requires a lot of random numbers and therefore requires a large amount of random numbers, resulting in a heavy load on the calculation process, it takes a long time to calculate and requires the use of a high-performance calculation means, so that a maintenance plan is made. The cost of In particular, when Monte Carlo simulation is used for the calculation, since the error is inversely proportional to the square root of the number of samples, the convergence time is very slow. The problem cannot be ignored.

  Therefore, in view of the problems of the prior art, the present invention can determine the probability of the vicinity of an important state in the maintenance plan, such as whether inspection is necessary, in more detail than before, and without performing a large amount of calculations. The purpose is to provide a maintenance planning method.

In the present invention for solving the above problems, and classifying means for grouping and classifying parts and units constituting the plant or product, an inspection means for determining the distribution of the damage state of the components and units L *, plant Ya Storage means for storing a plurality of data of damage state distributions L ′ _n of parts and units constituting the product, and determining the degree of coincidence between the damage state distribution L * and the plurality of L ′ _n , A selection means for selecting L ′ _n having the highest degree of coincidence with L * and a damage state calculation means for calculating the damage state of parts and units are used. a maintenance plan method of selecting the components and units for performing repairs, parts and units constituting the plant or product for each component or unit of similar structure by said classifying means A first step of grouping and extracting a part or a part of each of the groups classified by the classification means by the inspection means to obtain a damage state distribution L *; stores the parts and units plurality of data distribution L _'n made by damaged state based on empirical data of the storage means for, said the distribution of intact condition L * which has been determined by the testing means by said selecting means 'to determine the _n respective degree of coincidence, L * a highest degree of coincidence L' multiple L acquired from the storage means and the second step of selecting _n, the in the second step by said calculation means based on the selected the L _'n in the selection means, a component or unit of the same group withdrawn by said checking means in said first step and said first The dispersion of the damage state of the components and units that have not been removed by the inspecting means modeled from actual data in step, and a third step of calculating the modeled in the simulation, the calculation in the previous SL third step Based on the result of the damage state of the part or unit calculated by the means, the part or unit for inspecting / repairing the part or unit calculated by the calculating means in the first step is selected. .

In the second step, the damage state distribution L * inspected in the first step is compared with a plurality of data of the damage state distribution L ′ _n created based on the demonstration data of the parts and units constituting the plant or product. By determining the degree of coincidence and selecting L ′ _n having the highest degree of coincidence with L *, it is possible to determine which distribution of the proof data is closest to L *. Here, L′ n can be determined based on, for example, the lot of the part or unit, product data, delivery quality document, use place, use condition, and the like.
In the third step, based on L ′ _n selected in the second step, damage to parts or units that are in the same group as the part or unit extracted in the first step and not extracted in the first step. By calculating the state, it is possible to calculate and determine the damage state even for parts and units that have not undergone sampling inspection. That is, a part of many parts is extracted, and the entire distribution state can be estimated from the extracted state.

Further, the selecting means in the second step is characterized rows Ukoto degree of coincidence judgment of L * and L _'n using Kolmogorov over Smirnov test.
The Kolmogorov-Smirnov test is a kind of hypothesis test in statistics, and is effective for examining whether two populations are different based on a finite number of samples.

In addition, the maintenance planning system further obtains a weighted probability density function w weighting the damaged state that requires inspection / repair based on the verification data on the damaged state of the parts and units constituting the plant or product. weighting the probability density function determining means, 0 or more but 1 or less to generate a uniform random number u _n or ultra uniform random number u * n _(hereinafter referred to as a uniform random number, etc.), the weighting probability the uniform random number, etc. The weighted probability variable obtaining means for obtaining a random variable z _n corresponding to the uniform random number or the like, and the plant or product is configured as the cumulative probability of the cumulative distribution function of the weighted probability density function w obtained by the density function determining means. The data relating to the probability of the damage state relating to the parts or units to be stored is stored in the second storage means, and the data relating to the probability of the damage state is stored by the weighting probability calculation device. A weighting probability calculation device for calculating a weighting probability of a damage state of a part or unit constituting the plant or product by Monte Carlo simulation using the obtained z _n as a random number, and weighting based on the weighted probability density function w The damage probability calculating means for calculating a damage value of a part or unit constituting a plant or product from the weight of the damage value obtained by the restoration value W and the weighted probability calculation device. with which the distribution L of the intact condition _'n, the distribution L of the damage probability calculating means the intact condition by calculating based on the probability of the state of the computed units by' and calculates the _n .
This makes it possible to calculate the distribution L ′ _n obtained in detail for the important parts in the maintenance plan without performing a large amount of calculations, so that L * and L ′ _n for the vicinity of the important parts in the maintenance plan can be calculated. Comparison accuracy is improved.

  As described above, the maintenance planning method of the present invention has the advantage that the probability of the vicinity of an important state in the maintenance plan, such as whether inspection is necessary, can be obtained in more detail than before, and without performing a large amount of computation. There is.

The maintenance plan method in Example 1 is shown with the block diagram. It is the flowchart which showed the flow which calculates | requires the weighted probability density function w. FIG. 3A is a graph showing an example of the weighted probability variable w, and FIG. 3B is a graph showing the cumulative distribution function of the weighted probability variable w shown in FIG. FIG. 4A is a graph showing an example of a probability density function of a damage state related to parts and units constituting a plant or a product, and FIG. 4B is an accumulation of probability density functions shown in FIG. It is a graph showing a distribution function. The maintenance plan method in Example 2 is shown with the block diagram. It is a conceptual diagram of L' _n accommodated in the memory | storage means.

  Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the drawings. However, the dimensions, materials, shapes, relative arrangements, and the like of the components described in this embodiment are not intended to limit the scope of the present invention unless otherwise specified, but are merely illustrative examples. Not too much.

FIG. 1 is a block diagram illustrating a maintenance planning method according to the first embodiment. 1 is a weighted probability density function determining means for determining a weighted probability density function w that gives a weight to a desired damage state based on demonstration data on the damage state of parts and units constituting the plant or product. Based on a weighted probability density function w, a weighted random variable generator that generates a weighted random variable z _n that is a random number weighted to the desired damage state, 3 is the weighted random variable z _n as a random number, Weighting for calculating the weighting probability of the damage state of the parts or units constituting the plant or product weighted to the desired damage state by the Monte Carlo simulation from the data regarding the damage state probability of the plant or product previously stored in the database 6 4 is based on the weighted probability density function w. Damage probability calculation means for calculating a restoration value W for restoring weighting, and calculating a damage state probability of a part or unit constituting a plant or product from the restoration value W and the weighting probability of the damage state; It is a maintenance planning means that calculates and calculates the damage status of parts and units that make up the plant and product from the probability of damage status of the components and units, and selects the parts and units to be inspected and repaired based on the calculation results .
Hereinafter, it will be described in detail along the flow of the block diagram shown in FIG.

  First, the weighted probability density function determining means 1 determines the weighted probability density function w. FIG. 2 is a flowchart showing a flow for obtaining the weighted probability density function w. In FIG. 2, first, a temporary weighted probability density function w ′ is appropriately set in step S11. In the next step S12, a probability density function I is obtained from the provisional weighted probability density function w ′ set in step S11 by using Monte Carlo simulation.

In the next step S13, the probability variable at the peak of the probability density function I is compared with a probability variable important in the maintenance plan to determine whether or not they match. Here, the coincidence may be a complete coincidence, but an appropriate range may be provided and the coincidence may be determined within the range.
If YES is determined in the step S13, the process proceeds to a step S14. In the step S14, the provisional weighting probability density function w ′ at the determination time is determined as the weighting density function w. If NO is determined in step S13, the process proceeds to step S15 described later.
Here, a random variable important in the maintenance plan is a random variable corresponding to a damage state of a part or unit constituting a plant or product that needs to be examined in detail. Depending on the product, part or unit and the purpose of use, for example, if the part is cracked, the crack depth that requires inspection, the crack depth that requires repair, the crack depth that leads to the occurrence of an accident, etc. Determined based on empirical data. Damage conditions that need to be examined in particular in this way are determined based on demonstration data, so a lot of demonstration data is collected, for example, crack depth that requires inspection, crack depth that requires repair, etc. Is preferably used.

Next, in step S15, it is determined whether or not the peak random variable of the probability density function I is larger than the random variable important in the maintenance plan. If YES is determined in step S15, the probability variable of the probability density function I has a large probability variable. Therefore, in step S16, the temporary weighting variable w ′ is changed so that the probability variable of the peak of the probability density function I decreases. Return to S12. On the other hand, if NO is determined in step S15, the peak probability variable of the probability density function I is small. Therefore, in step S17, the temporary weighted probability variable w ′ is changed so that the peak probability variable of the probability density function I decreases. The process returns to step S12 again.
Thus, the weighting density function w is determined by repeating Steps S11 to S17.

Then weighting random variable generator 2 generates a weighted random variable z _n. FIG. 3A is a graph showing an example of the weighted random variable w, where the vertical axis represents the probability and the horizontal axis represents the random variable. FIG. 3B is a graph showing the cumulative distribution function of the weighted probability variable w shown in FIG. 3A. The vertical axis represents the cumulative probability, and the horizontal axis represents the random variable.

In the weighted random variable generator 2, first, the weighted probability density function w determined by the weighted probability density function determining means 1 as shown in FIG. 3A is accumulated as shown in FIG. 3B. Convert to distribution function. Then generates a 0 or more and 1 or less of n uniform random number u _n or ultra uniform random number u * n _(hereinafter represented by uniform random numbers or the like as called u _n), shown respectively in FIG. 3 (B) To the cumulative probability of the cumulative distribution function. FIG. 3B shows only u _i which is one of _n uni. Thus the random variable z _n corresponding to the u _n which fit the cumulative probability obtained from the cumulative distribution function. FIG. 3B shows only a random variable z _i corresponding to u _i .
Z _n obtained in this way is defined as a weighted random variable.

Next, the weighting probability calculating means 3 calculates and calculates the weighting probability. The weighted random variable z _n generated by the weighted random variable generating device 2 is used as a substitute for a random number, and the probability variable x _n (w _{n) is determined} from the probability of damage states related to parts and units constituting the plant and product stored in the database 6. ) And calculating and calculating the weighting probability of the damage state of the parts and units constituting the plant or product by Monte Carlo simulation.

Here, the random variable x _n (w _n ) is generated as follows. FIG. 4A is a graph showing an example of a probability density function of damage states related to parts and units constituting a plant or a product, where the vertical axis represents probability and the horizontal axis represents random variables. FIG. 4B is a graph showing the cumulative distribution function of the probability density function shown in FIG. 4A, where the vertical axis represents the cumulative probability and the horizontal axis represents the random variable. Each of the weighted random variables z _n is applied to the cumulative probability of the cumulative distribution function shown in FIG. FIG. 4B shows only z _i which is one of _n z _n . In this way, the random variable x _n (w _n ) can be generated by determining the random variable corresponding to z _n fitted to the cumulative probability from the cumulative distribution function.

In addition, here, the probability of the damage state regarding the parts and units constituting the plant or product stored in the database 6 is, for example, the probability of occurrence of damage, the probability of occurrence of the damage that has occurred, the probability of detecting the damage that has occurred, etc. Of which one or more probabilities are used.
Further, the weighting probability of the damage state of the parts or units constituting the plant or product calculated by the weighting probability calculating means 3 is the weighting probability that requires inspection, the weighting probability that requires repair, or an accident. The weighting probability is given.

However, the weighted probability of the damage state of the parts and units constituting the plant and product obtained by the weighted probability calculation means 3 by the Monte Carlo simulation using the weighted random variable z _n which is a biased variable in this way as a substitute for the random number. Does not represent the probability of damage of the parts and units constituting the plant or product. Therefore, the damage state probability calculation device 4 calculates the probability of damage states of the parts and units constituting the plant and product by restoring the weights from the weighted damage state probabilities of the parts and units constituting the plant and product. calculate.

  The damage state probability calculating device 4 restores the weights as follows and calculates and calculates the damage state probability of the parts and units constituting the plant and product. Here, as an example of the probability of damage states of parts and units constituting the plant or product, an example of calculating the probability of occurrence of an accident by calculating the depth of cracks occurring in the parts will be given.

また、ｎ個の重み付け確率変数ｚ_ｎを数２で表すことができる。

この時、全体（ｘ_１、ｘ_２、・・・ｘ_ｎ）の重みを復元するための復元値をＷとすると、Ｗは数３で表すことができる。

First, n random variables to be used in the Monte Carlo simulation can be expressed by Equation 1.

Further, n weighting probability variables z _n can be expressed by Equation 2.

At this time, if the restoration value for restoring the weight of the whole (x ₁ , x ₂ ,... X _n ) is W, W can be expressed by Equation 3.

亀裂深さＬ_０以上で事故が発生する、つまり亀裂深さＬ_０が事故発生に至る境界深さである場合、事故が発生するシミュレーション結果数Ω^Ｘと、事故が発生しないシミュレーション結果数Ω^０はそれぞれ数５及び数６を用いて求めることができる。

従って、事故発生確率Ｐ^Ｘは前記Ω^Ｘ及びΩ^０を復元値Ｗで復元して数７を用いて求めることができる。

Moreover, the crack depth L of the simulation result can be expressed by Equation 4,

When an accident occurs at a crack depth L ₀ or more, that is, when the crack depth L ₀ is a boundary depth leading to the occurrence of the accident, the number of simulation results Ω ^X where the accident occurs and the number of simulation results Ω ⁰ where no accident occurs Can be obtained using Equation 5 and Equation 6, respectively.

Therefore, the accident occurrence probability P ^X can be obtained by using Equation 7 after restoring the Ω ^X and Ω ⁰ with the restoration value W.

  Next, the maintenance plan means 5 calculates the damage state probability calculation means 4 and calculates the damage state of the parts and units constituting the plant and product from the damage state probability of the parts and units constituting the plant and product. And a part or unit to be inspected / repaired is selected based on the calculation result.

  The maintenance planning method of the first embodiment can be applied to any part or unit as long as it is a maintenance plan for parts and units constituting a plant or product. It is preferable to use a part or unit that is connected to a part or unit that causes a large recovery cost if damage occurs because of the connection, and the repair criteria can be adjusted. The repair criteria can be adjusted, for example, when there are multiple inspection methods, and when repairs are made according to the crack length found in the inspection, it is possible to determine the criteria for the crack length that the repair is deemed necessary. .

  As such a preferable specific example, for example, crack inspection of a tail cylinder of an incinerator of a gas turbine can be given. In the inspection of cracks in the incinerator of a gas turbine, how long is the crack depending on the inspection method, and what is the error between the actual crack length and the crack length detected in the inspection? The basic data of the relationship is organized, and the criteria for the depth of cracks detected in inspections that require repair are determined. Therefore, it is effective to perform simulations on a plurality of inspection methods and repair methods, and to select effective inspection / repair measures by giving weights in the vicinity of an important state in terms of maintenance management whether or not serious damage will occur.

FIG. 5 is a block diagram showing the maintenance planning method according to the second embodiment. 11 is a classification means for classifying parts and units constituting a plant or product into parts and units having similar structures and grouping them, and 12 is a damaged state by extracting a part of the parts and units for each of the classified groups. Is an inspection means for inspecting the distribution L * of the gas, and reference numeral 13 denotes a plurality of damage states created based on the damage state distribution L * and verification data of parts and units constituting the plant and product stored in the storage means 14. Is a selection means for determining the degree of coincidence with the distribution L ′ _{n of} L * _n and selecting L ′ _n having the highest degree of coincidence with L *, and 15 is based on L ′ _n selected by the selection means 14. Computation means for computing the damage state of the parts and units that are in the same group as the parts and units extracted by the inspection means 12 and not extracted by the inspection means 12; A means for selecting the components and units to inspect and repair based on the wound state.
Hereinafter, a detailed description will be given along the flow of the block diagram shown in FIG.

  First, the classification means 11 classifies the parts and units constituting the plant and products into parts and units having similar structures.

Next, the inspection unit 12 extracts a part or part of each unit classified by the classification unit 11 and obtains a damage state distribution L * by inspection. The number of samples to be extracted is not particularly limited, but if it is too large, the cost of extraction and inspection will increase, and if it is too small, L' _n having the highest degree of coincidence with L *, which will be described later, cannot be selected accurately. Since the management accuracy is lowered, it is preferable to determine the number of sampling according to the required cost and the maintenance management accuracy.

Next, a plurality of damages created by the selection means 13 based on the verification data of the unit of the parts constituting the plant and the product are stored in the storage means 14 with the damage state distribution L * obtained by the inspection means 12. The degree of coincidence is determined by comparison with the state distribution L ′ _n . Kolmogorov-Smirnov test, which is a kind of hypothesis test, is used to determine the degree of coincidence.

Here, the damage state distribution L ′ _n created based on the verification data of the parts and units constituting the plant and product stored in the storage means 14 will be described. FIG. 6 is a conceptual diagram of L ′ _n stored in the storage means 4. As shown in FIG. 6, M types of data (1) to (M) are stored in the storage unit 14, and each has different damage occurrence time probability distributions and occurrence frequencies. And wherein (1) ~ (M) of the M types of data, each said damage occurring timing probability and frequency related parameter _{(m i, σ i) (} i integer from 1 to M) and the parameter _{(m i,} σ _i ) having a distribution function L ′ (m _i , σ _i ) corresponding to a damaged state, and L ′ (m _i , σ _i ) corresponds to L ′ _n .

Moreover, the damage occurrence probability distribution and occurrence frequency shown in FIG. 6 are determined on the basis of demonstration data. For example, in the case of a product part, the part design data, delivery quality data, production lot, use condition Etc.
Further, in FIG. 6, in order to obtain the parameters (m _i , σ _i ) and the damage state distribution L ′ (m _i , σ _i ), the damage occurrence time probability distribution and the occurrence frequency are used. Alternatively, for example, a probability distribution of damage progress may be added to obtain parameters and damage state distributions from three or more types of data.
Further, the distribution of the parameters and the damage state may be obtained based on the damage state probability of the parts and units constituting the plant or product obtained by the damage state probability calculation means 4 shown in FIG.

Next, the calculation means 15 is in the same group as the parts and units extracted by the inspection means 12 based on L ′ _n selected by the selection means 14 and is damaged by the inspection means 12 and has not been extracted. The state is calculated by calculation, and the maintenance plan means 16 selects a part or unit to be inspected / repaired based on the damage state of the part or unit calculated by the calculation means.

The maintenance planning method of the second embodiment can be used for any part or unit as long as it is a maintenance plan for parts and units constituting a plant or product. A hole is mentioned. Usually, there are dozens of cooling holes in the incinerator tail tube, and cracks may occur from each hole. When the maintenance planning method of the second embodiment is not applied, 100% inspection will be performed, but the ease of crack generation and growth varies depending on the state where the incinerator is placed, that is, the positional relationship between the incinerator and the cooling hole. There is a track record. Therefore, by preparing a plurality of L ′ _n corresponding to the results and inspecting some cracks in the cooling holes, it is possible to calculate and calculate cracks in other cooling holes. In addition, since it is also known what kind of positional relationship between the incinerator and the cooling holes, cracks are likely to grow, so by inspecting the cooling holes where cracks are likely to grow, the cracks in the remaining cooling holes Can model from the actual data what kind of dispersion corresponds to the inspection result. By performing simulation using the modeled model, the growth rate to the depth corresponding to the crack damage is compared according to the sampling rate, and the simulation is used for selecting the optimal rate, selecting the interest rate, and selecting the sampling site. It becomes possible.

  Probability in the vicinity of important states in maintenance plans, such as whether inspection is necessary, can be determined in more detail than before, and without requiring a large amount of computation, so it can be used for maintenance management plans for parts and units of plants and products. Can be used.

DESCRIPTION OF SYMBOLS 1 Weighted probability density determination means 2 Weighted random variable generator 3 Weighted probability calculation means 4 Damage state probability calculation means 5 Maintenance plan means 11 Classification means 12 Inspection means 13 Selection means 14 Storage means 15 Calculation means 16 Maintenance plan means

Claims (3)

Classification means for classifying and grouping parts and units that make up a plant or product,
An inspection means for obtaining a distribution L * of the damage state of the component or unit;
Storage means for storing a plurality of data of damage state distributions L ′ _n of parts and units constituting the plant or product;
'To determine the respective degree of coincidence _n, the highest degree of coincidence L and L *' distribution L * and the plurality of L of the intact condition and selecting means for selecting the _n,
Damage state calculation means for calculating the damage state of parts and units;
A maintenance planning method for selecting parts and units to be inspected and repaired using a maintenance planning system that performs the selection means by simulation,
Said classification means grouped by classifying parts and units constituting the plant or product for each component or unit of similar structure by, some of the parts and units for each of the groups classified by the classifying means by said inspection means A first step to extract the damage state distribution L *;
Stores the plurality of data distributions L _'n damage state were prepared based on the empirical data of the parts or units constituting the plant or product in the storage means, the damage state which has been determined by the testing means by said selecting means A second step of determining the degree of coincidence between the distribution L * and each of the plurality of L ′ _n acquired from the storage means, and selecting L ′ _n having the highest degree of coincidence with L *;
Based on L ′ _n selected by the selection means in the second step by the calculation means, the parts and units extracted by the inspection means in the first step are in the same group and the first A third step of modeling the distribution of damage states of parts and units not extracted by the inspection means in the step from actual data, and calculating the modeling by simulation ,
Based on the results of the damage state before Symbol third parts and units which are calculated by the computing means in step, component Ya to inspect and repair of the computed components and units by the computing means in said first step A maintenance planning method characterized by selecting a unit. 2. The maintenance plan determination method according to claim 1, wherein, in the second step, the selection means determines the degree of coincidence between L * and L′ _n using a Kolmogorovius Mirnov test. The maintenance planning system further includes:
A weighted probability density function determining means for obtaining a weighted probability density function w weighted to a damaged state that needs to be inspected / repaired based on proof data on a damaged state of a part or unit constituting a plant or a product;
A uniform random number u _n of 0 or more and 1 or less or a super-uniform random number u * _n (hereinafter referred to as a uniform random number or the like) is generated, and the uniform random number or the like is obtained by the weighted probability density function determining means. A weighted probability variable obtaining unit that obtains a cumulative probability of a cumulative distribution function of the weighted probability density function w and obtains a random variable z _n corresponding to the uniform random number, and the like;
Data relating to the probability of damage state relating to the parts and units constituting the plant or product is stored in the second storage means, and the calculated z _n is used as a random number from the data relating to the probability of damage state by a weighted probability arithmetic unit. A weighting probability calculation device for calculating the weighting probability of the damage state of the parts and units constituting the plant and product by simulation;
The weighting based on the probability density function w determined weighting reconstruction value W,該復original value W and the probability of damage to the state of the parts and units constituting the plant or product from the weighting probability of the obtained intact condition by said weighting probability calculation unit The damage probability calculating means for calculating
The distribution L damaged state _'n, the distribution L of the damage probability calculating unit state the intact condition by calculating a probability based on the of the computed units by' claim 1, characterized in that to calculate the _n or 2. The maintenance plan method according to 2.

Images