GNSS无线电掩星探测能提供全球覆盖的地球电离层^[1-3]和中性大气^[4-8]剖面，同时具有全天候、长期稳定、高精度和高垂直分辨率等优点，对数值天气预报^[9]和气候^[10]研究具有重要作用，在天文、气象、空间以及国防领域具有广泛的应用前景^[11].最初的GNSS掩星探测计划是1995年的GPS/MET计划^[12]，该计划由美国实施，在低轨卫星上安装了一台GPS掩星探测器，对GNSS掩星测量地球大气的理论进行了验证.GPS/MET在轨运行了两年时间，获得了很多宝贵的观测资料.科学家们利用GPS/MET提供的掩星数据，成功地反演得到了中性大气0~60km高度的大气参数剖面^[13-15].结果证明了GNSS无线电掩星测量地球大气理论的可行性，从此以后，世界各国开展了许多GNSS无线电掩星探测项目^[16]，掩星探测地球大气技术取得了极大进展.

美国的GPS系统是至今为止最完善也最稳定的GNSS系统，因此，在已有的掩星探测项目中，使用的GNSS信号来源主要为GPS卫星信号.然而，近年来不少国家和地区也在积极研制自己的卫星导航系统.俄罗斯从1978年10月开始筹备发射自己的全球导航卫星系统（GLONASS），该系统由21颗工作卫星和3颗在轨备用卫星组成，均匀分布在3个轨道平面上.欧空局（ESA）也在筹建民用导航卫星系统（GALILEO），到2011年底已经发射了4颗实验卫星.然而，GLONASS和GALILEO系统的建设都由于资金问题有所延迟.中国的北斗卫星导航系统（COMPASS）是我国正在实施的自主研发、独立运行的全球卫星导航系统^[17]，具有开放性、自主性、兼容性、渐进性的特点，其建设按照“先区域、后全球”的总体思路分步实施，在2012年建立包含14颗北斗卫星的中国及周边区域导航系统，并将在2020年全面建成全球覆盖的北斗卫星导航系统^[18].截止2011年底，我国已成功发射了十余颗北斗导航卫星，基本完成了中国区域的覆盖.

GNSS掩星事件的数量及切点分布与GNSS系统的空间卫星星座有着密切关系.在不同的GNSS卫星星座下，需要进行大量的仿真以选择合适的LEO轨道参数，用以获取尽量多的掩星数量和全球覆盖并均匀分布的掩星切点.国内外学者在这方面的工作针对GPS卫星星座的较多.赵世军等^[19]在一定的大气折射指数分布模型和现有GPS星座的假设条件下，利用射线追踪法模拟了掩星事件，并讨论了LEO卫星的轨道倾角和高度对掩星事件发生的数量和分布的影响.徐晓华等^[20]分析了对于单颗LEO卫星，GPS掩星事件的分布和数量随着LEO轨道参数包括轨道升交角距、升交点赤经、轨道高度和倾角而变化的规律，并采用数值模拟方法，对不同卫星数和不同星座参数的玫瑰型星座在一天内所观测的掩星事件的数量和分布进行了比较分析^[21].杜晓勇等^[22]通过仿真计算，定量地讨论了轨道参数对LEO-LEO掩星事件数量及分布的影响.

与GPS卫星星座由6个轨道的MEO卫星所构成不同，中国的北斗卫星导航系统空间星座由地球静止轨道^[23]（GEO）、倾斜地球同步轨道（IGSO）和中圆地球轨道（MEO）三种轨道的卫星组成，并且北斗的MEO卫星轨道分布与GPS的MEO卫星轨道分布有所不同.现阶段北斗导航系统已基本完成由14颗星组成的区域导航系统的建设，在中国区域内具有较好的可见性和定位精度以及良好的可用性^[24].与GPS相比，北斗系统可以方便的提供双频或三频载波相位^[25]，有利于掩星反演.如何高效地使用北斗卫星星座进行GNSS掩星探测值得深入研究.国内外文献中还未发现描述利用北斗卫星星座，特别是利用北斗GEO与IGSO卫星进行掩星探测的研究.

本文对全球范围内的北斗掩星事件进行了模拟，并在不同的LEO轨道高度、倾角、升交点赤经、近地点角距条件下，对北斗掩星事件的数量和分布进行了仿真分析.仿真中针对北斗GEO、IGSO和MEO卫星特点，对其掩星事件特性分别进行了研究.结论对利用北斗卫星导航系统进行掩星探测具有参考价值.

本文中采用的北斗导航卫星星座构成参考文献[26]，文献中采用了总共35颗卫星的北斗星座配置，包括5颗GEO、3颗IGSO、27颗MEO.考虑到最新的北斗接口文件^[17]中包括5颗倾斜地球同步轨道（IGSO）卫星，因此在本文涉及的仿真中将北斗卫星总数扩展为37，包括5颗地球静止轨道（GEO）、5颗倾斜地球同步轨道（IGSO）卫星和27颗中圆地球轨道（MEO）卫星.5颗北斗GEO卫星的轨道位置分别为58.75°E、80°E、110.5°E、140°E和160°E.5颗北斗IGSO卫星的倾角为55°，交叉点经度为118°E，同一时刻的平近点角间隔72°.27颗北斗MEO卫星分布于升交点赤经分别为0°、120°、240°，倾角为55°的三个轨道平面，每个轨道上的第一颗卫星在仿真起始时刻的平近点角分别为0°、15°、30°，其余卫星平近点角依次增加45°.此外，3个轨道面上备用卫星的平近点角分别为10°、55°和105°^[26].北斗的MEO卫星倾角与GPS倾角接近，星下点轨迹也较为相似.北斗GEO和IGSO卫星的星下点轨迹见图 1，其中IGSO星下点轨迹呈8字型，且5颗IGSO卫星轨迹重合.

图 1

Fig. 1

图 1 北斗GEO和IGSO星下点轨迹 Fig. 1 Subsatellite point track of Compass GEO and IGSO

LEO卫星的轨道参数设置见表 1.在对某一参数进行研究时，则把其它参数设置为默认值.由于是近圆轨道，轨道偏心率被设置为一很小的常数.仿真时首先利用卫星开普勒轨道参数，计算不同时间北斗卫星、LEO卫星的位置和速度.然后计算两颗卫星和地球表面的几何关系，计算地心到两卫星连线切点距离和切点高度，同时计算出北斗卫星相对于LEO卫星的仰角和方位角.当满足如下条件时，认为发生一次掩星事件：

表 1(Table 1)

表 1 LEO卫星的轨道参数设置 Table 1 LEO satellite orbit parameters settings

表 1 LEO卫星的轨道参数设置 Table 1 LEO satellite orbit parameters settings

（1）切点在北斗卫星和LEO卫星连线之间；

（2）切点高度在0~100km之间；

（3）北斗卫星相对于LEO的方位角在天线接收范围之内（前后向掩星天线极限功率波束宽都为40°）.

分别在不同的LEO卫星轨道的近地点角距和升交点赤经条件下，计算一天中北斗掩星事件数量和分布的变化.从仿真结果可以看出，近地点角距和升交点赤经对北斗卫星的掩星事件切点的经度及纬度分布的变化和掩星事件数量和分布影响并不明显.

不同近地点角距条件下对掩星事件数量的仿真结果见图 2a，由图中可以看出，一天内的北斗掩星总数随LEO卫星的近地点角距变化不大，平均掩星次数约为600次，变化范围为±10次.其中上升掩星和下降掩星的数量都在300次左右.而不同的北斗卫星造成的掩星事件中，MEO卫星掩星事件数量最多，并稳定在445次左右，而GEO和IGSO掩星事件数量都稳定在80次左右.说明在不同的LEO近地点角距条件下，北斗的三种轨道卫星掩星事件数量都保持稳定.图 2b显示了一天内北斗掩星事件数量随LEO升交点赤经变化的规律.从图中可以看出，MEO掩星事件数量变化具有周期性，在LEO卫星的升交点赤经在60°、180°和300°时，无论是上升还是下降的MEO掩星事件数量都位于波谷，而MEO掩星数量的波峰位于波谷的左右15°.联系仿真中采用的MEO三个轨道的升交点赤经设置，可以发现在LEO轨道与MEO轨道的升交点赤经互补时，掩星事件数量达到最低，实际进行北斗掩星探测时应避免这种情况.北斗GEO掩星事件数量基本稳定不变.而北斗IGSO掩星事件数量变化不大.

图 2

Fig. 2

图 2 不同近地点角距和升交点赤经条件下的北斗掩星事件数量 （a）不同近地点角距；（b）不同升交点赤经. Fig. 2 Compass occultation number variation with LEO′ sargument of perigee and LEO′s right ascension of ascending node

按照表 1中的参数设置，将LEO轨道高度设置为在300~1500km间，以50km为阶梯变化，对北斗掩星事件进行了模拟计算.图 3显示了一天内不同LEO轨道高度条件下的北斗掩星数量.从图中可以看出，不论是GEO、IGSO还是MEO卫星，其掩星事件数量都随LEO轨道高度呈下降趋势.而总的北斗掩星事件数量从300km轨道的690次下降到1500km轨道的490次左右.

图 3

Fig. 3

图 3 不同LEO轨道高度条件下北斗掩星事件数量 Fig. 3 Compass occultation number variation with LEO′s orbit height

另一方面，在不同的LEO轨道高度条件下，北斗掩星事件切点的地理分布也稍受影响.图 4分别显示了LEO轨道高度为300km和1500km时的掩星事件切点分布.对两个高度的掩星事件切点的经度和纬度分布进行比较，可以发现不同轨道高度下掩星事件切点的经度分布变化不大；而轨道高度越低，掩星事件切点的纬度分布越集中，轨道高度变高，掩星事件切点的纬度分布趋于分散.且LEO轨道高度对GEO卫星掩星事件的影响相对较大，其切点分布变化较明显.

图 4

Fig. 4

图 4 不同LEO轨道高度条件下北斗掩星事件切点分布 （a）300km轨道高度；（b）1500km轨道高度. Fig. 4 Compass occultation latitudinal and longitudinal distribution with various LEO orbit height

设置LEO卫星轨道倾角为0°~90°，步长5°，其它默认参数设置见表 1，计算一天内北斗掩星事件随卫星轨道倾角的变化.从仿真中发现，在不同的轨道倾角条件下，北斗掩星事件数量和分布的变化较大。

北斗掩星事件数量随轨道倾角变化的曲线见图 5.图中可以看出，北斗MEO掩星事件数量随轨道倾角提高而增加，直到85°倾角附近开始小幅下降.而10°~20°倾角之间北斗MEO掩星事件数量有一个急剧增加的过程.北斗GEO掩星事件数量在倾角为30°以下时保持在115次左右，此后掩星事件数量随倾角增加而减小.IGSO掩星事件数量基本稳定.总的掩星事件次数的峰值在20°~30°倾角附近.

图 5

Fig. 5

图 5 LEO轨道倾角对北斗掩星数量的影响 Fig. 5 Compass occultation number variation with LEO′s orbit inclinations

而不同的的LEO轨道倾角条件下，北斗掩星事件的切点分布也变化很大.图 6和图 7分别显示了LEO轨道倾角为30°和90°时的北斗掩星事件分布.倾角为30°时的北斗掩星切点基本集中在北纬30°与南纬30°之间，北斗GEO掩星切点都分布在在固定经度点上，而北斗IGSO掩星切点基本集中在0~50°E和130~180°W两片区域.当LEO倾角为90°时，北斗掩星事件切点在高纬度地区的分布明显增加，可以探测到地球极地区域，但赤道区域的掩星事件较为稀少.与倾角30°时相比，北斗GEO掩星切点分布有明显变化，基本集中在极地附近，且纬度分布不像30°时一样集中.北斗IGSO掩星切点在LEO倾角为90°时变得较为发散.

图 6

Fig. 6

图 6 LEO轨道倾角为30°时北斗掩星切点分布 （a）GEO掩星事件及切点分布；（b）IGSO掩星事件及切点分布；（c）MEO掩星事件及切点分布；（d）掩星事件及切点分布. Fig. 6 Global distribution of Compass occultation (LEO orbit inclination is 30 degree)

图 7

Fig. 7

图 7 LEO轨道倾角为90°时北斗掩星切点分布 （a）GEO掩星事件及切点分布；（b）IGSO掩星事件及切点分布；（c）MEO掩星事件及切点分布；（d）掩星事件及切点分布. Fig. 7 Global distribution of Compass occultation (LEO orbit inclination is 90 degree)

