FR 2478321 A1 19810918 FR 8005836 A 19800314 1 2478321 La présente invention concerne un système de me- sure de vitesse d'un fluide dans les trois directions de l'espa- ce à partir de phénomènes turbulents et en particulier dans l'air et dans l'eau. Elle peut s'appliquer à tout système de télédé- tection acoustique et en particulier dans les buts d'études météorologiques et de surveillance des pollutions atmosphéri- ques et dans les buts d'études océanographiques. Le principe d'une telle mesure est connu il consiste à émettre un train d'ondes sonores, à une fréquence donnée; ce signal émis,est diffusé à distance par la turbulence thermique à petite échelle et une faible partie de l'onde dif- fusée est renvoyée en direction du récepteur. Cette turbulence est transportée par le fluide de telle sorte que le signal dif- fusé est affecté de la variation apparente de fréquence, dite effet Doppler, qui correspond à la vitesse du fluide dans la zone diffusante. L'émission est réalisée par une antenne direc- tive généralement du type parabolique, au foyer de laquelle est situé un générateur d'ondes sonores. La réception est réalisée par le même type d'antenne, le signal diffusé étant capté par un transducteur quelconque, tel qu'un microphone, situé au voi- sinage du foyer de l'antenne. Parmi les stations mettant en oeuvre ce princi- pe on distingue deux types qui l'un et l'autre sont améliorés par l'utilisation du système selon l'invention - le système de mesure du type bistatique o l'antenne récep- trice est distincte de l'antenne émettrice. - le système de mesure du type monostatique o une même antenne fonctionne successivement comme émetteur et comme récepteur. Les stations du type monostatique utilisent généralement trois antennes; elles sont dites alors du type "trimonostatique", o chaque antenne, généralement parabolique, est repérée par le site et l'azimuth de son axe d'émission. L'ensemble des trois antennes détermine alors un système de coordonnées précis dans lequel sont repérées les trois compo- santes du vecteur de vitesse du fluide, chaque composante étant mesurée selon l'axe des antennes. La valeur d'une composante est donnée après analyse du signal capté, par la mesure du "glisse- ment" de fréquence qui affecte la fréquence de l'onde émise si Fo désigne cel3M--ci, la fréquence de l'onde diffusée sera 2 247832 1 F0 +AFd o la "composante Doppler" àFd., ou glissement de fréquence, est donnée par:Fd V d = 2 x - F0 C o o C et V désignent respectivement la célérité du son dans le fluide et la composante de la vitesse à mesurer. Les problèmes essentiels posés par une telle me- sure consistent à identifier la composante Doppler rélle compte- tenu, d'une part, du faible niveau réfléchi et d'autre part, de la présence de bruit parasite, au voisinage de la fréquence d'émission, qui va "biaiser" la valeur brute du glissement de fréquence: la valeur mesurée AFm étant alors inférieure à la valeur de la composante Doppler réelle AF Les techniques d'exploitation d'un tel signal sont très variées. Parmi celles-ci notons la méthode de l'ana- lyse spectrale par transformée de Fourier rapide, la technique des filtres adjacents nombreux à bandes passantes fines et cou- vrant l'ensemble de la bande passante utile, la technique du passage à zéro, la technique du "phase lock"... Mais ces métho- des ne donnent pas satisfaction parce qu'elles sont généralement lourdes, imprécises et peu fiables. En particulier, la méthode d'analyse spectrale nécessite un échantillonnage fin d'o un temps de calcul important difficilement compatible avec une me- sure en temps réel. La méthode des filtres adjacents est peu fiable et d'autre part les filtres successifs sont difficiles à mettre en oeuvre; quant à la technique du "passage à zéro", bien que simple d'emploi, elle ne donnait jusqu'à présent que des résultats médiocres car fortement faussée par la présence de bruits de fond. La présente invention remédie à ces inconvénients en constituant un système de traitement du signal à partir d'une technique connue d'exploitation du signal combinée à un ensemble de tests ou de critères de validation guj permet de valider ou de fréquence d'invalider en temps réel la mesure'fidu signal et d'afficher immédiatement les trois composantes de la vitesse. Sur la base de cette technique on opère une me- sure sur un nombre prédéterminé de périodes. Cette mesure est validée ou invalidée à partir de critères qui reposent sur le franchissement de seuils dans des intervalles de temps limités. Ces critères sont effectués d'une part, à l'échelle de chaque période, sur les intervalles de temps séparant le franchissement 3 2478321 soit de deux seuils soit du même seuil deux fois de suite dans le même sens, d'autre part sur la différence de temps de fran- chissement de seuils entre le début et la fin du nombre de pé- riodes prédéterminé. - Le résultat brut obtenu au niveau de chaque mesure est toutefois fortement biaisé par la présence de bruits de fond. Etant donné que l'on connait les caractéristiques du signal utile (signal à bande étroite) et du bruit de fond (signal à bande large dont la densité spectrale est sensiblement linéaire au voisinage de la bande considérée) on peut calculer la loi statistique de la différence entre la valeur de la mesure brute et de la valeur réelle. Un procédé nouveau permettant l'estimation des termes intervenant dans la correction qui en résulte, est une autre caractéristique importante de l'invention. Les termes intervenant dans cette correction sont en particulier les va- leurs efficaces du signal et du bruit de fond dans des bandes de fréquences voisines de la fréquence d'émission. Plus préci- sément la correction permettant d'obtenir le glissement réel de fréquence consiste en l'application d'un facteur sur le glisse- ment brut mesuré, égal au rapport des puissances efficaces du bruit et du signal. Dans un mode de réalisation préféré le signal global capté par l'antenne est transformé par un transducteur en un signal électrique analogique qui après un conditionnement préalable, est dérivé dans deux circuits en parallèle comportant chacun un système de filtrage associé à des fonctions électro- niques classiques. Le premier système de filtrage appelé "filtre de travail" est centré sur la fréquence d'émission. Sa largeur de relative aF/F est de l'ordre de 10o; cette bande est dénom- mée bande de travail". La mesure élémentaire de fréquence s'ef- fectue à partir du signal prélevé en sortie de "filtre de tra- vail selon la technique classique du passage à zéro. Le signal est soumis à un comptage de périodes au moyen d'une base de temps précise dont le déclenchement s'opère après franchisse- ment du seuil bas. Celui-ci est déterminé au préalable relati- vement à l'écart-type du bruit de fond électronique de l'appa- reillage. Il participe en outre aux critères de validation selon l'invention. Chaque période du signal est validéeen fonction de 4 2478321 l'intervalle de temps maximal séparant le franchissement entre un le seuil bas etkseuil haut; ce dernier est déterminé relative- ment à l'écart-type du bruit de fond acoustique; sa valeur est avantageusement pilotée en temps réel par la mesure simultanée du niveau de bruit acoustique. Cet intervalle de temps de fran- chissement doit correspondre, à une certaine tolérance près, à une fraction de la période du signal d'émission. Deux autres critères possibles utilisent de façon analogue les intervalles de temps qui doivent correspon- dre aux franchissements successifs d'un même type de seuil, soit du seuil bas, soit du seuil haut. Comme le précédent ces intervalles peuvent être prédéterminés et affichés dans une mémoire. Enfin un quatrième critère est l'intervalle de temps que ne doit pas dépasser, en tenant compte d'une tolérance due à l'effet Doppler, un nombre déterminé de périodes du si- gnal rétrodiffusé par exemple 5 périodes. La figure 5 est une représentation des critères de validation appliqué à la mesure de lafréquence du signal un des tests de validation sur une période étant représenté ici par 16, le test sur un nombre prédéterminé de périodes étant représenté par 22. En parallèle à cette mesure de fréquence, le même signal, en sortie de ce "filtre de travail", est dérivé dans une branche comportant en série amplificateur redresseur et intégrateur et qui donne une grandeur analogique fonction de travail de la valeur efficace V F(S) du signal S dans la bande passante; et de celle VF (B) du bruit de fond (B) dans la même bande avec: VF (B+S) = V P(B±S); en appelant P la puissance efficace, mesurée par le carré de la valeur efficace V correspondante. F -------------------------------------- Avec l'hypothèse justifiée que le bruit de fond acoustique B et le signal utile S (porteur de la composante Doppler) sont pra- tiquement décorrellés, on a approximativement P (B+S) = P (B) + P(S) Le second système de filtrage équipé du "filtre d'estimation de bruit" donne une grandeur analogique propor- tionnelle à la valeur efficace du signal dans la bande passante du filtre dite "bande d'estimation du bruit" soit dans le cas général: VF (B'+S')=V P(B'+S') 2478321 Avec la même hypothèse que précédemment on a appro- ximativement: P (B'+S') = P(B') + P (S') avec P(S')= puissance efficace du signal reçu dans la bande "d'estimation de bruit". P(B') = puissance efficace du bruit de fond acoustique reçu dans la bande "d'estimation de bruit". Ces grandeurs sont en relation directe avec celles mesurées par"le filtre de travail" car on a constaté (fig. 2) une évolution sensiblement linéaire décroissante du bruit de fond aérien, relativement à la fréquence, dans une bande de fréquence suffisamment large au voisinage de la fréquence d'é- mission. En pratique cette évolution se traduit par VF (B') # XVF (B) d'o P (B') AP (B) Xétant voisin du rapport des largeurs des bandes passantes respectives du"filtre d'estimation de bruit" et du "filtre de travail". De même les valeurs efficaces VF (SI) et VF (S) mesurées dans les bandes d'estimation de bruit et de travail sont dans un rapport T dont les valeurs dépendent comme suit du choix de la bande d'estimation de bruit &*) dans une première version la bande passante du "filtre d'estimation de bruit" recouvre la bande passante"de travail" en étant plus large que celle-ci. Dans ce cas: T _ d (/I >O dans le cas général) ) dans une seconde version la bande passante "d'estimation de bruit" est sensiblement ad- jacente à celle "de travail". Dans ce cas: /t.Ocar P (S') = O t) dans unetroisième version la bande passante'W'estimation de bruit" est constituée par deux bandes de fréquence restreintes situées au voisinage et de part et d'autre de la bande passante "de travail". Dans ce cas on a: "iO car P(S') = O Ces informations numérisées par un système clas- sique de conversion analogique/digitale sont transférées, le moment voulu, dans un microcalculateur. Celui-ci efectue une évaluation du rapport P (B)/P (S) dans la bande de fréquence "de travail" par une résolution classique du système linéaire P (B) + P (S) = P (B+S) 6 2478321 A 2 p (B) +p p (S) = p (B'+S') Compte tenu des hypothèses précédentes et en par- ticulier des caractéristiques du signal diffusé, schématisées sur la figure 4, à savoir signal à bande étroite ff L 1/i' # Do situé dans une bande passante principale:2,aFtelle que 4,x Le calcul du glissement réel de fréquence du signal utile ou composante Doppler se fait alors en première approximation selon l'algorithme: AFd = tFm (1 + P(B)) la vitesse se calculant suivant l'expression classique préci- tée:Fd 2V F C o On notera que la "bande passante"de travail" 2 AiF est choisie de telle manière qu'elle permette de mesurer toutes les vitesses de vent situées dans une plage fixée AV soit telle que: àE > * 2 AV F0 C L'invention sera précisée à l'aide du complément de description qui va suivre, qui se réfère aux dessins annexés et se rapporte plus particulièrement à un exemple de mise en oeuvre de l'invention destiné à mesurer à la vitesse du vent. - La figure 1 représente un diagramme spatio temporel de pro- pagation d'une impulsion acoustique. - la figure 2 est une représentation d'une évolution type du bruit de fond acoustique et électronique en fonction de la fréquence sur laquelle est située la position relative des bandes passantes du "filtre de travail" et du "filtre d'estima- tion de bruit". - La figure 3 est un bloc diagramme de l'ensemble du système de télédétection acoustique. - Les figures 4 et 5 ont été citées plus haut. - La figure 6 est un schéma fonctionnel de la chaîne de trai- tement de signal. Le diagramme spatio-temporel de la figure 1 correspond à l'émission d'un pinceau d'ondes sonores 1 à une 7 2478321 fréquence déterminée très stable; celles-ci sont diffusées par la turbulence qui affecte les différentes couches de fluide at- teintes par le signal. La fréquence d'émission utilisée est choisie en fonction de la distance R d'exploration retenue, ainsi que du fluide mesuré. Elle peut être par exemple de 2000 Hz dans l'atmosphère et 10000 Hz dans l'eau. La résolution spatiale ou profondeur d'exploration, est déterminée par la durée 20de l'impulsion sonore d'amplitude A. Une faible partie de l'onde émise 1 est diffusée en direction de l'antenne 8 et est captée après un intervalle de temps correspondant au temps requis par les ondes sonores 1 pour parcourir la somme des distances telles que R s'éparant l'émission 1 et l'antenne 8 de la zone diffusante 2. Les signaux diffusés à différentes distances, et ici en particulier à différentes altitudes s'échelonnent dans le temps, au niveau de la réception, selon un ordre corres- pondant aux distances croissantes des zones turbulentes ren- contrées, à partir d'un instant: 2 R fixé par le microcalcu- lateur, sur un intervalle de temps----2 e donné, comme repré- senté sur le diagramme. Sur la figure 2 est représentée une évolution type de la densité spectrale du bruit de fond acoustique 4-a et électronique 46 en fonction de la fréquence. Sur cette même figure on a représenté la position de la "bande de travail" 18 centrée sensiblement sur la fréquence s'émission F et en pointillé les différentes positions possibles de la "bande d'estimation de bruit" 17. Ces deux bandes de fréquences 18 et 17 sont complémentaires et couvrent une bande de fréquence quasi conti- nue dans la zone d'évolution sensiblement linéaire 19 du bruit de fond acoustique. La figure 3 représente le bloc-diagramme d'un système de télédétection acoustique. La transmission du signal acoustique est réalisée par une chaîne d'émission comprenant à l'amont un générateur 6 de fréquence électronique stable du type oscillateur à cristal par exemple. Chaque train d'ondes émis 3 est initialisé par un signal issu du microcalculateur et sa durée &est déterminée par une logique d'émission in- tégrée dans le système de génération de fréquence 6. Le signal ainsi délivré à une fréquence donnée très stable, est amplifié 8 ?478321 dans un amplificateur de puissance 7 avant d'exciter une chambre de compression, située au foyer d'une antenne 8, qui génère les ondes sonores 1. A la réception, le signal acoustique diffusé 9 est transformé par un capteur, tel qu'un microphone, ou par la membrane de la même chambre de compression, en un signal électrique analogique qui est soumis à une préamplification dans un transformateur préamplificateur 10. Ce signal analogique global, porteur du bruit de fond et de la composante Doppler, est ensuite analysé dans une unité de traitement analogique 11 qui donne en continu des informations analogiques sur les carac- téristiques du signal 9 capté, telles que: fréquence et valeur efficace selon des modalités qui font partie de.l'invention et seront détaillées plus loin. Ces informations sont transfor- mées par un convertisseur analogique/nojrn tr 12 et stockées dans une mémoire "tampon" du microcalculateur 5 avant d'être transférées, le moment voulu, dans l'unité de traitement numé- rique de ce microcalculateur 5 qui effectuera, sur ces valeurs, un certain nombre d'opérations élémentaires programmées. Les différents ordres de commande du microcal- culateur 5, à l'ensemble du système de mesure sont symbolisés par des liaisons 14 et en particulier les signaux de chrono- métrage. La figure 6 est le schéma fonctionnel de l'en- semble de la chaîne de traitement du signal selon l'invention. Ce signal en sortie de l'antenne 8 est dérivé dans deux circuits en parallèle 15a et 15b qui comportent respectivement, disposés en série: un système de filtrage 18 et 17, un amplificateur 23a, 23b, un redresseur 21a et 21b et un intégrateur 30a et 30b. Dans le circuit 15a la bande passante de largeur L1du filtre 18 appelé "filtre de travail" est centrée sur la fréquence d'émission F0. Dans le circuit 15b le filtre 17 applelé "filtre d'estimation de bruit" est représenté dans les trois versions possibles, selon que sa bande passante 20 est o) centrée sur la fréquence d'émission F avec une largeur de bande L2 plus grande que L1 adjacente à la bande passante Ru "filtre de travail" 18 dans cette version, il existe une autre configuration possible, non figurée, qui serait symétrique, par rapport à la bande 9 2478321 passante du "filtre de travail" 18a de celle qui est représentée. ) constituée de deux bandes disjointes 20a et 20b situées symétriquement par rapport à la bande passante du "filtre de travail" 18a et adjacentes à celle-ci. Le signal électrique analogique en sortie de l'amplificateur 23a, compris dans la bande de fréquences: F0 + F du "filtre de travail" 18 est dérivé dans deux branches en parallèle 24a et 24c. Dans la branche de dérivation 24a, le signal électrique ainsi filtré passe successivement dans un redresseur 21a et un intégrateur 30a. En sortie des intégrateurs 30a et 30b on dispose de grandeurs proportionnelles aux valeurs efficaces du signal reçu, dans les bandes passantes respectives du "filtre de travail" et du "filtre d'estimation de bruit" 17. Ces valeurs analogiques sont transformées, dans le convertisseur analogiquewImItrlq 12, en valeurs numériques pour chacun des niveaux ou tranches d'altitude analysés, puis stockées dans les adresses correspondantes d'une mémoire tampon du microcalculateur 5 avant d'être transférées, le moment voulu dans l'unité de traitement numérique de ce microcalculateur 5. En parallèle à cette mesure, effectuée sur la branche de dérivation 24a, le signal issu du "filtre de travail" 18 est analysé période par période dans la branche de dérivation 24c. Le signal ainsi filtré dans la bande: F0 tF porteur de la composante Doppler est dérivé dans deux lignes 25d et 25e comportant chacune un comparateur de seuil ou de niveau 26d et 26e. Lorsque l'amplitude du signal dépasse la valeur affichée par le comparateur de seuil bas S.B. 26e celui-ci déclenche un compteur de périodes 28 par l'intermédiaire d'un compteur d'im- pulsionSintégré dans un système de chronométrage 27 à double entrée qui comprend une logique de comptage. La seconde entrée du système de chronométrage 27 est en relation avec le comparateur à seuil haut S.H. 26d qui renseigne celui-ci sur la poursuite o l'arrêt du comptage des impulsions selon la valeur de l'amplitude du signal par 26d rapport au seuil haut S.H.4prédéterminé. Le nombre de périodes validées au moyen des décrits critères+précédemment, ainsi que le nombre d'impulsions corres- pondantes sont transformés en informations binaires par le 2478321 convertisseur analogique numérique 12 et stockés dans la mémoire tampon avant d'être transférés dans l'unité de calcul de micro- calculateur 5.' On dispose ainsi de la totalité des informations donnant à la fois le nombre de périodes valides successives ainsi que l'intervalle de temps correspondant à leur établisse- ment, la valeur efficace du signal porteur de la composante Doppler et du signal issu du filtre d'estimation de bruit. les informations issues du "filtre d'estimation de bruit" 17 soit: P (B' + S') sont respectivement identifiées dans le microcalculateur 5 à: VP =T+ >q1J P(B) dans la configuration 17o Dû.\JP P(B) dans la configuration 17- avec _ L2 Là Le transfert de ces informations vers le micro- calculateur 5 est symbolisé par la liaison 29. Le microcalcu- lateur 5 peut alors effectuer les calculs sur ces valeurs selon les algorithmes décrits plus haut. îl 2478321 R E V E N D I C A T I O N S 1. Système de mesure à distance des vitesses re- latives au sein d'un fluide en mouvement, en particulier de l'eau et de l'air et comportant une émission de trains d'ondes acoustiques selon au moins une direction déterminée: - une réception des ondes diffusées en particulier par la turbulence à distance de l'émission. - une mesure de la fréquence des ondes reçues. - une détermination du glissement de fréquence, égal à la différence des fréquences émises et reçues. - une évaluation de la composante de vitesse du fluide par la formule de DOPPLER. - caractérisé par une mesure de fréquence globale dans une bande passante dite "de travail" sensiblement centrée sur la fréquence d'émission et par une estimation simultanée du rapport entre les puissances du bruit de fond et du signal utile captés ce qui permet de corriger cette mesure globale. 2. Système selon la revendication 1 o la mesure de fréquence globale est réalisée selon la méthode classique du passage à zéro, caractérisé par une validation période par période de ladite mesure de fréquence. 3. Système selon la revendication 2 caracérisé en ce que la validation période par période est conditionnée par un intervalle de temps qui sépare le franchissement suc- cessif d'un seuil bas et d'un seuil haut. 4. Système selon l'une des revendications 2 ou 3 caractérisé en ce que la validation période par période est conditionnée par un intervalle de temps qui sépare le franchis- sement successif d'un même seuil dans le même sens. 5.Système selon l'une des revendications 3 ou 4 caractérisé en ce que la valeur du seuil bas est de l'ordre de grandeur de l'écart-type du bruit de fond électronique. 6. Système selon l'une des revendications 3 à 5 caractérisé en ce que la valeur du seuil haut est de l'ordre de grandeur de l'écart-type du bruit de fond acoustique dans la bande "de travail". 7. Système selon l'une des revendications 1 à 6 caractérisé par une validation de la mesure de fréquence, conditionnée par un intervalle de temps séparant le franchisse- ment d'un même seuil, dans le même sens, sur un nombre global fixé de périodes. ?478321 8. Système selon l'une des revendications 1 à 7 caractérisé en ce que l'estimation du rapport entre les puissan- ces du bruit et du signal captésdans la bande "de travail" est faite par l'intermédiaire d'une mesure simultanée dans uoebande de fréquence complémentaire, dite"d'estimation du bruit". 9. Système selon la revendication 8 caractérisé en ce que la bande "d'estimation de bruit" recouvre la bande "de travail", en étant plus large que celle-ci. 10. Système selon la revendication 8 caractérisé en ce que la bande"d'estimation de bruit" est adjacente à la bande "de travail" ou à son proche voisinage. 11. Système selon la revendication 10 caractérisé en ce que la bande "d'estimation de bruit" se compose de deux bandes passantes adjacentes à la bande "de travail", ou à son voisinage de part et d'autre de celle-ci. 12. Système selon l'une des revendications 8 à 11 caractérisé en ce que la valeur efficace des ondes reçues dans - la bande passante "d'estimation de bruit" est identifiée au pro- duit de la valeur efficace du bruit de fond acoustique dans la bande passante "de travail" par un facteur voisin du rapport des largeurs des bandes passantes "d'estimation de bruit" et "de travail" 13. Système selon l'une des revendications 1 à 12 caractérisé en ce que la correction permettant d'obtenir le glissement réel de fréquence, à partir de la valeur du glisse- ment brut, consiste en l'application d'un facteur multiplicatif, sur le glissement brut mesuré,égal au rapport des puissances efficaces du bruit et du signal dans la "bande de travail". 14. Appareil pour la mise en oeuvre du système selon l'une des revendications 1 à 13 caractérisé par une chaîne de traitement du signal comportant les éléments suivants: - une antenne réceptrice captant les ondes diffusées, équipée d'un transducteur délivrant un signal électrique, associé à un préamplificateur - une unité de traitement de ce signal, comportant, dans chacun de deux circuits en parallèle, un système de filtrage associé en série à des éléments tels que amplificateur, redresseur et intégrateur, le premier circuit équipé d'un filtre passe- bande centré sur la fréquence d'émission, appelé "filtre de travail", le second circuit équipé d'au moins un filtre passe- 13 2478321 bande appelé "filtre d'estimation de bruit" - un estimateur des puissances transmises à travers chacun des circuits - un système de comptage des périodes du signal, sur une déri- vation du premier circuit, " de travail" - une mémoire tampon dans laquelle ces données sont stockées ces - une unité de traitement numérique de'données calculant la va- leur corrigée de la vitesse du fluide à distance de l'appareil. 15. Appareil selon la revendication 14 rattachée à l'une des revendications 2 à 13 caractérisé en ce que le système de comptage de périodes est réalisé par un système de chronométrage associé à: - un comparateur de niveau à seuil bas qui déclence le chrono- métrage - un comparateur de niveau à seuil haut qui renseigne le système de chronométrage - au moins une mémoire de validation conservant la valeur de l'intervalle de temps pendant lequel doit s'effectuer le chronomètrage entre deux comparaisons de niveau. 16. Appareil selon la revendication 14 carac- térisé en ce que la bande passante du "filtre d'estimation de bruit" couvre avec la bande passante du "filtre de travail" une bande de fréquence quasi continue à l'intérieur des limites de laquelle le rapport des largeurs de bande de ces deux fil- > 25 tresest voisin du rapport de la densité spectrale du bruit dans les bandes considérées; la bande passante du "filtre de travail" sensiblement centrée sur la fréquence d'émission F couvrantâ l'intérieur de ces limites une plage de fréquences 2 a F telle que A F > 2A V F C o tAV désigne la plage des vitesses de vent que l'on veut mesurer.