Физика сеизмичких таласа и S−P интервал
Лоцирање епицентра земљотреса заснива се на чињеници да се различити сеизмички таласи крећу кроз Земљину кору различитим брзинама. Примарни или компресиони таласи (P-таласи) су најбржи и први стижу до сеизмолошке станице. Секундарни или попречни таласи (S-таласи) путују спорије и стижу касније. Како се удаљеност од извора земљотреса повећава, временски размак између доласка ове две групе таласа — познат као S−P интервал — постаје све већи.
Под претпоставком константне брзине простирања у једноличном моделу коре, овај временски интервал омогућава директну процену удаљености од станице до извора. Однос између времена и растојања дефинисан је формулом:
Апроксимација растојања: D = Δt ÷ (1/Vs − 1/Vp)
У овој једначини, Δt представља унесени S−P интервал, док су Vp и Vs брзине P-таласа и S-таласа које дефинише корисник. Овај прорачун претпоставља плитак извор земљотреса и занемарује сложене факторе као што су закривљене путање зрака кроз дубље слојеве Земље.
Геометрија тријангулације на сфери
Подаци са једне сеизмолошке станице пружају само информацију о удаљености, што значи да се епицентар може налазити било где на кружници чији је полупречник једнак израчунатом растојању. Када имамо две станице, њихови кругови се пресецају у највише две тачке. Тек са три станице добијамо трећи круг који омогућава једнозначно одређивање положаја епицентра.
За прорачуне на закривљеној површини, Земља се моделује као сфера са константним полупречником:
Површинско растојање: d = R × хаверсин централног угла, са R = 6371 km
У реалним условима, услед грешака у мерењу времена и локалних варијација у брзини таласа, кругови се ретко пресецају у једној савршеној тачки. Због тога се користи математички поступак оптимизације који тражи најбоље решење за све укључене станице:
Уклапање епицентра: минимизовати Σ(dᵢ − Dᵢ)² за све станице
Овај метод најмањих квадрата проналази координате које минимизују разлике између сферних растојања на моделу Земље (dᵢ) и растојања изведених из S−P интервала (Dᵢ) за све станице истовремено.
Унос података и ограничења модела
Калкулатор омогућава унос података за најмање 3 и највише 8 сеизмолошких станица. За сваку станицу потребно је дефинисати следеће параметре:
- Назив станице
- Географска ширина: Вредност мора бити у опсегу од −90° до 90°.
- Географска дужина: Вредност мора бити у опсегу од −180° до 180°.
- S−P интервал: Вредност мора бити већа од 0 и не већа од 2000 секунди.
У оквиру секције Модел пропагације, корисник подешава физичке параметре средине:
- Брзина P-таласа: Мора бити већа од 0 и не већа од 20 km/s, и мора бити већа од брзине S-таласа.
- Брзина S-таласа: Мора бити већа од 0 и не већа од 20 km/s, и мора бити мања од брзине P-таласа.
- Несигурност времена (1σ): Мора бити већа од 0 и не већа од 60 секунди.
Корисник може брзо тестирати рад алата кликом на дугме Учитај пример за Токио, које аутоматски попуњава поља тестним подацима, или ресетовати све уносе кликом на дугме Очисти.
Анализа резултата и одступања
Након покретања прорачуна, алат приказује детаљне резултате у неколико секција. У главном блоку Решење епицентра приказују се израчунате координате под ознаком Процењени епицентар. Квалитет добијеног решења се оцењује кроз неколико категорија:
- Квалитет уклапања растојања: Приказује се као блиско уклапање растојања, мешовито уклапање растојања или лоше уклапање растојања.
- Квалитет геометрије станица: Описује просторни распоред станица око епицентра и може бити јака геометрија станица, умерена геометрија станица или слаба геометрија станица.
- Статистички индикатори: Приказују се вредности за RMS резидуал растојања, Највећи резидуал растојања, Приближна несигурност позиције 1σ и Несигурност растојања 1σ по станици.
Табела под називом Уклапање станица садржи следеће колоне за сваку појединачну станицу:
- Станица
- S−P растојање (израчунати полупречник на основу времена доласка)
- Сферно растојање (стварна удаљеност од уклопљеног епицентра на сфери)
- Резидуал (разлика између сферног и S−P растојања)
Поред табеле, Шематски приказ геометрије станица визуелизује међусобни положај станица и епицентра, уз напомену: „Само релативна геометрија — без мапе, граница или података о земљотресима у реалном времену.”
Математички извод и замена вредности
У одељку Формула и замена приказан је детаљан ток прорачуна са конкретним вредностима које је корисник унео:
- Фактор растојања: Приказује се израчунати фактор конверзије времена у удаљеност:
Фактор растојања = vp × vs ÷ (vp − vs) = factor km по секунди S−P - Прорачун по станицама: За сваку станицу се приказује замена вредности:
station: interval s × factor km/s = distance km - Најбоље уклапање: Приказују се коначне координате и средње квадратна грешка:
Најбоље сферно уклапање = latitude°, longitude°; RMS резидуал = rms km - Пропагација грешке: Приказује се како се несигурност мерења времена преноси на просторну несигурност:
Несигурност времена 1σ од timing s даје несигурност растојања 1σ од distance km; геометрија станица ово преноси на око position km на главној хоризонталној оси.
Правила валидације и грешке при уносу
Алат врши строгу проверу свих унетих вредности како би се обезбедила стабилност нумеричког алгоритма. У случају неправилних уноса, приказују се следеће поруке:
- Ако број станица није у дозвољеном опсегу:
Користите између три и осам станица. - Ако унесена вредност није број:
where: „token“ није број. - Ако је географска ширина ван опсега:
Станица station: географска ширина мора бити од −90° до 90°. - Ако је географска дужина ван опсега:
Станица station: географска дужина мора бити од −180° до 180°. - Ако је временски интервал ван граница:
Станица station: унесите S−P интервал већи од 0 и не већи од 2000 секунди. - Ако су брзине таласа ван опсега:
field мора бити већа од 0 и не већа од 20 km/s. - Ако је редослед брзина погрешан:
Брзина P-таласа мора бити већа од брзине S-таласа. - Ако је несигурност времена ван граница:
Несигурност времена мора бити већа од 0 и не већа од 60 секунди. - Ако две станице имају идентичне координате:
Станице station и other имају исте координате. Користите различите локације станица. - Ако су станице постављене у линију или преблизу једна другој:
Ове станице не ограничавају стабилну дводимензионалну локацију. Додајте станицу ван тренутне линије или групе. - У случају прекорачења опсега током рачунања:
Вредност или међурезултат прелази подржани опсег бројева.
Сигурност података и обрада
Сви прорачуни се обављају локално, директно у веб прегледачу корисника. Координате станица, S−P временски интервали и сви добијени резултати остају искључиво на локалном уређају и не отпремају се на спољне сервере.
Често постављана питања (FAQ)
Да ли је S−P растојање тачно површинско растојање?
Не. Једначина са константном брзином процењује заједничку путању путовања на основу разлике у времену доласка. Третирање тог полупречника као сферног површинског растојања подразумева плитак извор и једноставну структуру брзина. Стваран ток времена путовања зависи од дубине, удаљености и материјала кроз који сваки зрак пролази.
Шта значи позитиван или негативан резидуал?
Резидуал је уклопљено сферно растојање минус растојање процењено на основу S−P. Позитивна вредност означава да је уклопљена тачка удаљенија од те станице него њен израчунати S−P полупречник. Негативна вредност значи да је уклопљена тачка ближа станици. Велики мешовити резидуали показују да се кругови појединачних станица не пресецају у једној тачки при изабраним брзинама таласа.
Зашто су потребне три станице?
Једно S−P мерење дефинише само круг могућих локација, без смера. Два круга се најчешће пресецају у две тачке. Трећа станица је неопходна да би се одредио јединствен епицентар, док додатне станице (четврта и свака следећа) омогућавају процену грешке кроз резидуале. Такође, станице које су правилно распоређене око извора пружају знатно стабилније решење од оних које леже у правој линији.
Шта је укључено у несигурност позиције?
Она преноси појединачну вредност несигурности времена 1σ кроз фактор S−P растојања и геометрију станица у близини уклопљене тачке. Не укључује нетачан модел брзине, дубину жаришта, погрешно идентификоване фазе, одступање сата станице нити тродимензионалну структуру Земље, тако да не представља потпуну тврдњу о тачности.