Расплавы, 2023, № 2, стр. 113-121

Распределение радиоактивных элементов при шлаковом переплаве конструкционных материалов оболочек ТВЭЛов ВВЭР

О. Н. Будин a*, И. В. Кузнецов a, М. Ю. Каленова a, С. А. Красиков bc, А. С. Щепин a

a Ведущий научно-исследовательский институт химической технологии
Москва, Россия

b Институт металлургии Уральское отделение РАН
Екатеринбург, Россия

c Уральский государственный горный университет
Екатеринбург, Россия

* E-mail: o.n.budin@gmail.com

Поступила в редакцию 06.07.2022
После доработки 07.11.2022
Принята к публикации 20.11.2022

Полный текст (PDF)

Аннотация

На радиохимическом заводе ФГУП “ПО “Маяк” в процессе переработки отработавшего ядерного топлива (ОЯТ) образуются около 170 т/год металлических радиоактивных отходов (МРАО), преимущественно представленных фрагментами конструкционных материалов (КМ) оболочек ТВЭЛов, содержащих остаточные количества радионуклидов после растворения топлива, и отработавших тепловыделяющих сборок (ОТВС), отправляемых на складирование. Достичь компактизации и дезактивации МРАО возможно способом, основанном на шлаковом переплаве в индукционной печи с холодным тиглем. Для установления распределения актинидов и продуктов деления (ПД) проведено термодинамическое моделирование процесса шлакового переплава КМ оболочек ТВЭЛов и ОТВС реакторной установки ВВЭР-1000 и экспериментальная верификация полученных данных. Показано наиболее вероятное распределение актинидов и ПД по продуктам плавки. Большая часть кюрия и америция концентрируются в металле – 99 и 94 мас. % соответственно. Максимальное извлечение урана в шлаковую фазу в виде диоксида UO2 составляет до 40 мас. %. Распределение плутония в интервале температур 1500–2000°С по шлаковой (в виде оксидов PuO и PuO1.61) и металлической фазам происходит практически в равном соотношении. Установлено, что в составе газоаэрозольного потока, представленного ПД, преобладает до 99.78 мас. % цезия, содержание европия и америция составляет 0.05 и 0.17 мас. % соответственно.

Ключевые слова: термодинамическое моделирование, шлаковый переплав, оболочки ТВЭЛов, флюс, плутоний, уран

ВВЕДЕНИЕ

Переработка отработавшего ядерного топлива на радиохимических производствах РФ приводит к накоплению МРАО, отнесенных ко 2-му классу, обусловленному наведенной активностью и содержанием следовых количеств топливных материалов (остаточного растворителя/промывочного раствора, загрязняющего поверхность, и ПД, диффундирующих в толщу металла) [15]. На площадке ФГУП “ПО “Маяк” существующая практика обращения с ними подразумевает пересыпание противопожарным наполнителем в массовом соотношении 3 : 1 и складирование навалом [5]. Хранение материалов с нефиксированным загрязнением создает угрозу выхода радионуклидов в окружающую среду, а избыточный объем некондиционированного продукта ведет к неэффективному использованию существующих хранилищ.

Кандидатным способом переработки МРАО является шлаковый переплав в индукционной печи с холодным тиглем (ИПХТ), обеспечивающий максимальную степень компактизации, минимизацию вторичных радиоактивных отходов (РАО) в виде отработанных тигельных материалов, изложниц и прочих, возможность возвращения и повторного использования актинидов и ПД, извлекаемых в шлак [69]. При этом оставшаяся часть делящихся материалов, содержавшихся в хрупких поверхностных оксидных пленках толщиной до 100 мкм, восстанавливается, равномерно распределяется и фиксируется в устойчивой металлической матрице [10].

В настоящей работе проведено исследование распределения актинидов и ПД по продуктам плавки с использованием термодинамического (ТД) моделирования в зависимости от состава металлической фракции и температуры процесса с последующей экспериментальной верификацией части расчетных данных.

ТЕРМОДИНАМИЧЕСКОЕ МОДЕЛИРОВАНИЕ

Объектом исследований являлись оболочки ТВЭЛов ОТВС ВВЭР-1000 блока № 1 Балаковской АЭС с выгоранием 47.28 ГВт · сут/тU. Удельные содержания актинидов и ПД, полученных в 2021 г. специалистами АО “Радиевый институт им. В.Г. Хлопина” в рамках реализации проекта Единого отраслевого тематического плана Госкорпорации “Росатом”, и принятых в качестве исходных данных для ТД моделирования, представлены в табл. 1.

