胡世士，袁飛虎，黎國亮，魏巍，李偉

(1.中國葛洲壩集團水泥有限公司，湖北 武漢 430000；2.葛洲壩宜城水泥有限公司，湖北 襄陽市 441400；3.武漢理正工程科技有限公司，湖北 武漢 430000)

1 礦床地質(zhì)概況

宜城市胡咀-馬頭山礦區(qū)大地構(gòu)造位置處于揚子準地臺東北緣大巴山臺褶帶東段，屬大洪山斷褶束襄陽—三陽店復式背斜之客店坡-板凳崗倒轉(zhuǎn)向斜的西段。礦區(qū)分布地層主要有志留系、二疊系、第4系，其中二疊系下統(tǒng)棲霞組灰?guī)r段第1層、第2層和二疊系下統(tǒng)茅口組第2段為水泥用灰?guī)r含礦層位。礦區(qū)分為兩個礦段，北部為胡咀礦段，南部為馬頭山礦段，兩個礦段相隔400～600m，中間為第4系覆蓋的巖溶洼地。

胡咀礦段構(gòu)造較為簡單，呈一向北東傾斜的單斜構(gòu)造，地層傾向為25°～45°，傾角為34°～43°。馬頭山礦段為一軸向北西的向斜構(gòu)造，南翼地層傾向為7°～70°，傾角為15°～40°，多在30°左右。北翼地層傾向南西180°～280°，傾角為15°～43°。礦區(qū)二疊系下統(tǒng)棲霞組灰?guī)r段第1層（P1q2?1）、第2層（P1q2?2）及二疊系下統(tǒng)茅口組第2層（P1m2）共3個層位為水泥用灰?guī)r含礦層。胡咀礦段礦體（層）總長約2000m，出露標高為130～196m；馬頭山礦段礦體（層）總長約700m，出露標高為110～203m，兩礦段礦體地表出露寬度為40～60m，厚度為30～50m。礦體形態(tài)簡單，呈層狀產(chǎn)出，由于受斷層影響，各礦層沿走向及傾向均不連續(xù)，厚度變化較大。

2 理論方法與數(shù)值模擬

2.1 距離冪次反比法原理

距離冪次反比法（Inverse Distance Weighted,IDW）是最常用的空間內(nèi)插方法之一，通過已知的鉆孔巖樣品位數(shù)據(jù)來計算空間任意一點的礦石品位，不同的品位數(shù)據(jù)對同一點的影響不同，與待計算點距離越近的品位數(shù)據(jù)對該點品位值的影響越大[1?2]。計算公式為：

式中，g是估計值；gi是第i個樣本；di是距離；p是距離的冪次，它的大小顯著影響著估值的結(jié)果。

受品位變化影響，距離冪次反比法的不同冪次對應不同的適用范圍和估值效果。當品位變化較穩(wěn)定時，可取冪次為2，一般適用于Cu、Fe等元素的估值；當品位變化較大時，可取冪次為3，一般適用于Au、Ag的估值。

2.2 實體模型的建立

2.2.1 地表模型構(gòu)建

依據(jù)地形圖等高線標注信息，利用DIMINE軟件對等高線賦值，然后利用賦好高程的等高線進行地形建模，得到的礦區(qū)地形地貌實體模型見圖1。

圖1 礦區(qū)地形地貌

2.2.2 礦體建模

本文采用非標準網(wǎng)格法進行礦體建模。該方法是將整體建模區(qū)域分成單個單元內(nèi)的建模。在每個單元格的建模后，把所有單元格模型組合在一起構(gòu)建生產(chǎn)礦體模型，如圖2所示。

圖2 ?375～?430m東部生產(chǎn)礦體模型

3 統(tǒng)計分析與儲量計算

3.1 地質(zhì)數(shù)據(jù)庫的建立與統(tǒng)計分析

3.1.1 建立地探數(shù)據(jù)庫

