殷長春, 張博, 劉云鶴, 任秀艷， 齊彥福， 裴易峰，邱長凱， 黃鑫， 黃威， 繆佳佳, 蔡晶

吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院， 長春 130026

?

航空電磁勘查技術(shù)發(fā)展現(xiàn)狀及展望

殷長春, 張博*, 劉云鶴, 任秀艷， 齊彥福， 裴易峰，邱長凱， 黃鑫， 黃威， 繆佳佳, 蔡晶

吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院， 長春 130026

航空電磁作為一種高效的地球物理勘查技術(shù)手段，其發(fā)展在國外(加拿大、澳大利亞等國家)已十分成熟.然而，在我國該項技術(shù)仍處于發(fā)展當中，在國內(nèi)目前尚未形成具有實際探測能力的航空電磁系統(tǒng)和解釋手段.這一現(xiàn)狀嚴重影響了我國對地形地質(zhì)條件復(fù)雜區(qū)域(比如廣大西部地區(qū))礦產(chǎn)資源勘查的需求.本文旨在通過系統(tǒng)介紹航空電磁勘查技術(shù)中的基礎(chǔ)理論、關(guān)鍵技術(shù)、儀器系統(tǒng)、數(shù)據(jù)處理、解釋及應(yīng)用，并對未來我國航空電磁勘查技術(shù)的發(fā)展提出建議，使讀者了解該技術(shù)未來發(fā)展方向和研究熱點，以期該項技術(shù)在我國得到快速發(fā)展并獲得廣泛應(yīng)用.

礦產(chǎn)資源勘查； 航空電磁勘查技術(shù)； 頻域和時域航空電磁系統(tǒng)； 電磁數(shù)據(jù)處理； 數(shù)據(jù)成像和解釋

1 引言

航空電磁(AEM)是利用飛機平臺搭載電磁勘探設(shè)備的地球物理勘查技術(shù).由于該方法采用飛行平臺作為搭載裝置，無需地面人員接近勘查作業(yè)區(qū)，特別適合高山、沙漠、湖泊沼澤、森林覆蓋等地形復(fù)雜地區(qū).國外航空電磁勘查自20世紀五六十年代開始發(fā)展，目前已形成時間域和頻率域，主動源和被動源，固定翼和直升機吊艙平臺的系列航空電磁勘查系統(tǒng).應(yīng)用領(lǐng)域由傳統(tǒng)的金屬礦產(chǎn)和油氣資源勘查，發(fā)展到包括環(huán)境工程、地下水和地?zé)豳Y源、海洋地形調(diào)查、極地研究等應(yīng)用領(lǐng)域.國內(nèi)七八十年代曾發(fā)展過航空電磁系統(tǒng)，曾研制出了固定翼頻率域和時間域航空電磁系統(tǒng).90年代由于地質(zhì)行業(yè)蕭條航空電磁系統(tǒng)的開發(fā)基本處于停滯狀態(tài).近年，隨著國家經(jīng)濟高速發(fā)展，對礦產(chǎn)資源的需求急劇增加，而國內(nèi)地形地質(zhì)條件良好地區(qū)的礦產(chǎn)資源已基本勘查完畢.礦產(chǎn)勘查的目標相應(yīng)地轉(zhuǎn)向占我國國土面積近三分之二的廣大西部地形條件復(fù)雜地區(qū)(高山、沙漠)和東部大面積森林覆蓋、湖泊沼澤地區(qū).這些地區(qū)由于地面人員無法接近，地球物理觀測程度較低.人們把眼光逐漸轉(zhuǎn)向航空地球物理勘查，特別是航空電磁勘查技術(shù).然而，由于國外的高技術(shù)封鎖和壟斷，航空電磁關(guān)鍵技術(shù)(比如超大功率發(fā)射、大動態(tài)范圍多分量采集、吊艙動態(tài)監(jiān)測和控制技術(shù)、飛行平臺背景場去除等)無法取得突破性進展，使得國內(nèi)航空電磁技術(shù)發(fā)展一直停留在地面原理樣機階段，無法實現(xiàn)生產(chǎn)飛行.本文通過對航空電磁勘查技術(shù)及國內(nèi)外發(fā)展現(xiàn)狀、航空電磁理論、數(shù)據(jù)處理和反演解釋等做出系統(tǒng)介紹，并結(jié)合應(yīng)用實例，對航空電磁勘查原理與應(yīng)用進行分析，以期為我國航空電磁勘查關(guān)鍵技術(shù)的突破指明方向，從而推動我國航空電磁勘查向?qū)嵱没较虬l(fā)展，早日實現(xiàn)我國地形復(fù)雜地區(qū)礦產(chǎn)資源勘查.本文在第2節(jié)對航空電磁勘查技術(shù)進行總體介紹；第3節(jié)介紹國內(nèi)外航空電磁勘查技術(shù)發(fā)展現(xiàn)狀；第4節(jié)介紹航空電磁正反演理論；第5節(jié)介紹航空電磁勘查技術(shù)應(yīng)用；最后，對我國航空電磁技術(shù)發(fā)展做出展望.

2 航空電磁勘查技術(shù)

航空電磁勘查技術(shù)涉及到的基本概念較多，闡述這些概念是了解這種勘探方法的前提.下面我們對這些方法技術(shù)涉及到的基本概念做出介紹.

2.1 航空電磁勘查技術(shù)基本原理

航空電磁采用發(fā)射線圈發(fā)射電磁信號，通過接收線圈接收經(jīng)空氣傳播的一次場信號和經(jīng)地下介質(zhì)感應(yīng)產(chǎn)生的二次場信號.當?shù)叵虏淮嬖陔娦圆痪鶆蝮w時，接收線圈接收到的二次場信號為地下半空間的感應(yīng)信號(圖1a)；當?shù)叵麓嬖陔娦援惓ｓw時，異常體也參與了電磁感應(yīng)，導(dǎo)致二次場信號與不含異常體時出現(xiàn)了差異(圖1b).通過分析由地下介質(zhì)產(chǎn)生的電磁場，可以得到地下電性結(jié)構(gòu)分布信息.航空電磁通常分為頻率域和時間域兩種方法.頻率域航空電磁法利用發(fā)射線圈發(fā)射諧變電磁信號，接收線圈接收相同頻率的電磁信號；而時間域航空電磁法則利用發(fā)射線圈發(fā)射脈沖信號，接收線圈接收隨時間變化的瞬變信號.頻率域航空電磁根據(jù)發(fā)射和接收機的關(guān)系分為：水平共面(HCP)，直立共軸(VCA)，直立共面(VCP)和交叉裝置；而時間域航空電磁通常采用共中心或分離裝置.

2.2 發(fā)射波形

時間域航空電磁采用的發(fā)射波形種類很多，可大致分為階躍、脈沖、方波、半正弦波、三角波和梯形波，不同發(fā)射波形的激勵響應(yīng)強度有所不同.根據(jù)Liu(1998)，對斷電(off-time)電磁響應(yīng)來說，在發(fā)射信號寬度和強度相同的情況下，方波電磁響應(yīng)強于其他發(fā)射波形的電磁響應(yīng)且方波越寬，電磁響應(yīng)強度越大；多數(shù)波形存在最優(yōu)波形寬度(對應(yīng)最大響應(yīng)幅值的波形寬度)，半正弦波的最優(yōu)波形寬度為目標體時間常數(shù)的2.5倍.陳曙東等(2012)對自由空間中導(dǎo)電回線電磁響應(yīng)研究發(fā)現(xiàn)，當目標體時間常數(shù)較大時，不同發(fā)射波形的目標體電磁響應(yīng)幅值之比近似等于波形與時間坐標軸所構(gòu)成區(qū)域的面積之比.對供電(on-time)電磁響應(yīng)，當目標體時間常數(shù)較小時，方波電磁響應(yīng)強于其他發(fā)射波形；方波電磁響應(yīng)隨目標體時間常數(shù)的減小而增加，而其他波形隨目標體時間常數(shù)的減小趨于穩(wěn)定值；梯形波電磁響應(yīng)比方波小，但大于三角波和半正弦波.隨著時間常數(shù)的增加，方波激勵響應(yīng)持續(xù)減小，而半正弦波和三角波的激勵響應(yīng)先減小，而后趨于不變，并在時間常數(shù)較大時強度超過方波.因此，利用三角波和半正弦波的on-time信號易于探測地下低阻目標體(陳曙東等，2012).

從關(guān)于趨膚深度的討論我們知道，通過發(fā)射不同頻率的電磁信號，我們可以勘探不同深度的目標體，實現(xiàn)航空電磁的測深目標.同樣，對于時間域航空電磁系統(tǒng)，由于早期道信號高頻成分豐富，主要反映淺層信息；而晚期道信號中低頻成分豐富，主要反映深部地電信息.通過觀測全時時間域電磁信號，我們可以研究地下介質(zhì)從淺到深的電性分布特征.

2.3 On-time和off-time

如前文所述，時間域航空電磁發(fā)射脈沖信號，通過接收線圈接收瞬變電磁信號.發(fā)射裝置供電和斷電期間分別稱為on-time(供電) 和off-time(斷電)，對應(yīng)接收線圈采集的數(shù)據(jù)稱為on-time和off-time數(shù)據(jù).on-time數(shù)據(jù)包含了一次場和二次場信號，且一次場比二次場大得多；off-time數(shù)據(jù)是在發(fā)射線圈斷電情況下采集的，因此只含有地下介質(zhì)感應(yīng)的二次場.目前航空電磁實際解釋中，off-time數(shù)據(jù)應(yīng)用較多.隨著航空電磁數(shù)據(jù)處理技術(shù)發(fā)展，人們對on-time數(shù)據(jù)的重視程度逐漸增加(殷長春等，2015).必須指出的是飛機作為一個飛行平臺是一個大的良導(dǎo)體，在電磁信號中由飛機產(chǎn)生的背景場必須予以去除.

圖1 航空電磁勘查技術(shù)原理(a) 不存在異常體；(b) 存在異常體(Greg Hodges,2011,個人通訊).Fig.1 Airborne electromagnetic technology(a) Without abnormal body; (b) With abnormal body (Greg Hodges, 2011, Personal communication).

2.4 Fiducial、PPM、bucking技術(shù)

航空電磁數(shù)據(jù)采樣間隔通常為0.1 s，稱之為一個Fiducial(簡稱Fid).由于頻率域航空電磁一次場和二次場同時存在，且一次場比二次場信號強得多，通常系統(tǒng)通過bucking技術(shù)分離出一次場和二次場，同時將二次場和一次場作比值，并乘以106，稱為ppm.然而，對于時間域航空電磁，由于主要應(yīng)用off-time數(shù)據(jù)，或者在on-time數(shù)據(jù)中剔除一次場，所以觀測信號直接使用磁感應(yīng)信號的單位(T/s，pT/s，fT/s等).航空電磁接收線圈接收到的信號中，能夠反映地下介質(zhì)信息的僅二次場信號，因此在對采集數(shù)據(jù)進行處理之前必須去除一次場.頻率域航空電磁采用bucking技術(shù)，而在時間域中通常采用互相關(guān)技術(shù)予以去除.

2.5 全區(qū)視電阻率概念

視電阻率是地下電性特征的綜合反映.它定義為可以替換地下不均勻介質(zhì)的均勻半空間電阻率，如果該均勻半空間對于相同的發(fā)射接收條件，在接收機處產(chǎn)生相同的電磁信號.電磁勘查中傳統(tǒng)的早期和晚期視電阻率難以給出連續(xù)而且準確的電性參數(shù)，為此人們提出了全區(qū)視電阻率的概念，即全時段/全區(qū)域視電阻率.殷長春和樸華榮(1991)研究了電偶源地面電磁全區(qū)視電阻率定義.他們發(fā)現(xiàn)由于磁場是均勻半空間電阻率的單值函數(shù)，使用磁場定義全區(qū)視電阻率優(yōu)于磁場時間導(dǎo)數(shù)，并利用多項式擬合電阻率關(guān)于磁場的反函數(shù)定義了全區(qū)視電阻率.強建科等(2010)研究了適用于航空電磁多分量觀測全時域電阻率計算方法，張瑩瑩(2013)使用牛頓迭代方法進行電偶源地空系統(tǒng)全域視電阻率定義，戚志鵬等(2014)對大定源瞬變電磁三分量全域視電阻率進行了定義.他們的視電阻率定義同樣是基于航空電磁響應(yīng)對于半空間電阻率的單值性，并通過對磁場關(guān)于半空間電阻率函數(shù)進行泰勒展開和線性近似，由假設(shè)的初始視電阻率值利用迭代方法計算全區(qū)視電阻率.

