Фізика сейсмічних хвиль та інтервал S−P
Під час землетрусу вивільнена енергія поширюється крізь надра Землі у вигляді сейсмічних хвиль різних типів. Першими до реєструвальних станцій доходять первинні хвилі (P-хвилі), які є поздовжніми хвилями стиснення та розтягнення. За ними слідують вторинні хвилі (S-хвилі) — поперечні хвилі зсуву, що мають меншу швидкість поширення. Оскільки P-хвилі рухаються швидше за S-хвилі, часова затримка між їхнім приходом збільшується пропорційно до відстані від джерела землетрусу до сейсмічної станції.
Ця різниця в часі реєстрації фаз називається інтервалом S−P. Використовуючи спрощену модель однорідної кори зі сталими швидкостями хвиль, зв’язок між часовим інтервалом та відстанню описується рівнянням:
Δ t = D × (1 / Vₛ - 1 / Vₚ)
Шляхом математичного перетворення отримуємо формулу для обчислення відстані від станції до джерела:
D = Δ t ÷ (1 / Vₛ - 1 / Vₚ)
Де:
- D — відстань до джерела землетрусу;
- Δ t — інтервал часу між приходом хвиль S та P;
- Vₚ — швидкість P-хвилі;
- Vₛ — швидкість S-хвилі.
Геометрія тріангуляції на сфері
Для визначення географічного положення епіцентру використовується метод тріангуляції. Одне вимірювання інтервалу S−P на поодинокій станції дає лише радіус віддалення, окреслюючи коло можливих положень джерела на поверхні Землі. Дані двох станцій обмежують область пошуку до двох точок перетину їхніх кіл. Для однозначного визначення координат епіцентру на площині чи сфері необхідні спостереження щонайменше трьох станцій.
Калькулятор моделює Землю як сферу зі сталим радіусом R = 6371 км. Розрахунок відстані між точками на сфері базується на тригонометричних співвідношеннях:
Відстань на поверхні: d = R × гаверсинус центрального кута, де R = 6371 km
У реальних умовах через похибки вимірювань та неоднорідність середовища кола відстаней від трьох або більше станцій рідко перетинаються в одній точці. Тому для пошуку оптимальних координат застосовується математичне узгодження методом найменших квадратів. Алгоритм мінімізує суму квадратів різниць між теоретичними сферичними відстанями та відстанями, обчисленими за інтервалами S−P для всіх станцій:
Узгодження епіцентру: мінімізація Σ(dᵢ - Dᵢ)² для всіх станцій
Вхідні параметри та обмеження моделі
Для проведення розрахунків користувач має задати параметри середовища та дані спостережень для групи від 3 до 8 станцій.
Параметри моделі поширення:
- Швидкість P-хвилі: має бути більше 0 і не більше 20 km/s. Також вона обов'язково повинна перевищувати швидкість S-хвилі.
- Швидкість S-хвилі: має бути більше 0 і не більше 20 km/s, а її значення має бути меншим за швидкість P-хвилі.
- Похибка часу (1σ): випадкова похибка визначення моменту приходу фаз, що має бути більше 0 і не більше 60 секунд.
Дані спостережень станцій:
- Назва станції (текстовий ідентифікатор).
- Широта: географічна координата в межах від −90° до 90°.
- Довгота: географічна координата в межах від −180° до 180°.
- Інтервал S−P: тривалість затримки в секундах, що має бути більше 0 і не більше 2000 секунд.
Для швидкого ознайомлення з роботою алгоритму можна натиснути кнопку Завантажити приклад Токіо, яка автоматично заповнить поля тестовими даними, або очистити форму за допомогою кнопки Очистити.
Аналіз результатів та оцінка похибок
Після обробки введених даних калькулятор надає детальний звіт, що містить оцінку якості розв'язку та геометричні характеристики мережі спостережень.
Якість узгодження відстаней та геометрії
Точність визначення епіцентру залежить не лише від точності вимірювання часу, а й від взаємного розташування станцій відносно джерела події. Інструмент класифікує ці параметри за такими категоріями:
- Якість узгодження відстаней: набуває значень «тісне узгодження відстаней», «змішане узгодження відстаней» або «слабке узгодження відстаней».
- Якість геометрії станцій: класифікується як «сильна геометрія станцій», «помірна геометрія станцій» або «слабка геометрія станцій».
Статистичні показники та нев'язки
Для кожної станції у таблиці результатів розраховується індивідуальна нев'язка — різниця між сферичною відстанню до підібраного епіцентру та радіусом, обчисленим за часом S−P. Додатково виводяться інтегральні показники:
- Середньоквадратична нев’язка відстані (RMS) — узагальнена міра відхилення розв'язку.
- Найбільша нев’язка відстані серед усіх задіяних станцій.
- Приблизна похибка положення 1σ та Похибка відстані 1σ на станцію, які показують просторову невизначеність розв'язку на основі заданої похибки часу.
Обробка помилок та валідація даних
Калькулятор містить жорсткі правила перевірки введених даних для запобігання математичній нестабільності:
- Кількість станцій: якщо спробувати виконати розрахунок для кількості станцій поза межами діапазону 3–8, з'явиться повідомлення:
Використовуйте від трьох до восьми станцій.Швидкість P-хвилі має бути більшою за швидкість S-хвилі.- Якщо станції розташовані на одній лінії або занадто тісно:
Ці станції не забезпечують стабільного двовимірного визначення місця. Додайте станцію осторонь від поточної лінії або скупчення.
Конфіденційність обробки даних
Усі обчислення виконуються безпосередньо у веб-браузері користувача. Введені координати станцій, часові інтервали та результати розрахунків не передаються на зовнішні сервери й залишаються на локальному пристрої.
Часті запитання (FAQ)
Чи є відстань за S−P точною відстанню на поверхні?
Ні. Рівняння зі сталою швидкістю оцінює спільний шлях проходження за різницею часу приходу хвиль. Розгляд цього радіуса як сферичної відстані на поверхні припускає мілке джерело та просту структуру швидкостей. Реальні годографи залежать від глибини, відстані та середовища, через яке проходить кожен промінь.
Що означає позитивна чи негативна нев’язка?
Нев’язка — це узгоджена сферична відстань мінус відстань, обчислена за S−P. Позитивне значення означає, що узгоджена точка розташована далі від цієї станції, ніж її радіус за S−P; негативне — що вона ближче. Великі різноспрямовані нев’язки показують, що кола не мають одного чіткого перетину за вибраних швидкостей.
Чому потрібні три станції?
Одне вимірювання S−P дає коло відстані, а не напрямок. Два кола можуть перетинатися у двох точках. Третя станція зазвичай дозволяє зробити однозначний вибір, тоді як четверта та наступні станції роблять розбіжність видимою через нев’язки. Станції, розподілені навколо джерела, обмежують як широту, так і довготу краще, ніж станції вздовж однієї лінії.
Що враховує похибка визначення положення?
Вона поширює єдине значення похибки часу 1σ через коефіцієнт відстані S−P та геометрію станцій поблизу узгодженої точки. Вона не враховує неправильну швидкісну модель, глибину фокусу, помилково визначені фази, зміщення годинника або тривимірну структуру Землі, тому не є повною оцінкою точності.