Difference between revisions of "Liquid loading"
(→Nomenclature) |
(→See also) |
||
(101 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 1: | ||
__TOC__ | __TOC__ | ||
− | + | == Brief == | |
− | [[Liquid loading]] is a | + | [[Liquid loading]] is a phenomenon when the gas phase does't provide sufficient transport energy to lift the liquids out of the well. |
− | === Math & Physics | + | In '''1969''' Turner et al. published an empirical correlation defining the [[Liquid loading]] gas velocity<ref name=Turner/>. |
+ | |||
+ | [[File: Liquid Loading.png|x400px|Liquid Loading|right]] | ||
+ | |||
+ | == Math & Physics == | ||
The minimum gas velocity to remove the liquid equation: | The minimum gas velocity to remove the liquid equation: | ||
− | :<math> v_g = 1.593 \sigma^{1/4} \frac{({\rho_L-\rho_g})^{1/4}}{\rho_g^{1/2}}</math> | + | :<math> v_g = 1.593\ \sigma^{1/4}\ \frac{({\rho_L-\rho_g})^{1/4}}{\rho_g^{1/2}}</math><ref name= Lea/> (1) |
+ | |||
+ | ''Note that original paper<ref name=Turner/> uses 20.4 constant which is good if the surface tension units are in 'lbf/ft'. Also 20.4 is calculated by adjusting the theoretical equation constant by 20 percent. For the reference PROSPER uses 2.04. PQplot uses 1.593.'' | ||
The minimum gas rate to remove the liquid equation: | The minimum gas rate to remove the liquid equation: | ||
− | :<math> q_g = 3. | + | :<math> q_g = 3.067\ \frac{P\ v_g\ A}{T\ z} </math> (2) |
+ | |||
+ | == PQplot workflow == | ||
+ | In order to calculate the liquid loading point in [[PQplot]]:<br> | ||
+ | For every point in the [[VLP]] curve:<br> | ||
+ | 1. Critical gas velocity at the wellhead is calculated by eq (1)<br> | ||
+ | 2. Wellhead flowing gas velocity is calculated | ||
+ | |||
+ | :<math> v_{g wellhead} = \frac{q_g \times 10^6}{86400 A}\ \frac{14.7}{P}\ \frac{T}{520}\ \frac{z}{1}</math> | ||
+ | |||
+ | 3. Liquid loading flag is 1 if: | ||
+ | |||
+ | :<math> v_{g wellhead} < v_{g critical}</math> | ||
+ | |||
+ | 4. Critical gas rate is the last [[VLP]] curve point at which liquid loading flag is 1. | ||
+ | |||
+ | == Discussion == | ||
+ | To avoid the [[Liquid loading]] the gas velocity should be above the [[Liquid loading]] velocity. | ||
+ | |||
+ | The higher the gas rate the higher the gas velocity. | ||
− | === Nomenclature | + | The lower the wellhead flowing pressure the higher the gas rate. |
+ | |||
+ | The bigger the tubing ID the higher the gas rate. | ||
+ | |||
+ | In case when the gas rate is limited by the [[Reservoirs|Reservoir]] deliverability smaller tubing ID will increase the gas velocity. | ||
+ | |||
+ | ==Example. Calculating critical gas velocity and rate == | ||
+ | Example 3-2 from <ref name= Lea/> | ||
+ | ===Input data=== | ||
+ | :<math> P = 400 </math>, flowing wellhead pressure, psia | ||
+ | :<math> T =120 </math>, flowing wellhead temperature, °F | ||
+ | :<math> \rho_L =67</math>, liquid density, lbm/ft3 | ||
+ | :<math> \sigma = 60 </math>, water surface tension, dyne/cm | ||
+ | :<math> SG_g = 0.6 </math>, gas specific gravity | ||
+ | :<math> z =0.9 </math>, gas compressibility factor | ||
+ | Production string = 2-3/8 inch tubing with 1.995 in ID, A = .0217 ft2<br> | ||
+ | Water is the produced liquid | ||
+ | |||
+ | ===Solution=== | ||
+ | Calculate the gas density: | ||
+ | :<math> \rho_g = \frac{28.967\ SG_g\ p}{z\ 10.732\ T_R} = \frac{28.967\ 0.6\ 400}{0.9\ 10.732\ (120+460)}=1.24</math>, lbm/ft2 | ||
+ | |||
+ | |||
+ | Calculate the critical gas velocity form (1): | ||
+ | :<math> v_g = 1.593\ \sigma^{1/4}\ \frac{({\rho_L-\rho_g})^{1/4}}{\rho_g^{1/2}} = 1.593\ 60^{1/4}\ \frac{({67-1.24})^{1/4}}{1.24^{1/2}}=11.3</math>, ft/sec | ||
+ | |||
+ | |||
+ | Calculate the critical gas rate form (2): | ||
+ | :<math> q_g = 3.067\ \frac{P\ v_g\ A}{T\ z} = 3.067\ \frac{400\ 11.3\ 0.0217}{(120+460)\ 0.9}=0.578 </math>, MMscfd | ||
+ | |||
+ | ===Solution with PQplot=== | ||
+ | [[PQplot]] calculates the following critical gas rate: | ||
+ | :<math> q_g = 0.56 </math>, MMscfd | ||
+ | |||
+ | [[File:Liquid loading example 3-1.png|1000px| Liquid loading example 3-1]]<br> | ||
+ | [https://www.pengtools.com/pqPlot?paramsToken=c97c4f3b499cf59b82d61c374d3ab21a Liquid loading example 3-1 (Public model)] | ||
+ | |||
+ | == Nomenclature == | ||
:<math> A </math> = flow area, ft^2 | :<math> A </math> = flow area, ft^2 | ||
− | :<math> P </math> = flowing pressure, psia | + | :<math> P </math> = flowing wellhead pressure, psia |
:<math> q_g </math> = gas rate, MMscf/d | :<math> q_g </math> = gas rate, MMscf/d | ||
:<math> \rho_g </math> = gas density, lbm/ft3 | :<math> \rho_g </math> = gas density, lbm/ft3 | ||
:<math> \rho_L </math> = liquid density, lbm/ft3 | :<math> \rho_L </math> = liquid density, lbm/ft3 | ||
− | :<math> \sigma </math> = surface tension, dyne/cm | + | :<math> \sigma </math> = surface tension, dyne/cm (ref values: 60 - water, 20 - condensate) <ref name=Turner/> |
− | :<math> T </math> = flowing temperature, °R | + | :<math> SG_g </math> = gas specific gravity |
+ | :<math> T </math> = flowing wellhead temperature, °R | ||
:<math> v_g </math> = gas velocity, ft/sec | :<math> v_g </math> = gas velocity, ft/sec | ||
− | :<math> z </math> = gas compressibility factor at flowing P & T | + | :<math> z </math> = gas compressibility factor at flowing P & T, dimensionless |
+ | |||
+ | == See also == | ||
+ | [[Gray|Gray correlation]] | ||
+ | |||
+ | == References == | ||
+ | <references> | ||
+ | |||
+ | <ref name=Turner>{{cite journal | ||
+ | |last1= Turner |first1=R. G. | ||
+ | |last2= Hubbard |first2=M. G. | ||
+ | |last3= Dukler |first2=A. E. | ||
+ | |title=Analysis and Prediction of Minimum Flow Rate for the Continuous Removal of Liquids from Gas Wells | ||
+ | |journal=Journal of Petroleum Technology | ||
+ | |number=SPE-2198-PA | ||
+ | |date=Nov 1969 | ||
+ | |pages=1475–1482 | ||
+ | |url=https://www.scribd.com/doc/269398353/Friction-Factors-for-Pipe-Flow-MoodyLFpaper1944 | ||
+ | |url-access=registration | ||
+ | }}</ref> | ||
+ | |||
+ | <ref name=Lea>{{cite book | ||
+ | |last1= Lea |first1=J. F. | ||
+ | |last2= Nickens |first2=H. V. | ||
+ | |last3= Wells |first2=M. | ||
+ | |title=Gas Well Deliquification | ||
+ | |publisher=Energy Tribune Publishing Inc | ||
+ | |date=2011 | ||
+ | |url=https://www.sciencedirect.com/book/9780750682800/gas-well-deliquification | ||
+ | |url-access=subscription | ||
+ | }}</ref> | ||
+ | |||
+ | |||
+ | Lea, J. F., Lea Jr., J. F., Nickens, H. V. (2011). Gas Well Deliquification. Нидерланды: Elsevier Science. | ||
+ | </references> | ||
− | + | [[Category:pengtools]] | |
− | + | [[Category:PQplot]] | |
− | + | {{#seo: | |
− | + | |title=Liquid loading Turner velocity | |
+ | |titlemode= replace | ||
+ | |keywords=Liquid loading, gas flow, Turner velocity, Trammel, Coleman, petroleum engineering | ||
+ | |description=Liquid loading is a phenomenon when the gas phase does't provide sufficient velocity to lift the liquids out of the gas well. | ||
+ | }} |
Latest revision as of 10:45, 8 May 2024
Contents
Brief
Liquid loading is a phenomenon when the gas phase does't provide sufficient transport energy to lift the liquids out of the well.
In 1969 Turner et al. published an empirical correlation defining the Liquid loading gas velocity[1].
Math & Physics
The minimum gas velocity to remove the liquid equation:
- [2] (1)
Note that original paper[1] uses 20.4 constant which is good if the surface tension units are in 'lbf/ft'. Also 20.4 is calculated by adjusting the theoretical equation constant by 20 percent. For the reference PROSPER uses 2.04. PQplot uses 1.593.
The minimum gas rate to remove the liquid equation:
- (2)
PQplot workflow
In order to calculate the liquid loading point in PQplot:
For every point in the VLP curve:
1. Critical gas velocity at the wellhead is calculated by eq (1)
2. Wellhead flowing gas velocity is calculated
3. Liquid loading flag is 1 if:
4. Critical gas rate is the last VLP curve point at which liquid loading flag is 1.
Discussion
To avoid the Liquid loading the gas velocity should be above the Liquid loading velocity.
The higher the gas rate the higher the gas velocity.
The lower the wellhead flowing pressure the higher the gas rate.
The bigger the tubing ID the higher the gas rate.
In case when the gas rate is limited by the Reservoir deliverability smaller tubing ID will increase the gas velocity.
Example. Calculating critical gas velocity and rate
Example 3-2 from [2]
Input data
- , flowing wellhead pressure, psia
- , flowing wellhead temperature, °F
- , liquid density, lbm/ft3
- , water surface tension, dyne/cm
- , gas specific gravity
- , gas compressibility factor
Production string = 2-3/8 inch tubing with 1.995 in ID, A = .0217 ft2
Water is the produced liquid
Solution
Calculate the gas density:
- , lbm/ft2
Calculate the critical gas velocity form (1):
- , ft/sec
Calculate the critical gas rate form (2):
- , MMscfd
Solution with PQplot
PQplot calculates the following critical gas rate:
- , MMscfd
Liquid loading example 3-1 (Public model)
Nomenclature
- = flow area, ft^2
- = flowing wellhead pressure, psia
- = gas rate, MMscf/d
- = gas density, lbm/ft3
- = liquid density, lbm/ft3
- = surface tension, dyne/cm (ref values: 60 - water, 20 - condensate) [1]
- = gas specific gravity
- = flowing wellhead temperature, °R
- = gas velocity, ft/sec
- = gas compressibility factor at flowing P & T, dimensionless
See also
References
- ↑ 1.0 1.1 1.2 Turner, R. G.; Hubbard, A. E.; Dukler (Nov 1969). "Analysis and Prediction of Minimum Flow Rate for the Continuous Removal of Liquids from Gas Wells". Journal of Petroleum Technology (SPE-2198-PA): 1475–1482.
- ↑ 2.0 2.1 Lea, J. F.; Nickens, M.; Wells (2011). Gas Well Deliquification. Energy Tribune Publishing Inc.