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

книги из ГПНТБ / Скуба, В. Н. Исследование устойчивости горных выработок в условиях многолетней мерзлоты

.pdf
Скачиваний:
4
Добавлен:
20.10.2023
Размер:
7.98 Mб
Скачать

и условии ограниченности на бесконечности lim Т2(t, R) = Ге Д-*-оо

принять координату протаивания i?, равной нулю, и прирав­ нять результат температуре плавления Гпл. Тогда получаем

Tm = tB — {t„ Те) ехр (7г2х2г') erfc (А К ^ т 7).

(Ш.35)

Отсюда

ехр (Д2к2т') erfc (А J/"х2т') = !В_ .Г ПЛ< 1-

(III.36)

1 е

 

В момент т = т ' задача становится двухфазной. Для ее ре­ шения воспользуемся квазистационарным методом (Лыков, 1967), применение которого в данном случае допустимо, так как решение уравнения теплопроводности при граничном усло­ вии первого рода является асимптотой для решения при граничном условии третьего рода. Для наших расчетов потоком из мерзлой зоны можно пренебречь, тем более что до­ пускаемая при этом ошибка идет в расчетный запас. В соответ­ ствии с используемым методом распределение температуры в талой зоне представим в виде

T1= c1B-j-C2.

(III.37)

Для определения произвольных постоянных сг и С2 имеем

(111.34) и условие

 

Г1= Г ПЛ

при R = R T (т).

(III.38)

Выполняя необходимые преобразования, находим

 

 

1

 

 

Тг

1 +АД: -[(ЛЯА

Т-ал) h (tB ^пд)^].

(III. 39)

Уравнение движения границы оттаивания определяем из условия Стефана, которое для данного случая имеет вид

■,V1 ш \ R=RT “

<IR„

~ЗГ-

Подставляя (III. 39) в (III. 40), получаем

d*r

М р в -Упл)

1

dr

LnV*wn

* 1 + hRT

При начальном условии

(III.40)

(III.41)

Дт (^l)

(III. 42)

42

интегрирование (III.41) дает

2Xifc2 (t„ — T ^ )

(ш .43)

(i + м гт)2 = i + — z^ _ j e L (т _ Т]).

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

1 -f 2fтд —тJ иw=(1 + hR д)2,

где

к =

с* (*в — Тпл)

М ^ в -Г д л )

^2 '

«2ЬлР2^л

Если р корень уравнения

ехр (р2) erf с (р) = -?— ^ 5 =

т

 

i- i® -

то из (III. 36) имеем

(111.44)

(111.45)

(111.46)

Т' ___ L

(111.47)

~

ftW

 

Функция exp(P2)erfc(p)

затабулирована

(Карслоу, Егер,

 

 

Т

 

 

1964). Зависимость Р о т ---- -- показана на рис. 15.

Подставляя (III.47) в (III.44), после несложных преобра­

зований

получаем

 

 

 

кг (2/стди 2— R l) -

2R Rh - 2/ф 2 = 0 .

 

 

 

 

(III.48)

 

 

Отсюда,

полагая

 

 

 

 

F„ = ^ .

(III.49)

Рис. 15. Зависимость р от

 

 

 

находим

 

 

Т

Т

 

 

- t ± -

при _Л > 4,6; р =

 

 

 

Т

h = —

f + V l +

2йрг (2kF0 - i )

 

 

1 -1 1 -

■Кд

2kF„ — l

 

 

 

 

(III.50)

 

я

43

В общем случае

 

Л =

(III.51)

Для примера рассчитаем толщину однослойной изоляции из пенополиуретана, не допускающей протаивания более чем на 3 м за 105 ч (10 лет), полагая температуру стенки выработки равной температуре воздушного потока. Для этого случая

^i = l,l ккал/м-ч- град; а= оо ; Я*= 0 ,0 3 ккал/м-ч*град; /?д= 3 м;

тд= 105ч;

Х2 = 2 , 5- 103м2/ч;

р=2200 кг/м3; 1/л= 80

ккал/кг;

tB = 2° С;

Те= —2°С.

 

 

 

Из уравнения (III.46)

определяем

р= 0,75. По

формуле

(III. 45)

и (III.49) находим &=Q,1 и

F0=21,8. Из

(III.50)

определяем <^=0,1631/м. Из (III. 51) находим

 

 

б = Ь -

- 4 - = 0,167

м.

 

 

Aj

tl

 

 

Выражение (III.51) при 2kFo— 1= 0 обращается в бесконеч­ ность, т. е. толщина изоляции равна нулю и граничное условие третьего рода переходит в условие первого рода. Это означает, что условие протаивания не больше чем на 7?д метров за тд часов будет выполнено и без изоляции.

Если в плоскости 7?д, тд построить кривую 2kFo—i или

ЬлР^л

(III.52)

m

2*1 (*в-Гщ.)

 

то для точек, лежащих выше этой кривой, изоляция требуется, а для точек, лежащих ниже, она не нужна. Для данных нашего

