Влияние нарушенности горных пород на их свойства и состояние Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

News of the Ural State Mining University 4 (2017)

УДК 622.278 РО! 10.21440/2307-2091-2017-4-62-65

ВЛИЯНИЕ НАРУШЕННОСТИ ГОРНЫХ ПОРОД НА ИХ СВОЙСТВА И СОСТОЯНИЕ

О. Г. Латышев, О. О. Казак

The influence of the geological material disturbance on its properties and condition

O. G. Latyshev, O. O. Kazak

The properties of the geological material and rock masses are to a degree defined by their fracture structure, or the disturbance degree. The natural massifs possess the property of hierarchy (M.A. Sadovsky), which consists in the structure's similarity on different levels of the hierarchy. In accordance with all that, this article deals with the problem of the regularity of the formation of the strength properties as well as deformative characteristics on the level of minute cracks, rock samples and block massifs. The minute cracks survey with the help of the luminescent method of crack detection is the basis of the theoretical treatment of the issues of destructibility of the rock from the energetic viewpoint. In the scale of mineral granules sizes, strength properties as well as deformative characteristics of the rocks are exponentially connected with the samples' damages on a microlevel. The quantity of the common porosity may serve as a measure of a fracture degree on a representative samples' level. The statistical dependencies of a given quantity with the basic mechanical characteristics of the rocks have been defined. The properties of the rock mass are defined by its sectional structure and geometrical characteristics of the partings in the scales, which are comparable to the mine working sizes. The sizes of rock blocks are estimated by the Weibull distribution. Considering the intrinsic fractures from the fractal positions allows to define the contact rigidity of their edges and to predict the modulus of the massif deformation. The strength passport of the rock mass is defined with the structural attenuation factor, which is adequately evaluated by the variation of the strength of the adjoining rocks and also by the modulus of fracture array. The regularities defined create the basis for projecting the technology of the processes of rock destruction. The forecast of the strength property and deformation characteristics of the rock mass are the basis of the calculation of stability of underground workings with the finite-element method.

Keywords: rocks and massifs; fracture structure; levels of hierarchy; strength; deformability.

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

Ключевые слова: горные породы и массивы; трещинная структура; уровни иерархии; прочность; деформируемость.

Теория и практика свидетельствуют о том, что свойства и состояние горных пород в первую очередь определяются степенью их нарушенности. Применительно к задачам горного производства основным показателем нарушенности является трещиноватость пород. В соответствии с теорией академика М. А. Садовского [1, 2] горные породы и их массивы представляют собой самоподобную иерархическую структуру, образованную трещинами. Свойства и состояние пород на любом уровне иерархии . определяется величиной критерия [3]: К = Я/Ь, где Я. - среднее расстояние между трещинами; Ь. - средний размер трещины. Тогда для прогноза характеристик горных пород и их массивов следует определить их свойства на различных уровнях иерархии. В данной работе исследуются три уровня: микротрещины, нарушенность на уровне образца и трещинная структура массива.

Для изучения нарушенности пород на микроуровне на кафедре шахтного строительства УГГУ разработан метод люминесцентной дефектоскопии [4, 5]. На рис. 1 представлены элек-

тронные фотографии образца (туф) на разных стадиях нагруже-ния горной породы.

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