礦山前期已經(jīng)整理了地質(zhì)勘探、生產(chǎn)勘探以及深部勘探的地質(zhì)數(shù)據(jù)，根據(jù)數(shù)據(jù)特點按地質(zhì)數(shù)據(jù)庫建庫要求，將數(shù)據(jù)整理成開口表、測斜表、樣品表和巖性表，進行鉆孔校驗和生成。在孔口、測斜、樣品、巖性表4個文件經(jīng)邏輯校驗無誤后，為了保證地質(zhì)數(shù)據(jù)準確性及與原數(shù)據(jù)的一致性，對所有鉆孔進行空間檢查、分析和校正。

在鉆孔數(shù)據(jù)庫空間檢查、分析和校正完成后，完成了鉆孔數(shù)據(jù)庫的建立。

3.1.2 原始樣統(tǒng)計分析

胡咀-馬頭山水泥原料礦主要元素為CaO和MgO，為有效了解礦區(qū)品位分布，對礦床所有樣品的CaO和MgO品位進行統(tǒng)計分析，如圖3、圖4所示。其中，CaO品位最大值為55.68，品位平均值為49.76；MgO品位最大值為9.34，品位平均值為1.15。

3.2 地質(zhì)數(shù)據(jù)庫的樣長組合與統(tǒng)計分析

3.2.1 樣長組合

通過勘探數(shù)據(jù)庫中的勘探數(shù)據(jù)可以對塊段模型中單元區(qū)塊的各種參數(shù)進行估值，同時也可以計算礦床儲量。根據(jù)地質(zhì)統(tǒng)計學原理[3]，為了確保各種參數(shù)的無偏估計，所有樣本數(shù)據(jù)都應處于相同的負荷下，即同類型參數(shù)的地質(zhì)樣品段的長度應 相同。

圖3 全礦區(qū)CaO品位分布直方圖

圖4 全礦區(qū)MgO品位分布直方圖

運用地質(zhì)統(tǒng)計學原理，對過濾后的鉆孔數(shù)據(jù)庫樣長進行統(tǒng)計，可看出原始樣品的樣長平均值為5.72m，如圖5所示，因此組合樣長取5.7m。

圖5 地質(zhì)數(shù)據(jù)庫的樣長的直方圖統(tǒng)計分析直方圖

樣品組合的計算公式為：

式中，GC為組合樣參數(shù)值；Li為第i個樣品的長度；LC為組合樣的長度；m為參與組合樣計算的樣品數(shù)。經(jīng)過組合后的樣長數(shù)據(jù)庫如圖6所示。

圖6 樣長組合數(shù)據(jù)庫

3.2.2 組合樣統(tǒng)計分析

對“礦體內(nèi)部樣長組合數(shù)據(jù)庫”中的CaO和MgO含量進行統(tǒng)計分析，統(tǒng)計參數(shù)見表1，組合后CaO和MgO品位統(tǒng)計直方圖如圖7、圖8所示。

表1 CaO和MgO組合后品位基本統(tǒng)計參數(shù)

圖7 組合后CaO品位統(tǒng)計直方圖

圖8 組合后MgO品位統(tǒng)計直方圖

由圖7、圖8可以發(fā)現(xiàn)，CaO和MgO品位的平均值在組合前后基本一致。組合樣的品位數(shù)據(jù)是估計推算塊段模型和礦床儲量計算的基礎(chǔ)，因此，確定組合樣長是建模的重中之重。在組合樣計算時，需仔細確定組合樣長，根據(jù)理論和多次實驗得出，取原始組合樣長平均值是比較合理的。

3.3 塊段3D模型可視化與品位估計

塊段3D模型是礦山儲量計算及品位估計的基礎(chǔ)，其思路是按照一定的標準將礦床細分為大小一致的單元塊，然后根據(jù)已知品位的樣品單元塊，利用空間插值的方法對整個礦床的單元塊進行品位推估，然后在估值的基礎(chǔ)上統(tǒng)計儲量，為實現(xiàn)礦山資源儲量的動態(tài)管理采切開拓設(shè)計、開采評價等提供數(shù)據(jù)支持。

