Failure Occurrence Probability Assessment of Claus Furnace Package in Sulfur Recovery Unit Using Fault Tree Analysis Based on Conventional and Fuzzy-Based Approach

authors:

avatar Abbas Bakbaki 1 , avatar Nader Nabhani 1 , avatar Bagher Anvaripour 1 , * , avatar Gholamabbas Shirali 2

Department of Safety and Protection Engineering, Abadan Faculty of Petroleum Engineering, Petroleum University of Technology, Abadan, IR Iran
Department of Occupational Health Engineering, School of Health, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran

how to cite: Bakbaki A, Nabhani N, Anvaripour B, Shirali G. Failure Occurrence Probability Assessment of Claus Furnace Package in Sulfur Recovery Unit Using Fault Tree Analysis Based on Conventional and Fuzzy-Based Approach. Jundishapur J Health Sci. 2018;10(2):e58162. https://doi.org/10.5812/jjhs.58162.

Abstract

Background:

If shutdown scenario of burner of the sulfur recovery unit takes place, toxic release, fire, and explosion accident can easily occur. Therefore, it is essential to assess the basic causes of burner shutdown scenario. Fault Tree Analysis (FTA) could be used to assess the occurrence probabilities of burner shutdown scenario and its basic causes/events. The failure occurrence probability of these Basic Events (BE) are often not available in lack of data and uncertain conditions.

Objectives:

This study was done to provide a comprehensive approach for analyzing and calculating BE failure occurrence probability affecting the shutdown scenario, using combined fuzzy logic, expert judgment, and FTA.

Methods:

The study was carried out from June to December, 2016. In this study, a fuzzy-based approach based on expert judgment was proposed to calculate the occurrence probability of burner shutdown scenario in lack of data and uncertain conditions. The brainstorming and FMEA, and HAZOP study were first used to identify fault events of the fault tree. Then, based on these methods, the fault tree was constructed. Subsequently, the failure occurrence probabilities of BE and shutdown scenario were calculated using the fuzzy-based approach and a conventional approach. Finally, the Fussell-Vesely importance analysis was used to rank the BE in FTA.

Results:

Results showed that the occurrence probability of shutdown scenario was 4.76E-04 per year. Since the failure occurrence probabilities of some BEs were not available, using failure probability functions in the conventional approach cannot provide failure occurrence probabilities of those BE. Therefore, the occurrence probability of shutdown scenario based on the conventional approach was not available. Based on the Fussell-Vesely importance measure analysis, it was determined that the blower failure while running, air pre-heater blockage, and shut-off valve fail close were 3 major causes of “burner shutdown scenario”.

Conclusions:

The fuzzy-based approach could derive a failure occurrence probability of BE based on expert subjective judgments using the failure possibility distributions (FPDs) in lack of data conditions. This study overcame the weaknesses of the conventional approach, to calculate BE failure occurrence probabilities via fuzzy logic.

1. Background

Oil refineries, especially sulfur recovery units (SRU), are exposed to various types of hazardous events and risks in daily operations, including process risks, risks due to human error (HE), toxicity and mechanical failures, material release, and fire and explosion risks. Thus, it is very important to identify hazardous events and risks, perform risk assessments, and take proper interventions to mitigate hazardous events and risks (1, 2).

Risk assessment with FTA is a method for calculating the occurrence probability of hazardous events or system failure. The FTA is an up-down and deductive approach, by which an undesirable hazardous event, referred to as the top event (TE), is graphically decomposed to possible intermediate and basic causes in increasing detail to determine the causes or combinations of causes of TE. Fault tree provides a graphical approach for determining how failures could occur both qualitatively and quantitatively (3).

Given that the precise and relevant failure rate of BEs, normally required for the quantitative phase of conventional FTA, are not available for process industries, expert judgments are often suggested to assess failure occurrence probability of BEs. Unfortunately, expert judgments are afflicted with qualitative and subjective uncertainties (4, 5). The concept of fuzzy set or fuzzy logic was first proposed by Zadeh (1965) to solve the subjective and qualitative uncertainties of real engineering problems. Fuzzy logic uses possibility distributions to represent expert judgments (6-8).

