7.

FLUJO DE GAS BAJO CONDICIONES ESTABLES

Nomenclatura
Simbolo
A
C
D
d
E
e

f
g

rea
Constante de conversin de unidades
Dimetro
dimetro en unidades prcticas
Eficiencia de flujo
Rugosidad absoluta
Rugosidad relativa
Factor de friccin
Aceleracin de la gravedad

h
L
l
M
m
NRe
P
Pb
q
q
qb
qsc

entalpa especfica
Longitud
Longitud en unidades prcticas
Peso molecular
Masa
Nmero de Reynolds
Presin
Presin base
Tasa de flujo
Calor
Tasa de flujo a condiciones base
Tasa de flujo a condiciones base en unidades
prcticas
Constante universal de los gases
Entropa
Tiempo
Temperatura
Velocidad, flujo volumtrico
Volumen especfico
Trabajo
Trabajo especfico
Distancia
Factor de compresibilidad de los gases
altura

R
S
t
T
u
v
W
w
x
Z
z

Dimensiones
L2
L

L
adimensional
adimensional
32,2 ft/s2
9.8 m/s2
L2/t2

adimensional
M/Lt2
L3/t
ML2/t2
L3/t

ML2/t2T

L/t
L /m
mL2/t2
L2/t2
L
adimensional
L
3

Introduccin
La tubera es el medio de transporte ms comn para llevar el gas de un sitio a otro, aunque
cuando ya las distancias son demasiado grandes en algunos casos se tiene que transportar el
gas en tanqueros (metaneros). El gas puede viajar por tuberas a presiones altas o bajas, en
el primer caso por ejemplo cuando se tienen redes de recoleccin, gasoductos o redes de
inyeccin en sistemas de bombeo neumtico, en el segundo caso, por ejemplo, cuando se
tienen redes de distribucin de gas.

7.1 Ecuacin General para Flujo de Gas en Tuberas (2, 12)

El problema de flujo en tuberas para cualquier fluido se puede analizar partiendo de la
primera ley de la termodinmica o ley conservacin de energa, la cual se puede expresar en
unidades de energa por libra masa como:
1
h u 2 gz q w
(7.1)
2
donde:
h
:
Cambio en entalpia
2
u
:
Cambio en el cuadrado de velocidad
z
:
Cambio en altura
q
:
Calor entregado o recibido por el fluido
w
:
Trabajo realizado por o sobre el fluido
Para el caso de flujo de gas, se pueden hacer las siguientes suposiciones:
u2= 0
z= 0
T= Constante (flujo isotrmico)
w= 0
o sea que la ecuacin (7.1) queda como:
h q

(7.2)

Pero de acuerdo con la termodinmica:

h Ts vP
o en forma diferencial:
dh Tds vdP
y
Tds dq dLw
o sea que llevando estas tres ltimas expresiones a la ecuacin (7.2) se tiene:
(7.3)
vdP dLw 0
donde dLw se conoce como las prdidas irreversibles de energa ocasionadas, por ejemplo,
por friccin, y v y P son volumen especfico y presin respectivamente .

Se tienen varias expresiones para calcular dLw, una de las ms conocidas es la ecuacin de
Moody la cual tiene la siguiente forma:
dL W fu 2

(7.4)
dL
2D
donde:
u
:
velocidad del fluido,
L
:
longitud a travs de la cual ocurre las prdidas de energa,
D
:
Dimetro de la tubera
Llevando la ecuacin (7.4) a la ecuacin (7.3) se tiene:
fu 2dL
vdP
0
(7.5)
2D
pero volumen especfico es el inverso de densidad y de acuerdo con la expresin para
calcular la densidad de un gas:
ZRT
v
PM
donde P y T son las condiciones de presin y temperatura a las que se encuentra el gas, M
es el peso molecular del mismo , R la constante universal de los gases cuyo valor depende
de las unidades usadas para las variables de la ecuacin de estado y Z es el factor de
compresibilidad..
Adems:
q
u
A
y si se expresa q en trminos de volumen medido a condiciones base se tiene:
P T Z
q qb b
P Tb Zb
donde el subndice b se refiere a condiciones base y qb es la tasa de flujo de gas medida a
condiciones base; o sea que entonces velocidad queda como:
P T Z 1
u qb b
P Tb Zb A
P T Z 4
u qb b
P Tb Zb D2
y llevando las expresiones anteriores a la ecuacin (7.5), recordando que M=29g y Zb =1,
se tiene:
2
2
2
2
ZRT
f 4
2 Pb T Z
dP

*q

5 dL
b

P
*29

g
b P D

Simplificando y separando variables, se puede escribir:

f q P
0,8106 5 b b
D Tb

R 1
dL
29 g

p1

PdP
ZT
p2

y despus de integrar y despejar qb

0,5
R

Tb D5 P1 PdP
qb

29*0,8106 Pb f g L P2 ZT
la cual , finalmente se puede dejar como

0,5

0,5

T D5 P1 PdP
(7.6)
qb C b

Pb f g L 2 ZT
donde C es una constante que depende de las unidades usadas para las variables. Cuando se
usan unidades absolutas del sistema internacional (SI), la constante vale 18,8 y cuando se
usan unidades absolutas del sistema ingls la constante vale 46,1

La ecuacin (7.6) es la forma general de la ecuacin de flujo para gas en tuberas

suponiendo flujo horizontal y en estado estable. Se conoce como la ecuacin de Clinedinst
y para aplicarla se requiere resolver el integral el cual se transforma en trminos de la
presin y temperatura seudorreducida y su valor se puede obtener de tablas existentes en la
literatura.(ver referencia 1)
Una forma ms comn de la ecuacin general par flujo de gas en tuberas se obtiene de la
ecuacin (7.6) tomando las variables Z y T como valores promedios y constantes , lo cual
permite sacarlas del integral y adems efectuar analticamente ste ; la ecuacin queda de la
siguiente forma:
0,5

T d5 *(P12 P22 )
(7.7)
q b Cg b

Pb f g LZT
donde Z y T son los valores promedios de Z y T tomados como constantes, Cg es una
constante que depende de las unidades de las variables. La tabla 27 muestra valores de Cg
para diferentes grupos de unidades en los sistemas SI e Ingls.

Con respecto a la ecuacin (7.7) se debe hacer claridad sobre los factores Z y f. Z es el
factor de compresibilidad calculado a condiciones promedias, P y T . Con respecto aT,
normalmente se considera flujo isotrmico y por tanto T = T = temperatura de flujo; en
cuanto a la presin sta s vara ampliamente y por lo tanto se deben proponer formas de
calcular un valor promedio, y entre las ms conocidas se tienen
P P
P 1 2
(7.8)
2
3
3
2 P1 P2
(7.9)

P 2
3 P1 P2 2

Tabla 27. Valores de la constante Cg en la ecuacin (53) para diferentes grupos de unidades en los sistemas
ingls y S.I.

Sistema
Ingls

Sistema
Mtrico
(S.I)

Grupo
qh
P
Unidades
1
PC/s lb/pie.S
2
2
PC/hr
3
MPC/D
lpca
4
KPC/D
lpca
lpca
1
m3/s
Pa
3
2
m /hr
kPa
3
m3/hr
kPa
3
4
m /D
kPa
5
m3/D
bars

R
R
R
R

pie
pulg.
pulg.
pulg.

pie
milla
milla
pie

32,64
3,23
7,75210-5
5,633

K
K
K
K
K

m
cms
mm
mm
mm

