張 敏 楊 陽 鄧中飛 陳邦學(xué) 朱慧秦

（新疆維吾爾自治區(qū)地質(zhì)礦產(chǎn)勘查開發(fā)局第十一地質(zhì)大隊 昌吉 831100）

當(dāng)前在地質(zhì)勘查、礦山開發(fā)及地質(zhì)研究等方面加快礦業(yè)信息化和三維可視化是業(yè)界已達成的普遍共識，以AutoCAD、Mapgis等二維軟件為主的情況也正在逐漸變化，三維軟件逐漸在地礦行業(yè)普及應(yīng)用。同時以數(shù)學(xué)統(tǒng)計為核心的地質(zhì)統(tǒng)計學(xué)儲量計算方法，將逐漸發(fā)展成為主流的儲量計算方法，以手工計算為主的傳統(tǒng)方法逐漸被地質(zhì)統(tǒng)計方法所替代[1]。礦區(qū)儲量計算也從二維向三維逐漸發(fā)展，將礦床的三維建模與儲量計算結(jié)合起來是當(dāng)前礦產(chǎn)資源儲量計算新的手段，它能夠更加直觀、準(zhǔn)確的對礦產(chǎn)儲量進行計算。

1 建立礦區(qū)數(shù)據(jù)庫及地形模型

數(shù)據(jù)庫的建立是為三維建模提供數(shù)據(jù)基礎(chǔ)，其中記錄了大量原始地質(zhì)數(shù)據(jù)信息。以鉆孔數(shù)據(jù)庫為例，其主要由工程定位、樣品化驗數(shù)據(jù)、工程編錄數(shù)據(jù)等文件組成（表1）。數(shù)據(jù)庫的構(gòu)建同時實現(xiàn)了對礦山地質(zhì)資料數(shù)據(jù)的動態(tài)查詢和管理，能夠三維顯示地質(zhì)信息、進行基本統(tǒng)計分析、地質(zhì)解譯等[2]。

表1 鉆孔數(shù)據(jù)庫各文件示例表

數(shù)字地形模型（DTM）就是地形表面形態(tài)屬性信息的數(shù)字表達，是帶有空間位置特征和地形屬性特征的數(shù)字描述。軟件中采用TIN技術(shù)構(gòu)建地表模型，TIN模型可以根據(jù)地形的具體特點和復(fù)雜程度而確定采樣點的密度和位置，避免數(shù)據(jù)的冗余，并相對準(zhǔn)確的描述出地形的結(jié)構(gòu)和一些局部地區(qū)，同時又便于進行地形計算和分析，能夠較好地反映實際地形信息[3]。將矢量化的地形圖導(dǎo)入到軟件中進行高程賦值后，生成數(shù)字地形模型。

2 三維礦床建模

三維礦床建模是為了將礦體、斷層等在三維空間可視化，同時也是對數(shù)據(jù)的有效集成和信息挖掘的最佳方法。將技術(shù)人員對所研究區(qū)域各地質(zhì)體的認(rèn)識有效的體現(xiàn)在所建立的模型中，使所建立的模型最接近于實際地質(zhì)環(huán)境，從而為后期礦塊品位插值、儲量估算等奠定基礎(chǔ)。3DMine軟件提供了多種方法進行三維地質(zhì)建模，如自動建模、基于原始勘探線剖面的建模、基于鉆孔數(shù)據(jù)庫的建模等。

2.1 基于勘探線剖面的建模

在勘探線剖面圖的基礎(chǔ)上，整理提取出礦體邊界線、礦體編號、坐標(biāo)網(wǎng)、地表界線等有用信息，進行優(yōu)化、糾錯等處理，去除不利于建模的無關(guān)信息，經(jīng)過多次坐標(biāo)變換將二維剖面數(shù)據(jù)轉(zhuǎn)換到三維空間中正確的位置，使各剖面圖的位置與勘探錢相互對應(yīng)（圖1）。軟件中設(shè)置合適的三角網(wǎng)參數(shù)，在相鄰的勘探線剖面圖中選擇同一礦體的閉合線，依次連接成三角網(wǎng)。在實際的建模過程中會有礦體的分支復(fù)合較多、對應(yīng)較難的現(xiàn)象，此時在三維礦體的連接、外推及分支復(fù)合的處理上，可創(chuàng)建多個“輔助線”和“分區(qū)線”對礦體形態(tài)進行控制，最終形成完整的礦體模型。

