Surpac 在金屬礦山三維建模中的應用研究

任志遠 靳 洋 熊章華 熊開治

(1.重慶市地勘局208地質隊，重慶 400700;2.重慶一三六地質隊，重慶 401147;3.重慶市江津區國土局，重慶 402200)

1 礦床地質概況

本次研究的鋅銅礦體屬似層狀矽卡巖型鋅銅礦床，共有7個主要礦體，平均控制長685 m，寬831 m～1 747 m，平均厚度5.56 m。礦床中最主要的有用金屬元素為Zn，占探獲資源總量的49%，平均品位為4.23%，其次為 Cu，Ag，平均品位分別為0.17%，320×10－6，伴生礦產有硫、砷、銦、鎘、鎵等，含量達到綜合利用的要求。

2 三維地質數據庫的創建

2.1 原始資料的處理

本次收集了礦山勘探工程數據，按照Surpac軟件格式要求對鉆孔數據進行處理:1)刪掉無品位信息的鉆孔;2)刪掉礦體外的鉆孔資料;3)將所有鉆孔名的第一個字符設置為大寫;4)將向下鉆孔的傾角設置為負值;5)修正因人為因素造成的錯誤數據;6)保留本次研究的鋅(Zn)、銅(Cu)、銀(Ag)三種元素的品位數據值，其他的刪除;7)因缺少圍巖巖性資料，最后將數據表分成3個表，分別為鉆孔表(Collar)、測斜表(Survey)和品位分析表(Sample)。通過處理后，得到44個鉆孔數據資料，739組測斜數據和2 473個樣品化驗數據，并保存為.csv或.txt格式。

2.2 鉆孔數據庫的創建

通過Surpac導入Zn，Cu，Ag三種元素數據源，創建地質數據庫.ddb，建立了鉆孔表(Collar)、測斜表(Survey)和品位分析表(Sample)，各表儲存的信息如圖1所示。

圖1 鉆孔數據庫

2.3 數據庫后期處理

鉆孔數據庫是礦床三維數字建模的基礎，后續實體模型、品位模型的創建及統計分析、儲量計算等都離不開該數據庫。另外，該數據庫還具有數據編輯、審核、更新、查詢、統計分析以及三維可視化顯示等操作。

3 礦體數學模型的構建及研究

3.1 樣品的統計分析

本次采用普通克立格法進行樣品的統計分析，根據勘探線產生組合點，取樣長為1 m，組合樣品點后，利用Surpac軟件的基礎統計窗口對樣品品位的離散值進行平均值、方差值等統計，生成Zn，Cu，Ag三種元素的品位統計分析直方圖，并根據統計報告得到該礦區Zn，Cu，Ag三種元素的統計分析結果(見表1)。

表1 樣品統計分析結果

3.2 變異函數統計分析及驗證

1)特高品位的處理方法。由表1得Zn，Cu，Ag變化系數均大于150%，屬不均勻礦床，應對特高品位處理，一般特高品位的最低界限值是其平均品位的15倍，故本次Zn品位高于23.43(156.2% ×15)的品位為特高品位，以 23.43代替;Cu品位高于26.13(174.2% ×15)的樣品品位視為特高品位，以26.13代替;Ag品位高于44.67(297.8% ×15)的樣品品位視為特高品位，以44.67代替。由鉆孔數據可知，Cu品位最大值為3.8，沒有特高品位，而 Zn品位最大值為24.366，Ag品位最大值為1 044.74，均屬于特高品位，應對Zn，Ag進行特高品位處理。通過Surpac軟件編輯菜單設置運算字段后，得到處理后的Zn，Ag品位統計分析直方圖并生產統計分析結果(見表2)。

表2 特高品位處理后的統計分析結果

對比表1和表2可知，由于Zn品位中大于特高品位的樣品很少，經過處理后統計直方圖變化不大，而Ag的變化比較大，雖然統計后的均值、方差、標準離差、偏度等變化較大，但表2較為合理。

2)變異函數的計算和統計分析。變異函數計算是品位估值計算的基礎，一般沿礦體走向、傾向和垂向進行實驗變異函數計算，通過變異函數建模窗口設置參數并輸入礦體走向、傾向和垂向基礎數據，運行后得到Zn，Cu，Ag的變異函數散點，經擬合后可得理論變異函數曲線，這些參數將為克立格賦值提供依據。

3)變異函數驗證。變異函數驗證選擇方差驗證法，判斷條件如下:a.計算的殘差均值要基本接近于0;b.實際方差與理論方差間的誤差要求小于15%;c.計算的兩標準的殘差參數要基本接近95%;d.繪出的兩類直方圖要基本服從正態分布。進行驗證后的結果均服從這4個條件時，說明擬合的數據是比較正確的。

