Известия РАН. Физика атмосферы и океана, 2021, T. 57, № 5, стр. 526-538

Зависимость антропогенного потока тепла от температуры воздуха (на примере Санкт-Петербурга)

А. С. Гинзбург a*, И. А. Евсиков b, В. А. Фролькис bcd

a Институт физики атмосферы им. А.М. Обухова РАН
119017 Москва, Пыжевский пер., 3, Россия

b Санкт-Петербургский государственный архитектурно-строительный университет
190005 Санкт-Петербург, 2-я Красноармейская ул., 4, Россия

c Главная геофизическая обсерватория им. А.И. Воейкова
194021 Санкт-Петербург, ул. Карбышева, 7, Россия

d Санкт-Петербургский государственный экономический университет
191023 Санкт-Петербург, ул. Садовая, 21, Россия

* E-mail: gin@ifaran.ru

Поступила в редакцию 25.05.2021
После доработки 07.06.2021
Принята к публикации 09.06.2021

Полный текст (PDF)

Аннотация

Антропогенный поток тепла, связанный с расходом энергии на отопление зданий, оказывает существенное влияние на тепловой баланс урбанизированных территорий и интенсивность городского острова тепла. Энергопотребление городского хозяйства сильно зависит от метеорологических условий, их годового и суточного хода. Особенно это касается той части энергопотребления, которая расходуется на поддержание комфортной температуры внутри помещений. В статье используется открытая веб-картографическая платформа “OpenStreetMap”, которая позволяет провести инвентаризацию потоков антропогенного тепла по разности внутренней и внешней температур воздуха на основе строительных нормы и теплофизических свойств ограждающих конструкций. Сформулирована простая нестационарная модель взаимодействия городского пограничного слоя атмосферы с городским покрывающим слоем, содержащим антропогенные источники тепла.

Ключевые слова: энергопотребление городского хозяйства, теплофизические свойства зданий, антропогенный поток тепла, городской покрывающий слой, характеристики отопительного периода, энергобалансовые модели

ВВЕДЕНИЕ

Воздействию антропогенной эмиссии тепла на атмосферу урбанизированных территорий в последнее время посвящено большое количество публикаций, обзор которых можно найти, например, в статье А.С. Гинзбурга и П.Ф. Демченко (Известия РАН. Физика атмосферы и океана, 2019, № 6). Главной особенностью атмосферы городской агломерации является наличие переходного слоя между поверхностью и верхней границей городских строений – городского покрывающего слоя, в котором происходит тепловое загрязнение городской атмосферы, влияющее на структуру и свойства всего атмосферного пограничного слоя. Выделяемая в атмосферу часть потребляемой энергии, отнесенная к единице площади подстилающей поверхности, называется антропогенным потоком тепла (АПТ) – QF и состоит из потоков тепла от зданий, транспорта и метаболизма живых организмов [1]. На урбанизированных территориях основной вклад в АПТ дает энергопотребление зданий. Состояние современных исследований климата больших городов и роли антропогенного потока тепла подробно описано в [2], а для Москвы в [3].

Энергия, затраченная на поддержание температуры внутри зданий в комфортном диапазоне, поступает в окружающее пространство в виде АПТ. Многочисленные эмпирические оценки и численные эксперименты показали, что АПТ может увеличивать температуру воздуха в городском покрывающем слое, который расположен между поверхностью и огибающей крыш зданий, и давать существенный вклад в интенсивность городского острова тепла, который определяется разностью температур между городскими и прилегающими сельскими территориями. Результаты таких исследований для городов России описаны в работах [412] и других публикациях. Зимой интенсивность городского острова тепла на 50–80% определяется наличием антропогенного потока тепла [8].

Доля АПТ, которая связана с энергией, потребляемой внутри зданий, существенно зависит от температуры воздуха снаружи помещений. Из-за влияния температуры атмосферного воздуха на энергопотребление и обратного влияния АПТ на температуру приземного воздуха в режиме отопления возникает отрицательная обратная связь между температурой воздуха и антропогенным потоком тепла. Наиболее современным методом исследования такой обратной связи является использование моделей мезомасштабной циркуляции атмосферы (WRF, COSMO и др.) совместно с моделями городской застройки и городского покрывающего слоя и моделями энергетики зданий.

Воздействие АПТ на атмосферу во многом определяется тем, как этот поток тепла распределяется между компонентами теплового баланса поверхности и городского покрывающего слоя. Часть АПТ передается теплопередачей из помещений в атмосферу через стены и крыши зданий. Если действие АПТ при моделировании атмосферного пограничного слоя ограничивается фиксацией температуры внутри зданий, то учитывается только вклад теплового потока зданий. Эту часть выделяют в отдельное слагаемое – тепловой поток зданий. Он модифицирует поток тепла через подстилающую поверхность, который в уравнении энергетического баланса учитывается как поток теплового накопления и отвечает за тепловую инерцию элементов ландшафта.

Часть АПТ поступает в атмосферу и (частично) в грунт и водные резервуары в виде тепловых стоков от работы систем нагревания. Эта часть АПТ называется потоком теплового стока (WHF – waste heat flux). В случае поступления непосредственно в атмосферу, эта часть АПТ дает вклад в атмосферный поток явного и скрытого тепла.

