CN114002717A - 一种基于投入使用风险的gso轨位选取方法 - Google Patents

一种基于投入使用风险的gso轨位选取方法 Download PDF

Info

Publication number
CN114002717A
CN114002717A CN202111278443.3A CN202111278443A CN114002717A CN 114002717 A CN114002717 A CN 114002717A CN 202111278443 A CN202111278443 A CN 202111278443A CN 114002717 A CN114002717 A CN 114002717A
Authority
CN
China
Prior art keywords
satellite
rail
candidate
frequency band
longitude
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111278443.3A
Other languages
English (en)
Inventor
张周
桑玮
胡向晖
王彤彤
许左宏
王利利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin (binhai) Intelligence Military-Civil Integration Innovation Center
Original Assignee
Tianjin (binhai) Intelligence Military-Civil Integration Innovation Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin (binhai) Intelligence Military-Civil Integration Innovation Center filed Critical Tianjin (binhai) Intelligence Military-Civil Integration Innovation Center
Priority to CN202111278443.3A priority Critical patent/CN114002717A/zh
Publication of CN114002717A publication Critical patent/CN114002717A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/37—Hardware or software details of the signal processing chain
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明公开了一种基于投入使用风险的GSO轨位选取方法，步骤如下：建立频段划分表，根据用户实际需求，选择一到多个频段类型；确定候选弧段，用户输入国家名称或在地图手绘区域，确定目标覆盖区域的经纬度；用户确定多个使用频段，并对每个频段分别进行如下操作：在目标弧段内，查找使用该频段且已登记入国家主表的轨位资源，记录到已使用资源表；在目标弧段内，查找使用该频段且已具有C资料的轨位资源，记录到未使用资源表；在目标弧段上，建立候选轨位表，依次对轨位进行判断并排序；根据排序规则，得到候选轨位序列，依据用户需求个数按顺序依次选取多个轨位。本发明详细的给出了轨位选取的依据，靠性强、操作简便，具有较强的工程实现性。

Description

