УДК 621.187.11
Расчет концентраций ионных примесей в предельно разбавленных водных растворах типа конденсата и питательной воды паровых котлов (рб > 10 МПа)
А.Б. Ларин, А.Я. Сорокина ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,
г. Иваново, Российская федерация E-mail: [email protected]
Авторское резюме
Состояние вопроса: Расчет концентраций ионных примесей питательной воды энергетических котлов по показателям приборов-анализаторов давно интересует специалистов, однако практическая реализация расчетных методов стала возможной лишь с развитием систем химико-технологического мониторинга. В условиях сверхчистой питательной воды котлов современных энергоблоков наиболее достоверным автоматически измеряемым параметром является удельная электропроводность охлажденных проб. На базе измерений электропроводности прямой и Н-катионированной пробы возможен расчет ряда важных нормируемых параметров. В связи с этим необходима разработка расчетной методики определения рН, концентраций аммиака и хлоридов. Материалы и методы: Исходными данными (материалами) являются нормативные документы, определяющие требования к уровню концентраций примесей, и результаты химического анализа питательной воды различных энергоблоков ТЭС. В качестве метода принят расчетный метод, разработанный на кафедре ХХТЭ ИГЭУ и основанный на теории растворов электролитов. Особенностью метода является использование лишь измерений электропроводности охлажденных проб и оригинальной расчетной методики, применимой в широком диапазоне соотношений анионного состава примесей питательной воды.
Результаты: Разработана методика и алгоритм расчета рН и концентраций аммиака, хлоридов, гидрокарбонатов в питательной воде при изменении соотношения концентраций анионов.
Выводы: Представленная разработка (методика) позволяет по измеренным и приведенным к 25 оС значениям удельной электропроводности (% и %н) охлажденных проб питательной воды энергетических котлов производить косвенные (расчетные) измерения концентраций хлоридов, гидрокарбонатов, аммиака и значения рН и может использоваться для совершенствования систем химического контроля и управления водным режимом энергетических котлов ТЭС. Разработанная методика применима в системах химико-технологического мониторинга энергетических котлов различных типов ТЭС. Практическая значимость исследования состоит в реализации методики для определения рН сверхчистых сред, что подтверждено патентом на изобретение.
Ключевые слова: химический контроль водного режима на ТЭС, измерение электропроводности водных растворов, расчетная методика нормируемых параметров рабочей среды.
Calculation of ionic impurities concentrations in extremely diluted aqueous solutions such as condensate and feed water in steam boilers (pb > 10 MPa)
A.B. Larin, A.Ya. Sorokina, Ivanovo State Power Engineering University, Ivanovo, Russian Federation E-mail: [email protected]
Abstract
Background: The problem of calculating feed water ionic impurity concentrations in power boilers by analyzer devices indicators has long been in the focus of researchers' attention but it only became possible to practically implement the calculation methods when the chemical technology monitoring began to develop. Specific electric conductivity of cooled samples is the most reliable automatically measured parameter of modern power plant boilers operating in ultrapure feed water conditions. It is possible to calculate a number of important rate parameters by measuring direct and cation electric conductivity. All this makes it important to develop a calculation method for determining pH, ammonia and chloride concentrations. Materials and methods: The initial data (materials) are the regulatory documents setting the requirements for impurity concentrations and the results of chemical analysis of feed water of different power generating units of heat power plants. As the research method we used the calculation method developed by the Department of Chemistry and Chemical Technologies in Power Engineering of Ivanovo State Power Engineering University based on electrolyte solution theory. The distinguishing feature of this method is that it only uses measurements of electric conductivity of cooled samples and an original calculation method applicable to a wide range of anion composition ratios of feed water impurities.
Results: We have developed a method and an algorithm of calculating pH and concentrations of ammonia, chlorides, hydrocarbonates in feed water at anion concentration alterations.
Conclusions: The proposed method allows us to make indirect (design) measurements of chloride, hydrocarbonate, ammonia and pH concentrations based on the measured and reduced to 25 оС electric conductivity (% and %н) values of power boiler feed water cooled samples. It can be used to improve systems of water regime chemical monitoring and control in heat power plant power boilers. The developed method can also be used in chemical technology monitoring systems of power boilers of different heat power plant types. Its practical applicability consists in the ability to determine pH of ultra-pure mediums, which is confirmed by a patent of invention.
Key words: chemical monitoring of heat power plant water regime, water solution electric conductivity measurement, calculation method of working medium rate parameters.
DOI: 10.17588/2072-2672.2017.2.013-017
Введение. Расчет концентраций ионных примесей питательной воды барабанных и прямоточных энергетических котлов СВД и СКД по измерениям электропроводности и рН охлажденных проб может выполняться решением системы уравнений, полученных путем преобразования обобщенной математической модели ионных равновесий [1]. Промышленные испытания показали необходимость совершенствования расчетной методики и метода химического контроля в силу следующих причин:
1. В сверхчистой воде с удельной электропроводностью менее 0,3 мкСм/см усиливается неустойчивость в измерениях рН, что заставляет искать другой, отличный от потенциометрического, метод измерения рН.
2. Н-катионированная проба питательной воды представляет собой сложную систему предельно разбавленного раствора смеси электролитов HCl, NaCl, NaHCO3, H2CO3, равновесное состояние которой определяется как составом и концентрациями примесей питательной воды, так и условиями работы Н-катионитовой колонки.
3. Принятые в методике [1] допущения не удовлетворяют точности косвенных измерений с переходом к питательной воде котлов-утилизаторов энергоблоков ПГУ (рб < 7,0 МПа) в связи с недостаточной декарбонизацией питательной воды в атмосферном деаэраторе тепловой схемы блока.
Методы исследования. Предлагаемая методика расчета ионных примесей питательной воды основана лишь на измерениях удельной электропроводности охлажденных прямой и Н-катионированной проб питательной воды энергетических котлов, включая котлы-утилизаторы ПГУ. Некоторые частные решения математической модели ионных равновесий с учетом названных факторов получены с участием авторов и представлены в литературе [2, 3]. Ниже приведены результаты теоретического анализа ионных равновесий в питательной воде и представлен расчетный алгоритм, применимый для широкого круга задач автоматического химического контроля качества водного теплоносителя современных энергоблоков ТЭС и АЭС. Подобный подход давно интересовал специалистов как в России [4, 5], так и за рубежом [6, 7], однако практическая реализация расчетных методов стала возможной лишь с развитием систем химико-технологического мониторинга [2].
Преобразование обобщенной математической модели ионных равновесий [2] определяет концентрацию хлорид-ионов (мкмоль/дм3) в питательной воде в виде выражения
[ОГ ]
1000 •
Хн
426,2 • а (1 + П) +126,5 • (1-а )(1 + n) - n • 31,9
(1)
где а1 = 0,1-0,99 - параметр, характеризующий эффективность удаления катионов натрия на Н-катионитовой колонке*; п = [НСО3-]Н/[СГ]Н = = 0,1-2,0 - отношение концентрации гидрокарбонатов и хлоридов в фильтрате Н-колонки, определяемое условиями ее работы. Идентичность параметров «а» и «п» позволяет объединить их в один с расширением функциональной нагрузки. Т.е. принимаем а1 = 1, но учитываем его влияние на вычисление [СГ] эмпирическим коэффициентом «п». Тогда выражение (1) можно записать в виде 1000-Хн
[Ol- ] = ■
(2)
426,2 - (1 + п) - п - 31,9'
где хн - удельная электропроводность Н-катио-нированной пробы при 25 оС; п = 0,1-2,0 - эмпирический параметр, определяемый условиями работы Н-колонки и зависящий от соотношения гидрокарбонатов и хлоридов.
Лабораторные и промышленные исследования позволяют определить значение параметра «п» в следующих диапазонах:
- для прямоточных котлов СКД п = 0,1-0,5;
- для барабанных котлов (рб > 10 МПа) и котлов-утилизаторов ПГУ (рб > 7 МПа) п = 0,5-1,0;
- для барабанных котлов-утилизаторов ПГУ с рб < 7 МПа и глубоко обессоленной воды п = 0,9-1,5; котлов-утилизаторов ПГУ с рб > 7 МПа п = 0,5-1,0.
Изменение параметра «п» в отмеченных диапазонах отвечает условию оперативного определения концентраций примесей в пределах 10 %-го отклонения расчетных и аналитических измерений. В условиях отсутствия практических методик аналитических измерений хлоридов, гидрокарбонатов, аммиака в сверхчистых средах (хн < 0,3 мкСм/см) предложенный расчетный метод является наиболее приемлемым. Достоверность методики может быть оценена с помощью использования высокочувствительных методов анализа: метода ионной жидкостной хроматографии [8] или других методов накопления (концентрирования) ионных примесей, например, в результате упаривания питательной воды до состоя-
параметр ai определяется из выражения 1000•
[H+ ]н = ai
V + V
При ai = 1 ионный обмен на
Н-колонке идет на 100 %, тогда [H+ ]н = [H+ ^ = [Cl-
ния котловой воды барабанных котлов и испарительных установок [2].
Конкретные значения эмпирического параметра «п» определяются наладкой измерительной системы на действующих котлах и могут быть приняты одинаковыми для типового оборудования. Концентрация хлоридов (мкмоль/дм3) при заданных «л» имеет простое выражение:
[С1-] = 2,15 • хн (при п = 0,1); [а ] = 1,55 • хн (при п = 0,55); [С1-] = 1,22 • хн (при п = 1,0); [а ] = 0,98 • хн (при п = 1,5). Учитывая баланс форм углекислоты до и после Н-катионитовой колонки [2], можно записать выражение концентрации гидрокарбонат-ионов [НС03-] в виде
0,45 + 2,34-Хн
[НСО3 ] = п - [О!- ]
(3)
0,45 + 10-рН+6 ' Для питательной воды различных типов энергетических котлов значения рН рекомендуется поддерживать дозировкой аммиака на уровне 8,0 и выше, кроме нейтрально-кислородного водного режима, почти повсеместно замененного на кислородно-аммиачный. Тогда знаменатель выражения (3) можно представить в виде (0,45+а), где а < 0,01 для рН > 8, что позволяет игнорировать значение «а» в рамках принятых допущений.
Тогда выражение концентрации гидрокарбонат-ионов (мкмоль/дм3) принимает простой вид: [НСОз-] = (1+5,2 хн) п[СГ]. (4)
Подстановкой значений «п» и [С1] получим:
• при п = 0,1 [НСОз-] = 0,215хн+1,118хн2;
• при п = 0,55 [НСОз-] = 0,853хн+4,43хн2;
• при п = 1,0 [НСОз-] = 1,22хн+6,34хн2;
• при п = 1,5 [НСОз-] = 1,47хн+7,64хн2.
Значение рН пробы питательной воды можно определить, рассчитывая концентрацию гидроксильных ионов [ОН-] путем совместного решения уравнений электропроводности (5) и электронейтральности (6) при измеренных значениях удельной электропроводности прямой (х) и Н-катионированной (хн) проб:
1000х = хма+ [Мв+ ] + хМн+ [МН+ ] + v [Н+ ] +
4 Н (5)
+хсГ [С!- ] + хнсо- [НСО- ] + хон- [ОН- ],
[N8+ ] + [Ы Н+ ] + [Н+ ] = [С!- ] + [НСО3 ] + [ОН- ], (6) где [N8+ ], [ЫН+ ] и др. - концентрации катионов
з
и анионов в питательной воде, мкмоль/дм ;
водности при бесконечном разбавлении, Ом-1-см2-моль-1.
Как отмечено выше, при дозировке аммиака значение рН питательной воды больше или равно 8,0 (рН > 8,0). В этом случае слагаемыми Хн+[Н+ ] и [Н+ ] в уравнениях (5) и (6)
можно пренебречь ввиду малости их значений. Учитывая соотношение концентраций №+ и ЫН4+ не менее чем 1:5, можно принять
Хыа [N8+ ] + Хын+ [ЫН+ ] = ([N8+ ] + [ЫН+ ]) - 70 ,
тогда решением уравнений (5) и (6) относительно концентрации гидроксильных ионов ОН-, получим
[ОН-] = 3,73-х - (0,546 + 0,427 п + 2,219пхн) [СГ].
(7)
Из (7) значение рН рассчитывается по выражению
рН = 6 - 1д
0,01 [ОН- ]
(8)
(10)
^N8+ , Х
'N1-1+
и др. - подвижности этих ионов,
равные значениям эквивалентной электропро-
Концентрация аммонийных ионов [ИН4+] выражается из уравнения (6) с учетом принятых допущений в виде
] = [С!-] + [НСО- ] + [ОН- ] - [N8+ ]. (9)
Для питательной воды принятых к рассмотрению энергетических котлов с хн < 0,3 мкСм/см можно записать [С!-] + [НСО3] + [ОН-] - [N8+] »
И [С!-] + [НСО-] + [ОН-], т.е. пренебречь в сумме слагаемым [№+]. Тогда концентрация ^Н+ ] (мкмоль/дмз) получится из выражения
[1\1Н4+]=3,73-х+[СГ] (0,454+0,573п+2,98пхН).
(11)
Концентрация аммиака [ИНз], мкг/дм , определяемая суммой равновесных концентраций [ЫН4+]+[ЫН4ОН], будет равна
[ЫНз]=[ЫН4+](17+0,904[ОН-]). (12)
Результаты исследования. Результаты расчета по представленной методике значений рН и концентрации аммиака [ЫНз] при удельной электропроводности пробы 2,5 мкСм/см приведены на рис. 1, 2. Анализ полученной расчетной зависимости (рис. 1) показывает, что при увеличении удельной электропроводности Н-катионированной пробы от 0,1 до 1,0 мкСм/см снижение рН находится в пределах 0,2 единиц при всех значениях «п», тогда как увеличение хн до 1,5 мкСм/см влечет резкое уменьшение рН с 8,8 до 8,15 уже при увеличении «п» от 0,1 до 1,0 и далее - до 7,35 при значении п = 1,5. Изменение концентрации аммиака при этом имеет более равномерный характер, увеличиваясь от 240 мкг/дмз при ХН = 0,1 мкСм/см (и любых «п») до 364 мкг/дмз при хн = 1,5 мкСм/см и п = 1,5.
Рис. 1. Расчетная зависимость рН предельно разбавленного водного раствора с удельной электропроводностью 2,5 мкСм/см
Рис. 2. Расчетная зависимость концентрации аммиака предельно разбавленного водного раствора с удельной электропроводностью 2,5 мкСм/см
Подобные расчеты, выполненные при значении удельной электропроводности пробы Х = 5 мкСм/см, дают аналогичные с рис. 1, 2 результаты, увеличивая концентрацию аммиака и уменьшая диапазон изменений рН.
Обсуждение и выводы. Проведенные расчеты показывают, что для сверхчистых сред питательной воды энергоблоков ТЭС провал значений рН в коррозионно опасную область (рН < 8,0) может происходить даже в условиях существенных дозировок аммиака (более 250 мкг/дм3) при увеличении концентрации углекислоты, например, вследствие ухудшения работы деаэратора питательной воды.
Использование предложенной методики для расчета рН и концентраций ионных примесей питательной воды энергоблоков различных ТЭС подтверждает возможность ее применения для анализа качества ВХР и сделанный выше вывод (см. таблицу).
Таким образом, принятая методика позволяет по измеренным и приведенным к 25 оС значениям удельной электропроводности (х и ХН) охлажденных проб питательной воды энергетических котлов производить косвенные (расчетные) измерения концентраций хлоридов, гидрокарбонатов, аммиака и значения рН, что обеспечивает возможность автоматического химического контроля и управления водным
режимом современных энергоблоков в широком диапазоне изменения параметров и концентраций ионов, в том числе, связанной углекислоты.
Список литературы
1. Бушуев Е.Н. Математическое моделирование ионных равновесий водного теплоносителя // Теплоэнергетика. - 2009. - № 7. - С. 13-18.
2. Ларин Б.М., Ларин А.Б., Колегов А.В.
Измерения электропроводности и рН в системах мониторинга водного режима ТЭС. - Иваново, 2014. -С. 57-161.
3. Ларин Б.М., Ларин А.Б. Современное состояние химического контроля качества водного теплоносителя на ТЭС // Теплоэнергетика. - 2016. - № 5. -С. 70-74.
4. Мостофин А.А. Уточнение показаний кондуктометров с предвключенными Н-катионитовыми фильтрами // Электрические станции. - 1974. - № 1. -С. 79-81.
5. Кострикин Ю.М., Коровин В.А., Рубчинская С.М.
Влияние повышения температуры пробы на значение рН и удельную электрическую проводимость // Теплоэнергетика. - 1982. - № 1. - С. 76.
6. Bellows J.C., Weaver K.L. An on-line Steam Cycle Chemistry diagnostic System // Philadelphia. USA. ASME IEEE Power Generation Conference. - 1988. - С. 34-40.
7. Emory H. Hill, Robert D. Bartholomew. Rigorous calculation of sodium-to-phosphate mole rations for phosphate treatment programs / Power plant chemistry. - 2006. - 8 (9). -Р. 526-536.
8. Воронов В.Н., Петрова Т.И. Проблемы организации водно-химических режимов на тепловых электростанциях // Теплоэнергетика. - 2002. - № 7. - С. 2-6.
Результаты измерений и расчета нормируемых и диагностических показателей качества питательной воды
Название Измеренные данные Расчетные данные
объекта n x, мкСм/см хн, мкСм/см рН [NH3], мкг/дм3 рН [СП, з мкг/дм [HCO3-], мкг/дм [NH3], мкг/дм3
ГТЭС Терешково 1,0 3,77 0,221 9,13 9,14 0,27 0,58 426,13
1,5 9,13 0,22 0,70 427,17
Петрозаводская 0,55 4,6 0,225 9,22 587 9,23 0,35 0,42 564,90
ТЭЦ 1,0 9,22 0,27 0,59 566,52
Северо-Западная 0,55 8,04 0,19 9,47 1400 9,47 0,30 0,32 1328,72
ТЭЦ 1,0 9,47 0,23 0,46 1330,24
0,55 0,196 0,187 7,62 7,64 0,29 0,30 18,16
Конаковская ГРЭС 0,1 7,62 0,23 0,43 18,98
0,1 0,453 0,212 8,16 8,15 0,46 0,09 35,64
0,55 8,13 0,33 0,38 37,49
References
1. Bushuev, E.N. Matematicheskoe modelirovanie ionnykh ravnovesiy vodnogo teplonositelya [Mathematical modeling of ionic balances of the water heat carrier]. Teploenergetika, 2009, no. 7, pp. 13-18.
2. Larin, B.M., Larin, A.B., Kolegov, A.V. Izmereniya elektroprovodnosti i pH v sistemakh monitoringa vodnogo rezhima TES [Measurements of conductivity and pH in water regime monitoring systems of thermal power plants]. Ivanovo, 2014, pp. 57-161.
3. Larin, B.M., Larin, A.B. Sovremennoe sostoyanie khimicheskogo kontrolya kachestva vodnogo teplonositelya na TES [The current state of chemical quality control of the water heat carrier at thermal power plants]. Teploenergetika, 2016, no. 5, pp. 70-74.
4. Mostofin, A.A. Utochnenie pokazaniy konduktome-trov s predvklyuchennymi H-kationitovymi fil'trami [Calibration
of conductivity meters with H-cation prefilters]. Elektricheskie stantsii, 1974, no. 1, pp. 79-81.
5. Kostrikin, Yu.M., Korovin, V.A., Rubchinskaya, S.M. Vliyanie povysheniya temperatury proby na znachenie pH i udel'nuyu elektricheskuyu provodimost' [Influence of test temperature increase on pH and specific electric conductivity]. Teploenergetika, 1982, no. 1, p. 76.
6. Bellows, J.C., Weaver, K.L. An on-line Steam Cycle Chemistry diagnostic System. Philadelphia. USA. ASME IEEE Power Generation Conference, 1988, pp. 34-40.
7. Emory, H. Hill, Robert, D. Bartholomew. Rigorous calculation of sodium-to-phosphate mole rations for phosphate treatment programs. Power plant chemistry, 2006, 8 (9), pp. 526-536.
8. Voronov, V.N., Petrova, T.I. Problemy organizatsii vodno-khimicheskikh rezhimov na teplovykh elektrostantsiyakh [Problems of water and chemical regimes organization at thermal power plants]. Teploenergetika, 2002, no. 7, pp. 2-6.
Ларин Андрей Борисович,
ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,
кандидат технических наук, доцент кафедры химии и химических технологий в энергетике,
e-mail: [email protected]
Larin Andrei Borisovich,
Ivanovo State Power Engineering University,
Candidate of Engineering Sciences (PhD), Associate Professor of the Department of Chemistry and Chemical Technologies in Power Engineering, e-mail: [email protected]
Сорокина Анастасия Ярославовна,
ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,
инженер кафедры химии и химических технологий в энергетике,
e-mail: [email protected]
Sorokina Anastasia Yaroslavovna,
Ivanovo State Power Engineering University,
Engineer of the Department of Chemistry and Chemical Technologies in Power Engineering, e-mail: [email protected]