鄭 偉 許 石 孫海龍 范玉海

（西安煤航遙感信息有限公司,陜西 西安 710199）

航空磁法測(cè)量具有效率高、周期短、不受地表?xiàng)l件約束、不受地形影響、可在不同高度測(cè)量及能大面積作業(yè)的優(yōu)點(diǎn),這些優(yōu)點(diǎn)使航空磁法廣泛應(yīng)用于地質(zhì)區(qū)域調(diào)查、金屬探礦、尋找油氣藏、推斷地球深部構(gòu)造等領(lǐng)域。航空磁法測(cè)量時(shí)按照預(yù)設(shè)的測(cè)線飛行,其測(cè)量方式?jīng)Q定了相鄰測(cè)線間存在著視水平漂移,需要進(jìn)行調(diào)平處理,調(diào)平的質(zhì)量直接決定了航磁資料的質(zhì)量,對(duì)數(shù)據(jù)處理和解釋具有重要的意義[1-3]。

1 調(diào)平方法研究

調(diào)平誤差可分為系統(tǒng)誤差和系統(tǒng)噪聲兩大類(lèi)。系統(tǒng)誤差源包括地球日變磁場(chǎng)的變化、飛行器和探測(cè)器之間的滯后誤差、探頭的方向差等,系統(tǒng)誤差可以通過(guò)簡(jiǎn)單的加減算法校正；系統(tǒng)噪聲是校正系統(tǒng)誤差以后所保留的誤差,它通常是由數(shù)據(jù)中的梯度和高度誤差引起的,該誤差通過(guò)調(diào)平方法進(jìn)行校正[4-6]。筆者利用Geosoft 公司Oasis montaj 軟件附帶的航磁資料為例對(duì)航磁調(diào)平的方法進(jìn)行了研究（原始未處理磁場(chǎng)圖如圖1）,后應(yīng)用于山陽(yáng)地區(qū)航磁資料處理中,為地質(zhì)找礦提供了磁法依據(jù)。

圖1 實(shí)驗(yàn)數(shù)據(jù)原始未處理磁場(chǎng)圖

1.1 預(yù)處理

航磁調(diào)平的預(yù)處理包括日變磁場(chǎng)改正、正常場(chǎng)改正及飛行海拔高度校正。航磁預(yù)處理與地磁處理方法基本相同,分別根據(jù)日變觀測(cè)得到的日變數(shù)據(jù)、國(guó)際地磁參考模型及測(cè)量飛行出現(xiàn)幾率最多的某一海拔高度進(jìn)行逐點(diǎn)校正。

1.2 滯后校正及飛行方向差校正

飛機(jī)上GPS 與磁探頭并不在同一位置,使得測(cè)量存在滯后偏差,需對(duì)數(shù)據(jù)進(jìn)行滯后校正。校正方法為利用反向重復(fù)線峰值之間的距離的一半與平均地速的比值求取滯后時(shí)間,計(jì)算公式如公式（1）。

式中：t-滯后校正時(shí)間；δS-磁異常的位置差；V-飛機(jī)每秒的平均地速。

通過(guò)對(duì)比圖2 中A、B 兩幅圖（為對(duì)比效果明顯截取部分區(qū)域作對(duì)比）可以發(fā)現(xiàn)通過(guò)滯后校正使得相鄰測(cè)線峰值趨勢(shì)平滑,“V”型誤差消失。

圖2 截取的磁場(chǎng)分布對(duì)比圖

磁探頭在不同方向通過(guò)同一點(diǎn)時(shí)測(cè)量值會(huì)有一定變化,因此需要進(jìn)行飛行方向差校正。根據(jù)不同方向通過(guò)同一點(diǎn)的飛行實(shí)驗(yàn)數(shù)據(jù)完成飛行方向差校正。通過(guò)對(duì)比圖2 中A、B 兩幅圖可以發(fā)現(xiàn)校正后圖中條帶性誤差得到很大的壓制,但仍存在一些系統(tǒng)噪聲。

1.3 常規(guī)調(diào)平

調(diào)平方法發(fā)展迅速,按照數(shù)據(jù)調(diào)平和網(wǎng)格化的順序可分為基于測(cè)線數(shù)據(jù)的調(diào)平和基于網(wǎng)格化數(shù)據(jù)的調(diào)平兩大類(lèi)。調(diào)平的基本原理為以最小二乘法為基礎(chǔ),使測(cè)線和切割線在其交點(diǎn)處的平方和最小。大致流程如下[7]：

1.3.1 求取切割線修正值

ΔTci為磁場(chǎng)平靜區(qū)第i 條切割線與所有和它相交的測(cè)線的交點(diǎn)的差值的算術(shù)平均值,n 為此切割線上與它所相交測(cè)線交點(diǎn)的數(shù)量,Tci為切割線原場(chǎng)值,Tj為交點(diǎn)處測(cè)線上的場(chǎng)值,以上場(chǎng)值的單位都是nT。

1.3.2 對(duì)切割線進(jìn)行修正

Tcij為切割線改正后的值,Tc0ij為切割線未修正前的原場(chǎng)值,單位為nT。

