A Novel Underlying Algorithm for Reducing Uncertainty in Process Industry Risk Assessment

: Normal fuzzy fault tree is a classic model in the field of process industry risk assessment, and it can provide reliable prior knowledge for machine learning. However, it is difficult to adapt the traditional approximate calculation method to highly nonlinear problems, and this may introduce model uncertainty. To solve this problem, this study proposes an accurate calculation algo-rithm. In the proposed algorithm, first, an exact α -cut set of normal fuzzy fault tree is derived according to the exact calculation formula of normal fuzzy numbers and in combination with the cut-set theorem. Subsequently, the relationship between the membership function and the exact cut set is derived based on the representation theorem. Finally, according to the previous derivation, the coordinates of the point on the exact membership curve are found within the range of x from 0 to 1. Based on this, an accurate membership graph is drawn, the membership curve is evenly divided with the area enclosed by the x -axis, and the fuzzy median is obtained. The results of the two chemical accident cases demonstrate that the proposed algorithm has a strong ability to handle uncertainty and can significantly reduce the uncertainty of the process industry risk assessment results. The results also reveal that the superiority of the accurate calculation algorithms becomes more obvious when the mean failure probability of basic events is larger or the accident tree is more complex. This study provides a high-accuracy underlying algorithm for process industry risk assessment, and it is of great value for improving system security.


Introduction
The advent of the artificial intelligence era has promoted the development of process industry risk assessment [1,2].As a classic method in the field of process industry risk assessment, the fuzzy fault tree has been widely applied in various fields [3,4].
However, Yazdi et al. [5] pointed out that the existence of some factors leads to uncertainty in fuzzy fault tree risk assessment results, including the complexity of the object being evaluated, the insufficiency of the required data, the hesitation factors in the decision-making process, and the approximations and simplifications in the modeling process.Furthermore, uncertainties are transmitted from the lowest level of the fault tree to the highest level of the fault tree, seriously affecting the credibility of the fuzzy fault tree risk assessment results.Therefore, uncertainty has become a key issue that needs to be discussed in the fuzzy fault tree risk assessment process [6].
In order to better handle uncertainty in the case of insufficient failure data, the Bayesian theorem is widely used as an analytical tool to update probability distributions.Soltanali et al. [7] established a fuzzy fault tree-Bayesian model to optimize complex automobile manufacturing equipment and improve the accuracy of prediction based on two types of uncertainty factors: insufficient data and system complexity.Kaushik et al. [8] further proposed a similarity aggregation method based on intuitionistic fuzzy fault tree and Bayesian network, which promptly update the failure probability of new evidence.Zhao et al. [9] used fuzzy sets to calculate the probability of basic events and transformed the fault tree into a Bayesian network using Noisy-OR gate, solving the problem of the absolute description of the relationship between basic events and intermediate events.Ding et al. [10] combined FMEA, FTA, and Bayesian analysis to fill the gap in risk assessments of preheating and exclusion systems under uncertain conditions, expressing the polymorphism of events and the uncertainty logic relationship between events.The combination of fuzzy fault tree and Bayesian network can allow for the calculation of the prior probability of each failure based on historical data, thus being able to handle objective uncertainty.
To resolve the uncertainty of information input in fuzzy fault tree analysis, evidence theory and expert opinion are also good choices.Pang et al. [11] analyzed the uncertainty of systems based on fuzzy fault trees and evidence theory, overcoming the uncertainty of imprecise evidence.Zaib et al. [12] used expert opinions to address the uncertainty of human error events.However, while these obtain the probability of basic events, uncertainty related to knowledge is likely to occur.Therefore, the subjectivity of expert opinions has become an inevitable issue in the handling of uncertainty in fuzzy fault trees.Celik et al. [13] used fuzzy set theory to address the subjective uncertainty in expert heuristic processes, emphasizing that the use of linguistic terms is an appropriate choice for dealing with the uncertainty of basic event probabilities.Purba et al. [14] used the FPFTA method to quantify the uncertainty of experts' overall understanding.
For some complex systems, experts often cannot provide accurate membership functions for system component failures.In order to overcome the subjective uncertainty of experts in fuzzy fault tree risk assessment, Kaushik and Kumar [15] developed a new IF importance measurement method that effectively measures the uncertainty of basic events by providing cut intervals in an -intuitionistic fuzzy environment.Hu et al [16] considered the hesitant factors in uncertainty quantification and proposed an IF importance measurement method that is beneficial for identifying key basic events.Lu et al. [17] proposed an improved similarity aggregation method based on the butterfly optimization algorithm to address the uncertainty problem of basic event failure probability in threedimensional EVAC systems and improve the accuracy of aggregation results.Garg et al. [18] proposed a fuzzy risk assessment method based on the granularity Z-number to improve the utilization of information and the degree of uncertainty in management decisions.Feng et al. [19] proposed a fuzzy importance sampling method to improve the computational efficiency of system failures with fuzzy uncertainty.
From the above analysis, it can be seen that the fuzzy fault tree analysis method produces two uncertainties in the process of outputting risk assessment results [20]: one is the uncertainty of the input quantity (parameter), which is related to the true value of the input x, and the other is structural model uncertainty, that is uncertainty related to model errors, which is related to fuzzy fault tree membership models and approximate calculations.
The above literature solves the uncertainty of information input in the fuzzy fault tree model from different perspectives.To a certain extent, the uncertainty of expert understanding in fuzzy fault tree risk assessment is reduced, and the uncertainty of model input parameters in the process of fuzzy fault tree risk assessment is improved.
However, there is currently little research on the uncertainty of the fuzzy fault tree's own structural model.The uncertainty of the fuzzy fault tree's own structural model stems from the gap between knowledge and practice.It may take the form of an insufficient understanding of phenomena that already occur in the system or the simplification of complex phenomena leading to "inappropriate" assumptions about the structure of the model.
The purpose of this study is to propose a new algorithm based on the cut-set theorem and the representation theorem, with the aim of solving the problems of difficulty in making accurate calculations and determining uncertainty in approximate calculations during the normal fuzzy fault tree risk assessment process.The formula derivation and proof can provide theoretical support for the implementation of the algorithm.The method of drawing a precise membership graph can help to quickly obtain accurate risk assessment results.This will reduce the uncertainty of the process industry risk assessment and improve the accuracy of risk assessment.It will provide a high-precision underlying algorithm for process industry risk assessment.

