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

книги / Сдвижение горных пород.-1

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

I

I

I

Масш т абы : _L расстояния ?!??£?* T

сдвиж ения

*fl

TtTtfi

,

* ft rtrt ffl

f

ft t*rn<

 

t

ft

T ift]

 

t

f t

f t

f t

tf t f t

//

Рис. 2.

Поле векторов смещения, вычислен­

ное

при oi= 0° (Подмосковный бассейн)

Масштабы;

расстояниямОЮ20904950*

*■■■■| |

cà&uce/w.г М О М Ц И Ю Ю Оаи

Рис. За. Поле векторов смещения, вычисленное при cL= 45° (Дон­ басс)

У

Рис. 36. Векторное поле прио(.= 15‘

V *

8оу

(9)

 

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

 

 

 

 

l - х р

( 10)

 

 

 

 

ОЖ

 

 

 

 

где % и t] — соответственно горизон­

Рис. 4.

Схема-

линейного

тальная

и вертикальная

компонен­

распределения

составляю­

ты смещения. Коэффициент К в этом

щей X

по вертикали

частном

случае является

расстояни­

ем от поверхности до «нейтрального слоя», в точках которого v —

= 0 и ? = 0.

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

II

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

15 +12=о.

(П )

0Х tty

 

Из уравнений(11) и (Ю)следует уравнение относительно вертикаль­ ной компоненты:

(12)

Уравнение (12) показывает, цго вертикальные перемещения

под­

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

Уравнение (10) подобно уравнению

92w

(13)

Z3 * 2’

 

которое имеет место в прогибающей упругой и пластической

балке. Деформация t x

пропорциональна кривизне

нейтраль­

ной оси и коэффициент

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

расстоянию

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

Коэффициент К представляет функцию от глубины', т.е. функ­ цию от у= Н. Экспериментальные данные показывают, для ус­ ловий горизонтального залегания в точках поверхности

К=0,2Н.

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

K = ty= 0,2y,

и, следовательно,

(14)

Т ч = А у $ 1

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

ривать как среду,

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

среду пластическую.

и г] весьма малые по сравнению с размера­

Если величины ^

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

указанной, дает

и упругое решение этой задачи.

 

Как известно, решение уравнения

(14) есть:

 

v = -

nv

S ± * \ + &

( l z *

(15)

Y

VA y /

 

у

 

где S — половина ширины выработки,

m — мощность выработки,

А= const « 0,2.

[3]

подробно развил

решение уравнения

(14),

Р.А. Муллер

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

ми о % и tj достаточно хорошая.

Ill

Уравнение (14) для случая негоризонтального залегания (ot^O) должно несколько усложниться:

В нашей

работе [4]

было показано, что для

условий о£ ^ О

получается

.уравнение:

 

 

 

 

1* = К <»

JÜn

9v

(16)

 

 

 

} х г + B W H

 

Здесь коэффициент В(у) отражает влияние угла падения пласта. Нами было выяснено [4] его значение. При горизонтальном за­ легании геометрическое место центров,т. е. середин кривых проги­ бов, представляет вертикальную линию, проходящую через сере­ дину выработки. При ci # 0 такая линия должна была бы быть нормалью, проходящей через середину. Но сечения, параллель­ ные пласту, прогибаются по нормали и середины этих линий про­ гиба смещаются от нормали, проходящей через середину выра­ ботки. Коэффициент В(у) входит в уравнение кривой, представ­ ляющей геометрическое место середин кривых прогибов по норма­ ли для сечений, параллельных пласту (рис. 5, линия 00, ).

Решение уравнения (16) будет:

V

х - Ш +s,

д , »

- 0 { y ) - s o

14Гу

}

4 4

У £у

(17)

 

М 1

где-ц— величина смещений

по нормали

к пласту, ш — вынима­

емая мощность пласта, S — половина ширины выработки,

 

А=К(у) = const,

 

0(ÿ) = j B(y)dy,

 

 

 

х-в(у)+&

о

 

 

 

 

 

х - 6(y) + s

 

 

сИ

 

VF У

 

о

 

 

 

 

 

 

 

 

 

Рис. 5. Схема отклонения от нормали центров кривых изгибов

— табличный интеграл Гаусса. При малых углах падения (о^<

<30°) линию 00, можно, в первом приближении, принимать прямой, представляющей геометрическое

место точек, в которых h b = 0. Тогда, как показано в [4], дх

®(y) = x=ytg<ç .

( 18)

Таким образом,

 

Ъ(у) = tg<p ..

(19)

Из рис. 6 следует, что ф = Л — 90-f-oL Угол А — есть угол накло­

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

непосредственных

наблюдений

и приблизительно равен 90+

I

2

fi

п

*

> гДе 1* и Р — углы сдвижения, тоже достаточно изучен­

ные. Поэтому можно написать:

go + - 1 ^ 1 - 9 0 +оС=

-f 6L

Рис. 6. Схема расчета кривых смещений

Рис. 8. Кривые прогиба по нормали, вычисленные по формуле (17)

На рис./ для сравнения приведены рассчитанная и фактичес­ кая кривые оседания поверхности, а на рис.8 — кривые прогибов по нормали, рассчитанные по формуле (17).

Для пространственного случая при ol= Ü имеем [4J :

+ Ь )

Решение этого уравнения, как это выполняется в курсах мате­ матической физики [5], применительно к нашей задаче будет:

(a-s)z4-(v-z)2

1

2Х(У)У

 

И110 2яК{у)у е

&Ad,v.

(21)

Для условий прямоугольной выработки размерами

2S1 и 2S2

при начале координат в центре выработки будет:

 

