Резервуар для ливневых стоков

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

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

Резервуар для ливневых стоков. Общий вид

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

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

ПоказательМесяцы
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 метров составляет минус 0,5 ⁰С.

Резервуар для ливневых стоков. Геология

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

Наименование
показателя
ОбозначениеЕд.измеренияНасыпной грунтИГЭ №141100.
Суглинок
слабольдистый незасоленый
ИГЭ №140300.
Суглинок мягкопластичный
ИГЭ №141000Э.
Суглинок нельдистый
ИГЭ №141100Э.
Суглинок
слабольдистый незасоленый
Физические свойства
Плотность грунта в сухом состоянииRdfкг/м³17401450148018401500
Суммарная влажностьWtotд.е.0,180,3060,2830,1490,286
Степень засоленностиDsal%-0,0490,0570,0590,051
Число пластичностиIpд.е.-0,1220,1320,1020,123
Влажность на границе раскатыванияWpд.е.-0,2030,2050,1810,207
Теплофизические свойства
Температура начала замерзанияTbf⁰С-0,15-0,12-0,2-0,17-0,15
Теплопроводность талого грунтаλthВт/(м*К)2,321,591,531,551,59
Теплопроводность мерзлого грунтаλfВт/(м*К)2,731,811,761,761,81
Объемная теплоемкость талого грунтаCthМДж/м³2,633,183,112,793,16
Объемная теплоемкость мерзлого грунтаCfМДж/м³1,932,352,352,362,38
Коэффициенты кривой незамерзшей воды
Ww(t)= A + B / (C - t),
где t - температура, ⁰С
А-0,00130,07420,07810,06050,0761
B-0,00160,110,1040,16570,1085
C-0,44520,58520,48511,70060,5444
Теплота фазового переходаQfМДж/м³-----

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

Резервуар для ливневых стоков представляет собой железобетонную конструкцию в виде параллелепипеда. Под монолитным днищем выполнена бетонная подготовка толщиной 100 мм из бетона класса B15.

Грунты основания используются по I принципу, т.е. с сохранением мерзлого состояния в процессе строительства и эксплуатации.

Температурная стабилизация

Технические решения по температурной стабилизации предусматривают сооружение теплоизоляционного экрана под днищем резервуара.

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

Прогнозное моделирование осуществлялось с учетом требований к порядку составления прогноза изменения температурного режима грунтов, РСН 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,2990,3500,3840,3670,158----0,0450,1410,226
Плотность снега, кг/м³144153159166386----98115131
Теплопроводность снега, Вт/(м*К)0,2070,2180,2260,2360,672----0,1570,1740,191
Коэффициент теплопередачи, Вт/(м²*К)0,6440,5860,5570,6063,10311,45010,2509,53010,0102,5981,0940,776
Расчетная область

Модель расчетной области представляет собой параллелепипед с размерами 70 х 60 х 50 метров.

Расчетная область показана на рисунке ниже.

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

Резервуар для ливневых стоков. 3D Модель

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

Тепловое влияние резервуара моделировалось граничными условиями 3-го рода. Предполагается что слив дождевого стока в резервуар осуществляется в весенний и осенний период в течении 15 суток, а в летний период по 5 суток каждого месяца. В период слива дождевого стока в резервуар коэффициент теплоотдачи на внутренней поверхности резервуара принимается равным 20 Вт/(м²*К). В период отсутствия дождевого стока коэффициент теплоотдачи рассчитывается как коэффициент теплопередачи, т.е. как обратная величина суммы термических сопротивлений снежного покрова, сопротивления теплоотдачи ко внутреннему пространству резервуара и сопротивления теплоотдачи от внутреннего пространства к стенке. График коэффициента теплопередачи приведен на рисунке.

График теплоотдачи резервуара

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

График тепловыделения бетона B15

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

Заполнение ячеек расчетной сетки материалами приведено на рисунке ниже.

Резервуар для ливневых стоков. Материалы
Результаты расчета

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

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

Трехмерное температурное поле ММГ в основании резервуара на 15.09.2003

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

Температурное поле грунтов на 15 мая 1-го года (начало эксплуатации).

Температура ММГ в основании резервура на 15.05.2000

Температурное поле грунтов на 15 октября 1-го года.

Температура ММГ в основании резервуара на 15.10.2000

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

Температура ММГ в основании резервуара на 15.10.2001

Температурное поле грунтов на 15 октября 3-го года.

Температура ММГ в основании резервуара на 15.10.2002

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

Глубина, мТемпература грунта по термометрической скважине, ⁰С
На начало расчета
(15 мая)
На конец летних периодов (15 октября)
1 год2 год3 год
1,0-0,50-0,01-0,02-0,02
2,0-0,50-0,04-0,04-0,04
3,0-0,50-0,12-0,12-0,12
4,0-0,50-0,19-0,23-0,25
5,0-0,54-0,31-0,37-0,40
6,0-0,46-0,42-0,49-0,55
7,0-0,41-0,43-0,48-0,53
8,0-0,41-0,42-0,45-0,49
9,0-0,42-0,43-0,44-0,47
10,0-0,44-0,45-0,45-0,47
11,0-0,47-0,47-0,47-0,48
12,0-0,49-0,49-0,49-0,49
13,0-0,51-0,51-0,50-0,51
14,0-0,53-0,53-0,52-0,52
15,0-0,54-0,54-0,54-0,54

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

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

Выводы

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

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

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

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

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

1 × 3 =