1.3.3 對(duì)測(cè)線進(jìn)行修正。

ΔTj與ΔTci數(shù)值相同,T0ij為測(cè)線未修正前的原場(chǎng)值,Tij為一次調(diào)平后得到的場(chǎng)值,若進(jìn)行反復(fù)迭代,可以使結(jié)果更為可靠,提高精度,達(dá)到設(shè)計(jì)的總精度標(biāo)準(zhǔn)。

1.4 微調(diào)平

經(jīng)過(guò)常規(guī)調(diào)平后的數(shù)據(jù)資料仍會(huì)存在一些水平差,需要經(jīng)過(guò)微調(diào)平處理。微調(diào)平采用頻率域和空間域組合濾波,在無(wú)切割線的條件下通過(guò)分離地質(zhì)信息與噪聲實(shí)現(xiàn)航磁資料的調(diào)平,在數(shù)據(jù)處理中起到了重要的作用。微調(diào)平步驟如下：

1.4.1 條帶狀干擾的提取

條帶干擾提取沿測(cè)線方向采用低通濾波、垂直測(cè)線方向采用高通濾波的方法。可采用巴特沃斯濾波結(jié)合余弦方向?yàn)V波結(jié)合的方法提取條帶干擾。

1.4.2 水平調(diào)整

將分離出的條帶干擾從測(cè)線中去除即可獲得微調(diào)平處理后的結(jié)果,若測(cè)線間仍存在一定的水平差,可重復(fù)微調(diào)平過(guò)程直至取得合理的效果。

通過(guò)常規(guī)調(diào)平和微調(diào)平的結(jié)果與原始數(shù)據(jù)對(duì)比（圖3）發(fā)現(xiàn),調(diào)平處理在保持原始資料分辨率的同時(shí),沒(méi)有移動(dòng)或歪曲有意義的地質(zhì)異常,去除了系統(tǒng)噪聲和系統(tǒng)誤差,提高了資料的信噪比。

圖3 磁場(chǎng)分布對(duì)比圖

2 調(diào)平方法在山陽(yáng)航磁資料處理中的應(yīng)用

山陽(yáng)地區(qū)位于秦嶺東西復(fù)雜構(gòu)造帶東段中部,區(qū)內(nèi)地層主要由泥盆系及石炭系以淺海相、濱海相的碎屑巖占優(yōu)勢(shì)的淺變質(zhì)巖系構(gòu)成區(qū)域地質(zhì)總體。區(qū)內(nèi)褶皺斷裂發(fā)育,中酸性小巖體廣泛分布,圍巖接觸變質(zhì)普遍,與中酸性小巖體有關(guān)的以Cu,Mo 為主的礦產(chǎn)及受泥盆系青石埡組層控的菱鐵多金屬礦產(chǎn)比較豐富,與小巖體及泥盆系地層關(guān)系密切的Au 異常廣泛分布[8-11]。

成礦類(lèi)型集中地表現(xiàn)為矽卡巖型和熱液型,其控礦因素主要為巖漿巖體、成礦圍巖與地質(zhì)構(gòu)造。在空間位置上,已發(fā)現(xiàn)的礦體或礦化體距離巖體一般不超過(guò)200m,且已知巖體周邊因圍巖或構(gòu)造強(qiáng)度不同,其礦化強(qiáng)度及規(guī)模雖不盡相同,但都發(fā)生了不同程度的礦化,說(shuō)明礦（化）體與巖體之間已是一種“有巖有礦”的相互對(duì)應(yīng)關(guān)系。巖漿巖體已然成為區(qū)內(nèi)礦體尋找最有力的標(biāo)志。因此利用磁法資料圈定巖體對(duì)進(jìn)一步找礦有很大幫助。筆者根據(jù)研究的調(diào)平方法對(duì)山陽(yáng)地區(qū)航磁資料進(jìn)行了處理,圈定了8 處巖體、分析出5 處斷裂（如圖4）,為該區(qū)找礦提供了磁法依據(jù)。

圖4 山陽(yáng)地區(qū)航磁等值線分布及巖體、斷裂推斷圖

3 結(jié)論

本文通過(guò)分析航磁資料的調(diào)平原理和應(yīng)用,以實(shí)際航磁資料處理為例,指出了調(diào)平處理的關(guān)鍵性步驟和技術(shù)手段（滯后校正、飛行方向差校正、常規(guī)調(diào)平及微調(diào)平）,討論了每項(xiàng)處理對(duì)結(jié)果的影響：其中滯后校正可以使相鄰測(cè)線峰值趨勢(shì)平滑,消除"V"型誤差,方向差校正可以壓制條帶性干擾,常規(guī)調(diào)平及微調(diào)平使數(shù)據(jù)的系統(tǒng)噪聲進(jìn)一步壓制,這些使我們深入理解了調(diào)平對(duì)航磁資料處理的實(shí)際意義。之后筆者把研究的技術(shù)方法應(yīng)用到山陽(yáng)航磁資料的處理中,圈定了幾處巖體,為尋找該區(qū)多金屬礦產(chǎn)提供了方向。