Jak działa sejsmograf i co można wyczytać z wykresu trzęsienia ziemi?

Dlaczego wykres z sejsmografu bywa lepszą „relacją na żywo” niż relacje w mediach

Gdy na lokalnej grupie mieszkańców mnożą się wpisy „zatrzęsło?”, ktoś zwykle wrzuca zrzut sejsmogramu z najbliższej stacji i sprawa się klaruje w kilka minut. Wykres trzęsienia ziemi jest bezstronny: rejestruje wstrząsy, gdy inne źródła jeszcze milczą. Problem zaczyna się, gdy próbujemy samodzielnie zrozumieć, co oznaczają pierwsze ząbki, szerokie wahania i nagłe „wyjście poza skalę”. Dobrze ustawiony sejsmograf pokaże dużo: moment nadejścia fal P i S, odległość do epicentrum, siłę wstrząsu i różnicę między trzęsieniem a eksplozją. Źle zinterpretowany wykres podsunie błędny wniosek, np. że „mieliśmy duże trzęsienie”, gdy instrument się tylko nasycił lub zadziałał hałas miejski.

Poniżej przewodnik: jak działa sejsmograf i co – krok po kroku – da się odczytać z wykresu trzęsienia ziemi, a także kiedy takie wnioski są wiarygodne, a kiedy lepiej się wstrzymać.

Krótkie pytania użytkownika i decyzje, które pomogą uniknąć błędów

Najczęstsze pytania praktyczne

• Czy to, co widzę na sejsmogramie, to faktycznie trzęsienie ziemi, czy tylko hałas?
• Jak rozpoznać na wykresie fale P i S oraz fale powierzchniowe?
• Czy da się z jednego wykresu oszacować odległość do epicentrum i siłę wstrząsu?
• Kiedy amplituda na wykresie przekłada się na magnitudę, a kiedy nie wolno tego robić?
• Jakie filtry i skale czasowe włączyć, aby nie „zgubić” kluczowych cech zdarzenia?
• Jak odróżnić wstrząs górniczy lub eksplozję od naturalnego trzęsienia?
• Jaki typ sejsmografu ma sens: czuły szerokopasmowy czy akcelerograf odporny na silne drgania?

Decyzyjna lista kontrolna przed interpretacją wykresu

• Sprawdź zegar: czy czas na wykresie jest zsynchronizowany (UTC lub lokalny, bez przesunięć)?
• Sprawdź kanał: pionowy (Z) pomoże złapać fale P, poziome (N/E lub X/Y) dobrze pokażą fale S i powierzchniowe.
• Oceń tło hałasu: czy to dzień w mieście czy noc w cichym miejscu? Jaki jest poziom szumu poza zdarzeniem?
• Wyklucz nasycenie: płaskie „ścięte” wierzchołki sygnału oznaczają clipping i dyskwalifikują amplitudową analizę.
• Zidentyfikuj P i S: różnica czasów S–P to wstępna miara odległości (z zastrzeżeniami regionalnymi).
• Dobierz filtr: dla fal powierzchniowych zwykle pasmo niższych częstotliwości (np. 0,02–0,2 Hz), dla lokalnych P i S – wyższe (np. 1–10 Hz).
• Podeprzyj się więcej niż jednym kanałem i – jeśli to możliwe – więcej niż jedną stacją.

Sejsmograf w praktyce: z czego się składa i co naprawdę mierzy

Rdzeń urządzenia: czujnik, który „stoi w miejscu”

Sejsmograf opiera się na zasadzie bezwładności: masa sejsmetryczna chce pozostać w spoczynku, gdy grunt się porusza. Urządzenie mierzy różnicę ruchu między masą a obudową związaną z podłożem. W nowoczesnych stacjach masa ta jest sprzężona elektronicznie, co poszerza pasmo i stabilizuje wskazania. Zależnie od konstrukcji rejestrujemy przemieszczenie, prędkość cząsteczkową lub przyspieszenie. Informacja ta trafia do przetwornika A/C i dalej do rejestratora.

Typy czujników: szerokopasmowy kontra akcelerograf

Co do zasady rozróżniamy dwie grupy instrumentów:

• Sejsmometry szerokopasmowe – czułe na bardzo małe wibracje i bardzo szerokie pasmo (od fal długookresowych po lokalne wysokie częstotliwości). Świetne do trzęsień regionalnych i odległych (telesejsmicznych), słabiej znoszą nasycenie przy bardzo silnych ruchach blisko epicentrum.
• Akcelerografy – zaprojektowane do rejestracji silnych drgań bez nasycania; mierzą przyspieszenie. Stosowane w inżynierii trzęsień ziemi i przy silnych wstrząsach lokalnych.

Tor pomiarowy: od ruchu do pliku

Kluczowe elementy toru to: czujnik, wstępne wzmocnienie, filtr antyaliasingowy, przetwornik ADC (zwykle 24-bit), zegar (GPS lub inny system synchronizacji), rejestrator danych i transmisja do centrum. Każdy z tych komponentów może zafałszować odczyty, jeśli jest źle dobrany lub skalibrowany. Stąd w profesjonalnych stacjach dba się o codzienną kontrolę opóźnienia zegara, szumy własne i odpowiedź częstotliwościową instrumentu.

Jak działa sejsmograf: od bezwładności masy do zrozumiałego wykresu

Zasada bezwładności i odpowiedź instrumentu

Masa sejsmometru jest mechanicznie odsprzęgnięta od obudowy. Gdy grunt przyspiesza, na wykresie widać proporcjonalny sygnał: w zależności od konstrukcji mierzona jest prędkość, przemieszczenie lub przyspieszenie. Odpowiedź instrumentu (tzw. response) mówi, jak urządzenie przenosi różne częstotliwości. Fachowo rekonstruuje się „prawdziwy” ruch gruntu przez dekonwolucję odpowiedzi instrumentu. Bez tego porównywanie amplitud z różnych stacji bywa złudne.

Skala czasu i filtracja

