МИНИСТЕРСТВО ЗДРАВООХРАНЕНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

ФЕДЕРАЛЬНОЕ УПРАВЛЕНИЕ МЕДИКО-БИОЛОГИЧЕСКИХ И ЭКСТРЕМАЛЬНЫХ ПРОБЛЕМ

2.6.1. Ионизирующее излучение, радиационная безопасность

РАСЧЕТ ОЖИДАЕМЫХ ЭФФЕКТИВНЫХ ДОЗ ВНУТРЕННЕГО ОБЛУЧЕНИЯ ПЕРСОНАЛА ПО РЕЗУЛЬТАТАМ ИЗМЕРЕНИЙ АКТИВНОСТИ РАДИОНУКЛИДОВ В БИОПРОБАХ С ИСПОЛЬЗОВАНИЕМ КОМПЬЮТЕРНОЙ ПРОГРАММЫ ММК-01

Методика выполнения расчетов

МВР 2.6.1.60-2002

 

Содержание

Предисловие

1. Методика выполнения расчетов МВР 2.6.1.60-2002 «Расчет ожидаемых эффективных доз внутреннего облучения персонала по результатам измерений активности радионуклидов в биопробах с использованием компьютерной программы ММК-01» разработаны авторским коллективом специалистов.

Исполнители: к.т.н. Молоканов А.А. (ГНЦ - ИБФ),

к.м.н. Антипин Е.Б. (ФУ «Медбиоэкстрем»),

д.т. н. Бадьин В.И (БФАИИЧ ГЦ ГСЭН),

к.ф - м.н. Крючков В.П. (ГНЦ- ИБФ),

зав. лаб. Цовьянов А.Г. (ГНЦ - ИБФ),

к.т.н. Ермилов А.П (ЦМИИ ГП ВНИИФТРИ),

к.т.н. Антропов СЮ (ЦМИИ ГП ВНИИФТРИ).

2. Утверждены и введены в действие Федеральным управлением медико-биологических и экстремальных проблем (Федеральное Управление «Медбиоэкстрем») при Минздраве России 4 ноября 2002 г.

3. В настоящих методических указаниях реализованы требования законов Российской Федерации:

«О санитарно-эпидемиологическом благополучии населения» ФЗ-52 от 30.03.1999 г.;

«О радиационной безопасности населения» ФЗ-3 от 09.01.1996 г.;

«Об использовании атомной энергии» ФЗ-170 от 21.11.1995 г.;

«Об обеспечении единства измерений» 487-1 от 27.04.1993 г.;

«Об информации, информатизации и защите информации» ФЗ-24 от 20.02.1995 г.;

«Об обеспечении единства измерений» № 4871-1 от 27.04.1993 г.;

«О стандартизации» № 5154-1 от 10.06.1993 г.

4. Введены впервые.

Согласованы с Директором Центра метрологии ионизирующих излучений ГП «ВНИИФТРИ» В.П. Ярыной 18 мая 2001 г.

Утверждены Заместителем Главного государственного санитарного врача по объектам и территориям, обслуживаемым Федеральным управлением «Медбиоэкстрем» О.И. Шамовым 04 ноября 2002 г.

1. Назначение и область применения

1.1. Методика выполнения расчетов (далее МВР) совместно с компьютерной программой ММК-01 (далее ММК-01) предназначена для обеспечения требований Норм радиационной безопасности (НРБ-99, п.7) и последующих Методических указаний (МУ 2.6.1.16-2000 и МУ 2.6.1.26-2000) к контролю уровней облучения персонала от источников внутреннего облучения.

1.2. Данная МВР совместно с ММК-01 обеспечивает определение индивидуальных ожидаемых эффективных доз (ОЭД) внутреннего облучения персонала на основе измерений, проведенных в рамках индивидуального дозиметрического контроля (ИДК).

1.3. Областью применения данной МВР является расчет индивидуальных ОЭД, обусловленных ингаляционным поступлением в организм радионуклидов в стандартных условиях облучения, в котором в качестве исходных данных используются результаты индивидуальных измерений активности радионуклидов (³Н, ⁵⁹Fe, ⁵⁷Co, ⁵⁸Co, ⁶⁰Co, ⁸⁵Sr, ⁸⁹Sr, ⁹⁰Sr, ¹⁰⁶Ru, ¹²⁵I, ¹²⁹I, ¹³¹I, ¹³⁴CS, ¹³⁷CS, ²²⁶Ra, ²²⁸Ra, ²²⁸Th, ²³²Th, ²³⁴U, ²³⁵U, ²³⁸U, ²³⁷Np, ²³⁸Pu, ²³⁹Pu, ²⁴⁰Pu, ²⁴¹Am, ²⁴²Cm, ²⁴⁴Cm, ²⁵²Cf) в биопробах (моча, кал) человека.

1.4. ММК-01 на основе результатов измерений активностей из рабочей базы данных программно-аппаратурного комплекса «Прогресс» или из собственной базы данных проводит расчеты индивидуальных ожидаемых эффективных доз облучения работников за любой календарный год и выдает информацию по уровню внутреннего облучения работника в виде протокола и графических иллюстраций

