-
中国地质灾害频发,北起内蒙古,南至海南、广东,东起浙江、山东,西至西藏、新疆,超过22个省份不同程度面临着滑坡灾害的威胁,给人民生命财产安全造成了严重危害[1]。为深入了解滑坡地质灾害变形的内在规律,需要对滑坡灾害的位移形变进行有效监测。全球导航卫星系统(global navigation satellite system,GNSS)技术以全天时观测、高精度、实时性等优势在大坝、桥梁、边坡等形变监测领域得到了广泛应用,是地表形变监测的重要手段之一[2-5]。文献[6]利用全球定位系统(global positioning system,GPS)技术开展了实时滑坡监测研究,结果表明水平定位精度达到15 mm,高程精度达到20 mm;文献[7]利用连续运行参考站(continuously operating reference stations,CORS)观测数据,对意大利北部某区域进行了沉降监测,水平和垂直精度都达到毫米级;文献[8]研究了GPS用于滑坡监测的尺度与方位基准问题,指出监测时需要顾及基准点不稳定导致的基准变化;文献[9]针对单基准站的稳定性问题,提出了多基准站动态实时相对定位变形监测方法,采用法方程域融合处理方法,其定位精度和稳定性均优于坐标域融合处理方法。
在传统的GNSS形变监测方法中,一般选取在待监测区域附近稳定的区域自建基准站作为监测参考点[7-9]。因此针对每个滑坡地质灾害隐患点,均需要单独建设GNSS基准站,并且根据地质灾害隐患点区域大小和地质特性不同,需要建设一个甚至多个GNSS基准站,导致监测系统成本增加,不利于GNSS形变监测技术的大范围应用推广。因此,国内外学者提出利用CORS参考站网进行形变监测[10-13],文献[14]提出了虚拟基站+自适应解算的技术路线,采用服务端固定非差模糊度后生成虚拟观测值,监测端采用实时动态(real-time kinematic,RTK)测量技术动态与静态相结合的解算模式估计监测点位移变化,结果表明可以获得厘米级的定位精度,然而该方法仅通过单基线解算估计确定待监测点的位置坐标,各监测点之间没有组成三角网,无法利用三角网约束信息进行平差提高数据处理精度,造成了有效信息的浪费,不利于形变监测精度的进一步提高。本文提出一种基于虚拟参考站(virtual reference station,VRS)的GNSS滑坡形变监测方法,首先利用现有CORS系统生成VRS,并将VRS与各监测点组成Delaunay三角网,进行基线解算及网平差处理,估计得到待监测点的三维位置坐标,然后确定各监测点的位移变化数据,最后利用中国云南省屏边苗族自治县(下简称屏边县)的实测数据进行算法实验和结果分析。
-
本文提出的基于VRS的GNSS滑坡变形监测方法的数据处理流程如图 1所示,可以分为两个部分:
1)服务端VRS观测数据生成。服务端利用CORS基准站观测数据,首先固定双差宽巷(wide-lane,WL)和窄巷(narrow-lane,NL)整数模糊度,计算高精度站星几何距离;然后解算得到双差对流层和电离层延迟,建立区域对流层、电离层延迟模型;最后内插生成VRS观测值数据。
2)监测端基线解算及网平差处理。首先利用监测站坐标组成Delaunay三角网;然后逐基线进行解算,得到每条基线的向量;其次依据三角形闭合差为零的准则,对基线向量解算结果进行检验;最后将检验通过的基线解算结果进行网平差处理,得到最终的监测点坐标向量。
-
与网络RTK数据处理方法和流程类似,服务端先利用CORS基准站观测数据固定双差模糊度,考虑基准站间距达到20~70 km,一般采用无电离层组合观测值消除电离层延迟的影响[15-16],观测方程为:
]]>∇ Δ P a , b , I F i , j = ∇ Δ ρ a , b i , j + m a , b i , j × Δ T a , b + ε ( ∇ Δ P a , b , I F i , j ) ∇ Δ L a , b , I F i , j = ∇ Δ ρ a , b i , j + m a , b i , j × Δ T a , b + λ I F × ∇ Δ N a , b , I F i , j + ε ( ∇ Δ L a , b , I F i , j ) 式中,
表示双差算子; 、 为卫星编号, 为测站编号; 分别为伪距和载波相位无电离层组合观测值; 为双差站星距离; 为对流层投影系数; 为站间单差天顶对流层延迟; 为无电离层观测值波长; 表示双差无电离层组合模糊度; 、 分别为伪距和载波观测值的噪声。 由于无电离层组合观测值不具有整数特性,故将其分解为WL和NL模糊度进行固定[17-19],其中WL模糊度可采用MW(Melbourne-Wübbena)组合观测值平滑取整进行固定[19],固定方程为:
]]>∇ Δ N a , b , N L i , j = f 1 + f 2 f 1 ∇ Δ N a , b , I F i , j - f 2 f 1 - f 2 ∇ Δ N a , b , W L i , j ∇ Δ N a , b , W L i , j = i n t < ∇ Δ L a , b , M W i , j > 式中,
分别表示双差WL和NL模糊度; 分别表示L1和L2观测值的频率;int表示取整操作; 表示双差MW组合观测值; 表示历元间平滑滤波。 模糊度固定后,首先建立区域大气误差模型,然后根据待监测区域范围和监测站的概略坐标,得到需要虚拟的参考站数量及其坐标,最后利用线性内插模型(linear interpolation model,LIM)[20-21]插值得到大气误差改正模型,同时虚拟生成对应的VRS观测值数据[22]。
-
监测端数据处理时,将服务端数据处理生成的VRS观测值与监测站观测数据组成双差观测值。VRS虽然与待监测点相距较近,一般仅几十到几百米,但考虑VRS观测值实际由距离待监测点几十千米的CORS站计算得到,因此,双差观测值仍需要考虑对流层和电离层残余误差的影响[23-24],因此本文采用L1、L2双频非组合模型,基线解算观测方程为:
]]>∇ Δ P a , b , 1 i , j = ∇ Δ ρ a , b i , j + ∇ Δ I a , b i , j + m a , b i , j × Δ T a , b + ε ( ∇ Δ P a , b , 1 i , j ) ∇ Δ P a , b , 2 i , j = ∇ Δ ρ a , b i , j + μ × ∇ Δ I a , b i , j + m a , b i , j × Δ T a , b + ε ( ∇ Δ P a , b , 2 i , j ) ∇ Δ L a , b , 1 i , j = ∇ Δ ρ a , b i , j - ∇ Δ I a , b i , j + m a , b i , j × Δ T a , b + λ 1 ∇ Δ N a , b , 1 i , j + ε ( ∇ Δ L a , b , 1 i , j ) ∇ Δ L a , b , 2 i , j = ∇ Δ ρ a , b i , j - μ × ∇ Δ I a , b i , j + m a , b i , j × Δ T a , b + λ 2 ∇ Δ N a , b , 2 i , j + ε ( ∇ Δ L a , b , 2 i , j ) 式中,
分别表示L1和L2双差伪距观测值; 分别表示L1和L2双差载波相位观测值; 表示L2相对L1观测值的电离层比例系数; 表示L1观测值上双差电离层延迟; 分别表示L1和L2观测值的波长; 分别表示L1和L2观测值的双差模糊度; 、 、 表示观测值噪声。 为方便数据处理,本文采用静态逐次滤波解算式(3)中的待估计参数,包括基线向量、站间单差天顶对流层延迟、双差斜路径电离层延迟以及双差模糊度参数,其中基线向量在一次解算中认为保持不变,站间单差对流层延迟采用随机游走模型,而双差电离层斜路径延迟参数认为历元间不相关。站间单差对流层延迟和双差电离层斜路径延迟可利用服务端提供的对流层和电离层模型的精度作为先验信息进行约束。基线解算后,为进一步提高数据处理精度,本文将VRS与各监测站组成Delaunay三角网,添加三角网约束信息进行网平差,约束方程为:
∇ Δ X a , b + ∇ Δ X V R S , a + ∇ Δ X b , V R S = 0 ]]>式中,
、 分别表示VRS与Delaunay三角网另外两个监测站a、b的双差基线向量; 表示监测站a、b的双差基线向量。 -
屏边县地处云南省南部、红河哈尼族彝族自治州东南部,受东南海洋暖湿气流的影响,境内湿润多雨,地势北高南低,由北向南倾斜,由于河流的切割,地貌形成了“四河三山六面坡”的结构,高山横亘连绵,V型谷较多,地形极其复杂,是全州地质灾害隐患点较多的县(市)之一。2020年全县共有地质灾害隐患点139处,包括滑坡点94处,崩塌点30处,泥石流沟14条,地面塌陷1处,受地质灾害隐患威胁群众2 762户16 046人,潜在经济损失达5.2亿元。
本文选取其中30处滑坡灾害隐患点进行实验分析,如按常规的一点一基站方式,共需布设30个基站,若要实现全县所有隐患点监测则需要建设139个基站,将造成极大的成本压力。而红河哈尼族彝族自治州已建成了全域126个CORS基准站,平均站间距约20 km,部分CORS站点分布如图 2所示。充分利用已有的CORS基准站生成VRS代替实体基准站,进行滑坡形变监测应用,可大大降低监测成本。此外,VRS利用多基准站虚拟生成,消除了基准点不稳定导致基准变化的影响。
-
电离层延迟和对流层延迟模型精度是影响监测端数据处理精度的主要误差,本文选取基线长度分别为20 km、40 km和80 km以内的观测数据(分别对应图 2中红、黄、蓝三角网),分析电离层延迟和对流层延迟模型精度。利用多个基准站数据固定模糊度后,分离得到区域电离层和对流层延迟模型,采用LIM内插得到虚拟站与基准站之间双差电离层和对流层延迟,同时利用监测站实测数据计算某一卫星对的双差电离层和对流层延迟,比较两者的差异以分析电离层和对流层延迟模型的精度。
图 3和图 4分别为某一卫星对双差电离层和对流层模型的内插值和计算值的时间序列。由图 3和图 4可以看出,对不同的基线长度,经过网络RTK内插模型计算的双差电离层和对流层延迟与实际计算的双差电离层和对流层延迟具有较好的一致性,可以较好地改正电离层和对流层延迟的影响。
图 3 不同基线长度双差电离层建模精度分析
Figure 3. Accuracy Analysis of Double-Difference Ionospheric Modeling with Different Baseline Lengths
图 4 不同基线长度双差对流层建模精度分析
Figure 4. Accuracy Analysis of Double-Difference Tropospheric Modeling with Different Baseline Lengths
图 5和图 6分别为所有卫星对双差电离层和对流层延迟的内插值和计算值的差异概率分布。由图 5和图 6可以看出,双差电离层和对流层延迟内插值和计算值的差异均呈现正态分布的特性,且随着基线长度的增加,内插值和计算值差异逐渐增大,说明随着站间距离的增加,电离层和对流层的内插模型的精度逐渐降低。当站间距小于40 km时,虽然双差电离层和对流层延迟模型的标准差(standard deviation,STD)仅能达到厘米级,但是差异的均值小于1 cm。而滑坡形变监测数据处理一般采用事后解算模式,通过多历元观测数据可以有效削弱内插模型厘米级噪声的影响,实现高精度数据处理。
-
考虑双差电离层和对流层延迟模型的内插精度影响,本文将CORS基准站站间距设置为小于40 km,服务端数据处理对应生成30个VRS观测数据,与90个监测站数据进行事后模拟,通过基线解算及网平差处理获得监测站的三维坐标变化。根据30个测站一周的观测数据,利用不同观测时长(2 h、4 h、6 h、8 h、10 h、12 h)数据的解算结果,统计N、E、U方向的均方根(root mean square,RMS)变化,如图 7所示。
图 7 不同观测时长数据的解算结果精度
Figure 7. Accuracy of Data Processing Results with Different Observation Duration
由图 7可以看出,数据处理解算结果的精度与观测时段长短密切相关,2 h时长观测数据的解算结果仅能达到厘米级精度,在95%置信区间内,U方向的RMS可以达到25 mm;随着观测时长增加,解算精度逐渐提高,当观测时长达到8 h时,在95%置信区间内,N、E方向的RMS控制在5 mm以内、U方向的RMS控制在10 mm以内。
图 8为网平差前后所有基线的RMS均值。由图 8可以看出,经过网平差处理,基线解算结果在N、E、U方向的RMS均值分别从3.2 mm、3.8 mm和5.8 mm降至2.9 mm、3.0 mm和5.1 mm,精度分别提高了9%、14%和12%。选取某一监测站连续15天的8 h观测数据,采用本文方法和传统方法处理结果的时间序列如图 9所示。由图 9可以看出,传统形变监测方法由于基准站与监测站相距较近,其解算结果在N、E、U方向的精度分别达到0.51 mm、0.54 mm和1.26 mm;本文方法的解算结果在N、E、U方向的精度分别为3 mm、3 mm和6 mm,虽然解算精度不及传统方法,但也能够满足毫米级的形变监测应用需求。
-
本文提出了一种基于VRS的GNSS变形监测方法,直接利用现有CORS观测数据虚拟生成VRS,代替实体参考站实现高精度变形监测。利用中国云南省屏边县约30处滑坡监测点进行实验分析,结果表明,当CORS基准站间距小于40 km时,双差对流层和电离层延迟模型内插精度优于10 mm,利用观测时长大于8 h的数据,可实现平面精度优于5 mm、高程精度优于10 mm的形变监测。
Landslide Deformation Monitoring Method and Performance Analysis Based on GNSS Virtual Reference Station
-
摘要: 全球导航卫星系统(global navigation satellite system,GNSS)技术在大坝、桥梁、滑坡等形变监测领域应用广泛。针对传统监测方法需要架设基准站、导致监测成本增加的问题,提出一种基于虚拟参考站的GNSS变形监测方法,利用连续运行参考站(continuously operating reference stations,CORS)观测数据,生成虚拟参考站(virtual reference station,VRS)代替实体参考站,实现高精度变形监测。利用中国云南省屏边苗族自治县约30处滑坡监测点进行实验分析,结果表明,当CORS基准站间距小于40 km时,双差对流层和电离层延迟模型内插精度优于10 mm,利用观测时长大于8 h的数据,可以实现平面精度优于5 mm、高程精度优于10 mm的形变监测。Abstract:
Objectives Global navigation satellite system (GNSS) technology is widely used in deformation monitoring such as dams, bridges, and landslides. In traditional monitoring methods, the reference stations need to be built, which leads to the increased monitoring cost. Methods This paper proposes a GNSS deformation monitoring method based on virtual reference station (VRS), which uses the existing continuously operating reference stations (CORS) observation data to generate VRS instead of physical reference station to achieve high accuracy deformation monitoring. Firstly, the existing CORS observation data are used to obtain the atmospheric error correction models and generate VRS. Secondly, the baselines between VRS and each monitoring point are calculated. In order to further improve the accuracy of data processing, Delaunay triangulation networks are formed, and the baseline vector solution results are checked according to the criterion that the triangle closure error is zero. Finally, the checked baseline solution results are adjusted to obtain the final three-dimensional position coordinates of the monitoring points. Results This paper uses about 30 monitoring points in Pingbian Miao Autonomous County, Yunnan Province for experimental analysis. The results show that the interpolation accuracis of the double-difference tropospheric and ionospheric delay are better than 10 mm if the distance between CORS reference stations is less than 40 km. The observation data longer than 8 hours can realize deformation monitoring with horizontal accuracy better than 5 mm and vertical accuracy better than 10 mm. Conclusions The proposed method can realize the practical application of landslide monitoring. In the future, we can study how to improve the atmospheric error model and further improve the accuracy of data processing. -
-
[1] 殷跃平. 中国地质灾害减灾战略初步研究[J]. 中国地质灾害与防治学报, 2004, 15(2): 1-8 doi: 10.3969/j.issn.1003-8035.2004.02.001 Yin Yueping. Initial Study on the Hazard-Relief Strategy of Geological Hazard in China[J]. The Chinese Journal of Geological Hazard and Control, 2004, 15(2): 1-8 doi: 10.3969/j.issn.1003-8035.2004.02.001 [2] 黄声享, 尹晖, 蒋征. 变形监测数据处理[M]. 武汉: 武汉大学出版社, 2003 Huang Shengxiang, Yin Hui, Jiang Zheng. Data Processing of Deformity Monitoring[M]. Wuhan: Wuhan University Press, 2003 [3] 谭行, 蒋健, 陈炜, 等. 分布式GNSS地质灾害实时监测系统设计[J]. 测绘地理信息, 2020, 45(6): 118-123 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202006025.htm Tan Hang, Jiang Jian, Chen Wei, et al. Design of Distributed GNSS Geological Hazards Real-Time Monitoring System[J]. Journal of Geomatics, 2020, 45(6): 118-123 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202006025.htm [4] Zumberge J F, Heflin M B, Jefferson D C, et al. Precise Point Positioning for the Efficient and Robust Analysis of GPS Data from Large Networks[J]. Journal of Geophysical Research: Solid Earth, 1997, 102(B3): 5005-5017 doi: 10.1029/96JB03860 [5] 王利, 张勤, 范丽红, 等. 北斗/GPS融合静态相对定位用于高精度地面沉降监测的试验与结果分析[J]. 工程地质学报, 2015, 23(1): 119-125 https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201501022.htm Wang Li, Zhang Qin, Fan Lihong, et al. Experiment and Results of High Precision Land Subsidence Monitoring Using Fused BDS/GPS Data and Static Relative Positioning[J]. Journal of Engineering Geology, 2015, 23(1): 119-125 https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201501022.htm [6] Gili J A, Corominas J, Rius J. Using Global Positioning System Techniques in Landslide Monitoring [J]. Engineering Geology, 2000, 55(3): 167-192 doi: 10.1016/S0013-7952(99)00127-1 [7] Rawat M S, Joshi V, Rawat B S, et al. Landslide Movement Monitoring Using GPS Technology: A Case Study of Bakthang Landslide, Gangtok, East Sikkim, India[J]. Journal of Development and Agricultural Economics, 2011, 3: 194-200 [8] 张勤. GPS监测滑坡形变的基准研究[J]. 西安工程学院学报, 2001, 23(4): 69-71 https://www.cnki.com.cn/Article/CJFDTOTAL-XAGX200104014.htm Zhang Qin. The Datum Research of the Slope Deformation Monitor by GPS[J]. Journal of Xi'an Engineering University, 2001, 23(4): 69-71 https://www.cnki.com.cn/Article/CJFDTOTAL-XAGX200104014.htm [9] 唐洪林. 多基准站Multi-GNSS实时动态相对定位及在变形监测中的应用[D]. 西安: 长安大学, 2020 Tang Honglin. Multi-GNSS Real Time Kinematic Relative Positioning of Multiple Base Stations and Its Application in Deformation Monitoring[D]. Xi'an: Chang'an University, 2020 [10] D'Oreye N, Arjona A, Bagalwa M, et al. A Permanent Geodetic GNSS Network to Monitor Ground Deformation in the Virunga Volcanic Province[C]// Active Volcanism & Continental Rifting (AVCOR), Gisenyi, Rwanda, 2013 [11] Yusup A, Othman R, Musliman I A, et al. ISKANDARnet CORS Network Integrity Monitoring [J]. Jurnal Teknologi, 2014, 71(4): 11-19 [12] 隆华平. GPS连续运行参考站网络在变形监测中的应用研究[D]. 武汉: 中国地质大学(武汉), 2005 Long Huaping. Application of GPS Continuous Operation Reference Station Network in Deformation Monitoring[D]. Wuhan: China University of Geosciences (Wuhan), 2005 [13] 刘新华, 尚俊娜, 施浒立. 基于GNSS-RTK技术的双层滑坡监测系统设计与实现[J]. 电子技术应用, 2020, 46(11): 1-7 https://www.cnki.com.cn/Article/CJFDTOTAL-DZJY202011002.htm Liu Xinhua, Shang Junna, Shi Huli. Design and Implementation of a Double-Layer Landslide Monitoring System Based on GNSS-RTK Technology[J]. Application of Electronic Technique, 2020, 46(11): 1- 7 https://www.cnki.com.cn/Article/CJFDTOTAL-DZJY202011002.htm [14] 张鸣之, 湛兵, 赵文祎, 等. 基于虚拟参考站技术的滑坡高精度位移监测系统设计与实践[J]. 中国地质灾害与防治学报, 2020, 31(6): 54-59 https://www.cnki.com.cn/Article/CJFDTOTAL-ZGDH202006007.htm Zhang Mingzhi, Zhan Bing, Zhao Wenyi, et al. Design and Practice of High Precision Landslide Displacement Monitoring System Based on VRS[J]. The Chinese Journal of Geological Hazard and Control, 2020, 31(6): 54-59 https://www.cnki.com.cn/Article/CJFDTOTAL-ZGDH202006007.htm [15] 省天琛, 张生鹏, 赵利江, 等. 低成本单频精密单点定位在城市导航定位中的应用[J]. 测绘地理信息, 2022, 47(2): 23-25 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202202006.htm Sheng Tianchen, Zhang Shengpeng, Zhao Lijiang, et al. Application of Low-Cost Single-Frequency Precise Point Positioning in Urban Navigation and Positioning[J]. Journal of Geomatics, 2022, 47 (2): 23-25 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202202006.htm [16] Yuan Y B, Mi X L, Zhang B C. Initial Assessment of Single- and Dual-Frequency BDS-3 RTK Positioning[J]. Satellite Navigation, 2020, 1: 31 [17] Kozlov D, Tkachenko M. Centimeter-Level, RealTime Kinematic Positioning with GPS+GLONASS C/A Receivers[J]. Navigation, 1998, 45(2): 137-147 [18] Ge M, Gendt G, Rothacher M, et al. Resolution of GPS Carrier-Phase Ambiguities in Precise Point Positioning(PPP)with Daily Observations[J]. Journal of Geodesy, 2008, 82(7): 389-399 [19] Geng J H, Shi C, Ge M R, et al. Improving the Estimation of Fractional-Cycle Biases for Ambiguity Resolution in Precise Point Positioning[J]. Journal of Geodesy, 2012, 86(8): 579-589 [20] 赵传华, 秘金钟, 党亚民, 等. 电离层内插模型的分析[J]. 测绘科学技术学报, 2013, 30(2): 140-143 https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC201302009.htm Zhao Chuanhua, Bei Jinzhong, Dang Yamin, et al. Analysis on Ionosphere Interpolation Models[J]. Journal of Geomatics Science and Technology, 2013, 30(2): 140-143 https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC201302009.htm [21] 宋伟伟, 何成鹏, 辜声峰. 不同纬度区域电离层增强PPP-RTK性能分析[J]. 武汉大学学报·信息科学版, 2021, 46(12): 1832-1842 doi: 10.13203/j.whugis20210243 Song Weiwei, He Chengpeng, Gu Shengfeng. Performance Analysis of Ionospheric Enhanced PPPRTK in Different Latitudes[J]. Geomatics and Information Science of Wuhan University, 2021, 46 (12): 1832-1842 doi: 10.13203/j.whugis20210243 [22] 唐卫明. 大范围长距离GNSS网络RTK技术研究及软件实现[D]. 武汉: 武汉大学, 2006 Tang Weiming. Research on Techniques of Large Area and Long Range GNSS Network RTK and Developing Network RTK Software[D]. Wuhan: Wuhan University, 2006 [23] Uzel T, Eren K, Gulal E, et al. Monitoring the Tectonic Plate Movements in Turkey Based on the National Continuous GNSS Network[J]. Arabian Journal of Geosciences, 2013, 6(9): 3573-3580 [24] Li B F, Teunissen P J G. GNSS Antenna ArrayAided CORS Ambiguity Resolution[J]. Journal of Geodesy, 2014, 88(4): 363-376 -