FR 3081316 A1 20191129 FR 1800493 A 20180522 6895 - DD18588 1 SYSTEME AUTOMATISE DE CONTROLE DE LA GLYCEMIE D'UN PATIENT Domaine La présente demande concerne le domaine des systèmes automatisés de contrôle de glycémie, et vise plus particulièrement, dans un tel système, la détermination d'un coefficient représentatif de la sensibilité à l'insuline du 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) déposée le 21 septembre 2016, et dans la demande de brevet français N°1756960 (B15860/DD18588) 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 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 à commande prédictive, dans lesquels la régulation de la dose d'insuline administrée tient compte d'une prédiction de 316895 - DD18588 2 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. Plus généralement, de nombreux systèmes automatisés de 5 contrôle de glycémie tiennent compte de l'historique de glycémie, de l'historique de prise de repas, et de l'historique d'injection d'insuline d'un patient pour déterminer des doses d'insuline à administrer au patient en vue de maintenir sa glycémie dans une plage souhaitée. 10 Dans les systèmes automatisés de contrôle de glycémie, un paramètre qui joue un rôle essentiel dans la détermination des doses d'insuline à administrer au patient est le coefficient de sensibilité à l'insuline du patient, aussi appelée coefficient de sensibilité de compensation, ou ratio de compensation, c'est-à- 15 dire la quantité d'insuline nécessaire pour faire baisser la glycémie de un gramme par litre (en UI/g/1 - où UI désigne une unité internationale d'insuline, soit l'équivalent biologique d'environ 0,0347 mg d'insuline humaine). Un problème qui se pose est que le coefficient de 20 sensibilité à l'insuline peut varier de façon significative d'un patient à un autre, ou, chez un même patient, en fonction des conditions dans lesquelles se trouve le patient, et notamment en fonction de l'activité physique du patient. En pratique, les systèmes automatisés de contrôle de glycémie connus se basent sur un coefficient de sensibilité à l'insuline fixe, par exemple fourni par le diabétologue du patient. Il n'est donc pas tenu compte des variations à court terme du coefficient. Il en résulte que les quantités d'insuline injectées au patient sont parfois inadaptées, entrainant un risque d'hyperglycémie ou d'hypoglycémie. 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 un capteur de glycémie, un dispositif de mesure d'une activité B16895 - DD18588 3 physique du patient, et une unité de traitement et de contrôle, dans lequel : l'unité de traitement et de contrôle est configurée pour générer, par convolution d'un signal PA fourni par le dispositif 5 de mesure d'une activité physique du patient avec une fonction mathématique décroissante H, un signal IPA représentatif de l'influence de l'activité physique du patient sur sa sensibilité à l'insuline ; et l'unité de traitement et de contrôle est configurée pour 10 calculer, à partir d'un premier modèle mathématique fCR spécifique au patient et en tenant compte du signal IPA et d'une unique valeur de glycémie Gr mesurée par le capteur, un coefficient CR représentatif de la sensibilité à l'insuline du patient. Selon un mode de réalisation, le système comporte en 15 outre un dispositif d'injection d'insuline, et l'unité de traitement et de contrôle est configurée pour commander le dispositif d'injection d'insuline en tenant compte du coefficient CR. Selon un mode de réalisation, l'unité de traitement et 20 de contrôle est configurée pour prédire, à partir d'un deuxième modèle mathématique, 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. Selon un mode de réalisation, le premier modèle 25 mathématique est une fonction d'équation CR = fcR(IPA,Gr) = a x IPAb + c x Grd + e où a, b, c, d et e sont des paramètres spécifiques au patient. Selon un mode de réalisation, le premier modèle 30 mathématique est une fonction d'équation CR = Gr x fcR(IPA) = Gr X ( a x IPAb + c) où a, b et c sont des paramètres spécifiques au patient. Selon un mode de réalisation, l'unité de traitement et de contrôle est configurée pour mettre en oeuvre une étape de 316895 - DD18588 4 calibration automatique du premier modèle foR en tenant compte d'un historique de glycémie mesurée par le capteur, d'un historique d'insuline injectée au patient, d'un historique d'ingestion de glucose par le patient, et d'un historique du 5 signal d'activité physique PA du patient, sur une période d'observation passée. Selon un mode de réalisation, l'unité de traitement et de contrôle est configurée pour, lors de l'étape de calibration automatique, mesurer une pluralité de valeurs du coefficient de 10 sensibilité à l'insuline réel CRr du patient lors d'une pluralité d'évènements de mesure compris dans la période d'observation passée. Selon un mode de réalisation, chaque évènement de mesure correspond à une plage temporelle continue allant d'un instant 15 initial tinit à un instant final tfinal, répondant aux critères suivants : - l'instant tinit se situe dans une phase d'hyperglycémie, c'est-à-dire une phase dans laquelle la glycémie du patient est supérieure à un seuil 20 prédéterminé ; - un bolus de correction, c'est-à-dire une dose d'insuline a été administré au patient après le début de la phase d'hyperglycémie et avant l'instant tinit, en vue de limiter la durée de la phase 25 d'hyperglycémie ; - la glycémie du patient décroit de façon continue entre l'instant initial tinit et l'instant final tfinal ; et - aucune ingestion de glucose par le patient n'a été 30 effectuée entre l'instant tinit - Ti et l'instant tfinal, où Tj est une durée de jeun prédéterminée. Selon un mode de réalisation, l'unité de traitement et de contrôle est configurée pour, lors de l'étape de calibration automatique, déterminer le premier modèle mathématique fCR par B16895 - DD18588 5 régression à partir de ladite pluralité de valeurs du coefficient de sensibilité à l'insuline réel CRr. Selon un mode de réalisation, la fonction H est une fonction exponentielle décroissante. 5 Selon un mode de réalisation, le dispositif de mesure comprend une interface utilisateur par laquelle le patient déclare les activités physiques qu'il réalise. Selon un mode de réalisation, le dispositif de mesure comprend un ou plusieurs capteurs adaptés à mesurer des grandeurs 10 représentatives de l'activité physique du patient. Selon un mode de réalisation, le dispositif de mesure comprend un capteur de mouvement et/ou un capteur de rythme cardiaque. Brève description des dessins 15 Ces caractéristiques et 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 20 de blocs, un exemple d'un système automatisé de régulation de la glycémie d'un patient selon un mode de réalisation ; 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 ; 25 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 ; les figures 4A et 4B sont des diagrammes représentant l'évolution, en fonction du temps, respectivement d'un signal PA 30 représentatif de l'activité physique du patient et d'un signal IPA, représentatif de l'influence de l'activité physique du patient sur son coefficient de sensibilité à l'insuline ; les figures 5A et 5B sont des diagrammes représentant l'évolution du coefficient de sensibilité à l'insuline d'un 35 patient en fonction de sa glycémie et du signal IPA ; B16895 - DD18588 6 la figure 6 est une représentation simplifiée d'un modèle mathématique utilisé dans le système de la figure 1 pour déterminer le coefficient de sensibilité à l'insuline du patient ; et 5 la figure 7 est un diagramme illustrant un exemple d'un procédé pouvant être mis en oeuvre par le système de la figure 1 pour mettre à jour le modèle mathématique de la figure 6. Description détaillée De mêmes éléments ont été désignés par de mêmes 10 références dans les différentes figures et, de plus, les diverses figures ne sont pas tracées à l'échelle. Par souci de clarté, 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 15 et de contrôle des systèmes décrits 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. En outre, le dispositif de mesure de glycémie et le dispositif d'injection d'insuline des 20 systèmes décrits 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. Sauf précision contraire, les expressions "approximativement", "sensiblement", et "de l'ordre de" signifient à 10 % près, de 25 préférence à 5 % près. 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) 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 316895 - DD18588 7 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. Dans cet exemple, l'unité de traitement 25 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. L'unité de traitement et de contrôle 105 est adaptée à 30 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. Pour cela, l'unité de traitement et de contrôle 105 comprend un 35 circuit de calcul numérique (non détaillé), comprenant par exemple B16895 - DD18588 8 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 smartphone configuré pour mettre en oeuvre un procédé de 5 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 l'évolution future de sa glycémie en fonction du temps. Plus 10 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 mathématique, par exemple un modèle physiologique décrivant l'assimilation de l'insuline par le corps du patient et son impact 15 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 et d'un coefficient représentatif 20 de la sensibilité à l'insuline du patient, l'unité de traitement et de contrôle 105 détermine les doses d'insuline qu'il convient 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 mathématique) du patient reste dans des limites 25 acceptables, et en particulier pour limiter les risques d'hyperglycémie ou d'hypoglycémie. La figure 2 est une représentation simplifiée d'un modèle mathématique 201 (MPC) utilisé dans le système de la figure 1 pour prédire l'évolution future de la glycémie du patient. Sur 30 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 quantité d'insuline injectée au patient ; B16895 - DD18588 9 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 si fournissant un signal G(t) représentatif 5 de l'évolution, en fonction du temps t, de la glycémie estimée du patient. Le modèle mathématique 201 est par exemple un modèle physiologique. A titre d'exemple, le modèle 201 est un modèle physiologique compartimental comportant, outre les variables 10 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 un système d'équations différentielles 15 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 la figure 2 par un vecteur [INIT] appliqué 20 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é "Nonlinear model predictive control of glucose concentration in subjects with type 1 diabetes" de 25 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, 2002). Plus généralement, tout autre modèle physiologique 30 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é. D'autres paramètres, appelés ci-après paramètres temps-35 dépendants, sont en revanche susceptibles d'évoluer dans le temps. B16895 - DD18588 10 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 les 1 à 20 minutes, par exemple toutes les 5 minutes, pour 5 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 de préférence 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 physiquement les paramètres temps-dépendants du système 10 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 oeuvre par le système de la figure 1. 15 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 oeuvre un procédé de ré-estimation des paramètres temps-dépendants 20 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 précédant l'étape de calibration. Plus particulièrement, lors de 25 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 cette période), et compare la courbe de glycémie estimée par le 30 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 représentative de l'erreur entre la courbe de glycémie estimée 35 par le modèle et la courbe de glycémie réelle pendant la période B16895 - DD18588 11 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 5 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 : t, m A1T z Icr(t)- G(t)I2 t=to-AT où t est une variable temps discrétisée, t0-AT correspond à 10 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 réelle mesurée par le capteur 101 pendant la période [to-AT, -bo], 15 et G est la courbe de glycémie estimée à partir du modèle pendant la période [t0-AT, -Loo]. 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 passée. L'algorithme de recherche de paramètres optimaux utilisé 20 lors de cette étape n'est pas détaillé dans la présente demande, les modes de réalisatioC 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 fonction de coût. 25 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 comportement du patient à partir du modèle. Pour définir les états 30 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 calibration du modèle, le patient se trouvait dans un état B16895 - DD18588 12 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 à 5 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 10 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 des variables d'état sont alors déterminés de la même façon que 15 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 20 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 [to, 25 teTpred] de durée T- pred , 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 30 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 [-Co, to+T- pred] . A l'issue de cette étape, l'unité de traitement et de contrôle 105 peut B16895 - DD18588 13 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 de 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 Selon un aspect d'un mode de réalisation, lors de l'étape 305, l'unité de traitement et de contrôle 105 estime le coefficient de sensibilité à l'insuline CR du patient en tenant compte de l'activité physique du patient. Pour cela, comme cela apparaît sur la figure 1, le système comprend un dispositif 107 15 de mesure d'une activité physique du patient. Le dispositif 107 est relié à l'unité de traitement et de contrôle 105, par exemple par liaison filaire ou par liaison radio (sans fil), et communique à l'unité de traitement et de contrôle 105 un signal PA(t) représentatif de l'évolution, en fonction du temps t, d'une 20 activité physique du patient. Le dispositif 107 est par exemple un dispositif mobile porté par le patient tout au long de la journée et/ou de la nuit. A titre d'exemple, le dispositif 107 est une simple interface utilisateur (non détaillée) par laquelle le patient 25 déclare les activités physiques qu'il réalise. A titre d'exemple, le signal PA(t) correspond à un niveau d'intensité physique déclaré par le patient par le biais du dispositif 107. Le dispositif 107 est par exemple muni d'un clavier permettant à l'utilisateur de saisir, sur une échelle de 0 à N, où N est un 30 entier strictement positif, par exemple égal à 3, le niveau d'intensité de l'activité physique qu'il est en train de réaliser, la valeur 0 correspondant à une activité physique considérée comme nulle ou négligeable, et la valeur N correspondant au niveau d'intensité d'activité physique maximal du patient. B16895 - DD18588 14 A titre de variante, le dispositif 107 comprend un ou plusieurs capteurs adaptés à mesurer des grandeurs représentatives de l'activité physique du patient. A titre d'exemple, le dispositif 107 comprend au moins un capteur de mouvement (non 5 détaillé sur la figure 1), par exemple un accéléromètre. Le dispositif 107 peut en outre comprendre un capteur du rythme cardiaque du patient (non détaillé sur la figure 1). Dans ce cas, le signal PA(t) est par exemple un signal représentatif de la dépense énergétique du patient, calculé à partir des données de 10 sortie du ou des capteurs du dispositif 107, par exemple tel que décrit dans l'article intitulé "Prior automatic posture and activity identification improves physical activity energy expenditure prediction from hip-worn triaxial acceleromêtry" de M. Garnotel et al. (Journal of Applied Physiology (1985). 2017 15 Nov 30), ou dans l'article intitulé "An original piecewise model for computing energy expenditure from accelerometer and heart rate signals" de H. Romêro-Ugalde et al. (Physiological measuremênt, 2017 Jul 28;38(8):1599-1615). A titre de variante, le signal PA(t) peut être une 20 combinaison d'un signal mesuré au moyen d'un ou plusieurs capteurs du dispositif 107, et d'un signal d'intensité d'activité physique déclaré par le patient au moyen d'une interface utilisateur du dispositif 107. Selon un aspect d'un mode de réalisation, l'unité de 25 traitement et de contrôle 105 est configurée pour calculer, à partir du signal PA(t) représentatif de l'évolution, en fonction du temps t, de l'activité physique du patient, un signal IPA(t) représentatif de l'évolution, en fonction du temps t, de l'influence de l'activité physique passée ou en cours sur le 30 coefficient de sensibilité à l'insuline du patient. Plus particulièrement, le signal IPA(t) est généré par convolution du signal PA(t) avec une fonction de décroissance H(t). La fonction H(t) est par exemple une fonction exponentielle décroissante. Plus généralement, toute autre fonction décroissante représentative de 35 la diminution de l'influence d'une activité physique sur le B16895 - DD18588 15 coefficient de sensibilité à l'insuline du patient peut être utilisée, par exemple une fonction linéaire décroissante, une phase de décroissance d'une fonction parabolique ou hyperbolique, etc. La fonction de décroissance H(t) peut prendre en compte 5 l'intensité et la durée de l'activité physique, ainsi qu'un temps caractéristique d'extinction Text, c'est-à-dire une durée au-delà de laquelle l'influence de l'activité physique passée sur le coefficient de sensibilité à l'insuline est considérée comme nulle ou négligeable. 10 A titre d'exemple, la fonction IPA(t) est définie comme suit : I P A(t) = (PA * H)(t) avec : 1 t H (t) =-e-7 15 où T est une variable, en min-1, fonction de la durée et de l'intensité des activités physiques passées. En pratique, le signal PA(t) et la fonction H(t) peuvent être échantillonnés à une période d'échantillonnage T, par exemple de l'ordre de 1 minute. La variable T peut être calculée à chaque 20 instant t, sur une fenêtre temporelle glissante [t-Text ; t]. A titre d'exemple, la variable T est définie comme suit : Next T = 1(k) x (1 "ext k=1 Next étant un entier supérieur à 1 définissant la durée de la période d'extinction Text, telle que Text =Next X T, et I 25 étant un vecteur de Next valeurs représentant l'intensité de l'activité physique effectuée par le patient sur la période [tText ; t], échantillonnée à la période d'échantillonnage T. A titre d'exemple, la période d'extinction Text est comprise entre B16895 - DD18588 16 12 et 72 heures, par exemple de l'ordre de 48 heures, soit Next = 2880 pour T = 1 minute. Les figures 4A et 4B sont des diagrammes représentant, à titre d'exemple illustratif, l'évolution, en fonction du temps, 5 respectivement du signal PA (figure 4A) représentatif de l'intensité d'une activité physique effectuée par le patient, et du signal IPA (figure 4B) représentatif de l'influence de l'activité physique du patient sur son coefficient de sensibilité à l'insuline. Dans cet exemple, le signal IPA est calculé par 10 l'unité de traitement et de contrôle à partir de la fonction de décroissance exponentielle H définie ci-dessus. Selon un aspect d'un mode de réalisation, lors de l'étape 305 du procédé de la figure 3, l'unité de traitement et de contrôle 105 estime le coefficient de sensibilité à l'insuline 15 OR du patient à partir d'une unique valeur du signal IPA et d'une unique valeur de glycémie mesurée par le capteur 101, en se basant sur un modèle mathématique prédéterminé. Plus particulièrement, lors de l'étape 305, l'unité de traitement et de contrôle 105 calcule le coefficient de sensibilité à l'insuline OR estimé du 20 patient à partir d'une fonction mathématique prédéterminée fcR telle que CR =Gr(t)xfcR(IPA(0), où Gr(t) est la valeur de glycémie réelle du patient mesurée par le capteur 101 à un instant t courant, par exemple à l'instant t=to de début de la période de prédiction considérée à l'étape 303, IPA(t) étant la valeur du 25 signal IPA à l'instant t. Lors de l'étape 305, l'unité de traitement et de contrôle 105 détermine alors les doses futures d'insuline à administrer au patient en tenant compte du coefficient de sensibilité OR ainsi calculé. Les inventeurs ont montré qu'il existe, pour un patient 30 donné, une corrélation forte entre l'évolution temporelle du signal IPA du patient, et l'évolution temporelle du coefficient de sensibilité à l'insuline du patient, normalisé par rapport à la glycémie du patient. Les inventeurs ont notamment montré que l'ajustement, en temps réel, du coefficient de sensibilité à 35 l'insuline du patient en fonction de sa glycémie instantanée et B16895 - DD18588 17 du signal IPA, permet de déterminer avec une meilleure précision les doses d'insuline futures à administrer au patient et ainsi limiter les risques d'hyperglycémie ou d'hypoglycémie. Les figures 5A et 5B sont des diagrammes représentant 5 respectivement, pour deux patients distincts, l'évolution du rapport du coefficient de sensibilité à l'insuline du patient divisé par sa glycémie (en ordonnée, en UI/g/l/g/l), en fonction de la valeur de la variable IPA du patient (en abscisse). Chaque diagramme comprend une pluralité de points 501 correspondant 10 chacun à une mesure du rapport Rr = CRr/Gr entre le coefficient de sensibilité à l'insuline réel CRr du patient, et à une valeur correspondante (i.e. corrélée temporellement) de la glycémie réelle Gr du patient. Le coefficient de sensibilité à l'insuline réel CRr peut 15 être mesuré par toute méthode connue de mesure du coefficient de sensibilité à l'insuline d'un patient, par exemple par des méthodes du type décrit dans les demandes de brevet US2010/0198520, US2013/0211220 et W02017/040927. Dans un mode de réalisation préféré, le coefficient de 20 sensibilité à l'insuline réel CRr du patient est déterminé à partir de l'historique de glycémie du patient (par exemple mesuré par le capteur 101 dans le système de la figure 1), de son historique de prise de repas, et de son historique d'injection d'insuline, selon la méthode suivante. 25 A partir de l'historique de données du patient, on identifie des évènements de mesure, c'est-à-dire des plages de temps pendant lesquelles le coefficient de sensibilité à l'insuline est isolé, c'est-à-dire pendant lesquelles on observe une diminution de la glycémie du patient liée à l'administration 30 d'insuline. A titre d'exemple, les évènements sélectionnés sont des plages temporelles continues allant d'un instant initial tinit à un instant final tfinal, répondant aux critères suivants : - l'instant tinit se situe dans une phase d'hyperglycémie, c'est-à-dire une phase dans laquelle B16895 - DD18588 18 la glycémie du patient est supérieure à un seuil prédéterminé, par exemple de l'ordre de 1,40 g/1 ; - un bolus de correction, c'est-à-dire une dose d'insuline supplémentaire a été administré au patient 5 après le début de la phase d'hyperglycémie et avant l'instant tinit, en vue de limiter la durée de la phase d'hyperglycémie ; - la glycémie du patient décroit de façon continue entre l'instant initial tinit et l'instant final tfinal ; 10 - aucune ingestion de glucose par le patient n'a été effectuée entre l'instant tinit - Ti et l'instant tfinal, où Ti est une durée de jeun prédéterminée, par exemple supérieure ou égale à lh et de préférence supérieure ou égale à 2h. 15 A titre d'exemple, l'instant tinit correspond au pic de glycémie de la phase d'hyperglycémie. L'instant tfinal correspond par exemple à un instant de stabilisation ou de remontée glycémique suivant la phase d'hyperglycémie, ou encore à une perturbation telle qu'un repas ou une ingestion de glucose. 20 Pour chaque évènement identifié, le coefficient de sensibilité à l'insuline réel CRr du patient est calculé comme suit : CRr - AI/AG, où AI désigne la quantité d'insuline consommée pendant 25 l'évènement et AG désigne la différence entre la glycémie réelle du patient à l'instant tinit de début de l'évènement et la glycémie réelle du patient à l'instant tfinal de fin de l'évènement. La quantité d'insuline AI consommée pendant l'évènement peut par exemple être calculée en tenant compte des doses d'insuline 30 administrées avant et pendant l'évènement, et de la cinétique d'absorption de l'insuline par le corps. A titre d'exemple, la quantité d'insuline AI consommée pendant l'évènement correspond à la différence entre la quantité d'insuline embarquée ("insulin on board" en anglais) du patient, c'est-à-dire la quantité d'insuline 35 encore active (c'est-à-dire encore susceptible d'avoir un effet B16895 - DD18588 19 sur la glycémie) à l'instant tinit de début de l'évènement et la quantité d'insuline embarquée à l'instant tfinal de fin de l'évènement. La détermination de la quantité d'insuline embarquée du patient aux instants tinit et tfinal peut être réalisée par 5 toute méthode connue de détermination de la quantité d'insuline embarquée d'un patient. A titre d'exemple, la détermination de la quantité d'insuline embarquée du patient à un instant t peut être calculée par convolution, sur une période allant d'un instant antérieur à l'instant t jusqu'à l'instant t, d'une courbe 10 représentative de l'évolution, en fonction du temps, de la quantité d'insuline injectée au patient avant l'instant t, et d'une fonction fiog représentative de la cinétique de consommation de l'insuline par le corps, par exemple la fonction t-i fros(t) = (1+ x e 15 où t est la variable temps discrétisée et TI est une constante de temps de durée prédéterminée, par exemple comprise entre 40 et 60 minutes, par exemple de l'ordre de 47 minutes. Pour chaque évènement, la valeur de glycémie réelle du patient retenue pour le calcul du ratio Rr est par exemple la 20 valeur de glycémie réelle Gr(t ,-init) à l'instant tinit de début de l'évènement. Ainsi, pour chaque évènement identifié, le ratio Rr est calculé comme suit : Rr = CRr/Gr Pour chaque évènement, la valeur retenue du signal IPA 25 est par exemple la valeur IPA(tinit) à l'instant tinit de début de l'évènement. Pour chaque patient, pour définir une fonction ou un modèle mathématique fCR spécifique au patient, un nombre d'évènements Nbev relativement élevé est d'abord identifié dans 30 les données d'historique du patient, et, pour chaque évènement, une valeur du ratio Rr et une valeur du signal IPA associée sont mesurées. A titre d'exemple, le nombre d'évènements Nbev utilisé pour définir la fonction foR est compris entre 20 et 100, par (tt)- 316895 - DD18588 20 exemple entre 30 et 60, par exemple de l'ordre de 40. La fonction fcR est ensuite déterminée par régression à partir des Nbev points de mesure spécifiques au patient (les points 501 des figures 5A et 5B). A titre d'exemple, pour chaque patient, la fonction fcR 5 est obtenue par régression à partir des Nbev points de mesure 501 obtenus pour le patient. Pour chaque patient, la fonction fcR est alors une fonction d'équation : fcR (IPA) = a x IPA' + c où a, b et c sont des paramètres spécifiques au patient, 10 le paramètre b correspondant à l'ordre du modèle. A titre d'exemple, le paramètre b est fixé égal à 1, le modèle étant alors un modèle linéaire. Sur chacun des diagrammes des figures 5A et 5B, une droite 503 représente la fonction fcR déterminée pour le patient, 15 reliant la variable IPA du patient au ratio estimé R entre son coefficient de sensibilité à l'insuline et sa glycémie. A titre d'exemple, dans le système de la figure 1, la fonction fcR peut être déterminée de façon automatisée par l'unité de traitement et de contrôle 105, lors d'une phase de calibration. 20 L'unité de contrôle et de traitement 105 peut de plus être configurée pour mettre à jour de façon automatisée, par exemple périodiquement, les paramètres de la fonction fcR, pour tenir compte des nouvelles données d'historique enregistrées par le système de régulation au fur et à mesure de son utilisation 25 par le patient. Dans un mode de réalisation préféré, le nombre Nbev d'évènements pris en compte à chaque mise à jour des paramètres de la fonction fcR reste constant. Autrement dit, à chaque fois qu'un nouvel évènement est pris en compte pour la mise à jour des 30 paramètres de la fonction fcR, un évènement antérieur, par exemple l'évènement le plus ancien, est exclu du modèle, ce qui permet que le modèle ne se fige pas et puisse évoluer au cours du temps. 316895 - DD18588 21 La figure 6 est une représentation simplifiée du modèle mathématique fcR utilisé dans le système de la figure 1 pour déterminer le coefficient de sensibilité à l'insuline du patient. Sur la figure 6, le modèle est représenté sous la forme d'un bloc 5 de traitement comportant une entrée e3 sur laquelle est appliquée un signal Gr(t) représentatif de la glycémie réelle du patient, mesurée par le capteur 101, à un instant de mesure t, une entrée e4 sur laquelle est appliqué le signal IPA(t) calculé par l'unité de traitement et de contrôle 105, et une sortie s2 fournissant un 10 signal CR(t) représentatif du coefficient de sensibilité à l'insuline estimé du patient à l'instant t, tel que : CR(t) = Gr (t) x R(t) = Gr (t) X fcR(IPA(t)) La figure 7 est un diagramme illustrant un exemple d'un procédé pouvant être mis en oeuvre par le système de la figure 1 15 pour mettre à jour le modèle mathématique fCR pour tenir compte d'un nouvel historique de données du patient. Le procédé de la figure 7 comprend une étape 701 au cours de laquelle l'unité de traitement et de contrôle acquière et mémorise les données de glycémie réelle Gr mesurée par le 20 capteur 101, les données d'insuline effectivement injectée par le dispositif 103, les données d'ingestion de glucose par le patient, ainsi que les données d'activité physique, sur une période d'observation. La période d'observation considérée à l'étape 701 est choisie pour comprendre au moins un évènement permettant une 25 mesure du coefficient de sensibilité réel CRr du patient. Le procédé de la figure 7 comprend en outre, après l'étape 701, une étape 703 de calcul du ratio Rr = CRr/Gr du patient et de la valeur correspondante du signal IPA du patient (par exemple la valeur du signal IPA au début de l'évènement sur 30 lequel est basé le calcul du ratio Rr=CRr/Gr), à partir des données acquises pendant la période d'observation. Le procédé de la figure 7 comprend de plus, après l'étape 703, une étape 705 de mise à jour du modèle fcR en tenant compte B16895 - DD18588 22 de la nouvelle valeur du ratio Rr et de la valeur correspondante du signal IPA déterminées à l'étape 703. Des modes de réalisation particuliers ont été décrits. Diverses variantes et modifications apparaîtront à l'homme de 5 l'art. En particulier, on a décrit ci-dessus des exemples de réalisation dans lesquels un modèle mathématique fCR est utilisé pour calculer, à partir du signal IPA, une première valeur R représentative du ratio Rr=CRr/Gr du patient, et, par 10 multiplication de cette première valeur avec une valeur de glycémie Gr mesurée par le capteur 101, le coefficient de sensibilité à l'insuline estimé CR du patient (selon la formule susmentionnée CR(t) = Gr (t) X R(t) = Gr (t) X fc R (I P A(t))) . A titre de variante, le modèle mathématique f -CR peut 15 être un modèle deux variables d'entrée permettant de calculer directement le coefficient estimé CR à partir du signal IPA et d'une valeur de glycémie Gr mesurée par le capteur 101, selon la formule suivante : CR(t) =q.cl?(IRAM,GrOW 20 Pour chaque patient, la fonction fCR est alors une fonction d'équation : fcR(IP A, Gr) = a x IP Ab + c x Gra + e où a, b, c, d et e sont des paramètres spécifiques au patient, les paramètres b et d correspondant aux ordres du modèle. 25 A titre d'exemple, les paramètres b et d sont fixés égaux à 1, le modèle étant alors un modèle linéaire. Par ailleurs, les modes de réalisation décrits ne se limitent pas à l'exemple particulier de système de régulation de glycémie décrit en relation avec les figures 1 à 3, à savoir un 30 système à commande prédictive utilisant un modèle mathématique pour prédire l'évolution future de la glycémie du patient et ajuster en conséquence les doses d'insuline à administrer au patient. Plus généralement, la méthode proposée d'estimation en B16895 - DD18588 23 temps réel du coefficient de sensibilité à l'insuline OR du patient, à partir du signal IPA et d'une unique valeur de glycémie mesurée sur le patient, peut être mise en oeuvre dans tout système de régulation de glycémie pouvant tirer profit d'une estimation 5 in-situ et en temps réel du coefficient de sensibilité à l'insuline REVENDICATIONS1. Système automatisé de contrôle de la glycémie d'un patient, comportant un capteur de glycémie (101), un dispositif (107) de mesure d'une activité physique du patient, et une unité de traitement et de contrôle (105), dans lequel : l'unité de traitement et de contrôle (105) est configurée pour générer, par convolution d'un signal PA fourni par le dispositif de mesure d'une activité physique du patient avec une fonction mathématique décroissante H, un signal IPA représentatif de l'influence de l'activité physique du patient 10 sur sa sensibilité à l'insuline ; et l'unité de traitement et de contrôle (105) est configurée pour calculer, à partir d'un premier modèle mathématique fCR spécifique au patient et en tenant compte du signal IPA et d'une unique valeur de glycémie Gr mesurée par le 15 capteur (101), un coefficient CR représentatif de la sensibilité à l'insuline du patient. 2. Système selon la revendication 1, comportant en outre un dispositif d'injection d'insuline (103), dans lequel l'unité de traitement et de contrôle (105) est configurée pour 20 commander le dispositif d'injection d'insuline (103) en tenant compte du coefficient CR. 3. Système selon la revendication 2, dans lequel l'unité de traitement et de contrôle (105) est configurée pour prédire, à partir d'un deuxième modèle mathématique (201), 25 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. 4. Système selon l'une quelconque des revendications 1 à 3, dans le premier modèle mathématique est une fonction 30 d'équation CR = fcR(IPA,Gr) = a x IPAb + c x Grd + e B16895 - DD18588 25 où a, b, c, d et e sont des paramètres spécifiques au patient. 5. Système selon l'une quelconque des revendications 1 à 3, dans lequel le premier modèle mathématique est une fonction 5 d'équation CR = Gr X fcR(IPA) = Gr X ( a x IPAb + c) où a, b et c sont des paramètres spécifiques au patient. 6. Système selon l'une quelconque des revendications 1 à 5, dans lequel l'unité de traitement et de contrôle (105) est 10 configurée pour mettre en oeuvre une étape de calibration automatique du premier modèle foR en tenant compte d'un historique de glycémie mesurée par le capteur (101), d'un historique d'insuline injectée au patient, d'un historique d'ingestion de glucose par le patient, et d'un historique du signal d'activité 15 physique PA du patient, sur une période d'observation passée. 7. Système selon la revendication 6, dans lequel l'unité de traitement et de contrôle (105) est configurée pour, lors de l'étape de calibration automatique, mesurer une pluralité de valeurs du coefficient de sensibilité à l'insuline réel CRr du 20 patient lors d'une pluralité d'évènements de mesure compris dans la période d'observation passée. 8. Système selon la revendication 7, dans lequel chaque évènement de mesure correspond à une plage temporelle continue allant d'un instant initial tinit à un instant final tfinal, 25 répondant aux critères suivants : - l'instant tinit se situe dans une phase d'hyperglycémie, c'est-à-dire une phase dans laquelle la glycémie du patient est supérieure à un seuil prédéterminé ; 30 - un bolus de correction, c'est-à-dire une dose d'insuline a été administré au patient après le début de la phase d'hyperglycémie et avant l'instant tinit, en vue de limiter la durée de la phase d'hyperglycémie ; B16895 - DD18588 26 - la glycémie du patient décroit de façon continue entre l'instant initial tinit et l'instant final tfinal ; et - aucune ingestion de glucose par le patient n'a été 5 effectuée entre l'instant tinit - Tj et l'instant tfinalf où Tj est une durée de jeun prédéterminée. 9. Système selon la revendication 7 ou 8, dans lequel l'unité de traitement et de contrôle (105) est configurée pour, lors de l'étape de calibration automatique, déterminer le premier modèle mathématique fcR par régression à partir de ladite pluralité de valeurs du coefficient de sensibilité à l'insuline réel CRr. 10. Système selon l'une quelconque des revendications 1 à 9, dans lequel la fonction H est une fonction exponentielle 15 décroissante. 11. Système selon l'une quelconque des revendications 1 à 10, dans lequel le dispositif de mesure (107) comprend une interface utilisateur par laquelle le patient déclare les activités physiques qu'il réalise. 20 12. Système selon l'une quelconque des revendications 1 à 11, dans laquelle le dispositif de mesure (107) comprend un ou plusieurs capteurs adaptés à mesurer des grandeurs représentatives de l'activité physique du patient. 13. Système selon la revendication 12, dans laquelle le 25 dispositif de mesure (107) comprend un capteur de mouvement et/ou un capteur de rythme cardiaque.