2. Нормативные ссылки

2.1. В настоящих Методических указаниях использованы основные положения следующих руководящих документов:

• СП-2.6.1.758-99. Нормы радиационной безопасности (НРБ-99). Гигиенические нормативы. М.: Центр санитарно-эпидемиологического нормирования, гигиенической сертификации и экспертизы Минздрава России, 1999.

• МУ 2.6.1.16-2000. Определение индивидуальных эффективных и эквивалентных доз и организация контроля профессионального облучения в контролируемых условиях обращения с источниками излучения. Общие требования.

• МУ 2.6.1.26-2000 Дозиметрический контроль профессионального внутреннего облучения. Общие требования.

3. Термины, определения и условные обозначения

3.1 Термины и определения.

Аэрозоль - дисперсная система с газообразной средой и с твердой, жидкой или смешанной дисперсной фазой.

Аэрозоль радиоактивный - аэрозоль, в дисперсную фазу которого входят радионуклиды.

Аэродинамический диаметр частицы аэрозоля - это диаметр частицы с плотностью, равной 1 г/см³, имеющей ту же конечную скорость оседания в воздухе при нормальных условиях, что и данная частица.

Активностный медианный аэродинамический диаметр (AMAD) - такое значение аэродинамического диаметра частиц дисперсной фазы радиоактивного аэрозоля, при котором 50 % активности указанного аэрозоля приходится на частицы, имеющие диаметр меньше, чем AMAD, а 50 % - на частицы, имеющие аэродинамический диаметр больше, чем AMAD.

Доза индивидуальная эффективная (эквивалентная в органе или ткани) - эффективная доза (эквивалентная доза в органе или ткани), которая была бы получена стандартным работником, если бы он находился в тех же производственных условиях и выполнял те же работы с источником, что и данный индивид. Значение индивидуальной дозы приписывается индивиду по результатам дозиметрического контроля.

Доза эффективная ожидаемая при внутреннем облучении (Е(τ)) - величина, используемая как мера риска возникновения отдаленных последствий внутреннего облучения всего тела человека и отдельных его органов и тканей с учетом их радиочувствительности. Она представляет сумму произведений ожидаемой эквивалентной дозы в органах и тканях на соответствующие взвешивающие коэффициенты:

                                                                             (5)

Значение τ принимают для персонала равным 50 лет. Единица ожидаемой эффективной дозы - зиверт (Зв).

Доза эффективная (эквивалентная) годовая - сумма эффективной (эквивалентной) дозы внешнего облучения, полученной за календарный год, и ожидаемой эффективной (эквивалентной) дозы внутреннего облучения, обусловленной поступлением в организм радионуклидов за этот же год. Единица годовой эффективной дозы - зиверт (Зв).

Контроль дозиметрический индивидуальный (ИДК) - контроль облучения персонала, заключающийся в определении индивидуальных доз облучения работника на основании результатов индивидуальных измерений характеристик облучения тела или отдельных органов каждого работника.

Облучение внутреннее - облучение органов и тканей человека в результате поступления радионуклидов в организм человека.

Период контроля - промежуток времени между последовательными индивидуальными измерениями характеристик облучения каждого работника при проведении ИДК.

Поступление радионуклида ингаляционное - активность радионуклида, которая проникает в организм через органы дыхания.

Работник стандартный - воображаемый человек, обладающий биологическими и физическими свойствами, присущими среднестатистическому здоровому взрослому человеку. Свойства стандартного работника включают:

• антропометрические характеристики тела, отдельных органов и тканей человека;

• характеристики физиологических показателей человека;

• параметры биокинетики химических элементов в органах и тканях человека, рекомендованные МКРЗ и использованные при определении значений допустимых уровней облучения, установленных в Нормах.

Стандартная модель - модель определения индивидуальных доз, используемая в рамках индивидуального дозиметрического контроля и включающая:

• определение индивидуальных доз облучения по результатам индивидуальных систематических измерений характеристик облучения работника;

• распространение на объект контроля стандартных условий облучения.

Специальная модель - модель определения индивидуальных доз, используемая в рамках индивидуального дозиметрического контроля при определенных условиях (Е(τ) > У_Д) и включающая:

• определение индивидуальных доз облучения по результатам индивидуальных систематических измерений характеристик облучения работника;

• учет реальных условий облучения при интерпретации результатов индивидуальных систематических измерений характеристик облучения работника.

Тип химического соединения при ингаляции - категория дисперсной фазы радиоактивного аэрозоля в классификации по скорости перехода радионуклида из легких в кровь, установленной в дозиметрической модели органов дыхания МКРЗ:

• тип «М» (медленно растворимые соединения): при растворении в легких человека веществ, отнесенных к этому типу, наблюдается компонента активности радионуклида, поступающая в кровь со скоростью 0,0001 сут^-1;

• тип «П» (соединения, растворимые с промежуточной скоростью): при растворении в легких человека веществ, отнесенных к этому типу, основная активность радионуклида поступает в кровь со скоростью 0,005 сут^-1;

