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: Difference between revisions

PetroWiki
Jump to navigation Jump to search
No edit summary
 
Line 1: Line 1:
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''').  
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''').<br/><gallery widths="300px" heights="200px">
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 319 Image 0001.png|'''Fig. 1—Fluid-flow system (courtesy of AMEC Paragon).'''
File:Vol3 Page 319 Image 0001.png|'''Fig. 1—Fluid-flow system (courtesy of AMEC Paragon).'''
</gallery>
</gallery>


==Piping design==
== Piping design ==
The minimum basic parameters that are required to design the piping system include, but are not limited to, the following.  
 
<br>
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 characteristics and physical properties of the fluid.
* The pressure, temperature, and elevation at Point A.
*The desired mass-flow rate (or volume) of the fluid to be transported.
* The pressure, temperature, and elevation at Point B.
*The pressure, temperature, and elevation at Point A.
* The distance between Point A and Point B (or length the fluid must travel) and equivalent length (pressure losses) introduced by valves and fittings.
*The pressure, temperature, and elevation at Point B.
<br>
*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.  
 
<br/>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.<br/><br/>The equation is stated as<br/>[[File:Vol3 page 319 eq 001.PNG|RTENOTITLE]] '''(Eq. 1)'''<br/>where


==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.
<br>
<br>
The equation is stated as
<br>
[[File:Vol3_page_319_eq_001.PNG‎]] '''(Eq. 1)'''
<br>
where
{|
{|
|''Z''
|=
|elevation head, ft,
|-
|-
|''P''  
| ''Z''
|=  
| =
|pressure, psi,  
| elevation head, ft,
|-
| ''P''
| =
| pressure, psi,
|-
|-
|''ρ''  
| ''ρ''
|=  
| =
|density, lbm/ft<sup>3</sup>,  
| density, lbm/ft<sup>3</sup>,
|-
|-
|''V''  
| ''V''
|=  
| =
|velocity, ft/sec,  
| velocity, ft/sec,
|-
|-
|''g''  
| ''g''
|=  
| =
|gravitational constant, ft/sec<sup>2</sup>,  
| gravitational constant, ft/sec<sup>2</sup>,
|-
|-
|and  
| and
|-
|-
|''H''<sub>''L''</sub>  
| ''H''<sub>''L''</sub>
|=  
| =
|head loss, ft.  
| head loss, ft.
|}
|}
<br>
 
'''Fig. 2''' presents a simplified graphic illustration of the Bernoulli equation.  
<br/>'''Fig. 2''' presents a simplified graphic illustration of the Bernoulli equation.<br/><br/><gallery widths="300px" heights="200px">
<br>
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 320 Image 0001.png|'''Fig. 2—Sketch four Bernoulli equation (courtesy of AMEC Paragon).'''
File:Vol3 Page 320 Image 0001.png|'''Fig. 2—Sketch four Bernoulli equation (courtesy of AMEC Paragon).'''
</gallery>
</gallery><br/>Darcy’s equation further expresses head loss as<br/>[[File:Vol3 page 319 eq 002.PNG|RTENOTITLE]] '''(Eq. 2)'''<br/>and<br/>[[File:Vol3 page 319 eq 003.PNG|RTENOTITLE]] '''(Eq. 3)'''<br/>where
<br>
 
Darcy’s equation further expresses head loss as
<br>
[[File:Vol3_page_319_eq_002.PNG]] '''(Eq. 2)'''
<br>
and
<br>
[[File:Vol3_page_319_eq_003.PNG]] '''(Eq. 3)'''
<br>
where
{|
{|
|''H''<sub>''L''</sub>
|=
|head loss, ft,
|-
|-
|''f''  
| ''H''<sub>''L''</sub>
|=  
| =
|Moody friction factor, dimensionless,  
| head loss, ft,
|-
| ''f''
| =
| Moody friction factor, dimensionless,
|-
|-
|''L''  
| ''L''
|=  
| =
|pipe length, ft,  
| pipe length, ft,
|-
|-
|''D''  
| ''D''
|=  
| =
|pipe diameter, ft,  
| pipe diameter, ft,
|-
|-
|''V''  
| ''V''
|=  
| =
|velocity, ft/sec,  
| velocity, ft/sec,
|-
|-
|''g''  
| ''g''
|=  
| =
|gravitational constant ft/sec<sup>2</sup>,  
| gravitational constant ft/sec<sup>2</sup>,
|-
|-
|Δ''P''  
| Δ''P''
|=  
| =
|pressure drop, psi,  
| pressure drop, psi,
|-
|-
|''ρ''  
| ''ρ''
|=  
| =
|density, lbm/ft<sup>3</sup>,  
| density, lbm/ft<sup>3</sup>,
|-
|-
|and  
| and
|-
|-
|''d''  
| ''d''
|=  
| =
|pipe inside diameter, in.  
| pipe inside diameter, in.
|}
|}


== Reynolds number and Moody friction factor ==
== 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
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<br/>[[File:Vol3 page 320 eq 001.PNG|RTENOTITLE]] '''(Eq. 4)'''<br/>where
<br>
 
[[File:Vol3_page_320_eq_001.PNG]] '''(Eq. 4)'''
<br>
where
{|
{|
|''ρ''
|=
|density, lbm/ft<sup>3</sup>,
|-
|-
|''D''  
| ''ρ''
|=  
| =
|pipe internal diameter, ft,  
| density, lbm/ft<sup>3</sup>,
|-
| ''D''
| =
| pipe internal diameter, ft,
|-
|-
|''V''  
| ''V''
|=  
| =
|flow velocity, ft/sec,  
| flow velocity, ft/sec,
|-
|-
|and  
| and
|-
|-
|''μ''  
| ''μ''
|=  
| =
|viscosity, lbm/ft-sec.  
| viscosity, lbm/ft-sec.
|}
|}
<br>
 
The Reynolds number for liquids can be expressed as
<br/>The Reynolds number for liquids can be expressed as<br/>[[File:Vol3 page 320 eq 002.PNG|RTENOTITLE]] '''(Eq. 5)'''<br/>where
<br>
 
[[File:Vol3_page_320_eq_002.PNG]] '''(Eq. 5)'''
<br>
where
{|
{|
|''μ''
|=
|viscosity, cp,
|-
|-
|''d''  
| ''μ''
|=  
| =
|pipe inside diameter, in.,  
| viscosity, cp,
|-
| ''d''
| =
| pipe inside diameter, in.,
|-
|-
|''SG''  
| ''SG''
|=  
| =
|specific gravity of liquid relative to water (water = 1),  
| specific gravity of liquid relative to water (water = 1),
|-
|-
|''Q''<sub>l</sub>  
| ''Q''<sub>l</sub>
|=  
| =
|liquid-flow rate, B/D,  
| liquid-flow rate, B/D,
|-
|-
|and  
| and
|-
|-
|''V''  
| ''V''
|=  
| =
|velocity, ft/sec.  
| velocity, ft/sec.
|}
|}
<br>
 