m
m
km
km
km

13,36
0,4786
4,78610-5
1,1510-3
1,1510-3

o tambin se puede calcular Z de

Z Z2
Z 1
2

(7.10)

La ecuacin (7.9) se obtiene de la siguiente manera (2):

La presin promedia tambin se puede plantear de la siguiente forma: Supongamos que se
tuviera un grfico del comportamiento de la presin con la distancia, o sea:

De acuerdo con el Teorema del valor medio:

L

P(x)dx

(7.11)

X2

o sea que si se tuviera una expresin para P(x) se podra calcular la integral P entre 0 y X2

Una expresin para P(x) se puede obtener as: Supongamos un punto cualquiera x entre 0 y
X2; en tal punto se puede establecer aplicando la ecuacin general de flujo en tuberas
(ecuacin (7.7)):

T 1
q b Cg b
Pb ZT g f

T 1
Cg b
Pb ZT g f

1/2

1/2

P12 P 2 5
d

x

1/2

P 2 P22 5
d

X 2 x

1/2

De estas dos expresiones se puede obtener:

P12 P 2 P 2 P22
, donde L = X2

x
Lx
1/2
2 x 2
2
P P1 P1 P2
L

(7.12)

Llevando la ecuacin (7.12) a la (7.11) se tiene:

0.5
L
x 2
2
2
3/2 L
0 P1 L P1 P2 dx 2
1
1
2 x 2

P
L 2
P P1 P22
2 1
L
3 P1 P2
L
0 L

3/2
2 1
P 2 P12 P22 P12
2
2 1

3 P1 P2

2 P13 P23
3 P12 P22

(7.9)

1
1
se conoce como factor de transmisin.
f
f
El factor de friccin f, depende de un parmetro conocido como nmero de Reynolds y de
la rugosidad de la tubera, aunque cuando el nmero de Reynolds es grande f, depende
solamente de la rugosidad de la tubera.
El factor

El nmero de Reynolds se define como:

uD
N Re

donde:

:
Densidad del fluido
u
:
Velocidad del fluido
D
:
Dimetro de la tubera

:
Viscosidad del fluido

El trmino u se conoce como flujo msico que es tasa msica por unidad de rea y se
representa por G, o sea que la expresin anterior quedara como:
GD q g *D
N Re

A
y cuando se da la tasa de flujo en trminos de qb y se reemplaza por su definicin a partir
de la ecuacin de estado de los gases , se tiene
g
P
N Re C b q b
(7.13)
Tb D
donde C es una constante que depende de las unidades usadas para las variables y las
variables conservan las definiciones que se han venido dando en el texto. Cuando se usan
unidades absolutas del sistema ingls el valor de C es 7.42*10-4 y cuando se usan unidades
absolutas del sistema SI el valor de C es 0.0044, y cuando se usan las siguientes unidades
de campo del sistema ingls: P en Lpc., q en KPCN/D, en cp. , y D en pulgadas, la
constante es 711.8
La expresin para NRe dada por la ecuacin (7.13) es para el caso particular en el que el
conducto sea circular. La expresin general para Nre es
ud e
N Re
(7.14)

donde de es el dimetro equivalente y es igual a 4 veces el radio hidrulico definido por

Area de Flujo
Rh
Perimetro Mojado
Para el caso de un conducto circular:
AT = d2/4
2r = d = Permetro mojado.
2
d
4 d
Rh
d
4
de = 4Rh = d
para el caso de flujo anular:
d o2 di2
rea de Flujo:
4
Permetro Mojado : do di
d o2 d i2

Rh

4 1 d o di d o di
d o di
4
d o di

d o di

(7.15)

4
de do di 4*R h

De acuerdo con el valor del nmero de Reynolds se define el rgimen de flujo y de acuerdo
con ste ltimo se tienen expresiones para calcular f, de la siguiente manera (2).

Flujo laminar (NRe < 2.000)

1
0,125N 0Re,5
f

(7.16)

Flujo crtico (2.000 < NRe < 4.000)

1
1, 4142N 0,15
Re
f
1,16

Flujo de transicin 4.000 <N Re < 200

e
1
9,34
1,14 2log

f
D N Re f
1,16

Flujo turbulento N Re > 200

1
e
1,142log
D
f

(7.17)

(7.18)

(7.19)

En las expresiones anteriores D es el dimetro de la tubera y e es la rugosidad absoluta o

sea la magnitud de las irregularidades que se presentan en la superficie de la tubera.
Ikoku (8), plantea la siguientes expresiones para calcular el factor de friccin considerando
que solo hay dos regmenes de flujo: laminar y turbulento. La expresin para el flujo
laminar es la ecuacin (7.16), pero cuando se tiene flujo turbulento la expresin para f
depende de si se trata de tuberas suaves o rugosas.
Para tuberas suaves:
f = 0,0056 + 0,5N-0,32
Re

(7.20)

La ecuacin (7.20) se conoce como ecuacin de Drew - Koo y McAdams y se aplica para
NRe entre 3103 y 3106.
Cuando se trata de tuberas rugosas se utiliza la ecuacin de Colebrook and White:
2e
1
18, 7
1, 74 2log

f
D N Re * f

(7.21)

La ecuacin (7.21) tiene la desventaja de que no es explcita en f y por tanto para hallarlo se
debe recurrir a ensayo y error.
Una ecuacin posterior, de (1976), la de Colebrook es de (1939), es la de Jain(9)

e 21,25
1,14 2log 0,9
f
D N Re

(7.22)

Las ecuaciones (7.21) y (7.22) son las ms usadas para obtener f.

En el caso de flujo de gas se acostumbra generalmente hablar, para valores de N Re > 2000,
de flujo parcialmente turbulento y flujo totalmente turbulento(10). Se habla de flujo
parcialmente turbulento cuando cerca a la pared de la tubera hay una zona de fluido donde
an permanece el flujo laminar y en la parte central de la tubera hay flujo turbulento;
cuando desaparece la zona de flujo laminar se habla de flujo totalmente turbulento.
En general, en la zona de flujo parcialmente turbulento el factor de friccin depende del
nmero de Reynolds y en la zona totalmente turbulento depende de la rugosidad relativa de
la tubera y para un valor dado de esta variable se mantiene constante con el nmero
Reynolds.
El rgimen parcialmente turbulento est muy asociado con tuberas lisas (de baja
1
rugosidad) y para obtener
se plantea una expresin conocida como ley de flujo en
f
tuberas lisas, dada por
N
1
4 log Re
(7.23)
f
1
f
Por otro lado el flujo totalmente turbulento est asociado con tuberas rugosas (rugosidad
alta) y el factor de transmisin 1 se calcula con una expresin conocida como ley de
f
flujo en tuberas rugosas, dada por
1
4 log 3,7
(7.24)

f
donde es la rugosidad relativa de la tubera.
Para aplicar la ecuacin (7.23) o (7.24) es necesario definir si se tiene flujo parcialmente
turbulento o totalmente turbulento y aunque existen procedimientos para ello(10), un criterio
aproximado es que para valores de NRe entre 4.000 y 500.000 se puede considerar flujo
parcialmente turbulento y para valores de NRe mayores de 500.000 se considera flujo
totalmente turbulento.

