区域成矿预测中化探元素异常下限优化探讨.pdfVIP

第九届全国矿床会议论文集 区域成矿预测中化探元素异常下限优化探讨 白万成 (武警黄金指挥部，北京100055) 2007年，“全国区域地球化学勘查数据库”向公众开放，从根本上解决了区域化探数据共享和大区域 数据准备问题，为在成矿省、成矿带范围内开展矿床统计预测提供了可能。 众所周知，20世纪80年代开始的全国地球化学扫面计划历经20多年才得以完成，期闻不同单位、 不同时期在不同地区以标准图幅为单元开展工作。当时确定元素异常的原则，是以工作图幅(或几个联 测图幅)为数据集合经统计确定异常下限，圈定元素异常。但是，在成矿省、成矿带范围内使用这些数 据时，有两个问题需要先行解决：其一，不同单位、不同时期测量的系统误差，即“数据台阶”问题； 其二，不同地质单元元素含量背景值的差异造成的异常失真问题。恰当地使用归一化方法(刘大文， 2004)能较好地解决上述问题。归一化方法适用于调整由于不同地区背景值的差异造成的局部信息强化 或弱化，其实质是将不同子区的背景值调整到全区数据平均值水平，而保持各子区内数据的相对变化水平 不变。以下讨论所涉及的区域化探数据均指经过上述归一化处理后的数据，称为元素的“视含量”。 勘查地球化学中，“异常”是一个相对概念，它服务于应用目的。尽管对于异常下限的确定原则和方 法，学术界进行了长期研究，有的已经被行业标准推荐使用，但对于小比例尺、大区域的化探数据，异常 下限的确定仍有进一步探讨的必要。笔者认为，作为特定矿种(如金矿)的找矿指示标志，只有那些具 有显著指示意义的元素和分布范围才是真正的“异常”，因此，这里就涉及一个“异常下限优化”问题。 即对一个特定矿种，在一种能够指示其存在的元素中，优选出一个异常下限，使其圈定的范围能更好地指 示该矿种矿床的空间位置。 “信息量法”是矿床统计预测中的常用方法之一。该方法的各个参数均有明确的地质意义，易于为地 质专家理解和接受。这里探索将其用于研究、优化元素异常下限。 N，N ’ oi·o 种矿床存在的信息量。s为研究区单元总数，Ⅳ为含矿单元数，5i为该标志状态出现的单元数，Ⅳi为该标 志状态与矿床同时出现的单元数。以化探元素异常为例，，^M越大，表明该元素异常对B矿种找矿的指示 意义越好，反之如果该数值太小，说明这种元素异常对B矿种找矿没有指示意义。 研究信息量计算公式可以发现，对一个区域而言，含矿单元数Ⅳ和单元总数5是已知的、不变的， 对于某种元素异常歹，其出现的单元数si随着异常下限的不同而变化，下限提高则面积缩小，单元数Si 减小，反之亦然。同时由于异常面积变化，异常单元内的含矿单元数川一般情况下也相应发生变化(但 通常不是成比例变化)，由此造成信息量值的变化。这种异常下限变化导致．si、Ⅳi变化进而影响信息量的 变化，为我们利用信息量法优化元素异常下限奠定了基础。 异常下限优化的基本途径是：以元素视含量均值与标准差之和为基准，在其上下一定范围(不超过2 倍标准差区间)内取不少于5组数值作为待定异常下限进行分布统计，计算其找矿信息量并进行对比。 统筹考虑信息量(厶『)、异常单元数(．sf)和异常单元中的含矿单元数(川)，选定“优化的”异常下限。 一般情况下应考虑以下几点： (1)要求Li大于0．2。如果用各待定异常下限计算的厶，均小于0．2，则认为该元素对金矿找矿不具有 指示意义，不用其作为指示标志。 (2)要求si不大于有效单元数的20％。s，太大，即异常范围过大，异常也就失去了指示意义。 ·582· 白万成：区域成矿预测中化探元素异常下限优化探讨 (3)要求M不小于含矿单元总数的20％。从信息量计算公式可知，如果异常区中含矿单元很少而％ 值较大，则一定意味着5彳艮小，si太小，就失去了在区域上扩大找矿空间的作用。 以上优化异常下限的三项原则可以简称为“三个0．2原则”。 现以全国金矿资源潜力评价工作中化探数据的处理为例，说明在成矿省范围内优化异常下限的方法。 区域化探数据来自“全国区域地球化学勘查数据库”，使用多元地学数据管理与分析系统(GeoExpl)，提