Kalkulator triangulacji epicentrum trzęsienia ziemi

Fizyka fal sejsmicznych i interwał S−P

Lokalizowanie epicentrum trzęsienia ziemi opiera się na różnicy prędkości rozchodzenia się fal sejsmicznych w skorupie ziemskiej. Podczas nagłego przesunięcia mas skalnych w ognisku (hipocentrum) generowane są różne rodzaje fal, z których najważniejsze dla celów lokalizacji są fale podłużne P (ang. primary lub compressional) oraz fale poprzeczne S (ang. secondary lub shear).

Fale P poruszają się szybciej niż fale S. W rezultacie, w miarę oddalania się od źródła wstrząsu, odstęp czasowy między przyjściem fali P a fali S stale rośnie. Mierząc ten odstęp – określany jako interwał S−P (Δt) – na stacji sejsmologicznej, można wyznaczyć odległość stacji od epicentrum. Pojedyncza stacja rejestrująca określa jedynie promień okręgu, na którym leży źródło. Dopiero zestawienie danych z wielu stacji pozwala na jednoznaczne wskazanie położenia geograficznego wstrząsu.

Podstawowe równanie wiążące różnicę czasu przyjścia fal z odległością (D) przy założeniu stałych prędkości fal P (Vₚ) i S (Vₛ) ma postać:

Δt = D × (1/Vs − 1/Vp)

Po przekształceniu pozwala ono bezpośrednio obliczyć odległość stacji od źródła:

D = Δt ÷ (1/Vs − 1/Vp)

Geometria triangulacji i dopasowanie sferyczne

Tradycyjna metoda triangulacji na płaszczyźnie zakłada, że okręgi narysowane wokół trzech stacji o promieniach odpowiadających wyznaczonym odległościom przetną się dokładnie w jednym punkcie. W rzeczywistości, z powodu błędów pomiarowych, uproszczeń modeli prędkości oraz sferycznego kształtu Ziemi, okręgi te rzadko przecinają się idealnie w jednym miejscu.

Kalkulator przyjmuje model Ziemi jako kulę o stałym promieniu R = 6371 km. Odległość powierzchniowa (d) między współrzędnymi stacji a potencjalnym epicentrum obliczana jest przy użyciu wzoru haversine dla kąta środkowego:

d = R × kąt środkowy haversine, gdzie R = 6371 km

Aby znaleźć najbardziej prawdopodobne położenie epicentrum przy użyciu danych z od 3 do 8 stacji, narzędzie stosuje metodę najmniejszych kwadratów. Algorytm minimalizuje sumę kwadratów różnic pomiędzy sferyczną odległością geometryczną (d) a odległością obliczoną z interwału S−P (D) dla wszystkich stacji:

Dopasowanie epicentrum: minimalizacja Σ(dᵢ − Dᵢ)² dla wszystkich stacji

Dzięki temu podejściu, zamiast polegać na idealnym przecięciu trzech okręgów, kalkulator znajduje optymalne matematycznie rozwiązanie dla całej sieci stacji.

Parametry wejściowe i ograniczenia modelu

Do poprawnego obliczenia położenia epicentrum wymagane jest wprowadzenie danych z minimum 3 i maksimum 8 stacji. Wprowadzenie mniejszej lub większej liczby stacji spowoduje wyświetlenie komunikatu o błędzie: Użyj od trzech do ośmiu stacji.

Dla każdej stacji użytkownik określa następujące parametry:

• Szerokość geograficzna: Wartość liczbowa w przedziale od −90° do 90°. Przekroczenie tego zakresu wywołuje błąd: Stacja 1: szerokość geograficzna musi mieścić się w przedziale od −90° do 90°.
• Długość geograficzna: Wartość liczbowa w przedziale od −180° do 180°. Przekroczenie tego zakresu wywołuje błąd: Stacja 1: długość geograficzna musi mieścić się w przedziale od −180° do 180°.
• Interwał S−P: Wartość większa niż 0 i nie większa niż 2000 sekund. Nieprawidłowa wartość generuje błąd: Stacja 1: wpisz interwał S−P większy niż 0 i nie większy niż 2000 sekund.

W sekcji Model propagacji definiuje się parametry prędkości fal oraz niepewności pomiarowych:

