Calculo de La Densidad Del Gas Natural Basado en Su Gravedad Especifica

UNIVERSIDAD MAYOR DE SAN ANDRES
FACULTAD DE INGENIERIA
INGENIERIA PETROLERA
PET-230
PROGRAMACION APLICADA
“CALCULO DE LA DENSIDAD DEL GAS

NATURAL BASADO EN SU GRAVEDAD
ESPECIFICA”
ALUMNA: LLAPACU BALTAZAR MARY ISABEL
DOCENTE: ING HERMAS HERRERA CALLEJAS

Universidad Mayor de San Andrés Programación Aplicada
Ingeniería Petrolera PET – 230
CALCULO DE LA DENSIDAD DEL GAS NATURAL
BASADO EN SU GRAVEDAD ESPECIFICA
1. INTRODUCCIÓN.
Generalmente el gas natural es una mezcla de gases de hidrocarburos y no hidrocarburos. Los

hidrocarburos gaseosos que normalmente se encuentran en el gas natural son metano, etano,
propano, butano, pentano y pequeñas porciones de hexano y superiores. Los gases no
hidrocarburos (o impurezas) incluyen el dióxido de carbono, sulfuro de hidrógeno y nitrógeno.
El conocimiento de las relaciones presión-volumen-temperatura (PVT) y las demás
propiedades físicas y químicas de los gases es esencial para resolver problemas en ingeniería de
reservorios del gas natural. Estas propiedades incluyen:
 Peso molecular aparente, Ma
 Gravedad específica, γg
 Factor de compresibilidad, z
 Densidad, ρg
Comportamiento de los Gases Reales.
Al tratar los gases a muy baja presión, la relación de gases ideales es una herramienta conveniente y
generalmente satisfactoria. A presiones más elevadas, el uso de la ecuación de estado de los gases
ideales puede conducirnos a errores tan grandes como del orden de 500%, comparados a errores de 2
a 3% a presión atmosférica.
Básicamente la desviación de los gases reales de las condiciones de la ley de gas ideal se incrementa
a medida que crece la presión y temperatura y varía ampliamente con la composición del gas. Los
gases reales se comportan de modo diferente a los gases ideales. La razón para esto es que la ley de
los gases perfectos fue derivada bajo la suposición de que el volumen de las moléculas es
insignificante y que no hay atracción o repulsión molecular entre ellas. Este no es el caso para los
gases reales.
Se han desarrollado numerosas ecuaciones de estado con el intento de correlacionar las variables de
presión-volumen-temperatura para gases reales con datos experimentales. Con el fin de expresar una
relación más exacta entre las variables p, V y T, un factor de corrección llamado factor de
compresibilidad de los gases, factor de desviación de los gases, o simplemente factor Z debe
introducirse en la ecuación de estado de los gases ideales para tomar en cuenta el alejamiento de los
gases de la idealidad. La ecuación tiene la siguiente forma:
pV= znRT
Donde p = presión absoluta, psia
V = volumen, ft3
T = temperatura absoluta, °R
n = numero de moles de gas, lb-mole
R = la constante universal de los gases que, par alas unidades anteriores, tiene el valor de
10.73 psia ft3/lb-mole °R
Significado del valor de Z.
Univ. Mary Isabel Llapacu Baltazar 1

