автореферат диссертации по разработке полезных ископаемых, 05.15.02, диссертация на тему:Научные основы расчета конструктивных параметров систем подземной разработки руд с учетом структуры массива и порядка ведения горных работ

ВВЕДЕНИЕ

1. ИЗУЧЕННОСТЬ ВОПРОСА И ЗАДАЧИ ИССЛЕДОВАНИЙ.

1.1. Особенности ведения горных работ на больших глубинах.

1.2 Методы управления устойчивостью элементов систем разработки

1.3. Модели массива горных пород и методы моделирования его напряженно-деформированного состояния.

1.4. Цель и задачи исследований.

2. МЕТОД РАСЧЕТА НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ НЕОДНОРОДНОГО МАССИВА ГОРНЫХ ПОРОД

2.1. Основные положения метода конечных элементов применительно к расчету напряженно-деформированного состояния скальных массивов.

2.2. Основные численные процедуры метода конечных элементов

2.3. Комплекс программ для расчета НДС скальных массивов с учетом неоднородности деформационных свойств и первоначального напряженного состояния.

Выводы.

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

3.1. Понятие структуры массива, ее основные параметры и методы изучения.

3.2. Физико-механические характеристики трещин и способы их учета при расчетах НДС и устойчивости подземных конструкций

3.3. Контакт-элемент для моделирования трещин при решении двумерных задач.

3.4. Контакт-элемент для моделирования трещин при решении объемных задач.

Выводы.

4. ИЗУЧЕНИЕ И УЧЕТ НАПРЯЖЕННОГО СОСТОЯНИЯ МАССИВА В ГЕОМЕХАНИЧЕСКИХ РАСЧЕТАХ

4.1. Методы измерения напряжений

4.2. Влияние неоднородности деформационных свойств массива на исходное поле напряжений.

4.3. Влияние структуры массива на поле первоначальных напряжений

4.4. Методический подход к измерению напряженного состояния скального массива.

Выводы.

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

5.1. Учет порядка изменения конфигурации выработанного пространства

5.2. Учет последовательности очистных работ с закладкой выработанного пространства при расчетах НДС массива скальных пород

5.3. Учет условий контактов закладки с вмещающими породами . 194 Выводы.

6. РЕАЛИЗАЦИЯ РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЙ.

6.1. Геомеханическое обеспечение горных работ Высокогорского

ГОКа.

6.1.1. Оптимизация параметров компенсационных камер на шахте "Магнетитовая".

6.1.2. Моделирование фильтрационных свойств массива при заполнении шламами Главного карьера.

6.2. Геомеханическое обеспечение подземных горных работ на шахтах Гороблагодатского рудоуправления.

6.2.1. Оценка условий отработки запасов нового горизонта -240 м шахты "Южная".

6.2.2. Изучение причин разрушения потолочинных целиков на шахте "Валуевская".

6.3. Геомеханическое обеспечение подземных горных работ на шахтах Донского ГОКа.

6.3.1. Мероприятия по повышению устойчивости выработок на шахте "Молодежная".

6.3.2. Определение напряженного состояния массива на сопряжении ствола с квершлагом на шахте" Центральная".

6.4. Определение несущей способности предохранительного целика Гологорского карьера на шахте "Магнезитовая"

6.5. Выбор параметров камерной системы разработки с обрушением налегающих пород на шахте "Сидеритовая" Бакальского РУ

6.6. Оценка запаса устойчивости элементов камерно-столбовой системы разработки на шахте "Центральная" Джездинского РУ

6.7. Выбор параметров этажно-камерной системы разработки с обрушением налегающих пород на шахте "Западный Каражал"

Атасуйского ГОКа.

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

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

В основу диссертации положены результаты более 20 плановых научно-исследовательских хоздоговорных и госбюджетных работ, выполнявшихся с участием автора в лаборатории горного давления ИГД УрО РАН (ИГД МЧМ СССР) и на кафедре разработки рудных месторождений УГГГА. В последние годы работа выполнялась при поддержке фонда фундаментальных исследований Минвуза РФ (конкурс грантов 1997 г., грант "Разработка методики расчета напряженно - деформированного состояния массива горных пород с учетом стадийности развития горных работ") и Российского фонда фундаментальных исследований (грант 99-05-64263 "Влияние иерархии и ориентации трещин на динамические явления в массиве скальных пород при разработке месторождений полезных ископаемых").

Объект исследований - массив скальных пород в зоне влияния выработанного пространства.

Предмет исследований - методы расчета напряженно - деформированного состояния (НДС) и устойчивости элементов систем подземной разработки.

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

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

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

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

2. Разработать методы расчета НДС массива с учетом развития горных работ.

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

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

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

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

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

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

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

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

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

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

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

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

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

3. Разработаны контакт - элементы для моделирования протяженных трещин при решении плоских и объемных задач деформируемого твердого тела.

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

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

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

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

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

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

Апробация работы. Основные результаты работы докладывались на Всесоюзном семинаре "Численные методы в горном деле" (Апатиты, 1986), территориальной научно-технической конференции "Методы оценки напряженного состояния массивов горных пород при разработке полезных ископаемых Урала" (Свердловск, 1987), международных конференциях по механике горных пород (1988 г. - Австрия и 1990 г. - Португалия), Международном симпозиуме "Проблемы безопасности при эксплуатации месторождений полезных ископаемых в зонах градопромыш-ленных агломераций" (1995 г - г. Москва - г. Пермь), Всероссийских конференциях "Управление напряженно-деформированным состоянием массива горных пород при открытой и подземной разработке месторождений полезных ископаемых" (1994 г., 1996 г. - г. Екатеринбург), Международной конференции "Экологические проблемы горного производства, переработка и размещение отходов" (1995 г. - г. Москва) и постоянной действующем региональном семинаре по механике горных пород при ИГДУрО РАН (1999 г.).

Публикации. Основные положения диссертации опубликованы в 24 печатных трудах.

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

Выводы

1. Разработанный метод расчета НДС с учетом последовательности выемки запасов позволяет определить условия возникновения, а также размеры и форму областей разрушения массива (зон неупругой деформации) и выбрать оптимальный порядок отработки запасов.

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

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

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

6. РЕАЛИЗАЦИЯ РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЙ

Разработанные метод и комплекс программ для расчета НДС скальных массивов широко применялись автором при геомеханическом обеспечении подземных горных работ, в основном на рудниках черных металлов Урала и Казахстана. Разработанный комплекс программ также используется при проведении практических занятий по курсу Теомеханика" и спецкурсу "Расчет НДС линейно- и нелинейно-упругих сред" для студентов специализации "Подземная разработка рудных месторождений" Уральской горно-геологической академии. Кроме того, в 1997 г комплекс программ был передан лаборатории горного давления ИГД УрО РАН.