Es importante adems observar que para calcular f usando la ecuacin (7.23) se debe
recurrir a un proceso de ensayo y error, y para aplicar la ecuacin (7.24) se requiere
conocer ; las referencias (2), (9) y (10) muestran criterios para determinar este valor.
Serghides (10) partiendo de la ecuacin de Colebrook y aplicando un mtodo numrico
iterativo para su solucin obtuvo una expresin para calcular directamente el factor de
friccin f.
La ecuacin de Colebrook usada por Serghides fu
e
1
2,51
D

(7.25)
2 log

3, 7 N Re f
f

y llego a la siguiente expresin de tres parmetros, que es vlida para NRe > 2100 y
cualquier valor de rugosidad relativa (/D):
2

B A

f A

2B

donde,
e
12
A 2 log D
3, 7 N Re

e
2,51A
B 2log D
N Re
3, 7

e
2,51B
D

C 2 log

N Re
3, 7

(7.26)

(7.27)

(7.28)

(7.29)

Serghides adems hizo un estudio comparativo de la ecuacin (7.25) y otras 7 ecuaciones

comunes usadas para calcular f y encontr que la ecuacin propuesta por l arrojaba
desviaciones promedia y mxima de 0.0002% y 0.0023% respectivamente; mientras las
dems ecuaciones la de mejor comportamiento mostraba para estas mismas variables
valores de 0,027% y 0,138% respectivamente.
Finalmente, f se puede obtener de grficos existentes en la literatura como los que se
muestran en las figuras 82 y 83, pero se debe tener claro que cuando se vaya a determinar f
de grficos (Moody o Fanning) es necesario saber de cual grfico se trata pues fMoody = 4
fFanning; el criterio de aclaracin es que en el grfico de Moody f = 0.064 cuando N Re =
1.000 y en el de Fanning f = 0,016 cuando NRe = 1.000.
Cuando no se conoce la rugosidad de la tubera se recomienda los siguientes valores para
las rugosidades absolutas en pulgadas:
Tubera de Produccin nueva: t = 0,0006

10

Figura 82. Factor de Friccin de Moody.

11

Figura 83. Factor de Friccin de Fanning

Lnea de Flujo:
Tubera Galvanizada
Tubera Recubierta

t = 0,0007
t = 0,006
t = 0,01 0,1

Ejemplo 7.1
Un gas de g = 0.7 y g = 0.2 cp fluye a 30 pies/s a travs de una tubera de 6 pulgadas de
dimetro. La tubera es de acero comercial y es nueva.
Calcular f por:
Ecuacin de Colebrook
Ecuacin de Jain
Ecuacin de Serghides

12

Grficos

Solucin
Calculo de qg (Suponiendo flujo a condiciones normales)
86.400
q g u *A = 30*(d 2 / 4)
1.000
2
6 86.400
= 30*
4 12 1.000
508,94 KPCN / D
Clculo de Rugosidad
Para una tubera nueva de acero comercial se toma como rugosidad absoluta 0.00015 pies;
o sea que:
e 0,00015

0,0003
6
D
12

Clculo por Mtodo de Jain

e 21,25
1
1,14 2 * Log 0,9
f
D N Re
El nmero de Reynolds, usando ecuacin (7.13), es:
g
P
N Re 711,8 b *q g
Tb
D
14, 7
0, 7
711 *
*508,94*
520
0, 02*6
59671,6 Flujo Turbulento
1
21, 25

1,14 2log 3*104

59.671, 60,9
f

f 0, 0212

Mtodo de Colebrook:
Supongamos inicialmente para f un valor 0,021, y el valor calculado ser, usando
ecuacin (65):

18, 7
1, 74 2 log 2 e

D N Re 0, 021

= 0,0213

13

Mtodo de Serghides:
Usando las ecuaciones (7.25) - (7.28) se tiene:
e
12
A 2log D
7,1
3, 7 N Re

2,51A
B 2log D
6,8439
N Re
3, 7

e
2,51B
C 2log D
6,8666
N Re
3, 7

B A 2 0,02122
f A
C 2B A

7.2. Ecuaciones Prcticas para Flujo de Gas (12)

Observando la ecuacin general para el flujo de gas en tuberas ecuacin (7.7) y el
1
procedimiento para calcular f
se encuentra que la aplicacin de esta ecuacin es un
f
proceso de ensayo y error (si se desconoce qb y/o d), por lo tanto para simplificar un poco
esta ecuacin se han propuesto diferentes alternativas para manejar el trmino relativo al
factor de friccin.
El trmino 1/ f de la ecuacin (7.7) se ha llamado factor de transmisin y de acuerdo a la
forma como se considere este factor se tienen tres tipos de ecuaciones as:
Factor de transmisin constante
Factor de transmisin en funcin del dimetro de la tubera
Factor de transmisin en funcin del nmero de Reynolds
El primer caso ser para cuando se tiene flujo laminar, el segundo caso, como ya se
mencion , ser aplicable al rgimen totalmente turbulento y el tercero al rgimen
parcialmente turbulento.
Esto ha dado origen a varias ecuaciones para analizar el flujo de gas en tuberas, y se vern
a continuacin algunas de las ms usadas en las diferentes situaciones donde es comn
trabajar con flujo de gas: redes de recoleccin y gasoductos, redes de distribucin y redes
domiciliarias.

14

7.2.1. Ecuacin de Weymouth.

Se obtiene reemplazando el factor de transmisin en la ecuacin (7.7) por la siguiente
expresin:
1
C*d1/6
(7.30)
f
donde d es el dimetro de la tubera y C es una constante que depende de las unidades de d.
Cuando d est en pies, C = 8.46 y cuando d est en metros, C = 10.31
Llevando la ecuacin (7.30) a la ecuacin (7.7) se obtiene la siguiente expresin general
para la ecuacin de Weymouth:
T d16/3
q b Cw b

Pb g ZTL

0,5

2
1

P22

0,5

(7.31)

donde Cw es la constante de la ecuacin de Weymouth y las dems variables ya fueron

definidas en la ecuacin (7.7).
El valor de Cw depende de las unidades usadas para las variables . La tabla 28 muestra el
valor para esta constante cuando se usan diferentes grupos de unidades en los sistemas
ingls e internacional.
De acuerdo con la ecuacin (7.30) la ecuacin de Weymouth se aplica a flujo totalmente
turbulento y generalmente se recomienda para presiones medias ( entre 100 y 500
lpca.(690.5 y 3452 Kpa.)) y dimetros medios (entre 6 y 12 pulgs. (152 y 305
mm)),situacin que normalmente se da en redes de distribucin primarias.
De acuerdo con la ecuacin (7.30) la ecuacin de Weymouth se aplica a flujo totalmente
turbulento y generalmente se recomienda para presiones medias (entre 100 y 500 lpca.
(690.5 y 3452 Kpa.)) y dimetros medios (entre 6 y 12 pulgs. (152 y 305 mm)), situacin
que normalmente se da en redes de distribucin primarias.
Tabla 28. Constante de la Ecuacin de Weymouth para Diferentes Grupos de Unidades en los Sistemas Ingls
y S.I. de Unidades.

Grupo
qh
P
Unidades
1
PC/s
Lbm/pie.s2
2
MPC/D
lpca
Sistema Ingls
3
KPC/D
lpca
4
PC/hr.
Lpca.
1
m3/s
Pa
3
2
m /hr
kPa
Sistema
3
m3/hr
kPa
Internacional,
3
4
m /D
kPa
S.I
5
m3/D
bars.

15

Cw

R
R
R
R
K
K
K
K
K

pie
millas
millas
millas
m
m
kms
kms
kms

pie
pulgs.
pulgs.
pulg.
m
cms
mm
mm
mm