一种基于投入使用风险的GSO轨位选取方法
技术领域
本发明涉及卫星频率轨道技术领域，特别是一种基于投入使用风险的GSO轨位选取方法。
背景技术
随着卫星技术的不断成熟以及人们对通信需求的爆炸增长，各国对频率轨道资源的抢占非常激烈，不仅申报迅速且申报数量较大，为了提升轨位选取效率，量化轨位选取技术点，并降低后期协调压力具有非常重要的意义。
现有的轨位选取来自于专家的经验分析，根据历史经验和官方信息判断周边轨位的使用情况以及后期协调难度，该方法完全依赖人工，难以大规模执行，不利于空间资源的储备。
发明内容
本发明的目的在于提供一种可靠性强、操作简便、易于实现的基于投入使用风险的GSO轨位选取方法。
实现本发明目的的技术解决方案为：一种基于投入使用风险的GSO轨位选取方，包括以下步骤：
步骤1、建立频段划分表，根据用户实际需求，选择一到多个频段类型；
步骤2、确定候选弧段，用户输入国家名称或在地图手绘区域，确定目标覆盖区域的经纬度；
步骤3、用户确定多个使用频段，并对每个频段分别进行步骤4～步骤8的操作；
步骤4、在目标弧段内，查找使用该频段且已登记入国家主表的轨位资源，记录到已使用资源表；
步骤5、在目标弧段内，查找使用该频段且已具有C资料的轨位资源，记录到未使用资源表；
步骤6、在目标弧段上，建立候选轨位表，依次对轨位进行判断并排序；
步骤7、根据排序规则，得到候选轨位序列，依据用户需求个数按顺序依次选取多个轨位。
本发明与现有技术相比，其显著优点为：(1)详细清楚的给出了轨位选取的依据，减轻了干扰计算辅助判断后期协调工作的难度；(2)利于空间资源的储备，靠性强、操作简便，能够大规模执行，具有较强的工程可实现性。
附图说明
图1是基于投入使用风险的GSO轨位选取方法的流程图。
图2是SNL查询页面示意图。
图3是SNL查询结果页面示意图。
图4是卫星网络详情示意图。
图5是BIU页面示意图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
结合图1，本发明基于投入使用风险的GSO轨位选取方法，包括以下步骤：
步骤1、建立频段划分表，根据用户实际需求，选择一到多个频段类型；
步骤2、确定候选弧段，用户输入国家名称或在地图手绘区域，确定目标覆盖区域的经纬度；
步骤3、用户确定多个使用频段，并对每个频段分别进行步骤4～步骤7的操作；
步骤4、在目标弧段内，查找使用该频段且已登记入国家主表的轨位资源，记录到已使用资源表；
步骤5、在目标弧段内，查找使用该频段且已具有C资料的轨位资源，记录到未使用资源表；
步骤6、在目标弧段上，建立候选轨位表，依次对轨位进行判断并排序；
步骤7、根据排序规则，得到候选轨位序列，依据用户需求个数按顺序依次选取多个轨位。
作为一种具体实施方式，步骤1中建立频段划分表，如表1所示。根据用户实际需求，选择其中的一到多个频段类型。
表1频段划分表
Figure BDA0003330417530000021
Figure BDA0003330417530000031
作为一种具体实施方式，步骤2中所述确定候选弧段，用户输入国家(地区)名称或在地图手绘区域，确定目标覆盖区域的经纬度，具体如下：
步骤2.1、判断该区域的经纬度是否能由一颗卫星完全覆盖；
①判断该区域的纬度是否在[-81°18′12″,81°18′12″]内，在则继续判断，不在则提示纬度过高，GSO卫星无法覆盖该纬度地区；
②判断该区域的经度跨度是否不大于2×81°18′12″，即当目标区域均在东经或西经时，满足|θ2-θ1|≤2×81°18′12″，当目标区域跨越0度经线时，满足|θ2+θ1|≤2×81°18′12″，当目标区域跨越180度经线时，满足360°-|θ2+θ1|≤2×81°18′12″，满足则判定能完全覆盖，否则提示需要多颗GSO卫星通信；
步骤2.2、判定该目标区域能由一颗GSO卫星覆盖后，确定该区域的最左侧经、纬度Left_long、Left_lati，以及最右侧经、纬度Right_long、Right_lati，卫星覆盖该区域时的最左、最右侧经度值通过计算获得，地球站天线俯仰角、地球站经度和卫星经度之间的计算关系如下：
Figure BDA0003330417530000032
其中，E表示地球站天线的俯仰角，α表示卫星与地球站的经度值差，B1表示地球站的纬度，r表示地球半径，R表示卫星到地心的距离；
设置地球站天线的俯仰角E为10°，r/R＝0.1511；将每个国家的纬度值代入上式得到α，arc_from＝Left_long+α，arc_to＝Right_long－α，由此得到卫星的目标弧段，即[arc_from，arc_to]。
作为一种具体实施方式，步骤3中确定多个使用频段，用户需具体到频段类型，系统进一步获取对应的上下行的值，用户选择其中的一到多个频段类型后，针对每个频段类型，都将进行步骤4～步骤7的操作。
作为一种具体实施方式，步骤4所述在目标弧段内，查找使用该频段且已登记入国家主表的轨位资源，记录到已使用资源表，具体如下：
在目标弧段内，查找使用该频段且已登记入国家主表的轨位资源，将对应的卫星网络名称和ID number记录到已使用资源表，已使用资源表包括4列，每一列分别为序号、唯一标识、卫星网络名称、国家，如表2所示，将ID number存入ntc_id字段，由于ID number可能存在多个，因此ID number以数组形式存储。
步骤4.1、查找已具有N资料的轨位资源；
如图2所示，在ITU提供的SNL网址通过SNL查询页面查找N资料，具体如下，在frequency行输入频段最低值和频段最高值，emission行选择All，longitude行输入经度范围，即目标弧段的最左侧和最右侧经度值，space or Earth行选择Geostationary，submission reason行选择notification，之后点击select。经过查找后，得到一个列表如图3所示，将列表中的卫星网络名称SATELLITE_NAME和国家ADM/ORG存到已使用空间资源表。之后依次点击每一行的卫星网络名称，得到列表如图4所示，将图4中的ID number列中不重复的值记录到已使用空间资源表的唯一标识中，以数组形式存放。
表2已使用空间资源表
序号 唯一标识 卫星网络名称 国家
步骤4.2、查找待投入使用的轨位资源；
如图5所示，在BIU页面输入相应值查找C资料，查找过程与步骤4.1类似，区别在于图2中submission reason行选择coordniation，得到如图3所示的列表后，记录卫星名称SATELLE NAME和国家ADM/ORG，将其代入BIU网址，如图5所示为该网址页面，在列表中查询是否具有相同卫星名称和国家，当查询到时，判断该条记录对应的Status是否为C，当为C时，点击该网络名称得到如图4所示内容，将图4中的ID number列中不重复的值记录到已使用空间资源表的唯一标识中，以数组形式存放。
作为一种具体实施方式，步骤5所述在目标弧段内，查找使用该频段且已具有C资料的轨位资源，将对应的卫星网络名称和ID number记录到已占用空间资源表。与步骤4.2相同，区别在于在BIU网址中，选择Status为N的资料存储。
表3已占用资源表
序号 唯一标识 卫星网络名称 国家
作为一种具体实施方式，步骤6所述在目标弧段上，建立候选轨位表，如表4所示，依次对轨位进行判断并排序，其中候选轨位表包括4列，每一列分别为候选轨位、周边已占用轨位间隔、ΔT/T、待协调数量，ΔT表示由于干扰造成的受扰卫星等效噪声温度增加量，T表示等效卫星链路噪声温度，依次进行如下操作：
步骤6.1、删除已使用资源表和未使用资源表的轨位；
步骤6.2、以目标弧段最左侧为初始值，目标弧段最右侧为终点值，以0.2度为间隔，依次假定当前轨位为候选轨位，将候选轨位与自身周边3度范围内已使用和已占用轨位进行干扰计算，将计算结果代入表4候选轨位表的ΔT/T列，同时计算候选轨位和与之进行干扰计算轨位的间隔，将当前轨位与周边已占用轨位依次计算，并以数组形式存储；周边轨位均计算完毕后，筛选出对应ΔT/T数组中的最大值，当最大值不大于6％时，令待协调数量为0，否则令待协调数量为ΔT/T数组中不小于6％的数量；
表4候选轨位表
候选轨位 周边已占用轨位间隔 ΔT/T 待协调数量
进一步地，所述将候选轨位与自身周边3度范围内已使用和已占用轨位进行干扰计算，具体如下：
依次读取已使用资源表和已占用资源表，读取对应的ntc_id值，在srs数据库中查找该ntc_id的资料，之后进行如下操作：
(1)判断干扰和受干扰系统方向是否相同：
当相同时，计算由于干扰造成的受扰卫星等效噪声温度增加量ΔT，令ΔT＝ΔTs+ΔTe，其中ΔTs表示干扰发射下卫星S接收天线的输出，
Figure BDA0003330417530000061
ΔTe表示干扰发射下地球站接收天线的输出，
Figure BDA0003330417530000062
具体地，pe′表示传输到发射地球站的每Hz最大功率密度，g1′(θt)表示地球站在卫星S′方向的发射天线增益，g2(δe′)表示从卫星S到发射地球站方向上的接收天线增益，k表示玻尔兹曼常数1.38×10-23J/K，lu表示地球站到卫星S的上行链路自由空间损耗，g3(ηe)表示在ηe方向上卫星S的传输天线增益，g4(θt)表示在卫星S′方向上地球站的接收天线增益，ld表示卫星S到地球站的下行链路自由空间传输损耗，ps′表示传输到卫星天线的每Hz最大功率密度；
当相反时，进行地面站和地面站之间的协调处理；
Figure BDA0003330417530000063
ΔT＝γΔTs，卫星干扰辐射引起的等效噪声温度的增加
Figure BDA0003330417530000064
其中，g3′(ηs)表示在ηs方向上卫星S的传输天线增益，g2(δs′)表示从卫星S到卫星S′方向上的接收天线增益，ls表示从卫星S′到卫星S的自由空间损耗，γ表示从卫星S接收天线的输出到地球站接收天线的输出的传输增益，γ′表示从卫星S′接收天线的输出到地球站接收天线的输出的传输增益，ΔTs′表示表示干扰发射下卫星S′接收天线的输出，ps表示传输到卫星S天线的每Hz最大功率密度，g3′(ηs′)表示在ηs方向上卫星S′的传输天线增益，g2(δs)表示从卫星S′到卫星S方向上的接收天线增益；。
(2)当网络的主管部门通知极化情况时，可以进一步考虑极化隔离系数，则当干扰和受干扰系统方向相同时，
Figure BDA0003330417530000065
方向相反时，
Figure BDA0003330417530000066
其中，Yu、Yd、Yss均由表5给出，Yd表示下行链路上的隔离系数，Yu表示上行链路的隔离系数，Yss表示星间链路的隔离系数；
表5极化隔离系数表
Figure BDA0003330417530000071
(3)计算ΔT/T，当ΔT/T＞6％时，两个网络之间需要协调，否则不需要协调。
步骤6.3、将候选轨位分为两类，第一类是待协调数量为0的候选轨位，第二类是剩余候选轨位，令第一类候选轨位的优先级均高于第二类候选轨位的优先级；
将第一类轨位进行如下判断：
用户确定候选频段，即潜在可能开发的频段，之后将该频段代入步骤4中，查找该频段周边的轨位占有情况，判断该频段周边轨位的占有情况，记录资料总数，将候选轨位按照资料总数，由小到大进行排列，作为第一类轨位的排序。
将第二类轨位进行如下判断：
将轨位按照周边3度范围内待协调传统运营商轨位的数量按从小到大进行排序，传统运营商轨位如表6所示。当存在数量相同的情况时，比较与传统运营商轨位的ΔT/T最大值，该值较小者优先级更高。
表6国际传统运营商轨位表
operator sat_name status lifetime frequency orbit risk_interval_deg
INMARSAT Inmarsat GX5 moving 16+ Ka 11°E 1.5
measat Measat 3 active 15 C,Ku 91.5°E
Skynet Skynet 5B active 15+ X,UHF 24.5°E
Eutelsat Eutelsat 3B active 15 Ku,C,Ka 3°E 1.5
MUOS MUOS-1 active 15 Ka 177°W
WGS WGS 10 active 14 X,Ka 53°W 1.5
ViaSat Inc ViaSat 1 active 12 Ka 115°W 1.5
Intelsat Intelsat 24 retired 11 Ku 31E
JSAT JCSAT-2B active 15 C,Ku 154°E
作为一种具体实施方式，步骤7、根据步骤7中的排序规则，得到候选轨位序列，根据用户需求个数按顺序依次选取多个轨位。
实施例
步骤1、建立频段划分表，如表1.1所示。
表1.1频段划分表
Figure BDA0003330417530000081
步骤2、确定候选弧段。用户选择覆盖地区为中国。
1、经过查找得到，中国的纬度在[-81°18′12″,81°18′12″]内，并且经度跨度不大于2×81°18′12″，因此可以由一颗GSO卫星完全覆盖。
2、根据目标区域的最左侧和最右侧经度，计算目标弧段，计算公式如下：
Figure BDA0003330417530000082
其中，E表示地球站天线的俯仰角，α表示卫星与地球站的经度值差，B1表示地球站的纬度，r表示地球半径，R表示卫星到地心的距离；
设置地球站天线的俯仰角E为10°，r/R＝0.1511；将目标区域的纬度值代入上式得到α，arc_from＝Left_long+α，arc_to＝Right_long－α。由此得到卫星的目标弧段，即[73.39，139.1]。
步骤3、确定多个使用频段，用户根据实际需求，选择频段类型为标准C和标准Ku，即3700-4200MHz、5925-6425MHz、12.2-12.75GHz、14.0-14.45GHz。将标准C和标准Ku频段分别进行步骤3到步骤6的计算。
步骤4、在目标弧段内，查找使用标准C和标准Ku且已登记入国家主表的轨位资源和待投入使用的轨位资源，将对应的卫星网络名称和ID number记录到已使用资源表，如表1.2和表1.3所示，由于数据项较多，该表只展示查询到的部分内容。
表1.2标准C频段已使用空间资源表
序号 唯一标识 卫星网络名称 国家
1 90504231 INSAT-2(74) IND
2 101540108 APSTAR-76E CHN
3 94540047 THAICOM-A2B THA
表1.3标准Ku频段已使用空间资源表
序号 唯一标识 卫星网络名称 国家
1 93540116 INTERBELAR-2 RUS/IK
2 98543166 INTERSPUTNIK-75E-Q RUS/IK
3 101540108 APSTAR-76E CHN
步骤5、在目标弧段内，查找使用使用标准C和标准Ku且已具有C资料的轨位资源，将对应的卫星网络名称和ID number记录到已占用空间资源表，如表1.4和表1.5所示。
表1.4标准C频段已占用资源表
序号 唯一标识 卫星网络名称 国家
1 120520041 IK-74.5E RUS/IK
2 119520006 HISPASAT-92A E
3 100543885 COMPASS-80E CHN
表1.5标准Ku频段已占用资源表
序号 唯一标识 卫星网络名称 国家
1 114540193 GSAT-NS IND
2 118520285 QATARSAT-G5-77E QAT
3 120520146 MADAR-77E UAE
步骤6、在目标弧段上，建立候选轨位表，如表1.6所示，构建过程如下：
1、删除已使用资源表和未使用资源表的轨位；
2、以目标弧段最左侧为初始值，目标弧段最右侧为终点值，以0.2度为间隔，依次假定当前轨位为目标轨位，与周边3度范围内已使用和已占用轨位进行干扰计算，将计算结果代入表4的ΔT/T列，同时计算二者间隔，将当前轨位与周边已占用轨位依次计算，并以数组形式存储。周边轨位均计算完毕后，筛选出对应ΔT/T数组中的最大值，当最大值不大于6％时，令待协调数量为0，否则令待协调数量为ΔT/T数组中不小于6％的数量。
表1.6候选轨位表
候选轨位 周边已占用轨位间隔 ΔT/T 待协调数量
96.59 [0.39,0.41…,3.41] [0,6％,…,2％] 72
96.79 [0.21,0.21,…,3.29] [4％,5％,…,3％] 50
86.79 [0.21,0.21,…,3.29] [3％,5％,…,0] 66
86.19 [0.19,0.19,…,3.31] [6％,1％,…,1％] 45
96.39 [0.19,0.29,…,3.39] [3％,8％,…,3％] 79
86.39 [0.11,0.11,…,3.41] [2％,6％,…,3％] 31
114.39 [0.11,0.11,…,3.41] [2％,8％,…,4％] 72
86.59 [0.09,0.09,…,3.41] [9％,7％,…,4％] 61
114.19 [0.09,0.19,…,3.31] [7％,3％,…,5％] 67
96.99 [0.01,…,0.79,…,3.49] [3％,4％,…,2％] 21
86.99 [0.01,0.01,…,3.49] [2％,3％,…,5％] 30
96.19 [0.01,0.09,…,3.31] [7％,4％,…,6％] 75
其中，干扰计算过程如下：
依次读取已使用资源表和已占用资源表，读取对应的ntc_id值，在srs数据库中查找该ntc_id的资料，之后通过计算进行如下操作：
(1)判断干扰和受干扰系统方向是否相同。
当相同时，计算由于干扰造成的受扰卫星等效噪声温度增加量ΔT，令ΔT＝ΔTs+ΔTe，其中ΔTs表示干扰发射下卫星接收天线的输出，
Figure BDA0003330417530000101
ΔTe表示干扰发射下地球站接收天线的输出，
Figure BDA0003330417530000111
具体地，pe′表示传输到发射地球站的每Hz最大功率密度，g1′(θt)表示地球站在卫星S′方向的发射天线增益，g2(δe′)表示从卫星S到发射地球站方向上的接收天线增益，k表示玻尔兹曼常数(1.38×10-23J/K)，lu表示地球站到卫星S的上行链路自由空间损耗，ld表示卫星S到地球站的下行链路自由空间传输损耗，ps′表示传输到卫星天线的每Hz最大功率密度。
当相反时，将通过类似于用于地面站和地面站之间协调的协调程序来处理，与地心角θg有关。卫星接收天线输出的噪声温度增量
Figure BDA0003330417530000112
等效链路噪声温度的表观增量为ΔT＝γΔTs，卫星干扰辐射引起的等效噪声温度的增加
Figure BDA0003330417530000113
(2)当网络的主管部门通知极化情况时，可以进一步考虑极化隔离系数，则当干扰和受干扰系统方向相同时，
Figure BDA0003330417530000114
方向相反时，
Figure BDA0003330417530000115
其中，Yu、Yd、Yss均由下表给出，Yd表示下行链路上的隔离系数，Yu表示上行链路的隔离系数，Yss表示星间链路的隔离系数，如表1.7所示。
表1.7极化隔离系数表
Figure BDA0003330417530000116
(3)计算ΔT/T，当ΔT/T＞6％时，证明两个网络之间需要协调，否则不需要协调。
3、将候选轨位分为两类，第一类是待协调数量为0的候选轨位，第二类是剩余候选轨位，令第一类候选轨位的优先级均高于第二类候选轨位的优先级。
将第一类轨位进行如下判断：
用户确定候选频段，即潜在可能开发的频段，之后将该频段代入步骤4中，查找该频段周边的轨位占有情况，判断该频段周边轨位的占有情况，记录资料总数，将候选轨位按照资料总数，由小到大进行排列，作为第一类轨位的排序。
将第二类轨位进行如下判断：
将轨位按照周边3度范围内待协调传统运营商轨位的数量按从小到大进行排序，传统运营商轨位如表1.8所示。当存在数量相同的情况时，比较与传统运营商轨位的ΔT/T最大值，该值较小者优先级更高。
表1.8国际传统运营商轨位表
operator sat_name status lifetime frequency orbit risk_interval_deg
INMARSAT Inmarsat GX5 moving 16+ Ka 11°E 1.5
measat Measat 3 active 15 C,Ku 91.5°E
Skynet Skynet 5B active 15+ X,UHF 24.5°E
Eutelsat Eutelsat 3B active 15 Ku,C,Ka 3°E 1.5
MUOS MUOS-1 active 15 Ka 177°W
WGS WGS 10 active 14 X,Ka 53°W 1.5
ViaSat Inc ViaSat 1 active 12 Ka 115°W 1.5
Intelsat Intelsat 24 retired 11 Ku 31E
JSAT JCSAT-2B active 15 C,Ku 154°E
经排序后，得到的轨位顺序为：
[86.79,86.19,96.59,114.19,96.99,114.39,86.39,86.59,96.79,96.19,96.39,86.99,86.39]。
步骤7、根据步骤6中的排序规则，得到候选轨位序列，根据用户需求个数按顺序依次选取多个轨位。
综上所述，本发明详细清楚的给出了轨位选取的依据，减轻了干扰计算辅助判断后期协调工作的难度；利于空间资源的储备，靠性强、操作简便，能够大规模执行，具有较强的工程可实现性。