The Reynolds number for gases can be expressed as
<br/>The Reynolds number for gases can be expressed as<br/>[[File:Vol3 page 321 eq 001.PNG|RTENOTITLE]] '''(Eq. 6)'''<br/>where
<br>
 
[[File:Vol3_page_321_eq_001.PNG]] '''(Eq. 6)'''
<br>
where
{|
{|
|''μ''
|=
|viscosity, cp,
|-
|-
|''d''  
| ''μ''
|=  
| =
|pipe inside diameter, in.,  
| viscosity, cp,
|-
| ''d''
| =
| pipe inside diameter, in.,
|-
|-
|''S''  
| ''S''
|=  
| =
|specific gravity of gas at standard conditions relative to air (molecular weight divided by 29),  
| specific gravity of gas at standard conditions relative to air (molecular weight divided by 29),
|-
|-
|and  
| and
|-
|-
|''Q''<sub>''g''</sub>  
| ''Q''<sub>''g''</sub>
|=  
| =
|gas-flow rate, MMscf/D.  
| gas-flow rate, MMscf/D.
|}
|}
<br>
 
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.
<br/>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.<br/><br/><gallery widths="300px" heights="200px">
<br>
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 322 Image 0001.png|'''Fig. 3—Friction-factor chart (courtesy of AMEC Paragon).'''
File:Vol3 Page 322 Image 0001.png|'''Fig. 3—Friction-factor chart (courtesy of AMEC Paragon).'''


File:Vol3 Page 321 Image 0001.png|'''Table 1'''
File:Vol3 Page 321 Image 0001.png|'''Table 1'''
</gallery>
</gallery><br/>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<br/>[[File:Vol3 page 321 eq 002.PNG|RTENOTITLE]] '''(Eq. 7)'''<br/>where
<br>
 
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
<br>
[[File:Vol3_page_321_eq_002.PNG]] '''(Eq. 7)'''
<br>
where
{|
{|
|=
|kinematic viscosity, centistokes,
|-
|-
|ϕ  
| γ
|=  
| =
|absolute viscosity, cp,  
| kinematic viscosity, centistokes,
|-
| ϕ
| =
| absolute viscosity, cp,
|-
|-
|and  
| and
|-
|-
|''SG''  
| ''SG''
|=  
| =
|specific gravity.  
| specific gravity.
|}
|}
<br>
 
<br>
<br/><br/><gallery widths="300px" heights="200px">
<gallery widths=300px heights=200px>
File:Vol3 Page 324 Image 0001.png|'''Fig. 4—Standard viscosity/temperature charts for liquid petroleum products (courtesy of ASTM).'''
File:Vol3 Page 324 Image 0001.png|'''Fig. 4—Standard viscosity/temperature charts for liquid petroleum products (courtesy of ASTM).'''


Line 228: Line 198:
File:Vol3 Page 327 Image 0001.png|'''Fig. 6—Hydrocarbon-gas viscosity vs. temperature (courtesy Western Supply Co.).'''
File:Vol3 Page 327 Image 0001.png|'''Fig. 6—Hydrocarbon-gas viscosity vs. temperature (courtesy Western Supply Co.).'''
</gallery>
</gallery>
<br>


==Pressure drop for liquid flow==
== Pressure drop for liquid flow ==


===General equation===  
=== General equation ===
 
'''Eq. 3''' can be expressed in terms of pipe inside diameter (ID) as stated next.<br/>[[File:Vol3 page 323 eq 001.PNG|RTENOTITLE]] '''(Eq. 8)'''<br/>where


'''Eq. 3''' can be expressed in terms of pipe inside diameter (ID) as stated next.
<br>
[[File:Vol3_page_323_eq_001.PNG]] '''(Eq. 8)'''
<br>
where
{|
{|
|''d''
|=
|pipe inside diameter, in.,
|-
|-
|''f''  
| ''d''
|=  
| =
|Moody friction factor, dimensionless,  
| pipe inside diameter, in.,
|-
|-
|''L''  
| ''f''
|=  
| =
|length of pipe, ft,  
| Moody friction factor, dimensionless,
|-
|-
|''Q''<sub>l</sub>
| ''L''
|=  
| =
|liquid flow rate, B/D,  
| length of pipe, ft,
|-
|-
|''SG''  
| ''Q''<sub>l</sub>
|=  
| =
|specific gravity of liquid relative to water,  
| liquid flow rate, B/D,
|-
|-
|and
| ''SG''
| =
| specific gravity of liquid relative to water,
|-
|-
|Δ''P''  
| and
|=  
|-
|pressure drop, psi (total pressure drop).  
| Δ''P''
| =
| pressure drop, psi (total pressure drop).
|}
|}
<br>
 
===Hazen Williams equation===
 
The Hazen-Williams equation, which is applicable only for water in turbulent flow at 60°F, expresses head loss as
 
<br>
=== Hazen Williams equation ===
[[File:Vol3_page_323_eq_002.PNG]] '''(Eq. 9)'''
 
<br>
The Hazen-Williams equation, which is applicable only for water in turbulent flow at 60°F, expresses head loss as<br/>[[File:Vol3 page 323 eq 002.PNG|RTENOTITLE]] '''(Eq. 9)'''<br/>where
where
 
{|
{|
|''H''<sub>''L''</sub>
|=
|head loss because of friction, ft,
|-
|-
|''L''  
| ''H''<sub>''L''</sub>
|=  
| =
|pipe length, ft,  
| head loss because of friction, ft,
|-
|-
|''C''  
| ''L''
|=  
| =
|friction factor constant, dimensionless ('''Table 2'''),  
| pipe length, ft,
|-
|-
|''d''  
| ''C''
|=  
| =
|pipe inside diameter, in.,  
| friction factor constant, dimensionless ('''Table 2'''),
|-
|-
|''Q''<sub>l</sub>
| ''d''
|=  
| =
|liquid flow rate, B/D,  
| pipe inside diameter, in.,
|-
|-
|and
| ''Q''<sub>l</sub>
| =
| liquid flow rate, B/D,
|-
|-
|gpm  
| and
|=  
|-
|liquid flow rate, gal/min.  
| gpm
| =
| liquid flow rate, gal/min.
|}
|}
<br>
 
<br>
<br/><br/><gallery widths="300px" heights="200px">
<gallery widths=300px heights=200px>
File:Vol3 Page 327 Image 0001.png|'''Table 2'''
File:Vol3 Page 327 Image 0001.png|'''Table 2'''
</gallery>
</gallery><br/>Pressure drop can be calculated from<br/>[[File:Vol3 page 323 eq 003.PNG|RTENOTITLE]] '''(Eq. 10)'''
<br>
 
Pressure drop can be calculated from
== Pressure drop for gas flow ==
<br>
 
[[File:Vol3_page_323_eq_003.PNG]] '''(Eq. 10)'''
=== General equation ===
<br>


==Pressure drop for gas flow==
The general equation for calculating gas flow is stated as<br/>[[File:Vol3 page 323 eq 004.PNG|RTENOTITLE]] '''(Eq. 11)'''<br/>where


