FR 3109217 A1 20211015 FR 2003527 A 20200408 Méthode de calcul d’un flux d’au moins un gaz émis par une source dans l’atmosphère, procédé de mesure, système et nécessaire associés Cette méthode comporte les étapes suivantes : - récupération de données de teneurs en au moins un gaz, mesurées dans l’atmosphère à l’écart de la source (14) suivant une pluralité de lignes parallèles à une première direction ; - intégration des teneurs relevées sur chaque ligne suivant la première direction pour obtenir une teneur globale intégrée sur chaque ligne ; - intégration du produit des teneurs globales intégrées sur chaque ligne par une vitesse de vent présent sur la ligne suivant une deuxième direction perpendiculaire à la première direction, pour obtenir un flux brut de gaz ; - détermination du flux de gaz émis par la source (14), à partir du flux brut de gaz. Figure pour l'abrégé : figure 2 Méthode de calcul d’un flux d’au moins un gaz émis par une source dans l’atmosphère, procédé de mesure, système et nécessaire associés La présente invention concerne une méthode de calcul d’un flux d’au moins un gaz émis par une source dans l’atmosphère, mis en œuvre par un système de calcul. Les gaz à mesurer sont notamment des gaz à effet de serre tel que le méthane ou le dioxyde de carbone. 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, compte tenu du vent, le flux de gaz engendré par la source se disperse et se propage dans l’atmosphère sous forme d’un panache. Il est donc généralement difficile et peu précis de mesurer des émissions émises par une source ponctuelle lorsqu’on s’éloigne de la source. Un but de l’invention est de disposer d’une méthode permettant de calculer le flux d’au moins un gaz émis par une source dans l’atmosphère, en particulier d’un gaz à effet de serre, la méthode ne nécessitant pas de données prises au plus près de la source, tout en étant précise et simple à mettre en œuvre. A cet effet, l’invention a pour objet une méthode du type précité, comportant les étapes suivantes : - récupération de données de teneurs en au moins un gaz, mesurées dans l’atmosphère à l’écart de la source suivant une pluralité de lignes parallèles à une première direction ; - intégration des teneurs relevées sur chaque ligne suivant la première direction pour obtenir une teneur globale intégrée sur chaque ligne ; - intégration du produit des teneurs globales intégrées sur chaque ligne par une vitesse de vent présent sur la ligne suivant une deuxième direction perpendiculaire à la première direction, pour obtenir un flux brut de gaz ; - détermination du flux de gaz émis par la source, à partir du flux brut de gaz. La méthode selon l’invention peut comprendre une ou plusieurs des caractéristiques suivantes, prises isolément ou suivant toute combinaison technique possible : - la méthode comprend, entre les étapes d’intégration, une étape d’interpolation d’une courbe de teneurs globales intégrées en fonction d’une coordonnée suivant la deuxième direction, à partir des teneurs globales intégrées calculées ; - l’interpolation est réalisée par une interpolation cubique, en particulier par une interpolation cubique par partie ; - la méthode comporte une étape préalable de récupération de données de vitesses de vent présent sur chaque ligne ; - la méthode comporte une étape préalable de détermination d’un vent moyen commun à toutes les lignes ; - la methode comprend, après l’étape de récupération des données de teneurs, le calcul d’un fond continu de gaz présent dans l’atmosphère et le traitement des données de teneurs pour éliminer le fond continu ; - la détermination du fond continu comprend pour chaque ligne, le calcul d’une valeur moyenne de teneurs mesurée sur la ligne, puis l’élimination des teneurs situées au-dessus de la valeur moyenne, et la répétition des étapes précédentes jusqu’à ce que la différence entre deux valeurs moyennes successives soit inférieure à un seuil de convergence ; - la première direction est horizontale, la deuxième direction étant verticale ; - les flux d’au moins deux gaz émis par la source sont calculés. 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 teneurs en au moins un gaz dans l’atmosphère à l’écart de la source suivant une pluralité de lignes parallèles à une première direction par vol d’un drone muni d’un ensemble de mesure de données représentatives des teneurs en au moins un gaz ; - transfert des 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é de mesure selon l’invention peut comprendre l’une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toute combinaison techniquement possible : - le procédé 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. - le procédé comprend une étape de mesure de la vitesse du vent. L’invention a également pour objet un système de calcul d’un flux d’au moins un gaz émis par une source dans l’atmosphère, comportant : - un module d’obtention de données de teneurs en au moins un gaz, mesurées dans l’atmosphère à l’écart de la source suivant une pluralité de lignes parallèles à une première direction ; - un module d’intégration des teneurs relevées sur chaque ligne suivant la première direction pour obtenir une teneur globale intégrée sur chaque ligne ; - un module d’intégration du produit des teneurs globales intégrées sur chaque ligne par une vitesse de vent présent sur la ligne suivant une deuxième direction perpendiculaire à la première direction, pour obtenir un flux brut de gaz ; - un module de détermination du flux de gaz émis par la source, à partir du flux brut de gaz. L’invention a aussi pour objet un nécessaire de mesure des émissions d’au moins un gaz par une source dans l’atmosphère comportant : - un drone, apte à voler dans l’atmosphère à l’écart de la source suivant une pluralité de lignes parallèles à une première direction ; le drone étant apte à mesurer des données représentatives des teneurs en au moins un gaz suivant chaque ligne parallèle à la première direction ; - un système de calcul tel que défini plus haut, propre à recevoir les données représentatives mesurées par le drone. Le nécessaire selon l’invention peut comprendre l’une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toute combinaison techniquement possible : - le drone est apte à mesurer en continu des données représentatives des teneurs d’au moins deux gaz présents dans l’atmosphère. 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 lesquelles : - 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 détaillée du panache résultant de l’émission par la source dans des conditions de vent établi ; - la figure 4 est 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 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 d’une courbe reliant les teneurs intégrées obtenues à partir de plusieurs lignes en fonction de l’altitude, et l’interpolation faite entre ces teneurs ; - la figure 7 est une vue d’une estimation de la vitesse du vent en fonction de l’altitude, pouvant être utilisée dans la mise en œuvre du procédé selon l’invention ; - la figure 8 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 gaz émis par une source dans l’atmosphère 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. De préférence, l’émission en au moins deux gaz présents dans l’atmosphère est mesurée par le procédé selon l’invention. Les gaz sont de préférence le méthane et le dioxyde de carbone. 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. L’installation 12 comporte au moins une source 14 émettant des gaz dont la teneur est mesurée. Dans l’exemple représenté sur la figure 3, la source 14 est une torche. Elle é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 8, dans le cas où le vent V est plus faible, la zone de surélévation 18 est plus haute et la zone de propagation 20 s’étend de manière inclinée par rapport à l’horizontale. 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 au moins un gaz, de préférence en au moins deux gaz, à une pluralité de positions dans l’atmosphère, à l’écart de la source 14. 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 ou de chaque gaz émis par la source 14 dans l’atmosphère, à partir des données représentatives de teneurs en chaque gaz dans l’atmosphère mesurées par le drone 22. Le drone 22 est propre à réaliser les mesures nécessaires à la collecte de données représentatives des teneurs en au moins un gaz présent 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 d’une trajectoire 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 une trajectoire prédéfinie, ou alternativement, à être piloté manuellement par un opérateur, pour réaliser 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 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 au moins un capteur propre à réaliser des mesures de données représentatives de teneurs en au moins un gaz présent 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 au moins deux 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 à 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 des données représentatives de teneurs en au moins un gaz, pour calculer, sur chaque ligne parallèle 50, des teneurs successives en au moins un gaz le long de la ligne parallèle 50, dont un exemple est visible sur la figure 5. La mémoire 68 contient en outre un module 72 d’intégration des teneurs sur chaque ligne 50 suivant la première direction D1, pour obtenir une teneur globale intégrée TGI sur chaque ligne 50. La mémoire 68 contient en outre un module 74 d’interpolation d’une courbe 75 de teneurs globales intégrées TGI le long de la deuxième direction D2 transversale à la première direction D1 (voir figure 6), à partir des teneurs globales intégrées TGI calculées sur chaque ligne 50 et un module 76 d’intégration du produit de la teneur globale intégrée TGI selon la première direction par une vitesse de vent V, l’intégration étant réalisée suivant la deuxième direction D2, pour obtenir un flux brut Qb de gaz circulant dans le panache 16. La mémoire 68 contient également un module 78 de détermination d’un flux Qg de gaz émis par la source 14 par correction du flux de gaz brut Qb en fonction de la conformation du panache pour obtenir. Le module d’obtention et de traitement 70 est propre à recevoir les données représentatives des teneurs mesurées en au moins un gaz, de préférence en au moins deux 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. Il est propre à transformer les données représentatives mesurées, en teneurs 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 T 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 illustré sur la figure 5. 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 de teneur sur chaque courbe 71, sur la base d’un seuil S prédéterminé d’occurrence d’un pic, puis à supprimer les pics 71A 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, 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 représentant les teneurs T en fonction de la position X sur chaque ligne 50. Le module d’intégration 72 est propre à intégrer la courbe 71 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 TGI sur chaque ligne 50, par l’équation suivante : 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, une teneur globale intégrée TGI (Z) est obtenue. Le module d’interpolation 74 est propre à interpoler, à partir des teneurs globales intégrées TGI(Z) sur chaque ligne 50, associée à leurs coordonnées Z suivant la deuxième direction, une courbe 75 continue de teneurs globales intégrées TGI en fonction de la coordonnée Z suivant la deuxième direction D2, comme illustré par la figure 6. La courbe interpolée 75 est par exemple obtenue par une interpolation cubique, en particulier une interpolation cubique par partie jusqu’à convergence. Le module d’intégration 76 est propre à intégrer le produit de la vitesse du vent V(Z) mesurée ou obtenue à chaque coordonnée Z le long de la deuxième direction D2 avec la teneur globale intégrée TGI(Z) correspondant à cette coordonnée, obtenue à partir de la courbe interpolée 75, pour obtenir un flux brut Qb passant par le plan de mesure Pm suivant la formule : Q où Zmin et Z max sont les coordonnées minimales et maximales le long de la deuxième direction D2 pour lesquelles une ligne 50 de mesures a été obtenue. Dans l’exemple illustré par la figure 2, la vitesse du vent V(Z) est prise comme une vitesse du vent moyenne constante selon la deuxième coordonnée Z. En variante, une courbe de vent V(Z) en fonction de la deuxième coordonnée suivant la deuxième direction Z est établie, tel qu’illustrée par la figure 7, et la vitesse du vent V(Z) à chaque deuxième coordonnée Z est utilisée pour calculer le produit avec la teneur globale intégrée TGI(Z) à la deuxième coordonnée Z et réaliser l’intégration. Le module d’intégration 76 est ainsi apte à obtenir un flux total brut Qb de chaque gaz mesuré passant à travers le plan de mesure Pm visible sur la figure 3. Puis, le module de détermination 78 est propre à corriger la valeur de flux total brut Qb mesuré pour tenir compte de la conformation du panache 16. Par exemple, dans le cas où le plan de mesure Pm est vertical, un angle d’inclinaison α de la direction de flux dans le panache 16 au niveau du plan de mesure Pm est calculé, en fonction d’une valeur de la hauteur de la zone de surélévation 18, et d’une forme supposée de panache dans la zone de transport 20, calculée en fonction du vent. Un flux total de gaz Qt passant à travers un plan Pp perpendiculaire à la direction de flux est alors calculé sur la base du flux brut Qb calculé par le module d’intégration 76 et de l’angle d’inclinaison α déterminé, par exemple en supposant que la section du panache est circulaire perpendiculairement au flux. Selon le principe de conservation de la masse, le flux de gaz Qg émis par la source 14 est alors égal au flux passant par le plan Pp. Un procédé de mesure va maintenant être décrit. Initialement, le drone 22 est mis en vol pour effectuer une trajectoire 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 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 données représentatives des teneurs en au moins un gaz, de préférence en au moins deux 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 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 données représentatives des teneurs mesurées en au moins un gaz, de préférence en au moins deux 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. Il transforme les données représentatives en teneurs 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 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 la figure 4. 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 la courbe 71 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 TGI sur chaque ligne 50, par l’équation exposée ci-dessus : 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, une teneur globale intégrée TGI (Z) est obtenue. Le module d’interpolation 74 interpole ensuite, à partir des teneurs globales intégrées TGI sur chaque ligne 50, associée à leurs coordonnées Z suivant la deuxième direction, une courbe 75 continue de teneurs globales intégrées TGI en fonction de la coordonnée Z suivant la deuxième direction D2, comme illustré par la figure 5. La courbe interpolée 75 est par exemple obtenue par une interpolation cubique, en particulier une interpolation cubique par partie jusqu’à convergence. Le module d’intégration 76 intègre alors le produit de la vitesse du vent V(Z) mesurée ou obtenue à chaque coordonnée Z le long de la deuxième direction D2 avec la teneur globale intégrée TGI(Z) correspondant à cette coordonnée (Z), obtenue à partir de la courbe interpolée 75, pour obtenir un flux brut Qb passant à travers le plan de mesure suivant la formule : Q Dans l’exemple illustré par la figure 2, la vitesse du vent V(Z) est prise comme une vitesse du vent moyenne constante selon la deuxième coordonnée Z. En variante, une courbe de vent V(Z) en fonction de la deuxième coordonnée suivant la deuxième direction Z est établie, tel qu’illustrée par la figure 7, et la vitesse du vent V(Z) à chaque deuxième coordonnée Z est utilisée pour calculer le produit avec la teneur globale intégrée TGI(Z) à la deuxième coordonnée Z et réaliser l’intégration. Le module d’intégration 76 obtient ainsi un flux total brut Qb de chaque gaz mesuré passant à travers le plan de mesure Pm. Puis, le module de correction 78 corrige la valeur de flux total brut Qb mesuré pour tenir compte de la conformation du panache 16, comme décrit plus haut. Par exemple, dans le cas où le plan de mesure Pm est vertical, un angle d’inclinaison α de la direction de flux dans le panache 16 au niveau du plan de mesure Pm ainsi qu’une vitesse d’évolution du panache 16 au niveau du plan de mesure Pm sont calculés en fonction d’une valeur de la hauteur de la zone de surélévation 18, et d’une forme supposée de panache dans la zone de transport 20. Un flux total de gaz Qg passant à travers un plan perpendiculaire Pp à la direction de flux est alors calculé sur la base du flux brut Qb calculé par le module d’intégration 76, de l’angle d’inclinaison α déterminé et de la vitesse calculée d’évolution du panache 16, par exemple en supposant que la section du panache est circulaire perpendiculairement au flux. La vitesse d’évolution est calculée le long de l’axe du panache et correspond à la norme du vecteur vitesse du panache qui est perpendiculaire au plan Pp. Selon le principe de conservation de la masse, le flux Qg de gaz émis par la source 14 est alors égal au flux passant par ce plan. 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, le calcul est simple et efficace pour obtenir une détermination précise du flux émis par la source 14. Cette méthode 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 autre variante, la vitesse du vent est mesurée directement au droit de l’installation 12, par exemple par un dispositif de mesure 100 de type LIDAR VENT, tel qu’illustré par la figure 2. Méthode de calcul d’un flux d’au moins un gaz émis par une source (14) dans l’atmosphère, mis en œuvre par un système de calcul (24), comportant les étapes suivantes : - récupération de données de teneurs en au moins un gaz, mesurées dans l’atmosphère à l’écart de la source (14) suivant une pluralité de lignes (50) parallèles à une première direction (D1) ; - intégration des teneurs relevées sur chaque ligne (50) suivant la première direction (D1) pour obtenir une teneur globale intégrée sur chaque ligne (50) ; - intégration du produit des teneurs globales intégrées sur chaque ligne (50) par une vitesse de vent présent sur la ligne (50) suivant une deuxième direction (D2) perpendiculaire à la première direction (D1), pour obtenir un flux brut de gaz ; - détermination du flux de gaz émis par la source (14), à partir du flux brut de gaz. Méthode selon la revendication 1, comprenant, entre les étapes d’intégration, une étape d’interpolation d’une courbe (75) de teneurs globales intégrées en fonction d’une coordonnée suivant la deuxième direction (D2), à partir des teneurs globales intégrées calculées. Méthode selon la revendication 2, dans laquelle l’interpolation est réalisée par une interpolation cubique, en particulier par une interpolation cubique par partie. Méthode selon l’une quelconque des revendications précédentes comportant une étape préalable de récupération de données de vitesses de vent présent sur chaque ligne (50). Méthode selon l’une quelconque des revendications précédentes comportant une étape préalable de détermination d’un vent moyen commun à toutes les lignes (50). Méthode selon l’une quelconque des revendications précédentes comprenant, après l’étape de récupération des données de teneurs, le calcul d’un fond continu de gaz présent dans l’atmosphère et le traitement des données de teneurs pour éliminer le fond continu. Méthode selon la revendication 6, dans lequel la détermination du fond continu comprend pour chaque ligne (50), le calcul d’une valeur moyenne de teneurs mesurée sur la ligne (50), puis l’élimination des teneurs situées au-dessus de la valeur moyenne, et la répétition des étapes précédentes jusqu’à ce que la différence entre deux valeurs moyennes successives soit inférieure à un seuil de convergence. Méthode selon l’une quelconque des revendications précédentes, dans laquelle la première direction (D1) est horizontale, la deuxième direction (D2) étant verticale. Méthode selon l’une quelconque des revendications précédentes, dans laquelle les flux d’au moins deux gaz émis par la source (14) sont calculés. Procédé de mesure des émissions d’une source (14) dans l’atmosphère comprenant les étapes suivantes : - collecte de teneurs en au moins un gaz dans l’atmosphère à l’écart de la source (14) suivant une pluralité de lignes (50) parallèles à une première direction (D1) par vol d’un drone (22) muni d’un ensemble (34) de mesure de données représentatives des teneurs en au moins un gaz ; - transfert des 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 10, 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 (22) étant effectué en fonction de la configuration de panache prédéterminée. Procédé selon l’une quelconque des revendications 10 à 11, comprenant une étape de mesure de la vitesse du vent. Système de calcul (24) d’un flux d’au moins un gaz émis par une source (14) dans l’atmosphère, comportant - un module (70) d’obtention de données de teneurs en au moins un gaz, mesurées dans l’atmosphère à l’écart de la source (14) suivant une pluralité de lignes (50) parallèles à une première direction (D1) ; - un module (72) d’intégration des teneurs relevées sur chaque ligne (50) suivant la première direction (D1) pour obtenir une teneur globale intégrée sur chaque ligne (50) ; - un module (76) d’intégration du produit des teneurs globales intégrées sur chaque ligne (50) par une vitesse de vent présent sur la ligne (50) suivant une deuxième direction (D2) perpendiculaire à la première direction (D1), pour obtenir un flux brut de gaz ; - un module (78) de détermination du flux de gaz émis par la source (14), à partir du flux brut de gaz. Nécessaire (10) de mesure des émissions d’au moins un gaz par 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 pluralité de lignes (50) parallèles à une première direction (D1) ; le drone (22) étant apte à mesurer des données représentatives des teneurs en au moins un gaz suivant chaque ligne (50) parallèle à la première direction (D1) ; - un système de calcul (24) selon la revendication 13, propre à recevoir les données représentatives mesurées par le drone (22). Nécessaire (10) selon la revendication 14, dans lequel le drone (22) est apte à mesurer en continu des données représentatives des teneurs d’au moins deux gaz présents dans l’atmosphère.