Система термостабилизации грунтов оснований
резервуарного парка V=2000 м³

Рассмотрим выполнение прогнозного моделирования температурного поля грунтов основания вертикальных резервуаров V=2000 м3. Рассматриваются два технических решений по термостабилизации грунтов оснований:

  • проветриваемое подполье;
  • проветриваемое подполье совместно с одиночными сезонно-действующими охлаждающими устройствами (термостабилизаторами).

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

Резервуарный парк V=2000. Общий вид

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

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

ПоказательМесяцы
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 ⁰С.

Резервуарный парк V=2000. Геология
Свойства грунтов приведены в таблице ниже.
Наименование
показателя
ОбозначениеЕд.измеренияНасыпной грунтИГЭ №140200.
Суглинок
ИГЭ №140300.
Суглинок
ИГЭ №141100Э.
Суглинок
Физические свойства
Плотность грунта в сухом состоянииRdfкг/м31740170017001700
Суммарная влажностьWtotд.е.0.180.2470.2760.256
Степень засоленностиDsal%----
Число пластичностиIpд.е.-0.070.070.07
Влажность на границе раскатыванияWpд.е.-0.170.170.17
Теплофизические свойства
Температура начала замерзанияTbf⁰С-0.15-0.21-0.21-0.21
Теплопроводность талого грунтаλthВт/(м*К)2.321.801.801.80
Теплопроводность мерзлого грунтаλfВт/(м*К)2.732.022.022.02
Объемная теплоемкость талого грунтаCthМДж/м32.633.243.453.30
Объемная теплоемкость мерзлого грунтаCfМДж/м31.932.492.592.52
Коэффициенты кривой незамерзшей воды
Ww(t)= A + B / (C - t),
где t - температура, ⁰С
A-0.00130.04780.04780.0478
B-0.00160.05710.05710.0571
C-0.44520.64250.64250.6425
Теплота фазового переходаQfМДж/м3----

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

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

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

Проектирование термостабилизации грунтов оснований ведется по двум вариантам:

  • Вариант 1. Проветриваемое подполье.
    Данный вариант предусматривает использование только проветриваемого подполья для термостабилизации грунтов оснований.
  • Вариант 2. Проветриваемое подполье + вертикальные термостабилизаторы.
    Данный вариант предусматривает совместное использование проветриваемого подполья с установкой одиночных вертикальных

    сезоннодействующих охлаждающих устройств (термостабилизаторов).

Схема расположения свай и термостабилизаторов в основании резервуарного парка V=2000 м3

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

Прогнозное моделирование температурного поля грунтов осуществлялось численными методами с учетом требований к порядку составления прогноза изменения температурного режима грунтов РСН 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.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
Теплоотдача, Вт/(м2*К)9.539.7710.2510.7311.4511.4510.259.5310.0110.259.539.29
Высота снега, м0.1580.1850.2030.1940.084----0.0240.0750.119
Плотность снега, кг/м3118123127134307----94103111
Теплопроводность снега, Вт/(м*К)0.1780.1830.1860.1930.483----0.1530.1620.170
Коэффициент теплопередачи, Вт/(м2*К)1.0050.8980.8430.9113.84311.4510.259.5310.013.9521.7711.237
Вариант 1. Проветриваемое подполье
Расчетная область

Расчетная область для численного моделирования температурного поля грунтов представляет собой параллелепипед с размерами 80 х 60 х 50 метров. Расчетная область показана на следующем рисунке.

3D модель для расчета термостабилизации грунтов резервуарного парка V=2000 м3

Схема расположения свай и строительных конструкций в модели приведена на рисунке.

Резервуарный парк V=2000. без термостабилизаторов. Схема свай

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

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

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

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

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

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

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

Трехмерное температурное поле в основаниирезервуарного парка V=2000 м3

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

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

Отепляющее воздействие свай на многолетнемерзлые грунты оснований резервуарного парка V=2000 м3. В начальный момент времени.

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

Отепляющее воздействие свай на многолетнемерзлые грунты оснований резервуарного парка V=2000 м3. Осень первого года эксплуатации.

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

