化探异常下限计算

《化探异常下限计算》由会员分享，可在线阅读，更多相关《化探异常下限计算（11页珍藏版）》请在装配图网上搜索。

1、化探-异常下限-计算方法大全及详解化探异常下限计算方法大全及详解谭亲平地球化学研究所目录1. 传统方法，均值加标准差12. 直方图解法23. 概率格纸图解法34. 多重分形法。65. 85%累计频率法。7小结8传统方法，均值加标准差在excel中用过函数，求均值，求标准差，先对数据中的极大/极小值进行剔除，大于/小于三倍标准差的剔除掉，直到无剔除点。然后用均值加2倍标准差求异常下限。KIABDE1-0. 61979均值AVERAGE(AA6923432-0. 5S503标准差STDEV (A:A)0.6039033-0. 56S64异常下限R倍标淮差D1+2*D2.Z 060154-0. 55

2、204真值FOWER(10,D3j_114. 8555-0. 53766-0. 5220S7-0522SS1图，D列中的函数，E列中的结果。图一中的化探数据的异常下限114.86.。直方图解法图KHistogram= 1EQ0 54528Std. &ev.= ,337SeiD2OM = 1.023双峰间谷底处（0.7）为异常下限。求真值得5.所以，异常下限位5。用 Ca=Co+2*a =0.16+2*0.665=1.49，求得为真值为 31。图2首先，做频率直方图，（图1的数据是某化探区数据）含量频率分布图上呈现双峰曲线，左边是背景部分，右边是异常部分，图2另一个化探区的数据，是单峰曲线，在频

3、率极大值的0.6倍处画一条平行直线，与曲线一侧相交，其横坐标长度即为a。概率格纸图解法.图3,图3是概率格纸。发现纵坐标（累计频率）是不均匀的。把样本值小于或等于某个样本n,的数据频率累加，即得到小于或等于n,的累积频率。概率格纸用excel能轻松的做出来。制造方法如下。99一 乎 3. 09W2323 6. 1日。如如0L G44 监霸 4. 7350859L油83S5 如蹄2. 32&S4T9 5. 41&5S022. 32&甜?5. 4165502氐妁的色M白住诵狷If .tSTTO 觐率r心 分费 w七升茜数 w 5花HE5 3, 6】战花字 第0. 5244Q05 3. G14E32

4、80D. 212 3. 93105S5 2&0. 5di212 3. 53105350T调疝网21整翁2i的2&穿涉29S31巽责3i5 1叵亟西11 -3.0*02321 -2.S3634S Or7E30S4+1 -2.3?6348 O.7630S44 5 -1.64 4954 j. 4 453187 255 -1.644954 I. 445ST8T10 -1.2*1552 1.808&S07 10 -1.231552 1.8036307 却 -0.BH5S1 2t 2486111 20 -fl.G4ie21 2. 2436111 30 -0.52:440 2,565831S 2530 -0

5、.52440140 -0,25334750 -1.39E-160 0. 25S34T1 3. 343579460 0. 253-3471 3. 34357947Q 70图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.09023

6、2（C5），（处理化探数据的时候，加的也是相同的数）。即输入公式D5” =C5 $C$5.。E列为x值，根据实际化探数据，设定最大和最小值。我们这里随便设为0、25。然后画，“带直线的散点图”，一条线作为一组数，共15组。选择数据源图表数宠区域图例项凉列）数括区域因太复杂而无法显示。如果选择新的数据区域，则将替换系列宙格中 的所有系列；编辑敷据系列系列名称如:Y轴系列值翌匕 =1确定? IX取消添加15次。即可。然后设置每一条线的“数据系列格式”，数据标记选项为“无”，线条为“黑色”，线粗为“ 1磅”。然后把纵坐标删掉，在图表中人为添加每条线纵坐标标签：觐 -!阍咬本框凹13垂亘文本框&。加上

7、一条线的标签后，通过复制，粘贴，可快速完成，改好数值，调整位置即为图3.。现在就是应用了将图1中的化探数据，应用如下。横坐标的间隔为0.1.。求累计频率。将图1中的化探数据，用概率格纸求解如下。图5发现是有两条斜率不等的直线所综合形成的曲线。与图1中的双峰曲线一致。左边的直线反映背景，右边的曲线反映异常。应用多重母体分解法，以拐点为界，左侧背景占60%,右侧异常 占40%。将换算成单一背景母体的累计频率=背景部分每个点的累计概率*100/60。异常母体累计频率=（异常部分每个点的累计频率-60）*100/40.。再分别绘图，如上图。因此背景部分累计 频率97.7处的横坐标即为异常下限（0.76

8、）真值为5.6。发现和直方图解法求得的相近。多重分形法。多重分形法将背景与矿化异常的形成认为是两个相互独立的过程，它们分别满足不同的幕指数分布。目前利用分形技术进行地球化学异常下限确定的方法主要有（含量）周长法、（含量）面 积法、（含量）距离法、（含量）频数法等，（含量）求和法，这里采用（含量）求和法进行讲解。设分形求和模型:NQkq-DC i0），式中C为元素含量，又称特征尺度，k为比例常数（k0），D为一般分维数，N（C）为当元素含量为2时所有大于等于2的元素含量的和数。分形求和模型两 边分别取对数得到一元线性回归模型:logN（C.）=-Dlog（C.）+log（C.），用最小二乘法求出

9、斜率D的估计量，即为分维数，其散点大致分布在两段直线上，采用分段拟合分别求出两段线性方程，两段直线的交点为背景与异常的分界点，即异常下限值。图6图1中的化探数据在图6中，异常下限为1.6 (40)。85%累计频率法。把样本值小于或等于某个样本n,的数据频率累加，即得到小于或等于n,的累积频率。一般使用累积频率85%时的样本值作为异常下限。累计频率求解方法如下氏11。050. 220. 10-33：0. 15G. 44房0.550.-250. 660. 3.Q 770. 35c. e8Qt 4。9g0. 451100. 51.1110. 551.21.2耻6L 3130. 651. 4140, 71. 5150. 751.6160. 81. 7170. 851. 818顷91. 9190. 95201，假如有一组数据，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%累计频率法比较符合实际。工作中按实际情况选择方法