一种基于斜返电离图和层析成像的联合反演方法
【技术领域】
[0001]本发明涉及一种基于斜返电离图和层析成像的联合反演方法，属于智能控制领域。
【背景技术】
[0002]电离层是电波传播的介质，电离层探测系统工作时电离层是无线电信号的传输媒质，电离层中电子浓度的分布状态和起伏变化会会造成电离层折射指数的改变，从而对回波信号产生影响，导致回波信号相位扰动，回波信号多普勒谱展宽，信噪比降低。因此从广义上讲，掌握电离层电子浓度的分布状态与变化规律，有利于提高无线电通信、授时、定位、导航、遥测和预警等系统的精度。
[0003]在实际电离层层析问题中，由于地面GPS接收机布站不均且数量有限、电子密度重建所需的数据缺失比较严重、以及观测仰角有限等问题，从而导致基于GPS数据的电离层层析反演通常是一个秩亏的问题，也就是在电离层层析过程中有一些格网点没有任何GPS射线穿过，得不到任何观测信息来对初始值进行修正。因此能够将不同探测手段获得的电离层数据资料都利用起来进行联合反演和优化反演，对于电离层结构的准确认识和提前预测具备极强的科学意义和应用价值。
[0004]在传统的基于GPS数据的电离层层析成像实验中，由于待重建区域中的地面接收机数量比较少并且通常分布不均匀，再加上GPS射线高度角的有限性使得垂直射线和接近垂直的射线占绝大多数，只有与待反演区域上下边界都落在区域内的射线才能用来层析反演，因此很难获得待反演区域电离层内部结构的全面信息，所以将待反演电离层空间进行离散化之后，有些像素并没有GPS射线通过，即没有观测信息进行修正，那么该像素的电子密度值与初始值保持一致，没有任何改进，这就意味着单纯基于GPS数据的电离层层析成像的电子密度反演结果与电离层电子密度的真实分布状态普遍存在一定的误差。
[0005]另一方面，在电离层层析成像中迭代初始值很多情况下是由国际参考电离层模型给出的，但是因为IRI模型是一种统计经验模型，它只能反映出平静电离层状况下的平均状态，特别是当空间天气异常或存在扰动时，IRI模型提供的电子密度值与电离层真实值存在着比较大的差异。正因为这些问题的存在，重建得到的电离层电子密度分布图一般会有一定程度的失真现象，图像上直观的表现为相邻像素点间的电子密度值不连续，往往存在一些突变。
[0006]有鉴于此，有必要提供一种基于斜返电离图和层析成像的联合反演方法，以克服上述问题。
【发明内容】
[0007]本发明的目的是:为了克服现有技术存在的问题，本发明提出一种基于斜返电离图和层析成像的联合反演方法。
[0008]本发明所采用的技术方案是:一种基于斜返电离图和层析成像的联合反演方法，其特征在于，对层析模型各个像素点电子密度初始值，先利用斜返反演进行修正，增加先验信息，然后再进行层析反演，层析反演基于像素基模型，利用GPS观测数据求解各条射线的电离层TEC值;然后根据卫星信号射线求出系数矩阵;最后反求电子密度。
[0009]如上所述的基于斜返电离图和层析成像的联合反演方法，其特征在于，所述联合反演方法的具体步骤如下:
(I)确定电离层反演区域，并根据电子密度分布规律对待反演区域进行相应的网格划分；
(2 )利用IRI2012模型获得各格网点的电子密度初始值，并把初始值代入第一步建立的像素基电离层模型中；
(3)对斜返电离图进行反演以获得电离层电子密度的实时分布，把所得的反演结果代入电离层模型从而对各格网点的电子密度初始值进行初步修正；
(4)利用GPS观测数据求出在基于像素指标函数的电离层层析反演中待反演区域范围内每条GPS射线对应的绝对TEC值；
(5)根据WGS-84大地坐标系计算出GPS卫星的瞬时位置，然后根据已知的地面站坐标求解出每条GPS射线在每个格网内截距的长度，从而得出系数矩阵；
(6)利用MART算法进行迭代求解，重建反演区域内的电子密度分布。
[0010]如上所述的基于斜返电离图和层析成像的联合反演方法，其特征在于，步骤(I)中，划分的格网间隔在经度和玮度方向上都取为1°，根据电离层电子密度的分布特性，高度方向上的间隔在500km以下取为10km，500?100km取为50km，100km以上取为200km，因此，待反演区域被划分成共12826个网格点。
[0011]如上所述的基于斜返电离图和层析成像的联合反演方法，其特征在于，步骤(4)中，每条GPS射线上的TEC值等于该条射线穿过的相应像素内的截距与该像素内的电子密度值的乘积之和，将范围扩大，m条射线的情况下，TEC的值可用方程组表示，进一步写为矩阵形式，为通过GPS观测求得的m条信号线穿过电离层区域的TEC值所构成的列向量；由于采用像素指标函数，那么每个格网内的电子密度值在某一时刻可假设是常数并均匀分布的，确定了电离层反演区域和反演时刻之后，需要进行层析成像中很重要的工作一一层析方程的建立:首先统计出完整穿过待反演区域的GPS射线数目，即该GPS射线与电离层反演区域最上和最下两个球面的交点均落在反演区域内，然后按照一定的规则对所有的格网点及GPS射线进行排序，为建立层析方程及下一步的迭代求解做准备。
[0012]如上所述的基于斜返电离图和层析成像的联合反演方法，其特征在于，步骤(5)中，为了满足电离层层析成像的精度要求，在计算GPS卫星瞬时位置时考虑地球自转的影响，忽略极移和岁差的影响。
[0013]如上所述的基于斜返电离图和层析成像的联合反演方法，其特征在于，步骤(6)中，每一步修正对应一次TEC测量，即针对一个方程进行，每m步迭代为一轮迭代，m为总的GPS射线数目，在迭代过程中修正的依据是第k步迭代产生的电离层电子密度对应的斜距TEC与实际观测求出的斜距TEC之比，从而对电离层电子密度逐步做出相应的修正，使图像趋近于最终解。
[0014]本发明的有益效果是:本发明的方法是一种将斜返扫频电离图与电离层层析成像相结合的优化算法，通过VFSA算法对扫频电离图的前沿进行反演得到实时电离层参数来构建电离层模型，然后利用GPS数据进行电离层层析成像，从而将GPS数据与斜返数据进行同化，使之更加接近电离层实际状况，这也是本发明的创新之处。与传统电离层层析方法相比该算法精度有了一定的提高，与垂测观测数据符合的更好。在观测数据不足情况下，该联合反演算法能够克服传统的电离层层析方法过分依赖经验模型的缺点。并且，本发明是首次提出将斜向返回探测数据与GPS数据结合在一起的联合反演算法，有效地改善了 GPS数据不充足的问题。
【附图说明】
[0015]图1是本发明的基于斜返电离图和层析成像的联合反演方法的方法流程图。
【具体实施方式】
[0016]为了更好地理解本发明，下面结合实施例进一步阐明本发明的内容，但本发明的内容不仅仅局限于下面的实施例。本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样在本申请所列权利要求书限定范围之内。
[0017]在本实施例中，如图1所示，本发明提出了一种将斜返扫频电离图与电离层层析成像相结合的联合反演方法，总的思路是:在IRI2012模型给出层析模型各个像素点电子密度初始值之后，先利用斜返反演结果对其进行修正，增加先验信息，然后再进行层析反演，从而克服那些没有任何GPS射线穿过的像素对经验电离层模型给出的电子密度迭代初始值的依赖。实验结果表明，该方法能够提高电子密度的重建精度，是两种反演方法的改进和融入口 ο
[0018]鉴于像素基电离层层析具有简单、方便的特点，本发明所进行的电离层层析反演基于像素基模型