República Bolivariana de Venezuela

Ministerio del Poder Popular para la Educación Superior

Universidad de los Andes
Facultad de Ingeniería – Escuela de Ingeniería Química
Departamento de Química Industrial
Caracterización del crudo pesado venezolano Boscan
Empresa Petrol CYBSA: Sarali Quintero, Álvaro Ramírez, Jenny Picón, Carlos Torres, Balmore
Zambrano
Debido a la compleja composición de los crudos, es necesario realizar una caracterización del
mismo con el fin de conocer un poco más a fondo sus propiedades químicas y físicas, las cuales
son de gran importancia no solo para saber en qué tipo de refinería puede llevarse a cabo su
tratamiento, sino también para su comercio. Ahora bien, con el objetivo de realizar la
caracterización del crudo Boscan, por medio de la página de John Jechura se pudo obtener los
datos del estudio del crudo Boscán 1 representados en la tabla 1, mientras que los cálculos
expuesto en la presente entrega están fundamentados con el capítulo tres del libro Taher A.
Alsahhaf, Amal Sayed Elkilani, M.A. Fahim - Fundamentals of Petroleum Refining.
Tabla 1. Propiedades del crudo

Gravedad API 10,9

Azufre (%Peso) 5,5
Punto de inflamación (ºC) 57
Densidad (g/mL) @ 15 ºC 0,993
Punto de fluidez (ºC) 10
Viscosidad dinámica (mPa.so Sc) @ 15 ºC 485500
Viscosidad de Saybolt (SUS) 90000

Tabla 2. Distribución del Punto de Ebullición

Temperatura (ºC) %Peso

120 1
140 1
160 2
180 2
200 3
250 6,6
300 9,9
350 15
400 21
450 28
500 36
550 46
600 55
 650 64
700 73

Con los datos expresados en la tabla 2 es posible obtener la curva que pertenece a la
norma ASTM D 2887 Método de prueba estándar para la distribución del rango de temperatura de
ebullición de fracciones de petróleo por Cromatografía de gases, la cual expresa la relación entre
una distribución de temperaturas de ebullición versus el porcentaje en peso del destilado.
Obteniendo así la curva expuesta en la figura 1 realizando una simulación de destilación con el
software Aspen Hysys.

Figura 1. Curva ASTM D 2887

Es posible convertir la curva ASTM D 2887 a la curva ASTM D 86 por medio del software
antes mencionado con el objetivo de encontrar de manera más precisa los valores de las
temperaturas de ebullición con respecto al porcentaje de volumen de destilado, obteniendo de tal
forma la figura 2.
Figura 2. Curva

Ahora bien, de la curva ASTM D 86 es posible obtener los siguientes valores:

Tabla 3. Distribución de Temperatura versus el porcentaje en volumen de destilado

%Volumen Temperatura (ºC) Temperatura (ºF)

0 150 302
10 330 626
30 485 905
50 565 1049
70 660 1220
90 745 1373
95 850 1562

Con los datos de la tabla 3 es posible hacer la conversión de la curva ASTM D 86 a la curva
de TBP, la cual muestra información más detallada sobre la volatilidad del crudo o sus fracciones 2.
Dicha conversión es posible mediante la ecuación propuesta por Riazi y Daubert (1980) 3.

TBP=a ( ASTMD 86)b Ecuación 1

Teniendo en cuenta que los factores de a y b son constantes que varían según el
porcentaje volumétrico que se tiene de destilado y la temperatura está expresada en grados
Rankine, como puede verse en la tabla 4 extraída del capítulo tres del libro Taher A. Alsahhaf,
Amal Sayed Elkilani, M.A. Fahim - Fundamentals of Petroleum Refining 4.
 Tabla 4. Constantes de Riazi y Daubert para la conversión de la curva ASTM D 86 a TBP

%Volumen a b
0 0,9167 1,0019
10 0,5277 1,0900
30 0,7429 1,0425
50 0,8920 1,0176
70 0,8705 1,0226
90 0,9490 1,0110
95 0,8008 1,0355

Sustituyendo los valores en la ecuación 1, se obtienen los siguientes datos:

Tabla 5. Resultados obtenidos de la conversión de la curva ASTM D 86 a TBP

%Vol a b T1 (ºC) T1(ºR) T2(ºR) T2(ºC)

