ФЕДЕРАЛЬНЫЙ НАДЗОР РОССИИ ПО ЯДЕРНОЙ

И РАДИАЦИОННОЙ БЕЗОПАСНОСТИ

ПОСТАНОВЛЕНИЕ

от 29 декабря 1998 г. N 3

ОБ УТВЕРЖДЕНИИ И ВВЕДЕНИИ В ДЕЙСТВИЕ

РУКОВОДСТВА ПО БЕЗОПАСНОСТИ РБ-006-98 "ОПРЕДЕЛЕНИЕ ИСХОДНЫХ

СЕЙСМИЧЕСКИХ КОЛЕБАНИЙ ГРУНТА ДЛЯ ПРОЕКТНЫХ ОСНОВ"

Федеральный надзор России по ядерной и радиационной безопасности постановляет:

Утвердить и ввести в действие с 1 июля 1999 г. руководство по безопасности РБ-006-98 "Определение исходных сейсмических колебаний грунта для проектных основ".

Начальник

Госатомнадзора России

Ю.Г.ВИШНЕВСКИЙ

Утверждено

Постановлением

Госатомнадзора России

от 29 декабря 1998 г. N 3

Введено в действие

с 1 июля 1999 года

РУКОВОДСТВО ПО БЕЗОПАСНОСТИ

ОПРЕДЕЛЕНИЕ ИСХОДНЫХ СЕЙСМИЧЕСКИХ КОЛЕБАНИЙ ГРУНТА

ДЛЯ ПРОЕКТНЫХ ОСНОВ

РБ-006-98

Руководство содержит общие положения и рекомендации по определению исходных сейсмических колебаний грунта для проектных основ.

СПИСОК СОКРАЩЕНИЙ

ВОЗ - возможный очаг землетрясений

ДСР - детальное сейсмическое районирование

МАГАТЭ - Международное агентство по атомной энергии

МРЗ - максимальное расчетное землетрясение

ОИАЭ - объект использования атомной энергии

ПЗ - проектное землетрясение

СНиП - строительные нормы и правила

СР - сейсмическое районирование

ОСНОВНЫЕ ТЕРМИНЫ И ОПРЕДЕЛЕНИЯ

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

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

Гипоцентр макросейсмический - точка очага, соответствующая максимальной плотности выделившейся сейсмической энергии; определяется по макросейсмическим наблюдениям.

Землетрясение максимальное расчетное (МРЗ) - землетрясение, вызывающее на площадке строительства сотрясение максимальной интенсивности за период 10000 лет.

Землетрясение проектное (ПЗ) - землетрясение, вызывающее на площадке строительства сотрясение максимальной интенсивности за период 100 лет.

Интенсивность (J) сейсмическая - интегральная макросейсмическая мера сейсмического воздействия, определяемая по статистике повреждений эталонных зданий и сооружений, реакции предметов, реакции людей, параметрам движения грунта, параметрам очага землетрясения и расстоянию до пункта наблюдения, изменениям на дневной поверхности.

    Колебание   сейсмическое   грунта   для    проектных   основ -

                                                  ..

сейсмограммы  и спектры  реакции  (акселерограммы X,  велосиграммы

.

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

площадки  воздействиям  при  ПЗ  и  МРЗ  для  50%  и  84%  уровней

обеспеченности.

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

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

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

Продолжительность сейсмических колебаний - время, в течение которого амплитуды колебаний превышают фоновые колебания более чем на 10%.

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

Расстояние гипоцентральное - расстояние от гипоцентра до пункта наблюдения.

Расстояние кратчайшее до поверхности разлома - широко распространенная мера расстояния. Как показала практика, это дает минимальную дисперсию эмпирических данных.

Расстояние эпицентральное - расстояние от эпицентра до пункта наблюдения.

Сейсмичность рассеянная или фоновая - фоновыми землетрясениями считаются землетрясения невысоких магнитуд (обычно с М <= 4), которые трудно связать с известными сейсмогенерирующими структурами. При консервативном подходе считается, что фоновые землетрясения могут произойти в любой точке рассматриваемой территории. Повторяемость землетрясений, определенная для всей территории, в среднем в 5 раз ниже реальной. Таким образом, оценка сейсмической опасности, связанной с фоновыми землетрясениями, сводится к установлению наличия эпицентров микроземлетрясений в непосредственной близости от площадки.

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

    Спектр реакции  (ответа)  -  совокупность  абсолютных значений

                                                     ..

максимальных  амплитуд,  соответственно,  ускорений (у), скоростей

 .

(у),  смещений  (у)  линейно-упругой  системы осцилляторов с одной

степенью  свободы в зависимости от собственной частоты осциллятора

и  затухания 1%, 2%, 5%, 10% критического. Типовой спектр реакции,

приведенный  в  нормах, является огибающей множества нормированных

по  уровню  спектров  и поэтому существенно шире спектра колебаний

отдельно взятого землетрясения.

Спектры нормальные - имеют логарифмическую ширину, близкую к средней: Rs = 0,5 - 0,8 и составляют 1/3 всей выборки мировых данных.

Спектры узкополосные - имеют ширину менее Rs = 0,5 и составляют 1/3 всей выборки мировых данных.

Спектры широкополосные - имеют ширину более Rs = 0,8 и составляют 1/3 всей выборки мировых данных.

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

Ускорение нулевого периода - ускорение спектра ответа в асимптотической (твердотельной) области спектра, которая обычно лежит в диапазоне частот более 33 Гц.

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

Ширина спектра (Rs) - разница логарифмов частот правого и левого склонов спектра на уровне 0,5 от максимального. Логарифмическая ширина спектра слабо зависит от магнитуды землетрясения и в среднем близка к Rs = 0,65. Это свойство уменьшает погрешности экстраполяции слабых воздействий в область сильных.

Эпицентр - проекция точки гипоцентра землетрясения на земную поверхность; характеризуется географическими координатами.

1. НАЗНАЧЕНИЕ И ОБЛАСТЬ ПРИМЕНЕНИЯ

1.1. Руководство по безопасности "Определение исходных сейсмических колебаний грунта для проектных основ" (далее - Руководство) предназначено для определения параметров сейсмических колебаний грунта для проектных основ при обосновании сейсмостойкости ОИАЭ.

1.2. Цель Руководства - изложить возможные подходы к определению исходных сейсмических колебаний грунта для проектных основ, дать оценку приемлемости и установить приоритеты использования этих подходов в приложениях к ОИАЭ различных типов и конкретным целям анализа их безопасности.

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

1.4. Руководство может применяться для анализа действующих, проектируемых и сооружаемых ОИАЭ.

2. СЕЙСМИЧЕСКИЕ КОЛЕБАНИЯ ГРУНТА. МЕТОДЫ И ПОДХОДЫ

2.1. Сейсмические воздействия подразделяются на стандартные и локальные.

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

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

2.2. Сейсмические колебания грунта на площадке зависят от следующих основных факторов:

- положения активных разломов и их параметров (длина, глубина заложения, направление движения, скорость движения);

- положения зон ВОЗ и их параметров (максимальная магнитуда, глубина очага, механизм очага, параметры сейсмического режима);

- удаления площадки от центра активного разлома или зоны ВОЗ;

- характеристики затухания интенсивности сейсмических волн и изменения спектрального состава колебаний на пути распространения колебаний от потенциального очага землетрясения до площадки;

- сейсмических характеристик грунтовых условий площадки (скорости распространения поперечных сейсмических волн, их коэффициентов демпфирования, плотности и мощности слоев грунта).

2.3. Для определения сейсмических воздействий допускается использовать любой из перечисленных ниже методов (подходов) или их комбинаций, которые можно объединить в три основные группы:

I. Методы, использующие записи сильных землетрясений максимального расчетного уровня, имевших место на площадке (подход 1), или имеющиеся аналоговые записи сильных землетрясений (подход 2).

II. Методы, основанные на моделях разлома:

- теоретический метод (подход 3);

- полуэмпирический метод (подход 4).

III. Методы, использующие стандартные спектры:

- методы синтезирования (моделирования, генерации) расчетных акселерограмм и спектров действия с установленными оценками параметров движений грунта при расчетных воздействиях во временной или (и) спектральной форме (подход 5).

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

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

а) следует отдавать предпочтение подходу 1, использующему записи сильных движений от землетрясений на площадке максимального расчетного уровня, поскольку они наиболее удовлетворяют реальной площадке;

б) применение полуэмпирического метода предпочтительно тогда, когда записей сильных движений нет, но имеются данные о параметрах разлома и о распределении скоростей между разломом и площадкой. Использование подхода 4 позволяет получать достаточно надежные результаты;

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

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

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

3. ОБЩИЕ ПОЛОЖЕНИЯ ДЛЯ ОПРЕДЕЛЕНИЯ СЕЙСМИЧЕСКИХ

КОЛЕБАНИЙ ГРУНТА ДЛЯ ПРОЕКТНЫХ ОСНОВ

3.1. Схема алгоритма определения исходных сейсмических колебаний грунта для проектных основ приведена на рис. 1.

┌──────────┐

│ Площадка │

└─────┬────┘

      \/

┌─────────────┐

│  Расчетное  │

│землетрясение│

└──────┬──────┘

       \/

┌─────────────────────┐    ┌─────────────┐    ┌────────────────┐    ┌───────────────────┐

│   Записи сильных    │    │  Параметры  │    │    Параметры   │    │Магнитуда, эпицент-│

│движений на площадке,│нет │  разломов,  │нет │разломов, записи│нет │ральное расстояние,│нет

│   соответствующих   ├───>│распределение├───>│малых сотрясений├───>│грунтовые условия  ├────────┐

│  уровню расчетных   │    │  скоростей  │    │   на площадки  │    │на площадке        │        │

│    землетрясений    │    │             │    │                │    │                   │        │

└──────────┬──────────┘    └──────┬──────┘    └────────┬───────┘    └──────┬────────────┘        │

           │ Да                   │ Да                 │ Да                │ Да                  │

           \/                     \/                   \/                  \/                    \/

┌─────────────────────┐   ┌─────────────┐   ┌────────────────┐   ┌──────────────────┐   ┌──────────────┐

│Использование записей│   │Теоретический│   │Полуэмпирический│   │   Использование  │   │Использование │

│сильных сотрясений на│   │   подход,   │   │    подход,     │   │    стандартных   │   │других записей│

│площадке, соответст- │   │ основанный  │   │ основанный на  │   │    спектров и    │   │    сильных   │

│вующих уровню расчет-│   │ на моделях  │   │моделях разломов│   │продолжительностей│   │  сотрясений  │

│ного землетрясения   │   │  разломов   │   │                │   │                  │   │              │

└──────────────┬──────┘   └───────┬─────┘   └───────┬────────┘   └───────┬──────────┘   └──────┬───────┘

               │                  └────────┬────────┘                    │                     │

               ├───────────────────────────│─────────────────────────────│─────────────────────┘

     ┌─────────┴────┐            ┌─────────┴────────┐           ┌────────┴────┐

     │    Метод,    │            │      Метод,      │           │   Метод,    │

     │ использующий │            │  основанный на   │           │использующий │

     │записи сильных│            │ моделях разломов │           │ стандартные │

     │  сотрясений  │            │                  │           │   спектры   │

     └──────────────┘            └──────────────────┘           └─────────────┘

Рис. 1. Схема алгоритма определения исходных

сейсмических воздействий

В справочных Приложениях 1, 2, 3 к Руководству приводятся соответственно рекомендации, одноименные названиям справочных приложений:

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

- оценка и генерация расчетных сейсмических воздействий методами, основанными на модели разлома;

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

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

I. Определение параметров сейсмического колебания грунта на свободной поверхности грунта (разделы 3 - 5 Руководства);

II. Определение параметров расчетного сейсмического колебания грунта для отметки коренной породы, перекрытой толщей относительно более слабых грунтов (раздел 6 Руководства).

3.3. Решение задачи I этапа предполагает определение:

1) сейсмических колебаний грунта в свободном поле на площадке размещения объекта;

2) частотной характеристики колебаний грунта;

3) максимальных амплитуд (ускорение, скорость, перемещение) колебаний грунта;

4) акселерограмм и спектров ответа.

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

3.5. Если для целей анализа сейсмостойкости требуется использовать два уровня расчетных землетрясений ПЗ и МРЗ, то сейсмические колебания грунта определяются в каждом случае индивидуально и непосредственно для каждого уровня. В этом случае им соответствуют различные наборы колебаний.

Допускается принимать уровень сейсмических воздействий для ПЗ в два раза ниже, чем для МРЗ.

3.6. Исходные сейсмические колебания грунта для проектных основ должны быть получены для свободной поверхности и коренной породы с учетом конкретных сейсмотектонических грунтовых условий площадки.

3.7. Должны быть определены величины параметров сейсмических колебаний грунта и получены оценки их неопределенности. Спектр реакции должен быть определен для затухания 5%.

3.8. Спектр реакции рекомендуется пересчитать для целей анализа сейсмостойкости ОИАЭ для значений затуханий 1%, 2% и 10% с учетом влияния весовых и инерционных характеристик сооружения, размещенного на грунте.