Fuzzy sets and fuzzy membership functions have been applied for calculating the occurrence probability of hazardous events where relevant failure data of BE are not available. The FTA using fuzzy membership functions was first suggested by Tanaka et al. (1983), as a comprehensive method to solve problems, where exact failure occurrence probabilities of BE are not available (9). Another study of FTA, based on fuzzy sets, was performed by Shi et al. (2014). They suggested an algorithm using Fuzzy FTA (FFTA) and Fuzzy Analytical Hierarchy Process (FAHP) in order to perform risk assessment during a fire and explosion accident for steel oil storage tanks in petroleum industries (5). Rajakarunakaran et al. (2015) showed that expert judgments in the form of qualitative FPDs are appropriate for FTA when relevant failure rate data are not available (6). In addition to the above studies, Mohsendokht (2017) established a model to assess the uranium hexafluoride release from a uranium conversion plant by using fuzzy set theory, possibility concept, and the conventional approach (4). All the studies mentioned above, represent the ability of expert subjective judgment based on fuzzy set theory for solving the lack of data problem in risk and safety assessment studies.

The application of fuzzy-based approach overcomes the limitation of the risk assessment by conventional FTA of the Claus furnace package of SRU by taking expert subjective judgments to assess process failures occurrence probabilities.

2. Objectives

This research aimed at assessing failure occurrence probability of hazardous events based on fuzzy logic, expert judgment, and FTA in lack of data conditions.

3. Methods

The study was carried out at a petroleum refinery during June and December 2016. In the current study, FTA according to fuzzy-based and conventional approaches was applied to assess the failure of burner shutdown in the SRU. The overall framework of the study is shown in Figure 1.

Proposed Framework Based on Fuzzy-Based and Conventional Approaches
Proposed Framework Based on Fuzzy-Based and Conventional Approaches

3.1. Fault Tree Construction Based on Risk Identification

In this step, a critical hazardous event in the SRU is first initiated, called the TE of the fault tree. To construct the fault tree, it is essential to identify hazardous fault events i.e., BEs, intermediate events (IEs), and TE. The logic gates, such as “AND” and “OR” gates, could be used to connect these fault events in the fault tree diagram. The “OR” gate shows that the output event (i.e. TE or IE) occurs if any of the input events (i.e. BEs or IEs) occur. The “AND” gate shows that the output event (i.e. TE or IE) occurs only when all the input events (i.e. BEs or IEs) occur at the same time. The goal of risk identification is to enable identification of all hazardous fault events associated with SRU at different hazardous event levels, e.g. from TE level to IE and BE levels with a view to construct a fault tree diagram. Moreover, construction of a fault tree helps find its cut sets (CSs) while its quantitative analysis calculates the failure occurrence probability of TE using its BEs’ failure occurrence probabilities. A cut set is a set of BEs, whose simultaneous existence involves the occurrence of the TE, and which is independent of the occurrence of the other BEs of the fault tree (3, 6).

3.2. Failure Probabilities of BEs Calculation Based on the Conventional Approach

In the conventional approach, the failure rates of components (or BE) are required to calculate the failure occurrence probabilities of BEs. In this approach, failure probability functions can be used to transform the failure rates of BEs to their failure occurrence probabilities. There are 4 major types of component failure probability functions for calculating the failure occurrence probabilities of BE, which include constant failure probability (failure per demand), non-repairable components, repairable components with revealed failure mode, and repairable components with the periodically proof test (3, 4). For non-repairable components, the failure occurrence probability of BEs can be calculated using Equation 1.

Where λ is the failure rate of BE (number of failures per hour); t is mission time.

Moreover, the BE failure occurrence probability of repairable components could be calculated using Equation 2.

Where μ is the repair rate of BE (number of repairs per hour).

3.3. Failure Probabilities of BEs Calculation Through the Fuzzy Logic Approach

The steps of fuzzy logic approach or fuzzy-based approach are described as follows.