_ (7i-x)%(\>-z)z

 

 

e

«K(y>y

cl À d v.

(22)

 

 

Интегрируя, получим в символах Гаусса:

где К (ч) = Ау.

Совершенно другим путем профессор Е.Литвинишин получил для сдвигающейся толщи уравнение параболического типа, но бо­ лее полное:

(24)

Профессор Литвинишин и его ученики весьма подробно ис­ следовали решение, уравнения (24) применительно к процессу сдвижения пород.

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

I V

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

За функцию влияния принимается экспонента в форме ана­ логичной гауссовской экспоненте, как это имеем в формуле (22). Если функция влияния выбрана достаточно хорошо, то расчет сдвижений поверхности сводится к суммированию (интегриро­ ванию) этого влияния по всей площади выработки, как это де­ лает формула (22).

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

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

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

V

 

Можно утверждать, что параболические уравнения

(16) и

(24) отражают главные свойства механизма сдвижения

пород.

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

влияния в форме

аналогичной

экспоненте

Гаусса, дают

наибо­

лее практически удовлетворительные результаты.

 

 

 

 

 

ЛИТЕРАТУРА

 

 

 

 

1.

С. Г.А в ер ш ин . , Горные работы

под сооружениями

и

водоема­

ми. М., Углетехиздат,

1954.

 

 

 

 

 

 

 

2.

В.В.С о к о л о в с к и й.

Статика

сыпучей

среды. Технико-теорети­

ческое издательство. М., 1954.

 

 

 

 

 

 

 

3.

Р.А.Муллер. Влияние

горных

выработок

на

деформацию

земной

поверхности. М., Углетехиздат,

1958.

 

 

 

 

 

 

4.

С.Г.Авершин.

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

влиянием подземных разработок. Л., ВНИМИ, 1960.

 

 

 

5.

С.Л.С о б о л е в.

Уравнения математической

физики.

М

Гостех-

издат,

1954.

 

 

 

 

 

 

 

 

п докт°Ра МаРтоша на VI совещании IV группы М еждуна­ родного бюро по механике горных пород в г. Кракове в 1969 г.

ТРУДЫ ВСЕСОЮЗНОГО НАУЧНО-ИССЛЕДОВАТЕЛЬСКОГО ИНСТИТУТА

ГОРНОЙ ГЕОМЕХАНИКИ И МАРКШЕЙДЕРСКОГО ДЕЛА (ВНИМИ)

Сб. 89

1973 г_

Канд.технмаук ИЛ.Петухов

ВЛИЯНИЕ НАНОСОВ НА СДВИЖЕНИЕ ГОРНЫХ ПОРОД ПО НАПЛАСТОВАНИЮ

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

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

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

положенного

в сторону

восстания

(ABCD,

 

рис.1),

определяется

условиями специального пре- -

 

 

 

 

 

 

дельного

равновесия

[2,3]

 

À

У

*

t a

и может быть записана урав­

 

 

 

 

 

 

 

нением

= Ntgp'+ K'L-f

н

 

V

\

/

SOT

! Çtf

Рт

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 2Рбз+Рбо+Рн.

О)

 

 

 

 

 

 

где Ntgp1—трение по кон­

/

х 7п

 

 

 

 

тактам

слоев;

K L — сцеп-

 

 

 

 

ление

Tid

контактам

слоев;

 

т

 

 

 

 

2P(j3—боковое

защемление

 

 

 

 

 

 

призмы

 

сползания;

Pg0—

 

 

рис. i

 

 

боковой

 

отпор,

возникаю-

 

 

 

 

щий при прогибе слоев.в об­

отпор прогибающихся

наносов.

ласти сдвижения;

Рн — боковой

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

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

пользуя гипотезу проф. Г.И.Кузнецова о блокообразовании при

сдвижении горных пород

над

выработанным

пространством

[4] :

 

 

Т=

 

1И±\

 

 

(2).

 

 

~

2Ь/’

 

 

 

 

 

 

 

 

где

6с;нс— временное

сопротивление сжатию

горной породы

од­

ного

слоя; h — мощность

слоя; к — коэффициент,

учитывающий

отличие фактического

распределения напряжений

от принято­

го; Г}0— величина прогиба слоя.

при

наличии

пер­

Из выражения (2)

можно

установить, что

воначального распора, что всегда имеет место в массиве горных пород, боковой распор исчезает при прогибах, равных 0,2—0,5 мощности прогибающегося слоя.

Анализ геологоразведочных данных показал, что массив гор­ ных пород большинства угольных месторождений на 70—85% состоит из слоев мощностью менее 2 ж, а мощность разрабаты­ ваемых пластов более 1 м. Следовательно, при вертикальных сдвижениях более 0,2—0,5 м боковой отпор прогибающихся сло­ ев будет незначительным.

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

Для определения силы бокового распора прогибающихся

наносов можно принять ту же

гипотезу проф. Г.И.Кузнецова

[5], тогда

 

 

Т„ =

СОК 1т-н

(3)

 

 

Силы бокового распора наносов могут оказывать сопротив­ ление сдвигу призмы до тех пор, пока сцепление (Сн) и трение

(Рн) по контакту коренных пород с наносами будут больше сил бокового распора наносов

 

« W b P „ T.

(4)

V н вс*§?*;

P i= hHft, tgpH,

(5)

где Селе — сопротивление сжатию наносов; hH— мощность

нано­

сов; Нв — глубина верхней границы горных работ; уи — объемный , вес наносов; рн — угол трения по контакту коренных пород

снаносами.

Втех случаях, когда неравенство (4) не будет соблюдаться, коренные породы сдвигающейся призмы ABCD (рис. i) будут

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