6.1. Геомеханическое обеспечение горных работ Высоко-гоского ГОКа

6.1.1. Оптимизация параметров компенсационных камер на шахте "Магнетитовая"

Запасы блока № 21 шахты "Магнетитовая" Высокогорского ГОКа отрабатываются системой этажного принудительного блокового обрушения руды на компенсационную камеру. Рудное тело представляет собой штокообразную наклонную залежь. Размеры залежи в плане на гор. -450 м составляют примерно 200м х 200 м. Прочность образцов руды и пород колеблется в пределах 100 - 170 МПа. До горизонта -370 м (отметка дневной поверхности +220м) сложностей с обеспечением устойчивости элементов систем разработки не возникало. При отрабтке запасов первой панели горизонта -450 м были зарегистрированы большие объемы самообрушения стенок и кровли компенсационной камеры, приведшие к потере вееров скважин и части доставочных выработок.

Проведенный комплекс натурных экспериментов по определению напряженного состояния [148, 149] показал, что напряженное состояние массива по сравнению с гор. -370м изменилось незначительно, величина вертикальных напряжений составила 40 МПа, а горизонтальных - 2025 МПа. При этом было отмечено увеличение интенсивности трещино-ватости рудного массива по сравнению с вышележащим горизонтом.

Величины напряжений на контуре камеры, найденные с учетом объемности составляют 20-25 МПа сжатия по целику между камерой и зоной обрушения и 1-2 МПа растяжения в кровле камеры. Изменение размеров камер в разумных пределах существенных изменений в распределение напряжений не вносит.

Первоочередными рекомендациями было уменьшение размеров компенсационных камер в плане с оформлением у них шатровой кровли. Эта мера, хоть и незначительно, увеличивает прочность обнажения массива, а также сокращает сроки формирования камер. В результате ВГОК уменьшил размеры компенсационных камер с 50м х 15м до (20-25)м х (12-15)м. В дальнейшем, на основании результатов решения объемных линейно-упругих задач теории упругости, было предложено оформлять камеры в виде песочных часов в вертикальном сечении (рис. 6.1). Использование наклонных плоскостей приводит к разгрузке центральной части камер от вертикальных напряжений, что позволяет увеличить размеры таких камер по сравнению с камерами прямоугольной формы даже в неустойчивых массивах. Отработка такой камеры в опытной панели показала применимость такого технологического решения.

Компенсационная камера для условий шахты "Магнетитовая"

Проведенные расчеты также показали, что около выработок доставки формируется зона неупругих деформаций мощностью до 1,5 м. Поэтому шахте было рекомендовано увеличить несущую способность крепей откаточных и скреперных выработок. От крепления набрызг-бетоном шахта перешла к использованию арочной крепи из СВП-27 с плотностью 0,6 - 0,7 арки/м и забутовкой закрепного пространства. Скреперные выработки в зоне неустойчивых пород, примыкающей к висячему боку залежи, стали крепить бетоном.

Описанный комплекс мероприятий позволил вести отработку запасов гор. -470 м без больших осложнений.

6.1.2. Моделирование фильтрационных свойств массива при заполнении шламами Главного карьера

В проекте реконструкции Черемшанского хвостохранилища было предусмотрено складирование хвостов МОФ в выработанное пространство Главного карьера. Борта и дно карьера сложены известняками, имеющими сложное блоковое строение. В пределах контура карьера за-картировано 5 крупных крутопадающих тектонических зон дробления мощностью 20-30 м. Кроме того, известняки расчленены 3 основными системами трещин, имеющими раскрытие от 3-5 до 200 мм. При инженерно-геологической съемке в известняках было выявлено несколько карстовых каналов сечением до 0,4 м. И, наконец, массив подработан шахтой "Магнетитовая", причем восточный борт захвачен зоной обрушения, образовавшейся при отработке Западно-Ревдинской залежи. Провальные воронки выходят в борт карьера выше гор. +110 м, и в настоящее время полностью перекрыты отвалом из полускальных и глинистых пород.

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

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

Граничные условия для моделирования принимались по результатам измерения напряженного состояния, полученными ИГД МЧМ СССР [151]. Скальный массив моделировался линейно-упругой средой с модулем деформации 5000 МПа и коэффициентом Пуассона 0,3. Породы, заполняющие зону обрушения, и отходы обогатительной фабрики аппроксимировались идеально сыпучей средой. Изучение изменения раскрытия трещин велось на основе решения линейно-упругих плоских задач [167, 172 - 175]. Результаты расчетов показаны на рис. 6.2.

Изменение раскрытия трещин в борту Главного карьера а)

200 400 600 а - раскрытие трещин в незаполненном карьере; б - раскрытие трещин в заполненном карьере; положение противофильтрационных экранов;

-2" - раскрытие трещин (мм)

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

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

6.2. Геомеханическое обеспечение подземных горных работ на шахтах Гороблагодатского рудоуправления

6.2.1. Оценка условий отработки запасов нового горизонта -240 м шахты "Южная"

В связи с доработкой запасов этажа -80 - 160 м перед шахтой "Южная" встал вопрос об организации очистных работ на нижележащем горизонте -240 м. Целью работы являлось прогнозирование геомеханических характеристик массива и разработка технологических рекомендаций по совершенствованию отработки блока № 1 нового горизонта.

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

Выполненный комплекс исследований включал в себя измерения НДС массива и оценку его деформационных свойств экспериментально-аналитическим методом и методом сейсмометрии, а также расчет НДС массива для различных вариантов отработки, как блока, так и запасов этажа в целом [41, 42]. Результаты натурных экспериментов показали, что величина первого инварианта тензора напряжений с глубиной значительно увеличивается (рост отношения скоростей продольных и поперечных упругих волн составил 50 - 65%). Кроме того, экспериментально было выявлено формирование зон неупругих деформаций мощностью до 4 м вокруг откаточных выработок. Отметим, что при этом очистные работы на горизонте не велись.

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

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

6.2.2. Изучение причин разрушения потолочинных целиков на шахте "Валуевская"

Целью данной работы являлось установление причин разрушения потолочинного целика на шахте "Валуевская" и разработка рекомендаций по обеспечению его устойчивости [176]. Шахтой отрабатывалось наклонное рудное тело мощностью около 15 м. Камерные запасы извлекались нормально. Однако при этом целики сильно деформировались, скважины в них терялись, а иногда происходило их частичное разрушение (без отбойки руды выпускалось до 20-80% запасов целика при отработке камеры).

