Mechanizmy rozkładu naprężeń zależne od rozmiaru modelu numerycznego i parametrów zrobów w eksploatacji wielowarstwowych pokładów węgla

Dlaczego naprężenia podziemne są ważne dla bezpieczeństwa wydobycia

Współczesne wydobycie węgla coraz częściej prowadzone jest głęboko pod ziemią, gdzie kilka warstw węgla zalega jedna nad drugą. Gdy jedna z tych warstw jest eksploatowana, otaczająca ją skała nie pozostaje nieruchoma: wygina się, pęka i przekierowuje swoje wewnętrzne siły. Jeśli tych przesunięć sił, czyli naprężeń, nie rozumie się poprawnie, może dojść do zawałów stropu, awarii filarów i nagłego uwolnienia gazów. W badaniu postawiono pozornie proste pytanie o duże implikacje dla bezpieczeństwa: na ile nasze komputerowe modele tych naprężeń zależą od tego, jak zarysujemy pudełko modelu i jak odwzorujemy puste, zasypane przestrzenie pozostawione po wydobyciu?

Wgląd w ułożone warstwy węgla

Naukowcy skupili się na chińskim rejonie wydobywczym, gdzie pięć udostępnialnych pokładów węgla leży stosunkowo blisko siebie. W miarę postępu wydobycia pustki powstające w tych warstwach — zwane zrobami — układają się jedna nad drugą, rozdzielone jedynie ławicami skalnymi o różnej wytrzymałości. Aby zbadać, jak zmieniają się naprężenia w takim układzie, zespół użył FLAC3D, powszechnie stosowanego programu symulacyjnego w inżynierii górniczej. Zbudowali dwa warianty podziemnego świata: wąski, cienki model tylko nieco szerszy niż pojedynczy panel ścianowy oraz duży model pełnej szerokości, sięgający znacznie dalej w bok. Następnie zasymulowali realistyczną sekwencję wybierania paneli w różnych pokładach, śledząc, jak ciężar nadkładu przekazywany jest na pozostały węgiel i skały w miarę tworzenia kolejnych pustek.

Jak rozmiar modelu zmienia obraz

Można by oczekiwać, że mniejszy model, z sztucznymi granicami bocznymi bliżej obszaru wydobycia, zniekształci pole naprężeń — i tak się dzieje, ale w specyficzny sposób. Model cienki ma tendencję do pokazywania silniejszego narastania naprężeń na krawędziach świeżo wydobytego panelu, szczególnie we wczesnych etapach, gdy wydobyto tylko jeden albo dwa pokłady. Ponieważ boki modelu nie mogą poruszać się swobodnie, zachowują się jak sztywne ściany, zmuszając naprężenia do koncentrowania się przy krawędziach zrobów. W większym modelu naprężenia rozkładają się łagodniej, a kontury wyglądają bardziej realistycznie. Jednak po wydobyciu trzech lub więcej pokładów różnica w wartościach szczytowych naprężeń między modelem cienkim a dużym maleje. Co kluczowe, oba modele lokują szczyty naprężeń prawie w tych samych miejscach wzdłuż pokładów węgla — rozmiar modelu zmienia głównie wysokość tych szczytów, a nie ich położenie.

Co naprawdę ma znaczenie we wnętrzu zrobów

Zdecydowanie większa różnica pojawiła się, gdy zespół zmienił sposób odwzorowania samego zrobu. W jednym wariancie zrob traktowano jako prawdziwą pustkę — tzw. model Null — nieoddającą żadnego oporu, więc naprężenia kumulują się głównie przy jego krawędziach. W drugim wariancie, nazwanym Double-Yield, zasypisko traktowano jako luźny, ale podatny na zagęszczenie materiał, który stopniowo przejmuje część obciążenia. W tym bardziej realistycznym ustawieniu naprężenia nie skupiają się tylko przy krawędziach zrobu; część obciążenia przejmowana jest przez zagęszczający się rumosz, a następnie przenoszona ku górze do nadkładu skalnego. W miarę eksploatacji kolejnych pokładów i nakładania się zrobów Double-Yield oddaje, jak naprężenia mogą narastać wewnątrz zasypanych stref i przechodzić przez nie, podczas gdy model Null pozostawia duże, nierealistyczne obszary niemal zerowego naprężenia. Wybór modelu zrobu znacząco przesuwa miejsca występowania szczytów naprężeń wzdłuż pokładów węgla — dużo bardziej niż jakakolwiek zmiana zewnętrznych rozmiarów siatki numerycznej.

Rola zapadającego się stropu

Badanie objęło także wpływ kąta, pod jakim strop zapada się do zrobu, na zachowanie naprężeń w modelu Double-Yield. Testując kilka kątów zapadania, autorzy stwierdzili, że bardziej strome i rozleglejsze zapadanie stropu prowadzi do silniejszego zagęszczenia rozkruszonej skały i lepszego kontaktu między fragmentami. W efekcie strefa zasypana przejmuje większą część ciężaru nadkładu, naprężenie wewnątrz zrobu rośnie, a główne koncentracje naprężeń przesuwają się ku górze, do warstw stropowych nad eksploatowanym panelem, zamiast pozostawać ściśle skupione przy krawędziach panelu. Takie zachowanie lepiej pokrywa się z obserwacjami terenowymi niż prostsze założenie pustki i podkreśla znaczenie kalibracji parametrów zrobu na podstawie rzeczywistych pomiarów odkształcalności zasypisk pod ziemią.

Co to oznacza dla bezpieczniejszej eksploatacji wielopokładowej

Mówiąc wprost, badanie pokazuje, że przy eksploatacji wielu warstw węgla poprawne odwzorowanie zrobu w modelu jest ważniejsze niż narysowanie bardzo dużego pudełka modelu. Wąski model nadal może przewidzieć, gdzie powstaną niebezpieczne punkty koncentracji naprężeń, pod warunkiem że inżynierowie uwzględnią, iż może on przeszacowywać intensywność tych szczytów. Jednak zastosowanie realistycznego, podatnego na zagęszczenie modelu zrobu — dostrojonego do warunków danego obiektu i kąta zapadania stropu — jest kluczowe, aby uchwycić, jak naprężenia przemieszczają się przez nałożone na siebie zasypane strefy i do pozostających pokładów węgla. Te wskazówki pomagają projektantom kopalń wybierać efektywne, a jednocześnie wiarygodne symulacje, poprawiając rozmieszczenie filarów, podpór i wyrobisk, tak by niewidzialne siły w głębokich złożach węgla były kontrolowane zanim przekształcą się w katastrofy.

Cytowanie: Wang, N., Paneiro, G.A., Li, Y. et al. Mechanisms of stress distribution influenced by numerical model size and goaf parameters in multi-coal seam mining. Sci Rep 16, 11137 (2026). https://doi.org/10.1038/s41598-026-42013-0

Słowa kluczowe: eksploatacja wielopokładowa węgla, zagęszczanie zrobów, modelowanie numeryczne, redystrybucja naprężeń, stabilność stropu kopalni