Materiały konferencyjne SEP 2021

R – odległość hipocentralna, obliczona według zależności = √ 2 + ℎ 2 , gdzie r jest odległością epicentralną a h to średnia głębokość wstrząsów, odpowiadająca poziomowi eksploatacji w badanym rejonie, równa 950 m, S k – czynnik określający lokalne warunki geologiczne (k – numer stanowiska, k=1÷3), wartość ta wynosi 0 jeśli nie uwzględnia się na nim amplifikacji lub 1 jeśli estymujemy współczynnik amplifikacji. α 0 , α 1 , α 2 , α 3 , δ k – estymowane współczynniki relacji tłumienia Po przeprowadzeniu analizy regresji dokonuje się weryfikacji uzyskanego modelu. Sprawdza się między innymi znaki uzyskanych współczynników oraz ich zgodność z założeniami teoretycznymi, które posłużyły do budowy modelu. Ponieważ relacja tłumienia (1) opisuje zmienność szczytowych amplitud drgań sejsmicznych, współczynnik związany z logarytmem energii powinien być dodatni (amplituda drgań wzrasta ze wzrostem energii sejsmicznej), a współczynniki związane z odległością i jej logarytmem powinny być ujemne (amplituda drgań maleje ze wzrostem odległości). Za pomocą odpowiednich parametrów i testów statystycznych sprawdzana jest istotność uzyskanych współczynników. Do oceny jakości otrzymanych modeli służy analiza wartości resztowych, czyli różnic pomiędzy wartościami empirycznymi a teoretycznymi, wynikającymi z oszacowanej funkcji regresji. Reszty powinny mieć rozkład zgodny z rozkładem Gaussa. Nie jest to jednak zawsze konieczne i w przypadku dużej liczebności zebranych obserwacji wystarczającym jest, aby rozkład wartości resztowych nie odbiegał skrajnie od rozkładu normalnego a weryfikacja normalności rozkładu reszt może sprowadzić się do wizualnej analizy histogramu reszt lub wykresu normalności. Zebrane dane można także poddać testowi statystycznemu na normalność rozkładu. Rozkład reszt może być zaburzony przez tzw. wartości skrajne, związane z błędami pomiarowymi lub błędami przetwarzania. Dla poprawy modelu regresyjnego z danych usuwa się wartości skrajne i ponownie przeprowadza się analizę regresji. W celu poprawy dopasowania modelu do danych pomiarowych można usunąć do 20% obserwacji. Do analizy regresji danych sejsmometrycznych z obszaru KWK Sośnica wykorzystano model, w którym usunięto w obliczeniach współczynnik α 3 , związany z tłumieniem nieelastycznym z uwagi na niemający fizycznego sensu dodatni znak tego parametru obliczony według modelu (1). W związku z powyższym do dalszych obliczeń wykorzystano model opisany równaniem: = 0 + 1 + 2 + (2) gdzie przyjęto oznaczenia tak jak we wzorze (1) 4. ANALIZA STATYSTYCZNA RELACJI TŁUMIENIA Do estymacji współczynników regresyjnych dla relacji tłumienia prędkości PGV Hmax i przyspieszeń PGA 10Hz drgań wykorzystano odpowiednio 739 i 730 rejestracji. Współczynniki determinacji wyniosły 0,85 i 0,82, co oznacza, że przyjęte modele tłumienia prędkości i przyspieszeń drgań w 85 % i 82 % wyjaśniają zmienność danych pomiarowych. Estymatory błędów standardowych wynoszą odpowiednio 0,16 i 0,17. Testem Fishera sprawdzono czy obliczone modele regresyjne różnią się istotnie statystycznie od modelu o wartości stałej. Zweryfikowano hipotezę zerową H 0 , zakładając że wszystkie współczynniki regresji oprócz składowej stałej są zerem. Wyliczona wartość p dla tego testu zarówno dla prędkości i przyspieszeń drgań wyniosła mniej niż 10 -4 , co pozwala odrzucić hipotezę zerową.