Modelowanie postępującej awarii stromych zboczy skalnych z użyciem połączonej metody elementów skończonych i dyskretnych

Dlaczego chwiejne ściany górskie mają znaczenie

Wysokie, skalne ściany wąwozów mogą wydawać się solidne, lecz w odpowiednich warunkach nagle się załamują, zrzucając miliony ton skał w dół. Dla ludzi mieszkających poniżej dużych zapór lub podróżujących górskimi drogami przewidywanie, kiedy takie zbocza mogą zawieść, to kwestia bezpieczeństwa, planowania i kosztów. W tym badaniu zaprezentowano nowy sposób cyfrowego "testowania wytrzymałości" stromych zboczy, który pozwala inżynierom zobaczyć nie tylko czy zbocze jest niebezpieczne, ale dokładnie jak i kiedy może się rozpaść.

Ukryte słabości w wyniosłych klifach

Zbocza gór rzadko składają się z jednorodnej, bezskazowej skały. Warstwy piaskowca i łupka, stare spękania oraz strefy wietrzenia tworzą mozaikę pasm mocnych i słabych. W zachodnich Chinach, zwłaszcza wzdłuż głęboko wciętych kanionów górnego biegu rzeki Lancang, projekty hydrotechniczne są wznoszone bezpośrednio na takich stromych ścianach. Choć te zbocza mogą odkształcać się powoli pod własnym ciężarem, mogą ulec katastrofalnej awarii, gdy pęknięcia połączą się wzdłuż wrażliwych warstw. Tradycyjne narzędzia inżynierskie potrafią oszacować ogólny margines bezpieczeństwa, ale nie oddają w pełni, jak drobne rysy zaczynają się u podnóża zbocza, rozrastają w górę i w końcu uwalniają poruszającą się masę skał.

Ograniczenia starszych narzędzi prognostycznych

Inżynierowie od dawna korzystają z trzech głównych typów metod komputerowych do oceny stabilności zboczy. Pierwsza, metoda równowagi granicznej, równoważy siły napędowe i oporowe wzdłuż założonej powierzchni osuwania; jest szybka, lecz wymaga wcześniejszego zgadywania kształtu tej powierzchni i opisuje jedynie moment awarii. Druga, metoda elementów skończonych, śledzi odkształcenia skały ciągłej, ale słabo radzi sobie z nagłym pojawieniem się i wzrostem pęknięć. Trzecia, metoda elementów dyskretnych, przedstawia zbocze jako oddzielne bloki, które mogą przesuwać się i zderzać, lecz nie opisuje wcześniejszego etapu, kiedy skała zachowuje się jako masa ciągła. Żadna z tych metod samodzielnie nie potrafi płynnie przeprowadzić symulacji od pozornie nienaruszonej skały do rozrzuconych głazów na dnie doliny.

Cyfrowe zbocze, które łamie się realistycznie

Autorzy łączą zalety tych podejść w ramach zwanej metodą elementów skończonych–dyskretnych (FDEM) i łączą ją ze strategią "zwiększania grawitacji". W ich wirtualnym zboczu skała jest przedstawiona jako wiele niewielkich, solidnych kawałków spajanych wzdłuż niewidocznych spoin. W miarę wzrostu symulowanego naprężenia te spoiwa mogą słabnąć i pękać, przekształcając ciągłą skałę w oddzielne bloki, które ślizgają się i zderzają. Zamiast zgadywać powierzchnię awarii model stopniowo zwiększa efektywne działanie grawitacji, aż zbocze wykazuje nagły skok przemieszczeń i energii kinetycznej. Poziom grawitacji w tym punkcie daje współczynnik bezpieczeństwa, a ewoluujący wzór pęknięć i ruchu pokazuje dokładnie, jak przebiega awaria.

Przetestowanie metody w realnym kanionie

Aby sprawdzić, czy cyfrowe zbocze zachowuje się jak w rzeczywistości, zespół zamodelował 796‑metrowe zbocze skalne obok elektrowni wodnej na górnym biegu rzeki Lancang. Zbudowali uproszczony przekrój uwzględniający różne strefy wietrzenia i istotną uskokę, a następnie użyli zmierzonych właściwości skały do ustalenia parametrów modelu. W miarę stopniowego zwiększania grawitacji symulacja odtworzyła realistyczną sekwencję: najpierw w najbardziej wietrzejącej skale u podnóża pojawiły się drobne rysy; rysy te połączyły się w ciągłe słabe pasmo; potem duża masa skalna ześlizgnęła się wzdłuż tego pasa, rozpadła na bloki i zatrzymała się jako szeroki usyp w dolinie. Przebyte odległości, głębokości spękań i końcowa szerokość rumowiska w dużym stopniu zgadzały się z obserwacjami terenowymi, a różnice zwykle nie przekraczały 20 procent.

Porównanie ze standardowymi metodami

Naukowcy porównali swoje nowe podejście z powszechnie stosowanymi metodami równowagi granicznej i elementów dyskretnych, obliczając, jak stabilne wydawało się to samo zbocze w każdym przypadku. Wszystkie trzy podejścia dały zbliżone współczynniki bezpieczeństwa i zgadzały się co do ogólnego położenia głównej powierzchni osuwania wzdłuż słabo zwięzłej, silnie wietrzejącej warstwy. Kluczowa różnica polegała na tym, że nowe podejście FDEM–grawitacja robiło znacznie więcej niż potwierdzało jedynie, czy zbocze jest marginalnie stabilne. Ujawniało też moment wystąpienia pierwszego ruchu, rolę uskoku w wywoływaniu wtórnych zawałów z wyższych partii oraz sposób, w jaki bloki oddziaływały i fragmentowały się podczas zjazdu — szczegóły istotne przy projektowaniu wzmocnień i planowaniu systemów monitoringu.

Znaczenie dla bezpieczniejszych inwestycji górskich

Autorzy dochodzą do wniosku, że połączone ramy modelowania mogą wiarygodnie śledzić pełną historię osuwiska, od pierwszych włoskowatych pęknięć do końcowej sterty gruzu. Dla zbocza na rzece Lancang obliczony współczynnik bezpieczeństwa sugeruje, że ściana skalna jest tylko nieznacznie stabilna w warunkach obecnych, co oznacza konieczność dodatkowych wzmocnień i ostrożnych wykopów. W szerszym ujęciu metoda daje inżynierom narzędzie do wskazania, gdzie najwcześniej może się pojawić uszkodzenie, gdzie umieścić kotwy i instrumenty monitorujące oraz jak przyszłe czynniki wyzwalające, takie jak trzęsienia ziemi czy intensywne opady, mogą zmienić sytuację. Chociaż obecne symulacje są dwuwymiarowe i wymagające obliczeniowo, rozszerzenie ich na trzeci wymiar i dodanie realistycznych czynników wyzwalających może uczynić je potężnym elementem długoterminowego zarządzania bezpieczeństwem w stromych obszarach górskich.

Cytowanie: Xu, J., Deng, Z., Feng, Y. et al. Modeling progressive failure in steep rock slopes using the combined finite-discrete element method. Sci Rep 16, 11180 (2026). https://doi.org/10.1038/s41598-026-40966-w

Słowa kluczowe: osuwiska, stabilność zboczy skalnych, modelowanie numeryczne, elektrownie wodne, grawitacyjne uszkodzenie