Step 1: Expert weighting evaluation

The first step of the fuzzy-based approach is to select experts and expert weighting evaluation. To perform the FTA based on expert judgments, a group of experts based on Cooke et al. (2008) indicators need to be assembled and selected (10). In some cases, the credibility of experts based on various criteria, such as title (professional position), age, experience level, and educational qualification, is widely different. Therefore, this study considered the relative weighting factor of each expert from 0 to 1. The expert criteria and their corresponding scores are shown in Table 1. For expert j, the relative weighting factor could be calculated using Equation 3 (6, 11):

Table 1.

Expert Criteria and Weighting Scores

RowConstitutionClassificationScore
1TitleChief Engineer5
Supervisor4
Engineer3
Technician2
Operator1
2ExperienceGreater than 305
20 - 304
10 - 193
6 - 92
Less than 61
3Educational qualificationPHD5
Master4
Bachelor3
HND2
School level1
4AgeGreater than 504
4
40 - 503
30 - 392
Less than 301

Where pj is the sum of scores of expert j and n is the number of experts.

Step 2: Select the FPD and their corresponding fuzzy membership functions

In this study, the FPD was used to describe the relationship between an expert subjective judgment (linguistic variable) and its fuzzy membership function value, which is applied in the expert subjective evaluation process (4, 5). This FPD is used to represent failure possibilities of BEs. Based on Miller’s study, the proposed FPD includes 6 linguistic variables, namely, very low (VL), low (L), fairly low (FL), fairly high (FH), High (H), and very high (VH). Based on Chen-Hwang scales (12, 13), the fuzzy triangular membership functions μà (r1, r2, r3) for VL, L, FL, FH, H, and VH are presented in Equations 4 to 9.

Step 3: Expert subjective evaluation process

The purpose of this evaluation was to gather the qualitative failure possibility of a BE from a group of experts using interviews. Using the linguistic variables defined in Equations 4 to 9, each expert in the risk assessment team provided their subjective judgment of the failure possibility of each selected BE.

Step 4: Aggregating the fuzzy failure possibilities obtained

This step was performed to aggregate different failure possibilities obtained from experts. The sum-product method could be used to aggregate different fuzzy failure possibilities of a BE and it is expressed by the following Equation (5):

Where μAggregated (x) is the aggregated fuzzy failure possibilities score (AFFPS), νj represents a relative weighting factor of expert j and μij is the fuzzy triangular membership function i assigned by expert j. k is the number of fuzzy membership functions and n is the number of experts.

Step 5: De-fuzzifying the aggregated fuzzy failure possibilities obtained

The next step of fuzzy logic approach is to de-fuzzify aggregated fuzzy failure possibilities score (AFFPS) to Crisp Failure Possibilities Score (CFPS). Area De-fuzzification technique (ADT) method is the defuzzification process of the AFFPS to a CFPS. The ADT method for a fuzzy triangular membership function μA ̃(r1, r2, r3) is described as follows (14):

Where x* is a defuzzified output equal to CFPS.

Step 6: Converting the crisp failure possibilities to failure probabilities of BE

To generate failure occurrence probabilities of BE, the CFPS values must be converted to Failure Probability (FP) values. The conversion function could be done by Onisawa logarithmic function, which is described as follows (11, 14):

3.4. Calculating the Occurrence Probability of TE

In the current study, based on the failure occurrence probability of BEs, the occurrence probability of a CS was calculated using Equation 13. Also, Equation 14 was used to calculate the occurrence probability of the TE (5, 6).

Where Q(TE) is the occurrence probability of the TE, Q(CSi) is the occurrence probability of cut set i, n is the number of CSs, and m is the number of BEs in the largest CS.

3.5. BEs and CSs Importance Measure Analysis

In this study, importance measure (IM) analysis was performed to assess the contribution of each BE to the TE occurrence and to identify the critical weak points of the burner H-901. The importance measure of BEs can be calculated using the Fussell-Vesely (F-V) equation (3, 4):

