謝徽 高幫飛 李寒濱 張書琛 陳偉康

摘要：常福龍金礦床為受剪切帶控制的低溫?zé)嵋好}型金礦床。由于礦化的三維連續(xù)性較差，傳統(tǒng)手工圈連礦體難以準(zhǔn)確刻畫礦化的真實分布。隱式地質(zhì)建模具有速度快、可重復(fù)性強(qiáng)、以數(shù)學(xué)為依據(jù)、人為干擾少等特點，尤其適合處理樣品數(shù)據(jù)稠密且表面形態(tài)復(fù)雜的問題。通過引入徑向基函數(shù)（RBF），對常福龍金礦床開展了三維隱式地質(zhì)建模與品位估值，旨在探討如何更好地解決復(fù)雜不連續(xù)礦體的圈定和估值問題。研究表明，在適當(dāng)?shù)募s束條件下，隱式建模可以給出合理的地質(zhì)模型。RBF估值與地質(zhì)統(tǒng)計學(xué)方法的估值準(zhǔn)確度相當(dāng)，而前者呈現(xiàn)出的礦化結(jié)構(gòu)更為清晰。研究工作對于促進(jìn)隱式地質(zhì)建模技術(shù)在中國資源量估算領(lǐng)域的應(yīng)用有著一定的啟示意義。

關(guān)鍵詞：隱式地質(zhì)建模;徑向基函數(shù);品位估值;資源量估算;地質(zhì)統(tǒng)計學(xué)

引 言

常福龍金礦床位于華北地臺北緣大青山南麓，為受剪切帶控制的低溫?zé)嵋好}型金礦床[1-2]。礦體產(chǎn)于北西向—北西西向剪切帶中，礦化以浸染狀、細(xì)脈狀、網(wǎng)脈狀、角礫狀為主，單個礦脈寬度從幾厘米到幾米不等，其三維連續(xù)性整體較差。

目前，傳統(tǒng)礦體的圈定主要依據(jù)相鄰見礦工程及礦脈或容礦構(gòu)造裂隙產(chǎn)狀“就礦連礦”。實踐證明，傳統(tǒng)礦體圈定和資源量估算給出的結(jié)果難以很好指導(dǎo)礦山生產(chǎn)。隱式地質(zhì)建模（Implicit Geological Modelling）是近20年來發(fā)展起來的一項建模新技術(shù)，為當(dāng)前國內(nèi)外地質(zhì)建模的熱點和前沿[3-16]。其原理是基于隱函數(shù)（Implicit Function）對離散的三維點數(shù)據(jù)插值而實現(xiàn)表面重建（Surface Reconstruction）。隱式地質(zhì)建模具有速度快、可重復(fù)性強(qiáng)、以數(shù)學(xué)為依據(jù)、人為干擾少等特點[17]，尤其適合處理樣品點數(shù)據(jù)稠密且表面形態(tài)復(fù)雜的問題[18-20]。

本文嘗試采用徑向基函數(shù)（Radial Basis Functions，RBF）開展常福龍金礦床品位建模和估值，并將估值結(jié)果與傳統(tǒng)IDW和Kriging法進(jìn)行比較。其目的是探討及利用隱式建模技術(shù)，解決復(fù)雜不連續(xù)礦化體圈定和估值問題，也以此為契機(jī)與同行交流，為隱式建模技術(shù)在地質(zhì)和礦業(yè)領(lǐng)域的應(yīng)用提供借鑒。

1 地質(zhì)背景

常福龍金礦床位于華北地臺北緣大青山大佘太—固陽—武川區(qū)域東西向深大斷裂的次級斷裂常福龍溝—德勝營—大石槽斷裂內(nèi)[2]。成礦帶劃分上屬于華北地臺北緣金成礦帶的中段，除常福龍金礦床外，該金成礦帶上還分布有烏拉山、東伙房、卯獨慶、白乃廟等幾十處金礦床（點）[21-22]。東西向展布的大青山復(fù)式褶皺及其伴生逆沖斷裂構(gòu)成了區(qū)域構(gòu)造主體格架[23]。其中，逆沖斷裂總長幾十千米，寬幾十米到幾百米。

礦區(qū)內(nèi)F1斷裂走向300°（見圖1-a）），整體南傾，傾角30°～60°，其上盤為太古界烏拉山群片麻巖、大理巖，下盤為元古界渣爾泰群碳質(zhì)板巖、板巖和砂巖（見圖1-b））;其下盤由多條NWW向斷裂組成的剪切帶控制蝕變和礦化[2]。

