葉慶森

(核工業(yè)二三○研究所，長沙 410007)

0 前言

礦床原生暈指示元素垂向分帶規(guī)律的研究有助于正確判斷礦床(體)剝蝕程度，尋找隱伏礦床(體)，是構造地球化學原生疊加暈找礦和地球化學模式找礦的基礎。元素地球化學活動性的差異使得不同指示元素的原生暈在空間上呈現(xiàn)序次分布(組分分帶)。在二十世紀七十年代前蘇聯(lián)學者格里戈良提出指示元素垂向分帶指數計算法[1]的基礎上，國內化探界學者先后提出了多種計算方法，主要有廣義襯值法、濃集重心法、濃集指數法、比重指數法、含量梯度法等[2-9]，為元素分帶序列研究方法的不斷發(fā)展和完善做出了貢獻。下面就上述幾種主要的計算方法予以簡單介紹，并指出各自的特點。

1 分帶指數法

分帶指數法[1]是一種傳統(tǒng)的原生暈指示元素分帶序列計算方法，由前蘇聯(lián)學者格里戈良提出。它是根據指示元素分帶指數的最大值所在截面(中段)位置由淺至深將指示元素順序排列。

1.1 計算方法

以A礦床為例，其計算步驟可概括如下：

(1)計算各截面有關指示元素的線金屬量(表1)，并按元素逐一標準化至同一數量級。

(2)計算不同截面各指示元素的分帶指數值D(表2)。某指示元素的分帶指數最大值(Dmax)所在的截面即為該指示元素在分帶序列中的位置，由此就可初步排出它們之間自上而下的分帶序列為：Au、Ag、Pb-W-Sn、Cu、Co-Bi。如果某截面同時存在兩個或多個元素Dmax，則需進行下一步計算。

表1 A礦床8個指示元素在不同截面的線金屬量(m·10-6)

表2 A礦床8個指示元素的分帶指數值

(3)計算指示元素的分帶指數梯度值ΔG，比較ΔG的大小即可確定各元素在垂向分帶序列中的確切位置。ΔG大者在上，小者在下，設共有p個截面參與計算。

(1)

式(1)中D1、Dm、Dp分別是當兩個或多個指示元素的Dmax同在最上截面(如地表)、同在中部截面、同在最下截面時相應元素的Dmax值；Dl為某指示元素在l截面的分帶指數值。

從公式(1)可以看出，當指示元素的異常不連續(xù)即某截面無異常時，就必須人為地賦給Dl一定的數值，否則線金屬量Dl為“0”，就無法計算出ΔG的值。

如A礦床中Au、Ag、Pb和Sn、Cu、Co的Dmax分別同在最上截面(地表)和中部中段2，同時Sn、Cu、Co在某些截面出現(xiàn)了異常不連續(xù)(線金屬量為“0”)的情況。為計算ΔG，將指示元素標準化線金屬量為“0”的相應位置賦值0.001。計算結果表明，Au、Ag、Pb和Sn、Cu、Co自上而下的確切垂向分帶序列分別為Ag-Au-Pb和Co-Cu-Sn。

因此本例確切的分帶序列為：Ag-Au-Pb-W-Co-Cu-Sn-Bi。

需要指出，如果在表1中分別刪去Bi和Sn，利用分帶指數法分別重新計算剩余7個元素的分帶指數(表3、表4)和必要的ΔG，可分別得出垂向分帶序列Ag-Pb-Co-Cu -W-Au-Sn和Ag-Au-Pb-W-Co-Cu-Bi。

比較表2(8個元素)和表3(據表1中去掉Bi所在列計算)可以看出，兩者的垂向分帶序列發(fā)生了明顯變化，即由于參與計算的指示元素的個數不同，分帶序列發(fā)生了變更，如根據表2(8個元素)，Au排在第2位，但根據表3，Au卻排到第6位。同樣比較表3和表4(分別據表1中去掉Bi和Sn所在列計算)可見，兩者的分帶序列也有明顯變化，如依照表3，Au排到第6位，但依照表4，Au卻排在第2位了。

1.2 特點

縱觀上述計算過程可知，格里戈良分帶指數計算法雖然面世最早，也曾一度視為“經典”，但其不足之處是顯而易見的，如：

(1)該方法的計算過程比較復雜。

(2)當某指示元素在某中段異常不連續(xù)，線金屬量為“0”時，必須人為地賦給一個很小的非“0”值，否則垂向分帶序列無法計算。然而這個很小的非“0”值可能對整個計算結果產生較大的影響，這顯然是不客觀的。