2.6 Footprint

航空電磁勘探數(shù)據(jù)量龐大，導(dǎo)致航空電磁數(shù)據(jù)處理解釋耗時驚人.幸運的是，航空電磁測量系統(tǒng)是緊湊系統(tǒng)，對于單個測點來說“敏感區(qū)域”遠小于整個測區(qū).這個“敏感區(qū)域”就是航空電磁系統(tǒng)的“footprint”.航空電磁中footprint的第一個明確定義由Liu和Becker(1990)給出，此后很多學(xué)者對footprint進行了研究(Kovacs et al.，1995；Reid and Vrbancich，2004; Reid et al.,2006；Beamish，2003).Yin等(2014)將footprint定義為均勻?qū)щ姲肟臻g中的最小區(qū)域，該區(qū)域邊界上感應(yīng)電流在接收機處產(chǎn)生相同的二次場，同時該區(qū)域中地下感應(yīng)電流在接收機處產(chǎn)生的二次場占整個半空間產(chǎn)生的二次場的90%.研究發(fā)現(xiàn)，航空電磁系統(tǒng)footprint主要取決于飛行高度.通常情況下，footprint為飛行高度的3～4倍.利用footprint的概念，我們在從事航空電磁數(shù)據(jù)反演時，可先基于footprint大小對測區(qū)進行單元劃分，在完成各測區(qū)單元的電磁數(shù)據(jù)反演后，集成得到測區(qū)總的電性分布.這為復(fù)雜二、三維模型反演節(jié)省大量計算成本.

2.7 良導(dǎo)和高阻體頻率域和時間域電磁信號特征

導(dǎo)電介質(zhì)中頻率域電磁信號在小感應(yīng)數(shù)(對應(yīng)于高阻或低頻)時，電磁信號的實部比虛部小；隨著感應(yīng)數(shù)的增大，電磁信號的實部逐漸增大直至達到飽和，而虛分量經(jīng)過極大值后逐漸衰減為零.在感應(yīng)數(shù)很大時(對應(yīng)于低阻或高頻)，電磁信號的實分量比虛分量大得多.相比之下，時間域電磁信號對于高阻體在斷電早期很強，然而隨時間衰減很快；而對于良導(dǎo)體，斷電早期電磁信號幅值相對較小，但隨時間衰減緩慢.因此，良導(dǎo)體對應(yīng)的晚期道信號較強.利用時域/頻域電磁信號的變化特征，我們能從航空電磁觀測數(shù)據(jù)中有效識別地下高阻和良導(dǎo)體.

圖2給出了頻域和時域航空電磁信號特征及對高阻和低阻體識別的應(yīng)用實例.圖中A處地下存在高阻體，而C處為良導(dǎo)體.從圖可以看出，高阻介質(zhì)(A處)的時間域電磁信號在早期較強，晚期電磁響應(yīng)衰減較快，表現(xiàn)為各時間道信號非常稀疏.相應(yīng)地，頻率域信號在各頻段虛部均大于實部；良導(dǎo)體(C處)的時間域電磁信號在早期相對較弱，在晚期衰減較慢，表現(xiàn)為晚期各時間道信號比較密集.通過分析時域和頻域航空電磁數(shù)據(jù)的特征，可以定性識別地下高阻和低阻體.

2.8 地下介質(zhì)磁性、介電常數(shù)、各向異性對航空電磁響應(yīng)的影響數(shù).利用低頻電磁信號求解地下介質(zhì)的磁化率，利用高頻電磁信號求解地下介質(zhì)的介電常數(shù)，而電阻率則可利用所有頻點的電磁信號進行求解.

Huang 和Fraser(2001)討論地下介質(zhì)磁性和介電常數(shù)對航空電磁響應(yīng)的影響特征，提出了利用航空電磁數(shù)據(jù)求解地下介質(zhì)的電磁特征和介電參

圖2 良導(dǎo)體C和高阻體A航空電磁響應(yīng)特征(Greg Hodges, 2011,個人通訊)Fig.2 AEM responses for conductive body C and resistive body A (Greg Hodges, 2011, Personal communication)

航空電磁勘探信號對地下介質(zhì)各向異性也相當敏感.Yin和Fraser(2004a)研究一維各向異性介質(zhì)對航空電磁響應(yīng)的影響特征，得出如下結(jié)論：對電磁信號實分量，介質(zhì)各向異性的影響隨頻率的增加而變大，各向異性影響的最大值出現(xiàn)在高頻段；對信號虛分量，各向異性影響隨頻率增大先增加，到達某一頻率后開始變小，各向異性影響的最大值出現(xiàn)在中間頻段.對于VCA裝置，當飛行方向改變時，各向異性對電磁響應(yīng)的影響特征發(fā)生明顯變化；而對于HCP裝置，各向異性對電磁信號的影響不受飛行方向的影響.利用VCA裝置電磁響應(yīng)隨飛行方向的變化特征，我們可以有效地識別地下各向異性的主電導(dǎo)率方向.

3 航空電磁勘查技術(shù)國內(nèi)外發(fā)展現(xiàn)狀

航空電磁勘查技術(shù)已經(jīng)有大約70年的發(fā)展歷史， 經(jīng)歷了從簡單的單頻系統(tǒng)逐步發(fā)展成為多頻電磁系統(tǒng)、時間域系統(tǒng)、直升機和固定翼系統(tǒng)等.本節(jié)首先介紹航空電磁在國外的發(fā)展現(xiàn)狀，特別對目前具有國際領(lǐng)先水平的航空電磁系統(tǒng)和相關(guān)技術(shù)進行介紹，然后再回顧國內(nèi)航空電磁的發(fā)展歷史和現(xiàn)狀.

3.1 國外航空電磁發(fā)展現(xiàn)狀

1948年，固定翼航空電磁系統(tǒng)Stanmac-McPhar在加拿大試飛成功，這標志著第一個航空電磁勘探系統(tǒng)的誕生(Fountain，1998).1954年，航空電磁勘探方法在加拿大New Braunswick省發(fā)現(xiàn)了Health Steele 礦床，這一發(fā)現(xiàn)極大地促進了航空電磁勘查技術(shù)和系統(tǒng)的發(fā)展.1955年，第一個吊艙式硬支架直升機航電系統(tǒng)誕生.此后，Tony Barringer開發(fā)了第一個時間域固定翼INPUT系統(tǒng).該系統(tǒng)后經(jīng)歷不斷更新，一直沿用至今. Slichter(1955)闡明使用天然場代替人工場源可以大幅度提高航空電磁系統(tǒng)的探測深度.基于這一理論，第一個被動源航空電磁系統(tǒng)AFMAG于1958年在加拿大研制成功.

在隨后的幾十年里，航空電磁系統(tǒng)得到不斷的優(yōu)化和改進.1965年，INPUT系統(tǒng)升級成為Mark V INPUT系統(tǒng)，增大了發(fā)射磁矩，進而大幅度增加了系統(tǒng)的有效探測深度.1967年，一個新的多頻吊艙式航空電磁系統(tǒng)F-400被成功研發(fā).與之前的系統(tǒng)相比，該系統(tǒng)采用水平偶極作為發(fā)射源，這種裝置對飛機的改造最小，便于在不同飛行裝置上移植.1969年，人們通過改造地面Turam系統(tǒng)研發(fā)出Turair半航空系統(tǒng).該系統(tǒng)通過利用地面大線圈發(fā)射增大了有效勘探深度.1970年，單頻多分量接收DIGHEM I型航空電磁系統(tǒng)研發(fā)成功.該系統(tǒng)采用X軸方向發(fā)射，接收X、Y、Z三分量.此外，Questor對Mark VI INPUT系統(tǒng)進行升級，由此生產(chǎn)了兩個新的勘探平臺——Skyvan(1971)和Trislander(1973)系統(tǒng).同時，1970年前后印度開發(fā)了可以接收B和dB/dt的固定翼時間域航空電磁系統(tǒng)，然而該系統(tǒng)只做了試驗飛行.此后，McPhar對F-400做出進一步改進，將其升級為多頻電磁系統(tǒng)并安裝在固定翼和直升機飛行平臺上.1976年，固定翼COTRAN和Smelting EM-30系統(tǒng)研發(fā)成功.COTRAN系統(tǒng)采用了INPUT系統(tǒng)的基本結(jié)構(gòu)，發(fā)射波形為方波，接收X和Z兩個分量，然而該系統(tǒng)沒有實現(xiàn)商業(yè)飛行.Smelting EM-30系統(tǒng)通過采用大收發(fā)距的共軸裝置達到較大的勘探深度.70年代末期，大收發(fā)距頻率域固定翼航空電磁系統(tǒng)基本被淘汰，固定翼航空電磁以時間域系統(tǒng)為主.直升機系統(tǒng)則朝著多頻發(fā)射、多分量接收方向發(fā)展.1976年，在DIGHEM I的基礎(chǔ)上，DIGHEM II成功研發(fā)，并通過系統(tǒng)改進實現(xiàn)雙頻發(fā)射(900 Hz和3600 Hz)和多分量接收.

20世紀80年代中期，由于各大公司削減礦產(chǎn)勘查預(yù)算，航空電磁勘查技術(shù)走向蕭條.然而，航空電磁系統(tǒng)的研發(fā)并未因此停滯不前.1983年，Geoterrex引進了一個新的INPUT航空電磁勘探平臺CASA 212，并于1985年研發(fā)成功針對這一平臺，且與INPUT系統(tǒng)具有相似結(jié)構(gòu)的GEOTEM系統(tǒng).相比于固定翼航空電磁系統(tǒng)受到的沖擊，基于直升機平臺的航空電磁系統(tǒng)得到迅速發(fā)展，很多新的直升機航空電磁系統(tǒng)研發(fā)成功.1982年，Questor公司成功研發(fā)直升機INPUT系統(tǒng)，以期在地形復(fù)雜地區(qū)獲得好的飛行數(shù)據(jù)；加州大學(xué)伯克利分校成功研發(fā)UNICOIL直升機系統(tǒng)(UNICOIL cryogenic helicopter system).該系統(tǒng)采用超導(dǎo)UNICOIL作為發(fā)射和接收線圈，工作頻率為40 Hz，具有較大的探測深度.然而，該系統(tǒng)由于需要大量液氦而飛機存儲空間有限，缺乏實用價值.到20世紀90年代，航空電磁系統(tǒng)主要發(fā)展方向為：1)多頻、高分辨率的直升機吊艙系統(tǒng)；2)低頻固定翼時間域系統(tǒng).

經(jīng)歷了航空電磁勘查技術(shù)的蕭條，20世紀90年代礦產(chǎn)資源勘探開始回暖，給航空電磁帶來了新活力.固定翼時間域航空電磁系統(tǒng)基本實現(xiàn)了三分量接收、可選擇發(fā)射波形寬度和頻率、接收B和dB/dt信號等功能；而直升機吊艙式頻率域系統(tǒng)基本實現(xiàn)了多頻發(fā)射和接收，且發(fā)射和接收線圈多安裝在6～9 m的吊艙中.其中，90年代早期，World Geoscience開發(fā)了SALTMAP系統(tǒng).該系統(tǒng)是第一個專門為近地表電阻率成像而設(shè)計的固定翼時間域航電系統(tǒng).同時，Geotech公司引入HUMMING BIRD 直升機航電系統(tǒng).該系統(tǒng)發(fā)射和接收線圈安裝在4.5 m的吊艙中，以增加系統(tǒng)的靈活性；而Elliott Geophysics International 公司研發(fā)出一種新型半航空系統(tǒng)FLAIRTEM.該系統(tǒng)本質(zhì)上是TURAIR系統(tǒng)的時間域版本.90年代晚期，航空電磁系統(tǒng)研發(fā)如雨后春筍般蓬勃興起.其中，直升機系統(tǒng)包括1997年研發(fā)成功的Anglo ExplorHEM、 AERODAT HELITEM和HELI QUESTEM等，1998年研發(fā)成功的THEM和DIGHEM RES BIRD，1999年研發(fā)成功的AeroTEM、 NewTEM和SIAL PHOENIX等；固定翼系統(tǒng)包括1997年研發(fā)成功的GEOTEM， 1998年研發(fā)成功的MEGATEM和QUESTEM，1999年研發(fā)成功的GEOTECH HAWK系統(tǒng)(Thomson et al., 2007).

