Determinación Experimental de Permeabilidades Relativas en Dos Fases Utilizando Imágenes de Tomografía y Análisis Por Simulación Numérica.
Determinación Experimental de Permeabilidades Relativas en Dos Fases Utilizando Imágenes de Tomografía y Análisis Por Simulación Numérica.
Determinación Experimental de Permeabilidades Relativas en Dos Fases Utilizando Imágenes de Tomografía y Análisis Por Simulación Numérica.
Director
Codirector (Opcional)
Bucaramanga
2019
OBJETIVO
OBJETIVOS ESPECIFICOS
Ajustar software de simulación existente para que realice el ajuste automático de los
muy representativa, que recree de una manera análoga el comportamiento in-situ de los fluidos
Por consiguiente, el uso de la tomografía computarizada es una gran herramienta debido a que
durante el desarrollo de los experimentos; esto ayuda a tener una mejor estimación de las
Como consecuencia de este trabajo grado se espera dejar una base experimental que sirva de
FLUJO MULTIFASICO
De forma sencilla el flujo multifásico es el flujo simultaneo de diferentes fases (solido, liquido,
gas). En cuanto se refiere a flujo multifásico en el área de petróleo quiere decir, flujo entre
Tipos de fluidos
𝜕𝑉
=0
𝜕𝑝
𝑝 𝑉
𝑑𝑉
−𝑐 ∫ 𝑑𝑝 = ∫
𝑃𝑟𝑒𝑓 𝑉𝑟𝑒𝑓 𝑉
𝑉
𝑒 𝑐(𝑝𝑟𝑒𝑓 −𝑝) =
𝑉𝑟𝑒𝑓
𝑝=presión, psia
1 1 𝜕𝑧
𝑐𝑔 = − ( )
𝑝 𝑧 𝜕𝑝 𝑇
Regímenes de flujo
Flujo de estado estacionario. El flujo de estado estable se reconoce si la presión en cada lugar
del yacimiento permanece constante, lo que significa, no hay cambio con el tiempo.
Esto se presenta cuando el yacimiento está lleno y además es soportado por un fuerte acuífero.
Flujo en estado inestable. También conocido como flujo transitorio, hace referencia a que la
contante.
Flujo en estado pseudo estable- En este régimen de flujo el cambio de la presión con respecto
Ley de Darcy. La Ley de Darcy es esencial para el movimiento de los fluidos en medios
porosos.
La ecuación fue desarrollada en 1856 por Henry Darcy, expresa que la velocidad de un fluido
𝑞 1 𝑑𝑝
𝑣= 𝛼 (− )
𝐴 𝜇 𝑑𝑥
𝑞 𝑘 𝑑𝑝
𝑣= = (− )
𝐴 𝜇 𝑑𝑥
𝑑𝑝
= Gradiente de presión, atm/cm
𝑑𝑥
PERMEABILIDAD RELATIVA
Múltiples estudios han concluido que la permeabilidad efectiva de un fluido es una función de la
La permeabilidad absoluta es una característica del medio poroso que cuantifica la capacidad de
flujo de una roca, cuando hay presencia de un solo fluido. En presencia de dos o mas fluidos que
fluyen al mismo tiempo, la permeabilidad relativa de cada fase es una relación entre la
𝑘𝑒𝑜
𝑘𝑟𝑜 =
𝑘
𝑘𝑒𝑤
𝑘𝑟𝑤 =
𝑘
𝑘𝑒𝑔
𝑘𝑟𝑔 =
𝑘
Donde
𝑘= Permeabilidad absoluta
siguiente forma, primero se satura el núcleo con agua (fase mojante), luego se desplaza el agua
con aceite (fase no mojante), hasta una saturación residual de agua. Este proceso se puede
Proceso de imbibición. Este proceso realizado en laboratorio parte saturando el núcleo con
agua (fase mojante), después se desplaza el agua hasta saturación irreducible, desplazando aceite.
Este proceso es de drenaje que sirve para determinar las saturaciones originales que se
encontraban en yacimiento. La fase mojante (agua) se reinyecta, este proceso es conocido como
imbibición.
En ocasiones los datos de permeabilidad relativa no pueden ser obtenidas de las muestras reales
de reservorio y por tanto se requiere otras alternativas para dicho fin. Se han implantado varios
𝑆𝑜
𝑆𝑜∗ =
1 − 𝑆𝑤𝑐
𝑆𝑤 − 𝑆𝑤𝑐
𝑆𝑤∗ =
1 − 𝑆𝑤𝑐
𝑆𝑔
𝑆𝑔∗ =
1 − 𝑆𝑤𝑐
Donde
Correlación de Wyllie y Gardner. Wyllie y Gardner (1958) observaron una relación lineal
agua 𝑆𝑤∗
Sistema aceite-agua
𝑆𝑤∗
𝑘𝑟𝑤 = (𝑆𝑤∗ )2 − 𝑘𝑟𝑜 ( )
1 − 𝑆𝑤∗
Sistema gas-aceite
𝑆𝑜∗
𝑘𝑟𝑜 = 𝑆𝑜∗ − 𝑘𝑟𝑔 ( )
1 − 𝑆𝑜∗
Correlación de Torcaso y Wyllie. Torcaso y Wyllie (1958) propusieron una expresión para
(𝑆𝑜∗ )4
𝑘𝑟𝑜 = 𝑘𝑟𝑔 ( )
(1 − 𝑆𝑜∗ )2 (1 − (𝑆𝑜∗ )2 )
Correlación de Pirson. Pirson (1958), las ecuaciones son aplicadas para rocas mojadas por
agua.
Imbibición
2
𝑆𝑤 − 𝑆𝑤𝑐
(𝑘𝑟 )𝑛𝑜 𝑚𝑜𝑗𝑎𝑛𝑡𝑒(𝑛𝑤) = (1 − ( ))
1 − 𝑆𝑤𝑐 − 𝑆𝑛𝑤
Drenaje
0.5
(𝑘𝑟 )𝑛𝑜 𝑚𝑜𝑗𝑎𝑛𝑡𝑒(𝑛𝑤) = (1 − 𝑆𝑤∗ )(1 − (𝑆𝑤∗ )0.25 √𝑆𝑤 )
𝑆𝑤 =Saturación de agua
Método de Corey
Sistema gas-petróleo
𝑘𝑟𝑜 = (1 − 𝑆𝑔∗ )4
Sistema agua-petróleo
1 − 𝑆𝑤 4
𝑘𝑟𝑜 = ( )
1 − 𝑆𝑤𝑐
𝑆𝑤 − 𝑆𝑤 4
𝑘𝑟𝑤 =( )
1 − 𝑆𝑤𝑐
Este tipo de expresiones son generalmente utilizadas en simulaciones numéricas. Las expresiones
Sistema aceite-agua
1 − 𝑆𝑤 − 𝑆𝑜𝑟𝑤 𝑛𝑜
𝑘𝑟𝑜 = (𝑘𝑟𝑜 )𝑆𝑤𝑐 = ( )
1 − 𝑆𝑤𝑐 − 𝑆𝑜𝑟𝑤
𝑛𝑤
𝑆𝑤 − 𝑆𝑜𝑟𝑤
𝑘𝑟𝑤 = (𝑘𝑟𝑤 )𝑆𝑜𝑟𝑤 =( )
1 − 𝑆𝑤𝑐 − 𝑆𝑜𝑟𝑤
1 − 𝑆𝑤 − 𝑆𝑜𝑟𝑤 𝑛𝑝
𝑃𝑐𝑤𝑜 = (𝑃𝑐 )𝑆𝑤𝑐 = ( )
1 − 𝑆𝑤𝑐 − 𝑆𝑜𝑟𝑤
𝑛𝑔𝑜
1 − 𝑆𝑔 − 𝑆𝑙𝑐
𝑘𝑟𝑜 = (𝑘𝑟𝑜 )𝑆𝑔𝑐 =( )
1 − 𝑆𝑔 − 𝑆𝑙𝑐
𝑛𝑔
𝑆𝑔 − 𝑆𝑔𝑐
𝑘𝑟𝑔 = (𝑘𝑟𝑔 )𝑆𝑤𝑐 =( )
1 − 𝑆𝑙𝑐 − 𝑆𝑔𝑐
𝑛𝑝𝑔
𝑆𝑔 − 𝑆𝑔𝑐
𝑃𝑐𝑔𝑜 = (𝑃𝑐 )𝑆𝑙𝑐 =( )
1 − 𝑆𝑙𝑐 − 𝑆𝑔𝑐
usualmente mediante el método de mínimos cuadrados con el fin de concordar con datos
TECNICAS DE MEDICIÓN
Estado inestable
Se trata de una prueba de desplazamiento dinámico, y existen tres formas de determinar las
promedio del fluido. Sus cálculos son simples solo basta con aplicar la ley de Darcy de
promedio.
saturación de fluido efluente, para desplazar el fluido a una velocidad de inyección constante. La
Método Toth et al. (2002), se emplea para desplazamiento constante de la tasa de inyección y
presión constante. Permite la estimación directa con base en los datos de desplazamiento. Ofrece
Estado estable
Es una condición de equilibro donde la caída de presión como la tasa de flujo me mantienen
constantes en el espacio del medio poroso y en el tiempo, esto es necesario para que las
saturaciones se conserven constantes. En una prueba de flujo bifásico, los fluidos inyectados
deben alcanzar condiciones estables de manera independiente. Una vez se consiga este equilibro
𝑄𝑤 ∗ 𝜇𝑤 𝐿
𝑘𝑒𝑤 = = 𝑓(𝑆𝑤 )
𝐴 ∗ ∆𝑝𝑤
𝑄𝑜 ∗ 𝜇𝑜 𝐿
𝑘𝑒𝑜 = = 𝑓(𝑆𝑤 )
𝐴 ∗ ∆𝑝𝑜
ESTIMACION DE SATURACIONES A PARTIR DE INFORMACION DE
TOMOGRAFIA
Son varios los autores que se han referido al tema sobre la obtención de saturaciones in-situ a
partir de los experimentos de inyección de fluidos, esto desde hace más de treinta años. Este
proceso consiste, en que la roca por medio de la naturaleza de sus fluidos verá afectada sus poros
y sus coeficientes de atenuación; también es importante saber que la atenuación en los materiales
se relaciona con la densidad, numero atómico, lo cual se puede aprovechar para “dopar” los
fluidos que intervienen en el experimento y mezclando estos con compuesto que tengan un
elevado número atómico, para producir un aumento en la atenuación, distinguiendo que la lectura
de las imágenes sea más practica y permitiendo tener clara la presencia del fluido y cuantificando
partir de diferentes cortes para entender su estructura interna, por esto la TAC aprovecha la
diferencia que existe en la atenuación de los diferentes materiales sobre los rayos- X que lo
atraviesan, esto nos permite que el resultado sean representaciones graficas de los diferentes
coeficientes de atenuación, que se pueden discretizar en unidades mas pequeñas y las cuales
denominamos voxels.
Podemos encontrar una expresión para determinar la saturación en el trabajo de (Wellington &
Vinegar, 1987) . Esta ecuación (1) se usa para sistemas bifásicos y aplica en imágenes voxel a
fenómenos: el efecto fotoeléctrico, que busca usar la densidad y el numero atómico del material y
la energía del fotón, y el efecto Compton que depende de solo de la densidad y la energía (Hunt,
P., Engler, P., and Bajsarowicz C. 1988). Observe que del coeficiente de atenuación de un vóxel,
Para realizar este proceso de estimación de saturaciones nos hemos apoyado del Sistema
este consiste en un tomógrafo multicorte de 64 corte por rotación (GE Optima 660) acoplado
psi).
este consiste en que la salmuera se dopa con yoduro de potasio, permitiendo que aumente su
atenuación y múltiples escaneos se hacen con diferentes intervalos de tiempo de prueba, para así
poder tomar los datos de avance de los fluidos y estimar las saturaciones.
Es posible observar de forma más clara las características físicas y los cambios que sufre la roca
de forma gráfica; las propiedades físicas de la roca usada pueden intervenir y afectar el
experimento, por esto es importante que esa tenga una forma homogénea o no la podríamos usar
Cabe resaltar que este método utilizado no es el único que existe para determinar las saturaciones;
también por medio de las densidades encontradas en la roca se pueden realizar imágenes CT, es
el caso de la universidad de Calgary que hace uso de esta metodología y es una de las alternativas
cuando las condiciones operacionales no nos permiten el uso práctico del primer método
planteado
METODOLOGIA
El equipo TAC debe ser un escáner médico de energía graduable, de manera que se pueda
realizar escaneo a energía alta (superior a 100 kV) y energía baja (Inferior a 100 kV) (Wellington,
presiones determinadas, de tal forma que dé como respuesta imágenes de tomografía durante el
proceso de desplazamiento.
Para mejor contraste visual, normalmente se utilizan dopantes, que simplemente son trazadores
que se agregan a los fluidos los cuales provocan un aumento en la atenuación y por tanto se hace
Los perfiles de saturación serán obtenidos basándose en las imágenes procedentes de tomografía,
escaneo.
Es necesario definir un área de interés para que el software (Image-J) analice la información que
compete a dicha región además de pedir al software un perfil de valor promedio del número CT,
y luego de procesar las imágenes se tiene una escala de colores donde se puede evidenciar el
frente de avance del aceite, con los perfiles obtenidos se procede a obtener las saturaciones.
Ajustar software de simulación existente para que realice el ajuste automático de los
El ajuste se realizará con los datos que se obtuvieron en los experimentos de desplazamiento,
dicho ajuste se llevará a cabo en una herramienta implementada en Matlab, la cual permite
Se utilizará el algoritmo propuesto por Ortiz, Carrillo y Rojas (2017) para determinar las
(pp. 331-456).
283-329).
Akin, S., & Kovscek, A. R. (2003). Computed tomography in petroleum engineering research.
Berg, S., Armstrong, R., Ott, H., Georgiadis, A., Klapp, S. A., Schwing, A., … Enzmann, F.
(2014). Multiphase Flow in Porous Rock Imaged Under Dynamic Flow Conditions with Fast X-
Johnson, E. F., Bossler, D. P., & Naumann, V. O. (s. f.). Calculation of Relative Permeability
London, M., Cameron, S. M., Donald, J., & Wassmuth, F. R. (2014). Waterflooding Experiments
with X-ray CT Imaging. SPE Heavy Oil Conference-Canada. Presentado en SPE Heavy Oil
Moghadasi, L., Guadagnini, A., Inzoli, F., Colapietro, D., Bartosek, M., & Renna, D. (2015).
Mohanty, K. K., & Salter, S. J. (s. f.). Multiphase Flow in Porous Media: II. Pore-Level
experimental para la estimación de permeabilidades relativas en dos y tres fases por medio de
Ortiz Meneses, A. F., Herrera Otero, E., & Santos, N. (2017). Estimación de saturaciones in-situ
Ortiz Meneses, A. F., Plata Chaves, J. M., Herrera Otero, E., & Santos, N. (2017).
Sedaghat, M. H., Gerke, K., Azizmohammadi, S., & Matthai, S. K. (2016). Simulation-based
Sharma, B. C., Brigham, W. E., & Castanier, L. M. (1997). CT imaging techniques for two-phase
DOE/BC/14994--T1, 501560).
Tiab, D., & Donaldson E. C. (2014). Relative Permeability. En Petrophysics (pp. 826-834)
Torabi, F., Mosavat, N., & Zarivnyy, O. (2016). Predicting heavy oil/water relative permeability
Toth, J., Bodi, T., & Szucs, P. (s. f.). Determining Relative Permeability from Unsteady-State
Withjack, E. M., & Devier, C. (s. f.). The Role of X-Ray Computed Tomography in Core
Analysis. 12.
Zhang, D., Li, S., Jiao, S., Shang, Y., & Dong, M. (2019). Relative permeability of three
immiscible fluids in random porous media determined by the lattice Boltzmann method.