Beregner til triangulering af jordskælvsepicenter
Beregneren til triangulering af jordskælvsepicenter er et pædagogisk værktøj, der gør det muligt at estimere epicentret for et jordskælv på en sfærisk jordmodel ved hjælp af observationer fra tre til otte målestationer. Ved at indtaste koordinater og S−P-bølgeankomstintervaller for hver station, sammen med brugerdefinerede bølgehastigheder og tidsmæssig usikkerhed, beregner værktøjet den bedst tilpassede placering af epicentret. Værktøjet anvender en mindste kvadraters metode til at minimere afvigelserne i afstand på tværs af alle stationer. Det viser det estimerede epicenter, styrken af stationsgeometrien samt de enkelte stationers residualværdier for at illustrere, hvor godt observationerne stemmer overens.
Alle indtastede stationskoordinater, S−P-tider og beregningsresultater bliver på denne enhed og uploades ikke til eksterne servere.
Fysikken bag seismiske bølger og S−P-intervallet
Når et jordskælv finder sted, frigives der energi i form af seismiske bølger, som udbreder sig gennem Jorden. De to primære kropsbølger, der registreres på seismiske målestationer, er P-waves (primære eller kompressionsbølger) og S-waves (sekundære eller forskydningsbølger).
P-waves bevæger sig hurtigere end S-waves. Derfor vil tidsintervallet mellem deres ankomster til en målestation vokse i takt med, at afstanden fra jordskælvets kilde øges. Under antagelsen om konstante bølgehastigheder i en ensartet jordskorpe opfylder forskellen i ankomsttid ligningen:
Δt = D × (1/Vs − 1/Vp)
Hvor:
- Δt er S−P-intervallet i sekunder.
- D er distancen fra stationen til epicentret.
- Vp er P-wave-hastighed.
- Vs er S-wave-hastighed.
Ved at omforme denne ligning kan distancen fra en given målestation til epicentret beregnes direkte som:
D = Δt ÷ (1/Vs − 1/Vp)
Denne formel danner grundlaget for at bestemme den radius, inden for hvilken jordskælvet må have fundet sted i forhold til den enkelte station.
Geometrisk triangulering og mindste kvadraters metode
En enkelt målestation giver kun en distance, hvilket definerer en cirkel på Jordens overflade, hvor epicentret kan ligge. To stationer danner to cirkler, som typisk skærer hinanden to steder, hvilket indsnævrer de mulige placeringer til to punkter. For at bestemme et unikt epicenter er det nødvendigt med mindste tre stationer.
I den virkelige verden vil målinger og tidsregistreringer altid indeholde små unøjagtigheder, hvilket betyder, at cirklerne fra tre eller flere stationer sjældent skærer hinanden i ét præcist punkt. For at løse dette problem anvender værktøjet en mindste kvadraters optimering. Metoden minimerer summen af de kvadrerede afstandsafvigelser på tværs af alle aktive stationer:
minimer Σ(dᵢ − Dᵢ)² på tværs af alle stationer
Her repræsenterer dᵢ den sfæriske overfladeafstand fra det estimerede epicenter til stationen, mens Dᵢ repræsenterer den distance, der er beregnet ud fra stationens S−P-interval. Jordkloden modelleres her som en perfekt kugle med en konstant radius R på 6.371 km. Overfladeafstanden beregnes ved hjælp af haversinus-formlen:
d = R × haversinus-centralvinkel, med R = 6371 km
Indtastningsparametre og valideringsregler
For at udføre beregningen skal brugeren angive en række specifikke inputværdier. Værktøjet validerer alle indtastninger i henhold til følgende grænser og regler:
Stationsobservationer
Der skal tilføjes mellem 3 og 8 målestationer. Hvis der forsøges beregnet med færre end 3 eller flere end 8 stationer, vises fejlmeddelelsen: Brug mellem tre og otte stationer.. For hver station angives:
- Stationsnavn
- Breddegrad: Skal være en numerisk værdi fra −90° til 90°. Ved overtrædelse vises en fejlmeddelelse, hvor stationsnavnet indsættes dynamisk, f.eks.:
Station 1: Breddegrad skal være fra −90° til 90°.. - Længdegrad: Skal være en numerisk værdi fra −180° til 180°. Ved overtrædelse vises en fejlmeddelelse, hvor stationsnavnet indsættes dynamisk, f.eks.:
Station 1: Længdegrad skal være fra −180° til 180°.. - S−P-interval: Skal være et tal større end 0 og højst 2.000 sekunder. Ved overtrædelse vises en fejlmeddelelse, hvor stationsnavnet indsættes dynamisk, f.eks.:
Station 1: Indtast et S−P-interval, der er større end 0 og højst 2000 sekunder..
Udbredelsesmodel
- P-wave-hastighed: Skal være større end 0 og højst 20 km/s.
- S-wave-hastighed: Skal være større end 0 og højst 20 km/s.
- Tidsmæssig usikkerhed (1σ): Skal være større end 0 og højst 60 sekunder. Ved overtrædelse vises:
Tidsusikkerheden skal være større end 0 og højst 60 sekunder..
Yderligere valideringsregler
- Ikke-numeriske værdier: Hvis et felt indeholder ugyldige tegn, vises en fejlmeddelelse, der angiver placeringen og den indtastede værdi som dynamiske parametre:
input: ”abc” er ikke et tal.. - Hastighedsforhold: P-wave-hastigheden skal altid være højere end S-wave-hastigheden. Hvis dette ikke er tilfældet, vises:
P-wave-hastighed skal være større end S-wave-hastighed.. Hvis en af hastighederne er uden for grænserne, vises en fejlmeddelelse, hvor feltnavnet indsættes dynamisk:Vₚ/Vₛ skal være større end 0 og højst 20 km/s.. - Identiske koordinater: To stationer må ikke have nøjagtig samme placering. Hvis det sker, vises en fejlmeddelelse, hvor de berørte stationsnavne indsættes dynamisk:
Stationerne 1 og 2 har de samme koordinater. Brug unikke stationsplaceringer.. - Ustabil geometri: Hvis stationerne ligger på en næsten ret linje eller i en meget tæt klynge, kan en todimensionel placering ikke bestemmes stabilt. Værktøjet afbryder beregningen med meddelelsen:
Disse stationer afgrænser ikke en stabil todimensionel placering. Tilføj en station væk fra den nuværende linje eller klynge.. - Numerisk overløb: Hvis beregningen overskrider de tekniske grænser, vises:
En værdi eller et mellemresultat overskrider det understøttede talområde..
Brugere kan klikke på Indlæs Tokyo-eksempel for at udfylde værktøjet med prædefinerede eksempeldata, eller klikke på Ryd for at nulstille alle felter.
Fortolkning af resultater og usikkerhed
Når beregningen er gennemført, præsenterer værktøjet en række resultater, der beskriver løsningens kvalitet og nøjagtighed:
- Estimeret epicenter: De beregnede koordinater for jordskælvets epicenter.
- Distancetilpasning: Kvaliteten af afstandsmatchet kategoriseres som enten tæt distancetilpasning, blandet distancetilpasning eller dårlig distancetilpasning.
- Stationsgeometri: Kvaliteten af stationernes geografiske fordeling i forhold til epicentret kategoriseres som stærk stationsgeometri, moderat stationsgeometri eller svag stationsgeometri.
- Residualværdier: For hver station beregnes en residualværdi, som er den tilpassede sfæriske distance minus den distance, der er udledt af S−P-intervallet. En positiv residualværdi betyder, at det beregnede epicenter ligger længere væk fra stationen end den radius, S−P-intervallet foreskriver. En negativ residualværdi betyder, at det ligger tættere på.
- Fejlforplantning: Værktøjet beregner, hvordan den tidsmæssige usikkerhed (1σ) forplanter sig gennem stationsgeometrien til en rumlig usikkerhed langs den primære horisontale akse. Dette rapporteres som en omtrentlig positions-1σ i kilometer. Denne usikkerhed inkluderer udelukkende den tidsmæssige usikkerhed og tager ikke højde for systematiske fejl som forkerte hastighedsmodeller eller faseidentifikationsfejl.
En skematisk fremstilling viser den relative geometri af stationerne og det beregnede epicenter. Denne grafik ledsages af teksten: "Kun relativ geometri – intet baggrundskort, grænser eller live jordskælvsdata.".
Pædagogisk model vs. virkelige seismiske kataloger
Dette værktøj er designet som en undervisningsmodel og anvender en række forenklinger af hensyn til overskuelighed og beregningshastighed. Det er vigtigt at forstå forskellene mellem denne model og professionelle seismiske lokaliseringer:
| Egenskab | Denne pædagogiske model | Professionelle seismiske kataloger |
|---|---|---|
| Jordens form | Perfekt sfærisk (R = 6.371 km) | Ellipsoidisk (f.eks. WGS84) |
| Bølgehastighed | Konstant og ensartet i alle retninger | Variabel, baseret på lagdelte eller 3D-hastighedsmodeller |
| Kilde-dybde | Antages at være en lav kilde (overfladenær) | Beregner både epicenter og hypocenterdybde (fokusdybde) |
| Bølgeudbredelse | Lineære overfladedistancer | Krumme bølgeveje gennem Jordens indre |
Metoden og formlerne i værktøjet er afstemt med pædagogiske ressourcer fra USGS, herunder USGS: Videnskaben om jordskælv, USGS: Seismiske bølger og lokalisering af jordskælv og USGS: Rejsetider for jordskælv.
Ofte stillede spørgsmål (FAQ)
Hvorfor er der brug for tre stationer?
Én S−P-måling giver en distancecirkel, ikke en retning. To cirkler kan mødes to steder. En tredje station afgør normalt valget, mens en fjerde eller efterfølgende station gør uoverensstemmelsen synlig via residualer. Stationer, der er spredt rundt om kilden, afgrænser både breddegrad og længdegrad bedre end stationer på en enkelt linje.
Er S−P-distancen den nøjagtige overfladedistance?
Nej. Ligningen med konstant hastighed estimerer en fælles udbredelsesvej ud fra forskellen i ankomsttid. At behandle denne radius som sfærisk overfladedistance forudsætter en lav kilde og en simpel hastighedsstruktur. Reelle løbetidskurver afhænger af dybde, distance og det materiale, som hver bølge passerer igennem.
Hvad betyder et positivt eller negativt residual?
Residualet er den tilpassede sfæriske distance minus den distance, der er beregnet ud fra S−P. Positiv betyder, at det tilpassede punkt er længere væk fra den pågældende station end dens S−P-radius; negativ betyder, at det er tættere på. Store, blandede residualer viser, at cirklerne ikke har ét rent skæringspunkt under de valgte hastigheder.
Hvad er inkluderet i positionsusikkerheden?
Den propagerer den enkelte tidsmæssige 1σ-værdi gennem S−P-distancefaktoren og stationsgeometrien nær det tilpassede punkt. Den inkluderer ikke en forkert hastighedsmodel, fokusdybde, fejlbehæftede faser, ur-afvigelser eller tredimensionel jordstruktur, så det er ikke en komplet nøjagtighedserklæring.