Cuaderno de Ejercicios de Propiedades de Los Fluidos Petroleros
Cuaderno de Ejercicios de Propiedades de Los Fluidos Petroleros
FACULTAD DE INGENIERÍA
Cuaderno de Ejercicios de
Propiedades de los Fluidos
Petroleros
MATERIAL DIDÁCTICO
Que para obtener el título de
Ingeniero Petrolero
PRESENTA
M. en I.
Eduardo Maldonado Cruz
A Dios, porque es por su amor y ayuda que hoy se ha podido concluir este gran logro. Gracias
Señor, concédeme ser una profesionista que ejerza con excelencia y ética y que este y mis demás
logros lleguen a agradar tu corazón.
A mi querido esposo Víctor Hugo Vásquez, gracias por haber sido un gran amigo y apoyo
para terminar este trabajo, Te Amo.
A mi pequeño hijo Hugo Abraham Vásquez, te amo mi niño, tu sola presencia ha motivado
cada esfuerzo.
A mis padres Justino López y Reina González, gracias por su apoyo durante y después de
haber concluido la carrera, son los mejores padres del mundo, los amo.
A mis hermanos Virginia López, Kevin López, Iván López, por su gran apoyo durante la
carrera y en esta recta final, los quiero mucho.
A mi suegra Alejandra Hernández, por haberme y habernos apoyado como familia, sobre todo
en estos últimos días, es una gran persona, mujer y ejemplo a seguir.
A mi director el Ingeniero Eduardo Maldonado Cruz, por todo su apoyo, ayuda, paciencia y
conocimientos que compartió conmigo, ha sido un honor para mí realizar este trabajo bajo su
dirección, muchas gracias por todo.
A mis sinodales Ing. Hector Erick Gallardo, Ing. Patricia Vega, Ing. Agustín Velazco e Ing.
Israel Castro, gracias por su apoyo, ayuda y tiempo para la conclusión de este trabajo.
3
ÍNDICE
ÍNDICE ......................................................................................................................
4
CAPÍTULO 1. .............................................................................................................
17
INTRODUCCIÓN .......................................................................................................
17
CAPÍTULO 2. .............................................................................................................
20
CAPÍTULO 3. .............................................................................................................
37
4
COMPORTAMIENTO DE FASES. ............................................................................
37
3.1 Definiciones y tipos de diagramas de fase. .......................................................
37
3.2 Sistemas de un solo componente. ....................................................................
38
3.2.1 Punto de burbuja. .......................................................................................
40
3.2.2 Punto de rocío. ...........................................................................................
40
3.2.3 Envolvente de saturación. ..........................................................................
40
3.2.4 Diagrama de fase de densidad-temperatura para un componente puro.....
41
3.2.5 Ley de los diámetros rectilíneos. ................................................................
42
3.2.6 Ecuación de Rackett (1911) .......................................................................
42
3.2.7 Ecuación de Rackett modificada. ...............................................................
43
3.2.8 Métodos para determinar la presión de vapor para una sustancia pura. ....
43
3.2.8.1 Método de Clausius-Clapeyron. ..............................................................
43
3.3 Sistemas binarios. .............................................................................................
46
3.4 Sistemas ternarios. ...........................................................................................
47
3.5 Sistemas multicomponentes de hidrocarburos. ................................................
48
3.6 EJERCICIOS RESUELTOS. .............................................................................
50
5
CAPÍTULO 4. .............................................................................................................
71
CAPÍTULO 5. .............................................................................................................
95
7
5.1.4 Calor latente de vaporización ∆Hm. ............................................................
99
5.2 Comportamiento de líquidos ideales. ..............................................................
100
5.3 Comportamiento de líquidos reales. ...............................................................
101
5.3.1 Coeficiente de viscosidad. ........................................................................
102
5.3.2 Dinámica de fluidos reales. .......................................................................
102
5.4 EJERCICIOS RESUELTOS ............................................................................
105
CAPÍTULO 6. ...........................................................................................................
111
9
7.2.2 Coeficiente de compresibilidad isotérmico del gas Cg. .............................
147
7.2.3 Coeficiente de viscosidad del gas(μg). .....................................................
149
7.3 Propiedades del aceite negro. ........................................................................
154
7.3.1 Densidad relativa del aceiteγo. .................................................................
154
7.3.2 Densidad en °API. ....................................................................................
154
7.3.3 Factor de volumen de formación de aceite BO. ........................................
155
7.3.4 Relación gas en solución aceite (Rs). .......................................................
155
7.3.5 Factor de volumen total de la formación Bt. ..............................................
156
7.3.6 Coeficiente de compresibilidad isotérmica del aceite cO. .........................
157
7.3.7 Coeficiente de viscosidad del aceite μo. ...................................................
159
7.4 Propiedades del agua de formación. ........................................................... 160
7.4.1 Composición. ............................................................................................
160
7.4.2 Factor de volumen del agua de
formación. ...............................................160
7.4.3 Densidad del agua de formación. .............................................................
161
7.4.4 Relación de solubilidad de gas en el agua. ..............................................
161
7.4.5 Coeficiente de compresibilidad isotérmica del agua. ................................
162
7.4.6 Coeficiente de viscosidad del agua de formación. ....................................
163
7.5 EJERCICIOS RESUELTOS ............................................................................
164
10
CAPÍTULO 8 ............................................................................................................ 175
11
CORRELACIONES PARA ESTIMAR LAS PROPIEDADES PVT DE LOS
HIDROCARBUROS. ....................................................................................................
191
9.1 Importancia de los modelos empíricos. ...........................................................
191
9.2 Correlaciones para estimar las propiedades del gas. .....................................
192
9.2.1 Presión y temperatura pseudocríticas. .....................................................
192
9.2.2 Factor de desviación Z. .............................................................................
193
9.2.3 Factor de volumen de gas. .......................................................................194
9.2.4 Compresibilidad del gas. ..........................................................................
194
9.2.5 Viscosidad del gas. ................................................................................... 195
9.3 Correlaciones para estimar las propiedades del aceite. ................................. 195
9.3.1 Corrección de la densidad relativa del gas por las condiciones del
separador.
............................................................................................................................... 195
9.3.2 Presión de burbuja: ...................................................................................
196
9.3.3 Relación gas-aceite. .................................................................................
196
9.3.4 Factor volumen del aceite. ........................................................................
197
9.3.5 Compresibilidad isotérmica del aceite. .....................................................
197
9.3.6 Viscosidad del aceite. ...............................................................................
198
9.4 Correlaciones para estimar las propiedades del agua. ...................................
199
9.4.1 Factor de volumen del agua. ....................................................................
199
12
9.4.2 Compresibilidad del agua. ........................................................................
200
9.4.3 Viscosidad del agua. .................................................................................
200
9.5 Proceso de ajuste de las correlaciones con datos de campo. ........................
201
9.5.1 Inclusión del gas seco. .............................................................................
201
9.6 EJERCICIOS RESUELTOS ............................................................................
202
13
10.3 Composiciones características de los fluidos
petroleros. ..............................210
10.3.1 Aceites negros. .......................................................................................
211
10.3.2 Aceites volátiles. .....................................................................................
212
10.3.3 Condensados de gas. .............................................................................
213
10.3.4 Gases húmedos. .................................................................................... 214
10.3.5 Gas seco ................................................................................................ 214
10.3.6 Criterios de clasificación. ........................................................................
215
10.4 Diagramas de fase característicos de los fluidos petroleros. ........................
216
10.4.1 Diagramas de fase de presión–temperatura para aceites negros. .........
216
10.4.2 Diagramas de fase de presión-temperatura para aceites volátiles. ........
217
10.4.3 Diagrama de fase de presión-temperatura para gas y condensados
(retrógrado). ........................................................................................................... 218
10.4.4 Diagrama de fase de presión-temperatura para gases húmedos. ..........
219
10.4.5 Diagrama de fase presión-temperatura para gas seco. ..........................
220
BIBLIOGRAFÍA .......................................................................................................
235
APENDICE A. ....................................................................................................... 237
Solución de una ecuación cúbica de estado. ........................................................ 237
Apéndice B ........................................................................................................... 239
Propuesta de código para resolver ecuaciones de estado cúbicas ......................
239
15
Listado de Figuras
FIGURA 3.1: DIAGRAMA DE FASES P-T PARA UNA SUSTANCIA PURA. ............................................................................................
38 FIGURA 3.2: DIAGRAMA DE FASE 𝑷 − 𝑽𝒆 PARA UNA SUSTANCIA
PURA. ...................................................................................... 40 FIGURA 3.3: DIAGRAMA DE FASES 𝒑
− 𝑻. ............................................................................................................................. 41
16
FIGURA
3.4: PRESIÓN DE VAPOR PARA UN COMPONENTE PURO CON LA ECUACIÓN CLAUSIUS-CLAPEYRON. ......................................... 44
FIGURA 3.5 CARTA DE COX PARA CALCULAR PRESIONES DE VAPOR DE PARAFINAS NORMALES (GPSA, 1987) ......................................
45
FIGURA 3.6: DIAGRAMA PRESIÓN- TEMPERATURA-COMPOSICIÓN PARA UN SISTEMA DE DOS COMPONENTES. ....................................
47
FIGURA 3.7: DIAGRAMA TERNARIO DE COMPOSICIÓN PARA UNA MEZCLA DE TRES COMPONENTES. ....................................................
48
FIGURA 3.8: ENVOLVENTE DE FASES P-T DE UNA MEZCLA MULTICOMPONENTE (AHMED, 2007). .....................................................
49
FIGURA 3.9: COMPORTAMIENTO PRESIÓN TEMPERATURA DE UNA MEZCLA DE METANO Y N-HEPTANO. .............................................
50
FIGURA 3.10: INTERSECCIONES DEL ETANO, PENTANO Y NONANO SOBRE LAS CARTAS DE COX ..........................................................
52
FIGURA 3.11: DIAGRAMA DE PRESIÓN COMPOSICIÓN PARA UNA MEZCLA DE DOS COMPONENTES CON A-METANO Y B-ETANO. (MCCAIN,
1990). .................................................................................................................................................................. 53
FIGURA 3.12 SOLUCIÓN GRÁFICA DEL EJERCICIO 2.2 ................................................................................................................
54 FIGURA 3.13: SOLUCIÓN GRÁFICA DEL EJERCICIO
2.3. .............................................................................................................. 56
FIGURA 3.14 DIAGRAMA TEMPERATURA-COMPOSICIÓN PARA UNA MEZCLA BINARIA DE METANO Y ETANO (MC CAIN, 1990) ................
57
FIGURA 3.15: SOLUCIÓN DEL EJERCICIO 3.4 ...........................................................................................................................
58 FIGURA 3.16: DIAGRAMA TERNARIO DE METANO, PROPANO Y N-PENTANO (A 𝟓𝟎𝟎 𝐩𝐬𝐢𝐚 Y
𝟏𝟔𝟎°𝐅). .............................................. 60 FIGURA 3.17 UBICACIÓN DE LOS PUNTOS 1, 2 Y 3 DEL EJERCICIO
3.5. ......................................................................................... 61 FIGURA 3.18: COMPOSICIONES DEL GAS EN
EQUILIBRO. ............................................................................................................ 62 FIGURA 3.19: COMPOSICIONES DEL GAS
EN EQUILIBRIO. ........................................................................................................... 62
FIGURA 3.20: DIAGRAMA TERNARIO DE UNA MEZCLA DE METANO, PROPANO Y N-PENTANO A CONDICIONES DE 𝟏𝟔𝟎°𝐅 Y 𝟏𝟓𝟎𝟎 𝐩𝐬𝐢𝐚. .
65
FIGURA 3.21: UBICACIÓN DE LOS PUNTOS 1, 2 Y 3 DEL EJERCICIO 3.7. ........................................................................................
66 FIGURA 3.22: COMPOSICIONES DE GAS DEL EJERCICIO
3.7. ....................................................................................................... 66 FIGURA 3.23: COMPOSICIONES DEL LÍQUIDO DEL
EJERCICIO 3.7. ................................................................................................ 67
FIGURA 4.1: CORRELACIÓN DE STANDING Y KATZ PARA LA OBTENCIÓN DEL FACTOR “Z” (H. WHITSON, 1994) ...................................
78
17
FIGURA
FIGURA 5.6: FLUIDO EN MOVIMIENTO A TRAVÉS DE UNA SECCIÓN DE TUBERÍA CON DIFERENTE DIÁMETRO . .......................................
104
FIGURA 5.7 SOLUCIÓN GRÁFICA DEL EJERCICIO 5.2. ...............................................................................................................
107
5.8: DEL EJEMPLO 1, REPRESENTA LA DISTANCIA ENTRE DOS TUBOS MANOMÉTRICOS Y LA DIFERENCIA DE ALTURAS DEBIDO A LA
CAÍDA DE PRESIÓN DE UN FLUIDO LAMINAR VISCOSO. ....................................................................................................
109
FIGURA 7.1. COMPORTAMIENTO DEL 𝑩𝒈, CON RESPECTO A LA PRESIÓN DEL YACIMIENTO Y A TEMPERATURA CONSTANTE. ..................
147
FIGURA 7.2: COMPORTAMIENTO ISOTÉRMICO DEL 𝒄𝒈 EN FUNCIÓN DE LA PRESIÓN DEL YACIMIENTO. ..............................................
148
FIGURA 7.3: COMPRESIBILIDAD PSEUDOREDUCIDA PARA GASES NATURALES (MCCAIN, 1990). ......................................................
149
FIGURA 7.4: COMPORTAMIENTO DE LA VISCOSIDAD CON RESPECTO A LA PRESIÓN DEL YACIMIENTO Y A DIFERENTES TEMPERATURAS ......
150
FIGURA 7.5: VISCOSIDAD DE GASES HIDROCARBUROS A PRESIÓN ATMOSFÉRICA (MC CAIN, 1990). ................................................
151
18
FIGURA
FIGURA 7.6: VISCOSIDAD DE GASES HIDROCARBUROS A PRESIÓN ATMOSFÉRICA CONSIDERANDO LA PRESENCIA DE 𝑯𝟐𝑺, 𝑵𝟐 Y 𝑪𝑶𝟐 (MC
CAIN, 1990) ........................................................................................................................................................
152
FIGURA 7.7 RELACIÓN DE VISCOSIDADES PARA GASES NATURALES CON DENSIDADES RELATIVAS DE 0.56 A 0.9 .................................
152
FIGURA 7.8 RELACIÓN DE VISCOSIDADES PARA GASES NATURALES CON DENSIDADES RELATIVAS DE 0.9 A 1.2 ....................................
153
FIGURA 7.9 RELACIÓN DE VISCOSIDADES PARA GASES NATURALES CON DENSIDADES RELATIVAS DE 1.2 A 1.5. ...................................
153
FIGURA 7.10 RELACIÓN DE VISCOSIDADES PARA GASES NATURALES CON DENSIDADES RELATIVAS DE 1.5 A 1.7 ..................................
154
FIGURA 7.11 COMPORTAMIENTO DEL 𝑩𝑶: EN UN ACEITE NEGRO Y EN FUNCIÓN DE LA PRESIÓN DEL YACIMIENTO. ............................
155
FIGURA 7.12: COMPORTAMIENTO DEL 𝑹𝑺 A TEMPERATURA CONSTANTE Y EN FUNCIÓN DE LA PRESIÓN DEL YACIMIENTO PARA UN ACEITE
NEGRO. ................................................................................................................................................................
156
FIGURA 7.13: COMPORTAMIENTO DEL FACTOR DE VOLUMEN TOTAL DE FORMACIÓN CON RESPECTO A LA PRESIÓN (DANESH, 1998). ....
157
FIGURA 7.14: COMPORTAMIENTO DEL 𝒄𝒐 EN FUNCIÓN DE LA PRESIÓN Y A TEMPERATURA CONSTANTE, CUANDO LA PRESIÓN DEL
YACIMIENTO ESTÁ POR ENCIMA DE LA PRESIÓN DE BURBUJA. ...........................................................................................
158
FIGURA 7.15: COMPORTAMIENTO DEL 𝑪𝑶, EN FUNCIÓN DE LA PRESIÓN Y A TEMPERATURA CONSTANTE, PARA PRESIONES DE YACIMIENTO
MENORES A LA PRESIÓN DE BURBUJA. .........................................................................................................................
159
FIGURA 7.16: COMPORTAMIENTO DE LA VISCOSIDAD DEL ACEITE CON RESPECTO A LA PRESIÓN Y A TEMPERATURA CONSTANTE . ............
159
FIGURA 7.17: COMPORTAMIENTO DEL 𝑩𝒘 CON RESPECTO A LA PRESIÓN. .................................................................................
161
FIGURA 7.18: COMPORTAMIENTO DEL 𝑹𝒔𝒘 CON RESPECTO A LA PRESIÓN DEL YACIMIENTO. ........................................................
162
FIGURA 7.19: COMPORTAMIENTO TÍPICO DE LA VISCOSIDAD CONTRA LA PRESIÓN A TEMPERATURA CONSTANTE. ...............................
163
19
FIGURA
FIGURA 8.1: OBTENCIÓN DE UNA MUESTRA DE FONDO, A) CÁMARA DE MUESTREO B) DISTRIBUCIÓN TÍPICA DE PRESIONES (DAKE, 1978).
..........................................................................................................................................................................
179 FIGURA 8.2: OBTENCIÓN DE UNA MUESTRA EN SUPERFICIE (DAKE,
1978). ............................................................................... 180
FIGURA 9.10: DIAGRAMA DE FASE DE PRESIÓN CONTRA TEMPERATURA PARA UN YACIMIENTO DE ACEITE VOLÁTIL. .............................
218
FIGURA 10.3: DIAGRAMA DE FASE DE UN GAS DE CONDENSACIÓN RETROGRADA (CLARK, 1969). ...................................................
219
FIGURA 10.4: DIAGRAMA DE FASE DE PRESIÓN CONTRA TEMPERATURA PARA UN YACIMIENTO DE GAS. HÚMEDO (CLARK, 1969). ........
220
FIGURA 10.5: DIAGRAMA DE FASES PARA UN YACIMIENTO DE GAS SECO (CLARK, 1969). ..............................................................
220
FIGURA 12.1 DIAGRAMA DE FASES PARA UN HIDRATO DE GAS PURO, G (GAS), H, (HIDRATO), L1 (AGUA LÍQUIDA), L2 (HIDRATO EN ESTADO
LÍQUIDO). ............................................................................................................................................................. 229
FIGURA 12.2 DIAGRAMA DE FASES PARA UN HIDRATO CON UNA MEZCLA DE GASES. .....................................................................
230 FIGURA 12.3 ESTRUCTURA DE LAS CAVIDADES DE LOS
HIDRATOS. ............................................................................................. 231
20
FIGURA
21
Listado de tablas
22
TABLA 6.16 SOLUCIÓN DEL EJERCICIO 6.9 (CONTINUACIÓN). ................................................................................................... 144
TABLA 7.1 CATIONES Y ANIONES PRESENTES EN EL AGUA DE FORMACIÓN ................................................................................... 160
TABLA 7.2 COMPOSICIÓN DE LA MEZCLA DEL EJERCICIO 7.5 ..................................................................................................... 169
TABLA 7.3 CÁLCULOS DEL EJERCICIO 7.5 .............................................................................................................................. 170
TABLA 7.4 COMPOSICIÓN DEL GAS DEL EJERCICIO 7.7 ............................................................................................................. 171
TABLA 7.5 CALCULO DE LA PRESIÓN Y TEMPERATURA PSEUDOCRÍTICAS DEL EJERCICIO 7.7 ............................................................. 172
TABLA 8.1 RESULTADO DE LA PRUEBA DIFERENCIAL PARA EL EJERCICIO 8.1. ................................................................................ 187
TABLA 8.2 RESULTADOS DE LA PRUEBA POR SEPARACIÓN DIFERENCIAL DEL EJERCICIO 8.2. ............................................................. 188
TABLA 8.3 CÁLCULOS DEL EJERCICIO 8.2. ............................................................................................................................. 189
TABLA 8.4 VALORES DE BO PARA EL EJERCICIO 8.3. ............................................................................................................... 190
TABLA 8.5 VALORES DE RS PARA EL EJERCICIO 8.4. ................................................................................................................ 190
TABLA 9.1 CONTANTES PARA LAS ESTIMACIONES DE TC Y PC. ................................................................................................... 192
TABLA 9.2 CONSTANTES PAR LA ESTIMACIÓN DE LA 𝒑𝒃. ......................................................................................................... 196
TABLA 9.3 CONSTANTES PARA LA ESTIMACIÓN DE BO. ............................................................................................................ 197
TABLA 10.1 COMPOSICIONES DE LOS FLUIDOS PETROLEROS ..................................................................................................... 210
TABLA 10.2 CRITERIOS PARA LA CLASIFICACIÓN DE HIDROCARBUROS SEGÚN MCCAIN. .................................................................. 215
TABLA 10.3 CRITERIOS DE CLASIFICACIÓN DE LOS FLUIDOS PETROLEROS SEGÚN MÉNDEZ. .............................................................. 215
TABLA 10.4 COMPOSICIONES DE LOS FLUIDOS PETROLEROS ..................................................................................................... 215
TABLA 11.1 COMPOSICIÓN DE LA MEZCLA DEL EJERCICIO 11.1 ................................................................................................. 226
TABLA 11.2 SOLUCIÓN DEL EJERCICIO 11.1 .......................................................................................................................... 226
TABLA 11.3 RESULTADOS DEL EJERCICIO 11.2 ...................................................................................................................... 227
23
Capítulo 1.
Introducción
Los fluidos petroleros se encuentran constituidos por componentes hidrocarburos y
no hidrocarburos, sus propiedades no son constantes y dependen de las condiciones de
presión, temperatura y composición del sistema, estas propiedades son muy
importantes puesto que dominarán los esquemas de explotación que se seguirán a fin
de optimizar la explotación de los mismos.
Las propiedades de los fluidos petroleros y su comportamiento son la base para las
demás disciplinas de la ingeniería petrolera, pues dominarán el comportamiento del
yacimiento y serán fundamentales para la interpretación y diseño de pruebas que
pueden proveer mucha información del yacimiento.
Estos apuntes son un esfuerzo por generar material que puede ser usado por los
alumnos de Ingeniería Petrolera que cursan el nuevo plan de estudios, además de
describir y repasar los temas fundamentales se presentan una serie de ejercicios que
pueden ayudar a repasar los conocimientos aprendidos en el aula de clase, estos
apuntes tienen además el objetivo de unificar los temas que deben de estudiarse a fin
de homogeneizar el conocimiento entre los alumnos.
En el capítulo 2, se describen generalidades básicas para el resto de los capítulos
como por ejemplo los sistemas de unidades: internacional, inglés, etc., expresiones con
24
inconsistencia de unidades y algunos conceptos básicos en la ingeniería de
yacimientos.
En el capítulo 3, se habla de sistema en una y varias fases, con uno, dos y tres
componentes, y se presentan los diagramas de fase de los mismos. Se obtienen
presiones de vapor de los componentes hidrocarburos con las Cartas de Cox, los
factores de compresibilidad “z” y se realizan cálculos de las composiciones de los
componentes hidrocarburos en una mezcla, dados diferentes diagramas de
composición - temperatura y composición – presión.
Para el capítulo 4, se estudiaron los gases ideales y reales, sus propiedades,
comportamiento, analizando ecuaciones básicas como la de estado entre otras, se
realizan conversiones fundamentales como es la de fracción mol a fracción peso, etc.
También se analizó el comportamiento de los mismos en mezclas y cuando contenían
heptanos y más pesados.
Posteriormente para el capítulo 5, se estudian las propiedades y comportamiento de
los líquidos ideales y reales, de igual manera se estudian cuando estos permanecen
estáticos y cuando están en movimiento.
En el capítulo 6 se analizan las ecuaciones de estado más comunes como las de Van
der Waals, Peng Robinson, entre otras, además de reglas de mezclado para encontrar
las propiedades de los líquidos y gases cuando están formando parte de una mezcla, el
objetivo general es encontrar el factor “z”, que nos ayudará principalmente a calcular la
presión de vapor de la sustancia.
Por otra parte en el capítulo 7, se habla acerca de las propiedades de los fluidos que
conforman los yacimientos, como son el gas, al aceite y el agua connata, entre estas se
mencionan el factor de volumen, la densidad, el coeficiente de viscosidad, la relación de
solubilidad entre otras.
En el capítulo 8 se habla acerca de las diferentes pruebas PVT que se aplican a las
muestras de hidrocarburos para saber con más exactitud ciertas propiedades de los
mismos, son de suma importancia debido a que de sus valores depende el diseño
correcto de la logística de la extracción de dichos hidrocarburos que a su vez también
está basada bajo la norma API RP-44 aquí mencionada.
Así mismo en el capítulo 9 se estudian las diferentes correlaciones que proporcionan
los valores de las propiedades de los hidrocarburos, aunque existen diversidad de ellas
25
en la literatura, las más usuales y también tratadas en este capítulo son la de Vázquez y
Beggs.
En el Capítulo 10 se habla acerca de las clasificaciones de diferentes autores que
hacen a los yacimientos de acuerdo a los fluidos que contienen.
En el capítulo 11 se habla del equilibro líquido – vapor, es de interés para la
ingeniería, se observan las variables principales que determinan este estado como la
temperatura, la presión y la composición del sistema.
Finalmente en el capítulo 12, se habla de los hidratos de gas, de sus propiedades y
comportamiento cuando están en una o varias fases, de los tipos de hidratos que hay y
se hace un modelamiento de la formación de hidratos.
Capítulo 2.
26
Fisico-química y termodinámica de hidrocarburos.
conversión, uno muy común en el cálculo de las propiedades de los fluidos petroleros
27
es el valor de la constate de gravitacional 𝑔𝑐 = 32.174 𝑙𝑏𝑚 − 𝑓𝑡
𝑙𝑏𝑓− 𝑠 2 , este valor no
En este sistema es común utilizar dos tipos de libras, las libras-masa (𝑙𝑏𝑚) para
medir masa y las libras-fuerza (𝑙𝑏𝑓) para medir fuerza, por lo que se trata de un sistema
de unidades inconsistente debido a la necesidad de realizar conversiones al operar las
magnitudes.
2.1.3 Sistema métrico de unidades.
Este incluye a los sistemas MKS, CGS y algunas unidades del sistema inglés, por
sus unidades de longitud, masa y tiempo, algunas de ellas están indicadas en la Tabla
2.3.
Tabla 2.3 Unidades del sistema métrico
MAGNITUD Sistema MKS Sistema CGS Sistema Inglés LPS
Longitud metro (m) Centímetro (cm) Pie (ft)
Masa Kilogramo (kg) Gramo (gr) Libra (lb)
Tiempo Segundo (s) Segundo (s) Segundo (s)
Materia kg- kmol gr-mol lb-mol
Temperatura °Celcius (°C) °Celcius (°C) °Fahrenheit (°F)
Para este sistema la materia únicamente se mide en unidades de masa, no se usan
las (𝑙𝑏𝑓), únicamente las (𝑙𝑏𝑚).
• Presión hidrostática
: presión ( ) 𝑙𝑏
: gravedad 𝑙𝑏 𝑓𝑡 𝑘𝑔 𝑚 : altura ( )( )
𝑃 𝑃𝑎 , ( 2 ) 𝑔 ( 2 ), ( 2) ℎ 𝑓𝑡 𝑚
𝑝𝑔 𝑎𝑏𝑠 𝑠 𝑠
𝑓𝑡3 𝑚3
Peso específico
𝜌𝑔 ......................................................
(2.2)
𝑔𝑐
Esfuerzo cortante
𝜇 𝑑𝑣 ................................................... (2.3)
29
𝛾=
𝜏= ( )
𝑔𝑐 𝑑𝑦
Donde
• Energía cinética
𝑚𝑣2 ............................ .............. (2.4)
𝐸𝑐 =
2𝑔𝑐
Donde
𝐸𝑐: Energía cinética (𝐽), (𝑓𝑡 − 𝑙𝑏𝑓) 𝑚: masa (𝑙𝑏𝑚), (𝑘𝑔)
• Energía potencial
𝑚𝑜𝑙
Mol Es el peso molecular de una sustancia -------- 𝑚𝑜𝑙
Presión Medida de desorden ( por el choque 𝐹 𝑙𝑏𝑓2, 𝑚𝑁2 entre las
moléculas, paredes del 𝑃=𝐴 𝑝𝑔 recipiente, etc.) en una
sustancia
Presión Es la suma de la presión manométrica 𝑃𝑎𝑏𝑠 = 𝑃𝑚𝑎𝑛 + 𝑃𝑎𝑏𝑠 𝑙𝑏𝑓2, 𝑚𝑁2
31
manométrica 𝑝𝑔
𝐶𝑜 = − 𝑉𝑜 (
1 𝛿𝑉𝑜 ) 𝑇
𝛿𝑃
32
Propiedades Propiedades que no dependen de la ____
Intensivas cantidad de materia en una sustancia o ______
cuerpo (temperatura, densidad,
viscosidad, etc).
Ejercicio 2.1
Completar la Tabla 2.6 de equivalencias.
Tabla 2.6 del Ejercicio 2.1 de equivalencias
Kg Slug 𝒍𝒃𝒎
1.- 4
2.- 10
3.- 8
Solución:
1 𝑙𝑏𝑚 1000 𝑔𝑟
10 𝑠𝑙𝑢𝑔 = 0.14𝑘𝑔 ( )( ) = 𝟎. 𝟑𝟎𝟖 𝒍𝒃𝒎
454 𝑔𝑟 1 𝑘𝑔
3.-
33
Ejercicio 2.2
Obtener la masa del gas propano que se encuentra confinado en el cilindro de
un émbolo a 𝟔𝟎°𝑪 y 𝟏𝟐𝟎 𝒂𝒕𝒎 y 𝟖𝟎 𝒍𝒕𝒔 de volumen. Exprese el resultado en 𝒍𝒃, 𝒌𝒈 y
𝒐𝒛. Considere el peso molecular del propano, 𝑪𝟑𝑯𝟖, igual a 𝟒𝟒. 𝟎𝟗𝟕 𝒍𝒃𝒎/𝒍𝒃𝒎 −
𝒎𝒐𝒍, además de que se comporta como un gas ideal y la constante universal de
los
𝒍𝒃 𝒇𝒕𝟑
𝒑𝒈
𝟐𝒂𝒃𝒔 gases 𝑹 = 𝟏𝟎. 𝟕𝟑𝟐.
𝒍𝒃𝒎 𝒎𝒐𝒍 °𝑹
Solución
𝑝𝑉 ∗ 𝑃𝑀
𝑚=
𝑅𝑇
El propano puede ser considerado como un gas ideal por lo que 𝑧 = 1, luego se
realizan las conversiones de unidades necesarias y se sustituyen en la ecuación
anterior.
Para la temperatura
°𝐹 = 1.8(°𝐶) + 32 = 1.8(60) + 32 =140
°𝑅 = °𝐹 + 459.69 = 140 + 459.69 = 𝟓𝟗𝟗. 𝟔𝟗 °𝐑. Para
la presión
𝑙𝑏
𝑝𝑔
14.695 2𝑎𝑏𝑠 𝒍𝒃
Para el volumen
𝑉 = 80 𝑙𝑡𝑠 ( 1𝑚3 ) (3.28 𝑓𝑡)3 = 𝟐. 𝟖𝟐𝟑 𝒇𝒕𝟑.
1 000 𝑙𝑡𝑠 1𝑚
Sustituyendo los valores para determinar la masa.
34
𝑙𝑏𝑚
1763.4 𝑝𝑔𝑙𝑏2𝑎𝑏𝑠 ∗ 2.823 𝑓𝑡3 ∗ 44.097 𝑙𝑏𝑚 𝑚𝑜𝑙
𝑚 == 𝟑𝟒. 𝟏𝟎𝟖𝟓 𝒍𝒃𝒎
𝑙𝑏 3
2
𝑎𝑏𝑠 𝑓𝑡
𝑝𝑔
10.732 ∗ 599.69 °R
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
Para las conversiones de la masa en 𝑙𝑏𝑚 a 𝑜𝑧 y 𝑘𝑔
16 𝑜𝑧
34.1085 𝑙𝑏𝑚 ( ) = 𝟓𝟒𝟓. 𝟕𝟑𝟔 𝒐𝒛
1 𝑙𝑏
0.4535𝑘𝑔
Ejercicio 2.3
En una plataforma petrolera se tienen 𝟐 tuberías concéntricas de 𝟏𝟐 𝒇𝒕 de altura
𝒑𝒈𝟐
𝒇𝒕
y la de la menor 𝟑 𝒑𝒈𝟐.
a) Calcular la presión hidrostática en el fondo de la tubería de base menor,
resultado en 𝒎𝒎𝑯𝒈.
Solución:
35
las unidades de 𝑙𝑏𝑚 a 𝑙𝑏𝑓.
𝑙𝑏 𝑓𝑡
𝜌𝑔ℎ (72 𝑓𝑡3)(32.174𝑠2)(12 𝑓𝑡) 𝒍𝒃𝒇
b) Como ambas tuberías están llenas a la misma altura su 𝑃ℎ es igual solo debe
presión atmosférica.
𝑇𝑢𝑏𝑚𝑎𝑦𝑜𝑟 = 𝟏. 𝟎𝟒 𝒂𝒕𝒎.
Ejercicio 2.4
𝒎
𝒌𝒈
𝟑 ). Indicar los resultados en 𝒇𝒕𝟑 , 𝒄𝒎𝟑 y 𝒈𝒂𝒍.
36
𝒍𝒃𝒎 𝒈𝒓 𝒍𝒃𝒎
Solución:
𝑚3 1 𝑓𝑡 3
0.4535 𝑘𝑔 𝑓𝑡3
1.26 𝑥10−3 ( ) ( ) = 0.0202
𝑘𝑔 0.3046𝑚 1 𝑙𝑏𝑚 𝑙𝑏𝑚
𝑚3 100 𝑐𝑚 3 1 𝑘𝑔 𝑐𝑚3
1.26 𝑥10−3 ( )( )=
1.26 𝑘𝑔 1𝑚 1000 𝑔𝑟 𝑔𝑟
𝑚3 0.4535 𝑘𝑔 1 𝑔𝑎𝑙 𝑔𝑎𝑙
1.26 𝑥10−3 𝑘𝑔 ( 1 𝑙𝑏𝑚 )
810𝑚3
Transformando unidades
−3
𝑚3 ( 1 𝑓𝑡 )3 (0.4535 𝑘𝑔) = 0.0197 𝑓𝑡3
1.23 𝑥10
𝑘𝑔 0.3046𝑚 1 𝑙𝑏𝑚 𝑙𝑏𝑚
𝑚3 100 𝑐𝑚 3 1 𝑘𝑔 𝑐𝑚3
1.23 𝑥10−3 ( )( )=
1.23 𝑘𝑔 1𝑚 1000 𝑔𝑟 𝑔𝑟
1.23 𝑥10−3 𝑚 𝑘𝑔3 (0.45351 𝑙𝑏𝑚 𝑘𝑔) (0 .0037851 𝑔𝑎𝑙 𝑚3) = 0.1473 𝑙𝑏𝑚 𝑔𝑎𝑙.
Ejercicio 2.5
Completar la Tabla 2.7.
37
Tabla 2.7 Ejercicio para completar
P m g
4 lb 𝒇𝒕
28.4
𝒔𝟐
𝒎
386.4 𝒍𝒃𝒇 9.81
𝒔𝟐
Solución:
Primera fila:
𝑓𝑡
𝑃 = 𝑚𝑔 = 4 𝑙𝑏 × 28.4 𝑠2 = 113.6 𝑙𝑏𝑓
Segunda fila, es necesario despejar el valor de la masa de la ecuación de peso y
s
𝑚 3.28 𝑓𝑡 𝑓𝑡
9.81 2 () = 32.1768 2,
𝑠 1𝑚 𝑠
sustituyendo
𝑃 386.4 𝑙𝑏
𝑚= = = 𝟏𝟐 𝒍𝒃.
𝑔 𝑓𝑡
32.1768 𝑠2
su energía
𝒔
cinética, expresar el resultado en 𝒇𝒕 y 𝒍𝒃𝒇.
Solución:
38
Sustituyendo en la ecuación 2.4 con el factor 𝑔𝑐:
𝑚𝑣2
𝐸𝑐 =
2𝑔𝑐
𝑓𝑡
(12 125 𝑙𝑏𝑚) (1800 𝑠 )
𝐸𝑐 =
𝑙𝑏𝑚 − 𝑓𝑡
2 (32.174 𝑙𝑏𝑓 − 𝑠2 )
Ejercicio 2.7
Solución:
𝑑 = 𝑣 = 88 𝑐𝑚3 = 𝟎. 𝟕𝟔𝟏𝟑 𝒄𝒎 𝟑 .
Realizando las conversiones para expresar los resultados en las unidades requeridas
3
𝑔𝑟 1 𝑘𝑔 100 𝑐𝑚 𝑘𝑔
0.7613 𝑐𝑚3 (1000 𝑔𝑟) ( 1
𝑚) = 761.3 𝑚 3 ,
3 3
𝑔𝑟 1 𝑙𝑏 100 𝑐𝑚 0.003785 𝑚 𝑙𝑏
0.7613 𝑐𝑚 (
3 453.593 𝑔𝑟 ) ( 1𝑚) ( 1 𝑔𝑎𝑙 ) = 635 𝑔𝑎𝑙 ,
3
𝑔𝑟 1 𝑙𝑏 30.46 𝑐𝑚 𝑙𝑏
0.7613 𝑐𝑚3 (453.593 𝑔𝑟) ( 1 𝑓𝑡 ) = 47.43 𝑓𝑡 3 ,
𝑚
1 3
−3
𝑚3. volumen
3 𝑙𝑖𝑡 () = 3 𝑥 10
1 000 𝑙𝑖𝑡
39
3
𝑔𝑟 𝑙𝑏 5.6146 𝑓𝑡 1 𝑠𝑙𝑢𝑔 𝑠𝑙𝑢𝑔
0.7613 𝑐𝑚 3 = 47.43 𝑓𝑡 3 (
) (32.17 𝑙𝑏𝑚) = 8.2774 𝑏𝑙 .
Ejercicio 2.8
Un tanque de 𝟑 𝒍𝒊𝒕𝒓𝒐𝒔 contiene oxígeno a 𝟐𝟎°𝑪 y 𝟐𝟓 𝒙 𝟏𝟎𝟓 𝑷𝒂 de presión.
Calcular la masa del oxígeno en el tanque, considere que la masa molecular del
oxígeno es de 𝟑𝟐 𝒌𝒈
y la presión atmosférica es de 𝟏 𝒙 𝟏𝟎𝟓𝑷𝒂. Expresar el
resultado en 𝒌𝒈.
𝑲𝒎𝒐𝒍
Solución:
𝑝𝑉 ∗ 𝑃𝑀
𝑚= .
𝑅𝑇
Sustituyendo:
𝑘𝑔
(26 𝑥 105𝑃𝑎)(3 𝑥 10−3 𝑚3) (32 )
𝐾𝑚𝑜𝑙
𝑚 = = 𝟎. 𝟏𝟎 𝒌𝒈.
𝐽
(8314 ) (293 𝐾)
𝐾𝑚𝑜𝑙 𝐾
Ejercicio 2.9
40
Determine el volumen de 𝟖 𝐠 de helio (𝐌 = 𝟒 𝐤𝐠/𝐦𝐨𝐥) a 𝟏𝟓 °𝐂 y 𝟒𝟖𝟎 𝐦𝐦𝐇𝐠.
Solución:
se mueve un tramo
𝒔
Solución:
41
𝐹𝑓 = = 0.014 𝑁.
Para la 𝑙𝑏𝑓:
1 𝑙𝑏𝑓 −𝟑
0.014 𝑁 ( 𝒍𝒃𝒇.
) = 𝟑. 𝟏𝟒𝟕𝟑 𝒙 𝟏𝟎
4.4482 𝑁
Ejercicio 2.11
Un objeto situado a 𝟏𝟎 𝒎 de altura y cuya masa es de 𝟑𝟐𝟎 𝒍𝒊𝒃𝒓𝒂𝒔 fue elevado
hasta una altura de 𝟖𝟎 𝒎.
a) Calcular la energía potencial de la posición inicial del objeto.
posición de objeto.
Solución:
𝑚𝑔ℎ
𝐸𝑝 =
𝑔𝑐
𝑓𝑡
(320 𝑙𝑏𝑚)(32.174 𝑠2)(32.82 𝑓𝑡)
𝐸𝑝 = 𝑙𝑏𝑚 − 𝑓𝑡
32.174 𝑙𝑏𝑓 − 𝑠2
𝑓𝑡
𝐸𝑝2 = 𝑚𝑔𝑔 ℎ = (320 𝑙𝑏𝑚)(32.174𝑙𝑏𝑚𝑠2−)(262𝑓𝑡.63 𝑓𝑡) = 84 041.6 𝑙𝑏𝑓 𝑓𝑡.
42
𝑐 32.174 𝑙𝑏𝑓 − 𝑠2
Ejercicio 2.12.
A 𝟑𝟐°𝑪 y 𝟐. 𝟏𝟎 𝒌𝒈
𝟐 , el volumen específico de cierto gas es de 𝟎. 𝟕𝟏 𝒎 𝟑
.
Determinar
𝒄𝒎 𝒌𝒈
temperatura en °𝐾.
𝑅𝑇
𝑘𝑔 𝑚3
(2.10𝑥 104
𝑃 𝑃 𝑉𝑠 𝑚2) (0.71 𝑘𝑔) 𝑚
𝑅= = = = 68.8 .
𝑤𝑇 𝑇 305 °𝐾 °𝐾
Para la densidad:
1
𝑤 𝑣𝑠
1 𝑚4 𝑼𝑻𝑴
𝜌 = 𝑔 = 𝑔 = 𝑣𝑠𝑔 = 0.71 𝑚𝑘𝑔3 ∗ 9.81 𝑠𝑚2 = 0.1436 𝑘𝑔𝑠2 =
𝟎. 𝟏𝟒𝟑𝟔 𝒎𝟑 .
de presión,
𝒑𝒈 𝒂𝒃𝒔
43
considerando que el butano tiene un factor de compresibilidad 𝒛 = 𝟎. 𝟓 y una 𝑴
= 𝟓𝟖. 𝟏𝟐 𝒈𝒓/𝒎𝒐𝒍, calcular la masa y el volumen específico del butano en 𝒌𝒈 y en
𝒍𝒃.
Solución:
Realizando conversiones.
De temperatura:
270°𝐾 𝑎 °𝑅
°𝐶 = °𝐾 − 273.16 = 290 − 273.16 = 16.84°𝐶 °𝐹 =
1.8(°𝐶) + 32 = 1.8(16.84°𝐶) + 32 = 62.312°𝐹
°𝑅 = °𝐹 + 459.69 = 62.312°𝐹 + 459.69 = 522.00°𝑅.
De volumen:
18 𝑔𝑎𝑙 𝑎 𝑙𝑏𝑚,
𝑙𝑏𝑚
( 1200 𝑙𝑏2𝑎𝑏𝑠)(2.4𝑓𝑡3) (0.1281 𝑙𝑏𝑚
𝑚𝑜𝑙) 𝑝𝑔
𝑚 = = 0.1317 𝑙𝑏𝑚.
2
𝑙𝑏 3 𝑎𝑏𝑠
𝑓𝑡 𝑝𝑔
(0.5) (10.732 ) (522.002 °𝑅)
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
44
Para los kg:
0.4535 𝑘𝑔
𝑚 = 0.1317 𝑙𝑏𝑚 ( ) = 𝟎. 𝟎𝟓𝟗𝟕 𝒌𝒈.
1 𝑙𝑏𝑚
Capítulo 3.
Comportamiento de fases.
Sistema: Parte del espacio aislada por fronteras reales o hipotéticas, cuyo objeto
principal es el estudio del mismo así como las interacciones entre los elementos que lo
conforman.
Fase: Es una parte en un sistema con consistencia homogénea y distinta a las
demás fracciones que lo conforman, por lo que está separada por superficies
fronterizas.
Las diferentes fases se identifican por su volumen o densidad.
Comportamiento de fase: Es el cambio que existe entre las fases al modificar sus
condiciones de temperatura y presión para las cuales pueden existir.
Diagrama de fase: Gráficas obtenidas de manera experimental y matemática del
comportamiento de las fases y los límites físicos para la existencia de cada una de
ellas.
45
Regla de la fase (Gibbs): Indica el máximo número posible de fases en equilibrio
que puedan coexistir y el número de componentes presentes, aplica solamente a un
sistema en equilibrio estable.
𝑭=𝑪−𝑷∓𝟐 .............................................. (2.1)
Donde:
𝐹: Grados de libertad
𝑃: Número de fases
Nombre del Masa Volumen Presión Temperatura Valor Factor compuesto molecular crítico
crítica Crítica °𝑭 de 𝒁𝑹𝑨 acéntrico
ω
(𝒍𝒃𝒎 ) (𝒍𝒃𝒎 𝟑 ) ( 𝒍𝒃𝟐 )
𝒎𝒐𝒍 𝒇𝒕 𝒑𝒈 𝒂𝒃𝒔
48
Figura 3.2: Diagrama de fase 𝑷 − 𝑽𝒆 para una sustancia pura.
La Fig. 3.2 muestra que mientras la temperatura incrementa, la presión disminuye
hasta desaparecer en el punto crítico, este se encuentra sobre una isoterma,
especialmente para un componente puro. El punto crítico se obtiene con las
ecuaciones 3.2 y 3.3,
( 𝜕𝑝) = 0, . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (3.2)
𝜕𝑉 𝑇𝑐
( 𝜕𝑝2
2
) = 0. ........... ................................... (3.3)
𝜕𝑉 𝑇𝑐
49
Figura 3.3: Diagrama de fases 𝒑 − 𝑻.
50
𝜌= 𝑙
[1+(1−𝑡𝑟)
2/7
]
𝑅𝑇𝑐𝑍𝑐
𝑀 𝑝𝑐𝑉𝑐 ................................................ (3.6)
𝑍𝑐 = .
𝑅 𝑇𝑐
𝑇 ................................................ (3.7)
𝑇𝑟 = .
𝑇𝐶
Donde:
.
𝜌𝑙: densidad del líquido saturado (𝑙𝑏𝑚𝑓𝑡3 ). 𝑀: Peso molecular (𝑙𝑏𝑚 𝑙𝑏𝑚
). 𝑃𝑐:
𝑚𝑜𝑙
𝑓𝑡
3
𝜌𝑙 = [1+(1−𝑡𝑟)2/7]
51
𝑅𝑇𝑐𝑍𝑅𝐴
3.2.8 Métodos para determinar la presión de vapor para una sustancia pura.
Donde:
𝑇1,2: Condiciones diferentes de temperatura.
𝐿𝑣: Calor de vaporización de una mol de líquido.
𝑃1,2: Condiciones diferentes de presión.
52
Figura 3.4: Presión de vapor para un componente puro con la ecuación Clausius-Clapeyron.
53
Figura 3.5 Carta de Cox para calcular presiones de vapor de parafinas normales (GPSA,
1987)
54
𝑝𝑣 = 𝑝𝑐𝐸𝑋𝑃 (𝐴 + 𝜔𝐵). ... .................................... (3.11)
Donde: 𝐴 = 5.92714 −
𝑇𝑟
6.09648 𝑟 𝑟
6 .............. (3.12) −
1.2886 ln ( )
15.6875 (3.13) 𝑇
𝐵 = 15.2518 −
+ 0.16934 𝑇−, 13.472 𝑙𝑛(𝑇𝑟) + 0.4357𝑇𝑟6 y ...............
𝑇𝑟
𝑃𝑣 ......................................... (2.14)
𝜔 ∶ −𝐿𝑜𝑔 (
𝑃𝑐) − 1.
= 0.7
Los valores del factor acéntrico para algunas sustancias se encuentran en la Tabla
3.3.
Los cambios de fase para un sistema de dos componentes pueden ser mostrados
de mejor manera mediante diagramas presión, temperatura y composición, la Fig. 3.6
es una representación tridimensional de los anteriores para un sistema de dos
componentes.
55
Figura 3.6: Diagrama Presión- Temperatura-Composición para un sistema de dos
componentes.
56
Figura 3.7: Diagrama ternario de composición para una mezcla de tres componentes.
Donde:
1.- Es el 100% mol del componente puro B.
2. Es una mezcla de dos componentes de 30% mol del componente A y 70% mol del
componente C.
3. Representa una mezcla de tres componentes (50% mol del componente A, 30%
mol del componente B y 20% mol del componente C).
7. La composición en el punto 7 es 50% mol del componente B, 15% mol del
componente A y 35% mol del componente C.
El componente puro B no se encuentra presente.
La composición se presenta en fracción mol o por ciento mol, manteniéndose la
presión y temperatura del sistema constantes.
57
Figura 3.8: Envolvente de fases p-T de una mezcla multicomponente (Ahmed, 2007).
Se observa que a una temperatura dada, las presiones de burbuja y de rocío son
distintas y forman una envolvente de fases en el diagrama p-T, contrario a un
componente puro. Las líneas de calidad se definen como las curvas que contienen a
todos los puntos a los que existen las dos fases unidas de la mezcla en equilibrio y a
ciertas proporciones referidas al líquido, por lo que de acuerdo a lo anterior, las curvas
de burbuja y de rocío son las líneas de calidad del 100 y 0 %, respectivamente. La
forma de la envolvente de fases depende de la concentración y tipos de componentes
en la mezcla. En la Fig. 3.9 se observa a la cricondenbara en la línea aa y se define
como la máxima presión a la cual dos fases pueden coexistir en equilibrio, por encima
de esta, no es posible formar una nueva fase la línea bb se denomina cricondenterma
es decir la máxima temperatura a la cual dos fases pueden coexistir en equilibrio, por
arriba de la cricondenterma no es posible formar una nueva fase. Estas líneas coinciden
con la máxima presión y temperatura a las coexisten las fases a condiciones de
equilibrio.
58
Figura 3.9: Comportamiento presión temperatura de una mezcla de metano y n-heptano.
Ejercicio 3.1
Etano 78 400
i-Pentano 140 120
Nonano 200 10
Solución:
Se requiere conocer la presión de vapor de cada componente, misma que puede ser
obtenida de la Carta de Cox de la Fig. 3.5, sobre la misma se selecciona el compuesto
y se interseca con una línea isotérmica, el valor de esta en presión determinará el
estado de agregación.
La Fig. 3.11 muestra los resultados de las intersecciones de cada componente
a) Etano
b) Pentano
líquido. c) Nonano
60
presión es de 10 (𝒑𝒈𝒍𝒃𝟐𝒂𝒃𝒔), como la 𝑃 > 𝑃𝑣, por lo tanto el nonano se encuentra en
estado líquido.
Figura 3.10: Intersecciones del etano, pentano y nonano sobre las cartas de Cox
Ejercicio 3.2
Calcular las composiciones y cantidades de gas y líquido formado cuando 𝟏𝟎
𝒍𝒃𝒎 −mol de una mezcla formada de 𝟑𝟓 % mol del componente A (metano) y
𝟔𝟓 % del componente B (etano) se lleva a un estado de equilibrio a −𝟑𝟐 °𝑭 y
𝒍𝒃
𝟕𝟓𝟎 𝟐 . Emplear el diagrama de fases presión-composición de la Fig. 3.11.
𝒑𝒈 𝒂𝒃𝒔
61
Figura 3.11: Diagrama de presión composición para una mezcla de dos componentes con A-
Metano y B-Etano. (McCain, 1990).
Solución:
curva de los puntos de burbuja este será el punto 2 y finalmente otra línea que toque a
la curva de los puntos de rocío y este será el punto 3. Se leen las composiciónes
62
trazando líneas perpendiculares a la línea formada anteriormente en los puntos 2 y 3 y
que intersequen el eje de composición, el punto 2 corresponde a la composición del
metano y el punto 3 a la del etano.
El procedimiento descrito anteriormente se encuentra graficamente en la Fig. 3.13.
Con la longitud de las líneas de unión se calculan las cantidades de líquido y gas:
12̅̅̅ 35 − 24 11 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑔𝑎𝑠
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑔𝑎𝑠 = = = = 0.57 ( )
23̅̅ 43 − 24 19 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙
13̅̅ 43 − 35 8 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑙í𝑞𝑢𝑖𝑑𝑜
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙í𝑞𝑢𝑖𝑑𝑜 = = = = 0.42 ( )
̅ ̅
23̅ 43 − 24 19 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙 Cantidad de gas =
(0.57)(10 𝑙𝑏𝑚 𝑚𝑜𝑙) = 5.7 𝑙𝑏𝑚 𝑚𝑜𝑙.
Cantidad de líquido= (0.42)(10 𝑙𝑏𝑚 𝑚𝑜𝑙) = 4.2 𝑙𝑏𝑚 𝑚𝑜𝑙.
Ejercicio 3.3
Calcular las composiciones y cantidades de gas y líquido formado cuando 𝟕
𝒍𝒃𝒎 − 𝒎𝒐𝒍 de una mezcla formada de 𝟖𝟓 % mol del componente A (metano) y
𝟏𝟓 % del componente B (etano) se lleva a un estado de equilibrio a −𝟏𝟓𝟎 °𝑭 y
𝒍𝒃
𝟐𝟎𝟎 𝟐 . Emplear el diagrama de fases presión-composición de la Fig. 3.11.
𝒑𝒈 𝒂𝒃𝒔
Solución:
traza una
𝑝𝑔 𝑎𝑏𝑠
línea horizontal que toque a la curva de los puntos de burbuja y otra línea que toque a la
curva de los puntos de rocío y estos son los puntos 2 y 3 respectivamente. Del gráfico
se obtienen las composiciónes trazando líneas perpendiculares a la línea formada
anteriormente en los puntos 2 y 3 que intersequen el eje de composición, el punto 2
corresponde a la composición del metano y el punto 3 a la del etano.
El procedimiento descrito anteriormente se encuentra graficamente en la Fig. 3.13.
64
Figura 3.13: Solución gráfica del ejercicio 2.3.
Con la longitud de las líneas de unión se calculan las cantidades de líquido y gas:
𝑙𝑏𝑚−𝑚𝑜𝑙 𝑔𝑎𝑠
̅12
̅̅̅ 85−56.5 28.5
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑔𝑎𝑠
̅̅̅ =
̅23 == 0.703
97−56.5
= 40.5(
). ̅13
̅̅̅ 97−85
̅23
̅̅̅ = 97−56.5 𝑙𝑏𝑚−𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙
𝑙𝑏𝑚−𝑚𝑜𝑙 𝑙í𝑞𝑢𝑖𝑑𝑜
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙í𝑞𝑢𝑖𝑑𝑜 = == 0.29 ( ).
𝑙𝑏𝑚−𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙
65
Ejercicio 3.4
Calcular las composiciones y las cantidades de líquido y gas formadas en 𝟗 𝒍𝒃𝒎
− 𝒎𝒐𝒍 de una mezcla binaria de gases con 𝟒𝟑% del componente A (metano) y 𝟓𝟑%
del componente B (etano) a condiciones de equilibrio de 𝟓𝟎𝟎 𝒑𝒔𝒊𝒂 y −𝟑𝟎°𝑭 a partir
del diagrama temperatura-composición de la Fig. 3.14.
Solución:
66
las curvas de puntos de burbuja y rocío para conocer las composiciones de las fases
líquida y vapor.
La representación gráfica de la solución se muestra en la Fig. 3.15.
67
Con la longitud de las líneas de unión se calculan las cantidades de líquido y gas.
12̅̅̅ 43 − 29 14 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑔𝑎𝑠
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑔𝑎𝑠 = = = = 0.36 ( ).
23̅̅ 67 − 29 38 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙
13̅̅ 67 − 43 24 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑙í𝑞𝑢𝑖𝑑𝑜
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑔𝑎𝑠 = = = = 0.63 ( ).
̅ ̅
23̅ 67 − 29 38 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙 Cantidad de gas =
(0.36)(7 𝑙𝑏𝑚 𝑚𝑜𝑙) = 𝟐. 𝟓𝟐 𝒍𝒃𝒎 𝒎𝒐𝒍 .
Cantidad de líquido= (0.63)(7 𝑙𝑏𝑚 𝑚𝑜𝑙) = 𝟒. 𝟒𝟏 𝒍𝒃𝒎 𝒎𝒐𝒍.
Ejercicio 3.5
Calcular las composiciones y las cantidades de líquido y gas formadas en
𝟖 𝒍𝒃𝒎 − 𝒎𝒐𝒍 de una mezcla de gases con 𝟔𝟎% de metano, 𝟐𝟑% de propano y 𝟏𝟕%
de n-pentano a condiciones de equilibrio de 𝟓𝟎𝟎 𝒑𝒔𝒊𝒂 y 𝟏𝟔𝟎°𝑭 a partir del
diagrama ternario de la Fig. 3.16.
Solución:
68
Figura 3.16: Diagrama ternario de metano, propano y n-pentano (a 𝟓𝟎𝟎 𝐩𝐬𝐢𝐚 y 𝟏𝟔𝟎°𝐅).
1.- Se ubican sobre el diagrama ternario las composiciones del problema, primero el
63% del metano, después el 23% del propano y finalmente 17% del n-pentano, el punto
de intersección de los tres será el punto 1.
2.- Una vez ubicado el punto 1, se traza una línea que toque a las curvas de puntos
de rocío y puntos de burbuja respectivamente, estos son los puntos 2 y 3. La
representación gráfica de los pasos 1 y 2 de la solución se muestran en la Fig. 3.17.
69
Figura 3.17 Ubicación de los puntos 1, 2 y 3 del ejercicio 3.5.
3.- Luego sobre el punto 2, se trazan líneas en dirección de cada uno de los
compuestos de la mezcla y se obtiene de cada lado del diagrama las composiciones del
gas en equilibro de cada componente.
4.- Del diagrama ternario se obtienen las composiciones para el punto 3, las
composiciones leídas corresponden al líquido en equilibro del sistema. La
representación gráfica de los pasos 3 y 4 se presentan en las Figs. 3.18 y 3.19.
70
Figura 3.18: Composiciones del gas en equilibro.
Con la longitud de las líneas de unión se calculan las cantidades de líquido y gas.
13̅̅̅ 0.66 − 0.63 0.03 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑔𝑎𝑠
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑔𝑎𝑠 = = = = 0.05 ( ).
23̅̅ 0.66 − 0.12 0.54 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙
12̅̅̅ 0.63 − 0.12 0.51 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑙í𝑞𝑢𝑖𝑑𝑜
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙í𝑞𝑢𝑖𝑑𝑜 = = = = 0.94 ( ).
̅ ̅
23̅ 0.66 − 0.12 0.54 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙 Cantidad de gas =
(0.05)(8 𝑙𝑏𝑚 𝑚𝑜𝑙) = 𝟎. 𝟒 𝒍𝒃𝒎 𝒎𝒐𝒍 .
Cantidad de líquido= (0.94)(8 𝑙𝑏𝑚 𝑚𝑜𝑙) = 𝟕. 𝟓𝟐 𝒍𝒃𝒎 𝒎𝒐𝒍.
Ejercicio 3.6
Se tienen 𝟑. 𝟓 𝒎𝒐𝒍 de una mezcla de gases de 𝟏𝟓% de metano, 𝟔𝟓% de
propano y 𝟐𝟎% de n-pentano a condiciones de equilibrio de 𝟓𝟎𝟎 𝒑𝒔𝒊𝒂 y 𝟏𝟔𝟎° 𝑭.
Calcular las composiciones y cantidades de líquido y gas formados, usar el
diagrama de la Fig.
3.16.
Solución:
Posteriormente sobre el punto 2 se trazan líneas que intersequen cada lado del
diagrama, para leer las composiciones del gas en equilibrio y lo mismo con el punto 3,
estas son las lecturas de las composiciones del líquido en equilibrio.
Con los valores de la Tabla 3.7 se determinan las cantidades de líquido y gas de
esta muestra:
La longitud de las líneas de unión proporcionan los datos para calcular las cantidades
de líquido y gas:
13̅̅ 0.30 − 0.15 0.15 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑔𝑎𝑠
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑔𝑎𝑠 = = = = 0.75 ( ).
23̅̅ 0.30 − 0.1 0.2 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙
12̅̅̅ 0.15 − 0.1 0.05 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑙í𝑞𝑢𝑖𝑑𝑜
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙í𝑞𝑢𝑖𝑑𝑜 = = = = 0.25 ( ).
̅ ̅
23̅ 0.30 − 0.1 0.2 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙 Cantidad de gas =
(0.75)(3.5 𝑙𝑏𝑚 𝑚𝑜𝑙) = 𝟐. 𝟔𝟐𝟓 𝒍𝒃𝒎 𝒎𝒐𝒍 .
Cantidad de líquido= (0.25)(3.5 𝑙𝑏𝑚 𝑚𝑜𝑙) = 𝟎. 𝟖𝟕𝟓 𝒍𝒃𝒎 𝒎𝒐𝒍.
Ejercicio 3.7
Se tienen 𝟗 𝒎𝒐𝒍 de una mezcla de gases de 𝟓𝟎% de metano, 𝟏𝟎% de propano y
𝟒𝟎% de n-pentano a condiciones de equilibrio de 𝟏𝟓𝟎𝟎 𝒑𝒔𝒊𝒂 y 𝟏𝟔𝟎° 𝑭. Calcular las
composiciones y cantidades de líquido y gas formados, usar el diagrama de la
Fig.
3.20.
73
Figura 3.20: Diagrama ternario de una mezcla de metano, propano y n-pentano a
condiciones de 𝟏𝟔𝟎°𝐅 y 𝟏𝟓𝟎𝟎 𝐩𝐬𝐢𝐚.
Solución:
74
Figura 3.21: Ubicación de los puntos 1, 2 y 3 del ejercicio 3.7.
75
Figura 3.23: Composiciones del líquido del ejercicio 3.7.
Se registran las lecturas en la Tabla 3.8 y se calculan las cantidades de líquido y gas
en equilibrio formados.
Tabla 3.8 Composición del líquido y gas del ejercicio 3.7.
Component Composición de gas Composición de líquido
e (%mol) (%mol)
Metano 82 40
Propano 8 10
n- 10 50
Pentano
TOTAL 100 100
Con la longitud de las líneas de unión se calculan las cantidades de líquido y gas:
21̅̅ 0.5 − 0.4 0.1 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑔𝑎𝑠
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑔𝑎𝑠 = = = = 0.23 ( ).
23̅̅̅ 0.82 − 0.4 0.42 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙
13̅̅ 0.82 − 0.5 0.32 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑙í𝑞𝑢𝑖𝑑𝑜
𝐹𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙í𝑞𝑢𝑖𝑑𝑜 = = = = 0.76 ( ).
̅
23̅ 0.82 − 0.4 0.42 𝑙𝑏𝑚 − 𝑚𝑜𝑙 𝑚𝑒𝑧𝑐𝑙𝑎 𝑡𝑜𝑡𝑎𝑙 Cantidad de gas =
(0.23)(9 𝑙𝑏𝑚 𝑚𝑜𝑙) = 𝟐. 𝟎𝟕 𝒍𝒃𝒎 𝒎𝒐𝒍 .
Cantidad de líquido= (0.76)(9 𝑙𝑏𝑚 𝑚𝑜𝑙) = 𝟔. 𝟖𝟒 𝒍𝒃𝒎 𝒎𝒐𝒍.
76
Ejercicio 3.8
En un recipiente se tiene n-butano a 𝟐𝟓𝟎 °𝑲, calcular la densidad de dicho gas
utilizando:
a) La correlación de Rackett.
Solución:
(550.6 2𝑎𝑏𝑠)(58.1203𝑚𝑜𝑙)(0.0703𝑙𝑏𝑚)
𝑝𝑔
𝑍𝑐 = 𝑙𝑏 3 = 0.2739.
2
𝑎𝑏𝑠𝑓𝑡
𝑝𝑔
(10.732 )(765.32°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑇𝑟 = 𝑇 𝑐 = 765.31°𝑅 = 0.5879.
77
Se obtiene el valor de 𝑍𝑅𝐴 para el n-butano de la Tabla 3.1 que luego se sustituye
en la ecuación 3.2.
𝑍𝑅𝐴 = 0.2730
Y el valor de 𝑇𝑟 es el que se obtuvo del inciso anterior,
𝑙𝑏𝑚 𝑙𝑏
𝑀 𝑝𝑐 (58.12030 𝑙𝑏𝑚 𝑚𝑜𝑙 )(550.60 𝑝𝑔 𝑎𝑏𝑠 )
2 𝒍𝒃𝒎 𝜌𝑙
2 𝜌𝐿 = 𝑙𝑏
2
= [1+(1−𝑡 𝑟 )7 ] 2 𝑓𝑡 3 [1+(1−0.5879 )7 ] =
( 10.732 ) (761.31°𝑅)(0.2730)
𝟑𝟗. 𝟎𝟗
𝒇𝒕𝟑 .
𝑝𝑔 𝑎𝑏𝑠
𝑍
𝑅𝑇𝑐 𝑅𝐴 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
Ejercicio 3.9
Con la correlación de Lee y Kesler calcular la presión de vapor del metano a
𝟑𝟏𝟓 °𝑲.
Solución:
𝑇𝑟 = 𝑇 𝐶 = 576.36 °𝑅 = 0.9837.
Ejercicio 3.10
Solución:
Ejercicio 3.11
Completar la Tabla 3.9 para determinar el número de grados de libertad.
Solución:
79
Tabla 3.10 Resultados.
80
Capítulo 4.
Comportamiento de los gases
4.1 Definición de una ecuación de estado y su utilidad.
Donde:
𝑓𝑡
𝑙𝑏 𝑝𝑔2𝑙𝑏𝑎𝑏𝑠 3
𝑃 = presión ( )
2 𝑅 = Constante de gases ideales= 10.73 (10.732 )
𝑝𝑔 𝑎𝑏𝑠 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
81
𝑉 = volumen (𝑓𝑡3) 𝑛 =número de moles 𝑇 = temperatura (°𝐹)
Su uso está limitado a un intervalo pequeño de presiones y temperaturas porque las
condiciones anteriormente descritas no siempre se pueden lograr.
Esta ecuación también describe el comportamiento de gases reales, con ciertas
limitaciones. El comportamiento ideal de los gases puede corregirse de las siguientes
dos formas:
1. Mediante un factor multiplicativo denominado factor de desviación 𝑍.
2. En forma aditiva para considerar los efectos moleculares.
82
𝑚 𝑝𝑀 .............................................. (4.6)
𝜌𝑔 = = .
𝑉 𝑅𝑇
4.3 Mezclas de gases ideales.
La Tabla 4.1 resume algunas propiedades físicas de los gases ideales y se expresan
en función del peso molecular aparente, volumen estándar, densidad, volumen
específico y densidad relativa.
Tabla 4.1 Propiedades físicas para mezclas de gases ideales.
de gas. ..(3.10)
=
Densidad Relación entre 𝛾𝑔 𝜌𝜌𝑎𝑖𝑟𝑒
𝑔𝑎𝑠
. . .(3.11)
relativa la densidad de un 𝑀
𝑎
83
Para entender el comportamiento de una mezcla de gases es necesario estudiar las
leyes que estudian sus propiedades físicas de cada gas y se les conoce como reglas de
mezclado, estas relacionan las propiedades de cada gas.
La suma de las presiones parciales de los gases será igual a la presión total.
𝑃𝑚𝑒𝑧𝑐𝑙𝑎 𝑑𝑒 𝑔𝑎𝑠𝑒𝑠 = 𝑃𝐶1 + 𝑃𝐶2 + 𝑃𝐶3 + ⋯ 𝑃𝐶𝑛, ........................ (4.13)
Dice que el volumen total de una mezcla de gases es igual a la suma de los
volúmenes de los componentes puros de la mezcla siempre y cuando se encuentren a
la misma temperatura y presión.
𝑣𝑗 𝑉𝑗𝑉 𝑉𝑗 , .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (4.17)
esta ecuación es equivalente a la ecuación 4.11, por lo que para un gas ideal, la
fracción volumen y la fracción mol de un componente C son iguales.
𝑦𝑗 = 𝑛 𝑗. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (4.18) 𝑛
Donde:
𝑣𝐶 = Fracción volumen de un componente en la fase gas.
𝑉𝐶 = Volumen ocupado por un componente.
𝑉 = Volumen total de la mezcla de gases.
84
𝑦𝐶= Fracción mol de un componente en la mezcla de gases.
𝑛𝐶 = número de moles del componente C.
𝑛 = Moles totales de la mezcla de gases.
4.3.4 Fracción peso.
𝑤𝑗 𝑚𝑗 𝑚𝑗 . .......................................... (4.19)
Donde
𝑤𝑗 = Fracción peso del componente j en la fase gas.
Donde:
𝑀𝑗 = Masa molecular del componente j.
Por lo que la masa total del sistema queda expresada de la siguiente manera:
𝑛
𝑚 𝑗
𝑚𝐶 , . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (4.21)
85
𝑚𝑗
...............................................
𝑀𝑤𝑗 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (4.24)
𝑦𝑗 𝑗
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (4.26)
4.4 Comportamiento de los gases reales.
Donde :
𝑉𝑟𝑒𝑎𝑙= Volumen del gas real (𝑓𝑡3) 𝑉𝑖𝑑𝑒𝑎𝑙= Volumen del gas ideal (𝑓𝑡3)
Los gases reales difieren su comportamiento de los gases ideales con el incremento
de los valores de presión y temperatura y los cambios en la composición a la que son
sometidos, por lo que 𝑧 es un factor de corrección del gas.
4.4.1 Densidad de un gas real
86
4.4.2 Factor de compresibilidad (𝒛)
Figura 4.124: Correlación de Standing y Katz para la obtención del factor “z” (H. Whitson,
1994)
87
Establece que los gases reales se comportan de manera similar al ser analizados
desde la presión, temperatura y volumen reducidos, es decir expresan cada variable en
función de los valores críticos.
Presión reducida:
𝑃 ................................................... (4.30)
𝑃𝑟 = ,
𝑃𝐶
temperatura reducida:
𝑇 .................................................. (4.31)
𝑇𝑟 = ,
𝑇𝑐
volumen reducido:
𝑉 .................................................... (4.32)
𝑉𝑟 = .
𝑉𝑐
Donde
𝑃𝑟, 𝑇𝑟 y 𝑉𝑟 son la presión, temperatura y volumen reducidos, respectivamente.
𝑃𝑐, 𝑇𝑐 y 𝑉𝑐 son la presión, temperatura y volumen críticos, respectivamente.
𝑗=1 temperatura
pseudocrítica:
𝑛
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (4.34) 𝑇𝑝𝑐 = ∑ 𝑦𝑗𝑇𝑐𝑗, .
𝑗=1
88
Donde:
𝑃𝑝𝑐: Presión pseudocrítica (𝑝𝑠𝑖𝑎).
Presión pseudoreducida:
temperatura pseudoreducida:
donde
𝑃𝑝𝑟: Presión pseudoreducida
más pesados.
89
Figura 4.2: Propiedades pseudocríticas del heptanos y más pesados.
amargos.
90
Método de Wichert-Aziz
Donde:
𝑇𝑝𝑐, : Temperatura pseudocrítica corregida (°𝑅) 𝑦𝐻2𝑆: Fracción mol del
𝐻2𝑆
91
Figura 4.3: Factor de ajuste de la temperatura pseudocrítica para gases amargos (Mc Cain,
1990).
Ejercicio 4.1
92
En un recipiente se tienen 𝟏𝟓 𝒍𝒊𝒕𝒓𝒐𝒔 de una gas a un presión de 𝟏. 𝟓 𝒂𝒕𝒎. ¿Cuál
será el volumen que ocupará este misma masa de gas si manteniéndose la
temperatura constante se lleva a una presión de 𝟏. 𝟖 𝒂𝒕𝒎?
Solución:
Solución:
Solución:
𝑉1 𝑉2
= ,
𝑛1 𝑛2
93
𝑉1𝑛2 (2 𝑙)(0.025 𝑚𝑜𝑙)
𝑛1 = == 𝟎. 𝟏 𝒎𝒐𝒍.
𝑉2 0.5𝑙
Ejercicio 4.4
Una esfera elástica está llena de nitrógeno y tiene un radio de 𝟐𝟎 𝒄𝒎. Si tiene
una temperatura de 𝟐𝟒 °𝑪 y una presión interna de 𝟏. 𝟎𝟐 𝒂𝒕𝒎. Determinar el
número de moles de nitrógeno en la esfera y la masa de nitrógeno necesarias
para inflar la esfera a estos valores.
Solución:
Ejercicio 4.5
Solución:
La masa molecular se obtiene de Tabla 3.1 del capítulo 3, de donde M = 72.015 𝑙𝑏𝑚
.
𝑚𝑜𝑙
𝑝𝑔2𝑎𝑏𝑠 𝑓𝑡
(10.732 ) (610 𝑅)
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
ideales con 𝟓𝟓% de metano, 𝟐𝟓% de etano y 𝟐𝟎% de propano, calcular la presión
parcial ejercida por cada componente dentro de la mezcla.
Solución:
Se observa que al ser una mezcla de gases, cada gas ejercerá una presión parcial y
es la suma de estos la que proporciona la presión total de la mezcla, por lo que se
emplea la ley de las presiones parciales de Dalton de la ecuación 4.15:
𝑃𝐶 = 𝑋𝐶 𝑃𝑚𝑒𝑧𝑐𝑙𝑎
𝑙𝑏 𝒍𝒃
𝑃𝑚𝑒𝑡𝑎𝑛𝑜 = (. 55) (1300 𝑝𝑔 2𝑎𝑏𝑠) = 𝟕𝟏𝟓 𝒑𝒈 𝟐𝒂𝒃𝒔.
95
𝑙𝑏 𝒍𝒃
𝑃𝑒𝑡𝑎𝑛𝑜 = (. 25) (1300 𝑝𝑔 2𝑎𝑏𝑠) = 𝟑𝟐𝟓 𝒑𝒈 𝟐𝒂𝒃𝒔.
𝑙𝑏 𝒍𝒃
𝑃𝑝𝑟𝑜𝑝𝑎𝑛𝑜 = (. 20) (1300 𝑝𝑔 2𝑎𝑏𝑠) = 𝟐𝟔𝟎 𝒑𝒈 𝟐𝒂𝒃𝒔.
𝒍𝒃
𝑃𝑡𝑜𝑡𝑎𝑙 = 𝟏𝟑𝟎𝟎 𝒑𝒈 𝟐𝒂𝒃𝒔.
Ejercicio 4.7
Se tiene una mezcla de gases con 𝟐𝟓% de metano, 𝟏𝟓% de etano, 𝟏𝟕% de
propano, 𝟏𝟑% de n-butano y 𝟑𝟎% de n-pentano, calcular su composición en
fracción mol.
Solución:
Para calcular la fracción mol se emplea la ecuación 4.26, los pesos moleculares se
obtienen de la Tabla 3.1 del capítulo 3, la Tabla 4.2 resume las operaciones y
resultados.
Tabla 4.2 Solución del ejercicio 4.7.
Ejercicio 4.8
Una mezcla de gases está conformada por 𝟑𝟓% de metano, 𝟐𝟓% de etano, 𝟏𝟕%
de propano y 𝟐𝟑% de n-butano. Calcular su densidad relativa.
Solución:
96
La densidad relativa se calcula con 𝛾𝑔 = 28 𝑀 𝑎
.96 que es la ecuación 4.12, para 𝑀𝑎
se emplea la ecuación 4.8, 𝑀𝑎 = ∑𝑛𝑐=1 𝑦𝑗𝑀𝑗, la Tabla 4.3 resume las operaciones
realizadas.
Sustituyendo:
𝑀𝑎 33.9961
𝛾𝑔 = == 𝟏. 𝟏𝟕𝟒.
28.96 28.96
Ejercicio 4.9
Solución:
Para calcular la presión pseudocrítica, se emplea la ecuación 4.33 𝑃𝑝𝑐 = ∑𝑛𝑗=1 𝑦𝑗𝑃𝑐𝑗,
las fracciones mol, se tienen del problema anterior y las propiedades críticas se
obtienen de la Tabla 3.1 del capítulo 3. La Tabla 4.4 muestra el desarrollo de los
cálculos para la solución.
Tabla 4.4 Solución de la presión pseudocrítica del ejercicio 4.9.
97
Componente Fracción Presiones 𝒚𝒋 𝑷𝒄𝒋
mol ( 𝒚𝒋) críticas (psia)
Metano 0.35 660.4 231.14
Etano 0.25 706.5 176.625
Propano 0.17 616 104.72
n-Butano 0.23 550.6 126.638
𝑃𝑝𝑐 = 639.123
Para la te mperatura pseudocrítica se usa la ecuación 4.34 𝑇𝑝𝑐 = ∑𝑛𝑗=1 𝑦𝑗 𝑇𝑐𝑗 , las
Ejercicio 4.10
Para una mezcla de gases con 𝟓𝟐% de metano, 𝟏𝟖% de etano, 𝟏𝟏% de propano
y 𝟏𝟗% de n-butano, calcular su factor de compresibilidad 𝒛, si está contenida en
𝟐𝟓, 𝟓𝟎𝟎 𝒇𝒕𝟑, con una presión de 𝟒 𝟓𝟎𝟎 𝒑𝒔𝒊𝒂 y a 𝟑𝟏𝟎°𝑭. Utilizar la Fig. 4.1.
Solución:
A partir de las propiedades pseudoreducidas se puede utilizar la Fig. 4.1, por lo que
se calculan las propiedades pseudocríticas para entonces sustituir en las ecuaciones
4.35 y 4.36. La Tablas 4.6 y 4.7 muestran los cálculos necesarios para obtener las
propiedades pseudocríticas.
Tabla 4.6 Solución de la presión pseudocrítica del ejercicio 4.10.
Componente Fracción Presiones 𝒚𝒋𝑷𝒄𝒋
mol (𝒚𝒋) críticas (psia)
Metano 0.52 660.4 343.408
Etano 0.18 706.5 127.17
Propano 0.11 616 67.76
n -Butano 0.19 550.6 104.614
98
𝑃𝑝𝑐 = 642.952
99
Figura 4.4: Solución del ejercicio 4.10.
Ejercicio 4.11
Se tiene una mezcla de gases en la que se incluyen heptanos y más pesados
(densidad relativa = 𝟎. 𝟗𝟓, M = 𝟏𝟕𝟎 𝒍𝒃𝒎/𝒍𝒃𝒎 − 𝒎𝒐𝒍), la composición de esta se
muestra en la Tabla 3.8 a 𝟏𝟐𝟎𝟎 𝒑𝒔𝒊𝒂 y 𝟐𝟎𝟎°𝑭. Obtener la presión y temperatura
pseudoreducidas y su factor de compresibilidad 𝒁, usar la Fig. 4.1.
100
Se calculan las propiedades pseudocríticas de cada componente, solo para el caso
de los heptanos y más pesados, estás se obtienen de la Fig. 4.2, la solución se muestra
en la Fig. 4.5.
101
𝑇 659.69°𝑅
𝑇𝑝𝑟 = == 1.10.
𝑇𝑐 598.98
Con estos valores se obtiene el factor de compresibilidad 𝑧 de la Fig. 4.1, donde 𝒛 =
𝟎. 𝟑𝟕.
Ejercicio 4.12
Se tiene una mezcla de gases cuya composición se muestra en la Tabla 4.10,
calcular su densidad relativa.
𝑃𝑗
𝑃𝑗 = 𝑦𝑗 𝑃𝑚𝑒𝑧𝑐𝑙𝑎 𝑦𝑗 = 𝑃 𝑚𝑒𝑧𝑐𝑙𝑎
102
Propano 1.2 0.04 44.1 2.03
i-Butano 0.9 0.03 58.123 2.00
n-Butano 0.7 0.02 58.123 1.56
i-Pentano 0.3 0.01 72.15 0.83
n-Pentano 0.2 0.007 72.15 0.55
26.04 𝑀𝑎 = 22.88
𝑃𝑚𝑒𝑧𝑐𝑙𝑎 =
Ejercicio 4.13
Para una mezcla de gases cuya composición se muestra en la Tabla 4.12,
calcular sus valores de presión y temperatura pseudoreducidas con el método de
Wichert-Aziz. Considerar para los heptanos y más pesados 𝜸𝒈 = 𝟎. 𝟖𝟎𝟕 𝒚 𝑴 =
𝒍𝒃𝒎
𝟏𝟒𝟐 .
𝒍𝒃𝒎−𝒎𝒐𝒍
𝐵 = 𝑦𝐻2𝑆 = 0.011,
donde
𝑉0 y 𝑃0, son condiciones iniciales de volumen y presión.
𝐶 es el promedio de los coeficientes de compresión desde 𝑝 hasta 𝑝0.
Es la presión ejercida por el vapor sobre un líquido cuando está en equilibrio, la Fig.
5.1 representa el comportamiento de la presión de vapor con respecto a la temperatura.
105
Figura 25: Comportamiento de la presión de vapor de un líquido con respecto a la
temperatura.
106
Figura 5.2 Logaritmo de la presión de vapor contra el reciproco de la temperatura absoluta
para n-Hexano (Burcik, 1979).
donde
𝑑𝑃°
: razón de cambio de la presión de vapor con respecto a la temperatura,
𝑑𝑇
Δ𝐻 : calor de vaporización para una cantidad de líquido,
107
𝑇: temperatura absoluta y
Δ𝑉: cambio de volumen del líquido desde la fase liquida hasta la fase gaseosa(𝑉𝑔
− 𝑉𝑙).
Si se considera que el vapor es un gas perfecto, la ecuación Clausius-Clapeyron se
simplifica como:
𝑑𝑃
= Δ𝐻𝑚,
°
................................................ (5.3)
𝑑𝑇 T𝑉𝑔 todas las variables deben estar calculadas para una mol de líquido,
donde:
Como el volumen molar del líquido es mucho muy pequeño con respecto al volumen
del gas, se considera despreciable.
Para el cálculo de 𝑉𝑔 se tiene:
O bien
𝑑 ln 𝑃°
𝑚 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Δ𝐻 𝑃°
(5.6)
=
𝑑 ln 𝑇 R𝑇2
donde
𝐶: Constante de integración.
108
Esta ecuación expresa claramente que el logaritmo de la presión de vapor es una
función lineal del inverso de la temperatura 1, por lo que la pendiente será − ∆𝐻𝑚. Si
𝑇 𝑅
1 ). ....................................... (5.8)
°
𝑙𝑛 𝑃1 ° = − ∆𝐻𝑚 ( 1 −
𝑃2 𝑅 𝑇1 𝑇2
También conocido como calor de vaporización, representa la energía que debe ser
aplicada para evaporar una mol de un líquido y no depende de la temperatura a la cual
se lleva a cabo la vaporización. Sin embargo este valor no varía mucho en un rango
pequeño de temperaturas, por lo que se puede considerar como constante.
109
Etano, C2 H6 6344 -128.2 19.1
Promedio 20
Se observa que los cocientes de los calores latentes entre la temperaturas del punto
de ebullición son muy similares, si se realiza un promedio de estos valores se obtiene
19.7 que redondeado es 20, a esto se le conoce como Regla de Trouton, funciona para
varios tipos de líquidos, en la literatura el valor de esta regla es de 21, su aplicación
principal es para calcular el calor latente de un líquido cuando se conoce su punto de
ebullición según Burcik (1979).
5.2 Comportamiento de líquidos ideales.
Los líquidos no soportan esfuerzos cortantes cuando están en equilibrio, sin embargo
los soportan entre capas adyacentes en un movimiento relativo, debido a que tienen
viscosidad. Un fluido ideal es aquel que no presenta esfuerzos cortantes, incluso
cuando están en movimiento, por lo que las fuerzas superficiales que sobre este ejerzan
son debidas a la presión, además de que no tiene viscosidad y por lo tanto no pierden
energía mecánica. Si un líquido ideal fluye estacionariamente (a una velocidad
constante) a través de una tubería horizontal estrecha y de sección transversal
constante, la presión será constante a lo largo de la tubería.
110
Figura 5.3: Perfil de velocidades de un fluido ideal y un fluido real.
donde:
𝑄: Caudal(𝑓𝑡 3
).
𝑠
paredes.
fluido.
111
Figura 5.4: Caída de presión que experimenta un líquido en movimiento.
5.3.1 Coeficiente de viscosidad.
donde:
𝐹: Fuerza de desplazamiento de la capa superior 𝑍: distancia entre capas
𝑣: velocidad de la capa superior 𝐴: Área de la capa
112
por lo que la ley de Poiseuille queda definida como:
113
Figura 5.6: Fluido en movimiento a través de una sección de tubería con diferente
diámetro.
Dado que el volumen 𝑉 considerado, ingresa con una velocidad 𝑣, el sistema recibe
una energía cinética y en virtud de su posición también recibe una energía potencial.
Por lo que la energía mecánica total del sistema es:
Dado que en los extremos izquierdo (1) y derecho (2) no hay acumulación ni salida
de fluido y que el fluido es incompresible, la energía mecánica que ingresa por la parte
izquierda es igual a la que egresa por la parte derecha.
114
5.4 EJERCICIOS RESUELTOS
Ejercicio 5.1.
El calor latente de un hidrocarburo es de 𝟓𝟑𝟔𝟎 𝒄𝒂𝒍, su presión de vapor a 𝟕° 𝑪
es de 𝟏𝟎𝟎𝟎 𝒎𝒎 𝑯𝒈, ¿cuál será su presión de vapor en psia a 𝟐𝟎°𝑪?.
Solución:
𝑐𝑎𝑙
∆𝐻𝑚 = 5,360 𝑐𝑎𝑙 𝑅 = 1.987 (𝑚𝑜𝑙 °𝐾 )
115
𝑃1° 5,360 𝑐𝑎𝑙 1 1
𝑙𝑛 = (− ) = 0.427,
1000 1.987 ( 𝑐𝑎𝑙 ) 293°𝐾 280°𝐾
𝑚𝑜𝑙°𝐾
𝑃1° °
= 1. = 1530 𝑚𝑚 𝐻𝑔.
53 → 𝑃1
1000
En 𝑝𝑠𝑖𝑎:
°
=
𝟐𝟗. 𝟔 𝒑𝒈𝒍𝒃𝟐𝒂𝒃𝒔. 𝑃1
Ejercicio 5.2
Calcular la presión de vapor en 𝒑𝒔𝒊𝒂 para el n-Hexano a 𝟑𝟎°𝑪.
Solución:
De la Fig. 5.2, se observa que la temperatura debe estar en °𝐾, por lo que 30°𝐶
equivale a 300°𝐾, ahora se obtiene el inverso de este para poder ubicarlo sobre la
gráfica y esto es igual a 0.003. La Fig. 5.7 muestra la solución grafica del ejercicio.
116
Figura 5.7 Solución gráfica del ejercicio 5.2.
Solución:
𝑃° = 6.36 𝑝𝑔 𝑙𝑏
2𝑎𝑏𝑠 𝑇1 = 75°𝐹 = 534.67°𝑅 𝑇2 = 90°𝐶 = 549.67°𝑅
𝑙𝑏 3
𝐵𝑇𝑈 𝑝𝑔2𝑎𝑏𝑠𝑓𝑡
𝑃2° ∆𝐻𝑚 1 1
𝑙𝑛 𝑃° = 𝑅 (𝑇1 − 𝑇2)
1
𝐵𝑇𝑈
𝑃° 8450 ( )
𝑙𝑛 1 = 𝑙𝑏 𝑚𝑜𝑙 ( 1 − 1 ) = 0.0401
𝑙𝑏 𝑙𝑏 3 534.67°𝑅 549.67°𝑅
6.36 𝑝𝑔2𝑎𝑏𝑠 𝑝𝑔2𝑎𝑏𝑠 𝑓𝑡
10.732 ( )
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑃1°
= 1.0409
𝑙𝑏
6.36
𝑝𝑔2𝑎𝑏𝑠
°
= 1.0409 (6.36 𝑝𝑔𝑙𝑏2𝑎𝑏𝑠 ) = 𝟔. 𝟔𝟐𝟎𝟏
𝒑𝒈𝒍𝒃𝟐𝒂𝒃𝒔 .
𝑃1
117
Ejercicio 5.4
La presión de vapor del etanol es de 𝟏𝟑𝟓. 𝟑 𝒎𝒎𝑯𝒈 a 𝟒𝟎°𝑪 y 𝟓𝟒𝟐. 𝟓 𝒎𝒎𝑯𝒈 a
𝟕𝟎°𝑪.
Calcular el calor molar de vaporización y la presión del etanol a 𝟓𝟎°𝑪.
Solución:
𝑃1° ∆𝐻𝑚 1 1
𝑙𝑛 𝑃° = 𝑅 (𝑇1 − 𝑇2)
2
𝑅 𝑙𝑛 𝑃1°°
∆𝐻𝑚 = 1 𝑃12
( − )
𝑇1 𝑇2
581 49
Considerando 𝑃°1 = 135,3 𝑚𝑚𝐻𝑔 = 2. 𝑝𝑔
𝑙𝑏
2𝑎𝑏𝑠 𝑃°2 = 542.5 𝑚𝑚𝐻𝑔 = 10. 𝑝𝑔
𝑙𝑏
2𝑎𝑏𝑠
Conociendo ∆𝐻𝑚, la presión de vapor 𝑃°3 a 50°𝐶 (𝑇3) se obtiene utilizando como datos
𝑃°1 y 𝑇1 así como despejando y sustituyendo en la ecuación 5.8.
𝑙𝑏1
𝑃3° = 2.581 𝑝𝑔 2𝑎𝑏𝑠 𝑒− (−563.67) = 𝟒. 𝟐𝟒 𝒑𝒈 𝒍𝒃𝟐𝒂𝒃𝒔
Ejercicio 5.5
118
Por una tubería de 𝟏/𝟖 de pulgada (𝟎. 𝟑𝟏𝟕𝟓𝒄𝒎) de diámetro pasa aceite de
𝟐𝟎°𝑪 y densidad de 𝟎. 𝟖 𝒄𝒎 𝒈𝒓
, descargando a la atmósfera con un gasto de 𝟎. 𝟏
𝟑
𝒎𝒍
𝒔 . Para medir la caída de presión en la tubería se colocan dos tubos
h
30 cm
Figura 5.8: Del ejemplo 1, representa la distancia entre dos tubos manométricos y la
diferencia de alturas debido a la caída de presión de un fluido laminar viscoso.
Solución:
𝑄 0.1𝑥10−6 𝑚3
𝑠
𝑣 = 𝐴 = 7.92𝑥10−6𝑚 2 = 1.26𝑥10−2 𝑚 𝑠 = 1.26 𝑐𝑚 𝑠
Ejercicio 5.6
119
Para el ejercicio 5.5, calcule la caída de presión en 𝒇𝒕 de altura
equivalentes entre los dos tubos manométricos.
Solución:
La caída de presión entre los dos puntos de la tubería está dada por:
∆𝑃 = 8 𝜋𝜂𝐿𝑟4 𝑄 = 8(30𝑥𝜋10(0−.00158753 𝑁.
𝑠⁄𝑚2𝑚)𝑥)40.30𝑚 (0.1𝑥10−6 𝑚3⁄𝑠) = 360𝑃𝑎.
La diferencia de altura debida entre los dos tubos manométricos se obtiene con:
∆𝑃 360𝑃𝑎
ℎ= = = 0.045 𝑚 = 𝟎. 𝟎𝟏𝟒𝟕𝒇𝒕.
𝜌𝑔 800 𝑘𝑔 𝑚
( 𝑚3 )(9.8 𝑠2)
𝒈𝒓
𝟎. 𝟖 𝒄𝒎 𝟑, 𝒎𝟏 = 𝟒𝟎𝟎 𝒈𝒓, considerando despreciable la masa de aire en las
burbujas, calcule:
b) La densidad de la mezcla.
Solución:
𝑉1 = 𝑚𝜌11 = 6001 𝑐𝑚𝑔𝑟 𝑔𝑟3 = 600 𝑐𝑚3, 𝑉2 = 𝑚𝜌22 = 0400.8 𝑐𝑚 𝑔𝑟𝑔𝑟3 = 500 𝑐𝑚3,
120
despejando 𝑉𝑏 = 𝑉𝑀 − 𝑉1 − 𝑉2 entonces:
𝑉𝑏 = 1200𝑐𝑚3 − 600 𝑐𝑚3 − 500 𝑐𝑚3 = 𝟏𝟎𝟎 𝒄𝒎𝟑.
b) Para la densidad de la mezcla:
𝑚𝑚 𝑚1 + 𝑚2 + 𝑚𝑏 600 𝑔𝑟 + 400𝑔𝑟 + 0 𝑔𝑟 𝒈𝒓
𝜌𝑀 = 𝑉 = 𝑉 = 1200𝑐𝑚3
= 𝟎. 𝟖𝟑 𝒄𝒎 𝟑.
𝑚 𝑚
Capítulo 6.
Ecuaciones de estado cúbicas.
121
Las ecuaciones de estado cúbicas tienen como principal objetivo corregir la ley de los
gases ideales para que pueda ser aplicable para los gases reales, entre las principales
se encuentran la que se describen a continuación.
6.1 Ecuación de estado de van der Waals. (Ecuación de estado de dos constantes
ó ecuación cúbica).
La ecuación de estado de van der Waals corrige estas suposiciones para explicar el
comportamiento de los gases reales, explica teóricamente que las moléculas de gas
ocupan un volumen significativo al formar parte de un volumen contenido en un
recipiente y que existen fuerzas de atracción entre las moléculas. Limita su aplicación a
presiones bajas.
La ecuación 6.1 corresponde a la ecuación de estado de van der Waals:
𝑎 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.1) (𝑝 + 𝑉 𝑀2 ) (𝑉𝑀 − 𝑏) =
𝑅𝑇,
donde:
𝑎,𝑏 : Constantes, dependen del tipo de gas.
𝑙𝑏
3
). 𝑅: Constante de los gases 10.732( 𝑝𝑔
2𝑎𝑏𝑠𝑓𝑡3 .
𝑉𝑀: Volumen molar real (𝑓𝑡 𝑙𝑏𝑚−𝑚𝑜𝑙−°𝑅)
𝑅𝑇 .......................................... (6.3)
𝑃𝑟𝑒𝑝𝑢𝑙𝑠𝑖ó𝑛 = ,
𝑉𝑀 − 𝑏
122
𝑃𝑎𝑡𝑟𝑎𝑐𝑐𝑖ó𝑛 = 𝑉 𝑎𝑀2. ............................................. (6.4)
La ecuación 6.5 representa la ecuación de van der Waals considerando los puntos
críticos de la sustancia:
𝑎 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.5) (𝑝𝑐 + 𝑉 𝑀𝑐2 ) (𝑉𝑀𝑐 − 𝑏)
= 𝑅𝑇𝐶, .
donde:
27 𝑅2𝑇𝑐2 ................................................ (6.6)
𝑎= ,
64𝑝𝑐
de igual manera las ecuaciones 6.9 y 6.10, permiten obtener los valores de las
constantes a y b:
, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.9) 𝑏 = Ω𝑏
2 2
𝑎 = Ω𝑎 𝑅 𝑃𝑇𝑐𝑐
, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.10)
2
𝑅𝑃 𝑇𝑐𝑐
donde:
Ω𝑎 = 0.421875 Ω𝑏 = 0.125 .
123
3
− (𝑏 + 𝑅𝑇) 𝑉𝑀2 + 𝑎 𝑉
𝑀 − (𝑎𝑏) = 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.11)
𝑉𝑀
𝑃 𝑝 𝑝
donde:
𝑏𝑝 ................................................... (6.14)
𝐵= .
𝑅𝑇
de vapor), en este caso la raíz más grande corresponde a la fase vapor y la raíz
124
𝑎 (6.15)
[𝑝 + 𝑇 ] (𝑉𝑀 − 𝑏) = 𝑅𝑇, . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1/2𝑉𝑀(𝑉𝑀+𝑏)
𝐴=
donde:
𝑅2𝑇2.5
𝑅2𝑇𝑐2.5 y . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.16)
𝑏𝑝 ................................................. (6.20)
𝑎=𝐵=Ω𝑎 .
𝑅𝑇
𝑃𝑐
𝑅 𝑇𝑐 ............................................... (6.17)
𝑏 = Ω𝑏 .
𝑃𝑐
de donde
𝑎𝑝 y............................................... (6.19)
𝑛 2 .......................................... (6.21)
𝑎𝑚 ,
𝑗=1
𝑛 .......................................... (6.22)
𝑏𝑚 = [∑ 𝑥𝑗𝑏𝑗] .
𝑗=1
Para mezclas de hidrocarburos en fase gas, están representadas por las ecuaciones
6.23 y 6.24:
𝑛 ......................................... (6.24)
𝑏𝑚 = [∑ 𝑦𝑗𝑏𝑗] .
𝑗=1 donde:
125
𝑥𝑗: Fracción mol del componente 𝑗 de la fase líquida.
𝐴= 𝑎𝑚𝑝. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.25)
𝑅2𝑇2.5
𝑎𝑚 𝑝
................................................. (6.26)
𝐵=
𝑅𝑇
[𝑝 + 𝑉 ] (𝑉𝑀 − 𝑏) = 𝑅𝑇, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
𝑀(𝑉𝑎𝑀𝑇 +𝑏) (6.27)
𝑎𝑇 = 𝑎𝑐𝛼. ................................................ (6.28)
(6.30)
𝑚 = 0.480 + 1.574𝜔 − 1.76𝜔2, .................................
.
𝜔 = −(log 𝑝𝑣𝑟 + 1). .......................................... (6.31)
donde :
𝑝𝑣𝑟: Presión de vapor reducida a una 𝑇𝑟 = 0.7 Para
las constantes 𝑎𝑐 𝑦 𝑏.
126
𝑅2𝑇𝑐2 .............................................. (6.32)
𝑎𝑐 = Ω𝑎 ,
𝑃𝑐
𝑅𝑇𝐶 ............................................... (6.33)
𝑏 = Ω𝑏 .
𝑃𝑐
Donde:
𝑏𝑝 ................................................... (6.35)
𝐵= .
𝑅𝑇
6.4 Ecuación de estado de Peng-Robinson.
Para m:
Sí ω ≤ 0.49 𝑚 = 0 .7464 + 1.54226𝜔 − 0.26992𝜔2 . . . . . . . . . . . (6.37)
Peng-Robison.
Están dadas por las ecuaciones 6.42 y 6.43 para la fase líquida:
128
𝑎𝑇 = ∑ ∑ [𝑥𝑖𝑥𝑗(𝑎𝑇𝑖𝑎𝑇𝑗) (1 − 𝛿𝑖𝑗)], ........................... 1/2 (6.42)
(6.43)
𝑏 = ∑(𝑥
𝑇𝑖𝑗 𝑗𝑏𝑗). .............................................. 𝑖𝑗
𝑇𝑖 𝑇𝑗
donde:
(6.46)
𝑎 = (1 − 𝛿 )(𝑎 𝑎 ).
.................................. 1/2
c) Se forma una matriz simétrica con los coeficientes de interacción binaria 𝛿𝑖𝑗 = 𝛿𝑗𝑖.
6.5 Otras ecuaciones de estado cúbicas.
Esta ecuación es una modificación a la ecuación de Van der Waals y toma en cuenta
la dependencia entre las fuerzas de atracción y la temperatura, es más exacta a
presiones y temperaturas bajas:
Dieterici razonó que para acercar una molécula a las paredes del recipiente es
necesario vencer la presión de las moléculas que están cerca de esta, explica que la
presión cerca de esta área debe ser menor debido a que debe haber menos moléculas
129
que en la parte central del recipiente y establece la relación entre las presiones con la
ecuación 6.48.
(6.48)
𝑝 = 𝑝𝑖𝑒(−𝑅𝑇), .............................................. 𝐴
donde
𝐴= 𝑎
. .. ............................................. (6.49)
𝑉𝑚
Es útil para cálculos con presiones altas, usa cinco constantes además de la
constante 𝑅 y está representada por la ecuación 6.52.
𝑐 ............................... (6.52)
𝑅𝑇 [1 − ( )]
3 𝐴
𝑃= 𝑉 𝑉 2 𝑉𝑇
donde:
(𝑉 + 𝐵) − 2.
𝑎 ............................................ (6.53)
𝐴 = 𝐴0 (1
− ),
𝑉
𝑏 ............................................ (6.54)
𝐵 = 𝐵0 (1 − ).
𝑣
130
Los valores de las constantes 𝐴0, 𝐵0, 𝑎, 𝑏 y c, tienen un valor en particular
dependiendo del gas en cuestión.
𝑅𝑇 𝐴 .. ............................... (6.55)
𝑉𝑀 = ( + 𝐵) (1 − 𝐸) − .
𝑝 𝑅𝑇
Donde
0
𝑅𝑇
𝑏𝑝 . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.57)
𝐵 = 𝐵0 (1 − ),
𝑅𝑇
𝐴 = 𝐴 (1 − 𝑎𝑝 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.56) ),
𝐸 = 𝑐𝑝
......... ......................................... (6.58)
.
131
𝑅𝑇4
Los valores de las constantes son los mismos de la ecuación original, estas fueron
obtenidas de manera experimental para gases puros, para mezclas de gases se usarán
las siguientes expresiones:
𝑗=𝑛 ............................................. (6.59)
𝑎 = ∑ 𝑦𝑗𝑎𝑗,
𝑗=1
𝑗=𝑛 . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.60)
𝑏 = ∑ 𝑦𝑗𝑏𝑗,
𝑗=1
𝑗=𝑛 ................................................ (6.61)
𝑐 = ∑ 𝑦𝑗𝑐𝑗.
𝑗=1
Para 𝐴0 y 𝐵0
𝑗=𝑛
2
...................................... (6.62)
𝐴0 = (∑ 𝑦𝑗𝐴0𝑗1/2) ,
𝑗=1
𝑗=𝑛
.......................................... (6.63)
𝐵0 = ∑ 𝑦𝑗𝐵0𝑗.
𝑗=1
132
6.6 EJERCICIOS RESUELTOS
Ejercicio 6.1
Para el n-Pentano calcular el valor de las constantes 𝒂 y 𝒃 de van der Waals.
Solución:
Para obtener 𝑎 y 𝑏, se usan las ecuaciones 6.6 y 6.7, los valores de las propiedades
críticas para el n-pentano se obtienen de la Tabla 3.1.
2
𝑙𝑏 3
2 2
𝑎𝑏𝑠𝑓𝑡
𝑝𝑔
𝑙𝑏 3
2
𝑎𝑏𝑠𝑓𝑡
𝑝𝑔
(10.732 )(845.49°𝑅)
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑙𝑏 𝟑
𝑏 = 8𝑅𝑝𝑇𝑐 = = 𝟐. 𝟑𝟐 𝒍𝒃𝒎 𝒇𝒕− 𝒎𝒐𝒍.
𝑐 8(488.6 𝑝𝑔2𝑎𝑏𝑠)
Ejercicio 6.2
Calcular la densidad a la que se encuentra una muestra de metano puro, cuya
temperatura es de 𝟏𝟐𝟎°𝑭 y la presión es de 5500 psia.
Solución:
Se calculan los valores a y b, con las ecuaciones 6.9 y 6.10, las propiedades críticas
se obtienen de la Tabla 3.1.
2
𝑙𝑏 3
𝑝𝑔
(10.732 2𝑎𝑏𝑠𝑓𝑡 2
) (576.36 °𝑅)
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
133
𝑎 = Ω𝑎 𝑅 2 𝑇 𝑐2
𝑃 𝑐 = (0.421875) 666.40𝑝𝑔2𝑙𝑏𝑎𝑏𝑠 = 8594.7576
2 2
𝑙𝑏𝑎𝑏𝑠𝑓𝑡3
𝑝𝑔
(10.732 )(576.36 °𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑏 = Ω𝑏 𝑅 𝑇 𝑐2
𝑃 𝑐 = (0.125) 666.40𝑝𝑔2𝑙𝑏𝑎𝑏𝑠 = 0.6911422
A partir de los valores obtenidos, se calculan los valores de 𝐴 y 𝐵 con las ecuaciones
6.13 y 6.14.
𝑙𝑏
𝑎𝑝 (8594.7576)(5500𝑝𝑔2 𝑎𝑏𝑠)
𝑅2 𝑇 2
= 2
𝐴 𝑙𝑏
2 𝑓𝑡 3 == 1.2205 y
( 10.732 ) (579.69 °𝑅)2 𝑝𝑔 𝑎𝑏𝑠
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑏𝑝 (0.691142)(5500𝑝𝑔2𝑙𝑏𝑎𝑏𝑠)
𝐵 = 𝑅𝑇 = 𝑝𝑔2𝑙𝑏𝑎𝑏𝑠 3 = 0.61069 .
𝑓𝑡
(10.732 )(579.69 °𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑄 = 3𝑎2−𝑎12 = 0.118426,
9
𝐽 = 9𝑎1𝑎2−27𝑎3−2𝑎13 = 0.19978.
54
𝑍1 = 1.1157 . Y
calculando la densidad:
𝑝𝑀
𝜌= ,
𝑧𝑅𝑇
𝑙𝑏𝑚 𝑙𝑏
(5500 𝑝𝑔 2𝑎𝑏𝑠 )(16.043 𝑚𝑜𝑙 )
𝒍𝒃𝒎 𝑙𝑏
2 𝑓𝑡 3
𝜌=
= 𝟏𝟐. 𝟕𝟏𝟐𝟏 𝒇𝒕𝟑 .
𝑝𝑔 𝑎𝑏𝑠
1.1157(10.732 )(579.69 °𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
Ejercicio 6.3
Con la ecuación de Redlich-Kwong, calcular la densidad de una muestra de
Solución:
2
𝑙𝑏
𝑓𝑡 3
𝑅2𝑇𝑐2.5 𝑝𝑔2 𝑎𝑏𝑠
( 10.732 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅) (549.61°𝑅)2.5
𝑎 = Ω𝑎 = (0.42747)
= 494016.786,
𝑙𝑏
𝑃𝑐 (706.5 )
𝑝𝑔2 𝑎𝑏𝑠
𝑙𝑏
𝑓𝑡 3
𝑝𝑔2 𝑎𝑏𝑠
( 10.732 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅) (549.61°𝑅)
𝑅 𝑇𝑐 𝑙𝑏
𝑏 = Ω𝑏 = (0.08664) (706.5 )
𝑝𝑔2 𝑎𝑏𝑠
= 0.7233.
𝑃𝑐
135
𝑙𝑏
(493504.78 ) (200 )
𝑎𝑝 𝑝𝑔2 𝑎𝑏𝑠
𝐴 = 2 2.5 = 2 = 0.0686,
𝑅 𝑇 𝑙𝑏
𝑓𝑡 3
𝑝𝑔2 𝑎𝑏𝑠
( 10.732 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅) (689.69°𝑅)2.5
𝑏𝑝 (0.7233)(200𝑝𝑔2𝑙𝑏𝑎𝑏𝑠)
𝐵 = 𝑅𝑇 = 2𝑙𝑏𝑎𝑏𝑠 3 = 0.0195.
𝑓𝑡 𝑝𝑔
(10.732 )(689.69°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝜃 ,
sustituyendo este valor en la ecuaciones A.6, A.7 y A.8:
𝑍 ,
𝑍 ,
𝑍 .
Ahora se considera 𝑧1 para la fase gaseosa y 𝑧3 para la fase líquida.
Con 𝑧1 para la fase gaseosa:
𝑙𝑏 𝑙𝑏𝑚
𝑝𝑀 (3200𝑝𝑔2𝑎𝑏𝑠)(40.784𝑚𝑜𝑙) 𝒍𝒃𝒎
Ejercicio 6.4
Se tiene una mezcla de hidrocarburos líquidos, cuya composición se muestra en
𝑻𝒄 = 𝟒𝟔𝟓°𝑹.
𝑎𝑗 = Ω𝑎 𝑃𝑐 𝑏𝑗 = Ω𝑏 𝑃𝑐𝑐
137
Componente Fracción Presión Temperatura Temperatura mol crítica crítica (°F)
crítica (°R)
Metano 0.36 666.4 116.67 586.26
Etano 0.21 706.5 89.92 559.51
Propano 0.16 616 206.06 675.65
i-Butano 0.09 527.9 274.46 744.05
n-Butano 0.0581 550.6 305.62 775.21
i-Pentano 0.0029 490.4 389.1 858.69
n-Pentano 0.008 488.6 385.8 855.39
Hexano 0.0045 436.9 453.6 923.19
Heptano 0.0835 385 1350 1819.59
Total 1
138
𝑙𝑏
𝑏𝑝 (1.1998)(3200 𝑝𝑔2 𝑎𝑏𝑠)
𝐵 = 𝑅𝑇 = 2𝑙𝑏𝑎𝑏𝑠𝑓𝑡3 = 0.55515.
𝑝𝑔
(10.732 )(669.79°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝜃 .
Sustituyendo este valor en las ecuaciones A.6, A.7 y A.8:
𝑍 .
De estas raíces se considera 𝑍1 para calcular la densidad de la fase gaseosa por ser
el factor de compresibilidad más grande y se considera 𝑍3 para la densidad de la fase
líquida, la masa aparente se determina mediante la ecuación 4.8:
𝑛
𝑀𝑎 = ∑ 𝑦𝑗𝑀𝑗 .
𝑗=1
M 𝑦𝑗 𝑀
𝑙𝑏𝑚
Componente
𝑦𝑗
𝑙𝑏𝑚 𝑚𝑜𝑙
Metano 0.36 16.04 5.774
139
Etano 0.21 30.07 6.315
Propano 0.16 44.1 7.056
i-Butano 0.09 58.123 5.231
n-Butano 0.0581 58.123 3.377
i-Pentano 0.0029 72.15 0.209
n-Pentano 0.008 72.15 0.577
Hexano 0.0045 86.177 0.388
Heptano 0.0835 142 11.857
𝑀𝑎 = 40.784
Para la fase gaseosa:
𝑙𝑏 𝑙𝑏𝑚
Ejercicio 6.5
Un recipiente totalmente cerrado contiene butano en fase líquida y gaseosa a
𝟏𝟓𝟎°𝑭. Calcular la densidad de cada fase utilizando la ecuación cúbica de estado
de Soave-Redlich-Kwong.
Solución:
Es necesario conocer los valores críticos de presión y temperatura, así como el factor
acéntrico para el propano, estos de la Tabla 3.1 del capítulo 3, de donde se obtienen los
siguientes valores:
𝑃𝑐 = 616, 𝑇𝑐 = 206.06°𝐹 y 𝜔 = 0.1522 .
Calculando la temperatura reducida con la ecuación 3.7:
𝑇 610°𝑅
𝑇= == 0.9158 .
140
𝑇𝐶 666.06°𝑅
Se determina 𝑚, con la ecuación 6.30:
𝑚 = 0.480 + 1.574𝜔 − 1.76 = 0.480 + 1.574(0.1522) − 1.76(0.1522) 2 = 0.6787 , se
calcula el valor de 𝛼 con la ecuación 6.29:
𝛼 = [1 + 𝑚(1 − 𝑇𝑟1/2)]2 = [1 + 0.6787(1 − (0.9158) 1/2)]2 = 1.0592 ,
6.32 y 6.33:
2
𝑙𝑏
𝑓𝑡 3
𝑝𝑔2 𝑎𝑏𝑠
( 10.732 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅 ) (666.06°𝑅)2
𝑅2
𝑇𝑐2
𝑙𝑏
𝑎𝑐 = Ω𝑎 (616 ) =
𝑝𝑔2 𝑎𝑏𝑠
(0.42748)
𝑙𝑏
𝑃𝑐 = 35 458.72 .
𝑝𝑔2 𝑎𝑏𝑠𝑓𝑡 6 (𝑙𝑏𝑚 𝑚𝑜𝑙)2
𝑙𝑏
𝑓𝑡 3
𝑝𝑔2 𝑎𝑏𝑠
( 10.732 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅) (666.06°𝑅)
𝑙𝑏
(616 2 )
𝑅𝑇𝐶𝑓𝑡3 𝑏 = Ω𝑏
= (0.08664) =
1.0053 .
𝑃𝑐 𝑙𝑏𝑚 𝑚𝑜𝑙
𝑝𝑔 𝑎𝑏𝑠
Con el valor obtenido de 𝑎𝑐, se procede a calcular 𝑎𝑇 , esto con la ecuación 6.28:
𝑙𝑏
𝑎𝑇 = 𝑎𝑐𝛼 = (35 458.72 𝑝𝑔 2 𝑎𝑏𝑠𝑓𝑡6(𝑙𝑏𝑚 𝑚𝑜𝑙)2) (1.0592)
𝑙𝑏
= 37 557.87 𝑝𝑔 2𝑎𝑏𝑠𝑓𝑡6(𝑙𝑏𝑚 𝑚𝑜𝑙)2 .
141
Y con este valor, ahora se calculan las constantes 𝐴 y 𝐵 y esto con las ecuaciones
6.34 y 6.35. La presión se obtiene de las cartas de Cox de la Fig. 3.5 de donde se lee
( )
𝑙𝑏
𝑏𝑝 (1.0053) (480 𝑝𝑔2𝑎𝑏𝑠)
𝐵= = = 0.0737 ,
𝑅𝑇 𝑙𝑏2𝑎𝑏𝑠 𝑓𝑡3
𝑝𝑔
(10.732 ) (610°𝑅)
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑝𝑔
(0.8089)(10.732 )(610°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
142
Ejercicio 6.6
Se tiene una mezcla de hidrocarburos cuya composición se muestra en la tabla
5.6, existe en equilibrio en dos fases a 𝟐𝟖𝟎𝟎 𝒍𝒃/𝒑𝒈𝟐 𝒂𝒃𝒔 y 𝟏𝟖𝟎°𝑭. Calcular las
densidades de las fases vapor y líquido empleando la ecuación de Soave-
RedlichKwong, SRK. Considerar que no existe interacción binaria entre los
componentes de la mezcla.
Tabla 6.6 Composición de la mezcla del ejercicio 6.6.
xj yj
Componente
(fracción mol líquido) (fracción mol vapor)
Solución:
𝑅2𝑇𝑐2 𝑅𝑇
𝑎 = Ω𝑎 𝑃𝑐 𝑏 = Ω𝑏
𝐶
𝑃𝑐 Ω𝑎 =0.42748 Ω𝑏=0.08664
El término 𝑎𝑇𝑖𝑗 del metano líquido recordando que no hay interacción binaria entre los
Los 𝑎𝑇𝑖𝑗 de las fracciones restantes tanto de la fase líquida como la fase vapor se
obtienen con el mismo procedimiento. La Tabla 6.8, resume los resultados de estos, así
como los resultados del parámetro 𝑏.
145
C3 H8 36,466.62 3,580.08 5,098.33 0.0804 0.19102 n-Butano,
n-C4 H10 59,717.44 4,008.70 1,030.15 0.0905 0.03879 n-Pentano,
n-C5 H12 88,869.64 10,479.08 2,094.47 0.2406 0.08018
n-Hexano, n-
C6 H14 124,892.80 12,422.68 1,489.76 0.2917 0.05833
Heptano,
C7 H16+ 405,776.11 14,927.89 895.10 0.3785 0.03785
Luego se calculan los términos A y B son las ecuaciones 6.40 y 6.41, esto para cada
fase.
Para la fase líquida:
lb
146
Para la fase vapor:
lb
lb
)
𝐵 𝑅𝑇 𝑅𝑇 2𝑙𝑏𝑎𝑏𝑠𝑓𝑡3 .
𝑝𝑔
(10.732 )(640°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑀𝑎 = ∑𝑛𝑗=1 𝑦𝑗𝑀𝑗,
Fase líquida:
𝑙𝑏 𝑙𝑏𝑚
Fase vapor:
𝑙𝑏 𝑙𝑏𝑚
𝑝𝑀 (2800𝑝𝑔2𝑎𝑏𝑠)(32.34𝑚𝑜𝑙) 𝒍𝒃𝒎
𝜌𝑔 = 𝑧𝑅𝑇 = 𝑝𝑔
2𝑙𝑏𝑎𝑏𝑠𝑓𝑡3 = 𝟑𝟐. 𝟗𝟔𝟕𝟓 𝒇𝒕𝟑 . (0.3999)(10.732
)(640°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
Ejercicio 6.7
Para una mezcla de hidrocarburos que existe en equilibrio en dos fases a
𝟐𝟓𝟎𝟎 𝒍𝒃/𝒑𝒈𝟐𝒂𝒃𝒔 y 𝟏𝟓𝟎°𝑭, cuya composición se muestra en la Tabla 6.9, calcular
las densidades de las fases vapor y líquido usando la ecuación de Peng-
Robinson.
Considerar que no existe interacción binaria entre los componentes de la mezcla.
Tabla 6.9 Composición de la mezcla del ejercicio 6.7.
Componente xj yj
(fracción mol líquido) (fracción mol vapor)
Solución:
148
Se deben conocer los valores de la temperatura y presión crítica, así como los
factores acéntricos de cada componente, estos se leen de la tabla del ejercicio anterior,
después se comienza a calcular 𝑇𝑟, m, α, 𝑎𝑐 y b, con las ecuaciones 3.7, 6.37, 6.29,
6.32 y 6.33. 𝑇𝑟 = 𝑇 𝑇𝐶 𝑚 = 0.7464 + 1.522𝜔 − 0.26992𝜔2 𝛼 = [1 + 𝑚(1 − 𝑇𝑟1/2)]2
𝑎𝑐 = Ω𝑎 𝑅 2 𝑇𝑐2
𝑃 𝑏 = Ω𝑏 𝑅 𝑇 𝐶
𝑃 𝑐 Ω𝑎 =0.457264 Ω𝑏=0.077796
𝑐
Se calcula 𝑎𝑇𝑖 con la ecuación 6.28. Luego 𝑎𝑇 y 𝑏 para la fase líquida con las
ecuaciones 6.42, 6.43 y para la fase vapor con las ecuaciones, 6.44 y 6.45
respectivamente.
1/2
𝑎𝑇𝑖𝑗 = ∑ ∑ [𝑥𝑖𝑥𝑗(𝑎𝑇𝑖𝑎𝑇𝑗) (1 − 𝛿𝑖𝑗)] 𝑏 = ∑(𝑥𝑗𝑏𝑗)
1/2
𝑎𝑇𝑖𝑗 = ∑ ∑ [𝑦𝑖𝑦𝑗(𝑎𝑇𝑖𝑎𝑇𝑗) (1 − 𝛿𝑖𝑗)] 𝑏 = ∑(𝑦𝑗𝑏𝑗)
149
Tabla 6.11 Calculo de constantes del ejercicio 6.7 (continuación).
LIQUIDO VAPOR
Ahora para cada fase se calculan los términos 𝐴 y 𝐵 son las ecuaciones 6.34 y 6.35:
Para la fase líquida:
lb
lb
(∑
)𝐵 𝑅𝑇 𝑅𝑇 2𝑙𝑏𝑎𝑏𝑠𝑓𝑡3 .
𝑝𝑔
(10.732 )(610°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
150
lb
)
𝐵 𝑅𝑇 𝑅𝑇 2𝑙𝑏𝑎𝑏𝑠𝑓𝑡3 .
𝑝𝑔
(10.732 )(610°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
𝑚𝑜𝑙.
𝑙𝑏𝑚
𝑙𝑏 𝑙𝑏𝑚
𝑝𝑔
𝜌𝑔 = 𝑧𝑅𝑇𝑝𝑀 = (2800𝑝𝑔2𝑎𝑏𝑠 )2𝑙𝑏(𝑎𝑏𝑠31.𝑓𝑡5𝑚𝑜𝑙3 ) = 87.22 𝒍𝒃𝒎 𝒇𝒕𝟑 .
151
(0.6243)(10.732 )(610°𝑅) 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
de gases a 𝟏𝟓𝟎 °𝑭, cuya composición se muestra en la Tabla 6.13 y que está
contenida en un recipiente de 𝟎. 𝟎𝟗 𝒇𝒕𝟑. Considerar que hay interacción binaria
entre los componentes de la mezcla y son: 0.02 entre metano y n-butano, 0.01
entre etano y n-butano, 0.0 metano y etano. Usar las ecuaciones de Soave-
RedlichKwong.
Tabla 6.12 Composición de la mezcla del ejercicio 6.8
Componente 𝒚𝒋
Fracción mol de vapor
Metano, C1H4 0.5
Etano, C2 H6 0.35
n-Butano, n-C4 H10 0.15
Solución:
𝑎𝑇
[𝑝 + (𝑉𝑀 + 𝑏)] (𝑉𝑀 − 𝑏) = 𝑅𝑇 ,
𝑉𝑀
por lo que es necesario calcular el volumen molar:
0.09 𝑓𝑡3 𝑓𝑡3
𝑉𝑀 = 𝑙𝑏𝑚= 150.845 𝑙𝑏𝑚 .
0.00568
𝑙𝑏𝑚 𝑚𝑜𝑙 𝑙𝑏𝑚 𝑚𝑜𝑙
Ahora se calcula 𝑎𝑇 con la ecuación 6.28 𝑎𝑇 = 𝑎𝑐𝛼, para cada fracción por lo que se
deben conocer la propiedades críticas de cada componente y posteriormente obtener
𝑇𝑟, m, α, 𝑎𝑐 y b, con las ecuaciones 3.7, 6.29, 6.30, 6.32 y 6.33 respectivamente.
𝑎𝑐 = Ω𝑎 𝑅 2 𝑇𝑐2
𝑃 𝑏 = Ω𝑏 𝑅 𝑇 𝐶
𝑃 𝑐 Ω𝑎 =0.42748 Ω𝑏=0.08664
𝑐
152
La Tabla 6.13 resume las operaciones realizadas:
𝑐1𝐻4 666.4 343.33 0.01 1.78 0.496 0.27 8,708.74 0.479 6,068.48
𝑐2𝐻6 706.5 549.93 0.09 1.11 0.632 0.93 21,074.79 0.7238 19,681.29
𝑛 − 𝑐4𝐻10 550.6 765.62 0.19 0.8 0.786 1.18 52,415.21 1.2929 61,640.44
𝑎𝑇𝑖𝑗 1
𝑏𝑖 = ∑(𝑦𝑗𝑏𝑗) ,
𝑓𝑡3
𝑏𝑖 = ∑(𝑦𝑗𝑏𝑗) = (0.5 ∗ 0.479) + (0.35 ∗ 0.7238) + (0.15 ∗ 1.2929) = 0.68678 .
𝑙𝑏𝑚 𝑚𝑜𝑙
Despejando p de la ecuación 6.27 y sustituyendo:
𝑅𝑇 𝑎𝑇
𝑝= − ,
(𝑉𝑀 − 𝑏) 𝑉𝑀(𝑉𝑀 + 𝑏)
𝑓𝑡 𝑙𝑏
𝑝𝑔 𝑎𝑏𝑠
( 10.732 𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅 ) (610°𝑅) 𝑝𝑔 2 𝑎𝑏𝑠 𝑓𝑡 6
15,603.74 (𝑙𝑏𝑚 𝑚𝑜𝑙 )2
3 𝑓𝑡 3
− 3 3 𝑓𝑡 3
( 150.845 −0.68678 ) 150.845 ( 150.845 +0.68678 )
𝑙𝑏
32
𝑝=,
𝑓𝑡 𝑓𝑡 𝑓𝑡
𝑙𝑏𝑚 𝑙𝑏𝑚 𝑚𝑜𝑙 𝑙𝑏𝑚 𝑙𝑏𝑚 𝑙𝑏𝑚 𝑚𝑜𝑙
𝑙𝑏𝑚 𝑚𝑜𝑙 𝑙𝑏𝑚 𝑚𝑜𝑙 𝑙𝑏𝑚 𝑚𝑜𝑙
𝑙𝑏
𝑝 = 42.91 𝑝𝑔 2𝑎𝑏𝑠.
Ejercicio 6.9
Calcular el volumen molar de una mezcla de gases cuya composición se
muestra
Primero se deben obtener los valores de las constantes a, b y c con las ecuaciones
6.59, 6.60 y 6.61 respectivamente. Los resultados de los cálculos realizados se
resumen en la Tabla 6.15.
154
Tabla 6.15 Solución de ejercicio 6.9
Componente Fracción 𝑎𝑗 𝒚𝒋𝒂𝒋 𝑏𝑗 𝒚𝒋𝒃𝒋 𝒄 𝒙𝟏𝟎−𝟔 𝒚𝒋𝒄𝒋
Mol
Metano 0.75 0.297 0.223 -0.2542 -0.191 11.986 8.990
Etano 0.15 0.939 0.141 0.3068 0.046 84.08 12.612
Propano 0.1 1.173 0.117 0.6877 0.069 112.1 11.210
𝐴 𝑦𝑗𝐴0𝑗1/2)2 𝐵 𝑦𝑗𝐵0𝑗
Tabla 6.16 Solución del ejercicio 6.9 (continuación).
(10.732 )(580°𝑅)
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
( )
𝑏𝑝 (−0.076)(280 𝑝𝑔2𝑙𝑏𝑎𝑏𝑠)
(580°𝑅)
𝑙𝑏𝑚 𝑚𝑜𝑙 °𝑅
( )
155
=
𝑐𝑝 (32.818𝑥106)(280 𝑝𝑔2𝑙𝑏𝑎𝑏𝑠) 𝐸 𝑅𝑇4 =
( )
𝒇𝒕𝟑
𝑉𝑀 = 𝟐𝟑 .
𝒍𝒃𝒎 𝒎𝒐𝒍
Capítulo 7.
156
Propiedades físicas de los fluidos de los yacimientos
petroleros
7.1.2 Aceite
El gas puede encontrarse disuelto dentro del aceite o libre dependiendo de las
condiciones de presión del yacimiento, su composición dependerá de las condiciones
de presión y temperatura a las cuales se formó, su composición puede ser determinada
mediante un análisis cromatográfico.
7.2 Propiedades del gas seco.
donde:
3
𝑓𝑡 @ 𝑐.𝑦.
𝐵𝑔: Factor de volumen del gas (𝑓𝑡3@ 𝑐.𝑒.)
Figura 7.1. Comportamiento del 𝑩𝒈, con respecto a la presión del yacimiento y a
temperatura constante.
Es el cambio fraccional del volumen del gas por unidad de presión, su expresión
matemática se puede expresar en función del volumen (𝑉), del volumen molar (𝑉𝑀) y del
volumen específico (𝑣):
158
𝜕𝑉 ............................................. (7.2)
𝑐𝑔 = − ( ) ,
𝑉 𝜕𝑝 𝑇
𝜕𝑉𝑀 ........................................... (7.3)
𝑔 𝑣
𝑐𝑔 = − () ,
𝜕𝑝 𝑇
𝑉𝑀 𝜕𝑝 𝑇
𝜕𝑣 ........................................... (7.4)
Otra 𝑐 =− ( ).
variante
de la expresión para 𝐶𝑔 con los términos de la ecuación de los gases reales es:
159
𝑝𝑝𝑟 𝑧 𝜕𝑝𝑝𝑟 𝑇𝑡𝑟
donde
160
𝑣: viscosidad cinemática(𝑐𝑒𝑛𝑡𝑖𝑠𝑡𝑜𝑘𝑒𝑠, 𝑐𝑠)
El comportamiento de la viscosidad en función de la presión se representa en la Fig.
7.4, de la cual se concluye que cuando la temperatura incrementa, se tienen dos casos;
a presiones bajas la viscosidad del gas se incrementa y a presiones altas la viscosidad
del gas disminuye.
𝑦𝑗𝑀𝑗
Es necesario conocer la composición de la mezcla, así como las viscosidades de
cada gas que la conforma a una presión y temperatura establecidas, de igual manera se
pueden emplear correlaciones gráficas, la Fig. 7.5 representa la presión de ciertos
gases hidrocarburos a presión atmosférica.
161
Figura 7.5: Viscosidad de gases hidrocarburos a presión atmosférica (Mc Cain, 1990).
por la viscosidad del gas (a una atmosfera) y el resultado será la viscosidad del gas en
cuestión a una presión alta.
162
Figura 7.6: Viscosidad de gases hidrocarburos a presión atmosférica considerando la
presencia de 𝑯𝟐𝑺, 𝑵𝟐 y 𝑪𝑶𝟐 (Mc Cain, 1990)
Figura 7.7 Relación de viscosidades para gases naturales con densidades relativas de 0.56 a
0.9
163
Figura 7.8 Relación de viscosidades para gases naturales con densidades relativas de 0.9 a
1.2
Figura 7.9 Relación de viscosidades para gases naturales con densidades relativas de 1.2 a
1.5.
164
Figura 7.10 Relación de viscosidades para gases naturales con densidades relativas de 1.5 a
1.7
Es la relación entre la densidad del aceite y la densidad del agua, ambas a mismas
condiciones de temperatura y presión:
𝜸 𝒐 = 𝜌𝑜
.......... ........................................ (7.13)
, 𝜌𝑤
donde
165
Es el volumen de aceite a condiciones de yacimiento necesario para obtener un barril
de aceite a condiciones atmosféricas:
𝐵𝑂 = 𝑉 𝑉𝑂𝑂 @
@ 𝑐. 𝑦.
𝑐 .𝑒 . , . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (7.15)
donde
𝐵𝑙𝑠 𝑑𝑒 𝑎𝑐𝑒𝑖𝑡𝑒 @ 𝑦
𝐵𝑂: Factor de volumen de formación (𝐵𝑙𝑠 𝑑𝑒 𝑎𝑐𝑒𝑖𝑡𝑒 @ 𝑒)
Figura 7.11 Comportamiento del 𝑩𝑶: en un aceite negro y en función de la presión del
yacimiento.
𝑅𝑆 = 𝑉 𝑔
@ 𝑝 𝑦 𝑇 @ 𝑐.𝑒
𝑉 𝑜 @ 𝑒. . . . . .. .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
(7.16) Donde:
3
𝑓𝑡 @ 𝑐. 𝑒.
𝑅𝑆: Relación de solubilidad (𝐵𝑙𝑠 @ 𝑐. 𝑒.)
166
𝑉𝑔 @ 𝑐. 𝑒.: Volumen del gas producido en la superficie @ estándar
...................... (7.18)
𝐵𝑔 = 𝑉 @ 𝑦 . 𝐵𝑙𝑠
𝑉𝑔𝑔 @ 𝑒. ( 𝐵𝑙𝑠
𝑑𝑒
𝑑𝑒
𝑔𝑎𝑠 @ 𝑦 .
𝑔𝑎𝑠 @ 𝑒. ), . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (7.19)
𝑒 . ( 𝐵𝑙𝑠
𝑅𝑠𝑏 = 𝑉 @
𝑉𝑔𝑜 @ 𝑒 . 𝐵𝑙𝑠 𝑎𝑐𝑒𝑖𝑡𝑒 ). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
𝑔𝑎𝑠 @ 𝑒 .
@ 𝑒. (7.20)
167
El término 𝑅𝑠𝑏 − 𝑅𝑠 es el volumen de gas libre que existe en el yacimiento. La Fig.
7.13 muestra el comportamiento del Bt con respecto a la presión.
Figura 7.13: Comportamiento del factor de volumen total de formación con respecto a la
presión (Danesh, 1998).
168
Figura 7.14: Comportamiento del 𝒄𝒐 en función de la presión y a temperatura constante,
cuando la presión del yacimiento está por encima de la presión de burbuja.
Cuando se tienen dos fases, esta propiedad puede ser determinada mediante 7.26:
𝑐𝑜 = − 𝐵 1𝑜 [(𝜕𝜕𝑝𝐵𝑂)𝑇 − 𝐵𝑔 (𝜕 .
𝜕𝑝𝑅𝑠)𝑇] ................................ (7.26)
169
7.3.7 Coeficiente de viscosidad del aceite 𝝁𝒐.
Se observa que:
a) Para presiones mayores a la de burbuja, la viscosidad aumenta conforme la
presión aumenta.
7.4.1 Composición.
Se conoce como salmuera, por los sólidos disueltos en ella (Na Cl en su mayoría), la
cantidad de sólidos presentes en la salmuera va de 10,000 hasta 300,000 ppm (gr
sólidos por cada 106 gr de agua de formación).
171
𝑚. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (7.28) 𝜌𝑤 = 𝑉 .
Representa el volumen de gas que puede ser disuelto por cada volumen de agua a
condiciones indicadas de presión y temperatura pero medidos a condiciones estándar.
Figura 7.18: Comportamiento del 𝑹𝒔𝒘 con respecto a la presión del yacimiento.
si:
• 𝑝𝑦 ≥ 𝑝𝑏, entonces 𝑅𝑠𝑤 es constante.
172
• 𝑝𝑦 ≤ 𝑝𝑏, entonces 𝑅𝑠𝑤 decrece.
𝐶𝑤 = − 𝐵 1𝑤 [(𝜕𝜕𝑝𝐵𝑤)𝑇 − 𝐵𝑔 (𝜕𝜕𝑝𝑅𝑠𝑤)𝑇].
.............................. (7.33)
7.4.6 Coeficiente de viscosidad del agua de formación.
173
Figura 7.20: Comportamiento de la viscosidad contra la temperatura.
Ejercicio 7.1
Solución:
Asumiendo que tienen un mayor comportamiento como gas ideal, el término 1 (𝜕𝑧) es
𝑧 𝜕𝑃 𝑇
igual a cero, por lo tanto la ecuación 7.4 puede ser expresada como:
1
−6𝑝𝑠𝑖−1.
𝑐𝑔 = 2000 𝑝𝑠𝑖 = 500 𝑥10
sustituyendo
𝑇𝑝𝑐 = 168 + 325(. 72) − 12.5(. 72)2 = 395.5°𝑅 .
175
Figura 7.21: Estimación del factor z de la figura de Standing y Katz
Fig. 4.1:
𝜕𝑧
( ) = −0.022 .
𝜕𝑝𝑝𝑟 𝑇𝑡𝑟=1.52
176
Ejercicio 7.2
el cuál tiene una presión promedio de 𝟐, 𝟎𝟎𝟎 𝒑𝒔𝒊𝒂 y una temperatura de 𝟏𝟐𝟎°𝑭, la
densidad relativa de este gas es de 𝟎. 𝟕𝟐. Calcular el gasto de gas en pies cúbicos
estándar por día (scf/día).
Solución:
Para iniciar se calculan las propiedades críticas con las ecuaciones 4.33 y 4.34, de
donde se obtiene: 𝑇𝑝𝑐 = 395.5 °𝑅, 𝑃𝑝𝑐 = 668.4 𝑝𝑠𝑖𝑎, luego se calculan las propiedades
El factor 𝑍 se estima con la Fig. 4.1 de Standing y Katz, de donde se obtiene 𝑍 = 0.78
𝑓𝑡3 1 𝑓𝑡3 1
𝑉. @ 𝑐. 𝑠.
15,000 @ 𝑐. 𝑦. ( ) = 15,000 @ 𝑐. 𝑠. ( ) ( )
𝑑í𝑎 𝐵𝑔 (𝑉𝑉..@@ 𝑐𝑐..𝑦𝑠..) 𝑑í𝑎 𝐵𝑔 𝑉. @ 𝑐. 𝑦.
15,0000 𝑓𝑡3 𝑄𝑔 = 6
.4𝑥10−3 𝑑í𝑎
@ 𝑐. 𝑠.
𝒔𝒄𝒇
𝑄𝑔 = 𝟐. 𝟑𝟒𝟑𝟕𝟓𝑴𝑴 𝒅í𝒂 .
Ejercicio 7.3
177
Calcular la densidad específica y la densidad en grados API de un sistema de
estándar.
Solución:
Ejercicio 7.4.
Calcular el coeficiente de compresibilidad isotérmico de un gas seco con una
densidad relativa de 𝟎. 𝟗 a una temperatura y presión de yacimiento de 𝟐𝟓𝟎° 𝑭 y
𝒍𝒃
𝟐, 𝟎𝟎𝟎 𝟐 .
𝒑𝒈 𝒎𝒂𝒏
Solución:
178
Figura 7.22: Propiedades pseudocríticas para el gas del problema 7.4.
De donde se obtiene: 𝑇𝑝𝑐 = 442 °𝑅, 𝑝𝑝𝑐 = 639 𝑝𝑠𝑖, posteriormente se calculan las
𝑐𝑝𝑟𝑇𝑝𝑟 = .54.
179
Despejando 𝑐𝑝𝑟 se tiene:
0.54 0.54
Ejercicio 7.5
En cierto pozo se produjo un gas cuya composición se muestra en la Tabla 7.2,
a 𝟐𝟓𝟎°𝑭 y a presión atmosférica. Calcular para dicha mezcla de gases su
viscosidad.
Tabla 7.2 Composición de la mezcla del ejercicio 7.5
Componente Fracción mol
𝒚𝒋
Metano 0.7
Etano 0.16
Propano 0.08
Butano 0.06
Solución:
180
Sustituyendo se tiene:
𝜇𝑔 = = 𝟎. 𝟎𝟏𝟏𝟑 𝒄𝒑.
Ejercicio 7.6.
Del pozo Kuna en México se extrajo una muestra de gas seco a una presión de
𝒍𝒃
𝟏𝟗𝟎𝟎 𝟐 , una temperatura de 𝟐𝟐𝟎°𝑭 y una densidad relativa de 𝟏. 𝟏. Calcular la
𝒑𝒈 𝒂𝒃𝒔
viscosidad de la muestra.
Solución:
𝑙𝑏𝑚
𝑀𝑎 = (𝑀𝑎𝑖𝑟𝑒)(𝛾𝑔) = (28.96)(0.9) = 26.064 (𝑙𝑏𝑚 ).
𝑚𝑜𝑙
mediante la Fig. 7.6, se determina la viscosidad para dicho gas de donde se tiene 𝜇𝑔
= 0.114, las propiedades pseudocríticas para dicho gas se calculan con la Fig. 7.22 de
donde se obtiene 𝑇𝑝𝑐 = 465 °𝑅, 𝑃𝑝𝑐 = 610 𝑝𝑠𝑖𝑎, luego se calculan las propiedades
3.11. Debido a que la densidad relativa del gas es de 1.1 se emplea la Fig. 7.8 para
encontrar una relación de viscosidades.
181
Figura 7.24: Relación de viscosidades para el ejercicio 7.6.
de donde se obtiene 𝜇
𝜇 𝑔
𝑔 1 = 2, por lo tanto 𝜇𝑔 = 2𝜇𝑔 = 2(0.114) = 𝟎. 𝟐𝟐𝟖 𝒄𝒑.
𝟐𝟐𝟎°𝑭 y cuya
𝒑𝒈 𝒂𝒃𝒔
b) Un comportamiento ideal.
182
a) Si se considera que tiene un comportamiento como gas real se calculan las
4.36
Fig. 7.3 se obtiene 𝐶𝑝𝑟𝑇𝑝𝑟 = 2.3 y despejando 𝐶𝑝𝑟 del resultado obtenido 𝐶𝑝𝑟 = 𝑇 2 .3
𝑝𝑟 =
= 2.19,
𝑝𝑐 627.22 (𝑝𝑔2𝑎𝑏𝑠)
Ejercicio 7.8
Estimar el gas en solución de un aceite cuya densidad relativa es de 𝟎. 𝟖 y se
encuentra a 𝟐𝟓𝟎𝟎 𝒑𝒔𝒊𝒂 de presión.
Solución:
184
Ejercicio 7.9
Se tiene una muestra de un aceite al que se le realizó un análisis PVT, entre los
datos obtenidos están 𝜸𝒈 = 𝟎. 𝟕𝟓, 𝑹𝒔 = 𝟓𝟎𝟐. 𝟗, 𝑹𝒔𝒃 = 𝟕𝟐𝟑 y 𝑩𝒐 = 𝟏. 𝟏𝟗. Calcular el factor de
volumen total (𝑩𝒕), considerando que la muestra se encuentra a 𝟏𝟎𝟎°𝑭 y
a 𝟏𝟕𝟎𝟎 𝒑𝒈𝒍𝒃𝟐𝒂𝒃𝒔.
Solución:
Ejercicio 7.10
En un laboratorio se trabajó con una muestra de 𝟓𝟔. 𝟓𝒇𝒕𝟑 de cierto aceite en una
celda a 𝟐𝟖𝟎𝟎 𝒑𝒈 𝒍𝒃
𝟐𝒂𝒃𝒔 a 𝟏𝟕𝟎°𝑭. En un segundo momento la presión se redujo a
𝒍𝒃
𝟑
. Calcular el
185
𝟏𝟖𝟎𝟎 𝒑𝒈𝟐𝒂𝒃𝒔 y se observó que el volumen incrementó a 𝟓𝟕. 𝟗 𝒇𝒕 coeficiente de
Solución:
𝑐𝑜(𝑝1−𝑝2),
=𝑒
𝑉1
𝑉2
−𝑙𝑛 = 𝑐𝑜(𝑝1 − 𝑝2)
𝑉1
2
𝑉
−𝑙𝑛 𝑐𝑜 = (
𝑝1−𝑉𝑝1
2) ,
sustituyendo
3
57.9𝑓𝑡
𝑙𝑛 56.5𝑓𝑡
3 −𝟓𝒑𝒔𝒊−𝟏.
𝑐𝑜 = − (2800𝑝𝑠𝑖−1800𝑝𝑠𝑖) = 𝟐. 𝟒𝟒 𝒙𝟏𝟎
Capítulo 8
Experimentos PVT
Los experimentos PVT son de gran importancia para la caracterización de los fluidos
y para entender el comportamiento termodinámico de los mismos al cambiar sus
condiciones de presión, temperatura y composición. Para poder medir las propiedades
de los fluidos en el laboratorio, se necesitan conocer dos datos: la presión estática
inicial y la temperatura del yacimiento, estos valores se obtienen con mediciones en el
pozo que se obtienen antes de comenzar la explotación.
186
8.1 Métodos de toma de muestras de los fluidos del yacimiento.
Es necesario obtener una muestra al inicio de explotación del pozo, esto bajo un
programa y tratando de asegurar que la proporción de gas y petróleo sea la misma que
la proporción que existía en el yacimiento antes de comenzar la explotación.
El programa de muestreo tiene por objetivo tomar una muestra del fluido original del
yacimiento, es recomendable adquirirla antes de que la presión de fondo fluyente sea
menor que la presión de saturación, con el fin de asegurar que el fluido se encontrará
en una sola fase y la composición será parecida a la original del yacimiento.
Si se sabe que el yacimiento en su estado original se encuentra bajo la presión de
saturación, es necesario esperar para perforar otros pozos.
Para la selección de métodos de acondicionamiento de pozo y estrategias de
obtención de las muestras, debe considerarse el fluido a producir y su comportamiento
al cambiar sus condiciones de presión, temperatura y volumen, debido a que una vez
que se alcanza la presión de saturación es muy poco probable que la composición
original del fluido sea reproducida en el laboratorio, en especial si se trata de un
yacimiento de gas y condensado, debido a que la composición puede ser muy diferente
del fluido original, esto por los fenómenos de condensación retrograda.
Un aspecto como es la localización del pozo seleccionado para el muestreo debe de
ser hecha cuidadosamente, la composición y las propiedades físicas de los fluidos del
yacimiento pueden variar a lo largo del área del yacimiento, estos pozos deben tener
valores de RGA y densidades relativas parecidas.
En yacimientos con grandes espesores es posible encontrar variaciones grandes en
la composición al cambiar la profundidad, también es importante tomar en
consideración las condiciones del fluido en las vecindades del pozo, esto se debe tomar
en consideración para poder acondicionar el pozo, una buena muestra del yacimiento
puede ser tomada siempre y cuando no se haya alcanzado la presión de saturación,
otra forma de obtener una buena muestra de fluido es cuando se tiene una condición de
flujo pseudo-estacionario, en esta condición el fluido que entra al pozo es muy parecida
al fluido original del yacimiento. Otras complicaciones que pueden aparecer es la
187
producción de dos zonas diferentes, o en pozos desviados, donde solo se permite el
muestreo de los fluidos pesados del yacimiento.
Otras consideraciones que deben de tomarse cuando se va a realizar el muestreo son:
1. No producción de agua
2. Mismos rangos de RGA y densidad del aceite
3. Alto índice de productividad.
4. Flujo estacionario sin cabeceo
Es una buena idea no tomar muestras de pozos cercanos al contacto agua aceite.
8.1.2 Métodos de muestreo.
Consiste en tomar muestras de gas y aceite del separador, con medidas adecuadas
del gasto de estas dos fases y con el pozo estabilizado. Algunas de sus ventajas son:
1. Colecta la muestra directamente de la formación.
2. La muestra no está afectada por la segregación en el pozo.
3. La muestra mantiene las condiciones originales.
8.1.2.3 Muestreo en la cabeza del pozo.
Si se conoce que el fluido esta siempre en una sola fase, este puede ser un método
efectivo y con buenos resultados, algunas ventajas son:
1. Menos costoso y complejo.
2. No se pierde producción.
3. Es el método preferido para yacimientos de gas y condensado.
8.1.3 Preparación el pozo para el muestreo.
188
Antes de un muestreo se debe de asegurar que la presión se mantiene por arriba de
la presión de saturación, esto para acondicionar el pozo, cuyo objetivo es el de
reemplazar los fluidos no representativos cercanos al pozo, esto se hace desplazando
el fluido cercano al pozo con fluido de partes más lejanas del yacimiento, puede
requerirse mucho tiempo para lograr acondicionar el pozo.
Consideraciones respecto al tipo de fluido del yacimiento y condiciones de operación.
Antes de tomar una muestra es necesario conocer
• tipo de fluido está siendo producido, condiciones de
189
Figura 8.1: Obtención de una muestra de fondo, a) cámara de muestreo b) distribución
típica de presiones (DAKE, 1978).
Es un buen método para casi todos los tipos de fluidos, consiste en tomar muestras
de gas y aceite a la salida del separador, con el pozo previa y propiamente
acondicionado, el objetivo es obtener fluidos que sean representativos.
190
Figura 8.2: Obtención de una muestra en superficie (DAKE, 1978).
La muestra obtenida de fluido se transfiere a una celda PVT, donde se realiza una
expansión a la temperatura del yacimiento, la presión se controla mediante una bomba
191
de desplazamiento con mercurio, se mide con un manómetro y se registra. Este
experimento se realiza para todos los tipos de fluidos. Consiste básicamente en llevar la
celda PVT a una presión mayor que la presión inicial del yacimiento, y reducir la presión
en una serie de etapas sin extraer gas de la celda PVT. Esta prueba tiene como objetivo
determinar la presión de saturación a temperatura del yacimiento, y determinar las
propiedades PVT del gas y del aceite a presiones mayores que la presión de
saturación.
Dicho proceso se ilustra en la Fig. 8.3.
Volumen relativo.
Factor de compresibilidad 𝑍 del gas.
Saturación de líquido.
Viscosidad.
Densidad.
192
Compresibilidad de fluidos arriba de psat.
Este experimento se realiza para todos los aceites negros. Consiste básicamente en
que a partir de la presión de saturación y a una temperatura de yacimiento constante se
baja la presión en la celda PVT en una serie de etapas y en cada una de ellas se extrae
todo el gas, el cual se mide y se determina su composición y sus propiedades. La Fig.
8.4 muestra dicho procedimiento.
Esta prueba tiene como objetivo determinar las propiedades del gas y del aceite a
presiones debajo de su presión de saturación.
La densidad y la viscosidad el aceite deberán medirse en todo el proceso.
Datos medidos:
Factor de volumen de aceite.
193
Algunos de estos datos en ocasiones es necesario calcularlos, entonces requerimos
como base conocer: el volumen retirado de la celda (𝑣𝑔 en condiciones estándares), el
volumen de aceite en la celda (𝑣𝑂) y el volumen de aceite residual (𝑣𝑜𝑟).
Con los datos anteriores podemos conocer 𝐵𝑔, 𝐵𝑜𝑑 y 𝑅𝑠𝑑, el subíndice 𝑑 en los últimos
parámetros indica que estos han sido obtenidos de una expansión diferencial.
El factor de volumen para el gas puede obtenerse como:
𝑅𝐵 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (8.2) 𝐵𝑔 = 𝑉𝑔 [𝑆𝑇𝐵],
𝑣 𝑔
𝑉 𝑜𝑟
𝑉𝑜𝑟
𝐵𝑔 = 𝑉𝑣𝑜𝑟𝑔 5.6151[𝑆𝐶𝐹𝑅𝐵 ].
𝑉𝑜𝑟
En unidades de campo. Donde 𝑉𝑔 es el volumen total de gas, el factor de desviación
𝑍 se calcula como:
𝑝𝑇𝑠𝑐𝐵𝑔 𝑝𝐵𝑔 ..................................... (8.4)
𝑧= = 35.37 .
𝑝𝑠𝑐𝑇 𝑇
Se conoce como 𝑅𝑠𝑖𝑑 a la relación gas-petróleo disuelta inicial determinada por medio
del ensayo diferencial. Puede estimarse conociendo la totalidad del gas liberado desde
la presión de burbuja hasta la de atmosférica por:
donde
𝐹 . .............................................. (8.7)
𝑉𝑜𝑟
194
𝐹𝑚á𝑥: es el valor más grande calculado de 𝐹, es decir a condiciones estándar.
El parámetro 𝐹 representa la cantidad acumulada de gas liberado a condiciones
estándar y relativa al volumen de aceite residual. En cada etapa de presión, el gas
remanente disuelto en el petróleo es igual al gas que estaba disuelto inicialmente
menos el gas liberado por la disminución de presión.
Este experimento se realiza para todos los yacimientos de tipo volátil y de gas y
condensado. Consiste básicamente en que a partir de la presión de saturación y a una
temperatura de yacimiento constante, se baja la presión en la celda PVT en una serie
de etapas y en cada una de ellas se extrae solamente una parte del gas, tomando como
referencia el volumen de hidrocarburos que se tiene a la presión de saturación. La Fig.
8.5 muestra el proceso de dicha prueba.
Esta prueba tiene como objetivo determinar las propiedades del gas y del aceite a
presiones debajo de su presión de saturación.
% Mol extraído.
195
Saturación de líquido.
Densidad de líquido.
Factor de encogimiento
196
Densidad relativa del gas
Cada etapa
La prueba de viscosidad del aceite se realiza a presiones por arriba y por debajo de
la presión de saturación en un viscosímetro de tipo canica, el cual ha sido calibrado
previamente.
Los valores de viscosidad a presiones por arriba y por debajo de la presión de
saturación deben de ser las mismas que se tiene en el experimento de separación
diferencial o de agotamiento a volumen constante.
En este equipo solo se mide la viscosidad del aceite, la viscosidad del gas se calcula
mediante correlaciones.
Ejercicio 8.1
197
Se realizó un experimento de expansión en equilibrio a una muestra de aceite,
los resultados se muestran en la Tabla 8.1. Calcular la presión de burbuja de
dicha muestra.
Tabla 8.1 Resultado de la prueba diferencial para el ejercicio 8.1.
Presión (psia) 𝒗𝒕/𝒗𝒐𝒕 Presión (psia) 𝒗𝒕/𝒗𝒐𝒕
4000 0.9856 2781 1.0000
3800 0.9880 2600 1.0302
3600 0.9903 2400 1.0550
3400 0.9928 2200 1.0852
3200 0.9948 2000 1.1139
3000 0.9971
Solución:
La presión de burbuja será aquella en la que la relación 𝑣𝑡 = 1, esto es a 2781 psia.
𝑣𝑜𝑡
Ejercicio 8.2
A una muestra de aceite, se realizó un experimento de expansión en equilibrio,
los resultados se muestran en la Tabla 8.2. Calcular
a) Los valores del factor del volumen del gas 𝑩𝒈.
Tabla 8.2 Resultados de la prueba por separación diferencial del ejercicio 8.2.
Presión (psia) 𝒗𝒕/𝒗𝒐𝒓 𝒗𝒈/𝒗𝒐𝒓 𝑽𝒈/𝒗𝒐𝒓
2781 1.2823 0.0000 0.0000
2600 1.2697 0.0476 7.6581
2400 1.2560 0.0491 7.3018
2200 1.2423 0.0509 6.9457
2000 1.2288 0.0547 6.7676
1800 1.2153 0.0595 6.5895
198
1600 1.2021 0.0657 6.4114
1400 1.1889 0.0738 6.2333
1200 1.1760 0.0872 6.2328
1000 1.1632 0.1031 6.0552
800 1.1507 0.1310 6.0548
600 1.1385 0.1726 5.8771
400 1.1267 0.2638 5.8763
200 1.1154 0.5215 5.6991
14.7 1.1060 6.7752 5.3428
S.C. 1.0000 0.0000 0.0000
Solución:
𝑉𝑔
𝐵𝑔 = 𝑉𝑣𝑜𝑟𝑔5.6151[𝑆𝐶𝐹𝑅𝐵 ].
𝑉𝑜𝑟
b) Para calcular 𝑅𝑠𝑑 es necesario agregar otra columna a la tabla para obtener los
𝐹 ,
𝑉𝑜𝑟
Ejercicio 8.3
Solución:
200
𝑉𝑂@ 𝑒. = 𝑉𝑂𝐵 @𝑜 𝑦.
𝟑. =
3
5501.053 𝑐𝑚
= 𝟓𝟐𝟐. 𝟑𝟏 𝒄𝒎
Ejercicio 8.4
Para la muestra del ejercicio anterior, también se registraron valores de 𝑹𝒔,
estos se muestran en la Tabla 8.5. Calcula el volumen de gas producido a
condiciones estándar.
Tabla 8.5 Valores de Rs para el ejercicio 8.4.
Prueba 𝑹𝒔
𝒎𝟑
(𝒎 𝟑𝒔𝒕𝒅)
Diferencial 15.49
Separador 13.45
Solución:
𝑉𝑔 @ 𝑐. 𝑒 = 𝑅𝑆 ∗ 𝑉𝑜@𝑐. 𝑒.
Para el experimento de separación diferencial:
𝑉𝑔 @ 𝑐. 𝑒 = 𝑅𝑆 ∗ 𝑉𝑜@𝑐. 𝑒. = 15.49 ∗ 5.3776𝑥10−4𝑚3 = 𝟖. 𝟑𝟐𝟗𝟗 𝒙𝟏𝟎−𝟑𝒎𝟑.
Capítulo 9.
Correlaciones para estimar las propiedades PVT de los
hidrocarburos.
201
9.1 Importancia de los modelos empíricos.
202
𝑩𝟑 -37.5 -11.1 -3.6
Donde
𝑃𝑝𝑐∗: Presión pseudocrítica corregida (psia).
𝑓(𝜌𝑝𝑟) = 𝑎𝜌𝑝𝑟6 + 𝑏𝜌𝑝𝑟3 + 𝑐𝜌𝑝𝑟2 + 𝑑𝜌𝑝𝑟 + 𝑒𝜌𝑝𝑟3 (1 + 𝑓𝜌𝑝𝑟2 )𝑒𝑥𝑝[−𝑓𝜌𝑝𝑟2 ] ..... (8.9)
−𝑔
𝑓 = 0.6845 . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (9.16)
𝑔 = 0.27𝑝𝑝𝑟 . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (9.17)
Normalmente dicho algoritmo converge rápidamente a no más de cinco iteraciones
para una aproximación aceptable como |𝜌𝑝𝑟𝑘+1 − 𝜌𝑝𝑟𝑘 | < 0.01𝜌𝑝𝑟𝑘 .
𝐵𝑔 = 0.02828 𝑧𝑇 (
𝑝 𝑆𝑇𝐵𝑅𝐵 ). ...................................... (9.19)
O bien en unidades inglesas se utilizan también los miles de pies cúbicos (𝑀𝑆𝐶𝐹)
para expresar el volumen de gas en condiciones estándares. Como 1𝐵 = 5.615 𝐶𝐹,
resulta:
204
𝑧𝑇 𝑅𝐵 (9.20)
𝐵𝑔 = 5.03 (
𝑝 𝑀𝑆𝐶𝐹 ). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
𝜕
Con 𝑇 en °𝑅 y 𝑝 en 𝑝𝑠𝑖𝑎.
9.2.4 Compresibilidad del gas.
Donde:
𝜕𝑧 1 5 2 2 2 .. (9.22)
𝜌𝑝𝑟 = 𝜌𝑝𝑟𝑇𝑝𝑟 [5𝑎𝜌𝑝𝑟 + 2𝑏𝜌𝑝𝑟 + 𝑐𝜌𝑝𝑟 + 2𝑒𝜌𝑝𝑟(1 + 𝑓𝜌𝑝𝑟 − 𝑓2𝜌𝑝𝑟4 )
𝑒𝑥𝑝(−𝑓𝜌𝑝𝑟2 )].
Puede ser estimada con la correlación de Lee (1966) representada por la ecuación
9.23:
𝜇 = 10−4𝐾𝑒𝑥𝑝[𝑋𝜌𝑌]. ....................................... (9.23)
separador.
Donde:
𝛾𝑔𝑠: Densidad relativa a obtenerse si la presión de separación fuera de 100 psig.
donde
206
API ≤30° API >30°
A 0.0362 0.0178
B 1.0937 1.1870
C 25.724 23.9310
Se puede calcular con la ecuación 9.29, solo aplica para presiones menores a la de
burbuja.
𝐴𝑃𝐼 ...................................... (9.29)
𝑅𝑠 = 𝑝𝐵𝐴𝛾𝑔𝑠𝑒𝑥𝑝 [ 𝐶].
𝑇
Esta correlación de Vázquez y Beggs funciona muy bien para crudos de API>15.
Para crudos más pesados, se comporta mejor la correlación de Standing (1977) y
puede ser utilizada para estimar 𝑝𝑏o 𝑅𝑠𝑖, y está representada por la ecuación 9.30.
𝑝𝑏 + 25.48 1.205
. . . . . . . . . . . . . . . . . . . . . . . . . . . (9.30) 𝑅𝑠𝑖 = 𝛾𝑔
(18.2𝑥100.00091𝑇−0.0125𝐴𝑃𝐼) .
Donde:
𝑅𝑠 Relación de solubilidad en SCF/STB.
𝑇: Temperatura en °𝐹.
𝑃: Presión en 𝑝𝑠𝑖𝑎.
9.3.4 Factor volumen del aceite.
207
𝑅𝑠 se mide en (𝑆𝑇𝐵 𝑆𝐶𝐹
), la temperatura en (°𝐹) y 𝐵𝑂 en (𝑆𝑇𝐵 𝑅𝐵
). Los valores de las
Se usa una correlación distinta dependiendo de si la presión del yacimiento está por
encima o por debajo de la presión de burbuja.
Para 𝑝 ≤ 𝑝𝑏
Propuestas por McCain, Rollins y Villena (1988) son válidas hasta 330 °𝐹 y 5300 𝑝𝑠𝑖.
𝑙𝑛(𝑐𝑜) = −7.633 − 1.497𝑙𝑛(𝑝) + 1.115𝑙𝑛(𝑇) + 0.533𝑙𝑛(𝐴𝑃𝐼) + 0.184𝑙𝑛(𝑅 𝑠𝑖) . . . (9.33)
Los resultados obtenidos presentan un error menor de 10% con presiones mayores
de 500 𝑝𝑠𝑖𝑎 y si la presión es menor que esta, el error es menor al 20%.
Cuando se conoce la presión de burbuja la correlación a usar es mejor y se presenta
en la ecuación 9.34:
𝑙𝑛(𝑐𝑜) = −7.573 − 1.450𝑙𝑛(𝑝) − 0.383𝑙𝑛(𝑝𝑏) + 1.402𝑙𝑛(𝑇)
+ 0.256𝑙𝑛(𝐴𝑃𝐼) + 0.449𝑙𝑛(𝑅𝑠𝑖) . . . . . . . . . . . . . . . (9.34)
Estas correlaciones tienen valides hasta 5300 𝑝𝑠𝑖 y 330°𝐹. Donde 𝑇 (°𝑅) y 𝑝(𝑝𝑠𝑖𝑎).
Para 𝑝 > 𝑝𝑏
Se utiliza la siguiente ecuación propuesta por Vázquez y Beggs (1980):
1 ........... (9.35)
208
𝑐𝑜 = 10 𝑝 (−1433 + 5𝑅𝑠𝑖 + 17.2𝑇 − 1180𝛾𝑔𝑠 + 12.61𝐴𝑃𝐼) ,
5
En donde 𝑇(°𝐹) y 𝜇𝑜𝐷(cp) y se restringe a 5 < 𝐴𝑃𝐼 < 58 y 60 < 𝑇 < 175.
Después se considerará el efecto del gas en solución a la presión requerida, con la
correlación de Beggs y Robinson (1975) mostrada en la ecuación 9.37:
donde:
𝐴 = 10.715(𝑅𝑠 + 100)−0.515 , ................................... (9.38)
donde
𝐵 = 2.6𝑝1.187𝑒𝑥𝑝(−11.513 − 8.98𝑥10−5𝑝). ...................... (9.41)
209
9.4 Correlaciones para estimar las propiedades del agua.
donde:
𝐵𝑤(𝑅𝐵/𝑆𝑇𝐵) y 𝑇(°𝐹).
Para agua sin gas las constantes son:
𝑎 = 0.9947 + 5.8𝑥10−6𝑇 + 1.02𝑥10−6𝑇2, ....................... (9.43)
210
Con 𝑝 en (𝑝𝑠𝑖).
Es necesario corregir la salinidad, esto con la correlación:
Cuyas unidades de 𝜇𝑤(𝑐𝑝), 𝑇𝑅(°𝑅) y las correcciones por salinidad, 𝑆𝑐2 y 𝑆𝑝. Se
calculan con:
𝑆𝑐2 = 1 − 0.00187𝑁𝑎 𝐶𝑙0.5 + 0.000218𝑁𝑎𝐶𝑙2.5 +
(𝑇0.5 − 0.0135𝑇)(0.00276𝑁𝑎𝐶𝑙 − 0.000344𝑁𝑎𝐶𝑙 1.5) , ............. (9.53)
.
𝑆𝑝 = 1 + 3.5𝑥10−12𝑝2(𝑇 − 40), ............................... (9.54)
211
tenerse en cuenta que la temperatura del separador no debe ser mayor a 140°𝐹 para la
aplicación de la correlación.
Ejercicio 9.1
Usando las correlaciones correspondientes a las propiedades del gas, calcule
para una presión de 𝟑 𝟐𝟓𝟕 𝒑𝒔𝒊𝒂:
a) El factor de volumen del gas.
Solución:
4.34 de donde obtiene: 𝑇𝑝𝑐 = 392.45 °𝑅, 𝑃𝑝𝑐 = 668.75 𝑝𝑠𝑖𝑎 , ahora se calculan los
valores pseudoreducidos con las ecuaciones 4.35 y 4.36 y se obtiene 𝑇𝑝𝑟 = 1.61,
𝑃𝑝𝑟 = 4.87.
212
Luego se obtiene el valor de 𝑍 mediante la correlación gráfica de Standing y Katz de
la Fig. 4.1, de donde se obtiene 𝑍 = 0.85. Este valor se sustituye en la ecuación 9.18 y
resulta:
𝑧𝑇 0.85(633°𝑅) 𝑹𝑩 𝑹𝑩
𝐵𝑔 = 0.02828 𝑝= 0.02828 (3,257 2𝑙𝑏𝑎𝑏𝑠) = 𝟎. 𝟎𝟎𝟎𝟒𝟔𝟕 (𝑺𝑻𝑩 ) = 𝟎. 𝟎𝟎𝟎𝟖𝟑𝟑 (𝑺𝑪𝑭 ).
𝑝𝑔
0.27𝑝𝑝𝑟
𝑧= ,
𝜌𝑝𝑟𝑇𝑝𝑟
0.27𝑝𝑝
𝑟
Ahora se usa la ecuación 9.22 para calcular 𝜕𝑧 . Los valores de las constantes
𝜕𝑝𝑝𝑟
se calculan con las ecuaciones 9.11, 9.12, 9.13, 9.14, 9.15 y 9.16.
𝑐𝑔 = 𝟐. 𝟓𝟏 × 𝟏𝟎−𝟒𝒑𝒔𝒊−𝟏.
213
Ejercicio 9.2
Solución:
La viscosidad del gas se calcula con ecuación 9.23, por lo que se necesita conocer
el peso molecular del gas y la densidad del gas a la presión y temperatura en que se
𝑔𝑟
𝑀 = 𝛾𝑔𝑀𝑎𝑖𝑟𝑒 = (0.71) (28.9 𝑚𝑜𝑙 ) = 20.5 𝑔/𝑚𝑜𝑙.
por lo tanto:
Ejercicio 9.3
Solución:
La densidad relativa del gas se corrige por las condiciones del separador con la
−5
(27) (544.67°𝑅) 𝑙𝑜𝑔 (150 𝑝𝑠𝑖𝑎)] = 0.721
𝛾𝑔𝑠 = 0.71 [1 + 5.912𝑥10
114.7
Es necesario incluir en la relación gas-aceite inicial el gas seco en el tanque,
sumando al 𝑅𝑠𝑖 proporcionado la corrección propuesta por la correlación 9.55. Pero
primero se necesita la densidad relativa del aceite, que puede ser determinada como:
141.5
°𝐴𝑃𝐼 = − 1.31.5,
𝛾𝑂
despejando
141.5 141.5
𝛾𝑂 = ° 𝐴𝑃𝐼+131.5 = 27+131.5 =
0.893,
sustituyendo en la ecuación 9.52:
log 𝑅𝑆𝑇 = 0.3818 − 5.506𝑙𝑜𝑔𝛾𝑜 + 2.902𝑙𝑜𝑔𝛾𝑔 + 1.327𝑙𝑜𝑔𝑝𝑠𝑒𝑝 − 0.7355𝑙𝑜𝑔𝑇𝑠𝑒𝑝
215
log 𝑅𝑆𝑇 = 0.3818 − 5.506𝑙𝑜𝑔(0.893) + 2.902𝑙𝑜𝑔(0.71) + 1.327𝑙𝑜𝑔(150𝑝𝑠𝑖𝑎)
− 0.7355𝑙𝑜𝑔(85°𝐹) = 1.689
10𝑙𝑜𝑔𝑅𝑆𝑇 = 101.689
Ejercicio 9.4
Solución:
Para la relación gas-aceite inicial se observa que la presión de 3 257 𝑝𝑠𝑖𝑎, es menor
que la presión de burbuja, entonces la relación gas-aceite disuelta se calcula con la
ecuación 9.29, las constantes son las de la Tabla 9.2, resultando:
Para el factor de volumen del aceite, de igual manera a presiones menores que la de
burbuja, el factor de volumen del aceite se calcula con la ecuación 9.31, los valores de
las constantes se muestran en la Tabla 9.3 y se obtiene:
216
(−1.811𝑥10−8) (543.73 𝑆𝐶𝐹
) (632.67°𝐹 − 60) (
) = 𝟏. 𝟐𝟖𝟕 𝑹𝑩
.
𝑆𝑇𝐵 0. 𝑺𝑻𝑩
Ejercicio 9.5
Calcular la compresibilidad isotérmica del aceite.
Solución:
𝑆𝐶𝐹
0.256𝑙𝑛(27) + 0.449𝑙𝑛 (618.90 ) = −9.67.
𝑆𝑇𝐵
𝑒𝑙𝑛(𝑐𝑜) = 𝑒−9.67
𝑐𝑜 = 𝟔. 𝟑 × 𝟏𝟎−𝟕𝒑𝒔𝒊−𝟏.
Ejercicio 9.6
Calcular la viscosidad del aceite.
Solución:
Calculando la viscosidad del petróleo sin gas disuelto esto con la ecuación 9.36
log log (𝜇𝑜𝐷 + 1) = 1. 8653 − 0.025086𝐴𝑃𝐼 − 0.5644𝑙𝑜𝑔𝑇
217
𝐴 = 10.715(𝑅𝑠 + 100)−0.515 = 10.715(543.73𝑆𝐶𝐹/𝑆𝑇𝐵 + 100) −0.515 = 0.3832,
la ecuación 9.37.
Capítulo 10.
Clasificación de yacimientos petroleros de acuerdo al
tipo de fluidos.
esta última clasificación toma especial importancia puesto que determinará la estrategia
Los yacimientos que existen de acuerdo al tipo de fluido son de aceite negro, aceite
volátil, gas y condensado, gas húmedo y gas seco, esta clasificación fue hecha por
existen moléculas largas, pesadas y no volátiles, así como otros diferentes compuestos
químicos.
Su RGA inicial está entre los 3,300 a 5,000 (ft3 de gas a c.e./Bl de aceite a c.e), por
lo que puede producir hasta un 35% de líquido del volumen del yacimiento.
temperatura, así como la localización del punto crítico del yacimiento en el diagrama
Se clasifican en
de dos fases o yacimiento con casquete de gas, en el cual el gas es la fase vapor
Existen otras clasificaciones que se son basadas en las propiedades que presentan
la composición
en:
Aceites negros, y
Aceites volátiles
del sistema. Una clasificación propuesta por McCain de acuerdo los diagramas de fase
220
Gas y condensado,
Gas húmedo, y
Gas seco.
mantienen constantes.
221
Su densidad relativa es menor a 45°𝐴𝑃𝐼 y desciende hasta un valor mínimo,
verdoso o café.
3
𝑓𝑡 𝑑𝑒 𝑔𝑎𝑠 @ 𝑐.𝑒.
3,300 .
𝐵𝑙 𝑑𝑒 𝑎𝑐𝑒𝑖𝑡𝑒 @ 𝑐.𝑒.
iniciar a descender.
En esta etapa la densidad relativa del aceite en °API decrece ligeramente hasta
222
Cuando la 𝑝𝑦 es menor que la 𝑝𝑏, la densidad relativa del aceite, es igual o
C7+, se encuentran en un rango entre el 12.5 y 20% 𝑚𝑜𝑙. La línea del 12.5% es
Presentan un factor inicial del volumen del aceite de la formación, 𝐵𝑜𝑖, mayor a
𝐵𝑙𝑠 𝑑𝑒 𝑎𝑐𝑒𝑖𝑡𝑒 @ 𝑐.𝑦.
2.0 .
𝐵𝑙𝑠 𝑑𝑒 𝑎𝑐𝑒𝑖𝑡𝑒 @ 𝑐.𝑒.
por debajo de la
223
Los líquidos (condensados) en la superficie presentan un color ligeramente café,
encontrando el fluido del yacimiento en fase vapor (gas) y presenta dos puntos
de rocío.
de la vida productiva del yacimiento. Producen 𝑅𝐺𝐴𝑖 mayores a 15,𝐵𝑙000 𝑑𝑒 𝑓𝑡 𝑎𝑐𝑒𝑖𝑡𝑒3 𝑑𝑒 𝑔𝑎𝑠 @ 𝑐@.𝑒
.𝑐.𝑒..
yacimiento.
se presenta.
224
Tienen 𝑅𝐺𝐴𝑖 por arriba de los 100 𝐵𝑙,000 𝑑𝑒 𝑓𝑡 𝑎𝑐𝑒𝑖𝑡𝑒3 𝑑𝑒 𝑔𝑎𝑠 𝑎 𝑐 .𝑎𝑒. 𝑐.𝑒. .
yacimiento.
a 0.7% 𝑚𝑜𝑙.
10.3.6 Criterios de clasificación.
Las siguientes tablas presentan los criterios para la clasificación de los yacimientos
de acuerdo a diferentes autores.
McCain utiliza los criterios de clasificación mostrados en la Tabla 10.2.
Tabla 10.2 Criterios para la clasificación de hidrocarburos según McCain.
Contenido, 𝑪𝟐 − 𝑪𝟔 , % > 35 26 − 35 18 − 26 1 − 18
Mol.
presentan una fase líquida. Conforme la explotación de los fluidos del yacimiento
y vapor).
observa:
decir, el aceite en el yacimiento contiene demasiado gas disuelto más del que
226
puede disolver. Una caída en la presión del yacimiento liberará gas formando
en el yacimiento, se inicia la liberación del gas disuelto en el aceite para formar las dos
227
aceite a condiciones del yacimiento.
Figura 9.10: Diagrama de fase de presión contra temperatura para un yacimiento de aceite
volátil.
(retrógrado).
De donde se observa:
• El punto crítico se localiza más hacia la izquierda baja del diagrama, resultado de
pesados.
presión inicia la condensación de líquido a partir del gas formando una cantidad
228
Figura 10.3: Diagrama de fase de un gas de condensación retrograda (Clark, 1969).
Figura 10.4: Diagrama de fase de presión contra temperatura para un yacimiento de gas.
húmedo (Clark, 1969).
229
Donde:
gas seco.
Figura 10.5: Diagrama de fases para un yacimiento de gas seco (Clark, 1969).
De donde se observa:
230
• Las trayectorias de producción tanto a condiciones de yacimiento como de
Capítulo 11.
Equilibrio líquido-vapor.
11.1 Definiciones.
231
donde:
𝐾𝑖: Constante de equilibrio del componente 𝑖.
𝑥𝑖: Fracción mol del componente 𝑖 en fase vapor.
𝑦𝑖: fracción mol del componente 𝑖 en la fase líquida.
ideal se mezclan.
Las soluciones ideales no existen, solamente los gases a bajas presiones y los
líquidos de series homologas se aproximan a una solución ideal a bajas presiones.
(11.3)
𝑧𝑛𝑡 = 𝑥𝑖𝑛𝐿 + 𝑦𝑖𝑛𝑉. ...................................... ..
Utilizando la relación:
𝑝 (11.4)
232
𝑦𝑖 = 𝑃𝑣𝑖 𝑥𝑖, .............................................
o también
∑ 𝑦𝑖 = ∑ 𝑧𝑛 𝑡
𝑝 = 1. .................................... (11.5)
𝑛𝑣+𝑝𝑣𝑖𝑛𝐿
ideales.
Con el fin de predecir las constantes de equilibrio de una solución ideal, la ley de
Raoult establece una presión parcial (𝑝𝑖)de un componente en un sistema
multicomponente es el producto de su fracción molar en la fase líquida y la presión de
vapor del componente. 𝑝𝑖 = 𝑝𝑣𝑖 × 𝑥𝑖.
............................................. (11.8)
233
𝑝𝑖 = 𝑝 × 𝑦𝑖. ............................................... (11.9)
o bien:
𝑦𝑖 𝑝𝑣𝑖
. ............................................. (11.11)
En general a presiones inferiores a los 100 𝑝𝑠𝑖𝑎, las leyes de Raoult y Dalton proveen
una forma simple para estimar las constantes de equilibrio. En otro caso, los cálculos
deben realizarse considerando una solución real y las constantes son una función de la
presión, temperatura y composición.
𝐾𝑖 = 𝐾(𝑝, 𝑇, 𝑧). .......................................... (11.12)
Dado que no siempre se cuenta con ecuaciones que representen correctamente a los
fluidos, se emplean correlaciones par estimar la constante de equilibrio a bajas y
moderadas presiones y temperaturas. Para estas correlaciones, se asume que son
independientes de la composición de la mezcla de fluidos. Para una mezcla
multicomponente, los valores de las constantes de equilibrio tenderán a converger a 1,
𝐾𝑖 = 1 para todos los componentes a una cierta presión conocida como la presión de
convergencia.
Donde:
𝑇𝑐 y 𝑝𝑐 son las temperaturas y presiones críticas del componente 𝑖 y 𝜔 es el factor
acéntrico de cada componente. Sin embargo existen tres factores que limitan el uso de
las constantes de equilibrios calculadas con las leyes de Raoult-Dalton y son las
siguientes:
234
a) La mezcla de gases se supone que se comporta como un gas ideal, sin
embargo
este comportamiento se establece solamente a bajas presiones.
b) La ley de Raoult solamente se aplica a soluciones ideales y la mayoría de
las
soluciones se separan de este comportamiento.
c) La mayoría de los componentes gaseosos puros que existen en los
hidrocarburos, se encuentran a presiones mayores que la crítica, por lo cual una
presión verdadera para un componente puro no existe.
11.4 Metodologías para determinar las condiciones de equilibrio de una mezcla.
Los cálculos flash son una parte integral de los procesos realizados en ingeniería de
yacimientos y de procesos, estos son requeridos cuando se desea conocer la cantidad
de moles de líquido y gas coexistiendo en un yacimiento a unas condiciones de presión
y temperatura dadas. Estos cálculos se llevan a cabo para determinar la composición
de las fases de hidrocarburos existentes. Dada la composición de un sistema de
hidrocarburos a una presión y temperatura específicas, los cálculos flash se llevan a
cabo para determinar las moles de la fase gaseosa, 𝑛𝑉, las moles de la fase líquida, 𝑛𝐿,
la composición de la fase líquida, 𝑥𝑖, y la composición de la fase gaseosa, 𝑦𝑖. Los
cálculos computacionales son iterativos y se calcula de la siguiente manera:
∑
𝑓(𝑛𝑣) = 𝑛
𝑧 (𝐾 − 1 )
𝑣 (𝑖𝐾 𝑖−𝑖 1 )+ 1 = 0, ................................ (11.14)
la cual se puede resolver para el número de moles de la fase vapor mediante alguna
técnica iterativa, una vez obtenido el valor de 𝑛𝑣 se puede emplear el método Newton
Raphson como:
(𝑛𝑣)𝑛𝑒𝑤 = 𝑛𝑣 − 𝑓 𝑓( 𝑛 )
.
′ ( 𝑛𝑣𝑣 ) ..................................... (11.15)
235
11.6 EJERCICIOS RESUELTOS
Ejercicio 11.1
Calcular la presión en el punto de burbuja a 𝟐𝟎𝟎 ° 𝑭 de una mezcla cuya
composición se muestra en la Tabla 11.1.
Tabla 11.1 Composición de la mezcla del ejercicio 11.1
Componente Composición
𝑪𝟑 0.56
𝒏𝑪𝟒 0.35
𝒏𝑪𝟓 0.9
Total 1.0
Solución:
𝑝𝑏 = ∑ 𝑧𝑗𝑝𝑣𝑗
𝑗=1
Ejercicio 11.2
236
Empleando la ecuación para cálculo de puntos de rocío, determinar la presión
en cuestión para la mezcla del ejercicio anterior.
Solución:
Primero deben determinarse las presiones de vapor, las cuales pueden ser obtenidas
mediante la carta de Cox de la Fig. 3.5 o mediante alguna correlación para la
temperatura deseada, para posteriormente emplear la ecuación 11.7, los resultados se
muestran en la Tabla 11.3.
𝑛
𝑧𝑗
𝑝𝑑 = 1/ ∑ ( )
𝑝𝑣𝑗 𝑗=1
𝑝𝑑 = 0 .003647 = 𝟐𝟕𝟒. 𝟏𝟖
𝒑𝒈 𝟐𝒂𝒃𝒔 .
Capítulo 12.
Hidratos de gas.
237
Una vez que el agua está cercana a su punto de congelamiento se comienzan a
formar pequeñas cavidades dentro de los enlaces, estos enlaces no serán estables a
menos que las cavidades recién formadas sean llenadas con gas. Estas nuevas
estructuras estabilizadas mediante moléculas de gas, son llamadas hidratos de gas.
Su apariencia es similar a la de la nieve o el hielo, sin embargo permanecen en
estado sólido hasta temperaturas tan altas como 35ºC. Ejemplos de gases que pueden
estabilizar las moléculas del agua son el metano (𝐶𝐻4) y el dioxido de carbono (𝐶𝑂2), la
temperatura de formación de los hidratos incrementa con la presión. Durante los
procesos de transporte de los hidrocarburos hasta las plantas de procesamiento la
caida de temperatura puede ser considerable y esto puede promover la presencia de
hidratos de gas que llevan al taponamiento de las líneas de producción, debido a ello es
necesario considerar la formación de hidratos cuando se está realizando el diseño de
las líneas de transporte (Pedersen, 2007).
La Fig. 12.1 muestra el comportamiento de fases de un hidrato de gas puro, la línea
AB representa las condiciones de equilibrio entre los hidratos en forma gaseosa y hielo.
la línea BC representa la línea de equilibrio entre el gas y el agua en forma líquida,
mientras que la línea CD representa los puntos de equilibrio entre los hidratos de gas en
estado líquido y el agua en forma líquida.
Cuando se trata de una mezcla de gases es posible que la curva de los hidratos pase
a través de la envolvente de fases de la mezcla de gases como es mostrado en la Fig.
12.2 a diferencia del componente puro, en la línea CD se tiene la coexistencia de tres
fases en equilibrio.
238
Figura 12.1 Diagrama de fases para un hidrato de gas puro, G (gas), H, (hidrato), L1 (agua
líquida), L2 (hidrato en estado líquido).
Figura 12.2 Diagrama de fases para un hidrato con una mezcla de gases.
Algunas moléculas pueden entrar dentro de cualquier estructura como es el caso del
metano, sin embargo moléculas tales como el Iso-Butano, sólo podrá entrar en las
cavidades grandes de la estructura II.
12.2 Modelamiento de la formación de hidratos.
240
1. Agua líquida o hielo (𝛼) → Estructra de hidrato vacia (𝛽)
El primer término del lado derecho, representa el efecto de estabilización del hidrato
causado por la adsorción de moléculas de gas, esta diferencia de potencial químici
puede ser calculado mediante la expresión 12.2:
𝑁𝐶𝐴𝑉 𝑁
𝐶𝑘𝑖𝑓𝑘
𝑌𝑘𝑖 = . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (12.3)
1 + ∑𝑁
𝑗 =1
𝐶𝑗𝑖𝑓𝑗
Wilson, sin embargo las estimaciones de los coeficientes de fugacidad serán para
todas las fases excepto para los sólidos puros y los hidratos, el valor es 1.0, con
equilibrio. La fase gas se asume como un gas ideal y se asume que la fase aceite
fugacidades estimadas en 1.
3. Una vez que se tienen las composiciones, se recalculan las fugacidades 𝑓𝑘𝑚𝑖𝑥,𝑘 =
1, 2, … 𝑁
𝜃 = ln 𝑓𝑤𝐻 − ln 𝑓𝑤𝑚𝑖𝑥
𝑁𝐶𝐴𝑉 𝑁
𝐶
= ∑ 𝑣𝑖 ln (1 − ∑ 1 + ∑𝑘𝑖𝑁𝑓𝑘𝐶𝑗𝑖𝑓𝑗) + ln 𝑓𝑤𝐻 ..... (12.4)
242
𝑖𝑘=1
𝑗=1
− ln 𝑓𝑤 .
Los nuevos estimados para los hidratos son calculados a partir de las
fugacidades de la mezcla mediante la ecuación 12.5:
ln 𝑓𝑤𝐻 = ln 𝑓𝑤𝑚𝑖𝑥 + 𝜃, 𝑘
. . . . . . . . . . . . . . . . . . . . . . . . . . . (12.5) = 1, 2, … , 𝑁 .
5. Una vez estimadas las composiciones de los componentes que no forman parte de
12.6.
𝑁𝐶𝐴𝑉
𝑥
𝑁 ;𝑘
𝑥 𝑤𝑘 = ∑ 𝑣𝑖 1 + ∑
𝐶𝑗𝑘𝑖=1
𝑓𝑘𝐶𝑗𝑖𝑓𝑗 . . . . . . . . . . . . . . . . . . . (12.6)
𝑖=1
= 1, 2, … , 𝑁
6. Basados en el comportamiento de las fugacidades, se puede estimar una nueva
7. De no converger, se regresa a 3.
243
Bibliografía
1. AHMED, Tarek, “Equations of State and PVT Analysis: Applications for Improved
Maracaibo, 1996.
3. BEGGS, H. D., and Robinson, J. R., “Estimating the Viscocity of Crude Oil
electrónica,
https://alojamientos.uva.es/guia_docente/uploads/2011/469/45750/1/Estructura_Fo
https://www.academia.edu/6228571/II.GASES_REALES (descargado el 05 de
8. BUECHE, Frederick J., “Física General”, Décima Edición, Mc Graw Hill, México,
2007.
10. CRAFT, B.C. Hawkins, M. and Terry, R.E., “Applied Petroleum Reservoir
12. CUAUTLI, Teresa, “Apuntes de Propiedades de los fluidos”, Tesis para obtener el
14. DANESH, Alí, “PVT and Phase Behavior of Petroleum Reservoir Fluids” Gulf
15. H. WHITSON, Curtis, “Fluid Properties”, Petroleum Engineering Data Book, Mayo,
1994.
245
1986.
APENDICE A.
si 𝐷 > 0, la ecuación solo tiene una raíz real y se obtiene con la expresión:
𝑍 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A.5) o
si 𝐷 < 0, la ecuación tiene tres raíces reales y se calculan con las expresiones
siguientes:
𝑍 , .................................... (A.6)
𝑍 , . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A.7) y
𝑍 . ............................... (A.8)
Donde
si 𝐷 = 0, la ecuación tiene tres raíces reales y dos de ellas son iguales y se calculan
1/3 − 𝑎1
................ ............................. (A.10)
𝑍1 = 2𝐽
3
y
𝑍2 = 𝑍3 = −𝐽1/3 − 𝑎31. ...................................... (A.11)
Apéndice B
class EOS
{
public static void EOSSolver()
{
//mixcomp();
coeficientes();
mixcoef();
EOSpar();
roots(); if
(v.root[3] == 0)
{
v.factorzL = Math.Min(v.root[0], Math.Min(v.root[1], v.root[2]));
}
else
{
v.factorzL = v.root[3];
}
mixcoefy();
EOSpar();
247
roots(); if
(v.root[3] == 0)
{
v.factorzV = Math.Max(v.root[0], Math.Min(v.root[1], v.root[2]));
}
else
{
v.factorzV = v.root[3];
}
}
public static void mixcoef()
{
v.a1 = 0;
v.B1 = 0;
for (int i = 0; i < v.contador; i++)
{ for (int j = 0; j <
v.contador; j++)
{
v.a1 = v.a1 + v.x[i] * v.x[j] * Math.Sqrt(v.arraya[i] * v.arraya[j] * v.aw[i] *
v.aw[j]);
}
v.B1 = v.B1 + v.x[i] * v.arrayb[i];
}
}
public static void mixcoefy()
{
v.a1 = 0;
v.B1 = 0; for (int i = 0; i
< v.contador; i++)
248
{ for (int j = 0; j <
v.contador; j++)
{
v.a1 = v.a1 + v.y[i] * v.y[j] * Math.Sqrt(v.arraya[i] * v.arraya[j] * v.aw[i] *
v.aw[j]);
}
v.B1 = v.B1 + v.y[i] * v.arrayb[i];
}
}
public static void coeficientes()
{
250
}
public static void EOSpar()
{
double delta1= 0, delta2 = 0;
v.A = v.a1 * v.presion * 1 / v.R / v.R / v.Tyac / v.Tyac;
v.B = v.B1 * v.presion / v.R / v.Tyac;
if
(v.global==3)
{ delta1 =2;
delta2 = -1;
} if
(v.global == 2)
{
delta1 = 1;
delta2 = 0;
} if
(v.global == 1)
{
delta1 = 1;
delta2 = 0;
}
//v.alpha = v.B - 1;
//v.betha = v.A - (2 * v.B) - 3 * Math.Pow(v.B, 2);
//v.gamma = -(v.A * v.B) + Math.Pow(v.B, 2) + Math.Pow(v.B, 3);
v.alpha = -(1+v.B-delta1*v.B);
v.betha = v.A + (delta2) * Math.Pow(v.B, 2) -delta1*v.B- (delta1) *
Math.Pow(v.B, 2);
v.gamma = -(v.A*v.B+delta2*Math.Pow(v.B,2)-delta2*Math.Pow(v.B,3));
}
public static void roots()
251
{ for (int i = 0; i <
4; i++)
{
v.root[i] = 0;
}
v.p = v.betha - Math.Pow(v.alpha, 2) / 3;
v.q = v.gamma + 2 * Math.Pow(v.alpha / 3, 3) - v.alpha * v.betha / 3;
v.disc = Math.Pow(v.q / 2, 2) + Math.Pow(v.p / 3, 3);
if (v.disc <
0)
{
v.m = 2 * Math.Sqrt(-v.p / 3);
v.o = 3 * v.q / v.p / v.m;
v.tetha = Math.Acos(v.o) / 3;
for (int i = 0; i <= 2; i++)
{
v.root[i] = v.m * Math.Cos(v.tetha + (i * 2f / 3f) * Math.PI) - v.alpha / 3;
}
}
else
{
}
}
}
}
253