• Prędkość fali P: Musi być większa niż 0 i nie większa niż 20 km/s.
• Prędkość fali S: Musi być większa niż 0 i nie większa niż 20 km/s.
• Prędkość fali P musi być większa od prędkości fali S.
• Niepewność czasu (1σ): Musi być większa niż 0 i nie większa niż 60 sekund.

Dla niepewności czasu limit określony jest komunikatem: Niepewność czasu musi być większa niż 0 i nie większa niż 60 sekund. W przypadku wprowadzenia wartości nieliczbowej w dowolnym polu, kalkulator wyświetli błąd: input: „abc” nie jest liczbą.

Użytkownik może skorzystać z przycisku Załaduj przykład z Tokio, aby automatycznie wypełnić formularz przykładowymi danymi, lub użyć przycisku Wyczyść, aby zresetować wszystkie pola.

Analiza wyników i interpretacja odchyłek

Po wykonaniu obliczeń w sekcji Rozwiązanie epicentrum prezentowane są następujące wyniki:

• Szacowane epicentrum: Wyznaczone współrzędne szerokości i długości geograficznej.
• Dopasowanie odległości: Jakość dopasowania opisywana jest jako bliskie dopasowanie odległości, mieszane dopasowanie odległości lub słabe dopasowanie odległości.
• Geometria stacji: Ocena układu przestrzennego stacji, oznaczana jako silna geometria stacji, umiarkowana geometria stacji lub słaba geometria stacji.
• Odchyłka RMS odległości (Root Mean Square) oraz Największa odchyłka odległości.
• Przybliżona niepewność pozycji 1σ oraz Niepewność odległości 1σ na stację.

W sekcji Dopasowanie stacji znajduje się tabela z kolumnami: Stacja, Odległość S−P, Odległość sferyczna oraz Odchyłka. Odchyłka (residual) to różnica między dopasowaną odległością sferyczną a odległością wynikającą z interwału S−P. Dodatnia odchyłka oznacza, że obliczone epicentrum leży dalej od stacji niż promień S−P, natomiast ujemna – że leży bliżej.

Wizualną reprezentację układu stacji przedstawia Schemat geometrii stacji opatrzony podpisem: „Tylko geometria względna — brak mapy, granic czy rzeczywistych danych o trzęsieniach ziemi.”

W sekcji Metoda, data i zakres wyświetlane są komunikaty o statusie obliczeń:

• Edytuj obserwacje stacji, aby ponownie dopasować epicentrum.

Metodologia obliczeń opiera się na standardach USGS. Informacja o weryfikacji metody wskazuje datę 6 sierpnia 2026.

Często zadawane pytania (FAQ)

Czy odległość S−P jest dokładną odległością powierzchniową? Nie. Równanie ze stałą prędkością szacuje wspólną drogę przebiegu fali na podstawie różnicy czasu przyjścia. Traktowanie tego promienia jako sferycznej odległości powierzchniowej zakłada płytkie źródło i prostą strukturę prędkości. Rzeczywiste krzywe czasu przebiegu zależą od głębokości, odległości i materiału, przez który przechodzi każdy promień.

Dlaczego potrzebne są trzy stacje? Jeden pomiar S−P daje okrąg odległości, a nie kierunek. Dwa okręgi mogą przecinać się w dwóch miejscach. Trzecia stacja zazwyczaj pozwala rozstrzygnąć ten wybór, natomiast czwarta i kolejne stacje uwidaczniają niedopasowania poprzez odchyłki. Stacje rozmieszczone wokół źródła lepiej określają zarówno szerokość, jak i długość geograficzną niż stacje ułożone w jednej linii.

Co oznacza dodatnia lub ujemna odchyłka? Odchyłka to dopasowana odległość sferyczna minus odległość wyznaczona z S−P. Wartość dodatnia oznacza, że dopasowany punkt leży dalej od danej stacji niż jej promień S−P; wartość ujemna oznacza, że leży bliżej. Duże, zróżnicowane odchyłki pokazują, że okręgi nie przecinają się czysto w jednym punkcie przy wybranych prędkościach.

Co jest uwzględnione w niepewności pozycji? Przenosi ona pojedynczą wartość niepewności czasu 1σ przez współczynnik odległości S−P oraz geometrię stacji w pobliżu dopasowanego punktu. Nie uwzględnia ona błędnego modelu prędkości, głębokości ogniska, błędnie zidentyfikowanych faz, błędu zegara stacji ani trójwymiarowej struktury Ziemi, więc nie stanowi pełnego określenia dokładności.