КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ

Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.

 

ОБЩИЕ СВЕДЕНИЯ


Номер 22-27-00415

НазваниеЭволюция волн цунами, вызванных землетрясением: роль динамических эффектов и фазовой дисперсии

РуководительСеменцов Кирилл Александрович, Кандидат физико-математических наук

Организация финансирования, регион Федеральное государственное бюджетное образовательное учреждение высшего образования «Московский государственный университет имени M.В.Ломоносова», г Москва

Период выполнения при поддержке РНФ 2022 г. - 2023 г. 

Конкурс№64 - Конкурс 2021 года «Проведение фундаментальных научных исследований и поисковых научных исследований малыми отдельными научными группами».

Область знания, основной код классификатора 07 - Науки о Земле, 07-504 - Механизмы возникновения и процессы распространения цунами, предвестники цунами

Ключевые словацунами, землетрясения, численное моделирование, сигма координаты, гравитационные волны, потенциальная теория

Код ГРНТИ37.25.21


СтатусЗакрыт досрочно


 

ИНФОРМАЦИЯ ИЗ ЗАЯВКИ


Аннотация
Существующие численные модели волн цунами обладают самой различной степенью сложности: от сложнейших кодов типа fully coupled, рассматривающих океан и подлежащее упругое полупространство как единую систему, до простейших кодов, решающих линейные уравнения теории длинных волн, позволяющих рассчитать трансокеанское распространение цунами даже на персональном компьютере. Выбор той или иной модели позволяет учесть одни физические факторы и пренебречь другими, а сравнительный анализ результатов моделирования позволяет оценить важность учета каждого из них. Качество описания явления напрямую зависит от точности физических представлений, положенных в основу модели. Причем важно помнить, что если речь идет об оперативном прогнозе цунами, то необходимо соблюдать разумный баланс между точностью расчетов и их скоростью. В ходе работы над проектом планируется численное воспроизведение волн цунами с помощью двух моделей, построенных в рамках потенциальной теории сжимаемой и несжимаемой жидкости. Эти модели превосходят двумерные вертикально-интегрированные коды, но уступают кодам типа fully coupled и кодам на основе уравнения Навье-Стокса. Преимущество исследуемых в проекте моделей над вертикально-интегрированными состоит в том, что, во-первых, они позволяют адекватно учесть сглаживающий эффект водного слоя и не требуют применения искусственных сглаживающих алгоритмов при вычислении начального возвышения поверхности воды в очаге цунами по заданной деформации дна, а во-вторых, корректно описывают фазовую дисперсию гравитационных волн. Результаты расчетов, полученные с помощью указанных моделей, будут сравниваться с реальными данными наблюдений, а также с расчетами, выполненными в рамках теории мелкой воды. Последнее сравнение представляет большой практический интерес, так как теория мелкой воды представляет собой наименее затратный вычислительно способ моделирования цунами, поэтому очень важно адекватно оценивать, насколько существенны те факторы, которыми мы пренебрегаем в длинноволновой модели ради достижения максимальной скорости счета. В последние годы в связи со значительным увеличением количества донных обсерваторий широко распространено использование донных датчиков давления в качестве сейсмометров. Принципиальная возможность такого использования основана на линейной связи между вертикальным ускорением дна и вариациями придонного давления, наблюдающейся в частотном диапазоне «вынужденные колебания». В рамках работы над проектом методом численного моделирования планируется выявить и конкретизировать условия, при которых эта связь нарушается. Модель, построенная в рамках потенциальной теории сжимаемой жидкости, позволит это сделать. Также с помощью обеих моделей планируется воспроизведение ряда реальных событий начала XXI века. В ходе этого воспроизведения мы планируем сравнивать между собой сами модели (по скорости счета и его точности), различные способы задания информации об очаге землетрясения (статический и динамический), различные типы очага землетрясения (Centroid Moment Tensor, Finite Fault Model), а также сопоставлять результаты моделирования с результатами расчетов в рамках теории длинных волн и записями станций DART, представляющими собой наиболее доступный, богатый и полный массив данных по цунами XXI века. Еще одно планируемое направление работ – получение зависимостей энергетических характеристик гравитационных волн, возбуждаемых поверхностными сейсмическими волнами в океане и водоемах суши, от параметров сейсмической волны и топографии дна. Интерес к механизмам генерации гравитационных волн подобного типа обусловлен тем, что их амплитуда в отличие от цунами определяется преимущественно амплитудой сейсмических колебаний дна, а не величиной остаточного смещения, то есть генерация таких волн – сугубо динамический эффект. Актуальность исследования роли динамических эффектов при генерации цунами крайне высока. Об этом могут свидетельствовать, например, катастрофические последствия события Тохоку 11 марта 2011 года, одной из причин которых явилась недооценка энергии цунами, вызванная, по-видимому, неучетом динамики горизонтальных движений дна.