Si el valor de Z es igual a 1 esto indica que el gas se comporta como ideal. Si el valor de Z es mayor
o menor que 1 el gas se comporta como un gas real. Mientras mas grande sea la desviación del valor
de Z con respecto a 1, mayor es la desviación del comportamiento respecto al comportamiento ideal
del gas.
Las relaciones de cartas dadas por Standing – Katz para sistemas de condensado de gas y gas,
pueden ser descritas en una nueva ecuación.
Esta ecuación solo necesita ingresar los valores de presión y temperatura, acompañadas con las
fracciones molares del gas, para poder encontrar al factor de desviación z.
Naturalmente, no se trata de una sola ecuación empírica, sino que es el resultado de los análisis
hechos por Yarborough – Hall, de las cartas analizadas por Standing – Katz.
Lo que se quiere realizar es un método iterativo, para el cálculo de z, para lo cual se basaron en la
ecuación de estado.
La facilidad de este método se debe a que los cálculos los realiza el computador, con lo cual, además
de ser un método rápido, es también un método que determina el factor de desviación del gas con
gran aproximación; debido a que realiza correcciones por la presencia de impurezas en el gas (ácido
sulfhídrico, anhídrido carbónico) como los más importantes en cuanto a impurezas.
Densidad
Debido a que la densidad se define como la masa por unidad de volumen de la sustancia,
la ecuación 5.3 puede reordenarse para estimar la densidad del gas a cualquier presión y
temperatura:
m pM
ρ g= =
V RT 1
donde ρg = densidad del gas, lb/ft3
Debe resaltarse que lb se refiere a lbs masa en cualquiera de las posteriores discusiones
de densidad.
Gravedad específica del gas
La gravedad específica se define como la relación de la densidad del gas a la densidad del
aire. Ambas densidades se miden o expresan a la misma presión y temperatura. Normalmente la
presión estándar Psc y la temperatura estándar Tsc se usan al definir la gravedad específica del gas:
ρg
γ g=
ρaire 2
Asumiendo que el comportamiento tanto de la mezcla de gas como del aire se describe por la
ecuación de gas ideal, la gravedad específica puede entonces expresarse como:
p sc M a
RT sc
γ g=
p sc M aire
RT sc
o

Ma Ma
γ g= =
M aire 28 . 96 3
donde γg = gravedad específica del gas
ρaire = densidad del aire
Maire = peso molecular aparente del aire = 28.96
Ma = peso molecular aparente del gas
Psc = presión estándar, psia
Tsc = temperatura estándar, °R
1. JUSTIFICACIÓN.
Conociendo la densidad de un gas, podemos determinar si un gas es más denso que el aire, o no.
Este cálculo nos permite de preveer si el gas se acumula en proximidad de la superficie por el
contrario
Si es menos denso se acumulara en la parte superior.
El gas natural es más ligero que el aire ( su densidad relativa es 0.61,aire =11.0) y a pesare de sus
altos niveles de inflamabilidad y explosividad las fugas o emisiones se disipan rápidamente en las
acpas superiores de la atmosfera, dificultando la formación de mezclas explosivas en el aire.
2. OBJETIVOS.
2.1. OBJETIVO GENERAL
Determinar la densidad del gas natural basado en su gravedad especifica.
2.2. OBJETIVOS ESPECÍFICOS
Determinar el factor z de los gases, mediante el método de Hall – Yarborough

Aplicar los conceptos iteración de programación, para resolver el método.
Aplicar el método de Newton – Raspón para resolver ecuaciones no lineales.
3. MARCO TEORICO.
Los primeros estudios fisicoquímicos realizados a partir del siglo 17 involucraron el estudio de los
gases dentro de un reservorio, debido a que éstos responden en forma más dramática a cambios en el
ambiente que los líquidos y los sólidos. La mayoría de estos primeros estudios estuvieron
focalizados en las variaciones de presión, temperaturas y volumen de una determinada porción de
gas (relaciones p-V-T).
Es asi que se establezca en forma breve las diferencias entre gas ideal y un gas real
 Para un gas ideal la variable "z" siempre vale uno, en cambio para un gas real, "z" tiene que
valer diferente que uno.

 La ecuación de estado para un gas ideal, prescinde de la variable "z" ya que esta para un gas
ideal, vale uno. Y para un gas real, ya que esta variable tiene que ser diferente de uno, así
que la formula queda de esta forma: pV=znRT.
 La ecuación de van der Waals se diferencia de las de los gases ideales por la presencia de dos
términos de corrección; uno corrige el volumen, el otro modifica la presión.
 Los gases reales, a presiones y temperaturas cercanas a las ambientales, actúan como gases
