Surpac在超大型多金屬礦三維地質建模中的應用*

王利東，李海泉，周云滿，和中華

(1.云南黃金礦業集團股份有限公司， 云南 昆明 650224；2.SURPAC 軟件中國辦事處， 北京 100043)

1 Surpac建模和資源儲量估算

根據礦區內工程取樣控制情況，創建礦區三維地質數據庫和礦體三維實體模型，將礦體劃分成若干個小長方體塊體組成的礦體塊體模型；根據樣品品位統計分析結果，符合正態分布的礦體采用克立格法對各個小塊進行品位估值，不符合正態分布的礦體采用距離冪次反比法對各個小塊體品位進行估值；得出各小塊體的品位值，將每小塊體賦給按礦區礦石類型或回歸分析計算得到的比重值，每個長方體或立方體小塊體的體積已知，軟件可以算出各小塊的礦石量、金屬量。

由此可見，建立三維地質數據庫是基礎，創建礦體實體模型是關鍵，建立礦體塊體模型并對各個小塊體進行品位估值是目的。

1.1 距離冪次反比法基本原理

塊平均品位是用參與估值樣品中心點到待估塊質心距離的冪次的倒數為系數求得的加權平均品位。估值公式為：

不難發現，距離冪次反比法估值時，周圍工程或樣品距塊質心越遠，其對估值結果的影響就越小，因此，參與估值的工程數和各個工程中參與估值的樣品數都需作一定的限定，一般工程數為距待估塊最近的3～10個，每個工程中參與樣品數為距待估塊質心最近的3～5件，以保證待估塊品位結果的可靠性和合理性。

1.2 克立格法估值基本原理

克立格法是一種求最優、線性、無偏內插值估計量的方法。它是在考慮了信息樣品的形狀、大小及其與待估塊段相互間的空間分布位置等幾何特征以及品位的空間結構之后，為了達到線性、無偏和最小估計方差的估計，而對每一樣品值分別賦予一定的權系數，最后進行加權平均來估計塊段品位的方法。是一種加權滑動平均法，通過把礦體劃分成許多小塊段(待估塊段)，根據待估塊段周圍有限鄰域內的信息逐塊估計，組合起來就得到了整個礦體的估計量。

克立格方程組如下：

它是一個含n+1個未知數(n個λα和一個拉格朗日因子μ)的n+1個方程組。

克立格方差由下式求得：

上述公式或方程的求解程序均已固化在Surpac軟件中，可以實現快速估值。

2 北衙金礦Surpac軟件建模應用情況

北衙金礦位于滇西北鶴慶縣城172°方向、平距47.5 km處。礦區開采歷史悠久，2012年該礦床探明金儲量已達200 t，屬超大型礦床。

2011年在北衙金礦運用Surpac軟件對萬硐山礦段內各主要礦體建模，用克立格法和距離冪次反比法估算了資源儲量，提交了云南省鶴慶縣北衙鐵金礦資源儲量核實地質報告，并通過評審備案。該礦段經2011～2012年進一步加密勘探，礦體中淺部有了一定的變化，在深部發現新礦體，因此，2013年再次采用Surpac軟件，重新建立了北衙金礦萬硐山礦段地質數據庫,對各種地質要素進行了地質解譯,創建了礦體實體模型和塊體模型,計算了金元素變異函數,用距離冪次反比法和克立格法重新估算了28個主要礦體的資源儲量；2013年建模剖面地質解譯和創建礦體實體模型時，礦體外推與國內地質礦產勘查規范要求實現了扣合，使得建出的礦體實體模型符合Surpac軟件建模估算資源儲量的基本原理和中國國內固體礦產質勘查報告編寫規范要求。

2.1 國內規范礦體外推一般要求

兩相鄰工程，一個見礦，一個未見礦時，礦體的圈連稱為有限外推。據中華人民共和國國土資源部頒布的地質礦產行業標準《巖金礦地質勘查規范》(DZ/T0205-2002)，巖金礦體圈定時的有限外推一般按下述原則進行。

(1) 當礦體長度與厚度存在正相關關系并經過足夠的統計資料證實時，可以根據見礦工程控制的實際厚度，按照比例外推；

