CN106646527A - 基于三测站数据特征点的电离层扰动传播测量方法及系统 - Google Patents
基于三测站数据特征点的电离层扰动传播测量方法及系统 Download PDFInfo
- Publication number
- CN106646527A CN106646527A CN201610874353.3A CN201610874353A CN106646527A CN 106646527 A CN106646527 A CN 106646527A CN 201610874353 A CN201610874353 A CN 201610874353A CN 106646527 A CN106646527 A CN 106646527A
- Authority
- CN
- China
- Prior art keywords
- propagation
- points
- time
- disturbance
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
Abstract
本发明公开了一种基于三测站数据特征点的电离层扰动传播测量方法及系统,包括建立局部坐标系,选取监测数据特征点,得到异常到达三个点的时间,得到异常扰动的传播时间差;求异常传播路径的斜率和扰动传播方向投影点坐标,求异常传播速度和方位角,得到测量结果。本发明利用测站穿刺点精确的几何关系,在电离层扰动时间较短的情况下,仅用少量特征点即可准确地自动解算出电离层波动传播的速度大小和方向。对于研究日地空间环境、电离层、磁层耦合、电离层行扰传播形态等具有重要科学意义和应用价值。
Description
技术领域
本发明涉及电离层波动传播分析技术领域,特别涉及一种基于三测站空间关系以及监测数据特征点的电离层扰动传播测量技术方案。
背景技术
地球高层大气由于受到太严辐射作用,原子、分子发生电离形成大量自由电子,从而在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轴正方向起,依顺时针方向到异常传播方向线之间的水平夹角)αAD、αAE可由D点坐标和E点坐标求出,最终以A点为原点的求出的异常传播方位角为
同理,可以将B、C设为原点,计算异常传播的速度和方位角,实现方式相同。
采用以上方案进行实验,采用中国大陆构造环境监测网络西藏地区三个GPS测站的数据。研究对象为2015年尼泊尔地震引起的同震电离层扰动。
步骤1:建立局部坐标系
选取图2中三个GPS测站xzar、xzzf、xzrk中任意一个测站的穿刺点为原点,记为点A,其他两个点记为A点和B点;经线方向为y轴,纬线方向为x轴,建立局部坐标系,同时计算其余两个GPS站的穿刺点在该局部坐标系中的坐标。
步骤2:选取监测数据特征点,计算扰动的传播时间。
从图3中三个监测站的二阶差分时间序列可以看出,三个箭头位置均出现了极小值,且三个曲线在该区域前后波形非常相似,因此选取该时域内最小值作为异常到达三个测站的时间,得到异常到达A、B、C三个点的时间分别为tA、tB、tC,则计算时间差dtAB=tB-tA,dtAC=tC-tA。
步骤3:求异常传播路径的斜率k和B、C在扰动传播方向投影点D、E坐标
步骤4:求异常传播速度和方位角
求得k后,D点和E点坐标都在步骤3已求出来,则异常传播速度大小为
异常传播方位角可由AD或AE求出,最终
最终计算结果如图4所示,V~=930m/s,α(Azi)~=56°,图中曲线为该时间段内三个GPS测站穿刺点的轨迹,箭头代表电离层扰动的方向。
具体实施时,本发明所提供方法可基于软件技术实现自动运行流程,也可采用模块化方式实现相应系统。
本发明实施例相应提供一种基于三测站数据特征点的电离层扰动传播测量系统,包括以下模块,
初始化模块,用于建立局部坐标系,包括根据已知的三个测站穿刺点的坐标,建立以其中一点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);
测量输出模块,用于求异常传播速度和方位角。
各模块具体实现可参见相应步骤,本发明不予赘述。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的。因此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (4)
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.一种基于三测站数据特征点的电离层扰动传播测量系统,其特征在于:包括以下模块,
初始化模块,用于建立局部坐标系,包括根据已知的三个测站穿刺点的坐标,建立以其中一点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)相应的异常传播方位角。
4.根据权利要求3所述基于三测站数据特征点的电离层扰动传播测量系统,其特征在于:时差提取模块中,选取监测数据特征点的实现方式为,从三个监测站的二阶差分时间序列取极小值或极大值,若三个曲线在该区域前后波形非常相似,选取该时域内极小值或极大值作为异常到达三个测站的时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610874353.3A CN106646527B (zh) | 2016-09-30 | 2016-09-30 | 基于三测站数据特征点的电离层扰动传播测量方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610874353.3A CN106646527B (zh) | 2016-09-30 | 2016-09-30 | 基于三测站数据特征点的电离层扰动传播测量方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646527A true CN106646527A (zh) | 2017-05-10 |
CN106646527B CN106646527B (zh) | 2019-01-04 |
Family
ID=58854675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610874353.3A Expired - Fee Related CN106646527B (zh) | 2016-09-30 | 2016-09-30 | 基于三测站数据特征点的电离层扰动传播测量方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646527B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111158021A (zh) * | 2018-11-08 | 2020-05-15 | 千寻位置网络有限公司 | 电离层行扰的估计方法及系统、预警终端 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002288774A (ja) * | 2001-03-27 | 2002-10-04 | Furuno Electric Co Ltd | 観測データ回収システム |
CN101893714A (zh) * | 2010-07-09 | 2010-11-24 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统广播电离层时延修正方法 |
CN102323598A (zh) * | 2011-07-29 | 2012-01-18 | 中国气象局北京城市气象研究所 | 一种电离层残差扰动量的检测方法、装置及系统 |
CN102928850A (zh) * | 2012-11-22 | 2013-02-13 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
JP2015138021A (ja) * | 2014-01-25 | 2015-07-30 | 国立研究開発法人 電子航法研究所 | 衛星航法システムにおける測位誤差の補正方法及びその装置。 |
-
2016
- 2016-09-30 CN CN201610874353.3A patent/CN106646527B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002288774A (ja) * | 2001-03-27 | 2002-10-04 | Furuno Electric Co Ltd | 観測データ回収システム |
CN101893714A (zh) * | 2010-07-09 | 2010-11-24 | 中国科学院测量与地球物理研究所 | 全球卫星导航系统广播电离层时延修正方法 |
CN102323598A (zh) * | 2011-07-29 | 2012-01-18 | 中国气象局北京城市气象研究所 | 一种电离层残差扰动量的检测方法、装置及系统 |
CN102928850A (zh) * | 2012-11-22 | 2013-02-13 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
JP2015138021A (ja) * | 2014-01-25 | 2015-07-30 | 国立研究開発法人 電子航法研究所 | 衛星航法システムにおける測位誤差の補正方法及びその装置。 |
Non-Patent Citations (1)
Title |
---|
姚宜斌等: "2015年尼泊尔地震的震前电离层异常探测", 《测绘学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111158021A (zh) * | 2018-11-08 | 2020-05-15 | 千寻位置网络有限公司 | 电离层行扰的估计方法及系统、预警终端 |
Also Published As
Publication number | Publication date |
---|---|
CN106646527B (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105136054B (zh) | 基于地面三维激光扫描的构筑物精细变形监测方法及系统 | |
CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
CN105759311B (zh) | 一种近实时地震震源位置定位方法 | |
CN105241465B (zh) | 一种道路更新的方法 | |
CN105158760B (zh) | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 | |
CN102591343B (zh) | 基于两行根数的卫星轨道维持控制方法 | |
CN102879820B (zh) | 基于三角网格的三维表层模型构建方法 | |
CN110007357B (zh) | 一种航空tem和航空mt联合反演方法 | |
CN104729486A (zh) | 一种基于似大地水准面精化的无验潮水深测量方法 | |
CN106291725A (zh) | 一种快速反演地下地质体空间位置的方法 | |
CN105093319A (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
CN113850908A (zh) | 一种考虑路径延长因素的地闪回击定位数据的优化方法 | |
Kanaujia et al. | 1D velocity structure and characteristics of contemporary local seismicity around the Tehri region, Garhwal Himalaya | |
CN104992054A (zh) | 基于时间序列二维化的电离层垂直总电子含量预报方法 | |
CN106556877A (zh) | 一种地磁通化方法及装置 | |
CN106874690B (zh) | 考量水平传播因素的地形重力波拖曳参数化方法 | |
Munro et al. | Divergence of radio rays in the ionosphere | |
CN101266153B (zh) | 测绘工程类陀螺全站仪精度评定方法 | |
CN105573963B (zh) | 一种电离层水平不均匀结构重构方法 | |
Solomentsev et al. | Three-dimensional assimilation model of the ionosphere for the European region | |
CN113960532A (zh) | 一种基于假想源的二次定位计算的微地震定位方法 | |
CN206459611U (zh) | 一种基于北斗卫星的超高层建筑高精度变形监测系统 | |
CN106646527A (zh) | 基于三测站数据特征点的电离层扰动传播测量方法及系统 | |
CN103745118B (zh) | 一种基于磁偶极子等效源法的地磁异常数据网格化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190104 Termination date: 20210930 |