星载等离子体探测仪回波的模拟研究

肖 雄，姜春华，杨国斌，赵正予,2

（1. 武汉大学 电子信息学院，武汉 430072；2. 哈尔滨工业大学（深圳），深圳 518055）

引 言

等离子体空间环境是空间环境的重要组成部分，在太阳辐射的作用下，该区域存在大量能够自由运动的电子和离子，其中，自由电子能够影响无线电波的传播特性。空间等离子体环境一般包括电离层、等离子体层和磁层。而等离子体层和磁层在电离层的上边界，并向外延伸至磁层顶，是连接深空等离子体环境与地球空间环境的过渡地区，主要由完全电离的电子、离子组成，中性大气成分非常少，带电离子的运动完全受地磁场控制。由于等离子体空间环境中的自由电子会影响在其中传播的无线电波的幅度和相位，在无线电深空探测和通信中，如果等离子体环境中存在电子密度的不规则体结构，那么这些不规则体结构将会严重影响深空探测和通信的性能。因此，空间等离子体环境的研究将为深空探测和通信工程的顺利实施提供重要的参考依据[1-2]。

随着现代电子技术的发展，空间等离子体环境的探测技术也逐渐发展起来，人类对于电离层的基本特征也有了丰富的认识，特别是随着卫星技术的发展，能够从电离层顶部对电离层进行观测研究。20世纪六七十年代，美国国家航空航天局（National Aeronautics and Space Administration，NASA）发射卫星“云雀1号”（Alouette-1）、“云雀2号”（Alouette-2）和ISIS-1、ISIS-2，搭载星载电离层顶部探测仪，用于探测顶部电离层。对于等离子体层和磁层观测，主要依赖于卫星的就位观测或无源等离子体探测手段，如WAVES系统[3]、PWI （Plasma Waves Instrument）系统[4]、范阿伦探针（Van Allen Probe）[5]等。搭载在IMAGE（Imager for Magnetopause-to-Aurora Global Exploration）[6]卫星上的RPI（Radio Plasma Imager）[7-8]是第一套主动发射电磁波探测等离子体层和磁层空间环境的星载探测系统。RPI是低功耗有源多普勒雷达系统，探测频率范围为3 kHz～3 MHz，RPI向自由空间中各个方向发射雷达电磁波，然后通过一副x-y轴为300 m、z轴为20 m的三正交天线接收三路回波信号[9]，可以从回波幅度、相位、延迟时间、频率、极化、多普勒频移和回波方向的数据中获得自由空间中反射点处等离子体的位置和特性。利用三根天线正交组成的深空天线组阵能更好地完成对深空等离子体的探测，更好地降低天线组阵的增益损失[10]。

随着卫星观测数据的增加，对等离子层的探索与模拟也愈发深入。Green等[11]通过对Angerami和Thomas[12]的扩散平衡模型、Kimura[13]的电离层和等离子层模型、Aikyo和Ondoh[14]的等离子层顶模型和Roelof和Sibeck[15]的磁层顶模型进行组合，构造背景等离子体环境，进一步完成射线追踪模拟。Huang等[16]通过对IMAGE卫星的实测回波数据进行分析得到了一种基于实测数据的等离子体层经验模型。

由于卫星就位观测模式对于研究等离子体环境的空间特征分布存在一定的缺陷，为进一步研究等离子体层和磁层中等离子体的物理过程及其和底部电离层之间的耦合关系，需要提供更为丰富的等离子体结构分布。星载等离子体探测仪可以获取电波传播路径上的电子密度分布特征，进一步丰富等离子体层和磁层的观测数据。为了从等离子体回波图中反演得到真实环境中的电子密度分布，首先需要对等离子体回波图有一定基本了解，并对比电子密度不同结构和等离子体回波图之间的关系。

星载等离子体探测仪的观测数据目前较少，主要结合等离子体层模型和射线追踪方法来研究探测仪回波特征结构。由于等离子体层的模型对于探测仪回波模拟有决定性作用，如果模型精度较差，则回波特征结构也和实际探测存在较大误差。因此本文采用较高精度的等离子体层模型来研究星载等离子体探测仪回波的结构特征。首先利用全球尺度等离子体模型（Global Core Plasma Model，GCPM）构造空间等离子体环境，然后利用基于冷等离子体的射线追踪方法对无线电波在等离子体环境中的传播进行模拟，同时加入无线电波在自由空间中传播的损耗，以此来仿真星载RPI的回波特征。同时，将RPI模拟回波和IMAGE卫星的实测结果进行了对比，初步结果表明，在中纬度地区和正北极区进行模拟传播时，群路径和传播损耗随着工作频率的增加而增加，整体趋势呈现抛物线型，与IMAGE卫星在极区地区探测的实测结果类似，这为后续从等离子体回波图中反演电子密度的分布提供了理论依据。

