Файл: Павлов, А. В. Искусственное оттаивание мерзлых пород теплом солнечной радиации при разработке россыпей.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 19.10.2024

Просмотров: 67

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

s (t ) = и т ,

где х — коэффициент пропорциональности.

Подставим значение s (t )

в условие Стефана:

 

2 1 +

х х

 

Решая это уравнение,

получим

 

/

"

KQw

т =

i,t -I КЧу,

*4,

 

"kitn 5 - 1

 

 

(IV.28)

(IV.29)

В критериальной форме выражение (IV.29) запишем в виде

где

 

Кг = ек ' - К х -

1,

(IV.30)

 

3 »

 

 

*2g№T

 

 

К,

к.

 

 

 

 

 

 

 

k2tn

 

Зависимость (IV.30) можно с большой точностью аппрокси­

мировать степенной функцией

 

 

 

 

 

 

к 1 = ак$п,

(iv.3i)

где Й=1,18,

п = 2,6.

 

 

его

значение, определим %

Подставляя вместо критерия К \

с учетом линейного изменения во времени термоизоляции на

по­

верхности:

 

QJ$nMn

 

 

 

 

 

 

(IV.32)

 

 

Kilo

 

 

 

 

 

 

 

 

Найдем отношение глубин промерзания под снегом с посто-

янным во

времени

термическим

 

/-

XX \

сопротивлением I s = —

и изменяющимся линейно во

времени

(s [t ] = k t ):

 

J,

 

Q X4W

V

•s +

,а 2^**дт

 

 

где

К ™ ( — 0 ,5 - / о , 2 5 + К1

m — = 0,62.

133


Для наблюдаемых в природе величин К 2 отношение изменяется от 1,17 до 1,22. В дальнейших расчетах примем его равным 1,2 и получим простую формулу для определения глубины сезонного промерзания горной породы под снегом при средней годовой температуре to, близкой к 0°

6 = 1,2 ( - S +

+

(JLV.33)

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

6 = 1 , 2 ^ - в (IV.34)

Для вычисления 8 в формулах (IV.25), (IV.34) необходимо знать температуру поверхности пород tu под слоем изоляции. Методика ее определения дана в следующем разделе.

Формулы (IV.25) и (IV.34) пригодны для расчета лишь максимальной величины £. Чтобы найти изменение g во вре­ мени, при интегрировании условия Стефана необходимо ввести переменные краевые условия.

Изменение температуры в зоне 3 с глубиной учитывается

тем, что величина q(\, т)

в соответствии с экспериментальными

данными считается пропорциональной не т,

а т-)-тн:

0(6.*) = ^3^0

1 + 0,6

1

(IV.35)

 

К т + тн ’

где параметр тн при всех т начиная с первого дня изменяется (по данным наблюдений в районе Загорска) от 6 до 10 дней (среднее значение тн принимается равным 8 дням).

Формула (IV.35) дает хорошие результаты при расчетах <?(6. *), совпадающие с данными наблюдений в природных условиях (см. рис. 27). Подставляя это значение q(\, т) в урав­ нение Стефана и интегрируя его с учетом произвольного изме­ нения во времени температуры поверхности tu (или темпера­ туры воздуха tS3), получаем окончательные формулы для опре­ деления | (температура замерзания породы принята равной ta,

а не 0°, как в формулах (IV.25), (IV.34).

Глубина промерзания: для породы под снегом

134


1 = 1

s2 + ___ 9__________

e У а т( т +

t h)

 

4w

 

 

 

(1 = 1,2),

 

(IV.36)

для оголенной породы

 

 

5 =

2Ан { ( ^ - д *

 

(IV.37)

------ 5— -------------------- 8 У а т ( т +

Тн).

Глубина

протаивания

 

 

еУоы(т + тн). (IV.38)

Н. И. Вотякова (1972) провела оценку точности формулы (1У.27) по сравнению с классическим решением Стефана, ко­ торое получено в виде трансцендентного уравнения и считается точным при постоянных краевых условиях. В решении Сте­ фана из всех факторов внешнего теплообмена входит только тем­ пература поверхности tu. Поэтому в качестве верхнего гранич­ ного условия при проверке и была задана эта величина. Таким образом, проводилась сравнительная оценка точности учета теплооттока в мерзлую толщу при протаивании пород, а также приближенной аппроксимации температурного поля в зоне протаивания. Тем же путем оценивалась точность применяе­ мых в практике формул Бергрена, Хакимова, Докучаева и Балобаева. Приводим некоторые результаты расчетов (табл. 49). Интервал значений средней годовой температуры пород был принят равным минус 1,2 — минус 12,4°.

 

 

 

 

 

 

Т а б л и ц а 49

Относительная

погрешность приближенных

формул

расчета глубины

 

 

протаивания пород, %

 

 

Автор

 

 

t o. °с

 

 

 

Время

 

 

 

 

 

 

— 1.2

-2,5

-5,0

- 7 ,4

- 9 ,9

-12 ,4

расчета

 

 

 

Бергрен

6,8

4,1

2,3

3,1

3,4

4,1

0,22

Хакимов

- 1 ,2

- 2 ,6

0,5

2,3

3,4

7,6

0,44

Докучаев

1,5

0,3

3,8

5,4

7,8

11,0

0,67

Павлов

1,4

0,5

0,8

0,8

1,2

4,1

0,30

Балобаев

1,5

—1,3

—0,6

- 2 ,1

—3,3

- 5 ,0

0,36

П р и м е ч а н и е . Время, необходимое для вычисления глубины протаивания пород из решения Стефана, принято за единицу.

135


Расчеты проводились при следующих исходных данных: порода супесчаная; объемный вес скелета 1500 кг/м3; весовая влажность и>с=20% ; период протаивапия т=3696 ч; средняя лет­ няя температура поверхности ta= 12,4°; объемная теплоемкость талой породы Ст=580, мерзлой См= 440 ккал/м3-град; коэф­

фициент

теплопроводности талой породы ?ц,=1,34, мерзлой

Я„=1,52

ккал/м-ч-град.

Выбранные формулы дают неодинаковое приближение к ре­ шению Стефана (его точность из трансцендентного уравнения, принятого за эталон, была оценена равной 0,2%). Наименьшую погрешность в учете теплооттока в мерзлую толщу дают фор­ мулы Павлова и Балобаева; у Бергрена она систематическая, но не превышает 7%. Несколько большая погрешность в форму­ лах Хакимова и Докучаева.

В действительности погрешности в расчете глубин промер- зания-протаивания могут намного превышать значения, при­ веденные в табл. 49. Это объясняется в основном следующими причинами:

а) исходные расчетные параметры (теплофизические харак­ теристики, влажность и пр.) назначаются, как правило, со зна­ чительной погрешностью, нередко превышающей 10% (Пав­ лов, 1973);

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

другого

поверхностного

покрова). Некоторые формулы (на­

пример,

Бергрена) не позволяют учесть этот важный фактор;

в) теплопередача при

промерзании-протаивании рассматри­

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

2. Промерзание пород под термоизоляционными покрытиями

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

Ц » , 2^ в з ( т - тз)

a J '

+ ^ w

е |/ а 3(т — т3)*.

*Все обозначения сохранены прежними.

136


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

нозах

глубин

промерзания процессы тепло- и

влагообмена

с атмосферой

можно

опускать, что

равносильно

допущению

R = L E = 0, а = о э . Тогда выражение для £ получит следующий

вид:

 

 

 

 

 

I =

-*из +п/ 4 +

— в! (Т~ Тз) -

е У а , ( т - т 3), (IV.39)

у

где

(IV.40)

лм

т8 — время запаздывания промерзания пород под слоем изо­ ляции, индекс «из» относится к изоляции.

Вспомогательный параметр е находится по табл. 48 в зави-

 

 

 

 

 

v

с *т°

и

V

 

^ м<п t *

симости от двух критериев: Лс = ------

л а =

—- —-

 

 

 

 

 

 

 

 

 

Ч-Ю

 

 

 

aj"w

 

 

температура породы под слоем изоляции).

 

 

 

 

 