• тип «Б» (быстро растворимые соединения): при растворении в легких человека веществ, отнесенных к этому типу, основная активность радионуклида поступает в кровь со скоростью 100 сут^-1.

Уровень введения индивидуального дозиметрического контроля (У_ВК) - такое значение годовой эффективной (ожидаемой) дозы, при действительном или предполагаемом превышении которого определение соответствующих доз следует проводить с помощью индивидуального дозиметрического контроля облучения работника.

Уровень действия (У_Д) - такое значение дозы, при действительном или предполагаемом превышении которого следует провести мероприятия по улучшению радиационной обстановки на рабочем месте, а при расчете дозы следует применять специальную модель определения индивидуальных доз.

Уровень исследования (У_И) - такое значение дозы, полученной в течение периода контроля, при превышении которого следует провести исследование причин повышения дозы и при необходимости провести мероприятия по улучшению радиационной обстановки на рабочем месте.

Условия облучения стандартные - определенные в Нормах для целей нормирования техногенного облучения условия воздействия техногенных источников на человека, которые характеризуются следующими параметрами:

• объемом вдыхаемого воздуха, с которым радионуклид может поступить в органы дыхания персонала группы А на протяжение календарного года: V_ПЕРС = 2,4 · 10³ м³;

• классификацией дисперсной фазы радиоактивного аэрозоля по скорости перехода радионуклида из легких в кровь согласно п. 8.3 Норм. Распределение химических соединений элементов по типам при ингаляции приведено в приложении П-3 к Нормам;

• логарифмически нормальным распределением активности по размерам частиц дисперсной фазы аэрозоля с активностным медианным аэродинамическим диаметром (АМАД), равным 1 мкм и стандартным геометрическим отклонением, равным 2,5.

Функция выведения радионуклида - функция, описывающая активность в суточном количестве мочи или кала человека, в зависимости от времени, прошедшем после однократного единичного поступления радионуклида в организм.

3.2. Условные обозначения и определения.

t_i - дата i-го измерения активности радионуклида; (t_i-1 - t_i) - i-й период контроля;

А_i - активность радионуклида в суточном количестве мочи или кала человека, измеренная в момент времени t_i;

ΔА_i - абсолютная погрешность измерения активности радионуклида в суточном количестве мочи или кала человека;

E_i - ожидаемая эффективная доза (ОЭД) на момент времени t_i;

П_i - суммарное ингаляционное поступление радионуклида на момент времени t_i;

I_j(t) - случайная функция поступления радионуклида за j-ый период контроля (скорость ингаляционного поступления радионуклида в момент времени t, Бк/с);

R(t) - функция выведения радионуклида из организма человека для рассматриваемых условий облучения;

ε_R - дозовый коэффициент перехода от поступления радионуклида к ОЭД для рассматриваемых условий облучения;

f - параметр, характеризующий вариабельность коэффициента (функции) выведения радионуклида с мочой или калом, которая наблюдается у реального человека;

АМАД - активностный медианный аэродинамический диаметр;

ИДК - индивидуальный дозиметрический контроль; МВИ - методика выполнения измерений; ОЭД - ожидаемая эффективная доза.

4. Общие положения

4.1. МВР совместно с ММК-01 реализует стандартную модель определения индивидуальных доз в системе ДК внутреннего облучения персонала:

Как видно из рисунка, система ДК внутреннего облучения персонала состоит из 4-х частей - измерительной, расчетной, организационной, информационной и 2-х уровней - группового дозиметрического контроля (ГДК), осуществляющего контроль доз по результатам измерений загрязнения воздуха помещений (рабочих мест), и индивидуального дозиметрического контроля (ИДК), осуществляющего контроль доз по результатам измерений активности радионуклидов в человеке. В свою очередь ИДК делится на 3 подуровня - «стандартный», «специальный» и «индивидуальный». Разбиение на уровни и подуровни служит для оптимизации затрат на определение величины индивидуальной дозы. Последовательность действий по определению индивидуальных доз персонала отражается в Регламенте ДК.

Расчетная часть позволяет связать измеряемые величины (обеспечиваемые измерительной частью) с результирующей индивидуальной ОЭД. Она позволяет решать задачу оптимизации системы ДК исходя из возможностей имеющейся в наличии измерительной части и системой установленных допустимых уровней и неопределенностей оценки индивидуальных доз, являясь, таким образом, инструментом для разработки Регламента ДК внутреннего облучения персонала предприятия.

Стандартная модель определения индивидуальных доз в системе ДК внутреннего облучения персонала является базовой моделью расчета, достаточной для разработки Регламента ИДК.

4.2. Стандартная модель, согласно МУ 2.6.1.16-2000 (п. 10.1), используется в рамках ИДК персонала в контролируемых условиях (НРБ-99, п. 3), когда значение индивидуальной дозы (полученное с использованием стандартной модели) не превышает уровня действия (УД).

4.3. Стандартная модель включает:

• определение индивидуальных доз облучения по результатам индивидуальных систематических измерений активности радионуклидов в моче или кале человека;

