賀耀文

（金川集團股份有限公司，甘肅 金昌 737100）

1 引言

隨著開采資源量的加大和找礦難度的增大，越來越顯示出礦產(chǎn)資源的預測和評價的重要性。通過對礦產(chǎn)資源的預測和評價，可提出礦山遠景，為今后找礦工作提供動力和指導建議。目前所采用的預測方法較多，諸如德爾菲法、概率法、邏輯信息法、回歸分析法、估算法、齊波夫定律等。其中齊波夫定律的地質(zhì)意義早在1975年就為N. J. Rowlands和D. Sampey所發(fā)現(xiàn)，他們論證了在相似成礦地質(zhì)背景下，同一礦區(qū)各個礦體的地質(zhì)資源儲量預測適用于齊波夫定律［1］，后來，國內(nèi)外許多專家學者也運用齊波夫定律進行了找礦潛力預測，并取得了良好效果。

本文試用齊波夫定律對甘肅省金昌市白家嘴子銅鎳礦床找礦潛力進行初步預測，探討齊波夫定律對白家嘴子礦區(qū)資源預測中應用的可能性，同時其預測結(jié)果將為白家嘴子礦區(qū)提供有用的找礦信息。

2 齊波夫定律

齊波夫定律是美國哈佛大學教授齊波夫于1949年提出的離散型概率分布數(shù)學模型。后來，人們發(fā)現(xiàn)這種分布規(guī)律廣泛存在于自然科學和社會科學中，且該方法簡單、有效，很快在世界各地被推廣使用［2］。

齊波夫定律的含義是：“如果有一組隨機數(shù)，將其從大到小排序后，如果最大數(shù)是次大數(shù)的2倍，是第3大數(shù)的3倍，……是第N 大數(shù)的N倍”［3］，則數(shù)組服從于齊波夫分布律。”換句話說，齊波夫定律就是指這樣的一個數(shù)列，其倒數(shù)比等于自然數(shù)列，這是帕雷托分布的極限情況［4］。數(shù)學表達式為：

式中：F為研究對象取值（礦體鎳金屬量）；R為等級值或秩，一般用自然數(shù)1，2，n表示；K為常數(shù)。齊波夫定律就是：

由（2）式可見，當F1=K，即秩R=1時研究對象取值等于K，所以我們只要求得最大值F1或K值，則其他各級的值將分別為：K/2，K/3，…，K/n。

首先，通過收集前期地質(zhì)勘查資料，可以統(tǒng)計出礦區(qū)內(nèi)已知礦體的資源量，計算出各個已知礦體的秩次和齊波夫預測秩，然后，計算齊波夫常數(shù)，即求出礦區(qū)內(nèi)潛在的最大礦體的資源量，最后，計算出礦區(qū)的總資源量（Q），即Q=F1+F2+……Fn，用資源總量減已知資源量就可以知道礦區(qū)內(nèi)潛在資源量。

3 白家嘴子礦區(qū)礦體地質(zhì)特征及資源現(xiàn)狀

白家嘴子礦區(qū)巖體長約6500m，寬20～500m，延深數(shù)百米至千余米。巖體東西兩端被第四系覆蓋，中部露出地表，上部已遭剝蝕，揭去覆蓋，巖體基巖面積約1.34km2。巖體走向310°、傾向南西,傾角50～80°，巖體被北東東向壓扭性斷層錯段，從西向東分為4段，依次編號為Ⅲ、Ⅰ、Ⅱ、Ⅳ個礦區(qū)（見圖1）。

圖1 白家嘴子銅鎳礦礦區(qū)示意圖

白家嘴子礦區(qū)整合采礦權范圍內(nèi)，主要礦體為Ⅱ-1、Ⅱ-2、Ⅰ-24、Ⅳ-1、Ⅲ-1號礦體，這5大主礦體地質(zhì)勘探時期的鎳金屬量累計約546萬t（表1），占全礦區(qū)的90.55%［5］。目前，開發(fā)該礦的礦業(yè)公司保有礦產(chǎn)資源儲量急劇減少，礦區(qū)深邊部地質(zhì)勘查工作也未取得較大成果。隨著該礦業(yè)公司的不斷發(fā)展，對銅鎳資源的需求將不斷增長，礦產(chǎn)矛盾日趨凸現(xiàn)，資源危機迫在眉睫。面對這一現(xiàn)狀，須通過地質(zhì)找礦，力爭探明新的銅、鎳資源量，以緩解該礦業(yè)公司快速發(fā)展帶來的保有資源銳減的矛盾，為該礦業(yè)公司長盛不衰發(fā)展提供資源保障。

4 白家嘴子礦區(qū)找礦潛力預測

