ESCUELA POLITÉCNICA NACIONAL

CARRERA DE INGENIERÍA EN PETRÓLEOS

INGENIERÍA DE GAS NATURAL

Dr. Johnny Zambrano

2. PROPIEDADES DE LOS GASES

La habilidad para calcular la afluencia de un sistema de producción de gas, incluso el reservorio y el

sistema de tuberías, requiere del conocimiento de muchas propiedades de gas a varias presiones y
temperaturas. Si el gas natural está en contacto con líquidos, como condensado o agua, el efecto de los
líquidos en las propiedades de los gases debe ser evaluado.

Este capítulo presenta los mejores y más amplios métodos usados para realizar los cálculos necesarios.
Alguna de la información presentados en este capítulo sólo se usarán en cálculos de reservorios y otros
solamente se usarán en el diseño de sistemas de tuberías en el capítulo siguiente.

Se realizarán problemas ejemplos y se presentarán gráficos para las correlaciones empíricas. La

aplicación de algunos de los métodos requiere de una computadora, y las correlaciones son incluidas en el
texto dado el caso.

GASES IDEALES

El entendimiento de la conducta de los gases con respecto a cambios de presión y temperatura se hacen
justificando la primera consideración de la conducta de los gases a condiciones cercanas a las condiciones
estándares de presión y temperatura; es decir:

p = 14.7 psia = 101.325 kPa (SPE usa 100 kPa)

T = 60 ºF = 520 ºR = 288.72 ºK (SPE usa 288 ºK)

Se dice que el gas se comporta como ideal a estas condiciones, y la mayoría de los trabajos iniciales con
gases se realizaron a estas condiciones. Un gas se define como ideal bajo las siguientes consideraciones:
(1) el volumen ocupado por las moléculas es pequeño comparado al volumen de gas total; (2) todas las
colisiones moleculares son elásticas; y (3) no existen fuerzas atractivas o repulsivas entre las moléculas.

La base para describir la conducta del gas ideal viene de la combinación de algunas de las así llamadas
leyes de los gases propuestas por experimentos iniciales.

Leyes Iniciales de los Gases

La Ley de Boyle. Boyle observó experimentalmente que el volumen de un gas ideal es inversamente
proporcional a la presión para un peso o masa de gas dado cuando la temperatura es constante. Esto
puede expresarse como:

1
Vα o pV =constante
p
Ley de Charles. Trabajando con gases a bajas presiones, Charles observó que el volumen ocupado por
una masa fija de gas es directamente proporcional a su temperatura absoluta, o

V
VαT o = constante
T
La Ley de Avogadro. La ley de Avogadro establece que bajo las mismas condiciones de temperatura y
presiona, volúmenes iguales de todo los gases ideales contienen el mismo número de moléculas. Esto es
equivalente a la declaración de que a una temperatura dada y presión un peso molecular de cualquier gas
ideal ocupa el mismo volumen que el peso molecular de otro gas ideal. Se ha mostrado que hay 2.73 x
1026 moléculas/lb-mol de gas ideal y que un peso molecular en libras de cualquier gas ideal a 60 ºF y 14.7
psia ocupa un volumen de 379.4 ft3. Una mol de un material es la cantidad de ese material cuya masa, en
el sistema de unidades seleccionado es numéricamente igual al peso molecular. Esto significa ese una
mol de cualquier gas ideal, es decir, 2.73 x 10 26 moléculas de cualquier gas, ocupará el mismo volumen a
una presión y temperatura dada.
Ley de Gas Ideal

Pueden combinarse las tres leyes de gases descritas previamente para expresar una relación entre la
presión, volumen, y temperatura; llamada la Ley de Gas Ideal.

Para combinar la Ley de Charles y la Ley de Boyle para describir la conducta de un gas ideal cuando se
cambian temperatura y presión, se asume una masa dada de gas cuyo volumen es V1 a presión p1 y
temperatura T1, e imaginan el proceso siguiente a través del cual el gas alcanza el volumen V2 a presión p2
y temperatura T2:

( V 1 a p 1 ,T 1 ) T =constante
⃗paso1 ( V a p 2 , T 1 )
1

( V a p2 ,T 1 ) P =constante
⃗paso2 ( V 2 a p2 , T 2 )
2

En el primer paso la presión se cambia de un valor de p1 a un valor de p2 mientras se mantiene la

temperatura constante. Esto cambia el volumen de V1 a V. En el paso 2, la presión se mantiene constante a
un valor de p2, y la temperatura se cambia de un valor de T1 a un valor de T2.

El cambio en el volumen del gas durante el primer paso puede describirse a través del uso de la Ley de
Boyle donde la cantidad de gas y temperatura se mantienen constantes. Así

p1 V 1
p1 V 1 =p 2 V o V = (2-1)
p2

donde V representa el volumen a presión p2 y temperatura T1. La Ley de Charles aplica subsecuentemente
al cambio en el volumen del gas durante el segundo paso donde la presión y la cantidad de gas se
mantienen constantes. Por consiguiente

V V2 V 2T1
= o V= (2-2)
T 1 T2 T2

La eliminación del volumen, V, entre las Ecuaciones 2-1 y 2-2 dan

p1 V 1 V 2 T 1
=
p2 T2

p1 V 1 V 2 p 2
= (2-3)
T1 T2

Así, para una cantidad dada de gas, pV/T = constante. La constante se designa con el símbolo R cuando
la cantidad de gas es igual a un peso molecular. Es decir,
 pV M
=R
T

donde VM representa el volumen de un peso molecular del gas a p y T.

Para mostrar que R es el mismo para cualquier gas, se invoca la Ley de Avogadro. En forma simbólica, la
Ley de estado
V MA =V MB

donde VMA representa el volumen de un peso molecular de un gas A y VMB representa el volumen de un
peso molecular de un gas B, ambos a presión p, y temperatura T. Esto implica que

pV MA R T
=R A o V MA = A
T p
y

pV MB R T
=R B o V MB = B
T p
donde RA representa el gas constante para un gas A y RB representa el gas constante para un gas B. La
combinación de las ecuaciones anteriores revela que

R AT RB T
= o R A=R B
p p
Así, la constante R es la misma para todo los gases ideales y se llama Constante Universal de los Gases.
Por consiguiente, la ecuación de estado para uno peso molecular de cualquier gas ideal es

pV M =RT ( 2−4 )

Para n moles de gas ideal esta ecuación se transforma en

pV =nRT (2-5)
donde V es el volumen total de n moles de gas en la temperatura, T, y presión, p. Donde n es la masa de
gas dividida para el peso molecular, la ecuación puede escribirse como

m
pV = RT
M
o, donde m/V es la densidad de gas,

pM
ρ= (2-6)
RT
Esta expresión es conocida por varios nombres como la ley de gas ideal, la ley de gas general, o la ley de
gas perfecto. Esta ecuación tiene valor práctico limitado por cuanto no se conoce la conducta del gas
como un gas ideal; sin embargo, la ecuación describe la conducta de la mayoría de los gases reales a
presión baja y dan bases para desarrollar las ecuaciones de estado que más adecuadamente describen la
conducta de gases reales a presiones elevadas.
El valor numérico de la constante R depende de las unidades usadas para expresar la temperatura, presión
y volumen.

