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

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

Рис. 170.

К выводу формулы для расчета временных коэффициентов z:

J — кривая оседаний на конечную стадию; 2 — при движущемся забое

наблюдений за развитием во времени сдвижения точки земной поверхности под которой проводится очистная выработка, вычисляя отношения измерен­ ных через определенные интервалы времени (например, ежеквартально) про­ межуточных значений оседания vZWH (см. кривую 2 на рис. 170) к конечному значению оседания у2К0Н, устанавливаемому впоследствии, когда очистная выработка достигнет размеров площади полной подработки [см. уравнение (311)]. Эти отношения

Zt = ^2 ДИН ^ 2 КОН

(Зо7 )

представляют собой значения факторов времени I, II и т. д. кварталов после начала воздействия подработки на данную точку земной поверхности. Так, например, для примера, приведенного в табл. 19, фактор времени для II квар­ тала составляет zu = 13 80 = 0,163. Таким образом, чтобы определить оседа­ ние на конец II квартала, достаточно будет вычисленное для данной точки зем­ ной поверхности конечное значение оседания vzкон = 80 см умножить на фактор времени zu = 0,16, т. е. vZfltHH= 80*0,16 = 13 см. При этом методе расчета предполагается, что очистной забой в процессе своего движения проходит рас­ стояние, равное диаметру всей площади полной подработки (см. рис. 170)у

Т А Б Л И Ц А 19

 

 

 

 

 

 

Вычисленное оседание в конечной

Временной

Участок отработки

Измеренное осе­

стадии процесса (z = l)

vz= aMet, см

 

 

коэффици­

(квартал)

дание v2 ДнН, см

 

 

квартальное

суммарное

ент 2

 

 

 

I

4

8

8

0,50

II

13

14

22

0,75

III

35,2

35

57

0,90

IV

54

15

72

1,00

V

69

8

80

1,00

Вторичное осе­

11

0

 

 

дание

и что скорость подвигания работ и глубина разработки соответствуют тем условиям наблюдений, для которых был вычислен временной коэффиИ#ент* Можно определять временной коэффициент и как отношение к наиболь­ шему возможному оседанию точки uz = aMet (кривая 1 на рис. 170) яри отра­ ботке площади, соответствующей положению забоя в момент времени t, т* е*

 

Н =

Vz дин

Уг,

 

 

так

что

для

примера, приведенного в табл. 19, будем иметь zn = 13

(8 ~

+ 14) =

0,6.

Таким

образом, при этом методе расчета сначала должно

быть

вычислено

конечное

оседание рассматриваемой точки для положения

забоя

на

конец

II

квартала vz = аМец = 1-0,22 = 22 см, произведение которого

на временной коэффициент zn = 0,6 даст оседание на конец II квартала, Рав" ное v2 = 22-0,6 = 13 см. Вычисленные по любой из приведенных формул

[(337) и (338)] временные коэффициенты не могут быть применены для условий, отличных от тех, при которых они определялись.

Оседания для промежуточных стадий процесса сдвижения могут быть вы­ числены по кривой оседания без учета замедления, т. е. при z = 1, если изве­ стна кривая замедления оседания, как показано на рис. 166 для забоя, ДР°" ходящего весь круг площади полной подработки. Найденные при помощи инте­ грационной сетки коэффициенты влияния для точки Р составляют eKfm = 0,537 при выходе забоя из круга площади полной подработки (vZK01l = 51,3 см) и et = = 0,4 или 74,5% екон при положении забоя, соответствующем моменту времени t, для которого производится расчет (t = 2 года). Если теперь из значения et-> выраженного в процентах, вычесть отсчитанную по диаграмме величину замед­ ления, равную 37,5%, то получим коэффициент влияния для промежуточной стадии процесса, равный 74,5—37,5 = 37%. Таким образом, оседание в момент

времени

t в данном случае составит у2ДИН =

51,3-0,37

= 19 см. Этот же резуль­

тат даст

расчет по формуле (338) при et =

vJaM (%),

если вместо того, чтобы

учитывать замедление оседания, определять временные коэффициенты Из соот­

ношения zt = у2Д„и et (здесь 37,5 74,5 = 0,5, смотри пунктирную линию на рис. 166):

дин== z ион ■—0,745 •0,50 •51,3 = 19 см.

