地球化学异常提取的自适应衬值滤波法-物探与化探.PDFVIP

　第３５卷第４期 物　探　与　化　探 Ｖｏｌ．３５，Ｎｏ．４　 　２０１１年８月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ａｕｇ．，２０１１　 地球化学异常提取的自适应衬值滤波法 １ １，２ 金俊杰 ，陈建国 （１．中国地质大学数学地质遥感研究所，湖北 武汉　４３００７４；２．中国地质大学地质过程与矿产资 源国家重点实验室，湖北 武汉　４３００７４） 摘 要：地球化学数据分布，特别是微量元素的频率分布，是正向偏斜的，其异常下限的计算也是基于对数正态分 布。当定性分析元素的分布状态时，大部分样本分布状态的结果是可以接受的，但是当定量计算异常下限时，一定 要结合具体样本的实际分布情况。引入ｇｈ分布函数对样本数据进行拟合，根据所得到的参数确定异常下限，将异 常下限与标准离差分割开来。结合子区中位数衬值滤波法，对个旧及其外围成矿区的化探数据进行处理，同时自 适应的调节滤波窗口大小，提取化探异常，从而对预测靶区进行圈定。 关键词：自适应；衬值滤波；ｇｈ分布；异常下限；个旧 中图分类号：Ｐ６３２　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２０１１）０４－０５２６－０６ 　　随着信息处理技术的发展，尤其是稳健统计学 了优秀的描述性统计和科学数据可视化分析方法， 在地球化学中的应用，化探数据的处理方法逐渐增 史长义等人引进了ＥＤＡ技术，以滤波的方式处理化 多，如趋势面法、归一化法、滑动平均法、小波分析方 探数据，通过ＥＤＡ计算异常下限，再以衬值的形式 法、衬度异常法、分区标准化方法、分区背景校正法、 圈定异常范围。ＥＤＡ技术中的一套统计参数是相 子区中位数衬值滤波法、局部奇异性分析方法 应滤波框中数据集的中位数（均值），上、下极值点 ［１］ 等 。这些方法围绕的核心内容是：怎样划定异常 和上、下四分点。异常下限的公式为 下限，从而圈定异常范围。传统方法确定异常下限 Ｆ ＝Ｑ ＋１．５×Ｓ， （１） ｕ ｕ ｈ 是基于经典统计学，认为元素服从正态或对数正态 其中，Ｆ为异常下限，Ｑ 为上四分点，Ｓ为内散度 ｕ ｕ ｈ 分布，微量元素频率分布偏斜的原因是样品自身是 （类似于标准离差）。这就避免了特异值对数据处 混合体造成［２］。从２０世纪６０年代起，勘查地球化 理的影响，只保留中心的背景场，所以只需根据数据 学数据均按照均值加２或３倍标准差来计算异常下 本身来确定异常，并不需要对数据进行预处理。具 限。然而针对实际情况，人们越来越认识到它的局 体方法是：将全区进行网格化，采用２个方形窗口套 限性，即背景总体与异常总体常常是重合的。关于 合取值，中心窗口为ｎ×ｎ个单位格子（内窗口），外 异常下限的圈定，前人也有很多研究，如概率图解 部环形窗口为ｍ×ｍ个单位格子（外窗口）；以内窗 ［３］ ［４］ 法、小波多尺度法 、含量排列法 、ＰＣＡ、ＥＤＡ、多 口中数据的中位数作为中心格子的异常值（Ｃ ）， Ｗ，ｍ ［５－６］ 重分形法 等，这些方法都伴有将背景曲面化的 外窗口中数据的中位数作为中心格子的背景值，利 过程。ｇｈ分布是