3.9. Набор акселерограмм или акселерограмма для проектных основ должен быть определен в соответствии с требованиями и критериями раздела 5 Руководства.

3.10. Должны быть определены две ортогональные горизонтальные и одна вертикальная компоненты колебаний грунта.

4. ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ СЕЙСМИЧЕСКИХ КОЛЕБАНИЙ ГРУНТА

4.1. Общие положения

4.1.1. Балльность площадки для естественных и техногенно-измененных условий определяется макросейсмической интенсивностью. Рекомендации по ее определению даны в справочном Приложении 4 к Руководству.

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

4.1.3. В случае представления сейсмических колебаний грунта для проектных основ с помощью спектра ответа его форма может быть как стандартной (не зависящей от площадки), так и характерной (зависящей от площадки). Могут быть использованы как широкополосные спектры, отражающие степень неопределенности доминирующих частот исходного сейсмического колебания, так и узкополосные спектры. При использовании узкополосных спектров центральная частота должна сдвигаться (расширяться) на 10% в обе стороны для учета неопределенности положения доминирующих частот. Узкополосные спектры могут быть рекомендованы для конкретных объектов с конкретной геологией при наличии необходимой сейсмологической информации о площадке. Широкополосные спектры следует использовать при разработке типовых проектов, которые могут привязываться к самым различным грунтовым условиям.

4.1.4. Стандартные и широкополосные спектры ответа должны задаваться для различных типов грунтов в виде семейства кривых для набора значений затухания в 1%, 2%, 5%, 10%.

    4.1.5. Следует определять   резонансный   период  спектра  Т ,

                                                                0

который    совпадает   с видимым   периодом   Т,   соответствующим

максимальной амплитуде  колебаний  грунта. Он может быть  вычислен

согласно рекомендациям Приложения 4 к Руководству.

4.1.6. Форму спектра ответа в двойном логарифмическом масштабе в первом приближении можно считать симметричной относительно резонансного периода.

4.2. Использование стандартных спектров ответа

4.2.1. Выбор спектра ответа (его формы) для проектных основ зависит от стадии разработки проекта (типовой проект, для конкретной площадки и т.д.) и от степени полноты сейсмологической информации, имеющейся для площадки строительства ОИАЭ. На стадии разработки типового проекта допускается использовать форму стандартных спектров ответа.

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

4.3. Горизонтальная компонента

4.3.1. В качестве стандартного спектра ответа горизонтальных компонент на свободной поверхности для интенсивности землетрясения 9 баллов по шкале MSK-64 допускается принимать спектр, приведенный на рис. 2, при условии обоснования его пригодности для данного случая применения.

                       График не приводится.

    ┌────────┬───────────────────────────────────────────────┐

    │Чаcтота,│              2                                │

    │   Гц   │Ускорение, м/с , при относительном затухании, %│

    │        ├───────────┬───────────┬───────────┬───────────┤

    │        │     1     │     2     │    5      │    10     │

    ├────────┼───────────┼───────────┼───────────┼───────────┤

    │    1   │    6,0    │    5,0    │    4,0    │    3,0    │

    ├────────┼───────────┼───────────┼───────────┼───────────┤

    │    2   │   26,0    │   20,0    │   13,0    │   10,0    │

    ├────────┼───────────┼───────────┼───────────┼───────────┤

    │   10   │   26,0    │   20,0    │   13,0    │   10,0    │

    ├────────┼───────────┼───────────┼───────────┼───────────┤

    │   30   │    5,0    │    5,0    │    5,0    │    5,0    │

    └────────┴───────────┴───────────┴───────────┴───────────┘

   Рис. 2. Стандартные спектры ответа горизонтальных компонент

     на свободной поверхности для интенсивности землетрясения

                        9 баллов по MSK-64

4.3.2. Спектры ответа для обеих горизонтальных компонент в этом случае могут приниматься идентичными.

4.4. Вертикальная компонента

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

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

Таблица 1

              СООТНОШЕНИЕ МЕЖДУ АМПЛИТУДАМИ УСКОРЕНИЙ

                        ..                   ..

        НА ВЕРТИКАЛЬНОЙ X V И ГОРИЗОНТАЛЬНОЙ X Н КОМПОНЕНТАХ

┌─────────────┬───────┬──────┬──────┬──────┬──────┬──────┬───────┐

│..           │       │      │      │      │      │      │       │

│Х Н, см/кв. с│250    │300   │350   │400   │450   │500   │550    │

│..           │       │      │      │      │      │      │       │

│X V, см/кв. с│125    │155   │185   │225   │275   │335   │400    │

├─────────────┼───────┼──────┼──────┼──────┼──────┼──────┼───────┤

│..           │       │      │      │      │      │      │       │

│Х Н, см/кв. с│600    │650   │700   │750   │800   │850   │900    │

│..           │       │      │      │      │      │      │       │

│X V, см/кв. с│470    │545   │625   │710   │800   │895   │1000   │

└─────────────┴───────┴──────┴──────┴──────┴──────┴──────┴───────┘

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

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

4.4.2. Отношение уровней двух горизонтальных компонент подчиняется логнормальному закону с нулевым средним значением и стандартным отклонением 0,16 ед. логарифма (если взять отношение более интенсивной ко второй горизонтальной компоненте, то отношение амплитуд в 16% случаев не превышает 1,15, в 50% случаев не превышает 1,65).

5. КРИТЕРИИ И МЕТОДЫ ОПРЕДЕЛЕНИЯ АКСЕЛЕРОГРАММ

5.1. Общие положения

5.1.1. Для синтезирования акселерограмм может быть использован обобщенный спектр реакции (медианные или 84% обеспеченности), полученный в результате статистической обработки набора аналоговых акселерограмм, подобранных для сходных сейсмотектонических и грунтовых условий.

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

5.1.2. В случае задания сейсмического колебания грунта в виде спектра ответа акселерограммы должны удовлетворять критериям, содержащимся в разделах 5.2 и 5.3 Руководства.

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

5.2.1. В дополнение к другим требованиям, приведенным ниже, акселерограммы должны быть отобраны, модифицированы либо получены численными методами таким образом, чтобы их временные параметры (длительность акселерограмм, огибающая колебаний) и амплитудные параметры (пиковое ускорение, пиковая скорость, пиковые перемещения) соответствовали тем, которые определены для площадки по ее макросейсмическим характеристикам.

5.2.2. Огибающая колебаний - сглаженная функция закономерного изменения пиковых амплитуд во времени. На рис. 3 приведена форма амплитуды огибающих колебаний.

                      График не приводится.

                   ┌─────────┬───────┬───────┐

                   │Магнитуда│Т  / Т │Т  / Т │

                   │         │ a    c│ b    c│

                   ├─────────┼───────┼───────┤

                   │    8    │ 0,08  │ 0,46  │

                   ├─────────┼───────┼───────┤

                   │    7    │ 0,12  │ 0,50  │

                   ├─────────┼───────┼───────┤

                   │    6    │ 0,16  │ 0,54  │

                   └─────────┴───────┴───────┘

           (0,31М-0,774)

    Т  = 10

     с

    Т  - зависящая от магнитуды  длительность колебаний в секунду,

     с

при Т > Т  амплитуда воздействия не превышает 1/10 максимальной амплитуды;

         с

    Т  - время нарастания амплитуд акселерограммы;

     a

    Т  - время начала спада амплитуд акселерограммы;

     b

    М - магнитуда.

           Рис. 3. Форма амплитуды огибающих колебаний

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

5.2.4. Пиковое ускорение акселерограммы определяется как:

а) максимальное значение ускорений;

б) абсолютное значение векторной суммы 2 горизонтальных и вертикальной компонент.

5.3. Критерии синтезирования акселерограмм

Акселерограммы должны удовлетворять следующим критериям:

5.3.1. Математическое среднее ускорений нулевого периода индивидуальных акселерограмм должно быть больше или равно значению пикового ускорения; каждая последующая частота должна отстоять от предыдущей на величину интервала, равную 10% предыдущей частоты.

5.3.2. В области частот от 0,5 до 33 Гц математическое среднее отношений значений спектра ответа (вычисленного на основе индивидуальных акселерограмм) и спектра ответа для проектных основ (отношения должны быть вычислены для всех частот, указанных в табл. 1) должно быть <= 1.

5.3.3. Ни одна точка математического среднего спектров ответа, вычисленных для индивидуальных акселерограмм, не должна лежать ниже 10% проектного спектра ответа. Значения спектров ответа должны быть рассчитаны с достаточно малым шагом по частоте. Рекомендуемые частоты и интервалы приведены в табл. 2.

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

Таблица 2

ЧАСТОТЫ, РЕКОМЕНДУЕМЫЕ ДЛЯ ВЫЧИСЛЕНИЯ СПЕКТРОВ ОТВЕТА

┌────────────────────────────────┬───────────────────────────────┐

│     Частотный диапазон, Гц     │          Приращение, Гц       │

├────────────────────────────────┼───────────────────────────────┤

│0,5 - 3,0                       │0,10                           │

├────────────────────────────────┼───────────────────────────────┤

│3,0 - 3,6                       │0,15                           │

├────────────────────────────────┼───────────────────────────────┤

│3,6 - 5,0                       │0,20                           │

├────────────────────────────────┼───────────────────────────────┤

│5,0 - 8,0                       │0,25                           │

├────────────────────────────────┼───────────────────────────────┤

│8,0 - 15,0                      │0,50                           │

├────────────────────────────────┼───────────────────────────────┤

│15,0 - 18,0                     │1,0                            │

├────────────────────────────────┼───────────────────────────────┤

│18,0 - 22,0                     │2,0                            │

├────────────────────────────────┼───────────────────────────────┤

│22,0 - 34,0                     │3,0                            │

└────────────────────────────────┴───────────────────────────────┘

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

    Статистическая независимость двух акселерограмм a (t) и  a (t)

                                                     1        2

подтверждается вычислением коэффициента корреляции:

                    Е (а (t) - m ) (a (t) - m )

                        1       1    2       2

             ро   = ---------------------------,

               12               *      *

                           сигма  сигма

                                1      2

    где:

    Е - математическое ожидание;

    m , m  - средние значения a (t) и a (t);

     1   2                     1       2

         *       *

    сигма , сигма  - стандартные отклонения.

         1       2

Две акселерограммы считаются статистически независимыми, если абсолютное значение коэффициента корреляции не превышает 0,3.

6. РАСЧЕТНЫЕ СЕЙСМИЧЕСКИЕ КОЛЕБАНИЯ ГРУНТА ДЛЯ ОТМЕТКИ

КОРЕННОЙ ПОРОДЫ

6.1. Сейсмические колебания грунта для отметки коренной породы определяются на основе анализа взаимодействия грунта и сооружения.

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

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

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

6.5. Вычисления, определяющие колебания на границах в соответствии с п. п. 6.2 и 6.3 Руководства, должны выполняться с использованием математических моделей и процедур, совместных с используемыми при проведении анализа системы "грунт-сооружение".

6.6. В случае частично заглубленных сооружений должны быть проведены расчеты изменения амплитудного и частотного составов сейсмических колебаний по глубине расположения фундамента. Амплитуды ускорений спектра ответа, вычисленного в свободном полупространстве на глубине фундамента, не должны лежать ниже 60% соответствующего проектного спектра ответа на поверхности грунта.

7. ОПРЕДЕЛЕНИЕ ХАРАКТЕРИСТИК ПАРАМЕТРОВ СЕЙСМИЧЕСКИХ

КОЛЕБАНИЙ ГРУНТА ЗАДАННОЙ ВЕРОЯТНОСТИ ПРЕВЫШЕНИЯ

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

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

Приложение 1

(справочное)

к Руководству "Определение

исходных сейсмических колебаний

грунта для проектных основ"

ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ИСХОДНЫХ СЕЙСМИЧЕСКИХ КОЛЕБАНИЙ

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

Этот метод применяется в двух случаях (см. рис. 1 Руководства):

