Summary

Imagerie Reservoir État Pore échelle des phases fluides multiples Utilisation de microtomographie à rayons X

Published: February 25, 2015
doi:

Summary

We present a methodology for the imaging of multiple fluid phases at reservoir conditions by the use of x-ray microtomography. We show some representative results of capillary trapping in a carbonate rock sample.

Abstract

Microtomographie à rayons X a été utilisé pour l'image, à une résolution de 6,6 um, l'arrangement pores échelle de résiduelle ganglions de dioxyde de carbone dans l'espace des pores d'une roche carbonatée à des pressions et des températures représentant de formations typiques utilisés pour stockage du CO 2. Équilibre chimique entre le CO 2, de la saumure et de roche phases a été maintenu en utilisant un réacteur à haute température haute pression, reproduisant les conditions loin du site d'injection. L'écoulement de fluide a été contrôlé à l'aide de pompes haute pression de la seringue à haute température. Pour maintenir représentatives des conditions in situ dans le scanner micro-CT une fibre de carbone haute pression micro-CT coreholder a été utilisé. Échange diffusif CO 2 à travers le manchon de confinement de l'espace des pores de la roche au fluide de confinement a été empêchée par entourant le noyau avec une triple enveloppe de feuille d'aluminium. Reconstruit contraste de saumure a été modélisé en utilisant une source de rayons X polychromatique, et la composition de la saumure wcomme choisi pour maximiser le contraste trois de phase entre les deux fluides et de la roche. Conduites d'écoulement flexibles ont été utilisés pour réduire les forces de l'échantillon lors de l'acquisition de l'image, ce qui pourrait causer le mouvement non désirée d'échantillon, une lacune majeure dans les techniques précédentes. Un thermocouple interne, placé directement adjacent au noyau de roche, associée à une enveloppe souple de chauffage externe et un contrôleur PID est utilisé pour maintenir une température constante à l'intérieur de la cellule d'écoulement. Des quantités importantes de CO 2 ont été pris au piège, avec une saturation résiduelle de 0,203 ± 0,013, et la taille des volumes plus grands ganglions obéissent distributions de droit de la puissance, conformément à la théorie de la percolation.

Introduction

Captage et stockage du carbone est le processus où le CO 2 est capté en provenance des grandes sources ponctuelles et stocké dans des roches poreuses, déplaçant saumures résidents de sorte qu'il reste dans le sous-sol pendant des centaines de milliers d'années 1. Le CO 2 réside dans le sous-sol sous forme de phase dense supercritique (SCCO 2), avec des propriétés radicalement différentes de CO 2 dans les conditions ambiantes. Il existe quatre principaux mécanismes par lesquels SCCO 2 pourrait être immobilisé dans le sous-sol: stratigraphique, la solubilité, minérale et le piégeage résiduel. Piégeage stratigraphique est où le CO 2 est maintenu en dessous de roches joint imperméable; solubilité piégeage est où le CO 2 se dissout dans la saumure résident entourant le CO 2 injecté 2-4; piégeage minéral est là phases minérales carbonatées sont précipités dans la roche 5; et le piégeage résiduel ou capillaire est où le CO 2 est maintenu par les forces de surfacegouttelettes de minuscules (des ganglions) dans l'espace des pores de la roche 6. Cela peut se produire soit naturellement, par la migration du panache de CO 2 7-9, ou peut être induite par l'injection de chasse saumures 10. Afin de comprendre les processus qui régissent l'écoulement et de piégeage de ce CO 2 dans le sous-sol d'une nouvelle série d'expériences doit être menée, en exploitant de nouvelles avancées dans la technologie pour mieux comprendre la physique fondamentale associée à flux multi-phase.

Microtomographie à rayons X a développé une technique au cours des 25 ans à partir de premières tentatives de visualiser les deux échantillons géologiques sec 11 et plusieurs phases fluides 12 à la méthode principale pour l'imagerie non-invasive des carottes de roche, à la fois à des fins de modélisation et expérimentale mise en œuvre 13-15. Parce que la microtomographie est non-invasive, il a la capacité d'étudier les systèmes dans des conditions représentatives, ce qui est particulièrement attractive pour le système CO 2 -brine-rock, que le comportement de l'écoulement polyphasique du SCCO 2 dépend fortement sur ​​les propriétés thermo-physiques, telles que la tension interfaciale et l'angle de contact, qui sont à leur tour une forte fonction des conditions de système tels que la température, la pression et la salinité de 16 à 18. Dans un tel système complexe, avec un tel ensemble vaste et mal compris de variables interdépendantes, des expériences utilisant des structures de pores idéalisées 19 ou fluides analogiques 20,21 peuvent ne pas être applicable à l'écoulement des processus dans la subsurface. Imagerie multiples fluides à des conditions représentatives d'un CO 2 formation d'injection prospective a, cependant, est resté un défi 22. Dans cette étude, nous présentons une méthodologie pour l'examen du comportement multi-fluide aux conditions du réservoir, en se concentrant sur ​​l'examen des capillaire piégeage 23,24. Cela comprendra la conception d'une stratégie d'imagerie, l'assemblage de la cellule fluide, l'injection stratégie et traitement de l'image suivante.

