Расчет системы термостабилизации склада тарного хранения масла

Рассмотрим выполнение прогнозного моделирования температурного поля грунтов основания склада тарного хранения масла. Для термостабилизации грунтов оснований применяется тепловой экран совместно с одиночными сезонно-действующими охлаждающими устройствами (термостабилизаторами).

Примечание: Пример подготовлен с использованием материалов ПАО «ВНИПИгаздобыча». Строительные решения, климатические параметры, геологические и геокриологические условия, последовательность строительных работ изменены для целей примера расчета и не соответствуют реальному объекту.

Склад. Вид общий

Климатические характеристики

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

ПоказательМесяцы
IIIIIIIVVVIVIIVIIIIXXXIXII
Температура в-ха, ⁰С-23.7-22.9-14.4-8.2 0.0 9.9 15.7 12.0 5.7-4.2-15.2-20.9
Скорость ветра, м/с3.03.13.33.53.83.83.33.03.23.33.02.9
Высота снега, м0.530.620.680.650.28----0.080.250.40

Геокриологические условия

Геокриологические условия площадки характеризуются сплошным развитием многолетнемерзлых грунтов, мерзлота сливающегося типа. Температура грунтов на глубине 10 метров составляет минус 1,0 ⁰С.

Склад. Геология
Свойства грунтов приведены в таблице ниже.
Наименование
показателя
ОбозначениеЕд.измеренияНасыпной грунтИГЭ №140200.
Суглинок тугопластичный
ИГЭ №130100.
Глина полутвердая
ИГЭ №141000Э.
Суглинок мерзлый
Физические свойства
Плотность грунта в сухом состоянииRdfкг/м31740147017101830
Суммарная влажностьWtotд.е.0,180,2880,1890,152
Степень засоленностиDsal%-0,0680,0710,08
Число пластичностиIpд.е.-0,1180,110,105
Влажность на границе раскатыванияWpд.е.-0,20,1960,188
Теплофизические свойства
Температура начала замерзанияTbf⁰С-0,15-0,34-0,4-0,47
Теплопроводность талого грунтаλthВт/(м*К)2,321,571,531,54
Теплопроводность мерзлого грунтаλfВт/(м*К)2,731,81,761,76
Объемная теплоемкость талого грунтаCthМДж/м32,633,12,892,8
Объемная теплоемкость мерзлого грунтаCfМДж/м31,932,322,352,37
Коэффициенты кривой незамерзшей воды
Ww(t)= A + B / (C - t),
где t - температура, ⁰С
А-0,00130,07360,07040,0631
B-0,00160,10260,11930,1818
C-0,44520,36190,60531,5792
Теплота фазового переходаQfМДж/м3----

Строительные решения

Грунты основания используются по I принципу (СП 25.13330.2012), т.е. с сохранением мерзлого состояния в процессе строительства и эксплуатации. Склад тарного хранения масла установлен на свайный фундамент с полами по грунту. В различных частях здания предусматривается поддержание различной температуры воздуха.

Склад. Температуры помещений

Температуры воздуха поддерживаются в холодный период времени. В теплый период температура приравнивается атмосферной. График температуры в помещении приведен на следующем рисунке.

Склад. График температуры

Коэффициент теплоотдачи внутри помещения принят согласно табл.4 СП 50.13330.2012 «Тепловая защита зданий», таблица значений условий теплообмена приведена ниже.

ПоказательМесяцы
IIIIIIIVVVIVIIVIIIIXXXIXII
Температура в-ха, ⁰С10,010,010,010,010,010,015,712,010,010,010,010,0
Теплоотдача, Вт/(м2*К)8,78,78,78,78,78,78,78,78,78,78,78,7

Сваи фундамента погружаются в грунт буроопускным способом. Бурится скважина диаметром большим чем погружаемая свая. Скважина заполняется цементно-песчаным раствором. Укладываемый раствор должен иметь положительную температуру (зимой подогреть до 20 ⁰С). Затрубное пространство сваи на высоту сезонно-талого слоя заполняется сухим непучинистым грунтом. Внутренняя полость сваи заполняется цементно-песчаным раствором.

