王玉祥，常 河，賈福聚，蔡家馭，蘇志宏

(昆明理工大學(xué)國土資源工程學(xué)院，云南 昆明 650093)

0 引 言

礦產(chǎn)資源儲量估算是礦產(chǎn)勘查工作中最重要的部分，也是礦產(chǎn)勘查工作成果的總結(jié)，所以在資源儲量估算過程中，依據(jù)礦區(qū)礦體真實(shí)賦存情況，在保證數(shù)據(jù)精度的情況下，選擇合理的估算方法進(jìn)行資源儲量估算尤為重要[1-3]。過去多年，傳統(tǒng)幾何學(xué)方法在資源儲量估算中廣泛應(yīng)用，且受到自然資源部的認(rèn)可，但隨著礦業(yè)技術(shù)及數(shù)字化時(shí)代的技術(shù)革新，通過地質(zhì)統(tǒng)計(jì)學(xué)方法進(jìn)行更真實(shí)的礦山模擬，從而進(jìn)行資源儲量估算成為了現(xiàn)實(shí)。在新技術(shù)推廣與技術(shù)迭代過程中，可靠性至關(guān)重要，驗(yàn)證其準(zhǔn)確性，是一個(gè)必要階段，所以使用原先受大眾認(rèn)可的傳統(tǒng)幾何學(xué)方法與地質(zhì)統(tǒng)計(jì)學(xué)方法進(jìn)行相互對比分析驗(yàn)證，更能直觀體現(xiàn)新技術(shù)的可靠性，同時(shí)再與實(shí)際開采數(shù)據(jù)相互印證，體現(xiàn)新方法的優(yōu)越性與準(zhǔn)確性。

本文以云南蒙自白牛廠銀多金屬礦區(qū)對門山礦段為例進(jìn)行資源儲量估算，其位于滇東南多金屬成礦區(qū)中部，礦區(qū)礦床Ag、Pb、Zn、Sn總儲量達(dá)到了大型-超大型規(guī)模，具有較高的經(jīng)濟(jì)價(jià)值[4-6]。目前礦山已進(jìn)行數(shù)十年的開采及礦床勘探工作，本文綜合歸納近年來地質(zhì)及探采數(shù)據(jù)，重新核實(shí)礦區(qū)資源儲量，分別以傳統(tǒng)的水平投影地質(zhì)塊段法和基于3DMine軟件的地質(zhì)統(tǒng)計(jì)學(xué)法進(jìn)行資源儲量估算對比分析。

1 地質(zhì)概況

1.1 礦區(qū)地質(zhì)

云南蒙自白牛廠銀多金屬礦區(qū)位于華南褶皺系-滇東南臺褶帶-文山-富寧褶皺束-薄竹山拱褶北西部，白牛廠北東向短軸背斜北東傾伏端，西接揚(yáng)子準(zhǔn)地臺，南連越北古陸。礦區(qū)由咪尾礦段、白羊礦段、穿心洞礦段、對門山礦段及阿尾礦段組成，區(qū)域內(nèi)古生界地層較發(fā)育，由老到新分別為寒武系、下奧陶統(tǒng)、泥盆系、石炭系、二疊系，中生界三疊系地層僅分布于礦區(qū)外圍，缺失中奧陶統(tǒng)、上奧陶統(tǒng)、志留系、侏羅系、白堊系、古近系及新近系[7]。寒武系主要發(fā)育在白牛廠礦區(qū)及其附近，是區(qū)內(nèi)重要的含礦地層，下統(tǒng)以絹云母粉砂質(zhì)板巖為主，頂部為核形石鮞狀灰?guī)r；中統(tǒng)主要為泥巖、板巖夾灰?guī)r和粉砂巖；上統(tǒng)由白云巖夾粉砂巖、泥巖組成，厚度大于2 120 m。區(qū)內(nèi)構(gòu)造較簡單，主體構(gòu)造線方向?yàn)楸睎|向，由一系列褶皺及壓扭性斷裂組成，其次為北西向及南北向斷裂構(gòu)造。其中巖漿活動較強(qiáng)烈，主要表現(xiàn)為印支早期基性巖漿噴發(fā)及燕山期酸性巖漿侵入，形成分布于東部的上二疊統(tǒng)峨眉山玄武巖，東南部的薄竹山花崗巖，白牛廠礦區(qū)的阿尾隱伏花崗巖及一些零星分布的輝綠巖、花崗斑巖、二長巖脈。