===General equation===
The general equation for calculating gas flow is stated as
<br>
[[File:Vol3_page_323_eq_004.PNG]] '''(Eq. 11)'''
<br>
where
{|
{|
|''w''
|=
|rate of flow, lbm/sec,
|-
|-
|''g''  
| ''w''
|=  
| =
|acceleration of gravity, 32.2 ft/sec<sup>2</sup>,  
| rate of flow, lbm/sec,
|-
| ''g''
| =
| acceleration of gravity, 32.2 ft/sec<sup>2</sup>,
|-
|-
|''A''  
| ''A''
|=  
| =
|cross-sectional area of pipe, ft<sup>2</sup>,  
| cross-sectional area of pipe, ft<sup>2</sup>,
|-
|-
|''V''<sub>1</sub>‘  
| ''V''<sub>1</sub>‘
|=  
| =
|specific volume of gas at upstream conditions, ft<sup>3</sup>/lbm,  
| specific volume of gas at upstream conditions, ft<sup>3</sup>/lbm,
|-
|-
|''f''  
| ''f''
|=  
| =
|friction factor, dimensionless,  
| friction factor, dimensionless,
|-
|-
|''L''  
| ''L''
|=  
| =
|length, ft,  
| length, ft,
|-
|-
|''D''  
| ''D''
|=  
| =
|diameter of the pipe, ft,  
| diameter of the pipe, ft,
|-
|-
|''P''<sub>1</sub>  
| ''P''<sub>1</sub>
|=  
| =
|upstream pressure, psia,  
| upstream pressure, psia,
|-
|-
|and  
| and
|-
|-
|''P''<sub>2</sub>  
| ''P''<sub>2</sub>
|=  
| =
|downstream pressure, psia.  
| downstream pressure, psia.
|}
|}
<br>
 
Assumptions: no work performed, steady-state flow, and ''f'' = constant as a function of the length.  
<br/>Assumptions: no work performed, steady-state flow, and ''f'' = constant as a function of the length.<br/>
<br>
 
<br>
=== Simplified equation ===
===Simplified equation===
 
For practical pipeline purposes, '''Eq. 11''' can be simplified to
For practical pipeline purposes, '''Eq. 11''' can be simplified to<br/>[[File:Vol3 page 327 eq 001.PNG|RTENOTITLE]] '''(Eq. 12)'''<br/>where
<br>
 
[[File:Vol3_page_327_eq_001.PNG]] '''(Eq. 12)'''
<br>
where
{|
{|
|''P''<sub>1</sub>
|=
|upstream pressure, psia,
|-
|-
|''P''<sub>2</sub>  
| ''P''<sub>1</sub>
|=  
| =
|downstream pressure, psia,  
| upstream pressure, psia,
|-
| ''P''<sub>2</sub>
| =
| downstream pressure, psia,
|-
|-
|''S''  
| ''S''
|=  
| =
|specific gravity of gas,  
| specific gravity of gas,
|-
|-
|''Q''<sub>''g''</sub>  
| ''Q''<sub>''g''</sub>
|=  
| =
|gas flow rate, MMscf/D,  
| gas flow rate, MMscf/D,
|-
|-
|''Z''  
| ''Z''
|=  
| =
|compressibility factor for gas, dimensionless,  
| compressibility factor for gas, dimensionless,
|-
|-
|''T''  
| ''T''
|=  
| =
|flowing temperature, °R,  
| flowing temperature, °R,
|-
|-
|''f''  
| ''f''
|=  
| =
|Moody friction factor, dimensionless,  
| Moody friction factor, dimensionless,
|-
|-
|''d''  
| ''d''
|=  
| =
|pipe ID, in.,  
| pipe ID, in.,
|-
|-
|and  
| and
|-
|-
|''L''  
| ''L''
|=  
| =
|length, ft.  
| length, ft.
|}
|}
<br>
 
The compressibility factor, Z, for natural gas can be found in '''Fig. 7'''.  
<br/>The compressibility factor, Z, for natural gas can be found in '''Fig. 7'''.<br/><br/><gallery widths="300px" heights="200px">
<br>
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 328 Image 0001.png|'''Fig. 7—Compressibility of low-molecular-weight natural gases (courtesy of Natl. Gas Processors Suppliers Assn.).'''
File:Vol3 Page 328 Image 0001.png|'''Fig. 7—Compressibility of low-molecular-weight natural gases (courtesy of Natl. Gas Processors Suppliers Assn.).'''
</gallery>
</gallery><br/>Three simplified derivative equations can be used to calculate gas flow in pipelines:
<br>
 
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.<br/>
 
==== Weymouth equation ====


* The Weymouth equation
This equation is used for high-Reynolds-number flows where the Moody friction factor is merely a function of relative roughness.<br/>[[File:Vol3 page 327 eq 002.png|RTENOTITLE]] '''(Eq. 13)'''<br/>where
* 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.
<br>
<br>
====Weymouth equation====
This equation is used for high-Reynolds-number flows where the Moody friction factor is merely a function of relative roughness.
<br>
[[File:Vol3_page_327_eq_002.png]] '''(Eq. 13)'''
<br>
where
{|
{|
|''Q''<sub>''g''</sub>
|=
|gas-flow rate, MMscf/D,
|-
|-
|''d''  
| ''Q''<sub>''g''</sub>
|=  
| =
|pipe inside diameter, in.,  
| gas-flow rate, MMscf/D,
|-
| ''d''
| =
| pipe inside diameter, in.,
|-
|-
|''P''<sub>1</sub>  
| ''P''<sub>1</sub>
|=  
| =
|upstream pressure, psia,  
| upstream pressure, psia,
|-
|-
|''P''<sub>2</sub>  
| ''P''<sub>2</sub>
|=  
| =
|downstream pressure, psia,  
| downstream pressure, psia,
|-
|-
|''L''  
| ''L''
|=  
| =
|length, ft,  
| length, ft,
|-
|-
|''T''<sub>1</sub>  
| ''T''<sub>1</sub>
|=  
| =
|temperature of gas at inlet, °R,  
| temperature of gas at inlet, °R,
|-
|-
|''S''  
| ''S''
|=  
| =
|specific gravity of gas,  
| specific gravity of gas,
|-
|-
|and  
| and
|-
|-
|''Z''  
| ''Z''
|=  
| =
|compressibility factor for gas, dimensionless.  
| compressibility factor for gas, dimensionless.
|}
|}
<br>


====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.
 
<br>
==== Panhandle equation ====
[[File:Vol3_page_330_eq_001.PNG]] '''(Eq. 14)'''
 
<br>
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.<br/>[[File:Vol3 page 330 eq 001.PNG|RTENOTITLE]] '''(Eq. 14)'''<br/>where
where
 