Claims (7)

1.一种基于投入使用风险的GSO轨位选取方法，其特征在于，包括以下步骤：
步骤1、建立频段划分表，根据用户实际需求，选择一到多个频段类型；
步骤2、确定候选弧段，用户输入国家名称或在地图手绘区域，确定目标覆盖区域的经纬度；
步骤3、用户确定多个使用频段，并对每个频段分别进行步骤4～步骤7的操作；
步骤4、在目标弧段内，查找使用该频段且已登记入国家主表的轨位资源，记录到已使用资源表；
步骤5、在目标弧段内，查找使用该频段且已具有C资料的轨位资源，记录到未使用资源表；
步骤6、在目标弧段上，建立候选轨位表，依次对轨位进行判断并排序；
步骤7、根据排序规则，得到候选轨位序列，依据用户需求个数按顺序依次选取多个轨位。
2.根据权利要求1所述的基于投入使用风险的GSO轨位选取方法，其特征在于，步骤1所述建立频段划分表，具体如表1所示：
表1频段划分表
Figure FDA0003330417520000011
Figure FDA0003330417520000021
3.根据权利要求1所述的基于投入使用风险的GSO轨位选取方法，其特征在于，步骤2所述确定候选弧段，用户输入国家名称或在地图手绘区域，确定目标覆盖区域的经纬度，具体如下：
步骤2.1、判断该区域的经纬度是否能由一颗卫星完全覆盖；
①判断该区域的纬度是否在[-81°18′12″,81°18′12″]内，在则继续判断，不在则提示纬度过高，GSO卫星无法覆盖该纬度地区；
②判断该区域的经度跨度是否不大于2×81°18′12″，即当目标区域均在东经或西经时，满足|θ2-θ1|≤2×81°18′12″，当目标区域跨越0度经线时，满足|θ2+θ1|≤2×81°18′12″，当目标区域跨越180度经线时，满足360°-|θ2+θ1|≤2×81°18′12″，满足则判定能完全覆盖，否则提示需要多颗GSO卫星通信；
步骤2.2、判定该目标区域能由一颗GSO卫星覆盖后，确定该区域的最左侧经、纬度Left_long、Left_lati，以及最右侧经、纬度Right_long、Right_lati，卫星覆盖该区域时的最左、最右侧经度值通过计算获得，地球站天线俯仰角、地球站经度和卫星经度之间的计算关系如下：
Figure FDA0003330417520000022
其中，E表示地球站天线的俯仰角，α表示卫星与地球站的经度值差，B1表示地球站的纬度，r表示地球半径，R表示卫星到地心的距离；
设置地球站天线的俯仰角E为10°，r/R＝0.1511；将每个国家的纬度值代入上式得到α，arc_from＝Left_long+α，arc_to＝Right_long－α，由此得到卫星的目标弧段，即[arc_from，arc_to]。
4.根据权利要求1所述的基于投入使用风险的GSO轨位选取方法，其特征在于，步骤4所述在目标弧段内，查找使用该频段且已登记入国家主表的轨位资源，记录到已使用资源表，具体如下：
在目标弧段内，查找使用该频段且已登记入国家主表的轨位资源，将对应的卫星网络名称和ID number记录到已使用资源表，已使用资源表包括4列，每一列分别为序号、唯一标识、卫星网络名称、国家，将ID number存入ntc_id字段，由于ID number可能存在多个，因此ID number以数组形式存储；
步骤4.1、查找已具有N资料的轨位资源；
在ITU提供的SNL网址通过SNL查询页面查找N资料，具体如下：在frequency行输入频段最低值和频段最高值，emission行选择All，longitude行输入经度范围，即目标弧段的最左侧和最右侧经度值，space or Earth行选择Geostationary，submission reason行选择notification，之后点击select；经过查找后，得到一个SNL查询结果列表，将列表中的卫星网络名称SATELLITE_NAME和国家ADM/ORG存到已使用空间资源表；之后依次点击每一行的卫星网络名称，得到卫星网络详情列表，将该列表中的ID number列中不重复的值记录到已使用空间资源表的唯一标识中，以数组形式存放；
步骤4.2、查找待投入使用的轨位资源；
在BIU页面输入相应值查找C资料，查找过程与步骤4.1区别在于，SNL查询页面submission reason行选择coordniation，得到SNL查询结果列表后，记录卫星名称SATELLENAME和国家ADM/ORG，代入BIU网址，在BIU页面列表中查询是否具有相同卫星名称和国家，当查询到时，判断该条记录对应的Status是否为C，当为C时，点击该网络名称得到卫星网络详情，将卫星网络详情的ID number列中不重复的值记录到已使用空间资源表的唯一标识中，以数组形式存放。
5.根据权利要求4所述的基于投入使用风险的GSO轨位选取方法，其特征在于，步骤5中在目标弧段内，查找使用该频段且已具有C资料的轨位资源，将对应的卫星网络名称和IDnumber记录到已占用空间资源表，与步骤4.2相同，区别在于在BIU网址中，选择Status为N的资料存储。
6.根据权利要求1所述的基于投入使用风险的GSO轨位选取方法，其特征在于，步骤6所述在目标弧段上，建立候选轨位表，依次对轨位进行判断并排序，其中候选轨位表包括4列，每一列分别为候选轨位、周边已占用轨位间隔、ΔT/T、待协调数量，ΔT表示由于干扰造成的受扰卫星等效噪声温度增加量，T表示等效卫星链路噪声温度，依次进行如下操作：
步骤6.1、删除已使用资源表和未使用资源表的轨位；
步骤6.2、以目标弧段最左侧为初始值，目标弧段最右侧为终点值，以0.2度为间隔，依次假定当前轨位为候选轨位，将候选轨位与自身周边3度范围内已使用和已占用轨位进行干扰计算，将计算结果代入候选轨位表的ΔT/T列，同时计算候选轨位和与之进行干扰计算轨位的间隔，将当前轨位与周边已占用轨位依次计算，并以数组形式存储；周边轨位均计算完毕后，筛选出对应ΔT/T数组中的最大值，当最大值不大于6％时，令待协调数量为0，否则令待协调数量为ΔT/T数组中不小于6％的数量；
步骤6.3、将候选轨位分为两类，第一类是待协调数量为0的候选轨位，第二类是剩余候选轨位，令第一类候选轨位的优先级均高于第二类候选轨位的优先级；
将第一类轨位进行如下判断：
用户确定候选频段，即潜在可能开发的频段，之后将该频段代入步骤4中，查找该频段周边的轨位占有情况，判断该频段周边轨位的占有情况，记录资料总数，将候选轨位按照资料总数，由小到大进行排列，作为第一类轨位的排序；
将第二类轨位进行如下判断：
将轨位按照周边3度范围内待协调传统运营商轨位的数量按从小到大进行排序，当存在数量相同的情况时，比较与传统运营商轨位的ΔT/T最大值，该值较小者优先级更高。
7.根据权利要求6所述的基于投入使用风险的GSO轨位选取方法，其特征在于，步骤6.2中所述将候选轨位与自身周边3度范围内已使用和已占用轨位进行干扰计算，具体如下：
(1)判断干扰和受干扰系统方向是否相同：
当相同时，计算由于干扰造成的受扰卫星等效噪声温度增加量ΔT，令ΔT＝ΔTs+ΔTe，其中ΔTs表示干扰发射下卫星S接收天线的输出，
Figure FDA0003330417520000041
ΔTe表示干扰发射下地球站接收天线的输出，
Figure FDA0003330417520000042
具体地，p′e表示传输到发射地球站的每Hz最大功率密度，g′1(θt)表示地球站在卫星S′方向的发射天线增益，g2(δe′)表示从卫星S到发射地球站方向上的接收天线增益，k表示玻尔兹曼常数1.38×10-23J/K，lu表示地球站到卫星S的上行链路自由空间损耗，g3(ηe)表示在ηe方向上卫星S的传输天线增益，g4(θt)表示在卫星S′方向上地球站的接收天线增益，ld表示卫星S到地球站的下行链路自由空间传输损耗，p′s表示传输到卫星天线的每Hz最大功率密度；
当相反时，进行地面站和地面站之间的协调处理；
Figure FDA0003330417520000051
ΔT＝γΔTs，卫星干扰辐射引起的等效噪声温度的增加
Figure FDA0003330417520000052
其中，g′3(ηs)表示在ηs方向上卫星S的传输天线增益，g2(δs′)表示从卫星S到卫星S′方向上的接收天线增益，ls表示从卫星S′到卫星S的自由空间损耗，γ表示从卫星S接收天线的输出到地球站接收天线的输出的传输增益，γ′表示从卫星S′接收天线的输出到地球站接收天线的输出的传输增益，ΔTs′表示表示干扰发射下卫星S′接收天线的输出，ps表示传输到卫星S天线的每Hz最大功率密度，g′3(ηs′)表示在ηs方向上卫星S′的传输天线增益，g2(δs)表示从卫星S′到卫星S方向上的接收天线增益；
(2)当网络的主管部门通知极化情况时，进一步考虑极化隔离系数，则当干扰和受干扰系统方向相同时
Figure FDA0003330417520000053
方向相反时
Figure FDA0003330417520000054
其中Yd表示下行链路上的隔离系数，Yu表示上行链路的隔离系数，Yss表示星间链路的隔离系数；
(3)计算ΔT/T，当ΔT/T＞6％时，两个网络之间需要协调，否则不需要协调。
CN202111278443.3A 2021-10-30 2021-10-30 一种基于投入使用风险的gso轨位选取方法 Pending CN114002717A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111278443.3A CN114002717A (zh) 2021-10-30 2021-10-30 一种基于投入使用风险的gso轨位选取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111278443.3A CN114002717A (zh) 2021-10-30 2021-10-30 一种基于投入使用风险的gso轨位选取方法

