Несогласное залегание слоев в верхней части геологического разреза – явление широко распространенное и давно замеченное в практике сейсмического микрорайонирования (СМР). Факторы, влияющие на формы рельефа подземных границ, принято делить как на эндогенные, протекающие внутри Земли, так и на экзогенные – внешние по отношению к грунтовому массиву и проявляющиеся через влияние внешних оболочек, окружающих твердую Землю. Другими словами, рельеф, в том числе и подземный, является результатом совместного действия эндогенных и экзогенных процессов. Для СМР вопросы учета влияния неровностей рельефа на сейсмическое волновое поле имеют важное значение, поскольку наиболее популярные компьютерные программы основаны на двух предположениях: плоскопараллельной модели грунтовых слоев и нормального падения сейсмических колебаний из нижнего полупространства. Справедливость второго предположения основана на том, что большинство землетрясений происходит на значительной глубине, где значения скоростей упругих волн значительно превышают соответствующие значения в грунтовых слоях. Это обстоятельство является причиной того, что угол выхода сейсмической радиации к дневной поверхности, в соответствии с законом Снеллиуса, близок к нулю. В отношении предположения о плоскопаралелльной модели грунтовой толщи такая определенность зачастую отсутствует, что делает сомнительным использование вышеупомянутых программ для правильной интерпретации влияния грунтового массива на соответствующее волновое поле. Как быть в таком случае? Конечно, идеальным было бы использовать усложненные программы обработки, адекватные реальной геологической обстановке, например с использованием метода конечных элементов. Но, во-первых, такие программы в практике СМР используются крайне редко, во-вторых, использование сложных и трудоемких программ обработки для решения частных вопросов СМР зачастую нерентабельно. Сказанное обусловливает поиск простых методических приемов, позволяющих получать приемлемые по точности и трудозатратам решения задач СМР с использованием предположений о плоскопаралелльности границ грунтовой толщи. Описанию одной из возможных методик учета неровностей подземного рельефа посвящена данная работа.

В статье представлены результаты сейсмического микрорайонирования на площадке строительства сосредоточенного объекта повышенной ответственности в г. Хабаровске. Исходная сейсмичность определялась по картам ОСР-2015-В. Выбор уровня карты ОСР был обусловлен ответственностью объекта согласно требованиям нормативных документов. Схема расположения объекта изысканий в районе г. Хабаровска приведена на рис. 1, из которого следует, что здание расположено на берегу р. Амур. В отношении назначения и ответственности объект относится к объектам, перечисленным в п. 2 табл. 3 нормативного документа [1]. Горизонтальные размеры площадки строительства составляют примерно 150 м, что позволяет классифицировать ее с точки зрения СМР как сосредоточенный объект [2]. Фундамент здания заглублен на 1,8 м от уровня дневной поверхности.

В соответствии с требованиями нормативных документов объект изысканий в г. Хабаровске относится к повышенному (I) уровню ответственности. Это означает, что сейсмичность территории, на которой предполагается расположить это здание, следует оценивать по карте общего сейсмического районирования ОСР-2015-В. Фрагмент карты ОСР-2015-В территории расположения площадки изысканий приведен на рис. 2.

На рис. 3 по вертикальной оси отложены значения сейсмической интенсивности, а по горизонтальной – логарифм значений соответствующего периода повторяемости Т. Жирная прямая линия показывает соответствующий график прямой наименьших квадратов отклонений от заданных значений интенсивности. В виде формулы эта зависимость выражается как ΔI = 1,513lgT + 1,6789. При подстановке Т = 1000 лет (уровень ответственности В) с округлением до 0,1 балла получим значение в баллах ΔI = 6,2. Это значение интенсивности, относящееся к средним грунтам 2-й категории, принимаем за исходную сейсмическую интенсивность площадки изысканий.

Определение приращений сейсмической интенсивности на площадке строительства объекта повышенной ответственности в данной работе основано на положениях нормативного документа СП 283.1325800.2016 [2]. В работе [3] изложен новый метод сейсмической микрорайонирования – метод решения прямой задачи СМР (МПЗ).

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

Расчетная мощность грунтового массива была определена в соответствии с нормативными требованиями: глубиной от дневной поверхности до грунта, с сейсмической жесткостью не менее 2000 тм-2 с-1.

Интенсивность исходных сейсмических воздействий относится к референтным грунтам с сейсмической жесткостью не менее R = 2000 тм-2 с-1. В случае, когда интенсивность сейсмических воздействий выражается в баллах, значения исходной сейсмичности берутся по карте ОСР с уменьшением на 1 балл.

Расчет приращения сейсмической интенсивности производится по формуле:

ΔI = 3,3lg (bR0/(R1+R0)),

где R – сейсмическая жесткость грунта;

индексы 0 и 1 показывают отношение их к референтному (в данном случае скальному) или исследуемому грунту;

b – максимум коэффициента динамичности.

Коэффициент 3,3 вместо 2,5, который приведен в нормах [2], использован, поскольку все расчеты по СМР, согласно документу [5], должны производиться в шкале МSK-64 в соответствии с данными по карте ОСР.

