Прогноз сдвижений и деформаций подрабатываемых скальных массивов рудных месторождений с учетом их структурно-тектонических особенностей тема диссертации и автореферата по ВАК РФ 25.00.20, кандидат технических наук Шустов, Денис Владимирович

• Шустов, Денис Владимирович
• кандидат технических наук
• 2012

Оглавление диссертации кандидат технических наук Шустов, Денис Владимирович

1.СОСТОЯНИЕ ВОПРОСА. ЦЕЛИ И ЗАДАЧИ ИССЛЕДОВАНИЯ.

1.1 Методические подходы к прогнозу параметров процесса сдвижения на рудных месторождениях.

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

1.3 Теоретические методы расчета сдвижений и деформаций.

1.4 Численные методы расчета сдвижений и деформацийподрабатываемых массивов.

1.5. Цели и задачи исследования.

2. ИНСТРУМЕНТАЛЬНЫЕ НАБЛЮДЕНИЯ ЗА СДВИЖЕНИЕМГОРНЫХ ПОРОД НА ТИШИНСКОМ И РИДДЕР-СОКОЛЬНОМ

МЕСТОРОЖДЕНИЯХ.

2.1. Сдвижение горных пород и земной поверхности при отработке Тишинского месторождения.

2.1.1. Краткая геологическая характеристика Тишинского месторождения.

2.1.2. Исходное напряженное состояние массива горных породТишинского месторождения.

2.1.3. Состояние наблюдательной станции в породах южного фланга месторождения.

2.1.4. Характер развития деформационных процессов по данным инструментальных наблюдений 1988-2009гг.

2.1.6. Развитие процессов сдвижения на северо-восточном фланге месторождения.

2.1.7. Развитие процессов сдвижения на западном флангеместорождения

2.1.8. Развитие концентраций горизонтальных деформаций на месторождении.

2.2. Сдвижение горных пород и земной поверхности в районе шахт №3 и "Вентиляционная".

2.2.¡.Краткая геологическая характеристика Риддер-Сокольного месторождения.

2.2.2. Развитие процессов сдвижения на месторождении.

ВЫВОДЫ.

3 ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ВЛИЯНИЯ СТРУКТУРНО-ТЕКТОНИЧЕСКИХ ОСОБЕННОСТЕЙ СКАЛЬНОГО МАССИВА НА ПАРАМЕТРЫ ПРОЦЕССА СДВИЖЕНИЯ.

3.1 Аналитическая часть модели. Выбор критерия разрушения.

3.2 Реализация модели в методе дискретных элементов.

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

ВЫВОДЫ.

4.1. Расчеты сдвижений и деформаций массива горных пород на южном фланге Тишинского свинцово-цинкового месторождения.

4.1.1. Моделирование процесса сдвижения на современном этапеотработки месторождения.

4.1.2. Прогноз параметров процесса сдвижения горных пород приотработке запасов до 20-го горизонта.

4.2. Моделирование развития процессов сдвижения горных пород в районе стволов шахт №3 и «Вентиляционная» Риддер-Сокольного месторождения

ВЫВОДЫ.

Введение диссертации (часть автореферата) на тему «Прогноз сдвижений и деформаций подрабатываемых скальных массивов рудных месторождений с учетом их структурно-тектонических особенностей»

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

Многолетние инструментальные наблюдения и теоретические исследования, проведенные в ведущих научно-исследовательских институтах, таких как ВНИМИ, УНИПРОМЕДЬ, ИГД МЧМ, ВИОГЕМ и др., позволили решить большой круг вопросов, связанных с охраной сооружений от подработки. Определены основные схемы развития процесса сдвижения на рудных месторождениях, выполнены исследования по прогнозированию углов сдвижения, установлены условия устойчивого состояния земной поверхности, для отдельных рудных месторождений имеются нормативные документы по охране сооружений от подработки.

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

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

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

Задачи исследований:

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