運行“地質統計—方差驗證”，輸入理論變異參數值，執行后獲得殘差直方圖及驗證統計分析報告(見表3)，由統計結果可知Zn，Ag方差驗證基本符合條件，說明Zn，Ag變異函數曲線擬合是正確的，而Cu的方差驗證有一定的誤差，這與Cu品位的分布特性有關，可多次擬合其變異曲線，直到基本符合要求。

表3 克立格殘差驗證結果統計表

4 礦體實體模型的建立與研究

礦體實體模型是Surpac三維模型的基礎，它的構建一般有三種方法:1)勘探線法;2)礦體的邊界線法;3)鉆孔數據法，一般綜合采用以上幾種方法，本次綜合采用了方法1)，3)，即通過鉆孔數據建立鉆孔三維模型，按相鄰勘探線間礦體的展布趨勢，大致確定礦體的范圍并連三角網，然后封閉后形成礦體實體，步驟如下:

1)利用建好的三維地質數據庫，在2D網格顯示圖上，執行“數據庫—剖面—定義剖面”操作，定義面前、面后及步距，大致沿著勘探線畫一條剖面線，完成定義剖面;2)在剖面方向上，根據礦體品位的數值(用不同顏色表示)及礦體的趨勢用線段連接大致圈定礦體，并保存為.str文件;3)執行“實體模型—創建三角網”連接段，最后封閉形成礦體實體模型(見圖2)。礦體三維實體模型建立后利用Surpac軟件可以進行以下操作:1)礦體三維可視化顯示;2)生成礦體表面積和體積報告;3)任意方向切制剖面;4)在模型上疊加地質數據庫數據后，可創建塊模型、品位模型;5)將鉆孔、坑道等的數據疊加后，可分析礦山工程布置情況。

圖2 礦體三維實體顯示

5 礦體品位塊體模型的創建與研究

為了表達礦體的品位分布情況，需建立品位塊體模型，模型的范圍覆蓋了礦山所有鉆孔控制的區域，結合礦山現有采礦方法、勘探網度及變異函數的分布特征等因素進行定義，確定參數后通過添加模型屬性值，選擇已建立的實體模型為定義約束對象，執行后調入塊體模型并顯示，同時加載定義的約束文件，即可顯示該礦體塊體模型(見圖3)。

圖3 礦體塊體模型

礦床品位模型建立后，利用Surpac軟件塊體模型部分的距離冪次反比法進行塊體估值計算，包括Zn，Cu，Ag的品位、金屬量以及塊體積和礦石量，還可用不同顏色、圖案表示各金屬元素品位的分布情況，也能生成任意方向的品位分布剖面圖，揭示礦體內部品位的分布變化情況。

6 礦體儲量的計算

根據創建的品位塊體模型，計算各金屬儲量，并生成各元素的品位—噸位曲線分別如圖4～圖6所示。

圖4 Zn品位—噸位圖

圖5 Cu品位—噸位圖

圖6 Ag品位—噸位圖

由圖4～圖6可知，礦體平均品位隨礦體邊界品位的增大而增大，而金屬量隨礦體邊界品位的增大而下降。由于礦床相關工業指標的任何變化(如邊界品位)都會導致礦體形態、平均品位、礦體儲量、礦山生產能力及服務年限等的變化，最終可導致礦山生產經濟效益等的變化，所以對礦床建模后可以選擇不同的方案來圈定礦體，并根據生產成本、生產規模、市場價格及投資費用等，綜合考慮資源的開發利用，最終選出適合礦山具體實情的礦體圈定方案，以便礦山企業能夠獲得更多的利潤。

7 結語

本次通過三維建模軟件Surpac創建了某礦山地質數據庫，建立了礦體三維實體模型、礦體品位塊體模型等，并采用普通克立格法對品位進行估值，根據估值結果計算礦體儲量。本次應用研究表明，利用Surpac軟件所建模型可靠，計算結果準確，是礦山企業進行礦床建模較好的開發工具，適合于在地質、測量、采礦等相關行業用于輔助礦山進行資源評估、采礦設計以及計劃編制等工作，能夠動態管理礦山的實時生產情況，具有廣泛的應用前景。

［1］Surpac Software International國際軟件公司.Surpac Vision軟件用戶使用手冊(第5版)［Z］.2003.

［2］張仁鐸.空間變異理論及應用［M］.北京:北京科技出版社，1993.

［3］陳建宏，古德生.礦業經濟學［M］.長沙:中南大學出版社，2007.