CC BY
26
- © Е.Н. Черных, Г.В. Шубин,
Б.Н. Заровняев, 2015
УДК 622.221
Е.Н. Черных, Г.В. Шубин, Б.Н. Заровняев
ОЦЕНКА ПАРАМЕТРОВ ТЕХНОГЕННЫХ ВИБРАЦИЙ НА КАРЬЕРЕ «УДАЧНЫЙ»
Для оценки влияния вибрационных воздействий от различного вида горно-транспортного оборудования на массив уступов и берм бортов карьера «Удачный», в течение двух сезонов проведен комплекс натурных измерений с использованием цифровых инженерно-сейсмометрических станции «Байкал-1» и «Байкал-12». Ключевые слова: массив бортов карьера, инженерно-сейсмометрические станциии «Байкал-1» и «Байкал-12», вибрационное воздействие, горнотранспортное оборудование, виброскорость, виброперемещения, виброускорения, землетрясения обычной продолжительности.
На современных карьерах применяются мощные высокопроизводительные механизмы: экскаваторы, погрузчики, автосамосвалы, буровые станки и вспомогательное оборудование с различными виброхарактеристиками.
При отработке глубоких карьеров происходят динамические воздействия от естественных и техногенных вибраций на борта карьера.
Постоянной динамической нагрузке подвергаются откосы уступов, расположенные вблизи транспортных коммуникаций и мест локальной концентрации работающего выемочно-по-грузочного, бурового и вспомогательного оборудования.
Таблица 1
Для получения практических значений вибрационных воздействий на массивы уступов и берм бортов карьера трубки «Удачная» в начале 2000-х годов в течение двух летних сезонов проведены натурные исследования вибропараметров сейсмостанцией «Бай-кал-12» по всему периметру карьера [1].
Для записи вибраций применялись цифровые инженерно-сейсмометрические станциии «Байкал-1» - на первом этапе работ, проведенном, и «Байкал-12» - на втором этапе. Основные технические характеристики приведены в табл. 1.
В комлекте с сейсмоприемниками СК-1П (собственный период 1 с, зату-
№ пп Параметры «Байкал-1» «Байкал-12»
1 Количество дифференциальных входов 12 14
2 Разрядность АЦП 16 24
3 Базовый коэффициент преобразования 1,6 мкв/дискрет 1,0 мкв/дискрет.
4 Программируемые коэффициенты усиления (по каждому каналу) - 0, 2, 4, 8
5 Объем внутренней энергонезависимой памяти - 4 Мбт.
6 Частота дискретизации, программируемая 100 до 5000 Гц
7 Напряжение питания 12 вольт 9-15 вольт пост. тока (не стаб.)
8 Потребляемая мощность не более 6 ватт не более 6 ватт
Технические характеристики инженерно-сейсмометрических станций
хание 0,5), станция «Байкал-1» позволяет регистрировать скорость колебаний в диапазоне частот 0,7-20 Гц, а станция «Байкал-12» - в диапазоне частот 0,7-450 Гц. Интервал регистрируемых станциями амплитуд составляет от 0,003 до 200 мкм/с.
При измерениях запись вибраций велась одновременно в двух пунктах наблюдения двумя трехкомпонентными сейсмоприемниками. Один из сейсмо-приемников устанавливался на верхней бровке, другой - на нижней. Ориентация производилась по трем направлениям: Х - поперек бровки, У - вдоль бровки, I - в вертикальной плоскости. В каждом пункте наблюдения производилась запись длительностью не менее 5 минут. Схема расположения пунктов наблюдений измерения вибраций в карьере трубки «Удачная» показана на рисунке. Пункты наблюдений были выбраны таким образом, чтобы охватить все разновидности инженерно-геологических условий и установить их реакцию на динамическое вибровоздействие от основных техногенных источников, действующих в карьере.
В число пунктов наблюдений включены такие участки в карьере, как места установки тросово-сетчатой завесы, дополнительный выезд, северный борт с вертикальными уступами и т.д., в которых отмечены деструктивные явления, а также те участки, в которых можно ожидать резкого возрастания динамических нагрузок в случае сильных сейсмических и техногенных воздействий. Всего измерения были выполнены в 33 (в первый сезон) и 32 (во второй сезон) пунктах. Расстояние между соседними пунктами наблюдения ориентировочно принималось в пределах 200 м.
В карьере трубки «Удачная» постоянным источником виброизлучений, за исключением периодов времени, когда производятся взрывы, являются работающие автосамосвалы, выемоч-
но-погрузочные машины и буровые станки, обеспечивающие непрерывность технологического цикла.
Все техногенные источники вибрации, действующие в пределах карьера, можно разделить на две группы:
1. источники, которые потребляют энергию из сети или сами ее вырабатывают;
2. транспорт: автомобильный, строительно-дорожные машины и механизмы, бульдозеры, грейдеры.
Особенностью источников вибраций, относящихся к первой группе, является то, что режим работы источников этого типа определяется режимом работы энергосети, посредством которой эти источники связаны между собой. Для спектров сейсмических колебаний, вызываемыми такими источниками, характерно наличие узких пиков с шириной, близкой к такой, которую дали бы гармонические колебания. Частота гармоник связана с частотой энергосети коэффициентом деления на 2п (где п - число пар полюсов). К источникам такого типа будут относиться машины и механизмы, работающие от сети, такие как электродвигатели и генераторы, компрессора, насосы.
К источникам второй группы можно отнести карьерные автосамосвалы, бульдозеры на колесном и гусеничном ходу, автогрейдеры, автопогрузчики и другой вспомогательный автомобильный транспорт.
Виброисточники первой группы возбуждают и передают грунту квазистационарные, гармонические вибрации непрерывно или в течение длительного времени. Для сейсмических колебаний, возбуждаемых такими источниками, выполняется условие:
'с * Т >> 1, (1)
где / и Тс - соответственно средняя частота и длительность посылки сигнала.
Виброисточники второй группы возбуждают и передают в грунт ударные
Блоки с разной степенью нарушенное™
О Аномально awcotio нарушенный блок (12-15 и более систем трещин)
Условные обозначения
11 I | Границы блоков. сформировавшиеся & обстановке растяжения
Границы блоков, сформировавшиеся в обстановке сжатия @ нарушенные блоки (1<М2 систем трещин)
\- Границы блоков, сформировавшиеся в обстановке сдвига
ергднг трушенные Блок* |В-)0 систем трещин)
О&пони с рзэпячмэи степенью нарушенное™ тектоническим* /Г"^ Слабо нарушенные 6по<и систем трещин 1 трещинами. Бьщеляннь»е А И. Мйлькиковым
Зоны трещиноватое™ О Аномально «мшенный бпок [менее 0 систем трещин|
Схема расположения пунктов измерения вибраций в карьере трубки «Удачная»
нагрузки и возбуждают сейсмические колебания импульсного типа. Некоторые источники можно отнести как к первой, так и ко второй группе. Например, экскаваторы с электроприводом будут возбуждать сейсмические колебания вибрационного и импульсного типов. Такие источники вибраций, как карьерные самосвалы в процессе движения будут возбуждать колебания виброимпульсного типа.
При измерении параметров вибраций может использоваться несколько режимов. В режиме виброметра измеряется виброперемещение колебательного процесса. В режиме велоси-метра - скорость колебательного процесса B(t) = S(t).
В режиме акселерометра определяется ускорение A(t) = S(t). Указанные параметры вибрации могут быть измерены в одной плоскости или могут быть определены их векторные значения.
Например, в режиме виброперемещений векторное значение будет определяться как:
S(t ) xyz =JS(t )l + S(t )2y + S(t )Z (2)
где S(t)xyz - вектор, S(t)x, S(t)y, S(t)z - амплитуды виброперемещений по компонентам Х, У, Z.
При измерении параметров вибраций используют их пиковое (максимальное) значение или среднеквадратичное. При анализе результатов измерений в дальнейшем, за исключением
специально оговоренных случаев, будем использовать среднеквадратичные значения величин колебательного процесса. Эта величина, как правило, используется в практике виброизмерений, так как по физическому смыслу она связана с мощностью переменных составляющих случайного процесса [2]. Среднеквадратичное значение колебательной величины за период называют действующим:
11
V т
-I-1
{[хЦ)]2 Л
(3)
По этой причине среднеквадратичное значение колебательной величины называют иногда эффективным значением. В дальнейшем, под уровнем колебаний будем подразумевать среднеквадратичное значение колебательной величины (виброперемещения, виброскорости, виброускорения). При анализе результатов измерений используем следующие обозначения:
5 , 5 , 5 - максимальные амплиту-
х' у' г ^
ды виброперемещений по компонентам Х, У, Z■
В , В , В - максимальные амплиту-
х' у' г ^ у
ды виброскоростей по компонентам Х, У, Z■
А , А , А - максимальные амплиту-
х' у> г
ды виброускорений по компонентам Х, У, Z■
5стх, 5сту, 5стг - среднеквадратичные амплитуды (уровни) виброперемещений по компонентам Х, У, Z■ Вст , Вст , Вст - среднеквадратич-
х' у1 г 1 1
ные амплитуды (уровни) виброскорости по компонентам Х, У, Z■
Астх, Асту, Астг - среднеквадратичные амплитуды (уровни) виброускорения по компонентам Х, У, Z■
5ст , 5ст , 5ст - векторные зна-
ху^ хуг' хуг 1
чения среднеквадратичных амплитуд виброперемещений по компонентам Х, У, ^
Вст , Вст , Вст - векторные значе-
хуг хуг хуг
ния среднеквадратичных амплитуд виброскоростей по компонентам Х, У, Z■
Аст , Аст , Аст - векторные значе-
хуг хуг хуг
ния среднеквадратичных амплитуд виброускорений по компонентам Х, У, Z.
Следует отметить, что в настоящее время нет общепризнанной методики, позволяющей определять степень влияния вибраций на устойчивость уступов и бортов карьеров. Если в оценке действия взрывов на массив имеется сравнительно большой экспериментальный материал, подкрепленный теоретическими разработками, из которых вытекает, что основным параметром, определяющим сейсмическое воздействие взрыва на горные породы и сооружения, является колебательная скорость, то в отношении вибраций такого определенного параметра нет. Поэтому, на данном этапе исследований динамического влияния вибраций и их применения для целей прогноза устойчивости разумно включать в анализ, наряду с виброскоростью, такие кинематические элементы как виброперемещения и виброускорения.
Измерения вибраций в карьере производилось в два этапа с интервалом около года [3]. Измерения на втором этапе были проведены в тех же пунктах наблюдения (ПН), что и на первом этапе, за исключением ПН-33, где по техническим причинам измерения не были выполнены. Измерения амплитуд производились по открытому каналу, при различных полосах пропускания, обусловленных амплитудно-частотными характеристиками станций «Байкал-1» и «Байкал-12». Это обстоятельство было учтено при анализе результатов измерений. Следует отметить, что основным параметром, характеризующим стационарный эргодический процесс, к которому относятся вибрации, являются среднеквадратичные амплитуды. Максимальные амплитуды характеризуют в основном импульсные сигналы, и дают представление о верхней границе динамического диапазона вибраций. Сравнивая амплитуды сейсми-
ческих сигналов в «тихих» и «шумных» пунктах наблюдения, можно оценить динамический диапазон сейсмических сигналов, регистрируемых в карьере трубки «Удачная». Диапазон изменения зависит от регистрируемого кинематического элемента: смещения, скорости или ускорения. Виброперемещения в пределах проведенной серии измерений изменяются в 34 раза, виброскорости в 291 раз, а степень изменения виброускорений достигает 832 раза. Этот результат объясняется тем, что каждый из указанных вибропараметров характеризует сейсмические колебания в различной спектральной области. При заданной полосе частот уровень виброперемещений определяется ее низкочастотной, виброскоростей -среднечастотной, а виброускорений -высокочастотной частью спектра регистрируемых колебаний. Высокочастотные колебания, которые определяют виброускрорение, вследствие затухания уменьшаются с увеличением расстояния до источника сильнее, чем низкочастотные, которые в большей степени определяют, например, виброперемещения. В качестве основного параметра, характеризующего вибрационный процесс, используем среднеквадратичные амплитуды, которые, как уже отмечалось, для краткости называют уровнем.
Максимальный уровень виброперемещений по измерениям в первой серии наблюдался в ПН-27б по компоненте Х, ствх = 0,24 мкм, виброскоростей по компоненте Х, стЬх = 0,23 мкм/с, виброускорений по компоненте У, стАу = 231,9 мкм/с2. Величина отношения уровня вибоперемещений по горизонтальным компонентам к уровню на вертикальной компоненте изменяется от 0,2 до 6,4. Максимальные значения этого отношения, характеризующего поляризацию колебаний, наблюдаются
в ПН-29, ст /ст = 4,7, ст /ст = 6,4, ми' вх вг ' ' вг ' '
нимальное значение в ПН-5, ствх/ствг = = 0,3, ст /ст = 0,2.
' ' ву вг '
По виброскорости диапазон отношений уровня колебаний по горизонтальным компонентам к вертикальной составляет от 0,35 до 3,64, по виброускорениям от 0,5 до 4,3. Максимальные значения отношений по виброскорости и виброускорению наблюдаются в ПН-29, ст/ст = 3,6, ст /ст =
' ' вг ' ' ву вг
= 3,5, ст /ст = 4,3, ст /ст = 3,9, ми' ' ах sz ' ' ау аг ' '
нимальные в ПН-5 по виброскорости, ст /ст = 0,5, ст /ст = 0,4 и по вибро-
вх sz ' ' ву- вг ' 1
ускорению в ПН-19, стаx/стsz = 0,8,
sx sz
ст /ст = 0,5.
ау аг '
Максимальный уровень виброперемещений по измерениям во второй серии наблюдался в ПН-27б по компоненте Х, ствх = 0,48 мкм, виброскоростей в ПН-28а по компоненте У, стЬх = 20,98 мкм/с, виброускорений в ПН-31б по компоненте У, ст =
Ау
= 525,8 мкм/с2. Величина отношения уровня вибоперемещений по горизонтальным компонентам к уровню на вертикальной компоненте изменяется от 0,4 до 3,2. Максимальные значения этого отношения наблюдаются в ПН-27б, ст /ст = 3,2 и в ПН-4б, ст /ст =
' вх вг ' ву вг
= 4,4, минимальное значение в ПН-20б, ст /ст = 0,36 и в ПН-21а, ст /ст =
вх вг ву вг
= 0,17. По виброскорости диапазон отношений уровня колебаний по горизонтальным компонентам к вертикальной составляет от 0,37 до 4,88, по виброускорениям от 0,23 до 4,76 Максимальные значения отношений по виброскорости наблюдаются в ПН-24а, ст /ст = 3,9 и в ПН-31б ст /ст =
' вх вг ' ву вг
= 4,76, по виброускорению в ПН-24а ст /ст = 3,9 и ПН-31б ст /ст = 4,76,
ах вг ау аг
минимальные в ПН-19а по виброскорости, ствх/ствг = 0,79 и в ПН 10б, ству/ствг = 0,37. По виброускорению в вуПН-19а, ст/ст = 0,5 и ПН-10б
ах вг
ст /ст = 0,23.
ау аг
Для анализа распределения амплитуд вибраций в карьере и уменьшения количества анализируемых параметров значения среднеквадратичных амплитуд по компонентам X, У, Z были пе-
реведены в векторные значения. Векторные значения вибропараметров приведены в табл. 2 и 3.
При анализе результатов измерений рассмотрены среднеквадратичные амплитуды (уровни) колебательных величин и их векторные значения.
Значения вибропараметра, превышающего средний уровень на величину более чем среднеквадратичное значение (стандарт). Такие значения вибропараметров, то есть превышающие среднее значение на величину более чем стандарт, выделены жирным шрифтом в табл. 2 и 3.
Максимальные амплитуды характеризуют, в основном, импульсные сигналы и дают представление о верхней границе динамического диапазона вибраций.
В качестве критерия повышенной вибрационной нагрузки приняты основании введенного критерия аномальные значения уровня вибропараметров по первой серии измерений наблюдались: по виброперемещениям - в ПН-1, ПН-4, ПН-29, виброскоростям -в ПН-1, ПН-8, ПН-27 ПН-29, ПН-32, виброускорениям - в ПН-8, ПН-27, ПН-31; ПН-32. По измерениям второй серии повышенные динамические нагрузки наблюдались в следующих пунктах наблюдения: по виброперемещениям - в ПН-1а, ПН-4а, ПН-10а, ПН-27а; ПН-28а, по виброскоростям -в ПН-1а, ПН-10а, ПН-28а, ПН-29а, ПН-30а, ПН-31а, по виброускорениям -в ПН-10а, ПН-27а, ПН-31а. В указанных пунктах повышенные динамические нагрузки на борта карьера могут
Таблица 2
Векторные значения вибропараметров по измерениям первого сезона
Наим. борта, абс. отм. № ПН 5охуг, мк м Во , хуг7 мк м/с Аохуг, мк м/с2 № ПН мк м Во , хуг7 мк м/с Аохуг, мк м/с2 Наим. борта, абс. отм.
Южный борт, +115:+70м 1 0,21 2,64 111,82 17 0,02 0,34 14,89 Северн. борт, гор.+160 м
2 0,05 0,86 51,08 18 0,02 0,59 34,43
3 0,06 1,13 70,06 19 0,05 0,69 31,59 Зап. борт, +70 м
Вост. борт. -35 м 4 0,16 1,66 63,19 20 0,09 1,13 53,48 Северн. борт, горизонт +70 м
Северный борт, горизонт -35 м 5 0,07 0,76 29,33 21 0,05 0,81 44,94
6 0,04 0,61 32,62 22 0,04 0,82 49,66
7 0,05 0,95 60,49 23 0,08 1,26 74,79
8 0,05 1,95 155,20 24 0,07 1,32 75,75
9 0,08 1,44 100,22 25 0,07 1,11 59,29
Юго-вост. борт, -35 м 10 0,04 0,86 55,80 26 0,08 1,21 66,27
Южн. борт, +130 м 11 0,06 1,21 71,83 27 0,20 2,76 207,65 Западн. борт, +70 : +55 м.
Западный борт, +160 : +145 м 12 0,05 1,22 94,14 28 0,07 0,90 42,72
13 0,04 0,60 32,07 29 0,21 2,21 105,85 Южный борт, Горизонт +155 : -20 м
14 0,05 1,02 67,57 30 0,05 1,05 66,81
15 0,03 1,02 51,26 31 0,04 1,66 126,73
16 0,03 0,76 51,58 32 0,05 2,18 196,93
33 0,09 1,64 105,16 Дно карьера гор. -155 м
Таблица 3
Векторные значения вибропараметров по измерениям второго сезон
Наим. борта, абс. отм. № ПН мк м Во , хуг7 мк м/с мк м/с2 № ПН мк м Bо , хуг7 мк м/с Aохyz, мк м/с2 Наим. борта, абс. отм.
Южный борт, +115 ++70м 1а 0,18 6,82 213,4 17а 0,04 1,01 48,7 Северн. борт, гор. +160 м
2а 0,08 2,45 84,5 18а 0,05 1,87 107,8
3а 0,04 1,09 50,3 19а 0,09 2,76 130,3 Зап. борт, +70 м
Вост. борт. -35 м 4а 0,27 4,78 174,9 20а 0,08 3,47 102,9 Северн. борт, горизонт +70 м
Северный борт, горизонт -35 м 5а 0,05 1,59 73,4 21а 0,11 2,54 111,1
6а 0,06 2,32 113,5 22а 0,09 4,28 132,0
7а 0,07 4,34 218,8 23а 0,08 3,15 114,7
8а 0,07 3,28 143,1 24а 0,07 3,75 182,9
9а 0,05 2,41 156,6 25а 0,08 2,02 74,5
Юго-вост. борт, -35 м 10а 0,39 14,98 508,7 26а 0,07 2,65 93,3
Южн.борт, +130 м 11а 0,06 2,03 72,3 27а 0,41 2,12 389,9 Западн. борт, +70 + +55 м
Западный борт, + 160 : +145 м 12а 0,08 2,46 71,5 28а 0,30 20,61 175,9
13а 0,04 1,41 52,8 29а 0,18 12,38 130,5 Южный борт, горизонт +155 + -20 м
14а 0,07 1,78 71,1 30а 0,05 10,57 189,2
15а 0,04 1,34 49,7 31а 0,03 14,09 490,3
16а 0,04 1,44 81,2 32а 0,02 0,93 40,9
служить дополнительным фактором, влияющим на величину и распределение естественных напряжений и привести к нарушению прочностных свойств массива.
Анализ полученных результатов вибрационного воздействия горного и транспортного оборудования позволил установить, что наиболее подвержены вибрационным нагрузкам в регистрируемые периоды времени были южный и западный борта карьера, во второй серии - Юго-восточный борт [4].
Для самых общих оценок можно воспользоваться эмпирической связью между максимальным ускорением колебаний, их длительностью и сейсмическим эффектом. Согласно этим данным увеличение длительности колебаний в п раз эквивалентно изменению сейсмического эффекта в раз [5]. Если воспользоваться этой, весьма приблизительной для наших условий оценкой,
то получается, что рассматриваемые длительные вибрационные воздействия примерно эквивалентны 2-3-балльным землетрясениям обычной продолжительности. Этот вывод, не претендуя на точность, тем не менее, показывает, что наблюдаемые динамические нагрузки, хотя и не способны вызвать мгновенные деформации, также могут быть одной из причин ведущих к катастрофическим последствиям и поэтому заслуживают самого тщательного изучения.
Таким образом, разработанная система мониторинга на карьере тр. «Удачная» позволяет осуществлять систематический контроль за значительной группой различных факторов, оказывающих определяющее значение на безопасность при ведении горных работ. По системе мониторинга, в отдельных случаях, можно прогнозировать и оперативно влиять на реально создающуюся неблагоприятную ситуацию на карьере.
СПИСОК ЛИТЕРАТУРЫ
1. Шмырко А.Н., Александров И.Н., Шубин Г.В. и др. Разработка системы мониторинга для безопасного ведения горных работ в условиях карьера «Удачный» // Горный информационно-аналитический бюллетень. -2004. - № 3 - С. 251-259.
2. Клюев В.В. Приборы и системы для измерения вибраций, шума и удара. Справочник в 2 томах. Т. 1. - М.: Машиностроение, 1978. - 447 с.
3. Черных Е.Н., Павленов В.А., Александров И.Н. Оценка устойчивости бортов карьера трубки «Удачная» / Материалы Между-
народной конференции 6-9 октября 2003 г. «Геодинамика и напряженное состояние недр земли». - Новосибирск, 2003. - С. 308-312.
4. Александров И.Н., Шмырко А.Н., Шубин Г.В., Кирюшин Д.И. Создание безопасных условий отработки сверхглубоких карьеров Якутии (на примере доработки карьера трубки «Удачная»). - Новосибирск: Наука, 2005. - 180 с.
5. Штейнберг В.В. и др. Методы оценки сейсмических воздействий // Вопросы инженерной сейсмологии. - 1993. - Вып. 34. -С. 5-94. ЕШ
КОРОТКО ОБ АВТОРАХ_
Черных Евгений Николаевич - кандидат технических наук, старший научный сотрудник, Институт Земной Коры, г. Иркутск,
Шубин Григорий Владимирович1 - кандидат технических наук, зав. кафедрой, Заровняев Борис Николаевич1 - доктор технических наук, директор, 1 Горный институт СВФУ.
UDC 622.221
ABOUT USE OF SCRAPER ON ALLUVIAL DEPOSITS
Chernih E.N., Candidate of Technical Sciences, Senior Research,
Institute of Earth Crust, 664033, Irkutsk, Russia,
Shubin G.V.1, Candidate of Technical Sciences, Head of Department,
Zarovnyaev B.N.1, Doctor of Technical Sciences, Director,
1 Mining Institute, North-Eastern Federal University named after M.K. Ammosov,
678015, Yakutsk, Russia.
To evaluate the effect of vibration effects on various types of mining transport equipment to an array of benches and berms pit walls «Lucky» for two seasons, a complex in-situ measurements using digital engineering seismic station «Baikal-1» and «Baikal-12».
Key words: An array of pit walls, engineering seismic station «Baikal-1» and «Baikal-12» vibration exposure, mining equipment, velocity, vibration displacement, vibration acceleration, the usual duration of the earthquake.
REFERENCES
1. Shmyrko A.N., Aleksandrov I.N., Shubin G.V. Gornyy informatsionno-analiticheskiy byulleten'. 2004, no 3, pp. 251-259.
2. Klyuev V.V. Pribory i sistemy dlya izmereniya vibratsiy, shuma i udara. Spravochnik v 2 tomakh, T. 1 (Instruments and systems for measuring vibrations, noise and impact. Reference book in two volumes, vol. 1), Moscow, Mashinostroenie, 1978, 447 p.
3. Chernykh E.N., Pavlenov V.A., Aleksandrov I.N. Materialy Mezhdunarodnoy konferentsii «Geodina-mika i napryazhennoe sostoyanie nedr zemli», 6-9 oktyabrya 2003 g. (Proceedings of International Conference on Geodynamics and Stress State of the Earth's Interior, 6-9 October 2003), Novosibirsk, 2003, pp. 308-312.
4. Aleksandrov I.N., Shmyrko A.N., Shubin G.V., Kiryushin D.I. Sozdanie bezopasnykh usloviy otrabotki sverkhglubokikh kar'erov Yakutii (na primere dorabotki karera trubki «Udachnaya») (Creation of safe mining conditions in super deep open pit mines in Yakutia (in terms of Udachnaya pipe development)), Novosibirsk, Nauka, 2005, 180 p.
5. Shteynberg V.V. Voprosy inzhenernoy seysmologii. 1993, issue 34, pp. 5-94.