Sejsmogram jest funkcją czasu. Skalę dobiera się do zdarzenia: kilkanaście sekund dla lokalnych wstrząsów, godziny dla telesejsmów. Filtry pasmowo-przepustowe eksponują interesujące pasma, ale mogą zniekształcać początki faz (tzw. efekt przebiegu filtra). Bezpieczną procedurą jest analiza surowego wykresu, a dopiero potem nakładanie łagodnych filtrów (np. Butterwortha 2–4 rzędu) i sprawdzanie stabilności wniosków.

Zegar i synchronizacja

Różnice w oznaczeniu czasu rzędu sekund decydują o precyzyjnym wyznaczeniu odległości do epicentrum przez różnice S–P. W praktyce używa się zegarów synchronizowanych GPS lub NTP. Jeśli stacja ma dryf zegara, wnioski o lokalizacji zdarzenia mogą być istotnie błędne.

Jak czytać wykres trzęsienia ziemi: fale P, S i powierzchniowe

Identyfikacja faz: od czego zacząć

Standardowo kanały oznacza się jako Z (pion), N (północ–południe), E (wschód–zachód). Krok po kroku:

1. Szukanie pierwszej energii na kanale Z – to zwykle fale P (podłużne); nadchodzą jako pierwsze, mają wyższe częstotliwości i mniejszą amplitudę.
2. Po czasie Δt pojawiają się fale S (poprzeczne), wyraźniejsze na kanałach poziomych, o większej amplitudzie i nieco niższej częstotliwości.
3. Jeszcze później widać fale powierzchniowe (Rayleigha i Love’a), z największymi amplitudami i najniższymi częstotliwościami, dominujące w odległych i silnych zdarzeniach.

Różnica czasów S–P i odległość do epicentrum

Różnica czasu nadejścia fal S i P – na tej samej stacji – pozwala oszacować odległość do epicentrum. Przy założeniach prędkości w skorupie (P ~ 6 km/s, S ~ 3,5 km/s) odległość D w kilometrach szacuje się orientacyjnie jako 8–9 razy różnica S–P w sekundach. Przykładowo, jeśli S–P = 10 s, odległość wynosi około 80–90 km. To przybliżenie; struktura ośrodka, lokalna geologia i głębokość ogniska wprowadzają odchylenia. Precyzyjna lokalizacja wymaga trzech lub więcej stacji i rozwiązania triangulacyjnego.

Polaryzacja pierwszego ruchu

Pierwsze wychylenie (do góry lub do dołu) bywa używane do oceny mechanizmu źródła. Zapis na kilku stacjach pozwala odróżnić ściskanie od rozciągania i zrekonstruować płaszczyzny uskoków (tzw. „beachball”). Dla amatorskiej interpretacji pojedynczej stacji to subtelna wskazówka, a nie dowód.

Co można wyczytać z sejsmogramu i gdzie kończą się możliwości

Informacje, które uzyskasz z jednej stacji

• Przybliżony czas zdarzenia w UTC/lokalnym i czasy nadejścia fal P, S.
• Szacunkową odległość do epicentrum z różnicy S–P.

Granice interpretacji amplitudy i skali

Duże wychylenie na wykresie nie oznacza automatycznie „silnego trzęsienia”. Co do zasady amplituda na sejsmogramie jest mieszanką trzech czynników: energii źródła, tłumienia na drodze (geometria + pochłanianie) oraz odpowiedzi instrumentu. Decyzja „silne/niezwykłe” ma sens dopiero po sprawdzeniu, że:

• instrument nie nasycił się (brak „ściętych” wierzchołków i nagłych płaskich plateau),
• zastosowana skala czasu nie wyolbrzymia pojedynczych cykli (np. zbyt duże powiększenie sekund),
• odpowiedź częstotliwościowa czujnika została uwzględniona lub dane zostały przeliczone do standardu (np. symulacja Wood–Anderson do ML),
• ścieżka propagacji nie jest wybitnie nietypowa (np. przejście przez basen osadowy silnie wzmacniający niskie częstotliwości).

Zakłócenia, które udają trzęsienie: co wykluczyć najpierw

W praktyce bywa tak, że krótkie, kontrastowe impulsy pochodzą z bliskiego otoczenia stacji. Zanim uznasz zdarzenie za sejsmiczne, wyklucz:

• prace budowlane i ruch ciężkich pojazdów – częste w godzinach dziennych, szczególnie na kanałach poziomych,
• wiatr i efekt przewiewania kabli/masztów – powolne, nieregularne wychylenia o niskich częstotliwościach,
• uderzenia drzwi, spadki przedmiotów w budynku – krótkie impulsy bez wyraźnej sekwencji P–S–powierzchniowe,
• mikrosejsmy oceanu – szerokie pasmo niskich częstotliwości (ok. 0,1–0,3 Hz) bez wyraźnych faz, rosną przy sztormie,
• przemagnesowanie lub błąd zasilania rejestratora – synchroniczne „skoki” na wszystkich kanałach bez struktury falowej.

Eksplozja a trzęsienie: kryteria wstępnego rozróżnienia

Odróżnienie nie zawsze jest możliwe z jednej stacji, ale istnieją wskaźniki, które pomagają wstępnie ocenić sytuację:

• Stosunek P/S – dla eksplozji bywa wyższy (silna P, słaba S), dla trzęsienia stosunek jest niższy i bardziej zrównoważony.
• Długość ogona (coda) – trzęsienia mają zwykle dłuższy coda decay, eksplozje krótszy i bardziej „impulsowy”.
• Głębokość i pora – płytkie, powtarzalne, w tych samych godzinach pracy kopalni/poligonu sugerują zdarzenia antropogeniczne.
• Widmo – eksplozje wykazują często bardziej wąskie pasmo i ostrzejszy początek; trzęsienia – bogatsze widmowo, z rozwiniętymi falami powierzchniowymi.

Decyzja „eksplozja” powinna być formułowana ostrożnie i raczej jako hipoteza, dopóki nie potwierdzą jej inne stacje lub źródła.