примера уравнение (II 1.52) перейдет в тд = 2* Ю3# |.

Время эффективной работы изоляции тэф (т. е. такой период, за который на контакте изоляции с породой температура при­ мет значение, близкое к граничному, например, Тг(тэф, 0)=61^в,

Т

\

 

подставив

(II 1.43) в

где - у ^ - < 0 < 1 ] можно определить,

(III.39)

с учетом (III.47) и R = 0. Тогда

 

 

 

t

т

 

(III.53)

 

Тг (г, 0) = t]

 

 

 

V 1 + 2А (™2Л2 — Р2)

 

Приравняв правую часть

(II 1.53)

при т = т вф

величине

44

®tB, получим для определения тЭф следующее выражение:

(III.54)

Формула (II1.54) позволяет определить эффективное время работы изоляции, толщина которой определяется выражением

(III.50).

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

Неравенство служит условием для определения времени, при котором решения осесимметричной и плоскопараллельной задач отличаются незначительно. Если неравенство (III.55) не выполняется, то задачу следует решить в осесимметричной постановке. Как и для плоскопараллельного случая, найдем время тд, за которое на границе Л =Ло температура станет рав­ ной Гпл. При этом, поскольку т'^Тд, время т7 можно опре­ делить из выражения (III.36), полученного для плоскопарал­ лельного случая. Такое допущение ведет к некоторому увели­ чению толщины изоляции, что пойдет в расчетный запас.

Для времени т ^ т ' в талой зоне

(т)т зададим ста­

ционарное распределение температуры

 

А хIn R -j-H.2.

(III.56)

Постоянные А х и А 2 найдем из граничных условий

 

(111.57)

TJ(Rt) = Тпл.

(111.58)

После этого формула (III.56) примет вид

 

 

(III.59)

Из условия Стефана

(III.60)

45

получим следующее уравнение для нахождения закона движе­ ния границы протаивания:

-

dx = Дтih 1пж +i t ) dR-

(ni-6i)

Введя безразмерные

величины

 

и

, х 2т

(III.62)

R ? = д2-; ^0 = —£; h = hR0,

0

^0

 

перепишем (III.61) в виде

 

khdF0' = (Д1пЛт-[ 1) R TdRT,

(И 1.63)

где к определяется по формуле (III.45).

 

Интегрируя (III.63) при начальном условии

 

где

Дт (^о) = 1,

(III.64)

 

 

получаем

2 (/’о — F l)h k = R l{h \n R ? 4- l) — 1.

(III.65)

Как и ранее, условием для выбора изоляции является протаивание не более чем на 7?д метров за тд часов. Безраз­ мерными координатами для данных величин служат

-

ЗД-Д

 

х2тд

R* —

a r - ’

F

(111.66)

 

Заметим, что в соответствии

с

(II 1.47)

 

 

 

 

(III.67)

Подставляя выражение

(III.66)

и (III.67) в (III.65), на­

ходим

 

 

 

 