276.1344
4.33*10-4
4.33*10-1
18.062
137.75
2.301
1.568*10-4
3.76*10-3
3.76*10-3

7.2.2. Ecuacin de Panhandle.

Esta ecuacin aplica para

una expresin en trminos del nmero de Reynolds, y

f
mientras que la ecuacin de Weymouth se aplica a condiciones de presiones medias y
dimetros medios, la de Panhandle se recomienda para presiones altas (mayores de 500
lpca. (3452 Kpa.)) y dimetros grandes (mayores de 12 pulgs. ( 305 mm.)) .
Se tiene dos versiones de la ecuacin de Panhandle conocidas como ecuaciones A y B de
Panhandle, como la B es la ms reciente se mostrar esta ltima.
La expresin propuesta por Panhandle para

1
f

es

1
16, 7* N Re 0,01961
(7.32)
f
donde NRe es el nmero de Reynolds que para el caso de flujo de gas en tuberas se puede
calcular de acuerdo con la ecuacin (7.13).
De acuerdo con la ecuacin (7.32), la ecuacin de Panhandle es aplicable al rgimen
parcialmente turbulento ya que f es funcin de Nre.
Llevando la ecuacin (7.13) a la (7.32) y luego esta a la ecuacin (7.7) se obtiene la
siguiente expresin como la ecuacin de Panhandle
0,51

T d 2,53 P 2 P22
1
q b Cp b 0,49 1
(7.33)

0,020
Pb g ZTL g
donde Cp es la constante de la ecuacin de Panhandle y su valor depende de las unidades
utilizadas para las variables. La tabla 29 muestra valores de Cp para diferentes grupos de
unidades en los sistemas ingls y S.I.
Tabla 29. Valores para Cp en la Ecuacin (7.33) para diferentes Grupos de Unidades en los Sistemas Ingls
y S.I. (mtrico).

Sistema
Ingls

Sistema
S.I

Grupo
qh
P
Unidades
1
PC/s Lbm./pie.s2
2
PC/hr
lpca
3
MPC/D
lpca
4
kPC/D
lpca
3
1
m /s
Pa
2
m3/hr
kPa
3
m3/hr
kPa
3
4
m /D
kPa
5
m3/D
bars

16

R
R
R
R
K
K
K
K
K

pies
millas
millas
millas
m
m
kms
kms
kms

pies lb/pie.seg
pulg.
cp.
pulg.
cp.
pulg.
cp.
m
Pa.s
cms
mPa.s
mm
mPa.s
mm
mPa.s
mm
mPa.s

C
535.2985
62.059
1.4910-3
1.49
222.88
9.2153
8.026*10-4
1.926*10-2
2.11*10-2

7.2.3. Ecuacin de Spitzglass

Esta ecuacin es aplicable a casos de cada de presin baja y presiones bajas (menores de
30 lpc.(207 Kpa.)) asociadas a tasas de flujo y tamaos de tubera bajos ( menores de 4
pulgs. (10.2 mm.)).
Aplica para

1
la siguiente expresin
f
0,5

1
100

(7.34)
f 1 C1 C *d
2
d

donde d es el dimetro de la tubera y C1 y C2 son constantes que dependen de las unidades

utilizadas para el dimetro. En el sistema ingls C1 = 3,6 y C2 = 0,03 cuando el dimetro
est en pulgadas y C1 = 0.3 y C2 = 0.036 cuando el dimetro est en pies. En el sistema S.I.
C1 = 91,44 y C2 = 1,18*10-3 cuando el dimetro est en mm, C1 = 9,144 y C2 = 1,18*10-2
cuando el dimetro est en centmetros y C1 = 0,09144 y C2 = 1,18 cuando el dimetro
est en metros.

Algunos autores presentan la ecuacin (7.34) con una constante de 88,5 en lugar de
100.(ref.10, pag.88)
Segn la ecuacin (7.34), la ecuacin de Spitzglass es aplicable al rgimen completamente
turbulento.
Spitzglass manipula la ecuacin (7.34) de la siguiente manera:
P12 P22 P1 P2 P1 P2

P1 P2 P
supone que
P2 0,9 P1 y por tanto (P1 +P2) 1,9P1, T = 520R en el sistema ingls y 298,8K en
sistema S.I., Z=1.
Y reemplaza la cada de presin por
P w *gh w
donde hw es la altura de una columna de agua que ejerce una presin igual a P, w es la
densidad del agua y g es la aceleracin debida a la gravedad.
Aplicando las modificaciones y suposiciones anteriores y llevando la ecuacin (7.34) a la
ecuacin (7.7) y recordando los valores para w y g en los sistemas ingls y S.I de unidades,
se obtiene la siguiente expresin general conocida como ecuacin de Spitzglass:

17

0,5

Tb
h w * P1 *d
(7.35)

q b Cs
C
P

1
b * L 1 C *d

2
g
d

donde,
qb, Tb, Pb, d, g, L y P1 tienen la misma definicin dada en la ecuacin (7.7) y hw es la cada
de presin (P1 - P2) dada en el equivalente a una columna de agua. Cs es la constante de la
ecuacin de Spitzglass y su valor depende de las unidades para las variables; C1 y C2 son
constantes, cuyos valores fueron dados dependiendo de las unidades usadas para el
dimetro.

La tabla 30 muestra valores para Cs, para diferentes grupos de unidades en los sistemas
ingls y S.I
Tabla 30. Valores de Cs en la Ecuacin (7.35) para diferentes Grupos de Unidades en los Sistemas Ingls y
S.I. (mtrico)

Grupo
qb
P
Unidades
1
PC/s
Lbm./pie.s2
Sistema
2
PC/hr
lpca
Ingls
3
PC/s
lpca
4
MPC/D
Lpca.
5
kPC/D
Lpca.
3
1
m /s
Pa
Sistema
2
m3/hr
kPa
S.I
3
m3/hr
kPa
4
m3/hr
bars
5
m3/hr
mbars

hw

R
R
R
R
R
K
K
K
K
K

pie
pulg.
pulg.
pulg.
pulg.
m
cm
mm
mm
mm

pie
millas
pie
millas
millas
m
m
kms.
m
m

pie
pulg.
pulg.
pulg.
pulg.
m
cm
mm
mm
mm

C
884.39
0.372
7.5110-3
8.93*10-6
8.93*10-3
1054.64
0.1201
3.8*10-6
1.201*10-5
3.8*10-4

7.2.4. Ecuacin de Mueller

Esta ecuacin es similar a la de Panhandle, porque considera el factor de friccin de Moody
(f) como una funcin del nmero de Reynolds; pero a diferencia de la ecuacin de
Panhandle, se aplica a situaciones de presin baja y dimetros pequeos.
La expresin propuesta por Mueller para el factor de transmisin es:
1
1, 675N Re 0,130
f
donde f es el factor de friccin de Moody y NRe es el nmero de Reynolds.

(7.36)

Recordando la expresin para el nmero de Reynolds dada por la ecuacin (7.13), llevando
esta a la ecuacin (7.36) y luego la expresin resultante a la ecuacin (7.7) se obtiene la
siguiente expresin para la ecuacin general de Mueller
18

T
d 2,725 P 2 P22
q b CM b 0,425 0,15 1

Pb g g ZTL

0,575

(7.37)

