固體礦產資源量估算規(guī)程 第 3 部分：地質統(tǒng)計學法

DZ/TXXXXX—XXXX

固體礦產資源量估算規(guī)程第3部分：地質統(tǒng)計學法

1范圍

本部分規(guī)定了地質統(tǒng)計學法估算固體礦產資源量的基本原理、適用條件及方法選擇、估算參數選擇、

數據準備、估算流程、資源量分類及估算成果。

本部分適用于勘查程度達到詳查階段及以上要求的固體礦產勘查、礦山設計、開發(fā)階段的資源量估

算工作。

2規(guī)范性引用文件

下列文件對于本文件的應用是必不可少的。凡是注日期的引用文件，僅所注日期的版本適用于本文

件。凡是不注日期的引用文件，其最新版本（包括所有的修改單）適用于本文件。

GB/T13908固體礦產地質勘查規(guī)范總則

GB/T17766固體礦產資源儲量分類

GB/T33444固體礦產勘查工作規(guī)范

DZ/T0033固體礦產勘查/礦山閉坑地質報告編寫規(guī)范

3術語和定義

下列術語和定義適用于本文件。

3.1

礦化域Domain

礦化域是基于地質控礦因素（如巖性、地層、構造、圍巖蝕變等）或礦化邊界或將二者綜合來考慮

圈連的，一個比較連續(xù)的礦化地質體。

3.2

區(qū)域化變量RegionalizedVariable

區(qū)域化變量是一種限制在一個特定的空間內（如一個礦體內或礦化域內）的變量，既有結構性（相

關性）又有隨機性的特征。如在特定的空間內的樣品品位之間是相關的，可以參與礦體品位的估算。超

出這個空間范圍的樣品，由于樣品品位之間相關性減弱，隨機性的成分增大，一般不參與礦體品位的估

算。

3.3

變異函數Variogram

1

DZ/TXXXXX—XXXX

變異函數是地質統(tǒng)計學的應用工具，它能夠反映區(qū)域化變量的空間變化特征—相關性和隨機性，反

映局部范圍和特定方向上地質特征的變化。變異函數一般用變異曲線來表示（見圖1）。

圖1變異函數球狀模型曲線示意圖

3.4

實驗變異函數ExperimentalVariogram

變異函數對任一給定礦床Ω是未知的，需要通過取樣值對之進行估計。設從礦床Ω中獲得一組樣

品，相距h的樣品對數為n(h)，那么變異函數γ(h)可以用下式估計：

2

1nh

hXzXzh........................（1）

ii

2nhi1

式中：

γ(h)--實驗變異函數值；

X(zi)--在zi處的樣品值；

X(zi+h)--在與zi相距h處的樣品值。

公式中含有1/2,又稱實驗半變異函數。實驗變異函數示意如圖2所示。

圖2實驗變異函數示意圖

2

DZ/TXXXXX—XXXX

3.5

理論變異函數TheoreticalVariogram

實驗變異函數由一組離散點組成，將實驗變異函數擬合為一個可以用數學解析式表達的數學模型即

理論變異函數（標準模型），如圖3所示。常見的變異函數的數學模型有以：球狀模型、指數模型、高

斯模型、空穴效應模型等。

圖3理論半變異函數擬合示意圖

3.6

塊金效應NuggetEffect

由于取樣分析過程中的誤差和微觀礦化作用的變化，兩個樣品的分析結果不會完全相同。因此，變

異函數值在原點附近實際上不等于零，這種現象稱為塊金效應。當距離R=0時，變異函數值不為0，變

異函數曲線和縱軸相交，這個值稱為塊金值(Nugget)，用C0表示。它表示距離R很小時兩點間品位的變

化，C0代表隨機的成分。

3.7

基臺值Sill

是參與變異函數計算數據的方差，用C表示，它反映某區(qū)域化變量在研究范圍內變異的強度。變異

函數曲線達到基臺后開始收斂，在基臺附近波動。C+C0組成總基臺。應用中主要看基臺值C和塊金常數

C0所占的比例，C代表相關的成分。

3.8

變程Range

3

DZ/TXXXXX—XXXX

用來表達區(qū)域化變量空間相關性的最大距離即是當變異函數曲線最初達到基臺時所對應的橫坐標

