不同型號磁通門磁力儀觀測資料質(zhì)量分析

胡秀娟，陳曉燕，王利兵，王 靜

(1.河北省地震局紅山基準(zhǔn)地震臺，河北 邢臺 054000；2.河北省地震局，河北 石家莊 050022)

0 引言

提供準(zhǔn)確可靠、連續(xù)完整的觀測資料是地磁臺站的首要任務(wù)。目前，磁通門磁力儀已成為地磁記錄的通用觀測儀器，廣泛應(yīng)用于地磁觀測臺站。磁通門磁力儀是矢量磁力儀，既可以測量磁場的方向，也可同時督視探頭溫度[1]。“十五”期間，地磁臺網(wǎng)中心為地磁臺站配備國產(chǎn)GM4型磁通門磁力儀；在2013年地磁臺網(wǎng)儀器更新改造項目期間，配備FGM-01型磁通門磁力儀；“十一五”背景場項目期間，配置GM4-XL型磁通門磁力儀。地磁臺站的相對記錄一般為2套儀器并行觀測，便于鑒別觀測設(shè)備故障和環(huán)境干擾，確保觀測數(shù)據(jù)的連續(xù)性和完整性。隨著磁通門磁力儀器的升級更新改造，觀測技術(shù)有明顯的提高。

磁通門磁力儀受探頭正交度、格值線性度、溫度系數(shù)、安裝定向準(zhǔn)確度等影響[2]。胡秀娟等利用本臺多套磁通門磁力儀實驗得到儀器定向誤差對地磁日變化記錄準(zhǔn)確度的影響[3]。

文章運用河北省3個地磁基準(zhǔn)臺站2016年的觀測數(shù)據(jù)，研究3種不同型號磁通門磁力儀受溫度及噪聲的影響情況。溫度對探頭內(nèi)鐵芯磁性材料、模擬電路電子元器件以及線圈性能產(chǎn)生影響，主要表現(xiàn)為：如果溫度漂移誤差較大，直接影響探頭各軸刻度因子和零偏值[4-5]，導(dǎo)致測量值與真實值之間的明顯偏差，這些觀測誤差在基線值中反映出來。為更好地對比分析，將陜西省乾陵基準(zhǔn)地震臺部分觀測數(shù)據(jù)一并引入。

研究結(jié)果可對臺站投入運行的磁通門類儀器進(jìn)行數(shù)據(jù)質(zhì)量評估，為改善觀測數(shù)據(jù)質(zhì)量提供參考。

1 臺站及觀測儀器概況

紅山基準(zhǔn)地震臺位于河北省隆堯縣山口鎮(zhèn)北部，始建于1966年邢臺地震后，地處太行山山前斷裂帶—隆堯、束鹿斷裂南端，臺基為震旦紀(jì)石英砂巖。

昌黎地震臺位于河北省昌黎縣城西10 km，龍家店鎮(zhèn)北, 始建于1969年，地處寧-昌斷裂北側(cè)，燕山山前沖積平原上，屬國家基本地磁臺。

涉縣地震臺位于河北省涉縣河南店鎮(zhèn)河二村，始建于1976年，地處涉縣盆地南部邊緣并靠近基巖出露的地方，是晉冀豫三省交界處重要的綜合觀測臺站。在中國地震背景場項目中，新建地磁相對記錄室和絕對觀測室，屬國家Ⅰ類地磁臺。

乾陵地震臺隸屬陜西省地震局，位于乾陵東乳峰南側(cè)，距縣城約8 km，始建于1966年，地處渭河斷陷盆地中段與陜北黃土高原的過渡地帶，臺基巖性為奧陶系石灰?guī)r，屬國家基本臺。

所選的4個臺站觀測環(huán)境良好，均可以排除環(huán)境對儀器的影響，且觀測數(shù)據(jù)在中國地震局的資料評比中均名列前茅，能為研究結(jié)果的可信度提供基本保障。每個臺站相對記錄室均采用地下式，同時配備兩套磁通門磁力儀與地磁絕對觀測儀，構(gòu)成完整的地磁數(shù)字化觀測系統(tǒng)。各儀器性能指標(biāo)參數(shù)如第6頁表1所示。

表1 地磁臺站儀器情況列表Table 1 Instruments in geomagnetic station

2 數(shù)據(jù)處理與分析

2.1 溫度的影響及計算方法

溫度系數(shù)是指當(dāng)磁通門磁力儀的輸出隨溫度線性變化時，與單位溫度變化對應(yīng)的測量值的變化。用于地磁臺站連續(xù)記錄的磁通門磁力儀一般要求溫度系數(shù)不大于1 nT/℃，由于各地磁臺站的觀測條件與實驗室有一定差距，加之儀器個體差異，在實際工作中儀器的各測向不一定都能達(dá)到這一指標(biāo)。

數(shù)據(jù)處理采用二元回歸分析方法。根據(jù)理論分析和實踐經(jīng)驗可知，影響基線值的主要因素是磁通門磁力儀的溫度系數(shù)和滑動量(滑動是隨時間單調(diào)增加或者減小的一種漂移，是時間的函數(shù))，基線值的回歸方程為：