Таблица 1.  

Содержание актинидов и ПД в оболочках ТВЭЛов

Элемент Форма нахождения Содержание, мкг/г оболочек Содержание, % от массы оболочек
U U3O8 180 0.02
Pu PuO2 7 7.78 · 10–4
Am AmO2 0.3 3.33 · 10–5
Cm Cm2O3 0.044 4.89 · 10–6
Cs Cs2O 6.2 6.89 · 10–4
Eu Eu2O3 0.014 1.56 · 10–6

По данным [11] и результатам анализа номенклатуры, образующихся МРАО при разборке ОТВС ВВЭР-1000, определены соотношения КМ, представленных нержавеющей сталью и циркониевым сплавом Э110 (табл. 2). Использование смеси фрагментов различного состава позволяет снизить температуру образующегося расплава.

Таблица 2.  

Составы имитаторов МРАО

Состав Массовое соотношение Описание МРАО
№ 1 93.3% Э110– 6.7% нержавеющая сталь КМ ТВЭЛов (состав, имитирующий смесь фрагментов оболочек ТВЭЛов с фиксирующими пружинами)
№ 2 82.2% Э110– 17.8% нержавеющая сталь КМ ОТВС (состав, имитирующий смесь фрагментов центральной части, головки и хвостовика ОТВС ВВЭР-1000)

Основной маркой нержавеющей стали, использующейся при производстве тепловыделяющей сборки ВВЭР, является 08Х18Н10Т либо 12Х18Н10Т [11]. В качестве металлической фазы при ТД расчетах принята марка нержавеющей стали 12Х18Н10Т, что обусловлено наличием соответствующих характеристик в базах данных.

В качестве флюса рассмотрена оксидная композиция состава, мас. %: 19.5CaO–2.5MgO–78Al2O3, не содержащая соединений кремния и бора, преимуществом которой является возможность гидрометаллургического извлечения актинидов, возвращаемых в ядерный топливный цикл [12].

Исследование распределения элементов по продуктам плавки при ТД моделировании взаимодействия КМ оболочек ТВЭЛов и ОТВС с флюсом выполнялось с использованием программы HSC-6.1 [13]. Работа программы основана на принципе минимизации свободной энергии Гиббса [12, 15] исследуемой замкнутой системы и позволяет рассчитать равновесный состав продуктов взаимодействия в многокомпонентных системах. Для проведения расчетов в программе используется простейшая модель идеальных растворов, где выполняется закон Рауля и коэффициент активности равен 1. В расчетах приняты следующие условия: масса сплава – 100 кг, исследуемый интервал температур 1000–2000°С, массовое соотношение металла к флюсу – 10 : 1. Интервал температур является наиболее интересным в части исследования вероятного протекания твердофазных и жидкофазных процессов. Верхняя граница обусловлена необходимостью расплавления флюса и его перегрева на 200–300°С выше температуры плавления. Учитывая сродство циркония к кислороду, принято, что процесс осуществляется в атмосфере 100% аргона при атмосферном давлении.

Результаты ТД моделирования распределения актинидов и ПД при взаимодействии металлической садки с флюсом представлены на рис. 1 и 2.

Рис. 1.

Распределение актинидов и ПД при взаимодействии сплава состава № 1 и оксидного флюса, мас. %: 19.5CaO–2.5MgO–78Al2O3. (а) Металлическая фаза; (б) шлаковая фаза; (в) газовая фаза.

Рис. 2.

Распределение актинидов и ПД при взаимодействии сплава состава № 2 и оксидного флюса, мас. %: 19.5CaO–2.5MgO–78Al2O3. (а) Металлическая фаза; (б) шлаковая фаза; (в) газовая фаза.

Установлено, что уран частично должен извлекаться в шлак (38–40 мас. %) вне зависимости от состава металлической фазы в виде UO2, а также образовывать соединения с фосфором (компонент стали 12X18Н10Т) – UP, UP2, U3P4, разлагающиеся при температурах свыше 1500°С.

Плутоний распределяется по шлаковой и металлической фазам практически в равном соотношении в интервале температур 1500–2000°С. При плавлении КМ ТВЭЛов и ОТВС извлечение плутония может происходить в виде оксидов PuO и PuO1.61.

