Materiały konferencyjne SEP 2018

Aby określić możliwości detekcji wstrząsów niskoenergetycznych dla analizowanej sieci sejsmologicznej dokonano analizy amplitud szumu sejsmicznego na kanałach pomiarowych w rejonie eksploatacji pokładu 209 ścianą 913 na wybranych losowo 108 cyfrowych sejsmo- gramach wstrząsów. Średnie wartości amplitud tych szumów dla poszczególnych kanałów przedstawiono w tabeli 4. Obliczono teoretyczne odległości detekcji wstrząsów niskoenergetycznych dla stanowisk pomiarowych, które przedstawiono w tabeli 5. Tabela 7. Estymowane maksymalne odległości [m] dla detekcji wstrząsów o założonym poziomie energetycznym Table 7. Estimated maximum distances [m] for tremors detection at the assumed energy level Poziom energii [J] Nr kanału/stanowiska 1 10 12 14 16 18 21 22 23 26 27 28 1  10 2 900 2250 1450 950 2700 2150 950 1400 3300 1800 1100 2100 2  10 2 1250 2950 1950 1300 3550 2850 1300 1900 4250 2400 1450 2800 3  10 2 1500 3450 2300 1550 4150 3300 1550 2250 4900 2850 1750 3250 4  10 2 1700 3850 2550 1750 4550 3700 1800 2500 5400 3200 2000 3650 5  10 2 1850 4200 2800 1900 4900 4000 1950 2750 5800 3450 2200 3950 1  10 3 2450 5300 2650 2550 6150 5100 2600 3600 7150 4450 2900 5000 Analiza danych wskazuje, że prowadzone obserwacje sejsmologiczne w rejonie ściany 913 charakteryzowały się błędami lokalizacji do 20m. Średni poziom szumów na poszczególnych kanałach pomiarowych wynosił 2,79  10 -7 [m/s]. Zastosowanie 64 kanałowej aparatury sejsmo- logicznej umożliwiło zwiększenie ilości stanowisk pomiarowych, a wykorzystanie mobilnych sond DLM2001 pozwoliło zmniejszyć odległości do obserwowanego rejonu. Dzięki wykorzy- staniu oddziałowych sieci obserwacyjnych wzrosły możliwości detekcji sieci w zakresie reje- stracji wstrząsów niskoenergetycznych, a zainstalowana sieć sejsmologiczna umożliwiała de- tekcję wstrząsów niskoenergetycznych od energii 1  10 2 J w sposób pozwalający na użycie tych zjawisk do obliczeń tomograficznych. 5. ANALIZA Analiza zmian pola prędkości została wykonana oprogramowaniem stworzonym i użyczonym przez prof. GIG dr. hab. inż. Adama Lurkę w wariancie tomografii pasywnej krzywoliniowej z wykorzystaniem algorytmu optymalizacji nieliniowej, którego dodatkową zaletą jest możli- wość uzyskania estymacji błędów wyznaczonych parametrów modelu [9]. Ze zbioru zareje- strowanych wstrząsów z rejonu ściany 913 w pokładzie 209 wyselekcjonowano zapisy, które pozwalały na przeprowadzenie obliczeń, a następnie określono dla nich czasy pierwszych wejść fali P. Pod uwagę brano zapisy, dla których można było określić co najmniej 5 czasów wejścia fali sejsmicznej . Do analizy wykorzystano 514 sejsmogramów wstrząsów zarejestro- wanych w okresie od 22.08.2013 do 2.10.2013 r. (5 cykli obliczeniowych). Obliczenia prze- prowadzono w 14 dniowych oknach czasowych z krokiem jednego tygodnia. Uzyskane, na podstawie zbiorów wstrząsów, rozkłady pola prędkości dla każdego cyklu obliczeniowego zostały skorelowane z aktywnością sejsmiczną w danym okresie i z sejsmicznością zarejestro-