Calculo Del Factor Z de Los Gases Por El

Universidad Mayor de San Andrs Programacin Aplicada
Ingeniera Petrolera PET 230
CALCULO DEL FACTOR Z DE LOS GASES POR EL METODO

HALL YARBOROUGH
I. INTRODUCCIN.
El factor de compresibilidad Z, es un factor de correccin, que se introduce en la ecuacin de estado

de gas ideal para modelar el comportamiento de los gases reales, los cuales se pueden comportar
como gases ideales para condiciones de baja presin y alta temperatura, tomando como referencia
los valores del punto crtico, es decir, si la temperatura es mucho ms alta que la del punto crtico, el
gas puede tomarse como ideal, y si la presin es mucho ms baja que la del punto crtico el gas
tambin se puede tomar como ideal.
La desviacin de un gas respecto de su comportamiento ideal se hace mayor cerca del punto crtico.
Remitindonos a la seccin de Gases Ideales tenemos:
Introduciendo el factor de correccin Z:
Por lo tanto:
El factor Z tambin se puede entender como:
Donde
Vactual: volumen especfico que se tiene del gas.
Videal: volumen especfico del gas tomado de la ec. de gas ideal.
Significado del valor de Z.

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 desviacin del valor
de Z con respecto a 1, mayor es la desviacin del comportamiento respecto al comportamiento ideal
del gas.
Univ. Henry Isaac Cruz Lucas 1

Las relaciones de cartas dadas por Standing Katz para sistemas de condensado de gas y gas,
pueden ser descritas en una nueva ecuacin.
Esta ecuacin solo necesita ingresar los valores de presin y temperatura, acompaadas con las
fracciones molares del gas, para poder encontrar al factor de desviacin z.
Naturalmente, no se trata de una sola ecuacin emprica, sino que es el resultado de los anlisis
hechos por Yarborough Hall, de las cartas analizadas por Standing Katz.
Lo que se quiere realizar es un mtodo iterativo, para el clculo de z, para lo cual se basaron en la
ecuacin de estado.
La facilidad de este mtodo se debe a que los clculos los realiza el computador, con lo cual, adems
de ser un mtodo rpido, es tambin un mtodo que determina el factor de desviacin del gas con
gran aproximacin; debido a que realiza correcciones por la presencia de impurezas en el gas (cido
sulfhdrico, anhdrido carbnico) como los ms importantes en cuanto a impurezas.
II. JUSTIFICACIN.
El factor de compresibilidad, z, de un gas natural es una medida de su desviacin respecto al

comportamiento del gas ideal. Su valor esta usualmente entre 0,8 y 1,2, pero puede estar tan bajo
como 0,3 y tan alto como 2,0. se usa en el calculo de las seudo-presiones del gas, y en conversiones
de volmenes a caudales de gas de condiciones estndar a condiciones de reservorio ( y viceversa).
Es a veces llamado factor de desviacin, y tambin factor de super-compresibilidad que no debe
confundirse con la compresibilidad (cambio de volumen por unidad de volumen por cambio
unitario en presin).
Se describe y aplica un algoritmo basado en el concepto de regiones de confianza en la optimizacin

de sistemas de produccin de gas conociendo propiedades pseudo-crticas necesarias, Sin embargo
estas propiedades pseudo-criticas, Ppc y Tpc, no representan las propiedades crticas actuales de la
mezcla de gas.
Estas propiedades son usadas como parmetros de correlacin dentro de la generacin de
propiedades de gas En orden de expresar una relacin mas exacta entre las variables p,V, y T, un
factor de correccin llamado factor de compresibilidad del gas. Bsicamente, representan la
magnitud de desviacin de los gases reales que aumentan con el incremento de la presin y
temperatura a diferencia de los gases ideales adems varan ampliamente con la composicin de los
gases, es as que los gases reales tienen diferente comportamiento que los gases ideales. El problema
de optimizacin se establece considerando que un sistema de produccin de gas dado debe producir
de manera ptima la cantidad suficiente para satisfacer una demanda
El sistema de produccin de gas se describe utilizando dos unidades conceptuales: pozo y factor de
compresibilidad. As, el algoritmo propuesto realiza la simulacin del proceso para determinar las
propiedades pseudo-criticas, peso molecular total o aparente, densidad total del reservorio, gravedad
especifica as como el factor de compresibilidad del gas para una topologa previamente fijada de
pozo. Los pozos se modelan rigurosamente a ciertos mtodos haciendo varias consideraciones del
mismo lo que resulta en un sistema de ecuaciones. Las propiedades termodinmicas son calculadas