Accurate Calculation Algorithms
To overcome the shortcomings of traditional approximate calculations, an accurate calculation algorithm is proposed in this section, and a theoretical proof process is provided.

Algorithm Principles
The principle of the accurate calculation algorithm is as follows: Firstly, convert the normal fuzzy numbers into classical sets based on the cut-set theorem.Secondly, discretize the normal fuzzy numbers based on the representation theorem, and draw accurate membership curves based on discrete data.Finally, calculate the area of the graph bounded by the membership curve and the x-axis, and further calculate the fuzzy normal median and fuzzy importance.

Theoretical Basis of the Algorithm
Normal fuzzy fault tree risk assessment is an effective method that can fully consider the opinions of multiple experts.Hence, m experts are invited to evaluate the possibility of basic event Xi, and then the failure probability Qi (i.e., Qi1, Qi2, ⋯ Qim) of Xi is obtained.Based on this, the mean value of the failure probability is calculated as qi, and the variance is σ i .According to Equation (A1), the membership function A  i (x) of Qi can be expressed as In production practice, the multiplication of failure probabilities represented by normal fuzzy numbers is too complicated, because the fault tree is complex and there are many basic events.Therefore, approximate calculation formulas are commonly used, as shown in Equation (2), However, approximate calculations can lead to uncertainty in the assessment results.Therefore, an accurate calculation algorithm is proposed in this paper.The derivation of the relevant formulas is carried out below in two cases:

A Normal Fuzzy Fault Tree with Only OR Gates
According to Equation (1), the normal fuzzy number A 1-Q i (x) of the probability of the non-occurrence of the basic event is obtained using According to Equations (A2) and (3), let A 1-Q i (x) = α, and the α-cut set of 1-Q i is calculated as Based on the fuzzy interval number operation rule, we have According to Equations ( 4) and ( 5), the exact α-cut set is obtained for the probability of basic events not occurring simultaneously, Based on the concept of bounded closed fuzzy numbers and Equation (A5), it can be seen that the membership function in which n basic events do not occur at the same time can be expressed as According to Equations (A6)-(A8), m A and n A are further derived (9) and then   (x) and R A (x) are derived based on Equations (A9) and (A10).