1) когда известны записи сильных землетрясений на площадке или рядом с ней, соответствующих уровню МРЗ или ПЗ;

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

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

   ┌───────────────────┐

   │Соответствующие МРЗ├────────────────────────────────────────────>│

   │  /акселерограммы  │                                      ┌ ─ ─ ─│─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ┐

 1 │  /велосиграммы    │  ┌─────────────┐                            │Параметры сейсмических воздействий

   │  /сейсмограммы    │  │  Механизмы  │                     │      │  ┌──────────────────────────┐    │

   │   на площадке     │  │   очагов    ├──┐                         ├─>│        Расчет            │

   └───────────────────┘  └─────────────┘  │                  │    ┌─┼─>│      _2                  ├──┐ │

                          ┌──────────────┐ │                    ┌──┼─┼─>│а   , а , V   , d   , Т   │  │

                          │  S(омега) -  │ │                  │ │  │ │  │ maх       maх   maх   maх│  │ │

                          │  частотная   │ │                    │  │ │  └──────────────────────────┘  │

                          │характеристика│ │                  │ │  │ │                                │ │

                          │сейсмического │ │                    │  │ │  ┌────────────────────────┐    │

                          │    канала    │ │                  │ │  │ ├─>│         Расчет         │    │ │

                          └───┬──────────┘ │                    │  ├─┼─>│      длительности      ├──┐ │

                              │        ┌───┘                  │ ├──┼─┼─>│     акселерограммы     │  │ │ │

                              \/       \/                       │  │ │  └────────────────────────┘  │ │

                           ┌─────────────┐  ┌────────────────┐│ │  │ │                              │ │ │

   ┌───────────────────┐   │   Решение   │  │Соответствующая │  │  │ │  ┌────────────────────────┐  │ │

   │Соответствующие МРЗ├──>│интегрального│  │МРЗ акселеро-   ├┼─┼─>│ ├─>│         Расчет         │  │ │ │

   │  /акселерограммы  │   │  уравнения  ├─>│грамма,         │  │  ├─┼─>│    спектров реакции    ├─>│ │

 2 │  /велосиграммы    │   │  I рода по  │  │пересчитанная   ││ ├──┼─┼─>│   SA(t), SV(t) SD(t)   │  │ │ │

   │  /сейсмограммы    │   │типу свертки │  │на площадку     │  │  │ │  └────────────────────────┘  │ │

   │   вне площадки    ├─┐ └─────────────┘  └────────────────┘│ │  │ │                              │ │ │

   └───────────────────┘ │       /\                             │  │ │  ┌────────────────────────┐  │ │

                         │       │                            │ │  │ ├─>│     Расчет спектра     │  │ │ │

                         │ ┌─────┴───────┐                      │  ├─┼─>│       мощности РS      │<─┘ │

   ┌────────────────────┐└>│    Расчет   │                    │ ├──┼─┼─>│                        │    │ │

   │Записи средних и    │  │ эмпирических│                      │  │ │  └────────────────────────┘    │

   │слабых землетрясений├─>│ передаточных│                    │ │  │ │                                │ │

   │на площадке и вне ее│  │   функций   │                      │  │ │  ┌────────────────────────┐    │

   └────────────────────┘  │    среды    │                    │ │  │ ├─>│       Расчет           │<───┤ │

                           └─────────────┘                      │  ├─┼─>│    интенсивности I     │    │

                                                                │  │ │  │                        │    │ │

   ┌────────────────────┐  ┌─────────────┐  ┌────────────────┐│ │  │ │  └────────────────────────┘    │

   │Мировой или региона-├─>│   Выборка   │  │Соответствующая │  │  │ │                                │ │

 3 │льный Банки данных  │  │акселерограмм├─>│МРЗ акселеро-   ├─>│  │ │  ┌────────────────────────┐    │

   │сильных движений    │  │ по магнитуде│  │грамма, подоб-  ││ │  │ └─>│   Расчет коэффициента  │    │ │

   └────────────────────┘  │и расстояниям│  │ранная к услови-│  │  └───>│    динамичности бета'  │<───┘

                           └─────────────┘  │ям площадки     ││ └──────>│                        │      │

                                 /\         └────────────────┘          └────────────────────────┘

                         ┌───────┴─────────┐                  └ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ┘

                         │Критерии отбора: │

                         │схожесть глубины │

                         │  и механизмов   │

                         │очагов, схожесть │

                         │грунтовых условий│

                         └─────────────────┘

Рис. 1.1. Блок-схема определения параметров сейсмических

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

1. ИСПОЛЬЗОВАНИЕ ЗАПИСЕЙ СИЛЬНЫХ ЗЕМЛЕТРЯСЕНИЙ

НА ПЛОЩАДКЕ

Если имеются записи сильных движений (землетрясений) на площадке (например, акселерограммы), соответствующие уровню МРЗ, то следует отдать предпочтение этому методу, поскольку в них (в этих записях) учтена вся информация, определяющая исходные сейсмические колебания грунта для проектных основ.

1.1. Акселерограмма паспортизирована следующими данными:

а) обозначения составляющих; б) масштаб времени; в) ускорений в виде отрезков, соответствующих по длине 1 с и 0,1 g; г) строчкой приводятся: балльность, дата землетрясения, эпицентральное расстояние в километрах и номер сейсмостанции (рис. 1.2 - не приводится).

Для обозначения трех составляющих колебаний принято, что буква В означает вертикальную составляющую, а буква Г - горизонтальную с указанием угла в градусах между направлением регистрируемого колебания и направлением к эпицентру (рис. 1.2).

Колебательный процесс является непериодическим колебанием с переменной амплитудой и периодом. Периодом считается удвоенный интервал между смежными нулевыми значениями ускорений.

1.1.1. Начальный отрезок времени обычно имеет относительно меньшие амплитуды. Он относится к продольным сейсмическим волнам. Его интервал тем длиннее, чем больше эпицентральное расстояние. При малых расстояниях от очага начальный отрезок составляет 1 - 3 с. Периоды колебаний в начальном участке относительно меньше.

1.1.2. Средний участок записи имеет наибольшие амплитуды ускорений. Он относится к поперечным волнам, осложненным вступлениями обменных и поверхностных волн. Периоды колебаний на среднем участке немного больше или такие же, как и на начальном участке. Переход от начального участка к среднему ясно выражен на записи.

1.1.3. Конечный участок отличается более длинными периодами. Амплитуда ускорений постепенно и нерегулярно уменьшается так, что трудно фиксировать окончание колебаний. Переход от среднего участка к конечному явно не выражен.

1.1.4. Общая продолжительность колебательного процесса не одинакова и тем больше, чем больше балльность и эпицентральное расстояние. Колебание продолжается приблизительно 10 - 40 с.

1.1.5. Число отклонений (амплитуд) на записи очень велико (более 100), что учитывается при спектральном анализе.

1.1.6. Вертикальная составляющая ускорения обычно имеет несколько меньшие амплитуды, чем горизонтальные (около 60 - 70%).

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

1.2. По акселерограммам можно определить:

- амплитудный уровень колебаний;

- продолжительность колебаний;

- спектральный состав колебаний.

1.2.1. Амплитудный уровень колебаний

По акселерограммам амплитудный уровень колебаний может быть задан:

- максимальной амплитудой;

- максимальным размахом колебаний;

- среднеквадратической амплитудой колебаний, равной корню квадратному из суммы амплитуд на представительном участке записи;

- спектральным уровнем на фиксированных периодах.

                                                             а

                                                              maх

    Для подавляющего большинства  акселерограмм  соотношение ----,

                                                              _2

                                                              а

                                                   _2

где а    -  амплитуды    максимального ускорения и а   -  величина

     mах

среднеквадратического   ускорения,   лежит  в  интервале  значений

2,41 - 2,47 [13].

    Следует иметь в виду, что с инженерной  точки  зрения величина

а    может  быть  не  представительна  в случае единичного выброса

 maх

большой  амплитуды  либо,  когда  период  колебаний с максимальной

амплитудой лежит вне  пределов  максимума частотной характеристики

колебаний сооружений. Вследствие  этих,  а  также некоторых других

факторов корреляции между а      и   интенсивностью сопровождаются

                           maх

значительным   разбросом,   близким   к величине среднего значения

(дельта  -> 100%). Однако в сочетании с длительностью максимальной

       а

фазы  величина а     х  тау  имеет  более   тесную  корреляцию   с

                mах

балльностью [14]. В качестве амплитудного  параметра  сейсмических

воздействий используются также амплитуды спектральных ускорений SA

колебаний осциллятора с 5% затуханием на   фиксированных  периодах

Т = 0,2 с  и  Т = 1,0 с,  описывающих  воздействия  в  коротко-  и

длиннопериодной областях спектра [15].

    1.2.2. Продолжительность колебаний

    Параметром   продолжительности   (длительности)   колебаний на

практике   используется   ширина  импульса  тау,  т.е.  промежуток

времени,   в   течение   которого   уровень   колебаний  превышает

а    / 2.   Схема  измерения амплитуды и ширины импульса  показана

 maх

на  рис. 1.3 (не приводится). Представлены случаи простых по форме

колебаний (рис.  1.3, "а") и запись волновой  группы,  разделенной

промежутком  времени  ДЕЛЬТА t, в  течение  которого а <= а    / 2

                                                           maх

(рис. 1.3, "б").  В  случае  сложных  колебаний  типа  (б) импульс

считается   единым  при  ДЕЛЬТА t <=  2 с.  Тогда  t = ДЕЛЬТА t  +

                                                               1

ДЕЛЬТА t  +  ДЕЛЬТА t .  В  противном   случае  (ДЕЛЬТА t  >  2 с)

                     2

рассматриваются два отдельных импульса.

1.2.3. Спектральный состав колебаний

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

- период, соответствующий максимальной амплитуде;

- спектр Фурье S(омега);

- спектры реакции, т.е. спектры ускорений SA(T), скорости SV(T) и смещений SD(T) осцилляторов с 5%-ным затуханием.

Выражения для спектров Фурье записываются в виде [16]:

                         +беск.        i омега t

             S(омега) = интеграл f(t) e          dt,         (1.1)

                         -беск.

    где:

    f(t) - колебательный  процесс;

    омега - круговая частота;

    S(омега) - комплексная   функция  круговой  частоты  полностью

                                                         i( )

определяется  амплитудным │S(омега)│ и фазовым спектром е        :

                                                         фи омега

                                   i( )

            S(омега) = │S(омега)│ е        .                 (1.2)

                                   фи омега

    В свою  очередь  амплитудный и  фазовый  спектры  определяются

через коэффициенты А(омега) и В(омега):

                            2           2        -1/2

             │S(омега)│ = [A (омега) + B (омега)]    ;

                                    В(омега)

                  фи(омега) = аrctg --------,                (1.3)

                                    А(омега)

    где:

                             t

              А(омега) = интеграл f(t) cos oмега t dt;

                             0

                             t

              B(омега) = интеграл f(t) sin oмега t dt.

                             0

    В инженерных оценках  большая часть операций  осуществляется с

использованием   модуля  спектра  |S(oмега)|,   т.е.  амплитудного

спектра Фурье.

    1.2.4. Спектры реакции

    Уравнение   движения  y(t)   линейного   осциллятора  с  одной

степенью свободы  и  затуханием  кси  при  движении основания  (t)

выражается соотношением:

      ..                  .           2         ..

      y(t) + 2 омега  кси y(t) + омега  y(t) = -x(t),        (1.4)

                    i                 i

    где:

    омега  -  частота    собственных    колебаний    незатухающего

         i

осциллятора.  Относительное  смещение  такого  осциллятора y может

быть записано в виде:

                                            -омега кси(t-тау)

                1             t     ..            i                             2 1/2

y = - -------------------- интеграл х(тау) е                  sin омега (1 - кси )    (t - тау) dтау. (1.5)

                     2 1/2    0

      омега  (1 - кси )

           i

                           .

    Относительная скорость y:

                      -омега (t-тау)

.       t     ..            i                          2 1/2

y = -интеграл х(тау) е               cos омега (1 - кси )    (t - тау) dтау +

        0

                                     -омега кси(t-тау)

             кси       t     ..            i                             2 1/2

    + ------------- интеграл х(тау) е                  sin омега (1 - кси )    (t - тау) dтау.        (1.6)

              2 1/2    0

      (1 - кси )

                          ..  ..

    Абсолютное ускорение (y + х):

          омега  (1 - 2 кси)                  -омега кси(t-тау)

 ..  ..        i                t     ..            i                             2 1/2

(y + х) = ------------------ интеграл х(тау) е                  sin омега (1 - кси )    (t - тау) dтау +

                    2 1/2       0

            (1 - кси )

                                   -омега кси(t-тау)

                     t     ..            i                             2 1/2

    + 2 омега кси интеграл х(тау) е                  соs омега (1 - кси )    (t - тау) dтау.          (1.7)

                     0

    Спектром  реакции  является  огибающая  максимальных  откликов

                                                             2 1/2

осцилляторов. При малых величинах затухания кси член (1 - кси )

и выражение для максимального относительного смещения  осциллятора

y = SD могут быть записаны в виде (см. 1.5):

                       -омега кси(t-тау)

         t     ..            i                             2 1/2

S  = [интеграл х(тау) е                  sin омега (1 - кси )    (t - тау) dтау]   .                  (1.8)

 D       0                                                                      maх

    Тогда для максимальной относительной скорости (см. 1.6):

                           SV = oмега  SD.                   (1.9)

                                     i

Спектры реакции строятся в двойном логарифмическом масштабе (рис. 1.4 - не приводится).

Они могут быть рассчитаны прямым методом решения уравнения (1.4) и пересчетом из одного спектра в другой. Например, спектр скорости, рассчитанный умножением на омега спектра смещений или делением на омега спектра ускорений, называется спектром псевдоскорости. Он отличается от определенного прямым способом спектра тем меньше, чем меньше затухание. Значения затухания большинства строительных конструкций лежат в интервале кси = 0,02 - 0,2 [17].

1.2.5. Коэффициент динамичности

    Динамический коэффициент равен [18]:

                                   _2

                           бета  = а    (t),                (1.10)

                               t    бета

    _2

