Emulsiones y Demulsificantes
Emulsiones y Demulsificantes
Emulsiones y Demulsificantes
INTRODUCCIN
una aproximacin mediante CFD multifsico del flujo dentro de un tanque tipo FWKO
(Abdulkadir y Hernandez-Perez, 2010). En este caso se calcula el comportamiento
estacionario del fluido en una malla bidimensional, utilizando un modelo Euleriano de dos
fases miscibles, que s pone en evidencia la notable influencia de las fuerzas boyantes.
MODELO MATEMTICO
En el trabajo de Barcel y colaboradores (Barcel et al., 2010) se detalla el desarrollo del
modelo de drift-flux isotrmico a partir del modelo de dos fluidos (Ishii, 1987), ms
completo pero tambin ms costoso desde el punto de vista computacional. El modelo de
drift-flux (Ishii y Hibiki, 2010), tiene como caracterstica principal el considerar a la mezcla
como un nico seudofluido y no como dos fases separadas, e introduce el concepto de
velocidad de deslizamiento para tener en cuenta la diferencia de velocidad entre las fases.
1 Ecuaciones bsicas
Llamaremos fase 1 o primaria al agua, y 2 o secundaria al petrleo. Utilizaremos los
correspondientes subndices para indicar campos de variables primarias y secundarias
tales como las velocidades v1 y v2 y las densidades 1 y 2. Definiremos tambin como
es usual en estos modelos las concentraciones volumtricas 1 y 2 para indicar las
fracciones del volumen ocupadas por agua y petrleo, respectivamente, en un dado punto
del espacio a un dado tiempo. El modelo de Drift-Flux utiliza los conceptos de flujo
volumtrico o velocidad del centro de volumen, definido como
j = 1v1 + 2v2, (1)
y de velocidad media o del centro de masa, definido como
vm = 11v1 + 22v2 m , (2)
donde
m = 11 + 22 (3)
es la densidad media. Definimos tambin la velocidad de deslizamiento vkj , o driftvelocity de la fase k, como la velocidad de dicha fase respecto del centro de volumen, o
expresado de otra manera, relativa al flujo volumtrico j. Es decir
vkj = vk j. (4)
Otras velocidades relativas tiles son las velocidades de las fases relativas al centro de
masa de la mezcla, definidas como
vkm = vk vm, (5)
para las que por definicin (Ishii y Hibiki, 2010) se cumple que
2. La velocidad de deslizamiento
Qued pendiente en el desarrollo anterior, entre otros detalles que veremos ms adelante,
la especificacin de la llamada velocidad de deslizamiento, v2j , utilizada en la Ec.(9). En
el presente modelo se ha optado por una relacin simple en la cual dicha velocidad
depende en cada punto del dominio de la concentracin local secundaria:
v2j = v2j (2). (13)
No se considera en este modelo ninguna dependencia de esta velocidad con la topologa
de la mezcla de fases. La determinacin de esta funcin debiera hacerse en cada caso
particular a travs de experimentos con los fluidos reales que se intenta simular. En el
presente trabajo se propone un modelo simple que intenta capturar la esencia del proceso
de separacin. El modelo propuesto es el siguiente:
La velocidad V0 debe ser interpretada como la velocidad de ascenso de una gota solitaria
de petrleo en un mar de agua, es decir, como la velocidad de Stokes de la gota. El
presente modelo Mecnica Computacional Vol XXXII, pgs. 1761-1783 (2013) equiere por
lo tanto definir un dimetro tpico de gota, a partir del cul se puede calcular este
parmetro V0. Respecto de los otros dos parmetros, debemos recordar que la fase
dispersa es en este caso el petrleo, si bien en el caso de los separadores de fases
estudiados en este trabajo el petrleo aparece en algunas situaciones en el papel de fase
continua, conteniendo agua no separable o emulsionada como fase dispersa. Este
fenmeno est representado en el presente modelo por el lmite 2max, por encima del
cual no es esperable separacin ulterior, ya que para este valor la velocidad de
deslizamiento es nula.
ANLISIS DIMENSIONAL
1. Descripcin de la herramienta
El anlisis dimensional es una de las tres herramientas fundamentales en el anlisis de
problemas fluidodinmicos, junto con el anlisis diferencial y el anlisis integral (o por
balances en volmenes de control). La tcnica de anlisis dimensional adquiere
relevancia sobre todo en la planificacin de experimentos y presentacin de resultados en
forma compacta, pero se utiliza con frecuencia en estudios tericos.
Esta tcnica se basa en que las variables involucradas en cualquier situacin fsica real
pueden ser agrupadas en un cierto nmero de grupos adimensionales (llamados nmeros ) los cuales son suficientes para caracterizar el fenmeno fsico. En general, la
cantidad de grupos adimensionales es menor a la cantidad original de variables
dimensionales. El anlisis dimensional permite adems relacionar los datos medidos en
un modelo experimental con la informacin requerida para el diseo de un prototipo a
escala real, ya que proporciona leyes de escala basadas en las similitudes geomtrica,
cinemtica y dinmica entre modelo y prototipo. Como ejemplo, podemos citar el conocido
problema de la prdida de carga por unidad de longitud, o gradiente de presin, en
tuberas horizontales en flujo completamente desarrollado
Dicho gradiente depende de 5 variables: 2 geomtricas (dimetro y rugosidad absoluta
de la tubera), 1 cinemtica (velocidad media del flujo), y 2 dinmicas (densidad y
viscosidad del fluido). El anlisis dimensional permite reducir el problema y expresarlo
como uno en el cul un gradiente de presin adimensional depende de slo dos variables
adimensionales: el conocido nmero de Reynolds y la rugosidad relativa respecto del
dimetro de la tubera. De esto modo, de 6 variables dimensionales pasamos a slo 3
nmeros , sin reducir en nada la generalidad del problema. Debe quedar claro sin
embargo que la tcnica no puede predecir qu variables son importantes ni permite
explicar el mecanismo involucrado en el proceso fsico, si no es con ayuda de las pruebas
experimentales.
MODELO NUMRICO
En este trabajo utilizamos una versin modificada del programa settlingFoam (Brennan,
2001), contenido dentro de la librera OpenFOAM R (http://www.openfoam.org). Dicho
software, de cdigo abierto permite resolver las ecuaciones mencionadas en el captulo
anterior mediante el mtodo de volmenes finitos, incluyendo modelos turbulentos para
simular las condiciones de transporte en este tipo de flujos. El cdigo utiliza el conocido
lazo tipo PISO para la resolucin del acople entre continuidad y cantidad de movimiento.
Las modificaciones aplicadas al cdigo settlingFoam implicaron la incorporacin de la
ecuacin de conservacin de energa, y un cambio sustancial en el modo en que se tratan
las viscosidades. Se ha llevado a cabo adems, un cambio respecto del tratamiento de la
restriccin de las concentraciones, discutido debajo, que permiti mantener la variable
concentracin de petr- leo dentro de sus lmites naturales 0 y 1 sin introducir
correcciones artificiales que terminaban produciendo prdidas y ganancias espreas de
masa. Las concentraciones k estn limitadas a tomar valores positivos y a cumplir la
restriccin algebraica X 2 k=1 k = 1. (24) Calculando 2 con la Ec.(9), y evaluando
posteriormente 1 con la Ec.(24), se asegura su cumplimiento. Resta sin embargo
restringir las soluciones posibles para 2 al rango [0, 1] para no generar soluciones no
realistas. Para esta aplicacin en particular, es normal adems encontrar situaciones en
donde parte del agua queda atrapada en la matriz continua de petrleo, y puede
considerarse inseparable a los fines prcticos de la simulacin. Debido a esto, la
concentracin 2 no alcanza nunca el valor 1, quedando limitada por debajo de una
concentracin mxima, 2max en nuestra notacin, cuyo valor debe ser determinado
para cada caso particular. Esto se puede modelar, como en nuestro caso, haciendo que
para esta concentracin o cualquiera superior se anule la velocidad de deslizamiento v2j
(ver Eq.(14)), obligando de esta manera a ambas fases a moverse en forma conjunta
impidiendo su separacin. En la implementacin original de Brennan (Brennan, 2001) la
restriccin sobre los valores de concentracin se impone a posteriori en cada iteracin
simplemente truncando por debajo y por Mecnica Computacional Vol XXXII, pgs. 17611783 (2013)
arriba el campo 2 para mantenerlo dentro del rango vlido. Esta solucin, que
aparentemente funcion para los casos mostrados por Brennan, se revel como no
adecuada para su aplicacin al rea petrleo, produciendo efectos espreos de creacin y
destruccin de masa de petrleo, sobre todo en las zonas de alta concentracin de ste.
Otros autores han propuesto soluciones matemticamente elegantes y elaboradas de este
problema notando que en general la ecuacin de transporte de la fase secundaria es una
ecuacin hiperblica con un flujo no convexo, lo que da lugar a ondas simples y
compuestas. El tipo de ondas depende de la ley elegida para la velocidad de
deslizamiento, y de las condiciones iniciales. Los esquemas elegidos debieran entonces
ser capaces de capturar choques, rarefacciones, y ondas compuestas. Se necesitan
esquemas de interpolacin de alto orden para mantener una precisin razonable sin
perder la capacidad de producir soluciones acotadas. Segn Marquez (Mrquez Damin,
2013), dos tcnicas posibles son los resolvedores Riemann-free y los esquemas TVD,
ambos con Flux Corrected Transport (FCT) limitados a un cierto mximo para evitar los
valores espreos de concentracin. Entre las muchas contribuciones presentes en la
literatura una de las ms importantes es el esquema KT de Kurganov y Tadmor
citepKurganov2000, capaz de usar los flujos en las caras obtenidos en el lazo PISO. Este
flujo es la parte convectiva del flujo total utilizado en la ecuacin de transporte de la fase
secundaria. Esta caracterstica lo hace adecuado para ser usado en un esquema de
volmenes finitos centrado en las celdas. Extensiones multidimensionales de esquemas
tipo KT se pueden definir escribiendo estos como esquemas MUSCL (Monotone
Upstream-centered Scheme for Conservation Laws) (Van Leer, 1979). Como ejemplo
reciente, un resolvedor Riemann-Free para esquemas de volmenes finitos centrados en
la celda ha sido propuesto por Marquez (Mrquez Damin, 2013), incluyendo una
extension multidimensional con flujos centrados basada en el trabajo de Greenshields et
al. (Greenshields, 2010). En este trabajo elegimos otro camino para obtener soluciones
acotadas y realistas, sin utilizar esquemas de alto orden. Se elimin el truncamiento
utilizado por Brennan y se modific la forma en que se discretiza por volmenes de control
la ecuacin de transporte, Ec.(9). La idea en que se basa la solucin es la siguiente.
Consideremos dos volmenes de control adyacentes, que comparten una dada cara. El
flujo de petrleo a travs de la cara tiene dos contribuciones, una convectiva y otra
difusiva, estando esta ltima representada en el modelo de drift-flux por la velocidad de
deslizamiento ya mencionada, la cul debe ser evaluada en esa cara. En el trabajo de
Brennan, las velocidades de deslizamiento se calculan en base a las concentraciones en
los volmenes de control adyacentes por esa cara, y se interpola linealmente esta
velocidad en la cara. A diferencia de esto, la propuesta presentada por Barcel y
colaboradores (Barcel et al., 2010) es utilizar la concentracin de petrleo del volumen
de control que se encuentra aguas abajo (downwind) respecto de esa cara para calcular
la velocidad de deslizamiento que se usar para evaluar los flujos de petrleo entre
volmenes de control vecinos.
De esta manera, por ejemplo, si el volumen de control aguas abajo ya tiene una
concentracin 2max la velocidad de deslizamiento utilizada para calcular el flujo
difusivo tendr valor nulo, impidiendo la llegada de ms petrleo a este volumen de
control por este mecanismo.
Emulsiones
Definiciones
1. [Fluidos de perforacin]
Una dispersin de un lquido inmiscible en otro mediante el uso de una sustancia qumica
que reduce la tensin interfacial entre los dos lquidos para lograr estabilidad. Se utilizan
dos tipos de emulsin como lodos: (1) emulsin de aceite en agua (o directa), conocida
como "lodo de emulsin" y (2) emulsin de agua en aceite (o inversa), conocida como
"lodo de emulsin inversa". La primera se clasifica como lodo a base de agua y la
segunda como lodo a base de aceite.
Tipo de dao en el cual hay una combinacin de dos o ms fluidos inmiscibles, entre ellos
gas, que no se separarn en componentes individuales. Se pueden formar emulsiones
cuando se mezclan los filtrados de fluidos o los fluidos inyectados con los fluidos del
yacimiento (por ejemplo, petrleo o salmuera), o cuando cambia el pH del fluido
producido, tal como ocurre despus de un tratamiento de acidificacin. La acidificacin
podra cambiar el pH de 6 o 7 a menos de 4. Las emulsiones se encuentran normalmente
en empaques de grava y disparos o perforaciones, o dentro de la formacin. La mayora
de las emulsiones se rompen fcilmente cuando se retira la fuente de energa que las
mantiene mezcladas. No obstante, algunos agentes estabilizadores naturales y artificiales,
tales como tensioactivos y pequeas partculas slidas, mantienen los fluidos
emulsificados. Los tensioactivos naturales, creados por bacterias o durante el proceso de
generacin de petrleo, pueden encontrarse en muchas aguas y petrleos crudos,
mientras que los tensioactivos artificiales son parte de muchos fluidos de perforacin,
terminacin o estimulacin. Entre los slidos ms comunes que estabilizan las emulsiones
se encuentran el sulfuro de hierro, la parafina, la arena, los sedimentos, la arcilla, el
asfalto, las incrustaciones y los productos de la corrosin. Las emulsiones se tratan,
generalmente, usando solventes mutuos.
3. [Petrleo pesado]
Dispersin de gotas de un lquido en otro lquido con el cual no es completamente
miscible. Las emulsiones pueden formarse en petrleos pesados que contengan una
cantidad significativa de asfaltenos. Los asfaltenos actan como tensioactivos con el agua
de tratamiento o de la formacin. Las gotas de la emulsin resultante tienen enlaces de
alta energa que crean una dispersin compacta de gotas que no se separan fcilmente.
Estas fuerzas que actan en la superficie pueden crear emulsiones de petrleo en agua o
de agua en petrleo. Dichas emulsiones requieren tratamientos de temperatura y
qumicos en el equipamiento de superficie para separarlas.
Desemulsificantes