礦區(qū)蝕變作用成礦早期主要為面狀硅化、綠泥石化，主成礦期以脈狀硅化和黃鐵絹英巖化為主，成礦晚期則疊加面狀或脈狀碳酸鹽化和高嶺土化[1]?？臻g分布上，F(xiàn)1斷裂下盤破碎帶中廣泛發(fā)育面狀硅化、綠泥石化及浸染狀黃鐵礦化，整個蝕變帶長約2 000 m、寬100～400 m。礦化地段多疊加細(xì)脈狀、網(wǎng)脈狀硅化和黃鐵絹英巖化;碳酸鹽化和高嶺土化主要分布在破碎帶兩側(cè)圍巖中。主成礦期的硅化以玉髓（見圖2-d）為主，呈皮殼狀、犬齒狀和梳狀分布在張性脈中，顯示低溫蝕變特征。

礦區(qū)內(nèi)礦體產(chǎn)出主要受北西西向剪切帶及伴生次級斷裂控制。主礦體由含礦構(gòu)造透鏡體和角礫巖組成，主礦體間更次級的斷裂和裂隙控制著硅質(zhì)脈狀、網(wǎng)脈狀和角礫狀礦化（見圖2-a、b、c）。從礦化空間分布來看，淺部礦化連續(xù)性好，礦體厚大，引張構(gòu)造和局部膨脹地段多;向深部，沿走向和傾向礦化連續(xù)性減弱，礦體變薄。這可能與張性斷裂-裂隙規(guī)模較大及在淺部較為密集和發(fā)育有關(guān)，強(qiáng)礦化發(fā)育在正花狀節(jié)理系花蒂之上的“金三角”區(qū)[24]。

整體而言，常福龍金礦床容礦剪切帶發(fā)育不成熟，流體作用以低溫蝕變?yōu)橹?，而缺少中高溫流體疊加改造，礦化較為分散[2]。圈連礦體時，主要依據(jù)相鄰工程見礦部位，參考礦脈及容礦構(gòu)造產(chǎn)狀“就礦連礦”。前期工作按照Ⅱ類+Ⅲ類勘查類型，以0.8 g/t為邊界品位共圈定了36條礦體（低品位礦體+工業(yè)礦體）。單個礦體長24～440 m，寬1.00～14.46 m。礦體走向290°～300°，傾向200°～210°，傾角（70°～90°）較陡。礦體呈脈狀、不規(guī)則透鏡體狀，局部呈現(xiàn)分支復(fù)合和尖滅再現(xiàn)特征。礦體的三維連續(xù)性較差，地表鉆探工程揭露礦體后，用坑道工程去追索和驗證礦體，有時難以達(dá)到預(yù)期效果。

2 勘查數(shù)據(jù)處理與分析

常福龍金礦床在勘查和開發(fā)階段施工了眾多探礦工程，積累了大量樣品數(shù)據(jù)，包括47個地表鉆孔、30個坑內(nèi)鉆探、5個中段坑道工程，總進(jìn)尺37 184 m，基本樣長1 m，累計基本分析樣品19 782件?；痉治鲈诔８｝埥鸬V化實驗室采用原子吸收光譜法完成，樣品檢出限為0.01 g/t，可以滿足數(shù)據(jù)統(tǒng)計分析的精度需求。為了避免樣品支集（Support）對建模的影響[25-26]，將鉆孔數(shù)據(jù)庫中部分0值及缺省值統(tǒng)一用檢出限的一半（0.005 g/t）來代替。處理后的樣品個數(shù)為22 497件，品位均值為0.21 g/t，最小值為0.005 g/t，最大值為148.29 g/t，標(biāo)準(zhǔn)差為1.88 g/t，變異系數(shù)（Cv）為9.06，偏度（Skewness）為43.1（見圖3-a））。樣品數(shù)據(jù)表現(xiàn)出很強(qiáng)的變化性。

由圖3-a）可知：曲線明顯的斷點為0.1 g/t和45.9 g/t，局部出現(xiàn)多個拐點或斷點，難以憑經(jīng)驗判斷樣品的特異值。本文采用分形方法[27-29]確定的特低值和特高值分別為0.09 g/t和19.0 g/t（見圖3-b）），將分形區(qū)間0.09～19.0 g/t作為后續(xù)品位估值的樣品區(qū)間。低于0.09 g/t的值將不會參與估值，而高于19.0 g/t的樣品值將由19.0 g/t來代替。ED1E4734-F179-45FD-8F81-33248B016491

3 隱式建模原理