表3 A礦床Au、Ag、W、Sn、Pb、Cu、Co的分帶指數值

表4 A礦床Au、Ag、W、Pb、Cu、Co、Bi的分帶指數值

(3)由于某指示元素的分帶指數值可隨參與計算的元素多寡、種類及其他元素線金屬量而變，從而影響其在分帶序列中的位置。由于各種條件的限制和影響，在實際工作中參與計算的指示元素種類和數量不可能且沒有必要完全統(tǒng)一。由此可見，由于參與計算的指示元素種類和數量的人為抉擇會使客觀存在的分帶規(guī)律發(fā)生改變。

綜上所述，該計算方法確實問題不少。

2 廣義襯值法

顧名思義，廣義襯值法[2]是從襯值或襯度(異常含量Ca與背景值Cb之比)的概念擴展而來。

2.1 計算方法

廣義襯值法的計算步驟如下：

(1)求指示元素在各截面的廣義襯值Km：

(2)

式(2)中Nm為m截面的樣品總數；Ci為指示元素在m截面樣品i的含量；Li為m截面樣品i的有效控制長度；Cb為指示元素的地球化學背景值。廣義襯值Km的意義相當于格里戈良法中分帶指數D。

(2)找出各指示元素廣義襯值最大值(Kmax)所在的截面位置，據此初步排出它們的垂向分帶序列。

(3)當某截面有兩個或多個指示元素的廣義襯值K同時為最大值時，則需按式(3)求廣義襯值梯度ΔGk(設有p個截面)。根據ΔGk的大小，可確定指示元素在垂向上的分帶序列，大者為上，小者為下。某指示元素的ΔGK定義如下：

(3)

式(3)中K1、Km、Kp分別是當兩個或多個指示元素的Kmax同在最上截面(如地表)、同在中部截面、同在最下截面時相應元素的Kmax值；Kl為某指示元素在l截面的廣義襯值。

根據該計算方法，可得到上述A礦床不同截面8個指示元素的廣義襯值(表5)和自上而下的垂向分帶序列(Ag-Au-Pb-W-Co-Cu-Sn-Bi)；即使表5中分別去掉Bi和Sn所在列后分別計算，得到的分帶序列為Ag-Au-Pb-W-Co-Cu-Sn和Ag-Au-Pb-W-Co-Cu-Bi，仍沒發(fā)生序次變化。

2.2 特點

(1)克服了格里戈良分帶指數計算法中當指示元素的異常在某截面不連續(xù)時，分帶指數梯度值ΔG無法直接計算的弊端。

(2)在廣義襯值法中，指示元素廣義襯值的計算是相互獨立的，因此廣義襯值不受其它元素的影響，即不會因參與計算的指示元素不同或多寡而使其在分帶序列中的位置發(fā)生變化。

(3)廣義襯值具有傳統(tǒng)襯值(即襯度)的概念，因而可揭示各指示元素在各截面的富集程度，從而能更客觀地反映指示元素分帶現(xiàn)象，地球化學意義明確。

(4)該方法數學意義非常明確，計算過程比格里戈良分帶指數法簡單得多，比如在廣義襯值法中，數據不需要標準化。

3 濃集重心法

該方法基于物理學原理。用剛性輕桿連結的n個質點的質心公式為

(4)

式(4)中mi為第i質點的質量；xi為第i質點在x軸上的坐標。

一般情況下，剛體的質心和重心相重合。如果把i截面的異?？闯梢粋€用輕桿連結的質點i，把指示元素的線金屬量(Mi)看成i點的質量(mi)，把標高(hi)看成i點的坐標(xi)，那么就可以求出該元素的濃集重心的空間位置，從而能夠定量研究元素的濃集重心在空間上的有序排列。

3.1 計算方法(設共有p個截面)

(1)計算各截面各指示元素的線金屬量(M)(也可以用異常含量和來代替)。

(2)參照上述質心公式，求出各指示元素的濃集重心(H)，根據H的大小排出由上而下的垂向分帶序列，大者在上，小者在下：

(5)

式中Ml為指示元素在l截面的線金屬量；hl為l截面異常標高。

3.2 特點

(1)在部分截面無異常的情況下，仍可計算分帶序列。

(2)計算出來的結果代表元素的濃集重心在空間上的分布位置，其數學意義和地球化學意義較明確。

(3)計算過程較簡便，且其間各指示元素相互獨立計算，即指示元素的濃集重心值不受其它元素的影響，因此不會因參與計算的指示元素不同或多寡改變元素在垂向分帶序列中的位置。

4 濃集指數法

4.1 計算方法