Вмещающие породы, представленные в висячем боку порфирита-ми и в лежачем боку сиенитами, имеют близкие прочностные и деформационные характеристики (модуль деформации вмещающих пород около 40000 МПа, коэффициент Пуассона равен 0,24 - 0,27). Магнетито-вые руды имеют пониженные физико-механические характеристики (модуль деформации около 9000 МПа, коэффициент Пуассона 0,3, прочность на одноосное сжатие 40-50 МПа). Рудное тело контролируется тектоническими разломами, представляющими собой зоны дробленых ми-лонитизированных пород мощностью около 0,4 м (рис. 6.3). Модуль деформации материала этих зон около 1000 МПа, коэффициент Пуассона 0,25, сцепление 0,5 МПа, угол внутреннего трения материала около 20°. Первоначальные вертикальные напряжения определяются весом налегающих пород, горизонтальные напряжения, действующие вкрест простирания рудного тела, превышают вертикальные на 10 МПа.

Положение (а) и деформация (б) потолочинного целика на ш. "Валуевская"

Рис. 6.3

Результаты моделирования НДС массива показывают, что отработка камеры вызывает смещение целика по тектоническим разломам в направлении восстания рудного тела на 10 - 16 см. При этом смещение пород висячего бока на 8-14 см больше, чем смещение целика относительно пород лежачего бока. Руда в целике находится практически в одноосном напряженном состоянии, величина сжимающих напряжений колеблется около 40-45 МПа, т.е. близка к пределу прочности рудного массива на одноосное сжатие.

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

6.3. Геомеханическое обеспечение подземных горных работ на шахтах Донского ГОКа

6.3.1. Мероприятия по повышению устойчивости выработок на шахте "Молодежная"

Шахта отрабатывает запасы хромитовых руд Кимперсайской группы месторождений. Массив вмещающих пород и руд разбит интенсивной крутопадающей трещиноватостью, имеющей хаотическое распределение по азимутам падения. Характерной чертой трещин является наличие низкопрочного заполнителя - серпофита и талькобрейнерита, имеющего прочность 0,2 - 0,5 МПа на одноосное сжатие и угол трения, равный 18 -20°. Мощность заполнителя от долей миллиметра до 2 - 3 см и более, в среднем - 1-2 мм. На основании анализа разрушений пород на обнажениях горизонтальных подсечек и вертикальных отрезных щелей, моделирования деформации контуров выработок вблизи выработанного пространства в качестве модели массива была выбрана среда с хрупким разрушением. Значения начальных модуля деформации и сцепления массива оценены как 2000 МПа и 1 МПа, величина остаточного сцепления близка к 0,1 МПа. Значения начального и остаточного углов внутреннего трения совпадают и равны 25°. Натурные наблюдения по скважинам с помощью РВП-27 и решение ряда нелинейно-упругих задач позволило установить факт формирования в окрестности выработанного пространства зоны неупругих деформаций большого (от 30 до 80 м и более) размеров (рис. 6.4). На основании полученных результатов шахте было рекомендовано увеличить расстояние между границей выработанного пространства и штреком лежачего бока, а также увеличить несущую способность крепи [177].

6.3.2. Определение напряженного состояния массива на сопряжении ствола с квершлагом на шахте "Центральная"

При проходке стволов шахты "Центральная" несколько раз возникали аварийные ситуации, выражавшиеся в обрушении пород из-за за-крепного пространства или в больших (до 10 - 15 см) деформациях уже пройденных и закрепленных участков стволов. Для уточнения параметров крепи и технологической схемы проходки стволов было необходимо уточнить напряженное состояние, формирующееся на сопряжении круглого ствола с горизонтальными выработками.

Распределение главных нормальных напряжений и положение ЗНД в окрестности выработанного пространства на ш."Молодежная"

Рис. 6.4

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

6.4. Определение несущей способности предохранительного целика Гологорского карьера на шахте "Магнезитовая"

Проект подземной разработки Саткинского месторождения магнезита предусматривает применение подэтажно-камерной системы с последующей закладкой пустот твердеющим материалом. В связи с отставанием строительства закладочного комплекса было принято решение вести отработку опытного блока открытыми камерами без их закладки с оставлением междукамерных целиков (МКЦ). Выемка МКЦ будет производиться после окончания строительства закладочного комплекса и заполнения первичных камер закладочными смесями. Размеры камер и МКЦ были приняты одинаковыми: высота камер 10 м, их ширина -8 м. Выемка магнезита должна осуществляться с применением самоходной техники.

Опытный блок располагается между ныне действующим Карагай-ским и уже отработанным Гологорским карьерами (рис. 6.5). Гологорский карьр заполнен вскрышными породами. Для предотвращения проникновения отвальных пород из Гологорского карьера в очистные камеры опытного блока был предусмотрен предохранительный целик из магнезита толщиной 30 м.

Рис. 6.5

Параметры целика при составлении проекта не были обоснованы. В ходе уточнения проектных решений возникла необходимость уточнения несущей способности целика. Моделирование напряжений в целике проводилось для двух стадий отработки запасов: первая стадия соответствовала моменту полной отработки всех камер первой очереди, вторая - моменту полной отработки запасов и заполнения выработанного пространства закладочным материалом [179]. Порядок выемки запасов вторичных камер не рассматривался. Результаты моделирования не только позволили оценить запас устойчивости целика, но и позволили обосновать уменьшение его толщины до 20 м, т.е. на 30%.

6.5. Выбор параметров камерной системы разработки с обрушением налегающих пород на шахте "Сидеритовая" Бакальского РУ

Блоком № 20 шахты "Сидеритовая" отрабатывалось слепое рудное тело, представляющее собой пластообразную залежь с углом падения 25-30° при максимальной мощности 70 м. Добываемые руды представлены сидеритом, лежачий бок залежи сложен доломитами. В висячем боку непосредственно над рудным телом залегает толща кварцево-глинистых сланцев, имеющих мощность до 70 м. Выше залегают кварциты, представляющие собой крепкий устойчивый массив. Необходимость применения камерной системы разработки обуславливается высокими требованиями, предъявляемыми к качеству руды: разубоживание не должно превышать 10%.

Первая очередь 20 блока отрабатывалась О-образной камерой с оставлением в центре столбчатого целика для поддержания кровли, сложенной неустойчивыми сланцами. Тем не менее, при отработке камеры происходило обрушение ее кровли, а засорение руды на конечной стадии отработки превысило 40%.

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

Средняя прочность на одноосное сжатие рудного массива составила 20 МПа, сланцев - 5 МПа. Величина горизонтальных напряжений, действующих по простиранию блока 12 МПа, вкрест простирания - 18 МПа. Вертикальные напряжения соответствовали весу налегающей толщи и на глубине 220 м составили 6 МПа. На основе решения объемных задач теории упругости шахте было рекомендовано отрабатывать 2-ю очередь блока двумя камерами, разделенными целиком. Ширина камер проектировалась равной 45 м, длина -110 м. Ширина целика -15 м, в центральном поперечном сечении - 30 м (рис. 6.6). В камере лежачего бока, высота которой более 55 м, предлагалось шатровое оформление кровли с оставлением рудного целика под неустойчивыми сланцами.

