摘要
目前铜多金属矿的钻孔分析数据多呈条带状分布,并且多金属之间相互影响,采用传统的协同克里格法对该类条带状分布的数据进行插值分析时,易出现端点数据权值偏大的现象,即条带状效应。对此,以协同克里格法为基础,提出了一种铜多金属矿储量估算插值方法,即权值校正的协同克里格插值方法。首先对原始钻孔数据进行统计分析、特异值处理以及样品组合处理,并对产生条带状效应的权值系数进行分析;然后采用拟牛顿法对协同克里格法插值得到的权值进行校正。以2012年西藏甲玛矿区勘探工程的铜多金属矿床钻孔品位数据为例进行储量估算,即以Cu品位为主区域化变量,Ag品位为协同区域化变量,计算了各自的试验变差函数和交差试验变差函数,并对协同克里格法以及权值校正的协同克里格法的插值效果进行了分析。结果表明,权值校正的协同克里格法插值效果较优,对于铜多金属矿的地质属性分析、储量估算具有一定的参考价值。
At present,the drilling analysis data of Cu-polymetallic ore present as a ribbon,with the interaction of multimetal. When using the classical Co-Kriging method to conduct interpolation analysis of the data distribute as a ribbon,it is apt to occur that the weight of endpoint data become too large,called as banded effect. In order to reduce the effect of banded effect in the process of data interpolation,based on the Co-Kriging method,a interpolation method for reserves estimation of the banned Cu-polymetallic deposit is proposed,that is the Co-Kriging method with calibration weights. Firstly,the original drilling data is processed by statistical analysis,specific values and composite samples,the weight coefficients with banned effect are analyzed in detail; then,the Quasi-Newton method is adopted to calibrate the weight values obtained by the Co-Kriging method. Taking the drilling grade data of 2012 Tibet Jiama copper polymetallic exploration engineering as an example,the reserves of the copper polymetallic deposit is estimated,in particular,Cu grade is regarded as the main regionalized variable and Ag grade is regarded as the synergy regionalized variable to calculate the experimental variation function and cross experimental variation function of Cu and Ag respectively. Besides that,the interpolation results of the Co-Kriging method and Co-Kriging method with calibration weights are analyzed. The research results show that the performance of the Co-Kriging method with calibration weights is superior to the classical Co-Kriging method,and it has some reference for geological attribute analysis and reserve estimation of Cu-polymetallic deposit.
引文
[1]黎枫佶.西藏墨竹工卡县甲玛铜多金属矿三维模型构建和资源量估算[D].成都:成都理工大学,2010.Li Fengji.Constructing 3D Models and Estimating Reserves of Jiama Copper Polymetallic Deposit,Mozhugongka County,Tibet[D].Chengdu:Chengdu University of Technology,2010.
[2]余牛奔,齐文涛,王立欢,等.基于3DMine软件的三维地质建模及储量估算——以新疆巴里坤矿区某井田为例[J].金属矿山,2015(3):138-141.Yu Niuben,Qi Wentao,Wang Lihuan,et al.Three-dimensional geological modeling and reserve estimation based on 3DMine software:taking a well field of Balikun mining area,Xinjiang as an example[J].Metal Mine,2015(3):138-141.
[3]程勖.地质统计学软件开发与应用[D].长春:吉林大学,2009.Cheng Xu.Development and Application of Geostatisties Software[D].Changchun:Jilin University,2009.
[4]唐攀,唐菊兴,唐晓倩,等.传统方法和地质统计学在矿产资源储量分类中的对比分析[J].金属矿山,2013(11):107-109.Tang Pan,Tang Juxing,Tang Xiaoqian,et al.The comparative research of traditional method and geostatistics method in mineral resource/reserve classification[J].Metal Mine,2013(11):107-109.
[5]余璨,李峰,张达兵,等.基于DMINE的红龙厂铜矿床地质建模与储量计算[J].金属矿山,2015(2):108-111.Yu Can,Li Feng,Zhang Dabing,et al.Geological modeling and calculation of the reserves of Honglongchang copper deposit based on DIMINE[J].Metal Mine,2015(2):108-111.
[6]李庆谋.多维分形克里格方法[J].地球科学进展,2005,20(2):248-256.Li Qingmou.Multi-fractal Krige interpolation method[J].Advances in Earth Science,2005,20(2):248-256.
[7]肖克炎,张晓华,王全明,等.应用改进的克里格法分离重力区域异常与局部异常[J].物探化探计算技术,1995,17(2):19-25.Xiao Keyan,Zhang Xiaohua,Wang Quanming et al.Separation of regional field and local field by improved Kriging method[J].Computing Techniques for Geophysical and Geochemical Exploration,1995,17(2):19-25.
[8]高美娟,朱庆忠,张淑华.利用贝叶斯-克里金估计技术进行储层参数预测[J].石油地球物理勘探,1999,34(4):390-397.Gao Meijuan,Zhu Qingzhong,Zhang Shuhua.Reservoir parameter prediction using Bayes-Kriging estimation technique[J].Oil Geophyical Prospecting,1999,34(4):390-397.
[9]李君,李少华,毛平,等.Kriging插值中条件数据点个数的选择[J].断块油气田,2010,17(3):277-279.Li Jun,Li Shaohua,Mao Ping et al.Choosing the number of conditional data in Kriging interpolation[J].Fault-Block Oil&Gas Feild2010,17(3):277-279.
[10]汪保,孙秦.改进的Kriging模型的可靠度计算[J].计算机仿真,2011,28(2):113-116.Wang Bao,Sun Qin.Structural reliability computation based on Kriging model[J].Computer Simulation,2011,28(2):113-116.
[11]孙卡,翁正平,李章林,等.排序扫描线算法在克里格估值中的应用[J].金属矿山,2009(7):73-76.Sun Ka,Weng Zhengping,Li Zhanglin,et al.Application of sequential scanning algorithm in Kriging prediction[J].Metal Mine,2009(7):73-76.
[12]牛文杰,朱大培,陈其明.滑动领域克里金插值法的改进[J].计算机辅助设计与图形学学报,2001,13(8):752-756.Niu Wenjie,Zhu Dapei,Chen Qiming.Improvement of moving neighborhood Kriging method[J].Journal of Computer-Aided Design&Compute Graphics,2001,13(8):752-756.
[13]李少华,王勇标,李君,等.克里金估值中条带效应的校正[J].石油地球物理勘探,2012,47(6):979-983.Li Shaohua,Wang Yongbiao,Li Jun,et al.Correcting the string effect of Kriging[J].Oil Geophyical Prospecting,2012,47(6):979-983.
[14]牛文杰,孟宪海,李吉刚.结合软数据的同位置协同克里金估值新方法[J].煤田地质与勘探,2011,39(2):13-17.Niu Wenjie,Meng Xianhai,Li Jigang.A new estimation method of collocated Co-Kring combined with soft data[J].Coal Geology&Exploration,2011,39(2):13-17.
[15]侯景儒,黄竞先.地质统计学的理论与方法[M].北京:地质出版社,1990.Hou Jingru,Huang Jingxian.Theory and Method of Geostatistics[M].Beijing:Geological Publishing House,1990.
[16]孙仁铎.空间变异理论及其应用[M].北京:科学出版社,2005.Sun Renduo.Spatial Variability Theory and Application[M].Beijing:Science Press,2005.
[17]孙红泉,康永尚,杜惠芝.实用地质统计学程序集[M].北京:地质出版社,1997.Sun Hongquan,Kang Yongshang,Du Huizhi.Utility of Geostatistical Procedures Set[M].Beijing:Geological Publishing House,1997.
[18]易君君,汪保,颜倩倩.基于新拟牛顿方程的优化算法设计及应用[J].宁波工程学院学报,2015,27(1):12-18.Yi Junjun,Wang Bao,Yan Qianqian.Design and application of optimization algorithm based on new Quasi-Newton equation[J].Journal of Ningbo University of Technology,2015,27(1):12-18.
[19]王小朋,刘翔峰.一种基于病态问题的修正牛顿法[J].河南科技大学学报:自然科学版,2015,36(1):86-91.Wang Xiaopeng,Liu Xiangfeng.A modified Newton method based on ill-conditioning problem[J].Journal of Henan University of Science and Technology:Natural Science,2015,36(1):86-91.
[20]张华军,赵金,罗慧.基于拟牛顿法的同时扰动随机逼近算法[J].华中科技大学学报:自然科学版,2014,42(9):1-4.Zhang Huajun,Zhao Jin,Luo Hui.Simultaneous perturbation stochastic approximation algorithm based on Quasi-Newton method[J].Journal of Huazhong University of Science&Technology:Natural Science Edition,2014,42(9):1-4.
