Abstract
Microtomografía de rayos X se utiliza para la imagen, a una resolución de 6,6 micras, la disposición de poros a escala de los ganglios dióxido de carbono residual en el espacio de poros de una roca de carbonato a presiones y temperaturas representante de formaciones típicos utilizados para el almacenamiento de CO2. El equilibrio químico entre las fases de CO 2, salmuera y roca se mantuvo usando un reactor de alta temperatura a alta presión, replicando condiciones lejos del sitio de la inyección. El flujo de fluido se controla mediante bombas de jeringa de alta temperatura de alta presión. Para mantener representativos condiciones in situ dentro del escáner micro-CT se utilizó una alta presión de fibra de carbono micro-CT coreholder. Difusiva intercambio CO 2 a través de la manga de confinar desde el espacio de poros de la roca para confinar el fluido fue impedido por que rodea el núcleo con un triple envoltura de papel de aluminio. Contraste salmuera reconstruida fue modelado usando una fuente de rayos x policromático, y la composición de salmuera wcomo elegido para maximizar el contraste de tres fases entre los dos fluidos y la roca. Líneas de flujo flexibles se utilizan para reducir las fuerzas de la muestra durante la adquisición de la imagen, que puede causar el movimiento indeseado de la muestra, una de las principales deficiencias en las técnicas anteriores. Un termopar interno, se coloca directamente adyacente al núcleo de roca, junto con una envoltura flexible de la calefacción externa y un controlador PID se utiliza para mantener una temperatura constante dentro de la celda de flujo. Las cantidades importantes de CO 2 fueron atrapados, con una saturación residual de 0,203 ± 0,013, y el tamaño de los ganglios de mayor volumen obedecen distribución de la ley de energía, en consonancia con la teoría de la percolación.
Introduction
Captura y almacenamiento de carbono es el proceso donde se captura el CO2 de grandes fuentes puntuales y se almacena en roca porosa, desplazando salmueras residentes por lo que permanece en el subsuelo durante cientos de miles de años 1. El CO 2 reside en el subsuelo como una densa fase supercrítica (CO2 supercrítico), con propiedades radicalmente diferentes al CO 2 en condiciones ambientales. Hay cuatro mecanismos principales por el cual SCCO 2 puede ser inmovilizado en el subsuelo: estratigráfica, solubilidad, mineral y la captura residual. Atrapamiento estratigráfica es donde se lleva a cabo de CO 2 por debajo sello rocas impermeables; atrapamiento solubilidad es donde CO2 se disuelve en la salmuera residente rodea el CO 2 inyectado 2-4; atrapamiento mineral es donde fases minerales de carbonato se precipitan en la roca 5; y la captura residual o capilar es donde CO 2 está en manos de las fuerzas de superficiegotas como diminutos (ganglios) en el poro-espacio de la roca 6. Esto puede ocurrir ya sea naturalmente, por la migración de la CO 2 penacho de 7-9, o puede ser inducida por la inyección de persecución salmueras 10. Con el fin de comprender los procesos que rigen el flujo y la captura de este CO2 en el subsuelo de una nueva serie de experimentos debe llevarse a cabo, el aprovechamiento de los nuevos avances en la tecnología para entender mejor la física fundamental asociada con flujo multifásico.
Microtomography de rayos X se ha desarrollado como una técnica en los últimos 25 años a partir de los primeros intentos de visualizar ambas muestras geológicas secas 11 y múltiples fases fluidas 12 para el método principal para la imagen no invasiva de núcleos de roca, tanto a efectos de modelado y para experimentación aplicación 13-15. Debido microtomography es no invasivo, que tiene la capacidad para estudiar los sistemas en condiciones representativas, que es particularmente atractive para el sistema -brine-rock CO 2, como el comportamiento del flujo multifase de CO2 supercrítico es altamente dependiente de las propiedades termofísicas, tales como la tensión interfacial y el ángulo de contacto, que a su vez son una fuerte función de las condiciones del sistema, tales como la temperatura, presión y 16-18 salinidad. En un sistema tan complejo, con un conjunto tan extenso y mal entendido de variables interdependientes, los experimentos que utilizan estructuras de poro idealizadas 19 o fluidos analógicas 20,21 no pueden ser aplicables a procesos de flujo en el subsuelo. Imaging múltiples fluidos en condiciones representativas de una formación de inyección prospectivo CO 2 ha, sin embargo, sigue siendo un reto 22. En este estudio nos planteamos una metodología para el estudio de la conducta de varios fluidos a condiciones de yacimiento, se centra en el examen de capilar atrapando 23,24. Esto incluirá el diseño de una estrategia de formación de imágenes, el montaje de la célula de líquidos, la inyección strategy y procesamiento de imagen posterior.
El examen experimental del comportamiento de flujo multifásico poro escala en sistemas de rocas reales se centra en la obtención de imágenes de núcleos de roca parcialmente saturados después de que tanto la inyección no humectante fase (drenaje) y la inyección fase humectante (imbibición). Estos líquidos se inyectan mediante la conexión de los núcleos de las bombas de inyección de fluidos utilizando líneas de flujo flexibles, mientras que confinar el núcleo usando un diseño coreholder tipo Hassler 25. Para imagen se utilizó con éxito la disposición in situ de CO2 supercrítico y salmuera, una novela y la configuración experimental altamente sensible, centrándose principalmente en el uso de un microscopio de alta resolución de rayos x 23,24,26. Los requisitos para la realización de experimentos a temperaturas y presiones elevadas son muy estrictas, y requieren de la evolución reciente de los dos servicios de tecnología de materiales y micro-CT. Los principales requisitos que se deben cumplir son que todo poseedor de núcleo / muestra necesita ser able para soportar altas temperaturas (HPHT) condiciones de alta presión sin dejar de ser lo suficientemente radiografía transparente para permitir la formación de imágenes eficaz. Instrumentos basados en el laboratorio imponen una restricción adicional, como el núcleo titular debe ser lo suficientemente pequeño de tal manera que la fuente de rayos X se puede colocar cerca de la muestra y que suficientemente grande ampliación de rayos x geométrica se puede lograr de manera que el espacio poroso es resuelto con eficacia. Aunque esta restricción se ha relajado un poco con la introducción de la óptica secundaria en las nuevas máquinas de micro-CT basado en laboratorio, no se ha eliminado por completo, especialmente si se desean tiempos de adquisición rápida, como aumentos ópticos más altos tienden a aumentar el tiempo necesario para adquirir imágenes.
Los experimentos con fluidos solubles proporcionan un desafío adicional al utilizar largos tiempos de adquisición, como CO 2 se difundirá a través de las porciones poliméricas de la Asamblea experimental, la reducción de la saturación de fluidos in-situ. Lall estos problemas hicieron que los tiempos de exploración de más de alrededor de 2 horas eran poco prácticos. Con el fin de mantener los tiempos de exploración por debajo de este requisito, especialmente estrictos para fuentes de proyección de laboratorio, el núcleo titular debe estar alrededor de 1 cm de diámetro. Un tamaño de coreholder más grande habría requerido el detector a ser mucho más lejos de la fuente para lograr el mismo aumento geométrico, reduciendo el flujo incidente de rayos x en el detector y por lo tanto el aumento de los tiempos de exposición requeridos de la proyección. La celda de flujo utilizado en estos experimentos se basó en un diseño de la célula Hassler tradicional, construido alrededor de un manguito de fibra de carbono, con un diseño de la manga similar a la utilizada por Iglauer et al 27, pero con dos alteraciones significativas:. 1) el compuesto de fibras de carbono utilizado en la fabricación de la manga se cambió a partir de fibras T700, con una rigidez de 230 GPa, a las fibras de M55, con una rigidez de 550 GPa. Esto no sólo reduce la cantidad de movimiento de la muestra durante la adquisición de la tomografía, pero también aumentó el máximo wojo presión de la célula de 20 MPa a 50 MPa. 2) El manguito se ha alargado de 212 mm a 262 mm para permitir que la fuente y el detector para estar tan cerca de la muestra como sea posible.
Una deficiencia experimental importante en el primer estudio que utiliza micro-CT para examinar CO 2 a condiciones de reservorio fue el uso de líneas de metal para controlar el flujo hacia y desde el núcleo porta-27. A medida que la muestra se hace girar con relación a las bombas, las líneas de flujo también necesitan ser rotados. Líneas de flujo Stiff pueden causar la muestra se mueva, lo que reduce la resolución efectiva de la imagen o hacer algo de o todo el conjunto de datos inutilizable. Para evitar que esto sustituimos todas las líneas de flujo cerca del escenario rotación con flexibles polieteretercetona tubo (PEEK). Estas líneas de flujo eran flexibles, proporcionando muy pequeñas fuerzas laterales (carga) al núcleo titular durante la adquisición. También nos apegamos las líneas de flujo para las válvulas conectadas a la etapa de la muestra, en lugar de fijar las líneas de flujo a la coreholder. Esto significaba que cualquier carga de línea de flujo existente se transmite directamente a la etapa, en lugar de a la muestra, lo que reduce la probabilidad de movimiento de la muestra. Una desventaja importante de usar el tubo de PEEK era que el CO 2 era capaz de difundirse lentamente a través de ella, a través de una escala de tiempo de alrededor de 24 hr. Esto significaba que el CO 2 salmuera saturada que queda en las líneas de flujo sería desaturar gradualmente.
Otra deficiencia experimental importante de estudios anteriores era el control inexacta de la temperatura. Esto puede afectar los resultados en un número de maneras. En primer lugar, la temperatura es un fuerte control sobre la tensión y ángulo de contacto interfacial 16-18. Por otra parte, la solubilidad de la roca tanto SCCO 2 y carbonato en salmuera es también dependiente 28 altamente temperatura. Control de la solubilidad es crítico, como cuando SCCO 2 se inyecta en un acuífero carbonato de solución salina se disolverá en la salmuera residente, formando un ácido carbónico altamente reactivo, que en turn empezar a disolver cualquier calcita presente. Cualquier imprecisión en el control de la solubilidad, por tanto, puede dar lugar a CO2 supercrítico disolución / exsolution o disolución sólida / precipitación.
Estudios previos 27 utilizan un fluido de confinamiento calienta para calentar el coreholder; sin embargo, esto era problemático. Tiene las desventajas asociadas con la dificultad de mantener con precisión una presión de confinamiento constante utilizando un suministro de agua de recirculación, lo que requiere baños de calentamiento adicionales para que el suministro. Además, este sistema sólo mantiene un control preciso de la temperatura en el punto del baño de calentamiento (no en el punto de soporte de núcleo, y el fluido de confinamiento se enfriaría entre el baño de agua y el soporte de núcleo). También requiere tanto una entrada y un puerto de salida para el fluido de confinamiento, aumentando el número de líneas de fluido unidos a la coreholder y así el aumento de carga de línea de flujo.
En lugar de utilizar un fluido de confinamiento climatizada, un cale flexibleschaqueta ng se utilizó para rodear el soporte de núcleo. Este método de calentamiento muy simple resultó en muy poca carga coreholder, y se dejó para el calentamiento preciso y exacto. Se utilizó una película de poliimida de calentamiento extremadamente delgada, con el fin de minimizar el tamaño de la muestra. La construcción de esta película se compone de un elemento grabado de cobre lámina de 0,0127 mm de espesor, encapsulada entre dos capas de 0,0508 mm película de poliimida. Los elementos de cobre presentes en la chaqueta no afectaron notablemente la calidad de imagen. La temperatura se midió usando un termopar que se sienta en el anillo de confinamiento de la célula. Se coloca en el exterior del manguito de confinamiento, lo más cerca posible al núcleo, lo que garantiza una lectura precisa, fiable y estable de la temperatura del fluido intersticial. La película termopar y la calefacción se conecta a un controlador a la medida Proporcional Integral Derivativo (PID), y las temperaturas se controla dentro de ± 1 ° C.
Para mantener un control completo over solubilidad inter-fase, y representan las condiciones actuales en el acuífero lejos del sitio de inyección, antes de la inyección de la salmuera se equilibró con SCCO 2 mezclando vigorosamente los dos fluidos junto con partículas pequeñas (1-2 mm) de la roca huésped en un reactor agitado y se calienta. Todos los componentes de contacto con el fluido dentro de este reactor están hechas de Hastelloy para minimizar la corrosión. El reactor contiene un tubo de inmersión filtrada para permitir el fluido más denso que ser extraído de la base del reactor (salmuera) y el fluido menos denso que ser extraído de la parte superior del reactor (SCCO 2). Bombas de jeringa de alta presión se utilizan para mantener el flujo de control de presión y en el espacio de poros de la roca y en el reactor, con una precisión de desplazamiento de 25,4 nl. El aparato experimental utilizado en este estudio se muestra en la Figura 1. La sal iónica utilizado para el experimento del que se extrajeron los resultados representativos fue yoduro de potasio (KI), ya que tiene un alto peso atómico y por lo que unacoeficiente de atenuación de rayos X de alta, por lo que es un agente de contraste eficaz. , Se requerirían menos atenuantes (sales tales como NaCl) o mezclas podrían ser utilizados sin embargo salinidades mayores para lograr la misma atenuación de rayos x.
Protocol
1. Imaging Estrategia Diseño
- Con el fin de predecir el rendimiento de imagen de diferentes opciones de soluto para la salmuera, calcular el espectro de rayos X de los rayos X incidentes 29-31. Incluir el impacto de la central titular, el conjunto del núcleo y los fluidos de confinamiento en el espectro de rayos x. Un espectro de rayos x ejemplo incidente utilizando un voltaje de aceleración de 80 kV y corriente de electrones de 87 mu se muestra en la Figura 1.
- Comparar este espectro a los factores de transmisión de la muestra que contiene diferentes formadores de poros fluidos. Simular los cambios en el factor de transmisión debido a cambios en el fluido intersticial usando la ley de Beer-Lambert, suponiendo una longitud efectiva óptica de las especies dentro de la muestra, y los coeficientes de atenuación de rayos x calculados (Figura 3) 32. Determinar el factor de transmisión global integrando sobre todas las energías de rayos X incidentes. Un ejemplo de los factores de transmisión efectivas resultantes para el matri rocax y el material de poros espacio, y los cambios en los factores de transmisión en relación con el caso cuando el poro-espacio está lleno de vacío se pueden ver en la Tabla 1.
- Elegir un soluto salmuera y concentración tal que el cambio en el factor de transmisión asociada con salmuera es aproximadamente la mitad del cambio en el factor de transmisión asociada con el sólido. Esto maximizará contraste trifásica en la imagen reconstruida. Pesar la cantidad necesaria de sal (7% (w / w) KI se utilizó para el experimento de la cual se tomaron los resultados representativos) y mezclar enérgicamente con agua desionizada.
- Alternativamente, si se desea una composición de salmuera específica, cambiar el voltaje de aceleración fuente de rayos x para cambiar el espectro de los rayos X incidentes.
2. Montaje de Equipos y Celular
- Montar el equipo como se muestra a continuación en la Figura 2. Uso PEEK líneas de flujo para reducir la carga de la muestra lateral en la celda de flujo.Pon a prueba cada conexión con cuidado que no haya fugas de fluidos.
- Coloque la salmuera, de la composición decidió durante los pasos 1.1 a 1.3, en la base del reactor. Envolver el calentador flexible alrededor de la celda de flujo.
- Construya accesorios terminales de metal. Retire el hilo del "final de un 1/16" 1/8 a 1/8 "accesorio reductor. A continuación, cortar pequeñas ranuras en la cara del "fin de la conexión 1/8 a distribuir el CO 2 inyectado sobre toda la cara del núcleo.
- Pasar el termopar de alta presión a través de las partes de extremo de metal de la celda de flujo y micro-sello usando ¼ "casquillos y la tuerca, por lo que la unión caliente del termopar se encuentra adyacente a la cara de entrada del núcleo, dentro del anillo de confinamiento de la célula.
- Perforar deseado de la muestra en un núcleo de 6,5 mm de diámetro y de 30 mm a 50 mm de longitud. Triturar hasta los confines de la plana central, para asegurar una buena conexión con las piezas de los extremos de metal. Envuelva este núcleo en papel de aluminio y coloque withina fluoro-polímero de la manga de elastómero.
- Conecte los extremos del manguito de elastómero para los accesorios de los extremos de metal. Añadir otra envoltura de papel de aluminio para el exterior de la manga de elastómero antes de colocar la junta caliente del termopar al lado del anillo de confinamiento de la célula y la adición de una envoltura final del papel de aluminio. Esto forma el conjunto del núcleo (Figura 4).
- Montar la célula de micro-flujo con el conjunto del núcleo sellada dentro de ella y conectar la célula para la etapa dentro del recinto micro-CT (Figura 5) mediante una abrazadera montada en la parte superior de la etapa de CT de rotación.
3. sistema de presurización
- Cierre todas las válvulas aparte de válvula 1, 2 y 3, según se define en la figura 3. Carga de CO 2 del cilindro en la bomba 1 y el reactor, a continuación, cierre la válvula 1. Levante lentamente la temperatura y la presión dentro del reactor a la deseada para el poro de fluido durante el experimento.
- Abrir la válvula 14 y cargar el fluido de confinamiento en la bomba 3. Cerrar la válvula de las válvulas 14. Abra 12 y 13. presurizar el anillo de confinamiento de la célula a al menos 10% mayor que la presión del fluido intersticial propuesto.
- Abra la válvula 11. Carga salmuera en la bomba 2. Cierre la válvula 11 y las válvulas abiertas 9, 8 y 6.
- Lentamente presurizar el espacio de poros de la roca hasta que está a la presión de poros de fluido deseado, llenando el espacio de poros de la muestra con salmuera que no ha sido equilibrada con CO2 supercrítico.
- 4. Abrir la válvula Flush más de 1.000 volúmenes de poro de salmuera equilibrada a través del núcleo de la recarga de la bomba 2 a un caudal constante. El volumen de poros se encuentra multiplicando el volumen del núcleo por la porosidad encontrado usando porosimetría de helio.
NOTA: Esto miscibly desplazar la salmuera un-equilibrada, garantizando 100% de saturación salmuera inicial y CREAcondiciones de corte en el núcleo similar a las condiciones del subsuelo en un acuífero en un punto ligeramente por delante de la parte delantera de un penacho de CO2 supercrítico.
4. Flujo de fluidos y de adquisición de imágenes
- Pasar a través de 10 volúmenes de poros (alrededor de 1 ml) de CO2 supercrítico a través del núcleo a velocidades de flujo muy bajas (1,67 × 10 -9 m 3 / seg), asegurando un bajo número capilar de alrededor de 10 -6. Continuamente tomar proyecciones 2D con el fin de medir con precisión el volumen inyectado total de observando el punto cuando SCCO 2 salmuera desplaza en el espacio de los poros.
- Pasar a través de 10 volúmenes de poros (alrededor de 1 ml) de salmuera equilibrada a través del núcleo a la misma velocidad de flujo baja, causando SCCO 2 a quedar atrapado como una fase residual en el poro-espacio.
- Después de los pasos 4,1 o 4,2, tomar las exploraciones de la muestra al drenaje imagen o imbibición, respectivamente. Utilice un tamaño de vóxel de tal manera que todo el diámetro del núcleo encaja dentro del campo de la view.
- Reconstruir las exploraciones mediante un programa de reconstrucción tomográfica. Para escanear toda la longitud del núcleo, mientras que retener un pequeño tamaño de vóxel, reconstruir volúmenes compuestos uniendo múltiples secciones superpuestas, adquiridos secuencialmente.
NOTA: Cada sección requiere alrededor de 400 proyecciones, tomando 15-20 minutos para adquirir, por lo que la exploración de un volumen completo compuesto tomó alrededor de 90 minutos.
5. Procesamiento de Imágenes y Segmentación
- Aplicar un borde medios no locales preservar filtro de 33,34 para el conjunto de datos y corregir las imágenes para cualquier endurecimiento del haz o suavizando los artefactos creados durante la reconstrucción de imágenes mediante el modelado de estos artefactos funciones gaussianas como radialmente simétricas 35.
- Segmento de los datos (gire la información de escala de grises en una representación binaria de la CO 2 dentro de la imagen) por el uso de un algoritmo de cuenca con una semilla generada a partir de un histograma 2D 36, el tratamiento de la CO
- Analizar esta imagen segmentada para encontrar tanto el número total de CO 2 voxels y también los tamaños de cada grupo conectado de CO 2 residual.
Representative Results
Los resultados de un único carbonato, Ketton piedra caliza, un oolita desde la parte superior de la piedra caliza de Lincolnshire miembro se analizaron en 3D con el fin de identificar y medir el volumen de cada ganglio único desconectado, que luego fue marcada (Figura 6). Todo el proceso se llevó a cabo dentro de los programas Avizo Fuego 8.0 y ImageJ 37.
Las imágenes parcialmente saturados segmentadas fueron analizados contando el número de voxels de residualmente atrapado CO2 supercrítico para encontrar la proporción del volumen ocupado por roca atrapado SCCO 2 - la capacidad de atrapar capilar. Esto entonces se puede convertir a una saturación residual (S r) dividiendo este valor por la porosidad tal como se obtiene usando porosimetría de helio. Un porcentaje importante de CO2 supercrítico quedaron atrapados como la saturación residual, con una saturación residual de 0,203 ± 0,013. Esto concuerda con los resultados encontrados en estudios previos utilizando micro-CT 23. Estudios núcleo mayor escala de captura residual en este tipo de roca mostraron una saturación residual inferior de 0,137 ± 0,012 38.
La entrada de salmuera en un núcleo saturado SCCO 2 es un proceso de imbibición, donde un fluido humectante (salmuera) invade cada poro, desplazando fluido no humectante (SCCO 2). En una roca fuertemente mojado esperamos que el agua llene áreas del espacio poroso en orden de tamaño 39,40, atrapando ganglios conexión durante el proceso llamado snap-off. Este proceso debe ser percolación como 41 lo pueden hacer predicciones acerca de la distribución del tamaño de los grupos aislados. El número n de grupos de volumen s (medida en voxels) debe escala como, donde τ es el exponente Fisher 42. Modelado de la red ha demostrado que en celosías regulares cúbicos tridimensionales el valor de este exponente es de alrededor de τ = 2,189 43. Una forma natural de la extracción de este exponente de datos realeses trazar la cantidad agrupada, como se define por Dias y Wilkinson 41.
que debe escalar a medida:
Esto entonces se representa en un gráfico log-log como una función de s (Figura 7), que muestra el comportamiento de ley de potencia para grandes ganglios, pero una escasa representación de la menor ganglios en comparación con el modelo de ley de potencia. El exponente se calcula mediante la exclusión de los ganglios menor que 10 5 voxels (aproximadamente el inicio del comportamiento de ley de potencia) y la realización de Levenberg-Marquardt regresión 44,45 residual usando un algoritmo de ajuste robusto menos absoluta 46,47. Esto se realizó usando un paquete de software comercial. El exponente Fisher para este sistema fue 2,287 ± 0,009, close al valor teórico de 2,189, lo que indica que la imbibición en este sistema es de hecho percolación similares. Más en general estos resultados confirman las conclusiones de los experimentos más grandes núcleo inundaciones 38,48,49 que SCCO2 actúa como la fase no humectante en carbonatos.
Figura 1. aparato experimental, mostrando las bombas, válvulas y reactores utilizados para controlar el flujo y el asentamiento de la coreholder dentro del recinto de micro-CT. Haga clic aquí para ver una versión más grande de esta figura.
Figura 2. El espectro de energía normalizada para los rayos X incidentes sobre el núcleo, fi se filtra a través de la coreholder, confinando manga y confinar fluido. Calculado utilizando SpekCALC 29-31.
Figura 3. Los coeficientes de atenuación lineal de diferentes fluidos y materiales de roca como una función de la energía fotónica.
Figura 4. Detalle del conjunto del núcleo, que muestra una envoltura de aluminio de triple alrededor del núcleo, evitando difusivo intercambio de CO 2 a través de la manga flouro-elastómero. Haga clic aquí para ver una versión más grande de esta figura.
0 / 52440fig5.jpg "/>
Figura 5. Detalle de la, aparatos de calefacción celda de flujo y la ubicación de la asamblea central dentro de la celda de flujo. El termopar debe colocarse lo más cerca posible a la cara de entrada del núcleo. Por favor, haga clic aquí para ver una versión más grande de esta figura .
Figura 6. Imagen del carbonato después del drenaje y la imbibición. (A) Una representación 3D del núcleo después del drenaje donde cada grupo no humectante fase se da un color diferente. (BF) Una representación 3D del núcleo después de cinco experimentos de imbibición, de color como se describe para (A). La amplia gama de colores indica una fase residual mal conectado. (G) Una sección transversal del núcleo después del drenaje. La fase más oscura es el CO2 supercrítico, la fase intermedia es la salmuera y la fase más ligero es el grano de la roca. (H) Una sección transversal del núcleo después de la imbibición. Por favor, haga clic aquí para ver una versión más grande de esta figura.
Figura 7. La distribución del tamaño de los ganglios residual se muestra en la Figura 6.
| Porosidad Material de relleno | Vacío | CO 2 | H 2 O | H2O - 7% en peso de NaCl | H2O - 7% en peso KI | Solid (CaCO 3) |
| Factor de transmisión | 0.25 | 0,247 | 0.243 | 0,242 | 0,224 | 0,202 |
| Cambio en el factor de transmisión en relación con la aspiradora | N / A | -0.003 | -0.007 | -0.008 | -0.026 | -0.048 |
Tabla 1. Resumen de los resultados (factor de transmisión y el cambio en el factor de transmisión en relación con el caso lleno de vacío) a partir de la simulación de las propiedades ópticas de rayos X del material de roca y de poros que llena el espacio fotografiado durante este estudio. Cada columna representa un material diferente llenando el espacio de poros de la roca en el coreholder.
Discussion
Los pasos más críticos para la formación de imágenes con éxito de fluidos multifásicos a presiones y temperaturas elevadas son: 1) El éxito en el aislamiento del fluido de poro del fluido confinar circundante; 2) el equilibrado eficaz de los fluidos y roca antes de la inyección; 3) control de la temperatura efectiva durante todo el experimento; y 4) la segmentación efectiva de las imágenes resultantes.
El uso de las envolturas de aluminio es crítico para el éxito en el aislamiento de la del fluido intersticial desde el fluido de confinamiento que rodea como en su ausencia intercambio difusivo a través de la manga es rápida, y la saturación dentro del núcleo no permanece constante durante la duración de la exploración. Este problema también puede ser evidente cuando el líquido permanece en las líneas de flujo de PEEK durante largos períodos de tiempo (> 2 h) antes de la inyección en el núcleo en el paso 4.1 y 4.2. Una vez más, CO 2 intercambios de modo difuso en todo el plástico, lo que provoca la salmuera para desaturar. Si esto desalmuera saturada se inyecta en el núcleo, la saturación en el núcleo disminuirá a medida que las agrupaciones residuales se disuelven por la salmuera inyectada.
Otros métodos para la equilibración de fluidos y rocas, incluyendo la recirculación de fluido 50, se han propuesto en la literatura. Estos métodos aumentan la complejidad de la configuración experimental, que a su vez han aumentado la cantidad de tiempo para cada experimento, lo que tendría, a su vez aumenta la probabilidad de que la salmuera en líneas de flujo habría desaturado difusiva.
El control eficaz de temperatura es esencial, y la presencia de un termopar dentro del anillo de confinamiento de la celda de flujo es crítico para este. Temperatura sólo se mide en un solo punto, lo que significa que puede haber algo de gradiente a través de la muestra, lo que conduce a un desequilibrio solubilidad y disolución o exsolution. Esto puede ser minimizado mediante la localización de la unión caliente del termopar lo más cerca posible a tél entrada de la cara del núcleo de roca.
La segmentación efectiva de las imágenes resultantes puede ser un verdadero reto con estos sistemas, como la segmentación de imágenes que contienen una saturación parcial de múltiples fluidos es mucho más difícil que la segmentación de imágenes en seco, por lo que el uso de un simple umbral universal de la escala de grises es insuficiente 51. El uso de la segmentación de cuencas no sólo da los resultados más fiables, en comparación con otros algoritmos en la literatura, pero también es la más efectiva para tratar con timbre y volumen parcial 35 artefactos.
Una de las limitaciones más importantes de esta técnica es que sólo se puede acceder al espacio de macro-poros de una roca. La microporosidad (en escalas más pequeñas que la resolución de la imagen) permanece inaccesible, y puede ser importante para el flujo multifase. Las resoluciones más altas revelan una mayor proporción de estas partes del espacio de los poros, sino que también corresponden a una disminución de la ésimoe campo de visión. La aplicabilidad de la técnica a un tipo de roca específico podría abordarse mediante la comparación de la resolución de la digitalización de la distribución del tamaño de garganta de poro obtenido usando el procedimiento independiente, como la presión capilar por inyección de mercurio.
Este método es una técnica líder para la formación de imágenes de poros escala de múltiples fluidos en condiciones de reservorio en los sistemas realistas, con las aplicaciones existentes, incluyendo un estudio comparativo transversal del capilar de captura 24 y la medición de ángulo de contacto 26, y el método es fácilmente aplicable a una amplia gama de sistemas porosos. El trabajo futuro podría estudiar, en la escala de poro, una amplia variedad de una sola fase y flujo multifásico en medios porosos problemas en condiciones representativas de los acuíferos subterráneos, los campos de petróleo y gas y otros sistemas geológicos profundos.
Acknowledgments
Agradecemos la financiación del Centro de Investigación de Carbonatos y Almacenamiento de Carbono Qatar (QCCSRC), proporcionado conjuntamente por Qatar Petroleum, Shell y Qatar Parque Científico y Tecnológico. También reconocemos la financiación del Consorcio Colegio Imperial en Pore-Scale Modelling.
Materials
| Name | Company | Catalog Number | Comments |
| High Pressure Syringe Pump | Teledyne ISCO | 1000D | |
| Parr Reactor | Parr Instrument Company | 4547A - hastelloy | |
| PEEK Tubing | Kinesis | 1560xL | |
| Potassium Iodide Salt | Sigma Aldrich | 30315-1KG | |
| Carbon Dioxide | BOC | CO2 - size E | |
| Thermocouple | Omega Engineering | KMTSS-IM300U-150 | |
| Kapton Flexible Heater | Omega Engineering | KH-112/10-P | |
| X-Ray Microscope | Zeiss | Versa XRM 500 | |
| Snoop Leak Detector | Swagelok | MS-SNOOP-8OZ | |
| Flouro-Elastomer Polymer (Viton) Sleeve | Fisher Scientific | 11572583 | |
| Micro-CT Coreholder | Airborne Composites | 262mm Coreholder | Constructed in conjunction with Imperial College |
| Tomographic program | Zeiss | XM-Reconstructor | |
| ImageJ - image processing | NIH | ImageJ | |
| Matlab | Mathworks | Matlab | Used for regression analysis |
| Avizo | FEI | Avizo |
References
- Special Report on Carbon Dioxide Capture and Storage. , Cambridge, U.K. (2005).
- Ennis-King, J., Paterson, L. Engineering aspects of geological sequestration of carbon dioxide. SPE 77809, Proceedings of the Asia Pacific Oil and Gas Conference and Exhibition. , (2002).
- Weir, G. J., White, S. P., Kissling, W. M. Reservoir storage and containment of greenhouse gases. Transport in Porous Media. 23 (1), 61-82 (1996).
- Lindeberg, E., Wessel-Berg, D. Vertical convection in an aquifer column under a gas cap of CO2. Energy Conversion & Management. 38, 229-234 (1997).
- Lin, H., Fujii, T., Takisawa, R., Takahashi, T., Hashida, T. Experimental evaluation of interactions in supercritical CO2/water/rock minerals system under geologic CO2 sequestration conditions. Journal of Materials Science. 43 (7), 2307-2315 (2007).
- Juanes, R., Spiteri, E. J., Orr, F. M. Jr, Blunt, M. J. Impact of relative permeability hysteresis on geological CO2 storage. Water Resources Research. 42 (12), (2006).
- Hesse, M., Orr, F. M. Jr, Tchelepi, H. A. Gravity Currents with Residual trapping. Journal of Fluid Mechanics. 661, 35-60 (2008).
- MacMinn, C., Szulczewski, M. L., Juanes, R. CO2 migration in saline aquifers. Part 1. Capillary trapping under slope and groundwater flow. Journal of Fluid Mechanics. 662, 329-351 (2010).
- MacMinn, C., Szulczewski, M. L., Juanes, R. CO2 migration in saline aquifers. Part 2. Capillary and solubility trapping. Journal of Fluid Mechanics. 688, 321-351 (2011).
- Qi, R., LaForce, T. C., Blunt, M. J. Design of carbon dioxide storage in aquifers. International Journal of Greenhouse Gas Control. 3 (2), 195-205 (2009).
- Flannery, B. P., Deckman, H. W., Roberge, W. G., D'Amico, K. L. Three-Dimensional X-Ray Microtomography. Science. 237 (4821), 1439-1444 (1987).
- Jatsi, J. K., Jesion, G., Feldkamp, L. Microscopic Imaging of Porous Media with X-Ray Computer Tomography. SPE 20495. SPE Formation Evaluation. , 189-193 (1993).
- Blunt, M. J., et al. Pore-scale imaging and modelling. Advances in Water Resources. 51, 197-216 (2013).
- Cnudde, V., Boone, M. N. High-resolution X-ray computed tomography in geosciences: A review of the current technology and applications. Earth-Science Reviews. 123, 1-17 (2013).
- Wildenschild, D., Sheppard, A. P. X-ray imaging and analysis techniques for quantifying pore-scale structure and processes in subsurface porous medium systems. Advances in Water Resources. 51, 217-246 (2013).
- Espinoza, D. N., Santamarina, J. C. Water-CO2-mineral systems: Interfacial tension, contact angle, and diffusion-Implications to CO2 geological storage. Water Resources Research. 46 (7), (2010).
- Li, X., Boek, E., Maitland, G., Trusler, J. P. M. Interfacial Tension of (Brines + CO2): (0.864 NaCl+0.136 KCl) at Temperatures between (298 and 448) K, Pressures between (2 and 50) MPa, and Total Molalities of (1 to 5) mol/kg. Journal of Chemical and Engineering Data. 57 (4), 1078-1088 (2012).
- Spiteri, E. J., Juanes, R., Blunt, M. J., Orr, F. M. Jr A New Model of Trapping and Relative Permeability Hysteresis for All Wettability Characteristics. SPE Journal. 13 (3), 277-288 (2008).
- Chaudhary, K., et al. Pore-scale trapping of supercritical CO2 and the role of grain wettability and shape. Geophysical Research Letters. 40, 1-5 (2013).
- Karpyn, Z. T., Piri, M., Singh, G. Experimental investigation of trapped oil clusters in a water-wet bead pack using X-ray microtomography. Water Resources Research. 46 (4), (2010).
- Wildenschild, D., Armstrong, R. T., Herring, A. L., Young, I. M., Carey, J. W. Exploring capillary trapping efficiency as a function of interfacial tension, viscosity, and flow rate. Energy Procedia. 4, 4945-4952 (2011).
- Silin, D., Tomutsa, L., Benson, S. M., Patzek, T. W. Microtomography and Pore-Scale Modeling of Two-Phase Fluid Distribution. Transport in Porous Media. 86, 495-515 (2011).
- Andrew, M. G., Bijeljic, B., Blunt, M. J. Pore-scale imaging of geological carbon dioxide storage under in situ conditions. Geophysical Research Letters. 40 (15), 3915-3918 (2013).
- Andrew, M. G., Bijeljic, B., Blunt, M. J. Pore-scale imaging of trapped supercritical carbon dioxide in sandstones and carbonates. International Journal of Greenhouse Gas Control. 22, 1-14 (2014).
- Method and Apparatus for Permeability Measurements. US Patent. Hassler, G. L. , 2,345,935 (1944).
- Andrew, M. G., Bijeljic, B., Blunt, M. J. Pore-scale contact angle measurements at reservoir conditions using X-Ray microtomography. Advances in Water Resources. 68, 24-31 (2014).
- Iglauer, S., Paluszny, A., Pentland, C. H., Blunt, M. J. Residual CO2 imaged with x-ray micro-tomography. Geophysical Research Letters. 38, (2011).
- Carroll, J. J., Slupsky, J. D., Mather, A. E. The Solubility of Carbon Dioxide in Water at Low Pressure. Journal of Physical and Chemical Reference Data. 20 (6), (1991).
- Poludniowski, G. Calculation of x-ray spectra emerging from an x-ray tube. Part II. X-ray production and filtration in x-ray targets. Medical Physics. 34 (2175), (2007).
- Poludniowski, G., Evans, P. M. Calculation of x-ray spectra emerging from an x-ray tube. Part I. Electron penetration characteristics in x-ray targets. Medical Physics. 34 (2161), (2007).
- Poludniowski, G., Landry, G., DeBlois, F., Evans, P. M., Verhaegen, F. SpekCalc: a program to calculate photon spectra from tungsten anode x-ray tubes. Physics in Medicine and Biology. 57 (19), (2009).
- Berger, M. J., et al. XCOM: Photon Cross Section Database (version 1.5). National Institute of Standards and Technology. , Gaithersburg, MD. (2007).
- Buades, A., Coll, B., Morel, J. -M. A non-local algorithm for image denoising. Computer Vision and Pattern Recognition. , (2005).
- Buades, A., Coll, B., Morel, J. -M. Nonlocal Image and Movie Denoising. International Journal of Computer Vision. 76 (2), 123-139 (2008).
- Schluter, S., Sheppard, A. P., Brown, K., Wildenschild, D. Image processing of multiphase images obtained via X-ray microtomography: A review. Water Resources Research. 50 (4), 3615-3639 (2014).
- Jones, A. C., et al. Assessment of bone ingrowth into porous biomaterials using MICRO-CT. Imaging Techniques for Biomaterials Characterization. 28 (15), 2491-2504 (2007).
- Schneider, C. A., Rasband, W. S., Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nature Methods. 9, 671-675 (2012).
- El-Maghraby, R. Measurements of CO2 trapping in Carbonate and Sandstone Rocks. , Imperial College, London. (2013).
- Roof, J. G. Snap-Off of Oil Droplets in Water-Wet Pores. SPE Journal. 10 (1), 85-90 (1970).
- Lenormand, R., Zarcone, C., Sarr, A. Mechanisms of the displacement of one fluid by another in a network of capillary ducts. Journal of Fluid Mechanics. 135, 337-353 (1983).
- Dias, M. M., Wilkinson, D. Percolation with trapping. Journal of Physics A: Mathematical and General. 19, 3131-3146 (1986).
- Fisher, M. E. The theory of equilibrium critical phenomena. Reports on Progress in Physics. 30 (2), 615 (1967).
- Lorenz, C. D., Ziff, R. M. Precise determination of the bond percolation thresholds and finite-size scaling corrections for the sc, fcc and bcc lattices. Physical Review E. 57 (1), 230-236 (1998).
- Levenberg, K. A method for the solution of certain non-linear problems in least squares. Quarterly Journal of Applied Mathmatics. 2, 164-168 (1944).
- Marquardt, D. W. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. Journal of the Society for Industrial and Applied Mathematics. 11 (2), 431-441 (1963).
- Holland, P. W., Welsch, R. E. Robust regression using iteratively reweighted least-squares. Communications in Statistics - Theory and Methods. 6 (9), 813-827 (1977).
- Huber, P. J. Robust Statistics. , John Wiley & Sons, Inc. New York, NY. (1981).
- Akbarabadi, M., Piri, M. Relative permeability hysteresis and capillary trapping characteristics of supercritical CO2/brine systems: An experimental study at reservoir conditions. Advances in Water Resources. 52, 190-206 (2013).
- Pentland, C. H., El-Maghraby, R., Iglauer, S., Blunt, M. J. Measurements of the capillary trapping of super-critical carbon dioxide in Berea sandstone. Geophysical Research Letters. 38, 4 (2011).
- El-Maghraby, R. M., Blunt, M. J. Residual CO2 Trapping in Indiana Limestone. Environmental Science and Technology. 47 (1), 227-233 (2013).
- Sheppard, A. P., Sok, R. M., Averdunk, H. Techniques for image enhancement and segmentation of tomographic images of porous materials. Physica A. 339 (1-2), 145-151 (2004).
