張 偉,黃美成,岳鵬軍
(1.北京睿思礦科科技有限公司,北京 100080;2.且末縣邦泰礦業(yè)投資有限公司,新疆 巴音郭楞 841900;3.中色地科藍天礦產(chǎn)(北京)有限公司,北京 100020)
在礦床研究和礦產(chǎn)資源量估算過程中,往往由于多期次的成礦作用和復(fù)雜多樣的控礦因素,造成成礦部位地質(zhì)構(gòu)造比較復(fù)雜,礦體形態(tài)也不規(guī)則,使原始數(shù)據(jù)的解譯工作變得十分困難,對礦床形態(tài)的理解也不能進行直觀的展示和分析。
針對上述問題,國際上多個礦業(yè)大國在20世紀中葉,開始著手編制了一系列的三維礦業(yè)軟件,如Surpac、Datamine、Micromine等大型專用三維礦業(yè)軟件,它們分別具有不同的優(yōu)缺點,其中Micromine軟件在完成數(shù)據(jù)采集和處理、三維建模、儲量估算等方面具快速和易操作性,因此其對地勘單位來說具有很強的實用性[1]。利用三維可視化,可以對礦床地質(zhì)特征、礦體空間分布特征進行直觀分析,同時實現(xiàn)相關(guān)數(shù)據(jù)的一體化管理,對提高礦產(chǎn)資源勘查過程中所需的分析判斷能力和找礦預(yù)測都有很大幫助。
資源量估算的基礎(chǔ)數(shù)據(jù)主要來源于鉆孔、探槽等原始數(shù)據(jù);將數(shù)據(jù)整合在一起,實現(xiàn)一體化管理便是數(shù)據(jù)庫的主要作用。關(guān)于鉆孔數(shù)據(jù),Micromine所使用的是自研的關(guān)系型數(shù)據(jù)庫(dhdb.格式),該數(shù)據(jù)庫僅作為定位表、測斜表、樣品分析表、巖性表等鉆孔數(shù)據(jù)的橋梁,本身并不具備儲存數(shù)據(jù)的功能,所以要求數(shù)據(jù)源是與之匹配的文件格式(DAT.格式)。這里我們可以借助Micromine的數(shù)據(jù)導(dǎo)入功能將EXCEL、CSV等常用數(shù)據(jù)格式轉(zhuǎn)化為所需的DAT格式,亦可從其他數(shù)據(jù)庫中導(dǎo)入[2]。
建立鉆孔數(shù)據(jù)庫前,需通過校驗來檢驗鉆孔數(shù)據(jù),通常原始數(shù)據(jù)會由于記錄、存放等各種原因造成數(shù)據(jù)的錯誤或缺失,鉆孔校驗功能可以很好的驗證鉆孔數(shù)據(jù)的合理性與準確性。
關(guān)于鉆孔數(shù)據(jù)的三維顯示,主要可分為顯示鉆孔軌跡、值、巖性三部分,這三部分基本滿足后續(xù)我們對礦體進行解譯和估算的需求。靈活使用三者的顯示設(shè)置,在實現(xiàn)鉆孔個性化顯示的同時,也可以達到不同的地質(zhì)目的,如顯示鉆孔軌跡中的鉆穿點功能,為繪制各類投影圖、切面圖提供了一個簡單、快捷的方式,此功能在探礦各個階段均可節(jié)省大量的時間。同時在Micromine 2022版中,還有多種其他設(shè)置可以使鉆孔顯示的內(nèi)容變得更加豐富、全面。
通過鉆孔的個性化顯示,結(jié)合剖面控制文件可以實現(xiàn)多個簡潔、直觀的礦體連接界面。剖面控制文件是通過勘探線或手工切割剖面來創(chuàng)建的、用于儲存一系列相應(yīng)空間位置的文件,其作用就是使計算機顯示特定位置和范圍內(nèi)的相關(guān)數(shù)據(jù)。
在地質(zhì)解譯、礦體圈連過程中,我們可以使用比例點功能與復(fù)制/移動線功能,準確的繪制1/2、1/3、2/3等滿足地質(zhì)勘查規(guī)范的輔助線和礦體界限并進行礦體外推。
構(gòu)建線框,首先需確定所要構(gòu)建線框的類型,如礦體、圍巖、DTM表面等;在建立工程文件夾時,Micromine針對三角片,自動創(chuàng)建有多個tridb.格式的三角片數(shù)據(jù)庫,用于分類管理線框模型,同時該類數(shù)據(jù)庫具有存儲數(shù)據(jù)的功能。關(guān)于線框的分支復(fù)合,Micromine采用的方法是,通過對構(gòu)建線框所用到的礦體界線進行編輯來實現(xiàn)的,將需分支、合并的礦體界線用分離即分割線功能將其分割為兩條或多條多段線,分別表示礦體分支或復(fù)合的位置,以達到針對不同分支分別構(gòu)建線框的目的[3]。在確定礦體延伸方位、傾角、距離等參數(shù)后,亦可采用對工程控制邊界處的礦體界線進行復(fù)制/移動線功能的操作,繪制礦體外推界線從而構(gòu)建外推線框。
數(shù)據(jù)處理,根據(jù)已構(gòu)建的線框,對鉆孔樣品等數(shù)據(jù)進行圈定,同時對樣品數(shù)據(jù)文件進行更新、修改。其主要目的是確定哪些樣品在礦體內(nèi)且要參與估值,同時對異常數(shù)據(jù)進行處理。
關(guān)于如何在三維空間內(nèi)挑選出,處于礦體線框內(nèi)的所有樣品,我們可以通過生成井中坐標功能,將所有鉆孔樣品線段壓縮為點,同時確定其三維坐標XYZ值;根據(jù)樣品點坐標,便能在空間上圈定哪些樣品在礦體線框內(nèi),并通過線框賦值給點的方式,在數(shù)據(jù)文件中反應(yīng)出相應(yīng)的樣品,通過增加過濾條件將滿足賦值要求的點提取出來[4]。針對這些樣品點,我們可以快速得到礦體的平均品位及與品位相關(guān)的數(shù)據(jù),如標準差、方差等,并進行初步分析;同時能獲得更為關(guān)鍵的、可以反映礦體品位均勻程度的變化系數(shù)或西舍爾系數(shù),從而判斷是否存在特高品位需要處理。
特高品位處理,國內(nèi)多采用6-8倍替換法,即使用礦體的平均品位替換超過平均品位6-8倍的樣品品位;國際上一般采用97.7%(或97.5%)分位數(shù)法,即采用礦體內(nèi)所有樣品品位的97.7%(或97.5%)分位數(shù)品位,替換超過該品位的樣品品位,同時可輔以QQ圖進行補充描述。
組合樣品,地質(zhì)統(tǒng)計學一般要求樣品等長組合。在數(shù)據(jù)處理的過程中,原則上要保持大多數(shù)樣品不變,這里可以采用統(tǒng)計直方圖的形式,定量的確定組合樣品的較合理的組合長度。地質(zhì)統(tǒng)計學之所以有這樣要求,一部分原因是在進行克里格法估值或者距離冪次反比法估值時,都是針對樣品點進行的,從而需要相同的長度來消除長度對樣品點品位的影響,以得到更合理、準確的估值。
建立塊體模型的目的在于將礦體分割成多個小塊體,分別計算各個小塊體的資源量,通過累計的方式計算整個礦體的資源量。
在建立空塊模型時,空塊的體積便是各個小塊體的固有屬性,計算它的資源量還需其體重及各個小塊體的平均品位,體重一般由化驗室給出或者是計算得出的品位與體重的函數(shù),而其平均品位就需要用到克里格法或距離冪次反比法進行估值。這里主要以距離冪次反比法為例。
距離冪次反比法的原理,簡單來說,就是在搜索范圍內(nèi),距離“待估塊體”越近的樣品點,其品位對“待估塊體”品位的影響就越大,所占權(quán)重也越大;冪次是控制所占權(quán)重大小的一個指數(shù),當冪次取的過大時,就意味著距離“待估塊體”最近的樣品點,對其估值具有決定性的作用,較遠的樣品點對估值影響微乎其微,所以一般在對貴金屬或者品位變化較大的礦體進行估值時,我們一般采用3次冪,其他一些礦種多采用2次冪,同時也可以根據(jù)實際情況進行論證,來確定具體的冪指數(shù);反比指得是一種常用的計算方式,這里就是指將距離取冪次后的倒數(shù),作為權(quán)系數(shù)進行加權(quán)平均。距離冪次反比法與克里格法在估值過程中,最主要的區(qū)別在于是否考慮礦體的各向異性和是否進行變異函數(shù)的計算,而當?shù)V體品位的均勻程度達到較均勻以上,即礦體內(nèi)所有樣品點品位的變異系數(shù)小于等于1.5時,可采用距離冪次反比法估值,且二者估值結(jié)果的差異通常不大,所以在國內(nèi),較為簡便的距離冪次反比法更受大家青睞。
在建立空塊模型時,小塊體的大小一般可以根據(jù)礦體規(guī)模的大小及工程間距的大小確定,可取工程間距的1/4、1/3、1/2,亦可根據(jù)實際情況、采礦工藝等因素選定;同時有兩種空塊模型可以建立,一種為次分塊模型(圖1),一種為塊系數(shù)模型(圖2)。兩種模型的區(qū)別在于礦體邊界處,空塊是否需要繼續(xù)分割為更小的塊體:次分塊可以設(shè)置繼續(xù)分割為多少小塊體,一般設(shè)置為2*2*2;而塊系數(shù)是保持礦體邊界處,空塊大小不變,以礦體占據(jù)空塊體積的比例,確定礦體邊界處空塊內(nèi)礦體的大小,同時自動生成塊系數(shù)字段。
圖1 次分塊模型
圖2 塊系數(shù)模型
確定搜索范圍,我們需要先建立一個搜索橢球體,其作用在于根據(jù)礦體的整體形態(tài),確定對“待估塊體”進行估值時,所用樣品點的方位與比例,即建立一個規(guī)則的、按比例縮小的礦體。一個橢球的空間形態(tài)可由6個參數(shù)確定,通俗的講可以由橢球的長、寬、厚,及它們的指向確定。構(gòu)建橢球體時,各個參數(shù)是比例關(guān)系,所以橢球可大可小,保持形態(tài)比例不變即可,但搜索范圍=給定的搜索半徑參數(shù)*各個軸的比例系數(shù),故而一般將主軸(即橢球長的方位或礦體走向)的比例系數(shù)設(shè)置為“1”,方便估值搜索范圍的確定。
根據(jù)搜索范圍,給空塊模型進行距離冪次反比法估值,通常估值不是一次性完成,我們需要根據(jù)資源量類型的劃定要求進行多次估值,一般由三次有效估值完成。
根據(jù)地質(zhì)統(tǒng)計學規(guī)范,在資源量類別的劃分過程中,應(yīng)根據(jù)搜索半徑、工程數(shù)、樣品數(shù)等因素對礦塊進行資源量分類,但由于規(guī)范中并沒有明確限定資源量類別的劃分參數(shù),所以可采用經(jīng)驗值來劃分。如在只考慮搜索半徑和工程個數(shù)的情況下采用:搜索范圍為小于1/2控制工程間距且參與估算的工程數(shù)大于等于3個孔時,可劃分為探明的資源量;搜索范圍為小于控制工程間距且參與估算的工程數(shù)大于等于2個孔時,可劃分為控制的資源量;搜索范圍為小于2倍控制工程間距且參與估算的工程數(shù)大于等于1個孔時,可劃分為推斷的資源量[5-9]??刂乒こ涕g距可根據(jù)實際情況放大到相應(yīng)規(guī)范要求的1.25倍左右,亦可通過相關(guān)論證確定最終的資源量分類參數(shù)。但類似的資源量分類標準,易產(chǎn)生“孤島”現(xiàn)象,即不同資源量類型的礦體在空間上相互穿插或包含,這時可以通過多段線賦值、線框賦值的方式,根據(jù)工程控制程度和地質(zhì)可靠程度進行人為干預(yù),采用我國的資源量分類規(guī)范進行分類。
生成塊體模型報告為資源量估算流程中出成果的一步。在生成塊體模型報告時,需通過設(shè)置默認體重或增加體重字段的方式給各個礦塊賦予相應(yīng)的體重值,結(jié)合各個礦塊的體積和先前通過克里格法或距離冪次反比法所得的各個礦塊的品位,就可以依據(jù)需求,設(shè)置不同元素、品位、高度、資源量類別等條件,生成不同的報告模式,實現(xiàn)對礦體資源量快速的、多樣的、全面的描述。這里我們需要注意的是,要根據(jù)不同的礦塊模型設(shè)置不同的參數(shù),如塊系數(shù)模型需要及時的勾選上塊系數(shù)字段,否則會造成估算資源量遠大于實際資源量的情況。同時我們可以在線框菜單欄下找到生成品位-噸位報告功能,通過線框和線框內(nèi)樣品點來直接估算礦體資源量,進行簡單的交叉驗證,通常兩者相差不大。
地質(zhì)工作中,一套理論的形成或者一次資源量的估算,都是經(jīng)過長時間對數(shù)據(jù)的分析、更新、探索而得出的。傳統(tǒng)的工作方法效率較低,且更新、修改的過程重復(fù)性較強,需要大量的人力、物力及時間。因此,三維地質(zhì)建模及其可視化、數(shù)據(jù)管理,成為了現(xiàn)代化礦業(yè)構(gòu)件中不可或缺的一部分。多年來我國Micromine用戶用實踐證明,通過Micromine軟件進行的資源量估算所得出的結(jié)果與傳統(tǒng)幾何法所計算出的資源量,結(jié)果基本吻合,其相對差在1%~15%之間;同時經(jīng)過大量專家的論證與評審,其功能、可靠性及報告通過率得到了保證。Micromine軟件作為大型的專業(yè)礦業(yè)軟件之一,對地勘單位在資源量估算、找礦研究等工作,可以提供一定的幫助。此外,對提升礦山的現(xiàn)代化建設(shè)和數(shù)字化管理也可以起到一定的作用。