Величину tn

определяем

приближенно, приравнивая зави-

симости

| =

1 /

/'2Х„<„т

И \

=

 

sH3 +

Г~с,

 

2 0 1 7 7

 

 

— —

] /

 

Н------ ----- :

 

 

 

у

 

"w

 

 

 

 

V

 

 

"w

 

 

 

 

еп-

-г1 =

1 + к п ~ у к п ( к п + 2),

 

 

(IV.41)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К .

 

* + •

 

 

 

 

 

( т > т ,) .

(IV.42)

 

 

 

 

 

 

 

 

 

Величина

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

по

табличным

данным

в

зави­

симости

от

значений критерия

Кп:

 

 

 

 

 

 

 

 

К

КП

 

 

 

 

§п

*п

 

ё„

*п

 

§п

 

0,00

1,00

0,10

0,64

0,5

 

0,38

1,4

0,22

4

 

0,10

 

0,01

0,87

0,15

0,58

 

0,6

 

0,35

1,6

 

0,20

5

 

0,08

 

0,02

0,82

0,20

0,54

0,7

 

0,32

1,8

 

0,18

6

 

0,07

 

0,03

0,78

0,25

0,50

0,8

 

0,30

2,0

 

0,17

8

 

0,06

 

0,04

0,75

0,30

0,47

0,9

 

0,29

2,5

 

0,15

10

 

0,05

 

0,05

0,73

0,35

0,44

1,0

 

0,27

3,0

 

0,13

15

 

0,03

 

0,07

0,69

0,40

0,42

 

1,2

 

0,24

3,5

 

0,12

20

 

0.00

 

Значение т3 определяется неодинаково в зависимости [от того, какие превращения происходят в слое изоляции.

137