礦區(qū)內(nèi)礦產(chǎn)豐富，主要有白牛廠礦區(qū)的銀、錫、鉛、鋅多金屬礦，茅山洞銻礦，黃鐵礦、菱鐵礦、無煙煤等。此外，在區(qū)外薄竹山花崗巖體接觸帶附近盛產(chǎn)鎢、錫、鐵、砷、鉛、鋅、銅、銀等礦產(chǎn)。

1.2 礦段地質(zhì)

本次研究的對門山礦段位于礦區(qū)中部，西起63勘探線，東至210勘探線，南自大尖坡，北至對門山以北F3斷裂地表出露附近，面積約2.795 km2。礦段出露地層以中寒武統(tǒng)、下寒武統(tǒng)為主，其次為下泥盆統(tǒng)，兩者間為低角度不整合接觸，缺失上寒武統(tǒng)、奧陶系和志留系。賦礦層位為中寒武統(tǒng)田蓬組頂部及中上部，巖性為灰-深灰色的一套中厚層白云巖與薄-中厚層粉砂巖、泥質(zhì)粉砂巖互層的巖石。構(gòu)造以北西向褶皺、斷裂為主，北東向及南北向斷裂、褶皺穿插其間組成基本構(gòu)造格架，其中北西向構(gòu)造是礦區(qū)主要的控礦及容礦構(gòu)造。巖漿巖以燕山期酸性巖為主，印支期基性巖不甚發(fā)育，二者皆為侵入巖，區(qū)內(nèi)沒有噴出巖[8-9]。對門山礦段主要發(fā)育晚期蝕變，常見晚期鐵錳白云石和次生泥質(zhì)斷續(xù)分布于礦體上覆地層中，在巖石化學(xué)成分上難以反映出它們的存在和變化，僅硫含量有從礦體向圍巖逐漸降低的趨勢(圖1)。

礦段內(nèi)礦石礦物成分復(fù)雜，礦物種類共計(jì)61種，其中原生礦物53種，次生礦物8種。主要金屬礦物有閃鋅礦、方鉛礦、黃鐵礦、鐵閃鋅礦；主要錫石礦物有錫石和黝錫礦；主要銀礦物有銀黝銅礦、黝銻銀礦。

1.3 礦體特征

對門山礦段共分為3個(gè)成礦帶(Ⅰ礦帶、Ⅴ礦帶和銅礦帶)，11個(gè)礦體(Ⅰ1礦體～Ⅰ4礦體、Ⅴ礦體、Ⅴ-2礦體、Ⅴ-3礦體、Ⅴ-6礦體、銅1礦體～銅3礦體)。由于其構(gòu)造環(huán)境、地層條件、巖體特征等地質(zhì)條件的特殊性，礦體連續(xù)較差，且分散的分布于巖體周圍，各個(gè)礦體相對獨(dú)立，由蝕變帶相互聯(lián)系。礦體厚度沿走向呈厚-薄-厚-薄，沿傾向呈薄-厚-薄的變化規(guī)律，除Sn以外的其他成礦元素品位變化穩(wěn)定程度均屬于穩(wěn)定(圖2)。

1-第四系；2-下泥盆統(tǒng)芭蕉箐組；3-下泥盆統(tǒng)坡腳組；4-下泥盆統(tǒng)坡松沖組；5-中寒武統(tǒng)龍哈組；6-中寒武統(tǒng)田蓬組；7-中寒武統(tǒng)大丫口組；8-下寒武統(tǒng)大寨組；9-下寒武統(tǒng)沖莊組；10-花崗斑巖；11-地質(zhì)界線；12-斷層及編號；13-礦段范圍；14-勘探線及編號；15-見礦鉆孔；16-未見礦鉆孔；17-礦體水平投影邊界線；18-水系；19-道路

1-中寒武統(tǒng)龍哈組；2-中寒武統(tǒng)田蓬組；3-下泥盆統(tǒng)坡松沖組；4-礦體及編號；5-鉆孔及編號；6-礦段范圍；7-斷層及編號；8-地質(zhì)界線

