Back to Publications
Process Engineering5 min read

Fault Tree (FTA) & Event Tree (ETA) Analysis in API Plants: Quantifying Reactor Runaway Risk

Kiran SeepanaSeptember 15, 20261251 Views
Go VIP (Watermark-Free)
Executive Summary & Scope

Fault Tree (FTA) & Event Tree (ETA) Analysis in API Plants: Quantifying Reactor Runaway Risk

Peer-Reviewed & PE Verified

ASME VIII • NFPA 68/69 • TEMA • ISO 9001 Alignment

This technical publication and associated design calculations have been reviewed for engineering consistency, unit integrity, and alignment with standard process design practices (Process Engineering).

# Fault Tree Analysis (FTA) & Event Tree Analysis (ETA) in API Manufacturing: Quantitative Reactor Runaway Probability & Minimal Cut Sets

When evaluating complex, multi-component safety systems in pharmaceutical API plants, process safety engineers require quantitative risk assessment (QRA) tools capable of identifying how multiple individual component failures combine to trigger a major industrial disaster.

  • Fault Tree Analysis (FTA) is a top-down, deductive analytical methodology that starts with an unwanted catastrophic outcome (Top Event) and works backward through Boolean logic gates to identify all combinations of basic equipment, instrument, and human failures that can cause it.
  • Event Tree Analysis (ETA) is a bottom-up, inductive analytical methodology that starts with a specific initiating event and maps forward through successive safeguard success or failure branches to quantify the frequencies of all potential end-state outcomes.

This technical guide details the Boolean algebra of FTA/ETA, the computation of Minimal Cut Sets (MCS), Common Cause Failure (CCF\text{CCF}) modeling, and a fully worked quantitative FTA/ETA case study for a batch reactor thermal runaway.


# 1. Fault Tree Analysis (FTA) Logic & Boolean Math

Fault trees model system vulnerabilities using standardized gate symbols:

          [ AND GATE ]                             [ OR GATE ]
     Output occurs ONLY if                   Output occurs if AT LEAST ONE
   ALL inputs occur together                      input occurs
           P = P1 x P2                        P = 1 - (1-P1)(1-P2)

# 1.1 Probability Math Formulas

# 1. AND Gate (Simultaneous Safeguard Failures)

For independent basic events E1,E2,…,EnE_1, E_2, \dots, E_n with failure probabilities P1,P2,…,PnP_1, P_2, \dots, P_n:

PAND=∏i=1nPi=P1×P2×⋯×PnP_{\text{AND}} = \prod_{i=1}^{n} P_i = P_1 \times P_2 \times \dots \times P_n

Key Insight: Connecting independent safeguards under an AND gate exponentially reduces top-event probability.

# 2. OR Gate (Single Point of Failure Path)

POR=1−∏i=1n(1−Pi)≈∑i=1nPi(for small Pi≪1)P_{\text{OR}} = 1 - \prod_{i=1}^{n} (1 - P_i) \approx \sum_{i=1}^{n} P_i \quad (\text{for small } P_i \ll 1)

Key Insight: Any basic event connected through an OR gate represents a single-point failure path that directly threatens plant safety.


# 2. Minimal Cut Sets (MCS): Finding Latent System Vulnerabilities

A Cut Set is any combination of basic events that, if they occur simultaneously, guarantees the occurrence of the Top Event.

A Minimal Cut Set (MCS) is a cut set from which no basic event can be removed without losing the guarantee of Top Event occurrence. It represents the smallest combination of component failures necessary and sufficient to cause a catastrophe.

Top Event Logic: T=MCS1∪MCS2∪⋯∪MCSk\text{Top Event Logic: } T = \text{MCS}_1 \cup \text{MCS}_2 \cup \dots \cup \text{MCS}_k
+----------------------------------------------------------------------------+
|                    MINIMAL CUT SET (MCS) ORDER RANKING                     |
|                                                                            |
| Order 1 Cut Set (Single Failure Path) ----> CATASTROPHIC VULNERABILITY     |
| Order 2 Cut Set (Dual Component Failures)-> Requires Active Redundancy     |
| Order 3 Cut Set (Triple Component Failures)-> High Safety Resilience       |
+----------------------------------------------------------------------------+

# 2.1 MOKOC Boolean Reduction Example

Using Boolean algebra rules (A⋅A=AA \cdot A = A, A+A⋅B=AA + A \cdot B = A):

Top Event T=E1⋅(E2+E1⋅E3)=E1⋅E2+E1⋅E1⋅E3=E1⋅E2+E1⋅E3\text{Top Event } T = E_1 \cdot (E_2 + E_1 \cdot E_3) = E_1 \cdot E_2 + E_1 \cdot E_1 \cdot E_3 = E_1 \cdot E_2 + E_1 \cdot E_3
Minimal Cut Sets: MCS1={E1,E2},MCS2={E1,E3}\text{Minimal Cut Sets: } \text{MCS}_1 = \{E_1, E_2\}, \quad \text{MCS}_2 = \{E_1, E_3\}

# 3. Common Cause Failure (CCF) & Beta-Factor Modeling

Independent redundancy can be compromised by shared environmental, maintenance, or calibration errors (Common Cause Failures).

The total failure rate λ\lambda of a redundant sensor pair is modeled using the β\beta-factor model:

λtotal=λs (Independent Failure)+λc (Common Cause Failure)\lambda_{\text{total}} = \lambda_s \text{ (Independent Failure)} + \lambda_c \text{ (Common Cause Failure)}
λc=β⋅λtotal\lambda_c = \beta \cdot \lambda_{\text{total}}
λs=(1−β)⋅λtotal\lambda_s = (1 - \beta) \cdot \lambda_{\text{total}}

Where β\beta typically ranges from 5%5\% to 10%10\% in industrial instrumentation.


