Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Сдвижение и дегазация пород и угольных пластов при очистных работах

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
8.07 Mб
Скачать

АКАДЕМИЯ НАУК УКРАИНСКОЙ ССР

ИНСТИТУТ г е о т е х н и ч е с к о й м е х а н и к и

В. И. Мякенький

СДВИЖЕНИЕ И ДЕГАЗАЦИЯ ПОРОД

УГОЛЬНЫХ ПЛАСТОВ ПРИ ОЧИСТНЫХ РАБОТАХ

КИЕВ « ПАУКОВА ДУМКА » 1У7&

МЭ9

6 П 1.2

УДК

6 22 .0 2:5 3 1

 

В монографии изложена результаты теоретических и эксперимен­

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

Рассчитана на научных и инженерно-технических работников, преподавателей и студентов вузов.

Под редакцией чл.-кор . АН УССР Ф.А.Абрамова

Рецензенты:

д -р техн.наук В.А.Бойко, канд.техн.наук Л.П.Роменский

Редакция тёхнической'литературы

М------ 196-75

М221(0 4 )-75 '

© Издательство "Наукова думка", 1975 г .

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

вышением

газоносности пластов и пород, газообильности добычных

участков,

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

Повышение газообильности выработок приводит к необходимости увели­

чения количества подаваемого воздуха в

шахту, что,в свою очередь, ведет

к увеличению сечений горных выработок,

а значит, и к повышению затрат

на их проведение, поддержание и проветривание.

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

В настоящее время наиболее

надежный способ отработки групп сбли­

женных выбросоопасных пластов -

отработка их после предварительной раз­

грузки путем выемки защитных пластов.

В.В Л'етеревенко-

Советскими

учеными С.Г.Авершиным, И.М.Петуховым,

вым, Г.И.Грицко»

Б.В.Власенко,

И.А.Фельдманом созданы

основы теории за ­

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

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

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

Работы советских ученых в области сдвижения горных пород С.Г.Авар­ ийна, Д.А.Козаковского, И.Г.Лисицы, М.В.Короткова, С.П.Колбенкова, Р.А.Муллера, А.Н.Медянцева, А.Г.Акимова, М.А. Иофиса, В.Н.Земисева и других позволили не только установить основные закономерности,но и полу­ чить математическое описание этого сложного явления.

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

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

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

динамики угольных пластов в процессе их сдвижения при ведении очистных

.работ выполнены для основных бассейнов нашей страны Г.Д.Лидиным и А.Г.АЙруни.

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

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

ных пластах,

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

с ко­

личественной

оценкой каждого источника метановыделения в общем газо ­

вом балансе.

 

нау­

ки и

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