COWAN等[17]首次提出隱式建模（Implicit Modelling）概念。傳統(tǒng)地質(zhì)建模方法是在2D剖面上手工數(shù)字化來定義礦巖邊界（見圖4-a、b），然后利用網(wǎng)格化方法創(chuàng)建3D表面[30]。隱式建模則是一種無網(wǎng)格化（Meshless）表面重建方法，采用體積函數(shù)（Volume Function）中提取的特定等值面（Isosurface）來刻畫目標(biāo)表面[19]。對于空間任一點P（x，y，z），其表達(dá)式為f（P）=C，其中f（P）為體積函數(shù)（即隱函數(shù)，Implicit Function），C為常數(shù)。為了構(gòu)建0值等值面，通常將落在面上的所有接觸點或節(jié)點均賦值為0，即f（P）=0;點P落入體積內(nèi)則賦為負(fù)值（f（P）<0，標(biāo)記為-ve），體積外的點賦為正值（f（P）>0，標(biāo)記為+ve）（見圖4-c）。利用隱函數(shù)估值算法可以自動獲得0值等值面形態(tài)（見圖4-d）。該方法的最大優(yōu)點是省去了大量手工數(shù)字化過程，直接根據(jù)巖石屬性或品位的數(shù)值化賦值結(jié)果來生成表面[17]，這一過程分別稱作巖性建模和品位建模。

徑向基函數(shù)（Radial Basis Function，RBF）是一種常用的隱函數(shù)，最早由HARDY[31]引入，用來處理地形數(shù)據(jù)。其基本表達(dá)式[32]為：

f（x）=∑Ni=1λi（‖x-xi‖）+p（x）（1）

式中：f（x）為RBF插值函數(shù);i=1，2，…，N;λi為線性組合權(quán)重系數(shù);‖x-xi‖為任何一點x到已知點xi的歐氏距離，xi為基函數(shù)的中心;p（x）為線性多項式。

因此，所有已知節(jié)點（node）的基函數(shù)進(jìn)行加權(quán)線性組合可以得到相應(yīng)的插值函數(shù)，該插值函數(shù)可以實現(xiàn)對任何曲面的逼近[33]。常用的基函數(shù)包括：薄板樣條（Thin-plate spline）φ（r）=r2lgr;高斯（Gaussian）φ（r）=e-cr2;多元二次曲面（Multi-Quadric）φ（r）=r2+c2 ;線性（Linear）φ（r）=r;立方（Cubic）φ（r）=r3（r=‖x-xi‖，c為常數(shù)）。不難看出，其函數(shù)僅與離中心點的距離有關(guān)。任何離中心點距離相等的點，其函數(shù)值相同，即為中心點球狀對稱，因而稱作徑向（Radial）。由于其簡單的形式和較高的逼近行為，RBF在過去20年成為了離散數(shù)據(jù)插值的有效工具[3-4，8，18，32]。國際通用三維礦業(yè)軟件如Leapfrog、Micromine、Minesight、Datamine等也采用RBF或改進(jìn)的RBF作為隱式建模的插值函數(shù)。

4 品位建模

常福龍金礦床礦化具有“面狀蝕變+細(xì)脈狀礦化”特征，這些高品位細(xì)脈的規(guī)模、產(chǎn)狀和分布在一定程度上決定了礦化強(qiáng)弱和工業(yè)礦體的分布。由Au品位大于或等于0.8 g/t樣品在三維空間的分布（見圖5）可知：工業(yè)礦化樣品在淺部及走向上連續(xù)性相對較好，而在傾向上變化較大，向深部樣品明顯變得稀疏和不連續(xù)。礦化的極不連續(xù)，使得傳統(tǒng)手工圈礦只能依賴于僅有的地質(zhì)產(chǎn)狀信息“就礦連礦”。

隱式建?？梢园堰@些難題交給數(shù)學(xué)模型。經(jīng)驗表明，只要約束條件合適，隱式建模就可以給出符合“地質(zhì)邏輯”的插值結(jié)果[3，7，13，34]。對常福龍金礦床原始樣品數(shù)據(jù)進(jìn)行全孔1 m等長組合作為品位建模的數(shù)據(jù)點。本次隱式建模的約束條件主要有：①用地形表面約束模型頂部界線;②設(shè)置插值橢球參數(shù)走向300°、傾角80°、傾伏角0°，長、中、短軸長度比為1∶0.1∶0.05;③插值函數(shù)選取RBF球狀（spheroidal）模型，設(shè)置基本變程為200 m，形狀參數(shù)為3，塊金值為0.15。根據(jù)研究需要，同時生成0.09 g/t（礦化域邊界）、0.8 g/t（低品位礦體邊界）和1.0 g/t（工業(yè)品位礦體邊界）的品位殼（grade shell）或線框（Wireframe，WF），分別記為WF_0.09、WF_0.8和WF_1.0，其插值精度分別為4 m、2 m和1 m。WF_0.8和WF_1.0的三維隱式品位模型見圖6，7勘探線三維隱式品位模型剖面見圖7。

由圖6、圖7可知，品位線框模型呈走向延長大于傾向延伸的“梭形”。這與礦化受張剪性構(gòu)造控制的特征相符，且與礦化在近地表富集向深部減弱特點一致（見圖1-b））。盡管也可以通過調(diào)整橢球中軸與長軸的比值來增強(qiáng)線框在傾向方向的延伸（可在上、下中段相應(yīng)位置均不見礦的情況下，延伸50～60 m，而按照國內(nèi)規(guī)范有限外推要求，這種情況下延伸應(yīng)在40 m以內(nèi)），但這樣明顯與實際坑道探礦情況相悖。因此，經(jīng)多次試驗后，堅持將橢球中軸設(shè)置成較小的數(shù)值。

