范小軍，陳 沖，王曉剛，陳 根

(1.四川省核工業(yè)地質(zhì)調(diào)查院，四川 成都 610066;2.中國地質(zhì)大學(xué)(武漢)資源學(xué)院，湖北 武漢430074)

西北半干旱地區(qū)巖屑地球化學(xué)異常提取方法的研究

范小軍1，2，陳 沖2，王曉剛1，陳 根1

(1.四川省核工業(yè)地質(zhì)調(diào)查院，四川 成都 610066;2.中國地質(zhì)大學(xué)(武漢)資源學(xué)院，湖北 武漢430074)

西北地區(qū)主要發(fā)育殘山丘陵地貌。風(fēng)成砂、急雨過后的黃土硬殼以及下伏膏巖層的普遍發(fā)育在一定程度上影響了地球化學(xué)勘查數(shù)據(jù)的代表性，也給異常的提取和解釋帶來了困難。針對這種特殊的自然景觀和地質(zhì)背景，以西北半干旱地區(qū)——阿拉善某地為例，介紹了工作區(qū)地質(zhì)特征，用傳統(tǒng)的統(tǒng)計學(xué)方法、趨勢面法、滑動平均法提取相關(guān)元素地球化學(xué)異常，在此基礎(chǔ)上結(jié)合地質(zhì)礦產(chǎn)資料進行綜合對比，探討了較為適合半干旱地區(qū)地球化學(xué)勘查數(shù)據(jù)處理的方法。

巖屑地球化學(xué)測量;異常提取;趨勢面法;滑動平均法;西北半干旱地區(qū)

0 引言

我國西北地區(qū)屬于半干旱氣候條件，嚴重缺水，地質(zhì)研究程度低。由于本區(qū)蒸發(fā)強烈，往往在短暫的暴雨后，會在山坡的下部、溝谷以及第四系中形成10～30 cm的堅硬黃土硬殼，第四系中甚至可厚達1 m以上。黃土殼中也有少量的細小巖屑與風(fēng)成砂，但是巖屑顆粒外被堅硬的黃土皮殼包裹，有些粒度與下部的巖屑粒度相差不大。另外，下伏的膏鹽層在工作區(qū)也比較發(fā)育，形成的地球化學(xué)障能夠?qū)Τ傻V及伴生元素異常的顯示起屏蔽作用。

在前人研究成果的基礎(chǔ)上，結(jié)合研究區(qū)實際地理景觀條件，認為研究區(qū)不適宜水系沉積物地球化學(xué)測量，也不適宜基巖地球化學(xué)測量。因此，采用巖屑地球化學(xué)測量方法。

1 區(qū)域地質(zhì)背景

從大地構(gòu)造位置上看，研究區(qū)地處塔里木、哈薩克斯坦、西伯利亞和華北4大板塊匯聚地帶，屬于天山—興蒙造山系(一級構(gòu)造單元)額濟納—北山弧盆系(二級構(gòu)造單元)公婆泉島弧(三級構(gòu)造單位)構(gòu)造的西部(潘桂棠等，2009)。北山地區(qū)的主體是哈薩克斯坦板塊的東延部分，其北界位于北緯42°以北的駱駝山—紅石山—黑鷹山一線，南界位于柳園—大奇山一帶(聶鳳軍等，2003)。

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

研究區(qū)出露的巖性主要有中基性巖漿巖(輝長巖體、閃長巖體)，酸性巖漿巖(花崗巖體)，變質(zhì)巖區(qū)(變質(zhì)熔巖、石英片巖、石英巖)，沉積巖區(qū)(灰?guī)r、砂巖、礫巖)，第四系洪沖積物(圖1)。

2.1 地層

研究區(qū)出露的地層主要有志留系圓包山組石英片巖和公婆泉組變質(zhì)火山巖、石炭系綠條山組長英角巖及石英砂巖和岌岌臺子組碳酸鹽建造;中生代地層主要有侏羅系龍鳳山組礫砂巖、赤金堡組砂礫巖和粉砂巖及泥巖;第四系。

2.2 構(gòu)造

區(qū)內(nèi)斷裂構(gòu)造主要為北西向、北東、東西向或近東西向，有高角度的正斷層、逆斷層或逆沖斷層，規(guī)模大，活動時間長，繼承性明顯，對區(qū)內(nèi)地層的分布起著控制作用。褶皺構(gòu)造復(fù)雜，與華力西期侵入體有密切關(guān)系的銅、鐵礦化，多分布在褶皺構(gòu)造轉(zhuǎn)折處的斷裂或接觸帶附近。

圖1 研究區(qū)簡略地質(zhì)圖

2.3 巖漿巖

區(qū)內(nèi)侵入巖分布廣泛，占圖幅面積的80%以上，包括基性巖、中性巖、中酸性巖及酸性巖類，其中以中酸性巖類分布最廣，呈巖基、巖枝及巖脈狀產(chǎn)生。侵入活動集中于華力西期?；鹕綆r為早石炭世芨芨臺子組(Cj)中基性火山巖，志留紀公婆泉組(Sg)中酸性火山碎屑巖。區(qū)內(nèi)脈巖發(fā)育，其中石英脈中可見孔雀石化，局部含有金礦化。

3 異常提取

風(fēng)力是破壞巖石和搬運風(fēng)化產(chǎn)物的最重要地質(zhì)外營力，化探異常均不同程度受到風(fēng)成砂的干擾，風(fēng)的吹蝕、分選搬運和風(fēng)積物的堆積與混入，不僅可改變和雜化異常分布，而且可使異常強度降低甚至完全消失。研究區(qū)具有獨特的半干旱荒漠區(qū)戈壁殘山景觀，風(fēng)化剝蝕嚴重，風(fēng)成砂、急雨過后包裹細小巖屑的黃土硬殼以及下伏膏巖層的普遍存在，導(dǎo)致元素背景值發(fā)生變化，影響地球化學(xué)數(shù)據(jù)的真實性和異常信息的提取。針對這種特殊的自然景觀和地質(zhì)背景，采用不同的異常提取方法進行探討和分析。

3.1 經(jīng)典的統(tǒng)計學(xué)方法

傳統(tǒng)的統(tǒng)計學(xué)方法要求數(shù)據(jù)服從正態(tài)分布或者對數(shù)正態(tài)分布，從數(shù)據(jù)統(tǒng)計特征分析可以看出，Cu、Au、Ni服從類似對數(shù)正態(tài)分布，Co、Cr、Mn 服從類似正態(tài)分布，因此要分別對數(shù)據(jù)進行相應(yīng)的處理。

對Cu、Au、Ni 3種元素進行取對數(shù)，再用對數(shù)平均值加減3倍標準差作為數(shù)據(jù)上下限對奇異值進行迭代替換(表1)。

對Co、Cr、Mn 3種元素?zé)o需進行取對數(shù)變換，直接用原始數(shù)據(jù)的平均值加減3倍標準差作為數(shù)據(jù)上下限對奇異值進行迭代替換(表2)。

表1 Cu、Au、Ni取對數(shù)—剔值后的特征值統(tǒng)計結(jié)果

表2 Co、Cr、Mn剔值后的特征值統(tǒng)計結(jié)果

從表1、表2可見，處理后的數(shù)據(jù)已基本符合正態(tài)分布，可以用于異常的圈定。

單元異常下限的計算公式：T=C0+Ks0(1)式中，T為異常下限;C0為背景值;s0為標準偏差;K為可靠性系數(shù)。

研究區(qū)內(nèi)風(fēng)化剝蝕較為嚴重，巖屑中各元素含量有所變化，并且項目屬于1∶5萬大面積化探，為了不損失一切可能的礦化信息，Co、Cr、Mn異常下限取值應(yīng)相對降低，信度系數(shù)應(yīng)提高，選擇信度α =0.1對應(yīng)的 K=1.65;Cu、Au、Ni是取對數(shù)后服從正態(tài)分布的，取對數(shù)的做法大大提高了異常下限(蔣敬業(yè)，2006)，因此，信度系數(shù)應(yīng)再提高，選擇K=1.5。

確定了單元素異常下限值T，Cu、Au、Ni的T值要返回(T0=power(10，T))為 T0值，再以 T0的1.0、1.5、2.0倍作為異常的外帶、中帶、內(nèi)帶下限值(表3)。

表3 各元素的外、中、內(nèi)帶值

Co、Cr、Mn 直接以 T 值的 1.0、1.5、2.0 倍作為異常的外、中、內(nèi)帶下限值。按照內(nèi)帶、中帶、外帶下限值做出了單元素的異常圖(圖2)。

由圖2明顯可見，異常主要與變質(zhì)火山熔巖和中基性巖關(guān)系最為密切，主要分布在公婆泉組(Sg)地層、輝長巖及圓包山組(Sy)地層，且與地質(zhì)體形態(tài)吻合較好;Cu、Au、Ni、Co、Cr、Mn 在尖山幅輝長巖區(qū)附近面積比較大，異常有明顯的濃度分帶;Cu異常主要分布在公婆泉組(Sg)地層中，與變質(zhì)火山熔巖的關(guān)系最為密切;Au異常除了主要分布在公婆泉組(Sg)地層中外，其次分布在黑云母花崗巖、赤金堡組(Kc)及綠條山組(Cl)地層中，與變質(zhì)火山熔巖、黑云母花崗巖、砂巖關(guān)系密切;提取的異常與Cu、Au礦化點有一定的吻合度，主要是與尖山輝長巖附近礦化點吻合。

圖2 傳統(tǒng)統(tǒng)計學(xué)方法單元素異常圖

異常形態(tài)不規(guī)整，面積大小不一，形態(tài)以鋸齒狀、長條狀、尖棱狀為主;研究區(qū)內(nèi)北西向深大斷裂是一條沖積溝，兩側(cè)高山巖屑物質(zhì)向下流動，提取的異常(如 Cu、Au、Ni、Mn)有沿深大斷裂分布的趨勢，說明傳統(tǒng)法不能抑制地形地貌對異常提取的干擾;雖然提取出的異常與礦化點有一定的吻合度，但是很多重要的礦化點部位都沒有提取到異常，比如尖山幅西南角弱背景區(qū)內(nèi)重要的Cu礦化點D10295(25.99%)、D0172(1.28%)及其周圍也沒有提取到異常;圈出的異常主要分布在公婆泉組高背景區(qū)的地質(zhì)體中，同時這意味著已漏掉其他低背景區(qū)的弱緩異常。

因此，傳統(tǒng)方法提取的異常有自身的局限性，由于不同的地質(zhì)單元一般具有不同的地球化學(xué)背景，在分布有多個地質(zhì)單元的地區(qū)開展化探工作時，如果不考慮地質(zhì)單元的地球化學(xué)背景差異，而在全區(qū)采用統(tǒng)一的異常下限，勢必會出現(xiàn)這樣的情況：若異常下限選擇過低，則會造成大面積的異常分布;反之，若異常下限過高，則會掩蓋低背景區(qū)的異常(趙榮軍，2006)，這就決定了提取的異常主要分布在含量相對較高的背景區(qū)內(nèi)。

3.2 基于Surfer軟件的趨勢面法

3.2.1 數(shù)據(jù)預(yù)處理 由于原始數(shù)據(jù)中個別異常點的特高值會對趨勢面的形態(tài)產(chǎn)生很大的影響，因此各元素數(shù)據(jù)中高于95%分位數(shù)的值全部用95%分位數(shù)替換。

3.2.2 異常下限計算 在Surfer軟件運用多項式回歸方法——三次曲面將替換后的數(shù)據(jù)進行網(wǎng)格化，再計算原始數(shù)據(jù)與網(wǎng)格化數(shù)據(jù)的殘差值，以殘差值均值加減3倍的標準差作為上下限對殘差值進行剔值處理，并對其進行統(tǒng)計特征分析(表4)。

表4 殘差值剔值數(shù)據(jù)特征統(tǒng)計

從表4可以看出，剔值后剩余值的分布模式接近正態(tài)分布，總體上噪音(隨機分量)對數(shù)據(jù)基本沒有影響，可以用傳統(tǒng)方法圈定異常，異常下限取值時，選擇K=1.5;求出異常下限值后，用該值的1.0、2.0、4.0倍作為異常外帶、中帶、內(nèi)帶的下限值圈定異常(圖3)。

結(jié)合圖2及圖3可見，趨勢面法圈出的異常與傳統(tǒng)法圈出的異常大體上類似，與礦化點吻合程度差不多。但是相對于傳統(tǒng)方法圈出的異常，趨勢面法圈出的異常有以下幾個特點。

(1)總體上異常面積縮小，如基東幅與石樁子井幅交界處的大面積Cu、Mn，高背景區(qū)的異常被削弱，異常形態(tài)相比而言變得較規(guī)整。

(2)異常的形態(tài)有所變化，產(chǎn)于高背景區(qū)的連續(xù)大面積異常被分解成斷續(xù)異常，如產(chǎn)于公婆泉組(Sg)中的 Cu、Mn異常。

(3)相對高背景區(qū)的大面積異常在消弱分解，西林陶勒幅和尖山幅低背景區(qū)中小異常的面積大小不變，有的甚至變大，如西林陶勒幅下部Co異常。

趨勢面法能有效地削弱高背景異常，突出低背景區(qū)的低緩異常，并且操作比較簡單，不需要對原始數(shù)據(jù)作過多的處理，在一定程度上也降低了處理數(shù)據(jù)過程中產(chǎn)生的隨機誤差。

3.3 滑動平均法

正如趨勢面描述的地球化學(xué)背景場一樣，滑動平均法也認為地球化學(xué)背景面是一個連續(xù)起伏變化的曲面，異常提取的關(guān)鍵也是從區(qū)域趨勢變化中尋找局部富集部位。

首先將化探數(shù)據(jù)網(wǎng)格化，使每一個網(wǎng)格中只放置1個值，若落在該網(wǎng)格中的是數(shù)個點，則要求均值，并將求得的均值再置入該格的中央。

圖3 趨勢面法單元素異常圖

然后，根據(jù)地質(zhì)情況在網(wǎng)格化數(shù)據(jù)的基礎(chǔ)上開“窗口”，“窗口”中心的趨勢值是根據(jù)該“窗口”范圍內(nèi)的各網(wǎng)格交點的含量、分布特點及數(shù)學(xué)模型計算求得?；瑒悠骄ㄊ菍ⅰ按翱凇贝笮》秶鷥?nèi)已網(wǎng)格化的數(shù)據(jù)的平均值作為“窗口”中心的趨勢值。計算完一個“窗口”后，將“窗口”分別進行由行到列的滑動，又依次求出整個區(qū)內(nèi)各滑動“窗口”中心的趨勢值?！按翱凇钡倪x擇要結(jié)合地質(zhì)情況和取樣分析數(shù)據(jù)的大致情況來進行。通常認為，大窗口指示區(qū)域的概率比較大，用來模擬地球化學(xué)背景的變化。而小窗口更偏重突出一些局部異常，用來模擬地球化學(xué)局部異常的變化。最后，用大、小“窗口”模擬的變化趨勢的差值來圈定異常。

具體的操作方法是先將數(shù)據(jù)網(wǎng)格化，經(jīng)過選擇對比，確定用3×3滑動平均網(wǎng)格化數(shù)據(jù)得到網(wǎng)格文件1，在3×3滑動平均網(wǎng)格化的基礎(chǔ)上進行21×21滑動平均，得到網(wǎng)格文件2，最后用網(wǎng)格文件1減去網(wǎng)格文件2得到殘差值。對殘差值以其平均值加減3倍標準差作為數(shù)據(jù)上下限對其奇異值進行迭代替換，對剔值后的殘差值進行統(tǒng)計(表5)。

表5 剔值后殘差值特征統(tǒng)計

從表5可以看出，處理后的數(shù)據(jù)服從類似正態(tài)分布，可以用傳統(tǒng)方法圈定異常，異常下限取值時，選擇K=1.5;求出異常下限值后，用該值的1.0、2.0、4.0倍作為異常外帶、中帶、內(nèi)帶的下限值，并根據(jù)下限值圈定單元素的異常圖(圖4)。

結(jié)合圖2、圖3及圖4可以看出，各種方法作出的異常圖對已知礦化點都有一定程度上的吻合，但很明顯可以看出來滑動平均法對已知礦化點吻合度是最好的，很多重要的礦化點都落到了綜合異常區(qū)域，比如尖山幅西南角弱背景區(qū)內(nèi)重要的Cu礦化點D10295(25.99%)、D0172(1.28%)周圍圈到了異常，最終它們將圈進綜合異常區(qū)域。

與趨勢面法相比，滑動平均法以小窗口擬合局部異常，以大窗口擬合背景，小窗口及其近鄰的數(shù)據(jù)都參與背景運算，可以排除局部變化對背景趨勢的影響，使得背景更加“干凈”，不僅可以突出強異常及壓抑背景噪音的影響，而且還可以強化低背景區(qū)弱小異常，因而它能發(fā)現(xiàn)很多趨勢面法沒有發(fā)現(xiàn)的新異常。

滑動平均法作出的異常圖圖面結(jié)構(gòu)均衡，異常形態(tài)比較規(guī)整，濃度分帶明顯，異常套合性較好，比如石樁子井幅東南角重要的Cu礦化點K6(1.28%)和采樣點1377-A2(0.2%)部位提取到Cu、Ni、Co、Cr、Mn 異常，在已知礦化點和未知礦化部位均有強異常;傳統(tǒng)方法得出的異常除了在已知礦化點及高背景值的公婆泉組變質(zhì)火山巖中有大面積的異常外，其他地方基本是一些面積小的弱異常;趨勢面法圈出的異常類似傳統(tǒng)法圈出的異常。

另外，滑動平均法是建立在非常規(guī)統(tǒng)計學(xué)基礎(chǔ)上的數(shù)據(jù)處理方法，不受原始數(shù)據(jù)分布規(guī)律的制約，可直接利用原始數(shù)據(jù)來識別異常，不需要對原始數(shù)據(jù)進行預(yù)處理，計算不受特高值的影響，這樣就比其他方法操作方便，人工干擾因素較少。

3.4 異常提取效果的對比評價

通過圖2—圖4的對比，就提取異常的效果而言，滑動平均法明顯優(yōu)于其他兩種方法，主要體現(xiàn)在以下8個方面。

圖4 滑動平均法單元素異常圖

(1)滑動平均法對背景曲面擬合要優(yōu)于其他兩種方法。由圖5可知，傳統(tǒng)方法背景曲面是統(tǒng)一大小的水平面，圈出的異常主要在高背景區(qū)的地層中;多項式回歸法擬合的背景面較粗略，擬合的背景面是傾斜的平面，在一定程度上能壓制高背景區(qū)強異常及突出低背景區(qū)的弱緩異常;滑動平均法擬合的背景是起伏變化的背景曲面，對背景的擬合最好，它能有效弱化高背景區(qū)的異常及強化低背景區(qū)的弱異常，圈出較多低背景區(qū)的弱緩異常。因此，滑動平均法在實際應(yīng)用中效果較好。

圖5 三種方法背景曲面擬合示意圖

(2)對于礦化點的吻合，滑動平均法在3種方法中最優(yōu)。在已知礦化點及其附近部位，滑動平均法圈出的單元素異常面積都要比兩種方法的大，特別是對石樁子井幅東南角的礦化點吻合最好，6種元素都圈出異常且套合很好?；瑒悠骄ㄈΤ龅木C合異常與大多數(shù)礦化點吻合，且其綜合異常圖與重要的礦化點如 D10295、K3、K6、D4043吻合;而傳統(tǒng)法和趨勢面法圈出的異常與礦化點沒有吻合或者僅一兩種元素吻合，其綜合異常圖也未與礦化點吻合。

(3)滑動平均法作出的異常圖較為規(guī)整，異常以圓形、長條狀和橢圓狀為主，傳統(tǒng)法和趨勢面法作出的異常圖以尖棱狀和長條狀為主。

(4)滑動平均法提取異常的背景區(qū)較干凈，基本沒有形態(tài)不規(guī)整的小異常。

(5)滑動平均法能很有效地弱化強背景區(qū)的大面積異常和強化低背景區(qū)的弱緩異常，圈出的異常圖圖面結(jié)構(gòu)較均衡，每個圖幅都有濃集分帶明顯的異常，每個元素異常圖中異常的分布位置大體相同且異常的套合及濃度分級較好;傳統(tǒng)方法圈出的異常中，Cu、Mn沿深大斷裂呈長條狀分布，Ni、Cr在西林陶勒幅幾乎沒有異常，這明顯與其他元素的異常圖很不協(xié)調(diào);趨勢面法圈出的異常也有類似的現(xiàn)象。

(6)滑動平均法提取的每個元素的異常面積適中，約占全區(qū)面積的15%，形態(tài)比較規(guī)整，異常的分帶較好，三級濃度分帶的異常比較多;傳統(tǒng)方法提取的異常除Cu、Au、Mn外，其他元素的異常面積都偏小，約占全區(qū)面積的8%，形態(tài)和分帶性也不好;趨勢面提取的異常面積都偏小，形態(tài)不規(guī)則，分帶性也不好。

(7)滑動平均法能有效抑制地形地貌對異常形態(tài)的干擾。北北西向大斷裂貫穿全區(qū)，斷裂中心是一條沖積溝，兩邊高山的巖屑向下流動。因此，傳統(tǒng)方法和趨勢面法圈出的異常都有沿著斷裂帶呈帶狀分布的特點，顯示出地形地貌控制異常形態(tài)的特點;滑動平均法圈出的異常就能抵制地形地貌的干擾，異常以圓狀、橢圓狀和長條狀為主，形態(tài)比較規(guī)整，濃度分帶比較好。

(8)滑動平均法不需要對原始數(shù)據(jù)進行預(yù)處理，人為干擾因素最少，而且操作方法簡便。

因此，滑動平均法較為適合西北半干旱地區(qū)異常信息的提取。

4 結(jié)論

通過對比3種常用的異常提取方法，探討了較為適合西北半干旱地區(qū)的異常提取方法，主要取得以下幾點認識和成果。

(1)每一種異常提取方法都是基于對地質(zhì)數(shù)據(jù)分布規(guī)律的認識。傳統(tǒng)方法認為是地質(zhì)數(shù)據(jù)可以擬合成正態(tài)分布模式，關(guān)鍵是將原始數(shù)據(jù)通過數(shù)學(xué)變換成正態(tài)分布模式，再求異常下限值?；瑒悠骄ê挖厔菝娣ㄕJ為數(shù)據(jù)在空間上的分布呈連續(xù)變化，趨勢面就是用多元函數(shù)模擬背景曲面的變化，再用原始分布與其差值求異常的分布?；瑒悠骄ㄓ眯〈翱谀M局部或噪音的變化，用大窗口模擬背景曲面的變化，再用兩者間的差值求異常的分布。

(2)滑動平均法較適合西北半干旱地區(qū)異常的提取，它能有效地弱化分解高背景區(qū)的異常，同時還能強化低背景區(qū)低緩的弱異常，并且提取出來的異常圖圖面結(jié)構(gòu)均衡，異常形態(tài)規(guī)整，能夠抑制地形地貌對異常的干擾，濃集分帶明顯，異常套合性較好，與地質(zhì)體和礦化點的吻合也較好。

馮偉華，張秦華，鮮鵬輝.2007.三層套合方差分析在區(qū)域化探中的應(yīng)用及計算機實現(xiàn)——以新疆某地1∶50 000區(qū)域化探為例［J］.物探化探計算技術(shù)，29(1)：44-47.

甘肅省地質(zhì)局.1972.石板井幅1∶20萬區(qū)域地質(zhì)調(diào)查K-47-(21)［R］.

韓東昱，龔慶杰，向運川.2004.區(qū)域化探數(shù)據(jù)處理的幾種分形方法［J］.地質(zhì)通報，23(7)：714-719.

蔣敬業(yè).2006.應(yīng)用地球化學(xué)［M］.武漢：中國地質(zhì)大學(xué)出版社.

李慧生，高德秀，杜明亮，等.1988.用滑動平均法確定地球化學(xué)異常［J］.數(shù)理統(tǒng)計與管理，(3)：32-34.

李隨民，姚書振，韓玉丑.2007.Surfer軟件中利用趨勢面方法圈定化探異常［J］.地質(zhì)與勘探，43(2)：72-75.

內(nèi)蒙古自治區(qū)地質(zhì)礦產(chǎn)局.1996.內(nèi)蒙古自治區(qū)巖石地層［M］.武漢：中國地質(zhì)大學(xué)出版社.

聶鳳軍，江思宏，白大明，等.2003.蒙甘新相鄰(北山)地區(qū)金銅礦床時空分布特征及成礦作用［J］.礦床地質(zhì)，22(3)：234-245.

潘桂棠，肖慶輝，陸松年，等.2009.中國大地構(gòu)造單元劃分［J］.中國地質(zhì)，26(1)：1-4.

謝學(xué)錦，邵躍，王學(xué)求.1999.走向21世紀礦產(chǎn)勘查地球化學(xué)［M］.北京：地質(zhì)出版社.

趙鵬大，胡旺亮，李紫金.1993.礦床統(tǒng)計預(yù)測［M］.北京：地質(zhì)出版社.

趙榮軍.2006.不同方法在欒川北部化探數(shù)據(jù)處理中的應(yīng)用［J］.地質(zhì)與勘探，42(3)：67-71.

張雨蓮，許榮科，陜亮，等.2009.北山公婆泉斑巖銅礦區(qū)域地質(zhì)地球化學(xué)異常結(jié)構(gòu)模式［J］.物探與化探，33(4)：406-410.

DZ/T 0011—1991.地球化學(xué)普查規(guī)范(比例尺1∶50 000)［S］.

Study on extraction method of rock geochemical anomaly in semi arid region in Northwest China

FAN Xiao-jun1，2，CHEN Chong2，WANG Xiao-gang1，CHEN Gen1

(1.Sichuan Institute of Nuclear Geology，Chengdu 610066，China;2.School of Resources，China University of Geology，Wuhan 430074，China)

Residual hills landforms were developed in Northwest China.The widespread development of loess crust and its underlying rock layers formed by aeolian sand and rain affected the representative of geochemical investigation data，brought difficulties for the extraction and explanation of the anomaly.In lights of the above unique natural landscape and geological settings，the authors presented the geological properties of a place in Alashan of Northwest China as a study case，extracted the related geochemical anomalies with conventional methods of statistical method，trend surface analysis and moving average method.And a comprehensive contrast was conducted in combination with the available information on geology and mineral resources，proper method of geochemical anomaly extraction in semi arid region of Northwest China was discussed in the text.

Rock geochemical anomaly survey;Anomaly extraction;Trend surface analysis;Moving average method;Semi arid region in Northwest China

P618

A

1674-3636(2012)01-0023-10

2011-09-29;編輯：陸李萍

范小軍(1986— )，男，助理工程師，資源勘查工程專業(yè)，主要從事地質(zhì)化探找礦工作，E-mail：364102412@qq.com

10.3969/j.issn.1674-3636.2012.01.23