Modélisation Monte-Carlo d'un accélérateur linéaire pour la prise en compte des densités pulmonaires dans le calcul de la dose absorbée en radiothérapie stéréotaxique

Monte-Carlo model of a linear accelerator for the absorbed dose computation of Stereotactic Radiotherapy in presence of very low lung densities

• Auteur / Autrice : Sara Beilla
• Direction : Xavier Franceries, Luc Simon
• Type : Thèse de doctorat
• Discipline(s) : Radio-physique et imagerie médicale
• Date : Soutenance le 27/09/2016
• Etablissement(s) : Toulouse 3
• Ecole(s) doctorale(s) : École doctorale Génie électrique, électronique, télécommunications et santé : du système au nanosystème (Toulouse)
• Partenaire(s) de recherche : Laboratoire : Centre de Recherche en Cancérologie de Toulouse (2011-....)

Mots clés

Mots clés contrôlés

Monte-Carlo, Méthode de

Radiothérapie

Poumon -- Cancer

Stéréotaxie

Mots clés libres

Monte-Carlo

Radiothérapie

Hétérogénéités

Algorithmes

Stéréotaxie pulmonaire

Monte Carlo

Radiotherapy

Heterogeneities

Algorithms

Lung Stereotaxic

Résumé

Le calcul de la distribution de dose en Radiothérapie externe se fait en routine clinique à l'aide de Systèmes de Planification de Traitement (TPS) commerciaux. Les algorithmes de calcul de ces TPS ont énormément progressé ces dernières années. Cependant ils sont basés sur des approximations qui restent acceptables pour la plupart des conditions cliniques mais qui montrent leurs limites dans certains cas notamment avec des petites tailles de champ d'irradiation et/ou des faibles densités massiques dans un milieu. Or ces deux conditions sont pourtant réunies dans le cadre de la radiothérapie stéréotaxique des tumeurs bronchiques. Si quelques études ont été réalisées pour des densités massiques classiques de poumon, aucune n'a été réalisée pour des densités pulmonaires très faibles comme par exemple lorsque le patient est traité en inspiration profonde (" Deep Inspiration BreathHold ", i.e. DIBH). Mes travaux de recherche de thèse proposent une étude du calcul de dose pour différentes densités massiques et différentes tailles de champ en se basant sur un modèle Monte-Carlo (MC). La première étape modélise un accélérateur de type TrueBeam(r) (Varian, Palo Alto, CA) en utilisant les données du constructeur. Le modèle est construit à l'aide de la plateforme GATE basée sur la librairie Geant4. Les éléments principaux de la tête de l'appareil sont modélisés. Les espaces de phases (fichiers de particules) fournis par le constructeur au format " .IAEAphsp " sont situés en amont des mâchoires. Pour valider ce modèle, une série de champs simples (3x3 à 20x20 cm2) dans un fantôme d'eau sont implémentés pour des faisceaux de photons de 6X FF (" Flattening Filter "), 6X FFF, 10X FF et 10X FFF (" Flattening Filter Free "). Les résultats (profils, rendements de dose) sont comparés à des mesures de référence obtenues dans une cuve d'eau : respectivement 99% et 97% des points de dose des rendements et des profils respectent les critères de gamma index de 2%-2mm. Une fois le modèle validé, nous avons réalisé une série de simulations pour des champs de petites tailles (3x3 à 8x8 cm2) avec des fantômes hétérogènes de formes simples, pour lesquels la mesure reste accessible. Pour cette dernière, ont été insérés des films radio-chromiques dans des fantômes composés de plaques de PMMA et de deux types de liège de densité 0,12 et 0,24 correspondant respectivement aux poumons en DIBH et en respiration libre. Les résultats du modèle MC pour les quatre énergies ont été confrontés aux mesures expérimentales et aux algorithmes AAA et Acuros (Varian). De façon générale l'algorithme AAA surestime la dose au sein de l'hétérogénéité pulmonaire pour les petites tailles de champ et les faibles densités massiques. Par exemple, pour un champ de 3x3 cm2 et une densité de 0,12 au sein de l'hétérogénéité, une surestimation de la dose absorbée dans le poumon de 16% est mise en évidence pour l'algorithme AAA. Enfin, le modèle est utilisé pour trois cas non mesurables : un objet-test numérique cylindrique hétérogène, des données tomodensitométriques d'un patient en DIBH pour un champ fixe et en arc-thérapie en condition de stéréotaxie pulmonaire. Les résultats ont démontré respectivement pour les études sur TDM une surestimation de la dose dans la tumeur de 7% et 5,4% et dans le poumon de 14% et 9,6% par AAA. D'un point de vue clinique, cela se traduit par un sous-dosage du patient et donc un risque de récidive.

For clinical routine in external Radiotherapy, dose computation is achieved using commercial Treatment Planning Systems (TPS). Since ten years, TPS algorithms have been improved. However they include approximations that are acceptable in most of the clinical cases but they show their limits in some particular conditions for example in presence of small fields and/or low mass y media. And these two conditions are found in the context of stereotactic body radiation therapy (SBRT) of lung tumor. Some studies were published for standard lung densities but none for very low y like in lung during Deep Inspiration Breath Hold (DIBH). This work is a study of dose computation based on a Monte Carlo (MC) model, for different field sizes and mass densities. The first step was to model a TrueBeam(r) linac (Varian, Palo Alto, CA) using data furnished by the manufacturer. This model is built using the Geant4-based GATE platform. The main compounds of the linac head are modeled. Space phase files (i.e. particles files) are furnished by Varian in "IAEAphsp" format and are integrated to the model above the main jaws. To validate this model, a set of simple fields (from 3x3 to 20x20 cm2) in a water phantom is implemented for different photon energies: 6FF, 6FFF, 10FF and 10FFF (FFF = "Flattening Filter Free"). Percentage depth dose (PDD) and lateral profiles are compared to reference measurement in a water tank: respectively 99% and 97% of all the points of these curves passed the Gamma Index test (2% 2mm). Once this validation was completed, a set of simulation was achieved with small field sizes (3x3 to 8x8 cm2) for simple heterogeneous phantoms for which the measurement was achievable. For this purpose, radiochromic films were inserted in phantoms made of PMMA slabs and two types of cork. Cork densities were 0.12 and 0.24 that correspond respectively to lungs during DIBH and free breathing. Results of the MC model for four energies are compared to experimental measurements and to AAA and Acuros Varian's algorithms. AAA algorithm overestimates the dose inside the lung heterogeneity for small field sizes and low density. As an example in the case a 3x3 cm2 field, inside the heterogeneity of density 0.12 an over estimation of 16% in the lung is observed for AAA. The model is finally used for three non-measurable cases: a cylindrical digital reference object and computerized tomography data of a patient during DIBH with a static and stereotactic arc field. Results showed respectively for CT studies an overestimation of dose in the tumor of 7% and 5.4% and in the lungs of 14% and 9.6% by AAA. From a clinical point of view, this means under-dosing the patient and thus a risk of recurrence.