2. Реализация модели среды с дилатансией, описывающей сдвижения горных пород по системам трещин в методе дискретных элементов.

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

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

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

Основные научные положения, выносимые на защиту:

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

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

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

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

Научная новизна работы:

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

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

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

Практическая значимость:

1. Экспериментально установлены угловые и линейные параметры процесса сдвижения горных пород и земной поверхности при разработке Тишинского и Риддер-Сокольного полиметаллических месторождений ОАО «КАЗЦИНК» Республики Казахстан.

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

3. Осуществлен прогноз сдвижений и деформаций массива пород при отработке рудной залежи Тишинского месторождения до 20-го горизонта.

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

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

Апробация работы. Основные результаты работы представлены на всероссийской конференции «Геомеханика в горном деле», проводимой в рамках III Уральского горнопромышленного форума (Екатеринбург, 2009 г.), международной конференции «Современные проблемы геомеханики, горного производства и недропользования» (Санкт-Петербург, 2009 г.), всероссийской конференции «Новые технологии в маркшейдерии и недропользовании» (Пермь, 2010 г.), научно-технических семинарах кафедры Маркшейдерского дела, геодезии и геоинформационных систем ПНИПУ.

Публикации. По результатам выполненных исследований опубликовано 6 работ, из них 5 в изданиях, рекомендованных ВАК.

Объем работы и ее структура. Диссертационная работа состоит из введения, четырех глав и заключения; содержит 139 страниц машинописного текста, включая 59 рисунков, 8 таблиц и список использованной литературы из 111 наименований.

Заключение диссертации по теме «Геомеханика, разрушение пород взрывом, рудничная аэрогазодинамика и горная теплофизика», Шустов, Денис Владимирович

выводы

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

Величины прироста оседаний земной поверхности в районе устья ствола «Вентиляционная» составит: при расчете методом дискретных элементов Дг| = 240мм, а горизонтальных сдвижений - = 670мм и при расчете методом конечных элементов Аг\ = 270мм, Д£, = 390мм. Однако наибольшую опасность представляют горизонтальные и вертикальные деформации массива горных пород в районе ствола. Результаты расчетов дают концентрацию как горизонтальных деформаций растяжения, так и вертикальных деформаций сжатия между 8 и 14 горизонтами, превышающую допустимые для крепи ствола значения. Так прирост горизонтальных деформаций растяжения на уровне 8-14гор. составит 0,4-0,7мм/м.Если в настоящее время горизонтальные деформации земной поверхности в районе ствола составляют 0,5-1,0мм/м, то такой прирост может привести к возникновению опасных деформаций крепи ствола. Данный вывод обязывает соответствующие службы Тишинского рудника внимательно отнестись к крепи и к состоянию ствола в целом на указанных глубинах. При появлении признаков деформации крепи и армировки ствола следует оборудовать наблюдательные станции для контроля горизонтальных и вертикальных деформаций.

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

2. Как показали сравнительные расчеты метод дискретных элементов дает более достоверные величины сдвижений и деформаций нежели метод конечных элементов при расчете подрабатываемых рудных месторождений. Средне-квадратическая ошибка расчета МДЭ в 2 - 5 раз меньше СКО расчета методом конечных элементов. Этим обосновывается третье защищаемое положение.

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

Максимальные деформации растяжения поверхности и массива горных пород в районе стволов при наличии щели составляют в пределах 1.5 мм/м. В районе щели наблюдается концентрация горизонтальных деформаций, которая на поверхности составляет свыше 4мм/м. Эти расчеты хорошо подтверждают разгрузочную функцию щели. Становится ясным, что при отказе от создания щели в 1991 г выемка запасов предохранительного целика однозначно привела бы к нарушению нормальной эксплуатации шахтных стволов и сооружений околоствольного комплекса. Т.е. щель выполняет свою функцию защиты подрабатываемых стволов шахт №3 и Вентиляционная и ее создание позволяет расконсервировать запасы руды из предохранительного целика.

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

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

