化探单元素异常统计与计算

刘先生的地质 新浪博客

（原创，转载必须注明出处）

一、前言

本人曾写过一记博客文章《化探数据处理方法与流程（手把手级）》，但不够详尽。本文尽量列出统计内含统计内容、计算公式、以及相应的Excel函数，方便各位人员的使用。

二、化探单元素异常统计内容

1、异常（图中圈定）

2、样品个数N

Excel函数：COUNTIF（条件计数）；

3、异常面积S（由图中圈出）

4、样品最大值Max

Excel函数：Max；

5、样品最小值Min

Excel函数：MIN；

6、算术平均值

公式： 

Excel函数：AVERAGE；

7、几何平均值

 公式：(式2-7)；

Excel函数：GEOMEAN；

8、中位数

公式：把这列数按从大到小(或从小到大）排列出来,若数据个数为奇数个,则最中间的一个数即中位数；若数据个数为偶数个,则最中间的两个数的平均数即中位数.

Excel函数：MEDIAN；

9、标准离差S₀

 公式： (式2-9)；

Excel函数：stdev.p；

10、变异系数Cv

公式：Cv=S₀/ AVERAGE (式2-10)；

11、异常下限T

公式：T= X+(2-3)*S₀   (式2-11)；

12、异常衬度Ac

公式：  (式2-12)；

13、异常规模Ad

 公式：(式2-13)；

14、异常NAP值

公式： (式2-14)；

15、浓集克拉克值(C)

计算公式： (式2-15)；

16、致矿系数(Z)

计算公式：Z＝Cv(全区)＋10×Cv(剔高值后)＋100×高值比例＋C   (式2-16)；

其中：高值是大于3倍的标准离差；C是浓集克拉克值。

17、偏度

公式： (式2-17)；

Excel函数：SKEW；

偏度临界值=2√(6/N)

18、峰度

公式： (式2-18)；

Excel函数：KURT；

峰度临界值=2√(24/N)

19、正态检验

公式： (式2-19)；

其中Xi：组中值或含量值；

fi：Xi所对应的频数；

F(X)：为经验累计频率；

F1(X)：为理论累计频率。

Excel函数：NORMDIST；

20、相关系数：

公式： (式2-20)；

其中d为对比序列的序差；n为对比序列的等级数。

Excel函数：CORREL。

三、综合分析

单元素异常统计的基础上进行综合分析包括：化探背景分析、聚类分析（R型与Q型）、因子分析、回归分析、对应分析、相关分析、特征分析、点聚图0-1法、隶属度和贴近度、信息量预测、原生晕元素分带序列、剥蚀度评价等。