進入21世紀，人們對航空電磁系統(tǒng)和勘查技術(shù)提出了更高要求.固定翼系統(tǒng)主要朝著具有更大發(fā)射功率和勘探深度的方向發(fā)展，而直升機系統(tǒng)主要朝著多線圈、寬頻帶、高分辨率和高精度方向發(fā)展.在此期間研發(fā)的航空電磁系統(tǒng)是目前各航空地球物理公司應(yīng)用的主力軍(表1)，主要有Fugro公司2001年研發(fā)的時間域固定翼MEGATEM II系統(tǒng)和2002年的頻率域直升機吊艙系統(tǒng)RESOLVE，2005年研發(fā)的時間域直升機吊艙系統(tǒng)HeliGEOTEM和2008年研發(fā)的HeliTEM系統(tǒng)； Geotech公司2004年研發(fā)的時間域直升機吊艙系統(tǒng)VTEM； Aeroquest公司2003年研發(fā)的時間域直升機吊艙系統(tǒng)AeroTEM，SkyTEM Surveys公司2003年研發(fā)的系列時間域直升機吊艙系統(tǒng)SkyTEM(Thomson et al., 2007).這些系統(tǒng)各具特色，承擔了全球超過90%航空地球物理勘查任務(wù).

3.2 國際主要航空電磁勘探公司

目前，國際上規(guī)模較大的航空電磁勘探公司有：加拿大Fugro(現(xiàn)為CGG)、Geotech、Aeroquest(現(xiàn)已被Geotech收購)和丹麥SkyTEM Surveys公司.

表1 直升機和固定翼航空電磁系統(tǒng)研發(fā)歷史(參考Thomson et al., 2007)
Table 1 Development history of helicopter and fixed wing AEM system (After Thomson et al., 2007)

這些公司擁有強大的技術(shù)實力和先進的勘探設(shè)備，占領(lǐng)了全球航空地球物理勘查市場.

1) CGG公司

CGG公司成立于1931年，現(xiàn)已發(fā)展成為世界領(lǐng)先的地球物理勘探公司之一.2013年，CGG Veritas收購了Fugro的地球科學(xué)部.航空地球物理勘探是CGG提供的全球技術(shù)服務(wù)之一.該公司擁有的航空電磁系統(tǒng)包括DIGHEM、RESOLVE、GEOTEM、MEGATEM、HeliGEOTEM、HeliTEM、MULTIIMPULSE等.這些系統(tǒng)涵蓋了直升機和固定翼、時間域和頻率域等，使CGG公司能夠出色地完成各種地質(zhì)條件和目標的勘探任務(wù).

DIGHEM和RESOLVE是全數(shù)字頻率域直升機系統(tǒng)，能夠提供先進的校準技術(shù)、較小的零漂、較低噪聲水平以及實時信號處理技術(shù).DIGHEM系統(tǒng)擁有三對水平共面和兩對直立共軸線圈、頻率覆蓋范圍大(900～56 kHz)，確保了該系統(tǒng)對地質(zhì)體有很大的靈敏度，適用于較淺目標體勘查(<120 m)，主要用于勘探淺部礦產(chǎn)資源、環(huán)境工程、地下水、海侵和極地研究等.RESOLVE系統(tǒng)包含五對水平共面(頻率400～140 kHz)和一對直立共軸線圈(頻率3300 Hz).該系統(tǒng)憑借著較大的頻率覆蓋范圍、較高水平和垂向分辨率、高采樣率和實時信號處理技術(shù)，廣泛應(yīng)用于構(gòu)造填圖、礦產(chǎn)、地?zé)岷偷叵滤Y源勘查.GEOTEM和MAGTEM系統(tǒng)是固定翼時間域航空電磁系統(tǒng).GEOTEM是裝載于雙引擎固定翼飛機CASA212上的時間域電磁勘查系統(tǒng)，發(fā)射磁矩可達1百萬Am2，發(fā)射基頻可調(diào)，可進行寬頻帶on-和off-time多分量觀測(Annan and Lockwood，1991).目前，該系統(tǒng)主要用于金屬礦、金剛石、鈾礦、油氣、地下水和其他地質(zhì)填圖.MEGATEM是基于GEOTEM系統(tǒng)研發(fā)的，采用四引擎固定翼飛機Dash-7作為飛行平臺(Smith et al., 2003).其發(fā)射磁矩超過2百萬Am2，是目前世界上發(fā)射功率和勘探深度最大的固定翼航空電磁系統(tǒng).它可進行寬頻帶on-和off-time多分量觀測.目前主要用于大測區(qū)深部目標勘查.HeliGEOTEM和HeliTEM系統(tǒng)是在GEOTEM基礎(chǔ)上發(fā)展起來的直升機吊艙式時間域航電系統(tǒng).HeliGEOTEM匯集了GeoTEM和MEGATEM系統(tǒng)的先進技術(shù)，具有更好的靈活性和更高的分辨率(Fountain et al.，2005).HeliTEM系統(tǒng)在HeliGEOTEM基礎(chǔ)上進行了較大改進.首先，由于采用銅線代替鋁管作為發(fā)射線圈，減輕發(fā)射線圈重量，增加發(fā)射磁矩，已達2百萬Am2；其次，由于將多分量接收機直接置于發(fā)射線圈上方(HeliGEOTEM系統(tǒng)采用獨立吊艙)，增加了系統(tǒng)對地下導(dǎo)電介質(zhì)的靈敏度；最后，HeliTEM系統(tǒng)可進行寬頻帶on-和off-time多分量觀測.由于結(jié)合直升機系統(tǒng)的高橫向分辨率和時間域系統(tǒng)大發(fā)射功率和大勘探深度的優(yōu)點，目前該系統(tǒng)廣泛應(yīng)用于礦產(chǎn)、油氣、地下水資源勘查和構(gòu)造填圖.由于所有CGG時間域航空電磁系統(tǒng)均采用三分量接收，采集到的二次場數(shù)據(jù)具有很低的噪聲水平，使得這些系統(tǒng)具有很高的探測靈敏度和很大的勘探深度.MULTIPULSE時間域航空電磁系統(tǒng)是由CGG現(xiàn)有的時間域電磁系統(tǒng)改進而成，可以使用直升機或固定翼兩種飛行平臺.該系統(tǒng)通過在前端先發(fā)射一個高能量的半正弦波，在末端再發(fā)射一個快速關(guān)斷的方波或梯形波，實現(xiàn)系統(tǒng)既能有很大的勘探深度，有利于深部目標體探測，又具有較高的近地表分辨率的特點(Chen et al.，2014).

基于眾多優(yōu)良的航空電磁勘探系統(tǒng)，CGG在全球范圍內(nèi)出色完成了各種航空電磁勘探任務(wù).在非洲南部, CGG公司使用DIGHTEM對煤田進行了勘查.2012年CGG公司使用HeliTEM在加拿大的Lalor湖對火山成因的硫化物礦床進行了勘探(Yang and Oldenburg，2013).在加拿大，該公司使用MULTIPULSE系統(tǒng)對Athabasca油砂進行了勘探.在Texas州西部，MEGATEM II系統(tǒng)被用于調(diào)查地下水鹽化度(Paine and Collins，2003).在德國北部海岸，RESOLVE系統(tǒng)對海侵情況進行了調(diào)查(Wiederhold et al.，2010).在澳大利亞西部，GEOTEM DEEP系統(tǒng)被用于鎳礦勘查(Wolfgram and Golden，2001).此外，CGG公司的航空電磁系統(tǒng)在尋找金礦、鈾礦及環(huán)境和工程勘探方面都取得較好的應(yīng)用效果.

2) Geotech公司

Geotech公司已有超過三十年的發(fā)展歷史.過去三十年中，Geotech在不斷完善已有的系統(tǒng)基礎(chǔ)上，積極開發(fā)新系統(tǒng)以提高服務(wù)質(zhì)量.1982年，Geotech開發(fā)了四頻直升機航電系統(tǒng)，并在加拿大內(nèi)陸水域進行勘探飛行；1983年，Geotech開發(fā)了第一個海水探測系統(tǒng)，并完成了NORDA水深勘查項目；1985年，Geotech開發(fā)了直升機航電系統(tǒng)，并于1987年為CRREL在阿拉斯加完成了冰層厚度勘查；1992年，Geotech開發(fā)了輕量級多頻數(shù)字式直升機航電系統(tǒng)；1998年，Geotech開發(fā)了可現(xiàn)場編程多頻固定翼航電系統(tǒng)；2000年，Geotech開發(fā)了全數(shù)字被動源電磁勘探系統(tǒng)；2002年，Geotech開發(fā)了全數(shù)字時間域航空電磁系統(tǒng)VTEM(Witherly et al.，2004)；2006年，Geotech使用ZTEM系統(tǒng)進行首次地球物理勘查飛行；2007年，Geotech第一次使用ZTEM系統(tǒng)進行商業(yè)性礦產(chǎn)勘查(Witherly and Sattel，2013)；2010年，該公司研發(fā)了固定翼ZTEM系統(tǒng)——FWZTEM.

Geotech公司提供的航空電磁勘查服務(wù)主要基于他們的VTEM系統(tǒng)(多功能時域電磁系統(tǒng))和ZTEM系統(tǒng).VTEM系統(tǒng)采用共中心、垂直偶極收發(fā)裝置以產(chǎn)生對稱的系統(tǒng)響應(yīng)(Witherly et al., 2004).這種收發(fā)裝置確保任何不對稱電磁響應(yīng)都是由地下異常體產(chǎn)生的，而非系統(tǒng)傾斜或是飛行方向?qū)е碌模瑥亩贡鎰e異常體位置和分析電磁數(shù)據(jù)更加容易.此外，VTEM系統(tǒng)采用了低噪聲接收線圈和大功率發(fā)射線圈，使得該系統(tǒng)具有很低的噪聲水平.由于具有以上優(yōu)點，VTEM系統(tǒng)深受顧客的青睞，到目前為止該系統(tǒng)已經(jīng)飛行超過200萬測線公里.ZTEM系統(tǒng)是在AFMAG系統(tǒng)基礎(chǔ)上研發(fā)的頻率域電磁勘查系統(tǒng).該系統(tǒng)可以搭載在直升機或固定翼兩種飛行平臺上.與其他商業(yè)電磁系統(tǒng)不同，該系統(tǒng)使用電離層電流或自然界產(chǎn)生的25～720 Hz的雷電信號作為激發(fā)場源(Witherly and Sattel, 2012).由于采用了特制的接收裝置和先進的信號處理技術(shù)，該系統(tǒng)擁有較低的噪聲、較高的分辨率和較大的勘探深度.ZTEM在過去四年時間里完成了超過25萬公里的勘探任務(wù).

利用VTEM和ZTEM系統(tǒng)，Geotech公司在全球完成了很多勘查任務(wù)(Witherly et al., 2004; Witherly and Sattel, 2012, 2013; Kaminski and Oldenburg, 2012).在坦桑尼亞Victoria地區(qū)，Geotech公司使用VTEM系統(tǒng)對金礦進行了勘查；在North Arizona，該系統(tǒng)被用于對老鈾礦區(qū)勘查；在加拿大Québec省，該系統(tǒng)被用于對Caber礦床進行勘查.在加拿大的Saskatchewan，Geotech公司使用ZTEM對鈾礦進行了勘查；在Nevada東南部，該系統(tǒng)被用于對金銀礦床進行勘查；2007年在加拿大安大略省，該系統(tǒng)被應(yīng)用于PGM-Cu-Ni礦藏勘探.