асж = ао ехР (-й - модуль упругости:

Есж = Ео ехр (-2й Щ.

Здесь индекс «0» обозначает свойства ненарушенной горной породы; й - средний размер минеральных зерен.

По мере нагружения горных пород наблюдаются рост трещин, их ветвление и слияние с образованием кластеров (рис. 1). На этом уровне выделить отдельные трещины с оценкой их концентрации не представляется возможным. Установлено, что в целом трещинная структура горных пород обладает фрактальными свойствами [5]. Мерой нарушенности пород может служить их кластерная размерность. Для ее определения на электронное изображение поверхности образца накладывается сетка с изменяющимся размером квадратной ячейки б. Для каждого типоразмера сетки определяется количество ячеек N , занятых изображением объекта (кластером). Фрактальная кластерная размерность определится соотношением: йк = 1п №/1п (1/ б.) + 1. Анализ экспериментальных данных показывает, что кластерная размерность закономерно возрастает с увеличением напряжений а в горной породе. Эта зависимость статистически надежно описывается уравнением параболы:

й = й п + к (а2 - Ь а),

к ко

где йк0 - начальная кластерная размерность ненагруженного образца; к, Ь - эмпирические коэффициенты, зависящие от типа горной породы. Для указанного на рис. 1 образца: йк0 = 2,48; к = = 8 • 10-5; Ь = 8,75.

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

■ fn

i&rtfi 4Y'.'Г ^

I ^ rf ' v > ? У

j rC * A. i n", \ jt VMr

4 iW ' vV \ A

- Л ' ^ iV)A ^ /л ^ЖЛР1

>-. - ,

" /у /г*

Ж 1 •'>• i ■ i

v ro J kj-1 • 4k I ik,-- У 7 s >

к ' ^ , -иг

Зависимость свойств горных пород от их пористости.

Уравнения связи

Параметры уравнений связи

Эффузивные породы

Бокситы

Р = (Р, ± Рк) Р, = 3,0 • 103 кг/м3

(1 - Р) к= 0,43

h = 0,72 ± 0,05

Е = Е0 (1 - аР)2 Е = 51 ГПа; а = 3,3;

h = 0,58 ± 0,07

р, = 3,15 • 103 кг/м3 к= 0,60 h = 0,93 ± 0,03

Е = 91 ГПа; а = 2,6; h = 0,71 ± 0,12

X = Х0 (1 - а,Р) а = 18,4 МПа; а, = 7,0; а = 19,7 МПа;

г = -0,49 ± 0,08 Ь = 0,90;

Л = -0,96 ± 0,02

X = Х0 Ь1""Р асж(0) = 126 МПа; а, = 7,0;

г = -0,47 ± 0,08 а = 171 ГПа;

Р = 2,3 ГПа; а, = 7,0; Ь = 0,92; г = -0,65 ± 0,09 л = -0,93 ± 0,03

Рк№) = 2,5 ГПа; Ь = 0,93; _Л = -0,91 ± 0,03_

Здесь р - объемная масса; Е - модуль упругости; Х - прочность горных пород при растяжении ар, сжатии а^ и контактная прочность Рк (г, г| -коэффициент корреляции и корреляционное отношение).

комплексных исследований скальных пород Урала [6] установлена система статистических уравнений связи свойств пород с величиной их пористости. В таблице дана сводка структуры и параметров этих уравнений для эффузивных пород и бокситов Североуральских месторождений.

Полученная система уравнений позволяет прогнозировать свойства пород с разной степенью нарушенности (пористости). Однако само определение пористости - достаточно трудоемкая операция, связанная с отбором и доставкой проб в специально оснащенную лабораторию. Для оперативного прогноза можно использовать зависимость пористости от скорости упругой волны С, которую можно достаточно просто измерить в производственных условиях: Р(%) = 14/С, где скорость волны измеряется в км/с.

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

F (L.) = 1 - exp

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

В качестве примера рассмотрим толщу базальтов месторождения Юбилейное (Башкортостан), залегающих на глубине 250-500 м. По данным скважинных исследований1, блочное строение массива определяется системами трещин с углом падения 20-45° и 70-80° к оси керна. Средний угол 0 = 30о. По результатам вероятностной оценки блочного строения параметры распределения Вейбулла для данных условий: Ь = 0,891; т = 2,31. Тогда вероятный размер породного блока составляет Ь. = 0,14 м.

Следующей задачей является прогноз свойств трещиноватого (блочного) массива. Наиболее значимыми с точки зрения горной технологии являются деформационные и прочностные характеристики массива. В соответствии с исследованиями К. В. Руппенейта [7] модуль деформации блочного массива определится соотношением:

E = -

E

1 +

J П (l - sin4©,)

СОБ 0

X

где Ь. - размер рассматриваемого блока; а - среднее расстояние

где Е0 - модуль упругости вмещающих пород; 0. - угол наклона системы трещин, определяющий средний размер блока Ь ; Хт -линейный размер рассматриваемой области массива (как правило, определяется размером выработки); г|,- - сводная характеристика трещиноватости пород, г|- = б/£Ь., б - средняя ширина зияния трещин; ^ - относительная площадь контактов берегов трещин.

Указанные параметры определяются на стадии инженерно-геологических исследований. Неоднозначной характеристикой здесь является величина относительной площади контактов трещин Отмечая это, К. В. Руппенейт предлагает в расчетах использовать некоторую среднюю величину ^ = 4 • 10-4. Однако совершенно очевидно, что эта величина будет существенно различной для разных по геометрии трещин. Кроме того, площадь контактов будет непрерывно меняться по мере деформации массива. Теоретические и экспериментальные исследования данного вопроса [8, 9] показали следующее.

При сближении берегов трещин деформируемых пород происходит смятие неровностей поверхностей соприкосновения и, как следствие, выполаживание их траекторий. Степень неровностей (изломанности) трещин однозначно оценивается их фрактальной размерностью й. Ее определение основано на законе Ричардсона [10], который устанавливает связь размера фрактального объекта Ь(б.) с шагом его измерения б.: Ь(б.) = а 5Д где р = 1 - й. Существо закона Ричардсона состоит в том,

L

1 Гуман О. М. Инженерно-геологическая и гидрогеологическая характеристика пород Юбилейного месторождения: отчет ООО «НПЦ Урапгеолпроект»: УГГУ, 2008. 411 с.

ТЕХНИЧЕСКИЕ НАУКИ

O. G. Latyshev, O. O. Kazak / News of the Ural State Mining University 4 (2017) 62-65

1,6

H

о

0

1

о.

<u

1,5

1,4

1,3

.0

™ 1,2

© 1,1

1,0

10 15

Номер шага сближения

20

25

Рисунок 2. Изменение фрактальной размерности противоположных берегов трещины при их сближении и деформировании.

что чем меньше масштаб измерения (раствор циркуля), тем все больше подробностей измеряемого объекта (трещины) входит в результат и тем больше его размер. Это в полной мере относится к траектории трещин как бесконечно изломанной, нигде не дифференцируемой кривой. Величина фрактальной размерности определяется углом наклона графика уравнения Ричардсона в логарифмических координатах: 1п[Цб.)] = б. 1п (1 - й) + 1п а. Тогда процесс деформации трещин будет сопровождаться снижением фрактальной размерности ее берегов [9, 11]. На рис. 2 в качестве примера приведены результаты статистического моделирования деформации трещины массива Екатеринбургского метрополитена. Видно, что по мере деформации берегов разница в их фрактальной размерности уменьшается и в пределе (при полном смыкании) становится равной.

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

