Petrleo

Modelo matemtico para la estimacin del porcentaje de asfaltenos en petrleos

crudos como funcin de la gravedad API, temperatura y viscosidad
Mathematical model for estimating the percentage of asphaltens in crude oil as a
function of the API gravity, temperature and viscosity

Modelo matemtico para a estimativa da percentagem de asfaltenos em petrleo como

funo da gravidade API, temperatura e viscosidade

Toms Marn1 Mara Albornett2 Alirio Romero3

Recibido: 4-3-17; Aprobado: 14-6-17
Resumen Abstract Resumo
Determinar el porcentaje de asfalte- Determine the asphaltenes per- Determinar a percentagem de asfal-
nos en petrleo es de gran impor- centage in oil is important by the tenos em petrleo de grande im-
tancia por los problemas que cau- problems that cause in the produc- portncia pelos problemas que cau-
san en los sistemas de produccin tion systems when precipitating and sam nos sistemas de produo ao
al precipitar y obstruir tuberas y obstruct pipes and facilities. The precipitar e obstruir encanamentos
equipos. El conocimiento del conte- knowledge of the content of e equipes. O conhecimento do con-
nido de asfaltenos es fundamental asphaltenes is fundamental to pre- tedo de asfaltenos fundamental
para predecir su comportamiento dict its behavior in the oil. We para predizer seu comportamento
en el petrleo. Se busc establecer sought to establish a mathematical no petrleo. Procurou-se estabele-
un modelo matemtico para estimar model to estimate the asphaltenes cer um modelo matemtico para esti-
el porcentaje de asfaltenos en mues- percentage in samples of crude oil mar a percentagem de asfaltenos
tras de crudo de campos petroleros of petroleum fields of the Monagas em mostras de cru de campos pe-
del Norte del estado Monagas, Ve- State North, Venezuela, taking as troleiros do Norte do Estado Mona-
nezuela, teniendo como variables la variables the gravity API, the tem- gas, Venezuela, tendo como vari-
gravedad API, la temperatura y la perature and the viscosity. There veis a gravidade API, a temperatura
viscosidad. Se realizaron medidas was realized measurement of API e a viscosidade. Realizaram-se me-
de API y viscosidad a diferentes tem- and viscosity to different tempera- dida de API e viscosidade a diferen-
peraturas (entre 23 y 40 C) a 27 tures (between 23 and 40 C) to 27 tes temperaturas (entre 23 e 40 C)
muestras de petrleo crudo y el por- samples of crude oil and the a 27 mostras de petrleo cru e a per-
centaje de asfaltenos se obtuvo por asphaltenes percentage it was ob- centagem de asfaltenos obteve-se
la norma ASTM D6560, para luego tained by the norm ASTM D6560, pela norma ASTM D6560, para de-
con el paquete estadstico Stat- for then with the statistical applica- pois com o pacote estatstico Stat-
graphics Centurion XVI determinar tion Statgraphics Centurion XVI to graphics Centurion XVI determinar
la ecuacin que relaciona la varia- determine the equation that relates a equao que relaciona a varivel
ble en estudio. Se plantearon siete the variable in study. Seven (7) em estudo. Propuseram-se sete (7)
(7) modelos matemticos y se utili- mathematical models appeared modelos matemticos e utilizaram-
zaron como parmetros el coefi- and there was used like parameters se como parmetros o coeficiente
ciente de correlacin R2, el error es- the coefficient of interrelation R2, de correlao R2, o erro regular m-
tndar medio absoluto y el error rela- the absolute average standard er- dio absoluto e o erro relativo mdio.
tivo promedio. Se obtuvo que el mo- ror and the average relative error. It Obteve-se que o modelo que me-
delo que mejor se ajust a los datos was obtained that the model that lhor se ajustou aos dados foi o mode-
fue el modelo exponencial en tem- better fitted to the information was lo exponencial em temperatura, du-
peratura, doble raz cuadrada en the exponential model in tempera- pla raiz quadrada em API e quadra-
API y cuadrado Y - logaritmo X en ture, double root squared in API and da Y - logaritmo X em viscosidade
2 2
viscosidad con R de 0,996, error es- squared Y - logarithmic X in viscos- com R de 0,996, erro regular mdio
2
tndar medio absoluto de 0,17 y ity with R of 0.996, absolute aver- absoluto de 0,17 e erro relativo m-
error relativo promedio de 5,1%. age standard error of 0.17 and aver- dia de 5,1%.
age relative error of 5.1 %.
Palabras clave/Keywords/Palavras-chave:
Asfaltenos, modelo matemtico, porcentaje, precipitacin, regresin, asphaltenes, mathematical model, percentage, precipitation,
regression, percentagem, precipitao, regresso.
1
IngQum. MSc. Departamento de Ingeniera de Petrleo. Escuela de Ingeniera y Ciencias Aplicadas. Universidad de Oriente,
Ncleo de Monagas (UDO). e-mail: tmarin@udo.edu.ve.
2
IngPet. UDO. e-mail: mariaalboragui@hotmail.com.
3
IngPet. UDO. e-mail: alirio_romero1992@hotmail.com.
GEOMINAS, Vol. 45, N 73, agosto 2017 73
T. Marn, M. Albornett, A. Romero

