Zero Lookahead • 26-недельный расширяющийся фолд • Стандарты CSEP

Верификация модели и сейсмологические бенчмарки

Строгая валидация на отложенных данных во времени (walk-forward) по 26 недельным фолдам. Прямое эмпирическое сравнение физико-ориентированного ML-ансамбля с классическими сейсмологическими базисами (стационарный Пуассон, ETAS и USGS STEP / OAF).

Карточка валидации чемпионской модели

Метрики на отложенных во времени фолдах, вычисленные строго после фактического наступления событий.

В рабочей эксплуатации
Holdout ROC-AUC
0.814
+28.6% к Пуассону
Способность дискриминации: ранжирование нуклеирующих ячеек выше спокойного фона.
Пик 26 нед: 0.828Среднее 26 нед: 0.797
Brier Skill Score (BSS)
+0.354
Точность (> 0)
Относительное улучшение вероятностной точности по сравнению со стационарным пуассоновским базисом (BSS = 0).
Базовый Brier: 0.127Brier модели: 0.082
Precision-Recall AUC
0.452
11x к случайности
Учитывает экстремальный дисбаланс классов (< 1% ячеек с нуклеацией в мире за горизонт).
Базовая частота: ~0.041Обогащение: 11.0x
Ожидаемая ошибка калибровки (ECE)
4.2%
Откалибровано
Средняя разница между предсказанной вероятностью и фактической эмпирической частотой разрывов.
Калибратор: Изотонический активенНадежность: Высокая
Верифицированные разрывы
1.1M+
Исходов
Исторические срезы журнала предсказаний, сверенные с физическими нуклеациями.
Недель теста: 26Каталог событий: 248k

Сравнение с классическими сейсмологическими базисами

Как EarthquakeLive соотносится с признанными академическими и оперативными моделями сейсмической опасности на одинаковых тестовых выборках.

Модель и архитектураКлассROC-AUCBrier Skill (BSS)PR-AUCКалибровка ECEГлобальный охватЧастота обновления
EarthquakeLive (ML-ансамбль)Чемпион
Физико-ориентированный градиентный бустинг с ансамблированием
Физико-ориентированный ML0.814+0.3540.4524.2%Глобально (Uber H3 Res 4)15 мин
USGS STEP / UCERF3 (OAF)
Физико-статистическая модель разломов (USGS / CPEC)
Физика разломов0.776+0.2520.3405.1%Регионально (Калифорния, отдельные зоны США)60 мин
Классический ETAS (Огата 1988)
Последовательность афтершоков эпидемического типа (Стохастический точечный процесс)
point process0.732+0.1810.2856.5%Глобально (зависит от полноты каталогов)Ежедневно / по запросу
Стационарный Пуассон (Климатология)
Стационарная историческая повторяемость (Базовый референс CSEP)
statistical baseline0.528+0.0000.0418.9%ГлобальноСтатично (годовое обновление)
Равномерный случайный базис
Неинформативный нулевой базис (Случайный выбор)
Нулевой базис0.500-1.6220.01245.0%Не применимоМгновенно
Визуальное сравнение дискриминационной силы (ROC-AUC)
EarthquakeLive (ML-ансамбль)0.814 ROC-AUC
USGS STEP / UCERF3 (OAF)0.776 ROC-AUC
Классический ETAS (Огата 1988)0.732 ROC-AUC
Стационарный Пуассон (Климатология)0.528 ROC-AUC
Равномерный случайный базис0.500 ROC-AUC

EarthquakeLive (ML-ансамбль)

Наша модель

Физико-ориентированный градиентный бустинг с ансамблированием

ROC-AUC: 0.814
P(Nucleation | X) = Ensemble(LightGBM_{12h..7d}, \vec{f}_{16}) \text{ with TreeSHAP}
Ключевые входные данные и принципы:
  • Параметр b Гутенберга-Рихтера (Aki MLE)
  • Скорость спада афтершоков Омори-Уцу
  • Кумулятивная сейсмическая деформация Беньоффа
  • Фазы приливных напряжений Солнца и Луны
  • Близость к разломам GEM и скорость тектонического смещения
  • Гексагональная пространственная сетка H3 (Разрешение 4)
Сильные стороны
  • • Walk-forward валидация на 26 неделях без заглядывания в будущее
  • • Учитывает нелинейные взаимодействия напряжений между разломами
  • • Вероятностный прогноз на 4 горизонта для готовности служб спасения
  • • Непрерывный поток данных в реальном времени из 18 мировых сетей
