Инженерная методика оперативной оценки сейсмовзрывного воздействия на законтурный массив в карьере
С.А. Козырев, Е.А. Власова, Е.А. Усачев
Горный институт Кольского научного центра Российской академии наук, г. Апатиты, Российская Федерация
Горная Промышленность №5 / 2024 стр.66-73
Резюме: С целью определения сейсмического воздействия массовых взрывов на непосредственно прилегающий к инициируемому блоку массив горных пород были проведены замеры скорости смещения грунта в ближней зоне действия взрыва на различных участках карьера «Железный» АО «Ковдорский ГОК». По степени реакции массива на динамические воздействия произведено районирование карьера по глубине. Выделены три области сейсмичности, для которых получены эмпирические зависимости изменения скорости смещения от приведенного расстояния и соответствующие коэффициенты сейсмичности и затухания. По ним рассчитаны размеры зоны наведенной трещиноватости относительно допустимой скорости смещения. Средний размер отдельности пород увеличивается с глубиной карьера: сначала резко (от 0,22 до 1,19 м) в интервале 0–150 м от дневной поверхности, затем постепенно (от 1,19 до 1,4 м) до глубины 450 м, а глубже стабилизируется на уровне 1,4–1,5 м, что и дает различие в размерах зон наведенной трещиноватости. Разработана инженерная методика оперативной оценки сейсмовзрывного воздействия на законтурный массив при производстве взрывных работ в карьере для расчета ширины зон заколов, трещинообразования в тыльную часть массива и зон проявления наведенной трещиноватости. Расчет на стадии проектирования позволяет подобрать параметры взрывания, обеспечивающие минимизацию сейсмического воздействия на законтурный массив. Для оперативного расчета положения границ разрушения в тыльную часть массива при взрыве скважин разработана компьютерная программа. Полученные расчетные значения хорошо согласуются с ранее полученными данными моделирования и практики.
Ключевые слова: карьер, взрывчатые вещества, скважинные заряды, массовый взрыв, сейсмическое действие взрывов, устойчивость уступов
Для цитирования: Козырев С.А., Власова Е.А., Усачев Е.А. Инженерная методика оперативной оценки сейсмовзрывного воздействия на законтурный массив в карьере. Горная промышленность. 2024;(5):66–73. https://doi.org/10.30686/1609-9192-2024-5-66-73
Информация о статье
Поступила в редакцию: 17.07.2024
Поступила после рецензирования: 27.08.2024
Принята к публикации: 16.09.2024
Информация об авторах
Козырев Сергей Александрович – доктор технических наук, главный научный сотрудник, заведующий лабораторией, Горный институт Кольского научного центра Российской академии наук, г. Апатиты, Российская Федерация; e-mail: Адрес электронной почты защищен от спам-ботов. Для просмотра адреса в вашем браузере должен быть включен Javascript.
Власова Елена Анатольевна – кандидат технических наук, старший научный сотрудник, Горный институт Кольского научного центра Российской академии наук, г. Апатиты, Российская Федерация; е-mail: Адрес электронной почты защищен от спам-ботов. Для просмотра адреса в вашем браузере должен быть включен Javascript.
Усачев Евгений Андреевич – ведущий технолог, Горный институт Кольского научного центра Российской академии наук, г. Апатиты, Российская Федерация; е-mail: Адрес электронной почты защищен от спам-ботов. Для просмотра адреса в вашем браузере должен быть включен Javascript.
Введение
Горнодобывающая промышленность в настоящее время претерпевает значительные изменения: повышаются требования к производительности технологических процессов, промышленной безопасности, охране окружающей среды. Взрывные работы – важнейший процесс горной технологии, от которого зависят последующие процессы добычи и переработки полезного ископаемого. Совершенствование технологии ведения взрывных работ для управления дробящим и сейсмическим воздействием по-прежнему актуально. Необходимость поддержания постоянного качества добываемой руды заставляет прибегать к помощи полевых экспериментов, подготовить и провести которые все труднее и труднее или даже невозможно. Численное моделирование становится важным инструментом для прогнозирования взаимодействия горной породы и взрывчатых веществ.
Особенности моделирования процесса разрушения горных пород взрывом
Под действием взрыва в массиве горных пород образуются зоны дробления (полного разрушения); заколов (перемещения); остаточных деформаций (мелких трещин) и микротрещин (сотрясения).
Зона дробления характеризуется полным разрушением горных пород на отдельные куски разных размеров. Параметры этой зоны зависят от размеров взрываемого блока и густоты расположения в нем взрывных скважин. Последующие зоны находятся за пределами блока в массиве.
Зона заколов начинается от бровки уступа и характеризуется интенсивной нарушенностью сплошности массива с образованием трещин (заколов), по которым смещаются отдельные блоки относительно друг друга. Прочность массива здесь снижается в десятки раз. Размеры этой зоны колеблются в пределах от нескольких метров до 10 и более. К зоне заколов примыкает зона остаточных деформаций, в которой формируются мелкие трещины с незначительными смещениями блоков горных пород. В результате сцепление пород уменьшается на 20–50%, силы трения по ослабленным поверхностям значительно снижаются. Ширина зоны может достигать десятков метров.
Зона сотрясений характеризуется развитием микротрещин и упругих колебаний, которые интенсивно затухают на расстоянии 65–80 м от места взрыва. Микротрещиноватость ослабляет массив, а упругие колебания при отдельных условиях могут оказывать отрицательное влияние на устойчивость откосов уступов, расположенных выше и ниже взрываемого блока.
Размеры зон взаимосвязаны и определяют границы структурного ослабления массива, в которых заметно изменяются его деформационно-прочностные характеристики и проницаемость. Размеры зависят от многих факторов, главные из них – количество одновременно врываемого взрывчатого вещества (ВВ), структурно-тектонические особенности массива, прочностные свойства горных пород и параметры буровзрывных работ.
В условиях крупномасштабных горных работ в карьере рудника «Железный» АО «Ковдорский ГОК» при разработке мероприятий по обеспечению длительной устойчивости уступов и бортов карьеров, особенно при постановке их в предельное положение, знание размеров этих зон позволит обосновать ширину приконтурных блоков и степень их нарушенности.
Прогнозирование результатов взрыва является очень сложной задачей. Математически она решается системами дифференциальных уравнений при наличии начальных и граничных условий. К настоящему времени выполнен большой объем исследований и разработано множество математических моделей разрушения, среди которых можно выделить модели: Е.И. Шемякина, В.В. Адушкина, В.Н. Родионова, В.А. Боровикова, И.Ф. Ванягина, М.Г. Менжулина, Н.Н. Казакова, Г.М. Крюкова, А.И. Шера и многих других. Единого представления о механизме разрушения среды до сих пор нет, что связано главным образом с многообразием, сложностью и скоротечностью явлений, сопровождающих взрыв [1–8]. Следует отметить, что указанные модели предназначены для расчета параметров волн напряжений в безграничном массиве при камуфлетном взрывании удлиненного заряда и не позволяют проследить всю динамику процесса при взрыве. Поэтому они не нашли широкого применения в практике ведения взрывных работ.
Решение задачи прогнозирования возможно численными методами. Вручную посчитать задачи не представляется возможным. Для расчета параметров взрывных работ и прогнозирования их результатов применяют численное моделирование с использованием компьютерных программ. Программное обеспечение способно выполнить моделирование взрыва, сделать расчет работы взрыва, прогнозировать смещение горной массы, гранулометрический состав, выход и количество негабаритов, оценить сейсмическое воздействие.
На сегодняшний день существует не так много верифицированных программных продуктов, позволяющих с высокой степенью достоверности моделировать процессы, происходящие при взрыве. Известным разработчиком программных решений является компания Ansys с программными продуктами LS-Dyna и Autodyn, которые позволяют моделировать весь процесс взрыва: детонацию заряда, расширение зарядной полости, истечение продуктов детонации, механическое взаимодействие продуктов детонации со средой, формирование и распространение ударных волн и волн напряжений в твердой среде, разрушение среды, сдвижение разрушенного материала и его фрагментацию 1. Компания MineExcellence (Австралия) – разработчик программных технологий для буровзрывных работ – также декларирует варианты решения проблем, возникающих во время производства взрывных работ: управление разлётом горной массы и гранулометрическим составом, определение сейсмического эффекта взрыва и направления распространения взрывных волн 2.
Из современных технологических решений можно выделить имитационную модель I-Blast Thierry Bernard Technologie (Франция) [9]. В модели используется численный метод расчёта – генетический алгоритм в сочетании с методом дискретных элементов. Задача решается в динамической постановке, результатом является модель взрыва в развале, которую можно оценить до реального взрыва.
Несмотря на широкие возможности программных продуктов они требуют больших навыков работы с ними и наличия компьютерных систем, что не всегда доступно промышленному сектору и не позволяет на стадии проектирования оценить все варианты взрывания. Работа по расчёту рациональных параметров БВР осложняется тем, что необходимо учитывать всё многообразие горно-геологических и горнотехнических условий. Поэтому для практиков необходимо иметь инженерные методики, позволяющие оперативно оценивать воздействие взрыва на массив горных пород в реальных условиях.
Для инженерных расчетов наиболее распространенной моделью по оценке зон разрушения вокруг удлиненного скважинного заряда является модель Holmberg and Persson (H&P) [10; 11], в которой учитывается геометрия расположения заряда относительно точки наблюдения, а оценка зон разрушения производится по критическим скоростям смещения, характерным для рассматриваемых горных пород (рис. 1).
Общая зависимость определения скорости смещения PPV исходя из модели Holmberg and Persson выглядит следующим образом:
(1)
где К – коэффициент грунтовых условий, зависящий от свойств пород в точке исследования и от типа ВВ; α – показатель степени затухания колебаний; Р – вместимость ВВ в погонном метре скважины, кг; R0 – расстояние от источника колебаний до исследуемой точки в плане, м; Н – длина заряда ВВ, м; XS – длина забойки, м; X0 – глубина расположения исследуемой точки, м.
Критическую скорость смещения PPVcrit рассчитывают по зависимости:
(2)
где σt – предел прочности горной породы на растяжение, МПа; VP – скорость распространения продольных волн, м/с; Е – модуль Юнга, ГПа.
О характере разрушения вокруг заряда судят по величине скорости смещения в рассматриваемой точке. При скорости смещения больше 4000 мм/с (8 PPVcrit ) считается зона раздавливания, при скорости от 1000 до 4000 мм/с (4 PPVcrit) – зона интенсивного дробления, при PPVcrit – граница зоны развития радиальных трещин.
Данная модель также предназначена для расчета скоростей смещения при взрыве в безграничном массиве и требует экспериментальным путем определять коэффициент грунтовых условий К и показатель степени затухания колебаний α для реальных условий.
Из отечественных моделей для инженерных расчетов наиболее известна модель М.И. Щуйфера и Е.А. Азарковича [12], в которой критериями разрушения являются критическая и допустимая скорости смещения в волне сжатия.
Отличительной особенностью данной модели является то, что она позволяет рассчитывать радиусы трещинообразования на уровне торцов и середины заряда, а также для точек у свободной поверхности заряда при различной его длине на основе экспериментальных данных, что соответствует реальным условиям взрывания. Зная размеры этих зон, можно прогнозировать линию отрыва и заколов, а следовательно, точнее подойти к обоснованию сетки скважин по сравнению с принятой формой зоны дробления для отдельного заряда в виде вертикального полуцилиндра. Поэтому нами за основу при определении ширины зоны трещинообразования и заколов в законтурном массиве принят методический подход, предложенный М.И. Шуйфером и А.Е. Азарковичем.
Методика расчета ширины зоны заколов и интенсивного разрушения
Радиус зоны заколов и трещинообразования от взрыва одиночного скважинного заряда определяется по преобразованной в [12] к единому виду зависимости:
(3)
где r – радиус зоны разрушения, м; С – коэффициент, учитывающий глубину рассматриваемой точки от взрываемого заряда: вглубь по оси заряда С = 0,7; на уровне торца заряда С = 1,0; на уровне середины заряда С = 1,2; на уровне верхней части заряда С = 1,1; на уровне дневной поверхности С = 1,3; d – диаметр заряда, м; Vкр – критическая скорость смещения в волне сжатия, при которой происходит разрушение массива с учетом категорий пород по трещиноватости и коэффициента структурного ослабления; α1 – показатель затухания скорости смещения от глубины расположения точки от заряда (на уровне торца и вглубь по оси заряда α1 = 0,125; на уровне середины заряда α1 = 0,25; на уровне верхнего торца заряда α1 = 0,375; на уровне дневной поверхности α1 = 0,8); α2 – коэффициент, учитывающий число свободных поверхностей (при двух поверхностях α2 = 1, при одной α2 = 1,5); ρ – плотность заряжания, кг/м3; ε – коэффициент, зависящий от вида ВВ: для граммонита 79/21 – 1, для эмулитов – 0,85.
Критическую скорость смещения в волне сжатия Vкр, при которой происходит разрушение массива, следует рассчитывать по зависимости, предложенной Ф.М. Ляховицким [13]:
(4)
где кр – отношение предела прочности горной породы на растяжение к пределу прочности на сжатие; Ср – скорость распространения продольных волн в массиве горных пород, м/с; кд – коэффициент динамичности (увеличение предела прочности горной породы на растяжение при динамических нагрузках. Для скальных горных пород кд = 1,5); кто – коэффициент ослабления, учитывающий степень раскрытия трещин для различных категорий пород по трещиноватости (табл. 1); μ – коэффициент Пуассона.
По зависимости (4) рассчитаны критические скорости смещения для горных пород Ковдорского месторождения (табл. 2) в зависимости от категории по трещиноватости с учетом коэффициента структурного ослабления (см. табл. 1). Прочностные и динамические характеристики пород приведены по данным лаборатории геомеханики Горного института КНЦ РАН.
Для оперативного расчета положения границ заколов и трещинообразования в тыльную часть массива при взрыве скважин разработана компьютерная программа. Порядок расчетов по программе.
1. Выбираются тип руды или породы, категория пород по трещиноватости и критическая скорость смещения.
2. Задаются диаметр скважин и плотность заряжания.
3. Производится расчет.
4. Программа автоматически отрисовывает положение границ заколов и трещинообразования в тыльную часть массива на различных высотных уровнях от взрываемого заряда (рис. 2).
Полученные расчетные значения хорошо согласуются с данными математического моделирования (рис. 3). Моделирование производилось в программном продукте Ansys Autodyn, с помощью которого нами уже решались задачи взрывного дела [14; 15].
Положение границы зоны заколов и трещинообразования (пунктирная линия) принималось таким, когда в этой зоне степень нарушенности массива достигала более 50% по шкале нарушений. Эту же линию можно принять и как линию отрыва, что хорошо согласуется с данными практики (рис. 4).
Методика расчета зон проявления наведенной трещиноватости в прилегающем массиве горных пород при производстве массовых взрывов
К зоне заколов примыкает зона наведенной трещиноватости. Размеры этих зон взаимосвязаны и определяют границы структурного ослабления массивов, в которых заметно изменяются его деформационно-прочностные характеристики. В результате этого сцепление пород уменьшается на 20–50%, силы трения по ослабленным поверхностям значительно снижаются. Ширина зоны может достигать нескольких десятков метров. При применении технологии щадящего взрывания при постановке бортов карьера в конечное положение ширина приконтурных блоков должна соответствовать радиусу зоны проявления наведенной трещиноватости в прилегающем массиве горных пород от взрывов зарядов дробления. Поэтому определение размера этой зоны является необходимым условием при постановке уступов в предельное положение. Существующие подходы к определению этих зон являются весьма ориентировочными и не учитывают глубину расположения взрываемого заряда, а зона разупрочнения принимается одинаковой по всей высоте уступа [7].
Основной причиной разупрочнения породного массива при взрыве является сейсмовзрывное воздействие, интенсивность которого принято оценивать величиной скорости смещения породы в рассматриваемой точке на определенном удалении от заряда.
Наиболее объективными являются результаты, полученные с помощью сейсмоизмерений в ближней зоне действия взрыва в различных зонах карьера. В качестве критерия оценки зоны наведенной трещиноватости в ближней зоне принимаются допустимые скорости смещения в зависимости от категории пород по трещиноватости во взрываемом блоке, которые, как правило, составляют половину от критической скорости (см. табл. 2).
С целью определения сейсмического воздействия массовых взрывов на непосредственно прилегающий к инициируемому блоку массив горных пород были проведены замеры скорости смещения грунта в ближней зоне действия взрыва на различных участках по глубине карьера. Карьер по глубине был разделен на три области. Первая область – верхняя часть карьера (+300 м / +100 м), вторая область – средняя часть карьера (+100 м / –100 м) и третья область – нижняя часть (–100 м / –300 м).
Для указанных областей были рассчитаны и построены зависимости скорости смещения от приведенного расстояния (рис. 5) и получены соответствующие коэффициенты сейсмичности и затухания для каждой области.
Далее рассчитаны зависимости для определения размера зоны наведенной трещиноватости RH.T.3 относительно допустимой скорости смещения. Полученные зависимости представлены в табл. 3.
На рис. 6 приведены размеры зон наведенной трещиноватости по трем областям, рассчитанные по вышеприведенной программе (см. рис. 2).
Если производится поскважинное взрывание, то в формулы для расчета зоны наведенной трещиноватости подставляется масса заряда ВВ в крайней от законтурного массива скважине.
При мгновенном взрывании группы зарядов расчет зоны наведенной трещиноватости рассчитывается по массе эквивалентного заряда, сейсмическое действие которого эквивалентно действию группы зарядов. Массу эквивалентного заряда можно определить по формуле 3 ,
где n – общее количество зарядов; Qi – масса i-гo заряда; rb – расстояние от охраняемого объекта до ближайшего заряда, м; ri – расстояние от охраняемого объекта до i-гo заряда, м.
На представленных расчетах размеров зон (рис. 7) видно, что наибольший размер зоны наведенной трещиноватости приурочен к верхней части карьера. В средней и нижней частях карьера зоны наведенной трещиноватости находятся примерно на одном уровне.
Более существенное разуплотнение прибортового массива в верхней части карьера связано с различной блочностью массива по глубине карьера. ОАО «ВИОГЕМ» [16] исследована динамика развития трещиноватости в плане и по глубине породного массива в карьере рудника «Железный» и показано, что средний размер отдельности пород увеличивается с глубиной карьера. Сначала резко (от 0,22 до 1,19 м) в интервале 0–150 м от дневной поверхности, затем постепенно (от 1,19 до 1,4 м) до глубины 450 м, а глубже стабилизируется на уровне 1,4–1,5 м, что и дает различие в размерах зон наведенной трещиноватости. Размер зоны наведенной трещиноватости также хорошо коррелирует с данными математического моделирования (см. рис. 7).
Заключение
Для определения ширины зоны трещинообразования и заколов в законтурном массиве применен методический подход М.И. Шуйфера и А.Е. Азарковича.
С целью определения сейсмического воздействия массовых взрывов на непосредственно прилегающий к инициируемому блоку массив горных пород были проведены замеры скорости смещения грунта в ближней зоне действия взрыва на различных участках карьера «Железный» АО «Ковдорский ГОК». По степени реакции массива на динамические воздействия произведено районирование карьера по глубине. Выделены три области сейсмичности, для которых получены эмпирические зависимости изменения скорости смещения от приведенного расстояния и соответствующие коэффициенты сейсмичности и затухания. По ним рассчитаны размеры зоны наведенной трещиноватости относительно допустимой скорости смещения.
Наибольший размер зоны наведенной трещиноватости приурочен к верхней части карьера. В средней и нижней частях размеры зон наведенной трещиноватости находятся примерно на одном уровне. Большее разуплотнение прибортового массива в верхней части карьера связано с различной блочностью массива – средний размер отдельности пород увеличивается с глубиной карьера.
Разработана инженерная методика оперативной оценки сейсмовзрывного воздействия на законтурный массив при производстве взрывных работ в карьере для расчета ширины зон заколов, трещинообразования в тыльную часть массива и зон проявления наведенной трещиноватости. Расчет на стадии проектирования позволяет подобрать параметры взрывания, обеспечивающие минимизацию сейсмического воздействия на законтурный массив.
Для оперативного расчета положения границ разрушения в тыльную часть массива при взрыве скважин разработана компьютерная программа. Полученные расчетные значения хорошо согласуются с ранее полученными данными моделирования и практики.
Список литературы
1. Шемякин Е.И. Сейсмовзрывные волны в процессе горного производства. М.: ННЦ ГП – ИГД им. А. А. Скочинского; 2004. 75 с.
2. Адушкин В.В., Спивак А.А. Геомеханика крупномасштабных взрывов. М.: Недра; 1993. 319 с.
3. Боровиков В.А., Ванягин И.Ф. Моделирование действия взрыва при разрушении горных пород. М.: Недра; 1990. 231 с.
4. Крюков Г.М. Феноменологическая квазистатическо-волновая теория деформирования и разрушения материалов взрывом зарядов промышленных ВВ. М.: Изд-во Моск. гос. гор. ун-та, 2003. 65 с.
5. Шер Е.Н., Черников А.Г. Расчет параметров радиальной системы трещин, образующейся при взрыве удлиненного заряда в хрупких горных породах. Фундаментальные и прикладные вопросы горных наук. 2015;2(2):299–303. Sher E.N., Chernikov A.G. Calculating parameters of radial crack system induced by elongated charge blasting in brittle rocks. Fundamentalnye i Prikladnye Voprosy Gornykh Nauk. 2015;2(2):299–303. (In Russ.)
6. Менжулин М.Г., Федосеев А.В. Определение размеров зон взрывного разрушения для трещиноватых горных пород с различными заполнителями трещин на примере Михайловского ГОКа. Записки Горного института. 2012;195:120– 123. Режим доступа: https://pmi.spmi.ru/pmi/article/view/6116 (дата обращения: 14.04.2024). Menchzhulin M.G., Fedoseev A.V. Amount of explosive destruction zones for cracked rocks with different crack filler for example Mikhailovsky GOK. Journal of Mining Institute. 2012;195:120–123. (In Russ.) Available at: https://pmi.spmi.ru/pmi/article/view/6116 (accessed: 14.04.2024).
7. Казаков Н.Н. Взрывная отбойка руд скважинными зарядами. М.: Недра; 1975. 192 с.
8. Боровиков В.А., Ванягин И.Ф. К расчету параметров волны напряжения при взрыве удлиненного заряда в горных породах. Взрывное дело. 1976;(76/33):74–85. Borovikov V.A., Vanyagin I.F. К расчету параметров волны напряжения при взрыве удлиненного заряда в горных породах. Explosion Technology. 1976;(76/33):74–85. (In Russ.)
9. Бернард Т. Новая модель для реалистичного моделирования взрыва в 3D. Глобус. Геология и бизнес. 2019;(4):112–127. Режим доступа: https://www.vnedra.ru/tehnologii/informacionnye-tekhnologii/novaya-model-dlya-realistichnogomodelirovaniyavzryva-v-3d-8604/ (дата обращения: 14.04.2024). Bernard T. A new model for realistic blast simulation in 3D. Globus. Geologiya i Biznes. 2019;(4):112–127. (In Russ.) Available at: https://www.vnedra.ru/tehnologii/informacionnye-tekhnologii/novaya-model-dlya-realistichnogo-modelirovaniya-vzryvav-3d-8604/ (accessed: 14.04.2024).
10. Holmberg R., Persson P.-A. The Swedish approach to contour blasting – proceedings of the 4th conference on explosives and blasting technique. In: Konya, C.J. (ed.), Society of Explosives Engineers, New Orleans, Lousiana, February 10–15, 1978, pp. 113–127.
11. Tannant D.D., Peterson J. Evolution of blasting practices at the EkatiTM diamond mine. In: 17th International Mining Congress and Exhibition of Turkey. 2001, pp. 297–304. Available at: https://api.maden.org.tr/uploads/portal/resimler/ekler/9eca5979ccbb752_ek.pdf (accessed: 14.04.2024)
12. Шуйфер М.И., Азаркович А.Е. Расчет размеров зоны трещинообразования при взрыве скважинных зарядов в скальном массиве. Взрывное дело. 1980;(82/39):191–209. Shuifer M.I., Azarkovich A.E. Calculation of the fracture zone sizes due to blasting of borehole charges in the rock mass. Explosion Technology. 1980;(82/39):191–209. (In Russ.)
13. Ляховицкий Ф.М. О соотношении упругих и прочностных свойств горных пород. В кн.: Геофизические исследования. М.: Изд-во МГУ, 1964. Вып. 1. С. 294–305.
14. Козырев С.А., Камянский В.Н., Аленичев И.А. Оценка взаимодействия скважинных зарядов при различных интервалах замедлений между ними. Взрывное дело. 2017;(117-74):60–75. Kozyrev S.A., Kamyansky V.N., Alenichev I.A. Estimation of interaction borehole charges at different delay time between them. Explosion Technology. 2017;(117-74):60–75. (In Russ.)
15. Козырев С.А., Камянский В.Н. Оценка влияния сейсмовзрывных нагрузок в ближней зоне взрыва. Горный информационно-аналитический бюллетень. 2017;(S23):316–324. https://doi.org/10.25018/0236-1493-2017-10-23-316-324 Kozyrev S.A., Kamyansky V.N. Estimation of the effect of seismic loads at the near zone of explosion. Mining Informational and Analytical Bulletin. 2017;(S23):316–324. (In Russ.) https://doi.org/10.25018/0236-1493-2017-10-23-316-324
16. Игнатенко И.М., Яницкий Е.Б., Дунаев В.А., Кабелко С.Г. Трещиноватость породного массива в карьере рудника «Железный» АО «Ковдорский ГОК». Горный журнал. 2019;(10):11–15. https://doi.org/10.17580/gzh.2019.10.01 Ignatenko I.M., Yanitsky E.B., Dunaev V.A., Kabelko S.G. Jointing of rock mass in open pit at the Zhelezny mine of the Kovdor Mining and Processing Plant. Gornyi Zhurnal. 2019;(10):11–15. (In Russ.) https://doi.org/10.17580/gzh.2019.10.01