Отработка запасов 2-й очереди блока (1,5 млн. тонн) продолжалась более 3 лет. Самообрушений при этом не наблюдалось. Экономический эффект в ценах 1988 г без учета обогатительного передела по общешахтным расходам составил более 200 тыс. руб. в год.

Схема отработки блока № 20

01=27 МПа

1-я очередь 2-я очередь б)

1 очередь (проект БРУ)

3-я очередь 4-я очередь

2 очередь (проект ИГД)

Рис. 6.6 а - план отработки блока N2 20; б - план отработки и форма камер при отработке 2-й очереди блока сидерит кварцевого инистые сланцы хххххххххххх: хх хх хх хх хх хх:

XX XX XX XX XX XX: ХХХ.ХХІІХ.ХХХХ.Х: доломит проект БРУ — — — проект ИГД очередность отработки

6.6. Оценка запаса устойчивости элементов камерно-столбовой системы разработки на шахте "Центральная" Джездинского РУ

Шахта отрабатывает пологопадающее рудное тело марганцевых руд со средней мощностью около 3 м. Кровля сложена песчаниками и конгломератами, подстилающие породы представлены гранитами. На некоторых участках на расстоянии 1,0 - 1,5 м от кровли пласта залегают глинистые отложения мощностью до 1 м, отделенные от марганцевых руд прослоями песчаника. Глубина залегания месторождения 100 - 200 м.

Для определения параметров крепи требовалось оценить запас устойчивости целиков и кровли камер в зависимости от расстояния до глинистых отложений и их мощности. Поверхность скольжения, которая может образоваться при вдавливании целиков в кровлю пласта, определяется решениями, полученными на основе метода предельного равновесия [61]. Коэффициент запаса устойчивости целиков по условиям их оснований определялся по этой поверхности методом горизонтальных сил с использованием расчетных величин напряжений. Для уточнения полученных результатов также был решен ряд нелинейно-упругих задач. -Результаты расчетов вошли составной частью в " Инструкцию по выбору и расчету параметров крепи." [180].

6.7. Выбор параметров этажно-камерной системы разработки с обрушением налегающих пород на шахте "Западный Каражал" Атасуйского ГОКа

Шахтой "Западный Каражал" отрабатывается железорудное месторождение, сложенное гематитовыми рудами. С увеличением глубины ведения горных работ участились случаи самопроизвольных обрушений кровли и стенок камер. В 1992 г. лаборатория горного дела ИГД УрО РАН разработала по просьбе Атасуйского ГОКа рекомендации по отработке блока № 5 на гор. +108 м. Первоначальный вариант проекта был составлен проектно-конструкторским отделом ГОКа совместно с ИГД Казахстана. По проекту запасы блока отрабатывались 11-ю различно ориентированными камерами, параметры которых были определены по аналогии с камерами вышележащего горизонта (рис. 6.7).

Комплекс работ по геомеханическому обоснованию отработки запасов, выполненный ИГД УрО РАН, включал измерения напряженного состояния массива методами щелевой разгрузки и электрометрии (метод ПЭЗ), моделирование НДС элементов систем разработки, предложенных в первоначальном проекте и разработку рекомендаций по порядку выемки запасов и параметрам камер и целиков. В результате шахте рекомендовано отказаться от междукамерного целика (МКЦ), разделяющего камеры 1 и 2 блока 5А, блок 5Б отрабатывать без оставления западного технологического целика, а также предусмотреть работы по принудительному обрушению кровли камеры 2 этого блока (рис. 6.7). Кроме того, было рекомендовано изменить положение МКЦ 1-2 в блоке 5В.

К настоящему времени блок N25 отработан практически полностью без существенных отклонений от проектных параметров.

План расположения камер блока №5 на гор. +108 м

Разрез по линии XV - XV

1Ш

Разрез по линии VII - VII

К,::

Разрез по линии XI - XI

Разрез по линии XIII - XIII

Рис. 6.7

- РУДЭ

- известняки кремнистые породы проект ГОКа —рекомендации ИГД

ЗАКЛЮЧЕНИЕ

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

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

1. Разработаны методология и комплекс программ для расчета напряженно-деформированного состояния массива с учетом его структуры, первоначального напряженного состояния, неоднородности физико-механических свойств пород и последовательности развития горных работ. Комплекс программ позволяет решать задачи в плоской и в объемной постановке. При этом практически снято ограничение по предельному количеству узлов.

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

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

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

5. Доказано, что необходим дифференцированный подход к заданию граничных условий при проведении расчетов НДС массива в окре

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

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

7. Разработано математическое описание параметров модели упрочняющейся среды для моделирования контакта закладки с вмещающими породами, а также блоков закладки из нетвердеющих смесей.

8. Результаты исследований применены проектными отделами рудников и институтом Уралгипроруда при решении практических задач на 9 шахтах, добывающих руды черных металлов. Экономический эффект по шахте "Сидеритовая" составил за 3 года 600000 руб. в ценах 1988 г.

Основные положения диссертации опубликованы 24 работах:

1. Бронников Д. М., Замесов Н.Ф., Богданов Г.И. Разработка руд на больших глубинах. М.: Недра, 1982. - 292 с.

2. Каплунов Д.Р., Левин Д.И., Александров A.A., Жегалин Я.А. Основные направления совершенствования техники и технологии разработки пологопадающих мощных и средней мощности месторождений крепки руд // там же. С. 89 - 100.

3. Курленя М.В., Штеле В.И., Шалауров В.А. Развитие технологии подземных горных работ. Новосибирск: Наука, 1985. - 184 с.

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

5. Черненко А.Р. Основные тенденции совершенствования технологии подземной добычи руд Кривбасса на больших глубинах // там же. -С. 86-89.

6. Славиковский О.В., Осинцев В.А. Подземная разработка рудных месторождений на Урале и основные направления совершенствования технологии добычи // Известия вузов. Горный журнал. 1997. - № 11 -12. -С. 90-117.

7. Зубков A.B. Проблемы, стоящие перед горным недроведением на современном этапе освоения недр // Проблемы геотехнологии и не-дроведения (Мельниковские чтения): Докл. международн. конф. Екатеринбург. - 1998. - Т. 1. - С. 3 - 5.

8. Яковлев В.Л., Волков Ю.В., Славиковский О.В. Основные направления физико-технических проблем освоения георесурсов подземной геотехнологией // там же. Т. 3. - С. 9-26.

