4 курс / Лучевая диагностика / Радиобиологическое_и_дозиметрическое_планирование_лучевой_и_радионуклидной
.pdfгде Mp – сигнал, измеренный ионизационной камерой в точке Р и скорректированный на влияние температуры, давления воздуха и рекомбинацию ионов.
Ионизационная камера должна иметь соответствующий колпачок, обеспечивающий электронное равновесие, и калибровочные
коэффициенты по экспозиционной дозе в воздухе Nх |
или по воз- |
душной камере в воздухе Nk. |
|
Шаг 1: Определение Хр в точке Р через |
|
Х p M p N x . |
(2.17) |
Шаг 2: Определение (Квоз)воз в точке Р через
(Квоз ) воз= 0,876·Хр . (2.18)
Возможно также непосредственное определение (Квоз)воз, если из-
вестно Nk:
(Квоз ) воз=Mp·Nk . (2.19)
Шаг 3: Определение ионизационной кермы в бесконечно малой массе m любого другого вещества (например воды), находящейся
в воздухе: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
en |
|
m |
|
|
|
|
|
|
|||||
(K |
m )воз |
(Kвоз )воз |
|
|
|
|
|
, |
(2.20) |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
воз |
|
|
|
|
|
|
|
|
m |
|
|
en |
|||||||
где |
|
|
– отношение массовых коэффициентов поглощения |
|||||
|
|
|
|
|||||
|
|
|
|
|
|
|||
|
|
|
воз |
|||||
энергии в веществе и в воздухе, усредненных по спектру пучка. |
||||||||
Шаг 4: |
|
Определение ионизационной кермы в сферической мас- |
се вещества, центрированной в точке Р и имеющей радиус rвещ для обеспечения электронного равновесия в точке Р:
Квещ воз К m воз к(rвещ ) , |
(2.21) |
где к(rвещ ) – поправочный множитель, учитывающий ослабление
фотонов в сферической массе вещества и аппроксимируемый формулой:
|
|
|
|
|
en |
|
|
|
|
к(r |
вещ |
) exp |
|
|
|
|
r |
, |
(2.22) |
|
|||||||||
|
|
|
|
|
вещ |
|
|
||
|
|
|
|
|
вещ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
131 |
|
|
|
где en вещ – массовый коэффициент поглощения энергии для
данного вещества и ρ – плотность вещества. Для воды, обычно выбираемой в качестве среды, к 0.985 для Со-60 и k = 1.0 для фотонов низких энергий.
Шаг 5: «Доза в небольшой массе вещества, находящегося в свободном пространстве».
D’вещ определяется из выражения:
D |
|
(K |
|
) |
|
|
0,876 |
|
en вещ |
|
X |
|
к(r |
) |
|
|
|
вещ |
воз |
|
|
|
|
|
р |
, (2.23) |
|||||||
вещ |
|
|
|
|
|
вещ |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
воз |
|
|
|
|
|
|
где β – коэффициент, равный 1,003; 1,001 и 1,0 для Со-60, Cs – 137 и рентгеновского излучения ниже 350 кВ.
Произведение 0,876 en / вещвоз обычно называют конверси-
онным фактором «рентген к грею» и обозначается fвещ. Для «дозы в небольшой массе вещества в воздухе», предполагаем β ≈ 1, тогда можно записать:
вещ |
|
вещ |
|
|
вещ |
|
. |
(2.24) |
D |
f |
|
Х |
|
K(r |
) |
|
|
2. Фантомные материалы
Дозиметрическое планирование лучевой терапии в соответствии с рекомендациями международных организаций должно проводиться с очень высокой точностью. Допустимая погрешность при расчете дозы не превышает ~ 3 %. Для достижения такой точности практически современные системы дозиметрического планирования опираются на экспериментальные измерения распределения поглощенной дозы в различных фантомах. Не менее важным обстоятельством, способствующим широкому распространению фантомов, является необходимость определения дозовых распределений внутри тела пациента. Современное состояние дозиметрической техники не позволяет, как правило, выполнить такие измерения непосредственно в теле человека.
Сегодня существует большое количество разных видов фантомов. Они изготовляются из разных материалов, имеют разные размеры и форму. Фантомы бывают гомогенные, гетерогенные и ан-
132
тропоморфные, последние моделируют тело типового человека в целом. Но наибольшее распространение получили фантомы из воды и близких к воде материалов (по эффективному атомному номеру и плотности электронов). В табл. 2.1 представлены физические свойства материалов, наиболее часто используемых в фантомах, а также физические свойства различных видов ткани человека.
|
|
|
|
|
Таблица 2.1 |
|
Физические свойства фантомных материалов |
|
|||
|
|
|
|
|
|
|
Химический |
Плотность, |
Эл./см3 |
|
Zэф |
Материал |
состав |
(г/см3) |
(х1023) |
|
(фотоэффект.) |
Воздух |
|
0,001293 |
3,01 |
|
7,64 |
Вода |
H2O |
1 |
3,34 |
|
7,42 |
Мышцы |
|
1 |
3,36 |
|
7,42 |
Жир |
|
0,916 |
3,48 |
|
5,92 |
Кость |
|
1,85 |
3,00 |
|
13,8 |
Полистерен |
(C8 H8)n |
1,03-1,05 |
3,24 |
|
5,69 |
Плексиглас |
(C5 O2 H8)n |
1,16-1,20 |
3,24 |
|
6,48 |
(perspex, licite) |
|
|
|
|
|
Полиэтилен |
(CH2)n |
0,92 |
3,44 |
|
6,16 |
Парафин |
CnH2n+2 |
0,87-0,91 |
3,44 |
|
5,42 |
Смесь D |
Параф.:60,8% |
0,99 |
3,41 |
|
7,05 |
|
Полиэт.:30,4% |
|
|
|
|
|
MgO:6,4% |
|
|
|
|
|
TiO2:2,4% |
|
|
|
|
Смесь М3 |
Параф.:100ч. |
1,06 |
3,34 |
|
7,05 |
|
MgO:29,06ч. |
|
|
|
|
|
CaCO3:0,94ч. |
|
|
|
|
Твердая вода |
Эпоксидно- |
1,00 |
3,34 |
|
7,35 |
|
резиновая смесь |
|
|
|
|
3. Процентная глубинная доза и ее свойства
3.1. Основные обозначения и определения
При расчете дозовых распределений и дозиметрическом планировании ЛТ используется довольно большое количество специаль-
133
ных понятий, терминов и обозначений. Имея это в виду, для большей наглядности и удобства приведем в табл. 2.2 и 2.3 предварительную сводку основных обозначений (по умолчанию, если в тексте не оговорены другие обозначения) и краткие определения этих величин.
|
|
|
Таблица 2.2 |
Обозначения геометрических величин, используемых в данной главе |
|||
|
|
|
|
Обозначение |
|
|
Значение обозначения |
|
|
|
|
d |
|
|
Глубина в ткани (фантоме) от поверхности, измеряемая па- |
|
|
|
раллельно оси пучка |
|
|
|
|
dmax (или dm) |
|
|
Глубина максимальной дозы |
dref |
|
|
Глубина ссылочной (опорной или рефересной) точки |
x |
|
|
Внеосевое расстояние от оси вращения коллиматора |
SAD |
|
|
Расстояние от источника до изоцентра (оси вращения гантри) |
SSD |
|
|
Расстояние от источника до облучаемой поверхности |
|
|
|
|
A (или r) |
|
|
Размер поля в изоцентре |
As (или rs) |
|
|
Размер поля на поверхности |
Am (или rm) |
|
|
Размер поля на глубине dmax |
Ad (или rd) |
|
|
Размер поля на глубине расчетной точки |
Aref (или rref) |
|
|
Ссылочный (опорный или рефересный) размер поля при ка- |
|
|
|
либровке (обычно 10х10 см2) |
ESQ |
|
|
Размер эквивалентного квадратного поля в расчетной точке |
|
|
|
|
|
|
|
Таблица 2.3 |
Сводка обозначений и кратких определений основных дозиметрических |
|||
|
|
величин, используемых в данной главе |
|
|
|
|
|
Обозначение |
|
Краткое определение |
|
|
|
|
|
W |
|
Наличие клинового фильтра |
|
D0 |
|
Калибровочная величина D/MU (MU – мониторная единица) в |
|
|
|
|
ссылочной точке для ссылочного размера поля (обычно 1 |
|
|
|
Гр/100 MU) |
PDD(d,As,SSD) или |
|
Процентная доза на глубине d в ткани для размера поля на |
|
P%(d,As,SSD) |
|
поверхности As и расстояния между источником и облучае- |
|
|
|
|
мой поверхностью SSD (в России часто применяется-РИП) |
TAR(d,Ad) |
|
Отношение ткань-воздух на глубине d для размера поля Ad в |
|
|
|
|
расчетной точке |
TMR(d,Ad) |
|
Отношение ткань-максимум на глубине d для размера поля Ad |
|
|
|
|
в расчетной точке |
|
|
134 |
|
Продолжение табл. 2.3 |
|
|
Обозначение |
Краткое определение |
TPR(d,Ad) |
Отношение ткань-фантом на глубине d для размера поля Ad в |
|
расчетной точке |
PSF(A) |
Пиковый фактор рассеяния для размера поля A, определенный |
|
в ткани на глубине dmax |
NPSF(A) |
Нормированный пиковый фактор рассеяния, определенный в |
|
ткани на глубине dmax |
Sp(dref,ESQ) |
Фактор рассеяния в фантоме для эквивалентного квадрата |
|
(ESQ) при dref, определенный в ткани на глубине dref |
Sc(A) |
Фактор рассеяния в коллиматоре для установленного поля А, |
|
называемый также "выходной фактор в воздухе" |
FOF(dref,A,ESQ) |
Выходной фактор поля для установленного поля А и эквива- |
|
лентного квадрата (ESQ) при при dref, определенный в ткани на |
|
глубине dref |
OAR(x,d,Ad) |
Внеосевое отношение для расстояния от оси пучка x на глубине |
|
d в ткани для размера поля Ad в расчетной точке |
WF(dref,Ad) |
Фактор клина: отношение выходных факторов с и без клина |
D(air,A,SAD) |
Поглощенная доза в воде (на MU) на CAX в воздухе на рас- |
|
стоянии SAD для размера поля в изоцентре А |
D(dref,Ad,SSD) |
Поглощенная доза в воде (на MU) на CAX в ткани на глубине |
|
dref на расстоянии источник-поверхность SSD для размера поля |
|
Ad в расчетной точке |
D(dmax,Ad,SSD) |
Поглощенная доза в воде (на MU) на CAX в ткани на глубине |
|
dmax на расстоянии источник-поверхность SSD для размера поля |
|
Ad в расчетной точке (=Dmax) |
D(d,As, SSD,W) |
Поглощенная доза в воде (на MU) на CAX в ткани на глубине d |
|
на расстоянии источник-поверхность SSD для размера поля As |
|
на поверхности с (возможно) клиновым фильтром |
D(d,Ad,SSD,W) |
Поглощенная доза в воде (на MU) на CAX в ткани на глубине d |
|
на расстоянии источник-поверхность SSD для размера поля Ad в |
|
расчетной точке с (возможно) клиновым фильтром |
D(x,d,Ad,SSD,W) |
Поглощенная доза в воде (на MU) на CAX точке на расстоянии |
|
x от CAX в ткани на глубине d на расстоянии источник- |
|
поверхность SSD для размера поля Ad в расчетной точке с (воз- |
|
можно) клиновым фильтром |
|
3.2. Процентная глубинная доза |
Процентная глубинная доза PDD (или P%) является одной из основополагающих величин в ЛТ. Она измеряется с высокой точ-
135
ностью в водных фантомах, и через нее определяются многие другие дозиметрические характеристики, используемые при дозиметрических расчетах в ЛТ. Процентная глубинная доза (или короче процентная доза) определяется как отношение поглощенной дозы на геометрической оси пучка в произвольной точке на глубине d к максимальной дозе на оси ( рис.2.2):
PDD(d, As , SSD) 100 D(d, As , SSD) / D(dmax , As , SSD), (2.25)
|
0 |
для E 0, 4 МэВ; |
где dmax |
|
для E 0, 4 МэВ. |
|
dm |
Рис. 2.2. Определение процентной глубинной дозы
На рис.2.3 приводятся типичные зависимости процентной дозы от глубины в водном фантоме для разных энергий источника или центрально-осевые дозовые распределения. Более детальные данные по процентным дозам для разных размеров полей и энергий источника даются в приложении. Для краткости эти распределения иногда называют глубинками. Форма этих кривых зависит от энергии источника, размера поля, расстояния от источника до поверхности фантома РИП (или в английской транскрипции SSD).
Характерной особенностью глубинных дозовых распределений является появление у этих кривых максимума для средних и высоких энергий излучения источника. Положение этого максимума dm
136
с ростом энергии смещается в сторону больших глубин, а район между поверхностью и dm называют областью накопления (англ. build up). Его появление связано с отсутствием электронного равновесия на малых глубинах (см. рис. 2.1). Как уже отмечалось в разделе 1, основная часть поглощенной дозы создается вторичными электронами, образующимися при взаимодействии фотонов с веществом. На малых же глубинах количество этих электронов недостаточно. Отметим одновременно, что глубинная зависимость кермы не имеет максимума.
|
100 |
|
10 MV |
25 MV |
|
|
|
|
|
||
Р% |
80 |
|
|
4 MV |
|
|
|
|
|
Co |
|
|
60 |
|
|
|
|
|
|
|
|
|
|
|
40 |
|
|
|
|
|
20 |
3.0 мм Cu |
|
|
|
|
|
|
|
|
|
|
5 |
10 |
15 |
20 |
d, см |
Рис. 2.3. Типичные глубинные дозовые распределения для разных энергий излу- |
|||||
чения источника |
|
|
|
|
3.3. Зависимость процентной глубинной дозы от размера и формы поля
Для полей с небольшими поперечными сечениями доза в водном фантоме на оси пучка почти целиком определяется первичным фотонным излучением. Под первичным здесь и далее понимаются фотоны, которые не испытали взаимодействия на пути от источника к точке расчета.
Сувеличением размера поля растет вклад в дозу от рассеянных
всреде фотонов. Так как эти фотоны имеют преимущественное направление вперед, особенно для высокоэнергетических фотонов, то относительный вклад в дозу от рассеянных фотонов растет с увеличением размера поля более быстро на больших глубинах.
137
В то же время на любой конкретной глубине (но равной или большей, чем dm) доза от первичных фотонов не меняется с изменением размера поля. Исключение представляют только поля с очень малыми поперечными размерами, так как в этом случае имеет место отсутствие уже поперечного (или бокового) электронного равновесия. Поэтому суммарная доза (от первичных и рассеянных фотонов) растет с увеличением размера поля быстрее, чем максимальная доза на оси пучка. Все это приводит к эффекту увеличения процентной дозы с увеличением размера поля (рис. 2.4.).
C увеличением энергии фотонов зависимость PDD от размера поля становится более слабой. Причиной такого эффекта является увеличение анизотропии рассеяния фотонов с ростом их энергии.
При постоянной площади поля процентная доза уменьшается с увеличением асимметрии поля независимо от энергии фотонов. Это происходит вследствие того, что с увеличением асимметрии увеличивается среднее расстояние между точками рассеяния фотонов и осью пучка. Как следствие – увеличивается вероятность поглощения рассеянных фотонов.
|
80 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
18MV |
|
|
|
70 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6MV |
|
|
|
60 |
|
|
|
|
|
|
|
|
|
P% |
50 |
|
|
|
|
|
|
Co-60 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
40 |
|
|
|
|
|
|
|
|
|
|
30 |
|
|
|
|
|
|
2mm Cu |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
204 |
6 |
8 |
10 |
12 |
14 |
16 |
18 |
20 |
|
. |
|
|
|
|
W,cm |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 2.4. Зависимость процентной дозы от стороны квадратного поля |
на глубине |
|||||||||
10 см для тормозного излучения разного спектра и радионуклида 60Со. Слой по- |
||||||||||
ловинного ослабления 2 мм Cu соответствует напряжению на трубке 240 кВ |
||||||||||
|
|
|
|
|
138 |
|
|
|
|
|
3.4. Переход от прямоугольных полей к эквивалентным квадратным полям
Облучение пациентов в лучевой терапии часто проводится полями с прямоугольным поперечным сечением. Так как прямоугольник имеет две степени свободы, накопление банка экспериментальных данных для таких полей затруднительно.
Удачным обстоятельством здесь является то, что с хорошей точностью возможен переход от прямоугольных к квадратным полям, имеющим одну степень свободы. Такой переход был предложен в работе [1] и осуществляется по принципу «Сохранение отношения S/P», где S – площадь; P – периметр.
Для прямоугольных полей с шириной a и длиной b
S |
|
а b |
(2.26) |
|
P |
2(a b) |
|||
|
|
Так как для квадратных полей a=b, то S/P=а/4, отсюда сторона
эквивалентного квадрата равна ESQ |
2 a b |
. Для эквивалентно- |
|
a b |
|||
|
|
||
го круглого поля радиус равен: |
|
|
rэкв |
4 |
|
S |
|
|
. |
(2.27) |
|
|
|
|
||||
|
|
|
|
||||
P |
|
|
|||||
|
|
|
|
|
|
Более точные значения ESQ были получены в работе [2] и приводятся в приложении.
3.5. Закон обратных квадратов
Излучение, испускаемое точечным изотропным источником, ослабляется в вакууме по закону, называемому "закон обратных квадратов". В соответствии с законом, если две расчетные точки расположены на расстояниях L1 и L2, то отношение доз в этих точках будет равняться:
D(L ) |
|
L |
2 |
|
1 |
|
2 |
. |
(2.28) |
|
|
|||
D(L2 ) |
|
L1 |
|
|
139
В случае облучательных установок таких, как ЛУЭ или аппараты с радионуклидом 60Со, источники не являются, строго говоря, истинно точечными и изотропными. Испускаемое ими фотонное излучение испытывает рассеяние в объеме источника (для 60Со) и мишени, в сглаживающем фильтре, в первичном коллиматоре, в коллимационных пластинах и, наконец, в воздухе. Тем не менее в достаточно широком интервале расстояний от источника (за исключением очень близких и далеких) этот закон выполняется с хорошей точностью.
Для более точных вычислений вводится понятие и определяется положение эффективного точечного источника. Для этого проводится измерение дозы (точнее выходного фактора (см. далее)) на оси пучка в воздухе на нескольких расстояниях от физического источника D(L). Полученные результаты наносятся на график в виде
зависимости y D(L0 ) / D(L) от расстояния (для ускорителей
обычно берется L0 = 100 см) и строится линейная регрессия (рис. 2.5). Пересечение аппроксимационной прямой с осью абцисс (y= 0) дает положение виртуального точечного источника.
3.6. Зависимость процентной глубинной дозы от расстояния источник – поверхность
При облучении РИП (в формулах, в основном, будет использоваться SSD) может меняться, поэтому необходимо уметь пересчитывать PDD, измеренную при конкретном значении РИП, к значениям PDD при произвольном РИП (рис.2.6). Будем считать, что в первом приближении вклад в дозу рассеянного излучения не меняется при изменении РИП. Тогда, учитывая экспоненциальное и геометрическое (по закону обратных квадратов) ослабление излучения, можно получить:
Р(d, A ,L ) 100 |
(L d |
m |
)2 |
|
e |
-μ(d -d |
|
) |
K |
|
; |
|
(2.29) |
||||||||
1 |
|
|
|
|
|
|
m |
|
|
|
|||||||||||
|
(L1 d )2 |
|
|
|
|
|
|
|
|||||||||||||
s |
1 |
|
|
|
|
|
|
|
|
|
|
s |
|
|
|
||||||
|
|
|
(L d |
m |
)2 |
|
|
- |
|
(d -d |
|
) |
|
|
|
|
|
||||
Р(d, A ,L ) 100 |
|
2 |
|
|
|
|
e |
|
|
m |
K |
|
, |
(2.30) |
|||||||
|
(L2 d )2 |
|
|
|
|
|
|
||||||||||||||
s |
1 |
|
|
|
|
|
|
|
|
|
|
|
s |
|
|
||||||
|
|
140 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|