ЗАКЛЮЧЕНИЕ

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

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

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

2. Для выяснения характера и степени влияния упругих и прочностных свойств массива и трещин на величины и характер сдвижений подрабатываемого массива был разработан и реализован в программе ЗЭЕС алгоритм учета деформирования по системам трещин в допредельной и запредельной стадиях и решена тестовая задача - расчет сдвижений и деформаций при создании строительного котлована в трещиноватом скальном массиве. Установлено, что для расчета сдвижений земной поверхности и массива горных пород следует использовать численную реализацию полной диаграммы деформирования массива по системам трещин с критерием разрушения Джагера, т.к. в нем используются традиционные характеристики прочности горных пород, получаемые при геомеханических испытаниях. Этими расчетами было подтверждено второе защищаемое положение.

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

4. На основе реализованного в программном комплексе ЗЭЕС алгоритма учета деформирования по системам трещин в допредельной и запредельной стадиях был выполнен прогноз сдвижений и деформаций массива горных пород и земной поверхности при отработке Тишинского месторождения до 20 горизонта. Установлено, что отработка залежи на Тишинском месторождении от 14 до 20 горизонтов не приведет к появлению деформаций ствола "Вентиляционный", которые приведут к нарушению его нормальной эксплуатации.

5. На основе реализованного в программном комплексе ЗБЕС алгоритма учета деформирования по системам трещин в допредельной и запредельной стадиях выполнены проверочные расчеты возможности управления сдвижениями массива созданием искусственной щели на рудных месторождениях. Выявлено, что данная мера эффективна для охраны ответственных сооружений, расположенных вблизи подрабатываемых бортов карьеров, а метод дискретных элементов пригоден для выполнения подобного вида расчетов.

Список литературы диссертационного исследования кандидат технических наук Шустов, Денис Владимирович, 2012 год

1. Авершин С.Г. Горные работы под сооружениями и водоемами.- М.: Уг-летехиздат, 1954. 323с.

2. Авершин С.Г. Некоторые свойства процессов сдвижения горных пород и вопросы расчета сдвижений// Труды ВНИМИ.- Л, 1961.-№43.- с.3-21.

3. Авершин С.Г. Расчет деформаций массива горных пород под влиянием подземных разработок.- Л.: ВНИМИ, i960.- 87с.

4. Акимов А.Г. Влияние трещиноватости массива на сдвижение горных пород при подземной разработке рудных месторождений // Труды ВНИМИ.-Л.,1961, с.32-45.

5. Акимов А.Г. Зависимость величин углов сдвижения от степени подра-ботанности земной поверхности // Труды ВНИМИ.-Л.,1964.-№>2. -с.99-104.

6. Акимов А.Г. Определение величины угловых параметров сдвижения на рудных месторождениях // Горный журнал.-1965.-.№2.-с. 63-67.

7. Акимов А.Г., Хакимов Х.Х. Обеспечение безопасной эксплуатации шахтных стволов. М.: Недра. 1988. - 216 с.

8. Андронова В. А. Анализ напряженно-деформированного состояния сложных регулярных стержневых конструкций, контактирующих с агрессивной средой, с использованием метода дискретных конечных элементов.: Автореф. дис. канд. тех. наук. Череповец, 1998. - 28с.

9. Ю.Ашихмин С.Г. Расчеты сдвижений и деформаций подрабатываемых скальных массивов рудных месторождений методом конечных элементов: Дис.канд.т.ех.наук./Перм.гос.тех.ун-т. -.Пермь., 1995 145.С.

10. П.Барях А. А., Стажевский С. Б., Тимофеев Е. А., Хан Г. Н. О деформированном состоянии породного массива над карстовыми пустотами. // ФТПРПИ, 2008, №6, С.3-12.

11. Барях A.A. деформирование и разрушение соляных пород и массивов: Автореф. дис.док.тех.наук Новосибирск, 1993 - 37 с.