Introduccin crudo. Los mtodos espectrom- Procesamiento de Hidrocarburos,

Los asfaltenos son los componen- tricos requieren de costosos equi- en total se tuvieron 27 muestras. A
tes ms pesados del petrleo cru- pos y su resultado depende de cur- cada muestra de petrleo crudo se
do, fsicamente son slidos de es- vas de calibracin, aunque tienen le determin la gravedad API por el
tructura amorfa, solubles en disol- la ventaja de ser mtodos rpidos, mtodo del Termohidrmetro, si-
ventes aromticos como el bence- su resultado est influenciado por guiendo el procedimiento estable-
no e insolubles en disolventes ali- las otras especies qumicas pre- cido en la norma ASTM D287 y se
fticos como el n-heptano. sentes en el petrleo que tambin corrigieron a temperatura estndar
Se han realizado diferentes estu- absorben luz visible. (60 F) a travs de las tablas dispo-
dios para determinar la forma en la En la bsqueda de metodologas nibles en la norma ASTM D1250.
que se encuentran los asfaltenos alternativas que permitan estimar Luego se determin el porcentaje
dentro del petrleo, siendo la ms con un mnimo margen de error el de asfaltenos de cada una de las
aceptada la teora coloidal pro- porcentaje de asfaltenos en el pe- muestra por precipitacin con n-
puesta por Kawanaka et al. (1989), trleo crudo, se han establecido heptano, de acuerdo a lo estableci-
la cual establece que los asfalte- modelos matemticos como el de do en la norma ASTM D6560. Para
nos se encuentran en suspensin Aliendres et al. (2013), el cual esti- la determinacin de las viscosida-
coloidal dentro del petrleo, esta- ma el porcentaje de asfaltenos en des de las muestras se utiliz un
bilizados por otros hidrocarburos funcin a la gravedad API, la tem- viscosmetro rotacional Brookfield
como las resinas. El sistema en peratura y el ndice de Refraccin, (ASTM D2196) y un bao termos-
suspensin es termodinmica- un modelo sencillo y que utiliza ttico en el cual se coloc la tem-
mente inestable por lo que cam- una mnima cantidad de muestra, peratura en 30, 35 y 40 C. Para la
bios en las variables presin, tem- sin embargo requiere del uso del xi- temperatura de 23 C, las viscosi-
peratura y composicin hace que leno como disolvente y se basa, al dades fueron determinadas sin el
la suspensin se rompa provocan- igual que el anlisis espectrofoto- uso del bao termosttico (tempe-
do la separacin de los asfaltenos mtrico en la capacidad de los as- ratura promedio del laboratorio) y
de la fase lquida, lo que produce faltenos para absorber luz visible, para la temperatura de 26 C se co-
depsitos de estos, que obstruyen teniendo ambos mtodos la misma locaron las muestras en bao de
tuberas y equipos. desventaja en cuanto a que no so- agua sin calentamiento (tempera-
El contenido de asfaltenos consti- lo los asfaltenos tienen esta pro- tura promedio del agua en el labo-
tuye un factor importante en la de- piedad dentro del petrleo. ratorio).
terminacin de los trayectos de pro- En base a lo anterior en la presente Luego de obtener los datos nece-
cesamiento y refinacin de un pe- investigacin se plantea el desa- sarios, se tabularon y se cargaron
trleo crudo (Akbarzadeh et al., rrollo de un modelo alternativo ba- en el programa estadstico Stat-
2007), por lo que es de suma im- sado en otra propiedad que es di- graphics Centurion XVI, para ini-
portancia su cuantificacin, pues rectamente proporcional al conte- cialmente establecer los modelos
de su cantidad depende la rapidez nido de asfaltenos como lo es la vis- matemticos que relacionan a la
con la cual ocurren los tapona- cosidad, pues al ser los asfaltenos variable dependiente (%ASF) con
mientos en los sistemas de pro- slidos suspendidos una mayor cada una de las variables indepen-
duccin y procesamiento del pe- cantidad de estos aumenta la re- dientes (API y VISC) por separado,
trleo. sistencia al flujo, pudindose rela- utilizando la funcin de regresin
La determinacin del porcentaje cionar entonces ambas variables. simple en su modalidad compara-
msico de asfaltenos en el petr- Para hacer el modelo ms confia- cin de modelos alternativos. Una
leo crudo se realiza comunmente ble se adicionaron, adems las va- vez obtenidos los modelos que
mediante dos metodologas espe- riables gravedad API para relacio- ms se ajustaron al comporta-
cficas, como lo son: el anlisis gra- nar el contenido de asfaltenos con miento de los datos, tomando co-
vimtrico establecido en normas la densidad y la temperatura, debi- mo parmetro de seleccin el coe-
como las ASTM D3279 y D6560 y do a que la viscosidad es depen- ficiente de determinacin R2 se pro-
anlisis espectrofotomtrico, el diente de esta variable. cedi a establecer siete (7) mode-
cual se basa en la relacin directa los matemticos combinando los
entre la absorcin de luz de los as- Metodologa anteriores y adicionando la varia-
faltenos y su cantidad (ASTM La investigacin se desarroll en el ble temperatura (T), segn lo esta-
D7996). Ambos mtodos tienen Laboratorio de Procesamiento de blecido por Aliendres et al. (2013);
ventajas y desventajas, los basa- Hidrocarburos del Departamento se utilizaron las funciones de re-
dos en gravimetra requieren de de Ingeniera de Petrleo de la Uni- gresin mltiple lineal y regresin
cierto tiempo de ensayo (alrededor versidad de Oriente, Ncleo de Mo- mltiple no lineal, obtenindose
de 9 horas) y medidas de masa y nagas. Para su desarrollo, se pro- los coeficientes de cada ecuacin.
volumen que conlleva a la acumu- cedi, en primer lugar, a la obten- Los parmetros de seleccin del
lacin de errores, adems del uso cin de las muestras de petrleo mejor modelo fueron el coeficiente
de disolventes como el h-heptano crudo de campos productores del de determinacin R2, el error es-
y el tolueno, sin embargo, permite rea norte del estado Monagas; tndar medio absoluto y el error re-
determinar el porcentaje de asfal- los mismos fueron identificados y lativo promedio entre los valores
tenos por separacin directa del almacenados en el laboratorio de experimentales y los predichos por
74 GEOMINAS, Vol. 45, N 73, agosto 2017
 Modelo matemtico para la estimacin del porcentaje de

