Бесплатный автореферат и диссертация по наукам о земле на тему
Исследование деформационных и реологических свойств горных пород и массивов для прогноза устойчивости подземных выработок с учетом их фрактальной геометрии
ВАК РФ 25.00.20, Геомеханика, разрушение пород взрывом, рудничная аэрогазодинамика и горная теплофизика

Автореферат диссертации по теме "Исследование деформационных и реологических свойств горных пород и массивов для прогноза устойчивости подземных выработок с учетом их фрактальной геометрии"

На правах рукописи

¿■/а.

005019723

Матвеев Алексей Александрович

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

Специальность 25.00.20 - «Геомеханика, разрушение горных пород, рудничная аэрогазодинамика и горная теплофизика»

Автореферат диссертации на соискание ученой степени кандидата технических наук

2 С Хи?

'-и й

Екатеринбург - 2012

005019723

Работа выполнена в ФГБОУ ВПО «Уральский государственный горный университет»

Научный руководитель - доктор технических наук, профессор

Латышев Олег Георгиевич.

Официальные оппоненты:

доктор технических наук, профессор Зотеев Олег Вадимович, Институт горного дела Уральского отделения РАН, заведующий лабораторией геомеханики; кандидат технических наук Аксенов Анатолий Аркадьевич, Уральский филиал научно - исследовательского института горной геомеханики и маркшейдерского дела (ВНИМИ), заведующий лабораторией.

Ведущая организация - ОАО «Научно - исследовательский и проектный институт обогащения и механической обработки полезных ископаемых» (ОАО «Уралмеханобр»).

Защита состоится 10 мая 2012 г. в 14.30 часов на заседании диссертационного совета Д 212.280.02 при ФГБОУ ВПО «Уральский государственный горный университет» по адресу: 620144, г. Екатеринбург, ул. Куйбышева, 30, 2-й учебный корпус, ауд. 2142.

С диссертацией можно ознакомиться в научной библиотеке ФГБОУ ВПО «Уральский государственный горный университет»

Автореферат разослан 05 апреля 2012 г.

Ученый секретарь диссертационного совета

доктор технических наук, профессор «-' ^ ' -"> — В. К. Багазеев

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы.

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

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

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

Работа выполнена в рамках исследований по следующей тематике: «Направленное изменение фрактальных характеристик, свойств и состояния пород поверхностно-активными веществами в процессах горного производства», в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 годы; «Теоретические основы прогноза прочности и устойчивости горных пород в подземных выработках с учетом фрактальных характеристик трещинной структуры» (задание Федерального агентства по образованию, 2011-2013 гг.).

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

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

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

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

Тема исследования соответствует паспорту специальности 25.00.20 (пп. 3,

5).

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

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

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

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

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

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

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

Защищаемые научные положения:

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

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

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

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

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

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

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

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

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

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

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

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

Реализация результатов работы. Методика прогнозирования устойчивости и разрушения горных пород в подземных выработках, включающая методы анализа лабораторных и скважинных исследований деформационных свойств массива, определения фрактальных характеристик трещин и подземных выработок, передана для использования в организации: Институт горного дела УрО РАН, ОАО «Уралгипротранс», ООО «Научно-производственное объединение VII У».

Теоретические результаты по анализу прочностных, деформационных и реологических свойств горных пород и массивов, оценке устойчивости пород в выработке используются при чтении лекций и проведении лабораторных работ по курсам: «Геомеханика», «Исследование процессов подземного строительства», «Физика горных пород».

Апробация работы

Основные результаты работы докладывались и обсуждались: на III Всероссийской молодежной научно-практической конференции по проблемам не-

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

Публикации

Основные положения диссертации опубликованы в 7 научных работах, из них 4 статьи - в ведущих рецензируемых научных изданиях.

Объем и структура работы

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

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

Глава 2 посвящена исследованию свойств горных пород Юбилейного месторождения и участков строящегося Екатеринбургского метрополитена в связи с характеристиками их трещиноватости.

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

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

Работа выполнялась в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 гг.; госбюджетной НИР «Теоретические основы прогноза прочности и устойчивости горных пород в подземных выработках с учетом фрактальных характеристик трещинной структуры», выполняемой по заданию Минобрнауки РФ (2011-2013 гг.).

Автор выражает благодарность и признательность научному руководителю, а также И. С. Осипову и В. В. Сынбулатову за дружеское участие и помощь в проведении исследований.

СОДЕРЖАНИЕ РАБОТЫ

Проблема разрушения и устойчивости горных пород в подземных выработках является ключевой при проектировании строительства и эксплуатации шахт и подземных сооружений. Базовыми параметрами, определяющими устойчивость выработок, являются свойства и состояние горных пород и породных массивов. Обширные исследования в области физики горных пород (JL И. Барон, Ж. С. Ержанов, О. Г. Латышев, Н. В. Мельников, Г. Я. Новик, М. М.

Протодьяконов, В. В. Ржевский, И. А. Турчанинов и др.) создаю^ фундамент для установления закономерностей, формирующих прочностные, деформационные и реологические характеристики пород.

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