12. И.Борщ-Компониец В.И., Макаров А.Б. Горное давление при отработке мощных пологих рудных залежей.- М.: Недра,1986.-271с.

13. Виттке В. Механика скальных пород. М.:Недра, 1988.-439с.

14. Временные правила охраны сооружений и природных объектов от вредного влияния подземных горных разработок месторождений руд цветных металлов с неизученным процессом сдвижения горных пород. -Л.: ВНИМИ. -1986. 76с.

15. Геомеханическая оценка состояния ответственных сооружений Тишинского месторождения на основе определения природного и техногенного полей напряжений. Отчет о НИР // ПермГТУ; Руководитель работы Кашников Ю.А., Пермь, 2003г. 79с.

16. Дорофеенко С. О. Моделирование сыпучих сред методом дискретных элементов: Дис. канд. физико-матем. наук/ Институт проблем химической физики РАН Черноголовка, 2008 114 с.

17. Ержанов Е.С. К теории сдвижения горных пород рудных месторождений // Труды ВНИМ,- Л., 1950. ЛВ2,- с.414-422.20.3емисев В.Н. Расчеты деформаций массива горных пород.- М.: Недра. -1973. 144с.

18. Инструкция по наблюдениям за сдвижением горных пород и земной поверхности при подземной разработке рудных месторождений.- М.:Недра, 1988. 112с.

19. Исследование и оценка динамических форм проявления горного давления на Тишинском месторождении ниже 10 горизонта: Отчет о НИР. Шифр темы 6-99-018. - Усть-Каменогорск, 1999. - 57 е. (.ВНИИцветмет).

20. Исследование проявлений горного давления на нижележащих горизонтах Тишинского месторождения и районирование горных пород по устойчивости с целью повышения безопасности отработки: Отчет о НИР/.ВНИИЦветмет, Рук. Крахин Н.С., Усть-Каменогорск, 1983.

21. Кашников Ю.А Деформирование породного массива в окрестности отрабатываемого рудного тела // Изв. вузов. Горный журнал.-1986,- №7.- с.43-51.

22. Кашников Ю.А. Анализ распределения напряжений и деформаций в подрабатываемых породах на основе расчета по программе "FEST-03".- Горный журнал.- 1993.-№12,- с.45-48.

23. Кашников Ю.А. О закономерностях сдвижения массива пород при разработке мощных крутопадающих рудных тел. В сб.: Совершенствование методов и.средств производства маркш.-.геод.работ. Л.:.ЛГИ.-1987. -.е. 16-19.

24. Кашников Ю.А. Сдвижение горных пород при подземной разработке мощных крутопадающих рудных тел сложной формы.: Дисс. канд. техн. наук. -Свердловск, 1982. 183с.

25. Кашников Ю.А., Ашихмин С.Г. Механика горных пород при разработке месторождений углеводородного сырья. М.: Недра, 2007. - 486с.

26. Кашников Ю.А., Кутовой С.Н. О параметрах призм сдвижения трещиноватого массива пород висячего бока рудных месторождений// Изв.вузов. Горный журнал.-1989.

27. Кашников Ю.А., Якушина Е.М., Ашихмин С. Г. Деформирование скального массива по системам трещин.//' Изв.вузов. Горный журнал. -1992. -№2.

28. Квочин В.А Управление сдвижением и удароопасностью массива горных пород при разработке железорудных месторождений Сибири на основе изучения их геодинамики: Дис.док.тех.наук. Новокузнецк,.2000. - 78 с.

29. Клишин В. И., Клишин С. В. Исследование процессов выпуска угля при отработке мощных пологих и крутых угольных пластов.// ФТПРПИ, 2010, №2, С.69-81

30. Клишин C.B. Применение метода дискретных элементов при анализе гравитационного движения гранулированного материала в сходящемся канале// Горный информационно-аналитический бюллетень.(научно-технический журнал), 2009, №12, С.273-277.

31. Кратч Г. Сдвижение горных пород и защита подрабатываемых сооружений.- М.: Недра, 1978. 494с.

