This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Aa
Sa
b
1
Icutntcm
cuttpatFo
h0
Process Safety and Environmental Protection 9 3 ( 2 0 1 5 ) 75–88
Contents lists available at ScienceDirect
Process Safety and Environmental Protection
journa l h om ep age: www.elsev ier .com/ locate /ps ep
n extension to Fuzzy Fault Tree Analysis (FFTA)pplication in petrochemical process industry
Sciences & Research Branch, Tehran Science and Research Branch University, Hesarak, Tehran, IranDepartment of Marine Engineering, Istanbul Technical University, Tuzla, 34940 Istanbul, Turkey
a b s t r a c t
Fault Tree Analysis (FTA) is an established technique in risk management associated with identified hazards specific
to focused fields. It is a comprehensive, structured and logical analysis method aimed at identifying and assessing
hazards of complex systems. To conduct a quantitative FTA, it is essential to have sufficient data. By considering the
fact that sufficient data is not always available, the FTA method can be adopted into the problems under fuzzy envi-
ronment, so called as Fuzzy Fault Tree Analysis (FFTA). This research extends FFTA methodology to petrochemical
process industry in which fire, explosion and toxic gas releases are recognized as potential hazards. Specifically, the
case study focuses on Deethanizer failure in petrochemical plant operations to demonstrate the proposed methodol-
ogy. Consequently, the study has provided theoretical and practical values to challenge with operational data shortage
n conventional FTA, the Failure Probabilities (FP) of systemomponents (i.e. Basic Events (BEs)) are treated as exact val-es. However, for many systems, it is very difficult to estimatehe precise failure rate or probabilities of individual compo-ents or failure events in the quantitative analysis of faultree structures. In other word, the crisp approach has diffi-ulty in conveying imprecision or vagueness nature in systemodeling (Liang and Wang, 1993; Lavasani et al., 2012).To remedy the gap about the mentioned inadequacy of the
onventional FTA, extensive research has been performed bysing fuzzy set theory. The pioneering work on this belongso Tanaka et al. (1983), which treated probabilities of BEs asrapezoidal fuzzy numbers, and applied the fuzzy extensionrinciple to determine the probability of Top Event (TE). Linnd Wang (1997) developed a hybrid method which can simul-aneously deal with probability and possibility measures in aTA. Sawer and Rao (1994) applied �-cuts to determine the FP
f the TE in mechanical systems modeling with Fuzzy Fault
Trees (FFTs). Cai et al. (1991) and Huang et al. (2004) adoptedpossibility theory to analyze FFTs. Dong and Yu (2005) appliedthe hybrid method to analyze FP of oil and gas transmissionpipeline. As another approach, Shu et al. (2006) used intuition-ist fuzzy methods to analyze FT on a printed circuit boardassembly.
Furthermore, Ping et al. (2007) used FFTA for assessing fail-ure of bridge construction. Toward marine accident analysisand prevention, Celik et al. (2010) proposed an investigationmodel based on FTA supported with fuzzy sets. Wang et al.(2013) employed FFTA for fire and explosion of crude oil tanks.Recently, Liu et al. (2014) used FTA in emergency responseplanning.
The main aim of this research is to extend FFTA methodol-ogy to petrochemical process industry. This section introducesthe existing applications of FFTA throughout the variousindustries. The steps of research methodology including iden-tifying BEs, obtaining FP of BEs with known failure rate, rating
. Lavasani).
Crisp Failure Possibility (CFP) of BEs into FP, calculating all
76 Process Safety and Environmental Protection 9 3 ( 2 0 1 5 ) 75–88
prop
Fig. 1 – Structure of
Minimal Cut sets (MCs) and FP of TE, and ranking of MCs areprovided in Section 2. In Section 3, a case study on an impor-tant section of petrochemical plant is demonstrated. The lastsection emphasis the highlights of the research.
2. Research methodology
In circumstances where the lack or incompleteness of dataexists, there is a need to incorporate expert judgment intothe FTA study. A framework proposed based on the fuzzy settheory with the FTA method is capable of quantifying thejudgment from experts who express opinions qualitatively.The new proposed framework is developed in eight differentstages. In the first stage, BEs with known failure rates isseparated from those BEs with a vague failure rate. Thesecond stage is to obtain the FPs of BEs with known failurerates. In the third stage, expert judgments are assignedto the vague BEs. These ratings are fuzzy numbers. Thefourth stage is an aggregation procedure. It is performed byaggregating of experts’ opinions. A defuzzification processwill then be adopted to transform the experts’ judgments(fuzzy possibility) to corresponding crisp possibility valuesby employing an appropriate algorithm. The sixth stage is toconvert crisp possibilities values to the FPs. MCs are identifiedand quantified in the seventh stage. In the last stage, rankingof all MCs can consequently be produced. Fig. 1 presents thestructure of proposed methodology.
2.1. Identifying BEs
As mentioned, the first step of the methodology is to iden-tify hazards with known failure rates from vague hazards.
Occurrence failure rate of some hazards are available from PDS(2010).
osed methodology.
2.2. Obtaining Failure Probability (FP) of BEs withknown failure rate
The foundation of a good analysis is the pedigree of failurerate or event probability data that is assigned to BEs. There arepredominantly three methods that can be used to determinethe occurrence probability of an event namely (Preyssl, 1995):
The statistical method involves the treatment of direct testof experience data and the calculation of the probabilities. Theextrapolation method involves the use of model predictionand similar condition or using standard reliability handbook.The expert judgment method involves direct estimation ofprobabilities by specialists.
A component is tested periodically with test interval. A fail-ure may occur at any time in the test interval, but the failureis only detected in a test. After a test/repair, the componentis assumed to be “as good as new”. This is a typical situationfor many safety-critical components, like sensors and safetyvalves. If an event failure of a kind which can be inspected,the component failure probability can be obtained from Eq. (1)(Spouge, 2000; Rausand and Hoyland, 2004).
P(t) = 12
�� (1)
where � is the component failure rate and � is the inspectioninterval.
Process Safety and Environmental Protection 9 3 ( 2 0 1 5 ) 75–88 77
If a component is of a kind which cannot be inspected.he component failure probability P, which is also called thenreliability, is determined from Eq. (2).
(t) = 1 − e−�t (2)
here � the component failure rate and t is the relevant timenterval. Based on the Maclaren series, the above equation for
can be obtained from Eq. (3) if �t � 0.1
(t) = 1 −(
1 + −�t
1!+ �2t2
2!+ �3t3
3!+ · · · + �ntn
n!
)∼= �t (3)
.3. Rating state
n this stage, experts express their opinions for each BE withespect to each subjective attribute. Expert elicitation is theynthesis of experts’ opinions of a subject where there isncertainty due to insufficient data because of physical con-traints or lack of resources. Experts’ elicitation is essentially
scientific consensus methodology. Expert elicitation gen-rally quantifies uncertainty. Examples of fields that haveontributed to probability elicitation are decision analysis,sychology, risk analysis, Bayesian statistics, mathematicsnd philosophy.
Quantification of subjective probabilities is employed in aumber of circumstances (Korta et al., 1996):
Evidence is incomplete because it cannot be reasonablyobtained.
Data exists only from analogous situations (one might knowthe solubility of one mineral and might use this informationto infer the solubility of another mineral).
There are conflicting models or data sources. Scaling up from experiments to target physical processes
is not direct (scaling of mean values is often much simplerthan rescaling the uncertainties).
Expert knowledge is influenced by individual perspec-ives and goals (Ford and Sterman, 1998). Therefore, completempartiality of expert knowledge is difficult to achieve. Anmportant consideration is the selection of heterogeneousxpert (e.g. both scientists and workers) or homogenous groupf experts (e.g. only scientists). The effect of difference in per-onal experience on expert judgment is assumed to be smallern homogenous group compared to a heterogeneous group. Aeterogeneous group of experts can have an advantage over aomogenous group through considering all possible opinions.
In this study, a heterogeneous group of experts is selectedor evaluating the probability of vague events. The weightingactors of experts are determined according to Table 1.
Rating of expert judgment can be carried out by linguisticerms which are used for soliciting expert opinions for eachE.
.4. Aggregating stage
ince each expert may have a different opinion accordingo his/her experiences and expertise in the relevant field, its necessary to aggregate expert’s opinion to reach a con-ensus. Different types of aggregation can be used (Hsu andhen, 1994; Aqlan and Ali, 2014). Aqlan and Ali (2014) pro-
osed a new method for aggregation of expert judgment bymploying triangle fuzzy numbers. As mentioned, Aqlan and
Ali (2014) model can just aggregate triangle fuzzy numberswhilst Hsu and Chen (1994) model is able to aggregate triangleand trapezoidal fuzzy numbers. Linguistic terms of this paperare mixture of triangle and trapezoidal fuzzy numbers. There-fore, Hsu and Chen (1994) method of aggregation is used inthis paper.
Hsu and Chen (1994) presented an algorithm to aggregatethe linguistic opinions of a homogenous/heterogeneous groupof experts. Suppose each expert, Ek (k = 1, 2, . . ., M) expresseshis/her opinion on a particular attribute against a specific con-text by a predefined set of linguistic variables. The linguisticterms can be converted into corresponding fuzzy numbers.The detailed algorithm is described as follows:
1. Calculate the degree of agreement (degree of similarity)Suv(Ru, Rv) of the opinions between each pair of expertsEu and Ev, where Suv(Ru, Rv). According to this approach,A = (a1, a2, a3, a4) and B = (a1, a2, a3, a4) are two standardtrapezoidal fuzzy numbers. Then the degree of similarityfunction of S, which is defined as:
S(A, B) = 1 − 14
4∑i=1
|ai − bi| (4)
where (A, B) ∈ [0, 1], the larger value of S(A, B), the greatersimilarity between two fuzzy numbers of A and B.
2. Calculate the Average Agreement (AA) degree AA(Eu) of theexperts.
AA(Eu) = 1M − 1
4∑u /= v
v = 1
S(Ru, Rv) (5)
3. Calculate the Relative Agreement (RA) degree, RA(Eu) of theexperts.
78 Process Safety and Environmental Protection 9 3 ( 2 0 1 5 ) 75–88
diag
Fig. 2 – Process flow
where ˇ(0 ≤ ≤ 1) is a relaxation factor of the proposedmethod. It shows the importance W(Eu) over RA(Eu). When
= 0 no importance has been given to the weight of anexpert and hence a homogenous group of experts is used.When = 1, the consensus degree of an expert is the sameas its importance weight. The consensus degree coefficientof each expert is good measure for evaluating the relativeworthiness of each expert’s opinion. It is responsibility ofthe decision maker to assign an appropriate value to ˇ.
5. Finally, the aggregated result of the experts’ judgment RAG,can be obtained as follows:
Defuzzification is the process of producing a quantifiableresult in fuzzy logic. Defuzzification problems emerge fromthe application of fuzzy control to the industrial processes(Zhao and Govind, 1991). The center of area defuzzificationtechnique is selected here. This technique was developed bySugeno in 1985 (Sugeno, 1999). This method can be expressedas:
∫
X∗ = �i(x)xdx∫
�i(x)dx(9)
ram of Deethanizer.
where X* is the defuzzified output, �i(x) is the aggregatedmembership function and x is the output variable. The aboveformula can be shown as follows for triangular and trape-zoidal fuzzy numbers. Defuzzification of fuzzy number A =(a1, a2, a3) is
X∗ =∫ a2
a1
x−a1a2−a1
xdx +∫ a3
a2
a3−xa3−a2
xdx∫ a2
a1
x−a1a2−a1
dx +∫ a3
a2
a3−xa3−a2
dx= 1
3(a1 + a2 + a3) (10)
Defuzzification of trapezoidal fuzzy number A =(a1, a2, a3, a4) can be obtained by Eq. (11).
X∗ =∫ a2
a1
x−a1a2−a1
xdx +∫ a3
a2xdx +
∫ a4
a3
a4−xa4−a3
xdx∫ a2
a1
x−a1a2−a1
dx +∫ a3
a2dx +
∫ a4
a3
a4−xa4−a3
dx
= 13
(a4 + a3)2 − a4a3 − (a1 + a2)2 + a1a2
(a4 + a3 − a1 − a2)(11)
2.6. Transforming CFP of BEs into FP
As aforementioned, there are data available for failure ratesof some events whilst the data associated with the othersare vague. There is inconsistency between FPs of certain BEs
and CFPs of vague events. This issue can be performed byusing Eq. (12). Onisawa (1988) has proposed a function which
Process
Safety
an
d En
viro
nm
enta
l Pro
tection
9
3
( 2
0 1
5 )
75–88
79
Table 2 – BEs and failure states of Deethanizer.
1.Trip of P-401
1.1 Power Failure
1.1.1 Power failure from source(Mobin)
1.1.2 Trip from Sub Station1.1.2.1 Breaker Failure1.1.2.2 Transformer Failure1.1.2.3 Human error to stop thepump
1.2 Equipment Failure
1.2.1 Poor Maintenance
1.2.1.1 Poor PM (planning andcontrol)
1.2.1.2 Bad/wrong installation1.2.1.2.1 Lack of supervision1.2.1.2 .2 Lack of competency ofmaintenance workers
1.2.1.3 Poor quality of equipmentparts
1.2.1.3.1 Procurement inadequacy1.2.1.3.2 Lack of supervision andinspection by asset integrity dept.
1.2.2 Blockage in suction of P-4011.2.2.1 icing due to dryersdeficiency1.2.2.2 Hydrate formation1.2.2.3 chocking of strainer
1.3 Human error to stop the pumpfrom DCS/Push bottom
1.3.1 Lack of knowledge (Training)1.3.2 Lack of skill (experience)1.3.3 Lack ofperception/carelessness
2. Failure FV 40071 A to close
2.1 Failure of Instrument Air (IA)2.1.1 Loss of IA from source(Mobin)2.1.2 Human error to close the IAvalves in plant/offsite battery limit
2.2. Failure/Error of FIC2.2.1 Equipment Failure2.2.2 Loose connections in IAnetwork
2.3 Human error in DCS to closethe FV
3. Failure of E-420
2.1 Failure of Instrument Air (IA)2.1.1 Loss of IA from source(Mobin)2.1.2 Human error to close the IAvalves in plant/offsite battery limit
4. Failure of TV-40075 to open 2.1 Failure of Instrument Air (IA)2.1.1 Loss of IA from source(Mobin)2.1.2 Human error to close the IAvalves in plant/offsite battery limit
80
Process
Safety
an
d En
viro
nm
enta
l Pro
tection
9
3
( 2
0 1
5 )
75–88
Table 2 (Continued)
5. Failure of PV-40102 to close4.1 Failure/Error of TICHuman error in DCS to open the FV2.1 Failure of Instrument Air (IA)Failure/Error of PIC2.3Human error in DCS to close theFVFailure of TV-40104
4.1.1 Equipment Failure2.1 Failure of Instrument Air (IA) 2.1.1 Loss of IA from source
(Mobin)2.1.1 Human error to close the IAvalves in plant/offsite battery limit
6.1.1 Failure/Error of TIC 6.1.1.1 Equipment Failure
6.2 Equipment Failure6.3 chocking due to polymerization
7.1 Failure of external body ofT-402
2.3 Human error in DCS to openthe FV6.2.1 Failure of FV 400916.2.2 6.2.3Failure of C-501Equipment failure due to corrosion
7.1.1 Bad manufacturing (Design)7.1.2 Poor welding
7.1.2.1 Poor Procedure7.1.2.2 Poor Supervision
7.1.3 Vibration7.1.4 corrosion 7.1.4.1 Internal
7.1.4.2 External7.2 Failure of supports 7.2.1 Bad manufacturing
7.2.2 Poor welding7.2.3 Vibration7.2.4 Lack of fire proofing in case offire in adjacent area
Process Safety and Environmental Protection 9 3 ( 2 0 1 5 ) 75–88 81
tatio
cbhta1
F
2
Bitc
P
s
2
OiiiSi
Fig. 3 – Fault tree represen
an be used for converting CFP to FP. This function is derivedy addressing some properties such as the proportionality ofuman sensation to the logarithmic value of a physical quan-
ity. The probability rate can be obtained from possibility rates follows (Onisawa, 1988, 1990, 1996; Onisawa and Nishiwaki,988; Lin and Wang, 1998):
P ={
1/10K, CFP /= 0
0, CFP = 0, K =
[(1 − CFP
CFP
)]1/3
× 2.301 (12)
.7. Calculating all MCs and FP of TE
y definition, a MC is a combination (intersection) of BEs lead-ng to the TE. The combination is a “minimal” combination inhat all failures are needed for the TE to occur. TE probabilityan be obtained from Eq. (13).
Where P(MCi) is the occurrence probability of minimal cutet i.
.8. Ranking of MCs
ne of the most important outputs of an FTA is the set ofmportance measures that are calculated for the TE. Suchmportance measures establish the significance for all the MCsn the FT in terms of their contributions to the TE probability.
ensitivity Analysis (SA) is the study of how the uncertainty
n the output of a mathematical model or system (numerical
n of Deethanizer failure.
or otherwise) can be apportioned to different sources ofuncertainty in its inputs. In other words, SA can be used asmethod for testing robustness of a model results in presenceof uncertainty. SA of FTs is estimated by Importance measurestechniques. The following probabilistic importance measuresare currently in use for FTA:
• Criticality importance• Risk Reduction Worth (RRW) or Top Decrease Sensitivity
(TDS)
SA can be useful for a range of purposes, including:
• Testing the robustness of the results of a model or systemin the presence of uncertainty.
• Increased understanding of the relationships between inputand output variables in a system or model.
• Uncertainty reduction: identifying model inputs that causesignificant uncertainty in the output.
• Searching for errors in the model (by encountering unex-pected relationships between inputs and outputs).
• Model simplification – fixing model inputs that have noeffect on the output, or identifying and removing redundantparts of the model structure.
Fussell-Vesely Importance (F-VI) is the contribution of theMC to the TE probability. F-VIs are determinable for every MCmodeled in the FT. This provides a numerical significance ofall FT elements and allows them to be prioritized. The F-VI iscalculated by summing all the causes (MCs) of the TE involving
the particular event. This measure has been applied to MCs todetermine the importance of individual MC. Where Qi(t) is the
82 Process Safety and Environmental Protection 9 3 ( 2 0 1 5 ) 75–88
Table 3 – Deethanizer BEs.
Deethanizer failure Fault tree ref BE failure rate
1 Power failure from source (Mobin) 1.1.1 Linguistic term2 Breaker failure 1.1.2.1 Failure rate3 Transformer failure 1.1.2.2 Failure rate4 Human error to stop the pump (P-401) 1.1.2.3 Linguistic term5 Poor PM (planning and control) (P-401) 1.2.1.1 Linguistic term6 Lack of supervision (P-401) 1.2.1.2.1 Linguistic term7 Lack of competency of maintenance workers 1.2.1.2.2 Linguistic term8 Procurement inadequacy 1.2.1.3.1 Linguistic term9 Lack of supervision and inspection by asset integrity dept. 1.2.1.3.2 Linguistic term
10 Icing due to dryers deficiency (Suction P-401) 1.2.2.1 Linguistic term11 Hydrate formation (Suction P-401) 1.2.2.2 Linguistic term12 Chocking of strainer (Suction P-401) 1.2.2.3 Linguistic term13 Lack of knowledge (Training) 1.3.1 Linguistic term14 Lack of skill (experience) 1.3.2 Linguistic term15 Lack of perception/carelessness 1.3.3 Linguistic term16 Loss of IA from source (Mobin) 2.1.1 Failure rate17 Human error to close the IA valves in plant/offsite battery limit 2.1.2 Linguistic term18 Equipment Failure (FIC of FV-40071A) 2.2.1 Failure rate19 Loose connections in IA network (FIC of FV-40071A) 2.2.2 Failure rate20 Human error in DCS to close the FV 2.3 Linguistic term21 Failure of C-501 due to poor PM 3.1 Failure rate22 Internal (inside erosion) (E-420) 3.2.1 Failure rate23 External (humidity) (E-420) 3.2.2 Failure rate24 Equipment failure (TIC of TV-40075) 4.1.1 Failure rate25 Loose connections in IA network (TIC of TV-40075) 4.1.2 Failure rate26 Failure/error of PIC (PV-40102) 5.1 Failure rate27 Equipment failure (TIC of TV-40104) 6.1.1.1 Failure rate28 Loose connections in IA network (TIC of TV-40104) 6.1.1.2 Failure rate29 Failure of FV-40091 6.2.1 Failure rate30 Failure of C-501 6.2.2 Failure rate31 Internal corrosion (E-422) 6.2.3.1 Failure rate32 External corrosion (E-422) 6.2.3.2 Failure rate33 Chocking due to polymerization (E-422) 6.3 Failure rate34 Bad manufacturing (design) (external body of T-402) 7.1.1 Failure rate35 Poor welding Procedure (external body of T-402) 7.1.2.1 Linguistic term36 Poor welding Supervision (external body of T-402) 7.1.2.2 Linguistic term37 Vibration (external body of T-402) 7.1.3 Failure rate38 Internal corrosion (external body of T-402) 7.1.4.1 Failure rate39 External corrosion (external body T-402) 7.1.4.2 Failure rate40 Bad manufacturing (supports) 7.2.1 Failure rate41 Poor welding (supports) 7.2.2 Linguistic term42 Vibration (supports) 7.2.3 Failure rate
(supp
43 Lack of fire proofing in case of fire in adjacent area
contribution of MC i to failure of the system, the importancemeasure can be quantified as follows (Modarres, 2006):
IFVi (t) = Qi(t)
QS(t)(14)
Qi(t) = probability of failure of MCi, QS(t) = probability of failureof TE due to all MCs.Risk Reduction Worth (RRW) measures thedecrease in the probability of the TE if a given MC is assurednot to occur. This importance measure can also be called the
Top Decrease Sensitivity (TDS) (Shu et al., 2006). RRW for a MCshows the decrease in the probability of the TE that would be
Table 4 – FP of the BEs with known failure rate.
BEs FP of BEs BEs FP of BEs
1.1.2.1 0.01 3.2.1 0.012
1.1.2.2 0.015 3.2.2 0.008
2.1.1 0.002 4.1.1 0.014
2.2.1 0.014 4.1.2 0.009
2.2.2 0.018 5.1 0.013
3.1 0.02 6.1.1.1 0.014
orts) 7.2.4 Linguistic term
obtained if the MC did not occur. Therefore, the RRW can becalculated by re-quantifying the FT with considering of theprobability of the given MC to 0. It thus measures the maxi-mum reduction in the TE probability. RRW can be used to testmodel robustness.
3. Case study
A petrochemical plant is built to manufacture Medium DensityPolyethylene (MDP) and High Density Polyethylene (HDP) for
the plastics processing industry. The facility has a capacityof 300,000 t/a based on 7920 h/y. The products are marketed
Process Safety and Environmental Protection 9 3 ( 2 0 1 5 ) 75–88 83
Fig. 4 – Chen and Hwang conversion scale.
Table 5 – Fuzzy number of conversion scale 6.
Linguistic terms Fuzzy sets
Very low (VL) (0,0,0.1,0.2)Low (L) (0.1,0.25,0.25,0.4)Medium (M) (0.3,0.5,0.5,0.7)High (H) (0.6,0.75,0.75,0.9)Very high (VH) (0.8,0.9,1,1)
up
SaD
3
ApiDttcpiswtao
Table 7 – Expert judgment on vague BEs.
BEs Expert judgment on vague BEs
E1 E2 E3
1.1.1 L M M1.1.2.3 VL L M1.2.1.1 L L M1.2.1.2.1 M H H1.2.1.2.2 M M M1.2.1.3.1 L L M1.2.1.3.2 M M H1.2.2.1 M M M1.2.2.2 L L L1.2.2.3 M H H1.3.1 M H M1.3.2 VL L M1.3.3 L VL L2.1.2 M L L2.3 L VL L7.1.2.1 L L M7.1.2.2 M H M7.2.2 L L VL7.2.4 L M L
nder the brand name LUPOLEN. For easy comprehension, thelant is divided into four main sections:
Plant Section 01: Purification (Area 700) and dosing of feed-stock (Area 100)Plant Section 02: Polymerization (Area 200) and powder hand-ling (Area 300 & 400)Plant Section 03: Extrusion and product handling (Area 500)Plant section 04: Granulate handling and logistic
Additionally, plant facilities are provided in Area 600.ection 02 is one of the most important plant sections. Deeth-nizer is one of the main equipment of this section. Therefore,eethanizer is selected as case study in this paper.
.1. Identifying BEs of Deethanizer
vent system on propylene side of the E-420 and E-422 isrovided to evacuate non-condensable components contained
n propylene that could accumulate. Bottom temperature ofeethanizer T-402 must be controlled carefully because at a
oo high temperature the heavier olefins in the bottom tendo form gums and polymers, thus fouling the bottom of theolumn and the reboiler. For this reason, a spare reboiler isrovided and when reboiling becomes inefficient due to foul-
ng, operator has to switch from the operating reboiler to thepare one. The reboiling control valve is on the steam inlet linehere it controls the condensing pressure and hence tempera-
ure. Steam is injected in order to get a maximum temperature
t Deethanizer bottom of about 82 ◦C. A piping pot is providedn the steam condensate line to make hydraulic guard.
Table 6 – Experts’ profiles and decision weights.
No of expert Title Servicetime (Year)
Edu
1 Senior 10–19 M2 Senior <5 P3 Junior 6–9 M
The overhead specification for the Deethanizer is thepropylene content in the hydrogenated C2 cut. Normal value is0.3 wt% (0.4 wt% max.). The liquid reflux flow rate is adjustedby FV-40071 to reach this specification. A polymerizationinhibitor is injected by the Polymerization Inhibitor InjectionPackage W-402 in the feed and at bottom of Deethanizer inorder to limit the fouling tendency. Based on Process FlowDiagram (PFD) which is shown in Fig. 2.
Considering the working principles and process flow dia-gram of Deethanizer unit, BEs associated with this section areidentified and listed in Table 2.
Fault tree representation of the Deethanizer failure is illus-trated in Fig. 3.
Considering the Deethanizer process flow diagram andrelated failure states, it is the next issue to determine whetherthe BEs are known or unknown.
S (E12) 0.75 AA (E1) 0.75S (E13) 0.75 AA (E2) 0.875S (E23) 1 AA (E3) 0.875RA (E1) 0.3 CC (E1) 0.3315RA (E2) 0.35 CC (E2) 0.3415RA (E3) 0.35 CC (E3) 0.327Weight of expert 1 (E 1) 0.363Weight of expert 2 (E 2) 0.333Weight of expert 3 (E 3) 0.304Aggregation for 1.1.1 0.2337 0.417125 0.417125 0.60055
3.2. Separating BEs with known failure rate from BEswith unknown failure rate
The elements of the FT are divided into failure probabilityanalysis of BEs with known probabilities of occurrence andsubjective linguistic evaluations of hazards with unknownprobability rate. 43 BEs are identified for Deethanizer failure.24 of them are BEs with known occurrence probabilities whilstthere are not historical data available for the other 19 BEs.Table 3 presents all of the BEs associated with the proposedFT.
3.3. Calculating FPs of BEs with known occurrenceprobabilities
As aforementioned, the foundation of a good analysis isthe pedigree of failure rate or event probability data that isassigned to BEs. Therefore, occurrence probabilities of haz-ard with known failure rate can be estimated by using Eqs.(1)–(3). For example, the rate of mechanical failure of homoge-nizer is 2.4 × 10−3 with 4 inspections in a year. Therefore, FP ofmechanical failure of homogenizer can be obtained by usingEq. (1) as follows:
1 −3 4 −3
FPmechanical failure of homogenizer =2
× 2.4 × 10 ×12
= 4 × 10
Table 9 – Aggregation calculation for each subjective BE.
FP of the BEs with known failure rates are calculated andpresented in Table 4.
3.4. Rating state
In the proposed method, a numerical approximation proposedby Chen and Hwang (1992) is used to convert linguistic term totheir corresponding fuzzy numbers. There are generic verbalterms in the system where scale 1 contains two verbal terms(linguistic terms) and scale 8 contains 13 verbal terms (linguis-tic terms). The typical estimate of a human memory capacity,is seven plus-minus two chunks, which means that the suit-able number for linguistic term selection for human beingsto make an appropriate judgment is between 5 and 9 (Miller,1956; Nicokis and Tsuda, 1985). Therefore, conversion scaleof 6 which contains 5 verbal terms is selected for performingthe subjective assessment of hazards with unknown failurerate. Fig. 4 presents the fuzzy linguistic scale that is used inthis paper to involve the judgments of experts with respect tohazards with unknown failure rate.
The given linguistic terms are in the form of both triangu-lar and trapezoidal fuzzy numbers. Table 5 represents all thefuzzy numbers in the form of trapezoidal numbers.
As mentioned, a heterogeneous group of experts is
employed to perform the judgment for the vague events.
Table 10 – Deffuzification process for all subjective BEs.
he weights of experts are not equal. Experts’ weights cane obtained based on their profiles and competencies. Inhis case, three experts are employed for performing theudgments. Table 6 shows the experts’ profiles and decisioneights.
Expert judgment on the BEs with unknown failure rates arellustrated by Table 7.
.5. Aggregation of BEs
n this stage, all the ratings are aggregated under each subjec-ive BE. As an example, detailed aggregation calculations for BEf 1.1.1 are given in Table 8. (Relaxation factor) is considered.5 in aggregation calculation of subjective BEs.
These calculations contain attribute based aggregation cal-ulation, such as average degree of agreement (AA), relativeegree of agreement of each expert (RA), etc. After the aggre-ation calculations, the results of all the BEs are presented in
The center of area deffuzification technique is employed tocalculate the deffuzification of all the subjective BEs. Table 10shows the result of subjective BEs deffuzification.
3.7. Converting CFP of BEs into FP
CFP of the subjective BEs can be transformed into the corre-
sponding FP by using Equation 12. Table 11 presents FP of allthe subjective BEs.
Fig. 5 – Result of sensitivity
.301607 0.0008 22
3.8. Calculating FP of TE
To quantify the probability of TE of the fault tree a proba-bility for each BE in the fault tree must be provided. TheseBE probabilities are then propagated upward to the TE usingthe Boolean relationships. In other words, conventional FTArules are employed for TE quantification. Therefore, all BEs areconsidered independent. The BE probabilities can be propa-gated upward using MCs. Table 12 presents the FP of all the
MCS. Furthermore, TE is obtained by using Eq. (13). The valueof FP of the TE is 0.3024 per year.
analysis for revised TE.
Process Safety and Environmental Protection 9 3 ( 2 0 1 5 ) 75–88 87
RW
3
Ti
FiedgampcghrMsMatrm
abuTItt
Rc
epFtmt
4
Aaptt
Fig. 6 – R
.9. Ranking of MCs
able 13 presents the ranking of MCs based on their calculatedmportance levels (Eq. (14)).
In a SA, an input data parameter, such as a componentP is changed, and the resulting change in the TE probabil-ty is determined. This is repeated for a set of changes usingither different values for the same parameter or changingifferent parameters, e.g. changing different FPs. Usually for aiven sensitivity evaluation, only one parameter is changed at
time. This is called a one-at-a-time sensitivity study. Thisethod is employed here to validate the sensitivity of the
roposed model. RRW is employed to perform SA. The RRWan be calculated by re-quantifying the MCs probability of theiven event set to 0. It is expected that eliminating of MCs thatave the highest contribution to the occurrence of TE shouldesult in reducing the occurrence rate of TE more than otherCs. Therefore, ranking of RRW values is expected to be the
ame as the initial ranking result of MCs. As shown in Table 14,Cs31 has the highest contribution to the TE occurrence prob-
bility. Therefore, the RRW value of MCs31 is 0.0164 which ishe highest as expected. It shows the ranking result whichemains the same as the previous calculation. The proposed
odel satisfies the aforementioned expectations.The fourth column of Table 14 shows the value of the TE
fter eliminating of MCs. Fig. 4 includes 40 bars; the red colorar shows the TE value which is 0.3024. All new TE values (Col-mn 5th of Table 14) are presented in blue bars. As shown inable 14, MC number 26 is the most critical MC of the system.f the value of MC 26 is reduce to zero, it would be expectedhe new TE value reduce more than others. Fig. 5 can confirmhe expectation.
FVI and RRW are employed for ranking of MCs in this paper.esults of FVI measure and RRW are shown in the 3rd and 6tholumns of Table 14.
As mentioned, one of the advantages of SA is to identifyrror in the model. Result of Table 14 and Fig. 5 show that theroposed model can produce robust outcomes. As obvious inig. 6, controlling the first 11 most critical MCs would reducehe TE rate from 0.3 to 0.15. It means that control of the deter-
ined MCs will ensure considerable safety improvement inhe Deethanizer section of the Arya Sasol Plant.
. Conclusion and discussions
s one of the heavy industry discipline, petrochemical plantsre required to implement effective and consistent safetylanning against potential hazards (i.e. fire, explosion and
oxic gas releases) in order to ensure sustainable produc-ion. This research focused on developing a FFTA methodology
results.
demonstrating with Deethanizer failure within petrochemicalplant operational concept. Consequently, the research high-lights the following points:
• A fuzzy methodology for FT evaluation seems to be analternative solution to overcome the weak points of the con-ventional approach.
• By using linguistic variables, it is possible to handle theambiguities in the expression of the occurrence of a BE. Inaddition, the state of each BE can be described in a moreflexible form, by using the concept of fuzzy sets.
• Instead of using CFP, FP is used to characterize the failureoccurrence of the system events. It can efficiently expressthe vagueness of the nature of system phenomena andinsufficient information. Further, regardless of the complex-ity of the system, it is also possible to identify which BE caninfluence system FP the most.
• The importance measure can provide useful information forimproving the safety performance of a system. F-VI mea-sure index assists the analyst in identifying the critical MCsin the system for reducing occurrence likelihood of a TE.
As further research, application of FFTA methodology toother critical processes in petrochemical plant can be con-ducted. Moreover, multi attribute decision making (MADM)techniques can be adopted into the proposed methodology tomake cost benefit analysis for controlling the determined MCs.
Acknowledgement
The authors gratefully acknowledge to HSE manager (HosseinCharkhand) of Arya Sasol Petrochemical Company (A.S.P.C.)for technical information support in the demonstration of thecase study on Deethanizer failure in petrochemical plant oper-ations.
References
Aqlan, F., Ali, M.E., 2014. Integrating lean principles and fuzzybow-tie for risk assessment in chemical industry. J. LossPrevent. Process Ind. 29, 39–48.
Cai, K.Y., Wen, C., Zhang, M., 1991. Fuzzy states as a basis for atheory of fuzzy reliability in the possibility context. Fuzzy SetsSyst. 42, 17–32.
Celik, M., Lavasani, S.M., Wang, J., 2010. A risk-based modellingapproach to enhance shipping accident investigation. Saf. Sci.48, 18–27.
Huang, H.Z., Tong, X., Zuo, M., 2004. Posbist fault tree analysis ofcoherent systems. Reliab. Eng. Syst. Saf. 84, 141–148.
Korta, J.P., Lee, M.P., Eisenberg, N., Dewispelare, A., 1996. BranchTechnical Position on the Use of Expert Elicitation in theHigh-level Radioactive Waste Program, 1st ed. U.S. NuclearRegulatory Commission, NUREG-1563.
Lavasani, S.M., Wang, J., Yang, Z., Finlay, J., 2012. Application ofMADM in a fuzzy environment for selecting the best barrierfor offshore wells. Exp. Syst. Appl. 39, 2466–2478.
Rausand, M., Hoyland, A., 2004. System Reliability Theory:Models, Statistical Methods and Applications, 1st ed. JohnWiley and Son, Inc., New Jersey, USA.
Sawer, J.P., Rao, S.S., 1994. Fault tree analysis of fuzzy mechanicalsystems. Microelectron. Reliab. 34, 653–667.
Shu, M.H., Cheng, C.H., Chang, J.R., 2006. Using intuitionisticfuzzy fault-tree analysis on printed circuit board assembly.Microelectron. Reliab. 46, 2139–2146.
Spouge, J., 2000. A Guide to Quantitative Risk Assessment forOffshore Installation, 2nd ed. CMPT publications, Aberdeen,UK.
Sugeno, M., 1999. Fuzzy Modelling and Control, 1st ed. CRC Press,Florida, USA.
Tanaka, H., Fan, L., Lai, S., Toguchi, K., 1983. Fault tree analysis byfuzzy probability. IEEE Trans. Reliab. 32, 453–457.
Wang, D., Zhang, P., Chen, L., 2013. Fuzzy fault tree analysis forfire and explosion of crude oil tanks. J. Loss Prevent. ProcessInd. 26, 1390–1398.