cada modelo planteado.

Modelos de regresin multivariables evaluados
Resultados En la tabla IV se muestran los modelos de regresin
Propiedades de las muestras de Petrleo Crudo multivariables que se establecieron para la relacin
En la tabla I se muestran los resultados obtenidos, %ASF vs API, VISC y T.
de acuerdo a lo establecido en la metodologa.
Modelos matemticos obtenidos a
Tabla I. Propiedades de las muestras de petrleo crudo. partir de regresin mltiple
Los resultados de los coeficientes
de cada uno de los modelos mate-
mticos evaluados, luego del uso
del programa Statgraphics Centu-
rion XVI, se muestran en la tabla V.

Resultados de parmetros esta-

dsticos obtenidos para cada mo-
delo
Los resultados de los parmetros to-
mados en consideracin para la se-
leccin del modelo matemtico se
muestran en las tablas VI a XI.

Discusin
En la tabla I se muestran los resulta-
dos de la caracterizacin de las 27
muestras de petrleo crudo perte-
necientes a los campos Furrial, Pun-
ta de Mata y Orocual del norte del
estado Monagas. Se tiene que el
22,22 % de la muestra corresponde
a crudos Livianos, los cuales son
provenientes del campo Carito en el
rea de Punta de Mata; 59,26 % de
la muestra es de crudos Medianos,
los cuales provienen del campo El
Furrial y el restante 18,52 % de la
muestra es de crudos Pesados del
campo Orocual. Se observa una
aparente relacin directa entre el
porcentaje de asfaltenos y la visco-
sidad de las muestras, es decir, que
a mayor viscosidad se tiene un ma-
yor porcentaje de asfaltenos, lo que concuerda con lo
Modelos de regresin de %ASF vs API y %ASF establecido por Ghanavati, et al. (2013), en que hay
vs VISC un incremento de la viscosidad al incrementar el con-
En las tablas II y III se muestran los modelos de regre- tenido de asfaltenos a una temperatura constante.
sin con mejor ajuste segn el coeficiente de determi- Las viscosidades de las muestras decrecen con la
nacin entre la variable dependiente y las dos varia- temperatura, siendo tambin consistente con Gha-
bles independientes de manera individual. navati et al. (2013), siendo mayores las variaciones
Tabla II. Modelos de regresin simple entre %ASF y las variables API. en los crudos pesados. La razn est
en la inmensa variedad de molcu-
las o de familias de molculas que
se encuentran en cualquier crudo, la
cual aumenta exponencialmente
cuando este es ms pesado (Len,
Tabla III. Modelos de regresin simple entre %ASF y VISC. 1998).
En las tablas II y III se observan los re-
sultados de los anlisis de regresin
entre el %ASF y las variables API y
VISC, en donde se muestra que no
existe una relacin lineal entre las
mismas, es decir, los modelos son de
GEOMINAS, agosto 2017 75
T. Marn, M. Albornett, A. Romero

