刘亚静 ，赵秀梅 ，2，甘德清 ，2

（1.华北理工大学 矿业工程学院，河北 唐山 063210；2.河北省矿业开发与安全工程实验室，河北 唐山 063210）

0 引言

我国钨资源储量和生产量都居于世界第一，随着对钨需求量的增加，大规模的矿山开采给生态环境带来巨大的影响[1-2]。为了保障矿区的可持续发展，准确模拟矿区地下岩层的分布形态[3-4]，李海港等[5]学者对钨矿进行三维建模可视化。由于离散钻探数据建立的曲面不能准确反映岩层的分布，需要确定精确的空间插值方法，利用最优插值方法，对矿区金属元素进行储量计算。

研究基于某矿区钻孔数据，使用GIS空间地质统计学和相关软件，分析数据特征[6-8]，通过交叉验证方法，得到最优插值方法，为矿体储量计算提供了依据。以期为该矿山地理信息领域的进一步研究提供参考。

1 研究区数据基础

研究矿井样本点经整理共110个样本点数据，将结果整理为EXCEL表格数据。为了将EXCEL数据转换为属性特征表，使用ArcGIS的添加x，y数据工具进行转换，利用ArcGIS定义投影工具，建立相应的数据系统[9-10]，得到钻孔的空间位置分布，如图1所示。

图1 矿区点分布Fig.1 Distribution of the mining area

2 矿区钻孔数据空间特征分析

2.1 地质统计学基础

克里金插值法又称空间插值，是一种基于变量理论和结构分析，对有限区域内的区域变量进行无偏最优估计的方法。采用已知样本点之间的距离和空间方位来反映空间相关性，对未知点进行线性无偏最优预计[11]。

克里金法插值公式如式（1）所示。

其中：λi为第i个样本点的未知权重；……