32. Крушатин Р.Ф. Исследование процессов сдвижения горных пород при отработке слепых залежей на уральских меднорудных месторождениях.: Авто-реф. Дисс.канд.техн.наук,- Свердловск, 1969.-. 19с

33. Кузнецов М.А., Кузнецова Е.И., Громов В.В. Прогнозирование характера и параметров процесса сдвижения горних пород на рудных месторождениях / Горный журнал.- 1977.-№.-9.с.65.-69.

34. Куликов В.В. Совместная и повторная разработка рудных месторождений. М.: Недра, 1971. - 328с.

35. Курленя М.В., Миренков В.Е. Методы расчета подземных сооружений.- Новосибирск, Наука, 1986. 211с.

36. Левкин Ю.М., Иофис И.М. Взаимодействие крепи выработки с подработанным породным массивом. // Геомеханика, 2003.

37. Лобанова Т.В. Разработнка научно-методических основ геомеханического обеспечения подземной отработки железорудных месторождений Сибири в геодинамически активном регионе.: Дис.док.тех. наук. -.Новосибирск, 2009 323 е.

38. Макаров А.Б. Практическая геомеханика. Пособие для горных инженеров М.: Горная книга, 2006. - 391с.

39. Мартынов Ю.И. Управление деформированием подрабатываемого массива горных пород глубокими щелями.-М.:Стройиздат, 1983.- 116с.

40. Мартынов Ю.И., Красовский А.Н. Расчет деформаций и сдвижений пород над горизонтальными выработками // Труды ВИОГЕМ, 1972, вып. 17. -с.24-29.

41. Модели скальных массивов и алгоритмы решения задач механики горных пород / Кузнецов СВ., Одинцов В.Н., Слоним М.А. и др.// Технико-экон. эффективн. комплексного освоения месторождений.-.: МПКОН АН СССР, 1987. -C.I68-181.

42. Муллер P.A. О некоторых диффренциальных и интегральных свойствах кривых оседаний и горизонтальных сдвижений мульды сдвижения // ФТПРПИ, 1981. №2.-с. 129-132.

43. Ноек Е, Bray J.W. Rock slope engineering.- Inst. Mining and Metallurgy, London.- 1981.

44. Ножин А.Ф., Зотеев В.Г. Расчет НДС трещиноватого скального массива вокруг горных выработок // ФТРПИ, 1978.-№5. с.9-14.

45. Правила охраны сооружений и объектов от вредного влияния подземных разработок на месторождениях руд цветных металлов. -Л: ВНИМИ (проект).-1991.

46. Правила охраны сооружений и природных объектов от вредного влияния подземных горных работ в Криворожском железорудном бассейне. -Л.: ВНИМИ, 1975.- 65с.

47. Правила охраны сооружений и природных объектов от вредного влияния подземных горных работ при разработке меднорудных месторождений Урала.- М.: Недра.- 1978,- 31с.

48. Правила охраны сооружений и природных объектов от вредного влияния подземных горных разработок на угольных месторождениях. М.:Недра, 1981.-288с.

49. Правила охраны сооружений и природных объектов от вредного влияния подземных разработок на месторождениях руд черных металлов Урала и Казахстана. Свердловск: ИГД МЧМ СССР.-1990.

50. Сашурин А.Д. Особенности процесса сдвижения горных пород в условиях тектонических напряжений железорудных месторождений // Горный журнал. 1980. - М. - с.47 - 49.

51. Сашурин А.Д. Сдвижение горных пород на рудниках черной металлургии. Екатеринбург: ИГД УрО РАН. 1999. - 268с.

52. Сашурин А.Д., Беркутов В.А., Храмцов Б.А. О мерах по охране стволов шахты "Северопесчанская" от подработки // Горный журнал. 1979. - .№5. -с.34-36.