A Normal Fuzzy Fault Tree with OR and AND Gates
Using Formula (A12) and the cut-set theorem, the exact cut set calculation formulas for   '(x) and   '(x) are derived.The derivation results are as follows: The formula for the cut set of membership function   '(x): ∑ ∏ (q iσ i ×√-lnα) The formula for the cut set of membership function   '(x): ∑ ∏ (q i +σ i ×√-lnα) The proof of the monotonic increasing property of the cut set of   '(x) is as follows: increases monotonically.The proof of the monotonic decreasing property of the cut set of   '(x) is as follows: Because d ∏ (q i +σ i ×√-lnα) we have thus, ∑ ∏ (q i +σ i ×√-lnα) Secondly, based on the derived Equations ( 14) and ( 15), as well as the monotonicity of the cut set, we obtain the coordinates of points on the exact membership curve.

Algorithm Execution Process
The specific execution process of the algorithm is shown in Figures 1 and 2. The algorithm is mainly divided into three steps.The first step is to calculate the cut-set interval under a given α.The second step is to discretize the result of multiplying the normal fuzzy numbers according to the cut-set interval and determine the ordered arrays.The third step is to draw a membership curve based on the ordered arrays and determine the fuzzy median value based on the area bounded by the membership curve and the x-axis.Firstly, the algorithm in Figure 2 needs to compute the expression for the failure probability of the top event.Secondly, if there are duplicate probabilities in the probability expression, then the duplicate probabilities need to be eliminated.Thirdly, an accurate cut-set expression for the probability of the top event needs to be expressed.The remaining steps, shown in Figures 1 and 2, are essentially the same.
The new algorithm adopts a data equal-width discretization method.It is necessary to reasonably select the number of discretization intervals based on the value ranges of α and x.This can avoid information loss or redundancy due to too few or too many intervals.

Empirical Analysis
To verify the effectiveness and feasibility of the algorithm, this study selects two types of fault trees for analysis, with the top events being electrostatic accumulation and liquid ammonia leakage.The failure probability of basic events in the fault tree is scored by experienced experts in the chemical industry.
This section includes the following content: Firstly, the execution process of expert scoring is introduced.Then, after obtaining the expert scoring results, the normal characteristics of the failure probability of basic events are verified, and the mean and variance of basic events are obtained.Finally, a new algorithm and approximate calculation method are used to draw a membership curve of the failure probability of top events and intermediate events; following this, the risk assessment results are compared and analyzed.

Expert Introduction and Scoring Rules
The introduction of experts is shown in Table 1, and the scoring rules for experts are shown in Table 2.

Failure Probability
Possibility of Occurrence α × 10 −1 very easy to occur α × 10 −2 easy to occur α × 10 −3 not easy to occur α × 10 −4 more difficult to occur α × 10 −5 extremely difficult to occur α × 10 −6  almost never happens In Table 2, α represents the membership degree of danger at the same probability level.The value of α is shown in Table 3.
Note that Tables 2 and 3 must be used jointly.

Case 1 Introduction
Case 1 selects electrostatic accumulation during oil storage as the top event for analysis [26] and analyses the electrostatic accumulation due to the relative movement of oil flow, filtration, mixing, spraying, sprinkling, flushing, perfusion, shaking, etc.The fault tree has three layers, including two intermediate events and eight basic events, as shown in Figure 3. Based on this, considering the fire source factor, an accident tree with an electrostatic spark explosion as the top event is drawn, as shown in Figure 4.
The nine basic events are as follows: X1-too fast refueling speed; X2-rough tank wall; X3-oil hitting wall; X4-personnel operation error; X5-friction between oil and air; X6-metal floating objects in oil; X7-the measuring tool does not meet the requirements; X8-short standing time; and X9-ignition source.

