1、化探-異常下限-計算方法大全及詳解譚親平地球化學研究所目錄1.傳統方法，均值加標準差. 12.直方圖解法. 23.概率格紙圖解法.34.多重分形法。. 65.85%累計頻率法。. 7小結. 8傳統方法，均值加標準差在excel中用過函數，求均值，求標準差，先對數據中的極大/極小值進行剔除，大于/小于三倍標準差的剔除掉，直到無剔除點。然后用均值加2倍標準差求異常下限。圖，D列中的函數，E列中的結果。圖一中的化探數據的異常下限114.86.。直方圖解法圖2首先，做頻率直方圖，（圖1的數據是某化探區數據）含量頻率分布圖上呈現雙峰曲線，左邊是背景部分，右邊是異常部分，雙峰間谷底處（0.7）為異常下限。

2、求真值得5.所以，異常下限位5。圖2另一個化探區的數據，是單峰曲線，在頻率極大值的0.6倍處畫一條平行直線，與曲線一側相交，其橫坐標長度即為。用Ca=Co+2*=0.16+2*0.665=1.49，求得為真值為31。概率格紙圖解法.圖3，圖3是概率格紙。發現縱坐標（累計頻率）是不均勻的。把樣本值小于或等于某個樣本ni的數據頻率累加，即得到小于或等于ni的累積頻率。概率格紙用excel能輕松的做出來。制造方法如下。圖4.圖4顯示了概率格紙的制造過程。原理就是把標準正態分布曲線投影到縱坐標上。首先確定縱坐標數值，如B列，0.1、1、5、10、20、30、40、50、60、70、80、90、95、9

3、9、99.9.。如果想要縱坐標線密一點，也可以插入更多的數。然后在C列中用NORMSINV函數，求對應頻率的分位數（如果把標準正態分布，正著放，分位數就是橫坐標）。這時的原點（0）在50%處，我們想要原點在0處，那么把C列的數統一加-03.090232（C5），-（處理化探數據的時候，加的也是相同的數）。即輸入公式”D5” =C5-$C$5。E列為x值，根據實際化探數據，設定最大和最小值。我們這里隨便設為0、25。然后畫，“帶直線的散點圖”，一條線作為一組數，共15組。添加15次。即可。然后設置每一條線的“數據系列格式”，數據標記選項為“無”，線條為“黑色”，線粗為“1磅”。然后把縱坐標刪掉，

4、在圖表中人為添加每條線縱坐標標簽：。加上一條線的標簽后，通過復制，粘貼，可快速完成，改好數值，調整位置即為圖3.。現在就是應用了將圖1中的化探數據，應用如下。橫坐標的間隔為0.1.。求累計頻率。將圖1中的化探數據，用概率格紙求解如下。圖5發現是有兩條斜率不等的直線所綜合形成的曲線。與圖1中的雙峰曲線一致。左邊的直線反映背景，右邊的曲線反映異常。應用多重母體分解法，以拐點為界，左側背景占60%，右側異常占40%。將換算成單一背景母體的累計頻率=背景部分每個點的累計概率*100/60。異常母體累計頻率=（異常部分每個點的累計頻率-60）*100/40.。再分別繪圖，如上圖。因此背景部分累計頻率97

5、.7處的橫坐標即為異常下限（0.76）真值為5.6。發現和直方圖解法求得的相近。多重分形法。多重分形法將背景與礦化異常的形成認為是兩個相互獨立的過程，它們分別滿足不同的冪指數分布。目前利用分形技術進行地球化學異常下限確定的方法主要有（含量)周長法、（含量)面積法、（含量)距離法、（含量)頻數法等, （含量)求和法，這里采用（含量)求和法進行講解。設分形求和模型:N(Ci)=kCi-D(i>0)，式中Ci為元素含量，又稱特征尺度，k為比例常數（k>0），D為一般分維數，N(Ci)為當元素含量為Ci時所有大于等于Ci的元素含量的和數。分形求和模型兩邊分別取對數得到一元線性回歸模型:lo

6、gN(Ci)=-Dlog(Ci)+log(Ci)，用最小二乘法求出斜率D的估計量，即為分維數，其散點大致分布在兩段直線上，采用分段擬合分別求出兩段線性方程，兩段直線的交點為背景與異常的分界點，即異常下限值。圖6圖1中的化探數據在圖6中，異常下限為1.6（40）。85%累計頻率法。把樣本值小于或等于某個樣本ni的數據頻率累加，即得到小于或等于ni的累積頻率。一般使用累積頻率85%時的樣本值作為異常下限。累計頻率求解方法如下，假如有一組數據，0.1、0.2、0.3、0.4。將他們排序，1、2、3、4、，將各自的序號除總數，如0.05=1/20選中0.85時的樣品值為異常下限。用此方法求得圖1中的化探數據異常下限為33.5.。小結表圖1中化探數據異常下限傳統法（均值加標準差）114.