donde, todas las variables con excepcin de g fueron definidas en la ecuacin (7.7) y g es
la viscosidad del gas; y CM es una constante que depende del sistema de unidades usado.
La tabla 31 muestra valores de CM para diferentes grupos de unidades en los sistemas ingls
y S.I.
Tabla 31. Valores de CM en la Ecuacin (7.37) para diferentes Grupos de Unidades en los Sistemas Ingls y
S.I. (mtrico).

Sistema
Ingls

Sistema
S.I

Grupo
qh
P
Unidades
1
PC/s
lbf/pie2
2
PC/hr
lpca
3
kPC/hr
lpca
4
MPC/D
lpca
5
PC/s
lb/pie.s2
1
m3/s
Pa
3
2
m /hr
kPa
3
m3/hr
kPa
4
m3/D
kPa
3
5
m /D
bars

R
pies
pies
R millasp pulgs.
R
ies
pulgs.
R millas pulgs
R
pies
pies
K
m
m
K
kms
cm
K
kms
mm
K
kms
mm
K
Kms
mm

lb/pie.s
CP
lb/pie.s
CP
lb/pie.s
Pa.s
mPa.s
mPa.s
mPa.s
mPa.s

CM
56,965
10,750
0,4966
2,5810-4
33,84
15,83
0,0302
5,7410-5
1,3810-3
2,7510-3

7.2.5. Ecuacin de Mueller para Baja Presin.

De las ecuaciones vistas, la de Spitzglass, como se plante, es una ecuacin para aplicar a
presiones bajas. Generalmente una presin baja est asociada a un tamao de tubera
pequeo.
Cualquiera de las ecuaciones vistas se puede modificar para llevarla a una ecuacin para
presiones bajas de la siguiente manera:
El trmino P12 P22 se puede escribir como:

2
1

P22 P1 P2 P1 P2

P1 P2

P1 P2
2
(7.38)
2P P1 P2
y el termino (P1 - P2) de la ecuacin (16) se puede expresar en trminos de la altura de una
columna de agua que ejerce una presin igual a (P1 - P2) as
(7.39)
P1 P2 w *g *h w
2

19

donde w es la densidad del agua, hw es la altura de la columna de agua que ejerce una
presin igual a (P1 - P2), g es la aceleracin de la gravedad cuyo valor en el sistema ingls
equivale a 32.2 pies/s2 y en el S.I. a 9.8 m/s2.
Recordando que en el sistema ingls w es 62.4 lbs/pie3 y en el S.I. 1000 Kg/m3, y llevando
la ecuacin (7.39) a la (7.38) se tiene
(7.40)
P12 P22 C*P*h w
donde C es una constante que en el sistema ingls vale 4018.56 y en el S.I. 19600 cuando
se tienen unidades absolutas para ambos sistemas.
La ecuacin (7.40) se puede llevar a cualquiera de las ecuaciones vistas y as se obtiene la
versin de dicha ecuacin para presiones bajas.
De esta manera, la ecuacin de Mueller para presiones bajas tiene la siguiente forma
general
0,575

Ph
T
d 2,725
(7.41)
q b C b 0,425 0,15 * w
Pb g g
TL

donde se ha asumido que Z = 1, por ser una ecuacin para presiones bajas, y C 'M es la
constante de la ecuacin de Mueller para presiones bajas cuyo valor depende de las
unidades usadas.
'
M

La tabla 32 muestra valores de C 'M para diferentes grupos de unidades en los sistemas
ingls y S.I.
Tabla 32. Valores de la Constante
Sistemas Ingls y S.I.

Sistema
Ingls

Sistema
S.I

Grupo
Unidades
1
2
3
4
1
2
3
4
5

qh

C 'M en la Ecuacin (7.41) para diferentes Grupos de Unidades en los

PC/s
lbf/pie2
PC/hr
lpca
PC/D
lpca
PC/s Lb/pie.s2
m3/s
Pa
m3/hr
kPa
m3/hr
kPa
m3/hr
mbars
m3/D
bars

hW

C'M

R
R
R
R
K
K
K
K
K

pies
pies
pies
pies
m
m
m
m
m

pies
pulgs
pulgs
pies
m
cms
mm
mm
mm

lb/pies
CP
CP
Lb/Pie.s
Pa.s
mPa.s
mPa.s
mPa.s
mPa.s

pies
pulg.
pulg.
pie.
m
cms
mm
mm
mm

911,803
326,2
7,83103
3997,74
4,65*103
0,63
3,1510-4
8,3810-4
1,06*10-3

20

7.3. Seleccin de la Ecuacin Adecuada

Como ya se menciono, antes de entrar a realizar clculos de flujo de gas en tuberas es
necesario definir el rgimen de flujo existente: laminar o turbulento y en este ltimo caso si
parcial o totalmente turbulento.
Cuando NRe < 2000 se tiene flujo laminar y en este caso

1
f

se calcula usando la ecuacin

(7.16) y este valor se lleva a la ecuacin (7.7).

Cuando NRe > 2.000 se puede hablar de flujo turbulento parcial o total. Aunque es posible
determinar si es parcial o totalmente turbulento, se considera que para NRe < 400000 se
tiene parcialmente turbulento y cuando NRe > 400000 se tiene flujo totalmente turbulento.
Una vez definido si es flujo parcial o totalmente turbulento, se podra usar la ecuacin
1
(7.23) en el primer caso y la (7.24) en el segundo caso para calcular
y luego usar la
f
ecuacin (7.7).
En lugar del procedimiento anterior se puede pensar en usar alguna de las ecuaciones
presentadas en el numeral 1.2, algunas de las cuales son para flujo parcialmente turbulento
(Panhandle, Mueller) y otras para flujo totalmente turbulento (Weymouth, Spitzglass).
En general se puede decir lo siguiente (3, 10, 12)
En redes de recoleccin y gasoductos (presiones altas y dimetros grandes) usar
Panhandle.
En redes de distribucin primarias (presiones medias y dimetros medios), usar
Weymouth para flujo totalmente turbulento y Mueller para totalmente turbulento.
En redes de distribucin secundarias (presiones bajas y dimetros pequeos) usar
Mueller de baja presin cuando se tiene flujo parcialmente turbulento o Spitzglass para
flujo totalmente turbulento.
En redes domiciliarias, los mismos criterios de las redes secundarias.

7.4. Eficiencia de Tubera (E)

Las ecuaciones planteadas para flujo de gas suponen que toda la tubera est ocupada por
gas; esto puede no ser cierto por la presencia de lquidos en el gas o la presencia de
escamas, slidos o productos de corrosin en la paredes de la tubera. Esto hace que el lado
derecho de las ecuaciones (7.7) y cualquiera otra de las ecuaciones para flujo de gas
presentadas en el numeral 1.2 se deba multiplicar por un factor E conocido como eficiencia
de la tubera y el cual es menor de 1. Un valor para E de 0.9 se considera satisfactorio y se
podr usar cuando el gas a transportar es seco y la tubera es nueva.

21

Normalmente y con el fin de mantener la eficiencia de las tuberas removiendo el lquido

condensado (hidrocarburos y agua) y los slidos depositados se envan raspadores
(conocidos como marramos, y el proceso se denomina marraneo (pigging)).
Suponiendo tuberas nuevas, se recomiendan los siguientes factores de eficiencia
dependiendo de tipo de gas que se est transportando:
Gas Seco E = 0.92
Gas de Cabeza de Pozo E = 0.77
Gas y Condensado E = 0.60
Cuando se tiene tuberas viejas, para los mismos fluidos mencionados antes se deben usar
valores menores.

