辛存林，廖建棟，董 凱，包小強(qiáng)，謝慶豪，趙 振，王 樂，季健健

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

?

·地球科學(xué)·

新疆喀拉尕依沙地區(qū)地球化學(xué)異常信息的提取

辛存林，廖建棟，董 凱，包小強(qiáng)，謝慶豪，趙 振，王 樂，季健健

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

通過泛克里金法和多重分形濾波技術(shù)相結(jié)合的方法,對(duì)新疆喀拉尕依沙地區(qū)成礦能量異常信息進(jìn)行分析,以此獲取該地區(qū)不同層次的致礦異常信息。研究結(jié)果表明:①新疆喀拉尕依沙地區(qū)成礦能量的最大變程為2.4 km,最大變程的方向?yàn)?10°,最大變程方向與研究區(qū)主斷裂F1走向基本一致；該地區(qū)成礦能量變化主要受到區(qū)域內(nèi)的主斷裂F1控制,其成礦能量含量的變化表現(xiàn)出較強(qiáng)的結(jié)構(gòu)性特征。②泛克里金法獲得的異常與低通濾波獲取的區(qū)域成礦能量的異?？傮w具有相似的空間分布特征,它們均反映了異常的分布趨勢(shì)和濃集特征,具有明顯的裂控-巖控特征。③多重分形濾波(S-A)法能夠提取多層次的致礦異常信息,其中包括了成礦能量的區(qū)域異常和隱伏背景中的弱異常和隱蔽異常,通過對(duì)背景和異常信息的關(guān)聯(lián)分析,得出有利于礦化的解釋,指示隱伏礦床的分布位置，同時(shí)為尋求未知礦床提供新的啟示。

成礦能量;泛克里金法;多重分形濾波;新疆喀拉尕依沙

隨著勘探地球化學(xué)的發(fā)展，為了更加準(zhǔn)確地提取異常信息和對(duì)異常進(jìn)行評(píng)價(jià),就要求原始數(shù)據(jù)測(cè)量點(diǎn)必須均勻分布。但是，由于取樣條件、地質(zhì)情況、人為等因素的制約,不可避免地出現(xiàn)數(shù)據(jù)局部集中或過度分散,從而引起元素空間異常信息的畸變。而成礦能量反映了全部所查明元素含量的綜合信息,并且可對(duì)區(qū)域內(nèi)集中或分散的元素自動(dòng)篩選[1],從而使提取的信息減少了外界環(huán)境的影響,使得研究結(jié)果更加科學(xué)、有效。

地球化學(xué)元素的空間變化不僅具有結(jié)構(gòu)性的變化特征,還具有隨機(jī)性的變化特征[2]。這種隨機(jī)性變化往往是由致礦異常因素所引起的,其尺度不變性特征可由分形和多重分形模型定量表征[2-3]。成礦能量是利用與研究區(qū)主攻元素正相關(guān)的某幾種元素的含量同區(qū)域上該元素的背景值的比值來反映元素能量的變化,所以也表現(xiàn)出分形和多重分形的特征。基于非線性模型的多重分形濾波技術(shù),根據(jù)區(qū)域地球化學(xué)數(shù)據(jù)在頻率域中具有不同的廣義自相似性特征將其分解為不同的空間模式,以此來識(shí)別或分解疊加地質(zhì)作用引起的與礦化有關(guān)的異常[4-7]。多重濾波技術(shù)是通過傅里葉變化把空間域轉(zhuǎn)變成時(shí)空域,通過構(gòu)建濾波器,在頻率域上運(yùn)用不同的濾波器,并將其轉(zhuǎn)換在空間域上保持或加強(qiáng)區(qū)域的信息,對(duì)低背景中的弱異常或者隱蔽異常進(jìn)行提取、識(shí)別和圈定,還可以將異常區(qū)域的背景和異常進(jìn)行分離。