Схема расположения свай склада тарного хранения масла
Температурная стабилизация грунтов оснований (ТСГ)

Технические решения по термостабилизации грунтов предусматривают совместное использование теплового экрана (теплоизоляция типа «Пеноплекс» толщиной 200 мм) с установкой одиночных наклонных сезонно-действующих охлаждающих устройств (термостабилизаторов).

Разрез задания с нанесением теплового экрана
Схема расположения сезоннодействующих охлаждающих устройств в основании склада

Прогноз температурного режима грунтов

Прогнозное моделирование температурного поля грунтов осуществлялось численными методами с учетом требований к порядку составления прогноза изменения температурного режима грунтов РСН 31-83 (п. 1.6., 2.35., 3.15.- 3.19., 4.16.-4.19.) и РСН 67-87. Расчеты выполнялись на программе «Борей 3D» (www.boreas3d.ru). Программа Борей 3D сертифицирована на предмет соответствия нормативных документов. Сертификат №RA.RU.АБ86.H01114.

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

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

Этапы выполнения строительства

Расчет системы термостабилизации грунтов выполнялся с поэтапным моделированием строительного периода:

  • Первый этап – возведение насыпи.
    Моделировалось изменение температурного режима грунтов при возведении насыпи. Принят наиболее критичный вариант для температур грунтов — возведение насыпи в конце летнего периода талым грунтом. Срок прогноза температурного поля грунтов составлял один календарный год. В качестве начального температурного поля приняты значения температур грунтов полученные при процедуре адаптации граничных условий.
  • Второй этап – период эксплуатации.
    Моделировалось тепловое воздействие от строительства и эксплуатации сооружения. Погружения свай и сезоннодействующих охлаждающих устройств (термостабилизаторов) принято на момент 15 мая. Температурное поле грунтов с первого этапа моделирования (возведения насыпи), принято в качестве начального температурного поля грунтов.
Адаптация условий теплообмена

Теплообмен на поверхности земли зависит от большого числа факторов: альбедо естественных или искусственных поверхностей, прямой и рассеянной солнечной радиации, конвективной (скорость ветра на поверхности земли), инфильтрации и испарения влаги, кондуктивной составляющей (теплопроводность и мощность растительного и снежного покровов) и т.д. Влияние некоторых из них на формирование теплового режима сложно корректно оценить. В связи с этим, расчеты производятся с использованием эффективных величин коэффициента теплопередачи. Методика расчета эффективной величины коэффициента теплопередачи (адаптации) реализована в ПО Борей 3D в автоматическом режиме. Методика адаптации будет приведена на сайте.

Значения адаптированных граничных условий приведены в таблице ниже.

ПоказательМесяцы
IIIIIIIVVVIVIIVIIIIXXXIXII
Температура в-ха, ⁰С-23,7-22,9-14,4-8,20,09,915,712,05,7-4,2-15,2-20,9
Скорость ветра, м/с3,03,13,33,53,83,83,33,03,23,33,02,9
Теплоотдача, Вт/(м2*К)9,539,7710,2510,7311,4511,4510,259,5310,0110,259,539,29
Высота снега, м0,2820,3300,3620,3460,149----0,0430,1330,213
Плотность снега, кг/м3141149155162377----98114128
Теплопроводность снега, Вт/(м*К)0,2030,2130,2210,2300,647----0,1570,1730,188
Коэффициент теплопередачи, Вт/(м2*К)0,6690,6070,5760,6273,15111,45010,2509,53010,0102,7091,1440,809
Расчетная область

Моделирование термостабилизации ведется в расчетной области, которая представляет собой параллелепипед с размерами 65 х 55 х 50 метров. Расчетная область грунтового массива с установленными сваями и термостабилизаторами показана на рисунке ниже.

3D модель для расчета термостабилизации вечномерзлых грунтов оснований склада тарного хранения масла

Схема расположения свай, сезоннодействующих охлаждающих устройств (СОУ, термостабилизаторов) и строительных конструкций в модели приведена на рисунке.

Схема расположения свай и вертикальных термостабилизаторов в модели для расчета температур многолетнемерзлых грунтов склада

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