7.5. Presiones de Trabajo Permisibles en Tuberas

Tericamente en un gasoducto las tuberas deben trabajar a presiones altas para tener una
mayor capacidad, pero en la prctica la presin mxima que puede soportar una tubera
depende de las propiedades mecnicas del acero y de las condiciones en que sta tenga que
trabajar. Para obtener una expresin que permita relacionar el espesor de la tubera con la
presin mxima de trabajo permisible se parte de la ecuacin general de esfuerzos de arco
que se puede deducir fcilmente y presenta la siguiente forma (1)
2ts
P
(7.42)
Do 2t
donde:
P
:
Presin mxima permisible,
t
:
Espesor de la tubera,
S
:
Lmite de proporcionalidad del acero,
Do
:
Dimetro externo de la tubera .
Con base en la ecuacin (7.42) el Instituto Americano de Normas Nacionales, ANSI por sus
iniciales en ingls, ha propuesto una serie de normas para seleccionar espesor de tuberas
dependiendo del fluido que transporta y las condiciones en que debe trabajar. Para el caso
de gasoductos la norma es la ANSI-31.8 que presenta la siguiente forma:
P *Do
t
(7.43)
2*S*F*E *T
donde:
F
: Es un factor de localizacin que depende del sitio por donde pase la lnea; es menor
de 1 y se alejar ms de este valor mientras mayor sea la posibilidad de prdidas
humanas en caso de accidente.
E
: Es un factor de eficiencia de las uniones, es menor de 1 y se aleja ms de ste valor
mientras menor sea la calidad de la soldadura y el control de calidad de sta.
T
: Es un factor de temperatura el cual es menor de 1 si la temperatura es mayor de
250F.(121.1C)

22

Otra norma propuesta por la ANSI para tuberas de gasoductos es la siguiente (2)
P *Do
t C
(7.44)
2(S*E P *T)
donde:
C:
Es una tolerancia por corrosin y por penetracin de rosca. y los dems trminos
tienen el mismo significado que en las ecuaciones (7.42) y (7.43).
En cuanto al factor de localizacin F, que aparece en la ecuacin (7.43) se debe agregar lo
siguiente(1):
La inclusin de este factor se debe a que los costos (humanos y econmicos) en el caso de
una falla de la tubera, varan dependiendo del sitio por donde esta pasa; pues unas veces la
tubera pasa por zonas despobladas donde si se presenta la falla los costos sern bajos, pero
otras veces puede pasar por zonas pobladas y en este caso una falla en la tubera ser
costosa y lamentable. Esto quiere decir que a lo largo de la tubera el valor de S
(resistencia a la cedencia) debe ser multiplicado por un factor de localizacin F, que puede
variar desde 0.72 en el caso ms favorable hasta 0.4 en los casos ms severos.
Para definir el valor de F se consideran cuatro tipos de localizaciones, de acuerdo con un
concepto conocido como ndice de una milla o ndice de las 10 millas. El ndice de una
milla se establece de la siguiente manera: Se toma una zona de una milla de longitud y
media de ancho y se cuentan los sitios donde habitan o hay posibilidad de que habiten o hay
posibilidad de que habiten personas, este nmero es el ndice de una milla. El ndice de 10
millas es el promedio de los ndices de 1 milla a lo largo de una longitud de 10 millas. Con
base en los valores del ndice de una milla y de 10 millas, se establecen las siguientes clases
de localizaciones:
Clase 1: Incluye baldios, zonas desrticas o montaosas, o planas en las que el ndice
de 10 millas es menor de 12 y el de 1 milla para cualquier seccin es menor de 20; para
este tipo de localizacin se toma una valor de 0.72.
Clase 2: Incluye zonas limtrofes alrededor de ciudades o pueblos, o reas cultivables o
industriales donde el ndice de 1 milla excede 20 o el de 10 millas excede 12. En este
caso se toma una factor de localizacin de 0.6.
Clase 3: Zonas subdivididas para propsitos residenciales o comerciales donde al
menos el 10% de los lotes estn sobre el derecho de va en la cual se va a tender la
tubera. Tambin incluye zonas en las que se encuentran edificaciones residenciales o
comerciales de tres pisos o menos. Aqu se toma F = 0.5.
Clase 4: Son localizaciones donde prevalecen edificios de ms de 4 pisos, hay trfico
pesado y denso y existen instalaciones en el subsuelo. Se le asigna un factor de 0.4.

7.6. Velocidad Permisible en Tuberas (2,11)

La velocidad de un fluido a travs de una tubera es algo que se debe analizar pues si sta es
demasiado alta se presentan problemas de desgaste de las paredes corrosin por friccin la
cual puede destruir la pelcula de inhibidor que en algunos casos protege la pared de la
23

tubera dejndola expuesta a la corrosin o si la friccin es muy alta se puede presentar

desgaste de la tubera por abrasin; por otra parte si la velocidad es demasiado baja se
pueden presentar problemas de depositacin de slidos y esto reduce el tamao de la
tubera.
Experimentalmente se ha encontrado que la mxima velocidad permisible de un gas en una
tubera para que no haya erosin se puede calcular de (2)
C
v c 0,5
(7.45)

donde:
vc
:
Velocidad erosional,

:
Densidad del fluido
C
:
Constante cuyo valor est entre 75 y 150 , y 366.3 y 732.6 cuando se usan
unidades absolutas de los sistemas ingls y SI de unidades respectivamente ; normalmente
se toma 100 y 488.
Recordando que v = q/A y la expresin para calcular densidad de un gas, la ecuacin (7.45)
se convierte en :
0, 5

(7.46)
(q e ) sc Cd
ZT
g

donde:
(qe)sc :
Es la tasa mxima permisible para evitar erosin de la tubera medida a
condiciones normales.
d
:
Dimetro de la tubera
P
:
Presin en la lnea
:
Gravedad especfica del gas
g
Z
:
Factor de compresibilidad a P y T
T
:
Temperatura de flujo
2

C es una constante que depende de las unidades usadas para las variables. Cuando se usan
unidades absolutas su valor es 24.82 y 17.72 para los sistemas ingles y SI de unidades
respectivamente ; y cuando se usan las unidades del grupo 4 para los mismos sistemas
ingls y SI que aparecen en la tabla 2 su valor es 1012.435 y 48.4 respectivamente.
Como en una tubera la presin vara desde P1 hasta P2 y el volumen de gas aumenta al
diminuir la presin , para aplicar la ecuacin (7.46) se debe usar la presin mnima a la que
se encontrar el gas en la tubera.

24

7.7. Determinacin del Dimetro de Tubera Requerido

Cuando se quiere transportar una cantidad dada de gas a travs de una distancia dada y con
una determinada caida de presin, el dimetro de la tubera requerida se calcula de la
siguiente manera:
Haciendo uso de las ecuaciones de flujo (Weymouth, Pan-Handle u otra) se determina
d, dimetro interno de la tubera.
Usando la ecuacin (7.46) se despeja d, el cual es el dimetro requerido para que no
haya velocidad erosional. El mayor de los dos dimetros calculados se selecciona como
el dimetro interno requerido.
Usando las ecuaciones (7.43) o (7.44) se encuentra el espesor de la tubera, recordando
que Do = DI + 2d y de sta forma se puede seleccionar el tamao comercial de tubera
requerida.

7.8. Denominacin de Tuberas (3)

