Методика расчета рН и концентраций ионных примесей питательной воды на ТЭС по измерениям удельной электропроводности Текст научной статьи по специальности «Химические технологии»

УДК 621.187.11

Методика расчета рН и концентраций ионных примесей питательной воды на ТЭС по измерениям удельной электропроводности

А.Б. Ларин, А.Я. Сорокина ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,

г. Иваново, Российская федерация E-mail: [email protected]

Авторское резюме

Состояние вопроса: Ужесточение норм, предъявляемых к качеству водного теплоносителя ТЭС, и развитие автоматических систем химического мониторинга за водно-химическим режимом требуют расширения диагностических возможностей систем мониторинга ТЭС. Использование наиболее точных измерений показателей качества теплоносителя, таких как удельная электрическая проводимость и рН, повышает надежность работы котельных установок и позволяет проводить диагностику состояния водного режима. Такое направление принято и в энергетике западных стран. В частности, фирмой Swan (Швейцария) разработаны и представлены на энергетическом рынке России автоматические анализаторы FAM Deltacon pH и AMI Deltacon Power. В связи с этим необходима разработка методик косвенных измерений и приборов автоматического химического контроля, не уступающих по техническому уровню лучшим мировым образцам.

Материалы и методы: Методики получены на основании глубоких исследований в области теории растворов электролитов, теории и практики ионных равновесий в водных растворах с использованием инновационных стендов кафедры ХХТЭ. Подобные методики и стенды впервые созданы в ИГЭУ и не имеют аналогов в других вузах России.

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

Ключевые слова: расчетная методика, питательная вода паровых котлов, показатель рН, математическая модель, удельная электропроводность, автоматический химический контроль, водно-химический режим ТЭС, нормы качества теплоносителя.

Methods of calculating pH and concentrations of feed water ionic impurities at heat power stations by measuring specific water conductivity

A.B. Larin, A.Ya. Sorokina Ivanovo State Power Engineering University, Ivanovo, Russian Federation E-mail: [email protected]

Abstract

Background: Tightening standards for the quality of thermal power plants water coolants and development of automatic chemical monitoring systems for water-chemical regime make it necessary to extend the diagnostic capacity of thermal power plant monitoring systems. Use of the most accurate measurement of coolant quality parameters, such as specific electrical conductivity and pH, increases boiler systems reliability and enables water regime diagnostics. This direction is adopted in power engineering in Western countries. In particular, the Swan firm (Switzerland) has developed and introduced to the Russian energy market its automatic analyzers «FAM Deltacon pH» and «AMI Deltacon Power». The purpose of our research is to develop techniques of indirect measurements and devices of automatic chemical control which are technologically not inferior to the best world samples.

Materials and methods: The developed techniques are the result of in-depth studies conducted in the field of the electrolyte solutions theory, the theory and practice of ionic balance in water solutions by using the innovative stands of the Department of Chemistry and Chemical Technologies in Power Engineering. Similar techniques and stands were for the first time designed at ISPEU and have no analogs in other higher education institutions of Russia. Results: The study deals with three variants of calculating indirect determination of pH and ammonia concentration based on specific conductivity measurements of refrigerated samples. The calculation techniques are based on specific solutions to the generalized mathematical model of ionic equilibria in conditions of good-quality feed water in boiler installations. The paper represents the results of laboratory research into and industrial testing of the new techniques of indirect determination of water quality standard values at heat power stations. An algorithm has been selected to calculate parameters of industrial systems of automatic chemical control.

Conclusions: The suggested techniques and algorithms have been tested in industrial environments at heat power plants and partly implemented in developing new devices and automatic chemical control systems. The scientific value of the results obtained consists in receiving and justification of the calculation techniques suitable for practical use at heat power plants.

Key words: calculation technique, feed water in steam boilers, pH indicator, mathematical model, specific conductivity, automatic chemical control, heat power plant water-chemical mode, coolant quality standards.