(2) 無規律可循時，一般按工程間距的1/2楔形外推或1/4平推(見圖1中a1、d3)；當邊部工程存在礦化現象(工程品位在邊界品位的1/2以上)時，則可按工程間距的2/3楔形外推或1/3平推(見圖1中a2、d2)；

(3) 見礦工程為米·百分值或米·克/噸值工程時，一般不得外推(見圖1中d4)；但對于薄脈型礦體，則可酌情外推。

當見礦工程之外沒有工程控制，或未見礦工程距離見礦工程較遠(距離大于相應地質可靠程度對應網度)時，礦體的圈連稱為無限外推。無限外推時，若礦體長度與厚度之間無規律可循，一般按相應地質可靠程度所對應網度的1/2楔形外推或1/4平推(見圖1中b1、b2、d1)。

相鄰兩工程，一個見工業礦，一個見低品位礦時，采用對角尖滅(見圖1中c1、c2)。

2.2 以往北衙金礦Surpac軟件建模礦體外推

據云南省鶴慶縣北衙鐵金礦資源儲量核實地質報告[3]，北衙金礦萬硐山礦段以往詳查，主要礦體工程控制網度已達80 m×80 m，局部已達40 m×40 m。

2.2.1 以往剖面解譯時礦體外推

在剖面礦體解譯方面，以往外推原則方法與2013年建模采用的是相同的。

相鄰兩個工程，一個見礦，另一個不見礦，先在未見礦工程中找到賦含所連礦體的礦化層或標志的中點作為對應點，然后從見礦工程向此對應點外推兩工程間距的二分之一得到礦體外推邊界點(見圖2中上、中下部)。

圖1 國內礦體外推示意

相鄰兩工程間距大于(333)類資源量網度要求的160 m時，或剖面邊部見礦工程外再無工程時，采用無限外推。先過邊部見礦工程和與其緊鄰的內側見礦工程切一輔助剖面，順著兩工程礦體的外延趨勢，從邊部見礦工程向外楔形外推控制的內蘊經濟資源量(332)工程間距的二分之一(40 m)得到礦體外推邊界點，如圖2中上部所示，或向外平推控制的內蘊經濟資源量(332)工程間距的四分之一(20 m)得到礦體外推邊界線，如圖2左側所示。

2.2.2 以往創建實體模型時礦體沿走向外推

(1) 邊緣剖面只有一個見礦工程時，沿相鄰剖面到邊緣剖面礦體走向變化趨勢，從邊緣剖面上見礦工程向外錐形外推控制的內蘊經濟資源量(332)工程間距的二分之一(40 m)得到外推點。

通識教育課程（必修部分）的比例為20%，學科專業知識課程的比例為35%，教師教育類的課程占到了45%的比重。其中，理論和實踐類課程各占62%和38%的比例。教師教育課程包括教學職業素養課、教育基本理論課、教學技術能力課和教師實踐能力等。教學計劃用兩年時間完成通識教育課程、學科專業課程、教師教育課程（教學職業素養課+教育基本理論課+教學技術能力課等），第五學期實施教師實踐能力課程?；蛘呱钊胄W課堂一線，進行教育實習（252學時），強化實踐教學技能，第六個學期進行畢業實習與創新創業實踐（252學時）。

(2) 邊緣剖面多于一個見礦工程時，用該剖面上礦體中心線或礦體外邊線，沿相鄰剖面到邊緣剖面礦體走向變化趨勢，從邊緣剖面上向外楔形外推控制的內蘊經濟資源量(332)工程間距的二分之一(40 m)得到外推邊界線或向外板推控制的內蘊經濟資源量(332)工程間距的四分之一(20 m)得到外推邊界面，如圖3(a)所示，外推部分的長度為邊緣剖面上見礦工程間的長度及剖面兩端外推的長度。

(3) 中間兩剖面間距大于(333)類資源量網度要求的160 m時，將該兩條剖面均視為邊緣剖面，按本小節①、②條分情況無限外推。

(4) 中間剖面間礦體邊界線：將相鄰各勘探線剖面解譯時得到的邊界線或點直接相連得到，建出的KT53礦體實體模型見圖3(a)。

