Jak napisać pracę magisterską z sejsmologii

Sejsmologia to dziedzina, w której rzadko zbierasz dane sam. Najczęściej pracujesz na zapisach z sieci sejsmicznych, które ktoś rejestrował latami. Twoja praca magisterska to analiza tych zapisów: filtrowanie, identyfikacja faz, liczenie, interpretacja. I tu zaczyna się problem, bo dane są ogromne, a pytanie badawcze trzeba zawęzić do czegoś, co obronisz przed komisją.

Spotykam to często. Student pobiera kilka miesięcy zapisów z IRIS albo ORFEUS, otwiera ObsPy i nagle topi się w gigabajtach sygnału. Albo odwrotnie: ma świetnie policzone magnitudy, ale nie potrafi wyjaśnić, co z nich wynika dla budowy skorupy w badanym regionie.

Ten artykuł pokaże Ci, jak ułożyć pracę magisterską z sejsmologii tak, żeby trzymała się kupy. Od wyboru tematu, przez metodologię, po rozdziały i błędy, które najczęściej wyłapuję podczas korekty.

Jedna uwaga na początek. Sejsmologia to dziedzina ilościowa i komisja oczekuje liczb, nie opisu. Jeśli boisz się statystyki i programowania, ta praca Cię tego nauczy, bo bez tego się nie obronisz. Ale to też dobra wiadomość: gdy raz opanujesz pipeline w Pythonie, reszta idzie szybko, a Twoje wyniki są powtarzalne i odporne na pytania w stylu „a skąd ta liczba”.

Wybór tematu i pytania badawcze

Sejsmologia dzieli się na kilka nurtów i od tego, który wybierzesz, zależy cała reszta. Inna metodologia obowiązuje przy analizie ryzyka sejsmicznego, inna przy tomografii, jeszcze inna przy badaniu mikrosejsmiczności indukowanej przez górnictwo. Wybierz jeden i trzymaj się go.

Obszary, które dobrze sprawdzają się w pracach magisterskich z sejsmologii:

Lokalizacja epicentrum i hipocentrum – wyznaczanie współrzędnych ognisk trzęsień na podstawie czasów przyjścia fal P i S z kilku stacji, ocena błędu lokalizacji, porównanie z katalogami referencyjnymi
Analiza magnitudy i sejsmiczności – obliczanie magnitud (lokalnej ML, momentu Mw), wyznaczanie parametru b w relacji Gutenberga-Richtera, analiza rozkładu sejsmiczności w czasie i przestrzeni
Tomografia sejsmiczna – rekonstrukcja struktury prędkościowej skorupy lub płaszcza na podstawie czasów przebiegu fal, inwersja, interpretacja anomalii prędkości
Sejsmiczność indukowana i wywołana – trzęsienia związane z górnictwem (np. GZW, LGOM), eksploatacją złóż, zatłaczaniem płynów; korelacja z działalnością człowieka
Ocena zagrożenia i ryzyka sejsmicznego – analiza PSHA (probabilistic seismic hazard assessment), krzywe zagrożenia, przyspieszenia gruntu (PGA), zastosowania inżynierskie
Analiza fal powierzchniowych – dyspersja fal Rayleigha i Love’a, wyznaczanie profili prędkości fal poprzecznych, metoda HVSR do oceny warunków gruntowych

Pytanie badawcze sformułuj precyzyjnie i tak, żeby dało się je zweryfikować. Nie: „Jak wyglądają trzęsienia ziemi w Polsce?”, ale: „Czy wartość parametru b relacji Gutenberga-Richtera dla wstrząsów indukowanych w rejonie Legnicko-Głogowskiego Okręgu Miedziowego różni się istotnie od wartości typowych dla sejsmiczności tektonicznej?”

Jedno pytanie główne, do trzech szczegółowych. Reszta poczeka na doktorat.

Metodologia pracy magisterskiej z sejsmologii

Metodologia musi być tak opisana, żeby inny sejsmolog pobrał te same dane, puścił ten sam pipeline i dostał Twoje wyniki. To kręgosłup całej pracy. Komisja zna te wymagania i sprawdza je punkt po punkcie.