{|
{|
|''E''
|=
|efficiency factor (new pipe: 1.0; good operating conditions: 0.95; average operating conditions: 0.85),
|-
|-
|''Q''<sub>''g''</sub>  
| ''E''
|=  
| =
|gas-flow rate, MMscf/D,  
| efficiency factor (new pipe: 1.0; good operating conditions: 0.95; average operating conditions: 0.85),
|-
| ''Q''<sub>''g''</sub>
| =
| gas-flow rate, MMscf/D,
|-
|-
|''d''  
| ''d''
|=  
| =
|pipe ID, in.,  
| pipe ID, in.,
|-
|-
|''P''<sub>1</sub>  
| ''P''<sub>1</sub>
|=  
| =
|upstream pressure, psia,  
| upstream pressure, psia,
|-
|-
|''P''<sub>2</sub>  
| ''P''<sub>2</sub>
|=  
| =
|downstream pressure, psia,  
| downstream pressure, psia,
|-
|-
|''L''<sub>''m''</sub>  
| ''L''<sub>''m''</sub>
|=  
| =
|length, miles,  
| length, miles,
|-
|-
|''T''<sub>1</sub>  
| ''T''<sub>1</sub>
|=  
| =
|temperature of gas at inlet, °R,  
| temperature of gas at inlet, °R,
|-
|-
|''S''  
| ''S''
|=  
| =
|specific gravity of gas,  
| specific gravity of gas,
|-
|-
|and  
| and
|-
|-
|''Z''  
| ''Z''
|=  
| =
|compressibility factor for gas, dimensionless.  
| compressibility factor for gas, dimensionless.
|}
|}
<br>
 
====Spitzglass equation====
 
<br>
 
[[File:Vol3_page_330_eq_002.PNG]] '''(Eq. 15)'''
==== Spitzglass equation ====
<br>
 
where
<br/>[[File:Vol3 page 330 eq 002.PNG|RTENOTITLE]] '''(Eq. 15)'''<br/>where
 
{|
{|
|''Q''<sub>''g''</sub>
|=
|gas-flow rate, MMscf/D,
|-
|-
|Δ''h''<sub>''W''</sub>  
| ''Q''<sub>''g''</sub>
|=  
| =
|pressure loss, inches of water,  
| gas-flow rate, MMscf/D,
|-
|-
|and
| Δ''h''<sub>''W''</sub>
| =
| pressure loss, inches of water,
|-
|-
|''d''  
| and
|=  
|-
|pipe ID, in.  
| ''d''
| =
| pipe ID, in.
|}
|}
<br>
 
Assumptions:  
<br/>Assumptions:
 
{|
{|
|''f''
|=
|(1+ 3.6/ d + 0.03 d ) (1/100),
|-
|-
|''T''  
| ''f''
|=  
| =
|520°R,  
| (1+ 3.6/ d + 0.03 d ) (1/100),
|-
| ''T''
| =
| 520°R,
|-
|-
|''P''<sub>1</sub>  
| ''P''<sub>1</sub>
|=  
| =
|15 psia,  
| 15 psia,
|-
|-
|''Z''  
| ''Z''
|=  
| =
|1.0,  
| 1.0,
|-
|-
|and  
| and
|-
|-
|Δ''P''  
| Δ''P''
|=  
| =
|< 10% of P 1 .  
| < 10% of P 1 .
|}
|}
<br>
===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.
<br>
<br>
====Simplified gas formula====
This formula is recommended for most general-use flow applications.
<br>
<br>
====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.
<br>
<br>
====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.
<br>
<br>
====Spitzglass equation====
The Spitzglass equation is recommended for low-pressure vent lines < 12 in. in diameter (Δ''P'' < 10% of ''P''<sub>1</sub>).
<br>
<br>
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.  
=== Application of the formulas ===
<br>
 
<br>
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.<br/>
The characteristics of horizontal, multiphase flow regimes are shown in '''Fig. 8'''. They can be described as follows:  
 
<br>
==== Simplified gas formula ====
* '''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.
This formula is recommended for most general-use flow applications.<br/>
* '''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.
==== Weymouth equation ====
* '''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.
The Weymouth equation is recommended for smaller-diameter pipe (generally, 12 in. and less). It is also recommended for shorter lengths of segments ( &lt; 20 miles) within production batteries and for branch gathering lines, medium- to high-pressure (+/–100 psig to &gt; 1,000 psig) applications, and a high Reynolds number.<br/>
<br>
 
<br>
==== Panhandle equation ====
<gallery widths=300px heights=200px>
 
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.<br/>
 
==== Spitzglass equation ====
 
The Spitzglass equation is recommended for low-pressure vent lines < 12 in. in diameter (Δ''P'' < 10% of ''P''<sub>1</sub>).<br/><br/>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.<br/><br/>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.
 
<br/><br/><gallery widths="300px" heights="200px">
File:Vol3 Page 331 Image 0001.png|'''Fig. 8—Two-phase-flow patterns in horizontal flow (courtesy of AMEC Paragon).'''
File:Vol3 Page 331 Image 0001.png|'''Fig. 8—Two-phase-flow patterns in horizontal flow (courtesy of AMEC Paragon).'''
</gallery>
</gallery><br/>'''Fig. 9'''<ref name="r1">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.</ref> 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.<br/><br/><gallery widths="300px" heights="200px">
<br>
'''Fig. 9'''<ref name="r1" /> 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.  
<br>
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 332 Image 0001.png|'''Fig. 9—Horizontal multiphase-flow map (after Griffith).'''<ref name="r1" />
File:Vol3 Page 332 Image 0001.png|'''Fig. 9—Horizontal multiphase-flow map (after Griffith).'''<ref name="r1" />
</gallery>
</gallery><br/>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.<br/><br/><gallery widths="300px" heights="200px">
<br>
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.  
<br>
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 333 Image 0001.png|'''Fig. 10—Two-phase-flow patterns in vertical flow (courtesy of AMEC Paragon).'''
File:Vol3 Page 333 Image 0001.png|'''Fig. 10—Two-phase-flow patterns in vertical flow (courtesy of AMEC Paragon).'''
</gallery>
</gallery>
<br>
 
====Bubble====
==== 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.  
 
<br>
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.<br/>
<br>
 
====Slug flow====
==== 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.  
 
<br>
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.<br/>
<br>
 
====Transition flow====
==== 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.  
 
<br>
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.<br/>
<br>
 
====Annular mist flow====
==== 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'''<ref name="r2" /> shows the various flow regimes that could be expected in vertical flow as a function of the superficial velocities of gas and liquid flow.  
 
