张忠龙 倪 喆 赵育飞

1）中国贵州 562400 黔西南民族职业技术学院

2）中国昆明 650224 云南省地震局

0 引言

地震发生受到构造运动控制（中国科学院地质研究所，1997），构造运动伴随着岩石圈磁场的变化，岩石圈磁场及其变化是地球深部物理过程的重要信息来源之一（丁鉴海等，1994）。地磁总强度数据通过日变通化、长期变改正处理后剥离得到岩石圈磁场，岩石圈磁场位场分离后得到地壳深部、浅表和居里面的能量加卸载信息，直接反映了地球内部各深度的物理过程和构造活动能量，比较真实可靠。但获取的地磁总强度数据是离散分布的点位，而实际采集的数据点位有限，难以直观反映岩石圈磁场变化，进而不利于判定岩石圈磁场异常区，需要对岩石圈磁场数据进行网格化处理。网格化数学模型的好坏不但影响到网格化岩石圈磁场数据的精度和可信程度，而且会进一步影响解释处理图件的质量、效果和可靠性。常用网格化数学模型有最小曲率、克里格、改进Shepard、反距离加权和径向基函数等（李振海等，2010）。文中将以上5 种数学模型应用于小江断裂带岩石圈磁场数据网格化，并采用实验方法对其精度进行比较分析。

1 数学模型基本原理

1.1 最小曲率

最小曲率法用公式表示为

1.2 克里格

克里格模型根据所给数据趋势，构造一个只与点间空间距离相关的半方差函数，从而求取各拟合点的权系数，据此进行数值内插（刁鑫鹏等，2012）。

对于区域化变量Z(xi,yi)，假设在n个位……