FR 3045989 A1 20170623 FR 1502640 A 20151218 1 Méthode d'allocation de ressources radio dans un système de communications par satellites défilants avec contrainte de niveau d'interférences vers un système géostationnaire L'invention concerne le domaine des systèmes de communications par satellites et en particulier les systèmes utilisant des satellites défilants, par exemple des satellites sur orbite basse ou « low earth orbit » en anglais. L'invention porte sur une méthode d'allocation de ressources radio dans un système de communications par satellites défilants opérant dans des bandes de fréquence réservées à des systèmes géostationnaires. L'invention concerne la prise en compte, dans la méthode d'allocation de ressources radio, du niveau d'interférences engendré par le système à satellites défilants vers un système de communications, quelconque, utilisant des satellites géostationnaires. La réglementation relative aux radiocommunications par satellites permet le déploiement de systèmes non géostationnaires dans des bandes de fréquence auparavant réservées aux systèmes basés sur des satellites 20 géostationnaires. Cependant, la coexistence de deux systèmes de communications par satellites exploitant, totalement ou partiellement, les mêmes bandes de fréquence pose un problème de gestion des interférences engendrées par un système vers un autre. En l'occurrence, le système géostationnaire est 25 susceptible d'être perturbé, voire brouillé, par les transmissions entre satellites et stations au sol appartenant au système non géostationnaire. Pour garantir une qualité de service donnée pour un système de communications géostationnaire, l'union internationale des télécommunications (UIT ou ITU en anglais) a fixé un niveau d'interférences 30 maximum autorisé impactant les systèmes de télécommunications géostationnaire dans leur ensemble et provenant des systèmes coexistant 3045989 2 sur les mêmes bandes de fréquence et utilisant des satellites non géostationnaires. Pour garantir qu'un système non géostationnaire n'engendre pas un niveau d'interférences supérieur au seuil autorisé, une solution consiste à 5 fixer, pour l'ensemble des liaisons entre les terminaux au sol et les satellites non géostationnaires, un seuil d'angle de séparation à respecter par rapport à l'arc géostationnaire. L'arc géostationnaire désigne la vue de l'orbite géostationnaire depuis la Terre. Un terminal du système non géostationnaire est autorisé à établir une 10 liaison de communication (montante, descendante ou bidirectionnelle) avec un satellite non géostationnaire si et seulement si son angle de séparation est inférieur au seuil d'angle de séparation fixé et ce indépendamment des autres liaisons établies par les autres terminaux. Ainsi, chaque terminal voit le nombre de satellites non géostationnaires avec lesquels il peut 15 potentiellement communiquer se réduire. En particulier, les systèmes non géostationnaires peuvent utiliser des ressources spectrales en bande Ka et les terminaux de ces systèmes peuvent opérer avec des angles d'élévation très faibles. De telles conditions entrainent l'exigence d'utilisation de puissances de transmission très importantes. De tels niveaux de puissance nécessitent alors l'utilisation d'un angle de séparation de valeur élevée pour ne pas brouiller les systèmes géostationnaires co-existants sur la même bande Ka. Typiquement, une valeur de seuil d'angle de séparation de l'ordre de 7 degrés est utilisée pour satisfaire l'ensemble de ces contraintes. Malheureusement, l'utilisation d'un angle de séparation avec une valeur fixe et élevée a pour conséquence que la couverture géographique assurée par le système non géostationnaire s'en trouve diminuée. Une solution simple permettant d'assurer un niveau de service et une couverture géographique suffisants pour le système non géostationnaire tout en respectant le seuil d'angle de séparation préconisé, consiste à augmenter le nombre de satellites non géostationnaires dans la constellation du système. Ainsi, on assure un nombre moyen de satellites disponibles pour 3045989 3 chaque terminal tout en interdisant les liaisons terminaux-satellites qui ne sont pas compatibles du seuil d'angle de séparation. Cependant, l'augmentation du nombre de satellites non 5 géostationnaires présente des inconvénients de surcoûts pour la conception du système global et également de sous-optimalité car l'ensemble des ressources de communication disponibles n'est pas pleinement utilisé du fait de l'empêchement lié à la coexistence avec les systèmes géostationnaires. 10 Un procédé statique décrit dans la demande de brevet US 2003/0073404 Al propose d'éviter les alignements des liens radio fréquence terre/espace entre un nouveau système non géostationnaire (système B) et d'autres systèmes géostationnaires existants (systèmes A) partageant la même bande de fréquence. 15 Ce procédé s'appuie sur la prise en compte de contraintes prédictives telles que les orbites des satellites et les caractéristiques radio- fréquence des stations et des satellites des systèmes A et B. Ce procédé statique est valable pour des systèmes opérant dans des bandes de fréquence basses comme la bande Ku. En effet, pour ce type de 20 systèmes, l'aspect géométrique permet à lui seul de planifier de façon statique et déterministe les liens radio fréquence terre/espace à utiliser par le système B sans provoquer d'interférences sur les autres systèmes A existants et tout en assurant la qualité de service attendue pour le système B. 25 Cette stratégie est justifiée dès lors que les atténuations atmosphériques sont négligeables (ce qui est le cas des systèmes exploitant des bandes de fréquence basses). Ce procédé statique n'est par contre plus approprié pour un nouveau système B opérant dans des bandes de fréquence hautes comme la bande 30 Ka. En effet, l'utilisation de faibles élévations combinées à d'importantes atténuations atmosphériques, comme celles induites par la pluie, impose un 3045989 4 surdimensionnement excessif en terme de puissance requise au niveau bord et au niveau sol de l'infrastructure du système B afin d'assurer la qualité de service attendue. Une variante du procédé ci-dessus consiste à imposer que les liens 5 radio fréquence terre/espace du nouveau système B ne puissent être mis en oeuvre que lorsque les satellites du système B sont vus par les stations du système B avec un angle d'élévation élevé. Ceci permet de réduire la dynamique des atténuations atmosphériques même dans des bandes de fréquence hautes comme la bande Ka. 10 Néanmoins, cette stratégie impose un surdimensionnement en terme de nombre de satellites à déployer par le système B afin d'assurer le taux de visibilité requis par ce dernier. Ce surdimensionnement a un fort impact sur le coût et le déploiement des infrastructures sol et bord du système B. 15 La présente invention propose de remédier aux inconvénients précités par une méthode d'allocation de ressources radio, dans un système de communications par satellites défilants, qui définit, de façon dynamique, un seuil d'angle de séparation à respecter pour chaque terminal en fonction d'un niveau global d'interférences engendré vers les systèmes géostationnaires. 20 L'invention permet de maximiser la couverture géographique en minimisant le nombre de satellites dans la constellation du système non géostationnaire tout en garantissant que le niveau d'interférences induit vers les systèmes géostationnaires exploitant des bandes de fréquences totalement ou partiellement identiques ne dépasse pas le seuil 25 réglementaire. L'invention a pour objet une méthode d'allocation de ressources radio pour l'établissement de communications par satellites dans un premier système de communications comprenant une constellation de satellites 30 défilants et un premier ensemble de terminaux, la méthode comprenant les étapes suivantes, pour chaque terminal du premier ensemble : 3045989 5 déterminer, un seuil d'angle de séparation distinct à respecter vis-à-vis d'une constellation de satellites géostationnaires, un angle de séparation étant défini comme l'angle topocentrique le plus faible sous lequel un satellite défilant donné et un satellite 5 géostationnaire quelconque de l'arc géostationnaire sont vus depuis un terminal donné du premier ensemble, autoriser le terminal à établir une liaison de communication avec un satellite défilant si l'angle de séparation associé à ladite liaison est supérieur ou égal au seuil d'angle de séparation. 10 Selon un aspect particulier de l'invention, le seuil d'angle de séparation distinct à respecter, pour chaque terminal du premier ensemble est déterminé de manière à ce que l'ensemble des liaisons de communications autorisées engendrent, sur un système de communications par satellites géostationnaires, un niveau de puissance interférente conforme à une limite 15 réglementaire. Selon un aspect particulier de l'invention, la détermination d'un seuil d'angle de séparation comprend : le calcul d'au moins un niveau de puissance interférente engendrée par l'ensemble des liaisons de communications entre 20 les terminaux du premier ensemble et les satellites défilants (vers un système de communications par satellites géostationnaires, la recherche itérative, pour chaque terminal du premier ensemble, de la valeur minimale du seuil d'angle de séparation permettant de maintenir ledit au moins un niveau de puissance interférente 25 calculé conforme à une limite réglementaire. Selon un aspect particulier de l'invention, la détermination d'un seuil d'angle de séparation comprend une phase d'initialisation consistant à attribuer une valeur initiale au seuil d'angle de séparation, pour chaque terminal, et à déclarer chaque terminal candidat à l'abaissement de son seuil 30 d'angle de séparation. 3045989 6 Selon un aspect particulier de l'invention, ladite recherche itérative comprend, pour chaque terminal et à chaque itération, l'association dudit terminal avec un satellite défilant respectant le seuil d'angle de séparation déterminé à l'itération courante et pour lequel l'angle de séparation est le 5 plus proche dudit seuil, les itérations étant poursuivies tant qu'au moins un terminal est candidat à l'abaissement de son seuil d'angle de séparation. Selon un aspect particulier de l'invention, les liaisons de communication entre les terminaux du premier ensemble et les satellites défilants sont bidirectionnelles ou unidirectionnelles montantes et ladite recherche itérative 10 comprend le calcul d'un premier niveau de puissance interférente engendrée par l'ensemble des liaisons de communications montantes vers l'arc géostationnaire pour un ensemble de positions supposées de satellites géostationnaires. Selon un aspect particulier de l'invention, le calcul dudit premier niveau 15 de puissance interférente est déterminé en choisissant, pour chaque satellite géostationnaire, une direction de pointage qui permet de maximiser la valeur dudit premier niveau de puissance interférente. Selon un aspect particulier de l'invention, ladite recherche itérative comprend la recherche d'au moins un satellite géostationnaire pour lequel le 20 premier niveau de puissance interférente calculé n'est pas conforme à la dite limite réglementaire. Selon un aspect particulier de l'invention, les liaisons de communications entre les terminaux du premier ensemble et les satellites défilants sont unidirectionnelles descendantes et ladite recherche itérative comprend le 25 calcul d'un second niveau de puissance interférente engendrée par l'ensemble des liaisons de communication descendantes vers un ensemble de positions supposées de terminaux appartenant à un système de communications par satellites géostationnaires. Selon un aspect particulier de l'invention, les liaisons de communication 30 entre les terminaux du premier ensemble et les satellites défilants sont bidirectionnelles et ladite recherche itérative comprend le calcul d'un second 3045989 7 niveau de puissance interférente engendrée par l'ensemble des liaisons de communication descendantes vers un ensemble de positions supposées de terminaux appartenant à un système de communications par satellites géostationnaires. 5 Selon un aspect particulier de l'invention, le calcul dudit second niveau de puissance interférente est déterminé en choisissant, pour chaque terminal appartenant à un système de communications par satellites géostationnaires, une direction de pointage qui permet de maximiser la valeur dudit second niveau de puissance interférente. 10 Selon un aspect particulier de l'invention, ladite recherche itérative comprend la recherche d'au moins un terminal appartenant à un système de communications par satellites géostationnaires pour lequel le second niveau de puissance interférente calculé n'est pas conforme audit seuil réglementaire. 15 Selon un aspect particulier de l'invention, ladite recherche itérative comprend, pour chaque premier ou second niveau de puissance interférente calculé non conforme audit seuil réglementaire, la recherche d'au moins un couple (terminal du premier ensemble, satellite défilant) dont la contribution, dans la valeur du premier ou du second niveau de puissance interférente, est 20 la plus élevée. Selon un aspect particulier de l'invention, ladite recherche itérative comprend, pour chaque terminal du couple (terminal du premier ensemble, satellite défilant) retenu, le rétablissement du seuil d'angle de séparation déterminé à l'itération précédente et le retrait dudit terminal de la liste des 25 candidats à l'abaissement du seuil d'angle de séparation. Selon un aspect particulier de l'invention ladite recherche itérative comprend, pour au moins un terminal du couple (terminal du premier ensemble, satellite défilant) retenu, l'adaptation des caractéristiques RF d'au moins une liaison impliquant ledit terminal et un satellite défilant, de sorte à 30 diminuer le niveau de puissance interférente de l'ensemble des liaisons de 3045989 8 communications autorisées sur un système de communications par satellites géostationnaires. Selon un aspect particulier de l'invention, la recherche itérative comprend une vérification de la conformité du premier niveau de puissance interférente 5 et/ou du second niveau de puissance interférente et, si la conformité est vérifiée, l'abaissement du seuil d'angle de séparation d'une valeur prédéterminée, pour chaque terminal candidat à l'abaissement de son seuil d'angle de séparation. Selon un aspect particulier de l'invention la recherche itérative est arrêtée 10 lorsqu'aucun terminal n'est plus candidat à l'abaissement de son seuil d'angle de séparation. L'invention a aussi pour objet un dispositif d'allocation de ressources radio pour l'établissement de communications par satellites comprenant des moyens, incluant au moins un processeur et une mémoire, configurés pour 15 mettre en oeuvre la méthode d'allocation de ressources radio selon l'invention. L'invention a aussi pour objet un programme d'ordinateur comportant des instructions pour l'exécution de la méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon l'invention, 20 lorsque le programme est exécuté par un processeur. L'invention a aussi pour objet un support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution de la méthode d'allocation de ressources radio pour l'établissement des communications par satellites selon l'invention, 25 lorsque le programme est exécuté par un processeur. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés qui représentent : 30 - La figure 1, un synoptique représentant le contexte d'application de l'invention, 3045989 9 La figure 2, un organigramme détaillant les étapes d'une méthode de conception d'un système de communications par satellites non géostationnaires selon l'invention, - La figure 3, un organigramme détaillant les étapes d'une méthode 5 d'allocation de ressources radio selon l'invention, - La figure 4, un organigramme détaillant les étapes d'une fonction de détermination dynamique d'un seuil d'angle de séparation, selon un premier mode de réalisation de l'invention, - La figure 5a, un schéma illustrant un exemple de calcul d'interférences 10 global relatif aux voies montantes, La figure 5b, un schéma illustrant un exemple de calcul d'interférences global relatif aux voies descendantes, La figure 6, un organigramme détaillant les étapes d'une fonction de détermination dynamique d'un seuil d'angle de séparation, selon un 15 deuxième mode de réalisation de l'invention appliqué aux liaisons montantes, - La figure 7, un organigramme détaillant les étapes d'une fonction de détermination dynamique d'un seuil d'angle de séparation, selon un deuxième mode de réalisation de l'invention appliqué aux liaisons 20 descendantes. La figure 1 représente schématiquement un exemple de scénario d'application de l'invention comprenant un premier système de communications par satellites défilants et un second système de 25 communications par satellites qui est un système géostationnaire. Le premier système de communications comprend une pluralité de satellites SATNGS01,SAT_NGS02 et un ensemble de terminaux au sol TER_NGS01,TER_NGS02. Les terminaux peuvent être fixes ou mobiles et sont aptes à communiquer avec chaque satellite du système selon une 30 liaison montante ou une liaison descendante ou une liaison bidirectionnelle LR1,LR2. De même, le second système de communications géostationnaire 3045989 10 comprend une pluralité de satellites géostationnaires et un ensemble de terminaux. Pour simplifier, sur la figure 1 ne sont représentés qu'un satellite SAT_GS0 et un terminal TER_GS0 constituant le second système de communications et communiquant l'un avec l'autre via une liaison par 5 satellite LRGEO. Un terminal peut être tout équipement de réseau capable d'échanger et de communiquer via une liaison sans fil avec un satellite. Il peut s'agir en particulier d'un ordinateur fixe ou portable, d'un téléphone fixe ou portable, d'un assistant personnel numérique, d'un serveur ou encore d'un modem 10 d'accès à l'Internet par satellite. Le premier système de communications par satellites et le second système de communications par satellites sont aptes à établir des communications dans des bandes de fréquence totalement ou partiellement identiques. Un exemple de bande de fréquences utilisée est la bande Ka. Un 15 autre exemple est la bande Ku. Le premier système de communications par satellites comprend en outre une station satellitaire ST_RRM comprenant un sous-système de contrôle chargé de gérer l'allocation de ressources de communication au sein du système. La présente invention peut être implémentée dans le sous- 20 système de contrôle. L'invention peut être implémentée à partir d'éléments matériel et/ou logiciel. Elle peut notamment être mise en oeuvre en tant que programme d'ordinateur comportant des instructions pour son exécution. Le programme d'ordinateur peut être enregistré sur un support d'enregistrement lisible par 25 un processeur. Le support peut être électronique, magnétique, optique ou électromagnétique. L'invention peut aussi être implémentée par un dispositif du type sous-système de contrôle embarqué dans une station satellite, le dispositif comprenant un processeur et une mémoire. Le processeur peut être un 30 processeur générique, un processeur spécifique, un circuit intégré propre à une application (connu aussi sous le nom anglais d'ASIC pour « Application- 3045989 11 Specific Integrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour « Field-Programmable Gate Array »). Le dispositif peut utiliser un ou plusieurs circuits électroniques dédiés 5 ou un circuit à usage général. La technique de l'invention peut se réaliser sur une machine de calcul reprogrammable (un processeur ou un micro contrôleur par exemple) exécutant un programme comprenant une séquence d'instructions, ou sur une machine de calcul dédiée (par exemple un ensemble de portes logiques comme un FPGA ou un ASIC, ou tout autre 10 module matériel). A titre d'exemple d'architecture matérielle adaptée à mettre en oeuvre l'invention, un dispositif selon l'invention peut comporter un bus de communication auquel sont reliés une unité centrale de traitement ou microprocesseur (CPU, acronyme de « Central Processing Unit » en 15 anglais), une mémoire morte (ROM, acronyme de « Read Only Memory » en anglais) pouvant comporter les programmes nécessaires à la mise en oeuvre de l'invention; une mémoire vive ou mémoire cache (RAM, acronyme de « Random Access Memory » en anglais) comportant des registres adaptés à enregistrer des variables et paramètres créés et modifiés au cours de 20 l'exécution des programmes précités; et une interface de communication ou E/S (I/O acronyme de « Input/ouput » en anglais) adaptée à transmettre et à recevoir des données. On introduit à présent une définition d'un angle, qui est identifié sur la 25 figure 1. Cet angle est utilisé par la suite pour expliciter les modes de réalisation de l'invention. On considère un angle at désigné par la suite par angle topocentrique. Pour rappel, un angle topocentrique est un angle sous lequel deux points donnés sont vus d'un point spécifique de la Terre. L'angle topocentrique 30 particulier at considéré dans le contexte de l'invention est un angle formé par la direction entre un terminal et un satellite du premier système et par la 3045989 12 direction entre ce même terminal et un satellite du système géostationnaire. Sur la figure 1, on a, à titre d'exemple, représenté l'angle topocentrique atil associé au terminal TER_NGS0 1 et aux satellites SAT_NGS01,SAT_GSO. On a représenté également l'angle topocentrique at22 associé au terminal 5 TER NGS02 et aux satellites SATNGS02,SAT_GSO. On définit également la notion d'angle de séparation, pour une liaison entre un terminal du premier système TER_NGS01,TER_NGS02 et un satellite défilant du premier système SAT_NGS01,SAT_NGS02 comme étant l'angle topocentrique le plus faible en faisant varier la position d'un 10 satellite géostationnaire sur l'ensemble des points possibles de l'arc géostationnaire. Autrement dit, l'angle de séparation est défini comme étant l'angle minimum sous lequel un satellite défilant de la première constellation et n'importe quel satellite géostationnaire sont vus depuis un terminal du 15 premier système de communications. On parle d'angle de séparation car cet angle définit la proximité entre la liaison reliant un terminal du premier système avec un satellite défilant de la première constellation et les satellites géostationnaires. Plus cet angle est faible, plus il y a de probabilité que cette liaison interfère sur les satellites 20 géostationnaires mais aussi sur les terminaux au sol des systèmes géostationnaires. En particulier, l'interférence générée sur les satellites géostationnaires provient du signal émis par un terminal de la première constellation vers l'espace, autrement dit sur la liaison montante Terre vers espace. Sur la figure 1, on a schématiquement représenté les interférences 25 INT_ M _1 générées par le terminal TER_NGS0 1 sur le satellite géostationnaire SAT_ GSO ainsi que les interférences INT_M_2 générées par le terminal TER_NGSO2 sur le satellite géostationnaire SAT_GSO. De même, l'interférence générée sur les terminaux au sol TER_GS0 appartenant à des systèmes géostationnaires provient du signal émis par un 30 satellite défilant du premier système de communications, autrement dit de la liaison descendante espace vers Terre. Sur la figure 1, on a identifié 3045989 13 également les interférences INT_D_1 générées par le satellite SAT_NGS01 vers le terminal TER_ GSO du système géostationnaire ainsi que les interférences INT_D_2 générées par le satellite SAT_NGSO2 vers ce même terminal TER GSO. _ 5 Comme indiqué en préambule, une solution permettant de respecter un niveau d'interférences réglementaire consiste à fixer un seuil pour cet angle de séparation, par exemple un seuil égal à 7°, en dessous duquel, la liaison entre un terminal et un satellite du premier système est interdite obligeant alors à sélectionner un autre satellite, pour l'établissement d'une 10 liaison avec le terminal, de sorte à respecter un angle de séparation supérieur au seuil fixé. Il convient de noter que le terme « angle de séparation » peut également être remplacé par « angle d'évitement » en français ou « avoidance angle » ou « separation angle » en anglais. 15 L'invention propose une méthode de configuration dynamique des seuils appliqués aux angles de séparation pour améliorer l'allocation de ressources de communication disponibles au sein du premier système de communications tout en garantissant un respect du niveau d'interférences 20 réglementaire engendré sur les systèmes de communications par satellites géostationnaires. La figure 2 schématise un organigramme détaillant les étapes d'une méthode de conception d'un système de communications par satellites selon 25 l'invention. Cette méthode a pour fonction globale la gestion statique de ressources radio du système. Une telle fonction est usuellement désignée par l'expression anglo-saxonne « static radio resources management ». Dans une première étape 201, un opérateur de satellites fournit une demande de besoins d'un système de communications qui comprend des 30 besoins en qualité de service, disponibilité globale, débit global, taux 3045989 14 d'erreurs ou autres caractéristiques générales de fonctionnement d'un tel système. Dans une deuxième étape 202, on évalue l'impact des conditions locales atmosphériques pour chaque terminal susceptible d'être intégré au 5 système. Cette évaluation peut se faire en caractérisant le niveau d'atténuation du signal reçu ou émis par un terminal en fonction de l'élévation du terminal. Dans une troisième étape 203, on effectue un calcul de bilan de liaison à la fois pour les liaisons montantes et les liaisons descendantes du 10 système. Le calcul du bilan de liaison prend notamment comme paramètre d'entrée le niveau maximum d'atténuation en fonction de l'élévation pour chaque terminal. Pour la liaison montante, le bilan de liaison peut être évalué par un calcul de la puissance isotrope rayonnée équivalente (PIRE) requise au niveau du terminal. Pour la liaison descendante, le bilan de liaison peut 15 être évalué par un calcul du rapport signal à bruit G/T requis au niveau du terminal. Dans une quatrième étape 204, on évalue les caractéristiques particulières des terminaux à déployer pour satisfaire le bilan de liaison. De telles caractéristiques peuvent notamment inclure le type d'antenne ou le 20 type d'amplificateur embarqué dans le terminal. En outre, des fonctionnalités supplémentaires de traitement de signal peuvent être incorporées dans les terminaux du système afin de lutter contre les atténuations atmosphériques du signal. Dans une cinquième étape 205, on vérifie que les besoins de 25 l'opérateur peuvent être satisfaits en prenant en compte les caractéristiques des terminaux sélectionnés à l'issue de l'étape 204 et en prenant en compte les plages d'atténuation du signal en fonction des plages d'élévation des différents terminaux. Si les besoins des connexions 200 ne peuvent être satisfaits, une 30 négociation est établie avec l'opérateur pour réduire les exigences en termes de qualité de service ou de débit. 3045989 15 Si au contraire, les besoins des connexions 200 sont satisfaits alors les besoins d'établissement de communications entre les terminaux et les satellites sont validés et on passe, par l'intermédiaire d'une étape 206 vers un processus d'allocation de ressources radio. L'étape 206 consiste à lister, 5 pour chaque terminal du système, ses caractéristiques radio-fréquence ainsi que le débit requis. La figure 3 détaille, sur un organigramme, les étapes de la méthode d'allocation de ressources selon l'invention. 10 Cette méthode a pour fonction globale une gestion dynamique de ressources radio du système ou « dynamic radio resources management » en anglais. En particulier, l'allocation de ressources radio est réalisée au fil du temps en prenant en compte les trajectoires dynamiques des satellites non géostationnaires. 15 La méthode débute à l'étape 206 qui analyse le nombre, le type, l'emplacement, les caractéristiques radio fréquence et le débit requis de chacun des terminaux du système déterminés à l'issue de la méthode de conception décrite à la figure 2. Dans une étape 208, on identifie les ressources radio disponibles pour 20 chaque terminal. Autrement dit, on identifie les satellites de la constellation du système qui sont en ligne de visée de chaque terminal et qui respectent un seuil d'angle de séparation déterminé, pour chaque terminal, via l'étape 207 décrite plus en détail ci après. Pour chaque satellite, l'identification des ressources radio peut également comprendre l'identification des porteuses 25 fréquentielles qui peuvent être utilisées pour communiquer et le nombre et le type des antennes embarquées à bord du satellite, en particulier la polarisation des antennes. A l'issue de l'étape 208, on obtient une liste de satellites et de ressources radio associées qui peuvent potentiellement être utilisées par 30 chaque terminal. 3045989 16 Dans une étape 209, on effectue, pour chaque terminal, un tri dans la liste des satellites obtenus afin d'établir un ordre de priorité pour l'allocation de ressources radio. Les priorités sont définies à partir d'un score attribué à chaque liaison, le score étant déterminé notamment en fonction d'un calcul 5 de bilan de liaison qui peut être fourni par l'étape 203 de la méthode de conception décrite à la figure 2. Dans une étape 210, on effectue ensuite une allocation de ressources radio pour chaque terminal en fonction du score de priorité attribué à l'étape 10 209. Les étapes 208 et 210 sont réalisées notamment à partir des contraintes opérationnelles du système 301. Dans une étape 211, on agrège toutes les liaisons radio actives à 15 partir des résultats de l'étape 210 d'allocation de ressources. Autrement dit, on identifie toutes les liaisons montantes et les liaisons descendantes entre terminaux et satellites afin d'avoir une vision globale des liaisons radio actives. 20 Dans une étape 207, on détermine ensuite dynamiquement, selon l'invention, pour chaque terminal, un seuil d'angle de séparation qui permet de respecter les besoins globaux en termes de débits requis tout en respectant les contraintes de niveau d'interférences sur les systèmes de communications par satellites géostationnaires opérant sur les mêmes 25 bandes de fréquence. Cette étape 207 prend notamment en entrée, les caractéristiques de la constellation des satellites du système, notamment les positions des satellites défilants, mais aussi les contraintes réglementaires de niveau d'interférences et les caractéristiques des systèmes géostationnaires que l'on souhaite préserver. 30 La détermination d'un seuil d'angle de séparation à respecter pour chaque terminal prend notamment en compte un calcul de niveau 3045989 17 d'interférences global du premier système vers l'arc géostationnaire et vers les stations sol. L'étape 207 fournit ensuite à l'étape 208, pour la prochaine itération 5 de la méthode, un seuil d'angle de séparation à respecter pour chaque terminal pour établir une liaison avec un satellite, de sorte à respecter les réglementations imposées pour le niveau d'interférences global. Les étapes 208,209,210 sont alors itérées en éliminant de la liste des satellites disponibles pour chaque terminal, ceux qui ne permettent pas 10 d'obtenir un angle de séparation supérieur ou égal au seuil d'angle de séparation. L'étape 210 d'allocation de ressources radio peut ainsi conduire à attribuer à un terminal d'autres satellites que ceux choisis initialement ou à la précédente itération. 15 Selon un aspect particulier de l'invention, les liaisons entre terminaux et satellites peuvent être bidirectionnelles ou unidirectionnelles. Un terminal peut être autorisé à communiquer avec deux satellites différents sur la voie montante et sur la voie descendante ou au contraire, les deux voies peuvent être liées en une seule liaison bidirectionnelle entre un terminal et un 20 satellite. Si les liaisons entre terminaux et satellites sont bidirectionnelles, alors l'étape 207 de la méthode selon l'invention consiste en la détermination d'un seuil d'angle de séparation associé à chaque terminal. Si, au contraire, un terminal peut communiquer sur une liaison 25 montante avec un satellite et sur une liaison descendante avec un autre satellite, alors l'étape 207 de la méthode est dupliquée pour déterminer un premier seuil d'angle de séparation applicable à chaque terminal pour les liaisons montantes et un second seuil d'angle de séparation applicable à chaque terminal pour les liaisons descendantes. 30 3045989 18 La méthode d'allocation de ressources radio selon l'invention décrite à la figure 3 peut être exécutée avec une cadence donnée ou en réaction à un évènement de la vie du système, par exemple l'entrée ou la sortie d'un terminal. 5 La figure 4 illustre un premier exemple de réalisation, selon l'invention, de la méthode de détermination dynamique d'un seuil d'angle de séparation pour chaque terminal, exécutée à l'étape 207. Selon ce premier exemple, les liaisons entre terminaux et satellites du premier système de communications 10 sont considérées bidirectionnelles, c'est à dire qu'un terminal communique avec un même satellite sur la liaison montante et sur la liaison descendante. Dans une première étape 401, les terminaux actifs du premier système de communication sont identifiés ainsi que les satellites défilants de la première constellation avec lesquels ils peuvent communiquer selon la liste établie à 15 l'issue de l'étape 208 de la méthode d'allocation de ressources radio décrite à la figure 3. Pour chaque terminal actif, on fixe une valeur prédéfinie initiale de seuil à respecter pour l'angle de séparation. Ce seuil peut être, par exemple fixé à une valeur de 7°. Dans une étape 402, tous les terminaux actifs sont 20 ensuite déclarés candidats pour l'abaissement potentiel de leur seuil d'angle de séparation. Cette déclaration peut se faire, par exemple, par le biais d'un vecteur ou d'une table dans laquelle sont listés les terminaux actifs identifiés par un identifiant numérique et une information binaire indiquant si le terminal est candidat ou non à l'abaissement potentiel de son seuil d'angle de 25 séparation. Dans une étape 403, on associe à chaque terminal du premier système, un satellite de la première constellation qui respecte le seuil d'angle de séparation fixé à l'étape précédente et qui correspond à l'angle de séparation le plus proche du seuil. A l'issue de l'étape 403, on obtient donc, 30 pour chaque terminal, un satellite associé de sorte que la liaison entre ce terminal et ce satellite respecte le seuil d'angle de séparation tout en 3045989 19 minimisant la valeur de l'angle de séparation. Dans certains cas, aucun satellite ne peut être retenu car aucun satellite ne permet de respecter le seuil d'angle de séparation initial fixé. Ce scénario est traité dans une variante de réalisation de l'invention décrite plus loin. 5 La méthode se poursuit ensuite avec les étapes 404 et 405 dans lesquelles un calcul de niveau d'interférences global pire cas est réalisé. La notion de calcul « pire cas » consiste à réaliser ce calcul pour un angle de séparation minimum, identifié à l'étape précédente 403, le plus proche du 10 seuil considéré pour chaque terminal. Selon le premier exemple de réalisation de la méthode de détermination d'un seuil d'angle de séparation décrit à la figure 4, un premier calcul de niveau d'interférences global pire cas pour les liaisons montantes entre terminaux et satellites du premier système de communication est réalisé à l'étape 404 et un second calcul de 15 niveau d'interférences global pire cas pour les liaisons descendantes entre satellites et terminaux du premier système de communication est réalisé à l'étape 405. Le premier et le second calcul de niveau d'interférences sont effectués en agrégeant les contributions de niveau d'interférences pour chaque terminal actif (pour la liaison montante) ou pour chaque satellite actif 20 (pour la liaison descendante). Le premier calcul de niveau d'interférences concerne les interférences générées par les émissions des terminaux actifs vers l'espace sur les satellites de l'arc géostationnaire. Il s'agit des interférences du type de celles référencées INT_M_1,INT_M_2 sur la figure 1. 25 Le second calcul de niveau d'interférences concerne les interférences générées par les émissions des satellites actifs vers la Terre sur les terminaux appartenant aux systèmes de communications par satellites géostationnaires. Il s'agit des interférences du type de celles référencées INT_D_1, INT_D_2 sur la figure 1. 30 304 5 9 89 20 Le niveau d'interférences global généré peut être calculé par le biais d'un calcul de puissance surfacique équivalente ou « equivalent power fluxdensity » en anglais. Le calcul de puissance surfacique équivalente peut, par exemple, être 5 effectué à l'aide de la méthode décrite dans la recommandation UIT-R S.1325-3 de l'Union Internationale des Télécommunications qui est reprise ici. Selon cette recommandation la puissance surfacique équivalente epfd est définie comme étant la somme des puissances surfaciques produites au 10 niveau d'une station de réception d'un système brouillé (située à la surface de la Terre ou en orbite, selon le cas) par toutes les stations d'émission d'un système brouilleur, compte tenu de la discrimination hors-axe d'une antenne de réception de référence que l'on suppose pointée vers sa direction nominale. 15 Dans le cas présent, le système brouilleur est le premier système de communications par satellites défilants et le système brouillé est un système de communications par satellites géostationnaires. La figure 5a représente, sur un schéma, un terminal TER_NGSO et deux 20 satellites SAT NGS01, SAT NGS02 du premier système de communications ainsi qu'un terminal TER_GSO et un satellite SAT_GS0 d'un système géostationnaire. Le paragraphe suivant explicite un exemple de calcul 404 de la puissance surfacique équivalente epfd_up calculée pour les interférences produites par les émissions du terminal TER_NGS0 vers le 25 satellite SAT GSO du système géostationnaire. Pour chaque position de l'arc géostationnaire, la puissance surfacique équivalente epfd_up pour les liaisons montantes (premier calcul de niveau d'interférence 404) peut être calculée par l'intermédiaire de la relation (1) 30 suivante : 3045989 21 =Na (1) epfd_up = 10 logio Eiop,,i0 Gt(sot,) Gr(s02.1 47c Ri2 Gr où: Na est le nombre de terminaux du premier système de communication visibles depuis un satellite géostationnaire (sur la figure 5a, un seul terminal 5 TER_NGS0 est représenté) i est l'indice du terminal considéré appartenant au premier système de communication, Pi est la puissance RF en entrée de l'antenne du terminal TER_NGSO considéré appartenant au premier système de communication (dBW), 10 G, (vu) est le gain de l'antenne d'émission du terminal TER_NGSO considéré appartenant au premier système de communication en direction d'un satellite géostationnaire SAT_GSO, Gr(e2i) est le gain de l'antenne de réception d'un satellite géostationnaire SAT_GSO en direction du ième terminal considéré du premier 15 système de communication TER_NGSO, G,,,, est le gain maximum de l'antenne d'un satellite géostationnaire SAT_GSO, (pli est l'angle entre l'axe de visée du terminal TER_NGS0 considéré (indice i) appartenant au premier système de communication et la direction 20 entre le terminal TER_NGSO et un satellite géostationnaire SAT_GS0 (angle représenté sur la figure 5a), (2i est l'angle entre l'axe de visée d'un satellite géostationnaire SAT_GS0 et la direction entre ce satellite SAT GSO et le ième terminal TER_NGSO considéré du premier système de communication (angle 25 représenté sur la figure 5a), Ri est la distance entre le terminal TER NGSO considéré du premier système de communication et le satellite géostationnaire SAT_GSO. 3045989 22 Le calcul de La puissance surfacique équivalente epfd_up pour les liaisons montantes est réalisé pour tous les points de l'arc géostationnaire avec un échantillonnage approprié. La direction de pointage du satellite géostationnaire SAT_GSO, autrement dit la valeur de l'angle cp2; est choisie 5 de sorte à maximiser la valeur du niveau d'interférence epfd_up. La direction de pointage du terminal du premier système TER_NGSO, autrement dit la valeur de l'angle (pi;, est donnée par la direction entre ce terminal TER_NGSO et le satellite SAT_NGS0 sélectionné à l'étape 403 (satellite qui permet de minimiser l'angle de séparation). 10 La valeur du seuil d'angle de séparation a une influence à la fois sur l'angle (pi; et sur l'angle cp2; puisque cette valeur définit l'éloignement relatif du satellite défilant SAT_NGSO par rapport au satellite géostationnaire SAT_GSO. En conséquence, la valeur du seuil d'angle de séparation a une influence sur le gain de l'antenne de réception d'un satellite géostationnaire 15 SAT_GS0 dans la direction de visée du terminal TER_NGS0 du premier système et sur le gain de l'antenne d'émission du terminal TER_NGS0 du premier système dans la direction de visée du satellite géostationnaire SAT_GSO. 20 La figure 5b représente, sur un schéma, deux terminaux TER_NGS01,TER_NGS02 et un satellite SAT_NGSO du premier système de communications ainsi qu'un terminal TER GSO et un satellite SATGS0 d'un système géostationnaire. Le paragraphe suivant explicite un exemple de calcul 405 de la puissance surfacique équivalente epfd_down calculée pour 25 les interférences produites par les émissions du satellite SAT_NGS0 vers le terminal TER GSO du système géostationnaire. De la même façon, pour chaque terminal TER_GS0 placé sur un point de la surface de la terre, la puissance surfacique équivalente epfd_down pour les 30 liaisons descendantes (second calcul de niveau d'interférence 405) peut être 3045989 23 calculée par l'intermédiaire de la relation (2), identique à la relation (1) mais appliquée aux variables suivantes : [i=±b1"°G't(e ) G' , (9'2 ) epfd down =10loglo 10, , j (2) .1.-1 4n R'',/ G' ,. . 5 Nb est le nombre de satellites de la première constellation visibles depuis un terminal appartenant au système géostationnaire (sur la figure 5b, on a représenté, à titre purement illustratif, un seul satellite SAT_NGSO), Fi est la puissance RF en entrée de l'antenne du satellite SAT_NGS0 considéré appartenant à la première constellation (dBW), 10 G', (q" est le gain de l'antenne d'émission du satellite SAT_NGSO considéré appartenant à la première constellation en direction d'un terminal TER_GS0 appartenant au système géostationnaire, G'r (ço12i ) est le gain de l'antenne de réception d'un terminal TER_GSO appartenant au système géostationnaire en direction du ième satellite 15 SAT_NGSO considéré de la première constellation, G'r X est le gain maximum de l'antenne d'un terminal TER GSO du système géostationnaire, (p'ii est l'angle entre l'axe de visée du satellite SAT_NGSO considéré (indice i) appartenant à la première constellation et la direction entre ce 20 satellite SATNGS0 et un terminal TER_GS0 appartenant au système géostationnaire (angle représenté sur la figure 5b), ('2i est l'angle entre l'axe de visée d'un terminal TER_GSO appartenant au système géostationnaire et la direction entre ce terminal TER_GS0 et le ième satellite SAT_NGS0 considéré de la première 25 constellation (angle représenté sur la figure 5b), R'i est la distance entre le satellite SAT_NGS0 considéré de la première constellation et le terminal TER_GS0 du système géostationnaire. 3045989 24 La relation (2) est appliquée pour un ensemble de positions supposées de terminaux d'un système géostationnaire selon un maillage prédéfini et approprié de la surface de la Terre. La direction de pointage du terminal TER_GSO appartenant au système géostationnaire, autrement dit la 5 valeur de l'angle (p'21 est choisie de sorte à maximiser la valeur du niveau d'interférence epfd_down. La direction de pointage du satellite SAT_NGSO, autrement dit la valeur de l'angle est la direction de pointage du satellite SAT_NGS0 sélectionné à l'étape 403 (satellite qui permet de minimiser l'angle de séparation) vers le terminal TER_NGS0 considéré. 10 La valeur du seuil d'angle de séparation a une influence à la fois sur l'angle (p'1; et sur l'angle ('2i puisque cette valeur définit l'éloignement relatif du satellite défilant SAT_NGSO par rapport au satellite géostationnaire SAT_GSO. En conséquence, la valeur du seuil d'angle de séparation a une influence sur le gain de l'antenne de réception d'un terminal TER_GSO 15 appartenant au système géostationnaire dans la direction de visée du satellite SAT_NGSO du premier système et sur le gain de l'antenne d'émission du satellite SAT_NGS0 du premier système dans la direction de visée du terminal TER_GSO. 20 La méthode se poursuit avec l'application d'une étape 406 consistant en un test de conformité des niveaux d'interférences calculés aux étapes 404 et 405 avec des limites spécifiées réglementairement, par exemple par l'Union Internationale des Télécommunications. Dans l'étape 406, on applique un premier test de conformité pour les 25 liaisons montantes consistant à comparer le premier niveau d'interférences global Epfd_up avec une première limite réglementaire. De même, on applique un second test de conformité pour les liaisons descendantes consistant à comparer le second niveau d'interférences global Epfd_down avec une seconde limite réglementaire. 30 Les tests de conformité de l'étape 406 peuvent consister plus généralement à évaluer si le premier niveau d'interférences global Epfd_up 3045989 25 et le second niveau d'interférences global Epfd_down sont acceptables ou non du point de vue du fonctionnement des systèmes géostationnaires. Si les deux tests de conformité réalisés à l'étape 406 sont conformes, 5 autrement dit si les deux niveaux d'interférences globaux sont chacun inférieurs aux limites réglementaires sur tout l'arc géostationnaire et pour toutes les hypothèses de positions de terminaux d'un système géostationnaire, alors la méthode se poursuit avec une étape 407 qui consiste à diminuer, d'une valeur prédéterminée, la valeur du seuil d'angle 10 de séparation pour chaque terminal candidat à l'abaissement de son seuil. A la suite de l'étape 407, la méthode redémarre à l'étape 403 pour une nouvelle itération. L'abaissement du seuil d'angle de séparation est ici possible car le niveau global d'interférences généré en fixant le seuil d'angle de séparation à 15 la valeur initiale établie à la première étape 401, ne dépasse pas les limites réglementaires. Si au contraire, l'un des deux tests de conformité réalisés à l'étape 406 n'est pas conforme, la méthode se poursuit avec les étapes 408 et/ou 20 409. Si le premier test de conformité, concernant les liaisons montantes, n'est pas conforme, autrement dit si le premier niveau d'interférences global Epfd_up est supérieur au premier seuil réglementaire pour au moins un point de l'arc géostationnaire, la méthode se poursuit avec une étape 408 qui 25 consiste à identifier les satellites géostationnaires, autrement dit les points de l'arc géostationnaire, qui subissent un niveau d'interférences non conforme. Si le second test de conformité, concernant les liaisons descendantes, n'est pas conforme, autrement dit si le second niveau d'interférences global Epfd_down est supérieur au second seuil réglementaire pour au moins une 30 position supposée d'un terminal d'un système géostationnaire, la méthode se poursuit avec une étape 409 qui consiste à identifier les terminaux des 3045989 26 systèmes géostationnaires, autrement dit les points du maillage de la Terre, qui subissent un niveau d'interférences non conforme. La méthode se poursuit avec l'étape 410 dans laquelle, pour chaque 5 non-conformité relevée, on identifie les couples (terminal, satellite) du premier système (tels que créés à l'issue de l'étape 403) dont la contribution, dans le calcul d'interférences global, est la plus élevée. Cette recherche est effectuée indépendamment pour le calcul de niveau d'interférences associé aux liaisons montantes epfd_up et pour le calcul de niveau d'interférences 10 associé aux liaisons descendantes epfd_down. Cette recherche peut se faire en évaluant séparément les termes des sommes mises en jeu dans les relations (1) et (2), autrement dit les niveaux d'interférences calculés pour chacun des Na terminaux (pour les liaisons montantes) ou chacun des Nb satellites (pour les liaisons descendantes) et en recherchant les couples 15 (terminal, satellite) qui produisent le niveau d'interférences le plus élevé pour chacun des deux tests de conformité. A l'issue de l'étape 410, on obtient donc une liste de terminaux du premier système qui présentent un seuil d'angle de séparation qui engendre, un niveau d'interférences non réglementaire sur les systèmes 20 géostationnaires. Dans une étape 411, pour chaque terminal de la liste identifiée à l'étape 410, on rétablit la valeur du seuil d'angle de séparation à la valeur fixée à l'itération précédente puis on retire le terminal de la liste des 25 terminaux candidats à l'abaissement de leur seuil d'angle de séparation. Dans une étape 412, on contrôle si la liste des terminaux candidats à l'abaissement de leur seuil d'angle de séparation est vide, si oui on passe à l'étape 413 qui consiste à mettre à jour, pour chaque terminal, la valeur du 30 seuil d'angle de séparation à la dernière valeur fixée. 3045989 27 Si cette liste 412 n'est pas vide, alors la méthode retourne à l'étape 403 pour une nouvelle itération. Selon une variante du premier mode de réalisation de l'invention tel 5 que décrit à la figure 4, l'étape 411 est modifiée pour prendre en compte des cas de figure où, lorsque l'on rétablit un seuil d'angle de séparation à sa valeur à l'itération précédente, on arrive à une situation où le terminal considéré ne peut plus communiquer avec aucun satellite sans enfreindre ce seuil. Autrement dit, aucun satellite n'est compatible de la valeur du seuil 10 d'angle de séparation pour le terminal considéré. Pour prendre en compte cette possibilité, l'étape 411 est modifiée de la façon suivante. Pour chaque terminal de la liste identifiée à l'étape 410, on rétablit la valeur du seuil d'angle de séparation à la valeur fixée à l'itération précédente si au moins un satellite est accessible, pour ce terminal, avec 15 l'ancienne valeur. Dans le cas contraire, la valeur du seuil d'angle de séparation est maintenue et le niveau d'interférences global du système est abaissé en adaptant les caractéristiques RF des liaisons impliquant le terminal et un satellite du système. L'adaptation peut concerner la densité de puissance du 20 signal émis par le terminal (liaison montante) ou du signal émis par le satellite (liaison descendante). L'adaptation peut également concerner le type d'antenne du satellite et/ou la fréquence des signaux échangés et/ou la polarisation de ces signaux. L'étape 411 se termine dans tous les cas par le retrait du terminal de 25 la liste des terminaux candidats à l'abaissement de leur seuil d'angle de séparation. Dans un deuxième mode de réalisation, selon l'invention, de la méthode de détermination dynamique d'un seuil d'angle de séparation pour 30 chaque terminal, exécutée à l'étape 207, les liaisons montantes et les liaisons descendantes sont considérées indépendantes. Autrement dit, un 3045989 28 même terminal peut communiquer avec un premier satellite selon une liaison montante et avec un second satellite selon une liaison descendante. Dans ce deuxième exemple, on considère séparément un premier seuil d'angle de séparation pour les liaisons montantes et un second seuil d'angle de 5 séparation pour les liaisons descendantes. La méthode décrite à la figure 4 est alors appliquée indépendamment pour les liaisons montantes et pour les liaisons descendantes. On a représenté les étapes de la méthode selon le deuxième mode de réalisation appliqué aux liaisons montantes sur la figure 6 et appliqué aux liaisons 10 descendantes sur la figure 7. La méthode décrite à la figure 6 pour les liaisons unidirectionnelles montantes est identique à celle décrite à la figure 4 pour les liaisons bidirectionnelles à l'exception du fait que les étapes 405 et 409, qui concernent le calcul d'interférences générées par les liaisons descendantes, 15 sont supprimées. La méthode décrite à la figure 7 pour les liaisons unidirectionnelles descendantes est identique à celle décrite à la figure 4 pour les liaisons bidirectionnelles à l'exception du fait que les étapes 404 et 408, qui concernent le calcul d'interférences générées par les liaisons montantes, 20 sont supprimées. Les deux exemples de réalisation de la méthode, selon l'invention, de détermination dynamique d'un seuil d'angle de séparation tels que décrits à l'appui de l'algorithme de la figure 4 sont donnés à titre illustratifs. Sans sortir 25 du cadre de l'invention, la méthode peut être mise en oeuvre par le biais d'autres implémentations. En particulier, le traitement itératif décrit à la figure 4 peut être remplacé par un algorithme consistant à calculer les niveaux d'interférence globaux Epfd_up,Epfd_down pour plusieurs seuils d'angles de séparation, pris dans une plage de variation. Par exemple, on peut calculer 30 les niveaux d'interférences globaux pour des seuils d'angles de séparation variant de 2° à 7° par pas de 0,5°. L'ensemble des résultats de calcul est 3045989 29 sauvegardé en mémoire puis une optimisation paramétrique est réalisée pour déterminer, pour chaque terminal actif, le seuil d'angle de séparation qui permet d'obtenir le niveau d'interférences global le plus faible ou qui permet d'obtenir un niveau d'interférences global respectant le seuil réglementaire. 5 REVENDICATIONS1. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites dans un premier système de communications comprenant une constellation de satellites défilants (SAT_NGS01,SAT_NGS02) et un premier ensemble de terminaux (TER_NGS01,TER_NGS02), la méthode comprenant les étapes suivantes, pour chaque terminal (TER_NGS01,TER_NGS02) du premier ensemble : déterminer (207), un seuil d'angle de séparation distinct à respecter vis-à-vis d'une constellation de satellites géostationnaires (SAT_GSO), un angle de séparation étant défini comme l'angle topocentrique le plus faible sous lequel un satellite défilant (SAT_NGS01,SAT_NGS02) donné et un satellite géostationnaire (SAT_GS0) quelconque de l'arc géostationnaire sont vus depuis un terminal (TER_NGS01,TER_NGS02) donné du premier ensemble, autoriser (210) le terminal (TER_NGS01,TERNGS02) à établir une liaison de communication avec un satellite défilant (SAT_NGS01,SAT_NGS02) si l'angle de séparation (atil, at22) associé à ladite liaison est supérieur ou égal au seuil d'angle de séparation. 2. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon la revendication 1 dans laquelle le seuil d'angle de séparation distinct à respecter, pour chaque terminal (TER_NGS01,TER_NGS02) du premier ensemble est déterminé de manière à ce que l'ensemble des liaisons de communications autorisées (210) engendrent, sur un système de communications par satellites géostationnaires, un niveau de puissance interférente conforme à une limite réglementaire. 3045989 31 3. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon l'une des revendications précédentes dans laquelle la détermination (207) d'un seuil d'angle de séparation 5 comprend : le calcul (404,505) d'au moins un niveau de puissance interférente engendrée par l'ensemble des liaisons de communications entre les terminaux (TER_NGS01,TER_NGS02) du premier ensemble et les satellites défilants (SAT_NGS01,SAT_NGS02) vers un 10 système de communications par satellites géostationnaires, la recherche itérative (406,407,408,409,410,411,412), pour chaque terminal (TER_NGS01,TER_NGS02) du premier ensemble, de la valeur minimale du seuil d'angle de séparation permettant de maintenir ledit au moins un niveau de puissance interférente 15 calculé conforme à une limite réglementaire. 4. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon la revendication 3 dans laquelle la détermination (207) d'un seuil d'angle de séparation comprend une phase 20 (401,402) d'initialisation consistant à attribuer (401) une valeur initiale au seuil d'angle de séparation, pour chaque terminal, et à déclarer (402) chaque terminal candidat à l'abaissement de son seuil d'angle de séparation. 25 5. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon la revendication 4 dans laquelle ladite recherche itérative comprend, pour chaque terminal et à chaque itération, l'association (403) dudit terminal avec un satellite défilant respectant le seuil d'angle de séparation déterminé à l'itération courante et pour lequel 30 l'angle de séparation est le plus proche dudit seuil, les itérations étant 3045989 32 poursuivies tant qu'au moins un terminal est candidat à l'abaissement de son seuil d'angle de séparation. 6. Méthode d'allocation de ressources radio pour l'établissement de 5 communications par satellites selon la revendication 5 dans laquelle les liaisons de communication entre les terminaux (TER_NGS01,TER_NGS02) du premier ensemble et les satellites défilants (SAT_NGS01 ,SAT_NGS02) sont bidirectionnelles ou unidirectionnelles montantes et ladite recherche itérative comprend le 10 calcul (404) d'un premier niveau de puissance interférente engendrée par l'ensemble des liaisons de communications montantes vers l'arc géostationnaire pour un ensemble de positions supposées de satellites géostationnaires. 15 7. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon la revendication 6 dans laquelle le calcul (404) dudit premier niveau de puissance interférente est déterminé en choisissant, pour chaque satellite géostationnaire, une direction de pointage qui permet de maximiser la valeur dudit premier niveau de 20 puissance interférente. 8. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon l'une des revendications 6 ou 7 dans laquelle ladite recherche itérative comprend la recherche (408) d'au 25 moins un satellite géostationnaire pour lequel le premier niveau de puissance interférente calculé n'est pas conforme à la dite limite réglementaire. 9. Méthode d'allocation de ressources radio pour l'établissement de 30 communications par satellites selon la revendication 5 dans laquelle les liaisons de communications entre les terminaux (TER_NGS01,TER_NGS02) du premier ensemble et les satellites défilants (SAT_NGS01,SAT_NGS02) sont unidirectionnelles 3045989 33 descendantes et ladite recherche itérative comprend le calcul (405) d'un second niveau de puissance interférente engendrée par l'ensemble des liaisons de communication descendantes vers un ensemble de positions supposées de terminaux (TER_GSO) appartenant à un système de 5 communications par satellites géostationnaires. 10. Méthode d'allocation de ressources de communication par satellite selon l'une des revendications 6 à 8 dans laquelle les liaisons de communication entre les terminaux (TER_ NGS01,TER_NGS02) du 10 premier ensemble et les satellites défilants (SAT_NGS01,SAT_NGS02) sont bidirectionnelles et ladite recherche itérative comprend le calcul (405) d'un second niveau de puissance interférente engendrée par l'ensemble des liaisons de communication descendantes vers un ensemble de positions supposées de terminaux (TER_GSO) appartenant 15 à un système de communications par satellites géostationnaires. 11. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon l'une des revendications 9 ou 10 dans laquelle le calcul (405) dudit second niveau de puissance 20 interférente est déterminé en choisissant, pour chaque terminal (TER_GS0) appartenant à un système de communications par satellites géostationnaires, une direction de pointage qui permet de maximiser la valeur dudit second niveau de puissance interférente. 25 12.Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon l'une des revendications 9 à 11 dans laquelle ladite recherche itérative comprend la recherche (409) d'au moins un terminal appartenant à un système de communications par satellites géostationnaires pour lequel le second niveau de puissance 30 interférente calculé n'est pas conforme audit seuil réglementaire. 304 5 9 89 34 13. Méthode d'allocation de ressources radio pour l'établissement de 5 communications par satellites selon l'une des revendications 8 à 12 dans laquelle ladite recherche itérative comprend, pour chaque premier ou second niveau de puissance interférente calculé non conforme audit seuil réglementaire, la recherche (410) d'au moins un couple (terminal du premier ensemble, satellite défilant) dont la contribution, dans la valeur du 10 premier ou du second niveau de puissance interférente, est la plus élevée. 14.Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon la revendication 13 dans laquelle 15 ladite recherche itérative comprend, pour chaque terminal du couple (terminal du premier ensemble, satellite défilant) retenu, le rétablissement (411) du seuil d'angle de séparation déterminé à l'itération précédente et le retrait dudit terminal de la liste des candidats à l'abaissement du seuil d'angle de séparation. 20 15. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon la revendication 13 dans laquelle ladite recherche itérative comprend, pour au moins un terminal du couple 25 (terminal du premier ensemble, satellite défilant) retenu, l'adaptation des caractéristiques RF d'au moins une liaison impliquant ledit terminal et un satellite défilant, de sorte à diminuer le niveau de puissance interférente de l'ensemble des liaisons de communications autorisées (210) sur un système de communications par satellites géostationnaires. 30 16. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon l'une des revendications 6 à 15 dans laquelle la recherche itérative comprend une vérification (406) de la conformité du premier niveau de puissance interférente et/ou du second 3045989 niveau de puissance interférente et, si la conformité est vérifiée, l'abaissement (407) du seuil d'angle de séparation d'une valeur prédéterminée, pour chaque terminal candidat à l'abaissement de son seuil d'angle de séparation. 5 17. Méthode d'allocation de ressources radio pour l'établissement de communications par satellites selon l'une des revendications 4 à 16 dans laquelle la recherche itérative est arrêtée lorsqu'aucun terminal n'est plus candidat à l'abaissement de son seuil d'angle de séparation. 10 18. Dispositif d'allocation de ressources radio pour l'établissement de communications par satellites comprenant des moyens, incluant au moins un processeur et une mémoire, configurés pour mettre en oeuvre la 15 méthode d'allocation de ressources radio selon l'une des revendications précédentes. 19. Programme d'ordinateur comportant des instructions pour l'exécution de la méthode d'allocation de ressources radio pour l'établissement de 20 communications par satellites selon l'une quelconque des revendications 1 à 17, lorsque le programme est exécuté par un processeur. 20.Support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution de 25 la méthode d'allocation de ressources radio pour l'établissement des communications par satellites selon l'une quelconque des revendications 1 à 17, lorsque le programme est exécuté par un processeur.