張 敏

（新疆維吾爾自治區(qū)地質礦產(chǎn)勘查開發(fā)局第十一地質大隊，新疆昌吉831100）

隨著計算技術的發(fā)展，出現(xiàn)以數(shù)學統(tǒng)計為核心的地質統(tǒng)計學儲量計算方法，并逐漸發(fā)展成為主流的儲量計算方法，以手工計算為主的傳統(tǒng)方法逐漸被地質統(tǒng)計方法所替代[1]。儲量計算也從二維向三維逐漸發(fā)展，將礦床的三維建模與儲量計算結合起來是當前礦產(chǎn)資源儲量計算新的手段，它能夠更加直觀、準確地對礦產(chǎn)資源儲量進行計算。本次選取新疆備戰(zhàn)—察漢烏蘇鐵礦L1 和L2 號礦體進行幾何法與地質統(tǒng)計學法資源量估算對比分析研究。

1 幾何法資源量估算

約21 世紀20 年代以前，國內地勘單位主要以MAPGIS6.7 為主作為地質信息數(shù)字化的平臺，前人已借助MAPGIS6.7平臺采用幾何法對備戰(zhàn)—察漢烏蘇鐵礦進行了人工儲量估算。本次在國內主流三維礦業(yè)軟件3Dmine 中用地質塊段法進行資源量估算試驗，完成礦段圈定、礦體投影、創(chuàng)建塊段及儲量報告等操作后，得出軟件與人工儲量計算的相對誤差為4.6%，相對誤差較小，該方法本身不存在誤差，推測軟件與人工計算引起誤差的原因是人為因素或投影產(chǎn)生的誤差等。利用3Dmine 軟件進行傳統(tǒng)的地質塊段法儲量計算較人工方法計算更加方便快捷，在數(shù)據(jù)量較大時更顯優(yōu)勢，同時能夠減少人為因素產(chǎn)生的誤差，提高技術人員的工作效率，便于后期的數(shù)據(jù)管理及應用等。

2 地質統(tǒng)計學法資源量估算

本次分析對比在三維礦體建模的基礎上運用地質統(tǒng)計學普通克里格法進行資源量估算?？死锔穹ㄊ窃诰C合考慮了樣品的大小、形狀、其待估區(qū)域相互之間的空間分布以及區(qū)域變量的空間結構信息后，為了達到最優(yōu)的、無偏的估計方差，對每個樣品設定相應的權重值，利用鄰近鉆孔或坑道的樣品的品位來對這些樣品中間的某個待估區(qū)域的品位進行估算的方法。同時基于地質統(tǒng)計學的儲量估算方法，考慮了礦石的空間品位變化，可提高資源儲量估算的準確性[2]。具體如下所述。

2.1 實驗半變差函數(shù)分析

地質統(tǒng)計學中變差函數(shù)是最基本與最重要的模擬工具，是計算區(qū)域化變量空間結構特征的函數(shù)[5]。通過計算變量得到實驗變差函數(shù)，擬合出理論變差函數(shù)，再將不同方向的擬合結果套合為一個函數(shù)，用于估值計算。半變差函數(shù)計算公式為：

式中：r(h)——實驗變差函數(shù)；h——步長；

n(h)——步長為h的樣品對數(shù)；

z(x+h)、z(x)——樣品數(shù)據(jù)值[3]。

礦化一般是具有一定的方向性的，備戰(zhàn)—察漢烏蘇鐵礦表現(xiàn)出空間上的各向異性（圖1），需要針對不同方向計算和分析其各自的半變差函數(shù)。根據(jù)礦體產(chǎn)狀計算三個方向的半變差函數(shù)，走向、傾向和厚度方向分別對應主軸、次軸和短軸，用組合樣長后的數(shù)據(jù)計算三個方向的半變差函數(shù)（圖2），最終選擇擬合理論變差函數(shù)的參數(shù)：塊金值、基臺值及變程。本次分析使用指數(shù)模型。

圖1 鐵礦結構各向異性橢球

