西舍尔法检验特异值处理合理性

——系列文章之二

刘先生的地质新浪博客

（转载必须注明出处）

一、前言

在《固体矿产资源量估算规程·第3部分：地质统计学法》中，附录3讲述了“西舍尔法”检验特异值处理合理性。

二、规范附录3

C.1在大多数的地质问题中，矿体的厚度或矿石的品位（地质变量x）不服从正态分布，但其对数值（即y=ln(x)）趋于正态分布，则称为地质变量服从对数正态分布。

C.2若地质变量不服从对数正态分布，但y=ln(x a)服从正态分布（其中a是一个常数）,则称地质变量服从三参数对数正态分布。

C.3服从对数正态分布或三参数对数正态分布的样品一般存在特异值，对此地质统计学有相应的处理方法。特异值处理的合理性可用西舍尔T估值检验（Sichel's T）。

C.4西舍尔T估值计算公式如下：

C.5西舍尔估值检验西舍尔估值的检验分为如下几个步骤：

a) 首先对已经处理特异值的样品值取自然对数；

b) 计算取自然对数后的样品品位平均值和方差；

c) 将对数平均值返回真值得到几何平均值；

d) 根据对数变换后的品位方差和样品个数，计算西舍尔系数；

e) 将几何平均值乘以西舍尔系数得到西舍尔T估值。

f) 当特异值处理后的样品的算术平均值小于并和西舍尔T估值接近时，判断特异值处理结果合理。

三、实例

国外某钼（Mo）矿，断裂控矿，脉状，中型。

参与资源量估算的样品数:1977件，最小值:0.003，最大值:4.820，平均值:0.356，变化系数:1.509。根据中国规范，取均值8倍作为特高品位（即当Mo品位大于2.5%时），特高品位样品共25件。处理方法：剔除样品，不能与计算。特高品位处理后样品1952件，最小值:0.003，最大值:2.456，平均值:0.336，变化系数:1.096。

利用三维软件进行统计（以3DMine为例），方法如下：

（1）打开三维软件，提取样品点

目的是：只统计参加资源量估算的样品，过滤掉其它样品。

（2）样品统计

操作方法：地质统计