0 0,9167 1,0019 150 761,67 707,081442 119,67
10 0,5277 1,09 330 1085,67 1074,72595 323,92
30 0,7429 1,0425 485 1364,67 1377,8333 492,31
50 0,892 1,0176 565 1508,67 1530,7433 577,26
70 0,875 1,0226 660 1679,67 1738,28939 692,57
90 0,949 1,011 745 1832,67 1889,05416 776,32
95 0,8008 1,0355 850 2021,67 2121,22651 905,31

Posterior a estos cálculos es posible graficar la curva TBP y compararla con la de ASTM D
86 como es ilustrado en
la figura 3.
 Figura 3. Curvas ASTM D 86 y Curva TBP

Los rangos de temperatura de ebullición del punto de quiebre de la curva pueden dar
información de los pseudocomponentes, los cuales son fracciones de crudo que no pueden ser
puntualmente identificados5, por tanto se agrupan y se usa la temperatura promedio del corte
NBP. La cantidad de estos dependen del rango de temperatura de ebullición del crudo.

Tabla 6. Número de cortes según los rangos de temperatura

Rango TBP Número de cortes

<37.8 ºC (100) ºF Componentes actuales
37.8-427 ºC (100-800) ºF 28
427-649 ºC (800-1200) ºF 8
649-871 ºC (1200-1600) ºF 4

Figura 4. TBP - NBP

Con la finalidad de continuar con la caracterización del crudo, es posible hallar el factor de
caracterización de Watson, Kuop el cual permitirá evaluar la parafinidad o aromaticidad de las
fracciones destiladas, pero antes es necesario hacer el cálculo del punto de ebullición de volumen
promedio (VABP), por tanto se sugieren las siguientes ecuaciones:
3
MeABP
Kuop= √ Ecuación 2
SG
MeABP=VABP−∆ Ecuación 3 Donde MeABP es la temperatura de ebullición promedio en
grados Rankine.

T 10+T 30 +T 50 +T 70+T 90
VABP= Ecuación 4
5
Todas las temperaturas deben estar expresadas en ºF
0.6667 0.333
ln ∆=−0.94402−0.00865 ( VABP−32 ) + 2.99791 SL Ecuación 5
 T 90−T 10
SL= Ecuación 6
90−10
141.5
ºAPI = −131.5 Ecuación 7
SG

Sustituyendo los valores en cada una de las ecuaciones, se tiene la siguiente tabla de
resultados:

Tabla 7. Resultados de los cálculos

SL 9.3375
∆ 89.7906
VABP (ºF) 1034.6
MeABP (ºR) 1404.4794
SG 0.9937
Kuop 11.2698

Teniendo el valor del Kuop es posible definir la clasificación del crudo Boscán por medio de
la siguiente tabla

Tabla 8. Clasificación de crudos de acuerdo al factor de caracterización Kuop

Hidrocarburos Kuop

Parafínicos >11.5

Nafténicos 10.5-11.5

Aromáticos <10.5

Como el Kuop obtenido se encuentra entre los valores de 10.5 y de 11.5 se puede
considerar al crudo Boscán como un crudo nafténico.
 Continuando
con la caracterización
del crudo, se procede a
calcular la gravedad
específica por medio de

pseudocomponentes, haciendo uso de los volúmenes expresados en la tabla 3 y de acuerdo a la

ecuación 8 es posible realizar una gráfica donde pueda observarse API-TBP versus el %Volumen de
destilado (ver figura 4).

API=−0,0002∗( %Vol )3 +0,0238∗( %Vol )2−( 1,7637∗%Vol ) +80,471 Ecuación 8

Figura 4. Gráfica API-TBP

El cálculo del peso molecular se realiza usando la ecuación de Pedersen, la cual expresa:

M =42.965 ⌊ exp ( 2.097 x 10−4 Tb−7.78712 SG+ 2.08476 X 10−3 TbSG ) Tb1.26007 SG 4.98308 ⌋
Ecuación 9

Teniendo en cuenta que las temperaturas deben estar expresadas en grados Kelvin y sería
el valor de la temperatura MeABP que es la temperatura media igual a 780.266 ºK y el valor de SG
conocido de 0.9937, se tiene el siguiente resultado:

M =474.6032
Para calcular la viscosidad del crudo se tiene una expresión dada por Pedersen, en la cual
se determina la viscosidad para fracciones liquidas de petróleo comprendidas en un rango de
(37.8-98.9) ºC que es equivalente a rango de (100-210) ºF. Se tienen dos variables, V 210 y V100 son
las viscosidades calculadas a 210 y 100 ºF respectivamente. Con los valores de Kuop de 11.2698 y
la gravedad API de 10.9, se calcula:

−4 2 8.0325 x 10−2 K +1.248

−3
logV 210=−0.463634−0.166532 ( API ) +5.13447 x 10 ( API ) −8.48995 x 10 K ( API )+
API +26.7
Ecuación 10

2 0.17161 K 2+ 10.9
log V 100 =4.39371−1.94733 K + 0.127690 K 2 +3.2629 x 10−4 ( API ) −1.18246 x 10−2 K ( API )+

Ecuación 11

Obteniendo unos valores de V210=28.8334 cSt y V100 = 4628.2546 cSt

Índice de Refracción es una importante propiedad física fundamental que puede ser usada
en conjunto con otras propiedades para caracterizar hidrocarburos puros y sus mezclas.
1
1+2 I
n= ( 1−I ) 2

Ecuación 12

I =a exp ( bTb+ cSG+dTbSG ) T eb S G f

Ecuación 13

Tabla 9. Constantes para la ecuación 13

Constantes Light fractions Heavy fractions

Rango de Peso molecular 70-300 300-600
Rango de Pto. de ebu. (ºF) 90-650 650-1000
a 2.266x10-2 2.341x10-2
b 3.905x10-4 6.464x10-4
c 2.468 5.144
d -5.704x10-4 -3.289x10-4
e 0.0572 -0.407
f -0.720 -3.333

Realizando los cálculos indicados, se obtiene los siguientes valores I= 0.3254 y n= 1.5643.

Como último paso en esta caracterización se deben determinar el tipo de composición

molecular que tiene el crudo, esto es para conocer específicamente el número de parafinas,
naftenos y aromáticos que posee el crudo, esto es usado para pesos moleculares mayores a 200, R
es la intercesión de la refractividad y VGC es una constante para viscosidad.

x p=a+ b ( R ) +c (VGC ) Ecuación 14

x n=d + e ( R ) +f ( VGC ) Ecuación 15

x a=g+ h ( R )+i ( VGC ) Ecuación 16

Tabla 10. Constantes para las ecuaciones 14, 15 y 16

Constantes Light fractions Heavy fractions

Rango de Peso molecular 70-200 200-600
a -13.359 2.5737
b 14.4591 1.0133
c -1.41344 -3.573
d 23.9825 -2.464
e -23.333 -3.6701
f -0.81517 -1.96312
g -9.6235 -4.0377
h -8.8739 2.6568
i -0.59827 1.60988

D
R=n− Ecuación 17 Hace referencia al corte de refracción
2
10 SG−1.0752 log ( V 100 −38 )
VGC 1= Ecuación 18
10−log ( V 100 −38 )

SG−0.24−0.022 log ( V 210 −35.5 )

VGC 2= Ecuación 19
0.755

Sustituyendo los valores y realizando los cálculos expuestos, se obtienen los siguientes resultados:

R= 1.0678 y VGC1 = 0.9466

Xp=0.2735

Xn= 0.4065 Lo cual tiene sentido, puesto que con el Kuop se concluyó que era un crudo Nafténicos.

Xa= 0.32

Es posible realizar los mismos cálculos para el ejemplo 2B4.1, página 228 del Manual API
sin Pass, Technical Data Book- Petroleum Refining. Por tanto se desea calcular el tipo de
distribución molecular de una fracción de petróleo de 0.9046 de gravedad específica, un índice de
refracción de 1.5002, una densidad líquida de 0.90, un MeABP de 798ºF y V100 de 336.

Con la ecuación 7 es posible conocer los grados API de esta fracción de crudo.

141.5
ºAPI = −131.5=24.923
0.9046
Empleando la ecuación 9 es posible determinar el peso molecular de la muestra,
recordando que la temperatura se debe sustituir en grados kelvin, por tanto se tiene un MeABP de
698.71ºK.
M =42.965 ⌊ exp ( 2.097 x 10−4 (698.71)−7.78712(0.9046)+2.08476 X 10−3 (697.71)(0.9046) ) 698.711.26007 0.904
Obteniendo un valor de M =377.305

Conocido el valor del peso molecular es posible determinar el valor del VGC (viscosity
gravity constant) por medio de la ecuación 18 antes expuesta.

10 (0.9046)−1.0752 log (336−38 )

VGC 1= =0.84851
10−log ( 336−38 )
Hallando la refractividad de la muestra por medio de la ecuación 17, se tiene:

D 0.90
R=n− =1.5002− =1.0502
2 2
Teniendo todos los datos necesario para emplear las ecuaciones 14, 15 y 16, se obtienen la
siguiente distribución molecular:

