УДК 622.261.004.6

Обоснование своевременности применения способов охраны при поэтапном поддержании горных выработок.

Кандидаты техн. наук Ю.А. Петренко, А.О. Новиков, инж. А.В. Захаренко (ДонГТУ)

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

Обследование состояния горных выработок по шахтам ГХК “Донуголь” и анализ проектных решений по их сооружению и поддержанию привели к выводу, что применение способов охраны осуществляется либо на стадии сооружения, либо при невозможности дальнейшей эксплуатации выработок без дополнительных мероприятий по повышению устойчивости вмещающего массива. Такой подход к решению вопроса не носит предупредительного характера и не позволяет максимально использовать несущую способность породного массива. В этой связи перспективным направлением повышения устойчивости горных выработок является их поэтапное поддержание. Проведенные шахтные инструментальные наблюдения за деформациями вмещающего выработку массива позволили установить, что основной критерий для своевременного применения того или иного способа охраны горных выработок критические смещения их контура, которые являются результатом формирования зоны неупругих деформаций вокруг выработки во времени.

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

Рассматривалась протяженная выработка радиусом R, пройденная на глубине H и поддерживаемая крепью с несущей способностью P0, работающей в режиме постоянного сопротивления. При этом принимались допущения, что породы, вмещающие выработку, однородны и изотропны, поле напряжений нетронутого массива гидростатично и равно g Н.

Возможность замены арочной и другой формы поперечного сечения выработок на круглую для аналитических исследований обоснована в работе К.В.Руппенейта [1]. Погрешность при этом не превышает 4 – 10%.

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

Распределение напряжений на границе невесомой полуплоскости принято равнокомпонентным: Рг = Рв = g Н. Исследования, выполненные в работах [2, 3], свидетельствуют о том, что на глубинах более 600 м величина коэффициента бокового распора l стремится к единице, а в работе [4] доказано, что при 0.7 £  l  <  1 погрешность в определении напряжений не превышает 15% по сравнению со случаем, когда l =1.

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

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

, (1)

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

Величина разрушающих касательных напряжений t раз определяется из диаграммы Мора

t раз = t д–[t ], (2)

где t д действующие в породе касательные напряжения, МПа,

t д = (s q - s r)cosj /2, (3)

где s q тангенциальное напряжение, МПа, s r радиальное напряжение, МПа, j угол внутреннего трения пород, град.; [t ] – допустимые касательные напряжения, МПа,

[t ] = ( C+ s r)sinj × cosj / (1–sinj ), (4)

где С удельная сила сцепления породы, МПа.

Учитывая, что начальное распределение напряжений после проведения выработки описывается уравнениями

s q ,s r} = , (5)

где Р отпор крепи массиву горных пород, МПа; r текущий радиус, м; Rв радиус выработки, м, и подставляя (3) и (4) в (2), получаем

t раз = (6)

Подставляя полученное выражение разрушающего напряжения в дифференциальное уравнение (1) и решая его относительно t с соблюдением следующих граничных условий: при t=0 ® y = 1 , при t=tраз ® y = 0, получаем выражение для определения времени образования вокруг выработки зоны разрушенных пород tраз радиусом r

tраз = , (7)

где D безразмерный коэффициент,

D = (1–sinj )/cosj .

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

Влияние перераспределения напряжений на время разрушения любого элемента массива можно учесть с помощью принципа линейного суммирования повреждений [7]. Сущность данного принципа сводится к тому, что если элемент массива горных пород находится в течение некоторого промежутка времени D tk = tk–tk–1 под действием постоянного напряжения t k (k = 1,2,¼ ,S) и в момент t = 0, y = 1, а в момент t = ts происходит разрушение (т.е. y  = 0), то величину потери сплошности за промежуток времени D tk можно определить по формуле

y n+1k–1–y n+1k=A (n+1)t nk× D tk=D tkk, (8)

где Тk время хрупкого разрушения элемента массива при постоянном напряжении,

Тk = [A (n+1)t nk]–1,

Складывая результаты последнего соотношения для всех k, получаем

. (9)

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

, (10)

где t1, t2, tk время разрушения зон с радиусом соответственно r1, r2, rk; T1k, T2k, Tkkвремя, необходимое для разрушения k-й зоны при напряжении, действующем при рассмотрении разрушения соответственно первой, второй и k-й зоны.

Согласно формуле (10) время разрушения k-го элемента массива, расположенного в окрестностях выработки, который до момента разрушения испытывал ступенчатое изменение напряжения, будет