Америций претерпевает “металлизацию” и частично может возгоняться в элементарном виде (до 6 мас. % при 2000°С) вне зависимости от состава металлической фазы, и менее 0.05 мас. % переходить в шлак как AmO2 и Am2O3. Аналогичное поведение наблюдается и для цезия, однако этому элементу свойственны большие упругости паров в тех же условиях. При 2000°С цезий до 89 мас. % переходит в газовую фазу, до 12 мас. % – в металлическую.

Поведение европия и кюрия практически идентично для обоих сплавов. До 40 мас. % европия распределяется в металлической фазе, около 32 мас. % ассимилируется шлаком как монооксид EuO, до 25.7 мас. % переходит в газовую фазу в элементарном виде. Кюрий полностью “металлизируется” и в шлак не извлекается, частично переходя в газообразное состояние (не более 0.02 мас. %).

При анализе состава газовой фазы установлено, что основным компонентом является цезий (около 99.78 мас. %), что закономерно вследствие относительно высокой упругости паров этого элемента [14].

ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ

Экспериментальная верификация проводилась при изучении взаимодействия имитаторов МРАО с флюсующей композицией в вакуумной печи “Вега-31” (ООО “ЭВТ”, Россия). Фрагменты сплавов заданного состава, полученных в инертной атмосфере на лабораторной электродуговой печи 5SA (Centorr Vacuum Industries, США), флюс и U3O8 загружались в тигли объемом 10 см3 из оксида циркония, модифицированного Y2O3. Радионуклидная составляющая эксперимента представлена только U3O8, что обусловлено крайне низким (менее 10–4% от массы металла) содержанием остальных актинидов и ПД, находящегося ниже предела обнаружения аналитических приборов. Количество материалов в каждом из экспериментов составляло: сплава 5–6 г; флюса – 10% от массы металла; U3O8 – 0.02% от массы металла. Процесс осуществлялся в атмосфере аргона высокой чистоты при температуре 1900°С и изотермической выдержке 60 мин. Скорость нагрева печи 15°С/мин. Заданная продолжительность обусловлена установлением равновесных содержаний радионуклидов в шлаке и металле [1518]. Снижение температуры с 1900°С до комнатной осуществлялось в свободном режиме с отключенными нагревателями.

Элементный состав флюсов и шлаков исследовался на энергодисперсионном рентгенофлюоресцентном спектрометре ARL PERFORM'X Sequential X-Ray (Thermo Fisher Scientific Inc., США). Содержание урана в металле определялось методом атомно-эмисионной спектрометрии с индуктивно связанной плазмой на спектрометре 700-ES (Agilent, Австралия).

На основе результатов элементного анализа проб рассчитывалось распределение урана по продуктам плавки (табл. 3), приведенное к 100 мас. %.

Таблица 3.  

Распределение урана в образцах, полученных в ходе эксперимента

Металлическая садка Распределение урана, мас. %
в шлаке в металле потери
82.8% Э110–17.2% 12Х18Н10Т 37.33 62.66 0.01
93.3% Э110–6.7% 12Х18Н10Т 36.11 64.87 0.02

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

Согласно результатам ТД моделирования вне зависимости от состава металлической фракции и температуры, поведение актинидов и ПД в процессе переплава практически идентичное. Из полученных данных следует, что большая часть актинидов (более 50 мас. %) в процессе переплава МРАО с оксидным флюсом концентрируется в металле. Переход актинидов и ПД в продукты переплава возможен в результате их взаимодействия с компонентами сплава Э110 (цирконий) и стали 12X18Н10Т (кремний, фосфор, углерод). В прямом направлении для интервала температур 1000–2000°С [13] возможно протекание реакций (табл. 4), характеризующихся отрицательными значениями энергии Гиббса.

Таблица 4.  

Реакции, характеризующиеся отрицательным значением энергии Гиббса

