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

книги / Сдвижение горных пород и защита подрабатываемых сооружений

..pdf
Скачиваний:
8
Добавлен:
13.11.2023
Размер:
24.46 Mб
Скачать

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

3.2.3.

Кровля как среда, свойства которой описываются реологической моделью Кельвина

До сих пор рассматривался процесс опускания кровли в очистной выработке в конечной стадии ее выемки, когда горные работы уже не ведутся, т. е. прини­ малось, что оседание происходит сразу же после воздействия нагрузки. Однако, как известно, процесс деформирования нагруженного слоя пород или поддер­ живающей его опоры завершается лишь по истечении некоторого времени. По­ этому следует рассматривать опирающиеся на закладку породы кровли как рео­ логическую среду, для которой зависимость деформаций от времени может быть описана той или иной реологической моделью, например моделью вязкоупругого тела Кельвина, представляющей его в виде сочетания упругой пру­ жины Гука с вязким поршневым элементом Ньютона (рис. 34). Напряжения в элементе Гука пропорциональны деформациям Е е, а в элементе Ньютона за­ висят от скорости деформирования r\ deldt, и, следовательно, общее напряже­ ние, возникающее в теле Кельвина при его деформировании

о = Ее-\-у) ((de/dt),

 

(29)

где ц — коэффициент

вязкости.

0, находящегося под действием

Для тела Кельвина с момента времени t =

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

S X

б

нагружение

Рис. 34.

Реологическая модель Кельвина (а), состоящая из соединения упругой пружины Е и вязкого поршневого элемента (вязкость ц), и отвечающие этой модели диаграммы деформаций

н в)

в

б

£

Рис. 35.

Распределение горного давле­ ния над зоной опорного дав­ ления и выработанным прост­ ранством при подвигании очи­

стного

забоя на расстоянии

х = vt

и изменение давления

в неподвижных точках

и Р

I — зона опорного давления;

II

выработанное пространство, запол­ ненное закладной

предельному значению е = о/Е, причем его изменение во времени описывается выражением

e = -|-(l-e-(£/Ti)<).

(30)

Однако в действительности при проведении очистной выработки давление пород кровли на пласт или на закладку в некоторой определенной точке не остается постоянным (рис. 35); возникающие сжимающие напряжения возра­ стают в зоне опорного давления, падают до нуля у очистного забоя и вновь возрастают над заложенным выработанным пространством до величины р = = yh, причем зависимость расстояния х между рассматриваемой точкой и за­ боем от времени t и скорости подвигания очистных работ v может быть выра­ жена уравнением х = L + vt для зоны опорного давления и уравнением х = = vt для выработанного пространства с закладкой. Слагаемое L, имеющее знак минус, по абсолютной величине равно ширине зоны опорного давления. Поэтому в преобразованное уравнение напряжений в теле Кельвина (29) дол­ жна быть подставлена переменная величина сжимающего напряжения б (я), изменяющаяся в зависимости от расстояния х между рассматриваемой точкой и забоем. Тогда для определения скорости деформирования в зоне опорного давления получим выражение

т) -^г р [ } + ( ^ “) / *e7n:cJ — Ее при x = L + vt9

(31)

где pj^l + ("Т7") /2em*J представляет собой переменное напряжение сг (я)

из уравнения (25).

Для определения скорости деформирования в поле закладки получим выражение

dz

