Week 14
Week 14
CH4051_Process Safety
Risk Analysis and Assessment
Prepared by:
For:
Ho Chi Minh City University of Technology
Learning Objectives
1. Probability theory.
2
Introduction
Hazards
Identification
Risk
Assessment
3
Probability Theory
4
Probability Theory
5
Probability Theory
POISSON DISTRIBUTION Probability R(t) component will not fail during (0,t)
This is called the reliability:
t t t t
6
Probability Theory
Example: A device is found to fail once every 2 years. What is the failure
rate, the failure probability and the reliability at the end of 1 year, and the
MTBF?
Answer:
8
Interaction
Look up the failure rates for these three components from Table 12-1.
Then compute the reliability and failure probability for each
component for 1 year time period.
(faults/yr) P=1-R
Control valve 0.60 0.55 0.45
Controller 0.29 0.75 0.25
DP cell 1.41 0.24 0.76
R = 0.10 = e −
= − ln ( 0.10 ) = 2.30 failures/year
11
(Un)revealed Failures
COMPONENT STATUS
Inspection
Unnoticed
Repaired failure Repaired
UNREVEALED FAILURE u 0 i 1
u 1
Unavailability ( r ) U= i Half of inspection interval
i 2
12
(Un)revealed Failures
13
(Un)revealed Failures
1 1
U = i = ( 2.30 / yr )( 0.083 yr ) = 0.0955
2 2
A = 1 − U = 1 − 0.0955 = 0.905
14
Probability Concept Application
Probability Theory
- Failures in Process Industries
Failure of a single
▪ Failure means any component or
release event caused single action
Single
by failure of Component
Failure
equipment/com Data for
ponents, or The Initiating failure rates
operational failure. event (mostly) are compiled
can be by industry
▪ Assessing risk is an attributed to
integral part of hazard the failures of Failures
management. component in resulting from
process several failures
▪ In general, risk industries: Multiple
Component and/or actions
represents the Failure
likelihood of the level Failure rates
of harm. determined
using FTA
15
Failures in Process
Probability Industries
Theory
Initiating Event (IE) Outcome
Event
(IE)
Failure of single
component
hazardous events
Catastrophic
16
TABLE 12-1
Failure Rate Data for Various Selected Process Components
Instrument Faults/Year
Controller 0.29
Control valve 0.60
Flow measurement (fluids) 1.14
Flow measurement (solids) 3.75
Flow switch 1.12
Gas–liquid chromatograph 30.6
Hand valve 0.13
Indicator lamp 0.044 Basic Fact:
Level measurement (liquids) 1.70
Level measurement (solids) 6.86 The more complex
Oxygen analyzer 5.65
pH meter 5.88 the device the higher
Pressure measurement 1.41
Pressure relief valve 0.022 the failure rate!
Pressure switch 0.14
Solenoid valve 0.42
Stepper motor 0.044
Strip chart recorder 0.22
Thermocouple temperature measurement 0.52
Thermometer temperature measurement 0.027
Valve positioner 0.44
17
Types of Failure Rate Data
Probability and Sources
Theory
18
Inspection
Probability&Theory
Testing
A test can be performed to estimate its failure rate. Ten identical components
are each tested until they either fail or reach 1000 hours, at which time the test
is terminated for that component. (The level of statistical confidence is not
considered in this example.) The results are as follows:
21
Example of Failure
Probability Rate Data
Theory
22
Other Sources – Failure
Probability Rate Data
Theory
Analysis Techniques
1. Frequency modelling techniques Used to estimate
2. Common-cause failure analysis frequencies or probabilities
3. Human reliability analysis from basic data. Typically
4. External events analysis used when detailed
Used historical data is not
available.
i. EVENT TREES
ii. FAULT TREES
24
Event Tree Analysis
Probability Theory
❑ An Event Tree:
• Begins with an Initiating Event
• Evaluates the impact of the Safety Function on the accident process
• Allows calculation of Failure
Probability for various scenarios
25
Generic Example
Probability Theory
26
Computational
Probability sequence
Theory
Example: Brake Fails – Barrier described negatively
Brake fails = 0.005 + 0.015 + 0.045 + 0.135 + 0.02 + 0.06 + 0.18 = 0.46
27
Computational
Probability sequence
Theory
Example: Brake Fails – Barrier described negatively
Failures and
Successes of
Initiating Event Various Defined
Various Intervening
(Cause) Final Outcomes
Safety
- these have an - These will have
Systems/Actions
associated associated
- These have an
frequency frequencies
average Probability
on Demand
29
Event Tree Analysis
Example
ID B (High Temp Alarm Alerts Operator)
0.01 failures/demand 1.We’ll call the initiating event A
and also note the occurrence per
year.
Success of
Safety 2.Draw a line from the initiating
Function B event to the first safety function (ID
A B) – a straight line up indicates the
Loss of coolant results for a success in the safety
(initiating event) Failure of function and a failure is represented
Safety by a line drawn down.
1 occurrence/year
Function B 3.We can assume the high temp
alarm will fail to alert the operator
1% of the time when in demand OR
0.01 failure/demand.(This is a
probability of failure on demand)
40
30
Event Tree Analysis
Example
Safety Function
ID B (High Temp Alarm Alerts Operator)
0.01 failures/demand
31
Event Tree Analysis
Example
• Initiating event: Loss of
Coolant
33
Evaluate Frequency of Loss Event
• Event tree analysis explores outcome • From the event tree the net failure
severities, probability, and risk frequency is the sum of the unsafe
• Is the risk tolerable? Can it be state frequencies
accepted? • Frequency = ADE +ABDE + ABCDE
• Consider additional safety functions = 0.025 failure/yr
to reduce risk (probability, severity) to = 1 failure every 40 years
tolerable levels.
• The corresponding risk is
• Select safety functions on the basis of considered too high, so the
✓ Effectiveness, operability frequency must be reduced.
✓ Cost, time • Add a high-temp reactor shutdown
✓ Feasibility (including schedule) system. Set the shutdown temp
above the alarm value to allow
operator to adjust coolant flow.
34
Potential for Reducing Frequency
▪ The event tree analysis shows that a dangerous runaway reaction will
occur on average 0.025 time per year, or once every 40 years. This is
considered too high for this installation.
▪ A possible solution is the inclusion of a high-temperature reactor
shutdown system. This control system would automatically shut down
the reactor in the event that the reactor temperature exceeds a fixed
value.
▪ The emergency shutdown temperature would be higher than the alarm
value to provide an opportunity for the operator to restore the coolant
flow.
35
Potential for Reducing Frequency
37
Fault Tree Analysis
DEDUCTIVE METHOD:
• Well-defined accident, top event of interest
• Possible causative hazards
• Quantitative: add failure data to fault tree
Minimal cut set – all possible ways to get to top event Logical functions
AND, OR, INHIBIT
Drawbacks • large trees computers & software
• never sure all failure modes are covered
• no partial failures
• accurate failure data required reference libraries
Additional Comment: The top event of a fault tree becomes the initiating
event of the event tree. This results in a BOW TIE diagram. 19
Fault Tree Symbols
20
Logical Connection Between Events
AND
C
The resulting output event requires the simultaneous
occurrence of all input event.
e.g., Event C will occur only if both events A and B occur
simultaneously, which is represented (for independent
A B
events, A, B) by
A•B=C AND-gate
OR
C
The resulting output event requires the occurrence of any
individual input event,
e.g., C will occur if either A or B occurs, which is
represented (when A, B each has low probability of A B
occurring) by
A+B=C OR-gate
21
Preliminary Steps - Fault Tree Construction
❑ Put the hazardous event (TOP EVENT) on the left-hand side (or top) of the
page.
❑ All immediate possible causes of the top event are identified and
placed next to it on the tree.
❑ If any of these events can cause the top event, they are joined to it by an
OR gate .
❑ If all are required before the top event occurs, they are
joined to it by an AND gate.
❑ If a combination of gates appears to be necessary at any point, then progress
is too rapid and a suitable intermediate stage should be sought.
❑ Describe these intermediate events in terms of earlier
events, using either AND or OR logic.
❑ Repeat this for each branch of the tree until no further detail is
necessary or possible.
21
Fault Tree Construction
21
Example: Fault tree model for hydride dust
explosion
21
Unit on Fault Tree Analysis and Rules
OR GATE rules :
can add P and P = unit of P
can add F and F = unit F
cannot add F and P =different unit
21
Example Fault Trees Analysis
22
Relationship between Fault Tree Analysis
and Event Tree Analysis
▪ Event trees
begin with an initiating event and work toward the top event
(induction).
▪ Fault trees
begin with a top event and work backward toward the
initiating events (deduction).
▪ Top events for Fault Trees are initiating events for Event Trees.
22
Advantages Of Fault Trees
ADVANTAGES:
23
Disadvantages Of Fault Trees
DISADVANTAGE:
• Fault tree only assumes ‘ hard ’ failures. Partial failures are not
considered.
24
Bow-tie Diagram
The top event from a fault tree becomes the initiating event for an
event tree.
25
Bow-tie Diagram
Incident
(Loss Event)
Initiating Events
(Causes)
Outcomes
Preventive Mitigative
Safeguards Safeguards
26
Bow-tie Diagram
Mitigative
Initiating Events Safeguards
(Causes) Probability Outcomes
of Ignition
Flash Fire
Proactive Probability of
Explosion
Safeguards
Vapor Cloud Explosion
Probability
Control Failure of Ignition
Preventive Building Explosion
Time at Safeguards Probability
Risk of Ignition
Human Error Fireball
Physical Explosion
Incident
Mechanical Failure (Loss Event) Chemical Exposure
Onsite Toxic
Each feasible path between an
Toxic Infiltration
initiating event and an outcome
represents a scenario with Offsite Toxic
applicable protective layers. 27
Bow-tie Diagram
1) Hazards
2) Initiating Events
3) Enabling Conditions
7) Mitigating Safeguards
8) Outcomes
28
9) Impacts
QRA: Quantitative Risk Assessment
Hazards
Identification Selection of
Release Incident
Selection of
Dispersion Model
Flammable Toxic
Flammable
Risk and/or Toxic?
Analysis
Selection of Selection of
Fire and Effect Model
Explosion Model
Mitigation
Factors
Risk
Assessment Consequence
Model
29
QRA: Quantitative Risk Assessment
• Define the initiating events and the incident sequence. For example, a
cooling water failure causes a runaway reaction that overpressures
the reactor vessel, causing the relief to open, discharging the reactor
contents.
• Use source models to estimate the discharge rate. For the reactor
example, this would require a source model to estimate the discharge
rate through the relief. (See Chapter 4.)
• Use a dispersion model to estimate the chemical concentrations
downwind of the release. (See Chapter 5.)
• Estimate the incident consequences for people, environment, and
property using effect models. (See Chapter 3.)
30
QRA: Quantitative Risk Assessment
• Estimate the potential incident frequencies using event trees and fault
trees.
• Estimate the risk by combining the consequences and frequencies.
• Combine the risk estimates for all the scenarios to estimate the
overall risk.
• Decide if the risk is tolerable. (See Sections 1-9 and 12-7 in 4th edition)
31
Type of Risk: Definitions
Many Hazards
Individual
Single Hazard
People
32
Voluntary and Involuntary Risk
Voluntary Risk – Risk that is consciously tolerated by someone seeking to
obtain the benefits of the activity that poses the risk.
Examples: Riding a car,
Riding a motorcycle,
Mountain climbing,
Skiing.
Involuntary Risk – Risk that is imposed on someone who does not directly
benefit from the activity that poses the risk.
Examples: Living in the vicinity of a chemical plant,
Riding a train,
Riding an airplane,
Visiting a mall. 33
QRA: Individual Risk – Risk Contours
10-6
10-5
Plant
10-4
Community
34
QRA: Individual Risk – Risk Contours
10-2
Broadly Unacceptable
10 -3
B
This is a plot of the
C
10-4 various societal risk
Fatalities/Year
D
10-5
E curves used throughout
10-6
10-9
1 10 100 1000
N, Fatalities
• A: United Kingdom’s Health Safety Executive (HSE) Criteria—Maximum tolerable societal risk
• B: Dutch—Maximum tolerable societal risk
• C: U.K. HSE—Negligible risk to workers and public
• D: New South Wales—Negligible societal risk
• E: Hong Kong—Acceptable societal risk 36
Example: F-N Curve
Use the data provided in Table 12-6 to draw an F-N curve.
1. The incident outcome case with the smallest number of fatalities is selected
first. This is case 5, which has one fatality. Case 2 is not selected because it has
zero fatalities.
2. All incident outcome cases with one or more fatalities are selected. These are
cases 1, 3, 4, and 5. The frequencies for these cases are added together to
create the plot shown below: 37
Example: F-N Curve
3. The case with the next highest number of fatalities is selected. This is
case 4, with three fatalities. The frequencies for all cases with three or
more fatalities are added together, as shown in Table 12-7.
4. This procedure is repeated for 6+, 13+, and greater than 13 fatalities.
The results are shown in the table below.
5. The data are plotted in figure below. The results are extrapolated to N =
1. The vertical lines are drawn at the actual number of fatalities. Note
that the results exceed many of the societal risk criteria in Figure 12-20.
38
Example: F-N Curve
39
Example: F-N Curve
-3
10
Frequency of N or More Fatalities, per year
-5
10
-6
10
1 10 100
Number of Fatalities, N
40
Definition of Risk Tolerance
• Risk Tolerance or acceptance is defined as “the maximum level of risk of a
particular technical process or activity that an individual or organization
accepts to acquire the benefits of the process or activity”
• Risk tolerance may also change with time as society, regulatory agencies
and individuals come to expect more from the chemical industry.
• As a consequence, a risk that was considered tolerable years ago may now
be deemed unacceptable. 41
QRA: Risk Tolerance Criteria
The recommended tolerability criteria based on the DOE risk criteria are as
follow (Environment Impact Assessment (EIA) Guidelines for Risk
Assessment, 2004):
The 1 x 10-6 fatalities / person per year individual risk contour should not
encompass involuntary recipients of industrial risks such as residential
area, schools, hospitals and places of continuous occupancy, etc.
The 1 x 10-5 fatalities / person per year individual risk contour should not
extend beyond industrial developments
The tolerability criteria based on the PETRONAS Risk Criteria as follow (PTS
60.2210 Quantitative Risk Assessment, 2006):
The annual risk greater than 1 x 10-3 fatalities / person per year represents
intolerable risk.
42
Risk Tolerance – Risk Matrix
Risk Matrix Likelihood
1. Select the severity from the highest box in either of columns 1, 2 or 3. Read the 4 5 6 7
Category and Safety Severity Level from the same row. LIKELY UNLIKELY IMPROBABLE IMPROBABLE.
BUT NOT
2. Select the likelihood from columns 4 thru 7. IMPOSSIBLE
3. Read the Risk Level from the intersection of the severity row and the likelihood Expected to
column. Expected to happen Expected to Not expected to
happen possibly happen possibly happen anywhere
TMEF: Target mitigated event frequency several times once over once in the in the division
TQ: Threshold Quantity over the life of the life of the division over the over the life of the
the plant. plant. life of the plant. plant
1 2 3 Safety 0 to 9 10 to 99
Human Health Fire, Explosion Chemical Severity Severity ≥ 100 years > 1000 years
years years
Impact Direct Cost in $ Impact Category Level
Public fatality 4
possible, Greater than Risk Level Risk Level Risk Level Risk Level
$10 MM ≥ 20x TQ CATASTROPHIC TMEF =
employee A A B C
1×10-6
fatalities likely
Severity
aLosttime injury (LTI): The injured worker is unable to perform regular job duties, takes time off for recovery, or is assigned modified work duties while
recovering.
bRecordable injury: Death, days away from work (DAW), restricted work or transfer to another job, medical treatment beyond first aid, or loss of
consciousness.
43
Table 1-15: Risk matrix for semi-quantitative classification of incidents.
Table 1-16 Threshold quantities (TQ) for a variety of chemicals. Source: AICHE/CCPS
2,000 kg = 4,400 lbm Ethyl acetate 200 kg = 440 lbm
Acrylamide Ethyl benzene Ammonia, anhydrous
Ammonium nitrate fertilizer Ethylenediamine Carbon monoxide
Amyl acetate Formic acid
Amyl nitrate Heptane 100 kg = 220 lbm
Bromobenzene Hexane Hydrogen bromide, anhydrous
Calcium oxide Methacrylic acid Hydrogen chloride, anhydrous
provided in Reference
1,000 kg = 2,200 lbm Carbon disulfide Nitric oxide, compressed
Acetic anhydride Cyclobutane Nitrogen trioxide
Acetone Diethyl ether or Ethyl ether Phosgene
Acetonitrile Ethane Phosphine
Aldol
Ammonium perchlorate
Ethylamine
Ethylene
Stibine
materials on course web page.
Aniline Furan
Arsenic Hydrazine, anhydrous
Barium Hydrogen, compressed
Benzene Lithium
Benzidine Methylamine, anhydrous
Butyraldehyde
Carbon tetrachloride
Potassium
Potassium cyanide
Each company customizes the
Coper chlorate Propylene oxide
Copper cyanide
Cycloheptane
Silane
Sodium
risk matrix for their operation.
Cycloheptene Sodium cyanide
Cyclohexene Sodium peroxide
Dioxane Trichlorosilane 44
Epichlorohydrin
Example – Risk Matrix
Solution: The Threshold Quantity (TQ) for acetone from the table is
1,000 kg. The release of 1,500 kg is 1.5 times the TQ. From Column 3
of the Risk Matrix – Chemical Impact - this is a MINOR severity
category. From the financial loss of $1,500,000, under column 2 of the
Risk Matrix – Fire, Explosion Direct Cost in $ - this is VERY SERIOUS.
45
Risk Tolerance – Risk Matrix
Risk Matrix Likelihood
1. Select the severity from the highest box in either of columns 1, 2 or 3. Read the 4 5 6 7
Category and Safety Severity Level from the same row. LIKELY UNLIKELY IMPROBABLE IMPROBABLE.
BUT NOT
2. Select the likelihood from columns 4 thru 7. IMPOSSIBLE
3. Read the Risk Level from the intersection of the severity row and the likelihood Expected to
column. Expected to happen Expected to Not expected to
happen possibly happen possibly happen anywhere
TMEF: Target mitigated event frequency several times once over once in the in the division
TQ: Threshold Quantity over the life of the life of the division over the over the life of the
the plant. plant. life of the plant. plant
1 2 3 Safety 0 to 9 10 to 99
Human Health Fire, Explosion Chemical Severity Severity ≥ 100 years > 1000 years
years years
Impact Direct Cost in $ Impact Category Level
Public fatality 4
possible, Greater than Risk Level Risk Level Risk Level Risk Level
$10 MM ≥ 20x TQ CATASTROPHIC TMEF =
employee A A B C
1×10-6
fatalities likely
Severity
aLosttime injury (LTI): The injured worker is unable to perform regular job duties, takes time off for recovery, or is assigned modified work duties while
recovering.
bRecordable injury: Death, days away from work (DAW), restricted work or transfer to another job, medical treatment beyond first aid, or loss of
consciousness.
46
Table 1-15: Risk matrix for semi-quantitative classification of incidents.
72