WF_0.09和WF_1.0線框內(nèi)樣品數(shù)據(jù)統(tǒng)計結(jié)果見圖8。由圖8可知：隨著圈礦邊界的提高，平均值和標(biāo)準(zhǔn)差呈明顯上升趨勢，而品位變化系數(shù)呈下降趨勢。以1.0 g/t為邊界，701個樣品數(shù)據(jù)的均值達(dá)到了3.89 g/t。另外，由于插值精度設(shè)置問題，仍有部分低于1.0 g/t的樣品被圈入線框，因而統(tǒng)計直方圖顯示出“多峰”分布特征（見圖8-b））。

5 品位估值

以上述隱式建模圈定的2個線框（WF_0.09和WF_1.0）作為Au品位估值的品位域。在對品位域開展進(jìn)一步邊界分析和變異函數(shù)分析基礎(chǔ)上進(jìn)行估值。

5.1 邊界分析

品位域邊界分析的目的是判斷采用域外一定范圍內(nèi)數(shù)據(jù)參與估值，還是僅采用域內(nèi)數(shù)據(jù)[35-37]。以0.09 g/t和1.0 g/t為邊界品位的品位域邊界分析見圖9。由圖9可知，品位域內(nèi)品位數(shù)據(jù)在接近邊界時，呈明顯降低趨勢，表現(xiàn)為“半軟”邊界特點?？紤]到品位域內(nèi)、外品位數(shù)值有顯著變化，同時品位建模時已經(jīng)將部分低于邊界品位（0.09 g/t和1.0 g/t）的樣品點圈入線框（見圖8），因此估值時采用硬邊界條件，即僅用線框內(nèi)樣品數(shù)據(jù)進(jìn)行估值。

5.2 變異函數(shù)分析ED1E4734-F179-45FD-8F81-33248B016491

以球狀模型對WF_0.09和WF_1.0線框內(nèi)樣品數(shù)據(jù)開展變異函數(shù)分析，結(jié)果見圖10。擬合結(jié)果如下：①WF_0.09線框內(nèi)，Au品位為0.09～19.0 g/t的樣品數(shù)據(jù)（去特異值后）設(shè)為標(biāo)準(zhǔn)差方差為2.18，設(shè)置標(biāo)準(zhǔn)化基臺值為1.0，則塊金值（C0）為0.18（即塊金效應(yīng)為18 %），長軸、中軸、短軸方向的變程（a）分別為89.85 m、44.75 m和3.23 m（見圖10-a）、b）、c））;②WF_1.0線框內(nèi)，Au品位為0.09～19.0 g/t的樣品數(shù)據(jù)（去特異值后）樣品方差為20.44，設(shè)置標(biāo)準(zhǔn)化基臺值為1.0，則塊金值為0.25（即塊金效應(yīng)為25 %），長軸、中軸、短軸方向的變程分別為179.30 m，54.21 m和19.49 m（見圖10-d）、e）、f））。此外，從實驗變異函數(shù)來看，3個方向上不同程度地出現(xiàn)了空穴效應(yīng)（Hole Effect），可能與統(tǒng)計方向上礦體三維延伸規(guī)模不大的礦化呈周期性分布有關(guān)。實驗變異函數(shù)結(jié)構(gòu)性總體較差，也反映了三維礦化分布的復(fù)雜性。

5.3 估值結(jié)果

按照8 m×4 m×2 m對WF_0.09和WF_1.0線框進(jìn)行塊劃分，建立塊模型（Block Model，BM），分別記作BM_0.09和BM_1.0。采用距離冪次反比（IDW）、普通克里格（OK）和徑向基函數(shù)（RBF）3種方法分別對塊模型進(jìn)行估值?？紤]到塊金效應(yīng)對IDW估值冪次的影響[35，38-39]，塊金效應(yīng)為18 %時，取冪次為3;塊金效應(yīng)為25 %時，取冪次為2。采用IDW和OK估值時，設(shè)置橢球估值樣品點數(shù)最少為4，最多為20。不同塊模型和插值方法的估值結(jié)果分別標(biāo)記為IDW_3_0.09、IDW_2_1.0、OK_0.09、OK_1.0、RBF_0.09和RBF_1.0。