Magnituda z pojedynczej stacji: kiedy oszacujesz, a kiedy lepiej poczekać

Przeliczenie amplitudy na magnitudę bywa kuszące, ale ma sens tylko w określonych warunkach.

• Kiedy ma sens:
  • dane są skalibrowane i znasz odpowiedź instrumentu,
  • zdarzenie jest lokalne/ regionalne i możesz zasymulować odpowiedź Wood–Anderson (ML),
  • odległość do epicentrum jest oszacowana (np. z S–P) i mieszcząca się w zakresie, dla którego obowiązują korekty tłumienia,
  • brak nasycenia i brak silnych efektów lokalnych (np. „site amplification”), albo potrafisz je skorygować.
• Kiedy nie ma sensu:
  • telesejsmy i zdarzenia bardzo odległe – jedna stacja bez modeli propagacji nie da wiarygodnej Mw,
  • silne wstrząsy w pobliżu – sejsmometr szerokopasmowy łatwo się nasyci; lepszy jest akcelerograf,
  • brak pewnego czasu i odległości – korekty geometrii i tłumienia będą przypadkowe,
  • jedyny dostępny kanał to mocno szumiący lub filtrowany wąskopasmowo zapis.

Dobór instrumentu do zadania: szerokopasmowy czy akcelerograf

Wybór czujnika to decyzja „kiedy tak / kiedy nie” w czystej postaci.

• Sejsmometr szerokopasmowy:
  • Wybierz, gdy: chcesz rejestrować telesejsmy i regionalne zdarzenia, zależy Ci na widmie od okresów długich do kilku–kilkunastu Hz, lokal jest spokojny sejsmicznie (niski poziom hałasu).
  • Uważaj, gdy: spodziewasz się silnych drgań w pobliżu – grozi clipping; rozważ parę z akcelerografem.
• Akcelerograf:
  • Wybierz, gdy: priorytetem są silne, bliskie wstrząsy i analiza przyspieszeń dla inżynierii trzęsień; instalacja w budynku/infrastrukturze.
  • Uważaj, gdy: interesują Cię słabe, dalekie sygnały – czułość może być niewystarczająca.

Szybka ścieżka decyzji: publikować wniosek czy poczekać?

Kilka prostych kroków porządkuje działanie i ogranicza ryzyko błędu.

1. Czas i synchronizacja: potwierdź poprawność zegara i spójność strefy czasowej.
2. Fazy: wskaż P na kanale Z i S na N/E; zapisz różnicę S–P.
3. Jedna kontra wiele stacji: sprawdź najbliższe stacje (jeśli dostępne). Brak korelacji – wstrzymaj się.
4. Nasycenie i filtry: odrzuć zapisy ze „ścięciem”; używaj łagodnych filtrów i sprawdź, czy wnioski nie zmieniają się po zmianie pasma.
5. Źródło antropogeniczne: porównaj z godzinami pracy kopalń/budów; powtarzalność w tym samym miejscu i porze to żółta flaga.
6. Magnituda: licz tylko wtedy, gdy spełniasz warunki z sekcji o ML; w pozostałych przypadkach odwołaj się do katalogów instytucji sejsmologicznych.

Jeśli po punktach 1–4 masz spójną sekwencję P–S–powierzchniowe, podobny sygnał na co najmniej dwóch stacjach i brak oznak nasycenia – komunikat „zarejestrowano wstrząs, odległość szacunkowa ~X km” jest co do zasady bezpieczny. Gdy którykolwiek warunek odpada, lepiej ograniczyć się do opisu: „zarejestrowano anomalię, trwają weryfikacje”.

Krótka praktyka: dwa scenariusze

Przykład 1: wieczór, stacja poza miastem. Na Z pojawia się wyraźne P, po 7 s na N/E – S, po kolejnych kilkunastu sekundach długookresowe fale powierzchniowe. Brak clippingu. S–P ~7 s → odległość ok. 60 km. Druga stacja 40 km dalej widzi podobny przebieg z przesunięciem. Wniosek: lokalne trzęsienie, ostrożna lokalizacja możliwa.

Przykład 2: południe, stacja w budynku. Krótki, ostry impuls bez wyraźnej sekwencji P–S, silniejszy na jednym kierunku poziomym, brak potwierdzenia na sąsiedniej stacji. Po sprawdzeniu – wózek widłowy na rampie obok. Wniosek: nie publikować tezy o wstrząsie.

Filtracja i prezentacja zapisu: decyzje, które zmieniają wniosek

Dobór filtra i skali wykresu potrafi zadecydować, czy zobaczysz czytelne P i S, czy jedynie „szum”. Zanim zmienisz pasmo, odpowiedz na dwa pytania: jakiej odległości się spodziewasz i czy chcesz wydobyć fazy, czy ocenić energię w określonym paśmie.

W praktyce zacznij od łagodnego pasma (np. 0,5–15 Hz dla regionu), zaznacz czasy P i S, a potem sprawdź, czy wnioski utrzymują się po zmianie filtra. Jeśli nie – decyzję odłóż do weryfikacji na surowych danych.

Parametry wykresu: skala czasu i zakres amplitudy

• Skala czasu:
  • Ustaw gęstszy podgląd (np. 30–60 s) do pickingu P/S;
  • rozszerz okno (5–10 min) przy ocenie fal powierzchniowych i coda. Przeskalowanie powinno być jawne – opisz zakres na osi.
• Zakres amplitudy:
  • Zaczynaj od auto‑scale, ale zawsze sprawdź manualnie, czy nie „ucinasz” szczytów lub nie chowasz cichych faz;
  • w razie wątpliwości porównaj z drugim kanałem (Z vs N/E) – niespójności sygnalizują nasycenie lub efekt lokalny.
• Opis toru:
  • Zaznacz filtr i jednostki (np. m/s, m/s², znormalizowane) bezpośrednio na wykresie; bez tego odbiorca nie porówna amplitud ani czasu narastania.

