Перспективы и темпы развития, безопасность населения, крупных сооружений гражданского и промышленного назначения Байкало-Монгольского региона непосредственно зависят от уровня сейсмической опасности исследуемой территории [New Catalog…, 1977], которая в основном оценивается в 7–9 баллов [Khilko et al., 1985]. В геоморфологическом отношении территория представляется как система межгорных понижений и горных сооружений. В геологическом плане район входит в систему Байкальской рифтовой зоны (БРЗ), ее юго-западный и южный фрагмент (рис. 1).

В структурном отношении юго-западный фрагмент включает систему Тункинских впадин и впадины Северной Монголии. Во фрагментах разломов, контролирующих впадины, зафиксированы зоны возникновения очагов землетрясений с магнитудой до 7.0 – это Мондинское землетрясение (1950 г.), Култукское землетрясение (2008 г.), Хубсугульское землетрясение (2021 г.), Быстринское и Кударинское землетрясения (2020 г.). Для них сейсмические сотрясения в эпицентре оцениваются в 8–9 баллов [Gileva et al., 2020; Tubanov et al., 2021; Emanov et al., 2022].

Южный фрагмент БРЗ включает котловину оз. Байкал и продолжение горного сооружения Хамар-Дабан до долины р. Селенги. Наличие активных в сейсмическом плане разломов и высокий ранг проявившихся здесь землетрясений с магнитудой до 7.5 соответствуют высокому показателю сейсмической опасности фрагмента [Chipizubov et al., 2003]. Присутствие сейсмических событий в настоящее время (рис. 1) указывает на то, что интенсивный сейсмический процесс продолжается.

В целом, оценка сейсмической опасности грунтовых слоев различного состава и состояния в рамках проблемы предсказания интенсивности землетрясения определяется следующими основными факторами: влиянием потенциальной сейсмичности зон возникновения очагов землетрясений (ВОЗ), особенностью трассы сейсмических волн, разломами и другими структурными неоднородностями. Инженерно-геологические условия на территории исследования существенно изменяют пространственную структуру волнового поля землетрясения. Достаточно отметить, что только различия в грунтовых условиях могут определить изменение амплитудного уровня колебания почти на порядок [Pavlov, 1988; Kalinina et al., 2017]. Наименьшие искажения в сейсмический сигнал вносят скальные грунты. По этой причине они наиболее выгодны как эталонная отсчетная точка для анализа пространственных вариаций интенсивности землетрясений. Но и само состояние скальных грунтов и процессы выветривания существенно ограничивают пределы точности, с которой могут быть описаны динамические характеристики землетрясений [Seismic Ground Properties, 1988].

С учетом сказанного выполненные исследования направлены на развитие методов изучения и прогноза сейсмических воздействий сильных землетрясений Байкало-Монгольского региона. Конкретной задачей является оценка влияния скоростных неоднородностей приповерхностных рыхлых отложений и глубинных слоев коренных «эталонных» грунтов, представленных сейсмогрунтовыми моделями, на прогнозируемые сейсмические воздействия сильных землетрясений региона.

Исследование основывалось на использовании уже известных методов и подходов, разработанных другими авторами и нами ранее в процессе оценок поведения грунтов различного состава и состояния при заданном уровне сейсмических воздействий [Dzhurik et al., 2015b, 2023]. Выбранные подходы к анализу поведения грунтовых слоев при сейсмических воздействиях излагаются в процессе обсуждения результатов работы, или на них даются соответствующие ссылки. Последовательность реализации поставленных задач связана с необходимостью задания расчетных акселерограмм, построения сейсмогрунтовых моделей и проведения теоретических расчетов.

Задание расчетных акселерограмм или моделирование сейсмических воздействий в настоящее время выполняется на основе использования двух подходов: моделирования воздействия для площадки строительства и моделирования воздействия для сооружения [Nazarov, 2012]. В нашем случае для территорий южной части Восточной Сибири и Северной Монголии реализовано формирование исходных сигналов для прогнозируемых сильных землетрясений, соответствующих принятой исходной сейсмичности, согласно картам общего сейсмического районирования (ОСР). В данном случае для решения поставленной задачи исходный сигнал сформирован для территорий с возможными 8-балльными сотрясениями для средних грунтов. Особенности задания исходного сигнала изложены в работе [Dzhurik et al., 2009], которая основана на использовании как местных записей землетрясений, так и записей мировой сети, достаточно близких к определенным параметрам зон ВОЗ исследуемого района. В развитии этого подхода нами использованы данные относительно сильных землетрясений, произошедших в последние годы (2020–2022) в наиболее опасных зонах ВОЗ, которые зарегистрированы сейсмическими станциями южной части Восточной Сибири (г. Иркутск) и Северной Монголии (г. Улан-Батор).