• распространение на объект контроля стандартных условий облучения.

4.4. Стандартные условия облучения характеризуются следующими параметрами:

• объемом вдыхаемого воздуха, с которым радионуклид может поступить в органы дыхания персонала группы А на протяжении календарного года: V_ПЕРС = 2,4 · 10³ м³;

• классификацией дисперсной фазы радиоактивного аэрозоля по скорости перехода радионуклида из легких в кровь согласно п. 8.3 НРБ-99 и распределением химических соединений элементов по типам при ингаляции согласно приложению П-3 к НРБ-99.

• логарифмически нормальным распределением активности по размерам частиц дисперсной фазы аэрозоля с активностным медианным аэродинамическим диаметром (АМАД), равным 1 мкм и стандартным геометрическим отклонением, равным 2,5.

4.5. Расчет доз проводят в предположении, что каждый человек соответствует характеристикам стандартного работника. При этом игнорируется возможное отличие действительной индивидуальной дозы от расчетной, обусловленное различием между характеристиками стандартного человека и биологическими характеристиками индивида.

4.6. Определение величины индивидуальной ОЭД (Е) и оценка неопределенности расчета дозы (U) проводится путем математического моделирования процесса поступления и выведения радионуклидов из организма стандартного человека и использования метода статистических испытаний (метода Монте-Карло) для нахождения распределения оценок величин Е и его параметров.

4.7. При определении величины индивидуальной ОЭД и оценке неопределенности расчета дозы учитывают:

• случайный характер поступления радионуклида в органы дыхания и неопределенность момента поступления радионуклида за период контроля;

• погрешность измерений активности радионуклидов в моче или кале человека согласно МВИ;

• неопределенность используемого биофизического метода измерения;

• величины активности радионуклида в моче или кале, обусловленные поступлением радионуклида в организм работника в прошлом.

5. Алгоритм расчета поступления и ОЭД и оценка их погрешности

5.1. Расчет поступления и ОЭД проводится с использованием программы ММК-01.

5.2. Исходными данными для расчета поступления радионуклида и ОЭД служат:

• результаты систематических измерений активности радионуклида в суточном количестве мочи или кала человека, А_mi;

• даты измерений, t_i;

• функция выведения радионуклида из организма человека для рассматриваемых условий облучения, R(t);

• дозовый коэффициент перехода от поступления радионуклида к ОЭД для рассматриваемых условий облучения, ε_R.

5.3. Оценка ингаляционного поступления радионуклида и ОЭД проводится в соответствии со следующей моделью:

5.3.1. Активность радионуклида в моче или кале человека обусловлена ингаляционным поступлением данного радионуклида за рассматриваемый и за все предыдущие периоды контроля и связана с поступлением соотношением:

                                           (5.1)

где i = 1, 2,..., n; n - число измерений; А_i - активность радионуклида в моче или кале человека в момент времени t_i; R(t) - функция выведения радионуклида с мочой или калом человека, приведенная в приложении 2; П_j - значение величины однократного ингаляционного поступления радионуклида за j-ый период контроля; I_j(t) - скорость ингаляционного поступления радионуклида за j-ый период контроля (промежуток времени между двумя последовательными измерениями [t_j-1 - t_j]), которая имеет следующий вид:

                                                                                      (5.2)

где τ_j - момент времени, в который произошло однократное ингаляционное поступление радионуклида за j-ый период контроля (нулевой период контроля либо начинается формально с нуля, либо с даты начала работы в условиях, при которых необходим ИДК, либо с даты «нулевого» измерения, если его проводят перед началом работы индивидуума).

5.3.2. ОЭД на момент времени t_i, E_i связана с поступлением П_i соотношением:

                                                                        (5.3)

Значения коэффициентов ε_R приведены в приложении 2.

5.4. Оценки величин ингаляционного поступления (П) и ОЭД (Е) проводят методом Монте-Карло. При этом, в каждом испытании для каждого периода контроля, j, следующим параметрам присваиваются случайные значения:

A_j - случайное, нормально-распределенное значение со средним, равным измеренному значению активности A_mi и дисперсией равной (ΔA_mi/2)²;

f_j - случайная переменная, имеющая логнормальное распределение с параметрами: x_g и σ_j, где x_g - геометрическое среднее и σ_j - геометрическое стандартное отклонение [5, 6];

τ_j - дата поступления является случайной величиной, равномерно распределенной в пределах (t_j-1, t_j).

5.5. Для каждого статистического испытания (с номером N) и присвоенных случайных значений A_jN и τ_jN в соответствии с заданной в п. 5.3 моделью рассчитываются значения величин П_jN; ∑П_jN и ОЭД (E_jN; ∑E_jN):

                                                              (5.4)

  (5.5)

                                                                                      (5.6)

                                                   (5.7)

где: N - номер статистического испытания.

