FR3138999A1 - Procédé d’estimation d’une concentration de gaz dégagé par un milieu - Google Patents
Procédé d’estimation d’une concentration de gaz dégagé par un milieu Download PDFInfo
- Publication number
- FR3138999A1 FR3138999A1 FR2208595A FR2208595A FR3138999A1 FR 3138999 A1 FR3138999 A1 FR 3138999A1 FR 2208595 A FR2208595 A FR 2208595A FR 2208595 A FR2208595 A FR 2208595A FR 3138999 A1 FR3138999 A1 FR 3138999A1
- Authority
- FR
- France
- Prior art keywords
- gas
- interest
- concentration
- medium
- collection chamber
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 239000012080 ambient air Substances 0.000 claims abstract description 14
- 239000008280 blood Substances 0.000 claims description 69
- 210000004369 blood Anatomy 0.000 claims description 69
- 238000009792 diffusion process Methods 0.000 claims description 36
- 239000003570 air Substances 0.000 claims description 35
- 238000005259 measurement Methods 0.000 claims description 34
- 230000008569 process Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 239000007788 liquid Substances 0.000 claims description 4
- 210000004204 blood vessel Anatomy 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000001939 inductive effect Effects 0.000 claims description 2
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 103
- 239000007789 gas Substances 0.000 description 76
- 229910002092 carbon dioxide Inorganic materials 0.000 description 55
- 239000002609 medium Substances 0.000 description 32
- 239000013598 vector Substances 0.000 description 31
- 239000011159 matrix material Substances 0.000 description 30
- 239000001569 carbon dioxide Substances 0.000 description 24
- 230000014509 gene expression Effects 0.000 description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 14
- 238000005070 sampling Methods 0.000 description 13
- 241000897276 Termes Species 0.000 description 8
- 230000003190 augmentative effect Effects 0.000 description 7
- 230000017531 blood circulation Effects 0.000 description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 230000005855 radiation Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- BVKZGUZCCUSVTD-UHFFFAOYSA-M Bicarbonate Chemical compound OC([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-M 0.000 description 4
- 239000010931 gold Substances 0.000 description 4
- 239000001257 hydrogen Substances 0.000 description 4
- 229910052739 hydrogen Inorganic materials 0.000 description 4
- -1 hydrogen ions Chemical class 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000011144 upstream manufacturing Methods 0.000 description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 description 3
- 230000004087 circulation Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 239000012528 membrane Substances 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 241000894007 species Species 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- CSCPPACGZOOCGX-UHFFFAOYSA-N Acetone Chemical compound CC(C)=O CSCPPACGZOOCGX-UHFFFAOYSA-N 0.000 description 2
- 208000010444 Acidosis Diseases 0.000 description 2
- 208000005223 Alkalosis Diseases 0.000 description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 206010020952 Hypocapnia Diseases 0.000 description 2
- RRHGJUQNOFWUDK-UHFFFAOYSA-N Isoprene Chemical compound CC(=C)C=C RRHGJUQNOFWUDK-UHFFFAOYSA-N 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- MWUXSHHQAYIFBG-UHFFFAOYSA-N Nitric oxide Chemical compound O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 2
- 230000007950 acidosis Effects 0.000 description 2
- 208000026545 acidosis disease Diseases 0.000 description 2
- 230000002340 alkalosis Effects 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 244000005700 microbiome Species 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 239000002002 slurry Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 208000025721 COVID-19 Diseases 0.000 description 1
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 208000006545 Chronic Obstructive Pulmonary Disease Diseases 0.000 description 1
- 235000008733 Citrus aurantifolia Nutrition 0.000 description 1
- 208000035473 Communicable disease Diseases 0.000 description 1
- 206010020591 Hypercapnia Diseases 0.000 description 1
- 235000011941 Tilia x europaea Nutrition 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 239000012491 analyte Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010241 blood sampling Methods 0.000 description 1
- 229910002091 carbon monoxide Inorganic materials 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 235000021183 entrée Nutrition 0.000 description 1
- 235000019441 ethanol Nutrition 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 239000001963 growth medium Substances 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 230000002209 hydrophobic effect Effects 0.000 description 1
- 208000000122 hyperventilation Diseases 0.000 description 1
- 230000000870 hyperventilation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 239000004571 lime Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000037081 physical activity Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007430 reference method Methods 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 208000023504 respiratory system disease Diseases 0.000 description 1
- 230000005070 ripening Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 150000003384 small molecules Chemical class 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14542—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/083—Measuring rate of metabolism by using breath test, e.g. measuring rate of oxygen consumption
- A61B5/0836—Measuring rate of CO2 production
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
- A61B5/14551—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1468—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using chemical or electrochemical methods, e.g. by polarographic means
- A61B5/1477—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using chemical or electrochemical methods, e.g. by polarographic means non-invasive
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/483—Physical analysis of biological material
- G01N33/487—Physical analysis of biological material of liquid biological material
- G01N33/49—Blood
- G01N33/4925—Blood measuring blood gas content, e.g. O2, CO2, HCO3
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Optics & Photonics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
Procédé d’estimation d’une teneur d’un gaz d’intérêt dans un milieu, à l’aide d’un dispositif de mesure, destiné à être disposé contre le milieu, le dispositif s'étendant entre une face de contact, destinée à être appliquée contre le milieu et une extrémité distale, le dispositif comportant une paroi latérale, s'étendant entre la face de contact et l'extrémité distale, le dispositif comportant : au niveau de la face de contact, au moins une ouverture d’admission, configurée pour collecter le gaz d'intérêt émis par le milieu, l'ouverture d’admission étant pratiquée à travers la face de contact; une chambre de mesure, comportant un capteur de gaz , le capteur de gaz étant configuré pour mesurer une concentration de gaz d’intérêt s’écoulant à travers la chambre de mesure;une chambre de collecte (30), reliée à la chambre de mesure, et délimitée par la paroi latérale, la chambre de collecte comportant au moins une ouverture latérale (34), ménagée à travers la face latérale, de façon à admettre de l’air ambiant dans la chambre de collecte
Description
Le domaine technique de l'invention est la mesure d'un gaz dégagé par un milieu à l’aide d’un dispositif compact. Le dispositif compact est appliqué contre le milieu. Le milieu peut être la peau d'un être vivant. Le procédé permet alors d’estimer une teneur d’un gaz dégagé par le milieu, Le gaz peut notamment être du dioxyde de carbone, pour des applications de capnométrie. Il s’agit alors d’estimer la teneur de dioxyde de carbone dissous dans le sang. Le milieu peut être un milieu liquide, par exemple de l’eau, ou du lisier. Le gaz peut être du dioxyde de carbone ou du méthane. Le milieu peut également être un milieu végétal, auquel cas le dispositif peut être utilisé pour étudier la respiration du milieu.
Certaines maladies respiratoires affectent les échanges gazeux entre le sang et l'air expiré. Le sang contient des gaz dissous, entre autres l'oxygène et le dioxyde de carbone (CO2), dont les pressions partielles respectives reflètent les échanges gazeux se produisant au niveau des poumons et des organes.
Pour évaluer la concentration de CO2dissous dans le sang, on peut avoir recours à un prélèvement sanguin. Il s'agit alors d'une méthode invasive, pouvant être douloureuse et délicate à appliquer, en particulier en néonatologie. De plus, elle ne peut être appliquée que de façon ponctuelle. En dépit de ces inconvénients, sa fiabilité est validée par le corps médical et elle constitue une méthode de référence. Une autre méthode consiste à estimer la teneur du sang en CO2de façon non invasive, en effectuant une mesure de la pression partielle de CO2diffusant à travers les tissus, et notamment la peau. Cette méthode est désignée par le terme de capnométrie transcutanée.
Dans le sang, le gaz carbonique est dissous suivant une espèce moléculaire (CO2), et suivant deux espèces ioniques : les ions carbonates CO3 2-et bicarbonates HCO3 -. L’équilibre entre ces espèces dépend du pH du sang, car les ions carbonates et bicarbonates sont en équilibre avec les ions hydrogène H+. Ainsi, la concentration des ions carbonates ou bicarbonates influence le pH du sang. Une augmentation de la concentration de CO2dissous (hypercapnie), se traduit par une augmentation de la quantité d’ions carbonates et bicarbonates, et, par équilibre, des ions hydrogène, ce qui engendre une diminution du pH du sang, ou acidose. Une augmentation de la concentration de CO2peut survenir lorsque l’élimination du CO2par les voies respiratoires est insuffisante, par exemple dans le cas de maladies pulmonaires obstructives chroniques, (COPD) ou de maladies infectieuses affectant les poumons, un exemple étant une infection due au COVID 19.
Inversement, une diminution de la concentration de CO2dissous (hypocapnie) abaisse la concentration d’ions hydrogène, ce qui entraîne une augmentation du pH, ou alcalose. Une hypocapnie peut être causée par exemple par une hyperventilation associée à une augmentation de la fréquence respiratoire. La survenue d’une acidose ou d’une alcalose peut avoir des conséquences sur le métabolisme. Ainsi, la concentration de CO2dans le sang est un paramètres vital important, qu’il convient de surveiller régulièrement pour certains patients à risque.
Le suivi de la concentration de CO2peut également concerner des patients en réanimation, ou chez des nouveaux-nés placés en incubateur. Cela peut également trouver des applications dans le suivi d’efforts, l’activité physique favorisant la production de gaz carbonique.
Ce type d'analyse basé sur des mesures transcutanées a été introduit dans les années 1980. Les mesures transcutanées permet d'effectuer un suivi continu, par exemple pour suivre les effets immédiats d'une prise en charge thérapeutique influençant la concentration du CO2.dans le sang. Elle peut également permettre de déterminer les instants auxquels une quantification plus précise, par prélèvement sanguin, est préférable. On comprend ainsi que les méthodes invasives et non invasives peuvent être complémentaires : l'une est précise est ponctuelle, tandis que l'autre peut être mise en œuvre en œuvre de façon continue pour un suivi longitudinal.
Un dispositif compact permettant la mesure du CO2transcutané a déjà été décrit dans WO2020/249466. Il s’agit d’un dispositif de mesure non invasif, porté par un utilisateur, pour estimer une concentration d'un gaz d'intérêt émis de façon transcutanée, le gaz d'intérêt pouvant par exemple être le dioxyde de carbone. La circulation du gaz dans le dispositif permet une collecte du gaz d'intérêt. Le gaz d’intérêt se propage à travers le dispositif par convection, sous l’effet d’une source de chaleur. Le dispositif comporte une face de contact, destinée à être appliquée contre l’utilisateur. Le dispositif comporte également une chambre de mesure, comportant un capteur du gaz d’intérêt, typiquement le CO2. Pour que l’augmentation de température se traduise par une convection, il est préférable que le la chambre de mesure soit disposée en dessus de la face de contact.
Les inventeurs proposent un perfectionnement du dispositif décrit dans WO2020/249466, visant à améliorer certaines performances, en particulier les performances de réponse temporelle. Les inventeurs proposent également un procédé permettant d’estimer, à l’aide du dispositif, la concentration d’un analyte dans le milieu face auquel le dispositif est disposé.
Un objet de l’invention est un procédé d’estimation d’une teneur d’un gaz d’intérêt dans un milieu, à l’aide d’un dispositif de mesure, destiné à être disposé contre le milieu, le dispositif de mesure s'étendant entre une face de contact, destinée à être appliquée contre le milieu et une extrémité distale, le dispositif de mesure comportant une paroi latérale, s'étendant entre la face de contact et l'extrémité distale, le dispositif de mesure comportant :
- au niveau de la face de contact, au moins une ouverture d’admission, configurée pour collecter le gaz d'intérêt émis par le milieu, l'ouverture d’admission étant pratiquée à travers la face de contact;
- une chambre de mesure, comportant un capteur de gaz, le capteur de gaz étant configuré pour mesurer une concentration de gaz d’intérêt s’écoulant à travers la chambre de mesure;
- une chambre de collecte, reliée à la chambre de mesure, et délimitée par une ouverture sur la paroi latérale, la chambre de collecte comportant au moins une ouverture latérale, ménagée à travers la face latérale, ou sur la paroi supérieure de la chambre de collecte de façon à admettre de l’air ambiant dans la chambre de collecte ;
le dispositif étant tel que :
- la chambre de mesure est disposée entre la face de contact et la chambre de collecte;
- le dispositif comporte un moyen d’entraînement, configuré pour entraîner l’air de la chambre de collecte vers une ouverture d’évacuation, l’entraînement d’air induisant un transport du gaz d’intérêt de la face de contact vers la chambre de collecte, à travers la chambre de mesure :
- le procédé comportant :
- a) mesure de la concentration de gaz d’intérêt dans la chambre de mesure ;
- b) à l’aide d’une unité de traitement, modélisation du transport du gaz d’intérêt, entre le milieu et la chambre de collecte, la modélisation comportant une prise en compte de la diffusion du gaz depuis le milieu jusqu’à la chambre de collecte à travers la face de contact et le dispositif, ainsi qu’une convection du gaz dans la cellule de collecte résultant de l’entraînement produit par le moyen d’entraînement ;
- c) à partir de la mesure résultant de l’étape a), et d’une prise en compte du modèle résultant de l’étape b), estimation de la teneur du gaz d’intérêt dans le milieu.
Selon un mode de réalisation,
- le dispositif est spatialement discrétisé selon un maillage spatial, définissant des points de maillage entre la face de contact et la chambre de collecte ;
- le modèle est un modèle spatio-temporel discrétisé, de façon à estimer une teneur en gaz d’intérêt en différents points du maillage, et en différents instants.
Le procédé peut être tel que l’étape b) comporte :
- modélisation de la diffusion du gaz d’intérêt à travers la chambre de mesure ;
- modélisation de la diffusion et de la convection du gaz d’intérêt dans la chambre de collecte.
L’étape b) comporte une modélisation du transport du gaz d’intérêt dans le milieu.
De préférence, l’étape c) met en œuvre un estimateur récursif. Il peut s’agir d’un estimateur récursif linéaire. L’estimateur récursif linéaire peut être un filtre de Kalman.
Selon une application, le milieu est la peau d’un utilisateur, la peau s’étendant en dessus d’un vaisseau sanguin. L’étape b) peut alors comporter :
- à partir de la concentration de gaz d'intérêt dans la chambre de mesure, résultant de l'étape a), estimation d'une concentration de gaz d'intérêt transcutané ;
- à partir de la concentration de gaz d'intérêt transcutané résultant de la sous-étape (i), estimation d'une concentration ou pression partielle de gaz d'intérêt dissous dans le sang de l'utilisateur.
Le procédé peut comporter une modélisation du transport du gaz d’intérêt émis par le milieu à travers le dispositif.
Le milieu à analyser peut être un milieu solide ou liquide.
Le moyen d’entraînement peut être un propulseur d’air ou un aspirateur d’air.
L'invention sera mieux comprise à la lecture de l'exposé des exemples de réalisation présentés, dans la suite de la description, en lien avec les figures listées ci-dessous.
La montre un premier mode de réalisation d’un dispositif selon l’invention.
La schématise les écoulements à l’intérieur du dispositif.
La schématise les principaux composants d’un capteur de gaz, disposé dans la chambre de mesure du dispositif. La vue est prise dans un plan perpendiculaire à un axe transversal.
La schématise les principaux composants d’un capteur de gaz, basé sur une absorption d’une lumière infra-rouge, et disposé dans la chambre de mesure du dispositif. La vue est prise dans un plan parallèle à l’axe transversal.
La schématise une décomposition du dispositif de mesure et du milieu à analyser en compartiments.
La schématise une discrétisation de chaque compartiment représenté sur la .
La illustre les principales étapes de mise en œuvre d’un procédé selon l’invention.
Les figures 6A, 6B et 6C montrent une estimation d’une concentration de CO2 en différents points d’échantillonnage lorsque la concentration dans le sang suit une variation en créneau.
La montre une estimation de la concentration de CO2 mesurée par le capteur de gaz du dispositif, ainsi qu’une application d’un bruit blanc sur ladite concentration. L’application du bruit blanc
Les figures 7A, 7B, 7C, 7D montrent une estimation d’une concentration de CO2 en différents points d’échantillonnage lorsque la concentration mesurée par le capteur suit la concentration bruitée simulée sur la .
La est une comparaison entre la concentration dans le sang en créneau, qui représente la réalité, et la concentration dans le sang estimée par l’estimateur récursif.
Les figures 1A et 1B sont des vues générales d'un exemple de dispositif 1 selon l'invention. Le dispositif 1 est destiné à être disposé au contact d’un milieu que l'on souhaite analyser. Dans l’exemple décrit, le milieu est de la peau S d'un utilisateur, humain ou animal. Le dispositif comporte un corps principal 2 ainsi qu'un élément de fixation 3, ce dernier étant, dans cet exemple, un bracelet. De façon alternative, le corps principal peut être disposé au contact du lobe d’une oreille, ou d’un doigt. Le support peut être une pince ou être intégré dans un casque auditif. De façon plus générale, le support est configuré pour maintenir le corps principal au contact du milieu à analyser.
Le dispositif est destiné à estimer une concentration d’un gaz d’intérêt émanant de la peau d’un utilisateur. Par gaz d'intérêt, il est entendu un gaz dont on souhaite déterminer une concentration dans un corps humain ou animal vivant, et plus particulièrement dans le sang. Dans l'exemple décrit ci-dessous, de façon non limitative, le gaz d'intérêt est le dioxyde de carbone, dont on cherche à estimer une teneur dans le sang de l'utilisateur. Selon d’autres possibilités, le gaz d’intérêt peut être, de façon non limitative, l’oxygène, l’alcool éthylique, le monoxyde de carbone, le méthane, le monoxyde d’azote, l’acétone ou l’isoprène, l’hydrogène, certaines drogues ou substances volatiles.
Le corps principal 2 comporte une face de contact 10, destinée à être apposée sur la peau S. La face de contact 10 est sensiblement planaire, au sens où elle s’étend parallèlement à un plan XY, certaines portions pouvant être inclinées par rapport au plan XY. Le corps principal comporte également une extrémité distale 4, opposée à la face de contact 10. La surface de contact 10 et l’extrémité distale 4 sont reliées l’une à l’autre par une face latérale 5, s’étendant autour d’un axe central Δ, parallèle à un axe transversal Z, ce dernier étant perpendiculaire au plan XY.
La face de contact 10 peut être formée par une membrane perméable au CO2ou une plaque comportant des ouvertures d’admission 12 permettant un passage du CO2. La face de contact peut comporter une membrane hydrophobe, de façon à éviter la diffusion de la vapeur d’eau à travers le dispositif 1. La face de contact peut être chauffée, ce qui permet d’augmenter le débit sanguin. Cela améliore la collecte.
Selon une possibilité, le dispositif comporte également un élément de chauffage, permettant de porter la face de contact 10, délimitant la chambre de mesure 20, à une température supérieure à 37°C, et de préférence comprise entre 40°C et 50°C, et de préférence entre 40°C et 45 °C, par exemple 42°C. L’élément de chauffage est par exemple une résistance ménagée sur la face de contact, produisant un chauffage par effet Joule. Une augmentation locale et modérée de la température, au voisinage de la peau, favorise en effet une augmentation du débit sanguin par dilatation des capillaires sanguins, ce qui augmente la diffusion d’un gaz d’intérêt transcutané, à travers la peau.
La face de contact 10 débouche sur une chambre de mesure 20. La fonction de la chambre de mesure est d’estimer une concentration de CO2du mélange gazeux circulant dans le corps principal 2, parallèlement, ou sensiblement parallèlement, à l’axe central Δ. A cette fin, la chambre de mesure 20 comporte un capteur de gaz 23. Plusieurs types de capteurs peuvent être utilisés à cette fin, par exemple des capteurs optiques ou des capteurs électrochimiques, ces derniers pouvant notamment être basés sur des oxydes métalliques (capteurs MOX), ou des capteurs photoacoustiques. Les inventeurs ont estimé qu’il était préférable d’utiliser un capteur optique, et plus précisément un capteur infra-rouge. Un tel capteur ne nécessite pas de maintenance particulière, et est particulièrement compact, ainsi que peu onéreux. De plus, un tel capteur est très spécifique pour caractériser des liaisons chimiques. Il convient à la détection de molécules de petites tailles, par exemple le dioxyde de carbone.
Le capteur de gaz est un capteur de type NDIR (Non Dispersive Infra Red). Ce type de capteur comporte une source de rayonnement infrarouge 24, émettant généralement dans une bande spectrale comprise entre 1 µm et 20 µm. Il comporte également au moins un photodétecteur de mesure 25, sensible au rayonnement infra-rouge. Le principe repose sur l’atténuation, par le gaz analysé, du rayonnement infra-rouge, émis par la source. La source infra-rouge 24 et le photodétecteur de mesure 25 forment le capteur de gaz 23. Le photodétecteur de mesure 25 est par exemple une thermopile. Un filtre 26 placé devant la thermopile détermine la longueur d’onde du rayonnement infrarouge mesuré par la thermopile. Différentes configurations du capteur de gaz sont décrites par la suite.
Le capteur de gaz 23 est configuré de façon à ce que le gaz d’intérêt transcutané, en l’occurrence le CO2, se propage entre la source de lumière et le photodétecteur, parallèlement à l’axe central Δ, ou sensiblement parallèlement à l’axe central Δ. Par sensiblement parallèlement, on entend parallèlement en admettant une tolérance angulaire inférieure à 30° ou 20° par rapport à la parallèle.
La source infra-rouge émet une lumière se propageant perpendiculairement à l’axe transversal Z. De préférence, la lumière émise par la source est un faisceau en forme de nappe lumineuse perpendiculaire à l’axe transversal Z. Selon l’axe transversal Z, l’épaisseur du faisceau est de préférence inférieure à 1 cm, et de préférence inférieure à 5 mm. De préférence, la nappe lumineuse s’étend sur au moins 50 %, voire au moins 80% de la section transversale de la face de contact 10, cette dernière étant comprise entre quelques cm² et 25 ou 30 cm². Par section transversale, on entend une surface perpendiculaire à l’axe transversal Z (ou à l’axe central Δ). La faible épaisseur du faisceau permet de diminuer le temps de réponse. La surface élevée de la section transversale du faisceau permet d’augmenter la quantité de gaz prélevée et donc la sensibilité de la mesure.
Afin de favoriser un écoulement selon l’axe transversal Z, à travers la chambre de mesure 20 la face de contact comporte une multitudes d’ouvertures d’admission 12, réparties selon une section transversale de quelques cm², par exemple entre 5 et 25 cm². Alternativement, la face de contact est une membrane poreuse, perméable au CO2. La chambre de mesure 20 débouche sur une chambre de collecte 30. L’interface entre la chambre de mesure 20 et la chambre de collecte 30 peut être formée par une ouverture de collecte large, ou par une pluralité d’ouvertures de collecte 22, ménagées dans une plaque de collecte 21 et réparties selon une section transversale de surface élevée, de façon similaire aux ouvertures d’admission 12. Chaque ouverture de collecte 22 débouche dans une chambre de collecte 30.
La chambre de collecte 30 comporte des ouvertures latérales 34, ménagées dans la face latérale 5 du dispositif 1. Les ouvertures latérales sont configurées pour permettre une admission d’air ambiant à l’intérieur du dispositif 1. L’air ambiant est l’air situé à l’extérieur du dispositif. Le flux d’air ambiant s’écoulant à travers les ouvertures latérales 34 est dirigé vers l’axe central Δ. Chaque ouverture latérale 34 est de préférence orientée perpendiculairement à l’axe central Δ, de façon à ce que l’air ambiant pénètre dans le dispositif 1 , à travers chaque ouverture latérale 34, selon une direction perpendiculaire, ou sensiblement perpendiculaire, à l’axe central Δ.
L’admission de l’air ambiant est obtenue par une pompe 41, agencée pour appliquer une dépression dans la chambre de collecte 30 par rapport à la pression ambiante. La pompe 41 peut être située dans un conduit d’évacuation 40, adjacent à la chambre de collecte 30. Outre l’admission d’air ambiant dans la chambre de collecte, la pompe permet l’écoulement du gaz d’intérêt à travers la chambre de mesure 10 jusqu’à la chambre de collecte 20, à travers les ouvertures de collecte 22. De préférence, le débit est ajusté de façon à ce que l’écoulement du gaz d’intérêt à travers la chambre de mesure 20 soit laminaire. Le débit d’air dans la chambre de collecte peut valoir par exemple entre 0,1 et 10 mL/min, suivant le volume de la chambre.
Dans la chambre de collecte 30, l’air ambiant, admis à travers les ouvertures latérales 34, se mélange avec le gaz d’intérêt, admis à travers les ouvertures de collecte 22. Le conduit d’évacuation 40 est également mis en dépression par la pompe 41, de façon à ce que l’air, mélangé au gaz d’intérêt, est évacué dans le conduit d’évacuation, jusqu’à une ouverture d’évacuation 42, formant l’extrémité distale 4 du dispositif. Dans l’exemple représenté sur les figures 1A et 1B, l’air se propage à travers l’ouverture d’évacuation 42 autour de l’axe central Δ. Selon une autre possibilité, l’air se propage à travers l’ouverture d’évacuation 42 selon un axe d’évacuation sécant, et notamment perpendiculaire, à l’axe central Δ.
Dans l’exemple représenté, les ouvertures latérales sont disposées sur deux faces opposées 51, 52de la paroi latérale. Cela favorise une symétrie de l’écoulement d’air dans le dispositif. Toutefois, un tel agencement symétrique des ouvertures latérales n’est pas nécessaire. Les ouvertures latérales 34 peuvent n’être ménagées que sur une face de la paroi latérale, le conduit d’évacuation étant alors disposé sur une face opposée de la paroi latérale.
Chaque ouverture latérale peut être associée à un filtre à CO2, par exemple un filtre comportant de la chaux, de façon à piéger le CO2présent dans l’air ambiant.
Le dispositif comporte une unité de traitement 50. L’unité de traitement 50 comporte des moyens de calcul, par exemple un microprocesseur ou un microcontrôleur, embarqué sur le dispositif, en étant relié par une liaison filaire ou sans fil à un téléphone portable ou à un calculateur de type PC. L'unité de traitement 50 est également configurée pour mettre en œuvre le procédé d'estimation de la teneur en CO2dissous du sang de l'utilisateur, à partir de la concentration de CO2résultant du capteur de gaz 23.
Les principes généraux régissant l’estimation de la teneur en CO2dissous dans le sang à partir de la concentration de CO2détectée dans le capteur de gaz ont été décrits dans WO2020/249466.
Le procédé comporte un calcul de la concentration du gaz carbonique dans la chambre de mesure à partir des mesures réalisées par les deux thermopiles à partir de modèles décrivant l’atténuation du rayonnement infrarouge aux deux longueurs d’onde associées à chacune des thermopiles (modèle de Beer Lambert, modèle linéaire quadratique). Ces modèles sont décrits dans WO2020/249466.
La schématise les courants fluidiques formés à l’intérieur du dispositif 1. Le CO2transcutané est admis dans le dispositif par les ouvertures d’admission 12 pratiquées dans la face de contact 10, et débouchant dans la chambre de mesure 20 (cf. flèche F1). Le CO2se propage à travers la chambre de mesure par diffusion, parallèlement à l’axe central, ce dernier étant parallèle l’axe transversal Z. (cf. flèches F2). L'air ambiant est admis à travers les ouvertures latérales 34 pratiquées dans la paroi latérale (cf. flèches F3). Le mélange gazeux comportant l'air ambiant et le CO2transcutané est réalisé dans la chambre de collecte 30. Ce mélange se propage vers l’ouverture d’évacuation 42 (cf. flèches F4) pour sortir du corps 2 du dispositif (cf. flèches F5)
Un aspect important du dispositif est le recours à la pompe 41, qui permet de former un courant de convection d’air au niveau de la chambre de collecte 30, pour évacuer le mélange air / gaz d’intérêt transcutané. Cela permet également d’accélérer la circulation d’air dans la chambre de collecte. Il en résulte un effet de pompage accélérant le transport du CO2du sang vers l’air, à travers la peau et la chambre de mesure. Il en résulte une diminution du temps de réponse du dispositif.
La chambre de mesure 20 comporte, de préférence, un capteur de température ainsi qu'un capteur d'humidité et un capteur de pression.
Le capteur de température sert à convertir la concentration mesurée en pression suivant la loi des gaz parfaits. A l’aide de la loi des gaz parfaits, nous en déduisons la pression partielle du gaz carbonique dans la chambre de mesure à la température :
Où est exprimée en (mol/m3),Rest la constante des gaz parfaits (0.0623637 m3.mmHg.K-1.mol-1), en K et en mmHg.
Le capteur de pression atmosphérique sert à calculer la pression relative exprimée en ppm (partie par million):
Où et sont exprimées dans la même unité, ici en mmHg.
Le capteur d’humidité sert à calculer la concentration de la vapeur d’eau dans la chambre de mesure. Un facteur d’atténuation du rayonnement important concerne la vapeur d’eau. La prise en compte de l’humidité nécessite de convertir la valeur d’humidité relative RHmesfournie par le capteur d’humidité en concentration de vapeur d’eau .L’humidité relative correspond à la pression partielle d’eau dans la chambre de mesure, notée , sur la pression de vapeur saturante à la température de la chambre de mesure.
La concentration en H2O dans la chambre de mesure peut être obtenue selon l’expression :
En combinant les deux équations précédentes, on obtient :
Où représente les valeurs de pression de vapeur d’eau saturante En fonction de la température ces valeurs peuvent être obtenues à partir de tables, où à partir d’une expression analytique, par exemple l’équation de Tetens
La concentration peut être utilisée pour estimer la concentration de CO2dans la chambre de mesure, en utilisant l’expression de Beer-Lambert
Où et sont les concentrations molaires, , , , les coefficients d’atténuation du gaz carbonique et de la vapeur d’eau, , les atténuations de l’air aux longueurs d’onde , .
Ainsi, la prise en compte de la concentration d’H2O permet d’améliorer la précision avec laquelle la concentration de CO2est mesurée, comme évoqué dans la demande de brevet WO2020/249466.
La montre la chambre de mesure dans un plan PXZparallèle à l’axe transversal Z. La chambre de mesure peut s’étendre entre deux grilles 21 et 21’ perméables au gaz d’intérêt. La grille 21 forme la paroi d’interface avec la chambre de collecte 30. Elle comporte des ouvertures de collecte, permettant l’écoulement du gaz d’intérêt de la chambre de mesure vers la chambre de collecte. De préférence, chaque grille est réfléchissante dans la bande spectrale d’émission de la source infra-rouge. La chambre de mesure peut comporter une grille réfléchissante 21’, similaire à la grille 21, au niveau de la face de contact 10. La présence de chaque grille réfléchissante permet d’augmenter la quantité de lumière se propageant à travers le gaz d’intérêt, vers chaque voie du photodétecteur. Elle permet également de limiter l’épaisseur de la nappe de lumière, selon l’axe transversal.
La distance d entre la grille 21 et la face de contact, ou entre les grilles 21 et 21’, est par exemple comprise entre 3 mm et 8 mm.
Les mesures résultant du capteur peuvent être bruitées, auquel cas un filtrage passe bas peut être appliqué. Le filtrage passe-bas peut être réalisé par exemple par un filtre à moyenne glissante dont la longueur est adaptée au niveau du bruit de mesure.
A partir de , la pression partielle de CO2 est déduite par :
A partir du capteur de pression atmosphérique de la chambre de mesure, la pression relative est calculée par :
Où est la pression atmosphérique dans la chambre de mesure.
A partir de la concentration de CO2 résultant du capteur de mesure, la concentration en gaz carbonique dans le sang peut être estimée en appliquant un modèle du transport de CO2 entre le sang et la chambre de mesure. On décrit par la suite un modèle de transport monodirectionnel, modélisant le transport de CO2, par diffusion et par convection.
La figure 3 schématise un exemple de dispositif apposé sur la peau S de l’utilisateur. Le sang occupe un compartiment sang B, entre les coordonnées et , et son épaisseur est . Le compartiment sang B est recouvert par de la peau S. La peau S s’étend entre les coordonnées et , et son épaisseur est . La peau est recouverte par le dispositif. Dans le dispositif de l’art antérieur, le dispositif comporte une superposition de la chambre de collecte et de la chambre de mesure. La chambre de collecte s’étend entre les coordonnées et , et son épaisseur est La chambre de mesure s’étend entre les coordonnées et , et son épaisseur est ; . Selon l’invention, le dispositif comporte une superposition de la chambre de mesure 20, de la chambre de collecte 30.
La diffusion est caractérisée par un coefficient de diffusion, qui dépend du milieu traversé, en l’occurrence le sang, la peau et l’air.
Dans la modélisation 1D, le transport du gaz carbonique est décrit uniquement selon l’axe central correspondant à l’axe Z . Le sang s’écoule, dans le compartiment sang, selon un courant de convection de vitesse moyenne Dans la peau, le CO2 migre uniquement par diffusion. Dans le dispositif, dans la cellule de mesure le CO2 migre uniquement par diffusion et dans la cellule de collecte le CO2 migre par diffusion et par convection. La convection résulte de la circulation d’air induite par la pompe dans la chambre de collecte. La convection d’air dans al chambre de collecte entraîne une diffusion à travers la face de contact et à travers la chambre de mesure. Elle est caractérisée par une vitesse moyenne selon l’axe Z.
En se basant sur une propagation du gaz d’intérêt selon une direction Z perpendiculaire à la surface de la peau (modèle monodimensionnel), on peut écrire dans tout milieu homogène ( _i-e de même constante de Henry) :
où :
En considérant un milieu homogène :
L’équation (10) devient
Dans l’expression (11), le terme correspond à la dynamique due àla convection. Le terme correspond à la dynamique due àde la diffusion. L’expression (11) correspond à une équation de convection – diffusion.
L’expression (11) modélise le transport du gaz d’intérêt selon l’axe Z à l’intérieur de chaque compartiment précédemment décrit (chambre de mesure, chambre de collecte).
Le modèle prend également en compte une équation d’observation :
Les conditions aux limites sont :
où correspond à la coordonnée, selon Z, de l’entrée du compartiment sang.
et
où correspond à la coordonnée, selon Z, de l’extrémité distale du dispositif.
On décrit à présent les conditions aux interfaces. Dans les équations qui suivent, les indices i et j désignent deux milieux successifs respectivement en amont et en aval de l’interface, l’amont et l’aval étant pris en compte dans le sens de propagation nette du CO2, c’est-à-dire orienté du compartiment sang vers la chambre de collecte.
Considérant que représente la coordonnée spatiale selon l’axe Z de l’interface entre le milieu i et le milieu j, est le flux gazeux sortant du milieu i et est le flux gazeux entrant dans le milieu j, par unité de surface. Pour une interface transparent au gaz nous avons :
avec :
et
Dans les expressions (16) et (17), et correspondent à la convection, tandis que et correspondent à la diffusion.
On considère qu’aux interfaces entre le sang et la peau, ou entre la peau et la chambre de mesure, ou entre la chambre de mesure et la chambre de collecte, la vitesse est nulle de part et d’autre de l’interface : . En effet, la convection n’intervient que dans la chambre de collecte. Dans la peau et dans la chambre de mesure, le CO2 se propage uniquement par diffusion.
En combinant (16) et (17), du fait de la continuité du flux, et en introduisant un facteur de transparence, on obtient alors :
Où correspond au facteur de transparence entre le compartiment i et le compartiment j. Le facteur de transparence correspond à un ratio entre la surface ouverte de l’interface, laissant passer le gaz d’intérêt, sur la surface totale de l’interface.
De part et d’autre de l’interface, la pression est identique :
Or
Où est la température du milieu i, et
Le symbole signifie « égalité qui sert de définition».
Pris
e en compte de l’écoulement du sang
dans le compa
r
timent sang
Dans le compartiment sang B : ,
Prise en compte de l’admission d’air dans la chambre de collecte.
Dans la chambre de collecte, l’admission d’air se traduit par un terme source correspondant au gaz carbonique présent dans l’air ambiant
Discrétisation du modèle.
Le dispositif et le milieu sur lequel il est appliqué ont été discrétisés
L’équation de convection – diffusion (11) peut être discrétisée spatialement comme représenté sur la . Un maillage spatial a été réalisé, générant 17 points d’échantillonnage équirépartis dans chaque compartiment, entre le compartiment sang B et la chambre de collecte. Chaque compartiment comporte 5 points de d’échantillonnage, dont 2 sont placés au niveau d’une interface inférieure ou supérieure. Au niveau d’une interface entre un compartiment aval et un compartiment amont, la concentration sortant dans le compartiment aval est référencée par un exposant « out » tandis que la concentration entrant dans le compartiment amont est référencée par un exposant « in ».
La concentration à l’intérieur du compartiment sang est discrétisée selon 3 points d’échantillonnage :
La concentration à l’intérieur du compartiment peau est discrétisée selon 3 points d’échantillonnage:
La concentration à l’intérieur de la chambre de mesure est discrétisée selon 3 points d’échantillonnage:
La concentration à l’intérieur et en sortie de la chambre de collecte est discrétisée selon 4 points d’échantillonnage:
Au niveau des interfaces entre les compartiments, les concentrations sont
Nous définissons des concentrations de CO2 aux interfaces :
- Au niveau de l’interface sang / peau :
- Au niveau de l’interface peau/cellule de mesure
- Au niveau de l’interface cellule de mesure / cellule de collecte
La concentration d’entrée est . Dans la formulation du problème direct, est connu. Lorsqu’on résout le problème inverse, c’est qui est recherché.
L’écart spatial entre chaque point du maillage d’un même compartiment est identique. Par la suite, l’écart entre chaque point du maillage est noté , l’indice i désignant le milieu : i = 1 pour le sang, i = 2 pour la peau, i = 3 pour la chambre de mesure et i = 4 pour la chambre de collecte.
Compte tenu de la discrétisation spatiale effectuée, l’équation de convection diffusion (expression 11) devient :
Dans cet exemple, : on utilise un schéma d’intégration dit implicite.
On définit un vecteur , de dimension (17, 1), qui comporte l’ensemble des points de maillage.
Le nombre total d’échantillons est tel que
où est le nombre d’échantillons au sein du compartiment .
Où est la hauteur du milieu et la distance entre deux points successifs dans le milieu est noté .
La position suivant z dans la colonne, à laquelle commence chaque milieu, est définie par . Cette position peut être calculée de manière générale par les relations suivantes:
On définit un vecteur par rapport aux concentrations évaluées dans différents points d’échantillonnage définis par le vecteur .
Dans cet exemple, on suppose que le CO2 est mesuré en sortie de la chambre de mesure :
Dans notre cas, nous avons : donc et
On définit un vecteur par rapport aux concentrations évaluées aux différents points d’état.
On définit un vecteur d’état augmenté de dimension en ajoutant, au vecteur la concentration en entrée du milieu sanguin:
L’expression (30) peut être écrite de façon matricielle :
Où :
Et sur les termes de rang ( , ) : (39)
En dehors des points de rang (q,q+1) et (q, q-1), les termes de ont une valeur nulle.
et contient le terme puit (i-e la sortie du sang au niveau de la peau dans le milieu sanguin).
- sur la diagonale principale (ligne et colonne de même rang
- au-dessus de la diagonale principale
- en-dessous de la diagonale principale
En dehors des points de rang (q,q+1) et (q, q-1), (q,q), les termes de ont une valeur nulle.
La matrice ,décrivant les conditions aux interfaces entre deux compartiments adjacents, de dimension [ ] ( = 13 dans cet exemple), est calculée comme suit :
Tous les termes de la matrice sont nuls, sauf :
- Les termes [1,1] et [1,2], qui correspondent au bord inférieur, à la coordonnée
- Les termes [(
- Les termes [(
- Les termes [(
- Les termes (
Les termes de la matrice différents de ceux précédemment explicités sont égaux à 0.
Chaque terme non nul de la matrice est tel que :
avec .
Chaque terme non nul de la matrice est tel que :
avec
et
et
La matrice ,décrivant les conditions aux interfaces entre deux compartiments adjacents, de dimension [ ] ( = 13 dans cet exemple), est calculée comme suit. Le seul terme non nul est le terme situé à la position (3, 2) correspondant à l’interface entre le sang et la peau
Le terme correspond à la sortie par convection du compartiment du sang
L’arrivée d’air dans la cellule de collecte est prise en compte dans la matrice de commande décrite par la suite.
Dans cet exemple, on met en œuvre un estimateur récursif pour estimer la concentration de CO2 dans le sang, en fonction du temps. Il s’agit ici d’un estimateur récursif linéaire, de type filtre de Kalman.
L’estimateur est basé sur le modèle d’espace état spatio-temporel précédemment décrit, auquel on ajoute une équation d’observation, qui traduit la relation entre chaque mesure et les états.
De façon connue avec ce type d’estimateur, la concentration, à un instant, est calculée de manière récursive, à partir d’une estimée de ladite concentration, résultant d’une itération précédente, et d’une mesure effectuée à l’instant t.
En partant de l’équation (36) :
On peut écrire :
Par la suite, on utilise les notations suivantes :
Le procédé comporte une phase d’initialisation, au cours de laquelle on utilise un vecteur d’état initial et une matrice de covariance de l’erreur à l’instant initial .
L’étape de correction est définie par les équations suivantes :
Le gain de Kalman est estimé par
Le gain de Kalman traduit les poids (i-e la confiance) donnés à la mesure bruitée et à son estimée
Estimation de la concentration à l’instant :
Estimation de la matrice de covariance de l’erreur d’estimation :
L’estimation de permet de connaître , ce dernier étant le terme de rang 1 du vecteur . A partir de , on peut estimer .
L’étape de prédiction est définie par les équations suivantes :
Prédiction de la concentration à l’instant k+1 :
Prédiction de la matrice de covariance de l’erreur
La variable que l’on cherche à estimer , se trouve dans le vecteur d’état à la position (terme de rang 1) Elle est liée à la variable d’état suivante par les termes de l’équation de transport. L’équation d’état qui lui est associée est :
Où est un paramètre caractérisant notre a priori sur le signal . Si notre a priori est que le signal en entrée est stationnaire, nous choisirons
L’expression (37) se traduit par :
où
Où la matrice de transition augmenté
e
est
Où pour le terme défini par la hauteur suivant l’axe z :
et
Pour un schéma implicite d’intégration dans le temps, l’équation d’état discrète s’écrit :
et
Où la matrice de dimension est définie par :
Où est la matrice identité de dimension
L’équation d’observation pour le système augmenté s’écrit :
Où :
Dans notre cas, comme l’observation est définie par une seule valeur, la matrice d’observation augmentée se réduit à un vecteur d’observation .
Où :
Le système d’équations est
L’algorithme de Kalman enchaine de façon récursive une étape de correction du vecteur d’état à un instant (ou ) et de prédiction du vecteur d’état
A partir de , la teneur du sang en CO2 peut être
La schématise les principales étapes de mise en œuvre d’un procédé selon l’invention.
Etape 100 : apposition du dispositif de mesure contre le milieu analysé et activation de la pompe.
Etape 110 : initialisation : Au cours de cette étape, on définit un vecteur d’état initial et une matrice de covariance de l’erreur à l’instant initial . L’instant de mesure k est ensuite incrémenté : k = 2.
Etape 120 : mesure de la concentration de CO2 dans la chambre de mesure. On obtient ainsi la grandeur mesurée .
Etape 130 : A partir du vecteur d’état initial , ou du vecteur d’état estimé lors d’une itération précédente ,et de la matrice de covariance de l’erreur à l’instant initial ,ou résultant d’une itération précédente , détermination du gain de Kalman à l’instant k (expression 41) et estimation du vecteur d’état (expression 42).
L’estimation de permet de connaître , ce dernier étant le terme de rang du vecteur : cela correspond à l’étape 135. A partir de , on peut estimer (expression 47) : étape 136.
Etape 140 : Estimation de la matrice de covariance de l’erreur d’estimation (expression 43).
Etape 150 : Prédiction du vecteur d’état à l’instant k+1.(expression 44) et prédiction de la matrice de covariance de l’erreur à l’instant suivant
Etape 160 : incrémentation de l’instant de mesure
Etape 170 : réitération des étapes 110 à 160 ou sortie de l’algorithme
Les étapes itératives 110 à 170 sont réitérées jusqu’à la sortie de l’algorithme. La sortie de l’algorithme peut être décidée par l’utilisateur, ou correspondre à un nombre prédéterminé d’itérations.
On a mis en œuvre le procédé précédemment décrit en prenant en compte les paramètres suivants :
- longueur (selon l’axe X) compartiment sang:
- largeur (selon l’axe Y) compartiment sang :
- hauteur compartiment sang (selon l’axe Z) :
- surface du compartiment sang au contact de la peau :
- volume du compartiment sang :
- Coefficient de Henry du gaz carbonique dans le sang :
- Coefficient de diffusion du gaz carbonique dans le sang (coefficient de l’eau) :
- Concentration initiale en CO2dans le sang :
- Débit sanguin :
- Vitesse transverse (selon x) :
- Vitesse axiale (selon z) :
- Vitesse transverse (selon x) :
- longueur compartiment peau(selon l’axe X):
- largeur compartiment peau (selon l’axe Y):
- hauteur compartiment peau (selon l’axe Z):
- surface de la peau en contact avec le dispositif :
- Coefficient de Henry du gaz carbonique dans la peau :
- Coefficient de diffusion du gaz carbonique dans la peau : :
- température de la chambre de mesure :
- longueur de la chambre de mesure (selon l’axe X):
- largeur de la chambre de mesure (selon l’axe Y):
- hauteur de la chambre de mesure (selon l’axe Z):
- Coefficient de Henry du gaz carbonique dans la chambre de mesure :
- Coefficient de diffusion du gaz carbonique dans la chambre de mesure :
- Concentration en Co2dans l’air entrant, on suppose qu’il y a un filtre
- Débit d’air convectif dans la chambre de collecte :
- Vitesse transversale (selon l’axe x) :
- Vitesse axiale (selon l’axe z) :
- Vitesse transversale (selon l’axe x) :
- surface de la face de contact entre la chambre de mesure et la chambre de collecte :
- température de la chambre de collecte :
- longueur de la chambre de collecte (selon l’axe X) :
- largeur de la chambre de collecte (selon l’axe Y):
- hauteur de la chambre de collecte (selon l’axe Z):
- Coefficient de Henry du gaz carbonique dans la chambre de collecte :
- Coefficient de diffusion du gaz carbonique dans la chambre de collecte :
- Facteur de transparence entre peau et la cellule de mesure : rpeau|mes=1 ;
- Facteur de transparence entre la cellule de mesure et la cellule de collecte : rmes|col=1
- Pas d’échantillonnage temporel :
- Variance du bruit de mesure : 10-6mol/m3
- Variance du bruit du modèle : 10-8mol/m3
- Coefficient du modèle autorégressif :
On a tout d’abord mis en œuvre un modèle direct, permettant d’estimer la concentration en chaque point d’échantillonnage à partir d’une concentration connue dans le sang. La concentration connue dans le sang suit un créneau. Les figures 6A, 6B, 6C décrivent l’évolution de la concentration en fonction du temps respectivement dans le compartiment sang, dans la peau, et dans la chambre de collecte. La montre l’évolution de la concentration au niveau du capteur de gaz, dans la chambre de mesure (courbe a). La courbe b de la montre la concentration au niveau du capteur de gaz de la courbe a, à laquelle on a ajouté un bruit gaussien pour simuler des mesures bruitées. Sur chacune des figures 6A à 6D, l’axe des abscisses correspond au temps et l’axe des ordonnées correspond à la concentration de CO2 estimée.
A partir des mesures bruitées simulées (courbe b de la ), on a mis en œuvre l’estimateur récursif précédemment décrit. Le vecteur d’état estimé à chaque instant de mesure permet d’obtenir une estimation de la concentration en différents points d’échantillonnage. Les figures 7A, 7B, 7C, 7D montrent l’évolution de la concentration de CO2 respectivement au niveau de la chambre de collecte, de la chambre de mesure, du compartiment peau et du sang. La montre la concentration dans le sang réelle (courbe a) et la concentration dans le sang résultant de l’estimation.
Sur chacune des figures 7A à 7E, l’axe des abscisses correspond au temps et l’axe des ordonnées correspond à la concentration de CO2 estimée.
Bien que décrite en lien avec une mesure de CO2transcutanée, l’invention peut être mise en œuvre dans d’autres applications, de façon à mesurer un gaz émis par un milieu, solide ou liquide. Le milieu peut notamment être :
- un végétal, par exemple un fruit, de façon à suivre le procédé de maturation ;
- du lisier, de façon à suivre l’émission de gaz, par exemple le méthane ;
- un milieu de culture d’organismes ou de microorganismes biologiques ;
- de l’eau, par exemple de l’eau douce ou de l’eau de mer, par exemple pour suivre la concentration de gaz ou réguler la concentration de CO2;
un sol, de façon à en étudier la respiration.
Claims (10)
- Procédé d’estimation d’une teneur d’un gaz d’intérêt dans un milieu, à l’aide d’un dispositif de mesure, destiné à être disposé contre le milieu, le dispositif de mesure s'étendant entre une face de contact, destinée à être appliquée contre le milieu et une extrémité distale, le dispositif de mesure comportant une paroi latérale, s'étendant entre la face de contact et l'extrémité distale, le dispositif de mesure comportant :
- au niveau de la face de contact (10), au moins une ouverture d’admission (12), configurée pour collecter le gaz d'intérêt émis par le milieu, l'ouverture d’admission étant pratiquée à travers la face de contact;
- une chambre de mesure (20), comportant un capteur de gaz (23), le capteur de gaz étant configuré pour mesurer une concentration de gaz d’intérêt s’écoulant à travers la chambre de mesure;
- une chambre de collecte (30), reliée à la chambre de mesure, et délimitée par une ouverture sur la paroi latérale, la chambre de collecte comportant au moins une ouverture latérale (34), ménagée à travers la face latérale, ou sur la paroi supérieure de la chambre de collecte de façon à admettre de l’air ambiant dans la chambre de collecte ;
- la chambre de mesure (20) est disposée entre la face de contact (10) et la chambre de collecte (30) ;
- le dispositif comporte un moyen d’entraînement (41), configuré pour entraîner l’air de la chambre de collecte vers une ouverture d’évacuation (42), l’entraînement d’air induisant un transport du gaz d’intérêt de la face de contact vers la chambre de collecte, à travers la chambre de mesure :
- le procédé comportant :
- a) mesure de la concentration de gaz d’intérêt dans la chambre de mesure ;
- b) à l’aide d’une unité de traitement, modélisation du transport du gaz d’intérêt, entre le milieu et la chambre de collecte, la modélisation comportant une prise en compte de la diffusion du gaz depuis le milieu jusqu’à la chambre de collecte à travers la face de contact et le dispositif, ainsi qu’une convection du gaz dans la cellule de collecte résultant de l’entraînement produit par le moyen d’entraînement ;
- c) à partir de la mesure résultant de l’étape a), et d’une prise en compte du modèle résultant de l’étape b), estimation de la teneur du gaz d’intérêt dans le milieu.
- Procédé selon la revendication 1, dans lequel :
- le dispositif est spatialement discrétisé selon un maillage spatial, définissant des points de maillage entre la face de contact et la chambre de collecte ;
- le modèle est un modèle spatio-temporel discrétisé, de façon à estimer une teneur en gaz d’intérêt en différents points du maillage, et en différents instants.
- Procédé selon l’une quelconque des revendications précédentes, dans lequel l’étape b) comporte :
- modélisation de la diffusion du gaz d’intérêt à travers la chambre de mesure ;
- modélisation de la diffusion et de la convection du gaz d’intérêt dans la chambre de collecte.
- Procédé selon l’une quelconque des revendications précédentes, dans lequel l’étape b) comporte une modélisation du transport du gaz d’intérêt dans le milieu.
- Procédé selon l’une quelconque des revendications précédentes, dans lequel l’étape c) met en œuvre un estimateur récursif.
- Procédé selon la revendication 5, dans lequel l’étape c) met en œuvre un estimateur récursif linéaire.
- Procédé selon la revendication 6, dans lequel l’estimateur récursif linéaire est un filtre de Kalman.
- Procédé selon l’une quelconque des revendications précédentes, dans lequel le milieu est la peau d’un utilisateur, la peau s’étendant en dessus d’un vaisseau sanguin, dans lequel l’étape b) comporte :
- à partir de la concentration de gaz d'intérêt dans la chambre de mesure, résultant de l'étape a), estimation d'une concentration de gaz d'intérêt transcutané ;
- à partir de la concentration de gaz d'intérêt transcutané résultant de la sous-étape (i), estimation d'une concentration ou pression partielle de gaz d'intérêt dissous dans le sang de l'utilisateur.
- Procédé selon l’une quelconque des revendications précédentes, comportant une modélisation du transport du gaz d’intérêt émis par le milieu à travers le dispositif.
- Procédé selon l’une quelconque des revendications précédentes, dans lequel le milieu à analyser est un milieu solide ou liquide.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR2208595A FR3138999A1 (fr) | 2022-08-28 | 2022-08-28 | Procédé d’estimation d’une concentration de gaz dégagé par un milieu |
FR2212730A FR3139000A1 (fr) | 2022-08-28 | 2022-12-03 | Procédé d’estimation d’une concentration de gaz dégagé par un milieu. |
PCT/EP2023/073458 WO2024046930A1 (fr) | 2022-08-28 | 2023-08-27 | Procédé d'estimation d'une concentration de gaz dégagé par un milieu |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR2208595 | 2022-08-28 | ||
FR2208595A FR3138999A1 (fr) | 2022-08-28 | 2022-08-28 | Procédé d’estimation d’une concentration de gaz dégagé par un milieu |
Publications (1)
Publication Number | Publication Date |
---|---|
FR3138999A1 true FR3138999A1 (fr) | 2024-03-01 |
Family
ID=85017967
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR2208595A Pending FR3138999A1 (fr) | 2022-08-28 | 2022-08-28 | Procédé d’estimation d’une concentration de gaz dégagé par un milieu |
FR2212730A Pending FR3139000A1 (fr) | 2022-08-28 | 2022-12-03 | Procédé d’estimation d’une concentration de gaz dégagé par un milieu. |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR2212730A Pending FR3139000A1 (fr) | 2022-08-28 | 2022-12-03 | Procédé d’estimation d’une concentration de gaz dégagé par un milieu. |
Country Status (1)
Country | Link |
---|---|
FR (2) | FR3138999A1 (fr) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010148692A (ja) * | 2008-12-25 | 2010-07-08 | National Cardiovascular Center | 表面ガス検知法および検知装置 |
WO2017023500A1 (fr) * | 2015-08-05 | 2017-02-09 | Vail Scientific, Llc | Appareil et procédé permettant la surveillance non-invasive d'oxyde nitrique et d'autres gaz dans le sang |
WO2020249466A1 (fr) | 2019-06-09 | 2020-12-17 | Commissariat à l'énergie atomique et aux énergies alternatives | Dispositif portable d'estimation de la pression partielle de gaz sanguin |
-
2022
- 2022-08-28 FR FR2208595A patent/FR3138999A1/fr active Pending
- 2022-12-03 FR FR2212730A patent/FR3139000A1/fr active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010148692A (ja) * | 2008-12-25 | 2010-07-08 | National Cardiovascular Center | 表面ガス検知法および検知装置 |
WO2017023500A1 (fr) * | 2015-08-05 | 2017-02-09 | Vail Scientific, Llc | Appareil et procédé permettant la surveillance non-invasive d'oxyde nitrique et d'autres gaz dans le sang |
WO2020249466A1 (fr) | 2019-06-09 | 2020-12-17 | Commissariat à l'énergie atomique et aux énergies alternatives | Dispositif portable d'estimation de la pression partielle de gaz sanguin |
Non-Patent Citations (1)
Title |
---|
"Abstracts", FUEL AND ENERGY ABSTRACTS, ELSEVIER, AMSTERDAM, NL, vol. 61, no. 6, 14 October 2020 (2020-10-14), pages 562 - 661, XP086295297, ISSN: 0140-6701, [retrieved on 20201014], DOI: 10.1016/J.FUELENEAB.2020.09.002 * |
Also Published As
Publication number | Publication date |
---|---|
FR3139000A1 (fr) | 2024-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6040578A (en) | Method and apparatus for multi-spectral analysis of organic blood analytes in noninvasive infrared spectroscopy | |
FR2905762A1 (fr) | Systeme et procede de detection et d'analyse d'au moins une substance volatile dans des echantillons d'haleine d'un sujet | |
EP3979911B1 (fr) | Dispositif portable d'estimation de la pression partielle de gaz sanguin | |
US7640140B2 (en) | Method of processing noninvasive spectra | |
US20150289782A1 (en) | Portable breath volatile organic compounds analyser and corresponding unit | |
US20030023170A1 (en) | Optically similar reference samples and related methods for multivariate calibration models used in optical spectroscopy | |
WO1997028437A9 (fr) | Procede et appareil d'analyse multispectrale en spectroscopie a infrarouge proche non vulnerante | |
JPH08136448A (ja) | 散乱吸収体内の散乱特性・吸収特性の測定方法及び装置 | |
JP7295167B2 (ja) | 全血試料内のパーセント酸素飽和度を測定する方法 | |
US11105739B2 (en) | Method and system for analyzing a sample using cavity ring-down spectroscopy, and a method for generating a predictive model | |
CN108872100B (zh) | 一种多次增强光谱高精度氨气检测装置及检测方法 | |
FR3138999A1 (fr) | Procédé d’estimation d’une concentration de gaz dégagé par un milieu | |
JP2004150984A (ja) | 近赤外分光による溶質濃度測定方法とその装置 | |
WO2024046930A1 (fr) | Procédé d'estimation d'une concentration de gaz dégagé par un milieu | |
FR3139001A1 (fr) | Dispositif portable d'estimation d’une concentration de gaz dégagé par un milieu. | |
JPH11178799A (ja) | 生体表層組織の分析方法及びその装置 | |
Koo | Measurement of blood analytes in turbid biological tissue using near-infrared Raman spectroscopy | |
Ageev et al. | Optoacoustic gas-analysis for diagnostics of biosystems | |
Kottmann et al. | New photoacoustic cell with diamond window for mid-infrared investigations on biological samples | |
JPWO2021260752A5 (fr) | ||
WO2024133307A1 (fr) | Systèmes d'analyse d'un échantillon biologique gazeux par spectroscopie térahertz dans le domaine temporel |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 2 |
|
PLSC | Publication of the preliminary search report |
Effective date: 20240301 |