The possibilities of using anisotropic models of rock mass to describe deformations of the mining area surface

• Autorzy: Wesołowski M.

Identyfikatory

DOI

10.1515/amsc-2016-0010

Warianty tytułu

PL

Możliwości zastosowania modeli anizotropowych górotworu do opisu deformacji powierzchni terenu górniczego

• Języki publikacji: EN

Abstrakty

EN

This paper presents the analysis of numerical modeling results of the influence of mining exploitation influence on the deformations of the area surface, with the use of an anisotropic elasto-plastic ubiquitous joint model and the transversely isotropic elastic model. The comparison of computer modeling results and geodetic measurements shows that with the use of transversely isotropic elastic model and ubiquitous joint model there is a possibility of the simultaneous description of both perpendicular and horizontal displacements of the area surface, caused by mining exploitation.

PL

Jak dotąd zasadniczym problemem związanym z modelowaniem numerycznym było zbudowanie modelu, który możliwie dokładnie opisywałby wszystkie etapy procesu deformacji górotworu. Niecki osiadania, uzyskiwane w wyniku modelowania numerycznego metodą elementów skończonych dla górotworu traktowanego jako ośrodek liniowo sprężysty izotropowy, znacznie różniły się od tych wyznaczanych z pomiarów geodezyjnych. Niecki te okazywały się zbyt rozległe i charakteryzowały się znacznie mniejszym (nawet 2,5-krotnie) nachyleniem zboczy w stosunku do niecek rzeczywistych. W celu wykazania, że stosowanie modeli numerycznych pozwala uzyskać w miarę dokładny opis rzeczywistych deformacji powierzchni terenu górniczego, przetestowano wiele modeli matematycznych, opisujących własności mechaniczne górotworu pod kątem zgodności z pomiarami geodezyjnymi. Testy te udowodniły, że w przypadku modelu liniowo sprężystego o transwersalnie izotropowej budowie warstw możliwe jest uzyskanie z obliczeń numerycznych w miarę dokładnego przybliżenia rzeczywistych deformacji powierzchni terenu górniczego (Tajduś, 2007; Białek et al., 2002; Mielimąka, 2009; Wesołowski et al., 2010). Pomimo poprawnego opisu deformacji powierzchni terenu stosowanie modelu transwersalnie izotropowego budzi wiele kontrowersji z uwagi na konieczność przyjęcia bardzo zróżnicowanych parametrów ośrodka w kierunku pionowym i poziomym. Pewne cechy podobieństwa do ośrodka transwersalnie izotropowego wykazuje model ubiquitous joint (Rys. 1). Odpowiednikiem kierunku prostopadłego do płaszczyzn izotropii modelu transwersalnie izotropowego może być tu kierunek prostopadły do płaszczyzn osłabienia. Model ubiquitous joint jest zatem anizotropowym modelem plastycznym zawierającym płaszczyzny osłabienia określonej orientacji. W modelu tym zaimplementowany został warunek wytrzymałościowy (uplastycznienia) Coulomba-Mohra. Do procesu modelowania numerycznego deformacji terenu górniczego, przeprowadzonego na potrzeby niniejszej pracy, wykorzystano program różnic skończonych FLAC (Itasca Consulting Group, Inc. 1992). Celem przeprowadzenia symulowanej komputerowej eksploatacji górniczej i określenia jej wpływu na deformacje powierzchni terenu zbudowano płaski model o wymiarach 2200 m × 913 m. Na głębokości 600 m (głębokość spągu pokładu) zamodelowano przeznaczony do eksploatacji pokład węgla o grubości 2 m. Schemat geometryczny modelu przedstawiono na Rys. 2. Wyniki symulacji komputerowej w zakresie opisu deformacji terenu górniczego porównane zostaną z przykładowymi pomiarami geodezyjnymi prowadzonymi na linii nr 100 podczas eksploatacji ścianowej pokładu 338/2 w KWK „Budryk” (Rys. 3). Parametry wytrzymałościowe oraz odkształceniowe warstw przyjęte zostały na podstawie literatury (Kidybiński, 1982; Prusek & Bock, 2008). Określając wartości parametrów płaszczyzn osłabienia posłużono się przypadkiem opisanym w pracy (Sainsbury et al., 2008). Zakres zmienności parametrów materiałowych warstw skalnych modelu ubiquitous joint oraz modelu transwersalnie izotropowego przyjęte do obliczeń przestawiono w tabelach 1 i 2. W pracy porównane zostały możliwości stosowania modelu transwersalnie izotropowego oraz modelu ubiquitous joint pod kątem zgodności opisu deformacji terenu górniczego wywołanych prowadzoną eksploatacją górniczą. W oparciu o wymienione powyżej modele górotworu przeprowadzona została symulacja komputerowa eksploatacji górniczej. Wyniki tych symulacji wykazały, że: 1. Przeprowadzony w ramach pracy cykl symulacji komputerowych wykazał, że zarówno dla modelu ubiquitous joint oraz modelu transwersalnie izotropowego istnieje możliwość jednoczesnego opisu zarówno pionowych, jak i poziomych ruchów górotworu, wywołanych eksploatacją górniczą. Uzyskanie bliskiego rzeczywistości opisu ruchów poziomych wymagało wprowadzenia płaszczyzn kontaktu (Interface) do modelowania połączeń międzywarstwowych. 2. Profil asymptotycznej niecki obniżeniowej ściany jest dla rozpatrywanych modeli asymetryczny względem wybranego pola. W profilach tych wartość maksymalnego nachylenia w rejonie krawędzi rozpoczynającej eksploatację jest nawet o kilkadziesiąt procent większa niż wartość maksymalnego nachylenia w rejonie krawędzi kończącej eksploatację. Podobne zależności dotyczą odkształceń poziomych. 3. Przedstawione wyniki symulacji komputerowych wskazują na to, że przy wykorzystaniu odpowiedniego ośrodka istnieje możliwość opisu kolejnych etapów deformacji terenu górniczego, w tym również wpływu kierunku prowadzenia eksploatacji na kształt profilu niecki obniżeniowej, kształtującej się nad postępującym frontem ścianowym. 4. Niewątpliwą zaletą modelu bazującego na ośrodku transwersalnie izotropowym jest stosunkowo mała ilość parametrów odkształceniowych koniecznych do obliczeń oraz możliwość łatwego dostosowania wyników obliczeń do wyników obserwacji geodezyjnych (Wesołowski, 2013). Uzupełnieniem prowadzonej analizy są rysunki 8 i 9 przedstawiające zasięg stref uplastycznienia rozpatrywanych układów modelu górotworu. Przeprowadzone w ramach pracy obliczenia komputerowe pokazały, że przy zastosowaniu numerycznych modeli górotworu opierającego się na sprężystym ośrodku transwersalnie izotropowym oraz anizotropowym modelu ubiquitous joint możliwy jest opis deformacji terenu górniczego jakościowo i ilościowo zgodny z obserwacjami geodezyjnymi.