<br>
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'''<ref name="r2">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.</ref> shows the various flow regimes that could be expected in vertical flow as a function of the superficial velocities of gas and liquid flow.<br/><br/><gallery widths="300px" heights="200px">
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 334 Image 0001.png|'''Fig. 11—Vertical-multiphase-flow map (after Taitel ''et al.'').'''<ref name="r2" />
File:Vol3 Page 334 Image 0001.png|'''Fig. 11—Vertical-multiphase-flow map (after Taitel ''et al.'').'''<ref name="r2" />
</gallery>
</gallery>
<br>


===Two phase pressure drop===
=== 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.
 
<br>
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.<br/><br/><gallery widths="300px" heights="200px">
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 335 Image 0001.png|'''Table 3'''
File:Vol3 Page 335 Image 0001.png|'''Table 3'''
</gallery>
</gallery>
<br>
 
===Simplified friction pressure drop approximation for two phase flow===
=== 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.
 
<br>
'''Eq. 16''' provides an approximate solution for friction pressure drop in two-phase-flow problems that meet the assumptions stated.<br/>[[File:Vol3 page 333 eq 001.PNG|RTENOTITLE]] '''(Eq. 16)'''<br/>where
[[File:Vol3_page_333_eq_001.PNG]] '''(Eq. 16)'''
 
<br>
where
{|
{|
|Δ''P''
|=
|friction pressure drop, psi,
|-
|-
|''f''  
| Δ''P''
|=  
| =
|Moody friction factor, dimensionless,  
| friction pressure drop, psi,
|-
|-
|''L''  
| ''f''
|=  
| =
|length, ft,  
| Moody friction factor, dimensionless,
|-
|-
|''W''  
| ''L''
|=  
| =
|rate of flow of mixture, lbm/hr,  
| length, ft,
|-
|-
|''ρ''<sub>''M''</sub>
| ''W''
|=  
| =
|density of the mixture, lbm/ft<sup>3</sup>,  
| rate of flow of mixture, lbm/hr,
|-
|-
|and
| ''ρ''<sub>''M''</sub>
| =
| density of the mixture, lbm/ft<sup>3</sup>,
|-
|-
|''d''  
| and
|=  
|-
|pipe ID, in.  
| ''d''
| =
| pipe ID, in.
|}
|}
<br>
 
The formula for rate of mixture flow is
<br/>The formula for rate of mixture flow is<br/>[[File:Vol3 page 333 eq 002.PNG|RTENOTITLE]] '''(Eq. 17)'''<br/>where
<br>
 
[[File:Vol3_page_333_eq_002.PNG]] '''(Eq. 17)'''
<br>
where
{|
{|
|''Q''<sub>''g''</sub>
|=
|gas-flow rate, MMscf/D,
|-
|-
|''Q''<sub>''L''</sub>  
| ''Q''<sub>''g''</sub>
|=  
| =
|liquid flow rate, B/D,  
| gas-flow rate, MMscf/D,
|-
| ''Q''<sub>''L''</sub>
| =
| liquid flow rate, B/D,
|-
|-
|''S''  
| ''S''
|=  
| =
|specific gravity of gas at standard conditions, lbm/ft<sup>3</sup> (air = 1),  
| specific gravity of gas at standard conditions, lbm/ft<sup>3</sup> (air = 1),
|-
|-
|and  
| and
|-
|-
|''SG''  
| ''SG''
|=  
| =
|specific gravity of liquid, relative to water, lbm/ft<sup>3</sup>.  
| specific gravity of liquid, relative to water, lbm/ft<sup>3</sup>.
|}
|}
<br>
 
The density of the mixture is given by
<br/>The density of the mixture is given by<br/>[[File:Vol3 page 334 eq 001.PNG|RTENOTITLE]] '''(Eq. 18)'''<br/>where
<br>
 
[[File:Vol3_page_334_eq_001.PNG]] '''(Eq. 18)'''
<br>
where
{|
{|
|''P''
|=
|operating pressure, psia,
|-
|-
|''R''  
| ''P''
|=  
| =
|gas/liquid ratio, ft<sup>3</sup>/bbl,  
| operating pressure, psia,
|-
| ''R''
| =
| gas/liquid ratio, ft<sup>3</sup>/bbl,
|-
|-
|''T''  
| ''T''
|=  
| =
|operating temperature, °R,  
| operating temperature, °R,
|-
|-
|''SG''  
| ''SG''
|=  
| =
|specific gravity of liquid, relative to water, lbm/ft<sup>3</sup>,  
| specific gravity of liquid, relative to water, lbm/ft<sup>3</sup>,
|-
|-
|''S''  
| ''S''
|=  
| =
|specific gravity of gas at standard conditions, lbm/ft<sup>3</sup> (air = 1),  
| specific gravity of gas at standard conditions, lbm/ft<sup>3</sup> (air = 1),
|-
|-
|and  
| and
|-
|-
|''Z''  
| ''Z''
|=  
| =
|gas compressibility factor, dimensionless.  
| gas compressibility factor, dimensionless.
|}
|}
<br>
The formula is applicable if the following conditions are met:
<br>
* Δ''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.


<br>
<br/>The formula is applicable if the following conditions are met:
===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:  
*Δ''P'' is less than 10% of the inlet pressure.
<br>
*Bubble or mist exists.
* In downhill lines, flow becomes stratified as liquid flows faster than gas.
*There are no elevation changes.
* The depth of the liquid layer adjusts to the static pressure head and is equal to the friction pressure drop.
*There is no irreversible energy transfer between phases.
* 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.
 
=== 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.
 
<br/>The pressure drop at low flow rates associated with an uphill elevation change may be approximated with '''Eq. 19'''.<br/>[[File:Vol3 page 335 eq 001.PNG|RTENOTITLE]] '''(Eq. 19)'''<br/>where


<br>
The pressure drop at low flow rates associated with an uphill elevation change may be approximated with '''Eq. 19'''.
<br>
[[File:Vol3_page_335_eq_001.PNG]] '''(Eq. 19)'''
<br>
where
{|
{|
|Δ''P''<sub>''Z''</sub>
|=
|pressure drop because of elevation increase in the segment, psi,
|-
|-
|''SG''  
| Δ''P''<sub>''Z''</sub>
|=  
| =
|specific gravity of the liquid in the segment, relative to water,  
| pressure drop because of elevation increase in the segment, psi,
|-
|-
|and
| ''SG''
| =
| specific gravity of the liquid in the segment, relative to water,
|-
|-
|Δ''Z''  
| and
|=  
|-
|increase in elevation for segment, ft.  
| Δ''Z''
| =
| increase in elevation for segment, ft.
|}
|}
<br>
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==
<br/>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.


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 ===


===Resistance coefficients===
The head loss in valves and fittings can be calculated with resistance coefficients as<br/>[[File:Vol3 page 336 eq 001.PNG|RTENOTITLE]] '''(Eq. 20)'''<br/>where


The head loss in valves and fittings can be calculated with resistance coefficients as
<br>
[[File:Vol3_page_336_eq_001.PNG]] '''(Eq. 20)'''
<br>
where
{|
{|
|''H''<sub>''L''</sub>
|=
|head loss, ft,
|-
|-
|''K''<sub>''r''</sub>  
| ''H''<sub>''L''</sub>
|=  
| =
|resistance coefficient, dimensionless,  
| head loss, ft,
|-
| ''K''<sub>''r''</sub>
| =
| resistance coefficient, dimensionless,
|-
|-
|''D''  
| ''D''
|=  
| =
|pipe ID, ft,  
| pipe ID, ft,
|-
|-
|and  
| and
|-
|-
|''V''  
| ''V''
|=  
| =
|velocity, ft/sec.  
| velocity, ft/sec.
|}
|}
<br>
 
The total head loss is the sum of all ''K''<sub>''r''</sub> ''V''<sup>2</sup>/2''g''.  
<br/>The total head loss is the sum of all ''K''<sub>''r''</sub> ''V''<sup>2</sup>/2''g''.<br/><br/>The resistance coefficients ''K''<sub>''r''</sub> 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. <ref name="r3">Crane Flow of Fluids, Technical Paper No. 410. 1976. New York City: Crane Manufacturing Co.</ref> The Natural Gas Processors Suppliers Assn. (NGPSA) ''Engineering Data Book''<ref name="r4">Engineering Data Book, ninth edition. 1972. Tulsa, Oklahoma: Natural Gas Processors Suppliers Assn.</ref> and Ingersoll-Rand’s ''Cameron Hydraulic Data Book''<ref name="r5">Westway, C.R. and Loomis,A.W. ed. 1979. Cameron Hydraulic Data Book, sixteenth edition. Woodcliff Lake, New Jersey: Ingersoll-Rand.</ref> are also good sources of references for the information. Some examples of resistance coefficients are listed in '''Tables 4 and 5'''.<br/><br/><gallery widths="300px" heights="200px">
<br>
<br>
The resistance coefficients ''K''<sub>''r''</sub> 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. <ref name="r3" /> The Natural Gas Processors Suppliers Assn. (NGPSA) ''Engineering Data Book''<ref name="r4" /> and Ingersoll-Rand’s ''Cameron Hydraulic Data Book''<ref name="r5" /> are also good sources of references for the information. Some examples of resistance coefficients are listed in '''Tables 4 and 5'''.
<br>
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 336 Image 0001.png|'''Table 4'''
File:Vol3 Page 336 Image 0001.png|'''Table 4'''


Line 825: Line 749:
File:Vol3 Page 340 Image 0001.png|'''Table 5 (Cont'd)'''
File:Vol3 Page 340 Image 0001.png|'''Table 5 (Cont'd)'''
</gallery>
</gallery>
<br>


===Flow coefficients===
=== Flow coefficients ===
 
The flow coefficient for liquids, ''C''<sub>''V''</sub>, 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<br/>[[File:Vol3 page 336 eq 002.PNG|RTENOTITLE]] '''(Eq. 21)'''<br/>In any fitting or valve with a known ''C''<sub>''V''</sub>, the pressure drop can be calculated for different conditions of flow and liquid properties with '''Eq. 22'''.<br/>[[File:Vol3 page 336 eq 003.PNG|RTENOTITLE]] '''(Eq. 22)'''<br/>where


The flow coefficient for liquids, ''C''<sub>''V''</sub>, 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
<br>
[[File:Vol3_page_336_eq_002.PNG]] '''(Eq. 21)'''
<br>
In any fitting or valve with a known ''C''<sub>''V''</sub>, the pressure drop can be calculated for different conditions of flow and liquid properties with '''Eq. 22'''.
<br>
[[File:Vol3_page_336_eq_003.PNG]] '''(Eq. 22)'''
<br>
where
{|
{|
|''Q''<sub>''L''</sub>
|=
|liquid-flow rate, B/D,
|-
|-
|and  
| ''Q''<sub>''L''</sub>
| =
| liquid-flow rate, B/D,
|-
| and
|-
|-
|''SG''  
| ''SG''
|=  
| =
|liquid specific gravity relative to water.  
| liquid specific gravity relative to water.
|}
|}
<br>
Again, the CV is published for most valves and fittings and can be found in Crane Flow of Fluids,<ref name="r3" /> Engineering Data Book,<ref name="r4" /> Cameron Hydraulic Data Book,<ref name="r5" /> as well as the manufacturer’s technical data.


===Equivalent lengths===
<br/>Again, the CV is published for most valves and fittings and can be found in Crane Flow of Fluids,<ref name="r3">Crane Flow of Fluids, Technical Paper No. 410. 1976. New York City: Crane Manufacturing Co.</ref> Engineering Data Book,<ref name="r4">Engineering Data Book, ninth edition. 1972. Tulsa, Oklahoma: Natural Gas Processors Suppliers Assn.</ref> Cameron Hydraulic Data Book,<ref name="r5">Westway, C.R. and Loomis,A.W. ed. 1979. Cameron Hydraulic Data Book, sixteenth edition. Woodcliff Lake, New Jersey: Ingersoll-Rand.</ref> 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.<br/><br/>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''<sub>''e''</sub>, can be determined from the resistance coefficient, ''K''<sub>''r''</sub>, and the flow coefficient, ''C''<sub>''V''</sub>, using the formulas given next.<br/>[[File:Vol3 page 341 eq 001.PNG|RTENOTITLE]] '''(Eq. 23)'''


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.
[[File:Vol3 page 341 eq 002.PNG|RTENOTITLE]] '''(Eq. 24)'''<br/>and<br/>[[File:Vol3 page 341 eq 003.PNG|RTENOTITLE]] '''(Eq. 25)'''<br/>where
<br>
<br>
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''<sub>''e''</sub>, can be determined from the resistance coefficient, ''K''<sub>''r''</sub>, and the flow coefficient, ''C''<sub>''V''</sub>, using the formulas given next.
<br>
[[File:Vol3_page_341_eq_001.PNG]] '''(Eq. 23)'''


[[File:Vol3_page_341_eq_002.PNG]] '''(Eq. 24)'''
<br>
and
<br>
[[File:Vol3_page_341_eq_003.PNG]] '''(Eq. 25)'''
<br>
where
{|
{|
|''K''<sub>''r''</sub>
|=
|resistance coefficient, dimensionless,
|-
|-
|''D''  
| ''K''<sub>''r''</sub>
|=  
| =
|diameter of the pipe, ft,  
| resistance coefficient, dimensionless,
|-
| ''D''
| =
| diameter of the pipe, ft,
|-
|-
|''f''  
| ''f''
|=  
| =
|Moody friction factor, dimensionless,  
| Moody friction factor, dimensionless,
|-
|-
|''d''  
| ''d''
|=  
| =
|pipe ID, in.,  
| pipe ID, in.,
|-
|-
|and  
| and
|-
|-
|''C''<sub>''V''</sub>  
| ''C''<sub>''V''</sub>
|=  
| =
|flow coefficient for liquids, dimensionless.  
| flow coefficient for liquids, dimensionless.
|}
|}
<br>
 
'''Table 6''' shows equivalent lengths of pipe for a variety of valves and fittings for a number of standard pipe sizes.
<br/>'''Table 6''' shows equivalent lengths of pipe for a variety of valves and fittings for a number of standard pipe sizes.<br/><br/><gallery widths="300px" heights="200px">
<br>
<br>
<gallery widths=300px heights=200px>
File:Vol3 Page 342 Image 0001.png|'''Table 6'''
File:Vol3 Page 342 Image 0001.png|'''Table 6'''
</gallery>
</gallery>


==Nomenclature==
== Nomenclature ==
 
{|
{|
|''Z''
|=
|elevation head, ft,
|-
|-
|''P''  
| ''Z''
|=  
| =
|pressure, psi,  
| elevation head, ft,
|-
| ''P''
| =
| pressure, psi,
|-
|-
|''ρ''  
| ''ρ''
|=  
| =
|density, lbm/ft<sup>3</sup>,  
| density, lbm/ft<sup>3</sup>,
|-
|-
|''V''  
| ''V''
|=  
| =
|velocity, ft/sec,  
| velocity, ft/sec,
|-
|-
|''g''  
| ''g''
|=  
| =
|gravitational constant, ft/sec<sup>2</sup>,  
| gravitational constant, ft/sec<sup>2</sup>,
|-
|-
|''H''<sub>''L''</sub>  
| ''H''<sub>''L''</sub>
|=  
| =
|head loss, ft.  
| head loss, ft.
|-
|-
|''f''  
| ''f''
|=  
| =
|Moody friction factor, dimensionless,  
| Moody friction factor, dimensionless,
|-
|-
|''L''  
| ''L''
|=  
| =
|pipe length, ft,  
| pipe length, ft,
|-
|-
|''D''  
| ''D''
|=  
| =
|pipe diameter, ft,  
| pipe diameter, ft,
|-
|-
|Δ''P''  
| Δ''P''
|=  
| =
|pressure drop, psi,  
| pressure drop, psi,
|-
|-
|''μ''  
| ''μ''
|=  
| =
|viscosity, lbm/ft-sec.
| viscosity, lbm/ft-sec.
|-
|-
|''SG''  
| ''SG''
|=  
| =
|specific gravity of liquid relative to water (water = 1),  
| specific gravity of liquid relative to water (water = 1),
|-
|-
|''Q''<sub>l</sub>  
| ''Q''<sub>l</sub>
|=  
| =
|liquid-flow rate, B/D,  
| liquid-flow rate, B/D,
|-
|-
|''S''  
| ''S''
|=  
| =
|specific gravity of gas at standard conditions relative to air (molecular weight divided by 29),  
| specific gravity of gas at standard conditions relative to air (molecular weight divided by 29),
|-
|-
|''Q''<sub>''g''</sub>  
| ''Q''<sub>''g''</sub>
|=  
| =
|gas-flow rate, MMscf/D.
| gas-flow rate, MMscf/D.
|-
|-
|γ  
| γ
|=  
| =
|kinematic viscosity, centistokes,  
| kinematic viscosity, centistokes,
|-
|-
|ϕ  
| ϕ
|=  
| =
|absolute viscosity, cp
| absolute viscosity, cp
|-
|-
|''Q''<sub>l</sub>  
| ''Q''<sub>l</sub>
|=  
| =
|liquid flow rate, B/D,  
| liquid flow rate, B/D,
|-
|-
|''w''  
| ''w''
|=  
| =
|rate of flow, lbm/sec
| rate of flow, lbm/sec
|-
|-
|''P''<sub>1</sub>  
| ''P''<sub>1</sub>
|=  
| =
|upstream pressure, psia
| upstream pressure, psia
|-
|-
|''P''<sub>2</sub>  
| ''P''<sub>2</sub>
|=  
| =
|downstream pressure, psia.  
| downstream pressure, psia.
|-
|-
|Δ''h''<sub>''W''</sub>  
| Δ''h''<sub>''W''</sub>
|=  
| =
|pressure loss, inches of water,  
| pressure loss, inches of water,
|-
|-
|''W''  
| ''W''
|=  
| =
|rate of flow of mixture, lbm/hr,  
| rate of flow of mixture, lbm/hr,
|-
|-
|''ρ''<sub>''M''</sub>  
| ''ρ''<sub>''M''</sub>
|=  
| =
|density of the mixture, lbm/ft<sup>3</sup>
| density of the mixture, lbm/ft<sup>3</sup>
|-
|-
|''P''  
| ''P''
|=  
| =
|operating pressure, psia,  
| operating pressure, psia,
|-
|-
|''R''  
| ''R''
|=  
| =
|gas/liquid ratio, ft<sup>3</sup>/bbl,  
| gas/liquid ratio, ft<sup>3</sup>/bbl,
|-
|-
|''T''  
| ''T''
|=  
| =
|operating temperature, °R,  
| operating temperature, °R,
|-
|-
|Δ''P''<sub>''Z''</sub>  
| Δ''P''<sub>''Z''</sub>
|=  
| =
|pressure drop because of elevation increase in the segment, psi,  
| pressure drop because of elevation increase in the segment, psi,
|-
|-
|Δ''Z''  
| Δ''Z''
|=  
| =
|increase in elevation for segment, ft.
| increase in elevation for segment, ft.
|-
|-
|''H''<sub>''L''</sub>  
| ''H''<sub>''L''</sub>
|=  
| =
|head loss, ft,  
| head loss, ft,
|-
|-
|''K''<sub>''r''</sub>  
| ''K''<sub>''r''</sub>
|=  
| =
|resistance coefficient, dimensionless
| resistance coefficient, dimensionless
|-
|-
|''C''<sub>''V''</sub>  
| ''C''<sub>''V''</sub>
|=  
| =
|flow coefficient for liquids, dimensionless.
| flow coefficient for liquids, dimensionless.
|-
|-
|''K''<sub>''r''</sub>  
| ''K''<sub>''r''</sub>
|=  
| =
|resistance coefficient, dimensionless,
| resistance coefficient, dimensionless,
|}
|}


==References==
== References ==
<references>
 
<ref name="r1">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.</ref>
<references />
<ref name="r2">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.</ref>
 
<ref name="r3">Crane Flow of Fluids, Technical Paper No. 410. 1976. New York City: Crane Manufacturing Co.</ref>
== Noteworthy papers in OnePetro ==
<ref name="r4">''Engineering Data Book'', ninth edition. 1972. Tulsa, Oklahoma: Natural Gas Processors Suppliers Assn.</ref>
<ref name="r5">Westway, C.R. and Loomis,A.W. ed. 1979. ''Cameron Hydraulic Data Book'', sixteenth edition. Woodcliff Lake, New Jersey: Ingersoll-Rand.</ref>
</references>


==Noteworthy papers in OnePetro==
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read


==External links==
== External links ==
 
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro


==See also==
== See also ==
[[Piping and pipeline systems]]
 
[[Piping_and_pipeline_systems|Piping and pipeline systems]]
 
[[Pipelines|Pipelines]]


[[Pipelines]]
[[Pipeline_pigging|Pipeline pigging]]


[[Pipeline pigging]]
[[Pipeline_design_consideration_and_standards|Pipeline design consideration and standards]]


[[Pipeline design consideration and standards]]
[[PEH:Piping_and_Pipelines]]


[[PEH:Piping and Pipelines]]
[[Category:5.3.2 Multiphase flow]] [[Category:4.2.1 Piping design and simulation]] [[Category:4.2 Pipelines, flowlines, and risers]]

Latest revision as of 10:38, 8 June 2015

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).

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
RTENOTITLE (Eq. 1)
where

