FR 3042038 A1 20170407 FR 1559328 A 20151001 1 PROCEDE D'OPTIMISATION DE LONGUEURS D'ONDE DE DETECTION POUR UNE DETECTION MULTI-GAZ DESCRIPTION DOMAINE TECHNIQUE ET ART ANTÉRIEUR L'invention concerne le domaine de la spectroscopie optique, et plus particulièrement celui de la détection d'espèces chimiques, en phase gazeuse ou liquide, par absorption spectrale. Un capteur de gaz optique, et plus particulièrement un capteur de type Non Dispersif Infra-Rouge (NDIR), permet de calculer la concentration d'un gaz en phase liquide ou vapeur. Il comporte une source infrarouge dont le rayonnement est dirigé à travers une cellule dans laquelle se trouve le gaz à détecter vers un détecteur qui est pourvu généralement d'un filtre spectral. La figure 1 représente schématiquement un tel capteur 10 comportant une source lumineuse infrarouge 12 disposée dans une cellule 14 munie d'une entrée 16 et d'une sortie 18 par lesquelles le gaz entre et sort de la cellule 14. Un rayonnement lumineux infrarouge émis par la source 12 est dirigé vers un détecteur 20 et passe à travers un filtre spectral 22 avant d'atteindre le détecteur 20. En variante, la source 12 peut être une source lumineuse laser et le capteur 10 ne comporte dans ce cas pas de filtre spectral en raison du caractère monochromatique de la lumière émise par une telle 2 0 source. La plupart des gaz possédant chacun une signature distincte dans le domaine moyen infrarouge (par exemple les longueurs d'ondes comprises entre environ 3 um et 15 um), une partie du rayonnement infrarouge émis par la source 12 dans le domaine MIR est absorbée par le gaz présent dans la cellule 14 et la baisse de l'intensité 25 du signal lumineux qui résulte de cette absorption est quantifiée par comparaison avec une voie de référence, non sujette à cette variation. Ce phénomène est décrit par la loi de BeerLambert, qui stipule que : /(À.) _ E (À.)/,c T(À) - lu (1) 3042038 2 avec Â la longueur d'onde considérée, T(Â) la transmittance à travers le gaz (c'est-à-dire sans prendre en compte le filtrage réalisé par le filtre spectral 22), I(X) l'intensité spectrale de sortie (mesurée par le détecteur), lo(X) l'intensité spectrale d'entrée (émise par la source), E(X) le coefficient d'absorption spectrale du gaz, L la longueur de 5 cavité traversée par le rayonnement infrarouge et C la concentration du gaz dans la cavité. La longueur L de la cavité est en général connue, tout comme le coefficient d'absorption spectrale E(X) du gaz qui peut être mesuré via une calibration préalable. L'intensité spectrale d'entrée 10(X) est plus difficile à mesurer et est souvent obtenue à partir d'une voie de référence. Cette voie de référence consiste généralement 10 en un second capteur pourvu d'un filtre optique situé hors de l'absorption du gaz à mesurer, par exemple disposé dans une cellule différente de celle dans laquelle le gaz est injecté, et comprenant une source émettant un faisceau lumineux à une longueur d'onde XR. Cette voie de référence est par exemple formée d'éléments similaires aux éléments 12, 14, 20 et 22 du capteur 10, la différence étant par exemple qu'aucun gaz n'est présent dans 15 la cellule de cette voie de référence. La transmittance T(X) de l'équation (1) est alors obtenue grâce à l'équation : 7, (À) = pRo (À.R) ( I A(À) _ 10- EMLC I Ao(À) J \I R(4)/ - avec IR (1R) absorption par le gaz dans cette voie de référence, et donc IA(À) = 410(À)10-EWLc donné par l'équation (1). Ici, l'indice R indique la voie de référence, l'indice A la voie active et l'indice 0 une mesure à concentration nulle (au niveau de l'entrée de la cellule). La calibration du rapport IRO(LR)/IAO(L) et la mesure du rapport /A(À)//R(1R) permet de calculer la transmittance T(X) à travers le gaz. En prenant le logarithme de cette transmittance T(X), il est possible de remonter à la concentration du gaz C : log T (À.) C = (3) EML Lorsque l'absorption du gaz est faible, il est possible de linéariser l'équation (1) et la concentration C est alors approximativement proportionnelle au facteur (1-T). (2) = /R0(4) puisque le faisceau lumineux ne subit aucune 3042038 3 Dans le cas où deux gaz différents sont destinés à être mesurés et discriminés l'un de l'autre, il est possible d'utiliser deux voies actives indépendantes, chacune possédant un filtre spectral positionné sur le spectre d'absorption de l'un des premier et second gaz (ou possédant chacune une source lumineuse laser émettant des 5 faisceaux lumineux à des longueurs d'onde différentes centrées vis-à-vis des spectres d'absorption des gaz). Chaque voie permet alors d'obtenir une information sur la concentration du gaz pour laquelle elle a été conçue. Cette approche n'est cependant valide que lorsqu'il n'y a pas de recouvrement entre les spectres d'absorption des deux gaz car sinon, il peut y avoir une 10 confusion entre les deux gaz détectés à une certaine longueur d'onde, et donc une erreur possible sur le résultat obtenu. Afin de résoudre cette indétermination, il est nécessaire de croiser les réponses des deux capteurs. La loi de Beer-Lambert selon l'équation (1) ci-dessus s'écrit alors pour deux gaz : 4)c2)L IT(À1) = 10- (E1(4)ci+E2( T(2.2) = (Ei (4)ci +E2(4)c2)L 15 avec Al correspondant à la longueur d'onde à laquelle la première mesure est réalisée (c'est-à-dire la longueur d'onde sur laquelle est positionné un premier filtre spectral du capteur dont la bande passante est de l'ordre de la largeur de la raie d'absorption du premier gaz, ou la longueur d'onde émise par une première source lumineuse laser du capteur), et Â2 correspondant à la longueur d'onde à laquelle la 20 deuxième mesure est réalisée (c'est-à-dire la longueur d'onde sur laquelle est positionné un deuxième filtre spectral du capteur dont la bande passante est de l'ordre de la largeur de la raie d'absorption du deuxième gaz, ou la longueur d'onde émise par une deuxième source lumineuse laser du capteur). Des équations (4), on en déduit deux équations à deux inconnues qui peuvent être résolues, pour déterminer les concentrations C1 et C2 des deux 25 gaz : c2 = L E2(22) 109 T(4)-E2 (4) 109 T(À-2) (5) { (À-1) E2 (À-2)-El (À-2) E2 (À-1) E1(4) 109 T(À-2)-E1 (À-2) 109 T(À-1) (À-1) E2 (À-2)-El (2L2) E2(À-1) (4) 3042038 4 Dans le cas où e1(X2) = E2(M) = 0 (pas de recouvrement entre les spectres d'absorption des deux gaz aux longueurs d'onde Al et Â2), les équations ci-dessus se simplifient en deux équations (3) indépendantes. Dans le cas contraire, une confusion entre les gaz détectés est possible. 5 Le cas ci-dessus peut s'appliquer à la détection de n gaz, avec n > 2. Un système de n équations à n inconnues, qui peut être résolu via diverses techniques, est alors obtenu. Comme indiqué dans le document de C.W. Brown et al., « Matrix representations and criteria for selecting analytical wavelengths for multicomponent spectroscopic analysis », Anal. Chem. 1982, 54(9), pp. 1472-1479, un tel système peut être 10 posé de façon matricielle, les équations (1) et (4) pour n gaz pouvant alors s'écrire : log 7T (À1)\ T(Â2) /£11 1) E1(12) E2 (11) - -- En(21)\ En(Â2) 7 LCi\ LC2 (6) \T (in)/ - E2(12) - - - \E1(Ân) - - En(4)/1 \LCn/I E2(Ân) Cette équation (6) peut s'écrire de façon plus condensée : log(T) = -(E)(LC) (7) Dans l'équation (7) ci-dessus, T et LC correspondent à des vecteurs et E 15 correspond à une matrice appelée matrice d'absorptivité. Au vu de l'équation (7), il apparait qu'obtenir le vecteur LC revient à une simple inversion de matrice, permettant donc d'obtenir finalement : (LC) = -(E-1) log(T) (8) En théorie, il suffit que la matrice E soit inversible pour garantir l'unicité 20 du vecteur de solution LC. Autrement dit, le passage de l'équation (7) à l'équation (8) est possible si les valeurs des colonnes de la matrice E sont indépendantes, c'est-à-dire s'il n'existe pas de combinaison linéaire pour passer des valeurs des coefficients d'absorption spectrale d'une colonne (relatives à l'un des gaz) aux valeurs des coefficients d'absorption spectrale des autres colonnes (relatives aux autres gaz) de la matrice E. Cependant, le 25 positionnement des mesures réalisées pour obtenir les valeurs du vecteur T, c'est-à-dire les valeurs )t,, influe fortement sur la précision du résultat obtenu. En effet, une erreur due à un bruit de mesure se propage sur le résultat qui correspond à la valeur de concentration affichée. Ainsi, en considérant le cas de deux mesures réalisées pour deux gaz aux 3042038 5 longueurs d'onde Al et X2, plus les ratios E1(À1)/E2(À1) et E2(À.2)/E1(À.2) sont élevés, plus la mesure est fiable. Dans les documents US 5 041 723, US 4 914 719, US 5 281 817 et US 6 307 201, différentes méthodes sont proposées pour choisir les valeurs des longueurs 5 d'onde ?t, sur lesquelles les filtres du capteur sont positionnés ou correspondant aux longueurs émises par les sources lumineuses laser du capteur. Ces approches sont toutefois relativement empiriques et ne peuvent pas être généralisées à n'importe quel nombre et type de gaz à détecter. EXPOSÉ DE L'INVENTION 10 Un but de la présente invention est de proposer un procédé permettant de déterminer des valeurs optimales de longueurs d'onde de détection, c'est-à-dire les longueurs d'onde sur lesquelles des filtres spectraux du capteur sont destinés à être positionnés ou les longueurs d'onde de faisceaux lumineux destinés à être émis par des sources lumineuses laser du capteur), pour réduire au maximum l'incidence des bruits de 15 mesure sur le résultat obtenu, et cela quel que soit le nombre n de gaz à détecter et quels que soient leurs spectres d'absorption, que ces spectres se recouvrent ou non. Pour cela, la présente invention propose un procédé d'optimisation de valeurs de n longueurs d'onde de détection d'un capteur de gaz optique apte à réaliser une détection de n gaz différents, comportant la mise en oeuvre des étapes suivantes : 20 a) calcul d'une valeur d'un déterminant d'une matrice d'absorptivité E dont des coefficients sont représentatifs de l'absorptivité spectrale de chacun des n gaz aux n longueurs d'onde de détection, l'étape a) étant répétée plusieurs fois en modifiant à chaque fois au moins une des valeurs desdites n longueurs d'onde de détection et telle que les valeurs desdites n longueurs d'onde de détection soient comprises dans une gamme de 25 valeurs pour laquelle l'absorptivité spectrale d'au moins un des n gaz est non nulle ; b) détermination des valeurs desdites n longueurs d'onde de détection pour lesquelles la valeur du déterminant de la matrice d'absorptivité E correspond à la valeur maximale parmi l'ensemble des valeurs calculées lors de l'étape a). 3042038 6 Ainsi, l'invention propose un procédé permettant, à l'aide d'un calcul d'erreur utilisant judicieusement la matrice d'absorptivité E, de choisir de façon optimale les longueurs d'onde de détection auxquelles les mesures sont destinées à être réalisées par le capteur afin de discriminer la présence de différents gaz lors d'une mesure 5 spectroscopique par absorption. Contrairement aux procédés de l'art antérieur se basant sur des méthodes essentiellement empiriques, le procédé selon l'invention peut s'appliquer à n'importe quel nombre de gaz et à n'importe quel type de gaz. De plus, ce procédé peut s'appliquer autant à un capteur de gaz optique de type NDIR ou multi-spectral. 10 Dans une première configuration, le procédé peut être tel que : - le capteur de gaz optique comporte au moins n sources lumineuses laser, - les coefficients de la matrice d'absorptivité E correspondent, ou sont uniformément proportionnels, aux coefficients d'absorption spectrale de chacun des n gaz 15 aux n longueurs d'onde de détection, et - les valeurs des n longueurs d'onde de détection déterminées à l'étape b) correspondent à des longueurs d'onde de faisceaux lumineux aptes à être émis par les n sources lumineuses laser. Dans une deuxième configuration, le procédé peut être tel que : 20 - le capteur de gaz optique comporte au moins une source lumineuse polychromatique et au moins n filtres spectraux chacun apte à permettre une détection d'un des n gaz, - les coefficients de la matrice d'absorptivité E correspondent, ou sont uniformément proportionnels, aux coefficients d'absorption spectrale moyens de chacun 25 des n gaz, pondérés ou non par des réponses spectrales des n filtres spectraux, et - les valeurs des n longueurs d'onde de détection déterminées à l'étape b) correspondent à des longueurs d'onde sur lesquelles les réponses spectrales des n filtres spectraux sont positionnées. Dans ce cas, les coefficients, notés É''./ (AI), de la matrice d'absorptivité E sont tels que : 3042038 7 f si MA i S (À) clÀ si (À i) = f Ai(À)S(À) dÀ. ou tels que : Si i) = f (À.) AiMS (À) clÀ avec ?t, correspondant aux n longueurs d'onde de détection, i étant tel 5 que 1 n ; EJ(X) correspondant aux coefficients d'absorption spectrale des n gaz, j étant tel que 1 j n ; A,(X) correspondant aux réponses spectrales des n filtres spectraux ; S(À) correspondant au spectre d'émission de la source lumineuse. 10 Dans ce cas, le procédé peut être tel que : - les n filtres spectraux sont choisis parmi un ensemble de m filtres spectraux, avec m > n dont un gabarit, une valeur de longueur d'onde sur laquelle la réponse spectrale est positionnée et une largeur de la réponse spectrale de chacun des filtres spectraux dudit ensemble de filtres sont prédéterminés, et 15 - les valeurs des n longueurs d'onde de détection sont choisies parmi les valeurs des longueurs d'onde sur lesquelles les réponses spectrales des filtres spectraux dudit ensemble de filtres sont positionnées. Ainsi, ce procédé peut être appliqué pour choisir, parmi plusieurs filtres prédéfinis, ceux les plus adaptés pour réduire les erreurs suivant le nombre et la nature des 20 gaz à détecter. En variante, le procédé peut être tel que : - les n filtres spectraux ont un même gabarit AM, et - les coefficients, notés (AI), de la matrice d'absorptivité E sont tels f AiMsi MS (À) clÀ si (À i) = f AiMS (À) dÀ. ou tels que : que : 25 3042038 8 èj(Ài) = f AiMsiMS (2)cIÀ avec ?t, correspondant aux n longueurs d'onde de détection, i étant tel que 1 i n ; EJ(X) correspondant aux coefficients d'absorption spectrale des n gaz, 5 j étant tel que 1 j n; A,(X) correspondant aux réponses spectrales des n filtres spectraux positionnées sur les longueurs d'onde A;; S(À) le spectre d'émission de la source lumineuse. Dans cette variante, le procédé sert à choisir, en utilisant pour la 10 détection de chacun des n gaz un filtre spectral d'un même gabarit et de même largeur de réponse spectrale, les positions des filtres les plus adaptées, c'est-à-dire les longueurs d'onde sur lesquelles les réponses spectrales des filtres spectraux sont positionnées, pour réduire les erreurs suivant le nombre et la nature des gaz à détecter. Dans ce cas, le procédé peut être tel que : 15 - les largeurs des réponses spectrales des n filtres spectraux sont ajustables, - l'itération de l'étape a) est répétée plusieurs fois en modifiant à chaque fois la largeur des réponses spectrales des n filtres spectraux, et - l'étape b) détermine les valeurs des n longueurs d'onde de détection 20 pour lesquelles la valeur du déterminant de la matrice d'absorptivité E correspond à la valeur maximale parmi l'ensemble des valeurs calculées lors des multiples itérations de l'étape a). Ainsi, le procédé permet d'optimiser également la largeur des réponses spectrales des n filtres spectraux utilisés pour réduire les erreurs en fonction du nombre et 25 de la nature des gaz à détecter. Le procédé peut être tel que : - le capteur de gaz optique ne comporte pas de voie de référence et comporte n+1 voies de mesure ; 3042038 9 - les coefficients de la matrice d'absorptivité E correspondent, ou sont uniformément proportionnels, aux coefficients d'absorption spectrale, moyens ou non, de chacun des n gaz, pondérés ou non par des réponses spectrales des n filtres spectraux, auxquels sont soustrait les coefficients d'absorption spectrale d'une des n+1 voies de 5 mesure du capteur de gaz optique servant de voie de référence. Ainsi, lorsque le capteur de gaz optique ne comporte pas de voie de référence permettant de réaliser une mesure d'absorption spectrale sans la présence de gaz dans cette voie, le procédé peut être tout de même mis en oeuvre en utilisant une des voies de mesure qui sert alors de voie de référence. La voie servant de référence peut 10 correspondre à n'importe laquelle des voies de mesure du capteur. Cela revient à choisir de façon optimale n+1 longueurs d'onde de détection de façon à maximiser le déterminant de la matrice d'absorptivité E. Les étapes a) et b) peuvent être mises en oeuvre par exemple via une méthode d'optimisation de type solveur d'équations par algorithme évolutionnaire. 15 BRÈVE DESCRIPTION DES DESSINS La présente invention sera mieux comprise à la lecture de la description d'exemples de réalisation donnés à titre purement indicatif et nullement limitatif en faisant référence aux dessins annexés sur lesquels : - la figure 1 représente un capteur de gaz optique de type NDIR ; 20 - la figure 2 représente les spectres d'absorption, ou coefficients d'absorption spectrale, de trois gaz destinés à être détectés par un capteur optimisé par la mise en oeuvre d'un procédé d'optimisation des valeurs des longueurs d'onde de détection du capteur, objet de la présente invention, selon mode de réalisation particulier ; - la figure 3 représente un gabarit de filtre du capteur optimisé par la 25 mise en oeuvre du procédé d'optimisation des valeurs des longueurs d'onde de détection, objet de la présente invention, selon le mode de réalisation particulier ; - la figure 4 représente les produits de convolution du gabarit du filtre représenté sur la figure 3 avec les spectres d'absorption représentés sur la figure 2 ; 3042038 10 - la figure 5 représente des longueurs d'onde de détection calculées lors de la mise en oeuvre du procédé selon le mode de réalisation particulier, et superposées aux produits de convolution de la figure 4 ; - la figure 6 représente des longueurs d'onde de détection calculées lors 5 de la mise en oeuvre du procédé selon le mode de réalisation particulier et superposées aux spectres d'absorption de la figure 2. Des parties identiques, similaires ou équivalentes des différentes figures décrites ci-après portent les mêmes références numériques de façon à faciliter le passage d'une figure à l'autre. 10 Les différentes parties représentées sur les figures ne le sont pas nécessairement selon une échelle uniforme, pour rendre les figures plus lisibles. Les différentes possibilités (variantes et modes de réalisation) doivent être comprises comme n'étant pas exclusives les unes des autres et peuvent se combiner entre elles. 15 EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Lors d'une mesure de concentrations de gaz par le capteur de gaz optique 10, telle que celle décrite précédemment, les principales sources d'erreurs proviennent de la source de rayonnement 12 et du détecteur 20 utilisés. En considérant que la source 12 émet une intensité spectrale I + alsour', avec alsour' correspondant au bruit de la source 20 12, et que le détecteur 20 rajoute un bruit 6kapteur' alors un signal I + 2 2 2 avec al = a/source + cricapteur est obtenu au final. Selon les cas, le terme du bruit issu de la source 12 ou du détecteur 20 prédomine, mais il est difficile de les dissocier. En pratique, il est plus facile de mesurer directement al comme l'écart type du signal mesuré par le capteur 10 en conditions fixes. En calculant la transmittance T = 1/10, il est possible de 25 définir un bruit GT tel que : (12 = (12 + (0-10)2 (9) Dans le cas d'un capteur 10 apte à réaliser la détection de n gaz, n transmittances T, sont mesurées, chacune avec une imprécision o-T, E [1 ...n]) donnée. 3042038 11 Ces valeurs se retrouvent sur la diagonale de la matrice de covariance ET associée, les valeurs hors diagonales étant nulles : /011 0 0 - ,TT2 u 0 2 ,T "Tn/ (ET) = 0 0 (10) L'erreur sur les transmittances créé une erreur sur les produits LC en se 5 propageant à travers l'équation (8) et mène à la formule suivante pour la matrice de covariance de LC : (ILc) = (E-1)(h)(IT) (MT (E-1)T Où h est la matrice jacobienne de la fonction f (fi (Ti, ... Tn) (-log Ti fn (Ti, ... Tn) -log Tn f (Ti, ... Tn) = définie par : (12) 10 La matrice jacobienne h de cette fonction f peut donc s'écrire selon l'équation : /1/T1 0 - - - 0 \ 1 0 1/T2 ... 0 0 0 - - - - 1/Tn/ If = In 10 (13) De plus, dans l'équation (11) ci-dessus, l'inverse de la matrice E peut s'exprimer en fonction de son déterminant et de sa comatrice par l'équation : 15 E-1 = 1 tCOM(E) (14) det(E) Les formules générales ci-dessus sont valides quel que soit le nombre n de gaz à détecter et peuvent être complétées s'il existe d'autres sources d'erreur. Plus spécifiquement, il est possible d'écrire la version étendue des variances de l'équation (11) dans le cas d'un capteur apte à réaliser la détection de deux gaz : 20 {CILCi2 = ( 2 _ aL C2 1 2 (ln 10( .EilE22-E21E12) [ T1 , ) 2 -21-Ti + ( T2 EliaT2) 1 2 2 [(E220-Ti) + (Eizer2)2 ln 10( .EilE22-E21E12) T1 T2 2 (15) ] 3042038 12 En utilisant les expressions de l'équation (11), il est possible d'évaluer précisément l'erreur sur toutes les concentrations à mesurer en fonction des valeurs de la matrice E correspondantes aux différentes raies d'absorption des gaz et des concentrations des différents gaz. Il est donc possible de sélectionner n longueurs d'onde autour 5 desquelles les n filtres doivent être positionnés (à partir de la matrice E correspondante) de façon à respecter un cahier des charges bien défini : toutes les erreurs en dessous d'une valeur seuil, une précision différente pour chaque gaz, etc. Du fait qu'un filtre correspond à un gabarit qui peut être translaté sur une plage de longueurs d'onde, le positionnement des filtres correspond au positionnement de son gabarit par rapport à une longueur d'onde 10 de référence. Cette longueur d'onde de référence peut changer selon la forme du gabarit du filtre. Par exemple, le positionnement d'un filtre gaussien, ou tout autre filtre dont le gabarit est symétrique par rapport à une longueur d'onde centrale, peut se faire en choisissant comme longueur d'onde de référence la longueur d'onde centrale de la réponse spectrale du filtre. Pour un filtre de type passe-haut ou passe-bas, la longueur d'onde de 15 référence peut correspondre par exemple à la longueur d'onde de coupure filtre. Cependant, dans le cas général, il est difficile d'optimiser le système avec les n valeurs diagonales de la matrice (ELS). Par exemple, dans les variances exprimées à l'équation (15), celles-ci comportent deux termes distincts : le premier est proportionnel à 2 l'inverse du déterminant au carré (correspondant aux termes ( \) 1 pour ln 10( .EllE22-E21E12) 20 GLC 2 1 et 6ic2 2)' et le second tient compte des erreurs en l'absence de gaz sur la transmittance et contient des cofacteurs de la matrice E au carré (correspondant aux 2 [(TE22aTi) 2 + (E12eT2)21 pour °ici2 et (E2leT1)2 + ( EllTeT2) termes 1 pour oic2 2). 1 T2 T1 2 Pour diminuer l'ensemble des erreurs, il est proposé de maximiser le déterminant dans les premiers termes mentionnés ci-dessus. Le déterminant de la matrice 25 d'absorptivité E est donc utilisé comme indicateur de performance pour la sélection des longueurs d'onde à utiliser pour la détection des gaz. Plus celui-ci est élevé en valeur absolue, plus l'erreur générale sur la détermination des concentrations est faible. 3042038 13 De manière générale, le procédé d'optimisation de valeurs de n longueurs d'onde de détection d'un capteur de gaz optique apte à réaliser une détection de n gaz différents comporte la mise en oeuvre des étapes suivantes : a) calcul de la valeur du déterminant de la matrice d'absorptivité E dont 5 des coefficients sont représentatifs de l'absorptivité spectrale de chacun des n gaz aux n longueurs d'onde de détection, l'étape a) étant répétée plusieurs fois en modifiant à chaque fois au moins une des valeurs desdites n longueurs d'onde de détection et telle que les valeurs desdites n longueurs d'onde de détection soient comprises dans une gamme de valeurs pour laquelle l'absorptivité spectrale d'au moins un des n gaz est non nulle ; 10 b) détermination des valeurs desdites n longueurs d'onde de détection pour lesquelles la valeur du déterminant de la matrice d'absorptivité E correspond à la valeur maximale parmi l'ensemble des valeurs calculées lors de l'étape a). L'ensemble des longueurs d'onde de détection potentielles sont donc évaluées en cherchant à maximiser ce déterminant. Pour un nombre de gaz n faible, par 15 exemple pour n 4, il est possible de recourir à une méthode exhaustive et tester toutes les longueurs d'onde de détection possibles. Pour un nombre de gaz n élevé, par exemple pour n 5, une méthode d'optimisation de type solveur d'équations peut être utilisée, correspondant par exemple à des algorithmes évolutionnaires tels que l'algorithme PSO, ou « Particle Swarm Optimization », décrit dans le document de J. Kennedy et al., « Particle 20 Swarm Optimization », Neural Networks, 1995. Proceedings., IEEE International Conference on, vol. 4, pp. 1942-1948, ou l'algorithme CMAES, ou « Covariance Matrix Adaptation Evolution Strategy », décrit dans le document de N. Hansen, « The CMA evolution strategy : A tutorial ». Eventuellement, une fois le déterminant de valeur maximale trouvé, il est 25 possible de calculer les valeurs exactes des erreurs sur les transmittances à partir de l'équation (11) afin de vérifier que les valeurs déterminées n'engendrent pas de problème de saturation. La méthode décrite ci-dessus donne de bons résultats pour un jeu de ?t, discrets, et est adaptée lorsque le capteur comporte des sources lumineuses laser qui sont 3042038 14 monochromatiques. Dans ce cas, le procédé permet d'identifier les longueurs d'onde ?t, optimales à émettre par les n sources lumineuses laser du capteur. Lorsqu'une source lumineuse polychromatique est utilisée avec des filtres de largeurs spectrales supérieures à celles des raies d'absorption des gaz à détecter, 5 au lieu d'utiliser directement les valeurs des coefficients d'absorption spectrale E(X) des gaz dans la matrice d'absorptivité E, les coefficients de la matrice d'absorptivité E sont choisis comme correspondant aux coefficients d'absorption spectrale moyens de chacun des n gaz pondérés par les réponses spectrales des n filtres spectraux. Ces coefficients sont calculés selon l'équation suivante : 10 IE./MA i (À) S (À) clÀ. ej (Ai) - f Ai(À.)S(À)dÀ. (16) avec 'Éi(À.,) le coefficient d'absorption spectrale moyen du gaz j par rapport à un filtre spectral dont la réponse spectrale type A, est positionnée sur la longueur d'onde )t,, et S(À) le spectre d'émission de la source lumineuse. Les indices i et j sont tels que 1 i n et 1 j n. Le terme EJ(X) correspond au coefficient d'absorption spectrale du 15 gaz j. Le calcul des coefficients 'Éi(À.,) est réalisé en intégrant la formule (16) ci-dessus sur tout le domaine où EJ(X) et Ai(Â) sont sensiblement non nuls, ce qui revient généralement à réaliser cette intégration sur le domaine correspondant à la largeur du filtre jusqu'à ce que Ai(Â) soit quasi nul. L'application de l'équation (16) ci-dessus correspond à l'application d'un 2 0 filtre d'un gabarit type donné (A), ce filtre étant répété pour différentes valeurs de ?t, pour simuler plusieurs positionnements de filtres. Dans un mode de réalisation particulier, lorsque les filtres à utiliser dans le capteur sont choisis parmi un ensemble donné de filtres dont les fonctions de transfert, ou réponses spectrales, sont définies (gabarit, largeur et longueur d'onde centrale définies), il est possible de calculer les valeurs Éii 25 correspondantes et de trouver le sous-ensemble de filtres qui permettent de maximiser la valeur du déterminant de la matrice d'absorptivité E. Les filtres les plus optimaux parmi ceux disponibles sont ainsi identifiés pour être utilisés dans le capteur. Lorsque les positions des filtres spectraux, c'est-à-dire les longueurs d'onde sur lesquelles les réponses spectrales des filtres sont positionnées, ne sont pas 3042038 15 définies et qu'elles peuvent être choisies, mais que les filtres spectraux ont un gabarit et une largeur définis (c'est-à-dire s'il est possible de positionner spectralement un filtre de réponse spectrale de forme A sur n'importe quelle longueur d'onde ) comme dans le cas de filtres interférentiels dont les paramètres de filtrage sont modifiés en modifiant les 5 paramètres des couches utilisées pour les réaliser), alors il est possible de construire la fonction Wi) en utilisant l'équation (16) précédente pour scanner toutes les valeurs possibles de positionnement ?t, dans la gamme d'intérêt. Dans le cas où, en plus de la position libre des filtres spectraux, les largeurs des réponses spectrales des filtres sont variables, l'opération décrite ci-dessus 10 peut être répétée pour chaque largeur possible de la réponse spectrale des filtres, et choisir les filtres optimums en faisant varier ces deux paramètres pour la création de la matrice E. Il est dans ce cas possible de réduire encore plus l'erreur en faisant correspondre les largeurs de filtres aux largeurs de raies d'absorption (par exemple en considérant la largeur à mi-hauteur des raies d'absorption). 15 Cependant, la possibilité d'avoir différentes largeurs de filtres va faire converger la recherche vers des filtres très étroits dont l'émissivité moyenne est plus importante. Dans le cas où la puissance spectrale délivrée par la source lumineuse est limitée (cas d'une source assimilable à un corps noir par exemple), cela signifie aussi réduire l'intensité spectrale du signal et augmenter l'erreur si le bruit est constant. 20 Lorsque les positions des filtres sont prédéfinies et ne sont pas modifiables, les coefficients de la matrice d'absorptivité E peuvent être calculés selon l'équation suivante : É.-"i (AI) = f si (À.) A i (À.) S (À.) dÀ (17) Ainsi, il est possible de facilement prendre en compte la quantité de flux 25 lumineux qui atteint le détecteur et qui est proportionnelle au terme f Ai(À.)S(À)dii., et trouver un optimum qui tienne compte de la quantité de flux lumineux qui atteint le détecteur, quelle que soit la largeur spectrale du filtre. Dans les cas décrits ci-dessus, l'utilisation des absorptivités moyennes pour former les coefficients de la matrice d'absorptivité E est préférable dans le cas de 3042038 16 faibles concentrations de gaz pour lesquelles le régime de l'exponentielle de la loi de BeerLambert est linéaire. En effet, la transmittance moyenne n'est pas la moyenne des transmittances, c'est-à-dire : (10-si(2)LCI, # 10-(si(2.))ALCi, (18) 5 De manière avantageuse, il est préférable de ne pas avoir un paramètre e(X).LC, qui concerne la transmittance du gaz seul, qui soit supérieur à environ 4 % pour garantir que la différence entre l'exponentielle et sa linéarisation reste inférieure à environ 5 %. Un exemple est décrit ci-dessous pour illustrer les méthodes 10 précédemment décrites. Cet exemple est mis en oeuvre pour déterminer les caractéristiques de trois filtres d'un capteur NDIR apte à discriminer trois gaz, ici trois hydrocarbures : le méthane CH4, l'éthane C2H6 et l'éthylène C2H4. Les spectres d'absorption, c'est-à-dire les coefficients d'absorption spectrale, de ces trois gaz sont représentés sur la figure 2. Cette 15 figure montre que ces spectres d'absorption se recouvrent de manière importante et que le choix des caractéristiques des filtres, notamment les longueurs d'onde de détection sur lesquelles les réponses spectrales des filtres sont destinées à être positionnées, n'est pas évident. Un filtre Lorentzien, dont le gabarit est représenté sur la figure 3, est par 20 exemple choisi dans cet exemple de mise en oeuvre. Ce gabarit, centré sur ?t, et de largeur r s'exprime selon l'équation : = 1 (19) rizÀ-À. D'autres types de filtres peuvent être utilisés, par exemple un filtre gaussien. 25 Les coefficients de la matrice d'absorptivité E sont calculés selon l'équation (16) ci-dessus. Les termes g'i(À.i) = f Ai(2.)sj(2i.)S(2)c12L de cette équation correspondent aux produits de convolution de la réponse du filtre avec le spectre d'absorption de chacun des trois gaz et sont représentés sur la figure 4. La référence 30 désigne le produit de convolution du spectre d'absorption du méthane avec le gabarit du 3042038 17 filtre représenté sur la figure 3, la référence 32 désigne le produit de convolution du spectre d'absorption de l'éthane avec le gabarit du filtre représenté sur la figure 3, et la référence 34 désigne le produit de convolution du spectre d'absorption de l'éthylène avec le gabarit du filtre représenté sur la figure 3. Une telle convolution a pour effet de lisser les spectres 5 d'absorption des gaz, ce qui permet de choisir une discrétisation spectrale plus grossière, et ainsi de réduire le temps de calcul nécessaire pour calculer le déterminant de valeur maximale de la matrice d'absorptivité E. Le déterminant de la matrice d'absorptivité E dont les coefficients correspondent aux absorptivités moyennes -Ei(J1) est calculé pour plusieurs positions 10 possibles des filtres. Les positions optimales des filtres, qui correspondent aux positions pour lesquelles la valeur de ce déterminant est maximisée, sont calculées et représentées sur la figure 5 dans laquelle ces positions sont superposées aux spectres d'absorption moyennés, et sur la figure 6 dans laquelle ces positions sont superposées aux spectres d'absorption initiaux des gaz. Ainsi, une première longueur d'onde de détection 36 sur 15 laquelle la réponse spectrale d'un premier des trois filtres destiné à détecter le méthane est calculée et est ici égale à environ 3,29 um, une deuxième longueur d'onde de détection 38 sur laquelle la réponse spectrale d'un deuxième des trois filtres destiné à détecter l'éthane est calculée et est ici égale à environ 3,35 um, et une troisième longueur d'onde de détection 40 sur laquelle la réponse spectrale d'un troisième des trois filtres destiné à 20 détecter l'éthylène est calculée et est ici égale à environ 3,18 um. Ces figures 5 et 6 illustrent bien que les positions optimales des filtres, c'est-à-dire les longueurs d'onde de détection sur lesquelles les réponses spectrales des filtres sont positionnées ne correspondent pas toujours aux valeurs maximales des spectres d'absorption des gaz à détecter. 25 Après le calcul de ces longueurs d'onde de détection, il est possible de calculer les erreurs sur les concentrations en utilisant l'équation (11) précédente à partir d'erreurs de transmission et de concentrations mesurées. Dans certains cas, il est impossible d'utiliser des filtres possédant une bonne réjection, c'est-à-dire coupant efficacement les longueurs d'onde hors de la bande 30 spectrale choisie. Dans le cas où il reste un fond d'absorption large bande, il est impossible 3042038 18 de créer une voie de référence qui soit insensible à la présence des gaz. Il est décrit ci-dessous le cas dans lequel le capteur ne comporte pas de voie de référence ou dans lequel l'une des voies de mesure est utilisée comme voie de référence. Dans une telle configuration, en considérant deux voies de mesure, la 5 transmittance est alors définie par l'équation suivante : V10/ ',./2.1 lo-Eic,irAc -Ei(4)Lc = 10_ T = (120)(4) =io REi (21)-ei (2,-))L Cl Ici, la deuxième voie est considérée comme la voie de référence, et £i(4) est l'absorptivité du premier gaz vue par le filtre associé à cette voie de référence. L'équation (20) est très proche de l'équation (2), la différence étant que le terme E1(À1) 10 est remplacé par le terme (E1(11) - £i(4)). Les étapes des procédés précédemment décrits peuvent donc s'appliquer en considérant, à la place des absorptivités de chacune des voies, les différences entre ces absorptivités et celle de la voie servant de référence. Lorsque n gaz sont considérés, n filtres de mesure et un filtre de référence indicé r sont utilisés. La matrice d'absorptivité E précédemment décrite dans l'équation (6) 15 est dans ce cas remplacée par la matrice suivante : 7(£1(L1) - El (Âr)) (E2(À1) - Ez(Âr)) ... (En(Âi) - En(Âr))\ (Ei (Â2 ) - El (Âr)) (E2(2L2) - Ez (Âr)) - (En (Â2 ) - En (Âr)) . . \ (Ei(Ân) - £1 (4)) (En (Â2) - Ez (Âr)) ... (En (Ân) - En (À,r))i Dans la matrice (21) ci-dessus, le vecteur ligne de l'absorptivité de la voie servant de référence est retranché à toutes les lignes de la matrice d'absorptivité initiale utilisée dans l'équation (6) précédente. Dans la méthode de calcul, cela revient à placer de 2 0 façon optimale n + 1 filtres de façon à maximiser le déterminant de la matrice (21). Il est à noter que les meilleurs résultats sont obtenus lorsqu'il est possible d'avoir une voie indépendante de tous les gaz, ce qui correspond à une voie de référence classique. Enfin, lorsque l'une des voies du capteur est utilisée comme voie de référence, le choix d'un filtre spécifique plutôt qu'un autre pour cette voie de référence n'a pas 25 d'influence sur le résultat obtenu. Dans les exemples précédemment décrits, le capteur utilisé correspond à un capteur comportant une cavité de longueur L connue, tel que le capteur NDIR 10. (20) (21) 3042038 19 Toutefois, ces procédés peuvent s'appliquer également à un imageur multi-spectral dans lequel un capteur est par exemple couplé à une roue à filtres tel que décrit dans le document « IR Gas Imaging in an Industrial Setting » de E. Naranjo et al., Proceedings of SPIE - The International Society for Optical Engineering (Impact Factor: 0.2). 04/2010, 5 vol. 7661, pp. 76610K-76610K-8, qui permet de sélectionner, via un choix de filtre, soit une voie active influencée par la présence du gaz, soit une voie de référence. La roue à filtres peut être remplacée par des absorbeurs sélectifs formés au niveau des pixels (résonateurs MIM, filtres pixelliques, etc.). Dans un tel système d'imagerie, le système optique est passif, et le 10 rayonnement assimilable à celui d'un corps noir situé en arrière-plan est utilisé comme source lumineuse. Le signal mesuré par un pixel filtré par une voie i est proportionnel à un flux reçu qui( est exprimé par l'équation suivante : (Pi = f [(PBB (TB)T (2-) + BB C) (1 - T(À))1ili(À)dÀ (22) avec ckBB la luminance du corps noir donnée par la loi de Planck, TB et Tc 15 respectivement les températures du fond et du gaz, et A,(X) la réponse spectrale du filtre i. En considérant la différence de flux entre deux pixels adjacents imageant deux zones de températures différentes TB' et TB2, un signal Ail), indépendant de la température du gaz est obtenu : J [oBB (712) (i) BB (711)1T (2-)Ai(1)CIÀ (23) 20 A partir du rapport entre le signal mesuré sur une voie active i et une voie de référence j, tel que T(À)Ai(À) = Ai(À) (le gaz n'absorbe pas dans la voie de référence), l'expression suivante proche de l'équation (2) est obtenue : f[OBB(TB2)-(PBB(TB1)1TMAiwcu. (24) ni] Cr rh (7, (1, (7, )1 n I J B2 J- Vvut/L Ai Les termes Â correspondent à des réponses spectrales moyennées sur À. 25 En calibrant le rapport À/24-, il est donc possible de remonter au produit longueur- concentration à l'aide de l'équation (3). REVENDICATIONS1. Procédé d'optimisation de valeurs de n longueurs d'onde de détection d'un capteur (10) de gaz optique apte à réaliser une détection de n gaz différents, comportant la mise en oeuvre des étapes suivantes : a) calcul d'une valeur d'un déterminant d'une matrice d'absorptivité E dont des coefficients sont représentatifs de l'absorptivité spectrale de chacun des n gaz aux n longueurs d'onde de détection, l'étape a) étant répétée plusieurs fois en modifiant à chaque fois au moins une des valeurs desdites n longueurs d'onde de détection et telle que les valeurs desdites n longueurs d'onde de détection soient comprises dans une gamme de valeurs pour laquelle l'absorptivité spectrale d'au moins un des n gaz est non nulle ; b) détermination des valeurs desdites n longueurs d'onde de détection pour lesquelles la valeur du déterminant de la matrice d'absorptivité E correspond à la valeur maximale parmi l'ensemble des valeurs calculées lors de l'étape a). 2. Procédé selon la revendication 1, dans lequel : - le capteur (10) de gaz optique comporte au moins n sources lumineuses laser, - les coefficients de la matrice d'absorptivité E correspondent, ou sont uniformément proportionnels, aux coefficients d'absorption spectrale de chacun des n gaz aux n longueurs d'onde de détection, et - les valeurs des n longueurs d'onde de détection déterminées à l'étape b) correspondent à des longueurs d'onde de faisceaux lumineux aptes à être émis par les n sources lumineuses laser. 3. Procédé selon la revendication 1, dans lequel : - le capteur (10) de gaz optique comporte au moins une source lumineuse polychromatique (12) et au moins n filtres spectraux (22) chacun apte à permettre une détection d'un des n gaz, 3042038 21 - les coefficients de la matrice d'absorptivité E correspondent, ou sont uniformément proportionnels, aux coefficients d'absorption spectrale moyens de chacun des n gaz, pondérés ou non par des réponses spectrales des n filtres spectraux (22), et - les valeurs des n longueurs d'onde de détection déterminées à l'étape 5 b) correspondent à des longueurs d'onde sur lesquelles les réponses spectrales des n filtres spectraux (22) sont positionnées. 4. Procédé selon la revendication 3, dans lequel les coefficients, notés É''./ (AI), de la matrice d'absorptivité E sont tels que : f si (À.) A i (À.) S (2)c1À 10 Ej f A i (À.) S (2)c1À ou tels que : Chi) = f (À.) AiMS (4c1À avec ?t, correspondant aux n longueurs d'onde de détection, i étant tel que 1 i n ; 15 EJ(X) correspondant aux coefficients d'absorption spectrale des n gaz, j étant tel que 1 j n ; A,(X) correspondant aux réponses spectrales des n filtres spectraux (22) ; S(À) correspondant au spectre d'émission de la source lumineuse (12). 20 5. Procédé selon la revendication 4, dans lequel : - les n filtres spectraux (22) sont choisis parmi un ensemble de m filtres spectraux, avec m > n, dont un gabarit, une valeur de longueur d'onde sur laquelle la réponse spectrale est positionnée et une largeur de la réponse spectrale de chacun des filtres spectraux dudit ensemble de filtres sont prédéterminés, et 25 - les valeurs des n longueurs d'onde de détection sont choisies parmi les valeurs des longueurs d'onde sur lesquelles les réponses spectrales des filtres spectraux dudit ensemble de filtres sont positionnées. 3042038 22 6. Procédé selon la revendication 3, dans lequel : - les n filtres spectraux (22) ont un même gabarit AM, et - les coefficients, notés (AI), de la matrice d'absorptivité E sont tels f A iM si MS (À) clÀ si (À ) = f A i(À.) S (À) clÀ ou tels que : -Si (À i) = f AiM si MS (À) clÀ avec ?t, correspondant aux n longueurs d'onde de détection, i étant tel que 1 i n ; 10 EJ(X) correspondant aux coefficients d'absorption spectrale des n gaz, j étant tel que 1 j n ; A,(X) correspondant aux réponses spectrales des n filtres spectraux (22) positionnées sur les longueurs d'onde A;; S(À) le spectre d'émission de la source lumineuse (12). 15 7. Procédé selon la revendication 6, dans lequel : - les largeurs des réponses spectrales des n filtres spectraux (22) sont ajusta bles, - l'itération de l'étape a) est répétée plusieurs fois en modifiant à 2 0 chaque fois la largeur des réponses spectrales des n filtres spectraux (22), et - l'étape b) détermine les valeurs des n longueurs d'onde de détection pour lesquelles la valeur du déterminant de la matrice d'absorptivité E correspond à la valeur maximale parmi l'ensemble des valeurs calculées lors des multiples itérations de l'étape a). 25 8. Procédé selon l'une des revendications précédentes, dans lequel : - le capteur (10) de gaz optique ne comporte pas de voie de référence et comporte n+1 voies de mesure ; que : 5 3042038 23 - les coefficients de la matrice d'absorptivité E correspondent, ou sont uniformément proportionnels, aux coefficients d'absorption spectrale, moyens ou non, de chacun des n gaz, pondérés ou non par des réponses spectrales des n filtres spectraux (22), auxquels sont soustrait les coefficients d'absorption spectrale d'une des n+1 voies de 5 mesure du capteur (10) de gaz optique servant de voie de référence. 9. Procédé selon l'une des revendications précédentes, dans lequel les étapes a) et b) sont mises en oeuvre via une méthode d'optimisation de type solveur d'équations par algorithme évolutionnaire.