(5) 主剖面間的加密剖面上沒有新增工程控制時，建實體前根據相鄰主剖面的礦體的解譯線進行加密解譯，即相鄰主剖面上對應的見礦部位相連，與加密的勘探線求交點，用這些交點相連得到加密線上外推的圖切礦體線，如圖4(c)中的藍色線。

(6) 礦體內“天窗”的扣除按剖面解譯的有限外推方法處理。

圖2剖面解譯礦體有限外推和無限外推截圖(據鶴慶北衙礦業有限公司，2013)

2.3 本次北衙金礦Surpac軟件建模礦體外推

據云南省鶴慶縣北衙鐵金礦資源儲量核實地質報告，經2011～2012年加密勘探，北衙金礦萬硐山礦段露采區礦體工程控制網度大部已達20～40 m×40 m，露采區外大部已達80 m×80 m。

2.3.1 本次剖面解譯時礦體外推

相鄰兩工程間距大于本次建模(333)類資源量網度要求的160 m時，或剖面邊部見礦工程外再無工程時，采用無限外推。先過邊部見礦工程和與其緊鄰的內側見礦工程切一輔助剖面，順著兩工程礦體的外延趨勢，從邊部見礦工程向外楔形外推本次控制的內蘊經濟資源量(332)工程間距的二分之一(40 m)或向外平推本次控制的內蘊經濟資源量(332)工程間距的四分之一(20 m)得到礦體外推邊界點(線)。如圖2左側所示。

2.3.2 本次創建實體模型時礦體外推

礦體走向外推遵循剖面解譯時礦體傾斜方向上的外推原則：

(1) 邊緣剖面只有一個見礦工程時，走向外推采用無限外推。沿相鄰剖面到邊緣剖面礦體走向變化趨勢，從邊緣剖面上見礦工程向外錐形外推本次控制的內蘊經濟資源量(332)工程間距的二分之一(40 m)得到外推點,見圖4(a)。

(2) 邊緣剖面多于一個見礦工程時，用該剖面上礦體中心線或礦體外邊線，沿相鄰剖面到邊緣剖面礦體走向變化趨勢，從邊緣剖面向外楔形外推控制的內蘊經濟資源量(332)工程間距的二分之一(40 m)得到外推邊界線，或向外板推控制的內蘊經濟資源量(332)工程間距的四分之一(20 m)得到外推邊界面，外推邊界線(面)的長度對應于邊緣剖面上見礦工程之間的長度，見圖3(b)中54線、82線處所示和圖4(b)所示。

(3) 中間兩剖面間距大于本次建模(333)類資源量網度要求的160 m時，將該兩條剖面均視為邊緣剖面，按本小節①、②條分情況無限外推,見圖4(c)。

(4) 中間剖面間的有限外推：先找到相鄰剖面上需向之外推的未見礦工程中賦含本礦體的礦化層或標志層，礦化層或標志層的中點即為有限外推的對應點，再從見礦工程向對應點楔形外推本次控制的內蘊經濟資源量(332)工程間距的二分之一(40 m)得到外推邊界點或向對應點平推本次控制的內蘊經濟資源量(332)工程間距的四分之一(20 m)得到外推邊界點。

這里還有一個前提條件，即見礦工程到對應點的距離不大于80 m，大于80 m時，則只楔形外推40 m或平推20 m。選擇向對應點楔形外推或是平推取決于鄰近剖面上的情況，如鄰近剖面上采用尖推，則也選用尖推，反之則選用平推，以保持礦體實體在這一段總體協調為準。

圖3北衙金礦KT53礦體外推對比圖(據鶴慶北衙礦業有限公司，2013)

(5) 主勘探線剖面間的加密剖面上沒有新增工程控制時，建實體前，根據相鄰主剖面的礦體的解譯線進行加密解譯：相鄰主剖面的對應的見礦部位相連，與加密的勘探線求交點，用這些交點相連得到加密線上外推的圖切礦體線，見圖4(c)中的藍色線。

(6) 礦體內“天窗”的扣除按中間剖面間的有限外推方法處理,見圖4d。

按上述原則建出的北衙金礦KT53礦體實體模型見圖3(b)。

3 Surpac軟件建模估算資源儲量注意事項

