雷顯權(quán)，莫曉東

(長(zhǎng)沙有色冶金設(shè)計(jì)研究院有限公司，　湖南 長(zhǎng)沙　410011)

?

基于DIMINE軟件的巖溶堆積型鋁土礦建模方法

雷顯權(quán)，莫曉東

(長(zhǎng)沙有色冶金設(shè)計(jì)研究院有限公司，湖南 長(zhǎng)沙410011)

摘要：介紹了利用DIMINE軟件構(gòu)建巖溶堆積型鋁土礦數(shù)字化模型的方法。數(shù)字化模型主要包括地質(zhì)數(shù)據(jù)庫(kù)、礦體實(shí)體模型和礦體塊段模型三部分。首先收集探礦工程資料，整理導(dǎo)入軟件形成地質(zhì)數(shù)據(jù)庫(kù)。根據(jù)地質(zhì)數(shù)據(jù)庫(kù)解譯礦層頂、底板，結(jié)合地形地質(zhì)圖上的礦體邊界，采用頂?shù)装彘]合法構(gòu)建礦體三維實(shí)體模型，并采用距離冪次反比法對(duì)塊段模型進(jìn)行屬性賦值，估算模型資源儲(chǔ)量，由此實(shí)現(xiàn)地質(zhì)資源數(shù)字化。

關(guān)鍵詞：數(shù)字化礦山；DIMINE軟件；巖溶堆積型鋁土礦；地質(zhì)建模

巖溶堆積型鋁土礦是原生沉積型鋁土礦經(jīng)巖溶風(fēng)化形成的次生鋁土礦床，主要賦存于巖溶洼地的第四系下更新統(tǒng)地層中，具有礦床規(guī)模大、礦石質(zhì)量?jī)?yōu)等特點(diǎn)[1]。該類型礦床主要分布于桂西北地區(qū)，近年在河南澠池、魯山等地也有發(fā)現(xiàn)[2]。由于礦體分布面廣、規(guī)模差別大、平面形態(tài)復(fù)雜，采用傳統(tǒng)方法進(jìn)行資源管理和制定礦山規(guī)劃難度大。近年來，數(shù)字化礦山技術(shù)發(fā)展迅速，為礦山實(shí)現(xiàn)信息化、自動(dòng)化、智能化管理提供了可能[3-4]。DIMINE軟件具有操作便利的特點(diǎn)[5]。因此，本文采用DIMINE軟件創(chuàng)建礦體模型，以實(shí)體模型為約束，以地質(zhì)數(shù)據(jù)庫(kù)為基礎(chǔ)對(duì)塊段模型進(jìn)行估值和資源儲(chǔ)量估算，由此實(shí)現(xiàn)地質(zhì)資源數(shù)字化，建模流程見圖1。

圖1　建模流程

1地質(zhì)數(shù)據(jù)庫(kù)

巖溶堆積型鋁土礦的探礦手段主要為淺井工程。系統(tǒng)收集地質(zhì)探礦報(bào)告等資料，按照DIMINE軟件對(duì)數(shù)據(jù)庫(kù)文件的格式要求，將原始工程數(shù)據(jù)整理為3個(gè)電子表格：孔口表、測(cè)斜表和樣品表(巖性表可以空缺)?？卓诒戆降V工程的地表位置和深度等信息，測(cè)斜表包含探礦工程在不同深度上的偏斜測(cè)量結(jié)果，樣品表包含探礦工程的采樣信息和分析化驗(yàn)結(jié)果。每個(gè)表格有特定的字段，每個(gè)字段有具體的類型、含義和要求。對(duì)各表的格式和內(nèi)容詳細(xì)說明如下。

(1) 孔口表字段說明(見表1)：“工程編號(hào)”為淺井編號(hào)，字符串型；“X坐標(biāo)”為淺井東坐標(biāo)，雙精度型，m；“Y坐標(biāo)”為淺井北坐標(biāo)，雙精度型，m；“Z坐標(biāo)”為井口高程，雙精度型，m；“井深”為淺井深度，雙精度型，m；“礦體號(hào)”為礦體編號(hào)，字符串型；“勘探線號(hào)”為勘探線編號(hào)，字符串型；“采場(chǎng)號(hào)”為采場(chǎng)編號(hào)，字符串型；“礦頂標(biāo)高”為礦體頂板標(biāo)高，雙精度型，m；“礦底標(biāo)高”為礦體底板標(biāo)高，雙精度型，m。

