王 樂， 辛存林， 趙 振， 廖建棟

(西北師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院，蘭州 730070)

新疆惠井子南山地區(qū)地球化學(xué)異常信息的識(shí)別與綜合評價(jià)

王 樂， 辛存林， 趙 振， 廖建棟

(西北師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院，蘭州 730070)

通過泛克里金法和多重分形濾波(S-A)法，對新疆惠井子南山地區(qū)成礦能量的綜合異常信息進(jìn)行識(shí)別、分析和評價(jià)，從而獲取不同層次的綜合異常信息及控礦因素。研究表明：①新疆惠井子南山地區(qū)成礦能量的最大變程為2.4 km,成礦能量連續(xù)性最佳方向?yàn)?20°，與研究區(qū)內(nèi)的主斷裂帶F2斷裂的走向基本一致，其成礦能量在變程范圍內(nèi)表現(xiàn)出較強(qiáng)的連續(xù)性和結(jié)構(gòu)性；②低通濾波獲取的區(qū)域成礦能量的異常在空間分布上與泛克里金法獲得的異常在空間上具有較強(qiáng)的相似性，其成礦能量巖控和裂控作用明顯；③多重分形濾波(S-A)法能夠提取多層次的致礦異常信息，其中包括了成礦能量的區(qū)域異常和隱伏異?；蛘弑尘爸械娜醍惓：碗[蔽異常。通過多重濾波方法提取的不同層次的異常，為下一步工作提供了新的研究靶區(qū)。

新疆惠井子； 成礦能量； 泛克里金法; 多重分形濾波法

0 引言

關(guān)于分析地球化學(xué)元素的空間特征、元素異常信息的分布位置和分布趨勢的研究，成為當(dāng)下地球化學(xué)研究的主要課題[1-4]。為此，地質(zhì)學(xué)家和地球化學(xué)家一直在探索一種既能表征地球化學(xué)元素含量空間分布規(guī)律，又能簡單、快速、準(zhǔn)確地確定元素異常下限的方法。許多學(xué)者發(fā)現(xiàn)地球化學(xué)元素含量在空間上具有一定的隨機(jī)性，除此之外，地球化學(xué)元素含量在空間上也具有一定的結(jié)構(gòu)性[5]。這種隨機(jī)性特征主要是由某些異常致礦因素導(dǎo)致的，其在尺度上是不變的，這種不變性可以用分形和多重分形模型來確定。

地質(zhì)統(tǒng)計(jì)學(xué)認(rèn)為地球化學(xué)元素的結(jié)構(gòu)性能夠反映其在空間上的自相關(guān)性特征， 這種特征在一定程度上也反映了地球化學(xué)元素在某種特定形成過程中所表現(xiàn)出的連續(xù)性，在這種過程中還有一些不明顯的地質(zhì)特征往往呈現(xiàn)出較強(qiáng)的隨機(jī)性[6]。因此，需要采用一些特定的數(shù)學(xué)模型來更加準(zhǔn)確地反映地球化學(xué)元素在空間上的某些分布變化，黃靜寧等[7]采用變差函數(shù)模型，在此基礎(chǔ)上結(jié)合泛克里金法，獲得了較為準(zhǔn)確地找礦信息[7]；Cheng[8-10]探索將非線性理論和方法應(yīng)用到深層次的致礦異常上，并且取得了不錯(cuò)的效果。與此同時(shí)，Cheng[11]又將多重分形的模型推廣到傅里葉能譜空間，提出廣義多重分形濾波方法，通過加強(qiáng)區(qū)域內(nèi)數(shù)據(jù)的局部結(jié)構(gòu)信息，從而達(dá)到對異常信息地識(shí)別、提取和確定，還可以將異常區(qū)域的背景和異常進(jìn)行分離。

筆者以新疆惠井子南山地區(qū)為研究對象，首先采用成礦能量法對該地區(qū)的地球化學(xué)元素異常信息進(jìn)行綜合分析，從而減少由于采樣間距的不均勻?qū)е碌臄?shù)據(jù)信息的過度分散或者集中。然后采用泛克里金法和多重分形濾波(S-A)法對研究區(qū)內(nèi)化探數(shù)據(jù)的空間變化特征進(jìn)行分析，提取研究區(qū)內(nèi)的致礦異常信息，并且結(jié)合研究區(qū)的地質(zhì)背景，對研究區(qū)的地球化學(xué)異常信息進(jìn)行評價(jià)。