con mtodos de aproximaciones de temperatura y presin pseudo-criticas de Standing Katz y

Brown. La evidencia de los resultados obtenidos muestra la bondad del algoritmo propuesto.
III. OBJETIVOS.
1. OBJETIVO GENERAL
Determinar el factor z de los gases, mediante el mtodo de Hall Yarborough
2. OBJETIVOS ESPECFICOS
Aplicar los conceptos iteracin de programacin, para resolver el mtodo.

Aplicar el mtodo de Newton Raspn para resolver ecuaciones no lineales.
IV. MARCO TEORICO.
Los primeros estudios fisicoqumicos realizados a partir del siglo 17 involucraron el estudio de los
gases dentro de un reservorio, debido a que stos responden en forma ms dramtica a cambios en el
ambiente que los lquidos y los slidos. La mayora de estos primeros estudios estuvieron
focalizados en las variaciones de presin, temperaturas y volumen de una determinada porcin 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 ecuacin 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 ecuacin de van der Waals se diferencia de las de los gases ideales por la presencia de dos
trminos de correccin; uno corrige el volumen, el otro modifica la presin.
Los gases reales, a presiones y temperaturas cercanas a las ambientales, actan 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 nocin de las desviaciones que sufren stos bajo determinadas
condiciones de temperatura, presin 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
(ntese 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 molculas s importa, no es despreciable
- Las fuerzas de interaccin entre molculas de los gases influye
Por fortuna, las ecuaciones termodinmicas de estado para gases reales tambin pueden utilizarse
para los fluidos supercrticos. Cada especie de gas puro (agua, dixido de carbono, oxgeno, metano
, etano ,propano ,etc.) tiene un nico punto crtico, sobre el cual desaparecen los lmites entre el
lquido y el vapor. Para el punto crtico se tiene que:
La manera obvia de cambiar a un gas en un lquido es enfriarlo a una temperatura por debajo de su
punto de ebullicin. Hay otra manera de condensar un gas para formar un lquido, sin embargo,
implica aumentar la presin en el gas. Los lquidos hierven a la temperatura a la cual la presin del
vapor es igual a la presin del lquido con sus alrededores. Elevar la presin en un gas por lo tanto
aumenta con eficacia el punto de ebullicin del lquido.
Suponga que tenemos vapor de agua (o vapor) en un envase cerrado a 120 y 1 atmsfera de
presin. Puesto que la temperatura del sistema est sobre el punto de ebullicin del agua, no hay
razn para que el vapor condense para formar un lquido. Nada sucede hasta que comprimimos
lentamente el envase de tal modo que se eleva la presin del gas hasta que la presin alcanza 2
atmsferas. En este punto, el sistema est en el punto de ebullicin en el que hierve el agua y algo
del gas se condensar para formar un lquido. Tan pronto como la presin en el gas exceda 2
atmsferas, la presin del vapor del agua a 120 no es lo bastante grande para que el lquido
hierva. El gas por lo tanto condensa para formar un lquido, como se indica en la figura
.
En teora, deberamos poder predecir la presin a la cual un gas condensa a una temperatura dada
consultando un diagrama de presin del vapor contra temperatura . En la prctica, cada compuesto
tiene una temperatura crtica (Tc). Si la temperatura del gas est sobre la temperatura crtica, el gas
no se puede condensar, sin importar la presin aplicada.
La existencia de una temperatura crtica fue descubierta por Thomas Andrews en 1869 mientras que
estudiaba el efecto de la temperatura y de la presin en el comportamiento del dixido de carbono.
Andrews encontr que l podra condensar el gas del CO2 en un lquido elevando la presin en el
gas, mientras l mantuvo la temperatura por debajo de 31.0 A esta temperatura se requiere de una
presin de 72.85 atmsferas para licuefarce el gas de CO2. Andrews encontr que era imposible
convertir gas de CO en lquido por encima de esta temperatura, no importaba que tanta presin se
aplicara.
Los gases no se pueden licuefarce a temperaturas por encima de la temperatura crtica porque en
este punto las caractersticas de los gases y de los lquidos son las mismas, y no hay base sobre la
cual distinguir entre los gases y los lquidos. La presin del vapor de un lquido a la temperatura
crtica se llama la presin crtica (Pc). La presin del vapor de un lquido nunca es ms grande que
esta presin crtica.
El gas se define como un fluido homogneo 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 pequeas porciones de hexano y superiores. Los gases no hidrocarburos (o impurezas)

incluyen el dixido de carbono, sulfuro de hidrgeno y nitrgeno.
El conocimiento de las relaciones presin-volumen-temperatura (PVT) y las dems propiedades
fsicas y qumicas de los gases es esencial para resolver problemas en ingeniera de reservorios del
gas natural. Estas propiedades incluyen:
Peso molecular aparente, Ma
Gravedad especfica, g
Factor de compresibilidad, z
Densidad, g
Volumen especfico, V
Coeficiente de compresibilidad isotrmica del gas, cg
Factor volumtrico de formacin del gas, Bg
Factor de expansin del gas, Eg
Viscosidad, g
Estas propiedades del gas pueden obtenerse de medidas directas en laboratorio o por prediccin de
expresiones matemticas generalizadas. Revisaremos las leyes que describen el comportamiento
volumtrico de los gases en trminos de presin y temperatura y correlaciones matemticas que son
ampliamente usadas en la determinacin de las propiedades fsicas de los gases naturales.
Comportamiento de los Gases Ideales.
La teora cintica de los gases postula que los gases son compuestos de un gran nmero de partculas
llamadas molculas. Para un gas ideal, se asume que el volumen de estas molculas es insignificante
comparado con el volumen total ocupado por el gas. Tambin se asume que estas molculas no
tienen fuerzas atractivas o repulsivas entre ellas, y que todas las colisiones de las molculas son
perfectamente elsticas.
Basado en la teora cintica de los gases anterior, la ecuacin matemtica denominada ecuacin de
estado puede derivarse para expresar la relacin existente entre la presin p, el volumen V y la
temperatura T para una determinada cantidad de moles de gas n. Esta relacin para gases perfectos
se llama la ley de los gases ideales y se expresa matemticamente por la siguiente ecuacin:
Donde
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
Comportamiento de los Gases Reales.
Al tratar los gases a muy baja presin, la relacin de gases ideales es una herramienta conveniente y
generalmente satisfactoria. A presiones ms elevadas, el uso de la ecuacin 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 presin atmosfrica.
Bsicamente la desviacin de los gases reales de las condiciones de la ley de gas ideal se incrementa
a medida que crece la presin y temperatura y vara ampliamente con la composicin del gas. Los
gases reales se comportan de modo diferente a los gases ideales. La razn para esto es que la ley de
los gases perfectos fue derivada bajo la suposicin de que el volumen de las molculas es
insignificante y que no hay atraccin o repulsin 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
presin-volumen-temperatura para gases reales con datos experimentales. Con el fin de expresar una
relacin ms exacta entre las variables p, V y T, un factor de correccin llamado factor de
compresibilidad de los gases, factor de desviacin de los gases, o simplemente factor Z debe
introducirse en la ecuacin de estado de los gases ideales para tomar en cuenta el alejamiento de los
gases de la idealidad. La ecuacin tiene la siguiente forma:
pV= znRT
Donde el factor de compresibilidad del gas z es una cantidad adimensional y se define como la
relacin del volumen actual de n-moles de gas a T y p al volumen ideal del mismo nmero 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 mayora de
los propsitos de ingeniera cuando son expresados en trminos de las dos siguientes propiedades
adimensionales:
Presin seudo reducida

Temperatura seudo reducida
Estos trminos adimensionales se definen por las siguientes expresiones:

p
p pr
p pc
T
T pr
T pc
Donde p = presin del sistema, psia
ppr = presin seudo-reducida, adimensional
T = Temperatura del sistema, R
Tpr = Temperatura seudo-reducida, adimensional
ppc, Tpc = Presin y Temperatura seudo-crticas, respectivamente, y definidas por las
siguientes relaciones:
p pc yi pci
i 1
T pc yi Tci
i 1

Debe destacarse que estas propiedades seudo crticas, como, ppc y Tpc, no representan las
propiedades crticas actuales de la mezcla de gas. Estas seudo propiedades se usan como parmetros
de correlacin en la generacin de propiedades del gas.
Basado en el concepto de propiedades seudo reducidas, Standing y Katz (1942) presentaron un
grfico generalizado del factor de compresibilidad del gas como se muestra en la Figura 1. El
grfico representa factores de compresibilidad de gas natural dulce como funcin de ppr y Tpr. Este
grfico es generalmente confiable para gas natural con menor cantidad de no hidrocarburos. Es una
de las correlaciones ms ampliamente aceptadas en la industria de petrleo y gas.
Calculo Directo de los Factores de Compresibilidad
Despus de 4 dcadas de existencia, el grfico del factor z de Standing-Katz es an ampliamente

usado como una fuente prctica de los factores de compresibilidad del gas natural. Como resultado,
hubo una necesidad aparente para una descripcin matemtica simple de la grfica. Varias
correlaciones empricas para el clculo del factor z se han desarrollado a travs de los aos. Las
siguientes tres correlaciones empricas se describen a continuacin:
Hall-Yarborough
Dranchuk-Abu-Kassem
Dranchuk-Purvis-Robinson
V. DESARROLLO.
El Mtodo de Hall-Yarborough
Hall y Yarborough (1973) presentaron una ecuacin de estado que representa exactamente la grfica
del factor z de Standing y Katz. La expresin propuesta se basa en las ecuaciones de estado de
Starling-Carnahan. Los coeficientes de la correlacin fueron determinados ajustndolos a los datos
tomados de la grfica del factor z de Standing y Katz. Hall y Yarborough propusieron la siguiente
formula matemtica:
0.06125 p pr t 1.2(1t ) 2
z e .(1)

Donde ppr = presin seudo reducida

t = reciproca de la temperatura seudo reducida, Tpc/T
Y = densidad reducida que puede obtenerse como solucin de la siguiente ecuacin:
2 3 4
F ( ) X 1 ( X 2) 2 ( X 3) X 4 0 (2)
(1 ) 3
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

X 3 90.7 * T 242.2 * T 2 42.4 * T 3

X 4 2.18 2.82 * T
La ecuacin anterior es una ecuacin no lineal y puede ser convenientemente resuelta para la
densidad reducida Y usando la tcnica de iteracin de Newton-Raphson.
El procedimiento computacional de resolver la ecuacin principal a cualquier presin seudo
reducida especificada ppr y temperatura Tpr se resume en los siguientes pasos:
Paso 1. Hacer una estimacin inicial del parmetro desconocido, Yk, donde k es un contador de
iteracin. Un intento inicial apropiado de Y se da por la siguiente relacin:
k 0.0125 p pr te 1.2(1t ) (3)

2
Paso 2. Sustituir este valor inicial en la ecuacin (2) y evaluar la funcin no lineal. A no ser que el
valor correcto de Y se haya seleccionado inicialmente, la ecuacin (2) tendr un valor un valor
diferente de cero para F(Y):
Paso 3. Una nueva estimacin mejorada de Y, como Yk+1, se calcula de la siguiente expresin:
f ( k )
k 1 k
f ' ( k )
Donde f(Yk) se obtiene evaluando la derivada de la ecuacin 5.38 a Yk, o:
1 4 4 2 4 3 4
f ' ( ) 2( X 2) ( X 3)( X 4) ( X 41)
(1 ) 4
Paso 4. Los pasos 23 se repiten n veces, hasta que el error, abs(Yk Yk+1), llegue a ser ms
pequeo que una tolerancia establecida, como 10-12:
Paso 5. El valor correcto de Y entonces se usa para evaluar la ecuacin (1) para el factor de
compresibilidad.
Hall y Yarborough resaltaron que el mtodo no se recomienda para su aplicacin si la temperatura

seudo reducida es menor a 1.
VI. 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

VII. 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 Combo1_Click()

Nc = Val(Text3.Text)
Adodc1.Refresh
Adodc1.Recordset.MoveFirst
For i = 1 To 20
If Not (Adodc1.Recordset.Fields("COMPONENTES") = Combo1.Text) Then
If Adodc1.Recordset.EOF Then
Else
Adodc1.Recordset.MoveNext
End If
Else
f=f+1
y(f) = Val(InputBox("Fraccion molar en porcentaje % y(" & Combo1.Text & ")="))
y(f) = y(f) / 100
s = s + y(f)
With Adodc1.Recordset
MSFlexGrid1.TextMatrix(f, 0) = f
MSFlexGrid1.TextMatrix(f, 1) = .Fields("COMPONENTES")
MSFlexGrid1.TextMatrix(f, 2) = Format(.Fields("PESO MOLECULAR"), "##,####0.0000")

pm(f) = .Fields("PESO MOLECULAR")

MSFlexGrid1.TextMatrix(f, 3) = Format(y(f), "##,####0.0000")
pc(f) = .Fields("PRESION CRITICA")
MSFlexGrid1.TextMatrix(f, 4) = Format(.Fields("PRESION CRITICA"), "##,####0.0000")
tc(f) = .Fields("TEMPERATURA CRITICA")
MSFlexGrid1.TextMatrix(f,5) = Format(.Fields("TEMPERATURA CRITICA"),"##,####0.0000")
End With
MSFlexGrid1.TextMatrix(f + 1, 3) = s
MSFlexGrid1.TextMatrix(f + 1, 2) = "Suma total"
Combo1.RemoveItem Combo1.ListIndex
If f = Nc Then
Combo1.Visible = False
Exit For
Command1.SetFocus
Else
Combo1.SetFocus
Exit For
End If
End If
Next i
End Sub
Private Sub Command1_Click()

Nc = Val(Text3.Text)
pcr = 0
tcr = 0
If s = 1 Then
For i = 1 To Nc
pcr = pcr + pc(i) * y(i)
Next i
Text4.Text = Format(pcr, "##,####0.0000")
Else
MsgBox "LA SUMA DE LAS FRACCIONES MOLARES NO SUMAN 1"
End If
If s = 1 Then
For i = 1 To Nc
tcr = tcr + tc(i) * y(i)
Next i
Text5.Text = Format(tcr, "##,####0.0000")
Else

MsgBox "LA SUMA DE LAS FRACCIONES MOLARES NO SUMAN 1"
End If
P = Val(Text1.Text)
T = Val(Text2.Text)
Psr = P / pcr
Tsr = T / tcr
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
Private Sub Command2_Click()

End
End Sub
Private Sub Form_Load()
Skin1.LoadSkin App.Path & "/GREEN.SKN"

Skin1.ApplySkin Me.hWnd
End Sub
VIII. DATOS DE PRUEBA.
El gas natural boliviano tiene la siguiente composicin, para una presin de 3800 psia y un
temperatura de 670 R, determinar z.
Componente % Molar
C1 87.832
C2 7.817
C3 1.357
nC4 0.079
iC4 0.055
nC5 0.016
iC5 0.015
N2 1.347
CO2 1.482
Sumatoria 100
IX. RESULTADOS.

Arrancamos el programa.
Introducimos datos.
Calculo de z.

X. CONCLUSIONES.
Se utilizaron satisfactoriamente las ecuaciones de iteracin de Newton Raphson., para

determinar la densidad reducida.
Este mtodo tiene sus limitaciones, para Tsr menores a 1, el mtodo arroja datos erroneos.
El diseo del algoritmo y la aplicacin de este fue todo un xito ya que se logro el proposito
esperado adems se obtuvo parmetros bsicos necesarios en campo.
La ejecucin del programa con un problema prctico de pozo gasifero demostr la veracidad del
funcionamiento del programa.
Se determino satisfactoriamente el valor de z, talvez con algunos errores muy pequeos, ya que
los valores de las presiones crticas y temperaturas crticas varan un poco.
XI. BIBLIOGRAFA.
Ingeniera de Reservorios Ing. Rolando Camargo Gallegos

Ingeniera de Reservorios I Ing. Hermas Herrera Callejas
BRADY.(1993) Qumica bsica. V ed. Noruega editores. Mjico.
http://www.tannerm.com/physical/vanderWaals/vander.htm
GLASTONE. Fundamentos de Fsicoqumica
XII. ANEXOS.

Figura 1 Grfico de factores de compresibilidad de Standing y Katz. (Cortesa de GPSA y GPA Engineering Data Book,
EO Edicin, 1987.)

Calculo Del Factor Z de Los Gases Por El

Cargado por

Copyright:

Formatos disponibles

Calculo Del Factor Z de Los Gases Por El

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 Del Factor Z de Los Gases Por El

Cargado por

Copyright:

Formatos disponibles

Universidad Mayor de San Andrs Programacin Aplicada

Ingeniera Petrolera PET 230

CALCULO DEL FACTOR Z DE LOS GASES POR EL METODO

El factor de compresibilidad Z, es un factor de correccin, que se introduce en la ecuacin de estado

Introduciendo el factor de correccin Z:

El factor Z tambin se puede entender como:

Significado del valor de Z.

Univ. Henry Isaac Cruz Lucas 1

El factor de compresibilidad, z, de un gas natural es una medida de su desviacin respecto al

Se describe y aplica un algoritmo basado en el concepto de regiones de confianza en la optimizacin

Univ. Henry Isaac Cruz Lucas 2

con mtodos de aproximaciones de temperatura y presin pseudo-criticas de Standing Katz y

Determinar el factor z de los gases, mediante el mtodo de Hall Yarborough

Aplicar los conceptos iteracin de programacin, para resolver el mtodo.

IV. MARCO TEORICO.

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

Univ. Henry Isaac Cruz Lucas 3

Univ. Henry Isaac Cruz Lucas 4

Univ. Henry Isaac Cruz Lucas 6

pentano y pequeas porciones de hexano y superiores. Los gases no hidrocarburos (o impurezas)

Comportamiento de los Gases Ideales.

Comportamiento de los Gases Reales.

Presin seudo reducida

Estos trminos adimensionales se definen por las siguientes expresiones:

Univ. Henry Isaac Cruz Lucas 8

Calculo Directo de los Factores de Compresibilidad

Despus de 4 dcadas de existencia, el grfico del factor z de Standing-Katz es an ampliamente

Donde ppr = presin seudo reducida

Univ. Henry Isaac Cruz Lucas 9

X 3 90.7 * T 242.2 * T 2 42.4 * T 3

k 0.0125 p pr te 1.2(1t ) (3)

Donde f(Yk) se obtiene evaluando la derivada de la ecuacin 5.38 a Yk, o:

Hall y Yarborough resaltaron que el mtodo no se recomienda para su aplicacin si la temperatura

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

Univ. Henry Isaac Cruz Lucas 10

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. Henry Isaac Cruz Lucas 11

VII. EL PROGRAMA EN CODIGO VISUAL BASIC.

Private Sub Combo1_Click()

Univ. Henry Isaac Cruz Lucas 12

pm(f) = .Fields("PESO MOLECULAR")

Private Sub Command1_Click()

Univ. Henry Isaac Cruz Lucas 13

MsgBox "LA SUMA DE LAS FRACCIONES MOLARES NO SUMAN 1"

Text6.Text = Format(Psr, "##,####0.0000")

If Tsr > 1 Then

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

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

Univ. Henry Isaac Cruz Lucas 14

MsgBox "Error, corrija datos"

Private Sub Command2_Click()

Private Sub Form_Load()