CN108303713A - 一种参考站稳定性的自动检验方法 - Google Patents
一种参考站稳定性的自动检验方法 Download PDFInfo
- Publication number
- CN108303713A CN108303713A CN201711371727.0A CN201711371727A CN108303713A CN 108303713 A CN108303713 A CN 108303713A CN 201711371727 A CN201711371727 A CN 201711371727A CN 108303713 A CN108303713 A CN 108303713A
- Authority
- CN
- China
- Prior art keywords
- station
- stability
- current
- stations
- euler
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000033001 locomotion Effects 0.000 claims abstract description 26
- 238000009826 distribution Methods 0.000 claims abstract description 14
- 239000013598 vector Substances 0.000 claims description 23
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 238000012544 monitoring process Methods 0.000 claims description 9
- 230000007774 longterm Effects 0.000 claims description 3
- 238000007689 inspection Methods 0.000 abstract description 4
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000010276 construction Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000010998 test method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 238000013112 stability test Methods 0.000 description 1
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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明为一种参考站稳定性的自动检验方法，参考站的稳定性与否直接决定着参考框架的精度和变形监测成果的可靠性。本发明提出了一种基于统计假设检验的欧拉运动模型，利用模型残差构造服从一定自由度为分布的标准统计量,通过置信区间来判断台站的稳定性，从而实现了参考站网稳定性的自动检验。该方法检验结果精确可靠、自动化程度高，可用于大规模观测台网的稳定性检测和稳定区域的精细划分。同时，为高精度参考框架与基准的建立提供了重要理论依据和技术支撑。
Description
技术领域
本发明涉及参考站网稳定技术领域，具体涉及一种参考站稳定性 的自动检验方法。
背景技术
参考框架是由空间大地测量众多台站三维坐标和速度场来维持 和实现。在控制网数据处理中，常常需要若干参考站作为框架站进行 约束，从而得到观测站的坐标。然而，由于个别参考站不可避免的观 测误差或可能受到的局部干扰或局部地壳形变等不稳定因素的影响, 从而导致同一个区域的多个台站会表现出不协调的位移，产生不自洽 性，使整个观测网平差结果错误。此外，大型建筑物变形监测(如桥 梁、大坝等)和地质灾害监测也需要以稳定的参考站为基准，看监测 点相对于参考站的变化量，进而评定该区域的地质稳定性。当这些参 考站点坐标的精度较低或含有粗差时，不仅会造成单位权方差估计值 不准确，并会致使监测获得的精密成果发生严重扭曲。可见，参考站 的稳定性与否直接影响着地球参考框架和变形监测成果的精度与可 靠性。因此，参考站稳定性分析检验在测量中有着至关重要的作用。 为了能够发现不稳定的参考站，多数是采用经验方法和先验知识,通 常布设多个参考站构成一个参考网，通过定期对参考网的复测、比较 两期或多期的坐标位移量以经验值进行判定相应的参考站是否稳定， 这样往往缺乏严密的理论依据，人为干预很大，分析结果缺乏可靠性。 而且当参考站受干扰因素较小时，或者多个台站都存在不稳定时，直 接从观测资料中逐一筛选、判断参考站稳定性具有很大难度。尤其是， 近年来，随着各区域观测网络的逐渐建设，台站数量在急剧增加，观 测资料日益丰富，迫切需要构建一种敏感度较高的模型，对大规模的 台站观测网络进行稳定性分析，根据分析结果选择一种平差模型、全 面甄别不稳定参考站、建立一个对变形分析比较有利的参考基准，从而提高参考站稳定性检验水平和自动化程度。
发明内容
为了解决上述技术问题，本发明提供了一种参考站稳定性的自动 检验方法，可高效地根据选取的稳定的台站可建立高精度参考框架基 准，确保了控制测量及变形监测成果的可靠性。
本发明要解决的技术问题是通过以下技术方案来实现的：一种参 考站稳定性的自动检验方法，其特征是：
步骤S1：建立所述参考站网运动模型，其中所述参考站包括多 个台站，如变形监测站；
步骤S2：根据欧拉运动理论，建立所述运动模型的欧拉矢量函 数；
步骤S3：采用最小二乘法求解所述欧拉矢量函数，并根据所述 欧拉矢量函数的解得到所述运动模型的模型残差，以及所述模型残差 的概率分布；
步骤S4：构造当前台站服从一定自由度的标准统计量，设置置 信区间，计算所述当前台站的所述台站统计量；
步骤S5：根据台站统计量是否落入置信区间检测所述台站统计 量的事件概率以判断所述当前台站的台站稳定性；
步骤S6：直至完成所有所述台站的稳定性判断。
本发明要解决的技术问题是通过以下技术方案进一步来实现的： 所述步骤5具体包括以下步骤：
步骤S51：检测所述台站统计量的事件概率，当所述事件概率为 正常概率时：判断当前台站为稳定台站；对下一台站，重复所述步骤 S4至所述步骤S5；
步骤S52：检测所述台站统计量的事件概率，当所述事件概率为 小概率时：判断所述当前台站为不稳定台站，剔除所述当前台站；对 下一台站，重复所述步骤S1至所述步骤S5。
本发明要解决的技术问题是通过以下技术方案进一步来实现的：
V＝Vh+Vq+ξ (1)
式中，V是空间大地测量观测三维速度矢量，可表示为Vh为 测站长期运动速度，Vq为局部地壳形变引起的速度，ξ为观测噪声引 起的速度。
本发明要解决的技术问题是通过以下技术方案进一步来实现的： 所述步骤S2中：设置欧拉运动理论公式为：
式中，r为地球半径，λ、φ分别为所述参考站的经度、纬度；其 中，为三维观测量，为已知量，用L表示；为所述台 站位置矢量，为已知量，用A表示；为欧拉三维矢量，为未知量， 用表示；即
本发明要解决的技术问题是通过以下技术方案进一步来实现的： 所述步骤S3中：
则对于所述参考站，利用最小二乘法求解算欧拉矢量函数为：
式中，P表示所述参考站的权重值；
进一步得到所述运动模型的所述模型残差：
V＝(A(ATPA)-1ATP-I)L (4)
当所述参考站由于观测误差、受到的局部干扰或局部变形,整个 观测站网会表现不协调的位移，出现不自洽性。令所述不协调的位移 为ε，则式(4)变为
此时，函数服从具有自由度的t分布。
本发明要解决的技术问题是通过以下技术方案进一步来实现的： 根据所述t分布，设置稳定限值其中所述稳定限值满足：
或
本发明要解决的技术问题是通过以下技术方案进一步来实现的： 所述步骤S4中：针对每一个所述台站，构造并计算所述台站统计量：
式中，σ0为所述运动模型中单位权中误差，为第k台站残差的协 因。
本发明要解决的技术问题是通过以下技术方案进一步来实现的： 所述步骤S51具体为：若所述当前台站的所述台站统计量满足 则所述当前台站为稳定台站；对下一台站，重复所述步 骤S4至所述步骤S5；
所述步骤S52：若所述当前台站的所述台站统计量满足则所 述当前台站为不稳定台站，剔除所述当前台站；对下一台站，重复所 述步骤S1至所述步骤S5。
本发明要解决的技术问题是通过以下技术方案进一步来实现的： 步骤S6之后，具有步骤S7：当完成所有所述台站的所述台站稳定性 的判断后，由于已经剔除所述不稳定台站，则所述参考站为稳定参考 站。
本发明有益技术效果为：
采用基于欧拉运动的参考站模型，利用t分布统计概率检验每个 台站的稳定性，提出不稳定性台站，从而实现了每个台站的稳定性检 验，进而实现了参考站网的稳定性的自动检验。该方法检验结果精确 可靠、自动化程度高，可用于大规模观测台网的稳定性检测和稳定区 域的精细划分，同时，为高精度参考框架与基准的建立提供了重要理 论依据和技术支撑。
附图说明
图1为本发明一具体实施例流程图；
图2为本发明一具体实施例采用的欧拉运动矢量示意图；
图3为本发明一具体实施例中的不稳定台站剔除前残差分布；
图4为本发明一具体实施例中的不稳定台站剔除后残差分布；
图5为本发明一具体实施例中的单位权中误差随检验次数的变 化。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明提出了一种参考站稳定性的自动检验方法，
首先，设置台站的速度：
V＝Vh+Vq+ξ (1)
式中，V是空间大地测量观测三维速度矢量，可表示为Vh为 测站长期运动速度，Vq为局部地壳形变引起的速度，ξ为观测噪声引 起的速度。
其次进行如下步骤：
步骤S1：建立所述参考站的系统模型，其中所述参考站包括多 个台站，如变形监测站；
步骤S2：根据欧拉运动理论，建立所述运动模型的欧拉矢量函 数；具体的，设置欧拉运动理论公式为：
式中，r为地球半径，λ、φ分别为所述参考站的经度、纬度；其 中，为三维观测量，为已知量，用L表示；为所述台 站位置矢量，为已知量，用A表示；为欧拉三维矢量，为未知量， 用表示；即
步骤S3：采用最小二乘法求解所述欧拉矢量函数，并根据所述 欧拉矢量函数的解得到所述运动模型的模型残差，以及所述模型残差 的概率分布；具体的，则对于所述参考站，利用最小二乘法求解算欧 拉矢量函数为：
式中，P表示所述参考站的权重值；
进一步得到所述运动模型的所述模型残差：
V＝(A(ATPA)-1ATP-I)L (4)
当所述参考站由于观测误差、受到的局部干扰或局部变形,整个 基准站网会表现出不协调的位移，令所述不协调的位移为ε，则式(4) 变为
此时，函数服从具有自由度的t分布。
进一步，根据所述t分布，设置稳定限值其中所述稳定限值满足：
或
其中α可以根据实际情况进行设置，比如α为0.05。
步骤S4：构造当前台站的台站统计量，计算所述当前台站的所 述台站统计量；具体的，针对每一个所述台站，构造并计算所述台站 统计量：
式中，σ0为所述运动模型中单位权中误差，为第k台站残差 的协因。
步骤S5：查表检测所述台站统计量的事件概率以判断所述当前 台站的台站稳定性；具体的，步骤S51：查表检测所述台站统计量的 事件概率，当所述事件概率为正常概率时：判断当前台站为稳定台站； 对下一台站，重复所述步骤S4至所述步骤S5；
更进一步，所述步骤S51具体为：若所述当前台站的所述台站统计量满足则所述当前台站为稳定台站；对下一台站，重 复所述步骤S4至所述步骤S5；
所述步骤S51：若所述当前台站的所述台站统计量满足则所述当前台站为不稳定台站，剔除所述当前台站；对下一台站，重 复所述步骤S1至所述步骤S5。
所述步骤S52：查表检测所述台站统计量的事件概率，当所述事 件概率为小概率时：判断所述当前台站为不稳定台站，剔除所述当前 台站；对下一台站，重复所述步骤S1至所述步骤S5
步骤S6：直至完成所有所述台站的所述台站稳定性判断。
步骤S7：当完成所有所述台站的所述台站稳定性的判断后，由 于已经剔除所述不稳定台站，则所述参考站为稳定参考站。
具体结合一实例分析。
采用179个台站组成的参考站，建立参考站系统，采用本发明所 提出的一种参考站稳定性的自动检验方法。经过9次试验，依次剔除 不通过的台站，最终确定69个测站作为稳定台站。其中，图3、图4 依次表示本实施例中不稳定台站剔除前残差分布、不稳定台站剔除后 残差分布。由图3、图4可以看出，不稳定台站剔除前N分量模型残 差值最大达到70mm，E分量模型残差最大达到25mm。剔除后N分 量模型残差基本分布在±0.4mm之内，E分量模型残差分布在±0.7mm之 内。可见，通过假设检验剔除不稳定台站后使模型拟合度得到了明显 改善，效果显著。
Claims (9)
1.一种参考站稳定性的自动检验方法，其特征在于：
步骤S1：建立所述参考站网运动模型，其中所述参考站包括多个台站，如变形监测站；
步骤S2：根据欧拉运动理论，建立所述运动模型的欧拉矢量函数；
步骤S3：采用最小二乘法求解所述欧拉矢量函数，并根据所述欧拉矢量函数的解得到所述运动模型的模型残差，以及所述模型残差的概率分布；
步骤S4：构造当前台站服从一定自由度的标准统计量，设置置信区间，计算所述当前台站的所述台站统计量；
步骤S5：根据台站统计量是否落入置信区间检测所述台站统计量的事件概率以判断所述当前台站的台站稳定性；
步骤S6：直至完成所有所述台站的稳定性判断。
2.根据权利要求1所述一种参考站稳定性的自动检验方法，其特征在于：所述步骤S5具体包括以下步骤：
步骤S51：检测所述台站统计量的事件概率，当所述事件概率为正常概率时：判断当前台站为稳定台站；对下一台站，重复所述步骤S4至所述步骤S5；
步骤S52：检测所述台站统计量的事件概率，当所述事件概率为小概率时：判断所述当前台站为不稳定台站，剔除所述当前台站；对下一台站，重复所述步骤S1至所述步骤S5。
3.根据权利要求1所述一种参考站稳定性的自动检验方法，其特征在于：所述步骤S1中：所述参考站速度可表示为：
V＝Vh+Vq+ξ (1)
式中，V是空间大地测量观测三维速度矢量或位移矢量，可表示为Vh为测站长期运动速度或位移量，Vq为局部地壳形变引起的速度或位移，ξ为观测噪声引起的速度或位移。
4.根据权利要求2所述一种参考站稳定性的自动检验方法，其特征在于：所述步骤S2中：设置欧拉运动理论公式为：
式中，r为地球半径，λ、φ分别为所述参考站的经度、纬度；其中，为三维观测量，为已知量，用L表示；为所述台站位置矢量，为已知量，用A表示；为欧拉三维矢量，为未知量，用表示；即
5.根据权利要求4所述一种参考站稳定性的自动检验方法，其特征在于：所述步骤S3中：
则对于所述参考站，利用最小二乘法求解算欧拉矢量函数为：
式中，P表示所述参考站的权重值；
进一步得到所述运动模型的所述模型残差：
V＝(A(ATPA)-1ATP-I)L (4)
当所述参考站由于观测误差、受到的局部干扰或局部变形,整个观测网的台站会表现不协调的位移，令所述不协调的位移为ε，则式(4)变为
此时，函数服从具有自由度的t分布。
6.根据权利要求5所述一种参考站稳定性的自动检验方法，其特征在于：根据所述t分布，设置稳定限值其中所述稳定限值满足：
或
7.根据权利要求6所述一种参考站稳定性的自动检验方法，其特征在于：所述步骤S4中：针对每一个所述台站，构造并计算所述台站统计量：
式中，σ0为所述运动模型中单位权中误差，即 为第k台站残差的协因数。
8.根据权利要求7所述一种参考站稳定性的自动检验方法，其特征在于：所述步骤S51具体为：若所述当前台站的所述台站统计量满足则所述当前台站为稳定台站；对下一台站，重复所述步骤S4至所述步骤S5；
所述步骤S52：若所述当前台站的所述台站统计量满足则所述当前台站为不稳定台站，剔除所述当前台站；对下一台站，重复所述步骤S1至所述步骤S5。
9.根据权利要求8所述一种参考站稳定性的自动检验方法，其特征在于：步骤S6之后，具有步骤S7：当完成所有所述台站的所述台站稳定性的判断后，由于已经剔除所述不稳定台站，则所述参考站为稳定参考站。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711371727.0A CN108303713A (zh) | 2017-12-19 | 2017-12-19 | 一种参考站稳定性的自动检验方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711371727.0A CN108303713A (zh) | 2017-12-19 | 2017-12-19 | 一种参考站稳定性的自动检验方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108303713A true CN108303713A (zh) | 2018-07-20 |
Family
ID=62870586
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711371727.0A Pending CN108303713A (zh) | 2017-12-19 | 2017-12-19 | 一种参考站稳定性的自动检验方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108303713A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109782275A (zh) * | 2019-03-14 | 2019-05-21 | 中国电建集团成都勘测设计研究院有限公司 | Gnss变形监测的参考点校验系统及方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101498785A (zh) * | 2009-03-16 | 2009-08-05 | 东南大学 | Gnss网络差分定位基准站网构造方法与动态更新方法 |
| CN105806208A (zh) * | 2016-03-11 | 2016-07-27 | 河南理工大学 | 一种基于gnss网形变化的形变异常检测方法 |
-
2017
- 2017-12-19 CN CN201711371727.0A patent/CN108303713A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101498785A (zh) * | 2009-03-16 | 2009-08-05 | 东南大学 | Gnss网络差分定位基准站网构造方法与动态更新方法 |
| CN105806208A (zh) * | 2016-03-11 | 2016-07-27 | 河南理工大学 | 一种基于gnss网形变化的形变异常检测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 徐克科,等: "《基于假设检验的刚性板块台站选取及全球板块运动模型的建立》", 《大地测量与地球动力学》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109782275A (zh) * | 2019-03-14 | 2019-05-21 | 中国电建集团成都勘测设计研究院有限公司 | Gnss变形监测的参考点校验系统及方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106855631B (zh) | 一种转发式gnss动态测量精度测试评估方法 | |
| Castaldo et al. | P‐RANSAC: An Integrity Monitoring Approach for GNSS Signal Degraded Scenario | |
| Teferle et al. | Using continuous GPS and absolute gravity to separate vertical land movements and changes in sea-level at tide-gauges in the UK | |
| CN108075826B (zh) | 一种实现卫星半物理仿真试验星地时间同步的方法及系统 | |
| CN104504247B (zh) | 一种gps双卫星故障的raim方法 | |
| CN105425248B (zh) | 单频gnss相位稳定性监测的高频逐历元相位差方法 | |
| JP2016153726A (ja) | 観測システム | |
| CN115453579A (zh) | 基于北斗ppp-rtk的合成星历a类故障完好性监测方法及装置 | |
| JP2017203650A (ja) | 観測システム | |
| CN109782275B (zh) | Gnss变形监测的参考点校验系统及方法 | |
| CN108303713A (zh) | 一种参考站稳定性的自动检验方法 | |
| CN105204047B (zh) | 一种卫星导航系统中观测量单个粗差的探测与修复方法 | |
| CN114609650A (zh) | 一种基于北斗全链条故障激励的完好性测试方法 | |
| Kowalczyk | The creation of a model of relative vertical crustal movements in the Polish territory on the basis of the data from active geodetic network EUPOS (ASG EUPOS) | |
| Ebeling | Ground-based deformation monitoring | |
| CN116990746A (zh) | 用于无线电监测测向系统及方法 | |
| CN107478209B (zh) | 超高层轴线控制网的检测方法 | |
| Kannemans | The Generalized Extreme Value statistical method to determine the GNSS integrity performance | |
| CN114114332B (zh) | 一种有效探测gnss基准站坐标时间序列不连续点的方法 | |
| Bayrak | Modelling the relationship between water level and vertical displacements on the Yamula Dam, Turkey | |
| CN114966779A (zh) | 一种基于北斗导航卫星监测站的山区峡谷定位方法及系统 | |
| CN112255647A (zh) | 弹道测量信息事后融合处理优化方法 | |
| CN110728007B (zh) | 一种基于模型特征的动态故障诊断方法 | |
| CN111336970A (zh) | 一种cors基准站发生位移的监测方法及系统 | |
| Imtiaz et al. | Design and implementation of receiver autonomous integrity monitoring algorithm on DSP for small UAV applications |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180720 |