1 研究區(qū)地質(zhì)特征

研究區(qū)位于卡拉麥里莫?dú)J烏拉板塊縫合帶南側(cè)(圖1)，屬于塔里木板塊博格達(dá)-哈爾里克晚古生代島弧[12-13]，隸屬于博格達(dá)-哈爾里克銅多金屬成礦帶。該成礦帶內(nèi)成礦地質(zhì)條件良好，具有很好地找礦潛力。出露的地層從老到新依次為古生界中石炭統(tǒng)居里得能組和上石炭統(tǒng)沙雷塞爾克組。研究區(qū)內(nèi)構(gòu)造作用和侵入作用強(qiáng)烈，構(gòu)造呈NNW方向分布在研究區(qū)中部，侵入體主要為石炭紀(jì)的安山玢巖，分布于研究區(qū)的中北部,區(qū)內(nèi)礦化蝕變作用明顯，巖石有不同程度的硅化、角巖化，局部的巖石比較破碎，常見褐鐵化現(xiàn)象。

圖1 東天山北部大地構(gòu)造單位略圖[14]Fig.1 The sketch map of the geotectonic units in the northern part of eastern Tianshan

2 工作方法

2.1 采樣點(diǎn)的布置

研究區(qū)采樣點(diǎn)的布置嚴(yán)格按照《地質(zhì)礦產(chǎn)勘查測量規(guī)范》(GB/T18341-2001)的要求，以100 m×250 m的距離為采樣間隔，取得巖屑化探分析樣693件，具體采樣點(diǎn)見圖2，但是由于有些地段地形的影響造成了部分采樣點(diǎn)的分散抑或集中，從而導(dǎo)致采樣點(diǎn)的不均勻分布。

圖2 研究區(qū)采樣點(diǎn)位圖Fig.2 Sampling points in the study area

2.2 樣品分析

本次樣品測試由核工業(yè)新疆理化分析測試中心完成，共完成測試樣品693件，分析了Au、Ag、Bi、Cu、Mg、Mo、Pb、Sb、Zn共9種元素，樣品分析方法由實(shí)驗(yàn)室根據(jù)自身實(shí)驗(yàn)儀器設(shè)備條件選擇確定成熟的分析方法，具體的分析方法見表1。其中Au、Ag、Bi、Mo、Pb元素的合格率達(dá)到100%，Cu、Mg元素的合格率為99.86%，Sb、Zn元素的合格率為99.52%。樣品的分析質(zhì)量可靠，能夠滿足工作地需要。

表1 各元素分析方法

3 研究區(qū)地球化學(xué)元素含量的統(tǒng)計(jì)

對研究區(qū)內(nèi)所采集的693件樣品中Au、Ag、Bi、Cu、Mg、Mo、Pb、Sb、Zn 9種微量元素統(tǒng)計(jì)分析，得到9種微量元素的平均含量、標(biāo)準(zhǔn)離差、元素分異系數(shù)和濃度克拉克值等參數(shù)(表2)。由表2可知，元素濃度克拉克值大于“1”的元素為Bi、Pb、Sb，說明這3種元素分布相對富集。變異系數(shù)大于“1”的元素有Mo、Sb、Au，說明這三種元素分布不均勻，分散程度高，活化和遷移作用強(qiáng)烈。從以上的分析中可以得出，Sb元素活動(dòng)性強(qiáng)(變異系數(shù)大于“1”)，濃度克拉克值高(濃度克拉克值接近于“1”)，這說明Sb元素在該地區(qū)相對富集，其分布不均勻，分散程度較高。在各種地質(zhì)作用和有利的地層條件下，極易富集成礦。因此，Sb元素作為研究區(qū)的主攻元素。

表2 惠井子地區(qū)地球化學(xué)元素含量統(tǒng)計(jì)表

注： Au元素的含量單位為10-9，其余元素的含量單位為10-6

4 泛克里金法對成礦能量綜合異常的確定

4.1 成礦能量的理論方法

成礦能量指的是成礦元素由原始分散狀態(tài)組合成富集狀態(tài)時(shí)所需的地球自然能。也就是說地球化學(xué)元素在成礦時(shí)元素富集或分散過程需要消耗的能量。具體到某個(gè)礦床就是利用礦區(qū)內(nèi)元素的含量值與元素的背景值的比值(克拉克濃度值)來反映能量地變化[14]。單位體積礦石(暈)上成礦能量(En)的計(jì)算公式可以用式(1)來表示[15]。

