趙 博

（河北省地球物理勘查院，河北 廊坊 065000）

地球化學(xué)異常是對(duì)礦產(chǎn)進(jìn)行勘測(cè)的重要方法，同時(shí)也是對(duì)礦區(qū)進(jìn)行評(píng)價(jià)的重要憑證和證據(jù)，其中的異常信息提取更是勘測(cè)過程中的關(guān)鍵部分。隨著大計(jì)算和人工智能的迅速崛起，針對(duì)異常提取的研究也有了很大的突破，方法也更加多樣化，如利用滑動(dòng)平均法來幫助勘察低背景區(qū)的異常情況，利用含量面積分形的方法來幫助確定異常的下限等等。因此，要根據(jù)不同的礦區(qū)實(shí)際地質(zhì)情況可以選擇不同的方法進(jìn)行有效的預(yù)測(cè)和勘察，同時(shí)還要進(jìn)行不同方法的比對(duì)，幫助提取更多的異常信息。

1 礦區(qū)地質(zhì)簡(jiǎn)析

圖1 地區(qū)主要結(jié)構(gòu)示意圖

金礦的分布位于準(zhǔn)格爾-哈薩克斯坦東北邊緣活動(dòng)帶附近，礦區(qū)內(nèi)斷裂分布走向差異較大，有NW 向斷裂、NNE 向斷裂以及近EW向斷裂，其中NW向斷裂占主要部分（圖1），而且礦區(qū)內(nèi)的礦產(chǎn)分布和規(guī)模大小受其影響[1]。礦區(qū)內(nèi)的礦石結(jié)構(gòu)主要有氧化礦和原生礦。其中氧化礦石金屬礦中褐鐵礦占多數(shù)，其次為黃鐵礦和毒砂。原生礦的主要成分為毒砂和黃鐵礦。而脈石礦物則有石英、方解石、正長石等。氧化礦的分布深度在地表下三十米，再往下大多是原生礦區(qū)域。

2 礦區(qū)化學(xué)元素的統(tǒng)計(jì)特征分析

2.1 元素組合特征

經(jīng)過對(duì)該金礦區(qū)的相關(guān)資料進(jìn)行研究之后，發(fā)現(xiàn)在其礦床元素中，Au 和Ag 與Pb 和Sn 與Mo 和As 等元素的襯度值大于一，同時(shí)它們的濃度克拉克值也明顯大于一，與其他元素對(duì)比來看相對(duì)較高，As 會(huì)低一點(diǎn)，而Cu，Pb，Zn 等元素則沒有表現(xiàn)出類似的異常，這就表明Au 是該礦區(qū)的主礦元素，經(jīng)過復(fù)雜的活化、變遷和富集作用之后，讓毒砂成為該礦區(qū)的主要硫化物[2]。此外，大量的As 元素，對(duì)于毒砂的富集也有著很好的促進(jìn)作用，使得該區(qū)域內(nèi)的毒砂含量豐富。

2.2 元素類聚

而經(jīng)過分析比對(duì)，不同于其他的元素與Au 的相似性差，W 和Mn 與Mo 和Hg 則顯示出一定的相似性，而與Pb、Zn、Ag、As 等元素則表現(xiàn)出相似度不高。W、Mn、Mo 和Au 則具有一些相同的特征，總結(jié)得出在該金礦區(qū)明顯存在著剝蝕的現(xiàn)象，就目前的開采情況來看，主要集中在原生礦區(qū)域的中下部分。Hg 和Au 存在相似特征，則表明在礦區(qū)的深部區(qū)域有隱伏礦體的存在。

2.3 因子分析