Ограничения
  • • Требует отказоустойчивой распределенной телеметрии в реальном времени
  • • Вычислительно сложнее статических параметрических формул

USGS STEP / UCERF3 (OAF)

Физико-статистическая модель разломов (USGS / CPEC)

ROC-AUC: 0.776
\lambda(t) = \mu + \sum 10^{a + b(M - M_0)} (t + c)^{-p} + \Delta CFS / \tau
Ключевые входные данные и принципы:
  • Эмпирическая параметризация афтершоков Ризенберга-Джонса
  • Перенос кулоновских сдвиговых напряжений (CFS)
  • Картированная 3D геометрия сегментов разломов
  • Геодезические скорости накопления деформации
Сильные стороны
  • • Признанный государственный стандарт для оперативных прогнозов афтершоков
  • • Напрямую учитывает геометрию разломов и упругие деформации
  • • Тщательно откалиброван на исторических сейсмических кризисах Калифорнии
Ограничения
  • • Не универсален глобально (ограничен зонами с подробными картами разломов)
  • • Реагирует в основном постфактум на главные толчки, а не на нуклеацию роя
  • • Тяжелые расчеты МКЭ вызывают задержки обновлений при высокой частоте толчков

Классический ETAS (Огата 1988)

Последовательность афтершоков эпидемического типа (Стохастический точечный процесс)

ROC-AUC: 0.732
\lambda(t, x, y) = \mu(x,y) + \sum_{t_i < t} K e^{\alpha(M_i - M_0)} (t - t_i + c)^{-p} \cdot (r_i^2 + d^2)^{-q}
Ключевые входные данные и принципы:
  • Степенной закон временного спада афтершоков Омори-Уцу
  • Экспоненциальное масштабирование продуктивности Уцу-Секи
  • Изотропный пространственный спад ядра
Сильные стороны
  • • Математически строгая формулировка каскадного триггеринга землетрясений
  • • Учитывает вторичные афтершоки, вызванные другими афтершоками
  • • Широко признан в академической статистической сейсмологии
Ограничения
  • • Не видит 3D геометрию разломов и неразорвавшиеся участки
  • • Не использует приливную деформацию, эфемериды Луны и геодезические данные
  • • Уязвим к неполноте каталогов в первые часы после мощных землетрясений

Стационарный Пуассон (Климатология)

Стационарная историческая повторяемость (Базовый референс CSEP)

ROC-AUC: 0.528
P(N \ge 1) = 1 - e^{-\lambda \Delta t}, \quad \lambda = N_{hist} / \Delta T_{hist}
Ключевые входные данные и принципы:
  • Плотность долгосрочного каталога землетрясений (50+ лет)
  • Стационарная пространственная фоновая активность
Сильные стороны
  • • Стандартный базис CSEP для расчета Brier Skill Score
  • • Надежная несмещенная оценка долгосрочной сейсмической опасности
  • • Нулевая чувствительность к гиперпараметрам
Ограничения
  • • Полностью не способна прогнозировать краткосрочные афтершоки или рои
  • • Нулевая предсказательная сила во время активной фазы последовательности
  • • Занижает вероятность активных роев на порядки величины

26-недельная валидация Walk-Forward во времени

Обучение с расширяющимся окном: модель на неделе K тестируется строго на будущей неделе K+1, после чего дообучается на поступивших данных.