(1)

式中：En表示由n個(gè)元素在成礦過程中所消耗的總能量,主要是由礦石中濃度克拉克值k較高元素決定；n為成礦、成暈時(shí)元素的個(gè)數(shù)；ki為成礦過程中第i個(gè)元素的克拉克濃度值。這種方法在描述綜合區(qū)域的化學(xué)特征、圈定遠(yuǎn)景區(qū)、確定含礦建造屬性以及估計(jì)礦化規(guī)模等方面有很好的效果[16]。通過該方法得到的成礦能量，只能在一定程度上反映元素遷移和富集到當(dāng)前狀態(tài)下所變化的相對值，不能夠反映在這個(gè)過程中成礦能量所消耗的絕對值,但能夠反映成礦能量在當(dāng)前狀態(tài)下的分布情況，對于表示某一地區(qū)的地球化學(xué)元素的綜合信息具有一定的實(shí)際意義[17]。

圖3 惠井子地區(qū)地質(zhì)略圖 Fig.3 Geological map of Huijingzi area

通過spss軟件對新疆鄯善縣惠井子南山693件樣品的9種元素?cái)?shù)據(jù)做標(biāo)準(zhǔn)化處理后，得到各相關(guān)系數(shù)矩陣如表3所示。在0.01信度下與主攻元素Sb正相關(guān)的元素為Au、Zn、Pb、Mo、Bi五種元素 ，根據(jù)這些元素的背景值，按成礦能量計(jì)算式(1)求出該研究區(qū)的成礦能量，以此表示研究區(qū)的綜合異常信息。

表3 微量元素的相關(guān)系數(shù)

注：n=1269，γ(0.05)=0.052，γ(0.01)=0.081

4.2 泛克里金法的原理

泛克里金法數(shù)據(jù)處理是區(qū)域化變量非平穩(wěn)條件下的一種最優(yōu)化無偏線性估值方法[18]。利用泛克里金法進(jìn)行數(shù)據(jù)處理，其優(yōu)勢主要有：①可以進(jìn)行無偏估計(jì)，也就是估計(jì)值與已知點(diǎn)的值一致；②能夠給出空間意義下的最小方差估計(jì)。泛克里金法通過變差函數(shù)和不同的漂移形式較為準(zhǔn)確地反映地質(zhì)變量在不同方向的變化程度，可以有效地度量區(qū)域化數(shù)據(jù)空間相關(guān)性和變異程度[7]。通過測定區(qū)域化變量可以有效的為泛克里金插值提供變異性的定量化支持。其表達(dá)式[19]為:

(2)

式中：h表示滯后距；N(h)表示在滯后距為h時(shí)所有樣本的個(gè)數(shù)。通常情況下變差函數(shù)需要用塊金值、方向、變程和基臺(tái)值這4個(gè)參數(shù)進(jìn)行表示。在進(jìn)行實(shí)驗(yàn)變差函數(shù)之后要從理論上選擇一個(gè)合適的變差函數(shù)模型。常見的變差函數(shù)模型有球狀模型、指數(shù)函數(shù)模型、高斯型、冪函數(shù)型等。其中球狀模型是應(yīng)用最廣泛的，從而保證泛克里金插值有解，達(dá)到最優(yōu)擬合。

4.3 研究區(qū)成礦能量統(tǒng)計(jì)特征分析

研究區(qū)的面積13 km2，根據(jù)《地質(zhì)礦產(chǎn)勘查測量規(guī)范》(GB/T18341-2001)的要求按照100 m×250 m的采樣網(wǎng)格采集693件樣品。筆者對Sb元素的原始值的均值加上3倍的方差剔除特異值，從而使計(jì)算出來的數(shù)據(jù)符合正態(tài)分布，得到成礦能量的統(tǒng)計(jì)參數(shù)見表4。研究區(qū)成礦能量的平均值為6.45，最高值505.8。

表4 成礦能量統(tǒng)計(jì)參數(shù)表與與變差擬合函數(shù)