Przykład z praktyki: rekord „cichy” w autoskali po włączeniu high‑pass 1 Hz nagle pokazuje wyraźne P i S z S–P ~5 s. Ten sam zapis w szerokim paśmie wyglądał jak mikrosejsmy. Decyzja o filtrze przesądziła o rozróżnieniu hałasu i zdarzenia lokalnego.

Triangulacja kilku stacji nieprofesjonalnych: kiedy próbować, a kiedy zrezygnować

Amatorska sieć potrafi dać poprawne położenie, ale tylko przy spełnieniu kilku warunków.

• Kiedy próbować:
  • masz co najmniej trzy stacje w promieniu do ~200 km, z zegarem GNSS lub znaną poprawką czasu,
  • sygnał P/S jest czytelny na osi Z i przynajmniej jednym kanale poziomym,
  • geometria nie jest skrajnie niekorzystna (stacje nie leżą prawie w linii prostej).
• Kiedy odpuścić:
  • różnice czasów między stacjami są porównywalne z niepewnością pickingu i synchronizacji,
  • obszar ma silne kontrasty prędkości (baseny osadowe obok masywów krystalicznych), a brak lokalnego modelu prędkości,
  • sygnał jest ewidentnie zdominowany przez hałas miejski lub wiatr – lepiej poczekać na katalog instytucjonalny.

Prosty workflow: wyznacz czasy P (i S, jeśli możliwe) na każdej stacji → z S–P oszacuj odległości → narysuj okręgi wokół stacji → sprawdź punkt/przecięcie pasów niepewności. Gdy okręgi się „rozjeżdżają” – przyczyną bywa pojedynczy błędny pick lub dryf zegara jednej stacji; w pierwszej kolejności weryfikuj te dwie rzeczy.

Korelacja kształtu fal: szybkie potwierdzenie, że to „to samo” zdarzenie

Dla stacji oddalonych o podobną odległość od epicentrum, wykresy po wyrównaniu na P powinny być podobne w pierwszych sekundach. Jeśli po filtracji w tym samym paśmie korelacja jest niska, a różni się nawet polaryzacja pierwszego ruchu – prawdopodobnie porównujesz inne zdarzenia lub zakłócenie z lokalnym hałasem.

Instalacja i serwis czujnika: decyzje terenowe, które „robią” jakość

• Lokalizacja:
  • Wybierz podłoże stabilne (skała, masywny beton) – zyskujesz lepszy stosunek sygnału do szumu.
  • Unikaj wnętrz budynków z ruchem ludzi/urządzeń oraz poddaszy – drgania konstrukcyjne dominują pasmo 1–20 Hz.
• Termika i wiatr:
  • Stabilizuj temperaturę (obudowa, zakopanie, izolacja) i eliminuj przewiewy – ograniczysz dryf i niskoczęstotliwościowe zakłócenia.
• Zasilanie i zegar:
  • Zapewnij podtrzymanie i synchronizację GNSS; bez pewnego czasu interpretacja faz traci sens decyzyjny.
• Kable i mocowanie:
  • Unieruchom kable, używaj przepustów z odciążeniem – „tańczące” przewody generują fałszywe niskie częstotliwości.

Krótki przykład: przeniesienie czujnika z piwnicy biurowca na płytę betonową poza budynkiem obniżyło tło w 2–8 Hz o kilka dB – faza S stała się rozróżnialna, co wcześniej było niemożliwe mimo tych samych filtrów.

Lista kontrolna przed udostępnieniem wykresu

• Na osi opisane: jednostki, filtr, okno czasowe i kanał(e).
• Brak oznak nasycenia i brak „płaściutkich” segmentów.
• Te same wnioski utrzymują się w co najmniej dwóch rozsądnych pasmach filtracji.
• Porównanie z sąsiednią stacją nie przeczy hipotezie (choćby z niższą rozdzielczością).

Najczęstszy błąd interpretacyjny i prosta kontra

Mylenie dużej amplitudy po filtracji z „silnym trzęsieniem”. Kontra jest prosta: zawsze zestaw wykres przefiltrowany z surowym, wskaż rodzaj filtra i sprawdź, czy relacja P/S i obecność fal powierzchniowych nie znikają po zmianie pasma. Jeśli znikają – decyzję odłóż, bo filtr mógł stworzyć pozór jednoznaczności.

Magnituda z pojedynczej stacji: kiedy wyliczać, a kiedy wstrzymać się z liczbami

Szacunek magnitudy kusi prostotą, ale bez kilku warunków łatwo o wynik pozornie precyzyjny. Kluczowe są: kalibracja toru (znana odpowiedź instrumentu), brak nasycenia oraz dobór właściwej definicji magnitudy do dystansu i pasma sygnału.

• Kiedy wyliczać:
  • tor jest skalibrowany (masz stałą czułości lub bieguny/zera) i pracujesz w jednostkach fizycznych po usunięciu odpowiedzi instrumentu,
  • brak clippingu na interesujących fragmentach P/S oraz stabilna linia bazowa,
  • zdarzenie lokalne–regionalne (kilka–kilkaset km), a sygnał mieści się w pasmach typowych dla ML/Md.
• Kiedy odpuścić:
  • nieznane wzmocnienie lub orientacja kanałów,
  • silny wpływ miejsca (mocne rezonanse 1–10 Hz bez możliwości korekty) albo bardzo nierówny szum tła,
  • telesejsm i wąskie pasmo – lokalne formuły ML przestaną działać; bez analizy widmowej i wielu stacji nie wyciągniesz wiarygodnego Mw.

W praktyce trzy podejścia mają sens:

