黃支剛 ,陳杜軍,宋 昊

(1.成都理工大學(xué) 地球化學(xué)系，成都 610059；2.西藏思明礦業(yè)股份有限公司， 拉薩 850000；3.內(nèi)蒙古地質(zhì)礦產(chǎn)勘查院，呼和浩特 010011)

0 引言

通過一次面積性的水系沉積物地球化學(xué)測量工作，可能發(fā)現(xiàn)一批異常，但并不是所有的異常都與礦化有關(guān)，由工業(yè)礦產(chǎn)引起的異常僅占少數(shù)，大多數(shù)異常由非礦化因素或分散礦化引起。在實(shí)際工作中，往往在工區(qū)內(nèi)圈出多個化探異常靶區(qū)，由于資金和人力資源等方面的限制，不可能對所有的異常進(jìn)行充分的研究。通過對異常進(jìn)行解釋推斷，初步篩選出與礦有關(guān)的異常，為了更好地安排異常檢查和驗(yàn)證工作，需要對異常的找礦遠(yuǎn)景進(jìn)行排序，達(dá)到快速縮小靶區(qū)，高效找到有經(jīng)濟(jì)意義礦床的目的[1]。

針對異常排序的方法有很多，各有側(cè)重點(diǎn)：有側(cè)重于地球化學(xué)指標(biāo)(如異常面積、異常強(qiáng)度、異常規(guī)模、元素組合特征等)；也有側(cè)重于成礦地質(zhì)條件(如巖漿巖、地層巖性、地質(zhì)構(gòu)造等[2])。排序時(shí)可以使用單一參數(shù)，也可以使用多個參數(shù)，但是往往一個工區(qū)內(nèi)，很多異常差異不大，使用單一的參數(shù)進(jìn)行排序往往效果不好，需要使用多個參數(shù)結(jié)合起來進(jìn)行排序[1]。使用不同參數(shù)排序所得的異常序列不同，顯然任意選擇其中一種序列均有不合理之處。因此綜合各種序列，給出一個定量決策模型，排出合理的序列就顯得十分必要。

為了對供選擇的方案集合U=(u1,u2,…,un)中的元素進(jìn)行排序，由m個專家組成的專家組M分別對U中的元素進(jìn)行排序，得到了m種意見，即V=(v1,v2,…,vm)，而這些意見都是模糊的，將這m種意見集中為一個較為合理的意見，稱為“模糊意見集中決策”[3]。

本次研究使用模糊意見集中決策方法，嘗試對化探異常的找礦前景給出合理的序列。

1 研究區(qū)地質(zhì)概況

研究區(qū)面積1 500 km2，地處班公湖-怒江成礦帶，區(qū)內(nèi)晚古生代以來地層發(fā)育較全，地質(zhì)構(gòu)造十分復(fù)雜，燕山期巖漿活動強(qiáng)烈，成礦條件非常有利(圖1)。研究區(qū)受到強(qiáng)烈的燕山期逆沖推覆作用及喜馬拉雅期右行剪切走滑作用，產(chǎn)生了大量近東西向的斷裂，控制了中酸性巖體的空間展布，其中區(qū)內(nèi)的許多基性巖脈及少量玄武巖體，可能與深斷裂相關(guān)。出露巖漿巖有侏羅紀(jì)的中細(xì)粒二長花崗巖、二長閃長巖、石英閃長巖，為淺成相和中淺成相。出露地層有二疊系的龍格組灰?guī)r、吞龍共巴組砂巖與板巖、曲地組砂巖等[4]。

圖1 研究區(qū)地質(zhì)略圖Fig．1 Geological sketch map of the study area

綜上所述，研究區(qū)成礦條件優(yōu)越，通過野外調(diào)查發(fā)現(xiàn)，銅礦化現(xiàn)象明顯，找礦前景較好，容易形成斑巖型和矽卡巖型兩種類型的銅礦床。因此，本次擬選用銅礦作為優(yōu)選礦種進(jìn)行異常排序。

2 1∶50 000水系沉積物測量

2.1 背景值和異常下限的確定

在研究區(qū)及其外圍進(jìn)行1∶50 000水系沉積物測量，采集水系沉積物的粒度為+20目～+80目，銅元素異常下限的確定采用計(jì)算法。

