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

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

со скоростью подвигания забоя v (м/сут) и радиусом площади полной подра­ ботки R (м) [407]:

 

dvzi

- (max) = 0,6

R - v

(364)

 

dt

 

 

dvхДИН

(max) = 0,25

R

v

(365}

 

dt

 

 

В СССР для условий Донецкого бассейна применяется формула

 

 

*ин (max) = 30

,

мм/сут,

(366)

где

v — подвигание

очистных работ, м/мес

[158].

 

Применяемые в ПНР формулы для расчета максимальных деформаций

земной

поверхности имеют следующий вид [202, 345]:

 

для

наклонов

 

 

 

 

 

^2 max

0,7-М ,0

аRМ

,

мм/м;

(367)

для

радиуса кривизны

 

 

 

 

 

 

Д2

 

(368)

 

Prmin =

650 -f- 720 а М

м;

для

деформаций сжатия

 

 

 

—esmax=— “TjT-

ИЛИ

О.бУг max, мм/м;

(369)

для деформаций растяжения

 

 

e s щ а х =

0,45

,

мм/м.

(370)

В формулах (367)—(370) вынимаемая мощность указывается в миллимет­ рах. Получаемые по этим формулам деформации (за исключением наклонов) больше, чем наблюдающиеся в Рурской области (табл. 20), однако еще большие деформации наблюдаются в Саарском бассейне1 [75, 234].

1 Первые формулы для приближенной оценки общего горизонтального сдвижения А«общ, по профилю мульды оседания были предложены в 1919 г. К. Леманом; они были выведены, исходя из предположения о поворотах породных блоков навстречу друг другу с обеих сторон выемочного поля — половина длины очистной выработки), причем оседание кровли в се­ редине очистной выработки принималось равным аМ. Эти формулы имели вид [227]:

оН

А$общ 1—2яД/ ^ ;

“Ь А$общ — 0,5 (—А^общ) •

Т А Б Л И Ц А

20

 

 

Условия

Формулы,

Определяемый

выведенные

подработки

аналитически

параметр

(полная или

для условий

сдвижения

неполная)

Рурской

 

 

области ФРГ

Оседание, мм

Полная

а М

 

Неполная

а М е

Наклон, мм/м

Полная

.

-

аМ

 

 

1 2

 

R

 

Неполная

0,7-

аМ

 

R

Сжатие, мм/м

Полная

 

 

аМ

0-25- Н

 

Неполная

0,5-

аМ

 

Н

Растяжение,

Полная

0,25-

аМ

Н

мм/м

 

 

 

 

 

Неполная

0,25-

аМ

 

Н

Минимальный

Полная

500-

R 2

а М

радиус кри­

 

 

 

визны, м

Неполная

 

 

/?2

 

250- а М

Формулы,

 

полученные

Место, для которого

эмпирически

для условий

определяются параметры

Саарского бас­

 

сейна ФРГ

 

аМ Середина мульды над цент­

аМ е ром очистной выработки

а А/ 1,9- Над контуром очистной вы­

I работки

Середина мульды над цент­ ром очистной выработки

1,3-

а А/

 

Н

 

0,65 -

а А/

Над серединой зоны опор­

Н

 

 

ного давления

12

250

а М

Над серединой очистной

 

выработки

(кривизна вогну­

 

 

тости) или над серединой зон ел

 

 

опорного

давления [(кривизна

 

 

выпуклости)

П р и м е ч а н и е : а — коэффициент оседания (0,5 при полной закладке и 0,9 при выемке с обру­ шением кровли); М —вынимаемая мощность (здесь в мм); е — коэффициент влияния, определяемый

при помощи интеграционной сетки; Н —глубина

разработки, м; я — радиус площади полной

попра-

ботки, м (Я = 0,7Н); г — расстояние по горизонтали от края мульды сдвижения до ее середины,

м (для

неполной подработки площадь выработки принята

равной 50% площади полной подработки).

 

В Голландии для оценки ожидаемых деформаций сжатия и растяжения

пользуются

приближенными формулами:

 

e smax l f 6 yz max^

\ мм/м;

(371)

max =

0»8l?z maxZ-1,

MM/M,

 

где I — длина области деформаций сжатия или растяжения (м) при выемке площади неполной подработки. Границей между зонами деформации сжатия и растяжения является линия, па которой оседание равно 50% максимального 1468].

 

В Остравском бассейне ЧССР для глубины разработки от 500 до 800 м

получены

следующие максимальные деформации:

)

для

наклонов

 

 

 

 

 

 

max =

4 , 4 ^

m ax,

мм/м;

(372)

для

деформаций

растяжения

 

 

% max = 0,'45l>zmax. мм/м;

(373)

для

деформации

сжатия

 

 

— ■e s n ia x = 0 , 6 3 y ; max, мм/м.

(374)

В этих формулах максимальное оседание выражается в метрах

[480].

 

Аналогичные формулы, применяемые в Великобритании, имеют следу­

ющий вид:

 

 

 

 

 

 

ДЛЯ

наклонов

 

 

 

 

 

 

-- Q/

vz max .

 

(375)

 

Uzщах —^

 

 