Publications (1)

Publication Number Publication Date
CN114002717A true CN114002717A (zh) 2022-02-01

Family

ID=79925749

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111278443.3A Pending CN114002717A (zh) 2021-10-30 2021-10-30 一种基于投入使用风险的gso轨位选取方法

Country Status (1)

Country Link
CN (1) CN114002717A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114826370A (zh) * 2022-03-16 2022-07-29 军事科学院系统工程研究院系统总体研究所 一种gso卫星的轨位选取方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020146982A1 (en) * 2001-04-04 2002-10-10 Mclain Christopher John Method and apparatus using transmit beam lobing for identifying an interfering mobile terminal
CN102523032A (zh) * 2011-12-08 2012-06-27 中国空间技术研究院 一种地球静止轨道卫星轨位确定方法
CN102520461A (zh) * 2011-12-08 2012-06-27 中国空间技术研究院 Ngso卫星地球探测系统对深空探测系统的干扰确定方法
US20150077288A1 (en) * 2013-09-17 2015-03-19 Accord Software & Systems Pvt Ltd Time To First Fix Optimization In A Satellite Navigation Receiver
US20150268349A1 (en) * 2014-03-21 2015-09-24 The Boeing Company Interference geolocation using a satellite constellation
WO2017198325A1 (de) * 2016-05-17 2017-11-23 Rheinmetall Air Defence Ag Antennenanordnung eines lenkflugkörpers mit einer radarantenne
CN112882071A (zh) * 2021-01-11 2021-06-01 北京新星宇航科技有限公司 一种静止轨道卫星网络轨位自动化选取方法
CN113162676A (zh) * 2021-03-26 2021-07-23 天津(滨海)人工智能军民融合创新中心 基于轨位多级联合风险的gso轨位效能评估方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020146982A1 (en) * 2001-04-04 2002-10-10 Mclain Christopher John Method and apparatus using transmit beam lobing for identifying an interfering mobile terminal
CN102523032A (zh) * 2011-12-08 2012-06-27 中国空间技术研究院 一种地球静止轨道卫星轨位确定方法
CN102520461A (zh) * 2011-12-08 2012-06-27 中国空间技术研究院 Ngso卫星地球探测系统对深空探测系统的干扰确定方法
US20150077288A1 (en) * 2013-09-17 2015-03-19 Accord Software & Systems Pvt Ltd Time To First Fix Optimization In A Satellite Navigation Receiver
US20150268349A1 (en) * 2014-03-21 2015-09-24 The Boeing Company Interference geolocation using a satellite constellation
WO2017198325A1 (de) * 2016-05-17 2017-11-23 Rheinmetall Air Defence Ag Antennenanordnung eines lenkflugkörpers mit einer radarantenne
CN112882071A (zh) * 2021-01-11 2021-06-01 北京新星宇航科技有限公司 一种静止轨道卫星网络轨位自动化选取方法
CN113162676A (zh) * 2021-03-26 2021-07-23 天津(滨海)人工智能军民融合创新中心 基于轨位多级联合风险的gso轨位效能评估方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈志鹏 等: "基于"一带一路"区域服务性能提升的北斗卫星导航系统组网方案分析", 《第九届中国卫星导航学术年会论文集——S04 卫星轨道与钟差》, 23 May 2018 (2018-05-23), pages 51 - 58 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114826370A (zh) * 2022-03-16 2022-07-29 军事科学院系统工程研究院系统总体研究所 一种gso卫星的轨位选取方法