x p=2.5737+1.0133 ( 1.0502 )−3.573 ( 0.84851 )=0.6061

x a=−4.0377+2.6568 ( 1.0502 )+ 1.60988 ( 0.84851 ) = 0.1185

x n=1− ( 0.6061+ 0.1185 )=0.2754

Ecuación generalizada para propiedades termo físicas: Realizando el cálculo del ejemplo
3.10 del capítulo 3 del libro Taher A. Alsahhaf, Amal Sayed Elkilani, M.A. Fahim - Fundamentals of
Petroleum Refining.
Riazi y Alsathaf en 1996 proponen un método que permite determinar una serie de
propiedades como el NBP, la densidad, el índice de refracción entre otros que expresan en la tabla
11. La ecuación general de la cual se hará uso es:

θ=θ∞−exp ( a−bM c ) Ecuación 20

θ Se define como cualquiera de las propiedades termo físicas que sean de interés, θ∞
(valor tabulado) es el límite evaluado para cuando se asume que M tiende al infinito, siendo M el
peso molecular, mientras que a, b y c son las constantes para cada propiedad.

Tabla 11. Constantes de las propiedades termo físicas.

θ Nombre θ∞ a b c
Tb Rango promedio de 1080 6.97996 0.01964 0.67
ebullición en ºK
SG Gravedad específica 1.07 3.56073 2.93886 0.1
d20 Densidad del líq a 20ºC 1.05 3.80258 3.12287 0.1
I Índice de Refracción 0.34 2.30884 2.96508 0.1
Tbr Tb reducida 1.20 -0.34742 0.02327 0.55
-Pc Presión crítica en bar 0 6.34492 0.72390 0.3
-dc Densidad crítica g/cm3 -0.22 -3.2201 0.00090 1
 -ω Factor acéntrico 0.30 -6.2520 -3.64457 0.1
σ Tensión superf dynes/cm 30.3 17.45018 9.70188 0.1
δ Solubilidad (cal/cm3)2 8.60 2.29195 0.54907 0.3

Ahora bien, para el ejemplo 3.10 se conoce el MeABP de 320ºC y la gravedad API de 34, se
emplea la ecuación 7 donde se define los grados API a 60ºF, con la intensión de conocer la
gravedad específica.

141.5
SG=
ºAPI +131.5
Sustituyendo los valores, se obtiene un valor de gravedad específica igual a 0.8549. Se
emplea la ecuación 9 para conocer el valor del peso molecular sustituyendo los valores conocidos
de gravedad específica y Tb= MeABP= 320ºC= 593.15ºK

M =42.965 ⌊ exp ¿
M =257.139
Teniendo el valor del peso molecular y conociendo los valores de los parámetros
expuestos en la tabla 11, es posible emplear la ecuación 20 con el fin de conocer todas las
propiedades termo físicas antes citadas.

θ θ∞ a b c Valor Obtenido
SG 1.07 3.56073 2.93886 0.1 0.8595
d20 1.05 3.80258 3.12287 0.1 0.8555
I 0.34 2.30884 2.96508 0.1 0.2825
Tbr 1.20 -0.34742 0.02327 0.55 0.3350
-Pc 0 6.34492 0.72390 0.3 12.416
-dc -0.22 -3.2201 0.00090 1 0.2517
-ω 0.30 -6.2520 -3.64457 0.1 0.8012
σ 30.3 17.45018 9.70188 0.1 28.565
δ 8.60 2.29195 0.54907 0.3 8.0566
Tabla 12. Resultados para el ejemplo 3.10

Es posible realizar los mismos cálculos para el crudo Boscán, con el fin de obtener una
caracterización del mismo más completa. Recordando que el valor del peso molecular del mismo
es igual 474.6032 , por consiguiente se obtienen los siguientes resultados:

θ θ∞ a b c Valor Obtenido
SG 1.07 3.56073 2.93886 0.1 0.9177
d20 1.05 3.80258 3.12287 0.1 0.9121
I 0.34 2.30884 2.96508 0.1 0.2985
Tbr 1.20 -0.34742 0.02327 0.55 0.8456
-Pc 0 6.34492 0.72390 0.3 5.7363
-dc -0.22 -3.2201 0.00090 1 0.2461
 -ω 0.30 -6.2520 -3.64457 0.1 1.3448
σ 30.3 17.45018 9.70188 0.1 29.704
δ 8.60 2.29195 0.54907 0.3 8.2975
Tabla 13. Propiedades termo físicas del crudo Boscán