Важнейшей характеристикой разрушаемости горных пород является паспорт прочности, выраженный огибающей предельных кругов напряжений. Используемые в настоящее время уравнения огибающей дают явно завышенные значения угла внутреннего трения ср. Так, рекомендуемое ГОСТ 21153.4-85 уравнение М. М. Протодьяконова (младшего) дает величину коэффициента внутреннего трения tg ф > 1 практически для всех скальных пород, что противоречит самой сути данного показателя. Свободной от указанного недостатка является концепция Г. Г. Литвинского, названная им «аналитической (дифференциальной) теорией прочности». Она описывает прочность структурно неоднородных и трещиноватых материалов, к которым в первую очередь относятся горные породы и их массивы. В соответствии с этой теорией уравнение огибающей принимается в виде:

т = т0

Ia

-+i , о)

'J-

где т0 - когезия сдвига, в классической интерпретации понимаемая как величина сцепления: Оо - соответствует прочности пород на растяжение. Неопределенной величиной здесь является показатель степени a (0<a< 1), который (по Г. Г. Литвинскому) оценивает долю сухого трения и зависит от площади истинного контакта минеральных зерен.

Относительную площадь истинного контакта минеральных зерен S0 (по отношению к площади разрушаемого образца S) может характеризовать показатель (О. Г. Латышев)

кс = SISo = то/тс, (2)

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

а - 0,7стс:

r=^(op+0,7ae)

2(ap+0,7осж)

- + 1

(3)

где К = 2ар - 2[ор(ар + асж)]1/2 + стсж; ар и стсж - прочность породы на одноосное растяжение и сжатие.

Непосредственными испытаниями горных пород установлена надежная зависимость (Я2 = 0,98) между коэффициентом а в формуле (1) и показателем кс:

а = 0,63+ 0,16 ln (Шс).

(4)

Таким образом, по описанной выше методике произведено построение паспортов прочности более 100 литотипов и определены статистические характеристики прочности принятых к изучению горных пород.

Анализ графиков деформации этих пород показывает следующее. В общем случае характер деформирования пород определяется тремя механизмами (стадиями). В начальной стадии нагружения наблюдается смыкание трещин и пор, что определяет отрицательную кривизну графика деформации (снижения темпа деформирования за счет уплотнения породы). После закрытия всех трещин наступает стадия чисто упругого деформирования, соответствующего закону Гука, - линейный участок графика деформации. Наконец, при достижении предела упругости наступает стадия пластической деформации, выраженная в положительной кривизне графиков деформации (опережающий рост деформации). Однако для многих пород надежное выделение этих стадий не представляется возможным в силу того, что эти механизмы действуют параллельно, и еще до достижения напряжениями предела упругости происходит развитие трещин, т. е. возникают необратимые изменения в горной породе.

Общий характер деформирования может быть описан уравнением

£ = £0/[1+Д/,о)], (5)

где Ео - начальный модуль упругости; F(/,a) - некоторая функция трещинова-тости, определяемая характеристиками трещинной структуры горных пород (/) при данном напряжений (<т). Зависимость этой функции от напряжения для широкой гаммы пород описывается квадратичным уравнением

F(/,a)==A(a2-ma), (6)

где к и m — параметры, определяемые характером трещиноватости горных пород. Для малотрещиноватых однородных пород установлено: к = 9'10"41/МПа; m = 5,8.

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

Е = Е0 [а + Ь (ст/Осж)]. (7)

где а и Ь — безразмерные эмпирические коэффициенты, зависящие от строения горных пород. Для указанных известняков а = 0,95; Ь~ 0,3 с коэффициентом корреляции г = 0,981 ± 0,012, что свидетельствует о достаточной статистической надежности данной зависимости.

Характер зависимости существенно меняется при увлажнении горных пород:

Е = Ео [А + с/(о/стсж)], (8)

где А - безразмерный коэффициент, характеризующий асимптоту графика, т. е. предельно возможное увеличение Е относительно модуля деформации; с — безразмерный параметр, характеризующий скорость роста Е с нагрузкой. Для изученных пород А = 1,72 и а = 0,12.

Для принятых к изучению пород определены их реологические характеристики: ползучесть, релаксация напряжений, длительная прочность. Для гор-

ных пород наиболее широко применяется теория наследственной ползучести (Ж. С. Ержанов). Центральная идея данной теории состоит в предположении, что ползучесть материала характеризуется определенной зависимостью деформации в данный момент времени от всей предыдущей истории его деформирования. Обобщенное уравнение реологии определяется видом ядра ползучести. На первом этапе производилось сравнение различных выражений ядра ползучести. Установлено, что наилучшее приближение к опытным данным обеспечивает степенное ядро ползучести (функция Абеля) в виде: ЦТ,т) = 8 (7 - т)"*1, где 5 и а — определяемые в эксперименте параметры ядра ползучести.

Результаты исследований горных пород в воздушно-сухом и увлажненном состоянии показали следующее. В зависимости от уровня приложенных напряжений горные породы проявляют затухающую или незатухающую ползучесть. Причем с увлажнением пород переход к незатухающей ползучести происходит при меньших по величине напряжениях. Количественно рост деформации во времени достаточно точно описывается уравнением Больцмана-Вольтерра со степенным ядром ползучести: е(г) = е0 + е0 [5/(1- а)] г1'", где Ео -мгновенная деформация (упругий прогиб балки). Однако финишные участки кривых затухающей ползучести данное уравнение заведомо не описывает. Эта стадия может быть описана гиперболической зависимостью

