Создание синтетических акселерограмм для решения инженерных и градостроительных задач (стр. 4 )

9

8

7

6

5

4

3

2

lg (A2d) 0.5

Рисунок 12. Соотношение между сейсмической интенсивностью и логарифмом энергии сейсмической волны lg (A2d)

I

10

9

8

7

6

5

4

lg W

Рисунок 13. Соотношение между сейсмической интенсивностью и центрами распределений мощности сейсмической волны

5. Построение банка региональных синтетических акселерограмм для объектов повышенной ответственности

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

При построении базы акселерограмм для объектов повышенной ответственности использован чисто эмпирический метод оценки параметров сейсмического движения грунта. Заметим, что степень повреждения объектов определяется не только вибрациями грунта, но и остаточными деформациями грунта, например, при выходе разлома на дневную поверхность. Уровень ускорений на обеих горизонтальных компонентах выбран одинаковым, вертикальная компонента принята равной 2/3 от величины на горизонтальной компоненты [30].

Уровни колебаний соответствуют шкале РСШ -10 (таблица 2).

В качестве исходных спектров для генерации акселерограмм используются обобщенные спектры из СНиП [30] (см. раздел 3.2.1).

Чтобы избежать порблем с высокочастотной частью нормативного спектра (см. раздел 3.2.1), для обобщенного спектра для частот превышающих 100 Гц было принято усиление равным единице, а в частотном интервале 100 Гц – 10 Гц спектр был аппроксимирован линейной функцией в двойном логарифмическом масштабе. В области низких частот низкочастотный склон спектра протянут до 4 с для грунтов 1-й и 2-й категорий, а для грунтов 3-й категории - до 8 с. Последнее весьма важно для объектов с большими собственными периодами.

Уровень кривой “бета” по действующим СНиП принимается для расчета синтетической акселерограммы равным β = 2.5. Нормативный спектр характеризует спектры множества землетрясений с различными магнитудами, механизмами очага, расстояниями до объекта и поэтому значительно шире любого отдельного реального спектра.

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

Обобщенные нормативные и скорректированные нами спектры из СНиП показаны на рисунках 14 и 15.

β

10

1

0.1

0

T, с

Рисунок 14. Обобщенный нормативный спектр (красная линия) для грунтов 1-й и 2-й категории и его модификация (черная линия)

b

10

1

0.1

0

T, с

Рисунок 15. Обобщенный нормативный спектр (красная линия) для грунтов 3-й категории и его модификация (черная линия)

Хотя в СНиП спектры для грунтов 1-й и 2-й категории одинаковы, продолжительность колебаний на этих грунтах неодинакова. В наших расчетах этот факт учитывается (см. раздел 3.6.6). При расстояниях в ближней зоне берется значение ширины импульса на границе этой зоны. В среднем оценка продолжительности колебаний составляет d = 5 с. Общая продолжительность записи реального землетрясения примерно в 5 раз превышает ширину импульса.

Карты ОСР-97 соответствуют грунтам второй категории. По картам ОСР и СНиП грунты подразделяются на три категории таким образом, что бы для грунтов первой категории оценка балльности снижалась на один балл, а для грунтов третьей категории повышалась на один балл. Разделение грунтов на три категории при этом вполне оправдано. Увеличение числа градаций грунта по сейсмическим свойствам не соответствовало бы точности оценки сейсмической опасности в баллах. Заметим, что расчет сейсмических воздействий в ускорениях грунта имеет примерно такую же точность. Если схема расчета конструкции того требует, вводить пониженные коэффициенты «бета» следовало бы на этом этапе. При использовании обобщенного спектра рассматриваются все группы грунтов, поскольку от категории грунта зависит ширина импульса.

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

Работа с банком синтетических акселерограмм для объектов повышенной ответственности

А) По карте ОСР-97-В определяем сейсмическую интенсивность I в интересующем нас пункте. Для многих населенных пунктов оценки интенсивности приводятся в объяснительной записке к карте, а также в СНиП.

Б) Если после названия населенного пункта стоит звездочка, значит там проводилось СМР. По результатам СМР определяется категория грунта. Более надежно определяется категория грунта по результатам СМР непосредственно на площадке. СНиП допускает для объектов повышенной ответственности определять категорию грунта по результатам инженерно-геологических исследований.

В) Определяется сейсмическая интенсивность на площадке в зависимости от категории грунта. Для этого из оценки интенсивности по карте ОСР-97-В вычитаем один балл для грунтов 1-й категории, прибавляем один балл для грунтов 3-й категории и оставляем оценку на карте без изменений для грунтов 2-й категории.

