Przestrzennie rozłożona charakteryzacja zwilżalności w ośrodkach porowatych

Dlaczego drobne styki cieczy mają znaczenie dla wielkich pytań energetycznych

W jaki sposób ropa, woda i gazy takie jak wodór czy dwutlenek węgla naprawdę przemieszczają się i zostają uwięzione wewnątrz skał głęboko pod ziemią? Odpowiedź zależy od „zwilżalności” – tego, jak bardzo powierzchnia skały preferuje jedną ciecz nad drugą – na skalach zbyt małych, by zobaczyć je gołym okiem. W artykule opisano nową metodę mapowania tej preferencji w trzech wymiarach, pora po porze, ujawniając ukryty patchwork zachowań, który silnie wpływa na podziemne składowanie CO₂, magazynowanie wodoru i wydobycie ropy, a także na wydajność zaawansowanych baterii i ogniw paliwowych.

Patrzenie we wnętrze skał w trzech wymiarach

Nowoczesna mikro-tomografia komputerowa rentgenowska (micro‑CT) pozwala naukowcom zajrzeć w trójwymiarowe wnętrze małych próbek skalnych i zobaczyć, gdzie różne ciecze zajmują przestrzeń porową. Z tych obrazów próbują wyznaczyć „kąt zwilżania” tam, gdzie powierzchnia fazy ciekłej spotyka powierzchnię stałą — prostą miarę geometryczną zwilżalności. W teorii kąt ten mówi, czy skała zachowuje się jako hydrofilowa, oleofilowa czy gdzieś pośrodku. W praktyce odczytanie tego kąta bezpośrednio z obrazów jest niezwykle trudne: dokładna linia styku skały, wody i oleju jest rozmyta na kilka pikseli, powierzchnie są chropowate, a komputerowa „segmentacja” faz nigdy nie jest idealna. Istniejące zautomatyzowane metody często wygładzają istotne szczegóły, przesuwają interfejsy o kilka pikseli i przez to zniekształcają kąty, zwłaszcza w skomplikowanych skałach lub układach mieszanej zwilżalności, gdzie zachowanie silnie różni się między porami.

Nowy sposób śledzenia krawędzi cieczy

Autorzy przedstawiają zautomatyzowany algorytm geometryczny, który omija najtrudniejszy aspekt problemu: precyzyjne wyznaczenie trójfazowej linii kontaktu. Zamiast mierzyć dokładnie na tej niejednoznacznej linii, metoda buduje szczegółowe siatki powierzchni granicznych skała–ciecz i ciecz–ciecz, a następnie ekstrapoluje lokalne kierunki powierzchni (normalne) z sąsiednich, lepiej rozdzielonych obszarów w stronę kontaktu. Te normalne są łączone, by obliczyć kąty zwilżania w wielu punktach wzdłuż każdej pętli styku. Przebieg pracy obejmuje staranne redukowanie szumu, solidną metodę segmentacji, która podąża za rzeczywistymi krawędziami intensywności na obrazach, oraz łagodne wygładzanie powierzchni, które usuwa pikselowe „schodkowe” artefakty bez kurczenia czy deformowania kształtów porów. Wbudowana kontrola jakości odrzuca pomiarowe odstępstwa, które ewidentnie odbiegają od lokalnego otoczenia, kosztem liczby na rzecz wiarygodności.

Testy metody na idealnych kształtach i rzeczywistych skałach

Aby sprawdzić dokładność, zespół najpierw zastosował algorytm do w pełni syntetycznych zestawów danych: cyfrowych kropli osadzonych na płaskich i zakrzywionych powierzchniach stałych, gdzie prawdziwy kąt zwilżania jest dokładnie znany. W szerokim zakresie kątów i rozdzielczości obrazów nowa metoda odtwarzała prawdziwe wartości z dokładnością około pięciu stopni, działając lepiej i bardziej konsekwentnie niż powszechnie używane narzędzia, zwłaszcza przy małych i dużych kątach, gdzie błędy zwykle rosną. Badacze przeszli następnie do rzeczywistych obrazów micro‑CT z kilku typów skał, w tym wapieni i piaskowców, zawierających układy olej–woda i wodór–woda w warunkach przepływu istotnych dla wydobycia ropy i podziemnego magazynowania gazów. Porównując wyniki automatyczne z pracochłonnymi, ręcznymi pomiarami kątów, wykazali, że ich metoda dobrze odpowiada ekspertyzie ludzkiej, unikając silnych uprzedzeń wprowadzanych przez starsze automatyczne techniki, które nadmiernie wygładzają interfejsy.

Ujawnianie ukrytych wzorców w zachowaniu skała–ciecz

Dysponując tysiącami wiarygodnych lokalnych pomiarów, autorzy zbudowali trójwymiarowe mapy kąta zwilżania w całej przestrzeni porowej. W skałach ogólnie hydrofilowych kąty są stosunkowo jednorodne i niskie, co potwierdza, że woda ma skłonność do przylegania do powierzchni stałych i przepływania przez wąskie narożniki, podczas gdy olej zajmuje centra porów. W „mieszanej” piaskowcu, zmienionym przez długą ekspozycję na ropę surową, średni kąt sugeruje jedynie łagodnie hydrofilowe zachowanie. Jednak mapa przestrzenna opowiada bogatszą historię: niemal 60% powierzchni pozostaje hydrofilowych, podczas gdy około 40% przechodzi w reżim pośredni, gdzie żadna z cieczy nie jest wyraźnie preferowana. Te pośrednie plamy to dokładnie miejsca, gdzie na obrazach obserwuje się nietypowe interfejsy w kształcie siodeł i złożone wzory napełniania, wyjaśniając zagadkowe kombinacje uwięzienia i przepływu, których nie da się uchwycić pojedynczym, uśrednionym kątem zwilżania.

Dlaczego to ma znaczenie dla energetyki i środowiska

Dla inżynierów próbujących przewidzieć, jak CO₂ lub wodór rozprzestrzenią się i pozostaną uwięzione pod ziemią, lub jak elektrolity i gazy przemieszczają się przez porowate warstwy ogniw paliwowych i baterii, znajomość jednego „średniego” parametru zwilżalności już nie wystarcza. Badanie pokazuje, że subtelne, przestrzennie zmienne zachowania zwilżania kontrolują, jak ciecze wnikają, zostają zablokowane lub omijają niektóre ścieżki. Nowy algorytm dostarcza map zwilżalności pora po porze wraz z przejrzystą oceną niepewności pomiaru, w ramach oprogramowania open‑source. Mówiąc przystępnie, zamienia on rozmazane obrazy rentgenowskie skał w szczegółowe „mapy preferencji” dla cieczy, oferując potężne narzędzie do projektowania bezpieczniejszych projektów składowania CO₂, stabilniejszych rezerwuarów wodoru, ulepszonych strategii wydobycia ropy oraz bardziej wydajnych urządzeń elektrochemicznych.

Cytowanie: Aljaberi, F., Belhaj, H., Foroughi, S. et al. Spatially distributed wettability characterization in porous media. Sci Rep 16, 12643 (2026). https://doi.org/10.1038/s41598-026-43688-1

Słowa kluczowe: zwilżalność, ośrodki porowate, kąt zwilżania, składowanie CO2, składowanie wodoru