9. Яковлев М.А., Жуков H.A. Технологические решения при подземной разработке рудных месторождений // Геомеханическое обоснование технологических решений при разработке руд подземным способом: Сб. науч. тр. / ИГД СО АН СССР. Новосибирск, 1984. - С. 3 -8.

10. Петухов И.М., Егоров П.В., Винокур Б.Ш. Предотвращение горных ударов на рудниках. М.: Недра, 1984. - 230 с.

11. Марков Г.А. Определение тектонических напряжений в верхней части земной коры // Методические и технические средства определения напряжений в горном массиве: Сб. науч. тр. Новосибирск: ИГД СО АН СССР, 1983.-С. 10-15.

12. Бич А.Я., Емельянов Б.И., Муратов A.A. Управление состояние массива горных пород. Владивосток: Дальневосточн. ун-тет, 1988. -284 с.

13. Шевяков Л.Д. О расчете прочных размеров и деформаций целиков. Изв. АН СССР. Отделение технических наук, 1941. - N 7, 8. - С. 313; N 9. - С. 43-58.

14. Слесарев В.Д. Механика горных пород. М.: Углетехиздат,1948. -303 с.

15. Методы определения размеров опорных целиков и потолочин. -М.: Академиздат, 1962. 200 с.

16. Фисенко Г.Л. Предельные состояния горных пород вокруг выработок,- М.: Недра, 1976. 272 с.

17. Методические указания по установлению размеров камер и целиков при камерной системе разработки руд цветных металлов / Нестеренко Г.Т., Скобцев Б.С., Палий В.Д. и др. Л.: ВНИМИ, 1972. - 21 с.

18. Ерофеев Н.П. Прогнозирование устойчивости горных выработок. Алма-Ата, 1977. - 80 с.

19. Разработка технологии и определение параметров разработки нижележащих горизонтов Соколовского подземного рудника: Отчет о НИР /ИГД МЧМ СССР. Рук. Влох Н. П., Шуплецов Ю. П. Свердловск, 1986. - 123 с.

20. Исследование и внедрение порядка разработки и параметров очистной выемки в блоках N 8 и N 20 ш. "Сидеритовая" БРУ : Отчет о НИР /ИГД МЧМ СССР. Рук. Шуплецов Ю. П. Свердловск, 1988. - 44 с.

21. Обобщение результатов исследований опытно-промышленных работ с обрушением на Соколовском подземном руднике: Отчет о НИР/ ИГД МЧМ СССР. Рук. Влох Н. П. Свердловск, 1989. - 132 с.

22. Руководство по определению безопасных объемов пустот, выбору и расчету изолированных сооружений при подземной разработке рудных месторождений Киргизии. Фрунзе: Илим, 1976. - 51 с.

23. Зубков A.B. Разработка методов управления горным давлением на основе решения трехмерных геомеханических задач: Автореф. дисс. докт. техн. наук / ИГД СО АН СССР. Новосибирск, 1990. -37 с.

24. Ефимов Ю.Н., Сапожников Л.В. Реализация метода конечных элементов для решения плоской и пространственной задач теории упругости// Известия ВНИИГ. Л.:Энергоатомиздат,1985. - Т. 186. - С. 3-6.

25. Зотеев О.В. Разработка методики математического моделирования напряженно-деформированного состояния массива горных пород: Дисс. канд. техн. наук / ИГД ММ СССР. Свердловск, 1991. - 126 с.

26. Опыт прогноза динамических явлений при производстве массовых взрывов / Курленя М.В., Еременко A.A., Никитин В.В. и др. Горный журнал. - 1994. - № 6. - С. 36 - 39.

27. Редькин В.А., Анохин А.Г., Селяев И.С. Выбор параметров ка-муфлетных взрывов для разгрузки удароопасных участков. Горный журнал. - 1990. - № 5. - С. 54-56.

28. Изыскание методов управления горным давлением при разработке мощных железорудных месторождение, склонных к горным ударам: Отчет о НИР / ИГД МЧМ СССР. Рук Зубков A.B. Свердловск,1986. -119с.

29. Баклашов И.В. Деформирование и разрушение породных массивов. М.: Недра, 1988. - 271 с.

30. Гудман Р. Механика скальных пород. -М.: Стройиздат, 1987. -232 с.

31. Иофис М.А. Научные основы управления деформационными и дегазационными процессами при разработке полезных ископаемых. -М.: ИПКОН АН СССР, 1984. 230 с.

32. Бурмин Г.М., Васильчиков Т.П., Губин К.И. Технология разработки удароопасных месторождений Горной Шории // Горное давление и технология подземной разработки руд на больших глубинах: Сб. науч. тр. М.: ИПКОН АН СССР, 1990. - С. 133-134.

33. Методические указания по определению углов наклона бортов, откосов уступов и отвалов строящихся и эксплуатируемых карьеров. -Л.: ВНИМИ,1972. 165 с.

34. Прочность и деформируемость горных пород / Карташов Ю.П., Матвеев Б.В., Михеев Г.В., Фадеев А.Б. М.: Недра, 1979. -268 с.

35. Протосеня А.Г., Ставрогин А.Н. Пластическое деформирование горных пород при переменных скоростях нагружения // Аналитические методы и применение ЭВМ в механике горных пород: Сб. науч. тр. Новосибирск: ИГД СО АН СССР, 1982. - С. 155-158.

36. Методология расчета горного давления / Кузнецов С.В, Одинцев

37. Методы и средства решения задач горной геомеханики / Кузнецов Г.Н., Ардашев К.А., Филатов H.A. и др. М.; Недра, 1987. - 248 с.

38. Шуплецов Ю.П. Прочность и деформируемость массивов скальных пород около выработанных пространств: Дисс. докт. техн. наук. -Екатеринбург, 1998. 293 с.

39. Розовский Л.Б. Введение в теорию геологического подобия и моделирования. М.: Недра, 1966.

40. Зотеев О.В. Математические модели трещиноватых неоднородных скальных массивов // Повышение эффективности буровзрывных работ: Сб. научн. тр./ ИГД МЧМ СССР. Свердловск, 1986. - Вып. 82.1. C. 87- 90.

41. Ухов С.Б. Скальные основания гидротехнических сооружений. -М.: Энергия, 1975. -263 с.

42. Виттке В. Механика горных пород. М.: Недра, 1990. -439 с.

43. Мерзляков В.П. Особенности анизатропии и механические свойства трещиноватых скальных массивов // Аналитические методы и применение ЭВМ в механике горных пород: Сб. научн. тр. / ИГД СО АН СССР. Новосибирск, 1982. - С. 32 - 33.

44. Айтматов И.Т. Геомеханика рудных месторождений Средней Азии. Фрунзе: Илим. - 1987. - 246 с.