Реакция Температурный интервал возможного протекания реакций, °C Значение изменения энергии Гиббса при Т = 1900°C (2000°C), кДж/моль
1 U3O8 + 4Zr = 3U + 4ZrO2 1000–2000 –641.9
2 PuO2 + Zr = Pu + ZrO2 1000–2000 –43.6
3 U3O8 + 8C = 3U + 8CO (г) 1800–2000 –253.0
4 U3O8 + 5C = 3UO + 5CO(г) 1500–2000 –357.9
5 U3O8 + 2C = 3UO2 + 2CO(г) 1000–2000 –503.7
6 2Cs2O = 4Cs(г) + O2(г) 1700–2000 –112.1
7 5Cs2O + 2P = 10Cs(г) + P2O5 1000–2000 –838.6
8 2Cs2O + P = 4Cs(г) + PO2 1800–2000 –61.1
9 Cs2O + C = 2Cs(г) + CO(г) 1000–2000 –356.2
10 PuO2 + C = PuO + CO(г) 1800–2000 –18.3
11 AmO2 + 2C = Am(г) + 2CO(г) 1900–2000 –7.7
12 AmO2 + 2C = Am + 2CO(г) 1900–2000 –25.3
13 2AmO2 + C = Am2O3 + CO(г) 1000–2000 –246.7
14 4AmO2 + C = 2Am2O3 + CO2(г) 1000–2000 –288.9
15 Cm2O3 + 3C = 2Cm + 3CO(г) 1800–2000 –135.2
16 Cm2O3 + 3Si = 2Cm + 3SiO(г) 1900–2000 –23.5
17 Eu2O3 + C = 2EuO + CO(г) 1700–2000 –45.5
18 Eu2O3 + 3C = 2Eu(г) + 3CO(г) 2000 (–20.7)
19 Eu2O3 + Si = 2EuO + SiO(г) 1900–2000 –8.3
20 2Eu2O3 + Si = 4EuO + SiO2 1500–2000 –8.2
21 AmO2 + Zr = Am + ZrO2 100–2000 –122.5
22 AmO2 + Zr = Am(г) + ZrO2 800–2000 –104.9
23 2Cm2O3 + 3Zr = 4Cm + 3ZrO2 600–2000 –562.1
24 2Cm2O3 + 3Zr = 4Cm(г) + 3ZrO2 1900–2000 –45.7
25 2Eu2O3 + 3Zr = 4Eu + 3ZrO2   100–2000 –105.2
26 2Eu2O3 + 3Zr = 4Eu(г) + 3ZrO2 1400–2000 –196.8

Как видно из табл. 4, “металлизация” америция возможна через восстановление AmO2 углеродом (11), (12), с учетом промежуточных стадий (13), (14). Аналогично может протекать восстановление кюрия (15), (16) и европия (17)–(20). Однако, более вероятно металлотермическое восстановление перечисленных металлов при взаимодействии их оксидов с цирконием (21)–(26), содержание которого в металлических отходах оболочек ТВЭЛов около 90 мас. % (табл. 2 и 3). Поэтому логично допустить вероятность протекания таких металлотермических реакций в процессе переплава, как для изначально присутствующих изотопов-загрязнителей (актинидов и некоторых продуктов деления), так и вносимых извне, образующихся при переработке ОЯТ РУ других компонентов (благородные металлы, Tc, Se, Te и др.).

Таким образом, на основе результатов ТД расчетов проведена предварительная оценка гипотетического распределения ПД и ДМ между фазами. Так как конкретные величины энергии Гиббса указывают лишь на возможность протекания реакций и никак не отражают скорости их протекания. Для подтверждения частично установленного взаимодействия ПД и ДМ со шлаком и сталью проведена их экспериментальная проверка. Выявление кинетических параметров вышеуказанных реакций затруднено в связи с крайне низким содержанием большинства радионуклидов в составе металлических оболочек ТВЭЛов.

Однако, сопоставляя экспериментальные значения распределения урана и результаты моделирования, можно заключить, что верификация ТД модели свидетельствует о согласии расчетных и экспериментальных данных и позволяет предположить, что распределение актинидов и ПД по металлической, шлаковой и газовой фазам описывается достаточно достоверно.

ЗАКЛЮЧЕНИЕ

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

Работа выполнена в рамках проекта ЕОТП – ТЦПМ-263 Госкорпорации “Росатом” (НИОКР “Разработка технологии обращения с циркониевыми оболочками ТВЭЛов, образующимися при переработке ОЯТ ВВЭР”). Работа выполнена с использованием оборудования ЦКП “Урал-М”.