1 GCPM模型

近年来，很多全球尺度的等离子体模型相继被开发，如基于卫星观测数据的经验模型IMAGE/RPI[16]、全球等离子体模型（Global Plasmasphere Ionosphere Density，GPID）[17]、标准等离子体模型（Standard Plasmasphere Ionosphere Model, SIM）[18]，以及GCPM[19]。在这些模型中，GCPM模型提供的电子密度及其导数是连续的，而这是等离子体中波传播数值模拟的一个必需条件[20]。同时，GCPM模型是一个开源的模型，基于以上条件考虑，本文最终采用GCPM模型（v2.4）来模拟等离子体空间环境。GCPM模型由多个区域模型集合而成，包括电离层模型、等离子体层、磁层槽区和极盖区等区域。该模型能够根据地磁活动指数Kp和太阳辐射（F10.7）以及时间等信息输出不同形态的等离子体环境图像。GCPM模型作为一个全球尺度的等离子体模型，能够根据不同的地磁环境合理地估计整个内磁层的热等离子体密度。该模型没有明显的外部边界，在白天可以延伸到外磁层区域，而晚上主要在地磁场的磁力线范围内适用。图1为2002年7月4日标准时间0点、3点、6点、9点4个不同时间段，Kp为1时的背景等离子体密度分布情况，可以看出图示背景等离子体密度存在明显不均匀分层，不同时刻的分布大致相同。图2为2002年7月4日标准时间0点时刻赤道区和极区位置电子密度随高度变化的剖面结构。

图1 2002年7月4日不同时间段背景等离子体密度分布Fig. 1 Distribution of background plasma density at different time periods on July 4th, 2002

图2 2002年7月4日0点时极区和赤道地区电子密度剖面图Fig. 2 Electron density profile at polar region and equatorial region at 0 UT on July 4th, 2002

2 基于冷等离子体的射线追踪方法

3 仿真结果分析

3.1 仿真条件

本文以2002年7月4日标准时间0点时刻的GCPM模型为背景等离子体环境进行电波传播的模拟。在模拟过程中，考虑两个不同地区即极区上空和中纬度地区上空的等离子体环境进行模拟。在图3中，发射仰角0°对应于x轴正方向，按逆时针方向增大，颜色深浅代表着射线发射仰角的大小，颜色由浅至深代表着发射仰角从小到大的变化。在极区上空的模拟中，设置卫星发射点位置为（3RE，8RE），射线群路径的步进为0.01RE，当返回的射线距离发射点位置小于0.1RE时，则判断该射线能够被卫星接收到。设置工作频率分别为0.1 MHz和1 MHz，射线的发射仰角以0.05°的步进从225°向245°变化，极区的射线传播路径图分别对应于图3（a）（0.1 MHz）和图3（b）（1 MHz）。从图3（a）和图3（b）中可知，由于极区部分存在等离子体密度的蝴蝶扇区分布，使得部分射线能返回到卫星的位置，而不是全部在极区上空反射远离发射点。图3（c）和图3（d）所示是在中纬度地区上空的位置，卫星发射点位置为（7.3RE，3.7RE），射线群路径的步进为0.01RE，同理，设置工作频率分别为0.1 MHz和1 MHz，射线的发射仰角以0.05°的间隔从210°向230°变化，中纬度地区射线传播路径图分别对应于图3（c）（0.1 MHz）和图3（d）（1 MHz）。

图3 处于不同发射起点的工作频率分别为0.1 MHz和1 MHz的射线0点时刻在GCPM中传播的寻常波的路径图Fig. 3 Path diagram of ordinary waves propagating in GCPM at 0 o'clock time of rays with working frequencies of 0.1 MHz and 1 MHz at different emission starting points respectively

为了获取RPI的等离子体回波图，本文将工作频率的起始频率设置为50 kHz，终止频率设置为1MHz，频率间隔为1 kHz，其它参数保持不变，群路径则由式（15）计算得到[24]， τ0表示起始射线和卫星的距离， τr表示射线反射点距离卫星的距离，在本文中τ0为0