РО1: 10.17588/2072-2672.2016.5.010-015

Состояние вопроса. Надежность, экономичность и безопасность работы энергоблоков ТЭС в значительной степени зависят от состояния водно-химического режима (ВХР),

регламентируемого правилами технической 1 2 эксплуатации и нормативными документами .

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

Для комплексного решения проблемы контроля состояния и ведения водно-химического режима (ВХР) энергетических котлов на ТЭС в последнее десятилетие широко используются системы химико-технологического мониторинга (СХТМ), основу которых составляют приборы АХК. Наиболее надежными приборами АХК в промышленной эксплуатации являются кондуктометры и рН-метры, диагностический потенциал которых используется не в полном объеме [1-3].

Материалы и методы. Методы математического моделирования на основе измерений величины рН, удельной электропроводности прямой и Н-катионированной пробы позволяют расширить спектр функциональных возможностей СХТМ, раскрыть приборный потенциал и, следовательно, достаточно эффективно выявлять недостатки и быстротекущие нарушения ВХР [4, 5].

1 СО 153-34.20.501-2003 (РД 34.20.501-95). Правила технической эксплуатации электрических станций и сетей Российской Федерации / Минэнерго России. -М.: СПО ОРГРЭС, 2003.

2 СТО 70238424.27.100.013-2009. Водоподготови-тельные установки и водно-химический режим ТЭС. Условие создания. Нормы и требования. Стандарт

организации. - М.: НП «ИнВЭЛ», 2009. - 93 с.

Таким образом, контроль состояния, диагностика нарушений и ведение ВХР в рамках СХТМ остаются одними из наиболее сложных задач, решение которых следует искать в направлении повышения информативности приборов и систем АХК на базе простых и надежных измерений, а также в разработке новых расчетных методик, адаптированных к условиям промышленной эксплуатации.

В настоящее время ведется практическая реализация алгоритмов и программ, созданных на основе частных решений обобщенной математической модели (ММ). Разрабатываемые алгоритмы и методики используются для создания новых систем мониторинга водного режима энергетических котлов ТЭС, в том числе блоков ПГУ. Появление на рынке высокотехнологичных приборов дает возможность измерения и расчета нескольких показателей качества при уменьшении количества устанавливаемых приборов на точку контроля [6]. Это направление отвечает мировой тенденции развития химического контроля на ТЭС и АЭС3.

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

Результаты исследования. Рассматриваются три варианта расчетной методики косвенного определения показателя рН и концентрации аммиака, основанные на измерениях удельной электропроводности охлажденных проб питательной воды, прямой (х) и Н-катионированной (хн). Способ основан на решении общей ММ ионных равновесий вод типа конденсата в граничных условиях качества питательной воды [7].

Обработку измеренных значений электропроводности проводят следующим образом. Приводят к температуре 25 оС, как требуют нормы контроля рабочей среды энергетических котлов [6], согласно следующим уравнениям: Х = Хизм /(1 + 0,02« - 25)); (1)

Хм =Хм,изм /(1 + 0,015(? - 25)). (2)

3 European standart EN 12952-12: 2003. Water - tube boilers and auxiliary installation - part 12: Requirement for boiler feed water and boiler water quality.

Проверяют достоверность результатов измерений на соответствие максимальному диапазону изменения величины удельной электропроводности, мкСм/см: 0,056 < х < 10; 0,056 < хН < 10.

Находят разность значений (х - хН) и определяют реакцию среды контролируемой пробы по следующему критерию:

- если х - хн > 0, то рН > 7 в присутствии подщелачивающего агента или без него;

- если х - хН < 0, то рН < 7 в присутствии подкисляющего агента.

Таким образом, определяют однозначную зависимость между удельной электропроводностью и величиной рН малобуферного предельно разбавленного водного раствора (рис. 1).

0,3

0,2

0,1

6,0

7,0

8,0

9,0