1) Ⅰ礦帶為銀、鉛、鋅、錫、銅等多組分同體共生礦，走向長約1.4 km，展布寬度30～800 m，處于63勘探線～210勘探線之間，賦存標(biāo)高為1 534～1 863 m，分布面積1.16 km2。Ⅰ礦帶中規(guī)模最大的礦體為Ⅰ1礦體，呈層狀連續(xù)產(chǎn)出，走向北西-南東向，傾向200°～240°，傾角5°～47°，鉆孔單工程見礦最厚為25.99 m，最薄為0.22 m,平均厚度為4.86 m。

2) Ⅴ礦帶沿走向長約1 000 m，水平展布寬度66～500 m，礦帶內(nèi)的礦體整體呈北西-南東向展布，傾向北東，礦體傾角較小，一般小于45°，深度越深傾角越緩，產(chǎn)狀基本與地層產(chǎn)狀一致，呈似層狀產(chǎn)出。礦帶內(nèi)大部分礦體均為單工程見礦，厚度品位較低，礦體厚度沿走向和傾向均較穩(wěn)定，各成礦元素變化穩(wěn)定程度均屬于較均勻。礦帶內(nèi)V礦體為主要礦體，走向北西-南東向，傾向55°～83°，傾角18°～55°，鉆孔單工程見礦最厚為20.83 m，最薄為0.47 m，平均厚度為5.18 m。

3) 銅礦帶中銀、鉛、鋅、錫組分含量較少，其沿走向方向長約660 m，水平展布寬度40～300 m，礦體傾角較小，呈似層狀產(chǎn)出。 銅1礦體為礦帶內(nèi)主要礦體，位于Ⅰ1礦體與Ⅰ2礦體之間，產(chǎn)狀與Ⅰ1礦體基本一致，呈似層狀、波狀產(chǎn)出，鉆孔單工程見礦最厚為19.26 m，最薄為2.59 m，平均真厚度為11.53 m。各礦體基本特征見表1。

表1 對門山礦段礦體特征Table 1 Ore body features in Duimenshan Mining Section

續(xù)表1

2 資源儲量估算過程

2.1 水平投影地質(zhì)塊段法

由于礦體主要受斷裂及其與背向斜虛脫空間復(fù)合控制，控礦斷裂破碎帶較連續(xù)完整，礦體呈似層狀、透鏡狀產(chǎn)出，平均傾角16°，產(chǎn)狀較為平緩，現(xiàn)有坑鉆工程分布相對較均勻，故采用水平投影地質(zhì)塊段法進(jìn)行儲量計(jì)算[10]。

水平投影地質(zhì)塊段法是在算數(shù)平均法基礎(chǔ)上改良而來，其原理是將圈出的礦體進(jìn)行水平投影，并在水平投影面上，按照工業(yè)類型，資源儲量級別等將礦體分為不同的理想形狀的板塊體塊段，然后分塊段算數(shù)平均計(jì)算儲量，累加后可得總儲量的傳統(tǒng)儲量計(jì)算方法[11]。

2.1.1 礦體圈定

本次礦體圈定以Pb、Zn為主混合圈定礦體，即在見礦單工程中從Pb、Zn任一組分等于或大于邊界品位的樣品一起圈入礦體，其他組分據(jù)其實(shí)際品位參加計(jì)算。礦體連接先連接斷層、地質(zhì)界線，再連接礦體；單斜礦體以直線連接，背向斜轉(zhuǎn)折部位用自然曲線連接；尚難利用礦與工業(yè)礦之間對角線連接，夾石尖滅于相鄰工程間距的1/2。

2.1.2 平均品位

單工程礦體平均品位由樣長與單樣品位加權(quán)平均求得，當(dāng)?shù)V體為多層礦時(shí)先用樣長加權(quán)計(jì)算礦層平均品位，再用礦層厚度加權(quán)計(jì)算單工程礦體平均品位。塊段礦體平均品位由塊段內(nèi)單工程礦體厚度加權(quán)平均求得，單礦體、礦段礦體平均品位用礦石量加權(quán)平均求得。

2.1.3 儲量估算

在以MapGIS軟件制作的1∶2 000資源量估算水平投影圖中用面積統(tǒng)計(jì)功能分別測定各塊段水平投影面積，以相應(yīng)塊段礦體平均傾角求得塊段的真面積，真面積與真厚度的乘積就是塊段的體積，與體重相乘可得塊段礦石量，最后對所有塊段進(jìn)行累加，得出礦段資源儲量。

2.2 三維建模距離冪次反比法資源儲量估算

