房江奇， 楊金政， 蔡文軍， 祝亞荔， 倪衛(wèi)沖

(1.核工業(yè)航測(cè)遙感中心， 石家莊 050002；2.中核集團(tuán)鈾資源地球物理勘查技術(shù)中心(重點(diǎn)實(shí)驗(yàn)室)， 石家莊 050002)

0 引言

航空磁法測(cè)量是將磁法探測(cè)儀器安裝在固定翼飛機(jī)、直升機(jī)或其他飛行器上，按一定高度沿測(cè)線飛行并記錄磁場(chǎng)變化，以解決地質(zhì)、地球物理問(wèn)題為目的的測(cè)量工作。該方法具有測(cè)量速度快、效率高，成本低，不受地域限制等優(yōu)勢(shì)，廣泛應(yīng)用于礦產(chǎn)勘查、油氣勘探、基礎(chǔ)地質(zhì)研究等領(lǐng)域(鐵、錳等礦產(chǎn)，能夠引起局部強(qiáng)磁異常，在磁法測(cè)量中可以被快速圈出勘查靶區(qū)[1-4])。另外，航空磁測(cè)按地表一定高度飛行測(cè)量，大大減弱了淺地表地質(zhì)體及地表磁性不均的干擾，能更清晰地反映深處的地質(zhì)情況，因而廣泛應(yīng)用于地球深部研究。

地球的磁場(chǎng)分為穩(wěn)定磁場(chǎng)和變化磁場(chǎng)。變化磁場(chǎng)起源于地球的外部，是疊加在地球基本磁場(chǎng)之上的各種短期變化磁場(chǎng)。按照起源的不同，變化磁場(chǎng)表現(xiàn)為兩大基本類型：①平靜變化；②干擾變化(磁擾)。平靜變化的主要成分是磁靜日變化，在中低緯度地區(qū)，日變幅平均在30 nT ～40 nT左右，在赤道地區(qū)，日變幅可達(dá)200 nT。以一天的變化來(lái)講，在不同的地區(qū)其形態(tài)變化特征均為白天變化大，而夜間比較平靜。地磁日變場(chǎng)是一個(gè)矢量場(chǎng)，在高精度航磁測(cè)量中，需對(duì)航空磁測(cè)數(shù)據(jù)進(jìn)行磁日變校正。

地磁場(chǎng)除了每天連續(xù)出現(xiàn)的周日變化外，還會(huì)發(fā)生無(wú)規(guī)律的突然變化。采用地磁觀測(cè)設(shè)備在地面對(duì)大地磁場(chǎng)進(jìn)行連續(xù)觀測(cè)中，短時(shí)間內(nèi)測(cè)量?jī)x器觀測(cè)到的數(shù)nT無(wú)規(guī)律的變化稱磁擾。突然出現(xiàn)的磁擾，變化劇烈，形態(tài)復(fù)雜。十分強(qiáng)烈的磁擾現(xiàn)象稱為磁暴，一般在全球同時(shí)發(fā)生。通常小磁暴較多，而大磁暴較少，有相當(dāng)數(shù)量的磁暴具有27天重復(fù)出現(xiàn)的特征。

依據(jù)中華人民共和國(guó)地質(zhì)礦產(chǎn)行業(yè)標(biāo)準(zhǔn)—航空磁測(cè)技術(shù)規(guī)范( DZ /T 0142-2010)(簡(jiǎn)稱“規(guī)范”)[5]，磁暴期間不得進(jìn)行航空磁測(cè)。當(dāng)進(jìn)行高精度航磁測(cè)量時(shí)，磁日變記錄連續(xù)出現(xiàn)梯度變化大于1 nT/min時(shí)，應(yīng)密切注意其變化；當(dāng)連續(xù)出現(xiàn)梯度大于5 nT/3min的非線性變化時(shí)，應(yīng)停止飛行或事后補(bǔ)飛。

磁暴的發(fā)生對(duì)于任何一項(xiàng)磁法測(cè)量勢(shì)必會(huì)有影響，如何確定影響大小，能否用數(shù)據(jù)處理手段減小影響，引起了廣大學(xué)者的重視。徐東禮等[6]對(duì)線性和非線性磁擾情況下的航空磁測(cè)數(shù)據(jù)進(jìn)行了影響分析和評(píng)價(jià)，認(rèn)為在嚴(yán)重磁擾情況下的航空磁測(cè)數(shù)據(jù)，經(jīng)過(guò)磁日變校正和其他各項(xiàng)修正后，可獲得與磁靜日重復(fù)觀測(cè)結(jié)果基本一致的航磁△T數(shù)據(jù)。

