Lattice-Boltzmann dla porowatych ośrodków: ponad 100 mln godzin GPU

Dlaczego drobne przestrzenie w skale mają znaczenie

Kiedy magazynujemy energię pod ziemią, wychwytujemy dwutlenek węgla lub eksploatujemy ogniwa paliwowe, sposób, w jaki dwa różne płyny przeciskają się przez drobne przestrzenie w skale, może zadecydować o powodzeniu technologii. Jednak obserwowanie tych płynów w rzeczywistych skałach jest wyjątkowo powolne i kosztowne. Niniejszy artykuł przedstawia ogromny, otwarty zestaw danych z symulacji komputerowych, który robi najbliższą rzecz: odtwarza, jak woda i płyn podobny do ropy przesuwają się względem siebie w realistycznych próbkach skał, wykorzystując ponad 100 milionów godzin pracy procesorów graficznych. Wynik to wspólne zasoby, z których każdy badacz może skorzystać, by testować pomysły dotyczące przepływów podziemnych bez potrzeby posiadania własnego superkomputera.

Badanie ruchu płynów w cyfrowych skałach

Autorzy skupiają się na „przepływie dwufazowym”, gdzie dwa niemieszające się płyny, na przykład woda i ropa, przemieszczają się wspólnie przez labirynt porów w skale. Kluczową wielkością dla inżynierów jest przepuszczalność względna, która mówi, jak łatwo każdy z płynów porusza się, gdy obecny jest drugi. Zwykle pełne pomiary wymagają tygodni starannych badań laboratoryjnych dla każdej próbki skały i zestawu warunków. Zamiast tego zespół wykorzystał wyspecjalizowane oprogramowanie symulacyjne, LBPM, aby obliczać przepływ bezpośrednio na trójwymiarowych obrazach prawdziwych skał. Te cyfrowe skały pochodziły z tomografii rentgenowskiej mikro-CT piaskowców i spieków szklanych, odwzorowując realistyczne kształty i rozmiary porów sięgające kilku mikrometrów.

Ogromny eksperyment wirtualny

Uruchamianie realistycznych symulacji w tak drobiazgowym wymiarze nadal jest bardzo kosztowne. Zespół wykorzystał superkomputery z tysiącami GPU, żeby przebadać warunki, które w laboratorium byłyby niepraktyczne. Zmieniali, jak silnie powierzchnia skały preferuje jeden płyn nad drugim (jej zwilżalność) oraz jak silnie płyny są przepychane przez ośrodek (opisane parametrem zwanym liczbą kapilarną). Dla czterech różnych materiałów porowatych wykonali zarówno protokoły przepływu „stacjonarnego”, jak i „niestacjonarnego”, naśladujące standardowe doświadczenia przepływowe stosowane w przemyśle naftowym. W sumie wygenerowali 50 pełnych krzywych przepuszczalności względnej oraz ponad 25 000 odrębnych konfiguracji płynów.

Widzieć kształty, nie tylko średnie

Powyżej średnich prędkości przepływu, symulacje śledzą szczegółowe kształty i powiązania płynów w czasie — informacje niemal niemożliwe do uzyskania na dużą skalę w eksperymentach. LBPM oddziela połączone drogi przepływu od izolowanych pęcherzy, czasem nazywanych gangliami, i mierzy ich objętość, pole powierzchni, krzywiznę, łączność, ciśnienie oraz ruch. Te wielkości są zapisywane na każdym kroku symulacji w prostych tabelach tekstowych, dzięki czemu użytkownicy mogą odtworzyć, jak powstają uwięzione kieszenie, jak się zrywają, łączą lub powoli opróżniają. Porównując różne typy skał i wzorce zwilżalności, zbiór danych pokazuje, jak subtelne zmiany preferencji powierzchni mogą przesunąć miejsca uwięzienia płynu i wpłynąć na jego przemieszczenie, co pomaga wyjaśnić obserwowane w laboratorium tendencje.

Wnioski o zwilżonej skale i uwięzionych płynach

Wykorzystując zestaw danych, autorzy weryfikują sensowność swoich symulacji, sprawdzając znane wzorce. Na przykład w miarę jak skała staje się bardziej „olejolubna”, ilość ropy pozostającej po przepłukaniu wodą zwykle maleje, co zgadza się z wcześniejszymi eksperymentami. W przypadkach bardziej hydrofilowych, przy słabo rozdzielonych wąskich gardłach, symulacje pokazują wczesne rozdrabnianie fazy niezwilżającej i niezwykle spłaszczone krzywe przepływu, ilustrując, jak rozdzielczość obrazu może wprowadzać w błąd. W innych skałach z lepiej odwzorowanymi porami przesunięcia krzywych przy zmianie zwilżalności zgadzają się z wcześniejszymi badaniami laboratoryjnymi. Miary łączności potwierdzają, że w pewnych pośrednich warunkach obie fazy pozostają silnie połączone, stan powiązany z złożonymi kształtami interfejsów i długo żyjącymi uwięzionymi skupiskami.

Wspólna podstawa dla przyszłych badań

Mówiąc prosto, artykuł dostarcza szczegółową mapę tego, jak dwa płyny mogą dzielić te same drobne przestrzenie w skale, w szerokim zakresie warunków, wszystko zakodowane w wielokrotnie użytecznych danych. Symulacje są porównane z najnowocześniejszymi eksperymentami obrazowania rentgenowskiego i zorganizowane tak, żeby inni mogli łatwo obliczać nowe krzywe podsumowujące, trenować modele uczenia maszynowego lub testować teorie dotyczące tego, jak szczegóły skali porów przekładają się na przepływ w dużej skali. Dla każdego, kto pracuje nad podziemnym magazynowaniem dwutlenku węgla, magazynowaniem wodoru, oczyszczaniem wód gruntowych lub komponentami ogniw paliwowych, ten otwarty zestaw danych oferuje potężne skróty: zamiast zaczynać od zera, można budować bezpośrednio na ponad 100 milionach godzin GPU starannie przygotowanych eksperymentów cyfrowych.

Cytowanie: Armstrong, R.T., Tavakkoli, O., Da Wang, Y. et al. Lattice-Boltzmann for Porous Media: 100M⁺ GPU Hours. Sci Data 13, 697 (2026). https://doi.org/10.1038/s41597-026-06823-1

Słowa kluczowe: porowate ośrodki, przepływ dwufazowy, digital rock, metoda sieci Boltzmanna, przepuszczalność względna