Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Климанов Дозиметрия ионизируюшчикх излучениы 2015

.pdf
Скачиваний:
64
Добавлен:
12.11.2022
Размер:
61.09 Mб
Скачать

 

q(t) =

 

I 0 f

[1 − exp(−λэфф t)] ;

(17.18)

 

 

λэфф

 

 

 

 

 

 

 

 

 

(при t

= ∞

q(t) =

I0 f

).

 

 

 

 

 

 

 

 

 

 

λэфф

 

Мощность дозы D(t) определяется соотношением:

 

D(t) =

I0 f

[1− exp(−λэфф t)] 1,6 1013 Eэфф ,Зв/с

(17.19)

λэфф m

 

 

 

 

 

 

 

 

С помощью соотношения (17.19) определяем предельнодопустимую концентрацию Av , Бк/м3, некоторого нуклида в возду-

хе, поступающего в дыхательную систему (и далее в щитовидную железу) со скоростью V, м3/с; ( I0 = Av V ).

Параметры задачи:

m = 0,02 кг – масса щитовидной железы; Eэфф = 0,23 МэВ/распад;

f = 0,23; λэфф = 1,06·10-6 1/с (Tэфф = 7,6 суток);

 

 

V = 7300 м3/год (годовое потребление воздуха);

 

 

D

 

= 0,020 Зв/год, (предельно допустимая мощность дозы).

 

доп

 

 

 

 

 

 

 

 

I0

=

 

P λ

эфф m

=

0, 02 1, 06 106

0, 02

=

f [1 − exp(−λэфф t)] 1, 6 1013 Eэфф

0, 23 0, 23 1, 6 1013

 

 

 

 

= 4,0·10-4 Бк/год.

 

 

 

 

 

 

 

 

 

A = I

0

/V = 4,0·10-4/7300 = 5,5 Бк/м3.

 

 

 

 

v

 

 

 

 

 

 

6.3. Двухкамерная модель

Диаграмма переходов нуклидов показана на рис. 17.8.

611

 

 

 

камера 1

 

камера 2

 

 

 

 

 

 

 

 

а2к1

q0,

 

q1, m1, c1, k1

 

 

q2, m2, c2, k2

(t

= 0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b1k2

 

 

 

 

aek1

 

 

 

 

 

 

 

 

 

afk1

 

 

b3k2

 

 

 

 

 

 

 

 

 

 

 

 

 

выведение

 

 

выведение

фиксация

 

 

 

 

Рис. 17.8.

Схема двухкамерной модели

Параметры и переменные варианта модели на рис. 17.8: к1 и к2 – скорость перехода активности из камер 1 и 2;

а2к1 доля активности, переходящей в единицу времени из камеры 1 в камеру 2;

b1к2 доля активности, переходящей в единицу времени из камеры 2 в камеру 1;

b3к2 доля активности, покидающая в единицу времени всю систему;

а2+ ае+ аf =1; b1+ b3=1.

Система дифференциальных уравнений для текущих значений активностей q1 и q2 в двух камерах имеет следующий вид:

dq1

= b k q

( k

 

+ λ

p

)q

;

 

 

 

 

 

 

dq2

1 2 2

1

 

1

 

(17.20)

= a k q

( k

2

+ λ

p

)q

2

.

 

 

dt

2 1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответствующие концентрации нуклидов в начальный момент времени t= 0:

q (0) = q

0

,

c (0) =

q0

= c

0

,

q

2

(0) = 0,

c

2

(0) = 0 ;

 

1

 

1

m1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зависимости концентраций в камерах 1 и 2 от времени имеют

вид:

 

 

 

 

 

 

 

 

 

 

k2

λб

 

λбt

 

k1 λб

 

( λб

+γ)t

 

c1(t) = c0exp(λpt)

 

 

e

 

+

 

e

 

 

;

γ

 

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

 

612

 

 

 

 

 

 

c

(t) = c

 

m1

a2k1 exp(λ

p

t) eλбt

eб+γ)t

;

(17.21)

 

2

0

 

m2

γ

 

 

 

 

 

 

 

 

 

 

 

 

где γ = (k

k

2

)2 + 4a b k k

1/ 2

, λ

б

= 0,5 (k + k

2

γ) .

1

 

2 1 1 2

 

 