Tabla IV. Modelos matemticos evaluados.

Tabla V. Modelos matemticos obtenidos.

Tabla VI. Resultados de los parmetros estadsticos para la seleccin. tipo no lineal, las correlaciones en-
tre las variables son fuertes, siendo
la variable API la de mayor influen-
cia sobre el %ASF. El comporta-
miento no lineal de la relacin
%ASF vs API tambin fue estableci-
do por Aliendres et al. (2013), cuyo
modelo asumi una relacin expo-
nencial directa entre estas varia-
bles.
Los diferentes modelos matemti-
cos establecidos se muestran en la
tabla IV, comenzando con un mode-
Tabla VII. Resultado de anlisis ANOVA de los errores absolutos medios. lo lineal mltiple,
que a pesar de
que la relacin
entre las varia-
bles no es lineal,
se tom como
modelo base; los
otros 6 modelos
Tabla VIII. Resultados del contraste mltiple de rangos para error se disearon tomando en cuenta la rela-
absoluto medio. cin no lineal entre las variables, sin em-
bargo, para el caso de la temperatura, se
estableci lineal en los modelos 2, 4 y 5,
ya que en el trabajo de Aliendres et al.
(2013) se estableci que la temperatura
tena una relacin lineal en la ecuacin.
Una vez obtenidos los 7 modelos mate-
mticos mostrados en la tabla V, se pro-
cedi a la seleccin del mejor, para lo
cual se tomaron en cuenta los parme-
tros establecidos en el estudio, mostra-
dos en la tabla VI. Se observa que los mo-
delos que mejor definen la variabilidad
del %ASF son los enumerados con 5, 6 y
7, por ser sus R2 los mayores, esto coincide con los menores errores tanto absolutos como relativos. Los an-
lisis ANOVA realizados a los errores tanto absolutos (Tabla VII) como relativos (Tabla IX), mostraron que exis-
76 GEOMINAS, Vol. 45, N 73, agosto 2017
 Modelo matemtico para la estimacin del porcentaje de

Tabla IX. Resultado de anlisis ANOVA de los errores relativos medios. te diferencia significati-
va entre ellos (P-Valor <
0,05), indicando que los
modelos tienen compor-
tamientos diferentes al
predecir el %ASF. Lo an-
terior queda demostra-
Tabla X. Resultados del contraste mltiple de rangos para do con el contraste ml-
tiple de rangos mostrado en las tablas VIII y X; en
la tabla VIII se observa que respecto a los errores
absolutos, existen tres grupos homogneos, des-
tacndose el primer grupo integrado por los mode-
los 5, 6 y 7, es decir, que respecto a los errores ab-
solutos no existe diferencia significativa entre los
resultados de estos tres modelos, siendo los de
menor error. Respecto al error relativo, se obser-
va en la tabla X, que el nico modelo que presen-
t un comportamiento no significativo fue el mode-
lo 1, entre los dems modelos no existe diferencia
significativa, destacndose nuevamente como
los de menor error los modelos 5, 6 y 7.
Para poder seleccionar un modelo de los 3 que mejor comportamiento mostraron, se realiz una validacin
externa, con 6 muestras adicionales de crudos medianos del campo Furrial (Tabla XI) para de esta forma se-
leccionar el modelo que mostrara una mejor prediccin del %ASF en muestras no utilizadas para la creacin
del modelo. Se observa que en funcin al error relativo porcentual medio, el modelo que mejor predijo el com-
portamiento de %ASF en las muestras fue el modelo 6, el cual fue de este modo el seleccionado en esta in-
vestigacin.