Where Q(CSi)j is the occurrence probability of the cut set i containing basic event j, Q(TE\BEj = 0) is the occurrence probability of TE by setting the occurrence probability of basic event j to 0, and nj is the total number of the CSs containing basic event j.

4. Results

4.1. Case Study Description

An SRU recovers sulfur that has to be removed from acid gas (hydrogen sulfide) in oil and gas refineries. In terms of sulfur removal efficiency and safety, the burner of the Claus furnace package is highly critical and important in the SRU. The duty of this stage is to add the quantity of oxygen (combustion process air) required for the Claus reaction and to remove the heat resulting from the reaction between the oxygen and the combustibles in the sour/acid gas. In the operation phase, the Fuel Gas (FLG) flow is needed for continuous firing operation in burner H-901. When the temperature is over the minimum set point (about 900°C), the process Acid Gas (ACG) can flow to the combustion chamber. An air blower is provided for combustion process air (PRA).

4.2. Fault Tree Construction

Based on the brainstorming approach, FMEA, and HAZOP, some important undesired TEs were identified, such as burner H-901 shutdown, gas release, fire, and explosion. In this research, there were 3 major reasons for selecting the shutdown scenario of the burner H-901 due to no flow of feeds (feed = process air + acid gas + fuel gas):

If shutdown scenario of the burner happens, toxic release, fire, and explosion accidents can easily occur.

The various existing shutdown factors for burner of the Claus furnace package process.

The destruction efficiency of the sour/acid gas is more than 65%.

At this stage, according to HAZOP and FMEA, the fault tree was constructed. The fault tree of the burner H-901 shutdown was constructed and is shown in Figure 2.

Fault Tree Structure of Burner H-901 Shutdown Failure
Fault Tree Structure of Burner H-901 Shutdown Failure

In qualitative analysis of FTA, the CSs of the fault tree for burner H-901 shutdown scenario were obtained by applying the Top-Down algorithm. In quantitative analysis, the BEs and CSs probabilities are very important for calculating the TE occurrence probability.

4.3. Calculating Failure Occurrence Probabilities of BEs Through the Fuzzy Logic Approach

Due to lack of failure data, an expert investigation for the failure occurrence probability of BEs of the burner H-901 shutdown was performed by interviewing a total of 4 experts, including a supervisor, chief engineer, process design, and operation personnel of SRU. First, using Equations 3, the relative weighting factors of 4 experts based on Table 1 was calculated and presented in Table 2. Next, using Equations 4 to 9 and interviews with experts, the failure possibilities for all BE were obtained in Table 3. Then, the AFFPSs for 32 BE were calculated using Equation 10 and presented in Table 3. Thereafter, the AFFPSs for all BEs in Table 3 were defuzzified using Equation 11 and presented in Table 4. Finally, using the Onisawa function in Equation 12, the failure occurrence probability for all BEs were calculated and presented in Table 4.

Table 2.

Relative Weighting Factors for 4 Experts

ExpertsTitleEducationExperienceAgeScoresWeight
E1SupervisorBachelor20 - 3040 - 50140.291667
E2Chief engineerMaster10 - 1930 - 39140.291667
E3TechnicianMaster6 - 930 - 39100.208333
E4EngineerMasterLess than 630 - 39100.208333
Sum = 48Sum = 1
Table 3.

Expert Subjective Evaluation and Fuzzy Aggregation Process