近幾年來(lái)，核工業(yè)航測(cè)遙感中心在黑龍江省和內(nèi)蒙古自治區(qū)開(kāi)展了大量的高精度航磁測(cè)量工作，并且在測(cè)量時(shí)，遇到磁暴事件發(fā)生。之后，依據(jù)“規(guī)范”要求，在磁靜日進(jìn)行了重新測(cè)量。筆者利用磁暴時(shí)航磁測(cè)量的數(shù)據(jù)和非磁暴日的重復(fù)測(cè)量數(shù)據(jù)，對(duì)比分析了這些磁暴時(shí)測(cè)量的歷史航磁數(shù)據(jù)質(zhì)量，對(duì)磁暴期間航測(cè)數(shù)據(jù)的利用具有指導(dǎo)意義。

1 航磁數(shù)據(jù)的比較

1.1 數(shù)據(jù)的獲取

總計(jì)選取了不同作業(yè)區(qū)、不同年度6條受磁暴影響的測(cè)線(圖1)，通過(guò)對(duì)磁日變數(shù)據(jù)每3min磁場(chǎng)強(qiáng)度最大變化統(tǒng)計(jì)分析，發(fā)現(xiàn)當(dāng)日有磁暴發(fā)生。對(duì)連續(xù)出現(xiàn)地磁場(chǎng)梯度大于5 nT/3 min的非線性變化時(shí)間段的測(cè)線數(shù)據(jù)，選擇磁靜日進(jìn)行了重新測(cè)量，第二次測(cè)量時(shí)的航磁數(shù)據(jù)和磁日變數(shù)據(jù)均符合“規(guī)范”的要求。

1.2 磁暴事件下的航磁數(shù)據(jù)分析

評(píng)價(jià)航磁數(shù)據(jù)質(zhì)量需要考慮以下幾個(gè)方面：

①航磁數(shù)據(jù)的動(dòng)態(tài)噪聲；②磁日變數(shù)據(jù)的噪聲；③航磁ΔT數(shù)據(jù)第一次測(cè)量和第二次測(cè)量的重現(xiàn)性。針對(duì)磁暴時(shí)間段測(cè)量的航磁數(shù)據(jù)，質(zhì)量采用兩種方法進(jìn)行評(píng)價(jià)分析：①內(nèi)符合精度法;②航磁數(shù)據(jù)四階差分法計(jì)算的動(dòng)態(tài)噪聲。

1.2.1 內(nèi)符合精度

內(nèi)符合精度法是一種通過(guò)第一次測(cè)量和第二次測(cè)量的數(shù)據(jù)計(jì)算重復(fù)測(cè)量精度的方法。計(jì)算一條測(cè)線兩次測(cè)量數(shù)據(jù)的內(nèi)符合精度均方誤差公式為[7-10]：

(1)

式中：n是重復(fù)線公共線段數(shù)據(jù)上的對(duì)應(yīng)點(diǎn)數(shù)量；δi是重復(fù)線公共線段各點(diǎn)測(cè)量值ΔTi與該點(diǎn)平均值ΔT的差。

依據(jù)航空磁測(cè)技術(shù)規(guī)范，航磁數(shù)據(jù)經(jīng)過(guò)了航磁數(shù)據(jù)和磁日變數(shù)據(jù)合并、飛行方向差校正、飛行海拔高度校正、地磁正常場(chǎng)校正、磁日變校正等處理。航磁測(cè)量時(shí)，航磁數(shù)據(jù)和坐標(biāo)數(shù)據(jù)采樣頻率為10 Hz。在基本保持和采樣點(diǎn)距一致的情況下，對(duì)坐標(biāo)數(shù)據(jù)和航磁ΔT數(shù)據(jù)進(jìn)行線性內(nèi)插，實(shí)際內(nèi)插點(diǎn)距在6.5 m～6.8 m之間。之后進(jìn)行了兩次測(cè)量公共線段提取及內(nèi)符合精度計(jì)算。