在計(jì)算過程中采用取樣間距0.2 km作為步長大小,角度容差為22.5°，帶寬為0.1 km，并用球狀模型進(jìn)行擬合變差函數(shù)。計(jì)算在0°、45°、90°、135°四個(gè)方向上的實(shí)驗(yàn)變差函數(shù)計(jì)算值點(diǎn)圖，最后得到擬合后的變差函數(shù)參數(shù)(表4)。由表4可以看出，新疆惠井子南山地區(qū)的成礦能量在其變程范圍之內(nèi)呈現(xiàn)較強(qiáng)的向異性；塊金效應(yīng)<25%，說明成礦能量在變程范圍內(nèi)具有較強(qiáng)的相關(guān)性；成礦能量連續(xù)性最佳方向(最大變程方向)為320°，與惠井子南山地區(qū)內(nèi)的主斷裂帶F2斷裂的走向基本一致，表明成礦能量受F2深大斷裂控制，F(xiàn)2斷裂可能是Sb的容礦構(gòu)造。在變差函數(shù)中成礦能量連續(xù)性變化的最大變程為2.4 km;由此表明，在研究區(qū)內(nèi)成礦能量具有較好的連續(xù)性。

4.4 應(yīng)用泛克里金法對地球化學(xué)成礦能量異常的提取

泛克里金法不僅能夠?qū)崿F(xiàn)線性無偏和最小估計(jì)方差的最優(yōu)估計(jì)，而且還可以給出指標(biāo)定量對其預(yù)測精度進(jìn)行評價(jià)[7]，與此同時(shí)，運(yùn)用泛克里金法繪制的地球化學(xué)圖能夠揭示成礦元素區(qū)域礦化異常特征以及影響成礦元素含量空間變化的因素[25]。通常情況下要求平均標(biāo)準(zhǔn)誤差較小，標(biāo)準(zhǔn)化均方根誤差接近為“1”，如表5所示，泛克里金法達(dá)到了較好的預(yù)測精度。 泛克里金法獲取的區(qū)域Sb異常見圖4。

表5 預(yù)測精度參數(shù)

圖4 泛克里金法獲取的區(qū)域成礦能量異常圖Fig.4 The regional obtain of the ore-forming energy by universal Kriging

5 應(yīng)用多重分形濾波技術(shù)對異常提取

5.1 多重分形濾波技術(shù)的基本原理

地球化學(xué)場的分布具有不均勻性和偶然性，除此之外，地殼的物質(zhì)結(jié)構(gòu)在各向異性上具有自相似的特征。多重分形濾波技術(shù)能夠從多重地球化學(xué)背景中提取有效的致礦異常信息，主要包括低背景值中的弱異常和隱蔽異常，為了能夠有效地從不均勻性和偶然性引起的離群值中提取地球化學(xué)場的異常，減少由于異常值造成的區(qū)域成礦背景的奇異現(xiàn)象[20]，需要對研究區(qū)的成礦能量進(jìn)行空間濾波處理。多重分形濾波(S-A)方法可以根據(jù)地球化學(xué)場分布的自相似性(尺度不變性)，采用濾波技術(shù)將異常從背景中分離出來[20-22]，解決了在復(fù)雜地質(zhì)背景和疊加異常的難題，具有減少變化背景影響、突出局部異常的功能?？梢岳媚茏V密度S與之間的函數(shù)關(guān)系來表示[7]：

A(≥S)∝S-β

(3)

其中：S為能譜密度；A(≥S)是能譜密度大于S的某一臨界值的區(qū)域面積；β為分形模型的指數(shù)系數(shù)(分維數(shù))。若S與A(≥S)服從指數(shù)關(guān)系，對其用最小二乘法繪制的雙對數(shù)圖(lnS與lnA)之間存在線性關(guān)系。用雙對數(shù)圖繪制的不同直線線段代表了不同的分形關(guān)系，兩條相交直線交點(diǎn)的橫坐標(biāo)就是譜能量密度值，通過這個(gè)值可以確定分形濾波過程中的閾值(臨界值)，進(jìn)而構(gòu)建不同的濾波器，從而實(shí)現(xiàn)對能譜密度的分離。通常情況下濾波器可被構(gòu)造為為高通、低通與帶通三種[7]分別如式(4)～式(6)所示：

(4)

(5)

(6)

已經(jīng)證明，波譜能量密度與波譜頻率成反比關(guān)系[23]，顯然，GA(w)的能譜密度低于GB(w)。從這個(gè)方面來看，GA(w)對應(yīng)的是高頻的部分，GB(w)對應(yīng)的是低頻部分。因此，GA(w)為高通濾波器 (異常濾波器)，GB(w)為低通濾波器(背景濾波器)，GC(w)只保留在區(qū)間(S1，S0)的能譜成分，因此，GC(w)為帶通濾波器(局部異常)，然后再利用Fouricr進(jìn)行逆變換，就得到分解后的背景值和異常值。