研究區(qū)的范圍依據(jù)水系沉積物測量所發(fā)現(xiàn)的異常圈定。首先計(jì)算在研究區(qū)內(nèi)銅元素的算術(shù)平均值和標(biāo)準(zhǔn)差；其次逐步剔除銅含量在算術(shù)平均值±3倍標(biāo)準(zhǔn)差之外的樣品，將未剔除的樣品作銅含量頻數(shù)直方圖，銅含量服從正態(tài)分布；然后計(jì)算研究區(qū)內(nèi)銅元素的算術(shù)平均值(Co)和標(biāo)準(zhǔn)差(σ)；最后利用公式(1)求出銅元素的異常下限(Ca)[5]。

Ca=Co+Kσ

(1)

其中Ca為異常下限值；Co為背景值；σ為標(biāo)準(zhǔn)差；K為試驗(yàn)常數(shù)，本研究區(qū)取2。由此確定了本區(qū)銅元素的異常下限為23×10-6。

2.2 銅異常的分布特征及初步篩選

研究區(qū)1∶50 000水系沉積物測量成果顯示，區(qū)內(nèi)的元素組合以銅鉛鋅為主。從銅元素來看，在研究區(qū)共圈選出13個異常區(qū)。本區(qū)圈選出的異常區(qū)，大小不一，形態(tài)不同；異常區(qū)的面積從數(shù)平方千米到數(shù)十平方千米不等，最大的可達(dá)24.3 km2；較大的異常區(qū)濃集中心相對較明顯，而小的異常區(qū)濃集中心強(qiáng)度相對較弱；元素組合相對簡單，銅鉛鋅異常套合較好。根據(jù)野外調(diào)查和室內(nèi)研究分析，剔除4個異常區(qū)，這4個異常區(qū)的異常面積較小，為數(shù)平方千米，異常強(qiáng)度低，元素組合較簡單，野外地質(zhì)調(diào)查也未發(fā)現(xiàn)有很好的礦化現(xiàn)象。剩下的9個異常區(qū)，分別對其進(jìn)行編號，為u1～u9。

3 異常排序的參數(shù)選擇

盡管化探是一種直接快速高效的找礦方法，實(shí)踐證明僅依靠化探來找礦效果并不理想，還需要與地質(zhì)條件進(jìn)行配合，才能對異常做出更符合實(shí)際的評價(jià)。常用于異常排序的參數(shù)有元素的異常強(qiáng)度、異常面積、異常規(guī)模、元素組合特征等，這些評價(jià)參數(shù)屬地球化學(xué)指標(biāo)。成礦遠(yuǎn)景不僅取決于地球化學(xué)指標(biāo)，不能僅憑異常的“高、大、全”來評價(jià)，因?yàn)楫惓５摹坝袩o、高低”受所選取的異常篩選準(zhǔn)則控制?！案?、大、全”異常不一定能成大礦，“弱小”異常不一定就成小礦，“無”異常也不一定就沒有礦，受礦化的剝蝕程度影響很大。因此，這就要求我們評價(jià)異常時(shí)還應(yīng)根據(jù)成礦地質(zhì)條件進(jìn)行綜合評價(jià)，為此，作者在進(jìn)行異常排序中引入“成礦地質(zhì)條件”參數(shù)。本研究確定選擇異常規(guī)模、元素組合特征、成礦地質(zhì)條件，作為異常排序的三個參數(shù)。

4 元素組合特征、成礦地質(zhì)條件的近似定量描述

由于元素組合特征、成礦地質(zhì)條件是描述性的，難以直接利用這兩個參數(shù)分別對異常進(jìn)行排序，需要把這兩個參數(shù)量對成礦遠(yuǎn)景反映的好壞程度轉(zhuǎn)化為數(shù)值來表示。這兩個參數(shù)的量化只能是近似的、模糊的，通過選擇使用“專家調(diào)查打分法”實(shí)現(xiàn)，即先由8位在該地區(qū)有多年找礦工作經(jīng)驗(yàn)的地質(zhì)工程師，對該地區(qū)9個異常的按元素組合異常及成礦地質(zhì)條件分別進(jìn)行打分，每個異常的每個參數(shù)的滿分為100分。將這些分值匯總并將統(tǒng)計(jì)結(jié)果交給各位專家研究，然后讓專家重新打分，將這次針對每個異常所打分值的平均值分別作為該異常的元素組合特征和“成礦地質(zhì)條件”的得分，分值越高表明該異常區(qū)的成礦潛力越大。

