资源描述：
《克里金法在七里坪多金属矿激电测深数据处理中的应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、文章编号:1672—7940(2008)06—0709—06克里金法在七里坪多金属矿激电测深数据处理中的应用徐顺强1,李清林1,余焱2,谢汝一1,陈安文2(1.中国地震局地球物理勘探中心,郑州450002;2.洛阳坤宇矿业有限公司,洛宁471700)摘要:高分辨率电法勘探已经成为勘探地球物理学界公认的发展方向,目前趋于成熟的是电阻率成像法。本文利用克里金(Kriging)插值方法对洛宁七里坪多金属矿激电测深数据进行处理,目的是得到高分辨率的解释结果,事实证明克里金插值在激电测深数据中能够取得较好的
2、效果。本文详细讨论了克里金法的原理、算法,以及在运用区域的模型设计,并对该方法在激电测深数据处理领域的发展做了有益的探索。关键词:克里金法;多金属矿;激电测深;数据处理中图分类号:P631文献标识码:A收稿日期:2008-07-30TheApplicationoftheKrigingMethodtotheIPSoundingDataProcessingoftheQilipingMulti-metalliferousOreXuShunqiang1,LiQinglin1,Yuyi2,XieRuyi1,
3、ChenAnwen2(1.GeophysicalExplorationCenterofChinaEarthquakeAdministration,Zhengzhou450002,China;2.LuoyangKunyuMiningIndustryLimitedCompany,Luoning471700,China)Abstract:Thehighresolutionelectricalprospectingmethodhasbeenrecognizedasthedirec2tionoftheres
4、earcharea.ThispaperprocessesthedataofQilipingmulti-metalliferousoreIPsoundingbasedontheKriginginterpolationmethod.Thepurposeistoobtaintheinter2pretationresultofthehighresolution.TheresultshowsthattheKriginginterpolationmethodisgoodintheIPsoundingdata.
5、Thispaperdiscussestheprinciple,calculationandmodeldesignoftheKrigingmethodindetailsandmakesagoodexplorationofthedevelop2mentofthemethodinthefieldofIPsoundingdataprocessing.Keywords:Kriging;themulti-metalliferousore;IPsounding;dataprocessing机技术的迅速发展,关于
6、激电测深数据处理模型趋于多样化,但是目前激电数据的二维反演主要针对高密度激电数据,而对常规的激电测深涉及较少。激电测深在异常比较明显、测点密度大的情况下,测量数据无需处理就可以直接勾绘拟断面图,根据拟断面图分析异常。而在异常效果不引言1电法勘探是勘探地球物理学中的主要学科之一,是电学领域中的一门应用科学。数十年来电法勘探这一学科发展迅速并趋于完善,随着计算作者简介:徐顺强(1982-),男,助理工程师,云南禄劝人,2006年毕业于中国科学技术大学,从事地震监测预报、矿产资源勘探等工作。E-mail
7、:sq-xu@126.com十分明显、测点少的情况下,对测量数据进行进一步处理就尤为重要。随着矿产开采业的发展,地表矿、浅部矿已基本探明,找矿工作逐渐转向寻找深部矿、隐伏矿,依靠地表地质现象结合测深数据寻找矿化带很明显已经不合时宜[1,2]。洛宁县七里坪金矿普查区,行政区划隶属洛宁县赵村乡,勘查区地理坐标:东经111°33′30″～111°37′30″,北纬34°10′00″～34°14′00″,面积22132km2。剔除划出的部分区域和周边局部区域,实际勘探区面积约15.9km2。设计剖面59条
8、,点距100m,线距200m。由于测区地形切割较大,实际点、线距远远大于设计的直线距离范围,这样导致实测剖面、实测测点之间出现了“盲区”,在需要的点上没有数据,可能一些激电测深能捕捉到的矿化带部分藏于其中,有些往往决定矿化带的走向,会导致激电测深成果存在较大的误差,这就需要以已知点为基点来进行插值规格化处理[7,8]。本文利用克里金插值方法对七里坪多金属矿激电测深数据进行处理,取得了较好的效果。,z(xn)为相应的观测值(图1)。区域化(x1),变量在x0处的值z3(x0)可采用一