5.2 應(yīng)用多重分形濾波法提取區(qū)域成礦能量異常

許多學(xué)者[23-24]通過研究發(fā)現(xiàn)地球化學(xué)元素的空間分布在空間域具有局部不均一性和各向異性，這些分布通常情況下服從多重分形的分布特征。由于成礦能量是利用礦區(qū)內(nèi)元素的含量值與元素的背景值的比值來反映區(qū)域成礦能量的變化，因此成礦能量也表現(xiàn)出一定的分形分布特征[17]。S-A法對數(shù)圖要求所有的值必須為正值，因此，要對計(jì)算過程中出現(xiàn)的負(fù)值進(jìn)行相應(yīng)的變換計(jì)算，從而滿足對數(shù)圖的需要[25]。

用S-A法獲取的新疆惠井子地區(qū)成礦能量的“能譜密度-面積”變換后的雙對數(shù)圖(圖5)，得到了不同的擬合線段D1、D2、D3。 由圖5可以確定，成礦能量擬合的3條直線線段的交點(diǎn)橫坐標(biāo)為：LgS1=1.23，LgS2=1.64，所對應(yīng)的成礦能量為S1=16.98，S2=43.65用以作為構(gòu)建濾波器的閾值。對于研究區(qū)而言SS0代表了區(qū)域異常。

圖5 成礦能量的lnA(>S)-lnS圖Fig.5 ln-ln plot showing relationshipsbetween power spectrum value S and area A(>S) on the ore-forming energy

6 結(jié)果與討論

通過對研究區(qū)內(nèi)693件樣品的9種微量元素的含量進(jìn)行分析可知，Sb元素可以作為研究區(qū)的主成礦元素。綜合該地區(qū)的地質(zhì)特征，利用泛克里金法獲得的區(qū)域成礦能量異常圖可知：成礦能量異常高值區(qū)主要集中在研究區(qū)的南部；從巖性上看其異常高值點(diǎn)主要分布在石炭紀(jì)安山玢巖(C2Jbtf)區(qū)。這充分的說明Sb元素集中分布在石炭紀(jì)的安山玢巖(C2Jbtf)區(qū),在構(gòu)造活動(dòng)較為強(qiáng)烈，熱液來源較為有利的情況下Sb元素極有可能形成礦化體。因此Sb元素可以作為該區(qū)域的主成礦元素，石炭紀(jì)安山玢巖(C2Jbtf)可以作為該區(qū)域?qū)ふ襍b礦化點(diǎn)的指示性巖體，進(jìn)一步的說明該地區(qū)Sb成礦受巖控的作用明顯。

通過泛克里金法獲得的區(qū)域成礦能量異常圖中可以看出，該區(qū)域的成礦能量異常在空間上的分布具有一定的不均勻性。其中，成礦能量高值的分布符合線性規(guī)律，從局部上看成礦能量分布比較密集。從走向上來看，與研究區(qū)內(nèi)主斷裂F2斷裂的走向一致，成北北西向分布。因此可以推測，該地區(qū)的成礦過程可能與研究區(qū)內(nèi)斷裂構(gòu)造活動(dòng)有關(guān)，具體來說，可能是華力西中期巖漿巖侵入到上石炭統(tǒng)居里得能組圍巖中，通過一系列的構(gòu)造熱液活動(dòng)萃取圍巖中的Sb元素，在有利的成礦部位富集。因此，研究區(qū)的成礦能量受裂控作用明顯。

通過對低通濾波法獲取的圖與泛克里金法獲取的區(qū)域異常圖對比分析可知: 通過低通濾波獲取的區(qū)域成礦能量的異常分布在空間分布上與泛克里金法獲得的異常總體上具有較強(qiáng)的相似性，低通濾波獲得的異常點(diǎn)明顯地要比泛克里金法獲得的異常點(diǎn)多，這說明低通濾波能夠通過低頻高能譜信號，屏蔽掉高頻低能譜的信號，從而更好地反映背景中的弱異常和隱蔽異常。泛克里金法可以獲得異常趨勢，低通濾波可以獲得詳細(xì)的異常點(diǎn)，二者所獲得的異常圖更加明顯的顯示了地球化學(xué)元素在研究區(qū)域上的分布趨勢和濃集特征，為下一步找礦工作提供了理論依據(jù)。