С

: k d„(1 - b d)>

f ov fi"

:(l- ш) 1

a

— + 1

V

л "

(l- Ш) <

+ 1

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

В исследованиях А. Н. Шашенко [13, 14] породный массив рассматривается как природная среда, нарушенная микродефектами различной природы и макродефектами в виде системы трещин, плоскостей ослабления и т. п. Для совместной оценки влияния данных факторов на прочность массива им предлагается следующая модель. Породный массив рассматривается как система плотно прилегающих друг к другу микроблоков, имеющих размеры и форму образцов, которые используются в стандартных испытаниях прочности. Те микроблоки, через которые проходит природная трещина массива, заведомо в лабораторных испытаниях не представлены и имеют нулевую прочность. Однако они содержатся в массиве и, следовательно, формируют генеральную совокупность. Тогда при организации классической выборки в статистике испытаний появятся нулевые величины. Это изменит величину дисперсии (коэффициент вариации) результатов, что и отразит влияние природной трещиноватости массива на его прочность. Вероятностный анализ данной модели позволил получить выражение коэффициента структурного ослабления массива [11]:

k =

: + b (l - V ) :(0,2/т+l) '

(1)

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

В частности, для указанной на рис. 2 трещины: ) = 1,468, к = 7,9 • 10-3, Ь = 0,68. Для любой произвольной трещины, учитывая линейность зависимости (1), для оценки параметров уравнения достаточно получить две точки (провести два измерения). Приведенный ранее анализ для блочного строения месторождения Юбилейное показал, что для толщи базальтов на глубине 250500 м модуль деформации Е = 0,9 ГПа.