值，用R表示。R反映了相關范圍，以待估塊為中心，R為半徑的范圍內的信息值對待估域進行估計。同

時，變程的大小反映了研究對象（如礦體）中某一區(qū)域化變量（如品位）的變化程度。

3.9

搜索橢球體SearchEllipsoid

搜索橢球體定義了對待估塊進行品位估算時，搜索的空間范圍和相關參數，用于品位估算時對樣品

的選擇。應用克里格法估算時應計算各個方向的變異函數，確定搜索橢球體主軸、次軸和短軸的方位并

根據變程求出軸比。

3.10

幾何異向性GeometricAnisotropy

在幾個方向上研究區(qū)域化變量時，當一種礦化現象在各個方向上性質不同時稱各向異性，根據變異

函數曲線，區(qū)域化變量在不同方向上表現出基臺值相同而變程不同時稱為幾何異向性。幾何異向性具有

相同的基臺值C而變程R不同（如圖4）。

圖4幾何異向性示意圖

3.11

帶狀異向性Zonalanisotropy

對應幾何異向性，根據變異函數曲線，區(qū)域化變量在不同方向上表現出基臺值不同，而不論變程相

同或不相同時稱為帶狀異向性（如圖5）。

4

DZ/TXXXXX—XXXX

圖5帶狀異向性示意圖

3.12

克里格法Kriging

是地質統(tǒng)計學中一種局部估計方法，以變異函數為主要工具，對區(qū)域化變量進行插值，求插值過程

中的最優(yōu)線性無偏估計量，再通過塊體約束計算資源量的方法（通常叫克里格法）。根據研究目的和條

件不同，分別有簡單克里格法（SimpleKriging）、普通克里格法（OrdinaryKriging）、對數正態(tài)克里

格法（LognormalOrdinaryKriging）和指示克里格法（IndicatorKriging）和泛克里格法（KrigingwithTrend）

等，其中應用最為廣泛的是普通克里格法。

3.13

距離冪次反比法InverseDistanceWeight

距離冪次反比法是一種與空間距離有關的插值方法,即在估計待估點的值時,按照距離越近權重值

越大的原則,利用已知點和待估點之間距離取冪次后的倒數為權重系數進行加權平均，再通過塊體約束

計算資源量的方法。

3.14

克里格方差EstimationVariance

采用克里格法估計時，估計誤差的方差稱為克里格方差，是衡量估計精度的度量。估計方差小證明

估計精度高。

4基本原理

4.1克里格法原理

5

DZ/TXXXXX—XXXX

應用克里格法對區(qū)域化變量進行局部估計時，將礦體劃分成許多相同或相似的長方體V，用在一定

范圍內的系列樣品值vα（α=1，2，…，n）為V估值，以估計方差最小、權系數之和等于1為條件，形

成克里格方程組：

n

Cv,vCv,V

1

n..........................（2）

1

1

式中：

Cv,v

--樣品對之間的協方差平均值；

Cv,V--樣品點與長方體中心之間的協方差平均值；

--樣品的權系數；

μ--拉格朗日因子。

n

以求解克里格方程組所得的權系數，采用加權平均的方法得到的值，即：vV。

V

1

4.2距離冪次反比法原理

距離冪次反比法是利用樣品點和待估塊中心之間距離取冪次后的倒數為權系數進行加權平均的方

法。以待估塊中心為圓心、以搜索半徑做圓（或橢圓），計算落入圓（或橢圓）內每一樣品與待估塊中

心的距離。計算公式：

nxn1

xi................................（3）

bmm

i1dii1di

式中：

xi--落入影響范圍的第i個樣品的值；

di--第i個樣品到待估塊中心的距離；

m--冪次；

xb--待估塊的估值。

5適用條件及方法選擇

5.1地質工作程度達到詳查及以上，樣品數據量和密度較大，足以計算出各個方向實驗變異函數，并

能通過擬合轉化為理論變異函數時采用克里格法。

5.2礦體形態(tài)為二維延展時，采用二維克里格法；礦體形態(tài)為三維延展時，采用三維克里格法。

5.3根據礦體結構：礦石自然類型、工業(yè)類型、工業(yè)品級和非礦夾石的形態(tài)、空間分布特征、種類和