j.[

,

 

 

для

растяжения

 

 

 

 

 

max

о 8

 

max

*

 

(376)

Для сжатия

 

 

 

 

 

 

 

fcs max “

О О

vz max

»

(377)

 

 

 

//

 

 

где максимальное оседание vZmax, которое в зависимости от величины выемоч­ ного поля меньше, чем полное оседание аМ , должно переводиться в милли­ метры. Деформации сжатия достигают критических значений при ширине выемочного поля 0,24Н и его длине 2R. Коэффициенты в первых двУх Фор­ мулах приведены для выемочного поля длиной 2R и шириной 0,5# [276].

Если принять в формулах для кривизны и наклонов R = 0,7Н и выра­ зить коэффициенты округленными числами, то можно построить наглядную диаграмму, по которой могут быть ориентировочно определены экстремальные значения параметров сдвижения земной поверхности для М = 1 м в зависимо­ сти от глубины разработки и способа выемки полезного ископаемого (рпс- 185). Для других значений вынимаемой мощности и коэффициентов полуяенные‘ результаты могут быть вычислены при условии, что их изменение проИсх°Дит по линейному закону, т. е., например, при М = 1,8 м полученные значения наклонов и горизонтальных деформаций нужно умножить, а значения радиуса кривизны — разделить на 1,8. Из этой диаграммы можно видеть, что с ростом глубины разработки деформации земной поверхности уменьшаются. Есл0 изве­ стны размеры области влияния на земной поверхности или радиус пЛ01ДаДи полной подработки R, то можно определить сдвижения для любых точек про­ филя по их расстоянию от точек, где сдвижения достигают экстремальных значений. Например, наклон в точке 1 на рис. 186 составляет vz

Рис. 185.

Зависимость экстремальных значений деформации земной поверхности над очистной выра­ боткой от глубины разработки, мощности пласта и системы разработки:

I — при выемке с закладкой выработанного пространства; II — при выемке с обрушением кровли

Рис. 186.

Схема расположения точек с экстремальными значениями деформаций земной поверхности относительно очистной выработки:

I, II, III, IV — линейное распределение соответственно растяжений, сжатий, кривизны, наклонов; V — кривая оседаний; VI и VII — растяжение соответственно при полной и неполной подработке

1 0 .6 .

Критические замечания по расчету сдвижений земной поверхности

Наука о сдвижении горных пород, опирающаяся на работы сравнительно небольшого круга исследователей-марШнейдеров, за несколько десятилетий XX в. сделала, несомненно, заметные успехи. Если не говорить об особо слож­ ных горнотехнических и горно-геологических условиях, создаваемых крутым падением пластов, наличием нарушенных зон и высокой концентрацией очист­ ных работ, то в настоящее время имеется возможность с помощью существу­ ющих методов расчета определять ожидаемые сдвижения земной поверхности с достаточной для производственных целей точностью (+10% при пологом

и±20% при крутом падении). Однако нельзя не признать, что еще ждут реше­ ния некоторые проблемы, которые или возникали вновь или же существовали

иранее, но лишь теперь стали актуальными в связи с возросшей степенью

подработанности породных массивов, ростом

глубины

разработок

до 1000 м

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

крупных

зданий и

сооружений

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

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

отказываться от решения этих специальных задач, необходимо внести ясность

в некоторые

вопросы, вытекающие

из следующих

критических

замечаний.

1.

До

настоящего времени

проведено

еще

слишком мало

сравнительных

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

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

2. Ни одним методом расчета, за исключением метода конечных элементов, до сих пор не поддавалось учесть фактическое местное распределение конвер­ генции в краевой зоне очистной выработки и явление вторичного оседания, вызванного влиянием ранее отработанной горной выработки, граничащей с рас­ четной выработкой (активизация сдвижения, вызванная влиянием старых горных работ). Правда, для методов расчета с помощью интеграционных сеток было предложено компенсировать влияние недостаточной или избыточной кон­ вергенции в краевой зоне введением различных значений коэффициента оседа­ ния для двух краевых полос (см. рис. 72) или добавлением к площади очистной выработки полосы шириной от 5 до 15 м [24, 87 J, однако из-за сложности практического осуществления этих рекомендаций их чаще всего не выполняют, -хотя локальное распределение величин конвергенции в выемочном поле как исходных данных, от которых зависит характер процесса сдвижения пород­ ного массива, должно вводиться в расчет с возможно большей точностью. Кроме того, часто не принимается в расчет влияние околоштрековых бутовых полос на процесс оседания (см. рис. И).

3.Чтобы иметь возможность с достаточной точностью определять ход раз­ вития процесса сдвижения во времени, необходимо вводить в расчет данные on изменении с течением времени таких факторов, как конвергенции в очистной выработке (см. рис. 171) или деформирования целиков (см. рис. 14), степени разрыхления и последующего уплотнения пород основной кровли в зоне обру­ шения (см. рис. 7), об уменьшении с ростом степени подработанности массива изгибной жесткости породных слоев, а также данные об изменении во времени деформированного состояния пород в средней и верхней части горного мас­ сива (см. рис. 173). Все это вероятно, можпо осуществить только теоретиче­ скими методами с многократным повторением расчета для различных сочетаний влияющих факторов.

4.Наблюдающаяся даже при горизонтальном залегании пласта асиммет­ рия профиля мульды оседания (см. рис. 165), явно связанная с направлением развития очистных работ, а также зависящая от степени подработанности и ориентировки систем трещиноватости породного массива [324] может быть учтена при расчете путем повышения «веса» ячеек интеграционной сетки в той

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

5.Дополнительное оседание породных слоев кровли, обусловленное возобновлением процесса сдвижения под влиянием соседних с очистной выра­ боткой старых (ранее отработанных) выработок, приводит к тому, что оседание земной поверхности над очистной выработкой оказывается больше, чем полу-