Brevet FR2737014A1 — publié le 19970124

SYSTEME ET PROCEDE DE MESURE DU CONTENU ELECTRONIQUE TOTAL DE LA IONOSPHERE DESCRIPTION Domaine technique La présente invention concerne un système et un procédé de mesure du contenu électronique total de la ionosphère. Etat de la technique antérieur Comme décrit dans la référence [l], la ionosphère est une série de couches ionisées excitées par les rayonnements solaires. Elle se caractérise par sa densité électronique N (nombre d'électrons par mètre cube) dépendant de l'altitude h et son contenu électronique total CET (nombre total d'électrons dans une colonne verticale de 1 m2 de section), qui peut être défini par les relations suivantes 10-16 e./m2 La ionosphère est un milieu dispersif pour les ondes électromagnétiques. Sa fonction de transfert Hiono(f) est telle que H(iono)(f) = A(iono)(f).e-j.#icno(f) Aiono(f) est le coefficient de transmission de la ionosphère. Ce coefficient dépend de la fréquence.Les transmissions ne sont plus très efficaces pour des fréquences inférieures à 100 MHz, en raison de réflexions ou rebonds des ondes électromagnétiques sur les couches ionisées. Le déphase #iono(f) de Hiono(f) , en première approximation, pour un trajet terre-espace est tel que #iono(f) = -2.#.40,3.CET 1 c.sin(Ei) f Ei est l'angle d'élévation de l'onde électromagnétique, à l'altitude du maximum de N (point subionosphérique). Le temps de propagation de groupe Igiono(f) de la ionosphère est donc #giono(f) 4o,3.cET . c.sin(Ei)T2 L'erreur de ADiono de mesure de pseudodistance se déduit de Igiono(f) : tDiono = c . Igiono(f) Le temps de propagation de phase Tcpiono(f) de la ionosphère est donc -40,3.CET 1 IÇiono(f) # c.sin (Ei) f L'erreur AViono de mesure de pseudovitesse radiale se déduit du temps de propagation de phase de la ionosphère, noté T(piono(f) AViono = d##iono(f) d#giono(f)# - 40,3.CET.dEi/dt - dt dt 40,3.dCET/dt c.sin(Ei).f2 Les systèmes de radio navigation localisation prétendant accéder à une bonne précision de localisation sont jusqu'à présent bifréquences. Supposons que les fréquences utilisées soient fl et f2, les mesures de distances correspondantes sont notées Dl et D2. On obtient un système de deux équations à deux inconnues Dl = d + ADiono(fI) D2 = d + ADiono(f2) d'où : CET = (D1-D2).sin(Ei).f12.f22 40,3.(f12 -f22) d = f2 . Dl - fi . D2 f22 - f12 Des équations similaires sont obtenues dans le cas de mesures de vitesse bifréquence. Parfois, certains systèmes de radio navigation localisation sont contraints par l'utilisation d'une seule fréquence, alors qu'une certaine précision de localisation est recherchée. Le système doit donc disposer d'une ionographie, afin de pouvoir corriger les mesures ; ce qu'il ne peut pas faire directement en général. Cette ionographie peut être effectuée à l'aide de toutes les mesures réalisées par le système lui-même. Cette opération est facilitée si le système dispose de références aux coordonnées parfaitement connues (citons le cas des balises d'orbitographie du système DORIS, qui pourtant est un système bifréquence). En ce cas, les fonctions de ionographie et d'orbitographie sont souvent couplées. Ainsi le système GPS, bien que bifréquence, réalise une ionographie qui est ensuite diffusée dans le message de référence des satellites, sous la forme de coefficients d'un modèle de ionosphère ajusté par des mesures de CET. Cette ionographie est alors utilisée par les récepteurs GPS mono fréquence placés au voisinage immédiat de la surface terrestre. Cet exemple montre qu'un système de radiolocation monofréquence peut également s'appuyer sur une ionographie réalisée par un système mondial bifréquence extérieur, tel que DORIS, PRARE, GPS ou GLONASS. Le modèle de ionosphère diffusé par la constellation GPS a été optimisé pour des trajets terre-espace. Il n'est donc pas toujours adapté au cas d'un récepteur GPS monofréquence placé en orbite basse. L'utilisation du filtrage de Kalman des mesures brutes permet de réaliser simultanément une orbitographie autonome ainsi qu'une ionographie autonome. De tels filtres de Kalman pour les récepteurs GPS C/A spatiaux sont actuellement à l'étude. Le vecteur d'état de ces filtres contient à la fois les paramètres d'orbite, les paramètres d'horloge ainsi que des paramètres de la ionosphérique. La référence [2] montre comment des mesures à une seule fréquence peuvent être utilisées pour réaliser une estimation en temps réel des retards de propagation dus à la ionosphère. La référence [3] décrit une estimation du retard ionosphérique absolu uniquement grâce à des mesures GPS à une fréquence. Le délai ionosphérique est obtenu en mesurant la divergence relative entre le code GPS et la porteuse GPS. La présente invention a pour objet de pallier les inconvénients des dispositifs de l'art antérieur, et donc de mesurer le retard ionosphérique en utilisant les émetteurs VHF en spectre étalé qui existent dans certains satellites, et notamment les petits satellites, permettant une liaison avec le sol. Exposé de l'invention La présente invention concerne un système de mesure du contenu électrique total de la ionosphère comprenant - à bord d'un satellite, un émetteur d'un signal cohérent en spectre étalé transmis en VHK - au sol, un dispositif de réception et de traitement du signal VEIF reçu permettant de mesurer le contenu électronique total de la ionosphère. Le segment bord peut être couplé avec une fonction d'initialisation autonome des liaisons spatiales directives ou/et une fonction de transmission de secours, prévue pour certains satellites. La présente invention concerne également un procédé de mesure du contenu électronique total de la ionosphère, à l'aide d'un signal en spectre étalé transmis à travers la ionosphère dans la bande VHF, ce procédé étant caractérisé en ce qu'on émet à partir d'un satellite un signal constitué d'un code pseudoaléatoire dans la bande VHF. On peut également émettre à partir d'une station sol un signal semblable au précédent, répété par un satellite dans la bande de VHF. Dans les deux cas, on reçoit le signal transmis par le satellite dans la station sol, les mesures disponibles dans la station (vitesse de groupe, vitesse de phase) étant utilisées pour déterminer le contenu électronique total de la ionosphère. Les mesures disponibles dans la station (distance, Doppler) sont corrigées des effets de la ionosphère à l'aide de la connaissance du contenu électronique total. Ces mesures corrigées sont utilisées pour déterminer l'orbite dudit satellite. Brève description des dessins - La figure l illustre le système de l'invention - la figure 2 illustre un exemple de segment bord classique - la figure 3 illustre une courbe d'orbitographie sur laquelle sont représentés les cercles de visibilité de stations d'orbitographie et les cercles de visibilité de zones de service. Exposé détaillé de modes de réalisation Comme représenté sur la figure l le segment bord 10 à l'intérieur d'une satellite 9 comprend un générateur il de signaux en spectre étalé en fréquence intermédiaire (par exemple, un générateur de code C/A GPS ou GLONASS en bande Ll) relié à une antenne 12 au travers d'un mélangeur 13, un oscillateur local OL 14 étant relié au générateur 11 et au mélangeur 13 dans le but d'assurer la cohérence entre le code et la porteuse du signal transmis. L'antenne 12 émet un signal en VHF qui, après traversée de la ionosphère 15, est reçu par l'antenne 16 d'une station sol 17.Cette station sol 17 comprend un récepteur de signaux en spectre étalé (par exemple GPS-C/A ou GLONASS-C/A 18, de bande L1) connecté à l'accès de l'antenne 16 au travers d'un amplificateur et d'un mélangeur 19, et en sortie à un ordinateur 20 lui-même relié à un réseau R de stations sols. Un oscillateur local OL 21 est relié au récepteur 18 et au mélangeur 19. le segment bord 10 et la station sol 17 peuvent avoir des fréquences OL ayant la même valeur. La figure 2 illustre un exemple classique de réalisation d'un segment bord 10. L'oscillateur 30 est relié à un générateur de code pseudo-aléatoire 31, lui-même relié à un sommateur 32, qui reçoit sur une autre entrée un signal provenant d'un générateur de message binaire 33 dont le signal d'entrée est référencé E. L'oscillateur 30 est également relié à un générateur de signal VHF 34. Un modulateur amplificateur 35 (par exemple modulateur MDP2 modulateur de phase à deux états) reçoit les sorties du sommateur 32 et du générateur 34. Un amplificateur 36 reçoit le signal de sortie de ce modulateur 35 et est relié à une antenne 37. Eventuellement un formateur de données 38 peut être positionné sur l'entrée du générateur 33. Un système de satellite de l'art connu permet de desservir plusieurs zones de service industrialisées contraignantes en terme d'environnement de brouillage radioélectrique, comme représenté sur la figure 3, les cercle de visibilité des stations d'orbitographie étant référencés 40 et les cercles de visibilité des zones de service étant référencés 41. La charge utile de transmission en VHF n'est pas sensée être utilisée en tant que telle en dehors du survol de ces zones de service. Un signal de localisation en spectre étalé bénéficiant du maximum de puissance disponible à bord est émis vers le sol lorsque le satellite n'est pas en visibilité d'une zone de service. Ce signal n'est donc pas perturbé par les signaux traités lors des survols des zones de service par la charge utile de transmission. Une ou plusieurs stations réduites de réception situées en dehors de ces zones de visibilité peuvent enregistrer les mesures de pseudovitesse et de pseudodistance réalisées à l'aide du signal reçu, pendant toute la durée du passage de chaque satellite. Ces stations retransmettent leurs mesures à un centre d'orbitographie, par l'intermédiaire d'une ligne de communication (filaire, spatiale ou mixte). Les sites des stations d'orbitographie sont choisis pour être peu brouillés. De plus, le signal de localisation utilise une liaison espace-terre, peu affectée par les brouilleurs. Le signal émis en dehors des périodes de visibilité des zones de service est choisi pour permettre la réalisation d'une orbitographie suffisamment précise pour représenter un poste minime dans le bilan d'erreur de localisation des terminaux utilisateurs situés dans les zones de service. Pour cela, on détermine précisément le retard et le Doppler ionosphérique altérant le signal reçu par les stations d'écoute déportées (ionographie associée aux stations d'orbitographie). Ce retard permet de corriger les mesures utilisées pour l'orbitographie du satellite. Lorsque cette orbitographie suffisamment précise est connue, il est possible de réaliser la ionographie associée aux zones de service à l'aide de terminaux de référence situés dans lesdites zones. En effet, si la position de ces terminaux de référence et des satellites est bien connue à tout instant, les mesures de type station-satelliteterminal-satellite-station permettent de déterminer la ionographie associée aux zones de service, tout comme les mesures de type station-satellite-station. La charge utile génère un pilote en spectre étalé en VHF. Ce signal est généré de façon cohérente à partir de l'oscillateur local (OL) bord. Il bénéficie de la chaîne d'amplification de la voie émission VHF, grâce à un commutateur. Afin d'écrire l'équation de mesure du contenu électronique total (CET) par le récepteur de la station d'orbitographie, on définit les notations suivantes fo = fréquence théorique de la porteuse du signal en spectre étalé (Hz) Ei = angle d'élévation ionosphérique (rad) c = vitesse de la lumière (m/s) Cr = glissement code/porteuse (m/s) Vgr = vitesse de groupe ou vitesse du code (m/s) Vçr = vitesse de phase ou vitesse de la porteuse (m/s) ti = instant de mesure(s) mg = mesure du temps de propagation de groupe, i.e. pseudodistance, dans le cas présent (m) ç = mesure de la phase de la porteuse du signal reçu (rad) af = écart type de la mesure de fréquence reçue (Hz) amg = écart type de la mesure du temps de propagation de groupe (m) On obtient Cr = Vgr - Vçr V mg(t2)-mg(tl) t2 - tl 1 c x) (t2) -#(t1) V#r = - x - x 2r = 2it fo t2-tl L'originalité du système de l'invention consiste à utiliser la bande de fréquence VHF bien plus adaptée que d'autres bandes (comme la bande L) pour réaliser les mesures de CET à l'aide d'un seule fréquence. Ces mesures peuvent être filtrées pour affiner l'estimation du CET. L'ensemble du segment bord, par exemple d'un petit satellite, peut être un émetteur VHF cohérent en spectre étalé standard. L'antenne d'émission et l'oscillateur peuvent être standard, eux aussi. L'oscillateur doit toutefois être suffisamment stable, de façon à ce que sa contribution à la mesure de Vçr soit négligeable. Le segment bord décrit précédemment peut être couplé avec une autre fonction de ces satellites, les petits satellites en particulier : l'initialisation autonome des liaison spatiales directives. Les petits satellites mettent en général en jeu deux types de liaisons spatiales - non directives (à relativement faible débit, en VHF par exemple) - directives (à relativement haut débit, en bande S par exemple). L'invention peut donc utiliser les besoins déjà existants de satellites, par exemple de petits satellites, qui utilisent un émetteur en spectre étalé en VHF pour transmettre des données vers le sol et permettre de faire des mesures de distance et/ou de vitesse dudit satellite en effectuant le traitement des signaux reçus au sol, pour effectuer une mesure du retard ionosphérique. On peut envisager un fonctionnement bidirectionnel, c'est-à-dire une liaison montante suivie d'une liaison descendante, en utilisant un répéteur transparent ou transpondeur à l'intérieur du satellite. Dans ce cas, la station sol utilisatrice est une station d'émission/réception. Cette station contient à la fois le dispositif 10 de génération et d'émission des signaux, ainsi que le segment sol de réception 17 décrit précédemment. Dans ce cas, les signaux générés par les oscillateurs locaux 14 et 21 sont générés de façon cohérente, à partir de la même référence de fréquence. Dans le cas d'un tel fonctionnement bidirectionnel, on peut également utiliser deux bandes de fréquence (par exemple bande UHF pour la liaison montante et bande VHK pour la liaison descendante). REFERENCES [1] "Techniques et technologies des véhicules spatiaux" (CNES, Cépaduès-Editions, tome 1, module 6, pages 628 à 633) [2] "Single Frequency Ionosphère Determination Using GPS" de M.L. Trethewey, I. Catchpole et A. Hansla (IONGPS-93 Proc., septembre 1993, pages 1373 à 1381) [3] "Estimation of Absolute Ionospheric Delay Exclusively Through Single- Frequency GPS Measurements" de C.E. Cohen, B. Peervan et B.W. Parkinson (ION GPS-92 Proc., septembre 1992, pages 325 à 330) REVENDICATIONS 1. Système de mesure du contenu électronique total de la ionosphère, caractérisé en ce qu'il comprend - à bord d'un satellite (9), un émetteur (11) d'un signal cohérent en spectre étalé transmis en VHF ; - au sol, un dispositif de réception (18) et de traitement du signal VHF reçu permettant de mesurer le contenu électrique total (CET) de la ionosphère à l'aide d'une seule fréquence. 2. Système selon la revendication 1, caractérisé en ce que le segment bord est couplé avec une fonction d'initialisation autonome des liaisons spatiales directives, ou/et une fonction de transmission de secours. 3. Procédé de mesure du contenu électronique total de la ionosphère, à l'aide d'un signal en spectre étalé transmis à travers la ionosphère dans la bande VHF, ce procédé étant caractérisé en ce que qu'on émet à partir d'un satellite un signal constitué d'un code pseudoaléatoire dans la bande VHF. 4. Procédé selon la revendication 3, caractérisé en ce qu'on émet à partir d'une station sol un signal semblable au précédent, répété par un satellite dans la bande VHF. 5. Procédé selon la revendication 3, caractérisé en ce qu'on reçoit le signal transmis par le satellite dans la station sol, les mesures disponibles dans la station (vitesse de groupe, vitesse de phase) étant utilisées pour déterminer le contenu électronique total de la ionosphère. 6. Procédé selon la revendication 3, caractérisé en ce que les mesures disponibles dans la station (distance, Doppler) sont corrigées des effets de la ionosphère, à l'aide de la connaissance du contenu électronique total, ces mesures corrigées étant utilisées pour déterminer l'orbite dudit satellite.