基于数据仿真验证三频信标TEC反演精度的方法及系统与流程

基于数据仿真验证三频信标tec反演精度的方法及系统
技术领域
1.本发明涉及tec测量技术领域，特别是涉及一种基于数据仿真验证三频信标tec反演精度的方法及系统。


背景技术：

2.三频信标技术是一种新型的tec测量技术，该技术是通过安装在低轨卫星上的三频信标发射机，向位于地面监测站的接收设备发射一组相位相干的无线电信号(vhf、uhf和l频段)，来获得大范围和高精度的电离层电子总含量、电子密度剖面等参数信息的一种新型高效的空间环境测量技术。
3.2018年，中国地震局成功发射了我国第一颗电磁试验卫星zh
‑
1号。zh
‑
1号上搭载有高精度磁强计，gnss掩星接收机，三频信标发射机的等载荷，开启了我国自主探测空间地球物理场和各类电磁扰动信号的先河
1.。自1964年阿拉斯加大地震有学者首次发现地震可以引起显著的电离层扰动起
2.，地震可以引发电离层扰动已经成为共识。近些年出现了许多获取高分辨率电离层参数的方法：掩星反演方法
3.，三频信标时延法
4.，三频信标多普勒频移法
5.等。
4.tec可以反应电离层的主要特性，并且与在电离层中传播的无线电电波的时间延迟与相位延迟密切相关。因此电离层tec信息可以广泛应用于导航定位及通信的电波信号修正等领域，同时我们也可以根据电波传播的时间延迟与相位延迟来反演电离层tec。
5.三频信标技术是一种可以利用三种频率电波信号的传播频移反演得到高精度tec的新兴技术。austen等(1988)最早进行电离层层析成像技术(cit)的研究，此后cit技术成为卫星信标和地面接收站配套使用的探测手段
1.。双频信标技术是指通过计算两个频率的差分时延或者多普勒频移，从而得出链路上的电子总含量(tec)。采用三频信标技术时，同样的方法可以得出链路上的tec，并且可以得到三组值。经相关的数学处理能够提高tec的精度及可靠性。美国海军实验室的bernhardt
6.等人在灵活运用数论知识的基础上，进行了严格的数学推导，得出了tec的算法。该算法一方面较大幅度提高了相位积分常数确定的精度，另一方面也使得相位积分常数的求解难度大大降低，对于三频信标技术的发展有着重大意义。美国首先将这一创新成果应用于军事和气象学的研究当中。2006年成功发射的cosmic卫星上就搭载了三频信标发射机tbb，该项目与我国台湾地区合作建立了低纬电离层探测链，对低纬地区的电离层形态和扰动进行研究，可以有效监测地震高发区域的电离层结构变化
7.。为此，中国地震局也进行了自己的星基电离层探测计划，2018年发射升空的张衡一号电磁试验卫星上也搭载有由中国电科22所研制的三频信标发射机。
6.但是由于一定的技术原因，该发射机未能正常工作。因此，如何基于我国目前的现状提出一种基于数据仿真验证三频信标tec反演精度的方法，对于中国地震局电磁试验卫星随后的工作具有重要的指导意义。
7.相关文献如下：
8.[1]张学民,申旭辉,赵庶凡,刘静,欧阳新艳,娄文宇,泽仁志玛,何建辉,钱庚.地
震电离层探测技术及其应用研究进展[j].地震学报,2016,38(03):356
‑
375.
[0009]
[2]leonard r s,barnes jr r a.observation of ionospheric disturbances following the alaska earthquake[j].journal of geophysical research,1965,70(5):1250
‑
1253.
[0010]
[3]林剑,吴云,刘经南.电离层gps掩星反演技术研究[j].地球物理学报,2009,52(08):1947
‑
1953.
[0011]
[4]吴健.用三频卫星信标测量电离层天气新方法[j].中国科学a辑,2000,30(z1):111
‑
114.doi:10.3321/j.issn:1006
‑
9232.2000.z1.029.
[0012]
[5]赵海生,许正文,吴健,王占阁,刘琨,李娜,陈金松.三频信标高精度tec测量新方法[j].空间科学学报,2011,31(02):201
‑
207.
[0013]
[6]paul a.bernhardt,carl l.siefring.new satellite
‑
based systems for ionospheric tomography and scintillation region imaging[j].paul a.bernhardt；carl l.siefring,2006,41(5).
[0014]
[7]赵海生.三频信标电离层扰动探测及层析成像关键技术研究[d].杭州电子科技大学2010.


技术实现要素：