它們的相互關系以及成礦期的構造，選擇相應的克里格法。

5.4根據區(qū)域化變量統(tǒng)計分析，結合變異函數的分布特征分別選擇如下方法：

a)當假設區(qū)域化變量的數學期望值為某一已知常數時，選用簡單克里格法;

b)當區(qū)域化變量滿足二階平穩(wěn)假設（或內蘊）且數學期望值未知時，選用普通克里格法;

6

DZ/TXXXXX—XXXX

c)當區(qū)域化變量不服從正態(tài)分布但服從對數正態(tài)分布時，選用對數正態(tài)克里格法;

d)對有多期礦化、區(qū)域化變量呈現多峰分布時，選用指示克里格法；

e)當區(qū)域化變量在大范圍內具有方向性的趨勢變化（即非平穩(wěn)性），或具有漂移存在時，選用泛

克里格法。

5.5當區(qū)域化變量變化系數滿足礦化在比較均勻的范圍內（對應礦化均勻或比較均勻的具體變化系數

值可查各分礦種規(guī)范，一般變化系數≤150%左右）可用距離冪次反比法或普通克里格法。

6估算參數選擇

6.1克里格法估算參數包括變異函數參數、搜索橢球體參數和塊體參數；距離冪次反比法估算參數包

括冪次、搜索橢球體參數和塊體參數。

6.2變異函數模型參數包括最佳基本滯后距、塊金值、基臺值、變程。

a)基本滯后距是計算變異函數時分隔樣品對的矢量的最小長度，一般由計算全向變異函數獲得。

當選定一個基本滯后距得到的全向實驗變異函數最容易擬合出理論變異函數時，即為最佳滯后

距。

b)塊金值反映出區(qū)域化變量在小尺度上的變異程度，表明變量的隨機成分。沿鉆孔方向以最小滯

后距計算變異函數時，所擬合的理論變異函數曲線與縱坐標的交點即為塊金值，見圖1。一般

情況下，一個估算域內只有一個塊金值。

c)基臺值表明了估算域內變量的相關程度，其數值由變量的方差減去塊金值取得。構成了變異函

數的總基臺，見圖1。

d)變程R表明了變量在估算域內自相關性存在的最大范圍，一般是理論變異函數曲線中最初達到

總基臺值時所對應的距離，見圖1。

6.3距離冪次反比法冪次參數選?。?/p>

a)冪次一般在1～3之間取值。變量的空間變化性越大、變化速度越快，則取值越大。

b)根據工程網度和樣品密度，密度大，取大值，反之取小值。

c)應通過交叉驗證冪次的合理性。

6.4搜索橢球體參數包括橢球體的方位角、傾伏角、傾角，主軸的搜索半徑、次軸的搜索半徑（或用

主軸與次軸的比值表示）、短軸的搜索半徑（或用主軸與短軸的比值表示）。若樣品點空間分布不均勻，

具有叢聚的情況，可將搜索橢球體分為多個扇區(qū)，每個扇區(qū)內可根據需要確定最少工程數、最少樣品數

和最多樣品數。

6.5常用的克里格方法及距離冪次反比法估算參數均包括橢球體參數、最多樣品數、最少樣品數、最

少工程數，變異函數參數只在各克里格方法中涉及，具體見參數匯總表（見表1）。

表1不同估算方法參數表

估值參數

方法

橢球體參數最多樣品數最少樣品數最少工程數變異函數參數附加特征參數

簡單克里格法√√√√√均值參數

普通克里格法√√√√√

指示克里格法√√√√√IK指示值

對數克里格法√√√√√對數變換參數

泛克里格法√√√√√漂移參數

距離冪次反比法√√√√冪次

7

DZ/TXXXXX—XXXX

6.6塊模型（也叫品位模型）的參數主要是設置塊模型方向和塊體尺寸（長、寬、高）以及進行次級

分塊，如圖6。塊模型屬性可記錄塊體所在的位置（即塊質心點坐標）、尺寸大小，還可以附加不同屬

性如礦巖類型、礦種、品位、體積質量、數學運算賦值的屬性、資源類別等。這些屬性均保存在塊模型