選取的6條測(cè)線，從圖1可以看出，測(cè)線上地磁場(chǎng)強(qiáng)度的變化各不相同。出現(xiàn)磁暴時(shí)間段，有時(shí)飛機(jī)正好位于地磁場(chǎng)平緩區(qū)，有時(shí)位于地磁場(chǎng)變化復(fù)雜區(qū)。對(duì)6條重飛測(cè)線進(jìn)行統(tǒng)計(jì)，研究其變化規(guī)律。通過(guò)計(jì)算(表1)可以看出，測(cè)線上非磁暴時(shí)間段所有重復(fù)測(cè)點(diǎn)參與計(jì)算時(shí)，內(nèi)符合精度在0.27nT～8.65 nT之間變化，不同地磁場(chǎng)梯度情況下，內(nèi)符合精度各不相同。剔除磁場(chǎng)梯度大于100 nT/km之后計(jì)算得出的內(nèi)符合精度在0.27 nT～3.53 nT之間(計(jì)算時(shí)磁暴時(shí)間段內(nèi)的數(shù)據(jù)沒(méi)有參加統(tǒng)計(jì))，內(nèi)符合精度明顯改善，可以看出磁場(chǎng)梯度變化大時(shí)，影響內(nèi)符合精度的計(jì)算結(jié)果。

圖1 航磁ΔT剖面圖Fig.1 Profile map of aeromagnetic ΔT(a)4971測(cè)線航磁ΔT剖面圖；(b)4130測(cè)線航磁ΔT剖面圖；(c)2770測(cè)線航磁ΔT剖面圖;(d)1720測(cè)線航磁ΔT剖面圖；(e)7360測(cè)線航磁ΔT剖面圖；(f)7370測(cè)線航磁ΔT剖面圖

通過(guò)分析磁暴時(shí)間段航磁ΔT數(shù)據(jù)(表2)可以看出，所有重復(fù)測(cè)點(diǎn)參與計(jì)算時(shí)內(nèi)符合精度在0.13 nT～8.08 nT之間變化，不同地磁場(chǎng)梯度情況下，內(nèi)符合精度也各不相同。剔除磁場(chǎng)梯度大于100 nT/km之后計(jì)算得出的內(nèi)符合精度在0.13 nT～2.78 nT之間，同樣可以明顯看出磁場(chǎng)梯度變化大距離為磁日變基站點(diǎn)距離測(cè)線最遠(yuǎn)點(diǎn)的距離；ε1為非磁暴時(shí)間段所有測(cè)點(diǎn)參加計(jì)算的內(nèi)符合精度；ε2為剔除磁場(chǎng)梯度大于100 nT/km之后計(jì)算的內(nèi)符合精度，點(diǎn)數(shù)為磁場(chǎng)梯度小于100 nT/km測(cè)點(diǎn)數(shù)；H1為測(cè)線第一次測(cè)量的平均離地高度；H2為測(cè)線第二次測(cè)量平均離地高度；Δ1為兩次測(cè)量的航磁ΔT平均值差。

表1 測(cè)線內(nèi)符合精度表

表2 磁暴時(shí)間段內(nèi)測(cè)線內(nèi)符合精度結(jié)果表

D為磁暴時(shí)間段內(nèi)的測(cè)線長(zhǎng)度，ε3為磁暴時(shí)間段內(nèi)的所有測(cè)點(diǎn)計(jì)算得出的內(nèi)符合精度；ε4為磁暴時(shí)間段內(nèi)的剔除磁場(chǎng)梯度大于100 nT/km之后計(jì)算的內(nèi)符合精度;ε5為磁暴時(shí)間段磁場(chǎng)平緩區(qū)內(nèi)測(cè)點(diǎn)計(jì)算的內(nèi)符合精度；Δ2為兩次測(cè)量的航磁ΔT平均值差。

時(shí)，影響航磁ΔT數(shù)據(jù)內(nèi)符合精度的計(jì)算結(jié)果。磁暴時(shí)間段內(nèi)磁場(chǎng)平緩區(qū)內(nèi)測(cè)點(diǎn)計(jì)算的內(nèi)符合精度在0.13 nT～1.08 nT之間，L1720和L2770均位于磁場(chǎng)變化復(fù)雜區(qū)，無(wú)法進(jìn)行磁場(chǎng)平緩區(qū)內(nèi)測(cè)點(diǎn)統(tǒng)計(jì)。