Г) По таблице 7 по величине полученной для площадки интенсивности I и категории грунта на строительной площадке определяется № синтетической акселерограммы, соответствующей рассматриваемой площадке. Расчетные акселерограммы в цифровом виде приводятся в Приложении 3 (таблицы 1 – 9).

Таблица 7

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

Категория грунта	Сейсмическая интенсивность на площадке
7	8	9
1	Acc71 (табл. 1 Приложения 3)	Acc 81 (табл. 4 Приложения 3)	Acc 91 (табл. 7 Приложения 3)
2	Acc72 (табл. 2 Приложения 3)	Acc 82 (табл. 5 Приложения 3)	Acc 92 (табл. 8 Приложения 3)
3	Acc73 (табл. 3 Приложения 3)	Acc 83 (табл. 6 Приложения 3)	Acc 93 (табл. 9 Приложения 3)

6. Оценка параметров сейсмического движения грунта при ДСР

На предыдущем этапе (см. раздел 3.3) были описаны параметры сейсмического движения грунта, которые определяют сейсмический эффект при ожидаемых землетрясениях. Оценить значения этих параметров по СНиП и картам ОСР для конкретной площадки строительства невозможно. Оценка параметров сейсмического движения грунта является одной из главных задач при ДСР. Работы по детальному сейсмическому районированию проводятся при проектировании и строительстве особо ответственных объектов. В таких случаях проведение ДСР экономически целесообразно и необходимо. При этом, при назначении сейсмическх воздействий появляется возможность использовать дополнительные данные.

Напомним, какие характеристики ожидаемых землетрясений необходимо оценить при ДСР:

- ожидаемая магнитуда землетрясения;

- тип подвижки в очаге;

- кратчайшее расстояние от площадки до поверхности разлома.

При инженерно-геологических изысканиях и сейсмическом микрорайонировании оцениваются:

- тип грунта на площадке;

- резонансные свойства грунтов.

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

-  пиковое ускорение грунта PGA,

-  преобладающий период колебаний T,

-  ширина импульса (продолжительность колебаний),

-  логарифмическая ширина спектра,

-  коэффициент динамического усиления.

По этим параметрам строится синтетическая акселерограмма. Для всех этих параметров известно среднеквадратичное отклонение.

6.1. Пиковое ускорение грунта (PGA)

Вблизи разлома выделяются три зоны с различным затуханием: очаговая, ближняя и дальняя.

В пределах очаговой зоны пиковое ускорение грунта PGA = PGA0 = const. Строго говоря, в этой зоне амплитуда даже несколько возрастает с расстоянием, но в этой зоне большую роль играют нелинейные эффекты, связанные с высоким уровнем колебаний, поэтому проще считать PGA0 в этой зоне постоянной. Эта зона протягивается от поверхности разрыва до расстояния R0, км, которое определяется следующим образом:

lg R0 = 0.33 MS - 1.51.

Значение PGA0 зависит от типа подвижки в очаге (таблица 8).

Вторая зона (ближняя) протягивается от R0 до R1:

lg R1 = 0.33MS - 0.61,

Таблица 8

Коэффициенты в эмпирических соотношениях затухания

PGA0, см/с2	Тип подвижки
1070	Поддвиг (зоны субдукции)
900	Взброс
757	Взбросо-сдвиг
637	Сдвиг (S)
536	Сбросо-сдвиг
450	Сброс

где R1 - граница между ближней и дальней зонами, на которой PGA = 170 см/с2 для всех типов подвижки в очаге. Влиянием категории грунта на величину PGA в очаговой и ближней зоне можно пренебречь. Статистически установлено, что на рыхлых грунтах в этих зонах амплитуда даже несколько меньше, чем на твердых.

PGA в ближней зоне описывается эмпирической формулой:

lg (PGA) = C (0.33MS - 0.61- lg R) + 2.23 ± 0.16,

170 см/с2 < PGA < PGA0, R0<R<R1,

где значения PGA0 приведены в таблице 6.

В дальней зоне (lg R > lg R1; PGA < 170 см/c2) PGA описывается следующим выражением:

lg (PGA) = 0.634 MS - 1.92 lg R + 1.076 ± 0.20.

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

Уровень PGA с доверительным интервалом задается выбором карты: ОСР-97-А (5% вероятность превышения в течение 50 лет) и ОСР-97-В (1% вероятность превышения в течение 50 лет). Заметим, что в СНиП задаются явно устаревшие оценки ускорения грунта, фактически по результатам, полученным до 1951 г. Кроме того, следует иметь в виду, что при заданной сейсмической интенсивности величина PGA при заданной интенсивности зависит еще и от ширины импульса, что в СНиП совершенно не учитывается. Фактически пиковое ускорение – единственный параметр сейсмического движения грунта, используемый в СНиП.

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

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