Рис. 1. Кривые, отвечающие значениям удельной электропроводности и рН при изменении концентрации водных растворов NaOH (С = 1-10-7 - 1-10-6 моль(экв)/дм3) и HCl (С = 1 -10-7 - 5-10-7 моль(экв)/дм3) для различных температур

Для чистой воды при температуре 25 оС удельная электропроводность связана с концентрацией ионов водорода, определяющих значение рН, следующим выражением: хн2О = 0,548 • Сн + , (3)

где хНгО измерена в мкСм/см, а концентрация

ионов водорода (С^ ) в мкмоль/дм .

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

1. Первый вариант расчетной методики предполагает анализ ионных равновесий питательной воды при отсутствии солевых примесей, когда в воде есть только угольная кислота

и аммиак. Рассматриваются равновесия как в фильтрате Н-катионитовой колонки, так и в пробе воды до Н-колонки.

Предположим, что в фильтрате Н-колон-ки может быть только угольная кислота. Тогда, согласно уравнению диссоциации по первой ступени

Н2СО3 ^ Н+ + НСО3- , имеем:

• уравнение, описывающее диссоциацию угольной кислоты в Н-фильтрате: С . • С

К, = 4,5 • 10 7 =■ Н НС°3'

(4)

12CO3

• уравнение электропроводности для Н-фильтрата:

Хн-103 = (Xlrh + ) - С ,+ = 394,3 - С . ;

лн v и+ hco" 'иг ' и+ '

'3

(5)

• уравнение электронеитральности для Н-фильтрата: С . = С

и+ hco:

3

(6)

(7)

(8) (9)

Из (4)-(6) следует: (^н+ )н = 2,536 хн • 10 6, моль/дм3;

(С )н = 2,536 хн • 10-6, моль/дм3;

(Сн2со3)н = 14,3 хН • 10-6, моль/дм3.

Обозначим (Снсо-)н + (Си2со3)н = а.

Запишем систему уравнений ионных равновесий для прямой пробы воды (до Н-ко-лонки) с учетом равенства суммы форм углекислоты до- и после Н-колонки:

• уравнение электронейтральности для анализируемой воды:

С . = С nh+

14 oh исо3 cog

• балансовое уравнение форм состояния

+ С..

hco"

+ 2С

co2

(1С)

углекислоты:

С — + Сго2_ = (C^CO)H + (^.о,)и = А;

(11)

• уравнения, описывающие диссоциацию слабых электролитов в анализируемой воде: а) гидрокарбоната:

K,, = 4,8 -10"11 =

11 СИ - ^O2"

С,

HCO"

б) воды: Kw = 1-1014 = С..+ - С

и+

OH"

(12)

(13)

• уравнение электропроводности для анализируемой воды:

х -10"3 = x + С + + x С + x С +

к nh+ nh+ oh" oh" hco" hco3

+2 -I 2 С 2

co2 co2

или

х-10"3 = 73,6 - cnи4 + 198,3 - +

+44,5 - ^3

при x, = x0,, ,

2 - 69,3 - С 2

co2

(14)