Źródła danych i charakterystyka stacji
Podajesz, skąd pobrałeś zapisy: konkretne sieci (np. IRIS DMC, ORFEUS, polska sieć PSN, EPOS), kody stacji, lokalizacje, typy sejsmometrów (szerokopasmowy, krótkookresowy, akcelerometr), częstotliwość próbkowania (np. 100 Hz), zakres dat. Format danych (miniSEED) i sposób dostępu (FDSNWS, protokół). Liczba stacji i ich rozmieszczenie geometryczne ma znaczenie dla błędu lokalizacji, więc to opisujesz.

Przetwarzanie wstępne sygnału
Usunięcie odpowiedzi instrumentu (deconvolution z funkcji przenoszenia), demeaning, detrending, filtracja (typ filtra, np. Butterworth, pasmo przepustowe w Hz, rząd filtra). Jeśli stosujesz okno czasowe wokół fazy, podajesz jego długość. Wszystko to jest powtarzalne tylko wtedy, gdy podasz parametry.

Identyfikacja i pikowanie faz
Jak wyznaczałeś czasy przyjścia fal P i S: ręcznie czy automatycznie (algorytmy STA/LTA, AR-AIC). Podajesz próg detekcji, długości okien STA i LTA. Przy pikowaniu ręcznym opisujesz kryteria. To jest najczęstsze źródło błędu w lokalizacji, więc komisja o to pyta.

Lokalizacja ognisk
Algorytm i program (HYPO71, HYPOINVERSE, NonLinLoc, czy moduł w ObsPy). Model prędkościowy (1D warstwowy, jaki, skąd zaczerpnięty, np. iasp91, ak135, lub model lokalny). Stosunek Vp/Vs. Sposób oceny niepewności lokalizacji (elipsa błędu, RMS rezyduów czasowych, gap azymutalny).

Wyznaczanie magnitudy
Typ magnitudy i wzór. Dla ML podajesz funkcję tłumienia i korekty stacyjne. Dla Mw opisujesz metodę: modelowanie widma przemieszczenia (model Bruna), wyznaczenie momentu sejsmicznego M0, częstotliwości narożnej. Podajesz, na ilu stacjach uśredniałeś i jaki był rozrzut między stacjami. Magnituda policzona z jednej stacji to wartość orientacyjna, nie wynik. Jeśli różne typy magnitud podajesz obok siebie, zaznaczasz wyraźnie, że to różne skale, i nie zestawiasz ich na jednym wykresie bez relacji przeliczeniowej.

Mechanizm ogniska (jeśli dotyczy)
Jeśli wyznaczasz mechanizm ogniska (fault plane solution), opisujesz metodę: polaryzacja pierwszych wejść fali P czy inwersja tensora momentu. Liczba stacji i rozkład azymutalny decydują o jednoznaczności rozwiązania. Podajesz program (np. FOCMEC, ISOLA) i sposób prezentacji (projekcja na dolną półkulę, beach ball). To pozwala powiązać sejsmiczność z konkretnym typem uskoku.

Analiza statystyczna i sejsmiczność
Relacja Gutenberga-Richtera: metoda wyznaczenia parametru b (najczęściej maximum likelihood wg Aki-Utsu, nie regresja liniowa), magnituda kompletności Mc i sposób jej wyznaczenia. Dla analiz czasowych: testy stacjonarności, identyfikacja roju lub sekwencji aftershocków (prawo Omoriego). Podajesz oprogramowanie i pakiety (ZMAP, ObsPy, Python z NumPy/SciPy, R).

Tomografia i inwersja (jeśli dotyczy)
Parametryzacja modelu (siatka węzłów, rozmiar komórek), metoda inwersji (LSQR, regularyzacja), test rozdzielczości (checkerboard test). Bez testu rozdzielczości wyniki tomograficzne są nie do obrony.

Struktura i rozdziały

Praca magisterska z sejsmologii ma zazwyczaj 60-90 stron. Typowy układ:

1. Wstęp – kontekst geologiczny i sejsmotektoniczny badanego regionu, przegląd dotychczasowych badań sejsmologicznych w tym obszarze, podstawy teoretyczne (rodzaje fal, mechanizm źródła) w zakresie niezbędnym dla pracy. Wskazujesz lukę: czego nie wiadomo i co Twoja praca dokłada.

