文档介绍:地球化学异常下限确定方法
在已知区做地化剖面:要求剖面较长,穿过矿化区(含蚀变区)和正常地层(背景),能区分含矿区和非矿区就可确定为下限。
2、概率格纸法:(可以不考虑野值)
以含量和频率作图
15%--负异常
50%--背景值
85%--X+δ(高背景)
98%-- ( X+2δ)异常下限
3、直方图法:(可以不考虑野值)
能分解出后期叠加的值就为异常下限
4、马氏距离法:(在计算时已考虑野值)
针对样本,实际为建立在多元素正态分布基础之上—多重样本的正态分布,超出椭球体时—异常样(如P3点)。
相似于因子得分的计算,最后为一个剔除异常样本时的计算值,实际计算出综合异常边界线。
当令m=1时,
上式化解为Xa=Xo±KS,这是我们较为熟悉的单元素(一维)计算异常下限常用公式。
该方法计算较为复杂:下面给出一个实例:
马氏距离(黑色虚线)圈定异常基本为两种以上元素异常的重合的部分。
上图中Hy-44与Hy-45综合异常中,由于As元素相连,传统方法无法分割。用该方法可分解为两个异常,后来实际查证中也证明:左边Hy-44为Au、Cu、Co的成矿,右边Hy-45为Au、Cu的成矿。
解决手工的随意性。
5、单元素计算法:(必须剔除野值)
Xa=Xo±KS
(Xa—异常下限,Xo—背景值,K—取值系数,S—标准离差)
从标准正态累积频率密度函数公式:
推断出当K=、2、3时,密度函数分别为95%,%,%
一般为计算方便,通常取K=2,这就是Xa=Xo±2S的来源。
通常应用时,用Xa=Xo±3S无限循环剔除,直到无剔除数据时,对于地球化学通常几百—上千的数据,基本保证数据为正态分布。则此时Xa=Xo+2S定为异常下限
为保证数据为正态分布,实际计算时先将数据转换为对数,此时由于数据离差变小,在剔除野值后,基本都能保证为正态分布。为进行下步计算处理有了理论保障。
6、数据排序法:(不考虑野值)
比较简单、实用
所有数据从小到大,按含量排序
做图(含量—纵坐标,1,2……n含量顺序序列----横坐标)
异常有明显的一斜率
但数据太多时不适合
7、累积频率法:(不考虑野值,在使用时为网格化数据)
目前较为普遍
元素含量高低分级,采用累频分级方式,分19级,
分级频率:
--2-3--8-15-25-40-60-75-85-92
--97-98---100(%)
异常85-90-95-100(%)和<15%
8、实际使用异常下限值的确定:
实际上各方法确定的异常下限都是可行的,关健是确定的这个值合不合理是值得商榷的。
在1:20万区域化探中,由于一般取水系沉积物,样品经过了充分的均一化,方差较小,数据基本为正态分布,剔除不了几个野值,此时计算下限与实际使用值变化不是很大(当然1:20万或1:25万由于区域较大,各分区中元素背景不一,异常下限是不同的,应该适当考