歐陽立勝

（廣東省地震局，廣東 廣州 510070）

分形濾波(S-A)技術(shù)在物探數(shù)據(jù)異常分解中的應(yīng)用

歐陽立勝

（廣東省地震局，廣東 廣州 510070）

分形濾波 (S-A)技術(shù)能夠有效地揭示地球物理場(chǎng)的空間相關(guān)性，對(duì)背景場(chǎng)以及異常場(chǎng)進(jìn)行分解。應(yīng)用分形濾波 (S-A)技術(shù)分解內(nèi)蒙古某銀礦區(qū)地面高精度磁測(cè)異常，結(jié)果表明：對(duì)研究區(qū)高精度磁異常進(jìn)行了有效分解，獲得了研究區(qū)背景場(chǎng)以及異常場(chǎng)特征，并剖析所獲得的異常場(chǎng)，分解異常8處，并據(jù)此結(jié)合地質(zhì)情況提出研究區(qū)的成礦有利區(qū)。

分形濾波；地面高精度磁測(cè)；異常分解

引言

在礦產(chǎn)勘查中，我們需要通過實(shí)施各種勘查技術(shù)來測(cè)量獲取所需的數(shù)據(jù),特別是通過物化探的方法獲取各種直接或者間接的成礦空間分布信息。然而我們所獲取的地球物理及地球化學(xué)數(shù)據(jù)往往是反映地表及地下多種地質(zhì)體和地質(zhì)過程的疊加結(jié)果，因此通過對(duì)地球物理及地球化學(xué)場(chǎng)的分解來研究有關(guān)的地質(zhì)體和地質(zhì)過程至關(guān)重要。同時(shí)，由于地質(zhì)過程往往具有多期次和空間相關(guān)性等特點(diǎn),所產(chǎn)生的地質(zhì)體或者相應(yīng)的場(chǎng)往往會(huì)呈現(xiàn)空間自相似性或統(tǒng)計(jì)自相似性。 這種自相似性可以幫助分解不同地質(zhì)過程所產(chǎn)生的同類地質(zhì)體［1］?；谶@一特征，文中選取內(nèi)蒙古西烏旗某銀礦床為研究對(duì)象，應(yīng)用分形濾波方法 （S-A）來分解研究區(qū)內(nèi)的物化探局部異常和背景場(chǎng)。研究表明,該方法對(duì)區(qū)分與礦有關(guān)的局部異常與背景場(chǎng)是很有效的。

1 S-A濾波方法原理

許多地質(zhì)過程包括礦化過程均可造成具有自相似性或統(tǒng)計(jì)自相似性的場(chǎng)-多重分形場(chǎng)。這種場(chǎng)的特點(diǎn)是場(chǎng)的形態(tài)變化具有與度量尺度的獨(dú)立性[1]。多重分形分析方法是一種用于研究具有自相似性或統(tǒng)計(jì)自相似性分布規(guī)律的有效方法，如果將多重分形場(chǎng) （測(cè)度）表示為f（ε），它表示在尺度為ε×ε范圍內(nèi)的某種測(cè)度。比如，一定范圍內(nèi)的地球化學(xué)元素含量;一定范圍內(nèi)的地球物理場(chǎng)的測(cè)度。這些測(cè)度與尺度呈指數(shù)關(guān)系:

α稱為奇異性指數(shù)，由不同α指數(shù)對(duì)應(yīng)的空間集合構(gòu)成多個(gè)空間鑲嵌的分形，稱為多重分形，其分形維數(shù)用維數(shù)譜函數(shù)來描述[2]。由于地質(zhì)作用的差異所形成的地質(zhì)體，其空間模式具有不同的廣義自相似性。它們?cè)陬l率域中也具有不同的廣義自相似性。在頻率域中對(duì)這些不同的廣義自相似性進(jìn)行識(shí)別，并采用濾波的方法選取部分頻率信息重新恢復(fù)空間模式，就可達(dá)到對(duì)異常場(chǎng)和背景場(chǎng)分離的目的。能譜密度S與能譜密度大于S的波數(shù)集合的面積A（波數(shù)的平方）之間的關(guān)系可以用下式來描述