表1和表2均統(tǒng)計(jì)了6條測(cè)線兩次測(cè)量ΔT的平均值差，非磁暴時(shí)間段和磁暴時(shí)間段內(nèi)兩次測(cè)量平均值變化分別在-0.85 nT～2.41 nT和-0.48 nT～2.90 nT之間。非磁暴時(shí)間段和磁暴時(shí)間段平均值差的變化(Δ1-Δ2)在-1.29 nT～0.26 nT變化，二者之間并沒(méi)有發(fā)生顯著變化。

內(nèi)符合精度的變化與磁場(chǎng)強(qiáng)度梯度的變化關(guān)系密切。以L4130為例(圖2)，內(nèi)符合精度采用分段計(jì)算的方式進(jìn)行統(tǒng)計(jì)，以150個(gè)測(cè)點(diǎn)為一組進(jìn)行滑動(dòng)統(tǒng)計(jì)，內(nèi)符合精度偏大的地方，往往出現(xiàn)在磁場(chǎng)強(qiáng)度梯度變化大的區(qū)域。從圖2中可以看出，磁場(chǎng)變化平緩區(qū)的內(nèi)符合精度遠(yuǎn)好于磁場(chǎng)變化復(fù)雜區(qū)計(jì)算的內(nèi)符合精度。

L2770磁暴時(shí)間段內(nèi)兩次測(cè)量ΔT的內(nèi)符合精度為2.70 nT(測(cè)線長(zhǎng)12.3 km)，主要是由于測(cè)線位于磁場(chǎng)變化劇烈區(qū)、兩次測(cè)量位置最大偏離30 m造成的。內(nèi)符合精度采用分段計(jì)算的方法，在一定范圍之內(nèi)最差為4.50 nT，位于磁場(chǎng)變化復(fù)雜區(qū)；最好為0.92 nT，位于磁場(chǎng)變化相對(duì)平緩區(qū)(測(cè)線長(zhǎng)6.4 km)。這說(shuō)明L2770測(cè)線磁暴時(shí)間段ΔT的內(nèi)符合精度變差，并不是由于磁暴引起。

1.2.2 航磁動(dòng)態(tài)噪聲

兩次測(cè)量的航磁數(shù)據(jù)進(jìn)行了動(dòng)態(tài)噪聲分析，采用動(dòng)態(tài)噪聲水平Si值來(lái)評(píng)價(jià)每條測(cè)線航磁測(cè)量數(shù)據(jù)質(zhì)量。選取采樣間隔為0.5 s的航磁數(shù)據(jù)，并且舍掉水平梯度大于600 nT/km異常上的測(cè)點(diǎn)值，按式(2)計(jì)算Si值：

(2)

圖2 4130測(cè)線航磁測(cè)量結(jié)果曲線圖(資料來(lái)源：核工業(yè)航測(cè)遙感中心實(shí)測(cè))Fig.2 Aeromagnetic survey curve of Line 4130(a)4130測(cè)線兩次測(cè)量航磁ΔT剖面圖；(b)4130測(cè)線兩次測(cè)量航磁ΔT內(nèi)符合精度曲線圖;(c)4130測(cè)線兩次測(cè)量飛行航跡偏差曲線圖；(d)4130測(cè)線兩次測(cè)量離地高度差值曲線圖;(e)4130測(cè)線兩次測(cè)量時(shí)磁日變變化曲線圖

表3列出了不同年度6條測(cè)線磁暴條件下和非磁暴條件下的航磁數(shù)據(jù)的動(dòng)態(tài)噪聲水平。從表3結(jié)果可以看出，非磁暴時(shí)間段兩次測(cè)量測(cè)線統(tǒng)計(jì)動(dòng)態(tài)噪聲最大差值為-0.105 7 nT；當(dāng)?shù)谝淮螠y(cè)量為磁暴時(shí)間段測(cè)量，第二次測(cè)量為非磁暴時(shí)間段測(cè)量時(shí)兩次測(cè)量測(cè)線統(tǒng)計(jì)動(dòng)態(tài)噪聲最大差值為0.112 47 nT。有3條測(cè)線(L1720、L2770、L7360)在有磁暴期間飛行時(shí)動(dòng)態(tài)噪聲相對(duì)偏大，但也出現(xiàn)兩條測(cè)線(L4971、L7370)動(dòng)態(tài)噪聲略微偏小的情況。這6條測(cè)線兩次測(cè)量航磁數(shù)據(jù)動(dòng)態(tài)噪聲水平?jīng)]有明顯變化。