ideales.
De esta simple ley se deducen las isotermas de un gas ideal:
El modelo de “gas ideal” permite definir un marco de referencia para estudiar el comportamiento de
los gases. En algunas ocasiones, podremos modelar los gases utilizando Leyes Ideales; sin embargo,
es de gran importancia tener una noción de las desviaciones que sufren éstos bajo determinadas
condiciones de temperatura, presión y volumen. Los gases naturales o reales presentan las siguientes
desviaciones del comportamiento ideal:
Esta desviaciones aparecen producto de la diferencia de volumen, por lo que definiremos el factor
de compresibilidad (Z), que corresponde a una medida de la “no-idealidad” en el comportamiento de
un gas:

Para un Gas Ideal, el factor de compresibilidad es unitario, mientras que para Gases Reales es mayor
o menor que 1. Ejemplos para el H2O, CO2 y O2 gaseosos:
“Estilizando” un poco el diagrama anterior, podemos generalizar las isotermas PV para un gas real
(nótese que a altas, T extremo superior derecho, el gas real se aproxima al comportamiento ideal):
Importante:
- Los gases reales pueden considerarse como ideales a bajas presiones y altas
temperaturas
- El volumen de las moléculas sí importa, no es despreciable
- Las fuerzas de interacción entre moléculas de los gases influye
Por fortuna, las ecuaciones termodinámicas de estado para gases reales también pueden utilizarse
para los fluidos supercríticos. Cada especie de gas puro (agua, dióxido de carbono, oxígeno,
metano , etano ,propano ,etc.) tiene un único punto crítico, sobre el cual desaparecen los límites
entre el líquido y el vapor. Para el punto crítico se tiene que:

Los gases no se pueden licuefarce a temperaturas por encima de la temperatura crítica porque en
este punto las características de los gases y de los líquidos son las mismas, y no hay base sobre la
cual distinguir entre los gases y los líquidos. La presión del vapor de un líquido a la temperatura
crítica se llama la presión crítica (Pc). La presión del vapor de un líquido nunca es más grande que
esta presión crítica.
El gas se define como un fluido homogéneo de baja viscosidad y densidad que no tiene un volumen
definido sino que se expande para llenar completamente el recipiente que lo contiene. Generalmente
el gas natural es una mezcla de gases de hidrocarburos y no hidrocarburos. Los hidrocarburos
gaseosos que normalmente se encuentran en el gas natural son metano, etano, propano, butano,
pentano y pequeñas porciones de hexano y superiores. Los gases no hidrocarburos (o impurezas)
incluyen el dióxido de carbono, sulfuro de hidrógeno y nitrógeno.
El conocimiento de las relaciones presión-volumen-temperatura (PVT) y las demás propiedades
físicas y químicas de los gases es esencial para resolver problemas en ingeniería de reservorios del
gas natural. Estas propiedades incluyen:
 Peso molecular aparente, Ma
 Gravedad específica, γg
 Factor de compresibilidad, z
 Densidad, ρg
 Volumen específico, V
 Coeficiente de compresibilidad isotérmica del gas, cg
 Factor volumétrico de formación del gas, Bg
 Factor de expansión del gas, Eg
 Viscosidad, μg
Estas propiedades del gas pueden obtenerse de medidas directas en laboratorio o por predicción de
expresiones matemáticas generalizadas. Revisaremos las leyes que describen el comportamiento
volumétrico de los gases en términos de presión y temperatura y correlaciones matemáticas que son
ampliamente usadas en la determinación de las propiedades físicas de los gases naturales.

Comportamiento de los Gases Ideales.
La teoría cinética de los gases postula que los gases son compuestos de un gran número de partículas
llamadas moléculas. Para un gas ideal, se asume que 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 fuerzas atractivas o repulsivas entre ellas, y que todas las colisiones de las moléculas son
perfectamente elásticas.
Basado en la teoría cinética de los gases anterior, la ecuación matemática denominada ecuación de
estado puede derivarse para expresar la relación existente entre la presión p, el volumen V y la
temperatura T para una determinada cantidad de moles de gas n. Esta relación para gases perfectos
se llama la ley de los gases ideales y se expresa matemáticamente por la siguiente ecuación:
pV = nRT
Donde p = presión absoluta, psia

