稳健数据分析方法在地球化学异常阈值计算中的应用研究.pdfVIP

第十二届全国数学地质与地学信息学术研讨会 稳健数据分析方法在地球化学异常阈值计算中的应用研究 周曙光u，周可法1，王金林1 (I、中国科学院新疆生态与地理研究所，新疆矿产资源研究中心，新疆，830011； 2、中国科学院大学北京io0049) 摘要：地球化学异常阈值的计算直接关系到地球化学异常制图结果。近而影响地球化学数据在矿产资源 勘查工作中的应用。本文对比了传统和稳健地球化学异常计算方法在化探异常阈值计算中的特点。结果显 示，传统方法在化探数据分析方面受到限制，而稳健方法可以作为更加可靠的地球化学异常阈值计算方法。 关键词：稳健数据分析；探索性数据分析；地球化学；异常阈值 0前言 传统地，为了确定地球化学异常，人们常常使用“平均值+n倍标准差”的方法来确定某 元素的异常下限。由于该方法简单易用，得到了广泛认可。然而，随着矿产资源勘查工作的 深入开展，明显致矿异常必将越来越少，提取地球化学异常也就更难。这对科研工作者提出 了新的要求。在同～数据源基础上，研发新的地球化学异常提取方法，发现有意义的弱小异 常将会成为一个热点方向。 对于化探异常信息提取，很多学者都做了有益的探索。为了解决单一异常下限方法确定 异常时常常会出现“有矿无异常或有异常无矿”的情况，李宝强【l】在西北地区利用滑动平均衬 值滤波法研究了地球化学异常的提取问题，认为滑动平均衬值滤波法是一种有效的数据处理 方法。子区中位数衬值滤波法被认为可以有效地避免地质背景对异常圈定的影响，解决低背 景区和高背景区弱小异常的识别问题，进行了一些有益的探索[2-4]。为了解决传统参数统计 方法依赖数据正态或对数正态分布的情况，探索性数据分析【5】被认为是一种稳健的数据分 析方法，可以直接分析原始地球化学数据【6】，该方法对数据分布两端25％的数据具有抵抗 力【7】，不会受到数据中极大和极小值的严重影响。在分形理论【8】提出以后，分形和多重分 形方法也在地球化学异常分离领域得到了快速发展，很多分形和多重分形技术在地球化学异 常分离领域得到广泛应用【9·13】。 本文主要对比研究探索性数据分析方法与常规方法在确定地球化学异常阈值方面的特 点，更深入了解地球化学数据的特点及如何选择更合适的地球化学异常阈值确定方法。 1地球化学数据分布 本文所用数据包含1444个样本数据，覆盖面积约7000km2，是中国区域地球化学调查 (ChinasNationalGeochemical MappingProject)数据的一部分，研究区位于新疆克拉玛依 129 第十二届全国数学地质与地学信息学术研讨会 附近．本研究以Ag元素为例进行探讨。 从图1可以看出，Ag元素浓度呈正偏态分布，Ag元素的正态Q-Q图显示，它不符合 正态分布。对Ag元素浓度数据进行对数转换(以自然对数为底)之后(图2)，数据趋向正 态分布，但是，通过对转换后数据的正态q-q图发现，数据仍然不符合正态分布。因此， 原始数据也不符合对数正态分布。 ∞一cIE∞∞一]c|c一_Io∞∞lI一匕∞30 Ag浓度(ppb) 图l Ag元素频数分布直方图和正态Q-Q图 ∞一cII-c∞∽一3 c|c一．Io∞∞|I—cm]o 图2对数转换后Ag元素频数分布直方图和正态Q-Q图 2探索性数据分析