Добавил:
kiopkiopkiop18@yandex.ru Вовсе не секретарь, но почту проверяю Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

4 курс / Лучевая диагностика / Радиобиологическое_и_дозиметрическое_планирование_лучевой_и_радионуклидной

.pdf
Скачиваний:
0
Добавлен:
24.03.2024
Размер:
13.93 Mб
Скачать

где 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