Artykuł metodologiczny

Dynamiczne obrazowanie stanu zbiornika w skali porowej reakcji w węglanach za pomocą szybkiej tomografii synchrotronowej

9.2K wyświetleń

⸱

DOI:

10.3791/53763

⸱

21 lutego 2017

 ,  ,  ,  ,  , 
1Department of Earth Science and Engineering, Imperial College London, 2Carl Zeiss X-Ray Microscopy, 3Diamond Manchester Imaging Branchline (I13-2), Diamond Lightsource

W tym artykule

Podsumowanie

Szybka tomografia synchrotronowa została użyta do dynamicznego zobrazowania rozpuszczania wapienia w obecności solanki nasyconej CO2 w warunkach zbiornikowych. Wykonano 100 skanów w rozdzielczości 6,1 μm przez okres 2 godzin.

Streszczenie

Trwałość składowania pod ziemią jest głównym problemem dla wychwytywania i składowania dwutlenku węgla. Pompowanie CO2 do zbiorników węglanowych może spowodować rozpuszczenie uszczelnień geologicznych i umożliwić ucieczkę CO2 . Jednak procesy rozpuszczania w warunkach zbiornikowych są słabo poznane. W związku z tym potrzebne są eksperymenty czasowo-rozdzielcze, aby obserwować i przewidywać naturę i szybkość rozpuszczania w skali porów. Szybka tomografia synchrotronowa to metoda wykonywania obrazów o wysokiej rozdzielczości w czasie złożonych struktur porowych znacznie szybciej niż tradycyjna μ-CT. Diamond Lightsource Pink Beam został użyty do dynamicznego obrazowania rozpuszczania wapienia w obecności solanki nasyconej CO2 w warunkach zbiornikowych. Wykonano 100 skanów w rozdzielczości 6,1 μm w ciągu 2 godzin. Obrazy podzielono na segmenty, a porowatość i przepuszczalność zmierzono za pomocą analizy obrazu i ekstrakcji sieciowej. Porowatość wzrastała równomiernie na całej długości próbki; Jednak tempo wzrostu zarówno porowatości, jak i przepuszczalności uległo spowolnieniu w późniejszych czasach.

Wprowadzenie

Głównym problemem związanym z wychwytywaniem i składowaniem dwutlenku węgla (CCS) jest długoterminowe bezpieczeństwo składowania1,2. Dwutlenek węgla, CO2, wstrzyknięty pod powierzchnię rozpuści się w solance gospodarza i utworzy kwas węglowy3,4,5. Ta kwaśna solanka może reagować i rozpuszczać otaczającą skałę, szczególnie jeśli skałą macierzystą jest wapień6. Rozpuszczanie może być korzystne i pozwolić na ciągłą przepuszczalność formacji7 i większą trwałość przechowywania8. Jednak integralność uszczelnień geologicznych może być zagrożona przez to rozpuszczenie i umożliwić migrację CO2 na powierzchnię9. Dokładne modelowanie predykcyjne trwałości magazynowania zależy zatem od pełnego zrozumienia rozpuszczania w systemie solanka-skała oraz rozkładu i szybkości ruchu płynu w podpowierzchniowym 10,11,12.

Jednak charakter i szybkość rozpuszczania w węglanach zależy zarówno od właściwości solanki13,14,15,16 oraz od skały macierzystej17. Szybkość rozpuszczania jest również silnie zależna od temperatury i ciśnienia solanki6, co sprawia, że rozwój eksperymentalnych technik pomiaru złożonych procesów zależnych od czasu w reprezentatywnych warunkach złoża ma kluczowe znaczenie.

Poprzednie eksperymenty wykazały, że szybkość reakcji w skali pola jest zazwyczaj o rzędy wielkości niższa niż w eksperymentalnych pomiarach reaktora wsadowego18,19. Wietrzenie, niejednorodność minerałów i niepełne mieszanie w niejednorodnym polu przepływu są możliwymi wyjaśnieniami tego zjawiska. Nie jest jednak możliwa ocena najważniejszych czynników bez bezpośredniej obserwacji ewoluującej przestrzeni porowej podczas reakcji. W związku z tym potrzebne są dynamiczne eksperymenty w skali porów, aby zapewnić zarówno wgląd w wzajemne zależności między transportem a reakcją, jak i zweryfikować modele predykcyjne.

Uznaną metodą eksperymentalną do badania procesów skalowania porów w zastosowaniach związanych z magazynowaniem węgla jest mikrotomografia rentgenowska (μCT)20,21. μ-CT ma kilka zalet: osiąga wysokie rozdzielczości przestrzenne do około 1 μm, jest nieinwazyjna i dostarcza trójwymiarowe obrazy. Rozpuszczanie wapienia badano w skali rdzenia (~cm)22 i stwierdzono, że reakcja skała-solanka zwiększa niejednorodność fizyczną. Aby pogłębić zrozumienie, w jaki sposób różne warunki transportu i reakcji zmieniają złożone struktury ciał stałych i porów, konieczne jest zmierzenie wywołanych reakcją zmian w geometrii przestrzeni porowej, topologii i przepływie w podpowierzchniowych systemach skalnych w temperaturach i ciśnieniach zbiorników oraz w wyższej rozdzielczości, w celu szczegółowego zbadania procesów w skali porów. W artykule opisano metodę badania reaktywnych procesów rozpuszczania w skałach o złożonej strukturze porów i skoncentrowano się na pomiarze czasu i szybkości reakcji zależnej od czasu i przestrzennie między solanką zakwaszoną CO2 a skałą wapienną w warunkach złożowych.