在不同的倾角下，对北斗掩星事件切点的经度和纬度分布进行统计的结果见图 8，图中显示了LEO轨道倾角I为0°、30°、60°、90°时的北斗掩星切点分布.从图中可以看出：北斗GEO掩星切点分布的峰值位置与倾角选择关系密切，一般切点分布集中在北纬I与南纬I（I为LEO轨道倾角）附近；在低倾角情况下掩星事件切点分布集中在几个固定经度点上，随倾角增大，经度分布变得分散.北斗IGSO掩星切点较均匀地分布在北纬I到南纬I之间；在低倾角情况下，北斗IGSO掩星切点的经度分布集中在150°W和45°E附近，在高倾角时，其经度分布也变得较为均匀.北斗MEO掩星切点的经度分布很均匀；纬度分布则较均匀地分布在北纬I到南纬I之间，但在倾角较大的情况下，赤道区域的掩星切点密度较低.总的北斗掩星事件切点的经度和纬度分布与MEO类似，但在LEO轨道倾角较低时，由于GEO掩星事件分布较为集中，会在固定的经度和纬度点造成尖峰.

图 8

Fig. 8

图 8 不同LEO轨道倾角条件下北斗掩星的空间分布 （a）轨道倾角0°；（b）轨道倾角30°；（c）轨道倾角60°；（d）轨道倾角90°. Fig. 8 Compass occultation latitudinal and longitudinal distribution with various LEO orbit inclinations

本文模拟并研究了北斗卫星在不同的LEO轨道参数条件下的掩星事件数量和分布，从仿真结果中可以看出北斗掩星事件具有如下特点：

（1）GEO掩星事件特点：在不同的LEO近地点角距和升交点赤经条件下，北斗GEO掩星事件数量变化不大；北斗GEO掩星事件数量随LEO轨道高度和LEO轨道倾角升高而下降.北斗GEO掩星事件切点的分布与LEO近地点角距和升交点赤经的关系不大；LEO轨道高度越低，GEO掩星事件切点的纬度分布越集中，轨道高度变高，GEO掩星事件切点的纬度分布趋于发散；GEO掩星事件切点分布与LEO轨道倾角关系密切，纬度分布的峰值在对应的南纬I和北纬I附近（I为LEO轨道倾角），经度分布在低倾角时很集中，在高倾角时有所发散.

（2）IGSO掩星事件特点：在不同的LEO近地点角距和升交点赤经条件下，北斗IGSO掩星事件数量变化不大；北斗IGSO掩星数量随LEO轨道高度升高而降低；随LEO倾角升高而稍有增加.北斗IGSO掩星事件切点的分布与LEO近地点角距、升交点赤经和轨道高度的关系不大；IGSO掩星事件切点较均匀地分布在南纬I到北纬I的纬度范围内，经度分布集中150°W和45°E附近，经纬度分布都随LEO轨道倾角增加而变得分散.

（3）MEO掩星事件特点：在不同的LEO近地点角距条件下，北斗MEO掩星事件数量变化不大；在LEO轨道与MEO轨道的升交点赤经度数互补时，MEO掩星事件数量达到最低；MEO掩星数量随LEO轨道提高而大幅降低；随轨道倾角提高而增加，直到85°倾角附近开始小幅下降.北斗MEO掩星事件切点的经度分布在全球范围内较为均匀；纬度分布集中在南纬I到北纬I的范围内，较为均匀，但在LEO倾角较高时，低纬度区域掩星事件数量较少.

结合三种轨道北斗卫星的特点，可以看出北斗掩星事件的数量和分布主要取决于数量最多的MEO卫星.但可以利用GEO和IGSO掩星事件的特点增加特定区域的掩星事件密度.例如图 7中，当LEO轨道倾角为90°时，北斗IGSO掩星事件切点在中国及周边地区上空较多，有利于提高此区域掩星探测密度，且此倾角条件下北斗GEO掩星事件切点集中于高纬度地区，有助于探测极地气候.在实践中应该根据具体需求选择合适的LEO轨道以尽量发挥IGSO和GEO卫星的区域探测作用.