• ML (lokalna): mierz amplitudę maksymalną na poziomych komponentach po przekształceniu do odpowiedzi Wood–Anderson i zastosuj tłumienie odległościowe właściwe dla regionu. Daje rozsądne wyniki dla odległości do kilkuset km, jeśli tor nie był nasycony.
• Md (czas trwania): licz czas zaniku cody powyżej progu. Przydaje się, gdy szczyty są spłaszczone lub szum wysoki. Wymaga stałego progu i tej samej filtracji dla porównań.
• PGV/PGA jako wskaźniki skutków: podaj prędkość/akcelerację szczytową wraz z pasmem i oknem. Można je wstępnie mapować na intensywność odczuć, ale tylko orientacyjnie i z zastrzeżeniem wpływu podłoża.

Przykład z praktyki: lokalne zdarzenie, S–P ~4 s, brak clippingu. Po dekonwolucji do prędkości i symulacji Wood–Anderson amplituda maksymalna na poziomych kanałach jest spójna; ML z dwóch stacji różni się o 0,2–0,3. Wniosek: można opublikować przedział magnitudy zamiast jednej cyfry.

Usuwanie odpowiedzi instrumentu i jednostki: decyzje przed pomiarem amplitudy

• Kiedy stosować dekonwolucję:
  • gdy chcesz porównywać amplitudy między stacjami lub wyliczyć ML/Md/PGV/PGA,
  • gdy potrzebujesz przemieszczenia (integracja prędkości) do oceny niskoczęstotliwościowych fal powierzchniowych.
• Jak to zrobić bezpiecznie:
  • użyj filtru wstępnego (pre‑filter), aby nie wzmacniać skrajów pasma,
  • stosuj filtry zero‑phase; podaj na wykresie docelową wielkość fizyczną (m/s, m/s², m) i pasmo,
  • sprawdź spójność kanałów: PGV/ML z N i E powinny być porównywalne (różnice rzędu dziesiątek procent są typowe, wielokrotność – to sygnał problemu).
• Kiedy zostać przy „counts”:
  • do szybkiej oceny obecności faz i względnych różnic między zdarzeniami na tej samej stacji,
  • gdy nie masz pewnych parametrów instrumentu – wtedy nie publikuj wartości fizycznych ani magnitudy.

Głębokość i charakter źródła: ile da się powiedzieć z kształtu fal

Różnica S–P daje głównie odległość epicentralną. Wnioski o głębokości wymagają ostrożności – bez modelu prędkości łatwo o błąd systematyczny.

• Wstępne przesłanki „płytkiego” źródła:
  • silne, długotrwałe fale powierzchniowe wyraźnie po fazach objętościowych,
  • duży stosunek amplitudy Rg/Rayleigh do S na dystansach lokalnych.
• Wskazówki „głębszego” ogniska:
  • wyraźne, twardsze P z mniejszą zawartością niskich częstotliwości i słabsze fale powierzchniowe,
  • dla telesejsm – obecność odbić pP/sP oddzielonych o kilka–kilkanaście sekund (wymaga dobrego SNR i filtracji niskoczęstotliwościowej).
• Gdy lepiej nie wnioskować:
  • silne efekty miejsca (basen osadowy, obiekty inżynieryjne) dominują zapis,
  • pojedyncza stacja i brak jednoznacznych fal powierzchniowych – opisz raczej „lokalne/telesejsmiczne” niż „płytkie/głębokie”.

Krótki przykład: dwa lokalne zdarzenia o podobnym S–P. W pierwszym fale powierzchniowe „niosą się” przez ponad minutę i dominują amplitudą – najpewniej płytkie. W drugim zapis kończy się szybko po S, a Rayleigh jest słaby – dopuszczalna hipoteza o większej głębokości, ale bez modelu prędkości pozostaje jako przypuszczenie.

Zakres deklaracji publicznej: co można powiedzieć „na pewno”, a co tylko warunkowo

• Można podać wprost:
  • czas przyjścia P/S na konkretnej stacji (z niepewnością pickingu),
  • oszacowanie odległości z S–P z podanym zakresem (np. ±20%),
  • obecność/nieobecność fal powierzchniowych w zadanym paśmie,
  • PGV/PGA lub ML/Md, jeśli tor skalibrowany i brak nasycenia (z zakresem zamiast jednej liczby).
• Lepiej ograniczyć się do opisu:
  • gdy istnieje ryzyko pomylenia zdarzenia z hałasem lokalnym,
  • gdy różne filtry prowadzą do sprzecznych wniosków co do faz,
  • gdy triangulacja opiera się na dwóch stacjach lub niespójnych zegarach.

Komunikując wstępne wyniki, trzy elementy porządkują przekaz: zakres niepewności, kontekst techniczny i język bez sugerowania pewności tam, gdzie jej nie ma. Zamiast „epicentrum było w X”, lepiej: „na stacji ABC zarejestrowano zdarzenie o cechach lokalnych; odległość oszacowana z S–P wynosi ok. 20–30 km od stacji (UTC 12:34:56), widoczne fale powierzchniowe”. Dodaj pasmo filtracji i jednostki, np. „prędkość gruntu, 0.5–10 Hz, zero‑phase”. Taka forma jasno rozdziela to, co stwierdzone, od tego, co dopiero wymaga sieciowego rozwiązania.

Sejsmogram czy hałas: decyzja rozpoznawcza w 60 sekund

Gdy masz tylko chwilę na rozstrzygnięcie, czy zapis to trzęsienie, czy lokalne zakłócenie, oprzyj się na kilku prostych wskaźnikach. Zestaw ich 2–3 i sprawdź spójność między kanałami.

Przykład z praktyki: krótki impuls o dużej amplitudzie, rozpoznany pierwotnie jako „P”, znika po lekkim zawężeniu pasma do 1–8 Hz i nie ma śladu na stacji 3 km dalej – decyzja: zakłócenie lokalne, nie publikować jako „lokalne trzęsienie”.

Wybuch czy trzęsienie: kiedy stawiać taką hipotezę

• Kiedy hipoteza „wybuch” ma sens:
  • dominacja wysokich częstotliwości i krótka coda w zapisie lokalnym,
  • duża amplituda P względem S (stosunek P/S nietypowo wysoki dla regionu),
  • silna fala Rg bardzo blisko (kilka–kilkanaście km) bez wyraźnej, późniejszej dominacji Rayleigh/Love.