53. Сашурин А.Д.,Беркутов В.А., Бекселев АЛИ. Расчет горизонтальных деформаций земной поверхности вокруг провалов на железорудных месторождениях." В кн.: Подземная разработка мощных рудных месторождений. -Свердловск, 1979. с.36-45.

54. Сдвижение горных пород на рудных месторождениях / Кузнецов М.А, Акимов А.Г., Кузьмин В.И. и др М.: Недра, 1971.-224с.

55. Смирнов А.Ф. О возможности предрасчета углов разрывов на рудных месторождениях // Труды ВНИМИ.-Л.,1964,№53.-с.239-251.

56. Смирнов А.Ф. Предрасчет углов разрывов для условий Кривбасса. // Горный журнал.- 1980.-Ш.-С. 52-54.

57. Современные проблемы механики скальных пород в энергетическом строительстве. М.:Энергоиздат, 1986. 302с.

58. Теоретические методы расчета сдвижений и деформаций земной поверхности, вызванных подземными горными работами.// Муллер P.A., Петухов И.А. и др. М.: ЦНИЗИуголь, 1977.- с.4-17.

59. Указания по охране сооружений и природных объектов от вредного влияния подземных горных разработок на полиметаллических месторождениях Лениногорского горно-обогатительного комплекса ОАО «.КАЗЦИНК.». Усть-Каменогорск. 2001г. - 54с.

60. Фадеев А.Б. Метод конечных элементов в геомеханике. М.:Недра, 1987.-221с.

61. Хан Г.Н. Исследование процесса обрушения смерзшейся породы методом дискретных элементов // Материалы 3-й Межд. конф. "Наукоемкие технологии добычи и переработки полезных ископаемых",.-.Новосибирск: Изд-во ИГД СО РАН, 2003. с.98-103.

62. Хан Г.Н. О несимметричном режиме разрушения массива горных пород в окрестности полости. // Физическая мезомеханика, 2008, №1, С.109-114.

63. Хрисчев Г., Мотов К., Георгиев К. Опазване на сьорьженията и обек-тите от вредното влияние на подземните минни работа.- София, Техника.-1978.— 268с.

64. Черный Г.И. Определение величин оседаний и деформаций земной поверхности при сдвижении пород в форме реологического течения // Изв. вузов. Горный журнал.- 1966.- №7.- с.3-9.

65. Черный Г.И. Устойчивость подрабатываемых бортов карьеров.- М.: Недра, 1980. 216с.

66. Шадрин А.Г. Теория и расчет сдвижений горных пород и земной поверхности.- Красноярск: изд-во Красноярского ун-та.-1990.-200с.

67. Юнаков Ю.Л. Исследование сдвижения горных пород и устойчивости горных выработок при комбинированном способе отработки месторождений: Автореф. Дисс канд.техн.наук.- Караганда.- 1980.

68. Andersson, J. and Dverstop, В., Conditional simulations of fluid flow in three-dimensional networks of discrete fractures. Water Resources Research, 1987;23(10): 1876-1886.

69. Ante Munjiza, The Combined Finite-Discrete Element Method Wiley, 2004.

70. Barbosa, R. and Ghaboussi, J., Discrete finite element method Proc. 1st US Conf. on Discrete Element Methods, Golden, Colorado, 1989.

71. Burman, В. C, A numerical approach to the mechanics of discontinua. Ph.D. Thesis, James Cook University of North Queensland, Townsville, Australia, 1971.

72. Byrne, R. J., Physical and numerical model in rock and soil-slope stability. Ph.D. Thesis, James Cook University of North Queensland, Townsville, Australia, 1974.

73. Chappel, B. A., The mechanics of blocky material. Ph.D. Thesis, Australia National University, Canberra, 1972.

74. Charlaix, E., Guy on, E. and Roux, S., Permeability of a random array of fractures of wisely varying apertures. Transport in Porous Media, 1987;2(1):31-43.