影響航磁數(shù)據(jù)動(dòng)態(tài)噪聲的因素包括：①儀器性能；②飛行離地高度；③地下磁性體的分布特征；④磁日變特征。因此第一次測(cè)量和第二次測(cè)量的結(jié)果不會(huì)完全一致，每條測(cè)線的航磁數(shù)據(jù)動(dòng)態(tài)噪聲也不相同。

2 非磁暴條件下測(cè)量的數(shù)據(jù)的比較

為了和磁暴時(shí)間段內(nèi)測(cè)量的數(shù)據(jù)精度進(jìn)行對(duì)比，同樣對(duì)均為非磁暴日的兩次測(cè)量進(jìn)行了內(nèi)符合精度分析和計(jì)算。抽選了2011年-2015年不同作業(yè)區(qū)的6條測(cè)線進(jìn)行分析(表4)，由表4可以看出，所有重復(fù)測(cè)點(diǎn)參加計(jì)算時(shí)，內(nèi)符合精度可在1.17 nT～8.48 nT之間變化，剔除磁場(chǎng)梯度大于100 nT/km之后的測(cè)點(diǎn)計(jì)算得出的內(nèi)符合精度在0.87 nT～2.83 nT之間，磁場(chǎng)梯度變化大時(shí)，影響內(nèi)符合精度的計(jì)算結(jié)果。采用磁場(chǎng)變化平緩區(qū)測(cè)點(diǎn)參加計(jì)算得出的兩次測(cè)量ΔT內(nèi)符合精度變化在0.41 nT～0.99 nT之間。采用平緩區(qū)測(cè)點(diǎn)進(jìn)行數(shù)據(jù)質(zhì)量評(píng)價(jià)，反映了設(shè)備的儀器狀態(tài)、性能指標(biāo)和當(dāng)時(shí)的地磁環(huán)境。

表3 重復(fù)線測(cè)量航磁數(shù)據(jù)動(dòng)態(tài)噪聲統(tǒng)計(jì)表

Si1、Si3為第一次測(cè)量原始航磁數(shù)據(jù)動(dòng)態(tài)噪聲；Si2、Si4第二次測(cè)量原始航磁數(shù)據(jù)動(dòng)態(tài)噪聲；Si1-Si2和Si3-Si4均為兩次測(cè)量動(dòng)態(tài)噪聲的差值

表4 非磁暴條件下重復(fù)線測(cè)量航磁ΔT內(nèi)符合精度表

D指磁日變基站點(diǎn)距離測(cè)線最遠(yuǎn)點(diǎn)的距離；ε1為所有測(cè)點(diǎn)計(jì)算得出的內(nèi)符合精度；ε2為剔除磁場(chǎng)梯度大于100 nT/km之后計(jì)算得出的內(nèi)符合精度；ε3為磁場(chǎng)變化平緩區(qū)測(cè)點(diǎn)參加計(jì)算得出的內(nèi)符合精度；H1為第一次測(cè)量時(shí)的平均離地高度；H2為第二次測(cè)量時(shí)的平均離地高度

3 磁日變數(shù)據(jù)對(duì)比分析

牡丹江和佳木斯磁日變基站相距258 km，2011年6月10日至10月8日期間，在牡丹江和佳木斯兩地開(kāi)展了多次磁日變測(cè)量，磁日變數(shù)據(jù)變化規(guī)律基本一致。但隨著測(cè)量時(shí)間的變化，二者之間的變化量略有差異。

以2011年牡丹江和佳木斯兩地磁日變基站測(cè)量結(jié)果為例(圖3)，7月19日06時(shí)19分03秒開(kāi)始測(cè)量，15時(shí)50分00秒結(jié)束測(cè)量，從測(cè)量結(jié)果來(lái)看，磁日變數(shù)據(jù)每3 min的最大變化小于3 nT，兩地磁日變數(shù)據(jù)的變化量在-4.25 nT～4.67 nT之間。對(duì)這個(gè)時(shí)間段內(nèi)兩地的磁場(chǎng)強(qiáng)度進(jìn)行相關(guān)性分析，相關(guān)系數(shù)為0.982 8。