Przeprowadzono kilka badań, które dotyczyły reakcji w złożonych węglanach23,24,25,26,27, ale ze względu na ograniczenia eksperymentalne lub obrazowe, zostały one albo ograniczone do obrazów przed i po reakcji, albo nie zostały ukończone w reprezentatywnych warunkach podpowierzchniowych. Menke i wsp.28 przeprowadził dynamiczne obrazowanie in situ reakcji pomiędzy solanką zakwaszoną CO2 a wapieniem kettońskim w skali porów przez okres kilku godzin oraz w temperaturze i ciśnieniu reprezentatywnym dla warstwy wodonośnej na głębokości około 1 km. Jednak Ketton jest stosunkowo jednorodną skałą o dużych ziarnach, którą łatwo zobrazować w bardzo krótkim czasie (~17 min) i przy niewielkiej liczbie projekcji (~400). Większość skał węglanowych ma złożone struktury porów, które wymagają wielu projekcji w celu dokładnego rozdzielenia, co może być bardzo czasochłonnym procesem przy użyciu tradycyjnej μ-CT - albo z monochromatyczną wiązką w źródle synchrotronowym, albo z laboratoryjnymi skanerami rentgenowskimi. W związku z tym potrzebna jest szybka metoda tomografii, aby dynamicznie obserwować zmiany w heterogenicznych węglanach wywołane reakcją.

Ilość czasu potrzebna do zobrazowania próbki jest kontrolowana przez strumień źródła promieniowania rentgenowskiego. Jedną z metod szybkiego skanowania jest użycie polichromatycznej wiązki źródła synchrotronowego20. Ta tak zwana "różowa wiązka" zapewnia o rzędy wielkości intensywniejsze światło niż źródła stacjonarne, a zatem obrazy mogą być wykonywane w skali kilkudziesięciosekundowej, a nie godzinnej. Undulator, który składa się z okresowej struktury magnesów dipolowych, wytwarza różową wiązkę. Wiązka elektronów jest zmuszona do ulegania oscylacjom podczas przechodzenia przez magnesy, w wyniku czego wypromieniowuje energię. Wytwarzana energia jest skoncentrowana w wąskich pasmach długości fal i jest bardzo intensywna. Lustra i filtry są następnie wykorzystywane do zawężania spektrum światła w celu dostosowania go do potrzeb eksperymentalnych. Lustra pochłaniają widmo o wysokiej energii, podczas gdy filtry pochłaniają niższe energie. Dzięki temu możliwe jest zawężenie widma do pożądanego pasma promieniowania przy użyciu tylko tych narzędzi.

Jednak korzystanie z tego intensywnego strumienia promieniowania rentgenowskiego nie jest pozbawione wyzwań. Promienie rentgenowskie o niższej energii widma różowej wiązki są pochłaniane przez próbkę w postaci ciepła. Może to zakłócać kontrolę temperatury aparatu in situ i powodować rozpuszczenie CO2 z roztworu20. Solanka nasycona CO2 jest bardzo wrażliwa zarówno na ciepło, jak i ciśnienie, dlatego niewielka zmiana równowagi termicznej może znacząco zmienić pH płynu in situ5. W związku z tym staranne elementy projektowe i kontrolne dla widma promieniowania rentgenowskiego muszą zostać włączone do wyposażenia linii wiązki przed obrazowaniem.

Szybka tomografia również generuje ogromne ilości danych z dużą szybkością. Ograniczenia związane z odczytem danych z kamery i późniejszym przechowywaniem stanowią poważne wyzwanie technologiczne. Niektórzy poradzili sobie z tym problemem, wykonując kilka kolejnych skanów i zapisując je w pamięci aparatu przed odczytaniem ich na zewnętrznych serwerach danych. Wymaga to jednak, aby eksperyment był stosunkowo krótki, ponieważ pamięć kamery może pomieścić tylko skończoną ilość danych. Binning danych w aparacie skraca również czas przesyłania, ponieważ zmniejsza ilość danych, które trzeba przesłać, ale może potencjalnie obniżyć jakość obrazów. Alternatywnie, dane mogą być przesyłane z aparatu po każdym skanowaniu przed rozpoczęciem następnego, co wydłuży całkowity czas między skanowaniami. W tym badaniu wykorzystano tę drugą metodę, przy czym każda akwizycja obrazu trwała ~45 sekund, a odczyt danych zajmował dodatkowe ~30 sekund.

Podczas wykonywania skanów z dużą szybkością, próbka musi obracać się znacznie szybciej niż przy tradycyjnym skanowaniu, dlatego potencjalne naprężenie kątowe na uchwycie rdzenia jest duże. Włókno węglowe, choć przezroczyste dla promieni rentgenowskich, jest elastyczne pod wpływem naprężeń. Jeśli próbka porusza się podczas akwizycji obrazu, może wystąpić rozmycie obrazu. Tuleja uchwytu rdzenia została zaprojektowana tak, aby była jak najkrótsza, aby złagodzić te potencjalne naprężenia. Dodatkowo na wszystkich elementach aparatury doświadczalnej w pobliżu stopnia zastosowano elastyczne rurki z polieteroeteroketonu (PEEK), tak aby stopień mógł się swobodnie obracać. Jedną z wad stosowania rurek PEEK jest to, że przepuszcza CO2 w dyfuzyjnych skalach czasowych. Płyn znajdujący się w przewodach przez długi czas będzie stopniowo ulegał desaturacji w ciągu około 24 godzin. Wszystkie przewody, które nie znajdowały się w pobliżu uchwytu rdzenia, zostały wykonane ze stali nierdzewnej, a płyn został wstępnie zrównoważony w energicznie zmieszanym reaktorze Hastelloy, podgrzanym i podgrzanym do warunków eksperymentalnych23,29,30.

Eksperymentalna aparatura jest przedstawiona na rysunku 1. Temperatura zbiornika jest utrzymywana w uchwycie rdzenia poprzez owinięcie zewnętrznej części tulei przezroczystą taśmą grzewczą rentgenowską i włożenie termopary przez otwór promieniowy ogniwa do płynu ograniczającego. Następnie regulator pochodnej proporcjonalnej całki (PID) regulował temperaturę z dokładnością do 1 °C. Warunki ciśnienia i przepływu utrzymywano za pomocą trzech wysokociśnieniowych pomp strzykawkowych, które osiągają prędkość przepływu 0,001 ml/min. Do eksperymentu użyto dwóch soli, wysoce absorbującej 25% wag. niereaktywnej solanki KI i nisko absorbującej 1% wag. KCl, 5% wag. Reaktywnej solanki NaCl. Różnica w tłumieniu sprawiła, że łatwo było dostrzec pojawienie się reaktywnej solanki w rdzeniu, dzięki czemu obliczenia objętości martwej stały się zbędne.