где а    (t) - амплитуда  (огибающая)  записи  ускорений колебаний

     бета

линейного осциллятора   с   заданной  частотой омега  и затуханием

                                                    0

кси    при  воздействии колебаний   x(t). В   инженерных  расчетах

наибольший    интерес    представляет    максимальное     значение

бета     = бета', которое  называется  коэффициентом  динамичности

    tmax

и характеризует  эффект  сейсмического воздействия  на сооружение.

Значение  бета'  зависит  от периода колебаний Т и затухания  кси.

Амплитудный    уровень,    т.е.    масштаб    коэффициента   бета,

устанавливается делением бета' на величину максимального ускорения

а    колебаний основания  осциллятора (ускорения колебаний почвы).

 mах

Так же используются кривые зависимости динамического  коэффициента

от периода колебаний бета(Т)  для фиксированных затуханий  (обычно

кси принимается равным 5%) [4]. Значительно   реже    используются

кривые бета по скоростям и  смещениям. В  сочетании  со  спектрами

реакции динамические коэффициенты по  ускорениям бета ,  скоростям

                                                     а

бета  и смещениям бета  определяются соотношениями:

    v                 d

бета  = SA / a   ; бета  = SV / V   ; бета  = SD / d   . (1.11)

    а         max      v         max      d         max

2. ИСПОЛЬЗОВАНИЕ АНАЛОГОВЫХ ЗАПИСЕЙ СИЛЬНЫХ ДВИЖЕНИЙ ГРУНТА

В случае применения аналоговых записей землетрясений для условий площадки размещения объекта эти записи необходимо пересчитать. Методика пересчета основана на теории распространения плоских волн в горизонтально-слоистых линейно-неупругих средах [20, 21].

    Расчет  акселерограммы  колебаний  для площадки a(t) по записи

u(t),   полученной  на  расположенной  вне  площадки  сейсмической

станции, сводится к решению интегрального уравнения  первого  рода

типа свертки [22]:

                  К(t - тау) а(тау) dтау = u(t),            (1.12)

    где  ядро  интеграла  K(t) описывает влияние строения среды на

пути  "очаг - станция",   "очаг - площадка",   а   также   влияние

регистрирующей аппаратуры. В спектральной области K(t) описывается

соотношением:

    К(омега) = F K(t) = эпсилон эта S(омега) ТЭТА(омега),   (1.13)

    где:

    F - оператор преобразования Фурье;

    S(омега) - частотная характеристика  сейсмографа;

    ТЭТА(омега) - отношение частотных характеристик среды на путях

распространения  волн  по трассе  "очаг - станция" G (омега)  и  к

                                                    i

площадке G (омега);

          2

    эпсилон - коэффициент, зависящий от типа регистратора,

         2

    oмега  - сейсмограф;

    i oмега - велосиграф;

    I - акселерограф;

    эта = ПСИ  / ПСИ ,   где  ПСИ  - направленность  излучения  из

             1      2            1

очага на станцию, ПСИ  - из очага на площадку.

                     2

    Если  механизм  очага  известен,   величина  эта  может   быть

рассчитана.  В противном  случае  величина  К(омега)  определяется

с  точностью  до  величины  этого  коэффициента.  Чтобы  уменьшить

влияние коэффициента эта,  следует выбирать записи как можно ближе

к площадке объекта.

Определение акселерограмм колебаний на площадке по записям на станции относится к классу некорректно поставленных задач, и поиск решения проводится в регуляризированном виде. Применяется специальная процедура регуляризации, заключающаяся в подборе коэффициентов с целью обеспечения регуляризированного решения уравнения (1.12). Математическое обеспечение в виде готовых вычислительных программ дается в работе [23]. Вопрос об отличиях в условиях распространения сейсмических волн сводится к определению значений комплексной функции ТЭТА (омега) выражения (1.11). Практический интерес представляет путь построения эмпирических передаточных функций среды с использованием землетрясений и взрывов. Эти передаточные характеристики строятся в виде отношений спектров сейсмических колебаний.

    Спектры записи  сейсмических  колебаний  в  пункте с номером i

могут быть приближенно представлены в виде:

          U (омега) = фи  W(a) G (омега) S (oмега),         (1.14)

           i            i       i         i

    где   W(oмега)  -  спектр    источника.   Отношение   спектров

синхронных записей одного сейсмического  события на  двух станциях

i  = 1,2 будет равно:

                 U (омега)                      S (омега)

                  1                              1

   R   (омега) = --------- = эта ТЭТА   (омега) ---------.  (1.15)

    1,2          U (омега)           1,2        S (омега)

                  2                              2

    Из этого выражения видно, что  значение  передаточной  функции

среды:

                            -1

              q   (t) = Re{F  [ТЭТА   (омега)]}             (1.16)

               1,2                 1,2

    может быть определено по записям слабых и  средних по величине

сейсмических событий  в  полосе  частот,  на  которой U (омега)  и

                                                       1

U (oмега) хорошо представлены.

 2

    Если в пункте с условным  номером  i = 1,  в  котором  удалось

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

зоны, и в пункте i = 2 на строительной площадке имеются синхронные

        (n)

записи U   (t)  нескольких  более   слабых  сейсмических   событий

        1

(n = 1, 2,...N) из  той же очаговой зоны, то формально связь между

соответствующими  записями  может  быть  выражена  в  виде системы

уравнений:

                   (n)           (n)            (n)

         интеграл R   (t - тау) U   (тау) dt = U   (t),     (1.17)

            0      1,2           2              1

    где спектр ядер R   (омега) имеет вид (1.15).  Решение каждого

                     1,2

Из уравнений системы (1.17) относительно q   (t) можно осуществить

                                          1,2

таким же образом, как  и решение  уравнения (1.12),  т.е. получить

его в виде:

                        f   (омега, альфа) U (омега) S (омега)

 альфа       1       -1  1,2                1         2

q     (t) = --- Re {F  {--------------------------------------}}, (1.18)

 1,2        эта                 U (омега) S (омега)

                                 2         1

где f   (омега, альфа) - регуляризирующая функция.  Учитывая,  что

     1,2

                          S (омега)

                           1                               2

     К(омега) = U (омега) ---------,  L(омега) = (К(омега))

                 2        S (омега)

                           2

                                 2р

    и  принимая  М(омега) = омега  ,  можно,   выбрав   конкретное

значение   параметра  р,   построить   регуляризирующую    функцию

f   (омега, альфа).

 1,2

    Отклонения    в    положениях    гипоцентров    землетрясений,

следовательно,  и  некоторые  различия  в  путях   распространения

сейсмических  колебаний,  эффекты  нелинейности  и  другие причины

                                                      (n)

могут  приводить  к тому,  что  отдельные  значения  q   (t) будут

                                                      1,2

несколько   отличаться  друг  от  друга.  Результирующее  значение

передаточной  функции будет определяться как среднее значение всех

реализаций:

                      _         1  N   (n)

                      q   (t) = - SUM q   (t).              (1.19)

                       1,2      N n=1  1,2

    При использовании записей инженерно-сейсмометрических  станций

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

таких случаях,  как  показано в [24], выполнять усреднение следует

                          (n)                             (m)

после  смещений  функций q   (t)  на отрезки  времени  тау   , при

                          1,2

которых наступают максимумы функции взаимной корреляции для каждой

               (n)

из реализаций q   (t) относительно одной из них.

               1,2

Описанный выше способ построения передаточных функций накладывает на них ряд ограничений. Функции могут быть определены только в классе абсолютно интегрируемых и будут описывать только воздействия деталей среды, сравниваемых по линейным размерам с длинами волн из частотного диапазона, в котором с достаточной точностью и детальностью удается получить спектры синхронных записей землетрясений из интересующей нас сейсмоактивной зоны. Точность описания воздействий различных деталей строения среды будет зависеть от того, насколько уровень полезного сигнала на соответствующих частотах будет превышать уровень помехи. Детально описанный выше алгоритм и пример его использования представлены в работе [21].

Близкая методика, основанная на тех же входных данных и позволяющая производить пересчет акселерограмм, полученных в одних грунтовых условиях, на площадки с другим близповерхностным разрезом, описана в работе [20].

    Задача  заключается  в  восстановлении сигнала,  подходящего к

подошве  слоистой  пачки,  и  получении акселерограмм на свободной

поверхности площадки. Для этого необходимо  знать спектр колебаний

на   свободной   поверхности  в  районе  станции  (опорная  точка)

F (омега) и спектральные характеристики  изучаемых  слоистых пачек

 v

Х (омега) в районе станции и Y (омега) на площадке. Входной сигнал

 v                            v

вычисляется по формуле:

                            N  F (омега)

        1     ДЕЛЬТА омега 2    v         i омега t

U (t) = -- Re    интеграл      --------- e          dомега. (1.20)

 0      пи          0          Х (омега)

                                v

    Соответствующая акселерограмма на площадке будет иметь вид:

                            N  F (омега)

         1    ДЕЛЬТА омега 2    v                  i омега t

U (t) = -- Re    интеграл      --------- Y(омега) e          dомега. (1.21)

 n      пи          0          Х (омега)

                                v

Спектральные характеристики среды рассчитываются по модифицированным формулам Томсона-Хаскелла для линейно-неупругой горизонтально-слоистой модели [25].

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

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

3. ПОДБОР АНАЛОГОВЫХ ЗАПИСЕЙ СИЛЬНЫХ ДВИЖЕНИЙ

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

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

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

На основании статистической обработки аналоговых спектров реакции вычисляются медианные и 84%-ной обеспеченности обобщенные спектры реакции грунта, которые соответствуют конкретным сейсмотектоническим и грунтовым условиям площадки размещения ОИАЭ. Кроме того, вычисляются наиболее вероятные максимальные ускорения и длительности, и оценивается степень их неопределенности.

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

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

Приложение 2

(справочное)

к Руководству "Определение

исходных сейсмических колебаний

грунта для проектных основ"

ОЦЕНКА И ГЕНЕРАЦИЯ РАСЧЕТНЫХ СЕЙСМИЧЕСКИХ

ВОЗДЕЙСТВИЙ МЕТОДАМИ, ОСНОВАННЫМИ НА МОДЕЛИ РАЗЛОМА

При оценке и генерации расчетных сейсмических воздействий (решение прямой задачи) необходимо:

1. Смоделировать процесс разрушения в разломе ("модель разлома").

2. Смоделировать излучение сейсмических волн в объем земной коры ("функции Грина").

3. Учесть локальные инженерно-геологические условия.

Задача может быть решена двумя методами: первый - чисто теоретический, где все явления описываются математически; второй - полуэмпирический, где отдельные части теории заменены экспериментальными данными.

1. ТЕОРЕТИЧЕСКОЕ МОДЕЛИРОВАНИЕ

При теоретических расчетах используются два подхода.

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

    Для расчетов короткопериодных движений грунта используются два

вида  стохастических  моделей  - барьерная модель А. Папагеоргиу и

                     2

К. Аки и модель омега  Т. Хэнкса и Р. Мак-Гуайра (рис. 2.1).

  ┌─ ─ ─ ─ ─ ─ ─ ─ ─ ┐                  ┌───────────────┐                                        ┌ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ── ┐

  │Детерминистическая│                  │ Теоретический │                                        │Параметры сейсмических воздействий

      модель очага            ┌────────>│    расчет     ├────────────────────────────────────────┼>│     ┌───────────────────┐      │

  │┌────────────────┐│        │         │"функций Грина"│                                        │ │     │             _2    │

   │Задание в источ-│         │         └───────────────┘                                          │     │Расчет а   , а ,   │      │

1 ││нике простран-  ││        │           /\           /\                                        │ │  ┌─>│        max        ├────┐

   │ственного и     ├┼────────┘           │            │                                           │  │  │V   , d   , T      │    │ │

  ││временного      ││              ┌─────┴─────┐ ┌────┴──────┐                                  │ │  │  │ max   max   max   │    │

   │распределения   │               │Локальные  │ │Передаточ- │                                    │  │  └───────────────────┘    │ │

  ││подвижек        ││              │грунтово-  │ │ные функции│                                  │ │  │  ┌───────────────────┐    │

   └────────────────┘               │геологичес-│ │среды      │                                    ├──┼─>│Расчет длительности├──┐ │ │

  └ ─ ─ ─ ─ ─ ─ ─ ─ ─┘              │кие условия│ │           │                                  │ │  ├─>│   акселерограммы  │  │ │

                                    └──────┬────┘ └────┬──────┘                                    │  │  └───────────────────┘  │ │ │

                                           │           │    ┌─ ─ ─ ─ ─ ─ ─ ─ ─ ─┐                │ │  │  ┌───────────────────┐  │ │

                                           │           │    │   Стохастический                     ├──┼─>│  Расчет спектров  │  │ │ │

                                           │           │           подход       │                │ │  │  │   реакции SA(t),  ├─>│ │

                                           │           │    │┌─────────────────┐   ┌────────────┐  │  ├─>│    SV(t), SD(t)   │  │ │ │

  ┌ ─ ─ ─ ─ ─ ─ ─ ─ ─┐ ┌────────────────┐  │           │  ┌─>│ Метод модельных ││  │Численная   ││ │  │  └───────────────────┘  │ │

     Стохастическая    │  5 параметров: │  │           │  │ ││  акселерограмм  ├┼┐ │модельная   │  │  │  ┌───────────────────┐  │ │ │

  │   модель очага   │ │1. длина;       │  │           │  │  │  (Монте-Карло)  │││ │акселеро-   ││ ├──┼─>│   Расчет спектра  │  │ │

   ┌────────────────┐  │2. ширина;      │  │           │  │ │└─────────────────┘ │ │грамма,     ├┼─┼─>│  │    мощности PS    │<─┘ │ │

  ││   Барьерная    ││ │3. максимальная │  \/          \/ │  ┌─────────────────┐││ │имитирующая ││ │  ├─>│                   │    │

   │  модель очага  ├┼>│   подвижка;    │ ┌──────────────┐│ ││Метод, основанный│ └>│реальную при│  │  │  └───────────────────┘    │ │

  ││А. Папагеоргиу и││ │4. скорость     ├>│  Амплитудный ├┼─>│    на теории    ││  │тех же      ││ │  │  ┌───────────────────┐    │

   │     К. Аки     │  │   вспарывания; │ │ спектр Фурье ││ ││    случайных    ├┼─>│средних     │  ├──┼─>│Расчет интенсив-   │<───┤ │

  │└────────────────┘│ │5. барьерный    │ └──────────────┘│  │    процессов    ││  │параметрах  ││ │  ├─>│ности I            │    │

