сотрудник
Северо-Восточный федеральный университет имени М.К. Аммосова (кафедра вычислительных технологий, Доцент)
сотрудник с 01.01.2019 по 01.01.2022
Якутск, Республика Саха (Якутия), Россия
Россия
Россия
ВАК 05.02.2000 Машиностроение и машиноведение
ВАК 05.17.00 Химическая технология
ВАК 05.23.00 Строительство и архитектура
УДК 60 Прикладные науки. Общие вопросы
ГРНТИ 81.01 Общие вопросы технических, прикладных наук и отраслей экономики
ОКСО 01.03.02 Прикладная математика и информатика
ББК 30 Техника и технические науки в целом
ТБК 50 Технические науки в целом
BISAC MAT007000 Differential Equations / General
Рассмотрен процесс оттаивания-промерзания вмещающего грунта подземного полиэтиленового трубопровода теплоснабжения с теплоизоляцией, эксплуатируемого в районах распространения многолетнемерзлых грунтов. Поставлена задача определения толщины теплоизоляции трубопровода и величины его заглубления на основе анализа температурного поля в грунте в течение многолетних циклов. Температурное поле в системе «теплоизолированная труба-грунт», изменяющееся во времени, предложено определять, решая уравнения теплопроводности в полярных координатах. При этом пучок труб рассматривается как одна труба эквивалентного радиуса. Такой подход позволил упростить условия сопряжения трубы, теплоизоляции, полиэтиленовой оболочки и вмещающего грунта, а также использовать простой численный метод решения уравнения – метод конечных разностей. Приведены результаты расчетов максимальной глубины оттаивания грунта от величины заглубления трубы теплоснабжения при различных толщинах теплоизоляции. Предложено определять оптимальную величину заглубления подземного трубопровода теплоснабжения из условия ежегодного установления нулевой изотермы вмещающего грунта на уровне границы деятельного слоя грунта к началу отопительного сезона. Возможность нахождения нулевой изотермы на таком уровне проиллюстрирована расчетами изотерм к началу отопительного сезона (сентябрь) при различных заглублениях и толщинах теплоизоляции. По результатам расчетов предложена методика расчетного определения толщины слоя теплоизоляции и величины заглубления трубопровода из условия минимизации функционала, характеризующего меру уклонения ординаты нулевой изотермы на левой границе области от уровня глубины деятельного слоя грунта в момент времени начала отопительного сезона.
трубопровод, теплоснабжение, многолетнемерзлые грунты, температура, глубина оттаивания, толщина теплоизоляции
Введение. В настоящее время в практику строительства тепловых сетей широко внедряются опытно-промышленные бесканальные варианты внутриквартальных предизолированных полиэтиленовых подземных трубопроводов теплоснабжения. В то же время, проектные решения прокладки таких трубопроводов недостаточно научно обоснованы. Это связано, прежде всего, с отсутствием методов определения таких параметров теплопровода, как толщина теплоизоляции и величина заглубления, которыми можно управлять тепловым режимом грунта в основании трубопровода. Управление тепловым режимом грунта имеет особую важность при проектировании тепловых сетей в регионах с многолетнемерзлыми грунтами.
Согласно нормативным документам расчет толщины теплоизоляции осуществляется по нормированной плотности теплового потока [1]. Такой расчет, при котором за температуру окружающей среды берется среднегодовая температура грунта на глубине заложения трубопровода, не пригоден для расчетов толщины теплоизоляции теплопроводов, прокладываемых в регионах распространения многолетнемерзлых грунтов. Кроме того, при таких расчетах толщина теплоизоляции не связана с величиной заглубления теплопроводов, что не соответствует реальному процессу теплового режима в многолетнемерзлых грунтах. При проектировании тепловых сетей в районах распространения многолетнемерзлых грунтов толщину слоя теплоизоляции необходимо определять в комплексе с величиной заглубления теплопровода, основываясь на анализе динамики границы оттаивания грунта. Отсутствие подобной методики расчета приводит к негативным последствиям при практическом использовании перспективных трубопроводов теплоснабжения подземной прокладки.
В условиях эксплуатации подземного теплопровода тепловое состояние грунта в течение годового цикла претерпевает значительные изменения. Это обусловлено не только изменением температуры окружающего воздуха, но и температурным режимом теплоносителя. На динамику температурного поля грунта существенное влияние имеют также толщина слоя теплоизоляции и глубина заложения трубопровода. В регионах распространения многолетнемерзлых грунтов прокладка теплопровода с использованием недостаточно обоснованных параметров может привести к ежегодному увеличению ореола оттаивания вокруг теплопровода и потере его устойчивости. Имеется достаточно большое количество научных работ, посвященных исследованию температурного поля в основании бесканальных полимерных трубопроводов [2–4]. В то же время актуальной проблеме сохранения многолетнемерзлых грунтов при воздействии подземного теплопровода уделяется недостаточное внимание.
В целях снижения теплопотерь подземные трубопроводы теплоснабжения прокладываются в пучке с другими трубопроводами, который можно рассматривать как один трубопровод [5]. Тогда температурное поле в многослойной трубе и грунте можно описать уравнением теплопроводности в полярных координатах, упрощающей условия сопряжения трубы, теплоизоляции, полиэтиленовой оболочки и вмещающего грунта.
Для учета теплоты фазового перехода воды в уравнении теплопроводности удобнее воспользоваться фронтовой моделью, согласно которой фазовое превращение происходит на границе мерзлого и талого грунта [6–9]. При численном решении уравнения теплопроводности, учитывающего фазовое превращение, трудность представляет определение границ твердой и жидкой фаз. Эффективным методом решения подобных уравнений является метод сквозного счета, предложенный в работе [10] и применяемый для решения многих прикладных задач [11–14]. Численные алгоритмы решения задачи определения динамики температурного поля грунта вокруг подземного полимерного трубопровода с теплоизоляцией рассматривались в работах [15, 16]. Результаты решения задач теплового режима в системе «грунт-труба» не применялись для решения задач проектирования теплопроводов, в частности определения толщины теплоизоляции и глубины заглубления.
Целью данной работы является разработка методики определения толщины теплоизоляции и величины заглубления подземного трубопровода теплоснабжения, прокладываемого в регионах с многолетнемерзлыми грунтами на основе анализа динамики нулевой изотермы (изотермы 0 °С) в области вмещающего грунта.
Постановка задачи. В рассматриваемой области (рис. 1) запишем уравнение теплопроводности в полярных координатах:
. (1)
На поверхности грунта граничное условие третьего рода запишем в виде:
(2)
В условии (2) эффективный коэффициент теплообмена αeff в зимний период вычисляется в зависимости от толщины снега, Tair– температура окружающего воздуха.
На удалении от теплопровода условие отсутствия теплового потока в горизонтальном направлении запишем в виде:
(3)
В виду симметричности температурного поля справедливы условия:
; (4)
На внутренней поверхности теплопровода происходит конвективный теплообмен с теплоносителем с температурой Тwater:
. (5)
На границах разнородных тел запишем условия идеального теплового контакта:
; (6)
, (7)
где i – номер слоя.
При t=0 зададим начальное условие:
. (8)
В уравнении (1) с начальным и граничными условиями использованы следующие обозначения: Т – температура; ρ(T), l(T) – эффективные (сглаженные) коэффициенты теплоемкости и теплопроводности; r, j – полярные координаты; t – время Поставленная задача решалась методом конечных разностей. Расчеты проводились с шагом по времени равным 3 часа. Проводились предварительные расчеты для определения начального распределения температуры следующим образом. За начальное время бралось первое число сентября. За начальное условие принималось однородное распределение температуры в рассматриваемой области с температурой -3 °С Проводились расчеты 3 летнего цикла, в результате которого распределение температуры в рассматриваемой области в начале сентября становится практически не зависящим от начального условия. Полученное температурное поле принималось за начальное условие при расчетах.
Рис. 1. Расчетная область предизолированного трубопровода и грунта:
Ω – область; Г1 – дневная поверхность; Г2, Г3, Г4, Г5 – границы области, j – угловая координата
Результаты расчетов. Обсуждения. Методику определения толщины теплоизоляции и величины заглубления подземного трубопровода теплоснабжения изложим на примере расчета температурного режима вмещающего грунта при воздействии трубопровода теплоснабжения с эквивалентным радиусом 81,5 мм.
В таблице 1 приведены значения теплофизических характеристик материалов, принятые в расчетах [17–18].
Таблица 1
Теплофизические характеристики грунта и материалов теплопровода
|
Параметр |
Труба из сшитого полиэтиленa |
Изоляция из пенополиуретана |
Оболочка из полиэтилена |
Талый грунт |
Мерзлый грунт |
|
Коэффициент |
0,35 |
0,05 |
0,42 |
1,4 |
1,5 |
|
Плотность, кг/м3 |
938 |
33 |
960 |
1700 |
1700 |
|
Удельная теплоемкость, Дж/(кг×К) |
2300 |
1800 |
1700 |
1850 |
1750 |
Примем, что температура теплоносителя
Twater в летний период (май-август) равна 10 °С, в отопительный период- меняется от 80 до 95 °С в зависимости от температуры окружающего воздуха. Влажность грунта равна 0,18. Коэффициент теплообмена α на внутренней поверхности трубы вычислялся по формуле [19] :
(9)
где ν – кинематическая вязкость воды, a – коэффициент температуропроводности воды, w – скорость течения воды в трубе; r1 – внутренний радиус трубы.
Эффективный коэффициент теплообмена αeff на дневной поверхности грунта определялся из формулы [7–9]:
, (10)
где αair – коэффициент теплообмена поверхности с окружающим воздухом; hsnow, lsnow – толщина и коэффициент теплопроводности снежного покрова.
В начальный момент времени расчета (1 апреля) температурное поле в грунте не известно. В качестве начального распределения температуры задавалось однородное температурное поле со значением -3 °С. Расчеты показали, что температурное поле в исследуемой области после нескольких годовых циклов практически не зависит от начального распределения температуры. В связи с этим далее приводятся результаты расчетов после трехлетнего цикла.
Изменение температуры окружающего воздуха во времени описывалось функцией среднесуточных температур, полученной статистической обработкой двадцатилетних температурных данных:
(11)
Заглубление тепловых сетей до верха оболочки при бесканальной прокладке полиэтиленовых трубопроводов согласно нормативным документам принималось не менее 0,7 м. Исследовалось влияние величины заглубления подземного трубопровода теплоснабжения на тепловое состояние вмещающего грунта. Варьировались величина заглубления трубопровода в интервале [70, 130] с шагом 20 см и толщина теплоизоляции в интервале [2, 5] с шагом 1 см. Расчетами с многолетними циклами показано, что, несмотря на минимальные температуры окружающего воздуха, максимальная глубина оттаивания достигается в январе. На рис. 2 представлены зависимости максимальной глубины оттаивания грунта от величины заглубления трубы теплоснабжения при различных толщинах теплоизоляции. Увеличение заглубления и уменьшение толщины теплоизоляции приводят к увеличению глубины оттаивания.
Рис. 2. Зависимости максимальных глубин оттаивания xот величины заглубления трубы Н при различных
толщинах теплоизоляции d
Расчеты показывают, что при всех рассматриваемых толщинах теплоизоляции и величинах заглубления нулевая изотерма расположена в области теплоизоляции, что свидетельствует об отсутствии замерзания теплоносителя. Такое поведение нулевой изотермы не позволяет при заданной толщине теплоизоляции величину необходимого заглубления определять из условия недопущения замораживания теплоносителя. Наиболее практичным является определение оптимального заглубления подземного трубопровода теплоснабжения из условия ежегодного установления нулевой изотермы вмещающего грунта на уровне границы деятельного слоя грунта к началу отопительного сезона. Такое заглубление будет оптимальным с точки зрения сохранения стабильной глубины оттаивания в окрестности трубопровода в течение продолжительного времени, поскольку процесс оттаивания грунта под воздействием теплопровода ежегодно будет повторяться с одного и того же уровня глубины. Увеличение толщины теплоизоляции выше расчетной при фиксированной величине заглубления приводит к смещению нулевой изотермы в основании трубопровода выше границы деятельного слоя, т.е. образованию под трубой валика из мерзлого грунта. Такая граница оттаивания грунта будет препятствовать миграции влаги в горизонтальном направлении и может привести к негативным последствиям.
Численными экспериментами, при которых варьировалась толщина теплоизоляции, для рассматриваемых климатических условий была установлена минимальная толщина теплоизоляции равная 3 см, обеспечивающая установление к началу отопительного сезона нулевой изотермы в грунте на уровне границы деятельного слоя (рис. 3–6). Координата глубины деятельного слоя из рисунков определяется по значению ординаты нулевой изотермы на расстоянии 8 м от начала координат.
К началу отопительного сезона нулевая изотерма в области вмещающего грунта восстанавливается до уровня границы деятельного слоя при величине заглубления трубопровода 70 см и толщине теплоизоляции 3 см. При заглублении трубопровода 90 см толщина теплоизоляции должна быть 4 см, при заглублении 110 см – 5 см. При заглублении H=130 см толщина теплоизоляции h должна быть более 5 см (рис. 6) и может быть рассчитана по формуле
.
|
А) |
Б) |
|
В) |
Г) |
Рис. 3. Изотермы температур вмещающих грунтов трубопровода теплоснабжения к началу отопительного
сезона (сентябрь) при заглублении 70 см при различных толщинах теплоизоляции d :
А – d=2; Б – d=3; В – d=4; Г – d= 5 см
|
А) |
Б) |
|
В) |
Г) |
Рис. 4. Изотермы температур вмещающих грунтов трубопровода теплоснабжения к началу отопительного
сезона (сентябрь) при заглублении 90 см при толщинах теплоизоляции d: А-Г соответствуют d рис. 3
|
А) |
Б) |
|
В) |
Г) |
Рис. 5. Изотермы температур вмещающих грунтов трубопровода теплоснабжения к началу отопительного
сезона (сентябрь) при заглублении 110 см при различных толщинах теплоизоляции d:
А–Г соответствуют d рис. 3
|
А) |
Б) |
|
В) |
Г) |
Рис. 6. Изотермы температур вмещающих грунтов трубопровода теплоснабжения к началу отопительного
сезона (сентябрь) при заглублении 130 см при различных толщинах теплоизоляции d :
А–Г соответствуют d рис. 3
Предлагаемая методика определения параметров подземного трубопровода теплоснабжения (величины заглубления Ни толщины слоя теплоизоляцииd) сводится к минимизации функционала:
, (12)
где y – уравнение в декартовых координатах, описывающее нулевую изотерму, зависящее от абсциссы, величины заглубления, толщины слоя теплоизоляции; ts – время начала отопительного сезона, l – ширина рассматриваемой области.
Функционал (12), характеризующий меру уклонения ординаты нулевой изотермы на левой границе области от уровня глубины деятельного слоя грунта в момент времени начала отопительного сезона. Минимизация функционала (12) реализуется методом последовательного анализа вариантов.
Заключение. Предлагаемая методика определения параметров трубопровода теплоснабжения из условия ежегодного установления нулевой изотермы вмещающего грунта к началу отопительного сезона на уровне границы деятельного слоя грунта, может применяться при проектировании подземных тепловых сетей из полиэтиленовых предизолированных труб в регионах распространения многолетнемерзлых грунтов.
1. СНиП 41-03-2003. Тепловая изоляция оборудования и трубопроводов. М.:2004.
2. Половников В.Ю., Хузеев В.А. Численный анализ влияния промерзания грунта в зоне прокладки на тепловые потери бесканальных теплопроводов // Инженерно-строительный журнал. 2013. №2. С. 16-24.
3. Размазин Г.А., Моисеев Б.В. Тепловое взаимодействие бесканальной прокладки теплопроводов с вечномерзлыми грунтами // Проблемы строительства, инженерного обеспечения и экологии городов / Сб. материалов I Международной научно-практической конференции. Пенза, ПДЗ, 2000. С.106-110.
4. Половников В.Ю. Численный анализ влияния нестационарности процессов теплопереноса в зонах размещения подземных теплопроводов на их тепловые режимы и тепловые потери // Известия Томского политехнического университета. Инжиниринг георесурсов. 2018. Т. 329. № 10. С. 76-84.
5. Пехович А.И., Жидких В.М. Расчеты теплового режима твердых тел. Л.: «Энергия», 1976. 352 с.
6. Andersland O.B, Ladanyi B. An introduction to frozen ground engineering. Springer Science & Business Media, 2013.
7. Слепцов В.И., Мордовской С.Д., Изаксон В.Ю. Математическое моделирование теплообменных процессов в многолетнемерзлых горных породах. Н.: Наука. Сибирская издательская фирма РАН, 1996.104 с.
8. Вабищевич П.Н., Варламов С.П., Васильев В.И., Васильева М.В., Степанов С.П. Численное моделирование температурного поля многолетнемерзлого грунтового основания железной дороги // Математическое моделирование. 2016. Т.28, №10. С. 110-124.
9. Васильев В.И, Максимов А.М., Петров Е.Е., Цыпкин Г.Г. Тепломассоперенос в промерзающих и протаивающих грунтах. М.: Наука, Физматлит, 1997. 224 с.
10. Самарский А.А., Моисеенко Б.Д. Экономичная схема сквозного счета для многомерной задачи Стефана // ЖВМ и МФ. 1965. Т.5.№5. С. 816-827.
11. Gornov V.F., Stepanov S.P., Vasilyeva M.V. etal. Mathematical modeling of heat transfer problems in the permafrost // AIP Conference Proceedings.2014. Vol. 1629. Pp. 424-431.
12. Pavlova N. V., Vabishchevich P. N., Vasilyeva M. V. Mathematical modeling of thermal stabilization of vertical wells on high performance computing systems // International Conference on Large-Scale Scientific Computing / Springer. 2013. Pp. 636-643.
13. Ascher U.M. Numerical methods for evolutionary differential equations. Society for Industrial Mathematics, 2008.
14. Vabishchevich P.N., Vasilyeva M.V., Pavlova N.V. Numerical simulation of thermal stabilization of filter soils // Mathematical Models and Computer Simulations. 2015. Vol. 7, No. 2. Pp. 154-164.
15. Акимов М.П., Мордовской С.Д., Старостин Н.П. Численный алгоритм для исследования влияния бесканального подземного трубопровода теплоснабжения на вечномерзлые грунты // Математические заметки ЯГУ. 2010. Т.17. №2. С.125-131.
16. Акимов М.П., Захаров П.Е., Матвеева О.И. Численное моделирование динамики температурного поля многолетнемерзлых грунтов при воздействии трубопроводов // Математические заметки СВФУ. 2014. Т. 21. №4. С. 61-70.
17. Гаврильев Р.И. Теплофизические свойства горных пород и напочвенных покровов криолитозоны. Н.: Изд-во СО РАН, 1998. 280 с.
18. Бабичев А.П., Бабушкина Н.А., Братковский А.М. и др. Физические величины. Справочник. М.: Энергоатомиздат, 1991. 1232 с.
19. Кутателадзе С.С. Основы теории теплообмена. Изд. 5-е. М.: Атомиздат, 1979. 416 с.
