Brevet FR2537722A1 — publié le 19840615

La présente invention concerne de manière générale la détection de gaz explosifs dans un milieu gazeux tel que l'air. De manière connue en soi, la détection d'un gaz explosif dans l'air est assurée par un appareil, usuelle- ment appelé explosimètre à filament catalytique, dans le- quel un filament, généralement de platine, est chauffé par effet Joule, c'est-à-dire par passage d'un courant électrique Le gaz explosif contenu dans l'air s'oxyde par catalyse au contact du filament, d'ou un échauffement supplé- mentaire de ce dernier L'élévation de température qui en résulte provoque une élévation de la résistance du fila- ment, dont la mesure permet d'avoir accès à la concentra- tion dudit gaz explosif dans l'air En pratique, on s'in- téresse au taux d'explosivité du gaz, c'est-à-dire au rapport de sa concentration par sa limite inférieure d'explo- sivité (LIE), c'est-à-dire la teneur en ce gaz au-dessus de laquelle il y a risque d'explosion; c'est pourquoi le résultat de la détection est usuellement présenté sous forme d'un pourcentage de la limite LIE précitée. A titre d'exemple, la Demanderesse a proposé dans son brevet FR-1 444 771 déposé le 4 mai 1965, un appareil de prélèvement et de dosage de gaz dans l'air qui est por- tatif Le dosage de gaz y est effectué de manière classique par mesure de la tension aux bornes de l'une des diagonales d'un pont résistif constitué par un filament détecteur et un filament compensateur montés en parallèle avec deux résistances dont l'une est avantageusement réglable, ce pont étant alimenté électriquement selon son autre dia- gonale. Dans un cadre plus général, la Demanderesse a également proposé, dans son brevet FR-1 577 448, déposé le 22 mai 1968, un procédé et un appareil de mesure d'une grandeur caractéristique d'un milieu gazeux, selon lesquels la valeur de la grandeur considérée est estimée à partir de la mesure d'une des données d'alimentation d'un filament détecteur lorsque la résistance de ce dernier est maintenue égale à celle d'un filament compensateur. A l'usage, de tels explosimètres se révèlent donner de très bons résultats pour le gaz explosif qui a servi à leur étalonnage, tandis que pour d'autres gaz, les taux de concentration estimés s'écartent très sensiblement des taux réels. Ces écarts résultent en partie des différences existant entre les gaz concernant notamment leur chaleur d'oxydation, leur température de combustion, leur conducti- bilité thermique, et surtout leur coefficient de diffusion dans l'air (l'hydrogène par exemple diffuse quatre fois plus vite que le méthane) Ces différences influent sur le renouvellement de gaz explosif auprès du filament, la qualité de la combustion (incomplète ou précoce selon la température du filament), le nombre de calories ainsi dispo- nibles, et leur évacuation dans le milieu gazeux environ- nant Il est à noter que ces écarts correspondent souvent à des sous-estimations des valeurs réelles, ce qui est défa- vorable du point de vue sécurité. L'invention vise à réduire ces écarts et a pour objet un procédé perfectionné de détermination du taux d'explosivité d'un milieu gazeux, ainsi qu'un dispositif pour sa mise en oeuvre, qui permettent d'échapper, au moins partiellement, à l'influence néfaste des sources d'écart précitées. L'invention propose à cet effet un procédé de détermination du taux d'Explosivité d'un milieu gazeux contenant au moins un gaz explosif, selon lequel on chauffe dans ce milieu un filament détecteur par passage de courant électrique, caractérisé en ce que: on allume le filament par connexion à une source d'alimentation; on régule l'alimentation du filament détecteur de manière à en maintenir la résistance constante; on prend en compte, en plusieurs instants successifs comptés à partir de l'allumage du filament, la valeur U(t) d'une grandeur U significative de l'alimentation dudit filament; on extrapole à partir desdites-valeurs U(t) une courbe exponentielle d'équation U(t) Uo exp (-k t) dont on identifie les paramètres U O et k; et on introduit les paramètres U O et k dans une fonction de traitement F adaptée à donner indépendamment. de la nature du gaz explosif présent dans ce milieu gazeux, à partir de ces paramètres, une estimation du taux dlexplosivité initial dudit milieu gazeux. Ce procédé de l'invention est fondé sur la décou- verte par la Demanderesse qu'il est possible, de manière surprenante, d'établir une fonction empirique simple,, permettant de relier les deux paramètres précités, U O coefficient préexponentiel, et k exposant, au taux d'ex- plosivité x O d'un mélange gazeux, qui soit sensiblement indépendante de la nature du ou des gaz en présence La précision sur le résultat dépend bien sûr de la précision avec laquelle a été établie ladite fonction empirique-; il faut cependant garder à l'esprit qu'une trop grande précision serait illusoire puisque la limite inférieure d'explosivité est une notion imprécise en elle-même, sa mesure dépendant notamment de la méthode et du matériel expérimentaux utilisés, et de paramètres tels que la tem- pérature. Dans le cas de mélanges de gaz explosifs, la mise en oeuvre de l'invention nécessite le choix préalable d'une hypothèse concernant la relation entre le taux d'explosivité du mélange, et ceux des gaz explosifs pris individuellement En première approximation, on admet généralement, conformément à la loi de Le Chatelier, qu'il suffit de faire la somme de ces derniers D'autres hypo- thèses peuvent néanmoins être proposées le cas échéant, sans sortir du cadre de l'invention. Selon l'invention, grâce à la mise en oeuvre d'une fonction F -udicieusezent choisie, le taux d'explosivité d'un mélange gazeux contenant un ou plusieurs gaz explosifs est estimé avec une précision sensiblement indépendante de la nature du ou des gaz explosifs présents. Il en résulte pour le procédé selon l'invention, et donc pour les dispositifs assurant sa mise en oeuvre, un caractère universel puisqu'ils sont adaptés, sans modification, à estimer le taux d'explosivité d'une grande variété de milieux gazeux. Il est à noter en outre que selon l'invention, on prend en compte l'écoulement du temps, et que l'on a accès au taux d'explosivité initial x O du milieu gazeux, c'est-à-dire avant toute perturbation liée aux opérations de mesure. Dans une forme préférée de l'invention, on met en oeuvre, en outre, un filament compensateur identique au filament détecteur, alimenté dans des conditions ana- logues, dans un milieu gazeux analogue à celui du fila- ment détecteur mais ne contenant aucun gaz explosif, et que la grandeur significative U de l'alimentation du fila- ment détecteur est la différence entre les tensions d'ali- mentatiôn des filaments détecteur et compensateur On élimine ainsi les fluctuations de l'environnement pouvant affecter les opérations de mesure et leur résultat. L'invention préconise tout particulièrement d'uti- liser comme fonction de traitement F le rapport, à un coefficient multiplicatif près, du paramètre préexponen- tiel U O de la fonction exponentielle extrapolée par l'expo- sant k de cette dernière Une telle fonction, extrêmement simple, ne nécessite la mise en oeuvre cue de moyens de calcul rudimentaires D'autres fonctions simples sont proposées par l'invention: elles diffèrent de la précédente par le remplacement du dénominateur par une fonction polynome de l'exposant du premier degré, du second voire du troisième degré en fonction de la précision recherchée, du nombre et de la nature des gaz que l'appareil peut avoir à doser. Par souci de simplification, l'invention préconise en outre de déterminer les paramètres de la loi exponen- tielle extrapolée en utilisant une expression approchée de cette dernière, c'est-à-dire son développement limité du premier ordre. Compte tenu de la disparité de comportement existant entre les gaz explosifs, notamment l'hydrogène vis-à-vis des autres hydrocarbures, l'invention préconise en outre d'introduire les paramètres U et k dans une fonction choisie en fonction de la position de l'exposant k vis-à-vis d'une valeur critique k Il se trouve en effet que la Demanderesse a découvert, dans le cas de l'hydro- gène, qui se comporte différemment des autres hydrocar- bures, une différence concomitante entre les exposants exponentiels auxquels ce gaz conduit, et ceux des autres hydrocarbures qui se révèlent rester à l'intérieur d'une plage étroite. Pour la mise en oeuvre du procédé précité, l'in- vention préconise un dispositif de détermination du taux d'explosivité d'un milieu gazeux, caractérisé en ce qu'il comporte. un filament détecteur disposé dans une cellule de mesure enfermant le milieux gazeux; un circuit régulateur comportant une source d'alimentation, adapté à alimenter électriquement ledit filament détecteur de manière à en maintenir constante la résistance; un convertisseur analogique-numérique commandé par un compteur-temporisateur; et une unité de traitement numérique adaptée à déterminer une estimation du taux d'explosivité du milieu gazeux. Compte tenu du caractère universel du procédé dont un tel dispositif assure la mise en oeuvre, et de la simplicité des éléments qui le composent, il résulte pour ces dispositifs une standardisation des séries de fabri- cation avec des coûts unitaires modérés conduisant à une large diffusion de ces dispositifs L'invention assure de ce fait un sensible accroissement de la sécurité dans tous les lieux o existent, ou risquent d'apparaître des gaz explosifs. D'autres objets, caractéristiques et avantages de l'invention ressortent de la description qui suit, donnée à titre d'exemple, en regard des dessins annexés sur lesquels: la figure 1 est un graphique indiquantpour différents hydrocarbures, les taux d'explosivité estimés par un explosimètre connu en fonction des taux réels; la figure 2 est un synoptique d'un dispositif selon l'invention; la figure 3 est une vue détaillée de la partié analogique de ce dispositif, la figure 4 est une vue schématique de la partie numérique de ce dispositif; la figure 5 est un graphique indiquant l'évolu- tion au cours du temps de la grandeur U; la figure 6 est un graphique corrélant, pour différents hydrocarbures et l'hydrogène, le paramètre UO et le taux d'explosivité réel dans l'air; la figure 7 est un graphique corrélant, pour divers hydrocarbures et l'hydrogène, le rapport U /k avec le taux d'explosivité réel dans l'air; la figure 8 est un diagramme indiquant le taux d'explosivité d'un mélange de méthane et de décane, estimés suivant un procédé selon l'invention, en fonction de la teneur en méthane; la figure 9 est un diagramme analogue au précé- dent, correspondant à un mélange d'hydrogène et de méthane; la figure 10 est un graphique corrélant, pour divers gaz explosifs, une grandeur caractéristique K 2 avec le coefficient de diffusion dans l'air D y la figure 11 est un graphique corrélant pour divers hydrocarbures, le rapport U/XO avec le paramètre k; et la figure 12 est un diagramme analogue à celui de la figure 9, estimé suivant un procédé amélioré selon l'invention. La figure 1 illustre les écarts entre les taux d'explosivité réels et mesurésî pour l'hydrogène et divers hydrocarbures, associés à un explosimètre de type connu, dont la référence est MSA DGE 2000 On retrouve que ces écarts correspondent presque toujours à une sous- estimation des taux réels, ce qui est préjudiciable à la sécurité. Sur ce graphique, ainsi que sur les suivants, les gaz explosifs sont désignés par des symboles dont la signification est précisée dans la table ci-dessous H Hydrogène C 1 Méthane C 2 Ethylène-Ethane C 3 Propane C 4 Butane C 5 Pentane C 7 Heptane-Acetone C 8 Octane C 1 Décane L'examen du graphique de la figure 1 révèle que les taux d'explosivité mesurés pour un hydrocarbure donné sont d'autant plus sous-estimés que cet hydrocarbure admet des molécules encombrantes, c'est-à-dire (en première approximation) qui diffusent lentement dans l'air vers le filament (et inversement) Ces écarts semblent donc pouvoir être attribués pour une grande part aux différences de coefficient de diffusion dans l'air existant entre ces gaz. L'hydrogène, qui n'est pas un hydrocarbure, se place, malgré sa faible masse moléculaire, entre le méthane et l'éthylène. On observe que, pour un taux d'explosivité x O réel donné, il y a un écart de près de 70 % entre les résul- tats obtenus pour les gaz extrêmes, le méthane et l'heptane. Cela explique pourquoi, à ce jour, on utilise autant d'ex- plosimètres qu'il y a de gaz explosifs à doser. Selon l'invention, on pallie cet inconvénient, au moins partiellement, à l'aide d'un dispositif dont la structure générale est présentée à la figure 2 Ce dernier se décompose en une partie analogique de mesure 1 suivie d'un convertisseur analogique-numérique CAN puis d'une unité de traitement numérique 2. L'unité de mesure 1 comporte un filament détecteur D disposé dans une cellule de mesure 10 enfermant une partie du milieu gazeux à doser Entre ses bornes, dont l'une est à la masse, est disposé un dispositif régulateur 11 adapté à réguler l'alimentation du filament D de manière à en maintenir constante la résistance malgré les réactions d'oxydation qu'il provoque, le cas échéant, dans le milieu gazeux enfermé dans la cellule 10. De manière avantageuse, l'unité de mesure 1 comporte en outre un filament compensateur C, identique au filament détecteurqui est placé dans un environnement gazeux, ana- logue au milieu gazeux considéré mais exempt de tout gaz explosif, enfermé dans une seconde cellule de mesure 12. Le filament compensateur participe à un montage identique à celui du filament détecteur, un régulateur 13, analogue au régulateur 11, en modulant l'alimentation de manière à en maintenir constante la résistance. Les bornes des filaments qui ne-sont pas à la masse sont reliées aux entrées d'un comparateur 14, adapté par exemple à mesurer la différence de tension existant entre ces bornes Cette comparaison des alimentations des deux filaments permet de supprimer l'influence de para- mètres expérimentaux susceptibles de fluctuer, notamment les variations de température extérieureàet l'usure du filament détecteur Sur ce dernier point, il est précisé que les filaments détecteur et compensateur sont avanta- geusement remplacés simultanément. La sortie du comparateur 14 est connectée à l'entrée du convertisseur analogique-numérique CAN suivi d'un élément de calcul 3 complété par un organe d'affi- chage 4. La structure de la partie analogique de mesure 1 est détaillée à la figure 3 Les régulateurs 11 et 13 ont la même structure La description qui suit est limitée au régulateur 11 Les éléments du régulateur 13 portent des références obtenues par addition de 10 à celles des éléments correspondants du régulateur 11. Le régulateur 11 du filament détecteur D comporte une barre d'alimentation 20 à potentiel constant ( 5 V dans l'exemple considéré) Une chaîne résistive formée de résis- tances R 1, R 2, R 3 montées en série est montée entre cette barre 20 et la masse Les résistances R 2 et R 3 sont montées 2537 t 22 entre le collecteur et l'émetteur d'un transistor T 1, formant élément d'allumage, dont la base est reliée au tra- vers d'une résistance R à une borne AF d'allumage de filament La borne commune des résistances R 1 et R 2 est à la base d'un transistor T 2 dont l'émetteur à son tour est à la base d'un second transistor T 3 Les collecteurs des transistors T 2 et T sont à la barre 20 tandis que 2 3 l'émetteur de T 3 est relié à la masse au travers d'un pont résistif auquel participe le filament D Ce pont comporte d'une part une résistance R 7 et D, et d'autre part une résistance R 6, en série avec une résistance variable Pl, et une résistance R 7 Un amplificateur opérationnel A 01 dont les entrées sont montées selon la diagonale trans- versale du pont résistif a sa sortie reliée à la jonction des résistances R 2 et R 3 La jonction R 6 R 7 est avanta- geusement reliée à la barre d'alimentation 20 au travers d'un circuit (C 1 + R 8) d'amortissement des oscillations dues à l'inertie thermique du filament. La borne AF à laquelle sont reliées les résistances analogues R 4 et R 14 est reliée à la barre 20 au travers d'une résistance R 9. Les bornes des filaments D et C reliées aux ampli- ficateurs opérationnels associés sont connectées au travers de résistances R 20 et R 21 aux entrées d'un amplificateur opérationnel de comparaison AO 3 dont la sortie est à l'entrée analogique du convertisseur analogique-numérique CAN Cette sortie de l'élément de comparaison A 03 est couplée, par C 3 et R 22 en parallèle, à son entrée venant du filament détecteur, tandis que l'entrée venant du filament compensateur est à la masse au travers d'une résistance R 23. Le convertisseur CAN est classiquement à sortie parallèle 12 bits Deux potentiomètres P 3 et P 4 en per- mettent le réglage en gain et le décalage en zéro-respecti- vement. L'unité de traitement numérique 2 est schématisée à la figure 4 Elle comporte essentiellement un micro- processeur MP entouré de ses périphériques, suivi de l'élément d'affichage 4. Le microprocesseur MP est relié au convertisseur CAN à travers un port d'entrée PE qui comporte en outre des bits pour commande par un interrupteur auxiliaire IA adapté à sélectionner, le cas échéant, une option de calcul numérique et pour un bouton-poussoir PA adapté à envoyer un signal d'interruption au microprocesseur, en plus de ceux associés à chaque conversion. Les périphériques du microprocesseur MP sont essentiellement une mémoire programme EPROM, une mémoire de données RAM, des décodeurs pour des sélections de cir- cuits et un compteur-temporisateur CT qui, entre autres fonctions, commande par exemple toutes les secondes, une conversion analogique-numérique (signal CC). Les résultats apparaissent sur les bits de sortie du port de sortie PS Ils sont pris en compte par des décodeurs(trois dans l'exemple considéré) qui commandent un nombre égal d'afficheurs LED 7 segments Le port PS comporte également des bits pour la commande des points décimaux, ainsi que la commande d'allumage des filaments (borne AF). Un poussoir de mesure PM commande la remise à zéro de tous les circuits et le démarrage du microprocesseur en début de programme. Les interconnexions entre les divers éléments sont faites de façon standard, à l'exception des bornes M 1 des ports d'entrée et de sortie qui sont reliées au poussoir PM. Au repos, le signal appliqué en AF rend les transistors T 1 et Til passants, court-circuitant les éléments d'alimentation T 2 3, T 12 et T 13: aucun courant ne traverse les filaments. Une pression sur le poussoir PM provoque l'appari- tion d'un signal en AF Les bases des transistors T 1 et T 9 i sont à un niveau bas de sorte que ces derniers sont bloqués, tandis que les transistors T 2 T 3 " T 12 et T 13 deviennent passants Dans le cas par exemple du filament détecteur, covamndés par le comparateur AO 1 en fonction de la différence de tension entre les bornes d'entrée de ces derniers, les transistors T 1 et T 2 régulent le courant dans le pont résistif de D de manière que la résistance du filament détecteur soit égale à (R 5 * R 7)/(P 1 + R 6), cette valeur étant ajustable par action sur le potentiomètre P 1 Le même raisonnement s'applique au filament compensateur On notera que ce système accélère le chauffage des filaments puisqu'il les maintient en surtension tant que la tempéra- ture voulue n'est pas atteinte. Les tensions existant aux bornes des filaments sont soustraites l'une de l'autre par le comparateur A 03, et leur différence est convertie, à instants répétés, par le convertisseur CAN qui assure ainsi, en des instants successifs, la prise en compte d'une grandeur significa- tive U de l'alimentation des filaments. L'unité de traitement numérique 3 prend en compte plusieurs valeurs successives U(t) de cette grandeur U, puis extrapole à partir de ces valeurs une fonction expo- nentielle d'équation U(t) U= exp (-k t), dont on identifie les paramètres U et k Ceux-ci sont alors introduits dans o une fonction de traitement définie empiriquement, qui donne à partir de ces paramètres une estimation du taux d Dexplo- sivité initial x O du mélange gazeux considéré. Un dispositif du genre précité, et le procédé qu'il met en oeuvre permettent d'échapper à la nécessité d'identifier le ou les gaz explosifs présents dans le milieu considéré dans la mesure ol la Demanderesse a découvert qu'il était possible de déterminer empiriquement une telle fonction de traitement xo = F (U 0, k) qui soit indépendante de la nature du ou des caz présents. La forme d'une telle fonction de traitement s'établit notamment sur la base des considérations théori- ques suivantes. Soit /G 7 la concentration d'un gaz G dilué dans un milieu tel que l'air, AH son enthalpie d'oxydation, et D son coefficient de diffusion dans le milieu gazeux considéré Il ressort d'un article de FIRTH, JONES et JONES (Combustion and Flame, 21 ( 1973), 303-311) qu'au contact d'un filament chauffé, le gaz G fournit à ce dernier une énergie Q qui, en régime permanent et à température cons- tante, en négligeant les phénomènes catharométriques, s'écrit: Q = K AH D / G 7 ( 1) o K est une constante caractéristique du filament. En passant de l'énergie Q à la puissance calori- fique P dégagée par la combustion du gaz G, et en intro- duisant la notion de taux d'explosivité x défini par: x = /G 7 / G 7 LIE ( 2) on obtient P = K 1 AH G 7 LIE D x ( 3) qui devient: P = K 1 K 2 Do x( 4) en posant K 2 = AH /-G/LIE ( 5) FIRTH, JONES et JONES ont en outre observé que le coeffi- cient K 2, qui est une grandeur caractéristique pour chaque gaz, varie peu d'un hydrocarbure à'un autre. En pratique, dans une cellule de mesure, le gaz s'épuise au cours du temps de sorte que le régime n'est pas permanent Compte tenu des faibles concentrations sous lesquelles sont présents les gaz à mesurer, on approxime la consommation de gaz à l'aide d'une cinétique du premier ordre On en déduit: x = x O exp (-kt) ( 6) x étant le taux d'explosivité en l'instant t = 0. En comparant l'expression de l'énergie calorifi- que fournie par le gaz pendant sa combustion obtenue par intégration de l'équation ( 4) entre O et l'infini, et son expression thermochimique, il découle que k est propor- tionnel au coefficient de diffusion Do. On vérifie que U est une grandeur expérimentale représentative de la puissance P qui admet pour expression approchée: U = U O exp (-k t) ( 6) avec Uo = 33 2 ' D xo ( 7) On en déduit que le rapport U /k ne dépend pas explicitement de D, d'o l'expression de x O indépendante du coefficient D U xo = F (UO, k) = K 4 K 2 k ( 8) K 4 étant une constante. Dans la mesure o l'on admet, en s'inspirant des observations précitées de FIRTH et al, que K 2 ne dépend pas de la nature des gaz, il s'avère que la fonction F, définie par l'équation ( 8), constitue une fonction de traitement au sens précisé plus haut. 1 Compte tenu de ce que l'équation ( 1) n'est valable qu'à température constante, l'équation de traitement déduite de l'équation ( 8) nécessite le maintien de la résistance des filaments à une valeur constante, d'o la présence des régulateurs 11 et 13 dans un dispositif selon l'invention. Par mesure de simplification, l'invention préco- nise de déterminer les paramètres préexponentiel U O et exposant k à partir du terme du premier ordre du développe- ment limité de l'équation exponentielle: U(t) = U O ( 1 k t) ( 9) La détermination des paramètres U O et k ne nécessite que la prise en compte par l'unité de traitement 3 de données correspondant à deux instants Ceuxci correspondent avantageusement à des durées de mesures qui correspondent à des rapports simples, par exemple du simple au double, ce qui simplifie considérablement les calculs. L'utilisation de données correspondant à plus de deux instants peut être préférée pour améliorer la précision sur U O et k à condition que les temps correspondants soient compatibles avec les hypothèses simplificatrices utilisées. L'invention préconise en fait la prise en compte des valeurs de la grandeur significative d'alimentation choisie en trois instants correspondant à des intervalles de temps satisfaisant à des rapports simples (par exemple ls, 2 S et 4 s); on traite successivement au moins deux des couples de valeurs ainsi obtenus (par exemple ls et 2 s 2 S et 4 s) et l'on compare les valeurs UO et k qui s'en déduisent La comparaison de ces dernières donne une esti- mation de la confiance à accorder à l'estimation finale du taux d'explosivité x 0, voire sert de base à l'obtention de valeurs corrigées. La figure 5 présente l'évolution du signal U en fonction du temps écoulé depuis l'allumage des filaments, pour un hydrocarbure donné, pour divers taux d'explosivité. On observe que les courbes, au-delà de ls environ, peuvent se décrire par une loi du type exponentiel ce qui justifie la relation ( 6) Un régime transitoire de chauffage de fila- ment et de stabilisation (entre 0,5 et ls environ), noté T, apparaît au début de ces courbes, rendant impossible toute mesure. Les figures 6 et 7 corrèlent les taux d'explosi- vité estimés par un dispositif selon l'invention pour deux fonctions de traitement simples (après prise en compte de U après 2 S et 4 s), avec les taux réels. La figure 6 correspond à une fonction de traite- mentsiitplifiée à l'extrême puisqu'elle conduit à estimer le taux d'explosivité à partir du seul paramètre U O; cette fonction présente cependant, par rapport à l'art antérieur, l'avantage de tenir compte de l'évolution du temps au cours des opérations de mesure Les dispersions entre courbes ex- trêmes restent cependant du même ordre de grandeur que celles de la figure 1. La mise en oeuvre comme fonction de traitement de la relation ( 8) conduit par contre à des résultats nettement plus sensibles puisque la figure 7 ne fait plus apparaître qu'un écart maximal de 30 % environ, soit pratiquement une réduction de moitié par rapport à la figure 6. Il est à noter que sur la figure 7, les courbes ne peuvent plus se classer comme sur la figure 1, en fonction des tailles moléculaires. La figure 7 comporte deux courbes associées à l'hydrogène Celle correspondant à H* résulte d'une cor- rection supplémentaire nécessitée par les caractéristiques particulières de l'hydrogène par rapport aux hydrocarbures. L'hydrogène présente en effet un coefficient K 2 cinq fois plus faible que pour les hydrocarbures mais un coefficient D plus élevé au traitement numérique des données une opération de compa- raison de l'exposant k de l'exponentielle à une valeur k o choisie entre la gamme de valeurs associée à l'hydrogène et celle associée aux hydrocarbures Si k est inférieur à k 0, on effectue les opérations associées à la relation ( 8), sinon on fait intervenir un coefficient multiplicatif spécifique de l'hydrogène. Les figures 8 et 9 correspondent au dosage, non plus d'un gaz seul, mais d'un mélange de deux gaz Les courbes qui y sont représentées résultent de l'application du traite- ment numérique qui vient d'être défini aux valeurs expéri- mentales indiquées dans l'article précité de FIRTH et al. Dans le cas des deux hydrocarbures les plus différents du point de vue diffusion (le méthane et le décane), on constate que la fonction de traitement de l'équation ( 8) conduit à une sous-estimation du taux d'à peine 20 % au maximm, ce qui est pleinement satisfaisant. La relation ( 8) s'avère ainsi particulièrement bien adaptée au dosage des hydrocarbures, qu'ils soient seuls ou en mélange. Lorsque l'hydrogène participe au mélange gazeux, il n'en est pas de même, ainsi qu'il apparaît à la figure 9. La mise en oeuvre de la correction précitée, spécifique à l'hydrogène, conduit à un segment vertical conduisant d'une sous-estimation du taux réel de 40 % environ à une sur- estimation de 60 %. La position horizontale de ce segment dépend de la valeur critique ko choisie La courbe du bas correspond à l'application de la relation ( 8) qui conduit pour l'hydro- gène à une nette sous-estimation de son taux d'explosivité (voir figure 7, courbe EH) La courbe haute correspond à une multiplication des valeurs de la courbe basse par le coefficient multiplicatif spécifique de l'hydrogène. Ces disparités dans les résultats proviennent en fait de ce que, ainsi que cela apparaît à l'examen de la figure 10, le coefficient 1 à l'autre, d'un gaz explosif à l'autre (H, Co). De manière à tenir compte de ces variations, l'invention préconise, lorsque cela apparaît utile, de modifier en conséquence la forme de la fonction de traitement mise en oeuvre pour le traitement de U 0 et k. L'invention propose ainsi une première forme amendée de la fonction de traitement U x = F' (U, k) = B ( 10) 0 O k+b dont la justification apparaît à la figureil,sur laquelle est reporté, en fonction de la valeur mesurée du paramètre k, le rapport du second paramètre U O par le taux réel d'explosivité x O pour plusieurs hydrocarbures -On peut tracer une droite,)passant sensiblement par ces points, mais celle-ci ne passe pas par l'origine comme le sous- entendait la retation ( 8), d'o l'intervention d'un coeffi- cient b dans la relation ( 10). Une autre fonction de traitement peut se déduire de la figure 10 précitée à partir de la droite qui y est * tracée, passant au plus près des points qui y sont reportés. En reportant l'expression de K 2 en fonction de D qui peut s'en déduire dans l'équation ( 7) et en éliminant D entre cette équation et l'expression de k, on définit une nouvelle fonction de traitement F" analogue à F', mais dont le dénominateur est un polynôme en k du second degré. Si l'on reporte sur la figure 11 d'autres points tels que l'hydrogène ou d'autres gaz explosifs différents des hydrocarbures, l'ensemble des points ainsi tracés se révèle pouvoir être décrit de façon convenable par une courbe-polynôme du troisième degré En reportant ce poly- nôme au dénominateur de la fonction de traitement, la Demanderesse a montré que cette nouvelle fonction de trai- tement conduit à une bonne approximation du taux d'explosivité d'un mélange contenant de l'hydrogène ainsi qu'en témoigne la figure 12 qui correspond au mélange méthane-hydrogène de la figure 9. Il s'avère ainsi que de nombreuses fonctions de traitement peuvent être Proposées dans le cadre de l'invention Ces fonctions, définies empiriquement ont bien sûr une forme d'autant plus complexe que l'on recherche une grande précision, qu'il y a de nombreux gaz en présence, et que le nombre de gaz susceptibles d'être dosés, seuls ou en mélange avec d'autres, est élevé. Dans le cas de la fonction de traitement définie par l'équation ( 8), l'invention préconise en outre d'éta- lonner le dispositif de dosage avec du méthane: ainsi que cela ressort de la figure 7, toutes les mesures seront ainsi sur-estimées, du moins dans le cas de gaz pratiquement purs ou presque, ce qui est favorable à la sécurité. Il va de soi que de nombreuses variantes de réalisation peuvent être proposées par l'homme de l'art sans sortir du cadre de l'invention, tant en ce qui concerne le filament (nature), la cellule de mesure (forme, volume), le dispositif régulateur, le comparateur, le convertisseur, que l'unité de traitement numérique et ses éléments d'affichage De même, les paramètres expéri- mentaux peuvent être choisis dans une large gamme de valeurs (température ou résistance de filament). L'invention a été décrite dans le cas d'un milieu gazeux à base d'air; elle s'applique bien str à tout autre environnement gazeux dont le taux d'explosivité doit rester limité. L'invention trouve ainsi un champ d'application très vaste, incluant les mines et les industries chixniquesp 1 'on veut se pré- server de toute fuite de gaz explosif Un dispositif selon l'invention peut alors être fixe, assurant automatiquement, en cas de dépassement d'un taux d'explosivité critique, le déclenchement de signaux d'alarme et de mesures d'urgence. REVENDICATIONS 1 Procédé de détermination du taux d'explosivité d'un milieu gazeux contenant au moins un gaz explosif, selon lequel on chauffe dans ce milieu un filament détecteur par passage de courant électrique, caractérisé en ce que on allume le filament par connexion à une source d'alimentation ( 20); on régule l'alimentation du filament détecteur (D) de manière à en maintenir la résistance constante; on prend en compte, en plusieurs instants successifs comptés à partir de l'allumage du filament, la valeur U(t) d'une grandeur U significative de l'alimenta- tion dudit filament; on extrapole à partir desdites valeurs U(t) une courbe exponentielle d'équation U(t) U O exp (-k t) dont on identifie les paramètres U O et k; et on introduit les paramètres U O et k dans une fonction de traitement F adaptée à donner, indépendamment de la nature du gaz explosif présent dans ce milieu gazeux, à partir de ces paramètres, une estimation du taux d'explosivité initial (x) dudit milieu gazeux. 2 Procédé selon la revendication 1, caractérisé en ce que, en outre: on allume, en même temps que le filament détecteur un filament compensateur (C), identique à ce dernierplacé dans un-environnement analogue à celui de ce dernier mais exempt de gaz explosifs, par connexion à une source d'alimentation ( 30); on régule l'alimentation de ce filament compen- sateur de manière à en maintenir constante la résistance; et on obtient la grandeur U par soustraction de la tension d'alimentation du filament compensateur à celle du filament détecteur. 3 Procédé selon la revendication 1 ou la revendication 2, caractérisé en ce que l'estimation du taux d'explosivité initial est obtenue par division de U par un polynôme en k. 4 Procédé selon la revendication 3, caractérisé en ce que l'estimation du taux d'explosivité initial est obtenue, à une constante multiplicatrice près, par division de U O par k. 5 Procédé selon l'une quelconque des revendica- tions 1 à 4, caractérisé en ce qu'on détermine les para- mètres UO et k à partir des valeurs U(t) en remplaçant l'équation de la courbe exponentielle par le terme du premier ordre de son développement limité. 6 Procédé selon l'une quelconque des-revendica- tions 1 à 5, caractérisé en ce que les temps associés à la prise en compte de deux valeurs successives U(t), comptés à partir de l'allumage du filament détecteur, sont dans un rapport 1 à 2. 7 Procédé selon la revendication 6, caractérisé en ce que des valeurs Uit) sont prises en compte en trois instants successifs. 8 Procédé selon l'une quelconque des revendica- tions 1 à 7, caractérisé en ce que l'on introduit les para- mètres UO et k dans l'une de plusieurs fonctions de traite- ment choisie en fonction de la position du paramètre k par rapport à une valeur critique ko. 9 Dispositif de détermination du taux d'explo- sivité d'un milieu gazeux, pour mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce qu'il comporte: un filament détecteur (D) disposé dans une cellule de mesure ( 10) enfermant le milieu gazeux; un circuit régulateur ( 11) comportant une source d'alimentation ( 20), adapté à alimenter électrique- ment ledit filament détecteur de manière à en maintenir constante la résistance; un couvertisseur analogique-numérique (CAN) commandé par un compteur-temporisateur (CT); et une unité de traitement numérique ( 12) adaptée à déterminer une estimation du taux d'explosivité du milieu gazeux. Dispositif selon la revendication 9, caracté- risé en ce qu'il comporte en outre: un filament compensateur (C) identique au filament détecteur, disposé dans une cellule de mesure ( 12) voisine de celle du filament détecteur mais ne contenant aucun gaz explosif; un circuit régulateur ( 13) analogue à celui associé au filament détecteur; et un élément de comparaison ( 14) adapté à fournir la différence entre les tensions d'alimentation des fila- ments détecteur et compensateur à ln trée du convertisseur analogique-numérique (CAN). 11 Dispositif selon la revendication 9 ou la revendication 10, caractérisé en ce que le régulateur ( 11, 13) définit avec le filament (D, C) un pont résiÈtif (P 1, R 5 R 6 R 7, D P 11, R 15, R 16, R 17, C) alimenté selon une diagonale au travers d'un élément d'alimentation (T 2, T 3 - T 12, T 13) commandé par un comparateur (A 01, A 011) disposé selon l'autre diagonale dudit pont résistif. 12 Dispositif selon la revendication 11, carac- térisé en ce que l'élément régulateur ( 11, 13) comporte en outre un élément d'amortissement des oscillations ther- miques du filament (R 8, C 1 R 18, C 11). 13 Dispositif selon l'une quelconque des reven- dications 9 à 12, caractérisé en ce que l'élément régula- teur ( 11, 13) est commandé au travers d'un élément d'allu- mage (T 1, T 4) depuis l'unité de traitement numérique. 14 Dispositif selon l'une quelconque des reven- dications 9 à 13, caractérisé en ce que l'unité de traite- ment numérique ( 2) admet plusieurs options de calcul donnant une estimation du taux d'explosivité du milieu gazeux en fonction des paramètres U O et k. Dispositif selon l'une quelconque des re- vendications 9 à 14, destiné au dosage d'hydrocarbures et d'hydrogène, caractérisé en ce que l'unité de traitement numérique contient en mémoire des données obtenues par étalonnage avec du méthane.