К - 1) + К К - ‘ )2 + 4*Р2 [ 4 + ( д | - 1) - R \ In R \

h =

 

i f f - l

-Л"1пДЙ

4№0+

(III.68)

46

 

Результаты расчета

па­

 

 

 

раметров теплоизоляции для

 

 

 

величины тд, равной 3; 10;

 

 

 

20

годам,

приведены

на

 

 

 

рис. 16 {1—3).

 

 

 

 

 

 

Вообще говоря, для осе­

 

 

 

симметричного случая выра­

 

 

 

жение для h отличается от

 

 

 

выражения (III.51) для плос­

 

 

 

копараллельного

случая.

 

 

 

Но

ввиду

малой

величины

 

 

 

отношения

26,

 

 

 

 

 

-гг— выражение

 

 

 

(II 1.51)

может быть исполь­

 

 

 

зовано и в

осесимметричном

 

 

 

случае.

 

 

 

 

Рис.

16. График для

определения

 

Рассчитаем толщину

од­

 

параметров изоляции

горных выра­

нослойной изоляции для вы­

боток в условиях многолетней мерз­

работки радиуса К0= 2 м при

лоты при £в=2,0°, Тпл= 0°, Те= —Ъ°.

условиях предыдущего при­

 

 

 

мера.

уравнения

(III.46) определяем р=0,75. По формулам

 

Из

(II 1.45)

и (III.66) находим &=0,1; # т=2,5; Fo=62,5. Из (111.68}

определяем /г=0,578. Из (III.51) и (III.68) находим

 

 

 

 

6 =

£ -

• Д?- =

0,094 м.

 

 

 

 

 

 

 

h

 

 

 

Сравнивая, видим, что при тех же условиях в плоскопарал­

 

 

 

 

лельном

случае

необходимая толщина изо­

 

 

 

 

ляции

в

1,77 раза больше, чем при рас­

 

 

 

 

смотрении цилиндрической выработки. Как

 

 

 

 

и для плоскопараллельного случая, в плос­

 

 

 

 

кости Кт,

F q м о ж н о построить кривую

 

 

 

 

^0

 

 

(III.69)

Рис. 17. Область применения изо­ ляции при Те=

= - 2 ° . Г„л = 0°,

разделяющую плоскость на две области

значений — R и Fq (рис. 17). В области* лежащей выше кривой, изоляция нужна, а в области, лежащей ниже, изоляция не требуется.

47

* * *

Сформулируем выводы по главе III.

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

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

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

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

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

шают необходимую толщину изоляции.

Г Л А В А I V

ТЕОРЕТИЧЕСКАЯ ОЦЕНКА УСТОЙЧИВОСТИ ГОРНЫХ ВЫРАБОТОК

ПРИ ПОЛОЖИТЕЛЬНОМ ТЕПЛОВОМ РЕЖИМЕ ШАХТ

§1. Влияние оттаивания пород вокруг горных выработок

на их устойчивость при слоистой трещиноватой кровле

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

48

большей степени, чем больше их льдистость. Это объясняется распадением значительных сил сцепления породных слоев и отдельностей по ледяным прослойкам, длительное сопротивле­ ние разрыву по которым составляет от 1,6 кг/см2 при темпера­ туре —2° до 3,6 кг/см2 при температуре —10°.

Исчезновение такого мощного «крепления» при оттаивании пород, безусловно, сказывается практически во всех случаях эксплуатации выработок на угольных шахтах. Об этом свиде­ тельствуют наблюдения за поведением незакрепленной и за­ крепленной кровли в горных выработках шахт, разрабаты­ вающих пласты угля в различных температурных зонах. Однако влияние положительной температуры на различные типы по­ род неодинаково. Глинистые породы, например, превращаются при оттаивании в несвязанную текущую массу, которая за­ полняет пустоты за крепью, гидростатически пригружает крепь и проникает в выработки через неплотности и щели в ней. Значительно в меньшей степени влияет положительная темпе­ ратура на коренные породы с известковым и кремнистым це­ ментом, имеющие относительно невысокую влажность. В дан­ ном случае определяющий фактор устойчивости обнажений — трещиноватость пород.

Интенсивно трещиноватые осадочные породы при оттаива­ нии превращаются в сыпучую массу. При крупноблочной сло­ истой кровле из-за наличия открытых (из-за вытаивания льда) трещин, поперечные размеры которых могут быть весьма зна­ чительными, а суммарная ширина в пределах пролета выра­ ботки такова, что породы имеют возможность опускаться на определенных участках, не испытывая при этом значительных сжимающих напряжений, схема работы породного слоя в прин­ ципе совпадает со схемами Г. Н. Кузнецова (1959, 1961) и А. А. Борисова (1948, 1962) и может быть представлена в виде трехшарнирной арки.

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

пролета выработки,

объемного

веса пород,

мощности

слоя,

величины смещения кровли и характера

прило­

жения

нагрузки.

 

 

 

Определение площади среднего шарнира, от которой в ко­

нечном

счете зависит

величина

горизонтального

распора,

по А. А. Борисову, производится эмпирическим путем. Г. Н. Куз­ нецов находит высоту области соприкосновения блоков исходя из их геометрии, что соответствует площади шарнира на еди­ ницу длины выработки. У обоих авторов в качестве исходной схемы взят случай, когда блоки, слагающие арку, плотно при­ жаты друг к другу и по мере опускания среднего шарнира об-

4 в. Н. Скуба

49

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

(рис. 18).

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

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

Рис. 18. Схема к расчету размера

Рис. 19. Схема для определе-

зон смятия пород в замке и пятах

ния стрелы прогиба трехшарнир-

трехшарнирной арки.

ной арки.

50

Угол поворота блока, определяющий период работы трех­ шарнирной арки, может меняться от 0° до некоторой величины, при которой проекция блока А С на горизонтальную ось МС'

будет равна

Тогда

ап = arcsin

L + k

— а.

(IV.2)

 

'■V^ 4

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

еесуществование невозможно.

Всвязи с тем, что блоки сминаются не только в замке, но и в пятах арки (точка А на рис. 18), величина области сопри­

косновения блоков в = EF) уменьшается вдвое при том же угле поворота блока. С учетом смятия блока в пятах арки ве­

личина области

соприкосновения

 

\ ^

h t + ‘- j sin (“ О+ а п) ~

- у - — ~2

сй

sin 2ап

(IV.3)

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

кСж= ~ .

(IV.4)

Выражение для определения величины области соприкосно­ вения блоков в замке с учетом их деформации при сжатии примет вид

] А | +

sin («0 +

«л) - -у- - \

Lo.сж

2Е

С6-

 

 

sin

2а^

(IV.5)

 

 

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

Данные о разрушающей нагрузке и области соприкосновения блоков в замке при определенно^ размере трещин приведены

4*

51

Соседние файлы в папке книги из ГПНТБ