北衙超大型金多金屬礦，主要礦體多達28個，礦體形態較為復雜，需估算資源量的元素有金、鐵(全鐵和磁性鐵)、銀、銅、硫、鉛、鋅等8種，且要分殘坡積礦、氧化礦、原生礦和不同礦石自然類型、證內證外、保有消耗、資源儲量類別等統計報告資源儲量，工作量極大，如用傳統方法完成此項工作，將耗費大量人力和時間。采用Surpac軟件建模進行資源儲量估算，利用該軟件強大的宏命令錄制編輯功能，將大量的數據計算、分析統計工作和資源儲量報告中一些重復性工作，實現了程序化、快速化完成的效果，三個月時間內就完成了地質建模、塊模型品位估值和資源儲量報告等工作。充分體現了該軟件在大型以上復雜礦區或礦山資源評價中的高效率優勢。

但需要注意的是，Surpac軟件建模估算的資源儲量，是礦體模型之內的部分，模型之外想估算資源儲量，不管把橢球體的搜索半徑放大到多少都是沒用的；另一方面，如果沒有嚴格遵循地質規律或國內地質規范，創建的礦體實體模型局部可能會超出礦體實際的空間形態，估算出的資源儲量就會偏大。因此，建模前就得把資源儲量需估算的范圍內的礦體模型都建全、建到位，并按礦區地質規律、國內地質勘查規范要求進行建模。

通過前面對北衙金礦前后兩次建模礦體外推方法原則的對比，不難發現2013年建立的實體模型更為合理和符合地質規律。

圖4本次建模礦體外推截圖(據鶴慶北衙礦業有限公司，2013)

以KT53為例，采用2011年外推方法與采用2013年外推方法建出的礦體實體模型見圖3(a)和圖3(b)。對比前面所述的礦體外推方法可見，以往使用Surpac軟件建模時，礦體走向外推沒遵循或沒完全遵循國內行業規范要求，而是將各相鄰剖面解譯時的邊界點(線)直接相連作為礦體實體模型的邊界。這樣處理，在實際勘探工程按工程網度近理想分布且都見礦時，影響和差別不大；但當實際勘探工程局部偏離勘探線或部分工程未見礦，或加密勘探線上只有少量或個別工程見礦時，建出的礦體實體模型就顯得不太合理或不符合地質規范要求，如圖3(a)中藍色線框區礦體實體是多建的部分，顯然會多算出一部分資源儲量，而紅色線框區則是少建的部分，會少算部分資源儲量。經估算，2013年外推法建模估算出的資源儲量較2011年外推法建模估算出的資源儲量見表1，對比數據顯示，金金屬總量減少了822千克，約26.66%，(332)類以上資源儲量增加了1031千克，約35.42%，(333)類資源量增加了209千克，約62.03%。影響較明顯。

表1 KT53礦體金資源儲量估算結果對比

注：A. 2013年外推法建模; B. 2011年外推法建模

4 結 論

Surpac軟件建模估算的資源儲量，礦體實體的建立非常重要，實體越接近礦體的實際空間形態，估算出來的資源儲量就越準確，而建好實體的關鍵是礦體外推方法原則。2013年北衙金礦萬硐山礦段Surpac軟件建模，針對剖面地質解譯和創建礦體實體模型時，礦體外推與國內地質礦產勘查規范要求相結合這一問題進行了改進，使得建出的礦體實體模型既符合Surpac軟件建模估算資源儲量的基本原理，又更趨于合理和符合地質規律，估算出來的資源儲量也較為準確。

參考文獻：

[1]候德義.找礦勘探地質學[M]. 北京：地質出版社，1984.

[2]鶴慶北衙礦業有限公司.云南省鶴慶縣北衙鐵金礦資源儲量核實地質報告[R].鶴慶縣:鶴慶北衙礦業有限公司,2011,2013.

[3]和中華，周云滿，和文訁，等.滇西北衙超大型金多金屬礦床成因類型及成礦規律[J].礦床地質，2013，32(2)：244-258.

[4]楊建宇，秦德先，康澤寧，等. 北衙金礦三維模型的建立及研究[J].有色金屬，2006，58(2)：17-20

[5]DZ/T0205-2002. 巖金礦地質勘查規范[S].

[6]張樹泉.地質統計學講義[M].北京：北京科技大學，2010:1-38.