(1)計算各指示元素在所有截面上含量之和的平均值CM：

表5 A礦床8個指示元素的廣義襯值

(6)

式中Cl是指示元素在l截面的含量，假設共有p個截面。

(2)計算各截面各指示元素含量與相應元素的CM之比值，并把該比值定義為該截面該指示元素含量的標準化值。

(3)計算同一截面所有指示元素含量標準化值之平均值CMB。

(4)計算某截面不同指示元素的濃集指數CI。某截面某指示元素的濃集指數CI等于該截面該元素含量標準化值與該截面CMB之比，并按該指示元素的濃集指數最大值CImax所在截面位置初步確定分帶序列。因為該比值越大，表明相應截面相應指示元素的濃集程度越高，故將該比值稱為濃集指數。

(5)當某截面有兩個或多個指示元素的濃集指數CI同時為最大值時，則按式(7)求濃集指數梯度ΔGCI(設有p個截面)。根據ΔGCI的大小，可確定出元素在垂向上的分帶序列，大者為上，小者為下。某指示元素的ΔGCI定義為式(7)：

(7)

式(7)中CI1、CIm、CIp分別是當兩個或多個指示元素的CImax同在最上截面(如地表)、同在中部截面、同在最下截面時相應元素的CImax值；CIl為某指示元素在l截面的濃集指數。

4.2 特點

(1)濃集指數方法與廣義襯值法有類似之處，克服了格里戈良分帶指數計算法中當指示元素的異常在某截面不連續(xù)時，分帶指數梯度值ΔG無法直接計算的弊端。

(2)原始數據標準化的結果是唯一的，為無量綱化處理，能消除地殼中各元素豐度差可能給計算帶來的影響。

(3)盡管各指示元素濃集指數的計算不相互獨立，濃集指數值可因參與計算的元素多寡或種類而變，但初步驗算表明，參與計算的指示元素多寡或種類不同不會影響指示元素在分帶序列中的位置。

(4)該方法的計算過程較復雜。

5 比重指數法

5.1 計算方法

(2)求指示元素i在l截面的比重指數Hli(即某指示元素在某截面的線金屬量占該指示元素在所有截面線金屬量總和之比例)：

(8)

式(8)中Mli為指示元素i在l截面的線金屬量，設共有m個截面(l=1，2，…，m)。

某指示元素i在某截面的比重指數Hi值大，則該元素在該截面所占的比重大。某指示元素i的比重指數最大值Himax所在的截面位置(自上而下)就是該元素在垂向分帶中的位置。

(3)當兩個或多個指示元素的Himax同在某截面時，可通過直接比較Himax的大小，來確定它們的具體位置：

當Himax同在最下截面時，大者為下。

當Himax同在某中部截面時，則相應指示元素i在該截面及其以上截面的Hi值之和∑M上大者或以下截面的Hi值之和∑M下小者為上。

5.2 特點

比重指數方法與廣義襯值法有類似之處：

(1)無需標準化計算，計算過程比較簡單。

(2)比重指數Hi值是通過同元素同量綱的關系來計算所占比重的，即指示元素Hi值的計算是相互獨立的，不會因參與計算的指示元素不同或多寡而使元素在分帶序列中的位置發(fā)生變化。

(3)該方法只需計算元素比重指數Hi，不必計算分帶指數梯度值ΔG，因而克服了當指示元素異常在某截面不連續(xù)時，ΔG無法直接計算的弊端。

6 含量梯度法

含量梯度法通過直接比較同一指示元素在各截面的含量大小來確定指示元素的垂向分帶序列。它與概率梯度法[10]基本相同，但比后者更簡便。

6.1 計算方法

(1)列出參與計算的所有指示元素在各截面中的含量。

挑取1～2環(huán)接入細菌斜面培養(yǎng)基上，37℃培養(yǎng)24 h后，用無菌生理鹽水洗下菌苔，制成菌懸液。采用麥氏比濁法，用無菌生理鹽水調節(jié)細菌含量為1×108cfu/ml。

(2)逐一比較同一指示元素在各截面中的含量，含量最大值Cmax所在的截面位置即為該元素在垂向分帶中的位置，據此得出基本的分帶序列。

(3)當兩個或多個指示元素的Cmax同在某截面時，則需比較含量梯度ΔGC的大小來確定指示元素在分帶序列中的確切位置。ΔGC大者為上，小者為下。

某指示元素的ΔGC定義為式(9)：

(9)

式中C1、Cm、Cp分別是當兩個或多個指示元素的Cmax同在最上截面(如地表)、同在中部截面、同在最下截面時相應元素的Cmax值；Cl為某指示元素在l截面中的含量。

