Fizica undelor seismice și intervalul S−P
Localizarea epicentrului unui cutremur se bazează pe diferența de viteză dintre diferitele tipuri de unde seismice generate în momentul producerii unei rupturi în scoarța terestră. Undele P (primare sau de compresiune) sunt cele mai rapide unde de volum și se propagă primele prin mediul terestru. Undele S (secundare sau de forfecare) sunt mai lente și ajung la stațiile de înregistrare după undele P. Această diferență de viteză face ca intervalul de timp dintre sosirea undei P și sosirea undei S (cunoscut ca intervalul S−P) să crească proporțional cu distanța parcursă de la sursa seismului.
O singură stație seismică ce înregistrează un anumit interval S−P poate determina doar distanța până la sursă, definind o rază de căutare sub forma unui cerc pe suprafața terestră. Pentru a localiza precis epicentrul, este necesară combinarea datelor de la mai multe stații prin tehnica triangulației.
Geometria triangulației și optimizarea prin cele mai mici pătrate
Atunci când se utilizează date de la o singură stație, localizarea este limitată la un cerc de distanță. Două stații seismice generează două cercuri care, în mod tipic, se intersectează în două puncte distincte. O a treia stație este esențială pentru a rezolva această ambiguitate și a indica un singur punct de intersecție comun.
În realitate, din cauza erorilor de citire a timpilor de sosire și a variațiilor din structura geologică, cercurile de distanță de la trei sau mai multe stații nu se intersectează niciodată perfect într-un singur punct. Din acest motiv, se utilizează o metodă de optimizare prin cele mai mici pătrate pentru a calcula o soluție de tip "cea mai bună potrivire". Această metodă matematică minimizează suma pătratelor diferențelor dintre distanțele sferice calculate de la epicentrul estimat la stații și distanțele teoretice deduse din intervalele S−P înregistrate.
Parametrii de intrare și configurarea modelului
Pentru a rula calculul de triangulație, utilizatorul trebuie să introducă date pentru un număr de minimum 3 și maximum 8 stații seismice. Pentru fiecare stație sunt necesare următoarele detalii:
- Nume stație
- Latitudine: Trebuie să fie o valoare numerică cuprinsă între −90° și 90°.
- Longitudine: Trebuie să fie o valoare numerică cuprinsă între −180° și 180°.
- Interval S−P: Timpul exprimat în secunde, care trebuie să fie mai mare de 0 și de cel mult 2000 secunde.
De asemenea, modelul de propagare necesită definirea următorilor parametri fizici:
- Viteza undei P: Trebuie să fie mai mare de 0 și de cel mult 20 km/s, depășind obligatoriu viteza undei S.
- Viteza undei S: Trebuie să fie mai mare de 0 și de cel mult 20 km/s, fiind obligatoriu mai mică decât viteza undei P.
- Incertitudine de timp (1σ): Valoarea abaterii standard a timpului de înregistrare, care trebuie să fie mai mare de 0 și de cel mult 60 de secunde.
Interfața oferă opțiunea Încărcare exemplu Tokyo pentru a popula automat câmpurile cu un set de date demonstrativ, precum și opțiunea Ștergere pentru a reseta toate intrările.
Formule matematice și aproximări sferice
Calculatorul utilizează un model simplificat în care Pământul este considerat o sferă perfectă cu o rază constantă R = 6371 km. Calculul distanțelor și optimizarea poziției se realizează pe baza următoarelor relații matematice:
- Aproximarea distanței: Distanța teoretică de la stație la sursă (D) este calculată direct din intervalul de timp S−P (Δ t) și vitezele undelor (Vₚ și Vₛ): „Aproximarea distanței: D = Δt ÷ (1/Vs − 1/Vp)”
- Distanța la suprafață: Distanța reală pe sferă (d) între coordonatele stației și epicentrul estimat se calculează folosind formula haversine pentru unghiul central: „Distanța la suprafață: d = R × unghiul central haversine, cu R = 6371 km”
- Potrivirea epicentrului: Algoritmul caută coordonatele geografice care minimizează suma reziduurilor pătratice pentru toate cele i stații active: „Potrivirea epicentrului: minimizarea Σ(dᵢ − Dᵢ)² pentru toate stațiile”
Această metodă didactică presupune un cutremur cu vatră de adâncime mică (superficial) și viteze constante ale undelor în întreaga crustă, omițând astfel adâncimea focală reală, căile curbate ale razelor seismice prin straturile interne ale Pământului și structurile tridimensionale de viteză.
Analiza rezultatelor și interpretarea erorilor
După rularea algoritmului, sistemul generează o serie de indicatori de calitate și date de analiză structurală:
- Rază calculată: Distanța estimată pentru fiecare stație pe baza timpului S−P.
- Soluție epicentru: Coordonatele geografice calculate pentru Epicentru estimat.
- Calitatea potrivirii: Evaluată vizual și etichetată drept potrivire strânsă a distanței, potrivire mixtă a distanței sau potrivire slabă a distanței.
- Calitatea geometriei: Evaluată ca geometrie puternică a stațiilor, geometrie moderată a stațiilor sau geometrie slabă a stațiilor.
- Reziduu de distanță RMS și Cel mai mare reziduu de distanță: Indică precizia globală a potrivirii matematice.
- Incertitudine poziție aproximativă 1σ și Distanță 1σ per stație: Reflectă propagarea erorii de timp introduse de utilizator prin geometria rețelei de stații.
Tabelul de potrivire a stațiilor (Potrivirea stațiilor) conține coloanele: Stație, Distanță S−P, Distanță sferică și Reziduu. Reziduul reprezintă diferența dintre distanța sferică potrivită și distanța dedusă din intervalul S−P. Un reziduu pozitiv indică faptul că epicentrul calculat se află mai departe de stație decât raza S−P, în timp ce un reziduu negativ arată că acesta este mai aproape.
Secțiunea Formulă și substituție detaliază pașii de calcul, afișând mesaje dinamice precum:
- Calculul factorului de distanță pe baza vitezelor undelor P și S introduse, exprimat în kilometri pe secundă de S−P.
- Distanța calculată pentru fiecare stație, înmulțind intervalul S−P înregistrat cu factorul de distanță determinat.
- Coordonatele celei mai bune potriviri sferice (latitudine și longitudine) și valoarea reziduului RMS în kilometri.
- Propagarea incertitudinii de timp 1σ introduse pentru a obține incertitudinea de distanță 1σ și transpunerea acesteia prin geometria stațiilor într-o incertitudine de poziție aproximativă pe axa orizontală majoră.
Validarea datelor și mesaje de eroare
Pentru a asigura stabilitatea matematică a modelului, calculatorul aplică reguli stricte de validare a datelor introduse:
- Număr de stații: Dacă se încearcă rularea cu mai puțin de 3 sau mai mult de 8 stații, se afișează eroarea:
Utilizați între trei și opt stații. - Valori non-numerice: Dacă un câmp conține caractere nepermise, se afișează:
input: „abc” nu este un număr. - Limite coordonate:
Stația 1: latitudinea trebuie să fie între −90° și 90°.Stația 1: longitudinea trebuie să fie între −180° și 180°.
- Interval S−P:
Stația 1: introduceți un interval S−P mai mare de 0 și de cel mult 2000 secunde. - Viteze unde:
Vₚ/Vₛ trebuie să fie mai mare de 0 și de cel mult 20 km/s.Viteza undei P trebuie să fie mai mare decât viteza undei S.
- Incertitudine de timp:
Incertitudinea de timp trebuie să fie mai mare de 0 și de cel mult 60 de secunde. - Coordonate duplicate:
Stațiile 1 și 2 au aceleași coordonate. Utilizați locații distincte pentru stații. - Geometrie instabilă: Dacă stațiile sunt coliniare sau extrem de grupate, algoritmul nu poate defini o poziție bidimensională stabilă și afișează:
Aceste stații nu constrâng o locație bidimensională stabilă. Adăugați o stație în afara liniei sau grupului actual. - Depășire numerică:
O valoare sau un rezultat intermediar depășește intervalul de numere acceptat.
Confidențialitate și procesare locală
Toate calculele matematice, coordonatele geografice introduse și timpii S−P asociați sunt procesați exclusiv local, în browserul web al utilizatorului. Datele nu sunt încărcate pe servere externe și nu părăsesc dispozitivul dumneavoastră.
Întrebări frecvente (FAQ)
De ce sunt necesare trei stații? O singură măsurătoare S−P oferă un cerc de distanță, nu o direcție. Două cercuri se pot întâlni în două locuri. O a treia stație rezolvă în mod normal alegerea, în timp ce o a patra sau următoarele stații fac vizibilă nepotrivirea prin intermediul reziduurilor. Stațiile răspândite în jurul sursei constrâng atât latitudinea, cât și longitudinea mai bine decât stațiile situate de-a lungul unei singure linii.
Este distanța S−P distanța exactă la suprafață? Nu. Ecuația cu viteză constantă estimează o cale de deplasare comună pe baza diferenței timpului de sosire. Tratarea acelei raze ca distanță sferică la suprafață presupune o sursă de adâncime mică și o structură simplă a vitezei. Curbele reale ale timpului de parcurs depind de adâncime, distanță și materialul traversat de fiecare rază.
Ce înseamnă un reziduu pozitiv sau negativ? Reziduul este distanța sferică potrivită minus distanța dedusă din S−P. Pozitiv înseamnă că punctul potrivit este mai departe de acea stație decât raza sa S−P; negativ înseamnă că este mai aproape. Reziduurile mari mixte arată că cercurile nu împart o singură intersecție curată la vitezele alese.
Ce este inclus în incertitudinea poziției? Aceasta propagă valoarea unică de timp 1σ prin factorul de distanță S−P și geometria stațiilor din apropierea punctului potrivit. Nu include un model de viteză incorect, adâncimea focală, fazele etichetate greșit, decalajul ceasului sau structura tridimensională a Pământului, deci nu reprezintă o afirmație completă de precizie.