GNSS网平差方法及装置制造方法及图纸

本申请公开了一种GNSS网平差方法及装置，该方法包括：根据已知点的速度参数和/或平面坐标，确定已知点的点位速度和方差-协方差阵；判断基线的观测历元是否相同；如果基线的观测历元不相同，则选取参考历元，计算在参考历元已知点的坐标；根据基线两端点的速度计算观测历元的基线两端点相对于参考历元中基线两端点的位移量；根据位移量修正对应基线的基线向量；根据基线的方差-协方差阵、基线两端点的方差-协方差阵、观测历元和参考历元确定修正后基线的方差-协方差阵；根据修正后的基线向量、修正后基线的方差-协方差阵、已知点的坐标和已知点的方差-协方差阵进行GNSS网平差。本申请可以消除由于点位位移引起的基线系统误差。请可以消除由于点位位移引起的基线系统误差。请可以消除由于点位位移引起的基线系统误差。

【技术实现步骤摘要】
GNSS网平差方法及装置


[0001]本申请涉及矿产勘探
，尤其涉及一种GNSS网平差方法及装置。

技术介绍

[0002]本部分旨在为权利要求书中陈述的本专利技术实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
[0003]随着大陆板块的运动，坐落在板块上的每个位置每年以不同的速度产生位移。以北京房山(BJFS)点为例，其点位移动速度为0.032米/年，拉萨(LHAZ)位移速度为0.049米/年，且与BJFS位移的方向并不一致。经过20年后BJFS坐标位移0.640米，LHAZ坐标位移0.980米。假设BJFS与LHAZ位移方向一致，则经过20年后会差生0.340米的差，如果位移方向相反，则会产生1.620米的差。
[0004]在工程应用领域，尤其是在大面积矿产勘探方面，必须布设全球导航定位系统(Global Navigation Satellite System，GNSS)网引入参考框架坐标，如ITRF框架，ITRF框架也是当前较为常用的参考框架。GNSS网中，利用基线来表征两点间的实时相对位置，本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种GNSS网平差方法，其特征在于，所述方法包括：获取基线参数、已知点的速度参数和/或平面坐标，所述基线参数包括基线的观测历元、基线向量、基线的方差-协方差阵和基线两端点的平面坐标；根据已知点的速度参数和/或平面坐标，确定已知点的点位速度和方差-协方差阵；判断基线的观测历元是否相同；如果基线的观测历元不相同，则从基线的观测历元中任意选取一个历元作为参考历元，根据已知点的点位速度、观测历元、参考历元和观测历元中已知点的坐标计算在参考历元已知点的坐标；根据基线两端点的速度计算观测历元的基线两端点相对于参考历元中基线两端点的位移量；根据位移量修正对应基线的基线向量；根据基线的方差-协方差阵、基线两端点的方差-协方差阵、观测历元和参考历元确定修正后基线的方差-协方差阵；根据修正后的基线向量、修正后基线的方差-协方差阵、参考历元已知点的坐标和已知点的方差-协方差阵进行GNSS网平差。2.根据权利要求1所述的方法，其特征在于，在判断基线的观测历元是否相同之后，所述方法还包括：如果基线的观测历元相同，则根据已知点的点位速度、观测历元、参考框架历元和在参考框架历元下已知点的坐标计算观测历元中已知点的坐标；根据观测历元中已知点的坐标、已知点的方差-协方差阵、基线向量和基线的方差-协方差阵进行GNSS网平差。3.根据权利要求2所述的方法，其特征在于，所述根据已知点的点位速度、观测历元、参考框架历元和在参考框架历元下已知点的坐标计算观测历元中已知点的坐标，包括：根据计算观测历元已知点的坐标X
P
；其中，t1表示观测历元，t0表示参考框架历元，V
P
表示已知点的点位速度，表示参考框架历元下已知点的坐标。4.根据权利要求1或2所述的方法，其特征在于，根据已知点的速度参数和/或平面坐标，确定已知点的点位速度和方差-协方差阵，包括：判断已知点是否为IGS跟踪站以及是否具有速度场模型；如果已知点是IGS跟踪站，则获取IGS提供的已知点的点位速度和方差-协方差阵；如果已知点非IGS跟踪站且具有速度场模型，则根据速度场模型中IGS跟踪站的站点速度、站点坐标和已知点的平面坐标确定已知点的点位速度和方差-协方差阵；如果已知点非IGS跟踪站且没有速度场模型，则根据IGS跟踪站的站点速度、站点坐标、方差-协方差阵和已知点的平面坐标，计算已知点的点位速度和方差-协方差阵。5.根据权利要求4所述的方法，其特征在于，当已知点非IGS跟踪站且具有速度场模型，根据速度场模型中IGS跟踪站的站点速度、站点坐标和已知点的平面坐标确定已知点的点位速度和方差-协方差阵，包括：获取速度场模型格网的四个格网点的速度，以及横坐标和纵坐标均最小、横坐标和纵坐标均最大的两个对角格网点的平面坐标；
根据计算已知点P的速度V
P
；根据计算已知点P的方差-协方差阵Q
P
；其中，V1、V2、V3、V4分别表示四个格网点P1、P2、P3、P4的坐标，(x1,y1)、(x2,y2)表示两个对角格网点的坐标，(x,y)表示已知点P的平面坐标，分别表示P1、P2、P3、P4的方差-协方差阵，协方差阵，6.根据权利要求4所述的方法，其特征在于，当已知点非IGS跟踪站且没有速度场模型时，根据IGS跟踪站的站点速度、站点坐标、方差-协方差阵和已知点的平面坐标，计算已知点的点位速度和方差-协方差阵，包括：获取IGS跟踪站的站点速度、站点坐标和方差-协方差阵；根据IGS跟踪站的站点坐标和已知点的平面坐标确定IGS跟踪站与已知点的距离，按照距离由近到远的顺序对IGS跟踪站进行排序，确定顺序排在前n个的IGS跟踪站；根据计算已知点P的速度V
P
；根据计算已知点P的方差-协方差阵Q
P
；其中，(x,y)表示已知点P的平面坐标，(x
i
,y
i
)表示第i个IGS跟踪站的站点坐标，V
i
表示第i个IGS跟踪站的站点速度；i＝1,2,...,n，n≥1，Q
i
表示第i个IGS跟踪站的方差-协方差阵。7.根据权利要求1所述的方法，其特征在于，根据已知点的点位速度、观测历元、参考历元和观测历元中已知点的坐标计算在参考历元已知点的坐标；根据基线两端点的速度计算观测历元的基线两端点相对于参考历元中基线两端点的位移量；根据位移量修正对应基线的基线向量；根据基线的方差-协方差阵、基线两端点的方差-协方差阵、观测历元和参考历元确定修正后基线的方差-协方差阵，包括：根据X'
P
＝(t
1-t2)V
P
+X
P
计算参考历元已知点的坐标X'
p
；根据Δb
AB
＝(t
2-t1)(V
B-V
A
)计算位移量Δb
AB
；根据b＝b
AB
+Δb
AB
计算修正后的基线向量b；根据计算修正后基线的方差-协方差阵Q
b
；其中，A、B表示基线的两端点，t2表示参考历元，t1表示观测历元，V
P
表示已知点的点位
速度，X
P
表示观测历元已知点的坐标，V
B
表示B点的速度，V
A
表示A点的速度，b
AB
表示修正前的基线向量，表示修正前的基线的方差-...

【专利技术属性】
技术研发人员：柴军兵，易昌华，殷学涛，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：