3DMine三維礦業(yè)工程軟件是由國內(nèi)創(chuàng)新設(shè)計(jì)的三維礦業(yè)軟件平臺，也是目前國內(nèi)礦業(yè)公司應(yīng)用最廣泛的軟件之一。本次三維建模儲量估算通過3DMine三維礦業(yè)工程軟件建立對門山礦段三維地質(zhì)模型和實(shí)測工程模型。同時(shí)，以地質(zhì)統(tǒng)計(jì)學(xué)原理為基礎(chǔ)，對對門山礦段礦體進(jìn)行品位估值與資源儲量估算。 常用的地質(zhì)統(tǒng)計(jì)學(xué)方法包括克里格法、距離冪次反比法以及最近距離法等，由于礦段內(nèi)探礦工程及樣品數(shù)量滿足統(tǒng)計(jì)學(xué)要求，且有用組分分布較均勻的特點(diǎn)，采用距離冪次反比法進(jìn)行資源儲量估算[11-12]。

2.2.1 地質(zhì)數(shù)據(jù)庫及樣品統(tǒng)計(jì)分析

對收集到的基礎(chǔ)地質(zhì)數(shù)據(jù)資料，分別按照探礦工程孔口信息、測斜信息、樣品化驗(yàn)信息、巖性信息進(jìn)行分類整理，并將整理好的文件導(dǎo)入3DMine軟件中，分別保存為孔口文件、測斜文件、樣品文件、巖性文件。通過軟件的數(shù)據(jù)校驗(yàn)功能對其空間邏輯關(guān)系進(jìn)行校驗(yàn)，校驗(yàn)無誤后便可成功創(chuàng)建地質(zhì)數(shù)據(jù)庫。同時(shí)，利用統(tǒng)計(jì)分析功能對Ag、Sn、Pb、Zn和Cu元素的品位特征、分布特征等進(jìn)行分析，為后續(xù)有用元素品位的空間品位賦值提供基礎(chǔ)(圖3)。

圖3 對門山礦段Ag、Zn、Pb、Sn樣品分布直方圖Fig.3 Distribution histogram of Ag,Zn,Pb and Sn samples in Duimenshan Mining Section

2.2.2 礦體模型

礦體模型是三維地質(zhì)模型中最關(guān)鍵、最重要的部分，不僅能充分反映礦體在三維空間的形態(tài)特征、礦體與礦體之間的關(guān)系，還可為后續(xù)塊段模型的建立、元素品位賦值和資源量估算等提供約束條件。

對門山礦段礦體模型的創(chuàng)建，主要依據(jù)礦區(qū)范圍內(nèi)8條勘探線及礦區(qū)外圍探礦工程，先分勘探線逐個(gè)圈定礦體邊界，之后連接相鄰勘探線間礦體邊界形成封閉的礦體實(shí)體模型(圖4)。

圖4 對門山礦段Ⅰ1礦體實(shí)體模型圖Fig.4 Entity model of Ⅰ1 ore body in Duimenshan Mining Section(注：圖中數(shù)字為勘探線編號)

由于白牛廠礦床成因復(fù)雜，成礦元素多，無法明確區(qū)分各個(gè)礦體的主成礦元素，礦床屬于異體共生礦床，且Ag、Pb、Zn、Sn、Cu為同體共生礦，其有用組分無法分采，考慮其復(fù)雜性，本次圈定礦體時(shí)選擇礦石量混圈和單元素單礦種單圈兩種方式。即在礦石量的估算過程中，按照四個(gè)成礦元素混圈方法進(jìn)行，即某個(gè)元素品位到達(dá)工業(yè)或邊界時(shí)，均圈入礦體計(jì)算礦石量，在金屬量估算中，見礦單工程中從對Ag、Sn、Pb、Zn、Cu任一組分等于或大于邊界品位的樣品單獨(dú)圈入礦體，視各個(gè)成礦元素為獨(dú)立礦種，不考慮共伴生關(guān)系，從而方便各個(gè)成礦元素金屬量估算工作。

2.2.3 塊段模型創(chuàng)建及品位推估

礦體模型只能反映礦體的空間形態(tài)，而無法表現(xiàn)出礦體局部品位和巖性及其分布情況，也無法估算，因此必須建立礦體的塊段模型。 塊段模型是品位估值和資源量計(jì)算的基礎(chǔ)，其基本思想是把礦體在三維空間內(nèi)按照一定的尺寸劃分成許多小塊段，然后對填滿整個(gè)礦體范圍內(nèi)的單元塊的品位根據(jù)已知的樣品進(jìn)行推估，并在此基礎(chǔ)上進(jìn)行資源量的計(jì)算[13]。

