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

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

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

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

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

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

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

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

0 引言

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

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

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

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

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

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

2 工作方法

2.1 采樣點的布置

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

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

2.2 樣品分析

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

表1 各元素分析方法

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

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

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

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

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

4.1 成礦能量的理論方法

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

(1)

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

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

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

表3 微量元素的相關系數

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

4.2 泛克里金法的原理

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

(2)

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

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

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

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

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

4.4 應用泛克里金法對地球化學成礦能量異常的提取

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

表5 預測精度參數

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

5 應用多重分形濾波技術對異常提取

5.1 多重分形濾波技術的基本原理

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

A(≥S)∝S-β

(3)

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

(4)

(5)

(6)

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

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

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

用S-A法獲取的新疆惠井子地區(qū)成礦能量的“能譜密度-面積”變換后的雙對數圖(圖5)，得到了不同的擬合線段D1、D2、D3。 由圖5可以確定，成礦能量擬合的3條直線線段的交點橫坐標為：LgS1=1.23，LgS2=1.64，所對應的成礦能量為S1=16.98，S2=43.65用以作為構建濾波器的閾值。對于研究區(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 結果與討論

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

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

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

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

7 結論

1)通過變差函數對該地區(qū)成礦能量的計算研究表明,新疆惠井子南山地區(qū)成礦能量的最大變程為2.4 km，成礦能量連續(xù)性最佳方向為320°，與該區(qū)域內主斷裂F2斷裂的走向一致，主要呈NNW向分布，其成礦能量可能與該斷裂有關，并且在其變程范圍內連續(xù)性較好，結構性較強。

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

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

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

[1] 王瑞廷, 毛景文.區(qū)域地球化學異常評價的現(xiàn)狀及其存在的問題[J].中國地質,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] 成秋明.空間自相似性與地球物理和地球化學場的分解方法[J].地球物理學進展,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] 王學求. 礦產勘查地球化學: 過去的成就和未來的挑戰(zhàn)[J].地學前緣,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] 謝學錦, 任天祥, 奚小環(huán), 等.中國區(qū)域化探全國掃面計劃卅年[J].地球學報，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] 成秋明.非線性成礦預測理論：多重分形奇異性-廣義自相似性-分形譜系模型與方法[J].地球科學-中國地質大學學報,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] 成秋明，張生元，左光仁，等.多重分形濾波方法和地球化學信息提取技術研究與進展[J].地學前緣，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異常信息提取與評價[J]．地球科學-中國地質大學學報，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]馬瑞士，王賜銀，葉尚夫，等.東天山構造格架及地殼演化[M].南京：南京大學出版社，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]馬瑞士，舒良樹，孫家齊.東天山構造演化與成礦[M].北京：地質出版社，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]劉訓，董凱，辛存林.新疆鄯善縣惠井子地區(qū)巖石銅異常暈特征及評價 [J].甘肅地質,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ū)地球化學特征及找礦前景預測[J].物探化探計算技術.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]楊利民,楊自安.成礦能量在化探綜合異常提取中的應用研究[J].礦產與地,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ū)地球化學異常信息的提取[J].西北大學學報(自然科學版)，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]王世稱，陳永良，夏立顯.綜合信息礦產預測理論與方法[M].北京：科學出版社，2000. WANG S C,CHEN Y L,XIA L X. Comprehensive information mineral prediction theory and method[M].Beijing:Science Press2000.(In Chinese)

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

[20]成秋明.空間模式的廣義自相似性分析與礦產資源評價[J].地球科學-中國地質大學學報，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]李慶謀，成秋明.分形奇異(特征)值分解方法與地球物理和地球化學異常重建[J]．地球科學-中國地質大學學報,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]成秋明．多維分形理論和地球化學元素分布規(guī)律[J]．地球科學-中國地質大學學報,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)計預測[M].北京:地質出版社, 1994. ZHAO P D,HU W L,LI Z J. Deposit statistical forecasting[M].Beijing:Geological publishing house,1994.(In Chinese)

[24]連永牢，胡天星，邵長來,等.黑龍江省遜克縣高松山金礦床微量元素地球化學[J].地質與資源，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]戴塔根，龔玲蘭，張起鉆.應用地球化學[M].長沙:中南大學出版社，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

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

王樂(1992-)，男，碩士，主要研究方向為地質礦產勘查，E-mail：1042796924@qq.com。

1001-1749(2017)02-0282-09

P 632

A

10.3969/j.issn.1001-1749.2017.02.20