BM_1.0塊模型的估值結(jié)果見圖11。由圖11可知：IDW、OK和RBF 3種方法估計的均值分別為3.41 g/t 、2.70 g/t和1.95 g/t，呈明顯下降趨勢。從塊模型（見圖11-a、c、e）也可以看出，IDW法高品位數(shù)據(jù)的連續(xù)性最好，其次為OK法，RBF高品位數(shù)據(jù)連續(xù)性最差。這一點也可以從塊品位的統(tǒng)計直方圖看出，IDW和OK品位數(shù)據(jù)呈現(xiàn)典型的“正態(tài)”分布特征（見圖11-b、d），而RBF低品位塊體占比較大（見圖11-f），因而平均品位較前二者明顯要低。從統(tǒng)計數(shù)據(jù)還可以看出，OK估值標(biāo)準(zhǔn)差和最大值都較IDW和RBF低，顯示出“平滑”作用特點。OK估值中部分負(fù)值塊體（最小值為-0.70 g/t），可能與屏蔽作用導(dǎo)致部分樣品點出現(xiàn)負(fù)權(quán)重有關(guān)[39-40]。

6 討 論

6.1 估值結(jié)果與組合樣品比較

為對比1 m等長組合樣品數(shù)據(jù)和塊模型結(jié)果，分臺階統(tǒng)計了大于工業(yè)品位1.0 g/t的平均品位（見圖12）。其中，1 m_組合為剔除了特高值（19.0 g/t）的樣品點數(shù)據(jù)。考慮到塊模型估值結(jié)果中均值最高在3.41 g/t，對應(yīng)m+3σ（m為平均值，σ為標(biāo)準(zhǔn)差）最大值為7.16 g/t，根據(jù)對礦床的整體認(rèn)識，有必要再對“特高值”進(jìn)行處理（傳統(tǒng)方法進(jìn)行特高值處理也是采用迭代法）。1 m_組合*為再次剔除了大于10 g/t數(shù)據(jù)后的樣品點：1858臺階剔除33個樣品點（Au平均品位為16 g/t）;1818臺階剔除16個樣品點（Au平均品位為15 g/t）;1778臺階沒有超過10 g/t樣品點，未進(jìn)行調(diào)整;1738臺階剔除4個樣品點（Au平均品位為16 g/t）;1698臺階剔除2個樣品點（Au平均品位為12 g/t）;1658臺階剔除1個樣品點（Au平均品位為14 g/t）。

由圖12可知：采用IDW法和OK法時，BM_1.0為塊模型時的估值平均品位較高（分別為3.46 g/t和2.82 g/t），而以BM_0.09為塊模型的估值平均品位最低（分別為1.78 g/t和1.72 g/t）;RBF估值結(jié)果居于中間位置，2種條件下的平均品位基本一致（分別為2.35 g/t和2.32 g/t）。換言之，估算域范圍對IDW和OK估值影響較大，以0.09 g/t圈礦（WF_0.09）導(dǎo)致大量低品位數(shù)據(jù)參與估值，大大降低了塊模型的均值。從不同臺階估值差異來看，從淺部1 858 m中段到深部1 658 m中段，不同方法表現(xiàn)出的均值差異呈加大趨勢。在1 858 m中段，礦化較為集中，礦體空間連續(xù)性較好，大于工業(yè)品位樣品數(shù)據(jù)相對較多，這種情況下OK_1.0和RBF_1.0的估值結(jié)果幾乎一致。向深部，隨著礦化連續(xù)性減弱，可供估值樣品點數(shù)據(jù)顯著減少，這種一致性不復(fù)存在。

與組合樣品的對比也可以反映不同方法的估值特點。IDW法在1 818 m中段和1 738 m中段，以及OK法在1 698 m中段和1 658 m中段的估值結(jié)果與1 m_組合樣品數(shù)據(jù)完全吻合;而1 858 m中段和1 778 m中段，IDW法和OK法都表現(xiàn)出對“極高值”和“極低值”的平滑作用，即高值低估和低值高估。RBF法在1 778 m中段與1 m_組合樣品數(shù)據(jù)最為接近。該中段組合樣品數(shù)據(jù)最高值為5 g/t，這表明品位中的高值數(shù)據(jù)對RBF估值的影響較大。RBF_0.09估值結(jié)果在不同臺階的品位變化趨勢與1 m_組合樣品基本一致。再次剔除特高值影響后的1 m_組合*樣品數(shù)據(jù)均值為2.40 g/t，與RBF法的估值結(jié)果較為接近。尤其是RBF_0.09，除了1 698 m中段與組合樣品數(shù)據(jù)有差異外，其余中段二者均值基本吻合。

6.2 估值結(jié)果與前人工作比較

為了進(jìn)一步說明RBF隱式建模和估值的效果，將塊模型估值結(jié)果與前人手工圈定的礦體進(jìn)行對比。低品位（0.8 g/t）和工業(yè)品位（1.0 g/t）1 m等長組合樣品數(shù)據(jù)在1 858 m中段的投影（視域范圍為上下各延伸4 m）見圖13，底圖為手工圈連的礦體中段平面圖。BM_0.09和BM_1.0塊模型見圖14，1 858 m中段進(jìn)行了IDW、OK和RBF 3種方法的估值結(jié)果對比，圖中僅顯示了低品位礦體和工業(yè)礦體的分布，即大于或等于0.8 g/t的塊模型估值結(jié)果。ED1E4734-F179-45FD-8F81-33248B016491

