Difference between revisions of "Hagedorn and Brown correlation"

From wiki.pengtools.com
Jump to: navigation, search
(Nomenclature)
(Nomenclature)
Line 100: Line 100:
 
:<math> \mu_L </math> = liquid viscosity, cp
 
:<math> \mu_L </math> = liquid viscosity, cp
 
:<math> \mu_g </math> = gas viscosity, cp
 
:<math> \mu_g </math> = gas viscosity, cp
 +
:<math> \mu_o </math> = oil viscosity, cp
 +
:<math> SG_o </math> = oil specific gravity, dimensionless
 +
:<math> SG_w </math> = water specific gravity, dimensionless
 +
:<math> SG_g </math> = gas specific gravity, dimensionless
 +
:<math> WOR </math> = water-oil ratio, bbl/bbl
 +
:<math> GLR </math> = gas-liquid ratio, scf/bbl
 +
:<math> R_s </math> = solution gas-oil ratio, scf/stb
 +
:<math> z </math> = gas compressibility factor, dimensionless
 +
:<math> T_R </math> = temperature, °R
 +
:<math> \sigma_L </math> = surface tension of liquid-air interface, dynes/cm
 +
:<math> \sigma_o </math> = surface tension of liquid-oil interface, dynes/cm
 +
:<math> \sigma_w </math> = surface tension of liquid-water interface, dynes/cm
  
 
== References ==
 
== References ==

Revision as of 12:46, 24 March 2017

Brief

Hagedorn and Brown is an empirical two-phase flow correlation published in 1965 [1].

It doesn't distinguish between the flow regimes.

The heart of the Hagedorn and Brown method is a correlation for the liquid holdup H_L[2].

Math & Physics

Following the law of conservation of energy the basic steady state flow equation is:

 144 \frac{\Delta p}{\Delta h} = \bar \rho_m + \frac{f q_L^2 M^2}{2.9652 \times 10^{11} D^5 \bar \rho_m} + \bar \rho_m \frac{\Delta{(\frac{v_m^2}{2g_c}})}{\Delta h}[1]

where

 \bar \rho_m = \rho_L H_L + \rho_g (1 - H_L)[1]

Colebrook–White [3] equation for the Darcy's friction factor:

 \frac{1}{\sqrt{f}}= -2 \log \left( \frac { \varepsilon} {3.7 D} + \frac {2.51} {\mathrm{Re} \sqrt{f}} \right)[4]

Reynolds two phase number:

 Re = 2.2 \times 10^{-2} \frac {q_L M}{D \mu_L^{H_L} \mu_g^{(1-H_L)}}[1]

Discussion

Why Hagedorn and Brown?

The of the consistently best correlations was found to be the empirical Hagedorn and Brown correlation
— Economides
.
Cry "Havoc" and let slip the dogs of war.
— William Shakespeare, Julius Caesar, act III, scene I


One of the consistently best correlations was found to be the empirical Hagedorn and Brown correlation. [2]

Flow Diagram

HB Block Diagram

Workflow

To find H_L calculate:

 M =SG_o\ 350.52\ \frac{1}{1+WOR}+SG_w\ 350.52\ \frac{WOR}{1+WOR}+SG_g\ 0.0764\ GLR[1]
 \rho_L= \frac{62.4\ SG_o + \frac{Rs\ 0.0764\ SG_g}{5.614}}{B_o} \frac{1}{1+WOR} + 62.4\ SG_w\ \frac{WOR}{1 + WOR}[5]
 \rho_g = \frac{28.967\ SG_g\ p}{z\ 10.732\ T_R} [5]
 \mu_L = \mu_o \frac{1}{1 + WOR} + \mu_w \frac{WOR}{1 + WOR}[5]
 \sigma_L = \sigma_o \frac{1}{1 + WOR} + \sigma_w \frac{WOR}{1 + WOR}[5]
 N_L = 0.15726\ \mu_L \sqrt[4]{\frac{1}{\rho_L \sigma_L^3}}[1]
 CN_L = 0.061\ N_L^3 - 0.0929\ N_L^2 + 0.0505\ N_L + 0.0019 [2]
 v_{SL} = \frac{5.615 q_L}{86400 A_p} \left ( B_o \frac{1}{1+WOR} + B_w \frac{WOR}{1+WOR} \right )[5]
 v_{SG} = \frac{q_L \left ( GLR-R_s \left( \frac{1}{1+WOR}\right) \right )}{86400 A_p}\ \frac{14.7}{p}\ \frac{T_K}{520}\ \frac{z}{1}[5]
 N_{LV} = 1.938\ v_{SL}\ \sqrt[4]{\frac{\rho_L}{\sigma_L}} [1]
 N_{GV} = 1.938\ v_{SG}\ \sqrt[4]{\frac{\rho_L}{\sigma_L}} [1]
 N_{D} = 120.872\ D \sqrt{\frac{\rho_L}{\sigma_L}} [1]
 H = \frac{N_{LV}}{N_{GV}^{0.575}}\  \left ( \frac{p}{14.7} \right )^{0.1} \frac{CN_L}{N_D} [2]
 \frac{H_L}{\psi} = \sqrt{\frac{0.0047+1123.32 H + 729489.64H^2}{1+1097.1566 H + 722153.97 H^2}} [6]
 B = \frac{N_{GV} N_{LV}^{0.38}}{N_{D}^{2.14}} [2]
 \psi = \begin{cases} 