2                      │   интервал     │         /\      │ │└─────────────────┘ ┌>│очага, среды│  │  │  └───────────────────┘    │ │

  │                  │ └────────────────┘         │       │  ┌─────────────────┐││ │и грунта    ││ │  │  ┌───────────────────┐    │

                                                  │       │ ││ Метод случайных │ │ └────────────┘  └──┼─>│Расчет коэффициента│    │ │

  │┌────────────────┐│ ┌────────────────┐         │       └─>│    вибраций     │││               │    │  │   динамичности    │<───┘

   │             2  │  │Модель очагового│         │         ││                 ├┼┘                    └─>│       бета'       │      │

  ││ Модель омега   ││ │ спектра Бруна  │         │          └─────────────────┘│                │       └───────────────────┘

   │очага Т. Хенкса ├┼>│М   ДЕЛЬТА сигма├─────────┘         └─ ─ ─ ─ ─ ─ ─ ─ ─ ─┘                └ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ─ ── ┘

  ││и Р. Мак-Гуайра ││ │ 0'             │

   └────────────────┘  └────────────────┘

  └ ─ ─ ─ ─ ─ ─ ─ ─ ─┘

Рис. 2.1. Блок-схема определения параметров сейсмических

колебаний грунта методами, основанными на "модели разлома"

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

    Барьерная  модель  описывается  пятью  основными  параметрами:

длиной   и   шириной   площадки   (очага), максимальной подвижкой,

скоростью вспарывания  и  барьерным   интервалом.  В  этой  модели

значение   высокочастотного  края  спектра F   зависит от размеров

                                            m

источника и других его параметров.

                  2

    В модели омега  значение f   определяется  эффектом  затухания

                              m

                                                          2

колебаний   в   приповерхностных   отложениях.    В  омега -модели

ускорения колебаний моделируются белым шумом в частотном интервале

f  - f  (где f  - угловая  частота  спектра  и  форма  спектра  по

 0    m       0

Дж. Бруну). Модель   характеризуется    величинами   момента   М ,

                                                                0

сброшенного напряжения ДЕЛЬТА сигма и f .  С  помощью  этой модели

                                       m

было    получено    хорошее    соответствие  экспериментальных   и

                              _2

теоретических значений а    и а    в интервале магнитуд 4,0 - 7,7.

                        max

Д. Бур  развил  обе  расчетные модели и  применил  метод случайных

колебаний, расширив диапазон магнитуд от 1,0 до 7,0. Использование

     2

омега -модели   для  расчетов   колебаний   от  очагов  с  М > 7,0

становится   некорректным   вследствие   нарушения  подобия,  т.е.

заметного  увеличения  коэффициента  формы K  =  L / W  в  области

                                            v

больших магнитуд.

    Для  расчетов  колебаний  от  очагов в виде барьерной модели и

     2

омега -модели применяются  два метода расчетов - метод Монте-Карло

и метод случайных колебаний.

    При использовании  обоих   методов   спектр   колебаний   R(f)

рассматривается в зависимости от следующих факторов:

               R(f) = С S(f) A(f) D(f) I(f),                 (2.1)

    где:

    С - масштабный коэффициент;

    S - спектр в источнике;

    D - затухание;

    I - реакция прибора в точке наблюдения;

                          R       F v

                           ТЭТА,j

                     С = --------------,                     (2.2)

                         4 пи ро  v   R

                                0  s0

    где:

    R       - функция направленности;

     ТЭТА,j

    F - эффект влияния свободной поверхности;

    v - часть энергии, приходящейся на горизонтальные  колебания;

    ро  и v   - плотность  и скорость  распространения  поперечных

      0    s0

волн в очаговой области;

    R - геометрическое расхождение волн.

                                            _

    Обычно принимается  F = 2  и  v = 1 / \/2.  На  расстояниях до

100 км для объемных волн  R = r (r - гипоцентральное  расстояние).

Для расстояний больше 100 км, где могут доминировать волны L ,

                                                            g

                                    _

                             R ~= \/r.

    Для  барьерной  модели  выражение  для  очагового спектра S(f)

[26]:

                           2        1/2

                        sin (пи f T)

 S(f) = {N [1 + (N - 1) ---------------]} S  (f) P(f, f   ), (2.3)

                                  2        0i          max

                          (пи f T)

    где:

    S  (f) -   спектр   индивидуального   субисточника   (круговая

     0i

трещина);

    N - общее число субисточников;

    Т - продолжительность вспарывания всего разлома,

                            ┌─              ─┐-1/2

                            │              2S│

               P(f, f   ) = │1 + (f / f   )  │    ,

                     max    │          max   │

                            └─              ─┘

                                   2

                         S  (f) = f  M  ,                    (2.4)

                          0i          0i

    где M   - сейсмический момент субисточника.

         0i

             2

    Для омега -модели:

                                          2

                 S(f) = M  / [1 + (f / f ) ],

                         0              0

                         6                       1/3

            f  = 4,9 х 10  v   (ДЕЛЬТА сигма / М )  ,        (2.5)

             0              s0                  0

   ДЕЛЬТА сигма   -  в  барах,  f   -  в  Гц, v    - в км/с и М  -

                                 0             s0              0

в дин. см.

    Учитывая нарушение принципа подобия для больших землетрясений,

необходимо   использовать   следующие   соотношения   для  спектра

источника [27]:

                                     1/2

           S(f) = M  / (1 + i f / f )   , f <= f

                   0               B            А

                        3/2                1/2

      S(f) = M  (f  / f)    / (1 + i f /f )   , f >= f ,     (2.6)

              0   a                      R            А

    где:

                   6      1/4                    1/3

      f  = 4,9 х 10  V   K    (ДЕЛЬТА сигма / М )   , <= М  ,

       A              sо                       0          0с

                       6      3/4                    1/3

          f  = 4,9 х 10  V   K    (ДЕЛЬТА сигма / М )   ;

           В              sо                       0

                6      1/4             1/3      -1/2

   f  = 4,9 х 10  V   K    ДЕЛЬТА сигма    М   М    , М  >= М  ,

    A              sо                       0с  0      0     0с

                     6      3/4                   1/3

        f  = 4,9 х 10  V   K    (ДЕЛЬТА сигма / М )  .       (2.7)

         В              sо                       0

    Здесь М   -  критическая  величина момента,  начиная с которой

           0с

нарушается  подобие.  Следует  отметить,  что в области сохранения

подобия для калифорнийских землетрясений K принималось равным 4.

Фактор усиления А(f) колебаний средой, т.е. так называемая специальная характеристика колебаний слоистой толщи, может быть рассчитана с применением метода Томсона-Хаскелла по программе Л.И. Ратниковой.

Иногда для определения А используется корректирующий фактор:

                         _________________

                       \/ро  v   / ро  v  ,

                           0  s0     i  si

    где   индекс  0  -  относится  к  области  очага,  i  -  точки

наблюдения.

    Фактор D(f) затухания с расстоянием записывается в виде:

               D(f) = exp[-пи f r / Q(f) V] P(f),            (2.8)

    где: Q - частотно  зависимый  фактор затухания в среде, разный

для разных регионов; P(f) - высокочастотный срез фильтра ;

                       P(f) = exp(-пи K  f).                 (2.9)

                                       0

    Наконец,  фильтр  I(f)  зависит  от  условий  регистрации. Для

стандартных акселерографов выражение для I(f):

                                   2

                                v f

               I(f) = --------------------------,           (2.10)

                        2    2

                      (f  - f ) - i (2 кси f f )

                             r                r

                                        2

    где f  = 25 Гц, кси = 0,6, V = (2 пи f ) .

         r                                r

    После   определения   спектра   R(f)   производится  следующая

процедура расчетов. Генератором случайных   чисел   (метод  Монте-

Карло)   генерируется  гауссовский  белый   шум  со спектром R(f),

фильтрованным  через  фильтр   коробчатой  формы  с  длительностью

                                        2

T   = 1 / f  + 0,05   для модели   омега   или 1 / f  + 0,05r  для

 W         о                                        А

модифицированной модели [28]. Здесь r - расстояние от очага.

    Далее производится обратное преобразование Фурье  во временную

область   и   получаются   временные   функции   a(t),  из которых

                               _2

оцениваются  величины  а    и  а .  Не  менее   20  генерированных

                        max

(расчетных) сигналов необходимы для хорошей оценки  а   .  Средний

                                                     max

по массиву реализаций спектр  лежит близко  к исходному расчетному

спектру.

    В   методе   случайных  вибраций  не  требуется  использование

генерации случайных  процессов. Спектр R(f) величины y (ускорение,

скорость,  смещение)  определяется  в  соответствии  с  принципом,

описанным выражением (2.1).

    Нулевой, второй и четвертый моменты  m , m , m  энергетической

                                          0   2   4

спектральной плотности рассчитываются по формуле:

                 1  омега        k       1/2

           m  = -- интеграл омега  / R(f)    dомега,        (2.11)

            k   пи    0

    где омега = 2 пи f.

                                _2

    Среднеквадратичная величина y   определяется соотношением:

                      _2             1/2

                      y   = (m  / T )   ,                   (2.12)

                              0    r

    где Т  (длительность колебаний) = Т     .

         r                             омега

    Для   определения   величины   спектра   реакции  используется

соотношение:

                                T          3

                                 0        v

               T  = T      + -------- (--------),           (2.13)

                r    омега   2 пи кси  v  + 1/3

                                        3

    где Т  и кси  -  период  собственных   колебаний  и  затухание

         0

осциллятора гамма = Т      / Т .

                     омега    0

    Математическое ожидание  Е(y   )  величины y    рассчитывается

                                max             max

с использованием точных либо асимптотических формул,  зависящих от

параметра эпсилон и числа экстремумов N, где:

                                          1/2

                    эпсилон = m  / (m  m )

                               2     0  4

                                      i+1  N

                           __     (-1)    C

                     _2   /пи  N           i        i

           E(y   ) = y  \/--- SUM ---------- эпсилон ,      (2.14)

              max          2  i=1      _

                                     \/l

         N

    где С  - биноминальный коэффициент, равный N! / l! (N - l).

         i

    Выражение   (2.14)    используется   для   малых    N  < 14,7,

                                                         i

эпсилон - 8.

    Для   больших   значений   N   предлагается   не   точное,   а

асимптотическое решение:

                  2           1/2                  1/2

       E(y   ) = y  {[2 ln(N)]    + гамма [2 ln(N)]   },    (2.15)

          max

    где в качестве гамма рекомендуется  константа  Эйлера,  равная

0,557216.

    После получения спектра  реакции с  помощью  метода  случайных

колебаний  задаются  стационарные  временные реализации с заданным

                            _2

спектром и величиной y    и y .

                      max

Указанные выше стохастические методы могут быть корректно использованы за пределами ближней зоны. Для ближней зоны была предложена модификация вышеописанного метода с использованием расчетного способа Монте-Карло. Однако в эту модель заложен ряд упрощающих предположений, в частности, предположение о постоянной по всему очагу скорости вспарывания. В то же время следует отметить большую чувствительность результатов расчетов к изменению скорости вспарывания [29].

Кроме того, применяется метод, основанный на представлении колебаний от большого очага как суммы колебаний от малых очагов - субисточников.

Сущность метода заключается в том, что эта субисточников совместно с временами начала их действия распределены случайно с однородной вероятностью в течение полного времени действия Т и излучаемые ими импульсы имеют масштабный коэффициент v.

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

гипоцентра.  На низких частотах излучение субисточников когерентно

с   уровнем   спектра,   пропорциональным   произведению   v  эта.

Высокочастотное   излучение   не  когерентно,  и  уровень  спектра

                    ___

