综合利用抗差拟合和Kriging逼近的地壳变形分析方法Crustal deformation analysis based on robust fitting of function model and Kriging statistics
赵丽华,王庆良,姚光飞
摘要(Abstract):
在地壳形变监测中,由于观测点相对较少,形变分析时常出现观测资料缺乏以及分布不均匀的情况,需要对某些未测点进行推估。考虑到同一区域变形一方面呈现一定的趋势性变化,另一方面又表现出局部的、随机的和不规则的特性,本文在基于可靠函数模型拟合的基础上,利用Kriging统计方法对区域未测点的变形量进行推估。顾及到观测数据中可能含有异常误差,不论是函数模型逼近,还是Kriging推估中半变异函数的拟合都将会受到影响,本文采用抗差估计的方法来抵制异常数据的影响。最后根据某实测地震网监测数据,通过实例验证该方法在未测点变形量估算中的有效性。
关键词(KeyWords): 地壳形变分析;函数模型拟合;Kriging;半变异函数;抗差估计
基金项目(Foundation): 国家自然科学基金资助项目(40774001,40841021,40902081);; 国家863项目(2007AA12Z331)
作者(Author): 赵丽华,王庆良,姚光飞
DOI: 10.16251/j.cnki.1009-2307.2010.01.010
参考文献(References):
- [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.