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

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

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

Номер 20-17-00180

НазваниеРазвитие сценарного подхода в задачах оценки сейсмической опасности и риска

РуководительШебалин Петр Николаевич, Доктор физико-математических наук

Организация финансирования, регион Федеральное государственное бюджетное учреждение науки Институт теории прогноза землетрясений и математической геофизики Российской академии наук (ИТПЗ РАН), г Москва

Период выполнения при поддержке РНФ

2020 г. - 2022 г.

, продлен на 2023 - 2024. Карточка проекта продления (ссылка)

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

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

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

Код ГРНТИ37.01.17

СтатусУспешно завершен

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

Аннотация
Актуальность задач точной оценки сейсмической опасности и сейсмического риска в современных условиях определяется двумя факторами. С одной стороны, опыт начала XXI столетия показывает, что наиболее разрушительные землетрясения часто происходят в областях существенной недооценки сейсмической опасности. Так, например, с момента публикации финальной карты Программы Глобальной Оценки Сейсмической Опасности (GSHAP) в 1999 году эффект каждого из более 60 землетрясений с магнитудой от 7.5 и выше был недооценен, причем в половине случаев вместо ожидаемых «легких» повреждений плохо построенных домов имели место «существенные», «значительные», или даже «тотальные» разрушения индустриально построенных сооружений, что суммарно определило невосполнимые потери более 700 тысяч человек погибших. Это приводит к тому, что на обновленных картах сейсмического районирования после подобных землетрясений повышается ожидаемая балльность в некоторой области, включающей его эпицентр. Так, например, 14 февраля 2013 г. на территории Республики Саха (Якутия) Российской Федерации произошло сильное землетрясение, которое явилось неожиданным для карт GSHAP, ОСР-97 и ОСР-2012. На перечисленных картах недооценка интенсивности сотрясений в эпицентральной зоне заведомо превышает 2 балла. В частности, интенсивность сотрясения в эпицентре была 8-9 баллов, а на карте ОСР-2012, характеризующей повторяемость превышения балльности 1 раз в 100 лет, эпицентр основного толчка 14 февраля 2013 г. находится в зоне 5.5-6 баллов (и даже на карте, характеризующей повторяемость превышения 1 раз в 2500 лет ожидаемое значение 7.5-8 баллов). Недооценка опасности вызвана либо занижением сейсмического потенциала в этой области (например, Олюторское землетрясение 2006, Катав-Ивановское землетрясение на Южном Урале 2018), либо недостаточным учетом локальных условий (например, землетрясение в Мичиоакан 1985, Мексика, или в Кобе 1995, Япония), либо обоими факторами (как это было в районе Спитакского землетрясения в Армении в 1989 г.). С другой стороны, в большинстве зон, отнесенных на картах ОСР к опасным, сильные воздействия не наблюдались со времени составления карт и, возможно, не будут наблюдены в течение многих десятков лет. Конечно, если речь идет о безопасности людей, об экологической безопасности, этот аргумент не должен быть определяющим при строительстве новых объектов или при укреплении старых. Это, очевидно, относится к 9-балльным зонам, в которых землетрясения могут приводить к обрушению зданий, или при строительстве атомных электростанций, плотин, хранилищ опасных отходов и т.д. Но в зонах умеренной опасности, в частности в 7- и 8-балльных зонах возникает вопрос о целесообразности огромных затрат на сейсмостойкость зданий и сооружений. Особенно актуально это становится в современных условиях, когда плановая продолжительность эксплуатации новых зданий и сооружений сокращается. Появление дешевых строительных материалов, материалов, допускающих повторное использование, совершенствование строительных технологий будет приводит к дальнейшему сокращению сроков эксплуатации зданий. В этих условиях локальное завышение ожидаемой балльности (интенсивности) фактически будет приводить к огромным экономическим потерям. Предварительные подсчеты показывают, что такие потери в среднем значительно превышают потери, связанные с занижением сейсмической опасности. Так, например, в Иркутской области основанные на оценках ОСР-97 затраты на сейсмостойкое строительство оказались фактически завышенными на два порядка. В настоящее время при составлении карт ОСР снижение балльности отдельных зон часто происходит не за счет более точных расчетов сейсмологов и геологов, а по чисто экономическим причинам, под давлением заинтересованных ведомств. Для сокращения экономических потерь от недооценки и переоценки сейсмической опасности необходимо развивать исследования по следующим основным направлениям: 1) совершенствование методов оценки сейсмического потенциала зон тектонических разломов (в современных условиях дополнительное внимание должно уделяться возможному изменению сейсмического потенциала в результате деятельности человека); 2) совершенствование методов определения возможного сейсмического воздействия землетрясений на здания и сооружения в зависимости от структуры очага землетрясения (с привлечением информации о региональной системе разломов) и его спектрально-энергетических характеристик, от характера затухания сейсмических волн с расстоянием от источника, от локальных грунтовых условий; 3) постепенный переход от оценок сейсмической опасности, определяющих ожидаемую максимальную балльность или пиковое ускорение грунта в заданном интервале времени, к расчетам сейсмического риска, определяющим ожидаемые экономические потери от воздействия землетрясений. Такие расчеты позволят в конкретных случаях оптимизировать расходы на повышение сейсмостойкости зданий при строительстве, а также задействовать систему страхования от воздействия землетрясений. Авторы заявленного проекта считают, что наиболее подходящим инструментом для расчетов по всем трем указанным направлениям является сценарный подход. В частности, под сценарным подходом авторы понимают многократное моделирование конкретных вероятных землетрясений, совокупность которых характеризует реальность. В частности, сейсмический потенциал какой-либо области при сценарном подходе представляется в виде синтетического каталога землетрясений, свойства которого совпадают со свойствами реального, а размер значительно больше и достаточен для статистически обоснованных оценок. Сценарный подход для расчета сейсмического воздействия состоит в построении синтетических сейсмограмм на основе детального моделирования очага потенциального землетрясения, с учетом затухания волн на пути до заданной точки. Для расчета сейсмического риска сценарный подход подразумевает многократное моделирование конкретных сценариев землетрясений (момент, очаг, воздействие), оценка экономических потерь для конкретных объектов или групп объектов в каждом сценарии и построение на этой основе статистических оценок рисков. Использование сценарного подхода предполагает сложные расчеты на ЭВМ, возможные лишь с использованием суперкомпьютеров. В том числе поэтому в России сценарный подход для оценки сейсмической опасности пока не применялся. В Европе, США и Японии (в последнее время и в Китае) сценарный подход реализуется и становится все более популярным, в частности, благодаря поддержке крупных перестраховочных компаний, таких как Swiss Re, Munich Re, SCOR. Вместе с тем рост популярности актуализировал противоречие между общепринятыми моделями сейсмичности и закономерностями, установленными по реальным данным. Разработка такой модели позволит устранить обозначенное противоречие и составит новизну проекта. Другая составляющая новизны заключается в развитии сценарного подхода для оценки сейсмических воздействий с учетом современных представлений об очаге сильного землетрясения.

