Estimation of An Event Occurrence for LOPA Studies. Randy Freeman S&PP Consulting Houston, TX

Similar documents
Statistical Intervals (One sample) (Chs )

We're in interested in Pr{three sixes when throwing a single dice 8 times}. => Y has a binomial distribution, or in official notation, Y ~ BIN(n,p).

We're in interested in Pr{three sixes when throwing a single dice 8 times}. => Y has a binomial distribution, or in official notation, Y ~ BIN(n,p).

Power of a hypothesis test

GEOMETRIC -discrete A discrete random variable R counts number of times needed before an event occurs

Hypothesis testing I. - In particular, we are talking about statistical hypotheses. [get everyone s finger length!] n =

MTH4451Test#2-Solutions Spring 2009

a. See the textbook for examples of proving logical equivalence using truth tables. b. There is a real number x for which f (x) < 0. (x 1) 2 > 0.

Applied Bayesian Statistics STAT 388/488

Preface Introduction to Statistics and Data Analysis Overview: Statistical Inference, Samples, Populations, and Experimental Design The Role of

Why Try Bayesian Methods? (Lecture 5)

EQ: How do I convert between standard form and scientific notation?

Name Solutions Linear Algebra; Test 3. Throughout the test simplify all answers except where stated otherwise.

EXAMPLE: INFERENCE FOR POISSON SAMPLING

CS5314 Randomized Algorithms. Lecture 15: Balls, Bins, Random Graphs (Hashing)

Key Words: Lifetime Data Analysis (LDA), Probability Density Function (PDF), Goodness of fit methods, Chi-square method.

P (E) = P (A 1 )P (A 2 )... P (A n ).

37.3. The Poisson Distribution. Introduction. Prerequisites. Learning Outcomes

Statistics. Lent Term 2015 Prof. Mark Thomson. 2: The Gaussian Limit

Saf-T-Graf Graphite Rupture Disks

Second-Order Homogeneous Linear Equations with Constant Coefficients

Probability Methods in Civil Engineering Prof. Rajib Maity Department of Civil Engineering Indian Institute of Technology, Kharagpur

Probability Methods in Civil Engineering Prof. Dr. Rajib Maity Department of Civil Engineering Indian Institution of Technology, Kharagpur

Week 1 Quantitative Analysis of Financial Markets Distributions A

CS 361: Probability & Statistics

(1) Introduction to Bayesian statistics

Why? 2.2. What Do You Already Know? 2.2. Goals 2.2. Building Mathematical Language 2.2. Key Concepts 2.2

Introduction to Probability, Fall 2013

Probability theory and inference statistics! Dr. Paola Grosso! SNE research group!! (preferred!)!!

Hypothesis Testing. ECE 3530 Spring Antonio Paiva

Percentage errors can ruin your day. Roland Martin, Stephan Kolassa June 27, 2012

Subject CS1 Actuarial Statistics 1 Core Principles

Learning Theory. Piyush Rai. CS5350/6350: Machine Learning. September 27, (CS5350/6350) Learning Theory September 27, / 14

Introduction to Survey Analysis!

Confidence intervals

Estimation of reliability parameters from Experimental data (Parte 2) Prof. Enrico Zio

Contingency Tables. Safety equipment in use Fatal Non-fatal Total. None 1, , ,128 Seat belt , ,878

Appendix A. Review of Basic Mathematical Operations. 22Introduction

Infinite Series. Copyright Cengage Learning. All rights reserved.

Learning Packet. Lesson 5b Solving Quadratic Equations THIS BOX FOR INSTRUCTOR GRADING USE ONLY

Regression, part II. I. What does it all mean? A) Notice that so far all we ve done is math.

Probability and (Bayesian) Data Analysis

Simulation. Where real stuff starts

24 = 5x. 4.8 = x. December 13, 2017

Extending the Number System

The t-distribution. Patrick Breheny. October 13. z tests The χ 2 -distribution The t-distribution Summary

BINOMIAL DISTRIBUTION

PVP STATISTICAL APPROACH OF STRESS CONCENTRATION FACTOR (SCF) WITHIN PIPELINE GIRTH WELDS ECA

A Short Introduction to the Theory of Statistical Learning. Kee Siong Ng, PhD