На территории г. Иркутска для записи землетрясений применялись цифровые инженерно-сейсмометрические станции «Байкал-7HR» с сейсмоприемниками А1638. Характеристики сейсмических каналов используемого комплекта аппаратуры позволяют регистрировать скорость колебаний в диапазоне частот 0.1–125.0 Гц, в диапазоне амплитуд ускорения от 10 мкм/с² до 5 м/с². Данная сейсмостанция установлена сотрудниками ИЗК СО РАН с целью инженерно-сейсмометрического мониторинга сооружений Иркутской ГЭС (52°14'18" N, 104°19'24" E). На территории г. Улан-Батора использовалась аппаратура «Kinemetrics FBA ES-T EpiSensor Accelerometer», которая регистрирует ускорения в м/с² в диапазоне частот меньше 0.01 Гц до 40 Гц без искажений. Сейсмостанция входит в состав глобальной сейсмографической сети (GSN), международный код – ULN (47°51'54.0" N, 107°03'11.0" E).

Для выбранных пунктов регистрации приемники установлены на бетонных основаниях ниже глубины залегания слоя сезонного промерзания, которые опираются, в свою очередь, на крупно- и среднезернистые юрские песчаники (г. Иркутск) и на гранитоиды мезозойского интрузивного комплекса (г. Улан-Батор). Обобщенные значения измеренных скоростей продольных волн (Vp и Vs) в указанных разновидностях грунтов в верхней зоне разреза до 10 м практически сравнимы и находятся в пределах от 1900 до 2200 м/с, а на глубине более 30 м достигают значений от 2800 до 3400 м/с. Скорости поперечных волн (Vs) меняются в интервале от 960 до 1990 м/с.

В первоначальный анализ спектрального состава зарегистрированных землетрясений, необходимого для выбора наиболее представительных записей, включено 17 трехкомпонентных записей землетрясений, зарегистрированных в г. Иркутске и г. Улан-Баторе (табл. 1). Их основные показатели следующие: энергетический класс 10–16, магнитуда 3.5–7.3, интенсивность сотрясения в баллах 2–5 (г. Иркутск) и до 1–2 баллов (г. Улан-Батор), эпицентральные расстояния в пределах 70–324 км (г. Иркутск) и 475–612 км (г. Улан-Батор).

Для первичной оценки расчетных спектров использовались записи землетрясений, зарегистрированных в Иркутске, как наиболее близких по эпицентральным расстояниям и, следовательно, наиболее представительных для выявления связи их параметров с основными параметрами очагов землетрясений. Отмечается, что регистрация землетрясений постоянной сейсмической станцией проводится в пределах города, с высоким уровнем помех. Это обстоятельство учитывалось при анализе спектров.

Расчетные спектры для отобранных акселерограмм до вычитания шума показаны на рис. 2. Они были вычислены с использованием окна длительностью 20.48 с (г. Иркутск) и 81.92 с (г. Улан-Батор) при шаге цифровки 0.01 с. Спектры шума представлены черной линией. Можно отметить общие особенности спектров землетрясений: это быстрое увеличение их значений до 1.0–1.5 Гц, их относительно пологая часть до 3–4 Гц, и далее наблюдается спад до 10 Гц. По интенсивности значения спектра шума закономерно увеличиваются с возрастанием частоты и по величине могут меняться от минимального значения в пределах 1–2 порядков. Проследить изменение максимумов спектров и их частот в зависимости от эпицентрального расстояния и магнитуды по представленным данным не представляется возможным, хотя по приведенным кривым после вычитания шума (табл. 1) отмечается, что для сравнимых по магнитуде землетрясений с увеличением эпицентрального расстояния величина S (значения спектра) снижается и, в связи с большим поглощением высокочастотных колебаний с увеличением эпицентрального расстояния, максимумы спектров могут несколько смещаться в сторону низких частот. Но и для отмеченной связи определенные коррективы может вносить глубина и механизм очагов. В целом, максимальные значения спектров меняются для горизонтальных компонент от 0.006 до 4.030 см/с и для вертикальных – от 0.005 до 2.940 см/с. Частоты основных максимумов спектров для трех компонент находятся в интервале 1.56–8.35 Гц. Общий анализ спектрального состава рассматриваемых землетрясений Байкало-Монгольского региона свидетельствует о том, что вполне обоснованно для формирования исходного сигнала достаточен выбор наиболее сильных и близких по эпицентральному расстоянию землетрясений, зарегистрированных в Улан-Баторе и Иркутске.

Спектры слабых и удаленных отобранных землетрясений по отношению к спектрам «шума» могут различаться в пределах до одного порядка, сильных и относительно сильных землетрясений – на 2–3 порядка (рис. 2). Этот показатель, кроме уже отмеченных ∆, K, М, имеет существенное значение и учитывался при выборе землетрясений для формирования исходного сигнала, в особенности для промышленных и административных центров, какими являются территории г. Иркутска и г. Улан-Батора.

Из приведенных данных (табл. 1) были отобраны пять записей сильных землетрясений, зарегистрированных одновременно в Иркутске и Улан-Баторе на грунтах I категории по сейсмическим свойствам. Отобранные землетрясения характеризуются длительностью колебаний до 20 с (г. Иркутск) и до 80 с (г. Улан-Батор) (рис. 3, а, б), интенсивностью в эпицентре 7.0–8.7 балла, магнитудой (М)>5, энергетическим классом (К)>14 и основными максимумами спектров на частотах 1.61–2.25 Гц (г. Иркутск) и 0.83–1.62 Гц (г. Улан-Батор). Отмечается, что по максимальным ускорениям записи отобранных землетрясений не превышают 10 см/с² для Иркутска и 1 см/с² для Улан-Батора (рис. 3; табл. 2).

Реализован способ, в котором учитываются местные условия формирования исходного сигнала. Основа этого способа опирается на теорию восстановления сигналов и изображений по фазовым характеристикам. В отличие от имеющихся подходов [Field, 2000; Boore, Thompson, 2014], по данным записей землетрясений региональной сети сейсмических станций задается амплитудный и фазовый спектры землетрясения с необходимыми параметрами (для задания амплитудного спектра можно также использовать мировые данные). В итоге формирование исходного сигнала сводится к ряду относительно простых процедур. Поскольку среда выступает формообразующим фактором очагового импульса, фазовый спектр так или иначе учитывает излучение очага и рассеивающие свойства неоднородностей среды.

Расчетные спектры (рис. 3, в, г) нормировались, и рассчитывался средний амплитудный спектр по пяти землетрясениям для NS- и Z-компоненты. Выбор NS-компоненты обусловлен тем, что в инженерной сейсмологии для оценки сейсмических воздействий принято использовать максимальную из двух горизонтальных компонент [Pavlenko, 2009; Aleshin, 2017]. Максимальное значение ускорений приходится на NS-компоненту землетряения № 3 (см. табл. 1). Расчет среднего спектра не является вполне корректной процедурой после отмеченных выше связей параметров спектров с параметрами очагов рассмотренных землетрясений (см. табл. 1), но при таком анализе средний расчетный спектр учитывает особенности колебаний из каждой зоны ВОЗ, для которой зарегистрировано землетрясение, по частотному составу, что, в первую очередь, сказывается на его расширении для выбранных уровней 0.5 и 0.7 от их максимума (табл. 2).

Фазовый спектр находился, по данным записи землетрясения № 1, как наиболее близкий по эпицентральному расстоянию и сравнимый с другими по интенсивности в г. Иркутске (рис. 3, а1), с помощью прямого преобразования Фурье [Jenkins, Watts, 1971]. Следующим этапом является Фурье-синтез расчетного среднего амплитудного спектра сильных землетрясений с полученным фазовым спектром [Dzhurik et al., 2015a]. Расчетные нормированные акселерограммы (рис. 3, д1–2) масштабировались по значениям максимальных амплитуд (рис. 3, е1–2), что не противоречит нормативным документам [Pavlov, 1988].

При таком подходе сформированный сигнал для рассматриваемых компонент отвечает параметрам очагов как близких, так и относительно далеких сильных землетрясений. Максимум спектра для NS-компоненты равен 2.15 Гц и 2.57 Гц – для Z-компоненты (рис. 3, ж, з; табл. 3), и по интенсивности колебаний поверхности, равной 100 см/с², для NS-компоненты и 61 см/с² – для Z-компоненты соответствует уровню сейсмичности для «эталонных» коренных пород, равному 7 баллам, т.е. исходный сигнал отвечает территориям южной части Восточной Сибири и Северной Монголии с исходной сейсмичностью, равной 8 баллам для средних грунтов (II категория по сейсмическим свойствам [Construction…, 2018]), и отражает особенности опасных зон ВОЗ.

В отношении построения сейсмогрунтовых моделей необходимо отметить, что при оценке сейсмических воздействий под приповерхностными неоднородностями понимается изменение скоростей сейсмических волн с глубиной и по площади в грунтовых слоях различного состава и состояния. Например, при сейсмическом микрорайонировании обязательная глубина исследований составляет 10–30 м [Pavlov, 1988], но для расчетов максимальных ускорений, преобладающих периодов и резонансных частот требуется ее увеличение вплоть до залегания «эталонных» коренных пород (приповерхностная скоростная неоднородность) [Dzhurik et al., 2023]. При этом не исключается и изучение глубинных неоднородностей как рыхлых, включая и кайнозойские отложения впадин Байкало-Монгольского региона, так и коренных пород верхней части земной коры практически до глубин возникновения сильных землетрясений в пределах Байкало-Монгольского регион (глубинная скоростная неоднородность). Это особенно существенно для оценки сейсмических воздействий близких сильных землетрясений для промышленных районов юга Сибири и Северной Монголии [Dzhurik et al., 2022а].

Для выявления степени влияния отмеченных приповерхностных и глубинных неоднородностей для рыхлых и коренных пород на сформированный исходный сигнал в первом приближении построено пять сейсмогрунтовых моделей (табл. 3). Для их построения использовались имеющиеся данные бурения, результаты сейсморазведочных и сейсмометрических измерений. Выбор по состоянию и общей мощности рыхлых отложений основывался на анализе проведенных ранее работ [Dzhurik et al., 2009] по обоснованию инженерно-сейсмологических условий строительства в регионе, отмеченных на схеме (см. рис. 1). Учитывались и полученные в последнее время обобщенные данные о скоростях сейсмических волн для наиболее распространенных типов грунтов исследуемой территории [Dzhurik et al., 2022b].

Первая модель соответствует, сверху до 30 м, «эталонным» коренным породам. Для юго-западной части Байкало-Монгольского региона отмечаются широкие пределы изменения скоростей сейсмических волн. В сильнотрещиноватых скальных породах диапазон изменения Vp составляет 1800–2300 м/с, наиболее вероятное их значение в этом состоянии равно 2000 м/с. В слаботрещиноватых, при их залегании до 30 м и более, пределы изменения скоростей от 2800 до 4000 м/с (табл. 3, м. 1). Для анализа влияния приповерхностных неоднородностей на сформированный исходный сигнал первая модель была дополнена сверху рыхлыми отложениями мощностью 10, 30, 90 и 210 м соответственно для моделей 2–5 (табл. 3, м. 2–5). Скорости Vp в рыхлых неводонасыщенных породах меняются от аномально низких 500–600 м/с в верхнем слое до 10 м, до 800–900 м/с – в нижележащих слоях. При этом Vs находится в пределах 240–450 м/с и отношение Vp/Vs – от 1.9 до 2.0. По сейсмическим свойствам это «средние грунты» II категории, к которым относится исходная сейсмичность по картам ОСР [Construction…, 2018].

Одним из основных параметров при расчетах сейсмических воздействий для построенных моделей является затухание сейсмических волн, которое характеризуется декрементом поглощения (D) и оценивается в зависимости от скоростей сейсмических волн согласно модели Гуревича [Gurevich, 1974]. В нашем случае для построенных сейсмогрунтовых моделей при диапазоне изменения Vp от 600 до 4000 м/с декремент поглощения находится в пределах от 0.240 до 0.028 (см. табл. 2) и при увеличении скорости до 7000 м/с уменьшается до 0.014.

Теоретические расчеты основывались на использовании методики и программы расчетов колебаний на поверхности и во внутренних точках среды [Ratnikova, 1973]:

S2(f)=U(f)·S1(f), (1)

где S1(f) – спектр колебаний сейсмической волны, падающей на пачку тонких слоев (рыхлых отложений), U(f) – частотная характеристика рыхлой толщи, S2(f) – спектр колебаний на свободной поверхности.

Необходимо отметить, что при расчетах сформированный на 8-балльные сейсмические воздействия (для средних грунтов) исходный сигнал использовался только для выявления сравнительной оценки его изменения для построенных моделей. При этом предполагалось, что вся энергия исходного сигнала уходит на упругое деформирование вышележащих слоев коренных и рыхлых пород.

Результаты расчетов показывают (рис. 4; табл. 4), что для модели коренных «эталонных» пород максимальное ускорение равно 117 см/с² для NS-компоненты и 42 см/с² – для Z-компоненты, преобладающая частота спектра, соответствующая его максимальным уровням (34.3 см/с – NS и 15.9 см/с – Z), равна 2.15 Гц. Резонансная частота слоя коренных пород мощностью 30 м равна 13 Гц, а уровень частотной характеристики не превышает 2.4 (табл. 4).

Дополненные рыхлыми грунтами модели следующим образом влияют на изменение основных параметров исходного сигнала. Согласно принятому авторами 8-балльному воздействию на рыхлые неводонасыщенные грунты (рис. 4, м. 2–5; табл. 4) максимальные ускорения при изменении мощности рыхлых отложений от 10 до 210 м увеличились и находятся в пределах 166–219 см/с² для NS-компоненты и 68–101 см/с² – для Z-компоненты.

По отношению к сейсмогрунтовой модели «эталона» (табл. 4) ускорение сейсмических колебаний увеличивается в 1.40–1.87 раза для NS- и в 1.61–2.40 раза – для Z-компоненты. Если оценивать по шкале балльности MSK-64, то это будет близко к интенсивности проявления сейсмичности в 8 баллов [Construction…, 2018].

Частоты основных максимумов спектров для неводонасыщенных грунтов для рассматриваемых мощностей находятся в интервале 1.17–4.64 Гц для NS- и в 2.05–2.78 Гц – для Z-компоненты. Но сам интервал частот на уровне 0.7·Smах несколько смещается в сторону низких частот и имеет тенденцию к уменьшению с увеличением мощности рыхлых отложений, что не всегда соблюдается из-за уровня и формы частотной характеристики: для моделей 4 и 5 имеется большее, чем для модели 3, количество резонансных пиков, растянутых в частотной области.

Показательными в этом отношении являются изменения параметров частотных характеристик. При увеличении общей мощности рыхлых отложений в интервале от 10 до 210 м для выбранного авторами эталона резонансная частота закономерно уменьшается для их неводонасыщенного состояния от 5.57 до 0.49 Гц (рис. 4, г; табл. 4). Для максимального значения уровня частотной характеристики такой закономерности в зависимости от общей мощности слоя не наблюдается.

Представленные особенности влияния приповерхностных неоднородностей рыхлых грунтов на исходный сейсмический сигнал рассматриваются авторами далее в сочетании с глубинными для коренных пород, также с учетом их скоростной неоднородности, вплоть до возможных глубин возникновения очагов землетрясений [Dzhurik, Drennov, 1991]. Последние для исследуемого региона регистрируются на глубине порядка 10 км [Golenetskii, 1997], и основная часть гипоцентров землетрясений расположена в пределах от 10 до 25 км [Radziminovich, 2022]. Для юго-западной части Байкало-Монгольского региона построена одна модель, соответствующая в пределах промышленных районов Иркутска и Улан-Батора распределению скоростей сейсмических волн до глубины 14 км. Для уточнения скоростного разреза использовались сведения, приведенные в работе [Zorin et al., 1994]. Для российской части трансекта (см. рис. 1) использовались данные глубинного сейсмического зондирования [Puzyrev, 1981], для монгольской – данные метода поверхностных волн по наблюдениям стационарных станций России и Монголии. По этим данным толщина земной коры варьируется в пределах 40–42 км и является типично континентальной вдоль всего трансекта (см. рис. 1). Другая особенность земной коры связана с постоянством средней скорости продольных волн, равной 6.4±0.1 км/с. В пределах Монголо-Амурской складчатой системы скорость сейсмических волн с поверхности фундамента до глубины 20–30 км увеличивается от 4–5 до 6.8 км/с. При этом по значениям Vp выделяются от четырех до шести слоев. На границе Мохо скорость достигает 8.1 км/с.

С учетом отмеченных данных о скоростях сейсмических волн построенная скоростная модель коренных пород (рис. 5, а, м. 1-6) сверху до 330 м представлена распределением с глубиной скоростей соответствующих модели коренных эталонных пород (см. табл. 3, м. 1). Ниже этого слоя модель 1 дополнена распределением скоростей Vp и Vs до возможных глубин гипоцентров близких землетрясений. Это распределение представлено в виде модели 6 (см. табл. 3). В итоге ниже 330 м до глубин 3–4 км скорость Vp увеличивается от 4 до 6 км/с [Dzhurik et al., 2009] и на глубине 14–15 км достигает значений 6.8–7.0 км/с.

Для сравнительного анализа влияния приповерхностных и глубинных неоднородностей на сформированный единый исходный сигнал (см. рис. 3) построенная модель (рис. 5, а, м. 1-6), в свою очередь, была дополнена сверху моделями рыхлых отложений мощностью 10, 30, 90 и 210 м. Сформированные таким способом модели 1-6, 2-6, 3-6, 4-6 и 5-6 представлены на рис. 5, а, а их цифровые значения сведены в табл. 3.

Отмечается, что при расчетах сформированный исходный сигнал (см. рис. 3) использовался для сравнительной оценки его изменения для построенных моделей и, следовательно, для относительного изменения параметров сейсмических воздействий для возможных территорий строительства Байкало-Монгольского региона, отвечающих построенным моделям. Предполагалось, что вся энергия исходного сигнала уходит на упругое деформирование вышележащих слоев коренных и рыхлых пород.

Результаты расчетов показывают (табл. 5), что для модели коренных «эталонных» пород максимальное ускорение равно 110 см/с² для NS- и 25 см/с – для Z-компоненты, преобладающая частота спектра, соответствующая их максимальному уровню (44.6 и 10.5 см/с), равна 1.76 и 2.34 Гц. Эти показатели свидетельствуют о преобладании низкочастотных колебаний коренных пород, представленных моделью 1-6, что подтверждается частотой основного максимума частотной характеристики, близкой к 1 Гц (табл. 5, м. 1-6).

Дополненная рыхлыми грунтами модель следующим образом влияет на изменение основных параметров исходного сигнала. Максимальные ускорения при изменении мощности рыхлых отложений от 10 до 210 м увеличились и находятся в пределах 146–207 см/с² для NS- и 29–48 см/с² –для Z-компоненты. Но определенная зависимость от мощности рыхлых грунтов по представленным данным начинает сказываться с 30 м, при этом Аmах снижается с 207 до 146 см/с² для NS-компоненты. Для Z-компоненты изменения незначительны (38–46 см/с²) и явного снижения не выявляется.

Такая же закономерность отмечается и для максимальных значений спектра. Это связано с затуханием сейсмических волн. Частота основного максимума спектра меняется незначительно – от 1.51 до 2.34 Гц для NS- и от 1.03 до 3.61 Гц – для Z-компоненты. Частотные характеристики в интервале изменения мощностей рыхлых отложений от 10 до 210 м для построенных моделей (рис. 5, г; табл. 5) показывают, что резонансная частота закономерно уменьшается c увеличением мощности рыхлых отложений от 4.83 до 0.49 Гц. В то же время наблюдается зависимость увеличения значений уровня частотной характеристики с увеличением общей мощности слоя рыхлых грунтов от 10 до 90 м. На 90 и 210 м значения сопоставимы.

Непосредственная оценка различия влияния слоев рыхлых неводонасыщенных грунтов для построенных сейсмогрунтовых моделей представлена для наглядности на рис. 6 в виде частотно зависимых функций U(f). Они получены согласно выражению:

U(f)=Ui(f)/U1(f), (2)

где Ui(f) – расчетные частотные характеристики для моделей 2-6 – 5-6 и U1(f) – расчетная частотная характеристика для модели эталонных коренных пород 1-6 (см. рис. 5, а).

Частоты основных максимумов относительных частотных характеристик несколько увеличились по сравнению с резонансными частотами для рассматриваемых мощностей рыхлых отложений (от 4.83 до 0.49 Гц), но не более чем на 1 Гц (от 5.8 до 0.5 Гц), и это различие уменьшается с увеличением общей мощности рыхлых слоев. Основные максимумы относительных характеристик также сдвигаются в сторону низких частот и меняются от 2.5 до 3.8, но их вклад от общего воздействия значительный. По отношению к максимумам частотных характеристик для построенных моделей, которые меняются от 3.1 до 10.2 (табл. 5), он может быть меньше в 1.24–2.68 раза.

Результаты проведенного анализа влияния неоднородностей грунтовых слоев верхней части земной коры на единый обоснованный исходный сигнал по данным теоретических расчетов приводят к выводу о снижении влияния неоднородностей верхней зоны разреза и увеличении влияния структурных неоднородностей на интенсивность сейсмических колебаний с уменьшением их частоты. Тем самым показана необходимость учета при районировании сейсмической опасности территорий приповерхностных и глубинных скоростных неоднородностей рыхлых отложений и коренных пород.

Дан общий анализ спектрального состава произошедших сильных землетрясений Байкало-Монгольского региона в 2020–2022 гг. по максимальным значениям спектров и соответствующим им частотам, в том числе и на различном уровне вероятности. Это позволило при использовании отобранных землетрясений, одновременно зарегистрированных в Иркутске и Улан-Баторе, по изложенной выше методике для грунтов первой категории впервые сформировать единый исходный сигнал с учетом вероятных очаговых зон для прогнозируемых близких и относительно далеких сильных землетрясений региона.

В соответствии с поставленными задачами для построения расчетных сейсмогрунтовых моделей использованы результаты сейсморазведочных и сейсмометрических измерений, выполненных авторами в последнее время, и данные, полученные ранее другими исследователями. Учитывались также имеющиеся обобщенные сведения об изменении скоростей сейсмических волн с глубиной для наиболее распространенных типов рыхлых неводонасыщенных грунтов до первых сотен метров и коренных пород до возможных глубин возникновения землетрясений. Построенные модели характеризуются мощностью слоев, изменением скоростей продольных и поперечных волн с глубиной, объемной массой и декрементом затухания.

В итоге при использовании сформированного исходного сигнала и построенных сейсмогрунтовых моделей выявлены особенности влияния скоростных неоднородностей приповерхностных и глубинных грунтовых слоев отдельно и в их сравнении, представленные авторами по основным показателям сейсмических воздействий: акселерограммам, спектрам, частотным характеристикам и их параметрам. Полученные расчетные величины сейсмических воздействий в первом приближении могут быть масштабированы для территорий, соответствующих построенным сейсмогрунтовым моделям, расположенным в пределах бассейна р. Селенги и на юго-западной части Байкальской сейсмической зоны, с исходной сейсмичностью, равной 7 и 9 баллам. Авторы также отмечают, что принятие сформированного исходного сигнала для всех 8-балльных территорий южной части Восточной Сибири и Северной Монголии является дискуссионным вопросом на данном этапе исследований, а сама методика может в дальнейшем совершенствоваться, в том числе и за счет накопления статистических данных по сильным землетрясениям региона.

Все авторы внесли эквивалентный вклад в подготовку рукописи, прочли и одобрили финальную версию перед публикацией.

All authors made an equivalent contribution to this article, read and approved the final manuscript.

Авторы заявляют об отсутствии конфликта интересов, связанного с этой рукописью.

The authors declare that they have no conflicts of interest relevant to this manuscript.