, (11)

где T1k , T2k , T3k ,¼ , Tkk расcчитываются по формуле (7), для чего необходимо знать параметры функции поврежденности А и n. Их значения определяли по результатам испытаний пород на длительную прочность, которые были получены И.Л.Черняком [8].

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

t раз = cosj (s разs ¥ ) / 2 . (12)

Учитывая то, что прочность массива пород в 28 раз меньше прочности в образце [9], и принимая для наших условий расчетную прочность массива в 2 раза меньше, чем в образце, формулу (12) преобразуем к виду

t раз = cosj (s разs ¥ ) / 4. (13)

Параметры функции поврежденности А и n, которые должны удовлетворять уравнению

tраз = [A(n+1) t nраз]–1, (14)

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

lgtраз=lg(n+1)A+nlgt раз. (15)

Таблица

Наименование породы

s ум,

Мпа

s ¥ ,

Мпа

А,

1/(сут× МПа)

n

Песчаник

140

100

0.0001511

1.07

Песчаный сланец

83

62

0.0003268

1.05

Глинистый сланец

43

32

0.00311

1.02

Для вычисления времени разрушения по формуле (7) необходимо также знать значение параметров С и Р.

Величина удельной силы сцепления С определяется из выражения

s ¥ /2 = (С+s ¥ /2)sinj , (16)

откуда

С = Кs ¥ (1–sinj )/(2 sinj ), (17)

где s ¥ предел длительной прочности пород при одноосном сжатии (в образце); К коэффициент ослабления, К=0,5.

Параметр Р в начальный момент времени является отпором крепи массиву горных пород. С развитием зоны разрушения вокруг выработки он представляет собой суммарный отпор крепи и разрушенного массива. По данным работы [10] его величина может быть найдена по формуле

, (18)

где Р0 отпор крепи массиву горных пород; r радиус зоны разрушенных пород; Rв радиус выработки; j угол внутреннего трения пород; Rз.р.п остаточная прочность пород в зоне разрушения, по данным работы [10] равна пределу прочности на растяжение для области ненарушенных пород.

Радиус зоны разрушенных пород определим в зависимости от граничных смещений. Для этого составим уравнение

p ( r2Rв2)Кр = p [ r2 –(RвUкр)2], (19)

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

Решая это уравнение, получим

. (20)

Средняя величина Кр может быть найдена по формуле

Кр = 1+0.667D Кр(Rв/r)2. (21)

Рекомендуется принимать D Кр = 0.2–0.25.

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

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

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

Библиографический список

1. Руппенейт К.В., Гомес Ц.О напряженном состоянии массива около подготовительной выработки. В кн.: Вопросы горного давления.– Новосибирск, СО АН СССР, 1962, вып. 12.

2. Рукин В.В., Руппенейт К.В. Механизм взаимодействия обделки напорных тоннелей с массивом горных пород. М.: Недра, 1969.

3. Либерман Ю.М. Естественное напряженное состояние массива горных пород. В сб.: Вопросы прочности подземных сооружений. Труды ВНИИСТ, 1962, вып. 12.

4. Баклашов И.В., Картозия Б.А. Механика подземных сооружений и конструкций крепи. М.: Недра, 1984.

5. Литвинский Г.Г., Бабиюк Г.В., Курман С.А. Моделирование на ЭВМ механизма разрушения пород вокруг выработки. В сб. научных трудов ДГМИ, Алчевск, 1996.

6. Литвинский Г.Г. Параметры хрупкого разрушения горных пород и условия возникновения горных ударов. В кн.: Разработка месторождений полезных ископаемых. Киев: Техніка, 1974, вып. 35.

7. Качанов Л.М. Основы механики разрушения. М.: Наука, 1974.

8. Черняк И.Л., Бурчаков А.С., Шехурдин В.И. Реологические свойства горных пород. М.: МИРГЭМ, 1966.

9. Литвинский Г.Г. Кинетика хрупкого разрушения породного массива в окрестностях горной выработки. ФТПРПИ. 1974, вып. 5.

10. Лютгендорф Г. Конвергенция горных выработок в упругом массиве при колебаниях горного давления. Глюкауф, 1975, №2.

11. Фисенко Г.Л. О взаимодействии крепи и пород в капитальных выработках. В кн.: Труды ВНИМИ.Л.; 1971, вып. 82.

Ó Ю.А. Петренко, А.О. Новиков, А.В. Захаренко, 1998