Теоретические основы химической технологии, 2021, T. 55, № 6, стр. 690-706
Сравнение сорбционной способности различных марок активированного угля по отношению к аргону, криптону и ксенону природного изотопного состава в статических условиях
Э. П. Магомедбеков a, *, А. О. Меркушкин a, А. В. Обручиков a, Д. А. Сахаров a
a Российский химико-технологический университет им. Д.И. Менделеева
Москва, Россия
* E-mail: eldar@muctr.ru
Поступила в редакцию 16.06.2021
После доработки 25.06.2021
Принята к публикации 12.07.2021
Аннотация
Определены константы Генри при адсорбции инертных газов природного изотопного состава (аргон, криптон, ксенон) на активированном угле марок АГ-3, ВСК-5, СКТ-3, NWC 12x40 в температурном интервале 15–80°С. Установлена положительная корреляция констант Генри с удельной поверхностью, объемом и поверхностью микропор, отнесенных к единице объема твердой фазы. Определены параметры температурной зависимости констант Генри и рассчитаны изостерические теплоты адсорбции инертных газов. Полученные экспериментальные результаты использованы для расчета стационарного режима работы радиохроматографической колонны системы спецгазоочистки атомных электростанций. Показано, что из всех исследованных марок наибольшие значения коэффициентов очистки воздушного потока от инертных радиоактивных газов обеспечивают кокосовые угли марок ВСК-5 и NWC 12x40, заметно превосходящие по этому показателю применяемый в настоящее время на большинстве российских АЭС торфяной уголь серно-калиевой активации марки СКТ-3. Уголь марки АГ-3 не позволяет достичь требуемой степени очистки.
ВВЕДЕНИЕ
Для всех технологических процессов и объектов ядерной энергетики характерной особенностью является присутствие источников радиационного риска, обусловленного выбросами и сбросами радиоактивности, которые при определенных условиях могут приводить к негативным воздействиям на человека и окружающую среду. Выбросы и сбросы могут быть как регламентными (постоянными или периодическими), находящимися под контролем эксплуатационного персонала, так и аварийными. В условиях нормальной эксплуатации атомных электростанций (АЭС) выбросы радионуклидов во внешнюю среду незначительны и состоят в основном из радионуклидов йода и инертных радиоактивных газов (ИРГ) [1, 2]. Проведенный авторами [3] анализ данных Европейской Комиссии по газовым выбросам АЭС Европы [4] показал, что доля ИРГ в суммарной активности выбросов составляет от 55 до 90%.
Согласно рекомендации Международного агентства по атомной энергии (МАГАТЭ) [5] при эксплуатации АЭС важным фактором воздействия на окружающую среду следует считать 12 радиоактивных изотопов инертных газов: 41Ar, 85Kr, 85mKr, 87Kr, 88Kr, 131mXe, 133Xe, 133mXe, 135Xe, 135mXe, 137Xe, 138Xe. Изотопы криптона и ксенона образуются при делении ядерного топлива; радиоактивный изотоп 41Ar образуется при активации нейтронами стабильного изотопа 40Ar, содержащегося в воздухе.
По данным работы [6], в России при нормальной эксплуатации АЭС 99% дозы облучения населения от выбросов ИРГ формируются следующими семью изотопами: 41Ar, 85mKr, 87Kr, 88Kr, 133Xe, 133mXe, 135Xe. Проведенная в работе [1] оценка радиационного воздействия на население от выбросов ИРГ на АЭС Европы выявила следующие изотопы, вносящие наиболее значительный вклад (указан в скобках) в формирование эффективной дозы: 41Ar (51%), 138Хе (17%), 88Kr (13.6%), 135Хе (7.2%), 135mХе (4.1%), 87Kr (3.4%), 133Хе (3.1%). Остальные изотопы ИРГ обеспечивают около 0.5% радиационного воздействия. Таким образом, очевидна необходимость снижения активности инертных радиоактивных газов в воздушных потоках АЭС перед их сбросом в атмосферу.
Для удаления ИРГ из воздушных потоков могут быть использованы различные технологии – криогенная дистилляция [7], экстракция фреонами [8], мембранная сепарация [9], адсорбция на высокопористых материалах [10]. Криогенная дистилляция – хорошо изученный процесс, широко применяемый для выделения стабильных изотопов криптона и ксенона из воздуха. Однако она требует специальной предварительной подготовки газа-носителя (в частности, удаления паров воды и оксидов азота), проводится при давлении, значительно превышающем атмосферное, что повышает риски, связанные с недостаточной герметичностью технологических контуров и возможностью протечек при извлечении радиоактивных изотопов инертных газов с высокой удельной активностью. Кроме того, существует опасность появления на этапе дистилляции взрывоопасных концентраций озона.
При экстракции смеси ИРГ жидким фреоном накопление коррозионноактивных продуктов радиолиза фреона накладывает дополнительные требования, предъявляемые к конструкционным материалам, используемым для изготовления оборудования, что неизбежно приведет к удорожанию процесса. Использование для выделения из воздуха и разделения радиоактивных криптона и ксенона при помощи каскада селективных органических мембран [9] осложняется вследствие радиационного повреждения последних, а также их высокой стоимости.
Одними из наиболее эффективных для удаления летучих радиоактивных примесей из воздуха являются адсорбционные процессы. Закономерности адсорбционных процессов достаточно хорошо изучены, а оборудование является типовым, надежным и недорогим в изготовлении и эксплуатации. Эффективными адсорбентами для извлечения из воздушных потоков инертных газов являются новые высокопористые материалы – углеродные нанотрубки, цеолиты (в частности, модифицированный морденит в H- и Ag-формах, демонстрирующий высокую сорбционную способность по отношению к криптону и ксенону соответственно [11, 12]), нанокомпозиционные сорбенты [13], а также новые перспективные пористые материалы с металлорганическим каркасом (metal-organic frameworks – MOFs) [10, 14, 15]. Следует отметить, что многие сорбенты на основе цеолитов довольно сложны в синтезе, кроме того, использование серебра для их модификации значительно повышает стоимость получаемых материалов. Основным препятствием для широкого промышленного использования сорбентов типа MOFs также является сложность и высокая стоимость их крупнотоннажного синтеза [10]. Старейшими, наиболее изученными и эффективными, а также сравнительно недорогими адсорбентами для инертных газов и летучих форм радиоиода являются активированные угли различных степеней и способов активации, которые могут применяться как в гранулированном виде, так и в составе композиционных материалов, обладающих меньшим гидравлическим сопротивлением при сопоставимой величине коэффициента очистки [16–18].
В настоящее время на АЭС функцию очистки радиоактивных технологических газовых сдувок оборудования при различных режимах эксплуатации энергоблоков от ИРГ, аэрозолей и радиоактивного йода-131 выполняет система спецгазоочистки (СГО), основное оборудование которой включает самоочищающиеся фильтры, цеолитовые фильтры и фильтры-адсорберы. Фильтры-адсорберы являются главным элементом системы СГО, где происходит высокоэффективная очистка газа от изотопов радиоактивных инертных газов (ксенон, криптон) и радиоактивных изотопов йода методом сорбции на активированном угле и выдержки [19].
Радиоактивные изотопы инертных газов при прохождении газового потока через фильтрующий слой адсорбируются частицами активированного угля и на некоторое время задерживаются на них. Таким образом, происходит медленное движение фронта радионуклидов по высоте колонны фильтра, где за счет естественного радиоактивного распада активность адсорбированных газов убывает вдоль слоя угля и должна достигать допустимого уровня к выходу из колонны.
Количество активированного угля, загружаемое в адсорберы системы СГО, должно соответствовать толщине слоя, обеспечивающей прохождение адсорбционных фронтов ИРГ через адсорбционные колонны за время, достаточное для их распада (за исключением долгоживущего 85Kr) до допустимого уровня суммарной объемной активности. Это время (время защитного действия) будет зависеть от природы адсорбтива, температуры, скорости газового потока и свойств адсорбента. Изменение этих параметров, а также замена адсорбента в адсорберах (например, вследствие перехода на активированный уголь другой марки или другого производителя) потребует подтверждения достаточности загруженного количества для обеспечения необходимой степени очистки газового потока от ИРГ.
Согласно действующему в России отраслевому стандарту [19] определение сорбционной способности используемого в системе СГО активированного угля по отношению к криптону и ксенону проводится в статических условиях. Для этого весовым методом получают 4 точки изотермы адсорбции при абсолютном давлении адсорбтива 0.5, 1, 5 и 10 кПа и температуре 25°С. Для заключения о пригодности испытуемого материала результаты сравнивают с результатами, полученными на стандартном образце угля. В настоящее время фильтры-адсорберы систем СГО на большинстве российских АЭС снаряжены торфяным углем сульфатно-калиевой активации типа СКТ-3 [20]. Детально адсорбция инертных газов на угле этой марки применительно к системе СГО исследована в [21, 22]. Однако, данная марка угля более не выпускается промышленностью, при этом актуальные данные по адсорбции инертных газов на коммерчески доступных современных марках активированных углей в условиях, близких к условиям работы системы СГО фактически отсутствуют. Многие современные публикации (например, [23–26]), ориентированы больше на процессы разделения и выделения инертных газов и описывают адсорбцию аргона, криптона и ксенона на различных высокопористых материалах при давлениях, значительно превышающих атмосферное. Данные, полученные для давлений до 1 атм при температурах вблизи комнатной, зачастую относятся к 1950–1970-м годам [27, 22 ].
Ключевой характеристикой, определяющей способность активированного угля удерживать адсорбат, является константа Генри (или коэффициент адсорбции). Именно от нее зависит скорость движения по радиохроматографической колонне фронта адсорбированного ИРГ.
Целью настоящей работы было получение изотерм адсорбции инертных газов (аргона, криптона и ксенона) природного изотопного состава на активированном угле различных марок в интервале температур 15–80°C, вычисление констант Генри и расчет эффективности удаления из газового потока основных дозообразующих изотопов аргона, криптона и ксенона в условиях, близких к реальным условиям работы системы СГО (в предположении о пренебрежимо малом влиянии бета-излучения сорбционное равновесие).
ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ
Изотермы адсорбции снимали на приборе Quadrasorb Kr/SI (Quantachrome Instruments). В качестве адсорбтивов были использованы инертные газы особой чистоты (99.9995%) – аргон, криптон и ксенон производства AKELA-N LTD. В качестве адсорбентов использовали активированные угли коммерчески доступных марок – ВСК-5, СКТ-3, АГ-3 (АО “ЭНПО Неорганика”) и NWC 12x40 (NWC Carbon). Перед каждым экспериментом адсорбенты подвергались дегазации в FloVac Degasser (Quantachrome Instruments) – нагревании в вакууме при 350°С до прекращения газовыделения (давление в нагретой вакуумированной ячейке изменялось менее чем на 0.001 мм рт. ст. в течение одной минуты). Перед отсоединением от системы дегазации и подключением к Quadrasorb ячейка с образцом заполнялась гелием под давлением ∼800 мм рт. ст. Навески дегазированных образцов составляли 8–12 г при адсорбции аргона и криптона и 1.0–2.0 г при адсорбции ксенона. Измерение свободного объема ячеек с адсорбентом проводилось автоматически перед каждым экспериментом с использованием гелия марки 6.0 (99.9999%). Термостатирование ячеек с адсорбентом на протяжении всего эксперимента осуществлялось с точностью ±0.1°С при помощи сосуда с водяной рубашкой, подключенного к жидкостному термостату LOIP LT 100.
Удельную поверхность образцов рассчитывали с использованием метода Брунауэра–Эммета–Теллера (БЭТ) по изотермам низкотемпературной (77.35 K) адсорбции азота. Площадь, занимаемую каждой адсорбированной молекулой азота, принимали равной 0.162 нм2 [28]. Для расчета характеристик пористой структуры по изотермам адсорбции азота использовали метод теории функционала плотности (QSDFT) [29, 30] (ядро модели для комбинации щелевидных и цилиндрических пор в углеродном адсорбенте). Математическую обработку изотерм адсорбции азота проводили с использованием специализированного программного обеспечения QuadraWin 7.1 (Quantachrome Instruments).
Измерение насыпной плотности с утряской проводили при помощи анализатора насыпной плотности Autotap (Quantochrome Instruments). Кажущуюся плотность частиц (отношение их массы к занимаемому ими объему) определяли методом гидростатического взвешивания в этаноле гранул угля, предварительно пропитанных расплавленным парафином, предотвращающим попадание спирта в поры.
Долю свободного объема в слое сорбента рассчитывали, исходя из значений насыпной и кажущейся плотности по уравнению
(1)
${{\varepsilon }} = \left( {1 - \frac{{{{{{\rho }}}_{{\text{b}}}}}}{{{{{{\rho }}}_{{\text{a}}}}}}} \right),$Удельный объем открытых пор и открытую пористость рассчитывали по уравнениям (2)–(3), измеряя давление гелия до и после напуска из калиброванной емкости (Quadrasorb) в ячейку известного объема с навеской образца:
где vо – удельный объем открытых пор, см3/г; Πo – открытая пористость образца, %; m – навеска образца, г; Vc – объем пустой ячейки для образца, см3; Vk – объем калиброванной емкости с гелием, см3; Pk – давление гелия в калиброванной емкости до напуска в ячейку с образцом, Па; Pс – давление гелия ячейке с образцом до напуска гелия из калиброванной емкости, Па; PHe – давление гелия в системе после напуска в ячейку с образцом, Па.РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
Характеристика исходных материалов. Перед адсорбционными экспериментами у исходных материалов были измерены насыпная плотность, кажущаяся плотность частиц, удельный объем открытых пор, открытая пористость и доля свободного объема в насыпном слое. Характеристики исходных материалов сведены в табл. 1
Таблица 1.
Характеристика | АГ-3 | ВСК-5 | СКТ-3 | NWC 12x40 |
---|---|---|---|---|
ρb, г/см3 | 0.499 | 0.438 | 0.414 | 0.539 |
ρa, г/см3 | 0.799 | 0.807 | 0.636 | 0.922 |
vo, см3/г | 0.802 | 0.774 | 1.133 | 0.590 |
Πo, % | 64.1 | 62.5 | 72.1 | 54.4 |
ε, % | 37 | 45 | 35 | 42 |
Характеристика пористой структуры исходных материалов. Для определения параметров пористой структуры исходных материалов были получены изотермы низкотемпературной (77.35 K) адсорбции и десорбции азота на образцах СКТ-3, ВСК-5, АГ-3 и NWC 12x40 (рис. 1). Все изотермы можно отнести к типу IVa по классификации Международного союза чистой и прикладной химии (ИЮПАК) [28], характерному для микро-мезопористых адсорбентов.
Изотерма адсорбции на образце NWC 12x40 ближе к типу 1b, характерному для сорбентов с преобладанием микропор и относительно малой внешней поверхностью. В области относительных давлений 0.45–1.0 на изотермах присутствует петля адсорбционного гистерезиса, причем ширина петли, измеренная как разность ординат десорбции и адсорбции, приблизительно постоянна во всей области гистерезиса. Такая характерная особенность, а также вид изотермы позволяют отнести ее к типу H4 по классификации ИЮПАК [28]. Петля адсорбционного гистерезиса наиболее выражена для образца СКТ-3, наименее – для NWC 12x40 и объясняется различными механизмами адсорбции и десорбции газа в мезопорах. Петля гистерезиса типа Н4 присуща системам с набором неоднородных капилляров, имеющих широкую часть полости и разнообразные по величине короткие узкие горла, и, по-видимому, отражает комбинацию микропор и мезопор, встречающуюся в некоторых типах активированных углей. Гистерезис на изотерме адсорбции для углей, имеющий форму Н4, означает, что уже при сравнительно небольшом давлении произошло заполнение микропор и дальнейшая адсорбция практически прекратилась. Наблюдающийся крутой подъем адсорбционной ветви изотермы вблизи Р/Р0 = 1 для образцов СКТ-3 и ВСК-5 может указывать, на то, что углеродная масса между полостями, имеющими размер мезо- и макропор, пронизана микропорами близкой к щелевидной форме [31, 32].
По изотермам адсорбции была рассчитана удельная поверхность адсорбентов с использованием метода БЭТ. Поскольку для микропористых адсорбентов выбор адекватного диапазона относительных давлений применимости метода БЭТ может быть затруднен, использовали рекомендации, приведенные в [28], позволяющие избежать субъективности при определении емкости монослоя: 1) значение энергетической константы C в уравнении БЭТ должно быть положительным; 2) применимость уравнения БЭТ должна быть ограничена диапазоном давлений, в котором член A(1 – P/P0) монотонно возрастает с увеличением P/P0; 3) Значение P/P0, соответствующее найденной емкости монослоя AM, должно лежать внутри выбранного диапазона. Все изотермы, построенные в координатах БЭТ в выбранных диапазонах относительных давлений, были линейны (r2 > 0.9999). Рассчитанные значения удельной поверхности представлены в табл. 2.
Таблица 2.
Характеристика | АГ-3 | ВСК-5 | СКТ-3 | NWC 12x40 |
---|---|---|---|---|
SBETm, м2/г | 750 | 1189 | 1471 | 1083 |
SBETv, м2/см3 | 599 | 960 | 936 | 999 |
vmm, см3/г | 0.205 | 0.399 | 0.472 | 0.362 |
vmv, см3/см3 | 0.164 | 0.322 | 0.300 | 0.334 |
vzm, см3/г | 0.178 | 0.091 | 0.232 | 0.077 |
vzv, см3/см3 | 0.142 | 0.073 | 0.148 | 0.071 |
vm, см3/г | 0.383 | 0.490 | 0.704 | 0.439 |
vv, см3/см3 | 0.306 | 0.395 | 0.448 | 0.405 |
Smm, м2/г | 475 | 1080 | 1173 | 917 |
Smv, м2/см3 | 380 | 872 | 746 | 845 |
Szm, м2/г | 189 | 92 | 177 | 99 |
Szv, м2/см3 | 151 | 74 | 113 | 91 |
Sm, м2/г | 664 | 1172 | 1350 | 1016 |
Sv, м2/см3 | 531 | 946 | 859 | 937 |
Расчет распределения пор по размерам с использованием модели QSDFT проводили по обеим ветвям изотермы (адсорбция + десорбция) для комбинации щелевидных и цилиндрических пор во всем исследованном интервале относительных давлений. Во всех случаях данный метод позволил достаточно хорошо описать экспериментальные изотермы – расхождение расчетных значений и экспериментальных данных были незначительны и составляли десятые и сотые доли процента. При этом попытка использования ядра QSDFT для пор другой формы (только щелевидных, цилиндрических, сферических, а также комбинации цилиндрических и сферических) приводило к существенно большим (в несколько раз) относительным отклонениям, особенно в области малых давлений (P/P0 < 5 × 10–3).
Интегральные кривые объема пор в зависимости от их ширины приведены на рис. 2. Каждая вертикальная “ступенька” на графиках рис. 2 соответствует порам определенной ширины, причем крутизна “ступеньки” говорит об узости соответствующего пика на дифференциальной кривой распределения пор по размерам – чем круче “ступенька”, тем более узким является распределение пор по размерам в этой области.
Суммарный объем микропор (d < 2 нм) и мезопор (2 нм < d < 50 нм), а также их удельная поверхность приведены в табл. 2. Анализ полученных данных выявил во всех образцах три группы пор. Наиболее значительный вклад в удельную поверхность (до 90%) и объем (до 80%) вносят микропоры с шириной менее 1.5 нм и модой ширины не более 0.7 нм (резкий подъем интегральных графиков удельного объема на рис. 2 до первой “ступеньки”). Наибольшим содержанием микропор характеризуются образцы ВСК-5 и NWC 12x40, наименьшим – АГ-3 (табл. 2). Резкий обрыв дифференциальных кривых распределения (рис. 3) со стороны левой границы диапазона (d ≈ 0.72 нм) при использовании ядра QSDFT может быть связан, с одной стороны, с недостатком экспериментальных точек для области P/P0 < 10–5, соответствующих заполнению узких микропор, а с другой – по-видимому, с особенностью реализации алгоритма решения интегрального уравнения адсорбции, не накладывающей требования равенства нулю объема пор на границах диапазона [30].
Две другие группы пор представлены узкими мезопорами с модами ширины 2.2–2.4 нм и 5.0–5.3 нм (рис. 3). Вторая группа мезопор наиболее выражена в образце СКТ-3, первая – в АГ-3. Для этого образца четкой границы между мезопорами первой и второй группы нет, можно говорить о широком распределении размеров мезопор в диапазоне 2.0–5.8 нм. Судя по наклону графиков на рис. 2 при d > 5 нм, в случае СКТ-3, АГ-3 и, в меньшей степени, ВСК-5, можно говорить также о третьей группе мезопор, не имеющей явно выраженной моды распределения. Иными словами, данные образцы содержат некоторое количество относительно широких мезопор, распределенных во всем диапазоне 5–50 нм. Относительный вклад мезо- и микропор в общую удельную поверхность и пористость образца приведен в табл. 2. Там же приведены значения объема микро- и мезопор, а также их поверхности, отнесенные к единице массы и единице объема твердой фазы.
Константы Генри и коэффициенты адсорбции аргона, криптона и ксенона. Для расчета констант Генри были проведены эксперименты по адсорбции аргона, криптона и ксенона на каждом адсорбенте в диапазоне температур 15–80°C при давлении адсорбтива от 0.76 до 700 мм рт. ст.
Известно, что изотермы адсорбции криптона и ксенона на углеродных адсорбентах имеют форму, близкую к изотерме Ленгмюра [22]. Однако, поскольку парциальные давления инертных газов (за исключением аргона) в воздушной смеси, направляемой в систему СГО, незначительны (вследствие невысокого содержания в воздухе стабильных изотопов и малого периода полураспада образующихся в реакторе радиоактивных изотопов), можно считать, что их адсорбция в области малых давлений подчиняется закону Генри. Для более точного определения констант Генри экспериментальные данные были обработаны с использованием уравнения Ленгмюра (4), достаточно точно описывающего начальный участок изотерм адсорбции и переходящего в уравнение Генри при малых давлениях адсорбтива:
Параметры уравнения (4) находили численно методом наименьших квадратов (МНК). Исходные данные изотерм с помощью уравнения Менделеева–Клапейрона предварительно преобразовывали так, чтобы концентрация адсорбтива в объемной и твердой фазах была приведена к одним и тем же единицам (мг/см3). Константы Генри, рассчитанные по таким преобразованным исходным данным, являются безразмерными и имеют физический смысл концентрационных констант равновесия адсорбции в области малых давлений.
Верхний порог давления (Pmax) при отборе точек для обработки по МНК выбирали таким, чтобы относительное среднеквадратичное отклонение результатов расчета от экспериментальных данных не превышало 1%. Этот порог давления уменьшался от аргона к ксенону и составлял 600, 50 и 10 мм рт. ст. для аргона, криптона и ксенона соответственно. Это связано с уменьшением давления насыщенного пара адсорбтива, вследствие чего отклонение изотермы от линейности для ксенона начинаются в области более низких давлений, чем для криптона и аргона, что затрудняет непосредственное определение константы Генри для ксенона по углу наклона начального участка изотермы. Обработка данных с использованием уравнений (4), (5) позволила обойти это затруднение. Схожий подход был использован в [33] при описании начальных участков изотерм адсорбции аргона, криптона и ксенона на цеолитах типа A. Рассчитанные значения констант Генри при разных температурах приведены в табл. 3.
Таблица 3.
T, °С | АГ-3 | ВСК-5 | СКТ-3 | NWC 12x40 | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Ar | Kr | Xe | Ar | Kr | Xe | Ar | Kr | Xe | Ar | Kr | Xe | |
15 | 4.47 | 29.5 | 469 | 9.84 | 70.3 | 1136 | 8.13 | – | 876 | 10.1 | 77.2 | 1366 |
20 | 4.08 | 25.2 | 396 | 8.97 | 61.8 | 1002 | 7.42 | 43.5 | 740 | 9.29 | 67.3 | 1125 |
25 | – | – | – | 8.17 | 53.0 | – | – | – | – | – | – | – |
30 | 3.44 | 19.5 | 293 | 7.47 | 46.8 | 689 | 6.14 | 33.3 | 548 | 7.76 | 52.3 | 796 |
35 | – | – | – | 6.85 | 41.1 | – | – | – | – | – | – | – |
40 | 2.93 | 15.0 | 205 | 6.29 | 36.4 | 504 | 5.28 | 26.9 | 383 | 6.58 | 40.0 | 600 |
50 | – | – | – | 5.39 | – | – | – | – | – | – | – | – |
60 | 2.18 | 9.84 | 119 | 4.71 | 22.9 | 266 | 3.98 | 17.5 | 234 | 4.89 | 25.7 | 326 |
80 | 1.69 | 6.71 | 65.4 | 3.55 | 16.0 | 151 | 3.08 | 12.0 | 119 | 3.72 | 17.3 | 179 |
Поскольку количество адсорбента лимитируется вместимостью адсорбера, на практике для расчета необходимого количества угля обычно используют приведенный к насыпному объему коэффициент адсорбции KA (м3/м3), связанный с константой Генри соотношением
(6)
${{K}_{{\text{A}}}} = {{K}_{{\text{H}}}}\frac{{{{{{\rho }}}_{{\text{b}}}}}}{{{{{{\rho }}}_{{\text{a}}}}}} = {{K}_{{\text{H}}}}\left( {1 - {{\varepsilon }}} \right).$По данным, представленным в табл. 1 и 3, были рассчитаны коэффициенты адсорбции инертных газов на различных марках активированного угля, учитывающие разницу в насыпной плотности и кажущейся плотности частиц (табл. 4). Они позволяют сопоставить эффективность задерживания ИРГ на различных марках активированного угля в одинаковых условиях.
Таблица 4.
T, °С | АГ-3 | ВСК-5 | СКТ-3 | NWC 12x40 | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Ar | Kr | Xe | Ar | Kr | Xe | Ar | Kr | Xe | Ar | Kr | Xe | |
15 | 2.79 | 18.4 | 293 | 5.34 | 38.1 | 617 | 5.3 | – | 570 | 5.91 | 45.2 | 799 |
20 | 2.55 | 15.8 | 247 | 4.87 | 33.5 | 544 | 4.8 | 28.3 | 482 | 5.43 | 39.4 | 658 |
25 | – | – | – | 4.43 | 28.8 | – | – | – | – | – | – | – |
30 | 2.15 | 12.2 | 183 | 4.05 | 25.4 | 374 | 4.0 | 21.7 | 357 | 4.54 | 30.6 | 466 |
35 | – | – | – | 3.72 | 22.3 | – | – | – | – | – | – | – |
40 | 1.83 | 9.36 | 128 | 3.41 | 19.8 | 274 | 3.4 | 17.5 | 249 | 3.85 | 23.4 | 351 |
50 | – | – | – | 2.93 | – | – | – | – | – | – | – | – |
60 | 1.36 | 6.15 | 74.3 | 2.56 | 12.4 | 144 | 2.6 | 11.4 | 152 | 2.86 | 15.0 | 191 |
80 | 1.06 | 4.19 | 40.8 | 1.93 | 8.7 | 82.0 | 2.0 | 7.83 | 77.5 | 2.18 | 10.1 | 105 |
Сравнение данных табл. 4 показало, что из четырех исследованных марок активированный уголь АГ-3 будет наименее эффективен при использовании в адсорбере системы СГО – по коэффициентам адсорбции всех трех инертных газов он значительно уступает остальным маркам угля. Коэффициенты адсорбции аргона для марок СКТ-3, ВСК-5 и NWC 12x40 близки между собой. Различие увеличивается при переходе к криптону и ксенону – для этих улавливания этих газов наиболее эффективным оказался уголь марки NWC 12x40, наименее – уголь СКТ-3. Однако обладающий наименьшим размером частиц NWC 12x40 будет создавать значительное гидравлическое сопротивление при работе адсорбера. В этой связи оптимальным выбором, возможно, будет уголь марки ВСК-5.
Связь констант Генри с параметрами пористой структуры адсорбентов. С практической точки зрения полезным было бы установить связь между константами Генри при адсорбции инертных газов на активированном угле с некоторыми общими параметрами пористой структуры и характеристиками сорбентов, измеряемых независимо. Это важно, поскольку при выборе сорбента не всегда есть возможность экспериментально исследовать его сорбционную способность к заданному адсорбтиву. Попытки установить такую связь предпринимались и ранее. Так, авторами [22] было установлено, что ряд по насыпной плотности углей марок СКТ соответствует ряду по коэффициенту адсорбции криптона и ксенона при малых парциальных давлениях. Теми же авторами было показано, что отсутствует корреляция между константами Генри и удельной поверхностью, измеренной по методу БЭТ или объема микропор, вычисленному в рамках теории Дубинина–Радушкевича. Исследование в настоящей работе активированных улей разных типов, существенно различающихся по своим свойствам, не подтвердило наличия эмпирической связи насыпной плотности с константами Генри. Так, обладающий наименьшей насыпной плотностью уголь СКТ-3 значительно превосходит по адсорбционной способности уголь АГ-3. Последний же, превосходя по насыпной плотности уголь ВСК-5, существенно уступает ему в значениях констант Генри для всех трех инертных газов.
Сопоставление табл. 2 и 3 показывает, что отнесенные к массе сорбента удельная поверхность по БЭТ (SBETm) и объем микропор (vmm), очевидно, также не могут служить критерием оценки сорбционной способности активированного угля по отношению к инертным газам. Так, обладающий максимальными удельной поверхностью и объемом микропор (1471 м2/г и 0.472 см3/г) уголь СКТ-3 имеет более низкие значения констант Генри для криптона и ксенона по сравнению с углями ВСК-5 и NWC 12x40 с удельной поверхностью 1189, 1083 м2/г и объемом микропор 0.399 и 0.362 см3/г соответственно. NWC 12x40, немного уступая ВСК-5 в удельной поверхности и объеме микропор, превосходит его в сорбционной способности криптона и ксенона.
Причина, по нашему мнению, может заключаться в том, что рассчитанные значения констант Генри относятся к объему твердой фазы (как отношение концентраций адсорбтива в подвижной и неподвижной фазах). Вследствие этого выстраивание корреляционных связей с удельными величинами, отнесенными к массе сорбента, не учитывающими кажущуюся плотность частиц, не будет корректным. В табл. 2 наряду с отнесенными к массе сорбента поверхностью и объемом пор приведены также значения этих показателей, пересчитанные на объем частиц.
Для выявления корреляционной связи нами были рассчитаны коэффициенты линейной корреляции Пирсона [34] констант Генри при 20°C с характеристиками пористой структуры адсорбентов, отнесенных к объему твердой фазы (табл. 5). Аналогичные расчеты, проведенные с константами Генри для других температур, дали близкие результаты.
Таблица 5.
Характеристика | Ar | Kr | Xe |
---|---|---|---|
SBETv, м2/см3 | 0.972 | 0.906 | 0.925 |
vmv, см3/см3 | 0.984 | 0.926 | 0.942 |
vzv , см3/см3 | –0.789 | –0.894 | –0.866 |
vv, см3/см3 | 0.744 | 0.593 | 0.633 |
Smv, м2/см3 | 0.989 | 0.938 | 0.946 |
Szv, м2/см3 | –0.960 | –0.941 | –0.932 |
Sv, м2/см3 | 0.987 | 0.930 | 0.941 |
Данные, представленные в табл. 5, показали наличие заметной положительной корреляционной зависимости между константами Генри и такими отнесeнными к объeму твeрдой фазы показателями, как удельная поверхность по БЭТ (SBETv), объем и поверхность микропор (vmv, Smv) для всех трех адсорбтивов. Отрицательная корреляция констант Генри с поверхностью и объемом мезопор (Szv, vzv) вызвана, очевидно, отрицательной корреляцией объема микропор с объeмом мезопор: при увеличении доли мезопор в образцах, объем и поверхность микропор уменьшается. Сильная положительная корреляция констант Генри с суммарным объемом (поверхностью) микро- и мезопор (vv, Sv) объясняется решающим вкладом в эти показатели именно микропор (см. табл. 2).
Известно, что параметры пористой структуры во многом зависят от способа активации. Было показано, что химическая активация гидроксидом и карбонатом калия дает более высокое содержание микропор, чем получается при использовании таких активирующих агентов, как ZnCl2 или H3PO4, применение которых увеличивает содержание мезопор в материале [35–37]. Таким образом, знание связи характеристик пористой структуры с константами Генри потенциально дает возможность целенаправленно проводить синтез адсорбентов с заданными свойствами. Следует отметить, что надежное выявление статистической взаимосвязи таких общих характеристик адсорбента, как удельная поверхность или объем микропор, требует обработки достаточной выборки экспериментальных данных. Этому вопросу будет посвящено отдельное исследование.
Параметры температурной зависимости констант Генри. Истинная константа равновесия связана с термодинамическими функциями системы и концентрационной константой Генри через соотношения
где ΔG – стандартное изменение энергии Гиббса при адсорбции; ΔH – стандартное изменение энтальпии при адсорбции; ΔS – стандартное изменение энтропии при адсорбции; K – истинная константа равновесия адсорбции; T – абсолютная температура; R – универсальная газовая постоянная, 8.31446 Дж/(моль K); f – множитель, включающий в себя коэффициенты активности адсорбтива в обеих фазах.Из уравнений (7) и (8) следует, что в некотором температурном интервале, где величины ΔH и ΔS можно считать не зависящими от температуры, а множитель f постоянным, логарифм константы Генри будет линейной функцией обратной температуры:
где a и b – коэффициенты уравнения линейной регрессии.Обработка данных из табл. 3 в соответствии с уравнением (9) показала, что в изученном температурном интервале (15–80°C) логарифм констант Генри линейно зависит от обратной температуры (рис. 4), что указывает также на линейность изостер адсорбции, построенных в координатах ln(P)–1/T, в этом температурном интервале. Линейность изостер, в свою очередь, позволяет при помощи уравнения Клапейрона–Клаузиуса, представленного в форме (10), рассчитать дифференциальную изостерическую теплоту адсорбции при допущении о применимости законов идеального газа [31, 38, 39]:
(10)
${{\left( {\frac{{\partial {\text{ln}}\left( P \right)}}{{\partial \left( {{1 \mathord{\left/ {\vphantom {1 T}} \right. \kern-0em} T}} \right)}}} \right)}_{{A = {\text{const}}}}} = - \frac{{{{q}_{{{\text{st}}}}}}}{R},$В табл. 6 приведены рассчитанные значения коэффициентов уравнения (9) при адсорбции аргона, криптона и ксенона на исследованных марках активированного угля, а также квадраты коэффициента линейной корреляции (r2) и значения дифференциальной изостерической теплоты адсорбции при достаточно малой степени заполнения поверхности адсорбентов, соответствующей условиям экспериментов.
Таблица 6.
Марка угля | Параметр | Ar | Kr | Xe |
---|---|---|---|---|
АГ-3 | a | –3.796 | –4.645 | –4.459 |
b | 1.525 × 103 | 2.309 × 103 | 3.064 × 103 | |
qst, кДж/моль | 12.68 | 19.20 | 25.48 | |
r2 | 0.9999 | 0.9995 | 0.9984 | |
ВСК-5 | a | –3.234 | –3.873 | –4.017 |
b | 1.590 × 103 | 2.340 × 103 | 3.197 × 103 | |
qst, кДж/моль | 13.22 | 19.46 | 26.58 | |
r2 | 0.9998 | 0.9994 | 0.9984 | |
СКТ-3 | a | –3.159 | –3.769 | –3.821 |
b | 1.512 × 103 | 2.209 × 103 | 3.061 × 103 | |
qst, кДж/моль | 12.57 | 18.37 | 25.45 | |
r2 | 0.9996 | 0.9998 | 0.9950 | |
NWC 12x40 | a | –3.120 | –3.841 | –3.676 |
b | 1.567 × 103 | 2.361 × 103 | 3.142 × 103 | |
qst, кДж/моль | 13.03 | 19.63 | 26.12 | |
r2 | 0.9999 | 0.9999 | 0.9986 |
Дифференциальные изостерические теплоты адсорбции мало различаются для исследованных адсорбентов и увеличивается в ряду адсорбатов Ar < Kr < Xe. Полученные значения в целом согласуются с данными других авторов [24] для активированных углей в том же интервале температур и давлений, и несколько превосходят (примерно на 2–4 кДж/моль) значения для цеолитов типа А [33] и MOFs на основе NiDOBDC [15]. Несколько более высокие значения (порядка 30 кДж/моль) приводятся в [22] для ксенона при его адсорбции из гелия на активированных углях марки СКТ.
Линейность изостер адсорбции инертных газах на микропористых сорбентах неоднократно отмечалась в других работах и, по-видимому, является общим свойством, характерным для такого рода систем [24]. Линейная аппроксимация изостер адсорбции (или зависимости логарифма константа Генри от обратной температуры) имеет важное практическое значение, позволяя рассчитывать адсорбционные равновесия в области параметров, выходящих за рамки экспериментальных исследований. Следует отметить, что параллельность прямых на рис. 4 дает возможность проводить сравнение сорбционной способности различных марок активированного угля на основании лишь одной-двух изотерм адсорбции.
Аналогично проведенным ранее расчетам для констант Генри нами были рассчитаны коэффициенты линейной корреляции Пирсона параметров уравнения (9) с удельной поверхностью по БЭТ, объемом и поверхностью микропор (табл. 7).
Таблица 7.
Характеристика | Параметр a | Параметр b | ||||
---|---|---|---|---|---|---|
Ar | Kr | Xe | Ar | Kr | Xe | |
SBETv, м2/см3 | 0.989 | 0.976 | 0.934 | 0.497 | 0.078 | 0.579 |
vmv, см3/см3 | 0.978 | 0.964 | 0.917 | 0.549 | 0.130 | 0.627 |
Smv, м2/см3 | 0.947 | 0.941 | 0.858 | 0.630 | 0.183 | 0.710 |
Анализ табл. 7 показал наличие заметной корреляционной связи энтропийного слагаемого уравнения (9) (параметр a) с такими общими характеристиками адсорбента, как удельная поверхность по БЭТ, объем и поверхность микропор, отнесенными к единице объема твердой фазы. Значимая их связь с параметром b, отражающим влияние энтальпийного фактора, напротив, выявлена не была.
Полученные экспериментальные данные по адсорбции инертных газов на различных марках активированного угля могут быть использованы для расчета стационарного режима работы радиохроматографической колонны, снаряженной активированным углем этих марок.
Расчет стационарного режима работы радиохроматографической колонны. Главной задачей радиохроматографической колонны в системе СГО является задержка поступающих в нее ИРГ на время, требуемое для их распада до приемлемых уровней активности. Поэтому ключевой характеристикой адсорбера, определяющей необходимое количество адсорбента, будут скорость движения фронта адсорбции по колонне и время защитного действия. Обе величины связаны с коэффициентом адсорбции через уравнение Н.А. Шилова [40]:
где τ – время защитного действия; Vb – насыпной объем адсорбента; w – объемный расход газового потока на входе в колонну; τ' – потеря времени защитного действия.Уравнение (11) с хорошим приближением может быть использовано для оценки минимально необходимого количества адсорбента при работе радиохроматографической колонны в стационарном изотермическом режиме, а также для расчета коэффициента очистки при известном объеме сорбента в адсорбере. Пользуясь законом радиоактивного распада и уравнением (11) при допущении τ $ \gg $ τ', можно связать с требуемым коэффициентом очистки от радионуклида в газовом потоке минимально необходимый объем адсорбента Vmin:
(12)
${{V}_{{{\text{min}}}}} = \frac{{w{{t}_{{1/2}}}}}{{{{K}_{{\text{A}}}}}}\frac{{\ln \left( {DF} \right)}}{{\ln \left( 2 \right)}},$Аналогичным образом можно выразить коэффициент очистки от конкретного радионуклида, задавая насыпной объем адсорбента Vb:
(13)
$DF = {\text{exp}}\left( {\frac{{{{V}_{{\text{b}}}}{{K}_{{\text{A}}}}}}{{w{{t}_{{1/2}}}}}\ln \left( 2 \right)} \right).$Поскольку в газовой смеси присутствуют изотопы разных элементов, необходимо знать вклад каждого в объемную активность исходной смеси. Суммарный коэффициент очистки может быть рассчитан в этом случае по уравнению
(14)
$\begin{gathered} D{{F}_{{{\Sigma }}}} = {{\left[ {\mathop \sum \limits_j \mathop \sum \limits_i {{{{\omega }}}_{{j,i}}}{\text{exp}}\left( { - \frac{{{\text{ln}}\left( 2 \right)}}{{{{t}_{{1/2}}}_{{j,i}}}}\frac{{{{V}_{{\text{b}}}}{{K}_{A}}_{j}}}{w}} \right)} \right]}^{{ - 1}}} = \\ = {{\left[ {\mathop \sum \limits_j \mathop \sum \limits_i \frac{{{{{{\omega }}}_{{j,i}}}}}{{D{{F}_{{j,i}}}}}} \right]}^{{ - 1}}}, \\ \end{gathered} $Величина Vmin не может быть выражена из уравнения (14) аналитически, для ее расчета необходимо применять численные методы решения трансцендентных уравнений.
Согласно [41] оптимальное время удержания изотопов ксенона и криптона, закладывающееся при проектировании радиохроматографической колонны для реакторов с водяным теплоносителем под давлением (PWR), составляет 40 и 2 сут соответственно. Рассчитанное с использованием уравнения Шилова минимально необходимое для этого количество активированного угля каждой из исследованных марок приведено в табл. 8 для разных температур эксплуатации адсорбционной колонны. Номинальный объемный расход очищаемого воздуха был принят равным 12 м3/ч [20, 41]. Для сравнения заметим, что на АЭС с реакторами этого типа российского производства система спецгазоочистки состоит из двух или трех взаимозаменяемых ниток, рабочий объем угольных адсорберов в каждой из которых составляет 20 или 40 м3 (экспортный вариант) [42].
Таблица 8.
T, °С | АГ-3 | ВСК-5 | СКТ-3 | NWC 12x40 |
---|---|---|---|---|
10 | 32 | 13 | 17 | 11 |
20 | 46 | 19 | 24 | 16 |
30 | 64 | 27 | 34 | 23 |
40 | 89 | 37 | 47 | 32 |
Минимально необходимые объемы угля АГ-3 значительно превосходят рабочий объем угольных адсорберов, что делает уголь этой марки непригодным для использования в системе СГО АЭС. Минимально необходимые для требуемого удержания изотопов криптона и ксенона объемы углей марок ВСК-5 и NWC 12x40 заметно ниже, чем марки СКТ-3 и близки между собой по этому показателю. Причем гидравлическое сопротивление колонны с углем ВСК-5 должно быть существенно ниже вследствие более крупных частиц угля этой марки.
Анализируя данные табл. 8, следует иметь в виду, что переход от Vmin к реальному объему адсорбента часто осуществляют, вводя в уравнение (12) коэффициент запаса (обычно 1.5 [21]), учитывающий неиспользованную статическую емкость (вклад τ'), а также флуктуации сорбционной способности сорбента от партии к партии.
Результаты расчета коэффициента очистки по уравнению (13) для отдельных нуклидов ИРГ при 20°С и насыпном объеме угля в адсорбере 20 м3 приведены в табл. 9, из которой видно, что короткоживущие изотопы ксенона (135Xe, 135mXe, 137Xe, 138Xe) и криптона (87Kr) полностью распадаются за время их движения по радиохроматографической колонне в штатном режиме эксплуатации (расчетный коэффициент очистки превышает 106) при использовании любой из исследованных марок активированного угля. Слабее всего удерживаются изотопы 41Ar (вследствие наименьших значений коэффициента адсорбции) и 131mXe (вследствие сравнительно большого периода полураспада). Активность изотопа 85Kr, очевидно, также не будет заметно снижаться из-за его большого периода полураспада 10.739 лет.
Таблица 9.
Нуклид | t1/2 [45] | АГ-3 | ВСК-5 | СКТ-3 | NWC 12x40 |
---|---|---|---|---|---|
20°C | 20°C | 20°C | 20°C | ||
41Ar | 109.6 мин | 5.1 | 23 | 21 | 30 |
85mKr | 4.48 ч | 61 | 5.7 × 103 | 1.4 × 103 | 2.5 × 104 |
87Kr | 76.3 мин | >106 | >106 | >106 | >106 |
88Kr | 2.825 ч | 6.9 × 102 | 8.9 × 105 | 9.8 × 104 | >106 |
131mXe | 11.84 сут | 2.8 | 9.0 | 7.3 | 15 |
133Xe | 5.248 сут | 10 | 1.4 × 102 | 88 | 4.4 × 102 |
133mXe | 2.198 сут | 2.5 × 102 | 1.4 × 105 | 4.4 × 104 | >106 |
135Xe | 9.14 ч | >106 | >106 | >106 | >106 |
135mXe | 15.29 мин | >106 | >106 | >106 | >106 |
137Xe | 3.818 мин | >106 | >106 | >106 | >106 |
138Xe | 14.14 мин | >106 | >106 | >106 | >106 |
Для расчета общего коэффициента очистки по уравнению (14) необходима информация об объемной активности каждого радионуклида ИРГ на входе в адсорбер либо о его вкладе в объемную активность. В доступной литературе такие сведения найти не удалось. В качестве оценочных значений для различных изотопов криптона и ксенона мы использовали расчетные данные об их удельной активности в теплоносителе первого контура реакторов типа PWR по данным [43] (табл. 10). В этой же работе проведено сопоставление полученных авторами расчетных значений с результатами непосредственных измерений.
Таблица 10.
Нуклид | 85Kr | 85mKr | 87Kr | 88Kr | 131mXe | 133Xe | 133mXe | 135Xe | 135mXe |
---|---|---|---|---|---|---|---|---|---|
нКи/г | 45.8 | 162 | 149 | 274 | 39.4 | 2600 | 70.2 | 820 | 125 |
Вклад 41Ar и 138Xe в объемную активность ИРГ мы рассчитывали по отношению их активности в криогенной азотной ловушке к суммарной активности изотопов 85mKr, 87Kr, 88Kr, 133Xe, 135Xe, полученной методом прямой гамма-спектрометрии, для газовых выбросов промышленного тяжеловодного реактора на ФГУП “ПО "Маяк” по данным [44]. Они составили 70% для 41Ar и 3.3% для 138Xe. Рассчитанные таким образом оценочные значения вклада отдельных радионуклидов ИРГ в их общую объемную активность (ω) приведены в табл. 11.
Таблица 11.
Нуклид | 41Ar | 85Kr | 85mKr | 87Kr | 88Kr | 131mXe | 133Xe | 133mXe | 135Xe | 135mXe | 138Xe |
---|---|---|---|---|---|---|---|---|---|---|---|
ω, % | 38.8 | 0.6 | 2.2 | 2.1 | 3.8 | 0.5 | 36.0 | 1.0 | 11.4 | 1.7 | 1.8 |
Следует отметить, что данные в табл. 11 являются усредненными значениями, поскольку сильно зависят от многих факторов, включая условия работы реакторной установки и продолжительность кампании. Замечено, что вклад 41Ar в активность ИРГ снижается, а изотопа 133Xe возрастает к окончанию кампании, что связано с увеличением числа микротрещин в оболочках тепловыделяющих элементов, через которые радиоактивные изотопы ксенона поступают в теплоноситель первого контура [44].
Рассчитанные по уравнению (14) значения общих коэффициентов очистки DFΣ для потока очищаемого воздуха 12 м3/ч при разных температурах и объемах угля в адсорбере приведены в табл. 12.
Таблица 12.
Сорбент | 10°C | 20°C | 30°C | 40°C | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
20 м3 | 30 м3 | 40 м3 | 20 м3 | 30 м3 | 40 м3 | 20 м3 | 30 м3 | 40 м3 | 20 м3 | 30 м3 | 40 м3 | |
АГ-3 | 14 | 33 | 68 | 8.3 | 19 | 39 | 5.6 | 11 | 22 | 4.1 | 7.3 | 13 |
ВСК-5 | 64 | 130 | 154 | 37 | 97 | 140 | 21 | 60 | 111 | 12 | 33 | 70 |
СКТ-3 | 58 | 124 | 152 | 34 | 90 | 136 | 19 | 54 | 104 | 11 | 30 | 64 |
NWC 12x40 | 78 | 140 | 156 | 49 | 113 | 148 | 29 | 79 | 128 | 17 | 48 | 95 |
Суммарный коэффициент очистки является интегральным показателем эффективности работы адсорбера в конкретных условиях и учитывает соотношение активности радионуклидов ИРГ во входящем газовом потоке, их периоды полураспада, а также коэффициенты адсорбции каждого инертного газа на данном адсорбенте. Для снижения активности суммы ИРГ в очищаемом воздушном потоке на два порядка величины, необходимо, чтобы температура эксплуатации адсорбера с активированным углем марок СКТ-3, ВСК-5 и NWC 12x40 не превышала 30°C при объеме угля в адсорбере не менее 30 м3 (для 30°С – не менее 40 м3). Активированный уголь марки АГ-3 не сможет обеспечить требуемою кратность снижения общей активности ИРГ в воздушном потоке.
ЗАКЛЮЧЕНИЕ
В результате проведенных исследований были определены константы Генри при адсорбции инертных газов (аргон, криптон, ксенон) природного изотопного состава на активированном угле марок АГ-3, ВСК-5, СКТ-3, NWC 12x40 в интервале температур 15–80°С. У всех образцов были также определены общие характеристики пористой структуры: удельная поверхность по БЭТ, объем и поверхность микро- и мезопор, распределение пор по размерам. Анализ полученных данных позволил выявить положительную корреляцию констант Генри с удельной поверхностью по БЭТ, объемом и поверхностью микропор, отнесенных к единице объема твердой фазы. Установленные закономерности потенциально дают возможность давать сравнительную оценку сорбционной способности различных марок активированного угля по отношению к инертным газам на основании общих характеристик пористой структуры, а также подбирать условия синтеза (активации) для получения углей с максимальной сорбционной способностью.
Определены параметры температурной зависимости констант Генри и рассчитаны изостерические теплоты адсорбции инертных газов. Линейность полученных зависимостей дает возможность прогнозировать значения констант Генри при температурах, выходящих за экспериментально исследованный диапазон. Более того, показано, что сравнение сорбционной способности различных марок активированного угля при разных температурах можно проводить на основании лишь одной-двух изотерм адсорбции.
На основании полученных экспериментальных данных был проведен расчет стационарного режима работы радиохроматографической колонны системы спецгазоочистки АЭС. Показано, что из всех исследованных марок, наибольшие значения коэффициентов очистки воздушного потока от суммы ИРГ обеспечивают кокосовые угли марок ВСК-5 и NWC 12x40, заметно превосходящие по этому показателю применяемый в настоящее время на большинстве российских АЭС торфяной уголь серно-калиевой активации марки СКТ-3. Уголь марки АГ-3 не позволяет достичь требуемой степени очистки. Более крупные частицы угля ВСК-5 по сравнению с NWC 12x40 и, как следствие, его меньшее гидравлическое сопротивление потоку воздуха делают его перспективным материалом для применения в адсорберах системы СГО АЭС.
Работа выполнена при поддержке Российского химико-технологического университета им. Д.И. Менделеева, проект 2020-008.
ОБОЗНАЧЕНИЯ
А | величина адсорбции при нормальных условиях, см3/г |
a | первый коэффициент уравнения линейной регрессии |
B | коэффициент, связанный с константой равновесия адсорбции |
b | второй коэффициент уравнения линейной регрессии |
DF | заданный коэффициент очистки |
d | ширина пор, нм |
f | множитель, включающий в себя коэффициенты активности адсорбтива в обеих фазах |
ΔG | стандартное изменение энергии Гиббса при адсорбции, Дж/моль |
ΔH | стандартное изменение энтальпии при адсорбции, Дж/моль |
K | константа адсорбции, см3/см3 |
m | масса образца, г |
P | давление, Па |
q | теплота адсорбции, Дж/моль |
r | коэффициент линейной корреляции Пирсона |
S | удельная поверхность, м2/г или м2/см3 |
ΔS | стандартное изменение энтропии при адсорбции, Дж/(моль K) |
T | температура, °С или K |
t1/2 | период полураспада радионуклида |
V | объем, см3 |
ν | удельный объем, см3/г или см3/см3 |
w | объемный расход газового потока на входе в колонну, м3/мин |
ε | доля свободного объема в слое адсорбента |
Π | пористость, % |
ρ | плотность адсорбента, г/см3 |
τ | время защитного действия, мин |
τ' | потеря времени защитного действия, мин |
ω | относительный вклад в общую исходную объемную активность |
ИНДЕКСЫ
0 | насыщение |
А | адсорбция |
a | кажущийся |
BETm | по БЭТ к массе |
BETv | по БЭТ к объему |
b | насыпной |
с | ячейка |
H | Генри |
i | изотоп |
j | радионуклид |
k | калиброванная емкость |
M | монослой |
m | микро- и мезопор к массе |
max | максимальное значение |
min | минимальное значение |
mm | микропор к массе |
mv | микропор к объему |
o | открытый |
st | изостерический |
v | микро- и мезопор к объему |
zm | мезопор к массе |
zv | мезопор к объему |
Список литературы
Пыркова А.А., Екидин А.А., Антонов К.Л. Поступление инертных радиоактивных газов в атмосферу при нормальной эксплуатации АЭС // Физика. Технологии. Инновации. Сборник материалов VI Международной молодежной научной конференции, посвященной 70-летию основания Физико-технологического института УрФУ. 2019. С. 279.
Merkushkin A.O., Obruchikov A.V., Magomedbekov E.P. Iodomethane Vapor Oxidation in Air Medium // Theor. Found. Chem. Eng. 2019. V. 53. № 4. P. 538. [Меркушкин А.О., Обручиков А.В., Магомедбеков Э.П. Окисление паров иодметана в воздушной среде // Теор. осн. хим. технол. 2019. Т. 53. № 4. С. 440.]
Пышкина М.Д. Определение основных дозообразующих нуклидов в выбросах АЭС PWR И ВВЭР // Биосферная совместимость: чел., рег., технол. 2017. № 2(18). С. 98.
RADD – European Commission RAdioactive Discharges Database (2021) Discharge reports. https://europa.eu/radd/nuclideDischargeOverview.dox?pageID=NuclideDischargeOverview. Cited 21 May 2021.
INPRO Methodology for Sustainability Assessment of Nuclear Energy Systems: Environmental Impact (IAEA Nuclear Energy Series № NG-T-3.15). Vienna: IAEA, 2016.
Екидин А.А., Васильев А.В., Васянович М.Е. Современные технологии управления воздействием на окружающую среду как инструмент соблюдения принципа ALARA // Биосферная совместимость: чел., рег., технол. 2017. № 2 (18). С. 67.
Zabaluev Y.V. Management of radionuclides from reprocessing plant gaseous effluents // IAEA Bull. 1979. V. 21. № 1. P. 23.
Monson P.R. Noble gas confinement for reactor fuel melting accidents // Proceedings of the 18th DOE Nuclear Airborne Waste Management and Air Cleaning Conference. Baltimore, Maryland, 1984. V. 1.
Stern S.A., Wang S.C. Permeation cascades for the separation of krypton and xenon from nuclear reactor atmospheres // AIChE J. 1980. V. 26. № 6. P. 891.
Nandanwar S.U., Coldsnow K., Utgikar V., Sabharwall P., Aston E.D. Capture of harmful radioactive contaminants from off-gas stream using porous solid sorbents for clean environment – A review // Chem. Eng. J. 2016. V. 306. P. 369.
Ianovski D., Munakata K., Kanjo S., Yokoyama Y., Koga A., Yamatsuki S., Tanaka K., Fukumatsu T., Nishikawa M., Igarashi Y. Adsorption of noble gases on H-mordenite // J. Nucl. Sci. Technol. 2002. V. 39. № 11. P. 1213.
Daniel C., Elbaraoui A., Aguado S., Springuel-Huet M.A., Nossov A., Fontaine J.P., Topin S., Taffary T., Deliere L., Schuurman Y., Farrusseng D. Xenon Capture on Silver-Loaded Zeolites: Characterization of Very Strong Adsorption Sites // J. Phys. Chem. C. 2013. V. 117. № 29. P. 15122.
Nandanwar S.U., Coldsnow K., Utgikar V., Sabharwall P., Aston E.D., Zhang Y. Synthesis and characterization of ETS-10: supported hollow carbon nano-polyhedrons nanosorbent for adsorption of krypton at near ambient temperatures // Adsorption. 2016. V. 22. P. 129.
Jian L., Thallapally P.K., Strachan D. Metal–Organic Frameworks for Removal of Xe and Kr from Nuclear Fuel Reprocessing Plants // Langmuir. 2012. V. 28. № 31. P. 11584.
Thallapally P.K., Grate J.W., Motkuri R.K. Facile xenon capture and release at room temperature using a metal–organic framework – A comparison with activated charcoal // Chem. Commun. 2012. V. 48. № 3. P. 347.
Магомедбеков Э.П., Меркушкин А.О., Обручиков А.В. Композиционный материал для сорбционной очистки воздуха от летучих форм радиоактивного иода. Пат. 2717818 РФ. 2020.
Obruchikov A.V., Magomedbekov E.P., Merkushkin A.O. Removal of radioactive methyliodide from the gas stream with a composite sorbent based on polyurethane foam // Nucl. Eng. Technol. 2020.V. 52. № 5. P. 1093.
Obruchikov A.V., Magomedbekov E.P., Merkushkin A.O. The composite sorption material for radioiodine trapping from air stream and the method for its preparation // J. Radioanal. Nucl. Chem. 2020. V. 324. P. 331.
Стандарт организации СТО 1.1.1.03.004.1633-2019. Определение допустимого количества фильтроциклов и оценка остаточного срока службы фильтрующего элемента фильтров адсорберов системы спецгазоочистки. М.: Концерн “Росэнергоатом”, 2019.
Технологические системы реакторного отделения ВВЭР-1000 с РУ В-320. Часть 2. Вспомогательные системы. М.: Министерство РФ по атомной энергии. Концерн “Росэнергоатом”. Балаковская атомная электростанция.
Очкин Д.В., Алёшин А.А., Казаков В.В., Голиков И.В. О применении активированного угля СКТ-3С для радиохроматографической очистки газообразных радиоактивных отходов. М.: ЦНИИатоминформ, 1984. С. 21.
Nakhutin I.E., Ochkin D.V., Tret’yak S.A., Dekalova A.N. The adsorption of krypton and xenon at low partial pressures on industrial samples of activated carbon // Sov. At. Energy. 1976. V. 40. № 3. P. 364. [Нахутин И.Е., Очкин Д.В., Третьяк С.А., Декалова А.Н. Исследование адсорбции криптона и ксенона при малых парциальных давлениях на промышленных образцах активного угля // Ат. энерг. 1976. Т. 40. № 3. С. 295.]
Malbrunot P., Vidal D., Vermesse J., Chahine R., Bose T.K. Adsorption Measurements of Argon, Neon, Krypton, Nitrogen, and Methane on Activated Carbon up to 650 MPa // Langmuir. 1992. V. 8. № 2. P. 577.
Потапов С.В., Фомкин А.А., Синицын В.А., Школин А.В. Адсорбция криптона на микропористых адсорбентах в области повышенных давлений // Физикохим. поверхн. защ. матер. 2010. Т. 46. № 6. С. 569.
Bazan R.E., Bastos-Neto M., Moeller A., Dreisbach F., Staudt R. Adsorption equilibria of O2, Ar, Kr and Xe on activated carbon and zeolites: single component and mixture data // Adsorption. 2011. V. 17. P. 371.
Кузнецова Т.А., Толмачёв А.М., Крюченкова Н.Г., Фирсов Д.А., Фомкин А.А. Термодинамика адсорбции криптона, ксенона, азота и кислорода на микропористом активном угле при температурах выше критических // Физикохим. поверхн. защ. матер. 2013. Т. 49. № 4. С. 339.
Kitani S., Takada J. Adsorption of Krypton and Xenon on Various Adsorbents // J. Nucl. Sci. Technol. 1965. V. 2. № 2. P. 51.
Thommes M., Kaneko K., Neimark V., Olivier P., Rodriguez-Reinoso F., Rouquerol J., S. W. Sing K. Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution (IUPAC Technical Report) // Pure Appl. Chem. 2015. V. 87. № 9–10. P. 1051.
Landers J., Gor G.Y., Neimark A.V. Density functional theory methods for characterization of porous materials // Colloids Surf., A. 2013. V. 437. P. 3.
Puziy A.M., Poddubnaya O.I., Gawdzik B., Sobiesiak M. Comparison of heterogeneous pore models QSDFT and 2D-NLDFT and computer programs ASiQwin and SAIEUS for calculation of pore size distribution // Adsorption. 2016. V. 22. P. 459.
Макаревич Н.А., Богданович Н.И. Теоретические основы адсорбции. Архангельск: САФУ, 2015.
Гаврилова Н.Н., Назаров В.В. Анализ пористой структуры на основе адсорбционных данных. М.: РХТУ им. Д.И. Менделеева, 2015.
Derrah R.I., Ruthven D.M. Sorption of the Inert Gases (Ar, Kr, and Xe) in Type A Zeolites // Can. J. Chem. 1975. V. 53. P. 996.
Гмурман В.Е. Теория вероятностей и математическая статистика. Учебное пособие для вузов. М.: Высшая школа, 2004.
Kwiatkowski M., Broniek E. An analysis of the porous structure of activated carbons obtained from hazelnut shells by various physical and chemical methods of activation // Colloids Surf., A. 2017. V. 529. P. 443.
Kwiatkowski M., Fierro V., Celzard A. Confrontation of various adsorption models for assessing the porous structure of activated carbons // Adsorption. 2019. V. 25. P. 1673.
Kwiatkowski M., Broniek E. An Evaluation of the Reliability of the Results Obtained by the LBET, QSDFT, BET, and DR Methods for the Analysis of the Porous Structure of Activated Carbons // Materials. 2020. V. 13. № 18. P. 3929.
Cimino R.T., Kowalczyk P., Ravikovitch P.I., Neimark A.V. Determination of Isosteric Heat of Adsorption by Quenched Solid Density Functional Theory // Langmuir. 2017. V. 33. № 8. P. 1769.
Nuhnena A., Janiak C. A practical guide to calculate the isosteric heat/enthalpy of adsorption via adsorption isotherms in metal–organic frameworks, MOFs // Dalton Trans. 2020. V. 49. P. 10295.
Матвейкин В.Г., Погонин В.А., Путин С.Б., Скворцов С.А. Математическое моделирование и управление процессом короткоцикловой безнагревной адсорбции. М.: Машиностроение-1, 2007.
Design of Off-Gas and Air Cleaning Systems at Nuclear Power Plants (Tech. Rep. Series). Vienna: IAEA, 1987.
Nakhutin I.E., Ochkin D.V., Pis’man B.Y., Babenko E.A., Krasikov A.N., Tret’yak S.A., Khlyustova T.M. A radiochromatographical system for scrubbing gas to remove radioactive contamination in atomic power plants with water-moderator water-cooled reactors // Sov. At. Energy. 1978. V. 44. № 3. P. 281. [Нахутин И.Е., Очкин Д.В., Письман Б.Я., Бабенко Б.Я., Красиков А.Н., Третьяк С.А., Хлюстова Т.М. Радиохроматографическая система очистки газов от радиоактивных загрязнений на АЭС с ВВЭР // Ат. энерг. 1978. Т. 44. № 3. С. 251.]
Mirza S.M., Iqbal J.M., Mirza N.M. Effect of flow rate transients on fission product activity in primary coolant of PWRs // Prog. Nucl. Energy. 2007. V. 49. № 2. P. 120.
Никтенко Е.И., Мурашова Е.Л. Состав выбросов инертных радиоактивных газов, образующихся в процессе эксплуатации промышленной реакторной установки на ФГУП “ПО "МАЯК” // Биосферная совместимость: чел., рег., технол. 2017. № 2 (18). С. 90.
IAEA Nuclear Data Services (2021) Live chart of nuclides: Nuclear structure and decay data. https://www-nds.iaea.org/relnsd/vcharthtml/VChartHTML.html. Cited 3 February 2021.
Дополнительные материалы отсутствуют.
Инструменты
Теоретические основы химической технологии