Conclusiones
El porcentaje de asfaltenos en crudos de los campos productores del norte de Monagas puede ser estimado
mediante el modelo matemtico propuesto, en funcin a la Gravedad API, la Viscosidad y la Temperatura, el
cual tiene la siguiente expresin:
ASF = e(1,38731+0,0102598*T)*((1,38344 + -0,146861*(API1/2))2)* (0,606354 + 0,198331*ln(VISC))2
para unos rangos de aplicacin de: API entre 11 y 37; Viscosidad entre 0,8 y 380 cP y temperatura entre 23 y
40 C; con un porcentaje de error relativo medio de 5,1 %.

Tabla XI. Resultados de validacin externa de los modelos 5, 6 y 7.

Referencias
Akbarzadeh K, Allenson, S, Creek, J, Jamaluddin,A, Marshall, A, Mullins, O, Solbakken, T. (2007, agosto).
Los asfaltenos: Problemticos pero ricos en potencial. Oileld Review, p. 443-458.
Aliendres, M, Rendn, M, Marn T. (2013). Modelo matemtico para la estimacin del porcentaje de asfalte-
nos en crudo en funcin de la gravedad API, la temperatura y el ndice de refraccin. Revista RITIUDO,
Vol 1 Nro 1, p. 52 62.
ASTM D287.(2012). Standard Test Method for API Gravity of Crude Petroleum and Petroleum Products (Hy-
drometer Method), USA: ASTM International, West Conshohocken, PA.
ASTM D1250. (2013). Standard Guide for Use of the Petroleum Measurement Tables, USA: ASTM Interna-
GEOMINAS, agosto 2017 77
T. Marn, M. Albornett, A. Romero

tional, West Conshohocken, Insolubles) in Crude Petro- Temperature on Viscosity of

PA. leum and Petroleum Prod- Iranian Heavy Crude Oil: Ex-
ASTM D2196. (2015). Standard ucts, USA: ASTM Interna- perimental and Modeling
Test Methods for Rheological tional, West Conshohocken, Study,. Revista Energy Fuels,
Properties of Non-Newtonian PA. Vol 27 Nro 12. p 72177232.
Materials by Rotational ASTM D7996 (2015). Standard Len, V. (1998). Nuevo enfoque so-
Viscometer, USA: ASTM In- Test Method for Measuring bre la visin molecular de un
ternational, West Visible Spectrum of crudo pesado. Revista Visin
Conshohocken, PA. Asphaltenes in Heavy Fuel Tecnolgica, Vol 5 Nro 2. p
ASTM D3279. (2012). Standard Oils and Crude Oils by Spec- 131 138.
Test Method for n-Heptane troscopy in a Microfluidic Plat- Kawanaka S, Leontaritis, K.J,
Insolubles, USA: ASTM Inter- form, USA: ASTM Interna- Park S.J, Mansoori, G.A.
n a t i o n a l , W e s t tional, West Conshohocken, (1989, julio). Thermodynamic
Conshohocken, PA. PA. and Colloidal Models of
ASTM D6560 (2012). Standard Ghanavati, M, Shojaei, M-J, Asphaltene Flocculation. Oil-
Test Method for Determina- Ramazani, A. (2013). Effects Field Chemistry. ACS Sym-
tion of Asphaltenes (Heptane of Asphaltene Content and posium Series, Vol. 396, p.

GEOMINAS

Fuente: https://www.emaze.com/@AOQCFZQR/Impactos-en-la-
Industria-Petrolera-copy1

Fuente: http://www.crezcamos.com/blog-11-m/86-conozca-el-impacto-
ambiental-de-los-derrames-de-petryoacuteleo.htm
78 GEOMINAS, Vol. 45, N 73, agosto 2017