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

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

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

В условиях Донецкого бассейна зависимость опускания кровли 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 п. В отечественной терминологии это понятие отсутствует (примеч. переводчика).