FR3060252A1 - Systemes de production de radionucleides - Google Patents

Systemes de production de radionucleides Download PDF

Info

Publication number
FR3060252A1
FR3060252A1 FR1771316A FR1771316A FR3060252A1 FR 3060252 A1 FR3060252 A1 FR 3060252A1 FR 1771316 A FR1771316 A FR 1771316A FR 1771316 A FR1771316 A FR 1771316A FR 3060252 A1 FR3060252 A1 FR 3060252A1
Authority
FR
France
Prior art keywords
time
constraints
radioactivity
objective function
distribution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
FR1771316A
Other languages
English (en)
Other versions
FR3060252B1 (fr
Inventor
Alexis Motto
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Medical Solutions USA Inc
Original Assignee
Siemens Medical Solutions USA Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens Medical Solutions USA Inc filed Critical Siemens Medical Solutions USA Inc
Publication of FR3060252A1 publication Critical patent/FR3060252A1/fr
Application granted granted Critical
Publication of FR3060252B1 publication Critical patent/FR3060252B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/67—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for remote operation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0633—Workflow analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/08—Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/083—Shipping
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/20—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the management or administration of healthcare resources or facilities, e.g. managing hospital staff or surgery rooms
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21G—CONVERSION OF CHEMICAL ELEMENTS; RADIOACTIVE SOURCES
    • G21G1/00—Arrangements for converting chemical elements by electromagnetic radiation, corpuscular radiation or particle bombardment, e.g. producing radioactive isotopes
    • G21G1/0005—Isotope delivery systems
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21G—CONVERSION OF CHEMICAL ELEMENTS; RADIOACTIVE SOURCES
    • G21G1/00—Arrangements for converting chemical elements by electromagnetic radiation, corpuscular radiation or particle bombardment, e.g. producing radioactive isotopes
    • G21G1/04—Arrangements for converting chemical elements by electromagnetic radiation, corpuscular radiation or particle bombardment, e.g. producing radioactive isotopes outside nuclear reactors or particle accelerators
    • G21G1/10—Arrangements for converting chemical elements by electromagnetic radiation, corpuscular radiation or particle bombardment, e.g. producing radioactive isotopes outside nuclear reactors or particle accelerators by bombardment with electrically charged particles
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Mathematical Physics (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Medical Informatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Marketing (AREA)
  • Public Health (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Biomedical Technology (AREA)
  • Tourism & Hospitality (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Data Mining & Analysis (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Game Theory and Decision Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • General Chemical & Material Sciences (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Educational Administration (AREA)
  • Nuclear Medicine (AREA)
  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)

Abstract

Système, qui comprend un cyclotron ; un module de chimie ; un système de dispensation de dose et un processeur 300 pour : déterminer une fonction d'objectif représentant des coûts de création des doses et de distribution des doses à des clients, la fonction d'objectif étant soumise à des contraintes de production, à des contraintes de distribution et à des contraintes de client ; donner une solution à la fonction d'objectif et produire un plan de production et un schéma de distribution.

Description

® RÉPUBLIQUE FRANÇAISE
INSTITUT NATIONAL DE LA PROPRIÉTÉ INDUSTRIELLE © N° de publication : 3 060 252 (à n’utiliser que pour les commandes de reproduction)
©) N° d’enregistrement national : 1771316
COURBEVOIE © Int Cl8 : H 05 H 13/00 (2017.01), G 01 R 33/64
DEMANDE DE BREVET D'INVENTION A1
©) Date de dépôt : 06.12.17. © Demandeur(s) : SIEMENS MEDICAL SOLUTIONS
© Priorité : 14.12.16 US 15378253. USA, INC. — US.
@ Inventeur(s) : MOTTO ALEXIS.
©) Date de mise à la disposition du public de la
demande : 15.06.18 Bulletin 18/24.
©) Liste des documents cités dans le rapport de
recherche préliminaire : Ce dernier n'a pas été
établi à la date de publication de la demande.
(© Références à d’autres documents nationaux ® Titulaire(s) : SIEMENS MEDICAL SOLUTIONS USA,
apparentés : INC..
©) Demande(s) d’extension : © Mandataire(s) : CABINET FLECHNER.
SYSTEMES DE PRODUCTION DE RADIONUCLEIDES.
FR 3 060 252 - A1
Système, qui comprend un cyclotron; un module de chimie; un système de dispensation de dose et un processeur 300 pour: déterminer une fonction d'objectif représentant des coûts de création des doses et de distribution des doses à des clients, la fonction d'objectif étant soumise à des contraintes de production, à des contraintes de distribution et à des contraintes de client; donner une solution à la fonction d'objectif et produire un plan de production et un schéma de distribution.
SYSTÈMES DE PRODUCTION DE RADIONUCLÉIDES
ARRIÈRE-PLAN [0001] L'imagerie médicale nucléaire utilise des composés radioactifs à bas niveau, ou radionucléides, pour produire des images de volumes internes d'un patient. Plus précisément, des radionucléides sont introduits dans un patient et absorbés par des organes, des os ou des tissus. Un scanner détecte des rayons gamma émis par les radionucléides absorbés et une image est produite sur la base des rayons gamma détectés.
[0002] Des doses de radionucléides sont produites habituellement par une radiopharmacie et distribuées à des clients (par exemple, un centre d'imagerie médicale nucléaire) suivant des termes définis à l'avance d'un contrat de client. Habituellement, un optimiseur de distribution détermine un schéma de distribution de toutes les doses reposant sur des paramètres de distribution connus et sur la radioactivité convenue et les temps de préparation de lot de dose et des doses individuelles de patient affectées à chaque lot. Ce qui sort de l'optimiseur de distribution permet de mobiliser des scénarios de production différents et des schémas d'affectation de dose pour déterminer l'impact sur le schéma de distribution et le coût total de la distribution.
[0003] L'opération de production ci-dessus est inefficace et il y a donc un besoin de systèmes procurant une intégration améliorée de la production et de la distribution des radionucléides.
L'invention a pour objet un système, caractérisé en ce qu'il comprend un cyclotron pour bombarder une cible, comprenant un premier matériau, par un faisceau de protons pour produire un deuxième matériau ; un système de module de chimie pour synthétiser le deuxième matériau en un marqueur radioactif de tomographie par émission de positrons ; un système de dispensation de dose pour préparer des doses du marqueur radioactif de tomographie par émission de positrons ; un dispositif d'entrée pour recevoir des valeurs de paramètre d'un opérateur et un processeur pour déterminer une fonction d'objectif représentant des coûts de création des doses et de distribution des doses à des clients, la fonction d'objectif étant soumise à des contraintes de production, à des contraintes de distribution et à des contraintes de client, les contraintes de production, les contraintes de distribution et les contraintes de client étant définies par des variables discrètes, des variables continues et des valeurs de paramètre ; recevoir les valeurs de paramètre du dispositif d'entrée ; donner une solution à la fonction d'objectif sur la base des valeurs paramètre reçues pour de des qui déterminer des valeurs variables discrètes et des variables continues, correspondent à des coûts appropriés de création des doses et de distribution des doses aux clients et produire un plan de production de marqueur radioactif de tomographie par émission de positrons et un schéma de distribution sur la base des valeurs des variables discrètes et des valeurs continues.
De préférence :
- Le processeur transforme, en outre, la fonction d'objectif en un modèle linéaire en nombre entier mixte.
- La transformation de la fonction d'objectif en un modèle linéaire en nombre entier mixte comprend la transformation de contraintes non-linéaires de la fonction d'objectif en des contraintes linéaires en nombre entier mixte stochastique à plusieurs degrés.
- Donner une solution à la fonction d'objectif comprend entrer le modèle linéaire en nombre entier mixte dans une application de résolution de programme en nombre entier mixte.
- Donner une solution à la fonction d'objectif repose sur une structure creuse du modèle linéaire en nombre entier mixte.
- Les contraintes de production comprennent une durée de bombardement, un temps de disponibilité du cyclotron, un temps d'indisponibilité du cyclotron, une radioactivité du premier matériau, un temps de déchargement de la cible, un temps de chargement de la cible, un temps de préparation du système de module de chimie, un temps de synthèse, un rendement en pourcentage et une radioactivité du marqueur radioactif de tomographie par émission de positrons à une fin de synthèse.
- Les contraintes de client comprennent une durée de distribution, une radioactivité ordonnée suivant les clients à un temps de distribution et une radioactivité totale de doses non-injectées à un site de client.
- Les contraintes de production comprennent une durée de bombardement, un temps de disponibilité du cyclotron, un temps d'indisponibilité du cyclotron, une radioactivité du premier matériau, un temps de déchargement de la cible, un temps de chargement de la cible, un temps de préparation du système de module de chimie, un temps de synthèse, un rendement en pourcentage et une radioactivité du marqueur radioactif de tomographie par émission de positrons à une fin de synthèse et dans lequel les contraintes de client comprennent une durée de distribution, une radioactivité ordonnée suivant les clients à un temps de distribution et une radioactivité totale de doses non-injectées à un site de client.
L'invention vise également un système comprenant une mémoire mémorisant des stades de traitement exécutables sur processeur et une unité de traitement pour exécuter les stades de traitement exécutables sur processeur pour déterminer une fonction d'objectif représentant des coûts de création des doses et de distribution des doses à des clients, la fonction d'objectif étant soumise à des contraintes de production, à des contraintes de distribution et à des contraintes de client, les contraintes de production, les contraintes de distribution et les contraintes de client étant définies par des variables discrètes, des variables continues et des valeurs de paramètre ; recevoir les valeurs de paramètre ; donner une solution à la fonction d'objectif sur la base d'un jeu de paramètres reçu pour déterminer des valeurs des variables discrètes et des variables continues, qui correspondent à des coûts appropriés de création des doses et de distribution des doses aux clients et produire un plan de production de marqueur radioactif de tomographie par émission de positrons et un schéma de distribution sur la base des valeurs des variables discrètes et des valeurs continues.
DESCRIPTION SUCCINCTE DES DESSINS [0004] L'agencement et l'usage de modes de réalisation apparaîtront facilement en considérant la description qui va suivre telle qu'illustrée dans les dessins annexés, dans lesquels des mêmes repères désignent des mêmes parties et dans lesquels :
[0005] la figure 1 illustre un système suivant certains modes de réalisation ;
[0006] la figure 2 illustre une opération de production 5 de radionucléide suivant certains modes de réalisation ;
[0007] la figure 3 est un diagramme de déroulement d'une opération suivant certains modes de réalisation ;
[0008] la figure 4 est un schéma synoptique fonctionnel d'un système suivant certains modes de réalisation ;
[0009] la figure 5 est un diagramme de déroulement d'une opération suivant certains modes de réalisation.
[0010] la figure 6 est un schéma synoptique fonctionnel d'un système suivant certains modes de réalisation et [0011] la figure 7 est un schéma synoptique d'un système informatique suivant certains modes de réalisation.
DESCRIPTION DÉTAILLÉE [0012] La description qui va suivre est donnée pour permettre à toute personne du métier de faire et d'utiliser les modes de réalisation décrits et indique le meilleur mode envisagé pour effectuer les modes de réalisation décrits. Mais, diverses modifications apparaîtront facilement à l'homme du métier.
[0013] La figure 1 illustre la production et la distribution de radionucléides suivant certains modes de réalisation. La production est régie par un plan de production et la distribution est régie par des schémas de distribution. Un plan de production est une collection de doses, qui sont affectées à un lot précis de produits.
Un lot individuel peut procurer suffisamment de produit pour jusqu'à un nombre fini de doses individuelles. Le nombre de doses de chaque lot dépend des localisations des clients et des ordres des clients. Le nombre de doses affectées à un lot dépend aussi de la dimension du lot, laquelle est régie à son tour par la quantité de radioactivité produite par les cyclotrons disponibles, par le rendement de la synthèse et par l'opération de contrôle de la qualité.
[0014] Un schéma de distribution est une collection de doses, qui sont affectées à un véhicule de distribution et qui sont envoyées aux lieux où se trouvent les clients. Les schémas de distribution sont spécifiques au site et reposent sur la géographie et sur des exigences locales de client. Les schémas de distribution peuvent varier en fonction du jour de la semaine et/ou des ordres saisonniers. En outre, les schémas de distribution peuvent nécessiter des changements de dernière minute pour faire face à des modifications des ordres (par exemple, des ordres nouveaux, des annulations ou une modification d'ordres antérieurs). D'autres contraintes associées à des schémas de distribution comprennent le nombre de véhicules de distribution qui sont disponibles, la capacité de chaque véhicule de distribution, etc.
[0015] Suivant certains modes de réalisation, un optimiseur de lot affecte des doses individuelles à des lots afin d'utiliser de la manière la plus efficace chaque lot. L'optimiseur de lot minimise simultanément les coûts associés à une production de lot, ainsi que les coûts associés à une distribution de dose. Les points de départ pour 1'optimisateur de lot sont une série de lots définis par l'utilisateur et une collection de doses individuelles au patient pour chaque lot, qui ont été affectées manuellement à chaque lot par l'utilisateur. Ce qui sort de l'optimiseur de lot est pareil à ce qui sort de l'optimiseur de distribution mentionné ci-dessus et en inclut également les affectations optimisées de dose à chaque lot, ainsi que la possibilité que l'utilisateur accepte les affectations optimisées.
[0016] Certains modes de réalisation prévoient un modèle de programmation stochastique à degré multiple, pour optimiser une chaîne de fourniture de radionucléides, qui s'accommode d'incertitudes de modèle, de contraintes opérationnelles douces et de contraintes physiques dures sur des variables de décision, et d'une approximation du modèle non-linéaire en un modèle linéaire en nombre entier mixte pour lequel il existe des dispositifs de résolution ouverts et dans le commerce. Les modèles sont écrits de manière à présenter leur structure creuse, qui peut être exploitée pour diminuer la complexité de calcul et l'encombrement d'une mémoire.
[0017] La figure 2 illustre une opération de production d'une radiopharmacie ou d'un autre vendeur de radionucléides. L'opération de production représente les blocs fonctionnels représentés à la figure 1 du cyclotron jusqu'à la dispensation de dose. Avant de passer en revue la figure 2, on définira ci-après les acronymes, paramètres, variables continues et variables discrètes mentionnés.
Acronymes
BG Groupe de lot.
BH Lot.
BOB Début du bombardement.
5 BOS début de la synthèse.
BR Branche.
CM Client.
CS Coffret.
CY Cyclotron.
10 DS Dose.
EOB Fin du bombardement.
EOU Fin du déchargement.
ETA Temps le plus tôt où la cible est disponible
MD Flacon multi-dose.
15 MS Site de fabrication.
ND Nœud.
PET Tomographie par émission de positrons.
PH Pharmacie.
QC Contrôle de qualité.
20 RN Radionucléide.
SY Système.
TN Temps.
TR Voyage.
UD Seringue à dose unitaire •
25 VH Véhicule.
Paramètres
TTBR “DliWfc Distance à parcourir entre l'origine et des nœuds de destination d'une branche f ;
TTBR “DUiWfc Durée de parcours entre 1 ' origine et des nœuds de destination d'une branche { ;
^NDcûjfc Nombre maximum de seringues à dose unitaire dans le coffret c ;
TTVH '‘•MRfùJk Vitesse kilométrique sous contrat avec un véhicule c
H-MAccûk Radioactivité maximum permise au site c du client ;
TTCY HSCCÙJl; Coût de démarrage du cyclotron c ;
^ETdWk Temps le plus tôt pour distribuer une dose d ;
^IAdiüfc Radioactivité d'injection d'une dose d ;
TTDS “iTdWfc Temps d'injection d'une dose d ;
TTDS 7TLTd(Jjk Dernier temps pour délivrer une dose d ;
TTND nSPnœk Redevance d'arrêt au site n ;
TTRN “DRrw/t Vitesse de décroissance d'un radionucléide r ;
7TSY “DTùjfc Echelon de temps discret ;
7TSY “FSüJfc Surcharge de carburant ;
7TSY reIVwk Nombre maximum de véhicules (internes) ;
SY n*Cœk Temps requis pour effectuer le contrôle de qualité ;
^DTtà),; Temps minimum d'inaccessibilité d'une cible t ;
4L, Temps minimum (ou le plus petit) de bombardement d'une cible t ;
Rendement de saturation pour une cible t ;
ttTG “UBtw,. Temps maximum (ou le plus grand) de bombardement
ttTG “VCtWfc d'une cible t ; Fonction de coût variable d'une cible t pour un
TTRN “ENridfc cyclotron c ; Enrichissement en radionucléide r.
Variables continues xRAbkù>k Radioactivité d'un lot b à un temps k ; xESbœk EOS d'un lot b ;
xRAc<j)k Radioactivité d'un cyclotron c à un temps k ; 5 xDTdœk Temps de distribution d'une dose d ;
DS xRA.dMk Radioactivité d'une dose d à EOB ;
<rDS ^ETdcdfc Fin d'un temps de bombardement EOB d'une dose d ;
„DS %RTd&)k Temps de préparation d'une dose d ;
TG xRAtku>k Radioactivité d'une cible t à un temps k ;
10 ~TG xBCtka>k Courant de faisceau appliqué à une cible t à un
temps k ;
<rTG xOFtka>k Nombre d'incréments temporels où la cible t n'a
pas été bombardée à la fin du temps k ;
rTG xONtka>k Nombre d'incréments temporels où la cible t a été
15 bombardée à la fin du temps k ;
<rTR ^CTtWfc Temps de fin du voyage t ;
tTR xDUiùjfc Durée du voyage t.
Variables discrètes
YvHcuwfc Egal à 1, si le coffret c est affecté au voyage t, ou sinon à 0 ;
CY ToNtfcw/; Egal à 1, si le cyclotron c est en fonctionnement au temps k, ou sinon à 0 ;
CY TsDtkWfc Egal à 1, si le cyclotron c est fermé au temps k, ou sinon à 0 ;
CY y STtkoifo Egal à 1, si le cyclotron c démarre au temps k, ou sinon à 0 ;
yc/\dcMk Egal à 1, si une dose d est affectée au coffret c, ou sinon à 0 ;
yCYdcwfc Egal à 1, si une dose d est affectée au cyclotron c, ou sinon à 0 ;
yETdkuk Egal à 1, si le EOB de la dose d est au temps k, ou sinon à 0 ;
SY yiVka>k Egal à 1, si tous les véhicules internes sont affectés à des voyages au temps k, ou sinon à 0 ; au début au temps k ;
yETtkù)k Egal à 1, si le bombardement se termine sur une cible c au début du temps k, ou sinon à 0 ;
yONtkù)k Egal à 1, si une cible t est en fonctionnement au temps k, ou sinon à 0 ;
ySOtkù)k Egal à 1, si une cible t est arrêtée au temps k, ou sinon à 0 ;
ySTtkWfc Egal à 1, si un bombardement a commencé sur la cible t au début du temps k, ou sinon à 0 ;
TR y PMtnka>k Egal à 1, si le voyage t arrive au nœud n au temps k, ou sinon à 0 ;
TR yBNtnk(iûk Egal à 1, si le voyage t arrive au nœud n au temps k ou sinon à 0 ;
VBEtlœk Egal à 1, si le voyage t traverse une branche T de
distribution, ou sinon à 0 ;
yvHtvwfc Egal à 1, si le voyage t est affecté à un véhicule d, ou sinon à 0 ;
yNDtnwfc Egal à 1, si le voyage t cesse au nœud n, ou sinon à 0 ;
[0018] De même, en ce qui concerne la notation utilisée ci-dessous, si x est une variable, x et x désignent sa limite inférieure et sa limite supérieure respectivement, si XY désignent une mnémonique à deux lettres de quelque
classe d'objet, JXY désigne le j eu d ' indice des objets
correspondants. Par exemple, ND est la mnémonique à deux
lettres de la classe de nœuds et 7nd désigne le jeu
d'indice des nœuds. Si i désigne un bord du réseau de distribution, o(-Z’) et d(-Z’) désignent le nœud d'origine et le nœud de destination respectivement de €. Si t désigne un voyage, o(t) et d(t) désignent le nœud d'origine et le nœud de destination respectivement de t. Le symbole /£M désigne le jeu d'indice de sites de clients fixes (c'està-dire k =1) et des sites de clients mobiles (c'est-àdire k = 2) . Le symbole 7™ désigne le jeu d'indice de temps accordé pour un voyage t, afin qu'il arrive au site c du client. Le symbole J^H désigne le jeu d'indice de véhicules internes (c'est-à-dire k = 1) et de véhicules extérieurs (c'est-à-dire k = 2) . Le symbole Ik s désigne le jeu de seringues à dose unique (c'est-à-dire k = 2) et de fioles à dose multiple (c'est-à-dire k = 2) . D'autres symboles ayant une portée plus étroite sont définis cidessous .
[0019] En revenant à la figure 2, la radioactivité d'un radionucléide au temps t, qt, est mis sous forme de modèle par l'équation différentielle du premier ordre suivante :
dqt = —Aqtdt, (1) [0020] dans laquelle λ désigne la constante de désintégration du radionucléide respectivement.
[0021] A partir de (1), on peut déterminer la radioactivité perdue par désintégration radioactive sur un intervalle de temps. Par exemple, sur la base de la radioactivité dans une seringue Qinj à dose unitaire au temps d'injection tINJ, la radioactivité à la fin de la synthèse qE0S est :
9inj = QeosexP(_T(tINJ - îeos)) (2 ) [0022] A partir de (2), on peut déduire qEOS de la manière suivante :
Qeos = qiNjeA(tlNJ“tEOs). (3) [0023] Opération de cyclotron [0024] Un cyclotron est un type d'accélérateur de particules ; c'est-à-dire un dispositif qui utilise des champs électromagnétiques pour propulser des particules chargées à de grandes vitesses et les contenir dans des faisceaux bien définis. Comme décrit à la figure 1, un cyclotron produit des radionucléides cibles ou des radioisotopes (par exemple, ion fluorure F) par bombardement par proton d'une cible contenant une matière première (par exemple, eau 18O) . Après déchargement des cibles bombardées d'un cyclotron, les cibles sont synthétisés chimiquement avec des molécules organiques (par exemple, du glucose) pour former un lot de marqueur radiopharmaceutique ou radioactif (par exemple, marqueur radioactif 18F-FDG). L'absorption de 18F-FDG par des tissus est un marqueur de l'absorption de glucose par le tissu, qui à son tour est en corrélation étroite avec certains types de métabolisme de tissu. Après avoir injecté 18F-FDG à un patient, un scanner PET peut former des images en deux dimensions ou en trois dimensions de la répartition de 18F-FDG au sein du corps.
[0025] Un équipement de fabrication peut comprendre un cyclotron ou plusieurs cyclotrons. Chaque cyclotron produit deux lignes de faisceau, de sorte qu'un cyclotron peut bombarder une cible ou plusieurs cibles. Il est ainsi possible d'utiliser n'importe quelle combinaison de cibles pour produire un lot de produit unique. Mais, si un cyclotron effectue un bombardement dual, les cibles peuvent être sur des lignes de faisceau différentes. Il n'est pas possible de bombarder simultanément deux cibles sur la même ligne de faisceau.
[0026] L'opération de chargement de cible de la figure 2 doit être achevée avant d'appliquer un faisceau à la cible. Si l'on utilise deux cibles dans l'opération de bombardement, l'opération de chargement de cible doit être achevée pour les deux cibles avant d'appliquer le faisceau à l'une des cibles. Les cibles peuvent être
chargées simultanément d' eau 18O. , Lorsque 1 ' on applique
en premier le faisceau à une cible, une opération
d'accord est nécessaire avant de pouvoir faire opérer le
faisceau à pleine puissance. Le début du bombardement (BOB) se produit après que le faisceau est à pleine puissance sur la cible. Si l'on utilise deux cibles dans l'opération de bombardement, le BOB peut être différent pour chaque cible. La radioactivité d'eau 180 et le temps nécessaire pour charger l'eau 180 peut différer pour chaque cible.
[0027] Le courant de faisceau est une mesure de flux de proton sur une cible et peut aller de 10 μΑ à 80 μΑ. La fin du bombardement (EOB) se produit lorsque le courant de faisceau est fixé à zéro. Si l'on utilise deux cibles dans l'opération de bombardement, le EOB sera le même pour les deux cibles. Le temps de bombardement est le laps de temps entre le BOB et le EOB. Le temps de bombardement peut aller de 5 minutes à 4 heures. Le rendement de saturation d'une cible représente la quantité théorique maximum d'ion fluorure 18F, qui peut être produite par une cible. Le rendement de saturation dépend de 1'énergie de faisceau du cyclotron et peut aller de 7 0 mCipA-1 à 2 60 mCipA-1.
[0028] La radioactivité de ion fluorure 18F produite dans une cible est donnée par :
qF = iBC aSY EN · (1 — e“ÂtB0MB) (4) ω/;(1 - ( 5 ) [0029] dans lequel qF désigne la radioactivité de l'ion fluorure 18F, iBC désigne le courant de faisceau, aSY désigne le rendement de saturation pour la cible, EN désigne l'enrichissement en l'eau 18O et tB0MB désigne le temps de bombardement.
[0030] Le temps de bombardement, le courant de faisceau, l'enrichissement en eau 18O et le rendement de saturation peuvent être différents pour chaque cible. Le temps de bombardement et le courant de faisceau peuvent être réglés lot par lot, afin de produire la radioactivité requise de l'ion fluorure 18F. Le rendement de saturation d'une cible ne change pas d'un lot à l'autre.
[0031] Les coûts de matériaux d'un lot dus à l'eau 18O augmentent avec la radioactivité de l'eau 18O utilisée dans la cible et avec le nombre de cibles utilisés dans le bombardement. Le coût dû à l'eau 18O dépend aussi de l'enrichissement. Le coût des matériaux n'est pas affecté par le courant de faisceau, le temps de bombardement ou le rendement de saturation. Le coût du travail augmente avec le temps de bombardement, en raison des coûts de travail indirects.
[0032] Suivant certains modes de réalisation, un optimiseur de lot détermine le nombre de cibles, le courant de faisceau, le BOB et le temps de bombardement pour chaque cible utilisée dans un lot.
[0033] Opération de déchargement de cible [0034] L'ion fluorure 18F produit par chaque cible est consolidé en une solution unique pendant l'opération de déchargement de la figure 2. Pour accomplir cela, chaque cible est déchargée séquentiellement dans un flacon de collecte (c'est-à-dire un flacon en v). Les cibles ne peuvent pas être déchargées simultanément. Le temps de déchargement peut être différent pour chaque cible. Le temps où le dernier déchargement de cible est achevé est défini comme la fin du déchargement (EOU). Après que le déchargement de la dernière cible est achevé, l'ion fluorure 18F total est donné par :
Qftotal = Qfi + qrie-^2 + qF3e~AtF3 + qF4e~^ ( 6 ) [0035] dans laquelle, Qftotal désigne la radioactivité totale de la désintégration de l'ion fluorure 18F corrigée à l'EOU, qF1 désigne la radioactivité de l'ion fluorure 18F produite dans le dernier déchargement de cible, qF2, qF3, qF4 représentent la radioactivité de l'ion fluorure 18F à la fin des déchargements de cible les plus tôt et tF2, tF3, tF4 représentent l'intervalle de temps entre un déchargement de cible le plus tôt et le EOU.
[0036] Opération de rechargement de cible en préparation d'un autre bombardement [0037] Après avoir achevé l'opération de déchargement, une cible peut être rechargée immédiatement en préparation d'un autre bombardement. Le temps nécessaire pour l'opération de rechargement est le même que le temps de l'opération initiale de chargement de cible. Après que la cible est rechargée, on peut commencer à faire arriver le faisceau sur la cible et le BOB suivant se produira après l'opération d'accord du faisceau. Cela représente le temps le plus tôt où la cible est disponible (ETA) pour le BOB suivant. Dans le cas d'un bombardement dual, le BOB suivant pour les deux cibles est le temps le dernier où l'une des cibles est disponible. On peut faire commencer le faisceau dès que les cibles pour le lot suivant ont été rechargées.
[0038] Opération de transfert de l'ion fluorure 18F [0039] Après que l'opération de déchargement de cible est achevée, la solution de l'ion fluorure 18F dans la fiole de collecte est transférée au module de chimie. Quel que soit le nombre de cibles, il y a seulement une opération de transfert pour chaque lot de produit. Un petit pourcentage de l'ion fluorure 18F total se perd pendant l'opération de transfert. Le temps de où l'opération de transfert de l'ion fluorure 18F est achevée est défini comme étant le début de la synthèse (BOS) . La radioactivité de l'ion fluorure 18F au BOS est donnée par :
Qbos = ^ftotai/I _ na)e ÂtTRAN (7) [0040] dans laquelle qB0S désigne la radioactivité de l'ion fluorure 18F au BOS, n désigne le nombre de cibles, a est le pourcentage en fraction perdu pendant l'opération de déchargement et tTRAN est l'intervalle de temps entre l'achèvement de l'opération de déchargement de la cible et le BOS.
[0041] Le pourcentage en fraction de la radioactivité perdue et le temps de transfert peuvent être généralisés comme étant les mêmes pour toutes les cibles à un équipement de fabrication.
[0042] Module de chimie [0043] Un équipement de fabrication peut avoir un module de chimie ou plusieurs modules de chimie. Un module de chimie peut être capable de traiter jusqu'à quatre lots de produit. Pour des modules de chimie ayant cette capacité, le module peut produire seulement des lots séquentiellement (c'est-à-dire qu'il n'est pas possible qu'un module de chimie prépare simultanément plus qu'un lot de produit).
[0044] Après que l'opération de transfert de l'ion fluorure 18F est achevée, un module de chimie transforme 1'ion fluorure 18F disponible en le médicament PET auquel on s'intéresse. Un pourcentage de l'ion fluorure 18F se perd pendant l'opération de chimie. Ce pourcentage est connu comme étant le rendement en pourcentage de l'opération de chimie. Le temps où l'opération de chimie est achevée est désigné comme étant la fin de la synthèse (EOS). La radioactivité du produit à 1’EOS est donnée par :
<7eos = qBosPYe“ÂtcHEM (8) [0045] dans laquelle <7eos désigne la radioactivité du produit à l'EOS, PY désigne le rendement en pourcentage de l'opération et tCHEM est l'intervalle de temps entre
1'EOS et le BOS (c'est-à-dire le temps de synthèse).
[0046] Chaque produit de médicament PET a un rendement en pourcentage unique et un temps de synthèse unique. Le rendement en pourcentage et le temps de synthèse peuvent varier à chaque équipement de fabrication. Le rendement en pourcentage et le temps de synthèse sont constants pour chaque lot de produit à un équipement de fabrication. Le coût des matériaux et du travail pour 1'opération de chimie comprend les coûts associés à la production et au test QC. Les matériaux et les coûts du travail sont constants pour chaque lot, mais sont différents pour chaque produit.
[0047] Nettoyage et préparation pour un autre lot [0048] Un module de chimie est associé à un temps de nettoyage, qui définit l'état jusqu'à ce que le module puisse être réutilisé pour un autre lot de produit. Le temps de nettoyage peut comprendre aussi un laps de temps nécessaire pour une désintégration radioactive de la radioactivité résiduelle, afin de permettre une intervention manuelle dans la préparation du lot suivant. Après expiration du temps de nettoyage, un module de chimie peut être utilisé pour préparer un autre lot de produit. Ce temps de nettoyage/préparation peut être unique pour chaque produit à un équipement de fabrication.
[0049] Test de contrôle de qualité [0050] Après que l'opération de chimie est achevée, le produit est échantillonné à des fins de test de contrôle de qualité (QC) . Cet échantillon est un pourcentage du produit provenant de l'opération de chimie. La radioactivité restante du produit est disponible pour une dispensation de dose. La radioactivité du produit disponible pour une dispensation de dose, désintégration corrigée de l'EOS, est donnée par :
Qdisp = (1 — ZOQeos (9) [0051] dans laquelle Qdisp désigne la radioactivité du produit disponible pour une dispensation de dose et β désigne le pourcentage du produit utilisé dans 1'échantillon.
[0052] Chaque médicament PET a un pourcentage d'échantillonnage unique et un temps d'échantillonnage unique. Le pourcentage d'échantillon et le temps d'échantillonnage peuvent varier à chaque équipement de fabrication. Le pourcentage d'échantillon et le temps d'échantillonnage sont constants pour chaque lot de produit à un équipement de fabrication.
[0053] Après que 1'échantillon a été retiré du produit, l'opération QC est effectuée et l'opération de dispensation de dose peut commencer sur le produit restant <7disp · Chaque médicament PET a un temps d'opération QC unique. Le temps d'opération QC peut varier à chaque équipement de fabrication. Le temps d'opération QC est constant pour chaque lot de produit à un équipement de fabrication.
[0054] Dispensation de dose [0055] L'opération de dispensation de dose peut être exécutée simultanément à l'opération QC. Les doses dispensées peuvent ne pas être libérées de l'équipement de fabrication jusqu'à ce que l'opération QC soit achevée et jusqu'à ce que le lot satisfasse tous les critères de libération QC. L'opération de dispensation de dose peut être utilisée pour préparer des flacons à dose multiple ou des seringues à dose unique en vue de leur répartition. La radioactivité agrégée, corrigée de la désintégration de toutes les doses dispensées ne peut pas excéder c/disp · Cette relation est donnée par :
Qdisp — Σ Qdose6 ^dose) (io) [0056] dans laquelle (/dose désigne la radioactivité du produit dans une dose individuelle au temps d’injection, toosE est l'intervalle de temps entre le temps d'injection et 1'EOS et la sommation s'effectue sur toutes les doses affectées au lot.
[0057] La figure 3 comprend un diagramme de déroulement d'une opération 300 suivant certains modes de réalisation. L'opération 300 peut être exécutée pour déterminer un plan de production et un schéma de distribution de radionucléides suivant un certain mode de réalisation. L'opération 300 et les autres opérations qui sont décrites ici peuvent être effectuées en utilisant toute combinaison appropriée de matériel, logiciel ou moyens manuels. Un logiciel incorporant ces opérations peut être mémorisé sur tout support tangible nontransitoire, y compris un disque fixe, un disque souple, un CD, un DVD, une mémoire flash ou une bande magnétique.
[0058] Des contraintes de production sont déterminées initialement à S310. Suivant certains modes de réalisation, les contraintes de production sont déterminées sur la base de paramètres, de variables continues et de variables discrètes comme énumérés cidessus . Les contraintes de production peuvent être mises sous forme de modèle pour supporter divers produits, chacun d'entre eux ayant une durée de demi-vie unique, le coût de fabrication et le temps de fabrication (par exemple, FDG, Amyvid, NaF, NH4, FLT et Lantheus). Suivant
certains modes de réalisation, les contraintes de
production peuvent comprendre des contraintes de
cyclotron, des contraintes de bloc de chimie, des
contraintes de bloc de contrôle de qualité et des
contraintes de dispensation de dose.
[0059] Contraintes de cyclotron [0060] Les contraintes de production peuvent comprendre des contraintes se rapportant directement au cyclotron ou aux cyclotrons utilisés pour bombarder la cible ou les cibles. Le nombre de cibles produites par campagne déterminent le débit et le coût d'un lot. Un cyclotron a des temps fixes minimum de disponibilité et d'indisponibilité, de sorte que la radioactivité de sortie d'un produit dépend de la durée de bombardement (c'est-à-dire l'intervalle de temps entre EOB et SOB).
[0061] Les contraintes (il) de production ci-dessous modélisent les contraintes de temps minimum de disponibilité et d'indisponibilité d'un cyclotron. La contrainte (lia) définit d'incréments temporels ou un cyclotron t a été en fonctionnement à la fin d'un intervalle de temps k. La contrainte (11b) définit le nombre d'incréments temporels où un cyclotron t n'a pas fonctionné à la fin d'un intervalle de temps k. Par définition, XoNtkœk est un nombre entier non-négatif et XoFtkcok est un nombre entier non-positif. Ainsi, le temps minimum d'indisponibilité est la valeur absolue de xOFtkœk · La contrainte (12a) définit la radioactivité d'une cible à partir de SOB et la contrainte (12b) indique que la radioactivité produite par un cyclotron est la somme de la radioactivité des cibles associées.
xONtko>k — (yONt(k-l)œk YsDt(k-l)Mk)xONt(k-l)œk + YoNtkcdk’ Vt ë 7TG, Vk 6 7™,Vü) G 7SC xoFtkœk = (1 - ySSt(k-i)cok - yïÎttk-^Ætk-Dc, + ySStkco, -1< vt e jtg, vk e j™, (x0Nt(k-l)Mk _ KUTtœk)(yONt(k-lj(,)k _ YoNtkcùk) θ’ Vt e Vk £ 7™, (x0Ft(k-l)œk ~ TÎDTtaik)(yONt(k-l)üjk _ YoNtkœk) -θ’ Vt G 7TG, Vk £ 7™, ySTtkcô, - yS TDtka>k = yjNtkœ, - ySNt(k-l)«k’ Vt £ 7TG, Vk £ 7™, yJw + ylotkrn, < 1’ vt £ 7tg, Vk e 7™,
..TG _ ,„TC» _TG „RN Γι _ „RN _SY „TG \] d'RAtkidk •t'BCtkülk Λ SYiWfc'* ENr(t)wfc L1 ^λΡ\, 71 DRr(t)wjfc n ) J (lia) (llb) (llc) (lld) (lie) (llf ) ,TG
Ittk„.,.TG ’-’WtkuJT, pvn/ RN SY TG λ ,,,ΤΟ Ί -AU /1 DRr(i)w;, /1 βΤω,. ·*- oF£àxua. ) J
Vieï^Vfcei™, (12a) xRAckœk = Zter/TG xRAtkf»)k’ Vc £ 7CY, (12b) xRAtkcok = xOFtktûkyONtkcûk + Ækco^Ntkc», ( 13 a ) [0062] Les contraintes suivantes équivalentes à (12a) peuvent être utilisées pour calculer xRAtkwk récursivement, la variable auxiliaire désignant la valeur de la fonction de désintégration au temps k :
xSLk = xtT(k-i)WkexP(- ^DRr(t)ü)k πϋΤω1<)(1 YsTtkœk) + ySTtkcùk (14a.) xRAtkwk = xBCtk(nk TÎSYt<»»k^ENr(t)ü3k(1 “ xüœk)YoNtkWk + (1 ~ Υ^ΝίΚω^ί^ (14b>
[0063] Contraintes du bloc de chimie [0064] Comme décrit ci-dessus, un bloc de chimie synthétise un lot de produit provenant des cibles produites au cyclotron. Un lot est donc une collection de cibles. La contrainte (15a) ci-dessous détermine la radioactivité d'un lot comme étant la somme d'activités sur toutes les cibles combinées. La contrainte (15b) assure que la radioactivité totale à un site de fabrication ne peut pas dépasser un seuil normalisé. Le nombre de boîtes de chimie est dépendant du site et dépendant du produit. Une unité de chimie a un cycle de fonctionnement minimum et le coût de la chimie est toujours constant par produit. Le temps nécessaire pour achever la chimie est constant par produit. Le plan de production peut comprendre une radioactivité recommandée à 1'EOS de chaque lot de produit.
[0065]
(15a) (15b) [0066] Contraintes du bloc de contrôle de qualité [0067] Le bloc de contrôle de qualité teste le produit radionucléide, ainsi que les matières premières et les fournitures utilisées pour la production du produit. Le bloc de contrôle de qualité soustrait une fraction, par exemple β, de la radioactivité d’un lot synthétisé. Pour simplifier les notations, à partir de maintenant, on écrira XRAbko>k P°ur (1 - β) XRAbkœk · [0068] Dispensation de dose [0069] Les lots qui passent avec succès les tests de contrôle de qualité sont dispensés en des flacons à dose multiple ordonnés suivant les clients ou en des seringues à dose unique en vue de la distribution. La radioactivité globale sur toutes les seringues de dose unitaire et tous les flacons à dose multiple dispensées à partir d’un lot donné est plus petite que la radioactivité du lot. La contrainte (16c) assure que le temps EOB d'une seringue u à dose unitaire coïncide avec un temps EOB du cyclotron c(u) affecté. La contrainte (16b) détermine le temps EOB de chaque dose et la contrainte (16b) assure que chaque dose est affectée à un lot.
(16a) (16b)
y.SDtCd)k<nk > YETdkœ,' Vd G 7DS, Vk G 7™, (160 xETdcok ~ Zke7™ yETdkü>kk' G 7DS, (16d) [0070] En retournant à l'opération 300, les contraintes de distribution sont déterminées à S320, après détermination des contraintes de production à S310. Les contraintes de distribution déterminées peuvent reposer aussi sur des paramètres, des variables continues et des variables discrètes.
[0071] Les contraintes de distribution suivant certains modes de réalisation relèvent de deux classes de coursiers : des coursiers extérieurs (par exemple, UPS, BDS, Carefree and MDS) et des coursiers internes. Par convention, chaque coursier peut servir un jeu donné d'emplacements de pharmacie et de client. Un véhicule sera associé à un itinéraire et distribuera à un nombre spécifié d'emplacements de client. Chaque véhicule a une charge utile maximum. Le temps maximum de EOS au temps d'injection est de 11 heures (par exemple, pour FDG). Les temps varieront suivant le produit. Chaque véhicule est autorisé à porter un niveau maximum de radioactivité dans chaque itinéraire (c'est-à-dire suivant sa licence de matériau radioactif (RAM)). Comme la radioactivité diminue en route, en raison de la décroissance radioactive, on compare ce niveau maximum au total de la radioactivité de chaque conteneur au début de l'itinéraire (par exemple, à l'équipement de production).
[0072] Les contraintes peuvent utiliser une moyenne (c'est-à-dire une moyenne de flacons par coffret) pour déterminer le nombre de coffrets à envoyer. Les coffrets doivent arriver 15 à 60 minutes avant le temps d'injection, le temps réel peut différer d'une convention à l'autre. Le modèle doit recommander la composition de chaque itinéraire sur la base de données de client, du temps de préparation de chaque cargaison et du temps de distribution de chaque cargaison. Dans le cas de certaines pharmacies, il est nécessaire que le véhicule revienne à la pharmacie pour décharger des coffrets vides et rapportér des clés. Il se peut que le véhicule doive arriver plus tôt à un emplacement de client (par exemple, pour satisfaire à des exigences de sécurité précises). Suivant l'emplacement, cette arrivée plus tôt peut varier entre 0 à X minutes et cela est pris en considération dans le calcul de la radioactivité des itinéraires les meilleurs et des temps de distribution. Certains modes de réalisation déterminent, en outre, des contraintes de distribution se rapportant au transport aérien.
[0073] La contrainte (17) ci-dessous se rapport à la cohérence de l'itinéraire, pour tout iG 7ND.
[0074]
' 1, si i = o(t), (17a)
- Σ EE, “ -1, si i = d(t), (17b)
fto(U = ,0, sinon (17c)
[0075] Un itinéraire n1 est pas défini par les tronçons
réels de route empruntés par un véhicule. En d'autres
termes un itinéraire peut avoir plus qu'une seule
possibilité de tronçon de route.
[0076] Un schéma de distribution est dit non-préemptif, si une fois qu'un véhicule a commencé un voyage, le véhicule doit achever le voyage avant de commencer un autre voyage. Un schéma est faisable seulement s'il est non-préemptif et chaque véhicule est affecté à au plus un voyage à un temps t.
[0077] La contrainte (18a) assure que chaque voyage est affecté à au plus un véhicule. La contrainte (18b) assure que le schéma de distribution n'est pas préemptif. Les contraintes (18c!
et assurent qu'un véhicule (18d) extérieur est affecté à un voyage seulement si un véhicule interne n'est pas disponible au temps k. En d'autres termes, les véhicules internes sont préférés.
ZveJVH YvHtvcûk — 1' £ 7TR, (18a)
TR TR > TR TR TR TR •''CTijW&.VDRijdw*. — •‘•CTifcWfci/DRffcdWfc 1 •‘'DUÎjWfci'DRfjrfw*.’ OU ,,,TR TR > TR TR , TR TR •''CTifcWn.yDRi^dwi,. — ••'crtjUi.yDMjdui,· r -MM,·ymdkdwi,.?
tkEîÎRy /tkyv eZVI1; (18b)
YlVkü>k “ ridef™ YvHtvaik- E 7TR, (18c)
YvHtvcok < YlVkœk< Vt e 7TR, Vv E 7™, (18d)
TR TR TR YANtnk<O]; ~ YBNtnk<jûk ~ YBNtn(k-l)a>k (18e)
[0078] Ensuite, à S330, on détermine des contraintes de client sur la base de paramètres, de variables continues et de variables discrètes. Il y a deux types de clients suivant certains modes de réalisation : (i) des clients fixes et (ii) des clients mobiles. Si 5™ et 32 M désignent les jeux d'indice de clients fixes et mobiles respectivement, la réunion de 7™ et de 7™ est 1® jen d'indice de tous les clients noté 7™. Un emplacement de client est fixé au temps où une dose est ordonnée. Un emplacement de client est ainsi associé à un nœud dans le réseau de la chaîne de fourniture. Conventionnellement, tous les clients doivent soumettre leur ordre de dose à 17 h pour une distribution le jour suivant. En pratique, des clients peuvent soumettre des ordres après 17 h pour une distribution le jour suivant ou même le même jour.
Zi „„SY \_DS <· „DS Avnr_^RN _SY Z_DS _ ,„DS V “iVwjJ'UArftvfc ~ 'EAcfeVfc ('XP[ 71 DR?’(d)wfc'‘DTWfe ('’lTf/u-'A,. KTdwjJ J < (1 + (19a) _DS < DS < _DS p 7DS πΕΤάω1< - XDTdtûk - KLTdœk< Va J ’ (19b)
Σ W* - 4L») - Ak - «j)«, «p[-4Sm- 4L»)) t - Ί ''+i · !,
Vceî^VfceZ :i9c) [0079] dans lesquelles u désigne la fonction Heaviside (ou stade unitaire). La contrainte (19a) assure que la différence absolue entre la radioactivité optimisée et ordonnée au temps d'injection ordonné d'une dose ne dépasse pas un seuil normalisé (par exemple 10%). La contrainte (19b) se rapporte à l'exigence du temps le plus tôt et le plus tard de distribution de la dose. La contrainte (19c) assure que la radioactivité de doses non-injectées à un site de client ne dépasse pas un seuil, comme déterminé par la licence RAM.
[0080] Les doses sont distribuées dans des récipients dénommés coffrets. Un coffret peut contenir 1, 2 ou 3 doses. Certains clients (par exemple, un personnel de centre d'imagerie) préfèrent des coffrets ayant une ou deux fioles, en raison de leur moindre poids et du personnel limité dont dispose le client. On fera l'hypothèse de la présence d'un nombre illimité de coffrets, qui seront utilisés pour distribuer les doses aux emplacements de client.
Ecejcs ycAdcWk = 1, Vd e JDS, (20a)
ZdC/DS yC DAdcWk *NDCWk, Vc £ 7CS, (20b)
Zvgjvh YvHcvcok = 1, Vc £ 7CS, (2 0c)
ΣνεΑΗ yvHcvœl( Æiccop Vc θ (20d) [0081] La contrainte (20a) assure que chaque dose est affectée à un coffret. La contrainte (20b) se rapporte au nombre maximum de doses qu'un coffret peut contenir. Chaque fiole à dose multiple est mise dans un coffret unique. La contrainte (20d) assure que chaque coffret est affecté à un voyage.
[0082] Les clients mobiles utilisent des coffrets de même poids, qui peuvent toutefois être d'une forme différente. Bien que ce ne soit pas habituellement le cas, il est possible que certaines restrictions de forme puissent être nécessaires. Les clients mobiles peuvent être traités comme des clients fixes. Chaque client mobile sera associé à un nombre d'envois sur la base de 1'emplacement où ils recevront des ordres.
[0083] On détermine une fonction d'objectif à S340 sur la base des contraintes de production, des contraintes de distribution et des contraintes de client, qui ont été déterminées. La fonction d'objectif consiste à minimiser les coûts totaux de production et de distribution Les coûts de distribution comprennent des d'arrêt, des coûts kilométriques, des surcharges de carburant et des péages. Les coûts de production englobent les coûts de mise en œuvre du cyclotron et les coûts de cible. Suivant certains modes de réalisation, en utilisant les contraintes déterminées comme décrit ci-dessus, une fonction d'objectif est la suivante :
escomptés. redevances [0084]
= E< | γ—Ά ,ç—\ i Zλi..... b·.·./? teb'P tel™ VH ^.BR TR MRWfc /l Dl&V/. »ORt(ùJk (21a)
+ Σ Σ Σ „ND TR SFWWfc (21b)
ί-εΧθΜ texTP r»exND
+ Σ Σ Σ r CY „TG I _TG f (21c)
fceï™ cexCY texTG )
[0085] Le terme dans (21a) représente le coût
kilométrique total, le terme dans (21b) représente la redevance d'arrêt total et le terme dans (21c) le coût total de production. Le coût total de production consiste en un coût fixe ou de démarrage (c'est-à-dire le premier terme dans (21c) ) et en un coût variable basé sur la radioactivité produite (c'est-à-dire le deuxième terme dans (21c)) .
[0086] Les péages sont inclus implicitement dans le modèle proposé par 1'intermédiaire des branches du réseau et de leurs coûts associés. On peut, par exemple, ajouter deux branches entre deux nœuds, une sans péage et 1'autre avec péage.
[0087] Le coût variable principal d'un stock de production est le coût du radionucléide 180. Le cyclotron est associé à un coût de mise en œuvre (à inclure dans le coût d'un lot) et les coûts de distribution sont : des redevances d'arrêt, des coûts kilométriques, des surcharges de carburant et des péages. Les coûts de véhicule interne seront calculés par l'intermédiaire des paramètres ci-dessus et représentés sous la forme d'un coursier. Il faut affecter des véhicules internes avant les véhicules extérieurs. La redevance d'arrêt est un coût contractuel fixe par client, par emplacement de distribution. Le coût kilométrique est une fonction linéaire déterminée par la vitesse associée et la distance réelle parcourue. Le taux kilométrique est déterminé par le coursier et indiqué dans la convention. Pour un itinéraire donné, la distance réelle parcourue par un véhicule de distribution est une variable aléatoire que l'on peut considérer comme une valeur moyenne sur la base de données antérieures. La surcharge de carburant dépend du coût de l'essence et est déterminée par le coût national en vigueur. La redevance d'arrêt est un coût fixe par emplacement de client à chaque arrêt. Le péage est un coût fixe par péage.
[0088] En conséquence, à S350, on détermine des valeurs de chaque variable discrète et continue de la fonction d'objectif, de manière à minimiser la fonction d'objectif sur la base de valeurs de paramètre données (c'est-à-dire connues). Les valeurs de variables déterminées constituent des variables d'un plan de production et d'un schéma de distribution.
[0089] La figure 4 est un schéma synoptique fonctionnel illustrant S350. suivant certains modes de réalisation. Comme représenté, la fonction 410 d'objectif comprend des contraintes de production, des contraintes de distribution et des contraintes de client. Chacune des contraintes est formulée sous la forme d'une combinaison ou de valeurs de paramètre, de variables continues et de variables discrètes. Afin de minimiser la fonction d'objectif (et de déterminer ainsi les valeurs de variables qui se traduisent par une minimisation de la fonction), le module 420 de minimisation évalue la fonction 410 d'objectif en utilisant des valeurs 430 de paramètres connues et des valeurs initiales des variables continues et discrètes des contraintes de production, des contraintes de distribution et des contraintes de client.
[0090] Ensuite, comme représenté à la figure 4, le module 420 modifie les valeurs de variables sur la base de l'évaluation de la fonction 410 d'objectif et réévalue la fonction 410 d'objectif sur la base des variables modifiées. L'opération se continue jusqu'à ce que le module 420 détermine que la fonction 410 d'objectif a été suffisamment minimisée. Des valeurs en cours des variables discrètes et continues sont alors utilisées pour produire le plan 440 de production et le schéma 450 de distribution.
[0091] Des systèmes pour mettre en oeuvre l'opération de minimisation du module 420 sont connus dans la technique. Des systèmes de ce genre peuvent présenter une complexité de calcul significative. C'est pourquoi, l'opération 500 de la figure 4 procure un système pour minimiser la fonction d'objectif et déterminer des valeurs de variable correspondantes d'une façon plus efficaces en informatique.
[0092] Les stades S510 à S540 peuvent se dérouler de la même façon que les stades S310 à S340 de l'opération 300 et c'est pourquoi on n'en répétera pas la description. Toutefois, on détermine un module linéaire en nombre entier mixte sur la base de la fonction d'objectif. La détermination du modèle linéaire en nombre entier mixte peut comprendre la transformation de contraintes nonlinéaires, telles que les contraintes de disponibilité, d'indisponibilité minimum et les contraintes de sortie de radioactivité maximum et minimum, en des contraintes linéaires en nombre entier mixte stochastique à degré multiple. En outre, on discrétise la contrainte de coût de démarrage non-linéaire en une fonction linéaire du temps de non-fonctionnement d'une cible. Un compteur linéaire de temps de non-fonctionnement est déduit en tant que sous-produit. Des fonctions de coût variable non convexe sont exprimées aussi sous la forme de fonctions linéaires en nombre entier mixte pièce par pièce.
[0093] Afin de linéariser des termes de produit, chaque monôme, dont le degré total par rapport à des variables continues est inférieur ou égal à un, est remplacé par un jeu de contraintes d'égalité ou d'inégalité linéaire en nombre intégré mixte. Par exemple, si x et y sont des variables scalaires continues, u, r et w sont des variables scalaires binaires et i, j, k, m et n des nombres entiers non-négatifs. Le degré total par rapport à des variables continues du monôme xlyiukvmwn est égal à i+j.
[0094] Par exemple, la contrainte en nombre entier mixte z = xy, où x est une variable continue limitée et y une variable binaire, est équivalente aux contraintes linéaires en nombre entier mixte suivantes :
y x < z < yx, (22a) (1 — y)x < x — z < (1 - y)x. (22b) [0095] La contrainte (22a) assure que z est égal à zéro si y est égal à zéro. La contrainte (22b) assure que z est égal à x si y est égal à 1. La contrainte en nombre entier mixte z = xy, où à la fois x et y sont d'une variable binaire, est équivalente aux contraintes linéaires en nombre entier mixte suivantes :
z < x, (23a) z<y, (23b) z > x + y - 1. (23c) [0096] Les contraintes (23a) te (23b) assurent que z est égal zéro si y est égal à zéro ou y est égal à zéro. La contrainte (23c) assure que est égal à 1 si à la fois x et y sont égaux à 0.
[0097] Les contraintes (11) en nombre entier mixte peuvent être transformées en la forme linéaire en nombre entier mixte équivalente suivante :
Σ,ΪΓ [l - ySNtkœJ = θ (24a)
k+K^tù)k- LT=k 1 ySiZk > πϊΖ,ΥίτΖ,. k = + 1,... ,|5™ | - + 1 (24b)
v|7™| ZjT=k [ySZMk - yJTtkcJ 0, k = |7™ | - + 2.....P™| (24c)
OÙ πυτΐωι( — min{|7™|, (TCuTtü)k xONtOtôk)yONtOwk}
YoNtkWk = 0 (25a) 1 [l - yc)NtT«kJ TtDTtcûkysDtk(1>k' k = TtS^tü)k + 1.....| J™| - ir$tù)k + 1 (25b)
Στ™' [1 ΥοΝίτω, YsDtkœ,] > 0, k = |J™| - + 2.....|7™ | (250 où îï^tWk = min{|J™|, (irj^ - x5^0Wk)(l - y5NtoWk)} [0098] La contrainte (24a) se rapporte aux contraintes de disponibilité minimum de chaque cible t qui a été en fonctionnement à une période 0 pendant une période de temps n'excédant pas leur disponibilité minimum. La contrainte (24b) se rapporte à la logique de disponibilité minimum pour tous les jeux de périodes consécutives de dimension ττυτίωι< - Les contraintes (24c) se rapportent à la logique de temps de disponibilité minimum pour les dernières périodes TTuTtwk “ 1 i c'est-à-dire si une cible t est démarrée dans l'une de ses heures, elle doit rester en fonctionnement jusqu'à la fin de l'horizon de planification. Les contraintes (25a) à (25c) , qui se rapportent à la logique d'indisponibilité minimum sont déduite de (24a) à (24c) par substitution de 1— yoNtk<ok' ySTtkœk' πυτίω1{ et xONtOcok Pour yoNtkœk> YDTtkœk’ πϋΤίω1< et xOFtO(ok respectivement.
[0099] En ce qui concerne la production de radioactivité et la désintégration, l’approximation linéaire en nombre entier mixte suivante peut être utilisée pour représenter des cas éventuels, dans lesquels un courant de faisceau variable peut être appliqué à une cible :
x^ - x£2tl«ok < (l - yS!Uk) + ^toky0 TGNtkü)k, vt e JTG, vk e J™, (26a) x^tko>k - x^t(k-i)<ok < π^1ωκ(ΐ - ysTXk) + ÎTRstWkysTtkMk· vt e vk e 3™ (26b) xRMkœk>0< Vt£7TG,Vke7™ (26c) [0100] Les contraintes de production mixte suivantes se rapportent à la faisibilité des états en fonctionnement et hors de fonctionnement et des transitions d'état (par exemple une cible ne peut pas être simultanément en fonctionnement et hors de fonctionnement).
YsTtkcôk - YSDtkWk - yONtkœk + Æ(k-l)Cok = θ- Vt £ 1TG, Vk G 7™, (27a)
yÎTtkcûk + yJntkcok - 1 < 0, vt e JTG, vk e J™. (27b)
[0101] Les contraintes suivantes définissent le nombre
d1 incréments temporels ou une cible a été en
fonctionnement :
TG < TG , i •'ΌΝ/Ά-ωγ. — ' i ΐ.ςγ, ' L (28a)
+ (XOîïtku,k +1)(1- î/oNtto*) æONi(fc-lM + (28b)
r TG ..... yTG TG < q AONtfcWA.i/ONtto*. — υ’ (28c)
> 0. (28d)
[0102] Les contraintes (28) définissent que le compteur
de temps en fonctionnement d'une cible t au temps k est
incrémenté de 1 si la cible t est en fonctionnement au temps k ; sinon, le compteur de temps est fixé à zéro.
[0103] Les contraintes suivantes définissent le nombre d'incréments temporels où une cible a été hors de fonctionnement :
UTG < tg , i — ^or-tlk-Putk + : ,,,TG , /-.TG , -, y. TG -> ,TG i ' oi'tfcu.’/, ' L'oit,/ru,/, t -Ua/onîA-wî,. — 1 ·*·>
,,TG _¥TG /-, _ ,,.TG \<n d‘OVtkuk '‘ΌΜΐίωιλ1 yoMltWk' — /•το > n «^OFiAwfc u!
(29a) (29b) (29c) (29d) [0104] Les contraintes (29) définissent que le compteur de temps de non-fonctionnement d’une cible t à un temps k est incrémenté de 1 si la cible t n’est pas en fonctionnement au temps k ; sinon le compteur de temps est fixé à zéro.
[0105] On résout le modèle linéaire en nombre entier mixte à S560. On effectue la résolution en se fondant sur des valeurs de paramètres connues pour déterminer des valeurs des variables de contraintes, qui se traduisent par un modèle linéaire en nombre entier mixte minimisé (c'est-à-dire minimisé à un degré acceptable). On utilise ensuite les variables pour déterminer un plan de production et un schéma de distribution, comme décrit cidessus. Avantageusement, une instance de la relaxation linéaire de modèle linéaire en nombre entier mixte peut être résolue en utilisant un procédé de programmation dynamique dual stochastique procuré par des dispositifs de résolution que l'on peut se procurer dans le commerce (par exemple CPLEX, Gurobi). Ces dispositifs de résolution peuvent avoir un effet de levier sur la structure creuse du modèle linéaire en nombre entier mixte pour procurer un calcul plus efficace, qui ne serait pas possible sinon.
[0106] La figure 6 est un diagramme synoptique fonctionnel illustrant l'opération 500 suivant certains modes de réalisation. Comme représenté, la fonction 610 d'objectif comprend des contraintes de production, des contraintes de distribution et des contraintes de client. On crée le modèle 620 linéaire en nombre entier mixte sur la base de la fonction 610 d'objectif à S550 et on l'entre dans le dispositif 630 de résolution à S560. Le dispositif 630 de résolution détermine ensuite les valeurs des variables discrètes et continues du modèle 620 linéaire en nombre entier mixte, qui minimise les coûts de production et de distribution. On utilise ensuite ces valeurs pour créer le plan 640 de production et de schéma 650 de distribution.
[0107] La figure 7 illustre un système 700 suivant certains modes de réalisation. Le système 700 peut exécuter n'importe laquelle des opérations décrites dans le présent mémoire suivant certains modes de réalisation. Le système 7 00 comprend un ou des unité (s) 710 de traitement et un dispositif 720 de mémorisation, chacun d'entre eux pouvant comprendre toute unité de traitement pour exécuter un code de programme exécutable par un processeur et tout système pour mémoriser une donnée électronique, qui est connue ou qui le devient.
[0108] Le dispositif 720 de mémorisation mémorise un code de programme exécutable par processeur d'un optimiseur. L'optimiseur peut déterminer un schéma de distribution et un plan de production sur la base de valeurs de paramètre décrites dans le présent mémoire. La détermination repose sur les contraintes de production, les contraintes de distribution, les contraintes de client et la fonction d'objectif mentionnées ci-dessus. Les valeurs de paramètre peuvent être mémorisées dans le dispositif 720 de mémorisation comme représenté et/ou peuvent être reçues d'un autre dispositif informatique ou d'un terminal 730 du système 700.
[0109] La détermination du plan de production et du schéma de distribution peut reposer sur un modèle linéaire en nombre entier mixte comme décrit ci-dessus à
propos de 1 'opération 500. Dans un cas de ce genre,
1'optimiseur peut utiliser un code de programme
exécutable par processeur d'un dispositif de résolution
MILP mémorisé également dans le dispositif 720 de mémorisation.
[0110] Le système 700 peut comprendre toute combinaison de matériel informatique, qui est connue ou qui le deviendra. Dans certains modes de réalisation, le système
700 est un système informatique de bureau, qui effectue toutes les déterminations décrites dans le présent mémoire. Dans d'autres modes de réalisation, le système 700 comprend une architecture client-serveur, dans laquelle le terminal 730 fonctionne principalement comme une interface utilisateur et les déterminations sont effectuées par un serveur éloigné. Le système 700 peut comprendre d'autres éléments, qui sont nécessaires pour son fonctionnement, ainsi que des éléments supplémentaires pour procurer des fonctions autres que celles décrites dans le présent mémoire.
[0111] L'homme du métier comprendra que diverses adaptations et modifications des modes de réalisation décrits ci-dessus peuvent être configurées sans sortir de l'esprit de la présente invention. Il va ainsi de soi que l'invention peut être mise en pratique d'une manière autre que décrite précisément dans la présent mémoire.
LÉGENDE DE LA FIGURE 2
Dose dispensing : Dispensation de doses
Quality control testing :
Test de contrôle de qualité
Chemistry Module de module : chimie
Cyclotron 1/2 : Cyclotron 1/2
Initialize cyclotron : Initialiser le cyclotron
Tune beam :
Accorder le faisceau
Load target 1/2/3/4 :
Charger la cible 1/2/3/4
Bombard target 1/2/3/4 : Bombarder la cible 1/2/3/4
Unload target 1/2/3/4 :
Décharger la cible 1/2/3/4
Reload target 1/2/3/4 :
Recharger la cible 1/2/3/4
Setup chemistry module :
Mettre en œuvre le module de chimie
Deliver collection vial : 25 Distribuer les flacons
Chemistry procèss :
Opération de chimie
Clearing & préparation process :
Opération de nettoyage et de préparation
Setup guality control process :
Mettre en œuvre l'opération de contrôle de qualité
Product sampling :
Echantillonnage de produit
Quality control testing :
Test de contrôle de qualité
Setup for dose dispensing :
Mise en œuvre de la dispensation de doses
Release :
Libération
Time :
Temps

Claims (20)

  1. REVENDICATIONS
    1. Système, caractérisé en ce qu'il comprend : un cyclotron pour bombarder une cible, comprenant un premier matériau, par un faisceau de protons pour produire un deuxième matériau ;
    un système de module de chimie pour synthétiser le deuxième matériau en un marqueur radioactif de tomographie par émission de positrons ;
    un système de dispensation de dose pour préparer des doses du marqueur radioactif de tomographie par émission de positrons ;
    un dispositif d'entrée pour recevoir des valeurs de paramètre d'un opérateur et un processeur pour :
    déterminer une fonction d'objectif représentant des coûts de création des doses et de distribution des doses à des clients, la fonction d'objectif étant soumise à des contraintes de production, à des contraintes de distribution et à des contraintes de client, les contraintes de production, les contraintes de distribution et les contraintes de client étant définies par des variables discrètes, des variables continues et des valeurs de paramètre ;
    recevoir les valeurs de paramètre du dispositif d'entrée ;
    donner une solution à la fonction d'objectif sur la base des valeurs de paramètre reçues pour déterminer des valeurs des variables discrètes et des variables continues, qui correspondent à des coûts appropriés de création des doses et de distribution des doses aux clients et produire un plan de production de marqueur radioactif de tomographie par émission de positrons et un schéma de distribution sur la base des valeurs des variables discrètes et des valeurs continues.
  2. 2. Système suivant la revendication 1, le processeur transformant, en outre :
    la fonction d'objectif en un modèle linéaire en nombre entier mixte.
  3. 3. Système suivant la revendication 2, dans lequel la transformation de la fonction d'objectif en un modèle linéaire en nombre entier mixte comprend la transformation de contraintes non-linéaires de la fonction d'objectif en des contraintes linéaires en nombre entier mixte stochastique à plusieurs degrés.
  4. 4. Système suivant la revendication 2, dans lequel donner une solution à la fonction d'objectif comprend entrer le modèle linéaire en nombre entier mixte dans une application de résolution de programme en nombre entier mixte.
  5. 5. Système suivant la revendication 2, dans lequel donner une solution à la fonction d'objectif repose sur une structure creuse du modèle linéaire en nombre entier mixte.
  6. 6. Système suivant la revendication 2, dans lequel les contraintes de production comprennent une durée de bombardement, un temps de disponibilité du cyclotron, un temps d'indisponibilité du cyclotron, une radioactivité du premier matériau, un temps de décharge de la cible, un temps de charge de la cible, un temps de préparation du système de module de chimie, un temps de synthèse, un rendement en pourcentage et une radioactivité du marqueur radioactif de tomographie par émission de positrons à une fin de synthèse.
  7. 7. Système suivant la revendication 6, dans lequel les contraintes de client comprennent une durée de distribution, une radioactivité ordonnée suivant les clients à un temps de distribution et une radioactivité totale de doses non-injectées à un site de client.
  8. 8. Système suivant la revendication 1, dans lequel les contraintes de production comprennent une durée de bombardement, un temps de disponibilité du cyclotron, un temps d'indisponibilité du cyclotron, une radioactivité du premier matériau, un temps de décharge de la cible, un temps de charge de la cible, un temps de préparation du système de module de chimie, un temps de synthèse, un rendement en pourcentage et une radioactivité du marqueur radioactif de tomographie par émission de positrons à une fin de synthèse et dans lequel les contraintes de client comprennent une durée de distribution, une radioactivité ordonnée suivant les clients à un temps de distribution et une radioactivité totale de doses non-injectées à un site de client.
  9. 9. Système comprenant :
    une mémoire préconisant des stades de traitement exécutables sur processeur et une unité de traitement pour exécuter les stades de traitement exécutables sur processeur pour :
    déterminer une fonction d'objectif représentant des coûts de création des doses et de distribution des doses à des clients, la fonction d'objectif étant soumise à des contraintes de production, à des contraintes de distribution et à des contraintes de client, les contraintes de production, les contraintes de distribution et les contraintes de client étant définies par des variables discrètes, des variables continues et des valeurs de paramètre ;
    recevoir les valeurs de paramètre ;
    donner une solution à la fonction d'objectif sur la base d'un jeu de paramètre reçu pour déterminer des valeurs des variables discrètes et des variables continues, qui correspondent à des coûts appropriés de création des doses et de distribution des doses aux clients et produire un plan de production de marqueur radioactif de tomographie par émission de positrons et un schéma de distribution sur la base des valeurs des variables discrètes et des valeurs continues.
  10. 10. Système suivant la revendication 9, l'unité de traitement exécutant, en outre, les stades de traitement exécutable sur processeur pour :
    transformer la fonction d'objectif en un modèle linéaire en nombre entier mixte.
  11. 11. Système suivant la revendication 10, dans lequel la transformation de la fonction d'objectif en un modèle linéaire en nombre entier mixte comprend la transformation de contraintes non-linéaires de la fonction d'objectif en des contraintes linéaires en nombre entier mixte stochastique à plusieurs degrés.
  12. 12. Système suivant la revendication 10, dans lequel donner une solution à la fonction d'objectif comprend entrer le modèle linéaire en nombre entier mixte dans une application de résolution de programme en nombre entier mixte.
  13. 13. Système suivant la revendication 10, dans lequel donner une solution à la fonction d'objectif repose sur une structure creuse du modèle linéaire en nombre entier mixte.
  14. 14. Système suivant la revendication 10, dans lequel les contraintes de production comprennent une durée de bombardement, un temps de disponibilité du cyclotron, un temps d'indisponibilité du cyclotron, une radioactivité du premier matériau, un temps de décharge de la cible, un temps de charge de la cible, un temps de préparation du système de module de chimie, un temps de synthèse, un rendement en pourcentage et une radioactivité du marqueur radioactif de tomographie par émission de positrons à une fin de synthèse, et dans lequel les contraintes de client comprennent une durée de distribution, une radioactivité ordonnée suivant les clients à un temps de distribution et une radioactivité totale de doses non-injectées à un site de client.
  15. 15. Procédé mis en œuvre par ordinateur comprenant :
    déterminer une fonction d'objectif représentant des coûts de création des doses et de distribution des doses à des clients, la fonction d'objectif étant soumise à des contraintes de production, à des contraintes de distribution et à des contraintes de client, les contraintes de production, les contraintes de distribution et les contraintes de client étant définies par des variables discrètes, des variables continues et des valeurs de paramètre ;
    recevoir un jeu de paramètre ;
    donner une solution à la fonction d'objectif sur la base du jeu de paramètre reçu pour déterminer des valeurs des variables discrètes et des variables continues, qui correspondent à des coûts appropriés de création des doses et de distribution des doses aux clients et produire un plan de production de marqueur radioactif de tomographie par émission de positrons et un schéma de distribution sur la base des valeurs des variables discrètes et des valeurs continues.
  16. 16. Procédé suivant la revendication 15, comprenant, en outre :
    la fonction d'objectif en un modèle linéaire en nombre entier mixte.
  17. 17. Procédé suivant la revendication 16, dans lequel la transformation de la fonction d'objectif en un modèle linéaire en nombre entier mixte comprend la transformation de contraintes non-linéaires de la fonction d'objectif en des contraintes linéaires en nombre entier mixte stochastique à plusieurs degrés.
  18. 18. Procédé suivant la revendication 17, dans lequel donner une solution à la fonction d'objectif repose sur une structure creuse du modèle linéaire en nombre entier mixte.
  19. 19. Procédé suivant la revendication 17, dans lequel les contraintes de production comprennent une durée de bombardement, un temps de disponibilité du cyclotron, un temps d'indisponibilité du cyclotron, une radioactivité du premier matériau, un temps de décharge de la cible, un temps de charge de la cible, un temps de préparation du système de module de chimie, un temps de synthèse, un rendement en pourcentage et une radioactivité du marqueur radioactif de tomographie par émission de positrons à une fin de synthèse, et dans lequel les contraintes de client comprennent une durée de distribution, une radioactivité ordonnée suivant les clients à un temps de distribution et une radioactivité totale de doses non-injectées à un site de client.
  20. 20. Système suivant la revendication 15, dans lequel les contraintes de production comprennent une durée de bombardement, un temps de disponibilité du cyclotron, un temps d'indisponibilité du cyclotron, une radioactivité du premier matériau, un temps de décharge de la cible, un temps de charge de la cible, un temps de préparation du système de module de chimie, un temps de synthèse, un rendement en pourcentage et une radioactivité du marqueur radioactif de tomographie par émission de positrons à une fin de synthèse, et dans lequel les contraintes de client comprennent une durée de distribution, une radioactivité ordonnée suivant les clients à un temps de distribution et une radioactivité totale de doses non-injectées à un site de client.
    1/7 ο
FR1771316A 2016-12-14 2017-12-06 Systemes de production de radionucleides Active FR3060252B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US15378253 2016-12-14
US15/378,253 US20180165420A1 (en) 2016-12-14 2016-12-14 Systems for Radionuclide Production

Publications (2)

Publication Number Publication Date
FR3060252A1 true FR3060252A1 (fr) 2018-06-15
FR3060252B1 FR3060252B1 (fr) 2020-12-11

Family

ID=60950113

Family Applications (1)

Application Number Title Priority Date Filing Date
FR1771316A Active FR3060252B1 (fr) 2016-12-14 2017-12-06 Systemes de production de radionucleides

Country Status (3)

Country Link
US (1) US20180165420A1 (fr)
FR (1) FR3060252B1 (fr)
GB (1) GB2558778A (fr)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003094107A2 (fr) * 2002-05-02 2003-11-13 Manugistics, Inc. Planification et ordonnancement de production fondes sur des contraintes
FR2867084A1 (fr) * 2004-03-02 2005-09-09 Gen Electric Systemes, procedes et appareils de perfusion de produits radiopharmaceutiques

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003094107A2 (fr) * 2002-05-02 2003-11-13 Manugistics, Inc. Planification et ordonnancement de production fondes sur des contraintes
FR2867084A1 (fr) * 2004-03-02 2005-09-09 Gen Electric Systemes, procedes et appareils de perfusion de produits radiopharmaceutiques

Also Published As

Publication number Publication date
GB201720029D0 (en) 2018-01-17
GB2558778A (en) 2018-07-18
US20180165420A1 (en) 2018-06-14
FR3060252B1 (fr) 2020-12-11

Similar Documents

Publication Publication Date Title
Tetzlaff et al. Identifying birth places of young isolated neutron stars
Asner et al. Averages of b-hadron, c-hadron, and tau-lepton Properties
Klein et al. Blind analysis in nuclear and particle physics
Kaspi et al. A massive radio pulsar binary in the Small Magellanic Cloud
Irrgang et al. Hypervelocity stars in the Gaia era-Runaway B stars beyond the velocity limit of classical ejection mechanisms
Rossi et al. Joint constraints on the Galactic dark matter halo and Galactic Centre from hypervelocity stars
Boubert et al. On the kinematics of a runaway Be star population
Southworth et al. A multiplicity study of transiting exoplanet host stars-II. Revised properties of transiting planetary systems with companions
Marchetti Gaia EDR3 in 6D: searching for unbound stars in the galaxy
Dupourqué et al. Constraints on the antistar fraction in the Solar System neighborhood from the 10-year Fermi Large Area Telescope gamma-ray source catalog
CN109377224A (zh) 区块链交易方法及装置、电子设备
CN111818095A (zh) 基于区块链的退税方法、装置及电子设备
CN110033370A (zh) 账户创建方法及装置、电子设备、存储介质
EP2738103A1 (fr) Procédé et système de mise à poste d&#39;un satellite
Yamaguchi et al. Wide post-common envelope binaries containing ultramassive white dwarfs: evidence for efficient envelope ejection in massive asymptotic giant branch stars
FR3060252A1 (fr) Systemes de production de radionucleides
Macias et al. Cosmic ray models of the ridge-like excess of gamma rays in the Galactic Centre
CN110009323A (zh) 区块链交易方法及装置、电子设备、存储介质
Fayolle et al. Decoupled and coupled moons’ ephemerides estimation strategies application to the JUICE mission
Unnikrishnan Correlation functions, Bell's inequalities and the fundamental conservation laws
Reed et al. Analysis Coordinator Report
Dubedout et al. A GRASP for real life inventory routing problem: application to bulk gas distribution
Van Leeuwen et al. Dynamic modelling of the Hipparcos attitude
Ragozzine et al. Kepler-9: The first multi-transiting system and the first transit timing variations
Boehnke et al. Xenon isotopes in the MORB source, not distinctive of early global degassing

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 3

PLSC Publication of the preliminary search report

Effective date: 20200508

PLFP Fee payment

Year of fee payment: 4

PLFP Fee payment

Year of fee payment: 5

PLFP Fee payment

Year of fee payment: 6