地質(zhì)勘查資料是對礦區(qū)中各礦體的各種地質(zhì)特性觀測的結(jié)果。齊波夫定律認為礦體是觀測對象的函數(shù)。人們對礦床（體）的認識深淺程度，反映了觀察者的認識和理解能力及其思維活動能力的特性。一個礦床（體）的賦存條件及其是否能被人們所發(fā)現(xiàn)，在很大程度上取決于觀察者對其認識的水平。齊波夫秩在理論上可以取到無窮大，齊波夫秩越大，預測值越趨近于0。地質(zhì)相關行業(yè)標準決定齊波夫取值范圍，即預測值不應小于最小工業(yè)礦床儲量。

4.1 主礦體的鎳金屬量統(tǒng)計

用地勘時期5個主礦體的鎳金屬量為基數(shù)進行預測，各礦體按鎳金屬量從大到小排序并賦秩（見表 1）。

表1 主礦體鎳金屬量排序表

4.2 求已知礦體的秩和齊波夫預測秩

根據(jù)各已知礦體與已知最大礦體儲量比值（Fx／F1）乘以自然數(shù)的積，建立最接近自然數(shù)的數(shù)列（見表2）。

由公式（1）和（2）計算：nF2／F1，nF3／F1，nF4／F1，nF5／F1。其中n為自然數(shù)（齊波夫秩）。把計算結(jié)果接近自然數(shù)1，2，… ，6的值分別記為序列I，II，III，IV，V列于表2中，并用小括號注明其齊波夫秩。然后計算表2中每一行即齊波夫序列的均值和標準離差，也一并列于表2中。選取最優(yōu)序列的原則就是序列的標準離差最小。比較表2中的各序列的均值和標準離差，明顯的看出第I列的標準離差最?。?.014）。因此可以確定白家嘴子礦區(qū)鎳金屬量最大的F1單元的齊波夫秩為1，其余各已知單元的齊波夫秩（R秩）分別為4、5、13、19（見表2第一列括號內(nèi)數(shù)字）。

表2 各系列的均值和標準離差

4.3 求最大礦體的金屬量

根據(jù)公式（1）計算常數(shù)K，也就是預測的最大礦儲量：

4.4 預測結(jié)果及可靠度分析

由公式（2）計算各級理論儲量Yi= K/Ri（Ri為相應的等級值）等級值取到116（等級值取116時，鎳金屬量為3萬t，達到鎳礦床最小可采規(guī)模），計算結(jié)果見表3。

表3 鎳金屬量預測表

齊波夫預測結(jié)果表明：白家嘴子礦區(qū)預測的總金屬量為1859萬t，隨著這些年對金川銅鎳礦床勘查工作的不斷提高，先后總計查明546萬t鎳金屬量，根據(jù)表3的預測結(jié)果，仍有111個礦體累計約1313萬t金屬量的找礦潛力。

表4 白家嘴子礦區(qū)鎳金屬量預測誤差表

圖2 白家嘴子礦區(qū)鎳金屬量齊波夫定律分布圖

由表4和圖2可知，已查明的5大礦體均已預測出來，所預測的最大礦體即Ⅱ-1礦體。5大礦體的預測值與實際值絕對誤差很小，擬合度很高，說明預測過程可靠，預測結(jié)果符合齊波夫定律。需要注意的是，應用齊波夫定律進行資源預測對區(qū)域性、戰(zhàn)略性找礦有參考價值，但它不能指出礦床產(chǎn)出的位置，所以該方法與礦體定位預測相結(jié)合是今后發(fā)展的方向［6］。

5 應用條件

（1）應用齊波夫定律進行礦產(chǎn)資源預測，要求礦區(qū)內(nèi)礦體是在相同地質(zhì)時期、相同地質(zhì)條件、相同地質(zhì)背景下形成的，且各礦體資源量滿足齊波夫分布定律。

（2）齊波夫定律雖然在使用上方便、簡單，但在數(shù)學推導上存在空白區(qū)，致使難以對誤差大小進行判斷［10］。

6 結(jié)論

（1）根據(jù)齊波夫定律預測得出白家嘴子礦區(qū)總資源量1859萬t，已控制資源量546萬t，尚有1313萬t潛在資源量，其資源前景可觀。

（2）白家嘴子礦區(qū)最大礦體的資源量為342萬t，與預測出的資源量348萬t十分接近，兩者可相互對照、參考，說明用齊波夫定律來預測資源量是可行的。

（3）應用齊波夫定律預測方便、簡單，所需已知數(shù)據(jù)的數(shù)量較少，尤其適合于一些地質(zhì)信息較少的地區(qū)和老礦區(qū)去發(fā)現(xiàn)新礦床。

本文試用齊律對甘肅金昌白家嘴子礦區(qū)進行了找礦前景預測，預測結(jié)果在理論上符合數(shù)學模型及推演法則。本文旨在起到拋磚引玉的作用，希望和更多的地質(zhì)同仁們一起對用數(shù)學模擬進行找礦預測的方法進行交流探討。