Expert Opinions on Case 1
Based on Tables 2 and 3, six experts evaluated the possibility of basic events, and the results are shown in Tables 4 and 5.  2, 4, and 5, the failure probability of the basic event is obtained, as shown in Table 6.After obtaining the failure probability of the basic event in Case 1, the K-S (Kolmogorov-Smirnov) test is performed.It is checked whether the probability of failure conforms to a normal distribution, and the test results are shown in Table 7.In Table 7, it can be seen that the exact significant consistency of Xi (i = 1~9) is greater than 0.05.So, the assumption that the probability of all basic events occurring follows a normal distribution is accepted at a significance level of 0.05.

Case 1: Drawing of Membership Curve
Based on the new algorithm (Figures 1 and 2), a failure probability membership curve of the top event is drawn considering the following situations: all events occur; except Xi (i = 1~9), all events occur; X1 and X4 occur simultaneously; and X2, X3, X5, X6, X7, and X8 occur simultaneously.Based on the membership curve, mTXTAe, mTXTBe, mTxie (i = 1~8), mTXM1e, and mTXM2e are determined.The drawing results are shown in Figures 5-16.
Subfigure (a) is a graph of the membership curves obtained by simulation using the new algorithm; subfigure (b) is a graph of the membership curves obtained by the traditional approximation calculation; and subfigure (c) is a comparison of the results obtained by the two calculation methods.The meanings of the letters in the subfigures refer to the explanations in Figure A1 in the Appendix A.4.

Case 1: Discussion and Analysis
According to Figures 5-16, firstly, the fuzzy median value of the accurate algorithm simulation and approximate calculation is obtained.Secondly, the absolute uncertainty reduction value of the fuzzy median value is computed.Finally, the relative uncertainty reduction value is calculated.The results are shown in Table 8.It can be seen in Table 8 that, after the simulation using accurate calculation algorithms, the uncertainty of both intermediate events and top events is reduced.The uncertainty reduction values of XM1 and XM2 are 4.75% and 5.75%, respectively.The uncertainty reduction value of the top event (electrostatic accumulation) is 7.53%.Considering all cases, the uncertainty reduction value is between 4.75% and 11.05%.
The above data illustrate that, by using approximate calculations, system uncertainty can be greatly propagated.After using the new algorithm, this uncertainty is significantly reduced.The uncertainty reduction value significantly increases from the intermediate layer to the top layer of the fault tree.
When the number of basic events occurring simultaneously is 2, 6, 7, or 8, the uncertainty of the fuzzy median value is significantly reduced.This conclusion can be reached by observing the fuzzy median of the top events and the intermediate events.
From the results, it can be seen that, in this case, the fault tree containing mixed gates has lower uncertainty than the accident tree containing only OR gates.This is because, after the transformation with AND gates, the number of elementary events that are multiplied becomes much less, and the probability of elementary events that are multiplied at the same time becomes much smaller.
In Figure 17, X i ′ indicates that basic event Xi does not occur.The trend of the changes in the fuzzy median obtained from the algorithm simulation and approximate calculation is basically consistent.However, the fuzzy median calculated using the accurate algorithm is significantly smaller.This shows that there is large uncertainty in the approximate calculation.According to Equation (A11) and the data in Table 8, the fuzzy importance of Xi (i = 1~8) is first calculated.This includes the traditional approximate calculation results and the improved accurate calculation results.Second, the absolute uncertainty reduction value of fuzzy importance is calculated.Finally, the fuzzy importance relative uncertainty reduction value is computed.The results are shown in Table 9.Based on the data in Table 9, a fuzzy importance comparison is drawn, as shown in Figure 18.It can be seen in Table 9 and Figure 19 that, after the simulation using the accurate algorithm, the uncertainty reduction effect of the fuzzy importance of X2, X4, X5, and X6 is the most obvious.The relative uncertainty reduction values are 50%, 46.56%, 51.29%, and 50.00%, respectively.

Case 2 Introduction
As shown in Figure 19, the top event of the fault tree is the leakage of liquid ammonia, which is divided into two layers.There are eight basic events: Y1-poor sealing at the connection point; Y2-poor flange sealing; Y3-poor valve sealing; Y4-poor sealing of the head manhole; Y5-corrosion perforation; Y6-damaged safety valve spring; Y7-incorrect selection of safety valve; and Y8-gasket was torn.