9月13日07時(shí)45分34秒開(kāi)始測(cè)量，12時(shí)21分58秒結(jié)束測(cè)量，從測(cè)量結(jié)果來(lái)看，兩地磁日變數(shù)據(jù)的變化量在-2.02nT～3.82 nT之間，最大變化達(dá)5.84 nT。對(duì)這個(gè)時(shí)間段內(nèi)兩地的磁場(chǎng)強(qiáng)度進(jìn)行相關(guān)性分析，相關(guān)系數(shù)為0.985 6。磁日變數(shù)據(jù)在11時(shí)39分44秒～11時(shí)42分42秒、12時(shí)05分16秒～12時(shí)06分08秒、12時(shí)07分51秒～12時(shí)09分36秒時(shí)間段每3 min磁場(chǎng)強(qiáng)度的最大變化大于5 nT，分別為7.27 nT、5.27 nT、5.91 nT，出現(xiàn)磁暴事件，但兩地觀測(cè)的磁日變數(shù)據(jù)的變化量并不大，分別為1.12 nT、1.43 nT、2.02 nT，均在正常變化范圍內(nèi)。對(duì)磁暴時(shí)間段內(nèi)兩地的磁場(chǎng)強(qiáng)度進(jìn)行相關(guān)性分析，相關(guān)系數(shù)為0.997 9。由圖3(b)可以看出，當(dāng)日發(fā)生磁暴時(shí)，兩地磁日變數(shù)據(jù)(磁場(chǎng)強(qiáng)度)的變化量并沒(méi)有增大。

因此，在一定距離范圍內(nèi)，磁日變基站選擇合適，兩地磁日變數(shù)據(jù)具有相同的變化規(guī)律，這也說(shuō)明當(dāng)磁暴發(fā)生時(shí)，通過(guò)磁日變校正可以消除由于磁場(chǎng)日變化對(duì)測(cè)量結(jié)果造成的影響。

圖3 佳木斯和牡丹江磁日變觀測(cè)站測(cè)量的磁日變數(shù)據(jù)對(duì)比圖(資料來(lái)源：核工業(yè)航測(cè)遙感中心實(shí)測(cè))Fig.3 Magnetic diurnal data comparison diagram of Jiamusi and Mudanjang base station(a)2011年7月19日測(cè)量的磁日變曲線及差值圖;(b)2011年9月13日測(cè)量的磁日變曲線及差值圖

4 重復(fù)線測(cè)量結(jié)果影響因素分析

影響因素主要包括以下幾個(gè)方面：

1)地形影響。兼顧航空放射性測(cè)量的航磁測(cè)量，采取沿地形起伏飛行，飛行離地高度通常在120 m左右。起伏飛行對(duì)磁異常測(cè)量結(jié)果影響明顯，其影響程度和起伏高差有很強(qiáng)的相關(guān)性[11]。

2)磁場(chǎng)變化劇烈程度。磁場(chǎng)梯度變化大的區(qū)域重復(fù)性相對(duì)較差，在此區(qū)域開(kāi)展重復(fù)性對(duì)比，會(huì)低估儀器的重復(fù)測(cè)量精度。

3)磁日變數(shù)據(jù)。選擇合適的磁日變觀測(cè)基站，并且基站到測(cè)區(qū)的距離盡可能近。詹志佳等[12]認(rèn)為,地磁總強(qiáng)度日變化的日變輻、形態(tài)與頻譜，在100 km～200 km的局部范圍內(nèi)基本相同，但在500 km大范圍內(nèi)則不同，在25°～40°緯度范圍內(nèi)，地磁總強(qiáng)度日變幅的緯度因子約1 nT/°～2 nT/°。雖然磁暴的產(chǎn)生是全球性的，具有同步性，受經(jīng)度差影響小，但受緯度差影響較大。地磁日變幅隨磁日變基站所處緯度的增加而增加，同緯度的磁日變站記錄的磁日變數(shù)據(jù)形態(tài)具有相近的變化規(guī)律[13-15]。牡丹江和佳木斯磁日變基站實(shí)測(cè)數(shù)據(jù)表明，在基站相距258 km時(shí)，磁日變數(shù)據(jù)變化規(guī)律是基本一致的。

4)離地飛行高度不同產(chǎn)生的影響。兩次測(cè)量結(jié)果對(duì)比時(shí),宜選取離地飛行高度接近的測(cè)線。

