Kalkulačka triangulácie epicentra zemetrasenia

Fyzikálne princípy seizmických vĺn a triangulácie

Pri vzniku zemetrasenia sa zo zdroja šíria rôzne druhy seizmických vĺn. Pre lokalizáciu epicentra sú kľúčové primárne kompresné vlny (P-vlny) a sekundárne strihové vlny (S-vlny). P-vlny sa pohybujú rýchlejšie ako S-vlny, čo znamená, že s rastúcou vzdialenosťou od zdroja sa časový interval medzi ich príchodmi na registračnú stanicu zväčšuje. Tento rozdiel v čase príchodu sa označuje ako interval S−P.

Jedna seizmická stanica dokáže z nameraného intervalu S−P určiť iba polomer kružnice, na ktorej sa zemetrasenie nachádza, nie však konkrétny smer. Dve stanice vymedzujú pretínajúce sa kružnice, ktoré sa zvyčajne pretínajú v dvoch bodoch. Na jednoznačné určenie polohy epicentra sú preto potrebné minimálne tri stanice.

Kalkulačka triangulácie epicentra zemetrasenia funguje ako bezplatný online výučbový nástroj, ktorý umožňuje odhadnúť epicentrum na sférickom modeli Zeme s konštantným polomerom R = 6 371 km. Výpočet využíva pozorovania z troch až ôsmich staníc, pričom zohľadňuje používateľom definované rýchlosti vĺn a časovú neistotu.

Vstupné parametre a obmedzenia modelu

Pre úspešný výpočet vyžaduje kalkulačka zadanie špecifických hodnôt pre každú stanicu a parametrov šírenia vĺn. Používateľ môže kliknúť na tlačidlo Načítať príklad Tokio pre automatické vyplnenie vzorových dát, alebo na Vymazať pre resetovanie všetkých polí.

Medzi povinné vstupy patria:

• Pozorovania staníc (minimálne 3 a maximálne 8 staníc):
  • Názov stanice
  • Zemepisná šírka: Musí byť v rozsahu od −90° do 90°.
  • Zemepisná dĺžka: Musí byť v rozsahu od −180° do 180°.
  • Interval S−P: Hodnota v sekundách, ktorá musí byť väčšia ako 0 a najviac 2 000 sekúnd.
• Model šírenia:
  • Rýchlosť P-vlny: Hodnota musí byť väčšia ako 0, maximálne 20 km/s a zároveň musí byť väčšia ako rýchlosť S-vlny.
  • Rýchlosť S-vlny: Hodnota musí byť väčšia ako 0, maximálne 20 km/s a zároveň musí byť menšia ako rýchlosť P-vlny.
  • Časová neistota (1σ): Hodnota musí byť väčšia ako 0 a najviac 60 sekúnd.

Spracovanie všetkých zadaných údajov prebieha lokálne v prehliadači používateľa. Súradnice staníc, časy S−P a výsledky zostávajú v tomto zariadení a neodosielajú sa na žiadne externé servery.

Matematické vyjadrenie a výpočet vzdialenosti

Nástroj aproximuje vzdialenosť k epicentru na základe predpokladu konštantnej rýchlosti vĺn v homogénnej zemskej kôre. Vzťah pre výpočet vzdialenosti vychádza z rovnice:

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

Po úprave dostávame vzorec pre odhad vzdialenosti (D) z časového rozdielu (Δt):

Aproximácia vzdialenosti: D = Δt ÷ (1/Vs − 1/Vp)

Pre prevod uhlových vzdialeností na sférickom povrchu Zeme sa používa haversine vzorec:

Povrchová vzdialenost: d = R × stredový uhol haversine, pričom R = 6371 km

Keďže reálne merania obsahujú nepresnosti, kružnice jednotlivých staníc sa málokedy pretnú v jednom presnom bode. Kalkulačka preto hľadá optimálne riešenie minimalizáciou štvorcov odchýlok medzi nameranými a modelovanými vzdialenosťami:

Prispôsobenie epicentra: minimalizovať Σ(dᵢ − Dᵢ)² cez všetky stanice

Interpretácia výsledkov a kvalita prispôsobenia

Po prebehnutí výpočtu nástroj zobrazí výsledné riešenie v sekcii Riešenie epicentra. Výstup zahŕňa tieto položky:

• Odhadované epicentrum: Výsledná zemepisná šírka a zemepisná dĺžka.
• Kvalita prispôsobenia vzdialenosti: Označená ako tesné prispôsobenie vzdialenosti, zmiešané prispôsobenie vzdialenosti alebo slabé prispôsobenie vzdialenosti.
• Kvalita geometrie staníc: Označená ako silná geometria staníc, stredná geometria staníc alebo slabá geometria staníc.
• RMS rezíduum vzdialenosti a Najväčšie rezíduum vzdialenosti.
• Približná neistota polohy 1σ a Neistota vzdialenosti 1σ na stanicu.

