1、河南嵩縣地區河南嵩縣地區金礦資源潛力估算金礦資源潛力估算龔慶杰中國地質大學（北京）地球化學教研室Tel:-mail: GeoMDIS & GeoExpl 培訓河南嵩縣地區金礦資源潛力估算河南嵩縣地區金礦資源潛力估算一、地質概況二、化探數據概況三、化探數據處理四、制作元素地球化學圖五、制作元素地球化學異常圖六、元素資源量估算七、元素地球化學組合異常圖八、剝蝕系數計算與資源量估算一、地質概況一、地質概況面積：面積：110*70km2地層：地層：Ar 太華群太華群地層：地層：Pt1 寬坪群寬坪群地層：地層：Pt21 熊耳群熊耳群地層：地層：Pt22 汝陽群汝陽群地

2、層：地層：Pt22 管道口群管道口群地層：地層：Pt3 欒川群欒川群地層：寒武系地層：寒武系地層：地層：K E N Q巖體：五丈山花崗巖體：五丈山花崗巖巖巖體：花山花崗巖巖體：花山花崗巖巖體：合峪花崗巖巖體：合峪花崗巖巖體：小型斑巖體巖體：小型斑巖體構造：近構造：近EW向斷裂向斷裂構造：近構造：近NW向斷裂向斷裂構造：近構造：近NE向斷裂向斷裂構造：近構造：近SN向斷裂向斷裂礦床：幾十處金礦礦床：幾十處金礦床床礦床：幾處鉬礦床礦床：幾處鉬礦床礦床：幾處鉛鋅礦床礦床：幾處鉛鋅礦床二、化探數據概況二、化探數據概況l1:20化探數據本實例：原始數據數據點： 1774點數據項：12項元素項：Ag As

3、 Sb Hg Cu Pb Zn Au W Sn Mo Bi三、化探數據處理三、化探數據處理l離散數據點位圖l離散數據分布檢驗l數據校正處理（濾波）離散數據點位圖離散數據點位圖取樣空間范圍，存在空白區取樣空間范圍，存在空白區規則網格取樣，取樣間距規則網格取樣，取樣間距2km分分布布檢檢驗驗非對數正態分布非對數正態分布存在奇異數據點（單點異常）存在奇異數據點（單點異常）表現在兩端（極低值與極高值）表現在兩端（極低值與極高值）分分布布檢檢驗驗AuAgAs分分布布檢檢驗驗SbHgCu分分布布檢檢驗驗PbZnW分分布布檢檢驗驗SnMoBi結論：存在奇異數據點（單點異常）結論：存在奇異數據點（單點異常）辦

4、法：濾波辦法：濾波濾波濾波刪除低端數據刪除低端數據l直接刪除低端奇異數據（置為空值）l刪除原則： 低于或等于檢出限 孤立數據點（從小到大排序時明顯偏小）按此原則本實例中刪除數據個數為：Au-1；Ag-6；Sb-2；Cu-1；As-0；Hg-0；W -1；Sn-5；Mo-1；Bi -6；Pb-0；Zn-0； l獲得無低端奇異點的濾波數據本實例：濾波刪低濾波濾波中位數中位數l中位數濾波范圍：5點即可（1.1倍網格間距）采用中位數插值方法獲得網格數據l將網格數據轉換為表數據l獲得無奇異點的濾波數據本實例：濾波中位三種數據的比較三種數據的比較原始數據刪低中位minmaxmax min min minm

5、axAg10 2680010424429 2400As2.11335152.152.518.8Sb0.122.369 0.174 0.203.18Hg109400110010600三種數據的比較三種數據的比較原始數據刪低中位minmaxmax min min minmaxCu5525111181010101Pb8.6 13279358.63 13.5 1242Zn 39.1610078 39.17 48.3615Au 0.025541110.300.342.7三種數據的比較三種數據的比較原始數據刪低中位minmax max min min minmaxW0.3236870.47 0.65 37

6、.6Sn0.736251.2-20.911Mo 0.09 61.2811 0.3370.4 25.4Bi0.02 28.117 0.06-1 0.05 5.55結論：中位數濾波數據穩健結論：中位數濾波數據穩健 利用中位數濾波數據利用中位數濾波數據濾濾波波分分布布檢檢驗驗非對數正態分布非對數正態分布無低端奇異數據點無低端奇異數據點AuAgAs濾濾波波分分布布檢檢驗驗非對數正態分布非對數正態分布存在高端奇異點存在高端奇異點SbHgCu濾濾波波分分布布檢檢驗驗PbZnW濾濾波波分分布布檢檢驗驗SnMoBi濾濾波波分分布布檢檢驗驗結論：仍存在奇異數據結論：仍存在奇異數據僅表現在高端（固有性質）僅表現在

7、高端（固有性質）四、制作元素地球化學圖四、制作元素地球化學圖l數據：中位數濾波網格數據l方法：累頻，19級l特點：與數據絕對大小無關l優點：每一元素均存在異常； 背景與異常對比清晰l缺點：出現頻率高的相同數據可能劃歸不同級元素地球化學圖元素地球化學圖 Au Ag As Sb元素地球化學圖元素地球化學圖 Hg Cu Pb Zn元素地球化學圖元素地球化學圖 W Sn Mo Bi五、制作元素地球化學異常圖五、制作元素地球化學異常圖l異常與背景的區分：定值（異常下限）；變值（趨勢面）l本實例采用定值（異常下限）的原因：地球化學景觀相同；數據范圍不大l確定異常下限的方法：累頻方式；均值標準差方式；分形方

