FR 2969768 A1 20120629 FR 1005078 A 20101223 Dispositif d'élimination de perturbations locales pour récepteur de référence de stations sol GNSS La présente invention concerne un dispositif pour éliminer les signaux parasites reçus sur les antennes des récepteurs de référence pour des stations sol GNSS (Global Navigation Satellite Systems). Elle s'applique plus particulièrement dans le cas des infrastructures sol pour les augmentations systèmes de type SBAS (Satellite Based Augmentation Systems), GBAS (Ground Based Augmentation Systems), ainsi que LAAS (Local Area Augmentation System). Les stations GNSS de référence, situées au sol à des positions connues a priori, contribuent à l'amélioration de la précision de positionnement des systèmes de navigation à bas de signaux GNSS. Malheureusement, les mesures fournies par ces stations peuvent être dégradées pour des raisons liées à l'environnement local à la station de référence. Notamment, les réflexions du signal satellite sur les structures situées dans l'environnement local de l'antenne de réception peuvent entraîner des trajets indirects multiples. II peut également contenir des sources d'interférences électromagnétiques, notamment les équipements Radio-Fréquence (RF) situées à proximité de la station. Les trajets indirects multiples et les interférences provoquent des erreurs sur les mesures de phase de code et de porteuse des différents signaux satellites GNSS utilisés. La robustesse des stations de référence aux trajets indirects multiples et aux interférences peut être assurée par utilisation d'une antenne de réception fixe de type FRPA (Fixed Radiated Pattern Antenna). Ces antennes fixes, qu'il s'agisse d'antennes de type « Choke-Ring » ou hélicoïdale, sont conçues dans le but de réaliser par simple filtrage spatial le difficile compromis entre : - la détection et la poursuite des signaux satellites dès les plus basses élévations, et ; - la réjection des trajets indirects multiples et des interférences, tous deux majoritairement situés aux basses élévations. Un inconvénient majeur des antennes fixes de type FRPA est non seulement que ce compromis est très difficile à réaliser à faible élévation, mais également que de telles antennes fixes ne permettent pas l'adaptation du récepteur à l'environnement local de chaque station. En effet, elles ont un diagramme de directivité fixe commun à toutes les stations, ce qui implique de fortes contraintes sur la spécification des sites d'installation. L'utilisation d'antennes fixes manque de souplesse face à la variété des environnements des sites d'installation. Des études sont en cours pour tenter d'améliorer la robustesse des antennes FRPA aux trajets indirects multiples et aux interférences, en les protégeant grâce à des structures de protection mécanique appelées IMLP (Interference + Multipath Local Protection). Ces protections permettent de mieux maîtriser les réflexions. Mais de telles protections présentent l'inconvénient majeur d'être encombrantes, typiquement un diamètre de 5 à 10 mètres et une hauteur de 2 à 3 mètres et d'être chères à cause des matériaux absorbants utilisés. La robustesse des stations de référence aux trajets indirects multiples et aux interférences peut également être assurée grâce à des traitements de filtrage fréquentiel et temporel réalisés au niveau du récepteur. Différentes techniques de filtrage sont généralement utilisées selon la nature de la perturbation. Mais un inconvénient majeur de ces techniques est qu'elles ne sont optimales que dans un domaine restreint d'hypothèses relatives à la nature de la perturbation. Du fait de leur spécialisation, elles nécessitent l'implémentation d'autant d'algorithmes dédiés, qui ne sont pas sans impact sur la qualité des mesures extraites, notamment sur la stabilité des biais de phases et sur la cohérence entre phase de code et phase de porteuse. La multiplication des algorithmes complique d'autant la complexité de la validation des performances de la station de référence. Le traitement spatial adaptatif d'antenne réseau de type CRPA 35 (Controlled Radiation Pattern Antenna) à formation de voie permet de s'adapter automatiquement et sans connaissance a priori de la configuration des sites d'installation. Mais ce type de traitement présente de nombreux inconvénients. D'une part, il implique l'utilisation d'une antenne réseau complexe et la mise en oeuvre des traitements associés. D'autre part, il conduit à des récepteurs eux aussi complexes et surtout sensibles aux défauts de calibrage des voies RF. Enfin, ce type de traitement permet la réjection des interférences, mais pas la réjection des trajets indirects multiples. L'invention a notamment pour but de rendre les stations de référence GNSS robustes aussi bien aux trajets indirects multiples qu'aux interférences liés aux sites d'installation desdites stations, ceci de manière auto-adaptable à l'environnement local de chaque site. Pour cela, elle propose de réaliser sur la voie de réception du signal d'intérêt une soustraction cohérente des sources perturbatrices, ces sources incluant les trajets indirects multiples et les interférences. Des voies de réception appelées par la suite « Voie de Référence des Perturbations » (VRP) sont créées, ces VRP permettant une isolation spatiale des trajets indirects multiples et des interférences. A cet effet, l'invention a pour objet un dispositif pour éliminer les signaux de perturbation reçus par une station GNSS de référence. Le dispositif comporte des moyens pour recevoir un signal d'intérêt émis par un satellite. Il comporte également des moyens pour recevoir les signaux de perturbation, lesdits moyens incluant des moyens pour recevoir lesdits signaux de perturbation isolés du signal d'intérêt. II comporte aussi des moyens pour soustraire les signaux de perturbation au signal d'intérêt, lesdits moyens incluant des moyens pour estimer la fonction de transfert différentielle W entre la voie de réception du signal d'intérêt et la voie de réception des signaux de perturbation, de manière à réaliser une soustraction cohérente desdits signaux. Par exemple, les signaux de perturbation peuvent inclure des interférences et des trajets indirects issus de réflections multiples du signal d'intérêt. Dans un mode de réalisation préférentiel, les moyens pour recevoir le signal d'intérêt peuvent inclure une antenne principale à diagramme de rayonnement à ouverture sensiblement omnidirectionnelle. Dans un mode de réalisation préférentiel, les moyens pour recevoir les signaux de perturbation peuvent inclure une antenne secondaire à diagramme de rayonnement directif à basses élévations. Avantageusement, les moyens pour recevoir les signaux de perturbation isolément du signal d'intérêt peuvent inclure des moyens pour maximiser le gain de l'antenne secondaire dans la direction des signaux de perturbation et des moyens pour minimiser le gain de l'antenne secondaire dans la direction du signal d'intérêt. Par exemple, l'antenne principale et l'antenne secondaire peuvent être les antennes d'une même station L AAS, l'antenne principale pouvant être asservie en boucle fermée de manière à poursuivre le signal d'intérêt, l'antenne secondaire pouvant être asservie en boucle ouverte à l'antenne principale, de manière à orthogonaliser le signal d'intérêt et les signaux de perturbation. Par exemple, la fonction de transfert différentielle W eut être estimée en calculant périodiquement W = R-'P, où R = EIY(k)YT (k)} désigne la matrice de covariance de la voie de réception des signaux de perturbation et P= E{X(k)YT (k)} désigne le vecteur d'intercorrélation entre les deux voies, X(k) et Y(k) désignant des vecteurs associés à des échantillons synchronisés du signal d'intérêt et des signaux de perturbation respectivement. Avantageusement, la soustraction peut être réalisée sur des signaux résultant du désétalement spectral par corrélation avec le code local, après compensation de la différence W desdites fonctions de transfert en calculant ,§'(k) = X(k) -WTY(k) , où S'(k) désigne une estimation d'un échantillon du signal d'intérêt après compensation des perturbations. Par exemple, la compensation de la différence W peut être réalisée au moyen d'un filtre FIR disposé sur la voie de réception des perturbations, les coefficients du filtre FIR étant ajustés périodiquement. Outre la simplicité du traitement de soustraction cohérente, un avantage clé de la présente invention est qu'elle ne nécessite pas de modéliser a priori la nature des sources de perturbation et des signaux interférents. Elle est aussi bien adaptée à la réflexion spéculaire qu'à la réflexion diffuse des trajets indirects multiples, aussi bien adaptée aux interférences bande étroite qu'aux interférences bande large, aussi bien adaptée aux ondes continues qu'aux ondes pulsées. Le même traitement s'applique indifféremment à tous types de perturbation et reste efficace sur tout type d'environnement local des stations de réception : il ne nécessite pas de contrôle a priori de ces environnements. D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit faite en regard de dessins annexés qui 15 représentent : la figure 1, par un schéma, un exemple de réalisation d'une station de référence LAAS selon l'invention; la figure 2, par un diagramme un exemple de combinaison de signaux selon l'invention; 20 - la figure 3, par un diagramme, un exemple de réalisation de l'invention appliquée à deux canaux de poursuite homologues des deux récepteurs d'une station LAAS. 25 Dans un mode de réalisation élémentaire couvrant une majorité des configurations des sites d'installation des stations de références GNSS, une antenne secondaire indépendante de l'antenne principale de réception des signaux d'intérêt émis par les satellites peut avantageusement être utilisée. La directivité de cette antenne secondaire peut permettre de recevoir 30 les signaux incidents à basses élévations, les basses élévations formant le principal secteur de réception des trajets indirects multiples et des interférences des émetteurs sol, de façon isolée des signaux d'intérêt. Le gain de cette antenne secondaire étant plus important sur les sources de perturbation à faible élévation que sur les signaux d'intérêt à forte élévation 35 situés en dehors de ce secteur, une VRP fournit une estimation des signaux perturbateurs qui permet, moyennant une modélisation adaptative des canaux de transmissions, une soustraction cohérente au signal fournit par l'antenne principale. La figure 1 illustre par un schéma un exemple de réalisation de la présente invention dans une station de référence LAAS de type IMLA (Integrated Multipath Limiting Antenna). Cet exemple de réalisation peut comporter une antenne principale HZA (High Zenith Antenna) pour la réception à hautes élévations de signaux d'intérêt X(t) émis par un satellite. Cette antenne HZA a un diagramme de rayonnement directif à ouverture sensiblement omnidirectionnelle comme illustré par la figure 1. Le présent exemple de réalisation peut également comporter une antenne secondaire MLA (Multipath Limiting Antenna) pour la réception de signaux perturbateurs Y«) à basses élévations, ces signaux perturbateurs Y(t) pouvant être émis par exemple par des radars aéroportés ou des réseaux de téléphonie mobile. Cette antenne MLA a également un diagramme de rayonnement directif. Il faut noter que X(t) est perturbé par Y(t), c'est-à-dire que X(t) inclut une composante perturbatrice que la présente invention se propose de compenser. Le traitement spatial de l'antenne MLA, notamment le pilotage de son diagramme de rayonnement, est optimisé de façon à assurer en continu la meilleure indépendance possible entre une VRP formée par l'antenne MLA et des signaux d'intérêt X(t) provenant d'un satellite. II s'agit notamment de maximiser le gain de l'antenne MLA dans la direction des signaux perturbateurs Y(t) et de minimiser le gain de l'antenne MLA dans la direction des signaux X(t). L'antenne MLA permet la formation de plusieurs VRPs, une VRP adaptée à chacun des satellites visibles. Ceci permet d'assurer, par orthogonalisation, la meilleure réjection possible entre les signaux satellites considérés et l'ensemble des sources de perturbation. La figure 2 illustre par un diagramme un exemple de combinaison selon l'invention des signaux X(t) disponibles sur la voie principale avec les signaux Y(t) disponibles sur la VRP. Il s'agit de réaliser une soustraction cohérente des sorties des deux antennes, c'est-à-dire entre la voie de mesure de l'antenne principale HZA et les VRP de l'antenne MLA, ceci sans calibrer au préalable les voies de traitement qui s'étendent du port d'entrée de chacune des deux antennes jusqu'à la sortie des signaux pour soustraction, ces voies de traitement induisant notamment un biais de phase de code et de porteuse. Pour cela, une estimation adaptative de la fonction de transfert différentielle W entre les deux voies de traitement peut être réalisée. Ainsi, si l'on note X(k) et Y(k) les vecteurs associés à des échantillons synchronisés du signal X(t) et du signal Y«) respectivement, 10 alors W peut se calculer de la manière suivante : W = R-'P où R = E{Y(k)YT (k)} désigne la matrice de covariance de la VRP et 15 P = E{X(k)YT (k)} désigne la matrice d'intercorrélation entre la VRP et la voie principale. La soustraction est réalisée après compensation de la différence W desdites fonctions de transfert de la manière suivante : 20 S(k) = X (k) - W T Y(k) où S(k) désigne une estimation de l'échantillon du signal d'intérêt après compensation des perturbations. 25 Cette compensation peut être effectuée par filtrage transverse appliqué sur la VRP grâce à un filtre numérique FIR (Finite Impulse Response). Cette méthode permet de respecter les contraintes de performances imposées aux stations de référence afin de respecter la phase 30 des signaux satellites, notamment en minimisant les variations de Temps de Propagation de Groupe (TPG) de la chaîne de traitement. Le traitement s'appliquant à l'amplitude et à la phase des signaux reçus, il peut être appliqué à cadence réduite, ce qui diminue le coût, après démodulation des signaux et désétalement spectral par corrélation avec le 35 code local, sur les signaux I&Q généralement échantillonnés à 50 Hz. Les signaux issus de la démodulation du code local sont ainsi débarrassés des biais provoqués par les trajets indirects multiples et les interférences, avant estimation du retard de groupe par le discriminateur de code et du retard de phase porteuse par le discriminateur de phase. Dans la station LAAS du présent exemple de réalisation, l'antenne principale HZA et l'antenne secondaire MLA ont des diagrammes de directivité complémentaires en élévation pour la poursuite des signaux satellites. C'est le cas de toutes les stations LAAS. Ceci permet, par couplage des sorties, d'assurer une sensibilité de réception homogène pour toutes les élévations. Cependant, cette complémentarité ne contribue en rien à l'amélioration de la robustesse aux perturbations de l'environnement. La présente invention propose avantageusement de coupler les traitements de poursuite des deux chaînes de réception. En effet, les récepteurs GNSS mettent classiquement en oeuvre des chaînes de poursuite de type boucle de code (DLL : Delay Lock Loop) et boucle de phase (PLL : Phase Lock Loop) en boucle fermée asservie aux écarts d'erreurs. Ceci suppose donc que les deux chaînes de traitement reçoivent les mêmes satellites avec une sensibilité suffisante pour assurer la continuité de chacune des poursuites. Or, la présente invention propose d'utiliser deux antennes afin justement de fournir des couvertures en élévation différentes pour pouvoir apprécier les perturbations sans être gêné par le signal d'intérêt : la présente invention compromet donc la poursuite simultanée des signaux satellites sur les deux récepteurs indépendants, puisque l'un des récepteurs ne dispose pas du signal d'intérêt. Cependant, grâce à l'augmentation du contraste, l'invention permet la poursuite de satellites même situés dans la zone de couverture en élévation de l'antenne secondaire. Ceci n'est pas possible avec les traitements classiques, qui ne permettent pas la séparation entre signal utile et perturbation dans cette zone. L'invention propose d'asservir la consigne de poursuite du récepteur de la VRP à celle du récepteur de la voie principale. Ainsi le récepteur de la VRP fonctionne en boucle ouverte. En l'absence de signal d'intérêt sur la VRP, l'invention permet, après démodulation, de collecter les signaux propres aux seules perturbations. Le récepteur de la voie principale fonctionne pour sa part classiquement, c'est-à-dire en boucle fermée asservie sur la phase de réception du signal d'intérêt venant du satellite. Ce signal, initialement dégradé par les perturbations locales visibles par l'antenne principale, est ensuite débarrassé selon l'invention par soustraction cohérente adaptative des perturbations estimées sur la VRP. La figure 3 illustre par un diagramme un exemple de réalisation de l'invention appliquée à deux canaux de poursuite homologues des deux récepteurs 10 et 20 de la station LAAS poursuivant le même satellite. Les deux récepteurs 10 et 20 sont partiellement synchronisés. Ils sont couplés au travers de la commande de leurs NCOs respectifs 15 et 25 (Numerically- Controlled Oscillator), en phase de code et de porteuse. Seul le récepteur maître 10 de l'antenne principale fonctionne en boucle fermée. Le récepteur esclave 20 de l'antenne secondaire fonctionne en boucle ouverte, en utilisant le même asservissement que celui de la voie principale. Il faut donc idéalement autant de voies couplées que de satellites visibles. La soustraction cohérente entre les deux voies s'effectue après égalisation des différences existantes entre les fonctions de transfert des deux voies, c'est-à- dire des différences entre les centres de phases des antennes, entre les biais de phase de code et de porteuse des canaux analogiques RF. L'égalisation est réalisée de façon adaptative et ne nécessite pas de calibrage préalable des antennes et des chaînes RF. En pratique, elle est réalisée sur la base de l'estimation de la matrice d'inter-corrélation entre les voies d'antenne. A partir de la matrice d'inter-corrélation, on estime la fonction de transfert différentielle des deux voies, qui permet ainsi de compenser la différence de marche des perturbations reçues. Cette compensation est réalisée au moyen d'un filtre FIR 26 dont les coefficients sont calculés en temps réel par le module 30 pour s'adapter à l'évolution des fonctions de transfert des deux chaînes de réception, selon la direction d'arrivée des satellites. Le filtre FIR est appliqué uniquement sur la voie de référence de perturbation, afin de ne pas perturber les mesures du signal satellite réalisées sur la chaîne principale. L'estimation des coefficients du filtre FIR est faite de façon adaptative par un module 30 à partir des signaux I&Q prélevés après corrélation par le code local, qui conserve l'information de phase du signal reçu et qui permet d'avoir un rapport signal sur bruit significatif. Les erreurs dues aux perturbations évoluant doucement dans le temps, essentiellement sous l'effet des évolutions des fonctions de transferts des antennes et des trajets multiples selon la direction des satellites, il est possible de rafraîchir les coefficients à un rythme inférieur à 1 Hz, et d'utiliser des signaux démodulés en bande de base et après désétalement spectral par corrélation avec le code local. Un filtre FIR à moins de 20 coefficients est suffisant pour décrire la fonction de transfert différentielle avec précision. Cependant, le nombre de coefficients pourra être dimensionné en fonction des performances visées. Dans un mode de réalisation amélioré, les VRP peuvent être adaptées à la configuration des sites d'installation des stations de référence, notamment aux sources connues de réflexion du signal et aux sources d'interférence. Il faut bien comprendre que le mode de réalisation à partir de deux antennes séparées décrit précédemment n'est pas limitatif, d'autres modes de réalisation pouvant être envisagés, notamment à partir d'une antenne réseau à agilité de faisceau. En effet, l'antenne MLA d'une station LAAS peut être constituée d'un ensemble d'antennes mises en réseau avec une apodisation fixe, c'est-à-dire un coefficient de pondération par antenne dans le réseau, cette apodisation permettant de garantir le diagramme souhaité. L'invention propose de rendre cette apodisation adaptative, afin de créer un « trou » de réception dans la direction de chacun des satellites GNSS d'intérêt, comme illustré par la figure 1 qui représente un trou dans le diagramme de l'antenne secondaire, ce trou étant dans la direction d'un satellite GPS émettant un signal d'intérêt représenté comme une sinusoïde. Pour cela, l'invention propose de créer un « zéro » de la fonction de directivité de l'antenne secondaire dans la direction de chacun des satellites. Cette apodisation adaptative peut également être mise à profit pour modifier le diagramme de l'antenne MLA afin de tenir compte de la spécificité du site d'installation de la station GNSS de référence, afin d'améliorer les performances du système. En effet, il est possible d'augmenter le gain dans la ou les directions d'arrivée des interférences, afin de mieux les estimer et donc mieux les annuler par soustraction. Ceci nécessite le calcul d'un jeu de coefficients pour chacun des satellites poursuivis, chaque jeu contenant autant de coefficients que d'antennes considérées dans le réseau MLA. Le récepteur d'une station GNSS de référence étant fixe, cela permet un rafraîchissement des coefficients lent puisque la vitesse angulaire de défilement des satellites est également faible. Un avantage de la présente invention est qu'elle permet directement une détection des interférences et donc facilite la surveillance du site de la station de référence. Un autre avantage de la présente invention est qu'elle peut aussi permettre d'optimiser le diagramme de réception de l'antenne secondaire afin de prendre en compte les contraintes spécifiques à chaque site d'installation des stations de référence, comme la direction des trajets indirects multiples et des interférences, en regard des signaux d'intérêt. Un autre avantage de la présente invention est que ses principes peuvent être étendus à tout dispositif muni de plusieurs voies VRP adaptées ou non à la configuration de chaque site de référence dont l'environnement a été caractérisé en termes de directions des sources d'interférences et/ou des sources de réflexion de trajets indirects multiples. Un autre avantage de la présente invention est que la VRP peut être utilisée pour assurer une surveillance indépendante des sites d'installation, en vue de prévenir l'apparition d'interférences accidentelles ou volontaires, ou encore de sources de réflexion imprévues. Un autre avantage de la présente invention est qu'elle peut être mise en oeuvre à faible coût dans une station LASS, en conservant 30 l'architecture existante de ses récepteurs. REVENDICATIONS1. Dispositif pour éliminer les signaux de perturbation reçus par une station GNSS de référence, le dispositif étant caractérisé en ce qu'il comporte : des moyens pour recevoir un signal d'intérêt émis par un satellite ; des moyens pour recevoir les signaux de perturbation, lesdits moyens incluant des moyens pour recevoir lesdits signaux de perturbation isolés du signal d'intérêt ; des moyens pour soustraire les signaux de perturbation au signal d'intérêt, lesdits moyens incluant des moyens pour estimer la fonction de transfert différentielle W entre la voie de réception du signal d'intérêt et la voie de réception des signaux de perturbation, de manière à réaliser une soustraction cohérente desdits signaux. 2. Dispositif selon la revendication 1, caractérisé en ce que les signaux de perturbation incluent des interférences et des trajets indirects issus de 15 réflections multiples du signal d'intérêt. 3. Dispositif selon la revendication 1, caractérisé en ce que les moyens pour recevoir le signal d'intérêt incluent une antenne principale à diagramme de rayonnement à ouverture omnidirectionnelle. 4. Dispositif selon la revendication 1, caractérisé en ce que les moyens pour recevoir les signaux de perturbation incluent une antenne secondaire à diagramme de rayonnement directif à basses élévations. 25 5. Dispositif selon la revendication 4, caractérisé en ce que les moyens pour recevoir lesdits signaux de perturbation isolément du signal d'intérêt incluent: des moyens pour maximiser le gain de l'antenne secondaire dans la direction des signaux de perturbation, et; 30 des moyens pour minimiser le gain de l'antenne secondaire dans la direction du signal d'intérêt. 6. Dispositif selon la revendication 5, caractérisé en ce que l'antenne principale et l'antenne secondaire sont les antennes d'une même station 20LAAS, l'antenne principale étant asservie en boucle fermée de manière à poursuivre le signal d'intérêt, l'antenne secondaire étant asservie en boucle ouverte à l'antenne principale, de manière à orthogonaliser le signal d'intérêt et les signaux de perturbation. 7. Dispositif selon la revendication 1, caractérisé en ce que la fonction de transfert différentielle W est estimée en calculant périodiquement W =R-'P , où R = *(k)YT (k)} désigne la matrice de covariance de la voie de réception des signaux de perturbation et P = E{X(k)YT (k)} désigne le vecteur d'intercorrélation entre les deux voies, X(k) et Y(k) désignant des vecteurs associés à des échantillons synchronisés du signal d'intérêt et des signaux de perturbation respectivement. 8. Dispositif selon la revendication 7, caractérisé en ce que la soustraction est réalisée sur des signaux résultant du désétalement spectral par corrélation avec le code local, après compensation de la différence W desdites fonctions de transfert en calculant S(k) = X(k)-WTY(k), où S'(k) désigne une estimation d'un échantillon du signal d'intérêt après compensation des perturbations. 9. Dispositif selon la revendication 8, caractérisé en ce que la compensation de la différence W est réalisée au moyen d'un filtre FIR disposé sur la voie de réception des perturbations, les coefficients du filtre FIR étant ajustés périodiquement.25