通過對實驗半變差函數(shù)和所擬合的理論半變差函數(shù)曲線圖進行分析，可以看出沿礦體厚度方向的半變差函數(shù)表現(xiàn)出一定周期性的跳動，說明礦體厚度方向元素的富集不均勻，礦體中有夾石存在。沿礦體傾向方向半變差函數(shù)的變程最大，說明礦體在傾向上有較好的連續(xù)性。Fe 元素的變化表現(xiàn)為各向異性，與實際情況相符。

2.2 建立空塊模型

通過建立空塊模型對礦體模型內部進行塊體分割，形成由n個體元模擬的區(qū)域，每個體元都有相應坐標并可以任意添加屬性，克里格法估值結果作為一個屬性值加入塊模型（圖3）。塊體單元的尺寸根據(jù)研究區(qū)地質復雜程度、礦體連續(xù)程度、勘探工程間距、研究尺度、計算機軟硬件條件等綜合確定。地質情況較復雜的區(qū)域塊體應盡量小，以保持實體模型邊界的精確度。

圖3 備戰(zhàn)—察漢烏蘇鐵礦空塊模型

2.3 交叉驗證

交叉驗證對所擬合的理論半變差函數(shù)的可靠性進行檢驗，其參數(shù)的可靠程度會直接影響礦床品位估算或儲量計算結果的精準性和可靠性，本次對比分析在估值計算前進行了交叉驗證。即從已知信息樣本總體中抽離某個點，用抽離點以外的其他信息結合半變差函數(shù)模型估計被抽離點，結果得到一系列估計值和誤差等參數(shù)，根據(jù)誤差判斷模型的準確性。對交叉驗證結果（表1）進行分析，誤差統(tǒng)計的均值為0.1377，誤差方差均值/克里格方差均值=56.3465/55.3810=1.017 接近于1。說明模型參數(shù)較可靠，可以用于后續(xù)計算。

表1 理論變異函數(shù)模型交叉驗證結果表

2.4 克里格估值

在塊體模型基礎上，設置估值參數(shù)、搜索參數(shù)、分區(qū)參數(shù)、孔約束等，用普通克里格法進行空間估值，在估值過程中通過改變主軸搜索半徑、分區(qū)設置、約束條件等參數(shù)進行了多次克里格估值（表2），確保所有的空塊均被賦予品位等屬性。

表2 普通克里格法估值參數(shù)表

3 資源量估算結果對比分析

對普通克里格法估值的塊體模型進行資源量報告，并與人工地質塊段法資源量計算結果進行比較，得出人工地質塊段法與普通克里格法兩種方法L1、L2礦段估算結果（表3）的相對誤差分別為0.73%、4.40%，相對誤差均小于5%，資源量估算結果的可靠性較強，同時提高了儲量估算的準確性。

表3 資源量估算結果對比表

兩種估算方法均有產(chǎn)生人為誤差的可能，如地質塊段法在制圖、投影、數(shù)據(jù)統(tǒng)計處理、計算等各方面產(chǎn)生的誤差，普通克里格法在三維建模、各類參數(shù)的選取等方面產(chǎn)生的人為誤差；兩種方法在體重值的選擇上可能產(chǎn)生誤差，地質塊段法在計算時不同塊段采用的體重值不同，而普通克里格法在儲量報告時采用的是平均體重值。誤差產(chǎn)生的原因是多方面、綜合性的，但兩種估算方法的相對誤差較小，估算結果可靠。

4 結論

（1）根據(jù)國內廣泛采用幾何方法進行資源儲量估算、提交報告的實際情況，以地質塊段法為例進行軟件和人工資源量計算的對比，結果顯示兩種方法的相對誤差小，且運用三維軟件計算資源儲量避免了人工計算繁瑣、人為因素多的缺陷，大大提高了技術人員的工作效率，并在實際生產(chǎn)中快速給出資源儲量報告。

（2）地質統(tǒng)計學普通克里格法能夠充分利用樣品信息，與幾何法相比較，在地質工作程度較高的詳查或勘探階段采用地質統(tǒng)計學法進行資源量估算具有明顯優(yōu)勢，通過三維地質建模實現(xiàn)了礦產(chǎn)資源儲量估算的三維可視化，促進了地質勘查工作由二維向二三維一體化轉變，同時使資源儲量估算成果和提交的報告便于與國際接軌。