8、式l本實例選用累頻方式（其它方式專門討論）累頻方式確定異常下限及濃度分帶累頻方式確定異常下限及濃度分帶l異常下限：85%l背景：0%85，顏色索引9l異常三級分帶：外帶：85%92% ，顏色索引272中帶：92%98% ，顏色索引274內帶：98%100% ，顏色索引278計算異常下限計算異常下限異常下限異常下限元素 異常下限 中帶下限 內帶下限Ag140180350As10.611.614.1Au2.94.812.5Bi0.550.771.9Cu313649Hg405070Mo1.92.88.1Pb65.994149.7Sb0.931.011.4Sn3.94.87.4W 4.986.712.

9、4Zn131148197濃度分帶濃度分帶元素地化異常圖元素地化異常圖 Au Ag As Sb元素地化異常圖元素地化異常圖 Hg Cu Pb Zn元素地化異常圖元素地化異常圖 W Sn Mo Bi六、元素資源量估算六、元素資源量估算1243567Au異常圖異常圖1243567異常序號異常名稱異常面積數據個數內帶均值豐度估計面金屬量1祈雨溝13.153 20.82501002窯溝36.839 17.26191693前河、店房39.1410 23.79484484牛頭溝、萑香洼 31.229 18.86772275上宮、虎溝11.06316192426康山7.722 25.52041047松里溝3

10、16.419747Au異常異常6個內帶個內帶豐度估計法資源量預測豐度估計法資源量預測l類比法思想l計算公式研究區元素質量=異常面積*深度*密度*異常均值當深度與密度一致時僅取決于異常面積和異常均值預測區控制區元素質量預測資源量成礦系數元素質量已知總儲量預測區控制區均值異常面積預測資源量成礦系數均值異常面積已知總儲量*面金屬量法資源量預測面金屬量法資源量預測l類比法思想l計算公式成礦元素質量=異常面積*單位面積面金屬量*深度*密度面金屬量=異常面積*異常超量=異常面積*(均值-下限)當深度與密度一致時僅取決于異常面積和異常超量預測區控制區超量異常面積預測資源量成礦系數超量異常面積已知總儲量*異常

11、序號異常名稱豐度估計面金屬量豐度估計資源量面金屬量資源量1祈雨溝25010030* 30* 2窯溝61916974 51 3前河、店房948448114 134 4牛頭溝、萑香洼67722781 68 5上宮、虎溝1924223 13 6康山20410424 31 7松里溝1974724 14 Au異常異常6個內帶的資源量預測個內帶的資源量預測剝蝕系數剝蝕系數l隨剝蝕深度不同，異常規模也不同。因此不需引入剝蝕系數，如豐度估計法、面金屬量法。l但異常在垂向上若具有對稱性，即礦體異常規模強，頭暈和尾暈弱，難以區分頭暈和尾暈。l直接導致類比時垂向深度相同，異常均值向下是增強或減弱，即垂向不均勻。l因

12、此引入：剝蝕系數七、元素地球化學組合異常圖七、元素地球化學組合異常圖l選定成礦指示元素成礦元素：Au待選指示元素：Ag As Sb Hg Cu Pb Zn W Sn Mo Bil元素地球化學組合異常圖異常空間套合Au Ag As Hg 組合異常圖組合異常圖Au Sb Cu 組合異常圖組合異常圖Au Pb Zn 組合異常圖組合異常圖Au W Sn 組合異常圖組合異常圖Au Mo Bi 組合異常圖組合異常圖八、剝蝕系數計算與資源量估算八、剝蝕系數計算與資源量估算l頭暈：As Sb Hgl礦體：Au Ag Cu Pb Znl尾暈：W Sn Mo Bil數據處理：均勻化（除以均值），消除量綱；頭暈、礦

13、體、尾暈各自累加；剝蝕系數=(體/體頭+尾/尾體+2尾/尾頭)/4剝蝕系數計算示例剝蝕系數計算示例礦體異常值暈比值剝蝕分帶 尾暈 礦體 頭暈 尾/尾頭 體/體頭 尾/尾體 系數無礦01200.33 0.00 0.08 礦頭1240.20.33 0.33 0.27 礦體2420.50.67 0.33 0.50 礦尾4210.80.67 0.67 0.73 無礦21011.00 0.67 0.92 GeoExpl軟件實現軟件實現l數據均勻化l多元素疊加分析l剝蝕系數分布檢驗l剝蝕系數網格化（最近點、1.1倍網格間距）l等值區與等值線l剝蝕系數與金異常空間對比分析金礦床與剝蝕系數金礦床與剝蝕系數考慮

14、剝蝕系數的資源量預測考慮剝蝕系數的資源量預測預測區控制區剝蝕系數元素質量預測資源量剝蝕系數元素質量已知總儲量/控制區控制區剝蝕系數元素質量已知總儲量控制區剝蝕系數元素質量已知總儲量成礦保存系數/異常序號異常名稱豐度估計面金屬量剝蝕系數豐度資源量面金屬量資源量1祈雨溝2501000.60230* 30* 2窯溝6191690.75559 40 3前河等9484480.582118 139 4牛頭溝等6772270.57286 72 5上宮等192420.45530 17 6康山2041040.38039 49 7松里溝197470.54426 16 Au異常異常6個內帶的資源量預測個內帶的資源量預測預測 資源量豐度豐度 估計估計面金 屬量 已探明序號 異常名