Прогнозирование землетрясений с помощью ETAS.inlabru

Основная информация

• ID статьи: 2510.13930
• Название: Earthquake Forecasting with ETAS.inlabru
• Автор: Ziwen Zhong
• Классификация: stat.AP (Статистика - Приложения)
• Дата публикации: 15 октября 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2510.13930v1

Аннотация

В данной работе исследуется прогнозирование землетрясений с использованием модели ETAS. Для решения проблемы высокой вычислительной сложности традиционных методов MCMC и сильной корреляции параметров предлагается применить байесовский вывод на основе метода INLA через пакет "inlabru". Статья проверяет эффективность метода через три эксперимента: (1) исследование влияния фиксированных параметров на апостериорные распределения других параметров, выявляя явное взаимное влияние параметров α и K; (2) сравнение распределений нормализованных интервалов между событиями для реальных и синтетических каталогов; (3) прогнозирование землетрясений в течение десяти недель после главного толчка с использованием событий перед главным толчком, оценивая точность прогноза с помощью N-теста и непрерывной оценки вероятности ранжирования (CRPS).

Исследовательский контекст и мотивация

Определение проблемы

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

Значимость

Точность прогнозирования афтершоков напрямую влияет на:

1. Предоставление пострадавшим большего времени для выживания
2. Лучшую защиту жизни спасателей
3. Оптимизацию распределения ресурсов при оказании помощи после землетрясения

Ограничения существующих методов

1. Оценка максимального правдоподобия (MLE): метод прямой и простой, но не позволяет получить неопределенность оценок параметров
2. Марковская цепь Монте-Карло (MCMC): как байесовский метод позволяет получить апостериорные распределения параметров, но имеет следующие проблемы:
   • Высокая вычислительная сложность
   • Высокая корреляция между параметрами модели
   • Сложность обеспечения хороших апостериорных распределений параметров

Исследовательская мотивация

Метод 'inlabru' на основе интегрированного вложенного приближения Лапласа (INLA) позволяет:

• Значительно сократить время подгонки модели через приближение апостериорного распределения
• Показывать отличные результаты, когда зависимости между параметрами ковариации и скрытыми переменными приблизительно линейны
• Внутренне оценивать гауссовское приближение, улучшая байесовский вывод

Основные вклады

1. Методологический вклад: применение метода inlabru на основе INLA к модели ETAS для прогнозирования землетрясений, решение проблемы вычислительной эффективности традиционных методов MCMC
2. Обнаружение взаимосвязей параметров: систематическое выявление через эксперименты взаимного влияния параметров в модели ETAS, особенно сильной корреляции параметров α и K
3. Повышение вычислительной эффективности: доказательство того, что фиксация α или K к истинным значениям может сократить время подгонки модели более чем вдвое
4. Проверка прогностической способности: установление рамки оценки прогнозов на основе реальных данных о землетрясениях, определение минимальных требований к данным для надежного прогнозирования
5. Теоретические выводы: уточнение условий, необходимых для получения хороших апостериорных распределений параметров, предоставление руководства для практического применения

Подробное описание метода

Определение задачи

Входные данные: историческая последовательность событий землетрясений, содержащая время возникновения t, местоположение s и магнитуду m Выходные данные: прогноз распределения вероятности событий землетрясений в будущем временном окне Ограничения: рамки модели ETAS на основе самовозбуждающегося точечного процесса

Архитектура модели

Основы модели ETAS

Модель ETAS основана на процессе Хокса, функция интенсивности которого:

λ_ETAS(t,m|H_t) = (μ + Σ K·e^(α(m_h-M_0))·(t-t_h/c + 1)^(-p))·π(m)

Где:

• μ: фоновая частота возникновения
• K: коэффициент производства, контролирующий интенсивность функции срабатывания
• α: масштабирование магнитуды, контролирующее влияние магнитуды исторических событий
• c: временное смещение
• p: параметр затухания афтершоков
• π(m): распределение магнитуды, подчиняющееся закону Гутенберга-Рихтера

Приближение логарифма правдоподобия

Используется разложение Тейлора и стратегия биннинга для приближения логарифма правдоподобия:

L(θ|H) = -Λ_0(T_1,T_2) - Σ Σ Λ_i(t_j^(bi), t_{j+1}^(bi)) + Σ log λ(t_i|H_{t_i})

Процесс байесовского вывода

