CN108957488B - 一种基于相位屏理论的电离层不规则体的漂移速度计算方法 - Google Patents
一种基于相位屏理论的电离层不规则体的漂移速度计算方法 Download PDFInfo
- Publication number
- CN108957488B CN108957488B CN201810901703.XA CN201810901703A CN108957488B CN 108957488 B CN108957488 B CN 108957488B CN 201810901703 A CN201810901703 A CN 201810901703A CN 108957488 B CN108957488 B CN 108957488B
- Authority
- CN
- China
- Prior art keywords
- ionospheric
- signal
- index
- change
- phase
- 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.)
- Active
Links
Images
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/14—Receivers specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及基于相位屏理论的电离层不规则体漂移速度计算方法,得相位变化与传播路径之间的关系;在接收机端可将相位的空间变化转换为相位变化与时间之间的关系,与电离层电子含量TEC的时间变化相关联;将电离层电子含量TEC的时间变化与ROTI指数时间变化率相联系;根据电离层幅值闪烁指数与信噪比,得出电离层不规则体在IGRF坐标系下的漂移速度以及其在x向,y方向上的分量与ROTI的时间变化率、电离层幅值闪烁指数、信噪比以及卫星仰角之间的函数式。本发明得到电离层不规则体的漂移速度与ROTI指数的时间变化率、电离层幅值闪烁指数S4、信噪比以及卫星仰角之间的关系,该方法可实现电离层不规则体及引起的电离层闪烁之间耦合机理的深入研究。
Description
技术领域
本发明涉及一种基于相位屏理论的电离层不规则体漂移速度计算方法,属于电离层领域。
背景技术
目前,卫星导航技术已成为世界各国在现代化战争中的重要武器组成部分。北斗卫星导航系统是中国自行研制、自主控制的卫星定位导航系统,已经实现亚太地区的定位、导航、授时、短报文等服务;BDS将在我国经济军事社会的发展中发挥重大作用。
卫星导航系统的空间部分大都位于20000公里以上,导航信号的传输都经过电离层,不可避免要受到电离层闪烁效应的影响。电离层闪烁影响无线电波的传播与电离层总电子含量TEC以及电离层不规则体相关。由于电离层电子浓度的不均匀,其中的相位折射指数通常小于1,会使穿过它的无线电波相比于自由空间存在一定的相位超前和相位色散;另外,电离层作为各向异性介质,其中传播的无线电波通常可分解为两个特征波模,两种波模具有不同的相速和传播路径,从而造成多普勒频移。所以,当无线电波穿过电离层时,其载噪比会减小,严重影响接收端信号的质量。同时电离层电子含量存在各种时空变化,即形成电离层不均匀体,其会受到太阳活动和地磁活动的影响,具有电离层暴、突然骚扰(SID)和电离层行扰(TID)等扰动现象,使电离层的形态更加复杂化。研究表明,电离层会对频率低于10GHz的无线电波传播产生明显影响,对于频率在L波段的全球卫星导航系统(GNSS)电离层的影响更大。
深入理解和认识电离层不规则体的形成和演化规律及其对卫星导航性能的影响,必须对其物理特征进行完备表征。电离层不规则体的磁纬向漂移速度定义在地磁坐标系下,且垂直于地磁子午向,是反映电离层不规则体空间尺度、演化速度及内部等离子体分布的重要物理特征,其大小及变化受到相关领域研究人员的广泛关注。目前研究人员(参见Priyadarshi,S.(2015).A review of ionospheric scintillation models.Surveys inGeophysics,36(2):295-324.
Carrano,C.S.,Groves,K.M.,Rino,C.M.,Doherty,P.H.(2016).A technique forinferring zonal irregularity drift from singlestation gnss measurements ofintensity(s4)and phase(σφ)scintillations.Radio Science,51(8):1263-1277.)主要通过地基GNSS观测网络、地面观测雷达及电离层测高仪来获得电离层不规则体的纬向漂移速度;美国波士顿学院的学者最新提出利用电离层闪烁监测仪获得的电离层闪烁指数和导航观测数据来反推电离层不规则体的形态,其方法具有一定的启示性,但其计算需要高频率(50Hz)的原始导航观测量,因此降低了方法的普适性。因此,设计并实现一种更简捷、普适和有效的电离层不规则体纬向漂移速度计算方法是本领域研究人员致力解决的难点之一。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种基于相位屏理论的电离层不规则体漂移速度计算方法,该方法在相位屏理论的基础上建立了电离层不规则体漂移速度与电离层幅值闪烁指数S4、信噪比以及电离层ROTI指数的时间变化率之间的映射关系,通过测量相关物理量来实现电离层不规则体纬向漂移速度的有效估计,并进一步解释了电离层幅度闪烁指数S4及电离层ROTI指数之间的相关性机理。本发明所需物理量可在具备电离层闪烁监测能力的GNSS闪烁监测接收机直接测量,为相关领域的研究提供方便可靠的纬向漂移速度计算方法。
本发明的技术方案如下:一种基于相位屏理论的电离层不规则体漂移速度计算方法,具体实现步骤为:
(1)根据一维相位屏理论模型,假设相位变化只发生在沿信号传播的x方向上,则信号沿传播路径上的相位变化可以表示为:
式中,I表示接收到信号的幅值强度,I0表示不存在闪烁现象时接收到的信号幅值,λ为信号的波长,φ表示信号的相位,s为信号的传播路径向量。
(2)由相位屏理论可知,在信号传播方向上的相位变化与该路径上的电离层总电子含量TEC的变化ΔNT(x)相关,即为:
φ(x)=-λreΔNT(r)
其中re=2.8×10-15为电子的半径。
(4)进一步公式可写为:
其中为电离层相位屏穿透点与卫星接收机之间的距离,Re为地球半径hipp为电离层穿透点高度,一般选取hipp=450km,v为在IGRF地磁参考坐标系下的电离层不规则体漂移速度,为卫星到观测站之间视线的天顶角,θ为观测站处穿透点的仰角,dROTI为ROTI指数的变化率。如图2所示。
(5)由此计算的电离层不规则体漂移速度为相对IGRF地磁参考坐标系的,如图2所示,在此坐标系下各地磁分量的关系如下式:
公式中Bx,By,Bz分别为x,y,z方向上的地磁强度。地磁偏角D和地磁倾角I定义如下:
因此电离层穿透点处观测的速度在IGRF地磁参考坐标系下可如下描述:
vipp=(vp,x+vi,x)x+(vp,y+vi,y)y
式中vp,x,vp,y分别为由于卫星在x(北向)和y(东向)方向运动而造成的电离层不规则体在x(北向)和y(东向)方向上的速度,vipp为电离层穿透点的速度,vi,x,vi,y分别为电离层不规则体在x(北向)和y(东向)方向上的漂移速度。
(6)进一步,从卫星到接收站的视线路径可以表示为:
k=-sin(δ)cos(αipp)x-sin(δ)sin(αipp)x+cos(δ)z
v=vipp·k。
本发明与现有技术相比的优点在于:
(1)相比传统的计算方法,本发明(图1所示)的优越性在于可以利用简单电离层闪烁监测仪的观测物理量来直接计算得到电离层不规则体的纬向漂移速度及时空分布,相比其它方法,具有简单、普适的优点,并具有较好的测量精度。如图1所示。
(2)相比传统的计算方法,本发明可有效研究不规则体纬向漂移速度的时间分布特征,利用长期观测获得电离层闪烁指数和GNSS导航观测量计算得到的电离层总电子含量,可实现电离层不规则体纬向漂移速度的时间特征和变化趋势,可获得指定区域范围内电离层不规则体的活动情况,其效果如图3所示。
(3)相比传统的计算方法,本发明建立了不规则体漂移速度与ROTI指数、S4指数、卫星导航信号接收信噪比之间的直观映射关系,有效解释了目前相关领域对ROTI与闪烁指数之间关系的一系列困惑,对电离层不规则体物理特征的完备表征具有较强的启发性。
附图说明
图1为本发明方法实现流程图;
图2为不规则体漂移速度计算空间示意图;
图3为不规则体漂移速度随月份及日小时的变化规律;
图4(a)分析了不规则体漂移速度与ROTI指数、仰角及S4指数之间的相关性;(b)分析了不规则体漂移速度与ROTI指数、信噪比α和S4之间的相关性。
具体实施方式
下面将结合附图及具体实施方式对本发明加以详细说明,需要指出的是,所描述的实施例仅旨在便于对本发明的理解,而不起任何限定作用。
本发明提出了一种基于相位屏理论的电离层不规则体漂移速度计算方法,该方法详细的描述了电离层不规则体漂移速度与卫星仰角、信噪比、幅值闪烁指数以及ROTI的变化率之间的关系。能够较好的描述电离层形态变化对闪烁产生的影响。
如图1所示,本发明具体实现步骤如下:
1、根据一维相位屏理论模型,假设相位变化只发生在沿信号传播的x方向上,可以得到信号的相位变化与传播路径的关系:
式中,I表示接收到信号的幅值强度,I0表示不存在闪烁现象时接收到的信号幅值,λ为信号的波长,φ表示信号的相位,r为信号的传播路径向量,s为一维相位屏理论中信号穿过电离层时的穿透点与卫星接收机之间的距离。
2、由相位屏理论可知,在信号传播方向上的相位变化与该路径上的电离层总电子含量TEC的变化ΔNT(x)相关,即相位变化为:
φ(x)=-λreΔNT(r)
其中re=2.8×10-15为电子的半径。
3、将相位在空间中的变化转换为随时间的变化并将ΔNT做近似处理ΔNT=ΔTEC,因此得到即相位变化与ROTI随时间的变化率相关。其中ΔNT为与信号传播路径相关的TEC变化;其中i表示不同的GNSS卫星,TECk+1和TECk分别表示在k+1时刻和k时刻的电离层总电子含量的值,Δtk=tk+1-tk为k+1时刻与k时刻之间的时间间隔;公式中的<ROT>为N个历元中ROT的平均值。
4、标准化接收信号的功率I0,即<I0>=I0,并认为<I0-I>=<I0>-<I>,其中I0为无闪烁时接收机接收到的信号功率,I为闪烁发生时接收机接收到的信号功率;
5、计算GNSS的信噪比与幅值闪烁指数最终得到其中为电离层相位屏穿透点与卫星接收机之间的距离,Re为地球半径hipp为电离层穿透点高度,一般选取hipp=450km,为卫星到观测站之间视线的天顶角,θ为观测站处穿透点的仰角,dROTI为ROTI指数的变化率,v为在IGRF地磁参考坐标系下的电离层不规则体漂移速度,dROTI为ROTI指数的变化率。
6、进一步根据电离层穿透点的速度计算不规则体漂移速度在IGRF地磁参考坐标系下各个方向上的分量。在此坐标系下各地磁分量的关系如下式公式中Bx,By,Bz分别为x,y,z方向上的地磁强度。地磁偏角D和地磁倾角I定义如下:因此电离层穿透点处观测的速度在IGRF地磁参考坐标系下可如下描述:vipp=(vp,x+vi,x)x+(vp,y+vi,y)y,式中分别为由于卫星在x(北向)和y(东向)方向运动而造成的电离层不规则体在x(北向)和y(东向)方向上的速度,vipp为电离层穿透点的速度,vi,x,vi,y分别为电离层不规则体在x(北向)和y(东向)方向上的漂移速度。
本发明提出了一种基于相位屏理论的电离层不规则体漂移速度计算方法,如图2所示,(a)图为电离层漂移空间示意图,R为地球半径,h为穿透点到地球表面的高度,其中P为电离层穿透点的位置(经纬高),R为接收机的位置(经纬高),I为地磁倾角,B地磁向量,k为卫星到接收站的视线路径,δ为卫星到观测站之间视线的天顶角,2β为接收机处穿透点的仰角,Vi表示不规则体表漂移速度,x、y、z分别为漂移速度在北向、东向、地向上的分量;(b)图为卫星到接收站的视线向量k在水平面上的分量;(c)图为卫星到接收站的视线向量k在竖直面上的分量,βp为卫星信号与水平面的夹角。
此方法详细的描述了电离层闪烁与信噪比、不规则漂移速度以及ROTI时间变化率之间的关系。
由仿真结果如图3和图4所示,图3计算了澳大利亚WEIP(左侧)与DWNI(右侧)观测站2011-2015年期间的S4(第1行)、ROTI(第2行)以及不规则体漂移速度(第3行)的时间分布图,图中可见,不规则体漂移速度与幅值闪烁指数S4以及ROTI指数存在相似的分布特征,随着S4与ROTI值的增大,漂移速度通常也会加速。图4中(a)图分析了不规则体漂移速度与ROTI、仰角以及S4之间的相关性,由图可知,在低仰角处更容易发现较大的漂移速度和较高的S4,(b)图分析了不规则体漂移速度与ROTI、信噪比α以及S4之间的相关性,可知在信噪比较低时,不规则体漂移速度较大更容易导致强烈的电离层幅值闪烁,这与实际情况相符合。
综述,本发明提出了一种利于电离层闪烁指数和电离层ROTI指数来直接测量电离层不规则体纬向漂移速度的方法,为电离层不规则体漂移速度的研究提供了较大的便利;该方法同时揭示了电离层闪烁指数、电离层ROTI指数、卫星信号信噪比和卫星仰角之间的直观映射关系,为进一步实现电离层不规则体物理特征的完备表征提供直接理论参考和技术依据。
以上所述仅为本发明的具体实施例,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于相位屏理论的电离层不规则体漂移速度计算方法,其特征在于,包括以下步骤:
A、根据一维相位屏理论模型,假设相位变化只发生在沿信号传播的x方向上,建立信号的相位变化与传播路径之间的关系式,具体表示为其中I表示闪烁发生时接收到信号的幅值强度,I0表示不存在闪烁现象时接收到的信号幅值强度,λ为信号的波长,φ表示信号的相位,r为信号的传播路径向量,s为一维相位屏理论中信号穿过电离层时的穿透点与卫星接收机之间的距离;
B、将上述信号相位变化与传播路径之间的关系式依据一维相位屏理论等效成电离层电子含量TEC的变化与传播路径之间的空间函数,在信号传播路径上的相位变化与该路径上的电离层电子含量TEC的变化ΔNT(x)相关,即相位变化为φ(x)=-λreΔNT(r),其中λ为信号的波长,re=2.8×10-15为电子的半径;
C、将ΔNT做近似处理ΔNT=ΔTEC,从而将TEC在空间中的变化转变为随时间的变化,根据一维相位屏理论将上述电离层电子含量变化与传播路径之间的空间函数重新建立为电离层电子含量变化的时间函数,具体表示为其中λ为信号的波长,s为一维相位屏理论中信号穿过电离层时的穿透点与卫星接收机之间的距离,re=2.8×10-15为电子的半径,v为在IGRF地磁参考坐标系下的电离层不规则体漂移速度;
D、根据ROTI指数的定义将C中电离层电子含量TEC变化的时间函数等效为ROTI指数的时间变化率指数,并将其定义为dROTI指数,即ΔNT为与信号传播路径相关的TEC变化;其中i表示不同的GNSS卫星,TECk+1和TECk分别表示在k+1时刻和k时刻的电离层电子含量的值,Δtk=tk+1-tk为k+1时刻与k时刻之间的时间间隔;公式中的<ROT>为N个历元中ROT的平均值,dROTI为ROTI指数的变化率;
E、将C中公式做近似处理,即认为<I0-I>=<I0>-<I>,其中I0为无闪烁时接收机接收到的信号功率,I为闪烁发生时接收机接收到的信号功率,并标准化接收信号功率I0,即<I0>=I0,得到公式其中λ为信号的波长,s为一维相位屏理论中信号穿过电离层时的穿透点与卫星接收机之间的距离,re=2.8×10-15为电子的半径,v为在IGRF地磁参考坐标系下的电离层不规则体漂移速度;
F、根据电离层幅值闪烁指数S4以及GNSS信噪比α的定义并替换上式中的I0与I,得到电离层不规则体在IGRF地磁参考坐标系下漂移速度与幅值闪烁指数S4和GNSS信噪比α以及ROTI指数时间变化率之间的关系式v为在IGRF地磁参考坐标系下的电离层不规则体漂移速度,dROTI为ROTI指数的变化率,λ为信号的波长;
G、根据电离层穿透点的速度计算公式可以将上式分解,最终分别得到在IGRF地磁参考坐标系下x方向,即北向和y方向,即东向上的电离层不规则体的漂移速度。
3.根据权利要求2所述的基于相位屏理论的电离层不规则体漂移速度计算方法,其特征在于:所述hipp=450km。
4.根据权利要求1所述的基于相位屏理论的电离层不规则体漂移速度计算方法,其特征在于:所述步骤G中,根据电离层穿透点的速度计算不规则体漂移速度在IGRF地磁参考坐标系下各个方向上的分量,在此坐标系下各地磁分量的关系如下式公式中Bx,By,Bz分别为x,y,z方向上的地磁强度,地磁偏角D和地磁倾角I定义如下:电离层穿透点处观测的速度在IGRF地磁参考坐标系下描述为:vipp=(vp,x+vi,x)x+(vp,y+vi,y)y,式中vp,x,vp,y分别为由于卫星在x和y方向运动而造成的电离层不规则体在x和y方向上的速度,vipp为电离层穿透点的速度,vi,x,vi,y分别为电离层不规则体x方向即北向和y方向即东向上的漂移速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810901703.XA CN108957488B (zh) | 2018-08-09 | 2018-08-09 | 一种基于相位屏理论的电离层不规则体的漂移速度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810901703.XA CN108957488B (zh) | 2018-08-09 | 2018-08-09 | 一种基于相位屏理论的电离层不规则体的漂移速度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108957488A CN108957488A (zh) | 2018-12-07 |
CN108957488B true CN108957488B (zh) | 2021-04-16 |
Family
ID=64467929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810901703.XA Active CN108957488B (zh) | 2018-08-09 | 2018-08-09 | 一种基于相位屏理论的电离层不规则体的漂移速度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108957488B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113253326B (zh) * | 2021-05-16 | 2021-12-21 | 中国矿业大学 | 基于测地型接收机的电离层不规则体漂移速度估计方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106134520B (zh) * | 2009-05-11 | 2014-05-28 | 中国电子科技集团公司第二十二研究所 | 中低纬地区电离层闪烁短期预报方法 |
CN103502844B (zh) * | 2011-03-25 | 2016-06-29 | 欧洲空间局 | 用于确定物体位置的方法、设备及系统 |
CN103675775B (zh) * | 2013-12-12 | 2016-01-20 | 北京理工大学 | 背景电离层对geo sar成像影响的分析方法 |
CN104933737B (zh) * | 2015-06-03 | 2016-04-27 | 北京航空航天大学 | 一种基于共轭梯度法的电离层层析成像混合反演方法 |
CN107561534B (zh) * | 2017-08-25 | 2020-08-04 | 北京航空航天大学 | 一种基于全极化高轨sar的电离层时变tec测量方法 |
-
2018
- 2018-08-09 CN CN201810901703.XA patent/CN108957488B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108957488A (zh) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9395443B2 (en) | Adaptive method for estimating the electron content of the ionosphere | |
Hoque et al. | Higher order ionospheric effects in precise GNSS positioning | |
Yang et al. | Global autonomous positioning in GNSS-challenged environments: A bioinspired strategy by polarization pattern | |
Nie et al. | An offshore real-time precise point positioning technique based on a single set of BeiDou short-message communication devices | |
CN109945870B (zh) | 一种伪距观测值和载噪比融合的伪卫星室内定位方法 | |
CN109975749A (zh) | 一种校正源存在条件下的短波单站直接定位方法 | |
CN106886039B (zh) | 基于城市三维地图的地面数字电视广播与导航卫星混合粒子滤波定位方法 | |
CN109613565A (zh) | 基于多星座gnss的电离层层析方法及系统 | |
CN114488228B (zh) | 一种适用于动态载体平台的gnss多路径误差削弱方法 | |
Chen et al. | Combination of ground-and space-based data to establish a global ionospheric grid model | |
JP2010112842A (ja) | Tecマップ及び受信機バイアスの作成及び測定方法及び装置 | |
CN108957488B (zh) | 一种基于相位屏理论的电离层不规则体的漂移速度计算方法 | |
Huang et al. | Analysis and improvement of ionospheric thin shell model used in SBAS for China region | |
CN111123345B (zh) | 一种基于gnss测量的经验电离层模型数据驱动方法 | |
CN110058270B (zh) | 一种基于钟差拟合的导航欺骗信号产生方法 | |
CN110146904B (zh) | 一种适用于区域电离层tec的精确建模方法 | |
Liu et al. | On the correlation between ROTI and S4 | |
Yoon et al. | Validation of ionospheric spatial decorrelation observed during equatorial plasma bubble events | |
CN113589392B (zh) | 基于数据仿真验证三频信标tec反演精度的方法及系统 | |
Mendoza et al. | A Method in Determining Ionospheric Total Electron Content Using GNSS Data for non-IGS Receiver Stations | |
Bo et al. | Research and experimental validation of the method for barometric altimeter aid GPS in challenged environment | |
Bishop et al. | GPS TEC and the plasmasphere: Some observations and uncertainties | |
CN118036338B (zh) | 一种基于掩星数据的电离层精细化三维重构方法 | |
Januszewski | How the ionosphere affects positioning solution using terrestrial and satellite navigation systems? | |
Pi | Comparisons of TEC Gradients Derived from Several GNSS-Based Approaches |
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 |