本发明涉及测量技术领域,特别涉及一种gnss(globalnavigationsatellitesystem,全球导航卫星系统)单站位移监测方法。
背景技术:
近年来,随着gnss接收机性能和解算技术的不断发展进步,gnss精确定位方法取得了系列进步,其应用领域也得到了极大拓展,gnss测量技术大致分为静态和动态测量两类。
静态测量需要观测时间长且同步观测接收机较多,数据处理需要专门软件进行解算、平差后才能应用于工作,主要用于高精度控制测量领域。
动态测量观测时间较短,现场可以快速得到测量结果,但因其精度较低主要在导航、测绘领域中放样、测图、变形监测等方面得到广泛应用,动态测量又大致分为以下几种方式:gnss单点定位(也叫绝对定位)、双差相对定位(rtk)、非差精密单点定位(ppp)。
gnss单点定位应用只需要一台gnss接收机即可测量快速测量定位,但只能达到米级定位精度。
双差相对定位(rtk)需要最少一台基站和至少一台流动站接收机,将基站架设在已知点上进行测量,测量时间几秒到几十秒,解算即可得到厘米级精度的结果,目前主要用于测绘领域精度要求较高的地方和工程测量施工放样。
非差精密单点定位(ppp),比静态测量观测时间短,比双差相对定位(rtk)观测时间长,定位精度可以达到厘米级精度,可用于建筑物进行了结构健康监测。
但是,以上三种定位方式都存在局限性:gnss单点定位因定位精度在米级,限制其一般只能用于导航领域,不能满足测图、放样等工作需要;双差相对定位(rtk)定位精度满足测图、放样等工作需要,但是需要基站作为参考站,所得到的动态定位结果是相对于参考站的位移情况,基站是必需的,要进行精确定位需要至少一台基站接收机和至少一台流动站接收机,要进行一系列基站设置、校正等准备工作后才能进行测量,因此增加了应用成本,作业效率低;非差精密单点定位(ppp)定位不需要基站,定位方式灵活,可用观测数据多,但是ppp定位模式需要接收精密星历,并且受到很多误差因素的影响,数据处理过程变得相对复杂,模糊度的快速固定也没有得到很好解决,与网络rtk相比,ppp定位精度较低,初始化时间更长,以至于至今都没有很好的商业接收机产品问世。
技术实现要素:
本发明的目的在于克服现有技术中所存在的现有gnss单点定位精度在米级,不能满足测图、放样等工作需要;rtk精度在厘米级,但是需要至少一台基站接收机和至少一台流动站接收机,要进行一系列基站设置、校正等准备工作后才能进行测量,因此增加了应用成本,作业效率低;ppp精度在厘米级,但是需要接收精密星历,并且受到很多误差因素的影响,数据处理过程变得相对复杂,模糊度的快速固定也没有得到很好解决,初始化时间较长的上述不足,提供一种gnss单站位移监测方法。
为了实现上述发明目的,本发明提供了以下技术方案:
一种gnss单站位移监测方法,将gnss接收机由a点移动至b点,a点观测时刻为t1,b点观测时刻为t2,移动过程中保持所述接收机对卫星的跟踪,还
包括在a点处和b点处分别建立该点处所述接收机与卫星相位观测之差的单差观测方程的步骤;
包括将a点处的单差观测方程和b点处的单差观测方程二次求差,得到双差观测方程的步骤;
包括所述接收机利用不同频率信号消除电离层因素对所述双差观测方程影响的步骤;
包括计算改正所述双差观测方程中卫星钟差变化的步骤;
包括将所述双差观测方程线性化,得到误差方程式的步骤;
包括通过最小二乘平差原理,获取精确的基线向量的步骤。
进一步的,所述gnss包括美国gps、俄罗斯glonass、欧盟galileo、中国beidou等卫星导航系统。
进一步的,所述卫星包括观测卫星s和参考卫星r,且参考卫星r的高度角大于观测卫星s的高度角。
进一步的,所述接收机在a点与观测卫星s的相位观测通过方程
其中,
那么,同理得出所述接收机在a点与参考卫星r的相位观测方程为
进一步的,所述接收机在a点与观测卫星s相对于参考卫星r的相位观测之差通过方程
其中,δ表示单差算子符号,
那么,同理得出所述接收机在b点与观测卫星s相对于参考卫星r的相位观测之差方程为
进一步的,b点处单差观测方程减去a点处的单差观测方程形成如下的所述双差观测方程
由于a点到b点的过程中,所述接收机都保持对卫星的跟踪,从a点到b点的距离和时间都较短,很多误差都是时空相关误差,可通过观测值历元间求差的方法得以消除,比如相对论效应、潮汐效应、轨道误差和对流层延迟等,且整周模糊度不会发生变化,均表示为
进一步的,由于电离层延迟变化较快,在数分钟内,即可发生较大变化,故需要使用双频进行消除,根据所述接收机的不同波长λ1和λ2,分别列出两个波长的所述双差观测方程:
进一步的,为了消除电离层
其中,if表示与电离层无关组合。
进一步的,通过二次多项式模型
其中,a1和a2分别是观测卫星s的钟速和钟漂,t0是卫星星历的参考时刻,t1是所述接收机位于a点的时刻,t2是所述接收机位于b点的时刻。
进一步的,将所述双差观测方程线性化,得到误差方程式包括以下步骤:
把已知项放在方程左边,把未知项放在方程右边,得到方程
将该方程线性化,得到误差方程式
其中,
进一步的,通过观测多颗卫星建立观测误差方程组,利用最小二乘平差原理,解算得到a、b两点基线向量的改正数。
进一步的,还包括利用卫星的不同高度角,对观测值进行加权的步骤,以进一步提高解算的基线向量的改正数的精度。
进一步的,根据a点和b点的初始坐标计算a、b两点的初始基线向量,再利用所述基线向量改正数,得到a点和b点精确的基线向量。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
运用本发明所述的一种gnss单站位移监测方法,适用于使用单个gnss接收机在短时间内测量空间距离较近的两点相对位置,依据对流层延迟等误差在短时间内较强相关性和没有周跳发生时模糊度固定不变的特点,提出gnss单站位移监测技术,研究对gnss单站数据时间序列信息进行处理,通过削弱或消除观测量中的时空相关误差,解算精确的单个接收机位移的基线向量,快速测量得到两点的相对位置关系;本发明技术不需要基站和通讯设备,不依赖接收精密星历,使用单个接收机即可达到厘米级定位精度,处理方法简单,节省设备购置成本和提高测量生产效率。
附图说明
图1为实施例中gnss接收机在t1时刻的a点移动至t2时刻的b点的示意图(本发明所述方法外业作业示意图);
图2为实施例中本发明所述方法外业作业流程图;
图3为实施例中本发明所述方法内业作业流程图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
实施例
如图1-3所示,本发明所述的一种基于美国gps的gnss单站位移监测方法,包括以下步骤:
步骤一、在a点设置gnss接收机,启动所述接收机进行数秒短时间的观测,a点与观测卫星s的相位观测可表示为:
其中,
同理可得a点与参考卫星r(参考卫星r的高度角大于观测卫星s的高度角)的相位观测方程为:
步骤二、对不同卫星间观测值求差,得到a点与观测卫星s相对于参考卫星r的观测方程为:
其中,δ表示单差算子符号,
步骤三、保持所述接收机对卫星的跟踪,将所述接收机迅速移动到b点(若移动过程中所述接收机与卫星断连,重新执行所述步骤一),进行数秒短时间的观测,则b点与观测卫星s相对于参考卫星r的观测方程为:
步骤四、由于a点到b点的过程中,所述接收机都保持对卫星的跟踪,从a点到b点的距离和时间都较短,很多误差都是时空相关误差,可通过观测值历元间求差的方法得以消除,比如相对论效应、潮汐效应、轨道误差和对流层延迟等,且整周模糊度不会发生变化,均表示为
其中,
步骤五、由于电离层延迟变化较快,在数分钟内,即可发生米级的较大变化,故需要使用双频进行消除,根据不同波长λ1和λ2,可分别列出双差观测方程:
步骤六、为了消除电离层
其中,if表示与电离层无关组合,由于gps信号基频f0=10.23mhz,f1=77*2f0,f2=60*2f0,因而
步骤七、卫星钟差的变化可通过二次多项式模型进行计算改正:
其中,a1和a2分别是观测卫星s的钟速和钟漂,t0是卫星星历的参考时刻,t1是所述接收机位于a点的时刻,t2是所述接收机位于b点的时刻;
步骤八、把已知项放在方程左边,把未知项放在方程右边,得到方程:
将该方程线性化,得到如下方程式:
其中,
步骤九、卫星到所述接收机几何距离的变化和卫星钟差的变化体现在相位观测值的变化上,通过观测多颗卫星建立观测误差方程组,利用卫星的不同高度角,对观测值进行加权,利用最小二乘平差原理,解算得到a、b两点基线向量的改正数,根据a点和b点的初始坐标计算a、b两点的初始基线向量,再利用所述基线向量改正数,由此得到a点和b点精确的基线向量:
δxab=xb0-xa0+dx,
δyab=ybo-yao+dy,
δzab=zb0-za0+dz。
运用本发明所述的一种gnss单站位移监测方法,通过编程实现其算法。在实际数据测试中,通过相邻两个点之间的相对位置测试和逐点之间(若干点两两之间)进行相对位置测试,通过本发明所叙述的方法得到的距离值和利用高精度全站仪实测结果进行比较分析,较差差值大部分能达到厘米级;本发明适用于使用单个gnss接收机在短时间内测量空间距离较近的两点相对位置,本发明技术不需要基站和通讯设备,不依赖接收精密星历,使用单个接收机即可达到厘米级定位精度,处理方法简单,节省设备购置成本和提高测量生产效率。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。