Expert Opinions on Case 2
Six experts evaluated the possibility of the basic events occurring, and the results are shown in Tables 10 and 11.  2, 10, and 11, the failure probability of the basic events is obtained, as shown in Table 12.
Table 12.The failure probability of the basic events.The normal distribution test results of Case 2 are shown in Table 13.It can be seen in Table 13 that the probability of Yi (i = 1~8) obeys the normal distribution.

Case 2: Drawing of Membership Curve
Using the new algorithms and traditional approximate calculation methods, a membership curve for the probability of the top event is drawn considering the following situations: all events occur; except for Yi (i = 1~8), all events occur; Y1, Y2, Y3, and Y4 occur simultaneously; and Y6 and Y7 occur simultaneously.According to the membership curve, mTXe, mTXie (i = 1~8), mTXM1e, and mTXM2e are determined, as shown in Figures 20-30.
Subfigure (a) is a graph of the membership curves obtained by simulation using the new algorithm; subfigure (b) is a graph of the membership curves obtained by the traditional approximation calculation; and subfigure (c) is a comparison of the results obtained by the two calculation methods.The meanings of the letters in the subfigures refer to the explanations in Figure A1 14.In Table 14, it can be seen that, when intermediate event YM2 occurs, the fuzzy importance obtained using the two calculation methods is completely consistent.This is because there are only two basic events that occur simultaneously, and the probability of occurrence is relatively small.
By comparing the reduction values of fuzzy importance uncertainty between the top event and the middle event, it can be seen that the new algorithm can reduce the uncertainty of propagation.Based on the data in Table 14, a fuzzy median uncertainty reduction effect diagram is drawn, as shown in Figure 31.In Figure 32, Y i ′ indicates that basic event Yi does not occur.It can be seen that the fuzzy median value obtained using the improved accurate calculation and the fuzzy median value obtained using the approximate calculation have the same change trend.When Y1 does not occur, the difference between the approximate calculation and the accurate calculation is the largest.
From the intermediate layer to the top layer of the fault tree, as the number of layers increases, the uncertainty reduction value significantly increases.In actual production practice, fault trees are often quite complex.Using new algorithms to handle complex fault trees will have a more significant effect on reducing uncertainty.
As for Case 1, the approximate and accurate fuzzy importance values were calculated, as shown in Table 15.According to Table 15, a fuzzy importance comparison chart was drawn, as shown in Figure 32.In Table 15 and Figure 32, it can be seen that the fuzzy importance of Yi (i = 1-8) significantly decreases.The changing trends of fuzzy importance obtained using the two calculation methods are consistent.The fuzzy importance relative uncertainty reduction effect of Y8 is the most obvious, with the relative uncertainty reduced by 10%.

Changing Rules of Uncertainty Reduction Effect
To further investigate the factors that affect the effectiveness of uncertainty reduction, Figures 33-35 are drawn for an analysis.By comparing Figures 33 and 34, it can be observed that the trend of the change in the mean of basic events is basically consistent with the trend of the change in the absolute uncertainty of fuzzy importance.The larger the mean of basic events, the more significant the reduction effect of fuzzy importance uncertainty after using optimization algorithms.The results of the risk assessments for electrostatic accumulation and the leakage of liquid ammonia conform to this pattern.In Figure 35, it can be seen that, with electrostatic accumulation as the top event, the relative fuzzy importance of the basic event decreases by a maximum of 51.29%.With the leakage of liquid ammonia as the top event, the relative fuzzy importance of the basic event decreases by a maximum of 14.95%.
When the top event is electrostatic accumulation, compared with the top event being the leakage of liquid ammonia, the relative uncertainty reduction effect of the fuzzy importance of the basic event is obvious.This is because the fault tree with electrostatic accumulation as the top event has a larger overall mean.Therefore, in different fault trees, the greater the overall mean value of basic events, the more obvious the effect of reducing fuzzy importance.
It can be seen that the uncertainty reduction effect of low-probability events in fault trees with a large overall mean is often greater than that of high-probability events in fault trees with a small overall mean.