• Kiedy raczej „tektonika”:
  • emergentny początek P i wyraźne, silniejsze S,
  • dłuższa coda z rozsądnie malejącą energią,
  • spójność z innymi stacjami i brak korelacji czasowej z harmonogramem strzelań w kopalniach.

W praktyce bywa różnie – pojedyncza stacja rzadko daje pewność. Jeśli po zmianie pasma P/S odwraca się lub zanika coda, odłóż kategoryzację do czasu weryfikacji sieciowej.

Back‑azymut z jednej stacji: kiedy wyznaczać kierunek źródła

Analiza polaryzacji i ruchu cząsteczkowego może z grubsza wskazać kierunek nadejścia fal, ale bywa wrażliwa na szum i orientację czujnika.

• Kiedy ma sens:
  • czytelne okno P o dobrym SNR na wszystkich kanałach,
  • znana orientacja N/E (błąd rzędu kilku stopni jest akceptowalny),
  • brak silnych rezonansów miejsca w wybranym paśmie (np. 1–8 Hz dla lokalnych).
• Jak przeprowadzić bezpiecznie (skrót procedury):
  1. wytnij krótkie okno P (np. 1–3 s dla lokalnych),
  2. zastosuj filtr pasmowy zero‑phase (np. 1–8 Hz) i usuń trend,
  3. policz macierz kowariancji kanałów N/E/Z i wektory własne,
  4. projekcja największej energii w płaszczyźnie poziomej daje szacowany back‑azymut,
  5. zweryfikuj na S (ruch eliptyczny w poziomie) i po rotacji N/E→R/T sprawdź, czy P jest głównie na R, a S na T.
• Kiedy odpuścić:
  • niski SNR lub nasycenie toru,
  • podejrzenie błędnej orientacji kanałów,
  • silne, wąskopasmowe wzbudzenia miejsca – wektory własne „gonią” rezonans, nie kierunek fali.

Krótki efekt decyzyjny: jeśli po rotacji na R/T energia P wyraźnie dominuje w R, a S w T w dwóch sąsiednich pasmach – kierunek jest użyteczny operacyjnie; w przeciwnym razie nie publikuj azymutu.

Widmo i spektrogram: kiedy wykres czasu to za mało

Widmo podpowiada, z jakich częstotliwości składa się sygnał; spektrogram pokazuje ich zmienność w czasie. To ułatwia odróżnienie fal objętościowych od powierzchniowych i hałasu.

• Wybieraj widmo/spektrogram, gdy:
  • chcesz ocenić zawartość wysokich częstotliwości (lokalne wstrząsy, strzały),
  • musisz rozseparować P/S od fal Rayleigh/Love,
  • różne filtry w dziedzinie czasu dają sprzeczne wnioski – widmo porządkuje obraz.
• Ustawienia, które zwykle działają:
  • okno 0.5–1 s dla lokalnych i 2–5 s dla telesejsm (kompromis czas/częstotliwość),
  • nakładanie 50–75% i skalowanie logarytmiczne amplitudy,
  • pre‑whitening lub łagodny pre‑emfaza, by nie „utopić” wyższych częstotliwości.