В состав работ, относящихся к СМР на площадках строительства, согласно нормативу [2], обязательно входит несколько таблиц:

– таблица физико-механических свойств грунтов. В ней для каждого инженерно-геологического элемента приведены основные данные – плотность, скорость сдвиговых волн, модули деформации, сцепление, сопротивление сжатию и т. п.;

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

Для инженерных изысканий на площадке строительства было пробурено 14 скважин глубиной до 39 м. План расположения скважин на площадке приведен на рис. 4.

По предварительным результатам камеральной обработки полевых работ и лабораторных исследований грунтов на рассматриваемом участке в соответствии с требованиями нормативных документов в разрезе на площадке изысканий выделено 13 инженерно-геологических элементов.

Основные физико-механические параметры грунта площадки изысканий, существенные для СМР, приведены в табл. 1. В соответствующем столбике дано краткое описание прослоя грунта, далее последовательно приводятся значения плотности ρ модуля общей деформации водонасыщенного грунта E, прочности на одноосное сжатие σ и скорости поперечных волн Vs. Скорость поперечных волн – важнейшая характеристика сейсмических свойств грунта. Обычно она определяется при сейсморазведочных работах.

В условиях слабой дифференциации свойств грунтового массива, частых переслаиваний, отсутствия резких сейсмических границ более надежные данные могут быть получены с использованием корреляционной связи параметров инженерно-геологических и сейсмических свойств каждого инженерно-геологического элемента [6].

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

Ввиду вышесказанного, наряду с непосредственным определением значений скоростей поперечных волн по данным сейсморазведки, авторы оценивали значения скоростей Vs, используя связи деформационных и упругих характеристик грунтов. В работе [6] приводятся графики связи динамического модуля упругости и модуля деформации для дисперсных грунтов.

Для дисперсных грунтов эта связь с достаточной для целей СМР степенью точности может быть выражена соотношением E = 20Eдеф, где Eдеф – модуль деформации, а E – динамический модуль упругости. Отсюда, используя выражение связи модуля упругости E и модуля сдвига G:

E = 2 (1 + σ) G,

где σ – коэффициент Пуассона.

Подставив G = ρVs2, можно получить выражение для Vs:

Подставляя в это выражение значения Едеф и ρ для соответствующих ИГЭ из табл. 1, можно получить значения Vs для каждого ИГЭ. Значения коэффициента Пуассона для дисперсных грунтов с большой степенью уверенности можно брать равными 0,3. Для скальных пород коэффициент Пуассона берется равным 0,25.

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

Значение скорости сдвиговых волн рыхлых грунтов было принято на основании сейсморазведочных определений и пересчета из значений модуля деформаций по штамповым испытаниям, а значение скорости Vs алевролитов вычислено из данных по прочности на сжатие с использованием соотношения:  где σсж – предел на одноосное сжатие. При этом если выразить σсж в Па, а ρ в кг/м 3, получим Vs в м/с. Отметим также, что в расчеты принималось значение σсж, относящееся к замоченным грунтам.

Перейдем к общей характеристике особенностей инженерно-геологического разреза на площадке изысканий. Анализ строения грунтовой толщи по скважинам показал, что всю совокупность скважин можно разделить на две группы. В одну можно отнести 9 скважин, которые вскрывают горизонт скальных пород ИГЭ 10, а в другую – 5 скважин, которые не вскрывают горизонт ИГЭ 10.

В качестве иллюстрации особенностей инженерно-геологического разреза по первой группе скважин на рис. 5 представлен разрез по линии II–II.

Как видим, граница раздела между слоями ИГЭ 8 и ИГЭ 10 выражена достаточно плавной линией, хотя в верхней части разреза отмечается чередование слабоконтрастных слоев и видны выклинивания различных слоев.

Разрез по линии IV–IV, представленный рис. 6, является иллюстрацией ситуации, когда высокоскоростной слой кремнистых сланцев ИГЭ 10 отсутствует. В целом картина примерно та же, что и на рис. 6: также заметны выклинивания слоев в верхней части разреза, но только отсутствует горизонт ИГЭ 10.

Согласно нормативному документу [2], если в 30-метровом массиве грунта отсутствует слой с сейсмической жесткостью более 2000 тм-2 с-1, следует использовать гипотетический слой с упомянутыми параметрами на глубине, большей 30 м. Без ограничения общности предположим, что данный слой располагается на глубине 40 м. После этих пояснений можно перейти к описанию построения моделей грунта и расчету соответствующих частотных характеристик.

Как было сказано выше, вся совокупность скважин, описывающих инженерно-геологические свойства грунта на площадке строительства, была разделена на две группы в зависимости от того, выделяется в разрезе слой достаточно жестких грунтов ИГЭ 10 или нет.

Сначала рассмотрим модели первой группы с жестким слоем ИГЭ 10. Глубина кровли этого слоя колеблется в пределах от 2,5 м (скважина 103) до глубины 25,2 м (скважина 234). Инженерно-геологические данные были использованы для построения соответствующих моделей грунтовой толщи с ИГЭ 10, представленных в табл. 2.