# 4. Worked Case Study: Batch Reactor Runaway FTA & ETA

# 4.1 Top Event Definition

Catastrophic Batch Reactor Rupture & Flammable Vapor Cloud Explosion (PTEP_{\text{TE}} per year).

                             [ TOP EVENT: Reactor Rupture ]
                                           |
                                      ( AND GATE )
                    +----------------------+----------------------+
                    |                                             |
      [ Reactor Overpressure ]                        [ Relief System Failure ]
            ( OR GATE )                                      ( OR GATE )
      +---------+---------+                            +---------+---------+
      |                   |                            |                   |
 [ Runaway ]    [ Loss of Cooling ]               [ PSV Blocked ]   [ Undersized Relief ]

# 4.2 Basic Event Probabilities (PiP_i per year)

  • E1E_1: Jacket Cooling Water Pump Trip (P1=1.0×10−1 /yrP_1 = 1.0 \times 10^{-1} \text{ /yr})
  • E2E_2: Temperature Controller Sensor Drift (P2=5.0×10−2 /yrP_2 = 5.0 \times 10^{-2} \text{ /yr})
  • E3E_3: Operator Fails to Respond to High TT Alarm (P3=1.0×10−1P_3 = 1.0 \times 10^{-1})
  • E4E_4: High-High Temp SIS Interlock Valve Fails to Close (P4=1.0×10−2P_4 = 1.0 \times 10^{-2})
  • E5E_5: Pressure Safety Valve (PSV) Isolated or Blocked (P5=1.0×10−2P_5 = 1.0 \times 10^{-2})

# 4.3 Quantitative FTA Calculation

  1. Heat Accumulation Gate (G1G_1, OR Gate):
P(G1)≈P1+P2=1.0×10−1+5.0×10−2=0.15 /yearP(G_1) \approx P_1 + P_2 = 1.0 \times 10^{-1} + 5.0 \times 10^{-2} = 0.15 \text{ /year}
  1. Uncontrolled Runaway Gate (G2G_2, AND Gate):
P(G2)=P(G1)×P3×P4=0.15×0.1×0.01=1.5×10−4 /yearP(G_2) = P(G_1) \times P_3 \times P_4 = 0.15 \times 0.1 \times 0.01 = 1.5 \times 10^{-4} \text{ /year}
  1. Top Event Gate (GTOPG_{\text{TOP}}, AND Gate with PSV Failure):
PTE=P(G2)×P5=1.5×10−4×0.01=1.5×10−6 /yearP_{\text{TE}} = P(G_2) \times P_5 = 1.5 \times 10^{-4} \times 0.01 = 1.5 \times 10^{-6} \text{ /year}

# 5. Event Tree Analysis (ETA) Forward Branching Matrix

Starting from the Initiating Event (E1E_1: Cooling Water Pump Trip, fIE=0.1 /yearf_{\text{IE}} = 0.1 \text{ /year}):

+---------------------------------------------------------------------------------------------------------+
| Initiating Event | High T Alarm (P3) | SIS Interlock (P4) | Emergency Quench | End-State Outcome        | Frequency (f/yr)
+---------------------------------------------------------------------------------------------------------+
| CW Pump Trip     | Success (0.90) ---+--------------------+------------------+-> 1. Controlled Safe Stop| 9.0 x 10^-2
| (f = 0.1/yr)     |                   |                    |                  |                          |
|                  | Fail (0.10) ------+-> Success (0.99) -+------------------+-> 2. Automated SIS Trip | 9.9 x 10^-3
|                  |                   |                    |                  |                          |
|                  |                   | Fail (0.01) -------+-> Success (0.90)-> 3. Quenched Safe Hold | 9.0 x 10^-5
|                  |                   |                    |                  |                          |
|                  |                   |                    +-> Fail (0.10) ---> 4. CATASTROPHIC EXPLOSION| 1.0 x 10^-5
+---------------------------------------------------------------------------------------------------------+

# 5.1 Outcome Frequency Summary

  • Outcome 1 (Controlled Stop): 9.0×10−2 /yr9.0 \times 10^{-2} \text{ /yr} (90.0%90.0\%)
  • Outcome 2 (SIS Automated Trip): 9.9×10−3 /yr9.9 \times 10^{-3} \text{ /yr} (9.9%9.9\%)
  • Outcome 3 (Manual Quench): 9.0×10−5 /yr9.0 \times 10^{-5} \text{ /yr} (0.09%0.09\%)
  • Outcome 4 (Catastrophic Explosion): 1.0×10−5 /yr1.0 \times 10^{-5} \text{ /yr} (0.01%0.01\%)

# 6. Applicable Engineering Standards & Codes

  • IEC 61025: Fault Tree Analysis (FTA).
  • ISO 17776: Petroleum and Natural Gas Industries - Offshore Production Installations - Guidelines on Tools and Techniques for Hazard Identification and Risk Assessment.
  • NUREG-0492: Fault Tree Handbook (U.S. Nuclear Regulatory Commission).
Process Safety & Engineering Publication Disclaimer
Terms of Service

This technical article is published strictly for educational, design modeling, and process engineering reference. Because chemical synthesis, solvent flammability envelopes, reaction calorimetry kinetics, and toxic off-gassing vary by plant configuration, all scale-up decisions, emergency relief sizing, and safety interlocks must comply with local statutory regulations (e.g., OSHA PSM 1910.119, PESO, ATEX/IECEx, NFPA) and undergo formal Process Hazard Analysis (PHA/HAZOP) prior to physical execution.

PharmaChemEng Technical Knowledgebase • www.pharmachemeng.com
Process EngineeringPharmaScale-Up
Comments (0)

Discussion

Please Log In to participate in the technical discussion.

No comments posted yet. Be the first to share your input!