Fyzikální principy seismických vln a triangulace
Při vzniku zemětřesení se z ohniska šíří seismické vlny různých typů. Nejdůležitější pro určení polohy epicentra jsou vlny P (primární, podélné) a vlny S (sekundární, příčné). Vzhledem k tomu, že P-vlny se pohybují horninami zemské kůry rychleji než S-vlny, časový odstup mezi jejich příchody na registrační stanici roste s rostoucí vzdáleností od zdroje. Tento časový rozdíl se označuje jako interval S−P.
Jedna seismická stanice dokáže z naměřeného intervalu S−P určit pouze radiální vzdálenost, ve které se zemětřesení odehrálo. Tato vzdálenost definuje kružnici možných poloh epicentra se středem v dané stanici. Protnutím kružnic ze dvou stanic získáme zpravidla dva možné body průsečíku. K jednoznačnému určení polohy epicentra jsou proto zapotřebí minimálně tři stanice. Pokud je k dispozici více stanic, kružnice se v důsledku reálných nepřesností v měření neprotnou v jediném bodě. V takovém případě se k nalezení nejvhodnějšího epicentra používá matematická optimalizace metodou nejmenších čtverců.
Vstupy a parametry modelu
Kalkulačka umožňuje zadat pozorování pro 3 až 8 seismických stanic. Pro každou stanici je nutné vyplnit následující údaje:
- Zeměpisná šířka: Musí být v rozsahu od −90° do 90°.
- Zeměpisná délka: Musí být v rozsahu od −180° do 180°.
- interval S−P: Časový rozdíl mezi příchodem P a S vln. Hodnota musí být větší než 0 a nejvýše 2000 sekund.
Dále se definují parametry modelu šíření vln v jednotné zemské kůře:
- Rychlost P-vlny: Musí být větší než 0 a nejvýše 20 km/s. Zároveň musí být tato rychlost vyšší než rychlost S-vlny.
- Rychlost S-vlny: Musí být větší než 0 a nejvýše 20 km/s. Zároveň musí být tato rychlost nižší než rychlost P-vlny.
- Nejistota časování (1σ): Směrodatná odchylka určení času příchodu vln. Musí být větší než 0 a nejvýše 60 sekund.
Pro rychlé otestování funkčnosti lze kliknout na tlačítko Načíst příklad Tokio, které vyplní kalkulačku vzorovými daty, nebo na tlačítko Vymazat pro resetování všech polí.
Matematický aparát a výpočty
Výpočet odhadu polohy epicentra probíhá na zjednodušeném sférickém modelu Země s konstantním poloměrem R = 6371 km. Výpočetní postup se skládá z několika kroků:
D = Δ t ÷ (1/Vₛ - 1/Vₚ)
Tento vztah předpokládá konstantní rychlosti vln a mělký zdroj zemětřesení.
d = R × haversine středového úhlu
minimalizovat Σ(dᵢ - Dᵢ)²
Reziduum pro každou stanici se definuje jako sférická vzdálenost minus vzdálenost S−P. Kladné reziduum znamená, že výsledné proložené epicentrum leží dále od stanice, než odpovídá naměřenému intervalu S−P. Záporné reziduum značí, že epicentrum je stanici blíže.
Interpretace výsledků a kvality řešení
Nástroj po dokončení optimalizace zobrazí výsledné souřadnice v sekci Odhadované epicentrum a vyhodnotí kvalitu řešení:
- Kvalita proložení vzdáleností: Výsledek je označen jako těsné přizpůsobení vzdálenosti, smíšené přizpůsobení vzdálenosti nebo špatné přizpůsobení vzdálenosti. Hodnotí se pomocí parametrů RMS reziduum vzdálenosti (střední kvadratická chyba) a Největší reziduum vzdálenosti.
- Kvalita geometrie stanic: Označuje se jako silná geometrie stanic, střední geometrie stanic nebo slabá geometrie stanic. Pokud jsou stanice rozmístěny rovnoměrně kolem epicentra, je geometrie silná. Pokud leží v úzkém shluku nebo v přímce, je geometrie slabá a poloha epicentra je v některém směru špatně definovaná.
- Schéma geometrie stanic: Zobrazuje relativní polohu stanic a epicentra. Schéma doprovází popisek: „Pouze relativní geometrie – bez mapy, hranic nebo živých dat o zemětřeseních.“
- Tabulka proložení stanic: Obsahuje sloupce Stanice, Vzdálenost S−P, Sférická vzdálenost a Reziduum.
V sekci Metoda, datum a rozsah se zobrazují stavové zprávy o úspěšnosti výpočtu:
- „Upravte pozorování stanic pro nové proložení epicentra.“
Pravidla validace a chybová hlášení
Kalkulačka provádí okamžitou kontrolu zadávaných dat. Při porušení pravidel se zobrazí příslušná chybová hlášení:
- Nedodržení limitu počtu stanic (3 až 8):
Použijte tři až osm stanic. - Chybný poměr rychlostí:
Rychlost P-vlny musí být větší než rychlost S-vlny. - Nejistota časování mimo rozsah:
Nejistota časování musí být větší než 0 a nejvýše 60 sekund. - Lineární nebo shlukové uspořádání stanic:
Tyto stanice nevymezují stabilní dvourozměrnou polohu. Přidejte stanici mimo stávající linii nebo shluk. - Přetečení rozsahu čísel:
Hodnota nebo mezivýsledek přesahuje podporovaný rozsah čísel.
Omezení modelu a ochrana soukromí
Tento nástroj slouží primárně pro výukové účely a demonstraci matematických principů lokalizace zemětřesení. Model předpokládá dokonale sférickou Zemi a konstantní rychlost šíření vln v homogenním prostředí. Zanedbává hloubku ohniska (hypocentra), zakřivení drah seismických paprsků v reálném zemském tělese a existenci rychlostně odlišných vrstev zemské kůry a pláště. Vypočtená nejistota polohy zohledňuje pouze uživatelem zadanou chybu měření času (1σ) a geometrické uspořádání stanic, nikoli systematické chyby rychlostního modelu či přístrojové odchylky.
Zpracování všech zadaných dat probíhá lokálně přímo ve vašem webovém prohlížeči. Souřadnice stanic, naměřené časy ani výsledky výpočtů se neodesílají na žádné externí servery.
Často kladené otázky (FAQ)
Je vzdálenost S−P přesnou povrchovou vzdáleností? Ne. Rovnice s konstantní rychlostí odhaduje společnou dráhu šíření z rozdílu časů příchodu. Považovat tento poloměr za sférickou povrchovou vzdálenost předpokládá mělký zdroj a jednoduchou strukturu rychlostí. Skutečné křivky doby šíření závisí na hloubce, vzdálenosti a materiálu, kterým každý paprsek prochází.
Proč jsou potřeba tři stanice? Jedno měření S−P dává kružnici vzdálenosti, nikoli směr. Dvě kružnice se mohou protnout na dvou místech. Třetí stanice obvykle tuto volbu vyřeší, zatímco čtvrtá nebo další stanice zviditelní nesoulad prostřednictvím reziduí. Stanice rozmístěné kolem zdroje vymezují zeměpisnou šířku i délku lépe než stanice v jedné linii.
Co znamená kladné nebo záporné reziduum? Reziduum je proložená sférická vzdálenost minus vzdálenost odvozená ze S−P. Kladná hodnota znamená, že proložený bod je od dané stanice dále než její poloměr S−P; záporná hodnota znamená, že je blíže. Velká smíšená rezidua ukazují, že kružnice při zvolených rychlostech nesdílejí jeden čistý průsečík.
Co je zahrnuto v nejistotě polohy? Propaguje jedinou hodnotu časování 1σ přes faktor vzdálenosti S−P a geometrii stanic v blízkosti proloženého bodu. Nezahrnuje nesprávný rychlostní model, hloubku ohniska, chybně určené fáze, odchylku hodin ani trojrozměrnou strukturu Země, takže se nejedná o kompletní vyjádření přesnosti.