V = volumen, ft3
n = numero de moles de gas, lb-mole
R = la constante universal de los gases que, par alas unidades anteriores, tiene el valor de
10.73 psia ft3/lb-mole °R
Comportamiento de los Gases Reales.
Básicamente la desviación de los gases reales de las condiciones de la ley de gas ideal se incrementa
a medida que crece la presión y temperatura y varía ampliamente con la composición del gas. Los
gases reales se comportan de modo diferente a los gases ideales. La razón para esto es que la ley de
los gases perfectos fue derivada bajo la suposición de que el volumen de las moléculas es
insignificante y que no hay atracción o repulsión molecular entre ellas. Este no es el caso para los
gases reales.
compresibilidad de los gases, factor de desviación de los gases, o simplemente factor Z debe
introducirse en la ecuación de estado de los gases ideales para tomar en cuenta el alejamiento de los
gases de la idealidad. La ecuación tiene la siguiente forma:
pV= znRT
Donde el factor de Comportamiento de los Gases Reales.

compresibilidad de los gases
, factor de desviación de los gases, o simplemente factor Z debe introducirse en la ecuación de
estado de los
gases ideales para tomar en cuenta el alejamiento de los gases de la idealidad. La ecuación tiene la
siguiente forma:
compresibilidad del gas z es una cantidad adimensional y se define como la relación del volumen
actual de n-moles de gas a T y p al volumen ideal del mismo número de moles a la misma T y p:
V V
z  actual 
Videal (nRT ) / p
El Factor z
Estudios de los factores de compresibilidad del gas natural de varias composiciones han mostrado
que los factores de compresibilidad pueden generalizarse con suficiente exactitud para la mayoría de
los propósitos de ingeniería cuando son expresados en términos de las dos siguientes propiedades
adimensionales:
 Presión seudo reducida

 Temperatura seudo reducida
Estos términos adimensionales se definen por las siguientes expresiones:

p
p pr 
p pc
T
T pr 
T pc
Donde p = presión del sistema, psia
ppr = presión seudo-reducida, adimensional
T = Temperatura del sistema, °R
Tpr = Temperatura seudo-reducida, adimensional
ppc, Tpc = Presión y Temperatura seudo-críticas.
Los parámetros pseudoreducidos pueden ser predecidos a través del valor de la gravedad
específica y las ecuaciones de Brown (1948), presentados en forma gráfica con una
aproximación conveniente de las propiedades pseudocríticas de los gases cuando solo
tenemos este dato disponible. Las expresiones matemáticas que pueden emplearse para
efectuar el cálculo son las siguientes:
ECUACIONES DE BROWN
En base a gravedad específica:
Sistemas de gas natural:

2
PPC =677+15. 0γ g −37 . 5γ g
2
T PC =168+325 γ g −12. 5 γ g
Donde:
Tpc =Temperatura pseudo crítica, °R
Ppc = Presión Pseudo crítica, Psia
g = Gravedad específica de mezcla de gas
Limitaciones:
Máximo.
5% N2
2% CO2
2% H2S
Efecto de los componentes no hidrocarburíferos en el factor z
Los gases naturales frecuentemente contienen componentes no hidrocarburíferos como el

nitrógeno, dióxido de carbono y sulfuro de hidrógeno. Podemos clasificar al gas natural en
amargo (si contiene elevadas cantidades de ácido sulfhídrico) o dulce (si el contenido de
ácido sulfhídrico es bajo). Ambos tipos de gas natural pueden o no contener nitrógeno,
dióxido de carbono o ambos.
La común ocurrencia de pequeños porcentajes de nitrógeno y dióxido de carbono, que en

parte es considerado en las ecuaciones anteriormente presentadas, llegando a considerar
como admisible una concentración inferior al 10% de estos componentes, una vez que
esta rebasa estos límites se presentan errores y consideramos esta presencia como
elevada, por lo que se hicieron modificaciones para efectuar el cálculo correspondiente.
Métodos de Ajuste para Mezclas con elevado contenido de compuestos No

Hidrocarburíferos
Existen dos métodos que han sido desarrollados para ajustar las propiedades
pseudocríticas del gas considerando la presencia de los components no hidrocarburíferos,
estos dos métodos son:
• Método de corrección de Wichert-Aziz
• Método de corrección de Carr-Kobayashi-Burrows
El Método de Corrección de Carr-Kobayashi-Burrows
Paso 1. Conociendo la gravedad específica del gas natural, calcular la presión y