1

 

7. Кинетика продуктов распада радона на фильтре

7.1. Кинетика продуктов распада радона на фильтре в процессе прокачки

Как отмечалось выше, воздух содержит радон с его продуктами распада. Основным способом определения концентраций радиоактивных аэрозолей в воздухе является прокачка воздуха через фильтры с последующим измерением активности, осажденной на фильтрах. Пусть смесь воздуха с аэрозолями прокачивается через фильтр со скоростью w; эффективность фильтра предполагается равной единице:η = 1. Накопление дочерних продуктов распада учитывается в соответствии со следующей частью цепочки распада

222Rn (см. также рис. 17.2):

 

 

 

α

β

β,α

 

RaA (218Po) → RaB (214Pb) → RaC (214Po)→ RaD (214Bi)→

~ 3,1 мин.

~ 26,8 мин.

164 мкс

~ 19,9 мин.

Обычно целью процедуры прокачки и последующего радиометрического определения значений активностей на фильтре является оценка удельной концентрации дочерних продуктов распада 222Rn в

воздухе. Если Nа, Nb ,Nc число атомов RaA, RaB, RaC (в виде

аэрозолей), осевших на фильтре за время t, и na, nb ,nc концентрации соответствующих атомов в воздухе, то система дифференциальных уравнений, определяющих кинетику атомов на фильтре, имеет вид:

613

 

 

 

 

 

 

 

dNa

= n

a

w λ

a

N

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(17.22)

 

 

 

 

 

 

 

dNb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= nb

w + λa N a λb Nb

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

= nc w + λb Nb λc N ,

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

где λa, λb, λc постоянные распада трех нуклидов.

 

 

При начальных условиях: t = 0, Nа, Nb, Nc

= 0 решение системы

(17.22) следующее:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

a

=

 

 

na w

 

F

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

a

 

 

 

 

 

 

 

 

a1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

b

=

 

na w

 

F

 

 

 

+

nba w

F

 

,

 

 

 

 

(17.23)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

b

 

 

 

 

 

 

 

b1

 

 

 

 

λ

b

 

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

c

=

na w

 

F

 

 

 

+

nb w

F

+

nc w

F

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

c

 

 

 

 

 

 

 

c1

 

 

 

 

λ

c

 

c2

 

 

 

λ

c

c3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Fa1 , Fb1 , Fb2 ,

Fc1 , Fc2 , Fc3 функции времени t и постоянных

распада λ :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fa1 = 1exp(λat),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fb1 = 1

 

λb

 

 

 

exp(λat) +

 

 

 

λa

 

 

 

 

exp( λbt),

 

 

 

 

 

 

 

λ

λ

 

 

 

 

λ

 

 

λ

a

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fb2 = 1exp(λbt),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fc1 = 1

 

 

 

λb λc

 

 

 

 

 

 

exp(λat) +

 

 

 

 

 

 

 

 

 

 

(17.24)

b λa ) (λc λa )

 

 

 

 

 

 

 

 

 

 

+

λa λc

 

 

 

 

 

 

exp(λbt)

 

+

 

 

 

 

 

 

λa λb

 

 

 

exp(λct) ,

 

b λc ) c λb)

 

c λa ) (λc λb)

 

Fc2 = 1

 

λc

 

 

 

exp(λbt) +

 

 

λb

 

 

 

exp(λct),

 

 

 

 

 

 

 

λ

λ

λ

c

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

Fc3 = 1exp(λct).

614

 

3 0

 

 

 

 

 

 

2 5

 

Q C

 

 

 

 

 

 

 

 

 

ед.

2 0

 

Q B

 

 

 

отн.

1 5

 

 

 

 

 

A,B,C,

1 0

 

 

 

 

 

Q

 

 

 

 

 

 

 

5

 

Q A

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

1 0 0

2 0 0

 

 

 

 

 

t, отн.ед.

 

 

 

 

Рис. 17.9. Зависимость величин числа атомов RaA, RaB и RaC на фильтре

 

 

от времени прокачки

 

 

 

 

Соответствующие величины активностей Qa,b,c = λ

N

a,b,c

;

 

 

 

a,b,c

 

 

 

Qa = na w Fa1,

 

(17.25)

 

Qb = na w Fb1 + nb w Fb2 ,

 

 