Ejemplo

Sea la presión expresada en pisa, el volumen en pies cúbicos, la temperatura en grados Rankin y moles en
libra-mol. La Ley de Avogadro determina que 1 lb-mol de cualquier gas ideal ocupa 379.4 ft 3 a 60 ºF y
14.7 psia. Por consiguiente,
pV ( 14 . 7 psia ) ( 379. 4 cu f t )
R= = =10. 73 psia cu ft / lb-mole ºR .
nT ( 1 lb-mole ) (520 ºR )

En la Tabla 2-1 se expresa los valores numéricos de R por el varios sistema de unidades.

Ejemplo 2-1

Calcular la masa de gas metano contenido a 1000 psia y 68ºF en un cilindro con volumen de 3.2 ft 3.
Asumir que el metano es un gas ideal.

Solución

Ejemplo 2-2

Calcular la densidad del metano a condiciones estándares.

Solución
Mezclas de gas ideales

El tratamiento anterior de la conducta de los gases sólo se aplica para singularizar los componentes del
gas. Como el ingeniero de gas raramente trabaja con gases puros, la conducta de una mezcla
multicomponente de gases se debe de tratar. Esto requiere la introducción de dos leyes de gas ideales
adicionales.

La Ley de Dalton. La Ley de Dalton establece que cada gas en una mezcla de gases ejerce una presión
igual a la que ejercería si ocupa el mismo volumen que la mezcla total. Esta presión se llama la presión
parcial. La presión total es la suma de las presiones parciales. Esta ley sólo es válida cuando la mezcla y
cada componente de la mezcla obedecen a la ley de gas ideal. Algunas veces se la denomina como la Ley
de presiones aditivas.

La presión parcial ejercida por cada componente de la mezcla de gas puede ser calculada usando la ley de
gas ideal. Considerar una mezcla que contiene nA moles del componente A, nB moles del componente B y
nC moles del componente C. La presión parcial ejercida por cada componente de la mezcla de gas puede
determinarse con la ecuación de gas ideal:

RT RT RT
p A=n A , p B=nB , pC =nC .
V V V
Según la Ley de Dalton, la presión total es la suma de las presiones parciales.

p= p A + p B +p C

RT RT RT
p=n A +n B + nC
V V V

RT
p=
V
∑ n j =VRT n .
La relación de la presión parcial del componente j, pj, a la presión total de la mezcla p es

pj nj nj
= = =y j, (2-7 )
p ∑nj n

donde yj es definida como la fracción molar del componente jth en la mezcla de gas. Por consiguiente, la
presión parcial de un componente de una mezcla de gas es el producto de su fracción molar por la presión
total.

La Ley de Amagat. La Ley de Amagat establece que el volumen total de una mezcla gaseosa es la suma
de los volúmenes que cada componente ocuparía a una presión y temperatura determinada. Los
volúmenes ocupados por los componentes individuales son conocidos como volúmenes parciales. Esta
ley sólo es correcta si la mezcla y cada uno de los componentes obedece la ley de gas ideal.

El volumen parcial ocupado por cada componente de una mezcla de gas que consiste en nA moles del
componente A, nB moles del componente B, y así sucesivamente, puede calcularse usando la ley de gas
ideal.

RT RT RT
V A=n A , V B=nB , V C =n C ...
p p p
Así, según Amagat, el volumen total es
V =V A +V B +V C

RT RT RT
V =n A +n B +nC
p p p

RT
V=
p
∑ n j = RT
p
n.

La relación del volumen parcial de componente j al volumen total de la mezcla es

RT
V j nj p nj
= = =y j. (2-8 )
V RT n
n
p
Esto implica que para un gas ideal la fracción del volumen es igual a la fracción molar.

Peso Molecular Aparente

Como una mezcla de gas está compuesta de moléculas de varios tamaños, no es estrictamente correcto
decir que una mezcla de gas tiene un peso molecular. Sin embargo, una mezcla de gas se comporta como
si fuera un gas puro con un peso molecular definitivo. Este peso molecular es conocido como un peso
molecular aparente y se define como

M a =∑ y j M j (2-9)

Ejemplo 2-3

El aire seco es una mezcla de gas consistente esencialmente de nitrógeno, oxígeno y pequeñas cantidades
de otros gases. Calcular el peso molecular aparente del aire dada su composición aproximada. El peso
molecular de cada componente puede ser encontrada en la Tabla 2-2.

Solución

Usualmente se considera un valor de 29.0 de exactitud adecuada para los cálculos de ingeniería.

Se define la gravedad específica de un gas como la relación de la densidad del gas a la densidad de aire
seco tomado a las condiciones estándares de temperatura y presión. Simbólicamente,
 ρg
γ g= (2-10 )
ρair
Asumiendo que la conducta del gas y aire pueden ser representada por la ley de gas ideal, la gravedad
específica puede darse como

Mg p
RT M M
γ g= = g = g,
M air p M air 29
RT

donde Mair es el peso molecular aparente del aire. Si el gas es una mezcla, esta ecuación se vuelve

Ma Ma
γ g= = ( 2-11)
M air 29

donde Ma es el peso molecular aparente de la mezcla de gas.

Ejemplo 2-4

Calcular la gravedad de un gas natural de la composición siguiente:

Solución
GASES REALES

Varias asunciones se han realizado para formular la ecuación de estado para los gases ideales. Puesto que
estas asunciones no son correctas para gases a las presiones y temperaturas que se desvían de ideal o de
las condiciones normales, se deben hacer correcciones para considerar la desviación de la conducta ideal.
El método para la corrección más ampliamente usado en la industria del petróleo es el factor de
compresibilidad del gas, comúnmente llamado el factor Z. Se define como una relación del volumen real
ocupado por una masa de gas a alguna presión y temperatura para el volumen del gas que ocuparía si se
comportara idealmente. Es decir,

V actual
Z= o V actual=ZV ideal
V ideal

La ecuación de estado es

V actual
pV ideal =nRT o p =nRT
Z

Por consiguiente, la ecuación de estado de cualquier gas se vuelve

pV =ZnRT (2-12)
donde, para un gas ideal, Z = 1.

El factor de compresibilidad varía con cambios en composición de gas, temperatura y presión. Puede ser
determinada experimentalmente. El resultado de determinaciones experimentales de factores de
compresibilidad normalmente se da gráficamente y usualmente toma la forma mostrada en figura 2-1. La
forma de la curva es consistente con el presente conocimiento de la conducta de gases; a presión muy
baja las moléculas están relativamente alejadas y las condiciones de conducta de gas ideal probablemente
están más a ser reunidas. A presión baja el factor de compresibilidad tiende a un valor de 1.0, lo que
indica que la conducta de gas de ideal ocurre de hecho. Los factores de compresibilidad para varios gases
del hidrocarburo se dan en las figuras 2-2, 2-3 y 2-4.
Ejemplo 2-5

Calcular la masa del gas metano contenido a 1000 psia y 68ºF en un cilindro con volumen de 3.2 ft 3.

Solución

Si se asume conducta ideal, la masa calculada podría haber sido:

m = 10.2 * 0.89 = 9.08 lbm.

Mezclas de gas reales

Las cartas del factor de compresibilidad están disponibles para la mayoría de los componentes simples de
gases de hidrocarburos ligeros, pero en práctica un componente simple de gas se encuentra raramente.
Para conseguir factores Z para las mezclas de gases naturales, se usa la ley de estados correspondientes.
Esta ley establece que la relación del valor de cualquier propiedad intensiva al valor de esa propiedad en
el estado crítico es relacionada a las relaciones de la temperatura y presión absoluta prevaleciente a la
temperatura y presión crítica para la misma función para todas las substancias similares. Esto significa
que todo los gases puros tienen el mismo factor Z a los mismos valores de presión y temperatura reducida
donde los valores reducidos se definen como
 T p
T r= , pr=
Tc pc

donde Tc y pc son la temperatura crítica y presionan para el gas, respectivamente. Los valores deben estar
en unidades absolutas.

Ejemplo 2-6

Calcular la densidad de etano a 900 psia y 110ºF. Las propiedades críticas, de la Tabla 2-2, son Tc = 550
ºR, pc = 708 psia, M = 30.1 lbm/lb-mol.

Solución

De Fig. 2-5, Z = 0.34, entonces:

Asumiendo conducta ideal, el valor calculado sería:

 = 13.03 * 0.34 = 4.43 lbm/ft3.

Se ha mostrado que la Ley de los Estados Correspondientes trabaja mejor para los gases de características
moleculares similares. Esto es ventajoso ya que la mayoría de los gases con los que trata el ingeniero de
petróleos están compuestos de moléculas de la misma clase de compuestos orgánicos conocidas como
hidrocarburos de parafina. La Ley de estados correspondientes se ha extendida para cubrir mezclas de
gases que están estrechamente relacionados químicamente. Puesto que es algo difícil de obtener el punto
crítico para las mezclas del multicomponentes, las cantidades de temperatura pseudocrítica y presión
pseudocrítica han sido concebidas. Estas cantidades se definen como

T pc=∑ y j T cj y p pc =∑ y j p cj (2-13)

Estas cantidades pseudocríticas se usan para las mezclas de gases en exactamente de la misma manera
como se usan las temperaturas críticas reales y las presiones críticas para gases puros. Debe entenderse,
sin embargo, que estas propiedades pseudocríticas simplemente se desarrollaron para usar en correlación
a los factores de compresibilidad y de ninguna manera relaciona a las propiedades críticas reales o
actuales de la mezcla de gases.

Ejemplo 2-7

Calcular la temperatura y presión pseudocrítica de la siguiente mezcla de gas natural. Usar las constantes
críticas dadas en la Tabla 2-2.

Tpc = 383ºR ppc = 667 psia

El factor de compresibilidad para gases naturales ha sido correlacionado usando propiedades

pseudocríticas y se han presentado en las figuras 2-6, 2-7, y 2-8. Los factores de compresibilidad son una
función de la composición así como la temperatura y presión. Se ha puntualizado que los componentes de
la mayoría de los gases naturales son hidrocarburos de la misma familia, y por consiguiente es posible
una correlación de este tipo.
Ejemplo 2-8

Calcular el volumen ocupado por 1 lb-mol de el gas natural dado en el ejemplo 2-7 a 120ºF y 1500 psia.

Solución

Mediante programas computacionales se puede determinar el factor Z. Uno de estos programas se

reproduce la figura 2-6.

En algunos casos la composición de un gas se dará en peso o en porcentaje de masa más que por
porcentaje molar. En este caso, la composición debe ser primero convertida a fracción molar o porcentaje
antes que las propiedades de la mezcla puedan ser calculadas. El ejemplo siguiente ilustra esta
conversión.

Ejemplo 2-9

Una mezcla de gas consiste de 50% de C 1, 30% de C2 y 20% de C3 en peso. Calcular el peso molecular
aparente y la gravedad específica de esta mezcla.

Solución

Asumir 100 lbm de gas como base.

Si la fracción de volumen se da condiciones diferentes a las estándares, la fracción de volumen debe

convertirse a una base de fracción molar, teniendo en cuenta la desviación de la conducta ideal. El
ejemplo siguiente ilustra este procedimiento.
Ejemplo 2-10

Para convertir de volumen a moles, asumir 100 ft 3 como base. Entonces usar:

La mayoría de los casos, la composición de un gas natural puede ser conocida y el peso molecular
aparente y las propiedades críticas pueden ser calculadas como previamente se describió.
Ocasionalmente, sin embargo, sólo la gravedad del gas puede ser conocida. También, es muy fácil de
medir la gravedad del gas en el campo. Si la composición es desconocida, o si los requerimientos de
exactitud no justifican cálculos más grandes, la figure 2-9 se puede usar para estimar las propiedades
pseudocríticas. Las propiedades pueden ser calculadas usando las ecuaciones siguientes

T pc=170 .5+307 .3 γ g (2-14 )

p pc =709 .6−58 .7 γ g (2-15)

Para fluidos condensados:

T pc=187 +330 γ g -71. 5 γ

g2

p pc =706−51. 7 γ g −11. 1 γ
g2
Ejemplo 2-11

Usando la correlación empírica (Fig. 2-6) para las propiedades empíricas, calcular TPC y ppc el gas en el
ejemplo 2-4. Comparar estos valores con los obtenidos usando la composición (ejemplo 2-7).

Solución

Los valores del ejemplo 2-7 fueron 383ºR y 667 psia.

FACTOR VOLUMÉTRICO DE FORMACIÓN DE GAS

En muchas operaciones que involucran producción de gas la tasa de flujo y las cantidades producidas son
medidas a las condiciones normales tales como scf/día o scf. Los cálculos de ingeniería de reservorios y
de línea de flujo requieren los volúmenes a condiciones in situ de presión y temperatura, y por
consiguiente se necesita un conveniente factor de la conversión de las condiciones estándares a las
condiciones in situ. Este factor de la conversión se llama el factor volumétrico de formación del gas y se
define como el volumen real ocupado por el gas a alguna presión y temperatura dividida por el volumen
que el gas ocuparía a las condiciones estándares. Eso es,

ZnRT
V p ,T p ZTp sc vol
B g= = =
V sc Z sc nRT sc Z sc T sc p std vol
psc

Usando Tsc = 520 ºR, psc = 14.7 psia y Zsc = 1, da

ZT (14 .7 ) ZT vol ft 3
B g= =0. 0283 , e. g . (2-16)
( 1 ) (520 ) p p std vol scf

A veces es conveniente expresar Bg en barriles de volumen, o

ZT bbls
B g =0 . 00504 (2-17)
p scf
Para las constantes dadas arriba, la presión está en psia y la temperatura en ºR. Para el SI (p = kPa, T =
ºK).

ZT (101 . 325 ) ZT
B g= =0 .351 (2-18)
( 1 ) ( 288. 72 ) p p

Ejemplo 2-12

Calcular el factor volumétrico de formación para un gas natural teniendo una gravedad de 0.7 a una
temperatura de 93ºC y una presión absoluta de 10343 kPa.

Solución
De la figura 2-6, Z = 0.88

CORRECCION POR IMPUREZAS

Los gases naturales frecuentemente contienen materiales diferentes a los hidrocarburos, como nitrógeno
(N2), dióxido del carbono (CO2) y sulfuro de hidrógeno (H2S). La presencia de estas impurezas afectan el
valor obtenido del gráfico del factor Z.

Un procedimiento por ajustar las propiedades críticas del gas fue propuesto por Wichert y Aziz en 1970.
Las propiedades críticas ajustadas entonces son usadas calculando las propiedades reducidas, y el factor
Z se obtiene de la figura 2-6.

El procedimiento por obtener el Z-factor para gases agrios es:

1. Determinar ppc y Tpc para el gas que usa la composición del gas de la figura 2-9 y
2. Calcule las propiedades críticas ajustadas,

T ' pc=T pc −ε ,

p pc T ' pc
p' pc =
T pc + ε ( B−B2 )

donde

ε=120 ( A 0 .9 − A 1. 6) +15 ( B0 . 5−B 4)

B= mole fraction H2 S ,

A= mole fraction CO2 +B

ε= correction factor, ºF
3. Calcular las propiedades reducidas usando las propiedades críticas corregidas.

T p
T r= , pr =
T ' pc p ' pc

4. Encontrar Z de figura 2-6 o usando la correlación del apéndice (subprograma).

Ejemplo 2-13

Un gas conteniendo 2.87% de CO2 y 23.27% de H2S tiene una presión crítica de 822 psia y una
temperatura crítica de 465 ºR. Encontrar el factor de compresibilidad del gas, Z, para p = 1000 psia y T =
100ºF.

Solución

De figura 2-6, Z = 0.78.

OTRAS ECUACIONES DE ESTADO

Una de las limitaciones en el uso de la ecuación de compresibilidad para describir la conducta de gases es
que el factor de compresibilidad no es constante, y por consiguiente no pueden hacerse manipulaciones
matemáticas directas pero pueden lograrse a través de técnicas gráficas o numéricas. La mayoría de las
otras ecuaciones de estado normalmente usadas se derivaron para factores de corrección, que corrigen la
ley de gas ideal para no idealmente, pueden asumirse constante, por consiguiente permiten usarse las
ecuaciones en cálculos matemáticos que involucran diferenciación o integración.

Se han propuesto muchas ecuaciones de estado para describir la conducta del gas y muchas
modificaciones y mejoras se han realizado. Se describirán sólo dos de las ecuaciones normalmente
usadas. Estas se usarán posteriormente en el cálculo de conducta de la fase.

Ecuación de Benedict-Webb-Rubin (BWR)

La ecuación presentada por Beattie y Bridgeman fue modificada y generó una ecuación con ocho
constantes empíricas.

RT ( )
Bo RT − A o −C o / T 2 ( bRT −a ) aα c γ γ
p=
VM
+
VM 2
+
VM 3
+ +
V M T2 V
6
M
1+
VM
3
(
EXP
VM 2 ) ( ) 2
( 2-19 )

Los parámetros Bo, Ao, Co, a, b, c, , y  son constantes para los compuestos puros y son funciones de
composiciones para las mezclas. Se dan las constantes para los compuestos puros en la tabla 2-3. Estas
constantes pueden combinarse para usar con mezclas de gases según las reglas de mezclas siguientes.
 2
A o= ∑ y j A
( )
oj 1/2

B o =∑ y j B oj

2
C o= ∑ y j C
( )
oj 1/2

3
a= ∑ y j a 1/3
( )
j

3
b= ∑ y j b
( 1/3 )
j

3
c= ∑ y j c
( )
j1/3

3
α= ∑ y j α
( )
j1/3

2
γ= ∑ y j γ
( 1/2 )
j

VM es el volumen molar en ft3/lb-molar.

Ecuación de Redlich-Kwong (RK)

La ecuación de Redlich-Kwong involucra sólo dos constantes empíricas en comparación a las ocho
requeridas en la ecuación de BWR. La ecuación original de RK es

RT a
p= − (2-20 )
V M −b V M ( V M +b ) T 0 .5

donde

Ca R 2 T 2. 5
c
a=
pc

Cb RT c
b=
pc
Para aplicar los cálculos con la ecuación de RK, especialmente para aplicación de mezclas, otras
constantes se han definido como

0. 5
CaT
( )
0 .5
a c 2. 5
A= 2 2. 5
R T ( ) =
pc T 2. 5

b C T
B= = b c
RT pc T

donde, para la temperatura en ºR y presión en psia, Ca = 0.42748, Cb = 0.08664.

La ecuación de RK posteriormente fue modificada por Soave para mejorar su exactitud cuando se aplicó
a las mezclas como sigue:

RT aα
p= − (2-21)
V M −b V M ( V M +b )

donde  es una función de temperatura. Usando las definiciones siguientes, un factor de compresibilidad
de gas puede ser calculado

C a α pT
c2
A=
pcT 2

C b pT c
B=
pc T

Z 3 −Z 2 +Z ( A−B−B2 ) −AB=0 (2-22 )

Las constantes A y se calculan B para las mezclas por

p 2
A m=C a
T2 (∑ ( yi T ci α i0 .5 / pci0 .5 )) (2-23)

p
B m=C b ∑ yi T ci / pci (2-24)
T
donde yi es la fracción molar del ith componente de la mezcla. La evaluación de  involucra el cálculo
del factor acéntrico, .

2
α= 1+S 1−T )) ,
( ( R 0. 5
(2-25)

S=0 . 480+1 . 574 ω−0 .176 ω 2 , (2-26)

log p c−1 .167

ω=C a T B
[ T c −T B
−1 . 0
] (2-27)
donde TB es el punto de ebullición del componente a 14.7 psia en ºR. También pueden obtenerse valores
de  de la Tabla 2-2.

Ejemplo 2-14

Usando la modificación de Soave para la ecuación RK, calcular el volumen molar de la siguiente mezcla
de gas a p = 250 psia, T = 100ºF.

Solución

Cálculos para C1:

De Tabla 2-2,  = 0.0104

COMPRESIBILIDAD DEL GAS

La compresibilidad isotérmica de un gas es la medida del cambio en volumen por unidad de volumen con
cambio de presión a temperatura constante. En forma de ecuación
 1 ∂V
C=− ( )
V ∂p T
(2-28)

No debe de confundirse con el factor de compresibilidad del gas, Z. La compresibilidad es requerida en

muchos ecuaciones de flujo de reservorios de gas y puede ser evaluada de la manera siguiente.

Compresibilidad del gas ideal

Para un gas ideal,

nRT
V=
p
, ( ∂V∂ p ) =− nRTp
T
2

y por consiguiente,

1 ∂V p nRT 1
C=−
V ∂p
=−
nRT
− 2 =
p p ( )
Compresibilidad del gas real

Para un gas real,

ZnRT
V=
p
y donde Z = f(p), debe ser incluido en la derivada

∂Z
p −Z
∂p
( ∂V∂ p ) =nRT
T p2
Sustituyendo estas expresiones en la ecuación que define compresibilidad da

− p nRT ∂ Z 1 1 ∂Z
C=
[ (
ZnRT p 2
p
∂p
−Z = −
p Z ∂p )] (2-29 )

La evaluación de C para gases reales requiere determinando cómo el factor Z varía con la presión a la
presión y temperatura de interés. Ya que la mayoría de las cartas y ecuaciones predice Z como una
función de presión y temperatura reducida, se ha definido la compresibilidad reducida como Cr = Cpc.
Esto puede ser expresado como una función de pr a un valor fijo de Tr por

1 1 ∂Z
Cr = −
( )
pr Z ∂ p r Tr

Valores de (Z/pr)Tr pueden obtenerse de la pendiente de una curva constante de Tr de la figura 2-6 al
factor Z de interés. Los valores de CrTr como una función de pr y Tr han sido presentados gráficamente
por Mattar y col.. en las figuras 2-10 y 2-11. También pueden calcularse los cambios de Z con la p usando
una expresión analítica calculando el factor Z a presión ligeramente sobre y por debajo de la presión de
interés. Es decir,
∂Z Z 1−Z 2
( ) (
∂ pr Tr
≈
p r 1− p r 2 )
Tr
Ejemplo 2-15

Calcular la compresibilidad de un gas a temperatura de 40ºF y una presión de 665 psia si Tc = 357 ºR, pc =
674 psia.

Solución

El problema también puede ser resuelto usando la aproximación de Z/pr:

VISCOSIDAD DEL GAS

La viscosidad de un fluido es una medida de la habilidad del fluido a fluir, o la relación de la fuerza del
esfuerzo (cisallamiento) a la tasa del esfuerzo (cisalla). La viscosidad normalmente se expresa en
centipoises o poises, pero puede convertirse a otras unidades para la compatibilidad de la unidad.

-2 -3 2
1 poise = 100 centipoises = 6.72 x 10 lbm/ft-seg = 2.09 x 10 lbf-seg/ft = 0.1 kg/m-seg
La viscosidad del gas es difícil medir experimentalmente y por propósitos de ingeniería puede ser
determinada con precisión de correlaciones empíricas. La correlación más ampliamente aceptada usada
en el pasado ha sido la de Carr y col. El método de Carr se presenta en las figuras 2-12 y 2-13. La
viscosidad a una atmósfera es obtenida primero de figura 2-12 y se corrige de impurezas diferentes a
hidrocarburos si es necesario. Esta viscosidad sólo es una función del peso molecular o gravedad y
temperatura. La corrección para presión se obtiene de figura 2-13 como una función de presión reducida
y temperatura.
Una expresión analítica para la viscosidad de gas de hidrocarburo fue presentada por Lee y col. en 1966.
La ecuación es

μg =K∗10−4 exp ( Xρgy ) ( 2-30 )

donde

( 9. 4 +0 . 02 M ) T 1 .5
K=
209+19 M +T

986
X =3. 5+ +0 . 01 M
T

y=2 . 4−0 . 2 X
En estas ecuaciones, T = ºR, g = cp, M = peso molecular y g = gm/cm3.

Ejemplo 2-16

Usando los métodos de Carr y Lee, calcular la viscosidad de un gas de gravedad específica 0.8 a una
presión de 2000 psia y una temperatura de 150ºF. El gas contiene 10% de H 2S y 10% de CO2.

Solución

Método de Carr

De figura 2-12, 1 = 0.0111 cp

Corrección por 10% de H2S = + 0.0003
Corrección por 10% de CO2 = + 0.0006

De figura 2-13, 1 = 1.61

Método de Lee

El método de Lee no incluye un método de corrección para impurezas que no son hidrocarburos; por lo
tanto, el valor obtenido podría ser para un gas de hidrocarburo puro. Sin embargo, si el factor Z usado en
el cálculo de la densidad del gas ha sido corregido, el método de Lee es válido para gases agrios.

A veces es conveniente expresar la viscosidad de un fluido como una viscosidad cinemática, . La

relación entre viscosidad y la viscosidad cinemática es

μ
ν= (2-31)
ρ
donde  es la densidad del fluido.
SISTEMAS GAS-AGUA

En muchos casos el diseño de ingeniería de operaciones de producción de gas involucraría gas natural en
contacto con agua. Esta agua puede ser agua connata del reservorio o agua producida de alguna otra zona.
En todo caso, es necesario determinar el agua contenida en el gas en el estado de vapor, el gas disuelto en
el agua, y bajo qué condiciones de presión y temperatura los hidratos de gas podrían ser formados.

La solubilidad de gas natural en agua

La solubilidad del gas natural en agua es muy baja a muchas presiones y temperaturas de interés en
ingeniería de producción de gas. Los factores primarios que afectan la cantidad de gas que podría ser
evolucionado desde el agua saturada con el gas depende de la presión, temperatura y contenido de sólidos
del agua. La relación se ilustra en figura 2-14. El factor de la corrección para salinidad puede calcularse
de

Rsw
=1−XY × 10 -4 (2-32 )
R swp

donde

Y = salinidad del agua, ppm

X = 3.471/T0.837
T = temperatura, ºF
Rsw = solubilidad del gas en salmuera
Rswp = solubilidad del gas en pura agua

Ejemplo 2-17

Calcular los ft3 de gas disuelto en una salmuera conteniendo 50000 ppm a una presión de 5000 psia y una
temperatura de 200ºF.
Solución

De figura 2-14, Rswp = 20 ft3/STB

La solubilidad del agua en gas natural

La cantidad de vapor de agua contenida en gas natural a varias temperaturas y presiones debe conocerse
en orden para relacionar la producción de agua de superficie a volúmenes equivalentes en el reservorio.
El contenido de agua también es importante en cálculos que involucran formación de hidratos. El
contenido de agua depende de la presión, temperatura y salinidad del agua. La figure 2-15 muestra la
masa o peso de agua fresca contenida en gas como una función de presión y temperatura. La ecuación
siguiente puede usarse para corregir la salinidad de agua.
 Ws
=1−2 . 87×10−8 Y 1. 266 (2-33)
W sp

donde

Ws = contenido de salmuera, lbm/MMscf,

Wsp = contenido de agua pura de figura 2-15,
Y = salinidad de agua, ppm.

Ejemplo 2-18

Determinar el contenido de agua en un gas natural en contacto con una salmuera de 100000 ppm a 3000
psia y 200ºF.

Solución

De figura 2-15, Wsp = 280 lbm/MMscf

Hidratos de gas

Los hidratos de gas son compuestos cristalinos, formados por la combinación química de gas natural y
agua bajo presión a temperaturas considerablemente sobre el punto de congelación del agua. En la
presencia de agua libre, los hidratos se formarán cuando la temperatura del gas está debajo de una cierta
temperatura, llamada la "temperatura del hidrato". La formación del hidrato a menudo es confundida con
la condensación, y la distinción entre los dos debe ser entendida claramente. La condensación del agua de
un gas natural bajo presión ocurre cuando la temperatura está a o por debajo del punto de rocío a esa
presión. El agua libre obtenida bajo tales condiciones es esencial para la formación de hidratos que podría
ocurrir a o por debajo de la temperatura del hidrato a la misma presión. De aquí, la temperatura del
hidrato podría ser debajo, y quizás igual que, pero nunca sobre la temperatura del punto de rocío.

Resolviendo problemas de producción de gas, es necesario definir, y por eso evita, condiciones que
promueven la formación de hidratos. Los hidratos pueden estrangular la sarta de flujo o producción,
líneas de superficie, y equipos de prueba de pozos. La formación de hidratos en la sarta de flujo podría
resultar en un valor más bajo para las presiones de cabeza de pozo medidas. En un dispositivo de
medición de tasa de flujo, la formación del hidrato podría resultar en tasas de flujo más bajas. La excesiva
formación del hidrato también puede bloquear completamente las líneas de flujo y equipo de superficie.

Condiciones que promueven la formación del hidrato son:

1. El gas a o por debajo de su punto de rocío del agua con agua "libre" presente,
2. baja temperatura, y
3. alta presión.

Una discusión extensa de formación del hidrato y prevención se da más adelante.

SISTEMAS DE CONDENSADO DE GAS

Como la exploración para gas natural se extiende a los horizontes más profundos, se descubren más
reservorios que contienen condensado de gas. El gas puede estar en fase gaseosa a condiciones de
reservorio inicial pero puede condensar para formar un poco de líquido en algún punto en el camino al
separador.
El diseño de ingeniería de estos sistemas de producción requieren de algún conocimiento de conducta de
fase. Esta sección trata del fenómeno de cambio de la fase y con métodos para calcular estos cambios de
fase.

Conducta de fase

La conducta de fase es simple para los sistemas de componentes simples pero se vuelve más complicado
según más componentes se agreguen al sistema. Una discusión del sistema más simple llevará a una
comprensión de los sistemas más complejos.

Fluido de componente simple. La conducta de fase de un fluido puede ser descrita por la determinación
de su respuesta a cambios de presión y de temperatura. En un líquido las moléculas están muy
íntimamente juntas, pero en un gas las moléculas están ampliamente separadas. Existen ciertas fuerzas
que tienden a confinar o dispersar las moléculas. Las fuerzas que confinan son principalmente la presión
y la atracción molecular. Las fuerzas de dispersión son la energía cinética y la repulsión molecular. Las
magnitudes relativas de las fuerzas de confinación y dispersión determinan si el fluido es un líquido o un
gas.

Un aumento en temperatura incrementa la energía cinética de las moléculas y por lo tanto las fuerzas de
dispersión, mientras que un aumento en la presión incrementa las fuerzas de confinamiento.

La conducta de un fluido de un solo componente puede ser ilustrada en un gráfico de presión versus
temperatura, denominado diagrama de fase. Una línea uniendo todos los puntos donde las fuerzas de
confinamiento y rechaza son iguales puede ser trazada. Esta es una línea de presión de vapor y se ilustra
en figura 2-16. En cualquier punto sobre la línea de presión de vapor el fluido está como líquido, mientras
que por debajo de la línea el fluido es un gas. En el punto crítico, Punto C, no hay ninguna distinción
entre el líquido y gas. Un cambio de fase a presión constante ocurre cuando la temperatura cambia, como
se ilustra por la línea 1-2. Un cambio de fase ocurrirá a temperatura constante según la presión se cambie
(línea 3-4).

Fluidos Multicomponentes. Cuando más de un fluido está presente, la diferencia en tamaño de la

molécula y la energía tiene una influencia en el cambio de fase. La figura 2-17 ilustra el diagrama de la
fase para una mezcla de dos componentes y la figura 2-18 muestra el diagrama de fase para una mezcla
de multicomponente de hidrocarburo típica. No hay ninguna transición exacta de un líquido a vapor o del
vapor a líquido, pero las moléculas pueden escapar desde el líquido o gas a presiones y temperaturas
diferentes debido a la atracción molecular. El sitio de todos los puntos donde la primera burbuja de gas
aparece en un líquido según las condiciones de presión y de temperatura sean cambiadas se llama la línea
de punto de burbuja. El sitio de todos los puntos donde la primera gota de líquido aparece según las
condiciones para un gas se cambia es la línea de punto de rocío. La presión más alta a la que un gas puede
existir se llama cricondenbárica mientras la temperatura más alta a la que un líquido puede existir es la
cricondentérmica.

Procesos separación. Los dos tipos de procesos de separación que pueden ocurrir según la presión o
temperatura en un fluido cambian se llaman diferencial y flash o equilibrio. Estos procesos se ilustran en
las figuras 2-19 y 2-20. En el proceso diferencial el gas es removido del contacto con el líquido cuando
este es separado, lo que significa que la composición de los líquidos restantes cambia constantemente. En
el proceso flash el gas permanece en contacto con el líquido, lo que significa que el proceso es a
composición constante. Ambos tipos de separación, o una mezcla de los dos, pueden ocurrir según el gas
es producido desde un reservorio hacia el separador.
Tipos de reservorio de gas

Los reservorios que producen gas natural pueden ser clasificados en esencialmente cuatro categorías.
Éstos son:

1. Gas seco - el fluido existe como un gas tanto en el reservorio como en el sistema de tubería. El único
líquido asociado con el gas desde un reservorio de gas seco es el agua. Un diagrama de fase de un
reservorio de gas seco se ilustra en la figura 2-21;

2. Gas húmedo - el fluido inicialmente existe como un gas en el reservorio y remanente en la fase gaseosa
según decline la presión a temperatura de reservorio. Sin embargo, produciéndose a la superficie, la
temperatura también decrece, causando condensación en el sistema de tubería y separador. Esta conducta
es ilustra en figura 2-22;

3. Condensado retrógrado de gas - el fluido existe como un gas a las condiciones iniciales del reservorio.
Según la presión del reservorio disminuya a temperatura de reservorio, la línea de punto de rocío es
cruzada y forma líquido en el reservorio. También se forma líquido en el sistema de tubería y separador.
La figura 2-23 describe la conducta de un sistema de gas retrógrado;

4. Gas asociado - muchos reservorios de petróleo existen a la presión de punto de burbuja del sistema de
fluido a condiciones iniciales. Puede producirse gas libre desde la capa de gas de semejante sistema. El
gas que inicialmente está disuelto en el petróleo también puede ser producido como gas libre en la
superficie. El diagrama de fase de semejante sistema dependerá de las propiedades del petróleo asociadas
con el gas. Un diagrama de fase ejemplo se muestra en figura 2-24.

El hecho que los hidrocarburos líquidos frecuentemente estén en contacto con el gas natural hace
imperativo que los métodos estén disponibles para calcular los volúmenes o masas de cada fase y también
la composición de cada fase existente a las varias condiciones de presión y temperatura.
Cálculos de la separación flash o de equilibrio

La conducta de fase de una solución ideal puede ser descrita combinando las leyes de Raoult y Dalton
para obtener la cantidad de gas y líquido que existen a las condiciones dadas. La Ley de Raoult establece
que la presión parcial de un componente en el gas es igual al producto de la fracción molar de ese
componente en el líquido multiplicado por la presión de vapor del componente puro. La Ley de Raoult
sólo es válida si el gas y las mezclas líquidas son soluciones ideales. La sentencia matemática de la Ley
de Raoult es

p j =x j p vj ,

donde pj representa la presión parcial de componente j en el gas en equilibrio con un líquido de

composición xj. La cantidad de pvj representa la presión de vapor que el componente puro j ejercerían a la
temperatura de interés.
La Ley de Dalton también puede usarse para calcular la presión parcial ejercida por un componente de
una mezcla de gas ideal:

pj= y j p
Las Leyes de Raoult y Dalton los dos representan la presión parcial de un componente en un gas. En el
caso de la Ley de Raoult, el gas debe estar en equilibrio con un líquido. Estas leyes pueden ser
combinadas eliminando presión parcial para dar una ecuación que relaciona las composiciones del gas y
fases líquidas en equilibrio a la presión Y temperatura a las que el equilibrio de gas-líquido existe, como

y j p=x j pvj

yj pvj
= (2-34)
xj p

Ecuación 2-34 puede usarse para calcular la relación de la fracción molar de un componente en el gas a
su fracción molar en el líquido. Para determinar valores de yj y xj, esta ecuación debe ser combinada con
otra ecuación relacionando estas dos cantidades. Semejante ecuación puede desarrollarse a través de la
consideración de un balance de material en el componente jth de la mezcla.

n representa el número total de moles en la mezcla,

nL representa el número total de moles en el líquido,
ng representa el número total de moles en el gas,
zj representa la fracción molar del componente jth en la mezcla total incluyendo las fases líquida y gas,
xj representa la fracción molar del componente jth en el líquido,
yj representa la fracción molar del componente jth en el gas,
zjn representa las moles del componente jth en la mezcla total,
xjnL representa las moles del componente jth en el líquido, y
yjng representa las moles del componente jth en el gas.

Un balance de material en el componente jth resulta en

z j n=x j n L + y j n g . (2-35 )

Combinación de ecuaciones 2-34 y 2-35 para eliminar yj resulta en

zjn
x j= .
pvj
nL + n
p g

Subsecuentemente por definición xj = 1,

z jn
∑ x j=∑ p vj
=1. (2-36)
nL+ ng
p
Las ecuaciones 2-34 y 2-35 también pueden combinarse para eliminar xj con el resultado,
 zjn
∑ y j= ∑ p
=1. (2-37 )
ng + nL
p vj

O las ecuaciones 2-36 o 2-37 pueden ser usadas para calcular las composiciones del gas y fases líquidas
de una mezcla en equilibrio. En cualquier caso una solución del ensayo y error se requiere. El cálculo se
simplifica si una mol de la mezcla total se toma como una base para que el ñL = nL / n, ñg = ng / n, y ñL +
ñg = 1. Esto resulta en la reducción de las ecuaciones 2-36 y 2-37 a

zj
∑ x j= ∑ p
=1 . (2-38 )

(
1+ n̄ g + vj −1
p )
y

zj
∑ y j= ∑ p
=1 . (2-39)
1+ n̄ L +
( )
p vj
−1

El uso de la ecuación 2-38 o 2-39 para predecir conducta de equilibrio gas-líquido severamente se
restringe desde que la Ley de Dalton se basa en la asunción de que el gas se comporta como una solución
ideal de gases ideales. Para los propósitos prácticos la asunción de gas ideal limita el uso de las
ecuaciones a presiones por debajo de aproximadamente 50 psia.

La ley de Raoult también es basada en la asunción de que los líquidos se comportan como una solución
ideal. La conducta de solución ideal sólo es aprovechada solamente si los componentes de la mezcla
líquida son químicamente y físicamente muy similares.

Otra restricción es que un compuesto puro no tiene una presión de vapor a temperaturas sobre su
temperatura crítica. Por lo tanto, el uso de estas ecuaciones se limita a temperaturas menores que la
temperatura crítica del compuesto más volátil en la mezcla. Por ejemplo, si el metano es un componente
de la mezcla, este método no puede usarse a temperaturas sobre -116ºF. Estas restricciones han llevado al
uso de correlaciones basado en observaciones experimentales. Estas correlaciones normalmente
involucran el uso de la relación de equilibrio K, define como

Yj
K j= .
xj

De ecuación 2-34, se ve que el k-valor es definido por la proporción de la presión de vapor del
componente por la presión del sistema, o

p vj Y j
= =K j .
p xj

Sustituyendo esta relación en las ecuaciones 2-38 y 2-39 resultados en

z
∑ x j=∑ 1+n̄ + jK −1 =1 . (2-40 )
g ( j )
 zj
∑ y j=∑ 1
=1 . (2-41)

( )
1+ n̄ L +
Kj
−1

La opción entre las ecuaciones 2-40 y 2-41 es arbitraria; cualquiera puede usarse con igual eficacia.
Cualquier ecuación requiere una solución de ensayo y error. Normalmente, los valores de ensayo
sucesivo de ñg o ñL son seleccionados hasta la suma sea igual a uno. Estas ecuaciones pueden usarse para
calcular la cantidad y composición de una mezcla por el procedimiento siguiente. Los datos requeridos
son la alimentación de la composición, la presión y temperatura en las que la separación tema lugar, y
factores K para cada componente a esta presión y temperatura.

1. Determinar el valor K apropiado para cada componente a p y T de separación.

2. Asumir un valor de ñg y calcular un valor supuesto de ñL de

n̄ L + n̄ g =1
3. Calcular la suma de las fracciones molares de todos los componentes usando el valor de ñg asumidos
en paso 2.

z
∑ x j=∑ 1+n̄ + jK −1
g ( j )
4. Si la suma calculara en paso 3 es igual a uno, la asunción para ñg fue correcta. Si no, asumir un nuevo
valor por el ñg y repetir el paso 3 hasta que la suma sea igual a uno.

Los valores calculados del xj en paso 3 representan la composición de las ñL moles de líquido resultante
de la separación. La composición de las ñg de gas puede calcularse de

y j =K j x j .

El paso 3 puede ser modificad para calcular la suma de las fracciones molares en la fase de gas
asumiendo un valor de ñL y usando la ecuación 2-41.

Ejemplo 2-19

Calcular la cantidad y composición de la fase gaseosa y líquida resultante de la separación flash del
líquido a una presión de 200 psia y una temperatura de 150ºF. Asumir que los gráficos para una presión
de convergencia de 5000 psia son aplicables.

Solución

Las ecuaciones previas también pueden usarse para determinar el punto de la burbuja y presiones de
punto de rocío de una mezcla a una temperatura dada. En el punto de la burbuja, la una cantidad de gas es
despreciable y

z j ≈x j ,
y

y j=x j k j=z j K j (2-42 )

∑ y j=∑ z j K j=1 .
A la presión de punto de rocío,

z j≈ y j ,

z
∑ x j=∑ Kj =1 . (2-43)
j

El punto de la burbuja o presiones de punto de rocío pueden ser calculados por ensayo y error asumiendo
varias presiones y determinando la presión que producirá K-valores que satisfacen Ecuaciones 2-42 o 2-
43.

Las relaciones siguientes ayudarán en la determinación de el punto de la burbuja o punto del rocío.

z
a. Si ∑ z j K j y ∑ Kj son > 1, p está entre pb y pd .
j

b. Si ∑ z j K j <1, p> p b .

z
c. Si ∑ k j <1, p< pd
j

Ejemplo 2-20

Determinar la presión de punto de burbuja de la siguiente mezcla a una temperatura de 150ºF.

Solución

De los ejemplos previos se conoce que la presión del punto de burbuja es mayor a 200 psia. Como una
primera aproximación, usar pb = 250 psia.

Ensayo 1 (pb = 250 psia)

La presión asumida es mayor que la presión del punto de burbuja desde que zjKj  1
Ensayo 2 (pb = 220 psia)

El hecho que zjKj > 1 implica que 220 psia está bajo la presión del punto de burbuja. Por interpolación
lineal, se selecciona la presión del punto de burbuja para un tercer ensayo de 230 psia.

Ensayo 3 (pb = 230 psia)

La presión del punto de burbuja es por lo tanto de 230 psia a una temperatura de 150ºF.

Determinación de las relaciones de equilibrio

Para aplicar los métodos descritos en la sección precedente, es necesario tener acceso a los K-valores para
cada componente como una función de presión y temperatura. El método más ampliamente usado de la
determinación de los K-valores es usar cartas experimentalmente determinadas publicadas por la NGPSA
(Engineering Data Book, Gas Processors Suppliers Asoc., Tulsa, 1972). El problema de usar estas cartas
es que debe conocerse la presión de convergencia, pk, para seleccionar la carta apropiada. En presiones
relativamente bajas los K-valores son muy independientes de la presión de convergencia, pero para
presiones de separación que se acercan presión de convergencia, ellos son muy sensibles a la presión de
convergencia. La determinación de la presión de convergencia apropiada a usar para un sistema dado es
iterativa. Un procedimiento sugerido es:

1. Estime un valor para la presión de la convergencia, p'k.

2. Leer los K-valores para cada componente de las cartas preparadas a p'k.

3. Realizar un cálculo rápido usando los K-valores del paso 2. Esto da la composición líquida.

4. Identificar el componente del hidrocarburo más ligero en la fase líquida que comprende al menos 0.001
de fracción molar; es decir, xj  0.001. Este normalmente será C1 o C2 para los sistemas de gas
condensado.

5. Calcular el peso promedio de temperatura crítica para los componentes más pesados restantes en la
fase líquida. Estos componentes más pesados restantes representan un componente pseudo-pesado en el
pseudo sistema binario.

6. Entrar en figura 2-25 a la temperatura crítica calculada en paso 5 y determina qué componente del
hidrocarburo actual tiene una temperatura crítica cerca al componente pseudo-pesado. Si el valor está
entre dos componentes reales, se requiere interpolación.

7. Trazar el sitio crítico entre el componente ligero y pseudo-pesado a la temperatura de separación y leer
la presión de convergencia pk en este punto. Si p'k  pk del paso 1, p'k = pk y se pasa a punto 2. Repetir
hasta que pk  p'k.
K-valores de las Ecuaciones de Estado

Las relaciones de equilibrio pueden definirse en términos de la fugacidad del componente en una mezcla.
La fugacidad, f, es una cantidad termodinámica definida en términos del cambio en energía libre pasando
de un estado a otro. Se ha mostrado que la fugacidad de un componente en una fase de una mezcla es
igual a la fugacidad de ese componente en la misma fase en el estado puro multiplicado por la fracción
molar del componente en la mezcla. Es decir, para las fases gas y líquido,

f gj ( mix )= y j f gj ( puro )

f Lj ( mix )=x j f Lj ( puro )

Si las fases están en equilibrio, las fugacidades de un componente en las fases líquido y de gas son
iguales, o

y j f gh=x j f Lj

yj f Lj
=K j = (2-44 )
xj f gj

Por consiguiente, los K-valores pueden aproximarse si la fugacidad de cada componente puede calcularse
a una presión y temperatura dada. Esto se puede lograr usando una ecuación de estado, como la ecuación
de Soave-Redlich-Kwong descrita antes. El coeficiente de fugacidad de un componente en una mezcla se
da por

fj bj a 0 .5 b
ln
px j
=
b
( Z−1 )−ln ( Z-B )−
A
B (
2 i0 . 5 − i ln 1+
a b )
B
Z ( ) ( 2-45 )

donde
 a0i . 5 α j0 . 5 T cj / p cj0. 5
=
a0 . 5 ∑ x j α 0. 5 p 0 . 5
j cj

b j T cj / p cj
=
b ∑ x j T cj / pcj

Los valores para Z, A, B, y  son calculados por las ecuaciones 2-22, 2-23, 2-24, y 2-25, respectivamente.
El xj representa la fracción molar del líquido o del gas.

El uso de una ecuación de estado para calcular la conducta de la fase involucra numerosos cálculos y
normalmente se aplica sólo mediante el uso de una computadora. Algunos programas simplificados se
han publicado para el uso en calculadoras de mano por Weber (Predicting Properties of Gas Mixtures,
Chemical Eng. May 19, 1980).

Ajuste de las Propiedades para las Mezclas de Condensado

En muchos casos no se medirán las propiedades totales del pozo o fluido del reservorio, pero se
conocerán las propiedades del gas y del condensado después de la separación. Muchos cálculos de los
reservorios y de los sistemas de tuberías requieren propiedades como la gravedad específica por el
reservorio o fluido del pozo. Los métodos por hacer los cálculos se dan debajo.

Gravedad específica de Mezclas. Formando un balance molal en el separador de efluentes de gas y

líquido, la ecuación siguiente puede ser derivada para calcular la gravedad de una mezcla:

γ g + 4584 γ o /R
γ gm = (2-46 )
1+ V o /R

donde

gm = gravedad específica de la mezcla (aire = 1),

g = gravedad específica del gas del separador,
o = gravedad específica del condensado,
Vo = volumen condensado vaporizando, y
R = relación de producción de gas-condensado.

La gravedad específica del condensado se relaciona a la gravedad API por

141. 5
γ o=
131. 5+API

Vo se obtiene de la figura 2-26 para los casos en los que o es medido a una presión de separador alta. Si
o es medida a las condiciones de tanque, el valor de Vo depende de la presión del separador y se obtiene
de la figura 2-27.
Ejemplo 2-21

Calcular la gravedad específica de un reservorio de gas dado:

g = 0.65 o de S.T. = 0.78 (50 ºAPI)

R = 30000 scf/STB
Presión del separador = 900 psig

Solución

De figura 2-27, Vo = 1110 scf/STB

La ecuación 2-46 puede ser usada para calcular la gravedad específica de un gas conteniendo agua si la
gravedad del gas seco y el contenido de agua son conocidos. El volumen de líquido vaporizando para el
agua es

Para calcular Rw en scf/STB de el contenido de agua en lbm/MMscf, usa