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

4

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

Глава 18. Дозиметрия в полях импульсного излучения

Некоторые источники ионизирующего излучения, используемые в различных областях радиационной физики, формируют поля импульсного излучения: источники ионизирующего излучения функционируют в импульсном режиме, т.е. за очень короткий интервал времени t происходит процесс выхода частиц (рис. 18.1). Линейные ускорители, бетатроны, синхротроны генерируют им-

пульсы излучения в интервале времен τ ~ 10-9 ÷ 10-6 с. Импульсы следуют периодически с частотой ν , сек-1, с интервалом времени

Т, сек, между импульсами, ν = T1 . В лучевой терапии применяются

ускорители с импульсами длительностью 3 ÷ 5 μс и частотой 200 ÷ 400 Гц. Существуют также устройства, генерирующие одиночные высокоинтенсивные импульсы электронов и тормозного излучения.

t

T τ

Рис. 18.1. Временные параметры импульсов излучения

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

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

621

1. Импульсный режим облучения ионизационной камеры

1.1. Соотношение Боуга при импульсном облучении

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

10-3 с) и с достаточно малой частотой ν = T1 , где Т время между

импульсами, достаточное для полного сбора ионов на электродах камеры. Основные этапы процессов в ионизационной камере после облучения её коротким импульсом излучения представлены на рис. 18.2; вариант 1 соответствует равномерному облучению газового объёма в момент облучения, после которого начинается процесс движения ионов в электрическом поле на электроды и рекомбинация; второй вариант (2) соответствует времени существования зоны перекрытия шириной Х, где происходит рекомбинация ионов, и двух зон с ионами одного знака; третий вариант (3) – исчезновение зоны перекрытия Х, завершение процесса рекомбинации, и сбора ионов разных знаков на электродах.

d

-++ -

++ - + +-

+-- + --

+-

-1+ -

-

 

Х

-

 

 

 

-

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

++ -

 

 

- -

+ +

-

+

+

 

 

 

- +

-

+

 

 

- -

+

-

+- +

 

 

-

+ +

 

 

2

 

 

 

 

 

3

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

622

Зависимость размера зоны перекрытия Х от времени t при равномерном движении ионов определяется соотношением:

X (t) = d (k

+ k

) U t

,

(18.1)

+

 

d

 

 

где d – межэлектродное расстояние, U – приложенное к камере напряжение, k+ , k– подвижности ионов.

Полное время Т существования зоны перекрытия Х, где происходит рекомбинация ионов после облучающего импульса:

T =

d2

 

[(k+ + k) U ] .

(18.2)

Суммарное количество прорекомбинировавших ионов Qα за время Т существования зоны перекрытия (одномерный случай, площадь электродов S = 1) определяется следующим интегралом:

Q

= T

α n2 (t) X (t)dt ,

(18.3)

α

 

 

 

 

0

 

 

где n(t) – объёмная концентрация ионов в момент времени t (в приближении nn+ = n ); α – коэффициент рекомбинации ионов.

Зависимость (после момента окончания облучения) уменьшения концентрации ионов от времени за счет рекомбинации определяет-

ся уравнением (при t = 0, n = n0 ):

n( t ) =

 

 

n0

 

,

1

+ n α t

с учетом которого

 

 

0

 

 

 

 

 

 

 

Qα = n0

 

d 1

1 ln(1 + η) ,

 

 

 

 

η

 

η= [ α n0 d 2 ] . (k+ + k) U

(18.4)

(18.5)

(18.6)

Так как f = (Q0 Qα ) / Q0 , ( Qα = n0 d ), итоговое соотношение для f (соотношение Боуга):

f =

1

ln(1 + η)

( lim

1

ln(1 + η) → 1).

(18.7)

 

η

 

η0

η

 

 

 

 

 

 

 

623

 

Целесообразно определить величину f относительно измеряемо-

го заряда Q = Q Q ; определим постоянную C =

