综合利用抗差拟合和Kriging逼近的地壳变形分析方法

英文篇名：Crustal deformation analysis based on robust fitting of function model and Kriging statistics
中文刊名：测绘科学
英文刊名：Science of Surveying and Mapping
作者：赵丽华 ; 王庆良 ; 姚光飞
英文作者：ZHAO Li-hua① ; WANG Qin-liang② ; YAO Guang-fei③(①School of Geology Engineering and Geomatics ; Chang'an University ; Xi'an 710054 ; China ; ②Second Monitoring Center ; China Earthquake Administration ; Xi'an 710054 ; China ; ③Northwest Research Institute of Engineering Investigations and Design ; Xi′an 710003 ; China)
中文关键词：地壳形变分析 ; 函数模型拟合 ; Kriging ; 半变异函数 ; 抗差估计
英文关键词：crustal deformation analysis ; function model fitting ; Kriging ; semivariogram function ; robust estimation
出版日期：2010-01-20
机构：长安大学地质工程与测绘工程学院;中国地震局第二监测中心;西北综合勘察设计研究院;
年：2010
期：01
出版单位：测绘科学

摘要

在地壳形变监测中,由于观测点相对较少,形变分析时常出现观测资料缺乏以及分布不均匀的情况,需要对某些未测点进行推估。考虑到同一区域变形一方面呈现一定的趋势性变化,另一方面又表现出局部的、随机的和不规则的特性,本文在基于可靠函数模型拟合的基础上,利用Kriging统计方法对区域未测点的变形量进行推估。顾及到观测数据中可能含有异常误差,不论是函数模型逼近,还是Kriging推估中半变异函数的拟合都将会受到影响,本文采用抗差估计的方法来抵制异常数据的影响。最后根据某实测地震网监测数据,通过实例验证该方法在未测点变形量估算中的有效性。
As the observation information is often insufficient or uneven in crustal deformation analyzing,the deformations of some unobserved points often need to be predicted.The deformation in one area,on one hand,shows a general tendency and on the other hand,is local,anomalous and stochastic,a model synthesized by functional models and statistical methods is proposed.Considering not only functional approximating but also semivariogram function fitting in Kriging statistics are impacted by measurement outliers,robust estimation method is applied in this paper.Finally,a crustal deformation monitoring network is calculated and analyzed,and the result proves the effectivity of the proposed method in this paper.

引文

[1]黄立人,陶本藻,赵承坤.多面函数拟合在地壳垂直运动中的应用[J].测绘学报,1993,22(1):25-32.
    [2]徐绍铨,李振洪,吴云孙.GPS高程拟合系统研究[J].武汉测绘科技大学学报,1999,24(4):336-340.
    [3]李景卫,高建.GPS高程拟合中多面函数及二次曲面函数的比较与分析[J].山东冶金,2006,28(3).
    [4]江在森,马宗晋,等.GPS初步结果揭示中国大陆水平应变场与构造变形[J].地球物理学报,2003,46(3).
    [5]江在森,牛安福,等.活动断裂带构造变形定量分析[J].地震学报,2005,27(6):610-619.
    [6]张希,江在森.考虑区域构造特征的最小二乘配置方法初步研究[J].中国地震,2001,17(4):403-407.
    [7]张希,江在森,等.借助最小二乘配置整体解算地壳视应变场[J].地壳形变与地震,1998,18(2).
    [8]杨元喜,刘念.重力异常的一种逼近方法[J].测绘学报,2001,30(3):192-196.
    [9]黄维彬.近代平差理论及其应用[M].北京:解放军出版社,1992:140-150.
    [10]孙付平,赵铭,宁津生,等.基于VLBI,SLR和GPS实测数据的现实板块运动模型[J].解放军测绘学院学报,1998,15(4):251-254.
    [11]杨元喜.抗差估计理论及其应用[M].北京:八一出版社,1993.
    [12]杨元喜.等价权原理——参数平差模型的抗差最小二乘解[J].测绘通报,1994,(6):33-35.
    [13]NOEL A,CRESSIE C.Statistics for Special Data[M].New York:John Wiley and Sons Inc,1993.