L'examen expérimentale du comportement d'écoulement multiphase dimension des pores dans les systèmes réels de roche se concentre sur la formation d'image de carottes de roche partiellement saturés après injection à la fois la phase non mouillante (drainage) et injection de mouillage de phase (imbibition). Ces fluides sont injectés en connectant les noyaux à des pompes d'injection de fluides utilisant des conduites d'écoulement flexibles, tout en confinant le noyau en utilisant une conception de coreholder Hassler de type 25. Pour l'image avec succès l'arrangement in situ de deux BCSC et de la saumure, un roman et dispositif expérimental très sensible a été utilisé, en se concentrant principalement sur ​​l'utilisation d'une haute résolution x-ray microscope 23,24,26. Les exigences relatives à la réalisation d'expériences à des températures et des pressions élevées sont très strictes, et nécessitent récents développements dans la technologie des matériaux et de micro-CT installations. Les exigences clés qui doivent être remplies sont que tout détenteur noyau / échantillon doit être able pour résister à haute pression haute température (HPHT) conditions tout en restant suffisamment x ray-transparente pour permettre l'imagerie efficace. Instruments fondés sur le Lab-imposent une contrainte supplémentaire, selon le porte-noyau doit être assez petit de sorte que la source de rayons X peut être placé à proximité de l'échantillon et que suffisamment grand géométrique grossissement x-ray peut être réalisé de telle sorte que l'espace des pores est efficacement résolu. Bien que cette contrainte a été quelque peu assouplies avec l'introduction de l'optique secondaires dans les nouvelles machines micro-CT effectuées en laboratoire, il n'a pas été complètement enlevé, surtout si les temps d'acquisition rapides sont souhaitées, comme les grossissements optiques plus élevés ont tendance à augmenter le temps nécessaire à l'acquisition images.

Des expériences avec des fluides solubles constituent un défi supplémentaire lors de l'utilisation de longs temps d'acquisition, sous forme de CO 2 va diffuser à travers les parties polymères du montage expérimental, réduire la saturation in situ fluide. Unll ces questions signifie que temps de balayage de plus de l'ordre de 2 h étaient inapplicables. Afin de réduire les temps d'analyse ci-dessous cette exigence, particulièrement sévères pour les sources à base de laboratoire, le porte-noyau doit être d'environ 1 cm de diamètre. Une taille de coreholder plus grand il aurait fallu le détecteur à être beaucoup plus loin de la source pour obtenir le même grossissement géométrique, ce qui réduit le flux incident de rayons X sur le détecteur et par conséquent de plus en plus nécessaire de projection des temps d'exposition. La cellule d'écoulement utilisée dans ces expériences a été basé sur une conception de cellule Hassler traditionnelle, construite autour d'un manchon en fibre de carbone, avec une conception de manchon similaire à celui utilisé par Iglauer et al 27, mais avec deux modifications importantes:. 1) le composite de fibres de carbone utilisé dans la fabrication du manchon a été modifié à partir de fibres T700, avec une rigidité de 230 GPa, de fibres M55, avec une rigidité de 550 GPa. Cela réduit non seulement la quantité de mouvement de l'échantillon pendant l'acquisition de tomographie, mais a également augmenté le maximum working pression de la cellule à partir de 20 MPa à 50 MPa. 2) Le manchon a été allongé de 212 mm à 262 mm pour permettre la source et le détecteur pour être aussi proche que possible de l'échantillon.