表1　孔口表格式

(2) 測(cè)斜表字段說明(見表2)：“工程編號(hào)”為淺井編號(hào)；“測(cè)斜深度”為測(cè)斜位置深度，雙精度型，m；“方位角”為測(cè)得方位角，雙精度型，(°)；“傾角”為測(cè)得傾角，雙精度型，(°)。

表2　測(cè)斜表格式

(3) 樣品表字段說明(見表3)：“工程編號(hào)”字段為淺井編號(hào)；“自”為采樣起始位置深度，雙精度型，m；“至”為采樣結(jié)束位置深度，雙精度型，m；“Al2O3”字段為Al2O3品位，雙精度型，%；“SiO2”為SiO2品位，雙精度型，%；“Fe2O3”為Fe2O3品位，雙精度型，%；“S”為硫品位，雙精度型，%；“A/S”為鋁硅比，雙精度型；“A/F”為鋁鐵比，雙精度型；“密度”為原礦體重，雙精度型，t/m3；“含礦率”為凈礦含礦率，雙精度型，kg/m3；“巖性”為樣品巖性，字符串型。

將整理好且校驗(yàn)無誤的孔口表、測(cè)斜表和樣品表導(dǎo)入DIMINE軟件生成“.dmd”格式數(shù)據(jù)庫(kù)，然后將“.dmd”格式數(shù)據(jù)庫(kù)轉(zhuǎn)換為“.dmg”格式數(shù)據(jù)庫(kù)，即為所建立的地質(zhì)數(shù)據(jù)庫(kù)。由于3個(gè)表格共享“工程編號(hào)”字段，數(shù)據(jù)庫(kù)各表之間便可實(shí)現(xiàn)數(shù)據(jù)對(duì)應(yīng)銜接。

表3　樣品表格式

2礦體實(shí)體模型

針對(duì)巖溶堆積型鋁土礦的特殊層狀礦體特征，采用頂、底板閉合法構(gòu)建礦體三維實(shí)體模型。具體步驟如下：

(1) 利用地質(zhì)數(shù)據(jù)庫(kù)分別創(chuàng)建見礦淺井的礦層頂板點(diǎn)和底板點(diǎn)；將地形地質(zhì)圖上的零點(diǎn)邊界線和灰?guī)r邊界線按規(guī)則連接，得到礦體邊界線；如果礦體存在天窗，還需提取地形地質(zhì)圖上的天窗邊界線。

(2) 利用實(shí)體建模模塊的區(qū)域DTM功能，選定礦體頂板點(diǎn)和礦體邊界線，以“距離冪”方式生成礦體頂板DTM面片。同理，選定礦體底板點(diǎn)和礦體邊界線生成礦體底板DTM面片。選定礦體頂板DTM面片，通過“提取開口線”命令生成礦體頂板邊界線，然后刪除礦體頂板DTM面片。同理，選定礦體底板DTM面片生成礦體底板邊界線，然后刪除礦體底板DTM面片。

(3) 利用線編輯模塊的定比例內(nèi)插線功能，選定礦體頂板邊界線和礦體底板邊界線，按0.5的比例生成一條多段線，即為新的礦體頂板邊界線，然后刪除之前的礦體頂板邊界線和礦體底板邊界線。利用線編輯模塊的移動(dòng)復(fù)制功能，選定新的礦體頂板邊界線，按XYZ分量模式以Z=-0.01復(fù)制生成一條多段線，即為新的礦體底板邊界線。

(4) 利用實(shí)體建模模塊的區(qū)域DTM功能，選定礦體頂板點(diǎn)和新的礦體頂板邊界線，以“曲線求交”方式生成新的礦體頂板DTM面片。同理，選定礦體底板點(diǎn)和新的礦體底板邊界線，生成新的礦體底板DTM面片。

(5) 利用實(shí)體建模模塊的線框功能，選定新的礦體頂板邊界線和新的礦體底板邊界線，以“最小周長(zhǎng)法”方式生成礦體邊界DTM面片。選定新的礦體頂板、底板和邊界DTM面片，通過“合并”命令生成礦體實(shí)體模型。

