книги / Сдвижение и дегазация пород и угольных пластов при очистных работах
..pdfХу X - координаты |
точки земной поверхности или |
массива горных по |
|||
род, в которой определяются параметры сдвижения, ы; |
у - скорость |
подви |
|||
гая ия очистного забоя |
лавы, ы/мес,; t - момент времени,в |
который оп |
|||
ределяются параметры |
сдвижения данной точки, |
мес. (начало |
отсчета |
вре |
|
мени выбрано с момента проведения разрезной |
печи); |
d^ - величина |
зави |
сания пород основной кровли у границ очистной выработки над движущимся забоем лавы, м; d - величина зависания пород основной кровли над не подвижной границей очистной выработки. (в данном случае над границей вы работки со стороны разрезной печи), м; Ÿ - табличная интегральная функ ция Гаусса.
Формула (1.9) представляет собой математическую модель динамической мульды сдвижения вдоль главного сечения для любого подрабатываемого слоя горных пород или для земной поверхности при условии полной подработки по оси у . Она позволяет определить величину оседания любой точки,нахо
дящейся в верхней полуплоскости хО% .для любого момента времени при
движении лавы со скоростью |
V . |
|
|
|
|
|
Для решения пространственной задачи рассмотрим разрабатываемый плаох |
||||||
угля мощностью т , |
залегающий на |
глубине f t |
от |
земной поверхности» |
За |
|
бой лавы движется со скоростью |
у от разрезной |
печи вдоль оси * . |
Спо |
|||
соб управления кровлей - полным обрушением. Как и в случав плоокой |
зада |
|||||
чи, над границами очистной выработки происходит |
зависание пород основной |
|||||
кровли. Исходя из |
соответствия |
кривой сдвижения над движущимся забоем |
||||
лавы интегральной |
функции Гаусса |
(см .р ас.!) |
в качестве функции влияния |
|||
принимаем предложенную наш |
ранее /3 7 / функцию |
(1*7) |
|
Определим величину оседания точки |
А , заданную координатами |
я , у х,ъ |
||||||
момент времени |
é |
. Точка А находится : на |
земной поверхности |
или на |
|
|||
горизонте z |
(в верхнем полупространстве) |
под действием отработки |
плас |
|||||
та мощностью |
т . Начало координат |
выберем на середине линии,ограничива |
||||||
ющей очистную выработку со-стороны |
разрезной печи (точка 0 на р и с .3 ) . |
|||||||
По аналогии |
с |
(1 .8) |
запишем выражение для определения величины |
оса- |
||||
дания точки |
А |
в |
момент времени t |
; |
|
|
|
|
|
|
|
vi-d. «-S, |
|
|
|
|
|
|
|
|
|
% ( л - * ) 2 (* -У )2 |
|
|
||
|
|
|
|
L* |
|
etA etУ. |
(1 .10) |
|
|
|
|
* v |
-с |
|
|
|
|
|
|
|
|
|
|
|
Сделав замену переменных и выполняя интегрирование,получим
х
|
|
|
|
|
( I . I I ) |
где |
. |
л |
, |
~ РасчетныЙ размер очистйой выработки по оси |
|
\ ~ ~ 2 — |
" |
||||
м; |
д2 - |
фактический |
размер очистной выработки по оси у , м. |
Это выра |
жение представляет собой математическую модель, динамической мульды сдви жения, образующейся на земной поверхности или на горизонте х в заданный
момент времени. Она дозволяет определить величину оседания точки о
координатами |
€, |
у , z |
земной поверхности или массива горних пород в |
момент времени |
t |
ори движении очистного забоя со скоростью Y . |
НО3
Рис.4 . Экспериментальные (спдошнне) и расчетные (штриховые) кривые сдвижений и деформаций земной поверхности по 7-й профильной линии шахты "Першотравнева” .
Результаты изложенных теоретических исследований хорошо согласуются с данными экспериментальных наблюдений за сдвижением земной поверхности (рис .4 -6).,
При подработке сооружений не редки случаи.когда они после ухода за боя лавы и затухания процесса сдвижения попадают в пределы плоского два мульды сдвижения. Производя расчет деформаций поверхности под такими сооружениями па момент отработки этой лавы и затухания процесса сдвиже ния, можно сделать вывод,что значения ожидаемых деформаций допустимы
дая данных сооружений. Однако в процесс в подработки по мере подхода ж удаления забоя лава сооружения могут претерпевать значительные де
формации. В связи с этим возникает эадача определения |
значений сдвиже |
||
ний и деформаций в заданной |
точке земной поверхности |
для любого момен |
|
та времени* |
Задача прогноза |
ожидаемых деформаций массива горных пород |
|
возникает и |
при решении вопросов охраны подземных сооружений, при веде- • |
||
нии горных |
работ вблизи затопленных ранее. выработок или под плывунами. |
Рис.5. Вертикальные разрезы |
вдоль |
1 ,2 ,3 и 4-й профильных |
линий реперов при отработке |
пласта |
и экспериментальные |
(сплошные) и расчетные (штриховые) |
кривые вертикальных |
|
сдвижений. |
|
С процессом сдвижения горных пород тесно связан и процесс естествен ной дегазации. Поэтому надежные методы расчета газообильности добычных участков могут быть получены только с учетом закономерностей процесса сдвижения горных пород. Подтверждением этого является работа Б.Г.Тарасо ва /5 9 /, где разработан метод определения пористости и проницаемости по род междуплас товой толщи по величине нормальных и касательных смещений горных пород и угля в зонах подработки и надработки. На основе этого ме тода получены дифференциальные уравнения, характеризующие процесс есте ственной дегазации междупластовой толщи и сближенного пласта.
Уравнение ( I , И ) позволяет перейти к определение величин наклонов, кривизны, горизонтальных сдвижений и деформаций для любого слоя горных
пород. |
Так,в точках земной |
поверхности или. заданном координатой z слое |
|
горных |
пород в плоскости |
a û z |
величины наклонов будут определяться |
по формуле |
|
|
У1
Кривизна в точке с координатой X в момент времени t определяет* ся так:
•(vt-X-x)#
Для нахождения величины горизонтальных сдвижений и деформаций вос пользуемся зависимостями,предложенными С.Г.Авершиным:
м ~ ’
где кп - некоторый коэффициент. |
Аналогичным(путем можно найти парамет |
|
ры сдвижения в любом сечении, параллельном оси * |
или у . |
|
Предложенный метод расчета хорошо согласуется |
с данными натурных |
наблюдений. Это следует,например,из сравнения показанных на рис.4 экспе риментальных и расчетных кривых сдвижений и деформаций земной поверхнос ти для второй профильной линии шахты "Першотравнева" в Западном Донбассе. При расчете сдвижений оседание каждого репера вычислялось как сумма осе даний от подработки каждой лавой в отдельности /377 .
Изложенный выше метод расчета параметров сдвижения земной поверхно сти позволяет также определять моменты времени, когда здания и сооруже ния будут претерпевать критические деформации,что позволит своевременно принимать соответствующие меры для их охраны.
п
3 . Определение параметров сдвижения подработанной толщи горных пород
Одним из основных параметров, характеризующих физико-механические свойства горных пород, а следовательно, предопределяющих характер процесса сдвижения при ведении очистных работ,являются граничные углы /6 ,3 (7 . Длину полумульды для земной поверхности,входящую в качестве основного параметра в полученные выше математические модели динамичес кой мульды сдвижения и в методы расчета с применением типовых кривых
/§ 7 , |
определяют по |
глубине разработки |
и граничным углам /5 4 / : |
|
|
|
|
( 1* 12) |
|
где /У - глубина разработки пласта, м; |
<Р- угол полных сдвижений, град; |
|||
â0 |
- граничный угод,град. |
|
|
|
|
При определении параметров сдвижения массива горных пород существую |
|||
щими аналитическими |
способами расчета /5 ,1 1 ,1 4 ,5 4 7 исходят из предпосы |
|||
лок, что зона сдвижения массива пород оконтурена граничными углами |
и ли |
|||
нией земной поверхности. Длина полумульды в слое, лежащем в кровле |
раз |
рабатываемого пласта, принимается равной нулю,что следует из уравнения (1.12) при Н*%*0, Однако практикой ведения очистных работ установлено,
(,мм
100
200
300
400
т
40 |
20 |
0 |
-20 |
-40 |
-60 |
-ВО |
-1001,м |
|
|
Положение забое лавы |
|
|
|
Рис.6. Экспериментальные (сплошные) и расчетные (штриховые) кривые вертикальных сдвижений глубинных реперов f i l и
fin (4-я восточная лава пласта г1, шахты "ЯсияовскаяГлубокая").
что зона сдвижения в подработанном массиве горных пород, расположенном вблизи разрабатываемого пласта,не может быть оконтурена граничными уг лами, она захватывает более обширные области впереди и позади забоя ла вы / о / . Последнее объясняется действием опорных нагрузок на пласт,его деформацией впереди лавы. Проведенные нами экспериментальные исследова ния по изучению процесса сдвижения земной поверхности и массива горных пород в условиях шахты "Ясиновская-Глубокая" /367 также подтвердили зго положение.
Сдвижение массива горных пород в условиях указанной шахты изучалось по профильным линиям, заложенным в I-м восточном вентиляционном штреке коренной -пявн пласта и находящемся на ПО м выше отрабатываемого пласта. Протяженность подземных профильных линий составила 120 и 140 м.
Вертикальные, разрезы профильных линий и план горных работ приведены на ри с.5 *7.
Jü
Рио.7 .Расположение реперов наблюдательных отанций при отработке пласта
16 (шахта |
"Ясиновская-Глубокая'', Донбасо). |
|
|
Рио.8. Схема расположения глубинных,реперов, |
заложенных с транспортного |
||
штрека 4-Й |
восточной лавы, пласта li |
(шахта |
"Яоиновокая-Глубокая , |
Донбасс). |
|
|
|
С целью исследования поведения массива и спутников пласта,располо женных вблизи очистной выработки, с транспортного штрека 4-й восточной
лавы пласта |
1$ в породах кровли, вмещающих спутники |
1$ • и |
лг; ,на |
|
расстоянии |
16 и 32 м от контура штрека и в почве выработки на глу |
|||
бине 9 м были заложены глубинные репера ftl, f t l |
и |
ftШ |
(рис.8). |
Величина и характер вертикальных смещений глубинных реперов изучался по
мере подвигания забоя лавы. Одновременно скважины, |
пробуренные под репе |
|||||
ра, использовались для замера скорости газовыделения |
из спутников |
I* |
||||
и ni] по мере подвигания |
забоя |
лавы. |
|
|
|
|
Анализ экспериментальных кривых вертикальных сдвижений по этим репе |
||||||
рам (см.рис.6) показывает,что |
качественно они подобны кривым сдвижения |
|||||
земной поверхности,которые как |
известно /3 ,2 6 /,могут |
быть описаны инте |
||||
гральной |
функцией Гаусса |
/ 5 / . |
|
|
|
|
Как |
отмечалось ранее, |
трансформирование указанных |
кривых в |
единич |
ную кривую и аппроксимация их интегральной кривой Гаусса дает возмож ность судить о достаточно высокой сходимости этих кривых. Это экснериментально подтверждает правомерность применения функции влияния /267 для описания процесса сдвижения массива горных пород при движущемся забое лавы.
Наблюдения 9а сдвижением глубинных реперов позволили также экспери ментально определить фактические длины полумудьд при движущемся забое
для горизонтов 16 и 32 м, которые |
оказались равными соответственно 120 |
и 130 м (ри с.9). Длины полумульд, |
рассчитанные по граничным углам для |
этих горизонтов, будут составлять |
всего лишь 15 ж 31. Естественно,что |
использование граничных углов для определения длины полумульды в масси ве горных пород и расчета по полученному значению 1г параметров сдви жения этого массива приведет к существенным ошибкам. В частности, длина зоны сдвижения вблизи пласта будет в 6-10 раз меньше фактической. Это обстоятельство необходимо учишвать не только при расчете сдвижений от дельных слоев, но' в при расчете дегазации сдвигающегося массива горных пород.
Рис.9 . График изменения длины полумульды при удалении от пласта.
Проведенные на шахте "Ясиновская-Глубокая" экспериментальные исследо вания позволили подучить информацию о параметрах сдвижения массивагорвш пород на различных расстояниях от разрабатываемого пласта,включая земную поверхность,а именно на расстояниях 16,32,110 и 307 м. В частности, по этим данным установлен линейный характер и вид уравнения, описывающего изменение длины полумульды с-удалением по нормали от разрабатываемого пласта. Для поверхности это уравнение имеет вид
L /fl(ctg< p
Уравнение длины полумульды в общем виде для любого горизонта можно записать
Lt *ka+ 6. |
(1.13) |
. |
» |
• |
Для определения к и b рассмотрим рис.9,согласно которому можно соста вить следующую систему уравнений:
|
297 - к SO? |
Ь |
|
откуда к*Ц$', Ь яПО. |
120 = к Гб * Ъ |
|
|
Тогда искомое уравнение (1.13) |
будет иметь вид |
||
|
LЛ* 4 6 x +no [м ] . |
(1.14) |
|
Таким образом, проведенные экспериментальные исследования в уело» |
|||
виях шахты "Ясииовская-Глубокая" |
позволили установить не только общий |
характер процесса сдвижения земной поверхности и массива горных пород, но и получить конкретное значение длины кривой сдвижения 1Ъ доя любого
горизонта. |
|
(I . I4) значений Lx ,а также |
Подстановка |
полученных по формуле |
|
других иоходных параметров в уравнение |
динамической мульды сдвижения |
|
( I . I I ) позволяет |
находить для любого момента времени величину вертикаль |
ных оседаний в любом горизонтальном слое пород,расположенном на расстоя
нии |
х |
от разрабатываемого |
пласта. |
|
||
|
Сравнение экспериментальных |
кривых оседания реперов профильных ли |
||||
ний, |
расположенных |
на поверхности |
и на расстоянии ПО, 32 и 16 м от плас |
|||
та,с |
рассчитанными |
по формуле |
(1 ,1 1 ) (см .рис.4 ,5 ,6 ) ,подтверждает правиль |
|||
ность полученных уравнений, а |
также их надежность. |
|||||
|
При выводе уравнения мульды сдвижения земной поверхности или произ |
|||||
вольного |
слоя горных пород интегрирование функции влияния производилось |
|||||
- не в |
пределах всей вынутой площади пласта, |
а с учетом величины зависания |
||||
U |
пород основной кровли у границ очистной выработки* |
|||||
|
Величина зависания пород кровли для земной поверхности определяется |
|||||
как |
функция глубины /6 ,1 8 /: |
|
где |
л - некоторый коэффици |
||
ент, завиоящий от физико-механических.овойотв породы. |
||||||
|
В работе /2 6 / получено следующее выражение для определения величины |
|||||
зависания основной кровли над неподвижной границей очистной выработки: |
||||||
|
|
|
|
|
|
(1.15) |
где |
V - граничный угод и угол полных сдвижений. В свою очередь,гра |
ничные углы характеризуют физико-механические свойства сдвигающейся тол щи горных пород.
Для движущегося забоя лавы в первом приближениивеличину зависания пород основной кровли над аабоем лавы необходимо рассматривать как функ цию скорости подвигания забоя и физико-механических свойств пород кров ли. В первом приближении ее можно принять примерно в два раза большей
величины зависания над неподвижным забоем. Эю необходимо учитывать при
расчетах сдвижений по формулам (1*9)» |
( I . I I ) . |
|
|
|
|
||||||
|
Так как |
величина |
зависания |
d |
является |
функцией расстояния от |
|||||
разрабатываемого пласта, |
запишем выражение для |
ее |
определения |
в общем |
|||||||
виде: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d - k x , |
|
|
|
|
(1 .16) |
|
где |
к - коэффициент, |
учитывающий физико-механические свойства |
пород; %- |
||||||||
расстояние |
от |
разрабатываемого пласта до земной |
поверхности |
или рас |
|||||||
сматриваемого слоя горных пород, м. |
|
|
|
|
|
||||||
|
Из формулы |
(I . I5) |
находим |
|
|
|
|
|
|
||
|
|
|
|
|
|
ciçV -etffd , |
|
|
|
|
|
|
|
|
|
|
н |
|
Z |
|
|
|
|
Тогда величина зависания пород основной крови для любого |
слоя, находяще |
||||||||||
гося |
на расстоянии £ |
от |
разрабатываемого пласта, |
определится |
так: |
||||||
|
|
|
|
|
, |
ctg4> -ctgâa |
|
|
|
( I .17) |
|
|
|
|
|
|
а* |
|
2 |
|
|
|
|
|
Следует |
отметить,что |
величина |
d |
- одна из |
основных характеристик |
|||||
сдвигающейся |
толщи горных пород.В |
частности, удвоенная величина зависа |
|||||||||
ния |
( 2 d ) |
|
характеризует минимальный пролет очистной |
выработки при |
|||||||
отходе от разрезной печи, при котором начинается процесс |
сдвижения зем |
||||||||||
ной поверхности. Как видно из уравнения (I . I7 ) |
.величина |
зависания для |
данного слоя является линейной функцией глубины или расстояния от разра батываемого пласта до данного слоя.
Величину зависания пород основной кровли можно рассматривать как по-, казатель, характеризующий затухание действия очистной выработки на сдви жение горных пород при удалении по аормали от разрабатываемого пласта.
В частности, |
для слоя пород, |
расположенного непосредственно над разра |
|
батываемым пластом, d=0, в |
сдвижении этого сдоя принимает участие вся |
||
ширина очистной выработки. |
При удалении от пласта на расстояние |
х |
|
величина d |
возрастает и в с движении принимает участие уже не вся |
ши |
рина выработки, а уменьшенная на величину d с |
каждой стороны. По форму |
|
ле ( I . I 7 ) можно рассчитать |
предельное расстояние по нормали,при котором |
|
затухает действие очистной |
выработки размером |
а : |
|
2d |
|
ctg ф - ctç d0
Так, для условий Донецкого бассейна 2,32 d , что качественно сог ласуется с данными моделирования на вальцмассе процесса сдвижения гор ных пород /4 8 / .