Ожидаемые результаты
Проект будет выполняться параллельно по трем составляющим в их взаимодействии: 1) развитие сценарного подхода для оценки сейсмического потенциала; 2) развитие сценарного подхода для оценки сейсмических воздействий от возможных землетрясений; 3) развитие сценарного подхода для вероятностных оценок экономического ущерба от землетрясений. В результате выполнения проекта будет разработана методика расчета сейсмической опасности в заданной точке с использованием сценарного подхода. По первому направлению будет разработана модель сейсмичности, в основе которой будет использована хорошо известная модель ETAS (Epidemic Type Aftershock Sequence) в ее изотропном пространственно-временном варианте [Zhuang et al., 2004]. В модели будет учтен недавно установленный грантозаявителями экспоненциальный закон продуктивности. В модели также будет учтена фрактальная природа пространственной неоднородности сейсмического режима с использованием локальных оценок параметров Общего закона подобия землетрясений, обобщающего классическое соотношение Гутенберга-Рихтера. Учет в модели новых закономерностей, подтвержденных по данным многолетних сейсмологических наблюдений на глобальном и региональном уровнях, позволяет обоснованно ожидать, что данный результат превысит мировой уровень. По второму направлению будет выполнено моделирование вероятного сейсмического воздействия, полученного с помощью функции источника сложной структуры в пространстве и во времени (что согласуется с современными представлениями об очаге сильного землетрясения), с применением современных вычислительных мощностей, для нескольких сценарных сильных землетрясений с учетом имеющихся сейсмологических, геологических, макросейсмических и геодинамических данных о системе региональных разломов, для районов прилегающих к Байкало-Амурской Магистрали. Моделирование сейсмического воздействия на основе комплексирования данных различных типов с целью построения различных сценариев соответствует последним тенденциям мировой практики оценивания сейсмической опасности. Это обосновывает соответствие ожидаемого результата мировому уровню. По третьему направлению для оценки сейсмических рисков для населения, городов и объектов инфраструктуры Иркутской области и Прибайкалья, включая территорию районов, прилегающих к Байкало-Амурской Магистрали и нефтепровода «Сила Сибири», будут использованы результаты расчетов, полученные в рамках первого и второго направлений. Ожидаемые оценки сейсмических рисков в районах проживания людей и размещения ответственных народнохозяйственных объектов определяют значимость результатов для экономики и социальной сферы. Подчеркнем, что универсальность разрабатываемого подхода позволит использовать ожидаемые результаты проекта для любых других сейсмоопасных территорий, что является важной составляющей решения общественно значимой задачи снижения рисков и ликвидации чрезвычайных ситуаций природного характера.

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