lg (AHмакс /AH2)= 0.223 – 0.008 MS – 0.021 lg R ± 0.111.

На основании обработки 2556 записей мира оценено отношение уровней двух случайно ориентированных горизонтальных компонент. Если взять отношение более интенсивной горизонтальной компоненты ко второй, то средняя оценка отношения двух горизонтальных компонент составит 1.28, отношение амплитуд в 16% случаев не превысит 1.15, и в 50% случаев не превысит 1.65 [25]. Для практического использования удобно уровень амплитуд на обеих компонентах считать одинаковыми.

6.3. Соотношение амплитуд на горизонтальной и вертикальной компонентах

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

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

Таблица 9

Соотношение между амплитудами ускорений на вертикальной V и горизонтальной H компонентах

AH, см/с2 AV, см/с2	250 125	300 155	350 185	400 225	450 275	500 335	550 400
AH, см/с2 AV, см/с2	600 470	650 545	700 625	750 710	800 800	850 895	900 1000

Это соотношение практически не зависит от магнитуды, типа подвижки в очаге и категории грунта.

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

Величина стандартного отклонения составляет 0.16 ед. десятичного логарифма.

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

6.4. Коэффициент динамического усиления

Эта величина, как безразмерная, согласно теории размерностей и подобия должна быть весьма стабильной. Действительно, в работе [25] показано, что этот коэффициент для реальных землетрясений существенно выше нормативных. Среднее значение коэффициента динамического усиления при 5% демпфировании по мировым данным без учета магнитуды, типа подвижки в очаге, расстояния и категории грунта

lg β = 0.56 ± 0.10; (β = 3.63).

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

Отметим, что значение этого коэффициента в обобщенных спектрах, приведенных в СНиП, существенно меньше, что объясняют особенностями поведения зданий в неупругой стадии. Однако в более строгих нормах для объектов ядерной энергетики значение этого коэффициента принимается равным β = 3.2.

6.5. Преобладающий период колебаний

В СНиП преобладающие периоды косвенно задаются обобщенным спектром реакции (спектром ответа). Однако это задание произведено только с учетом категории грунта и совершенно не учитывает более важные факторы – магнитуду землетрясения и расстояние, не говоря уж о типе подвижки в очаге землетрясения.

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

Преобладающий период колебаний в дальней зоне (R > 0.33MS - 0.61) по среднемировой зависимости равен:

lg T = 0.15 MS + 0.25 lg R - lg vp + C1 + C2 ± 0.20,

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

R - гипоцентральное расстояние;

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

C1 = - 0.2 для зон субдукции, - 0.10 для взбросов, 0.00 для сдвигов и 0.10 для сбросов;

C2 = - 1.11 - член, определяющий влияние неучтенных факторов.

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

lg T = 0.15 Ms + 0.25 lg R + C1 – 1.9 ± 0.20,

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

В ближней и очаговой зонах преобладающий период можно считать величиной постоянной, такой же, как на границе ближней и дальней зон (lg R =0.33 M – 0.61).

В работе [97] показано следующее:

- видимый период T хорошо коррелируется с преобладающим периодом TRS соответствующего спектра реакции:

lg TRS = lg T ± 0.10

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

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

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

Величина продолжительности колебаний d, с - есть интервал времени между первым и последним разом превышения огибающей амплитуд 0.5 максимального значения.

Величина d, с в дальней зоне (R > 0.33 MS -0.61) по среднемировой зависимости равна:

lg d = 0.15 MS+ 0.50 lg R, км + C1 + C2 + C3 ± 0.30,

где С1 = 0.25 для сбросов, 0.12 для сбросо-сдвигов, 0.00 для сдвигов, - 0.12 для взбросо-сдвигов и - 0.25 для взбросов, в том числе и в зонах субдукции; С2 = - 0.15 для грунтов 1-й категории, 0.00 для грунтов 2-й категории и 0.25 для грунтов 3-й категории.

Среднее значение коэффициента С3 = - 1.3. Обработка данных по местным землетрясениям позволяет снизить величину стандартного отклонения до 0.25 ед. десятичного логарифма.

В ближней и очаговой зонах продолжительность колебаний можно считать величиной постоянной, такой же, как на границе ближней и дальней зон
(lg R =0.33 M – 0.61).

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

По эмпирическим данным нами найдено соотношение:

I = 2.5 lg A, см/с2 + 1.25 lg d, с + 1.05

6.7. Логарифмическая форма спектра (ширина спектра)

