Estimation of the Chemical Device Failure Probability Based on Monte Carlo Simulation Algorithm
窦站,蒋军成,朱常龙,张明广,赵声萍
DOI: https://doi.org/10.3969/j.issn.1009-6094.2012.06.034
2012-01-01
Applied Mechanics and Materials
Abstract:In this paper,we would like to propose a new method based on Monte Carlo Simulation(MCS) theory for chemical device failure probability estimation,particularly that of the chemical facilities failure characteristic of domino effect.As is seen,in petroleum chemical works,it is often the case that chemical device failure often leads to domino effect,which might be not only affected by the internal factors,but also by external factors,such as domino effect brought about from the outside.Therefore,it is of great necessity to consider the domino effect when we come to the case of such kinds of device failure.Needless to say,chemical equipment accidents with domino effect usually take place in a complicated manner due to the diversified reasons or unexpected or expected causes.Great knowledge and vast experience are needed to make a sound assessment or estimation of such events based on the probability theory.Even if the situation was made clear,errors in calculation,such as round-off error,things would remain suspicious because the data used for such assessment may not be adequate in size and accuracy for the number of likely scenarios is too huge.Therefore,at the first stage(initialization),it is necessary to specify the data and information of the number of pieces of equipment(n),the failure rates(λi) for any equipment in isolation,escalation probabilities(Pji),number of iterations(N),time step(Δt) and the final time(tf),etc.It is also necessary to initialize a matrix,named F,which can show the failure probability of each component at each respective time step from 0.Then,for each time step and each component(i)(which can be chosen randomly,which should not be chosen in favor of the higher failure chance of any particular component,rather than done objectively);And,next,it is needed to work out the instantaneous probability of failure(Pi) and compare the probability with a randomly generated number(r1),in case that Pi r1 means that the component under study fails at the time step,whose related position in F can be turned to one.In this case,it is assumed that in each Δt as a component is likely to fail for the sake of being too simple.In the innermost loop component,j is to be checked it is likely to cause any domino effect,i.e.for each component(j) if Pji is greater than a randomly generated number(r3),when the Domino effect is prone to occur(based on the concept of MCS) and therefore its position in F can be changed to one.Such tracing steps should be repeated for N times to make up for the stochastic nature of MC algorithm.For a random number to be quoted,it is necessary to use the algorithm in two steps,a predefined function in the software named MATLAB.Since the random numbers involved here are by nature a set of variables with their values uniformly distributed between 0 and 1,different algorithms are to be used to generate such random numbers.And,as a result,we have made efforts in assessing the applicability and flexibility of the method while applying it for our estimation.From what is described above,it can be seen that the simulation method we have developed enjoys great advantages in comparison with the purely analytical probability methods.The major advantage of our method is independent of the accuracy of the results of the complexity of the system.In addition,when it is used,it is possible to work out the failure probability as a function of time.