Protokół

1. Projektowanie strategii obrazowania

  1. Oblicz widma promieniowania rentgenowskiego linii badawczej przy najwyższej energii i strumieniu różowej wiązki, aby przewidzieć wydajność obrazowania przy użyciu eksperymentalnej krzywej strojenia i pomiaru transmisji filtrów. Przykładowe widma promieniowania rentgenowskiego Diamond Lightsource I13-2 'Pink Beam' pokazano na rysunku 2.
  2. Ponieważ promienie rentgenowskie o niższej energii powodują nagrzewanie próbki i nie zwiększają kontrastu obrazowania, należy odfiltrować dolną część widma promieniowania rentgenowskiego tak, aby do obrazowania próbki używane były tylko promienie rentgenowskie o najwyższej energii. Wybierz filtry liniowe z materiałów, które pochłaniają światło o pożądanych niskich długościach fal odpowiednich dla źródła światła, obliczając teoretyczną transmisję filtra przy dostępnych długościach fal światła31,32. Tutaj użyj aluminium i złota jako linii wiązki przy tym źródle światła.
    1. Stosować filtr pasmowo-przepustowy składający się z kombinacji filtrów rentgenowskich jako górnoprzepustowych i lustra rentgenowskiego działającego w pobliżu kąta krytycznego jako filtra dolnoprzepustowego. W takim przypadku należy użyć zestawu filtrów węglowych pirolitycznych 0,2 mm i aluminiowych 0,2 mm, a do lustra paska pokrytego platyną pod kątem padania 1,15 mrad. Lustro odbija tylko światło poniżej 30 keV, a dodatkowe filtry liniowe zostały zainstalowane o średnicy 2 mm Al i 0,1 μm Au, które mają piki absorpcyjne odpowiednio na poziomie 13 i 22 keV, aby filtrować więcej promieniowania rentgenowskiego o niższej energii. Rysunek 3 przedstawia aparaturę do obrazowania linii badawczej.
  3. Wybierz scyntylator, który obficie mieni się na dostępnych liniach badawczych, częstotliwościach i strumieniu światła. W tym przypadku przesiewacz scyntylacyjny jest wykonany z wolframu kadmu o grubości 250 μm (CdWO4), który jest ułożony w stos z wolframem ołowiowym o grubości 750 μm (PbWO4). Następnie wybierz obiektyw i kamerę, które mają odpowiednie pole widzenia i rozdzielczość czasu przyciągania dla wymagań eksperymentalnych. W takim przypadku połącz obiektyw 1,25x o aperturze 0,04 z kamerą PCO EDGE 5,5 CMOS i użyj go do uchwycenia 4-milimetrowego pola widzenia z liczbą klatek na sekundę 0,001 s.
  4. Wybierz technikę "flyscan" do akwizycji obrazu, ponieważ ta metoda obracania stolika zmniejsza drgania próbki. Tradycyjna akwizycja wymaga, aby stolik zatrzymywał się na każdym kroku kątowym, wykonywał rzut, a następnie przechodził do następnego kąta. Akwizycja obrazu podczas tych dynamicznych tomografii została wykonana za pomocą "flyscanu", który wykonuje skan tomograficzny w miarę ruchu stolika i przyjmuje przyrost kątowy taki, że różnica między każdą kolejną projekcją jest niewielka. Metoda "flyscan" eliminuje niewielkie efekty drgań ruchu start-and-stop i szybciej zapewnia wyższą jakość obrazu.

2. Montaż sprzętu i komórki

  1. Załaduj rdzeń do ogniwa, przygotowując się do zalania rdzenia.
    1. Najpierw owiń rdzeń jedną warstwą folii aluminiowej i włóż do rękawa (np. Vitonu) (Rysunek 4).
    2. Przyciąć tuleję na wymiar tak, aby była o 2 mm krótsza niż łączna długość rdzenia i wewnętrznych łączników końcowych. Złączki końcowe to krajowe złączki rurowe 1/16" z gwintem rurowym (NPT) do złączek, które zostały obrobione do średnicy zewnętrznej 5 mm, podczas gdy tuleja ma średnicę wewnętrzną 4 mm.
    3. Naciągnij tuleję na złączki końcowe 5 mm, aby uzyskać szczelne uszczelnienie. Upewnij się, że między łącznikami końcowymi a rdzeniem nie ma żadnej przestrzeni, aby upewnić się, że ciśnienie ograniczające nie spowoduje nadmiernego ściśnięcia tulei i nie przyciśnie przepływu.
    4. Owiń złączki i tuleję dwiema dodatkowymi warstwami aluminium, zarówno aby zapobiec dyfuzji gazowego CO2 do płynu ograniczającego, jak i utrzymać tuleję na miejscu na złączkach i zapobiec łączeniu płynów ograniczających i porowych przez ścieżkę hydrauliczną.
    5. Złóż z powrotem uchwyt rdzenia, wsuwając rurkę i uszczelki z powrotem na miejsce, a następnie uszczelnij zaślepki i łączniki końcowe, wymieniając.
  2. Zamontuj uchwyt rdzenia na stoliku i podłącz przewody przepływowe i elektryczne.
  3. Przetestuj obrót stolika i upewnij się, że wszystkie przewody przepływowe i elektryczne mogą się swobodnie obracać w zakresie od -90° do 90°.
  4. Wykonaj suchy skan całego rdzenia przed rozpoczęciem eksperymentu.
    1. Zeskanuj rdzeń w zachodzących na siebie odcinkach o szerokości i długości około 4 mm. Skalibruj czas ekspozycji skanowania do średniej wartości zliczania około 15 000, co zapewnia wysoki stosunek sygnału do szumu bez nadmiernego nasycenia scyntylatora. Wykonaj każdy suchy skan z co najmniej 2,400 projekcjami, aby zachować kontrast fazowy i ostrość krawędzi.
    2. Wykonaj płaskie i ciemne zdjęcia scyntylatora, aby można było uwzględnić wszelkie uszkodzenia i szumy zewnętrzne podczas rekonstrukcji. Zrób płaskie powierzchnie, przesuwając uchwyt rdzenia poza pole view i robiąc zdjęcie samego scyntylatora z włączoną wiązką. Weź ciemne kolory tą samą metodą przy wyłączonej wiązce.