Значения плотности и скорости поперечных волн каждого инженерно-геологического элемента брались из табл. 1. Параметры нижнего полупространства определялись свойствами ИГЭ 10.

Частотные характеристики рассчитывались с применением программы NERA [7]. Входное воздействие осуществлялось коротким импульсом амплитудой 0,05g, что соответствует 6-балльному сотрясению. Выходной характеристикой программы NERA, имеющей значение для расчетов приращений сейсмической интенсивности, является максимум спектра реакции при 5 %-ном затухании, нормированного на значение при нулевом периоде, так называемом коэффициенте динамичности. Расчеты показали, что значение максимума коэффициента динамичности b изменяется в пределах от 1,9 (скважина 234) до 2,4 (скважины 201 и 212). Графики коэффициентов динамичности для моделей по скважинам 234 и 201 приведены на рис. 7 и 8. Максимум графика рис. 8 соответствует резонансу на толще грунтов между дневной поверхностью и кровлей слоя ИГЭ 8.

Теперь рассмотрим случай, когда в верхней части грунтового массива слой с повышенной сейсмической жесткостью не выделяется. Именно это имеет место для пяти скважин: 217, 223, 226, 228 и 237. В этом случае норматив [2] рекомендует расчетную мощность определять равной глубине расположения границы с сейсмической жесткостью не менее 2000 тм-2 с-1. Но как это сделать, не проводя бурения, в документе [2] не указано.

В книге [4] было экспериментально показано, что в случае постоянного значения сейсмической жесткости в массиве рыхлого грунта коэффициент динамичности практически не изменяется при увеличении мощности рыхлого массива. Но что будет в случае слоистости грунтового массива, на это упомянутый эксперимент не отвечает.

Проверим на модельном примере, как глубина кровли скального полупространства влияет на частотные характеристики массива, т. е. на величину коэффициента b и соответствующий период. На рис. 9 приведены два графика коэффициентов динамичности – для моделей с глубиной кровли скалистого фундамента 40 и 50 м. При этом средние сейсмические жесткости в верхнем рыхлом слое близки. Коэффициенты динамичности и периоды двух моделей различаются весьма мало, причем меньшей глубине соответствует большая величина и более короткий период.

Параметры моделей грунта по скважинам без слоя ИГЭ 10 приведены в табл. 3. Для нижнего полупространства на глубине 40 м от дневной поверхности приняты значения Vs = 800 м/с, а плотности ρ = 2,5 т/м 3.

Полученные значения коэффициента b для моделей без ИГЭ 10 находятся в пределах 2–2,7 (скважина 217). Можно отметить, что значения коэффициента b для этой группы моделей в целом заметно выше соответствующих значений с более высоким положением кровли скального грунта. Лишь только для одной модели из пяти отмечено низкое значение b = 2. Графики коэффициентов динамичности для скважин 217 и 223 приведены на рис. 10 и 11.

Острый локальный пик на спектре реакции по скважине 217 показывает, что в этом случае имеет место резонанс в верхней части грунтового массива. Точнее всего частоту резонанса можно определить по спектру Фурье, приведенному на рис. 12.

Частота резонанса грунтового массива скважины 217, определенная по максимуму спектра Фурье f = 5,87 Гц, суммарная мощность рыхлых слоев H = 11,6 м. Отсюда по формуле V = 4f × H, подставив значения H и f, можем получить значение скорости V = 270 м/c, что хорошо согласуется со значением средней скорости Vs для грунтового разреза скважины 217. Это, помимо прочего, подтверждает правильность выбора значений скорости Vs для каждого ИГЭ по данным инженерной геологии.

Таким образом, получены значения максимумов коэффициентов динамичности b всех моделей по скважинам на площадке строительства. Они находятся в пределах 1,9–2,7. Эти значения использовались в расчетах по определению приращения сейсмической интенсивности с учетом резонансных эффектов.

Для определения приращения сейсмической интенсивности за счет грунтовых условий была использована методика расчетов по формуле метода прямой задачи СМР.

Приращение балла при этом описывается формулой:

ΔI = 3,3lgbR0/(R1+R0),

где R – сейсмическая жесткость расчетной мощности грунта;

индексы 0 и 1 относятся соответственно к референтному и исследуемому грунту;

b – максимум коэффициента динамичности.

Расчетная мощность грунта соответствует мощности до границы с R > 2000 тм-2 с-1 в случае ее наличия в верхней 30-метровой толще разреза. В противном случае такая граница гипотетически определялась на глубине 40 м. Средние значения плотности, скорости поперечных волн и сейсмической жесткости верхней рыхлой части массива грунта рассчитывались по формулам: ρсред = Σρihi / Σhi; Vсред = Σhi /(Σhi/ vi); Rсред = ρсред × Vсред.

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

Результаты расчетов приращений сейсмической интенсивности и суммарной сейсмической интенсивности на площадке строительства приведены в табл. 4.

В результате проведенных исследований получены следующие основные результаты:

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