т] -JJ- = р ([ — е“пх) —Ег при x = vt,

где р (1 — е~пл:) представляет собой напряжение а (х) из уравнения (24).

В результате интегрирования выражения (31) получим для зоны опорного давления

ЕL-x

ек = К е * » + -§-

(33)

----+ 771V

Я

где

представляет собой зависящую от скорости подвигания очистных работ дефор­ мацию а (v)/E.

Интегрированием выражения (32) и добавлением величины предваритель­ ной конвергенции е*тах получим для поля закладки (индекс v обозначает закладку)

Е

X

 

е0 = Ке ч

• -f -£ - 1

(34)

где

также представляет собой зависящую от скорости подвигания очистных работ деформацию a (v)/E>

Значение входящего в выражения (33) и (34) множителя К получается из условия, что при t — 0, т. е. при х ^ L в зоне опорного давления (точка Р г)

и при х =

0 в поле закладки (точка Р 2), зависящая от времени деформация так­

же равна нулю, т. е. е = 0. Из равенства

 

 

Е

О

 

(35)

Ке ’> 0 = 0

 

 

 

следует, что при е° = 1 значение К равно нулю. Для точки с абсциссой х = L

при условии, что L имеет весьма большую

абсолютную величину и что е771* =

= emL

0, выражение (33) будет иметь

вид

 

 

 

 

(36)

Максимальная вертикальная деформация угольного пласта у очистного забоя (предварительная конвергенция) при скорости подвигания очистных работ v составляет

л\ _ Р

Ьmax-----jjr

Рис. 36.

Распределение вертикальных составляющих гор­ ного давления и деформаций или оседаний в очи­ стной выработке для различной скорости подвигания забоя в соответствии с формулами (33) и (34):

I— зона опорного давления; II — выработанное прост­

ранство, заполненное закладкой; я0 — прекращение очи­ стных работ

Для поля

закладки в

точке с абсциссой х = О (точка Р 2) из уравнения

(34) получим

 

 

щах»

 

(38)

а для точки с

абсциссой х =

оо

e v — “Ц - +

max*

(3 9 )

На рис. 36 показаны кривые изменения напряжений Я деформаций для различных скоростей подвигания очистных работ. ПлощадЯ и 4 2, заклю­ ченные между кривыми изменения напряжений и параллельной оси абсцисс прямой, соответствующей давлению от веса пород покрывающей толщи, для одинаковых интервалов времени неодинаковы, поскольку равновесие напря­ жений должно установиться только после того, как пройдет процесс деформиро­ вания, т. е. будут преодолены сопротивление трения и силь* Сцепления.

3.2.4.

Эмпирические формулы для расчета Конвергенции в выработанном пространстве

Во многих горнопромышленных районах при помощи корреляционного ана­ лиза исследовалась взаимосвязь результатов измерения опусканий кровли или конвергенции, выполненных в выемочных штреках и очи£Тцых выработках, с расстоянием от забоя, глубиной разработки, вынимаемой Мощностью пласта и скоростью подвигания очистных работ. При применении 9тцх эмпирических зависимостей было замечено, что оседания пород, наблюдаюхДИеся в выемочных штреках у границ выемочного поля, чаще всего бывают менв*Че, чем в главных сечениях очистной выработки по простиранию или падений0 пласта, а также не всегда свободны от влияния динамики проведения штрека* Кроме того, вели­ чина конвергенции с, которая, как известно, представляв^ собой величину сближения двух точек, расположенных по одной вертикали в кровле и почве выработки, и является суммой опускания кровли ски подняла почвы сп может

Рис. 37.

Зависимость конвергепции с и скорости кон­ вергенции cd с момента подработки наблюда­ тельной точки А от расстояния х между этой точкой и забоем:

I — призабойное

пространство; II

— выработанное

пространство; 1 — маркированные

наблюдательные

точки; 2 — стойка

крепи

 

С, мм

О

2 Ь 6 8 W X, м

U

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

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

с = Ylg (а+Ра)

(40)

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

Скорость конвергенции, т. е. приращение cd конвергенции, между двумя точками (мм/сут), достигает наибольшего значения на расстоянии от 2 до 35 м позади фронта очистных работ, в точке перегиба кривой конвергенции (рис. 37). Чем жестче кровля и чем выше скорость подвигания очистных работ, тем даль­ ше отодвигается от очистного забоя точка перегиба кривой конвергенции.

Для условий Верхне-Силезского каменноугольного бассейна опускание кровли впереди движущегося забоя на расстоянии х (в м) с достаточной точ­ ностью описывается уравнением

wn = и?0е~Ьл\

(41)

где w0 — опускание кровли (в мм) у кромки угольного

пласта; Ь — коэффи­

циент, зависящий от способа закладки (при выемке с

обрушением Ъ= 0,09

и при гидравлической закладке выработанного пространства b = 0,07). Опу­ скание кровли в поле закладки вычисляется по формуле

в которой М — вынимаемая мощность пласта, а а и Ъх — коэффициент*»!, характеризующие деформационные свойства пород кровли в призабойном пространстве и в заложенном выработанном пространстве.

В условиях Донецкого бассейна зависимость опускания кровли w в при­ забойном пространстве до края зоны обрушения или закладки от скорости подвигания очистных работ v (м/сут) выражается показательной функцией

«’приз =aM be$v,

(43)

т. е. опускание кровли в призабойном пространстве пропорционально вынимае­ мой мощности пласта М и ширине призабойного пространства Ь. Эмпирические коэффициенты а и 0 составляют, например, при выемке с обрушением и крепле­

нием призабойного пространства шагающей крепью: а = 74,2

и р = 0,13.

Во многих горнопромышленных районах Европы была установлена взаимо­

связь между с у т о ч н о й к о н в е р г е н ц и е й cd (м/сут),

измеренной

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

выражен

отношением c j v d, зависящим от глубины разработки

/г, коэффи­

циента оседания а и вынимаемой мощности пласта М :

 

tga =

= 0,2 у^/г (аМ)л/«.

(44)

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

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

Увеличение конвергенции с течением времени в условиях остановленного забоя, знание которой важно для расчета развития процесса сдвижения во вре­ мени (см. рис. 21), удалось наблюдать до сих пор в условиях Рурской области только на протяжении 60 ч. По результатам этих исследований установлено, что конвергенция пары наблюдательных точек, удаленных от забоя на расстоя­ ние х , в течение времени t возрастает в соответствии с показательной функцией

спрнз = A tB-\-CxD,

(45)

пока не достигнет статического конечного значения. Точные значения коэф­ фициентов и показателей для условий описываемых наблюдений установить не удалось — ориентировочно А = 0,17; В = 0,5, С = 0,5 и D = 0,6. Для статистического определения кривой конвергенции или опускания кровли на всем протяжении от зоны опорного'давления до середины выработанного про­ странства, зависящей от влияния большого числа различных факторов, необхо­

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

Однако для многих методов расчета сдвижения горных пород и земной поверхности достаточно заменить кривую опускания кровли ступенчатой ли­ нией, состоящей из отрезков длиной 25 или 50 м, выражая высоту этих ступе­ ней, т. е. разность оседаний на данном интервале, в процентах от максималь­ ного оседания (см. рис. 8).

3.3.

Расчет сдвижений горных пород

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

Т А Б Л II Ц А

Предмет

Эмпирические методы

Модель,

положенная

Искомые величины

расчета

п оси

ну метода

Конверген­

Формулы типовых кри­

1. Упругий прогиб бал­

Конвергепция в различ­

ция в вые­

вых опускания кровли

ки па податливом основа­

ных

точках

или

опуска­

мочном

(недостаток:

отсутствие

нии

 

 

 

несвязных

ние

кровли

и

поднятие

поле

универсальности)

2. Массив

почвы

 

 

 

 

 

 

 

пород (например,

песча­

 

 

 

 

 

 

 

 

 

ная модель)

 

 

тело

 

 

 

 

 

 

 

 

 

3. Вязкоупругое

 

 

 

 

 

 

 

 

 

(модель

Кельвина)

(недо­

 

 

 

 

 

 

 

 

 

статок: не учтена неупру-

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

но распределение нагруз­

 

 

 

 

 

 

Сдвижение

Метод

интеграцион­

ки по площади)

 

 

Вертикальные

 

и го­

1. Модели

сплошной

 

массива

ных сеток (недостаток: не­

среды для упругих и вяз-

ризонтальные

составля­

горных

известно влияние геоме-

коиластичных

материа­

ющие

сдвижения

точек

пород

ханических

характери­

лов;

Стохастическая сре­

горного

массива;

величи­

 

стик)

 

2.

ны в конечной стадш! про­

 

 

 

да;

 

 

 

конечных

цесса и в промежуточные

 

 

 

3. Метод

моменты

 

 

 

 

 

 

(упругих)

элементов (не­

 

 

 

 

 

 

 

 

 

достаток:

решения

полу­

 

 

 

 

 

 

 

 

 

чены только для плоской

 

 

 

 

 

 

 

 

 

задачи,

не

учитываются

 

 

 

 

 

 

 

 

 

разрывные деформации и

 

 

 

 

 

 

 

 

 

фактор времени,

слоисто­

 

 

 

 

 

 

 

 

 

сть

учитывается

 

только

 

 

 

 

 

 

 

 

 

в методе конечных элемен­

 

 

 

 

 

 

 

 

 

тов)

 

 

 

 

 

 

 

 

 

 

 

 

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

Методы расчета можно подразделить на эмпирические и основанные на мо­ делях массива горных пород (табл. 4). Искомыми величинами в этих методах являются вертикальные и горизонтальные составляющие сдвижения, из кото­ рых могут быть получены величины линейных деформаций пород в вертикаль­ ном направлении (см. рис. 19), а также горизонтальные растяжения и сжатия (см. рис. 23). Одним из перспективных методов расчета является метод, основан­ ный на применении палеток (интеграционных сеток). Из методов, основанных на моделях массива горных пород, можно назвать многочисленные теоретиче­ ские решения для зоны сдвижения, основанные на теории упругости (при по­ мощи сосредоточенных нагрузок, дифференциальных уравнений и конечных элементов). В ПНР разработаны методы, в которых массив горных пород рас­ сматривается как среда из несвязанных между собой частиц, движение которых подчиняется законам математической статистики, а в СССР — методы, в кото­ рых массив горных пород уподобляется вязкопластичному телу.

3.3.1.

Метод интеграционных сеток

Графический расчет сдвижений горных пород над очистными выработками при горизонтальном залегании разрабатываемого пласта, выполняемый при помощи интеграционных сеток (палеток), до сих пор остается в ФРГ единствен­ ным применяемым в широком масштабе методом расчета сдвижений. При этом методе условно принимают, что весь массив горных пород над рассматриваемой точкой снят, и ведут расчет так же, как для точек земной поверхности, — ин­ тегрированием дифференциалов оседаний или горизонтальных сдвижений, обусловленных выемкой элементарной площади очистной выработки в пределах области влияния [9]. Таким образом, метод, разработанный для расчета сдви­ жений земной поверхности, распространяется на любые горизонты горного мас­ сива, причем принятый для земной поверхности метод разбивки интеграционной сетки сохраняется и только размер этой сетки уменьшается в соответствии с высотой рассматриваемой точки породного массива над очистной выработкой (рис. 38).

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

 

 

;

Рис. 38.

пород

в

Схема к расчету оседания горных

точке Р при помощи интеграционной сетки:

Р 2

7 — V — зоны; 1 — земная поверхность;

2—расчет­

ный горизонт; з — очистная выработка

 

 

для которой производится расчет, вызывает опускание этой точки на 1/5 пол­ ного оседания. Для облегчения графического расчета каждая из этих кольце­ вых зон разбивается на восемь одинаковых секторов и, следовательно, вся пло­ щадь полной подработки разбивается на 40 равноценных частей (ячеек сетки), каждая из которых оценивается в 2,5% полного оседания.

Расчет для выемочного участка любой формы ведется следующим образом: вычерченная на листе кальки интеграционная сетка диаметром 2R = 2 h ctg у (площадь полной подработки) накладывается на план горных работ по разра­ батываемому пласту (горизонту) и ее центр совмещается с спроектированной (отвесно) на поверхность пласта точкой для которой производится расчет оседания. Затем подсчитывают число п ячеек сетки, оказавшихся в пределах контура отработанной части пласта и, умножив полученное число на цену

каждой

ячейки (2,5%), получают коэффициент влияния е по

формуле

е =

0,025тг,

(46)

откуда ожидаемое оседание точки Р породного массива

 

Vz — evг п.

(47)

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

в свою очередь, от его мощности. Только геометрические значения граничного угла у, высоты рассматриваемой точки h l и координат очистной выработки определяют положение точки Р 1? в которой частичное оседание v2 на вертикаль­ ной линии Р Р Гпереходит в полное оседание vzn. Деление сетки на зоны в зави­ симости от глубины, при котором могли бы учитываться стратиграфическое строение породного массива и давление от веса пород покрывающей толщи, еще не может быть выполнено из-за недостаточного количества эксперименталь­ ных данных.

Так называемая зона полного оседания в пределах построенной из точки Р области влияния, имеющей форму конуса с образующими, наклоненными под углом, равным граничному углу у, в которой горные породы теоретически должны быть свободны от деформаций сжатия или растяжения, образуется, по всей вероятности, только над серединой площади сверхполной подработки когда уплотнение закладки в очистной выработке достигает предельной вели­ чины. Существенные ошибки в результатах расчета оседаний с помощью инте­ грационных сеток могут быть вызваны неправильным выбором «ступеней осе­ дания» у контура очистной выработки [коэффициент а в формулах (10) и (47)], в особенности при малой расчетной высоте h, не превышающей 200 м.

3.3.2.

Массив горных пород как упругая среда

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

В несжимаемом упругом полупространстве (р, = 0,5), на поверхности кото­ рого приложены равномерно распределенная нагрузка и равномерное обжатие по площадке круга диаметром 2а (рис. 39), осесимметричные оседания на ради­ альном расстоянии г и вертикальном расстоянии h для равномерно распреде­ ленного давления [32] могут быть получены из выражения84

(48)

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

Соседние файлы в папке книги