Qra Version F
Qra Version F
) ( 0058 . 0 17 . 1 (2)
where r = puncture resistance of pipe;
R = pipeline radius;
w = pipeline wall thickness;
L
t
= length of excavator tooth;
w
t
= width of excavator tooth;
u
= ultimate tensile strength of pipe steel.
Equation (2) has been calibrated against tests in a recent joint industry project [Ref. 29]. It
was shown that modelling uncertainty can be taken into account by modifying equation (2) as
follows:
C w w L
w
R
r
u t t
+ +
]
]
]
) ( 0058 . 0 17 . 1 (3)
where C is an additive error term to be described by an appropriate probability distribution.
The limit state function for puncture is then
F-r = 0 (4)
To apply the puncture limit state equation, the following variables require to be described by
probability distributions:
Excavator mass
Modelling uncertainty C
Ultimate tensile strength
Bucket tooth length and width
Pipe radius and wall thickness
CORRIB FIELD DEVELOPMENT PROJECT
Onshore Pipeline Quantified Risk Assessment
g:\052377 - enterprisecorrib additional work\ctr 140 - qra update\02-f-3-835- rev 05 - qra.doc Page 51 of 58
A1.3 Gouge / Dent Model
Reference 30 indicates that 18% of damage events in Transco transmission pipelines
involved dents containing gouges while 82% were plain gouges.
Gouges situated in dents, or plain gouges, are assessed using a fracture mechanics
approach assuming that the gouge behaves as a crack.
Using the level 2 assessment methodology of PD6493
[Ref. 31], a failure assessment
diagram (FAD) can be developed in terms of the fracture parameters K
r
and S
r
. K
r
is a
measure of how close the pipe is to failure by brittle fracture and S
r
is a measure of how close
the pipe is to failure by plastic collapse. These parameters are defined as follows:
( ) ( ) {
c
b b m m
r
K
a a Y a Y
K
1
+
(5)
,
`
.
|
,
`
.
|
w
a
w M
a
S
y
a
m
f
n
r
1 15 . 1
1
(6)
where a = gouge depth;
m
= membrane stress
=
,
`
.
|
R
D
h
2
8 . 1 1 ; (7)
b
= bending stress
=
R
D
w
R
h
2
2 . 10 ; (8)
R = pipe radius;
D = dent depth;
h
= hoop stress;
m
= membrane stress;
M
a
= Folias factor
=
2
1
2
26 . 0 1
,
`
.
|
]
]
]
+
Rw
L
; (9)
L = defect axial length;
w = pipe wall thickness;
Y
m
= normalised stress intensity factor (pure membrane stress)
CORRIB FIELD DEVELOPMENT PROJECT
Onshore Pipeline Quantified Risk Assessment
g:\052377 - enterprisecorrib additional work\ctr 140 - qra update\02-f-3-835- rev 05 - qra.doc Page 52 of 58
=
4 3 2
4 . 30 7 . 21 6 . 10 23 . 0 12 . 1
,
`
.
|
+
,
`
.
|
,
`
.
|
+
,
`
.
|
w
a
w
a
w
a
w
a
; (10)
Y
b
= normalised stress intensity factor (pure bending stress)
=
4 3 2
0 . 14 1 . 13 32 . 7 39 . 1 12 . 1
,
`
.
|
+
,
`
.
|
,
`
.
|
+
,
`
.
|
w
a
w
a
w
a
w
a
; (11)
K
1c
= pipe material fracture toughness
= ( )
5 . 0
6 . 17
1000
]
]
]
v
C
A
E
; (12)
A = cross-sectional area of Charpy specimen;
E = Youngs modulus of pipe steel;
C
v
= Charpy energy.
The limiting value for K
r
is
crit
r
K and failure occurs when the following condition is reached
crit
r r
K K (13)
A failure assessment diagram (FAD) can be constructed by plotting
crit
r
K against S
r
(see
Figure 5.1).
2
1
2
2
sec ln
8
]
]
]
,
`
.
|
r r
crit
r
S S K
for S
r
< 1 (14)
0
crit
r
K for S
r
1 (15)
The failure of a gouge within a dent (i.e.
crit
r r
K K ) leads to the creation of a through wall
defect. The resulting through wall defect will result in a failure by either a leak or a rupture,
depending on the length of the gouge. If the length of the through wall defect is greater than a
critical length L
c
given by equation (16) [Ref. 32] then a rupture will occur. If the length is less
than this critical length, the through-wall defect will be stable and the failure will lead to a leak.
,
`
.
|
]
]
]
]
,
`
.
|
4 . 0
1
15 . 1
2
Rw
L
y
h
c
(16)
CORRIB FIELD DEVELOPMENT PROJECT
Onshore Pipeline Quantified Risk Assessment
g:\052377 - enterprisecorrib additional work\ctr 140 - qra update\02-f-3-835- rev 05 - qra.doc Page 53 of 58
FAIL
SAFE
0
0.2
0.4
0.6
0.8
1
1.2
0 0.2 0.4 0.6 0.8 1 1.2
S
r
K
r
c
r
i
t
Figure A1 - PD6493 Level 2 Failure Assessment Diagram (FAD)
Dent depth cannot be based on historical damage statistics because it is a function of pipeline
geometry, steel properties and operating pressure. The dent depth is therefore derived from
the estimated impact force using the equation given by Corder and Chatain
[Ref. 33] as
follows:
( )
38 . 2
4
1
80
15 . 1
4 . 1
49 . 0
,
`
.
|
,
`
.
|
+
w
R p
w
F
D
y
y
op
(17)
where D = dent depth
F = impact force
w = wall thickness
p
op
= internal pressure
R = mean radius
y
= yield stress
The impact force is estimated using the correlation with excavator mass that was described in
the sub-section relating to punctures.
CORRIB FIELD DEVELOPMENT PROJECT
Onshore Pipeline Quantified Risk Assessment
g:\052377 - enterprisecorrib additional work\ctr 140 - qra update\02-f-3-835- rev 05 - qra.doc Page 54 of 58
For the dent/gouge limit state function, the following variables need to be described by means
of probability distributions :
Wall thickness
Yield stress
Charpy energy
Gouge depth
Gouge length
Excavator mass
Internal pressure
In addition allowance has to be made for modelling uncertainty. This is done by multiplying
crit
r
K by a factor representing the modelling uncertainty, and this factor is itself described by a
probability distribution developed from test data.
A1.4 Solution of the Limit State Models
For the limit states considered in this report, the reliability calculations have been performed
using the Monte Carlo simulation technique. The method involves sampling at random from
the probability distributions representing the parameters in the limit state function in a large
number of numerical experiments. Each sample is used to check the limit state function. If the
limit state function is violated, the pipeline has failed. The numerical experiment is repeated
many times, each with a randomly chosen set of variables. If N iterations are undertaken and
n is the number of failures, the failure probability per damage event is given by:
N
n
P
f
(18)
The accuracy of the Monte Carlo simulation depends on the number of iterations in the
simulation. The method is rigorous and any desired accuracy can be obtained by performing
enough iterations.
Note that P
f
is the failure probability per damage event. This must then be multiplied by the
probability of the damage event occurring to obtain the overall failure probability (e.g. per km-
year).
For the external interference limit state, the probability of a puncture, or of a dent/gouge that
results in a leak or a rupture, is evaluated as follows:
i
p
p
P
N
n
P (19)
i
l
l
P
N
n
P (20)
i
r
r
P
N
n
P (21)
CORRIB FIELD DEVELOPMENT PROJECT
Onshore Pipeline Quantified Risk Assessment
g:\052377 - enterprisecorrib additional work\ctr 140 - qra update\02-f-3-835- rev 05 - qra.doc Page 55 of 58
where P
p
= probability of a puncture (per km-year);
P
l
= probability of a gouge/dent resulting in a leak (per km-year);
P
r
= probability of a gouge/dent resulting in a rupture (per km-year);
n
p
= number of iterations resulting in F>r;
n
l
= number of iterations resulting in
crit
r r
K K
and L < L
c
;
n
r
= number of iterations resulting in
crit
r r
K K
and L L
c
;
N = total number of iterations;
P
i
= probability of an external interference incident (per km-year).
A1.5 Distributions used in the Reliability Calculations
Table A1 summarises the probabilistic description of the parameters in the reliability
calculations. The steel properties are typical for modern pipeline steels and have been based
on mill certificates from previous projects.
The gouge lengths and depths are those that have been reported by British Gas Technology
as experienced in damage events to the national gas transmission system. The data indicate
that the depth and the length are uncorrelated and independent of the pipeline properties.
In Ref. 18, British Gas Technology indicated a frequency of external interference events of
1.86 x 10
-3
per km year as applying to typical 36-inch diameter national transmission system
pipelines. This has been conservatively assumed for the study location.
The distribution of excavator weights given in Figure A2 is from North American data [Ref.
34]. This distribution was considered slightly over-conservative for the area in which the
pipeline is located as heavy excavating equipment is unlikely to be used. In the model used,
therefore, the distribution was truncated to impose a maximum limit of 40 tonnes. This is
still considered to be a conservative assumption for the project area.
CORRIB FIELD DEVELOPMENT PROJECT
Onshore Pipeline Quantified Risk Assessment
g:\052377 - enterprisecorrib additional work\ctr 140 - qra update\02-f-3-835- rev 05 - qra.doc Page 56 of 58
Table A1 - Parameter distributions assumed in the models
Variable Distribution Type Mean COV Reference
Operating
pressure
Constant Design pressure 0
Pipe radius Constant Nominal value 0 Varies little
Wall thickness Normal Nominal (t
nom
) +/- 1.0mm Fab. Tol.
Pipe yield stress Lognormal 1.08 SMYS 0.04 Typical
Pipe ultimate
tensile strength
Normal 1.06 SMTS 0.03 Typical
Pipe Charpy
energy
Lognormal 130 J 0.40 Typical
Excavator tooth
length
Normal 62.5 mm 0.14 Ref. 28
Excavator tooth
width
Normal 4.5 mm 0.167 Ref. 28
Additive puncture
model error
Normal 0.883 kN 30.24 Ref. 29
Gouge depth Weibull Shape parameter = 0.73
Scale parameter = 0.98 mm
Ref. 18
Gouge length Offset logistic Shape parameter = 0.043
Scale parameter = 24.84 mm
Offset parameter = 30.13 mm
Ref. 18
Excavator mass Histogram See Figure B2 Ref. 34
Dent/gouge
multiplicative
model error
Lognormal 1.0 0.30 J P Kenny
assessment
using published
test data
CORRIB FIELD DEVELOPMENT PROJECT
Onshore Pipeline Quantified Risk Assessment
g:\052377 - enterprisecorrib additional work\ctr 140 - qra update\02-f-3-835- rev 05 - qra.doc Page 57 of 58
Figure A2 - Distribution of excavator weight according to Ref. 34 (US)
A2 References (Appendix A)
27 Spiekhout, J, 1995. A New Design Philosophy for Gas Transmission Pipelines Designing for
Gouge Resistance and Puncture Resistance
28 Driver R G and Zimmerman T J E. A Limit State Approach to the Design of Pipelines for
Mechanical Damage. Proceedings of the 17t h International Offshore and Arctic Engineering
Conference, OMAE 98-1017, Lisbon, Portugal, July 1998
29 Joint Industry Project to Develop Guidance for Limit State, Reliability and Risk Based Design
and Assessment of Onshore Pipelines, C-FER 2000
30 Example of the Application of Limit State, Reliability and Risk Based Design to the Uprating of
an Onshore Pipeline, A M Edwards, R J Espiner and A Francis, BG Technology, Report No.
R3125, Issue No.1, 6 August 1999, Commercial in confidence to Joint Industry Project
Sponsors
31 BSI, PD 6493 : 1991 Guidance on Methods for Assessing the Acceptability of Flaws in
Fusion Welded Structures
32 R J Espiner, British Gas Technology, private communication
33 Corder I and Chatain P, EPRG Recommendations for the Assessment of the Resistance of
Pipelines to External Damage, 1995
34 Reliability Based Planning of Inspection and Maintenance of Pipeline Integrity, Final Report
PR-244-9517, C-FER, published by American Gas Association, 1997
CORRIB FIELD DEVELOPMENT PROJECT
Onshore Pipeline Quantified Risk Assessment
g:\052377 - enterprisecorrib additional work\ctr 140 - qra update\02-f-3-835- rev 05 - qra.doc Page 58 of 58
This page intentionally left blank