3) Aeroquest公司

Aeroquest航空地球物理公司曾經(jīng)擁有包括十架固定翼飛機和超過十套時間域航電系統(tǒng).憑借其全球業(yè)務(wù)能力，Aeroquest在世界各地尋找礦產(chǎn)和油氣資源，提供準確的近地表成像結(jié)果.Aeroquest公司的時間域航空電磁勘探主要基于他們的直升機AeroTEM系統(tǒng)，該系統(tǒng)在全球范圍內(nèi)完成超過50萬測線公里的飛行觀測.通過不斷完善自己的系統(tǒng)，Aeroquest已經(jīng)開發(fā)了三代AeroTEM航電系統(tǒng)(AeroTEM II，AeroTEM IV，AeroTEM HD).AeroTEM系統(tǒng)具有以下特色：采用了聚焦footprint技術(shù)，可以提供地下導(dǎo)體的電性特征信息，圈定異常體范圍；能夠采集on-time和早期off-time數(shù)據(jù)，滿足弱導(dǎo)電地質(zhì)體勘探和電導(dǎo)率成像的需求；擁有大功率發(fā)射(在30 Hz時可達到1百萬Am2)，接收信號具有較高的信噪比，保證了較大的勘探深度(勘探深度可以達到500 m)；同時記錄X和Z分量數(shù)據(jù)，對垂直地質(zhì)體十分敏感.基于以上優(yōu)勢，AeroTEM時間域航空電磁系統(tǒng)被廣泛應(yīng)用于金屬礦和金剛石勘查、油砂、地下水調(diào)查等方面(Balch et al., 2003).

Aeroquest公司的頻率域航空電磁系統(tǒng)包括Impulse和GEM2A.系統(tǒng)具有如下特色：發(fā)射信號頻帶較寬(330～96 kHz)，具有很好的電導(dǎo)率成像功能，對弱導(dǎo)電性地質(zhì)體具有很高的探測靈敏度，勘探深度可達100 m.Aeroquest公司頻率域航空電磁勘探系統(tǒng)可提供視電阻率和視深度計算、離散電磁異常體篩選與分級、電磁異常體的數(shù)值模擬及系列成像產(chǎn)品、頻率域響應(yīng)譜分析等.Aeroquest頻率域航空電磁系統(tǒng)被成功應(yīng)用于金屬礦和金剛石勘探、覆蓋層厚度、地下管網(wǎng)、地質(zhì)填圖和工程勘察等.Aeroquest公司于2012年被Geotech公司收購.

4) SkyTEM Surveys公司

SkyTEM Surveys是國際領(lǐng)先的航空地球物理公司之一.該公司可以采集時間域電磁、磁和放射性數(shù)據(jù)(radiometric)，并提供先進的數(shù)據(jù)處理技術(shù).SkyTEM能夠為客戶提供不同的系統(tǒng)選擇.每一個系統(tǒng)都具有不同的特點，可以滿足不同勘查需求.SkyTEM的勘查系統(tǒng)主要有SkyTEM101、SkyTEM302、SkyTEM304、SkyTEM508、SkyTEM512.SkyTEM101系統(tǒng)專為近地表(50～120 m)電阻率精確成像設(shè)計.該系統(tǒng)采用了符合空氣動力學(xué)的框架作為載體，降低了系統(tǒng)的振動噪聲，可以采集off-time的超早期無偏差(bias-free data)數(shù)據(jù)(斷電后3 μs)，從而得到近地表高分辨率解釋結(jié)果. SkyTEM302系統(tǒng)是針對埋深150～200 m的地質(zhì)體成像設(shè)計的.由于具有快速的關(guān)斷時間(4 μs)、較低的噪聲水平及超早期時間道數(shù)據(jù)記錄，SkyTEM302是解決近地表問題、高阻地層問題的理想選擇.SkyTEM304從2004年已成功應(yīng)用于地下水、礦產(chǎn)資源、油氣勘探以及巖土工程等領(lǐng)域.該系統(tǒng)最初配置了X和Z分量電磁接收線圈.2007年，系統(tǒng)進行了更新升級，添加了與發(fā)射機同步的總磁場觀測技術(shù).該系統(tǒng)把所有的傳感器安置在承載框架上，使得這些傳感器與地表距離較近，從而確保了采集的數(shù)據(jù)具有較高的精度.SkyTEM508系統(tǒng)的開發(fā)旨在獲得具有較高信噪比的晚期道數(shù)據(jù).該系統(tǒng)與SkyTEM304在傳感器安置方面采取了相同措施，確保了采集的數(shù)據(jù)具有較高精度.該系統(tǒng)可以在飛行過程中進行遠程設(shè)置，在35 ms的時間道仍然可以采集高質(zhì)量off-time數(shù)據(jù).SkyTEM512是SkyTEM開發(fā)的深部航空電磁探測系統(tǒng).由于擁有750000 Am2的發(fā)射磁矩，該系統(tǒng)的勘探深度可達500 m.此外，該系統(tǒng)采用獨特的安置在剛性框架上的零耦合接收線圈，使其具有較高的信噪比，從而進一步提高了系統(tǒng)的勘探深度和分辨率.

3.3 國內(nèi)航空電磁發(fā)展現(xiàn)狀

我國航空電磁系統(tǒng)的發(fā)展起步較早，然而由于歷史原因目前技術(shù)還極不成熟.從20世紀50年代末開始，原地礦部物化探所開始研制長導(dǎo)線半航空電磁探測儀器，后轉(zhuǎn)入硬支架頻率域航電系統(tǒng).1970—1974年，國土資源部航空物探遙感中心(原航空物探大隊)開展過天電系統(tǒng)的研制，因空難事故中止；1976—1983年原長春地質(zhì)學(xué)院開展過固定翼時域(M-1)系統(tǒng)的研制和試飛，并在黑龍江物探隊和湖北航空物探隊應(yīng)用于生產(chǎn)飛行；1976年北京地質(zhì)儀器廠開展直升機時間域航電系統(tǒng)研制；1974—1980年，桂林冶金地質(zhì)研究所開展直升機時間域系統(tǒng)研制；1981—1983年原長春地質(zhì)學(xué)院在對M-1系統(tǒng)改進基礎(chǔ)上研制M-2型固定翼時間域系統(tǒng)，由于缺少經(jīng)費中途停止.20世紀90年代，由于地質(zhì)行業(yè)經(jīng)歷了轉(zhuǎn)型和不景氣，航空地球物理，特別是航空電磁勘查系統(tǒng)和技術(shù)研發(fā)處于停滯狀態(tài).進入21世紀，隨著國家經(jīng)濟高速發(fā)展和對能源和礦產(chǎn)資源的需求激增，地質(zhì)行業(yè)迎來春天！航空地球物理電磁勘查技術(shù)和儀器系統(tǒng)研發(fā)再度受到相關(guān)部門的高度重視.目前國內(nèi)自主研發(fā)的系統(tǒng)主要包括中國地質(zhì)科學(xué)院地球物理地球化學(xué)勘查所研發(fā)的固定翼三頻航空電磁系統(tǒng)HDY-402和國土資源航空物探遙感中心研發(fā)的CHTEM時間域直升機航空電磁勘查系統(tǒng).HDY-402系統(tǒng)采用直立共面發(fā)射和接收裝置，收發(fā)距為19.2 m，工作頻率為463 Hz、1563 Hz和8333 Hz，最大采樣率8次/秒，噪聲水平<20 ppm，零漂小于100 ppm/小時，勘探深度為100 m.由于受飛行高度影響，信號較弱，該系統(tǒng)應(yīng)用受到限制，主要用于淺部礦產(chǎn)和水資源勘查.相比之下，我國自主研發(fā)的時間域航空電磁系統(tǒng)由于技術(shù)原因，發(fā)射功率較小(～28萬Am2)，目前未能在我國礦產(chǎn)資源勘查中獲得應(yīng)用.近年，國內(nèi)相關(guān)單位從國外引進航電系統(tǒng)，主要有中國冶金地質(zhì)局從加拿大引進的TS150系統(tǒng)和中國科學(xué)院電子所從俄羅斯引進的Impulse_A5系統(tǒng).TS150時間域直升機系統(tǒng)最高發(fā)射偶極矩為155000 Am2，可以記錄全波形B和dB/dt，基頻25～150 Hz，可多分量數(shù)據(jù)采集.Impulse_A5時間域直升機系統(tǒng)最大發(fā)射磁矩1.6×105Am2，發(fā)射基頻為25 Hz和75 Hz，懸掛裝置距離地表高度30～50 m，最大勘探深度可達400 m.另外，航遙中心還曾于2010—2011年間從加拿大Aeroquest公司引進AeroTEM-IV時間域航電系統(tǒng)，然而由于該系統(tǒng)發(fā)射磁矩小(<40萬Am2)，勘探深度淺，至今未能在我國礦產(chǎn)資源勘查中發(fā)揮作用.

我國航空電磁勘探的另一個途徑是利用國外地球物理公司提供的飛行觀測服務(wù).這種途徑由于造價過于昂貴使用不多.國內(nèi)有少數(shù)單位目前從事半航空電磁勘查技術(shù)和系統(tǒng)研發(fā)，然而，由于該技術(shù)發(fā)射源與接收機距離較遠，存在嚴重的體積效應(yīng)，喪失了常規(guī)航空電磁勘查系統(tǒng)的緊湊性(compact)和高目標分辨率的特點，可以預(yù)計該系統(tǒng)在我國礦產(chǎn)資源勘查中將不會發(fā)揮太大作用.

航空電磁具有很長的發(fā)展歷史，國際上的航空電磁勘查系統(tǒng)相當完善，技術(shù)已經(jīng)十分成熟.然而，我國的航空電磁勘查系統(tǒng)的研發(fā)還處于起步階段，尚未成功研發(fā)出能夠用于地形復(fù)雜地區(qū)和深部礦產(chǎn)資源勘查的航空電磁探測系統(tǒng).為了滿足我國深部和廣大地形復(fù)雜地區(qū)找礦的需求，迫切需要研發(fā)適合我國地形地質(zhì)條件的航空電磁系統(tǒng).

4 航空電磁正反演理論

4.1 航空電磁一維正演模擬

頻率域航空電磁法中常用的觀測裝置包括水平共面(HCP)和垂直共軸(VCA).根據(jù)Yin和Hodges(2007)及Yin 等(2014)，HCP裝置的磁場垂直分量為

(1)

而對于VCA裝置，x方向的水平磁場分量為

(2)

為計算時間域航空電磁響應(yīng)，通常采用將頻率域航空電磁響應(yīng)轉(zhuǎn)換到時間域的方法.根據(jù)牛之璉(2007)，利用傅里葉變換可以得到時間域響應(yīng)

(3)

式中H(ω)是頻率域響應(yīng).(3)式可轉(zhuǎn)化為貝塞爾函數(shù)積分，并可利用漢克爾變換進行計算.殷長春等(2013)對航空電磁一維正演進行了詳細的介紹，而殷長春等(2015)對航空電磁探測能力進行了研究.

4.2 航空電磁二、三維正演方法

航空電磁二、三維正演方法主要包括有限元法、有限差分法、積分方程法.有限差分法由于方法簡單易于實現(xiàn)，在電磁正演計算中應(yīng)用較為廣泛(Newman and Alumbaugh, 1995; Streich，2009；Liu and Yin，2013，2014), 但該方法目前只能采用規(guī)則矩形/六面體網(wǎng)格對計算區(qū)域進行剖分，無法精確擬合地形和復(fù)雜異常體邊界(Sasaki and Nakazato, 2003).有限元法既可使用規(guī)則矩形/六面體網(wǎng)格對計算區(qū)域進行剖分(Sugeng, 1998；王若等，2006)，又可使用不規(guī)則的三角形/四面體網(wǎng)格對計算區(qū)域進行剖分 (Key and Weiss, 2006; Ren et al., 2013)，因此可以更好地適用于復(fù)雜模型的模擬.隨著計算機技術(shù)的發(fā)展，有限單元法在電磁正演模擬中得到了廣泛的應(yīng)用，正逐漸成為一種高效、通用的正演模擬手段.積分方程法也是一種常用的電磁響應(yīng)模擬手段(Farquharson and Oldenburg,2002；Abubakar et al.,2006；Kamm and Pedersen,2014).然而，該方法最終形成的求解矩陣是密實的，占用內(nèi)存過大，只可用于簡單模型電磁響應(yīng)的正演模擬.航空電磁系統(tǒng)具有多源性特征，計算每一個測點的正演響應(yīng)都需要求解一次方程.傳統(tǒng)的迭代求解大型線性方程組方法對每個測點需要單獨求解，計算效率很低.多波前分解算法只需進行一次分解即可通過改變右端源項求解出所有測點的響應(yīng)，大幅度提高正演計算效率，能夠滿足航空電磁正演模擬的需求.

