W kierunku trójwymiarowego dyskretnego modelowania sieci uskoków przy użyciu zintegrowanych wielowymiarowych danych z odsłonięć

Dlaczego podziemne pęknięcia mają znaczenie

Ukryte pod naszymi stopami sieci drobnych pęknięć w skale mogą decydować o tym, czy woda podziemna przepływa swobodnie, czy ropa i gaz mogą być wydobywane efektywnie, albo czy wstrzyknięty pod ziemię dwutlenek węgla pozostanie na miejscu. Jednak te sieci pęknięć trudno zobaczyć w trzech wymiarach, ponieważ odwierty i badania sejsmiczne dają tylko częściowe obrazy. W tym badaniu pokazano, jak szczegółowe zdjęcia skalnej ściany w Brazylii, w połączeniu z nowoczesnymi metodami obliczeniowymi i statystyką, można przekształcić w realistyczne modele 3D tych podziemnych pęknięć, poprawiając naszą zdolność przewidywania ruchu płynów w podłożu.

Naturalne okno do podłoża

Naukowcy skupili się na odsłonięciu Gaivota w basenie Potiguar w Brazylii, gdzie warstwy skał węglanowych formacji Jandaíra występują zarówno jako szeroka, niemal płaska powierzchnia, jak i stroma, pionowa ściana. Ta naturalna geometria daje dwa komplementarne widoki tego samego masywu skalnego: z góry, gdzie można mapować długie ślady pęknięć, oraz z boku, gdzie widoczne są całe płaszczyzny uskoków. Zespół wykorzystał zdjęcia z drona i fotogrametrię do zbudowania szczegółowego cyfrowego modelu 3D odsłonięcia, a następnie pieczołowicie zamapował ponad 1600 pęknięć powierzchniowych i prawie 500 płaszczyzn uskoków. Ponieważ podobne skały stanowią ważne złoża węglowodorów, miejsce to służy jako rzeczywiste laboratorium do zrozumienia, jak układy pęknięć kontrolują magazynowanie i przepływ płynów.

Od map pęknięć do trójwymiarowych rodzin szczelin

Przekształcenie tych obserwacji w użyteczny model 3D wymagało podziału pęknięć na sensowne „rodziny”, które mają podobne orientacje. Autorzy zastosowali metodę grupowania K-means, dostosowaną do danych leżących na sferze, aby pogrupować płaszczyzny uskoków 3D w cztery zestawy kierunków. Następnie przetestowali, jak zwarte i dobrze zrównoważone są te grupy, używając statystyki Fishera — narzędzia mierzącego, jak silnie kierunki skupiają się wokół średniej. Te rodziny kierunków posłużyły jako kręgosłup modelu: każda rodzina reprezentuje dominujący sposób, w jaki skała jest rozłamana, odzwierciedlając złożoną historię tektoniczną i krasową regionu.

Ujmowanie rozmiarów, wzorców i realistycznych objętości

Znajomość orientacji pęknięć to tylko połowa historii; modelarze potrzebują też realistycznych długości i rozstawów pęknięć. Dla pęknięć powierzchniowych zespół zbadał rozkład długości, testując kilka kandydackich form matematycznych, w tym prawa potęgowe i wykresy wykładnicze. Kluczowe parametry oszacowano za pomocą metody optymalizacji uczenia maszynowego zwanej losowym spadkiem gradientu (stochastic gradient descent). Większość rodzin pęknięć podążała za wzorcem prawa potęgowego, co oznacza, że jest wiele małych pęknięć i stopniowo coraz mniej dużych — cecha pękania rosnącego w sposób samopodobny, czyli fraktalny. Aby uniknąć budowania modeli zbyt małych, by były reprezentatywne, lub niepotrzebnie dużych, autorzy obliczyli także „reprezentatywną objętość elementarną” — minimalny rozmiar bloku, dla którego właściwości pęknięć, jak pole powierzchni na objętość, stają się stabilne. Ten krok zapewnił, że ich sześciany pęknięć 3D odzwierciedlają zachowanie uśrednione, a nie lokalne osobliwości.

Budowanie i testowanie syntetycznych światów pęknięć

Z orientacjami, rozkładami długości i reprezentatywną objętością w ręku, badacze wygenerowali dwa rodzaje modeli 3D pęknięć. Model pseudo-deterministyczny bezpośrednio włączał zamapowane płaszczyzny uskoków, przypisując długości losowane z dopasowanych rozkładów. Model w pełni stochastyczny tworzył nowe, losowe pęknięcia, które przestrzegały tych samych statystyk dla każdej rodziny i były dodawane aż do osiągnięcia docelowego pola powierzchni pęknięć na objętość. Oba modele następnie przecinano w wielu kierunkach, aby obliczyć standardowe miary pęknięć w 2D (długość na jednostkę pola) i 3D (pole na jednostkę objętości), a także stopień połączeń między pęknięciami. Porównanie wykazało, że modele syntetyczne wiernie odtwarzały intensywność i łączność pęknięć z rzeczywistego odsłonięcia, zwłaszcza gdy analizowano je rodzinami, a nie wszystkie naraz.

Łączenie wskazówek powierzchniowych z ukrytymi strukturami

Jednym z najbardziej praktycznych wniosków jest silne powiązanie między miarami pęknięć w 2D i 3D. Autorzy stwierdzili, że długość pęknięć na jednostkę pola mierzona na płaszczyźnie jest ściśle skorelowana z polem pęknięć na jednostkę objętości w otaczającej skale, z wartościami korelacji powyżej 0,9. Zaobserwowali też, że wraz ze wzrostem intensywności pęknięć łączność zwykle rośnie w podobny sposób, co sugeruje, że gęstsze sieci pęknięć zapewniają także bardziej ciągłe drogi przepływu płynów. Co ważne, te zależności wynikły z modeli opartych na rzeczywistych danych z odsłonięć, ale rozszerzonych na objętości znacznie wykraczające poza to, co jest bezpośrednio widoczne.

Co to oznacza dla wody, energii i magazynowania

Dla osób nietechnicznych kluczowy przekaz jest taki, że staranna integracja obrazów powierzchni o wysokiej rozdzielczości, geometrii odsłonięcia 3D i zaawansowanych narzędzi statystycznych pozwala przekształcić fragmentaryczne obserwacje pęknięć skalnych w solidne trójwymiarowe modele. Modele te pomagają przełożyć to, co geolodzy widzą na odsłoniętych klifach, na prognozy dotyczące ukrytych skał, które przechowują wodę pitną, ropę i gaz, czy wtryskiwane płyny odpadowe. Pokazując, że relatywnie dostępne pomiary 2D mogą wiarygodnie przewidywać właściwości i łączność pęknięć w 3D, praca ta oferuje praktyczny proces roboczy, który może poprawić projektowanie i bezpieczeństwo operacji podziemnych, od zarządzania wodami podziemnymi po produkcję energii i magazynowanie węgla.

Cytowanie: Racolte, G., Marques, A., Sales, V. et al. Towards three-dimensional discrete fracture network modeling using integrated multidimensional outcrop data. Sci Rep 16, 10087 (2026). https://doi.org/10.1038/s41598-026-37359-4

Słowa kluczowe: sieci uskoków, trójwymiarowe modelowanie geologiczne, złoża węglanowe, fotogrametria odsłonięć, przepływ płynów w podłożu