2. Cel pracy i hipotezy – jasno, maksymalnie strona.

3. Podstawy teoretyczne – propagacja fal sejsmicznych, równanie falowe w niezbędnym zakresie, definicje magnitud, relacja Gutenberga-Richtera. Tyle teorii, ile potrzeba do zrozumienia metod, nie podręcznik.

4. Dane i obszar badań – charakterystyka sieci, stacji, okresu obserwacji, mapa rozmieszczenia stacji (rycina z podziałką, obowiązkowo).

5. Metodyka – szczegółowo, jak opisano wyżej.

6. Wyniki – katalogi, mapy epicentrów, wykresy magnitud, profile tomograficzne. Ryciny i tabele z opisami. Tekst opisuje liczby słownie, nie ogranicza się do „jak widać na ryc. 5”.

7. Dyskusja – zestawiasz wyniki z innymi badaniami i katalogami (EMSC, ISC, lokalne). Wyjaśniasz rozbieżności, oceniasz ograniczenia (geometria sieci, model prędkościowy, kompletność katalogu). Sugerujesz dalsze kierunki.

8. Wnioski – 5-8 punktów, każdy odpowiada na postawioną hipotezę.

9. Literatura – minimum 40-60 pozycji; czasopisma: Bulletin of the Seismological Society of America, Geophysical Journal International, Journal of Geophysical Research, Acta Geophysica.

Załączniki – katalog zdarzeń, sejsmogramy przykładowe, tabele rezyduów, mapy w wyższej rozdzielczości.

Potrzebujesz pomocy z pracą dyplomową?

Wyślij fragment tekstu, bezpłatną wycenę otrzymasz w 24 h. Poprawki bez limitu w cenie usługi.

Wyślij do wyceny

Najczęstsze błędy

Brak testu rozdzielczości w tomografii – kolorowe mapy anomalii prędkości bez checkerboard testu to obrazek, nie wynik. Komisja pyta o rozdzielczość zawsze
Parametr b liczony regresją liniową – to klasyk. Standardem jest metoda największej wiarygodności (Aki-Utsu), bo regresja przeszacowuje b przy małych próbach. Recenzent geofizyk to wychwyci
Pominięcie magnitudy kompletności – analiza Gutenberga-Richtera bez wyznaczenia Mc jest niewiarygodna, bo poniżej progu kompletności brakuje słabszych zdarzeń i statystyka się sypie
Brak usunięcia odpowiedzi instrumentu – liczenie magnitudy z surowego zapisu w zliczeniach (counts), bez przeliczenia na prędkość lub przemieszczenie, zafałszowuje amplitudy
Lokalizacja bez oceny błędu – podanie współrzędnych epicentrum bez RMS, elipsy błędu i gap azymutalnego to połowa informacji. Komisja chce wiedzieć, jak pewna jest ta lokalizacja
Mylenie intensywności i magnitudy – magnituda to miara energii ze źródła (jedna na zdarzenie), intensywność to skutki w danym miejscu (skala MSK, EMS-98, różna w różnych punktach). Mieszanie tych pojęć to błąd merytoryczny
Model prędkościowy bez uzasadnienia – „użyto modelu ak135” bez słowa, dlaczego akurat ten i czy pasuje do lokalnej geologii. Dla wstrząsów lokalnych model globalny często jest nieadekwatny
Niespójne jednostki magnitudy – raz ML, raz Mw, raz mb, bez zaznaczenia, że to różne skale, których nie wolno mieszać w jednym wykresie bez konwersji
Wykresy bez skali i orientacji – mapa epicentrów bez podziałki, strzałki północy i odwzorowania. Sejsmogram bez osi czasu i amplitudy
Dyskusja w oderwaniu od tektoniki – policzone zdarzenia bez powiązania z budową geologiczną i znanymi uskokami. Sejsmologia bez tła geologicznego wisi w próżni
Brak opisu pipeline’u przetwarzania – skrypt w Pythonie zrobił całą robotę, ale w pracy nie ma ani słowa o kolejnych krokach, parametrach filtra i wersjach pakietów. Skoro analiza jest komputerowa, opisujesz ją tak, żeby dało się ją odtworzyć. Recenzent może poprosić o powtórzenie