3. Zwiększanie ciśnienia w systemie

  1. Załaduj solankę z chlorku potasu (KCl) o zawartości 1% wag. i chlorku sodu (NaCl) o stężeniu wagowym (NaCl) do zdemontowanego reaktora, wlewając płyn przez górną część zbiornika reaktora.
    1. Dodaj sproszkowaną skałę węglanową, aby uzyskać pożądaną kwasowość solanki. W tym przypadku nie dodano węglanu.
    2. Zmontuj reaktor, dokręcając i owijając go taśmą grzewczą, a następnie wkładając sondę temperatury do góry.
    3. Załadować CO2 do pompy wtryskowej, otwierając zawór 1 (V1 na rysunku 1).
    4. Zamknij zawór 1 i zwiększ ciśnienie w pompie wtryskowej do 100 barów.
    5. Otwórz zawór 2, aby zalać reaktorCO2. Podgrzać reaktor do 50 °C za pomocą owijki grzewczej sterowanej PID w połączeniu z sondą temperatury i ciągle mieszać za pomocą mieszadła porywającego napędzanego zewnętrznym silnikiem elektrycznym. Równoważyć solankę z CO2 w temperaturze 10 MPa i 50 °C przez okres od 2 do 6 godzin, aby upewnić się, że solanka jest całkowicie nasycona CO2 i że węglan jest całkowicie rozpuszczony.
  2. Przed podłączeniem uchwytu rdzenia należy całkowicie oczyścić układ z powietrza i ewentualnych środków strącających w przewodach z poprzednich eksperymentów. Aby to zrobić, połącz linie powyżej i poniżej uchwytu rdzenia, aby ominąć uchwyt rdzenia (U1 i U2).
    1. Załaduj wodę dejonizowaną (DI) do pompy odbiorczej przez zawór 11, ustawiając pompę odbiorczą na ponowne napełnienie.
    2. Otwórz zawory 7, 4 i 3 i użyj pompy odbiorczej w trybie stałego ciśnienia, aby przekierować wodę DI z powrotem przez system i wyprowadzić z zaworu 3 poniżej reaktora. Użyj około dziesięciu objętości systemu, aby upewnić się, że przewody są wolne od powietrza i wypłukane do czysta.
  3. Opróżnij pompę odbiorczą, a następnie załaduj solankę o stężeniu 25% wag. KI do pompy odbiorczej przez zawór 11 i załaduj wodę DI do pompy ograniczającej przez zawór 10.
    1. Zamknij zawór 10 i otwórz zawory 8 i 6. Użyj pompy ograniczającej, aby ograniczyć rdzeń do 2 MPa.
    2. Zamknij zawór 11 i zwiększ ciśnienie w pompie odbiorczej do 10 barów.
    3. Otwórz zawory 9, 7, 4 i 3 i wykorzystaj powstały spadek ciśnienia do przepuszczenia solanki domieszkowanej KI przez rdzeń.
    4. Stopniowo zwiększaj ciśnienie ograniczające i porowe, aż do ustalenia rozsądnego natężenia przepływu. Przepuścić około dwóch pełnych objętości solanki systemu przez rdzeń i spuścić płyn przez zawór 3 poniżej reaktora. W ten sposób całe powietrze jest usuwane z systemu, a rdzeń jest zalewany solanką o wysokim kontraście, co sprawia, że nadejście nieoczyszczonej solanki reaktywnej jest łatwe do zaobserwowania.
    5. Zamknij zawór 3 i stopniowo zwiększaj ciśnienie ograniczające i porowe, aż rdzeń zostanie ograniczony do 12 MPa, a ciśnienie porowe wyniesie 10 MPa. Włączyć regulator PID, aby doprowadzić rdzeń do 50 °C.
    6. Zatrzymaj pompę odbiorczą, zamknij zawór 3 i otwórz zawór 5 u podstawy reaktora, aby połączyć system reaktora z rdzeniem.