圖1 勘探線剖面建模過程示意圖

2.2 基于鉆孔數(shù)據(jù)庫的人機交互建模

在鉆孔數(shù)據(jù)庫的基礎(chǔ)上，經(jīng)過數(shù)據(jù)庫檢驗、特高品位處理、組合樣長、鉆孔設(shè)置及創(chuàng)建勘探線等步驟后，將一定容差范圍內(nèi)的鉆孔投影到切割的剖面上，根據(jù)邊界品位、鉆孔數(shù)據(jù)信息及地質(zhì)規(guī)律和建模原則對不同剖面的礦體進行解譯。解譯后根據(jù)相鄰剖面礦體賦存位置、地質(zhì)規(guī)律、外推和尖滅原則等建立礦體面模型，最后為保證模型的可靠性要進行實體驗證，對錯誤實體進行修正。礦體模型的創(chuàng)建流程復(fù)雜，在圈連過程中容易出現(xiàn)邏輯錯誤，如自相交、無效邊、開放邊等，因此必須對建立的模型進行驗證，驗證合格的實體模型才可以進行查詢、編輯、體積計算等操作。建立的礦體模型具備三維可視化功能，并可在其內(nèi)部任意設(shè)置剖面瀏覽，全面掌握礦體的形態(tài)、產(chǎn)狀及空間分布規(guī)律等（圖2），技術(shù)人員能夠更方便直觀的研究礦床的品位特征、分布規(guī)律及礦化與地質(zhì)環(huán)境的關(guān)系等。

圖2 切割狀態(tài)下三維礦床模型效果圖

2.3 自動建模

3Dmine提供一種在鉆孔數(shù)據(jù)庫基礎(chǔ)上的自動建模方法，通過設(shè)置邊界條件、網(wǎng)格、剔除厚度等參數(shù)，直接生成礦體模型。該自動建模的礦體模型是基于固定的算法生成，針對不同的礦床其礦體形態(tài)可能與專業(yè)技術(shù)人員的預(yù)期有一定的差異，且不同礦床的礦體形態(tài)與礦床成因和地質(zhì)環(huán)境密切相關(guān)，故因綜合考慮判斷自動建模模型的可靠性或僅做參考使用。另外，利用地質(zhì)統(tǒng)計學(xué)對礦體進行結(jié)構(gòu)各向異性分析，分析得出主軸/次軸、主軸/短軸、主軸方位角、主軸傾伏角、次軸傾角等參數(shù)值。將各向異性分析結(jié)果加入到自動建模過程中生成礦體模型（圖3）。通過自動建模結(jié)果的比較，認(rèn)為加入礦體結(jié)構(gòu)各向異性參數(shù)為指導(dǎo)的自動建模較直接自動建模的結(jié)果可靠性更強，更加接近礦體的實際真實形態(tài)。

圖3 結(jié)構(gòu)各向異性指導(dǎo)的自動建模效果圖

3 儲量估算及方法比較

3.1 地質(zhì)塊段法儲量計算

軟件中用地質(zhì)塊段法進行儲量估算試驗（圖4），完成礦段圈定、礦體投影、創(chuàng)建塊段及儲量報告等，得出軟件與人工儲量計算的相對誤差為4.6%，相對誤差較小，該方法本身不存在誤差，推測軟件與人工計算引起誤差的原因是人為因素或投影產(chǎn)生的誤差等。

圖4 地質(zhì)塊段法儲量計算三維塊段模

利用3Dmine軟件進行傳統(tǒng)的地質(zhì)塊段法儲量計算較人工方法計算更加方便快捷，在數(shù)據(jù)量較大時更顯優(yōu)勢，同時能夠減少人為因素產(chǎn)生的誤差，提高技術(shù)人員的工作效率，方便后期的數(shù)據(jù)管理及應(yīng)用等。

3.2 地質(zhì)統(tǒng)計學(xué)儲量估算

