CC BY
31
УДК 504.55.054:622(470.6)
ГЕОМЕХАНИЧЕСКИЕ ПРЕДПОСЫЛКИ СОХРАНЕНИЯ УСТОЙЧИВОСТИ ВЫРАБОТОК ПРИ РАЗРАБОТКЕ ВОДООБИЛЬНЫХ МЕСТОРОЖДЕНИЙ
Ю.В. Дмитрак, В.И. Голик, В.В. Венигор
Изложены результаты моделирования поведения водообильных массивов. Подтверждена коррелятивная зависимость между параметрами водонасыщения и напряжениями в структурных породных блоках и устойчивостью пройденных в них выработок. Показано, что устойчивость выработок увеличивается при уменьшении продолжительности контактов пород с влагой, а также, что при добыче руд вариантами с восходящей выемкой руд и закладкой пустот твердеющими смесями вероятность прорыва вод и ослабления устойчивости выработок уменьшается. Приведены количественные параметры указанной зависимости.
Ключевые слова: разработка, месторождение, водообильность, подземный способ, прорыв вод, геомеханика, массив, структурный блок, порода.
В эпоху научно-технического прогресса вопросы защиты окружающей среды от негативного влияния горных технологий стали государственным делом и частью межгосударственных отношений. Экологические показатели горного производства не улучшились, а воздействие на окружающую среду превышает восстановительные возможности природы [1].
Воздействие горного дела на литосферу характеризуется рамками земельных отводов и глубиной горных разработок. Горные технологии являются инструментом воздействия на окружающую среду, поэтому их совершенствование является приоритетным направлением инженерной мысли современности.
Ухудшение условий разработки является установившейся тенденцией горного производства [2]. Один из опасных элементов - разработка месторождений с невысокими коэффициентами фильтрации с прорывами воды в производственную зону. Даже при достаточной гидрогеологической информации с высокой достоверностью прогнозировать опасные участки массива трудно [3].
На Северо-Кавказских рудниках такие прорывы создавали опасность для работающих и недр. На Квайсинском руднике разжиженная горная масса прорывалась на глубину 1300 м. На Згидском и Фиагдонском рудниках прорывы насыщенных водой пород в очистные выработки происходили нередко.
На Сибирских железных рудниках в результате горных работ объемы прорвавшейся породы достигали 2000 м . На руднике им. 40-летия ВЛКСМ Лениногорского полиметаллического комбината обводненные породы распространялись по выработке со скоростью 1,5 м/с.
На руднике «Муфулира» над обрушенной зоной 20 млн тонн отвальных хвостов в полужидком состоянии прорвались в выработку менее чем за 15 мин, затопили половину рудника, унеся жизнь 89 человек. Аварийные прорывы обводненных пород определяют существенную долю потерь руды (табл. 1).
Таблица 1
Потери руды на Сибирских рудниках, тыс. т/год
Рудоуправление Всего потери При выпуске От аварий
Темир-Тауское 335,5 115,0 220,5
Шерегешское 923,3 297,8 625,5
Ташгольское 704,4 237,5 466,9
В результате освоения Михайловского месторождения КМА вокруг основных водозаборных узлов (г. Курск, г. Железногорск) сформировались региональные депрессионные воронки, определяющие дефицит питьевой воды нормативного качества.
Сложными инженерно-геологическими и гидрогеологическими условиями характеризуется Яковлевское месторождений железных руд, разрабатываемое под неосушенным высоконапорным водоносным горизонтом с оставлением водозащитного целика, сложенного водонеустойчи-выми рудами с потенциальной опасностью прорывов подземных вод для обеспечения безопасного ведения горных работ. Вскоре ту же задачу предстоит решать и на Коробковском месторождении КМА при освоении запасов следующего горизонта месторождения.
В Подмосковном бассейне из-за обводненности затрудняется эксплуатация месторождений. На 1 т добываемого угля откачивают около 10
3
м воды.
С прорывами воды в выработки связано много аварий, самая крупная из них произошла в 1880 г. в польской шахте «Renard» унесла жизни 200 человек.
В Подмосковном бассейне при освоении Никулинского и Бельцевского месторождений с мощным слоем обводненных известняков участились прорывы воды. Известные способы осушения и управления массивами не обеспечили безопасной и эффективной работы шахт. Высокая обводненность явилась одним из решающих факторов при ликвидации шахт Подмосковного бассейна.
Гидрогеологические осложнения оказывают комплексное воздействие на все процессы горного производства. Выбросы пород в выработки
вызывают необходимость капитального ремонта оборудования, эксплуатация электроустановок требует дополнительных мер по их защите и т.п.
Профилактика гидрогеологических осложнений при подработке водных объектов сводится к понижению уровня вод [4].
Широкое распространение получил способ понижения уровня подземных вод путем создания системы дренажных устройств, которая позволяет улавливать подземные воды дренажной системой и отводить их за пределы зоны. Ограждение горных выработок заключается в создании вокруг объекта искусственной преграды, оттеснение вод - в создании по контуру объекта или в его центре системы специальных устройств, вскрывающих водоносный горизонт скважиной или стволом. В водоносный горизонт, участвующий в обводнении, нагнетается сжатый воздух, который оттесняет воду из зоны ведения добычных работ.
При разработке месторождений со сложной гидрогеологией осушительные мероприятия не всегда увязывают в единый технологический цикл с добычей полезного ископаемого.
Профилактика обводнения эффективна не всегда, поэтому важное место уделяется геомеханическим факторам управления водонасыщенны-ми массивами.
Район Квайсинского рудного поля сложен юрскими, меловыми, палеогеновыми и четвертичными отложениями. В основании стратиграфического разреза залегает сорская свита, ядро антиклинали представлено глинистыми сланцами и слюдистыми песчаниками [5].
Квайсинское месторождение представлено эпитермальными рудами, локализованными в трещинах, жилах и межпластовых залежах. Мощность рудной зоны 0,4...6,0 м. Рудные тела имеют сложные, прерывистые формы, образующие в ряде случаев систему линзообразных тел с неравномерным распределением рудной массы.
Квайсинское месторождение относится к классу месторождений со сложными геологическими условиями в связи с крупной тектонической зоной и наличием карстовых полостей, вскрываемых в зоне аэрации на малой глубине (25.30 м) от поверхности, где они имеют форму вертикальных каналов переменного сечения. Подземные воды содержатся в четвертичных, палеогеновых, меловых и юрских образованиях, но главная роль принадлежит трещино-карстовым и трещинным подземных водам.
Разгрузка подземных вод осуществляется в подземные выработки рудника в виде отдельных и групповых водопроявлений и обводненных участков, приуроченных к зонам тектонических нарушений.
Для эксплуатации месторождений с высокой обводненностью необходима гидрогеологическая информация, которая может быть получена прямыми измерениями дебета воды на замерных станциях в выработках. На Квайсинском месторождении наиболее представительными оказались исследования водопритока в штольне "Капитальная", где сбрасывались на
поверхность рудничные воды, формирующиеся во всем месторождении.
3 3
Водоприток изменялся от 100 м /ч в январе до 375 м /ч в апреле - мае.
При разработке месторождения системами с обрушением руды и пород водообильность интенсифицирует развитие процесса сдвижения массива горных пород. В свою очередь, деформирование пород изменяет их проницаемость. Выработанное пространство может быть интерпретировано в виде полости шириной до 20 м, длиной до 700 м и глубиной до 700 м, заполненной дискретными минералами.
Деформации горного массива выходят на поверхность и являются причиной водонасыщенности.
По данным выполненных экспериментальных исследований модуль трещиноватости пород изменяется в среднем от 2...7 при влажности 5 % в сухих забоях до 27...31 в местах водопроявления. При влажности 20 % модуль трещиноватости в сравнении с сухими забоями увеличивается в 5 - 6 раз. Полученные данные позволяют конкретизировать зависимость модуля трещиноватости от влажности [6].
Раскрытие трещин при водонасыщении увеличивает дифференцирование массива на дискретные блоки, снижая его устойчивость. Поэтому важны принципы обеспечения устойчивости выработок на участках с различной обводненностью и структурной ослабленностью.
Величина пролета обнажений может быть найдена графическим построением в зависимости от угла падения рудного тела а.
^ 90° + ф 2к т
При а >-- Ь =-+ •
2 90° +ф sin а
т . 90°
ПАо Mtga + tg--)
90 + ф т v ь ь о ' m
при а <-- L =-2-+
2 90° + ф sin а
tga' tg—^
где h - высота выработанного пространства; ф - угол внутреннего трения пород; m - мощность рудного тела.
На величину угла обрушения влияют внешние факторы, в том числе водонасыщенность. Величина угла и зоны обрушения в результате горных работ на Квайсинском месторождении определена съемкой поверхности над выработанным пространством. Полученный фактический угол обрушения сравнивается с расчетными данными (табл. 1).
Комплекс исследований гидрогеологических условий разработки включает в себя определение гидрогеомеханической позиции массива пород: внешнего контура подработки водных объектов, границ и размеров зон водопроводящих трещин, обрушения и сдвига пород, а также фильтрационных процессов [7].
Таблица 1
Угол обрушения пород по альтернативным данным
Автор Угол обрушения
ВНИМИ о 85
Голик В.И. о 84
Казикаев Д.М. о 80
Фактически измеренный о 69
Наиболее распространенным способом определения прочностных
свойств является срез со сжатием. Исследования проведены в матрицах,
позволявших изменять величину нормальных и касательных усилий в
плоскости среза под углами 30, 45 и 60°. Для каждого угла среза найдена
величина нормальных и касательных напряжений в момент среза:
P . P
а =— sin а; т = — cos а, S' S'
где Р - предельная сжимающая нагрузка в момент разрушения образца; S - площадь среза; а - угол наклона плоскости среза к направлению внешнего усилия.
Исследованиями порфиритов лежачего бока рудного тела и известняков висячего бока отмечено уменьшение величины критических напряжений от 0 до 20 % при изменении влажности на 25.30 % вследствие проникновения воды в породу с ослаблением связей между структурными элементами породы, в результате чего прочность пород изменяется (табл. 2, 3).
Таблица 2
Изменение свойств порфиритов в зависимости от влажности
Влаж- Угол Нормальное Касательное Угол внутрен- Сила сцепле-
ность об- наклона напряжение, напряжение т, него трения, ния , МПа
разца, % матрицы, град МПа МПа град
0 30 19,0 32,8
0 45 37,2 37,2 30 23,0
0 60 67,4 38,9
10 30 15,2 26,3
10 45 31,0 31,0 25 18,0
10 60 58,3 33,7
20 30 12,1 20,9
20 45 26,0 26,0 22 15,0
20 60 47,2 27,3
Гипотезы горного давления исходят из положения о перераспределении первоначальных напряжений в массиве в результате техногенного вмешательства [8 - 11].
Скальным массивам больше других отвечает гипотеза горного давления М.М. Протодьяконова. Над выработанным пространством образуется свод естественного равновесия, а давит на кровлю выработки вес пород в пределах свода (рис. 1).
Таблица 3
Результаты измерения основных прочностных свойств известняков
в зависимости от влажности
Влажность Угол накло- Нормальное Касательное Угол внут- Сила сцеп-
образца, % на матрицы, напряжение, напряжение т, реннего ления , МПа
град МПа МПа трения, град
0 30 11,2 19,4
0 45 21,1 21,1 29 12,0
0 60 41,2 23,8
10 30 8,0 13,8
10 45 16,7 16,7 23 9,0
10 60 33,3 19,3
20 30 7,2 12,4
20 45 14,0 13,0 21 8,0
20 60 29,2 16,9
Рис. 1. Условие создания плоской кровли: 1 - граница свода естественного равновесия пород; 2 - отслоившиеся породы;
3 - несущий слой заклинившихся структурных блоков
По С.В. Ветрову реальные массивы состоят из отдельных блоков, ограниченных нарушениями. Они образуют устойчивые системы при самозаклинивании [12]. Поверхности свода естественного равновесия оказываются образованными гранями заклинивающихся и плотно прижатых друг к другу отдельных блоков.
Свод с максимальным пролетом и высотой, который может возникнуть при данной прочности пород и размерах структурных блоков,
Ь = 3
23
10Я
к 2 '
где ^ - горизонтальный размер структурных блоков, м; Я'- временное сопротивление пород одноосному сжатию, кг/см ; Н - глубина расположения основания свода, м; у - объемный вес пород, т/м ; к1 - коэффициент запаса устойчивости свода.
Исследованные гипотезы устойчивости выработок не учитывают сложной структуры месторождений, поэтому авторами предложено выражение для определения предельных пролетов обнажений для пород основной и непосредственной кровли, учитывающее место нахождения выработки в массиве относительно рудовмещающих структур. Допустимый по условию сохранения плоской формы пролет обнажений основной кровли [13 - 15]
Ь1 = 1,713
10Rd о22КО
У
где Ясж - предел прочности на сжатие, кг/см ; й02 - вертикальный размер структурного блока основной кровли, м; У0 - коэффициент устойчивости пород основной кровли; к2 - коэффициент запаса устойчивости плоского обнажения, 2,5 - 3,0; у - объемный вес пород, т/м .
Для сложноструктурных месторождений гидротермального типа характерно оруденение в середине части шва. Непосредственная кровля ведет себя как отслаивающаяся пачка - в ее структурных блоках не образуется замкнутого свода. Для непосредственной кровли устойчивый пролет
Ь = \Шп 2
10Я_,
^2 Уь
где - вертикальный размер структурных блоков пород непосредственной кровли, м; Ь - мощность отслаивающейся пачки непосредственной кровли, м.
Гипотеза С.В. Ветрова более точно воспроизводит процессы образования устойчивых пролетов, однако для месторождений со сложной гидрогеологией необходимо учесть влажность пород.
Взаимодействие воды с породой ослабляет структурные связи пород и уменьшает устойчивость массивов. В породах с коагуляционными и конденсационными связями происходят набухание и размокание, а с кристаллизационными связями снижение механической прочности. Движущиеся потоки воды, проникая по трещинам, способствуют их раскрытию, разрушая связи и вынося их потокам. Движущиеся потоки воды изменяют механическую прочность и устойчивость при воздействии фильтрационного потока. Поэтому целесообразно введение коэффициента ослабления, который бы учитывал
Ксл = /
А(М;); £»(„); £*.(„)
А Ск Фк
где Ь / 1с - отношение ширины или высоты выработки к размеру структурного элемента пород; свсж / стсж - отношение предела прочности на одноосное сжатие породы во влажном и сухом состоянии; Ст/ Ск - отношение удельных сцеплений в монолите породы и на контактах структурных элементов; с / срк - отношение угла наклона поверхности ослабления к направлению наибольшего нормального напряжения и угла трения на контакте элементов; - влажность пород,
Геомеханическая устойчивость зависит от взаимодействия полей напряжений. В исходном состоянии (нетронутый массив) напряжения оцениваются С1 + к с 23, в котором коэффициент к позволяет учесть фактор
сложноструктурности. Условие геомеханической сбалансированности массива описывается моделью
где с1 - вертикальная составляющая главных напряжений; с23 - горизонтальная результирующая составляющая главных напряжений; к - коэффициент, зависящий от степени искажения напряжений структурно -тектоническими условиями; ссж - напряжение в зоне влияния выработок; х1 - физико-механические свойства пород, х2 - рудного тела; х3 - влияние обводненности.
Потеря несущей способности слоя заклинившихся пород в основании свода естественного равновесия происходит при скалывании наиболее напряженных участков структурных блоков (рис. 2).
Рис. 2. Схема разрушения прочности слоя заклинившихся блоков
при скалывании пород в замке
Вода, омывая поверхность структурных блоков и проникая внутрь их по микротрещинам, уменьшает прочность пород и ускоряет разрушение породной конструкции, провоцируя разгрузку напряжений в своде при обрушении пород в выработку [16 - 18].
Управляющее воздействие на прочность системы заключается в использовании параметров, при которых жесткость перекрытия обеспечивается за счет остаточной прочности пород и напряжений бокового распора, которые могут быть увеличены укреплением пород.
Управление напряжениями в нижнем слое нарушенных пород состоит в реализации остаточной несущей способности пород ограничением
Косл = Жосл (.
С + кс23 < сс ж= /х(0X1; 0х2; dx3;...; dxn),
с Ж
величины пролетов раскройкой шахтного поля на геомеханически сбалансированные участки. В пределах геомеханически сбалансированных участков могут быть применены минимизированные по затратам труда и материалов технологии.
Исследование закономерностей геомеханической системы и поведения обводненных рудовмещающих массивов позволяет разработать геотехнологические решения, способные управлять состоянием водообиль-ных массивов в оптимальном режиме.
При разработке водообильных объектов опасна связь между очистным пространством и водным объектом [19]. Сущность предлагаемой технологии заключается в том, что отработка месторождения ведется в направлении снизу вверх системой с закладкой выработанного пространства твердеющей смесью. Выемка полезного ископаемого от воды осуществляется под защитой нетронутого рудного массива над фронтом очистных работ. Коэффициент фильтрации пород в целиковой части значительно ниже, чем в нарушенном сдвижением массиве, поэтому поверхностные и атмосферные воды не попадают на новые рабочие горизонты и не вызывают гидрогеологических осложнений.
При этом в начальный период освоения нового горизонта увеличение притока за счет срабатывания статических запасов нижних горизонтов возможно, но по мере проведения выработок приток воды стабилизируется.
Выемка полезного ископаемого производится горизонтальными слоями снизу вверх с одновременной закладкой выработанного пространства твердеющими смесями. Выемочный блок имеет длину до 60 м, ширина его ограничивается мощностью рудного тела, а высота блока равна высоте этажа.
Рекомендуемая технология разработки позволяет:
- управлять состоянием водообильных скальных массивов путем решения геомеханической задачи;
- минимизировать потери и разубоживание руды, поскольку системы с закладкой выработанного пространства являются наиболее ресурсосберегающими;
- сохранить земную поверхность и минимизировать ущерб окружающей среды;
- повысить производительность и безопасность работ.
Экономический ущерб на ликвидацию гидрогеологических осложнений, помимо затрат непосредственно на осушение месторождений,
У = У + З
У ликв.го У л 1 Зосуш. м 5
где Ушкв.го - ущерб на ликвидацию осложнений; Зосуш. м - затраты на осушительные мероприятия.
Ущерб на ликвидацию последствий выбросов водонасыщенных пород
Ц =УС, 226
где V - объем выбросов, м3; C - себестоимость ликвидации 1 м3выбросов, руб.
Экономический ущерб в социальной сфере складывается от потери трудоспособности при заболеваемости рабочих из-за плохих условий труда; ущерба от выплат из фонда социального страхования; ущерба от лечения потерявших трудоспособность и т.п.
Список литературы
1. Ivashchuk O., Ivashchuk Or. Automation and Intellectualization to Control the Ecological Situation in the Urbanized Territories // Proceedings of the 7th International Conference on Intelligent Data Acquisition and Advanced Computing Systems (IDAACS). Berlin, Germany, 2013. Vol. 2. Р. 814-820.
2. Голик В.И. Природоохранные технологии разработки рудных месторождений: учеб. пособие для вузов. М.: Инфра, 2014. 190 с.
3.Wu Qiang, Pang Wei, Dai Yingchun Vulnerability forecasting model based on coupling technique of GIS and ANN in floor groundwater inrush // Journal of China Coal Society. 2006. 31(3). Р. 314-319.
4.Wu Qiang, Zhang Zhilong, Zhang Shengyuan A new practical methodology of the coal floor water inrush evaluating // Journal of China Coal Society.2007. 32(11). Р. 1 121-1 126.
5. Габисов В.Ю. Определение путей проникновения подземных вод на отрабатываемые горизонты Квайсинского рудника // Сб. науч. тр. науч.-техн. конф. СКГМИ. Владикавказ, 1993. С. 69 - 70.
6. Голик В.И., Габисов В.Ю. Влияние природных и техногенных условий на сохранность окружающей среды при добыче полезных ископаемых // Изв. вузов. Горный журнал. 1995. № 7.
7. Дубинин Н.Г., Хромцов В.Ф., Шеховцов В.С. Предотвращение прорывов глинистых пород при разработке рудных месторождений // Новосибирск, 1989. 118 с.
8. Казикаев Д.М., Осипенко Ю.С. Разработка рудных месторождений под водными объектами. М.: Недра, 1989. 192 с.
9. Ржевский В.В., Новик Г.Я. Основы физики горных пород. М.: Недра, 1984. 359 с.
10. Борисов А.А. Механика горных пород и массивов. М.: Недра, 1980. 246 c.
11. Gattinoni P., Pizzarotti E. M., Scesi L. Engineering Geology for Underground Works // Springer, 2014. 312 p.
12. Ветров С.В. Допустимые размеры обнажений горных пород при подземной разработке руд. М.: Наука, 1975. 198 с.
13. Golik V., Komashchenko V., Morkun V. Geomechanical terms of use of the mill tailings for preparation // Metallurgical and Mining Industry. 2015. №. 4. Р. 321-324.
14. The theory and practice of rock massifs control in the ore mining / V.Golik, V.Komaschenko,V. Morkun,T.N. Gvozdkova // Metallurgical and Mining Industry. 2016. № 1. Р. 209-212.
15. Modelling of rock massifs tension at underground ore mining / V.Golik, V.Komaschenko, V.Morkun, O.Burdzieva // Metallurgical and Mining Industry. 2015. № 8. Р. 540-543
16. Griffith A. A. The theory of rupture // Proc. Ict. Int. Congr. Appl. Mech. Delft. 1924. P. 55-63.
17. Dold B., Weibel L. Biogeometallurgical pre-mining characterization of ore deposits: An approach to increase sustainability in the mining process // Environmental Science and Pollution Research. 2013. Vol. 20. No. 11. P. 7777 -7786.
18. Hencher S. Practical Engineering Geology. 2012. 450 p.
19. Бабец А. М., Лейзерович С. Г., Тиганова Ю. В. Вклад ОАО «НИИКМА» в решение проблем Курской магнитной аномалии // Горный журнал. 2014. №8. С.67-72.
Дмитрак Юрий Витальевич, д-р техн. наук, проф., ректор, v.i.golik a mail.ru, Россия, Владикавказ, Северо-Кавказский государственный технологический университет,
Голик Владимир Иванович, д-р техн. наук, проф., v. i.golikamail.ru, Россия, Владикавказ, Северо-Кавказский государственный технологический университет,
Вернигор Владимир Владимирович, асп., [email protected], Россия, Москва,
МИСиС
GEOMECHANICAL CONSERVATION PRINCIPLES SUSTAINABILITY OF DEVELOPMENTS AT MINING WATER DEPOSITS
Yu. V. Dmitrak, V.I. Golik, V. V. Venigor
The results of modeling the behavior of water bodies are described. The correlation between water saturation parameters and stresses in structural rock blocks and the stability of excavations passed in them was confirmed. It is shown that the stability of excavations increases with decreasing duration of contact of rocks with moisture. It is shown that when mining ores with variants with ascending excavation of ores and laying of voids with hardening mixtures, the probability of breakthrough of water and weakening of the stability of the excavations decreases and the quantitative parameters of this dependence are given.
Key words: development, deposit, water supply, underground method, water breakthrough, geomechanics, massif, structural block, rock.
Dmitrak Yuryi Vitalievich, doctor of technical sciences, professor, the Rector, v.i.golika mail. ru, Russia, Vladikavkaz, North-Caucasian State Technological University,
Golik Vladimir Ivanovich, doctor of technical sciences, professor, v. i. golikamail. ru, Russia, Vladikavkaz, North-Caucasian State Technological University,
Venigor Vladimir Vladimirovich, postgraduate, vovrsoros@yu. ru, Russia, Moscow,
MISiS
Reference
1. Ivashchuk O., Ivashchuk Or. Automation and Intellectualization to Control the Ecological Situation in the Urbanized Territories // Proceedings of the 7th International Conference on Intelligent Data Acquisition and Advanced Computing Systems (IDAACS). Berlin, Germany, 2013. Vol. 2. R. 814-820.
2. Golik V.I. Prirodoohrannye tehnologii razrabotki rudnyh mestorozhde-nij: ucheb. posobie dlja vuzov. M.: Infra. 2014. 190 s.
3. Wu Qiang, Pang Wei, Dai Yingchun Vulnerability forecasting model based on coupling technique of GIS and ANN in floor groundwater inrush // Journal of China Coal Society. 2006. 31(3). R. 314-319.
4. Wu Qiang, Zhang Zhilong, Zhang Shengyuan A new practical methodology of the coal floor water inrush evaluating // II - the vulnerable index method. Journal of China Coal Society.2007. 32(11). R. 1 121-1 126.
5. Gabisov V.Ju. Opredelenie putej proniknovenija podzemnyh vod na otraba-tyvaemye gorizonty Kvajsinskogo rudnika // Sb. nauch. tr. nauch.-tehn. konf. SKGMI/ Vladikavkaz. 1993. S. 69 - 70.
6. Golik V.I., Gabisov V.Ju. Vlijanie prirodnyh i tehnogennyh uslovij na sohrannost' okruzhajushhej sredy pri dobyche poleznyh iskopaemyh // Gornyj zhurnal. Izv. Vuzov, 1995. № 7.
7. Dubinin N.G., Hromcov V.F., Shehovcov V.S. Predotvrashhenie proryvov glin-istyh porod pri razrabotke rudnyh mestorozhdenij // Novosibirsk, 1989. 118 s.
8. Kazikaev D.M., Osipenko Ju.S. Razrabotka rudnyh mestorozhdenij pod vod-nymi ob#ektami. M.: Nedra, 1989. 192 s.
9. Rzhevskij V.V., Novik G.Ja. Osnovy fiziki gornyh porod. M.: Nedra, 1984. 359 s.
10. Borisov A.A. Mehanika gornyh porod i massivov. M.: Nedra, 1980. 246 c.
11. Gattinoni P., Pizzarotti E. M., Scesi L. Engineering Geology for Underground Works// Springer, 2014. 312 p.
12. Vetrov S.V. Dopustimye razmery obnazhenij gornyh porod pri podzemnoj razrabotke rud. M.: Nauka, 1975. 198 s.
13. Golik V., Komashchenko V., Morkun V. Geomechanical terms of use of the mill tailings for preparation // Metallurgical and Mining Industry, 2015. №. 4. R. 321-324.
14. The theory and practice of rock massifs control in the ore mining / V.Golik, V.Komaschenko,V. Morkun,T.N. Gvozdkova // Metallurgical and Mining Industry, 2016. № 1. R. 209-212.
15. Modelling of rock massifs tension at underground ore mining / V.Golik, V.Komaschenko, V.Morkun, O.Burdzieva // Metallurgical and Mining Industry, 2015. № 8. R. 540-543
16. Griffith A. A. The theory of rupture // Proc. Ict. Int. Congr. Appl. Mech. Delft. 1924. P. 55-63.
17. Dold B., Weibel L. Biogeometallurgical pre-mining characterization of ore deposits: An approach to increase sustainability in the mining process // Environmental Science and Pollution Research. 2013. Vol. 20. No. 11. P. 7777-7786.
18. Hencher S. Practical Engineering Geology // Spon Press, 2012. 450 p.
19. Babec A. M., Lejzerovich S. G., Tiganova Ju. V. Vklad OAO «NIIKMA» v res h-enie problem Kurskoj magnitnoj anomalii // Gornyj zhurnal, 2014. №8. S.67-72.