е = (/(ш + Ь) (9)

с пределом = На. Для изученных пород а = (5 - 6)103.

По результатам экспериментального изучения параметров ползучести определены функции релаксации напряжений: о(Г) = Ег0 + ^{¿е/Ж) - для затухающей и = Еео ехр(-£/г|7) - для незатухающей ползучести, где п - коэффициент динамической вязкости (ц = Е Г0), Па-с; /о - период релаксации.

Анализ результатов исследования ползучести показывает, что в ходе деформирования пород происходит перестройка их структуры, проявляющаяся в развитии трещин и подвижках минеральных зерен. В этой связи такие параметры как модуль упругости Е и коэффициент Пуассона V уже не являются константами породы, что затрудняет расчеты напряженно-деформированного состояния горных пород и массивов. Ю. М. Работнов показал, что при установлении и использовании временных функций данных показателей: Е({) = Ео/[ 1 + и у(/) = 0,5 - (0,5-у)/[1+/г(0] задачу линейной наследственности можно рассматривать как задачу теории упругости. Тогда задача сводится к определению функции Дг). В результате исследований нами определены параметры данной функции для всех принятых к изучению горных пород:

= ' (10) 1 - а

Процессы упругого деформирования и ползучести сопровождаются и определяются изменением трещинной структуры горных пород. Для изучения данного вопроса выполнена серия экспериментальных исследований. Образцы горных пород ступенчато нагружались на прессе с параллельным замером их деформации. После каждой ступени нагружения методом люминесцентной дефектоскопии получали цифровые фотографии трещинной структуры пород и по разработанной методике (О. Г. Латышев, И. С. Осипов, А. Н. Еремизин) опре-

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

Для изучения данного процесса использован аппарат фрактального кластерного анализа (Б. Мандельброт). Количественной мерой кластеров служит фрактальная размерность поверхности образца. Для ее определения использован метод покрытия. На изображение поверхности (фотоснимок) накладывается сетка с изменяющимся размером квадратной ячейки 5,. Для каждого типоразмера сетки определяется количество ячеек ЛГ/, занятых изображением объекта (кластером). Фрактальная кластерная размерность трещинной структуры определяется по углу наклона графика данной зависимости, построенного в двойных логарифмических координатах:

1пЛГ .

(И)

- + 1.

1п(1/5)

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

2,9 -

2,8

2,7

2,6

О 20 40 60 60

Напряжение, МПа

Рис. 1. Рост фрактальной размерности кластеров с увеличением нагрузки

Данную зависимость можно описать уравнением параболы вида:

с!^ + К (а2 - тп о), (12)

где й^о - фрактальная кластерная размерность поверхности исходной горной породы; та - эмпирические коэффициенты, отражающие темп роста фрактальной размерности по мере нагружения породы (ка) и кривизну графика зависимости (та). Для изученных образцов туфа ¿¡^ = 2,48; ка = 8'10'51/МПа и та = 8,75 (при измерении напряжений а в МПа). Коэффициент достоверности аппроксимации (коэффициент детерминации) Я2 = 0,995 свидетельствует о статистической надежности данной зависимости. Положительная кривизна графика отражает ускоряющийся рост дефектности пород с нагрузкой, что соответствует современным представлениям о накоплении повреждений в ходе разрушения горных пород.

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

Фрактальная кластерная размерность

Рис. 2. Зависимость модуля деформации от фрактальной размерности кластеров

горных пород

Данная зависимость описывается уравнением

£ = *(«/,- 2). (13)

Для указанной породы (туфа) к = 5,1 ГПа. Близкий к единице (г = 0,968) коэффициент корреляции свидетельствует о статистической значимости данной закономерности. Следует отметить важное обстоятельство. Полученные соотношения справедливы до момента начала разрушения горной породы, когда лавинообразный рост трещин сопровождается прогрессирующим пластическим течением. Здесь действуют другие закономерности.

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

Полученные результаты являются базой для прогноза устойчивости и разрушаемости горных пород в подземной выработке. Однако для обеспечения достоверности прогноза требуется оценка свойств и состояния породного массива. Из современных методов натурных испытаний породного массива наиболее разработанными в техническом и методическом плане являются скважин-ные исследования. Такие исследования проведены в массиве строящегося Екатеринбургского метрополитена с помощью прессиометра (дилатометр И7 096 швейцарской фирмы БокхреП АС). Произведены измерения 25 интервалов длиною по 1 м в пяти скважинах. Для построения графиков деформации и их анализа разработана соответствующая компьютерная процедура.

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

О 5 Ю 15

Деформация (Е-3)

Рис. 3. Выделяемые стадии деформирования породного массива

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

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

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

Е0 = ов/Ео- (14)

Следующим шагом является выделение необратимых (пластических) деформаций массива. Для этого на каждом шаге нагружения из общей деформации ео6щ вычитается ее упругая составляющая, определяемая модулем упругости Ем:

£пл(|)= Б0бщ(,) — а,/£м. ■ (15)

По вычисленным значениям ёШ1(0 производится построение графика необратимых деформаций массива и оценивается величина модуля пластичности Ет как тангенса угла наклона линии тренда данного графика. Аналогично методике исследования деформации образцов для породного массива определяется уравнение функции трещиноватости (6). Для различных участков массива величина т = 0,2, а коэффициент темпа закрытия трещин изменяется в пределах к = 0,5 - 1,5. По мере смыкания трещин (изменения функции трещиноватости) закономерно увеличивается модуль деформации массива. Его относительное

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

Обобщение и анализ исследований позволили установить надежные зависимости модулей общей деформации породного массива от модуля упругости слагающих его пород £0:

£даф = 0,17£0 (16)

с коэффициентом корреляции г = 0,91.

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

Паспорт прочности горных пород в массиве в соответствии с концепцией Г. Г. Литвинского определится на основе зависимости (1): т = т0ч*[а/уоо + l]a\ где показатель сплошности у может быть оценен соотношением:

у = £0(ГП)/£о(М), (17)

где £0(т) - модуль упругости горных пород, измеренный на образце; £0(м) - модуль упругости слагающих массив горных пород. Первый показатель определяется лабораторными испытаниями в начальной стадии нагружения образцов, методика которых изложена выше. Модуль упругости слагающих массив пород (Ео(м)) определяется по результатам скважинных испытаний (формула (14)).

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

На основе изучения законов распределения прочности горных пород и статистики его экстремальных (минимальных) значений установлено уравнение масштабного эффекта: as = om¡n + P/S03/2*. Здесь S0 - площадь обнажения пород (площадь горной выработки); amin = ц + p(-lnP)1" - показатель, характеризующий минимальное статистически вероятное значение частного результата прочности при S -» оо, соответствующее прочности неограниченного массива; ц - крайняя левая граница прочности, т. е. асимптота кривой масштабного эффекта; Р - параметр масштабного фактора, имеющий размерность напряжения; к -коэффициент однородности горной породы; Р~ принятый уровень значимости.

Паспорт прочности трещиноватого (блочного) массива описывается уравнением:

-jqAW

т = *о Л*™)

- + 1

Я* стрК кту

дулем трещиноватости породного массива JT как f(k^) = 1,08 JT

(18)

Здесь функция коэффициента структурного ослабления оценивается мо-

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

а + Ь (1 - Уа) *ств =-—. (19)

стр с (0,2УТ +1)

где Уа - коэффициент вариации прочности пород; а, Ь, с - эмпирические коэффициенты, Для условий Юбилейного месторождения а = 0,7, Ъ = 0,8, с = -1,25. По результатам исследований составлены прогнозные диаграммы прочности породного массива Юбилейного месторождения по глубине.

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

Все методы прогноза устойчивости и нагрузки на крепь подземных выработок разделяют на две группы: расчет по схеме заданной нагрузки и расчет по схеме заданной деформации. В первом случае при критических напряжениях предполагается разрушение пород на контуре выработки с образованием свода естественного равновесия (М. М. Протодьяконов, П. М. Цимбаревич, В. Д. Сле-сарев, Н. М. Покровский и др.). Второй подход (расчет по заданным деформациям) основан на оценке образования вокруг выработки зоны неупругих (запредельных) деформаций и использовании методов механики сплошной среды (А. Лабасс, К. В. Руппенейт, Ю. М. Либерман, А. Н. Ставрогин, А. Г. Протосеня и др.). При переходе к содержательному расчету неизбежно возникает задача выбора той или иной концептуальной схемы. Из постановки задачи и самой сути различных гипотез горного давления следует, что концепцию образования свода естественного равновесия следует применять при «малых» глубинах заложения выработки, а концепцию заданных деформаций - при «больших» глубинах. Однако конкретных (количественных) признаков «малых» и «больших» глубин в литературе не приводится. Это совершенно понятно, поскольку никакие жесткие рекомендации не могут учесть всего многообразия особенностей реальных массивов. Но для конкретных горно-геологических условий такие рекомендации могут быть получены.

В принципе эта задача так или иначе решается в любых исследованиях устойчивости. Анализ позволяет выделить следующие основные подходы. В соответствии с одним из них (А. П. Максимов) разрушение массива определяется чистым сдвигом, и выбор критической глубины основывается на критерии Кулона-Мора. Другой подход основан на критерии нормальных напряжений с учетом комплексной деформационной характеристики массива (Н. С. Булычев). В исследованиях И. В. Баклашова и О. В. Тимофеева учитывается технология проведения подземной выработки, в частности влияние буровзрывных работ. Общий подход к решению данной задачи предложен Г. Г. Литвинским. В качестве критерия им вводится безразмерный параметр - «локальный нормированный критерий разрушения (ЛНКР)»:

= (20)

где - функция действующего тензора напряжений в данной точке масси-

ва: ^(о,С,у) - функция прочности пород с учетом действующих напряжений. Однако конкретного выражения критерия Г. Г. Литвинским не приводится. Это вполне объяснимо, так как учесть все особенности разнообразных реальных горно-геологических условий в единой модели не представляется возможным. Поэтому для прогноза устойчивости горных пород в выработке необходима разработка конкретной формулы критерия, отражающей реальные условия проходки выработки. На основе функции (20) и методических подходов Н. С. Булычева и И. В. Баклашова нами предлагается следующая формула критерия:

У = . (21) КнКсК3Квасж

Конкретное наполнение данного расчетного соотношения производится по результатам проведенных в диссертации исследований. Здесь Кх - коэффициент, учитывающий отличие реального напряжения массива от гидростатического (у//), может быть достоверно определен только натурными измерениями НДС в конкретной точке массива. При отсутствии такой информации следует использовать рекомендации СНиП. Величина Ки = К1УКР учитывает возможное разупрочнение пород при увлажнении (К¡у) и влияние горного давления (Кр). Коэффициент Кс учитывает влияние масштабного фактора и трещинной структуры массива. Приведенные выше результаты исследований позволяют определить величину данных коэффициентов. В дополнительной оценке нуждаются показатели, отражающие деформационные характеристики массива особенности технологии проходки выработки (Кд) и концентрацию напряжений на ее контуре (Ка).

Коэффициент Н. С. Булычевым К$ предлагается определять по формуле: К3 = 1 + 1/эт ф (П/п ф - 1), где П£ = есж/е>, - характеристика пластичности горных пород, представляющая собой отношение предельной деформации при разрушении пород сжатием (есж) и ее упругой составляющей (еу). Величина П£ определится по деформационным характеристикам горных пород (уравнения (5) - (8)) или результатами скважинных исследований (см. рис. 4). Угол внутреннего трения ф при данном напряжении а = уН определится анализом паспорта прочности по изложенной выше методике.

В работах И. В. Баклашова и О. В. Тимофеева установлено, что изменение модуля деформации массива Е по мере удаления от контура выработки г за счет разгрузки массива и нарушения пород взрывными работами (образования зоны трещиноватости) можно учесть коэффициентом Кд = (1 - а г '"), где л= 1/^ (1 + О; а = 0,98". Параметр £,, зависящий от глубины зоны нарушения породного массива взрывными работами, определяется шахтными наблюдениями. На стадии проектирования величину рекомендуется определять по таблицам в зависимости от коэффициента крепости пород и глубины заложения выработки.

Показатель К„ можно сопоставить с коэффициентом концентрации напряжений на контуре выработки, величина которого определяется размерами

обнажения и характером неровностей ее контура. Наиболее распространенная методика оценки величины Ка заключается в следующем (И. В. Баклашов, К. В. Руппенейт). Контур выработки описывается системой параметрических уравнений, решение которой находится в комплексном виде (метод Колосова-Мусхелишвили). Для упрощения решения задачи строится конформное отображение выработки, т. е. такой линии контура, для которого математическое ожидание отклонений равно нулю (£С0 = 0). С учетом вероятностной природы отклонений контура выработки от проектного коэффициент концентрации напряжений вычисляется по формуле

Ка = о»/уН = 2[1 + q(P) 2я m Dm\, (22)

где m — среднее число пересечений линии реального контура выработки с ее проектным профилем на участках г, равных по длине эквивалентному радиусу выработки Rb', D - дисперсия величин отклонения контура выработки от проектного сечения Се, выраженная в долях Rb- Величина множителя q(P) зависит от принимаемого уровня надежности расчетов Р: при Р = 85 % ^(0,85) = 1,44; при Р = 95 % 9(0,95) = 1,96; при Р = 99 % д(0,99) = 2,66.

В соответствии с формулой (22), величина коэффициента концентрации напряжений определяется двумя характеристиками: дисперсией отклонений D и общим числом пересечения профилем выработки ее проектного (конформного) контура т. Эти показатели косвенно характеризуют степень изломанности линии контура выработки, но если дисперсия D является объективным показателем геометрии выработки, то величина m представляется малоинформативной. Более объективной характеристикой степени изломанности линии контура выработки является ее фрактальная размерность df. Для изучения данного вопроса разработана математическая имитационная модель (метод Монте-Карло) генерации контура горной выработки, основанная на теории итерационных функций и ее случайных реализациях. Корреляционный анализ результатов реализации модели показал, что фрактальная размерность практически точно соответствует данным характеристикам выработок. Полученное уравнение связи

df= 0,42 In {mDm) (23)

характеризуется коэффициентом достоверности аппроксимации, близким к единице (R2 = 0,985). С учетом этого уравнение оценки коэффициента концентрации напряжений можно представить в следующем виде:

/¡Q^l + IßgtfOÄJ*-1. (24)

Для проверки применимости (адекватности) данного уравнения исследованы 65 сечений горизонтальных выработок СУБРа. По результатам шахтных замеров произведено вычисление коэффициентов концентрации напряжений по классической формуле (22) и по предлагаемому уравнению (24), основанному на определении фрактальной размерности линии контура выработки. Результаты сравнения и анализа показали следующее (рис. 5). Фрактальная размерность ¿//достаточно точно соответствует коэффициенту концентрации напряжений Ка на контуре горных выработок. На рис. 5 точками отмечены значения Ка, вычисленные традиционным способом по формуле (22); линия соответствует уравнению (24) при q{0,95) = 1,96 и ß = 0,8. Коэффициент вариации составил 1,6 %.

3,5

3,4

X

X ■з 3,3

1 г 3.2

я X X 3,1

и ь 1 3

© X с а 2,9

л X

X •е- 2,8

•8» А 2.7

Ж 2,6

2,5

• Г—

< * / п

• > г

• -• • ! * ¡и— •

1.3 1.4 1,5 1.6 1,7 1,8 1,9 Фрактальная размерность

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

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

Определенные по результатам исследований параметры уравнения (21) позволили составить прогнозную диаграмму устойчивости горных пород в выработках по глубине Юбилейного месторождения (рис. 6). Сопоставление ее с результатами расчетов по критериям Н. С. Булычева и И. В. Баклашова показало их качественное соответствие.

По полученной диаграмме определены интервалы глубин, где породы устойчивы и выработки можно не крепить; крепление необходимо, и расчет следует производить по схеме заданной нагрузки; расчет крепи требуется производить с учетом реологии массива (расчет по схеме заданной деформации). По схеме с заданной нагрузкой на основании рекомендаций СНиП П-94-80 определены значения критерия устойчивости С и нагрузки на крепь ствола по условиям Юбилейного месторождения. При этом приведенные выше результаты исследований позволили конкретизировать многие коэффициенты СНиП, что повысило надежность прогноза.

Аналогичные расчеты выполнены с учетом изученных деформационных и реологических свойств горных пород. Определено давление на крепь выработок по методике К. В. Руппенейта. В соответствии с этой методикой упругие свойства трещиноватого породного массива зависят от геометрической характеристики систем трещин:

Л, = 5,//,- Ь, (25)

Рис. 6. Прогнозная диаграмма устойчивости горных пород в выработках по глубине Юбилейного месторождения

где 5, - средняя ширина раскрытия трещин /-Й системы; /, - среднее расстояние между трещинами; ^ - относительная площадь скальных контактов берегов трещин. Расстояние между трещинами отдельностей /,• и ширину их зияния 5, можно получить усреднением натурных замеров. Более сложной в определении является относительная площадь скальных контактов берегов трещин К. В. Руппенейтом рекомендуется принимать ^ = 3-10'4 как • некоторую константу. Однако совершенно очевидно, что эта величина будет существенно различной для реальных трещин массива. Более того, она будет меняться в ходе деформирования массива. Следовательно, данный вопрос нуждается в специальных исследованиях.

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

мерности средней линии трещины б//. Реализация модели позволила установить зависимость относительной площади контактов берегов трещины от ее фрактальной размерности (рис. 7).

6.00е-04

а

о

& 5.00е-04

о

л 4.00е-04

9

3

о | 3,00е-о4

к

X 2.00е-04

с

г

5 и 1.00е-04

X ь

о 0.00е+00

1,42 1,43 1,44 1,45 Фрактальная размерность

1,46

Рис. 7. Зависимость площади контактов берегов трещины от ее фрактальной размерности

В общем виде уравнение данной зависимости запишется в виде:

£• = £¿/0 (1-& (26) где фрактальная размерность средней линии исходной трещины; к и Ь -эмпирические коэффициенты, зависящие от геометрии трещины и соотношения нормальных и сдвигающих деформаций. Изучение реальных трещин массива показало, что по порядку величин: к = (5-15)-10"3, Ъ = 0,5-0,8.

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

' КаЧП

и = Аг0

ех]

м.

-1

(27)

где г0 - эквивалентный радиус выработки; кс - коэффициент структурного ослабления пород. При А = 26 данное уравнение точно соответствует номограмме для определения типовых смещений СНиП Н-94-80 (рис. 8), но имеет то преимущество, что здесь в явном виде учитываются размеры выработки, коэффициенты концентрации напряжений и структурного ослабления.

Результаты исследования реологических характеристик горных пород позволяют определить функцию ползучести Ф = 8 / "7(1 - а), где а и 5 - параметры ядра ползучести, и длительную прочность пород адл, которые являются основой прогнозирования устойчивости выработок по схеме заданной деформации (Ж. С. Ержанов). Методика расчета базируется на известном решении А. Н. Динника о НДС идеально упругого изотропного массива с заменой (в соответ-

ствии с идеей Ю. Н. Работного) упругих констант горных пород их временными операторами.

-(5-5

-К=20

И=30

15=50 И=60 Я=70

0 200 400 600 800 1000 1200 1400 1600 1800 Глубина, м

Рис. 8. Номограмма к определению типовых смещений

В частности, временная функция модуля упругости: £(/) = Е0-ехр [-ш 5 Г(1 - а) где со = (1 - а)1"0; Г() - гамма-функция. По аналогичным формулам определены зависящие от реологических свойств пород коэффициент Пуассона у(/), коэффициент горизонтального распора и произведен расчет напряженно-деформированного состояния пород вокруг выработки применительно к изученным горным породам: радиальные напряжения - стг = Х(/) (1 - Я2/г2) уН, тангенциальные напряжения - о9 = Я.(г) (1 + Я2/г2) уН и перемещения - Ц. = { Х(/) [1 + у(0]/[1 - Е(0} уН, где Л - радиус выработки, г — расстояние от ее оси.

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

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

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

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

7. По схеме с заданной нагрузкой на основании рекомендаций СНиП II-94-80 определены значения критерия устойчивости С и нагрузки на крепь ствола по условиям Юбилейного месторождения. При этом приведенные выше результаты исследований позволили конкретизировать многие коэффициенты СНиП, что повысило надежность прогноза.

8. На основе разработанной математической модели деформации трещины определены закономерности формирования упругих и реологических характеристик породного массива; выполнены расчеты устойчивости подземных выработок по схеме заданной деформации.

9. Методика прогнозирования устойчивости и разрушения горных пород в подземных выработках, включающая методы анализа лабораторных и скважинных исследований деформационных свойств массива, определения фрактальных характеристик трещин и подземных выработок, передана для использования в организации: Институт горного дела УрО РАН, ОАО «Уралгипрот-ранс», ООО «Научно-производственное объединение УТТУ». Теоретические результаты по анализу прочностных, деформационных и реологических свойств горных пород и массивов, оценке устойчивости пород в выработке используются при чтении лекций и проведении лабораторных работ по курсам: «Геомеханика», «Исследование процессов подземного строительства», «Физика .горных пород».

Опубликованные работы по теме диссертации.

Статьи, опубликованные в ведущих рецензируемых научных журнал лах, определенных Высшей аттестационной комиссией:

1. Латышев О. Г., Соколов В. В., Матвеев А. А. Оценка коэффициента концентрации напряжений по фрактальной размерности контура горной выработки // Изв. вузов. Горный журнал. - 2009. - №3. - С. 47-52.

2. Латышев О.Г., Сынбулатов В.В., Матвеев A.A. Анализ деформационных характеристик трещиноватого массива при строительстве Екатеринбургского метрополитена // Горный информационно-аналитический бюллетень. Вып. 6. - М.: Горная книга, 2010. - С. 97-103.

3. Прогноз деформационных характеристик трещиноватых горных пород и массивов / Латышев О. Г., Матвеев А. А., Мартюшов К. С., Еремизин А. Н. // Изв. вузов. Горный журнал. - 2011. -№ 7. - С. 92-97.

4. Матвеев А. А. Исследование параметров трещинной структуры горных пород и породных массивов // Изв. вузов. Горный журнал. - 2012. - № 2. - С.

5. Латышев О.Г., Соколов В.В., Матвеев A.A. Анализ паспорта прочности скальных горных пород // Проблемы недропользования: материалы III Всероссийской молодежной научно-практической конференции, 10-13 февраля 2009 г., г. Екатеринбург. - Екатеринбург: ИГД УрО РАН, 2009. - С. 287-296.

6. Латышев О.Г., Соколов В.В., Матвеев A.A. Оценка структурного ослабления горных пород в массиве // Добыча, обработка и применение природного камня: сб. науч. тр. - Магнитогорск: МГТУ, 2010. - С. 279-285.

7. Оценка деформационных характеристик породного массива для прогноза устойчивости горных выработок / Латышев О.Г., Матвеев A.A., Сынбулатов В.В., Строев Ю.М. II Труды III Международной конференции «Проектирование, строительство и эксплуатация комплексов подземных сооружений», 19-21 мая 2010 г., г. Екатеринбург. - Екатеринбург: УГГУ, 2010. -С. 223-228.

Подписано в печать Формат 60x84 1/16

Бумага писчая. Печать на ризографе

Печ. л. 1,0. Тираж 120 экз. Заказ 50.

Издательство Уральского государственного горного университета 620144, г. Екатеринбург, ул. Куйбышева, 30 Отпечатано с оригинал-макета в лаборатории множительной техники УГГУ

64-67.

Статьи, опубликованные в других изданиях:

Текст научной работыДиссертация по наукам о земле, кандидата технических наук, Матвеев, Алексей Александрович, Екатеринбург

61 12-5/2751

ФГБОУ ВПО «Уральский государственный горный университет»

На правах рукописи

(7

Матвеев Алексей Александрович

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

Специальность 25.00.20 - «Геомеханика, разрушение горных пород, рудничная

аэрогазодинамика и горная теплофизика»

Диссертация на соискание ученой степени кандидата технических наук Научный руководитель - д-р техн. наук, профессор Латышев О.Г.

Екатеринбург - 2012

ОГЛАВЛЕНИЕ

ВВЕДЕНИЕ.......................................................................................................................4

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

1.1. Требования нормативных документов..............................................................8

1.2. Теория и методы оценки устойчивости и разрушаемости горных пород в выработках......................................................................................................................18

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

1.4. Гипотезы горного давления и формирования нагрузки на крепь.................33

1.5. Элементы фрактальной геометрии...................................................................39

1.5.1. Введение фрактальной размерности.....................................................39

1.5.2 Виды фрактальных размерностей..........................................................41

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

2. ИССЛЕДОВАНИЕ СВОЙСТВ ГОРНЫХ ПОРОД.................................................46

2.1. Объекты и методы исследований.....................................................................46

2.1.1. Объекты и предмет исследований.........................................................46

2.1.2. Методики определения плотностных и водно-физических свойств. 49

2.1.3. Методика построения паспортов прочности горных пород...............50

2.1.4. Методика исследования деформационных характеристик................51

2.1.5. Методы реологических испытаний.......................................................53

2.2. Плотностные и прочностные свойства горных пород...................................55

2.2.1. Объемная масса горных пород..............................................................55

2.2.2. Прочностные и упругие свойства горных пород.................................57

2.2.3. Анализ паспорта прочности горных пород..........................................60

2.3. Деформационные характеристики горных пород...........................................64

2.3.1. Общий характер деформирования........................................................64

2.3.2. Определение модуля деформации.........................................................67

2.3.3. Влияние влажности на деформацию горных пород............................70

2.4. Реологические свойства горных пород............................................................72

2.4.1. Выбор уравнения ядра ползучести........................................................72

2.4.2. Результаты экспериментов.....................................................................77

2.4.3. Анализ результатов исследования реологических характеристик.... 80

2.5. Связь деформационных характеристик с трещиноватостью горных пород ..........................................................................................................................................83

2.5.1. Постановка и результаты экспериментальных исследований..........83

2.5.2. Анализ результатов исследования........................................................90

Выводы.......................................................................................................................95

3. ПРОГНОЗ СВОЙСТВ ПОРОДНОГО МАССИВА.................................................97

3.1. Общие положения..............................................................................................97

3.2. Скважинные исследования деформационных характеристик массива........98

3.2.1. Приборное обеспечение экспериментов...............................................98

3.2.2. Планирование и организация экспериментальных исследований.....99

3.2.3. Анализ результатов исследований......................................................102

3.3. Прогноз деформационных характеристик породного массива...................105

3.3.1. Стадии деформирования породного массива....................................105

3.3.2. Методика определения деформационных характеристик массива. 106

3.3.3. Закономерности деформирования породных массивов....................109

3.4. Оценка деформационных характеристик массива после его укрепления цементацией..................................................................................................................114

3.5. Паспорт прочности горных пород в массиве................................................116

3.6. Оценка свойств горных пород в массиве.......................................................118

3.6.1. Учет влияния горного давления и влажности на свойства пород... 118

3.6.2. Оценка масштабного фактора..............................................................120

3.6.3. Определение коэффициента структурного ослабления....................123

3.6.4. Прогноз прочности породного массива..............................................128

Выводы.....................................................................................................................129

4. ПРОГНОЗ УСТОЙЧИВОСТИ ПОДЗЕМНЫХ ВЫРАБОТОК............................131

4.1. Прогноз устойчивости по схеме заданной нагрузки для условий Юбилейного месторождения.......................................................................................131

4.1.1. Определение критерия и категории устойчивости горных пород... 131

4.1.2. Определение нагрузок на крепь вертикального ствола....................136

4.2. Использование фрактальной размерности линии контура подземной выработки для прогноза концентрации напряжений...............................................140

4.2.1. Методы оценки коэффициента концентрации напряжений на контуре выработки......................................................................................................................140

4.2.2. Фрактальная размерность как мера неровностей контура подземной выработки......................................................................................................................142

4.2.3. Расчет коэффициента концентрации напряжений по фрактальной размерности линии контура выработки.....................................................................145

4.3. Оценка деформационных характеристик трещиноватого породного массива ........................................................................................................................................146

4.4. Устойчивость горных пород в выработке с учетом деформации и реологии массива...........................................................................................................................153

4.4.1. Общие принципы расчета по схеме заданной деформации.............153

4.4.2. Прогноз устойчивости с учетом деформации породного массива.. 155

4.4.3. Прогноз устойчивости с учетом реологии породного массива........159

4.5. Критерии устойчивости горных пород в подземных выработках..............163

4.6. Реализация результатов исследований..........................................................170

Выводы.....................................................................................................................170

ЗАКЛЮЧЕНИЕ............................................................................................................172

БИБЛИОГРАФИЧЕСКИЙ СПИСОК.........................................................................174

ПРИЛОЖЕНИЕ 1 180

ПРИЛОЖЕНИЕ 2 199

ПРИЛОЖЕНИЕ 3 204

ВВЕДЕНИЕ

Актуальность темы.

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

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

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

Работа выполнена в рамках исследований по следующей тематике: «Направленное изменение фрактальных характеристик, свойств и состояния пород поверхностно-активными веществами в процессах горного производства», в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 годы; «Теоретические основы прогноза прочности и устойчивости горных пород в подземных выработках с учетом фрактальных характеристик трещинной структуры» (задание Федерального агентства по образованию, 2011-2013 гг.).

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

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

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

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

Тема исследования соответствует паспорту специальности 25.00.20 (пп. 3,

5).

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

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

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

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

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

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

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

Защищаемые научные положения:

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

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

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

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

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

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

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

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

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

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

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

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

Реализация результатов работы. Методика прогнозирования устойчивости и разрушения горных пород в подземных выработках, включающая методы анализа лабораторных и скважинных исследований деформационных свойств массива, определения фрактальных характеристик трещин и подземных выработок, передана для использования в организации: Институт горного дела УрО РАН, ОАО «Уралгипротранс», ООО «Научно-производственное объединение УГГУ».

Теоретические результаты по анализу прочностных, деформационных и реологических свойств горных пород и массивов, оценке устойчивости пород в выработке используются при чтении лекций и проведении лабораторных работ по курсам: «Геомеханика», «Исследование процессов подземного строительства», «Физика горных пород».

Апробация работы

Основные результаты работы докладывались и обсуждались: на III Всероссийской молодежной научно-практической конференции по проблемам не-

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

Публикации

Основные положения диссертации опубликованы в 7 научных работах, из них 4 статьи - в ведущих рецензируемых научных изданиях.

Объем и структура работы

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

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

Глава 2 посвящена исследованию свойств горных пород Юбилейного месторождения и участков строящегося Екатеринбургского метрополитена в связи с характеристиками их трещиноватости.

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

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

Работа выполнялась в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 гг.; госбюджетной НИР «Теоретические основы прогноза прочности и устойчивости горных пород в подземных выработках с учетом фрактальных характеристик трещинной структуры», выполняемой по заданию Минобрнауки �