文件中，既可以通過估值的方法賦值，也可以通過屬性名稱直接賦值，還可以進行相關數學運算賦值。

圖6塊模型示意圖

7數據準備

7.1建立數據庫，從地質數據庫或電子表格中提取參與資源量估算的基礎數據，至少包括勘查工程或

生產工程定位數據表（見表2）、測斜數據表（見表3）、樣品分析數據表（見表4）和地質巖性表（見

表5）。表中所列的字段是必不可少的。

表2工程定位數據表

鉆孔名稱東坐標（E）北坐標（N）高程（R）終孔深度（m）工程類別

表3測斜數據表

鉆孔名稱井深（m）方位角(度)傾角（度）

表4樣品分析數據表

鉆孔名稱樣品編號從（m）到（m）樣長（m）品位(分列)

表5地質巖性表

鉆孔名稱從（m）到（m）巖性代碼巖性描述

7.2對數據進行糾錯及完整性和邏輯性檢查。

7.3通過三維軟件對工程數據進行位置對比查看和空間關系分析。

8

DZ/TXXXXX—XXXX

8估算流程

8.1地質解譯建立礦化域或礦體模型

8.1.1將地質數據庫中的鉆孔、探槽和坑道等工程數據繪制成圖，作出剖面圖或中段平面圖，參見附

錄A。礦巖界線清晰時，在剖面圖或中段平面圖上根據品位結合地層、巖性、構造、以及礦體產狀等地

質特征進行地質解譯，圈定礦體或礦化域及夾石邊界線。

8.1.2根據礦化特征，礦巖界線呈過渡關系時，選取礦化品位值或低于一般工業(yè)指標圈連礦化域，采

用礦塊指標體系圈定礦體范圍，參見附錄B。

8.1.3剖面圖上礦體或礦化域邊界線的圈定和外推原則，參見本規(guī)程通則部分和有關的固體礦產勘查

工作規(guī)范。

8.1.4進行資源量估算通常做的模型包括：地形三維模型、礦體或礦化域模型、夾石模型、巖體模型

和地質構造模型等。

8.1.5采用礦體模型或礦化域模型作為估算范圍時，需要對模型進行質量檢查，確保模型的合理性。

8.2估算域的劃分

8.2.1當礦體被斷層、巖脈等線性構造切割或錯斷見圖7，使得構造兩側區(qū)域化變量統(tǒng)計特征存在明

顯差異，應以此線性構造為界劃分不同的區(qū)域、設置不同的估算參數進行估算。

圖7礦體被斷層錯斷示意圖

8.2.2當礦體出現褶皺，使得礦體的走向和傾角發(fā)生轉折見圖8，無法使用一個搜索橢球體，應將礦

體劃分成不同的估算域，分別設置不同的估算參數進行估算，或用動態(tài)橢球參數進行估值。

圖8礦體的產狀受褶皺影響發(fā)生改變示意圖

8.2.3當礦體內出現多個礦石類型或者有用元素局部富集見圖9、或者受巖性控制圖10，而使區(qū)域化

變量統(tǒng)計特征存在明顯差異，應將礦體按不同礦石類型或富集特征劃分成不同的區(qū)域，分別設置不同的

估算參數進行估算。

9

DZ/TXXXXX—XXXX

貧礦體

富礦體

圖9根據礦體的富集特征不同劃分估算域示意圖

圖10根據巖性不同劃分估算域示意圖

8.3樣品數據統(tǒng)計分析及特異值處理

8.3.1對樣品統(tǒng)計分析及特異值處理時，應以礦化域或礦體為單位進行。

8.3.2對區(qū)域化變量進行估值前，分別繪制區(qū)域化變量的分布直方圖、累計頻率分布曲線或概率圖等，

判斷其分布特征，并計算區(qū)域化變量的統(tǒng)計特征，如最小和最大值、樣品數、平均值、中值、方差、標

準差、偏度、峰度和變化系數等。

8.3.3識別和處理特異值可采用數理統(tǒng)計的方法：分位數法、估計鄰域法、影響系數法、概率曲線和

累計頻率分布曲線法等。特異值處理后的檢驗參見附錄C。

8.3.4如采用對數正態(tài)克里格法或指示克里格法進行變量估值時，不進行特異值的處理。