元素組合特征的量化，優(yōu)先選擇Cu、Pb、Zn作為量化的組合元素，如果異常區(qū)內(nèi)這三個元素均有異常、相關(guān)系數(shù)高、濃集中心明顯，異常位置套和好，則該異常的元素組合特征得分越高。

“成礦地質(zhì)條件”參數(shù)的量化，考慮的主要因素有區(qū)域地質(zhì)背景、成礦地質(zhì)環(huán)境(地質(zhì)構(gòu)造背景、地層巖性、控礦構(gòu)造、巖漿巖等)等，對成礦有利的因素，會提高“成礦地質(zhì)條件”參數(shù)的得分。

5 異常排序

根據(jù)異常規(guī)模(使用Cu面金屬量來度量)、元素組合特征、成礦地質(zhì)條件的量化數(shù)值，對研究區(qū)的9個異常分別進(jìn)行了排序，異常區(qū)的屬性特征及排序結(jié)果見表1。

表1 研究區(qū)銅異常的多參數(shù)排序表

Tab.1 Multi-parameters ranking list of Cu anomany in study area

異常編號異常規(guī)模(Cu面金屬量)序數(shù)元素組合特征序數(shù)成礦地質(zhì)特征(滿分100分)序數(shù)u1282.47626578u2255.28479607u3612.53881912u4232.49518549u5723.32823863u6435.25764834u7376.56587656u8761.61655745u9517.84862941

設(shè)這9個異常為U=﹛u1,u2,u3,u4,u5,u6,u7,u8,u9﹜，依據(jù)銅礦異常區(qū)的異常規(guī)模、元素組合特征和成礦地質(zhì)條件，分別對U中的9個異常區(qū)進(jìn)行排序，得到3種排法，即V=﹛v1,v2,v3﹜。如果第j個異常在第i個排法中位于第k位，令Bi(uj)=n-k，那么稱Bi(uj)為uj的Borda數(shù)[3]。此時(shí)U的所有元素可按Borda數(shù)的大小排序，此排序就是比較合理的。

評價(jià)異常的三個參數(shù)(異常規(guī)模、元素組合特征、成礦地質(zhì)條件)，在異常的整體評價(jià)中重要程度基本相同，因此給予相同的權(quán)重，而按異常的某特征或指標(biāo)排列出來的名次，對找礦而言其重要性是不同的，因此名次的權(quán)重也是不一樣的。名次降低，成礦前景降低，根據(jù)Borda數(shù)確定權(quán)重，權(quán)重值相應(yīng)下降，權(quán)重值見表2。使用模糊意見集中決策對研究區(qū)的異常作出排序，排序表見表2。

1)按表2所給的權(quán)重計(jì)算加權(quán)的Borda數(shù)。

B(u1)= 3×0.07+4×0.09+2×0.05=0.67;

B(u2)= 2×0.05+1×0.04+3×0.07=0.35;

B(u3)= 7×0.15+9×0.19+8×0.17=4.12;

B(u4)= 1×0.04+2×0.05+1×0.04=0.18;

B(u5)= 8×0.17+7×0.15+7×0.15=3.46;

B(u6)= 5×0.11+6×0.13+6×0.13=2.11;

B(u7)= 4×0.09+3×0.07+4×0.09=0.93;

B(u8)= 9×0.19+5×0.11+5×0.11=2.81;

B(u9)= 6×0.13+8×0.17+9×0.19=3.85;

2)按加權(quán)Borda數(shù)集中后的排序?yàn)椋?/p>

u3，u9，u5，u8，u6，u7，u1，u2，u4

6 不同異常排序方法的對比與討論

最先使用異常強(qiáng)度對異常排出的序列(u8,u5,u3…)進(jìn)行查證。結(jié)果表明8號異常，雖然異常面積大，但是野外調(diào)查發(fā)現(xiàn)該異常是由灰白色中細(xì)粒二長花崗巖體里的數(shù)個銅礦化破碎帶引起。礦化破碎帶規(guī)模較小，寬度0.3 m～1 m不等，延伸短(20 m～60 m)，由于這些銅礦化破碎帶距離比較近，計(jì)算機(jī)做異常等值線圖時(shí)，把這幾個銅異常圈選在一起構(gòu)成大異常，由于鉛鋅礦化在破碎帶中呈星散狀，在該異常區(qū)內(nèi)出現(xiàn)數(shù)個Pb、Zn小異常，由此推斷8號異常找礦遠(yuǎn)景小。