FAQ

Czy mogę napisać pracę magisterską z sejsmologii bez własnych pomiarów terenowych?

Tak, i tak powstaje większość prac. Sejsmologia opiera się na danych z otwartych sieci: IRIS DMC, ORFEUS, EIDA, EPOS udostępniają zapisy przez protokół FDSNWS za darmo. Analiza tych danych jest pełnoprawną pracą badawczą, pod warunkiem że pytanie badawcze i metoda do tego pasują. Opisujesz wtedy dokładnie pochodzenie danych w rozdziale metodycznym. Miałam pod korektą prace oparte wyłącznie na danych z repozytoriów i broniły się bez zarzutu.

Jakiego oprogramowania użyć do analizy sejsmologicznej?

Standardem akademickim jest ObsPy, czyli biblioteka Pythona do pracy z danymi sejsmicznymi: pobieranie, filtracja, pikowanie, podstawowe analizy. Do lokalizacji: HYPOINVERSE, NonLinLoc, HYPO71. Do analizy sejsmiczności: ZMAP. Do tomografii: pakiety inwersyjne zależne od metody. Wszystkie są darmowe. Podaj wersję każdego programu w metodyce.

Ile zdarzeń sejsmicznych potrzebuję, żeby analiza była wiarygodna?

To zależy od pytania. Do wyznaczenia parametru b relacji Gutenberga-Richtera metodą największej wiarygodności mówi się o minimum 50-100 zdarzeniach powyżej magnitudy kompletności, im więcej tym węższy przedział ufności. Do lokalizacji pojedynczego ogniska potrzebujesz zapisów z co najmniej 4 stacji, ale im więcej i im lepsza geometria, tym mniejszy błąd. Ustal minimalne n z promotorem zanim zaczniesz liczyć.

Czy potrzebuję silnej matematyki, żeby zrobić tomografię sejsmiczną?

Tomografia to problem odwrotny i wymaga zrozumienia podstaw algebry liniowej oraz regularyzacji. Nie musisz wyprowadzać wszystkiego od zera, bo używasz gotowych pakietów inwersyjnych. Musisz jednak rozumieć, co robi inwersja, dlaczego potrzebna jest regularyzacja i jak interpretować test rozdzielczości. Bez tego nie obronisz wyników. Jeśli matematyka inwersji Cię przerasta, rozważ temat z lokalizacji lub analizy sejsmiczności, gdzie aparat jest prostszy.

Podsumowanie

Praca magisterska z sejsmologii broni się szczegółami metodycznymi. Dane masz z sieci, narzędzia są darmowe, ale to, jak je przetworzyłeś i jak ocenisz niepewność wyników, decyduje o wartości naukowej. Najsłabszym punktem takich prac bywa dyskusja oderwana od geologii i wyniki podane bez oceny błędu. Zadbaj o jedno i drugie.

Jeśli zbliżasz się do obrony i chcesz mieć pewność, że praca jest dopięta – zapraszam na dobrzenapisane.pl. Koretuję prace z nauk o Ziemi regularnie: sprawdzam język, spójność terminologii (magnituda kontra intensywność, jednostki, oznaczenia faz), podpisy pod rycinami i opisy metod. Napisz, odpiszę w ciągu 24 godzin.

Dowiedz się więcej

Stein, S., Wysession, M. (2003). An Introduction to Seismology, Earthquakes, and Earth Structure. Blackwell Publishing
Lay, T., Wallace, T.C. (1995). Modern Global Seismology. Academic Press
Shearer, P.M. (2009). Introduction to Seismology. Cambridge University Press
IRIS Data Management Center: https://www.iris.edu
ORFEUS – Observatories and Research Facilities for European Seismology: https://www.orfeus-eu.org

Inne artykuły na blogu:
– Jak znaleźć literaturę do pracy dyplomowej
– Metodologia pracy magisterskiej – od czego zacząć
– Korekta pracy magisterskiej – co obejmuje i ile kosztuje