1. Использование начальных параметров θ_0 для получения приближения логарифма правдоподобия рядом Тейлора
2. Добавление логарифма априорного распределения, вычисление моды апостериорного распределения θ_1*
3. Обновление параметров через линейный поиск: θ* = αθ_0 + (1-α)θ_1*
4. Итерация до сходимости (разница менее 1%)

Технические инновации

1. Оптимизация стратегии биннинга: применение адаптивного биннинга с узкими бинами слева для повышения точности и широкими бинами справа для снижения вычислительной нагрузки
2. Обратное преобразование вероятности: преобразование выборок гауссовского распределения в выборки целевого априорного распределения
3. Линейное приближение: реализация эффективного приближения функции правдоподобия через разложение Тейлора

Экспериментальная установка

Наборы данных

1. Землетрясение Aquila: используется для генерации синтетических каталогов и анализа параметров
2. Землетрясение Amatrice (2016-2017, Италия): реальная последовательность землетрясений с тремя главными толчками
3. Синтетические каталоги: генерируются на основе параметров землетрясения Aquila для контролируемых экспериментов

Показатели оценки

1. N-тест (тест количества): оценка точности прогнозируемого количества событий

   δ_2 = |{N_j|N_j ≤ N_obs, j=1,...,m}|/m
   

2. Непрерывная оценка вероятности ранжирования (CRPS): одновременная оценка точности и точности

   S(F,N_true) = Σ(F(k) - I(N_true ≤ k))²
   

Методы сравнения

• Полная модель ETAS в сравнении с моделями с фиксированным одним параметром
• Сравнение моделей с различным количеством исторических главных толчков
• Сравнение распределений реальных и синтетических данных

Детали реализации

Априорные распределения:

• μ ~ Gamma(0.3, 0.6)
• K, α, c ~ U(0, 10)
• p ~ U(1, 10)

Начальные значения: μ₀=0.5, K₀=0.1, α₀=1, c₀=0.1, p₀=1.1

Результаты экспериментов

Основные результаты

Эксперимент 1: Эффект фиксирования параметров

Ключевые выводы:

• Параметры α и K показывают явное взаимное влияние: фиксация одного из них к истинному значению улучшает апостериорное распределение другого
• Между c и p также существует взаимное влияние
• μ относительно независим, фиксация других параметров не влияет на апостериорное распределение μ

Повышение вычислительной эффективности:

Эксперимент 2: Распределение интервалов между событиями

• Распределение нормализованных интервалов между событиями для реальных данных и синтетических каталогов в основном совпадает
• По сравнению с экспоненциальным распределением Exp(1), наблюдается больше событий с короткими и длинными интервалами, указывая на кластеризацию
• Увеличение μ приближает распределение к Exp(1) (однородный процесс Пуассона)
• Увеличение p значительно увеличивает количество событий с короткими интервалами

Эксперимент 3: Оценка прогностической способности

Результаты прогнозирования землетрясения Amatrice:

• При использовании второго главного толчка: во все десять недель наблюдается переоценка, но истинные значения находятся в пределах 95% доверительного интервала
• При использовании третьего главного толчка: значение δ₂ близко к 0.5, прогноз более точен
• CRPS показывает, что прогноз третьего главного толчка значительно лучше, чем второго

Абляционные эксперименты

Влияние количества исторических главных толчков:

1. Без главных толчков: апостериорные распределения всех параметров, кроме μ, крайне плохие
2. Один главный толчок: можно получить базовый пригодный прогноз, но с ограниченной точностью
3. Два главных толчка: значительное улучшение прогностической способности, особенно апостериорного распределения параметра α

Анализ конкретных случаев

Порядок сложности оценки параметров:

1. μ (самый простой): может быть оценен с небольшим количеством событий
2. c, p (средний): требуют больше событий землетрясений, могут быть оценены в спокойных каталогах
3. α, K (самый сложный): требуют по крайней мере одной полной последовательности главный толчок-афтершоки; больше главных толчков значительно улучшают апостериорное распределение α

Экспериментальные выводы

1. Корреляция параметров: существует сильная корреляция между α-K и c-p, требующая совместного рассмотрения
2. Требования к данным: надежное прогнозирование требует по крайней мере одной полной последовательности главный толчок-афтершоки
3. Улучшение прогноза: увеличение количества исторических главных толчков значительно повышает прогностическую способность
4. Оптимизация вычислений: разумное использование априорной информации может значительно повысить вычислительную эффективность

Связанные работы

Основные направления исследований