В расчете учтено тепловое влияние от свай при устройстве фундамента, как с подогревом бетона при бетонировании свай, так и с теплом, выделяемым при гидратации. График тепловыделения от цементно-песчаного раствора М150 приведен на рисунке.

График тепловыделения раствора М150

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

ПоказательМесяцы
IIIIIIIVVVIVIIVIIIIXXXIXII
Температура в-ха, ⁰С-23,7-22,9-14,4-8,20,09,915,712,05,7-4,2-15,2-20,9
Скорость ветра, м/с1,81,92,02,12,32,32,01,82,02,01,81,8
Коэф. теплоотдачи к поверхности конденсатора, Вт/(м2*К)17,2917,5418,0318,52-----18,0317,2917,01
Эффективный коэф. теплопередачи, отнесенный к наружной поверхности испарителя, Вт/(м2*К)18,5018,7619,2819,81-----19,2818,5018,19

Область исследования разбивается на прямоугольные элементы произвольных размеров. Разбивка производится горизонтальными и вертикальными плоскостями (разбивочными плоскостями), параллельными соответствующим границам области. В расчетной области шаг горизонтальной сетки изменялся в пределах от 0,1 до 1,0 м. Вертикальная разбивка осуществлялась в пределах от 0.1 до 3.0 м. Расчетная область состоит из 276х184х97 ячеек (4,9 млн. ячеек).

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

Результаты расчета работы системы термостабилизации вечномерзлых грунтов приведены на рисунках ниже.

Трехмерное температурное поле грунтов на 15 ноября 2-го года эксплуатации.

Расчетное трехмерное температурное поле грунтов в основании склада

Сечение (разрез) температурного поля грунтов оснований.

На 15 октября 0-го года (начало эксплуатации).

Температурное поле грунтов оснований склада на начальный момент времени

На 15 января 1-го года.

Разрез температурного поля вечномерзлых грунтов оснований склада

На 15 октября 1-го года.

Температурное поле многолетнемерзлых грунтов оснований склада в результате работы системы термостабилизации грунтов

На 15 октября 2-го года.

Температурное поле многолетнемерзлых грунтов оснований склада в результате работы наклонных термостаблизаторов

На 15 октября 3-го года.

Расчет термостабилизации грунтов оснований склада

Расчетные температуры по термометрической скважине т.скв.1 (показана на схеме расположения элементов модели), расположенной под зданием приведены в таблице.

Глубина, мТемпература грунта по термометрической скважине, ⁰С
На начало расчета
(15 мая)
На конец летних периодов (15 октября)
1 год2 год3 год
1,01,337,917,797,70
2,01,850,980,36-0,09
3,01,220,18-0,51-1,15
4,00,23-0,52-1,38-2,11
5,0-0,58-1,10-2,10-2,84
6,0-0,80-1,55-2,60-3,34
7,0-0,90-1,86-2,89-3,60
8,0-0,96-1,99-2,97-3,66
9,0-0,99-1,97-2,86-3,53
10,0-0,99-1,80-2,61-3,25
11,0-0,99-1,55-2,25-2,86
12,0-0,99-1,31-1,86-2,41
13,0-0,98-1,13-1,51-1,97
14,0-0,98-1,03-1,25-1,60
15,0-0,97-0,99-1,10-1,33

Время выполнения расчета модели, состоящей из 4,9 млн. ячеек на 30-ти летний период составляет:

  • На центральном процессоре (Intel i7-6700К) — 32 часа.
  • На видеокарте (GeForce GTX 1080) – 1 час 18 минут.

Выводы

Исходя из представленных расчетов можно сделать следующие выводы:

  • Выполнение отсыпки талым грунтом в конце летнего периода приводит к повышению температур грунтов в верхней части геокриологического разреза. Восстановление исходного температурного состояния происходит в течении 12-15 лет эксплуатации.
  • С помощью только наклонных СОУ, установленных под углом до 35°, не удается проморозить массив талого грунта, получившегося в результате отсыпки талым грунтом, в первые годы эксплуатации сооружения.

Файлы для загрузки

Файлы примера доступны по ссылке.

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *

десять − два =