基于局部无缝Delaunay三角网反距离加权法构建中国大陆速度场  

Construction of Velocity Field in Chinese Mainland Based on Local Seamless Delaunay Triangulation with Inverse Distance Weighting Method

作　　者：任营营 王解先 王虎[2] 连丽珍 侯阳飞 王永哲[4] REN Yingying;WANG Jiexian;WANG Hu;LIAN Lizhen;HOU Yangfei;WANG Yongzhe(College of Surveying and Geo-Informatics,Tongji University,Shanghai 200092,China;Chinese Academy of Surveying and Mapping,Beijing 100830,China;Shanghai Astronomical Observatory,Chinese Academy of Sciences,Shanghai 200030,China;Institute of Geophysics,China Earthquake Administration,Beijing 100081,China)

机构地区：[1]同济大学测绘与地理信息学院,上海200092 [2]中国测绘科学研究院,北京100830 [3]中国科学院上海天文台,上海200030 [4]中国地震局地球物理研究所,北京100081

出　　处：《武汉大学学报（信息科学版）》2021年第7期1071-1080,共10页Geomatics and Information Science of Wuhan University

基　　金：国家自然科学基金(41874042);国家自然科学基金青年科学基金(11903065)。

摘　　要：现有研究中,中国大陆速度场模型容易出现数据量小、连续性不佳、现势性不强等问题。基于近7年的中国地壳运动观测网络工程连续运行基准站GNSS(global navigation satellite system)观测数据,解算得到高精度陆态网基准站点的点位坐标和速度场,并利用提出的局部无缝Delaunay三角网反距离加权模型构建中国大陆格网速度场。相较现有的NNR-NUVEL1A(no net rotation Nubia velocity 1A)欧拉矢量模型与常用的速度场数值拟合模型,局部无缝Delaunay三角网反距离加权模型的精度最高,其水平速度场拟合精度优于1.5 mm/a;采用陆态网约1 800个区域站进行外部检核,结果表明,东、北方向速度差值的绝对值平均值分别为1.11 mm/a、0.90 mm/a,中误差分别为1.50 mm/a、1.35 mm/a。相较于整体三角网而言,局域三角网在大陆边缘地区的拟合精度更优,水平速度场的插值精度更高,平均可以提高约0.3 mm/a,且不易产生粗差。局部无缝Delaunay三角网反距离加权模型不仅考虑了邻近点的距离和方位信息,还可以刻画出更为精细的局部特征,同时克服了边缘地区整体三角网跨度过大以及二级块体边缘处三角网不连续的缺点。Objectives: With the development of modern millimeter level geodetic technology, velocity field is one of the indispensable basic elements in Geodetic work. Under the background that current velocity models have limitations such as small data, poor continuity and timeliness update, a method based on local seamless Delaunay triangulation with inverse distance weighting is presented to establish the velocity field model in Chinese mainland. Methods: Aiming at the limitations of the current velocity models, we obtain the high-precision coordinates and velocities of each CMONOC(crustal movement observational network of China) station, through processing the long-term CMONOC GNSS(global navigation satellite system)observation data of recent seven years by the GAMIT/GLOBK software. And we construct the velocity field of Chinese mainland grid with the inverse distance weighted model of locally seamless Delaunay triangulation. Compared with that calculated by NNR-NUVEL1 A(no net rotation Nubia velocity 1 A) Euler vector model and some other numerical fitting models, the velocity obtained by the new method in this research has the highest precision. Results: The fitting precision of horizontal velocity field constructed by this new model is better than 1.5 mm/a. By using more than 1 800 regional stations for external checking,it shows that the average of absolute value of velocity differences in East and North directions are 1.11 mm/a and 0.90 mm/a, respectively, and the corresponding mean square errors are 1.50 mm/a and 1.35 mm/a. In continental margin areas, the fitting precision of horizontal velocity filed by the method of local triangulation net is much higher than that from the whole triangulation net model. Furthermore, based on the method of local triangulation net, interpolation precision of horizontal velocity filed is higher, an improvement of about 0.3 mm/a on average, and the gross errors could be largely reduced. Conclusions: The new model,which is based on the method of local seamless Delaunay triangulation

关 键 词：局部无缝Delaunay三角网 反距离加权法 CGCS2000 速度场 数值拟合模型 精度检核 

分 类 号：P228[天文地球—大地测量学与测量工程]