α d2

 

и

[(k+ + k) U ]

0

α

 

 

 

 

 

 

 

 

 

 

 

 

 

переменную η= C

n

, где

n = Q /V

( n =

n

, Vk – объём камеры).

 

 

 

f

k

0

f

 

 

 

 

 

 

 

 

 

 

 

После подстановок этих величин в соотношение (18.7) это уравнение преобразуется к виду:

f = ln(1 +

C n)

f

 

;

(18. 8)

C n

 

 

 

 

f

 

 

решение (18.8) следующее:

 

 

 

 

n

 

 

 

 

 

f = C

 

 

 

 

.

 

(18.9)

exp(C n) 1

 

Из определения η следует соотношение:

 

 

η=

T

 

=

T

 

 

 

 

 

 

 

 

,

 

 

 

 

(18.10)

1

 

τα

 

 

 

 

 

α n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

ние f зависит от отношений времени существования области перекрытия Х и времени рекомбинации. Для полного восстановления состояния газа время завершения всех процессов τ ≈ 2Т.

Часть р освобожденных в процессе ионизации электронов до момента их прилипания к молекулам газа и образования отрицательных ионов может достичь положительного электрода, что нарушает сделанные выше предположения о характере протекающих в камере процессах. Учет вклада свободных электронов, имеющих большую подвижность, может быть включен в соотношение (18.7) для расчета эффективности собирания ионов следующим образом:

f =

1

ln(1

+ exp( p η) 1) .

(18.11)

 

η

 

p

 

В случае р → 0 уравнение (18.11) переходит в (18.7): 624

lim

1

ln(1

+ exp( p η) 1)

1

ln(1 + η) ).

(18.12)

p 0

η

 

p

η

 

 

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

 

0,9

 

 

 

 

 

p= 0,5

f (η)

0,6

 

p= 0,25

 

 

 

0,3

 

 

 

 

 

p= 0

 

0,1

1

10

η

Рис. 18.3. Зависимости эффективности собирания ионов при различных вкладах свободных электронов

В случае, когда значительное количество импульсов имеют место во время сбора ионов, используются аппроксимирующие условия непрерывного облучения с установлением усредненного по времени значения q q , входящего в соотношение Боуга. Если

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

Типичные ионизационные камеры применимы для измерения доз до максимальных значений ~10 мГр/имп; использование спе-

циальных камер позволяет увеличить верхний предел измерений до 100 мГр/имп. Основная проблема измерения ионизационными камерами больших значений доз в импульсе – учет поправок на эффекты переноса свободных электронов, рекомбинации и влияния объёмного заряда носителей на электрическое поле в камере. Ме-

625

тод определения этих поправок – расчет системы уравнений ионизационной камеры с учетом объёмного заряда ионов.

1.2. Экспериментальный метод определения эффективности собирания ионов в импульсном режиме

Как и в случае непрерывного облучения, целесообразным является экспериментальный способ определения эффективности собирания ионов; в случае импульсного облучения на основании приведенных выше соотношений возможна реализация следующей процедуры: проводятся измерения зарядов Q1 и Q2 для двух одинаковых импульсов излучения при различных напряжениях на камере U1 и U2. Соответствующие значения величин эффективности собирания ионов f1 и f2 определяются соотношениями:

 

 

 

 

f

= Q1

,

 

 

 

 

Q1

=

1

ln(1

+ η ) ;

 

