选取最小曲率、克里格、改进Shepard、反距离加权和径向基函数等5种网格化数学模型，对小江断裂地磁总强度加密区岩石圈磁场数据进行数据网格化，采用均方根预测误差和插值数据残差均方根等评价指标对网格化结果进行评价，结果表明，克里格插值与反距离加权插值法的精度最高。进一步比较克里格插值与反距离加权插值法的网格化图形质量，结果显示克里格插值网格化过程中兼顾了数据的平滑性和各实测点与待估点之间的空间位置关系，避免了系统误差，得出克里格插值更适用于岩石圈磁场数据网格化的结论。
The paper compared and discussed five commonly used interpolation methods. The interpolation methods were evaluated using the parameters such as the root mean square prediction error and the interpolation data root mean square for the lithospheric magnetic field data in the fault zone of Xiaojiang. The results showed that the Kriging interpolation method and inverse distance interpolation method give the highest accuracy. Further, the grid graphics of the two interpolation methods were evaluated. The result showed that Kriging interpolation takes into account the smoothness of the data and the spatial position relationship between the measured points and the points to be estimated during the Gridding process, avoiding systematic errors. It can be concluded that Kriging interpolation is more suitable for Gridding of lithospheric magnetic field data.
2021,42(3): 137-143 收稿日期：2020-11-14
DOI：10.3969/j.issn.1003-3246.2021.03.016
基金项目：2015年地震行业科研专项：中国综合地球物理场观测——大华北地区（项目编号：201508009）协作任务：大华北地区西南缘地磁观测研究
作者简介：张忠龙,男,硕士,工程师,国家注册测绘师,主要从事地震地磁数据处理、高精度GPS数据反演、GIS教学等工作
参考文献：
陈欢欢, 李星, 丁文秀. Surfer8.0等值线绘制中的十二种插值方法[J]. 工程地球物理学报, 2007, 4(1):52-57.
刁鑫鹏, 吴侃. 基于Kriging算法与曲面拟合的三维激光扫描点云数据插值研究[J]. 大地测量与地球动力学, 2012, 32(1):110-113.
丁鉴海, 卢振业, 黄雪香, 等. 地震地磁学[M]. 北京:地震出版社, 1994:97-98.
管泽霖, 李建成, 晁定波, 等. WZD94中国重力大地水准面研究[J]. 武汉测绘科技大学学报, 1994, 19(4):292-297.
李振海, 汪海洪. 重力数据格网化方法比较[J]. 大地测量与地球动力学, 2010, 30(1):140-144.
杨元德, 熊云琪, 王泽民, 等. 几种插值方法在构建南极DEM中的比较[J]. 大地测量与地球动力学, 2013, 33(5):63-66.
郑晓月, 康锟鹏. 一种基于加权平均的曲面拟合算法在电法勘探中的应用[J]. 河南理工大学学报(自然科学版), 2008, 27(1):106-108.
中国科学院地质研究所. 中国地震地质概论[M]. 北京:科学出版社, 1997:53-54.
Pan Guocheng. Geostatistical analysis of the structure of the theory and methods[J]. Global Geology, 1997, 16(3):70-82.
Smith W H F and Wessel P. Gridding with continuous curvature splines in tension[J]. Geophysics, 1990, 55(3):293-305.