Кривые для других скоростей подвигания очистных работ (найример, 0,5 м/сут, см. рис. 166) могут быть получены из кривых, построенных по дан­ ным наблюдений, если принять, что интервалы времени должны удлиняться пропорционально длительности воздействия подработки, причем вторичное оседание должно быть независимым от скорости подвигания забоя и одинако­ вым для всех кривых замедления. Преимущество этого метода расчета состоит в том, что можно обойтись без разбивки выемочного участка на части, отраба­ тываемые в определенные календарные сроки (а такой календарный план на стадии проектирования горных работ, как правило, неизвестен) и только впоследствии определить коэффициенты влияния, соответствующие выбранным участкам выемочного поля, имеющим произвольную величину [501*

10.3.2.

Временной коэффициент, отнесенный к единичной площади

Существует еще один подход к определению временного коэффициента, когда его относят не к постепенно возрастающей отрабатываемой площади, как эта было в предыдущем разделе [уравнения (337) и (338)], а к некоторой неизмен­ ной единичной площади, отрабатываемой за выбранный промежуток времени (например, за один квартал), и к степени влияния, соответствующей концу этого промежутка [87, 318]. Степень влияния выемки такой единичной (квар­ тальной) площади на точку земной поверхности в конце каждого интервала времени (I, II и т. д. кварталов) может выражаться, например, цифрами, при­ веденными в табл. 19* т. е. соответственно 50, 75, 90 и 100% полной величины коэффициента влияния. Для вновь добавляющихся единичных площадей / / , I I I и т. д. имеют силу, с учетом длительности их воздействия, такие же времен­

ные коэффициенты zt, как и для единичной площади I. Таким образом, в конце

III квартала оседание точки

Р (см. рис. 170), обусловленное влиянием I еди­

ничной площади, составит

0,9аМгь I I единичной площади — 0,75аМец и

I I I единичной площади — 0,5аМещ, так что общее оседание к концу III квар­ тала составит vZ3 = (90 X ej -f 0,75гц -f 0,5ет )аМ .

Закономерность образования формул для расчета сдвижений земной по­ верхности на промежуточных стадиях развития процесса видна из приводимых ниже уравнений, составленных для определения величин оседаний точки Р земной поверхности на конец I, II, III и т. д. кварталов после начала влияния подработки, если очистные работы ведутся непрерывно и на точку Р оказывают

влияние

и единичных площадей (z4 = 1):

 

 

vzi =

аМ (ziCj);

 

 

(

 

v22 = аМ (z2e1 +

z^ j);

 

 

 

VZ3йМ (Zs^j “Г

-f" 2LCJJJ),

 

 

(339)

Vz4

ClM (z4Cj “l" Z3£JJ “I" 22^III ~H sl^iv)’

 

 

 

vzb — ci M(z4ex zAen -f-

^oCjy-f~

 

 

v2 кон ~

(z4Cj -f- Z4CJJ-f- 24CJJJ -f- z4CjV -J- z4eу-p

~h z4eu)»

 

Процесс сдвижения земной поверхности в

данном случае

завершается

по прошествии и + (4—1) интервалов времени

(кварталов).

В матричном

выражении уравнения (339) принимают вид [212]

 

\vzAmx)= a M [Z ][e}.

 

 

(340)

Элементы этих или треугольных матриц располагаются по схеме уравнений (339), причем z4 = 1, и выемочное поле должно состоять из и единичных

' г д и н /

^1

иг2

vZ3

^24

^25

Vz^

Z2

Zl

 

 

z3

z2

zx

 

Z\

 

Z2

Zi

^4

ii»4

Z3

 

{ Z 4

- 4

s 4

Z 4

(341)

Число строк и столбцов здесь составляет и + (4—1). Коэффициенты влия­ ния et для отдельных квартальных выемочных площадей определяются обыч­ ным способом при помощи интеграционной сетки.

После рассмотрения этих основных соображений нужно еще показать, как этот общепринятый временной коэффициент zt, характеризующий измене­ ние степени воздействия площади очистной выработки на точку земной поверх­ ности через определенное число четвертьгодовых интервалов, прошедших после начала подработки, может быть выведен из хода развития процесса оседания этой точки земной поверхности и вычисленных коэффициентов влияния et. Обратимся снова к табл. 19, в которой приведены измеренные через кварталь­ ные интервалы величины оседаний точки Р, равные игДИН = v2l, и вычисленные при z = 1 составляющие конечного оседания vz для отдельных единичных

площадей (третий столбец табл, 19). Из уравнения (339), после некоторых преобразований, получим:

z, =

.

 

 

 

ехаМ

 

 

 

 

22 — (e\laM ) Z1 .

 

(342)

 

ехаМ

 

 

 

 

 

игя— (еш аМ) zi — (enaM) z2

 

 

егаМ

 

 

Подставив в эти выражения значения параметров, приведенные в табл. 19,

получим

временные

коэффициенты:

 

8 -— 0 5*

„ ^ *3 -1 4 -0-5

= 0J5.

 

-

g

 

^3 ”

35,2 — 35.0,5 — 14.0,75

— 0,9

и т. д.

 

8

 

 

 

В приведенном численном примере мы поступали так, как если бы отдель­ ные квартальные выемочные участки характеризовались одним и тем же коэф­ фициентом оседания а, определенным по средней конечной конвергенции во всем выемочном поле. В действительности же конвергенция в краевой зоне ши­ риной от 30 до 60 м оказывается значительно меньшей, чем средняя конвер­ генция по всему выемочному полю, как это можно видеть из графиков, при­ веденных на рис. 9—11. Поэтому в уравнение (342) должны подставляться

конечные

значения

конвергенции

at соответствующих выемочных участков

/, / / , I I I

и т. д.

(см. рис. 170),

если требуется определить временной коэф­

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

10.3.3.

Развитие во времени процесса оседания непосредственной кровли

На рис. 171 показан ход развития процесса оседания пород непосредственной кровли в зависимости от увеличения размеров очистной выработки. Аналогич­ ным образом происходит оседание основной кровли над выемочным полем, отрабатываемым с обрушением кровли. Можно видеть, что, например, оседание, н?!_2 на только что отработанном участке 2 заметно меньше оседания непосред­ ственной кровли wx_5 на участке 5 вскоре после отработки этой части пласта. Поэтому, чтобы правильно определить, с учетом фактора времени, степень влияния или конвергенцию et (t) отдельных участков в краевой зоне очистной выработки, нужно сопоставить опускание породного слоя непосредственной кровли на данный момент времени (t) не со средней или максимальной конвергенцией аМ всего выемочного поля, а только с конечным оседанием

а/М£соответствующего участка 7, 2, 3 и т. д. по

формуле

с,- (t)

u>i (0

%.

(343)

aiM