УДК 621.311.22
ПРИКЛАДНОЙ ПРОГРАММНЫЙ КОМПЛЕКС ДЛЯ ПРОЕКТИРОВАНИЯ, ОРГАНИЗАЦИИ ЭКСПЛУАТАЦИОННОГО КОНТРОЛЯ И НАЛАДКИ АТМОСФЕРНЫХ ДЕАЭРАТОРОВ
ЛЕДУХОВСКИЙ Г.В., канд. техн. наук, КОРОТКОВ А.А., НЕНАЕЗДНИКОВ А.Ю., инженеры
Приведено описание прикладного программного комплекса, обеспечивающего расчет параметров работы деаэратора атмосферного давления при известных конструктивных характеристиках струйных и барботажных устройств. Показана возможность и приведены примеры использования разработанного программного комплекса для проектирования, организации эксплуатационного контроля технического состояния и наладки атмосферных деаэраторов.
Ключевые слова: термический деаэратор, десорбция кислорода, десорбция диоксида углерода, программный комплекс, моделирование процессов деаэрации, струйный отсек, барботажное устройство.
АPPUCATЮN SUITE FOR DESIGNING, ORGANIZATION FIELD INSPECTION AND SETTING UP ATMOSPHERIC DEAERATORS
G.V. LEDUHOVSKY, Candidates of Engineering, A.A. KOROTKOV, Engineer, A.J. NENAEZDNIKOV, Engineer
The description of application suite providing calculation of atmosphere pressure deaerator's operational characteristics with certain design characteristics of bubbling and stream equipment is shown. The possibility and examples of developed bundled software's usage for designing, organizing field inspection of technical condition and setting up of atmosphere pressure deaerators are shown.
Keywords: steam deaerator, oxygen desorption, carbon dioxide desorption, bundled software, deaeration process simulation, jet compartment, bubbler.
Особенности конструктивного исполнения деаэраторов различных моделей существенно влияют на их эксплуатационные характеристики. Поэтому в настоящее время технологически оптимальные режимы работы таких деаэраторов устанавливаются исключительно в ходе их натурных испытаний. Проведение последних сопряжено с существенными затратами времени, трудовых и энергетических ресурсов. В связи с этим разработка компьютерных моделей деаэраторов, описывающих протекающие в них технологические процессы с требуемой точностью и прецизионностью, является актуальной задачей.
Для разработки математической модели деаэратора в целом необходимы модели отдельных струйных и барботажных элементов, включенных в технологическую схему деаэратора. Количество таких отдельных элементов
в деаэраторах современных конструкций может достигать десяти, а схема их соединения по потокам воды и греющего пара должна предусматривать рециркуляцию, объединение и разделение потоков (рис. 1). Опубликован ряд альтернативных математических моделей струйных и барботажных элементов [1-3], применение которых при одних и тех же исходных данных приводит к существенно различающимся результатам. По этой причине использованию тех или иных моделей предшествует проверка их адекватности на основе результатов натурных испытаний деаэраторов различных конструкций. Для ряда элементов или диапазонов работы этих эле-
ментов опубликованные математические модели отсутствуют.
11 4
Рис. 1. Пример деаэрационной колонки деаэратора современной конструкции (ДА-300м): 1 и 2 - штуцеры подвода исходной воды; 3 - смесительный колодец; 4 - верхний струйный отсек; 5 - нижний струйный отсек; 6 - барботажный колодец; 7 - барботаж-ный коллектор деаэрационной колонки; 8 - непровальный барботажный лист; 9 - опускные трубы; 10 - гидрозатвор барботажного колодца; 11 - пароперепускные трубы верхней тарелки; 12 - штуцер отвода выпара; стрелками показаны направления движения пара
Нами предложены критериальные уравнения, обеспечивающие расчет процессов нагрева воды и десорбции из нее кислорода в струйных отсеках деаэраторов при струйнокапельном гидродинамическом режиме работы этих отсеков [4]:
19(^2)= ^"^Н^0''83 * (1)
х рг-0,054 рг-0,832 к-0,290
lgI — 1 = 0,89101 — 1 Lap0264Sc-1K-0882 , (2)
H
0,964
С2
. p„.Wn2d с 2 v ж
где Lap = ж п—; Fr =—; Рг = -ж-;
gd
К =------r-----; Sc =^ж;
сж (t2 - 1 °ж
Wn =
W вх - W вых
Wn-------------Wn— при Wпвх/Wпвых > 1,7,
W вх 2,3 lg W—
Wneb'x
W = WnBX + Wneblx „вхп„вых.
при WpT/WpT1* < 1,7,
D вхn вх ex _ Dn_Un_
Wn =
a
W вых = Dn vvn
вых и вых ип
a
4ввх
Wж =~^ж-
nd прж
, L = 3Wж
Рвжх d 3
где t1 и t2 - температуры воды соответственно на входе и выходе струйного отсека, оС; tH -температура насыщения пара в отсеке, оС; С1 и С2 - концентрации растворенного в воде кислорода соответственно на входе и выходе струйного отсека, мкг/дм3 (проба при 25 оС); Lap, Fr, Pr, K и Sc - критерии соответственно Лапласа, Фруда, Прандтля, Кутателадзе и Шмидта; Н - высота отсека от нижней поверхности струеобразующей тарелки до борта нижней тарелки, м; d - диаметр отверстий струеобразующей тарелки, м; рж - средняя плотность воды, кг/м ; Wn - средняя скорость парового потока в отсеке, м/с; стж - средний коэффициент поверхностного натяжения воды, Н/м; W»,; - средняя скорость истечения воды из отверстий струеобразующей тарелки, м/с; g -ускорение свободного падения, м/с2; \ж - коэффициент кинематической вязкости воды, м /с; аж - коэффициент температуропроводности воды, м2/с; r - скрытая теплота парообразования, кДж/кг; сж - истинная изобарная теплоемкость воды, кДж/(кгоС); йж - коэффициент молекулярной диффузии кислорода в воде, м2/с; <Зжвх - массовый расход воды на входе в отсек, кг/с; n - число отверстий струеобразующей тарелки, шт.; Wnex и Wneblx - скорости парового потока соответственно во входном и
вх вых
выходном сечениях отсека, м/с; ип и ип -удельные объемы пара соответственно во входном и выходном сечениях отсека, м3/кг; авх, авых - площади проходного сечения соответственно на входе и выходе пара из отсека, рассчитанные через высоту отсека H, за вычетом площади, занимаемой условными струями воды постоянного по высоте диаметра d, м2; L - средняя высота зоны чисто струйного режима течения воды в отсеке, м.
После тщательной проверки адекватности модели, аналогичные (1), (2), могут быть использованы для составления модели деаэратора в целом. В разработанном нами программном комплексе использован следующий подход. На первом этапе выполняется анализ конструктивной схемы деаэратора с выделением отдельных струйных и барботажных элементов, элементов потерь давления, объединения и разделения потоков. Цель состоит в составлении расчетной схемы деаэратора. Например, для деаэратора ДА-300м, конструктивная схема колонки которого показана на рис. 1, оборудованного барботаж-ным устройством деаэраторного бака, расчетная схема включает 15 элементов (рис. 2).
Полученная таким образом расчетная схема вносится в программный комплекс в качестве исходных данных. Для этого необходимо указать тип каждого элемента (струйный кольцевой, струйный сегментный, барботажный лист, деление потока и т.п.), его конструктивные характеристики (по запросу) и задать направления движения потоков теплоносителей в виде матрицы коммутации.
По заданной расчетной схеме автоматически синтезируется математическая модель деаэратора в целом. Для этого используются процедура выбора алгоритма расчета того или иного элемента из базы и связи алгоритмов расчета отдельных элементов по входным и выходным параметрам потоков теплоносителей. В результате пользователю представляется экранная форма с полями для ввода значений параметров внешних потоков деаэратора (рис. 3). При этом имеется возможность просмотра всех результатов расчета по каждому элементу, а также ошибок, возникших в ходе расчета (рис. 4). Например, алгоритм расчета предусматривает контроль отсутствия перелива воды через борт струеобразующей тарелки, провала воды через отверстия барбо-тажного листа, отсутствия выпара барботажно-го устройства, срабатывания гидрозатворов.
В настоящее время программа предусматривает выполнение только прямого расчета: по значениям параметров исходной воды, направляемой на деаэрацию, греющего пара (основного и барботажного) и давления выпара рассчитываются параметры деаэрированной воды на входе и выходе каждого отдельного элемента. При этом проектный расчет деаэратора проводится следующим образом: за-
ст
Ж
Ж
ст
ж
дается конструктивная схема аппарата, после чего проверяется ее работоспособность и вычисляются характеристики эффективности технологического процесса.
ХОВ
\
2
14
Ч
8 г 1 1 1 3 1
\
1 7 | 4
і
ІДВ
6 А _ \
5 А
БП
Рис. 2. Расчетная схема деаэратора ДА-300м: 1 и 2 -соответственно первая и вторая по ходу движения воды зоны верхнего струйного отсека; 3 - нижний струйный отсек; 4 - барботажный лист; 5 - объем воды в деаэраторном баке; 6 - надводное пространство деаэраторного бака; 7, 8 - элементы разделения парового потока; 9, 11 - элементы потерь давления пара; 10 - элемент смешения паровых потоков; 12 - элемент разделения потоков воды; 13, 14 - элементы смешения потоков воды; 15 - охладитель выпара деаэратора; ХОВ - основной поток химически очищенной воды, подаваемой в деаэратор; ХОВ на ОВА - поток воды на охладитель выпара деаэратора; В - выпар деаэратора; КВ - конденсат выпара деаэратора за охладителем выпара; ОП и БП - потоки основного и барботажного пара, подаваемого в деаэратор; ДВ - отвод деаэрированной воды; сплошные линии - направление движения потоков воды; пунктирные линии - направление движения потоков пара
Далее приведем примеры использования программного комплекса для целей организации эксплуатационного контроля за техническим состоянием деаэраторов и их режимноналадочных испытаний.
Расчетная программа, разработанная для деаэратора ДСА-300 ст. № 1 деаэрацион-но-питательной установки теплосилового цеха ОАО «Северсталь», позволяет контролировать соответствие работы деаэратора характеристикам, полученным в ходе его режимноналадочных испытаний. Стабилизировав режим работы деаэратора и измерив при двухтрех различных режимах необходимые теплофизические и химические характеристики потоков воды и пара, выполняется расчет значений параметров деаэрированной воды. Сравнение результатов расчета и результатов измерения позволяет оценить соответствие деаэратора его режимной характеристике. Их систематическое расхождение указывает на изменение условий работы деаэратора, связанное с возникновением дефектов внутренних элементов. В настоящее время, к сожалению, не удается диагностировать, какой именно из внутренних элементов поврежден.
Имеется и первый положительный опыт использования описанной программы: с ее помощью в порядке оперативного контроля установлено возникновение дефекта барботажного коллектора деаэраторного бака в виде разгерметизации одного из его фланцев. При этом выявлено систематическое превышение фактического значения остаточной концентрации растворенного в деаэрированной воде кислорода над расчетным значением этого показателя (в среднем более чем на 30 %). Уменьшение температуры исходной воды на входе в деаэратор с 73 до 49 оС подтвердило вывод о наличии дефекта: в баке возникали гидроудары, чего раньше при той же гидравлической нагрузке не наблюдалось.
Использование разработанного программного комплекса существенно сокращает объем режимно-наладочных испытаний. Целью испытаний является получение основной режимной характеристики деаэратора в графическом виде - зависимости остаточной концентрации кислорода в деаэрированной воде от температуры исходной воды и гидравлической нагрузки деаэратора. Такая характеристика строится при постоянных значениях удельного расхода выпара деаэратора, давления пара в колонке, удельного расхода пара на барботаж, уровня воды в деаэраторном баке и других параметров. Ясно, что основная режимная характеристика в таком виде не отражает реальной эффективности деаэратора, поскольку не учитывает изменчивость многих факторов. Получение поправок, учитывающих эту изменчивость, в ходе натурных испытаний деаэратора практически невозможно, поскольку требует проведения большого числа опытов.
1
О Thermal deaerators. Design and Calculate
Файл Сервис База данных О программе
■ I
Щ Новый расчет таты расчета
Открыть Сохранить Сохранить как ...
Экспорт ...
® Выход
ІБарботажньїй пар xj|
G,t/4..............
t, *С..............
Р, кгс/см2.........
02, мкг/дмЗ........
Що, мкг-экв/дмЗ.. Х25, мкСм/см.......
3,4
194,5
0,613
0,11
2,0
у ок
G. мЗ/ч 300.4
t +С |59,5
02, мкг/дмЗ 6355
[До, мкг-экв/дмЗ... 1400
Х25. мкСм/см 22S
■У ОК.
Основной пар
У QK
t, *С 1192,5
Р. кгс/см2 10,192
02, мкг/дмЗ F"
[До, мкг-экв/дмЗ.... |0,11
Х25, мкСм/см |2,0
Общие данные
Атмосфера
Р, кПа.............Р' 1
Уровень
h, мм........
Выпар
1902
Р, кгс/см2..........0,190
У ок
Р, кгс/см2 0,192
ХОВ на вкоде в ОВА
t/C 19,5
ХОВ на выходе из ОВА
[/С 52,4
КВ на выходе из ОВА
мЗ/ч 1,684
С 81,1
Рис. 3. Основная экранная форма программного комплекса для расчета атмосферных деаэраторов
0 Thermal deaerators. Design and Calculate
Файл Сервис База данных О программе Исходные данные Результаты расчета |
^Lnjxj
Результаты расчета
Ошибки
-РАСЧЕТ ТЕПЛООБМЕНА: НСО-
Средняя температура насыщення пара в отсеке, оС 102,6755
Коэффициент температуропроводности воды на входе, м2/с Ш94Е-7
Коэффициент температуропроводности воды на выноде, м2/с Ш05Е-7
Средний коэффициент температуропроводности БОДЫ в струйном отсеке, м2/с Ш49Е-7
Удельный объем веды в состоянии насыщения при среднем давлении е струйном отсеке, мЗ/кг 0,00104545
Коэффициент поверхностного натяжения боды б состоянии насыщения при среднем давлении в струйном отсеке, Н/м
0.0583
Критерий Лапласа для воды в состоянии насыщения при среднем давлении в струйном отсеке 114,4940
Критерий Фруда для воды 2,3885 '
Коэффициент температуропроводности воды в состоянии насыщения при среднем давлении в струйном отсеке, м2/с
1,ББ77Е-7
Критерий Прандтля для воды в состоянии насыщения при среднем давлении в струйном отсеке 1,7186
Критерий Кутателадзе для воды в состоянии насыщения при среднем давлении в струйном отсеке 32,7112
Температура воды к конце струйной зоны в нижнем струйном отсеке (Модель: ИГЭУ], оС 102,6279
ПЕРЕГРЕВ основного пара, оС:
8,7В
ПЕРЕГРЕВ барботажного пара, оС:
4,45
ПЕРЕГРЕВ пара на входе в колонку, оС:
6,44
ПРЕДУПРЕЖДЕНИЕ: Теплофизическое состояние пара на вкоде в пароперепускные трубы соответствует состоянию влажного пара ПРЕДУПРЕЖДЕНИЕ: Теплофизическое состояние пара на вкоде в ВСО соответствует состоянию влажного пара
ПРЕДУПРЕЖДЕНИЕ: Пропуск пара через гидрозатвор барботажного колодца
Рис. 4. Экранная форма, содержащая детальные результаты и ошибки расчета
Полный комплекс таких поправок может быть получен расчетным путем с использованием разработанного программного комплекса. Для этого необходимо провести испытания деаэратора при 9-10 режимах. В большинстве случаев этого достаточно, чтобы «настроить» программный комплекс под конкретный деаэратор. После этого все интересующие характеристики определяются расчетным путем. Завершающим этапом является проверка адекватности полученной основной режимной характеристики с комплексом поправок, для чего
проводятся дополнительные опыты при работе деаэратора в заранее выбранных режимах.
Заключение
Разработанный прикладной программный комплекс обеспечивает расчет параметров работы деаэратора атмосферного давления при известных конструктивных характеристиках струйных и барботажных устройств, заданных теплотехнических и химических параметрах качества теплоносителей.
Программный комплекс может быть использован при проектировании струйно-барботажных деаэраторов атмосферного давления для проверки работоспособности конструкции деаэратора и предварительного расчета показателей его эффективности.
Показана возможность использования разработанного программного комплекса для целей эксплуатационного контроля за техническим состоянием атмосферных деаэраторов: по систематическому отклонению расчетных показателей работы деаэратора от опытных значений возможно диагностирование внутренних дефектов деаэратора.
Использование программного комплекса позволяет существенно сократить объем натурных режимно-наладочных испытаний и по-
высить прогностические возможности основной режимной характеристики деаэратора.
Список литературы
1. Оликер И.И. Термическая деаэрация воды в отопительно-производственных котельных и тепловых сетях. - Л.: Стройиздат, 1972.
2. Оликер И.И., Пермяков В.А. Термическая деаэрация воды на тепловых электростанциях. - Л.: Изд-во «Энергия», 1971.
3. Шарапов В.И., Цюра Д.В. Термические деаэраторы. - Ульяновск: Ульян. гос. техн. ун-т, 2003.
4. Ледуховский Г.В. Совершенствование технологии десорбции кислорода в струйно-барботажных деаэраторах атмосферного давления: Дис... канд. техн. наук. -Иваново , 2008.
Ледуховский Григорий Васильевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры тепловых электрических станций, [email protected]
Коротков Александр Александрович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
инженер кафедры тепловых электрических станций,
Ненаездников Александр Юрьевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
инженер кафедры тепловых электрических станций,