Для проектирования технологии горных работ и прогноза устойчивости выработок определяющее значение имеет прочность трещиноватого породного массива. Комплексной характеристикой разрушаемости пород является паспорт прочности. Применительно к трещиноватым массивам установлено [11], что наиболее адекватной оценкой служит уравнение паспорта прочности, основанное на теории дифференциальной прочности Г. Г. Литвинского [12]:

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

Анализ закономерностей распределения прочностных свойств горных пород для условий месторождения Юбилейное позволил получить значения данных коэффициентов: а = 0,7; Ь = 0,8; с = 1,25. В частности, на рис. 3 представлены паспорта прочности базальтов в образце и их толщи в массиве.

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

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

-20

20 40 60 80 Нормальные напряжения, МПа

100

120

-Образец

■Массив

Рисунок 3. Сопоставление паспортов прочности базальта в образце и массиве.

0

5

т

0

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

ЛИТЕРАТУРА

1. Садовский М. А., Голубева Т. В., Писаренко В. Ф., Шнирман М. Г. Характерные размеры горной породы и иерархические свойства сейсмичности // Изв. АН СССР. Физика Земли. 1984. № 2. С. 3-15.

2. Садовский М. А. Избранные труды: Геофизика и физика взрыва. М.: Наука, 2004. 440 с.

3. Журков С. Н., Куксенко В. С., Петров В. А., Савельев В. Н., Султанов У. О. О прогнозировании разрушения горных пород // Изв. АН СССР. Физика Земли. 1977. № 6. С. 11-18.

4. Латышев О. Г., Корнилков М. В. Направленное изменение фрактальных характеристик, свойств и состояния пород поверхностно-активными веществами в процессах горного производства. Екатеринбург: Изд-во УГГУ, 2016. 407 с.

5. Латышев О. Г., Корнилков М. В. Исследование трещинной структуры горных пород как фрактального объекта. Saarbrücken: LAP LAMBERT Academic Publishing, 2015. 156 с.

6. Латышев О. Г. Разрушение горных пород. М.: Теплотехник, 2007. 672 с.

7. Руппенейт К. В. Деформируемость массивов трещиноватых горных пород. М.: Недра, 1975. 223 с.

8. Латышев О. Г., Прищепа Д. В. Прогноз деформационных характеристик трещиноватого породного массива // Изв. вузов. Горный журнал. 2017. № 1. С. 80-86.

9. Латышев О. Г., Прищепа Д. В., Франц В. В. Статистическое моделирование природных трещин // Изв. вузов. Горный журнал. 2016. № 5. С. 38-45.

10. Мандельброт Б. Фрактальная геометрия природы: пер. с нем. М.: Изд-во ИКИ, 2002. 656 с.

11. Латышев О. Г. Неоднородность трещинной структуры и прочность горных пород // Изв. вузов. Горный журнал. 2014. № 6. С. 152-159.

12. Литвинский Г. Г. Аналитическая теория прочности горных пород и массивов. Донецк: Норд-Пресс, 2008. 207 с.

13. Шашенко А. Н., Тулуб С. Б., Сдвижкова Е. А. Некоторые задачи статистической геомеханики. Киев: Пульсари, 2002. 304 с.

14. Шашенко А. Н., Пустовойтенко В. П. Механика горных пород. Киев: Новий друк, 2003. 400 с.

REFERENCES

