calculo de dosis 3.2 pencil beam

1

Objetivos: Comprender la forma como se calcula la dosis empelando el método de Pencil Beam.

Calculo de Dosis3.2 Pencil Beam

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Dr. Willy H. GerberInstituto de Fisica

Universidad AustralValdivia, Chile

2

Absorción

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

z

Radiación ionizante al penetrar materia:

Fotones lateralesgenerados por Rayleigh

Scattering con generación de electrones en la distancia z

3

Absorción en un volumen dV

En tres dimensiones debemos considerar que la Intensidad decrece en función del radio:

R r

Muestra

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-2-4-Interaccion-rayos-gamma-materia-08.08

4

Absorción dependiente de la energía

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Si la constante de absorción es dependiente de la energía del fotón:

El flujo “extraído” por volumen para un espectro incidente Φ(E) será:

Si el tejido presenta un coeficiente de absorción variable

el flujo “extraído” por volumen será:

5

Modos de Scattering

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Pero que significa “extraer” del haz

Rayleigh

Compton

Fotoeléctrico

Campo de Núcleo Electrón e-

Positrón e+

Pares-núcleo

Campo de un electrón Electrón e-

Positrón e+

Pares-electrón

Fotón

Fotón + Electrón Electrón

6

Situación geométrica

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

El flujo absorbido en un punto es en consecuencia en la aproximación Pencil beam igual a:

En donde la propagación depende del coeficiente de absorción total:

y la absorción en el volumen dV del coeficiente parcial:

7

Aproximaciones Pencil Beam

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

z

1. Los fotones desviados por scattering Rayleigh solo se consideran en la reducción de la intensidad del haz pero no en lo que se refiere a sus consecuencias.

Aproximación Pencil BeamAt

enua

ción

Atenuación en AguaNIST

Para agua con densidad de 1 g/cm3 el coeficiente de atenuación en energías del orden de 1-6 MeV es 10-5 1/cm.Para anchos de 101 cm la perdida de energía por scattering Rayleigh es del orden de 10-4 o sea despreciable.

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

9

Aproximaciones Pencil Beam

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

z

2. Se asume que los electrones generados en el scattering tienen un camino libre despreciable depositando la energía en el mismo punto de la interacción.

Aproximación Pencil Beam

Stopping Power en AguaNIST

Para un stopping power de 2 MeV cm/g y una densidad de 1 g/cm3 y energías entre 1 y 6 MeV el camino recorrido es de algunos cm lo que constituye un ERROR NO DESPRECIABLE.

Para estudiar el error que se comete con esta suposición debemosanalizar el nivel de stopping power a lo que están expuestos los electrones:

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Aproximación Pencil Beam

Si comparamos el stopping power de tejido, hueso, aire (cavidades) y aguaobservamos que todos son similares y que el real problema es la densidad del medio

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Aproximación Pencil Beam

Por ello los principales errores del método Pencil Beam se referirán a problemas de localización de la dosis:

Aire (cavidad)

Tejido

Liquido (agua)Hueso

Dosis segúnPencil Beam

Haz

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Dosis realCavidades formanun caso extremo

13

Aproximaciones Pencil Beam

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Φ(0)

Φ(z)

z

3. Consideramos que los fotones que resultan del scattering son también absorbidos como energía sin estudiarse su propagación posterior.

14

Aproximaciones Pencil Beam

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Otro de los errores que se observan es la falta de build-up (en realidad el método supone que se logra en forma instantánea al ocurrir el scattering) o sea la zona en que los electrones alcanzan la situación de equilibrio.

Para el caso monocromático y sin un efecto de dilución por la propagación radial según Pencil Beam se tendría una dosis igual a

Sin efecto build-up

Profundidad

% d

e la

dos

is m

áxim

a

100%

0%

15

Aproximaciones Pencil Beam

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

La medición de dosis en profundidad muestra que de considerarse la entrega de energía sin propagación posterior existiría ante todo un problema en la modelación de la superficie:

Profundidad [cm]

Fotones en agua

Dosis

[% v

alor

máx

imo]

Build up

Nota: la mayoría de los “Pencil beam” comerciales muestran un comportamiento build-up lo que indica que técnicamente incluyen una convolusion simple.

16

Situación geométrica

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

La dosis se calcula sumando sobre todas las contribuciones que pueden irradiar el punto en estudio:

17

Situación geométrica

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

En este caso debemos sumar sobre las contribuciones definidas en el flujo de la superficie Φ(ρ,E):

ρ

r

R

ρ

18

Situación geométrica

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

De esta forma se puede calcular la dosis en un punto r como:

En donde Π es el kernel de la integral y en este caso corresponde a:

Lo que permite la estimación de la dosis en función de la ubicación y propiedades del tejido y de las características del haz incidente con los errores principalmente debido a no considerar la difusión de los electrones desde la zona del scattering.

19

Ejercicio con planilla

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

El calculo numérico requiere de estudiar como los fotones se propagan por el tejido.

Se subdividen los 20 cm de tejido en intervalos de 1 cm

Fotones “perdidos”por Rayleigh

Fotones que“generan” electrones

20

Ejercicio con planilla

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Una simulación simple se deja realizar mediante una planilla Excel:

Constantes deabsorción

Simulación

Dosis

21

Ejercicio con planilla

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Una simulación simple se deja realizar mediante una planilla Excel:

Constantes deabsorción

Simulación

Dosis

22

Ejercicio con planilla

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Caso 1: Radiación única por tejido de densidad uniforme (fotón monocromático)

Dosis [Gy] Densidad [kg/m3]

Profundidad [cm]

23

Ejercicio con planilla

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Caso 2: Con zona de mayor densidad

Dosis [Gy] Densidad [kg/m3]

Profundidad [cm]

24

Ejercicio con planilla

www.gphysics.net – UFRO-2008-Master-Fisica-Medica-3-2-Pencil-Beam-08.08

Caso 3: A lo largo de un haz que se cruza con dos similares

Dosis [Gy] Densidad [kg/m3]

Profundidad [cm]