圖13 1 858 m中段工業(yè)礦化樣品分布? 從塊模型反映的礦化結(jié)構(gòu)來看，對于3種方法而言，BM_0.09塊模型結(jié)構(gòu)性明顯較好，表現(xiàn)為礦脈的分支復(fù)合較為清晰，而BM_1.0的塊模型連續(xù)性較好，但破壞了礦化的結(jié)構(gòu)特征（見圖14）。此外，部分地段的塊模型因原始工業(yè)品位樣品點過少而未給出估值結(jié)果。對比基于BM_0.09塊模型的不同方法，RBF反映的礦化結(jié)構(gòu)最為清晰，且個別單工程見礦地段也在塊模型中顯示。這無疑對后續(xù)找礦勘查提供了思路和方向。

6.3 隱式建模的地質(zhì)約束

隱式地質(zhì)建模的最大特點是速度，而不是準(zhǔn)確性、無偏性、細(xì)節(jié)[20]?？焖俳？梢宰尩刭|(zhì)工程師將更多精力用在考慮模型的地質(zhì)約束上[17]。

傳統(tǒng)建模方法主要依賴于地質(zhì)工程師的認(rèn)識，隱式地質(zhì)建模亦是如此。無各向異性橢球約束和將橢球中軸與長軸

比值由0.1改為0.6（其他參數(shù)不變）的隱式建模結(jié)果見圖15、圖16。沒有各向異性約束情況下，RBF模型邊界多為“氣球”形狀，明顯失真;而改變橢球參數(shù)則直接增強(qiáng)了工業(yè)品位礦體在傾向上的連續(xù)性，與坑道探礦實際情況不符。

地質(zhì)體模型完全或部分使用隱式建模工具來構(gòu)建，并不意味著其比使用其他建模工具創(chuàng)建的地質(zhì)體模型更好或更差[20]。與傳統(tǒng)建模方法一樣，隱式建模也需要建模者的經(jīng)驗和必要地質(zhì)條件約束。然而，如何憑借地質(zhì)工程師的經(jīng)驗，準(zhǔn)確判斷礦化在走向和傾向的連續(xù)性，評價礦化的變異性，以及界定特高值等，都不是輕易可以回答的問題。需要建模者開展系統(tǒng)的野外勘查，還要與現(xiàn)場地質(zhì)工程師充分交流，才可能對某一礦床的礦化分布規(guī)律、品位變化乃至礦床成因有較為清晰的認(rèn)識。隱式建模也并非任何時候都能替代傳統(tǒng)建模方法，正如英國著名統(tǒng)計學(xué)家BOX[41]所說的“所有的模型都是錯誤的，但有些是有用的”。隱式建模提供了一個新的數(shù)學(xué)工具，可以讓地質(zhì)工程師在較短時間內(nèi)，根據(jù)自身對礦床的認(rèn)知，設(shè)置和調(diào)整約束參數(shù)，不斷地“試錯”，最終構(gòu)建出最具實際意義的模型。

7 結(jié) 語

地質(zhì)行業(yè)位處礦業(yè)開發(fā)的上游和工業(yè)鏈的頂端。地質(zhì)模型和資源量估算的可靠性在一定程度上決定了后續(xù)礦山開發(fā)的風(fēng)險和收益。因此，在礦山投資開發(fā)之前，投入足夠多的時間和精力，創(chuàng)建符合實際情況的地質(zhì)模型及開展可靠品位估值十分必要。針對內(nèi)蒙古常福龍金礦床礦化不連續(xù)、手工圈連難的實際情況，引入基于徑向基函數(shù)的隱式建模和插值方法，為復(fù)雜地質(zhì)問題提供了新的解決方案。

隱式建模真正引入地質(zhì)行業(yè)不超過20年，國內(nèi)也僅有少數(shù)學(xué)者在對這一課題進(jìn)行研究，目前尚屬于起步階段。然而，就像地質(zhì)統(tǒng)計學(xué)在國內(nèi)的推廣和使用一樣，可能還需要幾十年時間，才能寫入勘查規(guī)范，成為地質(zhì)行業(yè)的工作準(zhǔn)則。本文作為隱式建模技術(shù)應(yīng)用的一個嘗試，希望能夠起到借鑒作用，同時也期待與廣大地質(zhì)工作者一道，為推動中國資源儲量估算工作邁上新臺階貢獻(xiàn)出自己的力量。

[參 考 文 獻(xiàn)]

