匹配濾波方法在場源分離上的應用

張 翔，常樹帥，李兵海

（核工業航測遙感中心，石家莊 050002）

匹配濾波方法在場源分離上的應用

張 翔，常樹帥，李兵海

（核工業航測遙感中心，石家莊 050002）

介紹了匹配濾波方法的原理、應用條件、實現過程。用該方法對內蒙古白云鄂博礦區的主礦段和西礦段航磁異常進行了垂向位場的分離，利用徑向功率譜法對深、淺源場異常進行磁性體埋深反演，反演結果表明，該區在1 000m～1 300m之間可能還存在鐵礦體。驗證了匹配濾波方法在一定條件下，可以較為快速地進行位場分離，實現磁性體深度準確反演。

匹配濾波；功率譜法；白云鄂博礦區；場源分離；航磁異常

0 前言

匹配濾波方法是分離區域場和局部場的關鍵技術之一，對準確反演磁性體深度及異常解釋有著重要意義。疊加的地質體其磁異常的頻譜具有不同特點，①淺部地質體產生的異常比深部地質體產生的異常要強的多，具有高頻的特點；②深部磁性體產生的異常寬緩，具有低頻特點。這種差異使得從位場圖中濾掉近地表的高頻干擾得到深部的寬緩異常成為可能；反之也可以濾去深部低頻部分，達到突出淺部高頻異常的目的［1］。Spector［2］曾使用匹配濾波方法處理科迪雷拉山區的航磁數據，有效地消除了火山巖覆蓋的干擾，從而得到與成礦有關的深部火成巖引起的異常。賀日政等［3］利用匹配濾波方法重新解釋了青藏高原中西部航磁ΔT化極異常數據，推斷青藏高原內的中、上地殼的分界深度為19km。該深度青藏高原中、上地殼間殼內脆韌性轉換帶，此外劉青松等［4］利用兩次匹配濾波方法，對沙市實測布格重力異常進行處理，比較合理地分離出由基底起伏鹽系地層引起的異常。

1 方法原理

重磁異常經從空間域到頻率域的轉化后，為全頻段信號，它的功率譜和振幅譜為指數形式，因此它的徑向平均對數功率譜或振幅譜，可被不同斜率的直線所擬合，直線的斜率與其等效層的深度相關，在這種以直線擬合位場對數譜（功率譜或振幅譜）的模型下，匹配濾波技術可用來分離垂向疊加場［5］。

設觀測場頻譜由深部和淺部地質體異常疊加之和即

式中：T深、T淺分別為深部和淺部磁性體引起的磁場值，與磁化強度與磁性體的水品尺度有關，而與徑向頻率r無關；H和t1為深部磁性體的埋深和厚度，相應的h和t2為淺部磁性體的埋深和厚度。對于深部磁性體不妨認為t1→∞，有A深≈T深e－Hr，那么有

公式（4）即為深、淺部磁性體疊加磁場的總能譜表達式。

當r較大（高頻段）時，有

W1即為“濾波算子”，當要濾去淺部高頻成分時，只需將ΔT（u，v）乘以W1就可達到目的，同理，也可將E（r）寫成式（6）的形式：

這時當要濾去深部低頻成分時，將ΔT（u，v）乘以W2構成“匹配濾波”即可。

因此構建匹配濾波器的關鍵在于求出W1和W2兩個“濾波算子”，而W1和W2依賴于H、h、T深、T淺等參數。

由式（7）與式（8）采用圖解法（直線擬合法）確定所需的H、h、T深、T淺等參數。徑向頻率對數譜可被不同斜率的直線擬合，直線的斜率即為相應的等效層的深度。

上述過程可以一次在計算機上實現：對于算出的功率譜，利用可視化技術顯示對數功率譜曲線，并用鼠標畫出深源場與淺源場的回歸直線，計算其斜率和縱軸截距，即可構制出匹配濾波算子，進行垂向不同深度場源的分離。具體實現步驟如下：

1）利用傅立葉變換，由實測異常求頻譜。

2）由傅立葉變換的實部與虛部求對數功率譜。

3）根據對數功率譜曲線LnE（r）—r，求H、h、T深、T淺等參數。

4）把實測異常頻譜乘以相應濾波因子，得到淺源場或深源場的頻譜。

5）反傅立葉變換得到分離的淺源場與深源場。

2 應用實例

2．1 地質背景及航磁特征

