книги / Сдвижение и дегазация пород и угольных пластов при очистных работах
..pdfАКАДЕМИЯ НАУК УКРАИНСКОЙ ССР
ИНСТИТУТ г е о т е х н и ч е с к о й м е х а н и к и
В. И. Мякенький
СДВИЖЕНИЕ И ДЕГАЗАЦИЯ ПОРОД
УГОЛЬНЫХ ПЛАСТОВ ПРИ ОЧИСТНЫХ РАБОТАХ
КИЕВ « ПАУКОВА ДУМКА » 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 |
принимает значение |
'
За основание этого треугольника принимаем длину полумульды 1Х «рассчи тываемую для земной поверхности по граничным углам ^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 |
, м; |
* Вид функции влияния изменился за очет смещения ее графика для общего положения точки А относительно начала координат.