Advanced Herd Management Probabilities and distributions

Independence 1 2 P(H) = 1 4. On the other hand = P(F ) =

Econ 325: Introduction to Empirical Economics

Physics 403. Segev BenZvi. Credible Intervals, Confidence Intervals, and Limits. Department of Physics and Astronomy University of Rochester

CHAPTER 10 Comparing Two Populations or Groups

CHAPTER 10 Comparing Two Populations or Groups

Probability Method in Civil Engineering Prof. Dr. Rajib Maity Department of Civil Engineering Indian Institute of Technology, Kharagpur

5-7 The Pythagorean Theorem

Chapter 24. Comparing Means. Copyright 2010 Pearson Education, Inc.

Probability calculus and statistics

A Discussion of the Bayesian Approach

ANALYSIS OF INDEPENDENT PROTECTION LAYERS AND SAFETY INSTRUMENTED SYSTEM FOR OIL GAS SEPARATOR USING BAYESIAN METHODS

Bayesian Inference. STA 121: Regression Analysis Artin Armagan

Protean Instrument Dutchtown Road, Knoxville, TN TEL/FAX:

Algebra 1B notes and problems March 12, 2009 Factoring page 1

Math/Stat 352 Lecture 8

MATH 10 INTRODUCTORY STATISTICS

MAE 493G, CpE 493M, Mobile Robotics. 6. Basic Probability

I - Probability. What is Probability? the chance of an event occuring. 1classical probability. 2empirical probability. 3subjective probability

Conditional distributions

Contingency Tables. Contingency tables are used when we want to looking at two (or more) factors. Each factor might have two more or levels.

Lehmer Random Number Generators: Introduction

Simulation. Where real stuff starts

appstats27.notebook April 06, 2017

Math 562 Homework 1 August 29, 2006 Dr. Ron Sahoo

The Inductive Proof Template

Solutions to In-Class Problems Week 14, Mon.

Chapter 3. Estimation of p. 3.1 Point and Interval Estimates of p

2. I will provide two examples to contrast the two schools. Hopefully, in the end you will fall in love with the Bayesian approach.

3.4. The Binomial Probability Distribution

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 3

Math101, Sections 2 and 3, Spring 2008 Review Sheet for Exam #2:

BERTINORO 2 (JVW) Yet more probability Bayes' Theorem* Monte Carlo! *The Reverend Thomas Bayes

( )( b + c) = ab + ac, but it can also be ( )( a) = ba + ca. Let s use the distributive property on a couple of

Solving Quadratic & Higher Degree Equations

Principles of AI Planning

Statistical techniques for data analysis in Cosmology

THE POSTDOC VARIANT OF THE SECRETARY PROBLEM 1. INTRODUCTION.

(c) P(BC c ) = One point was lost for multiplying. If, however, you made this same mistake in (b) and (c) you lost the point only once.

Chapter 27 Summary Inferences for Regression

Math/Stats 425, Sec. 1, Fall 04: Introduction to Probability. Final Exam: Solutions

Slides for Data Mining by I. H. Witten and E. Frank

The Components of a Statistical Hypothesis Testing Problem

Homework 4 Solutions

Probability Theory and Statistics (EE/TE 3341) Homework 3 Solutions

Please bring the task to your first physics lesson and hand it to the teacher.

Lesson #33 Solving Incomplete Quadratics

Quantitative Understanding in Biology 1.7 Bayesian Methods

Confidence Intervals. Confidence interval for sample mean. Confidence interval for sample mean. Confidence interval for sample mean

Practice Exam 1. (A) (B) (C) (D) (E) You are given the following data on loss sizes:

p = q ˆ = 1 -ˆp = sample proportion of failures in a sample size of n x n Chapter 7 Estimates and Sample Sizes

Transcription:

Estimation of An Event Occurrence for LOPA Studies Randy Freeman S&PP Consulting Houston, TX 77041 713 408 0357 rafree@yahoo.com 1

Problem Your LOPA team members tell you that the initiating event of concern has never happened in the history of the process unit. No bad events in 20 years. What do you do? Ignore the event and move on It can never happen What else? 2

