FR 2952722 A1 20110520 FR 0958152 A 20091118 MESURE DE CARACTERISTIQUES VELOCIMETRIQUES OU VIBROMETRIQUES UTILISANT UN DISPOSITIF DE TYPE LIDAR A DETECTION HETERODYNE La présente invention concerne la mesure de caractéristiques vélocimétriques ou vibrométriques en utilisant un dispositif de type LIDAR à détection hétérodyne, aussi appelé LIDAR cohérent. Elle concerne un procédé ainsi qu'un dispositif pour effectuer une telle mesure. Les dispositifs de type LIDAR, pour «Llght Detection And Ranging» en anglais, sont très utilisés pour de nombreuses applications de détection et de mesure. Parmi ces applications, on peut citer la mesure de vitesses anémométriques, l'étude d'écoulements fluides, l'étude de phénomènes vibratoires, la mesure de vitesses de véhicules, des applications du domaine 1 o du transport aérien, etc. Leur fonctionnement consiste à émettre une onde optique cohérente en direction d'un volume-cible où est situé un objet ou un phénomène à étudier, et à recueillir une partie de cette onde qui est diffusée à partir du volume-cible. Dans un dispositif LIDAR cohérent, la partie recueillie de l'onde diffusée est soumise à une détection hétérodyne. La mesure de vitesse 15 est déduite d'un décalage Doppler qui est mesuré dans un signal de la détection hétérodyne. Toutefois, l'onde diffusée qui est recueillie peut provenir non seulement de l'objet qui est situé dans le volume-cible, vers lequel le dispositif LIDAR est pointé, mais aussi d'objets qui sont extérieurs au volume-cible et qui ont aussi 20 une capacité significative à diffuser l'onde optique. Par exemple, pour une mesure anémométrique, l'onde est diffusée par des gouttelettes ou des poussières qui se trouvent dans le volume-cible et qui sont entraînées par le vent, mais un nuage qui est situé au-delà du volume-cible peut avoir une capacité de diffusion optique très supérieure aux gouttelettes et poussières du 25 volume-cible. Dans ce cas, la présence du nuage à l'arrière-plan, avec éventuellement une vitesse de déplacement du nuage, peut gêner ou empêcher la détection de la partie d'onde qui est diffusée par le contenu du volume-cible. Une situation comparable est celle de la mesure de vitesse d'un 2952722 -2- écoulement fluide lorsqu'une paroi en arrière-plan produit elle-même une diffusion ou une réflexion importante de l'onde qui est émise par le dispositif LIDAR. Les dispositifs LIDAR pour lesquels le système d'émission de l'onde 5 optique est séparé d'un système de réception de l'onde diffusée permettent de sélectionner le volume-cible par rapport à l'arrière-plan. Les deux systèmes d'émission et de réception sont chacun orientés vers le volume-cible, si bien que celui-ci est sélectionné par rapport à l'arrière-plan, à l'intersection des directions de pointage respectives du système d'émission et du système de réception. Toutefois, de tels dispositifs LIDAR, dits bistatiques, nécessitent d'être installés avec précision pour sélectionner le volume-cible voulu, avec une distance de séparation entre les systèmes d'émission et de réception qui soit suffisante. Ils sont donc mal adaptés pour être utilisés à bord d'un véhicule ou d'un aéronef, notamment. Le plus souvent, ils peuvent difficilement être déplacés à cause de leur structure en plusieurs parties séparées. Enfin, les dispositifs LIDAR bistatiques sont mal appropriés pour réaliser une cartographie de la vitesse mesurée, en variant le volume-cible à l'intérieur d'un volume d'analyse beaucoup plus vaste. Les dispositifs LIDAR monostatiques, pour lesquels les systèmes d'émission et de réception d'onde sont réunis, ne présentent pas ces inconvénients. Autrement dit, ils peuvent facilement être déplacés et utilisés à bord de véhicules et d'aéronefs. Par contre, les mesures de vitesse qui sont obtenues avec de tels dispositifs monostatiques peuvent être perturbées et empêchées par un élément diffusant qui est situé dans la direction de pointage, en arrière ou en avant du volume-cible, même si l'éloignement du volume-cible est pris en compte dans la détection hétérodyne. En outre, un dispositif LIDAR monostatique permet de réaliser facilement une cartographie de vitesse mesurée, en variant la direction de pointage du dispositif pour balayer un champ d'analyse. Pour sélectionner la partie utile du signal de détection hétérodyne qui contient l'information de mesure, par rapport à la partie du même signal de détection hétérodyne qui ne provient pas de l'objet étudié, il est connu de 2952722 -3- répéter l'acquisition du signal hétérodyne, d'analyser séparément les signaux de détection hétérodyne qui sont obtenus à chaque répétition, et de cumuler tous les résultats de ces analyses. Une telle accumulation consiste le plus souvent à additionner les résultats des analyses spectrales. Elle permet de 5 séparer l'information de mesure qui concerne l'objet étudié lorsque cette information est constante, par rapport à un bruit d'arrière-plan qui est variable ou incohérent. Néanmoins, un tel traitement ne permet pas de distinguer entre deux objets qui sont situés dans la direction de pointage du dispositif LIDAR monostatique, et qui ont chacun une vitesse de déplacement sensiblement 10 constante. Enfin, il est aussi connu d'effectuer une détection synchrone de la partie d'onde qui est diffusée, pour la distinguer de contributions au signal détecté qui sont indépendantes de l'onde optique émise. Dans ces conditions, un premier but de l'invention est de fournir des 15 mesures de vitesse en utilisant un dispositif de type LIDAR monostatique à détection hétérodyne, qui ne présentent pas les inconvénients précédents. Un second but de l'invention est que de telles mesures puissent être réalisées en utilisant des dispositifs LIDAR existant avant la présente invention, en modifiant peu et facilement ces dispositifs et leur utilisation. 20 En particulier, l'invention a pour but de mesurer des vitesses avec une sélectivité du volume-cible par rapport à l'avant-plan et à l'arrière-plan dans la direction de pointage, qui soit améliorée. Notamment, une telle sélectivité améliorée est recherchée pour réaliser des cartographies de mesures. Enfin, un autre but de l'invention est de permettre des mesures avec un 25 dispositif LIDAR qui soit peu encombrant, et dont la masse n'est pas augmentée significativement par rapport aux dispositifs existants. Pour atteindre ces objectifs et d'autres, l'invention propose un procédé de mesure d'une caractéristique vélocimétrique ou vibrométrique qui utilise un dispositif de type LIDAR à détection hétérodyne, dans lequel un signal 30 d'émission d'onde optique est produit en direction d'un volume-cible à partir d'une tête optique du dispositif, et un signal rétrodiffusé est recueilli à travers cette même tête optique, puis est détecté par détection hétérodyne de façon à 2952722 -4- produire un signal de détection hétérodyne. Autrement dit, le dispositif LIDAR qui est utilisé est du type LIDAR cohérent monostatique. Le procédé comprend en outre les étapes suivantes pour un cycle de mesure : 5 /1/ produire une modulation d'une caractéristique de phase de l'onde optique dans le signal d'émission, telle qu'une modulation de fréquence ou une modulation de phase ; /2/ démoduler le signal de détection hétérodyne par rapport à la modulation de la caractéristique de phase, en compensant un retard de 1 o propagation du signal d'émission et du signal rétrodiffusé entre la tête optique du dispositif et le volume-cible ; /3/ combiner une analyse spectrale du signal de détection hétérodyne démodulé avec une accumulation pour des fenêtres temporelles successives de découpage, de façon à isoler une contribution au signal 15 de détection hétérodyne qui provient du volume-cible ; et /4/ obtenir un résultat de la mesure de caractéristique vélocimétrique ou vibrométrique, à partir d'une analyse d'effet Doppler de la contribution isolée à l'étape /3/, De plus, dans un procédé qui est conforme à l'invention : 20 - la modulation est obtenue en décalant la caractéristique de phase de l'onde optique pendant des créneaux temporels successifs de modulation («successive modulation time-slots» en anglais), conformément à un pas d'incrément fixe qui est multiplié par des facteurs multiplicatifs q respectivement affectés aux créneaux 25 temporels de modulation ; - les facteurs multiplicatifs q sont égaux à ap[n], où : n est un nombre entier qui est supérieur ou égal à quatre, et qui est constant pendant le cycle de mesure, a est un autre nombre entier, qui est aussi constant pendant le cycle 30 de mesure, 2952722 -5- p est un facteur de puissance entier qui est positif ou nul, strictement inférieur à n-1, et avec des valeurs variables qui sont affectées respectivement aux créneaux temporels de modulation, ap[n] désigne la réduction de ap modulo n, et 5 a est choisi tel que a[n] soit différent de zéro et de l'unité. Lorsque que la caractéristique de phase qui est décalée selon l'invention est la fréquence f de l'onde optique du signal d'émission, les décalages de celle-ci pendant les créneaux temporels successifs de modulation sont Af = Afo x q, Afo étant l'incrément fixe de fréquence. 1 o Alternativement, lorsque l'invention est appliquée à la phase initiale du signal d'émission, Lcp = ~cpo x q sont les décalages de la phase initiale de l'onde optique du signal d'émission, Acpo étant l'incrément fixe de phase initiale. Un procédé selon l'invention comporte donc les avantages qui sont procurés par un dispositif LIDAR monostatique à détection hétérodyne. 15 II comporte aussi les avantages qui sont procurés par l'accumulation de résultats d'analyse ou de détection, notamment pour séparer l'information qui est utile pour la mesure par rapport à un bruit variable et incohérent qui est aussi détecté pendant le cycle de mesure. En outre, l'application à la caractéristique de phase de l'onde optique, 20 de décalages qui varient d'un créneau temporel au suivant, améliore la résolution spatiale le long de la direction de pointage du dispositif monostatique. Typiquement, l'invention procure une résolution spatiale de mesure, le long de la direction de pointage, qui est de l'ordre de C x Tmoa / 2, où C est la vitesse de propagation de l'onde optique, et Tmoa est la durée 25 individuelle des créneaux temporels de modulation correspondant aux facteurs multiplicatifs q successifs. L'invention est encore améliorée dans sa facilité de mise en oeuvre et/ou dans sa capacité à améliorer la résolution spatiale qui est obtenue, en utilisant l'un et/ou l'autre des perfectionnements suivants : 30 - les facteurs multiplicatifs q des créneaux temporels successifs de modulation peuvent former une séquence qui est répétée -6-périodiquement pendant le cycle de mesure ; les facteurs multiplicatifs q successifs peuvent comprendre au moins trois valeurs différentes, et de préférence au moins cent valeurs différentes ; - le nombre n peut être un nombre premier («prime number» en anglais) ; et - le nombre entier a peut être une racine primitive («primitive root» en anglais) du groupe-quotient multiplicatif (Z/nZ)* («quotient-group (Z/nZ)* under multiplication» en anglais), où Z est l'anneau des 1 o nombres entiers relatifs, * indique que la valeur zéro est exclue, les facteurs multiplicatifs q étant alors n-1 valeurs distinctes pour p variant de 0 à n-2. Des mises en oeuvre de l'invention qui sont particulièrement préférées sont obtenues lorsque n est un nombre premier de l'ordre de quelques 15 centaines, et a est l'une des racines primitives du groupe-quotient multiplicatif (Z/nZ)*. Dans des premiers modes de mise en oeuvre de l'invention, la caractéristique de phase de l'onde optique qui est décalée peut être une fréquence de cette onde optique, de façon à réaliser une modulation de 20 fréquence. Dans ce cas, la durée individuelle des créneaux temporels de modulation est inférieure ou égale à la durée individuelle des fenêtres temporelles de découpage. L'étape /3/ comprend alors les sous-étapes suivantes : /3-1/ analyser spectralement le signal de détection hétérodyne démodulé 25 à l'intérieur d'une fenêtre temporelle de découpage, de façon à obtenir un résultat d'analyse spectrale pour cette fenêtre temporelle ; et /3-2/ répéter la sous-étape /3-1/ pour plusieurs fenêtres temporelles successives de découpage, et accumuler les résultats d'analyse spectrale qui ont été obtenus respectivement pour ces fenêtres 30 temporelles, la contribution au signal de détection hétérodyne qui provient du volume-cible étant prépondérante dans cette accumulation. 2952722 -7- Dans des seconds modes de mise en oeuvre de l'invention, la caractéristique de phase de l'onde optique qui est décalée peut être une phase initiale de cette onde optique, de façon à réaliser une modulation de phase. Selon une première possibilité pour une telle modulation de la phase, 5 la durée individuelle des créneaux temporels de modulation peut être inférieure ou égale à la durée individuelle des fenêtres temporelles de découpage. Dans ce cas, l'étape /3/ peut comprendre les sous-étapes suivantes : /3-1/ analyser spectralement le signal de détection hétérodyne démodulé à l'intérieur d'une fenêtre temporelle de découpage, de façon à obtenir 1 o un résultat d'analyse spectrale pour cette fenêtre temporelle ; et /3-2/ répéter la sous-étape /3-1/ pour plusieurs fenêtres temporelles successives de découpage, et accumuler les résultats d'analyse spectrale qui ont été obtenus respectivement pour ces fenêtres temporelles, la contribution au signal de détection hétérodyne qui 15 provient du volume-cible étant prépondérante dans cette accumulation. Selon une seconde possibilité pour la modulation de phase, la durée individuelle des créneaux temporels de modulation peut être supérieure ou égale à la durée individuelle des fenêtres temporelles de découpage. Dans ce cas, l'étape /3/ peut comprendre les sous-étapes suivantes : 20 /3-1/ accumuler des segments du signal de détection hétérodyne démodulé qui correspondent à des fenêtres temporelles successives de découpage, pendant les créneaux temporels successifs de modulation de sorte que la contribution au signal de détection hétérodyne qui provient du volume-cible soit accumulée de façon 25 constructive ; et /3-2/ analyser spectralement un résultat de l'accumulation. Enfin et d'une façon générale, l'invention est compatible avec l'utilisation d'un dispositif LIDAR qui est adapté pour produire un signal d'émission continu ou par impulsions. 30 L'invention propose aussi un dispositif LIDAR à détection hétérodyne, qui est du type monostatique et adapté pour effectuer des mesures 2952722 -8- vélocimétriques ou vibrométriques, et qui comprend en plus par rapport à des dispositifs existants : - un modulateur de phase, qui est disposé pour moduler au moins un signal de source de l'émission d'onde optique ; 5 - une unité de commande, qui est reliée à une entrée de commande du modulateur de phase, et qui est adaptée pour commander un fonctionnement de ce modulateur pour mettre en oeuvre un procédé tel que décrit précédemment ; et - des moyens de compensation d'un retard de propagation du signal 1 o d'émission et du signal rétrodiffusé entre la tête optique et le volume-cible. Un tel dispositif LIDAR selon l'invention peut être facilement obtenu en modifiant ou en programmant de façon appropriée un dispositif connu antérieurement. 15 D'autres particularités et avantages de la présente invention apparaîtront dans la description ci-après d'exemples de mise en oeuvre non limitatifs, en référence aux dessins annexés, dans lesquels : - les figures 1a et 2b montrent deux architectures de dispositifs LIDAR avec lesquelles l'invention peut être mise en oeuvre ; 20 - les figures 2a et 2b sont deux diagrammes qui montrent un résultat d'analyse de signal rétrodiffusé pour un premier mode de mise en oeuvre de l'invention ; et - les figures 3a et 3b correspondent respectivement aux figures 2a et 2b pour un second mode de mise en oeuvre de l'invention. 25 Dans les figures 1 a et 1 b, des éléments qui sont désignés par des références identiques sont eux-mêmes identiques ou ont des fonctions qui sont identiques. En outre, les références suivantes désignent les composants qui sont indiqués ci-après, et qui sont utilisés de la même façon que dans un dispositif LIDAR à détection hétérodyne, de type monostatique, tel que connu 30 avant la présente invention : -9- 100 : le dispositif LIDAR à détection hétérodyne dans son ensemble ; 1 un oscillateur laser, qui est adapté pour produire une onde optique OL; 2 : un séparateur optique, qui est disposé pour diviser l'onde optique OL en un signal de source d'émission SSE et un signal de référence SRef ; 3 un amplificateur optique, qui est adapté pour produire un signal d'émission SE à partir du signal de source d'émission SSE ; 11 une unité de mélange-détection («mixing and detection unit») ; 12 : une tête optique, qui est adaptée pour transmettre le signal d'émission SE en direction d'un volume-cible et pour recevoir un signal rétrodiffusé R ; et 20 : des moyens d'analyse. En outre, D désigne une direction de pointage du dispositif 100, qui est fixée par l'orientation de la tête optique 12. V est le volume-cible dans lequel se trouve l'objet dont une caractéristique vélocimétrique ou vibrométrique est recherchée. Par conséquent, la direction de pointage D est orientée vers le volume-cible V. Une source parasite de rétrodiffusion P peut être située à une distance quelconque de la tête optique 12 le long de la direction de pointage D, au delà du volume-cible V ou entre le volume-cible V et cette tête optique 12. Dans ces conditions, le signal rétrodiffusé R qui est recueilli par la tête optique 12 comprend une première contribution RV qui est produite par l'objet situé dans le volume-cible V, et au moins une seconde contribution RP qui est produite par la source parasite de rétrodiffusion P. Dans la suite, la contribution RV est appelée contribution utile, et la contribution RP est appelée contribution parasite. L'oscillateur laser 1, le séparateur optique 2 et l'amplificateur 3 peuvent appartenir à une unité laser d'émission 10, qui produit le signal d'émission SE. De façon connue et particulièrement avantageuse, chacun de ces composants peut être réalisé à partir d'au moins une fibre optique qui conduit l'onde optique 2952722 -10- OL ou les signaux SSE, SRef ou SE. De façon courante, l'onde optique OL et les signaux SSE, SRef et SE appartiennent à la bande de rayonnement électromagnétique infrarouge, pour laquelle la longueur d'onde est comprise entre 1,535 lm et 1,565 pm (micromètre). 5 L'unité de mélange-détection 11 est adaptée pour transmettre le signal d'émission SE à la tête optique 12, et pour produire un signal de détection hétérodyne SRF à partir du signal rétrodiffusé R qui est reçu par la tête optique 12 et du signal de référence Sref• Les moyens d'analyse 20 sont adaptés pour analyser le signal de 1 o détection hétérodyne SRF. D'une façon qui est couramment utilisée, ces moyens d'analyse 20 découpent le signal de détection hétérodyne SRF en segments de signal selon des fenêtres temporelles successives de découpage. Ils réalisent alors une combinaison d'opérations d'analyse spectrale et d'accumulation, ainsi qu'une analyse d'effet Doppler. 15 Le dispositif 100 peut comprendre aussi d'autres composants, d'une façon qui est usuelle pour des dispositifs LIDAR à détection hétérodyne. De tels composants ne sont pas repris dans la présente description, dans la mesure où ils n'ont pas de liaison directe avec l'objet de l'invention. En plus des composants précédents, un dispositif 100 selon l'invention 20 comprend : - un modulateur de phase 4, pour moduler au moins le signal de source d'émission SSE ; - une unité de commande 40, pour commander un fonctionnement du modulateur de phase 4 ; et 25 - des moyens de compensation d'un retard de propagation du signal d'émission SE et du signal rétrodiffusé R entre la tête optique 12 et le volume-cible V. Le modulateur de phase 4 peut comprendre une cellule de Pockels. Une telle cellule de Pockels produit un champ électrique variable à l'intérieur 30 d'un matériau actif qui est adapté pour modifier la phase d'une onde optique en fonction du champ électrique. De tels modulateurs de phase sont aussi bien -11 - connus de l'Homme du métier. Dans le dispositif 100 de la figure 1 a, le modulateur de phase 4 est disposé pour recevoir en entrée le signal de source d'émission SSE qui est produit par le séparateur optique 2, et pour transmettre en sortie le signal de source d'émission modulé à l'amplificateur optique 3. Autrement dit, le modulateur 4 est situé entre le séparateur 2 et l'amplificateur 3. Dans ce premier cas, seul le signal de source d'émission SSE comporte la modulation qui est produite par le modulateur 4. Le signal de détection hétérodyne SRF qui est produit par l'unité de mélange-détection 11 est alors démodulé avant d'être 1 o transmis aux moyens d'analyse 20. Pour cela, le dispositif comprend en outre des moyens de démodulation 5 qui sont agencés pour démoduler le signal de détection hétérodyne SRF conformément au fonctionnement du modulateur de phase 4. Les moyens de compensation du retard de propagation du signal d'émission SE et de la contribution utile RV sont agencés pour compenser ce 15 retard lors de la démodulation. Dans une réalisation avantageuse, les moyens de démodulation 5 peuvent être numériques et incorporer les moyens de compensation du retard. Dans le dispositif alternatif 100 de la figure 1 b, le modulateur de phase 4 est disposé pour recevoir en entrée l'onde optique OL qui est produite par 20 l'oscillateur laser 1, et pour transmettre en sortie cette onde optique modulée au séparateur optique 2. Autrement dit, le modulateur 4 est situé entre le l'oscillateur laser 1 et le séparateur 2. Dans ce second cas, le signal de source d'émission SSE et le signal de référence SRef sont modulés d'une façon identique. Les moyens de compensation du retard de propagation du signal 25 d'émission SE et de la contribution utile RV sont alors agencés pour retarder le signal de référence SRef qui est transmis à l'unité de mélange-détection 11. Par exemple, les moyens de compensation du retard peuvent comprendre une ligne à retard 6 qui est interposée sur la ligne de transmission du signal de référence SRef, entre le séparateur 2 et l'unité de mélange-détection 11. 30 Les deux dispositifs 100 des figures 1 a et 1 b sont équivalents vis-à-vis de l'invention, si bien que toutes les mises en oeuvre de l'invention qui sont décrites dans la suite peuvent utiliser indifféremment l'un ou l'autre de ces 2952722 -12- dispositifs. A titre d'illustration, le dispositif 100 de la figure 1 b est adopté dans la suite de cette description. Par ailleurs, l'Homme du métier sait aussi comment commander le fonctionnement du modulateur de phase 4 pour créer une modulation de 5 fréquence ou une modulation de phase de l'onde optique OL. Il est seulement rappelé que la phase d'une onde optique est de la forme 2ruf•t + cp, où f et cp sont respectivement la fréquence et la phase initiale de l'onde, et t le temps. La modulation de fréquence consiste à faire varier dans le temps la fréquence f, et la modulation de phase consiste à faire varier dans le temps la phase initiale (p. 1 o On décrit maintenant la modulation qui est introduite selon l'invention. Dans un premier temps, l'expression de caractéristique de phase désigne aussi bien la fréquence f que la phase initiale (p. La durée de détection du signal rétrodiffusé R, pour un même cycle de mesure, est divisée en créneaux temporels de modulation successifs. Cette 15 division en créneaux temporels est prévue pour la modulation du signal de source d'émission SSE. Elle est indépendante a priori du découpage de la même durée en fenêtres temporelles successives pour l'analyse du signal de détection hétérodyne SRF. La modulation est produite en décalant la caractéristique de phase 20 d'une façon qui est constante à l'intérieur de chaque créneau temporel, et qui varie d'un créneau temporel à un autre. Selon l'invention, ce décalage est multiple d'un incrément fixe, conformément à la formule suivante : A(Cphase) = A(Cphase)o x q (1) OÙ Cphase est la caractéristique de phase qui est modulée, A(Cphase) est le 25 décalage de cette caractéristique Cphase à l'intérieur de chaque créneau temporel de modulation, A(Cphase)o est l'incrément fixe et q le facteur multiplicatif qui est affecté à ce créneau temporel. En outre, selon l'invention, le facteur multiplicatif q est donné par la formule suivante : 30 q = ap[n] (2) où : n est un nombre entier qui est supérieur ou égal à quatre, 2952722 -13- a est un autre nombre entier tel que a[n] soit différent de zéro et de l'unité, et p est un facteur de puissance entier positif ou nul, strictement inférieur à n-1. 5 Les nombres n et a sont constants pendant chaque cycle de mesure. Ils peuvent être fixés définitivement par une programmation initiale de l'unité de commande 40. Le facteur de puissance p prend des valeurs qui varient pour des créneaux temporels différents. De préférence, il peut prendre toutes les valeurs 1 o entières de 0 à n-2, dans un ordre quelconque pour une série de créneaux temporels successifs. Toutefois, l'ordre croissant de ces valeurs de p : 0, puis 1, puis 2,... , jusqu'à n-2 est préféré. Eventuellement, seulement un nombre restreint de ces valeurs peuvent être utilisées pour le facteur de puissance p. Le facteur multiplicatif q et le décalage A(Cphase) varient donc en conséquence. 15 ap[n] désigne la réduction de ap modulo n. Autrement dit, le facteur multiplicatif q du décalage A(Cphase) de la caractéristique de phase est le reste d'une division euclidienne de ap par n. La valeur du facteur multiplicatif q pour chaque créneau temporel de la modulation est donc elle-même un nombre entier positif non nul. 20 D'une façon générale, le nombre entier a peut être négatif, mais il peut être sélectionné entre l'unité et n, ces deux valeurs étant exclues. Mathématiquement, les valeurs successives du facteur multiplicatif q sont les éléments d'un sous-groupe du groupe-quotient (Z/nZ)* doté de l'opération de multiplication. 25 De façon générale, un décalage A(Cphase) de la caractéristique de phase CphaSe qui a été introduit dans le signal d'émission SE, apparaît dans la contribution parasite RP avec un retard temporel, par rapport au même décalage A(Cphase) dans la contribution utile RV. Ce retard est dû à l'écart Ad entre les éloignements respectifs du volume-cible V et de la source parasite P, 30 selon la direction D. A l'inverse, un même décalage A(CphaSe) apparaît d'abord dans la contribution parasite RP si la source parasite P est située entre le dispositif 100 et le volume-cible V. Lorsque l'écart d'éloignement Ad est 2952722 -14- supérieur en valeur absolue à c x Tmod/2, où C est la vitesse de propagation du signal d'émission SE et Trima est la durée individuelle des créneaux temporels de modulation, la contribution parasite RP et le signal de référence SRef présentent entre eux, lors de la détection hétérodyne, des différences 5 successives de la caractéristique de phase Cphase qui sont variables. Il en résulte un étalement des valeurs de la caractéristique de phase Cphase qui sont relatives à la contribution parasite RP, dans l'accumulation qui est réalisée par les moyens d'analyse 20. A l'inverse, le retard de propagation de la contribution utile RV est compensé lors de la démodulation du signal de détection 1 o hétérodyne SRF (figure 1 a), ou par la ligne à retard 6 (figure 1 b). La contribution utile RV ne fait donc pas apparaître de décalage variable de la caractéristique de phase Cphase dans le signal de détection hétérodyne SRF qui est analysé, pendant tout le cycle de mesure. De cette façon, la contribution utile RV devient prépondérante dans l'accumulation qui est réalisée par les moyens 15 d'analyse 20 après une durée d'accumulation suffisante. Ainsi, la contribution utile RV peut être isolée dans le signal rétrodiffusé R, pour en déduire la vitesse du ou des objet(s) qui est (sont) situé(s) dans le volume-cible V. Cette ou ces vitesse(s) est (sont) déduite(s) en effectuant une analyse d'effet Doppler qui est bien connue de l'Homme du métier. 20 De préférence, n est un nombre premier suffisamment grand et a est une racine primitive, aussi désignée par générateur, du groupe-quotient (Z/nZ)* doté de la multiplication, dont la valeur nulle est exclue. Dans ce cas, le facteur multiplicatif q peut prendre successivement n-1 valeurs qui sont toutes différentes : a°[n]=l , al[n]=a, a2[n], ..., an-2[n]. Ces valeurs de q sont alors 25 suffisamment nombreuses pour produire un étalement efficace pour la contribution parasite RP. Cet étalement peut permettre de réduire le nombre de fenêtres temporelles successives pendant lesquelles la détection et l'analyse du signal rétrodiffusé R sont répétées, pour séparer la contribution utile RV de la contribution parasite RP. Cette séparation est efficace même si la 30 contribution parasite RP est plus beaucoup intense que la contribution utile RV. Par exemple, n peut être égal à 173 et a peut être égal à 11 ou 17, ou n peut être égal 317 et a peut encore être égal à 17. Pour les figures qui sont présentées dans la suite de cette description, et pour obtenir une illustration 2952722 -15- claire du principe de l'invention, n est pris égal à 7 et a égal à 3, 3 étant une racine primitive du groupe-quotient multiplicatif (Z/7Z)*. Les valeurs du facteur de puissance p et celles du facteur multiplicatif q forment deux séquences qui sont répétées avec une durée de période égale à 6 x Tmod. 5 En général, la durée Tmoa des créneaux temporels de la modulation est choisie de sorte que son inverse 1 /Tmoa appartienne au domaine des radiofréquences (RF). La détection hétérodyne qui est réalisée par l'unité de mélange-détection 11 comprend alors une étape de démodulation RF. Dans le premier mode de mise en oeuvre de l'invention qui est décrit 1 o maintenant, en référence aux figures 2a et 2b, la caractéristique de phase Cphase qui est modulée de la façon qui vient d'être décrite, est la fréquence f. Ainsi, l'invention introduit une modulation de fréquence pour l'onde optique du signal d'émission SE. Dans ce cas, la relation (1) devient : Af=Èfoxq (1') 15 où Afo est un incrément fixe de fréquence. Par exemple, Afo peut être de l'ordre de 10-8 x f, f étant la fréquence de l'onde optique OL qui est produite par l'oscillateur laser 1. Le diagramme de la figure 2a montre, en fonction du temps noté t, les décalages de fréquence respectifs du signal de référence SRef d'une part, de la 20 contribution utile RV au signal rétrodiffusé R qui est produite par le contenu du volume-cible V d'autre part, et de la contribution parasite RP au même signal rétrodiffusé R par ailleurs, tels que ces décalages sont synchronisés pour la détection hétérodyne. La fréquence de la contribution utile RV est modulée à chaque instant comme celle du signal de référence SRef, avec un écart AfDoppler 25 par rapport à celle-ci qui correspond à l'effet Doppler qui est provoqué par le déplacement du contenu du volume-cible V. A la différence, la fréquence de la contribution parasite RP est modulée comme le signal de référence SRef, mais avec un retard qui est égal au temps supplémentaire de propagation pour la contribution parasite RP. Ce temps supplémentaire de propagation est 30 2 x Ad/C, et est pris égal à 2•Tmod pour les figures 2a et 2b, à titre d'illustration. En outre, on a supposé pour ces figures que la contribution parasite RP ne présentait pas d'effet Doppler dû à un déplacement de la source parasite P. 2952722 -16- Néanmoins, un tel déplacement de la source parasite P ne modifie pas l'efficacité de l'invention pour isoler la contribution utile RV. Le diagramme de la figure 2b montre, encore en fonction du temps t, les fréquences de battement qui apparaissent lors de la détection hétérodyne 5 pour la contribution utile RV et pour la contribution parasite RP. Il correspond aux décalages de fréquence Af qui sont indiqués dans la figure 2a. Ces fréquences de battement correspondent à l'écart à chaque instant, entre les fréquences respectives de la contribution utile RV et du signal de référence SRef d'une part, et entre les fréquences respectives de la contribution parasite 1 o RP et du signal de référence SRef d'autre part. Elles peuvent être déterminées numériquement. Pour la contribution utile RV, la fréquence du battement est constante sur toute la durée du cycle de mesure, et est égale à Afpopp,er. Grâce à l'invention, la fréquence du battement pour la contribution parasite RP prend les valeurs -3 x Afo, -2 x Afo, - Afo, +Afo, +2 x Afo et +3 x Afo, dans un ordre qui 15 dépend de n, de a et de l'ordre des valeurs du facteur de puissance p. L'ordre des valeurs successives de la fréquence de battement pour la contribution parasite RP n'a aucune fonction d'information de mesure. Ainsi, lorsque le signal de détection hétérodyne SRF est analysé spectralement, la contribution parasite RP est dispersée sur les fréquences -3 x Afo, -2 x Afo, - Afo, +Afo, 20 +2 x Afo et +3 x Afo, alors que la contribution utile RV reste concentrée à Afoopp,er. Cette dernière apparaît donc avec une amplitude prépondérante dans le résultat de l'accumulation. Elle est isolée puis analysée de façon usuelle pour obtenir le résultat de mesure vélocimétrique ou vibrométrique. Ce mode de mise en oeuvre de l'invention par modulation de la 25 fréquence est particulièrement adapté pour mesurer une vitesse anémométrique. Il peut alors être mis en oeuvre à bord d'un aéronef, notamment pour fournir une mesure de référence de la vitesse du vent. En particulier, la vitesse du vent peut ainsi être mesurée à plus de 400 m (mètre) de l'aéronef, de sorte que l'écoulement d'air dans la portion-cible V, à l'origine 30 du résultat de la mesure, n'est pas perturbé par le déplacement de l'aéronef lui-même. Pour cela, la durée individuelle Tmoa des créneaux temporels de modulation peut être comprise entre 0,2 ps (microseconde) et 10 µs. Par exemple, Tmoa peut être prise égale à 2 µs, ce qui permet de séparer dans le 2952722 -17- signal rétrodiffusé R la contribution utile qui provient d'une distance sélectionnée par rapport à des contributions parasites qui proviennent d'autres distances, différentes de plus de 300 m par rapport à la distance sélectionnée. Dans une telle application de mesure anémométrique, les objets à 5 l'origine de la contribution utile RV sont des poussières ou des gouttelettes qui sont situées à l'intérieur du volume-cible V. Leur contribution RV au signal rétrodiffusé R peut alors être considérablement moins intense qu'une contribution parasite RP d'un nuage plus éloigné, mais qui est aussi situé dans la direction de pointage P. Grâce à l'invention, la contribution utile RV peut néanmoins être séparée et analysée pour en déduire une mesure précise de la vitesse du vent à la distance qui est sélectionnée à partir de l'aéronef. Dans les seconds modes de mise en oeuvre de l'invention qui sont décrits ci-dessous, la caractéristique de phase Cphase qui est modulée selon l'invention est la phase initiale (p. Ainsi, l'onde optique du signal d'émission SE est modulée en phase par l'invention, et la relation (1) devient : L p=A~Poxq (1") où A90 est l'incrément fixe de phase initiale. De préférence, cet incrément fixe A90 est pris égal à 2•rr/n, afin que la modulation produise un étalement maximal de la contribution parasite RP. Les figures 3a et 3b sont obtenues de la même façon que les figures 2a et 2b, de sorte que l'Homme du métier comprendra l'illustration correspondante de l'invention sans qu'il soit nécessaire de répéter le mode d'obtention de ces figures. De nouveau, la durée de chaque fenêtre temporelle du découpage d'analyse est plus longue que celle de chaque créneau temporel de la modulation, cette modulation portant maintenant sur la phase initiale (p. L'effet Doppler de la contribution utile RV apparaît maintenant comme un décalage de phase supplémentaire, qui croît linéairement et qui s'ajoute au décalage Ècp de la modulation de la phase initiale pour la contribution utile RV. AfDoppler est alors la pente de la variation de Lcp à l'intérieur de chaque créneau temporel de modulation. Afin d'obtenir une mesure constante de l'écart de phase qui est produit par l'effet Doppler, notamment à la fin de chaque fenêtre temporelle du découpage d'analyse, il est préférable que la durée individuelle 2952722 -18- Tmod des créneaux temporels de modulation soit un diviseur de la durée de la fenêtre temporelle du découpage d'analyse, qui est notée TdéC. Le signal de détection hétérodyne SRF est d'abord démodulé par rapport aux sauts de phase initiale, puis analysé spectralement. Lors de cette analyse spectrale, la 5 contribution utile RV n'est pas atténuée, alors que la contribution parasite RP est combinée avec une fonction de sinus cardinal au carré, dont la largeur est égale à 1 /Tmod. Ainsi, ce type de modulation de phase produit encore une dispersion spectrale, sélectivement pour la contribution parasite RP. Les résultats de cette analyse spectrale pour les fenêtres successives de 10 découpage temporel sont cumulés pour isoler la contribution utile RV, puis l'analyse d'effet Doppler fournit le résultat de la mesure vélocimétrique ou vibrométrique. Un tel mode de mise en oeuvre de l'invention par modulation de phase est particulièrement adapté pour mesurer en laboratoire une vitesse 15 d'écoulement. Elle permet en effet de distinguer efficacement la contribution de l'écoulement dans le signal rétrodiffusé R qui est recueilli, par rapport à une contribution de réflexion du signal d'émission qui est produite par une paroi fixe d'arrière-plan. Enfin, dans une variante de la mise en oeuvre de l'invention par 20 modulation de phase, la durée individuelle Tmod des créneaux temporels de modulation est supérieure à la durée individuelle TdéC de la fenêtre temporelle du découpage d'analyse. Le signal de détection hétérodyne SRF, démodulé par rapport aux sauts de phase initiale, est d'abord découpé en segments qui correspondent aux fenêtres temporelles. Ces segments sont ensuite 25 accumulés, c'est-à-dire additionnés les uns aux autres, sur toute la durée du cycle de mesure. Cette accumulation est constructive pour la contribution utile RV et destructive pour la contribution parasite RP. De préférence, la durée individuelle Trima des créneaux temporels de modulation est un multiple de la durée individuelle TdéC de la fenêtre temporelle du découpage d'analyse. 30 Autrement dit, Tmod = z x Tdéc, où z est un nombre entier strictement positif. En effet, la phase initiale est alors constante à l'intérieur de chaque fenêtre temporelle. Dans ce cas, l'accumulation de toutes les valeurs de la phase initiale cp pour la contribution parasites RP est égale à ûz, alors que la même 2952722 -19- accumulation pour la contribution utile RV est égale à z x (n-1). Pour cette variante de mise en oeuvre, l'invention produit plutôt une atténuation de la contribution parasite RP, plutôt qu'une dispersion spectrale de celle-ci. L'analyse spectrale puis l'analyse d'effet Doppler sont effectuées finalement 5 sur le résultat de l'accumulation, pour obtenir le résultat de la mesure vélocimétrique ou vibrométrique. Il est entendu que l'invention peut être appliquée en introduisant des adaptations par rapport aux modes de mise en oeuvre qui ont été décrits en détail ci-dessus. En particulier, la modulation qui est introduite par l'invention 10 peut être appliquée à une onde optique qui est continue comme cela a été décrit, ou bien qui est produite sous forme d'impulsions successives. En outre, l'invention peut être utilisée non seulement pour réaliser une mesure de vitesse du contenu du volume-cible, mais aussi pour mesurer d'autres caractéristiques vélocimétriques ou vibrométriques, telles qu'une largeur de distribution de 15 vitesse, une fréquence de vibration, une amplitude de vibration, etc. REVENDICATIONS1. Procédé de mesure d'une caractéristique vélocimétrique ou vibrométrique utilisant un dispositif de type LIDAR à détection hétérodyne (100), dans lequel un signal d'émission d'onde optique (SE) est produit en direction d'un volume-cible (V) à partir d'une tête optique dudit dispositif, et un signal rétrodiffusé (R) est recueilli par ladite tête optique, puis ledit signal rétrodiffusé est détecté par détection hétérodyne de façon à produire un signal de détection hétérodyne (SRF), ledit procédé comprenant les étapes suivantes pour un cycle de mesure : /1/ produire une modulation d'une caractéristique de phase de l'onde optique dans le signal d'émission (SE) ; /2/ démoduler le signal de détection hétérodyne (SRF) par rapport à la modulation de la caractéristique de phase, en compensant un retard de propagation du signal d'émission (SE) et du signal rétrodiffusé (R) entre la tête optique du dispositif (100) et le volume-cible (V) ; /3/ combiner une analyse spectrale du signal de détection hétérodyne (SRF) démodulé avec une accumulation pour des fenêtres temporelles successives de découpage (TdéC), de façon à isoler une contribution au signal de détection hétérodyne provenant du volume-cible (V) ; et /4/ obtenir un résultat de la mesure de caractéristique vélocimétrique ou vibrométrique, à partir d'une analyse d'effet Doppler de la contribution isolée à l'étape /3/, et dans lequel : - la modulation est obtenue en décalant la caractéristique de phase de l'onde optique pendant des créneaux temporels successifs de modulation (Tmod), conformément à un pas d'incrément fixe multiplié par des facteurs multiplicatifs q respectivement affectés aux dits créneaux temporels de modulation ; les facteurs multiplicatifs q sont égaux à ap[n], où : 2952722 - 21 - n est un nombre entier supérieur ou égal à quatre, constant pendant le cycle de mesure, a est un autre nombre entier, aussi constant pendant le cycle de mesure, 5 p est un facteur de puissance entier positif ou nul, strictement inférieur à n-1, et avec des valeurs variables respectivement affectées aux créneaux temporels de modulation, ap[n] désigne la réduction de ap modulo n, et a est choisi tel que a[n] soit différent de zéro et de l'unité. 10 2. Procédé selon la revendication 1, suivant lequel les facteurs multiplicatifs q des créneaux temporels successifs de modulation (Tmod) forment une séquence qui est répétée périodiquement pendant le cycle de mesure. 3. Procédé selon la revendication 1 ou 2, suivant lequel les facteurs multiplicatifs q successifs comprennent au moins trois valeurs différentes. 15 4. Procédé selon la revendication 3, suivant lequel les facteurs multiplicatifs q successifs comprennent au moins cent valeurs différentes. 5. Procédé selon l'une quelconque des revendications précédentes, suivant lequel le nombre n est premier. 6. Procédé selon la revendication 5, suivant lequel le nombre entier a 20 est une racine primitive d'un groupe-quotient multiplicatif (Z/nZ)*, où Z est l'anneau des nombres entiers relatifs, * indique que la valeur zéro est exclue, les facteurs multiplicatifs q étant n-1 valeurs distinctes pour p variant de 0 à n-2. 7. Procédé selon l'une quelconque des revendications 1 à 6, suivant 25 lequel la caractéristique de phase de l'onde optique qui est décalée est une fréquence (f) de ladite onde optique, de façon à réaliser une modulation de fréquence, et suivant lequel une durée individuelle des créneaux temporels de modulation (Tmod) est inférieure ou égale à une durée individuelle des fenêtres 2952722 - 22 - temporelles de découpage (Taéo), l'étape /3/ comprenant les sous-étapes suivantes : /3-1/ analyser spectralement le signal de détection hétérodyne (SRF) démodulé à l'intérieur d'une fenêtre temporelle de découpage (Taéo), de 5 façon à obtenir un résultat d'analyse spectrale pour ladite fenêtre temporelle ; et /3-2/ répéter la sous-étape /3-1/ pour plusieurs fenêtres temporelles successives de découpage, et accumuler les résultats d'analyse spectrale respectivement obtenus pour lesdites fenêtres temporelles, la 1 o contribution au signal de détection hétérodyne (SRF) provenant du volume-cible (V) étant prépondérante dans l'accumulation. 8. Procédé selon la revendication 7, utilisé pour mesurer une vitesse anémométrique. 9. Procédé selon la revendication 7 ou 8, suivant lequel la durée 15 individuelle des créneaux temporels de modulation (Tmoa) est comprise entre 0,2 ps et 10 µs. 10. Procédé selon l'une quelconque des revendications 7 à 9, mis en oeuvre à bord d'un aéronef. 11. Procédé selon l'une quelconque des revendications 1 à 6, suivant 20 lequel la caractéristique de phase de l'onde optique qui est décalée est une phase initiale (cp) de ladite onde optique, de façon à réaliser une modulation de phase. 12. Procédé selon la revendication 11, suivant lequel l'incrément fixe (A(po) de la phase initiale (cp) de l'onde optique est égal à 2•rr/n. 25 13. Procédé selon la revendication 11 ou 12, suivant lequel une durée individuelle des créneaux temporels de modulation (Tmod) est inférieure ou égale à une durée individuelle des fenêtres temporelles de découpage (Tdéc), l'étape /3/ comprenant les sous-étapes suivantes : 2952722 - 23 - /3-1/ analyser spectralement le signal de détection hétérodyne (SRF) démodulé à l'intérieur d'une fenêtre temporelle de découpage (Tdéc), de façon à obtenir un résultat d'analyse spectrale pour ladite fenêtre temporelle ; et 5 /3-2/ répéter la sous-étape /3-1/ pour plusieurs fenêtres temporelles successives de découpage, et accumuler les résultats d'analyse spectrale respectivement obtenus pour lesdites fenêtres temporelles, la contribution au signal de détection hétérodyne (SRF) provenant du volume-cible (V) étant prépondérante dans l'accumulation. 10 14. Procédé selon la revendication 13, suivant lequel la durée individuelle des créneaux temporels de modulation (Tmod) est un diviseur de la durée individuelle des fenêtres temporelles de découpage (Tdéo). 15. Procédé selon la revendication 11 ou 12, suivant lequel une durée individuelle des créneaux temporels de modulation (Tmod) est supérieure ou 15 égale à une durée individuelle des fenêtres temporelles de découpage (Tdéc), l'étape /3/ comprenant les sous-étapes suivantes : /3-1/ accumuler des segments du signal de détection hétérodyne (SRF) démodulé correspondant à des fenêtres temporelles successives de découpage (Tdéc), pendant les créneaux temporels successifs de 20 modulation (Tmod) de sorte que la contribution au signal de détection hétérodyne (SRF) provenant du volume-cible (V) soit accumulée de façon constructive ; et /3-2/ analyser spectralement un résultat de l'accumulation. 16. Procédé selon la revendication 15, suivant lequel la durée 25 individuelle des créneaux temporels de modulation (Tmod) est un multiple de la durée individuelle des fenêtres temporelles de découpage (Tdéo). 17. Dispositif de type LIDAR (100) adapté pour effectuer des mesures vélocimétriques ou vibrométriques, et comprenant : - un oscillateur laser (1), adapté pour produire une onde optique (OL) ; 2952722 - 24 - - un séparateur optique (2), disposé pour diviser l'onde optique en un signal de source d'émission (SSE) et un signal de référence (SRef) - un amplificateur optique (3), adapté pour produire un signal d'émission (SE) à partir du signal de source d'émission ; 5 - une tête optique (12), adaptée pour transmettre le signal d'émission en direction d'un volume-cible (V) et pour recevoir un signal rétrodiffusé (R) ; et - une unité de mélange-détection (11), adaptée pour transmettre le signal d'émission à la tête optique, et pour produire un signal de 1 o détection hétérodyne (SRF) à partir du signal rétrodiffusé reçu par ladite tête optique et du signal de référence ; et - des moyens d'analyse (20) du signal de détection hétérodyne, adaptés pour combiner une analyse spectrale dudit signal de détection hétérodyne avec une accumulation pour des fenêtres temporelles 15 successives de découpage (Tdéc), et pour effectuer une analyse d'effet Doppler, le dispositif étant caractérisé en ce qu'il comprend en outre : - un modulateur de phase (4) disposé pour moduler au moins le signal de source d'émission ; 20 - une unité de commande (40), reliée à une entrée de commande du modulateur de phase, et adaptée pour commander un fonctionnement dudit modulateur de phase pour mettre en oeuvre un procédé selon l'une quelconque des revendications précédentes ; et - des moyens de compensation d'un retard de propagation du signal 25 d'émission (SE) et du signal rétrodiffusé (R) entre la tête optique (12) et le volume-cible (V). 18. Dispositif selon la revendication 17, dans lequel le modulateur de phase (4) comprend une cellule de Pockels. 19. Dispositif selon la revendication 17 ou 18, dans lequel le modulateur de phase (4) est disposé pour recevoir en entrée le signal de source d'émission 2952722 - 25 - (SSE) produit par le séparateur optique (2), et pour transmettre en sortie ledit signal de source d'émission modulé à l'amplificateur optique (3), le dispositif comprenant en outre des moyens de démodulation (5) agencés pour démoduler le signal de détection hétérodyne (SRF) conformément au 5 fonctionnement du modulateur de phase (4), et les moyens de compensation de retard étant agencés pour compenser le retard de propagation du signal d'émission (SE) et du signal rétrodiffusé (R) lors de la démodulation. 20. Dispositif selon la revendication 17 ou 18, dans lequel le modulateur de phase (4) est disposé pour recevoir en entrée l'onde optique (OL) produite 1 o par l'oscillateur laser (1), et pour transmettre en sortie ladite onde optique modulée au séparateur optique (2), de sorte que le signal de source d'émission (SSE) et le signal de référence (SRef) sont modulés d'une façon identique, les moyens de compensation de retard (6) étant agencés pour retarder le signal de référence (SRef) transmis à l'unité de mélange-détection (11). 15