Физика сейсмических волн и интервал S−P
При возникновении землетрясения энергия высвобождается в виде сейсмических волн, которые распространяются сквозь толщу Земли. Первыми регистрируются продольные волны, или P-волны (первичные), которые совершают колебания вдоль направления распространения. За ними следуют поперечные волны, или S-волны (вторичные), колеблющиеся перпендикулярно направлению движения. Поскольку P-волны движутся быстрее, чем S-волны, разница во времени их прибытия на сейсмическую станцию увеличивается по мере удаления от источника.
Этот промежуток времени называется интервалом S−P. Измерив его, можно вычислить расстояние от станции до эпицентра, используя скорости распространения волн в земной коре. В рамках калькулятора применяется допущение о постоянной скорости волн, при котором разность времени прибытия удовлетворяет уравнению:
Δt = D × (1/Vs − 1/Vp)
Где Δt — интервал S−P, Vs — скорость S-волны, Vp — скорость P-волны, а D — расстояние до источника. Путем математического преобразования вычисляется расстояние от станции до эпицентра:
D = Δt ÷ (1/Vs − 1/Vp)
Геометрия триангуляции на сфере
Данные одной сейсмической станции позволяют определить только расстояние до источника, очерчивая на поверхности Земли окружность возможного положения эпицентра. Две станции дают две пересекающиеся окружности, что сужает область поиска до двух возможных точек. Для однозначного определения координат эпицентра требуются данные как минимум трех станций.
Калькулятор моделирует Землю как сферу с постоянным радиусом R = 6371 км. Расстояние по поверхности между точкой на сфере и станцией рассчитывается через гаверсинус центрального угла:
d = R × гаверсинус центрального угла, где R = 6371 km
Для расчета используются географические координаты станций: широта (в пределах от −90° до 90°) и долгота (в пределах от −180° до 180°).
Метод наименьших квадратов для подбора эпицентра
В реальных условиях из-за погрешностей измерений и неоднородности земной коры окружности расстояний от трех или более станций редко пересекаются в одной точке. Для нахождения оптимального решения применяется математический метод наименьших квадратов. Калькулятор минимизирует сумму квадратов разностей между расчетными сферическими расстояниями от станций до оцениваемой точки и расстояниями, полученными на основе интервалов S−P:
минимизация Σ(dᵢ − Dᵢ)² по всем станциям
Этот подход позволяет использовать избыточные данные (от 3 до 8 станций) для повышения точности локализации и оценки качества исходных наблюдений.
Входные параметры и ограничения модели
Для проведения расчета пользователю необходимо задать параметры модели распространения волн и данные наблюдений станций.
Параметры модели распространения
- Скорость P-волны: должна быть больше 0 и не более 20 km/s. Также она должна строго превышать скорость S-волны.
- Скорость S-волны: должна быть больше 0 и не более 20 km/s, а также быть меньше скорости P-волны.
- Погрешность времени (1σ): должна быть больше 0 и не более 60 секунд.
Данные наблюдений станций (от 3 до 8 станций)
- Широта: от −90° до 90°.
- Долгота: от −180° до 180°.
- Интервал S−P: больше 0 и не более 2000 секунд.
Для быстрой проверки работы калькулятора можно нажать кнопку Загрузить пример Токио, которая заполнит поля демонстрационными данными, или очистить форму с помощью кнопки Очистить.
Ограничения и ошибки ввода
Калькулятор выполняет строгую проверку вводимых данных и выводит соответствующие предупреждения:
- При попытке расчета с числом станций менее 3 или более 8:
Используйте от трех до восьми станций. - Нарушение соотношения скоростей:
Скорость P-волны должна быть больше скорости S-волны. - Погрешность времени вне диапазона:
Погрешность времени должна быть больше 0 и не более 60 секунд. - Нестабильное расположение станций (например, на одной линии):
Эти станции не позволяют определить стабильное двухмерное положение. Добавьте станцию в стороне от текущей линии или группы. - Превышение разрядности вычислений:
Значение или промежуточный результат превышает поддерживаемый диапазон чисел.
Анализ результатов и оценка погрешности
Приближение расстояния: D = Δt ÷ (1/Vs − 1/Vp)
Расстояние по поверхности: d = R × гаверсинус центрального угла, где R = 6371 km
Подбор эпицентра: минимизация Σ(dᵢ − Dᵢ)² по всем станциям
После выполнения расчета калькулятор выводит детальный отчет в разделе Решение для эпицентра:
- Оцененный эпицентр: рассчитанные координаты широты и долготы.
- Соответствие станций: таблица со столбцами Станция, Расстояние по S−P, Сферическое расстояние и Невязка. Невязка рассчитывается как сферическое расстояние минус расстояние по S−P. Положительная невязка означает, что подобранный эпицентр находится дальше от станции, чем расчетный радиус S−P, а отрицательная — ближе.
- Качество соответствия расстояний: хорошее соответствие расстояний, смешанное соответствие расстояний или плохое соответствие расстояний.
- Геометрия станций: сильная геометрия станций, умеренная геометрия станций или слабая геометрия станций.
- Статистические показатели: Среднеквадратичная невязка расстояния, Наибольшая невязка расстояния, Приблизительная погрешность положения 1σ и Погрешность расстояния 1σ на станцию.
Визуальная схема расположения станций отображается в блоке Схема геометрии станций с обязательной пометкой: «Только относительная геометрия — без карты, границ и реальных данных о землетрясениях.»
Конфиденциальность вычислений
Все расчеты выполняются непосредственно в веб-браузере на устройстве пользователя. Введенные координаты станций, интервалы времени и результаты вычислений не передаются на внешние серверы, что обеспечивает конфиденциальность обрабатываемых данных.
Часто задаваемые вопросы
Является ли расстояние по S−P точным расстоянием по поверхности?
Нет. Уравнение с постоянной скоростью оценивает общий путь прохождения волны по разности времени прибытия. Рассмотрение этого радиуса как сферического расстояния по поверхности предполагает неглубокий источник и простую структуру скоростей. Реальные годографы зависят от глубины, расстояния и среды, через которую проходит каждый луч.
Что означает положительная или отрицательная невязка?
Невязка — это разность между подобранным сферическим расстоянием и расстоянием, рассчитанным по S−P. Положительное значение означает, что подобранная точка находится дальше от этой станции, чем ее радиус по S−P; отрицательное — ближе. Большие разнонаправленные невязки показывают, что окружности не имеют одного четкого пересечения при выбранных скоростях.
Почему необходимы три станции?
Одно измерение S−P дает окружность расстояния, а не направление. Две окружности могут пересекаться в двух точках. Третья станция обычно позволяет сделать однозначный выбор, а четвертая или последующие станции делают несоответствие видимым через невязки. Станции, распределенные вокруг источника, ограничивают как широту, так и долготу лучше, чем станции, расположенные вдоль одной линии.
Что входит в погрешность определения положения?
Она проецирует единичное значение погрешности времени 1σ через коэффициент расстояния S−P и геометрию станций вблизи подобранной точки. Она не учитывает неверную модель скоростей, глубину очага, ошибки идентификации фаз, смещение часов станций или трехмерную структуру Земли, поэтому не является полной оценкой точности.