Помимо сезонного хода, взаимодействие городской атмосферы и энергопотребления подвержено и суточному циклу. Это во многом связано с суточным циклом поступающей к поверхности солнечной радиации. Также может давать вклад суточный цикл деловой активности. Численные эксперименты показывают, что максимум реакции температуры воздуха городского покрывающего слоя на АПТ достигается ночью. Аналогичное усиление парникового эффекта в относительно холодное время (ночью и зимой) имеет место при росте концентрации СО2 в атмосфере [13].

Эффект ночного увеличения реакции на АПТ по численным экспериментам продолжается иногда и в дневное время. Это связано, по-видимому, с комбинированным действием двух факторов. Во-первых, с сохранением дневных аномалий за счет потока теплового накопления. Во-вторых, с подавлением вертикального перемешивания при переходе к устойчивой стратификации при исчезновении солнечного нагрева поверхности, в частности уменьшения толщины атмосферного перемешанного слоя. Вместе с тем АПТ может приводить и к увеличению толщины перемешанного слоя, в частности, даже полному исчезновению ночной приземной инверсии.

Для оценки антропогенного потока тепла QF, как правило, используются три основных подхода: прямые измерения тепловых потоков; инвентаризация, которая заключается в суммировании всех потребителей энергии; дистанционные спутниковые измерения потоков теплового излучения с выделением антропогенных потоков по локальным метеорологическим данным. Эти подходы подробно описаны в [7].

В настоящей статье применяется подход, основанный на альтернативной инвентаризации, в рамках которой антропогенный поток тепла QF рассчитывается на основе строительных норм, теплофизических свойств ограждающих конструкций и разности между внутренней и внешней температурами воздуха на основе открытой веб-картографической платформы “OpenStreetMap” (OSM) [14]. Платформа OSM позволяет описать реальную трехмерную застройку района или города, указать теплофизические свойства ограждающих конструкций отдельных зданий (стен, окон, крыш, подвалов), температуры внутреннего и наружного воздуха и, таким образом, рассчитать антропогенный поток тепла QF.

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

В рамках проекта РФФИ № 20-05-00254 “Антропогенное воздействие и обратные связи в атмосферном пограничном слое урбанизированных территорий” разрабатывается модель покрывающего слоя городской атмосферы, предназначенная для исследования действия обратных связей и описания антропогенного воздействия на городской пограничный слой атмосферы в мезомасштабных климатических моделях. Основные уравнения такой модели, основанной на работах [1518] и аналогичных исследованиях, приведены в данной статье.

ГОРОДСКАЯ ЗАСТРОЙКА И ПОКРЫВАЮЩИЙ СЛОЙ АТМОСФЕРЫ

Городская застройка создает внутри пограничного слоя атмосферы так называемый городской покрывающий слой (UCL – urban canopy layer). Городской покрывающий слой (полог застройки) простирается от поверхности до средней высоты застройки. В рамках этого слоя на климатические характеристики особенно сильно влияют тепловые свойства зданий и поверхностей и возникающие из-за геометрии зданий и улиц потоки воздуха.

Рис. 1.

Городская застройка, покрывающий и пограничный слои атмосферы над урбанизированной территорией.

Высота городского покрывающего слоя значительно меньше высоты отдельных зданий, в различных работах для пригородов ее величина составляет менее 10 м, в жилых многоэтажных районах – около 20 м, а в малоэтажном историческом центре города – примерно 10–15 м [1921].

В данной работе городская застройка описывается с помощью трехмерного моделирования зданий на основе открытой веб-картографической платформы “OpenStreetMap” (OSM) [14]. Проект OSM является online-ресурсом, содержащим картографическую информацию о городской топографии, геометрии и материалах зданий, дорогах, зеленых насаждениях и т.п. База данных хранится на сервере и постоянно обновляется и пополняется благодаря сообществу картографов, ГИС-профессионалов и инженеров. Информацию с OSM платформы можно свободно экспортировать в формате файлов, имеющих xml-структуру. В базе данных OSM структурировано хранятся координаты точек, объединение которых позволяет построить контуры зданий, дорог, границ различных площадок, зеленых насаждений и т.д. Все точки рассматриваемой территории имеют свои уникальные идентификационные номера id. Здания определяются атрибутами, описывающими их характеристики. Ключевыми атрибутами являются: адрес, категория, высота, этажность, материалы стен и крыш. Количество атрибутов у разных зданий отличается: не у всех зданий имеется атрибут высоты, чаще указывается этажность, иногда есть информация об арках. Если отсутствует атрибут высоты, то используется атрибут, отвечающий за этажность здания, на основании которого рассчитывается высота здания (предполагая высоту этажа 3 м). В случае отсутствия обоих атрибутов высота здания может задаваться, например, по данным ГИС ЖКХ [22]. Для построения контуров зданий используют их атрибуты и соответствующие id точек. После определения всех атрибутов формируется контур здания и образуется 3-D объект.