Z = elevation head, ft,
P = pressure, psi,
ρ = density, lbm/ft3,
V = velocity, ft/sec,
g = gravitational constant, ft/sec2,
and
HL = head loss, ft.


Fig. 2 presents a simplified graphic illustration of the Bernoulli equation.


Darcy’s equation further expresses head loss as
RTENOTITLE (Eq. 2)
and
RTENOTITLE (Eq. 3)
where

HL = head loss, ft,
f = Moody friction factor, dimensionless,
L = pipe length, ft,
D = pipe diameter, ft,
V = velocity, ft/sec,
g = gravitational constant ft/sec2,
ΔP = pressure drop, psi,
ρ = density, lbm/ft3,
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
RTENOTITLE (Eq. 4)
where

ρ = density, lbm/ft3,
D = pipe internal diameter, ft,
V = flow velocity, ft/sec,
and
μ = viscosity, lbm/ft-sec.


The Reynolds number for liquids can be expressed as
RTENOTITLE (Eq. 5)
where

μ = viscosity, cp,
d = pipe inside diameter, in.,
SG = specific gravity of liquid relative to water (water = 1),
Ql = liquid-flow rate, B/D,
and
V = velocity, ft/sec.


The Reynolds number for gases can be expressed as
RTENOTITLE (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
Qg = 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
RTENOTITLE (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.
RTENOTITLE (Eq. 8)
where

d = pipe inside diameter, in.,
f = Moody friction factor, dimensionless,
L = length of pipe, ft,
Ql = 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
RTENOTITLE (Eq. 9)
where

HL = head loss because of friction, ft,
L = pipe length, ft,
C = friction factor constant, dimensionless (Table 2),
d = pipe inside diameter, in.,
Ql = liquid flow rate, B/D,
and
gpm = liquid flow rate, gal/min.




Pressure drop can be calculated from
RTENOTITLE (Eq. 10)

Pressure drop for gas flow

General equation

The general equation for calculating gas flow is stated as
RTENOTITLE (Eq. 11)
where

w = rate of flow, lbm/sec,
g = acceleration of gravity, 32.2 ft/sec2,
A = cross-sectional area of pipe, ft2,
V1 = specific volume of gas at upstream conditions, ft3/lbm,
f = friction factor, dimensionless,
L = length, ft,
D = diameter of the pipe, ft,
P1 = upstream pressure, psia,
and
P2 = 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
RTENOTITLE (Eq. 12)
where

P1 = upstream pressure, psia,
P2 = downstream pressure, psia,
S = specific gravity of gas,
Qg = 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.
RTENOTITLE (Eq. 13)
where

Qg = gas-flow rate, MMscf/D,
d = pipe inside diameter, in.,
P1 = upstream pressure, psia,
P2 = downstream pressure, psia,
L = length, ft,
T1 = 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.
RTENOTITLE (Eq. 14)
where

E = efficiency factor (new pipe: 1.0; good operating conditions: 0.95; average operating conditions: 0.85),
Qg = gas-flow rate, MMscf/D,
d = pipe ID, in.,
P1 = upstream pressure, psia,
P2 = downstream pressure, psia,
Lm = length, miles,
T1 = temperature of gas at inlet, °R,
S = specific gravity of gas,
and
Z = compressibility factor for gas, dimensionless.


Spitzglass equation


RTENOTITLE (Eq. 15)
where

Qg = gas-flow rate, MMscf/D,
ΔhW = pressure loss, inches of water,
and
d = pipe ID, in.


Assumptions:

f = (1+ 3.6/ d + 0.03 d ) (1/100),
T = 520°R,
P1 = 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 P1).

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.
RTENOTITLE (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/ft3,
and
d = pipe ID, in.


The formula for rate of mixture flow is
RTENOTITLE (Eq. 17)
where

Qg = gas-flow rate, MMscf/D,
QL = liquid flow rate, B/D,
S = specific gravity of gas at standard conditions, lbm/ft3 (air = 1),
and
SG = specific gravity of liquid, relative to water, lbm/ft3.


The density of the mixture is given by
RTENOTITLE (Eq. 18)
where

P = operating pressure, psia,
R = gas/liquid ratio, ft3/bbl,
T = operating temperature, °R,
SG = specific gravity of liquid, relative to water, lbm/ft3,
S = specific gravity of gas at standard conditions, lbm/ft3 (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.
RTENOTITLE (Eq. 19)
where

ΔPZ = 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
RTENOTITLE (Eq. 20)
where

HL = head loss, ft,
Kr = resistance coefficient, dimensionless,
D = pipe ID, ft,
and
V = velocity, ft/sec.


The total head loss is the sum of all Kr V2/2g.

The resistance coefficients Kr 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, CV, 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
RTENOTITLE (Eq. 21)
In any fitting or valve with a known CV, the pressure drop can be calculated for different conditions of flow and liquid properties with Eq. 22.
RTENOTITLE (Eq. 22)
where

QL = 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, Le, can be determined from the resistance coefficient, Kr, and the flow coefficient, CV, using the formulas given next.
RTENOTITLE (Eq. 23)

RTENOTITLE (Eq. 24)
and
RTENOTITLE (Eq. 25)
where

Kr = resistance coefficient, dimensionless,
D = diameter of the pipe, ft,
f = Moody friction factor, dimensionless,
d = pipe ID, in.,
and
CV = 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/ft3,
V = velocity, ft/sec,
g = gravitational constant, ft/sec2,
HL = 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),
Ql = liquid-flow rate, B/D,
S = specific gravity of gas at standard conditions relative to air (molecular weight divided by 29),
Qg = gas-flow rate, MMscf/D.
γ = kinematic viscosity, centistokes,
ϕ = absolute viscosity, cp
Ql = liquid flow rate, B/D,
w = rate of flow, lbm/sec
P1 = upstream pressure, psia
P2 = downstream pressure, psia.
ΔhW = pressure loss, inches of water,
W = rate of flow of mixture, lbm/hr,
ρM = density of the mixture, lbm/ft3
P = operating pressure, psia,
R = gas/liquid ratio, ft3/bbl,
T = operating temperature, °R,
ΔPZ = pressure drop because of elevation increase in the segment, psi,
ΔZ = increase in elevation for segment, ft.
HL = head loss, ft,
Kr = resistance coefficient, dimensionless
CV = flow coefficient for liquids, dimensionless.
Kr = resistance coefficient, dimensionless,

References

  1. 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. 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. 3.0 3.1 Crane Flow of Fluids, Technical Paper No. 410. 1976. New York City: Crane Manufacturing Co.
  4. 4.0 4.1 Engineering Data Book, ninth edition. 1972. Tulsa, Oklahoma: Natural Gas Processors Suppliers Assn.
  5. 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

See also

Piping and pipeline systems

Pipelines

Pipeline pigging

Pipeline design consideration and standards

PEH:Piping_and_Pipelines