Qc = na w Fc1 + nb w Fc2 + nc w Fc3

 

 

 

Зависимости Qa , Qb

и Qc от времени прокачки показаны на рис.

17.9.

 

 

 

 

 

 

При t →∞ значения функций F 1 и

 

 

 

 

Qa (∞) = na w

,

(17.26)

Qb (∞) = na w + nb w

Qc (∞) = na w + nb w + nc w.

 

,

 

С учетом распада RaC полная активность Q = Qa + Qb + Qc и

Q() = (4na + 3nb + 2nc ) w

(17.27)

615

 

При равновесии продуктов распада 222Rn в воздухе Qa = Qb = Qc

и

 

 

 

λa na

= λb nb = λc nc , откуда следует:

 

 

 

 

 

 

na

=

λb

;

na

=

λc

;

nb

=

λc

и na : nb : nс = 1:8,8 : 6,45.

 

nb

λa

nc

λa

nc

λb

 

 

 

 

 

 

 

При раздельном режиме измерении α или β активности из

схемы распада следует:

 

 

RaB

β активность,

RaA

α активность;

RaC

 

RaC

 

 

Qα = Qa + Qc ,

Qβ = Qb + Qc и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

Qα = na

(Fa1 + Fc1)+ nbFc2 + ncFc3

 

Qβ = [na (Fb1 + Fc1 )+ nb (Fb2 + Fc2 ) + nc Fc3 ] w .

 

При равновесии дочерних продуктов распада 222Rn

 

 

 

 

 

 

 

 

 

λa na = λb nb = λc nc = A0 ;

 

по измеренным величинам активностей Qαи Qβоднозначно опре-

деляется значение A0 и, соответственно, na, nb ,nc:

 

 

F

 

+

F

 

 

F

 

 

 

F

 

 

 

 

 

 

 

 

Q

=

 

a1

 

 

 

c1

+

 

c2

 

+

c3

 

w

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

λ

a

 

 

 

 

 

 

 

λ

b

 

 

λ

c

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

+ F

 

F

 

+ F

 

 

 

F

 

 

 

 

Q

=

 

b1

 

 

 

c1

+

 

 

c2

 

 

b2

 

+

c3

w

A .

(17.28)

 

λ

 

 

 

 

 

 

 

 

 

 

 

β

 

a

 

 

 

 

 

 

 

 

λ

b

 

 

 

λ

c

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из соотношений (17.28) может быть рассчитана зависимость

 

 

 

 

 

 

 

 

 

Qβ

 

= f (t) ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значение этого отношения при t →∞ и при условии равновесия:

Qβ

=

2na + 2nb + nc

=1,51 .

(17.29)

Q

 

 

2n

+ n

+ n

 

 

α

 

a

b

c

 

 

В крайнем случае неравновесия (присутствует только RaA)

616

 

Qβ

= 1.

 

 

 

Q

 

 

 

 

α

 

 

 

 

1 ,6

 

 

 

α

1 ,4

 

 

 

/Q

 

 

 

 

 

 

 

β

 

 

 

 

Q

 

 

 

 

 

1 ,2

 

 

 

 

1 ,0

 

 

 

 

0

1 0 0

 

2 0 0

 

 

t, мин

 

 

 

Рис. 17.10. Зависимость отношения

Qβ

от времени

 

 

 

Qα

 

7.2. Кинетика продуктов распада радона на

 

фильтре после прокачки

После окончания прокачки воздуха изменение количества ранее накопленных радионуклидов соответствует системе уравнений:

dN

 

 

 

 

 

 

 

 

 

 

 

 

a

 

= −λ

 

 

N

;

 

 

 

 

 

dt

a

 

 

 

 

 

 

 

 

a

 

 

 

 

(17.30)

dN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

= λ

 

N

λa N

;

 

dt

b

 

 

 

 

b

 

 

 

 

a

 

 

dN

 

 

 

 

 

 

 

 

 

 

 

 

c

= −λ

c

N

+ λ

b

N .

 

dt

 

 

 

 

 

c

 

 

b

Начальные условия при t′ = 0 :

Na(0) = Na (t), Nb(0) = Nb (t) Nc(0) = Nc (t) ,

где t – время прокачки. Решение системы (17.35) зависит от пара-

