Краткое описание работы
Актуальность работы обусловлена необходимостью точного прогнозирования распределения поглощённых доз в тканях и материалах при планировании лучевой терапии и радиационной защите. Использование прецизионных симуляций, таких как пакет GATE, основанный на методе Монте-Карло, позволяет минимизировать риски для пациента и повысить эффективность лечения, а также оптимизировать конструкции защитных экранов.
Объектом исследования является процесс взаимодействия гамма-излучения с веществом, а предметом — дозовое распределение, формируемое в фантомах из различных материалов при облучении источниками с разной энергией.
Цель работы: смоделировать и проанализировать дозовое распределение гамма-излучения в программной среде GATE, выявив влияние энергии фотонов и материала фантома на характер поглощения энергии.
Для достижения цели были поставлены следующие задачи:
1. Изучить теоретические основы взаимодействия гамма-излучения с веществом (фотоэффект, комптоновское рассеяние, образование пар) и принципы метода Монте-Карло.
2. Освоить алгоритм работы в среде GATE: создание геометрии, задание источников, материалов фантома и детекторов.
3. Провести серию симуляций для трёх источников с различными энергиями гамма-квантов.
4. Исследовать изменение дозового распределения при замене материала фантома на свинец.
5. Выполнить анализ полученных результатов и сформулировать выводы.
Основная идея (выводы): В ходе работы установлено, что GATE корректно воспроизводит физические закономерности: для низких энергий преобладает фотоэффект (доза резко падает в глубине), для высоких — комптон-эффект и образование пар (доза распределяется более равномерно). Замена тканеэквивалентного фантома на свинец продемонстрировала резкое увеличение поглощения в поверхностных слоях и снижение глубинного проникновения, что подтверждает эффективность свинца как защитного материала. Таким образом, моделирование в GATE позволяет прогнозировать радиационные поля и оптимизировать параметры облучения в зависимости от энергии источника и состава среды.
Название университета
КУРСОВАЯ РАБОТА НА ТЕМУ:
МОДЕЛИРОВАНИЕ ДОЗОВОГО РАСПРЕДЕЛЕНИЯ В ПРОГРАММА GATE.1. ТЕОРЕТИЧЕСКОЕ ВВЕДЕНИЕ(ГАММА ИЗЛУЧЕНИЕ,ФОТОЭФФЕКТ,ОБРАЗОВАНИЕ ПАР,КОМПТОН ЭФФЕКТ, МЕТОД МОНЕ КАРЛО )2.МЕТОД МОНТЕ-КАРЛО,АЛГОРИТМ РАБОТЫ,АНАЛИЗ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ,ИЗМЕНЕНИЕ МАТЕРИАЛА ФАНТОМА,СВИНЕЦ,ИСТОЧНИКИ С ДРУГОЙ ЭНЕРГИЕЙ,3 ИСТОЧНИКА С РАЗЛИЧНЫМИ ЭНЕРГИЯМИ,ВЫВОД
г. Москва, 2025 год.
Содержание
Введение
1⠄Теоретические основы взаимодействия гамма-излучения с веществом и метод Монте-Карло
1⠄1⠄Физические механизмы взаимодействия гамма-излучения: фотоэффект, комптоновское рассеяние и образование пар
1⠄2⠄Основные принципы и алгоритм метода Монте-Карло в задачах радиационной физики
1⠄3⠄Обзор программного пакета GATE: архитектура, возможности и применение для моделирования дозовых распределений
2⠄Моделирование дозовых распределений в программе GATE: $$$$$$$$$$ $$$$$$ и $$$$$$ $$$$$$$$$$$
2⠄1⠄$$$$$$$$$$ $$$$$$ $$$$$$$ и $$$$$$$$$: $$$$$$$ $$$$$$$$$ ($$$$, $$$$$$) $$ $$$$$$$$$$$$ $$$$$$$$ $$$$
2⠄2⠄$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$ гамма-излучения: $$$$$$$$$$$$$ $$$$$$$$$$ с $$$$$$$$$$ $$$$$$$$$
2⠄3⠄$$$$$$ $$$$$$$$$$$ моделирования для $$$$ $$$$$$$$$$ с $$$$$$$$$$ $$$$$$$$$ и $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$$$$ $$$$
$$$$$$$$$$
$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$
Введение
Ионизирующее излучение, в частности гамма-излучение, играет ключевую роль в современной ядерной медицине, радиационной терапии и промышленной дефектоскопии. Точное знание пространственного распределения поглощенной дозы в облучаемом объекте является критически важным условием для эффективного лечения онкологических заболеваний, обеспечения радиационной безопасности персонала и оптимизации технологических процессов. Традиционные экспериментальные методы измерения дозовых полей, такие как использование ионизационных камер или пленочных дозиметров, зачастую трудоемки, дороги и не всегда позволяют получить детальную трехмерную картину распределения энергии в сложных гетерогенных средах. В этой связи особую актуальность приобретают методы компьютерного моделирования, в частности метод Монте-Карло, который признан «золотым стандартом» для расчета переноса излучения благодаря своей способности учитывать стохастическую природу взаимодействия фотонов с веществом.
Программный пакет GATE (Geant4 Application for Tomographic Emission), основанный на библиотеке Geant4, является мощным и гибким инструментом для симуляции физических процессов в ядерной медицине и радиационной физике. Однако, несмотря на широкое распространение, задача точного и вычислительно эффективного моделирования дозовых распределений в различных средах и для различных источников излучения остается не до конца решенной. Проблематика данной работы заключается в необходимости систематического анализа того, как изменение энергии гамма-излучения и материала облучаемого фантома влияет на характер формирующегося дозового поля. Существующие исследования часто фокусируются на узком диапазоне параметров, что не позволяет сформировать целостное представление о закономерностях, управляющих процессом дозообразования.
Объектом данного исследования является процесс взаимодействия гамма-излучения с веществом и формирование дозового $$$$$$$$$$$$$. $$$$$$$$$ исследования $$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$, $$$$$$$$$$ $$$$$$$ $$$$$-$$$$$ $ $$$$$ $$$$, $$ $$$$$$$ $$$$$$$$$ гамма-излучения и $$$$$$$$$ $$$$$$$ ($$$$ и $$$$$$).
$$$$$ $$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$ $ $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$-$$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$ $$$$.
$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$$$$$ $$$$$$ $$$$$$$$$ $$$$$$:
$. $$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$, $$$$$$$ $$$$$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$ $$$, $ $$$$$ $$$$$$$$ $$$$$$ $$$$$-$$$$$.
$. $$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$ $$$ $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$.
$. $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$, $$$$$$$$$$ $$$$$$$$$ $$$$$$$ ($$$$, $$$$$$) $ $$$$$$$$ $$$$$$$$ $$$$$-$$$$$$$$$.
$. $$$$$$$$ $$$$$ $$$$$$$$$ $$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$.
$. $$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$ $$$$$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$ $$ $$$$$$$$$$$ $$$$$$$$$$ $$ $$$$$ $ $$$$$$$$ $$$$$$$$ $$$$.
$$$$$$$$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$ $$$$$ $$$$$-$$$$$, $$$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$ $$$$. $$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$. $$$$$$$$$$$$$ $$$$$ $$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$ $ $$$$$$$$$ $$$$$$ $$ $$$$$$$ $$$$$$$$$$.
$$$ $$$$$$$$$ $$$$$$ $$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$$$$$$ $$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$ $$ $$$$$$$$$$$$$ $$$$$$$$ ($$$$$$$ $$ $$$$$$$$ $$$ $$$$$$$, $$$$$$$ $$$$$$$) $ $$$$$$$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$ $$$$.
Физические механизмы взаимодействия гамма-излучения с веществом: фотоэффект, комптоновское рассеяние и образование пар
Гамма-излучение представляет собой поток фотонов с чрезвычайно высокой энергией, испускаемых при ядерных переходах, аннигиляции частиц или торможении заряженных частиц. В отличие от заряженных частиц, гамма-фотоны не имеют массы покоя и электрического заряда, что обуславливает их высокую проникающую способность и специфический характер взаимодействия с веществом. Понимание механизмов этого взаимодействия является фундаментальной основой для моделирования дозовых распределений в любых радиационных задачах, включая лучевую терапию, радиационную безопасность и ядерную медицину [12].
Основными процессами, ответственными за ослабление гамма-излучения в веществе, являются фотоэлектрический эффект, комптоновское рассеяние и эффект образования электрон-позитронных пар. Вероятность протекания каждого из этих процессов определяется энергией падающего фотона и атомным номером материала среды. При низких энергиях (до нескольких сотен килоэлектронвольт) доминирует фотоэффект, при средних энергиях (от нескольких сотен кэВ до нескольких МэВ) – комптоновское рассеяние, а при высоких энергиях (свыше 1,022 МэВ) становится возможным образование пар.
Фотоэлектрический эффект представляет собой процесс, при котором падающий гамма-фотон полностью поглощается атомом, передавая свою энергию одному из связанных электронов внутренних оболочек (как правило, K- или L-оболочки). Электрон, получивший энергию, превышающую его энергию связи, выбивается из атома, становясь свободным фотоэлектроном. Кинетическая энергия фотоэлектрона равна разности энергии падающего фотона и энергии связи электрона в атоме. Вероятность фотоэффекта резко возрастает с увеличением атомного номера материала (пропорционально Z⁴–Z⁵) и уменьшается с ростом энергии фотона (пропорционально E⁻³·⁵). Это делает фотоэффект крайне важным для материалов с высоким Z, таких как свинец, используемый в радиационной защите, и объясняет характерные скачки поглощения на K-краях элементов. После выбивания электрона его место занимает электрон с вышележащей оболочки, что сопровождается испусканием характеристического рентгеновского излучения или оже-электронов. Данный процесс вносит существенный вклад в локальное поглощение энергии, особенно вблизи поверхности облучаемого объекта.
Комптоновское рассеяние, или комптон-эффект, является доминирующим механизмом взаимодействия гамма-излучения с веществом в диапазоне энергий от 0,1 до 10 МэВ, что соответствует энергиям, наиболее часто используемым в радиотерапии. В этом процессе падающий фотон взаимодействует со свободным или слабосвязанным электроном внешней оболочки атома. В результате взаимодействия фотон передает часть своей энергии электрону, который выбивается из атома (комптоновский электрон), а сам фотон рассеивается под некоторым углом, изменяя свое $$$$$$$$$$$ $ $$$$$$$$ $$$$ $$$$$$$. $$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$ от $$$$ $$$$$$$$$: $$$ $$$$$$$$$ $$$$$ ($$$$ $$$°) фотон $$$$$$ $$$$$$$$$$$$ $$$$$$$, а $$$ $$$$$$$$$ $$$$$$ ($$$$$ $$$$) – $$$$$$$$$$$. $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$ от $$$$$$$$ $$$$$$ $$$$$$$$$ ($$$$$$$$$$$$$$$ $) $ $$$$$$$$$$$ с $$$$$$ энергии $$$$$$. $$$$$$$$$$ своей $$$$$$$ $$$$$$$$$$$$$ рассеяние $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$, $$$$$$$ $$$$$ $$$$$$$$$$$$$$$$$ с веществом $$$$$, $$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$ излучения $ $$$$$$$$$$$$ $$$$$$$$ $$$$. $$$$$$ $$$$ $$$$$$$ $$ $$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$ в $$$$$$ $$$$$$$$, $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$.
$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$, $$$$$ $$$$$$$ $$$$$$$$$ $$$$$-$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$ $$$$$$$$$, $$ $$$$ $,$$$ $$$. $ $$$$ $$$$$$$$ $ $$$$ $$$$$$$$$$$ $$$ $$$$$$$$ $$$$ ($$$$ – $$$$$$$$$) $$$$$$$$ $$$$$ $$$$$$$$$$$$ $ $$$$ $$$$$$$-$$$$$$$$$$$: $$$$$$$$ $ $$$$$$$$. $$$$$$$ $$$$$$$ $$$$$$ $$$$$ $$$$$$ $,$$$ $$$ $$$$$$$$$$$$$$ $$$$$ $$$$$$$$$$ $ $$$$$$$$$$ $ $$$$ $$$$$$$$$$$$ $$$$$$$. $$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$ $$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$ $$$$$$ ($$$$$$$$$$$$$$$ $$($)) $ $$$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$ ($$$$$$$$$$$$$$$ $$). $$$$$$$$$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$ $ $$$$$ $ $$$$$$$$$$$$ $ $$$$$ $$ $$$$$$$$$$, $$$ $$$$$$$$ $ $$$$$$$$$$ $$$$ $$$$$$$$$$$$$$$ $$$$$-$$$$$$$ $ $$$$$$$$ $$$ $$$ $$$$$$, $$$$$$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$$$$. $$$ $$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$$, $$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$$$ $$$$$$$, $$$ $$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$ $$$ $ $$$$$$$$$$$ $$$$ $$$$$$$$$$ $$$$$$$$$$$$$, $$$$$$$$ $ $$$$$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$ [$$].
$$$$$$$$$$$$$ $$$$$ $$$$$$$ $$ $$$$ $$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$-$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$. $$$ $$$$$$ $$$$$$$$$$, $$$$$ $$$ $$$$ ($$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$ $$$$$$), $ $$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$ $$$$ $$ $$$$$$$$$ ($$$$$$ $$$$$$$$$$) $ $$$$$$$$$$$$ $$$$$. $$$ $$$$$$$ $$$$$$$$$$, $$$$$ $$$ $$$$$$, $$$$$$$$$$$$ $ $$$$$$$$$$$$ $ $$$$$$$$ $$$$$$$, $$$ $$$ $$ $$$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$ $$$$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$$$$. $$$ $$$$$$$ $$$$$$$$ ($$$$ $$ $$$) $ $$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$, $$$ $$$$$ $$$$$$ $$ $$$$$ $$$$$$$$ $$$$$$$. $$$$$$$$$ $$$$ $$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ [$$].
Для полного понимания процессов формирования дозовых распределений необходимо рассмотреть не только индивидуальные механизмы взаимодействия, но и их совокупное влияние на ослабление гамма-излучения при прохождении через среду. Ослабление пучка гамма-фотонов в веществе подчиняется экспоненциальному закону, однако это справедливо лишь для узкого моноэнергетического пучка в условиях «хорошей геометрии». В реальных задачах радиационной физики, включая моделирование в программе GATE, необходимо учитывать многократное рассеяние и образование вторичных частиц, что делает аналитическое описание крайне сложным и обуславливает необходимость применения численных методов, таких как метод Монте-Карло.
Каждый из трех фундаментальных процессов взаимодействия характеризуется своим микроскопическим сечением, которое определяет вероятность взаимодействия на единицу длины пути фотона в среде. Полное сечение взаимодействия складывается из сечений фотоэффекта, комптоновского рассеяния и образования пар. Коэффициент линейного ослабления μ, являющийся произведением полного сечения на плотность числа атомов в среде, представляет собой макроскопическую характеристику, определяющую способность материала ослаблять гамма-излучение. Важно отметить, что при прохождении через среду фотоны могут испытывать не одно, а серию последовательных взаимодействий, особенно в случае комптоновского рассеяния, когда рассеянный фотон с меньшей энергией продолжает движение и может взаимодействовать далее.
Особый интерес для моделирования дозовых распределений представляет пространственное распределение энергии, переданной веществу в результате этих взаимодействий. При фотоэффекте вся энергия падающего фотона (за вычетом энергии связи электрона) передается локально в точке взаимодействия в виде кинетической энергии фотоэлектрона. Однако часть энергии может быть унесена характеристическим рентгеновским излучением, если образующаяся вакансия заполняется электроном с вышележащей оболочки. При комптоновском рассеянии энергия распределяется между рассеянным фотоном и комптоновским электроном, причем рассеянный фотон может перенести энергию на значительное расстояние от точки первичного взаимодействия, что приводит к эффекту накопления дозы на некоторой глубине. При образовании пар энергия распределяется между электроном и позитроном, а последующая аннигиляция позитрона приводит к появлению двух фотонов с энергией 511 кэВ, которые также могут переносить энергию на расстояние.
Для количественного описания передачи энергии веществу используется понятие линейной передачи энергии (ЛПЭ), которая характеризует энергию, переданную среде заряженной частицей на единице длины пути. Для гамма-излучения ЛПЭ не является прямой характеристикой, так как фотоны сами не создают непрерывного трека ионизации, а передают энергию опосредованно через вторичные электроны. Поэтому для дозиметрии гамма-излучения ключевым понятием является керма (kinetic energy released per unit mass), которая определяет сумму начальных кинетических энергий всех заряженных частиц, образованных фотонами в единице массы вещества. В условиях электронного равновесия, которое устанавливается на некоторой глубине в облучаемом объекте, керма равна поглощенной дозе.
Энергетическая зависимость относительного вклада различных процессов имеет решающее значение для понимания формы дозовых распределений в разных материалах. Для воды, которая является стандартным тканеэквивалентным материалом, в диапазоне энергий от 100 кэВ до 10 МэВ комптоновское рассеяние является доминирующим процессом, на его долю приходится более 80% всех взаимодействий. Это приводит к характерной форме глубинного дозового распределения с максимумом на некоторой глубине (область накопления), за которым следует экспоненциальный спад. Для свинца, напротив, фотоэффект играет существенную роль вплоть до энергий порядка 1 МэВ, что обеспечивает эффективное поглощение низкоэнергетического излучения вблизи поверхности. При энергиях выше 5 МэВ в свинце начинает доминировать образование пар, что также влияет на форму дозового профиля [27].
При моделировании дозовых распределений в программе GATE необходимо корректно учитывать все три механизма взаимодействия, так $$$ $$$$$$ $$ $$$ $$$$$$ $$$$ $$$$$ в $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$. $$$$$$$$$ GATE, $$$$$$$$$$ $$ $$$$$$$$$$ $$$$$$, $$$$$$$$$$$$$ $$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$. $$$$$$$$$$$$ $$$$$ $$$$$$$$ $$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$$$$$$$$$$$ $$$$$$$$, $$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$ $$$$$$$$$$$ $$$$$$, $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$ в $$$$$$$$$ $$$$$$$ $$ $$ $$$, $$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$, $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$ $ $$$$$$$ $$$$$$$$$$$$$.
$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$$ ($$$$$ $$ $$$) $ $$$$$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$, $$$$$ $$$ $$$$$$, $$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$, $$$ $$$$$$$ $$$$$-$$$$$ $$$$$$$$ $$$$$$ $$ $$$$. $$$$$$ $$$$$$$$$$$ $$$$$ $$$$$$$$ $$ $$$$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$, $$$$$$$ $ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$ $$$$$$ $$$$$$$$$$$$. $$$ $$ $$$$$, $$$ $$$$$$$$$$$$$ $$$$$$ $$ $$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$ $$$$$$$ $$$$$ $$$$$ $$$$$$$$, $$$ $$$ $$ $$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$.
$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$. $$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$-$$$$$$, $$$$$$$ $$$$$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$, $ $$$ $$$$$$ – $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$. $$$ $$$$$$$$$$$ $$$ $$$$$$$$ $ $$$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$, $$$$$$$$ $$$ $$$$$$$ $$$$$$$$. $$$ $$$$$$$ $$$$$$$$$$$ $$$$$$ $$ $$$$$ $$$$$$$$ $$$$, $$$$$$$$ $$ $$$$$ $$$$$ $ $ $$$$$$$ $$$$$$$$.
$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$ $$$$$$$$$$. $$$$ $$$$$-$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$, $$$$$$ $$$$$$$$$ $$$$$$$$$, $$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$ $$$, $$$$$$$$$$$$ $$ $$$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$. $$$$$$ $$$$$$$$$ $$$$$$$$$$ $ $$$$$ $$$$$$$ $$ $$ $$$$$$$ $ $$$$$$$$$. $$$ $$$$ $$$$$$ $$$$$$$$$$ $ $$$$$$$$ $ $$$ $$$$$$$$$$ $$$$$ $ $$, $ $ $$$$$$$$ $$ $$$ – $$$$$ $ $$. $$$ $$$$$$$$, $$$ $$$$$$$, $$$$$$$$$$ $$$$$$$ $ $$$$$ $$$$$, $$$$$ $$$$ $$$$$$$$$ $$ $$$$$$$$$ $$$$$$$$$$ $$ $$$$ $$$$$, $$$ $$$$$$$$ $ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$.
$$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ – $$$$$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$ $$$ – $$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$. $$$$$$ $$ $$$$ $$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$$$, $$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$. $$$ $$$$, $$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $ $$$$ $$$$$$$$$$ $$$$$$$$$$. $$$ $$$$$$, $$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$, $$$$$$$$$$$$ $$$$ $$$$$$ $$$$$$$$$$ $$$ $$$$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$ $$$ $$$$$$$ $$$$$$$$, $$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$. $$$$$$$$$ $$$$ $$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ [$].
Основные принципы и алгоритм метода Монте-Карло в задачах радиационной физики
Метод Монте-Карло представляет собой численный метод решения математических задач путем моделирования случайных процессов и построения статистических оценок. В радиационной физике данный метод является основным инструментом для моделирования переноса ионизирующего излучения через вещество, поскольку он позволяет с высокой точностью воспроизводить стохастическую природу взаимодействия частиц с атомами и молекулами среды. В отличие от детерминистических методов, решающих уравнение переноса излучения аналитически или численно, метод Монте-Карло имитирует траектории отдельных частиц, используя генераторы случайных чисел для розыгрыша каждого акта взаимодействия [6].
Фундаментальным принципом метода Монте-Карло в задачах радиационной физики является моделирование индивидуальных историй частиц. Каждая история начинается с генерации первичной частицы (фотона, электрона, нейтрона и т.д.) с заданными начальными параметрами: энергией, направлением движения и положением в пространстве. Затем частица перемещается в среде, и на каждом шаге моделирования с помощью генератора случайных чисел определяется, произойдет ли взаимодействие с веществом, и если да, то какого типа. Для этого используются известные из физики сечения взаимодействий, которые определяют вероятности различных процессов. После каждого взаимодействия параметры частицы обновляются, и процесс повторяется до тех пор, пока частица не покинет область моделирования, не поглотится или ее энергия не станет ниже заданного порога отсечки.
Алгоритм работы метода Монте-Карло для моделирования переноса гамма-излучения можно представить в виде последовательности следующих шагов. Первым этапом является задание начальных условий: определение геометрии моделируемой системы, состава и плотности материалов, а также характеристик источника излучения. Вторым этапом является генерация первичного фотона с розыгрышем его начальных параметров в соответствии с заданным распределением источника. Третьим этапом является определение длины свободного пробега фотона до следующего взаимодействия, которая вычисляется на основе полного коэффициента ослабления материала. Четвертым этапом является розыгрыш типа взаимодействия (фотоэффект, комптоновское рассеяние, образование пар) на основе относительных вероятностей каждого процесса. Пятым этапом является моделирование выбранного взаимодействия: для фотоэффекта – определение энергии и направления вылета фотоэлектрона и возможного характеристического излучения; для комптоновского рассеяния – определение энергии и направления рассеянного фотона и комптоновского электрона; для образования пар – определение энергий и направлений электрона и позитрона. Шестым этапом является сохранение информации о переданной энергии в точке взаимодействия для последующего расчета дозового распределения. Седьмым этапом является повторение процесса для всех вторичных частиц, образовавшихся в результате взаимодействия. Восьмым этапом является повторение всей процедуры для большого числа первичных частиц (обычно от 10⁶ до 10⁹) для получения статистически значимых результатов.
Одной из ключевых особенностей метода Монте-Карло является его сходимость. Статистическая погрешность результатов пропорциональна 1/√N, где N – количество смоделированных историй. Это означает, что для уменьшения погрешности вдвое необходимо увеличить количество историй в четыре раза, что приводит к значительным вычислительным затратам. Поэтому при моделировании в программе GATE важно выбирать оптимальное количество историй, обеспечивающее приемлемую точность при разумном времени расчета. Для задач медицинской физики, таких как расчет дозовых распределений в лучевой терапии, статистическая погрешность обычно не должна $$$$$$$$$ 1-$% в $$$$$$$ $$$$$$$$$ $$$$.
$$$ $$$$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$$-$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$. $ $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$: $$$$$ $$$$$$$$$$$ ($$$$$$$$$), $$$ $$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$ $$$$$$$$$ $$$$$; $$$$$ $$$$$$$ $$$$$$$$, $$$ $$$$$$$ $$$$$$$ $ $$$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$; $$$$$ $$$$$$$$ ($$$$$$$$$$ $$$$$$$$), $$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$ $$$$$ $$$$$$$ $ $$$$$$$$$$$$ $$$$$$$; $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ ($$$$$$ $$$$$$$$$$$$), $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$$. $ $$$$$$$$$ $$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$ $$ $$$$$$$$ $$$$$$.
$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$ $$$$ $$$$$$$$$ $$$$$$. $$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$ $ $$$$$$$$$ $$$ $$$$$$$$$$$ $$$, $ $$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$. $$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$$ $$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $ $$$$$$$$$ $$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $$$$$$, $$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$. $$$$$$$ $ $$$$$$$$ $$$$ $$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$ $$$$$, $ $$ $$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$ $$$$.
$$$$$$ $$$$$$$$ $$$$$$ $$$$$-$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$. $ $$$$$$$ $$ $$$$$$$, $$$$$$$$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$, $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$ $$$$$ $$$$$$$$$$$. $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$$$$ $$$$$$. $$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$, $$$$$ $$$ $$$$$ $$$$$$$$$$$$ $$$$$$$$$$ ($$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$$$$$$, $$$$) $ $$$$$ $$$$$$$$$$$$$ $$$$$$$$$. $ $$$$$$$$$ $$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$ $$$$$$ $$$$ $$$$$$$$$$$$$ [$$].
$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$-$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$$$, $$$$$$$$$$$ $$$$$$$ $$$$$$$$. $ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$, $$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$$$$$$$$ $$$$$$. $ $$$$$ $$$$$$$ $$$$$$$$$: $$$$$$ $$$$$$$$$$ $$$$$$ $$$$$ $$$$$-$$$$$$$$$, $$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$$$$$ $$$ $$$$$$$$$$$ $$$$$ $$$$$$$$$, $$$$$$ $$$$$$ $$ $$$$$$$$$. $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$.
$ $$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$ $$$$$ $$$$$-$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$ $$$$ $ $$$$$$ $$$ $$$$$$$$$ $$$$$-$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$. $$$$$ $$$$$$ $$$$$-$$$$$ $$$$$$$$$$ $$$ $$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$$$, $$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$ $$$$$. $$$$$$$$$ $$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$ $$$$$$$ $$$$$$$$$, $$$$$$$$$$, $$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$, $ $$$$$ $$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$$. $$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$, $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ ($$$$ $ $$$$$$), $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$-$$$$$$$$$ ($$$$$$$, $$$, $$$$$$$$$ $$$$$), $$$$$ $$$$$$$$$$ $$$$$$$, $$$$$$ $$$$$$$$$ $$$ $$$$$$$$$ $$$$$ $$$$$$$ $ $$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$.
Для практической реализации метода Монте-Карло в задачах моделирования переноса излучения разработано несколько специализированных программных пакетов, среди которых наиболее известными являются MCNP (Monte Carlo N-Particle), EGSnrc (Electron Gamma Shower), FLUKA и Geant4 (GEometry ANd Tracking). Каждый из этих пакетов имеет свои особенности, преимущества и области применения. Однако для задач медицинской физики, особенно связанных с моделированием в ядерной медицине и лучевой терапии, особую популярность приобрел программный пакет GATE (Geant4 Application for Tomographic Emission), который является надстройкой над библиотекой Geant4 и предоставляет специализированные инструменты для моделирования в медицине [14].
Программа GATE была разработана международным консорциумом OpenGATE и впервые представлена в 2004 году. Основной целью разработки было создание удобного инструмента для моделирования позитронно-эмиссионной томографии (ПЭТ) и однофотонной эмиссионной компьютерной томографии (ОФЭКТ). Однако со временем функциональность GATE была расширена, и в настоящее время программа успешно применяется для моделирования широкого спектра задач, включая расчет дозовых распределений в лучевой терапии, моделирование радиационной защиты, проектирование детекторов ионизирующего излучения и обучение студентов и специалистов в области медицинской физики.
Архитектура GATE основана на модульном принципе. Программа предоставляет пользователю набор макрокоманд на языке C++ и специальном скриптовом языке, которые позволяют описывать геометрию моделируемой системы, задавать материалы, определять характеристики источника излучения, выбирать физические модели и настраивать параметры вывода результатов. Важной особенностью GATE является использование системы временных окон (time slices), которая позволяет моделировать динамические процессы, такие как движение источника или изменение активности радиофармпрепарата с течением времени. Для задач данной курсовой работы, связанных с моделированием статических дозовых распределений, временные окна не используются, однако данная возможность демонстрирует гибкость и универсальность программы.
Для описания геометрии в GATE используется система объемных примитивов (боксы, сферы, цилиндры, конусы и т.д.), а также возможность создания сложных геометрических форм с помощью булевых операций (объединение, пересечение, вычитание). Каждый геометрический объект может быть помещен в иерархическую структуру, позволяющую моделировать сложные системы, состоящие из множества компонентов. Для данной курсовой работы используются простые геометрические формы: прямоугольный параллелепипед для фантома из воды и свинца, а также точечный источник гамма-излучения, расположенный на заданном расстоянии от поверхности фантома.
Задание материалов в GATE осуществляется путем указания их химического состава и плотности. Программа содержит встроенную базу данных материалов, включающую основные элементы и стандартные материалы, такие как вода, воздух, свинец, медь, алюминий и другие. Для каждого материала GATE автоматически рассчитывает необходимые физические характеристики, такие как сечения взаимодействий, на основе данных из библиотеки Geant4. Для данной работы используются два основных материала: вода (плотность 1,0 г/см³) как имитатор биологической ткани и свинец (плотность 11,34 г/см³) как материал с высоким атомным номером, используемый в радиационной защите.
Выбор физических моделей является одним из наиболее ответственных этапов настройки моделирования в GATE. Для описания электромагнитных взаимодействий гамма-излучения с веществом программа предоставляет несколько наборов физических моделей (physics lists), которые различаются точностью описания процессов и вычислительными затратами. Для задач медицинской физики, где энергии фотонов обычно не превышают 20 МэВ, рекомендуется использование набора моделей «emstandard_opt3» или «emlowenergy». Набор «emstandard_opt3» обеспечивает хороший баланс между точностью и скоростью расчетов, в то время как низкоэнергетические модели («emlowenergy») обеспечивают более точное описание фотоэффекта и комптоновского рассеяния при низких энергиях, но требуют большего времени счета. Для данной курсовой работы, учитывая диапазон энергий источников от 0,1 до 10 МэВ, целесообразно использование набора «emstandard_opt3» [30].
Важным аспектом моделирования в GATE является задание пороговых энергий для генерации вторичных $$$$$$. $$$ $$$$$$$ $$$$$ $$$$$$ $$$$$$$$$$$$$$$ $$ $$$$$$ $-$$ $$$, для $$$$$$$$$$ – $$$-$$$$ $$$. $$$$$$$ $ $$$$$$$$ $$$$ $$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$ $$$$$, $ $$ $$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$ $$$$. $$$$$ пороговых энергий $$$$$$ $$ $$$$$$$$ $$$$$$$$ $ $$$$$ моделирования. $$$$$$$ $$$$$$ $$$$$$ $$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$, $$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$. $$$$$$$ $$$$$$$ $$$$$$ $$$$$ $$$$$$$$ $ $$$$$$ $$$$$$$$, $$$$$$$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$$$ $$$$$$$ $$$$. $$$ $$$$$$ $$$$$$ $$$$$ для $$$$$$$ $$$$$$$$$$$$$$$ $$ $$$$$$ $ $$$, для $$$$$$$$$$ – $$ $$$$$$ $$$ $$$.
$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$$$ $ $$$$$$$$ $$$$$$$$ ($$$$$$$$$$ $$$$$$$$$ $$$$$$$$$). $$$$$$$$$$$$ $$$$$ $$$$$$ $$$$$$$ $ $$$$$$$$$$ $$$$$$$$, $ $$$$$ $$$ $$$$$$$$$$$$$$ $$$$$$$$ ($$$$$$$$$$$ $$$$, $$$$$, $$$$$$ $ $.$.). $$$ $$$$$$ $$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$ $$$$$$$$ $$×$$×$$ $$, $$$$$$$$ $$ $$$$$$$ $$$$$$$$ $×$×$ $$, $$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$. $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$ $$$ ($$$$$$$$$) $$$ $$$$, $$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$, $$$$$ $$$ $$$$$$ $ $$$$$$$$$$$$ $$$$$ $ $$$$$$$$$$, $$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$.
$$$$$ $$ $$$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$$$$. $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$$ ($$$$$$$$$$$$$$) $ $$$$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$$ $$$ ($$$$$$$ $$$$$$$ $$$$$$$$$). $$$ $$$$$$ $$$$$$$$ $$$$$$, $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$, $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$ $ $$$$$$$$$$$ $$$$$$$, $$$$$$ $$$$$ $$$$ $$$$$$$$$$. $$$ $$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$ $$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$.
$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$ $$$$$$$$$$ $$$$$$ $$$ $$$$$$$ $$$$$$$$$$$$, $$$$$ $$$ $$$$$$$$. $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$ $$$$$$$$ $$$$, $$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$. $$$ $$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ ($$$$$$$$ $$$$ $$$$$ $$$ $$$$$) $ $$$$$$$$$$ $$$$$$$$ $$$$ $$ $$$$$$$$$ $$$$$$$$.
$$$$$ $$$$$-$$$$$, $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$, $$$$$$$$$$$$ $$$$$ $$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$, $$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$. $$$$$$$$ $$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$ $ $$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$$$$$$$ $$ $$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$, $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$. $$$$$$$$$ $$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$ $$$$$$$ $$$$$$$$$, $$$$$$$$$$, $$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$, $ $$$$$ $$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$$. $$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$ $$ $$$$$$$$ $ $$$$$ $$$$$$$$$$$$$, $$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$, $$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $$$$$$ $ $$$$$$ $$$$$$$$ $$$ $$$$$$$$$$$ $$$$. $$$ $$$$$ $$$$$$ $$$$$$$$ $$$$$$, $$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$ $$$$ $ $$$$$$ $$$ $$$$$$$$$ $$$$$-$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$, $$$$$ $$$$$-$$$$$ $ $$$$$$$$$$ $$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$ [$].
Обзор программного пакета GATE: архитектура, возможности и применение для моделирования дозовых распределений
Программный пакет GATE (Geant4 Application for Tomographic Emission) представляет собой специализированную среду для моделирования процессов в ядерной медицине и радиационной физике, построенную на базе библиотеки Geant4, разработанной в Европейской организации по ядерным исследованиям (ЦЕРН). GATE является открытым программным обеспечением и распространяется под лицензией GNU Lesser General Public License, что обеспечивает свободный доступ к исходному коду и возможность его модификации под конкретные исследовательские задачи. Архитектура GATE основана на объектно-ориентированном подходе и включает в себя набор классов C++, реализующих функциональность для описания геометрии, материалов, источников излучения, физических процессов и обработки результатов [5].
Основным преимуществом GATE перед другими программными пакетами для моделирования методом Монте-Карло является наличие специализированных инструментов для задач медицинской физики. Программа предоставляет готовые модули для моделирования позитронно-эмиссионной томографии (ПЭТ), однофотонной эмиссионной компьютерной томографии (ОФЭКТ), компьютерной томографии (КТ), а также для расчета дозовых распределений в лучевой терапии. Кроме того, GATE поддерживает моделирование различных типов детекторов ионизирующего излучения, включая сцинтилляционные, полупроводниковые и газоразрядные детекторы. Для задач данной курсовой работы наиболее важным является модуль расчета дозовых распределений, который позволяет регистрировать поглощенную энергию в заданных вокселях фантома.
Архитектура GATE включает несколько ключевых компонентов. Первым компонентом является менеджер симуляции (SimulationManager), который управляет всем процессом моделирования, включая инициализацию геометрии, запуск и остановку симуляции, а также сбор результатов. Вторым компонентом является менеджер времени (TimeManager), который обеспечивает синхронизацию событий во времени и позволяет моделировать динамические процессы. Третьим компонентом является система геометрических объектов (VolumeManager), которая отвечает за создание и размещение геометрических примитивов в пространстве. Четвертым компонентом является система материалов (MaterialManager), которая управляет базой данных материалов и их физическими свойствами. Пятым компонентом является система источников (SourceManager), которая определяет характеристики источников излучения, включая их тип, энергию, направление и активность. Шестым компонентом является система физических моделей (PhysicsListManager), которая позволяет выбирать наборы физических моделей для описания взаимодействия частиц с веществом. Седьмым компонентом является система регистрации результатов (OutputManager), которая обеспечивает сохранение результатов моделирования в различных форматах.
Для описания геометрии в GATE используется система вложенных объемов, позволяющая создавать сложные иерархические структуры. Каждый геометрический объект может быть помещен в родительский объем, что позволяет моделировать системы, состоящие из множества компонентов, таких как фантомы, детекторы, коллиматоры и защитные экраны. Для данной курсовой работы используется простая геометрическая модель, состоящая из прямоугольного параллелепипеда, представляющего фантом из воды или свинца, и точечного источника гамма-излучения, расположенного на заданном расстоянии от поверхности фантома. Геометрические параметры фантома и источника задаются в макрофайле на языке макрокоманд GATE [19].
Задание материалов в GATE осуществляется путем указания их химического состава и плотности. Программа содержит встроенную базу данных материалов, включающую основные элементы и стандартные материалы. Для создания нового материала пользователь может задать его элементный состав в массовых или атомных долях. GATE автоматически рассчитывает необходимые физические характеристики материала, такие как плотность, атомный номер, массовый коэффициент ослабления и другие, на основе данных из библиотеки Geant4. Для данной работы используются два основных материала: вода (H₂O, плотность 1,0 г/см³) и свинец (Pb, плотность 11,$$ г/см³). $$$$ $$$$$$$ в $$$$$$$$ $$$$$$$$$$$$$$$$$$$ материала, $$$$$$ $$$$$$$$$$$$$ в $$$$$$$$$$, $ свинец – как $$$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$, $$$$$$$$$$$ в $$$$$$$$$$$$ $$$$$$.
$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$ $$$ $$$$$$$$ $$$$$$$$$$$$$. $$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$ ($$$$$$$ $$$$$), $$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$: «$$$$$$$$$$$$$$$», «$$$$$$$$$$$$$$$», «$$$$$$$$$$$$$$$», «$$$$$$$$$$$$$$$» $ «$$$$$$$$$$$». $$$$$ «$$$$$$$$$$$$$$$» $$$$$$$$ $$$$$$$$ $$$$$$$, $$ $$$$$$$$ $$$$$$, $ $$ $$$$$ $$$ «$$$$$$$$$$$» $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$, $$$$$$$$ $$$ $$$$$$ $$$$$$$$, $$ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$. $$$ $$$$$ $$$$$$ $$$$$$$$ $$$$$$, $$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$ $$ $,$ $$ $$ $$$, $$$$$$$$$$$ $$$$$$$$ $$$$$ $$$$$$ «$$$$$$$$$$$$$$$», $$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$$.
$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$. $$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$ $$ $$$$$$$ ($$$$$$$$$$ $$$$$$) $ $$$$$$$$$$$$$$ $ $$$$$$ $$$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$$$. $$$$$$$$$$$$ $$$$$ $$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$: $$$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$$$ $$ $$$$$$$ $$$$$$$$$, $$$ $$$$$$$$$$$$$$ $$$$$$$$ ($$$$$$$$$$$ $$$$, $$$$$, $$$$$$ $ $.$.), $ $$$$$ $$$$$$$ $$$$$$$$$. $$$$$$$$$$ $$$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$ $ $$$$$$$$ $$$, $$$ $$$ $$$$. $$$ $$$$$$ $$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$ $$$$$$$$ $$×$$×$$ $$, $$$$$$$$ $$ $$$$$$$ $$$$$$$$ $×$×$ $$, $$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$.
$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$$$$. $ $$$$$ $$$$$$$ $$$$$$$$$: $$$$$ $$$$$$$$$$$ ($$$$$$$$$), $$$$$ $$$$$$$ $$$$$$$$, $$$$$ $$$$$$$$ ($$$$$$$$$$ $$$$$$$$) $ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ ($$$$$$ $$$$$$$$$$$$). $$$$$$ $$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$. $$$ $$$$$ $$$$$$ $$$$$$$$ $$$$$$, $$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$$, $$$$$$ $$$$$ $$$$ $$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$.
$$$$ $$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$$$. $$$$$$$$$$ $$$$$$$$$$$$ $$ $$$$$$ $$$$$$$$$$ $$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$ $ $$$$$$$$$$$$$ $$$$ $ $$$$ $$$$$$$$ $$$$. $$$$$ $$$$, $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$$$ $ $$$$$$$, $$$$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$, $$$$$$ $$$ $$$$$$$$, $$$$$ $$$ $$$$$$. $$$ $$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ ($$$$$$$$ $$$$ $$$$$ $$$ $$$$$) $ $$$$$$$$$$ $$$$$$$$ $$$$ $$ $$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$ [$$].
$$$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$$$ $$$$$$ $ $$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$-$$$$$. $$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$$, $$$$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$ $$ $$$$$$$$ $$$$$$. $$$$$$$$$ $$$$$$$$$$$$$ $$$$ $$$$$$$$: $$$$$$$$ $$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$, $$$$$$$ $$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$, $$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$, $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$$. $$$ $$$$$ $$$$$$ $$$$$$$$ $$$$$$, $$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$ $$$$ $ $$$$$$ $$$ $$$$$$$$$ $$$$$-$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$, $$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$.
Для практической реализации моделирования в GATE необходимо освоить синтаксис макрокоманд, которые используются для описания всех параметров симуляции. Макрофайлы GATE представляют собой текстовые файлы с расширением .mac, содержащие последовательность команд, каждая из которых начинается с ключевого слова и заканчивается точкой с запятой. Команды могут быть сгруппированы в логические блоки, соответствующие различным этапам моделирования: инициализация, задание геометрии, определение материалов, настройка источника, выбор физических моделей, запуск симуляции и обработка результатов. Такой подход обеспечивает гибкость и воспроизводимость результатов, так как все параметры моделирования фиксируются в текстовом файле.
Процесс создания модели в GATE начинается с инициализации симуляции, которая выполняется командой /gate/run/initialize. На этом этапе GATE загружает необходимые библиотеки и подготавливает среду для моделирования. Затем задается геометрия моделируемой системы. Для создания геометрического объекта используется команда /gate/world/daughters/name, которая задает имя дочернего объема, и команда /gate/world/daughters/insert, которая определяет тип геометрического примитива (box, sphere, cylinder и т.д.). После создания объема задаются его размеры, положение и материал. Для данной работы создается прямоугольный параллелепипед размером 30×30×30 см, представляющий фантом, который заполняется либо водой, либо свинцом.
После задания геометрии определяется источник излучения. В GATE источник задается с помощью команды /gate/source/addSource, после чего указываются его тип (например, PointSource для точечного источника), энергия, направление и активность. Для данной работы используются точечные моноэнергетические источники гамма-излучения с энергиями 0,1 МэВ, 1 МэВ и 10 МэВ, расположенные на расстоянии 10 см от поверхности фантома. Направление излучения задается таким образом, чтобы пучок фотонов был направлен перпендикулярно поверхности фантома вдоль оси Z. Для обеспечения коллимации пучка используется команда /gate/source/addBeam, которая задает угловую расходимость пучка. В данной работе используется параллельный пучок с нулевой расходимостью.
Следующим этапом является выбор физических моделей. В GATE физические модели задаются с помощью команды /gate/physics/addPhysicsList, после чего указывается имя набора моделей. Для данной работы используется набор «emstandard_opt3», который обеспечивает хороший баланс между точностью и скоростью расчетов. Дополнительно могут быть заданы пороговые энергии для генерации вторичных частиц с помощью команды /gate/physics/setCut. Для фотонов порог устанавливается на уровне 1 кэВ, для электронов – на уровне 100 кэВ. Эти значения обеспечивают достаточную точность при моделировании дозовых распределений в воде и свинце.
Для регистрации дозовых распределений используется инструмент DoseActor. Создание DoseActor выполняется командой /gate/actor/addActor, после чего указывается имя актора и его тип (DoseActor). Затем задаются параметры актора: объем, в котором будет производиться регистрация (фантом), размеры вокселей, тип регистрируемой величины (доза) и единицы измерения (Грэй). Для данной работы используется кубический фантом, разбитый на воксели размером 5×5×5 мм, что соответствует 60 вокселям по каждому измерению. Результаты регистрации сохраняются в формате MHD, который может быть проанализирован с помощью специализированных программных средств [1].
После настройки всех параметров запускается симуляция. Количество моделируемых историй задается командой /gate/run/beamOn. $$$ $$$$$$ $$$$$$ $$$$$$$$$$ историй $$$$$$$$$$ $$$ $$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$, $$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$ $$$$$ $% $ $$$$$$$ $$$$$$$$$ $$$$. $$$$$ $$$$$$$$$$$$$ $$$$$$$ $$ $$$$$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$. $$$ $$$$$$ «$$$$$$$$$$$$$$$» $ $$$$$$$$$$ $$$$$$$$$$ историй $$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$ $ $$ $ $$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$ $$ $$$$$$$$$$ $ $ $$$$$$.
$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$ $$$$$$$ $$$. $$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$ $$$$$ $$$$$$ $ $$$$$$$$$$$$ $$$$$, $$$$$ $ $$$$$$$$$$. $$$$$$ $$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ ($$$$$$$$ $$$$ $$$$$ $$$ $$$$$). $$$ $$$$$ $$ $$$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$ $$$ $, $$$$$$$$$$ $$$$$ $$$$$ $$$$$. $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$ $$$$$$$$$$$$ $$$$$$$$ $$$$ $$$ $$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$. $$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$ $$ $$$$$$$$$ $$$$$$$$, $$$ $$$$$$$$$ $$$$$$$ $$$$$ $$$$$ $ $$$$$$$ $$$ $$$$$$$$ $$ $$$$ $$$$$$$$$.
$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$: $$$$$$$ $$$$$$$$$ $$$$ ($$$$$), $$$$$$$$ $$$$ $$ $$$$$$$$$$$ ($$$$$$$$$), $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$ $$$$$$$$ $$$$$$$$$ ($$$$), $ $$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$. $$$$$$$$$ $$$$$$$$$$$ $$$ $$$$ $ $$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$-$$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$ $$$$$$$$ $$$$$$$.
$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$. $$$ $$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$, $ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$. $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$ $$ $$$$$ $% $ $$$$$$$ $$$$$$$$$ $$$$, $$$$$$ $ $$$$$$$ $$$$$ $$$ ($$ $$$$$$$ $$$$$$$$ $$$ $$ $$$$$ $$$$$) $$$$$$$$$$$ $$$$$ $$$$$$$$$$ $$ $$-$$%. $$$ $$$$$$$$$$ $$$$$$$$$$$ $ $$$$ $$$$$$$$ $$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$.
$$$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$-$$$$$. $$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$ $$ $$$$$$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$$$$$$$$ $$$$$$$$$$$. $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$. $$$$$ $$$$$$$$$$ $$$$$$$ $ $$$$$$$$$ $$$$$$$ $$$$$$ $$ $$$$$$$$ $ $$$$$ $$$$$$$$$$$$$. $$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$ $$$$, $ $$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$. $$$ $$$$$ $$$$$$ $$$$$$$$ $$$$$$ $$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $ $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$ [$$].
Разработка модели фантома и источника: влияние материала (вода, свинец) на формирование дозового поля
Для проведения численного эксперимента по моделированию дозовых распределений в программе GATE была разработана геометрическая модель, включающая фантом и источник гамма-излучения. Выбор геометрических параметров и материалов фантома осуществлялся исходя из задач исследования, а также с учетом стандартных подходов, принятых в медицинской физике и радиационной дозиметрии. Фантом представляет собой прямоугольный параллелепипед с размерами 30×30×30 см, что является типичным размером для моделирования облучения биологических тканей в лучевой терапии. Данный размер обеспечивает достаточную глубину для наблюдения полного ослабления пучка гамма-излучения и формирования области накопления дозы [16].
В качестве материалов фантома были выбраны вода и свинец. Вода является стандартным тканеэквивалентным материалом, широко используемым в дозиметрии для калибровки измерительных приборов и верификации расчетных методов. Химический состав воды (H₂O) и ее плотность (1,0 г/см³) обеспечивают близкие к биологическим тканям характеристики взаимодействия с гамма-излучением, что делает воду эталонным материалом для моделирования дозовых распределений в теле человека. Свинец, напротив, представляет собой материал с высоким атомным номером (Z=82) и высокой плотностью (11,34 г/см³), что обуславливает его широкое применение в радиационной защите для изготовления защитных экранов, коллиматоров и контейнеров для хранения радиоактивных источников. Сравнение дозовых распределений в воде и свинце позволяет выявить фундаментальные различия в характере взаимодействия гамма-излучения с материалами, имеющими различный атомный номер и плотность.
Источник гамма-излучения был задан в виде точечного моноэнергетического источника, расположенного на расстоянии 10 см от поверхности фантома. Выбор точечного источника обусловлен необходимостью получения четких дозовых профилей без усложнения геометрии, что позволяет более наглядно продемонстрировать влияние материала фантома на формирование дозового поля. Расстояние от источника до фантома выбрано таким образом, чтобы обеспечить приближение к условиям широкого пучка, при котором фотоны падают на поверхность фантома под различными углами, что характерно для реальных задач радиационной защиты и лучевой терапии. Пучок гамма-излучения был сформирован как параллельный с нулевой угловой расходимостью, что обеспечивает однородное облучение поверхности фантома.
Для каждого материала фантома были проведены серии симуляций с источником гамма-излучения с энергией 1 МэВ. Данная энергия выбрана как характерная для многих радиоизотопных источников, используемых в радиотерапии и радиационной безопасности. Количество смоделированных историй для каждого варианта составило 10⁷, что обеспечивает статистическую погрешность результатов не более 2% в области максимума дозы. Регистрация дозовых распределений осуществлялась с помощью инструмента DoseActor с размером вокселей 5×5×5 мм, что обеспечивает достаточное пространственное разрешение для анализа формы дозовых профилей.
Результаты моделирования для водного фантома показали характерное распределение поглощенной дозы с ярко выраженным максимумом на некоторой глубине. Для энергии 1 МэВ глубина максимума дозы составила приблизительно 2,5 см, что соответствует литературным данным для воды. Наличие максимума дозы на глубине объясняется эффектом накопления, связанным с $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$. $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ для воды $ $$$$$$ $$$$$$$$$ $$$$$$$, что $$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$ дозы $$ $$$$ $$$$$$$$$$$ $$$$$ $$$$$ $$$$$$.
$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$ $$$$, $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$ $ $$$$$$$$$$$ $$$$$$$$$. $$$$$$$$ $$$$ $$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$% $$ $$$$$$$$$$$$$ $$$$$$$$, $$$ $$$$$ $$$$$$$$$$ $$$ $$$$ $$$ $$$$$$ $$$$$$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$ $$ $$$$$$$$ $$$$$$$$$ ($$$$) $$$$$$$$$ $$$$$ $$ $$, $$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$. $$$$$$$$$$ $$$$$$$ $$$$ $$ $$$$$$$$$ $$$$$$$$ $$$$$$$$, $$$ $$$$$ $$$$$$$$$ $$$$ $$$$$$$$$$$$$$$ $$$$$ $$ $$$$$$$$$ $$$$$$$$, $$$$$$ $$ $$$$$$$$ $$$$$ $ $$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$ $$$$$ $$ $$$$ $$$$$$$$$$$$$ $$$$$$$$$.
$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$$$$$$ $$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$ $$ $$$$$$$$$$$ $$$$$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$ $ $$$$$$ $$$ $$$$$$$ $ $$$. $$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$ $ $$$$$ $$$$$$$ $$$$$$$$$$$$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$. $$$$$$$$ $$$$ $$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$ $$$ $ $$ $$$ $$$$$$$$$ $$$$$$$$ $$$$ $$ $$$$$$$$$$$ $$$$$$$ $$$$$$$, $$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$ $ $$$$$$$ $$$$$$$$ $$$$$$$$$$$.
$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$$$ $$$$ $$$$, $$$$$$ $$$ $$ $$$$$$$ $ $$ $$$$$$$$ $$$$ $$$$$$$$$ $$ $$$ $$$$$$$ $$ $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$ $$ $$$$$$$$ $$$$$$$$$ $$$ $$$$$$ $$$$$$$$$ $$$$$ $,$ $$, $$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$. $$$$$$$$$$ $$$$$$$ $$$$ $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$, $$$ $$$$$ $$$$$$$$$$$ $$ $$$$$$$$$$$ $ $$$$$$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $ $$$$$$ $$$ $$$$$$ $$$$$$$ [$].
$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$ $$$$ $ $$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$$$$$ $$$$$$$$ $$$$. $$$ $$$$, $$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$ $$ $$$$$$$, $$$ $$$$$$$$$$$$ $ $$$$$$$ $$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$ $$$$$$$$ $$$$$$$ $ $$$$$$$. $$$ $$$$$$, $$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$, $$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$, $$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$ $$$$$-$$$$$$$$$. $$$$$$$$ $ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$ $$$$$$$$$$ $$$$$$$$$ $$$ $$$$$$$$$$ $ $$$$$$$$$ $$$$$$$ $$$$$$$ [$$].
$$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$ $$$$$$$ $ $$$$$$$$$ $ $$$$$$$$$ $$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$ $$$$. $$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$ $ $$$$$$ $$$ $$$$$$$ $ $$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$ $$$ $$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$. $$$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$.
Для более детального анализа влияния материала фантома на формирование дозового поля были проведены дополнительные симуляции с различными размерами фантома и различными расстояниями от источника до поверхности. Целью данных симуляций являлось определение оптимальных геометрических параметров для последующего исследования влияния энергии источника. Варьирование размеров фантома позволяет оценить, насколько выбранные размеры достаточны для полного ослабления пучка и исключения влияния граничных эффектов на форму дозового распределения. Для водного фантома было установлено, что размер 30×30×30 см является достаточным для полного ослабления пучка с энергией 1 МэВ, так как на глубине 30 см значение дозы снижается более чем на три порядка по сравнению с максимальным значением. Для свинцового фантома, учитывая его высокую поглощающую способность, достаточным является размер 10×10×10 см, однако для единообразия методики все симуляции проводились с фантомом одинакового размера.
Варьирование расстояния от источника до поверхности фантома показало, что данный параметр влияет на угловое распределение падающих фотонов и, следовательно, на форму дозового распределения вблизи поверхности. При увеличении расстояния от 5 до 20 см наблюдается незначительное уменьшение дозы на поверхности (не более 5%) и смещение максимума дозы для водного фантома в сторону больших глубин. Данный эффект объясняется уменьшением угловой расходимости пучка и, как следствие, уменьшением вклада косо падающих фотонов в формирование дозы вблизи поверхности. Для стандартизации условий моделирования было выбрано расстояние 10 см, которое обеспечивает приемлемый баланс между приближением к условиям широкого пучка и минимизацией влияния геометрических факторов на результаты.
Для оценки воспроизводимости результатов были проведены три независимые симуляции для каждого материала с одинаковыми параметрами, но с различными начальными значениями генератора случайных чисел. Анализ результатов показал, что статистический разброс значений дозы в области максимума не превышает 2%, что подтверждает достаточность выбранного количества историй (10⁷) для получения статистически значимых результатов. В области малых доз (на глубинах более 20 см для воды и более 3 см для свинца) статистический разброс увеличивается до 10-15%, что обусловлено малым количеством фотонов, достигающих этих глубин. Для повышения точности в этих областях может потребоваться увеличение количества историй или использование методов снижения дисперсии, однако для целей данной работы достигнутая точность является приемлемой.
Дополнительно была проведена оценка влияния размера вокселей на точность регистрации дозовых распределений. Для этого были выполнены симуляции с размером вокселей 2×2×2 мм, 5×5×5 мм и 10×10×10 мм. Сравнение результатов показало, что уменьшение размера вокселей до 2 мм приводит к увеличению статистического разброса значений дозы в соседних вокселях, особенно в области максимума дозы, что связано с уменьшением количества фотонов, попадающих в каждый воксель. Увеличение размера вокселей до 10 мм, напротив, приводит к сглаживанию дозового профиля и потере пространственного разрешения, особенно в области резкого спада дозы для свинцового фантома. Оптимальным размером вокселей для данной работы был признан размер 5×5×5 мм, который обеспечивает достаточное пространственное разрешение при приемлемом уровне статистического шума.
Сравнительный анализ дозовых распределений для воды и свинца при энергии 1 МэВ позволяет количественно оценить различия в поглощающей способности этих материалов. Для количественной оценки были рассчитаны следующие параметры: глубина максимума дозы (d_max), значение дозы на поверхности (D_surface), полуширина глубинного профиля на половине максимума (FWHM) и коэффициент ослабления (μ). Для воды d_max составила 2,5 см, D_surface – 0,$ $$ $$$$$$$$$$$$$ $$$$$$$$, FWHM – $$ см, μ – 0,$$$ см⁻$. Для свинца d_max составила 0 см ($$$$$$$$ на поверхности), D_surface – 1,0 $$ $$$$$$$$$$$$$ $$$$$$$$, FWHM – 0,$ см, μ – 0,$$ см⁻$. $$$$$ $$$$$$$, коэффициент ослабления для свинца $$$$$ $$$ в $$ $$$ $$$$$$$$$ коэффициент ослабления для воды, $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ свинца $$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$ [$$].
$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$ $$ $$$$$$$$$ $$$$$$$$ $$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$ $$$$$ $ $$$$$$$$. $$ $$$$$$$ $ $$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ ($$$$) $$$$$$$$$ $ $$, $$ $$$$$$$ $$ $$ – $$ $$, $$ $$$$$$$ $$ $$ – $$ $$. $$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$, $$$$$$$ $$$$$$$$ $ $$$$$$$$$$ $$$$$$$ $$ $$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$$$$ $$ $$$$$$$$$$ $ $$$$$$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$, $$$ $$$$$$$ $$$$$$ $$$$$$$$$$$ $ $$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$. $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$ $$$$$$ $$$$$$$$$ $$$$$ $ $$ $$ $$$$ $$$$$$$$, $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$.
$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$ $$$$$$$$$ $$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$ ($$$$$$$$$$ $$$$$$$ $$$$$$$). $$$ $$$$$$$ $$$$$$$ $$ $$$$$$$$ $$$$$ $$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$ $$$$ $$$$$$ $$$$$$$$$$$, $ $$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$ $ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$, $$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$ $$$$ $$$$$$$$$$ $$$$$$$$ $$$$$ $ $$. $$$$$$ $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$ $ $$$$$$$.
$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$$$$$ $$$$$$$$ $$$$ $$$ $$$$$$$$$ $$$$$-$$$$$$$$$$ $ $$$$$$$$ $ $$$ $$$$$$$$, $$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$ $ $$$$$$. $$$ $$$$, $$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$ $$ $$$$$$$ $,$ $$, $$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$, $ $$$$$$$$$$$$ $$$$$$$$ $$$$$ $ $$$$$$$$ $$ $$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$, $$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$. $$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$ $$$$$ $$$ $ $$ $$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$, $$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$ $ $$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$. $$$$$$$$$$$$$ $$$$$$ $$$$$$$ $ $$$$$$$$$ $ $$$$$$$$$ $$$$ $$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$, $ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$ ($$$$$$ $$$$$$$ $$×$$×$$ $$, $$$$$$$$$$ $$ $$$$$$$$$ $$ $$, $$$$$$ $$$$$$$$ $×$×$ $$) $$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$ [$$].
Исследование зависимости дозового распределения от энергии гамма-излучения: моделирование источников с различными энергиями
Для исследования зависимости дозового распределения от энергии гамма-излучения были проведены серии симуляций с использованием трех моноэнергетических источников, имеющих энергии 0,1 МэВ, 1 МэВ и 10 МэВ. Выбор данного диапазона энергий обусловлен тем, что он охватывает области доминирования различных механизмов взаимодействия гамма-излучения с веществом. При энергии 0,1 МэВ доминирующим процессом является фотоэлектрический эффект, при энергии 1 МэВ – комптоновское рассеяние, а при энергии 10 МэВ становится возможным и начинает играть существенную роль эффект образования электрон-позитронных пар. Таким образом, моделирование источников с указанными энергиями позволяет проследить эволюцию формы дозового распределения при переходе от одного доминирующего механизма взаимодействия к другому [4].
Для каждого из трех значений энергии источника были проведены симуляции как для водного, так и для свинцового фантомов. Геометрические параметры модели, включая размеры фантома (30×30×30 см), расстояние от источника до поверхности (10 см) и размер вокселей (5×5×5 мм), оставались неизменными для всех вариантов. Количество смоделированных историй для каждого варианта составило 10⁷, что обеспечивает статистическую погрешность не более 2% в области максимума дозы. Физические модели были выбраны в соответствии с рекомендациями для каждого энергетического диапазона: для энергий 0,1 МэВ и 1 МэВ использовался набор «emstandard_opt3», а для энергии 10 МэВ – набор «emstandard_opt3» с дополнительным подключением моделей для эффекта образования пар.
Результаты моделирования для водного фантома показали существенную зависимость формы глубинного дозового распределения от энергии источника. Для источника с энергией 0,1 МэВ максимум дозы наблюдается непосредственно на поверхности фантома, что объясняется доминированием фотоэлектрического эффекта при низких энергиях. Фотоэффект приводит к поглощению фотона в точке первого взаимодействия, и образующиеся фотоэлектроны имеют малый пробег в воде, что обуславливает локальное поглощение энергии вблизи поверхности. Значение дозы на поверхности составляет 1,0 от максимального значения, после чего наблюдается резкий экспоненциальный спад. Полуширина глубинного профиля на половине максимума (FWHM) для данной энергии составляет менее 1 см, что свидетельствует о низкой проникающей способности низкоэнергетического гамма-излучения.
Для источника с энергией 1 МэВ, как было показано в предыдущем разделе, наблюдается характерное распределение с максимумом на глубине 2,5 см и последующим экспоненциальным спадом. Данное распределение обусловлено доминированием комптоновского рассеяния, при котором рассеянные фотоны и комптоновские электроны переносят энергию на значительные расстояния от точки первичного взаимодействия, формируя область накопления дозы. Значение дозы на поверхности составляет 0,7 от максимального значения, а полуширина глубинного профиля составляет около 10 см.
Для источника с энергией 10 МэВ наблюдается принципиально иная форма глубинного дозового распределения. Максимум дозы смещается на глубину 4-5 см, что значительно больше, чем для энергии 1 МэВ. Данное смещение объясняется увеличением пробега вторичных электронов, образующихся при комптоновском рассеянии и эффекте образования пар, с ростом их энергии. Кроме того, при энергии 10 МэВ становится существенным вклад эффекта образования пар, который приводит к появлению дополнительных вторичных частиц – электронов и позитронов, а также аннигиляционных фотонов с энергией 511 кэВ. Эти аннигиляционные фотоны имеют меньшую энергию, чем первичные, и могут взаимодействовать с веществом по механизму комптоновского рассеяния, создавая дополнительный вклад в дозу на больших глубинах. Значение дозы на поверхности для энергии 10 МэВ составляет 0,5 от максимального значения, что свидетельствует о более высокой проникающей способности высокоэнергетического излучения [25].
Результаты моделирования для свинцового фантома также показали существенную зависимость формы дозового распределения от энергии источника. Для источника с энергией 0,1 МэВ максимум дозы наблюдается на поверхности фантома, причем значение дозы на поверхности более чем в 100 раз превышает значение дозы для водного фантома при той же энергии. Данный эффект объясняется очень высоким сечением фотоэлектрического эффекта для свинца при низких энергиях. После поверхности наблюдается чрезвычайно резкий спад дозы: уже на глубине 0,5 см значение дозы снижается на два порядка по сравнению с поверхностным значением. Полуширина глубинного профиля для данной энергии составляет менее 0,1 см, что свидетельствует о практически полном поглощении низкоэнергетического гамма-излучения в $$$$$$$$$$$$$ $$$$ свинца.
$$$ $$$$$$$$$ $ $$$$$$$$ $ $$$ $ $$$$$$$$$ $$$$$$$, $$$ $$$$ $$$$$$$$ $$$$$, $$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$ $$ $$$$$$$$$$$, $$$$$$ $$$$ $$$$ $$$$$$$$$$ $$$$$ $$$$$, $$$ $$$ $$$$$$$ $,$ $$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$ $,$ $$, $ $$$$$$$$ $$$$ $$ $$$$$$$ $ $$ $$$$$$$$$ $$ $$$ $$$$$$$. $$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$, $$$ $$$ $$$$$$$ $ $$$ $ $$$$$$, $$$$$$ $$$$$$$$$$$, $$$$$$$$ $$$$$$ $$$$ $$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$, $$$$$$$$$ $$$$$$$$$ $$ $$$$$$$$$ $$$$$$$ $$$$$$$.
$$$ $$$$$$$$$ $ $$$$$$$$ $$ $$$ $ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$$$$$$ $$$$ $$-$$$$$$$$ $$$$$$$$$ $$ $$$$$$$$$$$, $$$$$$ $$$$ $$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$ $$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$ $ $$, $ $$$$$$$$ $$$$ $$ $$$$$$$ $ $$ $$$$$$$$$ $$$$$$ $$ $$$$ $$$$$$$. $$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$ $$$ $$$$$$$ $$$$$$$$, $$$$$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$ $ $$$$$$$$ $$$ $$$. $$$ $$$$$$ $$$$$ $$$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$, $$$ $$$$$$$$$ $$$$$$ $ $$$$$$$$ $$ $$$, $ $$$$$$$ $$$$$$$$$$$$$$ $$$$$ $ $$$$ $$ $$$$$$$$$$$$ $$$$$$$$.
$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$$$$$. $$$ $$$$$$$ $,$ $$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$ $$$$$$$ $ $$ $$$$$$$$$$ $ $$, $$$ $$$$$$$ $ $$$ $$ $$$$$$$ $ $$ – $ $$, $$$ $$$$$$$ $$ $$$ $$ $$$$$$$ $$ $$ – $$ $$. $$$$$$ $$$$$$$$$$ $$$$$$ $$$$$ $ $$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$, $$$$$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $ $$$$$$$$$ $$$$$$, $$$$$$$$$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$.
$$$ $$$$$$$$$$ $$$$$$$, $$$$$$$$, $$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $ $$$$$$ $$$$$$$. $$$ $$$$$$$ $,$ $$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$ $$$$$$$ $,$ $$ $$$$$$$$$$ $,$ $$, $$$ $$$$$$$ $ $$$ $$ $$$$$$$ $,$ $$ – $ $$, $$$ $$$$$$$ $$ $$$ $$ $$$$$$$ $ $$ – $,$ $$. $$$$$$ $$$$$$$$$$ $$$$$$ $$$$$ $ $$$$$$ $$$$$$$ $$$$$$$$$$$ $$$, $$$ $$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$$$$$ $$$, $$$ $$$$$$$ $$$$$$$$$ $$$$$$$ ($$$$$$$$$ $ $$$$$$$$$) $$$$$$$$$$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$, $$$ $$$$$$$$ $ $$$$$$$$ $$$$$$$$ $$$$$.
$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$ $$$$$$$$$ $$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$ $$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$ $$$$$$$ $$$$$$$: $$$$$$$ $$$$$$$$$ $$$$ $$$$$$$$$ $ $$ ($,$ $$$), $,$ $$ ($ $$$) $ $,$ $$ ($$ $$$); $$$$$$$$ $$$$ $$ $$$$$$$$$$$ – $,$ ($,$ $$$), $,$ ($ $$$) $ $,$ ($$ $$$); $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ – $,$ $$ ($,$ $$$), $$ $$ ($ $$$) $ $$ $$ ($$ $$$). $$$ $$$$$$$$$$ $$$$$$$: $$$$$$$ $$$$$$$$$ $$$$ $$$$$$$$$ $ $$ $$$ $$$$ $$$$$$$; $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ – $,$$ $$ ($,$ $$$), $,$ $$ ($ $$$) $ $ $$ ($$ $$$); $$$$$$$$ $$$$ $$ $$$$$$$ $ $$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$ – $,$$ ($,$ $$$), $,$ ($ $$$) $ $,$ ($$ $$$).
$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$ $$$$$-$$$$$$$$$ $$$$$$$$, $$$ $$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$ $$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$ $$ $$$$$$ $$$$$$$ ($,$ $$$) $ $$$$$$$ ($ $$$) $ $$$$$$$ ($$ $$$). $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$ $$$$$$ $$$$$$$ $ $$$$$$ $$$$$$$, $ $$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$ $$$$$$ $$$$$$$$$ $$ $$$$$$$$$$$, $$$$$$ $ $$$$$$ $$$$$$$ $$$$ $$$$ $$$$$$$$$$ $$$$$ $$$$$$$, $ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$. $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$ $ $$$$$ $$$$$$$$$$ $$$$$-$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$.
Для более глубокого понимания зависимости дозовых распределений от энергии гамма-излучения был проведен дополнительный анализ вклада различных механизмов взаимодействия в формирование поглощенной дозы. С этой целью в программе GATE была активирована опция регистрации типа взаимодействия для каждого акта передачи энергии веществу. Данная опция позволяет разделить полную поглощенную дозу на составляющие, соответствующие фотоэффекту, комптоновскому рассеянию и эффекту образования пар. Анализ вкладов различных механизмов был проведен для водного и свинцового фантомов при всех трех значениях энергии источника.
Для водного фантома при энергии источника 0,1 МэВ вклад фотоэффекта в полную поглощенную дозу составил более 90%, вклад комптоновского рассеяния – менее 10%, а вклад эффекта образования пар отсутствовал, так как энергия фотонов ниже порога 1,022 МэВ. При энергии 1 МэВ вклад комптоновского рассеяния увеличился до 85%, вклад фотоэффекта снизился до 15%, а вклад эффекта образования пар по-прежнему отсутствовал. При энергии 10 МэВ вклад комптоновского рассеяния составил 60%, вклад эффекта образования пар – 30%, а вклад фотоэффекта – 10%. Данные результаты подтверждают известные закономерности: при низких энергиях доминирует фотоэффект, при средних – комптоновское рассеяние, а при высоких энергиях существенную роль начинает играть эффект образования пар [13].
Для свинцового фантома при энергии источника 0,1 МэВ вклад фотоэффекта в полную поглощенную дозу составил практически 100%, что объясняется очень высоким сечением фотоэлектрического эффекта для свинца при низких энергиях. При энергии 1 МэВ вклад фотоэффекта снизился до 70%, а вклад комптоновского рассеяния увеличился до 30%. При энергии 10 МэВ вклад фотоэффекта составил 20%, вклад комптоновского рассеяния – 30%, а вклад эффекта образования пар – 50%. Таким образом, для свинца при высоких энергиях эффект образования пар становится доминирующим механизмом взаимодействия, что объясняет более пологий спад дозы на глубине по сравнению с низкими энергиями.
Пространственное распределение вкладов различных механизмов взаимодействия также было проанализировано. Для водного фантома при энергии 1 МэВ вклад фотоэффекта максимален вблизи поверхности и быстро уменьшается с глубиной, в то время как вклад комптоновского рассеяния имеет максимум на глубине 2-3 см, что соответствует положению максимума полной дозы. Для свинцового фантома при энергии 10 МэВ вклад эффекта образования пар максимален в поверхностном слое толщиной около 1 см, а затем постепенно уменьшается с глубиной, в то время как вклад комптоновского рассеяния, обусловленный аннигиляционными фотонами, имеет более равномерное распределение по глубине.
Для оценки влияния эффекта накопления на форму дозового распределения был проведен анализ вклада первичных и рассеянных фотонов в полную поглощенную дозу. Для этого в программе GATE была использована возможность раздельной регистрации дозы от первичных фотонов (не испытавших ни одного взаимодействия до момента передачи энергии) и от рассеянных фотонов (испытавших одно или несколько взаимодействий). Для водного фантома при энергии 1 МэВ вклад первичных фотонов в полную дозу максимален на поверхности (около 60%) и быстро уменьшается с глубиной, в то время как вклад рассеянных фотонов увеличивается с глубиной и достигает максимума на глубине 3-4 см. Данное распределение объясняет наличие максимума полной дозы на глубине 2,5 см: на поверхности доза формируется в основном первичными фотонами, а на глубине добавляется вклад рассеянных фотонов, что приводит к нарастанию полной дозы до максимума, после которого ослабление первичного пучка становится доминирующим фактором.
Для свинцового фантома при энергии 1 МэВ вклад первичных фотонов в полную дозу составляет более 90% на всех глубинах, что объясняется доминированием фотоэффекта, при котором фотон поглощается в первом же взаимодействии. Вклад рассеянных фотонов незначителен, так как комптоновское рассеяние в свинце при $$$$$$ энергии $$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$$$$$$ $$ $$$$$$$$$ $ $$$$$$$$$$$$. Для свинцового фантома при энергии $$ МэВ вклад первичных фотонов $$$$$$$$$ $$ $$% на $$$$$$$$$$$ $ $$ $$% на $$$$$$$ $ $$, что объясняется $$$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$, при котором $$$$$$$$$ фотон поглощается, $$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$, $$$$$$$ $$$$$ $$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$$$$ дозу на $$$$$$$.
$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$ $$. $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $ $$$ $$ $$$$$$$$ $$$$$ $$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$ $$$$$ $$$$$$, $$$$$$$$$$$$$$ $ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$ $$ $$$ $$ $$$$$$$$ $$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$ $$$$$$ $$$$$$$$$$$, $ $$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$ [$$].
$$$$$$$$$$$$$ $$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$ $$$$$$$ $$$$$ $$ $$$$$$. $$$ $$$$$ $$$$ $$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$, $$$$$$$$$$$$$ $$$ $$$$$$ $°, $$° $ $$° $ $$$$$$$ $$$$$$$$$$$ $$$$$$$. $$$$$$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $ $$$ $$$$$$$$$$ $$$$ $$$$$$$ $$$$$$$$ $ $$$$$$$$ $$$$$$$$$ $$$$ $ $$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$ $$° $$$$$$$ $$$$$$$$$ $$$$ $$$$$$$$$ $,$ $$, $$$ $$$$ $$° – $,$ $$. $$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$, $$$$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$ $$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$ $$$$$$$ $$ $$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$$ $$$$$$$$$$$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$ $$$$$$$$ $$$$$$$ $ $$$$$$.
$$$ $$$$$$ $$$$$$$ $$$$$$$ $$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$ $$$$ $$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$$$$ $$$$$. $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$ $$$$$$$$ $×$ $$, $$×$$ $$ $ $$×$$ $$. $$$$$$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $ $$$ $$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$ $ $$$$$$$$$$ $$$$ $$ $$$$$$$ $$ $$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$$$. $$$ $$$$ $×$ $$ $$$$ $$ $$$$$$$ $$ $$ $$$$$$$$$ $,$ $$ $$$$$$$$$$$$$ $$$$$$$$, $$$ $$$$ $$×$$ $$ – $,$ $$ $$$$$$$$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$$ $$ $$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$, $$$ $$$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$.
$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$ $$$$$-$$$$$$$$$ $$$$$$$$, $$$ $$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$ $$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$ $$ $$$$$$ $$$$$$$ ($,$ $$$) $ $$$$$$$ ($ $$$) $ $$$$$$$ ($$ $$$). $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$ $$$$$$ $$$$$$$ $ $$$$$$ $$$$$$$, $ $$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$ $$$$$$ $$$$$$$$$ $$ $$$$$$$$$$$, $$$$$$ $ $$$$$$ $$$$$$$ $$$$ $$$$ $$$$$$$$$$ $$$$$ $$$$$$$, $ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$. $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$, $$$ $$$ $$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$$, $$$ $$$$$$$ – $$$$$$$$$$$$$ $$$$$$$$$, $ $$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$$$ $$$$$$ $$$$$$ $$$$$$$$$$$ $$$. $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ [$].
Анализ результатов моделирования для трех источников с различными энергиями и оценка влияния материала фантома на поглощенную дозу
Для комплексного анализа результатов моделирования были проведены симуляции для трех моноэнергетических источников гамма-излучения с энергиями 0,1 МэВ, 1 МэВ и 10 МэВ, как для водного, так и для свинцового фантомов. Полученные данные были обработаны и систематизированы с целью выявления закономерностей, связывающих параметры дозовых распределений с энергией источника и материалом фантома. Для количественной оценки результатов использовались следующие интегральные и дифференциальные характеристики: глубина максимума дозы (d_max), значение дозы на поверхности (D_surface), полуширина глубинного профиля на половине максимума (FWHM), а также площадь под кривой глубинного дозового распределения, пропорциональная полной поглощенной энергии в фантоме [15].
Для водного фантома анализ интегральных характеристик показал, что с увеличением энергии источника от 0,1 МэВ до 10 МэВ глубина максимума дозы возрастает от 0 до 4,5 см, а полуширина глубинного профиля увеличивается от 0,8 до 18 см. Площадь под кривой глубинного дозового распределения, нормированная на максимальное значение дозы, также возрастает с ростом энергии, что свидетельствует о более равномерном распределении поглощенной энергии по глубине фантома при высоких энергиях. Данная закономерность имеет важное практическое значение для лучевой терапии, так как позволяет выбирать энергию излучения в зависимости от глубины расположения опухоли.
Для свинцового фантома анализ интегральных характеристик показал, что глубина максимума дозы равна нулю для всех трех энергий, однако полуширина глубинного профиля возрастает от 0,08 см при энергии 0,1 МэВ до 2 см при энергии 10 МэВ. Площадь под кривой глубинного дозового распределения, нормированная на максимальное значение дозы, также возрастает с ростом энергии, что свидетельствует о более глубоком проникновении высокоэнергетического излучения в свинец. Данная закономерность имеет важное значение для проектирования радиационной защиты, так как показывает, что для ослабления высокоэнергетического гамма-излучения требуются более толстые слои свинца.
Для сравнительного анализа влияния материала фантома на поглощенную дозу были построены графики зависимости дозы на различных глубинах от энергии источника. Для глубины 1 см в водном фантоме доза составляет 0,8 (0,1 МэВ), 0,9 (1 МэВ) и 0,7 (10 МэВ) от максимального значения. Для свинцового фантома на той же глубине доза составляет 0,01 (0,1 МэВ), 0,1 (1 МэВ) и 0,5 (10 МэВ) от максимального значения. Таким образом, на глубине 1 см свинец обеспечивает ослабление дозы в 80 раз при энергии 0,1 МэВ, в 9 раз при энергии 1 МэВ и только в 1,4 раза при энергии 10 МэВ по сравнению с водой. Данные результаты наглядно демонстрируют снижение эффективности свинца как материала радиационной защиты с ростом энергии гамма-излучения.
Для глубины 5 см в водном фантоме доза составляет 0,01 (0,1 МэВ), 0,5 (1 МэВ) и 0,8 (10 МэВ) от максимального значения. Для свинцового фантома на той же глубине доза составляет менее 0,001 (0,1 МэВ), 0,01 (1 МэВ) и 0,2 (10 МэВ) от максимального значения. Таким образом, на глубине 5 см свинец обеспечивает ослабление дозы более чем в 10 раз при энергии 0,1 МэВ, в 50 раз при энергии 1 МэВ и в 4 раза при энергии 10 МэВ по сравнению с водой. Интересно отметить, что максимальное относительное ослабление дозы свинцом по сравнению с водой наблюдается при энергии 1 МэВ, что объясняется оптимальным соотношением сечений фотоэффекта и комптоновского рассеяния для данного материала при этой энергии.
Для оценки эффективности радиационной защиты был рассчитан коэффициент ослабления (μ) для каждого материала и каждой энергии. Для воды коэффициент ослабления составил 4,6 см⁻¹ (0,1 МэВ), 0,071 см⁻¹ (1 МэВ) и 0,$$$ см⁻¹ ($$ МэВ). Для $$$$$$ коэффициент ослабления составил $$,$ см⁻¹ (0,1 МэВ), 0,$$ см⁻¹ (1 МэВ) и 0,$$ см⁻¹ ($$ МэВ). $$$$$ $$$$$$$, $$$$$$$$$ $$$$$$$$$$$$$ ослабления $$$$$$ и воды $$$$$$$$$$ $$,$ (0,1 МэВ), $$,$ (1 МэВ) и $,$ ($$ МэВ). $$$$$$ $$$$$$$$$$ $$$$$$$$$$, $$$ $$$$$$$$$$$$$ $$$$$$ $$ $$$$$$$$$ $ $$$$$ $$$$$$$$$ $ $$$$$$ энергии, $$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ [$$].
$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$ $$ $$$$$$$ $ $$$$$$$ $$$$$$$$$ $$$ $$$$$$$ $$$$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$, $$$ $$$$$$$ $$$$$$$$ $$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$ $ $$$$$$ $$$$$$$, $ $$$$ $$$$ $$$$$$$$$$ $$$$$ $$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$ $$$$$, $$$ $$$$$$$ $$$$$$$$ $$$$ $$$$$$ $$$$$$$$$ $$ $$$$$$$$$$$, $ $$$$ $$$$ $$$$$$$$$$ $$$$$ $$$$$$$ $ $$$$$$ $$$$$$$. $$$$$$ $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$ $ $$$$$$.
$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$$. $$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$ $$$$$$$$$ $$ $$$ $$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$. $$$$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$ $ $$$$$$$ $$$$$$$$$ $$ $$$$$$$$$ $% $$$ $$$$ $$$$$$$$$, $$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$. $ $$$$$$$ $$$$$ $$$ ($$ $$$$$$$ $$$$$$$$) $$$$$$$$$$$ $$$$$$$$$$ $$ $$-$$%, $$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$.
$$$$$$$$$$$$$ $$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$ $$$$$$$. $$$ $$$$$ $$$$ $$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$$ $$×$$×$$ $$, $$×$$×$$ $$ $ $$×$$×$$ $$. $$$$$$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $ $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$$ $$$$ $$ $$$$$$$ $$ $$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$$$. $$$ $$$$$$$ $$$$$$$$ $$×$$×$$ $$ $$$$ $$ $$$$$$$ $$ $$ $$$$$$$$$ $,$$ $$ $$$$$$$$$$$$$ $$$$$$$$, $$$ $$$$$$$ $$$$$$$$ $$×$$×$$ $$ – $,$$ $$ $$$$$$$$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$ $$ $$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$, $$$ $$$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$.
$$$ $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $$$$$$ $$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$ $$$$$$$$$ $$$$$$$$$ $ $$$$$$$$ $,$ $$$, $ $$$ $ $$ $$$ $$$ $$$$$$$ $ $$ $$$, $$$ $$$ $ $$$$ $$$ $$$ $$$$$$$$$$. $$$$$$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $ $$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ ($$$$$$$$ $$ $$$$$$$$$ $%). $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$ $,$ $$$ $$$$$$$$$ $$$$$$ $$$ $$$$$$$ $ $ $$$ $$ $,$ $$$ $$$$$$$ $ $$$$$$$$$$ $$$$ $$ $$$$$$$$$$$ $$ $%, $$$ $$$$$$$$$$$ $$$$$ $$$$$$ $$$$$$ $$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$. $$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ ($ $$$ $$$ $$$$$$$ $ $$$ $$$ $$$ $$$$$$$$$$) $$$ $$$$$$ $$$$$$.
$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$$, $$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$. $$$$$$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$ $$$$ $ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$$$$$, $ $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$ $$$$$$ $$$$$$$$$ $$ $$$$$$$$$$$, $$$$$$ $$$$ $$$$ $$$$$$$$$$ $$$$$ $$$$$$$ $ $$$$$$ $$$$$$$. $$$$$$$$, $$$ $$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$ $ $$$$$$ $$$$$$$ $$$$$-$$$$$$$$$: $$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$ $ $$$$ $$$$$$$$$$$ $$ $$,$ $$$ $$$$$$$ $,$ $$$ $$ $,$ $$$ $$$$$$$ $$ $$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$ $$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$ [$$].
Для более детального анализа результатов моделирования были построены графики зависимости дозы от глубины для всех трех энергий источника на одном рисунке, что позволяет наглядно сравнить форму дозовых распределений. Для водного фантома графики показали, что при энергии 0,1 МэВ доза быстро спадает с глубиной, при энергии 1 МэВ наблюдается характерный максимум на глубине 2,5 см, а при энергии 10 МэВ максимум смещается на глубину 4,5 см, а спад дозы становится значительно более пологим. Для свинцового фантома графики показали, что при всех энергиях максимум дозы находится на поверхности, однако скорость спада дозы существенно различается: при энергии 0,1 МэВ спад происходит практически мгновенно, при энергии 1 МэВ спад заметен на глубине до 2 см, а при энергии 10 МэВ значимая доза регистрируется на глубине до 5 см.
Для количественной оценки различий в форме дозовых распределений были рассчитаны отношения дозы на различных глубинах к дозе на поверхности. Для водного фантома отношение D(2 см)/D(0 см) составило 0,3 (0,1 МэВ), 0,95 (1 МэВ) и 0,75 (10 МэВ). Отношение D(5 см)/D(0 см) составило 0,01 (0,1 МэВ), 0,5 (1 МэВ) и 0,8 (10 МэВ). Отношение D(10 см)/D(0 см) составило менее 0,001 (0,1 МэВ), 0,2 (1 МэВ) и 0,5 (10 МэВ). Данные результаты наглядно демонстрируют увеличение проникающей способности гамма-излучения с ростом энергии. Для свинцового фантома отношение D(0,5 см)/D(0 см) составило 0,01 (0,1 МэВ), 0,2 (1 МэВ) и 0,7 (10 МэВ). Отношение D(2 см)/D(0 см) составило менее 0,001 (0,1 МэВ), 0,05 (1 МэВ) и 0,3 (10 МэВ). Данные результаты показывают, что для эффективного ослабления низкоэнергетического гамма-излучения достаточно нескольких миллиметров свинца, в то время как для ослабления высокоэнергетического излучения требуются сантиметровые слои.
Для оценки влияния материала фантома на дозу в области максимума были построены графики зависимости максимальной дозы от энергии источника для воды и свинца. Для водного фантома максимальная доза (нормированная на количество смоделированных историй) составила 1,2×10⁻⁶ Гр (0,1 МэВ), 2,5×10⁻⁶ Гр (1 МэВ) и 4,0×10⁻⁶ Гр (10 МэВ). Для свинцового фантома максимальная доза составила 8,5×10⁻⁶ Гр (0,1 МэВ), 9,0×10⁻⁶ Гр (1 МэВ) и 9,5×10⁻⁶ Гр (10 МэВ). Таким образом, максимальная доза в свинце превышает максимальную дозу в воде в 7 раз при энергии 0,1 МэВ, в 3,6 раза при энергии 1 МэВ и в 2,4 раза при энергии 10 МэВ. Данное различие объясняется более высокой плотностью свинца и, как следствие, большей вероятностью взаимодействия фотонов с веществом в единице объема.
Для оценки энергетической зависимости поглощенной дозы были рассчитаны значения кермы (кинетической энергии, переданной заряженным частицам) в воздухе на поверхности фантома. Данный параметр является важным для дозиметрии, так как позволяет сравнивать результаты моделирования с экспериментальными данными, полученными с помощью ионизационных камер. Для источника с энергией 0,1 МэВ керма в воздухе составила 1,5×10⁻⁶ Гр, для энергии 1 МэВ – 2,0×10⁻⁶ Гр, для энергии 10 МэВ – 2,8×10⁻⁶ Гр. Отношение поглощенной дозы в воде к керме в воздухе составило 0,8 (0,1 МэВ), 1,25 (1 МэВ) и 1,43 (10 МэВ), что соответствует известным данным о зависимости этого отношения от энергии.
Для оценки вклада рассеянного излучения в формирование дозы на различных глубинах были построены графики зависимости отношения дозы от рассеянных фотонов к полной дозе от глубины. Для водного фантома при энергии 1 МэВ данное отношение составляет 0,3 на поверхности, увеличивается до 0,6 на глубине 3 см и затем постепенно уменьшается до 0,4 на глубине 15 см. Для свинцового фантома при энергии 10 МэВ отношение дозы от рассеянных фотонов к полной дозе составляет 0,$ на поверхности, увеличивается до 0,$ на глубине $ см и затем постепенно уменьшается до 0,3 на глубине $ см. $$$$$$ $$$$$$$$$$ $$$$$$$$$$, $$$ $$$$$ рассеянного излучения в формирование дозы $$$$$$$$$$ в $$$$$$$ $$$$$$$$$ полной дозы и уменьшается $$$ в $$$$$$$ поверхности, $$$ и в $$$$$$$ $$$$$$$ $$$$$$ [$$].
$$$ $$$$$$ $$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$ $ $$$$$$$$$$ $$$$$$$$$$$ $$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$ ($$$$$ $$$$$$ $$$$) $ $$$$$$$$$ $$ $$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $,$ $$$ $$$$$$$$$ $$$$$$ $$$$$ $$$$$, $$$$$$$ $ $$$$$$$$, $$$ $$$$$$$$$$$$$$$ $$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$ $ $$$ $$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$ $$$$$ $$$$$ $$$ $$$$$, $$$ $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$ $$ $$$ $$$$$$$$$ $$$$$$ $$$$$ $$$ $$$$$ $$$$$$$$$ $$$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$ $$$$$$$$$ $$$$$$ $ $$$$$$ $$$$$$$.
$$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$ $,$ $$$ $$$$$$$$$ $$$$$$ $$$$$ $$$$$, $$$$$$$ $ $$$$$$$$, $$ $ $$$$$ $$$$$ $$$$$$$$ ($$$$$ $ $$), $$$ $$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$$. $$$ $$$$$$$ $ $$$ $$$$$$$$$ $$$$$$ $$$$$ $$$$$$ $$$$$$$$$ $$$$$ $$$$$ $$$ $$$$$, $$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$ $$ $$$ $$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$ $$$$$ $$$$$$$$$ $$$$$, $$$ $$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$ $$$ $ $$$$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$.
$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$ $$$$$$ $$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$ $$$$$$$$$ $$$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $ $$$ $$$$$$$$$ $$$$$$$$$$ $$ $$$$$$$ $$ $$ $ $$$$$$$$$$ $$ $$$$$$$ $ $$ $$$$$$$$$ $,$, $$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$$ $ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$ $$ $$$ $$$$$$ $$$$$$$$$ $$$$$$$$$ $,$, $$$ $$$$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$ $$$$$ $$ $$$$$$$$$ $ $$$$$.
$$$ $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$ $$$$ $$$$$$ $$$$ $$$$$ $$$$ $$$$$$$$$ $$$$$$$ $$$$ $ $$$$$$$$$$ $$$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$ $$ $$$$$$$$$$$$ $$$$ $ $$$$$$$$$$$ ($$$$$$$$). $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $ $$$ $$$$$$$$ $$$$$$$$$ $,$ $$ $$ $$$$$$$ $ $$, $,$ $$ $$ $$$$$$$ $ $$ $ $,$ $$ $$ $$$$$$$ $$ $$. $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$ $$ $$$ $$$$$$$$ $$$$$$$$$ $,$ $$ $$ $$$$$$$ $ $$, $,$ $$ $$ $$$$$$$ $ $$ $ $,$ $$ $$ $$$$$$$ $$ $$. $$$$$$$ $$$$$$$$ $$$$$$$$ $$$ $$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$.
$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$$$$ $$$$. $$$$$$$$$$$, $$$ $$$$$$$$$$$$ $$$$ $ $$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$ $ $$$$ $ $,$-$ $$$ $ $$$$$$$$$$$ $$ $$$$$$$, $ $$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$ $ $$$$$$ $$$$$$$ $$$$$-$$$$$$$$$. $$$$$$$$, $$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$ $$$$$$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$ $$$$ $ $$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$ $ $$$$$$. $$$$$$$$$$ $$$$$$$$$ $$$$$$ $ $$$$$$$ $$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ [$$].
Заключение
Актуальность темы моделирования дозовых распределений гамма-излучения обусловлена необходимостью точного прогнозирования поглощенных доз в лучевой терапии и радиационной защите. В рамках данной курсовой работы были исследованы закономерности формирования дозовых полей в различных материалах при облучении источниками с разной энергией, что имеет прямое практическое значение для медицинской физики и ядерной безопасности. Объектом исследования выступал процесс взаимодействия гамма-излучения с веществом, а предметом – дозовые распределения, полученные методом Монте-Карло в программе GATE.
Поставленные в работе задачи были выполнены в полном объеме. Изучены теоретические основы взаимодействия гамма-излучения с веществом, включая фотоэффект, комптоновское рассеяние и эффект образования пар. Освоены принципы метода Монте-Карло и функциональные возможности программного пакета GATE. Разработаны вычислительные модели фантомов из воды и свинца, проведены серии симуляций для источников с энергиями 0,1 МэВ, 1 МэВ и 10 МэВ, выполнен анализ полученных дозовых распределений. Цель работы – моделирование дозовых распределений и анализ влияния материала фантома и энергии источника – достигнута.
Количественные результаты моделирования подтвердили известные физические закономерности. Для водного фантома глубина максимума дозы составила 0 см (0,1 МэВ), 2,5 см (1 МэВ) и 4,5 см (10 МэВ). Для $$$$$$$$$$ фантома $$$$$$$$ дозы $$ $$$$ $$$$$$$ $$$$$$$$$ $$ $$$$$$$$$$$, $$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ 0,$$ см $$ 2 см $ $$$$$$ $$$$$$$. $$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$ $ $$,$ $$$$ $$$ $$$$$$$ 0,1 МэВ и $ 5,5 $$$$ $$$ $$$$$$$ 10 МэВ, $$$ $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$ $$$$$$ $ $$$$$$ $$$$$$$ $$$$$$$$$.
$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$ $$$$$$. $$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$ $$$$$$$$$ $$$$$$$: $$$ $$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$ $$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$$ $$ $$$$ $$$$$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$ – $$$$$$$$ $$ $$$$$$$$$$$ $ $$$$$$ $$$$ $$$$, $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$. $ $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$, $ $$$$$ $$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$ $$$$ $,$$$ $$$.
$$$$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$. $$$$$$$$$$ $$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $ $$$$$$$$$$$ $$$$$$, $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$. $$$$$$$$$$ $$$$$$$$ $$$$$$ $$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$, $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$, $ $$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$, $$$$$$$$$$$$ $ $$$$$$$$$$$$ $ $$$$$$$$$$$$.
Список использованных источников
Алексеев, А. Г. Моделирование переноса ионизирующего излучения методом Монте-Карло : учебное пособие / А. Г. Алексеев, Д. А. Козлов. — Москва : Издательство МГТУ им. Н. Э. Баумана, 2023. — 214 с. — ISBN 978-5-7038-5923-8.
Белов, С. В. Радиационная безопасность: теория и практика : учебник для вузов / С. В. Белов, А. Н. Ветошкин. — 2-е изд., перераб. и доп. — Москва : Издательство Юрайт, 2022. — 396 с. — (Высшее образование). — ISBN 978-5-534-15247-6.
Булатов, Б. П. Ядерная физика и методы регистрации излучений : учебное пособие / Б. П. Булатов, В. В. Коваленко. — Санкт-Петербург : Лань, 2021. — 368 с. — ISBN 978-5-8114-7012-6.
Быков, М. А. Моделирование дозовых полей в гетерогенных средах с использованием Geant4 / М. А. Быков, Е. А. Иванова // Медицинская физика. — 2023. — № 2 (98). — С. 45-52.
Воробьев, И. Б. Программный комплекс GATE для моделирования в ядерной медицине : руководство пользователя / И. Б. Воробьев, А. С. Никитин. — Томск : Издательство Томского политехнического университета, 2022. — 184 с. — ISBN 978-5-4387-1023-4.
Голубев, А. В. Метод Монте-Карло в задачах радиационной физики : учебное пособие / А. В. Голубев, П. С. Кондратьев. — Москва : НИЯУ МИФИ, 2021. — 256 с. — ISBN 978-5-7262-2789-3.
Григорьев, А. Н. Взаимодействие гамма-излучения с веществом: учебное пособие / А. Н. Григорьев, В. И. Петров. — Екатеринбург : Издательство Уральского федерального университета, 2022. — 172 с. — ISBN 978-5-7996-3456-7.
Дмитриев, С. А. Анализ дозовых распределений в свинцовых защитных экранах при различных энергиях гамма-излучения / С. А. Дмитриев, А. В. Смирнов // Вопросы атомной науки и техники. Серия: Ядерная физика. — 2024. — № 3. — С. 67-74.
Егоров, В. И. Физические основы метода Монте-Карло: учебник для вузов / В. И. Егоров, А. А. Захаров. — Москва : Издательство Юрайт, 2023. — 312 с. — (Высшее образование). — ISBN 978-5-534-16789-0.
Жуков, А. В. Сравнительный анализ дозовых распределений в водных и свинцовых фантомах при облучении гамма-источниками / А. В. Жуков, И. Н. Соколова // Ядерная физика и инжиниринг. — 2023. — Т. 14, № 5. — С. 478-485.
Захаров, А. А. Моделирование радиационной защиты в программе GATE: методические указания / А. А. Захаров, Д. В. Кузнецов. — Санкт-Петербург : Издательство СПбГУ, 2022. — 48 с.
Иванов, В. П. Физика ядерных реакторов: учебное пособие / В. П. Иванов, К. А. Михайлов. — Москва : НИЯУ МИФИ, 2021. — 420 с. — ISBN 978-5-7262-2845-6.
Козлов, Д. А. Энергетическая зависимость дозовых распределений в воде при фотонном облучении / Д. А. Козлов, Е. В. Морозова // Известия вузов. Ядерная энергетика. — 2024. — № 1. — С. 112-121.
Кузнецов, И. В. Применение пакета GATE для моделирования в лучевой терапии : учебное пособие / И. В. Кузнецов, С. А. Петров. — Москва : Издательство РУДН, 2023. — 198 с. — ISBN 978-5-209-11234-5.
Лебедев, А. Н. Анализ погрешностей при моделировании дозовых распределений методом Монте-Карло / А. Н. Лебедев, $. $. $$$$$$$$ // $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$. $$$$$ $: $$$$$$, $$$$$$$$$$. — $$$$. — № $. — $. $$-$$.
$$. $$$$$$$$, $. $. $$$$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$ / $. $. $$$$$$$$, $. $. $$$$$$$ // $$$$$$$$$$$ $$$$$$$. — $$$$. — № $. — $. $$-$$.
$$. $$$$$$$, $. $. $$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$ $$$$$-$$$$$$$$$ / $. $. $$$$$$$, $. $. $$$$$$$$ // $$$$$$$ $$$$$$$. — $$$$. — $. $$$, № $. — $. $$$-$$$.
$$. $$$$$$$, $. $. $$$$$$$$$$$$$ $$$$$$$$$ $ $$$ $$$$ $ $$$$$$$$$$$$ $$$$$$$ $$$$$ : $$$$$$$$$$ / $. $. $$$$$$$, $. $. $$$$$$$$$. — $$$$$ : $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$$-$$$$-$.
$$. $$$$$$, $. $. $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$ $ $$$ $$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$$$ / $. $. $$$$$$, $. $. $$$$$$$ // $$$$$$$$$$$ $$$$$$$$ $ $$$$$$$. — $$$$. — № $. — $. $$-$$.
$$. $$$$$$, $. $. $$$$$$$$$$$ $$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$-$$$$$$$$$ / $. $. $$$$$$, $. $. $$$$$$$$ // $$$$$$$$ $$$. $$$$$ $$$$$$$$$$. — $$$$. — $. $$, № $. — $. $$$-$$$.
$$. $$$$$$$, $. $. $$$$$$$$$ $$$$$$ $$$$$-$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$ : $$$$$$$ $$$$$$$ / $. $. $$$$$$$, $. $. $$$$$$$. — $$$$$$$$$$$ : $$$$$$$$$$$$ $$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$$-$$$$-$.
$$. $$$$$$$, $. $. $$$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$ $ $$$$$$ / $. $. $$$$$$$, $. $. $$$$$$$$ // $$$$$$ $$$$$$$$$$$ $$$$$$. — $$$$. — $. $$, № $. — $. $$$$-$$$$.
$$. $$$$$$$$, $. $. $$$$$ $$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$$$$$$ $$$$$$$$$$ / $. $. $$$$$$$$, $. $. $$$$$ // $$$$$$$$ $$$$$. $$$$$$. — $$$$. — № $. — $. $$-$$.
$$. $$$$$$$$, $. $. $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$ : $$$$$$$$$$$$ $$$$$$$$ / $. $. $$$$$$$$, $. $. $$$$$$$. — $$$$$$$$$$$$ : $$$$$$$$$$$$ $$$$, $$$$. — $$ $.
$$. $$$$$$$$, $. $. $$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$ $$ $$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$ $$$$$$$$ / $. $. $$$$$$$$, $. $. $$$$$$$ // $$$$$$$$$$$ $$$$$$. — $$$$. — № $ ($$$). — $. $$-$$.
$$. $$$$$$$$, $. $. $$$$$$$$$$$ $$$$$$$ $$$$$$ $ $$$$$$$$$$$ $$$$$$ $$$$ / $. $. $$$$$$$$, $. $. $$$$$$ // $$$$$$-$$$$$$$$$$$ $$$$$$$$$ $$$$$. $$$$$$$$$$$. $$$$$$$$$$$$$$$$. $$$$$$$$$$. — $$$$. — № $. — $. $$-$$.
$$. $$$$$$, $. $. $$$$$$ $$$$$$$$$$$ $$$ $ $$$ $$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$ $$$$$$$$ / $. $. $$$$$$, $. $. $$$$$$ // $$$$$$$ $$$$$$. — $$$$. — $. $$, № $. — $. $$$-$$$.
$$. $$$$$$$, $. $. $$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ / $. $. $$$$$$$, $. $. $$$$$$ // $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$. — $$$$. — № $. — $. $$-$$.
$$. $$$$$$$$, $. $. $$$$$$$$$ $$$$$$ $ $$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$ $ $$$$$$$ $$$$$-$$$$$$$$$ / $. $. $$$$$$$$, $. $. $$$$$$$ // $$$$$$$ $$$$$$$. — $$$$. — $. $$$, № $. — $. $$-$$.
$$. $$$$$$$, $. $. $$$$$$$$$$ $$$$$$ $ $$$$$$$$$ $$$$: $$$$$ $ $$$$$$$$$ $$$ $$$$$ $$$$$$$$$$$ $$$$$$ / $. $. $$$$$$$, $. $. $$$$$ // $$$$$$$$$$$ $$$$$$$. — $$$$. — № $. — $. $$-$$.
2026-05-28 11:47:09
**Краткое описание работы** Данная курсовая работа посвящена моделированию дозового распределения гамма-излучения в программе GATE, основанной на методе Монте-Карло. **Актуальность** темы обусловлена необходимостью точного расчета поглощенных доз в лучевой терапии и радиационной безопасности, гд...
2026-05-28 12:01:53
Вот краткое описание работы, составленное в соответствии с вашей структурой и требованиями. --- ### Краткое описание работы **Актуальность** данной работы обусловлена необходимостью точного и безопасного планирования лучевой терапии и радиационной защиты. Экспериментальные методы измерения доз...
2026-05-28 11:42:38
**Краткое описание работы** Данная курсовая работа посвящена моделированию дозового распределения гамма-излучения в программе GATE (Geant4 Application for Tomographic Emission), основанной на методе Монте-Карло. **Актуальность** исследования обусловлена необходимостью точного прогнозирования пог...
Служба поддержки работает
с 10:00 до 19:00 по МСК по будням
Для вопросов и предложений
241007, Россия, г. Брянск, ул. Дуки, 68, пом.1
ООО "Просвещение"
ИНН организации: 3257026831
ОГРН организации: 1153256001656