王京貴 ：我國金資源總量的估算方法
- 作者：
- 安徽新天源建設咨詢有限公司
- 最后修訂：
- 2020-12-21 15:28:33
摘要：
一、概率法
1、中國科學技術大學黎彤先生(sheng)和冶金工(gong)業部地(di)質(zhi)局(ju)杜(du)春林先生(sheng)在《我(wo)國金礦(kuang)(kuang)資(zi)源潛力(li)的(de)估計和找礦(kuang)(kuang)問題的(de)建議(yi)》(1982年，內部資(zi)料) 中提出我(wo)國金礦(kuang)(kuang)的(de)總(zong)儲量(∑R) 為(wei)：
∑R = D'x F
其中 D' 為我國單位深度(1Km)陸殼內的金分布量，F 為我國金礦床的成礦率.根據我國陸地面積 9.64286x10 6 Km2，陸殼平均厚度 45 Km，陸殼內巖石平均密度 2.8 克/ cm3 以及中國陸殼(ke)內金的豐度 3.5 毫克 / 噸， 計算出我(wo)國陸殼(ke)內的金元素(su)總量(liang)(D).
D = 9.64286 x 106 公里2 x 45 Km x 2.8 克/厘米3 x 3.5 x 10-9 = 4.2 x 10-9 噸
再根據(ju)金(jin)礦山的開采深度 1 Km 范圍內，可(ke)計(ji)算我國在(zai) 1 Km 深度的陸(lu)殼中的金(jin)元(yuan)素量為:
D 4.2525 x 109 噸
D' = --- = ----------------- = 94.5 x 106 噸/公里(li)
α 45 Km
限(xian)于成礦率(F) 未(wei)知(zhi)， 將金的(de)總儲量(∑R) 分解成三項之和
∑R = R1 + R2 +R3
其中 R1 為我國歷代黃金產量所消耗的金礦儲量； R2 為我國現有的金礦儲量(金屬量)；R3 為我國今后可能獲(huo)得的金礦儲量(它是個未(wei)知數(shu)).
∑R = R1 + R2 + R3 = 3600 + 1700 + 12985 = 18285 噸(dun)
成礦率 F = ∑R / D' = 18285 / 94.5 x 106 = 1.934920635 x 10-4
2、東華(hua)理工大學王京貴(1984年，未曾發表)提(ti)出用概率(lv)法(fa)計算(suan)我國金礦資源總(zong)量的方法(fa)。
∑R = 中國陸地面積 x 計算(suan)深度 x 巖石平均密度 x 金的豐(feng)度 x 金的成(cheng)礦(kuang)概率
AR H D A P
其中(zhong)(zhong)的金(jin)成礦概率用金(jin)在地殼中(zhong)(zhong)的分配系(xi)數來近(jin)似代替.
金在地殼(ke)中(zhong)的分(fen)配系數(shu) = 地殼(ke)中(zhong)的金 / 地球中(zhong)的金
= 8400 x 108 / 48 x 1014
= 1.75 x 10-4
計算結果 ∑R = 9.6 x 106 x 1 x 2.8 x 3.5 x 10-9 x 0.000175
= 16464 噸
如果根據黎彤等提供的我國金資源數據 ∑R = 18285 噸及我國陸殼金總量(9450 萬噸)計算的成礦概率為 1.934920635 x 10-4.
我國金(jin)資源總量(liang)的計算(suan)結果為:
∑R = 9.6 X 106 X 1 X 2.8 X 3.5 X 10-9 X 1.934920635 X 10-4
= 18203.7 噸
與(yu)黎彤先生等(deng)的計算結果基本一致(zhi)。
3、應(ying)用(yong)"鑲嵌(qian)模型"估(gu)算每個地質(zhi)變(bian)量的礦產出概率(lv)，再應(ying)用(yong)
m5 (m'5)
P = ( Pm x ----- + Pm' x ------ ) / 2
155 155
其中 Pm 為產出(chu)概率(lv)(lv)(lv) < 0.65 的(de)(de) m 個(ge)地(di)質(zhi)變(bian)量(liang)的(de)(de)概率(lv)(lv)(lv)連乘， Pm' 為產出(chu)概率(lv)(lv)(lv)< 0.75 的(de)(de) m' 個(ge)地(di)質(zhi)變(bian)量(liang)的(de)(de)概率(lv)(lv)(lv)連乘。
二、概率分布(bu)法--吉波夫分布(bu)法
根據礦(kuang)床預(yu)測(ce)礦(kuang)田(tian)， 礦(kuang)田(tian)預(yu)測(ce)礦(kuang)帶(dai)， 礦(kuang)帶(dai)預(yu)測(ce)全國的思路進行估算。
1 1 1
∑R = Rmax ( 1 + --- + --- + ...... + ------- )
2 3 Rmax
其中(zhong) Rmax 是全國最大(da)礦帶(礦田(tian)或礦床)的(de)儲量(liang).
取膠東礦田資源量(liang) 700 噸(dun)(dun)， 則膠遼(liao)(liao)成礦帶(dai)的資源量(liang)為(wei)(wei) 700 x 7.12901= 4990.307 噸(dun)(dun). 以膠遼(liao)(liao)礦帶(dai) 5000 噸(dun)(dun)， 預(yu)測全國的資源量(liang)為(wei)(wei) 5000 x 9.09451= 45472.55 噸(dun)(dun)。
三、逆推法
計算公式為
m m
∑R = ∑Ri = t ∑SiMi
i=1 i=1
首先根據第 i 個有利成礦區的金礦資源量 Ri (本例中采用膠東)， 利用面積 Si 和成礦有利度 Mi， 可求(qiu)出分配系數 t 。再應用上式(shi)求(qiu)出金的資(zi)源總量。
對(dui)不同地區逆推的結果進行平均， 作為最終的資源總(zong)量。
四、富集系數法:
首先(xian)求(qiu)出已(yi)知礦(kuang)田(有利成(cheng)礦(kuang)區(qu))中的(de)富集系數。已(yi)知礦(kuang)田(有利成(cheng)礦(kuang)區(qu))越多結果越準確。最后求(qiu)其平均富集系數(F)， 則
; m
∑R = F x ∑Si
i=1
m
Si 為第 i 個有利成礦區的面積，∑Si 為全國有利成礦區的總面積。
i=1
由(you)于不(bu)同(tong)預(yu)測方法(fa)所圈定(ding)的(de)面(mian)積(ji)不(bu)同(tong)， 則可將不(bu)同(tong)方法(fa)計算得到(dao)的(de)資源量進行平(ping)均， 作為最終的(de)資源總量。
五、數量化模型Ⅰ法
根據15個變量(K1，J3，C3，C2，D3，D2，D1，S3，O3，O1，∈3，Z2，Z1，NW，γ5) 進(jin)行數量化模型Ⅰ法計算結果(guo)建立(li)如下回歸方(fang)程：
^
Y = 13.9595 + 24.4289 K1 + 22.1969 J3 - 14.9519 C3 + 11.2120 C2 +
73.7085 D3 - 83.4407 D2 + 20.0677 D1 - 22.0692 S3 - 84.8632 O3
+ 69.7806 O1 + 17.4474 ∈3 + 55.8213 Z2 - 70.5052 Z1 + 15.9616
NW + 24.4562 γ5
預測結果， 全(quan)國(guo)金資源(yuan)量為 37179 噸。
綜合上述計算結(jie)果，我國金資源量為：
概率法(fa)
1、中國科學技術大學黎彤(tong) 18285 噸
2、東華理(li)工(gong)大學王京貴(1984年) 16464 噸(dun)(dun) -- 18203.7 噸(dun)(dun)
3、應用"鑲嵌模型"估(gu)算
概(gai)率分(fen)布法：
吉波(bo)夫分(fen)布法(fa) 4990.307 噸 -- 45472.55 噸
逆推法：
數量化模(mo)型Ⅰ法 ; 37179 噸
------------------------------------------------
平均 19222.83噸 – 29785.06噸