Если измерять ширину спектра не в Гц, а в октавах, что более физично, придем к логарифмической ширине спектра реакции S [118]:

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

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

Практика показала, что величина S для сильных движений слабо зависит от магнитуды, что соответствует и теории размерностей и подобия, поскольку логарифмическая ширина спектра является величиной безразмерной. Коэффициент корреляции величины S с магнитудой не превышает 0.13. Таково же значение коэффициента корреляции величины S с расстоянием. Среднемировое значение логарифмической ширины спектра при 5%-ном демпфировании около двух октав или

S = 0.60 ± 0.22

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

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

S = 0.58 + C1 + C2 ± 0.20,

где C1 = - 0.1 для взбросов, 0.0 для сдвигов и 0.1 для сбросов,

C2 = - 0.08 для грунтов 1-й категории, 0.00 для грунтов 2-й категории и 0.08 для грунтов 3-й категории.

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

7. Создание банка региональных синтетических акселерограмм для особо ответственных объектов

Как известно, синтетические акселерограммы определяются следующими факторами (см. раздел 3.5):

- магнитуда землетрясения;

- глубина очага;

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

- тип подвижки в очаге;

- расстояние до поверхности разрыва;

- резонансные свойства грунтов.

Каждый из перечисленных параметров имеет несколько градаций. Только одному весьма важному параметру – расстоянию будет соответствовать множество различных акселерограмм. Создание банка акселерограмм, соответствующих комбинациям всех параметров невозможно. Тем более, что ДСР проводится только для особо важных объектов, для которых в обязательном порядке проводится и микрорайонирование. Отсюда правило, хорошо известное всем изыскателям:

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

По приведенным в разделах 3.3, 3.5 формулам оцениваем параметры движений грунта на площадке для землетрясений в различных зонах ВОЗ.

7.1. Построение локального спектра

Построение нормированного по амплитуде локального спектра начинается с выбора коэффициента динамического усиления b. Если расчет ведется точно по СНиП, то выбирается нормативный коэффициент, если желательно получение синтетической акселерограммы с параметрами, близкими к реальным, то выбираем b = 3.4, или хотя бы
b = 3.2, поскольку это значение используется в нормах для строительства объектов ядерной энергетики. Это значение откладываем на ожидаемом преобладающем периоде T0. Поскольку оценка преобладающего периода содержит случайную составляющую, считаем, что максимум PGA с заданным уровнем доверия попадет на полочку T0 ± ns. При n=1, преобладающий период ожидаемого землетрясения с вероятностью 67% попадает в интервал lg T0 ± 0.2. От концов этого интервала задаются склоны спектра, которые определяются величиной S. Среднее значение склонов пропорционально периоду. При 67% уровне доверия непревышения уклона спектра угол с вертикалью составит arctg7/6. В соответствии с построенным спектром строится синтетическая акселерограмма с использованием любой компьютерной программы. Если для строительной площадки представляют опасность несколько зон ВОЗ, и воздействия от различных не слишком сильно различаются, возможно построение единой синтетической акселерограммы.

Синтетические акселерограммы должны удовлетворять критериям пригодности, приведенным, например, в Нормах для объектов ядерной энергетики (см. раздел 3.2.2.1):

Коэффициент корреляции различных компонент не должен превышать 0.3.

Спектр реакции (спектр ответа) синтетической акселерограммы не должен отличаться от целевого спектра более, чем на 10%.

Для особо ответственных объектов для более точного учета грунтовых условий в обязательном порядке проводится микрорайонирование. Построенные по результатм ДСР локальные спектры учитывают среднюю для каждой категории грунтов реакцию. В результате СМР оценивается АЧХ грунтов на площадке, которая накладывается на потроенный локальный спектр.

7.2. Синтетическая акселерограмма, соответствующая рассеянной сейсмичности

Для фоновой сейсмичности расчет соответствующей акселерограммы приходится производить отдельно вследствие особых условий: нулевое эпицентральное расстояние и малая магнитуда. Оценки воздействий от фоновой сейсмичности всегда консервативны, поскольку сейсмогенерирующие структуры неизвестны и приходится учитывать возможность возникновения очага непосредственно под объектом на малой глубине. Ранее уже говорилось, что в таких условиях высоким ускорениям соответствует относительно невысокая интенсивность. Например, при афтершоке Ташкентского землетрясения с интенсивностью 7 баллов была получена акселерограмма с ускорением около 700 см/с2.

8. Оценки скорости колебаний и смещений грунта

Скорость и смещение колебаний можно определить по инструментальной шкале сейсмической интенсивности:

lg PGV, см/с = - 2.23 + 0.47 I.

lg PGD, см = - 4.26 + 0.68 I.

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