temperatura pseudocritica, aplicando las ecuaciones de Brown.
Paso 2. Ajustar los valores estimados de las propiedades pseudocriticas, empleando las
siguientes expresiones:

Donde
T’pc = temperatura pseudocritica, °R
Tpc = temperatura pseudocritica no ajustada, °R
yCO2 = fracción molar de CO2
yH2S = fracción molar de H2S
yN2 = fracción molar de N2
p’pc = presión pseudocritica ajustada, psia
ppc = presión pseudocritica no ajustada, psia
Paso 3. Usando los valores pseudocríticos ajustados calcular las propiedades
pseudoreducidos.
Paso 4. Calcular el factor Z de la gráfica de Standing Katz o por correlaciones.
Debe destacarse que estas propiedades pseudo críticas, como, p pc y Tpc, no representan las
propiedades críticas actuales de la mezcla de gas. Estas pseudo propiedades se usan como
parámetros de correlación en la generación de propiedades del gas.
4. DESARROLLO.
Paso 1. Conociendo la gravedad específica del gas natural, calcular la presión y

temperatura pseudocritica, aplicando las ecuaciones de Brown.
Paso 2. Usando los valores pseudocríticos ajustados calcular las propiedades
pseudoreducidos.
Paso 3. Calcular el factor Z de la gráfica de Standing Katz o por correlaciones.
El Método de Hall-Yarborough
Hall y Yarborough (1973) presentaron una ecuación de estado que representa exactamente la gráfica
del factor z de Standing y Katz. La expresión propuesta se basa en las ecuaciones de estado de
Starling-Carnahan. Los coeficientes de la correlación fueron determinados ajustándolos a los datos
tomados de la gráfica del factor z de Standing y Katz. Hall y Yarborough propusieron la siguiente
formula matemática:
0. 06125 p pr t [−1. 2( 1−t )2 ]

z= [ Υ
e ] …………….(5)
Donde ppr = presión seudo reducida

t = reciproca de la temperatura seudo reducida, Tpc/T
Y = densidad reducida que puede obtenerse como solución de la siguiente ecuación:
Υ +Υ 2 +Υ 3 +Υ 4
F( Υ )=X 1+ 3
−( X 2 )Υ 2 +( X 3 )Υ X 4 =0
( 1−Υ ) ……………(6)
Donde:
X 1 =−0 . 06125∗Psr +T *exp ((8−1. 2)∗(1−T )2 )
X 2 =14 . 76∗T−9 .76∗T 2 +4 .58∗T 3

2 3
X 3 =90 . 7∗T−242. 2∗T + 42. 4∗T
X 4 =2. 18+2 .82∗T
La ecuación anterior es una ecuación no lineal y puede ser convenientemente resuelta para la
densidad reducida Y usando la técnica de iteración de Newton-Raphson.
El procedimiento computacional de resolver la ecuación principal a cualquier presión seudo
reducida especificada ppr y temperatura Tpr se resume en los siguientes pasos:
Paso i. Hacer una estimación inicial del parámetro desconocido, Yk, donde k es un contador de
iteración. Un intento inicial apropiado de Y se da por la siguiente relación:
[ −1 . 2( 1−t )2 ]
Υ k =0 . 0125 p pr te ………(7)
Paso ii. Sustituir este valor inicial en la ecuación (6) y evaluar la función no lineal. A no ser que el
valor correcto de Y se haya seleccionado inicialmente, la ecuación (6) tendrá un valor un valor
diferente de cero para F(Y):
Paso iii.. Una nueva estimación mejorada de Y, como Yk+1, se calcula de la siguiente expresión:
k +1 k f (Υ k )
Υ =Υ −
f ' (Υ k )
Donde f’(Yk) se obtiene evaluando la derivada de la ecuación 5.38 a Yk, o:
2 3 4
1+4 Υ +4 Υ +4 Υ +Υ ( X 4−1 )
f ' (Υ )= 4
−2( X 2)Υ +( X 3 )( X 4 )Υ
(1−Υ )
Paso iv. Los pasos 2–3 se repiten n veces, hasta que el error, abs(Y k – Yk+1), llegue a ser más
pequeño que una tolerancia establecida, como 10-12:
Paso v. El valor correcto de Y entonces se usa para evaluar la ecuación (1) para el factor de
compresibilidad.
Hall y Yarborough resaltaron que el método no se recomienda para su aplicación si la temperatura