пропорционален  v \/эта.  В  рамках  модели Хэнкса  и Мак-Гуайра и

представления  подобия  спектров  по  К. Аки  и  Дж. Брюну  спектр

смещений  на  высоких  частотах  изменяется  по  закону f      .

                                                         -альфа

Тогда:

            ┌   ┐2альфа/дельта       ┌   ┐1-2альфа/дельта

            │М  │                    │М  │

            │ 0 │                    │ 0 │

      эта = │---│             ;  v = │---│               ,  (2.16)

            │М  │                    │М  │

            │ 0е│                    │ 0е│

            └   ┘                    └   ┘

                                                                2

    где М  - момент очага; М   - момент субисточника. Для  омега -

         0                  0е

модели альфа = 2 и дельта = 3 [3].

    Колебания  рассчитываются  как "функции Грина" и суммируются в

точке  наблюдения.  Операции расчета  "функции  Грина"  описаны  в

фундаментальной  работе  К. Аки и П. Ричардса [26].  Авторы метода

постулировали  некоторые  случайные распределения субисточников по

поверхности   площадки   разрыва  в  соответствии  с  предлагаемой

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

связано  с  предполагаемой  неоднородностью очага, и, естественно,

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

реальному   распределению   неоднородностей.   Так,   при  задании

распределения субисточников можно принимать во внимание результаты

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

подвижки   принимать  пропорциональной  длине  субразрыва.  Однако

следует иметь в виду, что в каждом конкретном случае распределение

параметров может отличаться от среднего или характерного по многим

случаям.   Так,   результаты   теоретических  и  экспериментальных

исследований  позволяют  сделать  вывод  о  тяготении максимальной

подвижки  D     к центру источника и о спаде значений D до нуля на

           max

краях  разрывов.  Однако в отдельных конкретных случаях отклонение

от  этой  тенденции  может  быть  довольно  значительным,  так как

определяется   пространственными  различиями  прочностных  свойств

горных пород.

Следует отметить перспективность разработанной технологии расчетов колебаний. Однако в настоящее время из-за слабого знания строения очагов реальных землетрясений описанные выше методы теоретических расчетов не находят широкого применения.

2. ПОЛУЭМПИРИЧЕСКИЙ МЕТОД

Полуэмпирический метод был предложен Харцелом. В этом методе в качестве "функций Грина" используются записи слабых землетрясений, что позволяет избежать необходимости вычисления этих функций теоретическими методами. Полуэмпирический метод является наиболее практичным при определении параметров сейсмических воздействий на основе модели разлома, но для его применения необходимы записи хотя бы слабых землетрясений. Ключевым вопросом применимости полуэмпирического метода является наличие в имеющихся записях составляющих, соответствующих периодам собственных колебаний проектируемых сооружений.

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

Полуэмпирический метод применялся и применяется в настоящее время в различных вариантах. Например, Канамори синтезировал волны Лява с периодами от 2 до 10 с с учетом разницы в сейсмических моментах между слабыми и сильными землетрясениями.

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

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

В России полуэмпирический метод в последнее время применен на Камчатке. Здесь ограниченность числа записей сильных движений компенсировалась использованием теоретических моделей, причем отдельные параметры этих моделей определены по слабым местным землетрясениям (рис. 2.2). С этой целью разработан простой приближенный метод решения прямой задачи прогноза параметров колебаний, в основу которого положен метод случайных функций; заданы параметры модели по имеющимся на тот момент и вновь полученным данным.

                               ┌──────────────────────────┐

               ┌─────────────┐ │        ┌───────────────┐ │

               │   Заданные  ├─┼─────┬─>│Расчет K  -    │ │

               │    М , R    │ │     │  │        Q      ├┐│

               │     W       │ │┌────┼─>│учет поглощения│││

               └─────────────┘ ││    │  └───────────────┘││

┌──────────┐   ┌─────────────┐ ││    │  ┌───────────────┐││   ┌──────────────┐    ┌───────────┐

│Опорный   │   │Опорный      │ ││    │  │Расчет K  -    ││└──>│Расчет спектра├─┬─>│   Счет    │

│спектр    │   │спектр Фурье │ ││    │  │        R      │└───>│ускорений FSA │ │  │ А  , V    │

│реакции   ├──>│ускорений для│ ││    ├─>│учет расхож-   ├────>│для заданных  │┌┼─>│  ск   ск  │

│для М  ,  │   │фиксированных├─┘│    │  │дения для      │     │M , R         │││  │           │

│     WO   │   │М  = М  ,    │  │ ┌──┼─>│протяженного   │┌───>│ W            │││  └───────────┘

│R         │┌─>│ W    WO     │  │ │  │  │источника      ││┌──>│              │││

│ 0        ││  │R = R        │  │ │  │  └───────────────┘││   └──────────────┘││  ┌───────────┐

└──────────┘│  │     0       │  │ │  │                   ││                   ││  │   Счет    │

            │  └─────────────┘  │ │  │  ┌───────────────┐││                   ├┼─>│А   , V    │

            │                   │ │  ├─>│Расчет K  -    │││                   ││  │ max   max │

┌──────────┐│   ┌─────────────┐ │ │  │  │        M      ├┘│                   ││  └───────────┘

│Расчет    ││   │ Параметры   │ │ │┌─┼─>│учет магнитуды │ │                   ││

│длитель-  ││   │ среды       ├─┘ ││ │  └───────────────┘ │                   ││  ┌───────────┐

│ности для ├┘   │эпсилон , Q ,│   ││ │                    │                   │├─>│   Счет    │

│М  , R    │    │       S   0 ├──┐││ │  ┌───────────────┐ │                   ├┼─>│  спектра  │

│ WO   0   │    │гамма , Т    │  │││ │  │Расчет K  -    │ │             D ────┼┼─>│ реакции RA│

└──────────┘    │     Q   100 │  │││ │  │        g      ├─┘                   ││  └───────────┘

                └─────────────┘  │││ │  │учет грунтовых │                     ││

                                 │││ │  │условий        │<───── g             ││  ┌───────────┐

               ┌────────────────┐│││ │  └───────────────┘                     │└─>│   Счет    │

               │Параметры закона││││ │                                        │   │  спектра  │

               │масштабирования ││││ │  ┌───────────────┐                     ├──>│мощности PS│

               │очагов:         ├┼┘│ ├─>│Расчет очаговой├─┐                   │   └───────────┘

               │связь L (M )    ││ │┌┼─>│длительности   │ │                   │

               │       S  W     ├┼─┘││  └───────────────┘ │  ┌───────────────┐│   ┌───────────┐

               │связь FSA(M )   ││  ││                    └─>│     Расчет    ││   │Счет       │

               │           W    ├┼──┘│  ┌───────────────┐    │  длительности ├┤   │интенсив-  │

               │связь T(M )     ││   │  │Расчет компо-  │ ┌─>│акселерограммы │└──>│ности I    │

               │         W      ││   └─>│ненты длитель- │ │  └───────────────┘    └───────────┘

               │параметр        ││      │ности, опреде- ├─┘

               │ДЕЛЬТА lg L     │└─────>│ляемой средой  │

               │           S    │       │               │

               └────────────────┘       └───────────────┘

Рис. 2.2. Блок-схема алгоритма

моделирования параметров колебаний грунта

    Первым этапом алгоритма   для   прогноза   параметров  сильных

сейсмических  колебаний грунта является расчет спектра Фурье (рис.

2.2).  Исходным  в  основном  варианте  алгоритма  очаговый спектр

М (f). Спектр Фурье FSA(f) рассчитывается по формуле:

 0

                         2

FSA(f / R, M ) = (2 пи f)  С х М (f / М ) х K  х K (f) х K (f) х K (f), (2.17)

            W                   0      W     R    Q       g       H

    где:

    С - масштабный коэффициент;

    M (f / M ) - очаговый спектр (в частности,  можно использовать

     0      W

модель   Бруна);

    K  -  член,     учитывающий     геометрическое     расхождение

     R

(сферическое)  с  поправкой  на  размер  источника  (зависящий  от

магнитуды);

    К (f) -  член,   учитывающий   потери   за  счет  поглощения и

     Q

рассеяния;

    K (f) -  поправка  к  спектру,  учитывающая  категорию  грунта

     g

(по СНиП: I, II или III;  соответственно:  скальный,  средний  или

мягкий грунт);

    К (f) - низкочастотный  фильтр,  определяющий срез  спектра на

     H

высоких частотах (параметр f   ).

                            max

    Вторым важным элементом расчета является  оценка  длительности

акселерограммы.  В  качестве  характеристики длительности в методе

случайных     колебаний     удобно    использовать    "эффективную

(прямоугольную)"  длительность Т  - длительность эквивалентного по

                                П

энергии  сигнала  с прямоугольной огибающей и среднеквадратической

амплитудой,   равной   среднеквадратической   амплитуде  исходного

сигнала  в  области  максимальных  колебаний.  Для  оценки  полной

длительности  в  системе  "очаг-среда-грунт"  удобно  использовать

"среднеквадратичную"  длительность  Т  ,  которая определяется как

                                     ск

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

комбинацией вкладов очага, среды и грунта:

                    2     2       2       2

                   T   = T     + T     + T    .             (2.18)

                    ск    ск,s    ск,m    ск,g

    Переход от Т   к Т   осуществляется по формуле:

                ск    П

                          Т  = К  х Т  .

                           П    Т    ск

    Коэффициент   К    может   принимать   значения   от   2  (для

                   Т

акселерограммы с гауссовой огибающей) до 3,5 (для акселерограммы с

прямоугольной огибающей).

    Расчет полной "среднеквадратичной" длительности осуществляется

по формуле (2.18) следующим образом:

    1) считается, что сигнал в очаге имеет прямоугольную огибающую

длительностью  T , тогда: T     = T  / 3,5. Длительность сигнала в

                S          ск,s    S

очаге   T   пропорциональна   размеру   очага:   T  = L (M ) / v .

         S                                        S    S  W     s

Коэффициент v  имеет смысл скорости разрыва в очаге;

             s

    2) "среднеквадратичная" длительность  импульсной реакции среды

оценивается по эмпирической формуле:

                    Т     = T    х (R / 100).               (2.19)

                     ск,m    100

    Параметр   T      (длительность   при   R = 100 км)   является

                100

настраиваемым входным параметром алгоритма.

Запись слабого землетрясения, полученную на скальном грунте, можно рассматривать как "функцию Грина" среды. Поэтому изучение параметров слабых землетрясений позволяет получать параметры среды для использования в расчетах параметров сильных движений.

    Зависимость  длительности  слабых  землетрясений от расстояния

можно использовать при расчете полной  длительности акселерограммы

по формулам  (2.18)  и  (2.19).   Осредненные   огибающие   слабых

землетрясений  используются  в  качестве огибающих "функции Грина"

среды  по  мощности  при  построении  модельных  акселерограмм. По

модельным        акселерограммам        далее       рассчитываются

"среднеквадратичная"   и  "прямоугольная"  длительности,  а  также

длительности   Т      (по   Аптикаеву)  и  Т    (по  Трифунаку)  и

                0,5                         90

определяются   коэффициенты   перехода   от   "среднеквадратичной"

длительности  к  длительностям  в  других  определениях. Огибающие

слабых  землетрясений  имеют  также  и самостоятельный интерес как

форма импульса рассеянных волн в реальной среде.

    Согласно определению "среднеквадратичной"   длительности,  при

расчете  Т    требуется  интегрирование  "хвоста"  записи (коды) с

          ск

        2

весом  t .  Это  приводит к существенной зависимости результата от

случайных  выбросов в коде, что нежелательно. Чтобы устранить этот

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

простой   метод,  предложенный  Раутиан  Т.Г.  (1981),  в  котором

считается,   что  кода   начинается   при  t = t  + k х (t  - t ).

                                                s         s    p

Значение  коэффициента  k = 2  было определено  по насыщению роста

значений Т     с ростом k.

          ск,m

    Существенный вклад в запись и  тем  самым  в результат расчета

"среднеквадратичной"   длительности  соответственно  вносит  также

микросейсмический  шум.  Характер искажений - фиктивное увеличение

Т    .  Зависимость  параметра  Т      от R можно аппроксимировать

 ск,m                            ск,m

линейной моделью вида: lg Т     = а + b lg R,  отдельно для каждой

                           ск,m

компоненты и каждого частотного интервала.

    При  построении  средних  огибающих допускается предположение,

что  форма  огибающей  стабильна,  то  есть  нормированная функция

огибающей  имеет  вид  F(t / T    )   и F(x)  одинакова  для  всех

                              ск,m

расстояний.   С   расстоянием   меняется   только масштаб времени,

задаваемый длительностью  Т    .   Такое  представление  позволяет

                           ск,m

определить огибающую  по  выборке  небольшого  объема как  среднее

нормированных огибающих, приведенных к единому расстоянию.

Приложение 3

(справочное)

к Руководству "Определение

исходных сейсмических колебаний

грунта для проектных основ"

ОЦЕНКА И ГЕНЕРАЦИЯ РАСЧЕТНЫХ СЕЙСМИЧЕСКИХ ВОЗДЕЙСТВИЙ

