Fluid Flow Manual PDF
Fluid Flow Manual PDF
Fluid Flow Manual PDF
March 1997
Manual sponsor:
For information or help regarding this manual, contact R.P. (Rob) Hohmann 242-2216
Printing History
Fluid Flow Manual
Second Edition First Revision Second Edition January 1990 October 1992 March 1997
The information in this Manual has been jointly developed by Chevron Corporation and its Operating Companies. The Manual has been written to assist Chevron personnel in their work; as such, it may be interpreted and used as seen fit by operating management. Copyright 1990, 1992, 1997 CHEVRON CORPORATION. All rights reserved. This document contains proprietary information for use by Chevron Corporation, its subsidiaries, and affiliates. All other uses require written permission.
March 1997
Chevron Corporation
Section Front Matter Table of Contents 50 100 200 300 400 500 600 700 800 900 1000 1100 Appendix A Appendix B Appendix C Appendix D Appendix E Appendix F Appendix G Appendix H Appendix I Index PC Disks Page
Date March 1997 March 1997 March 1997 January 1990 January 1990 October 1992 March 1997 March 1997 January 1990 January 1990 January 1990 January 1990 March 1997 March 1997 January 1990 January 1990 January 1990 March 1997 October 1992 March 1997 January 1990 January 1990 March 1997 March 1997 March 1997
Chevron Corporation
March 1997
Requesting Signature Send this completed form to: Document Control, Room 50-4328 Chevron Research and Technology Company 100 Chevron Way (P.O. Box 1627) Richmond, CA 94802
Date
Chevron Corporation
March 1997
Name Address
Phone
Chevron Corporation
March 1997
Fluid Flow
Manual Sponsor: R.P. Hohmann / Phone: (510) 242-2216 / E-mail : rpho@chevron.com
This document contains extensive hyperlinks to figures and cross-referenced sections. The pointer will change to a pointing finger when positioned over text which contains a link.
List of Current Pages 50 100 200 300 400 500 600 700 800 900 1000 1100 Using This Manual Introduction Static Pressure Acceleration Pressure Drop Friction Pressure Drop Fitting Pressure Drop Noncircular Conduits Open Channel Flow (Section not developed) Surge Pressure Pipeline Flow Fluid Properties Computer Programs 800-1 900-1 1000-1 1100-1 50-1 100-1 200-1 300-1 400-1 500-1 600-1
Appendices Appendix A Appendix B Appendix C Appendix D Appendix E Appendix F Appendix G Appendix H Appendix I
Conversion Tables Properties of Water Design Properties of Pipe PCFLOW Program PIPEFLOW-2 Program HOTOIL Program HOTOL* Program SURGE Program PCSURGE Program
Chevron Corporation
June 1993
50
Chevron Corporation
50-1
March 1997
Organization
This manual comprises Engineering Guidelines and appendices that address such concerns as: (1) designing piping to efficiently carry fluids, (2) determining open channel flow, (3) calculating surge pressures, (4) handling special pipeline problems, (5) fluid and pipe properties, (6) available computer programs.
Tabs
The colored tabs in the manual will help you find information quickly. White tabs are for table of contents, introduction, appendices, PC disks, index, and general purpose topics Blue tabs denote Engineering Guidelines Red tab marks a place for you to keep documents that are developed at your facility
March 1997
50-2
Chevron Corporation
Chevron Corporation
50-3
March 1997
Fig. 50-1
Fluid Flow Manual Quick-Reference Guide Task Fluid Flow Manual Sections
Learning Background Information Pressure drop calculations Pipeline friction heating Surge Open channel flow Computer programs 100, 200, 300, 400, 500 900 800 700 1100, Appendices D, E, F, G, H, I
Selecting the Best Computer Program Selection guide Detailed operation 1100 Appendices D, E, F, G, H, I
Calculating Flow Rates By PCFLOW program With flow charts With sophisticated programs 1100, Appendix D 400 1100, Appendices D, E, F, G, H, I
Finding Engineering Data Pipe dimensions Fluid properties Heat transfer properties Appendix C 1000 900
March 1997
50-4
Chevron Corporation
100 Introduction
Abstract
This section describes the scope of the Fluid Flow Manual and discusses its basic approach to fluid flow problems. Contents 110 120 130 140 150 Scope of the Fluid Flow Manual Basic Elements of Pressure Drop Importance of the Darcy-Weisbach Equation Nomenclature References Page 100-2 100-2 100-2 100-3 100-3
Chevron Corporation
100-1
January 1990
100 Introduction
The relationship between the three components of pressure drop may be expressed as follows: Psystem = Pelevation + Pacceleration + Pfriction
(Eq. 100-1)
These components of total system pressure drop are treated in Sections 200, 300, and 400, respectively, for simple cases. Special considerations are treated in the remaining sections. For example, Section 500 presents a method for approximating the combination of both acceleration and friction losses that occurs in valves, fittings, and pipe entrances.
(Eq. 100-2)
January 1990
100-2
Chevron Corporation
100 Introduction
D = pipe internal diameter, ft V = fluid velocity, ft/sec g = gravitational constant (32.17 ft/sec2) The Darcy-Weisbach equation defines the friction factor, f. Whenever possible the reader is encouraged to use this equation instead of the flow charts in Section 400. This equation is automated in the Incompressible Flow section of the PCFLOW program, which is provided on disk at the end of this manual.
140 Nomenclature
This manual does not contain a master list of nomenclature. Equation variables are defined following each equation.
150 References
The following selection of general references is supplemented by specific references in the applicable sections of the manual. 1. 2. 3. 4. 5. Fox, R. W., A. T. McDonald. Introduction to Fluid Mechanics. John Wiley & Sons, New York: 1978. Perry, R. H., C. H. Chilton. Chemical Engineers Handbook, Section 5. McGraw-Hill, New York: 1973. Streeter, V. L., E. B. Wylie. Fluid Mechanics. McGraw-Hill, New York. Engineering Data Book, Section 17. Gas Processors Association, Tulsa: 1987. Cameron Hydraulic Data. Ingersoll-Rand, Woodcliff Lake, N.J.: 1979.
Chevron Corporation
100-3
January 1990
Chevron Corporation
200-1
January 1990
where: Ps = static pressure, psi h = height of liquid column, ft = fluid density, lbm/cu ft
January 1990
200-2
Chevron Corporation
Equation 200-2 expresses head, in feet, equivalent to an arbitrary pressure, in psi: 144 H = P -------
(Eq. 200-2)
where: H = head, ft P = pressure, psi = fluid density, lbm/cu ft The conversion of head in feet to pressure in pounds per square inch for water at 60F is as follows: Ps = 0.433 h h = 2.31 Ps
Chevron Corporation
200-3
January 1990
Chevron Corporation
300-1
October 1992
(Eq. 300-1)
where: P = static pressure drop, psi = fluid density, lbm/cu ft V1 = upstream fluid velocity, ft/sec V2 = downstream fluid velocity, ft/sec Determination of acceleration pressure drop is particularly important when calculating the NPSHA of reciprocating pumps, to avoid cavitation. See Section 100 of the Pump Manual.
330 Discussion
Equations 300-1 and 300-2 describe acceleration loss at pipe entrances and reducers. Frictional losses (see Section 500) must be added to get the total loss for
October 1992
300-2
Chevron Corporation
this geometry. The fitting loss coefficients given in Section 500 for other types of valves and fittings (besides pipe entrances and reducers) take into account both acceleration and friction effects. During changes of phase (evaporation, flashing, and boiling), the velocity of a fluid must increase as the gas phase increases its mass flow rate. The pressure required to produce that acceleration is accurately described by Equations 300-1 and 300-2. The total pressure drop is the sum of the acceleration pressure drop and the flowing friction pressure drop. This friction loss can be difficult to calculate because the flow rates of the two phases are changing and, therefore, the friction pressure drop is changing as the fluid moves downstream. The static pressure that is converted to kinetic energy through the acceleration of a flowing fluid is theoretically recoverable as static pressure when the flow decelerates. However, since even carefully designed diffusers can recover only a fraction of the original static pressure, this recovery is not attempted in normal piping situations. In standard piping systems the kinetic energy of a flowing fluid is dissipated as turbulence at pipe exits and enlargements. Confusion on this point can arise because some authors attribute acceleration pressure loss not to the pipe entrance or reducer, but to the pipe exit or enlargement, where the potentially recoverable energy is finally lost. This gives some readers the false impression that there is a static pressure drop across pipe exits and enlargements. Static pressure drop produced by acceleration and friction effectsoccurs across pipe entrances and reducers, not their exits and enlargements.
Chevron Corporation
300-3
October 1992
Chevron Corporation
400-1
March 1997
March 1997
400-2
Chevron Corporation
(Eq. 400-1)
where: P = pressure drop, psi W = mass flow rate, lbm/hr = fluid density, lbm/ft3
Chevron Corporation
400-3
March 1997
where: Re = Reynolds number = absolute viscosity, cp There are no sharp divisions between the laminar, transition, and turbulent flow regimes. For design purposes, the recommended boundary between laminar and transition flow is Re = 1600. The recommended boundary between transition and turbulent flow is Re = 3400. These values provide relatively smooth transitions between regimes for calculated friction factors, and produce conservative results (tend to overpredict pressure drop) around the laminar-to-transition flow boundary. The friction factor for laminar flow (Re < 1600) can be derived analytically (without experimental components) to give: f = 64/Re
(Eq. 400-5)
The friction factor for transition flow (1600 < Re < 3400) cannot be predicted accurately. The following conservative value (overprediction) is recommended for most cases: f = 0.04
(Eq. 400-6)
The Moody Diagram (Figure 400-1) presents experimentally derived friction factors for turbulent flow (Re 3400). In turbulent flow the friction factor is a function of pipe roughness as well as the Reynolds number. At high Reynolds numbers the friction factor is a function of only relative roughness (absolute roughness/diameter). Figure 400-2 gives the relative roughness for various diameters and types of pipe.
March 1997
400-4
Chevron Corporation
Fig. 400-1
Moody Diagram
Chevron Corporation
400-5
March 1997
Fig. 400-2
Relative Pipe Roughness (/D) and Friction Factors (f) for Complete Turbulence
March 1997
400-6
Chevron Corporation
Many equations have been proposed to approximate the Moody Diagram friction factors. One of these is the Chen Equation (Equation 400-7), which is simple, accurate, and stable when used on small computers:
2 2 f = ------------------------------------------ 4 log 10 (A1-A2)
(Eq. 400-7)
where: --D A1 = --------------3.7065 5.0452 A2 = --------------- log 10 ( A3 ) Re --- D 7.149 0.8981 A3 = ---------------------- + ------------ Re 2.8257 = absolute pipe roughness, ft
1.1098
414 Correction Factor for Internal Roughness for Use With Flow Charts
Figures 400-4 through 400-13 are based on new steel pipe having an absolute roughness of 0.0018 inches. The effect of other values of roughness can be estimated by multiplying the pressure drop by a correction factor from Figure 400-3. Typical values of roughness, E, are as follows: Pipe Plastic Smooth Steel, New Galvanized Steel Absolute Roughness, 0.000005 ft 0.00015 ft 0.00042 ft
Chevron Corporation
400-7
March 1997
Pipe Cast Iron, Asphalted Transite Cast Iron, Uncoated, New Steel, Concrete Lined Concrete Riveted Steel
Correction factors in Figure 400-3 for 1 centistoke (cs) are typical of water or petroleum products ranging from 0.5 to 2.0 cs viscosity. Correction factors in Figure 400-3 for 10 centistoke (cs) are typical of crude oils or other liquids in the viscosity range from 5 to 20 cs. Correction factors are applicable for turbulent flow. No correction is required for laminar flow. The uncertainties in the transition range increase with roughness. The correction factors are the ratio fc / fo where: fc = friction factor from the Colebrook formula (Equation 400-12) fo = friction factor on which Figures 400-4 to 400-13 are based Where accurate performance data are required, pressure losses should be determined by test. If test measurements are not possible, the friction factor can be found with the Moody Diagram or calculated with the Chen Equation (Equation 400-7).
Fig. 400-3 Correction Factors for Internal Roughness (1 of 2)
Use for Viscosity = 1 Centistoke and Turbulent flow. No correction for laminar flow. Absolute Roughness, in. ID, in. Velocity, ft/sec 0.0010 0.0018 0.0050 0.0100 0.0300 0.1000 0.3000
3 5 10
3 5 10
10
3 5 10
March 1997
400-8
Chevron Corporation
Fig. 400-3
ID, in. 20
Velocity, ft/sec 3 5 10
50
3 5 10
Use for Viscosity = 10 centistoke and turbulent flow. Absolute Roughness, in. ID, in. Velocity, ft/sec 0.0010 0.0018 0.0050 0.0100 0.0300 0.1000 0.3000
3 5 10
3 5 10
10
3 5 10
20
3 5 10
50
3 5 10
Chevron Corporation
400-9
March 1997
Fig. 400-4
1-Inch PipeSchedule 40 (1 of 2) Correction Factor Table for 1 in. Pipe of Various Thicknesses
Schedule 5S 10S
Transition & Turbulent Flow 0.56 0.81 Fp 1.39 1.13 Fq 1.00 0.78 0.51 0.61 0.84 Fp 1.00 1.44 2.74
40 80 160
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 40. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 40. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity.
EXAMPLE 1: Given: Determine: Solution: Flow rate = 5 BPH; viscosity = 20 cs; specific gravity = 0.9; line size = 1 in. schedule 10S (ID = 1.097 in.) Pressure loss (psi/1000 ft) Enter flow chart at 5 BPH. Move across to viscosity of 20 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 14.3 psi/1000 ft. for 1.049-in. ID pipe. From correction table, for 1.097-in. ID pipe and laminar flow, Fp = 0.84 Pressure loss = (14.3) (0.84) = 12 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 16.0 psi/1000 ft; viscosity = 2 cs; specific gravity = 0.9; line size = 1 in. schedule 80 (ID = 0.957 in.) Flow rate (BPH) Enter flow chart at 16 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 2.0 cs in turbulent range. Move horizontally to flow rate of 10 BPH for 1.049-in. ID pipe. From correction table, for 0.957-in. ID pipe and turbulent flow, Fq = 0.78. Flow rate = (10) (0.78) = 7.8 BPH
March 1997
400-10
Chevron Corporation
Fig. 400-4
1-Inch PipeSchedule 40 (2 of 2)
Chevron Corporation
400-11
March 1997
Fig. 400-5
1-1/2-Inch PipeSchedule 40 (1 of 2) Correction Factor Table for 1-1/2 in. Pipe of Various Thicknesses
Schedule
Transition & Turbulent Flow Fp Fq 1.29 1.13 1.00 0.83 0.61 Fp 0.68 0.84 1.00 1.33 2.10
5S 10S 40 80 160
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 40. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 40. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate 12 BPH; viscosity = 20 cs; specific gravity = 0.9; line size = 1-1/2 in. schedule 10S (ID = 1.682 in.) Pressure loss (psi/1000 ft) Enter flow chart at 12 BPH. Move across to viscosity of 20 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 6.1 psi/1000 ft. for 1.610-in. ID pipe. From correction table, for 1.682-in. ID pipe and laminar flow, Fp = 0.84 Pressure loss = (6.1) (0.84) = 5.1 psi/1000 feet
EXAMPLE 2: Given: Determine: Solution: Pressure loss = 15.7 psi/1000 ft; viscosity = 5 cs; specific gravity = 0.9; line size = 1-1/2 in. schedule 80 (ID = 1.500 in.) Flow rate (BPH) Enter flow chart at 15.7 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 5.0 cs in turbulent range. Move horizontally to flow rate of 28 BPH for 1.610-in. ID pipe. From correction table, for 1.500-in. ID pipe and turbulent flow, Fq = 0.83. Flow rate = (28) (0.83) = 23.2 BPH
March 1997
400-12
Chevron Corporation
Fig. 400-5
1-1/2-Inch PipeSchedule 40 (2 of 2)
Chevron Corporation
400-13
March 1997
Fig. 400-6
2-Inch PipeSchedule 40 (1 of 2) Correction Factor Table for 2 in. Pipe of Various Thicknesses
Schedule
Transition & Turbulent Flow Fp Fq 1.25 1.12 1.00 0.84 0.58 Fp 0.72 0.84 1.00 1.29 2.24
5S 10S 40 80 160
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 40. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 40. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate = 24 BPH; viscosity = 60 cs; specific gravity = 0.9; line size = 2 in. schedule 10S (ID = 2.157 in.) Pressure loss (psi/1000 ft) Enter flow chart at 24 BPH. Move across to viscosity of 60 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 13.8 psi/1000 ft. for 2.067-in. ID pipe. From correction table, for 2.157-in. ID pipe and laminar flow, Fp = 0.84 Pressure loss = (13.8) (0.84) = 11.6 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 12.4 psi/1000 ft; viscosity = 5 cs; specific gravity = 0.9; line size = 2 in. Schedule 80 (ID = 1.939 in.) Flow rate (BPH) Enter flow chart at 12.4 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 5.0 cs in turbulent range. Move horizontally to flow rate of 48 BPH for 2.067-in. ID pipe. From correction table, for 1.939-in. ID pipe and turbulent flow, Fq = 0.84. Flow rate = (48) (0.84) = 40.3 BPH
March 1997
400-14
Chevron Corporation
Fig. 400-6
2-Inch PipeSchedule 40 (2 of 2)
Chevron Corporation
400-15
March 1997
Fig. 400-7
2-1/2-Inch PipeSchedule 40 (1 of 2) Correction Factor Table for 2-1/2 in. Pipe of Various Thicknesses
Schedule
Transition & Turbulent Flow Fp Fq 1.28 1.19 1.00 0.85 0.67 Fp 0.69 0.77 1.00 1.29 1.82
5S 10S 40 80 160
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 40. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 40. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate = 40 BPH; viscosity = 60 cs; specific gravity = 0.9; line size = 2-1/2 in. schedule 10S (ID = 2.635 in.) Pressure loss (psi/1000 ft) Enter flow chart at 40 BPH. Move across to viscosity of 60 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 11.2 psi/1000 ft. for 2.469-in. ID pipe. From correction table, for 2.635-in. ID pipe and laminar flow, Fp = 0.77 Pressure loss = (11.2) (0.77) = 8.6 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 12.8 psi/1000 ft; viscosity = 10 cs; specific gravity = 0.9; line size = 2-1/2 in. schedule 80 (ID = 2.323 in.) Flow rate (BPH) Enter flow chart at 12.8 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 10 cs in turbulent range. Move horizontally to flow rate of 74 BPH for 2.469-in. ID pipe. From correction table, for 2.323-in. ID pipe and turbulent flow, Fq = 0.85. Flow rate = (74) (0.85) = 63 BPH
March 1997
400-16
Chevron Corporation
Fig. 400-7
2-1/2-Inch PipeSchedule 40 (2 of 2)
Chevron Corporation
400-17
March 1997
Fig. 400-8
3-Inch PipeSchedule 40 (1 of 2) Correction Factor Table for 3 in. Pipe of Various Thicknesses
Schedule
Transition & Turbulent Flow Fp Fq 1.25 1.18 1.00 0.86 0.66 Fp 0.72 0.78 1.00 1.25 1.87
5S 10S 40 80 160
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 40. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 40. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate = 65 BPH; viscosity = 70 cs; specific gravity = 0.9; line size = 3 in. schedule 10S (ID = 3.260 in.) Pressure loss (psi/1000 ft) Enter flow chart at 65 BPH. Move across to viscosity of 70 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 8.8 psi/1000 ft. for 3.068-in. ID pipe. From correction table, for 3.260-in. ID pipe and laminar flow, Fp = 0.78 Pressure loss = (8.8) (0.78) = 6.9 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 10.7 psi/1000 ft; viscosity = 10 cs; specific gravity = 0.9; line size = 3 in. schedule 80 (ID = 2.900 in.) Flow rate (BPH) Enter flow chart at 10.7 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 10 cs in turbulent range. Move horizontally to flow rate of 120 BPH for 3.068-in. ID pipe. From correction table, for 2.900-in. ID pipe and turbulent flow, Fq = 0.86. Flow rate = (120) (0.86) = 103 BPH
March 1997
400-18
Chevron Corporation
Fig. 400-8
3-Inch PipeSchedule 40 (2 of 2)
Chevron Corporation
400-19
March 1997
Fig. 400-9
4-Inch PipeSchedule 40 (1 of 2) Correction Factor Table for 4 in. Pipe of Various Thicknesses
Schedule
Transition & Turbulent Flow Fp Fq 1.22 1.16 1.00 0.87 0.75 0.65 Fp 0.74 0.80 1.00 1.23 1.52 1.88
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 40. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 40. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate = 96 BPH; viscosity = 80 cs; specific gravity = 0.9; line size = 4 in. schedule 10S (ID = 4.260 in.) Pressure loss (psi/1000 ft) Enter flow chart at 96 BPH. Move across to viscosity of 80 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 5.05 psi/1000 ft. for 4.026-in. ID pipe. From correction table, for 4.260-in. ID pipe and laminar flow, Fp = 0.80 Pressure loss = (5.05) (0.80) = 4.04 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 7.9 psi/1000 ft; viscosity = 10 cs; specific gravity = 0.9; line size = 4 in. schedule 80 (ID = 3.826 in.) Flow rate (BPH) Enter flow chart at 7.9 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 10 cs in turbulent range. Move horizontally to flow rate of 210 BPH for 4.026-in. ID pipe. From correction table, for 3.826-in. ID pipe and turbulent flow, Fq = 0.87.
March 1997
400-20
Chevron Corporation
Fig. 400-9
4-Inch PipeSchedule 40 (2 of 2)
Chevron Corporation
400-21
March 1997
Fig. 400-10 6-Inch Pipe1/4-Inch Wall (1 of 2) Correction Factor Table for 6 in. Pipe of Various Thicknesses Schedule Inside Diameter, in. Transition & Turbulent Flow Fp 5S 10S 1/4 in. wall 40 80 120 160 6.407 6.357 6.125 6.065 5.761 5.501 5.189 .81 .84 1.00 1.05 1.34 1.67 2.21 Fq 1.13 1.11 1.00 .97 .85 .75 .64 Fp .84 .86 1.00 1.04 1.28 1.54 1.94 Laminar Flow Fq 1.20 1.16 1.00 .96 .78 .65 .52
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than 1/4-inch. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than 1/4-inch. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate = 400 BPH; viscosity = 150 cs; specific gravity = 0.9; line size = 6 in. schedule 40 (ID = 6.065 in.) Pressure loss (psi/1000 ft) Enter flow chart at 400 BPH. Move across to viscosity of 150 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 7.3 psi/1000 ft. for 6.125-in. ID pipe. From correction table, for 6.065-in. ID pipe and laminar flow, Fp = 1.04 Pressure loss = (7.3) (1.04) = 7.6 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 5.9 psi/1000 ft; viscosity = 20 cs; specific gravity = 0.9; line size = 6 in. schedule 80 (ID = 5.761 in.) Flow rate (BPH) Enter flow chart at 5.9 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 20 cs in turbulent range. Move horizontally to flow rate of 520 BPH for 6.125-in. ID pipe. From correction table, for 5.761-in. ID pipe and turbulent flow, Fq = 0.85.
March 1997
400-22
Chevron Corporation
Chevron Corporation
400-23
March 1997
Fig. 400-11 8-Inch Pipe1/4-Inch Wall (1 of 2) Correction Factor Table for 8 in. Pipe of Various Thicknesses Schedule 10S 20 40 60 80 120 160 Inside Diameter, in. 8.329 8.125 7.981 7.813 7.625 7.189 6.813 Transition & Turbulent Flow Fp 0.89 1.00 1.09 1.21 1.35 1.79 2.32 Fq 1.07 1.00 0.95 0.90 0.84 0.72 0.62 Fp 0.91 1.00 1.07 1.17 1.29 1.63 2.02 Laminar Flow Fq 1.10 1.00 0.93 0.86 0.78 0.61 0.49
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 20. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 20. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate = 840 BPH; viscosity = 250 cs; specific gravity = 0.9; line size = 8 in. schedule 40 (ID = 7.981 in.) Pressure loss (psi/1000 ft) Enter flow chart at 840 BPH. Move across to viscosity of 250 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 8.3 psi/1000 ft. for 8.125-in. ID pipe. From correction table, for 7.981-in. ID pipe and laminar flow, Fp = 1.07 Pressure loss = (8.3) (1.07) = 8.9 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 4.15 psi/1000 ft; viscosity = 20 cs; specific gravity = 0.9; line size = 8 in. schedule 80 (ID = 7.625 in.) Flow rate (BPH) Enter flow chart at 4.15 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 20 cs in turbulent range. Move horizontally to flow rate of 910 BPH for 8.125-in. ID pipe. From correction table, for 7.625-in. ID pipe and turbulent flow, Fq = 0.84.
March 1997
400-24
Chevron Corporation
Chevron Corporation
400-25
March 1997
Fig. 400-12 10-Inch Pipe1/4-Inch Wall (1 of 2) Correction Factor Table for 10 in. Pipe of Various Thicknesses Schedule Inside Diameter, in. Transition & Turbulent Flow Fp 10S 20S 40 60 80 120 160 10.420 10.250 10.020 9.750 9.564 9.064 8.500 0.92 1.00 1.11 1.27 1.39 1.80 2.44 Fq 1.05 1.00 0.94 0.87 0.83 0.72 0.60 Fp 0.94 1.00 1.10 1.22 1.32 1.64 2.11 Laminar Flow Fq 1.07 1.00 0.91 0.82 0.76 0.61 0.47
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 20S. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 20S. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate = 1100 BPH; viscosity = 300 cs; specific gravity = 0.9; line size = 10 in. schedule 40 (ID = 10.020 in.) Pressure loss (psi/1000 ft) Enter flow chart at 1100 BPH. Move across to viscosity of 300 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 5.15 psi/1000 ft. for 10.250-in. ID pipe. From correction table, for 10.020-in. ID pipe and laminar flow, Fp = 1.10 Pressure loss = (5.15) (1.10) = 5.7 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 3.70 psi/1000 ft; viscosity = 20 cs; specific gravity = 0.9; line size = 10 in. schedule 80 (ID = 9.564 in.) Flow rate (BPH) Enter flow chart at 3.70 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 20 cs in turbulent range. Move horizontally to flow rate of 1600 BPH for 10.250-in. ID pipe. From correction table, for 9.564-in. ID pipe and turbulent flow, Fq = 0.83.
March 1997
400-26
Chevron Corporation
Chevron Corporation
400-27
March 1997
Fig. 400-13 12-Inch Pipe1/4-Inch Wall (1 of 2) Correction Factor Table for 12 in. Pipe of Various Thicknesses Schedule Inside Diameter, in. Transition & Turbulent Flow Fp 10S 20S 40 80 120 160 12.390 12.250 11.938 11.376 10.750 10.126 0.95 1.00 1.13 1.42 1.86 2.48 Fq 1.03 1.00 0.93 0.82 0.70 0.60 Fp 0.96 1.00 1.11 1.34 1.69 2.14 Laminar Flow Fq 1.05 1.00 0.90 0.74 0.59 0.47
Fp = Pressure loss correction factor Fq = Flow rate correction factor NOTES: 1. Multiply pressure loss from flow chart by Fp for pressure loss with pipe walls other than Schedule 20S. 2. Multiply flow rate from flow chart by Fq to obtain flow rate with pipe walls other than Schedule 20S. 3. For SG1.0, multiply pressure loss at SG 1.0 by actual SG to obtain pressure loss. For known pressure loss, divide by SG, then enter chart at SG 1.0 to determine flow rate. 4. In laminar range, pressure loss is directly proportional to viscosity. To determine pressure losses for viscosities not shown, the ratio of a known viscosity to pressure loss at desired flow rate is applied to the actual viscosity. EXAMPLE 1: Given: Determine: Solution: Flow rate = 1500 BPH; viscosity = 400 cs; specific gravity = 0.9; line size = 12 in. schedule 40 (ID = 11.938 in.) Pressure loss (psi/1000 ft) Enter flow chart at 1500 BPH. Move across to viscosity of 400 cs. Move vertically to SG 0.9. Move diagonally to pressure loss of 4.55 psi/1000 ft. for 12.250-in. ID pipe. From correction table, for 11.938-in. ID pipe and laminar flow, Fp = 1.11 Pressure loss = (4.55) (1.11) = 5.05 psi/1000 feet EXAMPLE 2: Given: Determine: Solution: Pressure loss = 3.49 psi/1000 ft; viscosity = 20 cs; specific gravity = 0.9; line size = 12 in. schedule 80 (ID = 11.376 in.) Flow rate (BPH) Enter flow chart at 3.49 psi/1000 feet. Move diagonally to SG 0.9. Move vertically to viscosity of 20 cs in turbulent range. Move horizontally to flow rate of 2500 BPH for 12.250-in. ID pipe. From correction table, for 11.376-in. ID pipe and turbulent flow, Fq = 0.82.
March 1997
400-28
Chevron Corporation
Chevron Corporation
400-29
March 1997
Limitations
The method described here applies to isothermal gas-liquid flow, not to situations in which a phase change occurs; that is, constant gas-liquid ratios (by weight) are assumed. This method has not been verified for very long vertical piping (such as in oil wells) nor has the accuracy been established for horizontal piping more than 5-1/2 inches in diameter. In these cases the method should be used with caution, for vertical piping, PIPEFLOW-2 will yield better results. In addition, the limited experimental data available indicate that when the mixture velocity is less than 3 ft/sec the accuracy of the friction pressure drop calculations is very poor. This method is not fully applicable to flow of water-oil-gas (WOG) mixtures (socalled three-phase flow). This case requires the more powerful calculation methods of PIPEFLOW-2.
General References
Reference 1 (see Section 450) contains a more detailed discussion of two-phase flow. Reference 2 contains an extensive bibliography of two-phase literature.
The components of this equation, Pfriction, Pfittings, Pacceleration, and Pelevation are discussed in the following sections. The total pressure drop is calculated by evaluating each component individually and summing.
March 1997
400-30
Chevron Corporation
The five most widely used correlations are compared in Reference 3 using experimental data from a number of investigators. The data were carefully screened to eliminate unreliable measurements. The screened data, about 2600 points in all, cover pipe diameters from 1 to 5-1/2 inches and liquid viscosities from 1 to 20 centipoise. Of the five the most reliable correlation over this range of experimental conditions was the Lockhart-Martinelli correlation (see Reference 4). Another somewhat better correlation with the screened experimental data was achieved using similarity analysis (see Reference 5). This method is based on calculating a two-phase density, tp, and viscosity, tp, evaluated at the pipe entrance pressure and temperature and assumed constant for the friction and fitting pressure drop calculation, as follows: tp = l () + g (1.0 - )
(Eq. 400-9)
tp = l () + g (1.0 - )
(Eq. 400-10)
= liquid volume fraction at pipe entrance Equations 400-9 and 400-10 assume that both phases flow at the same velocity. The two-phase Reynolds number Retp is expressed as follows: 0.527W t V m D tp Re tp = -------------------- = -------------------D tp tp ----------- 1490
(Eq. 400-11)
where: Vm = velocity of mixture, ft/sec D = pipe inside diameter, ft Wt = mass flow rate of total fluid, lbm/hr For new steel pipe, two-phase Reynolds numbers should be used with the Moody diagram (Figure 400-1) to determine the friction factor f. If different pipe conditions exist or a more accurate determination is desired, the Colebrook formula (Equation 400-12) may be used.
Chevron Corporation
400-31
March 1997
where: = absolute pipe wall roughness, ft The need to proceed by trial and error is an inconvenience when using this equation for hand calculation, but a computer or Moody chart eliminates this problem. The equation reduces to the smooth tube equation when the wall roughness (left term in bracket) approaches zero or to Nikuradses Formula at high Reynolds numbers (when the right term in bracket approaches zero). The same absolute wall roughness, , should be used for both single-phase and two-phase flow calculations. The pressure drop due to friction may then be calculated as follows: L tp V m - P friction = f --- -------- ---------D 144 2g o L Wt f ------ --------5 D tp
(Eq. 400-13)
2 2
= 1.35 10
11
where: P = pressure drop, psi f = friction factor L = pipe length, ft go = gravitational constant (32.174 lbm ft/lbf sec2) Wl = flow rate of liquid, lbm/hr Wg = flow rate of gas, lbm/hr W t = Wl + W g This method of calculating friction pressure drop has the following characteristics: It reduces to the single-phase flow equations if the flow rate of either phase is zero. Except for the assumptions concerning two-phase density and viscosity (Equations 400-9 and 400-10), no empirical factors from two-phase flow data have been used. It is reasonably accurate for all flow patterns (see Section 427).
March 1997
400-32
Chevron Corporation
tp
(Eq. 400-14)
where: Z = compressibility factor T = temperature, R R = gas constant P1 = upstream pressure, psi P2 = downstream pressure, psi D = pipe inside diameter, ft
Chevron Corporation
400-33
March 1997
The actual flow density depends on how the liquid and gas are distributed in the pipe. The flow density in a short section of pipe of length L is given by Equation 400-16: LA g g + LA 1 1 tp = ---------------------------------------- = R g g + R 1 1 LA
(Eq. 400-16)
where: = actual density in pipe section Ag = area gas Al = area liquid Rg = gas volume fraction Rl = liquid holdup Rg is the fractional volume of the pipe filled with gas and Rl is the fractional volume of the pipe filled with liquid. Rl is called liquid holdup (see Equation 400-22). Because of the difference in velocity of the two phases, liquid holdup is greater downstream than at the entrance. Therefore, to calculate the actual flow density, the liquid holdup Rl has to be known along the pipe. The available correlations for liquid holdup were checked against experimental data from Reference 3. The correlation developed by Hughmark (see Reference 7, and below in this section, Liquid Holdup Correlation) was the best. The effects of bends and fittings on liquid holdup and, therefore, the flow density cannot be predicted at this time. Therefore, it is assumed that the same holdup correlation can be used even if the pipe contains bends and fittings. In two-phase flow, as in single-phase flow, the elevation head loss is expressed as follows: tp P elevation = -------- h 144
(Eq. 400-17)
where: h = static elevation, ft The flow density is calculated using Equation 400-16, where the gas density is evaluated at the average pressure. The elevation pressure drop term is included only in vertical upward flow. A conservative evaluation of acceleration pressure loss for vertical downward flow cannot take credit for the elevation pressure component in the downward section. Therefore, sections where the flow is downward should be treated as horizontal piping. No provisions have been made to handle inclined piping.
March 1997
400-34
Chevron Corporation
Figure 400-17 can be used to estimate the accuracy of a calculated friction pressure drop for any flow regime. For example, the calculated friction pressure drop for horizontal slug flow is within -18.0 to +12.0 percent of the actual value. Equation 400-18 may be restated as follows: P calc P exp = -----------------------% dev 1 + -------------100
(Eq. 400-19)
Chevron Corporation
400-35
March 1997
March 1997
400-36
Chevron Corporation
Fig. 400-15 Flow Pattern Map for Vertical Two-Phase Flow From Two Phase Slug Flow by Griffith & Wallis. Journal of Heat Transfer, Transactions of ASME Series C83 (Aug., 1961). Courtesy of ASME
Chevron Corporation
400-37
March 1997
Fig. 400-16 Calculated vs. Experimental Frictional Pressure DropHorizontal Flow Flow Regime Plug Stratified Wave Slug Annular Dispersed Bubble Range of Deviation (%) -22.3 to -2.3 -25.3 to +24.7 -21.0 to +39.0 -17.9 to +12.1 -59.2 to +15.8 -24.4 to +30.6 not given
Based on the range of deviation for horizontal slug flow, the actual value of a calculated pressure drop of 10 psi would be (approximately) between the following values: 10 P = --------------------------- = 12.2 psi 1 + ( 0.18 )
(Eq. 400-20)
A comparison between calculated and experimental friction pressure drop for vertical flow is not available.
March 1997
400-38
Chevron Corporation
Fig. 400-17 Correlation for the Flow Pattern Y From Holdup in Gas-Liquid Flow by G.A. Hughmark, Chemical Engineering Progress, Vol. 58, April, 1962, p. 62
Re Fr X = ---------------------1 where: Fr = Froude number = V2/Dg = liquid volume fraction at pipe entrance
-4
(Eq. 400-23)
Chevron Corporation
400-39
March 1997
D = diameter, ft g = gravitational constant (32.174 ft/sec2) The dimensionless numbers used in the variable X are shown in Equations 400-24 and 400-25. D Gm -----------------------------------( R11 + Rg g ) Re = ------------------------------------1490
(Eq. 400-24)
where: Gm = tpVm = mass velocity mixture (lbm/ft2-sec) Re Retp (from Equation 400-11) Vm ( ( Q1 + Q g ) A ) F r = ---------- = ---------------------------------------gD gD
2 2
(Eq. 400-25)
Q1 W1 1 = ----------------------------------- = ------------------W1 1 + Wg g Q 1 + Qg
(Eq. 400-26)
where: = specific volume The calculation procedure is to evaluate Re, Fr, and using Equations 400-24, 400-25, and 400-26. The variable X is then evaluated using Equation 400-23, and the flow parameter Y is determined from Figure 400-16. Using the flow parameter Y, the liquid holdup is found from Equation 400-22. An iterative calculation is required since the gas density used in Equation 400-22 is evaluated at the average pressure. The gas volume flow rate Qg used in Equations 400-25 and 400-26 is the inlet value evaluated using the inlet density. The actual flow density calculated using Equation 400-16 is then used to determine the elevation pressure drop in upwards vertical flow. The deviation between the calculated (Figure 400-17) and experimental (Reference 3) values of the liquid holdup varies by 25%. For vertical flow not as much experimental data are available. For the available data the deviation between experimental and calculated liquid holdup does not exceed 10 percent (see Section 450, Reference 7).
March 1997
400-40
Chevron Corporation
Of these options only COMFLOW and PIPEFLOW-2 consider heat transfer, and only PIPEFLOW-2 considers condensation. Condensation due to heat transfer is common in hot gas transmission and can significantly affect the friction pressure drop. Section 420 discusses two-phase flow pressure drop.
(Eq. 400-27)
where: Q = flow rate, SCFD To = standard absolute temperature, R Po = standard pressure, psia D = pipe ID, in. P1 = upstream pressure, psia P2 = downstream pressure, psia S = fluid specific gravity (air = 1) T = fluid absolute temperature, R L = length of pipeline, miles C1 through C7 = constants as shown in Figure 400-26
Chevron Corporation
400-41
March 1997
Fig. 400-18 Flow of Any Fluid1/2 to 6-Inch Pipe (1 of 2) Flow of Any Fluid (Gas, Vapor, Liquid) Nomenclature: M Mass flow rate, lb/hr P P1 & P2 P T Average line pressure, psia Initial and final pressure, psia Rate of pressure drop, psi/1000 ft Absolute temperature, F + 460 Viscosity, cp Specific volume of fluid, cu ft/sec Specific gravity of gas referred to dry air at standard conditions Line length, thousands of feet 2. Fluid, carbon dioxide gas; temperature = 70F; viscosity = 0.015 cp; inlet pressure = 25 psig; outlet pressure = 20 psig; line length = 600 ft; line size = 2 in. sch 40 = 2.067 in. ID. Average flowing pressure = 22.5 psig.
V G L
Examples: 1. Fluid, carbon dioxide gas; flow rate = 2500 lb/hr; temperature = 70F; viscosity = 0.015 cp; inlet pressure 15 psig; pipe size = 3 in. sch 40 = 3.068 in. ID. Assume average line pressure = 12 psig. Determine pressure loss. Miscellaneous Data: Specific volume of a perfect gas, V = 10.72 T / P(MOL. WT) 1 cubic foot per minute of gas at standard conditions = 4.58G pounds per hour 1 inch of water = 0.0361 pounds per square inch Notes: 1. The chart, strictly speaking, gives rate of pressure drop at a point in the pipe, but for a perfect gas will give the average rate of pressure drop if the specific volume at the average pressure is used. 2. Standard conditions = 60F and 14.73 psia
Chevron Corporation
400-43
March 1997
Fig. 400-19 Flow of Any Fluid2 to 48-Inch Pipe (1 of 2) Flow of Any Fluid (Gas, Vapor, Liquid) Nomenclature: M Mass flow rate, lb/hr P P1 & P2 P T Average line pressure, psia Initial and final pressure, psia Rate of pressure drop, psi/1000 ft Absolute temperature, F + 460 Viscosity, cp Specific volume of fluid, cu ft/sec Specific gravity of gas referred to dry air at standard conditions Line length, thousands of feet 2. Fluid, carbon dioxide gas; temperature = 70F; viscosity = 0.015 cp; inlet pressure = 25 psig; outlet pressure = 20 psig; line length = 600 ft; line size = 2 in. sch 40 = 2.067 in. ID. Average flowing pressure = 22.5 psig.
V G L
Examples: 1. Fluid, carbon dioxide gas; flow rate = 50 lb/hr; temperature = 70F; viscosity = 0.015 cp; inlet pressure 15 psig; pipe size = 3/4 in. sch 40 = 0.824 in. ID. Assume average line pressure = 12 psig. Determine pressure loss. Miscellaneous Data: Specific volume of a perfect gas, V = 10.72 T / P(MOL. WT) 1 cubic foot per minute of gas at standard conditions = 4.58G pounds per hour 1 inch of water = 0.0361 pounds per square inch Notes: 1. The chart, strictly speaking, gives rate of pressure drop at a point in the pipe, but for a perfect gas will give the average rate of pressure drop if the specific volume at the average pressure is used.
Chevron Corporation
400-45
March 1997
Fig. 400-20 Flow of Natural Gas1/2 to 6-Inch Pipe (1 of 2) Flow of Natural Gas Nomenclature: Q Flow rate, millions of standard cubic feet per day P P1 & P2 P L T D G Average line pressure, psia Initial and final pressure, psia Rate of pressure drop, psi/1000 ft Length of line, thousands of feet Temperature absolute degrees, F + 460 Internal diameter of pipe, in. Specific gravity of gas referred to dry air at standard conditions Viscosity, cp Supercompressibility factor 2. Line size = 4 in. sch 40 = 4.026 in. ID; inlet pressure = 100 psia; outlet pressure = 80 psia; specific gravity = 0.70.; viscosity = 0.012 cp; line length = 10,000 ft; temperature = 60F. Assume a perfect gas. Find flow rate.
Examples: 1. Flow quantity = 1.3 million cu ft/day; specific gravity = 0.70; viscosity = 0.012 cp; temperature = 200F; inlet pressure 2050 psia; pipe size = 2 in. sch 40; line length = 3000 ft. Find pressure drop. From GPSA (see Section 450, reference 14), supercompressibility factor = 0.882. Assume average pressure = 2000 psia.
Notes: 1. The chart, strictly speaking, gives rate of pressure drop at a point in the pipe, but will give the average rate of pressure drop if the specific volume at the average pressure is used. 2. Standard conditions = 60F and 14.73 psia
Chevron Corporation
400-47
March 1997
Fig. 400-21 Flow of Natural Gas2 to 48-Inch Pipe (1 of 2) Flow of Natural Gas Nomenclature: Q Flow rate, millions of standard cubic feet per day P P1 & P2 P L T D G Average line pressure, psia Initial and final pressure, psia Rate of pressure drop, psi/1000 ft Length of line, thousands of feet Temperature absolute degrees, F + 460 Internal diameter of pipe, in. Specific gravity of gas referred to dry air at standard conditions Viscosity, cp Supercompressibility factor 2. Line size = 4 in. sch 40 = 4.026 in. ID; inlet pressure = 100 psia; outlet pressure = 80 psia; specific gravity = 0.70; viscosity = 0.012 cp; line length = 10,000 ft; temperature = 60F. Assume a perfect gas. Find flow rate.
Examples: 1. Flow quantity = 1.3 million cu ft/day; specific gravity = 0.70; viscosity = 0.012 cp; temperature = 200F; inlet pressure 2050 psia; pipe size = 2 in. sch 40; line length = 3000 ft. Find pressure drop. From GPSA (see Section 450, reference 14), supercompressibility factor = 0.882. Assume average pressure = 2000 psia.
Notes: 1. The chart, strictly speaking, gives rate of pressure drop at a point in the pipe, but will give the average rate of pressure drop if the specific volume at the average pressure is used. 2. Standard conditions = 60F and 14.73 psia.
Chevron Corporation
400-49
March 1997
Fig. 400-22 Flow of Air1/2 to 6-Inch Pipe (1 of 2) Flow of Air Nomenclature: Q Flow rate, standard cubic feet per minute P P1 & P2 P L T D Average line pressure, psia Initial and final pressure, psia Rate of pressure drop, psi/1000 ft Length of line, thousands of feet Absolute temperature, F + 460 Internal diameter, in. Miscellaneous Data: 1 inch of water = 0.0361 psi 1 lb per hour = 13.2 cu ft of air per hour at standard conditions Notes: 1. The chart, strictly speaking, gives rate of pressure drop at a point in the pipe, but will give the average rate of pressure drop if the specific volume at the average pressure is used. 2. Standard conditions = 60F and 14.73 psia 2. Line size = 4 in. sch 40 = 4.026 in. ID. Length of line = 300 ft, inlet pressure = 100 psia, outlet pressure = 90 psia, temperature = 60 F.
Examples: 1. Line size = 1 in. sch 40 = 1.049 in. ID. Flow rate = 78 scfm, inlet pressure = 100 psia, temperature = 60F. Line length = 400 ft. Assume average pressure = 95 psia. From the chart.
Chevron Corporation
400-51
March 1997
Fig. 400-23 Flow of Air2 to 48-Inch Pipe (1 of 2) Flow of Air Nomenclature: Q Flow rate, standard cubic feet per minute P P1 & P2 P L T D Average line pressure, psia Initial and final pressure, psia Rate of pressure drop, psi/1000 ft Length of line, thousands of feet Absolute temperature, F + 460 Internal Diameter, in. Miscellaneous Data: 1 inch of water = 0.0361 psi 1 lb per hour = 13.2 cu ft of air per hour at standard conditions Notes: 1. The chart, strictly speaking, gives rate of pressure drop at a point in the pipe, but will give the average rate of pressure drop if the specific volume at the average pressure is used. 2. Standard conditions = 60F and 14.73 psia 2. Line size = 4 in. sch 40 = 4.026 in. ID. Length of line = 300 ft, inlet pressure = 100 psia, outlet pressure = 90 psia, temperature = 60 F.
Examples: 1. Line size = 4 in. sch 40 = 4.026 in. ID. Flow rate = 2650 scfm, inlet pressure = 100 psia, temperature = 60F. Line length = 400 ft. Assume average pressure = 95 psia. From the chart.
Chevron Corporation
400-53
March 1997
Fig. 400-24 Flow of Steam1/2 to 6-Inch Pipe (1 of 2) Table of Specific Volumes For Steam
Specific Volume (V) ABS PRESS. 1 in. Hg 1-1/2 2 2-1/2 3 2 psia 3 14.7 20 30 40 50 60 70 80 90 100 110 120 130 140 150 200 250 300 400 500 600 700 800 900 1000 1500 2000 SAT. TEMP, F 79.1 91.8 101.2 108.7 115.1 126.1 141.5 212.0 228.0 250.3 267.2 281.0 292.7 302.9 312.0 320.3 327.8 334.8 341.3 347.3 353.0 358.4 381.8 401.0 417.3 444.6 467.0 486.2 503.0 518.2 532.0 544.6 596.1 635.6 SAT. STEAM 653 445 340 275 232 174.0 118.9 26.8 20.1 13.74 10.50 8.51 7.17 6.20 5.47 4.89 4.43 4.04 3.72 3.45 3.22 3.01 2.28 1.841 1.541 1.160 .926 .768 .653 .565 .497 .442 .274 .1875 300F 922 615 461 369 308 226 150.5 30.5 22.4 14.82 11.04 8.78 7.26 400F 1142 695 521 417 348 256 170.4 34.6 25.4 16.89 12.62 10.06 8.35 7.13 6.22 5.50 4.93 4.47 4.08 3.75 3.46 3.22 2.36 500F 1163 776 582 466 388 286 190.3 38.8 28.4 18.92 14.16 11.30 9.40 8.04 7.02 6.22 5.58 5.06 4.63 4.26 3.95 3.68 2.72 2.15 1.765 1.283 .990 .792 600F 1287 857 844 515 430 315 210 42.8 31.5 20.9 15.68 12.53 10.42 8.92 7.79 6.92 6.22 5.64 5.16 4.76 4.41 4.11 3.06 2.42 2.00 1.474 1.156 .943 .790 .675 .584 .511 .279 700F 1408 939 705 563 470 345 230 46.9 34.5 22.9 17.19 13.74 11.44 9.79 8.56 7.60 6.83 6.20 5.68 5.24 4.86 4.53 3.38 2.68 2.22 1.647 1.301 1.069 .904 .779 .682 .604 .368 .247 800F 1530 1020 765 611 510 375 250 51.0 37.4 24.9 18.69 14.93 12.44 10.65 9.31 8.27 7.44 6.76 6.19 5.71 5.30 4.94 3.69 2.94 2.44 1.812 1.436 1.186 1.006 .872 .768 .684 .432 .305
Flow of Steam
Nomenclature: M Mass rate of flow, lbm/hr P P1 & P2 P V L Average Line pressure, psia Inlet and outlet pressures, psia Rate of pressure drop, psi/1000 ft Specific volume, cu ft/lbm Line length, thousands of feet
Notes: This chart, strictly speaking, gives rate of pressure drop at a point in the pipe, but for relatively small changes in pressure through the pipe will give the average rate of pressure drop if the specific volume at the average pressure is used. Example: Pipe size = 1 1/2 in. sch 40; inlet pressure = 151 psia; outlet pressure = 149 psia; pipe length = 487 ft; temperature = 500F. Find mass flow rate. Average pressure, P = (P1 + P2)/2 = 150 psia. From the table V = 3.68 cu ft per lbm. P= (P1 - P2)/L = (151 -149)/0.487 = 4.1 psi per 1000 ft from the chart, M = 405 lbm per hr.
Chevron Corporation
400-55
March 1997
Fig. 400-25 Flow of Steam2 to 48-Inch Pipe (1 of 2) Table of Specific Volumes for Steam
Specific Volume (V) ABS PRESS. 1 in. Hg 1-1/2 2 2-1/2 3 2 psia 3 14.7 20 30 40 50 60 70 80 90 100 110 120 130 140 150 200 250 300 400 500 600 700 800 900 1000 1500 2000 SAT. TEMP, F 79.1 91.8 101.2 108.7 115.1 126.1 141.5 212.0 228.0 250.3 267.2 281.0 292.7 302.9 312.0 320.3 327.8 334.8 341.3 347.3 353.0 358.4 381.8 401.0 417.3 444.6 467.0 486.2 503.0 518.2 532.0 544.6 596.1 635.6 SAT. STEAM 653 445 340 275 232 174.0 118.9 26.8 20.1 13.74 10.50 8.51 7.17 6.20 5.47 4.89 4.43 4.04 3.72 3.45 3.22 3.01 2.28 1.841 1.541 1.160 .926 .768 .653 .565 .497 .442 .274 .1875 300F 922 615 461 369 308 226 150.5 30.5 22.4 14.82 11.04 8.78 7.26 400F 1142 695 521 417 348 256 170.4 34.6 25.4 16.89 12.62 10.06 8.35 7.13 6.22 5.50 4.93 4.47 4.08 3.75 3.46 3.22 2.36 500F 1163 776 582 466 388 286 190.3 38.8 28.4 18.92 14.16 11.30 9.40 8.04 7.02 6.22 5.58 5.06 4.63 4.26 3.95 3.68 2.72 2.15 1.765 1.283 .990 .792 600F 1287 857 844 515 430 315 210 42.8 31.5 20.9 15.68 12.53 10.42 8.92 7.79 6.92 6.22 5.64 5.16 4.76 4.41 4.11 3.06 2.42 2.00 1.474 1.156 .943 .790 .675 .584 .511 .279 700F 1408 939 705 563 470 345 230 46.9 34.5 22.9 17.19 13.74 11.44 9.79 8.56 7.60 6.83 6.20 5.68 5.24 4.86 4.53 3.38 2.68 2.22 1.647 1.301 1.069 .904 .779 .682 .604 .368 .247 800F 1530 1020 765 611 510 375 250 51.0 37.4 24.9 18.69 14.93 12.44 10.65 9.31 8.27 7.44 6.76 6.19 5.71 5.30 4.94 3.69 2.94 2.44 1.812 1.436 1.186 1.006 .872 .768 .684 .432 .305
Flow of Steam
Nomenclature: M Mass rate of flow, lbm/hr P P1 & P2 P V L Average Line pressure, psia Inlet and outlet pressures, psia Rate of pressure drop, psi/1000 ft Specific volume, cu ft/lbm Line length, thousands of feet
Notes: This chart, strictly speaking, gives rate of pressure drop at a point in the pipe, but for relatively small changes in pressure through the pipe will give the average rate of pressure drop if the specific volume at the average pressure is used. Example: Pipe size = 6 in. 1/4 in. wall; inlet pressure = 151 psia; outlet pressure = 149 psia; pipe length = 487 ft; temperature = 500F. Find mass flow rate. Average pressure, P = (P1 + P2)/2 = 150 psia. From the table V = 3.68 cu ft per lbm. P= (P1 - P2)/L = (151 -149)/0.487 = 4.1 psi per 1000 ft from the chart, M = 13,800 lbm per hr.
Chevron Corporation
400-57
March 1997
This equation can be derived from basic pressure drop relations, but in the literature it is often presented in simplified form with certain empirical components. The two most widely accepted forms are the Weymouth Equation and the Panhandle Equation. The Weymouth Equation, in which friction is a function of the diameter, applies at high Reynolds numbers. The Panhandle Equation, in which friction is a function of the Reynolds number, applies at lower Reynolds numbers. The break point is defined as follows: Re = 9031D2.449
(Eq. 400-28)
where: D = inside diameter, in. The constants (C1 through C7) for the Weymouth and Panhandle equations are shown in Figure 400-25 both as presented in the literature and as derived without empirical components.
Fig. 400-26 Weymouth and Panhandle Equation Constant Equation Weymouth Weymouth Panhandle Panhandle where E Z = = = pipeline efficiency, ranging from 0.94 (new pipe) to 0.88 (old rough pipe) compressibility 0.41P r A 8 1 -------------- + ( 0.29 ) P r (From Source 4.) 4.04 Accurate within 10% if Pr<0.8 and Tr>1.0 T r or if Pr>0.8 and Tr>1.1 Tr16 T/Tc critical temperature, R operating temperature,R P/Pc critical pressure, psia operating pressure, psia C1 433.45 433.50 435.87 503.30 C2 Z 1 E 1 C3 1 1 1.0788 1 C4 2.667 2.667 2.6182 2.695 C5 1 Z 1 1 C6 1 1 0.8539 0.77 C7 0.5 0.5 0.5394 0.565 Source 1 2 3 2
A = Tr = Tc = T = Pr = Pc = P = Sources:
1. Natural Gas Processors Suppliers Association, Engineering Data Book, 1972. 2. Derived by W.A. Ebert, Chevron Engineering Department, 1984. 3. Baumeister and Marks, eds., Standard Handbook for Mechanical Engineers, McGraw-Hill, 1967. 4. Heat Transfer Research Inc., Computer Program Support Volume, pg. E1-47, 1976.
March 1997
400-58
Chevron Corporation
441 Assumptions
The charts in Figure 400-29 (and Lapples analysis) are based on the following assumptions: The friction factor (f) is constant along the length of the pipe. For the entire range of each chart, either the Perfect Gas Law applies or the compressibility factor (Z) and the ratio of specific heats (K) of the gases are constant. The charts are based on horizontal flow through constant cross-sectional area.
The design charts in Figure 400-29 are used when upstream conditions (usually static conditions within a vessel or reservoir) are known and either the discharge rate or downstream pressure are required for a given pipe size. In Figure 400-28, the typical problem is to determine mass flow rate G or pressure P2, given P0, T0, P3, L, and D. The velocity at Section 0 is assumed to be zero.
Chevron Corporation
400-59
March 1997
In Figure 400-29, flow rates are expressed as a ratio of the actual mass velocity, G, to a hypothetical maximum isothermal mass velocity through a nozzle, Gmax. Thus, it is first necessary to calculate Gmax from known conditions:
(Eq. 400-29)
where: G = mass velocity, lbm/ft2sec gc = conversion factor (32.17 lbm ft/lbf sec2) MW = molecular weight, lbm/mole e = base of natural logarithm (2.718) R = gas constant, 1546 ftlbf/lbmoleR T = absolute temperature, R, at location designated by subscript P = absolute pressure, lbf/ft2, at location designated by subscript V = specific volume, ft3/lbm, at location designated by subscript The friction factor (f) must also be established (see Section 410) prior to using the charts, although variations in f affect the answer very little. The initial value of f is usually assumed to be 0.0143 for gas flow.
March 1997
400-60
Chevron Corporation
Fig. 400-29 Design Charts for Gas Flow at High Pressure Drop (1 of 2) 1. To use design charts in Figure 400-29: a. Calculate an overall effective length L of straight pipe of diameter D, including equivalent length for valve and fitting losses (see Section 500). b. Assume a friction factor f for gas flow (usually assumed 0.0143) and calculate fL/D parameter. c. Calculate the hypothetical maximum mass velocity, Gmax, from g c P o 0.5 g c MW 0.5 lb m G max = -----------= P o --------------------------------ev eRT 2 o o ft sec d. Estimate K (ratio of specific heats) from Figure 400-27. e. Enter appropriate chart to determine P2/P0 or G/Gmax and solve for pressure P2 or mass flow. 2. Values for P2/P0 are valid only above the critical pressure ratio line which defines the point of sonic flow and maximum mass flow. Ratio P3/P0 is, however, valid over the entire range shown. Examples: Given: Air within a reservoir at 80F and 200 psig is discharging to the atmosphere through 20 feet of three-inch, schedule 40 pipe which includes two standard 90 long radius elbows. Determine: Discharge rate to the atmosphere Solution: 1. Calculate fL/D parameter (use consistent units) f = 0.0143; assumed D=3.068 in. = 0.256 ft L = 20 + L = 20 + (2)(0.256)(23) = 31.8 ft (see Section 500) 2. Calculate maximum mass velocity, Gmax To = 460 + 80 = 540R Po = (200 + 14.7)(144) = 30,900 lbf/ft2 MW = 29 lbm/mole lbm 0.5 ( 32.17 ) ( 29 ) G max = 30, 900 --------------------------------------------= 627 -----------------------( 2.718 ) ( 1546 ) ( 540 ) 2 ( sec ) ( ft ) 3. Find Flow Rate P3 = (14.7)(144) = 2120 lbf/ft2 P3/Po = 2120/30,900 = 0.0685 G/Gmax = 0.76 (K=1.4) G = (627)(0.76) = 476 lbm/(sec)(ft)2 A = (0.256)2/4 = 0.0515 ft2 Flow Rate = (G)(A) = (476)(0.515) = 24.6 lbm/sec
Chevron Corporation
400-61
March 1997
Fig. 400-29 Design Charts for Gas Flow at High Pressure Drop (2 of 2)
March 1997
400-62
Chevron Corporation
Chevron Corporation
400-63
March 1997
450 References
1. 2. 3. Scott, D. S. Properties of Concurrent Gas-Liquid Flow. Advances in Chemical Engineering, Vol. 4, p.199. New York: Academic Press, 1963. Gouse, W. S., Jr. An Index to the Two-Phase Gas-Liquid Flow Literature. MIT Report No. 9. MIT Press, 1966. Dukler, A. E., M. Wicks III, and R. G. Cleveland. Frictional Pressure Drop in Two-Phase Flow: A. A Comparison of Existing Correlations for Pressure Loss and Holdup. AIChE Journal 10 (1964), p. 38. Lockhart, R. W., and R. C. Martinelli. Proposed Correlation of Data for Isothermal Two-Phase, Two-Component Flow in Pipes. Chemical Engineering Progress, 45 (1949), p. 39. Dukler, A. E., M. Wicks III, and R. G. Cleveland. Frictional Pressure Drop in Two-Phase Flow: B. An Approach Through Similarity Analysis. AIChE Journal 10 (1964), p. 44. Streeter, V. L. Fluid Mechanics. 2nd Edition. New York: McGraw-Hill, 1958. Hughmark, G. A. Holdup in Gas-Liquid Flow. Chemical Engineering Progress Vol. 58 (April 1962), p. 62. Baker, O. Multiphase Flow in Pipelines. Oil and Gas Journal, 10 (Nov, 1958). Griffith, P., and G. B. Wallis. Two-Phase Slug Flow. Journal of Heat Transfer, Transactions of ASME Series C 83 (Aug 1961), p. 307.
4.
5.
6. 7. 8. 9.
10. California Research Corporation Standard Technical Books. California Research Corporation, Richmond, California, 1960. 11. Marks Mechanical Engineers Handbook. 6th Edition. New York: McGrawHill, 1958. 12. Perrys Chemical Engineers Handbook. 4th Edition. New York: McGraw-Hill, 1963. 13. Technical Data Book - Petroleum Refining. New York: American Petroleum Institute, Division of Refining, 1966. 14. S I Engineering Data Book. Tulsa: Gas Processors Suppliers Association, 1987. 15. Lapple, C.E. Isothermal and Adiabatic Flow of Compressible Fluids. Transactions of AIChE, Vol. 39 (1943), pp. 385-432. 16. Loeb, M. B. Graphical Solution of Compressible Fluid Flow Problems. NASA/Kennedy Space Center Document TR-256D, 1965. 17. Loeb, M. B. New Graphics for Solving Compressible Flow Problems. Chemical Engineering, Vol. 76, No. 11 (May 19, 1969). 18. Shapiro, A. H. The Dynamics and Thermodynamics of Compressible Fluid Flow, Vol. I. New York: The Ronald Press Company, 1953.
March 1997
400-64
Chevron Corporation
Chevron Corporation
500-1
March 1997
510 Introduction
Valves and fittings cause more energy loss than pipe of equal axial length. This loss may be relatively insignificant in long lines but, within process plants, it can be a major contributor to system losses. Losses at a change in section take two distinct forms, pressure loss and energy loss. At a well-rounded pipe entrance, there is a pressure loss due to the increase in velocity, but a negligible energy loss. At a pipe exit, pressure change is usually nominal and velocity energy is dissipated as turbulence. Head loss through a valve or fitting can be expressed in two ways: As the number of velocity heads lost As a length of straight pipe with a diameter and pressure drop equal to those of the valve or fitting
where: hf = head loss through the fitting, ft V = average velocity in the line, ft/sec K = constant for the fitting type g = gravitational constant (32.2 ft/sec2) Since V2/2g is the velocity head of the fluid, K is the number of velocity heads lost through the fitting. The average values of K for various valves and fittings are shown in Figure 500-2. K values for various entrance losses are given in Figure 500-1.
March 1997
500-2
Chevron Corporation
Fig. 500-1
Losses Through Entrances and Changes in Section Courtesy of Tube Turns Technologies
Chevron Corporation
500-3
March 1997
Fig. 500-2
March 1997
500-4
Chevron Corporation
where: h = friction head loss, ft f = friction factor L = length of fitting, ft D = diameter of fitting, ft Thus K equals f(L/D). By expressing this relation in terms of L, the loss through the fitting may be expressed as an equivalent length (L) of straight pipe (of the same diameter as the fitting). That is, K L K L = D --- or ---- = --D f f
(Eq. 500-3)
Figure 500-3 gives equivalent length of various sizes of valves and fittings. The L/D ratio provides equivalent length values in terms of diameters of straight pipe, so that one value can be applied to varying diameters of a valve or fitting. The equivalent length ratios shown in Figure 500-2 were calculated from the fittings K values (assumed constant for each type fitting) and a friction factor of 0.025 for liquids and 0.0143 for gases. Because the actual friction factor may differ appreciably, the equivalent length method should be used only for rough estimates or when the total equivalent length for valves and fittings is small compared to the length of straight pipe.
Chevron Corporation
500-5
March 1997
Fig. 500-3
March 1997
500-6
Chevron Corporation
550 Examples
Consult Figures 500-1 and 500-2 for the K and L/D values to be used in the following examples. In practice, manufacturers proprietary fitting loss data should be used whenever available. For elbows under three inches in diameter, increase loss by 30 percent. Add appropriate reducer losses to fitting losses for total loss through a reducing fitting, such as a tee with a reducing branch or a venturi pattern valve.
Solution
Calculate the pressure drops separately for each diameter. Diameter Ratio: DS/DL = 3.068 in./6.065 in. = 0.51 From Figure 500-1:
Fitting Reducer (friction) Reducer (acceleration) TOTAL K 0.16 0.56 0.72
Chevron Corporation
500-7
March 1997
Pipe Area = A = r2 1ft 6.065 2 2 A 1 = 3.14 ------------ in. ------------------2 2 144 in. = 0.201 ft2 1ft 3.068 2 2 A 2 = 3.14 ------------ in. ------------------2 2 144 in. = .051 ft2 Flow rate = Q lb 1 hr ft = 150, 000 ---- ------------------- -------------- -----hr 3600sec 62.4lb .75 = 0.89 ft3/sec Q Velocity = V = --A Q 0.89ft sec V 1 = ------ = ---------------------------2 A1 0.201ft = 4.44 ft/sec Q 0.89ft sec V 2 = ------ = ---------------------------2 A2 0.051ft =17.45 ft/sec
3 3 3 3 2
March 1997
500-8
Chevron Corporation
V One Velocity Head = -----2g V1 4.44 ft sec -------- = ------------------------------------2 2g 2 32.2 ft sec = 0.306 ft
2 2 2 2
V2 17.45 sq ft sec -------- = ----------------------------------------2 2g 2 32.2 ft sec = 4.73 ft Total pressure drop:
62.4 0.75 P 1 = ( 13.29 ) ( 0.306 ) -------------------------- = 1.322psi 144 62.4 0.75 P 2 = ( 0.72 ) ( 4.73 ) -------------------------- = 1.107psi 144 Total = 2.429 psi
(Eq. 500-4)
Chevron Corporation
500-9
March 1997
Solution
From Figures 500-1 and 500-2:
Fitting Gate Valve Elbows (6 x 13) Square Edged Entry (friction) Square Edged Entry (acceleration) TOTAL L/D 8 78 20 40 146
Equivalent length L = (L/D) (D) = (146) (6.065) = 855 in. = 71.2 ft Total equivalent line length 800 + L = 871.2 ft
560 References
1. 2. 3. 4. 5. 6. 7. Flow of Fluids Through Valves, Fittings, and Pipe. Crane Co., Crane Technical Paper No. 410-C, 1984. King, R.C., and S. Crocker. Piping Handbook, 5th Edition, McGraw-Hill, pp 167-181, 1967. Piping Engineering - Fluid Flow in Pipe. Tube Turns Bulletin No. 301, 1951. Simpson, L.L. Process Piping: Functional Design. Chemical Engineering, Vol. 76, No. 8, 1969. Perrys Chemical Engineers Handbook. 4th Edition, McGraw-Hill, New York, 1963. Beck, C. Laminar Flow Friction Losses Through Fittings, Bends, and Valves. Journal American Society Naval Engineers, vol. 56, p. 235-271, 1944. Kittredge, C.P., and D.S. Rowley. Resistance Coefficients for Laminar and Turbulent Flow Through One-Half-Inch Valves and Fittings. ASME Transactions, Vol. 79, pp 1759-1766, 1957.
March 1997
500-10
Chevron Corporation
Chevron Corporation
600-1
January 1990
610 Introduction
The calculation of pressure drop in noncircular conduits is handled differently for laminar and turbulent flow. Turbulent flow boundary layers are thin and relatively unaffected by proximity to the conduit walls. Laminar boundary layers, however, are thick, and the boundary layers from opposite walls often interact. Pressure drop for turbulent flow can be closely approximated based on the calculation of the conduits hydraulic diameter given in Section 620. For laminar flow, empirical data is needed to arrive at a reasonable approximation of the pressure drop for a specific case. Section 630 provides empirical correlations for the friction factor for rectangular and concentric annulus geometry. For this analysis, the transition between laminar and turbulent flow can be assumed to be at a Reynolds number of 2000.
VD h VD h Re = -------------- = ---------
(Eq. 600-2)
V dh = fL ----------------D h 2g
(Eq. 600-3)
where: f = friction factor from Figure 400-1, Moody Chart Dh = hydraulic diameter, ft Ax = cross-sectional flow area, ft2 Pw = wetted perimeter of channel, ft Re = Reynolds number = density, lbm/ft3 V = velocity, ft/sec
January 1990
600-2
Chevron Corporation
= kinematic viscosity, ft2/sec = viscosity, lb sec/ft2 f = friction factor dh = pressure drop in head loss, ft L = conduit length, ft g = gravitational constant, 32.17 ft/sec2
For concentric tube annulus choose C2 such that: a = radius of inner tube b = radius of outer tube Calculate the friction factor as follows: C2 f = -----Re
(Eq. 600-5)
Chevron Corporation
600-3
January 1990
Fig. 600-1 a/b(1) 0.0 0.025 0.05 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Laminar Flow Pressure Drop Constants C1(1) 96.0 94.0 90.0 85.2 76.8 70.4 65.6 62.8 60.0 58.8 57.6 57.2 56.8 C2(1) 64.0 81.2 86.0 89.2 92.4 93.6 94.8 95.2 95.6 96.0 96.0 96.0 96.0
640 References
Kays and Crawford, Convective Heat and Mass Transfer, McGraw- Hill, 1980.
January 1990
600-4
Chevron Corporation
Chevron Corporation
800-1
January 1990
810 Introduction
If a valve is closed rapidly in a line containing flowing liquid, the inertia of the flowing liquid will increase the pressure at the valve. This effect is called surge, and the increase in pressure is called surge pressure. Surge can cause extremely rapid changes in pressurerapid enough to cause the metallic percussions commonly called water hammer. The surge pressure wave will then propagate back up the line, and may cause mechanical damage. Water flowing at 10 ft/sec can generate a surge pressure rise of about 500 psi. Bulk modulus values for hydrocarbons are generally lower than for water, but surge pressures are still significant considerations in designing hydrocarbon piping systems. See Figures 800-1 through 800-3. This section provides a method for approximating the maximum surge pressure in a simple system. Because of nonlinear elements in the analysis, a more thorough calculation of surge pressure can be a complex problem. See Section 840 for sources providing more general solution techniques.
January 1990
800-2
Chevron Corporation
Fig. 800-1
Average Bulk Modulus for Crude Oil, Fuel Oil, Gas Oil, and Gasoline
Chevron Corporation
800-3
January 1990
Fig. 800-2
January 1990
800-4
Chevron Corporation
Fig. 800-3
Fig. 800-4
SurgeSimple Case
Chevron Corporation
800-5
January 1990
where: = speed of sound through liquid in pipe, ft/sec K = bulk modulus of liquid, psi. For hydrocarbon liquids, see Figures 800-2 and 800-3; for water, see Figure 800-4. = density of liquid, lbm/ft3 g = 32.2 ft lbm/sec2 lbf D = inside diameter of pipe, inches t = wall thickness of pipe, inches E = modulus of elasticity of pipe material, psi C = constant which depends on pipe fixity = 0.91 for line anchored against axial movement = 0.95 for unrestrained line A pressure disturbance generated at the valve will propagate back to the vessel and return to the valve in a propagation time equal to 2L/ (where L = line length between vessel and valve in feet). If the valve closing time (T) is less than 2L/, the surge pressure can be approximated by V 2L P = --------------- for T < ----- 144g
(Eq. 800-2)
where: P = surge pressure, psi V = total change in velocity, ft/sec T = valve closing time, sec 2L/ = propagation time, sec This solution is only an approximation tailored to this simple case. For example, this equation is not valid if the valve closing time is greater than 2L/. Section 840, references 1, 2, and 3, presents general techniques for calculating surge pressure accurately and in more complex situations.
January 1990
800-6
Chevron Corporation
840 References
1. 2. 3. Symposium on Surges in Pipelines, The Institution of Mechanical Engineers, Proceedings 1965-66, Vol. 180, Part 3E. Hydraulic Transients, Rich, G. R., Dover Publications, Inc., New York, 1963. Hydraulic Transients, Streeter, V. L., Wylie, E. B., McGraw-Hill, 1967.
Chevron Corporation
800-7
January 1990
Chevron Corporation
900-1
January 1990
910 Introduction
Section 900 addresses pipeline flow situations in which large temperature changes significantly affect fluid properties and flow characteristics. Other situations, involving typical liquids and gases at close to ambient temperatures or with small temperature changes, can be adequately addressed using the methods of Section 400. For long pipelines carrying fluids that require high pumping energy, the effects of friction heating should be investigated. Section 930 defines the relationship of temperature change to pumping energy for viscous fluids. Similarly, at high pressures, a pressure correction to viscosity may be necessary, as discussed in Section 940. Section 960 discusses heat transfer between the pipeline and its surroundings, including calculation of external heat transfer coefficients for pipelines in various ambient conditions. Section 950 identifies computer programs available for solving difficult temperature/flow problems over the length of a pipeline.
(Eq. 900-1)
January 1990
900-2
Chevron Corporation
BPD = oil flow rate, BPD r = oil density, lbm/ft3 Cp = oil specific heat, BTU/lbm F When Cp = 0.5 and = 58 the expression becomes:
High viscosities increase the ratio of HP to BPD and therefore increase the temperature change due to friction. How well the pipeline is insulated determines how much of the added heat will actually stay in the oil. The following factors determine the decrease in required pumping power: How much the viscosity is decreased by the increased temperature How sensitive the flow regime is to decreases in viscosity. Pressure drop in laminar flow is a stronger function of viscosity than in turbulent flow. Pressure drop in transition flow is not a function of viscosity
The effect of friction heating generally increases with: Flow rate Viscosity Insulation Line length
where:
(Eq. 900-3)
Chevron Corporation
900-3
January 1990
P = pressure, psig As with friction heating, the effect of increased viscosity on pumping energy requirements depends on the flow regimes sensitivity to viscosity. Since some pipelines have more than one flow regime (laminar, transition, or turbulent), the change in pumping requirements can be difficult to calculate without a computer program.
January 1990
900-4
Chevron Corporation
Buried Pipelines
The appropriate ambient temperature value for a buried line is the yearly average air temperature. The external heat transfer coefficient for buried pipe can be calculated as follows:
(Eq. 900-4)
where: ho = external heat transfer coefficient, Btu/hr ft2 F k = soil thermal conductivity, Btu/hr ft F D = pipe outside diameter, inches d = virtual pipe burial depth, inches = da + 12 k/ha where: da = actual pipe depth to center line, inches ha = ground to air heat transfer coefficient, Btu/hr ft2 F Ground-to-air heat transfer coefficients are typically 1 to 3 Btu/hr sq ft F for low to moderate winds. Soil thermal conductivity is mainly a function of moisture content. Typical values are between 0.2 and 1 Btu/ hr ft F. Figure 900-2 gives conductivity values for some soil, sand, and rock types, and other related materials. It also shows the relationship between soil density and thermal conductivity.
Chevron Corporation
900-5
January 1990
Ground moisture tends to migrate away from heated objects. Therefore, the actual soil conductivity around a buried hot (or warm) pipeline may vary with time and with distance from the line.
(Eq. 900-5)
where: ho = external heat transfer coefficient, Btu/hr ft2 F ha = ambient fluid heat transfer coefficient, Btu/hr ft2 F R1 = thermal resistance of layer 1, hr ft2 F/Btu R2 = thermal resistance of layer 2, hr ft2 F/Btu rop = outside radius of outermost layer on pipe, ft ro1 = outside radius of layer 1, ft ro2 = outside radius of layer 2, ft ri1 = inside radius of layer 1, ft ri2 = inside radius of layer 2, ft k1 = thermal conductivity of layer 1, Btu/hr ft F
January 1990
900-6
Chevron Corporation
k2 = thermal conductivity of layer 2, Btu/hr ft F ln = natural logarithm Thermal conductivity for subsea concrete (k1) coatings is about 0.5 Btu/hr ft F. Subsea ambient heat transfer coefficients (ha) are in the low one-hundreds for moderate currents. Using ha equal to 150 Btu/hr sq ft F should give acceptable accuracy. Figure 900-1 shows approximate ambient heat transfer coefficients for air.
970 References
1. Petroleum Refining, Technical Data Book. Washington, D.C.: American Petroleum Institute, 1970, pp 11-47.
Fig. 900-1
Chevron Corporation
900-7
January 1990
Soil Conductivity Chart (1 of 4) Conductivity, Btu/hr ft F Moisture Content dry 0.2 0.24 0.31 0.40 0.51 0.61 0.72 0.16 0.20 0.14 0.17 0.20 0.16 0.74 0.28 0.17 0.40 0.20 0.56 0.26 0.35 0.51 0.79 0.60 0.50 0.59 0.73 0.58 0.72 2% 4% 6% 8% 10% 12% 14% 20% 30% 1 2 2 2 2 2 1 2 2 2 3 Source
(dry density where reported) Soil Soil (80 lb/cu ft) Soil (90 lb/cu ft) Soil (100 lb/cu ft) Soil (110 lb/cu ft) Sandy soil Sand White sand, clean Yellow sand, clean Yellow sand and clay Clay Fine crushed quartz (100 lb/cu ft) (110 lb/cu ft) Crushed quartz (100 lb/cu ft) (110 lb/cu ft) (120 lb/cu ft) Graded Ottawa sand (100 lb/cu ft) (110 lb/cu ft) Fairbanks sand (100 lb/cu ft) (110 lb/cu ft) (120 lb/cu ft)
1.00 1.33
4 4
4 4 4
0.83 1.17
4 4
4 4 4
January 1990
900-8
Chevron Corporation
Soil Conductivity Chart (2 of 4) Conductivity, Btu/hr ft F Moisture Content dry 2% 4% 6% 8% 10% 12% 14% 20% 30% Source
(dry density where reported) Lowell sand (100 lb/cu ft) (110 lb/cu ft) Chena river gravel (110 lb/cu ft) (120 lb/cu ft) Crushed feldspar (100 lb/cu ft) (110 lb/cu ft) (120 lb/cu ft) Crushed granite (100 lb/cu ft) (110 lb/cu ft) (120 lb/cu ft) Dakota sandy loam (110 lb/cu ft) (120 lb/cu ft) Crushed trap rock (100 lb/cu ft) (110 lb/cu ft) (120 lb/cu ft) Ramsey sandy loam (100 lb/cu ft) (110 lb/cu ft)
4 4
0.75 1.08
4 4
4 4 4
4 4 4
4 4
4 4 4
4 4
Chevron Corporation
900-9
January 1990
Soil Conductivity Chart (3 of 4) Conductivity, Btu/hr ft F Moisture Content dry 2% 4% 6% 8% 10% 12% 14% 20% 30% Source
(dry density where reported) Northway fine sand (100 lb/cu ft) (110 lb/cu ft) Northway sand (100 lb/cu ft) (110 lb/cu ft) Healy clay (90 lb/cu ft) (100 lb/cu ft) (110 lb/cu ft) Fairbanks silt loam (90 lb/cu ft) (100 lb/cu ft) (110 lb/cu ft) Fairbanks silty clay loam (90 lb/cu ft) (100 lb/cu ft) (110 lb/cu ft) Northway silt loam (90 lb/cu ft) (100 lb/cu ft) (110 lb/cu ft) Iraq Desert Steppe Iraq Desert Sand
4 4
4 4
0.67 0.83
4 4 4
0.42
0.62 0.83
4 4 4
0.75
0.42
0.62 0.79
4 4 4
0.75
0.33
0.50 0.58
4 4 4 5 5
January 1990
900-10
Chevron Corporation
Soil Conductivity Chart (4 of 4) Conductivity, Btu/hr ft F Moisture Content dry 1.06 2% 4% 6% 8% 10% 12% 14% 20% 30% 5 Source
(dry density where reported) Abqaiq sand (in a limestone trench) Earth, coarse gravelly Concrete Common red brick Granite Limestone Marble Sandstone Material
0.30 0.54 0.36 1.73-3.98 1.26-1.33 2.07-2.94 1.83 Conductivity, Btu/hr ft F 0.25 0.35 0.65 1.00 6 6 6 0.70
3 1 2
6 Source 7 7 7 7
Soil (fairly dry, avg. California summer) Soil (wet weather, some drainage) Soil (heavy rains, but ground not flooded) Soil (marshy or constantly soaked) Sources
1. Krieth, F. Principles of Heat Transfer. New York: Harper & Row, 1973. 2. Flow of Hot Oil in Pipelines. Various sources. Chevron: discontinued. 3. Eckert, E. R. G., and R. M. Drake, Jr. Heat and Mass Transfer. New York: McGraw-Hill, 1959. 4. McAdams, W. H. Heat Transmission, 3rd ed. New York: McGraw-Hill, 1954. 5. 6.
Journal of the Institute of Petroleum, Vol. 36, No. 321, September, 1950.
Holman, J. P. Heat Transfer, 5th ed. New York: McGraw-Hill, 1981.
Chevron Corporation
900-11
January 1990
Chevron Corporation
1000-1
March 1997
1033 Constitutive Relationships 1034 Calculation of Flow Parameters 1035 Hydraulics Equations 1036 Computer Program HOTOIL 1037 Estimating Pipeline Restart Pressure Gradient 1038 Wax Deposition 1040 References 1041 Viscosity Conversion 1042 Viscosity Data 1043 Brine-in-Oil Emulsions 1000-54
March 1997
1000-2
Chevron Corporation
1010 Viscosity
Viscosity is a measure of the internal friction or resistance of a fluid to the relative motion of its parts. It may be regarded as the relationship between the force applied to a fluid and the rate of deformation produced in the fluid.
(Eq. 1000-1)
The coefficient is defined as the absolute, or dynamic, viscosity. Its metric system dimensions are as follows:
(Eq. 1000-2)
where: F = shearing force, dynes (dyne = gram cm/sec2) V/D = velocity gradient, sec-1 A = area of shear, cm2 Since the poise is a relatively large number, absolute viscosity is normally expressed in centipoise (0.01 poise). The English system expression for absolute viscosity is as follows:
(Eq. 1000-3)
or
(Eq. 1000-4)
Conversion factors for absolute viscosities are included as Figures 1000-1 through 1000-5.
Chevron Corporation
1000-3
March 1997
March 1997
1000-4
Chevron Corporation
Fig. 1000-2 Viscosity ConversionFrom Various Terms to Saybolt Universal (1 of 2) Courtesy of Hydraulic Institute
Chevron Corporation
1000-5
March 1997
Fig. 1000-2 Viscosity ConversionFrom Various Terms to Saybolt Universal (2 of 2) Courtesy of Hydraulic Institute
March 1997
1000-6
Chevron Corporation
Chevron Corporation
1000-7
March 1997
March 1997
1000-8
Chevron Corporation
Chevron Corporation
1000-9
March 1997
March 1997
1000-10
Chevron Corporation
Fig. 1000-4 Viscosity Conversion Equations Equation for Converting Kinematic Viscosity to Flow Times(1) 1 + EV T = DV + --------------------------------------------2 3 F + GV + HV + IV Unit Saybolt Universal(2) seconds at 100 F Saybolt Universal(2) seconds at 210 F Saybolt Furol seconds at 122 F Saybolt Furol seconds at 210 F Redwood No. 1 seconds at 140 F Redwood No. 2 seconds Engler degrees D 4.6324 4.6635 0.47170 0.47916 4.0984 0.40984 0.13158 E 0.03264 0.00677 0.0 0.0 0.0 0.0 0.0 where: T = flow time, seconds (or Engler degrees) V = kinematic viscosity, centistokes (CS) D,E,F,G,H,I = constants given below F 0.039302 0.039911 0.4895 0.3797 0.038014 0.38014 1.1326 G 0.02627 0.000938 -0.005213 0.0 0.001919 0.01919 0.01040 H 0.002397 0.000280 0.0000718 0.0001783 0.0000278 0.000278 0.00656 I 0.0000164 6 0.0000027 4 0.0 0.0 0.0000052 1 0.000521 0.0 Range (CS) >1.8 >1.8 >48 >48 >40 >73 >1.0
Equation for Converting Flow Times to Kinematic Viscosity(1) BT V = AT -------------3 T +C Unit Saybolt Universal seconds at 100 F Saybolt Universal seconds at 210 F Saybolt Furol seconds at 122 F Saybolt Furol seconds at 210 F Redwood No. 1 seconds at 140 F Redwood No. 2 seconds Engler degrees A 0.21587 0.21443 2.120 2.087 0.244 2.44 7.60 where: V T A,B,C = kinematic viscosity, centistokes (CS) = flow time, seconds (or Engler degrees) = constants given below
Range SUS >32 SUS >32 SFS >25 SFS >25 R1 >35 R2 >31 E >1.000
(1) Figure 1000-5 is based on these equations and should provide equal accuracy. (2) The following correction for Saybolt Universal seconds at other temperatures, is small and usually unnecessary:
Chevron Corporation
1000-11
March 1997
March 1997
1000-12
Chevron Corporation
where: = kinematic viscosity, stokes = absolute viscosity, poises = density, grams/cm3 The most commonly used unit of kinematic viscosity is the centistoke (cs = 0.01 stoke). The corresponding unit in English units is ft2/sec. To convert centistokes to English units, multiply by 1.08 x 10-5 to obtain ft2/sec.
Chevron Corporation
1000-13
March 1997
March 1997
1000-14
Chevron Corporation
Chevron Corporation
1000-15
March 1997
March 1997
1000-16
Chevron Corporation
Chevron Corporation
1000-17
March 1997
March 1997
1000-18
Chevron Corporation
Fig. 1000-11 Viscosity vs. Temperature, Refined ProductsAsphalts & Fuel Oils
Chevron Corporation
1000-19
March 1997
March 1997
1000-20
Chevron Corporation
Chevron Corporation
1000-21
March 1997
March 1997
1000-22
Chevron Corporation
The viscosities of petroleum products are generally specified in terms of time, and using one of the following viscometers: Saybolt Universal for lubricating oils, gas oils, and crude oils. Saybolt Universal viscosity is the time in seconds needed for the delivery of 60 cc of liquid from the Saybolt viscometer. Temperatures usually employed are 70F, 130F, and 210F. Saybolt Furol for crude residuum and heavy fuel oils. Saybolt Furol viscosity is the time in seconds needed for the delivery of 60 cc of liquid from the Saybolt viscometer equipped with a Furol outlet tube. The test is usually run at a temperature of 122F. The outflow time of the Furol instrument is approximately one-tenth that of the Universal.
Although Saybolt viscometers are used in petroleum product specification, capillary instruments are more widely used in making viscosity measurements. Capillary instruments give more precise and reproducible results which can be converted to centistokes or Saybolt seconds with the appropriate charts. Test procedures for kinematic viscosities and descriptions of several commercially available viscometers are given in ASTM Standard D-445-65. Equations and tables for converting kinematic viscosity to Saybolt Universal viscosity and Saybolt Furol viscosity are given in Figures 1000-4 and 1000-5.
Chevron Corporation
1000-23
March 1997
March 1997
1000-24
Chevron Corporation
Chevron Corporation
1000-25
March 1997
1. This chart can be used to estimate the viscosity of a blend of two oils where their viscosities and volume fractions for a given temperature are known. 2. Reasonably accurate answers can be expected only where the two oils being blended are of the same type. 3. Viscosity scales on the chart are based on Chevron Research Companys blending factor method in which the viscosity scale is drawn proportional to: ln [ l + ln (2)] Where equals kinematic viscosity in centistokes. Example: Given: Problem: Solution: A blend is to be made with 80% consisting of an oil with a 10,000 centistoke viscosity and 20% of an oil with a 100 centistoke viscosity. Determine the viscosity of the blend. 1. Connect the point for the higher viscosity oil (10,000 cs) on the left-hand scale with the point for the lower viscosity oil (100 cs) on the right hand-scale with a straight edge. 2. Read the viscosity of the blend at the percentage composition of the blend. For this example, the viscosity of the blend is 3,200 cs.
March 1997
1000-26
Chevron Corporation
Chevron Corporation
1000-27
March 1997
1. This chart may be used to estimate the order of magnitude of the effective viscosity of a brine-in-oil emulsion at a given temperature if the viscosity of the clean oil at that temperature and the percent brine in the emulsion are known. 2. To estimate an emulsion viscosity: Classify the emulsion as follows: Tight. Small brine particle sizemight occur in a wet well where a high gas-oil ratio results in extremely agitated flow. Medium. Brine particle size between that of tight and loose emulsionsmight occur if a pump is in poor condition or under moderate gas agitation. Loose. Large brine particle sizemight occur in production from a pumping well in which mechanical equipment is in good condition. Find the point on the appropriate curve for percent brine in emulsion and read the corresponding viscosity ratio. Calculate the viscosity of the emulsion by multiplying the viscosity of the clean oil by the viscosity ratio.
3. In general, emulsions are thixotropic fluids with appreciable shear strength. As the percent brine approaches the inversion point (often about 79%) the shear strength of the emulsion may be very high and the emulsion may take on a plastic-solid quality. 4. At the inversion point a brine-in-oil emulsion will invert to an oil-in-brine emulsion. Viscosity of the oil-inbrine emulsion will be approximately that of the brine. Example: Given: Problem: Solution: A tight emulsion of 40% brine in 12.3 API Inglewood crude is flowing at 100F. Determine the viscosity of the emulsion. From Figure 1000-17, viscosity ratio = 5.
From chart on Figure 1000-8, viscosity of clean oil = 34 cs. Viscosity of emulsion = (viscosity of clean) (viscosity ratio). Viscosity of emulsion = (34 cs) (5) = 170 cs.
March 1997
1000-28
Chevron Corporation
Chevron Corporation
1000-29
March 1997
March 1997
1000-30
Chevron Corporation
Chevron Corporation
1000-31
March 1997
1020 Gravity
In the petroleum industry the mass-volume relationship most commonly used for gases and liquids is gravity. Gravity can be expressed in several ways, for instance, specific gravity and degrees API. The density of a fluid, , is defined as its mass per unit volume. =m/v
(Eq. 1000-6)
The density of water at standard conditions of 14.7 psia and 39.4F (4C) is 62.43 1bm/cu ft or 1000 kg/m3. The standard conditions for defining water are not equivalent to the more commonly used, which are 14.696 psia and 60F. The specific weight, , of a unit volume of a material depends on its location in the earths gravity field. Thus: =g
(Eq. 1000-7)
Specific gravity (SG) is the ratio of the weight of a substance to the weight of an equal volume of water (or air, density 0.0763 lbm/cu ft) at standard conditions. Since specific gravity is related to standard conditions it does not change with location. Specific gravity may also be expressed as a ratio of density. API gravity is expressed in degrees and is related to specific gravity by the equation shown in Figure 1000-23. Other gravity scales exist and conversion tables are available. Since specific gravity is based on water at a standard temperature, substances which change volume with temperature have varying specific gravities. Figure 1000-20 shows the temperature relationship for oils. In Figure 1000-20: Specific gravities are referred to water at 60F The temperature correction to specific gravity is based on the ASTM-IP Petroleum Measurement Tables, Table 24
Figure 1000-19 gives API gravities and corresponding volume measurements for various refined petroleum products. Figures 1000-21 and 1000-22 show the thermal expansion of petroleum fractions. Figure 1000-21 provides a fast means of solving for resultant gravity after blending two oils of different gravities. The graph reproduces the equation on Figure 1000-23. Values obtained from Figure 1000-21 are approximate since no allowance has been made for shrinkage of the blend. Where more accurate data is required, a test should be made of the stocks in question. Specific gravities in Figure 1000-21 are referred to water at 60F. For use of the chart all gravity data should be at the same temperature.
March 1997
1000-32
Chevron Corporation
Chevron Corporation
1000-33
March 1997
March 1997
1000-34
Chevron Corporation
Chevron Corporation
1000-35
March 1997
March 1997
1000-36
Chevron Corporation
Fig. 1000-23 API and Baum Gravity Tables Courtesy of Goulds Pumps
Chevron Corporation
1000-37
March 1997
1021 Example
Find the gravity of a blend consisting of 21-1/2% 42API gravity oil and 78-1/2% 23.50API gravity oil. On Figure 1000-21, lay a straightedge connecting the specific gravity of the heavier oil at the top of the page (0.913) with the specific gravity of the lighter oil at the bottom (0.815). The intersection of this line with the composition of the blend gives the gravity of the blend: 0.892 specific gravity, or 27.1API gravity.
March 1997
1000-38
Chevron Corporation
during the drill stem test, keeping the separator and storage facilities hot, and taking the samples as late in the drill stem test as possible. Sample drums should be clean, preferably epoxy-coated, and (if possible) filled from flowing lines rather than from storage tanks. After the residual gas has been allowed to boil off, the drums should be sealed tightly for shipment to prevent loss of light ends and intrusion of water. In the Chevron Research laboratory the drums are heated and stirred to melt any wax that may have precipitated and deposited on the drum surfaces. This is done in a way that minimizes the loss of light ends. The crude is then put in smaller containers for ease of handling. The collection and handling process imposes a different history on the crude sample than would the production process. Before any testing is done, this history is erased, usually just before the test, by heating the sample enough to ensure that all the wax is back in solution. This temperature is often close to the reservoir temperature.
Characterization of Crude
Standard laboratory tests are normally used to characterize a crude sample. These tests usually include gravity, pour point, viscosity (in a capillary viscometer at higher temperatures), and water content. If the water content is higher than that expected from the production separators it is removed before the testing begins, because water can have a significant effect on the crude properties. Additional testing includes wax content, microscopic cloud point, specific heat capacity, distillation, and toluene and pentane insolubles. These tests are only done if required for a specific purpose.
Chevron Corporation
1000-39
March 1997
The temperature profile, T(t), is usually defined by the initial temperature, Ti, ambient temperature, Ta, and cooling parameter, , where T(t) = Ta + (Ti - Ta) e-t
(Eq. 1000-8)
If the time, t, is in hours then the parameter, , is in hours-1. As oil flowing in a pipeline cools, the shearing profile changes, increasing from the inlet to the outlet. This shearing (or turbulence) can be characterized by the power dissipation density (in horsepower/cubic foot). Although this change can be simulated in the laboratory, fluid properties are usually not very dependent on the amount of turbulence, and this parameter usually remains fixed in a test. In a typical test of consistency, an oil sample heated to Ti is placed in a device that simulates the shear. This device is immersed in a temperature-controlled bath and cooled according to the predefined temperature profile (Eq.1000-8) while being subjected to shear at a specified rate. When the oil temperature reaches the range of interest, the oil sample is transferred to a concentric cylinder viscometer that is also immersed in the temperature-controlled bath. Cooling continues while the sample is sheared in the viscometer. The temperature in the bath is controlled to within a fraction of a degree. This is an important feature of the laboratory procedure. Concentric cylinder viscometers come in various styles. The Ferranti model favored by Chevron Research is shown in Figure 1000-24. The oil sample is contained in the annular gap between the two cylinders. The outer cylinder is rotated at a controlled speed and the inner one is fixed to a torsion balance. The cylinder diameters, the annular gap between the cylinders, and the rotation speed of the outer cylinder define the shear strain rate imposed on the oil in the gap. The cylindrical area of the inner cylinder and the torque on the torsion balance define the resulting shear stress on the oil. The air bubble under the inner cylinder eliminates end effects. At various temperatures along the cooling profile, the rotational speed of the outer cup is varied to measure shear stress/shear strain rate pairs over a range of shear strain rates. In the Ferranti viscometer this is done at three shear strain rates. In other viscometers pairs can be measured at up to eight shear strain rates. During these shear scans the shear strain rate is held constant long enough to let transient or thixotropic effects decay. The resulting pairs are called equilibrium values. Chevron has found that these measurements can be quite misleading if the thixotropic effects often associated with waxy crude oils are not permitted to decay before the test values are recorded. The temperature range for these consistency tests is usually from the neighborhood of the pour point down to the lowest temperature of interest or to the point where the viscometer goes off scale. The use of the shear stress, shear strain rate pairs produced by these tests is discussed in the section on calculation of flow parameters.
March 1997
1000-40
Chevron Corporation
where:
Chevron Corporation
1000-41
March 1997
D = inside diameter of coil p = pressure difference across the coil L = length of the coil In general, the BAYS at any temperature Ta depends on the amount of static cooling T = Tt - Ta. A test program usually includes measuring the BAYS as a function of T. For many crudes, there is a value for T that gives a maximum value for BAYS.
Newtonian liquids have a certain constitutive relationship: the shear stress is proportional to the shear strain rate. The proportionality constant (or slope) of the Newtonian relationship is the absolute viscosity , a function of temperaturethe lower the temperature, the higher the slope. For Newtonian liquids we can express the relationship between shear stress, , and shear strain rate, s, as = s
(Eq. 1000-10)
March 1997
1000-42
Chevron Corporation
If the shear strain rate and the corresponding shear stress are known, then the absolute viscosity is = /s
(Eq. 1000-11)
The absolute viscosity of Newtonian fluids is a function of temperature only. Furthermore, Newtonian liquids will flow even when a very small shear stress is applied. In general, there is no reason for the constitutive relationship for liquids to be straight and to go through the origin. There are many such non-Newtonian relationships or flow models. In Figure 1000-26 the Power Law (or Pseudoplastic) Model goes through the origin but is curved. The Bingham Model is straight but does not go through the origin. The Herschel-Bulkley Model is neither straight nor does it go through the origin. The equations describing the relationship between shear stress, , and shear strain rate, s, for these models are as follows: Power Law Model: = sn
(Eq. 1000-12)
Bingham Model: = y + s
(Eq. 1000-13)
Herschel-Bulkley Model: = y + sn
(Eq. 1000-14)
For Equations 1000-13 and 1000-14 the parameter y is the intercept of the curve with the vertical axis. This intercept yield strength is not the same as the BAYS, and one cannot be estimated from the other. Often the BAYS, b, is larger than y. In addition, power law liquids have no intercept yield strength but can have a BAYS. Equations 1000-12 and 1000-13 are two-parameter models; to define the flow model one needs either (, n) or (, y). Equation 1000-14 requires three parameters to define the flow, (, n, and y). The Herschel-Bulkley model becomes the power law model when y = 0 and the Bingham model when n = 1. All three become the Newtonian model when y = 0 and n = 1. Heat loss and pressure drop can be calculated for any of these models. All of the flow parameters just discussed depend upon temperature. With nonNewtonian liquids one can also talk about the apparent viscosity, a = /s
(Eq. 1000-15)
Chevron Corporation
1000-43
March 1997
Apparent viscosity is a function both of temperature and shear strain rate; as the shear strain rate increases, the apparent viscosity decreases. This leads to the notion of shear thinning liquids. In contrast to the Newtonian and power law models, Bingham and Herschel-Bulkley liquids will not flow when a small shear stress is applied, but only when the applied shear stress exceeds the intercept value y. The study of flow models and the measurement of flow properties is the science of rheology. Often the flow properties and flow model for a liquid are called the rheology of the liquid. As we shall see, the selection of the flow model comes from an analysis of the laboratory data. When plotted as shear stress versus shear strain rate, the laboratory data suggest a model that gives the best fit of the data.
March 1997
1000-44
Chevron Corporation
The Newtonian curves on the right side of Figure 1000-27 are calculated with Equation 1000-10, s being the parameter of each curve. The circled points are typical data from a laboratory viscometer. The task is to draw smooth curves for each shear strain rate through these circled points, blending into the Newtonian curves on the right. These curves are also shown. Using the curves in Figure 1000-27, (,s) pairs can be read for any temperature. These pairs can be plotted on either linear or log-log graph paper. On linear graph paper, Bingham liquids will appear as a straight line. On log-log graph paper, Power Law liquids will appear as a straight line. The flow model selected should give the best representation of the data over the shear strain rate range of interest. Pairs from Figure 1000-27 at T = T1 are plotted in Figure 1000-28. Using least squares regression techniques, the parameters for Bingham and Power Law equations can be calculated. The resulting equations are Power Law: = 14.09s0.355 r2 = 0.998 Bingham: = 20.49 + .837s r2 = 0.973 The term r2 is a measure of how well the curve fits the data; the closer the value is to unity, the better the fit. For these points, the Power Law is a better fit. The same process can be applied at several temperatures. Using Figure 1000-27, three data pairs can be selected at any temperature and the parameters calculated for the selected flow model. These parameters can then be plotted against temperature as shown in Figure 1000-29. For this plot, the absolute viscosity-like term, , is converted to a kinematic viscosity-like term, , using the inversion of Equation 1000-16: (T) = (T)/SG(T)
(Eq. 1000-17)
Equation 1000-17 relates the data more easily to the higher temperature Newtonian measurements and is easier to use with HOTOIL. The term viscosity-like is used because, in the non-Newtonian range, the factor is not always given in units of viscosity. For the same reason the units of in Figure 1000-29 are stokes. The variable n in Figure 1000-29 is the second parameter for the selected flow model. See also Figure 1000-31. The curves in Figure 1000-29 characterize the constitutive relationship for the liquid over the temperature range of interest. These curves, in tabular form, will be used as input data for HOTOIL. Before discussing this program, we must discuss the equations needed for the hydraulics calculations.
Chevron Corporation
1000-45
March 1997
March 1997
1000-46
Chevron Corporation
Chevron Corporation
1000-47
March 1997
March 1997
1000-48
Chevron Corporation
(Eq. 1000-18)
In Figure 1000-30, the two p(V) equations are identical, using the Fanning friction factor, f, to relate the p term to the V term. For the laminar flow range, the relationship between f and the Reynolds number is also identical for Newtonian and NonNewtonian. The only difference is in the definition of the generalized Reynolds number which, for non-Newtonian, is expressed in terms of the generalized flow parameters, K and n. For the turbulent flow range, Figure 1000-30 shows the relationship between the friction factor and the Reynolds number. Note that the generalized non-Newtonian equation collapses to the Newtonian equation when n = 1 (When n = 1, K = ). Similarly, the generalized Reynolds number collapses to the Newtonian Reynolds number when n = 1. The equations in Figure 1000-30 for the turbulent flow range friction factor are for smooth tubes. A correlation relating f to Reg for rough tubes has not been published. To estimate the non-Newtonian rough tube turbulent friction factor, HOTOIL uses an approximation:
(Eq. 1000-19)
Chevron Corporation
1000-49
March 1997
p - V:
Reynolds No.:
p L D V f
= = = = =
pressure drop pipeline length pipeline inside diameter liquid velocity friction factor
d gc K, n
= = = =
liquid density gravity conversion factor liquid viscosity generalized flow parameters (see Figure 1000-31)
where: fr = rough-tube friction factor at Reg fs = smooth-tube friction factor at Reg from equation in Figure 1000-30 fnr = rough-tube friction factor for Newtonian liquids at Re = Reg fns = smooth-tube friction factor for Newtonian liquids at Re Newtonian friction factors are calculated with the Colebrook equation. The final element needed to make Figure 1000-30 complete is a criterion for establishing the transition from laminar to turbulent flow. This is provided in a paper by Dodge and Metzner and is reproduced in texts on non-Newtonian flow. The minimum Reynolds number for the turbulent flow regime can be approximated as follows: Ret = 4000 (1 - 0.25 n)
(Eq. 1000-20)
In the computer program the friction factor in the critical range between laminar and turbulent flow is assumed to be constant at fcr, where fcr is the turbulent friction factor at Ret. The laminar range is assumed to stop at Rel where
March 1997
1000-50
Chevron Corporation
(Eq. 1000-21)
This is the Reynolds number at the intersection of the laminar friction factor curve and the horizontal line f = fcr. The equations in Figure 1000-30 are valid for any liquid model. The liquid properties enter the equations through the terms K and n. For each of the models discussed above, the relationships between the generalized parameters K and n and the model-specific flow parameters , n, and y are shown in Figure 1000-31. The equations for the Herschel-Bulkley Model become the equations for the Power Law Model when y = 0 and for the Bingham Model when n = 1. Finally, the equations in Figure 1000-31 for the Bingham and Herschel-Bulkley models involve the term
(Eq. 1000-22)
where: tw = Dp/4L, the shear stress at the wall This presents some computational difficulties in that p is needed to calculate p. In the computer program this is done in iterative fashion, a value for w is assumed and then w is calculated. This process is repeated until the assumed value equals the calculated value.
Chevron Corporation
1000-51
March 1997
Fig. 1000-31 Generalized Flow Parameters Constitutive Equation Power law = sn Generalized Flow Parameters
Bingham = y + s
Herschel-Bulkley = y + sn
March 1997
1000-52
Chevron Corporation
(Eq. 1000-23)
where: b = BAYS Designs usually include a safety factor over this estimate. This model assumes that all of the gelled oil in the section will yield instantaneously at the pipe wall as a rigid body. This is rarely so. For instance, thermal shrinkage may result in autodestruction of the gel. Voids may be present that permit some flow before general yielding. Temperature gradients across the pipe radius may leave a slower cooling, liquid core in large diameter pipelines. Finally, pressure waves may generate much higher local pressure gradients as the wave moves down the pipeline. None of several attempts to model these mitigating effects has been successful. Fortunately, Equation 1000-23 usually overpredicts the pressure gradient needed. The peak value mentioned above may be used for the BAYS. If the design is capable of generating shear stresses significantly larger than this peak value, then the pipeline can be considered restartable. If more accurate estimates are needed, the plot of b versus T, along with the estimated temperature distribution along the pipeline, can be used to calculate the required pressure gradient distribution. After the initial yielding, the waxy structure continues to break down and, for a constant applied pressure gradient, the flow rate continues to increase. Estimates can be made of this breakdown process but the laboratory work needed to generate the flow properties and the engineering work needed to do the calculations are very time consuming. These estimates are not done routinely.
Mitigating Actions
There are several mitigating actions available if the handling of waxy crudes becomes troublesome. Lines can be kept hot with insulation or heaters. Non-waxy crudes, lighter oils (diluents or solvents), and flow-improving additives may be added. Additives can be effective, but are very crude-specific.
Chevron Corporation
1000-53
March 1997
One should not generalize about the behavior of waxy crudes. Each seems to be different, if not idiosyncratic, and must be considered a unique material. The design features selected and mitigating actions used will, consequently, be site specific.
1040 References
1041 Viscosity Conversion
1. ODonnell, R. J. Equations for Converting Different Viscosity Units. American Society for Testing and Materials, Materials Research and Standards, Vol. 9, No. 5. May, 1969. Standard Method for Conversion of Kinematic Viscosity to Saybolt Universal Viscosity or to Saybolt Furol Viscosity. American Society for Testing and Materials, ASTM Designation D2161-66, ASTM Standards, Vol. 17, 1969.
2.
7. 8.
March 1997
1000-54
Chevron Corporation
Chevron Corporation
1100-1
March 1997
1110 Introduction
Recently, PC-based fluid flow programs have proliferated, and this section does not attempt to cover them all or limit the user to those that are covered. Each program discussed below tends to be strongest in a particular area and best suited for a certain type of calculation. Users are encouraged to find the most productive program for the calculations they are doing, even if it is not mentioned in this manual. In addition to the brief discussions contained in this section, more extended treatment is provided for PCFLOW, PIPEFLOW-2, HOTOIL, HOTOL*, SURGE, and PCSURGE in Appendices D, E, F, G, H, and I respectively.
PIPEFLOW-2
PIPEFLOW-2 is a full featured, steady state pressure and flow simulation program. It handles single and multiphase piping systems and networks. It will model complex systems of pipelines, production and injection wells (including reservoir inflow to and from well-boxes), and surface facilities such as compressors and separators. Pressures, temperatures, flow rates, and fluid compositions can be calculated at any point in a network. The program is recommended for problems involving a network, multiphase flow, flashing liquids, or vertical flow, or for compositional modeling during flow calculations. For details on PIPEFLOW-2 program operation see Appendix E.
HOTOIL
HOTOIL is a PC based program that handles the full range of non-Newtonian fluid flow models as well as Newtonian flow. HOTOIL is as rigorous as HOTOL* in its heat transfer and Newtonian fluid property correlations. Friction heating and pressure correction to viscosity can be turned on or off during a run. Elevation profiles are accommodated. HOTOIL is very flexible, allowing inputs to be defined in English units, SI units, or any combination thereof. HOTOIL lets the user specify either the upstream or the downstream pressure in either psi or head. See
March 1997
1100-2
Chevron Corporation
Appendix F for details on HOTOIL program operation and ways to obtain a copy of the program.
HOTOL*
HOTOL* is useful for calculating pressure drop and heat transfer in pipelines where fluid temperature changes cause significant changes in fluid properties. The programs fluid property correlations assume the fluid is liquid phase hydrocarbon. The heat transfer routines assume the fluid is at or above ambient temperature. HOTOL* does not consider friction heating or pressure correction to viscosity, nor does it accommodate elevation profiles. HOTOL* is available on the mainframe engineering disk. Details on the use of the program are given in Appendix G. See also Section , HOTPIPE2.
SURGE
The SURGE program models surges (water hammer) for systems of considerable complexity in which the user may specify valves, pumps, and surge chambers in any arrangement. It can accommodate multiple inlets, outlets, loops, and dead ends, and simulate the normal function of relief valves and check valves. Other types of valves with specified characteristics may be opened or closed at specified times. Pumps may be started or stopped at specified times, or tripped if calculated pressure transients exceed a specified value. Discharge-pressure-controlled relief to pump suction can be simulated, including off-control operation (spill back valve closed or wide open). Initial flow in various branches may be finite or zero. Steel pipe, plastic pipe, concrete lined pipe, hoses, and combinations thereof can be modeled. SURGE does not model pressure pulsations from reciprocating pumps. SURGE is in the Chevron Research and Technology Company computer program library. Gaining access to it is explained in Appendix H. (SURGE was previously available only on Chevrons VM mainframe computer.) The program is also available on floppy disks from Rob Hohmann, CRTC Richmond, CTN 242-2216 for use on IBM compatible PCs. See also Section , PCSURGE and Section , HYDRESS.
PCSURGE
PCSURGE is a menu driven input and output system built around the original SURGE program. PCSURGE will run on an IBM compatible PC with a math coprocessor. The math co-processor is necessary because the SURGE calculations are performed in FORTRAN. PCSURGE allows quick, interactive creation and manipulation of SURGE input data files. It allows a user to load an input data file, edit and save the new data file, run SURGE calculations, and view the output on the screen or send the output to a file. The menu system was designed to be self-explanatory. Most users who are familiar with the SURGE program should be able to use PCSURGE with minimal reference to SURGE or PCSURGE documentation.
Chevron Corporation
1100-3
March 1997
For all but the simplest systems, it is important that the user start with a good sketch of the system to be modeled, with all branches and nodes numbered. The sketch should include geometry, elevations, initial flow rates, pressures, etc. Users of PCSURGE should refer to Appendix I for an thorough explanation of the programs operation, capabilities, inputs, and outputs and for instructions on obtaining a copy of the program.
HOTPIPE2
HOTPIPE2 is a modification of HOTOL*. It retains the rigorous fluid property correlations and heat transfer routines of HOTOL*, but also accommodates elevation profiles. The user can place pump stations and heater stations along the length of the pipeline, or instruct the program to automatically place pumps or heaters when specified pressures and temperatures are reached. Friction heating and pressure correction to viscosity can be turned on and off during a run, allowing results to be compared to programs which do not have those features. HOTPIPE2 runs on an IBM compatible personal computer. It is available from Rob Hohmann, CRTC Richmond, CTN 242-2216.
CPAS
The Chevron Production Analysis Simulator (CPAS), is a general-purpose multiphase simulator designed to analyze producing systems from the outer boundary of a reservoir to the outlet of a producing system. It can analyze flows across completions, up tubing strings, and through gas-lift valves and submersible pumps. CPAS is available on Houston VM HOVMB. To run CPAS type CPAS at the VM prompt. Contact Rob Hohmann, CRTC Richmond, CTN 242-2216, for current CPAS support information. CPAS was designed primarily to help engineers in producing organizations to evaluate and locate deficiencies in oil and gas producing systems, and to design production systems. A nodal analysis approach similar to the one used for analyzing centrifugal pump performance is used to evaluate pressure and rate relationships in oil and gas producing systems. With CPAS nodal analysis, engineers can determine the causes of low flow rates, evaluate the relative magnitude of pressure drops through completion zones and strings, identify bottlenecks, decide if a well needs an artificial-lift system, and determine incremental production rates. CPAS has proven to be a valuable tool in the design of more efficient, cost-effective production systems. The simulator can provide additional information, enabling engineers to make more informed decisions, and thereby increase production. Producing systems which can be evaluated by CPAS are: Reservoir Perforations Tubing, with all restrictions and artificial lift equipment Flowlines with all production facilities
The performances of reservoir, completion, well, and surface devices can be combined to calculate the operating point for naturally flowing or artificial-lift
March 1997
1100-4
Chevron Corporation
wells. Also, the effects of changing tubing size, flow-line size, submersible pump horsepower, or gas-lift gas rate in a production system can be evaluated and operating points calculated. Analysis Options. The primary analysis options in CPAS are: Production system performance (for tubing, flowline, gas lift, and submersible pump) Reservoir vs. production system performance (Inflow performance relationship (IPR) vs. tubing-surface intake) Surface vs. subsurface production system performance (tubing performance vs. flowline intake) Well completion performance Gas lift design
The first four options can be used for oil and gas single or multiphase producing systems. Production System Performance. This option provides pressure-flow rate relationships for a given production system. The effects of changes in one parameter on the system can be evaluated. For example, the effect of changing the tubing diameter on the pressure-flow rate relationship can be shown. A production system, as defined by the simulator, consists of one or a combination of the following devices: Surface devices Pipe Pumps Risers
The behavior of the production system is determined by calculating pressure losses for each device based on multiphase fluid-flow techniques. The total pressure loss is calculated from the wellbore to the outlet of the production system. The effect of each device on total pressure losses in the system can be evaluated with this option.
Chevron Corporation
1100-5
March 1997
The following device parameters can be varied within the same run:
Device Tubing Pipe Pumps Gas-lift valves Parameter Tubing inner diameter Pipe inner diameter Pump horsepower Gas-lift rate
The simulator contains various fluid flow correlations and fluid property calculations. The pressure loss for each device is calculated based on user-specified correlations. Reservoir vs. Production System Performance. A producing system consists of a reservoir and a production system. In this option, the effect of changes in one parameter in the producing system can be evaluated to determine whether the system is limited by the reservoir or by the production system (from wellbore to the outlet point). Reservoir performance can be described by an inflow performance curve (IPR) which is constructed from a wells productively index (PI). The following methods are available to construct an IPR curve: For oil Straight-line PI method (Muskat) Vogel method Fetkovich radial method (solution gas and water-drive reservoirs) Well test method For gas J method for known C and n Multipoint backpressure test
This option is also used to evaluate the behavior of a device in a production system (such as tubing intake, flowline intake, pump horsepower, and gas-lift rate). An operating point is calculated for every reservoir and production system curve. Surface vs. Subsurface Production System Performance. This option divides the production system into surface and subsurface production systems and calculates whether the entire production system is limited at the surface or subsurface. For example, a flowline that is too large or tubing that is too small can retard the production rate. This option is effective for evaluating bottlenecks in existing production systems. Well Completion Performance. This option determines the pressure drop across the completion interval for gravel-packed and open-hole completions. The effects of shot density, perforation diameter, and gravel-pack permeability on the production system can be evaluated.
March 1997
1100-6
Chevron Corporation
Gas Lift Design. The CPAS gas-lift option applies a universal design and spacing approach for all types of continuous flow gas-lift valves. It is used to design gas-lift strings based on user specified desired production rate, wellhead pressure, bottomhole conditions, etc. It calculates valve locations by depth with corresponding casing pressure and generates a plot of pressure and gas-lift valve locations by depth. Primary calculations include point of injection, gas-lift valve depths, fluid gradient, and design tubing wellhead pressure. Data Entry and Retrieval. The simulator features user-friendly input menus and an interactive output graphics system, and can be accessed with one command. The menus help users prepare their input files and run the program. The simulator is easy to use and can be successfully run by users with minimal computer knowledge. Output consists of a variety of reports and plots. Standard output includes a onepage summary report and a plot of the results. The following plots represent the available graphics: Production system curves (such as tubing intake) IPR vs. production system curves Subsurface vs. surface device curves Completion performance curves
Each plot contains a summary of input and output data and can be displayed on a terminal or routed to a printer.
HYDRESS
HYDRESS is a modification of the SURGE program. It was designed to handle fluid transients (surge pressures) in long flexible conduits. Its specific use is to model the time it takes high pressure and low pressure waves to travel from offshore platforms to subsea controls through rubber hose and stainless steel hydraulic control lines. HYDRESS is a mainframe FORTRAN program which resides on OELIB (Offshore Engineering Library) on HOVMB. To access HYDRESS type OELIB at the VM prompt.
ORIFICE
ORIFICE is a stand-alone flow meter design and analysis program. It is available on the VM system in the CRTC Technical Library. The full description and operating manual for ORIFICE are contained in Appendix A of the Instrumentation and Control Manual. The orifice design calculations performed by the program are suitable for flow meters in process plants and producing facilities where the orifice meters do not require accuracy better than plus or minus 2 percent. It is not suitable for custody transfer applications, where flow rate calculations must conform to American Gas Association Report No. 3 (AGA-3) Orifice Metering of Natural Gas (ANSI/API 2530).
Chevron Corporation
1100-7
March 1997
ORIFICE performs three kinds of calculations, as follows: Sizing new orifices given standard differential pressure and full scale flow rate Reranging existing orifices given orifice diameter and full-scale flow rate Calculating flow rate given differential pressure and orifice diameter
Three kinds of orifice plates are included in the program, as follows: Square edge orifice plates for clean liquid, gas, and low-velocity vapor (steam) flows in 2-inch or larger pipe ASME small bore orifice plates (for liquid vapor or gas) usually with 1/2 to 11/2-inch (nominal) pipe size, depending on the type of taps used Quadrant edge orifice (or quarter circle) plates for viscous liquids only; often used where the pipe Reynolds number is low (usually below 10,000)
Conic edge and integral orifices are not covered in the ORIFICE program. The program covers four types of pressure taps: Flange taps. One inch upstream and 1 inch downstream of the orifice plate Radius taps. One D (pipe inside diameter) upstream and 1/2 D downstream Corner taps. Front and rear face of orifice plate Pipe taps. Two-and-one-half D upstream and 8 D downstream
ORIFICE handles single phase flow for three kinds of fluid flows and their engineering units: Liquid. Flow rate units can be in barrels per day (bpd), barrels per hour (bph), gallons per minute (gpm) or gallons per hour (gph), all corrected to 60F Vapor (steam). Flow rate unit should be in pounds per hour (pph) Gas. Flow rate units should be in standard cubic feet per hour (scfh), i.e., at 60F and 14.73 psia, which are the base conditions used by AGA-3
The calculations in ORIFICE are based on the latest applicable industry standards and technical reports which contain certain limits. These limits are generalized to be incorporated in the program. ORIFICE may be accessed on the VM system as follows: 1. 2. 3. 4. Log on to HOVMA, HOVMB, or HOVMC Type in ATT T9ENGR Type in ORIFICE when you see the first menu for program choice in the Technical Library (ORIFICE is listed under VM programs) Select screen version desired
March 1997
1100-8
Chevron Corporation
COMFLOW
COMFLOW solves gas pipeline network flow problems using compressible flow equations. The network may consist of any number of upstream locations feeding one downstream location. The program accounts for heat transfer effects. Compressors and let-down valves may be specified. The program is written in FORTRAN and runs on the VM mainframe. To obtain a copy of the program contact Rob Hohmann, CRTC Richmond, CTN 242-2216.
Chevron Corporation
1100-9
March 1997
Page A-2
Chevron Corporation
A-1
January 1990
Appendix A
A1.0
Conversion Factors
The table of conversion factors is based on the Gould Pump Manual table, with some additions. They are not meant to include all possible units or conversions. The exact conversion factors may in some cases disagree with other authorities. English measuresunless otherwise designated, are those used in the United States, and the units of weight and mass are avoirdupois units. Gallondesignates the U.S. gallon. To convert into the imperial gallon, multiply the U.S. gallon by 0.83267. Likewise, the word ton designates a short ton, 2,000 pounds. Properties of waterit freezes at 32F, and is at its maximum density at 39.2F. In the multipliers using the properties of water, calculations are based on water at 39.2F in a vacuum, weighing 62.427 pounds per cubic foot, or 8.345 pounds per U.S. gallon.
A1.1
January 1990
A-2
Chevron Corporation
Appendix A
Multiply Acres Acres Acres Acres Acre-feet Acre-feet Acre-feet Atmospheres Atmospheres Atmospheres Atmospheres Atmospheres Atmospheres
By 43,560 4047 1.562 x 103 4840 43,560 325,851 1233.48 76.0 29.92 33.90 10.332 14.70 1.058
To Obtain Square feet Square meters Square miles Square yards Cubic feet Gallons Cubic meters Cm of mercury Inches of mercury Feet of water Kg/sq meter Lb/sq inch Tons/sq ft
Barrels-oil Barrels-Beer Barrels-Whiskey Barrels/Day-oil Bags or sacks-cement Board feet British thermal units British thermal units British thermal units British thermal units British thermal units British thermal units Btu/min Btu/min Btu/min Btu/min
42 31 45 0.02917 94 144 sq in. x 1 in. 0.2520 777.6 3.927 x 104 107.5 2.928 x 104 1.0551 x 103 12.96 0.02356 0.01757 17.57
Gallons-oil Gallons-Beer Gallons-Whiskey Gallons/Min-oil Pounds-cement Cubic inches Kilogram-calories Foot-lb Horsepower-hr Kilogram-meters Kilowatt-hr Joules Foot-lb/sec Horsepower Kilowatts Watts
1 0.01 0.01
Chevron Corporation
A-3
January 1990
Appendix A
Multiply Centimeters Centimeters Centimeters Centimeters of Mercury Centimeters of Mercury Centimeters of Mercury Centimeters of Mercury Centimeters of Mercury Centimeters/sec Centimeters/sec Centimeters/sec Centimeters/sec Centimeters/sec Centimeters/sec Cm/sec/sec Cubic centimeters Cubic centimeters Cubic centimeters Cubic centimeters Cubic centimeters Cubic centimeters Cubic centimeters Cubic centimeters Cubic feet Cubic feet Cubic feet Cubic feet Cubic feet Cubic feet Cubic feet Cubic feet Cubic feet/min Cubic feet/min Cubic feet/min
By 0.3937 0.01 10 0.01316 0.4461 136.0 27.85 0.1934 1.969 0.03281 0.036 0.6 0.02237 3.728 x 10-4 0.03281 3.531 x 10-5 6.102 x 10-2 10-4 1.308 x 10-4 2.642 x 10-4 9.999 x 10-4 2.113 x 10-3 1.057 x 10-3 2.832 x 10-4 1728 0.02832 0.03704 7.48052 28.32 59.84 29.92 472.0 0.1247 0.4719
To Obtain Inches Meters Millimeters Atmospheres Feet of water Kg/sq meter Lb/sq ft Lb/sq inch Feet/min Feet/sec Kilometers/hr Meters/min Miles/hr Miles/min Feet/sec/sec Cubic feet Cubic inches Cubic meters Cubic yards Gallons Liters Pints (liq) Quarts (liq) Cubic cm Cubic inches Cubic meters Cubic yards Gallons Liters Pints (liq) Quarts (liq) Cubic cm/sec Gallons/sec Liters/sec
January 1990
A-4
Chevron Corporation
Appendix A
Multiply Cubic feet/min Cubic feet/sec Cubic feet/sec Cubic inches Cubic inches Cubic inches Cubic inches Cubic inches Cubic inches Cubic inches Cubic inches Cubic meters Cubic meters Cubic meters Cubic meters Cubic meters Cubic meters Cubic meters Cubic meters Cubic meters/hr Cubic yards Cubic yards Cubic yards Cubic yards Cubic yards Cubic yards Cubic yards Cubic yards Cubic yards/min Cubic yards/min Cubic yards/min
By 62.43 0.646317 448.831 16.39 5.787 x 10-4 1.639 x 10-5 2.143 x 10-5 4.329 x 10-3 1.639 x 10-2 0.03463 0.01732 106 35.31 61023 1.308 264.2 999.97 2113 1057 4.40 764,554.86 27 46.656 0.7646 202.0 764.5 1616 807.9 0.45 3.366 12.74
To Obtain Pounds of water/min Millions gals/day Gallons/min Cubic centimeters Cubic feet Cubic meters Cubic yards Gallons Liters Pints (liq) Quarts (liq) Cubic centimeters Cubic feet Cubic inches Cubic yards Gallons Liters Pints (liq) Quarts (liq) Gallons/min Cubic centimeters Cubic feet Cubic inches Cubic meters Gallons Liters Pints (liq) Quarts (liq) Cubic feet/sec Gallons/sec Liters/sec
Decigrams Deciliters
0.1 0.1
Grams Liters
Chevron Corporation
A-5
January 1990
Appendix A
Multiply Decimeters Degrees (angle) Degrees (angle) Degrees (angle) Degrees/sec Degrees/sec Degrees/sec Dekagrams Dekaliters Dekameters Drams Drams Drams Dynes
By 0.1 60 0.01745 3600 0.01745 0.1667 0.002778 10 10 10 27.34375 0.0625 1.771845 1 x 105
To Obtain Meters Minutes Radians Seconds Radians/sec Revolutions/min Revolutions/sec Grams Liters Meters Grains Ounces Grams Newtons
Fathoms Feet Feet Feet Feet Feet of water Feet of water Feet of water Feet of water Feet of water Feet/min Feet/min Feet/min Feet/min Feet/min Feet/sec Feet/sec Feet/sec Feet/sec
6 30.48 12 0.3048 1/3 0.0295 0.8826 304.8 62.43 0.4335 0.5080 0.01667 0.01829 0.3048 0.01136 30.48 1.097 0.5924 18.29
Feet Centimeters Inches Meters Yards Atmospheres Inches of mercury Kg/sq meter Lb/sq ft Lb/sq inch Centimeters/sec Feet/sec Kilometers/hr Meters/min Miles/hr Centimeters/sec Kilometers/hr Knots Meters/min
January 1990
A-6
Chevron Corporation
Appendix A
Multiply Feet/sec Feet/sec Feet/sec/sec Feet/sec/sec Foot-pounds Foot-pounds Foot-pounds Foot-pounds Foot-pounds Foot-pounds/min Foot-pounds/min Foot-pounds/min Foot-pounds/min Foot-pounds/min Foot-pounds/sec Foot-pounds/sec Foot-pounds/sec Foot-pounds/sec
By 0.6818 0.01136 30.48 0.3048 1.286 x 10-3 5.050 x 10-7 3.240 x 10-4 0.1383 3.766 x 10-7 2.140 x 10-5 0.01667 3.030 x 10-5 5.393 x 10-3 2.280 x 10-5 7.704 x 10-2 1.818 x 10-3 1.941 x 10-2 1.356 x 10-3
To Obtain Miles/hr Miles/min Cm/sec/sec Meters/sec/sec British thermal units Horsepower-hr Kilogram-calories Kilogram-meters Kilowatt-hours Btu/sec Foot-pounds/sec Horsepower Gm-calories/sec Kilowatts Btu/min Horsepower Kg-calories/min Kilowatts
Gallons Gallons Gallons Gallons Gallons Gallons Gallons Gallons Gallons-Imperial Gallons-U.S. Gallons water Gallons/min Gallons/min Gallons/min Grains (troy)
3785 0.1337 231 3.785 x 10-3 4.951 x 10-3 3.785 8 4 1.20095 0.83267 8.345 2.228 x 10-3 0.06308 8.0208 0.06480
Cubic centimeters Cubic feet Cubic inches Cubic meters Cubic yards Liters Pints (liq) Quarts (liq) U.S. gallons Imperial Gallons Pounds of water Cubic feet/sec Liters/sec Cu ft/hr Grams
Chevron Corporation
A-7
January 1990
Appendix A
Multiply Grains (troy) Grains (troy) Grains/U.S. gal Grains/U.S. gal Grains/Imp gal Grams Grams Grams Grams Grams Grams Grams Grams/cm Grams/cu cm Grams/cu cm Grams/liter Grams/liter Grams/liter Grams/liter
By 0.04167 2.0833 x 10-3 17.118 142.86 14.254 980.7 15.43 .001 1000 0.03527 0.03215 2.205 x 10-3 5.600 x 10-3 62.43 0.03613 58.416 8.345 0.06242 1000
To Obtain Pennyweights (troy) Ounces (troy) Parts/million Lb/million gal Parts/million Dynes Grains Kilograms Milligrams Ounces Ounces (troy) Pounds Pounds/inch Pounds/cubic foot Pounds/cubic inch Grains/gal Pounds/1000 gals Pounds/cubic foot Parts/million
Hectares Hectares Hectograms Hectoliters Hectometers Hectowatts Horsepower Horsepower Horsepower Horsepower Horsepower Horsepower Horsepower Horsepower (boiler)
2.471 1.076 x 105 100 100 100 100 42.44 33,000 550 1.014 10.547 0.7457 745.7 33,493
Acres Square feet Grams Liters Meters Watts Btu/min Foot-lb/min Foot-lb/sec Horsepower (metric) Kg-calories/min Kilowatts Watts Btu/hr
January 1990
A-8
Chevron Corporation
Appendix A
Inches Inches of mercury Inches of mercury Inches of mercury Inches of mercury Inches of mercury (32F) Inches of water Inches of water Inches of water Inches of water Inches of water Inches of water
2.540 0.03342 1.133 345.3 70.73 0.491 0.002458 0.07355 25.40 0.578 5.202 0.03613 0.9478 x 10-3
Centimeters Atmospheres Feet of water Kg/sq meter Lb/sq foot Lb/sq inch Atmospheres Inches of mercury Kg/sq meter Ounces/sq inch Lb/sq foot Lb/sq inch
Joules
Btu
Kilograms Kilograms Kilograms Kilograms Kilograms (force) Kilograms-cal/sec Kilograms-cal/sec Kilograms-cal/sec Kilograms-cal/sec Kilogram-cal/min Kilogram-cal/min Kilogram-cal/min
980.665 2.205 1.102 x 10-3 103 9.80665 3.968 3086 5.6145 4186.7 3085.9 0.09351 69.733
Dynes Pounds Tons (short) Grams Newtons Btu/sec Foot-lb/sec Horsepower Watts Foot-lb/min Horsepower Watts
Chevron Corporation
A-9
January 1990
Appendix A
Multiply Kg/meter Kg/sq meter Kg/sq meter Kg/sq meter Kg/sq meter Kg/sq meter Kg/sq millimeter Kiloliters Kilometers Kilometers Kilometers Kilometers Kilometers Kilometers/hr Kilometers/hr Kilometers/hr Kilometers/hr Kilometers/hr Kilometers/hr Km/hr/sec Km/hr/sec Km/hr/sec Kilopascals Kilowatts Kilowatts Kilowatts Kilowatts Kilowatts Kilowatts Kilowatt-hours Kilowatt-hours Kilowatt-hours Kilowatt-hours Kilowatt-hours
By 0.6720 9.678 x 10-5 3.281 x 10-3 2.896 x 10-3 0.2048 1.422 x 10-3 106 103 105 3281 103 0.6214 1094 27.78 54.68 0.9113 .5399 16.67 0.6214 27.78 0.9113 0.2778 6.895 56.907 4.425 x 104 737.6 1.341 14.34 103 3414.4 2.655 x 106 1.341 860.4 3.671 x 105
To Obtain Lb/foot Atmospheres Feet of water Inches of mercury Lb/sq foot Lb/sq inch Kg/sq meter Liters Centimeters Feet Meters Miles Yards Centimeters/sec Feet/min Feet/sec Knots Meters/min Miles/hr Cm/sec/sec Ft/sec/sec Meters/sec/sec Pounds/sq inch Btu/min Foot-lb/min Foot-lb/sec Horsepower Kg-calories/min Watts Btu Foot-lb Horsepower-hr Kilogram-calories Kilogram-meters
January 1990
A-10
Chevron Corporation
Appendix A
Multiply
By 103 0.03531 61.02 10-3 1.308 x 10-3 0.2642 2.113 1.057 5.886 x 10-4 4.403 x 10-3 Length (ft.)
To Obtain
Liters Liters Liters Liters Liters Liters Liters Liters Liters/min Liters/min Lumber Width (in.) x Thickness (in.)/12
Cubic centimeters Cubic feet Cubic inches Cubic meters Cubic yards Gallons Pints (liq) Quarts (liq) Cubic ft/sec Gal/sec Board feet
Meters Meters Meters Meters Meters Meters Meters/min Meters/min Meters/min Meters/min Meters/min Meters/sec Meters/sec Meters/sec Meters/sec Meters/sec Meters/sec Microns Miles Miles
100 3.281 39.37 10-3 103 1.094 1.667 3.281 0.05468 0.06 0.03728 196.8 3.281 3.6 0.06 2.287 0.03728 10-6 1.609 x 105 5280
Centimeters Feet Inches Kilometers Millimeters Yards Centimeters/sec Feet/min Feet/sec Kilometers/hr Miles/hr Feet/min Feet/sec Kilometers/hr Kilometers/min Miles/hr Miles/min Meters Centimeters Feet
Chevron Corporation
A-11
January 1990
Appendix A
Multiply Miles Miles Miles/hr Miles/hr Miles/hr Miles/hr Miles/hr Miles/hr Miles/min Miles/min Miles/min Miles/min Milliers Milligrams Milliliters Millimeters Millimeters Milligrams/liter Million gals/day Miners inches Minutes (angle)
By 1.609 1760 44.70 88 1.467 1.609 0.8689 26.82 2682 88 1.609 60 103 10-3 10-3 0.1 0.03937 1 1.54723 1.5 2.909 x 10-4 1 x 10-5 0.22481 0.10197
To Obtain Kilometers Yards Centimeters/sec Feet/min Feet/sec Kilometers/hr Knots Meters/min Centimeters/sec Feet/sec Kilometers/min Miles/hr Kilograms Grams Liters Centimeters Inches Parts/million Cubic ft/sec Cubic ft/min Radians
Drams Grains Pounds Grams Ounces (troy) Tons (long) Tons (metric) Grains
January 1990
A-12
Chevron Corporation
Appendix A
Multiply Ounces (troy) Ounces (troy) Ounces (troy) Ounces (troy) Ounces (fluid) Ounces (fluid) Ounces/sq inch
To Obtain Pennyweights (troy) Pounds (troy) Grams Ounces (avoir) Cubic inches Liters Lb/sq inch
Parts/million Parts/million Parts/million Pennyweights (troy) Pennyweights (troy) Pennyweights (troy) Pennyweights (troy) Pounds Pounds Pounds Pounds Pounds Pounds Pounds Pounds (troy) Pounds (troy) Pounds (troy) Pounds (troy) Pounds (troy) Pounds (troy) Pounds (troy) Pounds (troy) Pounds (troy) Pounds of water Pounds of water Pounds of water
0.0584 0.07015 8.345 24 1.55517 0.05 4.1667 x 10-3 16 256 7000 0.0005 453.5924 1.21528 14.5833 5760 240 12 373.2417 0.822857 13.1657 3.6735 x 10-4 4.1143 x 10-4 3.7324 x 10-4 0.01602 27.68 0.1198
Grains/U.S. gal Grains/Imp gal Lb/million gal Grains Grams Ounces (troy) Pounds (troy) Ounces Drams Grains Tons (short) Grams Pounds (troy) Ounces (troy) Grains Pennyweights (troy) Ounces (troy) Grams Pounds (avoir) Ounces (avoir) Tons (long) Tons (short) Tons (metric) Cubic feet Cubic inches Gallons
Chevron Corporation
A-13
January 1990
Appendix A
Multiply Pounds of water/min Pounds/cubic foot Pounds/cubic foot Pounds/cubic foot Pounds/cubic inch Pounds/cubic inch Pounds/cubic inch Pounds/foot Pounds/inch Pounds (force) Pounds/sq foot Pounds/sq foot Pounds/sq foot Pounds/sq inch Pounds/sq inch Pounds/sq inch Pounds/sq inch Pounds/sq inch
By 2.670 x 10-4 0.01602 16.02 5.787 x 10-4 27.68 2.768 x 104 1728 1.488 1152 4.44822 0.01602 4.882 6.944 x 10-3 0.06804 2.307 2.036 703.1 0.146
To Obtain Cubic ft/sec Grams/cubic cm Kg/cubic meters Lb/cubic inch Grams/cubic cm Kg/cubic meter Lb/cubic foot Kg/meter Grams/cm Newtons Feet of water Kg/sq meter Pounds/sq inch Atmospheres Feet of water Inches of mercury Kg/sq meter Kilopascals
Quadrants (angle) Quadrants (angle) Quadrants (angle) Quarts (dry) Quarts (liq) Quintal, Argentine Quintal, Brazil Quintal, Castile, Peru Quintal, Chile Quintal, Mexico Quintal, metric Quires
90 5400 1.571 67.20 57.75 101.28 129.54 101.43 101.41 101.47 220.46 25
Degrees Minutes Radians Cubic inches Cubic inches Pounds Pounds Pounds Pounds Pounds Pounds Sheets
Radians Radians
57.30 3438
Degrees Minutes
January 1990
A-14
Chevron Corporation
Appendix A
Multiply Radians Radians/sec Radians/sec Radians/sec Radians/sec/sec Radians/sec/sec Reams Revolutions Revolutions Revolutions Revolutions/min Revolutions/min Revolutions/min Revolutions/min/min Revolutions/min/min Revolutions/sec Revolutions/sec Revolutions/sec Revolutions/sec/sec Revolutions/sec/sec
By 0.637 57.30 0.1592 9.549 573.0 0.1592 500 360 4 6.283 6 0.1047 0.01667 1.745 x 10-3 2.778 x 10-4 360 6.283 60 6.283 3600 4.848 x 10-4 1.076 x 10-3 0.1550 10-4 100 2.296 x 10-5 929.0 144 0.09290 3.587 x 10-4 1/9 8.0208 6.452
To Obtain Quadrants Degrees/sec Revolutions/sec Revolutions/min Rev/min/min Rev/sec/sec Sheets Degrees Quadrants Radians Degrees/sec Radians/sec Revolutions/sec Rads/sec/sec Revs/sec/sec Degrees/sec Radians/sec Revolutions/min Radians/sec/sec Revs/min/min
Seconds (angle) Square centimeters Square centimeters Square centimeters Square centimeters Square feet Square feet Square feet Square feet Square feet Square feet 1/Sq ft/gal/min Square inches
Radians Square feet Square inches Square meters Square millimeters Acres Square centimeters Square inches Square meters Square miles Square yards Overflow rate (ft/hr) Square centimeters
Chevron Corporation
A-15
January 1990
Appendix A
Multiply Square inches Square inches Square kilometers Square kilometers Square kilometers Square kilometers Square kilometers Square meters Square meters Square meters Square meters Square miles Square miles Square miles Square miles Square millimeters Square millimeters Square yards Square yards Square yards Square yards
By 6.944 x 10-3 645.2 247.1 10.76 x 106 106 0.3861 1.196 x 106 2.471 x 10-4 10.76 3.861 x 10-7 1.196 640 27.88 x 106 2.590 3.098 x 106 0.01 1.550 x 10-3 2.066 x 10-4 9 0.8361 3.228 x 10-7
To Obtain Square feet Square millimeters Acres Square feet Square meters Square miles Square yards Acres Square feet Square miles Square yards Acres Square feet Square kilometers Square yards Square centimeters Square inches Acres Square feet Square meters Square miles
Temp (C) + 273 Temp (C) + 17.78 Temp (F) + 460 Temp (F) - 32 Tons (long) Tons (long) Tons (long) Tons (metric) Tons (metric) Tons (short) Tons (short) Tons (short)
1 1.8 1 5/9 1016 2240 1.12000 103 2205 2000 32,000 907.1843
Abs. temp (C) Temp (F) Abs. temp (F) Temp (C) Kilograms Pounds Tons (short) Kilograms Pounds Pounds Ounces Kilograms
January 1990
A-16
Chevron Corporation
Appendix A
Multiply Tons (short) Tons (short) Tons (short) Tons (short) Tons of water/24 hrs Tons of water/24 hrs Tons of water/24 hrs
To Obtain Pounds (troy) Tons (long) Ounces (troy) Tons (metric) Pounds water/hr Gallons/min Cu ft/hr
Watts Watts Watts Watts Watts Watts Watt-hours Watt-hours Watt-hours Watt-hours Watt-hours Watt-hours
0.05686 44.25 0.7376 1.341 x 10-3 0.01434 10-3 3.414 2655 1.341 x 10-3 0.8604 367.1 10-3
Btu/min Foot-lb/min Foot-lb/sec Horsepower Kg-calories/min Kilowatts Btu Foot-lb Horsepower-hours Kilogram-calories Kilogram-meters Kilowatt-hours
91.44 3 36 0.9144
Chevron Corporation
A-17
January 1990
Appendix A
January 1990
A-18
Chevron Corporation
Chevron Corporation
B-1
January 1990
Appendix B
January 1990
B-2
Chevron Corporation
Appendix B
Chevron Corporation
B-3
January 1990
Appendix B
January 1990
B-4
Chevron Corporation
Appendix B
Chevron Corporation
B-5
January 1990
Bulletin TT 330
Chevron Corporation
C-1
January 1990
Appendix C
January 1990
C-2
Chevron Corporation
Appendix C
Chevron Corporation
C-3
January 1990
Appendix C
January 1990
C-4
Chevron Corporation
Appendix C
Chevron Corporation
C-5
January 1990
Appendix C
January 1990
C-6
Chevron Corporation
Appendix C
Chevron Corporation
C-7
January 1990
Appendix C
January 1990
C-8
Chevron Corporation
Appendix C
Chevron Corporation
C-9
January 1990
Appendix C
January 1990
C-10
Chevron Corporation
Appendix C
Chevron Corporation
C-11
January 1990
Appendix C
January 1990
C-12
Chevron Corporation
Appendix C
Chevron Corporation
C-13
January 1990
Abstract
This section discusses the computer program PCFLOW. The program automates many of the equations of fluid flow found in this manual. Contents D1.0 D2.0 D3.0 D4.0 D5.0 D5.1 D5.2 D6.0 D7.0 Description Of PCFLOW Programming Philosophy PCFLOW Limitations Accessing PCFLOW Installing PCFLOW on Your Computer Installing PCFLOW from the Technical Standards Website Installing PCFLOW from a Floppy Disk Running PCFLOW Program Support D-4 D-4 Page D-2 D-2 D-2 D-3 D-3
Chevron Corporation
D-1
March 1997
Appendix D
D1.0
Description Of PCFLOW
PCFLOW is a Windows program that runs in both Windows 3.1 and Windows 95. This new version of PCFLOW was created by converting the original BASIC language menu-driven program to an easy to use Visual BASIC compiled program. The PCFLOW program files include executable files, sample data, text files, and message files. PCFLOW automates the equations presented in the following sections of this manual: Section 410, Incompressible Flow Section 420, Two-phase Flow Section 430, Compressible Flow Section 440, Gas Flow At High Pressure Drop (Choked Flow) Section 500, Fitting Pressure Drop Section 600, Noncircular Conduits Section 820, Maximum Surge Pressure in a Simple Case
D2.0
Programming Philosophy
PCFLOW is easy to use. Its menu-driven screens are self-explanatory, and the equations used for calculations are the same as the ones presented in this manual. There is no user manual for PCFLOW. The user is as responsible for the input values when using PCFLOW as when using a calculator. No data checking is done by the program to constrain the users inputs to preconceived values. The only limit on user inputs is the inherent stability of the equations themselves in the Visual BASIC programming environment.
D3.0
PCFLOW Limitations
PCFLOW is not intended to replace, compete with, or discourage the use of commercially available piping design and fluid flow programs, which are often more useful in specific design cases. PCFLOW is not an all-purpose piping design program; it simply automates the equations for common fluid flow problems. PCFLOW is not designed to be a data base for fluid properties or piping data. Although PCFLOW provides typical input values in the help-screens and sample data files, the user should refer to this manual or other sources for input data. Fluid properties in the sample data files are from the VM mainframe program PPROP, which is also available in PC version.
March 1997
D-2
Chevron Corporation
Appendix D
D4.0
Accessing PCFLOW
PCFLOW is available three ways: 1. By accessing the Technical Standards Web site at http://techstds.rrc.chevron.com/tech_standards/ and downloading the program to your hard disk per the instructions given. By calling, e-mailing, or faxing Technical Standards customer service to request a disk from which the program can be installed to your hard disk. E-MAIL: TECHSTDS (use global address list) CTN: 242-7241 Fax CTN: 242-2157 By purchasing the Chevron Process Simulation Suite, a CD-ROM available from the Process Simulation Team at CRTC.
2.
3.
D5.0
6.
D5.2
Chevron Corporation
D-3
March 1997
Appendix D
D6.0
Running PCFLOW
To run PCFLOW, click the icon or the file called pcflow1.exe. A menu will appear on the screen and allow you to select the type of calculations you wish to perform. From this point, the program should be self-explanatory. Questions should be directed to the PCFLOW support contact (see Section D7.0) The following suggestions will help you make efficient use of PCFLOW: Modify sample data files rather than building input data from scratch. Use the HELP button to access help-screens that will assist in defining input data. Name your data files something other than DATA1 through DATA20 to preserve the sample data filenames.
D7.0
Program Support
If you have trouble running PCFLOW on your system or have questions about the program, please contact Yen Chen, CRTC, Richmond, CA, CTN-242-3437. If you have questions regarding the technical aspects of the program (i.e., correlations used, etc.), contact the Fluid Flow Manual sponsor Rob Hohmann, CTN 242-2216. Comments, questions, and suggestions concerning PCFLOW are welcome.
March 1997
D-4
Chevron Corporation
Contents E1.0 E2.0 E3.0 E4.0 E5.0 E6.0 E7.0 E8.0 E9.0 E10.0 E11.0 E12.0 E13.0 Introduction Data Preparation Using the Input Menu System Fluid Types Calculation Options General Input Requirements Pressure Loss Correlations Friction Calculations Flow Devices Fluid Property Correlations PVT Tables Output Reports Calibration Using PIPEFLOW-2
Page E-2 E-3 E-3 E-5 E-5 E-6 E-7 E-8 E-8 E-9 E-9 E-9 E-10
Chevron Corporation
E-1
October 1992
Appendix E
E1.0
Introduction
PIPEFLOW-2 is available on Houston VM HOVMB. Contact Rob Hohmann, CRTC Richmond, CTN 242-2216 for current support information. PIPEFLOW-2 includes a generalized, steady state pressure and flow simulator for single phase or multiphase fluid piping systems, and a computerized, mathematical description of fluid flow laws for producing and injection systems. Pressure and temperature losses and flow rates can be calculated at any point in the network for any combination of reservoir deliverability, fluid type, water cut and gas-oil ratio (GOR). PIPEFLOW-2 can also calculate two-phase steam flow behavior. The following discussion is only a brief introduction to the programs capabilities. A complete three-volume user manual is available from Chevron Exploration and Production Services Company. A multiday hands-on course is also available. PIPEFLOW-2 can handle problems ranging from a single pipe link to a highly complex network with hundreds of junction points. It can also link dynamically with Chevrons CRS-3D numerical reservoir simulator to simulate an entire oilfield system. Figure E-1 shows facilities that can be modeled using PIPEFLOW-2.
Fig. E-1
Optional output reports vary from concise summaries to detailed pressure profiles. A plotting postprocessor can display any calculated result. Most widely used pressure-flow correlations and fluid property techniques are included. New ones can easily be added. The program is dynamically dimensioned. It adjusts itself to the smallest possible increment of computer memory necessary for the given problem. This means that
October 1992
E-2
Chevron Corporation
Appendix E
the user with a small problem need not be aware that he is using a large, sophisticated system or pay a premium for more core than is necessary. It also means that problem sizes are limited only by available core.
E2.0
Data Preparation
The input menu system is recommended. It uses interactive, full screen forms to guide the user through problem setup. The system then automatically creates an input file and stores it on the users permanent disk. Each menu screen has a help function. Alternatively, the input forms guide the user in manually creating a dataset. Alternative approaches to preparing input to PIPEFLOW-2 are illustrated in Figure E-2.
Fig. E-2
E3.0
Chevron Corporation
E-3
October 1992
Appendix E
Fig. E-3
input menu system on the Houston VM, enter the command PIPEFLOW from the VM/CMS environment on HOVMB. After this command is entered a system identification logo appears, followed by the primary option menu. You may then create an input file or use any other feature (see Figure E-4).
Fig. E-4 Primary Options
October 1992
E-4
Chevron Corporation
Appendix E
E4.0
Fluid Types
Any fluid type currently available in PIPEFLOW-2 can be represented by the input menu system. These include the following: Black oil (multiphase oil, water, and gas) Condensate (multiphase condensate, water, and gas) Single phase gas Single phase liquid Steam (multiphase water and water vapor) Compositionally modeled fluids
E5.0
Calculation Options
The system has two primary options, single flow stream calculations and network calculations. A single flow stream is a combination of flow devices connected in a series. This combination makes up a flow link. The single stream option handles only one link at a time. A network consists of any branched or parallel combination of flow links.
Network Calculations
PIPEFLOW-2 network options solve for pressure and flow rates throughout a producing (or injecting) network, including both surface and wellbore flow devices. Some typical network applications are: Modeling an existing gathering system to help locate system bottlenecks Determining the effect of separator conditions on system capacity Determining the overall effect of planned system modifications Improving the design of new systems by rapidly evaluating responses under various conditions
Chevron Corporation
E-5
October 1992
Appendix E
Determining the combined effect of wellbore and surface facilities on reservoir performance using the PIPEFLOW-2/CRS interface Allocating known total system production to individual wells
Two network solution methods are available, the Newton-Raphson technique and a successive approximation method called DYNER (Dynamic Network Reduction). The Newton-Raphson method is normally used for networks with complex loops and crossovers. It solves a set of nonlinear differential equations that represents balanced flow conditions throughout the system using Newtons method of successive linear approximations. The DYNER network solution method is more straightforward (and generally faster) than Newton-Raphson. However, it is not capable of handling crossovers or loops and is best suited to tree-branched networks. DYNER calculates flow through each link on successive forward passes and pressure loss on corresponding backward passes through the network until the entire system is balanced.
Flow Tables
Solution time for multiphase network problems can often be reduced by first generating a set of pressure-flow tables for all flow links as functions of varying gas and water ratios. PIPEFLOW-2 can do this automatically. The program then interpolates among these tables to get finite-difference pressure-flow derivatives to use in balancing the network. The flow tables can be saved for additional runs, a significant saving on large studies.
E6.0
Since PIPEFLOW-2 network solution techniques are built around the concept of nodes and links, it will be helpful to define them. Node. A reference point for system pressures and fluid properties where links meet, begin, or end, or where flow can enter or leave the system
October 1992
E-6
Chevron Corporation
Appendix E
Link. A continuous-flow conduit between nodes. A link can be made up of any number of flow devices Device. An individual component of a link. A flow device is a segment of tubing or pipe of constant diameter. Several such devices might make up a tubing string for a well
Each link has a from and to node. PIPEFLOW-2 assumes the flow direction through the link to be positive from the from node to the to node. This, however, does not restrict the final solution of flow to that direction. It is helpful to think of nodes where flow enters or leaves the system as boundary nodes, and others as interior nodes. Node flow at interior nodes is zero by definition; node flow refers only to flow entering or leaving the system at that point and could be more precisely termed flow imbalance.
E7.0
Multiphase pressure drop calculations are very complex. They require prediction of liquid holdup and flow regime using complicated correlations and iterative calculations that move through each device in small increments. The most widely used pressure-loss-versus-flow-rate relationships are available in PIPEFLOW-2, as follows: Multiphase Vertical Orkiszewski Duns-Ros Hagedorn-Brown Beggs-Brill Beggs-Brill with Moody friction (recommended for initial trial run) Beggs-Brill no-slip Gray (gas condensate) Gray with Moody friction Aziz (friction)-Flanigan Angel-Welchon-Ross (water and free gas)
Chevron Corporation
E-7
October 1992
Appendix E
Multiphase chokes Fortunati (subcritical flow) Ashford (critical flow) Multiphase Horizontal Eaton Eaton-Flanigan Lockhart-Martinelli Beggs-Brill Beggs-Brill with Moody friction (recommended for initial trial run) Beggs-Brill no-slip Dukler Dukler-Flanigan Eaton (holdup)-Dukler
E8.0
Friction Calculations
Basic friction factor components for liquids and gases are calculated using a modified Wood-Colebrook iterative solution of the Moody diagram. The various pressure loss correlations combine the components in different ways. Some correlations have their own built-in friction calculations. For those that do not, PIPEFLOW-2 provides both the standard version and a modified version which substitutes the Moody-derived friction factors.
E9.0
Flow Devices
PIPEFLOW-2 can model the following devices: Well tubing Well annulus Surface pipelines Compressors Pumps Risers Chokes Well productivity and injectivity functions Gas lift valves Regulators User-defined P versus Q functions Check valves
October 1992
E-8
Chevron Corporation
Appendix E
Steam properties are calculated by subroutines that duplicate the Keenan and Keyes and ASME steam tables.
Chevron Corporation
E-9
October 1992
Appendix E
data. Examples are pressure and temperature profiles, gas lift analysis results, flow table generation results, and PVT property reports.
October 1992
E-10
Chevron Corporation
Abstract
This section discusses the computer program HOTOIL. The program calculates pressure drop and heat transfer in hydrocarbon pipelines and pipes where fluid temperature changes cause significant changes in fluid properties. HOTOIL is one of the few tools available that can perform pressure drop calculations for non-Newtonian fluids such as waxy crudes. Contents F1.0 F2.0 F3.0 F4.0 F5.0 F5.1 F5.2 F6.0 F7.0 F8.0 Description Of HOTOIL Programming Philosophy HOTOIL Limitations Accessing HOTOIL Installing HOTOIL on Your Computer Installing HOTOIL from the Technical Standards Website Running HOTOIL from a Floppy Disk Running HOTOIL Sample Input and Output Program Support F-4 F-4 F-5 Page F-2 F-2 F-2 F-3 F-3
Chevron Corporation
F-1
March 1997
Appendix F
F1.0
Description Of HOTOIL
HOTOIL is a compiled, Microsoft Quick Basic language program. It should run on any IBM-compatible personal computer. The HOTOIL files include the main program file, HOTOIL.EXE, and the following three sample data files: SAMPLEN.HOD - Sample data file for Newtonian fluid. SAMPLEP.HOD - Sample data file for the Power model of a non-Newtonian fluid. SAMPLEB.HOB - Sample data file for the Bingham model of a non-Newtonian fluid.
F2.0
Programming Philosophy
HOTOILs user interface is designed to be self explanatory and easy to use. This users guide is provided for completeness, but the instructions on the HOTOIL screen should be sufficient to lead users through the steps of creating a data file and running the HOTOIL calculations for most pipeline systems. HOTOIL uses the equations presented in Section 400, Friction Pressure Drop, of this manual for normal (Newtonian) fluids. The equations used by HOTOIL for modeling non-Newtonian fluids are presented in Section 1030 of this manual. HOTOIL breaks a pipeline into short calculation segments, each with a temperature drop of no more than 2 degrees F. Each segment is treated as isothermal. The total pressure drop equals the sum of the pressure drops in the individual segments. The program also calculates the internal film heat transfer coefficient used in the temperature drop calculations. Correlations are included for the laminar and turbulent heat transfer coefficients for both Newtonian and generalized non-Newtonian liquids. The effects of friction heating and of pressure on viscosity are optional and can be switched on and off when editing the data. HOTOIL assumes that pressure affects only the viscosity-like parameter of non-Newtonian liquids.
F3.0
HOTOIL Limitations
As its name implies, HOTOIL is useful only for single-phase liquid hydrocarbon cases. Several internal fluid property correlations for liquid hydrocarbons would introduce errors for systems with non-hydrocarbon fluids. HOTOIL will accept up to 15 sections of line with different geometry, elevations, external heat transfer coefficients, etc. However, these sections must be connected in series. The program cannot handle piping systems with parallel lines or branching networks. HOTOIL is a steady-state program. It does not address fluid transients associated with opening and closing valves, or starting and stopping pumps.
March 1997
F-2
Chevron Corporation
Appendix F
F4.0
Accessing HOTOIL
HOTOIL is available three ways: 1. By accessing the Technical Standards Web site at http://techstds.rrc.chevron.com/tech_standards/ and downloading the program to your hard disk per the instructions given. By calling, e-mailing, or faxing Technical Standards customer service to request a floppy disk from which the program can be installed to your hard disk. E-MAIL: TECHSTDS (use global address list) CTN: 242-7241 Fax CTN: 242-2157 By purchasing the Chevron Process Simulation Suite, a CD-ROM available from the Process Simulation Team at CRTC.
2.
3.
F5.0
6.
F5.2
Chevron Corporation
F-3
March 1997
Appendix F
F6.0
Running HOTOIL
The following suggestions will help you make efficient use of HOTOIL: Modify sample data files rather than building input data from scratch. Use the F9 key to access help-screens that will provide explanations for the topics where the cursor is currently placed. Give your data files new names to preserve the sample data files. Questions should be directed to the HOTOIL support contact (see Section F8.0) HOTOIL is operated from a single screen. This screen displays all the data necessary to define a flowing pipeline. Large tables are handled by scrolling sections of the screen. One scrolling section displays three lines of a pipeline geometry data table. The other scrolling section displays three lines of fluid properties as a function of temperature. The user controls the operation of HOTOIL with function keys. These keys are defined at the bottom of the screen. The functions that are available at any moment are shown with a dark background. The function key definitions are shown below. F1 - LOAD a data file. You will be asked to point to the file name. F2 - SAVE the current data file. You will be asked to give the file name. F3 - EDIT or change the loaded data file. F4 - ZERO or erase the values in the loaded data file to start afresh. F5 - Run the program with the data shown. F6 - STOP running the program. (This will not exit from HOTOIL.) F7 - INSert a new line into Configuration or Viscosity table. F8 - DELete a line from Configuration or Viscosity table or DELete a file. F9 - HELP information on using HOTOIL. The help is context sensitive. Alt-F9 - HELP information on HOTOIL output. F10 - EXIT the program and return to DOS Alt-F10 - Changes the monitor display to color or monochrome as needed. Monochrome may give better results on some color monitors.
F7.0
March 1997
F-4
Chevron Corporation
Appendix F
Fig. F-1
Fig. F-2
! 12 5 585 6 43/ 4 743 !; ! 12 5 495 6 764 4 5074 !; ! 12 5 /48 6 9/0 4 54:6 !;
56/000
600000
, ! $( , . = 51 61 > 41 31? /1 71 :1
81
501@
F8.0
Program Support
If you have trouble running HOTOIL on your system or have questions about the program inputs or outputs, please contact Rob Hohmann, CRTC Richmond, CTN-242-2216. Comments and suggestions concerning the HOTOIL program are also welcomed.
Chevron Corporation
F-5
March 1997
Abstract
This section explains the use of the Chevron Research and Technology Company computer program HOTOL*. This program calculates pressure drop, outlet temperature, and heat loss vs. flow rate in a hot oil pipeline. Contents G1.0 G2.0 G2.1 G2.2 G3.0 G3.1 G3.2 G4.0 G5.0 G6.0 G7.0 G8.0 G8.1 G8.2 G8.3 G8.4 G8.5 G9.0 Introduction Method of Calculation Establishing Zones Calculating Pressure Drop Input Data Format Number of Sections in Series (Line 1) Section Variables (Line 1+n) Data Input and Program Execution Program Output Limitations Revision History Equations Used In HOTOL* Basic Relations Property Variations Flow Regime Boundaries and Friction Factors Derivation of Iteration Equations Corrections to Isothermal Pressure Drop References G-17 G-5 G-5 G-7 G-7 G-7 G-3 Page G-2 G-2
Chevron Corporation
G-1
January 1990
Appendix G
G1.0
Introduction
If fluid temperature is constant along a pipe, then pressure drop is a function of the liquid properties, flow rate, and line size (see Section 410). In a hot oil pipeline, however, the temperature of the liquid is higher than the ambient temperature, and is not constant, but decreases as the oil flows through the pipe. As the temperature of the oil changes, its liquid properties change, particularly the viscosity. The HOTOL* computer program evaluates these temperature changes and accounts for changes in the following liquid properties of the oil: specific gravity, density, bulk modulus, specific heat, thermal conductivity, and kinematic viscosity. HOTOL* does not consider friction heating, or pressure correction to viscosity, nor does it accommodate elevation profiles. To be analyzed by HOTOL* a pipeline must consist either of one section or of sections connected in series. Each section is described by its length, outside diameter, wall thickness, internal roughness, ambient temperature, external heat transfer coefficient, number of identical parallel lines, and whether or not the fluid is mixed at the section inlet (see Section G3.0).
G2.0
Method of Calculation
HOTOL* defines each relevant liquid property as a function of temperature. These correlations are shown in Section G8.2. This makes it possible to express the Reynolds number as a function of temperature and pipeline geometry for any given flow rate. The program can then define the following flow regimes in terms of the Reynolds numbers at the boundaries of the regimes: Laminar flow Transition flow Turbulent flow Hyperturbulent flow
Both the Reynolds number for the turbulent/hyperturbulent boundary and the friction factor in the hyperturbulent regime are functions of internal roughness and inside diameter. Flow regime boundaries and friction factors are defined in Section G8.3. Given the pipeline data, the system inlet temperature and flow rate, and using the flow regime boundaries defined in Section G8.3, HOTOL* begins the calculation. The following terms are used in this discussion: Section. A length of pipe over which all relevant geometrical and heat transfer variables are constant. Sections are described thoroughly by the input data Zone. A subdivision of a sectionmade by the programthat has a single flow regime
January 1990
G-2
Chevron Corporation
Appendix G
G2.1
Establishing Zones
First, HOTOL* computes the temperature at each flow regime boundary within each section. It then divides the pipeline into a number of zones. HOTOL* establishes the end of one zone and the beginning of the next where either of the following occurs: The temperature of the liquid decreases enough to cause a change in the flow regime The end of the current section is reached
Within each zone the pipeline section variables remain constant; the Reynolds number changes, but the flow regime remains the same. Once established, the pipeline zones are the basic subdivisions on which all subsequent calculations are performed. When a pipeline zone is terminated by a change in flow regime, an iterative heat balance (see Section G8.4) is used to generate the length of the zone. A similar iteration is used to calculate the zone outlet temperature when the end of the zone is also the end of the section. In both cases the outlet temperature for one zone is the inlet temperature for the next.
G2.2
G3.0
This is followed by one line for each section (n = total number of sections): Lines 1+n; (n=1 to 4). Number of lines in parallel (1, 1, 1, 1); flow is mixed at the inlet of this section? ('yes', 'no', 'no', 'no'); length in feet (4000, 18600, 18600, 4000); outside diameter in inches (21, 21, 21, 21); wall thickness in
Chevron Corporation
G-3
January 1990
Appendix G
Fig. G-1
inches (0.375, 0.375, 0.375, 0.375); roughness in inches (0.0018, 0.0018, 0.0018, 0.0018); external heat transfer coefficient in Btu/hr/sq. ft./F (0.25, 13.3, 0.432, 0.25); ambient temperature in degrees F (60, 48, 48, 60) This is followed by one line for the entire pipeline (n = total number of sections): Line 2+n; (n=4). Name of fluid (18-character maximum)('FUTURE A-960'), gravity (15.9), UOP characterization factor (11.7), number of temperature/viscosity data pairs(2)
This is followed by one line for each data pair (m = number of data pairs): Line 2+n+m; (n=4, m=1 to 2). Temperature (100, 210), kinematic viscosity (centistokes)(854, 35.5)
This is followed by one line for the entire pipeline: Line 3+n+m; (n=4, m=2). Temperature of fluid at inlet (100) Line 4+n+m. First flow rate (60F bbl/hr) (500), last flow rate (2000), increment (1500)
G3.1
G3.2
January 1990
G-4
Chevron Corporation
Appendix G
Flow is mixed at the inlet of this section (yes/no). Consider the flow to be mixed at: The entrance to the first section A division into parallel lines A junction of parallel lines An increase or decrease of inside diameter by a factor of 1.2 or more
Internal roughness. A value of 0.0018 inches (0.00015 feet) is common for steel pipe External Heat Transfer Coefficient. For buried and submarine lines, refer to Section 900.
G4.0
G5.0
Program Output
HOTOL* output (see Figure G-2) includes a summary of some of the input data. The results of the computer calculations are appropriately labeled. For each flow rate, the program prints out the following: Liquid temperature at the outlet of each section Pressure drop (psi) and heat loss rate (million Btu/Hr) for each individual section Cumulative pressure drop and heat loss rate for the pipeline, up to the outlet of the section in question
Chevron Corporation
G-5
January 1990
Appendix G
Fig. G-2
January 1990
G-6
Chevron Corporation
Appendix G
G6.0
Limitations
HOTOL* applies only to liquid phase hydrocarbon pipelines because of the fluid property vs. temperature correlations used HOTOL* uses internal heat transfer coefficient correlations that are valid for liquid cooling only. If the inlet temperature to the pipeline system is below the ambient temperature the program will not run. If the liquid should cool to below the ambient temperature at any point in the pipeline, the program will preemptively assume that the liquid temperature and ambient temperature are equal, and continue the calculations on an isothermal basis for as long as necessary. The pressure drops calculated on this basis will be too low. In this case, a warning to this effect will appear in the output (see sample output) HOTOL* is capable of handling pipelines consisting of fifty series sections or less in its standard form. If the number input is greater than fifty, the program will print an error message. HOTOL* can be extended to handle any number of series sections. If greater capacity is needed contact CRTCs computer systems analyst HOTOL* cannot handle branch lines or loop lines other than multiple identical parallel lines The HOTOL* analysis is valid only for Newtonian fluids. It is not applicable if the fluid in question is cooled below its pour point temperature
G7.0
Revision History
HOTOL* was originally prepared by J. C. Sebastian (May 1968). The program logic was revised and expanded by W. A. Ebert (September 1976). G. H. Horne subsequently modified the program for multiple sections in series (November, 1977). The documentation (Design Practice A-871-3) was rewritten by G. W. Williams (October 1978) and revised by G. H. Horne (February 1979). Design Practice A-871-3 was republished as Appendix G of the Fluid Flow Manual in October 1988.
G8.0
(Eq. G-1)
Chevron Corporation
G-7
January 1990
Appendix G
where: L = length of zone, ft f = friction factor = liquid density, lbm/ft3 Di = inside diameter, ft V = average liquid, velocity ft/sec P = pressure drop, lbf/ft2 g = conversion factor, ft lbm/sec2 lbf
where: Q = rate of heat loss, Btu/sec w = mass flow rate, lbm/sec Cp = specific heat at average section bulk temperature, Btu/lbm/F T = change in bulk temperature from inlet to outlet, F
G8.2
Property Variations
The variations of fluid properties with temperature are calculated as follows:
(Eq. G-3)
Density
Per Section G9.0, reference 1. = D1 + D2(T - 60) + D3(T - 60)2
(Eq. G-4)
January 1990
G-8
Chevron Corporation
Appendix G
D2 = 62.4 [(8 x 10-4)SG60 - (1.05 x 10-3)] D3 = 62.4 [(6 x 10-7)SG60 - (4.9 x 10-7)]
Bulk Modulus
(Eq. G-5)
Specific Heat
Per Section G9.0, reference 1. Cp = Co (C1 + C2T)
(Eq. G-6)
where: Cp = specific heat, Btu/lbm/F T = temperature, F Co = 0.055 Kw + 0.35 C1 = 0.681 - 0.308 SG60 C2 = 0.000815 - 0.000306 SG60 Kw = Watson (or UOP) characterization factor=(Tb)1/3/SG60 Tb = normal (average) boiling point, R SG60 = specific gravity at 60
Thermal Conductivity
k = K1 + K 2 T
(Eq. G-7)
where: k = thermal conductivity, Btu/Hr F ft T = temperature, F K1 = 0.08934 - (6.7407 x 10-4) (API) + (5.77 x 10-6) (API)2 K2 = -2.976 x 10-5 + (1.104 x 10-6)(API) - (2.201 x 10-8) (API)2
Chevron Corporation
G-9
January 1990
Appendix G
Kinematic Viscosity
Input pairs of kinematic viscosity and temperature, (, T), are fitted to a curve of the form: ln[1 + ln(2)] = A - B ln(T + 460)
(Eq. G-8)
(Eq. G-9)
(Eq. G-10)
G8.3
January 1990
G-10
Chevron Corporation
Appendix G
(Eq. G-11)
(Eq. G-12)
G8.4
(Eq. G-13)
integrating:
(Eq. G-14)
Chevron Corporation
G-11
January 1990
Appendix G
Fig. G-3
= = = = = = = = = = = = = = = = = = = = = = = = =
length of zone distance from inlet of zone outside diameter of pipe inside diameter of pipe internal heat transfer coefficient external heat transfer coefficient overall heat transfer coefficient mass flow rate bulk temperature at position X bulk temperature at inlet of zone bulk temperature at outlet of zone bulk temperature averaged over the zone length wall temperature averaged over the zone length ambient temperature DiL = zone surface area distance from last point of mixing Reynolds number based on tb Prandtl number based on tb Grashof number with fluid properties evaluated at tb, t = (tb-tw) Prandtl number based on tw Grashof number with fluid properties evaluated at tw, t = (tb-tw) dynamic viscosity based on tb dynamic viscosity based on tw fluid velocity Nusselt number
January 1990
G-12
Chevron Corporation
Appendix G
(Eq. G-15)
(Eq. G-16)
(Eq. G-17)
where:
(Eq. G-18)
(Eq. G-19)
Calculation of tb
(Eq. G-20):
Chevron Corporation
G-13
January 1990
Appendix G
(Eq. G-21)
(Eq. G-22)
(Eq. G-23)
where: For a laminar regime with Reb 2000, the Nusselt number is the greatest of either:
(Eq. G-24)
January 1990
G-14
Chevron Corporation
Appendix G
For the transition zone 2000 < Reb < 104: Nu = * Nu2000 + (1 - *) Nu10,000
(Eq. G-26)
where: * = 1.25 - Reb/8000 Nu2000 = Nu @ Reb of 2000 per Eq. G-24 or 25 Nu10,000 = Nu @ Reb of 10,000 per Eq. G-27 For turbulent and hyperturbulent regimes with Reb 104:
(Eq. G-27)
For any zone, an iterative manipulation of Equations G-16 through G-23 can be used in either of two ways: If tbo (the bulk temperature at the outlet) is given, then L (the length of the zone) can be found If L is given, then tbo can be found
G8.5
where correction factors , , and Z are defined as follows: Radial Property Variation () (Correction for Mw/Mb not = 1)
Chevron Corporation
G-15
January 1990
Appendix G
(Eq. G-29)
Natural Convection ()
(Eq. G-30)
(Eq. G-31)
but
(Eq. G-32)
January 1990
G-16
Chevron Corporation
Appendix G
so
(Eq. G-33)
(Eq. G-34)
so
(Eq. G-36)
or
(Eq. G-37)
G9.0
References
1. Heat Transfer Research Inc. Design Manual.
Chevron Corporation
G-17
January 1990
Chevron Corporation
H-1
January 1990
Appendix H
January 1990
H-2
Chevron Corporation
Appendix H
Chevron Corporation
H-3
January 1990
Appendix H
January 1990
H-4
Chevron Corporation
Appendix H
Chevron Corporation
H-5
January 1990
Appendix H
January 1990
H-6
Chevron Corporation
Appendix H
Chevron Corporation
H-7
January 1990
Appendix H
January 1990
H-8
Chevron Corporation
Appendix H
Chevron Corporation
H-9
January 1990
Appendix H
January 1990
H-10
Chevron Corporation
Appendix H
Chevron Corporation
H-11
January 1990
Appendix H
January 1990
H-12
Chevron Corporation
Appendix H
Chevron Corporation
H-13
January 1990
Appendix H
January 1990
H-14
Chevron Corporation
Appendix H
Chevron Corporation
H-15
January 1990
Appendix H
January 1990
H-16
Chevron Corporation
Appendix H
Chevron Corporation
H-17
January 1990
Appendix H
January 1990
H-18
Chevron Corporation
Appendix H
Chevron Corporation
H-19
January 1990
Appendix H
January 1990
H-20
Chevron Corporation
Appendix H
Chevron Corporation
H-21
January 1990
Appendix H
January 1990
H-22
Chevron Corporation
Appendix H
Chevron Corporation
H-23
January 1990
Appendix H
January 1990
H-24
Chevron Corporation
Appendix H
Chevron Corporation
H-25
January 1990
Appendix H
January 1990
H-26
Chevron Corporation
Appendix H
Chevron Corporation
H-27
January 1990
Appendix H
January 1990
H-28
Chevron Corporation
Appendix H
Chevron Corporation
H-29
January 1990
Appendix H
January 1990
H-30
Chevron Corporation
Appendix H
Chevron Corporation
H-31
January 1990
Appendix H
January 1990
H-32
Chevron Corporation
Appendix H
Chevron Corporation
H-33
January 1990
Appendix H
January 1990
H-34
Chevron Corporation
Appendix H
Chevron Corporation
H-35
January 1990
Appendix H
January 1990
H-36
Chevron Corporation
Appendix H
Chevron Corporation
H-37
January 1990
Appendix H
January 1990
H-38
Chevron Corporation
Appendix H
Chevron Corporation
H-39
January 1990
Appendix H
January 1990
H-40
Chevron Corporation
Appendix H
Chevron Corporation
H-41
January 1990
Appendix H
January 1990
H-42
Chevron Corporation
Appendix H
Chevron Corporation
H-43
January 1990
Appendix H
January 1990
H-44
Chevron Corporation
Appendix H
Chevron Corporation
H-45
January 1990
Appendix H
January 1990
H-46
Chevron Corporation
Appendix H
Chevron Corporation
H-47
January 1990
Appendix H
January 1990
H-48
Chevron Corporation
Appendix H
Chevron Corporation
H-49
January 1990
Appendix H
January 1990
H-50
Chevron Corporation
Appendix H
Chevron Corporation
H-51
January 1990
Appendix H
January 1990
H-52
Chevron Corporation
Appendix H
Chevron Corporation
H-53
January 1990
Appendix H
January 1990
H-54
Chevron Corporation
Appendix H
Chevron Corporation
H-55
January 1990
Appendix I.
PCSURGE Program
Abstract
This document describes the computer program PCSURGE and gives guidance on its use. PCSURGE is a menu driven input and output system built around the original SURGE program, and is adapted for IBM PC compatible machines with a math co-processor. SURGE was previously available only on Chevrons VM mainframe computer. This guide to PCSURGE is not intended to be a complete description of the SURGE program. The SURGE program is described in detail in Appendix H of this manual. Users of both SURGE and PCSURGE should refer to Appendix H for an explanation of the programs capabilities, inputs, and outputs. Contents I1.0 I2.0 I3.0 I3.1 I3.2 I4.0 I5.0 I6.0 I6.1 I6.2 I7.0 General Description Accessing PCSURGE Installing PCSURGE on Your Computer Installing PCSURGE from the Technical Standards Website Running PCSURGE from a Floppy Disk Running PCSURGE Note on Out of Environment Space Error Description of Menu Options Main Menu Edit Data Menu Program Support I-7 I-3 I-3 I-4 Page I-2 I-2 I-2
Chevron Corporation
I-1
March 1997
Appendix I
I1.0
General Description
PCSURGE was created to allow quick, interactive creation and manipulation of SURGE input data files. It allows a user to load an input data file, edit the data, save the new data file, run the SURGE calculations, and view the output on the screen or send the output to a file, all from within the menu-driven environment. The menu system was designed to be self-explanatory. Most users who are familiar with the SURGE program should be able to use PCSURGE with minimal reference to the documentation for SURGE in Appendix H or the description of PCSURGEs features presented below. For all but the simplest systems, it is important that the user start with a good sketch of the system to be modeled, with all branches and nodes numbered. (Branches are pipe lengths of constant diameter. Nodes are pumps, tanks, junctions between branches, etc.) The sketch should include geometry, elevations, initial flow rates, pressures, etc. For a complete list of the program inputs, see the Description of Menu Options, Section I6.0, below.
I2.0
Accessing PCSURGE
PCSURGE is available three ways: 1. By accessing the Technical Standards Web site at http://techstds.rrc.chevron.com/tech_standards/ and downloading the program to your hard disk per the instructions given. By calling, e-mailing, or faxing Technical Standards customer service to request a disk from which the program can be installed to your hard disk. E-MAIL: TECHSTDS (use global address list) CTN: 242-7241 Fax CTN: 242-2157 By purchasing the Chevron Process Simulation Suite, a CD-ROM available from the Process Simulation Team at CRTC.
2.
3.
I3.0
March 1997
I-2
Chevron Corporation
Appendix I
5.
After you have saved pcsrge.exe to your hard disk, double-click the filename from Windows Explorer (Windows 95) or File Manager (Windows 3.x) to decompress the program files. After the program files have been decompressed, double-click on the file called pcsurge.exe to run the program.
6.
I3.2
I4.0
Running PCSURGE
To run PCSURGE it is necessary to have the files named PCSURGE.EXE and SURGE51.EXE in the same sub-directory on your hard disk or on a non-writeprotected floppy disk. Data files for use by PCSURGE must be in the same subdirectory. PCSURGE can be run without data files, since data can be entered directly into the program. However, it is often easier to load a data file, such as SAMPLE.DAT provided with PCSURGE, and then edit that data. To start PCSURGE type PCSURGE at the DOS prompt. Use the up and down arrow keys to highlight the desired operations from the Main Menu. Press ENTER to choose the highlighted option.
I5.0
Chevron Corporation
I-3
March 1997
Appendix I
I6.0
The following is a discussion of each of these options: Run Surge with Sample Data. When this option is chosen, the sample data file, SAMPLE.DAT, will be loaded and run. The user will see the SURGE output displayed on the screen while the program runs. The run will take about 30 seconds on a PC with a 386 central processor. The run can be terminated at any time with CTRL-BREAK (hold down the CONTROL key while pressing BREAK). Show Data Files. Choosing this option will cause all the files in the directory with the extension .DAT to be displayed on the screen. No checking is done to be sure they are actually SURGE data files. Retrieve Data File. This option allows the user to enter the name of a data file which will be retrieved and placed in memory for editing or running. If a data file has already been retrieved, then the name of that data file is displayed. If the user types a file name that doesnt exist in the directory, a message is printed on the screen and the user is instructed to try again. The user can also return to the Show Data Files option to review the data files that can be retrieved. A message on the Retrieve Data File screen tells the user that the extension .DAT does not need to be typed. It will be added automatically. If any other extension is typed, it will be ignored. Only files with the extension .DAT can be retrieved in PCSURGE. Edit Data. The Edit Data option of the Main Menu brings up an Edit Menu which gives the user five options for editing sections of the data file. The details of these options are discussed in Section I6.2. The sixth option on the Edit Menu is Return to Main Menu. Save Data. This option allows the user to type a file name. A message on the Save Data File screen tells the user that the extension .DAT does not need to be typed. It will be added automatically. If any other extension is typed, it will be ignored. Run Surge Program. Choosing this option will display the following choices for running the calculations with the data currently in memory: Run SURGE. Print output on SCREEN
March 1997
I-4
Chevron Corporation
Appendix I
Run SURGE. Print output to FILE: sample.OUT Save Current Data Retrieve Data File Return to Main Menu
The first two of these options show that the output from the program can either be routed to the screen or to a file. The user will probably find it useful to initially route the output to the screen while the data file is being developed and modified. Once the output and results are in the most useful format, the output can be saved in a disk file (with extension .OUT) for printing, or for manipulation in a word processor before printing. If there is not enough data in memory to run the calculations when the run option is chosen, a message to that effect will be displayed. The user can then continue editing the data in memory or retrieve another data file. Exit to DOS. The last option on the Main Menu is Exit to DOS. PCSURGE keeps track of the sequence in which data is edited and saved. If the data in memory has not been saved since it was edited, or has not been given a file name, the user will be given another opportunity to save that data before the final exit to DOS.
I6.2
Title. This option gives the user the opportunity to enter a 50 character title which will be printed on the program output for documentation purposes. The user can, however, ignore this option since no title is required. Fluid Properties. Values for the following three liquid properties must be entered to allow proper execution of the SURGE calculations. Bulk Modulus, psi (See Section 800 of this manual for typical values.) Density, lbm/ft3 Vapor pressure, psia
Boundary Conditions. The Boundary Condition Menu allows the user to add or delete boundary conditions. Boundary conditions are added to the data file one by one from the following list which is displayed on the screen. System Inlet: Reservoir or Pump Initially Open Simple Closing Outlet Valve Two Branch Junction, Optional Relief Valve Three Branch Diverging Junction
Chevron Corporation
I-5
March 1997
Appendix I
Three Branch Converging Junction Dead End or Constant Head Outlet Head Discontinuity: Pump or Head Loss Surge Chamber Initially Open Outlet Valve, Stepped Close Initially Open Interior Valve Pump w/ Pres. Controlled Relief to Suction Initially Closed Outlet Valve/Stepped Close Initially Closed Simple Close Interior Valve Initially Shutdown Interior Pump
After a specific boundary condition has been chosen from the list above, the user must enter the specific parameters needed to model that boundary condition. PCSURGE presents the user with an abbreviated explanation of these parameters. If there is any question about a boundary condition parameter the user can refer to the full explanation of each value in Appendix H of this manual. Each boundary condition and branch is numbered. The user specifies the incoming and outgoing branch numbers associated with each boundary condition. Branches. Like boundary conditions, branches are added and deleted one at a time. After a branch is added, the user must define the following parameters for each branch. A more thorough explanation of each of these values is given in Appendix H of this manual. Outside Diameter of Pipe, in. Thickness of Pipe Wall, in. Length of Pipe Section, ft. Flow Rate in Pipe Section, lbm/hr. Elevation at Upstream end of Branch, ft. Elevation at Downstream end of Branch, ft. Psig at Upstream end of branch, not including elevation head Psig at Downstream end of branch, not including elevation head Thickness of Concrete Lining Outside of Pipe, in. Pipe Elastic Modulus, psi. (Enter 0 for default for steel.)
Program Control. This option of the Edit Menu allows the user to control how the SURGE calculations proceed and how the output is formatted. A new user usually needs to experiment with the options for formatting printouts for various branches in the system being modeled. Fortunately, PCSURGEs menu interface allows a user to quickly go back and forth between making changes in the data file and seeing the results of those changes displayed on the screen. Total Run Time, sec. Number of nodes in shortest branch(es) (2 to 100) No. of the Branch for which summary results will be printed End of Branch for summary results printout, 1=inlet, 0=outlet
March 1997
I-6
Chevron Corporation
Appendix I
Plot P vs. T (Press vs. Time) for above point? 1=yes, 0=no Time between printouts, sec. Plot P vs. T + Print P, Head, Vel for points below? 1=yes, 0=no Branch No. of point 1 Position of point 1, 1=inlet, 0=outlet Branch No. of point 2 Position of point 2, 1=inlet, 0=outlet Branch No. of point 3 Position of point 3, 1=inlet, 0=outlet Branch No. of point 4 Position of point 4, 1=inlet, 0=outlet Branch No. of point 5 Position of point 5, 1=inlet, 0=outlet Branch No. of point 6 Position of point 6, 1=inlet, 0=outlet
It can also be useful for a new user to examine the program control parameters in the sample data file, SAMPLE.DAT, and observe the resulting output format by running the sample data. A sample data file and corresponding output file are shown in Appendix H, starting on page H-32, of this manual.
I7.0
Program Support
If you have trouble running PCSURGE on your system or have questions about the program inputs or outputs, please contact Rob Hohmann, CRTC Richmond, CTN242-2216. Comments and suggestions concerning the PCSURGE program are also welcomed.
Chevron Corporation
I-7
March 1997