CC BY
29
DOI: 10.21870/0131-3878-2020-29-2-32-48 УДК 621.039:614.876:504.055
Программный модуль радиологической защиты населения
Меняйло А.Н., Чекин С.Ю., Ловачёв С.С., Туманов К.А., Корело А.М., Иванов В.К.
МРНЦ им. А.Ф. Цыба - филиал ФГБУ «НМИЦ радиологии» Минздрава России, Обнинск
В статье представлено описание разработанного программного модуля радиологической защиты населения (ПМ РОЗА-Н). Данное программное обеспечение предназначено для расчёта радиационных рисков населения от возможного воздействия ядерных объектов, находящихся в районе проживания. Применяемая методология определяется современными рекомендациями Международной комиссии по радиологической защите (МКРЗ) и Научного комитета оОн по действию атомной радиации (НКДАР ООН). Определение доз внутреннего облучения основывается на дозовых коэффициентах, полученных по данным Американского агентства по защите окружающей среды (United States Environmental Protection Agency (EPA)). В работе рассмотрен пример расчёта с помощью ПМ РОЗА-Н пожизненного атрибутивного риска (LAR) онкосмертности для населения, проживающего вблизи предприятий опытного демонстрационного энергокомплекса российской ядерной энергетической системы нового поколения (проектное направление «Прорыв»). Расчёт демонстрирует существенное (десятикратное) снижение радиационного риска с увеличением возраста при облучении: от 2,1 х 10-7 для детей младше 10 лет до 1,96х10"8 для лиц возраста 70-79 лет. Усреднённая по возрасту и полу величина риска равна 1,00х10-7, что в 2 раза меньше оценки, получаемой с использованием номинального коэффициента риска НРБ-99/2009 (1,99х10-7). Новизна полученных результатов заключается в том, что впервые, благодаря разработанному программному модулю, радиационная безопасность населения, проживающего вблизи территории размещения ядерных объектов, может быть обоснована с учётом современных международных стандартов. Экономическая эффективность и значимость ПМ РОЗА-Н определяется тем, что он позволяет обосновать радиационную безопасность при выбросах любого радионук-лидного состава, обеспечивая оперативное планирование мероприятий по защите населения. При этом оценки риска могут производиться не только для населения как популяции, но и любого индивидуума, проживающего вблизи территории размещения ядерных объектов, что соответствует современной правоприменительной практике развитых стран.
Ключевые слова: радиационная безопасность, радиационный риск, радиобиологические эффекты, радиоиндуцированные онкологические заболевания, модели радиационного риска МКРЗ и НКДАР ООН, дозовые коэффициенты, путь поступления радионуклидов в организм, внешнее облучение, внутреннее облучение, радионуклиды, AMAD, ожидаемая эффективная доза, эквивалентная доза.
Введение
Программный модуль радиологической защиты населения (ПМ РОЗА-Н) предназначен для оценки возможного воздействия на население в терминах риска для здоровья человека от техногенного облучения. Для этого ПМ РОЗА-Н обеспечивает точечную оценку индивидуальных пожизненных радиационных рисков населения, проживающего вблизи ядерных объектов в следующих метриках: пожизненный атрибутивный риск онкозаболевания или онкосмерти LAR (от англ. lifetime attributable risk), риск радиационно-индуцированного случая REIC (от англ. risk of exposure induced cancer) или риск радиационно-индуцированной смерти REID (от англ. risk of exposure induced death), пожизненная атрибутивная доля LARF (от англ. lifetime attributable risk fraction), потерянные годы жизни или потерянные годы здоровой жизни YLL (от англ. years of life lost), потерянные годы жизни, приведённые на радиационно-индуцированный случай онкозаболевания или онкосмерти YLLRIC (от англ. years of life lost per radiation-induced cancer).
Меняйло А.Н.* - вед. научн. сотр., к.б.н.; Чекин С.Ю. - зав. лаб.; Ловачёв С.С. - мл. научн. сотр.; Туманов К.А. - зав. лаб., к.б.н.; Корело А.М. - ст. научн. сотр.; Иванов В.К. - зам. директора по научн. работе, Председатель РНКРЗ, чл.-корр. РАН. МРНЦ им. А.Ф. Цыба - филиал ФГБУ «НМИЦ радиологии» Минздрава России.
•Контакты: 249035, Калужская обл., Обнинск, ул. Королёва, 4. Тел.: (484) 399-32-81; e-mail: nrer@obninsk.com.
Программа может вычислять указанные радиационные риски для населения после внешнего облучения и, одновременно, после внутреннего облучения путём ингаляционного или пе-рорального поступления различных радионуклидов. Всего в ПМ РОЗА-Н можно провести расчёт рисков от внутреннего облучения для 746 различных изотопов, включая 25 благородных газов. Внешнее облучение от изотопов не учитывается. Пожизненные риски вычисляются для всех возможных возрастов (от 0 до 100 лет) в заданном в программе календарном году.
Риски рассчитываются на основе моделей годовых радиационных рисков Международной комиссии по радиологической защите (МКРЗ) [1] и частично (рак кости) Научного комитета ООН по действию атомной радиации (НКДАР ООН) [2]. Программа позволяет в качестве входных данных задавать показатели фоновой онкозаболеваемости и онкосмертности, показатели общей фоновой смертности, численность населения.
Программа предназначена для работы в среде операционных систем семейства Windows, начиная с Windows 7. Графический интерфейс пользователя ПМ РОЗА-Н представляет собой стандартный оконный графический интерфейс Windows.
ПМ РОЗА-Н был разработан при помощи объектно-ориентированного языка программирования Visual Basic .NET. Для функционирования ПМ РОЗА-Н требуется наличие установленной в системе программной платформы .NET Framework 4.7.
Новизна разработки заключается в том, что данный программный модуль впервые в атомной отрасли позволит быстро (за несколько минут) получать оценки радиационных рисков населения от радиоактивных выбросов любого радионуклидного состава, обеспечивая повышение противоаварийной готовности, оперативное планирование мероприятий по защите населения, а также возможность разработки уточнённых программ радиационного контроля, учитывающих оценки потенциальных радиационных рисков населения.
Настоящая работа выполнена в рамках проектного направления «Прорыв» Госкорпорации «Росатом».
Материалы и методы
Расчёт пожизненных радиационных рисков в ПМ РОЗА-Н основан на применении моделей риска, предложенных МКРЗ в Публикации 103 [1] и частично (рак кости) НКДАР ООН [2]. Публикация 103 МКРЗ [1] и отчёт НКДАР ООН [2] содержат модели интенсивностей радиационных рисков для различных органов и тканей человека, которые применяются для обоснования современной системы радиационной защиты.
Пусть X0 - показатель фоновой заболеваемости злокачественными новообразованиями (ЗНО) или смертности от них (отношение числа случаев за один год к полному числу лиц, находящихся под наблюдением за этот год). Воздействие радиации приводит к увеличению X0 на некоторую дополнительную величину SX. Согласно МКРЗ [1] радиогенный риск SX математически может быть представлен в виде двух форм: аддитивной и мультипликативной. В аддитивной форме EAR не зависит от фоновой заболеваемости, т.е. SX=EAR, где EAR - это избыточный абсолютный риск (от англ. excess absolute risk). В мультипликативной модели SX=X0ERR, где безразмерная величина ERR есть избыточный относительный риск (от англ. excess relative risk).
В работе [3] подробно даны формулы, по которым следует проводить расчёт SX согласно МКРЗ, но там не приводится модель расчёта EAR и ERR для рака кости, приведём её ниже.
Для рака кости использовалась модель НКДАР ООН [2], которая была приведена в соответствующем отчёте этой организации в табл. 52. Там даны как мультипликативная, так и аддитивная модели радиационно обусловленной заболеваемости раком кости. Формула (1) описывает мультипликативную модель, а формула (2) - аддитивную.
ERR НКДАР ( a , d ) = d 2 . 6,90379 . 10 7 . a ~4'472 , (1)
однокр. \ ' / ' ' v/
EAR НКДАР ( d ) = d 2 . 9,32940 . 10 -6 . (2)
однокр. V / ' v '
Здесь d- доза облучения в Гр; a - возраст, на который рассчитывается риск.
Два подхода в определении величины SX (аддитивный и мультипликативный) дают разный численный результат для популяций, отличных от тех, на которых строились модели радиационного риска. В этом случае МКРЗ [1] рекомендует вычислять 8Х как взвешенную сумму аддитивного и мультипликативного подходов со следующими весами мультипликативной модели: 0 - для молочной железы и костного мозга (лейкозы), 1 - для щитовидной железы, 0,3 - для лёгкого и 0,5 - для всех остальных локализаций.
Далее под EAR будем понимать величину 8Х, полученную с учётом такого взвешивания. НКДАР ООН 2] не даёт указаний как усреднять мультипликативную и аддитивную модель, но в ПМ РОЗА-Н мы применяем вес каждой модели, равный 0,5. Кроме того, в указанных выше моделях риска учитывается так называемый латентный период, т.е. минимальное время, через которое могут реализоваться радиационные риски. Для лейкозов латентный период принимается равным 2 года, а для солидных раков - 10 лет. При многократном облучении итоговый EAR в заданном возрасте будет суммой EAR, рассчитанных на этот возраст от каждого отдельного облучения.
Для вычисления пожизненных мер радиационного риска используются так называемые функция дожития (для вычисления рисков смертности от ЗНО) и функция здорового дожития (для вычисления рисков возникновения ЗНО).
Функция дожития S(s,e,a) характеризует вероятность для человека пола s и возраста e дожить до возраста а. Формула для вычисления S представлена ниже.
S (s,e,a) = exp
a > e
При этом в возрасте e и ранее функция дожития по определению принимается равной единице, то есть S (s ,e,a ) = 1 при a < e .
Здесь Л^бщ'смерт' (s,k ) - показатель общей фоновой смертности в год от всех причин для лиц пола s и в возрасте a.
Функция здорового дожития S'(s,c,e,a) характеризует вероятность для человека пола s и возраста e дожить до возраста а без возникновения ЗНО локализации с:
Г a"1 1
л / / \ 1 я общ . смерт./. \ . заб. t . ч „смерт./ . ч 1
S (s,c,e,a) = exp р (s,k ) + 20 (s'k'c)-^„ р (s,k,c )j (4)
a > e
a -1
z
k = e
общ . смерт .
( s' k )
(3)
Здесь лСмерт' - показатель смертности от ЗНО в год для заданного пола s, возраста a и
локализации с; Л™6' - показатель заболеваемости ЗНО в год для заданного пола s, возраста a и локализации с.
Функция дожития Sdos е( s ,c ,e,a ,g 1, ...,g n,d 1, ...,d n) учитывает увеличение смертности за счёт радиационного облучения. Она вычисляется аналогично функции дожития (3), но к показателю общей смертности добавляется EAR смертности от ЗНО:
Г a1 Я °6Щ ■ с""еРт ■ / g \ + ]
S dose ( s 'c'e'a '91.....3n'd1.....dn )= e*P 1-Е померк,' ). . .J. (5)
L k + EAR ( s,c'k'9i.....9n'di.....d n )J
a > e
Функция здорового дожития S ¿ose (s ,c,e,a ,g1 ,...,g n,d 1 ,...,d n) учитывает увеличение уровня заболеваемости ЗНО в зависимости от дозы облучения.
О. / - ч Г ^ Л°6щ-с"""""-(s,k ) + л"6-(S'k'C)-лсмерт(S'k'C) + !
SLse (S'C'e'a'g1'...'9n'd1'...'dn )= expl-£ 0 d3¡6. / ' ' 0 ( ' ) ° ( ' ' ) I . (6)
L к + EAR (s,c,k' g1'.. ' gn,di'...,dn) J
a > e
В (5) и (6) s - пол; с - локализация опухоли; n - количество облучений; gi,...,gn - возрасты на момент облучения; db...,dn - соответствующие эквивалентные дозы облучения, полученные органом или тканью с; е - текущий возраст человека, в котором он жив; EAR336' и EARCMepm' - избыточные абсолютные риски (с учётом взвешивания мультипликативной и аддитивной моделей) возникновения ЗНО и смерти от ЗНО соответственно, для заданного пола s, локализации опухоли с, на заданный возраст a, после многократного облучения. EARCMepm' получается из EAR336' путём умножения на так называемую долю летальности. Такой подход в расчёте рисков смерти - это рекомендация МКРЗ.
Пожизненный атрибутивный риск (LAR) характеризует количество избыточных над фоновым числом раковых заболеваний или смертей, которые могут произойти в течение всей последующий жизни после облучения группы лиц, у которых одинаковы модифицирующие риск факторы (пол, возраст при облучении и др.). LAR вычисляется путём суммирования EAR с весом вероятности здорового дожития (или просто дожития при расчёте рисков смерти) по всем возможным возрастам дожития, начиная от текущего возраста, в котором исследуемое лицо или группа лиц считаются живыми.
Учитывая зависимость EAR от пола, исследуемой локализации опухоли, возраста дожития, возраста при облучении и доз облучения представим общие формулы, по которым проводится расчёт LAR в ПМ РОЗА-Н. Формулы (7) и (8) отражает вычисления пожизненного атрибутивного риска возникновения ЗНО и смерти от ЗНО соответственно.
a
1 max í~!\
LAR 336 ( s 'C, e, gi,..., g n, d i,..., d n ) = -■ £ [S'(s,c,e,a )■ EAR за6' (s ,c, a, g ^ ..., g n,d ^ ...,d n )] , (7)
DDREF !|re L J
a
1 W aX /Q\
(s' C'e' g f'...' g n'd f'...' d n )= ■ £ [S (s^'a )■ EAR смерт ' (s' C'a' g,'...'g n'd,'...'d n)] . (8)
LAR
В (7) и (8) DDREF - коэффициент эффективности дозы и мощности дозы; S' - вероятность здорового дожития (без возникновения онкозаболевания локализации с) от возраста е до возраста a; S - вероятность дожития от возраста е до возраста a; amax - максимальный возраст дожития. МКРЗ в качестве amax принимает 100 лет.
Согласно МКРЗ [1] DDREF - это «Экспертно оцененный параметр, который объясняет обычно сниженную биологическую эффективность (на единицу дозы) радиационного воздействия малых доз и малых мощностей доз, если сравнивать её с эффективностью высоких доз и высоких мощностей доз». МКРЗ принимает величину DDREF, равной 2.
Риск радиационно-индуцированного случая (REIC) и риск радиационно-индуцированной смерти (REID) были введены НКДАР ООН [2]. REID согласно НКДАР ООН определяется как разность между уровнями смертности по данной причине в облучённой и необлучённой популяциях для заданного пола и возраста на момент облучения, становящегося дополнительной причиной смерти в облучённой популяции.
Риск REID полностью аналогичен риску LAR для смертности от рака. Различие заключается в том, что при вычислении LAR используется вероятность дожития S (3), а для вычисления REID - особая функция дожития Sdose (5), где учитывается уменьшение вероятности дожития за счёт смерти от ЗНО, вызванных облучением. Кроме того, НКДАР ООН [2] проводит нормировку REID на единицу дозы, но в ПМ РОЗА-Н, в целях сравнения с другими величинами пожизненного риска, эта нормировка не выполняется:
Sdose (s,c,e,a,g1,...,gn,d1,...,d п) ■ EARсмерт' (s,c,a,g1,...,gn,d1,...,dn)
REID (s ,c,e, g i,g п,d i, ...,d п) = £
a = e
Величина REIC (риск заболеваемости ЗНО) вычисляется аналогично REID (9):
a
m a
REIC (s ,c,e, g i,..., g n,d i, ...,d n )= £
Sdose ( S'C'e'a'gi'...'gn'd1'..'dn ) ■
■EAR заб' (s ,c, a ,g 1,..., g n ,d 1, ...,d n)
(9)
(10)
Величины пожизненного риска YLL и YLLRIC были введены в отчёте НКДАР ООН [2]. Они описывают радиационно-обусловленное уменьшение продолжительности жизни. YLL - потерянные годы жизни; YLLRIC - потерянные годы жизни, приведённые на радиационно-индуцированный случай онкосмерти.
Формула для вычисления YLL без нормировки на дозу согласно НКДАР ООН [2] выглядит следующим образом:
a
m ax
YLL (s ,c ,e,g i, ...,g n,d i, ■■■,d n )= £ [S (s,e,a)-^ (s ,c,e,a ,g,, ...,g n,d i, ...,d n)] . (11)
a = e
YLLRIC вычисляется по следующей формуле:
YLL (s,c,e,g.,...,g ,d,,...,d ) YLLRIC (s ,c,e,g ,...,g ,d1,...,d ) = -( ' ' ^ n i n) (12)
V ' ' '»i' n ' i' ' n> ni-ifW -i -i \
REID (s ,c,e,g i, ..^g n,d i, ...^ n) Величины, аналогичные YLL и YLLRIC, но только для оценки потерь лет здоровой жизни также вычисляются в ПМ РОЗА-Н.
a
max
YLL,nc ( s,c,e,gi.....g n,di.....dn )= £ [S '(s.c.e.a)-Sdose ^.c^^.g,.....gnd,.....dn)] . (13)
a = e
YLL (s.c.e.a .....a ,d .....d )
YLLRIC (s,c,e,g1.....g ,d1.....d )= -inc ( , , ^ n' i , n) . (14)
'nc n n REIC ^^^^i.....gn,di.....dn)
Для усреднения аддитивных величин пожизненного риска (LAR, LBR, YLL, REIC, REID) с весом численности населения требуется знать диапазон текущих возрастов, в которых произ-
a
max
водится усреднение и численность населения исследуемого региона в каждом текущем возрасте. Пусть, например, требуется вычислить средний LAR в диапазоне текущих возрастов e1 и e2. Пусть численность населения в заданном возрасте e для пола s обозначается как u(s, e). Тогда средний LAR в заданном диапазоне возрастов будет вычисляться по формуле (15):
LAR ( s,c,e1,e2,g1,...,gn,d1,...,dn )- —-1-• £ LAR ( s.c.e ,g „....g „d ^....d я )• и ( s,e ) . (15)
2: и (s.e ) e-e-
e -e,
Для производных величин риска (LARF, YLLRIC) среднее вычисляется по средним аддитивным величинам (LAR и LBR в случае LARF и YLL и REID (REIC) в случае YLLRIC).
ПМ РОЗА-Н позволяет вычислять пожизненные радиационные риски от внутреннего облучения после попадания заданных радионуклидов в организм человека. В этом случае вычисляются эквивалентные дозы облучения по заданным в программу ожидаемым эффективным дозам (ОЭД) или объёмной активности исследуемого радионуклида в воде или в воздухе. Методика определения эквивалентных доз по базе данных дозовых коэффициентов (Зв/Бк) была нами описана ранее [3]. В этой работе в качестве источника для базы данных дозовых коэффициентов использовалась распространяемая МКРЗ компьютерная программа «ICRP Database of Dose Coefficients: Workers and Members of the Public» [4]. В ПМ РОЗА-Н для построения базы данных дозовых коэффициентов используется другой источник: набор утилит под общим названием Dose and Risk Calculation Software (DCAL) [5]. Данный программный комплекс был разработан Американским агентством по защите окружающей среды (United States Environmental Protection Agency (EPA)).
Пакет утилит DCAL [5] предназначен для функционирования в среде операционной системы Windows. С помощью DCAL [5] были рассчитаны дозовые коэффициенты для 721 радионуклида. Для каждого радионуклида проводились расчёты для поступления этого радионуклида в организм человека с пищей, ингаляционному поступлению для десяти различных значений AMAD, мкм: (0,001; 0,003; 0,01; 0,03; 0,1; 0,3; 1,0; 3,0; 5,0; 10,0) и трёх вариантов растворимости: медленная (S), средняя (M) и быстрая (F), а также для ингаляционного поступления химических соединений для некоторых радионуклидов (табл. 3). Всего было проведено 22423 отдельных расчёта.
Дозовые коэффициенты для инертных газов невозможно получить в пакете DCAL [5]. Для этих радионуклидов применяется отдельный способ определения дозовых коэффициентов. В табл. C1 Публикации 119 МКРЗ [6] приведены мощности эффективных доз на единицу объёмной активности (Зв/сут на Бк/м ) для взрослых для ингаляционного поступления в организм инертных газов. Чтобы получить из этой таблицы дозовые коэффициенты (Зв/Бк) требуется умножить данные из таблицы на 365 дней и разделить объём вдыхаемого воздуха взрослым человеком в год (8100 м /год для взрослого человека по данным [7]).
Результаты
Разработанный ПМ РОЗА-Н состоит из двух основных окон: окно для ввода данных (рис. 1) и окно для вывода результатов вычислений (рис. 2). Пользователь может выбрать язык интерфейса программы (русский или английский), нажав на соответствующую кнопку (рис. 1). Для расчёта пожизненных рисков пользователь должен ввести все необходимые для проведения вычислений данные в окне ввода данных и нажать кнопку «Вычислить риски».
Выберите входные данные:
Тип расчета: -Пол:- -Добавить оолучение:—
# Заболеваемость # Мужской
О Смертность О Женский
РОЗА-Н
РАДИОЛОГИЧЕСКОЕ О&ЕСПЕЧЕНИЕ И ЗАЩИТА
Диапазон календарных лет Эффективная доза внешнего облучения. мЗв С-14: Пероральный; Ожидаемая эффективная доза, мЗв С-14; Медленный; 1.000; Ожидаемая эффективная доза, мЗв Н-3; Пероральный; Ожидаемая эффективная доза, мЗв Н-3; Медленный 1.000; Ожидаема я эффективна доза, мЗв
► 2020 1.1е-4 1.15е-3 7.4е-5 7.0Ge-4 2.4е-5
а 1
Рис. 1. Окно ввода входных данных ПМ РОЗА-Н. Введены данные для однократного облучения в 2020 г. по данным табл. 2.
Рис. 2. Окно с результатами вычислений ПМ РОЗА-Н. Показаны риски от внешнего облучения и поступления радионуклидов в организм в 2020 г. для женщин. Дозы облучения взяты из табл. 2.
Входные данные. Входные данные для ПМ РОЗА-Н делятся на два вида: те, что пользователь вводит в главном окне программы и те, что хранятся в виде файлов. Рассмотрим входные данные, которые пользователь вводит в главном окне программы (рис. 1). Вверху окна расположена кнопка «Редактировать входные данные». Нажатие на эту кнопку приводит к открытию окна проводника Windows на папке «Rates», расположенной в папке с исполняемым файлом ПМ РОЗА-Н. Здесь пользователь может отредактировать файлы, содержащие поло-
возрастные показатели фоновой онкозаболеваемости, онкосмертности, смертности от всех причин, половозрастную численность. Выбор нужных файлов осуществляется при помощи выпадающего списка с заголовком «Выберите входные данные». Этот список состоит из названий имён папок, которые содержатся в папке «Rates», расположенной в папке с исполняемым файлом ПМ РОЗА-Н. Данные о половозрастных показателях фоновой онкозаболеваемости, онкосмертности, смертности от всех причин и половозрастной численности будут считываться ПМ РОЗА-Н из выбранной пользователем папки.
Поле «Текущий календарный год» предназначено для ввода пользователем календарного года, который в ПМ РОЗА-Н будет считаться текущим, т.е. пожизненные риски будут рассчитываться из условия, что люди дожили до этого календарного года. По умолчанию в это поле заносится текущий календарный год.
Область «Тип расчёта» позволяет пользователю выбрать - какой тип риска нужно вычислять: онкозаболеваемость или онкосмертность.
Область «Пол» позволяет пользователю выбрать пол, для которого будут проводиться вычисления.
Область «Добавить облучение» содержит две кнопки «Внешнее» и «Внутреннее», которые позволяют пользователю добавить информацию о дозах внешнего и внутреннего облучения соответственно.
Кнопка «Вычислить риски» запускает расчёт радиационных рисков. Ниже кнопки «Вычислить риски» в главном окне программы расположена таблица, в которую пользователь вводит данные об облучении. Сначала пользователь должен добавить облучение: внутреннее или внешнее, нажав соответствующую кнопку. При этом, если в таблице не было строк, то появится одна строка. Для редактирования ячейки таблицы по ней нужно сделать двойной клик. В ячейки первого столбца пользователь должен ввести год облучения или диапазон календарных лет облучения через символ разделитель «-».
Если пользователь добавил внешнее облучение, то второй столбец должен называться «Эффективная доза внешнего облучения, мЗв». В ячейки этого столбца пользователь должен внести соответствующие данные. Столбцов для внешнего облучения в таблице будет столько, сколько раз пользователь нажмёт на кнопку «Внешнее» в разделе «Добавить облучение». Для удаления столбца с внешним облучением пользователю нужно нажать на заголовок столбца. В этом случае программа предложит удалить столбец.
Если пользователь добавил внутреннее облучение, то появится столбец с именем «Нажмите для выбора параметров внутреннего облучения». Для задания параметров внутреннего облучения пользователь должен кликнуть на этот заголовок таблицы. При этом появится окно, представленное на рис. 3.
В этом окне пользователь может выбрать радионуклид, от которого происходит внутреннее облучение, тип химического соединения этого радионуклида (возможные химические соединения представлены в табл. 1), путь поступления радионуклида в организм человека (ингаляция или поступление с пищей), в случае ингаляционного поступления - тип растворимости радионуклида и его AMAD, и дозу облучения в виде ожидаемой эффективной дозы в мЗв или (в случае ингаляционного поступления) объёмной активности радионуклида в воздухе или (в случае поступления с пищей) в виде активности на килограмм питьевой воды.
I Внутреннее аб... — X
Выберите радионуклид:
Выберите химическое соединение:
Выберите путь поступления радионуклида:
I Ингаляционный
Выберите тип растворимое™ (соединение):
Выберите АМАД:
1.000
Выберите размерность:
Удалить Подтвердить
Рис. 3. Окно задания параметров внутреннего облучения.
Таблица 1
Список радионуклидов ПМ РОЗА-Н, для которых имеется возможность выбора химического соединения при задании параметров внутреннего облучения
Радионуклид Тип соединения
C Двуокись углерода Моноксид углерода
H Органические соединения трития Элементарный водородный пар Тритиевая вода
Hg Пары ртути
I Йодистый метил Элементарный йод
Ni Карбонил никеля
Ru Тетроксид рутения
S Диоксид серы Сероуглерод
Te Теллуристый пар
После выбора параметров внутреннего облучения пользователь может нажать кнопку «Подтвердить», в этом случае в таблице в главном окне программы появится заголовок, содержащий в названии выбранные параметры. Или пользователь может нажать кнопку «Удалить» - это действие удалит в таблице соответствующий столбец для внутреннего облучения.
Пользователь может создать несколько строк для ввода данных. Для этого нужно кликнуть на строку, которая расположена в самом низу таблицы и выделена серым цветом. В каждой строке пользователь может вводить данные, соответствующие разным периодам облучения. Например, пусть в первой строке введён период облучения 2005-2009. Во второй строке пользователь может добавить данные об облучении, скажем, за 2013 г. и т.д. Для удаления строк следует нажать левой кнопкой мыши на поле левее строки, чтобы выделить эту строку и затем нажать на клавиатуре клавишу «Delete».
Рассмотрим входные данные, которые задаются в программу в виде файлов. Установочный пакет ПМ РОЗА-Н включает в себя следующие входные файлы.
1) Заранее подготовленные файлы с входными данными для Томской области и России в целом. Эти файлы содержат фоновые показатели онкозаболеваемости и онкосмертности на 2017 г. [8] по следующим локализациям: все солидные раки, желудок, толстая кишка, печень, лёгкое, молочная железа, яичник, мочевой пузырь, щитовидная железа, остальные солидные раки, лейкоз, а также численность населения и показатели фоновой смертности от всех причин.
Эти данные находятся в папке с именем «Rates», которая копируется при установке программы в каталог с исполняемым файлом. Данные для Томской области получены из «грубого» показателя, в предположении одинакового возрастного распределения показателей с такими данными по России в целом.
2) Файл NuclidesDataEPA.zip, который является zip-архивом, содержащим файлы со значениями дозовых коэффициентов для 721 радионуклида (нет благородных газов) по данным EPA [5]. Для каждого радионуклида представлена динамика дозовых коэффициентов от момента поступления радионуклида в организм с интервалом в один год. Дозовые коэффициенты даны для поступления радионуклида в организм человека с пищей, ингаляционному поступлению для 10 различных значений медианного по активности аэродинамического диаметра (AMAD), мкм: (0,001; 0,003; 0,01; 0,03; 0,1; 0,3; 1,0; 3,0; 5,0; 10,0) и трёх вариантов растворимости: медленная (S), средняя (M) и быстрая (F), а также для ингаляционного поступления химических соединений для некоторых радионуклидов.
3) Несколько файлов, содержащих варианты усреднения пожизненных радиационных рисков (дети и взрослые, по десятилетним и пятилетним интервалам). Это текстовые файлы, содержащие диапазоны возрастов для усреднения рисков. Они хранятся в папке «Avg», которая копируется при установке программы в каталог с исполняемым файлом.
Во время установки ПМ РОЗА-Н все указанные выше файлы копируются на компьютер пользователя.
Выходные данные. Результаты работы ПМ РОЗА-Н представляются в двух видах: на экране в окне программы и в виде файлов формата MS Excel.
В окне программы результаты вычислений пожизненных рисков представлены на двух вкладках: «Таблица» и «График». Вкладка «Таблица» содержит данные в табличном виде для каждого возраста в заданном календарном году (рис. 2). Столбцы содержат названия пожизненных рисков, значения которых были рассчитаны. Вкладка «График» позволяет пользователю отобразить результаты вычислений в виде графиков. Пользователь может выбрать нужные типы рисков (рис. 4).
Рис. 4. Окно ПМ РОЗА-Н, вкладка «График». Здесь отображаются графики выбранных пользователем радиационных рисков в зависимости от возраста в заданном календарном году. Показаны значения LAR от внешнего облучения и поступления радионуклидов в организм в 2020 г. для женщин. Дозы облучения взяты из табл. 2.
Программа позволяет проводить экспорт текущих результатов вычислений в файл формата MS Excel. Для этого пользователь должен нажать кнопку «Экспорт результатов в Excel». В этом случае будет запрошен путь для сохранения «.xlsx» файла. Для сохранения результатов в файл MS Excel на компьютере пользователя должно быть установлено соответствующее программное обеспечение. В выходном файле после экспорта будут содержаться значения пожизненных рисков. Название риска указано в заголовках столбцов. Кроме того, будут построены графики этих рисков от возраста в заданном календарном году. Примеры выходного файла представлены на рис. 5-7, на них даны листы «Входные данные», «Все возрасты» и «По десятилетним интервалам» соответственно.
А А В с
1 Текущий календарный год: 2020
2 Облучение: внешнее и внутреннее
3 4 Локализация: Дозы: Все ЗНО
5 Календарный год Эффективная доза внешнего облучения (мЗв) Внутреннее облучение
6 С-14, поступление с пищей, ожидаемая эффективная доза (мЗв)
7 2020 0.00011 0.00115
8
Рис. 5. Пример файла MS Excel, куда были экспортированы результаты вычислений.
Показан лист «Входные данные».
А А с Е F G H 1 J 1 К M N 0 Cl
1 Возраст в 2020 году LAR LBR LARF, % REID YLL, годы YLLRIC
3.00E-07
3 1 2.35Е-07 0.13106 0.00% 2.35 Е-07 4.13Е-06 17.60502
4 2 2.29Е-07 0.131178 0.00% 2.29 Е-07 3.96Е-06 17.3106
5 3 2.22Е-07 0.131296 0.00% 2.22Е-07 3.86Е-06 17.35676 2.50E-07
ft 4 2.17Е-07 0.131413 0.00% 2.17Е-07 3.77Е-06 17.40993
7 5 2.11Е-07 0.131531 0.00% 2.11 Е-07 3.69Е-06 17.46876 2.00E-07
3 6 2.05Е-07 0.131525 0.00% 2.05 Е-07 З.бЕ-Об 17.52007
9 7 1.97Е-07 0.131519 0.00% 1.97Е-07 3.41Е-06 17.30836
10 8 1.92Е-07 0.131513 0.00% 1.92Е-07 З.ЗЗЕ-Об 17.35149 Ч 1.50 Е-07
11 9 1.87Е-07 0.131507 0.00% 1.87Е-07 3.25Е-06 17.38262
1? 10 1.82Е-07 0.1315 0.00% 1.82Е-07 3.17Е-06 17.4144 1.00Е-07
13 11 1.78Е-07 0.131505 0.00% 1.78 Е-07 3.1Е-06 17.44581
14 12 1.70Е-07 0.131509 0.00% 1.70Е-07 2.94Е-06 17.25635
15 13 1.66Е-07 0.131513 0.00% 1.66Е-07 2.87Е-06 17.27634 5.00Е-08
16 14 1.62Е-07 0.131517 0.00% 1.62Е-07 2.79Е-06 17.275??
17 15 1.58Е-07 0.131521 0.00% 1.58Е-07 2.71Е-06 17.17006 0.00 Е+00
18 1& 1.54Е-07 0.131542 0.00% 1.54Е-07 2,53 Е-06 17.09866 SO so 100
19 17 1.50Е-07 0.131563 0.00% 1.50Е-07 2.53 Е-Об 16.88275 раст в 2020 году
?0 18 1.46Е-07 0.131584 0.00% 1.46Е-07 2.46Е-06 16.80119
21 19 1.43 Е-07 0.131605 0.00% 1.43 Е-07 2.38Е-06 16.70411
24 22 1.33 Е-07 0.131679 0.00% 1.33 Е-07 2.15 Е-06 16.11675
Рис. 6. Пример файла MS Excel, куда были экспортированы результаты вычислений. Показан лист «Все возрасты», а также разные виды пожизненных рисков после внешнего облучения и поступления радионуклидов в организм в 2020 г. для женщин. Дозы облучения взяты из табл. 2
А В С D Е F G H 1 J L M MOP 1 а
1 Диапазон возрастов в 2020 году LAR LBR LARF, % REID YLL, ГОДЫ YLLRIC
2 0-4 2.29Е-07 0.131178 0.00% 2.29Е-07 4.01Е-06 17.50952
3 5-9 1.98Е-07 0.131519 0.00% 1.98Е-07 3.45Е-06 17.40866
4 10-14 1.72Е-07 0.131509 0.00% 1.72Е-07 2.97Е-06 17.32709 \
5 15-19 1.50Е-07 0.131563 0.00% 1.50Е-07 2.54Е-06 16.93767 2.00E-07 --
6 20-24 1.33Е-07 0.131679 0.00% 1.33 Е-07 2.15Е-06 16.11631 \
7 25-29 1.18Е-07 0.131813 0.00% 1.18Е-07 1.81Е-06 15.361
8 30-34 1.06Е-07 0.131906 0.00% 1.06Е-07 1.53Е-06 14.52652 1.50E-07 \
9 35-39 9.55Е-08 0.131784 0.00% 9.55Е-08 1.29Е-06 13.51572 EÉ X.
10 40-44 9.47Е-08 0.131089 0.00% 9.47Е-08 1.15Е-06 12.12583
11 45-49 8.13Е-08 0.129394 0.00% 8.13Е-08 8.97Е-07 11.03907 1.00E-07 —».
12 50-54 6.90Е-08 0.126283 0.00% 6.90Е-08 6.82Е-07 9.884999
13 55-59 5.75Е-08 0.121211 0.00% 5.75Е-08 4.99Е-07 8.673302
14 60-64 4.64Е-08 0.113549 0.00% 4.64Е-08 3.45Е-07 7.42371 5.0QE-Q8
15 65-69 3.55Е-08 0.1028 0.00% 3.55Е-08 2.2Е-07 6.203462
16 70-74 2.48Е-08 0.088991 0.00% 2.48Е-08 1.26Е-07 5.071838 0 00E+00
17 18 75-79 80-100 1.54Е-08 5.81Е-09 0.072813 0.045615 0.00% 0.00% 1.54Е-08 5.81Е-09 6.41Е-08 1.74Е-08 4.166598 2.990812 ■ï- ,f ■
19 20 Диапазон возрастов в 2020 году
?1
Рис. 7. Пример файла MS Excel, куда были экспортированы результаты вычислений. Показан лист «По пятилетним интервалам», а также разные виды усреднённых пожизненных рисков после внешнего облучения и поступления радионуклидов в организм в 2020 г. для женщин.
Дозы облучения взяты из табл. 2.
Пример проведения расчёта при помощи РОЗА-Н. Рассмотрим пример расчёта с помощью ПМ РОЗА-Н величины LAR от однократного облучения дозой 5 мЗв. Усреднённые по полу и возрасту значения LAR смертности от ЗНО и LAR заболеваемости ЗНО равны 1,61 х10-4 и 3,09х10-4 соответственно. Оценка радиационного риска с использованием номинального коэффициента риска НРБ-99/2009 [7] 5,5х10-2 Зв-1 при дозе облучения 5 мЗв даёт величину риска 2,75х10-4, что в 1,71 раза выше рассчитанного ПМ РОЗА-Н среднего LAR смертности и в 1,12 раза ниже рассчитанного ПМ РОЗА-Н среднего LAR заболеваемости. Такие оценки LAR с помощью ПМ РОЗА-Н следует признать корректными, так как номинальный коэффициент риска НРБ-99/2009 [7] представляет собой удельный радиационный вред (на единицу дозы), который был получен путём специальной процедуры усреднения коэффициентов пожизненного радиационного риска онкосмертности и онкозаболеваемости [1].
Рассмотрим пример расчёта пожизненных радиационных рисков после годового облучения населения, обусловленного выбросами предприятий опытного демонстрационного энергокомплекса (ОДЭК).
По данным Е.В. Спирина [9] при одновременной работе всех предприятий ОДЭК годовая доза облучения человека от выбросов изотопов 3H и 14C равна 0,73 и 1,21 мкЗв соответственно, что составляет примерно половину суммарной дозы (3,62 мкЗв) от всех выбрасываемых радионуклидов. Эффективная доза обусловлена, в основном, их пероральным поступлением внутрь организма. Облучение за счёт ингаляционного поступления значительно меньше, доза от 3H и 14C составляет 0,024 и 0,074 мкЗв/год соответственно. Радионуклидный состав выбросов ОДЭК с соответствующими ожидаемыми эффективными дозами отдельно по радионуклидам приведён в табл. 2.
Таблица 2
Предполагаемый радионуклидный состав, пути поступления радионуклидов и эффективные дозы в результате нормальной работы ОДЭК
Радионуклид
Тип поступления радионуклида
ОЭД, мЗв/год
Внешнее облучение
14C 14C 3H 3H
238Pu 239Pu 240Pu 241Pu 242Pu 210Po 241Am 244Cm
Пероральное Ингаляционное Пероральное Ингаляционное Ингаляционное Ингаляционное Ингаляционное Ингаляционное Ингаляционное Пероральное Пероральное Пероральное
1,10E-04 1,15E-03 7,40E-05 7,06E-04 2,40E-05 7,21 E-05 1,54E-04 1,38E-04 1,25E-04 3,70E-07 1,70E-04 1,14E-04 7,86E-04
В расчёте примем, что для ингаляционного поступления изотопов углерода и водорода скорость сорбции будет медленная, а для изотопов плутония - быстрая. Для ингаляционного поступления примем AMAD равным 1.0. Результаты расчёта при помощи ПМ РОЗА-Н усреднённых по десятилетним возрастным группам пожизненных рисков возникновения ЗНО представлены в табл. 3.
Таблица 3
Пожизненные атрибутивные риски смерти от ЗНО для населения, проживающего
вблизи предприятий ОДЭК
Возрастная группа Мужчины Женщины Оба пола в среднем
0-9 1,76E-07 2,14E-07 1,96E-07
10-19 1,39E-07 1,61E-07 1,51E-07
20-29 1,04E-07 1,24E-07 1,15E-07
30-39 8,44E-08 1,01E-07 9,33E-08
40-49 6,79E-08 8,83E-08 7,88E-08
50-59 5,02E-08 6,27E-08 5,69E-08
60-69 3,44E-08 4,14E-08 3,82E-08
70-79 2,00E-08 1,92E-08 1,96E-08
80-89 1,18E-08 8,11E-09 9,82E-09
90-100 3,25E-09 8,50E-10 1,96E-09
Данные табл. 3 показывают, что при нормальной эксплуатации ОДЭК вклад в LAR населения от годового поступлении радионуклидов в организм за счёт предприятий ОДЭК является пренебрежимо малым (менее 10-6 [7]).
Усреднённый LAR по всем возрастам и полу в этом случае равен 9,28х10-8. Суммарная
-3
ОЭД в рассматриваемом примере равна 3,62х10- мЗв в год. Произведение этой дозы на коэффициент риска ЗНО из НРБ-99/2009 [7] 5,5х10- на мЗв даст значение пожизненного риска, равное 1,99х10-7. Это в 2,15 раза больше среднего риска, рассчитанного напрямую по моделям МКРЗ. Также точный расчёт позволяет увидеть сильную зависимость LAR от возраста при облучении. Из табл. 3 видно, что только риск для возрастной группы 0-9 лет примерно совпадает с риском, полученным через коэффициент риска из НРБ-99/2009, а уже в возрасте 30-39 лет фактический риск ниже в 2 раза.
Обсуждение
В разработанном программном модуле применяется метод оценки радиационных рисков на основе моделей, предложенных в основном МКРЗ [1]. В то же время существуют и другие международные модели радиационных рисков, предложенные, например, НКДАР ООН [2] или ВОЗ [10]. Так, НКДАР ООН [2] даёт модели рисков для таких локализаций опухолей, для которых не даёт модели МКРЗ [1], например кость, немеланомный рак кожи и мозг. В то же время НКДАР ООН [2] не даёт модели для локализации яичники.
В будущих версиях ПМ РОЗА-Н будет рассмотрена возможность проведения расчётов по моделям риска, отличным от моделей МКРЗ, а в текущей версии ПМ РОЗА-Н было принято решение учитывать влияние внутреннего облучения на развитие рака кости, поэтому для этой локализации использовалась модель, предложенная НКДАР ООН [2].
Ещё одним существенным отличием моделей, предложенных другими международными организациями, от моделей МКРЗ [1] является вопрос об использовании DDREF. Модели МКРЗ [1] являются менее консервативными, чем модели НКДАР ООН [2] и ВОЗ [10], в которых DDREF не используется.
Общей особенностью международных моделей радиационного риска является то, что они в основном идентифицировались по данным когорты выживших после атомной бомбардировки Хиросимы и Нагасаки в 1945 г. Для оценки применимости этих моделей для расчёта радиационных рисков для жителей Российской Федерации требуются дополнительные исследования.
Также следует отметить, что сами модели риска МКРЗ [1] не дают точного результата. Параметры этих моделей имеют неопределённости [11, 12]. Также на неопределённость результатов расчёта влияет неопределённость доз облучения и неопределённость фоновых значений онкозаболеваемости и онкосмертности. ПМ РОЗА-Н в настоящее время никак не учитывает эти неопределённости. В ПМ РОЗА-Н для расчётов применяются фоновые показатели за фиксированный календарный год. Однако при прогнозе на будущее такой подход вносит погрешность, т.к. фоновые показатели не являются статичными, а меняются во времени.
Одним из вкладов в неопределённость рисков от внутреннего облучения может являться то, что сами значения дозовых коэффициентов, полученные из программного комплекса DCAL [5], не являются точными. Этот фактор требует дополнительного исследования.
Заключение
В целях разработки программного обеспечения для расчёта радиационных рисков населения от воздействия ядерных объектов в замкнутом ядерном топливном цикле на основе быстрых реакторов разработан ПМ РОЗА-Н, предназначенный для оценки радиологической опасности в условиях нормальной работы и в аварийных ситуациях, с учётом современных международных рекомендаций и национальных требований к радиационной защите населения.
Разработанный программный модуль РОЗА-Н впервые в атомной отрасли реализует расчёт пожизненных радиационных рисков населения, выраженных шестью различными метриками риска: пожизненный атрибутивный риск (LAR); риск радиационно-индуцированного случая заболевания (REIC); риск радиационно-индуцированного случая смерти (REID); пожизненная атрибутивная доля (LARF); потерянные годы жизни и потерянные годы здоровой жизни (YLL); потерянные годы жизни, приведённые на радиационно-индуцированный случай онкозаболевания или онкосмерти (YLLRIC).
Оценка радиационного риска с использованием номинального коэффициента риска НРБ-99/2009 может значительно (до 2 раз) превышать оценку риска, рассчитанную непосредственно по моделям риска, рекомендуемым в настоящее время МКРЗ.
Практическая значимость ПМ РОЗА-Н определяется тем, что этот расчётный модуль позволяет получать оценки радиационных рисков любых половозрастных групп населения от выбросов любого радионуклидного состава, обеспечивая оперативное планирование и оптимизацию мероприятий по защите населения.
Литература
1. 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.
2. 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.
3. Меняйло А.Н., Чекин С.Ю., Кащеев В.В., Максютов М.А., Корело А.М., Туманов К.А., Пряхин Е.А., Ловачев С.С., Карпенко С.В., Кащеева П.В., Иванов В.К. Пожизненный радиационный риск в результате внешнего и внутреннего облучения: метод оценки //Радиация и риск. 2018. Т. 27, № 1. С. 8-21.
4. 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.2020).
5. DCAL Software and Resources. [Электронный ресурс]. URL: https://www.epa.gov/radiation/dcal-software-and-resources (дата обращения 13.04.2020).
6. ICRP, 2012. Compendium of dose coefficients based on ICRP Publication 60. ICRP Publication 119 //Ann. ICRP. 2012. V. 41 (Suppl.). P. 1-130.
7. Нормы радиационной безопасности (НРБ-99/2009). Санитарные правила и нормативы. СанПиН 2.6.1.2523-09. М.: Федеральный центр гигиены и эпидемиологии Роспотребнадзора, 2009. 100 с.
8. Злокачественные новообразования в России в 2017 году (заболеваемость и смертность) /Под ред. А.Д. Каприна, В.В. Старинского, Г.В. Петровой. М.: МНИОИ им. П.А. Герцена, 2018.
9. Спирин Е.В., Алексахин Р.М., Бажанов А.А. Структура дозы облучения населения при эксплуатации предприятий опытного демонстрационного энергокомплекса //Атомная энергия. Т. 124, №. 3. С. 169-173.
10. Health risk assessment from the nuclear accident after the 2011 Great East Japan earthquake and tsunami based on a preliminary dose estimation. Geneva: WHO, 2013. 172 p.
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 radiological protection of the public
Menyajlo A.N., Chekin S.Yu., Lovachev S.S., Tumanov K.A., Korelo A.M., Ivanov V.K
A. Tsyb MRRC, Obninsk
The article presents the newly-developed software module "PM ROSA-N" for radiological protection of the public. The softWare was designed for calculating radiation risks for people living close to nuclear power facilities and enterprises using radioactive materials or sources of ionizing radiation. When selecting the methodology authors followed current recommendations of the International Commission on Radiological Protection (ICRP) and the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR). Evaluation of internal radiation doses is based on dose factors derived from the data of the United States Environmental Protection Agency (EPA). The authors present the application of the "PM ROSA-N" for estimating lifetime attributable risk (LAR) of cancer mortality for the general public living near the facilities of the Pilot-demonstration energy complex, the new generation of the power system ("Proryv" (Breakthrough) Project). Estimates obtained with the use of the software module show considerable, tenfold, reduction of radiation risk in relation to age at radiation exposure; it varies from 2.1 x10"7 for children under 10 years to 1.96x10"8 in the age group of 70-79 years. The risk estimate averaged by age and sex is 1.00x10"7, it is two times less than the estimate obtained with the use of risk nominal coefficient 1.99x10-7, given in the Russian radiation safety standards (NRB-99/2009). The results obtained show that the newly developed software module the "PM ROSAN" is the new tool that ensures well-founded radiation safety of the public living close to nuclear power facilities with reference to the current international safety standards. The software module ROSA-N, that provides justification of safety and protective measures for the public in response to radionuclides emission, governs its economic efficiency and importance. The module can be used for estimating risk for any person living in the close proximity to a nuclear facility. It corresponds to the current enforcement practice used in the developed countries.
Key words: radiation safety, radiation risk, radiobiological effects, radio-induced oncological diseases, radiation risk models of ICRP and UNSCEAR, dose coefficients, route of radionuclide intake, external exposure, internal exposure, radionuclides, AMAD, expected effective dose equivalent dose.
References
1. 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.
2. 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.
3. Menyajlo A.N., Chekin S.Yu., Kashcheev V.V., Maksioutov IW^., 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 and internal exposure to radiation: method for estimating. Radiatsiya i risk - Radiation and Risk, 2018, vol. 27, no. 1, pp. 8-21. (In Russian).
4. 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.2020).
5. DCAL Software and Resources. Available at: https://www.epa.gov/radiation/dcal-software-and-resources (Accessed 13.04.2020).
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.; Korelo A.M. - Senior Researcher; Ivanov V.K. - Deputy Director, Chairman of RSCRP, Corresponding Member of RAS. A. Tsyb MRRC.
•Contacts: 4 Korolyov str., Obninsk, Kaluga region, Russia, 249035. Tel.: (484) 399-32-81; e-mail: nrer@obninsk.com.
6. ICRP, 2012. Compendium of dose coefficients based on ICRP Publication 60. ICRP Publication 119. Ann. ICRP, 2012, vol. 41 (Suppl.), pp. 1-130.
7. Norms of radiation safety (NRS-99/2009). Sanitary rules and regulations. SanPiN 2.6.1.2523-09. Moscow, Federal Center for Hygiene and Epidemiology of Rospotrebnadzor, 2009. 100 p. (In Russian).
8. 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).
9. Spirin E.V., Aleksakhin R.M., Bazhanov A.A. Composition of total radiation dose to public during the Pilot demonstration energy complex operation. Atomnaya energiya - Atomic Energy, 2018, vol. 124, no. 3, pp. 169-173. (In Russian).
10. Health risk assessment from the nuclear accident after the 2011 Great East Japan earthquake and tsunami based on a preliminary dose estimation. Geneva, WHO, 2013. 172 p.
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.