Аннотация результатов, полученных в 2020 году
Для целей корректного моделирования сейсмического процесса проведено исследование его иерархической структуры. В литературе сложилось представление о сейсмическом процессе, состоящем из основных толчков, форшоков и афтершоков. Однако такое представление усложнено наличием афтершоков у тех событий, которые сами являются афтершоками, а также наложением афтершоковых последовательностей от разных землетрясений. Таким образом, возможно проявление афтершоков афтершоков, афтершоков афтершоков афтершоков и тд. Такая иерархическая структура может характеризоваться уровнями иерархии, если считать, что каждое сейсмическое событие является афтершоком только одного предшествующего ему события. Для этого мы используем метод «ближайшего соседа» Заляпина-Бен-Зиона [Zaliapin, Ben-Zion, 2013]. Для каждого сейсмического события минимизируется функция «соседства» [Baiesi, Paczuski, M., 2004] по отношению ко всем предшествующим событиям и таким образом находится единственное событие-триггер. Функция «соседства» определяет степень взаимозависимости двух событий. Для рассматриваемого региона определяется пороговое значение этой функции, выше которого события можно считать независимыми. Если значение функции для «ближайшего соседа» оказывается выше порогового значения, то считается, что рассматриваемое событие не имеет триггера и тем самым является фоновым событием. В процессе подготовки к реализации данного проекта был установлен экспоненциальный закон продуктивности для землетрясений. В рамках проекта мы исследовали иерархическую структуру сейсмичности с учетом этого закона по данным мирового и 7 региональных каталогов землетрясений. Оказалось, что экспоненциальный характер распределения продуктивности сохраняется для разных уровней иерархии, при этом значение параметра экспоненциального распределения меняется незначительно. Была модифицирована функция «соседства» Байеси-Пачуцки, использующая точечную модель взаимодействия землетрясений – для пространственной компоненты функции введен масштабирующий коэффициент, зависящий от магнитуды главного события (эта зависимость установлена эмпирически). Аналогично работе [Zaliapin, Ben-Zion, 2013] с помощью модифицированной функции «соседства» было проведено сравнение реальной и модельной сейсмичности Калифорнии. Выявленные различия связанной (афтершоки) и фоновой сейсмичности оказались еще более контрастными, чем в цитированной работе. Полученные результаты представляются чрезвычайно перспективными для реализации целей проекта, поэтому в 2021 г. это направление исследований будет активно продолжено. Параметры сейсмического режима могут существенно варьировать вод воздействием различных факторов. Для целей проекта важно установить степень и характер влияния на сейсмичность режима флюидов в разломных структурах, а также влияние периодических воздействий. Для океанических землетрясений в зонах субдукции Камчатки и Новой Зеландии было установлено, что повторяемость афтершоков значительно возрастает (до двух раз) на максимумах как приливов, так и отливов. При этом для взбросовых землетрясений эффект сильнее на приливах, а для нормальных землетрясений на отливах. Результаты опубликованы в журнале Frontiers in Earth sciences [Shebalin, Baranov, 2020]. Установлена связь сейсмической активности в районе озера Байкал с сейшами на озере. Кроме того, с помощью непараметрической статистики критерия Кюйпера подтверждена связь повторяемости сейсмических событий с периодами лунных и солнечных циклов. Были исследованы временные вариации контрольного параметра Общего закона подобия землетрясений до и после сильных землетрясений. Показано, что коэффициенты ОЗПЗ в исследуемой области зависят от времени и коррелированы, а динамические изменения параметров в трех подобластях подчеркивают различие сейсмических режимов. Эти результаты исследования опубликованы в журнале Frontiers in Earth sciences [Nekrasova, Peresan, 2020]. Проведена модификация алгоритма прогноза области возможного возникновения сильных афтершоков для мониторинговой системы AFCAST, созданной в ИТПЗ РАН в рамках предшествующего проекта РНФ. Модифицированный алгоритм учитывает степенной характер убывания числа афтершоков с расстоянием от эпицентра главного толчка. Осуществлен перенос системы AFCAST на новый сайт ИТПЗ РАН в сети Интернет: https://www.itpz-ran.ru/ru/resultaty/maps-and-databases/afcast/. Разработана схема адаптированных расчетов, основанная на модели сейсмического процесса и построении на этой базе синтетического каталога землетрясений, свойства которого близко совпадают со свойствами реальной сейсмичности, а охватываемый период времени выбирается произвольно и может достигать геологических времен. Модель строится на основе оценки параметров Общего закона подобия и закона продуктивности для землетрясений. Синтетические каталоги позволяют рассчитать ожидаемые эффекты сейсмического воздействия на конкретный объект или группу объектов по известным законам затухания. Наиболее точные оценки дополнительно будут включать моделирование сейсмического излучения от предполагаемых протяженных очагов. Для территории Байкальской рифтовой системы моделировались очаги землетрясений с Mw=7-8 с учетом возможного отличия высокочастотной части спектров от очаговой модели Аки-Бруна («омега-квадрат»). В рамках проекта используется две модели протяженного очага. В первой модели стохастического моделирования [Mai, Beroza, 2002] распределение смещений в очаге рассматривается как случайный процесс с экспоненциальной автокорреляционной функцией или как гауссовский процесс, начинающийся в заданной точке вспарывания. Вторая модель подразумевает расчет для сложной очаговой структуры системы сейсмогенерирующих разломов [Käser M. et al., 2010]. В первый год выполения проекта моделирование проводилось по первой модели. В рамках проекта принято в печать 2 статьи (обе в журнале первого квартиля по SJR; на заключительной стадии подготовки, будут опубликованы в декабре 2020 г.), опубликованы тезисы 5 докладов на международных конференциях, один из докладов сделан в режиме онлайн (Генеральная ассамблея Европейского союза наук о Земле), 4 других доклада (включая один приглашенный пленарный) перенесены на март 2021 г. Информация о проекте опубликована в электронном журнале Вестник Отделения наук о Земле РАН и доступна по адресу https://onznews.wdcb.ru/dec20/modeli-sejsmicheskogo-protsessa-dolzhny-spasti-zhizni.html Литература 1. Baiesi M., Paczuski M. Scale–free networks of earthquakes and aftershocks // Phys. Rev. E. 2004. V. 69. No 6. doi 10.1103/PhysRevE.69.066106. 2. Mai P.M., Beroza G.C. Source scaling properties from finite-fault-rupture models // Bulletin of the Seismological Society of America, 2000, v. 90, №3, P. 604-615. 3. Käser M., Castro C., Hermann V., Pelties C. SeisSol–a software for seismic wave propagation simulations // in High Performance Computing in Science and Engineering, Garching/Munich 2009 (pp. 281-292). Springer, Berlin, Heidelberg. 4. Zaliapin I., Ben-Zion Y. Earthquake clusters in Southern California I: Identification and stability // Journal of Geophysical Research: Solid Earth. 2013. V. 118. No 6. P. 2847–2864. doi 10.1002/jgrb.50178