Что бы получить 3-D модель застройки региона, на визуальном языке программирования grasshopper написан скрипт, формирующий контуры зданий. При моделировании совокупности строений в смежных зданиях не учитываются примыкающие ограждающие конструкции, т.к. через них не происходит передача тепла.

В качестве примера рассматривается Васильевский остров (В.О.) Санкт-Петербурга, карта которого, полученная по данным OSM, приведена на рис. 2, и город Кудрово, расположенный в Ленинградской области, примыкающий к Санкт-Петербургу район новостроек на востоке и фактически являющийся частью Санкт-Петербургской агломерации.

Рис. 2.

Карта Васильевского острова (В.О) Санкт-Петербурга, построенная по данным OSM.

Для построения 3-D модели застройки рассматриваемая территория разбита на квадратные ячейки площадью 30 × 30 м2. На рис. 3 показана построенная 3-D модель сооружений Васильевского острова. На переднем плане 1-я линия и Большой проспект Васильевского острова, в правой верхней части видны линии границ р. Смоленки, а примыкающая к ней незастроенная территория – Смоленское кладбище, а также приведена 3-D модель застройки г. Кудрово. Этот рисунок иллюстрирует тонкую структуру городской застройки.

Рис. 3.

Трехмерная модель зданий по данным OSM В.О. Санкт-Петербурга (а) и г. Кудрово (б).

МЕТЕОРОЛОГИЧЕСКИЕ ХАРАКТЕРИСТИКИ ОТОПИТЕЛЬНОГО ПЕРИОДА

Поток тепла от зданий и сооружений состоит из потока теплопроводности, потока вентиляции помещений и потока радиационного теплообмена через окна. Для расчета коэффициента теплопередачи ограждающей конструкции (стены, окна, крыши и подвалы) необходимо задать ее параметры, т.е. материал, из которого изготовлена конструкция, ее теплофизические характеристики, толщину стены, степень остекления и т.д. Поскольку определение этих параметров для всех сооружений рассматриваемого региона является достаточно громоздкой задачей, то задаются значения теплофизических параметров для стен, окон, крыш и подвалов, которые применяются ко всем рассматриваемым сооружениям. В качестве значений таких параметров используются величины, определяемые на основе требований по тепловой защите зданий, приведенные в строительных нормах [2327]. Таким образом, рассматривается обратная задача, в которой сопротивление теплопередачи ограждающих конструкций вычисляется по значению “требуемого сопротивления теплопередачи”. В прямой задаче коэффициент теплопередачи ограждающей конструкции определяется ее конструктивными особенностями и материалом, из которого она сделана, и проверяется на соответствие “требуемому сопротивлению теплопередачи” [2325].

Согласно Сводам правил [23, 25] и ГОСТ [24], сопротивление теплопередачи ограждающей конструкции оценивается по показателю градусо-сутки отопительного периода (ГСОП)

(1)
${\text{ГСОП}} = \left( {{{t}_{{{\text{вн}}}}}--{{t}_{{{\text{от}}}}}} \right){{d}_{{{\text{от}}}}},$
где tвн – расчетная температура внутреннего воздуха здания, а tот – средняя температура атмосферного воздуха за время отопительного периода, °С, и dот – его продолжительность, сут.

Значение tот определяется как среднее по выборке среднесуточных температур, соответствующих отопительному периоду. В соответствии с СП 131.13330.2012 [23] для жилых и общественных зданий отопительный период начинается, когда среднесуточные температуры наружного воздуха не превышают 8°С, а для лечебно-профилактических, детских учреждений и домов-интернатов для престарелых не более 10°С.

Величина tвн – расчетная температура помещений, °С, принимаемая при расчете ограждающих конструкций групп зданий в соответствии с Г-ОСТом 30494–2011 [24]. Для жилых зданий, общежитий и гостиниц tвн должна быть в интервале 20–22°С, для лечебно-профилактических, дошкольных образовательных и общеобразовательных организаций, интернатов в интервале 16–21°С. Исходя из этого, в качестве расчетного значения принимается tвн = 20°С.

В Санкт-Петербурге период со средней суточной температурой воздуха, не превышающий 8°С, характеризуется tот = –1.3°С, dот = 213 сут, а со средней суточной температурой воздуха не более 10°С имеем tот = –0.4 °С, dот = 213 сут [23]. Для этих параметров величина ГСОП, согласно (1), будет равна 4537 и 4733 (°С сут) соответственно. Изменения средней температуры отопительного периода, его продолжительности и динамика ГСОП в Санкт-Петербурге за последние три десятилетия представлены на рис. 4–6.

Рис. 4.

Среднегодовая температура (°C) отопительного периода tот (1) в Санкт-Петербурге за последние 30 лет и ее тенденция (2) с угловым коэффициентом 0.42°C/10 лет.

Рис. 5.

Продолжительность (сутки) отопительного периода dот (1) в Санкт-Петербурге за последние 30 лет и его тенденция (2) с угловым коэффициентом 5.54 сут/10 лет.

Рис. 6.

Градусо-сутки (°С сут/год) отопительного сезона (1) в Санкт-Петербурге за последние 30 лет и их тенденция (2) с угловым коэффициентом –210.4°С сут/10 лет.

