FR 3077387 A1 20190802 FR 1850683 A 20180129 Procédé d'estimation d'une quantité de gaz Description DOMAINE TECHNIQUE Le domaine technique de l'invention est l'analyse d'un gaz par une méthode optique. ART ANTERIEUR Le recours à des méthodes optiques pour l'analyse d'un gaz est assez fréquent. Des dispositifs permettent de déterminer la composition d'un gaz en se basant sur le fait que les espèces composant un gaz présentent des propriétés spectrales d'absorption différentes les unes des autres. Ainsi, connaissant une bande spectrale d'absorption d'une espèce gazeuse, sa concentration peut être déterminée par une estimation de l'absorption de la lumière traversant le gaz, en utilisant la loi de Beer Lambert. Ce principe permet une estimation de la concentration d'une espèce gazeuse présente dans le milieu. La demande de brevet W02007064370 repose sur un tel principe. Ce document décrit un dispositif d'analyse d'un gaz comportant une source de lumière et un photodétecteur. Un filtre optique est disposé entre la source de lumière et le photodétecteur. Ce filtre présente des bandes passantes étroites, typiquement de largeur inférieure à 500 nm, correspondant aux bandes spectrales d'absorption des espèces gazeuses à détecter. Un tel filtre est désigné par le terme "multibandes". Mais une telle configuration présente certains inconvénients. Tout d'abord, la réalisation d'un filtre optique présentant des bandes passantes étroites, séparées les unes des autres, est coûteuse. En outre, un tel filtre peut présenter des difficultés techniques de fabrication. Ce document indique d'ailleurs que lorsque le nombre de bandes spectrales augmente, le recours à un tel filtre est difficilement envisageable. De plus, cette configuration n'est pas versatile et ne fonctionne que pour les espèces gazeuses correspondant au filtre optique réalisé. Il n'est pas possible de quantifier la concentration d'une espèce gazeuse dont la bande spectrale d'absorption ne correspond pas à une bande passante du filtre, sauf à changer le filtre. Par ailleurs, certaines espèces gazeuses présentent des bandes spectrales d'absorption se chevauchant. L'utilisation du filtre multibandes ne permet pas de les séparer et le dispositif ne permet pas de discriminer ces espèces, c'est-à-dire de déterminer une concentration de chacune de ces espèces. Le document US2009/213380 décrit un dispositif d'analyse de gaz, dans lequel un gaz est illuminé successivement selon différents spectres d'illumination, à l'aide d'une source de 2 lumière. A chaque illumination, le spectre d'illumination du gaz est modifié en interposant successivement différents filtres entre la source de lumière et le gaz à analyser. Un tel dispositif nécessite d'utiliser successivement plusieurs filtres, entre la source de lumière et le gaz, ce qui suppose des moyens mécaniques permettant une commutation entre les différents filtres. On 5 comprend qu'une telle configuration est complexe et n'est pas appropriée à un détecteur de gaz compact et peu onéreux. Elle suppose, d'autre part, de disposer d'un a priori sur les espèces gazeuses présentes dans le gaz, de façon à sélectionner les filtres mis en oeuvre. L'invention permet de surmonter ces inconvénients en proposant un dispositif et un procédé d'analyse d'un gaz simple, peu coûteux, et pouvant adresser différentes espèces gazeuses, y 10 compris les espèces gazeuses dont les bandes spectrales d'absorption se recouvrent. De plus, le dispositif est aisément adaptable à différents mélanges gazeux ; il ne nécessite pas la réalisation d'un filtre multibandes spécifique à chaque mélange gazeux étudié. Par ailleurs, le nombre d'espèces gazeuses présentes dans le mélange n'est pas limité par des contraintes de conception d'un filtre. En outre, l'invention ne nécessite pas des moyens de calculs puissants. En cela, elle 15 s'applique particulièrement à des capteurs de gaz nomades, dont la consommation en énergie doit être minimisée. EXPOSE DE L'INVENTION Un premier objet de l'invention est un procédé d'analyse d'un gaz, pour estimer une quantité de plusieurs espèces gazeuses différentes, dites espèces gazeuses d'intérêt, potentiellement 20 présentes dans le gaz analysé, chaque espèce gazeuse étant apte engendrer une absorption d'un rayonnement lumineux dans une bande spectrale d'absorption, le procédé comportant les étapes suivantes : a) disposition du gaz analysé entre une source de lumière et un photodétecteur de mesure, la source de lumière étant apte à émettre un rayonnement d'illumination se propageant 25 à travers le gaz, le photodétecteur de mesure étant apte à détecter un rayonnement transmis par le gaz et à délivrer un signal de détection représentatif d'une intensité du rayonnement détecté ; b) illumination du gaz analysé à l'aide de la source de lumière, le rayonnement d'illumination s'étendant selon une bande spectrale, dite bande spectrale 30 d'illumination, recouvrant la bande spectrale d'absorption d'une pluralité d'espèces gazeuses d'intérêt, à chaque illumination correspondant un spectre, dit spectre d'illumination, du rayonnement d'illumination ; 3 c) lors de l'illumination, détection, par le photodétecteur de mesure, d'une intensité d'un rayonnement transmis par le gaz analysé, dans au moins une bande spectrale de détection; d) réitération des étapes b) et c), en modifiant le spectre d'illumination du rayonnement 5 d'illumination; e) estimation d'une quantité de chacune des espèces gazeuses d'intérêt en fonction de l'intensité du rayonnement détecté lors de chaque illumination, chaque intensité étant associée à un spectre d'illumination; le procédé étant caractérisé en ce que l'étape e) comporte : 10 - une prise en compte, pour chaque spectre d'illumination, d'une fonction de réponse paramétrique, la fonction de réponse paramétrique étant définie par des paramètres, de telle sorte qu'à chaque spectre d'illumination sont associés des paramètres, les paramètres associés à deux spectres d'illumination distincts étant différents et préalablement définis ; 15 - une estimation d'une grandeur représentative de l'absorption du rayonnement d'illumination par le gaz, lors de chaque illumination, à l'aide de la fonction paramétrique. La bande spectrale de détection peut s'étendre entre une longueur d'onde minimale jusqu'à une longueur d'onde maximale. Elle peut notamment s'étendre continûment entre une longueur 20 d'onde minimale et une longueur d'onde maximale. Les étapes b) et c) peuvent être répétées selon un nombre d'itérations prédéterminé. A chaque spectre d'illumination peuvent être associés des paramètres dits spécifiques, définis pour chaque espèce gazeuse d'intérêt, au moins deux paramètres spécifiques respectivement définis pour deux espèces gazeuses d'intérêt distinctes étant différents les uns des autres. 25 L'étape e) peut comporter une estimation, lors de chaque illumination, d'une intensité, dite intensité non atténuée, qui serait mesurée par le photodétecteur de mesure en l'absence de gaz. La grandeur représentative de l'absorption du gaz peut alors être, pour chaque illumination, l'intensité mesurée lors de l'étape c), éventuellement normalisée par l'estimation de l'intensité non atténuée 30 L'étape e) peut comporter, pour chaque illumination, une comparaison entre : - la grandeur représentative de l'absorption du gaz estimée, pour chaque illumination, lors de l'étape e) ; - et une mesure de ladite grandeur, à l'aide de l'intensité mesurée lors de l'étape c). 4 Les quantités d'espèces gazeuses déterminées sont alors obtenues en fonction de chaque comparaison. On dispose d'une comparaison pour chaque spectre d'illumination. Chaque comparaison peut prendre la forme d'un écart quadratique, les quantités déterminées étant celles pour lesquelles une combinaison arithmétique de chaque comparaison se rapproche 5 d'une valeur prédéterminée. La combinaison arithmétique peut être une somme, les quantités déterminées étant celles pour lesquelles la somme est minimale. Il est avantageux que l'étape e) mette en oeuvre une méthode par moindres carrés non linéaire. Le procédé peut comporter également, lors de chaque illumination, une détection, par un photodétecteur de référence d'une intensité de référence d'un rayonnement de référence, 10 considéré comme n'ayant pas été absorbé par le gaz. L'étape e) comporte alors une prise en compte des intensités de référence respectivement détectées au cours de chaque illumination. La source de lumière peut être portée à différentes températures. Chaque spectre d'illumination correspond alors à une température de la source de lumière. Le procédé peut comporter, préalablement à l'étape e), une phase de calibration, effectuée pour 15 chaque spectre d'illumination, de façon à définir la fonction de réponse paramétrique associée à chaque spectre d'illumination, la phase de calibration comportant les étapes suivantes: cal-i) réalisation ou simulation des étapes a) à c) en considérant un gaz de calibration, comportant une quantité connue d'au moins une espèce gazeuse d'intérêt, et en considérant un spectre d'illumination, de façon à obtenir ou à calculer une intensité 20 détectée par le photodétecteur de mesure ou par un photodétecteur de calibration considéré comme représentatif du photodétecteur de mesure ; cal-ii) réitération de l'étape cal-i) en modifiant la quantité de l'espèce gazeuse d'intérêt ; cal-iii) prise en compte d'un modèle paramétrique de la fonction de réponse; cal-iv) à partir des intensités résultant de chaque étape cal-i) et du modèle paramétrique 25 prise en compte lors de l'étape cal-iii), détermination d'au moins un paramètre de la fonction de réponse. Ces paramètres sont associés au spectre d'illumination pris en compte lors de l'étape cal-i). Les étapes cal-i) à cal-iv) peuvent être effectuées en illuminant le gaz de calibration selon différents spectres d'illumination, de façon à obtenir des paramètres de la fonction de réponse 30 pour chaque spectre d'illumination. La phase de calibration peut être réalisée en utilisant une source de lumière de calibration et un photodétecteur de mesure de calibration, respectivement considérés comme similaires à la source de lumière et au photodétecteur de mesure. Certains paramètres peuvent être 5 déterminés par une calibration effectuée à l'aide d'un photodétecteur de calibration, d'autres paramètres étant déterminés au cours d'une calibration effectuée à l'aide du photodétecteur de mesure. Selon un mode de réalisation, la bande spectrale de détection recouvre des bandes spectrales 5 d'absorption de plusieurs espèces gazeuses d'intérêt. Un autre objet de l'invention est un dispositif d'analyse d'un gaz comportant : - une source de lumière, apte à émettre un rayonnement d'illumination, la source de lumière étant également apte à être modulée de telle sorte que le spectre du rayonnement d'illumination émis, dit spectre d'illumination, est variable dans une 10 même bande spectrale d'illumination ; - un photodétecteur de mesure, apte à détecter un rayonnement lumineux selon une bande spectrale de détection; - une enceinte, apte à recevoir le gaz, agencée de telle sorte que le photodétecteur de mesure soit configuré pour détecter un rayonnement lumineux ayant traversé le gaz ; 15 un processeur, configuré pour mettre en oeuvre l'étape e) d'un procédé selon le premier objet de l'invention, à partir du rayonnement détecté par le photodétecteur de mesure, suite à différentes illuminations du gaz, respectivement selon différents spectres d'illumination. D'autres avantages et caractéristiques ressortiront plus clairement de la description qui va suivre de modes particuliers de réalisation de l'invention, donnés à titre d'exemples non limitatifs, et 20 représentés sur les figures listées ci-dessous. FIGURES La figure 1 représente un exemple de dispositif d'analyse d'un gaz selon un mode de réalisation. La figure 2 montre l'évolution du spectre d'illumination d'une source de lumière de type corps noir en fonction de la température. 25 La figure 3 représente des spectres d'absorbance de différentes espèces gazeuses. La figure 4 illustre les principales étapes d'un procédé d'analyse d'un gaz selon l'invention. Les figures 5A et 5B représentent d'autres configurations de dispositifs d'analyse permettant une mise en oeuvre de l'invention. Les figures 6A à 6G montrent des résultats de simulations de fonctionnement d'un dispositif 30 d'analyse, le gaz analysé étant un mélange CO2 et CH4. EXPOSE DE MODES DE REALISATION PARTICULIERS 6 La figure 1 est un exemple d'un dispositif d'analyse 1 de gaz. Le gaz G est un mélange gazeux comportant une pluralité d'espèces gazeuses d'intérêt K étant un entier positif quantifiant le nombre d'espèces gazeuses différentes présentes dans le mélange. Le dispositif comporte une enceinte 10 définissant un espace interne à l'intérieur duquel se 5 trouvent : - une source de lumière 11, apte à émettre successivement un rayonnement d'illumination 12, selon différents spectres d'illumination Sj, de façon à illuminer un gaz G s'étendant dans l'espace interne ; - un photodétecteur de mesure 20, apte à détecter un rayonnement lumineux 14 10 transmis par le gaz G sous l'effet de son illumination par le rayonnement d'illumination 12; - un photodétecteur de référence 20ref, apte à détecter un rayonnement lumineux 14,-ef transmis par le gaz G dans une bande spectrale de référence 4X,-ef sous l'effet de son illumination par le rayonnement d'illumination 12. 15 L'enceinte peut comporter une admission 10,, et une évacuation 10ex permettant respectivement l'admission du gaz dans l'enceinte et son évacuation. L'admission et/ou l'évacuation du gaz peut être forcée ou spontanée. Ainsi, le gaz peut traverser l'enceinte par un simple effet de diffusion. Le photodétecteur de mesure 20 est associé à un filtre de mesure 18. Le photodétecteur de 20 référence 20ref est associé à un filtre de référence 18,-ef. Les fonctions du filtre de mesure 18 et du filtre de référence 18,-ef sont précisées dans la suite de la description, en lien avec l'équation (3). Le photodétecteur de mesure 20 détecte le rayonnement lumineux transmis 14 dans une bande spectrale de détection A20, cette dernière étant limitée par le filtre de mesure 18. Le photodétecteur de mesure 20 et le photodétecteur de référence 20ref peuvent être 25 spectralement non résolus. Il peut en particulier s'agir de photodiodes ou de thermopiles. Le photodétecteur de référence est optionnel. La source de lumière 11 est apte à émettre le rayonnement d'illumination 12 selon une pluralité de spectres d'illumination différents S1.. Sj.. Sj, différents les uns des autres, dans une même bande spectrale d'illumination 3,12 où J est un entier désignant le nombre de spectres. J j 30 est un entier indiçant chaque spectre d'illumination Sj, avec 1 j J. Par spectre d'illumination, on entend un histogramme de l'intensité /(À) d'un rayonnement lumineux en fonction de la longueur d'onde À. Le spectre de la source de lumière définit une bande spectrale 7 d'illumination Al2, cette dernière pouvant s'étendre entre le proche ultraviolet et l'infra-rouge moyen, entre 200 nm et 10 p.m, et le plus souvent entre 1 p.m et 10 p.m. A chaque illumination j par la source de lumière correspond un spectre d'illumination Sj. Le dispositif comporte un processeur 30, par exemple un microprocesseur. Ce dernier est 5 configuré pour recevoir un signal représentatif d'une intensité lumineuse détectée par le photodétecteur de mesure 20, et de l'éventuel photodétecteur de référence 20'f. Il permet d'exécuter les opérations de traitement décrites ci-après, de façon à déterminer une quantité de chaque espèce gazeuse recherchée. Le processeur 30 est relié à une mémoire 32 comportant des instructions permettant la mise en oeuvre de certaines étapes du procédé décrit ci-après. 10 Selon un mode de réalisation préféré, les différents spectres d'illumination Sj sont obtenus en modulant une température de la source de lumière 11. Cette dernière peut alors être une source de lumière de type corps noir ou corps gris. Selon un autre mode de réalisation, les différents spectres d'illumination sont obtenus en utilisant une source de lumière 11 regroupant une pluralité de sources de lumières élémentaires 111....11N, où N est un entier positif désignant le 15 nombre de sources de lumières élémentaires constituant la source. Les sources de lumière élémentaires peuvent être des diodes électroluminescentes dont l'éclairement peut être modulé de façon à former différents spectres d'illumination dans la bande spectrale d'illumination Al2. La figure 2 montre différents spectres d'illumination Sj d'une source de lumière 11, assimilée à 20 une source dite de type corps noir, en fonction d'une température Tj à laquelle est portée la source de lumière. Les spectres Sj suivent la loi de Planck, décrivant l'émission d'un corps noir 21/c hc 2 selon une expression de type L(À, Tj) = , (1) où Â.5 eÂk T -1 - I, (il, Tj) est la luminance spectrale, dépendant de la longueur d'onde À et de la température de surface Tj du corps noir, 25 - h est la constante de Planck, - k est la constante de Boltzmann, - c est la vitesse de la lumière dans le vide. Une telle source de lumière est par exemple obtenue par un simple filament chauffé à une température comprise par exemple entre 400°C et 800°C. Dans le cas de la figure 2, on a 30 représenté l'évolution de la luminance spectrale «À, Tj) dans une plage de température 8 comprise entre 600°C et 800 °C et selon une bande spectrale d'illumination 3,12 comprise entre 1p.met1011m. Ce type de source de lumière est particulièrement avantageux, car cela permet de moduler le spectre d'illumination Sj par une simple modulation de la température Tj de la source. A chaque 5 température Tj est associé un spectre d'illumination Sj. La modulation de la température de source peut être par exemple obtenue par une simple modulation d'un courant d'alimentation de la source de lumière. De plus, la modulation du spectre d'illumination est non linéaire. Par modulation non linéaire, on entend qu'en considérant respectivement un premier spectre d'illumination Sj et un deuxième spectre d'illumination Sj+i, la modulation est telle que le ratio 10 des intensités du spectre d'illumination entre une première longueur d'onde À1 et une deuxième longueur d'onde À2 est différent d'un spectre d'illumination à l'autre. Ceci peut être explicité par l'expression suivante : si (4) si-Ei (4) sj(A.2) sj+i(A.2) (2) où 5101.1) désigne la valeur de l'intensité, à une longueur d'onde .11, du rayonnement 15 d'illumination 12 de spectre Sj. La plupart des espèces gazeuses connues absorbent tout ou partie d'un rayonnement électromagnétique selon un spectre d'absorption présentant une ou plusieurs bandes spectrales d'absorption Ak, ainsi qu'un ou plusieurs pics d'absorption. L'absorption d'une espèce gazeuse Gk peut être définie selon la loi de Beer-Lambert, selon laquelle : 20 /(À) = Eray /0 (À, r ay)e-1.1k(A)xray (3) où - /(À) désigne l'intensité transmise par l'espèce gazeuse Gk à la longueur d'onde À, c'est-à-dire l'intensité du rayonnement détecté 14 à la longueur d'onde À; r ay) désigne l'intensité incidente à l'espèce gazeuse Gk à la longueur d'onde À, 25 c'est-à-dire l'intensité du rayonnement d'illumination 12 à la longueur d'onde À, selon une trajectoire ray prédéterminée. Cela correspond à l'intensité du rayonnement parvenant au photodétecteur de mesure 20, selon la trajectoire ray, sans atténuation par l'espèce gazeuse Gk. Il s'agit d'une intensité dite non atténuée. - 14(.1) désigne un coefficient d'absorption dit linéaire de l'espèce gazeuse Gk 30 considérée, à la longueur d'onde À ; - xray désigne la distance traversée par chaque rayonnement lumineux selon chaque trajectoire ray entre la source de lumière 11 et le photodétecteur de mesure 20. 9 Le filtre de mesure 18, auquel est associé le photodétecteur de mesure 20, permet de délimiter la bande spectrale de mesure du rayonnement 14 transmis par le gaz. Cela permet par exemple d'éviter que les mesures soient influencées par des bandes d'absorption de la vapeur d'eau, ces dernières étant situées autour de 2.6 p.m ainsi qu'entre 5.5 et 7.5 p.m. 5 Le filtre optique de référence 18ref définit la bande spectrale de référence 4Xref, correspondant à une plage de longueurs d'onde non absorbées, ou de façon négligeable, par le gaz G analysé. La bande spectrale de référence est par exemple centrée autour de la longueur d'onde de référence Xref 3.91 p.m. Le photodétecteur de référence 20ref, permet de mesurer une intensité de référence /J,'f. A partir de cette dernière, on peut estimer une intensité émise par la 10 source de lumière 11 dans la bande spectrale de référence 4Xref. La connaissance du spectre d'émission Sj permet, à partir de /J,'f, d'obtenir les intensités d'émission aux différentes longueurs d'onde de la bande spectrale d'illumination 412 et de la bande spectrale de détection 420. Ainsi, l'intensité de référence /J,'f permet d'estimer une intensité /0,/, non atténuée, qui serait mesurée par le photodétecteur de mesure 20 en l'absence d'atténuation par le gaz G dans 15 l'enceinte 10. Le recours à un photodétecteur de référence 20ref, pour estimer l'intensité non atténuée /0,/ mesurée en l'absence d'atténuation par le gaz, est utile lorsque l'intensité de la source de lumière varie dans le temps, notamment en raison d'un vieillissement de la source de lumière. Le coefficient d'absorption linéaire µk(À) peut être exprimé sous la forme d'un produit d'un 20 coefficient d'absorption linéaire molaire pk'(À) par une quantité ck de l'espèce gazeuse Gk. Le coefficient d'absorption linéaire molaire représente le coefficient d'absorption linéaire par mole de l'espèce gazeuse considérée. L'expression (3) devient : 1(À) = Eray 0(À,ray)e-14(A.)CkXraY (4) 25 Ainsi, / (À) gray - 0 .--, 10 (À, ray)hk(ck,À, Xray) (5) = / ray)e-14(A)ckxray = gray ik(À.)CkXray = e-1-1*(A)Xray (6), hk où hk (ck, À, = e-1-` désignant une fonction spectrale d'absorption de l'espèce gazeuse Gk lorsqu'un rayonnement de longueur d'onde À parcourt une distance Xray dans le gaz. Cette fonction représente une proportion du rayonnement d'illumination absorbé par l'espèce gazeuse Gk. 30 La figure 3 montre les spectres d'absorption, c'est-à-dire l'évolution du coefficient d'absorption linéaire molaire 4(À) de différentes espèces gazeuses, en fonction de la longueur d'onde. On 10 constate que ces spectres d'absorption présentent des pics d'absorption dans différentes bandes spectrales, dites bandes spectrales d'absorption Ak. Par pic d'absorption, on entend un maximum local de l'évolution de l'absorption en fonction de la longueur d'onde. Par bande spectrale d'absorption Ak, on entend une bande spectrale dans 5 laquelle l'absorption est considérée comme significative, c'est-à-dire supérieure à 50 Mol-1.cm 1, ou 100 Mol-1.cm-1 voire 200 Mol-lcm-1 lorsque l'absorption est quantifiée par le coefficient d'absorption linéaire molaire. Le CO2 dont on a représenté le dixième de la valeur du coefficient d'absorption linéaire molaire, présente une bande spectrale d'absorption entre 4.2 p.m et 4.4 ûm. Le N20 présente une première bande spectrale d'absorption comprise entre 4.3 p.m et 4.6 10 p.m, ainsi qu'une deuxième bande spectrale d'absorption entre 7.5 p.m et 8 p.m. Le CO comprend une bande spectrale d'absorption entre 4.5 p.m et 4.8 p.m. On constate que les bandes spectrales d'absorption du N20 et du CO se recouvrent partiellement. Selon un mode de réalisation de l'invention le spectre d'illumination .S'j peut s'étendre selon une bande spectrale d'illumination 3,12 comprenant les bandes spectrales d'absorption de 15 différentes d'espèces gazeuses recherchées. Ainsi, le gaz analysé est éclairé par un rayonnement d'illumination 12 selon une bande spectrale d'illumination 3,12 dont la largeur spectrale est de préférence supérieure à 500 nm, voire à 1 p.m, voire à 3 p.m. Par largeur spectrale, on entend une différence entre la longueur d'onde la plus élevée et la longueur d'onde la plus faible de la bande spectrale d'illumination de la source. Le photodétecteur de mesure 20 détecte le 20 rayonnement 14 transmis par le gaz G dans la bande spectrale de détection Am dont la largeur spectrale est de préférence supérieure à 500 nm, voire à 1 p.m, voire à 3 p.m. Selon un mode de réalisation, contrairement au dispositif décrit dans W02007064370, le rayonnement lumineux 14 transmis par l'échantillon est détecté par le photodétecteur de mesure 20 dans une bande spectrale de détection Am s'étendant continûment entre une 25 longueur d'onde minimale À20-min et une longueur d'onde maximale À20 -max - La bande spectrale de détection Am comporte la bande spectrale d'absorption Ak d'au moins deux espèces gazeuses d'intérêt Gk. De préférence, la largeur de la bande spectrale de détection Am s'étend sans discontinuité selon au moins 500 nm, voire 1 p.m, voire encore davantage, par exemple selon un intervalle spectral supérieur à 2µm ou à 5 p.m. Ainsi, il n'est pas nécessaire 30 d'interposer un filtre multibandes présentant différentes bandes passantes, dont chaque bande passante est adaptée à une espèce gazeuse Gk à détecter. Selon ce mode de réalisation, on peut, avec le même dispositif, déterminer les quantités d'espèces gazeuses dans différents mélanges, 11 en adaptant le traitement, par le processeur 30, des intensités détectées, en fonction des espèces gazeuses Gk susceptibles d'être présentes dans le gaz. Par exemple, la bande spectrale de détection Am s'étend entre 3 p.m et 8 p.m. Le même dispositif 1 peut alors permettre une détermination des concentrations de CO2, CH4 et CO dans un mélange de gaz naturel, ou une 5 détermination des concentrations de NO et NO2 dans un mélange de type NOx. Les espèces gazeuses d'intérêt sont déterminées a priori, en fonction de l'application et/ou du mélange gazeux analysé. Le terme continûment signifie que le signal détecté par le photodétecteur de mesure 20 comporte toutes les longueurs d'onde entre la longueur d'onde minimaleminet la longueur 10 d'onde maximale 2. -20-max évoquées ci-dessus. Un aspect important de l'invention est que le gaz analysé G est soumis à une pluralité J d'illuminations successives, respectivement selon différents rayonnements d'illumination 12, chacun présentant un spectre d'illumination Si adressant les bandes spectrales d'absorption des espèces gazeuses d'intérêt. De préférence, les rayonnements d'illumination s'étendent selon 15 une même bande spectrale d'illumination 3,12 mais leurs spectres respectifs Si, dans cette bande spectrale, fluctuent entre chaque illumination. Sur la figure 2, on a représenté différents spectres Si dans une même bande spectrale d'illumination 3,12 s'étendant entre 1 p.m et 10 p.m. A chaque illumination selon un spectre Sj, le photodétecteur de mesure 20 détecte un rayonnement 14 transmis par le gaz, dans la bande spectrale de détection Am précédemment 20 décrite. Le photodétecteur de mesure 20 délivre un signal dont l'intensité est représentative de l'intensité transmise par le gaz lorsqu'il est illuminé par le rayonnement d'illumination de spectre Si. Cette intensité peut être estimée par l'expression suivante : = fA. Eray IO,j (À, ray)fiLlhk(ck, À, Xray)(5À =gray loi (À,ray)fiLle-12kMckxraY 8.1 (7) 25 est une estimation de l'intensité du rayonnement d'illumination 14 transmis par le gaz G dans toutes les longueurs d'onde À de la bande spectrale de détection Am lorsque le gaz comporte une concentration ck d'espèces gazeuses Gk. Le photodétecteur de mesure 20 n'étant pas spectralement résolu, l'intensité correspond à une intensité dans la bande spectrale de détection A20. 30 Après discrétisation selon la longueur d'onde, = EA. Eray /0,; (À, ray) Ilk e-i-Lik(A)CkXray 8À(8) 12 où 8.1 correspond au pas de discrétisation en longueur d'onde. Si K désigne le nombre d'espèces gazeuses d'intérêt Gk dont on souhaite déterminer la quantité ck dans le gaz G, on effectue J séquences d'illumination / détection, de façon à obtenir/ signaux de détection, d'intensité îj telle qu'estimée par l'expression (8). Chaque signal de détection est 5 détecté, de préférence, dans la même bande spectrale de détection A20. De préférence, le nombre/ de spectres d'émission Sj différents les uns des autres est supérieur ou égal au nombre K d'espèces gazeuses d'intérêt :J > K. Selon une possibilité, par exemple décrite dans la demande de brevet FR1658243, les intensités mesurées I sont traitées par le microprocesseur 30 relié à la mémoire 32, dans laquelle sont 10 stockées différentes grandeurs, notamment les fonctions spectrales d'absorption hk(ck,À, x) de chaque espèce gazeuse d'intérêt. La détermination de la quantité ck de chaque espèce gazeuse d'intérêt Gk peut alors être obtenue par une minimisation d'une somme d'écarts quadratiques respectivement calculés lors de chaque illumination j. Chaque écart quadratique correspond à une comparaison entre une mesure de l'intensité 4 et son estimation Î. Les quantités c1 cK 15 recherchées sont obtenues selon l'expression : ... cK) = argmin 1(4 - îj)2 2 = argmin Ej (i - EA.Eray /oi (À,ray)flikc=le-14Mckxray 8À) (9) Une telle approche est fiable, mais l'algorithme de minimisation suppose le recours à un microprocesseur doté d'une puissance de calcul élevé, en particulier du fait des discrétisations 20 selon À et selon xray nécessaires pour l'estimation de Î. Une telle puissance de calcul peut être difficile à embarquer, en particulier dans des dispositifs nomades. De plus, elle s'accompagne d'une consommation électrique importante. Les inventeurs proposent une alternative, visant à obtenir une estimation îj par une fonction paramétrique, nécessitant moins de temps de calcul que l'estimation décrite en lien avec 25 l'expression (8). En particulier, la fonction paramétrique ne suppose pas une discrétisation selon la longueur d'onde À et selon la distance xray parcourue par le rayonnement dans le gaz. L'approche proposée est basée sur l'utilisation d'une fonction paramétrique fejejk,njk, de telle sorte que îj = (C1 --- = -4,1(1 - (1 + 11k=1 e-ai'kCk c) (10) 13 où : Î0,1 est une estimation de l'intensité incidente au gaz G dans la bande spectrale de détection A20. Il s'agit également de l'intensité qui serait mesurée par le photodétecteur de mesure 20 en l'absence d'absorption dans le gaz. Comme préalablement indiqué, 5 cette intensité peut être obtenue à partir de la mesure de l'intensité de référence /f,'f, lors de l'illumination selon le spectre d'illumination Sj, l'intensité de référence /f,'f étant mesurée à l'aide du photodétecteur de référence 20'f. - Of est un paramètre dépendant du spectre d'illumination S. Il s'agit d'un réel. aj,k est un paramètre dépendant du spectre d'illumination Sf et de l'espèce gazeuse 10 d'intérêt Gk. - 71.1,k est un paramètre dépendant du spectre d'illumination Sf et de l'espèce gazeuse d'intérêt Gk. feraj,k,nj,k est une fonction paramétrique, dite de réponse, dont les paramètres sont Of, aj,k et 71.1,k. Le paramètre Of ne dépend que du spectre d'illumination Sf, donc de la température Tf 15 lorsque la source de lumière 11 est modulée en température. Les paramètres aj,k et 71.1,k dépendent de chaque espèce gazeuse Gk et du spectre d'illumination S. Il s'agit de paramètres spécifiques, pouvant varier d'une espèce gazeuse à une autre. Pour un spectre d'illumination Sf, les paramètres O, aj,k et 71.1,k sont des constantes dans la bande spectrale de détection 420. Leurs valeurs respectives ne varient pas en fonction de la longueur d'onde. Cela permet de 20 minimiser la quantité d'information à stocker dans la mémoire 32. La fonction paramétrique de réponse fejejk,'jk permet de modéliser, en prenant en compte un nombre limité de paramètres, l'absorption du rayonnement d'illumination 12 par le gaz G. Le recours à la fonction paramétrique de réponse fejejk,njk permet une estimation simple et rapide de l'intensité îf, sans nécessiter une puissance de calcul importante. La fonction 25 paramétrique de réponse fejejk,njk est une approximation de la fonction d'atténuation, de telle sorte que : îj = îo,j 1.8 i,ai,k,ni,k(Ci - -- CK) EA.Eray /0,; (À, ray)ee-4(À)ckxr" 8/1 (11). Ainsi, à partir des intensités mesurées pour différents spectres d'illumination Si, on obtient les quantités c1 cK des espèces gazeuses recherchées par une minimisation de type : 14 (ci ... cK) = argmin - îj)2 2 = argmin Ei - /0,ifejej,k,'i,k (ci ... cK)) (12) Dans l'expression (12), l'intensité mesurée /j est une grandeur représentative de l'absorption du rayonnement d'illumination 12 par le gaz G. D'une façon générale, les quantités c1 ... cK des 5 espèces gazeuses recherchées peuvent être obtenues en comparant, lors de chaque illumination j, une grandeur mesurée, représentative de l'absorption du rayonnement d'illumination 12 par le gaz G, avec une estimation de cette grandeur. Les quantités c1 cK sont celles pour lesquelles une combinaison arithmétique de chaque comparaison prend une valeur prédéterminée. Dans l'exemple correspondant à l'expression (12), la comparaison est un écart quadratique, et la 10 combinaison est une somme. La valeur prédéterminée est alors égale à 0. La comparaison peut être également un ratio. La combinaison arithmétique peut comprendre une somme, et/ou un produit, et/ou une division. Dans l'expression (12), la grandeur représentative de l'absorption du rayonnement d'illumination par le gaz est l'intensité mesurée par le photodétecteur de mesure 20 lors de chaque illumination. Selon une variante, cette grandeur peut être l'intensité 15 If mesurée par le photodétecteur de mesure 20 normalisée par l'intensité non absorbée précédemment décrite, ou toute autre grandeur exprimant l'absorption du rayonnement d'illumination 12 par le gaz G. La fonction paramétrique de réponse fo jk,n jk peut être apprise sur la base d'expressions théoriques de la loi de Beer-Lambert, selon les expressions (7) ou (8). Elle peut également être 20 apprise au cours d'une phase d'apprentissage, durant laquelle on introduit un gaz de calibration comportant des quantités ck* connues d'espèces gazeuses d'intérêt, dites quantités de calibration. La phase de calibration consiste, pour chaque spectre d'illumination Sj, à effectuer plusieurs mesures en modifiant, à chaque mesure, les quantités ck* des espèces gazeuses d'intérêt dans le gaz de calibration. La phase de calibration est mise en oeuvre en utilisant le 25 photodétecteur de mesure ou un photodétecteur de calibration considéré comme représentatif du photodétecteur de mesure. En fonction des intensités /j(ck* ) mesurées par le photodétecteur utilisé lors de la calibration, dites intensités de calibration, on ajuste au moins un des paramètres O, aj,k,71,j,k de la fonction de réponse, après avoir pris en compte la forme de la fonction paramétrique. Dans l'exemple présenté ci-dessus, la forme de la fonction de 30 réponse est : 15 1.8 CI{) = (1 - OI (1 1-11( -a - kc nj,k (13) k=1e k La figure 4 résume les principales étapes d'un procédé selon l'invention. Les étapes 90 à 96 concernent la phase de calibration, permettant d'établir les paramètres Oraj,k,71,j,k de la fonction de réponse fejejk,njk. 5 Étape 90 : introduction d'un gaz de calibration comportant des quantités de calibration ck* connues d'espèces gazeuses d'intérêt, par exemple CH4 et CO2 et illumination du gaz de calibration selon un spectre de calibration Étape 92 : détection de l'intensité /i(ck* ) de calibration, et estimation, à partir d'une mesure de référence li,ref, de l'intensité /0,i incidente au gaz de calibration, cette dernière correspondant 10 également à l'intensité qui serait détectée par le photodétecteur utilisé lors de la calibration en l'absence de gaz de calibration. Les étapes 90 à 92 sont réitérées en modifiant le spectre d'illumination Si. Étape 94 : modification des quantités de calibration ck* et réitération des étapes 90 à 92. Étape 96 : lors de l'étape 96, on dispose, pour chaque spectre d'illumination Si, de valeurs 15 d'intensités /i(ck* ) correspondant aux quantités de calibration ck* respectivement utilisées. L'étape 96 est la définition de la fonction de réponse fejejk,nik. On considère une forme paramétrique prédéterminée, par exemple: fei,ai,k,ni,k (cl* ". cx* )-(1- (1 +inick=1 et on détermine, pour chaque spectre d'illumination Sj, le jeu de paramètres (0j,aj,k, 711,0 20 présentant une meilleure attache aux données. Il s'agit par exemple de minimiser les erreurs ) entre et fo. a n (Ci* Cg* ). , .d j,k Les paramètres aj,k ou 71,j,k peuvent être déterminés, pour différents spectres d'illumination Si, en utilisant un photodétecteur de calibration représentatif du photodétecteur de mesure 20. Le paramètre O correspond à un facteur de recalage. Il est de préférence déterminé sur chaque 25 photodétecteur de mesure 20, car il peut varier significativement entre deux photodétecteurs de mesure différents. Ainsi, les paramètres apc ou 71,j,k peuvent être communs à différents dispositifs d'analyse 1, le paramètre O étant ajusté sur chaque dispositif d'analyse 1. 16 La figure 6A montre un exemple de détermination des paramètres d'une fonction de réponse par ajustement selon des points expérimentaux obtenus par modélisation. Cette figure est décrite ultérieurement, en lien avec des modélisations expérimentales réalisées. La fonction de réponse fejejk,n jk étant définie, on peut procéder à une estimation des quantités 5 d'espèces gazeuses ck d'un gaz inconnu G.Pour cela, le gaz à analyser est introduit dans l'enceinte 10 et on procède aux étapes suisvantes Étape 100 : illumination du gaz selon un spectre d'illumination Sj et mesure de l'intensité Ij du rayonnement 14 transmis par le gaz, ainsi que de l'intensité de référence /'f Étape 200: Réitération de l'étape 100 en modifiant le spectre d'illumination Sj 10 Étape 300 : à partir des intensités de mesure des intensités de référence obtenues lors de chaque étape 100, et de la fonction de réponse fejejk,njk correspondant à chaque spectre d'illumination Si, estimation des quantités d'espèces gazeuses d'intérêt ck en utilisant l'expression (12). L'obtention des quantités ck par minimisation d'une somme d'écarts quadratiques, selon 15 l'expression (12), suppose le recours à un algorithme de minimisation suffisamment rapide et peu coûteux en mémoire. Pour chaque spectre d'illumination Sj, les fonctions de réponse fe pajk,nj,k sont non linéaires. Les inventeurs ont considéré qu'un algorithme de type moindres carrés non linéaire était adapté, comme décrit ultérieurement en lien avec des essais modélisés. Les figures 5A à 5B montrent différentes configurations d'un dispositif d'analyse 1 permettant 20 une mise en oeuvre de l'invention. Sur la figure 5A, on utilise un unique photodétecteur 20, lequel peut être successivement associé à un filtre de mesure 18 ou à un filtre de référence 18ref tels que précédemment décrits. Une telle configuration nécessite un moyen mécanique pour effectuer une commutation entre le filtre de mesure 18 et le filtre de référence 18ref. Ce moyen peut par exemple être une roue à filtre. Sur la figure 5B, on a représenté une configuration en 25 réflexion, selon laquelle le rayonnement traverse le gaz jusqu'à une surface réfléchissante 19, dans la bande spectrale d'illumination Al2, avant d'atteindre le photodétecteur de mesure 20 et le photodétecteur de référence 20ref. Les inventeurs ont simulé une mise en oeuvre de la méthode de mesure précédemment décrite, en simulant une détection d'une quantité de CH4 et d'une quantité de CO2 dans un gaz 30 comportant une concentration progressivement croissante de CH4 et une concentration fixe de 17 400 ppm de CO2. Une telle simulation peut être représentative de la survenue d'une fuite de CH4 dans de l'air, au cours de laquelle la quantité de CH4 augmente progressivement. La figure 6A représente un exemple de détermination des paramètres Of, af,k,71,j,k d'une fonction de réponse fejejk,,,jk. Pour déterminer ces paramètres, on a simulé une détection 5 d'une intensité de mesure If et d'une intensité référence /ref,f en faisant varier la concentration ck* en CH4 (axe des abscisses), et en CO2 (axe des ordonnées). Pour réaliser les simulations, on a considéré une source de lumière 11, de type corps noir, émettant un spectre Sf correspondant à une température Tf de 450°C. I' La figure 6A représente une évolution du ratio - calculé en fonction de différents couples de /0,f 10 concentrations en CH4 et en CO2. A partir des points obtenus par calcul, on a procédé à un ajustement (fit) de façon à obtenir les paramètres Of, aj,k,71,j,k de telle sorte que la valeur de la fonction de réponse, pour chaque concentration en CO2 et en CH4, corresponde aux valeurs du ratio calculé. Les valeurs If ont été calculées application de la loi de Beer-Lambert, en utilisant /0,/ les fonctions spectrales d'absorption hk des espèces gazeuses considérées. La valeur Îod est 15 connue par la simulation de la source de lumière. Dans cet exemple, l'ajustement a été effectué: - pour CH4, en se basant sur une concentration nulle de CO2, c'est-à-dire en considérant les valeurs du ratio le long d'une ligne horizontale d'ordonnée 0. considérant rod - Pour CO2, en se basant sur une concentration nulle de CH4, c'est-à-dire en les valeurs du ratio le long d'une ligne verticale d'abscisse 0. rod 20 On peut ainsi estimer les paramètres Of, aik,71,j,k de la fonction de réponse correspondant au spectre d'illumination S. . Cette opération est renouvelée pour l'ensemble des spectres d'illumination S. La figure 6B représente, pour différents couples de concentrations en CO2 et en CH4, l'erreur relative entre les valeurs calculées et les valeurs estimées. Les valeurs Î ont été estimées 10,j Ioj 25 en utilisant la fonction de réponse fejejk,'jk. Les estimations ont été effectuées pour chaque couple de concentrations en CO2 et en CH4. Dans la plage de concentrations considérées, l'erreur relative maximale est inférieure à 1.5%. 18 Pour mettre en oeuvre la minimisation décrite en lien avec l'expression (12), différents algorithmes ont été testés, l'objectif étant d'identifier l'algorithme présentant le meilleur compromis erreur/temps de convergence. 5 algorithmes ont été comparés : - algorithme de Nelder-Mead ; 5 - algorithme de Newton-Gradient conjugué, connu sous l'acronyme "Newton CG" ; - algorithme Broyden-Fletcher-Goldfarb-Shanno, connu sous l'acronyme "BFGS" ; - algorithme de type "brute", selon lequel on dispose de valeurs de l'intensité îj préalablement établies, lors de la phase de calibration, pour différents couples de concentrations discrètes en CO2 et en CH4. Lors de la phase de mesure, le couple de 10 concentrations sélectionné est celui pour lequel la valeur de l'intensité préalablement établie est la plus proche de l'intensité correspondant au cas testé. - un algorithme de moindres carrés non linéaires "leastsq". La figure 6C représente une évolution de l'erreur relative ECO2 affectant la détermination de la concentration de CO2 (axe des ordonnées) en fonction de la concentration de CH4 (axe des 15 abscisses), en mettant en oeuvre les cinq algorithmes décrits ci-dessus. La concentration de CO2 est fixe, et correspond à la concentration dans l'atmosphère, soit 400 ppm. La concentration de CH4, repérée sur l'axe des abscisses, représente une augmentation progressive de la concentration en CH4 résultant de la fuite. L'algorithme de Nelder-Mead conduit à des résultats moins satisfaisants que les quatre autres, dont les performances sont similaires et se 20 superposent à la ligne continue traçant la performance de l'algorithme des moindres carrés non linéaires. La figure 6D représente une évolution de l'erreur relative ECH4affectant la détermination de la concentration de CH4 (axe des ordonnées) en fonction de la concentration de CH4 (axe des abscisses), en mettant en oeuvre les 5 algorithmes décrits ci-dessus. On en tire les mêmes 25 conclusions qu'en lien avec la figure 6C. La figure 6E représente une évolution de la durée de convergence de l'algorithme de minimisation (axe des ordonnées), en fonction de la concentration en CH4 (axe des abscisses). On observe que l'algorithme des moindres carrés non linéaires mène à des temps de convergence plus faibles que les autres algorithmes testés. 30 On a simulé une détermination d'une concentration de CO2 et de CH4 en mettant en oeuvre un algorithme de minimisation de type moindres carrés non linéaires, dans un mélange gazeux dans 19 lequel la concentration de CO2 est de 400 ppm, la concentration de CH4 variant entre 0 et 25000 ppm. 6 configurations ont été testées : - 8 paliers de température T1...T8, variant entre 300 et 720 °C, régulièrement espacés par des incréments de 60°C (configuration a); 5 - 8 paliers températures T1...T8, variant entre 450 et 730 °C, régulièrement espacés par des incréments de 40° (configuration b); - 8 paliers températures T1...T8, variant entre 600 et 740 °C, régulièrement espacés par des incréments de 20° (configuration c); - 4 paliers températures T1...T4, variant entre 300 et 750 °C, régulièrement espacés par 10 des incréments de 150° (configuration d); - 4 paliers températures T1...T4, variant entre 450 et 750 °C, régulièrement espacés par des incréments de 100° (configuration e); - 4 paliers températures T1...T4, variant entre 600 et 750 °C, régulièrement espacés par des incréments de 50° (configuration f). 15 Les figures 6F et 6G montrent respectivement les évolutions des erreurs relatives ECO2 et ECH4 de quantification du CO2 et du CH4 selon les différentes configurations. On observe que l'invention est exploitable en combinant quatre températures, ce qui correspond à quatre spectres d'illumination différents. Les meilleures performances correspondent aux configurations combinant un nombre élevé de paliers de température et une importante plage 20 de variation de la température. Cependant, plus le nombre de paliers en température augmente, plus le temps de mesure augmente. Un compromis doit donc être défini, REVENDICATIONS1. Procédé d'analyse d'un gaz (G), pour estimer une quantité (ck) de plusieurs espèces gazeuses (Gk) différentes, dites espèces gazeuses d'intérêt, potentiellement présentes dans le gaz analysé, chaque espèce gazeuse étant apte engendrer une absorption d'un rayonnement lumineux dans une bande spectrale d'absorption (Ak), le procédé comportant les étapes suivantes : a) disposition du gaz analysé entre une source de lumière (11) et un photodétecteur de mesure (20), la source de lumière (11) étant apte à émettre un rayonnement d'illumination (12) se propageant à travers le gaz, le photodétecteur de mesure (20) étant apte à détecter un rayonnement (14) transmis par le gaz et à délivrer un signal de détection représentatif d'une intensité (Ii) du rayonnement détecté ; b) illumination du gaz analysé (G) à l'aide de la source de lumière (11), le rayonnement d'illumination (12) s'étendant selon une bande spectrale, dite bande spectrale d'illumination (Al2), recouvrant la bande spectrale d'absorption (Ak) d'une pluralité d'espèces gazeuses d'intérêt, à chaque illumination correspondant un spectre (S1), dit spectre d'illumination, du rayonnement d'illumination (12) ; c) lors de l'illumination, détection, par le photodétecteur de mesure (20), d'une intensité (Ii) d'un rayonnement transmis (14) par le gaz analysé, dans au moins une bande spectrale de détection (Am) d) réitération des étapes b) et c), en modifiant le spectre d'illumination (Si) du rayonnement d'illumination (12); e) détermination d'une quantité (ck) de chacune des espèces gazeuses d'intérêt en fonction de l'intensité (I]) du rayonnement détecté lors de chaque illumination, chaque intensité (I]) étant associée à un spectre d'illumination (Si); le procédé étant caractérisé en ce que l'étape e) comporte : une prise en compte, pour chaque spectre d'illumination, d'une fonction de réponse paramétrique, la fonction de réponse paramétrique (fej,'Rmik) étant définie par des paramètres (Oi, ai,k,ni,k), de telle sorte qu'à chaque spectre d'illumination sont associés des paramètres (Opai,k,ni,k), les paramètres associés à deux spectres d'illumination distincts (Si, Si_i) étant différents et préalablement définis ; 21 - une estimation d'une grandeur (/ î) représentative de l'absorption du gaz (G), lors de P chaque illumination, à l'aide de la fonction paramétrique (fej,'Rmik). 2. Procédé selon la revendication 1, dans lequel à chaque spectre d'illumination sont associés des paramètres dits spécifiques, (aj,k,71,j,k) définis pour chaque espèce gazeuse d'intérêt (ck), 5 au moins deux paramètres spécifiques respectivement définis pour deux espèces gazeuses d'intérêt distinctes étant différents les uns des autres. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape e) comporte une estimation, lors de chaque illumination, d'une intensité (4,j), dite intensité non atténuée, qui serait mesurée par le photodétecteur de mesure (20) en l'absence de gaz (G). 10 4. Procédé selon la revendication 3, dans lequel lors de l'étape e), la grandeur représentative de l'absorption du gaz (G) est, pour chaque illumination, l'intensité (Ii) mesurée lors de l'étape c), éventuellement normalisée par l'estimation de l'intensité non atténuée (/0,]). 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape e) comporte, pour chaque illumination, une comparaison entre : 15 - la grandeur (/_/) représentative de l'absorption du gaz estimée, pour chaque illumination, lors de l'étape e) ; - et une mesure (Ii) de ladite grandeur, à l'aide de l'intensité mesurée lors de l'étape c) ; les quantités déterminées (ck) étant obtenues en fonction de chaque comparaison. 6. Procédé selon la revendication 5 dans lequel chaque comparaison est un écart quadratique 20 ((îf - /1)2), les quantités déterminées étant celles pour lesquelles une combinaison arithmétique de chaque comparaison se rapproche d'une valeur prédéterminée. 7. Procédé selon la revendication 6, dans lequel la combinaison arithmétique est une somme, les quantités déterminées (ck) étant celles pour lesquelles la somme est minimale. 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape e) met 25 en oeuvre une méthode par moindres carrés non linéaire. 9. Procédé selon l'une quelconque des revendications précédentes, comportant également, lors de chaque illumination, une détection, par un photodétecteur de référence (20ref) d'une intensité de référence (/],'f) d'un rayonnement de référence (14ref), considéré comme n'ayant pas été absorbé par le gaz (G), l'étape e) comportant une prise en compte des intensités de 30 référence (/],'f) respectivement détectées au cours de chaque illumination. 22 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel chaque spectre d'illumination (Si) correspond à une température (Ti) de la source de lumière. 11. Procédé selon l'une quelconque des revendications précédentes, comportant, préalablement à l'étape e), une phase de calibration, effectuée pour chaque spectre 5 d'illumination (Si), de façon à définir la fonction de réponse paramétrique (fej,,,jk,'jk) associée au spectre d'illumination (Si), la phase de calibration comportant les étapes suivantes: cal-i) réalisation des étapes a) à c) en considérant un gaz de calibration, comportant une quantité connue d'au moins une espèce gazeuse d'intérêt (ck*), le gaz de calibration étant illuminé selon un spectre d'illumination (Si), de façon à obtenir une intensité 10 détectée par le photodétecteur de mesure (20) ou par un photodétecteur de calibration, considéré comme représentatif du photodétecteur de mesure; cal-ii) réitération de l'étape cal-i) en modifiant la quantité (ck*) de l'espèce gazeuse d'intérêt ; cal-iii) prise en compte d'un modèle paramétrique de la fonction de réponse (fejejkm,jk); 15 cal-iv) à partir des intensités résultant de chaque étape cal-i) et du modèle paramétrique prise en compte lors de l'étape cal-iii), détermination d'au moins un paramètre aj,k,nj,k) de la fonction de réponse. 12. Procédé selon l'une quelconque des revendications précédentes, dans lequel la bande spectrale de détection (420) recouvre des bandes spectrales d'absorption de plusieurs espèces 20 gazeuses d'intérêt (4k). 13. Dispositif (1) d'analyse d'un gaz (G) comportant : - une source de lumière (11), apte à émettre un rayonnement d'illumination (12), la source de lumière étant également apte à être modulée de telle sorte que le spectre du rayonnement d'illumination émis (Si), dit spectre d'illumination, est variable dans une 25 même bande spectrale d'illumination (Al2) - un photodétecteur de mesure (20), apte à détecter un rayonnement lumineux selon une bande spectrale de détection (Am); - une enceinte (10), apte à recevoir le gaz, agencée de telle sorte que le photodétecteur de mesure (20) soit configuré pour détecter un rayonnement lumineux (14) ayant 30 traversé le gaz ; - un processeur (30), configuré pour mettre en oeuvre l'étape e) d'un procédé selon l'une quelconque des revendications 1 à 12, à partir du rayonnement détecté par le 23 photodétecteur de mesure (20) suite à différentes illuminations du gaz (G), respectivement selon différents spectres d'illumination (Si).