新疆喀拉尕依沙地區(qū)地形條件復(fù)雜,在采集樣品時(shí)，樣品間距的不均一性導(dǎo)致研究區(qū)元素化探數(shù)據(jù)異常信息的不準(zhǔn)確。因此，本文采用成礦能量法綜合分析研究區(qū)元素異常,以減小樣品間距不均勻?qū)е碌臄?shù)據(jù)信息的不準(zhǔn)確性。通過泛克里金法和多重分形濾波技術(shù)分別對(duì)研究區(qū)成礦能量空間信息的變化進(jìn)行分析,獲取研究區(qū)成礦能量異常信息的空間變化特征及其控制因素,并對(duì)研究區(qū)綜合信息進(jìn)行弱異常和隱蔽異常信息的提取。最后，結(jié)合研究區(qū)地質(zhì)背景,對(duì)研究區(qū)異常信息進(jìn)行綜合評(píng)價(jià)。

1 研究區(qū)成礦地質(zhì)背景

東天山處于準(zhǔn)噶爾板塊與塔里木板塊接合處,在長(zhǎng)期的演化過程中經(jīng)歷了復(fù)雜的裂解和拼合,地質(zhì)構(gòu)造復(fù)雜,背景獨(dú)特,巖漿活動(dòng)頻繁,是中國(guó)主要的多金屬成礦帶,具有復(fù)雜的地質(zhì)成礦背景[8-9]。研究區(qū)位于東天山中段準(zhǔn)噶爾微板塊(Ⅱ)、北天山造山帶(Ⅵ)和東準(zhǔn)噶爾板塊(Ⅳ)的交匯處,地殼活動(dòng)活躍(見圖1)。區(qū)域上位于博格達(dá)銅、金成礦帶,該成礦帶成礦地質(zhì)條件良好,極具找礦潛力。

研究區(qū)出露的地層主要為古生界的石炭系以及中生界的三疊系和新生界的第四系。區(qū)內(nèi)巖漿作用強(qiáng)烈,主要為華力西中期輝綠玢巖(βμ2a4),侵位于中石炭統(tǒng)居里得能組地層中。區(qū)內(nèi)斷裂較為發(fā)育,呈NNW,NEE向展布,斷裂在走向上具波狀彎曲,穿插于石炭系地層內(nèi),對(duì)區(qū)內(nèi)成礦物質(zhì)的活化、遷移和富集成礦均有重要的控制作用(見圖2)。在沿?cái)嗔寻l(fā)育的一系列巖漿巖中,常見黃銅礦、輝銅礦、硅化、褐鐵礦化、絹云母化、孔雀綠泥石蝕變等礦化蝕變。

通過對(duì)研究區(qū)采集的1 269個(gè)樣品元素含量的均值、標(biāo)準(zhǔn)離差、變化系數(shù)和濃度克拉克值等參數(shù)進(jìn)行統(tǒng)計(jì)分析(表1)后表明，其中濃度克拉克值大于1的元素有Cu，Pb，Zn，Bi，這說明這4種元素具有相對(duì)富集和較強(qiáng)富集；變異系數(shù)大于1的元素有Ag，Cu，Mo，Sb,這說明這4種元素活化和遷移作用強(qiáng)烈。從上述特征可以看出,在該區(qū)地層中,既有同生聚集(濃度克拉克值大于1),又具強(qiáng)烈后生改造(變異系數(shù)大于l)的元素(Cu元素), 其他元素的濃度克拉克值和變異系數(shù)均較低,其活化和遷移作用也較為微弱。這說明Cu元素不僅在地層中相對(duì)富集,而且分布極不均勻,在有利的地層條件下,極有可能形成礦化體。由野外地質(zhì)調(diào)查可知，研究區(qū)分布有6條近東西向展布的綠泥石化蝕變帶,其普遍具孔雀綠泥石蝕變,巖石中可見斑點(diǎn)狀黃銅礦、輝銅礦等。綜上可知，Cu元素為本研究區(qū)成礦能量的主攻元素。

圖1 新疆北部銅、金礦床分布圖[10-11]Fig.1 Sketch map showing the major tectonic units and distribution of copper and gold deposits in North Xinjiang

表1 研究區(qū)微量元素含量特征統(tǒng)計(jì)表Tab.1 Statistical table of trace element in researched area

2 泛克里金法對(duì)成礦能量綜合異常的提取

2.1 成礦能量法的原理與計(jì)算

成礦能量是指成礦元素由分散狀態(tài)組合成富集體(礦體或礦石)時(shí)的地球自然能,也就是各種金屬成礦元素從最初的分散狀態(tài)到局部富集成礦、成暈所需要消耗的相應(yīng)的能量。具體到一個(gè)研究區(qū)就是利用與主攻元素正相關(guān)的某幾種元素的含量同區(qū)域上該元素的背景值的比值來反映能量的變化[12]。單位體積礦石(暈)上成礦能量消耗(En)的計(jì)算公式[13]為

(1)

式中:En表示由n個(gè)元素形成單位體積礦石或地球化學(xué)暈時(shí)所消耗的能量;n為成礦、成暈的元素?cái)?shù);ki為組成礦石或地球化學(xué)暈的第i個(gè)元素的克拉克濃度值。成礦能量只是相對(duì)地反映了元素遷移和富集增加或者減小的能量,不反映元素富集或者分散的絕對(duì)值和具體到那一時(shí)期地質(zhì)作用對(duì)能量分布的影響,僅僅是最終形成狀態(tài)的能量的分布。

圖2 研究區(qū)地質(zhì)略圖Fig.2 Geological map of researched area

通過Spss軟件對(duì)新疆喀拉尕依沙地區(qū)1 269個(gè)樣品的9種元素?cái)?shù)據(jù)作標(biāo)準(zhǔn)化處理后,計(jì)算得到各相關(guān)系數(shù)矩陣見表 2。與主攻元素Cu正相關(guān)的元素為Ag，Mo，Zn，Au， Mn， Sb，Pb (在1%信度下的顯著相關(guān)),根據(jù)這些元素的背景值,按成礦能量計(jì)算公式(1)求出該研究區(qū)的成礦能量,以此表示研究區(qū)的綜合異常。

2.2 泛克里金法的原理

泛克里金法主要是運(yùn)用變差函數(shù)和非平穩(wěn)隨機(jī)函數(shù)的協(xié)方差函數(shù)為已知條件下，是一種考慮漂移的無偏線性估計(jì)量的地質(zhì)統(tǒng)計(jì)方法。通過變差函數(shù)和不同方向的漂移形式可以準(zhǔn)確地反映地質(zhì)變量在不同方向的變化程度。

變差函數(shù)是地質(zhì)統(tǒng)計(jì)學(xué)中描述區(qū)域化變量空間結(jié)構(gòu)性和隨機(jī)性的基本工具。它通過測(cè)定區(qū)域化變量分隔等距離樣點(diǎn)間的差異來研究區(qū)域化變量的空間相關(guān)性和空間結(jié)構(gòu)。而基于變差函數(shù)的泛克里金法繪制的地球化學(xué)異常圖可以為區(qū)域找礦提供信息。其表達(dá)式[9]為

(2)

式(2)中，N(h)為滯后距為h時(shí)參加變差函數(shù)計(jì)算的樣品個(gè)數(shù),h為滯后距。變差函數(shù)由4個(gè)參數(shù)來表征:方向、塊金值、基臺(tái)值和變程。在做出實(shí)驗(yàn)變差函數(shù)之后,選擇一個(gè)最優(yōu)的理論變差函數(shù)模型，即球狀模型進(jìn)行最優(yōu)擬合，以保證克里格方程有解。

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

研究區(qū)占地面積29.25 km2。根據(jù)100 m×250 m的網(wǎng)度組合一個(gè)樣品,共采集了1 269個(gè)樣品。對(duì)這1 269個(gè)樣品的主成礦元素(Cu)的成礦能量參數(shù)進(jìn)行統(tǒng)計(jì)后可知(見表3)，研究區(qū)內(nèi)成礦能量的平均值為0.04,最高值為86.52。本文以原始數(shù)據(jù)的平均值加上3倍的標(biāo)準(zhǔn)差剔除特異值,使數(shù)據(jù)接近正態(tài)分布。對(duì)研究區(qū)特異值剔除之后,用實(shí)驗(yàn)變差函數(shù)計(jì)算成礦能量的變化特征。

在計(jì)算過程中采用兩倍的取樣間距0.2 km作為步長(zhǎng),角度容差22.5°,帶寬為0.1 km,用球狀模型擬合變差函數(shù)。首先計(jì)算其在0°，45°，90°，135°四個(gè)方向上的實(shí)驗(yàn)變差函數(shù)值點(diǎn)圖,擬合出的變差函數(shù)參數(shù)見表3所示。研究區(qū)的成礦能量在變程之內(nèi)有較強(qiáng)的向異性,其基臺(tái)值為塊金值的兩倍以上,這說明成礦能量在變程范圍內(nèi)具有較強(qiáng)的相關(guān)性；成礦能量連續(xù)性最佳方向?yàn)?10°,這與區(qū)域內(nèi)的主斷裂帶F1斷裂一致,這表明成礦能量受到北北西向斷裂控制,也與研究區(qū)內(nèi)華力西中期輝綠玢巖(βμ2a4)的分布較為一致。在變差函數(shù)中表現(xiàn)為成礦能量連續(xù)性變化的最大變程和最小變程分別為2.4 km和1.03 km,這表明成礦能量的變化在研究區(qū)巖漿巖的范圍之內(nèi)具有較好的連續(xù)性。同時(shí)，其塊金效應(yīng)(塊金值/基臺(tái)值)比較高，這表明成礦能量表現(xiàn)出更強(qiáng)的結(jié)構(gòu)性特征。

表2 微量元素的相關(guān)系數(shù)Tab.2 Correlation coefficient matrix of microelements

表3 成礦能量統(tǒng)計(jì)參數(shù)與數(shù)據(jù)變差函數(shù)擬合參數(shù)Tab.3 Statistical parameter and fitting semi-variogram parameters of the ore-forming energy concentration

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

根據(jù)以上參數(shù)擬合的變差函數(shù),采用泛克里金法對(duì)研究區(qū)進(jìn)行地球化學(xué)成礦能量異常識(shí)別分析,獲取區(qū)域成礦能量異常圖(見圖3)。

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

泛克里金法的優(yōu)勢(shì)不僅是線性無偏和最小估計(jì)方差的最優(yōu)估計(jì),還可以給出指標(biāo)定量評(píng)價(jià)其預(yù)測(cè)精度[14]。該方法一般要求平均誤差絕對(duì)值與標(biāo)準(zhǔn)化平均誤差接近于0,平均標(biāo)準(zhǔn)誤差較小,標(biāo)準(zhǔn)化均方根誤差接近為1。由表4可知,在本研究中，泛克里金法達(dá)到了較好的預(yù)測(cè)精度。

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

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

為了有效地從不均勻性、偶發(fā)性變化引起的離群值中提取地球化學(xué)場(chǎng)的異常,減少特高與特低的異常值造成的背景的奇異現(xiàn)象, 本研究對(duì)研究區(qū)成礦能量進(jìn)行空間濾波處理。多重分形濾波(S-A)法,是空間分析和頻譜分析的集成, 主要是利用頻域中明顯的廣義自相似性,根據(jù)場(chǎng)的能譜分布的不同,用濾波技術(shù)選取頻率信息,重新恢復(fù)空間模式，將異常從背景中分離出來[5-7]。簡(jiǎn)單來說，就是為了突顯某些成分并且壓制其他成分的影響。我們可以利用能譜密度S與A之間的函數(shù)關(guān)系表示

A(≥S)∝S-β，

(3)

其中:S為能譜密度,A(≥S)是能譜密度大于S的波數(shù)集合的面積,β為冪指數(shù)。對(duì)式(3)兩邊同時(shí)取對(duì)數(shù)，繪制雙對(duì)數(shù)圖(lnS與lnA)，并用最小二乘法對(duì)其進(jìn)行分段擬合成一條直線。不同段擬合的直線的交點(diǎn)對(duì)應(yīng)的橫坐標(biāo)值為能譜密度，用它來確定分形濾波器的閾值。利用這個(gè)閾值構(gòu)建不同的濾波器,從而達(dá)到對(duì)場(chǎng)的分離。濾波器一般分為高通、低通與帶通3種類型[9]，分別如下式所示:

(4)

(5)

(6)

由于能譜密度與波譜頻率成反比關(guān)系,所以GA(w)為高通濾波器，也稱為異常濾波器,它能濾掉能譜密度大于S0能譜成分,保留小于S0能譜成分。GB(w)為低通濾波器，也稱為背景濾波器,它能濾掉能譜密度小于S1能譜成分,保留大于S1能譜成分。GC(w)為帶通濾波器，也稱為局部異常,它只保留在區(qū)間(S1,S0)的能譜成分，再根據(jù)逆傅里葉變換就可以得到分解后的背景值和異常值的數(shù)據(jù)。

表4 預(yù)測(cè)精度參數(shù)Tab.4 Precision parameter of prediction

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

成礦能量是利用與主攻元素正相關(guān)的某幾種元素的含量同區(qū)域上該元素的背景值(區(qū)域背景值為一個(gè)定值)的比值來反映能量的變化的,所以也表現(xiàn)出分形分布的特征。

由于成礦能量計(jì)算出來有負(fù)值(成礦能量為負(fù)值只反映元素富集到現(xiàn)有程度能量增加或減少的相對(duì)值,而不反映元素富集或分散所消耗的能量的絕對(duì)值)不符合擬合能譜密度S(spectrum density)-面積A(area)雙對(duì)數(shù)圖的要求,故對(duì)成礦能量進(jìn)行變換,在其成礦能量上加一個(gè)正整數(shù)C,使得該地區(qū)成礦能量與正整數(shù)C之和Eni>0,使得成礦能量變換后符合指數(shù)模型的運(yùn)算要求[9]。

用S-A法繪制研究區(qū)成礦能量“能譜密度-面積”關(guān)系雙對(duì)數(shù)圖(見圖4)，用最小二乘法擬合不同斜率的直線，這些直線代表了不同的自相似性特征，它們通常對(duì)應(yīng)空間域中不同級(jí)次的地球化學(xué)異常模式。如圖4所示，由成礦能量雙對(duì)數(shù)圖擬合了3條斜率不同的直線，分別為：k1=-1.551，k2=-2.632，k3=-5.047，3條直線的交點(diǎn)橫坐標(biāo)為：S1=1.775，S0=3.044，其所對(duì)應(yīng)的成礦能量為r1=5.899，r0=20.982，用此可作為構(gòu)建濾波器的閾值。對(duì)于研究區(qū)成礦能量的分布而言，SS0代表了區(qū)域異常。

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

4 結(jié)果與討論

1)通過對(duì)研究區(qū)采集的1 269個(gè)樣品的元素含量進(jìn)行統(tǒng)計(jì),并且結(jié)合野外地質(zhì)調(diào)查綜合分析,確定Cu元素為研究區(qū)成礦能量的主攻元素。由研究區(qū)成礦巖層分布特征和元素成礦能量異常識(shí)別分析可知,成礦能量異常高值中心主要集中在華力西中期輝綠玢巖(βμ2a4)出露區(qū)及其附近。Cu元素礦化異常主要在華力西中期輝綠玢巖(βμ2a4)和中石炭統(tǒng)居里得能組:灰綠—灰褐色安山巖中。由泛克里金法獲取的成礦能量異常表明，該成礦能量的異常高值點(diǎn)與華力西中期輝綠玢巖(βμ2a4)分布吻合程度較高。這說明Cu元素不僅在華力西中期輝綠玢巖中相對(duì)富集,而且分布極不均勻,在有利的構(gòu)造熱液條件下,極有可能形成礦化體。華力西中期輝綠玢巖(βμ2a4)可作為該區(qū)域?qū)ふ褻u礦化點(diǎn)的指示性巖性特征,因此該地區(qū)銅成礦具有明顯的巖控特征。

2)從泛克里金法獲取的成礦能量異常分布和走向趨勢(shì)上來看,其異常具有明顯的濃度分帶特征,異常主要集中在研究區(qū)內(nèi)的最大斷裂F1的北側(cè),受到F1主斷裂和F3，F(xiàn)4次級(jí)斷裂的控制,呈北北西走向的帶狀分布,分布并不連續(xù)。由此可以推測(cè)，其成礦過程可能跟斷裂構(gòu)造有關(guān)的巖漿熱液作用具有一定的關(guān)聯(lián),可能是華力西中期巖漿巖侵位于中石炭統(tǒng)居里得能組圍巖中,通過熱液萃取圍巖中活化、遷移能力強(qiáng)的Cu元素到有利的成礦部位富集成礦。該異?？勺鳛閷ふ覠嵋盒豌~多金屬礦產(chǎn)的重要遠(yuǎn)景區(qū)。因此，該成礦過程既有明顯的巖控又具有明顯的裂控特征。

3)通過對(duì)低通濾波(見圖5A)獲取的異常與泛克里金法獲取的異常對(duì)比分析可知:① 低通濾波獲取的區(qū)域成礦能量的異常與泛克里金法獲得的異常總體具有相似的空間分布特征。② 泛克里金法所獲得的異常在區(qū)域上具有較明顯的分布趨勢(shì)和濃集特征,這種特征往往反映的是區(qū)域背景下的重大致礦地質(zhì)事件。③ 低通濾波獲得的異常點(diǎn)明顯比泛克里金法獲得的多，這說明低通濾波突出了背景中的弱異常和隱蔽異常。泛克里金法可以獲得異常點(diǎn)和異常趨勢(shì),低通濾波可以獲得詳細(xì)的異常點(diǎn)。它們兩者結(jié)合為下一步探礦工作指明了方向和依據(jù)。

4)通過帶通濾波獲得的局部異常(見圖5B)與泛克里金法和低通濾波獲得的異常(見圖5A)相比,異常范圍縮小并且突出,其局部異常出現(xiàn)了更強(qiáng)的方向性。這是因?yàn)橥ㄟ^帶通濾波對(duì)那些次生富集的不均勻性和偶發(fā)性引起的離群值與周圍的樣品相關(guān)性較差的數(shù)據(jù)已被屏蔽掉。而這些局部均勻的異常信息反映了元素在地質(zhì)過程中的運(yùn)移、富集規(guī)律,與特定的礦化過程有關(guān),為研究礦床的走向和尋找遠(yuǎn)景區(qū)礦床有一定的指示作用。

圖5 S-A法獲得的成礦能量異常圖Fig.5 The ore-forming energy S-A method

5 結(jié) 論

1)通過對(duì)成礦能量運(yùn)用變差函數(shù)的研究表明，新疆喀拉尕依沙地區(qū)成礦能量的最大變程為2.4 km,最大變程的方向?yàn)?10°。這表明該地區(qū)成礦能量變化主要受到區(qū)域內(nèi)的NNW主斷裂F1斷裂控制,具有較好的連續(xù)性。其比較高的塊金效應(yīng)(塊金值/基臺(tái)值)表明，它的成礦能量表現(xiàn)出較強(qiáng)的結(jié)構(gòu)性特征。

2)泛克里金法獲取的成礦能量異常反映了異常的分布趨勢(shì)和濃集特征,具有明顯的裂控和巖控特征。低通濾波獲取的區(qū)域成礦能量的異常與泛克里金法獲得的異常空間分布總體具有相似的空間分布特征,但是，低通濾波獲取的異常明顯比泛克里金法獲得的異常點(diǎn)多。

3)多重分形濾波(S-A)法能夠提取多層次的致礦異常信息,其中包括了成礦能量的區(qū)域異常和隱伏的或者背景中的弱異常和隱蔽異常。通過多重濾波方法提取的不同層次的異常為下一步工作提供了新的研究靶區(qū)。通過對(duì)背景和異常的異常信息的關(guān)聯(lián)分析,得出有利于礦化的解釋,指示隱伏礦床的分布位置，并為尋求未知礦床提供新的啟示。

[1] 楊利民,楊自安.成礦能量在化探綜合異常提取中的應(yīng)用研究[J].礦產(chǎn)與地,2008,22(1):74-77.

[2] 成秋明.非線性成礦預(yù)測(cè)理論:多重分形奇異性-廣義自相似性-分形譜系模型與方法[J].地球科學(xué)：中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2006,3(3):337-348.

[3] LIMA A, PLANT J A, DE VIVO B,et al. Interpolation methods for geochemical maps: Acomparative study using arsenic data from European streamwaters[J]. Geochem Explor Environ Anal, 2008, 8(1): 41-48.

[4] 成秋明, 張生元,左仁廣, 等. 多重分形濾波方法和地球化學(xué)信息提取技術(shù)研究與進(jìn)展[J]. 地學(xué)前緣,2009,16(2): 185-198.

[5] 成秋明.空間模式的廣義自相似性分析與礦產(chǎn)資源評(píng)價(jià)[J].地球科學(xué)：中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2004,29(6):733-743.

[6] 李慶謀,成秋明.分形奇異(特征)值分解方法與地球物理和地球化學(xué)異常重建[J].地球科學(xué)：中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2004,29(1):109-118.

[7] 成秋明.多維分形理論和地球化學(xué)元素分布規(guī)律[J].地球科學(xué)：中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2000,25(3):311-318.

[8] 崔彬, 和志軍, 趙磊, 等.新疆東天山中段金、銅成礦系統(tǒng)研究[J]. 地學(xué)前緣，2000,25(3):311-316.

[9] 黃靜寧,趙鵬大.東天山地區(qū)Cu-Au異常信息提取與評(píng)價(jià)[J].地球科學(xué)：中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2011,36(2):317-326.

[10] 秦克章.新疆北部中亞型造山帶與成礦作用[D].北京：中國(guó)科學(xué)院研究生院，2000.

[11] 韓春明,肖文交,崔彬,等.新疆北部晚古生代銅礦床主要類型和地質(zhì)特征[J].地質(zhì)學(xué)報(bào),2006,80(1):74-89.

[12] 楊利民,楊自安,羅鐵良,等.成礦能量在化探綜合異常提取中的應(yīng)用研究[J].礦產(chǎn)與地質(zhì),2008,22(1):74-77.

[13] 徐錫華.成礦能量-克立格雙模型在化探數(shù)據(jù)處理中的特殊功效[J].地質(zhì)找礦論叢,2002,17(1):52-56.

[14] 王景雷,孫景生,張寄陽(yáng),等.基于GIS和地質(zhì)統(tǒng)計(jì)學(xué)的作物需水量等值線圖[J].農(nóng)業(yè)工程學(xué)報(bào),2004,20(5):51-54.

(編 輯雷雁林)

Geochemical abnormal information extract in Kalagayisha Xinjiang

XIN Cun-lin, LIAO Jian-dong, DONG Kai, BAO Xiao-qiang, XIE Qing-hao, ZHAO Zhen, WANG Le, JI Jian-jian

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

Based on the universal Kriging method and Spectrum-Area(S-A) method, essential mineralization abnomal information in Kalagayisha areas, Xinjiang were analyzed, leadingto different levels of mine information.Results show that: ① Kalagayisha, Xinjiang mineralization can be one of the biggest variation range of 2.4 km, the direction of the maximum codomain is 310°, the largest variation direction basically accords with the strike of the main faults F1 in the study area, suggestting that the region mineralization can change is mainly controlled by the regional main faults F1 fault. Changes in concentrations of mineralization can show strong structural features. ② The anomaly obtained by the Universal Kriging Method and the regional ore-forming energy anomalies obtained by the low-pass filter have similar spatial distribution characteristics as a whole. They all reflect the trend and concentration distribution characteristics of abnormal information, with obvious crack and rock control characteristics. ③ The Spectrum-Area(S-A) method can extract the multi-level ore anomalies information, including the regional anomalies and mineralization can concealed in the background of weak anomalies. Through correlation analysis of the background and anomaly exception information, the mineralization explanation is reached, indicating the distribution of concealed deposit location and providing a new relation of seeking the unknown deposit.

the ore-forming energy; Universal Kriging Method; Spectrum-Area(S-A) method; Kalagayisha Xinjiang

2014-04-16

國(guó)家自然基金資助項(xiàng)目(41262001)；甘肅省科技支撐基金資助項(xiàng)目(1104FKCA116)；甘肅省自然科學(xué)基金資助項(xiàng)目 (1010RJZA020)；甘肅省高等學(xué)校研究生導(dǎo)師科研基金資助項(xiàng)目(1101-07)；蘭州市科技計(jì)劃基金資助項(xiàng)目(2011-2-89)

辛存林,男,甘肅秦安人,教授,從事地質(zhì)礦產(chǎn)與勘查研究。

P632

：ADOI：10.16152/j.cnki.xdxbzr.2015-02-020