4.3 航空電磁數(shù)據(jù)處理

航空電磁在野外實際勘探時，由于受到大氣流、雷電、飛機自身振動及速度不穩(wěn)定等因素的影響，測量數(shù)據(jù)包含大量噪聲.而對數(shù)據(jù)成像和反演有意義的是強度較弱的二次場信號，很小的噪聲對反演結(jié)果會產(chǎn)生很大影響.此外，隨著儀器科學(xué)的發(fā)展，航空電磁系統(tǒng)實現(xiàn)了多通道密集采樣，這在提高航空電磁勘探精度的同時，也使航空電磁數(shù)據(jù)量變得非常龐大，數(shù)據(jù)解釋成本極大提高.因此，通過數(shù)據(jù)預(yù)處理及處理技術(shù)，提高信噪比、改善數(shù)據(jù)質(zhì)量、同時適當精簡數(shù)據(jù)量，對航空電磁成像與反演解釋有重要意義.通常航空電磁數(shù)據(jù)預(yù)處理包括天電和運動噪聲去除、疊加和抽道等，而數(shù)據(jù)處理包括姿態(tài)校正、濾波和調(diào)平等.

天電是指在大氣層中存在大量的帶電電荷，由于雷電作用，這些電荷通過放電產(chǎn)生電磁輻射，然后經(jīng)過空氣和地面進行傳播.天電噪聲具有頻率高、幅值大的特點，在測量數(shù)據(jù)上表現(xiàn)出數(shù)據(jù)發(fā)生突跳.天電去除可采用Alpha-trim濾波實現(xiàn).

運動噪聲是飛機在飛行過程中，由于機身振動、飛機速度不均勻、以及大氣氣流變化等原因引起發(fā)射與接收線圈振動而產(chǎn)生的噪聲.去除系統(tǒng)運動噪聲主要采用陷波法及拉格朗日最優(yōu)化算法.

航空電磁系統(tǒng)在飛行過程中數(shù)據(jù)采集非常密集.如果在后續(xù)數(shù)據(jù)處理時對全部數(shù)據(jù)進行處理會花費很多時間.通過疊加和抽道，可以大幅度降低數(shù)據(jù)量，同時提高數(shù)據(jù)的信噪比.疊加通常是通過對設(shè)定窗口內(nèi)數(shù)據(jù)加權(quán)平均獲得，而數(shù)據(jù)在疊加后可按一定的時間完成抽道.

航空電磁數(shù)據(jù)采集過程中，發(fā)射和接收線圈由于受到氣候環(huán)境和飛行條件等因素的影響會發(fā)生擺動(Roll)、俯仰(Pitch)和旋轉(zhuǎn)(Yaw).這種線圈姿態(tài)的變化會對采集數(shù)據(jù)產(chǎn)生強烈的影響 (Fitterman and Yin，2004).為此，人們在采集數(shù)據(jù)時通常也記錄航空電磁系統(tǒng)的姿態(tài).Yin和Fraser(2004b)對姿態(tài)效應(yīng)的研究發(fā)現(xiàn)：根據(jù)系統(tǒng)姿態(tài)特征可計算出姿態(tài)校正因子，姿態(tài)效應(yīng)校正量可利用觀測值與校正因子的比值求得.王琦等(2013)、曲昕馨等(2014)在Yin和Fraser(2004b)研究的基礎(chǔ)上，進一步研究了航空電磁系統(tǒng)姿態(tài)效應(yīng)和位移效應(yīng)的綜合校正算法.

航空數(shù)據(jù)調(diào)平的主要目的是消除由于物理環(huán)境變化造成的系統(tǒng)誤差.在成像圖中零漂偏移表現(xiàn)為塊狀或者條帶狀分布，通過調(diào)平可以校正這些數(shù)據(jù)誤差.常用的調(diào)平方法包括條件濾波法、自動調(diào)平法、二維移動平均濾波自動調(diào)平方法以及基于線與線的航空頻率域測量數(shù)據(jù)相關(guān)性的調(diào)平方法(Huang and Fraser，1999；李文杰，2007；Huang, 2008).4.4 航空電磁數(shù)據(jù)成像與反演解釋

航空電磁反演和成像可從航空電磁數(shù)據(jù)中獲得地下電性分布信息.航空電磁數(shù)據(jù)成像本質(zhì)上不是反演，它是把電磁觀測數(shù)據(jù)轉(zhuǎn)換為某些中間參數(shù)，如視電導(dǎo)率和視深度等.這些參數(shù)表征地下介質(zhì)主要電性分布特征.由于沒有數(shù)據(jù)擬合，無需傳統(tǒng)反演算法的最優(yōu)化過程.成像算法速度快，能從海量航空電磁數(shù)據(jù)中提取出地下主要電性信息，適用于現(xiàn)場快速數(shù)據(jù)處理；同時它也可為復(fù)雜的航空電磁反演提供初始模型.頻率域常用的成像算法包括Sengpiel Section (Sengpiel，1988)和微分電阻率法 (Huang and Fraser，1996)；而時間域成像算法包括查表法、 CDI/CDT、EMFlow (Wolfgram and Karlik, 1995; Macnae et al., 1998) 等.

航空電磁數(shù)據(jù)反演是采用非線性最優(yōu)化算法實現(xiàn)數(shù)據(jù)擬合.其計算速度比成像慢，但能夠提供更為準確的地電結(jié)構(gòu)信息.航空電磁數(shù)據(jù)反演是通過不斷改變地下模型參數(shù)來減小模型響應(yīng)與實測數(shù)據(jù)的擬合差，當擬合差足夠小時，得到的模型即可認為是反映地下介質(zhì)電性分布的最佳模型.航空電磁數(shù)據(jù)反演分為一維反演和多維反演，一維反演速度快、技術(shù)成熟，在實際生產(chǎn)中被廣泛應(yīng)用；多維反演雖然能夠提供更為精細的地下介質(zhì)分布信息，但由于計算速度慢等原因，目前還處于理論研究階段，在實際生產(chǎn)中應(yīng)用較少.從搜索方法上說，航空電磁數(shù)據(jù)反演分為下降搜索和全球最小搜索法.下降搜索法主要是通過假設(shè)初始模型，在下降方向搜索新模型，以實現(xiàn)理論模型和實測數(shù)據(jù)的最優(yōu)擬合.由于采用下降搜索，反演迭代過程和反演結(jié)果與初始模型關(guān)系密切.當初始模型選擇合適時(離真實模型較近)，則反演很快收斂，反之當初始模型離真實模型較遠時，反演迭代容易陷于局部極小值.下降搜索法包括Marquardt(范翠松，2013)方法，LCI(Auken et al.，2000；Siemon，2009；蔡晶等，2014)，非線性共軛梯度 (Rodi and Mackie, 2001; Kelbert et al., 2008)，高斯-牛頓 (Daniel, 1973; Jin, 2008)等方法.全球最小搜索法通常先假設(shè)初始模型，并按一定規(guī)則搜索新模型，當模型結(jié)果和數(shù)據(jù)的擬合差下降時，則該模型被無條件接受，反之則附加一定條件(比如概率)予以接受.因此，反演迭代過程既可下降搜索，又容許一定程度的上升搜索.如此，反演模型在迭代過程中可以跳出局部極小值，達到全球極小值.全球極小搜索法與初始模型關(guān)系很小.該方法包括模擬金屬淬火Simulated Annealing(Yin and Hodges，2007)，Bayesian等算法(殷長春等，2014).

目前處理航空電磁二、三維反演問題常用的方法包括高斯牛頓方法(GN)、擬牛頓方法(QN)和非線性共軛梯度方法(NLCG).GN算法通常利用預(yù)處理的共軛梯度算法求解反演方程，每次的共軛梯度迭代中需要兩次額外的伴隨正演來計算靈敏度矩陣與向量的乘積.L-BFGS(標準的QN方法)和NLCG每次迭代都只需要計算一次梯度，即只需一次正演和一次伴隨正演.與GN方法相比，L-BFGS和NLCG每次迭代需要的計算量小，但收斂速度較慢(劉云鶴和殷長春，2013).GN方法由于需要精確計算海森矩陣，且在進行伴隨正演過程中需要存儲每個源的全域解向量，因此不適合多源發(fā)射的電磁反演.相比之下，L-BFGS和NLCG方法可利用近似方法計算海森矩陣，在確定搜索方向后僅需計算搜索步長即可得到模型更新量，在航空電磁二、三維反演中得到了廣泛的應(yīng)用.

航空電磁數(shù)據(jù)具有多源性和緊湊性的特點，為加快反演速度，直接分解求解方程和footprint等技術(shù)被應(yīng)用于航空電磁反演中.直接分解求解方程算法在模型不變而僅右端源項變化的情況下能一次性求解多個源對應(yīng)的響應(yīng)，大幅度提高多源問題的正演和伴隨正演速度，進而提高反演效率.由于航空電磁系統(tǒng)的影響范圍(footprint)較小，我們在對航空電磁數(shù)據(jù)進行反演時可以將一個大的區(qū)域劃分成很多個小的區(qū)域分別進行反演，從而減小反演模型的規(guī)模，加快反演速度.

航空電磁數(shù)據(jù)反演應(yīng)遵循由簡單到復(fù)雜的過程，通過利用成像技術(shù)獲得地下電性主要分布特征，并確定初始模型進行一維反演獲得地下電性分布的大致信息，最后根據(jù)一維反演結(jié)果確定復(fù)雜模型用于航空電磁數(shù)據(jù)二、三維反演.Yin等(2015)對航空電磁數(shù)據(jù)反演進行了詳細的介紹和評述，并給出了大量的模型算例.

5 航空電磁勘查技術(shù)應(yīng)用

航空電磁勘查技術(shù)應(yīng)用十分廣泛，主要包括礦產(chǎn)資源、油氣、環(huán)境與工程、地下水及地?zé)豳Y源勘查等.近年來在考古、核廢料處理和垃圾填埋、海洋地形調(diào)查、極地研究及軍事等方面也逐漸獲得應(yīng)用. 航空電磁勘查技術(shù)分為頻率域和時間域，機載平臺包括直升機和固定翼系統(tǒng).從20世紀70年代末，頻率域固定翼航電系統(tǒng)被逐漸淘汰，目前常用的系統(tǒng)主要為時間域固定翼系統(tǒng)和頻率域/時間域直升機系統(tǒng).時間域固定翼航電系統(tǒng)工作效率高、探測深度大，適合尋找埋藏較深、體積較大的目標體，但其分辨率較直升機系統(tǒng)低.直升機航電系統(tǒng)分辨率高，但飛行效率低，適合尋找勘探區(qū)較小的目標體.

5.1 航空電磁在礦產(chǎn)資源勘查中的應(yīng)用