Угловые коэффициенты, приведенные на рис. 4–5, носят исключительно информационный характер, а параметры уравнения регрессии на рис. 6 статистически значимы. Данные, приведенные на рис. 6, показывают, что скорость убывания значения ГСОП по Санкт-Петербургу за последние 30 лет увеличилась, что является еще одним индикатором потепления климата. За рассматриваемый период наблюдаются тенденции увеличения tот и уменьшения dот. Тренд величины ГСОП обусловлен изменением средней температуры атмосферного воздуха во время отопительного периода tот и его продолжительности dот.

ГСОП используется для потребности в отоплении в городах с централизованными системами отопления, начало и конец работы которых определяются решениями городских администраций. Для оценки полной потребности в отоплении помещений, в том числе и за пределами отопительного периода (когда, например, используются электрообогреватели, тепловые насосы и пр.), в различных климатических условиях используются так называемые градусо-сутки обогрева (heating degree-days – HDD). Этот подход был использован в работах [5, 11], что позволило проанализировать динамику потребности в теплоснабжении для отопления зданий в крупнейших российских городах, включая Санкт-Петербург.

Значение HDD рассчитывается как сумма ежедневных разностей между температурой комфорта внутри помещения (базовой температурой – tк)

(1a)
${\text{HDD}} = \sum\limits_{i = 1}^N {({{t}_{{{\text{к,}}i}}} - {{t}_{{{\text{нар,}}i}}})} ,$
где N – число дней со среднесуточной температурой наружного воздуха tнар ниже температуры, при которой возникает необходимость использовать отопление.

В России, Канаде, Японии и большинстве стран Европы tк принимается равной 18°C, в Великобритании – 15.5°C, а в США – 18.3°C (65°F) (см. [7]). Для климатических условий России величины HDD несколько больше значений ГСОП для тех же городов.

В городах России, Европы и Северной Америки с середины XX в. количество HDD и ГСОП заметно убывает. На рис. 7 и 8 в качестве примера приведены изменения HDD за последние десятилетия в среднем по городам Евросоюза и в Санкт-Петербурге.

Рис. 7.

HDD (°С сут/год) в среднем по городам Евросоюза с середины XX в. по данным, приведенным в [7] (1), его тенденция характеризуется угловым коэффициентом, равным –58.95 за период 1951–2014 (2), –81.50 за период 1967–2014 (3), –88.44 (°С сут/10 лет) за период 1991–2014 гг. (4).

Рис. 8.

HDD (1) в Санкт-Петербурге по [11] и их тенденция (2) с угловым коэффициентом –158.6°С сут/10 лет.

Тренды (тенденции) HDD, приведенные на рис. 7 и 8, статистически значимы, что демонстрирует существенное влияние современного потепления климата на потребность городского хозяйства Санкт-Петербурга в энергии для поддержания необходимого теплового режима помещений. Напомним, что величина HDD зависит от выбранного значения температуры комфорта (в [11] и на рис. 7 и 8 используется tк = 18°С), а значения ГСОП определяются в соответствии со Строительными правилами, в соответствии с которыми tвн = 20°С. Ниже используются значения ГСОП, полученные на основе Строительных правил.

Нормативные значения приведенного сопротивления теплопередаче $R_{i}^{{{\text{тр}}}}$ для разных элементов ограждающих конструкций, соответствующие ГСОП = 4000, равны: 2.8 м2 °С/Вт для стен (i = 1); 3.7 м2 °С/Вт для чердачных покрытий и перекрытий над неотапливаемыми подпольями и подвалами (i = 2); 0.63 м2 °С/Вт для светопрозрачных ограждающих конструкций (i = 3), кроме фонарей. Здесь i – номер категории ограждающей конструкции.

МЕТОДИКА ОЦЕНКИ АНТРОПОГЕННОГО ПОТОКА ТЕПЛА НА ОСНОВЕ ТЕПЛОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК ЗАСТРОЙКИ

Расчетная сетка, покрывающая рассматриваемую территорию, разбивает ее на ячейки и делит здания на отдельные части, по которым считается суммарный поток тепла QF внутри каждой j-й ячейки.

Расчет величины АПТ $Q_{F}^{i}$ от разных частей элементов ограждающих конструкций, попадающих в рассматриваемую j-ю ячейку, осуществляется по формуле:

(2)
$Q_{F}^{i}(j) = {{K}_{i}}({{t}_{{{\text{вн}}}}} - {{t}_{{{\text{нар}}}}}){{S}_{i}}(j),$
где ${{K}_{i}} = {1 \mathord{\left/ {\vphantom {1 {R_{i}^{{{\text{тр}}}}}}} \right. \kern-0em} {R_{i}^{{{\text{тр}}}}}}$ – коэффициент теплопередачи ограждающей конструкции i-й категории, tнар – температура наружного воздуха, °С, Si – суммарная площадь элементов ограждающей конструкции i-й категории в j-й ячейке, м2.

Окончательно, АПТ ${{Q}_{F}}(j)$ от частей ограждающих конструкций, удовлетворяющих нормативным значениям требуемого сопротивления теплопередаче $R_{i}^{{{\text{тр}}}}$ и попадающих в j-ю расчетную ячейку, площадью 30 × 30 м2, определяется соотношением