5.6. Используя значения поступлений и ОЭД за периоды контроля (П_jN; ∑П_jN и E_jN; ∑E_jN) рассчитывают годовые и соответствующие суммарные поступления и ОЭД, используя правило: годовое поступление (ОЭД) равно сумме полных годовых поступлений (ОЭД) за периоды контроля, полностью принадлежащие данному календарному году, и пропорциональных долей годовых поступлений (ОЭД) тех периодов контроля, которые принадлежат данному и последующему (предыдущему) календарному году.

5.7. В качестве результирующих значений поступлений и доз принимают средние (медианные) значения полученных распределений:

∑П_im(∑П_i50) = <∑П_iN>;    ∑E_im(∑E_i50) = <∑E_iN> ,                                  (5.8)

где i = 1, 2,... - порядковый номер календарного года

При этом,

если ∑П_i < ∑П_i-1, ∑П_i-2,..., ∑П_i-k, то величинам ∑П_i < ∑П_i-1, ∑П_i-2,..., ∑П_i-k присваивается «наилучшее» значение (∑П_i + ∑П_i-1 + ∑П_i-2 + ... ∑П_i-k)/(k + 1);

если ∑E_i < ∑E_i-1, ∑E_i-2,..., ∑E_i-k, то величинам ∑E_i < ∑E_i-1, ∑E_i-2,..., ∑E_i-k присваивается «наилучшее» значение (∑E_i + ∑E_i-1 + ∑E_i-2 + ... ∑E_i-k)/(k + 1).

5.8. Находят «наилучшие» средние (медианные) значения годовых поступлений и ОЭД П_im; E_im и П_i60; Е_i50 как разницу соответствующих «наилучших» суммарных значений ∑П_im(∑П_i50) и ∑E_im(∑E_i50).

5.9. В качестве погрешности расчета принимают разницу между верхней границей интервала (П₉₅, Е₉₅), в который попадают 95 % результатов П_iN и E_iN, полученных методом Монте-Карло, и средними значениями П_i и E_i.

5.10. Пример расчета поступления радионуклида и ОЭД с использованием приведенного выше алгоритма представлен в Приложении 3.

5.11. Используемые дозовые коэффициенты перехода от поступления радионуклидов к ОЭД ε_R и функции удержания радионуклидов в теле человека R(t) приведены в Приложении 2.

6. Оформление результата расчета индивидуальной дозы

4.8. Результат расчета оформляют протоколом в форме индивидуальной записи об облучении работника в течение всего контролируемого периода.

4.9. В записи об облучении работника отражаются:

• идентификационная информация об индивидууме и его профессиональной деятельности;

• индивидуальная ОЭД, полученная в течение текущего и всех предшествующих периодов контроля и календарных лет с указанием неопределенности ее расчета.

4.10. При оформлении протокола учитывают требования статьи 11 Федерального Закона «Об информации, информатизации и защите информации» № 24-ФЗ от 20.02.1995 г. о том, что индивидуальные записи об облучении работника относятся к категории конфиденциальной информации.

7. Библиографические данные

1. Нормы радиационной безопасности (НРБ-99). Гигиенические нормативы СП-2.6.1.758-99. М.: Центр санитарно-эпидемиологического нормирования, гигиенической сертификации и экспертизы Минздрава России, 1999.

2. Основные санитарные правила обеспечения радиационной безопасности (ОСПОРБ-99). Санитарные правила СП-2.6.1.799-99. М.: Центр санитарно-эпидемиологического нормирования, гигиенической сертификации и экспертизы Минздрава России, 2000.

3. МУ 2.6.1.16-2000. Определение индивидуальных эффективных и эквивалентных доз и организация контроля профессионального облучения в контролируемых условиях обращения с источниками излучения. Общие требования.

4. МУ 2.6.1.026-2000. Дозиметрический контроль профессионального внутреннего облучения. Общие требования.

5. Хохряков В.Ф. Обмен промышленных соединений плутония в организме человека. Моделирование транспорта и разработка косвенной дозиметрии. Дис. на соиск. уч. степ, докт. биол. наук. М., 1984, 322 с.

6. Бадьин В.И., Булдаков Л.А., Гастева Г.Н., Молоканов А.А., Мордашева В.В., Сорокин А.В. Исследование клинических эффектов у работников при хроническом поступлении плутония: 30-летнее наблюдение за критической группой лиц. IRPA Regional Congress on Radiation Protection in Central Europe, August 22 - 27, 1999, Budapest, Hungary.

7. L. Bertelly, A. Puerta, M.E. Wrenn and J.L. Lipsztein, Bioassay interpretation and dosivetry using specific absorption parameters for UO2 and U3O8. Radiation Protection Dosimetry Vol. 79, Nos 1 - 4, hh. 111 - 113 (1998) Nuclear Technology Publishing.

8. E. Ansoborlo, R.A. Guimette, M.D. Hoover et al, Application of in vitro dissolution tests to different Uranium compounds and comparison with in vivo data. Radiation Protection Dosimetry Vol. 79, Nos 1 - 4, hh. 33 - 37 (1998) Nuclear Technology Publishing

9. ICRP, Age-Dependent Doses to Members of the Public from Intake of Radionuclides: Part 1, ICRP Publication No. 56, Ann. ICRP 20 (2), Pergamon Press, Oxford (1989).