Available Methods to Handle Problem Assume event has happened once, M1 = 1/N Assume event has almost happened M2 = 1/2N Bayesian Results M3 = 1/3N M4 = 1/4N 3

Methods (cont) Poisson Arrivals of Failures At 90% limit M5 = 0.105/N = 1/10N At 10% limit M6 = 2.303/N 4

Methods (cont.) Chi Square Confidence Limit (95%) M7 = λ = χ α;2 / (2 n) = 5.991/(2 n) Chi Square Confidence Limit (50%) M8 = λ = χ α;2 / (2 n) = 1.3863/(2 n) 5

Methods (cont.) Binomial Distribution M9 = 1 [0.9] 1/n 6

Methods (cont.) Uniform Distribution Bayes Estimator M10 = 1/(N+2) 7

Methods (cont.) Hypothesis Test Normal distribution with 95% confidence M11 = ( n z + 2 α Zα = 1.6445 z 2 α ) 8

Methods (cont.) Explosive Initiation Test M12 = 1 [0.5] 1/N 9

1.00E+00 Best Guess Methods Legend Failure Rate, fail/yr. 1.00E-01 1.00E-02 M3 M1 M2 M3 M4 M8 M9 M10 1.00E-03 0 10 20 30 40 50 60 70 80 90 100 Nunber of Trials, N 10

1.00E+01 Confidence Limit Methods Failure Rate, fail/yr. 1.00E+00 1.00E-01 1.00E-02 M11 M12 M7 M6 Legend M5 M6 M7 M11 M12 M5 1.00E-03 0 10 20 30 40 50 60 70 80 90 100 Number of Trials, N 11

12

Recommendations for Basic LOPA Determine if the event of concern is physically possible. If not physically possible, delete from LOPA analysis Make sure you have a minimum of 10 years of data. If not use values from LOPA book Use M1 = 1/N for events found to be physically possible 13

Recommendations for Beyond Basic LOPA Reviews Use estimators M3, M8 or M12. M3 (Rule of 1/3N) is easier to remember and calculate. 14

Example Case 1 Data No occurrences in three years What frequency should be assigned? Answer Since the history is less than 10 years, use the frequency presented in LOPA book. 15

Example Case 2 Data One occurrence in ten years What frequency should be assigned? Answer Methods presented in this paper do not apply when an event has occurred. Start with a frequency of 1/10 and review value in LOPA book. 16

Example Case 3 Data No occurrences in ten years What frequency should be assigned? Answer Determine if history is valid Start with a frequency of 1/10 If LOPA at greater than order of magnitude significance, use 1/30 17

Example Case 4 Data No occurrences in 30 years What frequency should be assigned? Answer Determine if history is valid Start with a frequency of 1/10 If LOPA at greater than order of magnitude significance, use 1/90 18

Example Case 5 Data No occurrences in 20 years What frequency should be assigned? Answer Determine if history is valid No history of pump overpressure of PRV PRV set P = 150 psig Pump Deadhead P = 100 psig Not credible scenario Delete from LOPA Study 19

Can it really happen? We pressurize a vessel to 130 psig. The MAWP is 100 psig. What happens????? Does the vessel blow up???? 20

21

More Information Vessel in good condition per API 579 Fitness for Service Test Pressure = 1.5 * MAWP Test Pressure = 150 psig Piping is 2 inch schedule 40 with a catastrophic burst pressure of 7000 psig ANSI Class 150 Flanges with a MAOP of 285 psig at 100 F 22

ANSWER NOTHING HAPPENS! No rupture No weld tear No flange failure No release Nothing happens! Maybe an ASME code violation requiring an inspection for fitness for service 23

Lessons Learned Before starting detailed LOPA analysis verify that the scenario of interest can actually happen. Physics and Chemistry are powerful tools to sort out imaginary from real safety issues. 24

Conclusions There is no single correct answer Use an estimator that is consistent with the other assumptions in your LOPA study 25

Does the Audience Know of Another Estimator?? 26

More Details What to Do When Nothing Has Happened?, pp 204-211,Process Safety Progress, Vol. 30, No. 2, September 2011 Published online in Wiley Online Library (wileyonlinelibrary.com). DOI 10.1002/prs.10463, 2011 27

QUESTIONS???? 28

Have a Merry Christmas! 29