Ожидаемые результаты
1. Создание новой трехмерной вертикально-разрешающей численной модели волн цунами в рамках потенциальной теории несжимаемой жидкости. В отличие от двумерных вертикально-интегрированных моделей, данная модель позволит при задании начальных условий адекватно учесть сглаживающий эффект водного слоя без применения искусственных сглаживающих алгоритмов при вычислении начального возвышения поверхности воды в очаге цунами по заданной деформации дна, а во-вторых, корректно описать фазовую дисперсию гравитационных волн. 2. Формулировка рекомендаций по выбору оптимальной численной модели для прогноза цунами в каждом конкретном случае. Рекомендации будут основаны на сопоставлении результатов численного воспроизведения цунами, выполненного с помощью моделей различного уровня сложности в сочетании с данными об очаге землетрясения различной степени детальности. Также будет учтена точность расчета, его скорость и временная задержка в появлении в открытом доступе данных об очаге землетрясения. 3. Численное воспроизведение линейной связи между вертикальным ускорением дна и вариациями придонного давления в частотном диапазоне «вынужденные колебания». На основе результатов численного воспроизведение планируется выявление и конкретизация условий, при которых эта связь нарушается, и, как следствие, уточнение границ применимости недавно предложенного дистанционного метода проверки взаимной калибровки донных датчиков давления и донных сейсмометров. 4. Оценка важности учета динамики вспарывания разрыва в очаге землетрясения для корректного воспроизведения волн цунами. Оценка будет производиться на основе численного воспроизведения ряда сильных событий XXI века и позволит оптимизировать выбор модели очага землетрясения по уровню детальности. 5. Получение зависимостей энергетических характеристик гравитационных волн, возбуждаемых поверхностными сейсмическими волнами в океане и в водоемах суши, от параметров сейсмической волны и геометрии подводного склона. Приложение полученных закономерностей для исследований в области учета роли динамики движения дна при описании генерации цунами.


 

ОТЧЁТНЫЕ МАТЕРИАЛЫ