45. Баклашов И.В., Картозия Б.А. Механические процессы в породных массивах. М.: Недра, 1986. -272 с.

46. Беляев Е.В. Геомеханические преобразования и охрана подрабатываемых территорий. М.: Наука, 1988. - 191 с.

47. Мгалобелов Ю.Б. Исследование прочности береговых упоров арочных плотин учетом упругопластических деформаций // Сб. науч. тр. / Гидропроект им. С.Я. Жука. М., 1974. - № 33. - С. 18 - 24.

48. Выбор и обоснование расчетной схемы, разработка алгоритма решения задачи об устойчивости подготовительной выработки для различных типов забутованного материала: Отчет о НИР / МГИ. Рук. Ершов Л.В. № ГР 0183009931, инв. № 02830084447. М., 1983.-76 с.

49. Родионов В.Н., Сизов И.А., Цветков В.М. Основы геомеханики. -М.: Недра, 1986.-301 с.

50. Барях A.A., Константинова С.А., Асанов В.А. Деформирование соляных пород. Екатеринбург: УрО РАН, 1996. - 188 с

51. Турчанинов И.А., Иофис М.А., Каспарьян Э.В. Основы механики горных пород. Л.,: Недра, 1977.

52. Соколовский В.В. Статика сыпучей среды. М.: Физматгиз, 1960.-275 с.

53. Зотеев О.В., Осинцев В.А. Геомеханика: Учебное пособие. -Екатеринбург: УГГГА, 1997. 128 с.

54. Зенкевич О., Чанг И. Метод конечных элементов в теории сооружений и в механике сплошных сред. М.: Недра, 1974. - 240 с.

55. Зенкевич О., Морган К. Конечные элементы и аппроксимация. -М.: Мир, 1986.-318 с.

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

57. Сегерлинд Л. Применение метода конечных элементов. М.: Мир, 1979. -392 с.

58. Бенерджи П., Баттерфилд Р. Метод граничных элементов в прикладных науках. М.: Мир, 1984. - 494 с.

59. Крауч С., Старфилд А. Методы граничных элементов в механике твердого тела. М.: Мир, 1987. - 328 с.

60. Угодчиков А.Г., Хуторянский Н.М. Метод граничных элементов в механике деформируемого тела. Казань: КГУ, 1986. - 256 с.

61. Численные методы оценки устойчивости подземных сооружений: Сб. науч. тр. Апатиты, КФ АН СССР, 1988. - 111 с.

62. Юфин С.А., Харт Р.Д., Кюндалл П.А. Сравнительный анализ численных методов решения задач геомеханики // Энергетическое строительство. 1992. - № 7. С. 4 - 8.

63. Зотеев О.В., Осинцев В.А. Математическое описание объектов и процессов: Учебное пособие. Екатеринбург: УГГГА, 1998. - 124 с.

64. Румшинский Л.З. Математическая обработка результатов эксперимента. М.: Наука, 1971. - 192 с.

65. Зотеев О.В. Математическое моделирование неоднородных трещиноватых скальных массивов // Известия вузов. Горный журнал. Свердловск, 1986. 14 с. - Деп. в ВИНИТИ 22.12.86 N 6732 - В.

66. Зотеев О.В., Ножин А.Ф. Расчет напряженно-деформированного состояния скальных массивов горных пород // Устойчивость и технология формирования бортов и отвалов на глубоких карьерах: Сб. научн. тр./ИГД МЧМСССР, Свердловск, 1987. № 83. С. 112 - 117.

67. Демидович Б.П., Марон И.А. Основы вычислительной математики. М.: Физматгиз, 1966. - 664 с.

68. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. -М.: Наука, 1987. -600 с.

69. Икрамов Х.Д. Численное решение матричных уравнений. М.: Наука, 1984. 192 с.

70. Джордж А., Лю Д. Численное решение больших разреженных систем уравнений. М.: Мир, 1984. - 333 с.

71. Дэвис Дж. С. Статистический анализ данных в геологии: Пер. с англ. В 2 кн. М.: Недра, 1990 г.

72. Рац М.В., Чернышев С.Н. Трещиноватость и свойства трещиноватых горных пород. М.: Недра, 1970. - 164 с.

73. Шуплецов Ю.П. Об общности формирования очагов техногенных и природных землятрясений // Проблемы геодинамики, сейсмичности минерогении подвижных поясов и платформенных областей литосферы. Екатеринбург: Институт геофизики УрО РАН. - С. 236 - 238.

74. Методические указания по определению допустимых пролетов обнажений трещиноватых горных пород и размеров опорных целиков при подземной разработке рудных месторождений. М.: ИПКОН АН СССР, 1978.-92 с.

75. Методические указания по профилактике горных ударов с учетом геодинамики месторождений. Л.: ВНИМИ, 1983. - 118 с.

76. Геодинамическое районирование недр. Методические указания. -Л.: ВНИМИ, 1990.-127 с.

77. Маслов Н.Н., Павлищева Т.В., Фонарев П.А. Особенности внутреннего строения скальных массивов и их прочность на сдвиг по трещинам // Проблемы механики горных пород: Сб. науч. тр. М., Наука, 1987. -С. 231 -236.

78. Нейштадт Л.И., Пирогов И.А. Методы инженерно-геологического изучения трещиноватости горных пород. М.: Энергия, 1969. - 248 с.

79. Чернышев С.Н. Трещиноватость горных пород и ее влияние на устойчивость откосов. М.: Недра, 1984. - 111 с.

80. Вистелиус А.Б. Структурные диаграммы. М. - Л.: Изд-во АН СССР, 1958.- 158 с.

81. Газиев З.Г., Речицкий В.И. Вероятностная оценка надежности скальных массивов. М.: Стройиздат, 1985. - 104 с.

82. Зотеев В.Г. Выбор рациональной методики обработки результатов массовых замеров при изучении трещиноватости // Вопросы маркшейдерского дела. Белгород: ВИОГЕМ, 1962. - Ч. 2. - С. 60 - 67.

83. Зотеев В.Г., Морозов В.Н. Основы методики обработки массовых замеров трещиноватости на железорудных месторождениях различных генетических типов // Труды ИГД Минчермета СССР. Свердловск, 1975. - Вып. 47. С. 38 - 44.

84. Смирнов А.Ф. К методике обработки замеров трещиноватости // Труды ВНИМИ. Л., 1964. - Вып. 51. - С. 350.

85. Зотеев В.Г., Зотеев О.В., Костерова Т.К. Методика обработки массовых замеров трещиноватости на ЭВМ // Известия Уральского горного института. Сер.: Геология и геофизика. 1993. - Вып. 2. С. 28 - 34.

