协方差阵非负约束的赫尔默特方差分量估计  

Helmert variance component estimation with non-negative constraint of covariance matrix

作　　者：王乐洋[1,7] 赵雄 许文斌 汪驰升[5] 方楠 谢磊[6] WANG Leyang;ZHAO Xiong;XU Wenbin;WANG Chisheng;FANG Nan;XIE Lei(Faculty of Geomatics, East China University of Technology, Nanchang 330013, China;School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China;School of Geosciences and Info-Physics, Central South University, Changsha 410083, China;Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Ministry of Education, Changsha 410083, China;Key Laboratory of Geo-Environmental Monitoring of Great Bay Area, MNR, Shenzhen University, Shenzhen 518060, China;Department of Land Surveying and Geo-Informatics, The Hong Kong Polytechnic University, Hong Kong 999077, China;Key Laboratory of Mine Environmental Monitoring and Improving around Poyang Lake, Ministry of Natural Resources, Nanchang 330013, China)

机构地区：[1]东华理工大学测绘工程学院,江西南昌330013 [2]武汉大学测绘学院,湖北武汉430079 [3]中南大学地球科学与信息物理学院,湖南长沙410083 [4]有色金属成矿预测与地质环境监测教育部重点实验室,湖南长沙410083 [5]深圳大学自然资源部大湾区地理环境监测重点实验室,广东深圳518060 [6]香港理工大学土地测量及地理资讯学系,中国香港999077 [7]自然资源部环鄱阳湖区域矿山环境监测与治理重点实验室,江西南昌330013

出　　处：《测绘学报》2022年第3期351-360,共10页Acta Geodaetica et Cartographica Sinica

基　　金：国家自然科学基金(41874001,42174011,41664001,41804015);国家重点研发计划(2019YFC1509205)。

摘　　要：在多源数据联合反演同震滑动分布过程中,各类观测数据之间相对权比的确定是当今大地测量与地球物理联合反演领域的研究热点和难点。近年来,赫尔默特方差分量估计法因其可以同时估计虚拟观测数据及实际观测数据之间的比例因子的优势,在联合反演滑动分布领域中得到了广泛应用。然而,利用HVCE法确定多源数据联合反演相对权比过程中易出现负方差现象。针对以上问题,本文提出在HVCE法确定联合反演相对权比过程中添加线性不等式约束,利用带有线性不等式约束的平差模型进行求解(LC-HVCE),解决计算过程中出现负方差问题。为了验证LC-HVCE法的有效性,开展模拟地震试验,结果表明LC-HVCE法可以有效地解决HVCE法在确定联合反演相对权比过程中存在的不足,保证了同震滑动分布反演精度。蒙特卡洛模拟试验及拉奎拉实际地震反演结果进一步表明了LC-HVCE法在确定正则化参数与相对权比过程中较HVCE法的适用性与稳定性。How to determine the relative weight ratio in geodetic data inversion among various datasets is one of the challenging topics.In recent years,Helmert variance component estimation method(HVCE)has been widely used in the field of relative weights of joint inversion slip distributions due to its advantage of simultaneously estimating the scale factors between virtual and actual observations.However,the HVCE method sometimes suffers from the occurrence of negative variance.To overcome these problems,we propose linear inequality constraints based helmert variance component estimation method(LC-HVCE).We carry out synthetic experiment to verify the effectiveness and robustness of the LC-HVCE method and estimate the uncertainty of LC-HVCE method by Monte carlo method.We apply the proposed method to invert for the source parameters of the 2009 Mw6.3 LAquila earthquake(Italy).We find that LC-HVCE method is more applicable and stability than HVCE method in determining relative weight ratio and regularization parameter.

关 键 词：赫尔默特方差分量估计 相对权比 正则化参数 线性不等式约束 

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