研究區為著名的內蒙古白云鄂博稀土、鈮、鐵礦區的主礦和西礦段（圖1）。研究區地處內蒙古地軸北緣向內蒙古古生代地槽的過渡地帶，區內主要分布有中元古代的下白云鄂博群地層，為淺變質碎屑巖夾碳酸鹽巖及鐵礦，該層在研究區內分四個巖組，①哈拉霍疙特組（Jxh）；②比魯特組（Jxb）；③尖山組（Chj）；④都拉哈拉組（Chd）。其中尖山組（Chj）地層和都拉哈拉組（Chd）組地層為區內的含礦層位，巖性以板巖、石英巖、石英夾大理巖為主，局部見含鐵石英巖及灰巖。區內的侵入巖以二疊紀、石炭紀花崗巖為主，分布于礦床南北，其次是輝長巖類、閃長巖類等。研究區內發育多條斷裂，集中位于研究區東北部，走向主要為東西向、北西西向，地層產狀受這些斷裂構制的控制［6］。

由圖1可見，研究區中東部有一條東西向展布的強磁異常帶，長約14km，寬約5km，異常正負伴生，南正北負。異常幅值一般在5 000nT以上，最大值大于10 000nT。受地磁場斜磁化的影響，異常北側伴生有－2 600nT的負異常。該異常帶整體與白云鄂博稀土、鈮、鐵礦區主礦段和西礦段相對應。異常周圍為平靜的負背景場，場值在0nT～－200nT之間。

2．2 數據處理

分析的數據來源于核工業航測遙感中心承擔的“內蒙古包頭市－巴彥淖爾市一帶1∶50 000航磁、航放測量項目”［7］。數據處理流程見圖2。

圖1 研究區航磁ΔT異常圖及地質圖Fig．1 The map of the study area'sΔTaeromagnetic anomalies and geology

圖2 數據處理流程圖Fig．2 The flowchart of data processing

航磁ΔT值是地磁場相對變化量，是磁場總強度（航磁儀實測值）與正常地磁場總強度的模差。它反映的是在一定深度范圍內，地殼中沿水平和垂向分布的各種具磁性的地質體形成的磁異常的疊加場。

由于地磁場傾斜磁化的影響，可能造成正磁異常中心不是正好對應在磁性體地質體的正上方，而是沿傾斜磁化強度矢量水平投影的反方向上有不同程度的偏移錯動［8］，為消除傾斜磁化的這種影響，做了化極處理（圖3）。本區化極采用參數：磁傾角68°、磁偏角為－4．93°。

對比圖3（a）和圖3（b）發現，異常北側伴生的強負磁異常在化極后消失，異常中心位置向北移動200m～500m不等，異常南北兩側伴生弱負磁異常，正異常圈閉范圍與礦體邊界對應較好，表明經化極后較為有效地消除了地磁場斜磁化的影響。

對化極后的數據進行了匹配濾波處理，最后進行磁性體場源深度反演。

2．3 結果分析

對化極后的△T磁場計算了徑向平均對數功率譜并繪制了LnE（r）—r曲線圖（圖4）。從圖4中可以看出，曲線明顯地被分為5段，利用線性回歸分析方法對不同段曲線進行擬合（表1）。

由表1發現，頻率在0（1／m）～0．05× 10－2（1／m）（圖4中頻率小于f1）的低頻段所對應擬合直線段AB的斜率KAB＝－2．9，頻率在0．05× 10－2（1／m）～0．14×10－2（1／m）內（圖4中頻率在f1－f2之間）的中高頻段所對應擬合直線段BC的斜率KBC＝－0．57，二者相差5．1倍。由此表明，利用匹配濾波方法可以較好地分離不同深度的場源體。此外受白云鄂博礦區采礦和地表人文干擾的影響，在構造匹配濾波因子進行濾波前，應先濾除頻率范圍在0．14×10－2（1／m）～0．29×10－2（1／m）（圖4中頻率在f2－f3之間）高頻區段的頻譜響應。

圖3 航磁ΔT平面等值圖和化極平面等值圖Fig．3 The contour map ofΔTaeromagnetic anomalies and the pole contour map of ΔTaeromagnetic anomalies

表1 線性回歸分析擬合功率譜曲線結果表Tab．1 The results of fitting power spectral curve by linear regression analysis

分離后的深、淺源場的磁場如圖5所示，由圖5（a）可以看出，深源場磁場面貌簡單，幅值大幅降低，反映了基底的磁場特點，研究區中部有似“啞鈴狀”異常，推測為基底凸起區域；淺源場磁場圖中（圖5（b）），異常幅值變化不大，但在強磁異常南北兩側伴生－2 000nT～－1 000nT的強負磁異常，表明磁鐵礦體向下延伸有限。對異常進行反演計算，深源場場源平均埋深為3 000m，淺源場場源平均埋深為1 300m，目前白云鄂博主礦和西礦段探明的鐵礦體最大埋深分別為1 030m和855m［9］，推測在1 000m～1 300m仍存在鐵礦體，值得深入研究。

3 結論

采用匹配濾波方法對白云鄂博西礦段和東礦段進行垂向位場分離，取得了較好效果，對分離的淺源場，反演計算場源體的平均埋深，并與目前該礦區鐵礦最大埋深進行對比，推測在1 000m～1 300m還可能有鐵礦體。