morpas3.0軟件繪制出低通和帶通濾波獲取的成礦能量的地球化學(xué)異常分布如圖6所示。通過帶通濾波獲得的成礦能量局部異常(圖6(b))與泛克里金法(圖3)和低通濾波獲得的異常(圖6(a))相比較，帶通濾波獲取的成礦能量異常范圍縮小，濃集中心明顯減少且局部突出。究其原因主要是帶通濾波屏蔽和過濾掉那些不符合正態(tài)分布的異常值和離群值。剩余的異常信息在某種程度上能夠反映元素在地質(zhì)過程中的貧化和富集規(guī)律，這一切往往都與特定的礦化過程有關(guān)，從而為研究礦床的走向和尋找礦床的遠(yuǎn)景區(qū)提供一定的指示。

7 結(jié)論

1)通過變差函數(shù)對該地區(qū)成礦能量的計(jì)算研究表明,新疆惠井子南山地區(qū)成礦能量的最大變程為2.4 km，成礦能量連續(xù)性最佳方向?yàn)?20°，與該區(qū)域內(nèi)主斷裂F2斷裂的走向一致，主要呈NNW向分布，其成礦能量可能與該斷裂有關(guān)，并且在其變程范圍內(nèi)連續(xù)性較好，結(jié)構(gòu)性較強(qiáng)。

2)由泛克里金法獲取的成礦能量圖表明成礦能量的異常高值中心主要集中在中石炭統(tǒng)的安山玢巖(C2Jbtf)出露區(qū)域附近，因此中石炭統(tǒng)的安山玢巖可作為該區(qū)域?qū)ふ襍b礦化點(diǎn)的指示性巖性特征，其分布的走向與區(qū)域內(nèi)F2斷裂的走向一致，說明其成礦過程可能受斷裂構(gòu)造活動(dòng)影響，因此該區(qū)域成礦能量異常特征具有明顯的巖控和裂控特征。運(yùn)用多重分形濾波法的低通濾波獲取的區(qū)域成礦能量的異常分布在空間上與泛克里金法獲得的異常總體上具有較強(qiáng)的相似性，但是通過低通濾波獲得的異常點(diǎn)更多。

圖6 S-A法獲得的成礦能量異常圖Fig.6 The ore-forming energy S-A method(a)成礦能量區(qū)域異常圖(低通濾波)；(b)成礦能量的背景異常圖(帶通濾波)

3)通過對泛克里金法和多重分形濾波(S-A)法獲取的成礦能量異常圖可知：多重分形濾波(S-A)法能夠提取成礦能量的區(qū)域異常和隱伏異?；蛘弑尘爸械娜醍惓：碗[蔽異常，提供更多深層次的有效信息，能夠?yàn)檫M(jìn)一步的勘查工作提供有效準(zhǔn)確的靶區(qū)。

[1] 王瑞廷, 毛景文.區(qū)域地球化學(xué)異常評價(jià)的現(xiàn)狀及其存在的問題[J].中國地質(zhì),2005, 32(1): 168-175. WANG T R,MAO J W. Facts and problems of the evaluation Facts and problems of the evaluation of regional geochemical anomalies[J]. Chinese Geology, 2015,32(1):168-1759.(In Chinese)

[2] 成秋明.空間自相似性與地球物理和地球化學(xué)場的分解方法[J].地球物理學(xué)進(jìn)展,2001，16(2):8-17. CHENG Q M.Spatial Self-similarity and geophysical and geochemical anomaly decomposition [J], Progress in Geophysics,2001，16(2):8-17. (In Chinese)

[3] 王學(xué)求. 礦產(chǎn)勘查地球化學(xué): 過去的成就和未來的挑戰(zhàn)[J].地學(xué)前緣,2003,10(1): 239-248. WANG X Q. Exploration geochemistry: past Achievements and future challenges [J]. GeoscienceFrontiers,2003,10(1):239-248.(In Chinese)

[4] 謝學(xué)錦, 任天祥, 奚小環(huán), 等.中國區(qū)域化探全國掃面計(jì)劃卅年[J].地球?qū)W報(bào)，2009，30(6): 700-716. XIE X J,REN T X,XI X H, et al.The Implementation of the Regional Geochemistry-National Reconnaissance Program (RGNR) in China in the Past Thirty Years[J]. Journal of earth, 2009，30(6): 700-716. (In Chinese)