техники УССР,чл.-ков.АН УССР Ф.А.Абрамова,отдельные этапы (гл.1 и

Ш -

под руководством д-ра техн.наук И.Г.Лиоицы. Гл. 1,Ш,У, § 2-5

гл.1У

напиоаны кавд.техн.наук B .ll.l/лкеньким. Гл.П и § I глЛУ напиоаны канди­ датами технических наук В.И.Мякенькиы, Г.А. Шевелевым,В.С.Кулинячем.Те­ оретические исследования выполнены канц.техн.наук В.И.Мякеньким.

Авторы выражают глубокую благодарность д-ру техн.наук А.Т.Айруни, кандидатам технических наук-Г.Т.Василенко и А.В.Онищенко, инженерам Ю.А.Герасименко, Б.А.Бойко, Н.В.Молчанову, А.И.Воронкияу, Н.И.Меньшову за участие в выполнении отдельных этапов исследовании.

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

Г Л А В А I

СДВИЖЕНИЕ ГОРНЫХ ПОРОД И ЗЕМНОЙ ПОВЕРХНОСТИ ПРИ ДВИЖУЩЕМСЯ ЗАБОЕ ЛАВЫ

1 . Условия применимости интегральное функции Гауооа для описания процесса сдвижения горных пород

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

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

ющих трещин и пор, протекание физико-химических

процессов, связанных с

десорбцией газа

из пластов и передвижением

этого

га за в

сторону очист­

ной выработки.

Процесс десорбции метана из

угля

и пород

возникает вслед­

ствие

нарушения существующего равновесия системы уголь - газ (порода -

г а з ),

вызванного процессом сдвижения горных пород.

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

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

толщи пород,

окружающих очиотную выработку.

 

В результате исследований отечественных и зарубежных ученых установ­

лено, что для

описания мульды сдвижения земной поверхности при остановлен­

ном забое лавы и закончившемся процессе сдвижения горных пород может

быть

принята

интегральная

функция Гаусса /§ ,9 ,1 1 ,1 7 ,1 8 ,2 6 ,3 4 /. Кроме то­

го, в

работах

ркда ученых,

с различных позиций трактующих процесс сдвиже-

 

 

 

5

ння,

получены решения дифференциальных уравнений

для описания этого яв­

ления

в виде функции Гаусса 23 * 6 ,7 ,1 2 ,2 9 ,3 0 ,3 4 /.

Наряду с этим следует

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

ине могут изучаться в отрыве от нее.

Сцелы» изучения закономерностей сдвижения земной поверхности и мас­ сива-горных пород при движущемся забое лавы и математического описания

.8того процесса нами проанализированы данные наблюдений sa сдвижением 8емкой поверхности по 1-Й профильной линии реперов и площадной наблюда­ тельной станции, включающей пять профильных линий шахты "Першотравнева" в Западном Донбассе .выполненных кафедрой маркшейдерского дела Днепропет­ ровского горного института и наблюдений за сдвижением реперов четырех профильных линий,заложенных в массиве горных пород и на поверхности,а

Рис Л . Соответствие трансформированных кривых вертикальных сдвижений глубинных реперов интегральной кривой Гаусса:

I - 4 -я восточная лава пласта

/ / шахты "Ясиновекая-Глубокая" ком­

бината "Макеевуголь"; 2 - 5-я наблюдательная станция, I- я профильная

линия шахты "Першотравнева" ; 3 -

I- я и 2-я профильные линии, 2-й запад­

ный и восточный полевые откаточные штреки горизонта 915 м шахты

им.К.И.Поченкова в Донбассе.

 

также глубинных реперов,заложенных в скважины шахты "Яоиновокая-Глубо- кая" комбината "Макеевуголь" /2 6 /, выполненных сотрудниками шахты о уча­ стием авторов.

Для изучения поведения надрабатываемого выбросоопасного массива гор­ ных пород использованы данные' проведенных нами наблюдений за сдвижением

надрабатываемой толща горных пород по

1

в 2-й профильным линиям,зало­

женным соответственно во 2-м западном

в восточном полевых откаточных

штреках,

расположенных на

расстоянии

35-40 м ниже разрабатываемого

пласта

(горизонт 915

м) шахты им.Поченкова комбината "Макеевуголь"

в Донбассе.

Кривые вертикальных сдвижений реперов на различные моменты времени

(при различных положениях забоя лавы

относительно реперов линии)

транс­

формировались в

единичные кривые

(рис.1). При этом по оси абсцисс

от­

кладывалось отношение координаты

л

репера к

длине полумульды

I

по оси ординат

- отношение величины

оседания

£ репера' с координатой

лк величине максимального оседания £тах .

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

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

менно

предложена в

нашей стране Н.М.Ершовым / I I / и в Польше С.Кнотте

/1 7 /,

Основные ее

положения сводятся

к следующему.

 

 

Пусть имеется на поверхности или в массиве горных пород некоторая

точка

А

(рис.2),

на которую оказывает влияние разработка пласта мощ­

ностью

т .расположенного на глубине

Н .

Будем рассматривать сдвижения,

происходящие в верхнем полупространстве,

т .е . в полупространстве

с по­

ложительными значениями %. Конечное

оседание точки произойдет

при вы­

емке

части пласта

S'. Разобьем

отработанную часть пласта S на

а

эле­

ментарных площадок о площадью

Л Sj

каждая. Тогда выемка пласта в

пре­

делах à ô j

может вызвать

элементарное оседание

Ap .tоднако его

величина

будет зависеть не только от размера элементарной

площадки, но и от

удаления проекции точки

А

на плоскость пласта

от

элементарной

площад­

ки

A Si

 

, т .е .

гр .

-

dSj /

( л

t ÿ )

,где

J { x ,y ) -

некоторая

и з-

веатная функция координат (функция влияния).

 

 

 

 

 

 

 

 

 

Применив положение суперпозиции,

получим конечное

оседание

точки А

как предел, к которому отреыитоя сумма

/*/7

 

( х , у

;

при

стремле­

L AS- /

нии

AS-t

к нулю:

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

«

. о

 

 

Максимальное

 

 

 

 

 

р0

точки А произойдет

при

 

теоретическое оседание

условии выемки пласта на расстоянии +

оо

 

от проекции точки А на

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

 

 

 

 

 

 

 

 

 

 

 

+т»+то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

& * j

5 S t * y ) dxc(y-

 

 

 

 

U .2)

 

 

Для определения конечного вида функции влияния применим метод,пред­

ложенный

С.Кнотте

/1 7 /. Отметим,что

решения

задач, приведенных в

насто­

ящей работе,

относятся

к случаю

отработки горизонтально

и полого

зале­

гающего пласта и при способе управления кровлей полным обрушением.

 

Как отмечалось выше, кривые сдвижения земной поверхности и слоев

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

быть аппроксимированы интегральной функцией Гаусса (см .рис.1),

 

которая,

получается при

интегрировании функции Гаусса,

являющейся в нашем слу­

чае функцией

влияния.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно* раооте

/4 5 / функция Гаусса

имеет

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

■*'**

 

 

 

 

 

 

 

где л - абсцисса точки; А - некоторое вещественное число. Заметим,что

к 4 2* 2

-ор Этот интеграл для случая-плоской задачи соответствует (см. формулу(1.2))

максимальному значению оседания рд точки А при £ - 1 . Если величина мак­ симального оседания точки поверхности или слоя пород рд £ / , то она определяется из выражения

(1 .4)

Этому выражение соответствует функция влияния

 

«

h

 

U . 5 )

 

Я * )я?о ] р г е

.

 

Заменим кривую влияния f(x)

(см.рио.2) двумя отрезками прямых,обра­

зующих с осью

х равнобедренный треугольник. Площадь этого треуголь­

ника равна конечному оседанию

^ , а высота

-

максимальной ординате

кривой влияния

(1 .5 ),которая при условии z

°0

принимает значение

'

За основание этого треугольника принимаем длину полумульды «рассчи­ тываемую для земной поверхности по граничным углам ^54/. (Вопрос о рас­ чете длины полумульды для любого слоя подрабатываемого' массива горных пород будет рассмотрен ниже.) С учетом этого площадь треугольника,яв­ ляющаяся мерой конечного оседания точки, будет равна половине произве­ дения основания на высоту:

h l t

откуда

2 \Л Г

Подставляя значения h в функцию влияния (1 .5 ),

находим

-4 £ -Х *

(1. 6)

/(г )= Ц * -е т г *

Поверхность вращения этой кривой вокруг оси симметрии запишется в сле­ дующем виде:

J(x>y)=—f$ ~ e

z

( 1 . 7 )

ùz

 

 

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

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

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

Рассмотрим решение плоской

задачи, предусматривающей математическое

описание процесса сдвижения в

вертикальной плоскости

зсОх «перпендику­

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

 

Определим оседание точки А ,лежащей в верхней полуплоскости

л о г

(р » с .З ). Линию очистного забоя

ориентируем параллельно

оси Оу ,

Пусть

вдоль этой оси происходит полная подработка. Способ управления кровлей -

полним обрушением.

 

 

 

 

От точки Û вдоль оси

х

двинется забой лавы со скоростью

к • В

момент времени t забой

лавы

достигнет

точки С .пройдя путь

v t. Любо­

му моменту времени,в том числе

и моменту

t .соответствует определенная

вынутая площадь пласта. Каждая элементарная площадка указанной площади вызывает на поверхности элементарное оседание точки. Сумма действий каждой площадки в пределах всей вынутой части пласта вызывает конечное

оседание

точки

поверхности.

 

 

 

 

 

 

 

 

Выберем начало координат в точке О (на границе очистной выработки

со

стороны разрезной

печи),

ось

х

направим вдоль

движения забоя лавы.

В этом случае

согласно

(1 .1)

для точки А существует некоторая функция

/

(х)

характерная,

тем,

что

площадь,

ограниченная

кривой

/(& )

и отрезком ab

на

оси <£.является

мерой конечного

оседания этой точки

при выбмке пласта в пределах

ОС.

Отрезками

Оа и

ÔC

здесь

учтен фак­

тор зависания пород основной кровли

у границ очистной выработки. Причем

отрезок Оа характеризует величину

зависания пород основной кровли

над

остановленным и простоявшим некоторое время забоем лавы /1 8 ,3 4 /, а

ЪС-

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

ботах /2 3 ,3 4 /. £ геометрическом

смысле

отрезки

Оа и ÔO равны величи­

нам смещения точки перегиба в мульде

сдвижения относительно неподвиж­

ной границы очистной выработки или движущегося забоя лавы.

 

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

плоскости хОу,

можно определить из

следующего выражения

* :

 

 

 

 

J

е

Lt

4А .

 

 

(1 .8 )

Производя

 

а

 

 

 

 

 

по формуле

замену переменных и выполняя интегрирование

Ньютона -

Лейбница,

полу чаешГ ^

 

 

 

Jг

 

 

 

!(* ,* .Л

 

 

 

-

 

(1>9)

где р0 -

величина

максимального

оседания земной поверхности при пол­

ной подработке, м;

/^ “ Длина полумульды на горизонте

z

, м;

* Вид функции влияния изменился за очет смещения ее графика для общего положения точки А относительно начала координат.