seudo reducida es menor a 1.
Paso 4. Conociendo el factor de compresibilidad calcular la densidad del gas natural,

utilizando la siguiente ecuación:
28.96 P
ρ g= ×
10 73 Z ×T
Donde P = presión absoluta, psia
Z= factor de compresibilidad

5. ALGORITMO.
Diagrama de Flujo para Encontrar Z por Correlacion de Yarborough - Hall

inicio
Def X1,X2,X3,X4,Tsr,Psr,T,z D, ite, i,

sw,FD,DFD,Z
Leer Tsr, Psr
T = 1 / Tsr
X1 = (-0.03125) * Psr * T * Exp((-1.2) * (1 - T) ^ 2)

X2 = (14.76 * T) - (9.76 * T ^ 2) + (4.58 * T ^ 3)
X3 = (90.7 * T) - (242.2 * T ^ 2) + (42.4 * T ^ 3)
X4 = 2.18 + (2.82 * T)
ite = 100
D = 0.27 * Psr * T * Exp((-1.2) * (1 - T) ^ 2)
sw = False
i=1 To ite
FD = X1 + ((D + D ^ 2 + D ^ 3 - D ^ 4) / (1 - D) ^ 3) - (X2 * D ^ 2) + (X3 * D ^ X4)

DFD = ((1 + (4 * D) + (4 * D ^ 2) - (4 * D ^ 3) + D ^ 4) / (1 - D) ^ 4) - (2 * X2 * D) + (X3 * X4 * (D ^ (X4 - 1)))
D1 = D - (FD / DFD)
sw = True
Abs(D1 - D) <= (10 ^ (-6)) si
ite = i
no
no sw = True si
D = D1
Mensaje: Z = 0.06125 * ((Psr * T) / D)

"error, corrija datos" * Exp((-1.2) * ((1 - T) ^ 2))
i
YARBOROUGH = Z
fin

6. EL PROGRAMA EN CODIGO VISUAL BASIC.
Option Explicit
Dim f%, Nc%
Dim s#
Dim Ma#
Dim y(20) As Double
Dim pm(30) As Double
Dim tc(20) As Double
Dim pc(20), i As Double
Dim pcr As Double
Dim tcr As Double
Dim P, T As Double
Dim Z As Double, Sw As Boolean
Dim Tsr As Double, Psr As Double
Dim X1 As Double
Dim X2 As Double
Dim X3 As Double
Dim X4 As Double
Dim ite As Double
Dim D As Double
Dim FD As Double
Dim DFD As Double
Dim D1 As Double
Dim Yarborough As Double
Private Sub Command1_Click()

ge = Val(Text3.Text)
P = Val(Text1.Text)
T = Val(Text2.Text)
If ge < 1 Then
pc = 667 + 15 * ge - 37.5 * ge ^ 2
tc = 168 + 325 * ge - 12.5 * ge ^ 2
Text7.Text = pc
Text8.Text = tc
pr = P / pc
tr = T / tc
Text11.Text = pr
Text12.Text = tr
Text6.Text = Format(Psr, "##,####0.0000")
Text7.Text = Format(Tsr, "##,####0.0000")
If Tsr > 1 Then
T = 1 / Tsr
X1 = (-0.06125) * Psr * T * Exp((-1.2) * ((1 - T) ^ 2))

X2 = (14.76 * T) - (9.76 * T ^ 2) + (4.58 * T ^ 3)

X3 = (90.7 * T) - (242.2 * T ^ 2) + (42.4 * T ^ 3)
X4 = 2.18 + (2.82 * T)
ite = 100
Sw = False
D = 0.27 * Psr * T * Exp((-1.2) * (1 - T) ^ 2)
For i = 1 To ite
FD = X1 + ((D + D ^ 2 + D ^ 3 - D ^ 4) / (1 - D) ^ 3) - (X2 * D ^ 2) + (X3 * D ^ X4)

DFD = ((1 + (4 * D) + (4 * D ^ 2) - (4 * D ^ 3) + D ^ 4) / (1 - D) ^ 4) - (2 * X2 * D) + (X3 * X4 *
(D ^ (X4 - 1)))
D1 = D - (FD / DFD)
If Abs(D1 - D) <= (10 ^ -6) Then
Sw = True
ite = i
Exit For
Else
D = D1
End If
Next i
If Sw = True Then
Z = 0.06125 * ((Psr * T) / D) * Exp((-1.2) * ((1 - T) ^ 2))
Else
MsgBox "Error, corrija datos"
End If
Text8.Text = Z
Else
MsgBox "LA TEMPERATURA SEUDOREDUCIDA ES MENOR A 1"
End If
End Sub

ge = Val(Text3.Text)
P = Val(Text1.Text)
T = Val(Text2.Text)
dg = (2.7 * P * ge) / (Z * T)

Text14.Text = dg
End Sub

End
End Sub
Private Sub Text1_KeyPress(KeyAscii As Integer)

If InStr("0123456789." & Chr(8) & Chr(13), Chr(KeyAscii)) = 0 Then
KeyAscii = 0
End If
If KeyAscii = 13 Then
Text2.SetFocus
End If
End Sub

KeyAscii = 0
End If
If KeyAscii = 13 Then
Text3.SetFocus
End If
End Sub

KeyAscii = 0
End If
End Sub
7. DATOS DE PRUEBA.
El gas natural boliviano tiene una gravedad especifica de 0.6575, para una presión de 3800 psia y un
temperatura de 670 ºR, determinar su densidad.
8. RESULTADOS.
 Arrancamos el programa.

 Introducimos datos.
 Calculo de la densidad del gas natural.

9. CONCLUSIONES.
Se utilizaron satisfactoriamente las ecuaciones de iteración de Newton – Raphson., para

determinar el factor de compresibilidad z, y de esta manera llegar a calcular la densidad del gas
natural.
Este método tiene sus limitaciones, para Tsr menores a 1, el método arroja datos erroneos.
El diseño del algoritmo y la aplicación de este fue todo un éxito ya que se logro el proposito
esperado además se obtuvo parámetros básicos necesarios en campo.
La ejecución del programa con un problema práctico de pozo gasifero demostró la veracidad del
funcionamiento del programa.
10. BIBLIOGRAFÍA.
Ingeniería de Reservorios – Ing. Rolando Camargo Gallegos

Ingeniería de Reservorios I – Ing. Hermas Herrera Callejas
BRADY.(1993) Química básica. V ed. Noruega editores. Méjico.
http://www.tannerm.com/physical/vanderWaals/vander.htm
GLASTONE. Fundamentos de Físicoquímica

11. ANEXOS.
Figura 1 Gráfico de factores de compresibilidad de Standing y Katz. (Cortesía de GPSA y GPA Engineering Data Book,
EO Edición, 1987.)

Contenido
1. INTRODUCCIÓN.....................................................................................................................1
2. JUSTIFICACIÓN.......................................................................................................................3
3. OBJETIVOS...............................................................................................................................3
3.1. OBJETIVO GENERAL......................................................................................................3
3.2. OBJETIVOS ESPECÍFICOS.............................................................................................3
4. MARCO TEORICO...................................................................................................................3
5. DESARROLLO.......................................................................................................................10
6. ALGORITMO..........................................................................................................................11
7. EL PROGRAMA EN CODIGO VISUAL BASIC.................................................................13
8. DATOS DE PRUEBA.............................................................................................................15
9. RESULTADOS........................................................................................................................15
10. CONCLUSIONES................................................................................................................17
11. BIBLIOGRAFÍA..................................................................................................................17
12. ANEXOS..............................................................................................................................17

Calculo de La Densidad Del Gas Natural Basado en Su Gravedad Especifica

Cargado por

Copyright:

Formatos disponibles

