Brevet FR3069165A1 — publié le 20190125

SYSTEME AUTOMATISE DE REGULATION DE LA GLYCEMIE D'UN PATIENT Domaine La présente demande concerne le domaine des systèmes automatisés de régulation de glycémie, aussi appelés pancréas artificiels. Exposé de l'art antérieur Un pancréas artificiel est un système permettant de réguler automatiquement les apports en insuline d'un patient diabétique à partir de son historique de glycémie, de son historique de prise de repas, et de son historique d'injection d'insuline. On s'intéresse ici plus particulièrement aux systèmes de régulation de type MPC (de l'anglais "Model-based Predictive Control"), aussi appelés systèmes à commande prédictive, dans lesquels la régulation de la dose d'insuline administrée tient compte d'une prédiction de l'évolution future de la glycémie du patient, réalisée à partir d'un modèle physiologique décrivant l'assimilation de l'insuline par le corps du patient et son impact sur la glycémie du patient. Il serait souhaitable de pouvoir améliorer les 20 performances des pancréas artificiels à commande prédictive, et, plus particulièrement, de pouvoir améliorer la qualité de la prédiction de la glycémie future du patient, de façon à pouvoir 2 contrôler avec une plus grande pertinence les apports en insuline et limiter les risques de placer le patient en situation d'hyperglycémie ou d'hypoglycé mie. Résumé 5 Ainsi, un mode de réalisation prévoit un système automatisé de régulation de la glycémie d'un patient, comportant : un capteur de glycémie ; un dispositif d'injection d'insuline ; et une unité de traitement et de contrôle adaptée à 10 prédire, à partir d'un modèle physiologique, l'évolution future de la glycémie du patient sur une période de prédiction, et à commander le dispositif d'injection d'insuline en tenant compte de cette prédiction, dans lequel l'unité de traitement et de contrôle est 15 adaptée à : a) mettre en oeuvre une étape de calibration automatique du modèle physiologique en tenant compte de la glycémie mesurée par le capteur au cours d'une période d'observation passée ; b) à l'issue de l'étape de calibration, calculer au 20 moins un indicateur numérique représentatif de l'erreur entre la glycémie estimée à partir du modèle et la glycémie réelle mesurée par le capteur ; et c) ajuster la durée de la période de prédiction en tenant compte de la valeur de l'indicateur numérique. 25 Selon un mode de réalisation, l'indicateur numérique comprend l'écart quadratique moyen entre la glycémie estimée à partir du modèle et la glycémie réelle mesurée par le capteur pendant la période d'observation passée. Selon un mode de réalisation, l'indicateur numérique 30 comprend la différence entre la glycémie réelle mesurée par le capteur et la glycémie estimée par le modèle à un instant donné. Selon un mode de réalisation, l'indicateur numérique comprend la différence entre la dérivée de la glycémie réelle mesurée par le capteur et la dérivée de la glycémie estimée par 35 le modèle à un instant donné. 3 Selon un mode de réalisation, l'unité de traitement et de contrôle est configurée pour, à l'étape c), comparer la valeur de l'indicateur numérique à des premiers seuils, et sélectionner la durée de la période de prédiction parmi une pluralité de durées 5 prédéfinies en fonction du résultat de la comparaison. Selon un mode de réalisation, l'unité de traitement et de contrôle est en outre adaptée, après l'étape b), à déterminer, à partir de la valeur de l'indicateur numérique, si le modèle est suffisamment fiable pour pouvoir servir de base à la commande du 10 dispositif d'injection d'insuline, et, dans le cas contraire, à commander le dispositif d'injection d'insuline selon une méthode de substitution, sans tenir compte de la prédiction réalisée à partir du modèle. Selon un mode de réalisation, pour déterminer si le 15 modèle est suffisamment fiable, l'unité de traitement et de contrôle compare la valeur de l'indicateur numérique à un deuxième seuil. Selon un mode de réalisation, l'unité de traitement et de contrôle est adaptée à déterminer et ajuster automatiquement le deuxième seuil à partir de données passées mesurées sur le patient, de façon à ce que la commande du dispositif d'injection d'insuline soit basée sur les prédictions réalisées par le modèle au moins un certain pourcentage P du temps. Selon un mode de réalisation, la méthode de substitution 25 est une méthode de commande prédictive basée sur un modèle physiologique simplifié. Selon un mode de réalisation, la méthode de substitution consiste à commander le dispositif d'injection d'insuline pour délivrer des doses d'insuline préprogrammées correspondant à un 30 débit basal de référence prescrit au patient. Selon un mode de réalisation, la méthode de substitution consiste à commander le dispositif d'injection d'insuline pour administrer des doses d'insuline déterminées par l'unité de traitement et de contrôle en fonction du niveau courant de 4 glycémie mesuré par le capteur et/ou de la vitesse de variation de la glycémie mesurée par le capteur. Brève description des dessins Ces caractéristiques et avantages, ainsi que d'autres, 5 seront exposés en détail dans la description suivante de modes de réalisation particuliers faite à titre non limitatif en relation avec les figures jointes parmi lesquelles : la figure 1 représente de façon schématique, sous forme de blocs, un exemple d'un mode de réalisation d'un système 10 automatisé de régulation de la glycémie d'un patient ; la figure 2 est une représentation simplifiée d'un modèle physiologique utilisé dans le système de la figure 1 pour prédire l'évolution future de la glycémie du patient ; la figure 3 est un diagramme illustrant un exemple d'un 15 procédé automatisé de régulation de glycémie pouvant être mis en oeuvre par le système de la figure 1 ; et la figure 4 est un diagramme illustrant plus en détail un exemple d'un mode de réalisation d'un procédé automatisé de régulation de glycémie mis en oeuvre par le système de la figure 1. 20 Description détaillée De mêmes éléments ont été désignés par de mêmes références dans les différentes figures. Par souci de clarté, seuls les éléments qui sont utiles à la compréhension des modes de réalisation décrits ont été représentés et sont détaillés. En particulier, le dispositif de mesure de glycémie et le dispositif d'injection d'insuline du système de régulation décrit n'ont pas été détaillés, les modes de réalisation décrits étant compatibles avec tous ou la plupart des dispositifs de mesure de glycémie et d'injection d'insuline connus. De plus, la réalisation matérielle de l'unité de traitement et de contrôle du système de régulation décrit n'a pas été détaillée, la réalisation d'une telle unité de traitement et de contrôle étant à la portée de l'homme du métier à partir des indications fonctionnelles de la présente description. 5 La figure 1 représente de façon schématique, sous forme de blocs, un exemple d'un mode de réalisation d'un système automatisé de régulation de la glycémie d'un patient. Le système de la figure 1 comprend un capteur 101 (CG) 5 adapté à mesurer la glycémie du patient. En fonctionnement normal, le capteur 101 peut être positionné à demeure sur ou dans le corps du patient, par exemple à hauteur de son abdomen. Le capteur 101 est par exemple un capteur de type CGM (de l'anglais "Continuous Glucose Monitoring" - surveillance continue de glycémie), c'est-10 à-dire un capteur adapté à mesurer en continu (par exemple au moins une fois toutes les cinq minutes) la glycémie du patient. Le capteur 101 est par exemple un capteur de glycémie en sous-cutané. Le système de la figure 1 comprend en outre un dispositif 15 d'injection d'insuline 103 (PMP), par exemple un dispositif d'injection en sous-cutané. Le dispositif 103 est par exemple un dispositif d'injection automatique de type pompe à insuline, comportant un réservoir d'insuline relié à une aiguille d'injection implantée sous la peau du patient, la pompe pouvant 20 être commandée électriquement pour injecter automatiquement des doses d'insuline déterminées à des instants déterminés. En fonctionnement normal, le dispositif d'injection 103 peut être positionné à demeure dans ou sur le corps du patient, par exemple au niveau de son abdomen. 25 Le système de la figure 1 comprend en outre une unité de traitement et de contrôle 105 (CTRL) reliée d'une part au capteur de glycémie 101, par exemple par liaison filaire ou par liaison radio (sans fil), et d'autre part au dispositif d'injection 103, par exemple par liaison filaire ou radio. En 30 fonctionnement, l'unité de traitement et de contrôle 105 est adaptée à recevoir les données de glycémie du patient mesurées par le capteur 101, et à commander électriquement le dispositif 103 pour injecter au patient des doses d'insuline déterminées à des instants déterminés. Dans cet exemple, l'unité de traitement et de contrôle 105 est en outre adaptée à recevoir, par 6 l'intermédiaire d'une interface utilisateur non détaillée, des données cho(t) représentatives de l'évolution, en fonction du temps, de la quantité de glucose ingérée par le patient. L'interface utilisateur peut en outre être prévue pour permettre 5 de saisir des informations supplémentaires susceptibles de faciliter la régulation de la glycémie, par exemple des informations relatives à l'activité physique du patient, ou encore à son stress, ou toute autre information relative au métabolisme du patient, ou encore les types d'aliments ingérés par le patient 10 (gras ou non par exemple). L'unité de traitement et de contrôle 105 est adaptée à déterminer les doses d'insuline à injecter au patient en tenant compte notamment de l'historique de glycémie mesurée par le capteur 101, de l'historique d'insuline injectée par le dispositif 15 103, et de l'historique d'ingestion de glucose par le patient (ainsi que des éventuelles informations supplémentaires susmentionnées). Pour cela, l'unité de traitement et de contrôle 105 comprend un circuit de calcul numérique (non détaillé), comprenant par exemple un microprocesseur. L'unité de traitement 20 et de contrôle 105 est par exemple un dispositif mobile transpporté par le patient tout au long de la journée et/ou de la nuit, par exemple un dispositif de type smartphone configuré pour mettre en oeuvre un procédé de régulation du type décrit ci-après. Dans le mode de réalisation de la figure 1, l'unité de 25 traitement et de contrôle 105 est adaptée à déterminer la quantité d'insuline à administrer au patient en tenant compte d'une prédiction de l'évolution future de sa glycémie en fonction du temps. Plus particulièrement, l'unité de traitement et de contrôle 105 est adaptée, à partir de l'historique d'insuline injectée et 30 de l'historique de glucose ingéré (ainsi que des éventuelles informations supplémentaires susmentionnées), et en se basant sur un modèle physiologique décrivant l'assimilation de l'insuline par le corps du patient et son impact sur la glycémie, à déterminer une courbe représentative de l'évolution attendue de la glycémie 35 du patient en fonction du temps, sur une période à venir appelée 7 période de prédiction ou horizon de prédiction, par exemple une période de 1 à 10 heures. En tenant compte de cette courbe, l'unité de traitement et de contrôle 105 détermine les doses d'insuline qu'il conviendrait d'injecter au patient pendant la période de 5 prédiction à venir, pour que la glycémie réelle (par opposition à la glycémie estimée à partir du modèle physiologique) du patient reste dans des limites acceptables, et en particulier pour limiter les risques d'hyperglycémie ou d'hypoglycémie. Dans ce mode de fonctionnement, comme cela sera expliqué plus en détail ci-après, 10 les données de glycémie réelle mesurées par le capteur 101 sont utilisées principalement à des fins de calibration du modèle physiologique. La figure 2 est une représentation simplifiée d'un modèle physiologique MPC utilisé dans le système de la figure 1 15 pour prédire l'évolution future de la glycémie du patient. Sur la figure 2, le modèle est représenté sous la forme d'un bloc de traitement comportant : une entrée el sur laquelle est appliquée un signal i(t) représentatif de l'évolution, en fonction du temps t, de la 20 quantité d'insuline injectée au patient ; une entrée e2 sur laquelle est appliquée un signal cho(t) représentatif de l'évolution, en fonction du temps t, de la quantité de glucose ingérée par le patient ; et une sortie s fournissant un signal G(t) représentatif 25 de l'évolution, en fonction du temps t, de la glycémie du patient. Le modèle physiologique MPC est par exemple un modèle compartimental comportant, outre les variables d'entrée i(t) et cho(t) et la variable de sortie G(t), une pluralité de variables d'états correspondant à des variables physiologiques du patient, 30 évoluant en fonction du temps. L'évolution temporelle des variables d'état et de la variable de sortie G(t) est régie par un système d'équations différentielles comportant une pluralité de paramètres représentés sur la figure 2 par un vecteur [PARAM] appliqué sur une entrée pl du bloc MPC. La réponse du modèle 35 physiologique est en outre conditionnée par les états initiaux ou 8 valeurs initiales affectés aux variables d'état, représentés sur la figure 2 par un vecteur [INIT] appliqué sur une entrée p2 du bloc MPC. A titre d'exemple, le modèle physiologique MPC utilisé 5 dans le système de la figure 1 est le modèle dit de Hovorka, décrit dans l'article intitulé "Nonlinearmodel predictive control of glucose concentration in subjects with type 1 diabetes" de Roman Hovorka et al. (Physiol Meas. 2004;25:905-920), et dans l'article intitulé "Partitioning glucose distribution/transport, 10 disposal, and endogenous production during IVGTT", de Roman Hovorka et al. (Am J Physiol Endocrinol Metab 282: E992-E1007, 2002). Plus généralement, tout autre modèle physiologique décrivant l'assimilation de l'insuline par le corps d'un patient et son effet sur la glycémie du patient peut être utilisé, par 15 exemple le modèle dit de Cobelli, décrit dans l'article intitulé "A. System Model of Oral Glucose Absorption: Validation on Gold Standard Data", de Chiara Dalla Man et al. (IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 53, NO. 12, DECEMBER 2006). Parmi les paramètres du vecteur [PARAM], certains 20 peuvent être considérés comme constants pour un patient donné. D'autres paramètres, appelés ci-après paramètres temps-dépendants, sont en revanche susceptibles d'évoluer dans le temps. Du fait de cette variabilité de certains paramètres du système, il est en pratique nécessaire de ré-étalonner ou re-calibrer 25 régulièrement le modèle en cours d'utilisation, par exemple toutes les 1 à 20 minutes, par exemple toutes les 5 minutes, pour s'assurer que les prédictions du modèle restent pertinentes. Cette mise à jour du modèle, aussi appelée personnalisation du modèle, doit pouvoir être réalisée de façon automatique par le système de 30 la figure 1, c'est-à-dire sans qu'il soit nécessaire de mesurer physiquement les paramètres temps-dépendants du système sur le patient puis de les transmettre à l'unité de traitement et de contrôle 105. 9 La figure 3 est un diagramme illustrant un exemple d'un procédé automatisé de régulation de glycémie pouvant être mis en oeuvre par le système de la figure 1. Ce procédé comprend une étape 301 de re-calibration ou 5 mise à jour du modèle, pouvant par exemple être répétée à intervalles réguliers, par exemple toutes les 1 à 20 minutes. Lors de cette étape, l'unité de traitement et de contrôle 105 met en oeuvre un procédé de ré-estimation des paramètres temps-dépendants du modèle en tenant compte des données d'insuline effectivement 10 injectée par le dispositif 103 et des données de glycémie réelle mesurée par le capteur 101 pendant une période d'observation passée de durée AT, par exemple une période de 1 à 10 heures précédant l'étape de calibration. Plus particulièrement, lors de l'étape de calibration, l'unité de traitement et de contrôle 105 15 simule le comportement du patient sur la période d'observation passée à partir du modèle physiologique (en tenant compte des éventuelles ingestions de glucose et injections d'insuline pendant cette période), et compare la courbe de glycémie estimée par le modèle à la courbe de glycémie réelle mesurée par le capteur 20 pendant cette même période. L'unité de traitement et de contrôle 105 recherche alors, pour les paramètres temps-dépendants du modèle, un jeu de valeurs conduisant à minimiser une grandeur représentative de l'erreur entre la courbe de glycémie estimée par le modèle et la courbe de glycémie réelle pendant la période 25 d'observation. A titre d'exemple, l'unité de traitement et de contrôle recherche un jeu de paramètres conduisant à minimiser un indicateur m représentatif de l'aire entre la courbe de glycémie estimée par le modèle et la courbe de glycémie réelle pendant la période d'observation, aussi appelé écart quadratique moyen entre 30 la glycémie estimée et la glycémie réelle, par exemple défini comme suit : m= to Ig(t) - ,(t)12 t= to -AT 10 où t est la variable temps discrétisée, t0-AT correspond à l'instant de début de la phase d'observation passé, to correspond à l'instant de fin de la phase d'observation passée (correspondant par exemple à l'instant de début de l'étape de calibration du 5 modèle), g est la courbe d'évolution temporelle de la glycémie réelle mesurée par le capteur 101 pendant la période [t0-AT, -Lo], et § est la courbe de glycémie estimée à partir du modèle pendant la période [t0-AT, to]. A titre de variante, pour le calcul de l'écart quadratique moyen, la variable AT peut être remplacée par 10 le nombre de mesures réalisées pendant la période d'observation passée. L'algorithme de recherche de paramètres optimaux utilisé lors de cette étape n'est pas détaillé dans la présente demande, les modes de réalisation décrits étant compatibles avec les algorithmes usuels utilisés dans des domaines variés pour résoudre 15 des problèmes d'optimisation de paramètres par minimisation d'une fonction de coût. On notera que lors de l'étape 301, outre les paramètres temps-dépendants du modèle, l'unité de traitement et de contrôle 105 définit un vecteur [INIT] d'états initiaux (états à l'instant 20 t0-AT) des variables d'état du modèle, pour pouvoir simuler le comportement du patient à partir du modèle. Pour définir les états initiaux des variables d'état du modèle, une première possibilité consiste à faire l'hypothèse que, dans la période précédant la période d'observation [t0-AT, to] sur laquelle est basée la 25 calibration du modèle, le patient se trouvait dans un état stationnaire, avec un débit d'insuline injectée constant, et une prise alimentaire de glucose nulle. Sous cette hypothèse, toutes les dérivées du système d'équations différentielles peuvent être considérées comme nulles à l'instant initial t0-AT. Les valeurs à 30 l'instant t0-AT des variables d'état du système peuvent alors être calculées analytiquement. Pour améliorer l'initialisation, une autre possibilité consiste à faire les mêmes hypothèses que précédemment, mais en ajoutant la contrainte que la glycémie estimée à l'instant t0-AT soit égale à la glycémie réelle mesurée 35 par le capteur. Pour améliorer encore l'initialisation, une autre 11 possibilité est de considérer les états initiaux des variables d'état du modèle comme des variables aléatoires, au même titre que les paramètres temps-dépendants du modèle. Les états initiaux des variables d'état sont alors déterminés de la même façon que 5 les paramètres temps-dépendants du modèle, c'est-à-dire que l'unité de traitement et de contrôle 105 recherche un jeu de valeurs d'états initiaux [INIT] conduisant à minimiser une grandeur représentative de l'erreur entre la courbe de glycémie estimée par le modèle et la courbe de glycémie réelle pendant la 10 période d'observation passée. Le procédé de la figure 3 comprend en outre, après l'étape 301, une étape 303 de prédiction, par l'unité de traitement et de contrôle 105, de l'évolution temporelle de la glycémie du patient sur une période de prédiction à venir [L°, 15 teTpred] de durée Tp red , par exemple comprise entre 1 et 10 heures, à partir du modèle physiologique mis à jour à l'étape 301 et en tenant compte de l'historique d'insuline injectée au patient et de l'historique de glucose ingéré par le patient. Le procédé de la figure 3 comprend de plus, après l'étape 20 303, une étape 305 de détermination, par l'unité de traitement et de contrôle 105, en tenant compte de la courbe de glycémie future prédite à l'étape 303, des doses d'insuline à injecter au patient pendant la période de prédiction à venir [tO, to+Tp red] A l'issue de cette étape, l'unité de traitement et de contrôle 105 peut 25 programmer le dispositif d'injection 103 pour administrer les doses déterminées pendant la période de prédiction PL°, to+Tp redl- Les étapes 303 de prédiction de la glycémie et 305 et détermination des doses futures d'insuline à administrer peuvent par exemple être réitérées à chaque mise à jour du modèle 30 physiologique (c'est-à-dire après chaque itération de l'étape 301), à chaque nouvelle ingestion de glucose signalée par le patient, et/ou à chaque nouvelle administration d'une dose d'insuline par le dispositif d'injection 103. Dans le procédé susmentionné, la durée T- pred de la 35 période de prédiction de l'évolution future de la glycémie du 12 patient est un paramètre important, conditionnant la performance du système de régulation. Compte tenu de la dynamique relativement lente du système que l'on cherche à réguler, il serait souhaitable que la période de prédiction T- pred soit relativement longue, par 5 exemple de l'ordre de 4 heures ou plus, de façon à pouvoir anticiper et évaluer au mieux les besoins en insuline du patient. Toutefois, en pratique, les imperfections du modèle utilisé contraignent à limiter l'horizon de prédiction considéré. Selon un aspect d'un mode de réalisation, l'unité de 10 traitement et de contrôle 105 est adaptée, après chaque mise à jour du modèle physiologique (étape 301), à calculer un ou plusieurs indicateurs numériques représentatifs de la fiabilité du modèle mis à jour, et à ajuster la durée de prédiction Tp red en fonction de ces indicateurs. Plus particulièrement, si le 15 modèle mis à jour est jugé fiable, la durée de prédiction Tp red sera choisie relativement élevée, et, si le modèle est jugé peu fiable, la durée de prédiction T- pred sera choisie relativement faible. Par rapport à un système dans lequel la durée de prédiction Tpred 20 permet d'améliorer la qualité de la prédiction de la glycémie future du patient, et ainsi de contrôler avec une plus grande pertinence les apports en insuline. La figure 4 est un diagramme illustrant plus en détail un exemple d'un procédé automatisé de régulation de glycémie mis 25 en oeuvre par le système de la figure 1, dans lequel la durée de prédiction T- pred est ajustée en fonction d'une estimation de la fiabilité du modèle physiologique. Ce procédé comprend les mêmes étapes 301, 303 et 305 que dans l'exemple de la figure 3. Toutefois, le procédé de la figure 30 4 comprend en outre, après chaque étape 301 de mise à jour du modèle physiologique et avant la mise en oeuvre des étapes suivantes 303 de prédiction de la glycémie future du patient, et 305 de contrôle de la délivrance d'insuline à partir de la prédiction de glycémie, une étape 411 de calcul d'un ou plusieurs 35 indicateurs numériques de fiabilité du modèle mis à jour, et d'une est fixe, un avantage de ce mode de fonctionnement est qu'il 13 étape 413 d'ajustement de la durée de prédiction T- pred en fonction du ou des indicateurs de fiabilité calculés à l'étape 411. Lors de l'étape 411, l'unité de traitement et de contrôle 105 calcule un ou plusieurs indicateurs numériques 5 représentatifs de la fiabilité du modèle mis à jour à l'étape 301. A titre d'exemple, l'unité de traitement et de contrôle calcule trois indicateurs numériques de fiabilité MM, GD et SD. L'indicateur MM correspond à l'écart quadratique moyen entre la glycémie estimée à partir du modèle mis à jour et la courbe de 10 glycémie réelle mesurée par le capteur 101 pendant une période d'observation passée, par exemple une période de 1 à 10 heures précédant l'instant to, par exemple la période [t0-AT, to]. L'indicateur GD correspond à la différence entre la glycémie réelle mesurée par le capteur 101 et la glycémie estimée par le 15 modèle mis à jour à un instant donné, par exemple à l'instant to, et l'indicateur SD correspond à la différence entre la pente ou dérivée de la glycémie réelle mesurée par le capteur 101 et la pente ou dérivée de la glycémie estimée par le modèle mis à jour à un instant donné, par exemple à l'instant tO. 20 Lors de l'étape 413, l'unité de traitement et de contrôle détermine, à partir du ou des indicateurs numériques de fiabilité calculés à l'étape 411, la durée de prédiction Tp red à utiliser pour la mise en oeuvre de l'étape 303. A titre d'exemple, la durée de prédiction Tp red est choisie parmi n valeurs 25 prédéfinies Dl, Dn décroissantes, avec n entier supérieur ou égal à 2, en fonction de la valeur du ou des indicateurs numériques de fiabilité calculés à l'étape 411. A titre d'exemple, pour chacun des m indicateurs de fiabilité Ij calculé à l'étape 411, avec j entier allant de 1 à m et m entier supérieur ou égal à 1, 30 la valeur de l'indicateur est comparée à un jeu de n seuils prédéfinis S'il, ..., S'in de valeurs croissantes. L'unité de traitement et de contrôle 105 recherche alors le plus petit indice de seuil k tel que, pour chacun des m indicateurs Ij calculés à l'étape 411, la valeur de l'indicateur Ij soit inférieure au seuil 35 SIjk. L'horizon de prédiction T- pred est alors choisi égal à la 14 durée Dk. Plus généralement, selon l'objectif recherché, d'autres fonctions et/ou règles de décision permettant de déterminer la durée de prédiction T- pred à partir du ou des indicateurs de fiabilité calculés à l'étape 411 peuvent être implémentées. 5 Après l'étape 413, les étapes 303 et 305 peuvent être mises en oeuvre de façon similaire à ce qui a été décrit précédemment. On notera que dans certains cas, la fiabilité du modèle physiologique mis à jour à l'étape 301 peut être si faible qu'il 10 est préférable de cesser d'utiliser le modèle pour réguler la glycémie du patient. Dans l'exemple de la figure 4, l'unité de contrôle et de traitement 105 du système de régulation est en outre adaptée, après chaque mise à jour ou re-calibration du modèle physiologique 15 (étape 301), à partir du ou des indicateurs de fiabilité calculés à l'étape 411, à déterminer si le modèle mis à jour est suffisamment fiable pour pouvoir être utilisé pour réguler la glycémie du patient. Plus particulièrement, le procédé de la figure 4 20 comprend, entre les étapes 411 et 413, une étape 451 de vérification de la fiabilité du modèle mis à jour à l'étape 301. A titre d'exemple, la fiabilité du modèle peut être considérée comme suffisante par l'unité de traitement et de contrôle 105 lorsque les valeurs des indicateurs calculés à l'étape 411 sont 25 inférieures à des seuils prédéfinis, et insuffisante dans le cas contraire. A titre d'exemple, en reprenant les notations définies ci-dessus, la fiabilité du modèle peut être considérée comme suffisante par l'unité de traitement et de contrôle 105 lorsque pour chacun des m indicateurs de fiabilité Ij calculés à l'étape 30 411, la valeur de l'indicateur est inférieure au seuil SIjn correspondant, et insuffisante lorsque pour au moins l'un des indicateurs Ij, la valeur de l'indicateur est supérieure au seuil SIjn correspondant. Plus généralement, tout autre critère de qualité ou toute autre combinaison de critères de qualité peuvent 15 être utilisés à l'étape 451 pour déterminer si le modèle physiologique re-calibré à l'étape 301 est suffisamment fiable. Si le modèle physiologique est considéré comme suffisamment fiable à l'étape 451 (0), les étapes 413, 303 et 305 5 peuvent être mises en oeuvre de façon similaire à ce qui a été décrit précédemment, c'est-à-dire que l'unité de traitement et de contrôle 105 continue de se baser sur les prédictions réalisées par le modèle physiologique pour réguler l'administration de l'insuline au patient, en ajustant l'horizon de prédiction Tp red 10 en fonction du degré de fiabilité du modèle. Si le modèle physiologique est jugé insuffisamment fiable à l'étape 451 (N), l'unité de traitement et de contrôle 105 cesse d'utiliser ce modèle pour réguler l'administration de l'insuline au patient, et met en oeuvre une méthode de régulation 15 de substitution lors d'une étape 453. A titre d'exemple, lors de l'étape 453, l'unité de traitement et de contrôle 105 utilise un modèle physiologique simplifié, par exemple un modèle compartimentai comportant un nombre de variables d'état et un nombre de paramètres réduit par 20 rapport au modèle initial, pour prédire l'évolution de la glycémie du patient et réguler en conséquence l'injection d'insuline. A titre de variante, lors de l'étape 453, l'unité de traitement et de contrôle 105 cesse de mettre en oeuvre une commande prédictive, c'est-à-dire qu'elle cesse d'utiliser un 25 modèle physiologique pour prédire la glycémie future du patient et réguler en conséquence l'injection d'insuline. Dans ce cas, l'unité de traitement et de contrôle 105 commande par exemple le dispositif d'injection d'insuline 103 pour administrer des doses préprogrammées d'insuline, correspondant par exemple à un débit 30 basal de référence prescrit au patient. Alternativement, l'unité de traitement et de contrôle 105 utilise un algorithme de type matrice décisionnelle pour déterminer les doses d'insuline à administrer au patient, en fonction de divers paramètres observés tels que le niveau courant de glycémie mesuré par le capteur 101, 16 ou encore la vitesse de variation (ou pente) de la glycémie sur une période passée. Une telle méthode de substitution peut par exemple être utilisée pendant une période de temps prédéterminée. A l'issue de 5 cette période, les étapes 301 de calibration du modèle physiologique principal, 411 de calcul du ou des indicateurs de fiabilité du modèle physiologique principal, et 451 d'estimation de la qualité du modèle physiologique principal, peuvent être réitérées, pour, si la qualité du modèle physiologique principal 10 est jugée suffisante, réactiver l'utilisation du modèle principal pour réguler l'administration de l'insuline au patient. A titre d'exemple, les seuils utilisés à l'étape 451 pour déterminer si le modèle physiologique principal est suffisamment fiable pour pouvoir être utilisé sont choisis de 15 façon à maximiser la probabilité que le système de régulation fonctionne au moins un certain pourcentage P du temps, par exemple au moins 70% du temps, sur la base du modèle physiologique principal. Les seuils utilisés aux étapes 451 et 413 sont par 20 exemples déterminés sur la base d'une analyse de données passées mesurée sur un échantillon de plusieurs patients. A titre d'exemple, on peut rejouer l'algorithme de régulation pour une pluralité de patients sur un banc de test, et, pour chaque patient, à chaque mise à jour du modèle physiologique, pour chacune des n 25 valeurs possibles Dl, Dn de la durée de prédiction Tp red, calculer l'écart quadratique moyen, sur la période de prédiction Tpred, entre la glycémie estimée à partir du modèle mis à jour et la courbe de glycémie réelle mesurée par le capteur 101. A chaque mise à jour du modèle, on calcule en outre les m indicateurs de 30 fiabilité du modèle mis à jour Ii, Im. Ainsi, à chaque mise à jour du modèle, on dispose d'un jeu de m valeurs correspondant aux indicateurs de fiabilité du modèle tels que défini ci-dessus, et d'un jeu de n valeurs correspondant à des mesures effectives de fiabilité du modèle pour les n durées de prédiction Dl, ..., 35 Dn considérées. L'étude des corrélations entre les indicateurs de 17 fiabilité du modèle et les mesures effectives de fiabilité permet de déterminer les seuils à utiliser à l'étape 413 pour choisir la durée de la période de prédiction après chaque mise à jour du modèle, et/ou à l'étape 451 pour décider s'il convient ou non de 5 basculer sur une méthode de régulation de substitution. La détermination des seuils à partir des valeurs susmentionnées d'indicateurs de fiabilité et de mesures effectives de fiabilité peut être totalement ou partiellement automatisée. A titre de variante, les seuils utilisés aux étapes 451 10 et 413 sont déterminés de façon similaire à ce qui vient d'être décrit, mais uniquement sur la base des données passées mesurées sur le patient utilisateur du système, ce qui permet de personnaliser le fonctionnement du système de régulation. Dans ce cas, l'unité de traitement et de contrôle 105 peut être configurée 15 pour recalculer régulièrement les seuils utilisés à l'étape 413 et/ou à l'étape 451, REVENDICATIONS1. Système automatisé de régulation de la glycémie d'un patient, comportant : un capteur de glycémie (101) ; un dispositif d'injection d'insuline (103) ; et une unité de traitement et de contrôle (105) adaptée à prédire, à partir d'un modèle physiologique, l'évolution future de la glycémie du patient sur une période de prédiction, et à commander le dispositif d'injection d'insuline (103) en tenant compte de cette prédiction, dans lequel l'unité de traitement et de contrôle (105) est adaptée à : a) mettre en oeuvre une étape de calibration automatique du modèle physiologique en tenant compte de la glycémie mesurée par le capteur (101) au cours d'une période d'observation passée ; b) à l'issue de l'étape de calibration, calculer au moins un indicateur numérique représentatif de l'erreur entre la glycémie estimée à partir du modèle et la glycémie réelle mesurée par le capteur ; et c) ajuster la durée de la période de prédiction en tenant 20 compte de la valeur dudit au moins un indicateur numérique. 2. Système selon la revendication 1, dans lequel ledit au moins un indicateur numérique comprend l'écart quadratique moyen entre la glycémie estimée à partir du modèle et la glycémie réelle mesurée par le capteur (101) pendant la période 25 d'observation passée. 3. Système selon la revendication 1 ou 2, dans lequel ledit au moins un indicateur numérique comprend la différence entre la glycémie réelle mesurée par le capteur (101) et la glycémie estimée par le modèle à un instant donné. 30 4. Système selon l'une quelconque des revendications 1 à 3 dans lequel ledit au moins un indicateur numérique comprend la différence entre la dérivée de la glycémie réelle mesurée par le capteur (101) et la dérivée de la glycémie estimée par le modèle à un instant donné. 19 5. Système selon l'une quelconque des revendications 1 à 4, dans lequel l'unité de traitement et de contrôle (105) est configurée pour, à l'étape c), comparer la valeur dudit au moins un indicateur numérique à des premiers seuils, et sélectionner la 5 durée de la période de prédiction parmi une pluralité de durées prédéfinies en fonction du résultat de la comparaison. 6. Système selon l'une quelconque des revendications 1 à 5, dans lequel l'unité de traitement et de contrôle (105) est en outre adaptée, après l'étape b), à déterminer, à partir de la 10 valeur dudit au moins un indicateur numérique, si le modèle est suffisamment fiable pour pouvoir servir de base à la commande du dispositif d'injection d'insuline, et, dans le cas contraire, à commander le dispositif d'injection d'insuline (103) selon une méthode de substitution, sans tenir compte de la prédiction 15 réalisée à partir du modèle. 7. Système selon la revendication 6, dans lequel l'unité de traitement et de contrôle (105) est configurée pour, pour déterminer si le modèle est suffisamment fiable, comparer la valeur dudit au moins un indicateur numérique à un deuxième seuil. 8. Système selon la revendication 7, dans lequel l'unité de traitement et de contrôle (105) est adaptée à déterminer et ajuster automatiquement le deuxième seuil à partir de données passées mesurées sur le patient, de façon à ce que la commande du dispositif d'injection d'insuline (103) soit basée sur les prédictions réalisées par le modèle au moins un certain pourcentage du temps. 9. Système selon l'une quelconque des revendications 6 à 8, dans lequel la méthode de substitution est une méthode de commande prédictive basée sur un modèle physiologique simplifié. 10. Système selon l'une quelconque des revendications 6 à 8, dans lequel la méthode de substitution consiste à commander le dispositif d'injection d'insuline (103) pour délivrer des doses d'insuline préprogrammées correspondant à un débit basal de référence prescrit au patient. 20 11. Système selon l'une quelconque des revendications 6 à 8, dans lequel la méthode de substitution consiste à commander le dispositif d'injection d'insuline (103) pour administrer des doses d'insuline déterminées par l'unité de traitement et de 5 contrôle (105) en fonction du niveau courant de glycémie mesuré par le capteur (101) et/ou de la vitesse de variation de la glycémie mesurée par le capteur (101).