Słowa kluczowe

EN

rock mass; calculation model; subsidence trough; area deformations

PL

górotwór; model obliczeniowy; niecka obniżeń; deformacje terenu

• Wydawca: Instytut Mechaniki Górotworu PAN
• Czasopismo: Archives of Mining Sciences
• Rocznik: 2016
• Tom: Vol. 61, no. 1
• Strony: 125--136
• Opis fizyczny: Bibliogr. 17 poz., rys., tab., wykr.

Twórcy

autor

Wesołowski M.

• Institute of Mining, Faculty of Mining and Geology, Silesian University of Technology, 44-100 Gliwice, ul. Akademicka 2, Poland

Bibliografia

• [1] Białek J., 2013. Wpływ własności reologicznych górotworu oraz kierunku eksploatacji na kształt nieustalonych niecek obniżeniowych. Przegląd Górniczy, nr 8, 8-13.
• [2] Białek J., Mielimąka R., Wesołowski M., 2002. Ein linear, transversal-anisotropisches Gebirgsmodell zur Modellierung Abbaubedingter Gebirgsbewegungen. Schriftenreihe des Institutes fűr Markscheidewesen und Geodäsie der TU Bergakademie Freiberg, Heft, 1, 184-191.
• [3] Białek J., Wesołowski M., 2011. Problematyka numerycznego modelowania ruchów terenu górniczego na przykładzie eksploatacji pokładu 354 w KWK ”Chwałowice”. Prace Naukowe GIG, XI Dni miernictwa górniczego i ochrony terenów górniczych, Kwartalnik, nr 2/1.
• [4] Budryk W., Knothe S., 1950. Wpływ eksploatacji podziemnej na powierzchnie z punktu widzenia ochrony obiektów. Przegląd Górniczy, nr 11.
• [5] FLAC User’s Manual 1992. Itasca Consulting Group. Minneapolis.
• [6] Kidybiński A., 1982. Podstawy geotechniki kopalnianej. Wydawnictwo Śląsk, Katowice.
• [7] Kołodziejczyk P., Wesołowski M., 2010. The influence of deformational parameters of a numerical model on the subsidence basin profile for chosen working dept. Arch. Min. Sci., Vol. 55, No 4, p. 775-781.
• [8] Mielimąka, R., 2009. Wpływ kolejności i kierunku eksploatacji prowadzonej frontami ścianowymi na deformacje terenu górniczego. Wydawnictwo Politechniki Śląskiej, Gliwice.
• [9] Popiołek E., Ostrowski J., 1980. Próba ustalenia głównych przyczyn rozbieżności prognozowanych i obserwowanych poeksploatacyjnych wskaźników deformacji. Ochrona Terenów Górniczych, nr 58, Katowice.
• [10] Prusek S., Bock S., 2008. Assessment of rock mass stresses and deformations around mine workings based on threedimensional numerical modelling. Arch. Min. Sci., Vol. 53, No 3, p. 349-360.
• [11] Sainsbury B., Pierce M., Mas Ivars D., 2008. Simulation of rock mass strength anisotropy and scale effects using a ubiquitous joint rock mass (UJRM) model. Continuum and Distinct Element Numerical Modeling in Geo-Engineering, Detournay & Cundall (eds.) Paper: 06-02 Itasca Consulting Group, Inc., Minneapolis.
• [12] Tajduś K., 2007. Numeryczne określanie metodą elementów skończonych, wpływu eksploatacji podziemnej na powierzchnię terenu. Przegląd Górniczy, nr 5, 36-42.
• [13] Tajduś K., 2008. Określenie wartości parametrów odkształceniowych górotworu uwarstwionego w rejonie wpływów eksploatacji górniczej. Praca doktorska (niepublikowana), Akademia Górniczo-Hutnicza, Wydział Górnictwa i Geoinżynierii, Kraków.
• [14] Tajduś K., 2010. Determination of approximate value of GSI index for the disturbed rock mass layers in the area of Polish mines. Arch. Min. Sci., Vol. 55, No 4, p. 879-890.
• [15] Wesołowski M., Białek J., Kołodziejczyk P., Plewa F., 2010. Modelowanie wpływów eksploatacji górniczej przy wykorzystaniu modeli numerycznych. Wydawnictwo Politechniki Śląskiej, Gliwice.
• [16] Wesołowski M., 2001. Wybrane aspekty modelowania numerycznego ruchów górotworu pod wpływem eksploatacji podziemnej i jej oddziaływania na obiekty. Praca doktorska (niepublikowana), Politechnika Śląska, Wydział Górnictwa i Geologii, Gliwice.
• [17] Wesołowski M., 2013. Zastosowanie liniowego ośrodka transwersalnie izotropowego do modelowania deformacji terenu górniczego. Wydawnictwo Politechniki Śląskiej, Gliwice.

Uwagi

PL

Opracowanie ze środków MNiSW w ramach umowy 812/P-DUN/2016 na działalność upowszechniającą naukę