HB(i)=a0+a1T(i)+a2t(i)+ε，(i=1,2,…N),

(1)

式(1)中：LHB(i)為某分量某一次的基線值；T代表溫度；t代表時間；ε為剩余誤差；N為觀測次數(shù)[5]。一般情況下，可用全年剩余均方差衡量基線值質(zhì)量，如式(2)所示。

(2)

式中：m為自變量個數(shù)。

選取2016年紅山臺GM4型、昌黎臺FGM-01型、涉縣臺和乾陵臺GM4-XL型磁通門磁力儀的基線值數(shù)據(jù)和溫度數(shù)據(jù)，如圖1、圖2、第7頁圖3和圖4所示。

由圖1可見，紅山臺相對記錄室內(nèi)年溫差≤8.0 ℃。D和H分量基線值與溫度曲線呈明顯的正相關(guān)性，Z分量則呈明顯的負(fù)相關(guān)性；由圖2可見，昌黎臺相對記錄室內(nèi)年溫差≤8.0 ℃，D、H和Z分量基線值與溫度曲線均呈明顯的正相關(guān)性；由第7頁圖3可見，涉縣臺相對記錄室內(nèi)年溫差≤3.0 ℃，從外在形態(tài)上，不能看出D、H和Z分量基線值與溫度曲線的相關(guān)性；由第7頁圖4可見，乾陵臺相對記錄室內(nèi)年溫差≤6.0 ℃，D分量與溫度曲線呈明顯的正相關(guān)性，H和Z分量的相關(guān)性不明顯。

圖1 2016年紅山臺各分量基線值與溫度曲線圖Fig.1 Baseline value of each component and temperature at Hongshan station in 2016

圖2 2016年昌黎臺各分量基線值與溫度曲線圖Fig.2 Baseline value of each component and temperature at Changli station in 2016

圖3 2016年涉縣臺各分量基線值與溫度曲線圖Fig.3 Baseline value of each component and temperature at Shexian station in 2016

通過式1和式2對4個臺站4套儀器2016年的基線值數(shù)據(jù)進(jìn)行溫度系數(shù)、滑動量和剩余均方差的計算，結(jié)果如表2所示(溫度系數(shù)單位為nT/℃，滑動系數(shù)和剩余均方差單位為nT/a)。昌黎臺2月27日儀器死機(jī)重啟，乾陵臺7月28日更換UPS，分別造成H分量和Z分量基線值產(chǎn)生臺階，先將臺階修正，然后進(jìn)行計算。

圖4 2016年乾陵臺各分量基線值與溫度曲線圖Fig.4 Baseline value of each component and temperature at Qianling station in 2016

表2 各臺站儀器的計算參數(shù)Table 2 The calculation parameters of instruments at each station

由上表看出，溫度系數(shù)昌黎臺FGM-01儀D和Z分量的最大，分別為2.013 1 nT/℃和2.543 0 nT/℃，紅山臺GM4儀和涉縣臺GM4-XL儀H分量的溫度系數(shù)相當(dāng)。對于滑動系數(shù)，昌黎臺FGM-01儀Z分量的滑動系數(shù)最大，為-7.460 9 nT/a；其次是紅山臺GM4儀的H分量，為-4.112 7 nT/a。對于剩余均方差，昌黎臺FGM-01儀各分量的較紅山臺GM4儀、涉縣臺和乾陵臺GM4-XL儀稍大。總體來看，兩套GM4-XL儀的溫度系數(shù)和剩余均方差要小于FGM-01儀和GM4儀。

綜上所述，若基線值的變化除人為誤差外，主要受溫度和滑動影響，剩余標(biāo)準(zhǔn)差就是觀測誤差。通過表2也可以看出，3類儀器個別分量滑動系數(shù)較大，但是各分量的剩余均方差較小，說明剩余均方差與溫度系數(shù)和滑動系數(shù)的大小無明顯的相關(guān)性。盡管儀器各分量會受溫度變化的影響，但經(jīng)過處理可以消除，同樣保證了觀測資料的質(zhì)量[6]。

2.2 背景噪聲

目前，國家地磁臺網(wǎng)中心在做質(zhì)量評價時，引進(jìn)臺站參考背景噪聲指標(biāo)。地磁觀測參考背景噪聲是指固定臺站某套儀器對磁場變化響應(yīng)的靈敏程度[7]。

選取紅山臺、昌黎臺、涉縣臺和乾陵臺4個臺站2016年秒數(shù)據(jù)背景噪聲值(見表3)。數(shù)據(jù)資料由國家地磁臺網(wǎng)中心網(wǎng)站提供。

表3 2016年各臺站儀器秒數(shù)據(jù)背景噪聲值Table 3 Background noise values of second data at each station in 2016

由表3看出，河北省內(nèi)3個臺站的3種儀器中，涉縣臺GM4-XL儀每月各分量背景噪聲均小于紅山臺和昌黎臺各套儀器的。由于涉縣臺只配備GM4-XL儀，無法證明是儀器本身還是環(huán)境的原因所致。因此，將有該型號儀器的陜西乾陵臺背景噪聲一并引入。對于乾陵臺同一臺站同一觀測室內(nèi)不同型號的兩套磁通門磁力儀，GM4-XL儀各分量背景噪聲要明顯小于GM4儀，H分量尤為明顯，全年恒定值為0 nT。

