基于稳健统计学和EDA 技术的地球化学异常下限确定
2013-08-07NguyenTienThanh刘修国陈春亮DangThaiSon
Nguyen Tien Thanh ,刘修国 ,陈春亮 ,任 泽,Dang Thai Son
(1.中国地质大学 国际教育学院,武汉 430074;2.河内矿业与地质大学,越南 河内 130510;3.中国地质大学 信息工程学院,武汉 430074;4.中国地质大学 研究生院,武汉 430074;5.越南地质与矿产总局,越南 河内 130510)
0 前言
地球化学异常下限的确定是地球化学勘查的一个基本问题。二十世纪五十年代以来,地球化学家基本认为,地球化学场中的分布接近正态分布或者对数正态分布[1、2、4~9]。因此,使用传统统计方法确定地球化学异常下限时,首先要仔细研究地球化学数据,检验其是否服从正态分布或对数正态分布。为什么数据的分布情况对数据处理有这么重要的影响?因为在传统的统计方法中,相关分析、因子分析、判别分析和经典统计检验包括概率的计算,都是数据在服从正态分布假设的前提下进行的[3]。很多与地球化学相关的文献,都表明地球化学数据往往服从对数正态分布[4~7]。但新的研究表明,元素在地质体中的含量分布,并不局限于正态分布或对数正态分布[2、8、9],所以数据转换[10~13](如取对数、余弦、平方根等)和明显特异值(特高值和特低值)的剔除,直到数据服从正态或对数正态分布为止[1、2、9、13],是两个简单的解决方案。先计 算数据的平均值和标准偏差,然后对全部数据进行异常筛选和评价,一般是以平均值与二倍标准偏差之和作为地球化学异常下限。关键问题是数据转换并不能使数据更接近正态分布[14],忽略这一点,可能会导致错误的结果。至于明显异常点的剔除,数据处理的最终目的是确定异常点,剔除异常数据之前该如何确定哪些点是异常数据,哪些点不是?……
