Brevet FR3110701A1 — publié le 20211126

Méthode de calcul d’un flux d’au moins un premier gaz émis par une source dans l’atmosphère, utilisant un deuxième gaz traceur, procédé, système et nécessaire associés La méthode comporte les étapes suivantes : - obtention de premières données représentatives de teneurs en premier gaz, et de deuxièmes données représentatives de teneurs en un deuxième gaz traceur émis par la source (14) conjointement avec le premier gaz, - calcul d’au moins un coefficient de corrélation entre les teneurs en premier gaz et les teneurs en deuxième gaz à partir des premières données représentatives et des deuxièmes données représentatives ; - obtention d’un flux mesuré ou calculé de deuxième gaz émis par la source (14) ; - calcul d’un flux de premier gaz émis par la source (14), sur la base du flux mesuré ou calculé de deuxième gaz émis par la source (14) et du coefficient de corrélation. Figure pour l'abrégé : Figure 3 Méthode de calcul d’un flux d’au moins un premier gaz émis par une source dans l’atmosphère, utilisant un deuxième gaz traceur, procédé, système et nécessaire associés La présente invention concerne une méthode de calcul d’un flux d’au moins un premier gaz émis par une source dans l’atmosphère, mis en œuvre par un système de calcul, comprenant les étapes suivantes : - obtention de premières données représentatives de teneurs en premier gaz, mesurées dans l’atmosphère à l’écart de la source, suivant une trajectoire. Le premier gaz est notamment un gaz à effet de serre émis par la source lors d’une réaction chimique. La réaction chimique est par exemple une réaction de combustion d’un flux de méthane produisant notamment du dioxyde de carbone. Le premier gaz dont le flux est calculé est par exemple le méthane résiduel n’ayant pas subi de combustion. Les préoccupations sur la protection de l’environnement ont contribué au renforcement de la législation sur les émissions polluantes, notamment en Europe. De ce fait, les unités industrielles, telles que celles présentes dans l’industrie pétrolière ou chimique doivent s’adapter à des contraintes environnementales de plus en plus exigeantes. En particulier, des gaz à effet de serre sont émis lors d’opérations d’exploitation, de transport, de raffinage, et de dépôt des hydrocarbures. Ces émissions sont suivies par les exploitants et font régulièrement l’objet de mesures de réduction. Il est en particulier nécessaire de caractériser quelles sont les sources de gaz à effet de serre et quelles sont les quantités émises par ces sources pour s’assurer de leur contrôle et reporter les progrès effectués. Toutefois, l’identification des sources d’émission de gaz à effet de serre et la quantification des émissions diffuses et fugitives ne donnent pas entière satisfaction. En effet, les émissions sont très difficiles à mesurer, car souvent non canalisées, et potentiellement près de bassins ou de lacs ou à des endroits inaccessibles, par exemple en hauteur ou en milieu d’unité. Une difficulté majeure pour évaluer les émissions d’une source ponctuelle au sein d’une installation est souvent la difficulté, voire l’impossibilité de se rendre au plus près de la source pour mesurer le flux de gaz émis par la source dans l’atmosphère. En outre, la concentration de certains gaz émis est faible dans l’atmosphère, notamment lorsqu’il s’agit de résidus de combustion. Il est donc souvent difficile de mesurer de manière suffisamment précise, et ensuite de quantifier le flux de gaz émis, compte tenu des imprécisions de mesure. Néanmoins, le suivi de la quantité des gaz émis par la source à une faible concentration est souvent utile, notamment pour répondre aux obligations réglementaires, pour des raisons de sécurité ou/et pour optimiser le fonctionnement de l’installation. Un but de l’invention est donc de disposer d’une méthode permettant de calculer le flux d’au moins un premier gaz émis par une source, notamment lorsque le premier gaz présente une concentration réduite dans l’atmosphère, la méthode étant simple à mettre en œuvre, tout en étant précise. A cet effet, l’invention a pour objet une méthode du type précité, caractérisée par les étapes suivantes : - obtention de deuxièmes données représentatives de teneurs en un deuxième gaz traceur émis par la source conjointement avec le premier gaz, les deuxièmes données représentatives étant mesurées dans l’atmosphère à l’écart de la source, suivant la trajectoire, - calcul d’au moins un coefficient de corrélation entre les teneurs en premier gaz et les teneurs en deuxième gaz à partir des premières données représentatives et des deuxièmes données représentatives ; - obtention d’un flux mesuré ou calculé de deuxième gaz émis par la source ; - calcul d’un flux de premier gaz émis par la source, sur la base du flux mesuré ou calculé de deuxième gaz émis par la source et du coefficient de corrélation. La méthode selon l’invention peut comprendre l’une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toute combinaison techniquement possible : - le coefficient de corrélation est une constante représentative d’un rapport entre les teneurs en premier gaz le long de la trajectoire et les teneurs en deuxième gaz le long de la trajectoire ; - le calcul du coefficient de corrélation comporte l’intégration des premières données représentatives le long d’au moins une partie de la trajectoire pour obtenir une première teneur globale intégrée, l’intégration des deuxièmes données représentatives le long d’au moins une deuxième partie de la trajectoire pour obtenir une deuxième teneur globale intégrée, le coefficient de corrélation étant calculé à partir de la première teneur intégrée et la deuxième teneur globale intégrée ; - elle comprend l’intégration des premières données représentatives sur une pluralité de parties de trajectoire pour obtenir une première teneur globale intégrée correspondant à chaque partie de trajectoire, puis l’intégration des deuxièmes données représentatives sur la même pluralité de parties de trajectoire, le calcul pour chaque partie de trajectoire, d’un rapport de teneurs entre la première teneur globale intégrée et la deuxième teneur globale intégrée sur la partie de trajectoire, le coefficient de corrélation étant calculé comme la moyenne des rapports de teneurs calculés ; - la trajectoire comprend une pluralité de lignes parallèles ; - le calcul du flux de premier gaz émis par la source est effectué en outre sur la base du rapport entre la masse molaire du premier gaz et la masse molaire du deuxième gaz ; - les émissions de premier gaz et de deuxième gaz issues de la source résultent d’une réaction chimique, en particulier d’une combustion, le flux mesuré ou calculé de deuxième gaz étant calculé par un bilan de matière de la réaction chimique, en particulier par un bilan de combustion ; - le deuxième gaz est produit par la réaction chimique, un troisième gaz étant produit conjointement avec le deuxième gaz par la réaction chimique, le procédé comprenant l’obtention de troisièmes données représentatives de teneurs en troisième gaz, le flux mesuré ou calculé de deuxième gaz étant calculé en tenant compte des troisièmes données représentatives ; - le deuxième gaz est le dioxyde de carbone ; - le premier gaz est le méthane, le benzène, ou/et un composé organique volatile ; - le flux est calculé par l’équation Q1 = C x M1/M2 x Q2, où C est le coefficient de corrélation, M1 est la masse molaire du premier gaz, M2 est la masse molaire du deuxième gaz, et Q2 est le flux mesuré ou calculé du deuxième gaz émis par la source. L’invention a également pour objet un procédé de mesure des émissions d’une source dans l’atmosphère, comprenant les étapes suivantes : - collecte de premières données représentatives de teneurs en au moins un premier gaz dans l’atmosphère à l’écart de la source par un drone suivant une trajectoire ; - collecte conjointe par le drone de deuxièmes données représentatives de teneurs en un deuxième gaz traceur émis par la source conjointement avec le premier gaz, suivant la même trajectoire ; - transfert des premières et deuxièmes données représentatives collectées à un système de calcul ; - mise en œuvre par le système de calcul de la méthode de calcul telle que définie plus haut. Le procédé selon l’invention peut comprendre l’une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toute combinaison techniquement possible : - il comprend une étape préalable de détermination d’une direction du vent et/ou d’une configuration d’un panache d’émission en aval de la source, le vol du drone étant effectué en fonction de la configuration de panache prédéterminée. L’invention a également pour objet un système de calcul d’un flux d’au moins un premier gaz émis par une source dans l’atmosphère, comprenant : - un module d’obtention de données, propre à obtenir des premières données représentatives de teneurs en un premier gaz, mesurées dans l’atmosphère à l’écart de la source, suivant une trajectoire, le système étant caractérisé en ce que le module d’obtention de données est propre à obtenir des deuxièmes données représentatives de teneurs en un deuxième gaz traceur émis par la source conjointement avec le premier gaz, les deuxièmes données représentatives étant mesurées dans l’atmosphère à l’écart de la source, suivant la trajectoire, le système comportant : - un module de calcul d’au moins un coefficient de corrélation entre les teneurs en premier gaz et les teneurs en deuxième gaz à partir des premières données représentatives et des deuxièmes données représentatives ; - un module d’obtention d’un flux mesuré ou calculé de deuxième gaz émis par la source ; - un module de calcul d’un flux de premier gaz émis par la source, sur la base du flux mesuré ou calculé de deuxième gaz émis par la source et du coefficient de corrélation. L’invention a aussi pour objet un nécessaire de mesure des émissions d’une source dans l’atmosphère, comportant : - un drone, apte à voler dans l’atmosphère à l’écart de la source suivant une trajectoire ; le drone étant apte à mesurer le long de la trajectoire, des premières données représentatives de teneurs en au moins un premier gaz émis par la source et des deuxièmes données représentatives de teneurs en au moins un deuxième gaz traceur émis par la source conjointement au premier gaz ; - un système de calcul tel que défini plus haut, propre à recevoir les premières et deuxièmes données représentatives mesurées par le drone. L’invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d’exemple, et faite en se référant aux dessins annexés, sur lesquels : - la figure 1 est une vue schématique d’un premier nécessaire de mesure d’émissions selon l’invention ; - la figure 2 est une vue d’une source de gaz au sein d’une installation, et du panache émis par la source de gaz. - la figure 3 est une vue détaillée du panache résultant de l’émission par la source dans des conditions de vent établi ; - la figure 4 est une vue du plan de vol mis en œuvre par le drone du nécessaire de la figure 1 ; - la figure 5 est une vue des mesures de teneurs en au moins un premier gaz, réalisées sur une ligne horizontale lors de la mise en œuvre du plan de vol de la figure 4 ; - la figure 6 est une vue des mesures de teneurs en au moins un deuxième gaz traceur, réalisées sur une ligne horizontale lors de la mise en œuvre du plan de vol de la figure 4 ; - la figure 7 est une vue analogue à la figure 3, dans le cas d’une émission à faibles conditions de vent. Un nécessaire 10 de mesure d’émissions d’au moins un premier gaz émis par une source dans l’atmosphère, utilisant un deuxième gaz traceur est illustré schématiquement sur la figure 1. Le nécessaire 10 est destiné à mettre en œuvre un procédé de mesure d’émissions d’une installation 12 industrielle, représentée schématiquement sur la figure 2. L’installation 12 comporte au moins une source 14 émettant le premier gaz dont la teneur est mesurée et émettant le deuxième gaz traceur. Le premier gaz est de préférence un gaz présentant des émissions plus faibles que le deuxième gaz traceur. Le premier gaz et le deuxième gaz sont par exemple émis conjointement par la source 14, comme résultant d’une réaction chimique, en particulier d’une réaction de combustion. Par exemple, le deuxième gaz est un produit d’une réaction chimique utilisant le premier gaz comme réactif. Le premier gaz est alors un réactif résiduel n’ayant pas réagi lors de la mise en œuvre de la réaction chimique produisant le premier gaz. Avantageusement, comme illustré sur la figure 3, la source 14 est une torche mettant en œuvre une combustion d’un flux de méthane, qui constitue le premier gaz. Le deuxième gaz traceur est le dioxyde de carbone produit par la combustion du méthane au sein de la torche. Dans des variantes, d’autres gaz sont mesurables, tel que les gaz aromatiques, notamment le benzène ou encore le 1,3 Butadiène, le monoxyde de carbone, l’éthane et plus généralement les Composés Organiques Volatils. L’installation industrielle 12 est notamment une installation pétrolière, en particulier une installation d’exploitation, de transport, de raffinage, de traitement ou de dépôt d’hydrocarbures située en mer ou à terre. Dans l’exemple représenté sur la figure 3, la source 14 émet des gaz suivant un panache 16 qui se dégage à partir de la source 14 et qui se propage sous l’effet du vent V. Le panache 16 est entraîné par le vent V circulant dans l’atmosphère au voisinage de la source 14. Il présente avantageusement une zone 18 de surélévation du panache 16 sensiblement verticale et une zone 20 de propagation du panache, qui dans cet exemple est sensiblement horizontale. Dans l’exemple de la figure 7, dans le cas où le vent V est plus faible, la zone de propagation 18 est plus haute et la zone de propagation 20 s’étend de manière inclinée par rapport à l’horizontale. La teneur moyenne en premier gaz dans le panache 16 est inférieure à la teneur moyenne en deuxième gaz dans le panache 16, par exemple d’au moins un facteur 10, voire d’au moins un facteur 100. En référence à la figure 1, pour mettre en œuvre le procédé de mesure, le nécessaire de mesure 10 comporte un drone 22 de collecte de données représentatives de teneurs en premier gaz et en deuxième gaz, à une pluralité de positions dans l’atmosphère, à l’écart de la source 14, en suivant une trajectoire donnée 23, dont un exemple est donné sur la figure 4. Le nécessaire 10 comprend en outre un système de calcul 24, propre à mettre en œuvre une méthode de calcul d’un flux du premier gaz émis par la source 14 dans l’atmosphère, à partir des données représentatives de teneurs en premier gaz et en deuxième gaz dans l’atmosphère, mesurées par le drone 22 le long de la trajectoire 23. Le drone 22 est propre à réaliser les mesures nécessaires à la collecte de données représentatives des teneurs en premier gaz et en deuxième gaz présents dans le panache 16, à l’écart de la source 14. Il comporte un châssis 30, un ensemble de propulsion 32, propre à permettre le décollage du châssis 30 au-dessus et à l’écart du sol et son déplacement en vol dans l’atmosphère au-dessus du sol. Le drone 22 comporte en outre un ensemble de mesure 34, un ensemble 36 de commande de l’ensemble de mesure 34 et de préférence, un système de télétransmission 38. En référence à la figure 1, l’ensemble de propulsion 32 comporte une pluralité d’organes de propulsion 32A, qui sont ici des hélices entraînées en rotation par un moteur. L’ensemble de propulsion 32 comporte en outre une source d’énergie 32B formée ici par une batterie et un système 33 de localisation et de commande du déplacement du drone 22 dans l’atmosphère. Dans cet exemple, le drone 22 est un drone multi rotor à voilure tournante. Il est dépourvu d’ailes, sa portance étant assurée par l’ensemble de propulsion 32. Le drone 22 est par exemple un drone quadricoptère à voilure tournante, notamment un drone DJI M200 commercialisé par la société DJI. L’ensemble de propulsion 32 comporte une pluralité d’hélices rotatives autour d’axes sensiblement verticaux. Par « sensiblement verticaux », on entend généralement que les axes de rotation des hélices sont inclinés de moins de 30° par rapport à la verticale. Lorsque les moteurs des hélices sont alimentés électriquement par la batterie, les hélices sont entraînées en rotation autour de leur axe, entraînant un flux d’air descendant. Le système de localisation et de commande 33 comporte un capteur de position, notamment un GPS et/ou une centrale inertielle. Il comprend en outre une unité de commande, propre à piloter le déplacement du drone 22 le long de la trajectoire 23, pré-enregistrée avant le vol et chargée dans le système 33, ou de manière déportée et manuelle via une télécommande déportée. Le drone 22 est ainsi propre à suivre automatiquement la trajectoire prédéfinie 23, ou alternativement, à être piloté manuellement par un opérateur, pour réaliser la trajectoire prédéfinie 23 établie dans un plan de vol. De préférence, pour la mise en œuvre du procédé de mesure, le drone 22 est propre à effectuer une trajectoire 23 suivant un déplacement en échelle rampante, telle qu’illustré par la figure 4. Le drone 22 se déplace suivant une pluralité de lignes 50 parallèles à une première direction D1, avec un segment de connexion 52 entre chaque paire de lignes parallèles 50 adjacentes. Le segment de connexion 52 s’effectue suivant une deuxième direction D2 transverse à la première direction D1. Ici, la première direction D1 est une direction horizontale et la deuxième direction D2 est une direction verticale. Dans cet exemple, toutes les lignes parallèles 50 balayées par le drone 22 s’étendent sensiblement dans un même plan Pm de mesure vertical. L’étendue E1 des lignes 50 suivant la première direction D1 est choisie en fonction de la largeur du panache 16, pour balayer tout le panache 16. Cette étendue E1 est généralement supérieure à 20 m et est comprise entre 20 m et 500 m. La distance entre les lignes 50 est définie par une étendue E2 des segments de connexion 52 suivant la deuxième direction. Cette étendue E2 est par exemple supérieure à 1 m et notamment compris entre 1 m et 50 m. L’ensemble de mesure 34 comporte des capteurs propres à réaliser des mesures de données représentatives des teneurs en premier gaz et en deuxième gaz présents dans l’atmosphère, en une pluralité de points le long de chaque ligne 50. De préférence, les données représentatives des teneurs en premier gaz et en deuxième gaz sont collectées par l’ensemble de mesure 34 le long de chaque ligne 50. Les mesures sont réalisées en continu le long de la ligne 50. La fréquence de mesure de données représentatives en chaque teneur en gaz est par exemple supérieure à 1Hz et est notamment comprise 1 Hz et 100 Hz. Un exemple d’ensemble de mesure 34 est décrit dans la demande n° 20 03027 de la Demanderesse, déposée le 27 mars 2020 à l’Institut National de la Propriété Industrielle en France, dont le titre est « Drone de mesure de données représentatives de teneurs en au moins deux gaz présents dans l’atmosphère à l’écart du sol et procédé de mesure associé ». Le système de commande 33 comporte une unité de recueil de données qui comprend au moins une mémoire propre à stocker les données représentatives de chaque teneur en chaque gaz, associée à la position géographique le long de chaque ligne 50. L’unité de recueil de données est raccordée au système de télétransmission 38 pour permettre l’export des données vers le système de calcul 24, lors du vol du drone ou après le vol du drone. Le système de calcul 24 est ici situé au sol. Il comporte au moins un calculateur 60 et une interface homme machine comprenant un organe de commande 62 tel qu’un clavier, une souris et/ou un écran tactile, l’interface homme machine comportant aussi un afficheur 64, en particulier un écran. Le calculateur 60 comprend de manière connue au moins un processeur 66 et une mémoire 68 comportant des modules logiciels propres à être exécutés par le processeur 66 pour réaliser des fonctions. En variante, le calculateur 60 comporte des composants logiques programmables ou des circuits intégrés dédiés, destinés à réaliser les fonctions des modules qui seront décrits ci-après. En référence à la figure 1, la mémoire 68 contient un module 70 d’obtention et de traitement initial de premières données représentatives de teneurs T1 en premier gaz et de deuxièmes données représentatives de teneurs T2 en deuxième gaz, le long de chaque ligne 50 de la trajectoire 23. La mémoire 68 contient en outre un module 72 de calcul d’un coefficient de corrélation C entre les teneurs T1 en premier gaz et les teneurs T2 en deuxième gaz, à partir des premières données représentatives et des deuxièmes données représentatives. La mémoire 68 contient également un module d’obtention 74 d’un flux Q2 mesuré ou calculé de deuxième gaz émis par la source 14 et un module 76 de calcul d’un flux Q1 de premier gaz émis par la source, sur la base du flux mesuré ou calculé Q2 obtenu par le module d’obtention 74 et du coefficient de corrélation C. Le module d’obtention et de traitement 70 est propre à recevoir les premières données représentatives de teneurs T1 en premier gaz et les deuxièmes données représentatives de teneurs T2 en deuxième gaz, le long de chaque ligne 50, telles que mesurées par le drone 22 à chaque point de mesure, associées à la position géographique X du point de mesure le long de la ligne 50, chaque ligne 50 constituant une partie de la trajectoire 23. Il est propre à transformer les données représentatives mesurées, en teneurs T1, T2 en chacun des gaz, à chaque point de mesure X sur chaque ligne 50, sur la base d’une courbe de calibration associée à chaque gaz. Une courbe 71 de teneurs T1 en premier gaz en fonction d’une première coordonnée X le long de la ligne 50 suivant la direction D1 est ainsi obtenue, comme illustré sur la figure 5. Une courbe 75 de teneurs T2 en deuxième gaz en fonction d’une première coordonnée X le long de la ligne 50 suivant la direction D1 est ainsi obtenue, comme illustré sur la figure 6. Le module d’obtention et de traitement 70 est en outre propre à éventuellement filtrer les teneurs obtenues. Selon une première méthode, le module d’obtention et de traitement 70 est propre à détecter des pics 71A, 75A de teneur sur chaque courbe 71, 75 sur la base d’un seuil S prédéterminé d’occurrence d’un pic, puis à supprimer les pics 71A, 75A observés de la courbe obtenue pour obtenir une courbe de valeurs de fond en fonction de la première coordonnée X. Dans une variante, le module d’obtention et de traitement 70 est propre à mettre en œuvre un algorithme itératif dans lequel la valeur moyenne des teneurs le long de la ligne 50 est calculée, puis dans lequel toutes les teneurs se situant au-dessus de la valeur moyenne sont supprimées de la courbe 71, 75 puis à répéter les étapes de calcul de la valeur moyenne et de soustraction des teneurs se situant au-dessus de la valeur moyenne jusqu’à ce qu’un critère de convergence soit atteint. Le critère de convergence est par exemple que la différence entre les valeurs moyennes successives entre deux itérations soit inférieure à une valeur prédéterminée par exemple inférieure à 10%. Ainsi, un fond continu est déterminé et est soustrait de la courbe 71, 75 représentant les teneurs T en fonction de la position X sur chaque ligne 50. Le module de calcul 72 est propre à intégrer chaque courbe 71, 75 représentant les teneurs en chaque gaz le long de chaque ligne 50, suivant la première direction D1, sur toute la largeur de la ligne 50 pour obtenir une teneur globale intégrée TGI1, TGI2 sur chaque ligne 50, par les équations suivantes : où X min et X max sont les coordonnées géographiques caractérisant les bornes du panaches définies suivant l’étendue E1 du panache parallèlement à la première direction D1. Selon la première méthode de traitement des données par le module 70, l’intégrale de la courbe des valeurs de fond est également calculée et est soustraite de l’intégrale précédente. Selon la deuxième méthode, la courbe de valeurs de fond est soustraite de la courbe 70 des teneurs, avant intégration. Ainsi, pour chaque ligne 50 dans laquelle une mesure est intervenue, correspondant à une coordonnée Z suivant la deuxième direction D2, des teneurs globales intégrées TGI1 (Z), TGI2 (Z) correspondant à chaque gaz sont obtenues. Le module de calcul 72 est en outre propre à calculer, pour chaque ligne 50 dans laquelle une mesure est intervenue, correspondant à une coordonnée Z suivant la deuxième direction D2, le rapport R(Z) de la teneur globale intégrée TGI1(Z) en premier gaz sur la teneur globale intégrée TGI2(Z) en deuxième gaz. Le module de calcul 72 est également propre à calculer le coefficient de corrélation C comme étant égal à la moyenne, de préférence arithmétique, des rapports de teneurs R(Z) obtenus sur chaque ligne 50 définissant une partie de la trajectoire 23. Le coefficient de corrélation C correspond ainsi au rapport moyen des teneurs du premier gaz sur les teneurs en deuxième gaz sur chacune des parties de la trajectoire 23. Le module d’obtention 74 d’un flux mesuré ou calculé du deuxième gaz émis par la source est propre à obtenir ou à calculer un flux Q2 de deuxième gaz émis par la source 14 soit à partir d’une mesure de flux de deuxième gaz au niveau de la source 14, soit de préférence, à partir d’un calcul d’émission de deuxième gaz par la source, lorsque la source 14 émet le deuxième gaz lors d’une réaction chimique. Le calcul est de préférence réalisé par bilan de matière, sur la base des flux de réactifs fournis pour réaliser la réaction chimique, en supposant un rendement prédéterminé de la réaction chimique. De préférence, la réaction chimique est une combustion. Le deuxième gaz est alors du dioxyde de carbone émis par la source lors de la combustion. Un bilan de combustion est alors réalisé en fonction des flux d’alimentation dans l’installation 12 pour réaliser la combustion, en tenant compte d’un rendement de combustion prédéfini. Avantageusement, le bilan de combustion est réalisé en négligeant les émissions d’autres gaz produits par la réaction que le deuxième gaz, notamment les émissions de monoxyde de carbone dans le cas où le deuxième gaz est du dioxyde de carbone. Le calcul par bilan de matière est réalisé directement par le module d’obtention 74, ou est chargé à partir d’un système de calcul externe 77. Avantageusement, le bilan de matière, en particulier le bilan de combustion est réalisé à l’aide d’une feuille de calcul d’un tableur informatique, sur la base de la composition molaire et des flux (volumiques ou massiques) de réactifs introduits dans le réacteur permettent de réaliser la réaction. Un flux massique Q2 de deuxième gaz mesuré ou calculé est ainsi obtenu à l’aide du module d’obtention 74. Le module de calcul 76 du flux de premier gaz émis Q1 par la source est propre à calculer le flux massique Q1 en premier gaz émis par la source, en fonction du flux massique Q2 mesuré ou calculé par le module d’obtention 74, du rapport de la masse molaire M1 du premier gaz à la masse molaire M2 du deuxième gaz, et du coefficient de corrélation C déterminé par le module de calcul 72, à partir des mesures de teneurs en premier gaz et en deuxième gaz réalisées par le drone 22. De préférence, le flux de premier gaz émis par la source Q1 est obtenu par l’équation suivante : Q1 = C x M1/M2 x Q2 Eventuellement, ce calcul est répété pour une pluralité de vols du drone 22, chacun correspondant à une campagne de mesure distincte réalisée lors de l’émission de la source 14, puis le flux de premier gaz Q1 émis par la source 14 est moyenné sur l’ensemble des vols réalisés. Un procédé de mesure va maintenant être décrit. Initialement, le drone 22 est mis en vol pour effectuer une trajectoire 23 suivant un déplacement en échelle rampante dans un plan de mesure Pm, tel qu’illustré par la figure 4. Comme indiqué plus haut, le drone 22 se déplace suivant une pluralité de lignes 50 parallèles à une première direction D1, constituant des parties de la trajectoire 23, avec un segment de connexion 52 entre chaque paire de lignes parallèles 50 adjacentes, le segment de connexion 52 s’effectuant suivant une deuxième direction D2 transverse à la première direction D1. Les premières données représentatives de teneurs T1 en premier gaz et les deuxièmes données représentatives de teneurs T2 en deuxième gaz, sont collectées par l’ensemble de mesure 34 le long de chaque ligne 50. Les mesures sont réalisées en continu le long de la ligne 50. La mémoire de l’unité de recueil de données stocke les données représentatives de chaque teneur T1, T2 en chaque gaz, associées à la position géographique X le long de chaque ligne 50. Puis, lors du vol du drone 22 ou après le vol du drone 22, le système de télétransmission 38 exporte des données vers le système de calcul 24 au sol. Le module d’obtention et de traitement 70 reçoit les premières données représentatives de teneurs T1 en premier gaz et les deuxièmes données représentatives de teneurs T2 en deuxième gaz, le long de chaque ligne 50, telles que mesurées par le drone 22 à chaque point de mesure, associées à la position géographique X du point de mesure le long de la ligne 50, chaque ligne 50 constituant une partie de la trajectoire 23. Il transforme les données représentatives en teneurs T1, T2 en chacun des gaz, à chaque point de mesure sur chaque ligne 50, sur la base d’une courbe de calibration associée à chaque gaz. Pour chaque ligne 50, une courbe 71, 75 de teneur en chaque gaz en fonction d’une première coordonnée X le long de la ligne 50 suivant la direction D1 est ainsi obtenue, comme visible sur les figures 5 et 6. Le module d’obtention et de traitement 70 filtre éventuellement les teneurs obtenues par exemple selon la première méthode ou la deuxième méthode décrite plus haut. Puis, le module d’intégration 72 intègre chaque courbe 71, 75 représentant les teneurs en chaque gaz le long de chaque ligne 50, suivant la première direction D1, sur toute la largeur de la ligne 50 pour obtenir une teneur globale intégrée TGI1, TGI2 sur chaque ligne 50, par les équations suivantes : où X min et X max sont les coordonnées géographiques caractérisant les bornes du panaches définies suivant l’étendue E1 du panache parallèlement à la première direction D1. Selon la première méthode de traitement des données par le module 70, l’intégrale de la courbe des valeurs de fond est également calculée et est soustraite de l’intégrale précédente. Selon la deuxième méthode, la courbe de valeurs de fond est soustraite de la courbe 70 des teneurs, avant intégration. Ainsi pour chaque ligne 50 dans laquelle une mesure est intervenue, correspondant à une coordonnée Z suivant la deuxième direction D2, des teneurs globales intégrées TGI1 (Z), TGI2 (Z) correspondant à chaque gaz sont obtenues. Puis, le module de calcul 72 calcule, pour chaque ligne 50 dans laquelle une mesure est intervenue, correspondant à une coordonnée Z suivant la deuxième direction D2, le rapport R(Z) de la teneur globale intégrée TGI1(Z) en premier gaz sur la teneur globale intégrée TGI2(Z) en deuxième gaz. Le module de calcul 72 calcule ensuite le coefficient de corrélation C comme étant égal à la moyenne, de préférence arithmétique, des rapports de teneurs R(Z) obtenus sur chaque ligne 50 définissant une partie de la trajectoire 23. Le module d’obtention 74 d’un flux mesuré ou calculé du deuxième gaz émis par la source obtient ou calcule alors un flux Q2 de deuxième gaz émis par la source 14 soit à partir d’une mesure de flux de deuxième gaz au niveau de la source 14, soit de préférence à partir d’un calcul d’émission de deuxième gaz par la source, lorsque la source 14 émet le deuxième gaz lors d’une réaction chimique, comme décrit plus haut. Puis, le module de calcul 76 du flux de premier gaz émis Q1 par la source calcule le flux massique Q1 en premier gaz émis par la source, en fonction du flux massique Q2 mesuré ou calculé par le module d’obtention 74, du rapport de la masse molaire M1 du premier gaz à la masse molaire M2 du deuxième gaz, et du coefficient de corrélation C déterminé par le module de calcul 72, à partir des mesures de teneurs en premier gaz et en deuxième gaz réalisées par le drone 22. De préférence, le flux de premier gaz émis par la source Q1 est obtenu par l’équation suivante : Q1 = C x M1/M2 x Q2. Comme indiqué plus haut, ce calcul est éventuellement répété pour une pluralité de vols du drone 22, chacun correspondant à une campagne de mesure distincte réalisée lors de l’émission de la source 14, puis le flux de premier gaz Q1 émis par la source 14 est moyenné sur l’ensemble des vols réalisés. Le procédé de mesure selon l’invention est donc particulièrement simple à mettre en œuvre, puisqu’il nécessite une simple campagne de mesures à l’aide d’un drone 22 volant directement dans le panache 16, à l’écart de la source 14. Suite à cette campagne de mesures, un calcul simple et efficace pour obtenir une détermination précise du flux d’un premier gaz émis par la source 14, même si les teneurs en premier gaz mesurées dans le panache 16 sont faibles, grâce à la mesure simultanée de teneurs en un deuxième gaz traceur, dont le flux émis peut être mesuré ou calculé par bilan de matière. Ceci est particulièrement utile pour détecter des composés dont la teneur est faible dans l’atmosphère. Par exemple, le tableau ci-dessus illustre les teneurs en méthane émis par une source 14 constituée par une torche dans une installation pétrolière. Le débit total de gaz de torche et le débit de méthane permettent de calculer le flux Q2 de dioxyde de carbone en effectuant un bilan de combustion. La méthode selon l’invention est utilisée pour mesurer un coefficient de corrélation C entre le méthane et le dioxyde de carbone, lequel est utilisé pour déterminer le flux Q1 de méthane émis à partir du flux Q2 de dioxyde de carbone calculé par bilan de matière. Ceci permet de déterminer, à partir d’un débit moyen de méthane torché de 0,23 tonne/heure, une quantité de méthane résiduel égale à 0,7 g/h et donc une efficacité de 99% du point de vue de la destruction du méthane. Vol Débit gaz de torche (t/h) Débit de méthane fourni (t/h) Flux de dioxyde de carbone Q2 calculé par bilan de combustion (g/s) Coefficient de corrélation C Flux Q1 de méthane émis (g/s) 1 5,6 0,49 2286 0,00078 0,7 2 5,3 0,42 2121 0,00073 0,6 3 5,3 0,37 2056 0,00044 0,3 4 5,6 0,26 2040 0,00120 0,9 5 5,4 0,22 1949 0,00225 1,6 6 5,2 0,13 1883 0,00095 0,7 7 5,3 0,13 1895 0,00083 0,6 8 5,4 0,13 1911 0,00122 0,9 9 5,6 0,13 1959 0,00076 0,5 10 5,6 0,12 1911 0,00081 0,6 11 5,4 0,12 1741 0,00145 0,9 Moyenne 5,4 0,23 1977 0,00104 0,7 La méthode selon l’invention est apte à être mise en œuvre au voisinage d’installations industrielles 12 variées, même si ces installations sont inaccessibles ou/et nécessitent des précautions de sécurité. Les mesures peuvent être effectuées à moindre coût et de manière fréquente, ce qui permet notamment de suivre l’évolution des émissions engendrées par la source 14, et de s’assurer qu’elles sont sous contrôle ou qu’elles sont réduites. Dans une variante, représentée sur la figure 8, pour des conditions de vent peu fort, le panache 16 présente une configuration verticale et le plan de mesure Pm est un plan horizontal. La première direction D1 est alors une direction horizontale, et la deuxième direction D2 est une direction horizontale perpendiculaire à la première direction D1. Le procédé de mesure, incluant la méthode de mesure restent similaires à ceux décrits plus haut. Avantageusement, le procédé de mesure selon l’invention comporte une étape initiale de détermination de caractéristiques de conformation du panache 16, par exemple par mesure de la rose des vents s’appliquant sur la source 14 au moment de la campagne de mesure. Dans une variante, les teneurs en monoxyde de carbone sont mesurées dans le panache 16, par le drone 22. Ces teneurs sont intégrées comme décrit précédemment, et un coefficient de corrélation additionnel entre la teneur globale intégrée en monoxyde de carbone et la teneur globale intégrée en dioxyde de carbone est obtenu, pour corriger le flux Q2 de dioxyde de carbone calculé par bilan de combustion. Dans une autre variante, un coefficient de corrélation C(i) est calculé pour chaque partie de trajectoire 23, par exemple pour chaque ligne 50. Ensuite, un flux Q1(i) de premier gaz émis correspondant à chaque partie de trajectoire 23, notamment à chaque ligne 50 est obtenu par l’équation suivante : Q1(i) = C(i) x M1/M2 x Q2. Lorsque le flux Q1 de premier gaz émis par la source est stable au cours du temps il est alors possible de moyenner, notamment par une moyenne arithmétique, les flux Q1(i) de premier gaz émis correspondant à chaque partie de trajectoire 23, notamment à chaque ligne 50, pour obtenir un flux Q1 moyen. Au contraire, lorsque le flux Q1 de premier gaz évolue au cours du temps, il est possible de suivre l’évolution du flux Q1 au cours du temps en suivant les différents flux Q1(i) obtenus au cours du temps pour des parties successives de trajectoire 23 ou pour des vols successifs. Dans encore une autre variante, la trajectoire suivie par le drone 22 est différente de la trajectoire représentée. Par exemple, la trajectoire comprend une pluralité de lignes curvilignes et non rectilignes (par exemple en arc de cercle) parallèles les unes aux autres reliées par des segments de connexion. Ceci est avantageusement le cas lorsque les lignes sont horizontales. En variante encore, la trajectoire suivie par le drone 22 comporte par exemple une trajectoire en escargot, ou encore en zigzag permettant de couper le panache dans différentes sections perpendiculaires ou non à sa direction principale. Méthode de calcul d’un flux (Q1) d’au moins un premier gaz émis par une source (14) dans l’atmosphère, mis en œuvre par un système de calcul (24), comprenant les étapes suivantes : - obtention de premières données représentatives de teneurs (T1) en premier gaz, mesurées dans l’atmosphère à l’écart de la source (14), suivant une trajectoire (23), la méthode étant caractérisée par les étapes suivantes : - obtention de deuxièmes données représentatives de teneurs (T2) en un deuxième gaz traceur émis par la source (14) conjointement avec le premier gaz, les deuxièmes données représentatives étant mesurées dans l’atmosphère à l’écart de la source (14), suivant la trajectoire (23), - calcul d’au moins un coefficient de corrélation (C) entre les teneurs (T1) en premier gaz et les teneurs (T2) en deuxième gaz à partir des premières données représentatives et des deuxièmes données représentatives ; - obtention d’un flux mesuré ou calculé (Q2) de deuxième gaz émis par la source (14) ; - calcul d’un flux (Q1) de premier gaz émis par la source (14), sur la base du flux mesuré ou calculé (Q2) de deuxième gaz émis par la source (14) et du coefficient de corrélation (C). Méthode selon la revendication 1, dans laquelle le coefficient de corrélation (C) est une constante représentative d’un rapport entre les teneurs (T1) en premier gaz le long de la trajectoire (23) et les teneurs (T2) en deuxième gaz le long de la trajectoire (23). Méthode selon la revendication 2, dans laquelle le calcul du coefficient de corrélation (C) comporte l’intégration des premières données représentatives le long d’au moins une partie de la trajectoire (23) pour obtenir une première teneur globale intégrée, l’intégration des deuxièmes données représentatives le long d’au moins une deuxième partie de la trajectoire (23) pour obtenir une deuxième teneur globale intégrée, le coefficient de corrélation étant calculé à partir de la première teneur intégrée et la deuxième teneur globale intégrée. Méthode selon la revendication 3, comprenant l’intégration des premières données représentatives sur une pluralité de parties de trajectoire (23) pour obtenir une première teneur globale intégrée correspondant à chaque partie de trajectoire (23), puis l’intégration des deuxièmes données représentatives sur la même pluralité de parties de trajectoire (23), le calcul pour chaque partie de trajectoire, d’un rapport de teneurs entre la première teneur globale intégrée et la deuxième teneur globale intégrée sur la partie de trajectoire, le coefficient de corrélation (C) étant calculé comme la moyenne des rapports de teneurs calculés. Méthode selon l’une quelconque des revendications précédentes, dans laquelle la trajectoire (23) comprend une pluralité de lignes (50) parallèles. Méthode selon l’une quelconque des revendications précédentes, dans laquelle le calcul du flux (Q1) de premier gaz émis par la source (14) est effectué en outre sur la base du rapport entre la masse molaire (M1) du premier gaz et la masse molaire (M2) du deuxième gaz. Méthode selon l’une quelconque des revendications précédentes, dans laquelle les émissions de premier gaz et de deuxième gaz issues de la source (14) résultent d’une réaction chimique, en particulier d’une combustion, le flux mesuré ou calculé (Q2) de deuxième gaz étant calculé par un bilan de matière de la réaction chimique, en particulier par un bilan de combustion. Méthode selon la revendication 7, dans lequel le deuxième gaz est produit par la réaction chimique, un troisième gaz étant produit conjointement avec le deuxième gaz par la réaction chimique, le procédé comprenant l’obtention de troisièmes données représentatives de teneurs en troisième gaz, le flux mesuré ou calculé (Q2) de deuxième gaz étant calculé en tenant compte des troisièmes données représentatives. Méthode selon l’une quelconque des revendications précédentes, dans laquelle le deuxième gaz est le dioxyde de carbone. Méthode selon l’une quelconque des revendications précédentes, dans laquelle le premier gaz est le méthane, le benzène, ou/et un composé organique volatile. Méthode selon l’une quelconque des revendications précédentes, dans laquelle le flux (Q1) est calculé par l’équation Q1 = C x M1/M2 x Q2, où C est le coefficient de corrélation (C), M1 est la masse molaire du premier gaz, M2 est la masse molaire du deuxième gaz, et Q2 est le flux (Q2) mesuré ou calculé du deuxième gaz émis par la source (14). Procédé de mesure des émissions d’une source (14) dans l’atmosphère, comprenant les étapes suivantes : - collecte de premières données représentatives de teneurs en au moins un premier gaz dans l’atmosphère à l’écart de la source (14) par un drone (22) suivant une trajectoire (23) ; - collecte conjointe par le drone (22) de deuxièmes données représentatives de teneurs en un deuxième gaz traceur émis par la source (14) conjointement avec le premier gaz, suivant la même trajectoire (23) ; - transfert des premières et deuxièmes données représentatives collectées à un système de calcul (24); - mise en œuvre par le système de calcul (24) de la méthode de calcul selon l’une quelconque des revendications précédentes. Procédé selon la revendication 12, comprenant une étape préalable de détermination d’une direction du vent et/ou d’une configuration d’un panache d’émission en aval de la source (14), le vol du drone étant effectué en fonction de la configuration de panache prédéterminée. Système (24) de calcul d’un flux (Q1) d’au moins un premier gaz émis par une source (14) dans l’atmosphère, comprenant : - un module (70) d’obtention de données, propre à obtenir des premières données représentatives de teneurs (T1) en un premier gaz, mesurées dans l’atmosphère à l’écart de la source (14), suivant une trajectoire (23), le système étant caractérisé en ce que le module (70) d’obtention de données est propre à obtenir des deuxièmes données représentatives de teneurs (T2) en un deuxième gaz traceur émis par la source (14) conjointement avec le premier gaz, les deuxièmes données représentatives étant mesurées dans l’atmosphère à l’écart de la source (14), suivant la trajectoire (23), le système (24) comportant : - un module (72) de calcul d’au moins un coefficient de corrélation (C) entre les teneurs (T1) en premier gaz et les teneurs (T2) en deuxième gaz à partir des premières données représentatives et des deuxièmes données représentatives ; - un module (74) d’obtention d’un flux mesuré ou calculé (Q2) de deuxième gaz émis par la source (14) ; - un module (76) de calcul d’un flux (Q1) de premier gaz émis par la source (14), sur la base du flux mesuré ou calculé (Q2) de deuxième gaz émis par la source (14) et du coefficient de corrélation (C). Nécessaire (10) de mesure des émissions d’une source (14) dans l’atmosphère, comportant : - un drone (22), apte à voler dans l’atmosphère à l’écart de la source (14) suivant une trajectoire (23) ; le drone (22) étant apte à mesurer le long de la trajectoire (23), des premières données représentatives de teneurs en au moins un premier gaz émis par la source (24) et des deuxièmes données représentatives de teneurs en au moins un deuxième gaz traceur émis par la source (14) conjointement au premier gaz ; - un système de calcul (24) selon la revendication 14, propre à recevoir les premières et deuxièmes données représentatives mesurées par le drone (22).