• Uważaj, gdy:
  • za krótkie okno „rozmywa” niskie częstotliwości – fale powierzchniowe mogą wyglądać na słabsze,
  • za długie okno spłaszcza nagły początek P,
  • normalizacja „na okno” w spektrogramie (każde okno skalowane oddzielnie) maskuje realne zmiany energii w czasie,
  • brak korekty odpowiedzi instrumentu zaniża wysokie częstotliwości i zniekształca porównania między stacjami,

  Skalowanie wielkości: kiedy liczyć magnitudę z jednej stacji

  Wielkość zdarzenia da się wstępnie oszacować z pojedynczego zapisu, ale tylko przy spełnieniu kilku warunków technicznych. Najpierw zdecyduj, którą definicję magnitudy w ogóle możesz zastosować.

  Kiedy „tak” – oblicz wstępnie	Kiedy „nie” – wstrzymaj się
  kanały niesaturujące, znana i usunięta odpowiedź instrumentu (np. do prędkości gruntu)	clipping choćby na kilku próbkach lub niepewna kalibracja — amplitudy są bezużyteczne
  odległość i pasmo właściwe dla danej skali (ML: ~10–600 km, Md: lokalnie z codą, mb: telesejsmy)	bardzo blisko źródła (pole bliskie, komponent statyczny) lub bardzo daleko bez korekcji tłumienia
  rozsądny SNR i stabilna korekcja tłumienia z odległością (Q, geometria)	silne efekty miejsca bez możliwości korekty (basen osadowy, rezonanse konstrukcji)

  • ML (magnituda lokalna): zwykle bazuje na skorygowanej amplitudzie (mm) i odległości. Działa najlepiej dla sieci krótkookresowych w skali regionalnej. Wymaga regionalnych krzywych tłumienia; bez nich podawaj zakres (np. ±0.3–0.5).
  • Md (od cody): sumaryczna energia cody w ustalonym oknie i paśmie. Bywa stabilniejsza dla małych lokalnych zdarzeń w hałasie. Wrażliwa na wybór okna i filtr.
  • mb/mbLg: dla telesejsm i odległości regionalnych; bez korekcji widmowej i okna fazy P łatwo o błąd systematyczny.
  • Mw: co do zasady nie do uzyskania wiarygodnie z jednej krótkookresowej stacji bez modelowania widmowego i momentu sejsmicznego; jeśli już, to jako bardzo wstępne oszacowanie na bazie widma źródła i częstotliwości narożnej.

  Procedura skrótowa (ML/Md): usuń trend i średnią, zdekonwoluj odpowiedź instrumentu do prędkości, zastosuj pasmo odpowiadające skali (np. 1–10 Hz lokalnie), wybierz okno amplitudy (ML) lub cody (Md), skoryguj o odległość/tłumienie według przyjętego modelu regionalnego, podaj wartość wraz z zakresem niepewności. Jeśli dwie metody rozjeżdżają się o >0.5, nie publikuj liczby bez dodatkowej weryfikacji.

  Krótki przykład: lokalny wstrząs z czytelną codą; ML z amplitudy daje 2.4, Md z cody 2.1. Różnica akceptowalna przy hałasie tła i słabym Rg — komunikat: „M~2.1–2.4 (lokalne), metoda ML/Md, 1–10 Hz, prędkość gruntu, brak saturacji”.

  PGA/PGV a odczuwalność i szkody: kiedy ekstrapolować wnioski

  Warto rozdzielić dwie decyzje: co wolno powiedzieć o zapisie w punkcie pomiaru i co można uogólnić na okolicę. Co do zasady uogólnienia są ryzykowne bez sieci.

  • Kiedy dopuszczalne są ostrożne wnioski „u mnie na stacji”:
    • PGA/PGV policzone po korekcie instrumentu i bez nasycenia,
    • pasmo zgodne z praktyką inżynierską (np. 0.1–20 Hz) lub jasno wskazane inne,
    • brak silnych rezonansów miejsca w analizowanym paśmie (albo wyraźna adnotacja, że one dominują).
  • Kiedy nie ekstrapolować na obszar:
    • znaczna różnorodność warunków gruntowych wokół (doliny vs wychodnie skalne),
    • bliskie źródło z efektem kierunkowości (directivity) — sąsiednie miejscowości mogą doświadczyć skrajnie różnych wstrząsów,
    • brak potwierdzenia z innych stacji lub raportów makrosejsmicznych.

  Orientacyjnie PGV rzędu pojedynczych mm/s zwykle oznacza słabe odczucia, a wartości rzędu kilku cm/s bywają już wyraźnie odczuwalne. W praktyce bywa różnie — efekt miejsca potrafi zmienić obraz o rząd wielkości.

  Pole bliskie i efekty źródła: kiedy wykres wprowadza w błąd

  Bardzo blisko ogniska mechanika źródła i geometria fal mogą zdominować zapis i złamać typowe „reguły” P/S.

  • Scenariusze, w których klasyczne wskaźniki zawodzą:
    • skok statyczny (fling step) na accelerometrze lub długookresowy dryf na velocymetrze — amplitudy i PGV/PGA stają się mylące bez specjalnej obróbki,
    • węzły promieniowania — jedna z faz bywa nienaturalnie słaba, co podbija lub zaniża stosunek P/S,
    • silna bezpośredniość (rupture directivity) — coda krótka i „wysokoczęstotliwościowa” mimo tektonicznego charakteru.
  • Jak reagować decyzji:
    • porównaj sąsiednie stacje wzdłuż i w poprzek przypuszczalnego kierunku pękania,
    • sprawdź spójność w dwóch pasmach (np. 0.5–2 Hz i 2–8 Hz); jeśli wnioski zmieniają się o 180°, wstrzymaj kategoryzację,
    • oddziel komponent statyczny/ bardzo niskoczęstotliwościowy przed liczeniem PGV/PGA i magnitudy z cody.

  Lista kontrolna decyzji z pojedynczej stacji

  • czy instrument nie nasycił się i czy odpowiedź została poprawnie usunięta,
  • czy P i S są spójne na Z/N/E i w co najmniej dwóch pasmach,
  • czy odległość z S–P mieści się w zakresie wybranej skali magnitudy,
  • czy fale powierzchniowe nie dominują wskutek efektu miejsca,
  • czy sąsiednia stacja pokazuje choć zarys tego samego zdarzenia,
  • czy wynik jest stabilny po niewielkiej zmianie okien/filtrowania,
  • czy podajesz zakres niepewności zamiast jednej liczby,
  • czy nie próbujesz wnioskować ponad to, co wspiera sygnał (np. mechanizm źródła bez sieci).

  Automaty detekcji i progi: kiedy podkręcić, a kiedy przykręcić czułość

  Systemy STA/LTA i wielopasmowe wyzwalanie przyspieszają reakcję, ale konfiguracja to również decyzje o kompromisach.

  • Podnieś czułość (niższy próg), gdy:
    • sieć ma gęste pokrycie, a fałszywe alarmy są szybko weryfikowane,
    • celem jest rejestracja mikrosejsmiczności w spokojnym środowisku szumowym.
  • Obniż czułość (wyższy próg) lub dodaj warunki, gdy:
    • pojedyncza stacja w hałaśliwym otoczeniu generuje wiele fałszywych wzbudzeń,
    • brak synchronizacji czasu grozi błędną korelacją między stacjami,
    • pracują maszyny o znanych, wąskopasmowych częstotliwościach — włącz bramkowanie pasmowe lub test koherencji między kanałami.
  • włącz test stałej polarności między kanałami (np. czy Z jest zgodny fazowo z R dla P po rotacji); nagłe odwrócenia fazy między dniami to zwykle problem kablowy, nie „dziwne” zdarzenie.

Po każdej zmianie progów sprawdź skutki na historii danych: liczba wzbudzeń na dobę, odsetek fałszywych alarmów, czas do pierwszej próbki po progu i czytelność P. Jeśli znikają realne zdarzenia o znanej charakterystyce (np. kontrolne strzały w kopalni), czułość jest ustawiona zbyt agresywnie.

Synchronizacja czasu i odpowiedź instrumentu: punkt zerowy wiarygodności

Bez prawidłowego czasu i kalibracji amplitud każdy wykres jest tylko „ładnym obrazkiem”. W praktyce błąd o kilka sekund potrafi odwrócić wnioski o S–P, a brak odpowiedzi instrumentu przekłamuje amplitudy nawet o rząd wielkości.

• Czas:
  • zwykle GNSS/NTP w trybie ciągłym; dryf lokalnego zegara porównuj dobę do doby,
  • porównuj czasy przybycia z sąsiednimi stacjami (korelacja cody) — przesunięcia stałe są sygnałem problemu,
  • po utracie synchronizacji oznacz segment jako niepewny — automatyczne asocjacje w sieci bywają wówczas bezwartościowe.
