ch3 PDF
ch3 PDF
System Hydraulics
and Design
This chapter discusses the fundamentals of liquid pipeline hydraulics and the design
and operation (Chapter 5) of hydrocarbon liquid pipeline systems from a hydraulics
point of view. Pipeline system design is mainly concerned with line sizing, equipment
sizing and location, and flow capacity; while system operation is concerned with pipeline system or facility start-up and shut-down, product receipt and delivery, flow rate
changes, emergency shut-down, equipment failure, etc.
A proper pipeline system design requires a system approach taking into account
the following design disciplines:
Hydraulics
Mechanical design
Geo-technical design
Operations and maintenance design
These disciplines are closely interrelated because any decisions or changes in one
area of design directly affect or limit the options in another area. Through the hydraulic design, the pipeline route, pipe size, operating pressure and temperature and the
number of pump stations are determined. From a hydraulic design, mechanical designs
can be developed to meet the criteria of the design basis. The mechanical design is
dictated by the relevant codes and standards, resulting in pipe material selection and
specifications as well as burial depth requirements. Geo-technical design addresses
surface loads, water crossings, buoyancy control and geo-hazard management, which
can significantly affect the cost and safety, if the pipeline route traverses challenging
environments. The operation and maintenance consideration includes the necessary
control systems to operate the system within its design parameters, taking account of
the operating tasks while maintaining the functional integrity of the system.
The scope of this chapter includes the governing principles and equations of liquid pipeline hydraulics and their solutions in steady states. The design of any pipeline
system is based on various design factors such as flow profile over time and operating
pressures.
(3 1)
where
A = Cross sectional area of the pipe
The cross sectional area can change due to the changes in pressure and temperature:
A = A0 1 + cP ( P - P0 ) + cT (T - T0 )
(3 2)
where the subscript zero refers to base or standard conditions. cT is the coefficient for
thermal expansion of the pipe material and its effect on transients is negligibly small.
CP has a large effect on the acoustic speed of a pressure wave and is defined as:
cP =
1 D
1 - m2
E w
(3 3)
where
E = Youngs modulus of elasticity of the pipe
w = Pipe wall thickness
m = Poisson's ratio
The first term in the continuity equation represents the change of mass in a pipe
segment. It is often called line pack change. The line pack can be increased or decreased due to pressure and temperature changes. The line pack change is useful for
gas pipeline operation. It should be distinguished from the line fill volume, which is the
quantity of fluid contained in a pipeline. It is also a useful quantity for batch pipeline
operation. The second term represents the difference between mass flow into and out
of the pipe segment.
3.1.1.2 Momentum Equation
The momentum equation describes the motion of the fluid in the pipeline. It requires
fluid density and viscosity in addition to the pressures and flows. Applying Newtons
second law of motion to a fluid element together with the Darcy-Weisbach frictional
V
V P
h f rV | V |
+ rV
+
+ rg
+
= 0
t
x x
x
2D
(3 4)
where
r = Density of the fluid
V = Velocity of the fluid
P = Pressure on the fluid
h = Elevation of the pipe
g = Gravitational constant
f = Darcy-Weisbach friction factor
D = Inside diameter of the pipe
x = Distance along the pipe
t = Time
The first term is a force due to acceleration, and the second term a force due to
kinetic energy. These two terms are related to inertial force. The third term is a force
due to pressure difference between two points in a pipe segment. The fourth term is
a gravitational force, and the last term is a frictional force on the pipe wall, opposing
the flow.
The Darcy-Weisbach equation is used to calculate the pressure drop due to the
friction of fluid flow against the pipe wall. The friction pressure drop is linearly proportional to the fluid density and the friction factor, squarely proportional to fluid velocity, and inversely proportional to the pipe diameter. The friction pressure drop can be
expressed as follows:
f rV | V | 8 f rQ 2
= 2 5
2D
p D
(3 5)
In terms of flow rate, the frictional pressure drop is proportional to the square of
the flow rate and inversely proportional to the fifth power of the pipe diameter. Since
the frictional pressure drop and thus pipeline flow capacity depends highly on pipe
diameter, it is the most significant design parameter. The friction factor is related to
the energy losses resulting from fluid flow. It is a function of the Reynolds number and
pipe roughness. Depending on the Reynolds number, the type of pipe flow is classified
into three flow regimes: laminar flow, critical flow, and turbulent flow. Turbulent flow
can be further divided into partially turbulent, where the smooth pipe law applies, and
fully turbulent, where the rough pipe law applies.
The Reynolds number is dimensionless and the ratio of inertial forces to viscous
forces. It is defined by
Re =
| V | rD | V | D
=
m
(3 6)
where
m = dynamic viscosity (kg/m s)
n = m/r = kinematic viscosity (m2/s)
r = fluid density (kg/m3)
V = flow velocity, m/s
D = pipe inside diameter, m
k
1
2.51
= - 2 log
+
for Re 4,000
f
3.7D Re f
(3 7)
where k is the pipe roughness, D the pipe inside diameter, and Re is the Reynolds
number. For laminar flow, the friction factor is:
f =
64
for Re 2400
Re
(3 8)
The critical flow regime is defined between 2,400 < Re < 4,000, in which the flow
is unstable. Laminar flow is independent of pipe roughness, while partially turbulent
flow is dependent on Reynolds number and pipe roughness, and fully turbulent flow is
dependent only on relative roughness being independent of Reynolds number.
The Moody diagram, shown in Figure 3-1, relates the friction factor in terms of
Reynolds number and relative roughness.
The Colebrook-White equation is not easily solvable without a computer because
the friction factor appears on both right and left sides of the correlation. To facilitate an
(3 9)
for 106 < k/D < 102 and 5000 < Re < 108
Churchills formula
f = 8[(8/Re)12 + (A + B) 1.5]1/12
(3 10)
These equations correlate closely with friction factors on the Moody diagram. The
Fanning friction factor ff is occasionally used and related to the Darcy friction factor
as follows:
f = 4 ff
(3 11)
Other pressure drop equations, such as the Shell-MIT equations and Hazen illiams, are sometimes found in the literature. Since the Darcy-Weisbach equation
W
with the associated Darcy friction factor is most widely used in the petroleum pipeline
industry, it will be used throughout this book.
Most liquid hydrocarbon pipelines are operated in partially turbulent flow regimes,
with the exception of ethylene and ethane flow which may be in a fully turbulent regime and heavy crude which may be in a laminar flow regime.
3.1.1.3 Energy Equation
The energy equation accounts for the total energy of the fluid in and around the pipeline, requiring information regarding the flows, pressures, and fluid temperatures together with fluid properties and environmental variables, such as conductivity and
ground temperature.
4wrpCp T
T
r v v A
+ rvCv
+T
+
+
rCv +
D t
x
T x A x
f rv 2 | v | 4k dT
+
=0
2D
D dz
(3 12)
where
Cv = Specific heat of the fluid at constant volume
T = Temperature of the fluid
rp = Density of the pipe material
Cp = Heat capacity of the pipe material
k = Heat transfer coefficient
z = Distance from the pipe to its surroundings
The first term is the temperature change over time, the second is the rate of temperature change due to the net convection of fluid energy into the fluid element. The
third term describes the change rate due to expansion/compression of the fluid including the Joule-Thomson effect. The fourth term represents the heat flow to, or from,
r(P,T)V(P,T) = r(P0,T0)V(P0,T0)
or
r(P,T)Q(P,T) = r(P0,T0)Q(P0,T0)
(3 14)
This relationship is the basis of converting volume or flow rate from one pressure
and temperature condition to another including volume correction to base conditions.
Its application is illustrated with the following base design example (this example will
be extended further to a realistic design case):
Example: Base Case
A crude oil pipeline from CE to QU is 200 km long and is 20 in nominal diameter, with a
0.281 wall thickness. It is constructed of 5LX-56 electric resistance welded steel pipe. At
the injection point, crude oil of 32API gravity and ambient pressure enters the pipeline at
an initial flow rate of 18,000 m3/d at 15C. The average operating pressure and temperature
are 4000 kPag and 4C. Calculate the flow rate at the operating conditions.
Figure 3-2 illustrates this pipeline configuration, which will be used for the subsequent example problems in this chapter. CE is the initial injection station, QU is the
final delivery station, and TO a side stream delivery point.
Solution:
It is assumed that the API correction equation or equation of state (Refer to Chapter 2)
is applicable to convert the density at the base condition to the density at the operating
pressure and temperature.
Step 1. To determine the flow rate at the operating conditions, the crude density at
the same conditions should be determined.
The density equivalent to 32API gravity is calculated by applying the API
gravity and the specific gravity relationship, thus the specific gravity is g = 141.5/
(131.5 + API) = 0.8654, and the density is
dP
dh f rV | V |
+ rg
+
=0
dx
dx
2D
(3 15)
It can be assumed that the liquid density and velocity are constant between two
points along the pipeline. This assumption is valid as long as the distance between two
points is not long. Therefore, the pressure-flow equation can be obtained by integrating
the above steady-state momentum equation:
(3 16)
The left hand side is the pressure at the downstream point. The first term on the
right hand side is the pressure at the upstream point, the second the static pressure
or elevation head, and the third the friction head. The total pressure requirement in a
pipeline system consists of the following components:
Pressure changes due to elevation changes, depending only on the product density and difference between the elevations between two points on the pipeline;
Friction pressure drop due to flow rate or velocity, fluid density and viscosity,
and pipe diameter;
Pressure changes due to changes in pipe diameter and subsequent changes in
flow velocity.
For a given flow rate, the above pressure-flow equation allows us to calculate the
downstream pressure if the upstream pressure is known, and the upstream pressure if
the downstream pressure is known. Also, the flow rate can be calculated if the upstream
and downstream pressures are known.
If the static pressure term is set aside, the above equation can be arranged as
(3 17)
1880 J/kgC
9.5 cSt at 35C and 43.5 cSt and 4C
0C
Solution:
It is assumed that the viscosity of this product is Newtonian and that the density and
viscosity depend on temperature. The fluid density and viscosity are calculated at the
starting point temperature in the segment between two profile points. Let the inlet pressure be 8605 kPag, the same as for the isothermal case.
Step 1. Since the density and viscosity change with temperature, the temperature
relationships of density and viscosity need to be established to calculate these quantities as the temperature profile is calculated.
Step 2. Calculate the density and viscosity at the inlet conditions; r(35) = 851.0 kg/m3
and (16) = 9.5 cSt.
Step 3. Use the inlet temperature of the first segment to calculate the friction factor
of 0.0201 and the frictional pressure drop of 508 kPa.
Step 4. Calculate the temperature at the downstream point of the first segment.
The temperature increase due to the frictional pressure drop is 0.32C
To calculate the temperature drop due to conduction, the following values are
calculated iteratively:
the heat transfer coefficient, 0.324 W/m2C;
the log mean temperature, 34.1C;
the temperature drop at downstream temperature 2.1C;
hence the downstream temperature is 35 + 0.32 2.10 = 33.2C.
Step 5. Calculate the pressure and temperature at the other profile points by
repeating the above steps.
KMP
(km)
0
20
30
60
80
90
Elevation
(m)
Pressure
(kPag)
Temp (C)
KMP
(km)
Elevation
(m)
Pressure
(kPag)
Temp
(C)
30
55
45
30
70
100
8605
7889
7714
7060
6195
5674
35.0
33.2
32.4
30.0
28.6
27.9
100
130
150
160
180
200
130
100
60
110
150
130
5152
4586
4368
3666
2765
2364
27.3
25.4
24.3
23.9
22.8
21.9
It is expected that the total pressure requirement is lower than the pressure requirement under the isothermal assumption, because the operating temperature would
be higher and thus the values of density and viscosity are lower. Indeed, the delivery
pressure turns out to be much higher than the delivery pressure for the isothermal
case, and so the total pressure requirement is less by 2014 kPa. It is concluded that the
temperature effects have to be included in hydraulic calculations if the liquid injection
temperature is much higher than the ground temperature.
Note that the frictional pressure drop remains the same even though the elevation
changes.
Step 1. Use the same pressure gradient as obtained in the previous example.
Step 2. Calculate the pressures at the above profile points by adding the static
pressure difference to the frictional pressure drop.
KMP (km)
Elevation (m)
Pressure (kPag)
KMP (km)
30
55
45
30
70
100
8605
7650
7366
6384
5302
4676
100
130
150
160
180
200
0
20
30
60
80
90
Elevation (m)
Pressure (kPag)
130
100
60
110
150
130
4049
3196
2799
2001
919
350
Elevation (m)
KMP (km)
Elevation (m)
30
55
45
30
70
100
100
130
150
160
180
200
130
100
60
110
150
130
Liquid
Heat generation
Heat out
Insulation
Heat out
Some pipelines may be partially or wholly installed aboveground to save construction or maintenance cost. However, transmission pipelines are generally buried
in order to:
Minimize land use disturbance,
Provide longitudinal restraint along pipeline length,
Protect pipe from possible pipe material fatigue due to stress changes caused by
fluctuations in ambient temperature,
Minimize effects of changes in ambient temperature on fluid viscosity and
density,
Protect pipe from intentional or accidental damage, and
Use the pipeline right of way.
The temperature calculation from the energy equation is not simple even under a
steady-state condition. The steady-state energy equation can be derived by balancing
heat entering and leaving a pipe section, heat transferred from/to the pipe section, to/
from surrounding soil or ambient, and heat from friction. The heat balancing mechanism can be shown in Figure 3-5, and the heat balance is expressed as:
(3 18)
where
Hin = Heat entering a pipe section (w)
Hout = Heat leaving a pipe section (w)
Hcon = Heat transferred from/to the pipe section to/from surrounding soil or
ambient (w)
Hf = Heat from friction (w).
Described below is a temperature calculation procedure. Another method for calculating temperature profile is presented in Addendum 3.1, which includes a temperature calculation method for above-ground pipelines.
1. Assuming that the specific heat of the fluid remains constant at the entering
and leaving conditions, the heat entering and leaving a pipe section can be
expressed in terms of temperatures and engineering quantities as follows:
(3 19)
where
r = liquid density (kg/m3)
Q = flow rate (m3/hr)
Cp = specific heat of liquid, kJ/kg/C
Ti = temperature of liquid entering the pipe section, C
To = temperature of liquid leaving the pipe section, C
2. As the liquids flow through the pipe, the pipe pressure drops by friction, liquid
flows undergo an isenthalpic process, and as a result the pressure dissipated by friction becomes heat in the flowing fluid. The temperature of liquids rises in frictional
heating due to their volumetric properties as they are expanded in an isenthalpic
process. The effect of friction heating generally increases with flow rate, viscosity,
insulation, and line length. For large diameter pipelines and high flow rates, heat
generated by friction loss should be included in the temperature profile.
Heat of friction should be considered at high flow rates in large pipelines
to ensure that overheating does not occur. Pump stations operating on flow
control may experience increasing or decreasing discharge pressures as the
temperature of the fluid in the pipeline rises or falls after leaving the pump station. As the temperature increases, the fluid expands. As expansion continues
in the pipeline, the local pressure and volumetric flow rate increases. The heat
generated by frictional pressure drop is expressed as
Hw = q DPf = 0.278Q (DPf/Dx) L
(3 20)
where
Hw
= frictional heating, w
q
= liquid flow rate, m3/sec
DPf = frictional pressure drop, Pa
Q
= liquid flow rate, m3/hr
DPf/Dx = frictional pressure gradient, kPa/km
L
= length of the pipe section, km
3. Even though ground temperature along the pipeline is not normally measured
on a daily basis, it is an important parameter for designing a pipeline system.
Significant temperature changes can occur due to heat transfer through conduction between the liquid and surrounding soil. In describing the flow of heat
from pipeline to ground, Fouriers law of heat conduction is applied by taking
into account the heat transfer through pipe, insulation, and soil. The conduction
heat transfer can be expressed as:
Hcon = U A DTm = 2p DT L U DTm
where
U
A
DT
L
DTm
Tg
(3 21)
(3 22)
where
Rif = thermal resistance due to the boundary layer that builds up on the inside of
the pipe wall (m2C/w)
Rp = thermal resistance of the pipe wall (m2C/w)
Rins = thermal resistance of insulation (m2C/w)
Rs = thermal resistance due to the surrounding medium (m2C/w)
However, the heat transfer effects due to the boundary layer and pipe wall are
much smaller than those due to surrounding soil or insulation. Therefore, these two
terms are usually ignored, and only the last two terms are considered in the overall heat
transfer calculation.
Pipelines are not frequently insulated unless the fluid viscosity is so high that it can
be significantly reduced by heating the fluid. If the fluid such as heavy crude is heated,
certain parts of the pipeline are insulated. For an insulated pipe, the heat resistance can
be determined by,
Rins = (DT/kins) Ln(DT/D)
(3 23)
Ground
Insulation
Liquid
film
Steel
Pipe
Corrosion
coating
Outer
Jacket
(3 24)
The thermal conductivity is a measure of how easily heat conducts through the
material. It appears in Fourier's law of heat conduction. Generally, the thermal conductivity can be nearly constant over the temperature range normally encountered in
pipelines. Thermal conductivity is measured in units of W/(mC) (Table 3-2).
Certain portions of a pipeline may run above-ground, even for heated liquids, in
order to reduce the construction and other costs. Above-ground pipelines are usually
insulated. If the above-ground pipe length is long enough to affect the temperature profile, the heat transfer between the liquid and ambient air needs to be calculated. Since
the ambient air conditions can change significantly in a short time, their effects need
to be evaluated for design based on the average and worst conditions but are difficult
to assess for operation.
In heat transfer calculations, the log mean temperature difference between the
liquid in a pipe section and the surrounding soil is often used. This is because the fluid
Table 3-2. Thermal conductivity
Substance
0.450.70
0.851.05
1.902.25
0.350.50
0.700.85
1.051.55
0.020.05
(3 25)
where
Ti = temperature of liquid entering the pipe section (C)
T0 = temperature of liquid leaving the pipe section (C)
Tg = ground or surrounding medium temperature (C)
Therefore, the log mean temperature difference is determined by the equation:
DTm = Tm Tg = (Ti T0)/Ln[(Ti Tg)/(T0 Tg)]
(3 26)
Note that a log mean temperature is similar to a simple arithmetic average temperature for short pipe lengths over which the temperature is calculated, and that both
the log mean temperature and arithmetic average temperature contain the downstream
temperature that has to be calculated in the temperature profile calculation. Therefore,
an iterative technique is used to calculate either the log mean or arithmetic average
temperature and this can be easily implemented in software. A manual calculation
can also generate a reasonable temperature profile to the known upstream temperature
instead of using the log mean temperature.
Combining the above equations for temperature, we have
T0 = Ti + Pf /(rCp) Hcon/(rQCp)
(3 27)
Temperature
where
T0 = Outlet temperature (C)
Ti = Inlet temperature (C)
DPf = frictional pressure drop, Pa
r = density (kg/m3)
Q = flow rate (m3/sec)
Cp = specific heat (J/kg C)
T0
Temperature Profile
Ground Temperature
TG
Pipe Length
(3 28)
where
A = pipe surface area
dx = differential in distance
This equation can be integrated to obtain
Tx = Tg + [T0 Tg] exp[ (2p UDTX)/(rQCp)]
(3 29)
This equation shows that the temperature profile decays exponentially and that the
delivery temperature drops closer to the ground temperature. If the frictional heating
term is included, the overall temperature profile is elevated. The temperature equation
indicates that, assuming the ambient temperature is lower than the liquid temperature,
the liquid cools faster and its viscosity increases as flow rate decreases.
Note that the effect of friction heating increases with flow rates and viscosity because the frictional pressure drop is high. Therefore, a frictional heating term should be
included for the case of high flow rates and/or high viscosity liquid. Also, the calculation of a temperature profile is so complex and prone to error that it is beneficial to use
a computer software package to obtain quick and accurate results. Temperature-related
problems are more severe for larger pipelines because the conduction heat loss is proportional to pipeline surface area.
The surrounding environment is the key factor in the overall heat transfer coefficient, which is most critical in calculating the temperature profile along the pipeline.
Table 3-3 shows the range of overall heat transfer coefficients for an on-shore pipelines surrounding environment [14].
Example: Base Case Extension 3
The previous base case is extended to include the temperature profile by removing the
isothermal assumption. Assuming that the pipeline is not insulated, calculate the pressure and temperature profiles using the following data:
35C
4C
1.2 m
0.5 W/mC
U Value (W/m2 C)
0.853.69
0.280.85
1.704.54
0.571.14
3.978.52
0.571.15
1880 J/kgC
9.5 cSt at 35C and 43.5 cSt and 4C
0C
Solution:
It is assumed that the viscosity of this product is Newtonian and that the density and
viscosity depend on temperature. The fluid density and viscosity are calculated at the
starting point temperature in the segment between two profile points. Let the inlet pressure be 8605 kPag, the same as for the isothermal case.
Step 1. Since the density and viscosity change with temperature, the temperature
relationships of density and viscosity need to be established to calculate these quantities as the temperature profile is calculated.
Step 2. Calculate the density and viscosity at the inlet conditions; r(35) = 851.0 kg/m3
and (16) = 9.5 cSt.
Step 3. Use the inlet temperature of the first segment to calculate the friction factor
of 0.0201 and the frictional pressure drop of 508 kPa.
Step 4. Calculate the temperature at the downstream point of the first segment.
The temperature increase due to the frictional pressure drop is 0.32C
To calculate the temperature drop due to conduction, the following values are
calculated iteratively:
the heat transfer coefficient, 0.324 W/m2C;
the log mean temperature, 34.1C;
the temperature drop at downstream temperature 2.1C;
hence the downstream temperature is 35 + 0.32 2.10 = 33.2C.
Step 5. Calculate the pressure and temperature at the other profile points by
repeating the above steps.
KMP
(km)
0
20
30
60
80
90
Elevation
(m)
Pressure
(kPag)
Temp (C)
KMP
(km)
Elevation
(m)
Pressure
(kPag)
Temp
(C)
30
55
45
30
70
100
8605
7889
7714
7060
6195
5674
35.0
33.2
32.4
30.0
28.6
27.9
100
130
150
160
180
200
130
100
60
110
150
130
5152
4586
4368
3666
2765
2364
27.3
25.4
24.3
23.9
22.8
21.9
It is expected that the total pressure requirement is lower than the pressure requirement under the isothermal assumption, because the operating temperature would
be higher and thus the values of density and viscosity are lower. Indeed, the delivery
pressure turns out to be much higher than the delivery pressure for the isothermal
case, and so the total pressure requirement is less by 2014 kPa. It is concluded that the
temperature effects have to be included in hydraulic calculations if the liquid injection
temperature is much higher than the ground temperature.
Pipe grade
Pipe size
(in/mm)
Design pressure
(kPag)
Discharge
pressure at inlet
point (kPag)
Discharge pressure
at intermediate
station (kPag)
X65
X70
X70
22/558.8
20/508.0
22/558.8
8246
9765
8880
5066
7870
5066
5066
7870
5066
The capital cost due to the extra pumping power requirement for the 20
pipe is higher than the cost for the 22 pipe size, while the pipe cost for X70
with 20 diameter may cost less than the other two options. The extra capital
cost for the 22 line is more than 20% and is incurred by the extra pipe material
and construction expenses. However, the extra capital cost of the 22 diameter
pipe might be partly compensated by lower unit pumping cost. Assuming that
the annualized cost for the 20 pipe case is lowest, it is selected as the base
design.
The facilities such as the initiating pump station for the selected base design
would be designed to accommodate the capacity until the capacity increases
in the 10th year. In the 10th year, the additional facility increases include the
pumping capacity at the inlet point for the additional flow and an intermediate
pump station with the pumping capacity of 30,000 m3/d.
4. Develop alternative design cases and perform comparative studies against the
base design
Alternative 1: This alternative design is to use a pipe wall thickness larger
than 0.281 in order to increase the design pressure slightly higher than the
total pressure requirement. No intermediate pump station is required if the design pressure is slightly higher than the total pressure requirement. Note that
the required total pressure will be increased due to slightly smaller inside pipe
diameter.
The design pressure for the X70 22 pipe is lower than the total pressure requirement, which in turn is lower than the maximum operating pressure range.
The next largest nominal wall thickness is 0.312 or 7.92 mm, and its design
pressure is 9857 kPag or 1430 psig, but the required total pressure is 9914 kPag
for a flow rate of 30,000 m3/d. Therefore, the wall thickness is not sufficient to
meet the total pressure requirement without an intermediate pump station.
The next largest nominal wall thickness is 0.344 or 8.74 mm, which
can allow the design pressure to increase up to 10,870 kPag. For this wall
thickness, the required total pressure turns out to be 10,051 kPag. Since this
design pressure is higher than the required total pressure, no intermediate
pump station is required for the flow rate expected beyond the 10th year,
V = Q/A
(3 30)
where:
V = Liquid velocity
Q = Flow rate
A = pipe cross sectional area
It may be noted that there are a number of situations where selecting a pipe size based
on the optimum fluid velocity is not appropriate and a detailed analysis will be required
The pressure gradient or pressure drop per unit length of pipe is an important
measure for designing a safe and economic pipeline system. Since the liquid velocity
Suggested velocity
0.9 m/s to 2.0 m/s
1.5 m/s to 3.5 m/s
1.8 m/s to 4.0 m/s
2.4 m/s to 4.5 m/s
These figures approximate only but generally provide an economic pipeline and
piping design.
3.2.2.3.2 High Viscosity Fluids As the liquid viscosity increases above 10 cSt, the
suggested velocities are lower than those listed above. However, for high viscosity
liquids (i.e., these with viscosities approaching 1000 cSt and higher), pipeline and piping design would not be based purely on economic factors. For high viscosity liquids,
keeping the pressure drop to within acceptable limits is likely to be the key. It may be
noted that there are a number of situations where selecting a pipe size based on the
optimum fluid velocity is not appropriate and a detailed analysis will be required.
No pipeline systems can operate continuously for a full calendar year due to operational restrictions such as system maintenance or other reduced capacity operations.
The average daily flow is obtained by dividing the annual throughput by 365 (yearly
calendar days), and the actual maximum daily flow by the actual number of operating days. The ratio of operating days to calendar days is called load factor, so the load
factor can be defined as the average daily flow divided by the actual maximum daily
flow. Normally, the maximum daily flow is used for design in order to compensate for
the downtime. In the design procedure, a load factor of up to 95% is used for a simple
pipeline, while it may be as low as 85% for more complex systems or pipelines operated with expected large flow variations.
The minimum flow rate has to be defined for system design and operation, because
all equipment has maximum and minimum operational limits in capacity and efficiency.
For example, a pump can only operate within a flow bound between the maximum and
minimum capacity. In a highly mountainous terrain, slack flow conditions may occur
at low flows so that extra equipment specifications are required to operate the pipeline
safely. Refer to Section 3.3.3 for a detailed discussion of slack flow conditions.
Choice of operating pressures directly affects pipeline safety and operating requirements. The requirements include shipping capacity and volume demands, location
Pipe OD (mm/in)
Standard wall
thickness (mm/in.)
Weight (tons/km)
168.3 / 6.625
219.1 / 8.625
273.1 / 10.752
323.9 / 12.752
355.6 / 14
457.2 / 18
508.0 / 20
609.6 / 24
812.8 / 32
914.4 / 36
1066.8 / 42
1219.2 / 48
7.11 / 0.280
8.19 / 0.322
9.27 / 0.365
9.53 / 0.375
9.53 / 0.375
9.53 / 0.375
9.53 / 0.375
9.53 / 0.375
9.53 / 0.375
9.53 / 0.375
9.53 / 0.375
9.53 / 0.375
28.27
42.54
60.32
73.84
81.31
105.18
117.11
140.98
188.72
212.59
248.40
284.20
(3 31)
where
S = specified minimum yield strength (SMYS) of pipe, kPag, or psig
t = pipe wall thickness, mm or in.
D0 = outside pipe diameter, mm or in.
F = design factor or safety factor
L = location factor (L = 1 for liquid pipelines)
J = joint factor (to reflect the method of pipe joining generally taken to have a
value of 1)
T = temperature derating factor, to account for the effect of higher temperatures
on yield stress
The SMYS is a standard measure of the specified minimum yield strength for steel
pipe. Standards that are frequently used by the pipeline industry are API 5L: Specifications for Line Pipe, which includes API 5LX and 5LS. API 5LX specifies various
strength grades, ranging from Grade B, rated at 42,000 psig (289 MPag) to Grade
X120, rated at 120,000 psig (827 MPag), where the Grade X120 refers to the SMYS in
1000 psi. Pipes are manufactured to these specifications.
ASME B31.4 does not define the location factor. The design factor, F, specified
in ASME B31.4 is 0.72 for liquid pipelines regardless of the location of the pipeline,
while other codes such as CSA Z662 define the design factor differently depending on
the locations. The joint factor is 1 for all types of pipe manufactured to 5LX and 5LS
specifications. The temperature derating factor is generally taken as 1 for transmission
pipelines, because transmission lines are seldom operated beyond the temperature derating range. Several mechanical design aspects are discussed in the next chapter.
Effective pipe roughness is a pipe parameter that affects frictional pressure drop
and pipeline efficiency. It includes pipe roughness as well as other pressure loss terms
such as bends, welds, fittings, etc. It directly influences the friction factor of the fluid
flow; the larger the pipe roughness, the higher the frictional resistance. To reduce
roughness, pipes are internally coated or cleaned by pigging. Several examples of pipe
conditions and their corresponding roughness are listed in Table 3-5, showing also that
pipe roughness varies with pipe conditions.
3.2.2.4.2 Maximum Allowable Operating Pressure (MAOP) In the design of
pipelines and their components, the design engineer must ensure that the design pressure at any point along the pipeline is lower than or equal to the maximum design
pressure or maximum allowable operating pressure (MAOP). As discussed earlier,
the design pressure is proportional to pipe strength and the MAOP defines the maximum pressure permitted for steady-state pipeline operations which relates to the pipes
ability to withstand internal pressure. The MAOP is the sum of the pressure required
to overcome friction losses, static head pressure, and any required back pressure or
Table 3-5. Pipe roughness
Pipe conditions
Roughness (in.)
Roughness (mm)
0.00050.0008
0.00030.0005
0.00020.0003
0.00180.0020
0.01270.0200
0.00760.0127
0.00510.0076
0.04450.0508
Pump Capacity
Performance curves
Operating ranges (flow, pressure and temperature)
Pump efficiency
Cooler/heater requirement parameters
Station auxiliary equipment requirements and specifications
Energy/Power requirements and specifications
Driver requirements and specifications
Piping requirements and specifications
Mechanical factors
1. Pipe grade, pipe size or diameter, and wall thickness
2. Pipeline route and depth of cover
Capacity factors
1. Operating parameters
2. Station spacing and pumping costs
Reliability and safety factors
1. Valve spacing
2. Other valve-related costs
3.2.2.6.1 Pipe Grade, Size, and Wall Thickness It is critical to optimize the pipe
grade, diameter, and wall thickness to minimize the project cost. The pipe cost is based
on the grade, diameter, and wall thickness. For most pipeline systems, the pipe cost
is the highest material cost. In addition, these three factors have a direct effect on the
cost of installation. Pipeline economics begins with the selection of the pipe material.
Since pipe material for transmission lines is steel, it boils down to the selection of pipe
grade. Higher grade steels are more costly to produce and because of their chemical
composition require specific welding procedures. Nevertheless they do result in thinner pipe wall hence less steel tonnage, lower transportation costs, and reduced amounts
of welding. A case specific study is needed to determine if such steels are the optimal
solution to a given project.
One common economic decision is whether to construct a large line initially, or
put in a smaller line first and parallel it or add pumps at a later time. Once the need for
a pipeline system is recognized, the maximum pipe size is determined such that it can
be economically optimized. The larger the pipe size, the larger the carrying capacity
and the lower unit shipping costs. The pipeline capacity increases approximately by
5/2 power for a fixed pressure drop, but the pipe material cost increases significantly
and construction costs increase almost linearly as the size is increased.
The design pressure is directly proportional to pipe wall thickness for the same
grade and size. The larger the wall thickness for a given pipe size, the higher the design
pressure. The larger the wall thickness, the higher the pipe and construction costs.
Higher grade pipe requires thinner pipe wall for the same design pressure, resulting in
lower steel weight and reduced cost even though higher pipe grade costs more per ton.
Cost savings can also result from reduced construction costs.
3.2.2.6.2 Pipeline Route Both direct and indirect costs due to time delays have to
be taken into account in selecting a pipeline route. As noted earlier the costs of selecting a pipeline route are related to pipeline length, terrain features, intermediate supply
The conceptual design may include hydraulic and economic studies, which result
in overall system and financial requirements.
There are several types of pipeline system design; a new pipeline system, increasing the design capacity of an existing pipeline, and delivery from and/or injection to
other points outside the existing system. The increases in the capacity of the existing
system may require additional pump stations, a parallel line, or replacement of the
existing pipeline with a larger pipe size. Also, the route of the existing pipeline can be
moved due to significant supply/demand changes, or some existing pump stations may
be relocated to other sites to improve the operational efficiency and subsequently to
increase capacity.
After the conceptual design is approved, the pipeline system design is done to
achieve the minimum combined capital and operation costs. In the system planning
phase, the hydraulic and economic evaluation studies are performed in relative detail,
by taking into account the product properties and volumes to be transported, pipeline
route and terrain data, operating temperature ranges and possibly preliminary pressure
ranges, economic and financial data, and other factors such as environmental conditions and restrictions. Described below is a process for performing hydraulic and
economic studies:
1. Gather data
Receive the commitment from shippers for the proposed pipeline
Forecast the supply/demand volumes
Select a preliminary route for the pipeline
Solution:
Refer to the pressure-enthalpy diagram, Figure 3-18, which shows the phase behavior of ethane. The Pressure-Enthalpy diagrams show pressure on the vertical axis and
enthalpy on the horizontal axis. The diagrams are used in locating pipeline operating
points in terms of pressure and temperature and for designing control valves. Pipe
flow is almost an isenthalpic process, so the diagram shows a graph of the enthalpy
during various pressures and physical states. The critical point is defined at the critical pressure and critical temperature (point C in the figure), where the liquid phase
and vapor phase meet, and either phase cannot be distinguished. The rectangular
box in the diagram shows the operating pressure range of an ethane pipeline for an
operating temperature range (assuming that the operating temperature ranges from
0C to 30C (solid lines in the figure) and the pressure from 4500 kPa to 10,000kPa).
Since the operating temperature range is lower than the critical temperature, the
ethane in this operating condition remains in liquid phase. For different operating
temperatures, the operating pressure range should be different to avoid vaporization.
Density: 865.4 kg/m3 at 15C and 875.4 kg/m3 at the operating temperature
Viscosities at 4C: 43.5 cSt
Pipe roughness: 0.0457 mm
Delivery pressure: 350 kPag
Load factor: 90%
Wall thickness
(in/mm)
Velocity (m/s)
Wall thickness
(in/mm)
Velocity (m/s)
18/457.2
20/508.0
22/558.8
0.250/6.35
0.250/6.35
0.250/6.35
2.24
1.80
1.48
0.281/7.14
0.281/7.14
0.281/7.14
2.26
1.81
1.49
X70
Pipe size
(in/mm)
Wall thickness
(in/mm)
Design pressure
(psig/kPag)
Wall thickness
(in/mm)
Design pressure
(psig/kPag)
20/508.0
20/508.0
22/558.8
22/558.8
0.250/6.35
0.281/7.14
0.250/6.35
0.281/7.14
1180/8132
1315/9067
1064/7334
1196/8246
0.250/6.35
0.281/7.14
0.250/6.35
0.281/7.14
1260/8688
1417/9765
1145/7897
1288/8880
Since the operating pressure range is between 8100 kPag and 9500 kPag,
the design pressures far outside of the range are removed from further consideration. Therefore, the selected combinations for X65 pipe are 20 with wall
thicknesses of 0.25 and 0.281, 22 with wall thickness of 0.281, and those
for X70 are 20 with wall thicknesses of 0.250 and 0.281, and 22 with wall
thickness of 0.281. The allowable design pressure for the 20 with 0.281 wall
thickness exceeds the maximum operating pressure, but the combination is
selected for further consideration because it is within a tolerance level.
Pipe grade
Pipe size
(in/mm)
Wall thickness
(in/mm)
Design
pressure
(kPag)
Total
pressure
drop (kPag)
Total
pressure
req. (kPag)
X65
X65
X65
X70
X70
X70
20/508.0
20/508.0
22/558.8
20/508.0
20/508.0
22/558.8
0.250/6.35
0.281/6.35
0.281/7.14
0.250/7.14
0.281/7.14
0.281/7.14
8132
9067
8246
8688
9765
8880
7288
7288
4649
7400
7400
4649
7638
7638
4999
7750
7750
4999
All six combinations satisfy the total pressure requirements for 20,000 m3/d
flow. For the same combinations as above, calculate the required total pressure
drop and inlet pressure for 22,300 m3/d from the fourth year to the tenth year.
Pipe grade
Pipe size
(in/mm)
Wall thickness
(in/mm)
Design
pressure
(kPag)
Total
pressure
drop (kPag)
Total
pressure
req. (kPag)
X65
X65
X65
X70
X70
X70
20/508.0
20/508.0
22/558.8
20/508.0
20/508.0
22/558.8
0.250/6.35
0.281/6.35
0.281/7.14
0.250/7.14
0.281/7.14
0.281/7.14
8132
9067
8246
8688
9765
8880
8811
8946
5618
8811
8946
5618
9161
9296
5968
9161
9296
5968
Only the combinations of pipe size 22 with the wall thickness of 0.281 for
X65, and of the pipe size 20 with the wall thickness of 0.281 and the pipe size
22 with the wall thickness of 0.281 for X70 pipe, satisfy the pressure requirement with no intermediate pump station. It may not be cost-effective to install
and operate an intermediate pump station to accommodate a small amount of
the flow increase from the fourth year.
For the above three combinations, calculate the required total pressure and
inlet pressure for the flow rate of 30,000 m3/d from the tenth year on. It should
be noted that the pumping power requirement for the 20 pipe at the inlet point
is higher by 59% (8946/5618 = 1.59) than the power requirement for the 22
pipe size. Therefore, the pump units for the 20 pipe have to produce higher
head than those for the 22 pipe and thus their capital and operating costs are
higher. On the other hand, the required pressure for the 22 pipe is low for the
first 10 years, and so the facility usage would be limited unless further flow
increase is expected in earlier years.
Pipe grade
Pipe size
(in/mm)
Wall thickness
(in/mm)
Design
pressure
(kPag)
Total
pressure
drop (kPag)
Total
pressure
req. (kPag)
X65
X70
X70
22/558.8
20/508.0
22/558.8
0.281/7.14
0.281/7.14
0.281/7.14
8246
9765
8880
9432
15,039
9432
9782
15,389
9782
Pipe grade
Pipe size
(in/mm)
Design pressure
(kPag)
Discharge
pressure at inlet
point (kPag)
Discharge pressure
at intermediate
station (kPag)
X65
X70
X70
22/558.8
20/508.0
22/558.8
8246
9765
8880
5066
7870
5066
5066
7870
5066
The capital cost due to the extra pumping power requirement for the 20
pipe is higher than the cost for the 22 pipe size, while the pipe cost for X70
with 20 diameter may cost less than the other two options. The extra capital
cost for the 22 line is more than 20% and is incurred by the extra pipe material
and construction expenses. However, the extra capital cost of the 22 diameter
pipe might be partly compensated by lower unit pumping cost. Assuming that
the annualized cost for the 20 pipe case is lowest, it is selected as the base
design.
The facilities such as the initiating pump station for the selected base design
would be designed to accommodate the capacity until the capacity increases
in the 10th year. In the 10th year, the additional facility increases include the
pumping capacity at the inlet point for the additional flow and an intermediate
pump station with the pumping capacity of 30,000 m3/d.
4. Develop alternative design cases and perform comparative studies against the
base design
Alternative 1: This alternative design is to use a pipe wall thickness larger
than 0.281 in order to increase the design pressure slightly higher than the
total pressure requirement. No intermediate pump station is required if the design pressure is slightly higher than the total pressure requirement. Note that
the required total pressure will be increased due to slightly smaller inside pipe
diameter.
The design pressure for the X70 22 pipe is lower than the total pressure requirement, which in turn is lower than the maximum operating pressure range.
The next largest nominal wall thickness is 0.312 or 7.92 mm, and its design
pressure is 9857 kPag or 1430 psig, but the required total pressure is 9914 kPag
for a flow rate of 30,000 m3/d. Therefore, the wall thickness is not sufficient to
meet the total pressure requirement without an intermediate pump station.
The next largest nominal wall thickness is 0.344 or 8.74 mm, which
can allow the design pressure to increase up to 10,870 kPag. For this wall
thickness, the required total pressure turns out to be 10,051 kPag. Since this
design pressure is higher than the required total pressure, no intermediate
pump station is required for the flow rate expected beyond the 10th year,
Design
Design
pressure
Base
Design 2
Design 3
9765
10,149
9765
Inlet pressure
(9th year flow)
Inlet pressure
(beyond 10th year)
9720
6234
9300
Discharge pressure
at intermediate
station
18,133
11,499
9300
9242
6275
9300
The base design barely satisfies up to the seventh year transportation requirement without an intermediate pump station, while Alternative 3 provides more than the ninth year requirement within the design pressure. Also,
Alternative 3 fully utilizes the facility, but Alternative 2 does not. Still, both
alternatives need an intermediate pump station from the 10th year on, and
the pumping capacity at the inlet station has to be increased at the same time.
The intermediate pump station will be located at 100 km from the inlet station, because the pumping head at both stations is the same (the criteria for
locating pump stations are discussed in the next chapter). With the discharge
pressure of 9300 kPag, Alternative 3 has a higher flow capacity than needed.
In summary, Alternative 3 is selected as the best design under this flow condition, because:
With slightly more capital and operating costs than the base design, Alternative 3 offers more flexible operation,
If needed, the flow capacity can be increased significantly.
Sensitivity on a fast flow growth: The flow rate is expected to grow at the
yearly rate of 1000 m3/d from the first year on and to level off at 32,000 m3/d;
18,000m3/d or design flow rate of 20,000 m3/d in the first year, 19,000 m3/d or
design rate of 21,100 m3/d in the second year, etc. For these flow rate changes,
the required total pressures are calculated for these three cases:
Design
Design pressure
Inlet pressure
(kPag)
Year of pump
installation
Base
Design 2
Design 3
9765
10,149
9765
9436
9782
9440
3rd year
10th year
4th year
Alternative 2 does not require an intermediate pump station until the 10th
year, while the other two require it in 3rd and 4th year, respectively. Alternative
2 needs higher initial capital cost due to the higher pipe grade and larger pipe
size. However, Alternative 2, using a large pipe size, offers a better option in
terms of the operating cost for such a high flow growth rate.
Assume that the design factor of 0.72 is applicable and that the elevation profile is
flat and flow is isothermal. Figure 3-10 shows the configuration of this pipeline system.
Solution:
It is assumed that the Alternative 3 design has been used for the main line in anticipation of flow increase and the intermediate pump station has been operating.
Step 1. Determine the design pressure of the main and branch lines using the Barlow formula with the hoop stress limited to 72% of the SMYS
Pmain = 2 70,000 0.281 0.72/20 = 1416 psig = 9765 kPag
Pbranch = 2 52,000 0.25 0.72/12.75 = 1286 psig = 8868 kPag
Step 2. Calculate the required pressures at TO and the discharge pressure at the
intermediate pump station.
1. First calculate the pressure required at TO for the design flow rate that can meet
the branch line delivery pressure requirement within the design pressure limit
of 8868 kPag.
Q = Q1 + Q2
(3 32)
where Q represents the flow rate in the base conditions. Applying the common pressure
principle, we have
PB PC = P1 = P2
Elevation (m)
Kilometer
post (km)
Elevation (m)
30
55
45
30
70
100
110
130
150
160
180
200
100
300
770
425
150
130
Assume that the design factor of 0.72 is applicable and that the flow is isothermal.
Solution:
It is assumed that the elevation changes are gradual between two profile points, the peak
point pressure is kept at 350 kPag, and the minimum suction pressures are the same as
the delivery pressure. An intermediate pump station is located at KMP = 110km. Note
that the elevation changes in the first section between CE and the intermediate station
are mild, but the changes in the second section are significant.
Step 1. Calculate the pressures for the design flow rate of 30,000 m3/d at the above
profile points. The discharge pressure at CE is 9220 kPag so as to satisfy the minimum
suction pressure requirement, and the discharge pressure at the intermediate station is
9091 kPag so as to keep the peak point pressure at 350 kPag. The delivery pressure of
2067 kPag is obtained in order to keep the pipeline flow in a full flow condition. As a
result, the pressure profile is determined as shown in the table below.
KMP (km)
0
20
30
60
80
90
Elevation (m)
Pressure (kPag)
KMP (km)
Elevation (m)
Pressure (kPag)
30
55
45
30
70
100
9220
7502
6836
4708
2862
1854
110
130
150
160
180
200
100
300
770
425
150
130
350/9091
5876
350
2550
3400
2067
Elevation (m)
Pressure (kPag)
KMP (km)
Elevation (m)
Pressure (kPag)
30
55
45
30
70
100
1975
1575
1567
1415
886
536
110
130
150
160
180
200
100
300
770
425
150
130
350/6456
4558
350
3209
5375
5360
Product
Density at 4C
(kg/m3)
Viscosityat 4C
(cSt)
Batch size at 4C
(m3)
Vapor pressure
(kPa)
32API
35API
Condensate
875.4
857.6
705.0
43.5
21.0
0.7
20,000
15,000
25,000
10
15
95
It is assumed that these values are measured at the average operating temperature
of 4C. Design the batch pipeline including the delivery pressure.
Solution:
It is assumed that the base design is used; an intermediate pump station is located
100km downstream of CE and the design pressure is 9765 kPag.
Step 1. Determine the line fill volumes of the two sections of the pipeline. The
line fill volume is the volume of liquid contained in a segment of pipe, and is the pipe
volume in the ambient conditions, even though actual volume of liquid shrinks under
pressure. Addendum 3.4 discusses the effect of pressure and temperature on line fill
volume. A section is defined as the pipeline between two pump stations or between a
pump station and the delivery point. Therefore, the first section is defined from CE to
the intermediate station, where the second section starts, ending at QU.
Since the length of each section is the same, so is the line fill volume of each
section. Assuming that the pipe volume does not change in a pressurized condition, the line fill volume of each section becomes 19,150 m3.
Since this volume is smaller than the size of a 32API batch, the batch covers
thewhole section when fully lifted at CE or has passed the intermediate station.
Step 2. Select the product with the largest viscosity and the product with the highest vapor pressure among the three batch products.
The 32API batch has the highest viscosity among the three products, with a
viscosity of 43.5 cSt.
The condensate batch has the highest vapor pressure of 95 kPa or 6 kPag. Taking into account the minor pressure losses and transient effect, extra pressure of
400 kPa is added to the highest vapor pressure to get the minimum pressure of
400 kPag at the delivery and pump station.
Step 3. Determine the batch sequence. Here, the batch sequence is given below
without describing the sequencing method which is detailed in Chapter 5.
The batch sequence for minimizing the interfacial mixing is 32API 35API
Condensate, and the same sequence is repeated in the next batch cycle.
When these batches are placed in the pipeline, the batch line fill profile can be
shown in Figure 3-17.
Step 4. Calculate the pressure profile using the 32API properties.
Discharge pressure at CE = 7920 and suction pressure at the intermediate station = 400 kPag.
Discharge pressure at the intermediate station = 7920 and delivery pressure at
QU = 400 kPag.
Step 5. Determine the average pressure profile.
Calculate the pressure drops P32, P35 and Pcon for 32API, 35API and condensate, respectively, using the common minimum pressure of 400 kPag.
Calculate the total required pressure averaged over the weight of each batch
size:
Pavg = (20,000 P32 + 15,000 P35 + 25,000 Pcon)/(20,000 + 15,000 + 25,000)
This approach is acceptable if the load factor is low.
Vapor pressure
(kPa)
Thermal
expansion (/C)
Critical
pressure (kPa)
Critical
temperature (K)
I-Pentane
N-Pentane
I-Butane
N-Butane
Propane
Ethane
151
116
530
379
1370
6000 (*)
0.00160
0.00154
0.00216
0.00194
0.00280
0.015 (+)
3381
3370
3640
3798
4244
4872
460.4
469.7
407.8
425.1
369.8
305.3
Ethylene
9700 (*)
0.025 (+)
5040
282.3
(*) The vapor pressures and thermal expansions of these liquids are highly dependent on the pressure and temperature conditions. Therefore, a representative value does not have a definite meaning
for these products. These values are estimated by extrapolating measured values and are presented
for an illustrative purpose only.
(+) These values are estimated about 40C at 9000 kPa and presented for an illustrative purpose
only.
As shown in the table, the vapor pressures and thermal expansions of ethane and
ethylene are significantly higher than the other HVP liquids. Normally, these two products are transported in dense phase. For a hydrocarbon mixture, there is no clear line
dividing dense phase from the liquid phase or other single line dividing the dense
phase from the gas phase, but the dense phase lies between critical temperature and
cricondentherm if the pressure is above the cricondenbar. Phase change from denseto-liquid or vice versa is gradual. Ethane (C2H6), ethylene (C2H4), and carbon dioxide
(CO2), can be liquefied in pipelines at temperature and pressures even below the
critical point, and treated as liquids in transportation. Dense phase liquid is a highly
compressible liquid that shows properties of both liquid and gas; a density similar to
that of a liquid, but a viscosity similar to that of a gas. For liquid pipeline design and
operation, it is considered that the fluids are in dense phase if the pressure and temperature are around the critical pressure and critical temperature but above the vapor
pressure.
Because of the complex dependence of fluid properties on pressure and temperature in or near to the dense phase, pressure, and temperature should be determined
as accurate as possible and thus their calculations must be performed simultaneously
to achieve the desired accuracy. Reference [12] details the method of calculating
pressure and temperature in dense phase and identifies the following key design
parameters:
The critical point is not well defined nor are the properties near the critical
point. Therefore, one should try to avoid approaching the critical points too
closely.
Since most ethane or ethylene pipelines are operated in a fully turbulent flow
regime, the friction factor is independent of the Reynolds number and depends
only on the relative roughness of the pipe. Therefore, the accuracy of the pressure profile is sensitive to the values of the relative roughness. Sometimes,
other HVP products flow in a similar fully turbulent regime.
Both pressure and temperature profiles are relatively sensitive to the specified
value of the overall heat transfer coefficient, which in turn depends on soil
conductivity. The soil conductivity not only varies along the pipeline but also
changes frequently with moisture content.
The effect of the seasonal variation in the average soil temperature depends
on the difference between the fluid temperature and soil temperature. If the
Solution:
Refer to the pressure-enthalpy diagram, Figure 3-18, which shows the phase behavior of ethane. The Pressure-Enthalpy diagrams show pressure on the vertical axis and
enthalpy on the horizontal axis. The diagrams are used in locating pipeline operating
points in terms of pressure and temperature and for designing control valves. Pipe
flow is almost an isenthalpic process, so the diagram shows a graph of the enthalpy
during various pressures and physical states. The critical point is defined at the critical pressure and critical temperature (point C in the figure), where the liquid phase
and vapor phase meet, and either phase cannot be distinguished. The rectangular
box in the diagram shows the operating pressure range of an ethane pipeline for an
operating temperature range (assuming that the operating temperature ranges from
0C to 30C (solid lines in the figure) and the pressure from 4500 kPa to 10,000kPa).
Since the operating temperature range is lower than the critical temperature, the
ethane in this operating condition remains in liquid phase. For different operating
temperatures, the operating pressure range should be different to avoid vaporization.
Heavy crude is characterized by high density, high viscosity and high pour point,
and may contain a significant amount of wax and/or sulphur. Heavy crude may exhibit
non-Newtonian viscosity behavior at normal operating temperature ranges because its
pour point can be higher. It is known that the apparent viscosities of non-Newtonian
liquids are sensitive not only to temperature changes but also to the shear rate and cooling rates. Laboratory tests should be performed at the pipeline operating conditions to
determine the crudes viscosity types and behaviors in terms of the shear stress vs. shear
rate and yield stress over the operating temperature ranges including the pour points. The
types include Newtonian, dilatant, Bingham plastic, pseudoplastic, and thixotropic (timedependent) fluid, because heavy crudes show different fluid characteristics.
Another potential engineering problem in dealing with heavy crudes, and sometimes with light and intermediate crudes, is the significant presence of wax. A waxy
crude may exhibit Bingham plastic characteristics after gelling, requiring a finite shear
stress to initiate flow. Heavy and/or waxy crudes start developing a yield stress near
their pour point, which may require additional pressure to restart flow. It is known that
wax does not deposit in turbulent flow at high temperatures, certain parts of a pipeline
may have wax deposits, and wax deposits could have an insulating effect.
Table 3-7. Viscosity, temperature, and pour point [14]
Product
Specific
Gravity
Temperature
(C)
Viscosity
(cSt)
Temperature
(C)
Bitumen
Residuals
Crude
High wax
Diesel
Jet fuel
Gasoline
NGL
1.02
0.96
0.84
0.81
0.84
0.78
0.73
0.50
65
65
20
50
-1
-1
-1
-1
50,000
1000
11
7.4
2.8
2.2
0.8
0.23
120
120
50
60
27
27
27
38
0.2
55
32
13
35
(3 33)
where
T(t) = temperature at time t, (C)
T(0) = temperature at the time when the pipeline is shut down (C)
Tg = ground temperature (C)
DT = pipe outside diameter including insulation thickness for insulated pipe
t = time from start of static cooling, in second
d = inside pipe diameter, m
r = liquid density, kg/m3
Cp = heat capacity of the liquid averaged at the inlet and discharge temperatures,
kJ/kgC
If the heat conduction through the boundary layer and pipe wall is excluded
from the overall heat transfer coefficient, the calculated temperature would be lower
than the actual temperature, requiring somewhat lower restart pressure. If possible,
pipeline start-up or restart can be scheduled during periods of warmest ambient temperatures in order to avoid the difficult problems that may be encountered during
start-up.
3.3.7.4 Design Facilities
The effect of the yield stress of heavy crudes is non-trivial below pour point. When a
pipeline is shut in and thus the heavy crude cools down below the pour point, it requires
an extra pressure to put the crude in motion. This extra pressure requirement has to be
provided by a pump to initiate flow. The pressure required to initiate flow is sum of the
pressure differentials required to break the gel in each section of the pipeline. Since
yield strength is sensitive to temperature, the required pressure has to be determined on
each segment to reduce potential calculation error.
When starting up the pipeline after shut-in, the flow rate should be very low
to push the gelled crude gently. It is essential to establish the minimum flow rate
needed to be maintained during initial start-up, and it may be necessary to include
redundant provisions for emergency and planned shutdowns. In selecting a mainline pump, the maximum operating point should be satisfied as usual. If the minimum flow cannot be met by the mainline pump during an initiating period after
shut-in, special startup/restart pumps with the capability of high pressure and low
flow should be considered. Note that the performance of a centrifugal pump deteriorates for pumping high viscosity fluids, thus requiring a rerate of the pump perfor
mance. Refer to the next chapter for pump performance rerating for high viscosity
conditions.
Systems to consider would include standby pumps for displacing the crude oil
in the pipeline with water, and adding pour point depressant injection facilities. If
bitumen is blended with a diluent, a diluent blending and storage facility is required
at the lifting point. If the crude is heated for pumping, heaters have to be installed not
only at the initiating station but also intermediate pump stations. Depending on the
(3 34)
where
qh = heater duty required to heat the liquid to the discharge temperature, kJ/hr
Solution:
It is assumed that the pour point is low enough to transport the blended bitumen as a
Newtonian fluid and that the possible contents of diluent for winter condition are 45%,
40% and 35% and the contents for summer condition are 30%, 25% and 20%.
Density
(kg/m3)
Thermal
Expansion (C)
Minimum
Temperature (C)
862.4
8.27 104
321
34,800
40
878.8
7.95 10
15
330
33,600
35
895.2
7.66 104
23
342
32,400
30
911.6
7.37 104
32
330
31,200
25
20
928.0
944.4
7.12 104
41
50
324
339
30,000
28,800
45
6.89 104
Step 2. Calculate the pressure and temperature profiles. The maximum design
pressure for the selected pipe is 6950 kPag, so the blended bitumen is discharged at
6900 kPag.
The table below summarizes the injection and delivery temperatures and the delivery pressures for different amounts of diluent. The injection delivery temperatures are
determined in such a way that the viscosity at the delivery point is kept below 350 cSt.
The table shows the temperatures and their corresponding viscosities after the pipeline
is shut down for 120 hours.
Contents of
Diluent (%)
45
40
35
30
25
20
Injection
Temperature
(C)
Delivery
Temperature
(C)
Delivery
Pressure
(kPag)
Temperature
after 120 hours
(C)
Viscosity after
120 hours
(cSt)
15
34
52
62
83
104
7
15
23
32
41
50
532
1458
2019
2236
2774
3203
0.2
3.3
6.4
14.7
18.2
21.8
750
1264
2138
1917
2958
4359
The required injection temperature decreases as the amount of the diluent increases. As expected, the pressure requirements for the winter condition are higher
than those for the summer condition as a result of the operating temperatures in summer condition being much higher than in winter and the flow rates are lower for the
summer condition. Also, as the amount of diluent gets smaller, a heating facility has to
be installed to raise the injection temperature.
If this bitumen starts showing its non-Newtonian behavior about 2000 cSt, the
yield stress has to be measured in order to assess the requirement for extra pumping
facilities to dislodge the blended bitumen that was congealed during the 120 hours of
the shut-in period.
Step 3. Repeat the same calculations for the case where the pipeline is insulated
with 2 polyurethane insulation material. The insulation conductivity is 0.035W/mC.
Injection
Temperature
(C)
Delivery
Temperature
(C)
Delivery
Pressure
(kPag)
Temperature
after 120 hours
(C)
Viscosity after
120 hours
(cSt)
10
18
29
38
51
63
9.7
16
23
31
40
48
390
381
841
1139
1706
2109
5.6
9.8
15.4
23.8
30.4
36.5
410
610
774
742
853
997
Compared to the results of the un-insulated case, the insulated pipeline has the
following advantages over the un-insulated pipeline:
No extra pumping facilities are required even after shutting down for 120
hours,
There is less need for a heater because the required injection temperature is low,
or the heater duty is lower than the duty for the un-insulated pipeline even if a
heater is installed, assuming that the same amount of the diluent is mixed,
The diluent requirement is much smaller than the requirement for the uninsulated pipeline,
Restarting after the shut-down is much easier due to low viscosity.
Step 4. Finalize the pipeline system design
The 28 pipe with wall thickness of 0.35 and pipe grade X56 satisfies the pressure requirements for both winter and summer conditions.
The insulated pipeline can save both capital and operating costs by reducing
the requirements for extra facilities such as a heater and an extra pump to deal
with the congealed non-Newtonian crude.
The selection of the diluent requirement vs. heater installation is based on the
cost comparison of the diluent costs against the heater costs. If the pipeline is
insulated, 35% of diluent and 65% of bitumen blending can be sufficient in
winter. If the temperature of the blended bitumen is higher than 63C, a heater
is not required and thus 20% of diluent may satisfy the summer transportation
requirement.
For a simple pipeline system with relatively flat terrain, the criteria for locating
stations results in almost equal station spacing along the pipeline, and the number of
pump stations can be determined by dividing the total required pressure by the difference between the MAOP and the minimum pressure;
No. of stations = Total required pressure / (MAOP minimum suction pressure)
The procedure of locating pump stations is to start from the delivery pressure,
drawing the pressure gradient upstream to the intersection of the maximum design
pressure, which is superimposed on the elevation profile. If the discharge pressure of
the initiating station is smaller than the design pressure, then reduce the design pressure
and move the initial locations to further downstream locations. The same differential
pressure can be calculated by dividing the total pressure requirements by the number
of pump stations.
Example 1: Simple Pump Station Location
Refer to the design example described in Section 3.3.1. The total required pressure
is 15,389 kPag, maximum design pressure 9765 kPag, and minimum delivery pressure 350 kPag. It is assumed that the minimum suction pressure is 350 kPag. Since
the elevation profile is flat, the number of pump stations is obtained from the above
formula:
Therefore, the total number of pump stations required is 2; one at the initiating station and the other at an intermediate location. Applying the station location criteria, the
intermediate station is located at the mid-point of the pipeline as shown in Figure 3-20.
For a simple pipeline system with severe elevation changes, the station locations
can be determined by applying these criteria through trial and error on a graph. The
procedure of locating intermediate pump stations is as follows:
Step 1. Using the maximum design pressure as the discharge pressure at the initiating station, the first intermediate station is found at a location where the pressure
reaches the minimum suction pressure by drawing the pressure gradient on the elevation profile. In practice, a pressure allowance of 200 kPa to 300 kPa at the intake of the
pump station is required to account for the losses due to station piping, valves, fittings,
and other equipments.
PD = 7,780
0
Booster Pump
PS = 350
100 km
Main Line Pump
Distance
200 km
Step 2. Progressing downstream from the maximum design pressure at the intermediate station, the next intermediate station is located in the same way as above.
Repeat these steps until the minimum suction pressure of the last section is greater than
or equal to the delivery pressure.
Step 3. If the suction pressure is much greater than the delivery pressure, reduce
the discharge pressure equally at each pump station and then repeat the second step to
move the initial locations to upstream locations.
Step 4. If the discharge pressure has to be reduced significantly, the maximum design
pressure can be lowered by selecting lower grade pipe or thinner pipe wall thickness.
Figure 3-21 shows that the total pressure requirement is greater than the design
pressure. This pressure requirement can be met by installing an intermediate pump station or choosing a thicker pipe in the upstream segment where the required pressure is
Pressure
( kPag)
Head
(m)
Design Pressure
8,600
1,000
PD
PD
4,300
500
Ps
PB
PS
0
Distance (km)
200
The pipeline length is 350 km and the elevation profile is given below.
KMP
Elevation
0
10 m
50
250 m
150
250 m
190
250 m
230
250 m
290
310 m
320
460 m
350
10 m
Solution:
Step 1. Determine the design flow rate and the maximum design pressure.
Since the number of yearly operating days is 345 days, the load factor is
345/365 = 94.5%, and thus the design flow rate is 1175/0.945 = 1243 m3/hr, or
rounding up to 1250 m3/hr.
The design pressure is obtained by applying the Barlow formula and the design
factor for the X60 pipe grade; 8370 kPag.
Step 2. Calculate the pressure gradient.
Step 3. Determine initial station locations and calculate the pressures at the
locations.
Assuming that the discharge pressure at the initiating station is 8230 kPag, the
first pump station will be located approximately 53 km downstream with a suction pressure of 354 kPag.
Since the elevation difference is zero for about 180 km downstream of the
first intermediate station, the next pump station can be located with a similar
Pressure
(kPag)
8,600
8,116
8.105
PD
8,105
8,160
1,000
8,079
4,300
300
350
350
500
350
350
887kPag
0
50
100
150
200
250
300
350
PS
Distance (km)
Assume that the design factor of 0.72 is applicable and that the elevation profile is
flat and flow is isothermal.
Solution:
Step 1. Calculate the design pressure (MAOP) of the main and branch lines.
MAOP of the main line = (2 60,000 0.281 0.72/20) 6.895 = 8470 kPag
MAOP of the branch line = (2 52,000 0.250 0.72/16) 6.895 = 8067 kPag
Step 2. Calculate the pressure required at TO on the branch line side.
ADDENDA TO CHAPTER 3
A3.1 Temperature Calculation
Temperature has considerable influence on the design of pipelines and related facilities, including the establishment of facilities sizing and optimization, economic and
technical evaluation, etc.
Temperature and pressure influence all fluid properties. In fluid transmission pipelines, both pressure and temperature vary along the pipeline length. In long-distancetransmission pipelines traversing varied terrain, from permafrost regions to moderate
climate conditions, pipelines experience significant temperature changes. Temperature
change affects viscosity, density, and specific heat in liquid lines, particularly in crude
oil pipelines.
In any pipeline segment, the significant overall temperature change (DT) is due to
conduction and convection (DTc). However, there are other factors that affect the overall temperature change. These are (DTe) due to isentropic expansion caused by elevation change and due to isenthalpic expansion caused by friction (DTf) [17]. Therefore,
the overall temperature change in a pipeline segment is:
(A3 1)
The following illustrates a method of overall temperature change due to conduction and convection, DTc. For a pipeline (Figure A3-1) buried at a finite depth (ho)
with insulation, the following expression for computing fluid flow temperature To is
applicable, Holman [18].
Nomenclature:
Cp
= Fluid isobaric specific heat
Dp or D = Pipe outside diameter
2p DLki To - Tg
ln Di / Dp
(A32)
By introducing the shape factor, S, and rearranging the above equation, we have
To =
kg S
Tg
1+ a
kg
rQCp +
1+ a
rQCpTi +
(A33)
where
kg Di
ln
ki Dp
(A34)
2p DL
2
2h
2h
ln o + o - 1
Dp
Dp
(A35)
a=
and the shape factor is defined as
S=
rQCpU
(A36)
1 1
1
D 1
D
+
ln
+
ln i
p hf d 2kp d 2ki d
+ h D
a i
(A37)
In Eq. (A37), radiation heat losses are ignored as they are small at most normal pipeline operating temperatures. When the pipe is not insulated, the third term in
Eq.(A37) is reduced to zero and Di in the fourth term is set equal to D (i.e., outside
diameter of the pipe).
For above-ground pipeline, the film coefficient (ha) for air can be calculated from
the following equation recommended by Dittus and Boelter, Holman [18].
N u = 0.023 ( Re )
0.8
( Pr )n
(A38)
Where the Nusselt number, Nu, Reynolds number, Re, and Prantl number, Pr, are defined as follows:
h D
Nu = a
ka
rV D
Re = a a
ma
-m aCpa
Pr =
ka
0.4
(Vf )0.8
d 0.2
(A39)
k r V D
ha = C a a a i
Di m a
(A310)
Properties of air are provided elsewhere [18]. Values of constant C and exponent n
are dependent on the Reynolds number and are also given elsewhere [18].
For an offshore pipeline, ha can be calculated from Eq. (A3-10) with appropriate
values of Cp, r, k, and m for sea water, and knowing the current velocity.
The following expressions summarize the computation of DTe and DTf.
T V
T
DTe = - Es Ph =
=
Ph
P s Cp T p
(A311)
And
DTf = - JPf =
H
P
T
H
T Pf
P
=-
1 H
Pf
Cp P T
(A312)
Where Es and J are elevation sensitivity and Joule Thompson coefficients, respectively, Ph is pressure loss due to elevation change, and Pf is pressure loss in overcoming
friction. Es can be computed from graphs of pressure (P) and temperature (T) at constant entropy (s), and Pk can be calculated from graphs of enthalpy (H) versus pressure
(P) at constant temperature (T).
The sign of Joule Thompson coefficient J indicates whether fluid expansion
or compression will cause an increase or decrease in the temperature. As an example, in an expanding gas if J is positive, the gas will cool. A negative J in an
expanding gas indicates temperature rise, and is observed in expansion of some
specialgases,e.g., hydrogen. Methods for calculating Es and J are given elsewhere
[19].
The above procedure outlined above provides an accurate prediction of fluid flow
temperature under steady-state condition for buried and exposed pipelines. Sample
plots of temperature profiles for a liquid pipeline (carrying bitumen/condensate mixture) is shown in Figure A3-3.
Figure A3-3. T
emperature profile in a buried 12 pipeline transporting non-Newtonian
fluid[17]
rm =
(A313)
where
Ve = Erosional velocity, ft/sec
C = Constant defined as 75< C <150
P = Operating pressure, psia
Si = Liquid specific gravity (water = 1; use average specific gravity for hydrocarbon mixes such as interface) at standard conditions.
Ra = gas/liquid ratio, ft3/bbl at standard conditions
T = Operating Temperature oR
Sg = Gas specific gravity (air = 1) at standard conditions
Z = compressibility factor at the specified pressure and temperature,
Depending on gas/liquid ratio and the amount of entrainment of particulates, it is
a standard practice to limit the fluid velocity about 5 m/s. However, in two-phase flow
lines minimization of slugging in pipeline and separation equipment set the maximum
velocity. In this case the velocity is set to about 10 ft/sec (3 m/sec), ANSI/API RP-14E
[20]. This is particularly important in long lines with elevation changes.
(A314)
(A315)
where
DPm = Minor pressure loss
K = pressure loss coefficient (dimensionless).
V = velocity of liquid through valve, joint or fitting
r = liquid density
f = friction factor
D = pipe inside diameter
Le = equivalent length of straight pipe
The value of K is determined mainly by the flow geometry or the shape of the device. K is analogous to the term fL/D for a straight length in the Darcy equation.
It is more convenient to use K, because K is less dependent on the Reynolds
number and relative roughness. K and Le values are obtained from Reference [21].
Some representative values of K and Le are shown in Table A3-1. Since actual values
depend on the design and fabrication of the components, actual manufacturers data
shall be consulted in the final design.
Example of Equivalent Length Method: A piping system of a pump station is
140 m of NPS 24 pipe that has two 24 gate valves, three 24 ball valves, on swing
check valve, and two 90 standard elbows. Using the equivalent length concept, calculate the total pipe length of the station.
Solution:
Convert all valves and fittings in terms of 24 pipe as follows:
Types
Pipe entrance:
Sharp edged
Well-rounded
Pipe exit:
Sharp edged
Well-rounded
90 Bends:
Smooth threaded
Smooth flanged
Miter bend
Miter bend with vanes
Tee Branch flow:
Flanged
Threaded
Tee line flow:
Flanged
Threaded
0.5
0.03
1.0
1.0
0.9
0.3
1.1
0.2
1.0
2.0
Types
Valves (fully open):
Gate valve
Globe valve
Plug valve
Angle valve
Ball valve
Globe lift check valve
Angle lift check valve
Standard elbow:
45
90
Standard tee:
Flow through run
Flow through branch
Equivalent length
(Le/D)
(*)
8
340
18
150
3
600
55
16
30
20
60
0.2
0.9
(*) K values for valves vary with their type and size.
Total equivalent length of straight pipe and all fittings = 140 m + 112.78 m =
252.78 m
The pressure drop due to friction is calculated based on 252.78 m of pipe.
Gradual Enlargement and Reduction
Consider liquid flowing through a pipe of diameter D1. If the diameter enlarges or
reduces to D2, the pressure loss can be calculated as follows:
(A316)
where v1 and v2 are the velocity of the liquid in D1 and D2 pipes and A1 and A2 the areas.
The value of K depends on the diameter ratio and the different angle due to the enlargement or reduction. The figure below shows a gradual pipe enlarger and a reducer.
20
45
60
90
180
A2/A1 = 2.25
0.45
1.06
A2/A1 = 9
0.43
0.88
1.19
1.11
1.00
1.02
1.04
1.02
The loss coefficients for a reducer are given in the table below:
Table A3-3. Loss coefficients for reducer
Angle
20
45
60
90
180
A2/A1 = 0.5
0.05
0.06
0.06
0.12
0.26
A2/A1 = 0.25
0.03
0.07
0.07
0.17
0.41
Sudden Expansion: Minor loss for the sudden expansion of a pipe can be expressed as
DPm = r(1 D12/D22) 2 V12/2
(A317)
Sudden Contraction: Minor loss for the sudden contraction of a pipe can be
expressed as
DPm = r(1/Cc 1) 2 V 22/2
(A318)
0.0
0.585
0.2
0.632
0.4
0.659
0.6
0.712
0.8
0.813
1.0
1.0
Example: A tank open to the atmosphere is filled with 40API (specific gravity
of 0.825) oil to a height of 10 m from the bottom. A tap at the bottom of the tank is
opened, and oil flows out from the outlet. Assuming that the flow is steady and incompressible, determine the tank discharge pressure at point 2.
Solution:
Assume that as the oil is discharged out of the tank, the tank level drops so slowly that
the flow velocity at point 1 is negligibly small. The data for this example are:
Pipe = 16 OD 0.25 wt: inside diameter = 0.394 m
Pipe roughness = 0.0018 = 0.0457 mm
Crude viscosity = 3 cSt of 40API crude
10m
Oil
200m
16 OD
where f is the Darcy friction factor and L the pipe length, and K = Kentrance + Kexit are
the pressure loss coefficients for pipe entrance and exit, respectively.
Step 2. P1 is atmospheric pressures, v1 = 0, v2 = v, z1 = 10 m, and z2 = 0. Therefore,
the equation becomes
P2 = rgz1 rv2/2 frv2/2 L (Kentrance + Kexit)rv2/2
Static pressure at the tank outlet, rgz1 = 825 9.8 10/1000 = 80.9 kPa
Velocity pressure, rv2/2 = 825 (2.28) 2/2 = 2144 pa = 2.14 kPa
Pressure losses due to entrance and exit, (Kentrance + Kexit)rv2/2 = 1.5 2.14 =
3.21
Pressure at the pipe discharge point = 80.9 2.14 - 17.0 3.21 = 58.6 kPa
Step 5. Discuss the result of this example
The pressure losses due to pipe entrance and exit can be relatively significant for
low pressure system, whereas they are negligibly small for large pressure system such
as a transmission pipeline.
dVpipe = dVfluid
(A319)
By rearranging the expressions for both the bulk modulus and the coefficient of
volume thermal expansion of the fluid, it can be written that
dVfluid = BV dT - V dP / K
(A320)
where
B = coefficient of volumetric thermal expansion of the fluid.
K = bulk modulus of fluid
The determination of the incremental increase in volume of the pipe (dVpipe) is
dependent upon the degree to which the pipeline is free to move. Consider first the
case where a welded pipeline has unrestrained movement such that it is free to expand
in both the radial and longitudinal directions.
S - nSx DdP
eq = q
= 4tE (2 - n)
(A321)
S - nSq DdP
ex = x
= 4tE (1 - 2n)
(A322)
and
where
n = Poissons ratio
eq = Hoop strain (i.e., strain in radial direction)
ex = Axial strain (i.e., strain along the length of the pipeline)
L* = L(1 + ex)
(A323)
D* = D(1 + eq)
(A324)
The new volume of the pipe section is given by pL*(D*)2/4. Upon substitution
of the expressions for D* and L*, and ignoring the cross products of strain as being
negligible, the new pipe volume V* is given by:
V* =
pD 2 L
(1 + 2e q + e x)
4
(A325)
dVpipe =
VDdP
(5 - 4 n)
4tE
(A326)
dD = DadT
(A327)
dL = LadT
(A328)
and
dV = 3aV dT
(A329)
dP
DdP
V
(5 - 4 n) + 3adT = V BdT 4tE
(A330)
( B - 3a ) dT
1 D (5 - 4 n)
+
4tE
(A31)
B - 2a (1 + v ) dT
dP =
D 1 - v 2 / tE + 1/ K
(A332)
The effect of axial thermal stress on restrained pipelines capped at the ends can be
neglected [22]. Ignoring the temperature increase effects simplifies Equation (A3.4.14)
such that, for a restrained pipe, the theoretical volume required to pressure (squeeze)
the pipeline section is given by:
1
D
dV = V dP
1 - v2 +
tE
K
(A333)
Example:
To illustrate the effect of a change of temperature on pressure consider a pipe of outside
diameter 114.3 mm and wall thickness 3.18 mm transporting water under pressure.
Neglect the effect of axial thermal stress on a restrained pipeline Eq. (A332) can be
simplified to become [22]:
dP =
B - 2a
D
1 - v2 + C
Et
Figure A3-8. P
ressure change due to temperature change (medium water in 48 pipeline, no
flow condition)
E = 2.0 108
n = 0.3
C = 1/K = 4.86 107 (K= 2.057 106 kPa assumed for water at 10C)
Note that C (the reciprocal of the bulk modulus for liquid, in this case water) is
a function of the pressure and temperature. Assuming the temperature T is 10C then
upon substituting, B has the value 81.0728 106 and the pressure increase due to a
1C temperature change is:
dp =
REFERENCES
[1] American Petroleum Institute, 1984, Petroleum Liquid Volume Correction, API Publication
1101.
[2] NIST Standard Reference Database 4, Supertrap version 3, National Institute of Standards and
Technology, Gaithersburg, MD 20899, U.S.A.