Отепляющее воздействие свай на многолетнемерзлые грунты оснований резервуарного парка V=2000 м3. Осень пятого года эксплуатации.

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

Отепляющее воздействие свай на многолетнемерзлые грунты оснований резервуарного парка V=2000 м3. Десятый год эксплуатации.

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

Глубина, мТемпература грунта по термометрической скважине, ⁰С
На начало расчета
(15 октября)
На конец летних периодов (15 октября)
1 год5 год10 год
1.0-0.02-0.08-0.09-0.11
2.0-0.16-0.15-0.27-0.51
3.0-0.20-0.18-0.42-1.08
4.0-0.36-0.21-0.29-1.33
5.0-0.58-0.21-0.22-1.35
6.0-0.74-0.21-0.28-1.23
7.0-0.85-0.23-0.35-1.06
8.0-0.91-0.34-0.43-0.93
9.0-0.94-0.49-0.51-0.85
10.0-0.97-0.65-0.58-0.80
11.0-0.97-0.77-0.65-0.77
12.0-0.98-0.85-0.71-0.77
13.0-0.97-0.91-0.76-0.77
14.0-0.97-0.93-0.80-0.79
15.0-0.97-0.95-0.83-0.80

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

  • На центральном процессоре (Intel i7-6700R) — 15 часов 35 минут.
  • На видеокарте (GeForce GTX 1080) — 29 минут.
Вариант 2. Проветриваемое подполье + вертикальные термостабилизаторы (СОУ)
Расчетная область

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

Моделирование термостабилизации грунтов резервуарного парка V=2000 м3

Схема расположения свай и строительных конструкций в модели приведена на рисунке.

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

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

ПоказательМесяцы
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
Скорость ветра, м/с0.10.10.10.10.10.10.10.10.10.10.10.1
Коэф. теплоотдачи к поверхности конденсатора, Вт/(м2*К)10.2410.2410.2410.2410.24----10.2410.2410.24
Эффективный коэф. теплопередачи, отнесенный к наружной поверхности испарителя, Вт/(м2*К)9.279.279.279.279.27----9.279.279.27

Расчетная область состоит из 370х223х47 ячеек (3,9 млн. ячеек).

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

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

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

Трехмерное температурное поле грунтов в основании резервуарного парка V=2000 м3 с учетом работы системы термостабилизации грунтов

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

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

Результат работы вертикальных термостабилизаторов

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

Глубина, мТемпература грунта по термометрической скважине, ⁰С
На начало расчета
(15 октября)
На конец летних периодов (15 октября)
1 год5 год10 год
1.0-0.02-0.09-0.13-0.14
2.0-0.16-0.20-1.03-1.33
3.0-0.20-0.26-2.21-2.79
4.0-0.36-0.81-3.08-3.92
5.0-0.58-1.24-3.65-4.69
6.0-0.74-1.57-3.96-5.14
7.0-0.85-1.79-4.06-5.31
8.0-0.91-1.89-3.97-5.23
9.0-0.94-1.89-3.74-4.98
10.0-0.97-1.83-3.44-4.63
11.0-0.97-1.67-3.06-4.17
12.0-0.98-1.49-2.65-3.67
13.0-0.97-1.31-2.23-3.15
14.0-0.97-1.17-1.84-2.65
15.0-0.97-1.07-1.54-2.21

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

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

Выводы

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

  • Выполнение отсыпки талым грунтом в конце летнего периода приводит к повышению температур грунтов в верхней части геокриологического разреза. Восстановление исходного температурного состояния происходит в течении 12-15 лет эксплуатации.
  • Тепловыделения от свай при заполнении свай цементно-песчаными раствором и гидратации цемента в процессе твердения при водит к повышению температур грунтов. Восстановление исходного температурного состояния, в результате работы проветриваемого подполья, происходит в течении 5-6 лет эксплуатации.
  • В результате работы СОУ происходит промораживание грунтов и понижение температур грунтов. Компенсация теплового воздействия от свай фундамента и выполнение отсыпки талым грунтом происходи за 3 месяца работы СОУ. В целом, объем погруженных СОУ избыточен.

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

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

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

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

9 − два =