CC BY
36
ЭКОЛОГИЯ
J
УДК 53.072; 504.5
ОПРЕДЕЛЕНИЕ АГРЕГАТНОГО СОСТОЯНИЯ ПРОДУКТОВ ГИДРОЛИЗА ГЕКСАФТОРИДА УРАНА, ОСЕДАЮЩИХ НА КОЖУ ЧЕЛОВЕКА ПРИ АВАРИЙНОМ ВЫБРОСЕ
С.П. Бабенко1, А.В. Бадьин2
Московский государственый технический университет им. Н.Э. Баумана, Москва (e-mail: [email protected]);
2Московский государственный университет им. М.В. Ломоносова, Москва (e-mail: [email protected])
Рассмотрены модели загрязнения производственной среды при аварийном выбросе газообразного гексафторида урана в воздух рабочего помещения и для практического использования отобрана модель, которая лучше соответствует экспериментальным данным. Показано, что доза, получаемая человеком при перкутантном поступлении урана, определяется газовыми компонентами продуктов гидролиза гексафторида урана.
Ключевые слова: гексафторид урана, перкутанное поступление, математическая модель.
DETERMINATION OF STATE OF MATTER OF HYDROLYSIS PRODUCTS OF URANIUM HEXAFLUORIDE FALLING ON HUMAN SKIN WITH EMERGENCY DISCHARGE
S.P. Babenko1, A.V. Bad'in2
1Bauman Moscow State Technical University, Moscow (e-mail: [email protected]);
2Lomonosov Moscow State University, Moscow (e-mail: [email protected])
Two models of the industrial environment contamination with emergency discharge of gaseous uranium hexafluoride into the air of the working room are constructed. The model that best matches the experimental data is chosen for the practical use. It is shown that a dose acquired by a man in percutaneous absorption of uranium is determined by the gas components of the hydrolysis products of uranium hexafluoride.
Keywords: uranium hexafluoride, percutaneous absorption, mathematical model.
Гексафторид урана (UF6, ГФУ) является основным рабочим веществом в технологиях обогащения природного урана изотопом U235. С этими технологиями неразрывно связаны перспективы развития атомной энергетики, которой в настоящее время в России, как и во всех цивилизованных странах, уделяется очень большое внимание.
Известно, что гексафторид урана представляет большую токсическую и радиологическую опасность для человека [1]. В настоящей
работе рассматривается один из вопросов, которые нужно решить для обеспечения возможности расчета дозы, которую получает человек при перкутанном (через кожу) поступлении продуктов гидролиза иЕ6 в организм. Он заключается в определении агрегатного состояния (газы или аэрозоли), в котором продукты гидролиза ГФУ оседают на кожу человека. Для решения этой задачи проводилась верификация модели, построенной для расчета загрязнения производственной среды. Сравнивались рассчитанная поверхностная плотность массы урана р^и, кг/м2, осевшего на кожу человека и производственные поверхности к моменту времени г, и поверхностная плотность массы урана, осевшего на кожу человека и производственные поверхности к моменту времени г, полученная в эксперименте, моделирующем ситуацию аварийного выброса гексафторида урана на предприятии атомной промышленности [4, 6].
При построении модели загрязнения производственной среды использованы следующие положения. Рассматривается некоторая ограниченная область Ц в трехмерном координатном пространстве и некоторый промежуток времени (г0,Т). Вещество иЕ6 (газ) поступает в каждую точку области Ц и сразу же начинается процесс гидролиза, протекающий по схеме [2, 3]
№ + Н2О ^ иОР4 + 2НЕ,
(1)
иОР4 + Н2О ^ ИО2Р2 + 2НЕ.
Вещества ИО2Р2 и ИБ переходят из газообразного состояния в состояние аэрозоля (за счет процесса нуклеации), и за короткое время (Дг « 50 с [4]), в области Ц появляются следующие вещества: № (газ), ИОР4 (газ), Ш2Е2 (газ), ИБ (газ), Ш2Е2 (аэрозоль), ИБ (аэрозоль). Они выводятся из области Ц за счет воздухообмена (предполагается, что частицы каждого вещества исчезают в каждой точке области), диффузионно оседают на границу области Ц, дрейфуют в пространстве (под действием силы тяжести и силы сопротивления среды или вместе с потоком воздуха). Предполагается, что все аэрозольные частицы, содержащие молекулы одного из интересующих веществ, имеют одинаковый радиус, т.е. образуют монодисперсную систему аэрозолей.
Далее рассматриваются отдельно уравнения для урансодержащих и фторсодержащих веществ. Поэтому систему уравнений неразрывности для концентраций молекул веществ удобно записать в общем виде [7]
д N
—Щ = Дик - (Ик, grad щ) + ^ ак>тит + ^ (х, г), (2)
__т=1
к = 1,^, х е Ц, г е (Ьо,Т).
Здесь N — число веществ в составе газов и аэрозолей; пк (х, ¿) — концентрация молекул вещества с номером к в точке х в момент времени ¿; Бк — коэффициент диффузии частиц вещества с номером к (предполагается, что Бк не зависит от х и ¿); щ — скорость дрейфа частиц вещества с номером к (предполагается, что ьк не зависит от X и ¿); {ак,т}к — постоянные, описывающие процессы гидролиза,
нуклеации и воздухообмена; (х, ¿) — плотность мощности внешних источников вещества с номером к в точке х в момент времени ¿.
Для однозначного определения концентраций заданы следующие дополнительные условия, соответствующие рассматриваемым физическим ситуациям:
nk (x,to) = nk,0 (x), к = 1 ,N,x е Q; д
ak (x, t)— nk (x, t) + ßk (x, t) nk (x, t) = Rk (x, t), (3)
dn
к = х е дЦ, г е (¿о,т).
Здесь пк 0 (х) — концентрация молекул вещества с номером к в точке х в момент времени ¿0; П — единичная нормаль к поверхности дЦ, внешняя по отношению к области Ц; ак, вк, Як — некоторые коэффициенты, обсуждаемые далее. Первое из записанных равенств представляет собой начальное условие, второе — краевое условие третьего рода. Если Як = 0, то это условие отражает тот факт, что на границе области Ц частицы вещества с номером к частично поглощаются и частично отражаются. Если Як = 0, то это условие описывает поступление вещества с номером к внутрь области через ее границу. Коэффициенты ак, вк удовлетворяют условиям |ак | + |вк | = 0, аквк > 0 (эти условия связаны с теоремой единственности решения).
Если область Ц не является ограниченной, то в список дополнительных условий нужно включить условие регулярности решения на бесконечности.
Поскольку нет возможности найти явное решение общей системы уравнений (2), (3), переходим к рассмотрению частных случаев, в которых в соответствии с конкретными физическими ситуациями заданы конкретные исходные данные, а именно область Ц, временной промежуток (¿0, Т), коэффициенты диффузии Ок, векторы скорости щ, коэффициенты акт, плотности мощности внешних источников Гк (х, ¿), начальные концентрации пк,0 (х), коэффициенты ак (х, ¿), вк (х, ¿), Як (х, ¿). Для определения агрегатного состояния продуктов гидролиза ГФУ, в котором эти продукты оседают на кожу человека, достаточно рассмотреть два частных случая (модели) приведенных систем уравнений.
Первая модель процессов оседания описывает следующую ситуацию: оседание частиц происходит в полупространстве г > 0 на временном промежутке (0, есть диффузия газов, но нет макроскопического движения газов; есть (медленная) диффузия аэрозолей, нет макроскопического движения аэрозолей; есть гидролиз, есть нуклеа-ция газов, нет воздухообмена.
Полагаем, что
^(х,у,г,г) = ^(г, г) = ^к,оЯ (^ - г),
Пк,о(х, У, г) = Пк,о(г) = Пк,о, Ок(х, у, г) = «к = а,
вк(х,у,г) = вк = в, ^к(х,у,г) = 0
(здесь ^к о > 0, г1 > 0, ик о > 0; а, в — некоторые константы, Н — функция Хевисайда). Тогда систему уравнений неразрывности можно записать в виде
д д 2 N
—ик = Дк—-тик + V" а,к>тит + ^(г, г), дг дг2 ^ ' (4)
т=1 4 7
к = Т^, г е (0, г е (0,
nk|i=0 = nfeio(z), к = 1,N, z G (0, д
-«k^- nk + dz
= 0, к = 1,N, t G (0,
z=0
(5)
3С > 0 36 > 0 Ук = 1, N Vг е (0, V е (0,
(|ик(г,г)|< Се54) .
Применяя к системе (4), (5) преобразование Лапласа и решая получившуюся систему уравнений, находим изображение ик (г,р) функции ик(г, г), используя которое далее находим изображение ]к(г,р) =
дд = Дк д^ик (г,р) функции ;'к (г, г) = Д (г, г) (определяющей
плотность потока молекул вещества с номером к через нижнюю сторону плоскости г = г в момент времени ^ и изображение
1 Г
N1^ (г,р) = Б-^к (г,р) функции N1^ (г, г) = Б / ^т • ;'к (г,т) (опре-
Р Уо
деляющей число молекул вещества с номером к, прошедших через
нижнюю сторону участка плоскости г = г площади Б на временном промежутке [0, г]). Совершая обратное преобразование Лапласа, находим выражения для (0, г), ^ к(0, г).
Принято, что ^к о = 0 при к = , и1>о > 0, ик,о = 0 при к = 2, N (это означает, что рассматривается мгновенный выброс ГФУ в воздух рабочего помещения). При этом вычислялись масса урана
М1>к(0, £) = (0,£), осевшего на участок плоскости 2 = 0 пло-
щадью £ на временном промежутке [0, £] в составе вещества с номером k (здесь ш0,и = 238 • 1,67 • 10-27 кг — масса атома урана), поверхностная плотность массы урана рв,и,й = М1>к(0,£)/Б, осевшего на плоскость 2 = 0 на временном промежутке [0, £] в составе веще-
N
ства с номером к, и масса урана М1(0,Ь) = ^ М1>к(0,£), осевшего
к=1
на участок плоскости г = 0 площадью £ на временном промежутке м в составе всех веществ. При расчете массы урана, осевшего на кожу человека, принималось, что Б = 2 м2 [5]. Результаты расчетов, проведенных по полученным выражениям, представлены в табл. 1-3.
Таблица 1
Динамика оседания массы урана т на кожу человека в составе газов
(а = 0,01, в = 1)
t, с 10 30 50 120 240 360 600 5400 7200
m, мг 1 2,5 3,6 5,5 6,3 6,4 6,4 6,4 6,4
Из данных табл. 1 следует, что согласно первой модели примерно через 10 мин после выброса оседающая масса урана перестает зависеть от времени. Это означает, что газообразные продукты к моменту времени £ « 10 мин практически исчезают из атмосферы рабочего помещения.
Таблица 2
Динамика массы урана т, осевшего на кожу человека в составе всех урансодержащих продуктов гидролиза ГФУ и в составе каждого из них
(а = 0,01, в =1)
t, с m, кг, газы и аэрозоли m, кг, UF6 (газ) m, кг, UOF4 (газ) m, кг, UO2F2 (газ) m, кг, UO2F2 (аэрозоль)
10 1,02 • 10-6 3,03 • 10-7 6,28 • 10-7 8,77 • 10-8 3,44 • 10-12
30 2,52 • 10-6 3,11 • 10-7 1,52 • 10-6 6,88 • 10-7 9,72 • 10-11
50 3,61 • 10-6 3,11 • 10-7 1,90 • 10-6 1,40 • 10-6 3,69 • 10-10
120 5,53 • 10-6 3,11 • 10-7 2,18 • 10-6 3,03 • 10-6 2,49 • 10-9
240 6,31 • 10-6 3,11 • 10-7 2,20 • 10-6 3,79 • 10-6 7,97 • 10-9
360 6,44 • 10-6 3,11 • 10-7 2,20 • 10-6 3,91 • 10-6 1,40 • 10-8
600 6,47 • 10-6 3,11 • 10-7 2,20 • 10-6 3,94 • 10-6 2,62 • 10-8
5400 6,47 • 10-6 3,11 • 10-7 2,20 • 10-6 3,94 • 10-6 2,67 • 10-7
7200 6,47 • 10-6 3,11 • 10-7 2,20 • 10-6 3,94 • 10-6 3,56 • 10-7
Из данных табл. 2 следует, что почти одинаковую долю в поступлении урана имеют и0Р4 (газ) и и02Р2 (газ). На порядок меньшее
поступление обеспечивают иЕ6 (газ) и Ц"02Р2 (аэрозоль). Масса урана, осевшего в составе и02Р2 (аэрозоль), даже через 2 ч после выброса продолжает изменяться.
Из табл. 3 следует, что по экспериментальным данным масса урана, выброшенного в воздух помещения, хорошо совпадает с суммой следующих масс: осевшей к моменту времени г и находящейся в этот момент в воздухе. При этом рассчитанная осевшая масса на два порядка меньше выброшенной и состоит она в основном из массы урана, осевшего в составе газов. Это означает, что, в рамках первой модели загрязнения производственной среды, основная доля урана заключена в аэрозольных частицах, которые оседают очень медленно и слабо загрязняют производственные поверхности и кожу человека, что не соответствует экспериментальным данным. Такой результат позволяет заключить, что при описании движения частиц нужно учитывать их дрейф под действием силы тяжести и силы сопротивления среды.
Таблица 3
Расчетные и экспериментальные значения составляющих массы урана
Масса урана тпов, г, осевшего на производственные поверхности через 1,5 ч после выброса Масса урана Товозд, г, находящегося в воздухе через 1,5 ч после выброса Полная экспериментально зафиксированная масса урана т, г Масса урана т-выбр, г, выброшенного в воздух помещения
Теория (а = 0,01, в =1) Эксперимент
0,26 6,68 19,37 26,05 25,69
Вторая модель процессов оседания описывает следующую ситуацию: оседание частиц происходит в плоском слое 0 < г < Н (здесь Н > 0) на временном промежутке (0, есть диффузия газов, нет
макроскопического движения газов; нет диффузии аэрозолей, есть макроскопическое движение аэрозолей (аэрозольные частицы движутся под лействием силы тяжести и силы сопротивления среды); есть гидролиз и нуклеация газов, нет воздухообмена.
Для системы N веществ в составе газов (величины, относящиеся к этим веществам, имеют индекс к) и одного вещества в составе аэрозолей при задании (х,у, г, г) = 0, ик,о(х,у, г) = ик,о(г) = ик,о, а1,к(х,у,г) = 0, #1,к(х,у,г) = 0, а2,к(х,у,г) = 0, Й2,к(х,у,г) = 0; ^(х,у,г,г) = 0, ио(х, у, г) = 0, а1(х,у,г) = 0, Л1(х,у,г) = 0, а2(х,у,г) = 0, Д2(х,у,г) = 0, учитывая допущения об отсутствии диффузии аэрозоля, систему уравнений неразрывности можно записать в виде
d d2 —П = Dfe — n + dt dz2
N
+ У afe,mnm, k = 1,N, z G (0, h), t G (0, +ro); (6)
m=1
n|t=o = nfc,o(z), k = 1,N, z G (0, h);
|z=0 = 0, (7)
n|z=h = 0, k = T^N, t G (0, +ro);
d _ d N
(9)
—П = -у(г)—п + ^ Ьтпт, - С (0, й), г е (0, +го); (8)
т=1
п^=0 = 0, - е (0,й); <=л = 0, г е (0,
Здесь г — радиус аэрозольных частиц, -у(г) — модуль скорости дрейфа аэрозольных частиц радиуса г, {Ьт}т=т^ — коэффициенты, описывающие процесс нуклеации. Система уравнений для газов решалась методом разделения переменных, для аэрозолей — методом характеристик [7].
При расчетах рассматривался набор радиусов аэрозольных частиц UO2F2 в диапазоне от 5,567 • 10-8 до 1,391 • 10-5 м. Такой выбор определялся данными эксперимента, моделирующего аварийную ситуацию на предприятии атомной промышленности. В модельном эксперименте воздух помещения пропускался через следующие приборы, включенные последовательно: импактор с пятью каскадами, фильтр I, диффузионную батарею, фильтр II [4, 6]. Измерялась активность урана, осажденного на каждом элементе измерительной системы [6]. Знание параметров элементов измерительной системы позволило оценить диапазон радиусов аэрозольных частиц и02Б2 в воздухе помещения, из которого и выбирались радиусы частиц при изучении монодисперсной модели.
На рис. 1 приведены теоретические зависимости поверхностной плотности массы урана р^ и от времени оседания для аэрозолей при нескольких значениях радиуса аэрозольных частиц. Видно, что вначале для любого радиуса г, поверхностная плотность массы урана р^ и возрастает с течением времени. Затем, начиная с момента времени гнас,
величина р^ и, достигая значения 1,181 • 10-3 кг/м2, перестает зависеть от времени. Скорость возрастания р^и уменьшается с уменьшением радиуса частиц г.
Рис. 1. Зависимость поверхностной плотности массы урана рв,и • 104, кг/м2, от времени оседания Ь для аэрозольных частиц разных радиусов (г = 13,9, 5,24, 2,62, 1,05, 0,62 мкм, соответственно кривые 1, 2, 3, 4, 5; рmax = 1,81 • 10-3 кг/м2)
На рис. 2 приведены зависимости поверхностной плотности массы урана рзи,кг/м2, от времени оседания для газов и аэрозолей. Видно, что оседание урана в составе газов по второй модели происходит так же, как и по первой модели. В области насыщения значение рз,и для аэрозолей на два порядка больше, чем для газов. Горизонтальные линии на рис.2 соответствуют максимальному значению рз;и;шах, которого может достигать поверхностная плотность массы урана в соответствии с исходными экспериментальными данными. Верхняя горизонтальная линия соответствует максимальному значению рз;и;шах;а поверхностной плотности массы урана для частиц, оседающих под действием силы тяжести (аэрозолей) только на пол рабочего помещения. Используя исходные данные, нетрудно подсчитать, что рз;и;шах;а = 1,28 х х 10-3 кг/м2. Нижняя горизонтальная линия соответствует максимальному значению рз и шах г поверхностной плотности массы урана для
Рис. 2. Зависимость поверхностной плотности массы урана р$,и, кг/м2 времени £ для разных моделей оседания
от
газов, т.е. для частиц, оседающих за счет диффузии на пол, потолок и стены рабочего помещения. Используя исходные данные, нетрудно подсчитать, что рздша^г = 2,57 • 10-4 кг/м2.
Результаты расчета, проведенного по второй модели, находятся в хорошем качественном и количественном соответствии с максимальными значениями поверхностной плотности массы урана, соответствующими условиям эксперимента. Действительно, поверхностная плотность массы урана уже за время ^20 мин насыщается и для газов, и для аэрозолей. При этом величина р^и для аэрозолей в области насыщения хорошо совпадает с максимальным значением, соответствующим условиям эксперимента. Для урана в составе газов рассчитанное значение р^и в области насыщения много меньше максимально возможного. Это означает, что через несколько минут после выброса почти весь уран сосредоточен в аэрозольной компоненте.
На рис. 3 на кривой теоретической зависимости р^,и от времени оседания отмечены отдельные точки, соответствующие экспериментальным значениям р^и на фильтре Петрянова, на стеклах, на одежде и коже, измеренным через 2 ч после выброса ГФУ. Эти величины определены по непосредственно измеренным значениям числа а-частиц, испускаемых ураном за одну минуту с одного квадратного сантиметра поверхности.
Из сопоставления расчетных значений р^ и для газов и аэрозолей в области насыщения с экспериментальными точками следует, что оседание урана в составе продуктов гидролиза ГФУ на кожу человека определяется газовыми компонентами, что следует из сравнения рассчитанных и измеренных значений р^и- При этом
1) измеренное среднее по отдельным участкам кожи значение рв>и близко к рассчитанному для газов. Это значит, что перкутанное поступление урана в составе продуктов гидролиза ГФУ определяется газовыми компонентами;
2) значение ря,и на одежде определяется как газами, так и аэрозолями, хотя, по-видимому, из-за подвижности человека основная доля аэрозольных частиц быстро осыпается с одежды;
3) из рис. 3 следует, что почти вся осевшая на производственные поверхности масса урана определяется аэрозольными частицами;
4) даже на материале, который хорошо захватывает оседающие частицы, на полу, где оседание максимально, величина р^и остается меньше максимального значения. Это, по-видимому, связано с вторичным сдувом частиц [8] и с тем фактом, что в реальных условиях частицы дрейфуют не строго вертикально и часть частиц, которые согласно расчетам должны были осесть на пол, оседает на стены.
Полученный результат лег в основу метода расчета дозы, которую получает человек при перкутанном поступлении [7] урана.
Рис. 3. Зависимость поверхностной плотности массы урана рзи, кг/м2, от времени оседания для разных моделей
СПИСОК ЛИТЕРАТУРЫ
1. Гастева Г Н., Бадьин В. И., Молоканов А. А., Мордашева В. В. Клиническая токсикология химических соединений урана при хронической экспозиции // Радиационная медицина. В 2 т. Т. 2. - М. : ИздАт, 2001. -С. 369-389.
2. Гидролиз газообразного гексафторида урана в воздухе / В.И. Бадьин и др. // Отчет предприятия п/я М-5122 и п/я. В-2343. - ЦНИИАТОМИНФОРМ, 1975. - 76 с. - ИК № М 34559.
3. Экспрес с-методы измерения степени обогащения гексафторида урана и следовых количеств UF6 в HF атмосфере на основе диодных лазеров ближнего и среднего ИК диапазона // Г.Ю. Григорьев, А.И. Надеждинский, Ш.Ш. Набиев и др.: Препринт ИАЭ 6395/12 - М., 2006.
4. МирхайдаровА. Х.,БабенкоС. П., БадьинВ. И. Газоаэрозольная опасность соединений урана с фтором в аварийных ситуациях // Тез. докл. на междунар. конф. 24-26 апреля 2000 г. - М., СПб.: Гидрометеоиздат, 2000. -С. 229. - Доклад № 352.
5. Человек. Медико-биологические данные. Доклады рабочей группы Комитета II МКРЗ по условному человеку, публикация 60 / Пер. с англ. Ю.Д. Парфенова; Под ред. А.А. Моисеева. - М.: Медицина, 1977. - 496 с.
6. МирхайдаровА. Х. Исследование закономерностей загрязнения воздушной среды и поверхностей в производстве разделения изотопов урана: Дисс.... канд. тех. наук. - М.: 1978. - 217 с. - ИБФ МЗ; Архив ФМБА России.
7. Б а б е н к о С. П. Прогнозирование радиационного и токсического воздействия выбросов гексафторида урана методами математического моделирования: Дисс. ... д-ра техн. наук. - М.: 2008. - 375 с.
8. П а р х о м е н к о Г. М. К вопросу обоснования предельно допустимой загрязненности рабочих поверхностей при работе с а-излучателями // Защитные покрытия в атомной технике / Под ред. С.М. Городинского, В.Б. Тихомирова. -Госатомиздат, 1963. - С. 11-15.
Статья поступила в редакцию 21.02.2012