1. Sadovskiy M. A., Golubeva T. V., Pisarenko V. F., Shnirman M. G. 1984,

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

gtf.shs@ursmu.ru Ольга Олеговна Казак,

Уральский государственный горный университет Россия, Екатеринбург, ул. Куйбышева, 30

Kharakternye razmery gornoy porody i ierarkhicheskye svoystva seismichnosti [Characteristic dimensions of rock and hierarchical properties of seismicity]. Iz-vestiya Akademii Nauk SSSR. Fizika Zemli [USSR Academy of Sciences Bulletin. Physics of the Earth], no. 2, pp. 3-15.

2. Sadovskiy M. A. 2004, Izbrannye trudy: geofizika i fizika vzryva [Selected works: geophysics and physics of explosion]. Moscow, 440 p.

3. Zhurkov S. N., Kuksenko V. S., Petrov V. A., Savel'yev V. N., Sultanov U. O. 1977, O prognozirovanii razrushenya gornykh porod [On the prediction of rock destruction]. Izvestiya Akademii Nauk SSSR. Fizika Zemli [USSR Academy of Sciences Bulletin. Physics of the Earth], no. 6, pp. 11-18.

4. Latyshev O. G., Kornilkov M. V. 2016, Napravlennoye izmenenie fraktal'nykh kharakteristik, svoystv i sostoyaniya porod poverkhnostno-aktivnymi veshches-vami vprotsessakh gornogo proizvodstva [Directional change of the fractal characteristics, properties and state of rocks by means of the surfactants in the process of mining]. Ekaterinburg, 407 p.

5. Latyshev O. G., Kornilkov M. V. 2015, Issledovanie treshchinnoy struktury gornykh porod kak fraktal'nogo ob'ekta [The study of the fracture structure of the rocks as a fractal object]. Saarbrücken, 156 p.

6. Latyshev O. G. 2007, Razrushenie gornykh porod [Destruction of rocks]. Moscow, 672 p.

7. Ruppeneit K. V. 1975, Deformiruemost' massivov treshchinnovatykh gornykh porod [Deformability of fractured rocks]. Moscow, 223 p.

8. Latyshev O. G., Prishchepa D. V. 2017, Prognoz deformatsionnykh kharakter-istik treshchinnovatogo porodnogo massiva [Prediction of the deformation characteristics of the fractured rock massif]. Izvestiya vuzov. Gornyizhurnal[News of the Higher Institutions. Mining Journal], no. 1, pp. 80-86.

9. Latyshev O. G., Prishchepa D. V., Franz V. V. 2016, Statisticheskoe mode-lirovanye prirodnykh treshchin [Statistical modeling of natural cracks]. Izvestiya vuzov. Gornyi zhurnal [News of the Higher Institutions. Mining Journal]. # 5. Pp. 38-45.

10. Mandelbrot B. 2002, Fraktal'naya geometriya prirody [The Fractal Geometry of Nature]: tr. From German Moscow, 656 p.

11. Latyshev O. G. 2014, Neodnorodnost' treshchinnoy struktury i prochnost' gornykh porod [The heterogeneity of fracture structure and strength of rocks]. Izvestiya vuzov. Gornyi zhurnal [News of the Higher Institutions. Mining Journal], no. 6, pp. 152-159.

12. Litvinsky G. G. 2008, Analiticheskaya teoriya prochnosti gornykh porodimas-sivov [The analytical theory of the strength of rocks and massifs]. Donetsk, 207 p.

13. Shashenko A. N., Tulub S. B., Zdzvizhkova E. A. 2002, Nekotorye zadachistatis-ticheskoy geomekhaniki [Some issues of statistical geomechanics]. Kiev, 304 p.

14. Shashenko A. N., Pustovoytenko V. P. 2003, Mekhanika gornykh porod [Mechanics of rocks]. Kiev, 400 p.

Oleg Georgievich Latyshev,

gtf.shs@ursmu.ru Ol'ga Olegovna Kazak,

Ural State Mining University Ekaterinburg, Russia