6.2 特點

(1)含量梯度法與廣義襯值法、比重指數法等有類似之處，均通過同元素同量綱的關系來計算所占比重，即各指示元素含量梯度值的計算相互獨立，因而不會因參與計算的指示元素不同或多寡而使指示元素在分帶序列中的位置發(fā)生變化。

(2)該方法數據不需要標準化，計算過程更簡便。

(3)當指示元素的異常在某截面不連續(xù)時，不會造成分帶序列無法直接計算的現(xiàn)象。

(4)不足之處是該方法只考慮指示元素在各截面、各采樣點的含量，忽略了各采樣點指示元素含量的有效控制長度，即忽略了其含量的權重。因此為客觀全面起見，在確定各截面指示元素含量時，必須考慮各采樣點的有效控制長度。

7 結語

上述六種確定礦床原生暈指示元素垂向分帶序列的計算方法是近年來地質文獻中應用頻率較高的方法。作者認為：

(1)雖然廣義襯值法、濃集重心法、濃集指數法、比重指數法、含量梯度法一般都能克服格里戈良分帶指數計算法中當指示元素的異常在某截面不連續(xù)時，分帶指數梯度值ΔG無法直接計算的弊端，但由于濃集指數法中各指示元素濃集指數值的計算并非相互獨立，即同一礦床原生暈指示元素的濃集指數值可因參與計算的元素多寡或種類而變，而且其計算過程較復雜；含量梯度法僅考慮指示元素在各截面的含量，忽略了各采樣點指示元素含量的有效控制長度，因而該方法尚需完善和改進。由此可見，廣義襯值法、濃集重心法和比重指數法相對合理和客觀，其計算方法比較簡單，值得推廣。

(2)格里戈良分帶指數法雖然比較經典，但由于異常不連續(xù)，即某截面無異常時，其分帶指數梯度值ΔG無法直接計算而影響排序；而且某指示元素在分帶序列中的位置會因參與計算的指示元素不同或多寡而改變，因此不值得推廣。

(3)值得指出的是文獻[7]未正確介紹格里戈良分帶指數法(該方法并非將“各元素的含量換算成同一數量級次”)；文獻[8]介紹的數據標準化，未考慮到不同指示元素的含量級差大(如Au和Cu)時，無法實現(xiàn)標準化到同一數量級的情況。

(4)在計算礦床原生暈指示元素分帶序列前，必須查明控礦因素、綜合分析成礦地質條件、分解成礦作用期次及各期次指示元素組合，以便建立單一成礦期次的指示元素垂向分帶理想模型，合理解釋由于構造疊加暈引起的元素分帶序列“反常及反分帶現(xiàn)象”[11]。這對指導化探異常評價、提升礦產預測水平、提高模型找礦的效果意義重大。

參考文獻：

[1] 地質礦產部情報研究所．勘查地球化學手冊：第二卷[M]．北京：冶金工業(yè)出版杜，1990．

[2] 葉慶森，馮倡，滕潤球. MS巖體中段鈾礦床地球化學異常模式及其找礦意義[J]. 鈾礦地質，1988,4(6)：340-347.

[3] 樸壽成，連長云. 一種確定原生暈分帶序列的新方法——重心法[J]. 地質與勘探，1994,30(1)：63-65.

[4] 解慶林. 濃集指數法確定礦床原生暈元素軸向分帶序列[J]. 地質與勘探，1992,33(6)：55.

[5] 趙琦. 原生暈垂直分帶的元素比重指數計算法[J]. 物探與化探. 1989,13(2)：157-159.

[6] 張定源. 銀巖錫礦原生暈元素分帶序列計算方法研究[J]. 地質與勘探，1989,25(6)：45-49.

[7] 王建新，臧興運，郭秀峰，等. 格里戈良分帶指數法的改良[J]. 吉林大學學報：地球科學版， 2007,37(5)：884-888.

[8] 邢利琦，劉炳璋. 礦床原生地球化學暈分帶性研究[J]. 四川地質學報，2011,31(4)：489-495.

[9] 金獻革. 礦床原生暈元素空問分帶序列計算方法研究[J]. 科技與企業(yè)，2012(20): 274-275.

[10] 邱德同. 確定礦床原生暈指示元素分帶序列的新方法[J]. 地質與勘探，1989,30(8)：51-53.

[11] 李惠，禹斌，李德亮，等. 構造疊加暈找盲礦法及研究方法[J]. 地質與勘探，2013,49(1): 154-161.