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

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

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

Ху 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.

Рио.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. Естественно,что

использование граничных углов для определения длины полумульды в масси­ ве горных пород и расчета по полученному значению параметров сдви­ жения этого массива приведет к существенным ошибкам. В частности, длина зоны сдвижения вблизи пласта будет в 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)

Таким образом, проведенные экспериментальные исследования в уело»

виях шахты "Ясииовская-Глубокая"

позволили установить не только общий

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

горизонта.

 

(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 / .