基于三测站数据特征点的电离层扰动传播测量方法及系统技术方案

本发明专利技术公开了一种基于三测站数据特征点的电离层扰动传播测量方法及系统，包括建立局部坐标系，选取监测数据特征点，得到异常到达三个点的时间，得到异常扰动的传播时间差；求异常传播路径的斜率和扰动传播方向投影点坐标，求异常传播速度和方位角，得到测量结果。本发明专利技术利用测站穿刺点精确的几何关系，在电离层扰动时间较短的情况下，仅用少量特征点即可准确地自动解算出电离层波动传播的速度大小和方向。对于研究日地空间环境、电离层、磁层耦合、电离层行扰传播形态等具有重要科学意义和应用价值。

Ionospheric disturbance propagation measurement method and system based on data characteristic points of three stations

The invention discloses a three measuring ionospheric station data feature point measuring method and system based on disturbance propagation, including the establishment of local coordinate system, selecting feature points of monitoring data, get abnormal reached the three point time, spread time disturbances; coordinate slope and exception propagation path and disturbance propagation direction of projection. O, anomalous propagation velocity and azimuth angle measurement. The invention uses the precise geometric relation of the point of the measuring station, and when the ionospheric disturbance time is short, the velocity and direction of the ionospheric wave propagation can be accurately solved automatically with only a few characteristic points. It is of great scientific significance and practical value for the study of the solar terrestrial environment, ionosphere, magnetosphere coupling, ionospheric disturbances and propagation patterns.