[5] 成秋明.非線性成礦預(yù)測理論：多重分形奇異性-廣義自相似性-分形譜系模型與方法[J].地球科學(xué)-中國地質(zhì)大學(xué)學(xué)報(bào),2006,3(3):337-348. CHENG Q M.Singularity - Generalized Self-Similarity-FractalSpectrum(3S)Models[J].Earth Science-Journal of China University of Geosciences,2006,3(3):337-348.(In Chinese)

[6] 成秋明，張生元，左光仁，等.多重分形濾波方法和地球化學(xué)信息提取技術(shù)研究與進(jìn)展[J].地學(xué)前緣，2009,16(2):185-198. CHENG Q M,ZHANG S Y,ZUO G R, et al. Progress of multifractal filtering techniques and their applications in geochemical information extraction. [J]. Earth science Frontiers, 2009,16(2):185-198.(In Chinese)

[7] 黃靜寧，趙鵬大．東天山地區(qū)Cu-Au異常信息提取與評價(jià)[J]．地球科學(xué)-中國地質(zhì)大學(xué)學(xué)報(bào)，2011，36(2)：317-325. HUANG J N,ZHAO P D.Extraction and assessment of Cu-Au anomalies of the east Tianshan Ore-forming belt,northwestern China[J].Earth Science-Journal of China University of Geosciences , 2011，36(2)：317-325. (In Chinese)

[8] CHENG.Q.M.Multifratal modeling and spatial analysis with GIS:gold potential estimation in the Mitchell-Sulphurets area. northwestern British Columbia (Dissertation)[D].Ottawa,University of Ottawa,1994.

[9] CHENG.Q.M.Spatial and scaling modeling for geochemical anomaly separation[J]. Journal of Exploration Geochemistry,1999,65:175-194.

[10]CHENG.Q.M.A new model for quantifying anisotropic scale invariance and for decomposition of mixing patterns[J].Mathematical Geology,2004,36(3):345-360.

[11]CHENG.Q.M.Quantifying the generalized self-similarity of spatial patterns for mineral resources assessment[J]. Earth Science Journal of China University of Geosciences,2004,29(6):733-743.

[12]馬瑞士，王賜銀，葉尚夫，等.東天山構(gòu)造格架及地殼演化[M].南京：南京大學(xué)出版社，1993. MA R S,WANG C Y,YE S F, et, al. East Tianshan tectonic framework and crustal evolution[M].Nanjing: Nanjing university publishing house, 1993. (In Chinese)

[13]馬瑞士，舒良樹，孫家齊.東天山構(gòu)造演化與成礦[M].北京：地質(zhì)出版社，1997. MA R S,SHU l S,SUN J Q. East Tianshan tectonic framework and metallogeny [M]Beijing:Geological Publishing House,1997 .(In Chinese)

[14]劉訓(xùn)，董凱，辛存林.新疆鄯善縣惠井子地區(qū)巖石銅異常暈特征及評價(jià) [J].甘肅地質(zhì),2013,22(2):56-61. LIU X,DONG K,XIN C L.Characteristics of copper anomalifes halo in Huijingzi, Shanshan county,Xinjiang [J]. Gansu Geology, 2013,22(2):56-61. (In Chinese)