5)飛行航跡不一致產(chǎn)生的影響。兩次測(cè)量結(jié)果對(duì)比宜選取飛行航跡接近的測(cè)線。

6)儀器性能。

7)磁補(bǔ)償參數(shù)。

在實(shí)際工作中，磁場(chǎng)變化劇烈程度、離地飛行高度和飛行航跡是影響兩次測(cè)量的數(shù)據(jù)存在差異的主要原因。在排出上述影響后，兩次測(cè)量ΔT內(nèi)符合精度較大時(shí)，應(yīng)查明由什么原因引起。

5 結(jié)論

同一臺(tái)設(shè)備，在磁暴事件下測(cè)量的航磁數(shù)據(jù)，經(jīng)過(guò)磁日變校正、正常場(chǎng)校正、飛行高度校正等各項(xiàng)修正后，可獲得與在非磁暴日重復(fù)觀測(cè)基本一致的結(jié)果，平緩磁場(chǎng)區(qū)兩次測(cè)量結(jié)果對(duì)比好于磁場(chǎng)變化復(fù)雜地區(qū)，這與非磁暴條件下兩次測(cè)量的數(shù)據(jù)變化規(guī)律相同。

磁場(chǎng)變化復(fù)雜地區(qū)，測(cè)線兩次測(cè)量ΔT內(nèi)符合計(jì)算精度變化較大，不能真實(shí)反映數(shù)據(jù)的測(cè)量質(zhì)量和儀器性能。平緩磁場(chǎng)區(qū)磁暴時(shí)間段內(nèi)測(cè)點(diǎn)計(jì)算的內(nèi)符合精度在1 nT左右。

[1] 張津偉, 武力聰,楊春,等.幾種低空高精度航空磁測(cè)系統(tǒng)及找礦應(yīng)用分析[J].礦產(chǎn)與地質(zhì)，2014， 28(1): 124-128. ZHANG J W, WU L C, YANG C, et al. Several high resolution aeromagnetic systems surveyed at very low altitude and mine prospecting application analysis [J].Mineral Resources and Geology, 2014， 28(1): 124-128. (In Chinese)

[2] 史建民，陳行時(shí)，丁繼雙．黑龍江大興安嶺地區(qū)航磁異常特征及找礦效果[J]. 地質(zhì)與資源,2015,24(3): 271-276. SHI J M, CHEN X S, DING J S. Aeromagnetic anomaly characteristics and prospecting reflect of Daxinganling Region in Heilongjiang province [J]. Geology and Resources, 2015, 24(3): 271-276. (In Chinese)

[3] 崔志強(qiáng)，胥值禮，孟慶敏, 等．現(xiàn)行三類平臺(tái)航磁勘查系統(tǒng)特點(diǎn)及勘查效果評(píng)述[J]. 物探化探計(jì)算技術(shù),2015, 37(4): 437-443. CUI Z Q, XU Z L, MENG Q M, et al. Review on exploration effect and characteristics of aeromagnetic survey system based on current three types flying platform [J]. Computing Techniques for Geophysical and Geochemical Exploration, 2015, 37(4): 437-443. (In Chinese)

