You must log in to edit PetroWiki. Help with editing
Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information
Pressure drop evaluation along pipelines
The simplest way to convey a fluid, in a contained system from Point A to Point B, is by means of a conduit or pipe (Fig. 1).
Contents
Piping design
The minimum basic parameters that are required to design the piping system include, but are not limited to, the following.
- The characteristics and physical properties of the fluid.
- The desired mass-flow rate (or volume) of the fluid to be transported.
- The pressure, temperature, and elevation at Point A.
- The pressure, temperature, and elevation at Point B.
- The distance between Point A and Point B (or length the fluid must travel) and equivalent length (pressure losses) introduced by valves and fittings.
These basic parameters are needed to design a piping system. Assuming steady-state flow, there are a number of equations, which are based upon the general energy equation, that can be employed to design the piping system. The variables associated with the fluid (i.e., liquid, gas, or multiphase) affect the flow. This leads to the derivation and development of equations that are applicable to a particular fluid. Although piping systems and pipeline design can get complex, the vast majority of the design problems encountered by the engineer can be solved by the standard flow equations.
Bernoulli equation
The basic equation developed to represent steady-state fluid flow is the Bernoulli equation which assumes that total mechanical energy is conserved for steady, incompressible, inviscid, isothermal flow with no heat transfer or work done. These restrictive conditions can actually be representative of many physical systems.
The equation is stated as
(Eq. 1)
where
Z | = | elevation head, ft, |
P | = | pressure, psi, |
ρ | = | density, lbm/ft^{3}, |
V | = | velocity, ft/sec, |
g | = | gravitational constant, ft/sec^{2}, |
and | ||
H_{L} | = | head loss, ft. |
Fig. 2 presents a simplified graphic illustration of the Bernoulli equation.
Darcy’s equation further expresses head loss as
(Eq. 2)
and
(Eq. 3)
where
H_{L} | = | head loss, ft, |
f | = | Moody friction factor, dimensionless, |
L | = | pipe length, ft, |
D | = | pipe diameter, ft, |
V | = | velocity, ft/sec, |
g | = | gravitational constant ft/sec^{2}, |
ΔP | = | pressure drop, psi, |
ρ | = | density, lbm/ft^{3}, |
and | ||
d | = | pipe inside diameter, in. |
Reynolds number and Moody friction factor
The Reynolds number is a dimensionless parameter that is useful in characterizing the degree of turbulence in the flow regime and is needed to determine the Moody friction factor. It is expressed as
(Eq. 4)
where
ρ | = | density, lbm/ft^{3}, |
D | = | pipe internal diameter, ft, |
V | = | flow velocity, ft/sec, |
and | ||
μ | = | viscosity, lbm/ft-sec. |
The Reynolds number for liquids can be expressed as
(Eq. 5)
where
μ | = | viscosity, cp, |
d | = | pipe inside diameter, in., |
SG | = | specific gravity of liquid relative to water (water = 1), |
Q_{l} | = | liquid-flow rate, B/D, |
and | ||
V | = | velocity, ft/sec. |
The Reynolds number for gases can be expressed as
(Eq. 6)
where
μ | = | viscosity, cp, |
d | = | pipe inside diameter, in., |
S | = | specific gravity of gas at standard conditions relative to air (molecular weight divided by 29), |
and | ||
Q_{g} | = | gas-flow rate, MMscf/D. |
The Moody friction factor, f, expressed in the previous equations, is a function of the Reynolds number and the roughness of the internal surface of the pipe and is given by Fig. 3. The Moody friction factor is impacted by the characteristic of the flow in the pipe. For laminar flow, where Re is < 2,000, there is little mixing of the flowing fluid, and the flow velocity is parabolic; the Moody friction factor is expressed as f = 64/Re. For turbulent flow, where Re > 4,000, there is complete mixing of the flow, and the flow velocity has a uniform profile; f depends on Re and the relative roughness (Є/D). The relative roughness is the ratio of absolute roughness, Є, a measure of surface imperfections to the pipe internal diameter, D. Table 9.1 lists the absolute roughness for several types of pipe materials.
If the viscosity of the liquid is unknown, Fig. 4 can be used for the viscosity of crude oil, Fig. 5 for effective viscosity of crude-oil/water mixtures, and Fig. 6 for the viscosity of natural gas. In using some of these figures, the relationship between viscosity in centistokes and viscosity in centipoise must be used
(Eq. 7)
where
γ | = | kinematic viscosity, centistokes, |
ϕ | = | absolute viscosity, cp, |
and | ||
SG | = | specific gravity. |
Pressure drop for liquid flow
General equation
Eq. 3 can be expressed in terms of pipe inside diameter (ID) as stated next.
(Eq. 8)
where
d | = | pipe inside diameter, in., |
f | = | Moody friction factor, dimensionless, |
L | = | length of pipe, ft, |
Q_{l} | = | liquid flow rate, B/D, |
SG | = | specific gravity of liquid relative to water, |
and | ||
ΔP | = | pressure drop, psi (total pressure drop). |
Hazen Williams equation
The Hazen-Williams equation, which is applicable only for water in turbulent flow at 60°F, expresses head loss as
(Eq. 9)
where
H_{L} | = | head loss because of friction, ft, |
L | = | pipe length, ft, |
C | = | friction factor constant, dimensionless (Table 2), |
d | = | pipe inside diameter, in., |
Q_{l} | = | liquid flow rate, B/D, |
and | ||
gpm | = | liquid flow rate, gal/min. |
Pressure drop can be calculated from
(Eq. 10)
Pressure drop for gas flow
General equation
The general equation for calculating gas flow is stated as
(Eq. 11)
where
w | = | rate of flow, lbm/sec, |
g | = | acceleration of gravity, 32.2 ft/sec^{2}, |
A | = | cross-sectional area of pipe, ft^{2}, |
V_{1}‘ | = | specific volume of gas at upstream conditions, ft^{3}/lbm, |
f | = | friction factor, dimensionless, |
L | = | length, ft, |
D | = | diameter of the pipe, ft, |
P_{1} | = | upstream pressure, psia, |
and | ||
P_{2} | = | downstream pressure, psia. |
Assumptions: no work performed, steady-state flow, and f = constant as a function of the length.
Simplified equation
For practical pipeline purposes, Eq. 11 can be simplified to
(Eq. 12)
where
P_{1} | = | upstream pressure, psia, |
P_{2} | = | downstream pressure, psia, |
S | = | specific gravity of gas, |
Q_{g} | = | gas flow rate, MMscf/D, |
Z | = | compressibility factor for gas, dimensionless, |
T | = | flowing temperature, °R, |
f | = | Moody friction factor, dimensionless, |
d | = | pipe ID, in., |
and | ||
L | = | length, ft. |
The compressibility factor, Z, for natural gas can be found in Fig. 7.
Three simplified derivative equations can be used to calculate gas flow in pipelines:
- The Weymouth equation
- The Panhandle equation
- The Spitzglass equation
All three are effective, but the accuracy and applicability of each equation falls within certain ranges of flow and pipe diameter. The equations are stated next.
Weymouth equation
This equation is used for high-Reynolds-number flows where the Moody friction factor is merely a function of relative roughness.
(Eq. 13)
where
Q_{g} | = | gas-flow rate, MMscf/D, |
d | = | pipe inside diameter, in., |
P_{1} | = | upstream pressure, psia, |
P_{2} | = | downstream pressure, psia, |
L | = | length, ft, |
T_{1} | = | temperature of gas at inlet, °R, |
S | = | specific gravity of gas, |
and | ||
Z | = | compressibility factor for gas, dimensionless. |
Panhandle equation
This equation is used for moderate-Reynolds-number flows where the Moody friction factor is independent of relative roughness and is a function of Reynolds number to a negative power.
(Eq. 14)
where
E | = | efficiency factor (new pipe: 1.0; good operating conditions: 0.95; average operating conditions: 0.85), |
Q_{g} | = | gas-flow rate, MMscf/D, |
d | = | pipe ID, in., |
P_{1} | = | upstream pressure, psia, |
P_{2} | = | downstream pressure, psia, |
L_{m} | = | length, miles, |
T_{1} | = | temperature of gas at inlet, °R, |
S | = | specific gravity of gas, |
and | ||
Z | = | compressibility factor for gas, dimensionless. |
Spitzglass equation
Q_{g} | = | gas-flow rate, MMscf/D, |
Δh_{W} | = | pressure loss, inches of water, |
and | ||
d | = | pipe ID, in. |
Assumptions:
f | = | (1+ 3.6/ d + 0.03 d ) (1/100), |
T | = | 520°R, |
P_{1} | = | 15 psia, |
Z | = | 1.0, |
and | ||
ΔP | = | < 10% of P 1 . |
Application of the formulas
As previously discussed, there are certain conditions under which the various formulas are more applicable. A general guideline for application of the formulas is given next.
Simplified gas formula
This formula is recommended for most general-use flow applications.
Weymouth equation
The Weymouth equation is recommended for smaller-diameter pipe (generally, 12 in. and less). It is also recommended for shorter lengths of segments ( < 20 miles) within production batteries and for branch gathering lines, medium- to high-pressure (+/–100 psig to > 1,000 psig) applications, and a high Reynolds number.
Panhandle equation
This equation is recommended for larger-diameter pipe (12-in. diameter and greater). It is also recommended for long runs of pipe ( > 20 miles) such as cross-country transmission pipelines and for moderate Reynolds numbers.
Spitzglass equation
The Spitzglass equation is recommended for low-pressure vent lines < 12 in. in diameter (ΔP < 10% of P_{1}).
The petroleum engineer will find that the general gas equation and the Weymouth equation are very useful. The Weymouth equation is ideal for designing branch laterals and trunk lines in field gas-gathering systems.
Multiphase flow
Flow regimes
Fluid from the wellbore to the first piece of production equipment (separator) is generally two-phase liquid/gas flow.
The characteristics of horizontal, multiphase flow regimes are shown in Fig. 8. They can be described as follows:
- Bubble: Occurs at very low gas/liquid ratios where the gas forms bubbles that rise to the top of the pipe.
- Plug: Occurs at higher gas/liquid ratios where the gas bubbles form moderate-sized plugs.
- Stratified: As the gas/liquid ratios increase, plugs become longer until the gas and liquid flow in separate layers.
- Wavy: As the gas/liquid ratios increase further, the energy of the flowing gas stream causes waves in the flowing liquid.
- Slug: As the gas/liquid ratios continue to increase, the wave heights of the liquid increase until the crests contact the top of the pipe, creating liquid slugs.
- Spray: At extremely high gas/liquid ratios, the liquid is dispersed into the flowing-gas stream.
Fig. 9^{[1]} shows the various flow regimes that could be expected in horizontal flow as a function of the superficial velocities of gas and liquid flow. Superficial velocity is the velocity that would exist if the other phase was not present.
The multiphase flow in vertical and inclined pipe behaves somewhat differently from multiphase flow in horizontal pipe. The characteristics of the vertical flow regimes are shown in Fig. 10 and are described next.
Bubble
Where the gas/liquid ratios are small, the gas is present in the liquid in small, variable-diameter, randomly distributed bubbles. The liquid moves at a fairly uniform velocity while the bubbles move up through the liquid at differing velocities, which are dictated by the size of the bubbles. Except for the total composite-fluid density, the bubbles have little effect on the pressure gradient.
Slug flow
As the gas/liquid ratios continue to increase, the wave heights of the liquid increase until the crests contact the top of the pipe, creating liquid slugs.
Transition flow
The fluid changes from a continuous liquid phase to a continuous gas phase. The liquid slugs virtually disappear and are entrained in the gas phase. The effects of the liquid are still significant, but the effects of the gas phase are predominant.
Annular mist flow
The gas phase is continuous, and the bulk of the liquid is entrained within the gas. The liquid wets the pipe wall, but the effects of the liquid are minimal as the gas phase becomes the controlling factor. Fig. 11^{[2]} shows the various flow regimes that could be expected in vertical flow as a function of the superficial velocities of gas and liquid flow.
Two phase pressure drop
The calculation of pressure drop in two-phase flow is very complex and is based on empirical relationships to take into account the phase changes that occur because of pressure and temperature changes along the flow, the relative velocities of the phases, and complex effects of elevation changes. Table 3 lists several commercial programs that are available to model pressure drop. Because all are based to some extent on empirical relations, they are limited in accuracy to the data sets from which the relations were designed. It is not unusual for measured pressure drops in the field to differ by ± 20% from those calculated by any of these models.
Simplified friction pressure drop approximation for two phase flow
Eq. 16 provides an approximate solution for friction pressure drop in two-phase-flow problems that meet the assumptions stated.
(Eq. 16)
where
ΔP | = | friction pressure drop, psi, |
f | = | Moody friction factor, dimensionless, |
L | = | length, ft, |
W | = | rate of flow of mixture, lbm/hr, |
ρ_{M} | = | density of the mixture, lbm/ft^{3}, |
and | ||
d | = | pipe ID, in. |
The formula for rate of mixture flow is
(Eq. 17)
where
Q_{g} | = | gas-flow rate, MMscf/D, |
Q_{L} | = | liquid flow rate, B/D, |
S | = | specific gravity of gas at standard conditions, lbm/ft^{3} (air = 1), |
and | ||
SG | = | specific gravity of liquid, relative to water, lbm/ft^{3}. |
The density of the mixture is given by
(Eq. 18)
where
P | = | operating pressure, psia, |
R | = | gas/liquid ratio, ft^{3}/bbl, |
T | = | operating temperature, °R, |
SG | = | specific gravity of liquid, relative to water, lbm/ft^{3}, |
S | = | specific gravity of gas at standard conditions, lbm/ft^{3} (air = 1), |
and | ||
Z | = | gas compressibility factor, dimensionless. |
The formula is applicable if the following conditions are met:
- ΔP is less than 10% of the inlet pressure.
- Bubble or mist exists.
- There are no elevation changes.
- There is no irreversible energy transfer between phases.
Pressure Drop Because of Changes in Elevation
There are several notable characteristics associated with pressure drop because of elevation changes in two-phase flow. The flow characteristics associated with the elevation changes include:
- In downhill lines, flow becomes stratified as liquid flows faster than gas.
- The depth of the liquid layer adjusts to the static pressure head and is equal to the friction pressure drop.
- There is no pressure recovery in the downhill line.
- In low gas/liquid flow, the flow in uphill segments can be liquid "full" at low flow rates. Thus, at low flow rates, the total pressure drop is the sum of the pressure drops for all of the uphill runs.
- With increased gas flow, the total pressure drop may decrease as liquid is removed from uphill segments.
The pressure drop at low flow rates associated with an uphill elevation change may be approximated with Eq. 19.
(Eq. 19)
where
ΔP_{Z} | = | pressure drop because of elevation increase in the segment, psi, |
SG | = | specific gravity of the liquid in the segment, relative to water, |
and | ||
ΔZ | = | increase in elevation for segment, ft. |
The total pressure drop can then be approximated by the sum of the pressure drops for each uphill segment.
Pressure drop caused by valves and fittings
One of the most important parameters affecting pressure drop in piping systems is pressure loss in the fittings and valves, which is incorporated in the system. For piping systems within production facilities, the pressure drop through fittings and valves can be much greater than that through the straight run of pipe itself. In long pipeline systems, the pressure drop through fittings and valves can often be ignored.
Resistance coefficients
The head loss in valves and fittings can be calculated with resistance coefficients as
(Eq. 20)
where
H_{L} | = | head loss, ft, |
K_{r} | = | resistance coefficient, dimensionless, |
D | = | pipe ID, ft, |
and | ||
V | = | velocity, ft/sec. |
The total head loss is the sum of all K_{r} V^{2}/2g.
The resistance coefficients K_{r} for individual valves and fittings are found in tabular form in a number of industry publications. Most manufacturers publish tabular data for all sizes and configurations of their products. One of the best sources of data is the Crane Flow of Fluids, technical paper No. 410. ^{[3]} The Natural Gas Processors Suppliers Assn. (NGPSA) Engineering Data Book^{[4]} and Ingersoll-Rand’s Cameron Hydraulic Data Book^{[5]} are also good sources of references for the information. Some examples of resistance coefficients are listed in Tables 4 and 5.
Flow coefficients
The flow coefficient for liquids, C_{V}, is determined experimentally for each valve or fitting as the flow of water, in gal/min at 60°F for a pressure drop of 1 psi through the fitting. The relationship between flow and resistance coefficients can be expressed as
(Eq. 21)
In any fitting or valve with a known C_{V}, the pressure drop can be calculated for different conditions of flow and liquid properties with Eq. 22.
(Eq. 22)
where
Q_{L} | = | liquid-flow rate, B/D, |
and | ||
SG | = | liquid specific gravity relative to water. |
Again, the CV is published for most valves and fittings and can be found in Crane Flow of Fluids,^{[3]} Engineering Data Book,^{[4]} Cameron Hydraulic Data Book,^{[5]} as well as the manufacturer’s technical data.
Equivalent lengths
The head loss associated with valves and fittings can also be calculated by considering equivalent "lengths" of pipe segments for each valve and fitting. In other words, the calculated head loss caused by fluid passing through a gate valve is expressed as an additional length of pipe that is added to the actual length of pipe in calculating pressure drop.
All of the equivalent lengths caused by the valves and fittings within a pipe segment would be added together to compute the pressure drop for the pipe segment. The equivalent length, L_{e}, can be determined from the resistance coefficient, K_{r}, and the flow coefficient, C_{V}, using the formulas given next.
(Eq. 23)
K_{r} | = | resistance coefficient, dimensionless, |
D | = | diameter of the pipe, ft, |
f | = | Moody friction factor, dimensionless, |
d | = | pipe ID, in., |
and | ||
C_{V} | = | flow coefficient for liquids, dimensionless. |
Table 6 shows equivalent lengths of pipe for a variety of valves and fittings for a number of standard pipe sizes.
Nomenclature
Z | = | elevation head, ft, |
P | = | pressure, psi, |
ρ | = | density, lbm/ft^{3}, |
V | = | velocity, ft/sec, |
g | = | gravitational constant, ft/sec^{2}, |
H_{L} | = | head loss, ft. |
f | = | Moody friction factor, dimensionless, |
L | = | pipe length, ft, |
D | = | pipe diameter, ft, |
ΔP | = | pressure drop, psi, |
μ | = | viscosity, lbm/ft-sec. |
SG | = | specific gravity of liquid relative to water (water = 1), |
Q_{l} | = | liquid-flow rate, B/D, |
S | = | specific gravity of gas at standard conditions relative to air (molecular weight divided by 29), |
Q_{g} | = | gas-flow rate, MMscf/D. |
γ | = | kinematic viscosity, centistokes, |
ϕ | = | absolute viscosity, cp |
Q_{l} | = | liquid flow rate, B/D, |
w | = | rate of flow, lbm/sec |
P_{1} | = | upstream pressure, psia |
P_{2} | = | downstream pressure, psia. |
Δh_{W} | = | pressure loss, inches of water, |
W | = | rate of flow of mixture, lbm/hr, |
ρ_{M} | = | density of the mixture, lbm/ft^{3} |
P | = | operating pressure, psia, |
R | = | gas/liquid ratio, ft^{3}/bbl, |
T | = | operating temperature, °R, |
ΔP_{Z} | = | pressure drop because of elevation increase in the segment, psi, |
ΔZ | = | increase in elevation for segment, ft. |
H_{L} | = | head loss, ft, |
K_{r} | = | resistance coefficient, dimensionless |
C_{V} | = | flow coefficient for liquids, dimensionless. |
K_{r} | = | resistance coefficient, dimensionless, |
References
- ↑ ^{1.0} ^{1.1} Griffith, P. 1984. Multiphase Flow in Pipes. J Pet Technol 36 (3): 361-367. SPE-12895-PA. http://dx.doi.org/10.2118/12895-PA.
- ↑ ^{2.0} ^{2.1} Taitel, Y., Bornea, D., and Dukler, A.E. 1980. Modelling flow pattern transitions for steady upward gas-liquid flow in vertical tubes. AIChE J. 26 (3): 345-354. http://dx.doi.org/10.1002/aic.690260304.
- ↑ ^{3.0} ^{3.1} Crane Flow of Fluids, Technical Paper No. 410. 1976. New York City: Crane Manufacturing Co.
- ↑ ^{4.0} ^{4.1} Engineering Data Book, ninth edition. 1972. Tulsa, Oklahoma: Natural Gas Processors Suppliers Assn.
- ↑ ^{5.0} ^{5.1} Westway, C.R. and Loomis,A.W. ed. 1979. Cameron Hydraulic Data Book, sixteenth edition. Woodcliff Lake, New Jersey: Ingersoll-Rand.
Noteworthy papers in OnePetro
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read
External links
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro