FACULTAD INTEGRAL DEL NORTE

CARRERA DE INGENIERÍA PETROLERA

NOMBRE: LILIANA ZABALA JARE

REGISTRO: 212053264

DOCENTE: ING. EINAR MONTAÑO SUAREZ

FECHA: 16/3/2021
1. INTRODUCCIÓN

Los primeros investigadores en ingeniería de yacimientos se dieron cuenta que

para calcular los volúmenes de petróleo y gas in situ era necesario conocer:
Presión, volumen y temperatura; la presión de saturación o punto de burbujeo; la
cantidad total de gas disuelto en el petróleo como resultado de la liberación de su
gas disuelto.

El gas es el estado de agregación de la materia que no tiene forma ni volumen

propio, es decir, bajo ciertas condiciones de temperatura y presión permanece en
estado gaseoso. Sus moléculas se encuentran prácticamente libres, de modo que
son capaces de distribuirse por todo el espacio en el cual son contenidos.

2. OBJETIVO

Este documento está siendo realizado para ayudar a los estudiantes en cualquier
etapa de la carrera, con fines de instruir y de despejar dudas con respecto a la
compresibilidad del gas natural, los métodos que podemos usar y cómo se deben
usar.

3. DESARROLLO

El gas natural consiste principalmente de más del 70% de metano y proporciones

significativas de etano, propano, butano, pentano y pequeñas cantidades de
hexano, heptano y fracciones más pesadas. Habrá siempre alguna cantidad de
condensado y/o petróleo asociado como el gas.

3.1 COMPOSICIÓN DEL GAS

Composición típica:
Metano  CH4  usualmente > 70%
Etano  C2H6 2 a 10 %
Otros hidrocarburos:
Propano C3H8, Isobutano i  C4H10, Butano normal (n  C4H10), Isopentano i 
C5H12 , Pentano normal n  C5H12, Hexano C6H14, fracciones más pesadas
C7H+16.

Hidrocarburos cíclicos y aromáticos: ocasionalmente pueden ocurrir en pequeñas

proporciones.
Impurezas comunes: nitrógeno N2, dióxido de carbono CO2 y gas sulfhídrico
H2S

3.2 COMPORTAMIENTO DE LOS GASES

En la realidad no existen gases perfectos; sin embargo, muchos gases cerca de la

temperatura y presión atmosféricas se aproximan al comportamiento ideal. La
teoría cinética de los gases postula que el gas está compuesto de una cantidad
de partículas llamadas moléculas para un gas ideal, el volumen de estas
moléculas es insignificante comparado con el volumen total ocupado por el gas,
también se asume que estas moléculas no tienen una fuerza de atracción o
repulsión entre ellas y así se asume que todas las colisiones de las moléculas son
perfectamente elásticas.

3.3 LEY DE LOS GASES IDEALES

El gas ideal puede definirse como el gas cuyo volumen se reduce a la mitad
cuando la presión se dobla, y cuya presión se dobla si, manteniendo su volumen
constante, doblamos su temperatura absoluta. Estos no son más que enunciados
específicos de las Leyes de Boyle, Ley de Charles, Ley de Gay- Lussac y la Ley
de Avogadro.

3.6.1 LEY DE CHARLES

Corresponden a las transformaciones que experimenta un gas cuando la presión

es constante. Así tenemos que:

Eq. (1)

V1 V2
=
T 1 T2

PV=nRT Eq. (2)

Cuando la temperatura se acerca al cero absoluto, todos los gases tienden al

mismo comportamiento.
3.6.2 LEY DE BOYLE

Corresponde a las transformaciones que experimenta un gas cuando su

temperatura permanece constante.

P1V1 = P2V2 Eq. (3)

La curva que describe el gráfico P versus Volumen, corresponde a una isotérmica,

es decir a todos los puntos donde la temperatura es la misma.

3.6.3 LEY DE GAY – LUSSAC

Corresponde a las transformaciones que sufre un gas ideal cuando el volumen

permanece constante. Eq. (4)
P 1 P2
=
T1 T2
3.6.4 LEY DE AVOGADRO

Volúmenes iguales de distintas sustancias gaseosas, medidos en las mismas

condiciones de presión y temperatura, contienen el mismo número de partículas.

PV=nRT Eq. (5)

Donde:
P= Es la presión del gas
V= El volumen del gas
N= El número de moles
T= La temperatura del gas medida en kelvin
R= La constante de los gases ideales

3.4 DENSIDAD DEL GAS

La densidad del gas por definición es la relación entre las masas específicas del
gas y del aire, ambas medidas en las mismas condiciones de presión y
P2
temperatura, esto es: Eq. (6) g=
T2

Admitiéndose comportamiento de gas ideal, en la ecuación 5 el número de moles

n es la relación entre la masa de gas m y su masa molecular M.

Masa específica del aire es: Eq. (7)

P∗28,97
Paire=
RT
Densidad de un gas es: Eq. (8)
PM
RT
P x
P 28,47 M
γ g= g = =
Ṕ a RT 28,97
Donde:

g = Densidad del gas

M = Peso molecular, lbm/lb-mol
28,97 = Peso molecular del aire.

3.5 GASES REALES

El comportamiento de un gas real es diferente a un gas ideal, la razón para esto

es que la ley de los gases perfecto fue derivada bajo la asumición que el volumen
de moléculas es insignificante y no existe atracción o repulsión entre las
moléculas, lo cual no es el caso para gases reales.

De acuerdo con la ley definida por la ecuación 5 para las presiones y

temperaturas de trabajo. Para expresar de forma más real la relación entre las
variables P, V y T, un factor de corrección, denominado factor de compresibilidad
de gas Z, es introducido en la ecuación 5

PV=ZnRT Eq. (9)

Donde para un gas ideal Z=1.

El factor de compresibilidad varia con el cambio de presión y temperatura en la

composición del gas. Esto debe determinarse experimentalmente. Los resultados
de la determinación experimental del factor de compresibilidad son normalmente
dados gráficamente.

3.6 METODO DE OBTENCIÓN DEL FACTOR DE COMPRESIBILIDAD Z

El factor de compresibilidad Z, conocido también como el factor de compresión, es

la razón del volumen molar de un gas con relación al volumen molar de un gas
ideal a la misma temperatura y presión. Es una propiedad termodinámica útil para
modificar la Ley de los Gases Ideales para ajustarse al comportamiento de un gas
real. En general, la desviación del comportamiento ideal se vuelve más
significativa entre más cercano esté un gas a un cambio de fase, sea menor la
temperatura o más grande la presión. Eq. (10)

Vm
Z=
¿¿´ ¿ ¿

3.6.1 CORRELACIONES DE STANDING Y KATZ

La correlación de Standing y Katz favorece en función de presiones y
temperaturas reducidas y fue desarrollada con base en datos experimentales para
gases naturales sin impurezas. Su aplicación para gases ácidos requiere el uso
de factores de corrección para la presencia de CO 2 y H2S.

El procedimiento para la determinación de Z sigue los siguientes pasos:

Paso 1 Determinar las propiedades Pseudo Críticas

 Composición conocida.

De tabla 3 – 4 determinar el Peso Molecular, presión y temperatura Pseudo

Críticas para cada componente.
Obtener la Masa molecular aparente (Ma) de la mezcla. Con la siguiente
ecuación: Ma = yiMi Eq. (11)

Donde:

Ma = Peso molecular aparente de la mezcla.

yi = Fracción molar del componente i.
Mi = Peso molecular del componente i.

Obtener las coordenadas pseudo críticas. Estas pueden ser calculadas a

través de media ponderadas de las coordenadas críticas de cada
componente y su fracción molar en la mezcla: Eq. (12)
nc
Ppc=∑ YiPci
i=1
Eq. (13)
nc
Tpc=∑ YiTci
i=1
Donde:

Pci = Presión Pseudo Crítica del componente i.

Tci = Temperatura Pseudo Crítica del componente i.
Yci = Fracción molar del componente i.
nc = Número de componentes.

 Composición desconocida.
Con la densidad del gas conocida usar la figura 3.2 donde la presión y
temperatura pseudo críticas son dadas. O a través de las siguientes
correlaciones presentadas por Standing (1981):
Ppc = 677 – 15,0 g – 37,5 2g Eq. (14)

Tpc = 168 + 325 g – 12,5 2g Eq. (15)

 Para gas natural seco, es: Ppc = 706 – 51,7 g - 112g Eq. (16)

Tpc = 187 + 330g – 71,52g Eq. (17)

Paso 2 Correcciones de las propiedades Pseudo Críticas.

Debido a la presencia de gases que no son hidrocarburos, utilizamos el factor de

corrección de Wichert y Aziz dado por:

 = 120(A0,9 – A1,6) + 15(B0,5 – B4,0) Eq. (18)

Donde:
 = Factor de ajuste de las propiedades pseudo críticas.
A = Suma de las fracciones molares de H2S + CO2.
B = Fracción molar de H2S.

Calcule la temperatura pseudo crítica ajustada por la siguiente expresión:

Tpc = Tpc -  Eq. (19)
Y a la presión pseudo crítica ajustada por medio de:

P pc Tpc
Ppc=
T pc + B(1−B)
Eq. (20)

El parámetro  puede de ser también obtenido en la figura 4.4

Paso 3 Factor de compresibilidad Z

El factor de compresibilidad Z, de la mezcla gaseosa es obtenido de la carta de

Standing & Katz, figura 3.5 en función de las coordenadas pseudo reducidas, o
sea:
P
Ppr=
Ppc
Eq. (21)
T
Tpr=
Tpc
Eq. (22)
O, si hay presencia de contaminantes:

P
Ppr=
P ' pc
Eq. (23)
T
Tpr=
T ' pc

Eq. (24)
3.6.2 CORRELACIONES DE BRILL & BEGGS
Una modificación de las ecuaciones publicadas por Brill & Beggs (1974) favorece
valores del factor Z con precisión suficiente para la mayoría de los cálculos de
ingeniería:

1− A
Z= A + + CPDpr Eq. (25)
expB

Donde:

A = 1.39(Tpr – 0.92)0.5 – 0.36Tpr – 0.101 Eq. (26)

0.066
B= (0.62 – 0.23Tpr) Ppr +( −0.037)P2pr Eq. (27)
Tpr−0.86

C=0.132 – 0.32logTpr
D=antilog (0.3106 – 0.49Tpr + 0.1824T2pr) Eq. (28)

3.6.3 CORRELACIONES DE DRANCHUK, PURVIS Y ROBINSON

Este método se basa en la ecuación de estado desarrollada por Benedicto, Webb

and Rubin para representar el comportamiento de hidrocarburos leves. Utilizando
esa ecuación, conteniendo ocho constantes características de cada sustancia,
Dranchuk et al. Obtuvieron una ecuación explícita de la siguiente forma: Eq. (29)

A2 A3 A5 2 Pr 5 A 7 Pr 2
Z=1+( A 1+ + )Pr +( A 4 + )Pr + A5A6 + (1+AgPr2)EXP(-AgPr2)
Tr Tr 3 Tr Tr Tr 3
 Pr
Donde: Pr = 0,27
ZTr
Eq.(30)

Y las constantes de la correlación son dadas por:

A1 = 0.31506237 A5 = -0.61232032
A2 = -1.04670990 A6 = -0.10488813
A3 = -0.57832729 A7 = 0.68157001
A4 = 0.53530771 A8 = 0.68446549
3.6.4 CORRELACIONES DE HALL – YARBOROUGH

La ecuación de Hall- Yarborough, desarrollada usando la ecuación de estado de

Starling- Carnahan:
Eq. (31)
−1,2 2
0,06125 p Γt t e ( 1−t )
z=
y

Donde:
Pr = Presión pseudo reducida.
T = Reciproco, de temperatura reducida (Tpc/T).
Y = Densidad reducida con la cual puede ser obtenida como la solución de la
ecuación.

k Y + Y 2 +Y 3 −Y 4
-1,2(1-t)2
F = - 0,06125 Pprte + – (14,76t – 9,67t2 + 4,58t3) Y2 + (90,7t –
(1− y )3
242,2t2 + 42,4t3) Y(2,18+2,82t) = 0 Eq. (32)

Esta ecuación no lineal puede ser convenientemente resuelta usando las técnicas
simples de iteración de Newton – Raphson. Los pasos son:

 Tomar una estimación inicial de Yk, donde k es un contador de iteración.

Y1 = 0,001
 Sustituya este valor en la ecuación 30 sería un valor muy pequeño, el valor
de Fk diferente de cero.
 Utilizando la primera serie de expansión de Taylor, una mejor estimación
de Y puede ser determinada por la ecuación:
Fk
Yk+1= Yk - Eq. (33)
dFk /dY

Donde la expresión general para dF/dK puede ser obtenido como la

derivada de la ecuación 30.

1+ 4 Y + 4 Y 2 −4 Y 3 +Y 4
dF/dY = – (29,52t – 19,52T2 + 9,16t3)Y + (2,18 +
(1− y ) 4
2,82t)(90,7t - 242,2t2 + 42,4t3)Y(1,18+2,82t)


Iterando, usando la ecuación 30 y 31 antes que haya convergencia dentro
de una aproximación satisfactoria, Fk  0.
 Sustituya este valor correcto de Y en la ecuación 29 para determinar el
valor de Z.
3.7 CONCLUSIÓN

El método de Standing y Katz es un método confiable porque busca realizar un

ajuste numérico al método de Standing y Katz para obtener valores de Z por
medio de computadoras.
El método de Brill y Beggs no es válido para temperaturas seudorreducidas fuera
del rango de 1.2 a 2.4
Por el método de Dranchuk, Purvis y Robinson cuando se hace hace la
lecturación de los valores en las curvas de standing y Katz hay bastante error por
lo tanto no es un método confiable.
Al igual que el anterior este método de Hall y Yarborough también tiene errores al
momento de leer los valores en las curvas de Standing y Katz.

4.7 BIBLIOGRAFIA

[ CITATION Dav19 \l 16394 ]

[ CITATION Cra68 \l 16394 ]
[ CITATION Cam01 \l 16394 ]