8.4樣品等長度組合

8.4.1樣品等長度加權組合應限制在相應的礦體（層）內。其長度取值應為統(tǒng)計中主要的樣長值。

8.4.2組合后的樣段與原樣段的長度變量統(tǒng)計特征應保持最大程度上的一致。

8.4.3等長組合須在對原樣品進行特異值處理后進行。

8.5建立塊模型

8.5.1根據礦化域或礦體空間范圍、形態(tài)和產狀，采用一定長、寬、高尺寸的塊為單位，形成整個估

算區(qū)域的塊模型（可根據礦體產狀和形態(tài)創(chuàng)建次級分塊和旋轉塊模型）。

8.5.2影響塊模型尺寸的主要因素是礦體規(guī)模、勘查工程間距、開采方式（露采或地下開采）和采礦

工藝（礦房尺寸或露天開采臺階高度）、克里格效率等。其中克里格效率=（塊方差-克里格方差）/塊

方差，其值介于0和1之間。通過比對不同尺寸的塊模型的克里格效率值，克里格效率值最大時所對應

的塊模型尺寸是最優(yōu)的。

8.6建立變異函數模型

8.6.1確定最佳基本滯后距。

10

DZ/TXXXXX—XXXX

8.6.2確定塊金值。

8.6.3通過計算各方向的變異函數，確定估算域的最大連續(xù)性方向（主軸）、次連續(xù)性方向（次軸）

和最小連續(xù)性方向（短軸），三個方向在空間上應相互垂直。求得三個方向的實驗變異函數（見圖11），

采用適當的數學模型進行擬合，求出每一個方向的基臺值C、變程值R。

圖11三個方向的實驗變異函數示意圖

8.6.4變異函數結構的套合：

a)進行變異函數結構套合的目的是把不同距離上的變異性組合起來代表整個礦體的變異結構。每

個方向可以由不同的數學模型套合形成一個理論變異函數，但不同方向上數學模型的個數和類

型須一致。

b)進行估算前，需要將搜索橢球體三個軸向上的理論變異函數套合在一起。當三個軸的變異函數

結構具有相同的塊金值和基臺值而具有不同的變程值（見圖12、圖13和圖14），具有幾何異

向性的特征時，變異函數的套合方法是用搜索橢球體的主軸方向的數學模型來代表其他方向，

根據搜索橢球體三個軸向的變程值，再給出主軸與次軸及主軸與短軸的比值。

圖12主軸方向變異函數示意圖

11

DZ/TXXXXX—XXXX

圖13次軸方向變異函數示意圖

圖14短軸方向變異函數示意圖

8.7搜索橢球體的設置

8.7.1采用克里格法時，通過變異函數特征，確定搜索橢球體各軸向的產狀和半徑。

8.7.2確定搜索橢球體的產狀后，應將搜索橢球體與估算域或礦體進行疊加對比，觀察兩者的產狀是

否協調一致如圖15所示。

圖15搜索橢球體與礦體或估算域疊加示意圖

12

DZ/TXXXXX—XXXX

8.7.3確定搜索橢球體扇區(qū)數，并設置每個扇區(qū)的最少和最多樣品數，減少叢聚數據對估值的影響。

8.7.4采用距離冪次反比法時，搜索橢球體參數應依據礦體規(guī)模和產狀確定，還要結合勘查類型和工

程間距，也可以通過變異函數模型確定。

8.7.5采用可變橢球體進行估值，應使用實體模型對搜索橢球體參數進行約束。

8.8交叉驗證

8.8.1交叉驗證的目的是檢驗變異函數模型參數和搜索橢球體參數的合理性。

8.8.2交叉驗證的理想結果是真值與估計值的誤差均值趨近于“0”，誤差方差與克里格估計方差比值

趨近于“1”。

8.8.3若交叉驗證結果不理想，應修改或重新擬合理論變異函數并再次交叉驗證，直到取得一套合理

的變異函數參數和搜索橢球體參數為止。

8.9變量估值

8.9.1對每一個礦化域或礦體，可根據區(qū)域化變量的分布和統(tǒng)計特征，選擇普通克里格法或距離冪次