Conclusions
This study proposes a new algorithm that improves traditional normal fuzzy fault tree risk assessment results.First, the algorithm focuses on converting the precise membership curve into multiple cut-set intervals.Second, it further converts them into ordered arrays.Finally, an accurate membership graph is drawn based on this.Through this processing, the limitation where the number of basic events affects the risk assessment results of normal fuzzy fault tree is overcome.More importantly, shortcomings such as difficulty in obtaining accurate calculation results for complex systems and uncertainty in approximate calculations are solved.By combining two chemical accident cases for a verification analysis, the following conclusions are drawn: The new algorithm can significantly reduce the uncertainty of process industry risk assessment.The maximum relative uncertainty reduction in the fuzzy importance of basic events is 50.29%.Considering that accidents in production practice often cause a large amount of property damage, the reliability of the process industry risk assessment results will be greatly improved if the improved results are applied to engineering practice.
This study found that the new algorithm can play a greater role in the following three situations: first, when the average probability of basic events in a single fault tree is high; second, when the overall mean and variance of basic events in different fault trees are large; and, third, when the complexity of the fault tree increases and the number of layers increases.The new algorithm has good adaptability, and it can maintain stable operation and expected performance when the model parameters change, such as the number of basic events and the mean and variance of the probability of basic events.
The contribution of this study is that it can provide a high-precision underlying algorithm for process industry risk assessment, and it improves the reliability of risk assessment results for complex systems.Given the widespread application of fuzzy fault trees in the field of process industry risk assessment, the accurate calculation algorithm of normal fuzzy fault trees proposed in this paper, combined with other intelligent risk assessment algorithms such as neural networks and Bayesian analysis, is a direction worthy of in-depth research.

Figure 1 .
Figure 1.Accurate calculation algorithm for normal fuzzy fault tree (only including OR gates).

Figure 2 .
Figure 2. Accurate calculation algorithm for normal fuzzy fault tree (inclusion of OR gate and AND gate).

Figure 5 .
Figure 5.All basic events occurred (the top event is XTB).

Figure 6 .
Figure 6.Except for X1, all basic events occurred (the top event is XTB).

Figure 7 .
Figure 7. Except for X2, all basic events occurred (the top event is XTB).

Figure 8 .
Figure 8. Except for X3, all basic events occurred (the top event is XTB).

Figure 9 .
Figure 9. Except for X4, all basic events occurred (the top event is XTB).

Figure 10 .
Figure 10.Except for X5, all basic events occurred (the top event is XTB).

Figure 11 .
Figure 11.Except for X6, all basic events occurred (the top event is XTB).

Figure 12 .
Figure 12.Except for X7, all basic events occurred (the top event is XTB).

Figure 13 .
Figure 13.Except for X8, all basic events occurred (the top event is XTB).

Figure 16 .
Figure 16.All basic events occurred (the top event is XTA).

Figure 19 .
Figure 19.Fault tree of leakage of liquid ammonia.

Figure 20 .
Figure 20.All basic events occurred (the top event is YT).

Figure 21 .
Figure 21.Except for Y1, all basic events occurred (the top event is YT).

Figure 22 .
Figure 22.Except for Y2, all basic events occurred (the top event is YT).

Figure 23 .
Figure 23.Except for Y3, all basic events occurred (the top event is YT).

Figure 24 .
Figure 24.Except for Y4, all basic events occurred (the top event is YT).

Figure 25 .
Figure 25.Except for Y5, all basic events occurred (the top event is YT).

Figure 26 .
Figure 26.Except for Y6, all basic events occurred (the top event is YT).

Figure 27 .
Figure 27.Except for Y7, all basic events occurred (the top event is YT).

Figure 28 .
Figure 28.Except for Y8, all basic events occurred (the top event is YT).

3. 3 . 6 .
Case 2: Discussion and Analysis According to Figures 20-30, the fuzzy median values of the improved accurate calculation and traditional approximate calculation are obtained, as shown in Table

Figure 33 .
Figure 33.The relationship between the failure probability of basic events and the reduction in absolute uncertainty (Case 1).

Figure 34 .
Figure 34.The relationship between the failure probability of basic events and the reduction in absolute uncertainty (Case 2).

Figure 35 .
Figure 35.Comparison of uncertainty reduction effects between Case 1 and Case 2.

Table 3 .
The value and meaning of α.

Table 6 .
The failure probability of basic events.

Table 7 .
K-S calculation results of Case 1.

Table 13 .
Calculation results of Case 2.