Известия РАН. Серия физическая, 2020, T. 84, № 10, стр. 1509-1513
Расчетные исследования радиотоксичности облученного топлива ВВЭР-1200 и запасов активности при долгосрочном обращении с отработавшим ядерным топливом Белорусской АЭС
Ю. А. Корчева 1, *, Н. В. Горбачева 1, Н. Д. Кузьмина 1, Н. В. Кулич 1, А. М. Петровский 1
1 Государственное научное учреждение “Объединенный институт энергетических и ядерных исследований – Сосны” Национальной академии наук Беларуси
Минск, Республика Беларусь
* E-mail: julia.korchova@sosny.bas-net.by
Поступила в редакцию 11.05.2020
После доработки 02.06.2020
Принята к публикации 26.06.2020
Аннотация
Представлены результаты расчетов активностей продуктов деления и актинидов для средней глубины выгорания топлива Белорусской АЭС (55.6 МВт ⋅ сут ⋅ т–1 U). Получены темпы снижения суммарной активности отработавшего ядерного топлива (ОЯТ) в зависимости от времени выдержки. На основании проведенных расчетных исследований характеристик радиотоксичности ОЯТ ВВЭР-1200 определены перечни радиационно-опасных радионуклидов, определяющих потенциальную опасность на различных фазах обращения с ОЯТ Белорусской АЭС. Полученные результаты имеют значение для принятия решений по экологически приемлемому варианту обращения с отработавшим ядерным топливом, представленные в Стратегии обращения с ОЯТ Белорусской АЭС.
ВВЕДЕНИЕ
В настоящее время в Республике Беларусь завершается этап сооружения Белорусской АЭС и идет подготовка к вводу ее в эксплуатацию в составе двух энергоблоков АЭС-2006 с ВВЭР-1200. Объединенный институт энергетических и ядерных исследований – Сосны с участием АО “Техснабэкспорт”, как базового отраслевого центра российских технологий по обращению с отработавшим ядерным топливом (ОЯТ), разработал Проект Стратегии обращения с отработавшим ядерным топливом Белорусской АЭС (далее – Проект Стратегии), отражающий национальную политику в этой области. В соответствии с Законом Республики Беларусь “О государственной экологической экспертизе, стратегической экологической оценке и оценке воздействия на окружающую среду” [1] эксплуатирующей организацией проведена процедура по стратегической экологической оценке (далее – СЭО) Проекта Стратегии. Важным элементом процедуры по СЭО для обоснования экологически приемлемого варианта обращения с ОЯТ Белорусской АЭС стали разработка Экологического доклада по СЭО Проекта Стратеги и его обсуждение с общественностью.
Принимая во внимание современную мировую практику обращения с ОЯТ атомных электростанций, вариантами Проекта Стратегии предусматриваются следующие основные фазы [2, 3]:
– промежуточное хранение ОЯТ после выгрузки его из приреакторных бассейнов выдержки до отправки на переработку;
– переработка ОЯТ на предприятиях Российской Федерации;
– возврат и хранение продуктов переработки ОЯТ с учетом принципа радиационного эквивалента;
– захоронение продуктов переработки ОЯТ на территории Республики Беларусь.
По предварительным данным на Белорусской АЭС будет реализован четырехгодичный топливный цикл с одной перегрузкой в 12 мес, средняя расчетная глубина выгорания тепловыделяющих сборок (далее – ТВС) в установившемся режиме перегрузок ядерного топлива равна 55.6 МВт ⋅ сут ⋅ т–1 U. Общее количество ТВС с различной степенью обогащения, включая топливо нового типа обогащением 4.95% по 235U, выгружаемых за 60 лет эксплуатации двух энергоблоков Белорусской АЭС, составит ориентировочно Ntot = 5300 единиц, общая масса достигнет Mtot = 2500 т [3].
Целью настоящей работы являются прогнозные исследования запасов активности и характеристик радиотоксичности ОЯТ в объектах хранения и захоронения при долговременном обращении с ОЯТ Белорусской АЭС.
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ФОРМИРОВАНИЯ ЗАПАСОВ РАДИОАКТИВНОСТИ В ХРАНИЛИЩЕ ОЯТ АЭС
Для оценки запасов активности продуктов деления, в условиях имеющейся на текущий момент неопределенности данных по режиму облучения, техническим характеристикам отработавших тепловыделяющих сборок (далее – ОТВС), продолжительности фаз обращения с ОЯТ разработана вероятностная модель формирования запасов и прогнозирования темпов снижения радиоактивности ОЯТ в объектах ядерного топливного цикла Белорусской АЭС. Для разработки модели использована методология анализа дискретных вероятностных процессов [4]. Первичным событием моделируемого процесса является перемещение ОТВС из активной зоны в бассейн выдержки энергоблоков № 1 или 2, выгорание которой достигло требуемой глубины выгорания топлива. Принимая во внимание прогнозируемое количество ОТВС, отличающихся начальным обогащением k и глубиной выгорания $B_{n}^{k}$, недостаток на текущий момент информации по срокам выгрузки каждой ОТВС, дает нам основание рассматривать процедуру перемещения n-й ОТВС случайным процессом на интервале времени жизни АЭС $\left[ {{{t}_{A}},{{t}_{B}}} \right]$. Тогда параметр, характеризующий длительность выдержки ${\theta } = t - {{t}_{n}}$ n-й ОТВС и, соответственно, активность r – радионуклида $A_{n}^{r}({{t}_{p}}--{{t}_{n}})$ в момент времени, характерный для каждой фазы обращения t = tp, также является случайной величиной.
Суммарная активность основных долгоживущих радиационно-опасных продуктов деления и актиноидов для всей массы Mtot ОЯТ в момент t = tp определяется по следующей формуле:
(1)
$\begin{gathered} A\left( {{{t}_{p}}} \right) = \sum\limits_{n = 1}^N {\sum\limits_{r = 1}^R {A_{n}^{r}} } \left( {{{t}_{p}} - {{t}_{n}}} \right) = \\ = \sum\limits_{n = 1}^N {\sum\limits_{r = 1}^R {\left[ {A_{n}^{r}\left( {{{t}_{p}} - {{t}_{n}}} \right) \cdot \exp \left( { - \frac{{\ln 2}}{{T_{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}^{r}}}\left( {{{t}_{p}} - {{t}_{s}}} \right)} \right)} \right],} } \\ \end{gathered} $Задача оценки суммарной активности на каждой фазе обращения с ОЯТ сводится к простой, но трудоемкой процедуре суммирования, в которой параметр tn является случайной величиной, заданной на интервале [tA, tB] его возможных значений. С помощью процедур Монте-Карло, реализованных в разработанном для этих целей программном средстве имитационного моделирования CUB [5], реальная последовательность моментов выгрузки ОТВС моделируется равномерной выборкой случайных значений на интервале $\left[ {{{t}_{A}},{{t}_{B}}} \right]$. С помощью уравнения (1) строится выборка случайных значений активности $A\left( {{{\theta }_{n}}} \right)$. Размер выборки принимает значение N, равное суммарному числу ОТВС, наработанных на энергоблоках № 1 и 2 Белорусской АЭС за прогнозный срок эксплуатации.
ОЦЕНКА ЗАПАСОВ АКТИВНОСТИ НА РАЗНЫХ ФАЗАХ ОБРАЩЕНИЯ С ОЯТ БЕЛОРУССКОЙ АЭС
В качестве исходных данных использованы методические материалы документа РБ-093-14 “Радиационные и теплофизические характеристики отработавшего ядерного топлива водо-водяных энергетических реакторов и реакторов большой мощности канальных” по удельным активностям продуктов деления и актинидов в ОТВС с начальным обогащением 4.81% по 235U в зависимости от глубины выгорания в реакторе ВВЭР-1000 [6]. Для проведения расчетных исследований авторами подготовлены базы данных в среде Microsoft Excel по активностям продуктов деления и актинидов (в расчете на 1 ОТВС) для топлива с начальным обогащением 4.81% и средней глубиной выгорания топлива Белорусской АЭС (табл. 1). Оцененная неопределенность результатов расчета концентраций (удельных активностей) актинидов и продуктов деления при использовании аппроксимирующей зависимости лежит в диапазоне 1–64% [6], что коррелирует с неопределенностью оценки глубины выгорания топлива в ТВС расчетным либо экспериментальным методами. Выполненный авторами с помощью программного средства Serpent [7] численный расчет удельных активностей продуктов деления и актинидов для топлива с начальным обогащением 4.95% и выгоранием 55.6 МВт ⋅ сут ⋅ кг–1 U показал, что отличие от данных, представленных в табл. 1, составляет 7–15% [8], что для целей данной работы является вполне приемлемым.
Таблица 1.
Радионуклид | Активность на 1 ТВС, Бк ‧ ТВС–1 | T1/2, с | Радионуклид | Активность на 1 ТВС, Бк ‧ ТВС–1 | T1/2, с |
---|---|---|---|---|---|
144Ce | 2.53 ‧ 1016 | 2.46 ‧ 1007 | 93Zr | 3.19 ‧ 1010 | 4.83 ‧ 1013 |
134Cs | 4.73 ‧ 1015 | 6.51 ‧ 1007 | 125m Te | 4.11 ‧ 1013 | 5.01 ‧ 1006 |
135Cs | 1.17 ‧ 1010 | 7.26 ‧ 1013 | 237Np | 1.63 ‧ 1010 | 6.75 ‧ 1013 |
137Cs | 2.83 ‧ 1015 | 9.47 ‧ 1008 | 238Pu | 1.79 ‧ 1014 | 2.77 ‧ 1009 |
154Eu | 2.16 ‧ 1014 | 2.71 ‧ 1008 | 239Pu | 1.04 ‧ 1013 | 7.61 ‧ 1011 |
155Eu | 8.97 ‧ 1013 | 1.48 ‧ 1008 | 240Pu | 1.87 ‧ 1013 | 2.07 ‧ 1011 |
3H | 3.37 ‧ 1011 | 3.89 ‧ 1008 | 241Pu | 5.22 ‧ 1015 | 4.53 ‧ 1008 |
85Kr | 1.57 ‧ 1014 | 3.38 ‧ 1008 | 242Pu | 1.10 ‧ 1011 | 1.18 ‧ 1013 |
107Pd | 2.41 ‧ 1009 | 2.05 ‧ 1014 | 234U | 3.31 ‧ 1010 | 7.75 ‧ 1012 |
106Ru | 1.01 ‧ 1016 | 3.21 ‧ 1007 | 235U | 4.64 ‧ 1008 | 2.22 ‧ 1016 |
125Sb | 1.86 ‧ 1014 | 8.61 ‧ 1007 | 236U | 1.19 ‧ 1010 | 7.39 ‧ 1014 |
79Se | 9.87 ‧ 1009 | 1.04 ‧ 1012 | 238U | 8.48 ‧ 1009 | 1.41 ‧ 1017 |
151Sm | 9.77 ‧ 1012 | 2.84 ‧ 1009 | 241Am | 7.00 ‧ 1012 | 1.37 ‧ 1010 |
121m Sn | 3.03 ‧ 1011 | 1.74 ‧ 1009 | 243Am | 1.60 ‧ 1012 | 2.33 ‧ 1011 |
126Sn | 1.30 ‧ 1010 | 3.15 ‧ 1012 | 242Cm | 2.29 ‧ 1015 | 1.41 ‧ 1007 |
90Sr | 1.31 ‧ 1015 | 8.88 ‧ 1008 | 244Cm | 2.81 ‧ 1014 | 5.41 ‧ 1008 |
99Tc | 2.56 ‧ 1011 | 6.66 ‧ 1012 | 245Cm | 1.00 ‧ 1013 | 2.68 ‧ 1011 |
Моделирование запасов активности на каждой фазе обращения с ОЯТ по каждому из радиационно-опасных радионуклидов, а также суммарной активности, проведено с использованием программного средства имитационного моделирования CUB [5]. Результаты оценки суммарной активности продуктов деления и актинидов отработавшего топлива энергоблоков № 1 и 2 (5300 ед. ОТВС) в зависимости от времени выдержки приведены в табл. 2. Расчет показал, что после динамичного снижения в течение первых сотен лет удельная активность меняется медленно и через 1000 лет достигает уровня порядка 1.77 ⋅ 1017 Бк. Оцененная суммарная активность ОЯТ после выдержки 10 000 лет составит величину 3.88 ⋅ 1016 Бк.
РЕЗУЛЬТАТЫ РАСЧЕТОВ ХАРАКТЕРИСТИК РАДИОТОКСИЧНОСТИ ОЯТ ВВЭР-1200
Для оценки потенциальной опасности ОЯТ Белорусской АЭС на различных фазах обращения наряду с радиационными характеристиками, важно располагать сведениями по характеристикам радиотоксичности радионуклидов на различных временных фазах обращения [9, 10]. В Республике Беларусь для целей нормирования поступления радионуклидов в организм человека пероральным путем нормативно установлены уровни вмешательства по воде [11]. По определению радиотоксичность радионуклида по воде определяет объем воды, который необходим для разбавления активности i-го радионуклида, содержащегося в 1 т уранового топлива, для снижения до приемлемого уровня.
После выдержки в приреактором бассейне облученное топливо перемещается на хранение в промежуточное хранилище ОЯТ на срок не менее 50 лет, после которого, вследствие радиоактивного распада продуктов деления, радиационные характеристики ОЯТ снижаются до уровней, допускающих дальнейшее обращение на установках по переработке.
На рис. 1, 2 представлены темпы снижения радиотоксичности продуктов деления и актинидов на этапах хранения от 10 до 10 000 лет при долговременном хранении ОЯТ. Расчеты выполнены на основании данных, представленных в табл. 1.
На этапе промежуточного хранения топлива со средним выгоранием 55.6 МВт · сут · т–1 U радиотоксичность по воде определяется вкладом продуктов деления 106Ru, 144Ce, 90Sr, 137Cs, 134Сs, 147Pm, 154Eu. Радиотоксичность по воде продуктов деления за 100 лет снижается от 1018 до 1015 кг ⋅ т–1. Преобладающий вклад в суммарную радиотоксичность ОЯТ на временном этапе до 100 лет вносят изотопы актинидов 238Pu, 239Pu, 240Pu, 241Pu, 242Cm и 244Cm. Наблюдается рост удельной активности и, соответственно, радиотоксичности 241Am, что связано с его накоплением за счет распада материнского нуклида 241Pu. Радиотоксичность актинидов по воде снижается от 1019 до 1017 кг ⋅ т–1, что на два порядка превышает вклад в суммарную радиотоксичность продуктов деления. Оценка вкладов в суммарную радиотоксичность по воде ОЯТ показала, что на этапе хранения до 100 лет основной вклад дают радионуклиды 90Sr, 106Ru, 137Cs, 134Cs, 144Ce, 147Pm, 154Eu, 238Pu, 240Pu, 241Pu, 244Cm, а также 241Am.
На этапе долговременного хранения ОЯТ от 1000 до 10 тыс. лет потенциальная опасность определяется основными долгоживущими продуктами деления и актинидами, из которых основной вклад в радиотоксичность топлива вносит 99Tc, и в меньшей степени 93Zr, 79Se, 239Pu, 240Pu, 242Pu и 243Am. Радиотоксичность по воде продуктов деления снижается до 1013 кг ⋅ т–1. Радиотоксичность актинидов по воде снижается до 1015 кг ⋅ т–1. В целом, вклад актинидов на два порядка превышает вклад в суммарную радиотоксичность продуктов деления.
ЗАКЛЮЧЕНИЕ
Разработанная вероятностная модель, реализованная в виде программного средства CUB, обеспечивает прогнозирование запасов активности ОЯТ по перечню наиболее опасных продуктов деления и актиноидов на разных фазах обращения с ОЯТ АЭС с ВВЭР-1200.
Результаты моделирования показали, что вследствие радиоактивного распада продуктов деления на этапе хранения в промежуточном хранилище ОЯТ (в течение 100 лет) суммарная активность ОЯТ Белорусской АЭС снижается до уровня 4.87 ⋅ 1018 Бк. На этапе длительного хранения суммарная активность не превышает 1016 Бк.
На основании проведенных расчетных исследований характеристик радиотоксичности ОЯТ ВВЭР-1200 определены перечни радиационно-опасных радионуклидов, определяющих потенциальную опасность на различных фазах обращения с ОЯТ Белорусской АЭС.
Результаты исследований могут быть использованы в поддержку обоснования экологически приемлемого варианта обращения с ОЯТ Белорусской АЭС. Надежность оценок может быть повышена при получении детальных сведений по характеристикам ядерного топливного цикла Белорусской АЭС.
Список литературы
https://kodeksy-by.com/zakon_rb_o_gosudarstvennoj_ekologicheskoj_ekspertize.htm.
Status and trends in spent fuel and radioactive waste management. IAEA nuclear energy series № NW-T-1.14. Vienna, 2018.
Экологический доклад по стратегической экологической оценке проекта Стратегии обращения с отработавшим ядерным топливом Белорусской атомной электростанции. Минск, 2018.
Вентцель Е.С. Исследование операций. М.: Наука, 1972. С. 552.
Горбачева Н.В. // Докл. 6-й Межд. науч.-техн. конф. “Обеспечение безопасности АЭС с ВВЭР” (Москва, 2009). С. 020.
Радиационные и теплофизические характеристики отработавшего ядерного топлива водо-водяных энергетических реакторов и реакторов большой мощности канальных РБ-093-14 // ЯРБ. 2014. № 4(74). С. 49.
Leppaenen J. Serpent, continuous-energy Monte Carlo reactor physics burnup calculation code. VTT technical research center of Finland, 2013.
Petrovskii A.M., Rudak E.A., Korbut T.N., Kravchenko M.O. // J. Phys. Conf. Ser. 2018. V. 1133. Art № 012009.
A basic toxicity classification of radionuclides. Technical reports series № 15. Vienna: IAEA, 1963. P. 39.
Korchova J.A., Harbachova N.V., Kuzmina N.D. // J. Phys. Conf. Ser. 2018. V. 1133. Art № 012025.
Критерии оценки радиационного воздействия: гигиенический норматив. Нац. реестр правовых актов РБ. 8/26850, 2013.
Дополнительные материалы отсутствуют.
Инструменты
Известия РАН. Серия физическая