Публикации

1. Некрасова А.К., Пересан А. Uniﬁed Scaling Law for Earthquakes: space-time dependent assessment in Friuli-Venezia Giulia region Frontiers in Earth Sciences, - (год публикации - 2020) https://doi.org/10.3389/feart.2020.590724

2. Шебалин П.Н., Баранов А.А. Aftershock rate changes at different ocean tide heights Frontiers in Earth Scinces, - (год публикации - 2020) https://doi.org/10.3389/feart.2020.559624

3. Воробьева И.А., Шебалин П.Н., Баранов С.В. Deficit of large aftershocks and postseismic deformation Problems of Geocosmos–2020 Proceedings of the XIII International Conference and School, S33 (год публикации - 2021)

4. Кособоков В.Г., Щепалина П.Д., Некрасова А.К. Earthquake preferred days in the Lake Baikal and Yunnan-Sichuan Region EGU General Assembly 2020, Online, 4–8 May 2020, EGU General Assembly 2020, Online, 4–8 May 2020, EGU2020-720, https://doi.org/10.5194/egusphere-egu2020-720 (год публикации - 2020) https://doi.org/10.5194/egusphere-egu2020-720

5. Некрасова А.К., Кособоков В.Г. Unified Scaling Law for Earthquakes in assessment of seismic hazard and associated risks Problems of Geocosmos–2020 Proceedings of the XIII International Conference and School, S17 (год публикации - 2020)

