FR 3025892 A1 20160318 FR 1458746 A 20140916 1 DOMAINE TECHNIQUE GENERAL L'invention concerne le domaine technique de la surveillance radar et en particulier la surveillance de sites instables. En particulier l'invention est relative aux procédés et dispositifs pour la surveillance de sites instables rocheux. ETAT DE L'ART Par "site instable", on désigne ici et dans toute la suite du texte toute zone terrestre qui fait l'objet de modifications dans le temps, que cette zone soit naturelle ou non, que ces modifications correspondent à des mouvements de grande échelle (fissuration d'un sol par exemple) ou qu'elles soient localisées dans le temps ou en certaines parties de la zone (chutes de blocs rocheux par exemple). L'invention trouve notamment avantageusement application en volcanologie et sismologie. Elle trouve également avantageusement application pour la surveillance de zones de sol naturel présentant des risques importants de chutes de blocs rocheux. Un exemple d'utilisation sera pris avec les lignes de chemins de fer surplombées par des sites rocheux instables et sur lesquelles circulent notamment des trains à grande vitesse. Lorsqu'un rocher se détache du site et vient encombrer la voie, il existe un risque important d'endommagement de la locomotive et des rails. Il est ainsi souhaitable de savoir à chaque instant si la voie est dégagée. Des dispositifs en forme de filets maillés en fil conducteur de courant sont couramment observés aux abords des chemins de fer. Le maillage de ces derniers, typiquement de 15 cm, se casse lorsqu'un rocher supérieur à la dimension du maillage vient au contact du dispositif. Une coupure de courant est ainsi observée, ce qui permet d'alerter un opérateur quant à la présence d'un rocher supérieur à 15 cm à un 3025892 2 certain endroit sur les voies. La valeur de 15 cm est la taille maximale acceptable pour que la locomotive puisse pulvériser le rocher sans être endommagée. Un rocher plus grand impose une intervention pour dégager la voie. 5 Malheureusement quand un tel rocher est détecté trop tard et qu'un train se dirige vers la zone encombrée, un freinage d'urgence peut être déclenché afin de limiter les dommages de la locomotive. Mais un freinage d'urgence provoque une dégradation de la voie. Il est nécessaire de pouvoir détecter le plus tôt possible si un rocher 10 nécessitant une intervention d'un opérateur est sur la voie. En ce qui concerne l'obtention de la trajectoire, il existe d'autres dispositifs. Des distances-mètres optiques ont été utilisés en sismologie et volcanologie mais avec leur mire coûteuse et fragile, ils ne sont pas 15 adaptés pour être laissés à demeure sur les sites et en outre ils ne permettent que des mesures ponctuelles dans le temps et l'espace et ne donnent aucune information sur la taille du rocher. Le document US 4 300 121 propose une surveillance radar d'un site 20 rocheux instable en analysant les échos des ondes électromagnétiques émises par le radar et réfléchies par le rocher mais il ne permet pas de suivre l'évolution d'un point particulier ni d'avoir une estimation de la dimension. 25 Enfin, le brevet EP 0 811 854, à l'aide d'un radar également, permet d'établir la trajectoire d'un rocher mais il ne permet pas d'obtenir la taille de ce dernier. Or il existe un besoin d'identifier le plus rapidement possible si un rocher de taille supérieure à une dimension donnée présente un danger. 30 3025892 3 PRESENTATION DE L'INVENTION L'invention permet de pallier les inconvénients ci-dessus mentionnés et propose à cet effet un procédé de surveillance d'une cellule élémentaire 5 d'un site, au moyen d'un dispositif de surveillance comprenant un émetteur et un récepteur d'ondes électromagnétiques, le procédé comprenant les étapes suivantes mises en oeuvre dans une unité de traitement : a) émission vers la cellule élémentaire à surveiller, au moyen de 10 l'émetteur, d'au moins trois ondes électromagnétiques de fréquences respectives différentes et telles que la longueur d'onde associée à une première fréquence de battement définie par la différence entre une première et une deuxième desdites fréquences est supérieure à la taille radiale de la cellule élémentaire; 15 b) réception cohérente, au moyen du récepteur, des ondes réfléchies issues de la réflexion des ondes sur la cellule surveillée ; c) traitement des ondes réfléchies pour détecter un mouvement d'un objet dans la cellule élémentaire et estimer sa position dans la cellule élémentaire, ledit traitement mettant en outre en oeuvre, en fonction de 20 la longueur d'onde associée à une deuxième fréquence de battement définie par la différence entre la troisième des trois fréquences et l'une des deux autres, une caractérisation d'une dimension dudit objet. Un but de l'invention est ainsi de pouvoir caractériser la dimension d'un 25 rocher, en plus de connaître sa trajectoire. Il est ainsi possible d'estimer dès le départ de la chute du rocher, dans le cas des voies de chemins de fer, s'il faudra effectuer un freinage d'urgence ou non. Le procédé selon l'invention peut comprendre les caractéristiques 30 suivantes prises seules ou en combinaison : - ledit traitement est adapté pour caractériser la dimension dudit objet par rapport à une dimension seuil comprise entre la moitié 3025892 4 de la longueur d'onde associée à la deuxième fréquence de battement et deux fois celle-ci, - ledit traitement est adapté pour estimer la dimension dudit objet 5 si ladite dimension est comprise entre la moitié de la longueur d'onde associée à la deuxième fréquence de battement et deux fois celle-ci, - on émet, outre les trois ondes, une quatrième onde dont la 10 fréquence est telle que la fréquence de battement définie par la différence entre la troisième fréquence et cette quatrième fréquence est égale à la première fréquence de battement, la fréquence de battement définie comme la différence entre cette quatrième fréquence et l'une des deux premières étant égale à la 15 deuxième fréquence de battement, le procédé comprenant la réception de l'onde réfléchie par le récepteur après réflexion de l'onde émise, le traitement mettant en oeuvre une vérification de ce que le mouvement détecté est relatif à un unique objet, 20 - la vérification de ce que le mouvement détecté est relatif à un unique objet met en oeuvre une comparaison des corrélations entre les ondes réfléchies correspondant à la première fréquence de battement et/ou les ondes réfléchies correspondant à la deuxième fréquence de battement, 25 - le traitement des ondes réfléchies comprend les étapes suivantes: Cl) filtrage passe-haut des ondes réfléchies pour supprimer des composantes desdites ondes réfléchies sur au moins un objet(s) fixe(s) ou quasi-fixe(s) dans la cellule élémentaire, 30 C2) détermination de corrélations entre les ondes réfléchies pour en déduire la position de l'objet en mouvement ainsi que la dimension dudit objet en mouvement, 3025892 5 on détermine l'argument d'une corrélation entre deux ondes correspondant à la première fréquence de battement pour en déduire la position de l'objet en mouvement, 5 on détermine le module d'une corrélation entre deux ondes correspondant à la deuxième fréquence de battement pour en déduire la dimension de l'objet en mouvement, - la deuxième fréquence de battement est comprise entre 10 MHz 10 et 10 GHz, la longueur d'onde associée étant comprise entre 15 mm et 15 m, - la première fréquence de battement est comprise entre 0,1 MHz et 100 MHz, la taille radiale de la cellule élémentaire à surveiller 15 étant comprise entre 1m et 1km. L'invention propose aussi un dispositif de surveillance d'une cellule élémentaire comprenant : - un émetteur d'ondes électromagnétiques, configuré pour émettre 20 des ondes électromagnétiques, - un récepteur cohérent d'ondes électromagnétiques, configuré pour recevoir des ondes électromagnétiques réfléchies, une unité de traitement configurée pour mettre en oeuvre un procédé tel que décrit précédemment. 25 PRESENTATION DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, 30 et qui doit être lue en regard des dessins annexés, sur lesquels : - La figure 1 illustre un site instable à surveiller au moyen d'un procédé et un dispositif conformes à l'invention ; 3025892 6 - La figure 2 illustre en trois dimensions les formes géométriques de cellules élémentaires d'un site instable à surveiller au moyen d'un procédé et un dispositif conformes à l'invention ; - La figure 3 illustre des propriétés géométriques d'un objet en 5 mouvement dans le site à surveiller ; - La figure 4 illustre des étapes d'un procédé de surveillance d'un site conforme à l'invention selon un premier mode de réalisation ; - La figure 5 illustrent un axe sur lequel sont réparties des fréquences f1, f2 et f3 des ondes électromagnétiques impliquées 10 dans un procédé selon un premier mode de réalisation de l'invention ; - Les figures 6a et 6b illustrent des courbes représentant une dimension de l'objet en mouvement du site à surveiller en fonction du module du coefficient de corrélation entre des ondes 15 impliqués dans un procédé conforme à l'invention; - La figure 7 illustre un axe sur lequel sont réparties les fréquences f2, f3 et f4 des ondes électromagnétiques impliquées dans un procédé selon un second mode de réalisation de l'invention. 20 DESCRIPTION DETAILLEE La figure 1 illustre une surveillance d'un site comprenant un objet, typiquement un rocher 40 (ou bloc rocheux) que l'on surveille au moyen d'un dispositif 1 de surveillance. 25 Le dispositif 1 Le dispositif 1 comprend un émetteur d'ondes électromagnétiques 10, configuré pour émettre des ondes électromagnétiques, un récepteur 30 d'ondes électromagnétiques 20, configuré pour recevoir des ondes électromagnétiques réfléchies, et une unité de traitement 30. Afin de surveiller le site et notamment le rocher 40, on dispose le dispositif 1 à un point d'observation O de sorte qu'il puisse émettre des ondes 3025892 7 électromagnétiques vers le site. Il est néanmoins possible de placer l'émetteur 10 et le récepteur 20 à des endroits différents. L'émetteur d'ondes électromagnétiques 10 et le récepteur 20 forment généralement un radar monostatique. En particulier, l'émetteur 10 est 5 compatible d'une émission large bande. Il comporte une ou plusieurs antennes d'émission et de réception à réflecteur parabolique. Chaque antenne a un gain d'environ 30 dB et un angle d'ouverture à 3dB de l'ordre de 1,7°. Le champ éclairé par le radar à 1km a un diamètre d'environ 30m. Une lunette de visée facilite le pointage du radar vers le 10 site à surveiller. L'émetteur 10 peut émettre différentes formes d'ondes connues de l'homme du métier, par exemple multifréquence ou avec modulations linéaires de fréquence montantes et descendantes. 15 Les ondes électromagnétiques E1, E2 ou E3, de fréquences respectives f1, f2, f3 différentes, émises par l'émetteur 10 sont réfléchies par des éléments du site, notamment le rocher 40, pour revenir sous forme d'ondes réfléchies S1, S2 ou S3 au récepteur 20. L'unité de traitement 30 est configurée pour déclencher les émissions 20 d'ondes électromagnétiques et pour effectuer l'analyse des ondes reçues. Le site à surveiller 25 Comme illustré sur la figure 2, le site à surveiller est décomposé en un ensemble de cellules élémentaires S, chaque cellule élémentaire étant un volume défini par une largeur moyenne I, une hauteur moyenne h et une taille radiale r (ou profondeur). Les cellules élémentaires S sont typiquement des cases distance, selon la terminologie consacrée. La 30 largeur et la hauteur d'une cellule élémentaire S surveillée dépendent des angles d'émission utilisés par l'émetteur 10. La taille radiale r (ou profondeur) est liée à la fenêtre temporelle de réception du récepteur 20, i.e. au début de l'acquisition et à la durée d'acquisition. Les cellules 3025892 8 élémentaires ou cases distance peuvent se recouvrir partiellement. Les signaux issus des différentes cases distance S sont traités en parallèle. Pour bien localiser la position radiale de l'objet en mouvement à 5 l'intérieur de la case distance S, il faut que la profondeur r de la case distance S soit inférieure à la longueur d'onde associée au battement de la paire de fréquences utilisée (f2-f1 pour la figure 5). De cette façon, la phase de la fonction d'intercorrélation ne peut pas effectuer plus d'un tour (360°) à l'intérieur de la case distance S, ce qui assure l'unicité de 10 la valeur de la position radiale. L'objet surveillé En relation avec la figure 3 et comme mentionné précédemment, l'objet 15 40 surveillé est typiquement le rocher 40. Il peut néanmoins s'agir de n'importe quel élément pouvant se mettre en mouvement, le mouvement pouvant être nul, et dit diffractant, c'est-à-dire pouvant réfléchir les ondes électromagnétiques. Quand le rocher 40 se met en mouvement dans le site et donc dans une case distance S, le 20 mouvement étant typiquement une chute, il se met en rotation autour d'un axe de rotation. On peut ainsi lui définir une dimension 41 qui est la taille maximale du rocher perpendiculairement autour dudit axe de rotation. Cette dimension 41 est généralement la taille maximale du rocher 40, c'est-à-dire la distance maximale reliant deux points du 25 rocher 40. On définit la vitesse v du rocher 40 comme la vitesse de son centre de gravité et on définit la position radiale de l'objet 40 ou du rocher 40 comme la composante de la position de l'objet projetée sur l'axe entre le récepteur 20 et l'objet 40. Le suivi dans le temps de la position radiale 30 de l'objet 40 s'effectue sur l'ensemble des cases distance S et permet donc de savoir si l'objet se rapproche ou non du récepteur 20. 3025892 9 Procédé de détection A présent, un procédé de détection de mouvement d'objets diffractant dans un site à surveiller va être décrit en relation avec la figure 4. 5 Dans une première étape Pl, l'émetteur 10 émet en direction d'une case distance au moins trois ondes électromagnétiques E1, E2, E3, de fréquences respectives f1, f2, f3, différentes. Les ondes sont émises simultanément ou successivement, potentiellement sous forme d'une 10 onde complexe superposant les trois ondes E1, E2, E3. On définit une première fréquence de battement par la différencelf2-fil entre une première et une deuxième desdites fréquences. Cette première fréquence de battement est telle que la longueur d'onde 15 associée est supérieure ou égale à la taille radiale r de la case distance S utilisée. En effet, ainsi qu'il l'a été expliqué précédemment, à l'intérieur de la case distance S, il ne faut pas que le déphasage entre S1 et S2 varie de plus de 360°. Typiquement, la première fréquence de battement est comprise entre 20 0,1 MHz et 100 MHz, ce qui correspond à une taille radiale r de la profondeur de la case distance comprise entre 1m et 1km. Typiquement, la première fréquence de battement est de l'ordre de 5 MHz, soit une taille radiale r de la case distance S inférieure à 30m. 25 On définit une deuxième fréquence de battement 1f3-filou If3421définie par la différence entre la troisième des trois fréquences et l'une des deux autres. Cette deuxième fréquence de battement est telle que la longueur d'onde associée, appelée dimension prédéterminée DP, est caractéristique 30 d'une dimension seuil 42. Par prédéterminée, l'on entend qu'elle est déterminée au préalable. Dans la suite de la description, la longueur d'onde associée à la deuxième fréquence de battement sera toujours appelée distance prédéterminée. Cette dimension prédéterminée DP est 3025892 10 fixée de manière à centrer le domaine exploitable (région II de la figure 6b). Typiquement, la deuxième fréquence de battement est comprise entre 10 MHz et 10 GHz, ce qui permet d'apprécier la dimension de blocs 5 allant de 15 m à 1,5 cm. La figure 5 montre un exemple de répartition de ces fréquences selon l'invention. Dans une deuxième étape P2, le récepteur 20 reçoit les ondes réfléchies 10 Si, 52, S3. Dans une troisième étape P3, l'unité de traitement 30 analyse les ondes réfléchies S1, 52, S3 pour détecter un éventuel mouvement du rocher 40 du site ainsi que la dimension 41 dudit rocher 40. L'obtention des 15 positions radiales successives du rocher 40 permet de connaître sa trajectoire radiale dans case distance et in fine dans le site. Plus précisément, l'étape P3 comprend les sous-étapes suivantes : 20 La sous-étape P31 consiste à filtrer les ondes réfléchies pour supprimer les composantes des ondes qui sont caractéristiques d'objets fixes ou quasi-fixes. En effet, dans le site, se trouvent des objets fixes tels que des éléments rocheux stables ou quasi-fixes, tels que la végétation. Par quasi-fixe, on entend une vitesse de déplacement qui permet de 25 discriminer une chute d'objet de la végétation ou d'animaux, par exemple cette vitesse est inférieure à quelques mètres par seconde, typiquement 3m/s. Le filtrage se fait typiquement à l'aide d'un filtre passe-haut. Le filtrage permet ainsi d'éliminer le bruit. 30 La sous-étape P32 consiste à effectuer des corrélations entre les différentes ondes réfléchies Si, 52, S3 pour obtenir la position, puis la trajectoire, dans la case distance S, du rocher 40 et la dimension 41 du rocher 40. Plus précisément, l'argument d'une moyenne temporelle de 3025892 11 corrélation entre deux ondes correspondant la première fréquence de battement (ici Si et S2) permet de connaître les positions radiales successives et donc la trajectoire radiale de l'éventuel rocher 40 en mouvement et le module d'une moyenne temporelle de corrélation entre 5 deux ondes correspondant la deuxième fréquence de battement (ici Si et S3, S2 et S3) permet de caractériser la dimension 41 du rocher 40 en mouvement grâce à la dimension prédéterminé DP. On connaît la phase en évaluant Arg( ), où « Arg » est l'argument, « <> » représente une moyenne temporelle, Si et Si sont les 10 ondes réfléchis sous forme complexe et « * » symbolise le conjugué complexe. En effet, la position radiale du rocher 40 est obtenue par technique classique d'interférométrie dans le domaine fréquentiel. Lorsque l'on 15 émet les deux ondes électromagnétiques E1, E2 vers le rocher 40 en mouvement, la phase du produit de corrélation (correspondant au déphasage) mesuré entre les deux ondes de retour Si, S2 est proportionnelle à la première fréquence de battement f2-f1 et au déplacement radial du rocher 40 en mouvement. Ainsi, si l'on ne peut 20 garantir que la phase est comprise entre 0 et 360°, il existe une ambiguïté sur la position radiale du rocher 40. En ce qui concerne la caractérisation de la dimension 41 du rocher 40, la moyenne temporelle de la corrélation entre deux signaux correspondant 25 à la deuxième fréquence de battement est normalisée par la racine du produit des moyennes temporelles des autocorrélations de deux ondes Si, Si pour définir le coefficient C. Ainsi, par exemple pour Si et S3, le coefficient de corrélation C13 est défini par : C.1.3 = 30 Les figures 6a et 6b illustrent la courbe du module du coefficient de corrélation de deux ondes correspondant à la deuxième fréquence de 3025892 12 battement (C13, C23 par exemple) en fonction de la dimension 41 du rocher 40. Cette valeur est continue et décroissante. Le coefficient C13 tend vers 0 (ondes non corrélées) quand la dimension 41 du rocher 40 croit au-delà de la dimension prédéterminée DP et ce module tend vers 5 1 (ondes corrélées) quand la dimension 41 du rocher 40 décroît en deçà de la dimension prédéterminée DP. La courbe présente des effets de seuil en extrémité et une plage de décroissance manifeste se trouve entre DP/Y et X.DP, avec typiquement X=Y=2. Dans cette zone la courbe de la figure 6b est exploitable car il est nécessaire qu'il reste 10 une pente assez forte. Plus précisément, le module du coefficient de corrélation entre deux ondes correspondant à la deuxième fréquence de battement (S1 et S3/ S2 et S3 par exemple) est égal à environ 0,5 quand la dimension 41 du rocher 40 est sensiblement égale à la dimension prédéterminée DP. La 15 valeur du module du coefficient par rapport à 0,5 permet de savoir si la dimension 41 du rocher 40 est supérieure ou inférieure à la dimension prédéterminée DP. Il n'est néanmoins pas nécessaire que la dimension seuil 41 soit égale à la dimension prédéterminée DP. On définit trois intervalles dans la courbe 6b, le premier intervalle I 20 allant de 0 à DP/Y, le deuxième intervalle II allant de DP/Y à X.DP, le troisième intervalle III allant de X.DP à l'infini. Dans le premier intervalle I, le module du coefficient est quasi-constant et égal à 1, ce qui permet de conclure que l'objet 40 en mouvement est de dimension 41 largement inférieure à la dimension prédéterminée DP. 25 En revanche, si la dimension seuil 42 se situe dans le deuxième intervalle II, on peut dans ce cas estimer la dimension 41 de l'objet 40 en mouvement à partir du coefficient de corrélation et donc la comparer à la dimension seuil 42. Enfin, dans le troisième intervalle III, de la très faible corrélation, il est possible de conclure que la dimension 41 est 30 largement supérieure à DP. 3025892 13 Pour pouvoir obtenir une estimation de la dimension 41 de l'objet 40 et pour pouvoir comparer des dimensions 41 du même ordre de grandeur que la dimension seuil 42, il est nécessaire de fixer la dimension prédéterminée DP de sorte que la dimension seuil 42 se trouve dans le 5 deuxième intervalle II. Ainsi, d'une façon préférentielle, la dimension prédéterminée (et donc notamment la deuxième fréquence de battement) est choisie de sorte que la dimension seuil 42 se trouve dans le deuxième intervalle II. Par exemple, on peut choisir une dimension prédéterminée DP égale à la 10 dimension seuil 42 car c'est aux alentours de cette valeur que le module du coefficient offre le plus de sensibilité (pente maximale) pour évaluer la dimension 41. Dans l'exemple des voies de chemins de fer, il est nécessaire de pouvoir différencier des objets 40 du même ordre de grandeur (10 cm et 20 cm 15 si la dimension seuil 42 est fixée à 15 cm par exemple). Par erreurs d'incertitude, de mesure et de modèle, on entend que la valeur du module du coefficient de corrélation pour deux ondes correspondant à la deuxième fréquence de battement possède un degré 20 de précision limité dû au bruit sur les mesures et de modèle (un objet 40 possède une forme complexe et ne peut être simplement décrit pas la dimension 41, etc.). Mode de réalisation complémentaire 25 Selon un autre mode de réalisation, complémentaire au procédé précédent, durant l'étape Pl, on émet en outre une quatrième onde E4, à la fréquence f4 définie telle que 1f4-f31 = 1f2-fi1, et donc 1f4-f21 = 1f341.1. La fréquence f4 est représentée en figure 7. 30 Ainsi, les ondes réfléchies Si, S2, S3, S4 forment deux fois deux couples de fréquences : Si, S2 et S3, S4 ont la même première fréquence de battement et Si, S3 et S2, S4 ont la même deuxième fréquence de battement. 3025892 14 Grâce à l'onde S4, on dispose pour la première et la deuxième fréquences de battement d'un couple de ondes réfléchies, à la fois pour évaluer la trajectoire et pour évaluer la dimension 41 du rocher 40. En effet, dans le cas du rocher 40 en mouvement, les modules des 5 coefficients de corrélation des ondes ayant la deuxième fréquence de battement (C13 et C24) sont censés donner deux valeurs similaires puisqu'ils correspondent à la même deuxième fréquence de battement. Si ce n'est pas le cas, alors on peut considérer que le modèle utilisé n'est pas valide, ce qui peut signifier par exemple que plusieurs rochers 10 40 sont en mouvement dans la case distance S. De la même façon, l'analyse des phases des moyennes temporelles des corrélations entre deux ondes ayant la première fréquence de battement (Si, Sz et S3, S4) est censée donner deux valeurs similaires. Si ce n'est pas le cas, comme précédemment, le modèle d'un objet 15 unique en mouvement doit être remis en cause Dans une étape P4, comme décrit précédemment, les phase et modules des corrélations des couples d'ondes réfléchies ayant les mêmes fréquences de battement sont comparés pour évaluer l'unicité de l'objet 40 en mouvement. 20 Il est aussi possible d'émettre d'autres ondes de sorte que les fréquences de battement soient égales à la première ou à la deuxième fréquence de battement, pour multiplier les contrôles et ainsi savoir si les données sont fiables et s'il y a un seul ou plusieurs objets 40 en 25 mouvement. La surveillance l'ensemble des cases distances S d'un site dans une zone permet d'avoir l'ensemble de la trajectoire radiale de l'objet 40 sur toute 30 la zone. En particulier, dans le cas des lignes de chemins de fer, cela permet d'anticiper l'impact du rocher et la gestion de l'alerte. REVENDICATIONS1. Procédé de surveillance d'une cellule élémentaire (S) d'un site, au moyen d'un dispositif (1) de surveillance comprenant un émetteur (10) et un récepteur (20) d'ondes électromagnétiques, caractérisé par les étapes suivantes mises en oeuvre dans une unité (30) de traitement : a) émission vers la cellule élémentaire à surveiller (S), au moyen de l'émetteur (10), d'au moins trois ondes électromagnétiques (E1, Ez, E3) de fréquences respectives (f1, fz, f3) différentes et telles que la longueur d'onde associée à une première fréquence de battement définie par la différence (lfz-fil) entre une première et une deuxième desdites fréquences est supérieure à la taille radiale (r) de la cellule élémentaire (S); b) réception cohérente, au moyen du récepteur (20), des ondes réfléchies (Si, Sz, S3) issues de la réflexion des ondes (E1, E2, E3) sur la cellule surveillée ; c) traitement des ondes réfléchies (Si, Sz, S3) pour détecter un mouvement d'un objet (40) dans la cellule élémentaire (S) et estimer sa position dans la cellule élémentaire (S), ledit traitement mettant en outre en oeuvre, en fonction de la longueur d'onde associée à une deuxième fréquence de battement définie par la différence entre la troisième des trois fréquences et l'une des deux autres, une caractérisation d'une dimension (41) dudit objet (40). 2. Procédé selon la revendication précédente, dans lequel ledit traitement est adapté pour caractériser la dimension (41) dudit objet (40) par rapport à une dimension seuil (42) comprise entre la moitié de la longueur d'onde associée à la deuxième fréquence de battement et deux fois celle-ci. 3025892 16 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel dans lequel ledit traitement est adapté pour estimer la dimension (41) dudit objet (40) si ladite dimension (41) est comprise entre la moitié de la longueur d'onde associée à la 5 deuxième fréquence de battement et deux fois celle-ci. 4. Procédé selon l'une des revendications précédentes, dans lequel on émet, outre les trois ondes (E1, E2, E3), une quatrième onde (E4) dont la fréquence (f4) est telle que la fréquence de 10 battement (Wf4431) définie par la différence entre la troisième fréquence et cette quatrième fréquence est égale à la première fréquence de battement (1f2-fil), la fréquence de battement définie comme la différence entre cette quatrième fréquence et l'une des deux premières étant égale à la deuxième fréquence de 15 battement, le procédé comprenant la réception de l'onde réfléchie (S4) par le récepteur (20) après réflexion de l'onde émise (E4), le traitement mettant en oeuvre une vérification de ce que le mouvement détecté est relatif à un unique objet. 20 5. Procédé selon la revendication précédente, dans lequel la vérification de ce que le mouvement détecté est relatif à un unique objet met en oeuvre une comparaison des corrélations entre les ondes réfléchies correspondant à la première fréquence de battement (Si, Sz ; S3, S4) et/ou les ondes réfléchies 25 correspondant à la deuxième fréquence de battement (Si, S3 et Sz, S4). 6. Procédé selon l'une des revendications précédentes, dans lequel le traitement des ondes réfléchies comprend les étapes 30 suivantes : C1) Filtrage passe-haut des ondes réfléchies pour supprimer des composantes desdites ondes réfléchies sur 3025892 17 au moins un objet(s) fixe(s) ou quasi-fixe(s) dans la cellule élémentaire (S), C2) détermination de corrélations entre les ondes réfléchies pour en déduire la position de l'objet en mouvement ainsi 5 que la dimension (41) dudit objet en mouvement. 7. Procédé selon la revendication précédente, dans lequel : - on détermine l'argument d'une corrélation entre deux ondes correspondant à la première fréquence de battement 10 pour en déduire la position de l'objet (40) en mouvement, - on détermine le module d'une corrélation entre deux ondes correspondant à la deuxième fréquence de battement pour en déduire la dimension (41) de l'objet (40) en mouvement. 15 8. Procédé selon l'une des revendications précédentes, dans lequel la deuxième fréquence de battement est comprise entre 10 MHz et 10 GHz, la longueur d'onde associée (DP) étant comprise entre 15 mm et 15 m. 20 9. Procédé selon l'une des revendications précédentes, dans lequel la première fréquence de battement (If2-fil) est comprise entre 0,1 MHz et 100 MHz, la taille radiale (r) de la cellule élémentaire (S) à surveiller étant comprise entre 1m et 1km. 25 10.Dispositif (1) de surveillance d'une cellule élémentaire comprenant : - un émetteur d'ondes électromagnétiques (10), configuré pour émettre des ondes électromagnétiques, 30 - un récepteur cohérent d'ondes électromagnétiques (20), configuré pour recevoir des ondes électromagnétiques réfléchies, 3025892 18 - une unité de traitement configurée pour mettre en oeuvre un procédé selon l'une des revendications précédentes.