метров, λa , λb , λc , t , Nа(t), Nb(t), Nc(t): 617

N

 

N (t), λ

 

,t

;

 

 

 

 

 

 

 

 

 

 

 

 

a

(t ) = f

a

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

N

 

(t), N (t), λ

 

, λ

 

,t

;

 

 

 

 

 

(t ) = f

 

a

a

b

 

 

 

 

 

 

b

 

b

 

 

b

 

 

 

 

 

 

 

 

 

 

N

 

N

 

(t), N (t), N (t), λ

 

 

, λ , λ

 

,t

.

c

(t ) = f

c

a

a

c

 

 

 

 

 

b

 

 

c

 

 

 

 

b

 

 

Величины α- и β-активности через время выдержки t соответ-

ственно равны:

 

 

 

 

;

 

Qα

= λa Na

+ λc Nc

 

 

(17.31)

Qβ

= λa Na + λc Nc .

7.3. Фильтры для измерения концентрации аэрозолей

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

Аэрозольные частицы оседают в фильтре в результате:

столкновений с волокнами фильтра;

отклонения от первоначального направления движения, что увеличивает вероятность их захвата;

электростатического взаимодействия между фильтром и частицами.

При прокачке через фильтр толщиной х воздуха с концентрцией

аэрозолей А0, за фильтром их концентрация снизится до величины А(х); если аэрозоли одинакового размера, плотность фильтра однородна и скорость прокачки воздуха постоянна, то имеет место соотношение:

618

A(x) = A0 exp(k x) ,

(17.32)

где k – коэффициент фильтрации, см-1, зависящий от ряда факторов, основные из которых следующие:

свойства фильтра;

форма, размеры и плотность аэрозольных частиц;

агрегатное состояние частиц;

концентрация частиц в воздухе;

скорость прокачки воздуха.

Значение эффективности фильтра η определяется как

η =

A0 A(x)

= 1 exp(k x) ,

(17.33)

 

 

A0

 

а величина ε = 1-η – соответствующий коэффициент проскока. Определение эффективности фильтра может быть реализовано

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

h h

A0

A1

 

A2

 

 

 

 

Рис. 17.11. К определению эффективности фильтра

Для такой геометрии измерений одинаковые эффективности фиьтров η могут быть определены как

η=

A0 Ah1

=

c Q1

, η=

Ah1 Ah2

=

c Q2

,

(17.34)

A

A

A

A

 

 

 

 

 

 

 

0

0

 

h1

0

 

 

где Q1,2 – фиксированные в фильтрах активности, с – константа; из

этих отношений следует:

Ah1

=

Q2

и так как η = 1

Ah1

, то

 

A

Q

A

 

 

 

 

 

 

0

 

1

 

 

0

 

 

 

 

619

 

 

 

η= 1

Q2

.

(17.35)

 

 

Q

 

 

1

 

 

Величина проскока εn через n параллельных одинаковых фильтров рассчитывается по соотношению:

εn = ε1 ε2 εn

(17.36)

и, соответственно, ηn =1 εn .

Контрольные вопросы

1.Основные различия дозиметрии внешних потоков излучения

иинкорпорированных нуклидов.

2.Задачи дозиметрии инкорпорированных нуклидов.

3.Пути поступления нуклидов в организм человека и их классификация по степени опасности.

4.Приближенные оценки поглощенной энергии инкорпорированных нуклидов.

5.Зависимость интегральной дозы инкорпорированных нуклидов от времени.

6.Основные приближения, положенные в камерные модели.

7.Критерии опасности дозиметрических нагрузок инкорпорированных нуклидов.

8.Экспериментальный критерий равновесного состояния дочерних продуктов распада радона.

Список литературы

1.Иванов В. И. Курс дозиметрии. Учебник, 4е издание. М.: Энергоатомиздат, 1982.

2.Осанов Д.П. и др. Дозиметрия излучений инкорпорированных

радиоактивных веществ, М.: Атомиздат, 1970.

3.Л.С. Рузер. Радиоактивные аэрозоли. М.: Энергоатомиздат,

2001.

4.V.V. Kushin and V.V. Smirnov. Beta Dosimetry of Airborne Hot Particles from Chernobyl Fallout. Rad. Prot. Dos., 74, No.1/2, 1997, pp. 2737.

620

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]