礦產(chǎn)資源勘查是航空電磁技術(shù)最為廣泛的應(yīng)用領(lǐng)域.早在1954年，人們使用航空電磁勘查技術(shù)找到了鋅-鉛-銀-銅多金屬礦.1974年，人們使用Mark VI INPUT時間域航空電磁系統(tǒng)在加拿大魁北克省找到了鋅-銅-銀礦床(Reed，1981).使用的發(fā)射波形為半正弦波，發(fā)射時間為1.05 ms，接收時間為2.42 ms，垂直發(fā)射線圈距離地表400 ft，接收線圈拖曳在發(fā)射線圈后320 ft，垂直距離為220 ft.在澳大利亞西部和坦桑尼亞人們使用航空電磁系統(tǒng)發(fā)現(xiàn)Harmony、Maggie Hays North和Kabanga等三處硫化鎳礦床(Wolfgramand Golden，2001). 其中，Maggie Hays North礦床勘查使用了三個系統(tǒng)(GEOTEM Deep、QUESTEM 450和TEMPEST)進行勘探，系統(tǒng)操作頻率均為25 Hz.QUESTEM數(shù)據(jù)由于氣候原因?qū)е略肼曔^大未能參與解釋，但GEOTEM Deep和TEMPEST系統(tǒng)采集的數(shù)據(jù)給出的結(jié)果基本一致 (Petersand Buck，2000).Harmony和Kabanga礦床的勘查均使用了GEOTEM系統(tǒng)，但Harmony使用的是25 Hz和4 ms發(fā)射信號，而Kabanga使用的是75 Hz和1 ms發(fā)射信號.除以上應(yīng)用實例，航空電磁勘查技術(shù)在尋找金礦、鈾礦及金伯利巖(金剛石)等方面都獲得廣泛的應(yīng)用.2000年，F(xiàn)ugro公司使用MEGATEM系統(tǒng)對Matagami地區(qū)進行了飛行觀測.這次觀測直接發(fā)現(xiàn)了Perseverance礦床，并于一個月后在該礦床附近發(fā)現(xiàn)了Perseverance West和Equinox兩處新礦床.2005年，Geotech公司在Zambia Lusaka 銅金礦利用VTEM系統(tǒng)進行航空電磁勘查，數(shù)據(jù)解釋將目標礦體沿走向方向延長了4倍.2008年，Geotech公司在加拿大安省Northern Empire金礦進行時間域航空電磁勘查，揭示出位于金礦下部存在的良導(dǎo)體.2007和2010年Geotech公司利用VTEM和ZTEM系統(tǒng)在加拿大安省East Bull Lake礦區(qū)進行航空電磁觀測，其中VTEM成功查明位于兩個成礦帶交叉部位的目標體，并經(jīng)鉆探驗證；而ZTEM觀測數(shù)據(jù)解釋揭示了礦區(qū)800 m深的構(gòu)造特征.2012年，CGG公司使用HeliTEM在加拿大的Lalor湖對火山成因的硫化物礦床進行了勘探(Yang and Oldenburg，2013).雖然本次勘探目標深度超出了傳統(tǒng)航空電磁系統(tǒng)勘探范圍，但HeliTEM由于勘探深度大很好地確定了礦體的異常位置.5.2 航空電磁在油氣勘查中的應(yīng)用

地震方法通過圈定構(gòu)造實現(xiàn)油氣勘查.然而，地震方法難以確定構(gòu)造中的含油、含水特征，容易導(dǎo)致油氣開發(fā)中的干井問題，浪費資金.航空電磁勘探數(shù)據(jù)可以獲得地下介質(zhì)電阻率信息，在確定構(gòu)造前提下能夠進一步判斷構(gòu)造中是否儲藏油氣，降低干井率.2008年，Aeroquest公司使用AeroTEM IV平臺在莫桑比克進行了天然氣勘探(Pfaffhuber et al.，2009)，發(fā)射波形為75 Hz的三角波，通電時間和斷電時間分別為2.06和6 ms，發(fā)射線圈直徑為12 m，共五匝，接收線圈分別接收水平和垂直磁感應(yīng)分量.該航空電磁勘查項目為解譯測區(qū)構(gòu)造，識別天然氣運移通道提供了重要依據(jù).2004年，F(xiàn)ugro公司使用RESOLVE 直升機系統(tǒng)在美國蒙大拿州東北部的East Poplar油田進行了數(shù)據(jù)采集和成像(Smith et al.，2006).該系統(tǒng)的吊艙內(nèi)安置了六對線圈(五對為水平共面，一對為直立共軸)，發(fā)射接收頻率范圍為400～14 kHz，傳感器距離地面高度為30 m.該項目數(shù)據(jù)解釋結(jié)果為East Poplar油田區(qū)域地下鹽水分布的精細刻畫提供了有利技術(shù)參數(shù).

5.3 航空電磁在環(huán)境工程領(lǐng)域的應(yīng)用

航空電磁在環(huán)境工程方面的應(yīng)用主要包括壩基勘測、永凍層和極地冰層研究、地下管網(wǎng)調(diào)查、山體滑坡、海侵及淺海海底地形調(diào)查等.1999年，奧地利Vorarberg州發(fā)生嚴重的山體滑坡后，管理部門使用Geotech系統(tǒng)對該地區(qū)進行了地球物理勘查(Supper et al., 2008).該系統(tǒng)采用發(fā)射頻率為3200 Hz(共軸裝置)和7190 Hz(共面裝置)，系統(tǒng)同時還配有一個激光測高和兩個GPS裝置.此次勘查為該地區(qū)的山體滑坡風(fēng)險評估和制定有效的防治措施提供了重要的依據(jù).2003年，人們對Finland灣和Bothnia灣進行了海冰厚度探測(Haas，2004)，使用的裝置為Alfred Wegener研究所開發(fā)的直升機吊艙系統(tǒng).吊艙長3.5 m，重100 kg，懸掛在直升機下20 m，距離海面約15 m，系統(tǒng)采用3.6 kHz和112 kHz兩個操作頻率.數(shù)據(jù)處理結(jié)果表明，本次海冰厚度測量具有較高精度，測量的海冰厚度比已有資料給出的海冰厚度大很多，對資料更新和分析海冰變化情況具有重要的意義.2006年，SkyTEM公司在Jutland, Denmark西南部海岸進行了勘查，其目標是了解該地區(qū)海水污染對民眾生活的影響.該公司使用的SkyTEM系統(tǒng)發(fā)射線圈面積314 m2，最大發(fā)射磁矩120000 Am2，同時接收X和Z分量數(shù)據(jù).此次勘查通過反演得到了測區(qū)的電阻率剖面，揭示了該區(qū)的海水及地下水的運移特征.2008年至2009年，德國萊布尼茲研究所和聯(lián)邦地球科學(xué)與自然資源研究所(BGR)聯(lián)合對德國北部海岸的沿海含水地層和海侵情況進行調(diào)查(Wiederhold et al.，2010).本次調(diào)查采用BGR的頻率域系統(tǒng)RESOLVE和SkyTEM Aps的時間域系統(tǒng)SkyTEM，勘查結(jié)果揭示了這一區(qū)域淡水和海水的分布特征信息，為了解該地區(qū)海侵情況提供了依據(jù).2010年，F(xiàn)ugro公司使用RESOLVE系統(tǒng)對阿拉斯加附近Yukon平原的永凍層進行了勘查(Minsley et al.，2012).系統(tǒng)的操作頻率400～129 kHz，飛行的平均速度為30 m·s-1，飛行高度為30 m，測線長度共計900 km，覆蓋面積達300 km2.勘探數(shù)據(jù)及處理結(jié)果不僅提供了該地區(qū)永凍層的詳細信息，還揭示了該測區(qū)地下水分布信息.

5.4 航空電磁在地下水資源勘查中的應(yīng)用

地下水通常埋藏較淺，且其電阻率與圍巖有顯著差異.航空電磁法由于對淺層目標體具有很高的分辨率，是地下水資源勘查的理想選擇之一.目前用于地下水勘查的航空電磁系統(tǒng)主要有Fugro的頻域直升機RESOLVE系統(tǒng)、SkyTEM公司的時間域直升機系統(tǒng)等.1998年，SkyTEM Aps使用時間域電磁系統(tǒng)在馬來西亞Selangor州進行了地下水勘探.發(fā)射線圈面積為341 m2，最大發(fā)射偶極矩為150000 Am2，接收線圈接收X和Z分量數(shù)據(jù).此次勘查目的是在該地區(qū)尋找水源并確定打井位置，并為Kuala Lumpur、Putrajaya 和Selangor三個地區(qū)尋找可持續(xù)的水源.2001年，F(xiàn)ugro公司在美國Texas州西部的EI Paso和Sterling地區(qū)分別進行了162和372 km2的地下水鹽化度勘查(Paine and Collins，2003).本次勘查使用的裝置為MEGATEM II系統(tǒng)，飛行高度為120 m，三分量電磁接收器拖曳在發(fā)射線圈后131 m處，距離地面高度75 m，發(fā)射線圈發(fā)射頻率為30 Hz，發(fā)射電流為1330 A的正弦波，發(fā)射偶極矩為2100000 Am2. 此次勘探通過將航空電磁數(shù)據(jù)與地面電磁和井數(shù)據(jù)進行綜合分析，摸清Sterling地區(qū)地下水鹽堿化情況，并在EI Paso找到了良好的淡水水源.2005年，F(xiàn)ugro公司使用RESOLVE頻率域航空電磁系統(tǒng)對Chowilla河漫灘進行了地下水礦化調(diào)查(Tan et al.， 2009).吊艙距離地表30 m，使用五對水平共面(操作頻率為390～132 kHz)和一對直立共軸線圈對(操作頻率為3242 Hz).對該地區(qū)采集數(shù)據(jù)進行電導(dǎo)率成像和反演，揭示地下水鹽化情況和分布信息.2009年，SkyTEM在意大利Venice瀉湖對地表水、地下水轉(zhuǎn)換情況進行了調(diào)查(Viezzoli et al.，2010).使用裝置的發(fā)射偶極矩為200000 Am2，基頻為12.5 Hz，吊艙距離地表30～40 m.這次勘查在陸地和水交界位置采集的航電數(shù)據(jù)為確定地下水鹽化度提供了寶貴的信息，并為了解地表水和地下水之間的轉(zhuǎn)換關(guān)系提供了技術(shù)參數(shù).5.5 航空電磁在其他領(lǐng)域的應(yīng)用

航空電磁法因其工作效率高、對惡劣工作條件適應(yīng)性強，因此獲得更廣泛應(yīng)用.其領(lǐng)域已拓寬到農(nóng)業(yè)普查、核電站選址、尋找沉船、古墓探測等.1980年，航空電磁勘探在法國被應(yīng)用于核電站選址(Deletie and Lakshmanan，1986).他們使用的直升機系統(tǒng)的吊艙長9 m，內(nèi)部含有三個線圈對(垂直共面、水平共面、垂直共軸)，工作時吊艙距離地面高度30 m，工作頻率為375、900、3600和8000 Hz.這次勘查最終在法國找到了7處可供選擇的核電站建造地址.2001年，中國地質(zhì)科學(xué)院地球物理地球化學(xué)研究所研制的HDY-402三頻航空電磁系統(tǒng)在吉林省乾安地區(qū)開展水資源普查，對該地區(qū)農(nóng)業(yè)生態(tài)地質(zhì)進行評價.該系統(tǒng)收發(fā)裝置采用垂直共面線圈對，收發(fā)距為19.2 m，工作頻率分別為463、1563和8333 Hz，最高采樣頻率為8次/秒，噪聲水平小于20 ppm.本次調(diào)查查明了吉林省乾安地區(qū)土壤鹽漬化，為該區(qū)農(nóng)業(yè)生態(tài)的發(fā)展提供技術(shù)參數(shù).

6 結(jié)論與展望

我國現(xiàn)處于經(jīng)濟轉(zhuǎn)型和快速發(fā)展時期，對礦產(chǎn)資源的需求已經(jīng)成為制約國家經(jīng)濟發(fā)展的重要因素之一.國家近年不斷加大對礦產(chǎn)資源勘探開發(fā)力度，然而，我國地質(zhì)條件良好地區(qū)的礦產(chǎn)資源勘查工作已基本完成，勘查目標已轉(zhuǎn)移到地球深部和廣大西部無人區(qū).航空電磁勘查技術(shù)由于采用飛行平臺無需地面人員接近，勘查效率高，成本較低，特別適合地形地質(zhì)條件復(fù)雜地區(qū)資源勘查.其應(yīng)用領(lǐng)域涉及到礦產(chǎn)、環(huán)境工程、地下水及地?zé)豳Y源勘查等.因此，大力發(fā)展適應(yīng)于我國復(fù)雜地形地質(zhì)條件的航空電磁系統(tǒng)迫在眉睫.