[0015]
本发明要解决的技术问题是提供一种基于数据仿真验证三频信标tec反演精度的方法，使其能为三频信标tec反演精度提供一种有效的验证方法，证明三频信标反演高精度tec的精度与可行性。
[0016]
为解决上述技术问题，本发明采用如下技术方案：
[0017]
一种基于数据仿真验证三频信标tec反演精度的方法，包括：
[0018]
(1)利用站星距与电离层背景信息仿真三种频率信号的多普勒频移；
[0019]
(2)利用高斯分布仿真相位观测噪音；
[0020]
(3)基于多普勒频移与相位观测噪音获取仿真的三种频率信号的差分多普勒相移；
[0021]
(4)基于仿真的三种频率信号的差分多普勒相移数据反演电离层tec。
[0022]
作为本发明进一步地改进，所述步骤(1)为：
[0023]
利用已知的历元实时卫星坐标，测站坐标，得到历元站星距；然后利用站星距与iri模型电离层背景信息仿真三种频率信号的多普勒频移；
[0024]
相邻历元两信号链路之间的多普勒频移：
[0025][0026]
其中，δf为多普勒频移，右式第一项是卫星运动引起的多普勒频移，第二项是电离层引起的多普勒频移；第一项与信号频率f成正比，第二项与之成反比，利用差分消除第一项的影响；对相邻链路的距离s做差，得到δs，以时间间隔1s作为时间增量δt，则：
[0027][0028]
同样的：
[0029][0030]
分别得到三种频率信号的多普勒频移。
[0031]
进一步地，所述步骤(2)为：
[0032]
分别以1
°
，3
°
，6
°
的高斯分布的相位误差加入仿真的信号相位，模拟信号的接收误差；以实现利用高斯分布仿真相位观测噪音。
[0033]
进一步地，所述步骤(3)为：
[0034]
由多普勒频移和信号初始波长λ可以计算信号的波长变化量δλ：
[0035][0036]
基于站星距、波长变化量以及相位误差，可获取仿真差分多普勒相移。
[0037]
进一步地，所述步骤(4)为：
[0038]
所述tec的解算公式为：
[0039]
tec＝8.3165
×
10
16
[(δφ
13
×7‑
δφ
12
×
8)
mod(1)
+k]
[0040]
其中，δφ
12
与δφ
13
分别表示差分多普勒相移δφ
12
与δφ
13
的小数部分；δφ
12
表示1频与2频的差分多普勒相移；δφ
13
表示1频与3频的差分多普勒相移；mod(1)表示对1求余；k表示三频相位积分常数。
[0041]
进一步地，所述三频相位积分常数k的计算方法为：通过双站或者多站法，基于电离层局部球对称假设，通过双频相位，估计粗tec值；通过三频相位tec表达式，对积分常数进行计算。
[0042]
本发明还提供了一种基于数据仿真验证三频信标tec反演精度的系统，包括：一个或多个处理器；存储装置，用于存储一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现上述的基于数据仿真验证三频信标tec反演精度的方法。
[0043]
通过采用上述技术方案，本发明至少具有以下优点：
[0044]
本发明主要基于仿真技术来实现对电波在电离层中传播引起的多普勒频移的模拟，来实现利用三频信标技术反演高精度电离层tec，与传统的双频信标技术相比，三频技术通过将三个频率两两差分大大提高了tec相位积分常数，从而有效地降低了相位积分常数的求解误差，提升了电子密度测量精度，本发明基于仿真验证证明了三频信标技术反演高精度tec的精度与可行性，对于中国地震局电磁试验卫星随后的工作具有重要意义。
附图说明
[0045]
上述仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，以下结合附图与具体实施方式对本发明作进一步的详细说明。
[0046]
图1是四测站双频反演tec误差平均值图；
[0047]
图2是四测站双频反演tec误差rms图；
[0048]
图3是四测站三频反演tec误差平均值图；
[0049]
图4是四测站三频反演tec误差rms图；
[0050]
图5是十三个三频信标接收站分布示意图；
[0051]
图6是十三测站三频反演tec误差平均值；
[0052]
图7是十三测站三频反演tec误差rms。
具体实施方式
[0053]
本发明基于三频信标技术，利用仿真的方法来模拟电波信号在电离层中传播引起的多普勒频移，利用多普勒频移来解算相对tec，然后基于相对tec使用多站法得到双频tec，将精确的多普勒相移小数部分与粗测tec结合求解高精度tec。证明了三频信标技术反演高精度tec的精度与可行性，对于中国地震局电磁试验卫星随后的工作具有重要意义。
[0054]
下面对本发明进行详细展开描述：
[0055]
1、三频信标测量tec原理
[0056]
1.1双频测量tec原理
[0057]
由于电离层的色散特性，电离层对于不同频率信号的传播相位影响是不同的。同时发出的两组不同频率的信号，它们到达接收站的相位是不同的而且相位差中包含有传播路径上的电离层信息。因此我们可以利用两组信号的相位差来反演传播路径上的tec。
[0058]
一般采用高频信号(vhf与uhf)，近似的高频折射系数公式：
[0059][0060]
其中，q表示电波传播方向与磁场之间的夹角，f
b
与q值可以通过磁场模型求得，进而可以由折射系数n来反演电离层电子密度。在这里我们假设传播的电波信号和介质无弹性碰撞，并且磁场对其传播过程影响很小或者无影响，则上式公式可简化为：
[0061][0062]
其中，n为折射系数，f为信号频率，f
p
为电离层等离子体频率。
[0063][0064]
其中，n
e
为电子密度，e为电子电荷，m为电子质量，e0为自由空间介电常数。可得：
[0065][0066]
得：
[0067][0068]
假设信号传播至接收站的相位为φ
r
，信号发出的相位为φ
s
，它们之间的关系是：
[0069][0070]
其中，ω为角频率，l表示卫星至接收站的光学路径距离。
[0071]
两边对时间微分：
[0072][0073]
可得，多普勒频移δf：
[0074][0075]
其中，前一项表示卫星运动引起的多普勒频移，后一项表示在电离层介质引起的多普勒频移。可以看出，卫星运动的影响与信号频率成正比，电离层介质影响与信号频率成反比。因此，我们可以通过将两种频率信号线性组合差分的方式消除卫星运动的影响。在实验中我们采用三种频率信号：f1＝150.012mhz，f2＝400.032mhz，f3＝1066.752mhz，这三种频率信号由f
r
＝16.668mhz的基频信号产生。
[0076]
对双频信号的多普勒频移进行差分：
[0077][0078]
两边积分得：
[0079][0080]
φ(t)表示差分多普勒相移，φ0表示相位积分常数，表示电子密度沿传播路径的积分，即为tec。
[0081]
1.2三频信标测量tec原理
[0082]
三频信标技术是在双频的基础上额外增加一个频率信号。原理如上，将三种信号两两组合差分，我们可以得到三组多普勒频移来进行tec测量：
[0083][0084][0085]
[0086]
可得：
[0087][0088]
其中，m1,m2与m3互质，可通过数学关系
6.将上两式联立得到三频信标解算tec的公式：
[0089]
tec＝8.3165
×
10
16
[(δφ
13
×7‑
δφ
12
×
8)
mod(1)
+k]
[0090]
其中，δφ
12
与δφ
13
分别表示δφ
11
与δφ
13
对2π取余，即差分多普勒相位的小数部分；mod(1)表示对1求余。k表示三频相位积分常数，类似于双频相位积分常数，在确定精确的k值后，可以得到高精度tec。
[0091]
2、仿真过程与结果
[0092]
2.1仿真原理
[0093]
本发明主要利用仿真的方法来研究三频信标技术反演高精度电离层tec。为了仿真信号传播过程中的多普勒频移，我们需要已知测站坐标，各测站接收历元实时卫星坐标，对应实时卫星高度角，以及iri模型背景电离层tec信息。本发明初始建立了由四个测站组成的观测网，网型呈长方形分布，四个测站坐标如下：
[0094]
测站坐标1234lat(n)24
°
26
°
24
°
26
°
long(e)100
°
100
°
105
°
105
°
[0095]
这四个测站的观测历元间隔为1s，测站与卫星之间的背景电离层tec由iri模型给出，站点坐标与实时卫星坐标已知。
[0096]
首先，利用已知的历元实时卫星坐标，测站坐标，得到历元实时站星距。然后利用站星距与电离层背景信息等来仿真三种频率信号的多普勒频移。
[0097]
相邻历元两信号链路之间的多普勒频移：
[0098][0099]
其中，δf为多普勒频移，右式第一项是卫星运动引起的多普勒频移，第二项是电离层引起的多普勒频移。第一项与信号频率f成正比，第二项与之成反比，因此可以利用差分消除第一项的影响。对相邻链路的距离s做差，得到近似的δs，以时间间隔1s作为时间增量δt，则：
[0100][0101]
同样的：
[0102][0103]
可以分别得到三种频率信号的多普勒频移，由多普勒频移可以计算信号的波长变化量δλ：
[0104][0105]
利用站星距、波长变化量以及相位误差，可以计算得到多普勒相移：
[0106][0107]
按照之前所述，将三种频率信号多普勒相移进行两两差分可以消除卫星运动的影响：
[0108][0109][0110]
分别得到1频信号(150mhz)与2频信号(400mhz)的差分相移，1频信号(150mhz)与3频信号(1066mhz)的差分相移。利用差分相移可以求得相对tec值：
[0111][0112]
由于积分常数未知，所以求得的tec仅仅是历元间相对tec，仅仅代表相对变化。如果要知道绝对tec则需要对积分常数进行估计。
[0113]
在得到相对tec之后，需要进行多站法估计积分常数。我们使用四个测站，大致呈正方形分布。首先我们需要根据已知信息计算电离层穿刺点坐标(高度设为350km)，求得四个测站每个观测历元的穿刺点坐标(ipp)(郭文兴.基于卫星信标的tec分布和电子密度反演技术研究[d].西安电子科技大学,2017.)：
[0114][0115][0116][0117]
其中，φ是地心角，表示穿刺点ipp至地心的向径与接收站至地心的向径之间的夹角。θ表示卫星仰角，α表示接收站处的卫星方位角。r表示地球半径，h表示穿刺点高度，我们取350km的平均电离层高度。
[0118]
上面提到，我们可以由多普勒频移得到相对tec，在确定相位积分常数之后，才能求得绝对tec。一般情况下，采用leitinger1975年提出的两站法(leitinger r,schmidt g,tauriainen a.an evaluation method combining the differential doppler measurements from two stations that enables the calculation of the electron content of the ionosphere[j].revista cient
í
fica,1975,41(2):201
‑
213)，两站法又可以扩展成为多站法。在这里我们选取多站法来进行相位积分常数的估计。
[0119]
由于我们的反演需要三种频率信号同步观测的结果，因此首先对测量结果进行时间历元筛选，筛选条件是卫星高度角高于15
°
，并且选取四个测站共同时刻历元数据，共计
285个历元。多站法的中心思想是利用最小二乘原理来进行拟合：
[0120]
对于相距较近的两个测站我们可以任务同一时刻测量的垂直tec值应该相等。因此我们选取每个历元距离较近的两个测站组成一组方程，共计570个方程组成方程组来进行拟合。
[0121]
可以利用穿刺点仰角三角函数将路径倾斜tec转化为垂直tec：
[0122]
tec
v
＝tec
s
*cosχ
[0123]
tec
s1
*cosχ1‑
tec
s2
*cosx2＝0
[0124]
建立多站法方程组，求解方程组得到估计的每个测站初始绝对tec值，再加上由上一步我们得到的历元间相对tec值，可以得到双频绝对tec值。
[0125]
与双频信号不同，利用三频信号的计算原理是先利用双频计算tec粗值，再利用三频精测tec的相对值。将差分相位看做整数项与小数项，类似于gps观测中的整周模糊度问题，我们可以利用迭代的方法来求出整数项。三频信标的相位积分常数远远大于双频相位积分常数，比较好确定，再利用小数项的组合可以精确确定相对tec值，两者相加得到高精度绝对tec值。
[0126][0127]
其中与分别表示1频与2,3频的差分多普勒相位的小数部分，k为三频相位积分常数，是一整数，具体k的求解方法为：通过双站或者多站法，基于电离层局部球对称假设，通过双频相位，估计粗tec值；通过三频相位tec表达式，对积分常数进行计算。
[0128]
2.2程序实现及结果
[0129]
按照以上的原理和数据处理步骤，利用fortran语言编写算法程序，对四个测站的仿真数据进行解算。输入数据包括测站的站点坐标，历元时刻与对应的卫星坐标、卫星高度角以及iri模型的背景电离层信息。输出利用三频技术解算的电离层tec值，并可与先验模型背景电离层信息做对比，得到三频算法的解算精度。
[0130]
程序功能主要包含三大部分：
[0131]
1.利用已知数据仿真三频信标信号传播的多普勒频移/时延，得到历元间差分相对tec
[0132]
2.利用多站法计算粗的双频tec
[0133]
3.利用仿真相对tec小数部分测定高精度tec
[0134]
首先，利用双频技术将150mhz信号与400mhz信号进行组合，仿真得到双频信号的多普勒相移，然后利用多站法估计相位积分常数，将这两者代入计算公式即可得到双频tec。由于双频tec计算模糊度(即相位积分常数)较大，使用多站法估计的相位积分常不够精确，因此反演精度较差，故称为双频粗tec。
[0135]
图1、2中显示使用双频技术反演的四个测站tec与背景电离层tec差值的统计值。本发明使用平均值与均方根值rms来评定反演算法的精度，其中包括tec差值绝对值的统计值和相对值的统计值。可以看出，使用双频方法时四个测站反演tec的差值绝对平均值在0.7个tecu至2.5个tecu之间，差值相对平均值在3％至8.5％之间，差值绝对rms在0.7到2.7个tecu，差值相对rms在3％到9％之间。其中，第三测站的反演精度较高，差值绝对值可以达到1个tecu以下，其他站则在1.5到2.5个tecu。并且，高斯误差的大小对双频反演精度没有
明显的影响。这是由于双频相位积分常数较小，相位积分常数占有决定性的影响，因此相对多普勒相移的精度对双频反演精度影响不大。
[0136]
在得到双频粗tec以后，利用流程图中所给出的方法，对高精度的三频tec进行求解。上边原理提到过，三频tec之所以精度高，是因为其相位积分常数较大，相应的相位积分常数比较好求，精确的多普勒相移小数部分可以决定反演tec的精度。
[0137]
如图3、4，利用三频技术反演得到的tec精度很高。在不添加误差的情况下，反演得到的tec与背景tec差值平均值可以达到2.0e
‑
07到2.5e
‑
07，差值相对值可以达到4.8e
‑
7到8.6e
‑
7％，rms可以达到3e
‑
7个tecu左右，相对rms可以达到6e
‑
07到8.5e
‑
7％，可以看出利用三频技术反演tec技术的精度十分高，算法除去计算误差，反演得到的tec与背景tec一致。在实际观测中，在现有接收仪器的观测条件下，接收到的信号相位会有3
°
到6
°
的观测误差，相当于1/120周到1/60周的相位捕捉误差。电波信号是一种正弦波，传播一周有2π的相位变化，即观测误差会导致π/60至π/30的相位误差。因此我们分别1
°
，3
°
，6
°
的高斯分布的相位误差加入仿真的信号相位，模拟信号的接收误差，得到加入高斯白噪声误差的三频技术反演结果。如图3，反演的tec值与背景tec之差均在1个tecu以下，效果良好。
[0138]
2.3十三测站仿真计算结果
[0139]
中国地震局在我国西部地区总共部署了13个地面测站来接收三频信标信号，大致分布于西北到西南，主要是我国地震高发地区，希望利用高精度电离层信息来研究地震前兆及地震引起的电离层变化。这13个测站分别为：大理站，固原站，合作站，玛曲站，盐池站，玉溪站，昭通站，中卫站，天水站，红河站，普洱站，虾拉坨站，绵阳站。具体分布示意图如图5所示。
[0140]
与4个实验测站的处理方法相同，首先利用仿真多普勒频移的方法得到相对tec，然后利用精确的多普勒相移小数部分与估计的相位积分常数结合求得高精度的三频tec。得到13个测站的三频tec反演结果统计如图6、7所示，从图中可以看出：观测误差为3
°
时，tec计算的相对误差平均值在2.5％以内，相对误差rms在3％以内；观测误差为6
°
时，tec计算的相对误差平均值在4.5％以内，相对误差rms在6％以内。三频信标技术反演电离层tec的精度还是比较高的。
[0141]
综上所述，本发明利用仿真的方法模拟了电波信号穿过电离层所产生的多普勒频移，利用程序实现了利用三频信标多普勒来反演高精度电离层tec的技术。通过对比双频反演结果，发现三频信标技术对于提升电离层tec反演精度效果显著。并且观测误差对于三频信标技术的影响要显著大于对双频技术的影响，这主要是由于这是由于双频相位积分常数较小，而三频相位积分常数较之双频有60多倍的提升，所以观测量的精度直接影响了三频反演的精度。因此，三频技术的的推广使用也需要相关观测技术和仪器的支撑。目前，相关仪器的观测精度可以达到近似于高斯3
°
到6
°
之间，根据本发明的计算结果，是可以满足需要的。
[0142]
基于上述方法，本发明还提供了一种基于数据仿真验证三频信标tec反演精度的系统，包括：一个或多个处理器；存储装置，用于存储一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现上述的基于数据仿真验证三频信标tec反演精度的方法。由于系统硬件部分属于本领域的常规技术，在此不再展开赘述。
[0143]
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，本领域技术人员利用上述揭示的技术内容做出些许简单修改、等同变化或修饰，均落在本发明的保护范围内。