Brevet FR3079130A1 — publié le 20190927

B16770 - DD18479 1 SYSTEME DE PREDICTION DE LA GLYCENUE D'UN PATIENT Domaine La présente demande concerne le domaine des systèmes de contrôle de glycémie de façon générale, et vise plus particulièrement un système de prédiction de l'évolution future 5 de la glycémie d'un patient. Exposé de l'art antérieur On a déjà proposé, par exemple dans la demande de brevet français N°1658881 (B15018/DD16959) déposée le 21 septembre 2016, dans la demande de brevet français N°1658882 (B15267/ DD17175) 10 déposée le 21 septembre 2016, et dans la demande de brevet français N°1756960 (B15860/DD18479) déposée le 21 juillet 2017, des systèmes automatisés de régulation de glycémie, aussi appelés pancréas artificiels, permettant de réguler automatiquement les apports en insuline d'un patient diabétique à partir de son 15 historique de glycémie, de son historique de prise de repas, et de son historique d'injection d'insuline. Les systèmes de régulation décrits dans les demandes de brevet susmentionnées sont des systèmes de type MPC (de l'anglais "Model-based Predictive Control"), aussi appelés systèmes à 20 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 B16770 - DD18479 2 d'un modèle. Plus particulièrement, dans les systèmes décrits dans les demandes de brevet susmentionnées, le modèle utilisé pour prédire l'évolution future de la glycémie du patient est 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. De façon plus générale, de nombreuses applications automatisées de contrôle de glycémie utilisent des modèles physiologiques pour prédire l'évolution future de la glycémie d'un patient, et mettre en oeuvre, en tenant compte de cette prédiction, des actions visant à maintenir la glycémie du patient dans une plage souhaitée. Les modèles physiologiques connus pour la prédiction de l'évolution future de la glycémie d'un patient présentent toutefois certaines limitations. En particulier, dans certaines conditions, il arrive que les prédictions réalisées par les modèles physiologiques connus ne soient pas fiables. Ceci peut conduire à des erreurs lors du contrôle ou lors de régulation de la glycémie du patient, ce qui peut entrainer des risques pour le patient. Résumé Ainsi, un mode de réalisation prévoit un système automatisé de contrôle de la glycémie d'un patient, comportant une unité de traitement et de contrôle adaptée à prédire, à partir d'un premier modèle mathématique, l'évolution future de la glycémie G(t) du patient sur une période de prédiction, dans lequel le premier modèle prend pour entrées : une variable EE(t) représentative de l'évolution, en fonction du temps, de la dépense énergétique du patient ; une variable I0B(t) représentative de l'évolution, en 30 fonction du temps, de la quantité d'insuline embarquée du patient ; et une variable COB(t) représentative de l'évolution, en fonction du temps, de la quantité de glucose embarquée du patient, et dans lequel le premier modèle mathématique est un 35 modèle de type boîte noire, c'est-à-dire un modèle non B16770 - DD18479 3 physiologique défini par la seule observation des effets des variables d'entrée EE(t), I0B(t) et COB(t) sur la variable de sortie G(t), sans tenir compte des différents mécanismes physiologiques connus s'opérant dans le corps du patient. 5 Selon un mode de réalisation, le premier modèle mathématique est défini comme suit : G (t) = a. yG b. UEE + C. Ltiog + d. licog , où t est une variable temps discrétisée, a = [al, ..., ana] est un vecteur de paramètres de dimension na, avec na entier 10 supérieur ou égal à 1, b = [bl, bnb] est un vecteur de paramètres de dimension nb, avec nb entier supérieur ou égal à 1, c = [cl, ..., cnc] est un vecteur de paramètres de dimension nc, avec nc entier supérieur ou égal à 1, d = [dl, dnd] est un vecteur de paramètres de dimension nd, avec nd entier supérieur 15 ou égal à 1, yG = [G(t-1), G(t-na)] est un vecteur de régression de dimension na, uEE = [EE (t-nkl), EE (t-nkl-nb)] est un vecteur de régression de dimension nb, avec nid_ entier supérieur ou égal à 1, uIOB = ['OB(t-nk2), IOB (t-nk2-nc)] est un vecteur de régression de dimension nc, avec nk2 entier 20 supérieur ou égal à 1, et ucoB = [COB(t-nk3), COB(t-nk3-nd)] est un vecteur de régression de dimension nd, avec nk3 entier supérieur ou égal à 1. Selon un mode de réalisation, na est compris dans la plage [1,15], nb est compris dans la plage [1,15], nc est compris 25 dans la plage [1,15], nd est compris dans la plage [1,15], nid_ est compris dans la plage [1,15], nk2 est compris dans la plage [1,15], et nk3 est compris dans la plage [1,15]. Selon un mode de réalisation, en fonctionnement, l'unité de traitement et de contrôle reçoit un signal CHO(t) représentatif de l'évolution, en fonction du temps, de la quantité de glucose ingérée par le patient, et un signal I(t) représentatif de l'évolution, en fonction du temps, de la quantité d'insuline injectée au patient. B16770 - DD18479 4 Selon un mode de réalisation, l'unité de traitement et de contrôle est configurée pour calculer les variables d'entrée I0B(t) et COB(t) du premier modèle comme suit : I 0 B (t) Er° I(t- k) * (k) , et 5 COB(t) = Er, CHO(t - k)h (k) * -COB \--J où K est un entier supérieur à 1, k est un entier allant de 0 à K, hioB est une fonction d'action représentant l'évolution, en fonction du temps, de l'effet de l'insuline injectée sur l'absorption du glucose sanguin, et hcoB est une fonction d'action 10 représentant l'évolution, en fonction du temps, de l'effet du glucose ingéré sur la glycémie. Selon un mode de réalisation, les fonctions d'actions hIOB et hcoB sont définies comme suit : itioB (t) = [1. TIOB e rioB , et 15 hcoB(t) =[171.--ge 'COR, TCOB où TIOB et "cCOB sont des constantes de temps. Selon un mode de réalisation, le système comprend en outre un dispositif de mesure d'une activité physique du patient, la variable d'entrée EE(t) du premier modèle étant calculée par 20 l'unité de traitement et de contrôle à partir de données de sortie PA(t) du comprend 25 en outre un capteur du rythme cardiaque du patient. comprend Selon un mode de réalisation, le capteur de mouvements fournit un signal Scpm représentatif de mouvements effectués par le patient, et le capteur de rythme cardiaque fournit un signal SHR dispositif de mesure. Selon un mode de réalisation, le dispositif de mesure un capteur de mouvements du patient. Selon un mode de réalisation, le dispositif de mesure représentatif du rythme cardiaque du patient, et dans lequel B16770 - DD18479 5 la variable d'entrée EE(t) du premier modèle est calculée comme suit par l'unité de traitement et de contrôle : lai * SHR (t) + lel, SHR S1 EE(t) = a2 * SLc(t) + P2, SHR S2 où le signal SLc est une combinaison linéaire des 5 signaux SHR et Scmp, et les grandeurs a1,oc2, oc3, Pl, 132, 133, Si et S2 sont des paramètres du système. Selon un mode de réalisation, dans lequel l'unité de traitement et de contrôle est configurée pour prédire, à partir du premier modèle mathématique, l'évolution future de la glycémie 10 G(t) du patient sur une période de prédiction, déterminer si la glycémie G(t) prédite par le premier modèle reste ou non contenue dans une plage souhaitée, et déclencher une alarme si la glycémie G(t) prédite par le premier modèle sort de la plage souhaitée. Selon un mode de réalisation, le système comporte en 15 outre un capteur de glycémie, et un dispositif d'injection d'insuline, l'unité de traitement et de contrôle est adaptée à prédire l'évolution future de la glycémie du patient à partir d'un deuxième modèle mathématique et à commander le dispositif 20 d'injection d'insuline en tenant compte de cette prédiction, et l'unité de traitement et de contrôle étant configurée pour : a) mettre en oeuvre une étape de calibration automatique du deuxième modèle en tenant compte d'un historique de glycémie 25 mesurée par le capteur au cours d'une période d'observation passée ; b) à l'issue de l'étape de calibration, réaliser une première prédiction de l'évolution future de la glycémie du patient à partir du premier modèle et, sur la même période de 30 prédiction, une deuxième prédiction de l'évolution future de la glycémie du patient à partir du deuxième modèle ; B16770 - DD18479 6 c) calculer au moins un indicateur numérique représentatif de l'écart entre les première et deuxième prédictions ; et d) en fonction de la valeur de l'indicateur numérique, 5 commander le dispositif d'injection d'insuline en tenant compte de la deuxième prédiction ou commander le dispositif d'injection d'insuline sans tenir compte de la deuxième prédiction. Selon un mode de réalisation, le deuxième modèle est un modèle physiologique comportant un système d'équations 10 différentielles décrivant l'évolution d'une pluralité de variables d'état en fonction du temps. Selon un mode de réalisation, l'étape a) de calibration automatique du deuxième modèle comprend une étape d'estimation de paramètres du système d'équations différentielles par 15 minimisation d'une grandeur représentative de l'erreur, pendant une période d'observation passée, entre la glycémie estimée à partir du deuxième modèle physiologique et la glycémie mesurée par le capteur. Brève description des dessins 20 Ces caractéristiques et leurs avantages, ainsi que d'autres, 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 25 de blocs, un exemple d'un système 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 ; 30 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 ; la figure 4 représente de façon schématique, sous forme de blocs, un exemple d'un mode de réalisation d'un système 35 automatisé de prédiction de la glycémie d'un patient ; B16770 - DD18479 7 la figure 5 est une représentation simplifiée d'un modèle utilisé dans le système de la figure 4 pour prédire l'évolution future de la glycémie du patient ; la figure 6 est un diagramme illustrant un exemple d'un 5 procédé de construction du modèle de la figure 5 ; la figure 7 est un diagramme illustrant un exemple d'un procédé automatisé de contrôle de glycémie pouvant être mis en oeuvre par le système de la figure 4 ; et la figure 8 est un diagramme illustrant un autre exemple 10 d'un procédé automatisé de contrôle de glycémie pouvant être mis en oeuvre par le système de la figure 4. 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é, 15 seuls les éléments utiles à la compréhension des modes de réalisation décrits ont été représentés et sont détaillés. En particulier, la réalisation matérielle de l'unité de traitement et de contrôle des systèmes décrit n'a pas été détaillée, la réalisation d'une telle unité de traitement et de contrôle étant 20 à la portée de l'homme du métier à partir des indications fonctionnelles de la présente description. Sauf précision contraire, les expressions "approximativement", "sensiblement", et "de l'ordre de" signifient à 10 % près, de préférence à 5 % près. 25 La figure 1 représente de façon schématique, sous forme de blocs, un exemple d'un système automatisé de régulation de la glycémie d'un patient du type décrit dans les demandes de brevet français N°1658881, N°1658882 et N°1756960 susmentionnées. Le système de la figure 1 comprend un capteur 101 (CG) 30 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- à-dire un capteur adapté à mesurer en continu, par exemple de B16770 - DD18479 8 façon périodique (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 5 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 10 ê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. 15 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 20 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 (par exemple bolus et basal). Dans cet 25 exemple, l'unité de traitement et de contrôle 105 est en outre adaptée à recevoir, par 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. 30 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 103, et de l'historique d'ingestion de glucose par le patient. 35 Pour cela, l'unité de traitement et de contrôle 105 comprend un B16770 - DD18479 9 circuit de calcul numérique (non détaillé), comprenant par exemple un microprocesseur. L'unité de traitement et de contrôle 105 est par exemple un dispositif mobile porté par le patient tout au long de la journée et/ou de la nuit, par exemple un dispositif de type 5 smartphone configuré pour mettre en oeuvre un procédé de régulation du type décrit ci-après. Dans l'exemple de la figure 1, l'unité de 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 10 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 de l'historique de glucose ingéré, et en se basant sur un modèle physiologique décrivant l'assimilation de l'insuline par le corps 15 du patient et son impact sur la glycémie, à déterminer une courbe représentative de l'évolution attendue de la glycémie du patient en fonction du temps, sur une période à venir appelée 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 20 traitement et de contrôle 105 détermine les doses d'insuline qu'il conviendrait d'injecter au patient pendant la période de 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 25 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, 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. 30 La figure 2 est une représentation simplifiée d'un modèle physiologique 201 (PM) utilisé dans le système de la figure 1 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 : B16770 - DD18479 10 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 quantité d'insuline injectée au patient ; une entrée e2 sur laquelle est appliquée un signal 5 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 de l'évolution, en fonction du temps t, de la glycémie du patient. Le modèle physiologique 201 est par exemple un modèle 10 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'état correspondant à des variables physiologiques du patient, évoluant en fonction du temps. L'évolution temporelle des variables d'état et de la variable de sortie G(t) est régie par 15 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 201. La réponse du modèle physiologique est en outre conditionnée par les états initiaux ou valeurs initiales affectés aux variables d'état, représentés sur 20 la figure 2 par un vecteur [INIT] appliqué sur une entrée p2 du bloc 201. A titre d'exemple, le modèle physiologique 201 utilisé dans le système de la figure 1 est le modèle dit de Hovorka, décrit dans l'article intitulé "Nonlinearmodel predictive control 25 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, disposal, and endogenous production during IVGTT", de Roman Hovorka et al. (Am J Physiol Endocrinol Metab 282: E992-E1007, 30 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é. Parmi les paramètres du vecteur [PARAM], certains peuvent être considérés comme constants pour un patient donné. 35 D'autres paramètres, appelés ci-après paramètres temps- B16770 - DD18479 11 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 régulièrement le modèle en cours d'utilisation, par exemple toutes 5 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 la figure 1, c'est-à-dire sans qu'il soit nécessaire de mesurer 10 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. La figure 3 est un diagramme illustrant un exemple d'un procédé automatisé de régulation de glycémie pouvant être mis en 15 oeuvre par le système de la figure 1. Ce procédé comprend une étape 301 de re-calibration ou 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 20 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 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 25 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 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 30 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 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 35 représentative de l'erreur entre la courbe de glycémie estimée B16770 - DD18479 12 par le modèle et la courbe de glycémie réelle pendant la période 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 ou écart entre la courbe de 5 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 la glycémie estimée et la glycémie réelle, par exemple défini comme suit : 1 ,,to (eql) m = - AT Lit-4 -AT I Gr (t) - G (t) 1 2 10 où t est une 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 modèle), Gr est la courbe d'évolution temporelle de la glycémie 15 réelle mesurée par le capteur 101 pendant la période [t0-AT, -Lo], et G 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 le nombre de mesures réalisées pendant la période d'observation 20 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 des problèmes d'optimisation de paramètres par minimisation d'une 25 fonction de coût (avec ou sans contraintes). 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 t0-AT) des variables d'état du modèle, pour pouvoir simuler le 30 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 B16770 - DD18479 13 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 5 considérées comme nulles à l'instant initial t0-AT. Les valeurs à 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 10 estimée à l'instant t0-AT soit égale à la glycémie réelle mesurée par le capteur. Pour améliorer encore l'initialisation, une autre 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 15 des variables d'état sont alors déterminés de la même façon que 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 20 estimée par le modèle et la courbe de glycémie réelle pendant la 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 25 glycémie du patient sur une période de prédiction à venir [L°, 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. 30 Le procédé de la figure 3 comprend de plus, après l'étape 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+T- pred] . A l'issue 35 de cette étape, l'unité de traitement et de contrôle 105 peut B16770 - DD18479 14 programmer le dispositif d'injection 103 pour administrer les doses déterminées pendant la période de prédiction [tO, to+Tp red]Les étapes 303 de prédiction de la glycémie et 305 et détermination des doses futures d'insuline à administrer peuvent 5 par exemple être réitérées à chaque mise à jour du modèle 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. 10 Dans un système de régulation du type décrit en relation avec les figures 1 à 3, la fiabilité de la prédiction de l'évolution future de la glycémie est particulièrement importante pour déterminer correctement les doses d'insuline à administrer au patient, et par conséquent pour réguler correctement la 15 glycémie du patient. Plus généralement, dans de nombreuses applications de contrôle automatisé de glycémie, la prédiction de l'évolution future de la glycémie du patient joue un rôle important pour permettre le maintien de la glycémie du patient dans une plage 20 souhaitée (correspondant par exemple à une plage de normoglycémie). Comme indiqué précédemment, les modèles physiologiques connus pour la prédiction de l'évolution future de la glycémie d'un patient présentent toutefois certaines limitations. En 25 particulier, dans certaines conditions, il arrive que les prédictions réalisées par ces modèles ne soient pas fiables. Ceci peut conduire à des erreurs lors du contrôle ou lors de régulation de la glycémie du patient, ce qui peut entrainer des risques pour le patient. 30 Ainsi, il serait souhaitable de pouvoir disposer d'un système automatisé de prédiction de glycémie palliant tout ou partie des inconvénients des systèmes de prédiction de glycémie connus. Selon un aspect des modes de réalisation décrits, on 35 prévoit un système dans lequel la prédiction de l'évolution future B16770 - DD18479 15 de la glycémie du patient est réalisée en tenant compte non seulement de l'historique d'insuline injectée au patient et de l'historique de glucose ingéré par le patient, mais aussi d'une donnée représentative de l'activité physique du patient. 5 La figure 4 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 prédiction de la glycémie d'un patient. Le système de la figure 4 comprend une unité de traitement et de contrôle 401 (CTRL). En fonctionnement, l'unité 10 de traitement et de contrôle 401 est adaptée à recevoir 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, des données I(t) représentatives de l'évolution, en fonction du temps, de la quantité d'insuline injectée au patient, et des données PA(t) 15 représentatives de l'évolution, en fonction du temps, de l'activité physique du patient. Les données de glucose ingéré CHO(t) sont par exemple saisies par le patient par l'intermédiaire d'une interface utilisateur non détaillée. 20 Les données d'insuline injectée I(t) peuvent également être saisies par le patient par l'intermédiaire d'une interface utilisateur. A titre de variante, les données d'insuline injectées I(t) sont directement communiquées à l'unité de traitement et de contrôle 401 par un dispositif d'injection d'insuline. Dans 25 l'exemple de la figure 4, le système comprend un dispositif d'injection d'insuline 103 (PMP), par exemple identique ou similaire au dispositif 103 de la figure 1. Le dispositif d'injection d'insuline 103 est relié à l'unité de traitement et de contrôle 401, par exemple par liaison filaire ou par liaison 30 radio (sans fil), et communique à l'unité de traitement et de contrôle 401 les données d'insuline injectée I(t). A titre de variante, l'unité de traitement et de contrôle 105 est configurée pour commander électriquement le dispositif 103 pour injecter au patient des doses d'insuline déterminées à des instants 35 déterminés. Dans ce cas, les données d'insuline injectée I(t) sont B16770 - DD18479 16 connues de l'unité de traitement et de contrôle 401 et n'ont donc pas besoin de lui être transmises. Dans l'exemple de la figure 4, le système comprend en outre un dispositif 403 adapté à mesurer l'activité physique du 5 patient. Le dispositif 403 est reliée à l'unité de traitement et de contrôle 401, par exemple par liaison filaire ou par liaison radio (sans fil), et communique à l'unité de traitement et de contrôle 401 les données d'activité physique PA(t) du patient. Le dispositif 403 est par exemple un dispositif mobile porté par le 10 patient tout au long de la journée et/ou de la nuit. A titre d'exemple, le dispositif 403 comprend au moins un capteur de mouvements 405 (14VT), par exemple un accéléromètre. Le dispositif 403 peut en outre comprendre un capteur 407 (HR) du rythme cardiaque du patient. 15 Le système de la figure 4 peut de plus comporter un capteur 101 (CG), par exemple identique ou similaire au capteur 101 du système de la figure 1, adapté à mesurer la glycémie du patient. Le capteur 101 peut être relié, par exemple par liaison filaire ou par liaison radio (sans fil), à l'unité de traitement 20 et de contrôle 401. A titre d'exemple, en fonctionnement, l'unité de traitement et de contrôle 401 est adaptée à recevoir les données de glycémie du patient mesurées par le capteur 101. Dans le mode de réalisation de la figure 4, l'unité de traitement et de contrôle 401 est adaptée à prédire l'évolution 25 future de la glycémie du patient, et, éventuellement, à mettre en oeuvre une action de contrôle ou de régulation de la glycémie du patient, dans le but de maintenir la glycémie du patient dans une plage souhaitée. Plus particulièrement, l'unité de traitement et de contrôle 401 est adaptée, à partir de l'historique d'insuline 30 injectée, de l'historique de glucose ingéré, et de l'historique d'activité physique du patient, et en se basant sur un modèle mathématique, à déterminer une courbe représentative de l'évolution attendue de la glycémie du patient en fonction du temps, sur une période à venir. Pour cela, l'unité de traitement 35 et de contrôle 401 comprend un circuit de calcul numérique (non B16770 - DD18479 17 détaillé), comprenant par exemple un microprocesseur. L'unité de traitement et de contrôle 401 est par exemple un dispositif mobile porté 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 5 en oeuvre un procédé de prédiction et/ou de contrôle de la glycémie du patient. La figure 5 est une représentation simplifiée d'un modèle mathématique 501 (BLACK BOX) utilisé dans le système de la figure 4 pour prédire l'évolution future de la glycémie du 10 patient. Sur la figure 5, le modèle est représenté sous la forme d'un bloc de traitement comportant : une entrée el sur laquelle est appliqué un signal EE(t) représentatif de l'évolution, en fonction du temps t, de la dépense énergétique du patient, c'est-à-dire de l'effort physique 15 réalisé par le patient à l'instant t (à titre d'exemple, le signal EE(t) est représentatif de la quantité d'énergie dépensée à l'instant t par le patient pour permettre à son organisme de fonctionner correctement, d'assurer toutes ses fonctions physiques, et, le cas échéant, de réaliser un effort physique) ; 20 une entrée e2 sur laquelle est appliqué un signal I0B(t) représentatif de l'évolution, en fonction du temps t, de la quantité d'insuline embarquée ("insulin on board" en anglais) du patient, c'est-à-dire de la quantité d'insuline encore active (c'est-à-dire encore susceptible d'avoir un effet sur la glycémie) 25 à l'instant t dans le corps du patient ; une entrée e3 sur laquelle est appliqué un signal COB(t) représentatif de l'évolution, en fonction du temps t, de la quantité de glucose embarquée ("carbohydrate on board" en anglais) du patient, c'est-à-dire de la portion de la quantité de glucose 30 ingérée par le patient encore susceptible, à l'instant t, d'avoir un effet sur la glycémie du patient ; et une sortie s fournissant un signal G(t) représentatif de l'évolution, en fonction du temps t, de la glycémie du patient. Selon un aspect d'un mode de réalisation, le modèle 501 35 de la figure 5 est un modèle de type boîte noire, c'est-à-dire un B16770 - DD18479 18 modèle non physiologique, défini par la seule observation des effets des variables d'entrée EE(t), I0B(t) et COB(t) sur la sortie G(t), sans tenir compte des différents mécanismes physiologiques connus s'opérant dans le corps du patient et ayant 5 un effet sur la glycémie du patient. Le modèle 501 de la figure 5 est par exemple un modèle de type ABX, c'est-à-dire un modèle autorégressif à entrées externes. A titre d'exemple, le modèle 501 est défini par une équation unique, comme suit : G(t) = a. yG b. UEE C.Ltiog + d. licog (eq2) 10 où t est une variable temps discrétisée, a = [al, ..., ana] est un vecteur de paramètres de dimension na, avec na entier supérieur ou égal à 1 et de préférence supérieur ou égal à 2, b = [bl, bnb] est un vecteur de paramètres de dimension nb, avec nb entier supérieur ou égal à 1 et de préférence supérieur ou égal 15 à 2, c = [cl, cnc] est un vecteur de paramètres de dimension nc, avec nc entier supérieur ou égal à 1 et de préférence supérieur ou égal à 2, d = [dl, dnd] est un vecteur de paramètres de dimension nd, avec nd entier supérieur ou égal à 1 et de préférence supérieur ou égal à 2, yG = [G(t-1), G(t-na)] est un vecteur 20 de régression de dimension na, uEE = [EE (t-nkl-1), EE (t-nkl- nb)] est un vecteur de régression de dimension nb, avec nid_ entier supérieur ou égal à 1, ui0B = ['OB (t-nk2-1), IOB (t-nk2-nc)] est un vecteur de régression de dimension nc, avec nk2 entier supérieur ou égal à 1, et ucoB = [COB(t-nk3-1), COB(t-nk3- 25 nd)] est un vecteur de régression de dimension nd, avec nk3 entier supérieur ou égal à 1. Si l'on désigne par T la période d'échantillonnage du système, pour chacune des variables G, EE, IOB et COB du modèle, la valeur G(t), respectivement EE(t), respectivement I0B(t), 30 respectivement COB(t) correspond à la valeur de la variable à l'instant t, et la valeur G(t-j), respectivement EE(t-j), respectivement I0B(t-j), respectivement COB(t-j), avec j entier, correspond à la valeur de la variable à l'instant t-j*T. B16770 - DD18479 19 A titre d'exemple, la période d'échantillonnage du système est comprise dans la plage allant de 0,5 à 30 minutes, par exemple de l'ordre de 10 minutes. Dans cet exemple, les grandeurs na, nb, nc et nd sont 5 les ordres du modèle, et les grandeurs nki, nk2 et nk3 sont les délais du modèle. La variable d'entrée IOB du modèle peut être calculée à partir des données d'insuline injectée I(t) et d'une fonction d'action hioB représentant l'évolution, en fonction du temps, de 10 l'effet de l'insuline injectée sur l'absorption du glucose sanguin. A titre d'exemple, la variable d'entrée IOB du modèle est définie comme suit: 10B(t) = EL01(t - k) * hroB(k) (eq3) où K est un entier supérieur à 1 et k est un entier allant de 15 0 à K. La variable d'entrée COB du modèle peut être calculée à partir des données de glucose ingéré CHO(t) et d'une fonction d'action hcoB représentant l'évolution, en fonction du temps, de l'effet du glucose ingéré sur la glycémie (ou taux de glucose 20 sanguin). A titre d'exemple, la variable d'entrée COB du modèle est définie comme suit : COB(t) =ELoCHO(t - k)* hcoB(k) (eq4) Dans les équations eq3 et eq4 susmentionnées, la grandeur K définit la durée d de prise en compte de l'historique 25 d'insuline injectée ou de glucose ingéré pour le calcul des variables d'entrée I0B(t) et COB(t) du modèle, telle que d = K*T. A titre d'exemple, la durée d est comprise entre 500 et 2500 minutes, par exemple de l'ordre de 1500 minutes, par exemple égale à 1440 minutes. 30 Les fonctions d'action hioB et hcoB sont par exemple définies comme suit : r hroB(t) = [11- -le TIOB TIOB (eq5) B16770 - DD18479 20 (eq6) hcoB(t) = [1 +-t e TCOB TCOB OÙ TIOB est une constante de temps, par exemple comprise entre 5 et 120 minutes, par exemple de l'ordre de 50 minutes, et Où TCOB est une constante de temps, par exemple différente de TioB, 5 par exemple comprise entre 5 et 120 minutes, par exemple de l'ordre de 40 minutes. La variable d'entrée EE(t) du modèle peut être calculée à partir des signaux de sortie du capteur de mouvements 405 et du capteur de rythme cardiaque 407 (figure 4) du système. A titre 10 d'exemple, le capteur de mouvements 405 fournit un signal Scpm représentatif des mouvements effectués par le patient, par exemple représentatif d'un nombre de mouvements par minute effectués par le patient, par exemple représentatif du nombre de pas par minute effectués par le patient, et le capteur de rythme cardiaque 407 15 fournit un signal SHR représentatif du rythme cardiaque (nombre de pulsations par minute) du patient. Plus particulièrement, dans cet exemple, le signal SHR est égal à la différence entre le rythme cardiaque (nombre de pulsations par minute) du patient à l'instant t et le rythme cardiaque au repos du patient (qui peut 20 être une constante du système). La variable d'entrée EE(t) peut alors être calculée comme suit : lai * SHR(t) + lel, SHR Si EE(t) = a2 * SLc(t) + P2, SHR S2 où le signal SLc est une combinaison linéaire des 25 signaux SHR et Scmp, telle que SLc(t) = el*Scmp(t)+02 *SHR(t), et où les grandeurs al, a2, a3, pi, 132, 133, 01, 02, Si et S2 sont des paramètres du système. A titre d'exemple, les signaux EE(t), I0B(t) et COB(t) peuvent être calculés à une première période d'échantillonnage, 30 par exemple de l'ordre de 1 minute, puis échantillonnés à une période T plus élevée, par exemple de l'ordre de 10 minutes. B16770 - DD18479 21 A titre d'exemple, la valeur al est comprise entre 1 et 10, par exemple de l'ordre de 5,45, la valeur cc2 est comprise entre 100 et 500, par exemple de l'ordre de 256,1, la valeur a3 est comprise entre 10 et 200, par exemple de l'ordre de 86, la 5 valeur 01 est comprise entre -200 et 0, par exemple de l'ordre de -66,1, la valeur 02 est comprise entre -10 et 10, par exemple de l'ordre de -0,13, et la valeur E33 est comprise entre 10 et 200, par exemple de l'ordre de 82,4. A titre d'exemple, le seuil Si est compris entre 20 et 60, par exemple de l'ordre de 40, et le 10 seuil S2 est compris entre 0,1 et 1, par exemple de l'ordre de 0,5. Les valeurs 01 et 02 sont par exemple respectivement de l'ordre de 0,000084 et de l'ordre de 0,0086. Un exemple de procédé d'estimation de la dépense énergétique d'un individu, pouvant être utilisé pour calculer la variable d'entrée EE(t) du modèle, est 15 décrit plus en détail dans l'article intitulé "An original piecewise model for computing energy expenditure from accelerometer and heart rate signals" de H.M. Romero-Ugalde et al. (Physiological Measurement, vol. 38, no. 8, p. 1599, 2017). La figure 6 est un diagramme illustrant un exemple d'un 20 procédé de construction du modèle de la figure 5. Plus particulièrement, étant donnée la structure de modèle décrite ci-dessus et définie par l'équation eq2, la figure 6 illustre un exemple de procédé de dimensionnement du modèle, c'est-à-dire de détermination des ordres na, nb, nc et nd et des délais nki, nk2 25 et nk3 du modèle. Pour dimensionner le modèle, on utilise ici une base de données contenant, pour un ou plusieurs patients, un historique de l'activité physique PA(t) (par exemple déterminé à partir d'un historique des données de sortie d'un accéléromètre porté par le patient et d'un capteur du rythme cardiaque du 30 patient), des prises de glucose CHO(t), des prises d'insuline I(t), et de la glycémie réelle (mesurée par un capteur de glycémie) Gr(t) du patient. Le procédé de la figure 6 comprend une étape 601 au cours de laquelle on affecte, par exemple arbitrairement, une B16770 - DD18479 22 valeur à chacune des dimensions na, nb, nc, nd, nki, nk2 et nk3 du modèle. Le procédé de la figure 6 comprend ensuite une étape 603 de calcul des paramètres du modèle, à savoir les valeurs al, ..., 5 ana, bl, bnb, cl, ---, cnc, dl, ---, dnd- Lors de cette étape, la glycémie du patient est estimée à partir du modèle sur une période d'observation passée, en tenant compte de l'historique d'injection d'insuline, de prise de glucose, et d'activité physique, puis la courbe de glycémie G estimée par le modèle est 10 comparée à la courbe de glycémie réelle Gr mesurée par le capteur pendant cette même période. On recherche alors un jeu de paramètres al, ..., ana, bl, bnb, cl, cnc, dl, dnd 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 15 glycémie réelle pendant la période d'observation. A titre d'exemple, on recherche un jeu de paramètres conduisant à minimiser un indicateur 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 20 moyen entre la glycémie estimée et la glycémie réelle, par exemple identique ou similaire à l'indicateur m de l'équation eql susmentionné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 25 les algorithmes usuels utilisés dans des domaines variés pour résoudre des problèmes d'optimisation de paramètres par minimisation d'une fonction de coût, par exemple par la méthode des moindres carrés. On notera que dans le cas d'une base de données comportant les données d'historique de plusieurs patients, 30 l'étape 603 peut comprendre les étapes successives suivantes : - détermination d'une matrice de régression pour chaque patient à partir des données d'historique du patient ; - concaténation des matrices de régression des 35 différents patients ; et B16770 - DD18479 23 - détermination des paramètres optimaux du modèle conduisant à minimiser l'erreur entre la courbe de glycémie estimée par le modèle et la courbe de glycémie réelle. 5 Le procédé de la figure 6 comprend en outre, après l'étape 603, une fois les paramètres optimaux (pour les dimensions fixées à l'étape 601) du modèle déterminées, une étape 605 de détermination de l'erreur entre la courbe de glycémie estimée par le modèle et la courbe de glycémie réelle, sur la période 10 d'observation passée considérée. A titre d'exemple, lors de cette étape, on mémorise l'indicateur d'erreur calculé à l'étape 603 pour le jeu de paramètres considéré. Le procédé de la figure 6 est réitéré une pluralité de fois en modifiant à chaque fois une ou plusieurs dimensions (ordre 15 ou délai) du modèle. A titre d'exemple, l'ordre na est modifié de façon à balayer toutes les valeurs de la plage [1,15], l'ordre nb est modifié de façon à balayer toutes les valeurs de la plage [1,15], l'ordre nc est modifié de façon à balayer toutes les valeurs de la plage [1,15], l'ordre nd est modifié de façon à 20 balayer toutes les valeurs de la plage [1,15], le délai nid_ est modifié de façon à balayer toutes les valeurs de la plage [1,15], le délai nk2 est modifié de façon à balayer toutes les valeurs de la plage [1,15], et le délai nk3 est modifié de façon à balayer toutes les valeurs de la plage [1,15]. Ainsi, dans cet exemple, 25 la suite d'étapes 601, 602 et 603 est réitérée N = 15*15*15*15*15*15*15 = 170859375 fois. A l'issue des N itérations de la suite d'étapes 601, 602, 603, les dimensions et le jeu de paramètres correspondant, conduisant à la plus faible erreur entre la courbe de glycémie 30 estimée par le modèle et la courbe de glycémie réelle, sur la période d'observation passée, sont sélectionnés. A titre d'exemple, les dimensions finales du modèle sélectionnées à l'issue du procédé de la figure 6 sont les suivantes : na = 3, nb - 5, nc - 11, nd - 7, nkl= nk2 = nk3 - 2- B16770 - DD18479 24 Un avantage du modèle de type boîte noire décrit en relation avec la figure 5 est qu'il permet de réaliser des prédictions particulièrement fiables de l'évolution future de la glycémie du patient, ce que ne permettent pas les modèles 5 physiologiques connus. En effet, s'il est bien connu que l'insuline favorise l'absorption du glucose présent dans le sang par différents organes du corps, et notamment par les muscles, et que la sensibilité à l'insuline et l'absorption du glucose sanguin par les muscles augmentent en présence d'une activité physique, 10 la quantification des effets combinés de ces différentes variables est très difficile à réaliser. L'utilisation d'un modèle de type boîte noire ayant pour entrées la dépense énergétique, l'insuline embarquée, et le glucose embarqué, permet avantageusement d'obtenir des prédictions de glycémie particulièrement fiables, 15 notamment pendant et après une activité physique. Le système de prédiction de la figure 4 peut être utilisé pour de nombreuses applications. A titre d'exemple, le système de la figure 4 peut être utilisé pour mettre en oeuvre un système automatisé de régulation 20 de glycémie de type MPC tel que décrit en relation avec les figures 1 à 3, en remplaçant le modèle physiologique 201 de la figure 2 par le modèle boîte noire 501 de la figure 5. A titre de variante, le système de prédiction de la figure 4 peut être utilisé pour mettre en oeuvre un système 25 automatisé de régulation de glycémie de type MPC basé sur un modèle physiologique, tel que décrit en relation avec les figures 1 à 3, le modèle boîte noire 501 de la figure 5 étant utilisé en redondance du modèle physiologique 201, pour contrôler la fiabilité des prédictions réalisées à partir du modèle 30 physiologique. La figure 7 un diagramme illustrant un exemple d'un procédé automatisé de régulation de glycémie de type MPC basé sur un modèle physiologique, tel que décrit en relation avec les figures 1 à 3, dans lequel le modèle boîte noire 501 de la figure B16770 - DD18479 25 5 est utilisé pour contrôler la fiabilité des prédictions réalisées à partir du modèle physiologique. Le procédé de la figure 7 comprend une étape 701 de recalibration ou mise à jour du modèle physiologique, identique ou 5 similaire à l'étape 301 du procédé de la figure 3. Lors de cette étape, l'unité de traitement et de contrôle 401 met en oeuvre un procédé de ré-estimation des paramètres temps-dépendants du modèle physiologique en tenant compte des données d'insuline effectivement injectée par le dispositif 103 et des données de 10 glycémie réelle mesurée par le capteur 101 pendant une période d'observation passée. Plus particulièrement, lors de l'étape de calibration, l'unité de traitement et de contrôle 401 simule le comportement du patient sur la période d'observation passée à partir du modèle physiologique (en tenant compte des éventuelles 15 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 pendant cette même période. L'unité de traitement et de contrôle 401 recherche alors, pour les paramètres temps-dépendants du modèle, 20 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 d'observation. Lors de l'étape 701, outre les paramètres temps-dépendants du modèle physiologique, l'unité de traitement et de 25 contrôle 401 peut définir un vecteur d'états initiaux des variables d'état du modèle. Le procédé de la figure 7 comprend en outre, après l'étape 701, une étape 703 de prédiction, par l'unité de traitement et de contrôle 401, de l'évolution temporelle de la 30 glycémie du patient sur une période de prédiction à venir de durée Tpred Lors de cette étape, deux courbes de glycémie à venir G1 et G2 sont déterminées, la première (la courbe Gl) correspondant à la prédiction de glycémie réalisée à partir du modèle physiologique sur la période de prédiction à venir, de façon 35 identique ou similaire à ce qui a été décrit en relation avec la B16770 - DD18479 26 figure 3, et la deuxième (la courbe G2) correspondant à la prédiction de glycémie réalisée à partir du modèle boîte noire de la figure 5 sur la même période de prédiction à venir. Le procédé de la figure 7 comprend de plus, après l'étape 5 703, une étape 705 de comparaison entre la courbe de glycémie à venir G1 estimée à partir du modèle physiologique et la courbe de glycémie à venir G2 estimée à partir du modèle boîte noire. A titre d'exemple, lors de cette étape, l'unité de traitement et de contrôle 401 calcule un ou plusieurs indicateurs numériques 10 représentatifs de l'écart entre les courbes de glycémie prédites G1 et G2, par exemple l'écart quadratique moyen entre les courbes de glycémie prédites G1 et G2. Le procédé de la figure 7 comprend en outre, après l'étape 705, une étape 707 de vérification de la fiabilité du 15 modèle physiologique mis à jour à l'étape 701. Plus particulièrement, lors de cette étape, la fiabilité du modèle physiologique est considérée comme suffisante par l'unité de traitement et de contrôle 401 lorsque la(les) valeur(s) du ou des indicateurs d'écart entre la courbe de glycémie à venir G1 estimée 20 à partir du modèle physiologique et la courbe de glycémie à venir G2 estimée à partir du modèle boîte noire est(sont) inférieure(s) à un(des) seuil(s) prédéfini(s), et insuffisante dans le cas contraire. Si le modèle physiologique est considéré comme 25 suffisamment fiable à l'étape 707 (0), une étape 709 est mise en oeuvre, au cours de laquelle l'unité de traitement et de contrôle 401 détermine, en tenant compte de la courbe de glycémie future G1 estimée à partir du modèle physiologique, des doses d'insuline à injecter au patient pendant la période de prédiction à venir. A 30 l'issue de cette étape, l'unité de traitement et de contrôle 401 peut programmer le dispositif d'injection 103 pour administrer les doses déterminées pendant la période de prédiction à venir. Si le modèle physiologique est jugé insuffisamment fiable à l'étape 707 (N), l'unité de traitement et de contrôle 35 401 cesse d'utiliser le modèle physiologique pour réguler B16770 - DD18479 27 l'administration de l'insuline au patient, et met en oeuvre une méthode de régulation de substitution lors d'une étape 711. A titre d'exemple, lors de l'étape 711, l'unité de traitement et de contrôle 401 utilise la courbe de glycémie future G2 estimée à 5 partir du modèle boîte noire de la figure 5 pour déterminer les doses d'insuline à injecter au patient pendant la période de prédiction à venir. A titre de variante, lors de l'étape 711, l'unité de traitement et de contrôle 401 utilise un modèle physiologique simplifié, par exemple un modèle compartimentai 10 comportant un nombre de variables d'état et un nombre de paramètres réduits par 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 711, l'unité de traitement et de contrôle 401 cesse de mettre en oeuvre 15 une commande prédictive, c'est-à-dire qu'elle cesse d'utiliser un modèle mathématique 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 401 commande par exemple le dispositif d'injection d'insuline 103 pour administrer des doses 20 préprogrammées d'insuline, correspondant par exemple à un débit basal de référence prescrit au patient. Alternativement, l'unité de traitement et de contrôle 401 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 25 tels que le niveau courant de glycémie mesuré par le capteur 101, ou encore la vitesse de variation (ou pente) de la glycémie sur une période passée. Les étapes 701, 703, 705, 707 et, le cas échéant, 709 ou 711, sont par exemple réitérées périodiquement, par exemple 30 toutes les 1 à 20 minutes, et/ou à 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. B16770 - DD18479 28 La figure 8 est un diagramme illustrant un autre exemple d'un procédé automatisé de contrôle de glycémie pouvant être mis en oeuvre par le système de la figure 4. Dans cet exemple, l'unité de traitement et de contrôle 5 401 peut mettre en oeuvre un procédé de régulation de glycémie quelconque, par exemple un procédé de régulation en boucle fermée du type décrit en relation avec les figures 1 à 3, ou un procédé de régulation en boucle ouverte dans lequel le patient définit lui-même les doses d'insuline à injecter sur la base de son taux 10 de glycémie, de son historique d'injection d'insuline, et de son historique de prise de repas. Dans l'exemple de la figure 8, le modèle boîte noire de la figure 5 est utilisé par l'unité de traitement et de contrôle 401 pour prédire l'évolution future de la glycémie du patient, et 15 déclencher une alarme, par exemple une alarme sonore et/ou visuelle, lorsque le modèle prédit que la glycémie du patient va sortir d'une plage souhaitée, par exemple s'il existe un risque que le patient tombe en hypoglycémie. Le procédé de la figure 8 comprend une étape 801 de 20 prédiction, par l'unité de traitement et de contrôle 401, à partir du modèle boîte noire de la figure 5, de l'évolution temporelle de la glycémie du patient sur une période de prédiction à venir. Le procédé de la figure 8 comprend en outre, après l'étape 801, une étape 803 au cours de laquelle l'unité de 25 traitement et de contrôle 401 détermine si la glycémie à venir prédite par le modèle de la figure 5 présente ou non un risque pour le patient. Pour cela, l'unité de traitement et de contrôle 401 vérifie que la courbe de glycémie à venir prédite à l'étape 801 ne sort pas d'une plage souhaitée. 30 Si, à l'étape 803, la glycémie prédite est jugée acceptable (0), le procédé se termine lors d'une étape 805. Si en revanche, à l'étape 803, la glycémie prédite est jugée inacceptable ou dangereuse pour le patient (N), l'unité de traitement et de contrôle 401 déclenche une alarme lors d'une 35 étape 807, de façon à avertir le patient du risque d'hyperglycémie B16770 - DD18479 29 ou d'hypoglycémie, par exemple pour lui recommander de cesser toute activité physique et/ou de consommer du sucre en cas de risque d'hypoglycémie. Les étapes 801, 803 et, le cas échéant, 805 ou 807, sont 5 par exemple réitérées périodiquement. Des modes de réalisation particuliers ont été décrits. Diverses variantes et modifications apparaîtront à l'homme de l'art. En particulier, les modes de réalisation décrits ne se limitent pas aux exemples d'utilisation du système de prédiction 10 de glycémie de la figure 4 décrits en relation avec les figures 7 et 8. Plus généralement, le système de prédiction proposé, basé sur l'utilisation du modèle boîte noire de la figure 5, peut être utilisé pour toute application pouvant tirer profit d'une prédiction fiable de l'évolution future de la glycémie d'un 15 patient. De plus, les modes de réalisation décrits ne se limitent pas à l'exemple décrit ci-dessus dans lequel la variable d'entrée EE(t) du modèle est calculée en tenant compte d'un signal de sortie d'un capteur des mouvements du patient et d'un signal de 20 sortie d'un capteur du rythme cardiaque du patient. A titre de variante, la variable d'entrée EE(t) peut être calculée en tenant compte uniquement du signal de sortie du capteur de mouvements ou uniquement du signal de sortie du capteur de rythme cardiaque. Plus généralement, la variable d'entrée EE(t) peut être constituée 25 par tout signal ou combinaison de signaux représentatifs de l'activité physique du patient. A titre d'exemple, la variable d'entrée EE(t) peut être une combinaison d'un signal mesuré au moyen d'un capteur quelconque, par exemple le capteur de mouvements 405 et/ou le capteur de rythme cardiaque 407 de la 30 figure 4, et d'un signal d'intensité d'activité physique déclaré par le patient au moyen d'une interface utilisateur (non détaillée). A titre de variante, la variable d'entrée EE(t) peut être calculée uniquement à partir d'un signal d'intensité d'activité déclaré par le patient. B16770 - DD18479 30 Par ailleurs, les modes de réalisation décrits ne se limitent pas au cas particulier décrit ci-dessus dans lequel le modèle boîte noire utilisé pour prédire la glycémie future du patient est un modèle de type ARX. Plus généralement tout autre 5 type de modèle boîte noire peut être utilisé, par exemple un modèle de type ARMAX (de l'anglais "Auto Regressive Moving Average with eXternal inputs - modèle autorégressif à moyenne mobile avec des entrées externes), ou un modèle basé sur un réseau neuronal. De plus, outre les variables d'entrée EE(t), I0B(t) et 10 COB(t) décrites ci-dessus, le modèle de la figure 5 peut comprendre une variable d'entrée supplémentaire représentative de la variabilité du rythme cardiaque du patient, REVENDICATIONS1. Système automatisé de contrôle de la glycémie d'un patient, comportant une unité de traitement et de contrôle (401) adaptée à prédire, à partir d'un premier modèle mathématique (501), l'évolution future de la glycémie G(t) du patient sur une période de prédiction, dans lequel le premier modèle prend pour entrées : une variable EE(t) représentative de l'évolution, en fonction du temps, de la dépense énergétique du patient ; une variable I0B(t) représentative de l'évolution, en 10 fonction du temps, de la quantité d'insuline embarquée du patient ; et une variable COB(t) représentative de l'évolution, en fonction du temps, de la quantité de glucose embarquée du patient, dans lequel le premier modèle mathématique (501) est un 15 modèle de type boîte noire, c'est-à-dire un modèle non physiologique défini par la seule observation des effets des variables d'entrée EE(t), I0B(t) et COB(t) sur la variable de sortie G(t), sans tenir compte des différents mécanismes physiologiques connus s'opérant dans le corps du patient, 20 et dans lequel l'unité de traitement et de contrôle (401) est configurée pour prédire, à partir du premier modèle mathématique (501), l'évolution future de la glycémie G(t) du patient sur une période de prédiction, déterminer si la glycémie G(t) prédite par le premier modèle reste ou non contenue dans une 25 plage souhaitée, et déclencher une alarme si la glycémie G(t) prédite par le premier modèle sort de la plage souhaitée. 2. Système selon la revendication 1, dans lequel le premier modèle mathématique (501) est défini comme suit : G(t) = a.yG+ b.uEE + c.ujoE + 30 où t est une variable temps discrétisée, a = [al, - ana] est unvecteur de paramètres de dimension na, avec na entier supérieur ou égal à 1, b = [bl, bnb] est un vecteur de B16770 - DD18479 32 paramètres de dimension nb, avec nb entier supérieur ou égal à 1, c = [cl, cnc] est un vecteur de paramètres de dimension nc, avec nc entier supérieur ou égal à 1, d = [dl, dnd] est un vecteur de paramètres de dimension nd, avec nd entier supérieur 5 ou égal à 1, yG = [G(t-1), G(t-na)] est un vecteur de régression de dimension na, uEE = [EE(t -nkl), EE (t-nkl-nb)] est un vecteur de régression de dimension nb, avec nu entier supérieur ou égal à 1, umg [I0B(t -nk2), I0B(t -nk2-nc)] est un vecteur de régression de dimension nc, avec nk2 entier supérieur ou égal à 1, et ucoE =[COB(t-nk3), COB(t-nk3-nd)] est un vecteur de régression de dimension nd, avec nk3 entier supérieur ou égal à 1. 3. Système selon la revendication 2, dans lequel na est compris dans la plage [1,15], nb est compris dans la plage [1,15], nc est compris dans la plage [1,15], nd est compris dans la plage [1,15], nkl est compris dans la plage [1,15], nk2 est compris dans la plage [1,15], et nk3 est compris dans la plage [1,15]. 4. Système selon l'une quelconque des revendications 1 à 3, dans lequel, en fonctionnement, l'unité de traitement et de contrôle (401) reçoit un signal CHO(t) représentatif de l'évolution, en fonction du temps, de la quantité de glucose ingérée par le patient, et un signal I(t) représentatif de l'évolution, en fonction du temps, de la quantité d'insuline 25 injectée au patient. 5. Système selon la revendication 4, dans lequel l'unité de traitement et de contrôle (401) est configurée pour calculer les variables d'entrée I0B(t) et COB(t) du premier modèle (501) conne suit : 30 I 0 B (t) Elk(=0 i(t - k) * 1-110B(k) , et COB (t) = Elkc CHO (t - k) * hco B (k) , où K est un entier supérieur à 1, k est un entier allant de 0 à K, hioe est une fonction d'action représentant l'évolution, B16770 - DD18479 33 en fonction du temps, de l'effet de l'insuline injectée sur l'absorption du glucose sanguin, et hcoB est une fonction d'action représentant l'évolution, en fonction du temps, de l'effet du glucose ingéré sur la glycémie. 5 6. Système selon la revendication 5, dans lequel les fonctions d'actions hioB et hcoB sont définies comme suit : r (t) [1 + e 'JOB, et nos hCOB (t) = 1- e TCOB r TCOB où l'OB et tCOB sont des constantes de temps. 10 7. Système selon l'une quelconque des revendications 1 à 6, comprenant en outre un dispositif (403) de mesure d'une activité physique du patient, la variable d'entrée EE(t) du premier modèle (501) étant calculée par l'unité de traitement et de contrôle (401) à partir de données de sortie PA(t) du dispositif 15 de mesure (403). 8. Système selon la revendication 7, dans lequel le dispositif de mesure (403) comprend un capteur (405) de mouvements du patient. 9. Système selon la revendication 8, dans lequel le 20 dispositif de mesure (403) comprend en outre un capteur (407) du rythme cardiaque du patient. 10. Système selon la revendication 9, dans lequel le capteur de mouvements (405) fournit un signal Scpm représentatif de mouvements effectués par le patient, et le capteur de rythme cardiaque (407) fournit un signal SHR représentatif du rythme cardiaque du patient, et dans lequel la variable d'entrée EE(t) du premier modèle (501) est calculée comme suit par l'unité de traitement et de contrôle (401) : lai * SHR(t) + 161, SHR > S1 EE(t) = a2 * Syc(t) + )62, SHR 11. Système selon l'une quelconque des revendications 5 1 à 10, comportant en outre un capteur de glycémie (101), et un dispositif d'injection d'insuline (103), dans lequel l'unité de traitement et de contrôle (401) est adaptée à prédire l'évolution future de la glycémie du patient à partir d'un deuxième modèle mathématique (201) et à commander 10 le dispositif d'injection d'insuline (103) en tenant compte de cette prédiction, l'unité de traitement et de contrôle (401) étant configurée pour : a) mettre en oeuvre une étape de calibration automatique 15 du deuxième modèle (201) en tenant compte d'un historique de 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, réaliser une première prédiction de l'évolution future de la glycémie du 20 patient à partir du premier modèle (501) et, sur la même période de prédiction, une deuxième prédiction de l'évolution future de la glycémie du patient à partir du deuxième modèle (201) ; c) calculer au moins un indicateur numérique représentatif de l'écart entre les première et deuxième 25 prédictions ; et d) en fonction de la valeur de l'indicateur numérique, commander le dispositif d'injection d'insuline (103) en tenant compte de la deuxième prédiction ou commander le dispositif d'injection d'insuline (103) sans tenir compte de la deuxième 30 prédiction. 12. Système selon la revendication 11, dans lequel le deuxième modèle (201) est un modèle physiologique comportant un système d'équations différentielles décrivant l'évolution d'une pluralité de variables d'état en fonction du temps. B16770 - DD18479 35 13. Système selon la revendication 12, dans lequel l'étape a) de calibration automatique du deuxième modèle (201) comprend une étape d'estimation de paramètres du système d'équations différentielles par minimisation d'une grandeur 5 représentative de l'erreur, pendant la période d'observation passée, entre la glycémie estimée à partir du deuxième modèle physiologique et la glycémie mesurée par le capteur (101).