3.3.1 塊段三維模型的建立

三維塊體模型被廣泛應用于礦山儲量計算、品位估計以及露天礦山臺階爆破設(shè)計和開采境界優(yōu)化等。塊體的品位模型是將礦床完整地呈現(xiàn)在三維空間內(nèi)，按照一定的尺寸標準將其微分化，然后根據(jù)有限已知品位的樣品點對整個礦床的品位進行估值，同時，在此估值基礎(chǔ)上進行儲量的計算。

在DIMINE中打開計算礦山儲量的礦體文件，然后建立全包圍礦體的空塊段模型，如圖9所示。

圖9 生產(chǎn)礦體空塊段模型

靈活的約束可以實現(xiàn)工作中各種量的計算需要，比如可以進行采礦設(shè)計的礦巖量計算、貧化率計算和多種級別下礦量的統(tǒng)計等。距離冪估值、克立格估值、常量賦值、變量賦值、塊段法儲量計算等都需要用到塊段模型約束。在進行實體內(nèi)部約束之后的礦體塊段模型如圖10所示。

圖10 礦體約束后的塊段模型

3.3.2 品位推估

距離冪次反比法作為一種應用范圍比較廣的插值方法，其原理簡單易懂、計算結(jié)果準確可靠，本次研究品位空間插值采用此方法來完成[4]。

品位估值是利用已知樣品點對每一個單元塊進行品位的估算。根據(jù)樣品分布特點，估值主要分為3個部分：夾石估值、品位估值和估值驗證[5?6]。

（1）夾石估值。由于礦體內(nèi)部原始樣TFe品位為多峰分布，低品位礦石在礦床中占一定比例，在估值時首先需考慮低品位礦石的影響，即估出低品位礦石的影響范圍。夾石估值時用全部樣對夾石范圍進行估值。

（2）品位估值。在夾石類型定義完成后，以“礦體內(nèi)部+夾石外部”的約束，利用TFe≥20的樣品對每個礦體進行品位估值。估值時從“最小樣品搜索半徑”開始，不斷擴大搜索半徑直到所有的礦體都被估值。

（3）估值結(jié)果驗證。驗證塊段模型估值的正確性是一個非常重要的步驟，關(guān)系到后續(xù)儲量統(tǒng)計和資源科學管理，驗證估值結(jié)果一般有以下3種方法：

方法1：通過剖面對鉆孔化驗數(shù)據(jù)與塊段估值數(shù)據(jù)進行局部校驗。

方法2：運用塊段模型報告得出品位噸位曲線，通過實體模型體積乘比重得到的礦量與塊體模型的礦量比較，進而從全局上校驗塊體模型。

方法3：運用基本統(tǒng)計來比較塊段模型的平均品位和組合樣的平均品位。

本次采用方法3進行驗證，驗證對比結(jié)果見表2。

3.4 儲量計算

在品位估值完成后，利用DIMINE軟件分別完成了80m～206m生產(chǎn)臺階礦量的計算，結(jié)果見表3。

表2 礦區(qū)按中段估算資源品位統(tǒng)計

表3 礦區(qū)按中段估算資源儲量統(tǒng)計表

4 結(jié)論

采用大型三維礦業(yè)軟件DIMINE，先后完成了某礦區(qū)三維地表模型和三維地質(zhì)數(shù)據(jù)庫創(chuàng)建；基于地質(zhì)剖面和校正平面建立了生產(chǎn)礦體模型；運用距離冪次反比法對礦床CaO和MgO進行品位估值，并對估值數(shù)據(jù)按各分段分別進行了儲量計算和對比分析。經(jīng)過分析，所建模型和儲量統(tǒng)計結(jié)果均和傳統(tǒng)的地質(zhì)統(tǒng)計學分析結(jié)果相吻合。三維建模成果不僅實現(xiàn)了該礦山地表、礦體的三維可視化，還為礦山設(shè)計、計劃編排以及施工都提供了準確的數(shù)據(jù)支撐。