經(jīng)過對(duì)該礦區(qū)的地球化學(xué)元素分析比對(duì)發(fā)現(xiàn)，該區(qū)域內(nèi)的元素演變十分地復(fù)雜。因?yàn)榛鹕降囊蛩兀搮^(qū)域的巖體物質(zhì)形成受到火山物質(zhì)和周圍沉淀物的相互影響，導(dǎo)致物質(zhì)出現(xiàn)融合，而南明水組地層的水分則和大氣降水形成對(duì)流作用，演變成高溫?zé)嵋?，在溫度的作用下，不斷冷卻形成沉淀物，最后累積成礦床。首先，Au 元素的獨(dú)立性非常強(qiáng)，因?yàn)樵搮^(qū)域在經(jīng)過不斷的進(jìn)行活化、遷移等地層活動(dòng)，受到這些活動(dòng)的影響，從而出現(xiàn)大量的富集現(xiàn)象。這同樣給金元素積累成礦奠定了很好的基礎(chǔ)，同時(shí)因?yàn)檫@些地層活動(dòng)一直在進(jìn)行，且不斷地累積，導(dǎo)致了該區(qū)氧化礦的含量比原生礦的含量要高出很多。其次，因?yàn)榻鹪睾蚖、Cr、Bi、Mo 等元素存在一定的相似性的特征，表明了金元素的富集與W、Mo、Bi 元素有很大的關(guān)聯(lián)，為金元素富集成礦提供了非常好的條件。這些特征與礦區(qū)內(nèi)的硅化反應(yīng)，以及和鈉長石的蝕變，還有電氣石化的蝕變呈現(xiàn)出一致性。再者，砷元素的富集會(huì)經(jīng)歷三個(gè)以上的階段，金元素礦化主要出現(xiàn)在第二個(gè)階段，在金元素成礦之后，砷元素的富集規(guī)模會(huì)變得很大，導(dǎo)致其含量很高。雖然這個(gè)特征顯示砷元素的分布會(huì)影響著金元素的分布，但是這不是必備條件，因此不能把砷元素的含量分布當(dāng)作是尋找金礦的重要條件，根據(jù)其元素類聚分析來看，應(yīng)該著重于As，Sb，Bi，W，Mo，Hg等元素的集合，尋找它們的共同分布區(qū)域。在金元素的初步富集階段，呈現(xiàn)出與Mn，Ni，Co，As，Pb，Sn 元素有很強(qiáng)的相似性，根據(jù)這一特征可以得出，硫化物與金元素有很強(qiáng)的關(guān)聯(lián)性，有助于金元素的富集礦化。

3 礦區(qū)化學(xué)元素的空間分布

3.1 金元素分布在走向上的特征

在薩爾布拉克金礦區(qū)內(nèi)，呈現(xiàn)出NW-SE 向的分布狀態(tài)，呈現(xiàn)出整體分段礦化的現(xiàn)象。礦體展布為波狀，整體的礦體分布有抖傾狀的特點(diǎn)，在礦帶中的分散分布，其大小和品級(jí)都不盡相同。礦化較差的是東西側(cè)，存在明顯的蝕變礦化特征。同時(shí)，從所有的礦體元素分布來看，金元素的富集現(xiàn)象相對(duì)突出。

3.2 金元素分布在傾向上的特征

Au 元素和Hg，Sb，W，Bi，Pb，Ag，Mo，Sn 等元素整體表現(xiàn)出高氧化的現(xiàn)象，而在原生礦內(nèi)這種現(xiàn)象則比較少。在主斷裂的區(qū)域內(nèi)，Co，Mn，B，Cu，Ni，Zn 等元素的分布含量較少見，主要分布在破碎化程度高的氧化礦中，而在兩側(cè)區(qū)域的圍巖中，這些元素的含量則比較高。

3.3 金元素分布在垂向上的特征

金元素在地表垂直向下六十米富集現(xiàn)象明顯，在十米到四十米之間的位置，品級(jí)較高，且比較穩(wěn)定。同時(shí)，金元素的品級(jí)差異較大的現(xiàn)象，則表明其元素在礦區(qū)內(nèi)分布分散，這不利于對(duì)礦體的數(shù)量預(yù)測(cè)和開采。在垂直方向上，金元素有連續(xù)聚集的特征，從而顯現(xiàn)出多層礦化的特征。

4 地化異常的信息提取分析

近年來隨著礦質(zhì)勘測(cè)手段和技術(shù)的進(jìn)步，地球化學(xué)異常信息的提取方法也有了很大的進(jìn)步，針對(duì)不同的勘測(cè)要求，則有不同的信息提取方法進(jìn)行有針對(duì)性的應(yīng)用。本文主要針對(duì)傳統(tǒng)的統(tǒng)計(jì)方法和含量面積分形方法以及滑動(dòng)平均法這三種方法進(jìn)行闡述。

4.1 傳統(tǒng)的統(tǒng)計(jì)學(xué)方法

這種方式的關(guān)鍵在于數(shù)據(jù)要符合正態(tài)分布，亦或者要符合對(duì)數(shù)正態(tài)分布。主要利用Surfer 軟件對(duì)礦區(qū)內(nèi)的元素含量展開統(tǒng)計(jì)運(yùn)算。其對(duì)應(yīng)的特征值有最小值、最大值、平均值、中位數(shù)、標(biāo)準(zhǔn)離差、變化系數(shù)、偏度、峰度和K-S 吻合度。其中，檢測(cè)偏度則是為了檢查其正態(tài)分布的偏離程度，數(shù)值小于20 則表明符合正態(tài)分布，反之則符合對(duì)數(shù)正態(tài)分布。對(duì)于符合前者的數(shù)值，同時(shí)還要滿足大于平均值+2 倍的標(biāo)準(zhǔn)離差的概率要小于2.5%的條件。對(duì)于符合后者的數(shù)值，要滿足其值大于eμ+2σ的概率要小于2.5%的條件。然后進(jìn)行異常圈定，單一元素的異常下限檢測(cè)公式為T=C0+KS，測(cè)出T 的值后，再根據(jù)其值繪制異常圖。