где х0,( - предельная подвижность ионов при бесконечном разбавлении.

В уравнения (4)-(14) входят следующие величины:

• входные данные - показатели приборов АХК:

х, хН - измеряемые удельные электропроводности в анализируемой воде и Н-фильтрате соответственно, мкСм/см, приведенные к 25 оС;

• выходные данные - искомые концентрации ионов в анализируемой воде:

Сн,^н+ Снсо3Ссо2-Сон - концентрации соответствующих ионов в анализируемой воде, моль/дм ; (с^ )и,(^о- )н, (Сон- )н - концентрации соответствующих ионов в фильтрате Н-фильтра, моль/дм ;

• параметры (внутренние характеристики математической модели):

Кш, КI, К//, КМН4ОН - концентрационные константы ионных равновесий воды, углекислоты по первой и второй ступеням и аммиака; , х

Н ЫНд

и т.д. - предельные подвижности (или эквивалентные электропроводности) соответствующих ионов, См-см /г-экв.

Тогда расчетная система состоит из пяти уравнений (10)-(14), содержит пять неизвестных величин (не считая х и хн) и имеет одно решение. Решать можно методом подстановки или другим численным методом. Недостатком является большое различие численных значений определяемых величин, например СН+ = 110-9 при СМН4+ = 110-5, что может создавать значительные ошибки в определении рН или выводить решение за пределы допустимых значений, например получение отрицательных значений концентраций.

2. Второй вариант повторяет исходные данные первого варианта (угольная кислота и аммиак), но использует итерационный метод решения системы уравнений (10)-(14) с начальным заданием приближенного значения рН. При изменении рН в диапазоне от 8,5 до 9,5, что обычно имеет место для питательной воды энергетических котлов, задаемся средним значением рН1 = 9,0 лишь для преобразования уравнений (11) и (12).

Тогда из уравнения (12) получим

Тогда

С

К - Снсо-

со2

С

н+

4,8 -10 1-10-

11

С = 0,048 - С .

нсо- нсо-

Подставляя выражение Ссо2- в уравне

ние (11), получим: Снсо- (1+0,048)=А;

С

А

нсо:

- 1,048

А

С 2 = 0,048—— = 0,046 - А,

со2 '

1,048

а уравнение (10) будет иметь вид А

С^ = С + —— + 2 - 0,046 - А = 1,046 - А + С . N1-14 он 1,048 он-

Далее, подставляя СNн+, Снсо-, СсоГ в

уравнение (14), получим: х -10-3 = 73,6(1,046 - А + Сш-) +198,3 - Сон_ +

А

+44,5—А— + 2 - 69,3 - 0,046 - А, 1,048

или

х-10-3 = 77 - А + 73,6 - С__+198,3 - С__+ 48,84 - А,

' он- ' он-

или

С =Х- 103-125,84 - А = Х- 3,68-10-3-0,46 - А, он- 271,9 Х

3

моль/дм3, где х

Н2со3

(15)

измеряется в мкСм/см;

3

А = (Сн 2со)н + (Снссг )н , моль/дм .

Зная А и х, находим С и рН в виде

он-

рН2 = -1д

1-10

14

С

он-

Если расчетное значение рН отличается от принятого более чем на 0,1, что определяется нормативной точностью измерения рН, то повторяем расчет по методу 2 с заданием уточненного значения рН равного (рН1 + рН2)/2.

Далее находим концентрацию ЫН4+ по уравнению (10) и концентрацию ЫН3: СN^ = Сж+ -10-6 (1+5,32-10-4-Сон-)-17, мкг/дм3. (16)

Такой вариант привлекает простотой решения задачи с возможностью реализации в промышленных системах АХК.

3. Третий вариант (или оценочный метод) предполагает упрощение задачи путем исключения из исходных данных концентрации угольной кислоты. В абсолютно чистой воде растворен только ЫН4ОН. Тогда достаточно одного измерения удельной электропроводности (х) без Н-колонки.

Рассчитываем ионное равновесие: Ж4ОН ^ 1ЧН4+ + ОН-.

Уравнение диссоциации записывается следующим образом:

С....+ - С„

Кл = 1,76 -10-5 = _^н д р

(17)

NЩOH

Уравнение электропроводности имеет следующий вид:

Х -10-3 = (\,и+ + ^) Сми+ = 271,9 - С и+ (18)

или

См н+ = Сан- = 3,68 - Х -106, моль/дм3. (19)

Уравнение диссоциации воды можно представить в виде

9

Кт = 1-1014 = С^- Сон- . (20)

Из уравнений (17)^(20) следует:

рН = - 1д ^272 -10-9 ^ , (21)

где х - измеряется в мкСм/см,

и, подставляя выражение (19) в выражение

(16), получим

СNн = 62,56 - х +13,07-х2, мкг/дм3. (22)

В данном варианте расчетной методики будет получаться максимальное значение рН, так как при наличии углекислоты часть ионов ОН- будет расходоваться на связывание Н+ от диссоциации Н2СО3, и рН будет уменьшаться. Это может служить проверкой измеряемых значений рН в условиях сверхчистой питательной воды прямоточных котлов СКД и ССКП.

В таблице приведены результаты измеренных и расчетных значений по трем вариантам расчетной методики.

Выводы. При больших дозировках ЫН3 (больше 100 мкг/дм3) результаты расчета рН достоверны в пределах паспортной ошибки измерения прибора рН-метра, концентрации ЫН3 отличаются на 5-10 %.

Для предельно разбавленных водных растворов, типа питательной воды прямоточных котлов с хм < 0,2 мкСм/см и х < 1 мкСм/см, приближенные методы расчета, включая и вариант 1, могут давать отклонения значений рН от измеренных, превышающие паспортную ошибку прибора ± 0,05. Концентрация ЫН3 в этих случаях обычно не измеряется, так как в этих условиях химический анализ на ЫН3 не дает достоверных результатов. Тогда практически единственным методом остается расчетное определение ЫН3.

Результаты измеренных и расчетных значений рН и концентрации аммиака в питательной воде парового котла на энергетическом предприятии

Исходные данные Результаты измерений Результаты расчета

Энергетическое Вариант 1 Вариант 2 Вариант 3

предприятие x, мкСм/см xa мкСм/см рн Cnh3, мкг/дм3 рн смиэ , мкг/дм3 рн смиэ , мкг/дм3 рн смиэ , мкг/дм3

Петрозаводская 4,6 0,225 9,22 587 9,21 559,16 9,19 563,77 9,23 564,34

ТЭЦ 4,13 0,238 9,17 499 9,16 480,88 9,14 444,78 9,18 481,31

СевероЗападная ТЭЦ 8,04 0,19 9,47 1400 9,46 1303,22 9,46 1305,87 9,47 1347,85

Конаковская 0,196 0,187 7,62 - 7,44 21,56 7,47 21,96 7,86 12,76

ГРЭС 0,453 0,212 8,16 - 7,76 30,69 8,06 42,12 8,22 32,02

Костромская ГРЭС 0,686 0,108 8,0 - 8,37 52,95 8,37 53,00 8,4 49,07

Для питательной воды котлов СКД и ССКП возникает другая проблема, связанная с точностью измерений рН в сверхчистых средах. Сравнительный анализ измеренных значений х и рН (табл. 1) показывает, что измеренные значения х питательной воды блока Конаковской ГРЭС (х = 0,453 мкСм/см) меньше, чем аналогичная величина блока Костромской ГРЭС (х = 0,686 мкСм/см), а измеренное значение рН (рН = 8,16) больше аналогичного значения Костромской ГРЭС (рН = 8,0), чего не может быть. Поэтому расчетное значение рН = 8,37 более достоверное, по сравнению с измеренной величиной. Для условий питательной воды энергоблоков СКД и ССКП требуется разработка специальной методики калибровки рН-метров.

Представленные варианты расчетной методики базируются на штатных измерениях в системах химического контроля и не требуют сколько-нибудь существенных изменений и дополнений в системы АХК. Основу предложенных расчетных методик составляет сочетание частных решений обобщенной ММ ионных равновесий и измерений удельной электропроводности прямой и Н-катионированной

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

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

1. Воронов В.Н., Петрова Т.И. Проблемы организации водно-химических режимов на тепловых электростанциях // Теплоэнергетика. - 2002. - № 7. -С. 2-7.

2. Богачев А.Ф., Радин Ю.А., Герасименко О.Б. Особенности эксплуатации и повреждаемость котлов-утилизаторов бинарных парогазовых установок. - М.: Энергоатомиздат. - 2008. - С. 232.

3. Егошина О.В., Воронов В.П., Назаренко М.П. Современное состояние систем химико-технологического мониторинга на тепловых станциях на основе опыта МЭИ - ООО НИЦ «Элемент» // Теплоэнергетика. - 2013. - № 7.

4. Расчетный метод определения потенциально кислых веществ в питательной воде прямоточных котлов / Б.М. Ларин, Е.Н. Бушуев, А.Б. Ларин и др. // Теплоэнергетика. - 2008. - № 4. - С. 38-41.

5. Определение концентрации фосфатов в котловой воде путем измерения электропроводно-

сти / Б.М. Ларин, Е.Н. Бушуев, Ю.Ю. Тихомирова, С.В. Киет // Теплоэнергетика. - 2008. - № 7. - С. 21-27.

6. Ларин Б.М., Ларин А.Б., Колегов А.В. Измерения электропроводности и рН в системах мониторинга водного режима ТЭС. - Иваново, 2014. -332 с.

7. Бушуев Е.Н. Математическое моделирование ионных равновесий водного теплоносителя с использованием измерения электропроводности и рН // Теплоэнергетика. - 2009. - №7. - С. 13-18.

References

1. Voronov, V.N., Petrova, T.I. Problemy organi-zatsii vodno-khimicheskikh rezhimov na teplovykh elek-trostantsiyakh [Problems of organizing water-chemical modes at thermal power plants]. Teploenergetika, 2002, no. 7, pp. 2-7.

2. Bogachev, A.F., Radin, Yu.A., Gerasimenko, O.B. Osobennosti ekspluatatsii i povrezhdaemost' kotlov-utilizatorov binarnykh parogazovykh ustanovok [Operation and damageability of heat recovery boilers in combined-cycle power plants]. Moscow, Energoatomizdat, 2008, p. 232.

3. Egoshina, O.V., Voronov, V.P., Nazarenko, M.P. Sovremennoe sostoyanie sistem khimiko-tekhnologicheskogo monitoringa na teplovykh stant-siyakh na osnove opyta MEI - OOO NITs «Element» [State-of-the-art of chemical and technological monitor-

ing systems at heat power stations based on the experience of MEI - LLC Research Center Element]. Teploenergetika, 2013, no. 7.

4. Larin, B.M., Bushuev, E.N., Larin, A.B., Batti, M.K.L., Vinogradov, V.N., Kiet, S.V. Raschetnyy metod opredeleniya potentsial'no kislykh veshchestv v pitatel'noy vode pryamotochnykh kotlov [A calculation technique of identifying potentially sour substances in feedwater of non-circulation boilers]. Teploenergetika, 2008, no. 4, pp. 38-41.

5. Larin, B.M., Bushuev, E.N., Tikhomiro-va, Yu.Yu., Kiet, S.V. Opredelenie kontsentratsii fosfa-tov v kotlovoy vode putem izmereniya elektroprovod-nosti [Determining of phosphates concentration in boiler water through conductivity measurement]. Teploenergetika, 2008, no. 7, pp. 21-27.

6. Larin, B.M., Larin, A.B., Kolegov, A.V. Izmereniya elektroprovodnosti i рН v sistemakh monitoringa vodnogo rezhima TES [Measurements of conductivity and pH in heat power plant water mode monitoring systems]. Ivanovo, 2014. 332 p.

7. Bushuev, E.N. Matematicheskoe modelirova-nie ionnykh ravnovesiy vodnogo teplonositelya s is-pol'zovaniem izmereniya elektroprovodnosti i pH [Mathematical modeling of ionic balances of the water heat carrier by measuring conductivity and pH]. Teploenergetika, 2009, no. 7, pp. 13-18.

Ларин Андрей Борисович,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры химии и химических технологий в энергетике, e-mail: [email protected]

Сорокина Анастасия Ярославовна,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», инженер кафедры химии и химических технологий в энергетике, e-mail: [email protected]