將表3中GM4儀、GM4-XL儀和FGM-01儀12個月的背景噪聲取平均值，作圖觀察不同型號儀器背景噪聲的變化情況(見圖5、圖6和第9頁圖7)。

圖5 2016年各儀器D分量背景噪聲曲線圖Fig.5 Background noise curve of D component of each instrument in 2016

圖6 2016年各儀器H分量背景噪聲曲線圖Fig.6 Background noise curve of H component of each instrument in 2016

由圖5、圖6和圖7可見，GM4-XL儀相比GM4儀和FGM-01儀，D、H和Z分量的背景噪聲值均為最小，H分量最明顯。其中，D分量和H分量各月背景噪聲的平均值大多保持在0.03 nT左右，Z分量的平均值大多保持在0.05 nT左右；對于FGM-01型和GM4型儀，各分量背景噪聲值相當(dāng)。FGM-01儀D分量各月的背景噪聲保持在0.04 nT左右，相比GM4儀，只有6月和7月大于GM4儀；H分量各月的背景噪聲只在3月、11月和12月小于GM4儀；Z分量各月的背景噪聲保持在0.07 nT左右，全年各月份均大于GM4儀。

3 結(jié)論

圖7 2016年各儀器Z分量背景噪聲曲線圖Fig.7 Background noise curve of Z component of each instrument in 2016

基線值是相對記錄儀器和絕對觀測儀器工作狀態(tài)的綜合反映，是檢驗地磁觀測資料內(nèi)在質(zhì)量的主要標(biāo)準(zhǔn)。地磁臺站是通過基線值監(jiān)視和控制相對記錄儀器的溫度特性和漂移特性，以保證測值的正確性[8]。通過對4個地磁臺站3種不同型號磁通門磁力儀2016年基線值進(jìn)行分析發(fā)現(xiàn)，F(xiàn)GM-01儀D分量和Z分量的溫度系數(shù)值較GM4儀和GM4-XL儀稍大，GM4儀和GM4-XL儀各分量溫度系數(shù)值相當(dāng)。雖然3類儀器個別分量滑動系數(shù)較大，但各分量的剩余均方差較小，說明剩余均方差與溫度系數(shù)和滑動系數(shù)的大小無明顯的相關(guān)性，也表明觀測資料的可靠性。

通過對4個地磁臺站3種不同型號磁通門磁力儀2016年秒數(shù)據(jù)背景噪聲進(jìn)行分析表明，GM4-XL儀各分量的背景噪聲均小于GM4型和FGM-01型儀；FGM-01型和GM4型儀各分量的背景噪聲值相當(dāng)。

對于D分量，F(xiàn)GM-01儀背景噪聲要略小于GM4儀，H和Z分量，F(xiàn)GM-01儀背景噪聲要略大于GM4儀。對陜西乾陵臺同一臺站同一觀測室內(nèi)不同型號的兩套磁通門磁力儀的背景噪聲對比，進(jìn)一步排除了環(huán)境的影響，認(rèn)為GM4-XL儀背景噪聲小是儀器本身造成的。

綜上所述，地磁臺站在日常工作中，可進(jìn)行溫度系數(shù)、滑動系數(shù)及剩余均方差的計算，以便及時了解前兩項參數(shù)對觀測儀器的影響變化及基線值觀測質(zhì)量的好壞。通過對儀器背景噪聲的計算，同樣可以監(jiān)控儀器的工作狀態(tài)和觀測環(huán)境的變化。

參考文獻(xiàn)：

[1] 王曉美，騰云田，王喜珍，等.GM4型磁通門磁力儀性能檢測方法[J].地震地磁觀測與研究，2008，29(6):88-94.

[2] 張素琴，楊冬梅.地磁日變化記錄準(zhǔn)確度標(biāo)定結(jié)果分析[J].地震地磁觀測與研究，2011，32(5)：60-65.

[3] 胡秀娟，李西京，王靜，等.儀器定向誤差對地磁日變化記錄準(zhǔn)確度的影響研究[J].地震學(xué)報，2016，38(1):130-137.

[4] 李久春.三軸磁敏傳感器誤差分析與校正研究[J].測控技術(shù)，2008，27(8):11-16.

[5] 武修文，張松勇，薛析蒙，等.基于磁通門與加速度傳感器的三分量磁力儀研究[J].現(xiàn)代電子技術(shù)，2013，36(22):96-98.

[6] 王利森，張云昌，邊鵬飛，等.利用逐步回歸分析紅山地磁臺觀測數(shù)據(jù)質(zhì)量[J].華北地震科學(xué)，2013，31(3):54-56.

[7] 胡秀娟，張素琴，何宇飛，等.地磁觀測參考背景噪聲分析[J].地震地磁觀測與研究，2012，33(1):124-129.

[8] 高玉芬，周錦屏，趙永芬.地磁臺絕對觀測數(shù)據(jù)分析[J].地震地磁觀測與研究，1990，11(3)：7-10.