Summary
Synchrotron wurde schnell Tomographie dynamisch Bild Auflösung von Kalk in Gegenwart von CO 2 -gesättigten Sole bei Lagerstättenbedingungen verwendet. 100 Scans wurden bei einer Auflösung von 6,1 & mgr; m über einen Zeitraum von 2 h gebracht.
Abstract
Die unterirdische Speicherung Dauerhaftigkeit ist ein wichtiges Anliegen für die Kohlenstoffabscheidung und -speicherung. Pumping CO 2 in Karbonatvorkommen hat das Potenzial , geologische Dichtungen aufzulösen und ermöglichen CO 2 zu entkommen. Allerdings sind die Lösungsprozesse bei Reservoirbedingungen schlecht verstanden. So werden zeitaufgelöste Experimente erforderlich, um die Art und Geschwindigkeit der Auflösung an der Porenskala zu beobachten und vorherzusagen. Synchrotron schnelle Tomographie ist ein Verfahren zur Aufnahme hochauflösender zeitaufgelöste Bilder von komplexen Porenstrukturen wesentlich schneller als herkömmliche μ-CT. Der Diamond Light Source Rosa Strahl wurde dynamisch Bild Auflösung von Kalk in Gegenwart von CO 2 -gesättigten Sole bei Lagerstättenbedingungen verwendet. 100 Scans wurden bei einer Auflösung von 6,1 um über einen Zeitraum von 2 Stunden genommen. Die Bilder wurden segmentiert und die Porosität und Permeabilität wurden unter Verwendung von Bildanalyse und Netzwerk-Extraktion gemessen. Porosität erhöht gleichmäßig entlang der length der Probe; die Steigerungsrate der verlangsamte jedoch sowohl die Porosität und Durchlässigkeit zu späteren Zeitpunkten.
Introduction
Ein wichtiges Anliegen von Kohlenstoff - Abscheidung und Speicherung (CCS) ist eine langfristige Lagerung Sicherheit 1, 2. Kohlendioxid, CO 2, in den Untergrund injiziert wird in die Wirts Sole auflösen und Kohlensäure 3, 4, 5 bilden. Diese saure Sole hat das Potenzial , mit zu reagieren und das umgebende Gestein lösen, insbesondere wenn das Wirtsgestein ist Kalkstein 6. Die Auflösung kann für die weitere Bildung Permeabilität 7 und höhere Speicher Dauerhaftigkeit 8 günstig und lassen sein. Jedoch können geologische Dichtungsintegrität durch diese Auflösung beeinträchtigt werden und erlauben CO 2 an der Oberfläche 9 zu migrieren. Eine genaue prädiktive Modellierung von Speicher Dauerhaftigkeit ist somit abhängig von voll Verständnis Auflösung in der Sole-Felsen-System und die Verteilung unddie Rate der Flüssigkeitsbewegung in den Untergrund 10, 11, 12.
Jedoch ist die Natur und die Auflösungsgeschwindigkeit in Carbonate abhängig sowohl von den Eigenschaften der Sole 13, 14, 15, 16 und dem Wirtsgestein 17. Die Auflösungsraten sind auch stark abhängig von Soletemperatur und Druck 6, die Entwicklung experimenteller Techniken machen für lebenswichtige komplexe zeitabhängige Prozesse an repräsentativen Reservoirbedingungen zu messen.
Frühere Experimente haben beobachtet , dass Feld angelegte Reaktionsgeschwindigkeiten sind in der Regel um Größenordnungen niedriger als experimentelle Chargenreaktor Messungen 18, 19. Verwitterung, Mineral heterogenekeit und unvollständige Vermischung in einem heterogenen Strömungsfeld sind mögliche Erklärungen für dieses Phänomen. Es ist jedoch nicht möglich, die wichtigsten Faktoren ohne direkte Beobachtung des sich entwickelnden Porenraum während der Reaktion zu beurteilen. Somit sind dynamische Poren angelegte Experimente erforderlich, um sowohl die Einblicke in das Zusammenspiel von Transport und Reaktion bereitzustellen und Vorhersagemodelle zu validieren.
Eine etablierte experimentelle Methode zur Untersuchung von Porentechnischen Prozessen in Kohlenstoff - Speicheranwendungen ist die Röntgenmikrotomografie (μCT) 20, 21. μ-CT hat mehrere Vorteile: es wird eine hohe räumliche Auflösungen von bis zu etwa 1 & mgr; m, es ist nicht-invasiv und liefert dreidimensionale Bilder. Kalkstein Auflösung wurde am Kern (~ cm) untersuchten 22 skaliert und es wurde festgestellt , dass Rock-Sole Reaktion physikalische Heterogenität erhöht. Um das Verständnis für wie unterschiedlich transport und Reaktionsbedingungen zu verändern, die komplexen festen und Porenstrukturen ist es erforderlich, reaktionsinduzierte Veränderungen im Porenraum Geometrie, Topologie und fließen in unterirdischen Gesteins Systeme bei Reservoir Temperaturen und Drücken und bei einer höheren Auflösung zu messen, um im Detail zu untersuchen poren Maßstab Prozesse. Dieses Dokument beschreibt ein Verfahren zur reaktiven Lösungsprozesse in Fels mit komplexen Porenstrukturen zu studieren und konzentrieren sich auf die Zeit gemessen und ortsabhängigen Reaktionsgeschwindigkeit zwischen einem CO 2 -acidified Sole und Kalkstein - Felsen in Reservoirbedingungen.
Es wurden mehrere Studien, die in komplexen Carbonaten bei Reaktionsausgesehen haben , 23, 24, 25, 26, 27, aber durch experimentelle oder Abbildungs Constraints sie waren entweder beschränkt auf vor und nach der Reaktion wurden Bilder oder nicht abgeschlossenan repräsentativen Untergrundbedingungen. Menke et al. 28 wurde in situ Abbilden von Reaktion zwischen einem CO 2 -acidified Sole und Ketton Kalkstein an der Porenskala über einen Zeitraum von mehreren Stunden und bei einer Temperatur und Druck Vertreter eines Aquifers in etwa 1 km in Tiefe dynamischen durchgeführt. Jedoch ist eine relativ homogene Ketton Gestein mit großen Körnern, die dem Bild in sehr kurzer Zeit (~ 17 min) und mit wenigen Vorsprüngen einfach ist (~ 400). Die meisten Karbonatgestein haben komplexe Porenstrukturen, die viele Projektionen erfordern genau zu lösen, die eine sehr zeit traditionelle μ-CT mit intensiver Prozess sein kann - entweder mit einem monochromatischen Strahl mit einer Synchrotronquelle oder mit Benchtop-Röntgenscanner. Somit ist eine schnelle Methode der Tomographie erforderlich dynamisch reaktionsinduzierte Veränderungen in heterogenen Carbonate zu sehen.
Die Menge der Zeit es braucht, die dem Bild eine Probe durch den Fluss von th gesteuert wirde Röntgenquelle. Ein Verfahren zum Scannen ist schnell die polychromen Strahl einer Synchrotronquelle 20 zu verwenden. Diese so genannte 'Pink Beam liefert Größenordnungen intensiveres Licht als Benchtop-Quellen und daher können die Bilder auf der Zehner-Sekunden eher Zeitskalen als Stunde genommen werden. Ein Undulator, der von einer periodischen Struktur von Dipolmagneten besteht erzeugt das rosa Balken. Der Elektronenstrahl wird erzwungenen Schwingungen zu unterziehen, wie es die Magnete und als Folge durchläuft strahlt Energie. Die erzeugte Energie wird konzentriert Wellenlängenbändern zu verengen und ist sehr intensiv. Spiegel und Filter werden dann experimentellen Anforderungen zu verengen das Spektrum des Lichts anzupassen verwendet. Spiegel absorbieren die Hochenergiespektrum, während Filter die niederen Energien absorbieren. Es ist daher möglich, das Spektrum auf die gewünschte Strahlungsband nur mit diesen Werkzeugen zu verengen.
Allerdings Fluss mit dieser intensiven Röntgen ist nicht ohne Herausforderungen. Dasgeringere Energieröntgenstrahlen des rosa Strahlspektrum werden durch die Probe als Wärme absorbiert. Dies kann mit der Temperatursteuerung der in situ Vorrichtung stören und bewirken , dass CO 2 aus der Lösung 20 bis exsolve. CO 2 -gesättigten Sole ist sehr empfindlich sowohl gegen Hitze und Druck und damit eine geringe Änderung im thermischen Gleichgewicht signifikant den pH der in situ verändern Flüssigkeit 5. Damit eine sorgfältige Konstruktion und Steuerungselemente für die Röntgenspektrum muß vor der Bildgebung in die Strahllinie Ausrüstung eingebaut werden.
Fast Tomographie erzeugt auch eine große Menge an Daten mit einer hohen Rate. Die Grenzen der Daten von der Kamera und anschließender Lagerung bieten eine erhebliche technische Herausforderung ausgelesen. Einige haben überwinden dies durch mehrere aufeinanderfolgende Scans zu nehmen und sie auf dem Speicher der Kamera zu speichern, bevor sie auf externen Datenserver zu lesen. Dies erfordert jedoch, dass das Experiment sein relatively kurz wie Kameraspeicher kann nur eine begrenzte Menge an Daten zu halten. Binning der Daten auf der Kamera reduziert auch Übertragungszeit, da sie die Datenmenge reduziert, um übertragen zu werden, aber es hat das Potenzial, die Qualität der Bilder zu verringern. Alternativ können die Daten mit der Kamera werden nach jedem Scan vor dem Start des nächsten übertragen ab, die die gesamte Zeit zwischen den Scans zu erhöhen. In dieser Studie wurden die letztere Methode bei jeder Bildaufnahme unter ~ 45 Sekunden und Daten abgelesen zusätzlich ~ 30 s nimmt.
Bei Scans mit einer hohen Rate nehmen, muss die Probenbühne Spin viel schneller als mit herkömmlichen Scan und damit die mögliche Winkel Belastung der Kernhalter ist groß. Kohlefaser, während Röntgen transparent, flexibel, wenn betont. Wenn die Probe bewegt sich während der Bildaufnahme Bild Unschärfen auftreten. Die Haltehülse Kern wurde entwickelt, so kurz wie möglich sein, um diese potentiellen Belastungen zu mildern. Zusätzlich flexible polyether Ether-Keton (PEEK) Schlauch an allen Elementen der Versuchsapparatur der Nähe der Bühne verwendet wurde, so dass die Bühne drehen war frei. Ein Nachteil von PEEK Rohren ist , dass es durchlässig ist 2 auf diffusive Zeitskalen zu CO. Fluid in den Leitungen für lange Zeiträume befinden, werden allmählich entsättigt über einen Zeitraum von etwa 24 h werden. Alle Linien , die nicht in der Nähe der Kernhalter waren , wurden aus rostfreiem Stahl hergestellt und die Flüssigkeit war voräquilibriert in einer kräftig gemischt Hastelloy - Reaktor erhitzt und unter Druck zu experimentellen Bedingungen 23, 29, 30.
Die Versuchsapparatur ist in Abbildung 1 dargestellt. Reservoir Temperatur wird in dem Kernhalter gehalten, indem die Außenseite der Hülse in einer Röntgen transparent Heizband Umwickeln und ein Thermoelement durch die radiale Öffnung der Zelle und in das Sperrmedium eingeführt wird. Ein Proportional-Integral-Derivative (PID) Regler dann Temperatur geregelt innerhalb von 1 ° C. Druck und Strömungsbedingungen wurden beibehalten drei Hochdruckspritzenpumpen verwendet, die mit einer Flussrate von 0.001 ml / min genau sind. Zwei Salze wurden für das Experiment verwendet, ein hoch absorbierendes 25% wt KI unreaktiv Sole und einem niedrigen Absorbieren 1% wt KCl, 5% wt NaCl reaktive Sole. Der Unterschied in der Dämpfung machte es leicht, die Ankunft von reaktiven Sole im Kern machen Totvolumen Berechnungen unnötig zu sehen.
Protocol
1. Imaging-Strategie Entwurf
- Berechne den Röntgenspektren der Strahllinie am höchsten rosa Strahlenergie und Flussmittel, um die Abbildungsleistung unter Verwendung der experimentellen Kurve tuning und Messfilter Transmissionen vorherzusagen. Ein Beispiel für die Röntgenspektren der Diamond Light Source I13-2 'Pink Strahl' ist in Abbildung 2 dargestellt.
- Als untere energetische Röntgenstrahlen verursachen Erwärmung der Probe und ergibt nicht Bildgebungskontrast, auszufiltern, den unteren Teil des Röntgenspektrums derart, dass nur die höchste Energieröntgenstrahlen in Abbilden der Probe verwendet werden. Wähle in Netzfiltern von Materialien , die bei den gewünschten niedrigen Wellenlängen von Licht geeignet zu der Lichtquelle absorbieren , indem die theoretische Filterdurchlässigkeit bei den zur Verfügung stehenden Lichtwellenlängen - Berechnungs 31, 32. Hier verwenden Aluminium und Gold für die Strahllinie an dieser Lichtquelle.
- Verwenden, um ein Bandpassfilter, bestehend aus einem combination von Röntgenfilter als Hochpass und einem Röntgenspiegel nahe dem kritischen Winkel als Tiefpassfilter arbeitet. In diesem Fall verwenden Sie einen Satz von 0,2 mm Pyrokohlenstoff und 0,2 mm Aluminium-Filter und für den Spiegel einer platinbeschichteten Streifen unter einem Einfallswinkel von 1,15 mrad. Der Spiegel reflektiert nur Licht unterhalb von 30 keV und zusätzliche inline Filter haben von 2 mm Al und 0,1 & mgr; m Au, die Spitzen jeweils bei 13 und 22 keV haben Absorptions installiert mehrere der niedrigeren Energieröntgenstrahlen zu filtern. Figur 3 zeigt die Strahllinienabbildungsvorrichtung.
- Wählen Sie einen Szintillator, die reichlich an den Strahlrohren verfügbaren Lichtfrequenzen und Fluss funkelt. In diesem Fall wird die Scintillationsschirm von 250 um dicke Cadmium Wolframat (CdWO 4) hergestellt, die mit 750 & mgr; m-dicken Bleiwolframat (PbWO 4) gestapelt ist. Dann wählen Sie ein Objektiv und Kamera, die ein geeignetes Sichtfeld und Schnappen Zeitauflösung für die Experime habenntal Anforderungen. In diesem Fall Paar eine 1,25X Objektivlinse mit einer Apertur von 0,04 mit einem PCO EDGE 5.5 CMOS-Kamera und verwendet, um zu erfassen ein 4-mm Sichtfeld mit einer Bildrate von 0.001 s.
- Wählen Sie die 'FlyScan' Technik für die Bildaufnahme als diese Methode der Stufe Rotation reduziert Probe Vibration. Traditionelle Erwerb erfordert, dass die Bühne an jedem Winkelinkrement zu stoppen, einen Vorsprung nehmen, und dann zum nächsten Winkel bewegen. Die Bildaufnahme während dieser dynamischen Tomographien wurde mit einem 'FlyScan' erfolgen, die als die Stufe tomographischen Scan dauert bewegt und nimmt einen Winkelschritt, so dass die Differenz zwischen jedem aufeinanderfolgenden Projektions klein ist. Die 'FlyScan' Methode eliminiert die kleinen Schwingungseffekte des Start-und-Stop-Motion und bietet eine höhere Bildqualität schneller.
2. Montage von Anlagen und Zell
- Laden Sie den Kern in die Zelle in der Vorbereitung für Kernflut.
- Zuerst wickeln Sie den Kern in einer Schicht aus Aluminiumfolie und legen Sie in eine Hülse (zB Viton) (Abbildung 4).
- Schneiden Sie die Hülse der Größe, so dass es 2 mm kürzer als die Gesamtlänge des Kerns und der Innen Endstücken. Die Anschlussarten sind 1/16 "nationale Rohrgewinde (NPT) zu Anschlussverschraubungen, die bis 5 mm Außendurchmesser bearbeitet wurden, während die Hülse 4 mm Innendurchmesser.
- Strecken Sie die Hülse über die 5 mm Endstücken eine gute Abdichtung zu schaffen. Stellen Sie sicher, dass es keinen Raum zwischen den Endstücken ist und den Kern, um sicherzustellen, dass der Grenzdruck nicht über die Hülse nicht komprimieren und abschnüren fließen.
- Wickeln Sie die Armaturen und die Hülse in zwei zusätzliche Schichten aus Aluminium sowohl zu verhindern gasförmigem CO 2 diffundiert in das Sperrmedium und zu halten , die Hülse in Bezug auf den Armaturen und verhindern einen hydraulischen Weg von den Eingrenzungs verbinden und Porenflüssigkeiten.
- Setzen Sie den Kern holder wieder zusammen mit Schläuchen und Dichtungen wieder in Position schieben und Abdichten der Endkappen und Endstücken durch die Schrauben zu ersetzen.
- Scannen der Kern in sich überlappenden Teile um 4 mm in der Breite und Länge. Kalibrieren des Scan-Belichtungszeit auf einen durchschnittlichen Zählwert von etwa 15.000, die ohne Über Sättigen der Szintillator ein hohes Signal-Rausch-Verhältnis gewährleistet. Nehmen Sie jede trockene Scan mit mindestens 2.400 Projektionen Phasenkontrast und Kantenschärfe zu halten.
- Nehmen Sie flach und dunkle Bilder des Szintillators so dass Beschädigungen und Außengeräusche können bei der Rekonstruktion berücksichtigt werden. Nehmen Wohnungen durch den Kernhalter aus dem Blickfeld bewegtund Aufnehmen eines Bildes von nur dem Szintillator mit dem Strahl auf. Nehmen Sie darks die gleiche Methode mit dem Strahl aus verwenden.
3. System Druckerhoehung
- Legen Sie eine 1% Gew Kaliumchlorid (KCl) 5% Gew Natriumchlorid (NaCl) Sole in zerlegtem Reaktor durch Flüssigkeit in durch den oberen Teil des Reaktorbehälters zu gießen.
- In pulverisierte Karbonatfelsen die gewünschte Sole Säure zu erreichen. In diesem Fall wurde kein Carbonat zugesetzt.
- Montieren Sie den Reaktor durch Anziehen der Schrauben und mit Heizband Umpackzentrums und den Temperaturfühler in die obere einsetzen.
- Last CO 2 in die Einspritzpumpe durch Öffnen des Ventils 1 (V1 in Figur 1).
- Schließen Sie das Ventil 1 und unter Druck bis 100 bar Einspritzpumpe.
- Öffnen Sie Ventil 2 , den Reaktor mit CO 2 zu überfluten. Erhitzen des Reaktors auf 50 ° C, um eine PID-gesteuerte Heizung Packung in Kombination mit einem Temperaturfühler und con mitnuierlich rühren mit einem Rührer Mitreißen von einem externen Elektromotor angetrieben. Äquilibrieren die Sole mit CO 2 bei 10 MPa und 50 ° C für 2 bis 6 h , um sicherzustellen , daß die Sole vollständig mit CO 2 und das Carbonat gesättigt ist vollständig gelöst.
- Vor dem Kernhalter verbindet, spülen Sie das System vollständig von Luft und mögliche Ausfällungen in den Leitungen aus früheren Experimenten. Dazu verbinden die Linien oberhalb und unterhalb des Kernhalters den Kernhalter (U1 und U2) zu umgehen.
- Last deionisiertes (DI) Wasser in die Aufnahme Pumpe durch das Ventil 11 von der empfangenden Pumpeneinstellung nachzufüllen.
- Öffnen Sie die Ventile 7, 4 und 3 und die Empfangs Pumpe auf konstanten Druckmodus verwenden, um DI-Wasser rückwärts durch das System fahren und aus dem Ventil 3 unterhalb des Reaktors. Verwenden Sie etwa zehn Systemvolumen die Leitungen frei von Luft sind, um sicherzustellen, und sauber gespült.
- Leeren Empfangs Pumpe einnd dann laden eine Sole von 25% Gew KI in die Aufnahme Pumpe durch das Ventil 11 und Last DI-Wasser in die Eingrenzungs Pumpe durch das Ventil 10.
- Schließen Sie das Ventil 10 und geöffneten Ventilen 8 und 6. Verwenden Sie die Begrenzungs Pumpe den Kern auf 2 MPa zu beschränken.
- Schließen Sie das Ventil 11 und unter Druck bis 10 bar die Empfangspumpe.
- Öffnen Sie die Ventile 9, 7, 4, & 3 und verwenden Sie die resultierende Druckabfall KI zu fahren dotierten Sole durch den Kern.
- Schritt, um die einschließende und Porendruck schrittweise bis eine angemessene Strömungsgeschwindigkeit einstellt. Fahren Sie etwa zwei komplette System Volumina Sole durch den Kern und die Flüssigkeit ablaufen lassen durch das Ventil 3 unterhalb des Reaktors. Auf diese Weise wird die gesamte Luft aus dem System gespült, und der Kern mit hohem Kontrast Salzlösung geflutet, dass die Ankunft des undotierten reaktiven Sole leicht zu beobachten ist.
- Ventil 3 schließen und schrittweise die Eingrenzungs und Porendruck zu erhöhen, bis der Kern bei 12 MPa begrenzt ist und der Porendruck beträgt 10 MPa. Schalterauf den PID-Regler um den Kern 50 zu bringen, ° C.
- Stoppen die Empfangspumpe, Schließ-Ventil 3 und das offene Ventil 5 an der Basis des Reaktors das Reaktorsystem mit dem Kern zu verbinden.
4. Strömungs- und Image Acquisition
- Zentrieren Sie die Mitte des Kerns in das Sichtfeld und nehmen Sie 2-D-Projektionen kontinuierlich als der Kern überflutet wird Fortschritt der Kernflut zu verfolgen. Nehmen Sie 2-D-Projektionen durch auf dem rosafarbenen Strahl drehen und mit Hilfe der Kamerabilder zu nehmen, ohne die Bühne dreht. Ab 2-D-Projektionen vor reaktiven Fluidinjektion gibt eine klare vor Sole Bild, das später auf nachfolgende Sole gefüllt Bildern verglichen werden.
- Stellen Sie die Empfangspumpe somit bei der gewünschten Strömungsrate zum Nachfüllen Flüssigkeit aus dem Reaktor durch den Kern an den gewünschten Strömungsbedingungen ziehen , während die Einspritzpumpe verlassen Druck von dem vorderen Ende zu regulieren.
- Überwachen Sie die 2-D-Projektionenfür Änderungen in der Dämpfung, die die Ankunft von reaktiven Sole signalisieren. Wenn reaktive Sole ankommt, wird die Übertragung des Kerns zu erhöhen und die 2-D-Projektionen erhellt erheblich mehr Licht den Szintillator trifft und die dotierte Sole wird durch die stark röntgentransparente reaktiven Fluid verdrängt. Wenn es keine Dämpfungsdifferenz zwischen reaktiven und nichtreaktiven Salzlösung wird dann in Abhängigkeit auf der Strahllinienspektrum, Neustarten des Experiments aus Schritt 2.1 mit einer höheren Salzkonzentration KI oder einen anderen stark absorbierenden Salz verwendet, kann erforderlich sein.
- Stoppen Sie die 2-D-Scans und 3-D nehmen Tomographien nacheinander so schnell wie die Abbildungsvorrichtung ermöglicht. Verwenden Sie ca. 1.000 Projektionen pro Scan. Scannen der Kern nur 180 ° der Drehung unter Verwendung von (im Gegensatz zur traditionellen 360 °). Während weniger Rotationsgrade verringert das Signal-Rausch-Verhältnis verwenden, ist es schneller und Dehnung zu vermeiden hilft und die Strömung und elektrische Leitungen verheddern. Nehmen Sie 3-D-Scans bis either die Frist erreicht wird oder der Kern sieht ausreichend gelöst, dass es eine unmittelbare Gefahr der inneren strukturellen Kollaps ist (und damit sowohl der Grenzdruck und Zukunft gesamten Kern trocken Scandaten den Verlust verursacht).
- Nach dem letzten Scan genommen wird, drucklos , das System effizient weiter Umsetzen des Kerns zu vermeiden.
- Erste Station der Empfangspumpe. Dann schließen Ventil 5 Anschluss des Reaktors mit dem Rest des Systems.
- Schritt aus der Systemdruck nach unten die Eingrenzungs mit und Pumpen zu halten um 1 MPa mehr Druck auf das Sperrmedium zu empfangen.
- Wenn innerhalb von 1 MPa von Atmosphärendruck erreicht ist, öffnen Sie das einschließende und Empfangen von Pumpen Ventilen 10 und 11 und laufen in konstanten Flussmodus alle verbleibenden Flüssigkeit abzulassen.
- Fahren Sie den PID-Regler aus und öffnen Sie die 4-Wege-Union (U2) an der Spitze des Kernhalters alle verbleibenden Systemdruck zu lösen.
- Lösen Sie langsam die einschließende Liniewährend überschüssiges DI Fang Wasser mit saugfähigem Papier beschränken. Schließen Ventile 6 und 7 und trennen union 1 und die elektrischen Leitungen.
- Lösen Sie die Tischklemme und entfernen Sie von der Bühne den Kernhalter.
- die Kernanordnung aus dem Kernhalter vorsichtig entfernen und dann die Ärmel von den inneren Endarmaturen trennen. Sie nicht den Kern von der Hülse zu entfernen, da dies die zerbrechliche Kern reagiert beschädigen. Die Hülse bedeckten Kern in einem Becher voll von DI-Wasser alle potentiell reaktiven Sole zu verdünnen und alle Reaktion zu stoppen.
- Trocknen der gesamte Kern in einem 60 ° C Ofen für mindestens 12 h. Dann wieder montieren Kern auf der Bühne eine traditionelle Probenträger verwenden, und scannen sie wieder mit der gleichen Auflösung und Projektionen wie die anfängliche Trocken Scan.
5. Bildverarbeitung
- Korrigieren Sie die rekonstruierten Bilder für jede Strahlaufhärtung im Zusammenhang mit einem polychromen Strahl mit der Annahme, dass jeder betroffen sind radial symmetric Gauß - Funktionen 33.
- Filtern Sie die Bilder eine randerhaltenden Filters wie nichtlokalisierte Mitteln Signal zu erhöhen , um Rauschen 34, 35 (siehe Zusatz Datei ).
- Segment die trockenen Scanbilder einen Wendepunkt mit Segmentierung 36 - Algorithmus und definieren Sie die Samen als Fels und nichtig (siehe Zusatz Datei ).
- Nehmen Sie das erste Bild des Kerns mit reaktiven Sole und registrieren jedes nachfolgende Bild mit dem ersten Bild und resamplen es die Lanczos 37 Resampling - Methode mit dem ersten Bild als Referenz verwendet. Als ablaufenden Reaktion neigt Kanten zu verwischen, auf die Bilder Watershed Segmentierung ist nichtausreichend für eine genaue Segmentierung.
- Subtrahiert jedes reagierte Kern Bild aus dem ersten Bild, um das Differenzbild zu erhalten. Segment die Differenzbilder in Veränderung und keine Veränderung. Registrieren Sie den segmentierten trocken Scan zum ersten reaktiven Scan und dann subtrahieren aus dem segmentierten trockenen Scan die segmentierte Änderung des segmentierten zu erreichen , umgesetzt Bilder 38.
6. Modellierung
- Verwenden Sie die digitalisierten Bilder als Eingaben in eine entweder eine direkte Navier-Stokes - Strömungslöser 39, 40, oder ein Netzwerk - Extraktion Modell 41 (Abbildung 8) Permeabilitätsveränderungen charakterisieren und physikalische Einblicke in die Dynamik der Auflösung liefern.
Representative Results
Die Reaktion wurde zwischen Calcit und ungepufferte scCO 2 gesättigte Sole in 4 mm Durchmesser 1,2 cm langen Portland Karbonat Kern 42 abgebildet. Portland Karbonat ist ein relativ reines (<99%) Calcit oolite mit einem komplexen , heterogenen Porenstruktur 43. Der geringe Energieröntgenstrahlen wurden durch Hindurchleiten des Strahls durch 2 mm von Al und 0,1 & mgr; m Au filtriert. Ein CdWO 4 Szintillator mit einem 1,25X Objektivlinse und einer PCO EDGE Kamera wurden in der Detektoranordnung verwendet. Die trockenen Scans wurden mit 4000 Projektionen erworben, während die dynamischen Scans jeweils 1.000 Projektionen hatte. Gesamterfassungszeit betrug ~ 1 Minute 15 Sekunden pro Scan mit ~ 100 Scans über einen Zeitraum von 2 Stunden in Anspruch genommen.
Wiederaufbau und Artefaktentfernung wurde mit der Diamond Light Source proprietäre Software abgeschlossen. Jedes Bild besteht aus 2000 3 Voxel, which wurden dann binned Signal zu Rauschen in einem Bild führt bei einer Auflösung von 1000 3 Voxeln zu erhöhen von 6,1 um (Figur 5). Die Bilder wurden dann verarbeitet unter Verwendung der Bildverarbeitungsmodule in Avizo 8.1 und ImageJ Programme (siehe Zusatz File). Jedes Bild etwa 12 CPU Stunden und 3 Stunden GPU der Verarbeitung auf einem Computer mit einem 3,0-GHz-CPU und einem Tesla K20C GPU erforderlich.
Die segmentierten Bilder wurden als Zeitreihe für die Porosität Veränderungen analysiert, indem die Anzahl von Voxeln an Poren und Rock zählen. Während des Auflösens Porosität mit der Zeit zunimmt (Abbildung 6). Sichtkontrolle der segmentierten Bilder (Bild 7) zeigt das Vorhandensein eines Kanals in Strömungsrichtung. Wenn die Porosität als Funktion von Zeit und Entfernung von der Probe aufgetragen Einlass klar ist , dass ein Kanal in der ersten Stunde gebildet und verbreitert dann als das Experiment weiter (F ild 8).
Die segmentierten Bilder wurden dann als Eingabe in ein Netzwerk Extraktionsmodell verwendet , um Permeabilitätsänderungen (Abbildung 9) analysieren. Es wurde festgestellt, dass es einen starken Anstieg der Durchlässigkeit bei der ersten Stunde war, aber dann ist die Durchlässigkeit zu späteren Zeiten stabilisiert.
Abbildung 1. Die in situ - Versuchsapparatur. CO 2 wird durch die Einspritzpumpe und verwendet unter Druck stehende Salzlösung in den Reaktor zu äquilibrieren. Reactive Sole wird durch Kernanordnung durch die Empfangspumpe gezogen. Die Zelle wird von DI-Wasser in den Eingrenzungspumpe und erhitzt unter Verwendung Heizband durch ein Thermoelement in der Sperrflüssigkeit gesteuert beschränkt. Das experimentelle System ist miteinander verbunden unter Verwendung Schlauch und der Fluidstrom wird Ventile lenkbar Verwendung (V) und Unions (U).ps: //www.jove.com/files/ftp_upload/53763/53763fig1large.jpg "target =" _ blank "> Bitte hier klicken, um eine größere Version dieser Figur zu sehen.
Abbildung 2. Die Röntgenspektren der Diamond Light Source I-13I rosa Balken unter Verwendung berechnet sowohl die experimentelle Stimmkurve und theoretischen Spiegelreflexionsvermögen und Filterübertragung. Spiegel absorbieren Energien über 30 keV; Al und Au Filter absorbieren Energien unter 13 und 22 keV sind. Bitte klicken Sie hier , um eine größere Version dieser Figur zu sehen.
Abbildung 3. Die Abbildungsvorrichtung Strahllinie. Die Al und Au Blätter filtern, um die rosa Balken und die übrigen Röntgenstrahlen treffen the Kernanordnung. Ein Teil der Röntgenstrahlen werden von der Probe absorbiert, während der Rest durch die Probe passieren und auf den Szintillator, die im sichtbaren Spektrum fluoresziert. Dieses sichtbare Licht wird dann durch das Objektiv auf das CCD fokussiert, die das Licht in einem pixelierten Digitalbild übersetzt, wo die Pixelintensitätswert eine Funktion von der Anzahl der Röntgenstrahlen, die von dem Szintillator absorbiert werden. Bitte klicken Sie hier , um eine größere Version dieser Figur zu sehen.
Abbildung 4. Die Kernanordnung innerhalb des Kernhalters. PEEK Schlauch wird an den Innen Endarmaturen befestigt ist und durch die Stahl Endkappen gefädelt. Der Kern wird in Aluminiumfolie eingewickelt und in die Hülse eingeführt. Die Hülse wird dann über den Endarmaturen gestreckt, um eine wasserdichte zu schaffen seal und zwei weitere Schichten aus Aluminiumfolie zugegeben alles an seinem Platz zu halten und Gasdiffusion verhindern. Das Thermoelement ist an der Außenseite der Kernbaugruppe mit der äußeren Klebstoffschicht Aluminiumfolie befestigt. Abbildung von Menke et al modifiziert. 42. Bitte klicken Sie hier , um eine größere Version dieser Figur zu sehen.
Figur 5. Ein 2-D Schnitt des rekonstruierten Bildes vor (a) und nach (b) Auflösung. Die helleren Bereiche sind Getreide und die dunkleren Bereiche sind Pore. Unschärfen an den Rändern des Korns / Porengrenze kann in dem umgesetzten Teil des Porenraums (b) zu sehen ist. Bitte klicken Sie hier , um eine größere Version dieser Figur zu sehen.
Abbildung 6. Porösität aufgetragen mit der Zeit. Porosität erhöht sich linear mit einer geringen Abnahme in der Steigung in der zweiten Stunde der Auflösung. Abbildung von Menke et al modifiziert. 42. Bitte klicken Sie hier , um eine größere Version dieser Figur zu sehen.
Abbildung 7. Ein 3-D - Rendering der Änderung der Porosität nach 60 Minuten in das Experiment, in dem Grün die größte Veränderung der Porosität darstellt und rot am wenigsten. Eine klare poröse Kanal erstellt von Fluid-Feststoff chemische Reaktion in der Mitte des Kerns gesehen, wo die Auflösung am größten ist. Abbildung von Menke et al modifiziert."xref"> 42. Bitte klicken Sie hier , um eine größere Version dieser Figur zu sehen.
Abbildung 8. Profile der Porosität als Funktion der Entfernung von dem Probeneinlass. Porosität gleichförmig ist längs der Achse der Auflösung, aber die Auflösungsgeschwindigkeit ändert sich als Funktion der Zeit. Abbildung von Menke et al modifiziert. 42. Bitte klicken Sie hier , um eine größere Version dieser Figur zu sehen.
Figur 9 (A) ein Extraktionsnetz auf der segmentierten Bilder durchgeführt wird , 60 Minuten gezeigt, zeigtgroße Porenräume (Kugeln) und deren Verbindungen (Röhren). (B) Die berechnete Permeabilität gezeigt mit der Zeit mit einem starken Anstieg zwischen 40 und 60 Minuten als eine große Auflösung Kanal zu erhöhen hergestellt wird. Abbildung von Menke et al modifiziert. 42. Bitte klicken Sie hier , um eine größere Version dieser Figur zu sehen.
Discussion
Die wichtigsten Schritte für die dynamische Abbildung von Reaktions in heterogenen Porenstrukturen bei Lagerstättenbedingungen sind: 1) eine genaue Temperatursteuerung der Zelle innerhalb des rosa Strahls; 2) erfolgreich Kernhalter Stabilität auf einem sich schnell bewegenden Bühne; 3) eine effiziente Datenverarbeitung und Speichertechniken; und 4) wirksame Segmentierung von zeitaufgelösten Bilder.
Die Temperaturregelung ist von wesentlicher Bedeutung für Reservoir Zustand Bildgebung ein rosa Balken verwenden. Wenn die Temperatur über die Reaktortemperatur steigt, 2 CO wird im Porenraum exsolve und sowohl der pH - Wert der Salzlösung verändern und schaffen ganglia von überkritischem CO 2 in den Porenraum, der die Art der Auflösung 44 ändern könnte. Die Verwendung von Filtern, die niedrigere Energie Röntgenstrahlen zu absorbieren, ist von entscheidender Bedeutung für diese zusätzlichen Temperaturbelastung entfernt werden, die das Thermoelement und Heizungspackung ermöglicht extern Temperatur effektiv zu steuern. Allerdings senken Filter dieGesamtenergiedurchsatz des Strahls und somit müssen sparsam verwendet werden, um nicht wesentlich Gesamterfassungszeit zu erhöhen. Außerdem Filtertyp und Dicke müssen auf die spezifischen Energiewellenlängen und den Durchsatz des Strahllinie angepasst werden.
Der Kernhalter erfährt Rotations- und Vibrationsbelastungen während Tomographie Erwerb, dass die Kohlefasermanschette kann dazu führen, während der ersten Stufe Drehung und verwischen die Projektionen zu schütteln. Um dieses Potential zu minimieren, wird der Kernhalter als Kurz 6 cm Hülse für die Verwendung bei Synchrotrons. Diese Hülse würde für die Verwendung mit Tischscanner nicht günstig, da die Stahl Endarmaturen Minimierung der Quelle-Probe-Abstand und geometrische Vergrößerung hemmen würde. Jedoch mit einer parallelen Lichtquelle diese sind nicht Bedenken.
Jede tomographischen Scans in einer Serie erfolgt eine Größe haben über 20 GB was bedeutet, dass eine Reihe von 100 Scans 2 TB groß sein wird. Bei der Aufnahme viele Scans in einer Reihe sehr quickly sowohl das Instrument Bandbreite und Speicheroptionen bieten erhebliche Herausforderungen Managementdaten. Die experimentelle Bildgebungsgerät muss mit diesen Einschränkungen konzipiert werden, um voll und ganz die dynamische Bildgebung Potenzial schnell Tomographie realisieren. Die Datenübertragung von Engpässen müssen vor dem Start des Experiments identifiziert werden und die Technologie-Infrastruktur angepasst, so dass Themen wie Kamerageschwindigkeit ablesen, Übertragungsbandbreite und Speicherschreibgeschwindigkeit nicht Erfassungsgeschwindigkeit Potential hemmen.
Effektive Segmentierung der zeitaufgelösten Bilder der Auflösung stellt eine Herausforderung. Wenn ein tomographischen Scan in einem sich verändernden System genommen wird, kann die Ränder der Fest-Flüssig-Grenz verschwimmen. Diese Verwischung macht traditionelle Segmentierungstechniken wie Wendepunkt, der auf der Annahme arbeitet, dass die Grenzen werden die Regionen mit dem höchsten Dämpfungsverlauf, viel weniger erfolgreich. Um dies zu umgehen, das Differenzbild des unreacted und zur Reaktion gebracht Bilder berechnet, die ein Bild von nur Regionen des Wandels zur Verfügung stellt. Dieses Verfahren ermöglicht eine erfolgreiche Segmentierung der sich kontinuierlich ändernden Porenstruktur.
Synchrotron-Tomographie schnell mit einem Reservoir Skala Gerät gekoppelt ist eine leistungsfähige experimentelle Methode, die angepasst werden kann, eine Reihe von Anwendungen, einschließlich mehrphasigen Strömungsvorgänge, die Advektion-Dispersion und Transport in den chemisch heterogenen Medien zu erkunden. Allerdings ist die aktuelle Gerät zu einer Zeitauflösung in der Größenordnung von Sekunden, einphasigen Experimente und kleine Probengrößen beschränkt. Zukünftige Design-Upgrades können zusätzliche Pumpen sind für Drei-Phasen-Fähigkeiten, flux Erhöhung der Lage sein, größere Medien zu durchdringen, eine bessere Rekonstruktionstechniken, die für weniger Projektionen erlauben pro Scan und multivariate Ansätze zur Bildaufnahme und Segmentierung genommen werden, die weitere Informationen zu verbessern Tiefe, Breite und Genauigkeit.
Materials
| Name | Company | Catalog Number | Comments |
| NaCl salt | Sigma Aldrich | S7653-1KG | |
| KCl salt | Sigma Aldrich | P9333-1KG | |
| KI salt | Sigma Aldrich | 30315-1KG | |
| Core holder | Airbourne Composites | 110 mm Core holder | Constructed in conjunction with Imperial College |
| PEEK tubing | Kinesis | 1560xL | |
| Thermocouple | Omega Engineering | KMTSS-IM300U-150 | |
| Flexible Heating Tape | Omega Engineering | KH-112/10-P | |
| 1/16" Needle Valve | Hydrasun Ltd | MVE1002 | |
| High Pressure Syringe Pump | Teledyne ISCO | 1000D | |
| 600 mL Parr Reactor | Parr Instrument Company | 4547A - hastelloy | |
| CO2 Cylinder | BOC | CO2 - size E | |
| Viton | Fisher Scientific | 11572583 | |
| Aluminium Foil | Coroplast | 1510AWX | |
| ImageJ - image processing | NIH | ImageJ | |
| Matlab | Mathworks | Matlab | Used for data analysis |
| Avizo | FEI | Avizo | |
| Snoop Leak Detector | Swagelok | MS-SNOOP-8OZ |
References
- Herzog, H., Caldeira, K., Reilly, J. An issue of permanence: Assessing the effectiveness of temporary carbon storage. Clim. Change. 59, 293-310 (2003).
- Metz, B., Davidson, O., de Coninck, H., Loos, M., Meyer, L. Carbon Dioxide Capture and Storage. IPCC. , Cambridge University Press. 431 (2005).
- Langmuir, D., Hall, P., Drever, J. Aqueous Environmental Geochemistry. , Prentice Hall. (1997).
- Morse, J. W., Mackenzie, F. T. Geochemistry of sedimentary carbonates. , Elsevier. (1990).
- Peng, C., Crawshaw, J. P., Maitland, G. C., Trusler, J. M., Vega-Maza, D. The pH of CO 2-saturated water at temperatures between 308K and 423K at pressures up to 15MPa. J Supercrit Fluid. 82, 129-137 (2013).
- Peng, C., Crawshaw, J. P., Maitland, G. C., Trusler, J. M. Kinetics of calcite dissolution in CO 2-saturated water at temperatures between (323 and 373) K and pressures up to 13.8 MPa. Chem. Geol. 403, 74-85 (2015).
- Bachu, S., Nordbotten, J. M., Celia, M. A. Evaluation of the spread of acid gas plumes injected in deep saline aquifers in western Canada as an analogue to CO2 injection in continental sedimentary basins. Proceedings of 7th International Conference on Greenhouse Gas Control Technologies. 1, (2004).
- Bachu, S. Review of CO 2 storage efficiency in deep saline aquifers. Int J Greenh Gas Con. , (2015).
- Marland, G., Fruit, K., Sedjo, R. Accounting for sequestered carbon: the question of permanence. Environ Sci Policy. 4, 259-268 (2001).
- Daccord, G., Lenormand, R., Lietard, O. Chemical Dissolution of a Porous-Medium by a Reactive Fluid .1. Model for the Wormholing Phenomenon. Chem. Eng. Sci. 48, 169-178 (1993).
- Daccord, G., Lietard, O., Lenormand, R. Chemical Dissolution of a Porous-medium by a Reactive Fluid .2. Convection vs Reaction, Behavior Diagram. Chem. Eng. Sci. 48, 179-186 (1993).
- Maheshwari, P., Ratnakar, R., Kalia, N., Balakotaiah, V. 3-D simulation and analysis of reactive dissolution and wormhole formation in carbonate rocks. Chem. Eng. Sci. 90, 258-274 (2013).
- El-Maghraby, R., Pentland, C., Iglauer, S., Blunt, M. A fast method to equilibrate carbon dioxide with brine at high pressure and elevated temperature including solubility measurements. J Supercrit Fluid. 62, 55-59 (2012).
- Fredd, C., Fogler, S. Influence of Transport and Reaction on Wormhole Formations in Porous Media. AIChE. 44, (1998).
- Gharbi, O., Toth, A., Bijeljic, B., Boek, E., Blunt, M. PGE Seminar Series. , Imperial College London. (2013).
- Luquot, L., Gouze, P. Experimental determination of porosity and permeability changes induced by injection of CO2 into carbonate rocks. Chem. Geol. 265, 148-159 (2009).
- Cohen, C. E., Ding, D., Quintard, M., Bazin, B. From pore scale to wellbore scale: Impact of geometry on wormhole growth in carbonate acidization. Chem. Eng. Sci. 63, 3088-3099 (2008).
- Li, L., Peters, C. A., Celia, M. A. Upscaling geochemical reaction rates using pore-scale network modeling. Adv Water Resour. 29, 1351-1370 (2006).
- Swoboda-Colberg, N. G., Drever, J. I. Mineral dissolution rates in plot-scale field and laboratory experiments. Chem. Geol. 105, 51-69 (1993).
- Berg, S., et al. Real-time 3D imaging of Haines jumps in porous media flow. Proc. Natl. Acad. Sci. U.S.A. 110, 3755-3759 (2013).
- Blunt, M. J., et al.
- Ott, H., et al. Core-flood experiment for transport of reactive fluids in rocks. Rev. Sci. Instrum. 83, 084501 (2012).
- Gharbi, O. Fluid-Rock Interactions in Carbonates: Applications to CO2 storage. , Imperial College London. (2014).
- Noiriel, C., Gouze, P., Made, B. 3D analysis of geometry and flow changes in a limestone fracture during dissolution. J Hydrol. 486, 211-223 (2013).
- Hao, Y., Smith, M., Sholokhova, Y., Carroll, S. CO2-induced dissolution of low permeability carbonates. Part II: Numerical modeling of experiments. Adv Water Resour. 62, 388-408 (2013).
- Smith, M. M., Sholokhova, Y., Hao, Y., Carroll, S. A. CO2-induced dissolution of low permeability carbonates. Part I: Characterization and experiments. Adv Water Resour. 62, 370-387 (2013).
- Gouze, P., Luquot, L. X-ray microtomography characterization of porosity, permeability and reactive surface changes during dissolution. J. Contam. Hydrol. 120-121, 45-55 (2011).
- Menke, H. P., Bijeljic, B., Andrew, M. G., Blunt, M. J. Dynamic Three-Dimensional Pore-Scale Imaging of Reaction in a Carbonate at Reservoir Conditions. Environ. Sci. Technol. 49, 4407-4414 (2015).
- Andrew, M., Bijeljic, B., Blunt, M. Reservoir Condition Pore-scale Imaging of Multiple Fluid Phases Using X-ray Microtomography. J Vis Exp. , (2015).
- Andrew, M., Bijeljic, B., Blunt, M. J. Pore-scale imaging of trapped supercritical carbon dioxide in sandstones and carbonates. Int J Greenh Gas Con. 22, 1-14 (2014).
- Henke, B. L. Filter Transmission. , Available from: http://henke.lbl.gov/optical_constants/filter2.html (2015).
- Henke, B. L., Gullikson, E. M., Davis, J. C. X-ray interactions: photoabsorption, scattering, transmission, and reflection at E= 50-30,000 eV, Z= 1-92. At. Data Nucl. Data Tables. 54, 181-342 (1993).
- Schlüter, S., Sheppard, A., Brown, K., Wildenschild, D. Image processing of multiphase images obtained via X-ray microtomography: A review. Water Resour. Res. 50, 3615-3639 (2014).
- Buades, A., Coll, B., Morel, J. -M. A non-local algorithm for image denoising. Computer Vision and Pattern Recognition, 2005. CVPR 2005. IEEE Computer Society Conference on. 2, 60-65 (2005).
- Buades, A., Coll, B., Morel, J. -M. Nonlocal image and movie denoising. Int J Comput Vision. 76, 123-139 (2008).
- Sheppard, A. P., Sok, R. M., Averdunk, H. Techniques for image enhancement and segmentation of tomographic images of porous materials. Physica A. 339, 145-151 (2004).
- Lanczos, C. An iteration method for the solution of the eigenvalue problem of linear differential and integral operators. , United States Governm. Press Office. (1950).
- Andrew, M., Menke, H., Blunt, M. J., Bijeljic, B. The Imaging of Dynamic Multiphase Fluid Flow Using Synchrotron-Based X-ray Microtomography at Reservoir Conditions. Transport Porous Med. , 1-24 (2015).
- Raeini, A. Q., Blunt, M. J., Bijeljic, B. Modelling two-phase flow in porous media at the pore scale using the volume-of-fluid method. J. Comput. Phys. 231, 5653-5668 (2012).
- Bijeljic, B., Raeini, A., Mostaghimi, P., Blunt, M. J. Predictions of non-Fickian solute transport in different classes of porous media using direct simulation on pore-scale images. Phys Rev E. 87, 013011 (2013).
- Dong, H., Blunt, M. J. Pore-network extraction from micro-computerized-tomography images. Phys Rev E. 80, 036307 (2009).
- Menke, H. P., Bijeljic, B., Andrew, M. G., Blunt, M. J. Dynamic Pore-scale Imaging of Reactive Transport in Heterogeneous Carbonates at Reservoir Conditions. Energy Procedia. 63, 5503-5511 (2014).
- Bijeljic, B., Mostaghimi, P., Blunt, M. J. Signature of non-Fickian solute transport in complex heterogeneous porous media. Phys. Rev. Lett. 107, 204502 (2011).
- Ott, H., Oedai, S. Wormhole formation and compact dissolution in single-and two-phase CO2-brine injections. Geophys. Res. Lett. 42, 2270-2276 (2015).