(3)
$Q_{F}^{'}(j) = \sum\limits_{i = 1}^3 {Q_{F}^{i}} ,$

которое позволяет оценить среднее значение АПТ $Q_{F}^{'}$ для рассматриваемой территории

(4)
$Q_{F}^{'} = \frac{1}{N}\sum\limits_j {Q_{F}^{'}(j),} $
где N – количество ячеек расчетной сетки для изучаемого региона.

В результате транспортировки тепла к зданиям по магистралям, расположенным под землей, определенная его доля теряется в виде тепловых стоков $\Delta {{Q}_{{F,g}}}$. Следовательно, средний АПТ должен быть увеличен на соответствующую величину теплового стока,

(5)
${{Q}_{F}} = Q_{F}^{'} + \Delta {{Q}_{{F,g}}}.$

РАСЧЕТ АНТРОПОГЕННОГО ПОТОКА ТЕПЛА ДЛЯ КОНКРЕТНЫХ РАЙОНОВ САНКТ-ПЕТЕРБУРГА

Ниже приведены результаты моделирования АПТ для двух районов агломерации СПб: В.О., расположенный на западе, и г. Кудрово, примыкающий к СПб район новостроек на востоке. Эти районы выбраны с учетом преобладающей розы ветров, ориентированной с запада на восток. Расчеты проведены для температур наружного воздуха tнар в диапазоне от +8 до –30°С. В табл. 1 приведены основные геометрические параметры застройки.

Таблица 1.  

Геометрические параметры застройки Санкт-Петербурга

  Высота, м Площадь крыш, м2 Площадь стен, м2 Площадь окон, м2 Площадь
территории, м2
Доля крыш в общей площади
Среднее Медиана
В.О. 11.85 10 3 350 000 5 433 000 1 811 000 18 070 000 18.54%
Кудрово 19.01 9 693 095 1 731 975 577 325 4 517 000 15.34%

При составлении табл. 1 предполагается, что окна занимают порядка 25% фасада, а 75% фасада приходится на площадь стен. Заметим, что существенное расхождение по высоте среднего и медианного значений для Кудрово обусловлено большим количеством малоэтажных сооружений на его окраинах и высотной жилой застройкой (более 20 этажей). По современным оценкам в СПб при транспортировке тепла через городские магистрали теряется приблизительно 10% тепла, поэтому рассчитанные по территории средние значения АПТ, обусловленные теплопотерями зданий, увеличиваются на (см. (5)),

(6)
$\Delta {{Q}_{{F,g}}} = {{\beta }_{g}}Q_{F}^{'},$
где ${{\beta }_{g}} = 0.1$.

В табл. 2 даны оценки средних значений АПТ по рассматриваемым территориям для различных значений tнар с учетом тепловых потерь в магистралях (6).

Таблица 2.  

Средний по району АПТ QF (Вт/м2) от ограждающих конструкций при различных значениях tнар

tнар , °С +8 0 –5 –10 –15 –20 –25 –30
В.О. 4.63 7.71 9.64 11.56 13.50 15.42 17.35 19.27
Кудрово 5.34 8.90 11.12 13.34 15.57 17.79 20.01 22.23

Зависимость усредненного по региону АПТ от температуры наружного воздуха может быть описана в виде линейной формулы

(7)
${{Q}_{F}} = a + b{{t}_{{{\text{нар}}}}}.$

Для В.О. a = 7.711 Вт/м2, b = –0.385 Вт/(м 2 °С); для Кудрово a = 8.896 Вт/м2, b = –0.445 Вт/(м2 °С). В табл. 3 приведены оценки средних значений АПТ QF, рассчитанные только по массиву зданий городской застройки с учетом (5) для рассматриваемых территорий при различных значениях tнар (без учета тепловых потерь в магистралях).

Таблица 3.  

Средний по зданиям АПТ QF (Вт/м2) от ограждающих конструкций при различных значениях tнар

tнар, °С +8 0 –5 –10 –15 –20 –25 –30
В.О. 22.71 37.81 47.25 56.69 66.18 75.62 85.06 94.5
Кудрово 31.62 52.74 65.91 79.07 92.2 105.4 118.6 131.7

На рис. 9 и 10 приведены распределения АПТ для В.О. и Кудрово для температур наружного воздуха tнар, равных +8, –5 и –20°С.

Рис. 9.

Распределение АПТ, обусловленное городской застройкой, на В.О. СПб для температур наружного воздуха: а) 8°С, б) –5°С, в) –20°С.

Рис. 10.

Распределение АПТ (Вт/м2), обусловленное городской застройкой, в г. Кудрово для температур наружного воздуха: а) 8°С, б) –5°С, в) –20°С.

Для определения среднего за зимний период АПТ необходимо найти фактическое распределение температуры наружного воздуха. Используя данные ВНИИГМИ-МЦД [28] по среднесуточной температуре воздуха на метеорологической станции 26 063 – Санкт-Петербург за последние тридцать лет, оценена относительная частота дней с рассматриваемыми температурами в течение отопительного периода за интервал 1991–2020 гг., а также ежегодно.