МЕТОДОМ, ИСПОЛЬЗУЮЩИМ СТАНДАРТНЫЕ СПЕКТРЫ

    Единичная  акселерограмма с  точно  заданным  спектром реакции

может генерироваться с учетом [17]. Есть метод генерации расчетной

акселерограммы, спектр реакции которой SA(T, кси = 0,05)  по форме

соответствует кривой динамического коэффициента в действующем СНиП

11-7-81 [4]. Акселерограммы могут быть записаны в виде:

             0               na

            а    (t) = А(t) SUM B  sin(омега t + фи),        (3.1)

             бета           i=1  i

    где:

    A(t) -  огибающая,  а  фазовые  углы  фи   представляют  собой

                                            i

равномерно  распределенные  в  интервале  от  0 до 2 пи  случайные

величины.  Шаг  по  частоте  ДЕЛЬТА  омега определяется из условия

гладкости   спектра  реакции  как (омега  - омега   ) / омега    =

                                        i        i-1         i-1

= 0,06347. Расчеты проводятся в интервале периодов от 0,05 до 3 с.

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

                                       i,1

непосредственно    взятые    с   заданной   кривой   бета(Т)   для

соответствующих значений частот омега . Следующее приближение:

                                     i

                                   бета(Т)

                       B    = B    --------,                 (3.2)

                        i,2    i,1 бета (Т)

                                       1

    где:

    бета(Т) - кривые  СНиП,  а бета (T) - кривая,  полученная  при

                                   1

первом шаге итерации.

    Нужная точность  совпадения  расчетного   спектра  с  заданным

получается  после  ряда итераций. Применение программы на практике

показало,  что отличия расчетного спектра от заданного на периодах

Т > 0,05 с   могут   достигать   15%    после    10 итераций.    В

короткопериодной области вследствие наличия систематических ошибок

уровень расчетной кривой бета (Т) на периодах Т -> 0 не доходит до

                             i

значения бета(Т = 0) = 1.  Однако  близость  расчетных  и заданных

кривых  в  инженерном  интервале  периодов  допускает практическую

применимость метода.

                                              0

    Для  перехода  от базовой акселерограммы а    (Т)  к расчетной

                                              бета

          0

значения а    (t) следует умножить на коэффициент K:

          бета

                     ________

                     A  К  К

                      p  1  2            0

                K = ----------, a (t) = a    (t) К.          (3.3)

                        m        p       бета

                         кp

    Согласно  СниПа,  К  принимается  равным  0,1;  0,2;  0,4  для

расчетной сейсмичности 7, 8, 9 баллов соответственно;

    0,12 <= К  < 1 учитывает допускаемые повреждения сооружений;

             1

    0,85 <= m   <= 1,15    учитывает    повторяемость    расчетных

             кp

землетрясений.

Выбор расчетных спектральных кривых, уровня ускорений и длительности колебаний следует проводить с учетом необходимой обеспеченности (степени риска).

Приложение 4

(справочное)

к Руководству "Определение

исходных сейсмических колебаний

грунта для проектных основ"

ОПРЕДЕЛЕНИЕ БАЛЛЬНОСТИ ПЛОЩАДКИ

ДЛЯ ЕСТЕСТВЕННЫХ И ТЕХНОГЕННО-ИЗМЕНЕННЫХ УСЛОВИЙ

1. Балльность определяется макросейсмической интенсивностью I в баллах шкалы MSK-64.

2. При оценке I используются региональные соотношения между I, магнитудой М и расстоянием до очага r. Оценки I по эмпирическим соотношениям относятся к базовому (эталонному) грунту. Интенсивность I оценивается по формуле:

                    I = 1,5M - гамма lg r + c,               (4.1)

где гамма, с - коэффициенты,  характеризующие затухание балльности

в   зависимости   от   направления  сейсмического  излучения.  Для

направления  вкрест  структур гамма = -3,5 - 4,2, с = 3,5 - 4,4 (в

зависимости  от  региональных особенностей); для направления вдоль

структур  гамма  =  -3,0  -  3,4;  с  =  2,5  -  3,3.  Для средних

рекомендуется принимать гамма = -3,5, с = 3,0.

    Для М = 4,75 - 7,5 и r <= 150 км можно использовать формулу:

               I = (0,0075М - 0,08) r + 1,35М.               (4.2)

    3.  Изменение   интенсивности  ДЕЛЬТА I     за   счет   гидро-

                                           i

геологических условий оценивается по формуле:

                                                  -0,04h

        ДЕЛЬТА I  = 1,67 lg(ро  V  / ро  V ) + Y e      ,    (4.3)

                i             0  0     i  i

    где:

    ро - плотность грунта;

    V - скорость распространения волны сдвига;

    pо  V   и  pо  V  -   сейсмические  жесткости   эталонного (0)

      0  0       i  i

и исследуемого (i) грунта;

    h - глубина залегания уровня грунтовых вод;

    Y - коэффициент, учитывающий его влияние.

Приложение 5

(справочное)

к Руководству "Определение

исходных сейсмических колебаний

грунта для проектных основ"

ЭКСПРЕСС-МЕТОД ОЦЕНКИ КОЛИЧЕСТВЕННЫХ ХАРАКТЕРИСТИК

СЕЙСМИЧЕСКИХ КОЛЕБАНИЙ ГРУНТА ЗАДАННОЙ

ВЕРОЯТНОСТИ ПРЕВЫШЕНИЯ

1. Исходные данные для проведения экспресс оценок:

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

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

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

1.4. Выборки максимальных ускорений колебаний грунта при землетрясениях заданной интенсивности.

1.5. Динамические и спектральные особенности акселерограмм сильных землетрясений и их спектров реакции в зависимости от параметров очага землетрясения (магнитуда М), эпицентрального расстояния ДЕЛЬТА и грунтовых условий в пункте регистрации.

2. Расчет сейсмических колебаний грунта

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

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

                                               _____

    2.1.  Нормативное сейсмическое воздействие S (T),  учитываемое

                                                a

при   расчете   сейсмостойкости    зданий   и    сооружений     по

спектрально-динамической теории, определяется по формуле:

                      _____   __   _______

                      S (T) = a  х бета(Т),                  (5.1)

                       a       J

    где:

    Т  -  период  собственных колебаний зданий и  сооружений;

    __

    а  -  нормативное  ускорение;

     J

    _______

    бета(Т) - коэффициент динамичности.

                                                             _____

    2.2. Уровень  нормативного  сейсмического   воздействия  S (T)

                                                              a

определяется  нормативным  ускорением колебаний грунта  0,1; 0,2 и

0,4 g при  расчетной сейсмичности площадки соответственно 7, 8 и 9

баллов.   Влияние   повторяемости  сотрясений  в  нормах  частично

учитывается коэффициентом m  .

                           кp

    2.3.  Форма  нормативного  сейсмического воздействия  задается

кривой коэффициента динамичности бета(Т),  которая  не зависит  от

интенсивности   землетрясения   и   полностью  определяется  типом

грунтовых условий.

    Установленные   с   учетом   соотношения   (5.1)   нормативные

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

отличаться от спектров реакции реальных акселерограмм. Это отличие

обусловлено: разницей реальных и нормативных уровней ускорений при

землетрясении заданной интенсивности; отличием ординат приведенных

спектров  реакции  реальных  землетрясений  от  нормативной кривой

коэффициента  динамичности бета(Т).   Вероятность  Р__  превышения

                                                    aJ

                        __

нормативного  ускорения a  реальным a   при землетрясении заданной

                         J           PJ

интенсивности J ~= 0,35;  вероятность Р_______  превышения  уровня

                                       бета(Т)

кривой  коэффициента динамичности   бета(Т)  уровнем  приведенного

спектра  реакции  |S (T)|  реального  землетрясения примерно равна

                    a

0,5.

    2.4.  Полная  вероятность  P __     превышения    нормативного

                                SaJ

сейсмического    воздействия   при    заданной    вероятности   P

                                                                 J

землетрясения интенсивностью J в  первом  приближении  может  быть

оценена по формуле:

                    Р __ = Р  х Р__ х Р_______.              (5.2)

                     SaJ    J    aJ    бета(Т)

    2.5. Оценить годовую вероятность P   можно по среднему периоду

                                      J

повторения  сотрясений заданной интенсивности Т , предполагая, что

                                               J

характер распределения землетрясений во времени подчиняется закону

Пуассона:

                     P  = 1 - exp(-1 / Т ).                  (5.3)

                      J                 J

    При  достаточно  больших значениях Т  вместо соотношения (5.3)

                                        J

можно использовать более простое соотношение:

                            Р  = 1 / Т .

                             J        J

    В случае  необходимости  учета срока службы объекта t  формулу

                                                         0

(5.3) можно представить в виде:

                     Р  = 1 - exp(-t  / Т ).

                      J             0    J

    2.6.   Вероятности    превышения    нормативных   сейсмических

воздействий  для  зон   возникновения   землетрясений  с  периодом

повторения сотрясений 1 раз в 100, 1000 и 10000 лет соответственно

равны:

            100             -3    1000             -4

         Р ----- = 1,75 х 10  ; Р ----- = 1,75 х 10  ;

           S (Т)                  S (Т)

            a                      a

                         10000            -5

                       Р ----- = 1,75 х 10  .

                         S (Т)

                          a

    Вычисленный   уровень   вероятности   превышения   нормативных

сейсмических   воздействий   ниже  уровня  вероятности  проектного

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

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

вероятности событий Р     .

                     к.у.в

    2.7.    Для    определения   уровня   расчетных   сейсмических

воздействий,     вероятность     превышения     которых      будет

            -6    -1

Р      <= 10   год  , учитываются:

 к.у.в

    - сильнейшие землетрясения более редкой повторяемости;

    - рост уровня расчетных ускорений по сравнению с нормативными;

    - обобщенные кривые приведенных спектров  реакции, превышающих

                           _______

уровень нормативных кривых бета(Т).

    2.8.  Между  уровнями  максимальных ускорений колебаний грунта

a  , вероятностью превышения этого  ускорения  Р__  и  нормативным

 PJ                                             aJ

                     __

расчетным ускорением a  при заданной интенсивности землетрясения J

                      J

существует корреляционная зависимость [30] в виде:

                 __

           a   = a  [0,54 - 0,63 log(-log(1 - Р__))],        (5.4)

            PJ    J                            aJ

    где выражение в квадратных скобках  - коэффициент безопасности

К   ,  обеспечивающий   переход   от  нормативного   сейсмического

 без

            __

воздействия a   к воздействию  a   заданной вероятности превышения

             J                  PJ

Р__.

 aJ

    2.9. Значения коэффициента  безопасности К   ,  обеспечивающие

                                              без

заданный уровень вероятности превышения сейсмического воздействия,

даны   в   табл.   5.1,    где   приведены   значения   К      для

                                                         без

                  -6

Р _ ~= Р      = 10    при  Р        ~= 0,5  для различных периодов

 Sa     к.у.в               бета(Т)

повторения сотрясений Т   и сроков службы объекта t .

                       J                           0

Таблица 5.1

ЗНАЧЕНИЯ КОЭФФИЦИЕНТОВ БЕЗОПАСНОСТИ

┌────────┬───────────────────────────────────────────────────────┐

│        │                             -6    -1                  │

│   t ,  │Значения К    для Р      = 10   год   при бета  ~= 0,5 │

│    0   │          без      к.у.в                      m        │

│  годы  ├───────────────────────────────────────────────────────┤

│        │                         Т                             │

│        │                          J                            │

│        ├──────────┬───────────┬────────────────────────────────┤

│        │   100    │   1000    │       10000 (4545) <*>         │

├────────┼──────────┼───────────┼────────────────────────────────┤

│1       │   3,1    │    2,5    │              1,8               │

│        │  (1,8)   │   (1,2)   │             (0,6)              │

├────────┼──────────┼───────────┼────────────────────────────────┤

│30      │   4      │    3,4    │              2,8               │

├────────┼──────────┼───────────┼────────────────────────────────┤

│100     │   4,2    │    3,7    │              3,1               │

├────────┼──────────┼───────────┼────────────────────────────────┤

│1000    │   4,4    │    4,2    │              3,7               │

└────────┴──────────┴───────────┴────────────────────────────────┘

    ----------------------------

     <*>  В  скобках  приведены  значения  К   ,  вычисленные для

                                            без

                                   -4

 вероятности   превышения Р _ ~= 10    для   периода   повторения

                           Sa

 сотрясений 1 раз в 100, 1000 и 4545 лет.

2.10. Оценка формы спектров реакции грунта для конкретных сейсмотектонических и инженерно-геологических условий производится с учетом:

- спектров реакции аналоговых акселерограмм, рекомендованных для различных зон сейсмичности территории [31];

- спектров реакции калифорнийских и других землетрясений.

    Уровень максимальной реакции линейных осцилляторов  в  области

низких  частот  (менее 4 Гц) - огибающая S   (Т, дельта)  спектров

                                          ога

реакции  реальных  землетрясений,  нормированных  на  максимальное