表2 銅異常區(qū)的模糊意見集中決策排序表

9號、3號異常雖然異常強(qiáng)度不如8號異常大，但是野外調(diào)查發(fā)現(xiàn)，這兩個異常均位于一級水系口附近，由此導(dǎo)致的異常面積較小，相比8號異常，其異常規(guī)模也較小。9號異常產(chǎn)于灰白色中細(xì)?；◢忛W長巖巖榴與灰?guī)r的接觸帶上；而3號異常產(chǎn)于灰?guī)r的矽卡巖化破碎帶中，初步認(rèn)為受層間破碎帶控制。以上兩個異常均呈帶狀，9號異常礦化體寬約1 m～2 m，可見延伸100 m；3號異常礦化體受層間破碎帶控制，礦化體寬約0.5 m～1.2 m，可見延伸180 m。揀塊法單點(diǎn)樣品分析結(jié)果表明，9號樣品品位：Cu=23.4%、Pb=5.3%、Zn=2.9%;3號樣品品位Cu=17.61%、Pb=3.4%、Zn=2.2%。這兩個礦化體規(guī)模比8號大，異常元素組合也齊全，成礦遠(yuǎn)景比8號大。

由此看來，僅依據(jù)異常規(guī)模(異常面積、異常強(qiáng)度)，有時(shí)對異常排出的序列，不能真實(shí)客觀地反映出異常的找礦遠(yuǎn)景；而使用模糊意見集中決策方法排出的序列，因?yàn)榭紤]了影響成礦遠(yuǎn)景的多種因素，從多個角度來研究異常，綜合了得到的多種序列，給出了一個較合理的序列。即8號異常雖然強(qiáng)度大，但是元素組合及成礦地質(zhì)條件均較差，而把成礦遠(yuǎn)景較大的9號、3號異常排在8號異常的前面。

水系沉積物測量的結(jié)果，受礦化體剝蝕程度的影響很大，但是在異常查證階段要查明礦化體的剝蝕程度，還需大量人力、物力及時(shí)間，這是不合理的，因此本研究暫時(shí)未考慮剝蝕要素對化探異常的影響。通過所獲盡量多的野外及室內(nèi)資料，來推測異常的找礦遠(yuǎn)景，在異常查證階段是有實(shí)際意義的。

7 結(jié)論

經(jīng)野外查證證實(shí)，將元素異常組合、成礦地質(zhì)條件引入到評價(jià)參數(shù)中，使用模糊意見集中決策給出的異常序列，能克服只從單一角度評價(jià)異常的缺陷，可以從多角度研究異常，對異常做出更合理的評價(jià)。該方法簡單可行，是一種有效的異常排序方法，能更好地指導(dǎo)異常查證工作。

參考文獻(xiàn)：

[1] 楊小峰，劉長垠，張?zhí)┤?，? 地球化學(xué)找礦方法[M]. 北京：地質(zhì)出版社，2007.

[2] 彭梧山. 地球化學(xué)探礦[M]. 北京：地質(zhì)出版社，1992.

[3] 呂國安. 成礦區(qū)帶地球化學(xué)異常評價(jià)方法[M]. 北京：冶金工業(yè)出版社，2002.

[4] 謝季堅(jiān)、劉承平. 模糊數(shù)學(xué)及其應(yīng)用(第三版)[M]. 武漢：華中科技大學(xué)出版社，2006.

[5] 江西省地質(zhì)調(diào)查研究院. 中華人民共和國區(qū)域地質(zhì)調(diào)查報(bào)告1∶25萬羌多幅[R]. 2005.

[6] 羅先熔，文美蘭，歐陽菲，等. 勘查地球化學(xué)[M]. 北京：冶金工業(yè)出版社，2007.

[7] 朱華平，張德全. 區(qū)域化探異常的地球化學(xué)勘查評價(jià)方法技術(shù)進(jìn)展綜述[J]. 地質(zhì)與勘探，2003，39(3)：35-38.

[8] 袁見齊，朱上慶，翟裕生．礦床學(xué)[M]．北京：地質(zhì)出版社，1985