【技术实现步骤摘要】
基于三测站数据特征点的电离层扰动传播测量方法及系统
本专利技术涉及电离层波动传播分析
，特别涉及一种基于三测站空间关系以及监测数据特征点的电离层扰动传播测量技术方案。
技术介绍
地球高层大气由于受到太严辐射作用，原子、分子发生电离形成大量自由电子，从而在60～2000km高度范围内形成电离层。电离层扰动是指电离层结构偏离其常规形态的急剧变化，又称电离层骚扰。电离源的突变、非平衡态动力学过程、不稳定的磁流动力过程和某些人为因素等，都可引起电离层扰动。电离层扰动会对卫星通讯、GNSS定位等产生严重影响。而寻找电离层扰动传播源，研究其传播方向、速度、范围等，对于探索磁层、电离层、对流层之间相互作用，空间天气变化规律具有重要意义。目前国际上采用的电离层扰动传播速度提取方法主要有距离-时间图法、三测站干涉测量法等。其中距离-时间图法主要是在图上按照斜率计算，准确性较差，计算结果偏差较大；三测站干涉测量法需要较长的稳定扰动波形，且计算过程十分复杂，易受外界因素影响，在扰动波形较短时计算结果不准确。因此，本领域亟待更实用的技术方案出现。
技术实现思路
为了解决现有方法中存在的问题，本专利技术提供一种实现简单、结果精确的电离层扰动传播测量技术方案。本专利技术技术方案提供一种基于三测站数据特征点的电离层扰动传播测量方法，包括以下步骤，步骤1，建立局部坐标系，包括根据已知的三个测站穿刺点的坐标，建立以其中一点A为原点的局部坐标系，经线方向为y轴，纬线方向为x轴，计算其余两个测站穿刺点B、C在该局部坐标系中的坐标(xB,yB)、(xC,yC)；步骤2，选取监测数据特征点，得到异常到达A、B、C三个点的时间分别为tA、tB、tC，得到异常扰动从A到B和从A到C的传播时间，分别记为时间差dtAB＝tB-tA，dtAC＝tC-tA；步骤3，求异常传播路径的斜率k和B、C在扰动传播方向投影点D、E坐标(xD,yD)、(xE,yE)如下，步骤4：求异常传播速度和方位角，得到测量结果如下，异常传播速度大小为异常传播方位角为其中，αAD、αAE为D、E坐标(xD,yD)、(xE,yE)相应的异常传播方位角。而且，步骤2中，选取监测数据特征点的实现方式为，从三个监测站的二阶差分时间序列取极小值或极大值，若三个曲线在该区域前后波形非常相似，选取该时域内极小值或极大值作为异常到达三个测站的时间。本专利技术相应提供一种基于三测站数据特征点的电离层扰动传播测量系统，包括以下模块，初始化模块，用于建立局部坐标系，包括根据已知的三个测站穿刺点的坐标，建立以其中一点A为原点的局部坐标系，经线方向为y轴，纬线方向为x轴，计算其余两个测站穿刺点B、C在该局部坐标系中的坐标(xB,yB)、(xC,yC)；时差提取模块，用于选取监测数据特征点，得到异常到达A、B、C三个点的时间分别为tA、tB、tC，得到异常扰动从A到B和从A到C的传播时间，分别记为时间差dtAB＝tB-tA，dtAC＝tC-tA；投影模块，用于求异常传播路径的斜率k和B、C在扰动传播方向投影点D、E坐标(xD,yD)、(xE,yE)如下，测量输出模块，用于求异常传播速度和方位角，得到测量结果如下，异常传播速度大小为异常传播方位角为其中，αAD、αAE为D、E坐标(xD,yD)、(xE,yE)相应的异常传播方位角。而且，时差提取模块中，选取监测数据特征点的实现方式为，从三个监测站的二阶差分时间序列取极小值或极大值，若三个曲线在该区域前后波形非常相似，选取该时域内极小值或极大值作为异常到达三个测站的时间。本专利技术利用三个监测站的扰动出现的时间差以及三个监测站的空间几何关系提取电离层扰动在空间中传播速度的大小和方向。该技术方案利用测站穿刺点精确的几何关系，在电离层扰动时间较短的情况下，仅用少量特征点即可准确解算出电离层波动传播的速度大小和方向，并且在实验中得到了验证。本专利技术对于研究日地空间环境、电离层、磁层耦合、电离层行扰传播形态等具有重要科学意义和应用价值，具有重要的市场前景。附图说明图1是本专利技术实施例的技术方案原理图；图2是本专利技术实施例应用示例的三个监测站以及震中分布图；图3是本专利技术实施例应用示例的三个GPS站在电离层扰动期间的监测数据二阶差分时间序列图；图4是本专利技术实施例应用示例的按照三站特征点法计算得到的电离层扰动传播的速度和方向结果图。具体实施方式下面通过实施例，对本专利技术的技术方案作进一步具体的说明。本专利技术首先根据已知的三个测站穿刺点的坐标建立以其中一点为原点的局部坐标系，然后根据三个测站监测数据特征点出现的时间计算电离层波动到达三个测站之间的时间差；最后根据三个测站穿刺点几何关系以及波动到达三个测站穿刺点的时间差在局部坐标系中解算出电离层波动传播的速度大小和方向。参见图1，本专利技术实施例提供的技术方案包括如下步骤：步骤1：建立局部坐标系从GPS卫星到测站的卫星信号与电离层的交点称为穿刺点，三个测站在同一时刻会形成三个穿刺点。为了简化后续解算过程，本专利技术选取三个穿刺点中的任意一点作为原点，经线方向为y轴，纬线方向为x轴，建立局部坐标系。同时解算出其他两个点在局部坐标系中的坐标。图1中A、B、C三点为监测点，以A点为原点建立局部坐标系，B、C坐标为(xB,yB)、(xC,yC)。步骤2：计算扰动的传播时间假设电离层异常传播方向如图1中所示，异常平行向前推进。从B点向异常传播方向做垂线，交点为D点；从C点向异常传播方向做垂线，交点为E点。则异常到达B点的时间和到达D点的时间相同；同样，异常达到C点的时间和到达E点的时间相同。具体实施时，可选取监测数据特征点确定时间，优选的实现方式为，从三个监测站的二阶差分时间序列取极小值或极大值，若三个曲线在该区域前后波形非常相似，选取该时域内极小值或极大值作为异常到达三个测站的时间。A、B、C三个测站异常特征点出现时间已经确定，设异常到达A、B、C三个点的时间分别为tA、tB、tC，则异常从A到B和从A到C的传播时间分别为步骤3：求异常传播路径的斜率k和除原点以外其余两点在扰动传播方向投影点坐标设异常到达D、E三个点的时间分别为tD、tE，异常从A到D和从A到E的传播时间分别为dtAD、dtAE，传播距离为SAD、SAE。因为tB＝tD，tC＝tE，则dtAD＝dtAB，dtAE＝dtAC。设异常速度大小为V，方位角为α。则由于时间差dtAD和dtAE都有方向，即正负，能够表示异常传播的方向。而A点为原点，D、E在X轴的坐标xD、xE和在Y轴的坐标yD、yE也有方向，可得到上式中，dtAD＝dtAB，dtAE＝dtAC均为已知，设异常传播方向直线的斜率为k，则过C点和过B点的垂线斜率为-1/k。其余两点B、C电离层异常传播方向投影点分别为D点和E点，D点为直线AD和直线BD的交点，由AD和AB方程可得D点坐标为同理可得E点坐标为公式(3)、(5)、(6)联立最终求得，以A点为原点求出异常传播的斜率为并且根据(9)求得异常传播斜率后，即可将k回代到式(5)、(6)得到D、E坐标(xD,yD)、(xE,yE)。步骤4：求异常传播速度和方位角由于求得k后，D点和E点坐标都可以根据式(5)、(6)求出来，则异常传播速度大小为异常传播方位角(从Y轴正方向起，依顺时针方向本文档来自技高网...

【技术保护点】
一种基于三测站数据特征点的电离层扰动传播测量方法，其特征在于：包括以下步骤，步骤1，建立局部坐标系，包括根据已知的三个测站穿刺点的坐标，建立以其中一点A为原点的局部坐标系，经线方向为y轴，纬线方向为x轴，计算其余两个测站穿刺点B、C在该局部坐标系中的坐标(x

【技术特征摘要】
1.一种基于三测站数据特征点的电离层扰动传播测量方法，其特征在于：包括以下步骤，步骤1，建立局部坐标系，包括根据已知的三个测站穿刺点的坐标，建立以其中一点A为原点的局部坐标系，经线方向为y轴，纬线方向为x轴，计算其余两个测站穿刺点B、C在该局部坐标系中的坐标(xB,yB)、(xC,yC)；步骤2，选取监测数据特征点，得到异常到达A、B、C三个点的时间分别为tA、tB、tC，得到异常扰动从A到B和从A到C的传播时间，分别记为时间差dtAB＝tB-tA，dtAC＝tC-tA；步骤3，求异常传播路径的斜率k和B、C在扰动传播方向投影点D、E坐标(xD,yD)、(xE,yE)如下，步骤4：求异常传播速度和方位角，得到测量结果如下，异常传播速度大小为异常传播方位角为其中，αAD、αAE为D、E坐标(xD,yD)、(xE,yE)相应的异常传播方位角。2.根据权利要求1所述基于三测站数据特征点的电离层扰动传播测量方法，其特征在于：步骤2中，选取监测数据特征点的实现方式为，从三个监测站的二阶差分时间序列取极小值或极大值，若三个曲线在该区域前后波形非常相似，选取该时域内极小值或极大值作为异常到达三个测站的时间。3.一种基于三测站数...

【专利技术属性】
技术研发人员：翟长治，马符讯，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北,42