YACIMIENTOS II
PROPIEDADES DE LOS FLUIDOS
Docente: MSC. Ing. Adriangela Chiquinquirá Romero Sánchez
Especialista en Geología Petrolera / MSC. En Ing. De Gsas / MSC. En Gestión Ambiental
Email:
[email protected]PROPIEDADES DE LOS FLUIDOS
DE YACIMIENTO
Para evaluar un yacimiento
Petróleo en términos del
Agua rendimiento de
Gas natural hidrocarburos
Propiedades
físicas
En fondo del
En superficie
pozo
ANÁLISIS PVT
Los análisis PVT son un conjunto de pruebas que se hacen
en el laboratorio para determinar las propiedades de los
fluidos de un yacimiento y su variación con la presión, con la
finalidad de conocer el comportamiento de los mismos.
Los análisis PVT simulan en el laboratorio el agotamiento de
presión (depleción) de un yacimiento volumétrico e
isométrico midiendo exactamente los volúmenes de gas y
liquido separados en cada decremento de presión.
OBJETIVOS DE LOS ANÁLISIS PVT
Entre los objetivos de los análisis PVT, se encuentran:
1. Determinar ciertos parámetros del yacimiento y
condiciones actuales del pozo, para predecir el futuro
comportamiento del mismo.
2. Permite estimar nuevas reservas y diseñar la terminación
o completamiento del pozo más adecuado a las necesidades
del mismo.
3. Permite saber si existe suficientes hidrocarburos que
justifiquen la perforación de nuevos pozos o desarrollo de
nuevos campos.
4. Determinar la relación volumétrica de las fases desde
condiciones de yacimiento hasta las condiciones de
superficie, y también determinar la solubilidad del gas en el
crudo.
OBJETIVOS DE LOS ANÁLISIS PVT
Las principales propiedades que pueden ser obtenidas en
un estudio PVT son:
• Factor Volumétrico de formación del petróleo, Bo
• Relación Gas en Solución, Rs
• Coeficiente de compresibilidad isotérmica, Co
• Densidad del gas y del aceite, g y o
• Factor Volumétrico de formación del gas, Bg
• Factor de compresibilidad del gas, Z
• Viscosidad del gas y del aceite, µg y µo
• Composición del gas y del petróleo
DENSIDAD
La densidad se define como la cantidad de masa de fluido
por unidad de volumen a una temperatura y una presión
especifica.
Es una propiedad importante afecta los valores para el
calculo de volumen. Es necesario calcular la densidad del
fluido a condiciones de yacimiento ya que después de ser
sometido a cambios de presión y temperatura desde el
yacimiento a superficie es evidente que se disminuya el
volumen del fluido y se tome un valor de densidad incorrecto
.
𝑀 𝑃𝑀𝑎 UNIDADES:
ρg: es la densidad del fluido en lb/pie3
𝜌𝑜 = 𝜌𝑔 = Ma: es el peso molecular aparente lb/lbmol
𝑉 𝑍𝑅𝑇T: es la temperatura en °R
R: constante 10.73 lpca*ft^3/lbmol*ºR
P: presión en lpca.
GRAVEDAD ESPECÍFICA
La gravedad especifica se define como la razón de la
densidad del fluido a una densidad estándar ( En el
caso de ser aceite se toma como estándar la del
agua y si es gas se tomara la del aire). “Ambas
medidas a la misma temperatura y presión”.
“Es una propiedad adimensional”.
𝜌𝑜 𝜌𝑔
𝑦𝑜 = 𝑦𝑔 =
𝜌𝑎𝑔𝑢𝑎 𝜌𝑎𝑖𝑟𝑒
GRAVEDAD API
Es una medida de densidad , que en comparación con el agua,
precisa cuan pesado o liviano esta el petróleo
Extra liviano Liviano Mediano
>40 30-39.9 22-29.9
Pesado Extra pesado
10-21.9 <10
141.5
𝑦𝐴𝑃𝐼 = − 131.5
𝑦𝑜
PRESIÓN DE BURBUJA
La presión de Burbuja o también
llamada presión de Saturación (
Pb), se define como la presión a
la cual se forma la primera
burbuja de gas, al pasar de un
sistema de estado liquido a un
estado de dos fases (liquido-
vapor) e donde la fase liquida se
encuentra en equilibrio con una
cantidad mínima de gas libre.
RELACIÓN DE GAS DISUELTO EN
PETRÓLEO
La cantidad de gas disuelto en el
petróleo a condiciones de
yacimientos es llamado relación
de gas disuelto en el petróleo;
es decir, es la cantidad de gas
que se libera del petróleo a
medida que el mismo es
transportado desde el yacimiento
a la superficie. Esta relación es
definida en términos de cantidad
de gas y petróleo que llega a la
superficie durante su producción.
1 cm³ = 3.53E-5 ft³
1 cm³ = 6.29E-6 barriles
RELACIÓN DE GAS DISUELTO EN
PETRÓLEO
La relación de gas disuelto en el
petróleo (Rs), se define como la
cantidad de pies cúbicos normales
(PCN) de gas que puede disolverse en
un barril normal (BN) de petróleo ,
cuando ambos son llevados a las
condiciones de presión y temperatura
del yacimiento.
𝑃𝐶𝑁 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑔𝑎𝑠 𝑒𝑛 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛
𝑅𝑠 = =
𝐵𝑁 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑎𝑐𝑒𝑖𝑡𝑒
Típica gráfica de la solubilidad del gas en función de
presión a temperatura de yacimiento constante
Relación gas en solución – petróleo
Constante
hasta alcanzar
la Pb
Menos gas
permanece en
el líquido
Presión de yacimiento
En un yacimiento, la solubilidad del gas a presiones mayores a la presión de
burbujeo (𝑃𝑏 ) es constante e igual a 𝑅𝑠𝑖 , ya que a estas presiones no se desarrolla
gas en el espacio poroso y la mezcla como líquido se desplaza hacia el pozo.
Cuando la presión se reduce por debajo de la presión de burbujeo, se comienza a
formar gas en el yacimiento y menos gas permanece en el líquido. La diferencia
entre la solubilidad en el punto de burbujeo y la solubilidad a cualquier presión
por debajo del punto de burbujeo representa la cantidad de gas que ha sido
liberada en el espacio poroso.
FACTOR VOLUMÉTRICO DEL
PETRÓLEO
El volumen del petróleo que entra al tanque en superficie es menor al
volumen de petróleo que sale del yacimiento.
• Se define como el volumen que ocupa a condiciones de
yacimiento un barril normal de petróleo mas su gas en solución.
• También puede definirse como el cambio en el volumen que
experimenta la fase liquida al cambiar de condiciones de
yacimiento a condiciones en superficie como consecuencia de la
expansión liquida o liberación del gas en solución.
FACTOR VOLUMÉTRICO DEL
PETRÓLEO
𝑉𝑜 𝑃,𝑇
𝐵𝑜 =
𝑉𝑜 𝐶𝑁
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑝𝑒𝑡𝑟ó𝑙𝑒𝑜 + 𝑔𝑎𝑠 𝑑𝑖𝑠𝑢𝑒𝑙𝑡𝑜 𝑒𝑛 𝑒𝑙 𝑦𝑎𝑐𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑎 𝑐𝑜𝑛𝑑𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑑𝑒 𝑦𝑎𝑐𝑖𝑚𝑖𝑒𝑛𝑡𝑜
𝐵𝑜 =
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑝𝑒𝑡𝑟ó𝑙𝑒𝑜 𝑒𝑛 𝑒𝑙 𝑡𝑎𝑛𝑞𝑢𝑒 𝑑𝑒 𝑎𝑙𝑚𝑎𝑐𝑒𝑛𝑎𝑚𝑖𝑒𝑛𝑡𝑜 𝑎 𝑐𝑜𝑛𝑑𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑛𝑜𝑟𝑚𝑎𝑙𝑒𝑠
Donde 𝐵𝑜 es el factor volumétrico del petróleo en BY/BN, 𝑉𝑜 𝑃,𝑇 es el volumen de
petróleo más el gas disuelto en condiciones de presión y temperatura del
yacimiento en BY, 𝑉𝑜 𝐶𝑁 el volumen de petróleo medido en condiciones normales
BN.
FACTOR VOLUMÉTRICO DEL
PETRÓLEO
Líquido se
expande en el
yacimiento
El volumen de
líquido
disminuye por la
liberación de
gas
Factor Volumétrico del petróleo en función de presión
FACTOR VOLUMÉTRICO TOTAL
Es el volumen total del sistema dividido entre el volumen de petróleo a
condiciones de tanque
𝑉𝑡
𝐵𝑡 =
𝑉𝐶𝑁
Donde
𝑉𝑡 = 𝑉𝑜 + 𝑉𝑔
𝐵𝑡 también puede definirse como:
𝐵𝑡 = 𝐵𝑜 + 𝐵𝑔 𝑅𝑠𝑖 − 𝑅𝑠 𝐵𝑔 (BY/PCN) 𝐵𝑜 (BY/BN)
𝑅𝑠 (PCN/BN)
𝐵𝑜 y 𝐵𝑔 son los factores volumétricos del petróleo y gas respectivamente.
𝑅𝑠𝑖 (también 𝑅𝑠𝑏 ) es la relación gas disuelto en el petróleo inicial de
yacimiento. 𝑅𝑠 es la relación gas disuelto en el petróleo a la temperatura
del yacimientos y a la presión requerida.
Al comparar el factor volumétrico del petróleo con el factor volumétrico
total. Los dos valores son idénticos por encima de la presión en el punto
de burbuja ya que el gas no se ha desprendido en el yacimiento a estas
presiones. La diferencia entre las dos cantidades debajo de la presión en el
punto de burbuja refleja la emisión de gas en la formación. A medida que la
presión del yacimiento disminuye el volumen de gas liberado en la
formación incrementa y los dos factores volumétricos difieren.
Por encima de la presión de burbuja 𝑅𝑠𝑏 = 𝑅𝑠
𝐵𝑜 = 𝐵𝑡
VISCOSIDAD DEL PETRÓLEO
La viscosidad del petróleo, al igual que la viscosidad del
gas, es una medida de la resistencia al flujo.
Esta propiedad depende fuertemente de la temperatura
del yacimiento, la presión, la gravedad del petróleo, la
gravedad y solubilidad del gas.
VISCOSIDAD DEL PETRÓLEO
Generalmente la viscosidad del
líquido incrementa con el incremento
de la presión encima del punto de
burbuja. Por debajo de la presión de
burbuja, ambos gas y líquido están
presentes y un incremento en la
presión origina que los componentes
gaseosos ligeros presentes entren en
la fase líquida, disminuyendo la
viscosidad del líquido. Así, a mayor
gas disuelto en el petróleo del
yacimiento, menor es la viscosidad del Relación de la viscosidad del
petróleo en función de presión a
crudo y más fácilmente fluirá éste en temperatura constante
el yacimiento.
VISCOSIDAD DEL PETRÓLEO
De acuerdo con la presión la viscosidad del crudo se clasifica en:
• Viscosidad del petróleo muerto (μod): viscosidad a la presión
atmosférica (sin gas disuelto) y a la temperatura del yacimiento.
• Viscosidad del petróleo saturado (μob): viscosidad a la presión de
burbujeo y a la temperatura del yacimiento.
• Viscosidad del petróleo no saturado (μo): viscosidad a una presión por
encima del punto de burbujeo y a la temperatura de yacimiento.
COMPRESIBILIDAD ISOTÉRMICA
El coeficiente de compresibilidad es definido como el cambio
en el volumen debido al incremento de la presión por unidad
de volumen a una temperatura constante.
Presión mayor o igual al presión de burbuja (petróleo en una
sola fase con todo el gas disuelto en solución,
compresibilidad isotérmica subsaturada).Presión por debajo
de la presión de burbuja (el gas es liberado , compresibilidad
isotérmica saturada).
1 𝑑𝑣
𝑐= T
𝑉 𝑑𝑝
COMPRESIBILIDAD ISOTÉRMICA
Presión mayor o igual al presión de burbuja
Donde
Co :Compresibilidad del petróleo
Vo: Volumen de petróleo
Bo: Factor Volumétrico de Formación del petróleo
ρo : Densidad del petróleo
Presión por debajo de la presión de burbuja
Donde
Co :Compresibilidad del petróleo
Bg: Factor volumétrico de formación del gas
Bo: Factor Volumétrico de Formación del petróleo
Rs : relación gas disuelto
PROPIEDADES DE LOS FLUIDOS
DE YACIMIENTO
La información de las propiedades
Pruebas experimentales que se
realizan a muestras tomadas del
yacimiento.
físicas. Determinarlas en forma
aproximada aplicando
correlaciones derivadas
empíricamente a partir de datos
de campo en yacimientos con
fluidos similares.
PRESIÓN DE BURBUJA
La presión del punto de burbuja, se
determina en función de:
➢ Temperatura (T).
➢ Gravedad especifica del gas (yg).
➢ Gravedad API.
➢ Gas disuelto (Rs).
Las correlaciones que permiten determinar el punto de
burbujeo en función de estas variables son:
CORRELACIÓN DE STANDING
Standing Desarrollo una expresión matemática para
calcular el 𝑃𝑏 en función de la relación gas en solución
petróleo, gravedad especifica del gas, gravedad API y la
temperatura del yacimiento. 0.83
𝑃𝑏 = 18,2 ∗ 𝐹 − 1,4 𝑅𝑠𝑏 𝑎
𝐹= 10
𝛾𝑔
𝑎 = 0.00091 𝑇 − 0,0125 𝐴𝑃𝐼
EJEMPLO
Con los siguientes datos determine la presión de burbujeo para
las muestras, haciendo uso de la correlación de Standing:
Muestra T(ºF) Rs ºAPI Yg
s
1 250 751 47,1 0,851
2 220 768 40,7 0,855
Para calcular (Pb), hallamos el coeficiente a, y reemplazamos
en la ecuación así:
EJEMPLO
MUESTRA 1
𝑎1 = 0.00091 𝑇 − 0,0125 𝐴𝑃𝐼
𝑎1 = 0.00091 250 − 0,0125 47,1 = −0,36125
0.83
𝑅𝑠𝑏 𝑎
𝐹= 10
𝛾𝑔
0.83
751 −0,36125
𝐹= 10 = 121,25
0,851
𝑃𝑏 = 18,2 ∗ 𝐹 − 1,4
𝑃𝑏 = 18,2 ∗ 121,25 − 1,4 = 2181,31 𝑙𝑝𝑐
EJEMPLO
MUESTRA 2
𝑎1 = 0.00091 𝑇 − 0,0125 𝐴𝑃𝐼
𝑎1 = 0.00091 220 − 0,0125 40,7 = −0,30855
0.83
𝑅𝑠𝑏 𝑎
𝐹= 10
𝛾𝑔
0.83
768 −0,30855
𝐹= 10 = 138,92
0,855
𝑃𝑏 = 18,2 ∗ 𝐹 − 1,4
𝑃𝑏 = 18,2 ∗ 138,92 − 1,4 = 2502,88 𝑙𝑝𝑐
EJERCICIO
Con los siguientes datos, hallar la presión de burbujeo usando la
correlación de Vásquez y Beggs:
Muestra Rsb yg T °F ºAPI Tsep °F Psep
1 943 0,781 218 44,2 60 200
2 807 0,848 180 27,3 173 85
Para hallar Pb, se halla el coeficiente a y la gravedad del gas en
el separador y posteriormente se reemplaza en la ecuación.
CORRELACIÓN DE GLASO
1.7669+1.7447𝑙𝑜𝑔𝐹−0.30218 𝑙𝑜𝑔𝐹 2
𝑃𝑏 = 10
0.816
𝑅𝑠𝑏 𝑇 0.172
F= ∗
𝛾𝑔 𝐴𝑃𝐼 0.989
CORRELACIÓN DE GLASO
Adicionalmente, Glaso desarrolló las siguientes ecuaciones para corregir la
presión de burbujeo por efectos de gases no-hidrocarburos: 𝑁2 , 𝐶𝑂2 y/o 𝐻2 𝑆
Donde T en °F y 𝑦𝑁2 , 𝑦𝐶𝑂2 , 𝑦𝐻2𝑆 las fracciones molares de 𝑁2 , 𝐶𝑂2 y/o 𝐻2 𝑆 en los
gases de superficie
Jacobsen desarrolló la siguiente ecuación para corregir la presión de burbujeo
por efecto del 𝑁2
Donde T en °F y 𝑦𝑁2 , la fracción molar de 𝑁2 en los gases de superficie
CORRELACIÓN DE GLASO
Estas correcciones son simplemente factores de multiplicación aplicados a la
presión de burbujeo calculada, por lo tanto, cualquier factor de corrección puede
ser utilizado con cualquier correlación de 𝑃𝑏 , esto es:
De acuerdo a Ostermann y Cols., estas correcciones son necesarias, despreciar
su uso puede resultar en errores del 30% o más en la estimación de la presión
de burbujeo
𝑃𝑏𝑐 = 𝑃𝑏 ∗ 𝐶𝑁2 ∗ 𝐶𝐶𝑂2 ∗ 𝐶𝐻2𝑆
Para crudos volátiles, Glaso recomienda que el exponente de la temperatura se
cambie a 0.130.
EJERCICIO
Usando los siguientes datos y la correlación de Glaso, halle la
presión de burbujeo:
Muestra Rsb T ºAPI Yg
1 693 260 48,6 0,911
2 968 237 40,5 0,898
0,816
693 2600,172
𝐹1 =
0,911 48,60,989
F1= 12,54
𝑃𝑏 = 10 1.7669+1.7447 𝑙𝑜𝑔 12,54−0.30218𝑙𝑜𝑔 12,54 2
Pb1= 2082,53 lpca
0,816
968 2370,172
𝐹2 =
0,898 40,50,989
F2= 19,6
𝑃𝑏 = 10 1.7669+1.7447 𝑙𝑜𝑔 19,6−0.30218𝑙𝑜𝑔 19,6 2
Pb2= 3287,6 lpca
CORRELACIÓN DE LASATER
Para el cálculo de 𝑦𝑔 , se requiere conocer la razón gas disuelto −
petróleo a 𝑝 ≥ 𝑝𝑏 , 𝑅𝑠𝑏, la gravedad específica del petróleo, 𝛾�𝑜, y el peso
molecular efectivo del petróleo, 𝑀�𝑜 a condiciones normales. La
siguiente ecuación permite determinar 𝑦𝑔
𝑅𝑠𝑏
379,3
𝑦′𝑔∗ =
𝑅𝑠𝑏 350𝛾𝑜
+
379,3 𝑀𝑜
CORRELACIÓN DE LASATER
Para 𝛾′𝑔 ≤ 0,60 𝑃𝑏 𝛾𝑔
= 0,679𝑒𝑥𝑝 2,786𝛾′𝑔 − 0,323
𝑇
𝑃𝑏 𝛾𝑔 3,56
Para 𝛾′𝑔 > 0,60 = 8,26𝛾′𝑔 + 1,95
𝑇
Donde T en °R
Para 𝛾𝐴𝑃𝐼 ≤ 40 𝑀𝑜 = 630 − 10𝛾𝐴𝑃𝐼
Para 𝛾𝐴𝑃𝐼 > 40 𝑀𝑜 = 73110𝛾𝐴𝑃𝐼 −1,562
CORRELACIÓN DE AL-MARHOUN
0.715082 −1.87784 −3.1437 1.32657
𝑃𝑏 = 5.38088𝑥10−3 𝑅𝑠𝑏 𝛾𝑔 𝛾𝑜 𝑇
CORRELACIÓN DE
KARTOATMODJO, T. y SCHMIDT, Z.
𝑐4
𝑅𝑠𝑏
𝑃𝑏 = 𝐶3 𝛾𝐴𝑃𝐼
𝐶2
𝐶1 𝛾𝑔𝑐 10 𝑇+460
CORRELACIÓN DE
KARTOATMODJO, T. y SCHMIDT, Z.
CORRELACIÓN DE
KARTOATMODJO, T. y SCHMIDT, Z.
Al igual que Vázquez y Beggs, Kartoatmodjo y Schmidt desarrollaron la
siguiente correlación para corregir la gravedad específica del gas a una
presión de separación de 100 lpcm (114,7 lpca)
CORRELACIÓN DE
KARTOATMODJO, T. y SCHMIDT, Z.
CORRELACIÓN DE PETROSKY Y
FARSHAD
𝑃𝑏 = 112.727 𝐹 − 12.34
𝑅𝑠𝑏 0.5774 (4.561𝑥10−5𝑇 1.3911 −7.916𝑥10−4𝛾 1.5410 )
𝐹= 0.8439
10 𝐴𝑃𝐼
𝛾𝑔
GAS DISUELTO EN EL PETRÓLEO
Gas disuelto
El gas disuelto Rs, se calcula en función:
Presión
Temperatura
Gravedad API
Gravedad especifica
Gravedad especifica del gas
GAS DISUELTO EN EL PETRÓLEO
Los siguientes factores afectan la solubilidad del gas en el
petróleo:
• Presión: Al aumentar la presión aumenta Rs.
• Temperatura: Al aumentar la temperatura disminuye Rs.
• Gravedad del Petróleo: Al aumentar la gravedad API
aumenta Rs.
• Gravedad del gas: Al aumentar la gravedad especifica del
gas aumenta Rs.
CORRELACION STANDING
Standing propuso una correlación para determinar la solubilidad
del gas en función:
Presión
Gravedad API
Gravedad especifica del gas
Temperatura del sistema
La correlación fue desarrollada de un total de 105 puntos de datos
determinados de 22 mezclas de hidrocarburo y gas natural , de
crudo de California, esta correlación presenta un error de 10%
respecto a los valores determinados experimentalmente.
CORRELACIÓN DE STANDING
Correlación de Standing: Desarrollo una expresión matemática para calcular el 𝑅𝑠
en función de la presión, gravedad especifica, gravedad API y la temperatura del
sistema.
1,2048
𝑃
𝑅𝑠 = 𝛾𝑔 + 1,4 10𝑥
Donde
18,2
𝑥 = 0,0125𝐴𝑃𝐼 − 0,00091 𝑇
Donde
CORRELACIÓN DE STANDING
CORRELACIÓN VÁZQUEZ Y BEGGS
𝐴𝑃𝐼
0.0362𝛾𝑔𝑐 𝑃1.0937 𝑒𝑥𝑝 25.724 ; 𝐴𝑃𝐼 ≤ 30
𝑇 + 460
𝑅𝑠 =
𝐴𝑃𝐼
0.0178𝛾𝑔𝑐 𝑃1.187 𝑒𝑥𝑝 23.931 ; 𝐴𝑃𝐼 > 30
𝑇 + 460
𝑃𝑠𝑝
𝛾𝑔𝑐 = 𝛾𝑔 1 + 5.912𝑥10−5 𝐴𝑃𝐼𝑥𝑇𝑠𝑝 𝑙𝑜𝑔
114.7
FACTOR VOLUMÉTRICO DEL
PETRÓLEO
Correlación de Standing: Desarrollo una expresión matemática para calcular el 𝐵𝑜
en función de la relación gas en solución petróleo, gravedad especifica del gas,
gravedad específica del petróleo y la temperatura del yacimiento.
0,5 1,2
𝛾𝑔
𝐵𝑜 = 0,9759 + 0,000120 𝑅𝑠 + 1,25 𝑇
𝛾𝑜
FACTOR VOLUMÉTRICO DEL
PETRÓLEO
CORRELACIÓN VÁZQUEZ Y BEGGS
𝛾𝐴𝑃𝐼 𝛾𝐴𝑃𝐼
𝐵𝑜 = 1.0 + 𝐶1 𝑅𝑠𝑏 + 𝐶2 𝑇 − 60 + 𝐶3 𝑅𝑠𝑏 𝑇 − 60
𝛾𝑔𝑐 𝛾𝑔𝑐
Las correlaciones anteriores pueden ser utilizadas para estimar valores de 𝐵𝑜
a presiones por debajo de la presión de burbujeo, en este caso, se debe
utilizar 𝑅𝑠 a la presión deseada en vez de 𝑅𝑠𝑏
Los valores de 𝐵𝑜 puede determinarse a presiones mayores que la presión
del punto de burbujeo 𝑝 > 𝑝𝑏 mediante la siguiente ecuación:
𝐵𝑜 = 𝐵𝑜𝑏 𝑒𝑥𝑝 𝐶𝑜 𝑃𝑏 − 𝑃
FACTOR VOLUMÉTRICO TOTAL
CORRELACIÓN DE GLASO
CORRELACIÓN GLASO
8,0135𝑥10−2 +4,7257𝑥10−1 𝑙𝑜𝑔𝐹+1,7351𝑥10−1 𝑙𝑜𝑔𝐹 2
𝐵𝑡 = 10
𝑇 0,5 −1,1089 𝛾 2,9𝑥10−0,00027𝑅𝑠
𝐹 = 𝑅𝑠 𝑃 𝑜
𝛾𝑔 0,3
FACTOR VOLUMÉTRICO TOTAL
CORRELACIÓN DE AL-MARHOUN
𝐵𝑡 = 0,314693 + 0,106253𝑥10−4 𝐹 + 0,18883x10−10 𝐹 2
𝐹 = 𝑅𝑠 0.644516 𝛾𝑔 −1,07934 𝛾𝑜 0,724874 𝑃−0,76191 𝑇 2,00621
COMPRESIBILIDAD DEL
PETRÓLEO
CORRELACIÓN VÁZQUEZ Y BEGGS
Las siguientes correlaciones permiten determinar 𝐶𝑜 a presiones mayores que la
presión del punto de burbujeo 𝑝 > 𝑝𝑏
−1433 + 5𝑅𝑠𝑏 + 17,2 𝑇 − 1180𝛾𝑔𝑐 + 12,61 °𝐴𝑃𝐼
𝐶𝑜 =
𝑃𝑥105
𝑃𝑠𝑒𝑝
𝛾𝑔𝑐 = 𝛾𝑔 1 + 5,912 10−5 𝐴𝑃𝐼 𝑇𝑠𝑒𝑝 − 460 𝑙𝑜𝑔
114,7
Donde 𝛾𝑔𝑐 es la gravedad del gas a una presión de referencia del separador
(100lpcm), 𝛾𝑔 es la gravedad del gas en las condiciones actuales del separador, es
decir 𝑃𝑠𝑒𝑝 y 𝑇𝑠𝑒𝑝 , donde 𝑃𝑠𝑒𝑝 es la presión actual del separador en lpca y 𝑇𝑠𝑒𝑝 es la
temperatura actual del separador en °𝑅
CORRELACIÓN VÁZQUEZ Y BEGGS
CORRELACIÓN DE PETROSKY Y
FARSHAD
𝐶𝑜 = 1,705𝑥10−7 𝑅𝑠0.69357 𝛾𝑔0,1885 𝛾𝐴𝑃𝐼
0,3272 0,6729 −0,5906
𝑇 𝑃
COMPRESIBILIDAD DEL
PETRÓLEO
Cuando la presión del yacimiento esta por debajo de la presión en el punto de
burbuja, el volumen de líquido en el yacimiento decrece cuando la presión se
reduce. El cambio fraccional de volumen con respecto al cambio en la presión puede
ser representado por:
1 𝜕𝐵𝑜 𝜕𝑅𝑠
𝐶𝑜 = − − 𝐵𝑔
𝐵𝑜 𝜕𝑃 𝑇
𝜕𝑃 𝑇
Curva típica de Co en función de presión a Curva típica de Co en función de presión a
temperatura constante. (Presión por encima del temperatura constante. (Presión por debajo del
punto de burbuja) punto de burbuja)
COMPRESIBILIDAD DEL
PETRÓLEO
CORRELACIÓN DE McCAIN
La siguiente correlación de McCain, permite determinar coeficiente de compresibilidad a
presiones menores que la presión de burbuja:
𝑙𝑛𝐶𝑜 = −7,573 − 1,450ln(𝑃) − 0,383ln(𝑃𝑏 ) + 1,402ln(𝑇) + 0,256ln(𝐴𝑃𝐼) + 0,449ln(𝑅𝑠𝑏 )
Si no se dispone de la presión de burbuja, la siguiente correlación puede ser
utilizada:
𝑙𝑛𝐶𝑜 = −7,663 − 1,497ln(𝑃) − 1,115ln(𝑇) + 0,533ln(𝐴𝑃𝐼) + 0,184ln(𝑅𝑠𝑏 )
Si no se dispone de la presión de burbuja y de 𝑅𝑠𝑏 la siguiente correlación puede ser
utilizada:
𝑙𝑛𝐶𝑜 = −7,114 − 1,394ln(𝑃) − 0,98ln(𝑇) + 0,770ln(𝐴𝑃𝐼) + 0,446ln(𝛾𝑔 )
Donde 𝐶𝑜 el 𝑙𝑝𝑐 −1 , Presión en lpca y Temperatura en °𝑅
CORRELACIÓN DE McCAIN
COMPRESIBILIDAD DEL
PETRÓLEO
CORRELACIÓN DE STANDING
La siguiente correlación de Standing permite determinar coeficiente de
compresibilidad a presiones menores que la presión de burbuja:
1 𝜕𝑅𝑠 𝜕𝐵𝑜
𝐶𝑜 = − − 𝐵𝑔
𝐵𝑜 𝜕𝑃 𝑇
𝜕𝑅𝑠 𝑇
0.5 0.5 0,2
𝜕𝐵𝑜 𝛾𝑔 𝛾𝑔
= 14,4 ∗ 10−5 𝑅𝑠 + 1,25𝑇
𝜕𝑅𝑠 𝑇
𝛾𝑜 𝛾𝑜
𝜕𝑅𝑠 𝑅𝑠
=
𝜕𝑃 𝑇
0,83𝑃 + 21,15
CORRELACIÓN VASQUEZ-BEGGS
𝜕𝐵𝑜 𝛾𝐴𝑃𝐼
= 𝐶1 + 𝐶3 𝑇 − 60
𝜕𝑅𝑠 𝑇 𝛾𝑔𝑐
𝜕𝑅𝑠 𝑅𝑠
= 𝐶2
𝜕𝑃 𝑇
𝑃
CORRELACIÓN GLASO
0,526
𝜕𝐵𝑜 2 1 𝛾𝑔
= 2,91329 − 0,55366𝑙𝑜𝑔𝐹 10 −6,58511+2,91329𝑙𝑜𝑔𝐹−0,27683 𝑙𝑜𝑔𝐹
𝜕𝑅𝑠 𝑇
𝐹 𝛾𝑜
𝜕𝑅𝑠 𝑅𝑠
= 2,02777 0,5
𝜕𝑃 𝑇
𝑃 14,1811 − 3,3093𝑙𝑜𝑔𝑃
CORRELACIÓN AL-MARHOUN
𝜕𝐵𝑜 1,35556𝑥10−3 𝐹 + 4,723𝑥10−6 𝐹 2
=
𝜕𝑅𝑠 𝑇
𝑅𝑠
𝜕𝑅𝑠 𝑅𝑠
= 1,3984
𝜕𝑃 𝑇
𝑃
CORRELACIÓN PETROSKY Y FARSHAD
𝜕𝐵𝑜 𝛾𝑔 0,2914
= 83,313𝑥10−6 𝐹 2,0936 𝑅𝑠 −0,6262
𝜕𝑅𝑠 𝑇
𝛾𝑜 0,6265
𝜕𝑅𝑠 𝑅𝑠
= 1,73184
𝜕𝑃 𝑇
𝑃 + 139,05
VISCOSIDAD DEL PETRÓLEO
La estimación de la viscosidad del petróleo a presiones
iguales o por debajo del punto de burbuja se realiza en
dos pasos:
1. Se calcula la viscosidad del petróleo sin gas disuelto
(petróleo muerto), 𝜇𝑜𝑑 , a la temperatura del yacimiento.
2. Se ajusta la viscosidad del petróleo sin gas disuelto
(petróleo muerto), para tomar en cuenta el efecto de la
solubilidad del gas a la presión de interés.
A presiones mayores del punto de burbuja, es necesario
un tercer paso:
3. Ajustar la viscosidad al punto de burbuja para tomar en
cuenta la compresión y el grado de saturación del crudo
en el yacimiento
VISCOSIDAD DEL PETRÓLEO
MÉTODO PARA CALCULAR LA VISCOSIDAD DEL PETRÓLEO MUERTO
Beal, desarrollo una correlación gráfica para determinar la viscosidad del petróleo
muerto de acuerdo con la temperatura y la gravedad API del crudo. Posteriormente
Standig expresó esta correlación mediante la siguiente ecuación
𝑎
1,8 10 7 360
𝜇𝑜𝑑 = 0,32 +
𝐴𝑃𝐼 4,53 𝑇 + 200
𝑎 = 𝑎𝑛𝑡𝑖𝑙𝑜𝑔 0,43 + 8,33/𝐴𝑃𝐼
VISCOSIDAD DEL PETRÓLEO
MÉTODO PARA CALCULAR LA VISCOSIDAD DEL PETRÓLEO MUERTO
Correlación Beggs y Robinson, desarrollo una correlación basada en la correlación
gráfica de Beal para determinar la viscosidad del petróleo muerto.
−1,163
𝜇𝑜𝑑 = 10𝐴𝑇 − 1,0
𝐴 = 103,0324+0,02023𝐴𝑃𝐼
Donde 𝜇𝑜𝑑 es la viscosidad del petróleo muerto en cp medida a 14,7 lpca y
temperatura del yacimiento; y T, la temperatura en °𝐹
VISCOSIDAD DEL PETRÓLEO
MÉTODO PARA CALCULAR LA VISCOSIDAD DEL PETRÓLEO MUERTO
Correlación Glaso.
𝜇𝑜𝑑 = 3,141𝑥1010 𝑇 −3,444 (𝑙𝑜𝑔𝛾𝐴𝑃𝐼 )(10,313𝑙𝑜𝑔𝑇−36,447)
VISCOSIDAD DEL PETRÓLEO
MÉTODO PARA CALCULAR LA VISCOSIDAD DEL PETRÓLEOS SATURADOS
Correlación de Chew y Connally, desarrollo una correlación para ajusta 𝜇𝑜𝑑 de
acuerdo con la solubilidad del gas a la presión de saturación. Posteriormente Standig
expresó esta correlación mediante la siguiente ecuación
𝜇𝑜𝑏 = 10 𝑎 𝜇𝑜𝑑 𝑏
𝑎 = 𝑅𝑠 2,2 10−7 𝑅𝑠 − 7,4 10−4
0,68 0,25 0,062
𝑏= + 𝑑 +
10𝑐 10 10𝑒
𝑐 = 8,62 10−5 𝑅𝑠 𝑑 = 1,1 10−3 𝑅𝑠 e= 3,741 10−3 𝑅𝑠
Donde 𝜇𝑜𝑑 es la viscosidad del petróleo muerto en cp medida a 14,7 lpca y
temperatura del yacimiento; y 𝜇𝑜𝑏 , la viscosidad del petróleo en cp a la presión de
burbuja. y T, la temperatura en °𝑅
VISCOSIDAD DEL PETRÓLEO
Rango experimentales para usar la ecuación de Chew
VISCOSIDAD DEL PETRÓLEO
MÉTODO PARA CALCULAR LA VISCOSIDAD DEL PETRÓLEOS SATURADOS
Correlación de Beggs,
𝑏
𝜇𝑜𝑏 = 𝑎 𝜇𝑜𝑑
−0,515
𝑎 = 107,15 ∗ 𝑅𝑠 + 100
𝑏 = 5,44 𝑅𝑠 + 150 −0,338
VISCOSIDAD DEL PETRÓLEO
MÉTODO PARA CALCULAR LA VISCOSIDAD DEL PETRÓLEO NO SATURADOS
La viscosidad del petróleo por encima del punto de burbuja se estima, primero,
calculando la viscosidad al punto de burbuja y, luego, ajustándola para presiones
mayores.
Vásquez y Beggs propusieron la siguiente correlación para estimar dicha viscosidad:
𝑚
𝑃
𝜇𝑜 = 𝜇𝑜𝑏
𝑃𝑏
m= 2,6𝑃1,187 𝑒𝑥𝑝 −11,513 − 8,98𝑥10−5 𝑃
GASES REALES
• Los gases reales solo obedecen la ley de un gas ideal a bajas
presiones y temperaturas.
• En un gas real el volumen ocupado por las moléculas no es
despreciables.
• Existen fuerzas de atracción y repulsión entre las moléculas.
En el caso de la Industria petrolera, la ecuación general de los
gases puede ser utilizada dentro de ciertos límites para predecir el
comportamiento real de los gases, introduciendo un factor de
corrección o factor de desviación, así se tiene que:
Ecuación de compresibilidad
𝑃𝑉 = 𝑍𝑛𝑅𝑇 de estado o Ecuación para
gases reales.
La densidad de un gas real se expresa como:
𝑃𝑀
𝜌𝑔 =
PROFESORA: ING. YATNIELAH PIRELA
𝑧𝑅𝑇 YACIMIENTOS II
FACTOR DE COMPRESIBILIDAD (Z)
El factor de compresibilidad Z es determinado experimentalmente como la
razón del volumen real de n moles de un gas a determinada presión y
temperatura y su volumen ideal del mismo número de moles a la misma
presión y temperatura. Dicho valor no es una constante si no que varía con
los cambios en la composición de un gas, la temperatura y la presión.
𝑉𝑟𝑒𝑎𝑙
𝑧= Z gas ideal = 1
𝑉𝑖𝑑𝑒𝑎𝑙
El factor de compresibilidad, también es conocido como Factor de
supercompresibilidad, factor de desviación del gas o factor z
Gráfico típico del factor z, en
función de presión y a temperatura
constante
LEY DE ESTADO
CORRESPONDIENTE
La ley expresa que todos los gases puros
exhibirán un mismo comportamiento del factor
de compresibilidad z al mismo valor de presión
reducida y temperatura reducida.
𝑇 𝑃
𝑇𝑟 = 𝑃𝑟 =
𝑇𝑐 𝑃𝑐
DETERMINACIÓN DEL FACTOR DE
COMPRESIBILIDAD
MEZCLA DE GASES – PROPIEDADES SEUDOCRÍTICAS
REGLA DE KAY
Donde:
𝑦𝑗 = Fracción molar de j
𝑇𝑝𝑐 = 𝑦𝑗 𝑇𝑐𝑗 𝑃𝑝𝑐 = 𝑦𝑗 𝑃𝑐𝑗 𝑇𝑐𝑗 = Temperatura critica de
j
𝑗 𝑗 𝑃𝑐𝑗 = Presión critica de j
Estas propiedades seudocríticas fueron diseñadas
simplemente para su uso en la correlación de
propiedades físicas. Las propiedades seudocríticas no
son iguales a las propiedades críticas reales de una
mezcla de gases. Estas ecuaciones son conocidas
como Regla de mezcla de Kay.
DETERMINACIÓN DEL FACTOR DE
COMPRESIBILIDAD
MÉTODO STANDING Y KATZ
Basado en el concepto de propiedades seudoreducidas, Standing y Katz,
presentó un gráfico generalizado del factor de compresibilidad del gas. El
gráfico representa el factor de compresibilidad de gases dulces o con
pequeñas cantidades de componentes no hidrocarburos.
𝑇 𝑃
𝑇𝑝𝑟 = 𝑃𝑝𝑟 =
𝑇𝑝𝑐 𝑃𝑝𝑐
PASO 1: Determinar la presión y temperatura seudocritica.
Com Fracci Tempe Presió 𝒚𝒋 ∙ 𝑻𝒄𝒋 𝒚𝒋 ∙ 𝑷𝒄𝒋 𝑻𝒄𝒋 /𝑷𝒄𝒋 𝒚𝒋 . 𝑻𝒄𝒋 𝒚𝒋 .
𝒚𝒋 .
pone ón ratura n (𝑻𝒄𝒋 / 𝑷𝒄𝒋 𝑻𝒄𝒋
nte Molar Crítica Crítica 𝑷𝒄𝒋 ) 𝑷𝒄𝒋
𝒚𝒋 °R 𝒑𝒔𝒊𝒂 𝑷𝒄𝒋
𝑻𝒄𝒋
𝐶𝑂2 0,02 547,91 1071 10,96 21,42 0,512 0,010 0,014 0,335
𝑁2 0,01 227,49 493,1 2,27 4,93 0,461 0,005 0,007 0,102
𝐶1 0,85 343,33 666,4 291,83 566,44 0,515 0,438 0,610 11,305
𝐶2 0,04 549,92 706,5 22,00 28,26 0,778 0,031 0,035 0,827
𝐶3 0,03 666,06 616,0 19,98 18,48 1,081 0,032 0,031 0,805
i-𝐶4 0,03 734,46 527,9 22,03 15,84 1,391 0,042 0,035 0,959
n-𝐶4 0,02 765,62 550,6 15,31 11,01 1,390 0,028 0,024 0,653
1 384,38 666,38p 0,586 0,756 14,986
°𝑅=𝑻𝒑𝒄 sia=𝑷𝒑𝒄
Determinar la presión y temperatura seudocrítica por método de Kay
𝑇𝑝𝑐 = 𝑦𝑗 𝑇𝑐𝑗 = 384,38°𝑅 𝑃𝑝𝑐 = 𝑦𝑗 𝑃𝑐𝑗 = 666,38 𝑙𝑝𝑐𝑎
𝑗 𝑗
Determinar la presión y temperatura seudoreducida:
𝑇 639,69°𝑅 𝑃 3000𝑙𝑝𝑐𝑎
𝑇𝑝𝑟 = = = 1,66 𝑃𝑝𝑟 = = 4,50
𝑇𝑝𝑐 384,38°𝑅 𝑃𝑝𝑐 666,38 𝑙𝑝𝑐𝑎
𝑧 =0,85
Determinar la presión y temperatura seudocrítica por método de Stewart-Burkhadrt-Voo
2
𝑛 𝑛
𝑇𝑝𝑐 1 𝑇𝑐 2 𝑇𝑐 1 2 2
𝐽= = 𝑦𝑗 + 𝑦𝑗 = 0,586 + 0,756 = 0,576
𝑃𝑝𝑐 3 𝑃𝑐 𝑗
3 𝑃𝑐 𝑗
3 3
𝑗=1 𝑗=1
𝑛 𝐾2 14,986 2
𝑇𝑝𝑐 𝑇𝑐 𝑇𝑝𝑐 = = = 389,89°𝑅
𝐾= = 𝑦𝑗 = 14,986 𝐽 0,576
𝑃𝑝𝑐 𝑗=1
𝑃𝑐 𝑗
𝑇𝑝𝑐 389,89°𝑅
𝑃𝑝𝑐 = = = 676,90𝑙𝑝𝑐𝑎
𝐽 0.576
Determinar la presión y temperatura seudoreducida:
𝑇 639,69°𝑅 𝑃 3000𝑙𝑝𝑐𝑎
𝑇𝑝𝑟 = = = 1,64 𝑃𝑝𝑟 = = 4,43
𝑇𝑝𝑐 389,89°𝑅 𝑃𝑝𝑐 676,90 𝑙𝑝𝑐𝑎
𝑧 =0,84
PROPIEDADES SEUDOCRÍTICAS
DE UN GAS CUANDO NO SE
CONOCE SU COMPOSICIÓN
En caso de no conocer la composición del gas, las propiedades
seudocriticas pueden ser determinadas utilizando la gravedad especifica
del gas. Brown et al (1948) presento un método grafico para el calculo de
las propiedades seudocriticas. Standing (1977) expresó estas correlaciones
graficas en las siguientes ecuaciones matemáticas
Caso1: Sistemas de Gas Natural Caso 2: Sistemas de Gas Condensado
𝑇𝑝𝑐 = 168 + 325𝛾𝑔 − 12.5𝛾𝑔2 𝑇𝑝𝑐 = 187 + 330𝛾𝑔 − 71.5𝛾𝑔2
𝑃𝑝𝑐 = 677 + 15.0𝛾𝑔 − 37.5𝛾𝑔2 𝑃𝑝𝑐 = 706 − 51.7𝛾𝑔 − 11.1𝛾𝑔2
Donde:
𝑇𝑝𝑐 = Temperatura seudosrítica °𝑅
𝑃𝑝𝑐 = Presión seudocritica lpca
𝛾𝑔 = Gravedad especifica de la mezcla de gas
EJEMPLO
Calcular la densidad del ejercicio anterior, suponiendo que no se conoce la
composición de la mezcla del gas natural.
Calcular la gravedad especifica del gas
𝑙𝑏
𝑀𝑎 = 20,21
𝑙𝑏𝑚𝑜𝑙
𝑙𝑏
𝑀𝑎 20,21
𝛾𝑔 = = 𝑙𝑏𝑚𝑜𝑙 = 0,697
𝑀𝑎𝑖𝑟𝑒 𝑙𝑏
29
𝑙𝑏𝑚𝑜𝑙
Calcular la temperatura y presión seudocritica
𝑇𝑝𝑐 = 168 + 325𝛾𝑔 − 12.5𝛾𝑔2 = 168 + 325 0,697 − 12.5 0,697 2
= 388,45°𝑅
𝑃𝑝𝑐 = 677 + 15.0𝛾𝑔 − 37.5𝛾𝑔2 = 677 + 15.0 0,697 − 37.5 0,697 2
= 668,75lpca
Calcular la temperatura y presión seudoreducida
𝑇 639,69°𝑅 𝑃 3000 𝑙𝑝𝑐𝑎
𝑇𝑝𝑟 = = = 1,64 𝑃𝑝𝑟 = = = 4,48
𝑇𝑝𝑐 388,45°𝑅 𝑃𝑝𝑐 668,75 𝑙𝑝𝑐𝑎
Calcular el Factor de compresibilidad de la grafica de Standing y Katz
𝑧 =0,84
Calcular la densidad de la mezcla del gas
𝑚 𝑃𝑀𝑎
𝜌𝑔 = =
𝑉 𝑧𝑅𝑇
20,21𝑙𝑏
3000𝑙𝑝𝑐𝑎 ( )
𝜌𝑔 = 𝑙𝑏𝑚𝑜𝑙 = 10,51 lb/𝑝𝑖𝑒 3
𝑙𝑝𝑐𝑎 ∙ 𝑝𝑖𝑒 3
0,84 (10,73 )(639,69°𝑅)
°𝑅 ∙ 𝑙𝑏 − 𝑚𝑜𝑙
EFECTOS DE LOS COMPONENTES
NO HIDROCARBUROS EN EL
FACTOR Z
El gas natural usualmente contiene otros materiales adicionales a los
componentes hidrocarburos, como nitrógeno, dióxido de carbono y sulfuro
de hidrogeno. El gas natural se clasifica en gas dulce o gas ácido
dependiendo del contenido de sulfuro de hidrogeno. Ambos dulce o ácido
pueden contener nitrógeno, dióxido de carbono o ambos. Un gas
hidrocarburo se denomina gas acido si contiene un grano de sulfuro de
hidrogeno (𝐻2 𝑆) por 100𝑝𝑖𝑒 3
Concentraciones menores al 5 % de estos componentes no hidrocarburos
no afectarán gravemente la precisión. Sin embargo, errores en cálculos del
factor de compresibilidad mayores o iguales al 10 % pueden ocurrir a
mayores concentraciones de estos componentes no hidrocarburos en
mezclas de gases.
MÉTODO DE AJUSTE PARA
COMPONENTES NO HIDROCARBUROS
• Método de corrección de Wichert-Aziz
• Método de corrección Carr-Kobayashi-Burrows
Método de corrección Wichert-Aziz
Gases naturales que contienen 𝐻2 𝑆 y/o 𝐶𝑂2 frecuentemente muestran un
comportamiento del factor de compresibilidad diferente a los gases dulces.
Wichert y Aziz desarrollaron un procedimiento simple para calcular estas
diferencias. El método permite usar la grafica de Standing y Katz, usando
un factor de ajuste de temperatura seudocritica 𝜀 , el cual depende de la
concentración de 𝐻2 𝑆 y 𝐶𝑂2 en un gas ácido. Este factor de corrección es
luego usado para ajustar la presión y la temperatura seudocritica con las
siguientes expresiones:
′ =𝑇 −𝜀 Donde:
𝑇𝑝𝑐 𝑝𝑐 𝑇𝑝𝑐 =Temperatura seudocritica °𝑅
𝑃𝑝𝑐 =Presión seudocritica lpca
′
𝑃 𝑇 ′ 𝑇𝑝𝑐 =Temperatura seudocritica corregida °𝑅
′ = 𝑝𝑐 𝑝𝑐 ′
𝑃𝑝𝑐 𝑃𝑝𝑐 =Presión seudocritica corregida lpca
𝑇𝑝𝑐 + 𝐵(1 − 𝐵)𝜀 B= Fracción molar de 𝐻2 𝑆 en la mezcla de gas
𝜀= Factor de ajuste de temperatura seudocrítica
Factor de ajuste de temperatura seudocritica se expresa como:
𝜀 = 120 𝐴0.9 − 𝐴1.6 +15 𝐵0.5 − 𝐵4.0
Donde A es la suma de la fracción molar de 𝐻2 𝑆 y 𝐶𝑂2 en la mezcla de gas
𝐴 = 𝑦𝐻2𝑆 + 𝑦𝐶𝑂2
B= Fracción molar de 𝐻2 𝑆 en la mezcla de gas
Método de corrección Carr-Kobayashi-Burrows
Carr, Kobayashi y Burrows propusieron un método simple para ajustar las
propiedades seudocríticas del gas natural cuando componentes no
hidrocarburos están presentes. El método puede ser usado cuando la
composición del gas se desconoce.
′ = 𝑇 − 80𝑦
𝑇𝑝𝑐 𝑝𝑐 𝐶𝑂2 + 130𝑦𝐻2 𝑆 − 250𝑦𝑁2
′ = 𝑃 + 440𝑦
𝑃𝑝𝑐 𝑝𝑐 𝐶𝑂2 + 600𝑦𝐻2 𝑆 − 170𝑦𝑁2
Donde:
𝑇𝑝𝑐 =Temperatura seudocritica °𝑅
𝑃𝑝𝑐 =Presión seudocritica lpca
′
𝑇𝑝𝑐 =Temperatura seudocritica corregida °𝑅
′
𝑃𝑝𝑐 =Presión seudocritica corregida lpca
𝑦𝐶𝑂2 = Fracción molar de dióxido de carbono
𝑦𝐻2 𝑆 = Fracción molar de sulfuro de hidrogeno
𝑦𝑁2 =Fracción molar de Nitrógeno
Universidad de América
Propiedades de Fluidos del Yacimiento
LIMITACIONES AL USAR EL MÉTODO GRAFICO DE
STANDING KATZ
1. El gas debe ser rico en Metano (C1>90%)
2. El gas no debe tener impureza. El contenido de
20% de N2 produce un error del 4%. La presencia
de CO2 produce un error en el cálculo del factor Z
igual al valor del porcentaje de Co2 en la mezcla.
3. No presenta resultados a presiones y temperaturas
cercanas a la crítica.
4. No se recomienda su uso en el calculo de z a
presiones mayores de 10.000 lpca.
CORRELACIONES PARA EL CÁLCULO
DEL FACTOR Z
• Correlación de • Método de Papay, J
dranchuk. • Método de Sarem, A.
• Correlación de Brill, M.
j.p. y Beggs, H.D. • Método grafico de
• Correlación de Standing y Katz.
Yarborough y Hall.
Universidad de América
Propiedades de Fluidos del Yacimiento
MÉTODO DE PAPAY, J
La ecuación para el cálculo de z es la
siguiente
3,52𝑃𝑠𝑟 2
0,274𝑃𝑠𝑟
z = 1 − 0,9813𝑇 + 0,8157𝑇
10 𝑠𝑟 10 𝑠𝑟
Takacs, determinó que el error promedio de este método con
respecto a los valores leído en las curvas de Standing y Katz fue de
-4,873% para presiones y temperaturas seudorreducidas en el rango
de: 0,2 ≤ 𝑃𝑠𝑟 ≤ 15,0 y 1,2≤ 𝑇𝑠𝑟 ≤3,0. El método de Papay puede ser
utilizado cuando el gas natural contiene impurezas corrigiendo
𝑃𝑠𝑐 𝑦 𝑇𝑠𝑐 por el método de Wichert y Aziz.
Universidad de América
Propiedades de Fluidos del Yacimiento
MÉTODO DE SAREM, A. M.
Este método para determinar Z se basa en los Polinomios de
Legendre de grado 0 a 5. La ecuación básica de ajuste es:
5 5
𝑧= 𝐴𝑖𝑗 𝑃𝑖 𝑥 𝑃𝑗 𝑦 Ec. 1
𝑖=0 𝑗=0
donde: 2𝑃𝑠𝑟 − 15
𝑥=
14,8
2𝑇𝑠𝑟 − 4
𝑦=
1,9
Los Polinomios de Legendre de grado 0-5, 𝑃𝑖 y 𝑃𝑗 tienen las
siguientes formas:
𝑃0 𝑎 = 0,7071068
𝑃1 𝑎 = 1,224745𝑎
𝑃2 𝑎 = 0,7905695 3𝑎2 − 1
𝑃3 𝑎 = 0,9354145 5𝑎3 − 3𝑎
𝑃4 𝑎 = 0,265165 35𝑎4 − 30𝑎2 + 3
𝑃5 𝑎 = 0,293151 63𝑎5 − 70𝑎3 + 15𝑎
Donde “a” se reemplaza por “x” y “y” al efectuar la sumatoria de la
ecuación 1
Universidad de América
Propiedades de Fluidos del Yacimiento
La tabla presenta los valores del coeficiente 𝐴𝑖𝑗
Para los intervalos 0,1≤𝑃𝑠𝑟 ≤14,9 y 1,05 ≤ 𝑇𝑠𝑟 ≤2,95, el error del
método con respecto a los valores leídos en las curvas de
Standing y Katz fue menor de 0,4%. El método de Sarem puede
ser utilizado cuando el gas natural contiene impurezas
corrigiendo 𝑃𝑠𝑐 y 𝑇𝑠𝑐 por el método de Wichert y Aziz.
Universidad de América
Propiedades de Fluidos del Yacimiento
EJEMPLO DE
PROCEDIMIENTO
MÉTODO SAREM
X -0,255790394 Y -0,68723024
P0(X) 0,7071068 P0(Y) 0,7071068
P1(X) -0,313278006 P1(Y) -0,8416818
P2(X) -0,635391636 P2(Y) 0,3295538
P3(X) 0,63953446 P3(Y) 0,41050355
P4(X) 0,314742998 P4(Y) -0,89140241
P5(X) -0,801568425 P5(Y) 0,80738156
i j=0 j=1 j=2 j=3 j=4 j=5
0 2,1433504 0,0831762 -0,021467 -0,0008714 0,0042846 -0,0016595
1 0,3312352 -0,1340361 0,066881 -0,0271743 0,0088512 -0,002152
2 0,1057287 -0,0503937 0,0050925 0,0105513 -0,0073182 0,002696
3 0,052184 0,0443121 -0,0193294 0,0058973 0,0015367 -0,0028327
4 0,019704 -0,0263834 0,019262 -0,0115354 0,004291 -0,0081303
5 0,0053096 0,0089178 -0,0108948 0,0095594 -0,0060114 0,0031175
5 5
𝑧= 𝐴𝑖𝑗 𝑃𝑖 𝑥 𝑃𝑗 𝑦
𝑖=0 𝑗=0
MÉTODO DE HALL, K.R. Y YARBOROUGH,
L.
Hall y Yarborough basados en la ecuación de estado de
Starling-Carmahan desarrollaron las siguientes ecuaciones para
el cálculo de z: 2
0,06125𝑃𝑠𝑟 𝑡 𝑒𝑥𝑝 −1,2 1 − 𝑡 Ec. 2
𝑧=
𝑦
Donde:
𝑡 =Recíproco de la temperatura seudorreducida 𝑡 = 𝑇𝑠𝑐 /𝑇
𝑦 =Densidad reducida, la cual se obtiene a partir de la solución
de la siguiente ecuación:
𝑦 + 𝑦2 + 𝑦3 − 𝑦4 Ec. 3
−𝐴𝑝𝑠𝑟 + 3
− 𝐵𝑦 2 + 𝐶𝑦 𝐷 = 0
1−𝑦
2
𝐴 = 0,06125 𝑡 𝑒𝑥𝑝 −1,2 1 − 𝑡
𝐵 = 14,76𝑡 − 9,76𝑡 2 + 4,58𝑡 3
𝐶 = 90,7𝑡 − 242,2𝑡 2 + 42,2𝑡 3
𝐷 = 2,18 + 2,82𝑡
Takacs, determinó que el error promedio en el cálculo de z por este
método fue de -0,518% en comparación con el obtenido por el
método de Standing y Katz. Para obtener buenos resultados con
este método se recomienda que 𝑃𝑠𝑟 y 𝑇𝑠𝑟 estén entre los siguientes
rangos: 0,1 ≤ 𝑃𝑠𝑟 ≤ 24,0 y 1,2≤ 𝑇𝑠𝑟 ≤3,0. Para gases naturales con
impurezas se recomienda corregir 𝑃𝑠𝑐 𝑦 𝑇𝑠𝑐 por el método de Wichert
y Aziz.
Debido a que la ecuación anterior es no lineal, se requiere una
solución de ensayo y error para resolverla. Un método de
ensayo y error frecuentemente utilizado, es el método de
Newton-Raphson el cual utiliza el siguiente procedimiento
iterativo (aplicado a este caso):
1. Suponer un valor inicial de 𝑦1 y calcular 𝐹 𝑦1 , donde 𝐹 𝑦
es el término de la derecha de la ecuación 3.
2. Si 𝐹 𝑦1 ≅ 0 ó se encuentra dentro de una tolerancia
específica ±10−4 , se puede considerar que 𝑦1 es la
solución. En caso de lo contrario calcular un nuevo valor de
y utilizando la siguiente aproximación (Serie Taylor):
𝐹 𝑦1
𝑦2 = 𝑦1 − 𝑑𝐹 𝑦1 ,
𝑑𝑦
𝑑𝐹 𝑦1
donde se obtiene derivando la ecuación 3 con respecto y
𝑑𝑦
a 𝑇𝑠𝑟 constante, esto es:
𝑑𝐹 1 + 4𝑦 + 4𝑦 2 − 4𝑦 3 + 𝑦 4 𝐷−1
= − 2𝐵𝑦 + 𝐶𝐷𝑦
𝑑𝑦 1−𝑦 4
Universidad de América
Propiedades de Fluidos del Yacimiento
3. Hacer 𝑦1 = 𝑦2 y repetir el paso 2. Continuar el
procedimiento hasta obtener la solución.
4. Sustituir el valor correcto de y en la ecuación 2 para obtener
z
Universidad de América
Propiedades de Fluidos del Yacimiento
MÉTODO DE BRILL, J.P. Y BEGGS, H.D.
La ecuación para el cálculo de z es la siguiente:
1−𝐴
𝑧 =𝐴+ + 𝐶𝑃𝑠𝑟 𝐷
𝑒𝑥𝑝 𝐵
Donde:
0,5
𝐴 = 1,39 𝑇𝑠𝑟 − 0,92 − 0,36𝑇𝑠𝑟 − 0,10
0,066 2 +
0,32 6
𝐵 = 0,62 − 0,23𝑇𝑠𝑟 𝑃𝑠𝑟 + − 0,037 𝑃𝑠𝑟 𝑃𝑠𝑟
𝑇𝑠𝑟 − 0,86 109 𝑇𝑠𝑟−1
𝐶 = 0,132 − 0,32𝑙𝑜𝑔𝑇𝑠𝑟
2
𝐷 = 𝑎𝑛𝑡𝑖𝑙𝑜𝑔 0,3106 − 0,49𝑇𝑠𝑟 + 0,1824𝑇𝑠𝑟
Universidad de América
Propiedades de Fluidos del Yacimiento
La exactitud entre los valores calculados y leídos del grafico de
Standing y Katz fue verificado para valores de 0, 1, 2, 3, 4, 5, 7,
9, 11 y 13, y valores de 𝑇𝑠𝑟 de 1.2, 11.3, 1.4, 1.5, 1.7, 2.0 y 2.4,
los resultados fueron los siguientes:
El método no es válido para temperaturas seudoreducidas fuera
del rango de 1.2 a 2.4.
Para gases naturales con impurezas se recomienda corregir
𝑃𝑠𝑐 𝑦 𝑇𝑠𝑐 por el método de Wichert y Aziz.
Universidad de América
Propiedades de Fluidos del Yacimiento
MÉTODO DE DRANCHUK, P.M., PURVIS, R.A. Y
ROBINSON, D.B.
Este método es el resultado de un ajuste realizado a la
ecuación de estado de Benedict, Webb y Rubin (BWR), la cual
escribieron en la siguiente forma:
La densidad reducida 𝜌𝑟 , se calcula mediante la siguiente
expresión (obtenida de la ley de los gases):
𝜌 𝑧𝑐 𝑃𝑠𝑟 𝑃𝑠𝑟
𝜌𝑟 = = = 0,27 Ec. 4
𝜌𝑐 𝑧𝑇𝑠𝑟 𝑧𝑇𝑠𝑟
En esta ecuación se tomó el factor de compresibilidad del gas
en el punto crítico, 𝑧𝑐 = 0,27, considerado como un valor
apropiado para mezclas compuestas principalmente por
metano.
Universidad de América
Propiedades de Fluidos del Yacimiento
Además, utilizando datos de 1500 puntos, se determinaron los
siguientes valores para las constantes de 𝐴1 − 𝐴8 :
El método reprodujo los 1500 puntos de datos con error
absoluto de 0,54% y una desviación estándar de 0,00445. Los
rangos de aplicación son 0,2 ≤ 𝑃𝑠𝑟 ≤ 30,0 y 1,05≤ 𝑇𝑠𝑟 ≤3,0. Para
gases naturales con impurezas se recomienda corregir 𝑃𝑠𝑐 𝑦 𝑇𝑠𝑐
por el método de Wichert y Aziz.
Universidad de América
Propiedades de Fluidos del Yacimiento
Para resolver la formulación implícita de la densidad reducida,
se debe utilizar un procedimiento iterativo de ensayo y error. Se
utiliza el método iterativo de Newton-Raphson, se tiene:
𝐹 𝑍
𝑍2 = 𝑍1 − 𝑑𝐹 𝑍11 ,
𝑑𝑧
MÉTODO DE DRANCHUK, P.M., ABOU-KASSEM, J.H.
Este método es muy similar al anterior, sin embargo; Dranchuk
y Abou-Kassem utilizaron la ecuación de estado de Starling la
cual escribieron en la siguiente forma:
Universidad de América
Propiedades de Fluidos del Yacimiento
De nuevo utilizó la ecuación 4 para la densidad reducida y
mediante un procedimiento similar al método anterior y
utilizando los mismos 1500 puntos, determinaron los siguientes
valores para las constantes 𝐴1 − 𝐴11 :
El método reprodujo los 1500 puntos de datos con error
absoluto de 0,307% y una desviación estándar de 0,00378. Los
rangos de aplicación son: 0,2 ≤ 𝑃𝑠𝑟 ≤ 30,0 y 1,0 ≤ 𝑇𝑠𝑟 ≤3,0, y para
𝑃𝑠𝑟 < 1.0 con 0,7 < 𝑇𝑠𝑟 ≤ 1.0 , pero produce resultados
inaceptables en la región de 𝑇𝑠𝑟 = 1.0 y 𝑃𝑠𝑟 ≥ 1.0. Para gases
naturales con impurezas se recomienda corregir 𝑃𝑠𝑐 𝑦 𝑇𝑠𝑐 por el
método de Wichert y Aziz.
Universidad de América
Propiedades de Fluidos del Yacimiento
Para resolver la formulación implícita de la densidad reducida,
se debe utilizar un procedimiento iterativo de ensayo y error. Se
utiliza el método iterativo de Newton-Raphson, se tiene:
Universidad de América
Propiedades de Fluidos del Yacimiento
MÉTODO DE GOPAL, V.N.
Gopal desarrolló las siguientes ecuaciones para el cálculo del
factor z
FACTOR VOLUMÉTRICO DEL GAS
Es el volumen de gas a condiciones de yacimientos requerido para
producir una unidad volumétrica de gas a condiciones normales de
superficie.
Puede calcularse como el cociente entre el volumen ocupado por el
gas a temperatura y presiones de yacimiento, y el volumen
ocupado por el mismo peso de gas a condiciones normales.
𝑉𝑦
𝐵𝑔 =
𝑉𝑐𝑛
Factor volumétrico de gas
Típica gráfica del factor
de formación, 𝐵𝑔
volumétrico de gas de formación,
como función de la presión a
temperatura de yacimiento
constante
Presión de yacimiento
El volumen de n moles de un gas a condiciones de yacimiento
puede ser obtenido a partir de la ecuación de estado:
𝑧𝑛𝑅𝑇
𝑉𝑦 =
𝑃
Donde P y T representan la temperatura y presión de yacimiento.
Y el volumen del mismo número de moles de un gas a
condiciones normales, 𝑃𝑐𝑛 y 𝑇𝑐𝑛 , es:
𝑧𝑐𝑛 𝑛𝑅𝑇𝑐𝑛
𝑉𝑐𝑛 =
𝑃𝑐𝑛
Entonces, el factor volumétrico del gas de formación se puede
expresar de la siguiente manera:
𝑧𝑛𝑅𝑇
𝑉𝑦 𝑃 𝑧𝑇𝑃𝑐𝑛
𝐵𝑔 = = =
𝑉𝑐𝑛 𝑧𝑐𝑛 𝑛𝑅𝑇𝑐𝑛 𝑧𝑐𝑛 𝑇𝑐𝑛 𝑃
𝑃𝑐𝑛
Tomando en cuenta que las condiciones normales para 𝑇𝑐𝑛 = 520 R,
𝑃𝑐𝑛 = 14.65 𝑝𝑠𝑖𝑎 y 𝑍𝑐𝑛 = 1:
𝑧𝑇 14.65 𝑧𝑇 𝑃𝐶𝑌
𝐵𝑔 = = 0.0282 𝑃𝐶𝑁
1.0 520 𝑃 𝑃
También
𝑧𝑇 𝑃𝐶𝑌 𝐵𝑌 𝑧𝑇 𝐵𝑌
𝐵𝑔 = 0.0282 = 0,00502
𝑃 𝑃𝐶𝑁 5,615𝑃𝐶𝑌 𝑃 𝑃𝐶𝑁
T=°𝑅
P=lpca
COMPRESIBILIDAD DE LOS GASES
Se define como el cambio relativo en volumen por unidad de
cambio en presión, lo cual matemáticamente se expresa como:
1 𝜕𝑉
𝐶𝑔= −
𝑉 𝜕𝑃 𝑇
1 𝜕𝑉𝑚 𝐶𝑔 𝑙𝑝𝑐 −1
𝐶𝑔= −
𝑉𝑚 𝜕𝑃 𝑇
Se obtiene una ecuación generalizada:
1 1 𝜕𝑧
𝐶𝑔 = −
𝑃 𝑧 𝜕𝑃 𝑇
𝜕𝑧
Para un gas ideal, Z=1 y = 0 , la ecuación de la
𝜕𝑃 𝑇
compresibilidad del gas se expresaría como:
1
𝐶𝑔 =
𝑃
COMPRESIBILIDAD SEUDOREDUCIDA
La ley de estados correspondientes puede ser usada para reducir la
ecuación de compresibilidad del gas.
Sustituyendo la ecuación 5 en la ecuación 4:
𝑃 = 𝑃𝑠𝑐 𝑃𝑠𝑟 Ec. ….5
1 1 𝜕𝑧
𝐶𝑔 = −
𝑃𝑠𝑐 𝑃𝑠𝑟 𝑧 𝜕𝑃𝑠𝑐 𝑃𝑠𝑟 𝑇𝑠𝑟
Multiplicando la ecuación anterior por 𝑃𝑠𝑐 se obtiene:
1 1 𝜕𝑧
𝐶𝑔 𝑃𝑠𝑐 = − = 𝐶𝑠𝑟
𝑃𝑠𝑟 𝑧 𝜕𝑃𝑠𝑟 𝑇𝑠𝑟
El término 𝐶𝑠𝑟 es llamado Compresibilidad seudoreducida y puede ser
expresada como: Donde: 𝐶𝑔 = Compresibilidad del gas 𝑝𝑠𝑖 −1
𝐶𝑠𝑟 = Compresibilidad seudoreducida
𝐶𝑠𝑟 = 𝐶𝑔 𝑃𝑠𝑐 𝑃𝑠𝑐 = 𝑃𝑟𝑒𝑠𝑖ó𝑛 𝑆𝑒𝑢𝑑𝑜𝑐𝑟í𝑡𝑖𝑐𝑎 psi
COMPRESIBILIDAD SEUDOREDUCIDA
𝜕𝑧
Los Valores de Pueden calcularse de la pendiente de la
𝜕𝑃𝑠𝑟 𝑇𝑠𝑟 curva isotérmica de Tsr en la
Grafica generalizada Z vs Psr de Standing y Katz.
También puede calcularse a través de
correlaciones.
GOPAL, V.N.
VISCOSIDAD DE LOS GASES
Si no se conoce la composición del gas,
el siguiente grafico puede utilizarse para
obtener la viscosidad de la mezcla de
gases a presión de una atmósfera.
También se muestra la corrección del
valor de la viscosidad, que debe ser
aplicado para tomar en cuenta el efecto
de la presencia de componentes no
hidrocarburos.
Este método gráfico ha sido comparado
con los datos experimentales y es
bastante acertado para hidrocarburos
gaseosos, con errores menores a 2 %.
Cuando existe presencia de
componentes no hidrocarburos, causa
una desviación entre los resultados
gráficos y los experimentales mayores al
2%, aún después de realizar la
corrección.
VISCOSIDAD DEL GAS A ALTAS PRESIONES
Carr, Kobayashi y Burrows, desarrollaron una
correlación gráfica basada en el teorema de los
estados correspondientes para estimar la viscosidad
del gas natural en función de temperatura, presión y
gravedad específica.
VISCOSIDAD DEL GAS A ALTAS PRESIONES
El procedimiento es el siguiente:
1. Calcular la presión seudocritica, temperatura seudocritica y el peso
molecular aparente a partir de la gravedad especifica o de la
composición del gas. Corregir las propiedades críticas por las presencia
de componentes no hidrocarburos (𝐶𝑂2, 𝐻2 𝑆, 𝑁2 ), para concentraciones
mayores 5 % mol.
2. Obtener la viscosidad del gas natural a una atmósfera y a la temperatura
de interés (𝜇𝑔1 ). Esta viscosidad debe ser corregida por la presencia de
componentes no hidrocarburos. El efecto de los componentes no
hidrocarburos puede ser expresado matemáticamente por la siguiente
ecuación
𝜇𝑔1 = 𝜇𝑔1(𝑛𝑜 𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑎) + ∆𝜇 𝑁2 + ∆𝜇 𝐶𝑂2 + ∆𝜇 𝐻2 𝑆
Donde
𝜇𝑔1 = viscosidad del gas corregida a una atmósfera de presión y a la T del yacimientos, cp
𝜇𝑔1(𝑛𝑜 𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑎) = viscosidad del gas no corregida , cp
∆𝜇 𝑁2 = Corrección de la viscosidad por la presencia de 𝑁2
∆𝜇 𝐶𝑂2 = Corrección de la viscosidad por la presencia de 𝐶𝑂2,
∆𝜇 𝐻2 𝑆 = Corrección de la viscosidad por la presencia de 𝐻2 𝑆
Fundación Universidad de América
Propiedades de Fluidos del Yacimiento
3. Calcular la presión y temperatura seudoreducida.
4. Con la presión y temperatura seudoreducida se lee en el gráfico
la relación de viscosidad 𝜇𝑔 /𝜇𝑔1 . El término 𝜇𝑔 represente la
viscosidad del gas a la temperatura y presión requerida.
5. La viscosidad del gas 𝜇𝑔 , a temperatura y presión de interés es
calculada multiplicando la viscosidad a una atmósfera de presión,
𝜇𝑔1 , por la relación de viscosidad, 𝜇𝑔 /𝜇𝑔1 .
Relación de viscosidad 𝜇𝑔 /𝜇𝑔1 para gases
naturales
OTRO MÉTODO PARA EL CÁLCULO DE
VISCOSIDAD DEL GAS
Lee, González y Eakin, desarrollaron una correlación para estimar la
viscosidad del gas natural en función de temperatura del yacimiento,
densidad del gas y el peso molecular del gas. La ecuación esta
expresada como:
𝜌 𝑌 𝜌𝑔 𝑌
𝑔 𝑋
𝜇𝑔 = 10−4 𝐾𝑒𝑥𝑝 𝑋 𝜇𝑔 = 10−4 𝐾𝑒 62,4
62,4
Donde:
9,4 + 0,02𝑀𝑎 𝑇 1.5 𝜌𝑔 = Densidad del gas a presión y temperatura
𝐾= del yacimientos 𝑙𝑏/𝑝𝑖𝑒 3
209 + 19𝑀𝑎 + 𝑇 T= Temperatura del Yacimientos °𝑅
𝑀𝑎 = Peso molecular aparente de la mezcla de
986 gas
𝑋 = 3,5 + + 0,01𝑀𝑎
𝑇
Esta correlación puede predecir valores de viscosidad
con una desviación estándar de 2,7% y un máximo de
desviación de 8,99%. Esta correlación es menos
𝑌 = 2,4 − 0,2𝑋 precisa para gases con alta gravedad específica. Este
método no debe ser usado para gases ácidos.
BIBLIOGRAFÍA
• PARIS DE FERRER, Magdalena. Correlaciones para sistemas de petróleo. Fundamentos de
ingeniería de yacimientos. Maracaibo, Venezuela. Ediciones astro Data S.A, 2009. P. 123-127.
• ESCOBAR MACUALO, Freddy Humberto. Correlaciones para un sistema de petróleos. Fundamentos
de ingeniería de yacimientos. Neiva, Huila, Colombia. Editorial Universidad Surcolombiana. P. 145-
150.
• BÁNZER S, Carlos. Propiedades físicas del petróleo. Correlaciones numéricas PVT. Maracaibo,
Venezuela. Universidad del Zulia. Instituto de investigaciones petroleras, 1996. P. 48-59.
• AHMED, TAREK. PROPIEDADES DE LOS FLUIDOS. YACIMIENTOS II
• ARAMENDIZ J.J Y VELAZQUEZ M.A. CONSIDERACIONES Y PROCEDIMIENTOS PARA EL
ANÁLISIS PVT DEL CRUDO DE CAMPOS MADUROS. APLICACIÓN CAMPO COLORADO.
• CRAFT B.C. Y HAWKINS M.F. INGENIERÍA APLICADA DE YACIMIENTOS PETROLÍFEROS.
PRIMERA EDICIÓN. EDITORIAL TECNOS. MADRID, 1968