應用匹配濾波方法的物理前提是位場信號各頻段成分不相關，這可以從異常的徑向對數功率譜來判斷，不同頻率段擬合直線段斜率相差越大，濾波效果越好。

圖4 研究區航磁異常徑向平均對數振幅譜曲線圖Fig．4 The map of the study area's aeromagnetic anomalies average radial logarithmic amplitude spectrum curve

［1］曹禮剛，王文濤．基于匹配濾波算法的異常信息提取研究［J］．科技資訊，2011（17）：25－26．

CAO L G，WANG W T．The research of anomalies extraction base on the match filtering method［J］．Science and Technology News，2011（17）：25－26．（In Chinese）

［2］劉天佑，張勝業，潘玉玲．應用地球物理數據采集與處理［M］．北京：中國地質大學出版社，2006．

LIU T Y，ZHANG SH Y，PAN Y L，The date colletion and processing of applied geophysics［M］．Beijing：China university of geosciences publishing house，2006．（In Chinese）

圖5 經匹配濾波處理后的磁場圖Fig．5 The map of magnetic field after matched filtering process

［3］賀日政，高銳，鄭宏偉．青藏高原中西部航磁異常的匹配濾波分析與構造意義［J］．地球物理學報，2007，50（4）：1131－1140．

HE R ZH，GAO R，ZHENG H W．The match filtering analysis and structure significance for aeromagnetic anomalies of the Midwest of Qingzang Plateau［J］．Chinese Journal of Geophysice，2007，50（4）：1131－1140．（In Chinese）

［4］劉青松，王寶仁．應用多次匹配濾波技術進行垂向位場分離［J］．物探化探計算技術，1996，18（4）：279－286．

LIU Q S，WANG B R．The application of repeated match filtering method in vertical source separation［J］．Computing Techniques for Geophysical and Geochemical Exploration，1996，18（4）：279－286．（In Chinese）

［5］張培琴．位場異常頻譜與功率譜的統計分析［J］．物探化探計算技術，1991，13（1）：7－12．

ZHANG P Q．The Statistical analysis in spectrum and power spectrum of anomalies［J］．Computing Techniques for Geophysical and Geochemical Exploration，1991，1（1）：7－12．（In Chinese）

［6］內蒙古自治區地質礦產局，內蒙古自治區區域地質志［M］．北京：地質出版社，1991．

BUREAU OF GEOLOGY AND MINERAL RESOURCES OF INNER MONGOLIA．Regional geology of inner Mongolia［M］．Beijing：Geological publishing house，1991．（In Chinese）

［7］楊玉勤，杜化宇，李曉祿．內蒙古自治區巴彥淖爾市－包頭市一帶1：5萬航空磁法、伽瑪能譜測量報告［D］．核工業航測遙感中心，2012．

YANG Y Q，DU H Y，LI X L．The report of inner Mongolia Bayannaoer－Baotou area 1：50，000 aeromagnetic and gamma ray spectrometry survey［R］．Airbone Survey and Remote Sensing Center of Nuclear Industing，2012．（In Chinese）

［8］范正國，黃旭釗，熊盛青．磁測資料應用技術要求［M］．北京：地質出版社，2010．

FAN ZH G，HUANG X ZH，XIONG SH Q．The technology requirement for appling magnetic date［M］．Beijing：Geological publishing house，2010．（In Chinese）

［9］劉敬國，柳建勇．論白云鄂博鐵礦資源潛力與今后地質工作方向［J］．包鋼科技，2007，33（4）：9－14．

LIU J G，LIU J Y．The iron mine potential of Baiyunebo and outlook on the geological work［J］．Science and Technology of Baotou Steel Corporation，2007，33（4）：9－14．（In Chinese）

The application of match filtering method in source separation

ZHANG Xiang，CHANG Shu-shuai，LI Bing-hai
（Airborne Survey and Remote Sensing Center of Nuclear Industry，Shijiazhuang 050002，China）

The principle，application condition，implementation process of match filtering method was introduced in this paper，which was used in separating the vertical field of aeromagnetic anomalies in Inner Mongolia Baiyunebo mine＇s main ore section and west ore section．The radial logarithmic power spectrum method was used in magnetic anomaly depth inversion of deep－shallow source．The inversion results show that the area in 1 000to 1 300meters and possibly iron ore body．Under certain conditions，the matched filter can separate the magnetic field much quickly，and achieve the accurate depth of magnetic body inversion．

matched filter；power spectrum；Baiyunebo area；field source separation；aeromagnetic anomaly

P 631．2

A

10．3969／j．issn．1001-1749．2015．03．10

1001-1749（2015）03-0325-05

2014-08-15 改回日期：2014-11-12

內蒙古自治區項目（08－2－HK01）

張翔（1986－），男，工程師，主要從事航空物探數據處理和地球物理資料解釋工作，E-mail：zhangxiang19860626＠163．com。