2022年化探-异常下限-计算方法大全及详解.pdf

学习必备欢迎下载化探 -异常下限 -计算方法大全及详解谭亲平地球化学研究所目录1. 传统方法，均值加标准差. 12. 直方图解法 . 23. 概率格纸图解法 . 34. 多重分形法 . 65. 85% 累计频率法 . 7小结. 8传统方法，均值加标准差在 excel 中用过函数，求均值，求标准差，先对数据中的极大/极小值进行剔除，大于/小于三倍标准差的剔除掉，直到无剔除点然后用均值加2 倍标准差求异常下限图， D 列中的函数， E列中的结果图一中的化探数据的异常下限114.86.直方图解法精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 1 页，共 10 页 - - - - - - - - -学习必备欢迎下载精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 2 页，共 10 页 - - - - - - - - -学习必备欢迎下载图 2首先，做频率直方图，（图 1 的数据是某化探区数据）含量频率分布图上呈现双峰曲线，左边是背景部分，右边是异常部分，双峰间谷底处（0.7）为异常下限求真值得5.所以，异常下限位 5。

图 2 另一个化探区的数据，是单峰曲线，在频率极大值的0.6 倍处画一条平行直线，与曲线一侧相交，其横坐标长度即为用 Ca=Co+2*=0.16+2*0.665=1.49 ，求得为真值为31概率格纸图解法 .精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 3 页，共 10 页 - - - - - - - - -学习必备欢迎下载图 3，图 3 是概率格纸发现纵坐标（累计频率）是不均匀的把样本值小于或等于某个样本ni的数据频率累加，即得到小于或等于ni的累积频率概率格纸用excel能轻松的做出来制造方法如下精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 4 页，共 10 页 - - - - - - - - -学习必备欢迎下载图 4.图 4 显示了概率格纸的制造过程原理就是把标准正态分布曲线投影到纵坐标上首先确定纵坐标数值，如B 列， 0.1、1、5、10、20、30、40、50、60、70、80、90、95、99、 99.9. 如果想要纵坐标线密一点，也可以插入更多的数然后在 C 列中用 NORMSINV函数，求对应频率的分位数（如果把标准正态分布，正着放，分位数就是横坐标）。

这时的原点 （0） 在 50%处，我们想要原点在0 处，那么把 C列的数统一加-03.090232（C5） ，-（处理化探数据的时候，加的也是相同的数）即输入公式” D5” =C5-$C$5 E列为 x值，根据实际化探数据，设定最大和最小值我们这里随便设为0、25然后画，“带直线的散点图”，一条线作为一组数，共15 组精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 5 页，共 10 页 - - - - - - - - -学习必备欢迎下载添加 15 次然后设置每一条线的“数据系列格式”，数据标记选项为“无”，线条为“黑色”，线粗为“1 磅”然后把纵坐标删掉，在图表中人为添加每条线纵坐标标签：精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 6 页，共 10 页 - - - - - - - - -学习必备欢迎下载加上一条线的标签后，通过复制， 粘贴，可快速完成，改好数值，调整位置即为图3.现在就是应用了将图 1 中的化探数据，应用如下横坐标的间隔为0.1.求累计频率将图 1 中的化探数据，用概率格纸求解如下。

图 5发现是有两条斜率不等的直线所综合形成的曲线与图1 中的双峰曲线一致左边的直线反映背景，右边的曲线反映异常应用多重母体分解法，以拐点为界，左侧背景占60%，右侧异常占40%将换算成单一背景母体的累计频率=背景部分每个点的累计概率*100/60 异常母体累计频率=（异常部分每个点的累计频率-60）*100/40. 再分别绘图，精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 7 页，共 10 页 - - - - - - - - -学习必备欢迎下载如上图因此背景部分累计频率97.7 处的横坐标即为异常下限（0.76）真值为5.6发现和直方图解法求得的相近多重分形法多重分形法将背景与矿化异常的形成认为是两个相互独立的过程，它们分别满足不同的幂指数分布目前利用分形技术进行地球化学异常下限确定的方法主要有（含量)周长法、（含量 )面积法、 （含量 )距离法、 （含量 )频数法等 , （含量 )求和法，这里采用（含量)求和法进行讲解设分形求和模型:N(Ci)=kCi-D(i0)， 式中 Ci为元素含量， 又称特征尺度， k为比例常数（k0） ，D 为一般分维数， N(Ci)为当元素含量为Ci时所有大于等于Ci的元素含量的和数。

分形求和模型两边分别取对数得到一元线性回归模型:logN(Ci)=-Dlog(Ci)+log(Ci)，用最小二乘法求出斜率D 的估计量，即为分维数，其散点大致分布在两段直线上，采用分段拟合分别求出两段线性方程，两段直线的交点为背景与异常的分界点，即异常下限值精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 8 页，共 10 页 - - - - - - - - -学习必备欢迎下载图 6图 1 中的化探数据在图6 中，异常下限为1.6（40）85%累计频率法把样本值小于或等于某个样本ni的数据频率累加，即得到小于或等于ni的累积频率 一般使用累积频率85%时的样本值作为异常下限累计频率求解方法如下精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 9 页，共 10 页 - - - - - - - - -学习必备欢迎下载，假如有一组数据，0.1、0.2、0.3、0.4将他们排序，1、2、3、4、，将各自的序号除总数，如0.05=1/20 选中 0.85 时的样品值为异常下限用此方法求得图1 中的化探数据异常下限为33.5.。

小结表图 1 中化探数据异常下限传统法（均值加标准差）114.86直方图解法5概率格纸法5.6多重分形法4085%累计频率法33.5发现，不同的方法求的结果相差很大其中，绘制的异常图中，多重分形法，和85%累计频率法比较符合实际工作中按实际情况选择方法精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 10 页，共 10 页 - - - - - - - - -。