10. ICRP, Human Respiratory Tract Model for Radiological Protection, ICRP Publication No. 66, Ann. ICRP 24 (1 - 3), Elsevier Science, Oxford (1994).

11. ICRP, Age-Dependent Doses to Members of the Public from Intake of Radionuclides: Part 2, Ingestion Dose Coefficients, ICRP Publication No. 67, Ann. ICRP 23 (3/4), Elsevier Science, Oxford (1993).

12. ICRP, Age-Dependent Doses to Members of the Public from Intake of Radionuclides: Part 3, Ingestion Dose Coefficients, ICRP Publication No. 69, Ann. ICRP 25 (1), Elsevier Science, Oxford (1995).

13. ICRP, Age-Dependent Doses to Members of the Public from Intake of Radionuclides, Part 4, Inhalation Dose Coefficients, ICRP Publication No. 71, Ann. ICRP 25 (3 -4), Elsevier Science, Oxford (1995).

14. N. Ishigure. Database of Calculated Values of Retention and Excretion for Members of the Public Following Acute Intake of Radionuclides. Internal Dosimetry of Radionuclides. Occupational, Public and Medical Exposure. Proceedings of a Workshop, Oxford, United Kingdom, September 9 - 12 2002.

8. ПРИЛОЖЕНИЯ

8.1 Приложение 1.

Работа с программой ММК-01

1. Компьютерная программа ММК-01 предназначена для проведения расчетов индивидуальных доз персонала в соответствии с МБР. Она состоит из трех функциональных блоков:

• Блок формирования электронной базы данных результатов систематических измерений активности радионуклидов в биопробах человека;

• Блок проведения расчетов индивидуальных ожидаемых эффективных доз облучения работников за любой календарный год по результатам измерений, содержащихся в электронной базе данных;

• Блок выдачи информации по уровню внутреннего облучения работника в виде протокола и графических иллюстраций.

2. Блок формирования электронной базы данных осуществляет автоматическую запись необходимых исходных данных и результатов измерений в виде таблицы в формате Access и имеет форму:

Операции с электронной базой данных ИДК осуществляются при помощи:

• кнопок, расположенных в верхней части формы;

• алфавита, отображенного в левом окне с помощью «дерева» (если в базе данных имеется объект (человек) с фамилией, начинающейся с определенной буквы, значок [+] информирует об этом);

• вызова соответствующих таблиц (Персональные данные, Информация об измерениях).

База данных ИДК первоначально содержит тестовые примеры на все радионуклиды, на которые программа позволяет проводить расчеты индивидуальных ожидаемых эффективных доз облучения работников (о чем программа информирует значком [+] в левом окне).

3. Для вызова данных и проведения расчетов необходимо выбрать объект (человека), нажав указателем мыши на соответствующий [+] и соответствующий объект (фамилию). После этого нажать указателем мыши на вкладку «Информация об измерениях» - появится

таблица, отображающая даты измерений и тип измерения. Далее поместить указатель мыши на соответствующую дату (тип) измерения и нажать на правую кнопку - появится меню, позволяющее выбрать 4 операции (см. рисунок):

При выборе меню «Показать результаты текущего типа измерения» появится форма, отображающая указанные результаты измерений:

Она позволяет проводить следующие процедуры:

• Ввод данных для расчета.

• Редактирование (удаление) данных.

При нажатии кнопки «Ввести данные для расчета» появляется форма, содержащая таблицу с необходимыми для расчета доз исходными данными:

Программа позволяет проводить расчет «Дозы» и «Поступления» (для изотопов урана поступление рассчитывается в миллиграммах). При нажатии кнопки «Дозы» (или «Поступления») на форме появляется таблица с результатами расчета дозы (поступления):

В таблице представлены результаты расчета средних и медианных значений суммарной (E_s) и годовой (Е) ОЭД, верхние границы интервалов (Е₉₅), в которые попадают 95 % результатов расчетов E_iN, полученных методом Монте-Карло, а также относительная неопределенность (U/E) оценки значений величин Е₅₀. Воспроизводимость результатов (статистических) расчетов зависит от числа испытаний. Для изменения погрешности расчета (E_s) и (Е) предусмотрена возможность изменения числа испытаний (кнопка «Изменить число испытаний»), которая появляется, если нажать кнопку «Возврат», расположенную под таблицей результатов расчета.

При нажатии кнопки «Графическая иллюстрация результатов расчета» появляется форма, которая позволяет увидеть следующие результаты расчета в графическом изображении:

• среднее («наилучшее») значение ОЭД за каждый календарный год - годовая Е₅₀ (мЗв/год);

• среднее («наилучшее») значение суммарной ОЭД на окончание каждого календарного года - E_s50 (мЗв);

• верхнюю границу 95 %-го интервала частотного распределения годовых (Е₉₅) и суммарных (E_s95) ОЭД (мЗв).

При нажатии кнопки «Показать все возможные варианты расчета» появляется форма, которая позволяет дополнительно увидеть все возможные варианты расчета годовых ОЭД в графическом изображении:

4. Блок выдачи информации позволяет выводить результаты расчетов в виде протокола, в котором указывают:

• идентификационные данные работника;

• суммарную ОЭД внутреннего облучения работника (обусловленную ингаляционным поступлением данного радионуклида в организм) на текущий момент времени (а также на момент окончания любого предыдущего календарного года) и неопределенность ее расчета;

• годовые ОЭД внутреннего облучения работника (обусловленные ингаляционным поступлением данного радионуклида в организм) за последние 5 лет (или за все предшествующие годы работы).

8.2. Приложение 2

Дозовые коэффициенты ε_R и функции выведения R(t), используемые для расчета поступления радионуклидов и ОЭД

*⁾ Рекомендуемые в данной методике значения дозовых коэффициентов для соединений урана UO₂, U₃O₈ [7, 8].

Тритий

Функция выведения Н-3 с мочой при ингаляции паров тритированной воды (концентрация Н-3 в моче после однократного поступления 1 Бк трития в организм при ингаляции), Бк/л.

Значения параметров случайной переменной f: x_g = 1; σ_g = 1,5.

Железо

Функции выведения Fe-59 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 1,5.

Функции выведения Fe-59 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Кобальт

Функции выведения Со-57, 58, 60 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 1,5.

Функции выведения Со-57, 58, 60 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Стронций

Функции выведения Sr-85, 89, 90 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 1,5.

Функции выведения Sr-85, 89, 90 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Рутений

Функции выведения Ru-106 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 1,5.

Функции выведения Ru-106 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Йод

Функции выведения I-125, 129, 131 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 1,5.

Функции выведения I-125, 129, 131 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Цезий

Функции выведения Cs-134, 137 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 1,5.

Функции выведения Cs-134, 137 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Радий

Функции выведения Ra-226, 228 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Функции выведения Ra-226, 228 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 5,0.

Торий

Функции выведения Th-228, 232 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Функции выведения Th-228, 232 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 5,0.

Уран

Функции выведения U-234, 235, 238 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Функции выведения U-234, 235, 238 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 5,0.

Плутоний

Функции выведения Pu-238, 239, 240 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Функции выведения Pu-238, 239, 240 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 5,0.

Америций

Функции выведения Am-241 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Функции выведения Am-241 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 5,0.

Нептуний, кюрий и калифорний

Функции выведения Np-237, Cm-242, 244 и Cf-252 с мочой при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 2,0.

Функции выведения Np-237, Cm-242, 244 и Cf-252 с калом при ингаляции в стандартных условиях

Значения параметров случайной переменной f: x_g = 1; σ_g = 5,0.

8.3. Приложение 3

Пример расчета индивидуального поступления радионуклида и ОЭД

Рассмотрим процедуру расчета индивидуального поступления радионуклида и ОЭД на примере ингаляционного поступления плутония в организм персонала в течение 5 лет от начала работы (начала потенциального поступления).

Для простоты расчета (вручную) предположим, что поступление плутония в виде малорастворимых соединений (тип «М») аэрозоля с АМАД = 1 мкм происходило однократными порциями в середине каждого года в размере 1333 Бк каждая. При этом измерения плутония в моче проводили в конце каждого года (не обязательно точно), например, на 315, 730, 1095, 1460 и 1825-е сутки от начала работы. При таком поступлении активность плутония в суточном количестве мочи в соответствии с моделью (для среднего человека) должна иметь значения:

Однако в действительности, даже если человек, полностью соответствует средним параметрам (выведения), при определении радионуклида в моче имеются две составляющие неопределенности:

1. Статистическая погрешность определения активности радионуклида в измеряемом образце.

2. Неопределенность, связанная с биологической вариабельностью выведения радионуклида с мочой.

Результатом этих неопределенностей является то, что в приведенном примере ингаляционного поступления плутония реальные результаты измерений плутония в моче человека (соответствующего средним параметрам выведения) будут заметно отличаться от точных значений, например, они будут такими как показано на рисунке 1:

Рисунок 1.

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

1. А1 (t = 315) = (1,5 + 0,3) мБк/сутки

2. A2(t = 730) = (0,47 ± 0,1) мБк/сутки

3. A3(t = 1095) = (0,81 + 0,2) мБк/сутки

4. A4(t = 1460) = (3,9 ± 0,6) мБк/сутки

5. A5(t = 1825) = (3,3 ± 0,5) мБк/сутки

Расчет значений поступлений радионуклида за год и соответствующих ОЭД по приведенной выше стандартной модели состоит в следующем:

1. В каждом периоде контроля задают дату (число в сутках, t, от начала работы) однократного поступления радионуклида (которая является случайной величиной в соответствии с п. 5.4). В таблице 1 представлены интервалы времени между датой поступления (случайной величиной) и датой измерения, соответствующего j-ому периоду контроля (T1-t, T2-t и так далее).

Таблица 1.

2. Соответственно интервалам (T1-t, T2-t и т. д.) вычисляют значения функции выведения R(T1-t) и т. д.

Таблица 2.

