CC BY
8
DOI: 10.21870/0131-3878-2022-31-4-21-33 УДК 621.039.76
Программный модуль оценки достижения радиологической эквивалентности (РОЗА-РАО)
Меняйло А.Н., Чекин С.Ю., Ловачёв С.С., Туманов К.А., Иванов В.К.
МРНЦ им. А.Ф. Цыба - филиал ФГБУ «НМИЦ радиологии» Минздрава России, Обнинск
Основной проблемой широкомасштабного развития ядерной энергетики в мире является безопасное захоронение накапливаемых ею радиоактивных отходов (РАО). Опасность РАО на момент их образования для здоровья человека намного превышает опасность исходного природного уранового сырья, после использования которого эти отходы образовались. За счёт естественных процессов радиоактивного распада опасность РАО со временем уменьшается. Опасность отходов также может быть уменьшена за счёт искусственной трансмутации или сжигания радионуклидов, входящих в их состав. Для урановых ядерных топливных циклов естественным контрольным уровнем опасности РАО является уровень опасности исходного уранового сырья, по достижении которого могут приниматься решения об окончательном захоронении РАО. В статье представлен метод оценки времени достижения радиологической эквивалентности массы накопленных РАО и использованной при этом исходной массы природного уранового сырья. Представлено также описание программного модуля «Радиологическое обеспечение защиты - радиоактивные отходы» (ПМ РОзА-РАО), реализующего этот метод. Под радиологической эквивалентностью понимается равенство для населения соответствующих радиационных рисков. Применяемая методология определяется современными рекомендациями Международной комиссии по радиологической защите (МКРЗ) и Научного комитета оОн по действию атомной радиации (НКДАР ООН) по расчёту радиационных рисков человека от внутреннего облучения на основе эквивалентных доз в органах и тканях. Определение эквивалентных доз облучения в органах и тканях основывалось на дозовых коэффициентах радиоактивных веществ, предоставляемых для широкого использования Американским агентством по защите окружающей среды (United States Environmental Protection Agency, или EPA). Расчёты с помощью ПМ РОЗА-РАО показали, что в рамках рассмотренного наиболее вероятного сценария развития двухкомпонентной ядерной энергетики России для РАО, накопленных к 2100 г., время наступления радиологической эквивалентности составляет менее 100 лет (98,8 года). С помощью разработанного программного модуля впервые в атомной отрасли могут определяться параметры переработки облучённого ядерного топлива (состав и активности радионуклидов для контролируемого хранения, сжигания и трансмутации) и обращения с РАО, позволяющие минимизировать потенциальные радиационные риски населения, что соответствует современной правоприменительной практике развитых стран в области радиационной защиты.
Ключевые слова: радиоактивные отходы, природное урановое сырьё, радиационная эквивалентность, радиологическая эквивалентность, радиационная безопасность, радиационный риск, радиобиологические эффекты, злокачественные новообразования, модели радиационного риска МКРЗ и НКДАР ООН, дозовые коэффициенты, внутреннее облучение, радионуклиды, ожидаемая эффективная доза, эквивалентная доза.
Введение
Основной проблемой широкомасштабного развития ядерной энергетики в мире является безопасное захоронение накапливаемых ею радиоактивных отходов (РАО) [1]. Опасность РАО на момент их образования для здоровья человека намного превышает опасность исходного природного уранового сырья, после использования которого эти отходы образовались. За счёт естественных процессов радиоактивного распада опасность РАО со временем уменьшается. Опасность отходов также может быть уменьшена за счёт искусственной трансмутации или сжигания радионуклидов, входящих в их состав. Для урановых ядерных топливных циклов естественным
Меняйло А.Н.* - вед. науч. сотр., к.б.н.; Чекин С.Ю. - зав. лаб.; Ловачёв С.С. - мл. науч. сотр.; Туманов К.А. - зав. лаб., к.б.н. Иванов В.К. - науч. руководитель НРЭР, Председатель РНКРЗ, чл.-корр. РАН, д.т.н. МРНЦ им. А.Ф. Цыба - филиал ФГБУ «НМИЦ радиологии» Минздрава России.
•Контакты: 249035, Калужская обл., Обнинск, ул. Королёва, 4. Тел.: (484) 399-32-81; e-mail: [email protected].
контрольным уровнем опасности РАО является уровень опасности исходного уранового сырья, по достижении которого могут приниматься решения об окончательном захоронении РАО.
В статье представлено описание разработанного метода оценки времени достижения радиологической эквивалентности массы накопленных ядерной энергетикой РАО и использованной при этом исходной массы природного уранового сырья. Приведено также описание специализированного программного модуля «Радиологическое обеспечение защиты - радиоактивные отходы» (ПМ РОЗА-РАО), реализующего этот метод.
Под радиологической эквивалентностью в данном случае понимается равенство для населения соответствующих радиационных рисков. Программный модуль сравнивает пожизненный атрибутивный риск (LAR, от англ. lifetime attributable risk) заболеваемости злокачественными новообразованиями (ЗНО) после однократного поступления в организм человека основных дозооб-разующих радионуклидов, входящих в состав природного уранового сырья, и соответствующие величины LAR от радионуклидов, входящих в состав накопленных РАО после различных временных периодов их выдержки. Радиационный риск накопленных РАО уменьшается с увеличением длительности их выдержки и через некоторый период времени достигает уровня радиационного риска исходного уранового сырья. Этот период времени называется далее временем достижения радиологической эквивалентности.
Ранее контрольный уровень опасности РАО определялся по показателям так называемой «потенциальной биологической опасности» (ПБО), которая выражалась в единицах ожидаемой эффективной дозы (ОЭД) [2]. Период времени, в течение которого ПБО РАО достигает уровня ПБО исходного уранового сырья, называется далее временем достижения радиационной эквивалентности. Для рассмотренных в статье радионуклидных составов РАО время достижения радиологической эквивалентности меньше времени достижения радиационной эквивалентности.
С помощью ПМ РОЗА-РАО могут проводиться расчёты времени достижения радиологической и радиационной эквивалентности с природным урановым сырьём для произвольных смесей радионуклидов, составляющих РАО.
Материалы и методы
Величина LAR рассчитывалась на основе моделей интенсивности радиационных рисков Международной комиссии по радиологической защите (МКРЗ) [3], а для рака кости - на основе моделей Научного комитета ООН по действию атомной радиации (НКДАР ООН) [4]. Методика вычисления пожизненного атрибутивного риска от внутреннего облучения была изложена ранее [5, 6]. Величина LAR заболеваемости ЗНО, при условии дожития до текущего возраста е, вычисляется путём суммирования избыточного абсолютного риска EAR с весом функции здорового дожития по возрастам дожития a, начиная от текущего возраста е до возраста 100 лет.
LAR(s,c,e,gi.....gn,dí.....dn) = ■ T,a™[sXs,с,е,а)- EAR(s, c,a,gi.....gn,dí.....dn)]. (1)
Здесь DDREF - коэффициент эффективности дозы и мощности дозы; е - возраст на момент поступления радионуклида; S' - функция здорового дожития; amax - максимальный возраст дожития (МКРЗ [3] в качестве amax принимает равным 100 лет); s - пол; с - локализация опухоли; n - количество облучений (зависит от типа радионуклида); g-i,...,gn - возрасты при облучении; d1,...,dn - соответствующие эквивалентные дозы облучения, полученные органом или тканью с; EAR - избыточный абсолютный риск, вычисляемый по моделям МКРЗ [3].
Избыточный абсолютный риск - это приращение фонового показателя заболеваемости ЗНО или смертности от ЗНО после радиационного воздействия. Функция здорового дожития характеризует вероятность для человека заданного пола и возраста е дожить до заданного возраста a и не заболеть при этом ЗНО заданной локализации с.
Вычисление радиационных рисков производится по эквивалентным дозам, полученным из ОЭД. Прямой расчёт по ОЭД недопустим, т.к. эффективная доза не позволяет делать различий между отдельными органами.
Значения эффективных доз преобразуются по специальной методике в значения эквивалентных доз на органы с учётом накопления радионуклидов в организме и длительного их последующего вывода из него. Таким образом, при одном поступлении радионуклида в организм человека происходит его длительное облучение с определённым распределением эквивалентных доз.
Описание методики преобразования единичной ОЭД в динамику эквивалентных доз на органы и ткани было изложено в работе Меняйло и соавт. [5]. В этой работе в качестве базы данных дозовых коэффициентов использовалась распространяемая МКРЗ компьютерная программа «ICRP Database of Dose Coefficients: Workers and Members of the Public» [7]. В ПМ РОЗА-РАО для построения базы данных дозовых коэффициентов используется другой источник: набор утилит под общим названием Dose and Risk Calculation Software (DCAL) [8]. Данный программный комплекс был разработан Американским агентством по защите окружающей среды (United States Environmental Protection Agency, или EPA).
Пакет утилит DCAL [8] предназначен для функционирования в среде операционной системы Windows. С помощью него были рассчитаны дозовые коэффициенты для 721 радионуклида от поступления их в организм человека с пищей.
Для целей определения радиологической эквивалентности используется усреднённый пожизненный атрибутивный риск по всем возможным возрастам населения (0-100 лет).
Пусть численность населения в заданном возрасте е для пола s обозначается как u(s,e), тогда средний LAR в заданном диапазоне возрастов ei и ег будет вычисляться по формуле (2):
LAR(s, с, еъ е2,д1,..., дп, dt,..., dn) =-t^—— • LAR(s, c, e,g±,..., gn, dt,..., dn) • u(s, e). (2)
^e=ei u(s'e)
ПМ РОЗА-РАО вычисляет LAR заболеваемости ЗНО после однократного поступления в организм человека радионуклидного состава РАО в заданном времени выдержки и затем от поступления радионуклидного состава исходного природного уранового сырья. Полученные величины рисков усредняются по возрастам на момент поступления радионуклидов в организм и по полу. Такая процедура для рисков от РАО проводится для каждого времени выдержки. Для природного уранового сырья принимается неизменным во времени состав и дозы облучения.
В результате проведения вычислений ПМ РОЗА-РАО получает дискретные наборы данных вида: точка на оси времени выдержки РАО - суммарная (по радионуклидам) ОЭД и соответствующее значение LAR, а также вычисляется LAR для природного урана от ОЭД 1 мЗв.
Чтобы определить точку пересечения кривых РАО и природного уранового сырья требуется каким-либо способом соединить полученные точки данных. Для этого в ПМ РОЗА-РАО применяются два метода: линейная интерполяция (а-t + b) и экспоненциальная интерполяция (а • bl), где аи b - константы, а t- время выдержки РАО.
Коэффициенты а и b вычисляются по известным точкам данных. Пусть нужно соединить две точки данных с координатами на оси времени выдержки ti и t2. Времени ti соответствует риск ri, а
времени Î2 - риск ri. При этом ri>U>r2, где U - это уровень LAR природного уранового сырья. Тогда можно вывести следующие формулы, определяющие время, при котором значение риска от РАО будет равно значению риска от природного урана, т.е. время наступления радиологической эквивалентности:
для линейной интерполяции:
^S-^-tO + ti, О)
для экспоненциальной интерполяции:
t = ьм-ты - ( - ) +
Для определения времени достижения радиационной (по ОЭД) эквивалентности в формулах (3) и (4) риски следует заменить на соответствующие значения доз.
Из формул (3) и (4) видно, что время, через которое наступает радиологическая эквивалентность при линейной и экспоненциальной интерполяции точек данных, не зависит от абсолютных значений самой шкалы времени. То есть временной шкалой могут быть календарные годы выдержки, время выдержки после начала захоронения и т.д.
В результате алгоритм оценки времени достижения радиологической эквивалентности с использованием двух различных методик интерполяции входных точек данных заключается в следующем:
1. Рассчитать значения LAR для всех входных точек время-доза для РАО.
2. Определить константное значение U уровня риска для природного урана.
3. Определить интервал времени выдержки, в котором левая точка выше значения U, а правая - ниже.
4. По формулам (3) или (4) вычислить время наступления радиологической эквивалентности при линейной или экспоненциальной интерполяции.
Представленные выше методы оценки радиационных рисков и времени достижения радиологической эквивалентности реализует разработанный специализированный программный модуль (ПМ РОЗА-РАО). Он разработан при помощи объектно-ориентированного языка программирования Visual Basic .NET в среде разработки Microsoft Visual Studio Community 2019 и предназначен для работы в среде операционных систем семейства Windows, начиная с Windows 7. Для функционирования ПМ РОЗА-РАО требуется наличие установленной в системе программной платформы .NET Framework 4.7.2.
В качестве входных данных ПМ РОЗА-РАО принимает значения ПБО, выраженные в единицах ОЭД (Зв, мЗв и т.п.) отдельно по каждому радионуклиду РАО и в зависимости от времени выдержки РАО, а также соответствующие постоянные во времени значения ПБО для радионуклидов природного уранового сырья. В ПМ РОЗА-РАО принято, что состав природного уранового сырья постоянен по времени. В расчётах рисков используется ОЭД, полученная путём приведения ПБО РАО на 1 мЗв ПБО исходного природного уранового сырья. Иными словами, для целей расчёта принимается, что суммарная по радионуклидам ОЭД природного урана равна 1 мЗв, а для РАО высчитывается соответствующее значение. Под ОЭД понимается ожидаемая за жизнь эффективная доза от поступления радионуклидов в организм. Приписывается к году поступления этих радионуклидов.
Нормировку входных значений ПБО на 1 мЗв ПБО природного урана необходимо проделать, так как часть моделей годового радиационного риска квадратична по дозе (для ЗНО кости
и лейкозов). При вычислении времени наступления радиологической эквивалентности надо ориентироваться на очень малые дозы, сопоставимые с теми, которые могут быть получены отдельными индивидуумами. Программа позволяет в качестве входных данных задавать показатели фоновой онкозаболеваемости и онкосмертности, показатели общей фоновой смертности, численность населения. Результаты вычислений могут выводиться в экранное окно или в текстовый файл.
Результаты
Интерфейс пользователя ПМ РОЗА-РАО состоит из двух основных окон:
• основное окно (рис. 1), в котором пользователь может выбирать файлы входных данных, видеть результаты вычислений (время наступления радиационной и радиологической эквивалентности), сохранить результаты в выходной текстовый файл;
• окно (рис. 2), на котором представлен график, построенный по результатам вычислений.
fü РОЗА-РАО
I РАО .csv
Пр. урановое cbipbé.csv
Выберите входные фоновые показатели:
Время наступления радиологической (по рискам) и радиационной (по дозам) эквивапентностей, годы
По рискам
По дозам
Линейная Интерполяция: Экспоненциальная интерполяция:
Рис. 1. Главное окно ПМ РОЗА-РАО с результатами вычислений.
Для проведения расчётов времени наступления радиационной и радиологической эквивалентности пользователь сначала должен указать программе входные данные. Чтобы указать входной текстовый файл, содержащий ПБО РАО, нужно нажать на кнопку «Выбрать данные по РАО». Чтобы указать входной текстовый файл, содержащий ПБО исходного природного уранового сырья, нужно нажать на кнопку «Выбрать данные по пр. урану».
Далее пользователь может выбрать заранее подготовленные входные данные, содержащие фоновые показатели заболеваемости ЗНО и смертности от ЗНО, половозрастную структуру и показатели смертности от всех причин для исследуемой популяции. Нажатие на кнопку «Редактировать показатели» откроет окно проводника Windows, где пользователь сможет отредактировать такие входные данные или добавить новые.
После того, как все необходимые входные данные будут выбраны, пользователю станет доступна кнопка «Выполнить расчёты». Нажатие на эту кнопку запустит вычисления рисков LAR для всех времён выдержки РАО и для исходного природного уранового сырья.
По окончании вычислений внизу главного окна отобразятся полученные результаты: время наступления радиологической эквивалентности (по LAR) и время наступления радиационной эквивалентности (по суммарной ОЭД). Также станут доступны для нажатия кнопки «Сохранить в файл» и «Показать на графике».
Нажатие кнопки «Сохранить в файл» потребует от пользователя выбрать в соответствующем окне путь и имя текстового файла, куда ПМ РОЗА-РАО должен экспортировать результаты расчётов.
Нажатие кнопки «Показать на графике» откроет окно ПМ РОЗА-РАО (рис. 2), на котором будет отображён график рассчитанных значений LAR от РАО и исходного природного уранового сырья или значений суммарной (по радионуклидам) приведённой ОЭД РАО.
Ш РОЗА-РАО График - □ X
Логарифмический масштаб по оси X 3 Логарифмический масштаб по оси Y
Средний пожизненный атрибутивный риск (LAR) от однократного поступления с пищей в организм человека изотопов РАО и исходного природного урановго сырья Эффективная доза облучения получена путём приведения ПВО РАО на 1 мЗв ПВО природного уранового сырья
3 0.001
100 200 300 400 500 600 700 800 900 1000
Время выдержки, лет
Линейная интерполяция -Экспоненциальная интерпопяш«
Рис. 2. Окно ПМ РОЗА-РАО, отображающее на графике результаты вычислений.
Переключение между графиками осуществляется с помощью выпадающего списка. Также пользователь может переключать масштабы отображения данных по осям с линейного на логарифмический.
Входные данные. Входные данные для ПМ РОЗА-РАО делятся на три типа:
1. Данные о ПБО РАО в зависимости от времени выдержки и ПБО соответствующего количества исходного природного уранового сырья.
2. Данные о фоновых показателях заболеваемости ЗНО и смертности от ЗНО, половозрастной структуре и показателях смертности от всех причин для исследуемой популяции.
3. Данные о дозовых коэффициентах, которые необходимы, чтобы вычислять эквивалентные дозы на органы и ткани по известным значениям эффективных доз.
Данные первого типа подготавливаются пользователем в виде текстовых файлов с разделителями. В качестве разделителей может служить любой из следующих символов: табуляция, запятая, пробел или точка с запятой. Первая строка всегда содержит заголовки столбцов. Файл для входных данных по урановому сырью содержит два столбца: первый столбец - это радио-нуклидный состав, второй столбец - это соответствующее значение ОЭД. Файл для входных данных по РАО содержит следующие столбцы: первый столбец - это радионуклидный состав,
второй и последующий столбцы - это соответствующие значения ОЭД. Каждый отдельный столбец содержит значения ОЭД для одного значения времени выдержки. Значение времени выдержки должно указываться в первой строке как заголовок столбца. Пример открытых файлов входных данных этого типа представлен на рис. 3.
3 *з ■ □ = Пр, урановое сырьё,csv Excel
Файл Главная Вставка Разметк В Ът □
Д1 И Файл Главная Вставка 33метка страницы Формулы Данные Рецензирование В
А в с N1 4 = г f-
Радионуклид ПБО _1 1_
1
2 Rn-222 1.03 Е+10 А А B C D E 1 F 1 G 1 H i » 1
3 Ро-210 8.23Е+09 1 0 100 200 300 400 600 800 1000
4 РЬ-210 4.73Е+09 2 Ас-227 80400 96400 97000 97000 97000 96900 96800 96800
5 Ra-226 1.92Е+09 3 Ag-108m 52100 0 0 0 0 0 0 0
6 Th-230 1.41Е+09 4 Am-241 1.86E+10 1.87E+10 1.6E+10 1.36E+10 1.16E+10 8.41E+09 6.11E+09 4.43E+09
7 Ac-227 3.43 Е+08 5 Arn-242m 2.22E+09 1.36E+09 8.34E+08 5.1E+08 3.12E+08 1.17E+08 43700000 16300000
8 U-234 3.29Е+08 6 Am-243 1.16E+08 1.15E+08 1.14E+08 1.13E+08 1.12E+08 1.1E+08 1.08E+08 1.06E+08
9 U-238 3.01Е+08 7 Ba-133 90100 0 0 0 0 0 0 0
10 Ра-231 2.21Е+08 8 Bk-249 2930000 0 0 0 0 0 0 0
11 Ra-223 31193473 9 С-14 47000 46500 45900 45400 44800 43700 42700 41700
12 Th-234 22745241 10 Cd-109 73300 0 0 0 0 0 0 0
13 U-235 14621941 11 Cd-113m 9.81E+09 1.56E+08 553000 0 0 0 0 0
14 Bi-210 8935630 12 Се-144 5.62E+11 0 0 0 0 0 0 0
15 Th-227 2707767 13 Cf-249 160000 132000 108000 88800 72900 49100 33100 22300
16 Pb-214 958549.4 14 Cf-250 294000 0 0 0 0 0 0 0
17 Bi-214 752759.2 15 Cm-243 2.87E+08 27100000 2510000 231000 21400 0 0 0
18 Th-231 106144.5 16 Cm-244 4.79E+09 1.08E+08 2350000 51000 0 0 0 0
19 Pb-211 56321.55 17 Cm-245 5480000 5430000 5390000 5350000 5310000 5220000 5140000 5050000
20 18 Cm-246 691000 682000 672000 663000 653000 634000 616000 598000
Пр. урановое сырьё 19 Cs-134 6.56E+08 0 0 0 0 0 0 0
20 Cs-135 20000 20000 20000 20000 20000 20000 20000 20000
■я
i=i=| РАО <±>
I Готово Ш
Рис. 3. Пример открытых текстовых файлов, содержащих входные данные с ОЭД для ПМ РОЗА-РАО. Файлы открыты в программе MS Excel.
Фактическая единица измерения ОЭД (Зв, мЗв и т.д.) в этих файлах не имеет значения, т.к. для расчётов входные данные приводятся к суммарному значению ПБО для природного уранового сырья, равному 1 мЗв.
Данные о фоновых показателях заболеваемости ЗНО и смертности от ЗНО, половозрастной структуре и показателях смертности от всех причин для исследуемой популяции ПМ РОЗА-РАО берёт из входных файлов. Эти данные расположены в папке «Rates», расположенной там же, где и исполняемый файл ПМ РОЗА-РАО. Эта папка открывается в проводнике Windows, если пользователь нажмёт кнопку «Редактировать показатели» в главном окне ПМ РОЗА-РАО. Во время установки ПМ РОЗА-РАО в папку «Rates» помещаются заранее подготовленные данные о фоновых показателях для России в целом и Томской области, которые записываются в папках «Россия в целом» и «Томская область» соответственно. Эти данные получены из справочника [9]. Пользователь может добавлять другие папки в каталог «Rates» по своему усмотрению, но нужно соблюдать структуру файлов в этих папках.
Всего в каждой такой папке, соответствующей одной группировке фоновых показателей, должно содержаться 15 файлов. Файлы имеют расширение *.csv и являются обычными текстовыми файлами, в которых разделителем столбцов служит символ «;». Названия этих файлов следующие: «Все солидные раки.csv», «Желудок^та, «Костиста, «Лёгкие^та, «Лейкозы.csv», «Молочная железа.csv», «Мочевой пузырь.csv», «Остальные солидные раки^та, «Печень.csv», «Пищевод^та, «Толстый кишечник^та, «Щитовидная железа^та, «Яичник.csv», «Числен-ность.csv», «Смертность от всех причин.csv». Названия файлов соответствуют тем данным, которые в них содержатся. Так, файлы с названиями органов и тканей содержат фоновые показа-
тели заболеваемости ЗНО и смертности от ЗНО. Файл «Смертность от всех причин.csv» содержит фоновые показатели общей смертности от всех причин. Данные в этих файлах приведены на 100 000 человек. Файл «Численность.csv» содержит половозрастную численность населения.
При запуске ПМ РОЗА-РАО считывает названия папок из каталога «Rates» и добавляет соответствующие пункты в выпадающий список с заголовком «Выберите входные фоновые показатели:». Папку с описанными выше входными файлами пользователь должен выбрать в этом выпадающем списке перед тем, как нажать кнопку «Выполнить расчёты».
Для преобразования ОЭД в значения эквивалентных доз в органах и тканях ПМ РОЗА-РАО использует динамические ряды дозовых коэффициентов для поступления в организм человека различных радиоактивных веществ с пищей, которые могут быть получены через интерфейс программы EPA DCAL [8].
Выходные данные. Основные результаты работы ПМ РОЗА-РАО выводятся в главном окне программы (рис. 1, внизу) - это время наступления радиологической и радиационной эквивалентности и на графике (рис. 2). Здесь пользователь может переключаться между отображением значений ОЭД и LAR.
Эти же данные в табличном виде, а также время наступления радиационной и радиологической эквивалентности пользователь может сохранить как текстовый файл. Скриншот примера такого файла, открытого в программе «Блокнот», представлен на рис. 4.
Т] Пример выходного файлз.М- Блокнот — D
Файл Правка Формат Вид Справка
Суммарная (по радионуклидам) эффективная доза облучения от РАО, полученая путём приведения ПВО РАО на 1 мЗв ПВО природного уранового сырья
Время выдержки, лет Доза, мЗв
0 66.9511429047846
1Э0 1.80591029692944
200 1.2937667559872
300 0.953891273221637
400 0.786579731327523
600 0.597335642315448
800 0.492308624896472
1000 0.422088620909397
Средний пожизненный атрибутивный риск (LAR) от однократного поступления с пищей в организм человека изотопов РАО, расчитанный по дозам выше
Время выдержки, лет LAR
0 0. ,00278482376837883
100 2. ,84293493633016Е-05
200 1. . 93793641896731Е-05
300 1. ,30580408297115Е-05
4Q0 1. ,06432525631059Е-05
600 8. ,01999553847326Е-06
800 6. ,61519775348103Е-06
1000 5. ,69053719275005Е-06
от однократного поступления с пищей в организм человека изотопов природного уранового сырья в эффективной дозе 1 мЗв : 6.08598336322731Е-05
Время наступления радационной и радиологической эквивалентностей
Радиологическая эквивалентность (по ЬАЯ), лет Радиационная эквивалентность (по дозам), лет
Линейная интерполяция: 98.82 286.43
Экспоненциальная интерполяция: 83.40 284.51
Рис. 4. Пример выходного файла, открытого в программе «Блокнот».
Пример проведения расчёта при помощи РОЗА-РАО. В опубликованной ранее статье [10] решалась задача определения времени, через которое наступит равенство ОЭД РАО и ОЭД исходного природного сырья для смоделированного сценария развития в 21 веке ядерной энергетики России на основе АЭС с реакторами на тепловых и быстрых нейтронах. Исходными данными в этом сценарии для определения радиационной и радиологической эквивалентности являлись ОЭД от основных изотопов РАО и изотопов природного урана в различные периоды времени. Эти данные получены при моделировании сценария развития ядерной энергетики на тепловых реакторах и быстрых реакторах до 2100 г. и представлены в таблицах в той же работе [10].
В данной работе исходные данные из статьи [10] были подготовлены для проведения расчётов при помощи ПМ РОЗА-РАО (рис. 4). В результате вычислений было получено (рис. 1 и 4),
что время наступления радиологической (по LAR) эквивалентности равно 98,82 года по линейной интерполяции и 83,4 года по экспоненциальной интерполяции, время наступления радиационной (по ОЭД) эквивалентности равно 286,43 года по линейной интерполяции и 284,51 года по экспоненциальной интерполяции.
В самой статье [10] было показано, что равенство ОЭД для населения от РАО и природного урана (радиационная эквивалентность) достигается через 287 лет выдержки РАО, а равенство радиационных рисков (радиологическая эквивалентность) - через 99 лет выдержки. Таким образом, ПМ РОЗА-РАО проводит адекватные вычисления, поэтому результаты работы данного ПМ могут использоваться при обоснованиях методов обращения с РАО в проектном направлении «Прорыв».
Обсуждение
ПМ РОЗА-РАО рассчитывает время наступления радиологической эквивалентности на основе метрики пожизненного риска LAR для заболеваемости населения всеми локализациями ЗНО в совокупности, после однократного поступления радионуклидов в организм с пищей и водой. При этом LAR рассчитывается на основе моделей МКРЗ [3]. Результаты оценки времени наступления радиологической эквивалентности с использованием других метрик риска, вариантов облучения, моделей радиационного риска, источника дозовых коэффициентов и т.д. могут отличаться от той, что даёт ПМ РОЗА-РАО. Проведённые дополнительные расчёты, выполненные перед разработкой ПМ РОЗА-РАО, показали:
• Оценка времени наступления радиологической эквивалентности с использованием пожизненных метрик радиационного риска даёт в основном приблизительно одинаковые результаты независимо от пола, возраста и числа облучений, типа риска (онкозаболеваемость или он-космертность).
• Оценка времени наступления радиологической эквивалентности с использованием пожизненных радиационных рисков, посчитанных только для заданного органа или ткани, может существенно отличаться от такой оценки, полученной по рискам, посчитанным по всем органам в целом. Отличие может быть как в меньшую сторону (до 1,8 раза), так и в большую (до 7,7 раза).
• Сравнение оценок времени достижения радиологической эквивалентности, основанных на радиационных рисках, рассчитанных с использованием моделей радиационных рисков МКРЗ [3] и моделей радиационных рисков НКДАР ООН [5], показывает, что последние могут быть значительно (до 2,18 раза для метрик пожизненного риска и 2,4 раза для метрик годового риска) завышены.
• Использование разных источников для дозовых коэффициентов [7, 8] не оказывает существенного влияния на оценки времени достижения радиологической эквивалентности.
• Линейная интерполяция точек данных позволяет получить значительно более устойчивые оценки времени наступления радиологической эквивалентности. Но в то же время, эти оценки всегда выше, чем аналогичные, но посчитанные с использование экспоненциальной интерполяции.
Модели радиационного риска МКРЗ [3] идентифицировали по данным когорты выживших после атомной бомбардировки Хиросимы и Нагасаки в 1945 г. Для оценки применимости этих моделей для расчёта LAR для жителей Российской Федерации требуются дополнительные исследования.
Модели риска МКРЗ [3] не дают точного результата. Параметры этих моделей имеют неопределённости [11, 12]. Также на неопределённость результатов расчёта влияет неопределённость доз облучения, которая зависит, в частности, от неопределённости дозовых коэффициентов. Ещё одним фактором неопределённости в оценке LAR является неопределённость фоновых значений онкозаболеваемости и онкосмертности. ПМ РОЗА-РАО в настоящее время никак не учитывает эти неопределённости.
В ПМ РОЗА-РАО для расчётов применяются фоновые показатели за фиксированный календарный год. Однако при оценке радиологической эквивалентности, которая может наступить более чем через 100 лет, такой подход вносит погрешность, т.к. фоновые показатели не являются статичными, а меняются во времени.
Заключение
В целях разработки программного обеспечения для обоснования методов обращения с РАО в замкнутом ядерном топливном цикле на основе быстрых реакторов разработан ПМ РОЗА-РАО, предназначенный для оценки времени наступления радиационной и радиологической эквивалентности РАО и исходного природного уранового сырья.
Впервые в атомной отрасли с помощью разработанного ПМ РОЗА-РАО оценка радиологической состоятельности ОДЭК при переработке облучённого ядерного топлива и обращении с РАО может быть получена на основе национальных требований к радиационной защите населения с учётом выполнения современных международных рекомендаций и стандартов.
ПМ РОЗА-РАО позволяет определять параметры переработки облучённого ядерного топлива (состав и активности радионуклидов для контролируемого хранения, сжигания и трансмутации в быстрых реакторах) и обращения с РАО на основе минимизации потенциальных радиационных рисков населения, что соответствует современной правоприменительной практике развитых стран в области радиационной защиты.
Публикация подготовлена по результатам выполнения работ в Госкорпорации «Росатом» в рамках ПН «Прорыв».
Литература
1. Адамов Е.О., Ганев И.Х. Экологически безупречная ядерная энергетика. М.: НИКИЭТ им. Н.А. Доллежаля, 2007. 145 с.
2. Лопаткин А.В. Радиационно-эквивалентное обращение с РАО. Техническая справка 01.2017 НРРЭ. М., 2017. 21 с.
3. ICRP, 2007. The 2007 Recommendations of the International Commission on Radiological Protection. ICRP Publication 103 //Ann. ICRP. 2007. V. 37, N 2-4. P. 1-332.
4. United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR). Sources and effects of ionizing radiation. UNSCEAR 2006 Report Vol. I, Annex A: Epidemiological studies of radiation and cancer. New York: United Nation, 2008.
5. Меняйло А.Н., Чекин С.Ю., Кащеев В.В., Максютов М.А., Корело А.М., Туманов К.А., Пряхин Е.А., Ловачев С.С., Карпенко С.В., Кащеева П.В., Иванов В.К. Пожизненный радиационный риск в результате внешнего и внутреннего облучения: метод оценки //Радиация и риск. 2018. Т. 27, № 1. С. 8-21.
6. Меняйло А.Н., Чекин С.Ю., Ловачёв С.С., Туманов К.А., Корело А.М., Иванов В.К. Программный модуль радиологической защиты населения //Радиация и риск. 2020. Т. 29, № 2. С. 32-48.
7. ICRP Database of Dose Coefficients: Workers and Members of the Public; Ver. 3.0, official website. [Электронный ресурс]. URL: http://www.icrp.org/page.asp?id=402 (дата обращения 08.04.2021).
8. DCAL Software and Resources. [Электронный ресурс]. URL: https://www.epa.gov/radiation/dcal-software-and-resources (дата обращения 13.04.2021).
9. Злокачественные новообразования в России в 2017 году (заболеваемость и смертность) /под ред. А.Д. Каприна, В.В. Старинского, Г.В. Петровой. М.: МНИОИ им. П.А. Герцена, 2018.
10. Иванов В.К., Чекин С.Ю., Меняйло А.Н., Максютов М.А., Туманов К.А., Кащеева П.В., Ловачёв С.С., Адамов Е.О., Лопаткин А.В. Радиационная и радиологическая эквивалентность РАО при двухкомпо-нентной ядерной энергетике //Радиация и риск. 2019. Т. 28, № 1. С. 5-25.
11. Preston D.L., Kusumi S., Tomonaga M., Izumi S., Ron E., Kuramoto A., Kamada N., Dohy H., Matsui T., Nonaka H., Thompson D.E., Soda M., Mabuchi K. Cancer incidence in atomic bomb survivors. Part III: Leukemia, lymphoma and multiple myeloma, 1950-1987 //Radiat. Res. 1994. V. 137 (Suppl.). P. 68-97.
12. Preston D.L., Ron E., Tokuoka S., Funamoto S., Nishi N., Soda M., Mabuchi K., Kodama K. Solid cancer incidence in atomic bomb survivors: 1958-1998 //Radiat. Res. 2007. V. 168, N 1. P. 1-64.
Software module for assessing the achievement of radiological equivalence (ROZA-RAO)
Menyajlo A.N., Chekin S.Yu., Lovachev S.S., Tumanov K.A., Ivanov V.K
A. Tsyb MRRC, Obninsk
The main problem of the large-scale development of nuclear energy in the world is the safe disposal of accumulated radioactive waste (RW). The hazard of RW at the time of their formation to human health far exceeds the hazard of the source natural uranium ore, after the use of which these wastes were formed. Due to the natural processes of radioactive decay, the hazard of radioactive waste decreases over time. The hazard of waste can also be reduced by artificial transmutation or incineration of the radionuclides that make up RW composition. For uranium nuclear fuel cycles, the natural reference level of RW hazard is the hazard level of the uranium ore, upon reaching which decisions can be made on the final disposal of RW. The article presents a method for estimating the time of achieving radiological equivalence of the mass of accumulated RW of nuclear energy and the source mass of natural uranium raw materials used in this case. A description of the software module "Radiological support of protection - radioactive waste" (PM ROZA-RAO), which implements this method, is also presented. Radiological equivalence herein refers to the equality of the corresponding radiation risks for the population. The methodology used is determined by the current recommendations of the International Commission on Radiological Protection (ICRP) and the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) on the calculation of human radiation risks from internal exposure based on equivalent doses in organs and tissues. The determination of equivalent doses to organs and tissues was based on dose coefficients of radioactive substances provided for widespread use by the United States Environmental Protection Agency (EPA). Calculations by means of the PM ROZA-RAO showed that within the framework of the most likely scenario for the development of two-component nuclear energy system in Russia, for radioactive waste accumulated by 2100, the time of onset of radiological equivalence is less than 100 years (98.8 years). With the help of the developed software module, for the first time in the nuclear industry, the parameters of the processing of exposed nuclear fuel (the composition and activity of radionuclides for controlled storage, combustion and transmutation) and the management of RW can be determined, minimizing the potential radiation risks of the population, which corresponds to the best practice of developed countries in the field of radiation protection.
Key words: radioactive waste, natural uranium raw materials, radiation equivalence, radiological equivalence, radiation safety, radiation risk, radiobiological effects, malignant neoplasm, radiation risk models of ICRP and UNSCEAR, dose coefficients, internal exposure, radionuclides, expected effective dose, equivalent dose.
References
1. Adamov E.O., Ganev I.Kh. Ecologically impeccable nuclear energy. Moscow, N.A. Dollezhal NIKIET, 2007. 145 p. (In Russian).
2. Lopatkin A.V. Radiation-equivalent treatment of radioactive waste. Technical Reference 01.2017 NPRE. Moscow, 2017. 21 p. (In Russian).
3. ICRP, 2007. The 2007 Recommendations of the International Commission on Radiological Protection. ICRP Publication 103. Ann. ICRP, 2007, vol. 37, no. 2-4, pp. 1-332.
4. United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR). Sources and effects of ionizing radiation. UNSCEAR 2006 Report Vol. I, Annex A: Epidemiological studies of radiation and cancer. New York, United Nation, 2008.
5. Menyajlo A.N., Chekin S.Yu., Kashcheev V.V., Maksioutov Korelo A.M., Tumanov K.A., Pryakhin E.A., Lovachev S.S., Karpenko S.V., Kashcheeva P.V., Ivanov V.K. Lifetime attributable risks from external
Menyajlo A.N.* - Senior Researcher, C. Sc., Biol.; Chekin S.Yu. - Head of Lab.; Lovachev S.S. - Research Assistant; Tumanov K^. - Head of Lab., C. Sc., Biol.; Ivanov V.K. - Scientific Advisor of NRER, Chairman of RSCRP, Corr. Member of RAS, D. Sc., Tech. A. Tsyb MRRC. •Contacts: 4 Korolyov str., Obninsk, Kaluga region, Russia, 249035. Tel.: (484) 399-32-81; e-mail: [email protected].
and internal exposure to radiation: method for estimating. Radiatsiya i risk - Radiation and Risk, 2018, vol. 27, no. 1, pp. 8-21. (In Russian).
6. Menyajlo A.N., Chekin S.Yu., Lovachev S.S., Tumanov K.A., Korelo A.M., Ivanov V.K. Software module for radiological protection of the public. Radiatsiya i risk - Radiation and Risk, 2020, vol. 29, no. 2, pp. 32-48. (In Russian).
7. ICRP Database of Dose Coefficients: Workers and Members of the Public; Ver. 3.0, official website. Available at: http://www.icrp.org/page.asp?id=402 (Accessed 08.04.2021).
8. DCAL Software and Resources. Available at: https://www.epa.gov/radiation/dcal-software-and-resources (Accessed 13.04.2021).
9. Malignant neoplasms in Russia in 2017 (morbidity and mortality). Eds.: A.D. Kaprin, V.V. Starinskiy, G.V. Petrova. Moscow, P. Hertsen MORI, 2018. (In Russian).
10. Ivanov V.K., Chekin S.Yu., Menyajlo A.N., Maksioutov IWA., Tumanov K.A., Kashcheeva P.V., Lovachev S.S., Adamov E.O., Lopatkin A.V. Radiation and radiological equalities between natural uranium and radioactive waste in innovative two-component nuclear energy system. Radiatsiya i risk - Radiation and Risk, 2019, vol. 28, no. 1, pp. 5-25. (In Russian).
11. Preston D.L., Kusumi S., Tomonaga M., Izumi S., Ron E., Kuramoto A., Kamada N., Dohy H., Matsui T., Nonaka H., Thompson D.E., Soda M., Mabuchi K. Cancer incidence in atomic bomb survivors. Part III: Leukemia, lymphoma and multiple myeloma, 1950-1987. Radiat. Res., 1994, vol. 137 (Suppl.), pp. 68-97.
12. Preston D.L., Ron E., Tokuoka S., Funamoto S., Nishi N., Soda M., Mabuchi K., Kodama K. Solid cancer incidence in atomic bomb survivors: 1958-1998. Radiat. Res., 2007, vol. 168, no. 1, pp. 1-64.
