УДК 621.18:533.36
СОВЕРШЕНСТВОВАНИЕ СИСТЕМЫ ТЕПЛОВОЗДУХОСНАБЖЕНИЯ ГЛАВНОГО КОРПУСА ТЭС НА ОСНОВЕ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ
В.В. БУХМИРОВ, д-р техн. наук, Д.В. РАКУТИНА, канд. техн. наук, А.Ю. ГИЛЬМУТДИНОВ, асп.
Разработана математическая модель тепловоздушного режима главного корпуса Костромской ГРЭС с блоками 300 МВт. Приведено описание результатов вычислительного эксперимента. Предложены рекомендации по совершенствованию системы тепловоздухоснабжения главного корпуса Костромской ГРЭС.
Ключевые слова: температурное поле, система дифференциальных уравнений, математическая модель, уровень инфильтрации, повышение энергоэффективности.
IMPROVING HEAT AND AIR SUPPLY SYSTEM OF MAIN BUILDING AT STEAM POWER STATION ON BASIS OF MATHEMATIC SIMULATION
V.V. BUKHMIROV, Doctor of Engineering, D.V. RAKUTINA, Candidate of Engineering,
А.Уы. GILMUTDINOV, Post Graduate Student
The mathematical model of heat and air mode in the main building at Kostroma DPS with 300 MW power units is developed. The authors present the description of numerical experiment results. Recommendations for
improving the heat and air supply system of the main b
Key words: temperature field, differential equation energy efficiency.
Система тепловоздухоснабжения главного корпуса ТЭС представляет собой комплексный процесс аэродинамики воздушных потоков и теплообмена между тепловыделяющим оборудованием, приборами системы теплоснабжения, воздушной средой и внутренними ограждениями здания. Температурное поле внутри главного корпуса также зависит от температуры окружающей среды и от количества воздуха, забираемого из помещения на технологический процесс.
Математическое моделирование тепловоздушного режима главного корпуса ТЭС основано на решении системы дифференциальных уравнений движения (Навье-Стокса) и сохранения энергии текучей среды (Фурье-Кирхгофа). Система дифференциальных уравнений конвективного теплообмена воздушной среды дополняется уравнением неразрывности.
Дифференциальное уравнение переноса субстанции имеет вид д -
—(рФ) + div(p) = div(fФ gradФ) + Бф, (1)
dt
где р - плотность среды, кг/м3; W - вектор скорости, м/с; Гф - коэффициент переноса переменной Ф; Бф -источниковый член для переменной Ф.
Решение системы уравнений (1) с соответствующими условиями однозначности может быть выполнено на основе численных методов при помощи современных САЕ-систем (computer-aided engineering), к которым относятся такие программные продукты как Ansys, FlowVision, Star-CD, Fluent, Flow3D, Phoenics и др. Ниже предложена математическая модель тепловоздушного режима главного корпуса Костромской ГРЭС с блоками 300 МВт. Математическое моделирование выполнено в среде Phoenics.
Модель главного корпуса состоит из турбогенераторного, деаэраторного и парогенераторного отделений (рис. 1). Разделительная стенка между турбинным и парогенераторным отделениями отсутствует, однако система трубопроводов препятствует прямому перетоку воздуха в парогенераторное отделение. Для учета этого в математической модели между турбинным и парогенераторным отделениями установлены две
>uilding at Kostroma DPS are offered.
system, mathematical model, infiltration level, increasing the
виртуальные перегородки с условной проницаемостью 60 и 40 % соответственно.
Рис. 1. Трехмерная геометрическая модель главного корпуса в ПВК РЬюепюБ (1-10 - точки измерения температуры)
Математическая модель разработана для условий холодного времени года. Приток воздуха в помещение осуществляется через калориферы, встроенные в стены парогенераторного отделения и работающие «на просос», а также за счет инфильтрации через фрамуги световых проемов турбинного отделения. Удаление воздуха из помещения главного корпуса производится через воздухозаборные отверстия дутьевых вентиляторов.
Тепловыделения от оборудования и коэффициенты теплопередачи для наружных ограждений здания были заданы по проекту. Тепловую нагрузку приборов системы теплоснабжения главного корпуса принимали в соответствии с температурным графиком, утвержденным на станции. Расход воздуха, поступающего в главный корпус через калориферы парогенераторного отделения, задавали по данным измерений [1].
Размеры расчетной области: длина - 93 м, ширина - 21 м, высота - 55 м. Решение находили на адаптивной сетке с ее сгущением в местах наибольших градиентов температур. Количество узлов разностной сет-
ки составило 90593. Для решения была применена 1_УЕ1_-модель турбулентности, которая является разновидностью модели пути смешения Прандтля и предназначена для расчета течения жидкости в пространстве, загроможденном различными препятствиями [2].
Достоверность результатов расчета на математической модели проверена путем сопоставления с данными экспериментального исследования микроклимата в главном корпусе Костромской ГРЭС с блоками 300 МВт [1]. Для сравнения использованы средние значения температуры, измеренные в разных точках главного корпуса. Сопоставление результатов вычислительного эксперимента и опытных данных показало, что максимальная относительная погрешность расчета температуры не превышает 10 % (рис. 2).
Результаты расчетов на математической модели представлены в виде регрессионной зависимости температуры воздуха внутри помещения Твн от количества воздуха, забираемого дутьевыми вентиляторами из котельного отделения вдвпом, нагрузки системы теплоснабжения Отс и температуры наружного воздуха Тнв. Уравнения регрессии получены методом полного факторного эксперимента для десяти точек, которые характеризуют температуру в каждом условно выделен-
а)
Рис. 2. Расчетные и экспериментальные значения температуры отделение; ■ - экспериментальные данные; — - результаты
ном объеме помещения. Расположение точек внутри главного корпуса показано на рис. 1.
Расчетное исследование показало, что минимальная температура воздуха наблюдается в точке 1, что объясняется высокой инфильтрацией через оконные проемы турбинного отделения (рис. 3). Максимальная температура воздуха наблюдается перед воздухозаборными отверстиями дутьевых вентиляторов в точке 10. Наибольшее влияние нагрузка приборов системы теплоснабжения оказывает на температуру в турбинном отделении: при изменении Отс от 100 до 50 % температура в точках 1-5 (турбинное отделение) уменьшается на ~ 3-4 0С, а в точках 6-10 (котельное отделение) меняется незначительно (~ 1 0С).
В холодный период года фрамуги оконных проемов турбогенераторного отделения закрыты, и поэтому воздух поступает за счет инфильтрации, которая зависит от проницаемости 1. Проницаемость
оконных проемов, или удельная негерметичность 1, равна отношению суммарной площади неплотностей оконных проемов к суммарной площади поверхности.
б)
при ТНв = -16 0С: а — турбинное отделение; б — парогенераторное расчета
Твн, °С 30 25 20 15
10
10
30
точки измерения
а)
25
20
15
10
точки измерения
б)
Рис. 3. Температура воздуха внутри главного корпуса (Твн) при разной проницаемости оконных проемов турбинного отделения (1) (ввпом = 100 %, работают 4 блока из 8): а - Тнв = 0 0С; б - Тнв = - 30 0С; □ - 1 = 6,65 %; □ - 1 = 3 %
Адекватность математической модели тепловоздушного режима главного корпуса Костромской
ГРЭС подтверждается при 1 = 6,65 %, что соответствует уровню инфильтрации ~ 57 % при номинальной нагрузке блоков. Экспериментальное исследование показало, что инфильтрация в главном корпусе составляет 40-70 % от общего воздухообмена. Таким образом, значение инфильтрации, полученное при расчете на математической модели, полностью совпадает с результатами эксперимента.
На рис. 3 показано изменение температуры в главном корпусе при уменьшении проницаемости оконных проемов турбинного отделения с 6,65 до 3 %. Анализ диаграмм показывает, что при умень-
шении проницаемости до 3 % увеличение температуры воздуха в турбинном отделении составляет от 2 до 10 С, а в котельном отделении - от 0 до 3 0С при температуре наружного воздуха Тнв = 0 С. При Тнв = -30 С температура в турбинном отделении повышается на 17-25 0С, а в котельном отделении на 1-13 0С. Таким образом, для обеспечения допустимой температуры внутри главного корпуса тепловая нагрузка приборов системы теплоснабжения может быть значительно уменьшена за счет уменьшения инфильтрации.
На рис. 4 представлена зависимость
ОдвП°м = 1 (Твн, Отс) при двух значениях проницаемости оконных проемов турбинного отделения 1 .
5
0
5
5
0
5 10 15 20 25 Твн, 0С
Рис. 4. Зависимость Оде"°м = 1 (Твн, Отс) при Тнв = -30 0С и проницаемости оконных проемов турбинного отделения 6,65 % и 3 %: --------Отс = 100 %;-----------Ото = 90 %; --------Отс = 80 %; ^- Отс = 70 %; — — - Отс = 60 %; — — — - Отс = 50 %;
□ - 1 = 3 %, А - f = 6,65 %
Анализ графика показывает, что при Тне = -30 0С и
7 = 6,65 % количество дутьевого воздуха, которое можно забирать из помещения главного корпуса, не должно превышать 50 % при соблюдении требуемой СанПиН [3] температуры воздуха в рабочих зонах котельного, турбинного и деаэраторного отделений. Любое увеличение ЄдвП°м приведет к понижению температуры в рабочих зонах ниже допустимого значения Твн = 18 0С (заштрихованная область). Уменьшение проницаемости оконных проемов до 3 % позволит увеличить забор воздуха из помещения главного корпуса до 100 %, а также уменьшить нагрузку приборов системы теплоснабжения на 40 %.
Заключение
Разработанная математическая модель тепловоздушного режима главного корпуса ТЭС с блоками 300 МВт позволила реализовать принципиально новый подход к повышению энергоэффективности производственного помещения с тепловыделениями, основан-
ный на математическом моделировании при помощи современных САЕ-систем.
В результате проведенного исследования различных режимов функционирования системы тепло-воздухообмена главного корпуса ТЭС выявлен значительный резерв (до 40 %) уменьшения тепловой нагрузки приборов системы теплоснабжения главного корпуса за счет уплотнения оконных проемов турбинного отделения при соблюдении требуемых СанПиН температур в рабочих зонах.
Список литературы
1. Бухмиров В.В., Ракутина Д.В., Гильмутдинов А.Ю.
Экспериментальное исследование системы аэрации главного корпуса Костромской ГРЭС // Вестник ИГЭУ. - 2007. - Вып. 1. -С. 14-18.
2. Polis info system. - London: CHAM Ltd, 2004.
3. СанПиН 2.2.4.548-96 Гигиенические требования к микроклимату производственных помещений.
Бухмиров Вячеслав Викторович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор, зав. кафедрой теоретических основ теплотехники, телефон (4932) 26-97-78, e-mail: buhmirov@tot.ispu.ru
Ракутина Дарья Валерьевна,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры теоретических основ теплотехники, телефон (4932) 26-97-76.
Гильмутдинов Алексей Юрьевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», аспирант кафедры теоретических основ теплотехники, е-таН: awr@dsn.ru