Файл: Павлов, А. В. Искусственное оттаивание мерзлых пород теплом солнечной радиации при разработке россыпей.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