75. Cundall, P. A. and Hart, R. D., Numerical modelling of discontinua. Keynote lecture. Proc. 1st US Conf. on Discrete Element Methods, Mustoe, G. G. W., Henriksen, M. and Huttelmaier,.Golden, Colorado, USA, 1989.

76. Cundall, P. A. and Strack, O. D. L., A discrete numerical model for granular assemblies. Geotechnique, 1979a;29(l):47-65.

77. Cundall, P. A., A computer model for simulating progressive, large scale movements in blocky rock systems. Proc. Int. Symp. Rock Fracture, ISRM, Nancy, Vol. 1, PaperII-8, 1971.

78. Cundall, P. A., Adaptive density-scaling for time-explicit calculations Proc. 4th Int. Conf. on Numerical Methods in Geomechanics, pp. 23-26, 1982.

79. Dershowitz, W. S. and Einstein, H. H., Three-dimensional flow modelling in jointed rock masses. In: Herget, G. and Vongpaisal, S. (eds), Proc. of 6th Cong. ISRM, Montreal, Canada, Vol.1,.-87.92,.1987.

80. Elsworth, D., A hybrid boundary-element-finite element analysis procedure for fluid flow simulation in fractured rock masses. International Journal for Numerical and Analytical Methods in Geomechanics, 1986b; 10(6):569-584.

81. Endo, H. K., Long, J. C. S., Wilson, C. K. and Witherspoon, P. A., A model for investigating mechanical transport in fractured media. Water Resources Research, 1984;20(10): 1390-140.0.

82. Goodman, R. E., Methods of geological engineering in discontinuous rocks. West Publishing Company, San Francisco, CA, 1976.

83. Itasca (2007). Itasca Software Products 3DEC. Itasca Consulting Group Inc.: Minneapolis.

84. Jaeger G.C. Friction of Rocks and stability of rock slopes. Geotechnique 21/2(1971).95Jantson F. Die Syarmungsverteilimg in kreisformigen Schachtzonen unter Beruckaichtiguag elastichen Gebirgsver-haltens. Diss. Aachen, 1966.- 191s.

85. Jing L., Stephansson O. Fundamentals of Discrete Element Methods for Rock Engineering. Elsevier Science, 2002. - 450 c.

86. Jing, L., A review of techniques, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering. International Journal of Rock Mechanics and Mining.Sci.enc.es,.200.3;40(3):283-353.

87. Leichnitz W. Mechanische Eigenschaften von Felstrennflachen im direkten Scherversuch. Veroeff.des Instfuer Bodenmechanik und Felsmechanik der TH Karlsruhe, Heft 89.-1981.

88. Long, J. C. S., Gilmour, P. and Witherspoon, P. A., A model for steady fluid flow in random three dimensional networks of disc-shaped fractures. Water Resources Research, 1985;21(8):.l 10.5-1.115.

89. Pande, G., Beer, G. and Williams, J.R., Numerical Modeling in Rock Mechanics, John Wiley and Sons, 1990.101 .Robinson, P. C., Flow modelling in three dimensional fracture networks. UK AEA Harwell, AERER 11965, 1986.

90. Salamon M.D.G. Elastic modul of a stratified rock mass: Int. J. Rock Mech. Min. Sci., (5), 1968.- 519-527.

91. Salamon M.D.G. Rock mechanics of underground excavations: Proc. 3rd Congr. Int. Soc. Rock Mech. (1 ), 1974, Part B, Denver.-951-1099.

92. Shi, G., Discontinuous deformation analysis a new numerical model for statics and dynamics of block systems. Ph.D. Thesis, University of California, Berkeley, CA, 1988.

93. Smith, L. and Schwartz, F. W., An analysis of the influence of fracture geometry on mass transport in fractured media. Water Resources Research, 1984;20(9): 1241-1252.

94. Tsang, Y. W. and Tsang, C.-F., Channel model of flow through fractured media. Water Resources Research, 1987;22(3):467^179.

95. W.Wittke. Rock Mechanics. Teory and Applications with Case Histories. Sprlnger-Verlag.