用最小二乘法對(duì)logA～logS關(guān)系分段擬合求得對(duì)應(yīng)于不同的能譜密度范圍的冪指數(shù)β，根據(jù)不同的能譜密度范圍設(shè)計(jì)濾波器濾波，達(dá)到對(duì)場(chǎng)分離的目的[2～4]。目前該方法在Morpas3.0軟件上已經(jīng)得到實(shí)現(xiàn)。

2 研究區(qū)地質(zhì)概況

內(nèi)蒙古西烏旗某銀礦床屬中亞－蒙古地槽褶皺區(qū)、內(nèi)蒙古中部地槽褶皺系之內(nèi)蒙古大興安嶺地槽褶皺帶晚期的華力西地槽褶皺帶[5]。礦區(qū)出露地層為古生界下二疊統(tǒng)壽山溝組（P1s）砂巖、細(xì)砂巖、粉砂巖；中生界侏羅系上統(tǒng)瑪尼吐組 （J3mn）安山巖、安山玢巖和滿克頭鄂博組 （J3m）酸性凝灰?guī)r、集塊巖、含角礫凝灰?guī)r及沉凝灰?guī)r；第三系上新統(tǒng)五叉溝組 （N2Wc）黑色玄武巖和寶格達(dá)烏拉組 （N2b）棕色粉砂質(zhì)泥巖、紅色粘土和灰色砂礫巖；第四系 （Q4）沖洪積、沖坡積物及殘坡積碎石、風(fēng)成砂及亞砂土。

礦區(qū)斷裂發(fā)育，在主斷裂帶及兩側(cè)發(fā)育一系列的北西向、北東向及近南北向斷裂。在礦區(qū)內(nèi)以二采區(qū)為中心，形成一系列環(huán)狀、放射狀斷裂，這些斷裂為礦液的運(yùn)移和賦存提供了空間，斷裂構(gòu)造是本區(qū)的主要控礦因素。

礦區(qū)內(nèi)的巖漿巖主要為華力西晚期的超基性巖體和燕山期的巖株、巖墻和巖脈。礦區(qū)內(nèi)所見到的主要巖漿巖有：華力西晚期超基性巖，燕山晚期中性、酸性侵入體，花崗斑巖以及次流紋巖脈，花崗斑巖與閃長(zhǎng)玢巖脈等脈巖。

高精度磁測(cè)范圍為礦區(qū)西側(cè)，面積2.42 km2，區(qū)內(nèi)第四系覆蓋面積約占80%以上，地表可見有褐或褐紅色Fe－硅質(zhì)蜂窩狀石英脈出露于研究區(qū)的東側(cè)，地形為東高西低，實(shí)測(cè)研究區(qū)海拔高度為941.49～1 014.06 m，相對(duì)高差73 m。

3 基礎(chǔ)數(shù)據(jù)

本文選取物探數(shù)據(jù)為高精度磁測(cè)數(shù)據(jù) （ΔT），研究區(qū)測(cè)網(wǎng)布置為正南北向，線距50 m，點(diǎn)距10 m，總測(cè)點(diǎn)數(shù)5 292個(gè)，經(jīng)過一定的預(yù)處理以及校正，繪制ΔT等值線 (圖2)。

總磁場(chǎng)強(qiáng)度 （ΔT）極大值為392.42 nT，極小值－561.18 nT，從圖2中看，主要從異常區(qū)包括：圖幅北西，具有北東走向的長(zhǎng)條狀高值磁異常區(qū) （箭頭所示），其中包含幾個(gè)單峰異常，從磁異常特征分析推測(cè)為隱伏不深的花崗巖或花崗閃長(zhǎng)巖巖體引起；圖幅 （圖2）東南，包含強(qiáng)磁場(chǎng)異常，呈不規(guī)則狀半環(huán)形分布 （箭頭所示），環(huán)內(nèi)為負(fù)值異常場(chǎng)，依據(jù)礦區(qū)出露巖性層及研究區(qū)東部鉆孔揭示的巖層對(duì)比分析，推測(cè)為不規(guī)則狀蛇紋巖與砂巖接觸帶異常，上部負(fù)異常區(qū)輪廓清晰，推測(cè)為礦區(qū)含礦斷裂帶的西延伸部分，區(qū)內(nèi)地形較高地段可見氧化鐵帽分布于負(fù)磁異常區(qū)，為有利找礦勘查區(qū)。從上延2 km磁異常特征可見 （圖3），其具有向西南弧形轉(zhuǎn)折的趨勢(shì)，構(gòu)成外圍及深部找礦遠(yuǎn)景區(qū)。

