Back to Publications
Process Safety5 min read

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

Kiran SeepanaSeptember 6, 202612 Views
Executive Summary & Scope

A comprehensive guide on Fault Tree Analysis (FTA) and Event Tree Analysis (ETA) for pharmaceutical API reactors. Features Boolean algebra cut sets, quantitative runaway probability calculation, and SIL safety instrumented systems.

# 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=1i=1n(1Pi)i=1nPi(for small Pi1)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=MCS1MCS2MCSk\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 (AA=AA \cdot A = A, A+AB=AA + A \cdot B = A):

Top Event T=E1(E2+E1E3)=E1E2+E1E1E3=E1E2+E1E3\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×101 /yrP_1 = 1.0 \times 10^{-1} \text{ /yr})
  • E2E_2: Temperature Controller Sensor Drift (P2=5.0×102 /yrP_2 = 5.0 \times 10^{-2} \text{ /yr})
  • E3E_3: Operator Fails to Respond to High TT Alarm (P3=1.0×101P_3 = 1.0 \times 10^{-1})
  • E4E_4: High-High Temp SIS Interlock Valve Fails to Close (P4=1.0×102P_4 = 1.0 \times 10^{-2})
  • E5E_5: Pressure Safety Valve (PSV) Isolated or Blocked (P5=1.0×102P_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×101+5.0×102=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×104 /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×104×0.01=1.5×106 /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×102 /yr9.0 \times 10^{-2} \text{ /yr} (90.0%90.0\%)
  • Outcome 2 (SIS Automated Trip): 9.9×103 /yr9.9 \times 10^{-3} \text{ /yr} (9.9%9.9\%)
  • Outcome 3 (Manual Quench): 9.0×105 /yr9.0 \times 10^{-5} \text{ /yr} (0.09%0.09\%)
  • Outcome 4 (Catastrophic Explosion): 1.0×105 /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).
Fault Tree AnalysisFTAEvent Tree AnalysisETAMinimal Cut SetsProcess SafetyQuantitative Risk
Comments (0)

Discussion

Please Log In to participate in the technical discussion.

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