Category: Refining

  • Correlations for Conversion between True and Reid Vapor Pressures (TVP and RVP)

    Accurate measurement and prediction of crude oil and natural gas liquid (NGL) products vapor pressure are important for safe storage and transportation, custody transfer, minimizing vaporization losses and environmental protection. Vapor pressure specifications are typically stated in Reid Vapor Pressure (RVP) or/and True Vapor Pressure (TVP). In addition to the standard procedures for their measurements, there are rigorous and shortcut methods for their estimation and conversion.

    There are figures and monographs for conversion of RVP to TVP for NGLs (Natural Gas Liquids) and crude oil at a specified temperature. This tip will present simple correlations for conversion from RVP to TVP and vice versa at a specified temperature. The correlations are easy to use for hand or spreadsheet calculations. Figures generated using these correlations will be presented, too.

    Chapter 5 of reference [1] presents an excellent overview of TVP and RVP including their definitions, standard procedures for their measurements and diagrams for their conversion. The proceeding two paragraphs are extracted with minor revisions from reference [1].

    TVP is the actual vapor pressure of a liquid product at a specified temperature and is measured with a sample cylinder. TVP specifications must always be referenced to a temperature, which frequently falls between 30-50°C (86-122°F). TVP is difficult to measure and depends on the ratio of the vapor to liquid, V/L, in the measurement device. If V/L = 0, the vapor pressure is essentially equivalent to the bubble point of the mixture which is the highest vapor pressure value for the liquid. As V/L increases, i.e., a small amount of vapor exists at the point of measurement, the measured vapor pressure will decrease. The relationship between the measured TVP and V/L depends on the composition of the mixture. For “near pure” component mixtures, V/L has little effect on the measured vapor pressure. For mixtures with a large composition range, such as crude oil or condensate, the effect of V/L on the measured vapor pressure can be significant (See ASTM D 6377 – 10 for detail). Reference [1] lists the Standards used for TVP measurements.

    Because of the difficulty in accurately measuring TVP, an alternative method of measuring vapor pressure is frequently used. This is the RVP.  The RVP is a standard test set out in ISO 3007:1999, Petroleum products and crude petroleum – Determination of vapour pressure – Reid method. Another standard that applies to RVP is:  ASTM D323 – 08, Standard Test Method for Vapor Pressure of Petroleum Products (Reid Method). The RVP test is applicable for crude oils, condensates, and petroleum products such as gasoline (petrol) mixtures. A liquid sample is collected in the lower 20% chamber (see Fig 5.12 of reference [1]). The 80% chamber, which is filled with air (may also contain a small amount of water vapor) at atmospheric pressure.  Both chambers are cooled to 0°C (32°F). The 80% air chamber (at atmospheric pressure) is then connected to the 20% liquid chamber. The connecting valve is opened and the cylinder is heated in a water bath to 37.8°C (100°F). The pressure indicated on the gauge is the RVP.

    Reference [1] provides Figure 5.14 for conversion RVP to TVP for motor gasoline (petrol) and natural gasoline (C5+ NGLs) at various temperatures. Figure 5.15 of reference [1] is a nomograph that shows the approximate relationship between RVP and TVP for crude oil. It is commonly used for converting RVP to TVP at custody transfer points where the vapor pressure specification for the oil is a TVP, but the actual vapor pressure measurement is an RVP.

    Vazquez-Esparragoza et al. [2] present an algorithm to calculate RVP without performing the actual test. The algorithm, based on an air-and-water free model, uses the Gas Processors Association Soave-Redlich-Kwong [3] equation of state and assumes liquid and gas volumes are additive. This algorithm can be used to predict RVP of any hydrocarbon mixture of known composition. Since the calculations are iterative, it should be incorporated into a general purpose process simulator. Vazquez-Esparragoza et al. [2] reported good agreement between predicted and experimental values.

    Riazi et al. [4] presented a new correlation for predicting the RVP of gasoline and naphtha based on a TVP correlation. The input parameters for this correlation are the mid-boiling point, specific gravity, critical temperature, and critical pressure, where the critical properties may be estimated from the boiling point and specific gravity using available methods. They evaluated their proposed correlation with data collected on 50 gasoline samples from crude oils from around the world with API gravity ranges from 41 to 87, average boiling point ranges from 43 to 221 °C (110 to 430 °F) and RVP of 0.7 to 115 kPa (0.1–17 psi). The average error from their proposed correlation is about 6 kPa (0.88 psi) [4].

    Development of Model for Motor Gasoline and Natural Gasoline

    To develop the desired correlations for conversion of motor gasoline and natural gasoline RVP to TVP and vice versa, this tip generated 127 data points from Figure 5.14 of reference [1]. These data covered the full ranges of temperature, RVP and TVP of Figure 5.14.

    RVP to TVP: This tip used these data to determine the parameters to the Equations 1 through 3. The API 2517 originally reported the same forms of equations for crude oil.

    eq1a-3a

    TVP to RVP: Similarly this tip proposes the following equations for conversion from TVP to RVP.

    eq1b-3b

    where:

    T          = Temperature, °C (°F)

    RVP    = Reid Vapor Pressure, kPa (psi)

    TVP    = True Vapor Pressure, kPaa (psia)

    Note that the values of A1, A2, B1, and B2 are different in the above two sets of equations. The value of “C” is a function of the chosen units (SI versus FPS) and is consistent.

    Table 1 presents the optimized values of A1, A2, B1, B2, and C for three sets of data in FPS (Foot-Pound-Second) and SI (System International). The data set labeled “All” included 127 data points covering all of the data of combined motor gasoline and natural gasoline. The data set “Motor Gasoline’ and “Natural Gasoline’ cover 76 and 51 data points for motor gasoline and natural gasoline, respectively. Table 1 also presents the Average Absolute Percent Deviation (AAPD), the Maximum Absolute Percent Deviation (MAPD), Average Absolute Deviation (AAD), and the number of data points (NP) for each data set. The error analysis of Table 1 indicates that the accuracy of the proposed correlations is good. Their accuracy is as good as the quality of tabular data generated from Figure 5.14.

    Table 1. The optimized parameters for motor gasoline and natural gasoline

    tab1

    1AAPD           = Average Absolute Percent Deviation
    2MAPD           = Maximum Absolute Percent Deviation
    3AAD              = Average Absolute Deviation
    4NP                 = Number of data Points (NP)

    Figures 1a (SI) and 1b (FPS) present the tabular data (generated from Figure 5.14 of reference [1]) in the form of “legends”. This figure also presents the predicted TVP by the proposed correlations (Equations 1 through 3 and the corresponding parameters listed in Table 1) as a function of temperature and RVP in the form of “continuous lines” and “broken lines” for motor gasoline and natural gasoline, respectively.

    Figure 1a. TVP as a function of RVP and temperature for motor gasoline and natural gasoline (typical C5+ NGLs)
    Figure 1b. TVP as a function of RVP and temperature for motor gasoline and natural gasoline (typical C5+ NGLs)
    Figure 1b. TVP as a function of RVP and temperature for motor gasoline and natural gasoline (typical C5+ NGLs)

    Development of Model for Crude Oil

    To develop the desired correlations for conversion of crude oil RVP to TVP and vice versa, this tip generated 196 data points from Figure 5.15 of reference [1] using the correlations reported in API 2517. The API 2517 equations are the same as shown in Equations 1 through 3. These data covered the full ranges of temperature, RVP and TVP of Figure 5.15.

    RVP to TVP: Table 2 presents the API 2517 FPS parameters of A1, A2, B1, B2, and C for Equations 1a, 2a, and 3a (T, °F, RVP, psi and TVP, psia). This tip determined the SI (T, °C, RVP, kPa and TVP, kPaa) corresponding parameters. Table 2 also presents the SI parameters.

    Table 2. API 2517 parameters for crude oil TVP calculations
    (For use with Equations 1a, 2a and 3a)

    tab2

    Figures 2a (SI) and 2b (FPS) present the predicted TVP by Equations 1a, 2a and 3a and the corresponding parameters listed in Table 2 as a function of temperature and RVP for crude oil.

    TVP to RVP: To cover the full ranges of Figure 5.15 of reference [1], this tip proposes the following equations for conversion of TVP to RVP.

    fig1c-3c

    where:

    T          = Temperature, °C (°F)

    RVP    = Reid Vapor Pressure, kPa (psi)

    TVP    = True Vapor Pressure, kPaa (psia)

    Table 3 presents the optimized values of A1, A2, A3, B1, B2, B3, and C for RVP calculation in FPS and SI for the proposed Equations 1c through 3c. Table 3 also presents the Average Absolute Percent Deviation (AAPD), the Maximum Absolute Percent Deviation (MAPD), Average Absolute Deviation (AAD), and the Number of data Points (NP). The error analysis of Table 3 indicates that the accuracy of the proposed correlations is good and can be used for the estimation purposes.

    Table 3. The optimized parameters for crude oil TVP conversion to RVP

    (For use with Equations 1c, 2c & 3c)

    tab3

    1AAPD           = Average Absolute Percent Deviation

    2MAPD           = Maximum Absolute Percent Deviation

    3AAD              = Average Absolute Deviation

    4NP                 = Number of data Points (NP)

     

    Conclusions:

    For converting RVP to TVP of motor gasoline and natural gasoline, this tip presented simple correlations similar to the ones reported in API 2517 for crude oil RVP to TVP. This tip determined the correlation parameters by regressing the data generated from available diagrams. These correlations were also extended for easy conversion from TVP to RVP. Tables 1, 2, and 3 present the correlation parameters in SI and FPS system of units. Tables 1 and 3 presents the accuracy of the proposed correlations against the data generated from diagrams. Tables 1 and 3 indicate that the accuracy of these correlations is as good as the quality of the data in the original diagram and can be used for easy conversion of RVP to TVP or vice versa. These correlations are easy to use for hand or spreadsheet calculations and should be used for estimation purposes. For accurate measurements, correction procedures outlined in ASTM D 6377–10 and other guidelines should be consulted.  Several organizations are currently working to improve the accuracy of TVP estimation from RVP and/or VPCR(x) (ASTM D6377) measurement techniques. In all cases, Federal and State Laws and Regulations should be followed for safety and environmental protection.

    To learn more about similar cases and how to minimize operational problems, we suggest attending our G4 (Gas Conditioning and Processing), G5 (Advanced Applications in Gas Processing), and PF4 (Oil Production and Processing Facilities), courses.

    PetroSkills offers consulting expertise on this subject and many others. For more information about these services, visit our website at http://petroskills.com/consulting, or email us at consulting@PetroSkills.com.

    By: Dr. Mahmood Moshfeghian

    Figure 2a. TVP as a function of RVP and temperature for crude oil
    Figure 2a. TVP as a function of RVP and temperature for crude oil
    Figure 2b. TVP as a function of RVP and temperature for crude oil
    Figure 2b. TVP as a function of RVP and temperature for crude oil

    References:

    1. Campbell, J.M., Gas Conditioning and Processing, Volume 1: The Basic Principles, 9th Edition, 2nd Printing, Editors Hubbard, R. and Snow–McGregor, K., Campbell Petroleum Series, Norman, Oklahoma, 2014.
    2. Vazquez-Esparragoza, J.J., et al, “How to Estimate Reid Vapor Pressure (RVP) of Blends”, Bryan Research and Engineering, Inc, Bryan, Texas, 2015.
    1. Soave, Chem. Eng. Sci. 27, 1197-1203, 1972.
    2. Riazi, M.R., Albahri, T.A. and Alqattan, A.H., “Prediction of Reid Vapor Pressure of Petroleum Fuels”, Petroleum Science and Technology, 23: 75–86, 2005.
  • Adsorption Dehydration: Two-Tower vs Three-Tower System

    There are different process configurations for adsorption dehydration systems. The most common arrangements are two-tower and three-tower configurations. One can find the details of the adsorption dehydration process and the descriptions of equipment in Chapter 18 of John M. Campbell textbook [1]. Figures 1 and 2 present a simplified process flow diagram for two-tower and a three-tower configurations, respectively. These units can reduce the water content of a gas stream to less than 0.1 ppmv.  The gas industry normally uses adsorption dehydration units upstream of a liquefied natural gas (LNG) plant or a deep natural gas liquid (NGL) extraction plant where the gas temperature reduces to less than -160 °C (-256 °F) and -100 °C (-148 °F), respectively. Removal of water content to this very low level is essential to prevent freezing.

    Figure 1. A simplified process flow diagram for a two-tower adsorption dehydration system [1].
    Figure 1. A simplified process flow diagram for a two-tower adsorption dehydration system [1].
    In the two-tower system, while tower A is in the adsorption mode, tower B is regenerating. After tower A completes its adsorption cycle, it will switch to the regeneration mode and tower B starts its adsorption cycle. At any time one of the towers is adsorbing water while the other one is regenerating.

    In the three-tower configuration, at any time two towers (e.g. A and B) in staggered parallel adsorption while the third tower (e.g. C) is regenerating. In this configuration, half of the feed gas flow rate is going through tower A and the other half passes through tower B as shown in the cycle chart embedded in Figure 2.

    Figure 2. A simplified process flow diagram for a three-tower adsorption dehydration system [1].
    Figure 2. A simplified process flow diagram for a three-tower adsorption dehydration system [1].
    The PetroSkills’ May 2015 [2] and October 2015 tip of the month (TOTM) [3] discussed the efficient operation of molecular sieve dehydration units. Specifically, they discussed the benefits of standby time in the adsorption dehydration processes and impact of feed gas conditions.

    This month’s TOTM compares the required size of major equipment for the two-tower system with the three-tower system. The comparison considers the following parameters.

    1. Mass of Desiccant
    2. Bed Diameter
    3. Bed Height
    4. Regeneration Gas Rate
    5. Regeneration Heating Load
    6. Regeneration Cooling Load
    7. Regeneration Gas Heater Load
    8. Tower wall thickness and mass

    Table 1 presents descriptions of the two configurations.

    Table 1. Tower configurations

    tab1

    This discussion assumes a gas volume rate of 2.83×106 Sm3/d (100 MMSCFD) with externally insulated towers. It examines two feed conditions of:

    a. 30°C (86°F) and 6.207 MPaa (900 psia)

    • Estimated water content = 699 kg/106 Sm3 (43.7 lbm/MMSCF)
      • 2-Tower water load per tower per hour = 4 kg/h (182 lbm/hr)
        • Water load per tower per cycle = 989 kg (2185 lbm)
      • 3-Tower water load per tower per hour = 2 kg/h (91 lbm/hr)
        • Water load per tower per cycle = 662 kg (1456 lbm)

    b. 40°C (104°F) and 8.0 MPaa (1160 psia)

    • Estimated water content = 974 kg/106 Sm3 (61 lbm/MMSCF)
      • 2-Tower water load per tower per hour = 115 kg/h (254 lbm/hr)
        • Water load per tower per cycle = 1378 kg (3050 lbm)
      • 3-Tower water load per tower per hour = 5 kg/h (127 lbm/hr)
        • Water load per tower per cycle = 920 kg (2033 lbm)

    The other specified parameters are:

    1. Desiccant = Molecular Sieve Type 4A, 3.2 mm (1/8 in) diameter
    2. Desiccant loading capacity, ΔXnew = 19 wt% (mass of water/100 mass of desiccant)
    3. Desiccant life factor FL = 0.6 (Based on 3 years and average performance)
    4. Desiccant density = 705 kg/m3 (44 lbm/ft3)
    5. Desiccant heat capacity = 1.0 kJ/kg-°C (0.24 Btu/lbm-°F)
    6. Steel heat capacity = 0.5 kJ/kg-°C (0.12 Btu/lbm-°F)
    7. Feed gas relative density = 0.7
    8. Regeneration dry lean gas relative density = 0.59
    9. Regeneration gas pressure = 2.07 MPa (300 psia)
    10. Regeneration gas temperature to the heater = Feed gas temperature
    11. Regeneration gas temperature from heater = 288 °C (550 °F)
    12. Final bed regeneration temperature prior to cooling = 260 °C (500 °F).

    Calculation Results

    Based on the procedure and steps of Chapter 18 [1], this TOTM utilized a revised version of PetroSkills/Campbell GCAP software [4] to perform all of the calculations.

    Figure 3A shows variation of the required mass of desiccant per tower with the feed gas water content for two and three-tower systems.  As the feed water load increases, the required mass of desiccant increases for the specified and constant adsorption time. Note that the feed water load increases with increase in temperature and decrease in pressure. The water load per tower is a function of the feed water content, adsorption time and gas flow rate through the tower.

    Figure 3A. Mass of desiccant per tower vs the feed gas water content and number of towers

    The adsorption cycle times chosen for the two-tower and three-tower systems (12 hours and 16 hours, respectively) result in essentially the same total mass of water being adsorbed by both configurations. The total required mass of desiccant for the two-tower and three-tower systems are different because the length of mass transfer zone for the 3-tower system is lower than for the 2-tower system. Since desiccants are sold in 300 lbm (136.1 kg) increment, the mass of desiccant in each tower was rounded up to the next 300 lbm (136.1 kg).

    The calculated total mass of desiccants for the two-tower system are 19323 and 25855 kg (42600 and 57000 lbm) for the feed with the lower and higher water content, respectively. The corresponding desiccant masses for the three-tower system are 18779 and 25719 kg (41400 and 56700 lbm).  Similarly, Figure 3B presents the required mass of steel per tower as a function of feed gas water content and number of towers.

    Figure 3B. Mass of steel per tower vs the feed gas water content and number of towers
    Figure 3B. Mass of steel per tower vs the feed gas water content and number of towers

    The mass of the desiccant in the tower and pressure drop criteria establish the tower diameter and height of desiccant. This is a trial and error calculation. Figure 4 shows variation of the minimum required bed diameter with the feed gas water content and number of towers. Since the gas flow rate per tower of the two-tower system is two times higher than the gas flow rate per tower in the three-tower system, the diameter of each tower in the two-tower system must be larger to meet pressure drop criteria of less than 41 kPa (6 psi). Figure 4 indicates that the limiting factor is the pressure drop and the feed gas water content has small effect on the bed diameter. The calculated superficial gas velocity for all cases ranged from 0.10 to 0.15 m/s (20 to 30 ft/min).

    After determining the bed diameter, one can calculate the desiccant height from the mass of desiccant, bed diameter, and the desiccant density.  Figure 5 shows the variation of the minimum desiccant height with the feed gas water content and number of towers.  This figure indicates that feed with higher water content requires a taller bed and feed with lower water load requires shorter bed height. Since the diameter in two tower system was larger its height is shorter.

    Figure 4. Bed diameter vs the feed gas water content and number of towers
    Figure 4. Bed diameter vs the feed gas water content and number of towers
    Figure 5. Desiccant height vs the feed gas content and number of towers
    Figure 5. Desiccant height vs the feed gas content and number of towers

    Figure 6 shows the variation of the regeneration gas requirement with the feed gas water content while maintaining constant heating and cooling times.  This figure indicates that higher feed gas water content require more regeneration gas. However, the required regeneration gas rate is practically the same for the two-tower and three-tower systems under evaluation.

    Figure 6. Percent of % feed gas for regeneration vs the feed gas water content and number of towers
    Figure 6. Percent of % feed gas for regeneration vs the feed gas water content and number of towers

    Similarly, Figures 7, 8, and 9 show the variation of the required heating load, cooling load, and regeneration gas heater load, with the feed gas water content and number of towers. Since the towers in the two-tower system are larger than the towers in the three-tower system, the heating and cooling loads are larger for the towers in the two-tower systems (figures 7 and 8). However, Figure 9 indicates that the required regeneration heater loads are almost the same for the two-tower and three-tower systems.

    Figure 7. Heating load vs the feed gas water content and number of towers
    Figure 7. Heating load vs the feed gas water content and number of towers
    Figure 8. Cooling load vs the feed gas water content and number of towers
    Figure 8. Cooling load vs the feed gas water content and number of towers
    Figure 9. Regeneration gas heater load vs the feed gas water content and number of towers
    Figure 9. Regeneration gas heater load vs the feed gas water content and number of towers

    Summary:

    Water content of the feed gas affected by temperature, pressure and flow rate is the key factor in sizing and operation of adsorption dehydration system. Higher water load requires a larger size bed, higher heating and cooling loads and higher rate of regeneration gas.

    In the cases evaluated in this Tip of the Month, the tower diameters in the two-tower systems are larger but their heights are shorter. The mass of desiccants and mass of steel per tower in the two-tower system are larger than in the three-tower system. Therefore the heating and cooling loads are larger for the towers in the two-tower system. The regeneration gas heater loads are almost the same for two and three-tower system.

    To learn more about similar cases and how to minimize operational problems, we suggest attending our PF49 (Troubleshooting Oil and Gas Facilities), PF42 (Separation Equipment Selection and Sizing), G4 (Gas Conditioning and Processing), G5 (Gas Conditioning and Processing – Special), and PF4 (Oil Production and Processing Facilities), courses.

    PetroSkills offers consulting expertise on this subject and many others. For more information about these services, visit our website at http://petroskills.com/consulting, or email us at consulting@PetroSkills.com.

    Dr. Mahmood Moshfeghian

    References:

    1. Campbell, J.M., Gas Conditioning and Processing, Volume 2: The Equipment Modules, 9th Edition, 2nd Printing, Editors Hubbard, R. and Snow–McGregor, K., Campbell Petroleum Series, Norman, Oklahoma, 2014.
    2. Malino, H. M., http://www.jmcampbell.com/tip-of-the-month/2015/05/benefits-of-standby-time-in-adsorption-dehydration-process/
    3. Moshfeghian, M., http://www.jmcampbell.com/tip-of-the-month/2015/10/what-is-the-impact-of-feed-gas-conditions-on-the-adsorption-dehydration-system/
    4. GCAP Version 9.1.1, Gas Conditioning and Processing Software, Editor Moshfeghian, M., PetroSkills/Campbell, Norman, Oklahoma, 2015.
  • Gas-Liquid Separators Sizing Parameter

    Download PDF Version

    In the December 2014 tip of the month (TOTM) [1], we discussed troubleshooting of gas-liquid separators for removal of liquids from the gas stream leaving the separator. There are two methods for sizing gas-liquid separators: 1. Droplet settling theory method, 2. Souders-Brown approach. Historically the Souders-Brown equation has been employed as it can provide reasonable results and is easy to use, but has shortcomings in terms of quantifying separator performance.  References [2-4] provide the details on the droplet settling theory methods which can be used to more accurately quantify separator performance.  The Souders-Brown method is limited in that it is based on the average droplet size, but cannot quantify the amount of liquid droplets exiting the gas gravity section.

    In this TOTM, we will focus on the application of Souders-Brown approach in gas-liquid separators and present diagram, simple correlations and tables to estimate the Souders-Brown equation constant, KS (the so called sizing parameter). We will consider both vertical and horizontal gas-liquid separators. Knowing the actual gas flow rate through the vessel, one can use KS parameter to determine the maximum allowable gas velocity through the vessel and determine the required separator diameter. One can also use the appropriate value of KS to size the mist extractor in the vessel. The performance of a gas-liquid separator is highly dependent on the value of KS; therefore, the choice of appropriate  KS –values is important.

    Gas Gravity Separation Section

    The gas gravity separation section of a separator has two main functions:

    1. Reduction of entrained liquid load not removed by the inlet device
    2. Improvement / straightening of the gas velocity profile.

    Most mist extractors have limitations on the amount of entrained liquid droplets that can be efficiently removed from the gas, thus the importance of the gas gravity section to remove the liquids to an acceptable level upstream of the mist extractor.  This is particularly important for separators handling higher liquid loads. For scrubber applications with low liquid loadings,  the KS –values will be primarily dependent on the mist extractor type, and the gas gravity separation section becomes less important. For the higher liquid load applications, i.e. conventional separators, there are two approaches for sizing the gravity separation section to remove liquid droplets from the gas:

    1. The Souders-Brown approach (Ks Method)
    2. Droplet settling theory

    The Souders-Brown Approach

    If we consider a spherical liquid droplet with a diameter of DP in the gas phase two forces as shown in Figure 1 act on it. The drag force, FD, is exerted by flow of gas and gravity force, FG, is exerted by the weight of droplet. The drag force acts to entrain the liquid droplet while the gravity force acts to pull it down and separating it from the gas phase.

    Figure 1. Schematic of the forces acting on a liquid droplet in the gas phase [5]
    Figure 1. Schematic of the forces acting on a liquid droplet in the gas phase [5]
    Assuming plug flow with no eddies or disturbances, a single droplet and ignoring the end effect, at equilibrium (free fall or terminal velocity), these two forces are equal.

    Fd=Fg       (1)

    As presented in the Appendix, substitution of expressions for the drag and gravity forces in Equation 1, the maximum allowable gas velocity, VGmax, which prevents entrainment of liquid is obtained.

    eq2

    Equation 2 is called Souders-Brown [6] equation and KS is referred to as the design or sizing parameter. The terms ρG and ρL are the gas phase and liquid phase densities, respectively.

    Once the maximum gas velocity, VGmax, through the vessel is determined by Equation 2, one can calculate the minimum vessel diameter, Dmin by Equation 3.

    eq3

    Where:

    FG        = Fraction of cross section area available for gas flow (FG = 1 for vertical separators and is a function of liquid height for horizontal separators)

    qa         = Gas flow rate at the actual flowing condition

    The Design Parameter, KS

    The design parameter, KS, in the Souders-Brown equation is an empirical parameter and is a key factor for the sizing the gas-liquid separators vessel diameter as well as for determination of the mist extractor diameter. Its value depends on several factors including:

    • Pressure
    • Fluid properties (note temperature has a large impact on fluid properties)
    • Separator geometry
      • Vessel length and the liquid level (for horizontal separators)
    • Steadiness of flow
    • Inlet device design and performance
    • Relative amounts of gas and liquid
    • Most importantly – mist extractor type and design (e.g. mesh pad, vane pack, multi–cyclone)

    There are several sources that one can look up the KS –values for different applications.  In the following sections, we will discuss three sources.

    A. API 12 J

    The API 12J [7] recommends ranges of KS –values for vertical and horizontal gas-liquid separators. These values are presented in Table 1. The equivalent of API 12J for the North Sea region is NORSOK P-100.

    Table 1. API 12 J recommended range of KS –values for vertical and horizontal separators [7]

    tab1-2

    Per API 12J, “the maximum allowable superficial velocity, calculated form the above factors, is for separators normally having a wire mesh mist extractor. This rate should allow all liquid droplets larger than 10 microns to settle out of the gas. The maximum allowable superficial velocity or other design criteria should be considered for other type mist extractor. Mist extractor manufacturer’s recommended minimum distances upstream and downstream of the wire mesh between gas inlet and outlet nozzles should be provided for full utilization of the mist extractor. These values assume separators are equipped with standard mesh pad mist extractors” [7].

    B. Campbell Book

    The Ks method, Equation 2, is an empirical approach to estimate the maximum allowable gas velocity to achieve a desired droplet separation. For vertical separators with no mist extractor devices, Chap 11, Vol 2 of the Gas Conditioning and Processing book presents KS as a function of pressure and liquid droplet size [5]. This dependency of KS on pressure and droplet size is presented in Figure 2 [5]. Note for each droplet size a range of KS –values are given for a specified pressure.

    For horizontal separators, the sizing depends on (in addition to the droplet size, density of gas and liquid phases, and gas velocity) separator effective length, Le, and the depth available for gas flow, hG, (i.e. liquid level) in the separators.

    Figure 2. KS as a function of pressure and liquid droplet size for vertical separators with no mist extractor devices [5]
    Figure 2. KS as a function of pressure and liquid droplet size for vertical separators with no mist extractor devices [5]

    Sizing of the horizontal separators are more complicated. Referring to Figure 3, the effective Le may be defined in terms of separator actual length and diameter like Le=L-D. Therefore, the Souders-Brown parameter for horizontal separators, KSH, can be estimated in by Equation 4 in terms of KSV (read from Figure 2) for vertical separator [3].

    eq4

    If the calculated value of KSH by Equation 4 is greater than the maximum value of 0.7 ft/sec (0.21 m/s), it should be set equal to this maximum value.

    Figure 3. Schematic of a horizontal gas-liquid separator [5]
    Figure 3. Schematic of a horizontal gas-liquid separator [5]
    The horizontal separator sizing is a trial-and-error procedure. Normally, the Le/D and hg/D (or hL/D) are assumed and KSH, Vgmax, D are calculated by Equations 4, 2, and 3, respectively. The effective length and actual lengths are calculated by Equation 5.

    eq5

    Where:

    D         = Diameter

    FL        = Fraction of cross section area occupied by liquid (function of liquid height in horizontal separator)

    qL         = Liquid actual volume flow rate

    t           = Residence time per API 12J [7]

    If the calculated L/D is outside of its recommended range (normally 3 < L/D < 6), the liquid height in the vessel is changed and the calculations are repeated. For detail of calculations procedure refer to chapter 11 of reference [5].

    C. KS Correlations

    The curves for different droplet sizes shown in Figure 2 are fitted to a 3rd order polynomial (for droplet sizes of 100, 150, and 300 microns). The correlation is in the form of Equation 6 and its regressed coefficients a, b, c, and d are presented in Tables 2A and 2B for field (FPS) and System International (SI) units, respectively.

    eq6

    In Table 2, each droplet size in micron (µ) is preceded by letter L or U representing the lower and upper curve, respectively. The pressure is in psi and KS is in ft/sec for FPS (kPa and m/s in SI).

    The last row of Table 2 provides the average absolute percent deviation (AAPD) of the predicted KS by the proposed correlation from the corresponding values of Figure 2 values.

    Table 2A (FPS). Regressed coefficients for Equation 6 (P in psi and KS in ft/sec)
    Droplet Size: 100 – 300 microns

    tab2a

    Table 2B (SI). Regressed coefficients for Equation 6 (P in kPa and KS in m/s in SI).
    Droplet Size: 100 – 300 microns

    tab2b

    The two curves for 500 micron droplet size in Figure 2 were divided into 4 and 2 segments based on pressure range for the lower and upper curves, respectively. Each segment was fitted to a linear equation in the form of Equation 7 and its regressed coefficients e and f are presented in Tables 3A and 3B for FPS and SI units, respectively.

    eq7

    Table 3A (FPS). Regressed coefficients for Equation 7 (P in psi and KS in ft/sec)
    Droplet Size: 500 microns

    tab3a

    Table 3B (SI). Regressed coefficients for Equation 7 (P in kPa and KS in m/s in SI).
    Droplet Size: 500 microns

    tab3b-2

    D. Mist Extractors

    The mist extractor is the final gas cleaning device in a conventional separator.  The selection, and design to a large degree, determine the amount of liquid carryover remaining in the gas phase.  The most common types include wire mesh pads (“mesh pads”), vane-type (vane “packs”) and axial flow demisting cyclones.  Figure 4 shows the location and function of a typical mist extractor in a vertical separator.

    Mist extractor capacity is defined by the gas velocity at which re-entrainment of the liquid collected in the device becomes appreciable.  This is typically characterized by a KS –value, as shown in Equation 2.  Mesh pads are the most common type of mist extractors used in vertical separator applications.  The primary separation mechanism is liquid impingement onto the wires, followed by coalescence into droplets large enough to disengage from the mesh pad.  References [1-5] provide mesh pad examples.  Table 4 provides a summary of mesh pad characteristics and performance parameters.

    Figure 4. Typical mist extractor in a vertical separator [5]
    Figure 4. Typical mist extractor in a vertical separator [5]

    Table 4. Mesh pads KS and performance parameters [3, 5, 8]

    tab4

    Notes:

    • Flow direction is vertical (upflow).
    • Assume mesh pad KS – value decline with pressure as shown in Table 5. Table 5 was originally developed for mesh pads, but is used as an approximation for other mist extractor types. [9].
    • If liquid loads reaching the mesh pad exceed the values given in Table 4, assume capacity (KS) decreases by 10% per 42 L/min/m2 (1 gal/min/ft2). [2-4].
    • These parameters are approximate.

    Table 5. Mesh pad KS deration factors as a function of pressure [5]

    tab5

    Vane packs, like mesh pads, capture droplets primarily by inertial impaction.  The vane bend angles force the gas to change direction while the higher density liquid droplets tend to travel in a straight-line path, and impact the surface of the vane where they are collected and removed from the gas flow. Table 6 provides vane pack performance characteristics [3, 5, 8].

    In the case of demisting cyclones, the vendor should be consulted in regards to performance for the current operations of interest.

    Table 6. Typical vane-pack characteristics [3, 5, 8]

    tab6

    Notes:

    1. Assume vane-pack KS – value decline with pressure as shown in Table 5.
    2. If liquid loads reaching the vane pack exceed the values given in Table 2, assume capacity KS decreases by 10% per 42 L/min/m2 (1 gal/min/ft2). [2-4].
    3. These parameters are approximate only. The vane-pack manufacturer should be contacted for specific information.

    Conclusions: 

    We focused on the application of Souders-Brown Approach (SBA) in gas-liquid separators and presented diagram, simple correlations and tables to estimate the SBA design parameter, KS.

    • The SBA can provide reasonable results and is easy to use.
    • The SBA is limited in that it is based on the average droplet size, but cannot quantify the amount of liquid droplets exiting the gas gravity section and mist extractor section.
    • In a future TOTM we will discuss the droplet settling theory methods which can be used to more accurately quantify separator performance.
    • Sizing of three-phase gas-liquid hydrocarbon-liquid water separators are more complicated and will be discussed in another TOTM.

    To learn more about similar cases and how to minimize operational problems, we suggest attending our PF49 (Troubleshooting Oil and Gas Facilities), PF42 (Separation Equipment Selection and Sizing), G4 (Gas Conditioning and Processing), G5 (Gas Conditioning and Processing – Special), and PF4 (Oil Production and Processing Facilities), courses.

    PetroSkills offers consulting expertise on this subject and many others. For more information about these services, visit our website at http://petroskills.com/consulting, or email us at consulting@PetroSkills.com.

    Dr. Mahmood Moshfeghian

    References:

    1. Snow–McGregor, K., http://www.jmcampbell.com/tip-of-the-month/2014/12/troubleshooting-gas-liquid-separators-removal-of-liquids-from-the-gas/
    2. Bothamley, M., “Gas-Liquid Separators – Quantifying Separation Performance Part 1,” SPE Oil and Gas Facilities, pp. 22 – 29, Aug. 2013.
    3. Bothamley, M., “Gas-Liquid Separators – Quantifying Separation Performance Part 2,” SPE Oil and Gas Facilities, pp. 35 – 47, Oct. 2013.
    4. Bothamley, M., “Gas-Liquid Separators – Quantifying Separation Performance Part 3,” SPE Oil and Gas Facilities, pp. 34 – 47, Dec. 2013.
    5. Campbell, J.M., Gas Conditioning and Processing, Volume 2: The Equipment Modules, 9th Edition, 2nd Printing, Editors Hubbard, R. and Snow–McGregor, K., Campbell Petroleum Series, Norman, Oklahoma, 2014.
    6. Souders, M. and Brown, G. G., “Design of Fractionating Columns-Entrainment and Capacity,” Industrial and Engineering Chemistry, Volume 26, Issue 1, p 98-103, 1934.
    7. American Petroleum Institute, 12J, Specification for Oil and Gas Separators, 8th Edition, October, 2008.
    8. PF-49, Troubleshooting Oil and Gas Processing Facilities, Bothamley, M., 2014, © PetroSkills, LLC. All Rights reserved.
    9. Fabian, P., Cusack, R., Hennessey, P., Neuman, M., “Demystifying the Selection of Mist Eliminators, Part 1: The Basics,” Chem Eng 11 (11), pp. 148 – 156, 1993.

     

    Appendix

    Derivation of the Souders-Brown and Stokes’ Law Equations

    If we consider a spherical liquid droplet with a diameter of, DP, in the gas phase two forces as shown in Figure 1 act on it. The drag force, FD, is exerted by flow of gas and gravity force, FG, is exerted by weight of droplet. The drag force acts to entrain the liquid droplet while the gravity force acts to pull it down.

    Figure 1. Schematic of the forces acting on a liquid droplet in the gas phase [5]
    Figure 1. Schematic of the forces acting on a liquid droplet in the gas phase [5]

    At equilibrium, these two forces are equal.

    eqFdFg

    The drag force is expressed as:

    eqFd

    The droplet projected area, AP, is defined by:

    eqAp

    The gravity force, FG, is defined

    eqFg

    The volume of spherical droplet, VD, is calculated by

    eqVp

    Substitution of Equations 3 and 4 into Equation 1 and solving for the gas maximum velocity,

    eqVgmax

    For practical applications, the first term on the right hand side is replaced by KS

    eqKs

    Therefore, the maximum gas velocity which prevents entrainment of liquid is obtained.

    eqVgmaxKs

    Equation 6 is called Souder-Brown equation and KS is referred to as the design parameter.

    Where:

    AP        = Project area of droplet

    CD        = Drag coefficient

    g          = Acceleration of gravity

    gC        = Conversion factor

    V          = Gas velocity

    VP        = Volume of droplet

    ρG        = Gas density

    ρL         = Liquid density

    Once the maximum, VGmax, gas velocity through the vessel is determined by Equation 6, one can calculate the required minimum cross sectional area of vessel for gas flow by the following equation.

    eqAgmin

    Solving for the minimum vessel diameter, Dmin.

    eqDmin

    Where:

    FG        = Fraction of cross sectional area available for gas flow (FG = 1 for vertical separators and it is a function of liquid height for horizontal separators)

    qa         = Gas flow rate at the actual flowing condition

    The drag coefficient, CD, is a function of Reynolds number, Re=(DPG)/µG. For Stokes’ law Re <≈2

    eqCd

    Substitution of CD from Equation 9 into Equation 4 gives liquid droplet terminal velocity, VT, in the gas phase based on the Stokes’ law.

    eqVt

    Similarly, the terminal velocity for the other flow regimes like Intermediate and Newton can be derived based on their corresponding expressions for the drag coefficients [3].

    Download PDF Version

  • Effect of Relative Density (Specific Gravity) on the Saturated Water Content of Sweet Natural Gases

    In the past Tips of the Month (TOTM), we discussed the phase behavior and water content of lean sweet, sour natural gases and acid gases–water systems. Specifically, in the November 2007 [1], February 2014 [2], and September 2014 [3] Tips of the Month (TOTM), we discussed the phase behavior of water-saturated sour gases and acid gases. We also demonstrated the accuracy of shortcut and rigorous calculation methods. In the April 2015 TOTM [4] we introduced correlation and graphs to estimate the water content of sour gases.

    In this TOTM, we will study the effect of relative density (Specific Gravity, SG) on the saturated water content of sweet natural gases. The results of this study include the water content of sweet natural gases as a function of relative density in the range of 0.60 to 0.80.  Four temperatures of 4.4, 23.9, 37.8 and 149 °C (40, 75, 100, and 300 °F) were considered. For each temperature, the saturated water content was calculated for pressures of 1724, 3448, 6897, and 13 793 kPaa (250, 500, 1000 and 2000 psia).

    Water Content Calculation

    The dry gas compositions of the four mixtures studied in this study are presented in Table 1. The Soave-Redlich-Kwong equation of state (SRK EOS) [5] in the ProMax [6] was used to predict the water content of these gas mixtures at different pressures and temperatures. A simplified process diagram for this study is presented in Figure 1. The feed dry-gas at the specified pressure and temperature was saturated with water first and passed through a separator. The water content and specific gravity of vapor leaving the separator were recoded.

    Results and Discussion

    Figure 2 through 5 present the saturated water content of sweet natural gases as a function of the dry gas  relative density and pressures of 1724 to 13793 kPaa (250 to 2000 psia) for temperatures of 4.4, 23.9, 37.8 and 149 °C (40, 75, 100, and 300 °F), respectively.

    Table 1. The composition of relative density (Specific Gravity, SG) of the four gas mixture

    tab1

    Figure 1. Simplified process flow diagram
    Figure 1. Simplified process flow diagram
    Figure 2. Variation of saturated water content of sweet natural gas with the dry gas relative density and pressure at 4.4 °C (40 °F)
    Figure 2. Variation of saturated water content of sweet natural gas with the dry gas relative density and pressure at 4.4 °C (40 °F)
    Figure 3. Variation of saturated water content of sweet natural gas with the dry gas relative density and pressure at 23.9 °C (75 °F)
    Figure 3. Variation of saturated water content of sweet natural gas with the dry gas relative density and pressure at 23.9 °C (75 °F)
    Figure 4. Variation of saturated water content of sweet natural gas with the dry gas relative density and pressure at 37.8 °C (100 °F)
    Figure 4. Variation of saturated water content of sweet natural gas with the dry gas relative density and pressure at 37.8 °C (100 °F)
    Figure 5. Variation of saturated water content of sweet natural gas with the dry gas relative density and pressure at 149 °C (300 °F)
    Figure 5. Variation of saturated water content of sweet natural gas with the dry gas relative density and pressure at 149 °C (300 °F)

    While the results presented in the above diagrams are in agreement with Figure 6.1 of Reference [7], it is not in agreement with the suggested correction factor in the inset of Figure 20-3 of the GPSA data book [8] and the results of Reference [9]. For the gases with relative density of 0.6 up to 0.8, the GPSA figure indicates a correction factor of 1 to about 0.97 should be multiplied by the water content of saturated water content of sweet gas with relative density of 0.6. This study indicates that for the range of 0.6 to 0.8 relative density, the water content is not a function of the hydrocarbon gas composition. To verify the result of this study, the saturated water content of several gas mixtures were compared with the results of EzThermo [10] software. The SRK EOS in the EzThermo software was developed and regressed to predict the properties of sweet synthetic natural gas and natural gas compositions [11]. The same is also valid for the SRK EOS in ProMax. The comparison of these two software is presented in Table 2. This table indicates an excellent agreement between these two software.

    Table 2. Comparison of the ProMax [6] and EzThermo [10] predicted saturated water content results at 37.8 °C (100 °F)

    tab2

    Conclusions

    Figures 2 through 5 cover wide ranges of pressures and temperatures commonly encountered in the gas processing operation. The analysis of Figures 2 through 5 indicates that the gas relative density has minimal effect on the saturated water content of sweet natural gases. This conclusion is valid for the following ranges:

    1. Sweet gas relative density in the range of 0.6 to 0.8
    2. Temperature range of 4.4 to 149 °C (40 to 300 °F)
    3. Pressure range of 1724 to 13793 kPaa (250 to 2000 psia)

    Additional experimental water content data are being taken and analyzed by the GPA Research Committee to update the relative density correction factor as presented in [8].  From the data that are currently available, it appears that this correction is minimal if any for gas compositions, and temperature and pressure ranges that typically occur in the oil and gas facilities. This will be confirmed once the results from the GPA Research Committee have been published.

    To learn more about similar cases and how to minimize operational problems, we suggest attending our G4 (Gas Conditioning and Processing), G5 (Gas Conditioning and Processing – Special), G6 (Gas Treating and Sulfur Recovery), PF49 (Troubleshooting Oil and Gas Facilities),  and PF81 (CO2 Surface Facilities) courses.

    PetroSkills offers consulting expertise on this subject and many others. For more information about these services, visit our website at http://petroskills.com/consulting, or email us at consulting@PetroSkills.com.

    Dr. Mahmood Moshfeghian

    Reference:

    1. http://www.jmcampbell.com/tip-of-the-month/2014/09/lean-sweet-natural-gas-water-content-correlation/
    2. http://www.jmcampbell.com/tip-of-the-month/2007/11/water-sour-natural-gas-phase-behavior/
    3. http://www.jmcampbell.com/tip-of-the-month/2014/09/lean-sweet-natural-gas-water-content-correlation/
    4. http://www.jmcampbell.com/tip-of-the-month/2014/02/acid-gas-water-content/
    5. Soave, G., Chem. Eng. Sci., Vol. 27, pp. 1197-1203, 1972.
    6. ProMax 3.2, Bryan Research and Engineering, Inc., Bryan, Texas, 2014.
    7. Campbell, J.M., “Gas conditioning and Processing, Vol. 1: The Basic Principles”, 9th Edition, 2nd Printing, Editors Hubbard, R. and Snow–McGregor, K., Campbell Petroleum Series, Norman, Oklahoma, 2014.
    8. GPSA Engineering Data Book, Section 20, Volume 2, 13th Edition, Gas Processors and Suppliers Association, Tulsa, Oklahoma, 2012.
    9. Maddox, R.N., L.L. Lilly, M. Moshfeghian, and E. Elizondo, “Estimating Water Content of Sour Natural Gas Mixtures”, Laurence Reid Gas Conditioning Conference, Norman, OK, Mar., 1988.
    10. EzThermo, Moshfeghian, M. and R.N. Maddox, 2015.
    11. GPA Research Report RR-42, Predicting Synthetic Gas and Natural Gas Thermodynamic Properties Using a Modified Soave Redlich Kwong Equation of State, Oklahoma State University, August 1980
  • Effect of Chemical Additive on Crude Oil Pipeline Pressure Drop

    For transportation of crude oil, the pumping power requirement varies as the crude oil viscosity changes. Increasing °API or line average temperature reduces the crude oil viscosity. The reduction of viscosity results in higher Reynolds number, lower friction factor and in effect, lower pumping power requirements.

    In the March 2009 tip of the month (TOTM), procedures for calculation of friction losses in oil and gas pipelines were presented. The sensitivity of friction pressure drop with the wall roughness factor was also demonstrated. In the August 2009 TOTM, we also demonstrated the effect crude oil °API and the pipeline average temperature on the pumping requirement.

    In practical situations, an originating station takes crude out of storage and the midline stations taking suction from the upstream section of pipeline. Oil in the tank is often at ambient temperature, whereas once in the pipeline, the oil cools (or warms) to the same temperature as the ground.  In some parts of the world, the tank might be at +38 °C (+100 °F).  The first midline pumping station could operate at 18 °C (65 °F), and all subsequent pumping stations might operate at ground temperature, or notionally 9 °C (48 °F) with some seasonal variation. Therefore, a sound pipeline design should consider expected variation in crude oil viscosity which is normally a function of crude oil °API, and the line average temperature. In addition to °API and temperature, chemical additives may also affect crude oil pipeline pressure drop.

    To reduce pressure drop and increase pipeline capacity, oil industry has utilized drag reducing agents. Drag-reducing agents, or drag-reducing polymers, are additives in pipelines that reduce turbulence in a pipe. Usually used in petroleum pipelines, they increase the pipeline capacity by reducing turbulence and therefore allowing the oil to flow more efficiently [1]. In addition to drag reducing agents, another group of chemicals called “Incorporative Additives”, which reduces crude oil viscosity, may be used. Halloran presented a series of general reading articles on chemical additives [2-4].

    In this TOTM, we will demonstrate the effect of an incorporative additive on crude oil viscosity and consequently on pressure drop for crude oil pipeline transportation.

    Case Study: Part 1 – Viscosity Reduction

    The laboratory measured kinematic viscosity for different °API crude oil samples without and with “Incorporative Additive” at 50 °C (122 °F) reported by Oil Flux Americas [6] are shown in Table 1. The calculated density, absolute viscosity and percent reduction in viscosity for each oil sample at 50 °C (122 °F) are also shown in this table. As noted in this table, the lower °API (heavier oil), the greater the reduction in oil viscosity. The measured kinematic viscosities as a function of crude oil °API are shown in Figure 1. The absolute viscosity is calculated by multiplying the measured kinematic viscosity by density. The corresponding calculated absolute viscosities are also shown in Table 1 for crude oil samples “Without” and “With” additive, respectively.

    Table 1. Measured kinematic viscosity [6] and absolute viscosity for several crude oil samples at 50 °C (122 °F) without and with chemical additive.

    table1
    cSt = (mm)2/s cP = Poise/100 = Pa.s/1000 =kg/m-s/1000 = 0.000672 lbm/ft-sec * Used in the case studies
    Figure 1. Effect of chemical additive on crude oil absolute viscosity at 50°C (122 °F)
    Figure 1. Effect of chemical additive on crude oil absolute viscosity at 50°C (122 °F)

    The absolute viscosities (μ) at 50 °C (122 °F) are fitted to a quadratic equation as follows:

    equations1

    The absolute viscosities and the fitted correlations are shown in Figures 2 and 3 for crude oil samples “without” and “with” chemical additive, respectively.

    Case Study: Part 2 – Pressure Drop Calculations

    For a case study, we will consider a 55 km (34.18 miles) pipeline with an outside diameter of 406.4 mm (16 in) carrying crude oil with two separate flow rates of 7,950 and 15,900 m3/d (50,000 and 100,000 bbl/day). The wall thickness was estimated to be 5.7 mm (0.225 in). The wall roughness is 46 microns (0.0018 in) or a relative roughness (ε/D) of 0.0001. The procedures outlined in the March 2009 TOTM were used to calculate the line pressure drop due to friction. Since the objective is to study the effect of incorporative chemical additive, we will ignore elevation change.

    It is also assumed the line temperature is constant at 50 °C (122 °F). The change in pressure drop (ΔP) due to changes in crude oil viscosity for this case study will be calculated and presented in the following sections.

    Figure 2. Measured absolute viscosity at 50°C (122 °F) for crude oils without chemical
    Figure 2. Measured absolute viscosity at 50°C (122 °F) for crude oils without chemical
    Figure 3. Measured absolute viscosity at 50°C (122 °F) for crude oils with chemical
    Figure 3. Measured absolute viscosity at 50°C (122 °F) for crude oils with chemical

    Tables 2a and 2b show the calculated pressure drop for four different crude oils with varying viscosities in System International (SI) and field units (FPS – Foot, Pound and Second), respectively. The measured absolute viscosities were used to calculate pressure drops for all cases. The oil flow rate is 7,950 m3/d (50,000 bbl/day). Table 2 indicates that by using incorporative additives a reduction of up to 24% in pressure drop is achieved for this case study. The results in this table also indicate that the percent reduction in pressure drop (0.5% for the lightest oil) is not as high as the percent reduction in viscosity (3.7% for the lightest oil). Another observation is that the reduction in viscosity and consequently in pressure drop for light crudes oils is not as significant as for the heavy crudes.

    Table 2a. Pipeline pressure drop for four different crude oils without and with additive at 50 °C and oil flow rate of 7,950 m3/d

    table2a

    Table 2b. Pipeline pressure drop for four different crude oils without and with additive at 122 °F and oil flow rate of 50,000 bbl/day

    table2b

    Table 2c. Reynolds number and friction factor for the cases in Table 2 a and b

    table2c

    Similarly, for an oil flow rate of 15,900 m3/d (100,000 bbl/day), Tables 3a and 3b show the calculated pressure drop for the same four crude oils with varying viscosities. Tables 3a and 3b indicate that as flow rates are increased, less reduction in pressure drop is obtainable if the flow becomes turbulent. For the case of 16.4 °API, the reduction in pressure drop is 6.6% compared to 20.7 reduction when the flow rate was of 7,950 m3/d (50,000 bbl/day). The calculated Reynolds number, Moody friction factors for the cases of lower and higher oil flow rates are shown in Tables 2c and 3c, respectively.

    Table 3a. Pipeline pressure drop for four different crude oils without and with additive at 50 °C and oil flow rate of 15,900 m3/d

    table3a

    Table 3b. Pipeline pressure drop for four different crude oils without and with additives at 122 °F and oil flow rate of 100,000 bbl/day

    table3b

    Table 3c. Reynolds number and friction factor for the cases in Table 3 a and b

    table3c

    In order to show the impact of chemical on pipeline capacity for the same pressure drop, let’s consider the heavy crude oil with 12.7 °API. As shown in Table 2a and 2b for an oil flow rate of 7,950 m3/d (50,000 bbl/day) the pressure drop without chemical was 4.684 MPa (679 psia). For the same pressure drop and using the reduced viscosity due to addition of chemicals, the capacity increases to 10,472 m3/d (65,865 bbl/day). This is equivalent of 31% increase in pipeline capacity. Similarly, referring to Tables 3a and b for an oil flow rate of 15,900 m3/d (100,000 bbl/day) for the case of without chemical, the pressure drop was 9.367 MPa (1358 psia). The calculated capacity for the same pressure drop is 20943 m3/d (131,730 bbl/day). Again, a 31 % increase in pipeline capacity is observed.

    Conclusions:

    The following conclusions can be made based on this case study:

    1. The mechanisms of how drag reducing agents work are different from incorporative chemical additives. Incorporative chemical additives reduce viscosity.
    2. Utilizing incorporate chemical additives can reduce crude oil viscosity and consequently reduces the pipeline pressure drop significantly. For existing pipelines this means an increase in the capacity of the line and/or reduction in pump power requirement.
    3. The reduction of viscosity and pressure drop are more significant for heavier crude oils. As the oil gets lighter the effect of chemical additives is diminished. At lower temperatures the oil viscosity increases; therefore, the effect of chemical additives may become more significant for lighter crude oils, too.
    4. The percent reduction in pipeline pressure drop is not always as large as the percent reduction for viscosity.
    5. The incorporative chemical additives are most effective for laminar flow and/or heavier crude oils.
    6. A total cost analysis based on hydraulic design and chemical additives with consideration for HSE (health, safety, and environment) should be made for effective design and operation.

      To learn more about similar cases and how to minimize operational problems, we suggest attending our G4 (Gas Conditioning and Processing), PF4 (Oil Production and Processing Facilities), PL22 (Pipeline Systems Overview) and PL42 (Onshore Pipeline Facilities – Design, Construction and Operations), courses.

      PetroSkills offers consulting expertise on this subject and many others. For more information about these services, visit our website at http://petroskills.com/consulting, or email us at consulting@PetroSkills.com.

      Dr. Mahmood Moshfeghian

      References:

      1. Wikipedia, http://en.wikipedia.org/wiki/Drag_reducing_agent, 2015
      2. Halloran, M.D., “Taming Crude Behavior: Understanding production

        Additives – Part 1”, PennEnergy, Oil & Gas, September 22, 2014

      3. Halloran, M.D., “Taming Crude Behavior: Understanding production

        Additives – Part 2”, PennEnergy, Oil & Gas, September 24, 2014

      4. Halloran, M.D., “Taming Crude Behavior: Understanding production

        Additives – Part 3”, PennEnergy, Oil & Gas, September 26, 2014

      5. Halloran, M.D., “Incorporative Production Additives Lower HSE Concerns &

        Improve Processes”, Upstream Pumping-Wellhead Technology & Services, January/February 2015, http://upstreampumping.com/article/2015/incorporative- production-additives-lower-hse-concerns-improve-processes/

      6. Oil Flux Americas, LLC, www.oilfluxamericas.com, 2015

  • Benefits of Standby Time in Adsorption Dehydration Process

    Molecular sieves are used upstream of turboexpander units and LNG facilities to dehydrate natural gas to <0.1 ppmv water content.   In the natural gas industry, the molecular sieves employ heat to drive off the adsorbed water.  Figure 1 shows a typical flow schematic for a 2 tower system; Figure 2 shows a 3 tower system.

    Figure 1. Typical process flow diagram for a 2-tower adsorption dehydration system [1]
    Figure 1. Typical process flow diagram for a 2-tower adsorption dehydration system [1]
    The cyclical heating/cooling of the adsorbent results in a capacity (kg water/100 kg adsorbent; lbm water/100 lbm adsorbent) decline due to a gradual loss of crystalline structure and/or pore closure.  A more troublesome cause of capacity decline is contamination of the molecular sieves due to liquid carryover from the upstream separation equipment.

    Figure 3 shows a generic molecular sieve capacity decline curve.  A few important observations can be made from Figure 3:

    1. The life of the adsorbent is a function of the number of cycles, not the elapsed calendar time.
    2. The capacity decline is steep at the beginning but gradually flattens out. This assumes no step-change events such as NGL, glycol, and/or liquid amine carryover, bed support failure, etc.
    3. Shown in this figure are “Good”, “Average” and “Poor” curves that are a function of site specific factors.
    4. Locating one data point on Figure 3 from a performance test allows you to extrapolate the decline curve of the unit in question.

    If your regeneration circuit has excess capacity over the “normal design conditions”, i.e., a design factor, you have standby time.   This excess capacity allows you to reduce your online adsorption time and “turn the beds around” faster by regenerating the beds in a shorter cycle time.   When you are involved in the design of an adsorption unit, it is recommended to add 10 – 20% excess regeneration capacity.

    Because of the capacity decline curves flatten out, available standby time may be able to extend the life of a molecular sieve unit when your unit is operating on fixed cycle times.  Other operating options include: running each cycle to water breakthrough; and, reducing the cycle times in discreet steps throughout the life of the adsorbent.

     

    Figure 2. Typical process flow diagram for a 3-tower adsorption dehydration system [1]
    Figure 2. Typical process flow diagram for a 3-tower adsorption dehydration system [1]
    Figure 3. A generic molecular sieve decline curves [1]
    Figure 3. A generic molecular sieve decline curves [1]
    To illustrate the benefits of standby time, consider the following case study.   A natural gas processing plant has commissioned a new 3 tower molecular sieve dehydration unit to process 11.3 x 106 std m3/d (400 MMscfd) prior to flowing to a deep ethane recovery unit.  The unit is expected to run for 3 years before needing a recharge and the plant turnaround is based on this expectation.  The following assumptions are made:

    • 3 tower system (2 towers on adsorption, 1 on regeneration)
    • External Insulation
    • Tower ID = 2.9 m (9.5ft)
    • Each tower contains 24630 kg [54300 lbm] of Type 4A 4×8 mesh beads
    • Regeneration circuit capable of handling an extra 15% of flow
    • Unit is operated on fixed time cycles
    • No step-change events such as liquid carryover, poor flow distribution, etc.

    The design basis and molecular sieve design summary are shown in Tables 1 and 2.  The additional 15% of flow from the regeneration gas heater is well below the point at which bed lifting will occur.

    Table 1. Design basis for the case study
    Table 1. Design basis for the case study
    Table 2. Design Summary  for the Case Study
    Table 2. Design Summary for the Case Study

    The calculations presented here are valid for low pressure regeneration (less than 4100 kPaa (600 psia).  Using the concepts outlined in Chapter 18 of Gas Conditioning and Processing, Volume 2 [1]:  The Equipment Modules (9th Edition) we find a design life factor, FL, of 0.6 after 3 years (1 095 cycles) of operation at design conditions.  This point lies slightly above the “average” life curve as seen in Figure 4.

    After 12 months of operation, a Performance Test Run (PTR) is conducted.  The results are shown in Table 3.  The feed flow rate and temperature are slightly lower compared to the design values.  A water breakthrough time of 20.9 hours is recorded.  The FL is determined (using the concepts in Chapter 18) to be 0.68 after 365 cycles (one year of operation).  It is important and useful to understand the equation sequence of the concepts in Chapter 18, as shown by Equations 18.5 through 18.10 to arrive at the cited value for FL.   This data point is shown in Figure 5 and is seen to lie just below the generic “Average” curve.   Note that the slope of the curves are starting to flatten out.  Since the PTR FL is lower than the Design FL, the molecular sieves will experience water breakthrough if operated at design conditions in less than three years. Figure 6 shows the projected life factor, FL, after 3 years of service at design conditions.  If the capacity decline continues to follow the same trend as seen from the PTR, water breakthrough will occur after 750 cycles or just a little over 2 years from startup if operation continues at design conditions.  This is shown in Figure 7.

    Figure 4. Design condition life factor [1]
    Figure 4. Design condition life factor [1]

    Table 3. Results of Performance Test Run (PTR) after 12 months of operation

    tab3

    Figure 5. Performance test run (PTR) life factor [1]
    Figure 5. Performance test run (PTR) life factor [1]
    Figure 6. Projected life factor (red triangle) running at design conditions [1]
    Figure 6. Projected life factor (red triangle) running at design conditions [1]
    Figure 7. Projected life factor running at design conditions [1]
    Figure 7. Projected life factor running at design conditions [1]
                Because the unit has a regeneration circuit that can handle an additional 15% of flow, the complete regeneration cycle (heating, cooling, de- and re- pressurization) can be reduced to 7.0 hours.  This allows the beds to turn around faster.

    Using the reduced cycle time (the complete cycle time is now 21 hours vs the original 24 hours), we find an FL = 0.53.  This is because less water is being adsorbed per cycle.  This occurs at around the 1500 cycle mark as shown in Figure 8.

    If the plant elects to take advantage of the standby time and operate at reduced cycle time immediately following the PTR, the molecular sieves should last an additional 2.7 years, resulting in a total life of 3.7 years.   In this case, standby time will allow the unit to operate until the scheduled plant turnaround.

    Figure 8. Projected life factor (red triangle) if standby time is used [1]
    Figure 8. Projected life factor (red triangle) if standby time is used [1]
                We can draw the following conclusions from this case study:

    1. The methods presented allow the user to estimate the decline of their adsorbent based on only one performance test run for molecular sieve dehydrators using low pressure regeneration. This permits early formulation of a credible action plan.
    2. Site-specific factors will determine your unit’s decline curve. Consequently, conducting more than one performance test is highly recommended.  A poorly performing inlet separator, for example, could result in a unit exhibiting a more pronounced decline than indicated by the generic curves in Figure 3.
    3. Standby time offers a large degree of operating flexibility because the decline curves tend to level off; always try to build in standby time in any new molecular sieve design.
    4. Adsorption capacity is a function of the number of cycles, not calendar time.
    5. Install a good filter coalescer or filter separator upstream of your adsorption unit to keep the contaminants out of the system.

    The approach discussed in this Tip of the Month should help a facility engineer plan for the inevitable replacement of molecular sieves in their natural gas dehydration facility.

    To learn more about similar cases and how to minimize operational problems, we suggest attending our G4 (Gas Conditioning and Processing) and PF4 (Oil Production and Processing Facilities) courses.

    PetroSkills offers consulting expertise on this subject and many others. For more information about these services, visit our website at http://petroskills.com/consulting, or email us at consulting@PetroSkills.com.

    By: Harvey M. Malino

    Reference:

    1. Campbell, J.M., Gas Conditioning and Processing, Volume 2: The Equipment Modules, 9th Edition, 2nd Printing, Editors Hubbard, R. and Snow–McGregor, K., Campbell Petroleum Series, Norman, Oklahoma, 2014.
  • Process Safety and Low Oil Prices

    In this Tip of the Month, we reflect back to December 2008, and get a reminder from the United States Chemical Safety Board (CSB) to remain focused on process safety and accident prevention during this time of lower oil prices.

    During the economic downturn of 2008, oil prices dropped significantly. The latest drop in crude oil prices is similar. At that time, the CSB produced a video message asking companies to stay focused on process safety. That message is very relevant today.

    Process Safety and Low Oil Prices

    In the past, market conditions have occurred where oil prices have been low, such as we are experiencing today. Corporate cost cutting during these low oil price events have contributed to process safety incidents years later. In 2008, the United States Chemical Safety Board (CSB) Chairman John Bresland provided a reminder to oil companies that it is important to stay focused on process safety, even when prices are low. This was accomplished through a press release and a video safety message that is appropriate for this time [1].

    Low oil prices, combined with striking workers at US refineries increase the challenges faced by managers to insure that process safety is a core value of the organization.

    Containing overhead and operating costs during these market conditions may lead some to take shortcuts and make hasty decisions without considering all the process safety implications of these decisions. The attached press release and video safety message is as appropriate today as it was in 2008. This video message would be an excellent safety moment topic and hopefully will allow us to remain focused on process safety.

    Dec 22, 2008

    In First Video Safety Message, CSB Chairman John Bresland Calls for Industry to Remain Focused on Process Safety, Accident Prevention During Recession

    Washington, DC, December 22, 2008 – In his first video safety message, CSB Chairman John Bresland today said that chemical companies and refineries need to continue to invest in process safety and preventive maintenance, even as the economic downturn cuts into sales and profits.

    The four-minute video message was released on YouTube.com (http://www.youtube.com/safetymessages) and the text was posted on Blogger.com (http://safetymessages.blogspot.com).

    “My safety message for oil and chemical companies is clear: even during economic downturns, spending for needed process safety measures must be maintained,” Chairman Bresland stated in the message. He noted that the CSB investigation of the 2005 Texas City refinery disaster linked the accident to corporate spending decisions in the 1990s, when low oil prices triggered cutbacks in maintenance, training, and operator positions at the plant.

    “Unfortunately, around the country, refinery accidents continue to be a concern,” Chairman Bresland said, pointing to three major accidents that occurred at refineries in Texas this year, including a fire at a refinery in Tyler last month that fatally burned two workers and forced the refinery to shut down for months. “Today, as gasoline prices remain low, companies should weigh each decision to make sure that the safety of plant workers, contractors, and communities is protected.”

    Safety Messages are a new communication tool for the agency, consisting of short videos from the Chairman or the other board members. In the coming weeks and months, new messages will be released on a variety of current issues in chemical process safety.

    “I encourage all of our stakeholders to join the discussion on YouTube.com and Blogger.com and share their thoughts about the subject of these messages,” Chairman Bresland said. Comments and ideas for future Safety Messages can also be emailed to safetymessages@csb.gov.

    The CSB is an independent federal agency charged with investigating industrial chemical accidents. The agency’s board members are appointed by the president and confirmed by the Senate. CSB investigations look into all aspects of chemical accidents, including physical causes such as equipment failure as well as inadequacies in regulations, industry standards, and safety management systems.

    The Board does not issue citations or fines but does make safety recommendations to plants, industry organizations, labor groups, and regulatory agencies such as OSHA and EPA. Visit our website, www.csb.gov.

    For more information, contact Daniel Horowitz at (202) 261-7613 or Hillary Cohen at (202) 261-3601.

    To learn more about process safety, we suggest attending our PetroSkills HSE course, HS 45- Risk Based Process Safety Management or PS-2, Fundamentals of Process Safety To enhance process safety engineering skills we suggest any of the JMC foundation courses or our, PS 4 – Process Safety Engineering course.

    PetroSkills offers consulting expertise on this subject and many others. For more information about these services, visit our website at http://petroskills.com/consulting, or email us at consulting@PetroSkills.com.

    By: Clyde Young

    Reference:

    1. United States Chemical Safety Board, Press Release, December 22, 2008.
  • Troubleshooting Gas-Liquid Separators – Removal of Liquids from the Gas

    One of the most common problems in Oil and Gas Processing facilities is underperforming vapor / liquid separators. The most common types of gas-liquid separators are:

    • Slug catchers
      • Vessel / Finger-type
    • “Conventional” separators
      • Vertical / Horizontal
    • Scrubbers (i.e. Compressor Suction Scrubbers)
    • Gas “polishers”
      • Coalescing Filters / Filter Separators

    Underperforming separators generally result from either: 1. the wrong type of equipment was selected for the application, or 2. the correct type of equipment was selected, but the sizing methodology was inadequate. The type of separator required for an application depends largely on the gas-liquid ratio of the stream to be treated, and the flow variability of the process, as shown in Figure 1. As the flow variability increases with low to moderate gas-liquid ratios, the separator selection will move from a conventional separator to a slug catcher. For applications where there is a high gas-liquid ratio (i.e. very low liquid content), and the flow variability is moderately low, scrubbers and gas polishers would be the appropriate equipment selection depending upon the gas quality requirement for the treated stream.

    Figure 1. Gas-Liquid Separation Equipment Selection Map [1]
    Figure 1. Gas-Liquid Separation Equipment Selection Map [1]
    Figure 2. Troubleshooting Methodology [1]
    Figure 2. Troubleshooting Methodology [1]
    Unfortunately, once the equipment has been selected and installed, it is very costly to replace if the separator was not specified properly. Common separator performance issues are: too much liquid in the separated gas, inadequate slug/surge capacity, and too much gas in the separated liquid.

    This paper is focused on troubleshooting inadequate liquid removal from the gas for conventional separators (moderate to high liquid loads) and scrubbers (very low liquid loads).

    A troubleshooting methodology is provided in Figure 2 [Reference No]. The problem, in this case, is too much liquid in the separated gas stream. In order to effectively troubleshoot separator performance, it is required to understand the metrics of good performance, and the functions and analysis of the various components of the separation equipment.

    Typical performance metrics for separators are provided in Table 1. The specific performance requirements for a given separator are set by the sensitivity of the downstream process / equipment to the presence of liquids. For example, the product gas (sales gas) off of the cold separator in an NGL Extraction facility is sensitive to the presence of entrained liquids. The product gas can go off specification if there is too much carryover of liquids from the cold separator. On the other hand, the sensitivity of the downstream equipment from the facility inlet separator is much less, and the amount of liquids entrained in the gas is more tolerable.Table 1. Example separator performance metrics [1]

    tab1

    Separator Components

    The main components of a separator, shown in Figure 3, are the feed pipe, inlet device, gas gravity separation section, mist extractor and the liquid gravity separation section. The gas/liquid separator components will be briefly discussed in regards to their effects on gas/liquid separation performance. These effects need to be understood and quantified in order to troubleshoot separator operations, and to identify modifications that can be made to improve performance. The liquid gravity separation section will not be discussed.

    Figure 3. Parts of a Conventional Separator [2]
    Figure 3. Parts of a Conventional Separator [2]
    Inlet Feed Pipe

    The inlet feed pipe sizing and geometry is important as it is desired to keep the multiphase flow pattern “stabilized” in the piping to minimize the production of small liquid droplets, and liquid entrainment into the gas phase. Figure 4 [2] shows the effect of feed pipe velocity on liquid entrainment. Figure 5 [2] demonstrates how quickly the liquid entrainment increases once the entrainment inception point is reached.

    Figure 4. Effect of feed pipe velocity on liquid entrainment [2]
    Figure 4. Effect of feed pipe velocity on liquid entrainment [2]
    Figure 5. Example of liquid entrainment behavior in a gas-liquid system [2]
    Figure 5. Example of liquid entrainment behavior in a gas-liquid system [2]
    Some general guidelines for inlet piping to minimize liquid entrainment are:

    • Provide 10 diameters of straight pipe upstream of the inlet nozzle without valves, expansions/contractions or elbows.
    • If a valve is required, only use full port gate or ball valves.

    Inlet Device

    The main purpose of an inlet device is to improve separation performance. This is achieved by maximizing the amount of gas-liquid separation occurring in the feed pipe, minimizing droplet shearing, and optimizing the downstream velocity distributions of the separated phases into the separator. Schematics for inlet devices are shown in Figure 6. In large capacity, more critical separation applications, the vane-type and cyclonic inlet devices are commonly used. The simpler, and less expensive, impact (or diverter plates) are often used where the separation performance is less critical.

    Figure 6. Various separation equipment inlet devices [2]
    Figure 6. Various separation equipment inlet devices [2]
    Table 2 provides a comparison of the performance of various inlet devices.

    Table 2. Comparison of inlet devices [2]

    tab2

    The inlet momentum (ρmV2m – density*velocity2 of the mixture) of the feed stream is typically used to select and size inlet devices. Table 3 provides the suggested upper limits of inlet momentum values. For conditions where it is not practical to avoid higher feed pipe velocities, it must be recognized that failure to do so will result in higher entrainment loads, smaller droplet sizes, and more difficult separation conditions.

    Table 3. Inlet device ρV2 upper limits [3]

    tab3

    The quality of the flow distribution downstream of the inlet device is critical to mist extractor performance. Historically, tracer surveys have been used to provide an approximate indication of the continuous phase velocity within separators. In more recent years, the use of CFD (Computational Fluid Dynamics) has provided insight into the flow behavior of fluids, and has resulted in significant improvement in separator internals design. Separator performance is to a large degree dependent on the removal of droplets/ bubbles from the continuous phase. The efficiency of this removal is a function of the continuous phase velocity, thus the importance of understanding the factors that affect velocity profiles. Figure 7 provides an example of ideal versus actual gas velocity profiles within a separator.

    fig7
    Figure 7. Ideal and actual gas velocity profiles [3]
    Gas Gravity Separation Section

    The gas gravity separation section of a separator has two main functions: 1) reduction of entrained liquid load not removed by the inlet device, 2) improvement / straightening of the gas velocity profile.

    Most mist extractors have limitations on the amount of entrained liquid droplets that can be efficiently removed from the gas, thus the importance of the gas gravity section to remove the liquids to an acceptable level upstream of the mist extractor. This is particularly important for separators handling higher liquid loads. For scrubber applications with low liquid loadings, the Ks value will be primarily dependent on the mist extractor type, and the gas gravity separation section becomes less important.

    For the higher liquid load applications, there are two approaches for sizing the gravity separation section to remove liquid droplets from the gas: 1) Ks method, 2) Droplet settling theory.

    Historically the Ks method has been employed as it can provide reasonable results and is easy to use, but has shortcomings in terms of quantifying separator performance. References 3-5 provide the details on the droplet settling theory methods which can be used to more accurately quantify separator performance. The Ks approach is limited in that it only informs of the average droplet size, but cannot quantify the amount of liquid droplets exiting the gas gravity section.

    The Ks method (Equation 1) is an empirical approach to estimate the maximum allowable gas velocity to achieve a desired droplet separation.eq1

    Where:

    ρL             = liquid density kg/m3 (lbm/ft3)

    ρg         = gas density kg/m3 (lbm/ft3)

    Vgmax = maximum allowable gas velocity m/s (ft/sec)

    KS        = an empirical constant m/s (ft/sec)

    Figure 8 provides the relationship of Ks values for various droplet sizes and separator operating pressures for the gas gravity section. Typically, a Ks value is selected that will achieve removal of all entrained droplets larger than a specified target droplet diameter in the original design of the separator. For conventional separators the target droplet diameter is typically 150 microns, and for scrubbers the target droplet size should not exceed ~500 microns. This correlation can also be used to determine the performance of the gas gravity section based upon current operating conditions. The separator Ks value can be estimated from the actual velocity and fluid conditions, and the droplet size removed in the gas gravity section can be estimated from Figure 8.

    Figure 8. Ks vs. pressure and droplet size for empty vessels [2]
    Figure 8. Ks vs. pressure and droplet size for empty vessels [2]
    Mist Extraction Section

    The mist extractor is the final gas cleaning device in a conventional separator. The selection, and design to a large degree, determine the amount of liquid carryover remaining in the gas phase. The most common types include wire mesh pads (“mesh pads”), vane-type (vane “packs”) and axial flow demisting cyclones. Figure 9 shows the location and function of a typical mist extractor in a vertical separator.

    Figure 9. Typical mist extractor in a vertical separator [2]
    Figure 9. Typical mist extractor in a vertical separator [2]
    Mist extractor capacity is defined by the gas velocity at which re-entrainment of the liquid collected in the device becomes appreciable. This is typically characterized by a Ks value, as shown in Equation 1. Mesh pads are the most common type of mist extractors used in vertical separator applications. The primary separation mechanism is liquid impingement onto the wires, followed by coalescence into droplets large enough to disengage from the mesh pad. Figure 10 provides some mesh pad examples. Table 4 provides a summary of mesh pad characteristics and performance parameters.

    Figure 10. Mesh pad examples [1]
    Figure 10. Mesh pad examples [1]

    Table 4. Mesh pads summary of characteristics and performance parameters [1,4]

    tab4

    Notes:

    • Flow direction is vertical (upflow).
    • Assume mesh pad Ks values decline with pressure as shown in Table 5. Table 5 was originally developed for mesh pads, but is used as an approximation for other mist extractor types. [6].
    • If liquid loads reaching the mesh pad exceed the values given in Table 4, assume capacity (Ks) decreases by 10% per 42 L/min/m2 (10% per gal/min/ft2). [3-5].
    • These parameters are approximate.

    Table 5. Mesh pad Ks deration factors as a function of pressure [2]

    tab5

    Vane packs, like mesh pads, capture droplets primarily by inertial impaction. The vane bend angles force the gas to change direction while the higher density liquid droplets tend to travel in a straight-line path, and impact the surface of the vane where they are collected and removed from the gas flow. Figure 11 shows a schematic of a single-pocket vane mist extractor. Table 6 provides vane pack performance characteristics.

    Figure 11. Single-pocket vane schematic [2]
    Figure 11. Single-pocket vane schematic [2]

    Table 6. Typical vane-pack characteristics [1,4]

    tab6

    Notes:

    1. Assume vane-pack Ks values decline with pressure as shown in Table 5.
    2. If liquid loads reaching the vane pack exceed the values given in Table 4, assume capacity Ks decreases by 10% per 42 L/min/m2 (10% per gal/min/ft2). [3-5].
    3. These parameters are approximate only. The vane-pack manufacturer should be contacted for specific information.

    In the case of demisting cyclones, the vendor should be consulted in regards to performance for the current operations of interest.

    Troubleshooting

    When troubleshooting a separator, one needs to quantify the acceptable performance of the separator in terms of the amount of liquids in the separated gas. The separator physical condition and design is then assessed to determine the liquid removal capability of the separation equipment installed. Each separator component should be analyzed in terms of the current operating conditions versus the original design specifications.

    Table 7 provides a few common causes of inadequate liquid removal performance of a separator. The separator components that need to be reviewed are identified to determine the specific limitation. This table can serve as a road map for the calculations to work through when doing this type of analysis.

    Table 7. Common conditions resulting in inadequate separated gas quality [1]

    tab7

    There are numerous options available to improve the performance of a separator depending upon what the cause for the poor performance is. Depending upon the size and construction of the separator, it may be possible to retro-fit the separator internals. Another option may be modification of the inlet feed piping geometry and number of fittings upstream of the vessel if this is found to be less than ideal. The inlet device may be damaged, or in the bottom of the vessel. Higher efficiency inlet devices may be an option for consideration. Frequently, different mist extraction equipment can be selected to improve performance. For example, if the mist extractor Ks value is greater than the original design, a different mist extraction device could improve performance. The separator internals modifications may or may not be possible without welding on the vessel (which adds additional complications and cost to the project).

    The operating liquid levels should also be reviewed in terms of the distance of the normal operating liquid level in relation to the inlet feed device. If the liquid level is too high, the gas velocity from the inlet could be re-entraining liquids that were previously separated in the feed piping / inlet device. Unfortunately, in some cases the only way to improve performance is to cut rate (i.e. gas velocity) through the separator.

    To learn more about troubleshooting separators and other production equipment, we suggest attending our PF-49 (Troubleshooting Oil and Gas Processing Facilities), or PF-42 (Separation Equipment Selection and Sizing) for more details on the selection and specification of separators.

    By: Kindra Snow-McGregor

    Reference:

    1. PF-49, Troubleshooting Oil and Gas Processing Facilities, Bothamley, M., 2014, © PetroSkills, LLC. All Rights reserved.
    2. Campbell, J.M., Gas Conditioning and Processing, Volume 2: The Equipment Modules, 9th Edition, 2nd Printing, Editors Hubbard, R. and Snow–McGregor, K., Campbell Petroleum Series, Norman, Oklahoma, 2014.
    3. Bothamley, M. 2013. Gas-Liquid Separators – Quantifying Separation Performance Part 1. SPE Oil and Gas Facilities, Aug. (22 – 29).
    4. Bothamley, M. 2013. Gas-Liquid Separators – Quantifying Separation Performance Part 2. SPE Oil and Gas Facilities, Oct. (35 – 47)
    5. Bothamley, M. 2013. Gas-Liquid Separators – Quantifying Separation Performance Part 2. SPE Oil and Gas Facilities, Dec. (34 – 47)
    6. Fabian, P., Cusack, R., Hennessey, P., Neuman, M. 1993. Demystifying the Selection of Mist Eliminators, Part 1: The Basics. Chem Eng 11 (11): 148 – 156.

  • Impact of Gas-Oil Ratio (GOR) on Crude Oil Pressure Drop in Gathering Systems

    The use of multiphase flow systems is common practice in the oil and gas industry. Multiphase flow is often encountered in the well tubing, flow lines and gathering systems. For transport of oil and gas (and water) to downstream processing facilities the preference is normally a single pipeline in which both phases are transported simultaneously for economic reasons. Even in gas pipelines where the gas enters the line as a single phase fluid, condensation of liquids can occur due to pressure and temperature changes along the line.

    Modeling and simulation of a multiphase systems, even under steady-state conditions, is complex. There are a few tools designed specifically for modeling and analysis of complex multiphase systems such as PipePhase, PipeSim, OLGA, etc. [1].

    In the June 2008 Tip of the Month (TOTM), we demonstrated how general-purpose process simulation programs can be used to simulate gas dominated two-phase pipelines. In the August 2008 TOTM, we discussed the value of the simple Flanigan correlation and how it can be used to model and analyze the behavior of a wet gas transmission pipeline. The results of the Flanigan correlation were compared with more rigorous calculation methods for multiphase pipelines.

    In this TOTM, we will study the impact of gas-oil ratio (GOR) on pressure drop in crude oil gathering systems. Specifically, pressure drop along a gathering line for nominal pressures of 690, 3450, and 6900 kPag (100, 500, and 1000 psig) and nominal pipe size of 101.6 and 152.4 mm (4 and 6 inches) was calculated using multiphase rigorous method from commercial simulator. The calculated pressure drops are presented in graphical format as a function of the oil stock tank volume flow rate and GOR. Variation of thermo physical properties was considered.

    Case Study

    For the purpose of illustration, we considered a case study for transporting a crude oil of relative density of 0.852 (°API = 34.6) at stock tank condition combined with a gas with relative density of 0.751. The selected GORs were 0 (dead oil), 17.8, 356.5, and 891.3 Sm3 of gas/STm3 of oil (0, 100, 2000, and 5000 scf/STB). The compositions of oil and gas are presented in Table 1. The oil C6+ was characterized as 30 single carbon number (SCN) [2] ranging from SCN6 to SCN35 while the gas C6+ was characterized by 10 SCN ranging from SCN6 to SCN15. For details of the SCN components, see Table 3.2 on page 64 of reference [2]. The mole fraction of SCN components were determined by an exponential decay algorithm [3].

    Table 1. Feed composition at stock condition

    table1

    The following assumptions were made:

    1. Steady state conditions
    2. The line is 1.601 km (1 mile) long with nominal size of 101.6 and 152.4mm (4 and 6 inches), onshore buried line.
    3. Segment lengths and elevation changes are presented in Table 2. This elevation profile is considered to be approximately equivalent to “rolling” terrain.
    4. Pipeline inside surface roughness of 46 microns (0.046 mm, 0.0018 inch)
    5. Line nominal pressure 690, 3450, and 6900 kPag (100, 500, and 1000 psig)
    6. The feed enters the line at 15.6 ˚C and (60 ˚F)
    7. The ground/ambient temperature, is 15.6 ˚C and (60 ˚F)
    8. Water cut is 0 (no water in the feed).
    9. Overall heat transfer coefficients of 2.839 W/m2-˚C (0.5 Btu/hr-ft2-˚F), for onshore buried line (minor effect as inlet temperature = ambient ground temperature).
    10. Simulation software ProMax [4] and using the Soave-Redlich-Kwong (SRK) Equation of State [5] for vapor-liquid equilibrium and Beggs-Brill method for two-phase pressure drop calculation [6].

    Table 2. Line segment length and elevation change

    table2

    Results and Discussions:

    The two phase (oil and gas) flow through the gathering line was simulated by ProMax with SRK EOS for vapor-liquid equilibria and Beggs-Brill for two phase pressure drop calculations. Figures 1A and 1B present the calculated pressure drop per unit length as a function of oil stock tank volume rate and GOR for nominal line diameter of 101.6 mm (4 inches) at nominal line pressure of 690 kPag (100 psig) in SI (international) and FPS (Engineering) system of units, respectively. Figures 1A and 1B indicate that as the GOR increases from 0 to 891 Sm3/STm3 (0 to 5000 scf/STB), the pressure drop increases considerably. Consequently, as the GOR increases, the line capacity decreases.

    Figures 2A, 2B, 3A, and 3B present the results for the same line size but at nominal pressures of 3445 and 6900 kPag (500 and 1000 psig), respectively. Contrary to Figure 1, Figures 2 and 3 indicate that at these higher pressures as the GOR increases, the pressure drop decreases for low GOR value. However, for further increase of GOR the pressure drop increases considerably.

    Similar calculations were repeated for another line with nominal pipe size of 152.4 mm (6 inches) and the simulation results are presented in Figures 4 through 6. Figures 4 through 6 also demonstrate the same impact of GOR on the pressure drop, at higher pressures and low GOR, the pressure drop decreases. However, the impact of low GOR at higher pressures is less compared to the smaller line diameter.

    fig1a

    Figure 1A (SI). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 690 kPag for 101.6 mm pipe diameter

    fig1b

    Figure 1B (FPS). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 100 psig for 4 in pipe diameter

    fig2a

    Figure 2A (SI). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 3445 kPag for 101.6 mm pipe diameter

    fig2b

    Figure 2B (FPS). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 500 psig for 4 in pipe diameter

    fig3a

    Figure 3A (SI). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 6900 kPag for 101.6 mm pipe diameter

    fig3b

    Figure 3B (FPS). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 1000 psig for 4 in pipe diameter

    Conclusions

    The following conclusions can be made based on this case study:

    1. The GOR has a large impact on the capacity of crude oil gathering lines. In general as GOR increases the pressure drop increases which lowers the line capacity.
    2. At high pressures and low GOR, pressure drop is lower than the pressure drop for dead oil (solution gas is zero) because the viscosity of live oil is lower than viscosity of dead oil. This effect is bigger for the smaller line diameter.

    To learn more about similar cases and how to minimize operational problems, we suggest attending our PF 45 (Onshore Gas Gathering Systems: Design and Operation), G4 (Gas Conditioning and Processing), PF81 (CO2 Surface Facilities), PF4 (Oil Production and Processing Facilities), and PL4 (Fundamentals of Onshore and Offshore Pipeline Systems) courses.

    John M. Campbell Consulting (JMCC) offers consulting expertise on this subject and many others. For more information about the services JMCC provides, visit our website at www.jmcampbellconsulting.com, or email us at consulting@jmcampbell.com.

    By: Mahmood Moshfeghian

    Reference:

    1. Ellul, I. R., Saether, G. and Shippen, M. E., “The Modeling of Multiphase Systems under Steady-State and Transient Conditions – A Tutorial,” The Proceeding of Pipeline Simulation Interest Group, Paper PSIG 0403, Palm Spring, California, 2004.
    2. Campbell, J.M., Gas Conditioning and Processing, Volume 1: The Basic Principles, 9th Edition, 2nd Printing, Editors Hubbard, R. and Snow–McGregor, K., Campbell Petroleum Series, Norman, Oklahoma, 2014.
    3. Moshfeghian, M., Maddox, R.N., and A.H. Johannes, “Application of Exponential Decay Distribution of C6+ Cut for Lean Natural Gas Phase Envelope,” J. of Chem. Engr. Japan, Vol 39, No 4, pp.375-382 (2006)
    4. ProMax 3.2, Bryan Research and Engineering, Inc., Bryan, Texas, 2014.
    5. Soave, G., Eng. Sci. Vol. 27, No. 6, p. 1197, 1972.

    Brill, J. P., et al., “Analysis of Two-Phase Tests in Large-Diameter Flow Lines in Prudhoe Bay Field,” SPE Jour, p. 363-78, June 1981.

    fig4a

    Figure 4A (SI). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 690 kPag for 152.4 mm pipe diameter

    fig4b

    Figure 4B (FPS). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 100 psig for 6 in pipe diameter

    fig5a

    Figure 5A (SI). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 3445 kPag for 152.4 mm pipe diameter

    fig5b

    Figure 5B (FPS). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 500 psig for 6 in pipe diameter

    fig6a

    Figure 6A (SI). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 6900 kPag for 152.4 mm pipe diameter

    fig6b

    Figure 6B (FPS). Variation of pressure drop per unit length with oil stock tank volume rate and GOR at 1000 psig for 6 in pipe diameter

  • Lean Sweet Natural Gas Water Content Correlation

    In the October, November, December 2007 and February 2014 Tips of the Month (TOTM), we studied in detail the water phase behaviors of sweet and sour natural gases and acid gas systems. We also evaluated the accuracy of different methods for estimating the water content of sour natural gas and acid gas systems.

    The water vapor content of natural gases in equilibrium with water is commonly estimated from Figure 6.1 of Campbell book [1] or Figure 20.4 of Gas Processors and Suppliers Association, including corrections for the molecular weight (relative density) of gas and salinity of water [2].

    In this TOTM, we will present two new correlations for estimating the water content of lean and sweet natural gases. The performance of the proposed correlations will be compared with the rigorous simulation and shortcut method software and other correlations.

    Low Pressure System

    At low pressure conditions, less than 700 kPa (100 psia), the mole fraction of water in the gas phase can be estimated by dividing water vapor pressure, PV, at the specified temperature, T, by the system pressure, P. The vapor pressure of pure water, from 0 to 360, (32 to 680) can be calculated by the following relation [3].

    eq 1

    Where:

    eq 2

     

    The critical temperature, TC = 647.096 K and critical pressure, PC = 22.064 MPa, T in K, and PV in MPa, and

    a1 = −7.85951783,     a2 = 1.84408259,     a3 = −11.7866497,      a4 = 22.6807411,

    a5 = −15.9618719,     a6 = 1.80122502

    Knowing  one kmole of water = 18 kg (lbmole=18 lbm) and one kmole of gases occupy 23.64 Sm3 at standard condition of 15  and 101.3 kPa (one lbmole of gases occupy 379.5 SCF at standard condition of 60 and 14.7 psia), the water content is calculated by

    eq 3

    Moderate to High Pressure System

    For pressures higher than 700 kPa (100 psia), we propose a correlation similar to equation 6-276 in Chapter 6 of Standard Petroleum Handbook [4] as follows:

    eq 4

     

     

    Reference [4] presents the tabular values of A and B as a function of temperature. In this work, the temperature dependency of A and B in equation 3 is presented in the form of Gaussian Model.

    eq 5

    Where:

    table a

    A temperature range of -40 to 100 (-40 to 212) and pressure range of 6.8 to 680 atm have been considered. For the purpose of higher accuracy, the parameters in equations 4G are regressed and presented in Table 1 for 5 different temperature intervals, for SI (System International) and FPS (Foot-pound-Second) system of units. The temperature range in the 5th row in each system of units is used more frequently. The proposed correlations are suitable for spreadsheet calculations.

                Table 1. The parameters for equation 4G (Gaussian Model).

    table1

    Alternatively, the temperature dependency of parameters A and B in equation 3 can also be presented in the form of Polynomial Model.

    eq 6

     

     

    Similarly, the correlation parameters of equations 4P are presented in Table 2.

    The performance of the proposed correlation was evaluated against the water content of a lean sweet natural gas with a relative density (specific gravity) of 0.6 calculated by SRK equation of state of ProMax software [5]. The water saturator tool of ProMax was utilized. Figure 1A (SI) and Figure 1B (FPS) present the comparison between the proposed correlation (Eqs. 3 and 4G) results (solid curves) and the corresponding results from ProMax designated by geometric symbols.  Figure 1 indicates there is relatively good agreement between the proposed correlation and ProMax. Large deviations are observed for temperatures below -20 (-4). Almost in all conditions, except at very high temperature the proposed correlation is more conservative (over predicting) than the ProMax.

    Table 2. The parameters for equation 4P (Polynomial Model).

    table 2

    Evaluation of the Proposed Correlation

    The performance of the proposed correlation was also evaluated against the water content of a lean sweet natural gas predicted by GCAP  software [6].  The water content of GCAP is based on Figure 6.1 of Campbell book [1]. Figure 2A (SI) and Figure 2B (FPS) present the comparison between the proposed correlation (Eqs. 3 and 4G) results (solid curves) and the corresponding results from GCAP designated by geometric symbols.  Figure 2 indicates there is better agreement between the proposed correlation and GCAP compared to ProMax. Large deviations are observed for temperatures below -20 (-4). Almost in all conditions, except at very high temperature, the proposed correlation is slightly more conservative (over predicting) than the GCAP.

    Figures 1 and 2 indicate that in the temperature range of 0 to 70 (32 to 158 ) excellent agreement is observed between results of equation 3 and both software. Figure 2 also indicates that the GCAP results for 170 atm at very low temperature is inconsistent with other pressure results of GCAP.

    When equation 4P and the parameters presented in Table 2 were used instead of equation 4G, similar quality of results was obtained.

    fig1a

    Figure 1A (SI). Comparison of results between the proposed correlation (Eq 3) and ProMax

    fig1b

    Figure 1B (FPS). Comparison of results between the proposed correlation (Eq 3) and ProMax

    fig2a

    Figure 2A (SI). Comparison of results between the proposed correlation (Eq 3) and GCAP

    fig2b

    Figure 2B (FPS). Comparison of results between the proposed correlation (Eq 3) and GCAP

    Bukacek Correlation

    Bukacek [7] suggested a relatively simple correlation for the water content of lean sweet gas as follows:

    eq 7

    where T is in °F.

    This correlation is reported to be accurate for temperatures between  15.5 and 238°C (60 and 460°F) and for pressure from 0.105 to 69.97 MPa (15 to 10,000 psia). The pair of equations in this correlation is simple in appearance. The added complexity that is missing is that it requires an accurate estimate of the vapor pressure of pure water. In this study, we have used equation 1 for water vapor pressure.

    The performance of the proposed correlation was also evaluated against the Bukacek’s Correlation coupled with equation 1 for water vapor pressure. Figure 3A (SI) and Figure 3B (FPS) present the comparison between the proposed correlation (Eqs. 3 and 4G) results (solid curves) and the corresponding results from Bukacek’s correlation designated by geometric symbols.  Figure 3 indicates there is an excellent agreement between the proposed correlation and Bukacek’s correlation.

    Conclusions

    The following conclusions can be made regarding the proposed correlation:

    1. A relatively simple correlation for predicting the water content of lean sweet natural gas is presented that can be used for spreadsheet calculation.
    2. Based on Figures 1 and 2, the agreement between the predicted water content by the proposed correlation (Eqs. 3 and 4G or Eqs. 3 and 4P) and those predicted by ProMax and GCAP software is relatively good. The agreement deteriorates at temperatures below -20 (-4).
    3. In general, the estimated water content by the proposed correlation is conservative compared to ProMax and GCAP.
    4. A better agreement between the proposed correlation and GCAP compared to ProMax is observed.
    5. Excellent agreement is observed between the proposed and the Bukacek’s correlations.
    6. The GCAP results for pressure of 170 atm at low temperatures are inconsistent with respect to GCAP results at other pressure.
    7. The proposed correlation is easy and suitable for hand or spreadsheet calculations.

    To learn more about similar cases and how to minimize operational problems, we suggest attending our G6 (Gas Treating and Sulfur Recovery), G4 (Gas Conditioning and Processing), PF81 (CO2 Surface Facilities), PF4 (Oil Production and Processing Facilities), and PL4 (Fundamentals of Onshore and Offshore Pipeline Systems) courses.

    John M. Campbell Consulting (JMCC) offers consulting expertise on this subject and many others. For more information about the services JMCC provides, visit our website at www.jmcampbellconsulting.com, or email us at consulting@jmcampbell.com.

    By: Dr. Mahmood Moshfeghian

    Reference:

    1. Campbell, J.M., Gas Conditioning and Processing, Volume 1: The Basic Principles, 9th Edition, 2nd  Printing, Editors Hubbard, R. and Snow–McGregor, K., Campbell Petroleum Series, Norman, Oklahoma, 2014.
    2. GPSA Engineering Data Book, Section 20, Volume 2, 13th Edition, Gas Processors and Suppliers Association, Tulsa, Oklahoma, 2012.
    3. Wagner, W.  and Pruss, A.,  J. Phys. Chem. Reference Data, 22, 783–787, 1993.
    4. Standard Handbook of Petroleum, Natural Gas Engineering volume 2, Lyons, W. C., Editor, Gulf Professional Publishing, Houston, Texas, 1996
    5. ProMax 3.2, Bryan Research and Engineering, Inc., Bryan, Texas, 2014.
    6. GCAP 9.1, Gas Conditioning  and Processing, PetroSkills/Campbell, Norman, Oklahoma, 2014
    7. Bukacek, R.F., “Equilibrium Moisture Content of Natural Gases” Research Bulletin IGT, Chicago, vol 8, 198-200,  1959.

    fig3a

    Figure 3A (SI). Comparison of results between the proposed correlation (Eq 3) and Bukacek’s Correlation

    fig3afps

    Figure 3A (FPS). Comparison of results between the proposed correlation (Eq 3) and Bukacek’s Correlation