Para identificar adecuadamente la tubera que se requiere para trabajar en unas condiciones
dadas se deben especificar una serie de caractersticas relacionadas con tamao, espesor
tipo de acero y condiciones de trabajo. En general los parmetros a especificar son: (Ver
tabla 33, Norma ANSI-B31.3).
Tamao Nominal. Generalmente es un nmero entero, el ms prximo al dimetro
externo de la tubera cuando este es menor de 14 pulgadas (355mm.) o el mismo
dimetro externo cuando este es mayor de 14 pulgadas (355 m.m.).
Dimetro Externo. En el dimetro interno seleccionado para la tubera ms el espesor
requerido.
Espesor. Es el valor de t, obtenido por ejemplo de las ecuaciones (7.43) o (7.44).
Peso. En el peso de la tubera en lbs/pie. (Kgs/m).
Clase de Peso. Se habla de peso estndar (STD), extra - estndar (XS) y ultra - estndar
(XXS). Normalmente el peso STD, es el menor y el XXS es el mximo peso para una
tubera de un dimetro dado.
Cdigo de la Tubera. Para un tamao dado, es un nmero de dos o tres cifras,
terminando siempre en cero, con el cual se identifica el tipo de tubera.
Para seleccionar una tubera requerida, conociendo el tipo de acero, el dimetro interno, la
presin y la temperatura de trabajo, se va a la norma AINSI apropiada con esta informacin
y all se obtiene el dimetro nominal, el espesor, el peso, la clase de peso y el cdigo de la
tubera.

25

Tabla 33. Norma ANSI 31.3

26

Tabla 33. (Cont.)

27

Ejemplo 7.2 (Tomado de referencia (8))

Se tiene una tubera de 12.09 pulgadas de dimetro interno y 1000millas de longitudcon
una rugosidad absoluta de 0.0006 pulgadas para transportar un gas de g = 0.60 pulgadas
con P1 =400 Lpca. y P2 = 200 Lpca. y una temperatura de flujo promedia de 520R. Las
condiciones base son: Tb = 520R y Pb = 14,7 Lpca.
Se desea calcular la capacidad de la tubera en PC/hr.
Solucin
1. Utilizando la ecuacin (7.7)
a. Clculo de Propiedades del Gas , para g = 0,6
pPc 672 Lpca.
pTc 358R
P 400 200 / 2 300 Lpca.
pPr 300 / 672 0, 446
pTr 520 / 358 1, 453
Z 0,950
1 a 60F 0, 0103

a 300 Lpca y60F 1, 05

1
a 300 Lpca y60F 1, 05*0, 0103 0, 01082 cp.

Nmero de reynolds usando la ecuacin (7.13)

g Pb
qh
0, 6*14, 7 *q h
N Re 711,8
711,8
24*1.000 D Tb
520* 24*1.000*0, 01082*12, 09

2, 2016q h
donde qh es la tasa de flujo en PCN/hr
e 0, 0006

0, 00005
D 12, 09
b. Clculo de qh por ensayo y error:
Primer ensayo:
q h 100.000 PCN / hr.

N Re 2, 2*105
f 0, 0158

28

Usando la ecuacin (7.7) y el grupo 2 de unidades del sistema ingls de la tabla 25

se tiene
T 1
q h 3, 23 b
Pb f

0,5

P12 P22 D5

g ZTL

520 1
3, 23
14, 7 0, 0158

0,5

0,5

160.000 40.000 (12, 09)5

0, 6*0,95*520*100

0,5

929.560 PCN / hr.

Segundo ensayo:
q h 500.000 PCN / hr.

N Re 1,1*106
f 0, 012
q h 1'045.083 PCN / hr
Tercer ensayo:
q h 1'000.000 PCN / hr

N Re 2, 2*106
f 0, 012
q h 1'066.633PCN / hr.
2. Usando la ecuacin de Weymouth se tendra para el grupo de unidades de la tabla 26

7.9. Correcciones por Cambio de Altura (2), (5)

Cuando en una lnea de gasoducto se presenta territorio ondulado se deben hacer
correcciones por cambio de altura; para ello solo se requiere tener en cuenta los cambios de
altura y calcular el peso de una columna de gas.
El peso de una columna de gas es:
g *dh dP
PM *g
dh dP
ZRT
y suponiendo que h aumenta y la presin disminuye hacia arriba,
29 g *g
dh ZT *dP / P
R
g *29 g h
P2
dh ZT dP / P

0
P1
R
g *29 g h
P1
ZT dP / P ZT ln P1 / P2
P2
R

29

29 g h *g
P1
EXP
P2
ZRT
P1
,
P2
g * 29 g h
EXP

ZRT
donde h = h2 - h1; llamando
29 g h
gh
Sg
C
(7.47)
ZRT
ZT
donde C = 0,0342 y 0,01875 para los sistemas mtrico e ingls respectivamente, y como
1
P P1 P2 P1 1 S , entonces
e
e S 1
P P1 S
(7.48)
e
P1 en el punto donde se inicia el ascenso o el descenso y P2 el punto donde termina. Como
h es positivo hacia arriba, ser negativo hacia abajo y esto afectara el valor de S que puede
ser negativo o positivo; h es la diferencia entre h2 y h1.

La correccin total sera entonces

eSi 1
PT h P1 Si
e

(7.49)

En algunos casos es necesario corregir la capacidad de la tubera por cambio de altura, ya

que debido a estos cambios la tubera puede transportar mayor o menor cantidad de gas.
Una forma de calcular la verdadera capacidad de una tubera es desarrollando una ecuacin
para flujo inclinado partiendo a la ecuacin de balance de energa (ecuacin 7.1).
1
H v 2 h q w
2
y suponiendo que v2 0, w = 0 y reemplazando H por Tds + VdP, y Tds por q + dLw se
tiene:
fv 2 dL
VdP gdh
0
(7.50)
2d
Una relacin entre dh y dL se obtiene de la siguiente manera:
Supongamos que h y L tienen la siguiente orientacin:

30

donde h es la direccin vertical y L la direccin de flujo. De acuerdo con las orientaciones

propuestas, se puede escribir
h dh
, y por tanto
L dL
L
dL dh
h
y llevando la anterior expresin para dL a la ecuacin (7.50) queda

fv 2
0
VdP dh g
2dh / L
Ahora, en la ecuacin anterior se reemplaza v por la ecuacin
P T Z 4
u qb b
P Tb Zb D2
2

Pb T
f
4
VdP dh g
Z
qh
0

2d h / L P T b d 2

2
2
2

f
4 Pb 2 ZT

VdP dh g 5
qh
0

2d h / L Tb
P

y separando variables se tendra:

VdP
dh
2
2
f L 4 Pb ZT 2
g 5

qh
2d h Tb P
P1
h2
VdP

dh
2
P2

2
h1
f L 4 Pb ZT 2
g 5

qb
2d h Tb P
P1
ZT / PdP
Mg h 2

dh
2
P2
2
R h1
f L 4* Pb ZT 2
1

qb
2gd 5 h Tb P
C g h
donde C tiene los mismos valores que se plantean en la ecuacin (7.47).

31

(7.51)

La ecuacin (7.51) permite calcular la presin P1 en el fondo de una columna de gas cuando
la columna es esttica o cuando se presenta flujo. Cuando es esttica qb = 0 y la ecuacin
se convierte en:
P1 ZT
P2 P dP 0, 01875 g h
(7.52)
Cuando qb es diferente de cero se aplica la ecuacin (7.51).
Para aplicar la ecuacin (7.51) se requiere resolver la integral, la cual no es directa, por
tanto dependiendo del mtodo aplicado se tendrn diferentes soluciones(2). El mtodo ms
sencillo es el conocido como de Z y T constantes; en este caso la ecuacin (48) se
transforma en:
P1
dP / P
Mgh
P2 f *L / h 4 P 2 ZT 2 RZT
2
b
1
qb

