王 珂,錢建平,周永寧,任 東
(桂林理工大學地球科學學院,廣西 桂林 541000)
隨著科技的飛速發展,計算機在各個行業都得到廣泛應用,地質礦產領域亦不例外。計算機在地質領域最原始的功能是輔助地質人員更加便捷地制作各類圖件,以提高地質人員的工作效率。但是傳統的二維圖件表達各種地質體、地質現象之間的相互關系上具有一定的局限性,它不能很好地展現各種構造、各類地質體在三維空間真實形態及相互關系,而更多的是需要地質人員在腦中想象各種構造、地質體在三維中的真實形態,這樣假象出來的情況很容易和實際情況形成偏差,對實際工作的開展造成不良影響。隨著計算機技術的發展,三維圖像技術和計算機技術進一步的結合,構建三維可視化地質體模型的技術條件變得逐漸成熟,運用這一技術成果可以很好地解決上述問題。
構建三維地質體模型的技術就是指把二維地質圖和數據資料轉化為一標準格式,再在其中提取各種地質體的輪廓線,然后通過相關方法連接輪廓線使其形成形實體模型,并將其在三維空間中呈現出來,利用這種在三維空間中的地質模型分析各種地質體的大小、形態、相互關系及分布規律。這種三維地質體模型的構建為地質行業的發展提供了新的方法,為現代化礦山建設提供了新的方式。
三維地質建模(3D Geo-Modeling)主要是應用三維空間信息管理、現代計算機技術、地質解譯、空間分析與預測、地理信息技術和圖像處理及解譯技術,在三維空間模擬環境下分析和模擬地質體的技術。這個概念在上世紀九十年代由Simon.W.Houlding首先提出后,被廣泛地應用于礦產勘探、工程建設、國土資源、地理信息系統等多個領域(Simon W.Houlding,1994)[1]?,F在常見的三維地質建模多是以三維地質軟件為媒介來實現的,所以我們也可以說三維地質建模的發展過程也就是三維地質軟件的發展過程。這些軟件的發展依托于科技技術的發展,近年來在許多的地質領域的學者和科研人員共同努力下,三維地質方面的軟件有了很大的進展,例如3Dmine、Surpac和Dimine等軟件。三維地質建模是地質行業的新興課題,它特別適用于礦山地質研究,可以為礦區的生產規劃、采礦設計等服務,從而大大提高礦區的生產效率,為礦產資源的可持續開采做出一定的貢獻。
地質數據庫是現在常用的一種簡單有效的處理數據的工具,可以簡單直接的對數據進行檢索和處理。使用3Dmine建立數據庫,將收集的地質數據導入數據庫當中,通過3Dmine軟件將數字形式的勘察資料變換為三維圖形,供我們后期的應用與分析。
3DMine礦業軟件的數據庫格式是依據各種勘探工程編錄的格式而定的,其主要存儲的內容為:鉆孔空間參數數據(鉆孔孔口位置、鉆孔類型、鉆孔最大孔深等);鉆孔或坑道所揭露巖石性質特征;巖芯元素分析數據及地質剖面圖等。
3DMine礦業軟件使用的是Microsoft Access數據庫,用戶可以將鉆孔柱狀圖中的數據以文本或者Excel表格的形式導入其中。通過3DMine礦業軟件的轉換,這些文字形式的勘探資料就變成三維圖形的形態來管理及使用了。3DMine軟件的數據庫中有兩個必需的表格(定位表和測斜表)和若干非必需的表格(化驗表、巖性表等),其中定位表和測斜表是用以控制鉆孔在三維空間的位置和孔道軌跡,化驗表、巖性表等非必需性表格可以根據用戶的需求自己設定,本次研究共建立了四個表格如下:
2.1.1 定位表的作用
定位表的作用主要是錄入鉆孔等地質工程的坐標及相關的基本信息(表1)。其字段主要包含工程號、孔口坐標(X、Y、Z)、最大深度、軌跡類型及其它用戶自定義信息等,此外還有坑道、探槽等工程建立相關信息。

表1 定位表屬性結構
2.1.2 測斜表的作用
測斜表的作用主要是錄入工程號、鉆孔進尺、傾向方位和傾斜程度等相關信息(表2)。其字段包含工程號、深度、方位角、傾角等。

表2 測斜表屬性結構
2.1.3 化驗表的作用
化驗表的作用主要是錄入鉆孔或坑道所采金礦體樣的化驗分析結果等信息(表3),為組合樣生成以及之后的礦體品位賦值提供數值。字段主要含有工程號、從、至、金元素化驗結果值及其他信息等。

表3 化驗表屬性結構
2.1.4 巖性表的作用
巖性表的作用主要是錄入描述巖性及含礦程度等相關信息(表4)。字段主要含有工程號、從、至、巖性描述及巖性代碼信息等。