4. Przepływ płynów i akwizycja obrazu

  1. Wyśrodkuj środek rdzenia w polu widzenia i wykonuj projekcje 2D w sposób ciągły, gdy rdzeń jest zalewany, aby śledzić postęp zalewania rdzenia. Wykonuj projekcje 2D, włączając różową wiązkę i używając aparatu do robienia zdjęć bez obracania stolika. Rozpoczęcie projekcji 2D przed wstrzyknięciem płynu reaktywnego daje wyraźny obraz przed solanką, który zostanie później porównany z kolejnymi obrazami wypełnionymi solanką.
  2. Ustaw pompę odbiorczą tak, aby napełniała się z żądanym natężeniem przepływu, ciągnąc w ten sposób płyn z reaktora przez rdzeń w żądanych warunkach przepływu, pozostawiając pompę wtryskową w celu regulacji ciśnienia od przodu.
    1. Monitoruj projekcje 2D pod kątem zmian tłumienia, które sygnalizują nadejście reaktywnej solanki. Kiedy dotrze reaktywna solanka, transmisja rdzenia wzrośnie, a projekcje 2D znacznie pojaśnieją, gdy więcej światła uderzy w scyntylator, a domieszkowana solanka zostanie wyparta przez wysoce przezroczysty płyn reaktywny o wysokim promieniowaniu rentgenowskim. Jeżeli nie ma różnicy w tłumieniu między solanką reaktywną i niereaktywną, wówczas, w zależności od widma linii wiązki, może być wymagane ponowne rozpoczęcie eksperymentu od kroku 2.1 z wyższym stężeniem soli KI lub przy użyciu innej soli o wysokiej absorpcji.
    2. Zatrzymaj skany 2D i tomografie 3D kolejno tak szybko, jak pozwala na to aparatura do obrazowania. Użyj około 1 000 projekcji na skan. Skanuj rdzeń, obracając go tylko o 180° (w przeciwieństwie do tradycyjnego 360°). Użycie mniejszej liczby stopni obrotu zmniejsza stosunek sygnału do szumu, ale jest szybsze i pomaga uniknąć rozciągania i plątania przewodów przepływowych i elektrycznych. Wykonuj skany 3D, dopóki nie zostanie osiągnięty limit czasu lub rdzeń będzie wyglądał na wystarczająco rozpuszczony, aby istniało bezpośrednie niebezpieczeństwo wewnętrznego zapadnięcia się struktury (a tym samym spowodowania utraty zarówno ciśnienia ograniczającego, jak i przyszłych danych z suchego skanowania całego rdzenia).
  3. Po wykonaniu ostatniego skanowania należy skutecznie rozhermetyzować system, aby uniknąć dalszej reakcji rdzenia.
    1. Najpierw zatrzymaj pompę odbiorczą. Następnie zamknij zawór 5 łączący reaktor z resztą układu.
    2. Zmniejsz ciśnienie w układzie za pomocą pomp ograniczających i odbiorczych, utrzymując ciśnienie na płynie ograniczającym o około 1 MPa.
    3. Po osiągnięciu ciśnienia atmosferycznego w granicach 1 MPa otwórz pompę ograniczającą i odbiorczą za pomocą zaworów 10 i 11 i pracuj w trybie stałego przepływu, aby spuścić pozostały płyn.
    4. Wyłącz regulator PID i otwórz złącze 4-kierunkowe (U2) w górnej części uchwytu rdzenia, aby uwolnić pozostałe ciśnienie w systemie.
    5. Powoli poluzuj linię ograniczającą, jednocześnie łapiąc nadmiar wody ograniczającej DI za pomocą chłonnego papieru. Zamknij zawory 6 i 7 i odłącz złącze 1 oraz przewody elektryczne.
    6. Poluzuj stage clamp i wyjmij uchwyt rdzenia ze stolika.
  4. Ostrożnie wyjmij zespół rdzenia z uchwytu rdzenia, a następnie odłącz tuleję od wewnętrznych łączników końcowych. Nie wyjmuj rdzenia z tulei, ponieważ może to spowodować uszkodzenie delikatnego rdzenia, który przeszedł reakcję. Umieść rdzeń pokryty tuleją w zlewce pełnej wody demineralizowanej, aby rozcieńczyć potencjalnie reaktywną solankę i zatrzymać całą reakcję.
  5. Suszyć cały rdzeń w suszarce o temperaturze 60 °C przez co najmniej 12 godzin. Następnie ponownie zamontuj rdzeń na stoliku za pomocą tradycyjnego mocowania próbki i zeskanuj go ponownie w tej samej rozdzielczości i projekcjach, co początkowy skan na sucho.

5. Przetwarzanie obrazu

  1. Skoryguj zrekonstruowane obrazy pod kątem utwardzenia wiązki związanej z użyciem wiązki polichromatycznej, zakładając, że wszelkie oddziaływania są promieniowo symetrycznymi funkcjami Gaussa33.
  2. Filtruj obrazy za pomocą filtra zachowującego krawędzie, takiego jak niezlokalizowane sposoby zwiększania sygnału do szumu34,35 (patrz Plik uzupełniający).
  3. Podziel obrazy skanowania na sucho za pomocą algorytmu segmentacji zlewni36 i zdefiniuj ziarna jako skałę i pustkę (patrz Plik uzupełniający).
    1. Wykonaj pierwsze zdjęcie rdzenia z reaktywną solanką i zarejestruj każdy kolejny obraz na pierwszym obrazie, a następnie ponownie próbkuj go za pomocą metody ponownego próbkowania Lanczos37 z pierwszym obrazem jako odniesieniem. Ponieważ trwająca reakcja ma tendencję do rozmywania krawędzi, segmentacja działu wodnego na obrazach nie jest wystarczająca do dokładnej segmentacji.
    2. Odejmij każdy przereagowany obraz rdzeniowy od pierwszego obrazu, aby uzyskać obraz różnicy. Podziel obrazy różnic na zmiany i brak zmian. Zarejestruj podzielony na segmenty skan suchy do pierwszego skanowania reaktywnego, a następnie odejmij zmianę segmentacji od segmentowanego skanowania na sucho, aby uzyskać segmentowane obrazy przereagowane38.

6. Modelowanie

  1. Binaryzowane obrazy można wykorzystać jako dane wejściowe do bezpośredniego solvera przepływu Naviera-Stokesa39,40, lub modelu ekstrakcji sieciowej41 (Rysunek 8), aby scharakteryzować zmiany przepuszczalności i zapewnić fizyczny wgląd w dynamikę rozpuszczania.

Wyniki

Reakcję obrazowano pomiędzy kalcytem a nienasyconą solanką nasyconą scCO2 w rdzeniu węglanu portlandzkiego o średnicy 4 mm i długości 1,2 cm42. Węglan portlandzki jest stosunkowo czystym (<99%) oolitem kalcytowym o złożonej, heterogenicznej strukturze porów43. Promienie rentgenowskie niskiej energii filtrowano, przepuszczając wiązkę przez 2 mm Al i 0,1 µm Au. W układzie detektora zastosowano scyntylator CdWO4 z obiektywem 1,25X oraz kamerę PCO EDGE. Skanowanie na sucho wykonano przy użyciu 4 00 projekcji, natomiast skany dynamiczne obejmowały po 1 0 projekcji każdy. Całkowity czas akwizycji wynosił ~1 minutę 15 sekund na skan, przy czym w ciągu 2 godzin wykonano ~100 skanów.