Авторы выражают благодарность старшему научному сотруднику Е.М. Жилиной и младшим научным сотрудникам Т.В. Осинкиной и А.С. Русских лаборатории редких тугоплавких металлов Института металлургии УрО РАН за их вклад в выполнение исследований в данной работе.

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

  1. Piro M.H.A. (ed.). Advances in Nuclear Fuel Chemistry. Woodhead Publishing, 2020.

  2. Кащеев В.А., Шадрин А.Ю., Рыкованов Г.Н., Дырда Н.Д., Макеева И.Р., Хмельницкий Д.В., Алексеев П.Н. Объем радиоактивных отходов от переработки облученного ядерного топлива ВВЭР-1000 и варианты фракционирования //Атомная энергия. 2019. 127. № 2. С. 82–87.

  3. Kang K.H. Lee C.H., Jeon M.K., Han S.Y., Park G.I., Hwang S.-M. Characterization of cladding hull wastes from used nuclear fuels // Archives of Metallurgy and Materials. 2015. 60. № 2B. P. 1199–1203.

  4. Park G.I. Jeon M.K., Choi J.-H., Lee K.-R., Han S.Y., Kim I.T., Cho Y.-Z., Park H.-S. Recent progress in waste treatment technology for pyroprocessing at KAERI // Journal of Nuclear Fuel Cycle and Waste Technology (JNFCWT). 2019. 17. № 3. P. 279–298.

  5. Козлов П.В., Ремизов М.Б., Смелова Т.В., Сунцов Д.Ю., Шестоперов И.Н. Изучение процесса индукционно-шлаковой переплавки металлических ВАО от переработки ОЯТ в “холодном тигле” // Тезисы докладов. Тринадцатый Международный Уральский Семинар. “Радиационная физика металлов и сплавов”. Кыштым, 2019. С. 53.

  6. Кащеев В.А. Смелова Т.В., Мусатов Н.Д., Шестопёров И.Н., Сунцов Д.Ю., Тучкова А.И., Арсеенков Л.В. Оценка возможности извлечения делящихся материалов из конструкционных материалов ТВЭЛов при переработке СНУП ОЯТ // Вопросы атомной науки и техники. Сер. Материаловедение и новые материалы. 2016. № 4. С. 35–47.

  7. Готовчиков В.Т., Борзунов А.И., Середенко В.А. и др. Способ пирометаллургической переработки отходов, отработавших материалов и изделий. Дата публикации 2001.08.27. Патент РФ RU2172787.

  8. Каленова М.Ю., Кузнецов И.В., Щепин А.С. и др. Кондиционирование конструкционных материалов облученных ТВС методом индукционно-шлакового переплава в холодном тигле // Атомная энергия. 2018. 124. № 5. С. 273–278.

  9. Каленова М.Ю., Кузнецов И.В., Щепин А.С. Технология очистки конструкционных материалов ТВЭЛов методом индукционно-шлакового переплава в холодном тигле. текущее состояние и перспектива разработки // Вопросы атомной науки и техники. Серия: Материаловедение и новые материалы. 2017. 89. № 2. С. 71–80.

  10. Бычков С.И., Жирников Д.В., Алексеенко В.Н., Мацеля В.И. Регенерация металлических радиоактивных отходов радиохимических производств для целей повторного использования металлов и сплавов // Радиоактивные отходы. 2021. № 2. С. 33–38.

  11. Шмелев Драгунов Ю.Г., Денисов В.П., Васильченко И.Н. Активные зоны ВВЭР для атомных станций М.: ИКЦ “Академкнига”, 2004.

  12. Каленова М.Ю., Дмитриева А.В., Кузнецов И.В. и др., Очистка конструкционных материалов отработавших тепловыделяющих сборок от актинидных загрязнителей методом индукционно-шлаковой переплавки в холодном тигле // Тонкие химические технологии. 2016. 11. № 6. С. 83–90.

  13. Roine A. Outokumpu HSC Chemistry for Windows. Chemical Reaction and Equilibrium Software with Extensive Thermochemical Database. Pori: Outokumpu Research OY, 2006.

  14. Казенас Е.К., Цветков Ю.В., Термодинамика испарения оксидов. РАН, Ин-т металлургии и материаловедения им. А.А. Байкова. М.: URSS, 2008. С. 474.

  15. Лосицкий А.Ф., Ганза Н.А., Рождественский В.В. и др. Способ переработки металлических отходов, содержащих радионуклиды. Дата публикации 2000.11.20. Патент РФ RU 2 159 473.

  16. Константинов Е.А., Кижнеров Л.В., Кораблев Н.А. и др. Способ обработки металлических отходов, загрязненных радионуклидами. Дата публикации 2001.06.10. Патент РФ RU 2 168 780.

  17. Голубев А.А., Гудим Ю.А. Способ переработки металлических радиоактивных отходов и агрегат для его осуществления. Дата публикации 2009.01.09. Патент РФ RU 2 345 141.

  18. Петров Г.А., Суворов И.С., Соболев И.А. и др. Способ дезактивации радиоактивных металлических отходов. Дата публикации 2004.06.27. Патент РФ RU 2 231 843.

Дополнительные материалы отсутствуют.