未來我國航空電磁勘查技術(shù)發(fā)展面臨的技術(shù)難點主要包括：1)大功率多脈沖發(fā)射技術(shù)； 2)多分量大動態(tài)范圍電磁接收技術(shù)；3)on-和off-time電磁信號解釋技術(shù)，提高淺部地表和深部目標體探測能力；4)基于物理成因的航空電磁數(shù)據(jù)處理技術(shù)；5)快速三維電磁正、反演技術(shù).通過系統(tǒng)總結(jié)國內(nèi)外航空電磁勘查技術(shù)發(fā)展現(xiàn)狀，借鑒國外技術(shù)研究和系統(tǒng)開發(fā)的成功經(jīng)驗，以期突破國內(nèi)航空電磁發(fā)展的技術(shù)瓶頸，發(fā)展適合我國特殊地質(zhì)條件的航空電磁勘查技術(shù)和儀器系統(tǒng).

鑒于時間域和頻率域航空電磁系統(tǒng)各具特色，勘探目標和領(lǐng)域各不相同.建議在發(fā)展我國航空電磁勘查技術(shù)過程中同步研發(fā)時間域和頻率域電磁系統(tǒng).半航空系統(tǒng)由于采用了遠距離接地發(fā)射，接收信號存在嚴重體積效應(yīng)，喪失了航空電磁勘查技術(shù)的高分辨率優(yōu)勢，勘查精度無法與傳統(tǒng)全航空電磁技術(shù)相比，在實際礦產(chǎn)勘查中無法取代全航空系統(tǒng).無人機航電系統(tǒng)由于受到飛機本身特性的制約難以承擔大規(guī)模的勘探任務(wù).在發(fā)射裝置方面由于無人機系統(tǒng)承載重量小，導(dǎo)致其發(fā)射功率和勘探深度都十分有限；在數(shù)據(jù)采集方面，無人機飛行平臺穩(wěn)定性較差，采集的數(shù)據(jù)難以滿足精細勘查的要求；而在實際操作方面，無人機飛行平臺氣候適應(yīng)能力差，無法應(yīng)用于地形條件復(fù)雜的地區(qū)，同時續(xù)航時間有限.無人機和半航空電磁勘查系統(tǒng)由于其自身固有的弱點，難以擔當我國礦產(chǎn)資源勘查的主力軍，發(fā)展全航空電磁勘查系統(tǒng)應(yīng)該成為我國航空電磁勘查技術(shù)發(fā)展的主流方向.

Abubakar A, van den Berg P M, Habashy T M. 2006. An integral equation approach for 2.5-dimensional forward and inverse electromagnetic scattering.Geophys.J.Int., 165(3): 744-762.Annan A P, Lockwood R. 1991. An application of airborne GEOTEM in Australian conditions.ExplorationGeophysics, 22(1): 5-12.

Auken E, Sφrensen K I, Thomsen P. 2000. Lateral constrained inversion (LCI) of profile oriented data—the resistivity case. ∥ Proceedings of the EEGS-ES, Bochum, Germany, EL06.

Balch S L, Boyko W P, Paterson N R. 2003. The AeroTEM airborne electromagnetic system.TheLeadingEdge, 22: 562-566.Beamish D. 2003. Airborne EM footprints.GeophysicalProspecting, 51(1): 49-60.Cai J, Qi Y F, Yin C C. 2014. Weighted Laterally-constrained inversion of frequency-domain airborne EM data.ChineseJ.Geophys. (in Chinese), 57(3): 953-960, doi: 10.6038/cjg20140324.

Chen S D, Lin J, Zhang S. 2012. Effect of transmitter current waveform on TEM response.ChineseJ.Geophys. (in Chinese), 55(2): 709-716, doi: 10.6038/j.issn.0001-5733.2012.02.035.Chen T, Hodges G, Christensen A N, et al. 2014. Multipulse Airborne TEM Technology and Test Results Over Oil-sands. ∥ 76th EAGE Conference and Exhibition-Workshops.

Daniel J W. 1973. Newton′s method for nonlinear inequalities.NumerischeMathematik, 21(5): 381-387.

Deletie E, Lakshmanan J. 1986. Airborne resistivity surveying applied to nuclear power plant site investigation in France. ∥ Palacky G J ed. Airborne Resistivity Mapping: Geol. Surv. Canada, Paper 86-22: 145-152.

Farquharson C G, Oldenburg D W. 2002. An integral equation solution to the geophysical electromagnetic forward-modelling problem.MethodinGeochemistryandGeophysics, 35: 3-19.

Fan C S. 2013. Research on complex resistivity forward and inversion with finite element method and its application [Ph. D. thesis] (in Chinese). Changchun: Jilin University.

Fitterman D V, Yin C C. 2004. Effect of bird maneuver on frequency-domain helicopter EM response.Geophysics, 69(5): 1203-1215.Fountain D. 1998. Airborne electromagnetic systems-50 years of development.ExplorationGeophysics, 29(2): 1-11.

Fountain D, Smith R, Payne T, et al. 2005. A helicopter time-domain EM system applied to mineral exploration: system and data.FirstBreak, 23: 73-78.

Haas C. 2004. Airborne EM sea-ice thickness profiling over brackish Baltic sea water. ∥ 17th International Symposium on Ice Saint Petersburg, Russia. International Association of Hydraulic Engineering and Research, 12-17.

Huang H P, Fraser D C. 1996. The differential parameter method for multifrequency airborne resistivity mapping.Geophysics, 61(1): 100-109.

Huang H P, Fraser D C. 1999. Airborne resistivity data leveling.Geophysics, 64(2): 378-385.

Huang H P, Fraser D C. 2001. Mapping of the resistivity, susceptibility, and permittivity of the earth using a helicopter-borne electromagnetic system.Geophysics, 66(1): 148-156.

Huang H P. 2008. Airborne geophysical data leveling based on line to line correlations.Geophysics, 73(3): F83-F89.

Jin Q N. 2008. A convergence analysis of the iteratively regularized Gauss-Newton method under the Lipschitz condition.InverseProblem, 24: 045002, doi: 10.1088/0266-5611/24/4/045002.Kaminski V, Oldenburg D. 2012. The geophysical study of Drybones kimberlite using 3D Time Domain EM inversion and 3D ZTEM inversion algorithms. ∥ 22nd International Geophysical Conference and Exhibition, 26-29.Kamm J, Pedersen L B. 2014. Inversion of airborne tensor VLF data using integral equations.GeophysicalJournalInternational, 198(2): 775-794.Kelbert A, Egbert G D, Schultz A. 2008. Non-linear conjugate gradient inversion for global EM induction: resolution studies.Geophys.J.Int., 173(2): 365-381.

Key K, Weiss C. 2006. Adaptive finite-element modeling using unstructured grids: The 2D magnetotelluric example.Geophysics, 71(6): G291-G299.Kovacs A, Holladay J S, Bergeron C J J. 1995. The footprint/altitude ratio for helicopter electromagnetic sounding of sea ice thickness: comparison of theoretical and field estimates.Geophysics, 60(2): 374-380.

Li W J. 2007. An application of 2-D data leveling technique for frequency domain airborne electromagnetic survey.JournalofChengduUniversityofTechnology(ScienceandTechnologyEdition) (in Chinese), 34(4): 447-451.

Liu G M, Becker A. 1990. Two-dimensional mapping of sea-ice keels with airborne electromagnetics.Geophysics, 55(2): 239-248.

Liu G. 1998. Effect of transmitter current waveform on airborne TEM response.ExplorationGeophysics, 29(2): 35-41.

Liu Y H, Yin C C. 2013. 3D inversion for frequency-domain HEM data.ChineseJ.Geophys. (in Chinese), 56(12): 4278-4287, doi: 10.6038/cjg20131230.

Liu Y H, Yin C C. 2014. 3D anisotropic modeling for airborne EM systems using finite-difference method.JournalofAppliedGeophysics, 109: 186-194.

Macnae J, King A, Stolz N, et al. 1998. Fast AEM data processing and inversion.ExplorationGeophysics, 29(2): 163-169.

Minsley B J, Abraham J D, Smith B D, et al. 2012. Airborne electromagnetic imaging of discontinuous permafrost.GeophysicalResearchLetters, 39: L02503, doi: 10.1029/2011GL050079.Newman G A, Alumbaugh D L. 1995. Frequency-domain modeling of airborne electromagnetic responses using staggered finite differences.GeophysicalProspecting, 43(8): 1021-1042.

Niu Z L. 2007. Time Domain Electromagnetic Method Principle (in Chinese). Changsha: Central South University Press.

Paine J G, Collins E W. 2003. Applying airborne electromagnetic induction in groundwater salinization and resource studies, West Texas. ∥ Symposium on the Application of Geophysics to Engineering and Environmental Problems, 722-738.

Peters B, Buck P. 2000. The Maggie Hays and Emily Ann nickel deposits, Western Australia: A geophysical case history.ExplorationGeophysics, 31(2): 210-221.Pfaffhuber A A, Monstad S, Rudd J. 2009. Airborne electromagnetic hydrocarbon mapping in Mozambique.ExplorationGeophysics, 40(3): 237-245.Qi Z P, Zhi Q Q, Li X, et al. 2014. The definition of the full-zone apparent resistivity and the constrained inversion of the three components of fixed source TEM.GeophysicalandGeochemicalExploration(in Chinese), 38(4): 742-749.Qiang J K, Luo Y Z, Tang J T, et al. 2010. The algorithm of all-time apparent resistive for Airborne Transient Electromagnetic (ATEM) survey.ProgressinGeophysics(in Chinese), 25(5): 1657-1661, doi: 10.3969.issn.10042903.2010.05.018.

Qu X X, Li T L, Wang F. 2014. A study on the effect of bird′s attitude on helicopter EM data and the correction method.ChineseJ.Geophys. (in Chinese), 57(4): 1310-1321, doi: 10.6038/cjg20140428.

Reed L E. 1981. The airborne electromagnetic discovery of the Detour zinc-copper-silver deposit, northwestern Quebec.Geophysics, 46(9): 1278-1290. Reid J E, Vrbancich J. 2004. A comparison of the inductive limit footprints of airborne electromagnetic configurations.Geophysics, 69(5): 1229-1239.Reid J E, Pfaffling A, Vrbancich J. 2006. Airborne electromagnetic footprints in 1D earths.Geophysics, 71(2): G63-G72.

Ren Z, Kalscheuer T, Greenhalgh S, et al. 2013. A goal-oriented adaptive finite-element approach for plane wave 3-D electromagnetic modeling.GeophysicalJournalInternational, 194(2): 700-718.Rodi W, Mackie R L. 2001. Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion.Geophysics, 66(1): 174-187.Sasaki Y, Nakazato H. 2003. Topographic effects in frequency-domain helicopter-borne electromagnetic.ExplorationGeophysics, 34(2): 24-28.Sengpiel K P. 1988. Approximate inversion of airborne EM data form a multilayered ground.GeophysicalProspecting, 36(4): 446-459.

Siemon B, Auken E, Christiansen A V. 2009. Laterally constrained inversion of helicopter-borne frequency-domain electromagnetic data.JournalofAppliedGeophysics, 67(3): 259-268.

Slichter L B. 1955. Geophysics applied to prospecting for ores: Economic Geology. ∥ 50th Anniversary Volume.

Smith R, Fountain D, Allard M. 2003. The MEGATEM fixed-wing transient EM system applied to mineral exploration: a discovery case history.FirstBreak, 21(7): 73-77.

Smith B D, Thamke J N, Cain J M, et al. 2006. Helicopter Electromagnetic and Magnetic Survey Maps and Data, East Poplar oil field area, Fort Peck Indian Reservation, northeastern Montana, August 2004.Streich R. 2009. 3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: Direct solution and optimization for high accuracy.Geophysics, 74(5): F95-F105.

Sugeng F. 1998. Modelling the 3D TDEM response using the 3D full-domain finite-element method based on the hexahedral edge-element technique.ExplorationGeophysics, 29(4): 615-619.