這個(gè)方法是比較簡(jiǎn)單的計(jì)算異常下限的方法，優(yōu)點(diǎn)是容易操作，但是也有其弊端。因?yàn)樵谟?jì)算時(shí)要進(jìn)行有必要的踢值處理，這樣做就會(huì)破壞掉原有的數(shù)據(jù)，而且在K 值的選取上比較主觀，很難比較精確地反映出詳細(xì)的地質(zhì)情況。

4.2 含量面積分形方法

這個(gè)方法的關(guān)鍵是要大于或者等于某一條件的數(shù)，其關(guān)系模型為N(r)=Cr-D（r＞0）。其中，N 代表元素的含量，N(r)代表N 大于等于r 時(shí)圍成的面積，C 代表常數(shù)，D 代表分維數(shù)。為了得到公式中的D 值，要將(Nr1，Nr2，Nr3，...，Nrn)和(r1，r2，r3，...，rn)帶入到公式N(r)=Cr-D（r＞0）中，兩邊取對(duì)數(shù)，得到：logN(r)=Dlgr+lgC，然后得出D 的值。

這個(gè)方法相對(duì)于上面的統(tǒng)計(jì)方法，能夠客觀有效地反映出異常下限，而且不受到正態(tài)分布的限制和影響，是一種很好的異常分析方法。

4.3 滑動(dòng)平均法

該方法是根據(jù)具體的地質(zhì)圖，從而對(duì)區(qū)域進(jìn)行探測(cè)。首先要對(duì)采集的樣品進(jìn)行化學(xué)處理和分析，并對(duì)這些數(shù)據(jù)進(jìn)行網(wǎng)格化處理，每一個(gè)格子里面只能放一個(gè)值，如果在格子內(nèi)的是數(shù)個(gè)點(diǎn)，那么則要求取平均值，將其放置在格子內(nèi)。然后根據(jù)實(shí)際的地質(zhì)情況開窗口，窗口的選擇對(duì)于信息收集是非常重要的一步，過大則會(huì)漏掉一些異常較小的數(shù)據(jù)，過小則會(huì)受到取樣、分析過程的干擾。通常來說，大窗口主要應(yīng)用于地化背景變化的模擬，小窗口則應(yīng)用于局部變化中。

該方法的應(yīng)用是因?yàn)橹鲝垟?shù)據(jù)分布在空間上具有連續(xù)性的特點(diǎn)。基于這個(gè)特點(diǎn)，模擬出具有起伏的曲面，可以幫助突出低背景區(qū)異常較弱的信息，從而有針對(duì)性地找出符合要求的礦產(chǎn)區(qū)域。相對(duì)于前兩者來說，它圈定的異常范圍比較可控，不像含量面積方法的范圍那樣大，也不像統(tǒng)計(jì)法的那樣小。同時(shí)不需要對(duì)數(shù)據(jù)進(jìn)行一系列的處理，排除了干擾因素，是一種可控有效的異常提取方法。

5 結(jié)語

隨著經(jīng)濟(jì)和科技的指數(shù)型增長，我國的礦產(chǎn)消耗也隨著加劇。地化異常不但是進(jìn)行礦產(chǎn)數(shù)量估計(jì)的有效手段，同時(shí)也是尋找礦脈的重要方法。而隨著礦產(chǎn)開采的進(jìn)程推進(jìn)，礦產(chǎn)的探測(cè)技術(shù)和異常提取方面的研究也有了一定的突破，針對(duì)不同的開采要求，應(yīng)用不同的方法進(jìn)行地化數(shù)據(jù)分析。在薩爾布拉克金礦區(qū)的地質(zhì)分析基礎(chǔ)上，對(duì)元素特征和空間分布進(jìn)行了分析研究。在地化信息的提取方法中，傳統(tǒng)的統(tǒng)計(jì)方法比較簡(jiǎn)單，但受到的干擾較大，不夠精確。含量面積方法是一種很有效的異常分析方法，能夠客觀地反映出異常下限?；瑒?dòng)平均法則避免了前兩者的弊端，能夠可控的反映出地化異常，范圍可大可小，比較精確有效。