Аннотация результатов, полученных в 2022 году
В рамках потенциальной теории несжимаемой жидкости в приближении мгновенной деформации дна построена двумерная (0xz) численная модель, позволяющая рассчитать начальное возвышение поверхности воды в очаге цунами в бассейне переменной глубины, вызванное деформацией дна произвольной формы. Для тестирования численной модели получено и применено аналитическое решение задачи о начальном возвышении водной поверхности в бассейне с плоским наклонным дном, вызванном деформацией дна треугольной формы. С помощью построенной численной модели рассчитаны начальные возвышения поверхности воды во время Курильского землетрясения 13.01.2007 и Великого восточно-японского землетрясения Тохоку 11.03.2011 (вдоль произвольно выбранных разрезов). Полученные результаты использованы для проверки приближенного метода расчета начального возвышения, известного под названием фильтр Каджиуры. Путем сравнения численного решения с результатами применения фильтра Каджиуры установлено, что в качестве «средней глубины по источнику цунами» целесообразно брать «среднюю глубину области максимальных деформаций дна». Метод дистанционной проверки качества взаимной калибровки донных манометров и Z-акселерометров, основанный на теоретически предсказанной линейной связи между вертикальным ускорением дна и вариациями придонного давления, имеющей место в частотном диапазоне «вынужденные колебания», применен к шестнадцати станциям системы глубоководных обсерваторий DONET-1. Для применения указанного метода использованы записи, полученные станциями системы DONET-1 во время трех удаленных землетрясений 2017-2018 гг. Анализ батиметрии области постановки станций DONET-1 показал, что лишь восемь станций из шестнадцати расположены на достаточном для применения метода удалении от подводных склонов. Для этих восьми станций не обнаружено заметных нарушений взаимной калибровки датчиков. По данным глубоководных станций уровня моря DART и наземных барографов IRIS исследованы проявления взрывного извержения вулкана Хунга-Тонга-Хунга-Хаапай 15.01.2022 в тихоокеанском регионе. Установлено, что атмосферная волна Лэмба, вызванная взрывом вулкана, была сформирована в 04:27 UTC ± 4 мин и распространялась над акваторией Тихого океана со скоростью 312 ± 4 м/с. Показано, что первый отчетливо различимый импульс, зарегистрированный всеми глубоководными станциями уровня моря DART, представляет собой прямое проявление атмосферной волны Лэмба. Также показано, что волна Лэмба проявляется в вариациях придонного давления с усилением, при идентичности волновых форм. Установлено, что помимо проявлений атмосферной волны Лэмба, глубоководные станции уровня моря DART записали свободные гравитационные волны, возбужденные в океане атмосферной волной Лэмба при прохождении области резкого изменения глубин, и волны цунами, распространяющиеся непосредственно из области извержения вулкана. Выполнены теоретические оценки амплитуды поверхностных гравитационных волн, возбуждаемых в океане на скачке глубин бегущим атмосферным возмущением. Показано, что волна Лэмба с амплитудой 500 Па способна возбудить в области глубоководного желоба Тонга поверхностные гравитационные волны с амплитудой 0.5 м. Выполненные оценки согласуются с данными наблюдений. Проведено численное моделирование шести реальных крупных цунами-событий начала XXI века. Каждое событие моделировалось двумя способами: сценарий «пассивной генерации цунами» и сценарий «активной генерации». Максимальные высоты волн цунами, полученные в ходе этих расчетов, сравнивались между собой. Установлено, что для всех событий (за исключением катастрофического землетрясения Тохоку 11.03.2011) разница между максимальными высотами волн цунами существенно зависит от взаимной ориентации очага землетрясения, направления распространения деформации дна и точки вывода данных. Для события Тохоку 11.03.2011 амплитуда лидирующей волны цунами, рассчитанная в рамках «активного» сценария, по всем направлениям больше, чем амплитуда, рассчитанная в рамках «пассивного» сценария, причем разница составляет до 10%. Описан прототип автоматизированной системы оценки цунамигенности землетрясения “Tsunami Observer”, проанализированы и обобщены результаты работы системы за 3.5 года ее функционирования (28.01.2018 – 1.08.2022), проведено сопоставление расчетов, выполненных системой, с данными наблюдений. Также проведена оценка времени задержки в появлении информации о параметрах очага землетрясения в трех базах данных, с которыми взаимодействует “Tsunami Observer”: USGS, GEOFON, NIED. Установлено, что NIED предоставляет данные о фокальном механизме быстрее всего – за 7-8 минут после землетрясения, а служба GEOFON предоставляет данные несколько более оперативно, нежели USGS. Также было установлено, что для землетрясений с магнитудами Mw > 7 задержка редко превышает 1 час, а данные о Finite Fault Model появляются в оперативном доступе в лучшем случае через два часа после землетрясения. 1. М. А. Носов, К. А. Семенцов, С. В. Колесов, В. В. Прядун. Вулканогенное цунами 15.01.2022 по данным глубоководных станций DART. // Доклады Российской академии наук. Науки о Земле. — 2022. — Т. 507, № 1. — С. 61–66. 2. М. А. Носов, К. А. Семенцов, С. В. Колесов, В. В. Прядун. Проявления атмосферной волны Лэмба в вариациях придонного давления. // Вестник Московского университета. Серия 3: Физика, астрономия. — 2022. — № 6. — С. 1–8. 3. К. А. Семенцов, М. А. Носов. Расчет начального возвышения поверхности воды в очаге цунами в бассейне с произвольным рельефом дна. // Математическое моделирование. — 2023. — Т.35 (принята к печати 14.11.2022) 4. S. V. Kolesov, M. A. Nosov, K. A. Sementsov, A. V. Bolshakova, G. N. Nurislamova. Automatic Tsunami Hazard Assessment System «Tsunami Observer» // Geosciences. — 2022. — 12. (принята к печати 08.12.2022)

 

Публикации

1. Колесов С.В., Носов М.А., Семенцов К.А., Большакова А.В., Нурисламова Г.Н. Automatic Tsunami Hazard Assessment System «Tsunami Observer» Geosciences, - (год публикации - 2022)

2. Носов М.А., Семенцов К.А., Колесов С.В., Прядун В.В. Вулканогенное цунами 15.01.2022 по данным глубоководных станций DART Доклады Российской академии наук. Науки о Земле, том 507, № 1, с. 61–66 (год публикации - 2022) https://doi.org/10.31857/S2686739722601004

3. Носов М.А., Семенцов К.А., Колесов С.В., Прядун В.В. Проявления атмосферной волны Лэмба в вариациях придонного давления Вестник Московского Университета. Серия 3. ФИЗИКА. АСТРОНОМИЯ, № 6. С. 1–8 (год публикации - 2022)

4. Семенцов К.А., Носов М.А. Расчет начального возвышения поверхности воды в очаге цунами в бассейне с произвольным рельефом дна Математическое моделирование, - (год публикации - 2023)