表4 巖性表的屬性結構
以上這四張表,將鉆孔定位、測斜、取樣化驗、地質編錄等工作分門別類,使地質勘探等實際工作記錄起來更加條理清楚。雖然這四種表的內容以及格式等是有區別的,但是這四張表中都通過“工程號”這一屬性建立了聯系,形成了一個數據網絡。建立的地質數據庫包含定位表、測斜表、化驗表、巖性表這四種表格[2]。
在地表模型構建過程中,將礦區的AutoCAD地形地質圖導入3Dmine中作為底圖,并將礦區地表出露的地層巖性、地質構造等地質信息;鉆孔、探槽、建筑等工程信息水平投影到該底圖上的地質用圖。由于圖件上的線條高程不準確,因而圖中等高線進行處理如:坐標校正、等高線連接、線上冗余點清理等,將地形地質圖上一些不需要的信息進行刪除或隱藏處理,之后利用該軟件工具欄中線賦高程功能將所有等高線賦上其對應的高程值,最后利用DTM工具生成地表模型(圖1)。
三維地質體模型就是通過不規則的三角網將礦體的各個輪廓線連接起來,形成形態各異的密不透風的網。具體步驟是先將鉆孔柱狀圖、中段地質平面圖以及勘探線剖面圖等處理好,再加載到3DMine軟件中,然后延勘探線方向將鉆孔模型切剖面,再根據要求及實際情況圈閉連接所切剖面上的礦體品位分布線或者直接在勘探線剖面圖或中段地質平面圖上提取礦體輪廓線,最后通過閉合線到閉合線的連接三角網功能生成礦體實體模型(圖2)。在實際建模過程中,因為礦體的不規則性我們總會遇到一些存在分支、復合現象的礦體模型需要處理,所以,3DMine礦業軟件在連接礦體輪廓線時可以通過利用控制線以及分區線的方法來解決礦體形態不規則的問題從而使構建地礦體形態與真實情況更接。

圖1 地表模型示例圖

圖2 礦體模型示例圖
構建地質體實體模型的方法主要有以下三種:
(1)運用三維鉆孔模型及其上的三維品位數據信息圈定礦體剖面輪廓線從而建立地質體實體模型;
(2)利用礦區勘探線剖面圖直接在圖中提取地質體豎直輪廓線從而建立地質體實體模型;
(3)利用中段地質平面圖直接提取地質體水平輪廓線從而建立地質體實體模型。在實際建模過程中,我們需要依據礦區資料的多少選擇適合的建模方法。根據上述原則,對本研究區礦體實體模型創建的方法選擇具體情況具體分析。
礦體是在各種地質作用下形成的不規則的幾何體,我們利用3DMine礦業軟件構建礦體實體模型其實只是模擬了礦體的實際輪廓,但這個輪廓也只是一個表殼,沒有實質性的內容,因此我們創建塊體模型,塊體模型其實是利用若干個一定大小的長方體來填充這個不規則的幾何實體,然后通過累加所有充填的長方體體積就可以近似得到這個實體的體積,若再給每個長方體內部添加一些屬性如比重和品位等,我們就可以計算該礦體實體模型所含有的實際礦石量和內部礦的總重量等。
搜索塊體空間:應盡可能將整個礦體實體模型都包含在里面,有時候還要包含不可剔除的夾石。
定義塊體尺寸:一般情況下,主要依據礦體的類型、規模及開采方式等來定義其尺寸大小(王李管,2008)。塊體模型是若干大小相同的長方體累加而成的,但這種等大小的塊體累加不能將不規則的礦體實體模型全部充滿,在礦體實體模型邊緣還會留存小部分的空間未被充填,因此需要使用更次一級的子塊體充填實體模型邊緣未被充填的部分,這樣才能更貼近礦體實際輪廓從而保證儲量估算的準確性。二三級模塊按幾何級數進行尺寸分割,但尺寸不能設置太小,否則容易占用電腦內存導致運行緩慢。
選擇估值方法:一般情況下,主要依據礦床類型、鉆孔及坑探取樣化驗數量及分析結果等方面綜合考慮來選擇合適的品位估值方法。對于一般普查勘探層次的礦山而言,相關數據資料較少,適合選用距離冪次反比法或最近距離法。對于詳細勘探或已進行生產的礦山而言,各方面數據資料比較齊全,通常使用克里格法。
添加約束條件:就是限制塊體模型空間形狀,約束條件有以下五種:表面、實體、塊值、閉合線和區域文件,這些條件應用在不同的地方,目的是對塊體進行選擇,它是塊體估值過程中的重要環節。
塊體模型的屬性如同數據庫的字段名一樣會被用以存儲相應的內容,而品位屬性則是對塊體模型進行估值賦值。由于鉆孔和坑探采樣相對較少,已有的品位數據不可能包含每一個充填塊體所需的品位屬性,所以,我們采用估值方法插值來賦給所有塊體品位屬性,而估值方法一般采用距離冪次反比法及普通克里格法。
根據地質資料確定礦體的體重,將此數值通過塊體模塊中單一屬性賦值法賦給已成功估值的塊體就可以計算出該礦體的總的礦石量,下一步就是通過塊體模塊中的分類功能進行計算得到該礦體的金平均品位,最后通過以下公式就可計算出該區金的金屬總量(圖5)。
P=V×D×C
公式中P代表礦體中的金含量,V代表礦體的總體積,
D代表礦石的體重,C代表礦體中金的平均品位。
通過對礦體塊體模型中每個單元塊的累加平均統計計算可得礦體體積、質量和重量,同時也證明3Dmine軟件計算與傳統方法相比更加簡單、快捷。

表5 3Dmine礦體儲量計算報告示例
3Dmine礦業工程軟件是以礦區地質、測量和采礦為基礎,在實際使用過程中應用3Dmine礦業工程軟件利用地形地質圖、勘探線剖面圖、井下中段圖及鉆孔數據等資料分別對礦山的地表模型、三維礦體模型和塊體模型進行創建,利用三者創建出礦區的三維地質體綜合體模型,建立的三維地質體綜合模型可以直觀生動的表現出礦體的特征,以便于礦山后期的開發利用,又以三維地質體模型為基礎,利用地質統計學計算方法對建立的塊體進行估值,實現了在投影方向上進行三維塊體的劃分和儲量計算,解決了傳統計算儲量方法的不穩定和經常出現不可預測的誤差等問題,同時也彌補了單純利用地學統計法計算的工作量大的問題,提高了計算儲量的效率,創新了地質工作的傳統模式,對現代化礦山建設提供了指導作用。