Среднее значение АПТ по Лондону, Нью-Йорку и Москве (в границах 2011 г.) составляет 25.1, 21.1 и 55.9 Вт/м2 соответственно [5]. По оценкам, полученным по модели COSMO-CLM (см. [12]), АТП в разных районах Москвы находится в интервале 27–75 Вт/м2: в окрестности МКАД примерно 27 Вт/м2; на некоторых территориях, примыкающих к Садовому кольцу, в ~70–75 Вт/м2; на остальных территориях ~40–45 Вт/м2. Аналогичные оценки для Санкт-Петербурга дают на побережье Финского залива и на территориях, примыкающих к центру города, ~27–30 Вт/м2, а в центральных районах ~45–55 Вт/м2 [12]. Согласно [10] антропогенные потоки тепла за счет центрального отопления в Москве составляют примерно половину от общей величины АПТ. Аналогично в жилых районах Санкт-Петербурге антропогенные потоки тепла за счет центрального отопления оказываются меньше 15 Вт/м2. Эти оценки согласуются с предлагаемой методикой. Так, средние значения АПТ за отопительный период в интервале 1991–2020 гг., усредненные по территории, равны для В.О. и Кудрово равны соответственно 7.52 и 8.68 Вт/м2, а средние значения АПТ, усредненные по массиву зданий, расположенных на этих территориях, соответственно равны 36.88 и 51.40 Вт/м2.

Одна из причин таких значений АПТ на В.О. может быть связана с тем, что основным материалом стен в исторической застройке В.О. является кирпичная кладка, состоящая из 4 кирпичей для цокольного этажа или 2.5 кирпичей на уровне карниза здания [29], оказывается ниже современных требований к сопротивлению теплопередачи. Так, например, кирпичная кладка, состоящая из 3 кирпичей, обеспечивает сопротивлению теплопередачи, рассчитанное по СП 50.13330.2012 и равное 1.67 м2 °С/Вт, что на 60% меньше современных требований. На рис. 11 представлены среднегодовые значения АПТ, усредненные по территории (рис. 11а) и по массиву застройки (рис. 11б) для В.О. СПб и Кудрово за период 1991–2020 гг.

Рис. 11.

Среднегодовые значения АПТ: а) по территории и б) по массиву застройки за период 1991–2020 гг. для (1) В.О. СПб и (2) Кудрово; тенденция АПТ за период 2002–2020 гг. по территории (3, а) В.О. СПб и (4, а) Кудрово описывается угловым коэффициентом 0.44 и ‒0.51, а по массиву застройки ‒2.17 и –3.02 (Вт/м2)/10 лет соответственно.

Для более корректных оценок необходимо учитывать год постройки зданий и использовать соответствующие нормативы, которые за последние десятилетия существенно изменились. Такой подход описан, например, в [30], и его планируется использовать в дальнейшей работе.

ЭНЕРГОБАЛАНСОВАЯ МОДЕЛЬ ГОРОДСКОГО ПОКРЫВАЮЩЕГО СЛОЯ

Рассмотрим энергобалансовую модель городского покрывающего слоя (UCL – urban canopy layer). Уравнение притока тепла для случая пренебрежения горизонтальной адвекцией тепла внутри UCL имеет вид

(8)
$\rho {{c}_{p}}\frac{{\partial T(z)}}{{\partial t}} = {\text{div}}(Q{\kern 1pt} *(z) + {{Q}_{H}}(z) + {{Q}_{E}}(z)) + {{q}_{F}}(z),$
где Q* – эффективный интегральный радиационный поток, QH и QE – потоки явного и скрытого тепла, qF – антропогенный приток тепла, обусловленный городской застройкой, ρ – плотность и cp – теплоемкость воздуха.

Интегрируя (8) по вертикали от поверхности до высоты верхней границы UCL–hCL и пренебрегая горизонтальными изменениями, получаем уравнение притока тепла к UCL

