CC BY
20
БИОФИЗИКА И МЕДИЦИНСКАЯ ФИЗИКА
Определение энергетической зависимости коэффициента калибровки )д( наперстковой ионизационной камеры методом компьютерного
моделирования
A.B. Белоусов0, A.C. Осипов6
Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра физики ускорителей и радиационной медицины. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2. E-mail: а [email protected], ь [email protected] Статья поступила 26.04.2012, подписана в печать 25.07.2012.
Методом Монте-Карло моделируется энергетическая зависимость коэффициента калибровки наперстковой камеры по качеству пучка. Показано, что наибольшие отклонения наблюдаются в области энергий фотонов, соответствующих гигантскому дипольному резонансу на легких элементах. Расхождения могут достигать 10% при энергиях фотонов до 25 МэВ.
Ключевые слова: ионизационная камера наперсткового типа, метод Монте-Карло, коэффициент калибровки по качеству пучка.
УДК: 53.02. PACS: 87.53.Dq.
Введение
Ионизационные дозиметры, в частности ионизационные камеры наперсткового типа1, находят широкое применение в области радиационной медицины. Они широко используются при лучевой терапии фотонами тормозного излучения. При этом из-за специфики области применения на погрешности измерений накладываются жесткие требования. Суммарная погрешность доставки дозы, включая погрешности ее измерения, не должна превышать 5-7%.
Методику определения поглощенной дозы с использованием ионизационной камеры при облучении пучками фотонов высоких энергий водных фантомов подробно рассмотрел Hohfield [1]. Дополнительные исследования по этой теме и разработка методики проведены Andreo [2] и Rogers [3]. Процедура определения поглощенной дозы в воде, основанная на эталонах единицы поглощенной дозы в воде, реализована в национальных рекомендациях по дозиметрии [4-6] и включена в рекомендации МАГАТЭ [7]. Ионизационные камеры, в том числе наперсткового типа, калибруются по поглощенной дозе, создаваемой в воде пучками гамма-излучения 60 Со, коэффициент качества данного излучения обозначается Qo и является референсным при дальнейших измерениях.
Поглощенная доза в воде на опорной глубине zve{ для опорного пучка качества Qo дается выражением Dq0 = Mq0Nd q0 , где Mq0 — показания камеры в стандартных условиях, используемых в поверочной лаборатории; Nd,q0 — калибровочный коэффициент камеры по поглощенной дозе в воде, полученный из поверочной лаборатории. Если камера используется в пучке каче-
ства то значение поглощенной в воде дозы определяется выражением Оа = МдЫо^кд^, где коэффициент вносит поправку на различие между опорным пучком и используемым. Поправочный коэффициент на качество излучения определяется как отношение
калибровочных коэффициентов по поглощенной дозе в воде для пучков качества (2 и (¿о:
kQ,Qo =
Dq/Mq
DqJMQO ■
В идеале требуется непосредственно измерять поправочный коэффициент на качество излучения для каждой конкретной камеры при таком же качестве пучка, что и у пользователя. Однако данная процедура невозможна, так как требует наличия в поверочной лаборатории такого же, как и у конечного пользователя, клинического ускорителя. При отсутствии экспериментальных данных или невозможности непосредственного определения поправочные коэффициенты приходится рассчитывать теоретически. В большинстве случаев для пучков тормозного излучения с энергией более 5 МэВ полагают, что применима теория Брэг-га-Грея, тогда для поправочного коэффициента справедливо соотношение [2]
kQ,Qo =
QSwater)c2 (mwater)Q PQ (<SWater)c?o (®water)Q0 PQ0
которое включает следующие отношения (при качествах излучения С} и (¿о): тормозных способностей 5; средней энергии образования ионов ш; полных коэффициентов возмущения р, учитывающих все отклоне-
1 Наперстковыми (полостными) камерами называются маленькие камеры, ионизационный объем которых (газовая полость) окружен твердой стенкой (как правило, стенка является одним из электродов камеры). По измеренному ионизационному току наперстковые камеры позволяют определить дозу, поглощенную в материале стенки. Если стенки изготовлены из тканеэквивалентного материала, — дозу, поглощенную в ткани.
2 Обычно глубина 10 см в водном фантоме.
ния от идеальных условий. К сожалению, во-первых, коэффициенты возмущения р практически невозможно оценить теоретически, поскольку реальные условия измерений вне поверочных лабораторий уникальны. Во-вторых, в спектре тормозного излучения существенную часть составляет низкоэнергетическая компонента, для которой условия Брэгга-Грея нарушаются. В связи с этим пользователь практически всегда полагает, что поправочный коэффициент качества излучения равен единице во всем диапазоне энергий фотонов, используемых в лучевой терапии.
В диапазоне энергией гигантского дипольного резонанса (ГДР) продукты фотоядерных реакций, в первую очередь тяжелые заряженные частицы, также могут приводить к существенному нарушению условий Брег-га-Грея. В частности, из-за своего малого пробега они могут терять большую часть своей энергии в пределах чувствительного объема камеры. Целью работы является определение энергетической зависимости поправочного коэффициента на качество излучения методом компьютерного моделирования и определение максимальной погрешности, вносимой этим коэффициентом в результаты измерения поглощенной в водном фантоме дозы. При выполнении условий Брегга-Грея энергия, поглощенная в газовой полости ионизационной камеры, прямо пропорциональна энергии, поглощенной в окружающем полость веществе. Очень часто предполагается, что коэффициент пропорциональности р, представляющий собой среднее значение отношения тормозных способностей твердого вещества и газа, не зависит от энергии частиц. Тем не менее данная величина зависит от энергии фотонов, поскольку эффективный атомный номер вещества зависит от энергии. Кроме того, с появлением в результате фотоядерных реакций тяжелых вторичных частиц среднее отношение также меняется и зависит от доли энергии, выделившейся в результате взаимодействия тяжелых частиц. На все эти факторы влияет толщина надеваемых на камеру равновесных колпачков3, поэтому второй целью работы является определение энергетической зависимости коэффициента пропорциональности р при различных режимах облучения, т.е. при различных толщинах равновесных колпачков.
Используемые модели
Все вычисления проводятся с использованием пакета библиотек GEANT 4.9 (http://geant4.cern.ch/), который позволяет моделировать методом Монте-Карло прохождение ионизирующих излучений (в частности, фотонов) через различные материалы. В пакете предусмотрена возможность задания произвольных геометрических моделей произвольного элементного состава. В целях сокращения времени расчетов для данной работы выбран монохроматический мононаправленный источник фотонов, поперечное сечение которого представ-
ляет собой круг радиуса 10 см. В каждом отдельном вычислительном (в дальнейшем слово «вычислительный» будем опускать) эксперименте (расчет при заданных условиях) трассируется прохождение 108 первичных фотонов. В качестве фантома используется водяной куб размером 1x1x1 м. Проводятся две серии экспериментов: в первой ионизационная камера, имеющая сложный геометрический (рис. \,а) и элементный состав, помещается в водный фантом, в ее чувствительном объеме (газовой полости) фиксируются сгенерированные величины (проводятся «измерения») — количество актов ионизации и поглощенная энергия. Во второй серии «измерения» выполняются в водном фантоме, в той его области, где в предыдущей серии располагался чувствительный объем камеры (все составляющие камеры, первоначально имеющие различный элементный состав, замещаются водой), в дальнейшем независимо от эксперимента данный объем будет носить название «область интереса». В области интереса энергия обрезания4 Еси{ для фотонов и электронов равняется 10 эВ (минимальная энергия ионизации 13.6 эВ), вне области интереса энергия обрезания 10 кэВ (фотоны и электроны меньшей энергии интенсивно поглощаются в материалах стенки, слабо влияя на количество актов ионизации в области интереса). По аналогичным соображениям энергия обрезания для тяжелых частиц равна 10 эВ во всем объеме. Отметим, что
Рис. 1. Ионизационная камера наперсткового типа: а — вид изделия с равновесным колпачком; б — вид модели камеры в GEANT 4.9
° Для выполнения условий Брегга-Грея, при которых ионизационный ток камеры позволяет рассчитать поглощенную дозу в ее стенках, необходимо обеспечить условия, при которых толщина материала до газовой полости больше, чем пробег наиболее быстрых электронов, высвобожденных фотонным излучением в данном материале. Равновесные колпачки, надеваемые сверху на камеру, изменяют толщину ее стенок и позволяют измерять поглощенную дозу одной и той же камерой при различных энергиях фотонного излучения.
4 При значениях энергии ниже Еси\ частица дальше не трассируется, считается, что вся ее энергия выделяется в точке.
БИОФИЗИКА И МЕДИЦИНСКАЯ ФИЗИКА
71
программный пакет GEANT 4.9 не позволяет задавать различную энергию обрезания тяжелых частиц. Геометрическая модель ионизационной камеры, выполненная в программном пакете GEANT 4.9, представлена на рис. 1,6. В соответствии с рекомендациями [8-10] модель наперстковой камеры дополнительно окружается латунными колпачками различной толщины: в диапазоне энергий 2-6 МэВ толщина стенок равновесного колпачка 2 мм, в диапазоне энергий 6-15 МэВ толщина колпачка 4 мм, а при энергиях 15-25 МэВ — 6 мм.
Результаты и обсуждение
На рис. 2 представлен вклад продуктов фотоядерных реакций в полное количество актов ионизаций, произошедших в объеме интереса. Доля актов ионизаций, обусловленных продуктами фотоядерных реакций, в случае воды достигает максимального значения ~ 20% при энергии фотонов в районе 22 МэВ. Вид кривой схож с сечением фотоядерных реакций в области гигантского дипольного резонанса. В случае ионизационной камеры максимальное значение не превышает 7%, а кривая имеет сложную форму, что, по всей видимости, обусловлено сложным элементным составом окружающего область интереса пространства. В области энергий 15-20 МэВ значения кривых могут различаться во много раз (до 14.5).
Вклад, %
25
Вклад, % 4
20 15 10
ж
ф Вода ... ■ Камера . ♦—
♦ к
< ► V
1 1
♦ ♦ ■
**4 1** 1 ■
5000
10000 15000 Энергия, кэВ
20000 25000
Рис. 2. Вклад продуктов фотоядерных реакций в полное количество актов ионизации, произошедших в чувствительном объеме: ромбы — в случае водного фантома, квадраты — в объеме камеры
На рис. 3 представлен вклад продуктов фотоядерных реакций в дозу, поглощенную в области интереса. Вид кривых схож с кривыми вкладов в полное количество актов ионизации (рис. 2). В случае водного фантома вклад продуктов фотоядерных реакций достигает максимального значения ~ 3.5% в области энергий фотонов 22 МэВ. Если же объем интереса представляет собой чувствительный объем камеры, то вклад продуктов составляет ~ 1%. Из данных, представленных на рис. 2 и 3, видно, что нарушение условий Брэгга-Грея для тяжелых заряженных частиц — продуктов фотоядерных реакций — может привести к существенной зависимости результатов расчетов поглощенной дозы по калибровочному коэффициенту /гд0.
;
♦ ♦
♦ Вода ■ Каме
:ра
♦ ♦ 4 ►
4 ►
■
-----1 1 1
тт **1 1** М1 ♦ ■
5000
10000 15000 Энергия, кэВ
20000 25000
Рис. 3. Вклад продуктов фотоядерных реакций в поглощенную в чувствительном объеме дозу: ромбы — в случае водного фантома, квадраты — в объеме камеры
На рис. 4 показана зависимость поправочного коэффициента на качество излучения от энергии первичного пучка фотонов, полученная по результатам компьютерного моделирования. Кривая нормирована на значение коэффициента качества /гд0. Видно, что вплоть до энергий первичного фотонного пучка ~ 20 МэВ отклонение поправочного коэффициента на качество излучения от опорного значения не превышает 2%. В области энергий выше 20 МэВ характер зависимости изменяется, кривая плавно уменьшается и отклонение достигает ~ 10% при энергии 25 МэВ. Данный ход энергетической зависимости обусловлен увеличением вклада продуктов фотоядерных реакций в поглощенную объемом интереса дозу, что связано с изменением количества и энергии образовавшихся тяжелых частиц. Увеличение энергии, поглощенной объемом интереса в результате взаимодействия тяжелых частиц, приводит к изменению среднего отношения тормозных способностей всех частиц. Минимальное количество событий, зарегистрированных в области интереса, не менее 30000, относительная погрешность результатов в каждой точке не превышает 0.5%.
là1-0 PQ О
H
о о >3 I
A S M я"
g S, 0.5
CQ **
о
Он ^
в s
s * &
к
II1M * Éiè^
5000
10000 15000 Энергия, кэВ
20000 25000
Рис. 4. Энергетическая зависимость поправочного коэффициента на качество излучения
1а.
Ё и
■е
п £
10 8 6 4 2
25
а
20
1а.
I 15
■е
10
40
30
1о_
Ё и Я Я 15
•е -е-
п £
20
10
0.5 1.0 1.5 2.0 Энергия, МэВ
2.5
■ ■ ■
■ ■ ■ ■ 1 ■
в
70
60
, 50 1а.
Ё
8 40 я
т &
30
20
10
2 4 6
Энергия, МэВ
0 5 10 15 20 0 10 20 30
Энергия, МэВ Энергия, МэВ
Рис. 5. Энергетическая зависимость коэффициента р при различных толщинах равновесных колпачков: а — без колпачка; б — латунный колпачок толщиной 2 мм; в — латунный колпачок толщиной 4 мм; г — латунный
колпачок толщиной 6 мм
На рис. 5 представлено изменение коэффициента р в различных энергетических диапазонах при использовании латунных равновесных колпачков различной толщины. В диапазоне энергий 1-2 МэВ (без дополнительного равновесного колпачка) энергетическая зависимость крайне слабая, отклонения от среднего значения не превышают 7% (рис. 5,а). При энергиях 2-6 МэВ (толщина равновесного колпачка 2 мм) значение коэффициента скачкообразно изменяется, и кривая монотонно спадает по мере роста энергии фотонов. В данном энергетическом диапазоне значения коэффициента изменяются примерно в 2 раза (рис. 5,6). В диапазоне энергий 6-15 МэВ (толщина 4 мм) значение коэффициента снова скачкообразно изменяется примерно в 3 раза, кривая немонотонна, имеет множество максимумов и минимумов, что связано с различными вкладами фотоядерных взаимодействий в объеме интереса. Отклонения коэффициента ~р в данном энергетическом диапазоне не превышает 15% (рис. 5, в). При использовании колпачков толщиной 6 мм, что соответствует энергетическому диапазону 15-25 МэВ, зависимость выражена более ярко, как и прежде, наблюдается скачкообразное изменение примерно в 2 раза по сравнению с колпачком толщиной 4 мм. Значения в начале и в конце промежутка отличаются примерно
в 1.3 раза (рис. 5, г). Погрешность результатов вычислений в каждой точке не превышает 1%. Наблюдаемые скачкообразные изменения связаны в первую очередь с изменением эффективного атомного номера слоя вещества, окружающего газовую полость. Энергетические зависимости в исследуемых диапазонах обусловлены, во-первых, зависимостью эффективного атомного номера от энергии фотонов, а во-вторых, различными вкладами фотоядерных реакций.
Заключение
Поправочный коэффициент на качество излучения слабо зависит от энергии монохроматического пучка фотонов и практически не зависит от энергии вплоть до 20 МэВ (максимальные отклонения не превышают 2%). При дальнейшем увеличении энергии возрастает разница между поглощенной дозой, обусловленной продуктами фотоядерных реакций, и значения коэффициента плавно снижаются, уменьшаясь
примерно на 10% по сравнению со значением для опорного излучения 60Со.
Коэффициент ¿7 скачкообразно изменяется по мере изменения толщины равновесных колпачков, надетых на ионизационную камеру наперсткового типа, что связано главным образом с изменением среднего атомного
состава вещества, окружающего газовую полость. При отсутствии равновесного колпачка (энергия фотонов до 2 МэВ) данный коэффициент меняется незначительно (менее 7%), и можно полагать его постоянным. В интервале 2-6 МэВ значения коэффициента меняются примерно в 2 раза, в диапазоне энергий 6-15 МэВ — примерно в 1.15 и в промежутке 15—25 МэВ — примерно в 1.3 раза. Наиболее явно энергетические зависимости выражены при наличии колпачков толщиной 2 и 6 мм.
Список литературы
1. Hohfield К. 11 Dosimetry in Radiotherapy. Vol. 1. Vienna, 1998.
2. Andreo P. Ц Phys. Med. Biol. 1992. 37. P. 2189.
3. Rogers D.W. 11 Med. Phys. 1992. 19. P. 1227.
4. AAPM TG-51. // Med. Phys. 1999. 26. P. 1847.
5. Deutsches Institute fur Normung, Dosismessverfahren nach der Sondenmethode für Photonen- und Electronenstrahlung.
Teil 2: lonisationsdosimetrie. Deutsche Norm DIN 6800-2. Berlin, 1997.
6. Institute of Physical Science in Medicine. // Phys. Med. Biol. 1990. 35. P. 1355.
7. IAEA. The use of plane-parallel ionization chambers in high-energy electron and photon beams: An international code of practice for dosimetry. Technical Reports Series № 381. Vienna, 1997.
8. МАГАТЭ. Серия технических докладов № 398. Определение поглощенной дозы при дистанционной лучевой терапии. Международные практические рекомендации по дозиметрии, основанные на эталонах единицы поглощенной дозы в воде. Вена, 2004.
9. Dutreix A., Bjarngard В.Е., Bridier A. et al. 11 ESTRO Booklet N 3. Monitor Unit Calculation for High Energy Photon Beams. Physics for Clinical Radiotherapy. ESTRO, 1997.
10. Mijnheer В., Bridier A., Garibaldi С. 11 ESTRO Booklet N 6. Monitor Unit Calculation for High Energy Photon Beams. Practical Examples. Physics for Clinical Radiotherapy. ESTRO, 2001.
Definition of energy dependence of thimble ionisation chamber calibration factor kq by the method of computer modelling
A. V. Belousov , A. S. Osipov
Department of Accelerator Physics and Radiation Medicine, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.
E-mail: a [email protected], b [email protected].
In present work as a method Monte Carlo energy dependence of beam quality calibration factor of thimble ionization chamber is modelled. It is shown, that the greatest derivation are observed in area photons energy corresponding giant dipole resonance on easy elements. Divergences can reach 10% at photons energy to 25 MeV.
Keywords: thimble ionization chamber, Monte Carlo method, beam quality calibration factor. PACS: 87.53.Dq. Received 26 April 2012.
English version: Moscow University Physics Bulletin 1(2013).
Сведения об авторах
1. Белоусов Александр Витальевич — канд. физ.-мат. наук, ст. преподаватель; тел.: (495) 939-49-46, e-mail: [email protected].
2. Осипов Алексей Сергеевич — ассистент; тел.: (495) 939-49-46, e-mail: [email protected].
