Авторефераты по всем темам >>
Авторефераты по земле
На правах рукописи
рукописи АНДРИЕВСКИЙ Александр Порфирович
Физико-техническое обоснование параметров разрушения горного массива взрывом удлиненных зарядов
Специальность 25.00.20 - Геомеханика, разрушение горных пород, рудничная аэрогазодинамика и горная теплофизика
Автореферат диссертации на соискание ученой степени доктора технических наук
Новосибирск - 2009
Работа выполнена в Институте химии и химической технологии (ИХХТ СО РАН).
Научный консультант: доктор физико-математических наук Шер Евгений Николаевич
Официальные оппоненты: доктор технических наук, профессор Еременко Андрей Андреевич;
доктор технических наук, профессор Ермолаев Александр Иванович;
доктор технических наук, профессор Коростовенко Вячеслав Васильевич.
Ведущая организация - Московский государственный горный университет
Защита диссертации состоится л 11 декабря 2009 г. в 11 часов на заседании диссертационного совета Д 003.019.01 при Институте горного дела СО РАН (630091, г. Новосибирск, Красный проспект, 54)
С диссертацией можно ознакомиться в библиотеке ИГД СО РАН Автореферат разослан л____ ___________ 2009 г.
Ученый секретарь диссертационного совета доктор технических наук Попов Н.А.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы исследования. Взрывные работы являются основным методом отделения породы от массива в горном деле. Как отмечают ведущие ученые, наблюдается прямая зависимость уровня развития горной промышленности от уровня организации процессов разрушения горных пород при добыче полезных ископаемых. Объемы ведения взрывных работ в Советском Союзе считались наибольшими в мире. Потребление взрывчатых веществ (ВВ) для нужд горного производства приближалось к 2 млн. т в год. К середине 90-х годов прошлого века их потребление в России уменьшилось в два раза. В настоящее время с увеличением промышленного производства происходит рост потребления ВВ. При таких масштабах взрывных работ даже незначительное общее снижение удельного расхода ВВ и повышение выхода дробленной горной массы с 1 м скважины (шпура) приводит к большому экономическому эффекту в масштабах горнодобывающих отраслей России.
Во многих случаях неоправданный перерасход ВВ, а также низкие показатели выхода дробленной горной массы с 1 м скважины (шпура) обусловлены недостаточной теоретической обоснованностью применяемых методик расчета параметров БВР. Практическая проблема надежного обеспечения высоких показателей БВР лежит в створе важнейших задач горной науки о разрушении горных пород: определение механизма разрушения горной породы и изучение явлений, сопровождающих процессы разрушения и отделения горной породы от горного массива.
Таким образом, физико-техническое обоснование параметров разрушения породного массива удлиненными зарядами представляет актуальную проблему, имеющую важное научное и практическое значение.
Цель работы - получение закономерностей развития зон нарушенности и установление факторов, влияющих на процессы трещинообразования для физико-технического обоснования параметров взрывного разрушения горного массива удлиненными зарядами.
Идея работы заключается в использовании установленных закономерностей по определению формы, размеров, динамики развития зон разрушения и комплексно учитывающих физические свойства взрываемых массивов скальных пород, характеристики ВВ, геометрические размеры и места инициирования удлиненных зарядов.
Задачи исследований.
1. Проведение экспериментальных исследований параметров зон смятия и трещинообразования на реальных массивах скальных пород с последующей математической обработкой данных экспериментов для подтверждения и уточнения аналитических исследований.
2. Разработка теоретических моделей для расчета параметров зон смятия и трещинообразования при разрушении породного массива взрывом удлиненных зарядов, и согласование их с полученными экспериментальными данными.
3. Оценка влияния нарушенности массива на параметры зоны трещинообразования.
4. Определение влияния временных параметров процесса взрыва на развитие зон разрушения при различных физических свойствах взрываемых массивов скальных пород, характеристик ВВ, параметров и места инициирования удлиненных зарядов.
5. Разработка теоретических основ расчета рациональных параметров БВР при проведении горных выработок и скважинной отбойке полезных ископаемых.
Методы исследований.
Анализ и обобщение научно-технической информации, применение результатов теории сопротивления материалов, математическое моделирование, промышленные эксперименты и испытания, корреляционный и регрессионный анализ результатов экспериментов.
Научные положения, выносимые на защиту.
1. Определяющие параметры разрушения горного массива - радиус зоны смятия и предельно возможный радиус зоны трещинообразования, образующиеся при воздействии на монолитный скальный массив взрыва удлиненного заряда, описываются закономерностями в виде произведения степенных функций от диаметра заряда, плотности заряжания, скорости детонации ВВ, пределов прочности массива на сжатие и срез.
2. В трещиноватом массиве предельный радиус зоны трещинообразования увеличивается обратно пропорционально коэффициенту структурного ослабления, являющемуся функцией отношения предельно возможного радиуса зоны трещинообразования для монолитного массива к среднему расстоянию между трещинами.
3. Предельный размер зоны трещинообразования достигается при условии, что за время действия взрывного импульса волна напряжений в массиве успевает охватить всю зону, в противном случае размеры этой зоны ограничиваются расстоянием, пройденным продольной волной.
4. Разрушение массива по длине заряда происходит в объеме фигуры, состоящей из цилиндра с радиусом зоны трещинообразования и примыкающей к нему со стороны недозаряда (торца заряда, находящегося вблизи свободной поверхности) выпуклой половины тора с радиусом равным половине радиуса зоны трещинообразования. Граница воронки взрыва определяется линией пересечения половины тора с поверхностью обнажения, При длине недозаряда большей, чем половина радиуса зоны трещинообразования, происходит камуфлетное взрывание.
5. Установленные закономерности, подтвержденные представительными экспериментами, позволяют проектировать параметры БВР, обеспечивающие высокие техникоэкономические показатели взрывных работ, как при проходке горных выработок, так и при очистной выемке полезных ископаемых и в комплексе составляют физикотехническое обоснование параметров разрушения горного массива взрывом удлиненного заряда.
Достоверность научных положений, выводов и рекомендаций подтверждается представительным объемом промышленных экспериментов и испытаний, соответствием результатов теоретических исследований полученным данным опытно-промышленных испытаний, технико-экономическими показателями результатов внедрения.
Научная новизна работы:
- экспериментально определены размеры радиусов зон смятия и трещинообразования при взрыве удлиненных зарядов в различных горно-геологических и горнотехнических условиях;
- на основе применения решения статической задачи теории упругости при оценке размеров зон смятия и трещинообразования при взрывном воздействии на породный массив удлиненными зарядами получены зависимости, описывающие данные проведенных экспериментов с достаточной степенью точности;
- получена эмпирическая зависимость влияния естественной нарушенности массива на размер предельно возможного радиуса зоны разрушения массива;
- определено влияние времени действия взрывного импульса и скорости распространения продольной волны в массиве на развитие зоны трещинообразования;
- установлено, что образующаяся в массиве зона разрушения приближенно имеет форму цилиндра с примыкающей к нему со стороны торца заряда выпуклой половиной тора;
- получены аналитические зависимости определения минимальной длины заряда, при которой достигается полное развитие фигуры зоны разрушения в горном массиве для случаев прямого, обратного и центрального инициирования;
- разработан методологический подход определения схемы расположения зарядов и их параметров для проектирования взрывных работ при строительстве горных выработок и скважинной отбойки, комплексно учитывающий физические свойства горного массива, характеристики ВВ, размеры, форму и динамику образования зон разрушения.
ичный вклад автора состоит в: постановке задач и их решении, участии в промышленных экспериментах и испытаниях; анализе полученных результатов и выявлении зависимостей и ограничений, описывающих изучаемые явления; установлении научных и методологических основ физико-технического обоснования параметров разрушения горного массива удлиненными зарядами; разработке методических основ по определению параметров БВР при проведении подземных горных выработок и скважинной отбойки.
Практическое значение работы:
- разработана методика определения схемы расположения шпуровых зарядов и расчета их параметров для проектировании взрывных работ при проведении подземных горных выработок, позволяющая надежно обеспечивать снижение удельных расходов бурения и ВВ (от 15 % до 30 %) по сравнению с существующими их значениями при высоких значениях коэффициента использования шпуров (от 0,95 до 0,98);
- разработана методика определения схемы расположения скважинных зарядов и расчета их параметров при проектировании взрывных работ для открытых горных работ, исключающая необходимость перебура и обеспечивающая снижение удельных расходов бурения и ВВ на 15 % при снижении выхода негабарита до 1,7 % (абсолютных);
- предложен способ взрывного разрушения кристаллосодержащих горных пород удлиненными зарядами, позволяющий практически полностью исключить дробление добываемых кристаллов (Патент РФ № 2142610).
Реализация работы в промышленности. Разработанные методики и рекомендации работы внедрены и используются при проектировании скважинной отбойки на ОАО Севуралбокситруда, ОАО Боксит Тимана и при проектировании паспортов БВР для подземных горных выработок на ОАО Севуралбокситруда.
Апробация работы. Содержание работы и отдельные ее положения обсуждались на производственноЦтехнических совещаниях: концерна Совалюминий (г. Москва, 1990 г.), треста Бокситстрой (г. Североуральск, 1990 г.), НПО Сибруда (г. Новокузнецк, 1991 г.), горного отдела Комитета металлургии Российской Федерации (г. Москва, 1994 г.), ОАО Севуралбокситруда г. Североуральск 2004 г., ОАО Боксит Тимана г. Ухта 2006 г.; на Международной конференции по открытым горным, земляным и дорожным работам (г. Москва, 1995 г.); на II Международной конференции по буровзрывным работам (г. Москва, 1995 г.); на II Международной конференции по открытым горным работам (г. Москва, 1996 г.); на третьей Международной научно-технической конференции Современные технологии освоения минеральных ресурсов (г. Кемер, Турция, 2005г.); на четвертой Международной научно-технической конференции Современные технологии освоения минеральных ресурсов (г. Пекин, Китай 2006 г.).
Публикации. Основные положения диссертации изложены в 29 научных работах, в том числе 9 - в изданиях, рекомендованных ВАК РФ, в открытии, в монографии, в авторском свидетельстве СССР, в 6 патентах РФ.
Структура и объем работы. Диссертация состоит из введения, пяти глав, заключения, трех приложений, списка использованной литературы из 226 наименований, изложена на 351 страницах машинописного текста, содержит 36 таблиц и 60 рисунков.
Основное содержание работы
.
Изучением процесса разрушения горных пород взрывом занимались известные отечественные ученые: Ф.А. Баум, Е.Г. Баранов, А.Ф. Беляев, В.А. Боровиков, Д.М. Бронников, О.Е. Власов, С.С. Григорян, Г.П. Демидюк, М.Ф. Друкованый, Н.Г. Дубынин, Э.И. Ефремов, В.М. Кузнецов, Б.Н. Кутузов, Ф.И. Кучерявый, Н.В. Мельников, И.Э. Миндели, В.Н. Мосинец, Г.И. Покровский, В.В. Ржевский, В.Н. Родионов, М.А. Садовский, А.Ф. Суханов, А.Н. Ханукаев, Е.И. Шемякин, Е.Н. Шер и другие, которые заложили фундаментальные теоретические основы дальнейшего совершенствования методов расчета технологических параметров буровзрывных работ.
В первой главе приведен обзор состояния вопроса теории и практики взрывного разрушения горных пород.
Анализ существующих методик по расчету ЛНС при отбойке скальных пород на обнаженную поверхность показывает, что несмотря на их многофакторность они основаны на эмпирических зависимостях и работают в сравнительно узком спектре горно-технических условий. При значительном изменении условий проведения взрывных работ эти методики необходимо уточнять и корректировать.
Как известно, при взрыве заряда в массиве образуются три зоны: смятия (раздавливания), трещинообразования и упругих деформаций. В границах зон смятия и трещинообразования происходит разрушение массива, а их образование является общим при взрывном разрушении всех скальных массивов. Отсюда следует, что выявление закономерностей формирования этих зон может служить основой для разработки общих методологических основ определения параметров БВР в различных горно-геологических и горнотехнических условиях.
В главе 2 приводятся зависимости, позволяющие рассчитать размеры зон смятия и трещинообразования при взрыве удлиненного заряда.
На основании теории сопротивления материалов получена зависимость для определения размеров зоны смятия при взрывном нагружении скального массива удлиненными зарядами, согласно которой:
qD-0,Rсм = d = 0,3536d q0,5 D , (1) 8 где Rсм - радиус образующейся зоны смятия, м; d - диаметр заряжаемого шпура, м;
D - скорость детонации применяемого ВВ, м/с; q - плотность ВВ в заряде, кг/м3; - предел прочности пород на сжатие, Па.
Работоспособность полученной зависимости (1) проверялась в породах с различными физико-механическими свойствами. Всего были забурены 161 пара сближенных шпуров с различными расстояниями между ними. Диаметры шпуров в парах изменялись от 36 до 55 мм. В процессе экспериментов один из спаренных паралельных шпуров заряжали различными ВВ, отличающимися по скоростям детонации и плотности заряжания, а второй (компенсационный) не заряжался. Оценка размеров зоны смятия проводилась по характеру разрушения компенсационного шпура (учитывалось свойство пород в зоне смятия приобретать свойство текучести, т.е. должно происходить заполнение компенсационного шпура перемятой породой, если он находится в пределах радиуса зоны смятия). Данные, полученные в экспериментах, для установления степени достоверности зависимости (1), были обработаны методами статистической проверки гипотез.
Математическая модель связи имеет вид мультипликативной функции:
оп Rсмi = a0 dia1qia2 Dia3ia, (2) оп где Rсмi,di, qi, Di, i - соответственно радиус зоны смятия, диаметр заряда, плотность заряжания, скорость детонации ВВ и предел прочности пород на сжатие при i-м опытном взрыве; i - порядковый номер опытного взрыва; a0, Е, a4 - неизвестные параметры модели.
В результате регрессионного анализа мультипликативной функции (2) были получены оценки неизвестных параметров (a0,..,a4) и оценена степень тесноты множественной статистической связи (индекс корреляции R =0,88; расчетный критерий Фишера F(R) = 19,26, табличный критерий Фишера F(Rтабл) = 5,67).
Полученная по результатам опытных взрывов регрессия имеет вид оп 0,9777 -0,50Rсм = 0,3557d q0,4683 D1,051 (3) и хорошо согласуется с теоретической зависимостью (1).
Для определения размеров зоны трещинообразования при взрывном нагружении скального массива удлиненными зарядами на основании теории сопротивления материалов использовались две теории прочности: третья теория прочности и теория прочности О. Мора. При выводе аналитических зависимостей учитывалось влияние зоны смятия, в которой порода мелко измельчается, за счет чего происходит перераспределение напряжений в массиве. Окончательно зависимости имеют вид:
- по третьей теории прочности PВ rB -0,25 -0,RT = 1,4Rсм = 0,2102d q0,75 D1,5 ; (4) Rсм - по теории прочности О. Мора p PВ rB 1+ RT = Rсм = Rсм (5) p -0,25 -0, = 0,1486 1+ d q0,75 D1,5 , где RТ - радиус зоны трещинообразования, м; PВ - давление, развиваемое продуктами детонации ВВ, Па; rВ - радиус заряжаемого шпура; р - предел прочности пород на разрыв, Па; - предел прочности пород на срез, Па.
С целью проверки работоспособности полученных зависимостей (4), (5) в производственных условиях были забурены параллельные ряды шпуров различного диаметра (по три в ряду; 101 ряд) с переменной линией наименьшего сопротивления. Отбойка шпуров проводилась на параллельную относительно шпуров плоскость. После взрывов измерялись радиусы взрывных воронок.
Проведен сравнительный анализ между фактическими и теоретически возможными (4), (5) радиусами взрывных воронок. Из сравнительного анализа установлено, что при расчете радиуса зоны трещинообразования по зависимости (4) из всей выборки за пределы вышли 3 % данных, а по зависимости (5) вышли 30 % данных, поэтому в дальнейшем проводился анализ закономерности (4).
Для установления степени достоверности выведенной зависимости (4) по определению радиуса воронки взрыва (зоны трещинообразования) методами статистической проверки гипотез обработан статистический материал, полученный в экспериментах.
Как и при исследованиях зоны смятия, математическая модель связи представлена мультипликативной функцией, а полученная эмпирическая зависимость имеет вид:
-0,2499 -0,50RT = 0,1995 d1 q0,7503 D1,4998 . (6) В результате регрессионного анализа были получены оценки неизвестных параметров зависимости (а0,.., а5) и оценена степень тесноты множественной статистической связи (R = 0,891; F = 29,88; Fтабл = 3,99).
Поскольку F(R )= 29,88 > Fтабл,, а=0,005 = 3,99, то с вероятностью 99,5 % можно утверждать, что вариация экспериментальных радиусов воронки взрыва обусловлена на 79,3 % изменением уровней входных параметров (диаметра шпура, плотности заряжания, скорости детонации, пределов прочности пород на сжатие и срез).
Из вышеизложенного следует, что важнейшие параметры разрушения породного массива - радиус зоны смятия и предельно возможный радиус зоны трещинообразования, формирующиеся при воздействии на монолитный скальный массив взрывом удлиненного заряда, описываются закономерностями в виде произведения степенных функций от диаметра заряда, плотности заряжания, скорости детонации ВВ, пределов прочности массива на сжатие и срез.
Так как трещиноватостью обладает подавляющее большинство пород разрабатываемых месторождений, то возникает необходимость в ее изучении как фактора влияния на параметры взрывной отбойки.
Для расчета радиуса зоны разрушения трещиноватого массива при его взрывном нагружении удлиненными зарядами использовалось формула, аналогичная (4):
PB rB RPT = 1,4Rсм, (7) Rсм Кс экв где RPT - радиус зоны трещинообразования в нарушенном массиве, м; Кс - коэффициент структурного ослабления взрываемого трещиноватого массива; - эквивалентное наэкв пряжение разрушения пород взрывом, Па.
Из принятой схематизации трещиноватого массива следует, что при взрыве происходят процессы, аналогичные разрушению либо тонкостенных цилиндров, либо балок с одним или двумя защемленными концами.
Согласно теории сопротивления материалов расчет тонкостенных цилиндров ведется на разрыв, а балок по эквивалентному напряжению, равному пределу прочности материала на срез или на разрыв.
Таким образом, для расчета радиуса взрывной воронки получаем две зависимости:
PB rB RT RPT = 1,4Rсм = ; (8) Rсм Кс Kc ' PB rB RT ' RPT = 1,4Rсм =, (9) ' ' Rсм КC KC p ' ' где RT и RPT - соответственно радиусы зоны трещиообразования монолитного и нарушенного массива, определенные по разрывному разрушающему напряжению.
Объемы зоны развала, зоны дробления и их соотношения зависят от мощности применяемого ВВ, диаметра заряда и от прочностных характеристик массива. Если, например, расстояние между трещинами 0,5 м, то блочность массива 0,5х0,5 м окажет различное влияние на результаты взрывов скважинного и шпурового зарядов.
Величина взрывного импульса и прочностные характеристики массива аккумулированы в величине радиуса зоны трещинообразования, определенном для условий отбойки монолита. При этом влияние трещиноватости для условий взрывного разрушения нарушенного массива целесообразно определять в зависимости от количества трещин, приходящихся на радиус зоны разрушения, определенный для условий отбойки монолита.
С целью проверки данного предложения на шахтах Северо-уральского бокситового рудника была проведена серия экспериментов по шпуровой отбойке трещиноватого массива на обнаженную плоскость с изменением диаметров зарядов и типов взрывчатых веществ (взорвано 88 шпуров). Из формул (8), (9) можно получить выражение для Кс:
RT RT Kc = = ; (10) RPT R f ' ' RT RT ' Kc = =, (11) ' R RPT f где Rf - фактический радиус взрывной воронки, образующийся при разрушении трещиноватого массива.
' Для полученных значений KC и KC выполнен регрессионный анализ. Из регрессионного анализа установлено, что отношение (11) (в интервале от 0 до 10,2), осталось больше 1, что противоречит здравому смыслу, поэтому в дальнейшем рассматривалась схема расчета по пределу прочности на срез (8). Для выяснения вида функциональной связи коэффициента структурного ослабления Кc с трещиноватостью массива lm были рассмотрены три модели:
Kc1 = f (lm ) ; (12) Kc2 = f ; (13) lm RT Kc3 = f. (14) lm Коэффициенты корреляции между функцией и аргументом для всех моделей значимы и равны соответственно: 0,648, - 0,719, - 0,722.
Функция Кс = f(RT, lm) находилась методами регрессионного анализа по лучшим в совокупности характеристикам; наилучшей оказалась регрессия вида:
Kc = (15) RT 0,97 + 0,lm со следующими показателями: полная дисперсия 5,15; дисперсия регрессии 2,78; остаточная дисперсия 2,57; критерий Фишера 43,7 (табличный 13,1); коэффициент множественной корреляции 0,72.
Таким образом установлено: в трещиноватом массиве предельно возможный радиус зоны трещинообразования увеличивается обратно пропорционально коэффициенту структурного ослабления, являющемуся функцией отношения предельно возможного радиуса зоны трещинообразования для монолитного массива к среднему расстоянию между трещинами.
В пн. 2.2 работы исследовано влияние продолжительности действия взрывной нагрузки на размер зоны трещинообразования.
На результаты взрыва удлиненных зарядов существенное влияние оказывает длина заряда, а именно: увеличение длины заряда до какой-то определенной величины улучшает качество дробления разрушаемого массива и увеличивает объем образующейся взрывной воронки. Дальнейшее увеличение длины заряда никакого влияния на результаты взрыва не оказывают.
Исследования Э.И. Миндели показывают, что на результаты взрыва оказывает влияние не только длина, но и место инициирования колонкового заряда - прямое, центральное, обратное.
В зависимости от принятой схемы инициирования удлиненного заряда время действия взрывного импульса определяется по формулам:
lz 0,5lz lz t = ; tc = 2 ; to = 2 ;
(16) p D D D где tp, tc, to - соответственно время действия взрывного импульса при прямом, центральном и обратном инициировании, сек; lz - длина заряда, м.
Исходя из продолжительности действия взрывного импульса, волна напряжений в массиве распространяется по радиусу от точки инициирования:
lz Rtp = t cp = cp;
p (17) D lz Rtc = tc cp = c ;
p (18) D 2lz Rto = to c = c, (19) p p D где Rtp, Rtc, Rto - соответственно радиусы распространения волн напряжений в массиве за время действия взрывной нагрузки при прямом, центральном и обратном инициировании удлиненных зарядов, м; cp - скорость распространения продольной волны в массиве, м/с.
За основу расчета разрушения массива принята схема разрушения толстостенных цилиндров, у которых наружний диаметр равен бесконечности. Изменяя длину заряда, место его инициирования (а значит и длительность действия взрывного импульса) можно установить влияние продолжительности действия взрывного импульса на радиус образующейся взрывной воронки R. На практике могут возникнуть две ситуации:
f 1) времени действия взрывного импульса достаточно для охвата предельно возможной области разрушения для данных условий взрыва (рис. 1.);
2) времени действия взрывной нагрузки недостаточно для охвата предельной области разрушения для данных условий взрыва (рис. 2).
Наиболее наглядно возможные ситуации можно наблюдать при взрыве удлиненного заряда, расположенного перпендикулярно обнаженной плоскости.
Для определения зависимости радиуса взрывной воронки от длины заряда, схемы его инициирования, скорости детонации применяемого ВВ и физико-механических свойств горных пород на шахтах ОАО УСевуралбокситрудаФ был поставлен ряд экспериментов (69).
б в а а Рис.1 Схема образования воронки взрыва для случаев, когда Rti>RT при а - прямом, б - центральном, в - обратном инициировании зарядов (цифрой 1 обозначена точка инициирования заряда) а а в б Рис. 2 Схема образования воронки взрыва для случаев, когда Rti