Влияние влажности дутьевого воздуха на энергетическую эффективность и надежность жаротрубных водогрейных котлов

В современных условиях эксплуатации котельного оборудования на первый план выходит точность прогнозирования энергетической эффективности при переменных параметрах окружающей среды. Для инженера, эксплуатирующего котел в условиях переменного климата (например, в приморских регионах), важно понимать, как влажность воздуха «съедает» КПД. Одним из наименее изученных, но значимых факторов является переменное влагосодержание атмосферного воздуха (d_в), используемого в качестве окислителя. В данной работе на примере водогрейного котла Viessmann Vitomax 200-LW мощностью 2,1 МВт (тип M241) проведено математическое моделирование процессов теплообмена при изменении d_в в диапазоне от 5 до 25 г/кг сухого воздуха. Остальные исходные данные (температура воздуха 27°C, коэффициент избытка воздуха 1,12, температура топлива 16°C) оставались неизменными. Котёл работает на природном газе с низшей теплотой сгорания Q_p = 38838 кДж/м³ (при нормальных условиях). Теплоноситель – вода с параметрами 70/90°C.

Для реализации исследования была разработана динамическая поверочная модель, базирующаяся на нормативном методе теплового расчета котельных агрегатов, но дополненная алгоритмом итерационного подбора температур. В нормативной документации (например, ГОСТ Р 54825–2011) влажность воздуха при испытаниях принимается равной 10 г/кг, однако реальные условия эксплуатации могут существенно отличаться. В отличие от упрощенных инженерных методик, в расчетную модель были внесены уточняющие коэффициенты: в частности, для определения объема влажного воздуха и продуктов сгорания использован прецизионный коэффициент 0,00161, обеспечивающий корректный баланс влаги в газовом тракте.

Особенностью разработанного программного комплекса является использование VBA-макроса, который реализует метод последовательных приближений для замыкания теплового баланса в каждом элементе котла. Алгоритм в автоматическом режиме производит подбор температур за поверхностями нагрева до достижения невязки теплового баланса менее 2%, что позволяет исключить субъективную погрешность при анализе влияния влажности. Применение метода последовательных приближений позволило минимизировать невязку теплового баланса до уровня 0,24–0,25%, что существенно выше точности стандартных инженерных методик.

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

Результаты моделирования в диапазоне влажности d_в = 5…25 г/кг выявили специфическую динамику распределения тепловых нагрузок по тракту котла Vitomax 200-LW. Наиболее выраженная реакция на изменение влажности зафиксирована в радиационной части (топке) и поворотной камере.

С ростом влажности воздуха концентрация трехатомных газов (водяных паров) в продуктах сгорания увеличивается, что приводит к росту степени черноты газовой среды. Это интенсифицирует радиационный теплообмен, в результате чего температура газов на выходе из топки падает с 1233,24°C (при d_в = 5) до 1211,64°C (при d_в = 25). Суммарное падение температуры после радиационных поверхностей составило порядка 22°C.

Однако дальнейший ход температурного графика в конвективных пакетах демонстрирует эффект самокомпенсации. Несмотря на существенное снижение температуры газов на входе в трубную часть, температура на выходе из конвективных пакетов практически стабилизируется (изменение зафиксировано в 1°C). Этот эффект объясняется компенсационным ростом коэффициента теплоотдачи при увеличении скорости и объема влажных дымовых газов.

Увеличение влагосодержания ведет к росту суммарного объема продуктов сгорания. В условиях фиксированного сечения жаровых и дымогарных труб это неизбежно влечет за собой рост скорости газового потока. Так, в первом конвективном пакете скорость газов возрастает с 24,86 м/с до 25,52 м/с, во втором – с 13,92 до 14,30 м/с (при d_в = 5 и 25 г/кг соответственно). Рост скорости, согласно критериальным уравнениям, интенсифицирует конвективный теплообмен, однако одновременно увеличивает аэродинамическое сопротивление тракта.

Одновременно с ростом сопротивления наблюдается положительный эффект интенсификации теплообмена. Коэффициент конвективной теплоотдачи α_k в первом пучке труб увеличивается с 52,7 до 53,9 Вт/(м²*К). Рост коэффициента теплоотдачи происходит благодаря турбулизации потока (росту числа Рейнольдса) за счет более высокой плотности и скорости влажных газов. Именно этот фактор объясняет стабильность температуры газов на выходе из котла: падение температурного напора на входе в пучок компенсируется возросшей способностью среды передавать тепловую энергию.

Увеличение d_в влечет за собой рост балластной массы продуктов сгорания, что повышает их расчетную скорость в трубном пучке. Согласно критериальным уравнениям конвективного теплообмена, рост скорости потока увеличивает число Рейнольдса и, как следствие, коэффициент теплоотдачи k. Таким образом, интенсификация теплосъема в трубках полностью нивелирует снижение температурного напора на входе, обеспечивая стабильность температуры уходящих газов на уровне 185,71–186,79°C. Этот результат, к слову, полностью соответствует заявленному в паспорте на котел значению в 185°C, что в очередной раз доказывает правильность и точность расчета.

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

Результаты расчетов показывают, что при увеличении влагосодержания воздуха с 5 до 25 г/кг коэффициент полезного действия (КПД) котла снижается с 91,39% до 91,18% (рис. 1). Данное падение эффективности обусловлено ростом потерь тепла с уходящими газами (q₂), так как дополнительная влага, не участвуя в химической реакции горения, требует затрат энергии на нагрев. Это также подтверждается ростом часового расхода топлива со 212,84 до 213,33 м³/ч для поддержания заданной тепловой мощности.

Снижение КПД на 0,21% при изменении влажности на 20 г/кг может показаться незначительным в масштабе одного часа, однако в годовом исчислении для котла мощностью 2,1 МВт это ведет к перерасходу топлива в довольно значительном объеме, что напрямую отобразиться на технико-экономических результатах эксплуатации оборудования.

Рис. 1. Зависимость КПД котла от влагосодержания воздуха

Рост влажности дутьевого воздуха является прямой причиной снижения коэффициента полезного действия котла. Физика процесса обусловлена тем, что водяной пар, поступающий с воздухом, выступает в роли «балласта»: он не участвует в химической реакции горения, но требует затрат теплоты на нагрев до температуры уходящих газов. По результатам расчетов, потеря тепла с уходящими газами q₂ возрастает с 8,04% до 8,24% при увеличении влажности до 25 г/кг. Это влечет за собой необходимость повышения расхода топлива на 0,23% для поддержания заданной тепловой мощности, что подтверждает значимость учета климатического фактора при долгосрочном планировании энергопотребления.

Особое внимание следует уделить изменению температуры точки росы водяных паров (t_р). Согласно полученным данным, рост d_в приводит к повышению t_р с 57,37°C до 59,84°C, показано на (рис. 2). Однако для исследуемого котла Vitomax 200-LW сохраняется значительный запас надежности (более 100°C), что исключает риск низкотемпературной коррозии в исследуемом диапазоне.

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

Аэродинамический расчет зафиксировал рост сопротивления газового тракта на 4,45% (с 546,60 до 570,93 Па) как отображено на (рис. 3). Данный результат расчета также совпадает с заявленным в паспорте на агрегат значением в 550 Па. Этот фактор роста необходимо учитывать при подборе тягодутьевых машин для регионов с высокой влажностью, так как повышенное сопротивление может привести к выходу горелочного устройства за пределы рабочего диапазона по напору воздуха.

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

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