图4所示为在极区上空等离子体回波图的模拟结果，从图4中可知，在极区和在中纬度地区的回波图存在较大的差异，不再是一段平滑的抛物线。可以明显看出在部分频率范围内，群路径随工作频率的增长反而会减小。这是因为在极区地区，电子密度的分布会更加复杂，在极区上空存在着电子密度的凹槽。从图3（a）、图3（b）所示的射线传播路径图可以看出，群路径的变化即可以看作是射线反射点的变化，工作频率从0.1 MHz到0.3 MHz增长时，射线的反射点逐渐靠拢，对应着群路径也随之增大，但是由于极区地区电子密度分布复杂，存在电子密度凹槽区，所以反射点的位置也会随电子密度的变化而变化，表现出来的便是在0.3 MHz到0.65 MHz的频率区间内有随工作频率的增加，群路径反而减小的趋势，群路径在0.65 MHz以后继续随工作频率的增加而增大。图4所示为在中纬度地区上空模拟的结果，随着工作频率的增加，群路径和能量损耗也会随之增加。电波反射点越靠近，电子密度随空间的变化越小，所以群路径在工作频率达到0.8 MHz时趋于平稳状态。由于在中纬度地区的上空，电子密度随空间变化比较平稳，因此等离子体回波图是一段较为平滑的抛物线，如图5所示。通过中纬度地区和极区上空等离子体回波图的模拟，发现等离子体回波图能够很好地反映等离子环境中电子密度的分布特征，可以通过反演等离子体回波图来获取真实环境中电子密度随空间的分布结构特征。

图4 0点时刻极区上空的等离子体回波图的模拟结果Fig. 4 Simulation results of plasma echo patterns at polar regions at 0 UT

图5 0点时刻中纬度上空的等离子体回波图的模拟结果Fig. 5 Simulation results of plasma echo patterns at mid-latitude region regions at 0 UT

3.2 与在轨实测数据的对比研究

图6所示为IMAGE卫星RPI在轨实测回波图[25]，卫星处于北极上空，频率范围为50～300 kHz，横坐标代表频率变化，右侧纵坐标代表回波的群路径，颜色强度则是代表了相对幅度的大小。图7所示是处于相同位置的等离子体回波图的模拟结果，卫星发射点位置为(0,3.5RE)，频率为0.05～1 MHz。从图6可知，在自由空间中传播时，随着工作频率的增加，相对幅度越来越小，意味着回波损耗越来越大，该在轨实测等离子体回波图在50～300 kHz的频率范围内也是一段较为平滑的抛物线结构。在轨实测图中，频率为0.05MHz时对应的群路径约为0.8RE；频率为0.3 MHz时，群路径约为2.3RE。而在模拟结果图7中可以看出，频率为0.05 MHz时，计算出的群路径约为0.77RE；频率为0.3 MHz时，群路径约为2.0RE。在轨实测图的群路径变化趋势、变化范围以及回波损耗等特征都与图7的模拟结果类似。在轨实测图存在等离子体共振回波，而本文的模拟结果没有考虑共振回波。

图6 IMAGE卫星RPI实测回波图[25]Fig. 6 RPI measured echo map by IMAGE satellite[25]

图7 距离北极上空3.5RE处等离子体回波图模拟结果Fig. 7 Simulation results of plasma echo pattern at 3.5RE above the North Pole

4 结 论

本文利用等离子体模型GCPM模型（v2.4）作为背景模拟等离子体空间环境，然后通过射线追踪方法对无线电波在等离子体环境中的传播进行模拟，同时加入无线电波在自由空间传播的损耗，来仿真星载等离子体探测仪RPI的回波图。本文初步给出了不同地区（极区和中纬度地区）上空等离子体回波图的结构特征。从模拟仿真结果看，中纬度地区和正北极上空等离子体回波图的结构特征都和实测结果的回波类似，而从极区右上空接收的等离子体回波图则更为复杂，这是由于极区的电子密度分布更为复杂。模拟结果说明等离子体回波图能够很好地反映等离子体环境中电子密度在空间中的分布结构。根据本文的工作基础，后续将研究如何从等离子体回波图中反演得到电子密度随空间的分布，进一步为未来深空探测任务和通信工程的顺利实施提供参考。