Rekonstrukcja i usuwanie artefaktów zostały wykonane przy użyciu autorskiego oprogramowania Diamond Lightsource. Każdy obraz składał się z 203 wokseli, które następnie poddano binowaniu w celu zwiększenia stosunku sygnału do szumu, co poskutkowało uzyskaniem obrazu o rozmiarze 1003 wokseli przy rozdzielczości 6.1 µm (Rycina 5). Następnie obrazy przetworzono przy użyciu modułów przetwarzania obrazu w programach Avizo 8.1 i ImageJ (patrz Plik uzupełniający). Przetwarzanie każdego obrazu zajęło około 12 godzin CPU i 3 godziny GPU na komputerze z procesorem 3.0 GHz i kartą graficzną Tesla K20C GPU.

Zasegmentowane obrazy przeanalizowano jako szeregi czasowe w celu określenia zmian porowatości, licząc liczbę wokseli porów oraz skały. Podczas rozpuszczania porowatość wzrasta wraz z upływem czasu (Rysunek 6). Inspekcja wizualna zasegmentowanych obrazów (Rysunek 7) wykazuje obecność kanału w kierunku przepływu. Po wykreśleniu porowatości jako funkcji czasu oraz odległości od wlotu próbki wyraźnie widać, że kanał tworzy się w pierwszej godzinie, a następnie rozszerza się w miarę trwania eksperymentu (Rysunek 8).

Kolejno wykorzystano segmentowane obrazy jako dane wejściowe do modelu ekstrakcji sieci w celu analizy zmian przepuszczalności (Rysunek 9). Stwierdzono gwałtowny wzrost przepuszczalności w ciągu pierwszej godziny, po czym w późniejszym czasie przepuszczalność uległa stabilizacji.

Schemat układu przepływu cieczy z pompami, reaktorem i ogniwem Hasslera; wykorzystywany w eksperymentach z wtryskiwaniem CO2.
Rycina 1. in situ aparatura doświadczalna tlenek węgla(II)2 jest pod ciśnieniem wytwarzanym przez pompę wtryskową i służy do równoważenia solanki w reaktorze. Reaktywna solanka jest przeciągana przez zespół rdzenia za pomocą pompy odbiorczej. Komora jest ograniczona wodą dejonizowaną (DI) w pompie ograniczającej i podgrzewana za pomocą taśmy grzejnej sterowanej termoparą w płynie ograniczającym. System eksperymentalny jest połączony za pomocą przewodów, a przepływ płynu jest kierowany za pomocą zaworów (V) i złączek (U). Kliknij tutaj, aby wyświetlić powiększoną wersję tej figury.

Wykres strumienia fotonów w funkcji energii, analiza lustra Au-Al, piki absorpcji spektralnej.
Rysunek 2. Widma rentgenowskie wiązki różowej (pink beam) z Diamond Lightsource I-13I obliczone z wykorzystaniem zarówno eksperymentalnej krzywej strojenia, jak i teoretycznej odbijalności zwierciadeł oraz transmitancji filtrów. Lustra pochłaniają energie powyżej 30 keV; filtry z Al i Au pochłaniają energie odpowiednio poniżej 13 i 2 keV. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Schemat linii wiązki synchrotronowej; undulator, filtry, CCD; do obrazowania rentgenowskiego, analizy próbek.
Rysunek 3. Aparatura obrazująca linii wiązki. Arkusze Al i Au filtrują różową wiązkę, a pozostałe promienie rentgenowskie uderzają w zespół rdzenia. Część promieni rentgenowskich jest pochłaniana przez próbkę, natomiast pozostałe przechodzą przez nią i uderzają w scyntylator, który fluorescuje w widmie widzialnym. Światło widzialne jest następnie skupiane przez obiektyw na detektorze CCD, który przekształca to światło w zpikselowany obraz cyfrowy, w którym wartość natężenia piksela jest funkcją liczby promieni rentgenowskich pochłoniętych przez scyntylator. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Schemat głównego zestawu z termoparą, rurką PEEK i rdzeniem węglanowym do badań obrazowych.
Rycina 4. Zespół rdzenia wewnątrz uchwytu rdzenia. Przewody PEEK są przymocowane do wewnętrznych złączy końcowych i przewleczone przez stalowe zaślepki końcowe. Rdzeń zostaje owinięty folią aluminiową i włożony do tulei. Następnie tuleja jest naciągnięta na złącza końcowe, aby uzyskać wodoszczelną uszczelkę, a następnie dodaje się dwie dodatkowe warstwy folii aluminiowej, aby utrzymać wszystkie elementy na miejscu i zapobiec dyfuzji gazów. Termopara jest przymocowana do zewnętrznej strony zespołu rdzenia za pomocą zewnętrznej warstwy klejącej folii aluminiowej. Rysunek zmodyfikowany na podstawie pracy Menke. i wsp.42. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Obraz mikroskopowy, struktura porów skały osadowej, analiza porównawcza, wyniki badań geologicznych.
Rysunek 5. Przekrój 2D zrekonstruowanego obrazu przed (a) i po (b) rozpuszczeniu. Jasne obszary to ziarna, a ciemne obszary to pory. Rozmycie na krawędziach granicy ziarno/pora jest widoczne w przereagowanej części przestrzeni porowej (b). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Wykres porowatości w funkcji czasu ilustrujący trend wzrostowy, istotny dla badań nad analizą materiałową.
Rysunek 6. Porowatość w funkcji czasu. Porowatość zwiększa się liniowo z niewielkim spadkiem nachylenia w drugiej godzinie rozpuszczania. Rysunek zmodyfikowany na podstawie Menke i wsp.42. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Schemat analizy mikrostruktury gleby, przedstawiający rozkład cząstek na potrzeby badania składu gleby.
Rysunek 7. Trójwymiarowy rendering zmiany porowatości w 60. minucie eksperymentu; kolor zielony oznacza największą zmianę porowatości, a czerwony najmniejszą. W centrum rdzenia, gdzie rozpuszczanie jest najintensywniejsze, widoczny jest wyraźny porowaty kanał powstały w wyniku chemicznej reakcji ciało stałe-płyn. Rysunek zmodyfikowano na podstawie pracy Menke i wsp.42. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Wykres porowatości w funkcji odległości przedstawiający dane dla 0, 60 i 120 minut; analiza wykresu liniowego.
Rysunek 8. Profile porowatości w funkcji odległości od wlotu próbki. Porowatość jest jednorodna wzdłuż osi rozpuszczania, jednak szybkość rozpuszczania zmienia się w funkcji czasu. Rysunek zmodyfikowany na podstawie pracy Menke i wsp.42. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Symulacja sieci polimerowej i wykres przepuszczalności; analiza dyfuzji w czasie.
Rycina 9. (A) Ekstrakcja sieci przeprowadzona na obrazach po segmentacji jest przedstawiona dla 60 minuty i ukazuje duże przestrzenie porów (kule) oraz ich połączenia (rurki). (BWykazano, że obliczona przepuszczalność wzrasta w czasie, przy czym gwałtowny wzrost następuje między 40. a 60. minutą, co wiąże się z powstaniem szerokiego kanału rozpuszczania. Rysunek zmodyfikowano na podstawie pracy Menke. i inni42. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Dyskusja

Najbardziej krytycznymi etapami dynamicznego obrazowania reakcji w heterogenicznych strukturach porowych w warunkach zbiornika są: 1) dokładna kontrola temperatury komórki wewnątrz różowej wiązki; 2) udana stabilność uchwytu rdzenia na szybko poruszającej się scenie; 3) efektywne techniki przetwarzania i przechowywania danych; oraz 4) efektywna segmentacja obrazów rozdzielczych w czasie.

Kontrola temperatury jest niezbędna do obrazowania stanu zbiornika za pomocą różowej wiązki. Jeśli temperatura wzrośnie powyżej temperatury reaktora, CO2 rozpuści się w przestrzeni porów i zarówno zmieni pH solanki, jak i utworzy zwoje nadkrytycznego CO2 w przestrzeni porów, które mogą zmienić charakter rozpuszczania44. Zastosowanie filtrów do pochłaniania promieni rentgenowskich o niższej energii ma kluczowe znaczenie dla usunięcia tego dodatkowego stresu temperaturowego, co pozwala termoparze i owijce grzewczej skutecznie kontrolować temperaturę na zewnątrz. Jednak filtry obniżają całkowitą przepustowość energii wiązki i dlatego muszą być używane oszczędnie, aby nie wydłużyć znacznie całkowitego czasu akwizycji. Ponadto typ i grubość filtra muszą być dostosowane do określonych długości fal energii i przepustowości linii wiązki.

Uchwyt rdzenia jest poddawany naprężeniom obrotowym i wibracyjnym podczas akwizycji tomografii, które mogą powodować drgania tulei z włókna węglowego podczas obracania stolika i rozmycie projekcji. Aby zminimalizować ten potencjał, uchwyt rdzenia został zaprojektowany jako krótka tuleja o długości 6 cm do stosowania w synchrotronach. Tuleja ta nie byłaby korzystna do stosowania ze skanerami stołowymi, ponieważ stalowe łączniki końcowe hamowałyby minimalizację odległości źródło-próbka i powiększenie geometryczne. Jednak w przypadku równoległego źródła światła nie są to obawy.

Każdy skan tomograficzny wykonany w serii może mieć rozmiar ponad 20 GB, co oznacza, że seria 100 skanów będzie miała rozmiar 2 TB. W przypadku bardzo szybkiego wykonywania wielu skanów z rzędu, zarówno przepustowość urządzenia, jak i opcje pamięci masowej stanowią poważne wyzwanie w zakresie zarządzania danymi. Eksperymentalna aparatura do obrazowania musi być zaprojektowana z uwzględnieniem tych ograniczeń, tak aby w pełni wykorzystać potencjał obrazowania dynamicznego szybkiej tomografii. Przed rozpoczęciem eksperymentu należy zidentyfikować wąskie gardła w przesyłaniu danych, a infrastrukturę technologiczną dostosować tak, aby takie kwestie, jak szybkość odczytu z kamery, przepustowość transferu i prędkość zapisu w pamięci masowej, nie hamowały potencjału szybkości akwizycji.

Skuteczna segmentacja rozdzielczych w czasie obrazów rozpadu stanowi wyzwanie. Gdy skan tomograficzny jest wykonywany w zmieniającym się układzie, krawędzie granicy ciało stałe i ciecz mogą ulec rozmyciu. To rozmycie sprawia, że tradycyjne techniki segmentacji, takie jak zlewnia, która działa przy założeniu, że granice będą regionami o najwyższym gradiencie tłumienia, są znacznie mniej skuteczne. Aby obejść ten problem, obliczany jest obraz różnicowy obrazów nieprzereagowanych i przereagowanych, co daje obraz tylko obszarów zmian. Metoda ta pozwala na skuteczną segmentację stale zmieniającej się struktury porów.

Szybka tomografia synchrotronowa w połączeniu z aparaturą w skali zbiornikowej jest potężną metodą eksperymentalną, którą można dostosować do zbadania szeregu zastosowań, w tym procesów przepływu wielofazowego, adwekcji-dyspersji i transportu w chemicznie niejednorodnych ośrodkach. Jednak obecna aparatura jest ograniczona do rozdzielczości czasowej rzędu sekund, eksperymentów jednofazowych i małych rozmiarów próbek. Przyszłe ulepszenia konstrukcyjne mogą obejmować dodatkowe pompy do pracy trójfazowej, zwiększenie strumienia, aby móc penetrować większe media, lepsze techniki rekonstrukcji, które pozwalają na wykonanie mniejszej liczby projekcji na skan, oraz wielowymiarowe podejścia do pozyskiwania i segmentacji obrazu, które mogą jeszcze bardziej poprawić głębokość, szerokość i dokładność informacji.