塊段模型創(chuàng)建應(yīng)先進(jìn)行空白塊段的模型建立，空白塊段模型必須反映礦床的空間分布范圍，即平面范圍和深度范圍。本次建立的空白塊段模型以覆蓋所有礦體為原則，即平面范圍和深度范圍均能夠完全覆蓋所有礦體?；A(chǔ)塊尺寸的確定主要依據(jù)礦床規(guī)模、空間形態(tài)、勘探網(wǎng)度和擬采用的開采方式。在綜合白牛廠礦區(qū)礦體形態(tài)和特征基礎(chǔ)上，確定空白塊段模型的基礎(chǔ)尺寸為2 m×2 m×1 m，礦體內(nèi)部次級分塊尺寸為1 m×1 m×0.5 m。

采用的計(jì)算方法是距離冪次反比法，它是最常用的空間內(nèi)插方法之一，是一種與空間距離有關(guān)的插值方法，在計(jì)算插值點(diǎn)取值時(shí)按距離越近權(quán)重值越大的原則，用若干臨近點(diǎn)的線性加權(quán)來擬合估計(jì)點(diǎn)的值[14-15]。

搜索橢球體參數(shù)根據(jù)白牛廠礦床礦體形態(tài)特點(diǎn)，礦體走向較長，傾向次之，厚度最小，則橢球體三軸的半徑比為，東方向∶北方向∶高程方向=1∶1∶0.5。根據(jù)礦床勘探類型并結(jié)合礦山實(shí)際要求，以40 m、80 m、200 m的橢球搜索半徑對探明資源量、控制資源量、推斷資源量進(jìn)行資源儲量分類并進(jìn)行品位推估[16]。

3 資源儲量估算對比分析

3.1 計(jì)算方法及特點(diǎn)對比

水平投影地質(zhì)塊段法在儲量計(jì)算過程中，以現(xiàn)有坑鉆工程將礦體劃分為簡單的幾何形態(tài)塊段，從而進(jìn)行資源儲量計(jì)算，其計(jì)算公式見式(1)和式(2)。

Q=V×D

(1)

式中：Q為礦石量；V為體積；D為體重。

P=Q×C

(2)

式中：P為金屬量；Q為礦石量；C為品位。

上述算法優(yōu)點(diǎn)在于方法簡單，易于掌握，在計(jì)算過程中，作圖及計(jì)算都較為簡易便捷。但其缺點(diǎn)也較為明顯，若工作區(qū)礦體較為復(fù)雜，坑鉆工程較少，控制程度較低，則其計(jì)算結(jié)果誤差較大，且在礦體塊段劃分時(shí)，不能將礦體的空間信息考慮進(jìn)去，從而因塊段的不同劃分方式導(dǎo)致其資源儲量估算結(jié)果產(chǎn)生誤差。同時(shí)，在生產(chǎn)建設(shè)后期，其利用性差，隨著探采工程的加深，重新計(jì)算儲量時(shí)，需重繪儲量估算圖件，重新計(jì)算儲量[17]。

地質(zhì)統(tǒng)計(jì)學(xué)法資源儲量估算是一種新型資源儲量估算方法，它結(jié)合了三維礦業(yè)工程軟件與地質(zhì)學(xué)統(tǒng)計(jì)學(xué)原理進(jìn)行資源儲量計(jì)算，其計(jì)算公式見式(3)。

(3)

式中：g為估計(jì)值；gi為第n個(gè)樣本；Di為距離；p為距離的冪次。