圖1 研究區(qū)ΔT化極磁異常Fig.1 The magnetic pole anomaly（ΔT）in study area

圖2 研究區(qū)ΔT上延2 km趨勢(shì)圖Fig.2 The upward trend map in study area

4 S－A濾波方法異常提取

4.1 異常提取

在Morpas3.0系統(tǒng)下，應(yīng)用S－A分形濾波方法對(duì)高精度磁測(cè)數(shù)據(jù)進(jìn)行進(jìn)一步的異常提取，以盡可能的消除背景的差異對(duì)異常的影響。從圖3我們獲取了lnA（＞S）－lnS即面積－能譜密度雙對(duì)數(shù)圖關(guān)系圖 （圖4），橫坐標(biāo)代表能譜密度對(duì)數(shù)值變換，縱坐標(biāo)表示面積的對(duì)數(shù)值。該方法確保了S和A之間的冪律關(guān)系，不同斜率的直線通常對(duì)應(yīng)空間域中不同的地球物理背景和異常場(chǎng)分布模式，斜線段由最小二乘方法擬合，豎線表示不同指數(shù)關(guān)系的分界。 左邊斜線段反映低能譜分布范圍并形成異常濾波器，右邊斜線段反映高能譜分布范圍并形成背景濾波器。圖中，用最小二乘法模擬的右邊兩條線段的交點(diǎn)確立的閾值lnS0= 10.05。 S＜S0通常代表異常，S＞S0代表背景。

圖3 高精度磁測(cè)ΔT的lnA(＞S)－lnS圖Fig.3 The lnA(＞S)－lnS map of high－precision magnetic survey data ΔT

圖4 研究區(qū)ΔT由 S－A濾波獲取背景場(chǎng)Fig.4 The background field obtained by ΔT S－A filtering in study area

圖5 研究區(qū)ΔT由 S－A濾波獲取異常場(chǎng)Fig.5 The anomaly field obtained by ΔT S－A filtering in study area

通過S－A方法，并以閾值lnS0=10.05獲取高精度磁測(cè)空間域中的背景場(chǎng)和異常場(chǎng)分別如圖5和6，從圖5上看，研究區(qū)的磁測(cè)背景表現(xiàn)為以北東為軸向的一系列相間的正異常以及負(fù)異常，背景的最高值出現(xiàn)在南東，通過地質(zhì)檢測(cè)，該區(qū)域廣泛發(fā)育以蛇紋巖為主的超基性蝕變巖，因此造成較高的背景值。

4.2 異常評(píng)價(jià)

相較于背景場(chǎng)，對(duì)于異常場(chǎng)往往更為重視，因?yàn)楫惓?chǎng)是我們進(jìn)行礦產(chǎn)勘查以及遠(yuǎn)景區(qū)預(yù)測(cè)均的重要依據(jù)。通過S-A濾波方法分解異常場(chǎng)如圖6，圖6與圖2在某些范圍存在一定的相似，但也有區(qū)別，將其中的弱異常和隱蔽異常提取出來，并將各個(gè)異常編號(hào)，結(jié)合地質(zhì)情況，并進(jìn)行描述及評(píng)價(jià)。

Ⅰ號(hào)異常：位于研究區(qū)北西，由一串南北向的面積較小的單峰異常組合而成，對(duì)應(yīng)巖性為滿酸性流紋巖以及流紋質(zhì)火山角礫巖；

Ⅱ號(hào)異常：位于研究區(qū)北，由軸向?yàn)楸睎|的面積較大的雙峰異常組成，對(duì)應(yīng)的巖性推測(cè)為隱伏花崗巖，該異常有隱伏斷裂通過；