Súčasťou rozhrania je aj Schéma geometrie staníc, ktorá vizualizuje vzájomné priestorové usporiadanie staníc a odhadovaného epicentra. Táto schéma obsahuje upozornenie: „Iba relatívna geometria – bez mapy, hraníc alebo živých dát o zemetraseniach.“.

V sekcii Vzorec a dosadenie sa zobrazujú konkrétne kroky výpočtu vrátane odvodenia faktora vzdialenosti:

• Faktor vzdialenosti = vp × vs ÷ (vp − vs) = factor km na sekundu S−P
• Jednotlivé odvodené vzdialenosti pre každú stanicu v tvare: station: interval s × factor km/s = distance km
• Výsledné súradnice a RMS rezíduum: Najlepšie sférické prispôsobenie = latitude°, longitude°; RMS rezíduum = rms km
• Propagácia neistoty: Časová neistota 1σ s hodnotou timing s dáva neistotu vzdialenosti 1σ s hodnotou distance km; geometria staníc ju na hlavnej horizontálnej osi propaguje na približne position km.

Chybové stavy a validácia vstupov

Nástroj dôsledne kontroluje integritu zadávaných dát a geometrickú stabilitu siete staníc. Pri zistení nezrovnalostí zobrazuje nasledujúce chybové hlásenia:

• Nedostatočný alebo nadmerný počet staníc: Ak sa pokúsite spustiť výpočet s menej ako 3 alebo viac ako 8 stanicami, zobrazí sa: Použite tri až osem staníc.
• Chyba formátu: Ak zadaná hodnota nie je číselná: where: „token“ nie je číslo.
• Prekročenie rozsahov súradníc:
  • Stanica station: zemepisná šírka musí byť od −90° do 90°.
  • Stanica station: zemepisná dĺžka musí byť od −180° do 180°.
• Chybný interval: Stanica station: zadajte interval S−P väčší ako 0 a najviac 2000 sekúnd.
• Chybné rýchlosti vĺn: field musí byť väčšia ako 0 a nie viac ako 20 km/s. prípadne Rýchlosť P-vlny musí byť väčšia ako rýchlosť S-vlny.
• Chybná neistota: Časová neistota musí byť väčšia ako 0 a najviac 60 sekúnd.
• Duplicitné polohy: Stanice station a other majú rovnaké súradnice. Použite odlišné polohy staníc.
• Lineárne alebo zhlukové usporiadanie: Ak sú stanice usporiadané v priamke alebo tesnej skupine, výpočet nie je stabilný: Tieto stanice nevymedzujú stabilnú dvojrozmernú polohu. Pridajte stanicu mimo aktuálnej línie alebo zhluku.
• Pretečenie rozsahu: Hodnota alebo medzivýsledok presahuje podporovaný rozsah čísel.

Často kladené otázky (FAQ)

Prečo sú potrebné tri stanice?

Jedno meranie S−P poskytuje kružnicu vzdialenosti, nie smer. Dve kružnice sa môžu pretínať na dvoch miestach. Tretia stanica zvyčajne túto voľbu vyrieši, zatiaľ čo štvrtá alebo ďalšia stanica zviditeľní nesúlad prostredníctvom rezíduí. Stanice rozmiestnené okolo zdroja vymedzujú zemepisnú šírku aj dĺžku lepšie ako stanice pozdĺž jednej línie.

Je vzdialenosť S−P presnou povrchovou vzdialenosťou?

Nie. Rovnica s konštantnou rýchlosťou odhaduje spoločnú dráhu šírenia z rozdielu časov príchodu. Považovanie tohto polomeru za sférickú povrchovú vzdialenosť predpokladá plytké ohnisko a jednoduchú rýchlostnú štruktúru. Skutočné krivky hodochron závisia od hĺbky, vzdialenosti a materiálu, ktorým každý lúč prechádza.

Čo znamená kladné alebo záporné rezíduum?

Rezíduum je prispôsobená sférická vzdialenosť mínus vzdialenosť odvodená zo S−P. Kladná hodnota znamená, že prispôsobený bod je od danej stanice ďalej než jej polomer S−P; záporná hodnota znamená, že je bližšie. Veľké zmiešané rezíduá ukazujú, že kružnice pri zvolených rýchlostiach nemajú jeden čistý spoločný priesečník.

Čo je zahrnuté v neistote polohy?

Propaguje jedinú hodnotu časovej neistoty 1σ cez faktor vzdialenosti S−P a geometriu staníc v blízkosti prispôsobeného bodu. Nezahŕňa nesprávny rýchlostný model, hĺbku ohniska, chybne označené fázy, odchýlku hodín ani trojrozmernú štruktúru Zeme, takže nejde o úplné vyhlásenie o presnosti.