ускорение грунта а    (рис. 5.1 - не приводится) [32]:

                  max

                                                __

             S   (Т , дельта) = 2,55К       / \/Т ,          (5.5)

              ога  0                 дельта      0

    где К       = 8 / (дельта + 3,3) + 0,4.

         дельта

2.11. Выборочные оценки средних коэффициентов динамичности и стандартных отклонений, выполненные для различных значений собственных периодов линейных осцилляторов в диапазоне периодов 0,2 - 4 с, позволили установить:

- уменьшение среднего значения коэффициента динамичности с ростом собственного периода линейного осциллятора;

    - зависимость стандартного отклонения  сигма  от  собственного

периода  линейного осциллятора Т  и среднего значения коэффициента

                                0

              __

динамичности |Sа| для затухания 5%-ного критического в виде:

                                          __

                  сигма = (0,34 + 0,3Т ) |Sa|.               (5.6)

                                      0

2.12. Значения средних коэффициентов динамичности для диапазона периодов 0,2 - 4 с приведены на рис. 5.1 (кривая 6) и удовлетворительно согласуются с нормативными кривыми коэффициента динамичности бета (Т) для грунтов I, II и III категорий по сейсмическим свойствам (кривые 2, 3 и 4 на рис. 5.1 соответственно). На этом же рисунке приведена стандартная кривая нормированного на максимальное ускорение грунта спектра реакции (кривая 5), рекомендованного МАГАТЭ для горизонтальной составляющей при затухании, 5% [2].

    2.13.  Кривую обобщенного спектра реакции грунта в зависимости

от характерного периода Т    можно записать в виде [33]:

                         max

                                       Т

                    __                  i  2

                   |Sa| = exp[-(d log ----) ],               (5.7)

                                      T

                                       max

    где:

    Т  - период колебаний грунта;

     i

    d - показатель затухания.

    2.14. Влияние скорости распространения поперечных сейсмических

волн v , м/с, на показатель затухания оценивается по формуле:

      s

                     d ~= -0,2 + 0,8 log v .                 (5.8)

                                          s

    2.15. Для оценки формы обобщенного  спектра реакции грунта для

конкретных условий необходимо иметь сведения о характерном периоде

спектра   реакции   Т     и  скорости  распространения  поперечных

                     mах

сейсмических волн на площадке.

    2.16. Характерный  период спектра реакции [33] Т    получается

                                                    mах

по эмпирической зависимости:

                                                          3 ____

log Т   (М, ДЕЛЬТА, Т   ) = 0,08М + 0,16 log ДЕЛЬТА + log \/Т    - 1, (5.9)

     mах             рез                                     рез

    где: Т    -  резонансный  период  колебаний  грунта  площадки,

          рез

определяется экспериментально или расчетом:

                        Т    = 4h / v ,                     (5.10)

                         рез         s

    где:

    h - мощность слоя, (м);

    v  - скорость распространения поперечных  сейсмических  волн,

     s

м/с.

    2.17.    Уровень    прогнозируемого   спектра   реакции   S

                                                               пра

определяется    значением    коэффициента    безопасности    К   ,

                                                              без

                                 __

максимальным  ускорением грунта  а ,  уровнем  огибающей  спектров

                                  J

S   (Т, дельта),  устанавливаемым для характерного периода спектра

 ога

реакции Т    (при Т  = Т   )  по  формуле  (5.5),  и коэффициентом

         max       i    max

ослабления  К   ,  учитывающим  разницу  в  уровне спектров  S (Т)

             осл                                              а

реакции   жестких  рыхлых   грунтов   при  одинаковых  характерных

периодах:

                       __                         __

           S    = K    а  S   (Т          ) K    |S |,      (5.11)

            пра    без  J  ога  max,дельта   осл   а

    где К    = -1,7 + 0,9 log v  при 100 м/с < v  < 1000 м/с.

         осл                   s                s

    При v  >= 1000 м/с К    = 1, а при v  <= 100 м/с = 0,1.

         s              осл             s

                                                __

    2.18. Среднемаксимальные ускорения  грунта  а , см/кв. с,  при

                                                 J

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

формуле [31]:

              __          J

              а  = 0,8 х 2  при 0,1 с <= Т <= 0,5 с;        (5.12)

               J

              __   0,4   J

              а  = ---- 2  при 0,5 с < Т < 1,5 с.           (5.13)

               J    Т

    2.19. Среднемаксимальное ускорение   практически  совпадает  с

нормативным   расчетным   ускорением  при  землетрясении  заданной

интенсивности J при 0,1 с <= Т  <= 0,5 с.

                              i

2.20. Интенсивность землетрясения J при известных параметрах очага (магнитуда М, глубина очага Н) и эпицентральном расстоянии ДЕЛЬТА, коэффициентах уравнения микросейсмического поля а, b, с и грунтовых условиях на площадке можно оценить по формуле макросейсмического поля [19].

Приложение 6

(справочное)

к Руководству "Определение

исходных сейсмических колебаний

грунта для проектных основ"

СПИСОК ИСТОЧНИКОВ ИНФОРМАЦИИ

1. ПНАЭ Г-05-035-94 "Учет внешних воздействий природного и техногенного происхождения на ядерно- и радиационно-опасные объекты".

2. 50-SG-S1 Руководство по безопасности. "Учет землетрясений и связанных с ними явлений при выборе площадок для атомных станций". МАГАТЭ, 1981. 70 с.

3. Михайлова Н.Н. Прогноз сейсмических воздействий при учете местных условий. Автореферат диссертации канд. физ.-мат. наук. Москва: 1983. 24 с.

4. СНиП 11-7-81. Строительство в сейсмических районах. М.: Стройиздат, 1982. 49 с.

5. Сейсмическое районирование территории СССР. М.: Наука, 1980. 307 с.

6. Аптикаев Ф.Ф., Михайлова Н.Н. Форма спектра реакции в ускорениях для отдельного землетрясения. В кн.: Макросейсмические и инструментальные исследования сильных землетрясений. М.: Наука, 1985. С. 142 - 144.

7. Аптикаев Ф.Ф., Михайлова Н.Н. Некоторые корреляционные соотношения между параметрами сильных движений грунта // Бюллетень по инженерной сейсмологии. N 12. Ереван: издательство АН АрмССР, 1988. С. 48 - 56.

8. Аптикаев Ф.Ф., Шебалин Н.В. Уточнение корреляций между уровнем макросейсмического эффекта и динамическими параметрами движения грунта. В кн.: Исследования по сейсмической опасности. М.: Наука, 1988. С. 98 - 108.

9. Айзенберг Я.М., Залилов К.Ю. Методика генерирования расчетных ансамблей синтетических акселерограмм на основе региональной сейсмологической информации. В кн.: Детальные инженерно-сейсмологические исследования. М.: Наука, 1986. С. 61 - 74.

10. Аптикаев Ф.Ф., Роджан К., Фролова Н.И. Форма огибающей амплитуд ускорений на записях сильных движений. В кн.: Советско-американские работы по прогнозу землетрясений. Душанбе, М.: Дониш, 1979. Т. 2, кн. 2. С. 139 - 147.

11. Быков В.В. Цифровое моделирование в статистической радиотехнике. М.: Советское радио, 1971. 328 с.

12. Trifunac M.D. Lec. V. Routine computer prosessing of strong motion accelerograms. Earthquake Engeneering Research Laboratory, Report N EERL 73-03 Calif. Inst. of Tech, Pasadena, 1973. 360 p.

13. Vanmarke E., Lai Sh. Strong motion duration and RMS amplitude of earthquake records // Bull. Seismol. Soc. Amer. 1980. Vol. 70, N 4. P. 1293 - 1307.

14. Аптикаев Ф.Ф. Учет длительности колебаний при инструментальной оценке сейсмической интенсивности // Сейсмическая шкала и методы измерения сейсмической интенсивности. М.: Наука, 1975. С. 234 - 239.

15. Joyner W.B., Boore D.M. Measurement, characterization and prediction of strong ground motion // Proc. Earth. and Soil. Dyn. Div. ASCE. 1988. 11GT. P. 43 - 102.

16. Харкевич А. А. Спектры и анализ. М: Физматгиз, 1962. 234 с.

17. Ньюмарк Н., Розенблют Э. Основы сейсмоактивного строительства. М.: Стройиздат, 1980. С. 61 - 99.

18. Корчинский И.Л., Поляков С.В., Быховский В.А. и др. Основы проектирования зданий в сейсмических районах. М.: Стройиздат, 1961. С. 67 - 111.

19. Шебалин Н.В. Методы использования инженерно- сейсмологических данных при сейсмическом районировании. В кн.: Сейсмическое районирование СССР. М.: Наука, 1968. С. 95 - 111.

20. Ратникова Л.И., Сакс М.В., Кронрод Т.Л. К вопросу о пересчете акселерограмм сильных движений на различные грунтовые условия // Исследования по сейсмической опасности. М.: Наука, 1988. С. 81 - 93. (Вопр. инж. сейсмологии; Вып. 29).

21. Кендзера А.В., Скляр А.М., Роман А.А. и др. О возможности использования эмпирических передаточных функций среды при микрорайонировании территорий со сложным геологическим строением // Оценка эффекта сильных землетрясений. М.: Наука, 1989. С. 82 - 89. (Вопр. инж. сейсмологии; Вып. 30).

22. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука, 1979. 285 с.

23. Верлань А.Ф., Сизиков В.С. Интегральные уравнения: Методы, алгоритмы, программы. Киев: Наук. думка, 1986. 543 с.

24. Кендзера А.В., Стародуб Г.Р. Получение эмпирических частотных характеристик среды для прогнозирования воздействий глубокофокусных карпатских землетрясений // Современные геодинамические процессы и прогноз землетрясений. Киев: Наук. думка, 1987. С. 62 - 66.

25. Ратникова Л.И. Методы расчета сейсмических полей в тонкослоистых средах. М.: Наука, 1973. 123 с.

26. Аки К., Ричардс П. Количественная сейсмология. Т. 2. М.: 1983. С. 751 - 857.

27. Papageorgiou A., Aki К. A specific barrier model for quantitative description of inhomogeneous faulting and the prediction of strong ground motion. 2. Applications of the model // Ibid. P. 953 - 978.

28. Hanks Т., McGuire R.K. The character of high - frequency strong ground motion // Ibid. 1981. Vol. 71. P. 2071 - 2095.

29. Brune J. Tectonic stress and the spectra of seismic shear waves from earthquakes // J. Geophys. Res. 1970. Vol. 75. P. 4997 - 5009.

30. О нормах проектирования сейсмостойких атомных станций / В.П. Татарников, М.Л. Клоницкий, С.Г. Трушин, А.К. Беляничев. Энергетическое строительство, 1983, N 7. С. 67 - 69.

31. Медведев С.В. Определение балльности интенсивности землетрясения. В кн.: Землетрясения в СССР. М.: Наука, 1961. С. 103 - 125.

32. Бугаев Е.Г. Выбор ограниченного набора акселерограмм для проектирования унифицированной АЭС и типового оборудования. Сейсмостойкое строительство. Отечественный и зарубежный опыт, 1982. Вып. 9. С. 4 - 9.

33. Бугаев Е.Г. Оценка расчетного сейсмического воздействия заданной обеспеченности для особо ответственных объектов // Вопросы инженерной сейсмологии. 1984. Вып. 25. С. 43 - 50.

Приложение 7

(справочное)

к Руководству "Определение

исходных сейсмических колебаний

грунта для проектных основ"

ПЕРЕЧЕНЬ

МАТЕРИАЛОВ ПО ГЛУБИННОМУ СТРОЕНИЮ ЗЕМНЫХ НЕДР

1. Комплект карт по геологическому строению и физическому состоянию недр - раздел 2 в Геологическом атласе России 1:10000000, М. - С.-Пб., 1996.

2. Карта - Общее сейсмическое районирование территорий Российской Федерации - ОСР - 97, масштаб 1:8000000. Составлена в Объединенном институте физики земли РАН в 1997 г.

3. Карта активных разломов земной коры России, масштаб 1:8000000. Издана в Геологическом институте (ГИН) РАН в 1996 г.

4. Карта современных вертикальных движений земной коры территории СССР, масштаб 1:5000000 М. ГУГК СССР, 1998 г.

5. Карта ландшафтов геохимических условий территории России, масштаб 1:4000000. Издана в Институте геологии рудных месторождений, петрографии, минералогии и геохимии (ИГЕМ) РАН в 1997 г.

6. Материалы сейсмического зондирования недр России по геотраверсам. Работа выполнена в Институте ВНИИ Геофизика и в Госцентре ГЕОН Министерства природных ресурсов РФ.

7. Материалы по сейсмостратиграфии. Шлезингер А.Е. Региональная сейсмостратиграфия, М.: Научный мир, 1998. Работа проводилась в ГИН РАН и в Госцентре ГЕОН.

8. Материалы по сейсмической томографии земной коры, литосферы и мантии. Геотектоника N 4, 1997 и Геотектоника N 1, 1998, с. 5.

Директор НТЦ ЯРБ

Б.Г.ГОРДОН