Brevet FR2749652A1 — publié le 19971212

La présente invention a pour objet un procédé de mesure du temps de propagation d'un signal acoustique dans un écoulement de fluide entre un premier transducteur jouant le rôle d'émetteur et un second transducteur jouant le rôle de récepteur et situé à une distance déterminée du premier transducteur, le signal acoustique émis par le premier transducteur étant constitué par au moins une impulsion émise à une fréquence acoustique déterminée Fa et le signal acoustique reçu par le second transducteur comprenant une série d'oscillations caractéristiques dont l'amplitude est d'abord croissante sur plusieurs périodes, puis décroissante sur plusieurs périodes suivantes, I'enveloppe des oscillations caractéristiques présentant la forme d'un fuseau, le procédé consistant à échantillonner le signal acoustique reçu, à une fréquence d'échantillonnage Fe, à numériser le signal acoustique reçu échantillonné, et à rechercher, par analyse du signal acoustique reçu échantillonné et numérisé, le premier passage à zéro significatif des oscillations caractéristiques du signal acoustique reçu. II est connu depuis de nombreuses années de mesurer le débit d'un fluide (ou son volume) s'écoulant dans une conduite en utilisant la propagation des signaux acoustiques émis entre deux transducteurs acoustiques situés en des points espacés dans la direction de l'écoulement du fluide. En principe, un signal acoustique émis du premier transducteur vers le second transducteur est reçu par ce second transducteur et le temps de propagation Td de ce signal acoustique est mesuré. De même, on mesure le temps de propagation Tu d'un signal acoustique émis du second transducteur vers le premier transducteur après réception dudit signal par ce premier transducteur. Dans un compteur de fluide, le débit peut être obtenu en combinant une mesure des temps de propagation des deux signaux acoustiques émis entre les deux points dans des directions opposées avec une mesure des déphasages acoustiques induits dans chaque signal acoustique par la propagation de chacun des signaux acoustiques dans l'écoulement. La demande de brevet européen nO 0 426 309 décrit un exemple d'un tel système de mesure de débit, dans lequel le signal reçu est échantillonné et converti numériquement, la mesure du déphasage acoustique étant effectuée en réalisant une détection synchrone du signal numérisé. Lors de la mesure de la vitesse d'écoulement d'un gaz dans un compteur de gaz utilisant deux transducteurs ultrasonores, il est nécessaire de mesurer le temps de parcours de l'onde ultrasonore entre l'instant d'émission et l'instant de réception, dès lors que la célérité des ondes ultrasonores est dépendante de la nature du gaz. On a représenté sur la Figure 2 la forme d'un signal impulsionnel rectangulaire S1 de largeur T émis à un instant To par un premier transducteur ultrasonore disposé dans le flux d'écoulement d'un fluide en un premier point, et la forme du signal S2 constituant la réponse impulsionnelle reçue en écho à un instant T1 par un deuxième transducteur ultrasonore disposé dans le flux d'écoulement du fluide en un deuxième point distinct du premier point. Le signal acoustique S2 reçu par le second transducteur est constitué par une série d'oscillations caractéristiques Oc qui augmentent d'amplitude sur plusieurs périodes puis décroissent, I'enveloppe des oscillations caractéristiques présentant la forme d'un fuseau. Les oscillations caractéristiques Oc du signal S2 sont précédées et suivies d'oscillations parasites Op de faible amplitude. Pour déterminer l'instant T1 de début des oscillations caractéristiques, il convient de repérer le premier passage à zéro significatif des oscillations caractéristiques Oc du signal acoustique reçu S2. On a représenté sur la Figure 3, à une échelle plus grande que sur la Figure 2, un exemple de signal acoustique S2 reçu en écho à une impulsion rectangulaire S1 émise à une fréquence acoustique déterminée Fa. Pour déterminer le début des oscillations caractéristiques Oc, selon un procédé connu, on fixe une tension de seuil Vs, par rapport à laquelle on compare le niveau du signal acoustique reçu S2, la comparaison s'effectuant sur un signal numérique obtenu après échantillonnage du signal acoustique analogique reçu, à une fréquence d'échantillonnage Fe qui constitue par exemple un multiple de la fréquence acoustique Fa. Dans ce cas, on repère l'instant T2 auquel l'amplitude du signal reçu franchit la tension de seuil Vs, puis on identifie l'instant du passage à zéro précédent (ou suivant) que l'on considère comme l'instant T1 de départ des oscillations caractéristiques Oc du signal acoustique reçu S2. Une telle méthode de mesure peut conduire à des erreurs dès lors que les oscillations caractéristiques Oc du signal acoustique reçu S2 peuvent être plus ou moins amplifiées en fonction de la nature du gaz. Ainsi, on a représenté sur la Figure 4 une courbe S21 qui correspond à la forme d'un signal acoustique reçu pour de l'azote (N2) et une courbe S22 qui correspond à la forme d'un signal acoustique reçu pour un mélange de gaz carbonique et de méthane (CO2/CH4). On constate que la courbe S21 franchit la tension de seuil Vs à un instant T4 qui déclenche l'identification du passage à zéro précédent à un instant T3 qui est considéré à juste titre comme la marque du début des oscillations caractéristiques Oc. En revanche, on peut constater que la courbe S22, qui est en phase avec la courbe S21, franchit la tension de seuil Vs à un instant T6 postérieur à l'instant T4 et décalé par rapport à celui-ci de la valeur d'une période TR du signal reçu. L'instant T5 identifié alors pour la courbe S22 comme le point de passage à zéro précédant immédiatement le franchissement de seuil à l'instant T6, est alors considéré comme la marque du début des oscillations caractéristiques c de la courbe S22. Or, comme on peut le voir sur la Figure 4, la courbe S22 présente un lobe négatif qui ne vient qu'affleurer la valeur de la tension de seuil Vs sans atteindre ou franchir ce seuil. Par suite des atténuations ou amplifications du signal reçu S2, qui varient en fonction de la nature du gaz, la méthode classique de mesure du temps du premier passage à zéro des oscillations caractéristiques du signal reçu peut engendrer une erreur de plus ou moins une période qui fait diminuer de façon significative la précision des mesures. L'invention vise à remédier aux inconvénients précités et à permettre de diminuer la sensibilité du procédé de mesure du temps de passage à zéro d'un signal acoustique reçu, vis-à-vis des variations de formes d'onde de ce signal ainsi que vis-à-vis d'éventuelles perturbations extérieures susceptibles d'être détectées par une simple méthode de surveillance de franchissement de seuil, et qui pourraient engendrer des mesures erronées. Ces buts sont atteints, conformément à l'invention, grâce à un procédé de mesure du temps de propagation d'un signal acoustique dans un écoulement de fluide entre un premier transducteur jouant le rôle d'émetteur et un second transducteur jouant le rôle de récepteur et situé à une distance déterminée du premier transducteur, le signal acoustique émis par le premier transducteur étant constitué par au moins une impulsion émise à une fréquence acoustique déterminée Fa et le signal acoustique reçu par le second transducteur comprenant une série d'oscillations caractéristiques dont l'amplitude est d'abord croissante sur plusieurs périodes, puis décroissante sur plusieurs périodes suivantes, I'enveloppe des oscillations caractéristiques présentant la forme d'un fuseau, le procédé consistant à échantillonner le signal acoustique reçu, à une fréquence d'échantillonnage Fe, à numériser le signal acoustique reçu échantillonné, et à rechercher, par analyse du signal acoustique reçu échantillonné et numérisé, le premier passage à zéro significatif des oscillations caractéristiques du signal acoustique reçu, caractérisé en ce que, pour rechercher le premier passage à zéro significatif des oscillations caractéristiques du signal acoustique reçu, on définit au préalable une première période caractéristique idéale déterminant le premier passage par zéro des oscillations caractéristiques du signal acoustique reçu, en caractérisant cette période idéale par un rapport d'amplitudes théorique A entre les amplitudes maximales Pi- et Pi+ des deux lobes de cette période, on détermine pour chaque période du signal acoustique reçu échantillonné et numérisé, les amplitudes maximales P- et P+ des deux lobes de la période examinée, on compare le rapport de ces amplitudes P- et P+ au rapport d'amplitudes théorique correspondant de la période idéale puis, si le résultat de la comparaison est supérieur à une valeur de seuil Gs, on considère la période examinée comme une période parasite correspondant à du bruit tandis que, si le résultat de la comparaison est inférieur à cette valeur de seuil Gs, on considère la période examinée comme une période caractéristique et on détermine alors le passage par zéro entre deux lobes de cette période caractéristique, lequel passage à zéro est considéré comme le premier passage à zéro significatif des oscillations caractéristiques du signal acoustique reçu. Le rapport d'amplitudes théorique A entre les amplitudes maximales Pi- et Pi+ des deux lobes de ladite période idéale est déterminé au préalable à partir d'une moyenne, pour plusieurs gaz différents et à différents débits, du rapport entre les amplitudes maximales P- et P+ de la période caractéristique observée à partir de l'enregistrement de signaux acoustiques reçus. Le rapport d'amplitude A entre les amplitudes maximales Pi- et Pi+ des deux lobes de la première période caractéristique idéale d'un signal acoustique reçu est quasi constant en fonction de la nature des gaz et est indépendant d'un facteur de gain. Le choix d'un critère de comparaison basé sur ce rapport A permet de diminuer très fortement la dépendance de la mesure à la nature du gaz, et donc d'améliorer la précision et la fiabilité des mesures. Selon un premier mode de réalisation de l'invention, la comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitudes théorique A correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la façon suivante: Selon un deuxième mode de réalisation de l'invention, la comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitudes théorique A correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la façon suivante où z représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude maximale (P-) du lobe négatif de la période caractéristique examinée. Selon ce deuxième mode de réalisation, on limite les risques de pouvoir détecter des formes d'ondes ressemblant au critère de ressemblance ou gabarit G, qui seraient noyées dans du bruit. Pour réduire encore la probabilité de détecter des périodes ressemblantes de faible amplitude, selon un troisième mode de réalisation de l'invention, la comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitudes théorique correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la façon suivante: où z représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude maximale (P-) du lobe négatif de la période caractéristique examinée, et z' représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude maximale (P+) du lobe positif de la période caractéristique examinée. Afin d'éviter des erreurs dans la détermination des amplitudes maximales P- et P+ des deux lobes de chaque période examinée, erreurs qui peuvent être dues à un déphasage existant entre le signal acoustique reçu et le signal d'échantillonnage, on peut utiliser un algorithme d'approximation du maximum (en valeur absolue) à partir des points échantillonnés au voisinage de chaque maximum en vue d'obtenir une meilleure précision sur les valeurs P- et P+. Toujours pour éviter ces mêmes erreurs, lorsque le signal acoustique reçu est répétable, on effectue plusieurs mesures successives respectivement sur plusieurs signaux acoustiques consécutifs, en déphasant pour chaque signal acoustique le signal d'échantillonnage par rapport au signal acoustique considéré, afin d'obtenir plusieurs points échantillonnés décalés au voisinage du maximum. Selon une variante de réalisation permettant également de réduire la probabilité de détecter des périodes ressemblantes de faible amplitude, on effectue, pour chaque période examinée, deux ou trois calculs du critère de ressemblance G avec des valeurs du rapport d'amplitudes théorique différentes voisines du rapport d'amplitude théorique prédéterminé et ne s'écartant pas de celle-ci de plus de 10%, et l'on considère la période examinée comme une période caractéristique si le critère de ressemblance calculé G est dans tous les cas inférieur à la valeur de seuil Gs. Pour déterminer la valeur de seuil Gs du critère de ressemblance G, ainsi que la valeur du terme de décalage z, on trace au préalable, pour différents gaz et différents débits, en faisant varier le déphasage entre la fréquence d'échantillonnage Fe et le signal acoustique reçu et en appliquant un facteur de gain variant entre 0,45 et 1,5, un réseau de premières courbes représentant les valeurs maximums du critère de ressemblance G pour les périodes caractéristiques observées et un réseau de deuxièmes courbes représentant les valeurs minimums du critère de ressemblance G pour les périodes parasites observées, en fonction de différentes valeurs possibles du terme de décalage z, et on choisit la valeur de seuil Gs et celle du terme de décalage z en fonction de la marge de sécurité possible dans la zone intermédiaire entre le réseau de premières courbes et le réseau de deuxièmes courbes. A titre d'exemple, la valeur de seuil Gs peut être comprise entre 0,7 et 1,7 et la valeur du terme de décalage z peut être comprise entre 0,21 et 0,25 V. L'invention est avantageusement applicable à un procédé de mesure de la vitesse d'écoulement d'un fluide entre deux transducteurs disposés en des points de mesure espacés dans la direction d'un écoulement du fluide selon lequel la valeur de la vitesse d'écoulement à mesurer est obtenue en combinant une mesure des temps de propagation respectifs de chacun des deux signaux acoustiques émis entre les deux points dans des sens opposés par les transducteurs avec une mesure des déphasages acoustiques respectivement induits dans chaque signal acoustique par la propagation de chacun desdits signaux acoustiques dans l'écoulement, la mesure des temps de propagation respectifs de chacun des deux signaux acoustiques émis entre les deux points de mesure étant effectuée conformément au procédé selon la présente invention. D'autres caractéristiques et avantages de l'invention ressortiront de la description suivante de modes particuliers de réalisation, faite en référence aux dessins annexés, sur lesquels: - la Figure 1 est une représentation, sous la forme d'un schémabloc, de l'ensemble d'un dispositif de comptage de gaz connu auquel est applicable la présente invention - la Figure 2 est un diagramme montrant le principe d'une méthode connue de détermination du temps de propagation d'un signal acoustique dans un milieu fluide; - la Figure 3 est un diagramme agrandi d'une partie de la Figure 2 montrant le principe d'une méthode connue de détermination du temps de propagation d'un signal acoustique par comparaison du signal reçu avec une valeur de seuil; - la Figure 3a est un diagramme analogue à celui de la Figure 3 montrant la façon dont on opère sur un signal échantillonné, selon une méthode connue de détermination du temps de propagation d'un signal acoustique par passage à zéro; - la Figure 4 est un ensemble de deux diagrammes analogues à celui de la Figure 3 montrant les risques d'erreur dans une méthode connue de détermination du temps de propagation d'un signal acoustique faisant appel à une comparaison des signaux reçus avec une valeur de seuil; - la Figure 5 est un diagramme montrant le principe de mise en oeuvre du procédé selon l'invention pour la détermination du temps de propagation d'un signal acoustique dans un milieu fluide; et - la Figure 6 est un diagramme comportant des familles de courbes permettant de choisir des valeurs de paramètres utiles pour la mise en oeuvre du procédé selon l'invention. - la Figure 7 est un organigramme montrant un exemple des différentes étapes du procédé selon l'invention. Différents types de circuits électroniques peuvent être utilisés pour assurer un comptage de fluide par une méthode de mesure à ultrasons de la vitesse d'écoulement du fluide. La Figure 1 montre le schéma de principe d'un exemple de tels circuits électroniques qui sont associés à deux transducteurs ultrasonores 1, 2 disposés à distance l'un de l'autre dans une conduite 3 dans laquelle s'écoule un fluide tel qu'un gaz. Les deux transducteurs 1, 2 sont reliés à un bloc de commutation 4 qui comprend deux commutateurs 5, 6 et permet l'utilisation de chaque transducteur alternativement comme émetteur et récepteur. Un module d'émission 14 et un module de réception 17 sont reliés respectivement aux commutateurs 5, 6 du bloc de commutation 4. Le module d'émission 14 comprend un amplificateur opérationnel 16 et un convertisseur numérique-analogique 15. Le module de réception 17 comprend au moins un amplificateur 18 et un convertisseur analogiquenumérique 19 qui numérise et échantillonne simultanément le signal reçu. Une source d'énergie électrique 7 et un module 8 de gestion de l'alimentation électrique sont reliés notamment aux modules d'émission 14 et de réception 17 ainsi qu'au bloc de commutation 4 et à un microcontrôleur 10. Le microcontrôleur 10 comprend notamment une horloge à quartz 9, une unité arithmétique et logique, des circuits de mémoire vive et de mémoire morte, et peut coopérer avec des circuits d'affichage 13, une mémoire morte réinscriptible 12 et une liaison série 11 du type RS 232. Les compteurs de gaz destinés à équiper chaque abonné d'un réseau de distribution doivent être à la fois précis, fiables et le meilleur marché possible. Ces contraintes imposent d'éviter l'utilisation de composants coûteux, de mettre en oeuvre de petites piles d'alimentation de longue durée mais de capacité moyenne pour des raisons de sécurité et d'effectuer des calculs selon un procédé de mesure assurant à la fois une précision et une fiabilité élevées tout en restant suffisamment simple pour être économe en énergie. Le procédé de mesure selon l'invention qui sera explicité plus loin est ainsi avantageusement appliqué à un dispositif de mesure de la vitesse d'écoulement d'un fluide gazeux, permettant d'assurer le comptage de la consommation de ce fluide, en mettant en oeuvre deux transducteurs ultrasonores disposés à distance l'un de l'autre à l'intérieur d'une conduite, dans le sens d'écoulement du fluide. A titre d'exemple, les transducteurs ultrasonores peuvent fonctionner à une fréquence acoustique Fa de l'ordre de 40 kHz et la fréquence d'échantillonnage Fe des signaux acoustiques reçus est avantageusement de 320 kHz, c'est-à-dire huit fois plus élevée que la fréquence Fa. La fréquence d'échantillonnage Fe est avantageusement égale à un multiple de la fréquence acoustique Fa. Toutefois, la fréquence Fe peut également prendre d'autres valeurs. Lorsque les valeurs (multiples de Fa ou non multiples) de Fe sont trop basses, il est préférable d'utiliser par exemple un algorithme d'approximation au voisinage du maximum (en valeur absolue) de chaque période du signal acoustique, voire d'effectuer plusieurs mesures successives en déphasant à chaque nouvelle mesure le signal d'échantillonnage par rapport au signal acoustique considéré, ceci afin d'obtenir une plus grande précision sur les valeurs des amplitudes maximales (en valeur absolue). Ces méthodes seront détaillées ultérieurement. Lorsque les valeurs de Fe sont suffisamment élevées l'emploi de telles méthodes n'est pas nécessaire. Le signal acoustique émis par chaque transducteur est constitué par exemple par une impulsion émise à la fréquence acoustique déterminée Fa. Cette impulsion est par exemple rectangulaire. Plusieurs impulsions pourraient également être émises. D'une manière générale, le débit de fluide dans une conduite peut s'écrire: = = (StJ2) (Tu -Td)/Tu.Td où S désigne la section moyenne de passage offerte à l'écoulement entre les deux transducteurs acoustiques, L désigne la distance séparant les transducteurs, Td est le temps de propagation du signal acoustique émis dans le sens de l'écoulement par le premier transducteur (amont) jusqu'à la réception de ce signal acoustique par le second transducteur (aval), Tu est le temps de propagation du signal acoustique émis dans le sens contraire de l'écoulement par le second transducteur, jusqu'à la réception de ce signal acoustique par le premier transducteur (amont). Si la mesure de ce débit fait en outre intervenir une mesure des déphasages acoustiques induits dans chaque signal acoustique par la propagation de chacun des signaux acoustiques émis par les transducteurs dans l'écoulement, le débit de fluide dans une conduite peut s'écrire: De façon plus particulière, les temps de propagation Tu et Td peuvent s'écrire respectivement: Tu=L/(c-v) Td=L/(c+v) où c et v représentent respectivement la vitesse de propagation du signal acoustique et la vitesse du gaz. Pour une composition de gaz donnée, comme la vitesse de propagation c dépend principalement de la température qui varie peu lors des mesures de débit, le terme Tu.Td lui aussi varie peu et ce terme peut être recalculé avec une fréquence relativement faible. Toutefois, dès lors qu'il s'agit de la mesure absolue des temps de propagation, la précision obtenue pour cette mesure est directement reportée sur la précision de mesure du débit. Ainsi, il est essentiel d'éviter d'avoir un décalage sur cette mesure. Par exemple, si la précision recherchée sur le débit est de 1%, la précision sur la mesure des temps de propagation doit être voisine de 0,5%, ce qui correspond à une erreur inférieure à 2 As sur les temps de propagation avec les gaz courants et les longueurs usuelles L entre deux transducteurs (exemple: L = 150 mm et le gaz étant du méthane). Selon un procédé de mesure des temps de propagation Tu et Td, connu sous le nom de méthode de passage à zéro, on mesure les temps de propagation (Tu pz) et (Td pz) des signaux acoustiques émis par chacun des transducteurs jusqu'à la détection par passage à zéro et on en déduit les temps de propagation réels Tu et Td, avec Tu = (Tu pz - To) et Td = (Td pz - To) où To désigne une constante à retrancher pour obtenir le temps de propagation dans le gaz, To ne dépendant que de la nature des transducteurs et des circuits électroniques. La Figure 3a illustre la méthode de mesure par passage à zéro effectuée sur chaque signal ultrasonore reçu selon un procédé connu, en liaison avec un dispositif tel que celui de la Figure 1. Le signal reçu est amplifié par l'amplificateur 18 présentant un gain donné, puis converti en numérique et échantillonné simultanément dans le convertisseur 19. Dans l'exemple considéré, I'échantillonnage est effectué à une fréquence Fe par exemple égale à 8 Fa, soit 320 kHz si Fa vaut 40 kHz. Les instructions permettant d'effectuer les mesures, contenues dans une mémoire morte du microcontrôleur 10, permettent de comparer les valeurs de l'amplitude du signal reçu aux points échantillonnés et qui sont mémorisées dans une mémoire vive, à une valeur prédéterminée correspondant à un seuil Vs, et l'on recherche parmi les valeurs d'amplitude des points échantillonnés le premier point suivant A, qui a une valeur d'amplitude supérieure au seuil. A partir de ce point, on recherche le prochain passage par zéro de la courbe représentative du signal ultrasonore. Pour ce faire, on détermine deux points échantillonnés consécutifs B et C qui encadrent le point où l'amplitude de la courbe passe par la valeur zéro et l'on procède par interpolation linéaire entre B et C pour déterminer le point D et donc mesurer le temps de passage à zéro Tpz. Selon une variante de réalisation, illustrée sur les Figures 3 et 4 discutées précédemment, il est également possible de rechercher le dernier passage à zéro de la courbe représentative du signal ultrasonore avant que ladite courbe ne passe par le seuil. La présente invention vise à améliorer cette méthode de détermination d'un temps de propagation de signal par mesure du temps de passage à zéro. Comme on l'a mentionné précédemment en référence aux Figures 3 et 4, la simple comparaison entre les valeurs d'amplitude de tension aux points échantillonnés d'un signal acoustique reçu et une valeur de tension de seuil Vs peut conduire à des erreurs dans la détermination du temps de passage à zéro Tpz et donc dans le temps de propagation du signal acoustique. Si sur la Figure 4 le temps de passage à zéro T3 pour la courbe S21 relative au signal correspondant à un gaz constitué par de l'azote est correct, en revanche le temps de passage à zéro T5 pour la courbe S22 relative au signal correspondant à un gaz constitué par un mélange de dioxyde de carbone et de méthane, est incorrect et excessif. Le procédé selon l'invention permet d'éviter ce genre d'erreur et, d'une manière générale, de réduire la sensibilité du procédé de mesure du temps de passage à zéro à des variations de formes d'onde ainsi qu'à des éventuelles perturbations extérieures qui, si elles étaient détectées, engendreraient des erreurs de mesure. Selon l'invention, pour rechercher le premier passage à zéro significatif des oscillations caractéristiques 0c du signal acoustique reçu S2, on définit au préalable une première période caractéristique idéale, distincte des oscillations parasites de faible amplitude constituant du bruit, qui détermine le premier passage à zéro des oscillations caractéristiques Oc du signal reçu S2. Cette période idéale est caractérisée par un rapport d'amplitude A entre les amplitudes des valeurs maximales Pi- et Pi+ des lobes négatif et positif de cette période idéale. Lors de la réception d'un signal acoustique S2 à analyser, qui est échantillonné et numérisé de façon classique, on détermine, pour chaque période du signal acoustique reçu échantillonné et numérisé, les amplitudes des valeurs maximales P- et P+ des lobes négatif et positif de la période examinée. Pour cela, il n'est pas nécessaire d'examiner l'ensemble des points échantillonnés de la période, et la détermination des valeurs maximales (en valeur absolue) P- et P+ des lobes négatif et positif de la période examinée peut être faite par une recherche au simple voisinage de ces extremums. Cependant, il faut noter qu'il existe un déphasage entre le signal acoustique reçu et le signal d'échantillonnage. Ainsi, en fonction du déphasage existant, le point échantillonné représentera de façon plus ou moins fidèle les valeurs maximales (en valeur absolue) P- et P+ des lobes négatif et positif. Pour remédier à ce problème, il est possible d'utiliser un algorithme d'approximation du maximum (en valeur absolue) à partir des points échantillonnés au voisinage du maximum afin d'obtenir une estimation de la valeur réelle de ce maximum quel que soit le déphasage existant entre le signal acoustique reçu et le signal d'échantillonnage. Par exemple, on peut utiliser un algorithme d'approximation polynômiale connu en soi. Le problème précédemment évoqué pourrait également être résolu en augmentant la fréquence d'échantillonnage. Toutefois, lorsque le procédé selon l'invention est mis en oeuvre dans un compteur de fluide qui fonctionne sur piles, il est essentiel de réduire le plus possible la consommation énergétique du procédé et cette exigence n'est pas compatible avec une augmentation de la fréquence temps sur plusieurs mesures successives) à effectuer plusieurs mesures successives respectivement sur plusieurs signaux acoustiques consécutifs émis dans le même sens de propagation en déphasant le signal d'échantillonnage par rapport au signal acoustique lorsque l'on effectue une nouvelle mesure, afin d'obtenir après ces différentes mesures plusieurs points échantillonnés décalés au voisinage des amplitudes maximales P- et P+ des lobes négatifs et positifs. Par exemple, la première mesure est effectuée en échantillonnant normalement le signal acoustique reçu, puis la seconde mesure est effectuée sur le second signal acoustique en déphasant le signal d'échantillonnage par rapport audit second signal acoustique d'une valeur déterminée 1/2Fe. II est ainsi possible d'effectuer plusieurs déphasages successifs (de valeur 1/nue pour n déphasages) sur le deuxième signal acoustique et les signaux acoustiques suivants pour obtenir davantage de points échantillonnés au voisinage des maxima et ainsi avoir une précision encore plus grande sur ces amplitudes maximales P- et P+. II est également avantageux de combiner l'utilisation d'un algorithme d'approximation du maximum (en valeur absolue) avec cette solution afin d'accroître davantage encore la précision. On procède ensuite à une comparaison du rapport de ces amplitudes des valeurs maximales P- et P+ au rapport A correspondant de la période idéale. Si le résultat de la comparaison est supérieur à une valeur de seuil Gs, on considère la période examinée comme une période parasite correspondant à du bruit, et on continue l'analyse sur la période suivante. Si le résultat de la comparaison est inférieur à la valeur de seuil G5, on considère la période examinée comme une période caractéristique. On peut alors déterminer de façon fiable le passage par zéro associé à cette période examinée. Ce passage par zéro considéré comme le premier passage à zéro significatif des oscillations caractéristiques du signal acoustique reçu peut être avantageusement le passage par zéro situé entre les deux lobes positif et négatif de la période caractéristique. La Figure 5 montre une première oscillation caractéristique Oca de période AT avec un lobe positif dont l'amplitude présente une valeur maximale P+ et un lobe négatif dont l'amplitude présente une valeur maximale en valeur absolue P-. Le passage à zéro associé à cette période AT d'une première oscillation caractéristique Oca définit un temps de passage à zéro correspondant à un temps T7. Le rapport d'amplitude entre les amplitudes des valeurs maximales P- et P+ des deux lobes d'une période examinée du signal acoustique reçu constitue un critère de comparaison indépendant de la nature du gaz. Ainsi, le vieillissement de l'électronique peut entraîner des variations de gain qui n'affecteront pas ce rapport d'amplitude entre les amplitudes maximales P- et P+. La reconnaissance de forme du signal acoustique reçu peut ainsi s'effectuer d'une façon simple, qui limite la consommation d'énergie, puisque la comparaison s'effectue entre les seules amplitudes maximales P- et P+ des lobes et un rapport d'amplitudes théorique A entre les amplitudes maximales Pi- et Pi+ des deux lobes d'une période caractéristique idéale. Le rapport d'amplitudes théorique A est déterminé au préalable à partir d'une moyenne, pour plusieurs gaz différents et à différents débits, du rapport entre les amplitudes maximales P- et P+ de la période caractéristique observée à partir de l'enregistrement de signaux acoustiques reçus. Selon un mode particulier de réalisation, la comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitude correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la façon suivante: Le fait d'effectuer des comparaisons en prenant en compte des rapports d'amplitudes maximales des lobes d'une période du signal acoustique reçu, le rapport d'amplitudes théorique de référence A n'est pas affecté par la nature du gaz ou le gain de la chaîne d'amplification du signal reçu. Ceci contribue à faciliter les mesures. Toutefois, il pourrait se trouver, parmi les périodes parasites du signal qui constituent du bruit, des formes d'ondes ressemblantes à celle définie par le critère de comparaison G. Pour rendre la détection des périodes caractéristiques moins sensible au risque de prendre en compte des formes d'ondes ressemblantes noyées dans du bruit et qui ne constitueraient pas une première période ca des oscillations caractéristiques Oc du signal reçu, le critère de ressemblance peut être rendu légèrement dépendant de l'amplitude P-. Ainsi, selon un autre mode de réalisation particulier, a comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitudes correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la facon suivante où z représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude de la valeur maximale (P-) du lobe négatif de la période caractéristique examinée. Dans ce cas, le dénominateur du critère G est fortement diminué quand l'amplitude P- est voisine du terme de décalage z. Si le terme de décalage z présente une valeur légèrement supérieure à l'amplitude maximale du bruit (par exemple de l'ordre de 0,20 à 0,25 Volt), ceci permet d'effectuer plus rapidement une réjection des périodes parasites du signal reçu. Pour déterminer la valeur de Gs à laquelle doit être comparé le critère de ressemblance G, et la valeur du terme de décalage z, on trace un réseau de premières courbes représentant les valeurs maximum du critère de ressemblance G pour les périodes caractéristiques observées et un réseau de deuxièmes courbes représentant les valeurs minimum du critère de ressemblance G pour les périodes parasites observées, en fonction de différentes valeurs possibles du terme de décalage z, et on choisit la valeur de seuil Gs et celle du terme de décalage z, et on choisit la valeur de seuil Gs et celle du terme de décalage z en fonction de la marge de sécurité possible dans la zone intermédiaire entre le réseau de premières courbes et le réseau de deuxièmes courbes. A titre d'exemple, les familles de courbes 111 et 121 de la Figure 6 ont été obtenues en faisant varier le déphasage entre la fréquence d'échantillonnage et le signal et en appliquant un facteur de gain variant entre 0,45 et 1,5, les gaz pris en compte étant le méthane et un mélange de 20% de gaz carbonique et de 80% d'azote, les mesures ayant été faites à débit nul et à 7m3/heure. Les réseaux de courbes 111 et 121 permettent de choisir les valeurs du terme de décalage z et de la valeur de seuil Gs en fonction de la marge de sécurité voulue sur ces valeurs. Si l'on considère la Figure 6, on voit que l'on peut choisir pour G5 par exemple une valeur comprise entre environ 0,7 et 1,7 V et, pour le terme de décalage z, par exemple une valeur comprise entre 0,21 et 0,25 V. Si l'on choisit par exemple Gs = 2 et z = 0,22 V; la marge de sécurité sur z est comprise entre 0,20 et 0,25 V tandis que la marge de sécurité pour Gs est comprise entre 1 et 4. En pratique, avec cet exemple, une période examinée sera considérée comme caractéristique si G D'autres variantes de réalisation permettent d'accroître encore la sécurité et la fiabilité du procédé. Ainsi, selon un mode de réalisation particulier, la comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitudes correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la façon suivante où z représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude de la valeur maximale (P-) du lobe négatif de la période caractéristique examinée, et z' représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude de la valeur maximale (P+) du lobe positif de la période caractéristique examinée. Ceci permet de réduire la sensibilité à des périodes ressemblantes de faible amplitude et, sans augmenter sensiblement la complexité des calculs, de pouvoir diminuer le gain de la chaîne d'amplification jusqu'à des valeurs aussi faibles que 0,35, ce qui économise de l'énergie. Selon une autre variante, et quelle que soit la formule choisie pour le critère de ressemblance G, on effectue, pour chaque période examinée, deux ou trois calculs du critère de ressemblance G avec des valeurs du rapport d'amplitudes A différentes voisines du rapport d'amplitude théorique prédéterminé et ne s'écartant pas de celle-ci de plus de 10%, en ce que l'on considère la période examinée comme une période caractéristique si le critère de ressemblance calculé G est dans tous les cas inférieur à la valeur de seuil G5. La Figure 7 est un organigramme montrant un exemple de processus de mise en oeuvre du procédé selon l'invention. Si l'on considère un échantillonnage à une fréquence Fe, par exemple de 320 kHz, qui représente huit fois la fréquence acoustique Fa des transducteurs, le signal échantillonné comportera 8 points par période du signal reçu. Dans la mesure où le signal acoustique reçu présente la forme d'oscillations périodiques, on sait où se trouvent les minimums et maximums successifs, ce qui évite d'examiner l'ensemble des points échantillonnés. On commence donc, après la phase d'initialisation 201, à chercher au cours d'une étape 202, le prochain minimum local (d'amplitude P-). Au cours de cette étape 202, si l'on applique un critère de ressemblance G faisant intervenir un terme de décalage z, on ne retient que le premier minimum local dont l'amplitude P- est supérieure en valeur absolue au terme de décalage z. Quant ce premier minimum local a été trouvé, on passe à l'étape 203 qui consiste à chercher au voisinage d'une demi-période T/2 précédente, la valeur de l'amplitude maximale P+. A l'étape suivante 204, on calcule la valeur du critère de ressemblance G avec les valeurs précédemment trouvés du minimum Pet du maximum P+ de la période examinée qui présente au moins un minimum dont la valeur absolue est supérieure au terme de décalage z. A l'étape 205, on effectue un test de comparaison entre la valeur du critère G calculée et la valeur de seuil Gs. Si G 2 Gs pour la période examinée, on considère qu'il s'agit d'une période parasite, et il y a retour à l'étape 202 de recherche d'un prochain minimum local. Si G Le même processus sera renouvelé quelques instants plus tard pour une nouvelle mesure de détermination du temps de propagation d'un signal acoustique ultrasonore. REVENDICATIONS 1. Procédé de mesure du temps de propagation d'un signal acoustique dans un écoulement de fluide entre un premier transducteur (1) jouant le rôle d'émetteur et un second transducteur (2) jouant le rôle de récepteur et situé à une distance déterminée du premier transducteur (1), le signal acoustique émis par le premier transducteur (1) étant constitué par au moins une impulsion S1 émise à une fréquence acoustique déterminée Fa et le signal acoustique reçu par le second transducteur (2) comprenant une série d'oscillations caractéristiques dont l'amplitude est d'abord croissante sur plusieurs périodes, puis décroissante sur plusieurs périodes suivantes, I'enveloppe des oscillations caractéristiques présentant la forme d'un fuseau, le procédé consistant à échantillonner le signal acoustique reçu, à une fréquence d'échantillonnage Fe, à numériser le signal acoustique reçu échantillonné, et à rechercher, par analyse du signal acoustique reçu échantillonné et numérisé, le premier passage à zéro significatif des oscillations caractéristiques du signal acoustique reçu, caractérisé en ce que, pour rechercher le premier passage à zéro significatif des oscillations caractéristiques du signal acoustique reçu, on définit au préalable une première période caractéristique idéale déterminant le premier passage par zéro des oscillations caractéristiques du signal acoustique reçu, en caractérisant cette période idéale par un rapport d'amplitude théorique entre les amplitudes maximales Pi- et Pi+ des deux lobes de cette période, on détermine pour chaque période du signal acoustique reçu échantillonné et numérisé, les amplitudes maximales P- et P+ des deux lobes de la période examinée, on compare le rapport de ces amplitudes P- et P+ au rapport d'amplitudes théorique correspondant de la période idéale puis, si le résultat de la comparaison est supérieur à une valeur de seuil Gs, on considère la période examinée comme une période parasite correspondant à du bruit tandis que, si le résultat de la comparaison est inférieur à cette valeur de seuil Gsl on considère la période examinée comme une période caractéristique et on détermine alors le passage par zéro entre deux lobes de cette période caractéristique, lequel passage à zéro est considéré comme le premier passage à zéro significatif des oscillations caractéristiques du signal acoustique reçu. 2. Procédé selon la revendication 1, caractérisé en ce que le rapport d'amplitudes théorique A entre les amplitudes maximales Pi- et Pi+ des deux lobes de ladite période idéale est déterminé au préalable à partir d'une moyenne, pour plusieurs gaz différents et à différents débits, du rapport entre les amplitudes maximales P- et P+ de la période caractéristique observée à partir de l'enregistrement de signaux acoustiques reçus. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que la comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitudes théorique A correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la façon suivante: 4. Procédé selon la revendication 1 ou 2, caractérisé en ce que la comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitudes théorique A correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la façon suivante où z représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude maximale P- du lobe négatif de la période caractéristique examinée. 5. Procédé selon la revendication 1 ou 2, caractérisé en ce que la comparaison entre le rapport des amplitudes P- et P+ de la période examinée et le rapport d'amplitudes théorique A correspondant de la période idéale s'effectue en calculant un critère de ressemblance G qui constitue le résultat de la comparaison comparé à la valeur de seuil Gs, et qui est défini de la façon suivante où z représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude maximale P- du lobe négatif de la période caractéristique examinée, et z' représente un terme de décalage destiné à rendre le critère de ressemblance G légèrement dépendant de l'amplitude maximale P+ du lobe positif de la période caractéristique examinée. 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que, pour déterminer avec plus de précision les amplitudes maximales P- et P+ des deux lobes de chaque période examinée, on utilise un algorithme d'approximation du maximum (en valeur absolue) à partir des points échantillonnés au voisinage de chaque maximum. 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que, pour déterminer avec plus de précision les amplitudes maximales P- et P+ des deux lobes de chaque période examinée, on effectue plusieurs mesures successives respectivement sur plusieurs signaux acoustiques consécutifs en déphasant pour chaque signal acoustique le signal d'échantillonnage par rapport au signal acoustique considéré, afin d'obtenir après ces différentes mesures plusieurs points échantillonnés décalés au voisinage des amplitudes maximales. 8. Procédé selon l'une quelconque des revendications 3 à 7, caractérisé en ce que l'on effectue, pour chaque période examinée, deux ou trois calculs du critère de ressemblance G avec des valeurs du rapport d'amplitudes théorique A différentes voisines du rapport d'amplitudes théorique prédéterminé et ne s'écartant pas de celle-ci de plus de 10%, et en ce que l'on considère la période examinée comme une période caractéristique si le critère de ressemblance calculé G est dans tous les cas inférieur à la valeur de seuil Gs. 9. Procédé selon la revendication 4, caractérisé en ce que, pour déterminer la valeur de seuil G5 du critère de ressemblance G, ainsi que la valeur du terme de décalage z, on trace au préalable, pour différents gaz et différents débits, en faisant varier le déphasage entre la fréquence d'échantillonnage Fe et le signal acoustique reçu et en appliquant un facteur de gain variant entre 0,45 et 1,5, un réseau de premières courbes (111) représentant les valeurs maximums du critère de ressemblance G pour les périodes caractéristiques observées et un réseau de deuxièmes courbes (121) représentant les valeurs minimums du critère de ressemblance G pour les périodes parasites observées, en fonction de différentes valeurs possibles du terme de décalage z, et on choisit la valeur de seuil Gs et celle du terme de décalage z en fonction de la marge de sécurité possible dans la zone intermédiaire entre le réseau de premières courbes (111) et le réseau de deuxièmes courbes (121). 10. Procédé selon la revendication 4 ou 9, caractérisé en ce que la valeur de seuil Gs est comprise entre 0,7 et 1,7 et la valeur du terme de décalage z est comprise entre 0,21 et 0,25 V. 11. Procédé de mesure de la vitesse d'écoulement d'un fluide entre deux transducteurs disposés en des points de mesure espacés dans la direction d'un écoulement du fluide, selon lequel la valeur de la vitesse d'écoulement à mesurer est obtenue en combinant une mesure des temps de propagation respectifs de chacun des deux signaux acoustiques émis entre les deux points dans des sens opposés par les transducteurs, avec une mesure des déphasages acoustiques respectivement induits dans chaque signal acoustique par la propagation de chacun desdits signaux acoustiques dans l'écoulement, caractérisé en ce que la mesure des temps de propagation respectifs de chacun des deux signaux acoustiques émis entre les deux points de mesure est effectuée selon le procédé conforme à l'une quelconque des revendications 1 à 10.