反比法。其他的克里格法如對數正態(tài)克里格法和指示克里格法的實施步驟參見附錄D。

8.9.2采用普通克里格法估算時，對個別礦體規(guī)模較小且計算變異函數有困難的情況下，可采用距離

冪次反比法或其他方法進行估算。

8.9.3使用距離冪次反比法進行變量估值時，應根據工程間距、樣品數和估算域規(guī)模確定估值搜索半

徑；使用克里格法時，估值搜索半徑的選擇還須主要考慮變程值。

8.9.4對于采用雙指標體系圈定的礦體，最大搜索半徑可超過變程，確保礦體內所有塊體都能獲得相

應的估計值。采用礦塊指標體系時，估值的最大搜索半徑一般不超過變程值的4倍。

8.10變量估值結果驗證

8.10.1采用地質統(tǒng)計學法（含距離冪次反比法）對區(qū)域化變量進行估計后，應采用局部驗證或全局驗

證的方法進行檢驗。

8.10.2局部驗證的方法是在剖面圖或平面圖上比較參與估算的組合樣品值與相鄰塊體估計值之間的

相近程度、外推的距離和方向，判斷估值的可靠性。

8.10.3全局驗證的方法是沿某些方向按照一定的間距劃分一系列的區(qū)域，分別計算落在每一個區(qū)域中

參與估值的組合樣品平均值和塊模型估值的平均值，繪制成兩條曲線圖，通過兩條曲線的吻合度，判斷

估值的合理性和可靠性。

8.10.4可采用另一種適用的估算方法對主要礦體進行驗算，對比絕對誤差和相對誤差，判斷估值參數

的可靠性。

9資源量分類

9.1應根據搜索半徑、工程數、樣品數等因素對塊體進行資源量分類。對于“孤島”現象，即不同類

別在空間上相互穿插或包含，應根據工程控制程度和地質可靠程度進行人工干預。

9.2應在剖面圖上或平面圖上勾畫出不同分類輪廓界線并圈連成分類輪廓體，對塊模型的分類屬性進

行賦值。

9.3可采用國內的資源量分類規(guī)范進行分類。

10資源量估算成果

13

DZ/TXXXXX—XXXX

10.1資源量匯總參照本規(guī)程通則及資料性附錄E和國家相關規(guī)范執(zhí)行。

10.2地質統(tǒng)計學法提交資源量報告章節(jié)的編寫參見相關報告編寫規(guī)范。

14

DZ/TXXXXX—XXXX

AA

附錄A

（資料性附錄）

三維圖件的編制

A.1資源量計算的基本圖件

A.1.1按照DZ/T0078等規(guī)范要求，一般先制作地質勘查所需要的礦區(qū)（床）地形地質圖、勘查線剖面

圖、水平斷面圖、中段地質圖等。

A.1.2根據GB/T33444、DZ/T0079等規(guī)范要求，制作與資源量估算有關的礦體（層）頂底板等高線圖

（或含資源量估算圖）、資源量估算圖（水平或垂直縱投影圖）、砂礦或緩傾斜礦體（層）需要的頂底板

等高線和礦層等厚線圖（或含資源量估算圖）等。

A.2三維斷面圖件的制作

A.2.1在實際工作中，通常是按照二維平面圖件的方式完成相應的斷面圖，包括平面圖和剖面圖。這

些圖件的坐標均為標注值而不是實際的位置，比如，水平斷面圖中所有的標高均為0。在資源量估算中，

如果采用三維軟件按照實際剖面位置和工程坐標及測斜數據進行基本圖元（工程、礦段、構造等）的制

作，更有利于礦體的對應分析、塊段劃分和資源量估算。

A.2.2制作時，可以按照勘查線或水平中段分別進行切制，形成相應的斷面圖位置。再利用三維軟件

中工程的位置、軌跡和礦層信息，在一定的斷面范圍內，按照實際的位置保留工程的位置并形成點或者

線文件進行保存，從而形成了一個斷面、一定范圍內的三維斷面圖件的所有制圖元素（如圖A.1、A.2）。

圖A.1實際位置剖面圖

15

DZ/TXXXXX—XXXX

圖A.2實際高程的中段圖

A.3三維縱投影圖的制作

