基于Excel VBA在鋁土礦勘查中快速計算單工程化學樣平均品位、厚度的方法

朱曉冰 睢棟超 張超超

摘 要：固體礦產(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)勘查工作。