相對于水平投影地質(zhì)塊段法，其計(jì)算方法依托地質(zhì)三維模型進(jìn)行計(jì)算，而三維模型能精細(xì)表征礦區(qū)礦體和礦區(qū)工程。此方法運(yùn)用可視化思維將礦體一目了然地展示出來，更方便可視化操作，以及為后續(xù)的工程開發(fā)和勘探提供更直觀的支持。同時(shí)，在計(jì)算過程中，水平投影地質(zhì)塊段法將礦體塊段視為理想、品位均勻的板塊體進(jìn)行計(jì)算，且人為進(jìn)行塊段劃分時(shí)，易產(chǎn)生較大誤差，地質(zhì)統(tǒng)計(jì)學(xué)方法則根據(jù)已知的樣品對一定距離范圍內(nèi)的小又多的單元塊進(jìn)行品位估值，其距離范圍通過搜索半徑得到。相較于傳統(tǒng)資源儲量計(jì)算方法而言，該方法通過搜索橢球體半徑來進(jìn)行品位推估，更為合理和精準(zhǔn)，結(jié)果也更接近真實(shí)的礦體狀態(tài)。在礦山后續(xù)生產(chǎn)建設(shè)中，地質(zhì)三維模型可繼續(xù)完善使用，不必推倒重來，后續(xù)建設(shè)時(shí)可利用性強(qiáng)。但其操作較為繁瑣，在坑探工程保持一定數(shù)量基礎(chǔ)上，還需要考慮到參與礦體儲量估算的工程之間的相關(guān)性。計(jì)算中搜索橢球體半徑這一操作，也使得資源儲量類別界限不明顯，其參數(shù)設(shè)置顯著影響著估算結(jié)果的準(zhǔn)確性[18]。

3.2 水平投影面積對比

本文綜合考慮礦區(qū)復(fù)雜情況，以及兩種資源儲量估算方法對礦體的圈定原則的不同，在水平投影地質(zhì)塊段法中，以Pb、Zn元素為主，及兩者中任一元素大于邊界品位及圈入礦體；但在地質(zhì)統(tǒng)計(jì)學(xué)法資源儲量估算中，更充分地考慮到礦床的復(fù)雜性和真實(shí)性，采用Ag、Pb、Zn、Sn四種成礦元素混圈。無論選用何種方法進(jìn)行資源儲量計(jì)算，礦石量的計(jì)算就是礦體體積的計(jì)算，而在礦體產(chǎn)狀已知的情況下，其投影面積就成為了資源儲量計(jì)算中決定的因素。就本文兩種資源儲量計(jì)算方法而言，采用礦體水平投影面積對比更為直觀，礦區(qū)內(nèi)主要礦體水平投影面積對比見表2。

表2 水平投影面積對比Table 2 Horizontal projection area comparison table

經(jīng)對比，兩種資源儲量計(jì)算方法的礦體水平投影面積接近，相對誤差值較小，可認(rèn)為兩種計(jì)算方法所圈定礦體得形態(tài)大小基本保持一致。

3.3 品位對比

在資源儲量估算過程中，除了礦石量的估算，品位的推估也尤為重要。水平投影地質(zhì)塊段法應(yīng)用了加權(quán)平均數(shù)的計(jì)算方法，逐一由樣品加權(quán)到單工程加權(quán)，最后到塊段加權(quán)。地質(zhì)統(tǒng)計(jì)學(xué)法資源儲量估算則應(yīng)用空間差值法-距離冪次反比法進(jìn)行品位推估。兩者關(guān)于品位的計(jì)算方法大不相同，所以進(jìn)行礦區(qū)成礦金屬的品位對比，能相互應(yīng)照，反映兩種計(jì)算方法的真實(shí)性和準(zhǔn)確性，其結(jié)果對比見表3。 由表3可知，兩者相對誤差值較小，品位推估結(jié)果較吻合。

表3 礦段平均品位對比Table 3 Comparison of the average grade of ore section

3.4 礦石量及金屬量對比

資源儲量估算的最終目的是計(jì)算出礦石量和金屬量。水平投影地質(zhì)塊段法中礦石量的估算是將礦體視為多個(gè)理想板塊體進(jìn)行簡便化體積運(yùn)算，金屬量則是由各塊段的體積與塊段的金屬加權(quán)品位相乘得出。依托三維建模的地質(zhì)統(tǒng)計(jì)學(xué)法儲量估算則直接計(jì)算建模出的礦體體積完成礦石量估算，由于礦區(qū)礦體的復(fù)雜性，其金屬量的估算則是視Ag、Pb、Zn、Sn、Cu五種成礦元素為獨(dú)立礦種，單獨(dú)圈定來完成金屬量的估算。由表4可知，兩種方法的礦石量相對誤差小于5%，金屬量估算結(jié)果中僅Ag的相對誤差略大，其原因可能是礦段中Ag金屬量過少，其余四種金屬誤差均在可控范圍內(nèi)，可認(rèn)為兩種資源儲量估算結(jié)果可靠。