27170 B^3 - 317.52 B^2 + 0.5472 B + 0.9999,  &\mbox{if } B <= 0.025 \\
-533.33 B^2 + 58.524 B + 0.1171, & \mbox{if }B > 0.025 \\
2.5714 B +1.5962, & \mbox{if }B > 0.055
\end{cases} [6]
 H_L = \frac{H_L}{\psi} \times \psi[1]

Nomenclature

 p = pressure, psia
 h = depth, ft
 H_L = liquid holdup factor, dimensionless
 \bar \rho_m = average mixture destiny at flowing conditions, lbm/ft2
 f = friction factor, dimensionless
 q_L = total liquid production rate, bbl/d
 M = total mass of oil, water and gas associated with 1 bbl of liquid flowing into and out of the flow string, lbm/bbl
 D = pipe diameter, ft
 V_m = mixture velocity, ft/sec
 q_c = conversion constant equal to 32.174, lbmft / lbfsec2
 \rho_L = liquid destiny, lbm/ft2
 \rho_g = gas destiny, lbm/ft2
 \varepsilon = absolute roughness, ft
 Re = Reynolds number, dimensionless
 \mu_L = liquid viscosity, cp
 \mu_g = gas viscosity, cp
 \mu_o = oil viscosity, cp
 SG_o = oil specific gravity, dimensionless
 SG_w = water specific gravity, dimensionless
 SG_g = gas specific gravity, dimensionless
 WOR = water-oil ratio, bbl/bbl
 GLR = gas-liquid ratio, scf/bbl
 R_s = solution gas-oil ratio, scf/stb
 z = gas compressibility factor, dimensionless
 T_R = temperature, °R
 \sigma_L = surface tension of liquid-air interface, dynes/cm
 \sigma_o = surface tension of liquid-oil interface, dynes/cm
 \sigma_w = surface tension of liquid-water interface, dynes/cm

References

  1. Jump up to: 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 Hagedorn, A. R.; Brown, K. E. (1965). "Experimental study of pressure gradients occurring during continuous two-phase flow in small-diameter vertical conduits". Journal of Petroleum Technology. 17(04): 475–484. 
  2. Jump up to: 2.0 2.1 2.2 2.3 2.4 Economides, M.J.; Hill, A.D.; Economides, C.E.; Zhu, D. (2013). Petroleum Production Systems (2 ed.). Westford, Massachusetts: Prentice Hall. ISBN 978-0-13-703158-0. 
  3. Jump up Colebrook, C. F. (1938–1939). "Turbulent Flow in Pipes, With Particular Reference to the Transition Region Between the Smooth and Rough Pipe Laws"Paid subscription required. Journal of the Institution of Civil Engineers. London, England. 11: 133–156. 
  4. Jump up Moody, L. F. (1944). "Friction factors for pipe flow"Paid subscription required. Transactions of the ASME. 66 (8): 671–684. 
  5. Jump up to: 5.0 5.1 5.2 5.3 5.4 5.5 Lyons, W.C. (1996). Standard handbook of petroleum and natural gas engineering. 2. Houston, TX: Gulf Professional Publishing. ISBN 0-88415-643-5. 
  6. Jump up to: 6.0 6.1 Trina, S. (2010). An integrated horizontal and vertical flow simulation with application to wax precipitation (Master of Engineering Thesis). Canada: Memorial University of Newfoundland.