Ширина импульса (продолжительности колебаний) скоростей превышает продолжительность колебаний в ускорениях в среднем на 41%. Продолжительность смещений превышает продолжительность колебаний в ускорениях в среднем в 2.5 раза.

СПИСОК ЛИТЕРАТУРЫ

1.  Аптикаев колебания при землетрясениях и взрывах. М.: Наука, 19с.

2.  Aptikaev F. F. Instrumental seismic intensity scale// Proc. EE-21 C, Skopje-Ohrid, 2005. CD: Topic 2, pp. 1-9.

3.  , Аптикаев шкал типа MSK // Вычислительная сейсмология, вып. С. 210 – 253.

4.  Boore, D. M. (2005). SMSIM – FORTRAN programs for simulating ground motions from earthquakes: version 2.3 – a revision of OFR 96 – 80 –A. USGS.

5.  Seismological Research Letters, 1997. Vol. 68, No 1.

6.  Аптикаев сейсмических волн в поглощающей среде// Геофизика и математика. Материалы 1-й Всероссийской конференции. М.: ОИФЗ РАН. 1999. С

7.  , Эртелева сильных движений грунта в сейсмически активных зонах Ирана // Физика Земли, 41, № 2, 2005. С. 23-29.

8.  , , и др. Природные опасности России. Сейсмические опасности. М.: Крук, 20с.

9.  Харкевич и анализ. М.: Гос. Издательство физико-математической литературы. 19с.

10.  Safety Series No 50-SG-S1 (Rev.1). Earthquakes and Associated Topics in Relation to Nuclear Power Plant Sitting. IAEA, Vienna, 1991.61 pp.

, , и др. Природные опасности России. Сейсмические опасности. М.: Крук, 20с.

12.  Strong Motion Data Processing/SIMQKE-1. NISEE Research Software. Univ. of California, Berkeley. 1998. CD 32-EBA15.

13.  Saragoni G. R., Hart G. C. Nonstationary Analysis and Simulation of Earthquake Ground Motions. Univ. of California, Rpt. UCLA –ENG-72pp.

14.  СНиП II-7-81* Строительство в сейсмических районах. М.: Минстрой России, 2000.

15.  , Михайлова спектра реакции в ускорениях для отдельного землетрясения // Вопр. инж. сейсмол. - М.: Наука, 1985.- вып. 26. - С.

16.  , Эртелева задания обобщенного спектра реакции для проектных основ// Сейсмостойкое строительство, безопасность сооружений. 2001, № 4. С. 4-7.

17.  Нормы проектирования сейсмостойких атомных станций НП-031-01. М.: Госатомнадзор России, 20с.

18.  Уломов -детерминированная оценка сейсмических воздействий на основе карт ОСР-97 и сценарных землетрясений // Сейсмостойкое строительство. 2005. №4. С.60-69.

19.  Аптикаев записей сейсмических колебаний// Вопросы инженерной сейсмологии, выпуск 21, 1981. С. 3-8.

20.  О сейсмологической основе норм сейсмостойкого строительства в России. // Физика Земли, 2002, № 12.

21.  , Шебалин корреляций между уровнем макросейсмического эффекта и динамическими параметрами движения грунта //Исследования по сейсмической опасности. Вопросы инженерной сейсмологии. Вып. 29. М.: Наука, 1988. С. 98-108.

22.  Определение исходных сейсмических колебаний грунта для проектных основ. (РБ-006-98). Нормативный документ. –М.: НТЦ ЯРБ, 20

23.  Ambraseys N. N., Srbulov M. Attenuation of Earthquake-Induced Ground Displacements. Earthquake Engineering and Structural Dynamics, vol. 23, 1994.

24.  , Эртелева искусственных акселерограмм методом масштабирования реальных записей// Физика Земли. , 2002, № 7, с. 39-45.

25.  Mikhailova N. N., Aptikaev F. F. Some Correlation Relations between Parameters of Seismic Motions// J. of Earthquake Prediction ResearchVol. 5. - № 5. - Pp.

26.  Aptikaev F. F., Erteleva O. O., Sacks M. V. Properties of Response Spectra. XXIV Assembly of IUGG, 2007, Perugia, Italy. Session SS004, presentation № 000.

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

28.  Mc Guire R. K. and Barnhard T. P. The usefulness of ground motion duration in predicting the severity of seismic shaking. Preprint. 19p.

29.  , , Мокрушина сейсмических шкал семейства Меркалли. Вулканология и сейсмология. 2008. № 3. С. 1 – 5.

30.  Строительство в сейсмических районах. (Первая редакция) и пояснительная записка. М.: Госстрой России, 20с.