Calculo de La Densidad Del Gas Natural Basado en Su Gravedad Especifica

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Calculo de La Densidad Del Gas Natural Basado en Su Gravedad Especifica

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD MAYOR DE SAN ANDRES

“CALCULO DE LA DENSIDAD DEL GAS

DOCENTE: ING HERMAS HERRERA CALLEJAS

CALCULO DE LA DENSIDAD DEL GAS NATURAL

BASADO EN SU GRAVEDAD ESPECIFICA

Generalmente el gas natural es una mezcla de gases de hidrocarburos y no hidrocarburos. Los

Significado del valor de Z.

Univ. Mary Isabel Llapacu Baltazar 1

Univ. Mary Isabel Llapacu Baltazar 2

2.1. OBJETIVO GENERAL

Determinar la densidad del gas natural basado en su gravedad especifica.

2.2. OBJETIVOS ESPECÍFICOS

Determinar el factor z de los gases, mediante el método de Hall – Yarborough

Univ. Mary Isabel Llapacu Baltazar 3

De esta simple ley se deducen las isotermas de un gas ideal:

Univ. Mary Isabel Llapacu Baltazar 4

Univ. Mary Isabel Llapacu Baltazar 5

Univ. Mary Isabel Llapacu Baltazar 6

Comportamiento de los Gases Ideales.

Donde p = presión absoluta, psia

Comportamiento de los Gases Reales.

Donde el factor de Comportamiento de los Gases Reales.

Univ. Mary Isabel Llapacu Baltazar 7

 Presión seudo reducida

Estos términos adimensionales se definen por las siguientes expresiones:

Univ. Mary Isabel Llapacu Baltazar 8

Efecto de los componentes no hidrocarburíferos en el factor z

Los gases naturales frecuentemente contienen componentes no hidrocarburíferos como el

La común ocurrencia de pequeños porcentajes de nitrógeno y dióxido de carbono, que en

Métodos de Ajuste para Mezclas con elevado contenido de compuestos No

El Método de Corrección de Carr-Kobayashi-Burrows

Paso 1. Conociendo la gravedad específica del gas natural, calcular la presión y

Univ. Mary Isabel Llapacu Baltazar 9

Paso 1. Conociendo la gravedad específica del gas natural, calcular la presión y

0. 06125 p pr t [−1. 2( 1−t )2 ]

Donde ppr = presión seudo reducida

Univ. Mary Isabel Llapacu Baltazar 10

Donde f’(Yk) se obtiene evaluando la derivada de la ecuación 5.38 a Yk, o:

Hall y Yarborough resaltaron que el método no se recomienda para su aplicación si la temperatura

Paso 4. Conociendo el factor de compresibilidad calcular la densidad del gas natural,

Univ. Mary Isabel Llapacu Baltazar 11

Diagrama de Flujo para Encontrar Z por Correlacion de Yarborough - Hall

Univ. Mary Isabel Llapacu Baltazar 12

Def X1,X2,X3,X4,Tsr,Psr,T,z D, ite, i,

Leer Tsr, Psr

X1 = (-0.03125) * Psr * T * Exp((-1.2) * (1 - T) ^ 2)

FD = X1 + ((D + D ^ 2 + D ^ 3 - D ^ 4) / (1 - D) ^ 3) - (X2 * D ^ 2) + (X3 * D ^ X4)

Mensaje: Z = 0.06125 * ((Psr * T) / D)

Univ. Mary Isabel Llapacu Baltazar 13

6. EL PROGRAMA EN CODIGO VISUAL BASIC.

Private Sub Command1_Click()

If Tsr > 1 Then

X1 = (-0.06125) * Psr * T * Exp((-1.2) * ((1 - T) ^ 2))

Univ. Mary Isabel Llapacu Baltazar 14

X2 = (14.76 * T) - (9.76 * T ^ 2) + (4.58 * T ^ 3)

D = 0.27 * Psr * T * Exp((-1.2) * (1 - T) ^ 2)

FD = X1 + ((D + D ^ 2 + D ^ 3 - D ^ 4) / (1 - D) ^ 3) - (X2 * D ^ 2) + (X3 * D ^ X4)

If Abs(D1 - D) <= (10 ^ -6) Then