朱曉冰 睢棟超 張超超
摘 要:固體礦產(chǎn)資源儲量計算是地質(zhì)礦產(chǎn)評價中一項重要的工作,對于大型、特大型礦床中海量的化學樣數(shù)據(jù),如果使用傳統(tǒng)計算方式無疑工作量巨大、耗時超長、容易出錯,如果有關參數(shù)發(fā)生變化,還要重新檢查計算。本文基于Excel VBA為開發(fā)工具,針對鋁土礦及其共生礦產(chǎn),本文主要闡述計算單工程化學樣平均厚度、平均品位的計算原理、流程圖,并在實際工作中得到檢驗和應用,體現(xiàn)了高效率、高準確度、且靈活性強、對原始數(shù)據(jù)格式要求不高的特點。且對其他礦種資源儲量估算也有利用或借鑒的價值。
關鍵詞:資源儲量估算;化學樣;程序設計;Excel VBA;鋁土礦
0 引言
河南省澠池縣曹窯以西煤下鋁工作區(qū)位于河南省三門峽市陜州區(qū)--澠池縣一帶,總面積約51平方公里。本工作區(qū)從2011年至今,河南省有色金屬地質(zhì)礦產(chǎn)局第六地質(zhì)大隊先后提交了《河南省曹窯以西煤下鋁普查報告》、《河南省曹窯以西煤下鋁五門溝礦段詳查報告》、《河南省曹窯以西煤下鋁柴洼礦段詳查報告》、《河南省曹窯以西煤下鋁關底沃-扣門山礦段詳查報告》和《河南省曹窯以西煤下鋁青陽溝礦段詳查報告》等五個報告,累計施工鉆孔451個,進尺131190.51m,基本分析樣3551個。
另收集以往工作中及鄰區(qū)已施工的鉆孔75個,收集596個基本化學分析樣資料。本次工作和以往工作的鉆孔合計526個,基本化學樣4147個。
面對如此巨大的數(shù)據(jù)量,過去地勘單位估算資源量/儲量時,一般采用計算器來計算,計算周期長,計算過程復雜、繁瑣,效率極低、容易出錯,如果中間結(jié)果出錯或者需要修改,后面結(jié)果推倒重來,往往形成一個報告需要數(shù)月、甚至一兩年才能完成。報告完成人員工作強度大、任務重;隨著計算機開始逐漸普及,地勘單位人員開始使用計算機軟件如WPS、Excel等處理數(shù)據(jù),但水平不高,一般采用簡單的公式和函數(shù)計算,效率提升不多,計算過程仍然存在繁瑣、任務重、易出錯等缺點。因此,編寫一種當前滿足鋁土礦資源儲量報告估算需要,并且實現(xiàn)自動化,使計算過程工作高效、準確成為當務之急。
1 地質(zhì)特征
1.1 礦床地質(zhì)
本區(qū)鋁土礦含礦巖系為石炭系上統(tǒng)本溪組,鋁土礦礦體形態(tài)多呈層狀、似層狀或透鏡狀。礦床類型為沉積型,總體為特大型鋁土礦。本區(qū)無礦層露頭,均為隱伏礦體。
本區(qū)鋁土礦共生礦種有粘土礦(分上下層,個別報告還要求計算夾層粘土礦)、菱鐵礦、赤鐵礦、鐵礬土等。另有鎵、氧化鋰、輕稀土等伴生元素或有益組分。
化學分析樣品中鋁土礦基本分析項目為Al2O3、SiO2、Fe2O3、S、TiO2、燒失量(LOSS);共生礦種基本分析項目與鋁土礦基本相同,粘土礦另要求分析耐火度(SK)。鋁土礦及其共生礦種之間的區(qū)分除了以礦石的顏色、結(jié)構(gòu)、構(gòu)造特征區(qū)分外,還主要以各個分析項目品位的多少區(qū)分,如鋁土礦邊界品位要求Al2O3≥40%,A/S(鋁硅比)≥1.8;Ⅲ級硬質(zhì)粘土礦品位要求Al2O3≥40%,F(xiàn)e2O3≥3.5%。
該區(qū)鋁土礦體均為層狀或似層狀,傾角平緩(5~16度),探礦工程按勘探線布置,工程分布比較均勻,為便于生產(chǎn)利用,資源量估算選用地質(zhì)塊段法,在水平投影圖上進行。
1.2 單工程圈定原則
依據(jù)各單工程中樣品分析結(jié)果,凡單樣品位等于或大于邊界品位要求者均圈入礦體內(nèi);夾石樣品厚度小于相應夾石剔除厚度時并入相應礦體中。
2 程序設計
2.1 程序設計基礎
本人結(jié)合工作實際情況,經(jīng)過長期摸索,利用微軟公司的Excel VBA,編寫了計算鋁土礦單工程厚度、平均品位的程序,以及計算塊段平均厚度、平均品位的程序,用于該工作區(qū)的資源量估算工作。
Microsoft Excel是美國微軟公司編寫的一款電子表格軟件。界面直觀,有非常出色的計算功能和圖表工具,功能強大,普及率高;VBA是內(nèi)嵌于office軟件中的一個開發(fā)模塊,這個模塊提供程序自主開發(fā),語言基礎和Visual Basic相似。在Excel的“工具--宏--編輯器”可以打開。
2.2 原始數(shù)據(jù)格式要求
本程序?qū)υ紨?shù)據(jù)格式要求不高,只需要對每個鉆孔的原始化驗結(jié)果按從上到下的順序排列,不需要單獨提取含礦樣品。在Excel中表中A到V列對所有鉆孔樣品數(shù)據(jù)按表1格式整理。
2.3 程序設計原理及流程圖
為計算加權(quán)平均值(平均品位),為此特別編寫了一個加權(quán)平均自定義函數(shù)WeightedAverage,具體實現(xiàn)原理和代碼見本人的相關參考論文。
對于大批量工程的計算單工程平均品位、平均厚度的程序流程圖如圖1:
例如為了計算工程ZK4016的鋁土礦厚度和平均品位,依次對ZK4016的每個樣品進行比較,因為ZK4016的1號樣品為炭質(zhì)頁巖,非鋁土礦,忽略后比較下一個樣品;2號樣品為鋁土礦,記錄ZK4016的鋁土礦excel單元格地址為$L$657,然后依次比較3號樣品,非鋁土礦,忽略;4號樣品是鋁土礦,且與2號鋁土礦樣品的夾層真厚度0.85m(3號樣品的真厚度)≥0.8m(最大夾石剔除厚度),故夾石不能合并到鋁土礦中,此時記錄ZK4016的鋁土礦樣品地址為$L$657,$L$659,然后依次比較5至11號樣品,記錄ZK4016的鋁土礦地址為$L$657,$L$659, $L$661:$L$664。由于$L$657的真厚度為0.62m,≤最低可采厚度(0.80m),故應忽略,最終ZK4016的鋁土礦樣品地址為$L$659, $L$661:$L$664,然后依次計算ZK4016的樣品編號、各種厚度、單工程平均品位。真厚度的公式為“=SUM($G$659,$G$661:$G$664)”,A12O3的單工程平均品位公式為“=ROUND(weightedaverage($M$659,$M$661:$M$664,-5),2)”,其中參數(shù)“-5”為真厚度所在列(G列)相對于A12O3品位所在列(M列)的偏移列數(shù)值。ROUND為Excel中的四舍五入函數(shù)。
2.4 程序具體實現(xiàn)過程
為了在excel表格中記錄鋁土礦采樣位置、樣品編號、地址編寫了一個結(jié)構(gòu)體:
Public Type OreLayer
LayerFirstSampleRow&‘記錄每鋁土礦層的第一個樣品單元格所在行
LayerLastSampleRow&‘記錄每鋁土礦層的最后一個樣品單元格所在行
LayerSampleIDs$‘記錄單工程鋁土礦樣品編號(包含一層或多層)
LayerTopPosition#‘記錄每鋁土礦層的第一個樣品的采樣位置
LayerBottomPosition#‘記錄每鋁土礦層的最后一個樣品的采樣位置
LayerThickness#‘記錄每鋁土礦層的視厚度
LayerSampleAddress$‘記錄每個工程中鋁土礦樣品的單元格地址
LayerTrueThickness#‘記錄每鋁土礦層的真厚度
End Type
為了方便直觀顯示計算每個鉆孔的鋁土礦樣長、品位等信息,在excel的基本樣品分析結(jié)果表的W到AF列依次計算出鉆孔的樣號、視厚度、鉛直厚度、Al2O3(%)品位、SiO2(%)品位、Fe2O3(%)品位、S(%)品位、LOSS(%)等,并按單工程合并單元格。最終計算結(jié)果如表2(由于頁面太窄,為了方便顯示,C:V列隱藏)。
由于源代碼比較長,限于篇幅無法顯示,有興趣的讀者可以與本人聯(lián)系。
3 結(jié)論
所有代碼在Excel2016環(huán)境下全部編譯通過,并且在本礦區(qū)先后提交的五個報告經(jīng)過實際運用和檢驗,數(shù)據(jù)計算快速、準確、效率高,其中的關底沃--扣門山礦區(qū)的112個鉆孔的運算時間不到1秒,遠遠超過傳統(tǒng)的計算方法。
經(jīng)檢驗,本程序在計算超大量數(shù)據(jù)時快速、準確,體現(xiàn)了高效率、高準確度、靈活性強。得到了甲方單位和報告評審單位人員的認可和好評,值得進一步推廣應用。在貴金屬、有色金屬、非金屬計算資源儲量中也有利用和借鑒的價值。
參考文獻:
[1]何文平,董運如等.河南省澠池縣曹窯以西煤下鋁關底沃-扣門山礦段詳查報告[M].鄭州:河南省國土資源科學研究院,2008.
[2]朱曉冰.資源量、儲量估算中求加權(quán)平均值方法探討[J].華北國土資源,2015(5):69-72.
作者簡介:
朱曉冰(1979- ),男,2003年畢業(yè)于焦作工學院(現(xiàn)河南理工大學),工程師,現(xiàn)從事地質(zhì)礦產(chǎn)勘查工作。