gps信号对流层大气折射率的影响及分段模型.docxVIP

gps信号对流层大气折射率的影响及分段模型 当gps卫星信号通过地球外层传输到地面上时，它会受到流量层的影响。对流层对GPS信号的影响从天顶方向的2 m到地平高度角15°的25 m,即使通过模型修正,也只能修正其中的90%~95%;特别是湿大气部分,由于变化随机性强,很难用模型来估计。所以对流层延迟的改正精度一直制约着GPS测量精度的提高。然而由于对流层误差影响的复杂性,现有对流层模型的精度难以满足高精度基线解算的要求。 本文分析了不同高程段对流层大气折射率的变化情况,建立了新的大气折射模型,在此模型基础上进行分段积分,得到了站点的对流层改正模型,最后利用实际数据对新模型和传统模型进行了比较,验证了新模型的有效性。 1 流层大气折射建模方法分析 1.1 大气折射模型 电磁波在大气中传播时,由于受对流层的影响而使传播速度及传播路径发生变化,以小于光速的群速在传播,而传播路径也将产生弯曲(一般小于1 cm,可以不予考虑),所以研究对流层对电磁波传播的折射影响,首先就要研究对流层大气折射率。大气折射率随高度h呈线性变化: 式中,N(h)为高度h处的折射率;N0为地面折射率;ΔN为折射率梯度;h为从地面起算的高度。 在近地面1 km范围内,此线性模型与实际数据之间吻合较好;若超出该范围,则模型精度将随高度的增加逐渐下降。为了得到高、低空都适用的大气折射指数模型,平均大气折射率可用负指数模型来近似: 式中,c为指数衰减系数(1/m);h0为地面的海拔高度(m)。此模型对年平均或月平均大气折射率误差很小,一般小于5 N单位。 上述两种模型各有适用的区段,但都不能完全表达大气折射率的整个剖面。虽然在不同区段内大气折射率随高程呈线性或负指数变化,但干湿项折射率随高程的变化速度是不同的,所以要将它们分别考虑。 1.2 计算模型及计算参数 基于以上的理论分析,并结合多地区、多时段探空数据,作者提出了大气折射率分段模型。该模型充分考虑了大气折射率随高程变化的特点,将大气从测站地面到测站处的对流层顶分为3层:第1层是从测站地面到测站上空1 km处,该区段取线性模型;第2层是地面上空1 km处到高程12 km处,该区段取指数模型;第3层是从高程12 km处到测站处的对流层顶,该段取另一种指数模型。这样就可以更精确地描述整个剖面。 干分量折射率公式: 由于湿分量变化很大,且随机性强,目前尚没有可靠的模型来进行模拟,本文使用与干分量类似的公式来对湿分量折射率进行解算: 式中,下标d和w分别代表干项和湿项;N(h)为海拔高度h处的大气折射率(N单位);N(h0)为地面处的大气折射率(N单位);ΔN1为近地1 000m内的折射率梯度(1 N/m);α1为地面上空1 000m至海拔12 000 m的指数衰减率(1/m);α12为海拔12 000~60 000 m的指数衰减率。和分别为各区段均质大气高度的平均值,干项均质大气高可用下式计算: 式中,Rd为干气体常数;gm为平均重力;φ和h分别为用户纬度和高程(m);T为绝对气温(K)。 对于湿项均质大气高,按照φ划分为3个带: 因为大气折射指数的年变化很小,模型中各参数可用各地区的年平均和不同季节平均值来求得。 1.3 新模型的大气分辨率模型 本文利用成都大气探空站(高程544 m,北纬32.7°,东经103.83°)2004年7~9月每天2个时次(北京时间08时和20时)的大气探空数据来对新模型的精度进行验证。不同高度折射率计算公式为: 式中,P和T分别为规定层气压和温度;ew为水气压;A=6.11,a=7.5,b=237.3;td为规定层露点温度。 设H=h0,则用实测数据和新模型计算所得成都每天的大气折射率差值如图1所示。其差值中误差见表1。 从图1和表1中可以看出,新模型对干分量折射率的符合精度明显高于湿分量。对7~9月的解算结果进行比较可以看出,干分量符合精度较稳定,随季节、时段的变化较小;而湿分量在7、8月(湿度较大)时精度较差,9月(相对较干燥)时精度又有所提高。而且新模型的湿分量折射率在08时的精度比20时的精度要高,主要是受大气湿度的影响,这也符合了成都地区降水主要集中在夜间的事实。由新模型计算得出的干分量和湿分量折射率中误差均值分别为2.7 N和4.5 N,说明了新的大气折射率模型具有很高的精度。 湿分量对流层延迟为: 式中,hd为对流层外边缘的高度(m),hd=40 136+148.72(T-273.16);T为高程H处的绝对气温;hw=12 000 m为对流层湿项的协议高度。 3 实验与结果分析 3.1 实验数据和软件 本文实验所选用的GPS站点为中国境内部分IGS跟踪站(BJFS、KUNM、LHAS和WU-HN)和四川连续观测网中的成都站(CHDU)。所采用的数据为各站2004年7~9月