CC BY
150
ДОЗИМЕТРИЧЕСКОЕ ОБЕСПЕЧЕНИЕ РАДИОНУКЛИДНОЙ ДИАГНОСТИКИ И ТЕРАПИИ КОСТНЫХ МЕТАСТАЗОВ
Доля О.П., Матусевич Е.С., Клепов А.Н., Кураченко Ю.А.
Обнинский Государственный Технический Университет Атомной Энергетики Matusevich@ iate.obninsk.ru klyopov@obninsk.com
С целью создания системы дозиметрического обеспечения радионуклидной терапии (РНТ) диагностики и терапии пациентов с метастазами в кости разработаны математические модели, позволяющие определять активность радиофармпрепаратов (РФП) в метастазах и критических органах. Методы идентификации базируются на данных планарных сцинтиграфических и рентгенографических исследований пациентов и радиометрию экскреции РФП из организма. Расчеты распространения гамма- излучения в теле и системе регистрации проведены методом Монте-Карло, кинетика РФП - с применением камерных моделей, расчёты активности и поглощенных доз в очагах поражения части пациентов - с использованием воксельных моделей костных структур.
При множественных поражениях костей метастазами наиболее целесообразно использовать метод «внутреннего облучения», основанный на способности некоторых инкорпорированных радиофармпрепаратов избирательно накапливаться в метастатических очагах, и за счет облучения ß-частицами области нервных окончаний на границе здоровой и пораженной костных тканей (область периоста) снижать или полностью снимать болевой синдром. Этот метод и составляет существо паллиативной РНТ.
Зарубежный и отечественный клинический опыт показывает, что в большинстве случаев в той или иной мере паллиативный эффект достигается при уровнях активностей РФП, определяемых принятыми соглашениями (стандартными протоколами), при этом, как правило, поглощенные дозы оказываются сравнительно небольшими. Обычно при выборе уровней поглощенных доз в метастатических очагах учитывается только требования допустимого облучения критических внутренних органов и тканей пациента и безопасности окружающих его лиц. По существу, стандартные протоколы предписывают именно паллиативный режим РНТ с назначением достаточно невысоких уровней активностей РФП.
Относительно недавно обнаружена повышенная толерантность критических органов и тканей к облучению, что позволяет достигать высокого уровня очаговых доз и лучевого поражения метастазов. Например в [P. Anderson, 2001] сообщается о успешном использовании очень больших лечебных активностей (до 30 мКи/кг 153Sm) и, соответственно, больших очаговых доз при проведении РНТ пациентам с метастазами.
В случае применения высоких значений лечебных активностей резко возрастает роль дозиметрического обеспечения, которое включает:
- расчет поглощенных доз в очагах;
- проверку условий безопасного облучения критических структур пациента: красного костного мозга, остеогенного субстрата, крови, мочевого пузыря, почек, печени и др.;
- оценку и контроль уровня безопасного облучения контактирующих с пациентом лиц после выписки.
При разработке дозиметрического обеспечения РНТ диагностики и терапии костных метастазов на основе остеотропных РФП возникла необходимость обоснования и решения комплекса радиометрических и дозиметрических задач. С этой целью в клинике МРНЦ РАМН была проведена серия экспериментальных сцинтиграфических (В.В. Крылов, Г.А. Давыдов, Н.А. Олейник, О.Н. Спиченкова, авторы доклада) и рентгенологических исследований (В.В. Крылов, З.Н. Шавладзе,
А.В. Дементьев, авторы доклада) и радиометрических исследований (В.В. Крылов)
153
группы пациентов с костными метастазами, проходивших РНТ с Sm-оксабифором.
Авторами доклада был разработан ряд физико-математических моделей, обеспечивающих:
1. идентификацию величин накопленных активностей РФП в метастатических очагах ряда костей в процессе планарного сцинтиграфического обследования пациентов, что потребовало использования метода Монте-Карло для расчёта функции отклика системы регистрации и воксельных моделей костных структур;
2. идентификацию параметров кинетики накопления РФП в костных тканях и ряде критических структур на базе камерных моделей, что позволяет: - определить поглощенные дозы на костные ткани в целом и критические структуры (кровь, почки, мочевой пузырь);
3. обоснование амбулаторного варианта РНТ с 153Sm путем оценки уровня облучения окружающих лиц от остаточной активности в теле пациента и экскрециях;
Кроме того предложена формулировка модели идентификации параметров камерной модели с учетом ремоделирования здоровых костей, позволяющая определить эффективную площадь здорового эндоста и площадь поверхности костной ткани, пораженной MTS.
Разработано программное обеспечение на ФОРТРАНЕ, а также программные модули реализации расчетных воксельных моделей костных структур транспорта у-излучения применительно к программному комплексу MCNP.
Далее для иллюстрации возможностей разработанного дозиметрического обеспечения приведены некоторые результаты. О технологии их получения см. [13].
Результаты моделирования Монте-Карло функции отклика коллиматора на рассеянное у-излучения от источника в кости
Рис. 1. MCNP-модель расчета транспорта у- Рис. 2.
излучения в системе «коллиматор - тело
пациента».
Результаты дозиметрического аннализа РНТ (15^т-оксабифор) пациентов с костными метастазами по данным сцинтиграфического исследования
Разброс по группе пациентов вариаций поглощенных доз взятых по всем метастазам в организме
о ■I
?
о 5 о
140 120
100
1 1 5 7 9 11 И 15 17 19 21 23 25 27 29 Л 35 37 39 41 43 45 47
пациенты
□ нижняя граница вариации ■ верхняя граница вариации
поглощенных доз в метастазах по всему организму пациента --уровень минимальной терапевтической поглощенной дозы
Рис 3.
Расчет активностей и поглощенных доз в метастазах и здоровых костях пациентов с использованием воксельных моделей костных структур и
планарных сцинтиграмм
Пациент Ч. (рак предстательной железы, ведено 65 мКи).
1
( *.
X ■
А
I
Грудина
А
Ал
11,45 см
<->
у
Рис. 4. Срез томограммы позвонка Т^12.
Рис.5. Воксельная модель Т^12 и грудины.
Доза в очаге поражения (зоны 1 и 2 на рис.5) - 21 Гр. Дозы на критические структуры здоровой кости: на эндост - 0,6 Гр, на красный костный мозг - 0,14 Гр.
см
Решение задачи идентификации параметров кинетики РФП в костных тканях и критических структурах организма пациента
Ло
р
кровь /
кнр
г /
Н
почки
/
к
г /
и
Мочевой
пузырь /
Р)
л
В)
л лм (о
л лн (р
л ли (о
= -(к„Р + кмр + кНр + • Р(0 + крВ • В(0 + крм • М
= квр • Р(0 - (крв +1Г) • B(t)
= кмр • Р(0 - (крм + \) • М (t) = кнр • Р^)- (к + Лг) • Н(О
= к • н (О -л •и (О
л
Р(0) =1,
В(0) = М(0) = Н(0) = и(0) = 0 крм = 0 - ремоделирование костной ткани в пораженной кости отсутствует.
Пятикамерная модель динамики обращения РФП в организме пациента: Б - камера плазмы крови, B - камера здоровой кости, M - камера метастазированной кости, Н - камера почек, и - резервуар выведения активности (камера экскреции). кВР, кРВ, кмр, крм, кНР, к - параметры межкамерных коммуникаций (кинетические параметры).
Задача идентификации решается путем численной минимизация функционала невязки:
ф(()=Е ( - и(э ))2 + [ в(( )-ь;1ос ] • м()-шэ1ос ]},
i=1
где иЭ - экспериментальные значения функции экскреции, отнесенные к
13 3
введенной активности; Ь1ос, 1ос - экспериментальные значения локальных
функции накопления активностей (описываемых долевыми параметрами g1, g2) в соответствующих зонах интереса сцинтиграфического изображения - в здоровой и пораженной костных структурах в заданные моменты времени 13;
X = [кВР, квм, кн, крв, к, g1 g 2 2 - вектор идентифицируемых кинетических
параметров .
Результаты идентификации:
Площадь оставшейся здоровой трабекулярной поверхности скелета пациента Sb = 3,6 м2; более подробное обсуждение результатов в [2].
1) О.П. Доля, Е.С. Матусевич, А.Н. Клепов. «Математическое моделирование кинетики остеотропного радиофармпрепарата в организме пациентов с метастазами в кости» // Медицинская физика, №2, 2007г, С. 40-50.
2) О.П. Доля, А.Н. Клепов. «Дозиметрическая оценка паллиативной радионуклидной терапии и идентификация объемов пораженной метастазами костной ткани» // Медицинская физика, №1, 2008г.
3) О.П. Доля, А.Н. Клёпов, Ю.А. Кураченко, Е.С. Матусевич. «Моделирование Монте-Карло функции отклика гамма-камеры на гамма-излучение остеотропного радиофармпрепарата» // Медицинская физика, №2, 2008г.