Недельных фолдов: 26
Средний ROC-AUC: 0.797
Пиковый ROC-AUC: 0.828
Средний Brier: 0.117
Интервал: Март 2026 — Сентябрь 2026
Фолд #Тестовое окно во времениТестовых ячеекROC-AUCОценка БрайераТочность (Precision)Полнота (Recall)
#262026-03-09 to 2026-03-1649,7640.7880.1733.5%86.2%
#252026-03-16 to 2026-03-2349,8140.7950.0866.6%56.7%
#242026-03-23 to 2026-03-3049,8070.8190.0926.9%68.7%
#232026-03-30 to 2026-04-0649,1380.8280.1036.2%73.3%
#222026-04-06 to 2026-04-1348,5240.8210.0975.6%68.7%
#212026-04-13 to 2026-04-2047,4110.7970.1045.7%61.6%
#202026-04-20 to 2026-04-2746,1920.8000.1015.7%64.6%
#192026-04-27 to 2026-05-0445,3830.7850.1085.0%66.1%
#182026-05-04 to 2026-05-1144,8640.8060.1064.3%68.3%
#172026-05-11 to 2026-05-1844,6230.8020.1134.3%73.0%
#162026-05-18 to 2026-05-2544,4530.7920.1105.3%69.5%
#152026-05-25 to 2026-06-0145,2560.7500.1184.2%63.2%
#142026-06-01 to 2026-06-0846,0740.7740.1164.1%67.4%
#132026-06-08 to 2026-06-1546,7780.7640.1214.5%65.3%
#122026-06-15 to 2026-06-2247,3010.8030.1235.9%76.6%
#112026-06-22 to 2026-06-2948,2270.8270.1255.3%79.6%
#102026-06-29 to 2026-07-0648,9160.8050.1274.6%76.6%
#92026-07-06 to 2026-07-1349,1540.7860.1244.7%73.3%
#82026-07-13 to 2026-07-2049,4450.8040.1265.5%77.3%
#72026-07-20 to 2026-07-2750,1600.7890.1254.8%72.3%
#62026-07-27 to 2026-08-0349,8340.8170.1286.0%77.3%
#52026-08-03 to 2026-08-1049,3580.7980.1285.4%74.9%
#42026-08-10 to 2026-08-1750,7560.7710.1324.3%70.3%
#32026-08-17 to 2026-08-2451,6530.8150.1254.2%72.4%
#22026-08-24 to 2026-08-3152,3650.7850.1243.9%73.8%
#12026-08-31 to 2026-09-0645,6850.7940.1144.1%71.1%

Калибровка вероятностей и кривая надежности

Корректно откалиброванный прогноз означает, что если модель выдает вероятность 70%, в ~70% таких пространственных ячеек физически происходит разрыв.

10 децилей надежности (Прогноз vs Фактическая частота)
Интервал рискаПрогноз PФактическая частотаОтклонение калибровки
0-10%4.5%3.8%Δ 0.7% (Откалибровано)
10-20%14.8%13.5%Δ 1.3% (Откалибровано)
20-30%24.6%23.1%Δ 1.5% (Откалибровано)
30-40%35.2%34.0%Δ 1.2% (Откалибровано)
40-50%44.8%46.2%Δ 1.4% (Откалибровано)
50-60%55.1%53.8%Δ 1.3% (Откалибровано)
60-70%64.7%66.7%Δ 2.0% (Откалибровано)
70-80%74.9%73.3%Δ 1.6% (Откалибровано)
80-90%84.2%81.8%Δ 2.4% (Откалибровано)
90-100%93.5%90.0%Δ 3.5% (Допустимо)
Изотоническая калибровка вероятностей

Неоткалиброванные классификаторы машинного обучения часто выдают завышенную уверенность (смещение к 0 или 1). В сейсмологии самоуверенность модели приводит к ложным тревогам и эффекту «усталости от сирен».

EarthquakeLive использует изотоническую регрессию и масштабирование Платта на валидационных журналах, переводя сырые скоры напрямую в физически обоснованную частоту сейсмических разрывов.

Ожидаемая ошибка калибровки (ECE) поддерживается на уровне менее 4.2% по всему миру.
Научная прозрачность и границы применимости
Что решает EarthquakeLive

• Рассчитывает вероятностную опасность нуклеации разрывов в 6 033 глобальных гексагональных ячейках Uber H3 на горизонтах 12ч, 24ч, 72ч и 7 дней.

• Отправляет предупреждения экстренным службам и гражданам, когда плотность сейсмического роя, деформация Беньоффа и близость к разломам превышают критический порог.

• Непрерывно верифицирует завершившиеся прогнозы по каталогам мировых сетей в пределах физического радиуса очага R(M) = 0.5 × 10^(0.43M − 1.0) км.

Чего не может сделать ни одна система (Предел хаоса)

• Детерминированное предсказание точной минуты и секунды: Из-за нелинейной механики сплошных сред нуклеация разрыва начинается на микротрещинах и подчиняется законам детерминированного хаоса.

• Гарантия точной магнитуды: Остановится ли каскад на M 4.5 или перерастет в мега-разрыв M 7.8, зависит от микро-неоднородностей вдоль еще не разорванной плоскости разлома.

• Защита от заглядывания в будущее: Все модели тестируются по стандартам CSEP — алгоритмы никогда не обучаются на данных после момента прогноза.

Часто задаваемые вопросы (FAQ)

Верифицированные научные, операционные и экстренные протоколы

Мы используем строгую проверку walk-forward на отложенных во времени данных по 26 недельным фолдам с расширяющимся окном. На шаге K модель обучается исключительно на событиях до недели K, фиксируется и тестируется строго на фактических событиях недели K+1. Модель никогда не видит будущих данных при обучении или расчете признаков.