86. Введение в механику горных пород / под ред. Бока. М: Мир, 1983.-276 с.

87. Bandis S., Lumsden A.C., Barton N. Experimental studies of scale effects on the shear behaviour of rock joints // Int. J. Rock Mech. Sei. and Geomech. Abst. 1981. -Vol. 18.-P. 1 -21.

88. Никитин A.A., Старостина Т.Г., Рожкова Н.Ф. Методы оценки сопротивляемости сдвигу скальных массивов по сомкнутым трещинам // Извести ВНИИГ им, Б.Е. Веденеева. 1982. - Т. 156. С. 47-51.

89. Рекомендации по методике построения моделей скальных оснований плотин по параметрам сопротивляемости пород сдвигу: П 08 -83 / ВНИИГ. Л. - 1983. - 180 с.

90. Речицкий В.И., Эрлихман С.А. Современные методы определения прочности на сдвиг по трещинам // Геоэкология. 1997. - № 5. -С. 102-114.

91. Кочарян Г.Г., Кулюкин A.M. Исследование закономерностей обрушения подземных выработок в горном массиве блочной структуры при динамическом воздействии. Часть 2. О механических свойствах межблоковых промежутков // ФТПРПИ. 1994. - № 5. - С. 27 - 37

92. Власов А.Н., Рогозинский A.B., Ухов С.Б. Определение угла дилатансии в скальных породах при сдвиге по трещинам // Проблемы механики горных пород: Труды XI Российской конференции по механике горных пород. С. Пб., 1997. - С. 87 - 92.

93. Ngo D., Scordelis А. С. Finite element analysis in reinforced concrete Beam// J. of American concrete Inst. 1967. - Vol. 64. - N 3. - P. 152-163.

94. Goodman R.E., Tayler R.L., Brekke T.L. A model for the Mechanic of jointed rock // J. of the Soil Mechanic and Foundation Division. 1968. -Vol. 94. - N 3. - P. 637-659.

95. Речицкий В.И. Результаты экспериментальных исследований жесткости трещин в скальных породах // Геоэкология. 1998. - № 2. - С. 88-99.

96. Бенявски 3. Управление горным давлением. М.: Мир, 1990. -254 с.

97. Руппенейт К.В., Тарасова И.В. Способы учета влияния различных классов трещин на деформируемость массивов горных пород // Современные проблемы механики горных пород: Мат. 4 Всесоюзной конференции по механике горных пород. Л.: Наука, 1972. - С. 11 - 20.

98. Роза С.А., Зеленский Б.Д. Исследование механических свойств скальных оснований гидротехнических сооружений. М.: Энергия, 1967. -391с.

99. Глушко В.Т., Борисенко В.Г. Инженерно-геологические особенности железорудных месторождений. М.: Недра, 1978. - 253 с.

100. Кулатилаки П. Влияние трещин на прочность и деформатив-ность массивов скальных пород // Энергетическое строительство. -1992. № 8. С. 25-28.

101. Ким Д.Н. Влияние структуры на сдвиговую прочность массива и определение расчетных механических характеристик // Горное давление, сдвижение горных пород и методика маркшейдерских работ: Сб. науч. тр. Л.: ВНИМИ, 1969. - С. 568-585.

102. Протопопов И.И. Основные свойства и особенности напряженно-деформированного состояния блочных пород // Свойства горного массива и управление его состоянием: Сб. науч. тр. СПб: ВНИМИ, 1991. -С. 79-84.

103. Батугин С.А. Анизатропия массива горных пород. Новосибирск: Наука, 1988. - 86 с.

104. Brady В. H. G., Brown Е.Т. Rock Mechanics For Underground Mining. London: George Allen and Unwin, 1985. - 587 p.

105. Свобода Г., Штумволь M., Ито Ф. Численное моделирование повреждений в трещиноватых скальных породах с учетом шероховатости трещин // Энергетическое строительство. 1992. - № 12. С. 9-17.

106. Rowray A.L., Goodman R.E. Finite element analysid of crack initiation in a block rock model experiments // Rock Mechanics. 1979. - Vol. 4. -№4. - P. 203-233.

107. Трофимов B.A. Численное моделирование взаимодействия разлома и очистных выработок по мере развития горных работ // Физические основы прогнозирования разрушения горных пород: Сб. науч. тр. Фрунзе: Илим, 1985. - С. 71.

108. Назарова Л.А. Напряженное состояние наклонно-слоистого массива горных пород вокруг выработки // ФТПРПИ. 1985. - № 2. - С. 33-37.

109. Carol J., Alonso Е. A new jointed element for the analysis of fracture rock // 5th Congr. Int. Soc. Rock Mech. Melburn, Australia. - 1983.-P. 147-151.

110. Курленя M.B., Попов C.H. Теоретические основы определения напряжений в горных породах. Новосибирск: Наука, 1983. - 97 с.

111. Влох Н.П. Управление горным давлением на подземных рудниках. М.: Недра, 1994. - 208 с.

112. Пискарев B.K. Интерпретация деформаций керна при измерении напряжений методом разгрузки в трещиноватых породах // Методология измерения напряжений в массиве горных пород: Сб. науч. тр. / ИГД СО АН СССР. Новосибирск, 1978. - С. 12 - 16.

113. Влох Н.П., Зубков A.B., Феклистов Ю.Г. Совершенствование метода щелевой разгрузки // Диагностика напряженного состояния породных массивов: Сб. науч. тр. / ИГД СО АН СССР. Новосибирск, 1980. - С. 30 - 35.

114. Влох Н.П., Зубков A.B., Феклистов Ю.Г. Метод частичной разгрузки на большой базе // Диагностика напряженного состояния породных массивов: Сб. науч. тр. / ИГД СО АН СССР. Новосибирск, 1980. -С. 37 - 42.

115. Влох Н.П., Шуплецов Ю.П. Определение напряженного состояния массива горных пород по приращениям напряжений вокруг выработанных пространств // Напряженное состояние породных массивов: Сб. науч. тр. / ИГД СО АН СССР. Новосибирск, 1978. - С. 98 - 101.

116. Попов С.Н. Точная расчетная схема метода параллельных скважин // Диагностика напряженного состояния породных массивов: Сб. науч. тр. / ИГД СО АН СССР. Новосибирск, 1980. - С. 3 - 8.

117. Колмаков В.Д. Техника экспериментального измерения напряжений методом параллельных скважин в скальных породах // Исследования напряжений в горных породах: Сб. науч. тр. / ИГД СО АН СССР. -Новосибирск, 1985. С. 120 - 126.

118. Турчанинов И.А., Панин В.И., Кононенко B.C., Иванов В.И. и др. Определение напряжений в массиве скальных пород // Диагностика напряженного состояния породных массивов: Сб. науч. тр. / ИГД СО АН СССР. Новосибирск, 1980. - С. 42 - 48.