(9)
$\begin{gathered} \rho {{c}_{p}}{{h}_{{CL}}}\frac{{\partial {{T}_{{CL}}}}}{{\partial t}} = (Q_{{CL}}^{*} - Q_{0}^{*}) + ({{Q}_{{H,CL}}} - {{Q}_{{H,0}}}) + \\ + \,\,({{Q}_{{E,CL}}} - {{Q}_{{E,0}}}) + Q_{F}^{'}, \\ \end{gathered} $
где TCL – температура UCL; $Q_{F}^{'} = \int_0^{{{h}_{{CL}}}} {{{q}_{F}}(z)dz} $ – антропогенный поток тепла к UCL; индексами CL и 0 обозначены потоки на верхней и нижней границах UCL.

Граничное условие на постилающей поверхности получим, интегрируя (8) по толщине деятельного слоя,

(10)
${{c}_{g}}{{h}_{g}}\frac{{\partial {{T}_{g}}}}{{\partial t}} = Q_{0}^{*} + {{Q}_{{H,0}}} + {{Q}_{{E,0}}} + \Delta {{Q}_{{F,g}}},$
где cg, hg и Tg – теплоемкость, толщина и температура деятельного слоя подстилающей поверхности; $\Delta {{Q}_{{F,g}}}$ – АПТ, выделяющийся в поверхностном деятельном.

Заменяя температуру деятельного слоя Tg температурой поверхности T0 и объединяя (9) и (10), приходим к уравнению, определяющему эволюцию TCL.

(11)
$\begin{gathered} \rho {{c}_{p}}{{h}_{{CL}}}\frac{{\partial {{T}_{{CL}}}}}{{\partial t}} = Q_{{CL}}^{*} + {{Q}_{{H,CL}}} + {{Q}_{{E,CL}}} - \\ - \,\,{{c}_{g}}{{h}_{g}}\frac{{\partial {{T}_{0}}}}{{\partial t}} + Q_{F}^{'} + \Delta {{Q}_{{F,g}}}, \\ \end{gathered} $
где предполагается, что $\Delta {{Q}_{{F,g}}} = {{\beta }_{g}}Q_{F}^{'}$, βg – коэффициент, характеризующий потери на подземных коммуникациях.

Верхняя граница UCL состоит из поверхностей крыш и некоторой огибающей пространство между зданиями (более мелкие детали не рассматриваются), поэтому QF можно детализировать, явно вводя антропогенный поток тепла $\Delta {{Q}_{{F,CL}}} = {{\beta }_{{CL}}}Q_{F}^{'}$, проходящий через крыши в примыкающий приземный слой атмосферы, где βCL – коэффициент пропорциональности. Тогда

(12)
$Q_{F}^{'} = (1 - {{\beta }_{{CL}}})Q_{F}^{'} + {{\beta }_{{CL}}}Q_{F}^{'}.$

Здесь подразумевается, что антропогенный поток тепла с крыш сразу попадает в примыкающий приземный слой атмосферы. Коэффициент ${{\beta }_{{CL}}} < 1$ показывает часть АПТ, которая попадает в окружающую среду через крышу, непосредственно формируя вертикальный поток тепла на верхней границе UCL, а $(1 - {{\beta }_{{CL}}})$ определяет АПТ от боковых ограждающих конструкций непосредственно внутри острова тепла.

Таким образом, уравнение притока тепла к UCL можно представить как

(13)
$\begin{gathered} \rho {{c}_{p}}{{h}_{{CL}}}\frac{{\partial {{T}_{{CL}}}}}{{\partial t}} = Q_{{CL}}^{*} + {{Q}_{{H,CL}}} + {{Q}_{{E,CL}}} - \\ - \,\,{{c}_{g}}\frac{{\partial {{T}_{0}}}}{{\partial t}} + (1 - {{\beta }_{{CL}}})Q_{F}^{'} + {{\beta }_{{CL}}}Q_{F}^{'} + {{\beta }_{g}}Q_{F}^{'}, \\ \end{gathered} $
где полное значение АПТ определяется соотношением

(14)
${{Q}_{F}} = Q_{F}^{'} + \Delta Q_{{F,g}}^{'}.$

Коэффициент ${{\beta }_{{CL}}}$ можно оценить по табл. 2 при помощи формулы

(15)
${{\beta }_{{CL}}} = \frac{{{{K}_{2}}{{S}_{2}}}}{{{{K}_{1}}{{S}_{1}} + 2{{K}_{2}}{{S}_{2}} + {{K}_{3}}{{S}_{3}}}},$
где Ki – коэффициент теплопередачи i-й категории ограждающей конструкции, Si – их суммарные площади, i = 1 соответствует стенам, 2 – крышам и подвалам, 3 – окнам. Для В.О. СПб и Кудрово ${{\beta }_{{CL}}}$ равен 0.14 и 0.10 соответственно.

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

Работа выполнена при поддержке проекта РФФИ № 20-05-00254. Авторы благодарят рецензента за конструктивные замечания.

Список литературы

  1. Allen L., Lindberg F., Grimmond C.S.B. Global city scale urban anthropogenic heat flux: model and variability // Int. J. Climatol. 2011. V. 31. P. 1990–2005.

  2. Oke T.R., Mills G., Christen A., Voogt J.A. Urban Climates. Cambridge Univ. Press, 2017. 510 p.

  3. Кислов А.В. и др. Климат Москвы в условиях глобального потепления. Изд-во МГУ, 2017. 320 с.

  4. Гаврилова Ю.В., Смышляев С.П., Махура А.Г. Модельное исследование влияния крупных городов на региональные аспекты погоды // Ученые записки РГМУ. 2010. № 15. С. 107–116.

  5. Гинзбург А.С., Белова И.Н., Расплетина Н.В. Антропогенные потоки тепла в городских агломерациях // ДАН. 2011. Т. 439. № 2. С. 256–259.

  6. Александров Г.Г., Белова И.Н., Гинзбург А.С. Антропогенные потоки тепла в столичных агломерациях России и Китая // ДАН. 2014. Т. 457. № 1. С. 101–106.

  7. Гинзбург А.С., Решетарь О.А., Белова И.Н. Влияние климатических факторов на энергопотребление в отопительный сезон // Теплоэнергетика. 2016. № 9. С. 20–27.

  8. Varentsov M., Konstantinov P., Baklanov A., Esau I., Miles V., Davy R. Anthropogenic and natural drivers of a strong winter urban heat island in a typical Arctic city // Atmos. Chem. Phys. 2018. V. 18. P. 17573–1758.

  9. Ginzburg A.S., Dokukin S.A. Numerical modeling of anthropogenic heat flux impact on air temperature in Moscow in wintertime // IOP Conf. Series: Earth and Environmental Science. 2018. V. 211. https://doi.org/10.1088/1755-1315/211/1/012019

  10. Alexandrov G.G., Ginzburg A.S. Anthropogenic impact of Moscow district heating system on urban environment // Energy Procedia. 2018. V. 149. P. 161–169.

  11. Belova I.N., Ginzburg A.S., Krivenok L.A. Heating seasons length and degree days trends in Russian cities during last half century // Energy Procedia. 2018. V. 149. P. 373–379.

  12. Dokukin S.A., Ginzburg A.S. The influence of anthropogenic heat fluxes on the temperature and wind regimes of the Moscow and St. Petersburg regions // IOP Conf. Series: Earth and Environmental Science. 2020. V. 606. https://doi.org/10.1088/1755-1315/606/1/012010 (2020) 012010

  13. Демченко П.Ф., Голицын Г.С., Гинзбург А.С., Вельтищев Н.Н. Оценка суточного цикла парникового эффекта СО2 по одномерным моделям вертикальной структуры атмосферы // Изв. АН. Физика атмосферы и океана. 1994. Т. 30. № 5. С. 595–600.

  14. OpenStreetMap. – URL: https://www.openstreetmap.org (дата обращения: 29.03.2021).

  15. Salamanca F, Krpo A, Martilli A, Clappier A. A new building energy model coupled with an urban canopy parameterization for urban climate simulations – Part I. Formulation, verification, and a sensitive analysis of the model // Theor. Appl. Climatol. 2010. V. 99. P. 331–344.

  16. Bueno B., Roth M., Norford L., Li R. Computationally efficient prediction of canopy level urban air temperature at the neighborhood scale // Urban Climate. 2014. V. 9. P. 35–53.

  17. Kanda M., Kawai T., Kanega M., Moriwaki R., Narita K., Hagishima A. A simple energy balance model for regular building array // Boundary-Layer Meteorol. 2005. V. 116. P. 423-443.

  18. Ma S., Pitman A., Hart M., Evans J.P., Haghdad N., MacGill I. The impact of an urban canopy and anthropogenic heat fluxes on Sydney’s climate // Int. J. Climatol. 2017. V. 37. P. 255–270.

  19. Xiaodong He, Yuhuan Li, Xinran Wang, Liang Chen, Bu Yu, Yizhou Zhang, Shiguang Miao High-resolution dataset of urban canopy parameters for Beijing and its application to the integrated WRF/Urban modelling system // J. Cleaner Production. 2019. V. 208 P. 373–383.

  20. Balogh M., Kristof G. Fine scale simulation of turbulent flows in urban canopy layers // Quarterly J. of Hungarian Met. Service. 2010. V. 1. № 1–2. P. 135–148.

  21. Huq P., Carrillo A., White L.A., Redondo J., Dharmavaram S., Hanna S.R. The shear layer above and in urban canopies // J. Applied Meteorology and Climatology. 2007. V. 46. P. 368–376. https://doi.org/10.1175/JAM2469.1

  22. Государственная информационная система жилищно-коммунального хозяйства (ГИС ЖКХ). – URL: https://dom.gosuslugi.ru (дата обращения: 29.03.2021).

  23. Свод правил. Строительная климатология (СП 131.13330.2012). – URL: http://docs.cntd.ru/ document/1200095546 (дата обращения: 29.03.2021).

  24. Межгосударственный стандарт. Здания жилые и общественные. Параметры микроклимата в помещениях (ГОСТ 30494-2011). – URL: http://docs.cntd.ru/ document/1200095053 (дата обращения: 29.03.2021).

  25. Свод правил. Тепловая зашита зданий (СП 50.13330.2012). – URL: http://docs.cntd.ru/ document/1200095525 (дата обращения: 29.03.2021).

  26. Малявина Е.Г. Теплопотери здания: справочное пособие. М.: АВОК-ПРЕСС, 2007. 144 с.

  27. Богословский В.Н. Строительная теплофизика (теплофизические основы отопления, вентиляции и кондиционирования воздуха). М.: Высш. школа, 1982. 415 с.

  28. Единый государственный фонд данных ВНИИГМИ-МЦД. – URL: http://meteo.ru/data (дата обращения: 29.03.2021).

  29. Головина С.Г. Архитектурно-конструктивные особенности жилых зданий Санкт-Петербурга второй половины XVIII века // Градостроительство и архитектура. 2020. Т. 10. № 2. С. 71–77. https://doi.org/10.17673/Vestnik.2020.02.10

  30. Andrić I., Fournier J., Lacarrière B., Le Corre O., Ferrão P. The impact of global warming and building renovation measures on district heating system techno-economic parameters // Energy. 2018. V. 150. P. 926–937.

Дополнительные материалы отсутствуют.