2gd 5 Tb P
P1
PdP
Mgh

2
P2
f *L / h 4 Pb * ZT 2 RZT
(7.53)
P2

qb
2gd5 Tb
El integral de la ecuacin (7.53) es de la forma:
P1 PdP
P2 P 2 C 2
donde,
0,5

f *L / h 4 Pb ZT
2
C
q b
5

2gd Tb
2

el cual se puede resolver directamente como

1
ln P 2 C 2
2
por tanto la ecuacin (7.53) se puede llevar a
1 P12 C 2 Mg
h / ZT
ln

2 P22 C 2 R
P12 C2
Mg

EXP 2
h / ZT
2
2
P2 C
R

Mg h
y llamando
S se tiene:
R ZT

P12 C 2 2S
e
P22 C 2
Si se despeja P1 se tiene:

(7.54)

32

P12 P22e2S C2 e2S 1

(7.55)

La ecuacin (52) permite calcular la presin en el fondo de un pozo de gas.

Cuando el pozo est cerrado C2=0 y la ecuacin (7.55) se convierte en
P1 P2 eS

(7.56)

La Presin obtenida con la ecuacin (7.56) se conoce como presin de fondo esttica.
La presin obtenida con la ecuacin(7.55) se conoce como presin de fondo fluyente.
Para resolver la ecuacin (7.55) o (7.56) se debe recurrir a un proceso de ensayo y error
pues S y C involucran Z y para conocerlo se requiere la presin promedia.
Un procedimiento para calcular la presin en el fondo de un pozo de gas puede ser el
siguiente:
Suponer P1 P2 0, 25 P2 /100 h /100

Calcular P P1 P2 / 2

Calcular Z
Calcular S y C, este ltimo si se est calculando presin de fondo fluyente.
Calcular P1 con ecuaciones (7.55) o (7.56) segn el caso.
Comparar P1 supuesto con P1 calculado. Si no coinciden tomar el valor calculado como
el supuesto y repetir procedimiento.

El valor de P1 obtenido por el procedimiento anterior se puede mejorar si el proceso se hace

por tramos, o sea se divide la profundidad del pozo en un determinado nmero de tramos de
igual longitud y a cada tramo, empezando por el primero al cual se le conoce su presin en
la cabeza se le aplica el procedimiento anterior para encontrarle su respectiva P 1, la cual es
tambin la P2 del tramo siguiente al cual se le pasa entonces a determinar su P1 aplicando el
procedimiento descrito. As se continua hasta llegar al fondo del pozo.
Si se despeja ahora C2 de la ecuacin (7.54) se tiene:
C2 1 e2S P2 2e2S P12
C2 e2S 1 P12 P2e 2S

C
2

P12 P2 e 2S

2S

P12 P2 e 2S

e 2S 1

En la ecuacin anterior P1 es la presin mayor y P2 es la presin menor.

Reemplazando ahora C2 por su expresin equivalente se tiene:

33

4 Pb

P 2 P e2S
f h5
ZT q 2b 1 2S 2
2gd Tb
e 1

y despejando qb se tiene:
2g *d5 Tb
q

f L 4Pb
h
2
b

2g *d5 *h * Tb

4Pb

L*f

P12 P2e2S
2
2S
ZT 1 e
1

P12 P2e2S
2
2S
ZT 1 e
1

2SR Tb P12 P22 eS

2
S
4Pb L 2Mg ZT 1 e
2

2gd5

R
qb
29

0,5

Tb 2S P12 P22eS 5 1
d

4 Pb L g 1 eS
ZTf

0,5

(7.57)

La ecuacin (7.57) es similar a una de las formas de la ecuacin para flujo horizontal, si se
escribe de la siguiente forma
0,5
d5
R Tb
qb
29 4 Pb Le g ZTf
L e2s 1
donde Le
2S

T
d5
C b
Pb Le g ZTf

0,5

2
1

P22e2s

0,5

0,5
(7.58)
P12 P22e2s

donde C es una constante que depende de las unidades; cuando se usan absolutas su valor
es 13.3 y 32.55 para los sistemas S.I e Ingls, respectivamente.
Ahora si se aplica Weymouth a la ecuacin (7.58) y se usan las unidades del grupo 1 del
sistema ingls mostrado en la tabla 26 se tiene:
0,5

16
0,5
Tb d 3
(7.59)
P12 Pe2e 2S
q b 17, 44
Pb Le g ZT

y si se usan las unidades del grupo 2 del sistema ingls de la misma tabla la ecuacin (7.59)
se transforma en:
16
Tb d 3

q b 4,1856*10
Pb Le g ZT

0,5

2
1

P22e2s

0,5

34

(7.60)

Cuando el flujo es horizontal h=0 y por tanto S=0 y Le debe ser igual a L, tal como se
puede ver retomando la expresin para Le
L e2s 1
Le
2S
Al hacer S=0 la expresin anterior se vuelve indeterminada pero tomando el lmite cuando
S tiende a cero y aplicando la regla de Lhopital se tiene:
e2S 1
d / dS e2S 1
2e2S

L*lim
lim Le L*lim
L*lim
L
S0
S0
S0
S0
2S
d / dS 2S
2

35

Bibliografia
1. Arnold K. and Stewart M. Surface Production Operations (Vol. 2, Design of Gas Handling Systems and Facilities), Gulf Publishing Company, Houston - TX - USA,
1993.
2. Kumar S. Gas Production Engineering, Chap. 7, Gulf Publishing Company, Houston
- TX - USA, 1987
3. Arnold K. and Stewart, M. Surface Productions Operations (Vol. 1, Design of Oil
Handling Systems and Facilities), Gulf Publishing Company, Houston - TX - USA,
1986.
4. Gas Proccessors and Supliers Association (GPSA), Engineering Data Book(Vols. I
and II), Tulsa Okla - USA, Tenth Edition - 1987.
5. Campbell, J.M. Gas Conditioning and Processing, John M. Campbell and Co.
International Petroleum Institute Ltd., Norman - Okla - USA, 1968, Second Edition.
6. McCain Jr. W.D. Reservoir Fluids Properties Correlations - State of Art, SPERE
05/91, pag. 266.
7. Fourni, F.J.C. and Agostini J.P. Permeation Membranes Can Efficiently Replace
Conventional Gas Treatment Processes, JPT 06/87, pag. 707 - 712.
8. Ikoku, CH. W.Natural Gas Production Engineering, John Wiley and Sons, New York
N:Y:, USA, 1984
9. Beggs, H.D Gas Productions Operations, Oil and Gas Consultants International Inc.
Publications, Tulsa-Okla., USA, 1984
10. AGA, Gas Engineering and Operations Practices-Volume III, Distribution (Book D1, System Design), Arlington Va., USA, 1990.
11. Naranjo, A. A Manejo de Produccin en Campos de Gas, Universidad Nacional de
Colombia- Sede Medelln, Carrera Ingeniera de Petrleos, 1993.
12. Naranjo, A.A. Flujo de Gas bajo Condiciones Estables, Universidad Nacional de
Colombia Sede Medelln, Carrera Ingeniera de Petrleos, 1998.

36