[1] 陳志廣，高幫飛，黃榮偉，等.內(nèi)蒙古常福龍金礦床地質(zhì)特征、成礦模式及其地質(zhì)意義[J].黃金，2012，33（6）：11-16.

[2] 謝徽，高幫飛，陳志廣，等.內(nèi)蒙古常福龍剪切帶型金礦床形成的構(gòu)造控制[J].黃金，2015，36（1）：9-14.

[3] HILLIER M J，SCHETSELAAR E M，DEKEMP E A，et al.Three-dimensional modelling of geological surfaces using generalized interpolation with radial basis functions[J].Mathematical Geoscience，2014，46（8）：931-953.

[4] HILLIER M J，DE KEMP E A，SCHETSELAAR E M.Implicit 3-D modelling of geological surfaces with the generalized radial basis functions （GRBF） algorithm[R].Ottawa：Geological Survey of Canada，2017：1-15.

[5] VOLLGGER S A，CRUDEN A R，AILLERES L，et al.Regional dome evolution and its control on ore-grade distribution：insights from 3D implicit modelling of the Navachab gold deposit，Namibia[J].Ore Geology Reviews，2015，69：268-284.

[6] MARTIN R，BOISVERT J B.Iterative refinement of implicit boundary models for improved geological feature reproduction[J].Computers & Geosciences，2017，109：1-15.

[7] 郭甲騰，吳立新，周文輝.基于徑向基函數(shù)曲面的礦體隱式自動三維建模方法[J].煤炭學(xué)報，2016，41（8）：2 130-2 135.

[8] MAJDISOVA Z，SKALA V.Big geo data surface approximation using radial basis functions：a comparative study[J].Computers and Geosciences，2017，109：51-58.

[9] 鄒艷紅，黃望，陽寬達(dá)，等.基于楊赤中推估法空間插值的三維地質(zhì)隱式建模[J].地質(zhì)學(xué)刊，2017，41（3）：384-393.ED1E4734-F179-45FD-8F81-33248B016491

[10] 鄒艷紅，李高智，毛先成，等.基于隱函數(shù)曲面的三維斷層網(wǎng)絡(luò)建模與不確定性分析[J].地質(zhì)論評，2020，66（5）：1 349-1 360.

[11] WANG J M，ZHAO H，BI L，et al.Implicit 3D modeling of ore body from geological boreholes data using hermite radial basis functions[J].Minerals，2018，8（10）：443.

[12] WANG C，WANG G，LIU J J，et al.3D geochemical modeling for subsurface targets of Dashui Au deposit in Western Qinling（China）[J].Journal of Geochemical Exploration，2019，203：59-77.

[13] ZHONG D Y，WANG L G，BI L，et al.Implicit modeling of complex orebody with constraints of geological rules[J].Transactions of Nonferrous Metals Society of China，2019，29（11）：2 392-2 399.

[14] GUO J T，WU L X，ZHOU W H，et al.Section-constrained local geological interface dynamic updating method based on the HRBF surface[J].Journal of Structural Geology，2018，107：64-72.

[15] GUO J T，WANG J M，WU L X，et al.Explicit-implicit-integrated 3-D geological modelling approach：a case study of the Xianyan Demo-lition Volcano （Fujian，China）[J].Tectonophysics，2020，795：228648.

[16] GUO J T，WU L X，WANG J M，et al.Three-dimensional geological modelling and spatial analysis from geotechnical borehole data using an implicit surface and marching tetrahedra algorithm[J].Engineering Geology，2021，284：106047.

[17] COWAN E J，BEATSON R K，ROSS H J，et al.Practical implicit geological modelling[C]∥AusIMM.5th International Mining Geo-logy Conference.AusIMM，2003：89-99.

[18] STEWART M，DE LACEY J，HODKIEWICZ P F，et al.Grade estimation from Radial Basis Functions—How does it compare with conventional geostatistical estimation？[C]∥AusIMM.9th International Mining Geology Conference.Bendigo：AusIMM，2014：129-139.

[19] KNIGHT R H.Orebody solid modelling accuracy—A comparison of explicit and implicit modelling techniques using a practical example from the Hope Bay District，Nunavut，Canada[C]∥AusIMM.6th International Mining Geology Conference.Bendigo：AusIMM，2006：175-183.

[20] KENTWELL D J.Destroying the distinction between explicit and implicit geological modelling[J].The AusIMM Bulletin，2014，6：54-56.

[21] 孫德育，梁一鴻，張業(yè)明.內(nèi)蒙武川—固陽—大佘太韌性剪切帶與金礦[J].長春地質(zhì)學(xué)院學(xué)報，1990，20（4）：399-406，428.

[22] 鄭翻身，徐國權(quán)，馮貞，等.內(nèi)蒙古中部地區(qū)綠巖型金礦地質(zhì)特征及成礦遠(yuǎn)景預(yù)測[J].地質(zhì)學(xué)報，2005，79（2）：232-248.