119. Журин С.Н., Колесников В.И., Стрельцов В.И. Геомеханический мониторинг обводненных массивов. М.: НИА-Природа, 1997.- 188 с.

120. А. с. 1772594 СССР. Способ определения напряжений в исследуемом объекте / Зубков A.B., Феклистов Ю.Г., Крутиков A.B., Зотеев О.В. Опубл. в Б.И. - 1992. - № 40.

121. Методика геофизического прогноза удароопасности участков угольных пластов и рудных залежей. М.: Недра, 1980. - 88 с.

122. Л. М. Лазаркевич. Определение направлений действия главных нормальных напряжений методом электрометрии // Геофизические основы контроля напряжений в горных породах: Сб. науч. тр. Новосибирск, ИГД СО АН СССР, 1983. - С. 15 - 19

123. Проскуряков В.М., Бляхман A.C. Исследование напряженного состояния удароопасных пород сейсмическим методом // Геофизические основы контроля напряжений в горных породах: Сб. науч. тр. Новосибирск, ИГД СО АН СССР, 1983. - С. 52 - 55.

124. Кириченко Г.С., Малетин Л.В., Пахомов В.А. Способ определения напряженно-деформированного состояния закладочного массива //

125. Исследования проявления горного давления и технологии подземной разработки руд на больших глубинах: Сб. науч. тр. / ИПКОН АН СССР. -М., 1983. С. 15-30

126. Методы и средства контроля состояния и свойств горных пород в массиве/ Ватолин Е.С., Черняков А.Б., Рубан А.Д., Потапов A.M. М.: Недра, 1989. - 173 с.

127. Методические рекомендации по определению напряжений массива пород и оценке удароопасности горных выработок методом электрометрии. Свердловск: ИГД МЧМ СССР, 1986.-47 с.

128. Проведение экспериментально проверки и внедрение комплекса мер: предотвращающих горные удары на горнодобывающих предприятиях рудной промышленности: Отчет о НИР / ИГД МЧМ СССР. Рук. Липин Я.И. Свердловск, 1990.

129. Zoubkov A.V., Lipin Y.I. Stress state in the Urals crust and tectonic evolution // Rock Stress: Proc. Int. Symp. on Rock stress. Rotterdam: Balkema, 1997. - P. 403 - 408.

130. Кунин И.А., Соснина Э.Г. Напряженное состояние упругого слоя в среде // ФТПРПИ. 1973. - № 2. - С. 23 - 27.

131. Савченко С.Н., Козырев A.A., Мальцев В А. Напряженное состояние пород блочного строения // ФТПРПИ. 1994. - № 5. С. 38-47.

132. Серяков В.M. Об одном подходе к расчету напряженно-деформированного состояния массива горных пород в окрестности выработанного пространства // ФТПРПИ. 1997. - № 2. - С. 14-21.

133. Обоснование подготовки и отработки рудного тела № 2 Узель-гинского месторождения: Отчет о НИР / ИГД УрО РАН. Рук. Волков Ю.В., Сашурин А.Д. Екатеринбург, 1997. - 102 с.

134. Ветров C.B., Одинцев В.Н., Слоним М.Э. Сопоставление технологических схем разработки Яковлевского месторождения по фактору напряженно-деформированного состояния массива горных пород и закладки //тамже. С. 110-124.

135. Бушканец С.С., Скринник В.В. Экспериментальное изучение деформируемости галечниковых грунтов // Известия ВНИИГ Л.: Энер-гоатомиздат,1978. - Т. 122. - С. 19 - 23.

136. Лофицкий В.Н., Пахомов O.A. Результаты экспериментальных исследований деформируемости крупнообломочного грунта под большими нагрузками // Известия ВНИИГ Л.: Энергоатомиздат,1973. - Т. 102.-С. 137-144.

137. Нечипорович A.A., Рассказов Л.Н. Сопротивление крупнообломочных грунтов сдвигу // Гидротехническое строительство. 1969. - № 8.-С. 21 -26.

138. Павчич М.П., Пахомов O.A. Экспериментальные исследования деформируемости смесей аллювиальных грунтов под большим давлением // Гидротехническое строительство. 1979. - № 4. - С. 20 - 24.

139. Розанов H.H. Плотины из грунтовых материалов. М.: Строй-издат, 1983.-295 с.

140. Научное обеспечение опытно-промышленного заполнения Главного карьера Высокогорского ГОКа хвостами МОФ // Отчет о НИР / УГГГА. Рук. Зотеев О.В. Екатеринбург, 1995. - 36 с.

141. Курленя М.В., Серяков В.М. О методе расчета напряженно-деформированного состояния горных пород с учетом контакта кровли и почвы выработанного пространства. ФТПРПИ. - 1997. - № 2.

142. Курленя М.В., Серяков В.М. Геомеханическое состояние горного массива с учетом последовательности отработки // Проблемы геотехнологии и недроведения. Докл. междун. конф. (Мельниковские чтения). - Екатеринбург: ИГД УрО РАН, 1988. - Т. 3. - С. 74 - 79.

143. Совершенствование параметров конструктивных элементов систем разработки, опытно-промышленная проверка и внедрение их на шахтах Гороблагодатского рудоуправления // Отчет о НИР / ИГД ММ СССР. Рук. Зубков А.В О.В. Свердловск, 1990. - 110 с.

144. Определение напряжений приконтурного массива горных пород, в крепи и разработка технических решений по проходке стволов ш. "Центральная": Отчет о НИР / ИГД МЧМ СССР. Рук. Влох Н.П., Боликов В.Е. Свердловск, 1989. - 62 с.256

145. Инструкция по выбору и расчету параметров крепи горизонтальных горных выработок и управлению устойчивостью камер на шахте "Центральная" Джездинского РУ. ИГД ММ СССР. Свердловск, 1991. - 69 с.1. СПРАВКА

146. Разработанные О.В. Зотеевым метод и комплекс программ для моделирования напряженно-деформированного состояния (НДС) скальных массивов широко применялись при расчетах устойчивости элементов систем разработки на подземных рудниках Урала и Казахстана.

147. Разработанный О.В. Зотеевым метод моделирования НДС элементов систем разработки также использовался при обосновании параметров и порядка отработки запасов блока № 5 на шахте "Западный Кара-жал" АО "Атасуруда".

148. Переданный О.В. Зотеевым ИГД УрО РАН в 1997 г. комплекс программ для расчета НДС скальных массивов широко используется лабораторией горного давления при выполнении текущих НИР.

149. Ректор Уральской государственнойгорно-геологическо профессор1. И.В. Дементьев