BEsFailure Mode DescriptionJudgments of ExpertsAFFPS
X1No flow of fuel gas from upstreamVLLLVL(0.00,0.10,0.30)
X2PV-09717 fail closedFLFLHL(0.24,0.44,0.64)
X3FLG-4 pipe blockageLLLL(0.00,0.20,0.40)
X4FLG-2 pipe blockageLFHFLL(0.16,0.36,0.56)
X5Demister on V-908 blockage due to foulingLFLFLL(0.10,0.30,0.50)
X6XV-09105 fail closedFLFLHL(0.24,0.44,0.64)
X7XV-09107 fail closedFLFLHL(0.24,0.44,0.64)
X8FV-09115 fail closedFLFLHL(0.24,0.44,0.64)
X9Check valve on FLG-2 pipe pluggingLLLL(0.00,0.20,0.40)
X10Check valve on FLG-4 pipe pluggingLLLL(0.00,0.20,0.40)
X11Manual valve on FLG-2 pipe closure (HE)VLLFLFL(0.08,0.23,0.43)
X12Manual valve on FLG-2 plugging/blockageLLLL(0.00,0.20,0.40)
X13No flow of acid gas from upstreamFHHVHFH(0.54,0.74,0.90)
X14ACG-10 pipe BlockageLLFLVL(0.04,0.20,0.40)
X15ACG-12 pipe BlockageLLFLVL(0.04,0.20,0.40)
X16Demister on V-902 blockage due to foulingFHFHFLFL(0.32,0.52,0.72)
X17FV-09111 fail closedFHFHHL(0.36,0.56,0.76)
X18XV-09114 fail closedFHFHHL(0.36,0.56,0.76)
X19Manual valve on ACG-12 pipe closure (HE)VLLLL(0.00,0.14,0.34)
X20Manual valve on ACG-12 plugging/blockageLLLVL(0.00,0.15,0.36)
X21Air pre-heater (E-905) blockageVLLFLVL(0.04,0.14,0.34)
X22XV-09101 fail closedFLFLHL(0.24,0.44,0.64)
X23FV-09113 fail closedFLFLHL(0.24,0.44,0.64)
X24FV-09114 fail closedFLFLHL(0.24,0.44,0.64)
X25PRA-20 pipe redundant section blockageVLLVLVL(0.00,0.06,0.26)
X26PRA-8 pipe blockageLFLLL(0.06,0.26,0.46)
X27Blower C-901 fail while runningHHVHH(0.64,0.84,1.00)
X28PRA-20 pipe blockageVLLVLVL(0.00,0.06,0.26)
X29PRA-24 pipe blockageVLLVLVL(0.00,0.06,0.26)
X30Check valve on PRA-20 pipe pluggingVLLVLVL(0.00,0.06,0.26)
X31Manual valve on PRA-20 pipe closure (HE)VLVLVLVL(0.00,0.00,0.20)
X32Manual valve on PRA-20 plugging/blockageVLLVLVL(0.00,0.06,0.26)
Table 4.

BEs Failure Occurrence Probability and Their F-V Importance Measures

BEsCFPSFuzzy (FP)F-V IMRankConventional (FP)
X10.0227.54E-097.93E-09183.04E-06
X20.1142.76E-052.90E-0555.96E-05
X30.0338.54E-088.98E-08154.96E-09
X40.0868.78E-069.24E-0681.34E-08
X50.0672.85E-063.00E-06101.86E-05
X60.1142.76E-052.90E-0555.76E-05
X70.1142.76E-052.90E-0555.76E-05
X80.1142.76E-052.90E-0553.90E-05
X90.0338.54E-088.98E-08153.86E-05
X100.0338.54E-088.98E-08153.86E-05
X110.0551.12E-061.18E-06134.00E-06
X120.0338.54E-088.98E-08151.20E-05
X130.2122.71E-047.03E-054-
X140.0533.21E-078.34E-08162.67E-09
X150.0533.21E-078.34E-08161.11E-09
X160.1255.93E-051.54E-0571.86E-05
X170.1538.46E-052.20E-0561.29E-04
X180.1538.46E-052.20E-0566.97E-05
X190.0272.43E-086.32E-09204.00E-06
X200.0293.62E-089.39E-09191.20E-05
X210.0361.33E-072.80E-0434.99E-02
X220.1142.76E-055.80E-0223.92E-02
X230.1142.76E-051.66E-06112.67E-05
X240.1142.76E-051.60E-06122.67E-05
X250.0181.61E-099.64E-11211.34E-09
X260.0539.46E-075.49E-08171.77E-09
X270.2454.48E-049.42E-0117.96E-02
X280.0181.61E-093.38E-0692.60E-06
X290.0181.61E-093.38E-0692.18E-07
X300.0181.61E-093.38E-069-
X310.0115.35E-111.12E-0714-
X320.0181.61E-093.38E-069-

