DISPOSITIF ET METHODE DE REGULATION DE L ' INTENSITE D'UN FAISCEAU EXTRAIT D'UN ACCELERATEUR DE PARTICULES
Objet de l'invention
[0001] La présente invention se situe dans le domaine technique de la régulation de l'intensité d'un faisceau extrait d'un accélérateur de particules. [0002] La présente "• invention se rapporte à un dispositif destiné à la régulation rapide et précise de l'intensité d'un faisceau extrait d'un accélérateur de particules, et plus spécifiquement d'un cyclotron.
[0003] La présente invention se rapporte également à une méthode pour la régulation de l'intensité du faisceau extrait d'un accélérateur de particules.
[0004] La présente invention se rapporte enfin à l'utilisation de ce dispositif ou de cette méthode en protonthérapie et en particulier dans la technique de "Pencil Beam Scanning" .
Arrière-plan technologique et état de la technique [0005] Les cyclotrons sont des accélérateurs circulaires de particules qui sont utilisés pour accélérer des ions positifs ou négatifs jusqu'à des énergies de quelques MeV voire plus. Ce type d'appareils trouve des applications dans des domaines différents tels que l'industrie ou la médecine, plus précisément en radiothérapie pour la production de radio-isotopes ou en protonthérapie, en vue de traiter des tumeurs cancéreuses.
[0006] Les cyclotrons comprennent généralement cinq composants majeurs : la source d'ions qui génère les particules ionisées, le dispositif de confinement sous vide des particules ionisées, l' électroaimant qui produit le champ magnétique assurant le guidage des particules ionisées, le système accélérateur haute fréquence destiné à accélérer les particules ionisées, et le dispositif d'extraction permettant de dévier les particules ionisées de leur trajectoire d'accélération puis de les évacuer hors du cyclotron sous forme d'un faisceau à haute énergie cinétique. Ce faisceau est ensuite dirigé vers le volume cible.
[0007] Dans la source d'ions d'un cyclotron, les ions sont obtenus par ionisation, dans une enceinte fermée, d'un milieu gazeux constitué d'un ou plusieurs gaz, au moyen d'électrons fortement accélérés par résonance électronique cyclotronique sous l'action d'un champ magnétique haute fréquence injecté dans l'enceinte.
[0008] De tels cyclotrons peuvent être utilisés en protonthérapie. La protonthérapie vise à délivrer une dose élevée dans un volume cible à traiter bien défini tout en épargnant les tissus sains entourant le volume considéré. En comparaison à la radiothérapie conventionnelle (rayons X), les protons présentent l'avantage de déposer leur dose à une profondeur précise dépendant de l'énergie (pic de Bragg) . Plusieurs techniques pour distribuer la dose dans le volume cible sont connues .
[0009] La technique mise au point par Pedroni et décrite dans « The 200-Mev proton therapy project at the Paul Scherrer Insti tute : conceptual design and practical realization » MEDICAL PHYSICS, JAN. 1995, USA, vol .22, no . l , pages 37-53 , XP000505145 ISSN : 0094 -2405, consiste à découper le volume cible en volumes élémentaires appelés « voxels ». On dirige le faisceau vers un premier voxel, et
lorsque la dose prescrite est atteinte, on interrompt l'irradiation en déviant brusquement le faisceau au moyen d'un aimant rapide (fast kicking magnet) . On règle alors un aimant de balayage de manière à diriger le faisceau vers un voxel suivant, et on réintroduit le faisceau de manière . à irradier ce voxel suivant. Ce processus est répété jusqu'à irradiation du volume cible en entier. Un des inconvénients de ce procédé est que, en raison des interruptions et rétablissements successifs du faisceau entre deux voxels, le temps de traitement est important, et peut atteindre plusieurs minutes dans des conditions typiques.
[0010] La demande de brevet O00/40064 de la
Demanderesse décrit une technique améliorée, dite « pencil beam scanning », dans laquelle le faisceau ne doit pas être interrompu entre l'irradiation de chaque voxel individuel.
Le procédé décrit dans ce document consiste à déplacer le faisceau de manière continue de manière à "peindre" le volume cible couche après couche.
[0011] En opérant simultanément un déplacement du faisceau et une variation de l'intensité de ce faisceau, on parvient à conformer exactement la dose à délivrer au volume cible. La régulation de l'intensité du faisceau de protons est réalisée indirectement par une action sur le courant d'alimentation de la source d'ions. On utilise dans ce but un régulateur qui permet de réguler l'intensité du faisceau de protons. Toutefois, cette régulation n'est pas optimale.
[0012] Une autre technique utilisée en protonthérapie est la technique dite de « Double Diffusion ». Dans cette technique, la modulation de la profondeur d'irradiation (c'est-à-dire de l'énergie), est réalisée à l'aide d'une roue dite roue de modulation tournant à une vitesse de l'ordre de 600 tr/min. Les parties absorbantes de ce modulateur sont constituées d'un
matériau absorbant, tel que le graphite ou de lexan. Lors de la fabrication de ces roues de modulation, la modulation en profondeur obtenue est assez proche des prédictions . L'uniformité reste malgré tout en dehors des spécifications désirées. Pour atteindre les spécifications sur l'uniformité, plutôt que de ré-usiner les roues de modulation, il est moins onéreux d'utiliser une régulation de l'intensité du faisceau qui soit synchronisée sur la vitesse de rotation du modulateur d'énergie. La fonction de modulation est donc établie pour chaque modulateur d'énergie et est utilisée comme trajectoire fournie comme consigne au régulateur de l'intensité du faisceau. Une régulation rapide et précise de l'intensité du faisceau extrait d'un accélérateur de particules est donc nécessaire également dans les techniques de double diffusion utilisant une telle roue de modulation.
Buts de ' invention
[0013] La présente invention vise à fournir un dispositif et une méthode destinés à la régulation de l'intensité d'un faisceau extrait d'un accélérateur de particules, qui ne présente pas les inconvénients des procédés et dispositifs de l'état de la technique.
Résumé de 1 ' invention
[0014] La présente invention se rapporte à un dispositif de régulation de l'intensité du faisceau extrait d'un accélérateur de particules, tel qu'un cyclotron, utilisé par exemple pour la protonthérapie, lesdites particules étant générées à partir d'une source d'ions, caractérisé en ce qu' il comporte au moins :
- un comparateur déterminant un écart entre un signal digital représentatif de l'intensité du faisceau mesurée
à la sortie de l'accélérateur et une valeur de consigne de l'intensité du faisceau;
- un prédicteur de Smith, qui détermine, à partir de cet écart, une valeur corrigée de l'intensité de faisceau; - une table de correspondance inversée, fournissant, .à partir de la valeur corrigée de l'intensité de faisceau, une valeur de consigne pour l'alimentation du courant d'arc de la source d'ions. [0015] En outre, le dispositif selon l'invention peut comprendre un convertisseur analogique-digital convertissant le signal analogique directement - représentatif de l'intensité du faisceau mesuré à la sortie de l'accélérateur, et fournissant un signal digital. [0016] De préférence, le dispositif selon l'invention comprendra en outre :
- un filtre passe-bas filtrant le signal analogique directement représentatif de l'intensité du faisceau mesuré à la sortie de l'accélérateur, et fournissant un signal analogique filtré; - un régulateur à avance de phase échantillonnant ledit signal analogique filtré, compensant le retard de phase introduit par le filtre passe-bas, et fournissant un signal digital au comparateur.
[0017] Avantageusement, le dispositif de l'invention comporte des moyens de mise à jour du contenu de la table de correspondance inversée .
[0018] La fréquence d'échantillonnage est de préférence comprise entre 100 kHz et 200- kHz, et la fréquence de coupure du filtre passe-bas est de préférence comprise entre 2 et 6 kHz.
[0019] La présente invention concerne également une méthode de régulation, au moyen d'un dispositif de régulation digital fonctionnant à une fréquence d'échantillonnage donnée, de l'intensité du faisceau
extrait d'un accélérateur de particules, tel qu'un cyclotron, utilisé par exemple en protonthérapie, lesdites particules étant générées à partir d'une source d'ions, caractérisée en ce qu'elle comprend au moins les étapes suivantes :
- on mesure l'intensité du faisceau à la sortie de l'accélérateur de particules;
- on compare un signal digital représentatif de la mesure de l'intensité du faisceau avec la valeur de consigne de l'intensité du faisceau;
- on détermine, au moyen d'un prédicteur de Smith, une valeur corrigée de l'intensité de faisceau;
- on détermine, à partir de cette valeur corrigée de l'intensité de faisceau, au moyen d'une table de correspondance inversée, une valeur de consigne pour l'alimentation du courant d'arc de la source d'ions.
[0020] De préférence, dans la méthode selon l'invention, après la mesure de l'intensité du faisceau à la sortie de l'accélérateur de particules, on convertit le signal analogique directement représentatif de l'intensité du faisceau mesurée au moyen d'un convertisseur analogique digital pour obtenir un signal digital .
[0021] Selon une forme d'exécution de la méthode selon l'invention, - on filtre le signal analogique directement représentatif de l'intensité du faisceau mesurée au moyen d'un filtre passe-bas, donnant un signal analogique filtré;
- on échantillonne ledit signal filtré, et on compense le retard de phase introduit par le filtrage pour obtenir un signal digital.
[0022] Avantageusement, la correspondance entre une valeur pour l'alimentation du courant d'arc de la source d'ions et une valeur de l'intensité du faisceau mesurée à
la sortie de l'accélérateur est déterminée préalablement à la régulation.
[0023] Avantageusement, dans la correspondance entre une valeur de l'intensité du faisceau mesurée à la sortie de l'accélérateur et une valeur pour l'alimentation du courant d'arc de la source d'ions, les valeurs de l'alimentation du courant d'arc correspondant aux valeurs d'intensité de faisceau supérieures à une limite sont remplacées par la valeur d'alimentation du courant d'arc correspondant cette limite.
[0024] Enfin, la présente invention se rapporte également à l'utilisation du dispositif et de la méthode de l'invention en protonthérapie et en particulier dans les technique de "Pencil Beam Scanning" et de « double diffusion ».
Brève description des figures
[0025] La figure 1 représente un dispositif de régulation de l'intensité d'un faisceau extrait d'un accélérateur de particules suivant l'art antérieur.
[0026] La figure 2 représente la caractéristique du système, c'est-à-dire la correspondance entre une valeur 1^ pour l'alimentation du courant d'arc de la source d'ions et une valeur I„ de l'intensité du faisceau mesurée à la sortie de l'accélérateur.
[0027] La figure 3 représente un mode de réalisation d'un dispositif de régulation de l'intensité d'un faisceau extrait d'un accélérateur de particules selon l'invention. [0028] La figure 4 représente un second mode de réalisation d'un dispositif de régulation de l'intensité d'un faisceau extrait d'un accélérateur de particules suivant 1 ' invention.
Problèmes à la base de la présente invention [0029] En utilisant une régulation classique, par exemple PID, pour la mise en œuvre de la technique dite
"pencil beam scanning" telle que décrite dans la publication O00/40064 de la Demanderesse, on est confronté aux problèmes décrits ci-après.
[0030] Ainsi que le montre la figure 1, une valeur de consigne Ic de l'intensité du faisceau est fournie à un régulateur PID classique 10, qui détermine une valeur Ij^ du courant d'arc de la source d'ions 20. L'intensité du faisceau est mesurée au moyen d'une chambre d'ionisation 30, et le signal correspondant I est comparé à l'aide d'un comparateur "90 à la valeur de consigne le pour fournir un signal d'erreur ε. Selon la technique de balayage continu du faisceau, il est indispensable que l'intensité du faisceau varie de manière simultanée avec le déplacement, de manière à obtenir la conformité de la dose délivrée .
[0031] Un tel système présente les difficultés suivantes :
- un temps mort pur important est dû au temps de parcours important d'une particule entre son émission par la source d'ions 20 et sa sortie de la machine;
- la caractéristique du système, liant l'intensité du faisceau extrait de l'accélérateur de particules I à l'intensité du courant d'arc de la source d'ions I est fortement non linéaire, ainsi que le montre la figure 2;
- de plus, cette caractéristique peut varier au cours du temps, ainsi que le montrent les courbes en trait interrompu de la figure 2. Cette variation peut survenir rapidement en raison du chauffage ou du refroidissement du filament de la source d'ions lors de sa mise en
service. Elle peut également provenir du vieillissement du filament. Ces deux phénomènes conduisent à des variations de la caractéristique avec des constantes de temps très différentes; - le système est fortement bruité. L'intensité du faisceau généré par la source d'ions présente un bruit important, en particulier à la fréquence d'échantillonnage utilisée pour la mesure.
[0032] La régulation d'un tel système en utilisant les méthodes classiques de régulation telles que les techniques de feedforward, de feedback par action proportionnelle, intégrale et dérivée (PID) et de boucles en cascade a été évaluée. En raison du temps mort pur important, toutes ces méthodes donnent des réponses soit trop lentes, soit instables. Les méthodes classiques ne permettent pas non plus d'adresser le problème d'une caractéristique du système fluctuant en fonction du temps en utilisant une valeur moyenne des caractéristiques sur une période donnée, car les variations de gain d'une réponse à l'autre sont dans un rapport très important.
[0033] L'évolution de la caractéristique dépend de deux phénomènes bien découplés : le premier, de constante de temps brève, correspond au conditionnement de la source d'ions, c'est-à-dire à sa température. Le fonctionnement normal, continu ou intermittent à rapport cyclique élevé, chauffe la source d'ions rapidement. Ce temps d'établissement rapide en température pourrait permettre de travailler en boucle ouverte, c'est-à-dire sans tenir compte de la caractéristique réelle du système, en utilisant les méthodes classiques, pendant le temps de conditionnement. Cependant, ce compromis limite fortement l'utilisation d'une méthode classique en fonctionnement intermittent à rapport cyclique moyen, qui correspond souvent au mode de fonctionnement utilisé.
[0034] Le second phénomène, de constante de temps plus longue, est dû au vieillissement du filament et de la source d'ions elle-même. Cette évolution plus lente de la caractéristique pourrait donc donner lieu à l'utilisation d'une caractéristique moyenne du système. L'utilisation d'une caractéristique moyenne conduit cependant à une régulation soit trop lente, soit instable.
[0035] Il apparaît donc comme évident que les méthodes classiques de régulation ne peuvent pas résoudre de manière satisfaisante les problèmes du contrôle d'un tel système, c'est-à-dire un temps mort pur largement supérieur à la constante de temps principale du système (environ 4 fois) et une caractéristique non-linéaire évolutive et nécessitant une méthode de régulation adaptative. [0036] La régulation rapide et précise de l'intensité du faisceau extrait d'un accélérateur de particules se heurte donc à de nombreuses difficultés . Une telle régulation rapide et précise est cependant importante pour l'application de la technique du « pencil beam scanning ».
Description d'une forme d'exécution préférée de l'invention [0037] La présente invention se propose par conséquent de résoudre plus spécifiquement ce problême en utilisant selon une forme d'exécution préférée un dispositif de régulation 10 représenté à la figure 3 avec l'alimentation du courant d'arc de la source d'ions 20. La source d'ions produit un faisceau d'ions, qui est accéléré au cours de son parcours dans l'accélérateur, en est extrait, et traverse un dispositif de mesure 30 de l'intensité du faisceau à la sortie de l'accélérateur. Ce dispositif de mesure 30 peut être par exemple une chambre d' ionisation.
[0038] Le régulateur suivant l'invention à été appliqué pour un cyclotron présentant les caractéristiques exemplatives et non limitatives suivantes :
- énergie fixe: 235 MeV - temps mort pur : 60 μsec. Ce temps mort pur correspond au temps de parcours des ions dans l'accélérateur. Il correspond donc directement au temps nécessaire pour mesurer l'influence d'une modification de la consigne du courant d'arc de la source d'ions sur l'intensité du faisceau d'ions extraits de la machine
- constante de temps principale : 15 μs . Elle donne une indication sur le temps nécessaire à l'établissement, en boucle ouverte, de la réponse du système à une modification de consigne. - caractéristique du système fortement non-linéaire, ce qui conduit à une caractéristique en boucle ouverte correspondant quasiment à celle d'un système à dynamique hybride (tout ou rien) .
- évolution de la caractéristique au cours du temps. - signal mesuré très bruité. En effet, la source d'ions est instable, ce qui conduit un niveau de bruit très important pour l'intensité du faisceau après extraction.
Le rapport bruit/signal observé est de l'ordre de 150%.
En conséquence, lors d'une mise en œuvre digitale du régulateur, les fréquences d'échantillonnages retenues engendrent un rapport signal / bruit très faible.
[0039] Dans le dispositif de régulation de l'invention, représenté à la figure 3, les étapes suivantes sont réalisées : — la valeur de consigne de l'intensité du faisceau Ic est fournie sous la forme d'un signal analogique 0-10 V (10 V correspondant à une intensité du faisceau de 300 nA) ;
- l'intensité de faisceau est mesurée au moyen d'une chambre d'ionisation 30 et la mesure I„ M est fournie au dispositif de régulation 10 au moyen d'un signal analogique 0-15 μA (15 μA correspondant à une intensité du faisceau de 300 nA) ;
- ce signal analogique I est converti par un convertisseur
50 en un signal digital I ;
- ce signal I est comparé par le comparateur à la consigne
Ic pour fournir un signal d'erreur ε; - ce signal d'erreur ε est fourni au régulateur de type « Prédicteur de Smith » 80;
- la sortie IP du prédicteur de Smith 80 est alors fournie à l'entrée d'une table de correspondance inversée 40. La table de correspondance 40 fournit de manière numérique la relation non-linéaire entre le courant d'arc de la source d'ions I A, et l'intensité du faisceau I.M, d'ions extrait de l'accélérateur. Elle permet donc d'identifier la caractéristique non linéaire du système. La sortie de la table de correspondance inversée est convertie en un signal analogique de type 4-20 mA I qui est fourni par le dispositif de régulation 10 comme valeur de consigne pour l'alimentation du courant d'arc de la source d'ions. [0040] Des simulations montrent qu'un tel dispositif permet une bonne .régulation. Il est cependant sensible aux perturbations basse fréquence. Pour résoudre ce problème, on a développé une variante préférée du dispositif suivant l'invention, représenté à la figure 4. Dans ce dispositif 10, on introduit dans la contre-réaction un filtre passe-bas 60 et un régulateur à avance de phase 70. Le filtre 60 est par exemple un filtre passe-bas du premier ordre . La fréquence de coupure est de 4,5 kHz . Afin de
compenser le retard de phase introduit par le filtre, on utilise un régulateur à avance de phase 70 (dérivateur filtré), qui compense ce déphasage.
[0041] Tant le dispositif de la figure 3 que celui de la figure 4 comportent une table de correspondance inversée 40. Le contenu de cette table 40 est déterminé préalablement à chaque utilisation du dispositif de la manière suivante :
- le régulateur étant en boucle ouverte, la consigne du courant d'arc de la source d'ions 20 est portée progressivement de 0 à 20 mA sous la forme d'une rampe de 100 ms;
- l'intensité du faisceau est mesurée pour chacun des 4000 points échantillonnés; - la table obtenue est inversée, de manière à fournir une valeur correspondante du courant d'arc de la source d'ions I A, en fonction de l'intensité du faisceau IM„.
- Cette table inversée est chargée dans le dispositif de régulation 10. [0042] En pratique, cette opération est réalisée une douzaine de fois successivement. Ceci permet de s'assurer que les paramètres atteignent un palier, correspondant à la température de régime du filament. Afin d'éliminer le bruit, une moyenne des 4 dernières tables est calculée. Ces opérations, réalisées automatiquement, durent au maximum 1,5 s. Dans une variante de l'invention, les valeurs de I correspondant aux valeurs de I supérieures à une limite donnée sont remplacées par la valeur de I correspondant à cette limite. Les courbes de la figure 2 sont donc écrêtëes. Ceci est un élément de sécurité permettant de garantir que l'intensité du faisceau produit par l'accélérateur ne sera jamais supérieur à cette limite.
[0043] Le dispositif suivant l'invention est réalisé au moyen d'une carte électronique qui fait appel aux technologies digitales de type DSP (Digital Signal Processing) . [0044] La synthèse du prédicteur de Smith à été réalisée dans le domaine de Laplace et la discrétisation est fournie par la transformée en Z par la méthode de correspondance pôles-zéros. Un sur-échantillonnage aurait été adéquat pour éviter tout problème lié à la discrétisation mais les technologies DSP du moment ne nous ont pas permis de monter au-delà de 100 kHz.
[0045] La méthode de régulation selon la présente invention présente plusieurs avantages. Tout d'abord, elle permet une adaptation commandée, c'est-à-dire qu'elle demande un temps de calcul très petit en comparaison des méthodes modernes du contrôle adaptatif et permet un changement de structure très facile puisque l'identification est faite par construction d'une table de correspondance qu'il suffit alors d'inverser numériquement pour linéariser la caractéristique du système vu par le régulateur principal .
[0046] En outre, elle offre une flexibilité importante puisqu'elle pourrait trouver sa place dans la régulation précise, reproductible, robuste et performante de toute source d'ions équipant un cyclotron et ce de par l'avantage d'une régulation de type adaptatif permettant la ré-identification de la caractéristique du système lorsque celle-ci évolue dans le temps. Elle permet donc l'identification et la régulation d'un autre accélérateur que le cyclotron C235 pour laquelle cette régulation à été développée à la base.