(6) 如果存在天窗，則利用實(shí)體建模模塊的裁剪功能，選定礦體實(shí)體模型和天窗邊界線，以“封閉裁剪”方式將天窗部分的實(shí)體裁剪并刪除，保留部分即為挖掉天窗之后的礦體實(shí)體模型。如有多個(gè)天窗，則照此重復(fù)操作。

(7) 仔細(xì)檢查礦體實(shí)體模型，如果存在不封閉處或不合理的三角面片連接方式，則需手動(dòng)進(jìn)行調(diào)整，直到實(shí)體有效性檢測(cè)獲得通過。

上述步驟是構(gòu)建礦體實(shí)體模型的一般流程，所有礦體實(shí)體模型構(gòu)建方式基本相同，只是不同礦體的建模復(fù)雜程度可能存在差異。

3礦體塊段模型

將礦體實(shí)體模型導(dǎo)入DIMINE軟件視窗，利用地質(zhì)模塊的空塊功能，通過“從屏幕上自動(dòng)獲取范圍”的方式來創(chuàng)建空塊段模型，起點(diǎn)坐標(biāo)和延伸長(zhǎng)度可自動(dòng)獲取，基礎(chǔ)塊尺寸設(shè)為X=20、Y=20、Z=1。

對(duì)塊段模型添加“礦巖類型”、“體重”、“Al2O3”、“SiO2”、“Fe2O3”、“密度”、“含礦率”等屬性字段，利用礦體實(shí)體模型過濾地質(zhì)數(shù)據(jù)庫(kù)，將過濾后的地質(zhì)數(shù)據(jù)庫(kù)以組合樣長(zhǎng)為1進(jìn)行樣品組合，得到樣品組合文件。創(chuàng)建搜索橢球體，橢球長(zhǎng)半軸、次半軸、短半軸分別設(shè)為25、25、5，橢球方位角、傾伏角、傾角都為0°。

采用距離冪次反比法對(duì)塊段模型進(jìn)行估值。塊段模型選擇上述空塊段模型文件；估值約束設(shè)為礦體實(shí)體模型內(nèi)部，調(diào)整內(nèi)部級(jí)數(shù)使單元塊尺寸為X=20、Y=20、Z=1，外部級(jí)數(shù)比內(nèi)部級(jí)數(shù)大1；樣品文件選擇上述樣品組合文件；橢球體參數(shù)選擇上述搜索橢球體，其它估值參數(shù)為缺省值；估值變量勾選“Al2O3”、“SiO2”、“Fe2O3”、“密度”、“含礦率”等字段。

參數(shù)設(shè)置完畢，以單方案估值方式進(jìn)行估值。第一次估值結(jié)束，如果存在未估值塊，則需將搜索半徑擴(kuò)大一倍進(jìn)行第二次估值，以此類推，直至未估值塊數(shù)為0。每次重新估值時(shí)，需將原來的值覆蓋。估值完成后，可利用DIMINE軟件將塊段模型數(shù)據(jù)導(dǎo)出電子表格文件，以便進(jìn)行檢核和二次調(diào)用。

4資源儲(chǔ)量估算

利用地質(zhì)模塊的儲(chǔ)量計(jì)算功能進(jìn)行資源儲(chǔ)量估算。選擇需要估算儲(chǔ)量的礦體實(shí)體模型，輸入相應(yīng)的已完成估值的塊段模型，內(nèi)部級(jí)數(shù)和邊界級(jí)數(shù)設(shè)為與塊段估值時(shí)一致，其它參數(shù)分別設(shè)為：體重字段為“密度”，主統(tǒng)計(jì)字段為“Al2O3”，其它統(tǒng)計(jì)字段為“SiO2”、“Fe2O3”、“含礦率”。估算完成后可得到資源儲(chǔ)量統(tǒng)計(jì)表，表內(nèi)有詳細(xì)的礦體體積、體重、礦石量、品位、含礦率等數(shù)據(jù)。此外，還可通過設(shè)置高程區(qū)間、品位區(qū)間或其它分類字段來對(duì)資源儲(chǔ)量進(jìn)行分類統(tǒng)計(jì)。

