FR 2932277 A1 20091211 FR 0803169 A 20080606 PROCEDE DE PROTECTION D'UN UTILISATEUR DE RECEPTEUR DE RADIONAVIGATION VIS-A-VIS DE MESURES DE PSEUDO-DISTANCES ABERRANTES La présente invention se rapporte à un procédé de protection d'un utilisateur de récepteur de radio-navigation vis-à-vis de mesures de pseudo-distances aberrantes. L'erreur sur la position géographique indiquée par un récepteur de radio- navigation par satellite dépend des erreurs qui entachent les mesures de pseudodistances déterminées par l'algorithme qu'utilise ce récepteur et des erreurs contenues dans le message de navigation diffusé par le satellite. Afin de limiter ces erreurs et d'assurer la protection des utilisateurs, il faut disposer de moyens permettant d'identifier de telles erreurs pour les éliminer et de calculer une borne sur l'erreur de position en fonction des mesures disponibles, en tenant compte de spécifications d'intégrité et de continuité inhérentes au contexte d'utilisation (atterrissage d'un aéronef, etc). Ces moyens constituent la base d'un dispositif disposant d'une fonction RAIM (Receiver Autonomous Integrity Monitoring). Les récepteurs GNSS d'aviation civile actuels ne peuvent être utilisés sans une fonction RAIM. Les équipements actuels à fonction RAIM souffrent de deux problèmes : - ils sont complètement intégrés, ce qui signifie qu'il n'est pas possible de choisir séparément l'équipement qui acquiert le signal de navigation et celui qui calcule la position de l'équipement qui assure les fonctions d'intégrité, - ils sont basés sur des algorithmes de type moindre carrés qui sont déstabilisés par la présence de mesures erronées quelle que soit l'amplitude (même infinitésimale) de l'erreur qui affecte ces mesures, ce qui pose des problèmes de fiabilité des solutions de position proposées. La présente invention a pour objet un procédé de protection d'un utilisateur de récepteur de radionavigation vis-à-vis de mesures de pseudo-distances aberrantes, ce procédé pouvant être mis en oeuvre par un équipement indépendant de celui qui acquiert le signal de navigation et qui calcule la position de l'équipement assurant les fonctions d'intégrité, et qui ne soit pas déstabilisé par la présence de mesures erronées, quelle que soit l'amplitude des erreurs affectant ces mesures. Le procédé de protection conforme à l'invention est caractérisé, dans le cas le plus général, en ce que l'on pondère les N résidus de l'estimation du vecteur d'état lors de l'estimation de la solution de navigation, la pondération étant réalisée d'après la statistique estimée des erreurs de ces mesures, cette pondération pouvant être totale ou partielle. Si elle est totale, la prise en compte des résidus est totale, avec un facteur de pondération égal à 1 ou égal à 0 pour l'exclusion, et si elle est partielle, le facteur de pondération est compris entre 0 et 1. Le cas de la pondération totale couvre les méthodes robustes � RAIM-MSTD � et � RAIM-LTS �, et le cas de la pondération partielle couvre les méthodes robustes de type � RAIM-M-estimateurs �. Selon une autre caractéristique de l'invention, on sélectionne le sous-ensemble de N-h résidus de l'estimation du vecteur d'état de moindre variance, que l'on utilise cette variance, majorée d'un coefficient d'adaptation, ainsi que la moyenne du même sous-ensemble, pour centrer puis normaliser les N entrées, qui, une fois mises au carré, sont comparées à un seuil statistique et rejetées si elles sont supérieures à ce seuil, avec h=N-1 et/ou N-2, N étant le nombre de résidus. Le but de la comparaison au seuil statistique est de ne pas avoir à isoler systématiquement les h résidus de l'estimation du vecteur d'état. C'est la méthode � RAIM-MSTD �. Selon une autre caractéristique de l'invention, on sélectionne le sous-ensemble de N-h résidus de l'estimation du vecteur d'état tel que la somme de ces mesures au carré est minimale, on estime la moyenne et la variance de ce sous-ensemble, et on utilise cette moyenne et cette variance (cette dernière étant avantageusement multipliée par un coefficient d'adaptation supérieur à 1) pour centrer puis normaliser les N entrées, qui, une fois mises au carré, sont comparées à un seuil statistique et rejetées si elles sont supérieures à ce seuil. C'est la méthode � RAIM-LTS �. Selon une autre caractéristique de l'invention, les N résidus de l'estimation du vecteur d'état sont pondérés par une méthode itérative qui consiste, à chaque itération, à estimer la moyenne et la variance des résidus de l'estimation du vecteur d'état, à calculer des poids à l'aide d'une fonction de pondération dont les entrées sont les résidus centrés et normalisés par la variance, et à multiplier les résidus par ces poids, le processus itératif s'arrêtant lorsque la somme des carrés des différences entre deux résidus consécutifs est inférieure à un seuil, la moyenne et la variance (cette dernière majorée d'un coefficient d'adaptation) issues de ce processus étant utilisées pour centrer puis normaliser les N entrées, qui, une fois mises au carré, sont comparées à un seuil statistique et rejetées si elles sont supérieures à ce seuil. C'est la méthode � RAIM-M-estimateurs �. De façon détaillée, le procédé de protection conforme à l'invention est caractérisé en ce qu'il comporte les étapes suivantes de calcul de l'information d'intégrité par calcul des résidus de mesures de pseudo-distances obtenues à partir de la position géographique et du décalage d'horloge fournis par un récepteur primaire de radionavigation et à partir des mesures utilisées par le récepteur primaire pour obtenir ces valeurs : - Formation de tous les sous-ensembles de résidus de cardinal h=N-1 et/ou N-2, N étant le nombre de résidus, Calcul de l'écart-type de chaque sous-ensemble et détermination, 6n,in étant le plus petit écart-type et utilisé pour pondérer les résidus, hmin la valeur de cardinal pour laquelle amin a été obtenu, Ymi^ le vecteur de 20 cardinal hmin correspondant, et mmin la moyenne de Ymin Calcul d'un premier vecteur défini par : rv�. I = \/(Ymm'° )2 6m in les vecteurs rvecl sont triés par ordre croissant pour obtenir : - calcul du facteur f défini par : f = (x2) ' h m " ,1\ , où (, 2)�.,I) est l'inverse de la distribution du x2 à un degré de liberté, 25 Emin est alors re-pondéré de la façon suivante : trie rvec,l hmin) * trie ,,I f Emin calcul d'un nouveau vecteur de résidus : r e,,2 = 1).(7 ù mmin )2 6111111, 2 définition d'un seuil T tel que : T = (x 2) 1(P,1) , où P est une probabilité dont la valeur optimale doit être déterminée par simulation probabiliste, - Chaque élément de rvec,2 est comparé à T. Si rVec,, (i) > T , on déclare qu'une panne est détectée, et le satellite i est exclu. Selon une autre caractéristique de l'invention, on calcule une valeur d'horloge corrigée (11) lorsqu'une erreur dans les mesures d'entrée est détectée. La présente invention sera mieux comprise à la lecture de la description détaillée d'un mode de réalisation, pris à titre d'exemple non limitatif et illustré par le dessin annexé, sur lequel : - la figure unique est un diagramme simplifié des étapes successives de mise en oeuvre du procédé de l'invention. Sur le diagramme du dessin, on a symbolisé en 1 un récepteur primaire classique de radionavigation par satellite fournissant d'une part, sur une voie 2, la solution estimée par le récepteur primaire (position géographique et décalage d'horloge), et sur une autre voie 3 les mesures de pseudo-distances utilisées par le récepteur primaire pour obtenir la solution. Si les mesures transmises par le récepteur primaire ne sont pas prétraitées, il convient de leur faire subir un pré-traitement, connu en soi, les affranchissant des erreurs de propagation et de mesures, comme symbolisé par le rectangle 4 en traits interrompus. Le calcul de l'information d'intégrité comporte ensuite les étapes suivantes, à partir du calcul des résidus (5) de mesures de pseudo-distances (notés Yi, où i est l'indice du satellite, et N est le nombre de ces résidus) : 1. Formation de tous les sous-ensembles de résidus de cardinal h=N-1 et/ou N-2, 2. Calcul de l'écart-type de chaque sous-ensemble (6). Soit 6,,,;,, le plus petit écart-type. Celui-ci est pris comme référence (7). On définit ensuite : a. hm;,, la valeur de cardinal pour laquelle ami,, a été obtenu b. Ymir, le vecteur de cardinal hm;,, correspondant,30 . a,äir, est utilisé pour pondérer les résidus. On calcule un premier vecteur tel que : ~ (Y ù menin l / rive, = Emin c. mr,,ir, la moyenne de Ymir, ( rvec,l 4. Pour les différents satellites, les premiers vecteurs rvee,l sont triés par ordre croissant pour obtenir : On calcule par ailleurs le facteur f défini par : 5 fû(x2)-1 hmin N liberté., où (x2) ' (.,1) est l'inverse de la distribution du x2 à un degré de 5. ar,,ir, est alors repondéré de la façon suivante : 6m11,,2 = Ine ,_ min J 6. On calcule un nouveau vecteur de résidus : rvc ,c,2 = " (Y û menin )2 7. Le seuil T est défini par : T = (x2) �P,1), où P est une probabilité dont la valeur optimale doit être déterminée par simulation par une technique probabiliste de type � Monte-Carlo �. 8. Chaque élément de rVec,2 est comparé à T. Si rvec.2 (i) > T , on déclare qu'une panne est détectée, et le satellite i est exclu (8). 9. Une fois la détection des mesures éventuellement fautives effectuée, le calcul des rayons de protection (10) peut se faire de la façon habituelle compte tenu du nombre de mesures retenues. On notera que P, qui représente la probabilité d'exclure un satellite non- fautif, n'est cependant pas directement assimilable à la probabilité d'exclusion spécifiée par l'OACI. D'autre part, le test au point 8 est effectué indépendamment pour chaque satellite : il n'y a donc pas de détection sans exclusion. Les sorties de ce dispositif sont donc : 6min,2 - Une solution de position éventuellement corrigée (11), si une erreur dans les mesures d'entrée est détectée. Une valeur de rayon de protection (12) qui permet de garantir la solution de position du dispositif, incluant éventuellement une alarme. Le procédé de l'invention permet, de façon autonome de tout segment sol (donc implémentant une fonctionnalité RAIM) : - d'augmenter les performances d'un récepteur (dit � primaire �) du commerce sans fonction d'intégrité, - de détecter d'éventuelles erreurs contaminant des mesures en entrée du calcul de la position, par l'utilisation d'un algorithme statistique d'estimation robuste, c'est-à-dire réellement non affecté par les erreurs de mesures, et mettant en oeuvre un critère dynamique, - de calculer une correction robuste pour la position fournie par le récepteur 15 primaire en excluant une telle erreur si elle est détectée, - de calculer des bornes sur l'erreur de position en fonction des mesures disponibles, et en tenant compte de spécifications d'intégrité et de continuité inhérentes au contexte d'utilisation (par exemple en phase d'atterrissage d'un aéronef). Ces bornes sont des seuils qui ne doivent pas être dépassés plus d'une fois pour N 20 mesures de distances, avec N=105 à 10', par exemple, - des performances supérieures à celles de l'algorithme RAIM standard (utilisant la méthode des moindres carrés). Du fait qu'il n'existe pas actuellement de procédé permettant d'ajouter une fonction d'intégrité à un récepteur ne la fournissant pas d'origine, les dispositifs de 25 l'art antérieur, pour obtenir une performance équivalente d'intégrité/disponibilité, ne peuvent pas utiliser les algorithmes courants. Pour améliorer les performances, l'autre voie possible serait d'utiliser des traitements séquentiels ou filtrants, mais comme il est exclu de procéder à de tels traitements pour l'estimation du point (temps d'alarme d'intégrité trop long, qui est par exemple de l'ordre de 6 secondes 30 pour des applications civiles, mais peut être de l'ordre de plusieurs centaines de secondes avec un filtrage efficace), il n'y a pas d'autre solution totalement autonome pour assurer l'intégrité à ce niveau de performance. - Le procédé de l'invention permet d'optimiser le choix du matériel d'acquisition du signal de navigation indépendamment de celui de traitement RAIM. - Le procédé RAIM rendu robuste selon l'invention permet de fiabiliser l'estimation du point (position géographique), en améliorant les performances d'intégrité (détection et exclusion) par rapport à un RAIM standard. - Ce procédé apporte une capacité de détection et d'exclusion de mesures GNSS fautives optimisée par rapport au RAIM standard. En particulier, il permet d'améliorer le taux de disponibilité d'un système GNSS compatible avec les exigences aviation civile. En résumé, le procédé de l'invention offre des performances bien meilleures que celles des algorithmes standard tant en termes d'intégrité (détection) que de disponibilité (fausses alarmes, optimisation des rayons de protection) et peut être appliqué à un récepteur de radio-navigation ne possédant pas d'origine de protection RAIM rendue robuste contre les erreurs de mesure.20 REVENDICATIONS1 Procédé de protection d'un utilisateur de récepteur de radionavigation vis-à-vis de mesures de pseudo-distances aberrantes, caractérisé en ce que l'on pondère les N résidus de l'estimation du vecteur d'état lors de l'estimation de la solution de navigation, la pondération étant réalisée d'après la statistique estimée des erreurs de ces mesures, 2. Procédé selon la revendication 1, caractérisé en ce que la pondération est totale, avec un facteur de pondération égal à 1 ou O. 3. Procédé selon la revendication 1, caractérisé en ce que la pondération est partielle avec un facteur de pondération compris entre 0 et 1. 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce que l'on sélectionne le sous-ensemble de N-h résidus de l'estimation du vecteur d'état de moindre variance, que l'on utilise cette variance, majorée d'un coefficient d'adaptation, ainsi que la moyenne du même sous-ensemble, pour centrer puis normaliser les N entrées, qui, une fois mises au carré, sont comparées à un seuil statistique et rejetées si elles sont supérieures à ce seuil, avec h=N-1 et/ou N-2, N étant le nombre de résidus. 5. Procédé selon la revendication 1, caractérisé en ce que 1' on sélectionne le sous-ensemble de N-h résidus de l'estimation du vecteur d'état tel que la somme de ces mesures au carré est minimale, que l'on estime la moyenne et la variance de ce sous-ensemble, que l'on utilise cette moyenne et cette variance pour centrer puis normaliser les N entrées, qui, une fois mises au carré, sont comparées à un seuil statistique et rejetées si elles sont supérieures à ce seuil. 6. Procédé selon la revendication 5, caractérisé en ce que la variance est multipliée par un coefficient d'adaptation supérieur à 1. 7. Procédé selon la revendication 1, caractérisé en ce que les N résidus de l'estimation du vecteur d'état sont pondérés par une méthode itérative qui consiste, à chaque itération, à estimer la moyenne et la variance des résidus de l'estimation du vecteur d'état, à calculer des poids à l'aide d'une fonction depondération dont les entrées sont les résidus centrés et normalisés par la variance, et à multiplier les résidus par ces poids, le processus itératif s'arrêtant lorsque la somme des carrés des différences entre deux résidus consécutifs est inférieure à un seuil, la moyenne et la variance issues de ce processus étant utilisées pour centrer puis normaliser les N entrées, qui, une fois mises au carré, sont comparées à un seuil statistique et rejetées si supérieures. 8. Procédé selon la revendication 7, caractérisé en ce que la variance est multipliée par un coefficient d'adaptation supérieur à 1. 9. Procédé selon la revendication 1, caractérisé en ce qu'il comporte les étapes suivantes de calcul de l'information d'intégrité par calcul des résidus (5) de mesures de pseudo-distances obtenues à partir de la position géographique et du décalage d'horloge (2) fournis par un récepteur primaire de radionavigation et à partir des mesures (3) utilisées par le récepteur primaire pour obtenir ces valeurs : - Formation de tous les sous-ensembles de résidus de cardinal h=N-1 et/ou N-2, N étant le nombre de résidus, - Calcul de l'écart-type de chaque sous-ensemble (6) et détermination, Emin étant le plus petit écart-type et utilisé comme référence (7) pour pondérer les résidus, hmin la valeur de cardinal pour laquelle Emin a été obtenu, Ymin le vecteur de cardinal hmin correspondant, et mmin la moyenne de Ymin , - Calcul d'un premier vecteur défini par : rrec,, = '\/(Y û min ) °min - les vecteurs rvec,l sont triés par ordre croissant pour obtenir : r`r1e 2 )ùI - calcul du facteur f défini par : f = 1 , où 2 (.,1) est N l'inverse de la distribution du x2 à un degré de liberté, 25 - Emin est alors re-pondéré de la façon suivante : °min,2 mie h rvec,l min ) Emin f- calcul d'un nouveau vecteur de résidus : rvec 2 (Y ù m,nin ermin,2 définition d'un seuil T tel que : T = (x2) ' (P,l), où p est une probabilité dont la valeur optimale doit être déterminée par simulation probabiliste, - Chaque élément de rvec,2 est comparé à T. Si rVec 2 (i) > T , on déclare qu'une panne est détectée (8), et le satellite i est exclu. . 10. Procédé selon la revendication 9, caractérisé en ce que la simulation probabiliste utilisée pour la définition du seuil T est une méthode du type � Monte-Carlo �. 11. Procédé selon la revendication 9 ou 10, caractérisé en ce que lorsque les mesures 10 transmises par le récepteur primaire ne sont pas prétraitées, on leur fait subir un pré-traitement les affranchissant des erreurs de propagation et de mesures (4). 12. Procédé selon l'une des revendications 9 à 11, caractérisé en ce que l'on calcule une position corrigée (9) lorsqu'une erreur dans les mesures d'entrée est détectée. 13. Procédé selon l'une des revendications 9 à 12, caractérisé en ce que l'on calcule 15 une valeur d'horloge corrigée (11) lorsqu'une erreur dans les mesures d'entrée est détectée. 14. Procédé selon l'une des revendications 9 à 13 caractérisé en ce que l'on calcule des rayons de protection (10) relatifs aux mesures de façon habituelle. 15. Procédé selon l'une des revendications 9 à 14, caractérisé en ce que l'on émet 20 une alarme (12) lorsqu'une erreur dans les mesures d'entrée est détectée. 16. Procédé selon l'une des revendications 9 à 15, caractérisé en ce qu'il est appliqué à un récepteur de radio-navigation ne possédant pas d'origine de protection robuste contre les erreurs de mesure. 25