6. Шебалин П.Н. New quantitative characteristics of seismicity and their relationship with regional geodynamics Problems of Geocosmos–2020 Proceedings of the XIII International Conference and School, U1 (год публикации - 2021)

7. Шебалин П.Н., Баранов С.В. Earthquake productivity law Problems of Geocosmos–2020 Proceedings of the XIII International Conference and School, S32 (год публикации - 2021)

8. - Модели сейсмического процесса должны спасти жизни Вестник Отделения наук о Земле Российской академии наук, 04.12.2020 (год публикации - )

Аннотация результатов, полученных в 2021 году
Во второй год выполнения проекта в рамках задачи построения реалистичной модели сейсмического режима, исследования были направлены, главным образом, на построение модели сложного афтершокового процесса с учетом экспоненциального закона продуктивности землетрясений и конечных размеров очагов землетрясений, а также на разработку методики оценки параметров этой модели по реальным данным. Под продуктивностью землетрясения понимается число сейсмических событий в определенном диапазоне магнитуды, имеющих причинно-следственную связь с эти землетрясением. Недавно установленный экспоненциальный закон продуктивности землетрясений [Shebalin et al., 2020] является одним из ключевых свойств пространственно-временного группирования землетрясений, которое при моделировании сейсмического процесса применялось только его авторами. В рамках проекта был подтвержден закон продуктивности для широкого диапазона магнитуд (4 единицы магнитуды), а также проверен ряд гипотез с целью поиска интерпретации экспоненциального характера распределения продуктивности. Было показано, что продуктивность не является следствием закона Гутенберга-Рихтера. Разнообразие локальных условий в пространстве, например уровня фоновой сейсмичности, также не может объяснить экспоненциальный закон продуктивности. В рамках запланированных на 2021 г. исследований не удалось найти возможные варианты интерпретации экспоненциального закона продуктивности. Поиск будет продолжен в 2022 г. Для того, чтобы добиться близкого совпадения свойств моделируемого и реального сложного афтершокового процесса, необходимо, чтобы способ выделения связанных событий соответствовал той модели, которая будет использоваться для синтезирования каталога таких событий. Метод ближайшего соседа [Zaliapin, Ben-Zion, 2013] с использованием функции близости Баеси [Bayesi , Pascuzky, 2004], используемые для декластеризации - отделения связанных событий от независимых (фоновых) событий в каталоге землетрясений - обеспечивают это не полной мере. В рамках выполнения проекта в 2021 г. была исследована модификация функции Баеси, «симметричная» предполагаемому способу построения синтетического каталога связанных событий. Временная и пространственная компоненты были заменены на функции, обратные закону Омори-Утсу, и степенному закону убывания интенсивности афтершоков с расстоянием от гипоцентра основного толчка с фиксированной добавкой, моделирующей конечный размер очага землетрясения. Модифицированная функция не имеет нулевых значений, что исключает неоправданное отнесение к связанным событиям событий, совпадающих по времени или по положению гипоцентра, но значительно различающихся соответственно по пространству или по времени. Декластеризация каталогов землетрясений Японии и Прибайкалья с использованием двух вариантов функции близости дало близкие результаты; при использовании модифицированного варианта функции близости совпадающие по времени или гипоцентру события, далекие друг от друга по другому параметру не были идентифицированы как связанные. Была уточнена модель пространственного распределения афтершоков (степенное спадание с расстоянием от эпицентра основного толчка). На примере подробного каталога островной части Японии было установлено подобие этого распределения для разных магнитуд M основного толчка: при введении нормировки для расстояния от эпицентра основного толчка вида 10^(B*M), где B-коэффициент со значением около 0.5, распределения близко совпадают. Также была предложена и исследована модифицированная модель, учитывающая направление очага основного толчка [Баранов, Шебалин, 2017]. Оба варианта модели легко реализуются при моделировании сейсмического потока. На основе полученных в 2021 г. результатов был разработан алгоритм синтезирования последовательностей афтершоков в пространстве-времени-магнитуде. Разработана методика оценки параметров сложной модели афтершокового процесса по реальным данным. Разработанные в первый год выполнения проекта адаптированные методы расчета сейсмических рисков были применены для отдельных объектов Прибайкалья и Алтая-Саян. Для этих целей актуализированы и уточнены результаты анализа и картирования коэффициентов Общего закона подобия для землетрясений. Проведена оценка сейсмического риска протяженных линейных объектов Прибайкалья, а именно участков Байкало-Амурской и Транссибирской магистрали и газопровода Ковыкта – Саянск – Иркутск. Были продолжены работы по моделированию волновых форм для конкретных точек на поверхности (в том числе, для г. Иркутск). Для моделирования было выбрано Култукское землетрясение (27 августа 2008 г., MW=6.3, 76 км от ст. «Иркутск»). Пиковые амплитуды сценарных сейсмограмм были сопоставлены с реальными наблюденными данными о сильных движениях [Мельникова и др., 2014]. Для построения синтетических («сценарных») акселерограмм (от модельных очагов с М = 6, 6.5 и 7) впервые для Байкальского региона учитывалась структурная сложность источника за счет вариаций смещений (slip) в каждом из сегментов грида. Результаты моделирования показывают сопоставимые баллы (I = 6.2, 6.4, 6.7) непосредственно для территории города Иркутск, несмотря на разницу в магнитуде. В 2021 г. опубликовано или принято в печать 3 статьи, и представлено 5 докладов на виртуальных международных конференциях (ESC, EGU, IASPEI). Литература 1. Баранов С.В., Шебалин П.Н. О прогнозировании афтершоковой активности. 2. Оценка области распространения сильных афтершоков // Физика Земли. 2017. № 3. С. 43-61. DOI:10.7868/S0002333717020028 2. Shebalin P.N., Narteau C., Baranov S.V. Earthquake productivity law // Geophysical Journal International. 2020. V. 222. № 2. P. 1264-1269. DOI:10.1093/gji/ggaa252 3. Baiesi M., Paczuski M. Scale–free networks of earthquakes and aftershocks // Phys. Rev. E. 2004. V. 69. No 6. doi 10.1103/PhysRevE.69.066106. 4. Zaliapin, I. & Ben-Zion, Y., 2013. Earthquake clusters in southern California I: identification and stability, J. geophys. Res., 118(6), 2847–2864. 5. Мельникова В.И., Гилёва Н.А., Радзиминович Я.Б., Середкина А.И. Култукское землетрясение 27 августа 2008 г. с MW=6.3, I0=8-9 (Южный Байкал) // Землетрясения Северной Евразии, 2008 год. – Обнинск: ГС РАН, 2014. С. 386–407.