Supper R, R?mer A, Jochum B, et al. 2008. A complex geo-scientific strategy for landslide hazard mitigation-from airborne mapping to ground monitoring.AdvancesinGeosciences, 14: 195-200.

Tan K P, Munday T, Halas L, et al. 2009. Utilizing airborne electromagnetic data to map groundwater salinity and salt store at Chowilla, SA. ∥ ASEG Extended Abstracts 2009, 20th Geophysical Conference.

Thomson S, Fountain D, Watts T. 2007. Airborne Geophysics-Evolution and Revolution. ∥ Proceedings of Exploration 07, 5th Decennial International Conference on Mineral Exploration, 19-37.Viezzoli A, Tosi L, Teatini P, et al. 2010. Surface water-groundwater exchange in transitional coastal environments by airborne electromagnetics: The Venice Lagoon example.GeophysicalResearchLetters, 37: L01402, doi: 10.1029/2009GL041572.Wang Q, Lin J, Yu S B, et al. 2013. Study on influence and correction of coil attitude and bird swing for the fix-wing time-domain electromagnetic system.ChineseJ.Geophys. (in Chinese), 56(11): 3741-3750, doi: 10.6038/cjg20131115.Wang R, Wang M Y, Di Q Y. 2006. Electromagnetic modeling due to line source in frequency domain using finite element method.ChineseJ.Geophys. (in Chinese), 49(6): 1858-1866, doi: 10.3321/j.issn:0001-5733.2006.06.035.

Wiederhold H, Siemon B, Steuer A, et al. 2010. Coastal aquifers and saltwater intrusions in focus of airborne electromagnetic surveys in Northern Germany. ∥ 21st Salt Water Intrusion Meeting Azores (Portugal), June 2010.

Witherly K, Irvine R, Morrison E B. 2004. The Geotech VTEM Time Domain Helicopter EM System. ∥ SEG Int′l Exposition and 74th Annual Meeting *Denver, Colorado* 10-15 October 2004.

Witherly K, Sattel D. 2012. The application of ZTEM to porphyry copper-gold exploration. ∥ 22nd International Geophysical Conference and Exhibition, 26-29.

Witherly K, Sattel D. 2013. An assessment of ZTEM and time domain EM results over three mineral deposits. ∥ 23rd International Geophysical Conference and Exhibition, Melbourne, 1-4.Wolfgram P, Karlik G. 1995. Conductivity-depth transform of GEOTEM data.ExplorationGeophysics, 26(2-3): 179-185.

Wolfgram P, Golden H. 2001. Airborne EM applied to sulphide nickel-examples and analysis.ExplorationGeophysics, 32: 136-140.

Yang D K, Oldenburg D W. 2013. 3D conductivity models of Lalor Lake VMS deposit from ground loop and airborne EM data sets. ∥ 23rd International Geophysical Conference and Exhibition, Melbourne, Australia, 1-4.Yin C C, Piao H R. 1991. A study of the definition apparent resistivity in electromagnetic sounding.GeophysicalandGeochemicalExploration(in Chinese), 15(4): 290-299.

Yin C C, Fraser D C. 2004a. The effect of the electrical anisotropy on the response of helicopter-borne frequency-domain electromagnetic systems.GeophysicalProspecting, 52(5): 399-416.

Yin C C, Fraser D C. 2004b. Attitude corrections of helicopter EM data using a superposed dipole model.Geophysics, 69(2): 431-439.

Yin C C, Hodges G. 2007. Simulated annealing for airborne EM inversion.Geophysics, 72(4): F189-F195.

Yin C C, Huang W, Ben F. 2013. The full-time electromagnetic modeling for time-domain airborne electromagnetic systems.ChineseJ.Geophys. (in Chinese), 56(9): 3153-3162, doi: 10.6038/cjg20130928.

Yin C C, Huang X, Liu Y H, et al. 2014. Footprint for frequency-domain airborne electromagnetic systems.Geophysics, 79(6): E243-254.

Yin C C, Qi Y F, Liu Y H, et al. 2014. Trans-dimensional Bayesian inversion of frequency-domain airborne EN data.ChineseJ.Geophys. (in Chinese), 57(9): 2971-2980.

Yin C C, Ren X Y, Liu Y H, et al. 2015. Exploration capability of airborne TEM systems for typical targets.ChineseJ.Geophys. (in Chinese)， in Press.Yin C C, Ren X Y, Liu Y H, et al. 2015. Review on airborne electromagnetic inverse theory and applications.Geophysics, 80(4): W17-W31.

Zhang Y Y. 2013. Study on multi-component interpretation and full field apparent resistivity definition of semi-airborne transient electromagnetic method with electric dipole on the surface [Master′s thesis] (in Chinese). Xi′an: Chang′an University.

附中文參考文獻

蔡晶, 齊彥福, 殷長春. 2014. 頻率域航空電磁數(shù)據(jù)的加權(quán)橫向約束反演. 地球物理學(xué)報, 57(1): 953-960, doi: 10.6038/cjg20140324.

陳曙東, 林君, 張爽. 2012. 發(fā)射電流波形對瞬變電磁響應(yīng)的影響. 地球物理學(xué)報, 55(2): 709-716, doi: 10.6038/j.issn.0001-5733.2012.02.035.

范翠松. 2013. 基于有限元法的復(fù)電阻率正反演研究及應(yīng)用[博士論文]. 長春: 吉林大學(xué).

李文杰. 2007. 用于頻率域航空電磁數(shù)據(jù)的二維自動調(diào)平. 成都理工大學(xué)學(xué)報(自然科學(xué)版), 34(4): 447-451.

劉云鶴, 殷長春. 2013. 三維頻率域航空電磁反演研究. 地球物理學(xué)報, 56(12): 4278-4287, doi: 10.6038/cjg20131230.

牛之璉. 2007. 時間域電磁法原理. 長沙: 中南大學(xué)出版社.

戚志鵬, 智慶全, 李貅等. 2014. 大定源瞬變電磁三分量全域視電阻率定義與三分量聯(lián)合反演. 物探與化探, 38(4): 742-749.

強建科, 羅延鐘, 湯井田等. 2013. 航空瞬變電磁法的全時域視電阻率計算方法. 地球物理學(xué)進展, 25(5): 1657-1661, doi: 10.3969.issn.10042903.2010.05.018.

曲昕馨, 李桐林, 王飛. 2014. 直升機吊艙姿態(tài)變化對電磁場測量的影響規(guī)律及其校正方法研究. 地球物理學(xué)報, 57(4): 1310-1321, doi: 10.6038/cjg20140428.

王琦, 林軍, 于生寶等. 2013. 固定翼航空電磁系統(tǒng)的線圈姿態(tài)及吊艙擺動影響研究與校正. 地球物理學(xué)報, 56(11): 3741-3750, doi: 10.6038/cjg20131115.

王若, 王妙月, 底青云. 2006. 頻率域線源大地電磁法有限元正演模擬. 地球物理學(xué)報, 49(6): 1858-1866, doi: 10.3321/j.issn:0001-5733.2006.06.035.

殷長春, 樸華榮. 1991. 電磁測深法視電阻率定義問題的研究. 物探與化探, 15(4): 290-299.

殷長春, 黃威, 賁放. 2013. 時間域航空電磁系統(tǒng)瞬變?nèi)珪r響應(yīng)正演模擬. 地球物理學(xué)報, 56(9): 3153-3162, doi: 10.6038/cjg20130928.

殷長春, 齊彥福, 劉云鶴等. 2014. 頻率域航空電磁數(shù)據(jù)變維數(shù)貝葉斯反演研究. 地球物理學(xué)報, 57(9): 2971-2980.

殷長春, 任秀艷, 劉云鶴等. 2015. 時間域航空電磁法對深部良導(dǎo)體探測能力研究. 地球物理學(xué)報, 待發(fā)表.

張瑩瑩. 2013. 水平電偶源地空系統(tǒng)瞬變電磁法多分量解釋技術(shù)及全域視電阻率定義研究[碩士論文]. 西安: 長安大學(xué).

(本文編輯 何燕)

Review on airborne EM technology and developments

YIN Chang-Chun, ZHANG Bo*, LIU Yun-He, REN Xiu-Yan, QI Yan-Fu, PEI Yi-Feng, QIU Chang-Kai, HUANG Xin, HUANG Wei, MIAO Jia-Jia， CAI Jing

CollegeofGeo-explorationSciencesandTechnology,JilinUniversity,Changchun130026,China

As an effective and efficient geophysical tool, airborne EM (AEM) is specifically suitable for the exploration in areas of high mountains, desert, swamp, and forest. With the development of national economy of China, the demand for mineral resources increases sharply, geophysical explorations in areas with favorable geological conditions have been accomplished. The exploration target is switched to areas with unfavorable and complicated geological conditions, such as in Western China. Airborne EM works efficiently due to its moving platform of helicopter or fixed-wing aircraft, no human access is needed to the survey area. This technology is especially applicable for exploration in Western China with rugged mountains and deserts. In this paper, we review AEM technology systematically with the goal to make this technology quickly applicable in the mineral exploration in China.To make a comprehensive review, we introduce the AEM technology in following sequences. We first present basic AEM theory of on/off time, apparent resistivity and depth, footprint, etc. Then, we introduce the developments of AEM technology in Western countries, including CGG/Fugro, Geotech and Aeroquest in Canada, SkyTEM survey in Denmark, and we address specifically the status of this technology in China. Since the modeling and inversions are fundamental for AEM, we present 1D/2D/3D forward theory based on semi-analytical solutions, finite element, finite difference and integral equation methods. For the inversion theory, we follow the rule from simple to complex models by presenting 1D Marquardt, LCI, and 2/3D inversions by GS, NLCG and QN. Finally, we introduce successful applications of AEM technology in mineral, oil & gas, E&E, ground water, and natural hazard forest. We give our suggestions on future development of AEM technology in China.AEM technologies have been well developed in Western countries and are playing very important role in mineral, oil & gas, E&E, ground water explorations. However, AEM has not been well developed in China to the extent of practical use due to the fact that no breakthrough has been made on key technologies. Due to the fact that time-domain and frequency-domain AEM have different features and applications, we suggest to develop AEM technology in both frequency- and time-domain in China. The frequency-domain AEM is used in E&E, ground water exploration, while time-domain AEM is used in mineral, oil & gas in deep earth. Semi-airborne EM has serious problems of volume effect and thus is not recommended. Unmanned aerial vehicle (UAV) has a small load, making it difficult to transmit high-power signal. Besides, UAV doesn′t suit the complex climate and topography, so that it is not recommended, either.To explore the vast area in Western China with rugged mountains, desert, etc., we need to develop geophysical technologies based on moving platforms, like airborne EM. Considering that different AEM systems suit for different exploration targets, both frequency- and time-domain systems need to be developed. Both semi-airborne and UAV are not recommended for airborne EM due to their inherent defects, a full airborne EM technology is the best solution.

Mineral exploration; Airborne EM (AEM); Frequency-domain and time-domain AEM systems; AEM data processing; Imaging and data interpretation

國家自然科學(xué)基金項目(41274121)和國家青年基金項目(41404093)、國家重大科研裝備研究項目(ZDYZ2012-1-03和20130523MTEM05)、吉林大學(xué)研究生創(chuàng)新基金資助項目(2015053)聯(lián)合資助.

殷長春，男，1965年生，教授，國家“千人計劃”特聘專家，主要從事電磁勘探理論，特別是航空和海洋電磁方面的研究. E-mail: yinchangchun@jlu.edu.cn

*通訊作者張博，男，1989年生，博士生，主要從事航空電磁正反演研究. E-mail: em_zhangbo@163.com

10.6038/cjg20150804.

10.6038/cjg20150804

P631

2014-12-14，2015-06-29收修定稿

殷長春, 張博, 劉云鶴等.2015.航空電磁勘查技術(shù)發(fā)展現(xiàn)狀及展望.地球物理學(xué)報，58(8):2637-2653,

Yin C C, Zhang B, Liu Y H,et al. 2015. Review on airborne EM technology and developments.ChineseJ.Geophys. (in Chinese),58(8):2637-2653,doi:10.6038/cjg20150804.