Oświadczenia

Autorzy nie mają nic do ujawnienia.

Podziękowania

Z wdzięcznością dziękujemy za finansowanie z Qatar Carbonates and Carbon Storage Research Centre (QCCSRC), dostarczone wspólnie przez Qatar Petroleum, Shell i Qatar Science & Technology Park. Wyrażamy również wdzięczność za finansowanie i wsparcie zapewnione przez Diamond Lightsource i Uniwersytet w Manchesterze w I13 Imaging Branch.

Materiały

Lista materiałów użytych w tym artykule
NazwaFirmaNumer katalogowyKomentarze
Sól NaClSigma AldrichS7653-1KG
Sól KClSigma AldrichP9333-1KG
Sól KISigma Aldrich30315-1KG
Uchwyt rdzeniaAirbourne Composites110 mm Uchwyt rdzeniaZbudowany w połączeniu z
rurką PEEKImperial College Kinesis1560xL
TermoparaOmega EngineeringKMTSS-IM300U-150
Elastyczna taśma grzewczaOmega EngineeringKH-112/10-P
1/16 "Zawór iglicowyHydrasun LtdMVE1002
Wysokociśnieniowa pompa strzykawkowaTeledyne ISCO1000D
600 mL ReactorParr Instrument Company4547A - hastelloy
CO2 ButlaBOCCO2 - rozmiar E
VitonFisher Scientific11572583
Folia aluminiowaCoroplast1510AWX
ImageJ - przetwarzanie obrazuNIHImageJ
MatlabMathworksMatlabSłuży do analizy danych
AvizoFEIAvizo
Wykrywacz nieszczelności SnoopSwagelokMS-SNOOP-8OZ

Bibliografia

  1. Herzog, H., Caldeira, K., Reilly, J. An issue of permanence: Assessing the effectiveness of temporary carbon storage. Clim. Change. 59, 293-310 (2003).
  2. Metz, B., Davidson, O., de Coninck, H., Loos, M., Meyer, L. Carbon Dioxide Capture and Storage. IPCC. , Cambridge University Press. 431(2005).
  3. Langmuir, D., Hall, P., Drever, J. Aqueous Environmental Geochemistry. , Prentice Hall. (1997).
  4. Morse, J. W., Mackenzie, F. T. Geochemistry of sedimentary carbonates. , Elsevier. (1990).
  5. 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).
  6. 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).
  7. 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).
  8. Bachu, S. Review of CO 2 storage efficiency in deep saline aquifers. Int J Greenh Gas Con. , (2015).
  9. Marland, G., Fruit, K., Sedjo, R. Accounting for sequestered carbon: the question of permanence. Environ Sci Policy. 4, 259-268 (2001).
  10. 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).
  11. 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).
  12. 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).
  13. 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).
  14. Fredd, C., Fogler, S. Influence of Transport and Reaction on Wormhole Formations in Porous Media. AIChE. 44, (1998).
  15. Gharbi, O., Toth, A., Bijeljic, B., Boek, E., Blunt, M. PGE Seminar Series. , Imperial College London. (2013).
  16. 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).
  17. 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).
  18. Li, L., Peters, C. A., Celia, M. A. Upscaling geochemical reaction rates using pore-scale network modeling. Adv Water Resour. 29, 1351-1370 (2006).
  19. Swoboda-Colberg, N. G., Drever, J. I. Mineral dissolution rates in plot-scale field and laboratory experiments. Chem. Geol. 105, 51-69 (1993).
  20. 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).
  21. Blunt, M. J., et al. Pore-scale imaging and modelling. Adv Water Resour. 51, 197-216 (2013).
  22. Ott, H., et al. Core-flood experiment for transport of reactive fluids in rocks. Rev. Sci. Instrum. 83, 084501(2012).
  23. Gharbi, O. Fluid-Rock Interactions in Carbonates: Applications to CO2 storage. , Imperial College London. (2014).
  24. 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).
  25. 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).
  26. 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).
  27. 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).
  28. 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).
  29. Andrew, M., Bijeljic, B., Blunt, M. Reservoir Condition Pore-scale Imaging of Multiple Fluid Phases Using X-ray Microtomography. J Vis Exp. , (2015).
  30. 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).
  31. Henke, B. L. Filter Transmission. , Available from: http://henke.lbl.gov/optical_constants/filter2.html (2015).
  32. 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).
  33. 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).
  34. 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).
  35. Buades, A., Coll, B., Morel, J. -M. Nonlocal image and movie denoising. Int J Comput Vision. 76, 123-139 (2008).
  36. 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).
  37. Lanczos, C. An iteration method for the solution of the eigenvalue problem of linear differential and integral operators. , United States Governm. Press Office. (1950).
  38. 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).
  39. 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).
  40. 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).
  41. Dong, H., Blunt, M. J. Pore-network extraction from micro-computerized-tomography images. Phys Rev E. 80, 036307(2009).
  42. 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).
  43. Bijeljic, B., Mostaghimi, P., Blunt, M. J. Signature of non-Fickian solute transport in complex heterogeneous porous media. Phys. Rev. Lett. 107, 204502(2011).
  44. Ott, H., Oedai, S. Wormhole formation and compact dissolution in single-and two-phase CO2-brine injections. Geophys. Res. Lett. 42, 2270-2276 (2015).

Przedruki i uprawnienia

Poproś o pozwolenie na ponowne wykorzystanie tekstu lub ilustracji tego artykułu JoVE

Poproś o pozwolenie

Tagi

Obrazowanie zbiornik w w glanowychobrazowanie reakcji w skali por wsolanka nasycona CO2warunki z o owekalibracja wi zki promieni rentgenowskichuk ad do prz przep ukiwania rdzenianaliza porowato ci i przepuszczalno cimodel ekstrakcji siecidynamiczny interfejs p yn ska a

Powiązane artykuły