Публикации

1. Баранов С.В., Нарто К., Шебалин П.Н. Modeling and prediction of aftershock activity Surveys in Geophysics, - (год публикации - 2022)

2. Жанг Ю., Ву З., Романелли Ф., Вазари Ф., Джианг Ц., Гао С., Кособоков В., Панда Дж. Next-Generation EEW Empowered by NDSHA: From Concept to Implementation Geosciences, 11, 473, 21p. (год публикации - 2021) https://doi.org/10.3390/geosciences11110473

3. Некрасова А.К., Кособоков В.Г. Uniﬁed Scaling Law for Earthquakes: Lake Baikal Region A. Kosterov, N. Yu. Bobrov, E. Gordeev, E. Kulakov, E. Lyskova, I. Mironova (Eds). Problems of Geocosmos – 2020. Springer Proceedings in Earth and Environmental Sciences, Springer International Publishing, - (год публикации - 2022)

Аннотация результатов, полученных в 2022 году
В третий год выполнения проекта на сайте Института теории прогноза землетрясений и математической геофизики РАН, в котором выполнялся проект, был создан раздел, посвященный результатам работ по проекту https://www.itpz-ran.ru/ru/dejatelnost/current-research/project-rnf-20-17-00180/. В разделе размещены геоинформационные материалы, являющиеся результатом работ по проекту. Проведена подготовка калиброванного каталога землетрясений региона Алтай-Саяны-Прибайкалье с приведением магнитуды к единой шкале. Использованы источники данных как региональных сетей Единой геофизической службы РАН (Алтай-Саянский, Байкальский, Якутский регионы и регионы Приморья), так и данные мирового каталога ISC. Была проведена работа по объединению этих каталогов и идентификации дублей. С этой целью был использован метод на основе модифицированного метода ближайшего соседа вероятностной метрики в пространстве сетевых погрешностей определения эпицентров и времен сейсмических событий, разработанный ранее в рамках проекта РНФ 21-77-30010. В объединенном каталоге произведено удаление данных о взрывах. В результате объединения оказалось, что данные российских сетей и каталог ISC в значительной степени дополняют друг друга. Данные более чем о половине землетрясений содержатся либо только в каталоге ISC, либо только в российских региональных каталогах. Была проведена унификация магнитуды в объединенном каталоге на основе полученных регрессионных соотношений. Процедура создания калиброванного каталога детально описана в документе по ссылке: https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/calib_cat_report.pdf). Построенный калиброванный каталог землетрясений региона Алтай-Саяны-Прибайкалья доступен по ссылке: https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/S_Siberia_ForUpload.txt. По методу ближайшего соседа Заляпина-Бен-Зиона (2013) было проведено разделение каталога землетрясений на независимые (фоновые) и связанные события (афтершоки, форшоки, рои) (https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/S_Siberia_ForUpload_MainShocks.txt и https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/S_Siberia_ForUpload_Decluster.txt). Проведены локальные оценки представительной магнитуды каталога с помощью мультимасштабного метода [Vorobieva, et al., 2013] https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/Представительная_магнитуда.pdf. Важным результатом в 2022 г. была доработка метода оценки параметров Общего закона подобия землетрясений, обеспечивающего разную разрешающую способностью в местах с разной плотностью сейсмических событий. Идея метода состоит в том, что параметры подсчитываются в кругах фиксированного радиуса, но значения приписываются не центру круга, как это делается обычно, а к среднему положению землетрясений выборки. По новой методике проведена оценка параметров новой модели сейсмического режима для региона Алтай-Саяны-Прибайкалье. Построены карты вариаций параметра закона Гутенберга-Рихтера и продуктивности землетрясений https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/b.pdf, https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/a_3.5.pdf, https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/a_6.pdf, https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/a_8.pdf. и https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/Продуктивность.pdf. Разработана новая методика построения синтетического каталога землетрясений. Основная идея нашего подхода состоит в предположении о том, что закон Гутенберга-Рихтера выполняется на региональном уровне, а возможные локальные отклонения компенсируются при осреднении по всему региону. При моделировании это предположение учитывается таким образом, что времена и магнитуды событий в синтетическом каталоге определяются для всего региона независимо от местоположения в соответствии с региональным законом повторяемости, а их местоположение генерируется как векторная случайная величина с заданным распределением по пространству. Другое важное предположение состоит в том, что сильнейшие землетрясения происходят в ограниченном количестве зон в пределах рассматриваемого региона. В рамках проекта предполагается, что сильнейшие землетрясения приурочены к местам пересечений морфоструктурных линеаментов, что подтверждается фактическими данными. Разделение этих объектов на опасные и безопасные проводится на основе геолого-геофизической информации о них методами распознавания образов. При построении синтетического каталога возможное местоположение эпицентров сильнейших землетрясений ограничено окрестностями пересечений линеаментов, распознанных как опасные. Афтершоки обычно исключаются из рассмотрения, чтобы избежать неправомерного завышения сейсмической опасности в местах недавних сильных землетрясений. С другой стороны, афтершоки вносят существенный дополнительный вклад в общую сейсмичность. В нашем подходе вклад афтершоков учиывается с помощью стохастического моделирования афтершоковых процессов с помощью модели ETAS-e, являющейся модификацией широко известной модели ETAS с учетом экспоненциального закона продуктивности землетрясений, установленного ранее в рамках предшествующего проекта РНФ. По разработанной методике сгенерирован синтетический каталог землетрясений по региону Алтай-Саяны-Прибайкалье с учетом региональных значений закона повторяемости Гутенберга-Рихтера для определения магнитуды и времени фоновых событий, локальных значений для определения гипоцентров и модели ETAS-e для афтершоков https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/baikalSynthCat.txt. Каталог сгенерирован на условный период 40 000 лет. Для проверки соответствия моделей реальным данным разработана методика верификации модели сейсмического режима в виде синтетического каталога землетрясений и определены критерии качества модели. Методика применима как ретроспективно, так и для будущих периодов времени. Методика построена на функции правдоподобия – совместной вероятности конкретных реализаций землетрясений в течение периода теста T. Вероятность каждой из реализаций определяется из моделей. В общем случае модель представляется в элементах (ячейках) пространства и магнитуды на регулярной сетке. Мы используем шаг сетки 0.1 градуса по пространству и 0.1 единицы магнитуды. Фактические реализации модели в течение периода времени T в общем виде представляется с такой же дискретизацией. Совместная вероятность – это вероятность наблюдения w(i,j) землетрясений магнитуды i во всех пространственных элементах согласно модели. Предполагается, что значения в разных ячейках независимы. В этом случае совместное правдоподобие представляет собой произведение вероятностей, рассчитанных в каждой ячейке пространства и магнитуды. Проведена ретроспективная оценка качества сгенерированного синтетического каталога https://www.itpz-ran.ru/wp-content/rnf/rnf-20-17-00180/verification_results_syn_cat_Baikal.pdf. Результаты свидетельствуют о хорошем соответствии модели реальным данным. Были продолжены работы по моделированию волновых форм для сценарных землетрясений с протяженным очагом и M=6, 6.5, 7, 7.5, 8, на удалениях от 50 до 200 км, с помощью стохастического метода моделирования синтетических акселерограмм. Моделирование выполнялось не только в рамках классической модели Аки-Бруна, но и с помощью модели очагового спектра Аки-Бруна-Гусева, с двумя характерными частотами с различным скейлингом. Также по имеющимся литературным данным была уточнена использованная ранее модель ak135, а также учтено частотно-зависимое затухание, что позволило уточнить полученные ранее оценки в инженерном диапазоне частот 1-8 Гц. Опубликовано две статьи в журналах первого квартиля и принята в печать одна статья в журнале индексируемом РИНЦ, WOS, Scopus, опубликована статья в сборнике трудов международной конференции.

Публикации

1. Кособоков В.Г., Панца Дж. Seismic roulette: Hazards and risks Terra Nova, Terra Nova, 34(6), 475-494. (год публикации - 2022) https://doi.org/10.1111/ter.12617

2. Крушельницкий К.В., Маточкина С.Д, Селюцкая О.В. Spatio-temporal variations of earthquake productivity in Kamchatka, Kuril Islands and Japan Proceedings of the Third European Conference on Earthquake Engineering and Seismology - 3ECEES, c. 3949-3956 (год публикации - 2022)

3. Скоркина А.А. Modeling of Strong Ground Motion within the Baikal Rift Zone: the Irkutsk case Russian Journal of Earth Sciences (RJES), - (год публикации - 2023)

4. Шебалин П.Н., Баранов С.В, Воробьева И.А. Earthquake Productivity Law in a Wide Magnitude Range Frontiers in Earth Science, Frontiers in Earth Science, Т. 10. С. 881425. (год публикации - 2022) https://doi.org/10.3389/feart.2022.881425

Возможность практического использования результатов
не указано