(18.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Q

 

 

 

 

 

 

Q

 

 

η1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f2

=

Q2

,

 

 

 

 

 

Q2

=

1

 

ln(1 + η2 ) ,

 

(18.14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

Q

 

 

 

η2

 

 

 

 

 

 

 

 

где Q – полный заряд, созданный импульсом.

 

 

 

 

 

Отношение

 

 

 

f1

 

 

 

Q1

 

 

 

η2

 

ln(1+ η1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

=

 

 

;

 

 

 

(18.15)

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

f

2

 

 

 

 

 

 

η

 

 

ln(1+

η

2

)

 

 

 

 

 

 

η2

= U1

 

 

 

 

 

U1

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

при

и η

 

= η

 

соотношение (18.15) приводится к ви-

 

η

U

2

 

2

 

1

 

U 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ду:

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q1

 

 

U1

 

 

 

 

 

 

ln(1+η1)

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

,

(18.16)

 

 

 

 

 

 

 

 

 

 

ln(1+η U / U

 

)

 

 

 

 

 

 

 

Q

 

 

 

U

2

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

и отношение измеряемых величин Q1 является функцией η1. Вво-

Q2

дя обозначения C =

Q1

и U =

U1

и используя разложение функции

Q

 

 

U

2

 

 

2

 

 

 

 

 

 

 

 

626

f ( x ) = ln(1 + x) f ( x ) » x - x22 (диапазон значений f > 0,9 –

практически используемый диапазон), уравнение (18.16) преобразуется к виду:

 

 

 

 

 

 

 

η 2

 

 

 

 

 

 

 

η

-

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

C =U

 

 

 

1

 

2

 

,

(18.17)

 

 

 

 

 

 

U

η

-0,5

(U η )2

откуда

 

 

1

 

 

 

1

 

 

 

 

 

 

C -1

 

 

 

 

 

η

=2

 

 

 

.

 

 

(18.18)

 

 

 

 

 

 

 

1

 

 

 

C U -1

 

 

 

 

Известные значения С и U позволяют определить величину эффективности собирания ионов f1 из соотношения (18.13).

1.3. Метод определения эффективности собирания ионов в импульсном режиме облучения

Все вышеизложенное получено в предположении нулевой длительности импульса излучения. Если импульс имеет конечную длительность, необходимо учитывать процессы образования ионов, их рекомбинацию в течение импульса и после его окончания.

С этой целью решается нестационарная задача определения полного заряда ионов, собранных в плоской ионизационной камере, образованных ионизирующим импульсом длительностью τ0 при постоянном приложенном напряжении U0. Пренебрегается эффектами пространственного заряда и диффузией; подвижности положительных и отрицательных ионов предполагаются одинаковыми

kk+ = k . В этом случае система дифференциальных уравнений

в частных производных, моделирующая пространственновременной баланс разноименных ионов, имеет следующий вид (одномерный случай, (см. уравнения 11.17) в нормализованных переменных:

ν1 =1-ξ ν ν

2

-ν1

 

T

¢

1

 

 

l ,

(18.19)

 

 

 

 

 

 

 

 

 

 

 

ν2

=1-ξ ν

ν

 

+

ν2

 

T ¢

1

 

 

2

 

l

 

 

 

627

 

 

 

где q – скорость образования ионов в единице объёма, n1 иn2 – кон-

центрация положительных и отрицательных ионов, τ – время прохождения ионами расстояния d между электродами камеры; без-

размерные параметры: ν1 = qn1τ , ν2 = qn2τ , T ¢ = t/τ, l = x/d. Величина ξ включает основные параметры функционирования камеры:

ξ=α q τ

2

=

α q d

4

(18.20)

 

2 2

,

 

 

 

k U0

 

 

где α – коэффициент рекомбинации разноименных ионов.

В случае стационарного облучения камеры эффективность собирания ионов f определяется соотношением:

f = 1-ξò1

ν1ν2dl ,

(18.21)

0

 

 

где соответствующие значения ν1(T ¢ ,l) и ν2(T ¢ ,l) есть результат решения системы (18.19). Последнее соотношение следует из определения величины эффективности собирания ионов:

d

q α n1n2dx

f =

0

 

.

(18.22)

 

q

 

 

 

 

После подстановки в соотношение (18.22) соответствующих нормализованных переменных и с учетом определения ξ величина эффективности собирания f равна:

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q -ò

α ν1τq ν2

τqdx

 

 

 

 

 

f

=

0

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

(18.23)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=1-αqτ2

1

ν

ν

2

dl

=1-ξ

1

ν

ν

2

dl.

 

 

 

 

ò

1

 

 

 

 

ò

1

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

Величину эффективности собирания ионов в импульсном режи-

ме облучения с учетом длительности импульса T0¢ можно представить в следующем виде:

628

 

 

ξ

1

T ¢

 

1+T ¢

 

 

 

 

 

0

¢

 

0

¢ ¢ ¢

 

 

f

=1-

 

ò (ò ν1ν2dT

+

 

)dl ,

(18.24)

 

 

 

 

ò ν1ν2dT

 

 

T0¢ 0

0

 

 

T ¢

 

 

 

 

 

 

 

 

 

 

0

 

 

 

где ν1 и ν2 ̶значения нормализованных плотностей ионов во

время действия импульса T ¢ ;

ν¢

и

ν¢

̶ соответствующиезначе-

 

1

 

2

 

ния плотностей ионов после окончания импульса; время 1 + T0¢ в

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

Решения системы (18.19) в случае стационарного облучения ионизационной камеры для различных значений мощности дозы отражают отклонение формы пространственного распределения плотности ионов от линейной треугольной зависимости от координаты l (линейная зависимость является одним из основных приближений теории Боуга). На рис. 18.4 показаны распределения ионов при малой мощности дозы (ξ = 10-1) и большой (ξ = 10); первый случай соответствует условиям Боуга, которые нарушаются

при ξ = 10 (ξ ~ q). Нарушение линейности также имеет место в те-

чение действия импульса излучения.

 

1

1,0

 

 

 

 

 

ионов n

 

 

 

плотность

 

ξ = 1 0 − 1

 

0,5

 

 

Нормализованная

 

 

 

ξ = 1 0

 

0,0

 

 

 

 

 

 

0,0

0,5

1,0

 

 

Межэлектродное расстояние в ед. l

 

Рис. 18.4. Пространственное распределение плотности ионов для различных величин ξ

629

Система уравнений (18.19) учитывает также эффекты увеличения плотности ионов и рекомбинацию в течение действия импульса. Эти факторы, не учтенные при выводе соотношения Боуга, обусловливают различия в рассчитанных значениях эффективности собирания ионов f в непрерывном и импульсном режимах облучения ионизационных камер. На рис. 18.5 приведены результаты рас-

четов зависимости величины f от длительности импульса T0¢ (нормализованная величина) при различной мощности дозы (нормализованная величина ξT0¢); там же показаны результаты расчетов по

соотношению Боуга, которые для определенной мощности дозы зависят от произведения qτ0, т.е. от полного заряда образованных импульсом ионов. При коротких импульсах облучения в обоих случаях значения f различаются незначительно; при достаточно длительных импульсах необходим учет всех указанных выше факторов, что может быть реализовано при решении системы уравнений (18.19). Наиболее существенно это различие при больших длительностях импульса, сравнимых со временем собирания ионов, когда отсутствуют условия применимости соотношения Боуга.

 

 

 

 

 

 

 

 

 

 

ξ T

/ = 1 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0 , 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ T /

= 1 0 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0 , 0 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- 9

 

- 8

 

- 7

 

- 6

 

- 5

- 4

 

- 3

 

- 2

 

- 1

 

0

 

1

1 0

 

1 0

 

1 0

 

1 0

 

1 0

 

1 0T /

1 0

 

1 0

 

1 0

 

1 0

 

1 0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Рис. 18.5. Зависимость эффективности собирания f от нормализованной

 

длительности импульса T0¢ . Данные, соответствующие приближению Боуга,

показаны пунктирными линиями; сплошные линии – расчет по соотношению

 

 

 

 

 

 

 

 

(18.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

630

 

 

 

 

 

 

 

 

 

 

Соседние файлы в предмете Инженерные конструкции