4.4. Calculating Failure Occurrence Probabilities of BEs Through the Conventional Approach

OREDA (2002) and Lees (2004) are the main two sources used for assessment of the reliability of data for the conventional approach (3, 15). In the conventional approach, the failure occurrence probabilities of all BEs based on Equations 1 were calculated and presented in Table 4. Due to lack of failure data, the failure occurrence probability for some of the BE could not be calculated, such as the BEs X13, X30, X31, and X32.

4.5. Calculating Occurrence Probability of CSs and TE

Based on the fuzzy-based approach, the occurrence probabilities for 108 CSs were calculated using Equation 13. By assigning the occurrence probabilities of CSs, the occurrence probability of TE was 4.76 E-04 per year, which was calculated using Equation 14. In the conventional approach, it is not possible to calculate the TE and CSs occurrence probabilities, because of the lack of BEs failure data.

4.6. BEs importance Analysis

Equation 15 has been applied to rank the BE based on the F-V importance measure. Table 4 represents the F-V importance analysis results for 32 BE of the fault tree. Basic events F-V importance analysis revealed that process failures, such as X27 (blower C-901 fail while running), X22 (shut-off valve XV-09101 fail close), and X21 (air pre-heater E-905 blockage) play a very important role in burner shutdown occurrence.

5. Discussion

Results from Table 4 showed that failure probabilities of the 4 BE were not available in the conventional approach. These results are consistent with the findings of Shi et al. (2014) and Mohsendokht (2017) on the weakness of the conventional approach in calculating some of the BE in lack of data situations (4, 5). Moreover, the failure occurrence probabilities for some of the BE, such as the X3, X4, X14, X15, X25, X26, X28, and X29, were afflicted with high uncertainty levels. The main causes of this uncertainty were application of irrelevant and low confidence failure data to calculate the failure occurrence probability of BEs. Based on the findings mentioned above, the proposed fuzzy-based approach was required to calculate the failure occurrence probability of BE in this study.

Results of the studies of Omidvari et al. (2014) and Mohsendokht (2017) showed that the process HAZOP was useful for identifying process failures in the fault tree construction phase. However, it cannot be used to identify some of the hazardous events, such as human errors and some of the mechanical failures (4, 11). Therefore, the current study developed a combined HAZOP, brainstorming, and FMEA model to overcome the weakness of process HAZOP to identify the different hazardous events associated with TE. The results revealed that the proposed model is able to identify mechanical and process failures and human errors.

It was observed from importance analysis that process failures were ranked first in contributing to the TE occurrence in the current study, while results in Mohsendokht (2017) showed that a set of process failures and human errors caused TE occurrence (4). Results from Omidvari et al.’s study (2014) showed that a type of process failure, i.e. the obstruction of catalyst had the highest importance measure in fault tree (11). All findings mentioned above, revealed that process failures play a very important role in TEs occurrence for process hazardous scenarios.

Komal et al. (2015) proposed a new fuzzy approach to calculate the failure occurrence probability of BEs based on failure probability function fuzzy parameters in the reliability and risk studies. In their approach, Komal et al. (2015) developed a fuzzy method for computing various failure probability function fuzzy parameters, namely failure rate, repair rate, mean time to repair (MTTR), and expected number of failures (ENOF), for any process industry (16). In this study, fuzzy parameters were used to calculate the failure occurrence probabilities of BEs using the fuzzy lambda-tau method, while FPDs were used to calculate the BEs failure occurrence probability in the current research. Based on the findings mentioned above, in order to calculate the failure occurrence probability of BEs, a fuzzy-based approach using fuzzy failure and repair rate could be used in future studies.

6. Conclusions

In this study, lack of failure data for quantitative phase of conventional FTA was a major issue in the Claus section of SRU. To overcome lack of data, the fuzzy-based approach was applied to calculate the failure occurrence probability of BEs in this study. The main advantage of using FPDs in this study was that the proposed approach could directly express failure occurrence probability of each BE. Also, the results revealed that the proposed approach could eliminate the limitation of conventional FTA by using fuzzy logic.