Une lacune expérimentale majeur dans la première étude à utiliser micro-CT d'examiner CO 2 dans des conditions de réservoir a été l'utilisation de lignes métalliques pour contrôler le flux vers et depuis le porte-noyau 27. Lorsque l'échantillon est mis en rotation par rapport à la pompe, les lignes d'écoulement doivent également être mis en rotation. Lignes d'écoulement rigides peuvent provoquer l'échantillon de se déplacer, ce qui réduit la résolution de l'image effective ou faisant partie ou la totalité de l'ensemble de données inutilisable. Pour éviter cela, nous avons remplacé toutes les lignes de flux à proximité de la phase de rotation avec souple éther cétone de polyéther (PEEK) tubulure. Ces lignes de flux étaient flexible, en offrant de très petites forces latérales (charge) à la porte-noyau lors de l'acquisition. Nous avons également joint les lignes de flux aux vannes attachés à l'étape de l'échantillon, plutôt que de fixer les lignes de flux à l'coreholder. Cela signifie que ne importe quelle charge de ligne de flux existant a été transmise directement à l'étape, au lieu de l'échantillon, ce qui réduit la probabilité de mouvement de l'échantillon. Un inconvénient majeur de l'utilisation du tube de PEEK est que le CO 2 est capable de diffuser lentement à travers elle, sur une échelle de temps de l'ordre de 24 h. Cela signifie que le CO 2 de la saumure saturée à gauche dans les lignes d'écoulement se désaturer progressivement.

Un autre inconvénient majeur expérimentale des études antérieures était inexact contrôle de la température. Cela peut influer sur les résultats d'un certain nombre de façons. Tout d'abord, la température est un fort contrôle sur la tension et l'angle de contact interfaciale 16-18. En outre, la solubilité des deux SCCO 2 et de roches carbonatées en saumure est également très dépendante de la température 28. contrôle de la solubilité est critique, comme lorsque SCCO 2 est injecté dans un aquifère salin de carbonate il se dissoudra dans la saumure résident, formant un acide carbonique très réactif, ce qui, à turn commencer à dissoudre toute calcite présente. Toute inexactitude dans le contrôle de la solubilité peut donc conduire à SCCO 2 dissolution / exsolution ou solide dissolution / précipitation.

Des études antérieures 27 ont utilisé un fluide de confinement chauffée pour chauffer l'coreholder; mais ce ne était problématique. Il a les inconvénients liés à la difficulté de maintenir avec précision une pression de confinement constante à l'aide d'un approvisionnement en eau de recirculation, nécessitant bains de chauffage supplémentaires pour que l'offre. En outre, ce système maintient seulement un contrôle précis de la température au point du bain de chauffage (pas au point du support de noyau, et le fluide de confinement se refroidir entre le bain d'eau et le porte-noyau). Elle nécessite également la fois une entrée et un orifice de sortie pour le fluide de confinement, augmentant le nombre de conduites de fluide fixés au coreholder et la charge de la ligne d'écoulement de façon croissante.

Au lieu d'utiliser un fluide de confinement chauffée, d'un Infrastructures soupleveste ng a été utilisé pour entourer le support de noyau. Cette méthode de chauffage très simple abouti à très peu de charge de coreholder, et a permis pour le chauffage précis et exacts. Un film de polyimide de chauffage extrêmement mince a été utilisée, afin de minimiser la taille de l'échantillon. La construction de ce film est constitué d'un élément de gravure de cuivre de feuille 0,0127 mm d'épaisseur, encapsulée entre deux couches de film de polyimide de 0,0508 mm. Les éléments de cuivre présents dans l'enveloppe ne affectent pas sensiblement la qualité de l'image. La température a été mesurée en utilisant un thermocouple assis dans l'espace annulaire de confinement de la cellule. Il a été positionné sur le côté extérieur du manchon de confinement, aussi près que possible du noyau, en assurant une lecture précise, fiable et stable de la température du fluide des pores. Le film de thermocouple et le chauffage ont été connecté à un contrôleur sur mesure Proportionnelle Intégrale Dérivée (PID), et les températures ont été contrôlées à ± 1 ° C.

Pour maintenir complète vo de commandeer solubilité entre phases, et représentent les conditions présentes dans l'aquifère loin du site d'injection, avant l'injection de la saumure a été équilibrée avec SCCO 2 en mélangeant vigoureusement les deux fluides ensemble avec de petites particules (1 à 2 mm) de la roche d'accueil dans un réacteur agité et chauffé. Tous les composants en contact avec l'intérieur de ce réacteur sont faites de Hastelloy pour minimiser la corrosion. Le réacteur contient un tube plongeur pour permettre filtré pour le fluide dense pour être extraits de la base du réacteur (saumure) et le fluide moins dense pour être extraits de la partie supérieure du réacteur (SCCO 2). Pompes à seringue à haute pression ont été utilisés pour maintenir le flux de pression et de contrôle dans l'espace des pores de la roche et dans le réacteur, avec une précision de déplacement de 25,4 nl. L'appareil expérimental utilisé dans cette étude est illustré sur la figure 1. Le sel ionique utilisé pour l'expérience dont les résultats représentatifs ont été établis était iodure de potassium (KI), car elle a un poids atomique élevé et donc uncoefficient d'atténuation de rayons X haute, ce qui en fait un agent de contraste efficace. Moins d'atténuation (sels tels que NaCl) ou des mélanges peuvent être utilisés, cependant salinités plus grandes seront nécessaires pour atteindre le même atténuation de rayons X.

Protocol

1. Imagerie Stratégie de conception Afin de prédire la performance d'imagerie de choix de soluté différentes pour la saumure, calculer le spectre de rayons X de l'incident radiographies 29-31. Inclure l'impact de la porte-noyau, le noyau et les fluides de confinement sur le spectre x-ray. Un spectre de rayons X incident, par exemple en utilisant une tension d'accélération de 80 kV et le courant d'électrons de 87 uA est représenté sur la Figure 1. Comparer ce spectre pour les facteurs de transmission de l'échantillon contenant différents fluides de pores. Simuler l'évolution du facteur de transmission en raison de changements dans les pores du fluide en utilisant la loi de Beer-Lambert, en supposant une longueur optique effective de l'espèce dans l'échantillon, et calculées des coefficients d'atténuation de rayons X (figure 3) 32. Déterminer le facteur de transmission globale en intégrant sur toutes les énergies incidents x-ray. Un exemple des facteurs de transmission efficaces en résultent pour l'matri rockx et le matériau poreux de l'espace, et les variations des facteurs de transmission par rapport au cas où l'espace de pore est rempli de vide peut être vu dans le tableau 1. Choisir un soluté de la saumure et concentration telle que la variation du facteur de transmission associée avec de la saumure est d'environ la moitié de la variation du facteur de transmission associé avec le solide. Ceci permet d'optimiser le contraste de phase trois dans l'image reconstruite. Peser la quantité requise de sel (7% (p / p) de KI a été utilisé pour l'expérience dont les résultats représentatifs ont été effectuées) et mélanger vigoureusement avec de l'eau désionisée. Alternativement, si une composition de saumure spécifique est souhaitée, changer le x-ray tension d'accélération de la source pour changer le spectre de rayons X incidents. 2. Assemblée de l'Equipement et cellulaire Assemblez le matériel comme illustré ci-dessous à la figure 2. Utilisez PEEK flowlines de réduire charge latérale de l'échantillon sur la cellule d'écoulement.Testez chaque connexion attentivement les éventuelles fuites de fluide. Placez la saumure, a décidé de la composition au cours des étapes 01.01 à 01.03, à la base du réacteur. Envelopper le chauffage souple autour de la cellule d'écoulement. Construire embouts métalliques. Retirez le fil de la "fin d'un 1/16" 1/8 à 1/8 "raccord réducteur. Ensuite, couper les petites rainures dans le visage de la «fin du raccord 1/8 de distribuer le CO 2 injecté sur toute la face de l'âme. Passer le thermocouple haute pression à travers les parties d'extrémité métallique de la cellule de micro-écoulement et le joint à l'aide ¼ "viroles et l'écrou, de sorte que la soudure chaude du thermocouple se trouve adjacente à la face d'entrée du coeur, à l'intérieur de l'anneau de confinement de la cellule. Percer échantillon désiré dans un noyau de 6,5 mm de diamètre et de 30 mm à 50 mm de longueur. Meuler les extrémités de l'appartement de base, afin d'assurer une bonne connexion avec les embouts en métal. Enveloppez ce noyau dans du papier aluminium puis placez withina fluoro-polymère manchon en élastomère. Raccorder les extrémités du manchon en élastomère pour les raccords d'extrémité en métal. Ajouter une pellicule mince d'aluminium sur l'extérieur du manchon en élastomère avant de placer la soudure chaude du thermocouple à côté de l'anneau de confinement de la cellule et l'addition d'une enveloppe finale de feuille d'aluminium. Cela forme l'ensemble de noyau (figure 4). Assembler la cellule de micro-écoulement avec l'ensemble de noyau scellé à l'intérieur et connecter la cellule à l'étape à l'intérieur de l'enceinte micro-CT (Figure 5) à l'aide d'une pince montée sur le dessus de l'étage de rotation CT. 3. Système de pressurisation Fermez toutes les vannes en dehors de la vanne 1, 2 et 3, tel que défini dans la figure 3. Charge CO 2 du cylindre dans une pompe et le réacteur, puis fermer la vanne 1. élever lentement la température et la pression dans le réacteur à celle souhaitée pour le pores fluide pendant l'expérience. <li> Mélanger vigoureusement le réacteur pendant au moins 12 heures pour se assurer toutes les phases sont en équilibre chimique avant l'injection. Ouvrir le robinet 14 et chargez le fluide de confinement dans la pompe 3. Fermez la vanne 14. Ouvrir les vannes 12 et 13. Pressuriser le anneau de confinement de la cellule à au moins 10% supérieure à la pression de pore fluide proposé. 11. Ouvrez la vanne de charge de la saumure dans la pompe 2. Fermer le robinet 11 et vannes ouvertes 9, 8 et 6. Pressuriser lentement l'espace des pores de la roche jusqu'à ce qu'il est à la pression de fluide des pores désirée, le remplissage de l'espace des pores de l'échantillon avec de la saumure qui n'a pas été équilibrée avec SCCO 2. 4. Ouvrir la vanne de rinçage plus de 1000 volumes de pore de solution de sel équilibrée à travers le noyau de remplissage 2 à une pompe à débit constant. Le volume de pore est obtenue en multipliant le volume de base par la porosité trouvé avec porosimétrie hélium. NOTE: Ce sera miscibly déplacer la saumure non équilibrée, garantissant 100% initial saturation en saumure et creaconditions Ting dans le noyau semblable aux conditions souterraines dans un aquifère à un point légèrement en avant du front d'un panache SCCO 2. 4. Fluide Débit et Image Acquisition Passer à travers 10 volumes de pores (environ 1 ml) de SCCO 2 par l'intermédiaire du noyau à des débits très faibles (1,67 x 10 -9 m 3 / s), assurant un faible nombre capillaire de l'ordre de 10 -6. Prendre constamment projections 2D afin de mesurer avec précision le volume total injecté en observant le moment où SCCO 2 déplace la saumure dans l'espace des pores. Traversez 10 volumes de pores (environ 1 ml) de saumure équilibrée à travers le noyau en même faible débit, provoquant SCCO 2 pour devenir piégé comme une phase résiduelle dans l'espace des pores. Après les étapes 4.1 ou 4.2, prenez analyses de l'échantillon au drainage image ou imbibition respectivement. Utiliser une taille de voxel de telle sorte que tout le diamètre de l'âme se inscrit dans le domaine de la view. Reconstruire les analyses à l'aide d'un programme de reconstruction tomographique. Pour balayer toute la longueur du noyau tandis que conserver une faible taille de voxel, reconstruire volumes composites en cousant ensemble plusieurs sections qui se chevauchent, acquises de façon séquentielle. Remarque: chaque section requise autour des saillies 400, en 15 à 20 min d'acquérir, de sorte que le balayage d'un volume de composite entier a pris environ 90 min. 5. Traitement et Segmentation image Appliquez un bord de moyens non-locales préserver filtre 33,34 pour l'ensemble de données et de corriger les images pour tout durcissement du faisceau ou de ramollissement artefacts créés pendant la reconstruction de l'image en la modélisation de ces artefacts fonctions gaussiennes comme symétrie radiale 35. Segmenter les données (tournent les informations en niveaux de gris en une représentation binaire du CO 2 dans l'image) par l'utilisation d'un algorithme de bassin versant avec une graine généré en utilisant un histogramme 2D 36, le traitement de la CO <sub> 2 comme une phase et la saumure et la roche ainsi que l'autre phase. Analyser cette image segmentée pour trouver à la fois le nombre total de CO 2 voxels et aussi la taille de chaque cluster connecté de CO 2 résiduel.

Representative Results

Les résultats pour un seul carbonate, Ketton calcaire, un oolite de la partie supérieure Lincolnshire Limestone membres ont été analysés en 3D afin d'identifier et de mesurer le volume de chaque ganglion uniques déconnecté, qui a été ensuite marqué (figure 6). Tout le traitement a été menée dans les programmes Avizo incendie 8.0 et ImageJ 37. Les images segmentées partiellement saturés ont été analysés par comptage du nombre de voxels de façon résiduelle piégée SCCO 2 pour trouver la proportion du volume de roche occupé par piégé SCCO 2 – la capacité de piégeage capillaire. Cela peut ensuite être converti en une saturation résiduelle (S r) en divisant cette valeur par la porosité obtenue à l'aide d'hélium comme la porosimétrie. Des proportions importantes de SCCO 2 ont été piégés en saturation résiduelle, avec une saturation résiduelle de 0,203 ± 0,013. Ceci est en accord avec les résultats trouvés dans les études précédentes utilisant micro-CT 23. Des études plus importantes noyau échelle de piégeage résiduelle dans ce type de roche ont montré une saturation résiduelle inférieure de 0,137 ± 0,012 38. L'introduction de la saumure dans un noyau saturé SCCO 2 est un processus d'imbibition où un fluide de mouillage (eau glycolée) envahit chaque pore, déplacer un fluide non mouillant (SCCO 2). Dans une roche fortement humide l'eau nous nous attendons à l'eau pour remplir les zones de l'espace des pores par ordre de taille 39,40, le piégeage ganglions déconnecté au processus appelé sécable. Ce processus devrait être percolation comme 41 afin prédictions peuvent être faites au sujet de la distribution de la taille des groupes isolés. Le nombre n de groupes de volumes s (mesurée en voxels) doit se adaptent comme, où τ est l'exposant Fisher 42. Réseau de modélisation a montré que dans des réseaux réguliers cubiques en trois dimensions la valeur de cet exposant est autour de τ = 2,189 43. Une façon naturelle de l'extraction de cet exposant à partir des données réellesest de tracer la quantité binned, tel que défini par Dias et Wilkinson 41. qui devrait évoluer comme suit: Ce est alors tracée sur un terrain log-log en fonction de s (figure 7), montrant un comportement en loi de puissance pour les gros ganglions, mais une sous-représentation des plus petits ganglions par rapport au modèle de loi de puissance. L'exposant a été calculé en excluant les ganglions inférieure à 10 5 voxels (environ le début du comportement en loi de puissance) et la réalisation de Levenberg-Marquardt régression 44,45 résiduelle en utilisant un algorithme d'ajustement moins robuste absolue 46,47. Ceci a été réalisé en utilisant un logiciel commercial. L'exposant Fisher pour ce système était 2,287 ± 0,009, close à la valeur théorique de 2,189, ce qui indique que l'imbibition dans ce système est en effet comme percolation. Plus généralement ces résultats confirment les conclusions dans les grandes expériences core-inondation 38,48,49 que SCCO deux actes que la phase non mouillant dans les carbonates. Figure 1. appareil expérimental, montrant les pompes, vannes et réacteur utilisés pour contrôler le flux et le siège de la coreholder dans l'enceinte micro-CT. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure. La figure 2. Le spectre d'énergie normalisé pour les rayons X incidents sur le noyau, fi filtré à travers le coreholder, confinant manche et confinant fluide. Calculé selon SpekCALC 29-31. Figure 3. Les coefficients d'atténuation linéaires de différents fluides et de matériaux rocheux en fonction de l'énergie des photons. Figure 4. Détail de l'ensemble de noyau, montrant une pellicule d'aluminium triple autour du noyau, empêchant diffusion échange de CO 2 dans le manchon fluoro-élastomère. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure. 0 / 52440fig5.jpg "/> Figure 5. Détail de la cellule de débit, appareils de chauffage et de l'emplacement de l'ensemble de base au sein de la cellule d'écoulement. Le thermocouple doit être placé aussi près que possible de la face d'entrée de la base. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure . Figure 6. Image de carbonate après le drainage et l'imbibition. (A) de rendu en 3D de l'âme après drainage où chaque groupe de phase non mouillante est donnée à une couleur différente. (BF) un rendu 3D de l'âme après cinq expériences d'imbibition, comme décrit pour la couleur (A). La large gamme de couleurs indique une phase résiduelle mal connecté. (G) A-section transversale de l'âme après le drainage. La phase la plus sombre est le SCCO 2, la phase intermédiaire est la saumure et la phase la plus légère est le grain de la roche. (H) Une coupe transversale de l'âme après imbibition. Se il vous plaît cliquer ici pour voir une version plus grande de cette figure. Figure 7. La distribution de la taille des ganglions résiduelle représenté sur la figure 6. Matériau de remplissage Porosité Vide CO 2 H 2 O H 2 O – 7% en poids de NaCl H 2 O – 7% en poids KI Solide (CaCO 3) Facteur de transmission 0,25 0,247 0,243 0,242 0,224 0,202 Changement de facteur de transmission par rapport à l'aspirateur N / A -0,003 -0,007 -0,008 -0,026 -0,048 Tableau 1. Résumé des résultats (facteur et le changement de facteur de transmission par rapport au boîtier rempli de vide transmission) de simulation des propriétés optiques des rayons X de la matière de roche et des pores espace remplissage imagée au cours de cette étude. Chaque colonne représente un matériau différent remplir l'espace des pores de la roche dans le coreholder.

Discussion

Les étapes les plus critiques pour l'imagerie réussie de fluides polyphasiques, à des pressions et des températures élevées sont les suivantes: 1) L'isolement réussi du fluide interstitiel à partir du fluide de confinement environnante; 2) l'équilibrage efficace des fluides et de rock avant l'injection; 3) contrôle de la température effective pendant toute l'expérience; et 4) la segmentation effective des images résultantes.

Utilisation des enveloppes d'aluminium est essentielle pour l'isolement réussi de pore de fluide à partir du fluide de confinement entourant qu'en son absence échange diffusif à travers le manchon est rapide, et la saturation dans le noyau ne est pas constante pendant la durée de l'analyse. Ce problème peut également être évident lorsque le fluide reste dans les conduites de PEEK pour des périodes de temps prolongées (> 2 h) avant l'injection dans le noyau de l'étape 4.1 et 4.2. Une fois encore, échanges de CO 2 diffuse à travers le plastique, causant la saumure pour désaturer. Si cela dede la saumure saturée est injectée dans le noyau, la saturation dans le noyau diminue à mesure que les grappes résiduelles sont dissous par la saumure injectée.

D'autres procédés pour la mise en équilibre des fluides et des roches, dont recirculation de fluide 50, ont été proposées dans la littérature. Ces méthodes augmentent la complexité du montage expérimental, ce qui aurait augmenté la quantité de temps pour chaque expérience, ce qui aurait à son tour, augmente la probabilité que la saumure dans les conduites d'écoulement aurait diffusive désaturé.

Contrôle de la température efficace est essentielle, et la présence d'un thermocouple à l'intérieur de l'anneau de confinement de la cellule d'écoulement est critique pour ce produit. La température est mesurée seulement à un seul point, ce qui signifie qu'il peut y avoir gradient à travers l'échantillon, conduisant à un déséquilibre de la solubilité et la dissolution ou exsolution. Ceci peut être minimisé en plaçant la jonction chaude du thermocouple le plus près possible à tqu'il entrée face de l'âme du rock.

La segmentation efficace des images résultantes peut être un véritable défi avec ces systèmes, comme la segmentation d'images contenant une saturation partielle de plusieurs fluides est beaucoup plus difficile que la segmentation d'images sèches, donc l'utilisation de seuillage simple universelle de gris est 51 insuffisante. L'utilisation de la segmentation des bassins versants non seulement donne les résultats les plus fiables, par rapport à d'autres algorithmes de la littérature, mais est aussi la plus efficace à traiter avec sonnerie et de volume partielle artefacts 35.

L'une des limitations les plus significatifs de cette technique est qu'elle ne peut accéder à l'espace macro-pores d'une roche. La microporosité (sur des échelles plus petites que la résolution de l'image) reste inaccessible, et peut être important pour l'écoulement polyphasique. Les résolutions plus élevées révèlent une plus grande proportion de ces parties de l'espace des pores, mais correspondent aussi à une diminution de ee champ de vue. L'applicabilité de cette technique à un type de roche spécifique pourrait être traitée par la comparaison de la résolution de l'analyse de la distribution de taille des pores de la gorge obtenue par la méthode indépendante tels que la pression capillaire d'injection de mercure.

Cette méthode est une technique de premier plan pour l'imagerie pores échelle de plusieurs fluides aux conditions du réservoir dans les systèmes réalistes, avec des applications existantes, y compris une étude comparative transversale du capillaire piégeage 24 et la mesure de l'angle de contact 26, et la méthode est facilement applicable à une large gamme de systèmes poreux. Les travaux futurs pourraient étudier, à l'échelle des pores, une grande variété de monophasé et écoulement polyphasique dans poreuses problèmes des médias à des conditions représentatives des aquifères souterrains, des champs de pétrole et de gaz et d'autres systèmes géologiques profondes.

Disclosures

The authors have nothing to disclose.

Acknowledgements

Nous tenons à remercier le financement du Centre de recherches carbonates et stockage du carbone Qatar (QCCSRC), fourni conjointement par Qatar Petroleum, Shell, et le parc Qatar Science & Technology. Nous reconnaissons également le financement de l'Imperial College Consortium sur Pore Modélisation à l'échelle.

Materials

Name of Reagent/ Equipment Company Catalog Number Comments/Description
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

  1. . . Special Report on Carbon Dioxide Capture and Storage. , (2005).
  2. 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).
  3. 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).
  4. Lindeberg, E., Wessel-Berg, D. Vertical convection in an aquifer column under a gas cap of CO2. Energy Conversion & Management. 38, 229-234 (1997).
  5. 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).
  6. Juanes, R., Spiteri, E. J., Orr, F. M., Blunt, M. J. Impact of relative permeability hysteresis on geological CO2 storage. Water Resources Research. 42 (12), (2006).
  7. Hesse, M., Orr, F. M., Tchelepi, H. A. Gravity Currents with Residual trapping. Journal of Fluid Mechanics. 661, 35-60 (2008).
  8. 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).
  9. 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).
  10. 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).
  11. Flannery, B. P., Deckman, H. W., Roberge, W. G., D’Amico, K. L. Three-Dimensional X-Ray Microtomography. Science. 237 (4821), 1439-1444 (1987).
  12. 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).
  13. Blunt, M. J., et al. Pore-scale imaging and modelling. Advances in Water Resources. 51, 197-216 (2013).
  14. 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).
  15. 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).
  16. 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).
  17. 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).
  18. Spiteri, E. J., Juanes, R., Blunt, M. J., Orr, F. M. A New Model of Trapping and Relative Permeability Hysteresis for All Wettability Characteristics. SPE Journal. 13 (3), 277-288 (2008).
  19. 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).
  20. 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).
  21. 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).
  22. 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).
  23. 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).
  24. 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).
  25. Hassler, G. L. Method and Apparatus for Permeability Measurements. US Patent. , (1944).
  26. 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).
  27. Iglauer, S., Paluszny, A., Pentland, C. H., Blunt, M. J. Residual CO2 imaged with x-ray micro-tomography. Geophysical Research Letters. 38, (2011).
  28. 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).
  29. 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).
  30. 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).
  31. 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).
  32. Berger, M. J., et al. XCOM: Photon Cross Section Database (version 1.5). National Institute of Standards and Technology. , (2007).
  33. Buades, A., Coll, B., Morel, J. -. M. A non-local algorithm for image denoising. Computer Vision and Pattern Recognition. , (2005).
  34. Buades, A., Coll, B., Morel, J. -. M. Nonlocal Image and Movie Denoising. International Journal of Computer Vision. 76 (2), 123-139 (2008).
  35. 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).
  36. 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).
  37. Schneider, C. A., Rasband, W. S., Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nature Methods. 9, 671-675 (2012).
  38. El-Maghraby, R. . Measurements of CO2 trapping in Carbonate and Sandstone Rocks. , (2013).
  39. Roof, J. G. Snap-Off of Oil Droplets in Water-Wet Pores. SPE Journal. 10 (1), 85-90 (1970).
  40. 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).
  41. Dias, M. M., Wilkinson, D. Percolation with trapping. Journal of Physics A: Mathematical and General. 19, 3131-3146 (1986).
  42. Fisher, M. E. The theory of equilibrium critical phenomena. Reports on Progress in Physics. 30 (2), 615 (1967).
  43. 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).
  44. Levenberg, K. A method for the solution of certain non-linear problems in least squares. Quarterly Journal of Applied Mathmatics. 2, 164-168 (1944).
  45. 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).
  46. Holland, P. W., Welsch, R. E. Robust regression using iteratively reweighted least-squares. Communications in Statistics – Theory and Methods. 6 (9), 813-827 (1977).
  47. Huber, P. J. . Robust Statistics. , (1981).
  48. 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).
  49. 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).
  50. El-Maghraby, R. M., Blunt, M. J. Residual CO2 Trapping in Indiana Limestone. Environmental Science and Technology. 47 (1), 227-233 (2013).
  51. 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).

Play Video

Cite This Article
Andrew, M., Bijeljic, B., Blunt, M. Reservoir Condition Pore-scale Imaging of Multiple Fluid Phases Using X-ray Microtomography. J. Vis. Exp. (96), e52440, doi:10.3791/52440 (2015).

View Video