• Odpowiedź instrumentu:
  • trzymaj komplet plików kalibracyjnych i wersjonuj je; po serwisie przelicz widma szumu własnego,
  • dekonwolucja do prędkości lub przyspieszenia gruntu powinna być jawnie wskazana przy każdym wykresie i liczbie,
  • sprawdź liniowość: jeśli amplitudy bliskie zakresu pełnego powtarzają formę „spłaszczonych” grzbietów, zapis jest nasycony nawet bez oczywistego wykoszenia.

Co do zasady, każda decyzja ilościowa (magnituda, PGV/PGA, porównania między stacjami) bez uprzedniej korekty odpowiedzi instrumentu i kontroli czasu powinna być traktowana jako hipoteza robocza, nie wynik.

Jak czytać helikorder i typowy wykres sejsmogramu

Helikorder to dzienny zapis ułożony w wierszach; każda linia odpowiada określonemu przedziałowi czasu (np. 30 lub 60 minut). Zmiana koloru często sygnalizuje kolejną linię lub zakresu czasu, nie „siłę” zdarzenia. Skompresowane amplitude łatwo mylą oczy — duże zdarzenie może wyglądać jak „czarna taśma”, a lokalne, krótkie wstrząsy jak pojedyncze igły.

• Na co patrzeć najpierw:
  • czytelny początek P i późniejszy wzrost energii (S, fale powierzchniowe),
  • odstęp S–P i jego stabilność między komponentami — to najszybsza wskazówka odległości,
  • ślad nasycenia (płaskie „ścięte” grzbiety, nagłe przejścia do zakresu skrajnego),
  • wzory powtarzalne o stałej częstotliwości (maszyny, wiatr na maszcie) kontra pojedynczy impuls i coda,
  • czy zdarzenie nie „przechodzi” przez koniec wiersza — kontynuacja bywa na początku kolejnej linii,
  • różnice między Z a N/E: P zwykle dominuje na Z, S i fale powierzchniowe silniejsze na poziomych.

Jeśli helikorder jest mocno skompresowany, sięgnij po podgląd pojedynczego kanału z większą rozdzielczością. Ocena na podstawie „czarnej wstążki” prowadzi do nadinterpretacji — zamiast mocy sygnału często oglądamy kompresję graficzną.

Najczęstsze pomyłki na helikorderze

• zliczanie amplitud w milimetrach z ekranu bez znajomości skali i odpowiedzi instrumentu,
• mylenie lokalnego hałasu antropogenicznego (wąskie pasma, brak wyraźnej cody) z małym trzęsieniem,
• uznanie równoległych, rytmicznych impulsów za serię wstrząsów — często to „ripple-fire” zrobót strzałowych lub praca młota,
• przykładanie globalnych reguł (np. kolejności faz telesejsmicznych) do bliskich zdarzeń, gdzie dominuje pole bliskie i fale powierzchniowe.

Spektrogram jako szybka lupa

Sejsmogram mówi „kiedy” i „jak duże”, spektrogram podpowiada „w jakim paśmie”. W praktyce tandem tych dwóch widoków znacząco przyspiesza decyzje.

• Impulsy o szerokim paśmie, z szybkim spadkiem energii ku niskim częstotliwościom — typowe dla małych lokalnych wstrząsów.

Najczęściej zadawane pytania (FAQ)

Skąd wiem, czy to na sejsmogramie to trzęsienie ziemi, a nie hałas?

Co do zasady prawdziwe zdarzenie sejsmiczne ma kolejność faz: najpierw drobna, wyraźniejsza na kanale Z energia fal P, po chwili silniejsze fale S na kanałach poziomych, a następnie wolniejsze fale powierzchniowe. Hałas (ruch uliczny, wiatr, prace budowlane) bywa nieregularny, bez czytelnej sekwencji P–S i zwykle nie pojawia się jednocześnie na kilku stacjach.

Praktyczne kryteria rozróżnienia:

• Sprawdź kilka kanałów: P powinna wejść najpierw na Z, S dominować na N/E.
• Porównaj ze stacją w okolicy: jeśli obie widzą podobny przebieg w zbliżonym czasie – to nie jest lokalny hałas.
• Oceń tło: w dzień w mieście szum wysokoczęstotliwościowy bywa podwyższony; nocą sygnały są czystsze.
• Szukaj „clippingu”: płaskie, ścięte wierzchołki to nasycenie instrumentu – nie dowód „wielkiej siły”.

Przykład z praktyki: krótki, impulsowy pik tylko na jednym kanale poziomym w południe obok drogi to zwykle ciężarówka, nie trzęsienie.

Jak na sejsmogramie rozpoznać fale P i S oraz fale powierzchniowe?

Zacznij od kanałów: Z (pion) najlepiej pokazuje pierwsze nadejście P – sygnał bywa subtelny, wyższej częstotliwości i o mniejszej amplitudzie. Po czasie Δt pojawiają się fale S, zwykle wyraźniejsze na kanałach N i E, z większą amplitudą i nieco niższą częstotliwością.

Jeszcze później pojawiają się fale powierzchniowe (Rayleigha i Love’a). Mają największe amplitudy i najniższe częstotliwości; dominują przy silniejszych i dalszych zdarzeniach. Jeśli masz wątpliwość, nałóż łagodny filtr: 1–10 Hz ułatwia wyłapanie lokalnych P/S, a 0,02–0,2 Hz uwydatnia fale powierzchniowe.

Jak obliczyć odległość do epicentrum z różnicy czasów S–P?

Zmierz na tej samej stacji czas nadejścia P i S, a następnie policz Δt = t(S) − t(P). Przy typowych prędkościach w skorupie (P ~ 6 km/s, S ~ 3,5 km/s) odległość D szacuje się orientacyjnie jako 8–9 razy Δt w sekundach. Przykład: Δt = 10 s → D ≈ 80–90 km.