Similar Documents

Publication Publication Date Title
CN112882071B (zh) 一种静止轨道卫星网络轨位自动化选取方法
Vook et al. 5g new radio: Overview and performance
CN102523032B (zh) 一种地球静止轨道卫星轨位确定方法
CN107305246B (zh) 基于接收信号强度指示的定位方法和装置
CN113114397B (zh) 卫星移动通信系统天地同频共用集总干扰估计方法
CN114002717A (zh) 一种基于投入使用风险的gso轨位选取方法
CN109239740B (zh) 数据发送方法、装置、数据接收方法及装置
CN113627799B (zh) 基于军事效益的gso频轨资源效能评估方法
CN113630172A (zh) 基于空间资源数据库与轮廓线插值的上行干扰分析方法
CN114399652B (zh) 基于雷达图智能匹配的多指标体系geo卫星轨位选取方法
CN113162676B (zh) 基于轨位多级联合风险的gso轨位效能评估方法
Dadson et al. Computer prediction of field strength in the planning of radio systems
Verspieren et al. Innovative approach to data gathering in remote areas using constellations of store and forward communication cubesats
CN115967620A (zh) 一种自动实现地面站网设备定向匹配的方法与装置
CN113610363B (zh) 基于轨位投入使用风险的gso频轨资源效能评估方法
CN111491255B (zh) 一种回填经纬度定位扩量方法、系统及存储介质
Peña et al. Radiopropagation Simulations in 28 GHz, 38 GHz and 60 GHz with 128 Elements Massive MIMO Array for 5G Networks
CN114630362B (zh) 空间频谱资源对卫星通信系统下行链路的支持度计算方法
Shi et al. The interference analysis of satellite services in the 8025–8400 MHz frequency band
CN116248163B (zh) 一种面向大规模低轨星座的频率兼容性分析方法
CN113296743A (zh) 基于GSO C资料申报的SpaceCap软件映射方法
CN114665950B (zh) 空间频谱资源对卫星通信系统上行链路的支持度计算方法
Williamson Broadcasting by Satellite: Some Technical Considerations
Adetona et al. Television broadcasting transmission in Nigeria: application of technical data for re-planning
Saha Design of pan-india rural broadband network planning tool

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination