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

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

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

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

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

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

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

0 引言

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

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

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

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

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

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

2.1 地層

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

2.2 構(gòu)造

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

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

2.3 巖漿巖

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

3 異常提取

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

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

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

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

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

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

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

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

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

研究區(qū)內(nèi)風(fēng)化剝蝕較為嚴(yán)重，巖屑中各元素含量有所變化，并且項(xiàng)目屬于1∶5萬(wàn)大面積化探，為了不損失一切可能的礦化信息，Co、Cr、Mn異常下限取值應(yīng)相對(duì)降低，信度系數(shù)應(yīng)提高，選擇信度α =0.1對(duì)應(yīng)的 K=1.65;Cu、Au、Ni是取對(duì)數(shù)后服從正態(tài)分布的，取對(duì)數(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明顯可見(jiàn)，異常主要與變質(zhì)火山熔巖和中基性巖關(guān)系最為密切，主要分布在公婆泉組(Sg)地層、輝長(zhǎng)巖及圓包山組(Sy)地層，且與地質(zhì)體形態(tài)吻合較好;Cu、Au、Ni、Co、Cr、Mn 在尖山幅輝長(zhǎng)巖區(qū)附近面積比較大，異常有明顯的濃度分帶;Cu異常主要分布在公婆泉組(Sg)地層中，與變質(zhì)火山熔巖的關(guān)系最為密切;Au異常除了主要分布在公婆泉組(Sg)地層中外，其次分布在黑云母花崗巖、赤金堡組(Kc)及綠條山組(Cl)地層中，與變質(zhì)火山熔巖、黑云母花崗巖、砂巖關(guān)系密切;提取的異常與Cu、Au礦化點(diǎn)有一定的吻合度，主要是與尖山輝長(zhǎng)巖附近礦化點(diǎn)吻合。

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

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

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

3.2 基于Surfer軟件的趨勢(shì)面法

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

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

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

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

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

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

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

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

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

3.3 滑動(dòng)平均法

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

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

圖3 趨勢(shì)面法單元素異常圖

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

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

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

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

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

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

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

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

3.4 異常提取效果的對(duì)比評(píng)價(jià)

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

圖4 滑動(dòng)平均法單元素異常圖

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

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

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

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

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

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

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

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

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

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

4 結(jié)論

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

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

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

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

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

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

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

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

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

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

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

潘桂棠，肖慶輝，陸松年，等.2009.中國(guó)大地構(gòu)造單元?jiǎng)澐郑跩］.中國(guó)地質(zhì)，26(1)：1-4.

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

趙鵬大，胡旺亮，李紫金.1993.礦床統(tǒng)計(jì)預(yù)測(cè)［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