Ⅲ號(hào)異常：位于研究區(qū)北東，范圍較小，巖性為流紋巖與推測(cè)隱伏花崗巖接觸帶；

Ⅳ號(hào)異常：位于研究區(qū)中部，呈長(zhǎng)條形，對(duì)應(yīng)巖性為酸性流紋巖以及流紋質(zhì)火山角礫巖，靠近蛇紋巖東部,斷裂通過此區(qū)域；

Ⅴ號(hào)異常：由半環(huán)形異常組成，主要包含三個(gè)單峰異常區(qū)，該區(qū)域?qū)?yīng)巖性為華力西晚期蛇紋巖，該異常值較大；

Ⅵ 號(hào)異常：位于研究區(qū)東，地質(zhì)特征與Ⅴ號(hào)異常相似，兩者包圍一較為明顯的負(fù)異常區(qū)域；

Ⅶ號(hào)異常：帶狀，在巖性上恰好位于壽山溝組板巖砂巖與酸性凝灰?guī)r、集塊巖接觸帶；

Ⅷ號(hào)異常：位于研究區(qū)西部，孤立，面積較小，推測(cè)巖性為隱伏酸性巖。

結(jié)合本礦區(qū)成礦特點(diǎn)特別是超基性蛇紋巖與成礦的關(guān)系，認(rèn)為其中的Ⅳ、Ⅴ以及Ⅵ號(hào)異常區(qū)為研究區(qū)根據(jù)高精度磁測(cè)所推測(cè)的有利的成礦遠(yuǎn)景區(qū)。

5 結(jié)論

通過分形濾波S-A方法分析地球物理異常的富集特征，有利于識(shí)別與礦化相關(guān)的地球物理異常，有助于分析地球物理數(shù)據(jù)的空間自相似性，對(duì)區(qū)分不同的異常，獲取與礦化有關(guān)的信息具有獨(dú)到之處，將是一種識(shí)別非常規(guī)、難識(shí)別的礦藏的有效方法。

［1］成秋明.空間自相似性與地球物理和地球化學(xué)場(chǎng)的分解方法 ［J］.地球物理學(xué)進(jìn)展,2001,16（2）：8-17.

［2］成秋明.多重分形與地質(zhì)統(tǒng)計(jì)學(xué)方法用于勘查地球化學(xué)異?？臻g結(jié)構(gòu)和奇異性分析 ［J］.地球科學(xué)-中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2001,26（2）：161-166.

［3］李慶謀.分形奇異 (特征值)分解方法與地球物理和地球化學(xué)異常重建 ［J］.地球科學(xué)-中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2004,29（1）：110-118.

［4］李慶謀,劉少華.GIS環(huán)境下地球物理信號(hào)的奇異值分解、多維分形特征與應(yīng)用 ［J］.地球物理學(xué)進(jìn)展, 2003,18（1）：97-102.

［5］劉建明,張銳,張慶洲.大興安嶺地區(qū)的區(qū)域成礦特征 ［J］.地學(xué)前緣,2004,11（1）：269-277.

Abstract：Fractal filtering(S-A method)can reveal the spatial correlation of geophysical field and decompose the background field and anomaly field.Using fractal filtering(S-A method)to extract the geophysical anomalies of a silver deposits in Inner Mongolia,the results show that fractal filtering technology can decompose the high-precision magnetic anomaly effectively.Based on the characteristic of the background field and anomaly field in study area and analysis on the anomaly field,eight anomaly zones are extracted.Combining with the geological features,we give some suggestions on ore advantageous prospecting region in the study area.

Keywords：Fractal filtering(S-A method);Ground high-precision magnetic survey;Unusual decomposition

Application of Fractal Filtering in Unusual Decomposition of Geophysical Data

OUYANG Lisheng
(Earthquake Administration of Guangdong Province，Guangzhou 510070,China)

TU435

A

1001-8662（2010）04-0070-06

2010-07-20

歐陽立勝，男，1974年生，工程師.主要從事工程抗震研究. E-mail:1095705480@qq.com.