УДК 622.831.1
А.А. Неверов, С.А. Неверов
ПРОГНОЗНАЯ ОЦЕНКА НАПРЯЖЕННОГО СОСТОЯНИЯ ПОРОД ПРИ ОТРАБОТКЕ МАКМАЛЬСКОГО МЕСТОРОЖДЕНИЯ ЗОЛОТА
КРАТКАЯ ХАРАКТЕРИСТИКА МЕСТОРОЖДЕНИЯ
Макмальское золоторудное месторождение расположено в Кыргызской Республике на юго-западе гор Чаарташ. Поле месторождения сложено окварцованными, скарнированными известняками, скарнами, метасоматитами, относящимися к золотокварцевой формации малосульфидного типа. Рудная минерализация весьма разнообразна и представлена в основном пиритом, пирротином, арсенопиритом, кварцем, полевыми шпатами и пироксеном. Содержание золота в руде составляет 3.5 - 7.0 г/т.
Слагающие месторождение породы, представленные мраморами, гранитами и песчаниками обладают высокими прочностными и упругими свойствами с коэффициентом крепости по шкале проф. М.М Протодьяконова / = 14-19. Объемный вес руды - 2.73 т/м3, вмещающих пород - 2.7 т/м3. В морфологическом строении залежь имеет сложную форму, мощностью от 5 до 60 м. Рудная зона распространяется на глубину около 600 м с вертикальным падением.
В верхней границе месторождения локализованы две самостоятельных линзы - Северная и Южная (рис. 1), которые ниже гор. 2415 м образуют единую рудную залежь. Для визуализации
изображения на рис. 1б выработанное пространство карьера и открытых камер Северной, Южной линз представлено моделью твердого тела.
Северная линза имеет вертикальное падение и распространяется на глубину 250 м. В настоящее время она отработана до горизонта 2370 м. Южная линза имеет простую симметричную форму, параметры которой составляют: мощность 46 - 48 м, длина достигает 130 м.
Горно-складчатая область расположения месторождения предопределила его сложную структуру и неоднозначность геомеханической обстановки, которая характеризуется с одной стороны -наличием высокомодульных прочных и устойчивых пород, с другой - возможностью проявлений динамических явлений вследствие ведения работ под открытым выработанным пространством камер, нависающими породными консолями и при отработке разделительного рудо-породного целика. Отметим, что предприятие «Макмал» находится в 9-ти бальной сейсмогенерирующей зоне с повторяемостью сильных землетрясений 2 раза в 1000 лет.
В настоящее время на предприятии ведется отработка подкарьерных запасов весьма мощного Южного рудного тела. В соответствии с проектом месторождение вскрыто штольнями с высотою
а б
Рис. 1. Развитие горных работ на Макмальском золоторудном месторождении: а) общий вид поверхности месторождения; б) объемная твердотельная модель выработанного пространства (открытых камер и выработок) при отработке Южного рудного тела; 1 - карьер; 2 - целик; 3, 4 - соответственно выработанное пространство Северной и Южной линз; 5, 6 - породные консоли над отработанными Северной и Южной камерами
этажа 40 м и подэтажа - 20 м. Система разработки для выемки мощных рудных тел - этажное и по-дэтажное обрушение, для отработки рудных залежей мощностью до 10 м - подэтажных штреков.
На современном этапе развития горных работ технология выемки подкарьерных запасов Южного рудного тела была нарушена. В настоящее время отработка залежи ведется камерной системой разработки с камерно-целиковым порядком выемки и оставлением целиков. В результате под дном карьера образовалось выработанное пространство объемом над Северной и Южной линзами 211.2 и 94.2 тыс. м3 с нависающими породными консолями объемом 70.0 и 27.0 тыс. м3 соответственно (рис. 1). Последние уходят в борта карьера и имеют высоту от 20.0 до 50.0 м.
В центральной части рудного тела оставлен рудо-породный целик мощностью до 20 м, содержащий линзу забалансовой руды. Наличие этого целика и физико-механические свойства пород месторождения обусловили достаточно высокую устойчивость выработанному пространству и горным конструкциям в течение длительного периода времени (около 5 лет). Вместе с тем, опасность внезапных обрушений отдельностей нависших породных консолей, бортов открытых камер и оставленного рудо-породного целика, разделяющего Северную и Южную линзы, постоянно возрастает. Последнее обстоятельство создает аварийную ситуацию на предприятии и требует незамедлительной остановки горных работ до разработки специальных мероприятий, регламентирующих безопасную технологию добычи и порядок выемки запасов.
Целью выполненного исследования является геомеханическая оценка сложившейся горнотехнологической ситуации на Южном рудном теле Макмальского золоторудного месторождения и определение возможности его перспективного и безопасного освоения.
ПОСТАНОВКА ЗАДАЧИ
Моделирование механического поведения массива пород для выделенных условий проводилось в предположении объемного напряженного состояния методом конечных элементов (МКЭ) [1]. Руда и слагающие месторождение породы имитируется упругими, изотропными, линейно-деформируемыми, однородными материалами.
Недостаточная изученность исходного силового поля напряжений на месторождении предопределила постановку и решение задачи в двух вариантах.
1. Исходное напряженное состояние массива соответствует гипотезе Динника:
где - коэффициент бокового
распора, у - плотность пород, g - ускорение свободного падения, Н - расстояние до земной поверхности, ^ - коэффициент поперечных деформаций.
2. Исходное поле напряжений определяется тектоническими силами:
^ииии*
где Хд - коэффициент бокового давления.
И.Т. Айтматов для рудников Средней Азии предложил эмпирические зависимости, отражающие изменение величин горизонтальных напряжений с глубиной (Н) [2]:
' + сгЕ = аа + ад&ав,
<ТХ = 5.0 + = 2.0 +
В соответствии с этими данными, в качестве краевых условий второго варианта задачи, приняты следующие значения тектонического поля напряжений:
1^-2Л;.'#
На рис. 2а приведена расчетная схема и основные параметры к определению НДС массива. В центральной части модели заложена конфигурация целиков, открытых камер и консолей, которые отражают сложившуюся обстановку на предприятии. Глубина горных работ от дна карьера до горизонта +2370 м составляет 130-140 м.
На исследуемой базе приняты следующие граничные условия (рис. 2а):
На плоскости АВСБ -
х> 2}= г>=й
ЕР ОН- ' '
^ л, ?}=(.
БСОН- ' '
л ^=-чт
АВЕЕ-
Ту/0.$2)=тв/ах.5}=6
ВСОЕ-
[л/£рл л *}=луф
Я- */=*
АБНЕ-
# х}=1Г^О> # х}=С где оу, ах, а2 - вертикальная и горизонтальные нормальные компоненты тензора напряжений, тху, Тхг, Ху2 - касательные компоненты тензора напряжений, иу, их, и2 - вертикальная и горизонтальные компоненты вектора перемещений. Для учета веса горных пород окружающего массива в моделях задавалась объемная сила у?.
Исходные данные для расчетов: модуль Юнга Е = 60000 МПа, коэффициент Пуассона ц = 0.3, плотность соответственно породы и руды у= 2700 и 2730 кг/м3. Шаг дискретизации расчетной области в районе влияния открытых камер, рудо-
Рис. 2. Моделирование НДС массива пород: а) расчетная схема; б) дискретизация геомеханической модели
породного целика и консолей принят от 1.0 до 2.0 м, на внешнем контуре модели - до 12.0 м (рис. 2б). Общее количество элементов
, представленных 4-ех узловыми тетраэдрами и узлов объемной сетки, соответственно составило 2830000 и 557000.
Верхняя граница расчетной области соответствует земной поверхности. На границе контакта двух материалов (порода-руда) предполагается выполнение условия жесткого сцепления (равен-
„1 £ ство компонент вектора перемещении и,-
1 а 1 я ТСЧ. = Ин, — Т^,). Поверхности граничных
нагрузок удалены от области горных конструкций на расстояние, которое не оказывает влияние на природное поле напряжений (не менее чем на трех кратный размер наибольшего пролета выработанного пространства).
В качестве исследуемых областей принимались разрезы вдоль плоскостей: 0x7 - по горизонту + 2445 м; 0ух - поперек Северной и Южной камер (линия 9820); 0у7 - разрез соответственно
по рудо-породному целику и Северной камере (линии 10570 и 10590). Результаты расчетов представляются в виде картин изолиний напряжений. За максимальное главное напряжение принималось оь за минимальное - о2, растяжение - со знаком «минус».
АНАЛИЗ РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЙ
Особенности распределения полей напряжений 01, и о2 (МПа) в массиве пород в условиях гравитации приведены на рис. 3 - 6.
Обобщение результатов расчета данного варианта показало, что области наибольших значений напряжений 01 локализуются, в основном, в целике горизонта +2460 м и ниже, породных консолях и в бортах открытых камер (от 3 до 7-8 МПа).
Рудо-породный целик практически по всей высоте от гор. +2480 до гор. +2396 м испытывает двухосное сжатие (о2 = 0 МПа, максимальное и промежуточное напряжения являются положи-
Рис. 3. Характер распределения напряжений в массиве пород на гор. +2445м (I вариант):
а - а2; б -
Рис. 4. Характер распределения напряжений в массиве пород вдоль линии 9820 (I вариант):
а - а2; б - о1
а б
Рис. 5. Характер распределения напряжений в массиве пород вдоль линии 10570 (I вариант):
а - о2; б - о1
тельными). На участке целика в районе гор. +2445 м, где его размеры существенно снижаются, наблюдается незначительная зона одноосного растяжения о2 = -1 МПа. В остальных элементах конструкции наблюдается условие объемного сжатия. Ниже горизонта +2445 м борта камер со стороны целика находятся в состоянии близком к двухосному сжатию (о2 = 0, 01 = 1-2 МПа).
Из анализа следует, что наиболее нагруженными по напряжениям 01 и ттах являются участки породной консоли и стенок обнажений со стороны Северной камеры, однако уровень этих напряжений значительно ниже, чем пределы прочности вмещающих пород и руды на растяжение. Поэтому разрушений на этих участках наблюдаться не будет.
Аналогичная ситуация отмечается в районе Южной камеры (гор. +2445 м). В области Северной камеры уровень напряжений выше, чем в окрестностях Южной, что является следствием больших ее размеров. Породные консоли, как и целик, находятся в двухосном напряженном со-
стоянии, с постепенным переходом в борта карьера к объемному сжатию.
Концентрация сжимающих напряжений ох -зоны опорного давления - наблюдается в окрестности вертикальных границ открытых камер со стороны массива пород и в горных конструкциях (целиках, днищах, выработках) на глубоких горизонтах.
В целом, оценивая уровень действующих напряжений, в рассматриваемой задаче, отметим, что максимальными являются вертикально ориентированные напряжения и их абсолютные значения не являются критическими с позиции разрушения пород. Последнее обстоятельство не оказывает существенного влияния на устойчивость горных конструкций и безопасность ведения подземных работ.
На рис. 7 - 10 показано распределение действующих напряжений, сформированных отработкой залежи в условиях тектонического поля напряжений.
Рис. 6. Характер распределения напряжений в массиве пород вдоль линии 10590 (I вариант):
а - а2; б - сг;
Рис. 7. Напряжения в массиве пород на горизонте +2445м в условиях тектоники (II вариант):
а - 02; б - 01
В сравнении с предыдущей задачей, здесь на- наблюдаются в целике (01 = 9-12, ттах = 4-7 МПа)
блюдается заметный рост напряжений за счет тек- и консолях, как с северной (01 = 15-18, ттах = 6-9
тонических сил. Результаты анализа позволили МПа), так и с южной стороны (01 = 12-15, ттах =
установить места концентраций, где наибольшие 6.0-7.5 МПа).
абсолютные значения максимальных напряжений Наличие зон растяжение в целике 02 = 0 ^ -2
а
б
Рис. 9. Напряжения в массиве пород вдоль линии 10570 в условиях тектоники (II вариант):
а - 02; б - 01
б
Рис. 10. Напряжения в массиве пород вдоль линии 10590 в условиях тектоники (II вариант):
а - 02; б - 01
МПа свидетельствует о возможной потере его устойчивости. Однако полученные напряжения не являются запредельными, в виду высоких значений пределов прочности пород на сжатие и растяжение.
В породной консоли с северной стороны происходит постепенный переход от двухосного сжатия (а2 ~ 0 МПа) к объемному в направлении борта карьера. Поэтому консоль и борта карьера в районе открытых камер будут находиться в состоянии средней и ниже средней устойчивости.
Сопоставляя полученные результаты с напряженным состоянием в условиях гравитации, отметим, что качественная и количественная картины изменились в худшую сторону со значительными диапазонами различия. Так рост напряжений в целике составил более чем в 3 раза.
Наиболее благоприятная ситуация в рассматриваемом варианте характерна для стенок камер,
тия. В подобных случаях устойчивость пород очевидна.
Анализ результатов исследования НДС массива на Южном рудном теле месторождения при вариации параметров исходного силового поля земли показал, что величины напряжений, полученные в ходе расчетов, не превышают прочностных характеристик пород и не являются предельными, вследствие чего опасность обрушений элементов разработки снижается.
УСТОЙЧИВОСТЬ КОНСТРУКТИВНЫХ ЭЛЕМЕНТОВ ГОРНЫХ РАБОТ
Устойчивость элементов разработки оценивалась по двум критериям.
1. По критерию трещинообразования (статистический подход), в основе которого лежит длительная прочность и структурная нарушенность породы в массиве [3]:
гг, <1 К с с
которые испытывают состояние трехосного сжа-
Таблица 1. Прочностные свойства выделенных горных пород с учетом Кс
Тип породы Ф, град С, МПа асж, МПа ар, МПа
Мраморы слабоскарнированные (I - тип) 50 7.5 41.4 5.4
Мраморы скарнированные(11 - тип) 53 9.3 54.7 6.3
где 0д - напряжения в массиве, МПа; Кс - показатель учитывающий длительную нарушенность породы, д.ед.; 0сж - предел прочности породы на сжатие, МПа.
2. По критерию запаса прочности пород (Кзап), выраженному из уравнения прямой напряжений Мора (считается, что наиболее разрушительными являются сдвиговые напряжения):
т =
где т - предел прочности породы по касательным напряжениям, МПа; С - сцепление, МПа;0 - нормальное напряжение, действующее на площадке разрушения, МПа; ф - угол внутреннего трения.
Коэффициент запаса прочности пород [4]:
Луди------------------------21, Д. ед.
Ч-Ч
где 01 и 02 - соответственно максимальные и минимальные значения напряжений, действующих в массиве пород, МПа.
Условие Кзап < 1 соответствует процессу разрушения пород
Частичная неопределенность исходной информации предопределили рассмотрение показателя устойчивости для двух видов наименее прочных пород, слагающих Южное рудное тело (табл. 1) с учетом коэффициента структурного ослабления Кс = 0.6.
Результаты расчета устойчивости представлены в табл. 2.
В условиях гравитационного поля напряжений, учитывая. что массив пород испытывает значительно меньшие нагрузки, устойчивость горных конструкций не вызывает сомнений (Кзап составляет более 2.0). Вполне обоснованно предположить, что исследуемые элементы камер и целиков будут подвержены большим значениям горного давления при действии горизонтальных сил. Исходя из этого, к более детальной проработке была принята ситуация ведения горных работ в условиях прогнозной тектоники.
Консоли Северной и Южной камер при действии повышенных горизонтальных напряжений находятся в пограничном состоянии, устойчивость которых со временем потенциально снижается. Это обстоятельство может спровоцировать возможные внезапные обрушения пород. Зависшая консоль над Северной камерой характеризуется большей площадью обнажения по сравнению с Южной (8С = 1665 м2 против 8ю = 744 м2) и меньшей прочностью. Снижение устойчивости пород также наблюдается в бортах камер в интервале горизонтов 2460 - 2445 м.
Степень прочности рудо-породного целика определяется взаимным влиянием выработанных пространств Северной и Южной камер и характеризуется, на данном этапе ведения работ, как достаточно высокая.
Вместе с тем, дальнейшая разработка залежи без коренного изменения технологии очистной выемки может привести к полной потере устойчи-
Таблица 2. Коэффициенты запаса прочности горных пород в условиях тектонического поля напряжений
Наименование участка Кзап
Тип породы - I Тип породы - II
Борт Северной камеры:
гор. 2460 м 1.4 1.6
гор. 2445 м 1.3 1.5
гор. 2415 м 1.5 1.6
гор. 2396 м 1.4 1.5
гор. 2370 м 1.6 1.8
Борт Южной камеры:
гор. 2460 м 1.3 1.4
гор. 2445 м 1.3 1.4
гор. 2415 м 1.8 2.0
гор. 2396 м 1.7 1.9
гор. 2370 м 2.1 2.3
Консоль:
Северная 1.2 1.2
Южная 1.2 1.3
Целик:
гор. 2460 м 1.7 1.7
гор. 2445 м 1.2 1.2
гор. 2415 м 1.4 1.5
гор. 2396 м 1.7 1.9
гор. 2370 м 1.8 2.0
вости и обрушению целика в динамической форме. Коэффициент запаса прочности наиболее уязвимых участков массива не гарантирует сохранность выработанного пространства в случае значительных сейсмических событий в регионе, в том числе, массовых взрывов больших масштабов.
ПРЕДПОСЫЛКИ БЕЗОПАСНОГО
ОСВОЕНИЯ МЕСТОРОЖДЕНИЯ
Исходя из сложившихся горнотехнических условий на Южном рудном теле, добычные работы следует приостановить до погашения открытых камер закладкой и образования предохранительных породных подушек. Безопасную отработку руды, в оставленных целиках штолен № 7 и 6, возможно технически осуществить системой этажного обрушения с отбойкой руды на зажатую среду (рис. 11). Эта геотехнология достаточно широко известна в мировой практике.
Очистная выемка ведется сплошным фронтом. Руду отбивают вертикальными или слабонаклонными слоями взрыванием скважинных зарядов на обрушенную породу (в нашем случае, на сухую закладку). Выпуск отбитой руды ведут, через дуч-ки непосредственно на скреперный штрек и далее через полок грузится в вагоны. Очистное пространство не поддерживается и заполняется по мере выпуска руды налегающими породами. При всей простоте эта технология для достижения удовлетворительных показателей извлечения руды из недр обусловливает ряд требований к про-
п-р
Рис. 11. Технологическая схема отработки целиков на гор. 2445 - 2415 м системой этажного обруше ния: 1 - дучки; 2 - пудоппиемные воронки; 3 - отбитая, пуда: 4 - налегающие вмешаюшие породы; цессам отбойки и выпуска.
ЗАКЛЮЧЕНИЕ
Под влиянием временного фактора и атмосферного воздействия устойчивость сложившейся горной конструкции будет постоянно снижаться, создавая аварийность и возможную остановку дальнейшего освоения месторождения. Произведенные расчеты геомеханического состояния выработанного пространства по Южному рудному телу, на современном этапе добычи, позволяют оценить горнотехническую ситуацию на предприятии в меру безопасной, перспективной и пригодной для ведения горных работ с выполнением оп-
ределенных технологических ограничений и профилактических мероприятий, предусмотренных ЕПБ Кыргызской Республики.
Учитывая высокую сейсмичность региона, постоянные взрывные работы, очистную выемку на месторождении до выполнения требований, направленных на предупреждение негативных последствий неуправляемого поведения массива горных пород (высокая концентрация напряжений, воздушные удары при обрушении отдельностей обнажений и др.) и до образования предохранительных породных подушек следует остановить.
СПИСОК ЛИТЕРАТУРЫ
1. Зенкевич О. Метод конечных элементов в технике. - М.: Мир, 1975. - 538 с.
2. Айтматов И.Т. Методы и результаты изучения напряженного состояния массивов и создание на их основе эффективных способов управления горным давлением при подземной разработке руд / Ахматов В.И. и др. // ФТПРПИ.-1987.-№4.
3. Тектонические напряжения в земной коре и устойчивость горных выработок. / Турчанинов И. А., Марков Г.А., Иванов В.И., Козырев А.А. - Л., «Наука», 1978., 256 с.
4. КазикаевД.М. Геомеханика подземной разработки руд / Учебник для вузов, - М.: Изд-во МГГУ, 2005.
□ Авторы статьи:
Неверов Сергей Алексеевич
- канд.техн.наук, научн.сотр. ИГД СО РАН Тел. (383) 217-09-52, е-та^: пеуегоуаа 79@mail.ru
Неверов Александр Алексеевич
- канд.техн.наук, научн.сотр. ИГД СО РАН, е-таД: пеуегоуаа 79@mai1.ru