[15]李汶奎，楊員兵.云南省耿馬勐力簡鉛鋅礦區(qū)地球化學(xué)特征及找礦前景預(yù)測[J].物探化探計(jì)算技術(shù).2010，32(2):200-205. LI WQ,YANG X B.Geochemical Characteristics and metallogenic prognosis of Meng jian lead-zinc mine area in Geng Ma Yunnan Province[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2010，32(2):200-205.(In Chinese)

[16]楊利民,楊自安.成礦能量在化探綜合異常提取中的應(yīng)用研究[J].礦產(chǎn)與地,2008，22(1):74-77. YANG LM,YANG Z A. Study on the application of ore-forming energy to the extraction of geochemical composite anomalies[J]. Mineral Resources and Geology, 2008，22(1):74-77. (In Chinese)

[17]辛存林，廖建棟，董凱，等.新疆喀拉尕依沙地區(qū)地球化學(xué)異常信息的提取[J].西北大學(xué)學(xué)報(bào)(自然科學(xué)版)，2015，45(2)275-282. XIN C L,LIAO J D,DONG K,et al. Geochemical abnormal information extract in Kalagayisha Xinjiang[J].Journal of Northwest University(Natural Science Edition) 2015，45(2):275-282.(In Chinese)

[18]王世稱，陳永良，夏立顯.綜合信息礦產(chǎn)預(yù)測理論與方法[M].北京：科學(xué)出版社，2000. WANG S C,CHEN Y L,XIA L X. Comprehensive information mineral prediction theory and method[M].Beijing:Science Press2000.(In Chinese)

[19]趙鵬大.定量地學(xué)方法及應(yīng)用[M].北京:地質(zhì)出版社, 2004. ZHAO P D. Quantitative methods and application of science[M].Beijing:Geological Publishing House,2004.(In Chinese)

[20]成秋明.空間模式的廣義自相似性分析與礦產(chǎn)資源評價(jià)[J].地球科學(xué)-中國地質(zhì)大學(xué)學(xué)報(bào)，2004,29(6)：733-743. CHENG Q M. Quantifying the generalized self-similarity of spatial patterns for mineral resource assessment[J] Earth Science-Journal of China University of Geosciences. 2004,29(6)：733-743.(In Chinese)

[21]李慶謀，成秋明.分形奇異(特征)值分解方法與地球物理和地球化學(xué)異常重建[J]．地球科學(xué)-中國地質(zhì)大學(xué)學(xué)報(bào),2004,29(1):109-118. LI Q M,CHENG Q M. Fractal Singular-Value (egin-value) decomposition method for geophysical and geochemical anomaly reconstruction[J]. Earth Science-Journal of China University of Geosciences,2004,29(1):109-118. (In Chinese)

[22]成秋明．多維分形理論和地球化學(xué)元素分布規(guī)律[J]．地球科學(xué)-中國地質(zhì)大學(xué)學(xué)報(bào),2000，25(3)：311-318． CHENG Q M. Multifractal theory and geochemical element distribution pattern[J]．Earth Science-Journal of China University of Geosciences,2000，25(3)：311-318.(In Chinese)

[23]趙鵬大，胡旺亮，李紫金.礦床統(tǒng)計(jì)預(yù)測[M].北京:地質(zhì)出版社, 1994. ZHAO P D,HU W L,LI Z J. Deposit statistical forecasting[M].Beijing:Geological publishing house,1994.(In Chinese)

[24]連永牢，胡天星，邵長來,等.黑龍江省遜克縣高松山金礦床微量元素地球化學(xué)[J].地質(zhì)與資源，2010,19(4)：287-291. LIAN Y L,HU X,SHAO C L, et al. Trace element geochemistry of the Gaosongshan gold deposit in Xunke County HeilongJiang[J].Geology and Resources,2010,19(4):287-291.(In Chinese)

[25]戴塔根，龔玲蘭，張起鉆.應(yīng)用地球化學(xué)[M].長沙:中南大學(xué)出版社，2005. DAI T G , GONG L L,ZHANG QZ. Applied geochemistry[M].Changsha:CentralSouthUniversityPress,2005.(In Chinese)

Geochemical abnormal information recognition and comprehensive evaluation in Huijingzi south mountain area Xinjiang

WANG Le， XIN Cunlin， ZHAO Zhen， LIAO Jiandong

(College of Geography and Environmental Science，Northwest Normal University，Lanzhou 730070，China)

This paper aims to identify, analysis and evaluate the comprehensive abnormal information of mineralization energy in Nanshan mountain area of Huijingzi, Xinjiang province with the methods of universal Kriging and multifractal filtering(S-A), so as to get comprehensive abnormal information and ore-controlled factors at different levels. The results of this research are as follows: ①Huijingzi，Xinjiang mineralization can be one of the biggest variation range of 2. 4 km，the direction of the maximum codomain is 320°, basically according with the strike of the main faults F1 in the study area. Its mineralization energy within range can show strong structural features．②The mineralization energy's abnormal information in space distribution we got through low pass filtering method is fairly similar to what we got through universal Kriging. The former one has distinct features in rock and crack control. ③The spectrum-area (S-A) method can extract the multi-level ore anomalies information，including the regional anomalies and hidden anomalies. The mineralization can also concealed in the background of weak anomalies．All these information at different levels got through multiple filtering method is supplied as new target area for the next research.

Huijingzi Xinjiang; the ore-forming energy; universal Kriging method; spectrum-area( S-A) method

2016-03-14 改回日期：2016-06-27

國家自然科學(xué)基金(41262001);甘肅省科技支撐基金(1104FKCA116)

王樂(1992-)，男，碩士，主要研究方向?yàn)榈刭|(zhì)礦產(chǎn)勘查，E-mail：1042796924@qq.com。

1001-1749(2017)02-0282-09

P 632

A

10.3969/j.issn.1001-1749.2017.02.20