地質(zhì)統(tǒng)計學(xué)是以區(qū)域化變量理論為基礎(chǔ)，以變異函數(shù)為主要工具，研究那些在空間分布上既有隨機性，又有結(jié)構(gòu)性的自然現(xiàn)象的科學(xué)?；诘刭|(zhì)統(tǒng)計學(xué)的儲量估算方法，考慮了礦石的空間品位變化，可提高資源儲量估算的準(zhǔn)確性[4]。

實驗半變差函數(shù)分析。在地質(zhì)統(tǒng)計學(xué)中，變差函數(shù)是最基本與最重要的模擬工具，是計算區(qū)域化變量空間結(jié)構(gòu)特征的函數(shù)[5]。通過計算變量得到實驗變差函數(shù)，擬合出理論變差函數(shù)，再將不同方向的擬合結(jié)果套合為一個函數(shù)，用于估值計算。根據(jù)礦體產(chǎn)狀計算三個方向的半變差函數(shù)，走向、傾向和厚度方向分別對應(yīng)主軸、次軸和短軸，最終選擇擬合理論變差函數(shù)的參數(shù)：塊金值、基臺值及變程。本次分析使用指數(shù)模型。

克里格估值。估值前需建立礦體空塊模型和交叉驗證，對礦體模型內(nèi)部進行塊體分割，形成由n個體元模擬的區(qū)域，每個體元都有相應(yīng)坐標(biāo)并可以任意添加屬性，克里格法估值結(jié)果會作為一個屬性值加入塊模型。交叉驗證是對所擬合的理論半變差函數(shù)的可靠性進行檢驗，即從已知信息樣本總體中抽離某個點，用抽離點以外的其他信息結(jié)合半變差函數(shù)模型估計被抽離點，結(jié)果得到一系列估計值和誤差等參數(shù)，根據(jù)誤差判斷模型的準(zhǔn)確性。軟件中設(shè)置估值參數(shù)、搜索參數(shù)、分區(qū)參數(shù)、孔約束等，用普通克里格法進行空間估值，在估值過程中需通過改變主軸搜索半徑、分區(qū)設(shè)置、約束條件等參數(shù)進行了多次克里格估值，確保所有的空塊均被賦予品位等屬性。

資源儲量估算結(jié)果對比。對普通克里格法估值的塊體模型進行資源量報告，并與人工地質(zhì)塊段法儲量計算結(jié)果進行比較，得出人工地質(zhì)塊段法與普通克里格法兩種方法L1、L2礦段估算結(jié)果（表2）的相對誤差分別為0.73%、4.40%，相對誤差均小于5%，說明儲量估算結(jié)果的可靠性較強，同時提高了儲量估算的準(zhǔn)確性。不同的估算方法產(chǎn)生誤差是正常的，誤差產(chǎn)生的原因是多方面、綜合性的，但兩種估算方法的相對誤差較小，估算結(jié)果可靠。

表2 資源量估算結(jié)果對比表

4 結(jié)論

⑴通過三維礦床建模，改進了對礦區(qū)地質(zhì)數(shù)據(jù)的認(rèn)識和理解，有利于提高信息的利用率和空間分析能力，為地質(zhì)人員在三維空間中觀察、分析研究地質(zhì)現(xiàn)象以及空間展布特征等提供了一種全新的方法。

⑵根據(jù)國內(nèi)廣泛采用幾何方法進行儲量估算、提交報告的實際情況，以地質(zhì)塊段法為例進行軟件和人工儲量計算的對比，結(jié)果顯示兩種方法的相對誤差小，且運用三維軟件計算儲量避免了人工計算繁瑣、人為因素多的缺陷，大大提高了技術(shù)人員的工作效率，同時軟件可運用多種方法驗證對比，提高了儲量計算結(jié)果的可靠性，并在實際生產(chǎn)中快速給出儲量報告。

⑶基于3DMine軟件實現(xiàn)了礦區(qū)地質(zhì)建模及儲量估算的三維可視化，給出了基于三維可視化的礦床建模和儲量估算的流程方法，并運用多種方法快速有效地實現(xiàn)資源儲量可視化評價，促進礦區(qū)地勘工作由二維向二三維一體化轉(zhuǎn)變。