A.3.1礦體縱投影圖分為水平投影圖和垂直投影圖。水平投影圖可以投影到任意平面。但垂直縱投影

圖是可以實現與礦體走向相一致的三維投影圖（如圖A3）。

A.3.2制作縱投影圖的步驟：一是按照基線（垂直勘查線）方向作縱剖面，其走向與礦體走向一致，

但與礦體有一定的間距；其次將地形線、礦巖界線、探礦工程點、構造界線、塊段界線等投影到該剖面

上，并進行塊段劃分和資源量估算，同時需要保存這些文件。

圖A.3三維縱投影示意圖

16

DZ/TXXXXX—XXXX

BB

附錄B

（資料性附錄）

礦塊指標體系及圈礦原則

B.1礦塊指標體系一般根據地質礦化規(guī)律采用某一個品位（一般介于地質上的礦化品位與常用指標體

系中的邊界品位之間）圈出的一個比較完整的礦化域，在礦化域內按照一定的大小劃分估計品位的單元

塊，繼而對單元塊進行品位估值，再采用邊際品位界定單元塊是礦石還是廢石，然后統(tǒng)計資源量。

B.2圈礦原則

采用地質統(tǒng)計學法（含距離冪次反比法）進行礦化域（礦體）圈定應遵循的原則：

a)根據礦化特征（地層、巖性、構造、圍巖蝕變、礦石類型等），當礦巖界線呈過渡關系時，可

采用礦塊指標體系圈定礦化域。

b)礦化域的品位邊界可從樣品分布直方圖、累計分布曲線、概率圖中尋找礦化拐點確定，或通過

類比法確定，或采用邊界品位的1/2或根據尾礦品位值確定。礦化域品位邊界值，通常等于或

小于常用（雙）指標體系中邊界品位，其目的使所圈定的礦化域更完整連續(xù)。

c)在圈定礦化域時能剔除的夾石盡量剔除，如不剔除夾石，將對礦體的平均品位產生貧化，影響

品位的正態(tài)分布。

d)完成估算后，在塊模型中由大于或等于邊際品位的并且連續(xù)的、滿足開采條件的單元塊形成工

業(yè)礦體，然后統(tǒng)計資源量。

B.3當品位分布均勻，礦體與圍巖又是漸變關系時，采用礦塊指標體系的優(yōu)點更突出；相反條件下，

采用常用（雙）指標體系的有著顯著的優(yōu)點。

17

DZ/TXXXXX—XXXX

附錄C

（資料性附錄）

特異值處理合理性檢驗

C.1在大多數的地質問題中，礦體的厚度或礦石的品位（地質變量x）不服從正態(tài)分布，但其對數值（即

y=ln(x)）趨于正態(tài)分布，則稱為地質變量服從對數正態(tài)分布。

C.2若地質變量不服從對數正態(tài)分布，但y=ln(x+a)服從正態(tài)分布（其中a是一個常數）,則稱地質變量

服從三參數對數正態(tài)分布。

C.3服從對數正態(tài)分布或三參數對數正態(tài)分布的樣品一般存在特異值，對此地質統(tǒng)計學有相應的處理

方法。特異值處理的合理性可用西舍爾T估值檢驗（Sichel'sT）。

C.4西舍爾T估值計算公式如下:

*2

*lnZV

Zve（ne）................................（4）

其中系數：

22

12n14（n1）6

（）1

ne2e222!n1e223!(n1)(n3)e

式中：

*

Zv--西舍爾T估值；

*

elnZv--樣品的幾何平均值；

2

n(e)--西舍爾系數，為一個臺勞級數，一般取前三項即可滿足精度要求；

2

e--對數變換后的樣品方差；

n--樣品個數。

C.5西舍爾估值檢驗

西舍爾估值的檢驗分為如下幾個步驟：

a)首先對已經處理特異值的樣品值取自然對數；

b)計算取自然對數后的樣品品位平均值和方差；

c)將對數平均值返回真值得到幾何平均值；

d)根據對數變換后的品位方差和樣品個數，計算西舍爾系數；

e)將幾何平均值乘以西舍爾系數得到西舍爾T估值。

f)當特異值處理后的樣品的算術平均值小于并和西舍爾T估值接近時，判斷特異值處理結果合

理。

18