[4] 胥值禮，李軍峰，崔志強(qiáng)，等. 無(wú)人機(jī)航空磁測(cè)技術(shù)在多寶山整裝勘查區(qū)的應(yīng)用試驗(yàn)[J].物探化探計(jì)算技術(shù)，2016, 38(4): 501-506. XU Z L,, LI J F, CUI Z Q, et al. Application trial of UAV aeromagnetic survey technique in the Duobaoshan key geological exploration area[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2016, 38(4): 501-506. (In Chinese)

[5] 熊盛青，陳斌，趙百民，等． 中華人民共和國(guó)地質(zhì)礦產(chǎn)行業(yè)標(biāo)準(zhǔn)—航空磁測(cè)技術(shù)規(guī)范( DZ /T 0142-2010) [S]． 北京: 中國(guó)標(biāo)準(zhǔn)出版社,2010：3-8. XIONG S Q, CHEN B, ZHAO B M, et al. Geology and mineral resources industry standard of the People's Republic of China—Criterion of aeromagnetic survey ( DZ /T 0142-2010) [S]． Beijing: China Standard Publishing House, 2010:3-8. (In Chinese)．

[6] 徐東禮, 范正國(guó)，舒晴, 等. 兩種典型磁擾對(duì)航空磁測(cè)的影響[J]. 物探與化探, 2015,39(2): 362-365. XU D L，F(xiàn)AN Z G，SHU Q， et al． The influence of two typical kinds of magnetic disturbance on the airborne magnetic survey[J]．Geophysical and Geochemical Exploration，2015,39(2): 362-365. (In Chinese)

[7] 徐東禮，葉挺明，舒晴，等. 航磁重復(fù)線內(nèi)符合精度計(jì)算方法[J]. 物探與化探, 2016,40(1): 125-128. XU D L，YE T M，SHU Q，et al．The method of calculating internal accord accuracy for repeated lines in aeromagnetic survey[J]． Geophysical and Geochemical Exploration, 2016, 40(1): 125-128. (In Chinese)

[8] 徐東禮，駱遙，賈偉潔. 航磁重復(fù)線測(cè)量數(shù)據(jù)質(zhì)量評(píng)價(jià)方法研究[J]. 物探與化探, 2014, 38(4): 729-731. XU D L，LUO Y，JIA W J．A study of the evaluation method for data quality of aeromagnetic repeated linear measurement[J]． Geophysical and Geochemical Exploration, 2014, 38(4): 729-731. (In Chinese)

[9] 郭志宏, 熊盛青, 周堅(jiān)鑫,等.航空重力重復(fù)線測(cè)試數(shù)據(jù)質(zhì)量評(píng)價(jià)方法研究[J].地球物理學(xué)報(bào)，2008,51(5)：1538-1543. GUO Z H, XIONG S Q, ZHOU J X, et al. The research on quality evaluation method of test repeat lines in airborne gravity survey[J]. Chinese Journal of Geophysics，2008, 51(5)：1538-1543.(In Chinese)

[10]姜作喜, 張虹, 郭志宏. 航空重力測(cè)量?jī)?nèi)符合精度計(jì)算方法[J]. 物探與化探, 2010,34(5):672-676. JIANG Z X, ZHANG H, GUO Z H. The method for calculation of internal accord accuracy in airborne gravity survey[J]. Geophysical & Geochemical Exploration, 2010,34(5):672-676. (In Chinese)

[11]孟慶奎，朱宏偉，朱彥珍. 起伏飛行對(duì)航磁異常的影響研究[J]. 工程地球物理學(xué)報(bào), 2016,13(4):464-469. MENG Q K, ZHU H W, ZHU Y Z. The study of fluctuation flight on the influence of aeromagnetic anomaly [J].Chinese Journal of Engineering Geophysics, 2016, 13(4):464-469. (In Chinese)

[12]詹志佳，姚富鑫，M. JOHNSTON, 等. 中美地磁總強(qiáng)度日變化的分析比較[J]. 地震學(xué)報(bào), 1992, (14)1: 83-88． ZHAN Z J, YAO F X, M. JOHNSTON, et al. Analysis and comparison of the daily variation of total geomagnetic intensity about Sino-America [J].Acta Seismologica Sinica, 1992, (14)1: 83-88. (In Chinese)

[13]張敏，王喜珍，張文來(lái)，等．感應(yīng)式磁力儀磁暴數(shù)據(jù)分析[J]. 地球物理學(xué)進(jìn)展, 2014,29(4)：1966-1972. ZHANG M, WANG X Z, ZHNG W L, et al. Storm different in inductive magnetometer and fluxgate magnetometer [J].Progress in Geophysics, 2014,29(4) : 1966-1972. (In Chinese)

[14]林吉綏．廣州和粵東地區(qū)地磁日變的對(duì)比分析[J]. 熱帶海洋, 1984,3(1): 10-16. LIN J S. A comparative analysis of geomagnetic diurnal variation between Guangzhou and eastern Guangdong region [J].Tropic Oceanology, 1984, 3(1): 10-16. (In Chinese)

[15]徐行, 廖開(kāi)訓(xùn), 陳邦彥，等．多臺(tái)站地磁日變觀測(cè)數(shù)據(jù)對(duì)遠(yuǎn)海磁測(cè)精度的影響分析[J].海洋測(cè)繪, 2007,27(1): 38-40. XU X, LIAO K X, CHEN B Y, et al. The effect analysis of the observed diurnal magnetic variation from multi-stations on the accuracy of marine magnetic survey [J].Hydrographic Surveying and Charting, 2007, 27(1): 38-40. (In Chinese)