Acknowledgements

References

  • 1.

    Barua S, Gao X, Pasman H, Mannan MS. Bayesian network based dynamic operational risk assessment. J Loss Prevent Process Indust. 2016;41:399-410. https://doi.org/10.1016/j.jlp.2015.11.024.

  • 2.

    Bakbaki A, Nabhani N, Anvaripour B, Shirali G. Probabilistic risk assessment using fuzzy fault tree analysis based on two types of failure possibility distributions in process industries. J Occup Hyg Eng. 2017;4(2):41-52.

  • 3.

    Lees F. Lees' Loss prevention in the process industries: Hazard identification, assessment and control. Butterworth-Heinemann; 2012.

  • 4.

    Mohsendokht M. Risk assessment of uranium hexafluoride release from a uranium conversion facility by using a fuzzy approach. J Loss Prevent Process Indust. 2017;45:217-28. https://doi.org/10.1016/j.jlp.2017.01.004.

  • 5.

    Shi L, Shuai J, Xu K. Fuzzy fault tree assessment based on improved AHP for fire and explosion accidents for steel oil storage tanks. J Hazard Mater. 2014;278:529-38. [PubMed ID: 25010458]. https://doi.org/10.1016/j.jhazmat.2014.06.034.

  • 6.

    Rajakarunakaran S, Maniram Kumar A, Arumuga Prabhu V. Applications of fuzzy faulty tree analysis and expert elicitation for evaluation of risks in LPG refuelling station. J Loss Prevent Process Indust. 2015;33:109-23. https://doi.org/10.1016/j.jlp.2014.11.016.

  • 7.

    Zadeh LA. Fuzzy sets. Inf Control. 1965;8(3):338-53. https://doi.org/10.1016/s0019-9958(65)90241-x.

  • 8.

    Sa'idi E, Anvaripour B, Jaderi F, Nabhani N. Fuzzy risk modeling of process operations in the oil and gas refineries. J Loss Prevent Process Indust. 2014;30:63-73. https://doi.org/10.1016/j.jlp.2014.04.002.

  • 9.

    Tanaka H, Fan LT, Lai FS, Toguchi K. Fault-Tree Analysis by Fuzzy Probability. IEEE Transact Reliabil. 1983;R-32(5):453-7. https://doi.org/10.1109/tr.1983.5221727.

  • 10.

    Cooke RM, ElSaadany S, Huang X. On the performance of social network and likelihood-based expert weighting schemes. Reliabil Engin System Safety. 2008;93(5):745-56. https://doi.org/10.1016/j.ress.2007.03.017.

  • 11.

    Omidvari M, Lavasani SMR, Mirza S. Presenting of failure probability assessment pattern by FTA in Fuzzy logic (case study: Distillation tower unit of oil refinery process). J Chem Health Safety. 2014;21(6):14-22. https://doi.org/10.1016/j.jchas.2014.06.003.

  • 12.

    Miller GA. The magical number seven plus or minus two: some limits on our capacity for processing information. Psychol Rev. 1956;63(2):81-97. [PubMed ID: 13310704].

  • 13.

    Chen SJ, Hwang CL, Hwang FP. Fuzzy multiple attribute decision making(methods and applications. Berlin: Springer-Verlag; 2011.

  • 14.

    Purba JH, Lu JIE, Zhang G, Ruan DA. An Area Defuzzification Technique to Assess Nuclear Event Reliability Data from Failure Possibilities. Int J Comput Intelligence Applicat. 2012;11(4):1250022. https://doi.org/10.1142/s1469026812500228.

  • 15.

    SINTEF. OREDA. Offshore reliability data handbook. OREDA; 1984.

  • 16.

    Komal Chang D, Lee S. Fuzzy reliability analysis of dual-fuel steam turbine propulsion system in LNG carriers considering data uncertainty. J Nat Gas Sci Engin. 2015;23:148-64. https://doi.org/10.1016/j.jngse.2015.01.030.