3. Рассчитывают случайные значения измеренной активности А_j и коэффициента f_j (в соответствии с п. 5.4) для каждого периода контроля (j = 1, 2,..,5) и соответствующие поступления и ОЭД (за период контроля):

Таблица 3.

4. Рассчитывают годовые и соответствующие суммарные поступления и ОЭД по формулам (в соответствии с п. 5.6):

Таблица 4.

5. Повторяют процедуру расчета по пп. 1 - 4 многократно (в рассмотренном примере N = 500), получая массивы значений {П_iN}; {E_iN} и {∑П_iN }; {∑E_iN }. На рисунке 2 для иллюстрации показан массив значений {E_iN - серые линии}, полученный с помощью программы ММК-01:

Рисунок 2.

6. Находят средние значения годовых и соответствующих суммарных поступлений и ОЭД: П_i = <П_iN> и Е_i = <E_iN>; ∑П_i = <∑П_iN> и ∑Е_i = <∑E_iN>.

7. Находят медианные и верхние границы интервала, в который попадают 95 % результатов значений {П_iN}; {E_iN} и {∑П_iN}; {∑E_iN} годовых и соответствующих суммарных поступлений и ОЭД: П_i50; П_i95; ∑П_i50; ∑П_i95 и E_i50; E_i95; ∑E_i50; ∑E_i95. Значения E_i95 и ∑E_i95 также представлены на рис. 2 в таблице и в графическом виде (E_i95 - линия 1; ∑E_i95 - линия 2).

8. Находят «наилучшие» средние и медианные значения ∑П_im; ∑Е_im и ∑П_i50; ∑E_i50 в соответствии с п. 5.5:

∑Е_2m(= 24,3 мЗв) < ∑Е_2m(= 83,7 мЗв) → ∑Е_2m = (24,3 + 83,7)/2 = 54,0 мЗв

∑Е_3m(= 50,7 мЗв) < ∑Е_2m (= 54,0 мЗв) → ∑Е_3m = (24,3 + 83,7 + 50,7)/3 = 52,9 мЗв

∑Е_1m = ∑Е_2m = ∑Е_3m = 52,9мЗв

∑Е_4m(= 233,7 мЗв) > ∑Е_3m(= 52,9 мЗв) → ∑Е_4m = 233,7 мЗв

∑Е_5m (=194,4 мЗв) < ∑Е_4m (=233,7 мЗв) → ∑Е_5m = (194,4 +233,7)/2 = 214,05 мЗв

∑Е_4m = ∑Е_5m = 214,05 мЗв

∑Е₂₅₀ (= 19,7 мЗв) < ∑Е₂₅₀ 260(= 69,6 мЗв) → ∑Е₂₅₀ = (19,7 + 69,6)/2 = 44,7 мЗв

∑Е₃₅₀ (= 42,7 мЗв) < ∑Е₂₅₀(= 44,7 мЗв) → ∑Е₃₅₀ = (19,7 + 69,6 + 42,7)/3 = 44,0 мЗв

∑Е₁₅₀ = ∑Е₂₅₀ = ∑Е₃₅₀ = 44,0 мЗв

∑Е₄₅₀(= 182,4 мЗв) > ∑Е₃₅₀(= 44,0 мЗв) → ∑Е₄₅₀ = 182,4 мЗв

∑Е₅₅₀(=151,6 мЗв) < ∑Е₄₅₀(=182,4 мЗв) →∑Е₅₅₀ = (151,6 +182,4)/2 = 167,0 мЗв

∑Е₄₅₀ = ∑Е₅₅₀ = 167,0 мЗв

«Наилучшие» значения ∑E_im и ∑E_i50 представлены выше на рис. 2 в таблице, а также в графическом виде (∑E_i50 - линия 3).

9. Находят «наилучшие» средние и медианные значения годовых поступлений и ОЭД П_im; Е_im и П_i50; E_i50 как разницу соответствующих «наилучших» суммарных значений, которые представлены также на рис. 2 в таблице и в графическом виде (E_i50 - линия 4).

8.4. Приложение 4. (Обязательное) Список исполнителей

Исполнители:

Зав. кабинетом БФИ ГНЦ РФ «Институт биофизики»,

ведущий науч. сотр., к.т.н. …………………………………………………… А.А. Молоканов

Зав. отделом Федерального управления «Медбиоэкстрем», к.м.н. …………... Е.Б. Антипин

Зав. лабораторией БФАИИЧ ГЦ Госсанэпиднадзора

ФУ «Медбиоэкстрем», д.т.н. ……………………………………………………… В.И. Бадьин

Зав. лабораторией ГНЦ «Институт биофизики», к.ф.-м.н. …………………… В.П. Крючков

Зав. лабораторией ГНЦ «Институт биофизики». ……………………………... А.Г. Цовьянов

Вед. науч. сотр. ЦМИИ ГП «ВНИИФТРИ» Гостстандарта РФ, к.т.н. ………. А.П. Ермилов

Зав. лаб. ЦМИИ ГП «ВНИИФТРИ» Гостстандарта РФ, к.т.н. …………… С.Ю. Антропов