仔細(xì)檢核資源儲(chǔ)量統(tǒng)計(jì)表內(nèi)的各項(xiàng)數(shù)據(jù)，如有錯(cuò)誤或不合理之處，則需重新進(jìn)行資源儲(chǔ)量估算或塊段模型賦值，直至結(jié)果合理無誤。將通過模型估算得到的相關(guān)參數(shù)與地質(zhì)報(bào)告進(jìn)行對(duì)比，若偏差較大，則需仔細(xì)查找原因并作出調(diào)整，直到滿足規(guī)范要求。

5數(shù)據(jù)和模型更新

模型建好之后，如果有新增探礦工程，則需更新地質(zhì)數(shù)據(jù)庫(kù)、礦體實(shí)體模型和礦體塊段模型。將新的探礦工程資料整理成與原數(shù)據(jù)庫(kù)表格格式一致的孔口表、測(cè)斜表和樣品表，通過地質(zhì)模塊的數(shù)據(jù)庫(kù)更新功能將整理好的表格導(dǎo)入原地質(zhì)數(shù)據(jù)庫(kù)，便可得到新的地質(zhì)數(shù)據(jù)庫(kù)。更新礦體實(shí)體模型有2條途徑：如果新增工程比較少，可通過對(duì)原實(shí)體模型進(jìn)行局部調(diào)整以得到新的實(shí)體模型；如果新增工程很多，對(duì)原實(shí)體模型進(jìn)行局部調(diào)整的工作難度大，則需按照與原實(shí)體模型相同的建模流程來重新構(gòu)建新的實(shí)體模型。地質(zhì)數(shù)據(jù)庫(kù)和礦體實(shí)體模型更新之后，需對(duì)礦體塊段模型進(jìn)行重新賦值，賦值方法與原塊段模型賦值相同，只是在設(shè)置樣品文件和模型約束時(shí)，要選擇已更新的地質(zhì)數(shù)據(jù)庫(kù)的樣品組合文件和礦體實(shí)體模型。

6小結(jié)

介紹了利用DIMINE軟件構(gòu)建巖溶堆積型鋁土礦數(shù)字化模型的建模方法和一般流程。建模流程雖然簡(jiǎn)便，但實(shí)際操作過程比較繁瑣，特別是在構(gòu)建實(shí)體模型過程中總會(huì)出現(xiàn)一些意想不到的問題，需要不斷查找原因和進(jìn)行調(diào)整。整理地質(zhì)數(shù)據(jù)庫(kù)的基礎(chǔ)數(shù)據(jù)時(shí)要非常仔細(xì)，由于數(shù)據(jù)量非常大，很容易混淆或輸入錯(cuò)誤，要對(duì)數(shù)據(jù)表格進(jìn)行仔細(xì)校核。在對(duì)塊段模型進(jìn)行賦值時(shí)，要特別注意估值參數(shù)和約束條件不能設(shè)置和選擇錯(cuò)誤。總之，巖溶堆積型鋁土礦地質(zhì)建模是一項(xiàng)簡(jiǎn)單而又繁瑣的工作，建模過程中必須非常用心，才能快速、高效、準(zhǔn)確地將大量模型構(gòu)建完成。

參考文獻(xiàn)：

[1]祝瑞勤,奚小雙,吳塹虹,等.廣西平果堆積鋁土礦巖溶地貌演化與成礦作用[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版) ,2011,42(3):744-751.

[2]王盤喜,卞孝東,郭俊剛,等.河南巖溶堆積型鋁土礦的發(fā)現(xiàn)及其地質(zhì)特征[J].輕金屬,2013(12):1-3，52.

[3]潘冬.我國(guó)礦山數(shù)字化建設(shè)的探討[J].礦業(yè)研究與開發(fā),2006(S1):36-39.

[4]吳德政.數(shù)字化礦山現(xiàn)狀及發(fā)展展望[J].煤炭科學(xué)技術(shù),2014,42(9):17-21.

[5]蔣京名,王李管.DIMINE礦業(yè)軟件推動(dòng)我國(guó)數(shù)字化礦山發(fā)展[J].中國(guó)礦業(yè),2009,18(10):90-92.(收稿日期：2015-11-12)

作者簡(jiǎn)介：雷顯權(quán)(1985-)，男，工程師，從事地質(zhì)勘查工作，Email:xqlei_11@csu.edu.cn。