表4 礦石量、金屬量結(jié)果對比表Table 4 Comparison table of ore and metal content results

3.5 消耗量與采礦量對比

在保證新舊兩種資源儲量估算方法結(jié)果可靠的基礎(chǔ)上，進(jìn)行估算結(jié)果與采礦量的對比分析，以期驗(yàn)證兩種方法的精準(zhǔn)性。本次選用資源儲量估算的消耗量與礦山采礦量進(jìn)行對比，由于采礦過程中保留了較多的安全礦柱，對比時(shí)使用的采礦量數(shù)據(jù)實(shí)為選礦廠礦石量數(shù)據(jù)及礦柱相加的總量，其中礦柱重量為實(shí)測礦柱體積與礦石體重的乘積(表5)。由表5可知，兩種方法估算結(jié)果與采礦量的相對誤差均小于15%，但地質(zhì)統(tǒng)計(jì)學(xué)方法所估算的數(shù)值更接近真實(shí)數(shù)據(jù)，精確度更高。

表5 消耗量與采礦量對比表Table 5 Comparison table of consumption and mining volume

4 結(jié) 論

1) 本文綜合使用傳統(tǒng)的水平投影地質(zhì)塊段法和3DMine三維礦業(yè)工程軟件進(jìn)行地質(zhì)建模，并使用地質(zhì)統(tǒng)計(jì)學(xué)法中的距離冪次反比法對云南省蒙自白牛廠銀多金屬礦區(qū)對門山礦段進(jìn)行了礦區(qū)金屬品位估值及資源儲量估算。估算過程中，水平投影地質(zhì)塊段法相對簡便，操作簡單，其準(zhǔn)確結(jié)果和工程控制程度成正相關(guān)，但后續(xù)利用性差，僅能提供二維數(shù)據(jù)。而地質(zhì)統(tǒng)計(jì)學(xué)方法依托三維建模進(jìn)行運(yùn)算，其操作相對復(fù)雜，參數(shù)的設(shè)置對結(jié)果影響較大，但其可視化程度較高，可為礦山建設(shè)提供精度較高的三維數(shù)據(jù)，且只需完善數(shù)據(jù)，就能持續(xù)為后續(xù)的礦山建設(shè)提供支持。

2) 通過對兩種資源儲量估算方法的計(jì)算方法、礦體水平投影面積、成礦金屬平均品位、礦石量以及金屬量進(jìn)行詳細(xì)對比，其誤差基本在5%以內(nèi),Ag由于金屬量較小，計(jì)算誤差略大，但也小于7%，可認(rèn)為使用地質(zhì)統(tǒng)計(jì)學(xué)方法估算資源儲量較為準(zhǔn)確可靠。在保證地質(zhì)統(tǒng)計(jì)學(xué)方法資源儲量估算可靠性的前提下，再次進(jìn)行了與實(shí)際采礦量的分析對比，兩種方法與實(shí)際采礦量的相對誤差都小于15%，但地質(zhì)統(tǒng)計(jì)學(xué)方法所得結(jié)果明顯更接近實(shí)際采礦量，數(shù)據(jù)更為精準(zhǔn)。

3) 就目前兩種方法的對比分析來看，兩種資源儲量估算方法各有不足之處，但新老兩種方法之間的相互驗(yàn)證是一個(gè)必要環(huán)節(jié)，不僅可以驗(yàn)證地質(zhì)統(tǒng)計(jì)學(xué)方法的可靠性，更可以促進(jìn)其完善和推廣應(yīng)用。在本礦段資源儲量估算過程中，傳統(tǒng)幾何學(xué)方法雖然操作簡便且可以滿足生產(chǎn)要求，但在礦山的動態(tài)發(fā)展生產(chǎn)過程中，其可利用性差，三維可視性低，適用性唯一，遠(yuǎn)遠(yuǎn)不如地質(zhì)統(tǒng)計(jì)學(xué)方法，且在本文對比中，地質(zhì)統(tǒng)計(jì)學(xué)方法估算的資源儲量結(jié)果明顯更為接近真實(shí)數(shù)據(jù)。綜上，地質(zhì)統(tǒng)計(jì)學(xué)方法資源儲量估算優(yōu)于傳統(tǒng)幾何法資源儲量估算，更適合礦山后續(xù)的生產(chǎn)設(shè)計(jì)。