1. Развитие модели ETAS: от исходной модели Ogata (1988) к различным расширенным версиям
2. Методы оценки параметров: эволюция от MLE к MCMC и далее к методам на основе INLA
3. Оценка прогнозирования землетрясений: развитие показателей оценки, таких как N-тест и CRPS

Связь данной работы с другими исследованиями

• Основана на теории процессов Хокса
• Наследует основные идеи модели ETAS
• Проводит углубленное исследование на основе применения inlabru в работах Serafini et al. (2023) и Naylor et al. (2023)

Относительные преимущества

1. Значительно выше вычислительная эффективность по сравнению с традиционными методами MCMC
2. Систематический анализ взаимосвязей между параметрами
3. Предоставление практической рамки оценки прогнозов

Выводы и обсуждение

Основные выводы

1. Эффективность метода: метод inlabru показывает отличные результаты в модели ETAS, значительно повышая вычислительную эффективность
2. Взаимосвязь параметров: взаимная зависимость параметров α и K является ключевым фактором, влияющим на производительность модели
3. Прогностическая способность: модель обладает надежной способностью прогнозирования афтершоков, качество прогноза улучшается с увеличением исторических данных
4. Требования к данным: для надежного прогнозирования требуется по крайней мере одна полная последовательность главный толчок-афтершоки

Ограничения

1. Ограничения данных: синтетические каталоги содержат только один главный толчок, что ограничивает универсальность выводов
2. Определение интервала между событиями: рассматриваются только интервалы между соседними событиями, не учитывается внутренняя структура кластеров главных толчков
3. Вычислительные вызовы: время подгонки модели значительно увеличивается с увеличением количества главных толчков, сходимость затруднена
4. Недостаточный анализ характеристик: влияние других характеристик землетрясений, помимо количества главных толчков, на апостериорное распределение еще недостаточно изучено

Будущие направления

1. Стратегии обновления параметров: разработка методов инкрементального обучения, использование апостериорного распределения предыдущего главного толчка в качестве априорного для следующего
2. Эксперименты с несколькими главными толчками: повторение экспериментов на синтетических каталогах с большим количеством главных толчков
3. Анализ кластеризации событий: исследование распределения интервалов между событиями внутри и между кластерами главных толчков
4. Исследование влияния характеристик: систематический анализ влияния различных характеристик землетрясений на производительность модели

Глубокая оценка

Преимущества

1. Методологическая инновативность: успешное применение метода INLA к модели ETAS, решение практических вычислительных проблем
2. Полнота экспериментального дизайна: три эксперимента всесторонне оценивают производительность метода с разных углов
3. Углубленный анализ параметров: систематическое выявление взаимосвязей между параметрами модели ETAS
4. Высокая практическая ценность: предоставление четких руководств по требованиям к данным для практического применения прогнозирования землетрясений
5. Высокая достоверность результатов: использование реальных данных о землетрясениях для проверки, выводы имеют практическое значение

Недостатки

1. Недостаток теоретического анализа: отсутствие теоретического объяснения корреляции параметров
2. Ограниченный диапазон экспериментов: в основном основано на данных о землетрясениях в Италии, недостаточная географическая репрезентативность
3. Неполное сравнение методов: отсутствие подробного сравнения с другими современными байесовскими методами
4. Недостаточная количественная оценка неопределенности: анализ неопределенности прогноза еще недостаточно глубок

Влияние

1. Академический вклад: предоставление новой методологии для междисциплинарной области сейсмологии и статистики
2. Практическая ценность: прямое применение к практическому прогнозированию землетрясений и оценке сейсмического риска
3. Воспроизводимость: основано на открытом пакете R inlabru, обладает хорошей воспроизводимостью
4. Потенциал распространения: метод может быть расширен на другие приложения самовозбуждающихся точечных процессов

Применимые сценарии

1. Прогнозирование афтершоков: прогнозирование активности афтершоков в краткосрочной перспективе после землетрясения
2. Оценка сейсмического риска: оценка сейсмического риска на основе исторических данных
3. Планирование аварийного реагирования: предоставление научной поддержки принятия решений для оказания помощи после землетрясения
4. Сейсмологические исследования: изучение закономерностей пространственно-временного распределения сейсмической активности

Библиография

Основные цитируемые работы включают:

• Ogata (1988): исходное предложение модели ETAS
• Hawkes (1971): теоретическая основа процессов Хокса
• Rue et al. (2017): теоретическая рамка метода INLA
• Serafini et al. (2023): применение inlabru в процессах Хокса
• Naylor et al. (2023): байесовское моделирование модели ETAS

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