[23] 劉正宏，徐仲元，楊振升.內(nèi)蒙古大青山印支運動厘定[J].地質(zhì)論評，2003，49（5）：457-463.

[24] 汪勁草，王永東，程雄衛(wèi).正花狀節(jié)理系對內(nèi)蒙古常福龍金礦床的控制[J].桂林理工大學(xué)學(xué)報，2012，32（4）：447-454.

[25] 郭甲騰.地礦三維集成建模與空間分析方法及其應(yīng)用[D].沈陽：東北大學(xué)，2013.

[26] ZABRUSKY K.Wireframe construction in resource estimation：methods and practices[C]∥SME.SME Annual Meeting.Englewood：SME，2015：1-5.ED1E4734-F179-45FD-8F81-33248B016491

[27] BLENKINSOP T G.Cataclasis and processes of particle size reduction[J].Pure and Applied Geophysics，1991，136（1）：59-86.

[28] 高幫飛，陳志廣，孫剛，等.內(nèi)蒙古常福龍金礦床Au品位分布特征及其地質(zhì)意義[J].黃金，2012，33（9）：9-12.

[29] 高幫飛，謝徽，沈陽，等.一種確定資源儲量估算的品位域的方法：CN201611040814.3[P].2017-05-31.

[30] 何建春.地質(zhì)體三維空間形態(tài)隱式模擬的若干技術(shù)研究[D].長沙：中南大學(xué)，2011.

[31] HARDY R L.Multiquadric equations of topography and other irre-gular surfaces[J].Journal of Geophysical Research，1971，76（8）：1 905-1 915.

[32] CARR J C，BEATSON R K，CHERRIE J B，et al.Reconstruction and representation of 3D objects with radial basis functions[C]∥POCOK L.Proceedings of the 28th annual conference on computer graphics and interactive techniques.New York：Association for Computing Machinery，2001：67-76.

[33] 周少山.隱式曲面建模和體網(wǎng)格生成方法研究[D].武漢：武漢理工大學(xué)，2018.

[34] BASSON I J，THOMAS S A J，STOCH B，et al.The structural setting of mineralisation at Kolomela Mine，Northern Cape，South Africa，based on fully-constrained，implicit 3D modelling[J].Ore Geology Reviews，2018，95：306-324.

[35] 李紅兵，高幫飛.巖金礦床資源儲量估算品位域的確定及其意義[J].黃金，2021，42（4）：6-11.

[36] ABZALOV M Z.Applied mining geology[M].Berlin：Springer，2016.

[37] 高幫飛.資源量估算的邊界分析與礦化體圈定[J].地質(zhì)與勘探，2018，54（2）：415-425.

[38] COOMBES J.The art and science of resource estimation：a practical guide for geologists and engineers[M].Subiaco：Coombes Capability，2008.

[39] 高幫飛，李紅兵，張書琛，等.塊金效應(yīng)地質(zhì)意義及其對品位估值影響[J].黃金，2021，42（11）：6-13.

[40] BAAFI E Y，KIM Y C，SZIDAROVSZKY F.On nonnegative weights of linear kriging estimation[J].Mining Engineering，1986，38：437-442.

[41] BOX G E P.Science and statistics[J].Journal of the American Statistical Association，1976，71：791-799.

Implicit geological modeling and grade estimation

of the Changfulong Gold Deposit，Inner Mongolia

Xie Hui1，2，Gao Bangfei1，3，Li Hanbin1，3，Zhang Shuchen1，3，Chen Weikang1，3

（1.China Railway Resources Group Co.，Ltd.;

2.Beijing Technical Consulting Branch of China Railway Resources Group Co.，Ltd.;

3.La Miniere de Kalumbuwe Myunga SAS）

Abstract：Changfulong Gold Deposit is a low-temperature hydrothermal vein-type gold deposit controlled by shear zone.Due to the poor three-dimensional continuity of mineralization，it is difficult to accurately depict the true distribution of mineralization by the traditional method of hand-digitizing.Implicit geological modeling has the characteristics of fast speed，strong repeatability，mathematical and less man interference，so it is especially suitable to deal with the problems of dense sample data and complex surface morphology.In this paper，the radial basis function （RBF） is introduced to carry out three-dimensional implicit geological modeling and grade estimation of Changfulong Gold Deposit，intended to explore how to better solve the problem of delimitation and estimation of complex discontinuous mineralized bodies.The results show that implicit modeling can make a reasonable geological model under appropriate constraints.The accuracy of RBF estimates is similar to that of geostatistical methods，while the former shows a clearer mineralization structure.The research work in a way enlightens the promotion of implicit geological modeling technology application to the field of resource estimation in China.

Keywords：implicit geological modeling;radial basis function;grade estimation;resource estimation;geostatisticsED1E4734-F179-45FD-8F81-33248B016491