王秋根　吝利民　張聰聰　張慶揚

摘要文安臺 FHD-2B 質(zhì)子磁力儀噪聲值在2020—2021年出現(xiàn)地磁 F 分量觀測數(shù)據(jù)噪聲偏大。其噪聲與環(huán)境溫度存在高相關(guān)性，表現(xiàn)為冬季噪聲大，夏季噪聲小的現(xiàn)象。我們依次對磁房周圍環(huán)境、供電電壓、接地電阻、供電線路是否漏電、探頭是否漏油、工頻干擾等進行排查，未能發(fā)現(xiàn)明顯干擾源。通過研究 FHD 質(zhì)子磁力儀的原理發(fā)現(xiàn)，F(xiàn)HD 質(zhì)子磁力儀的配諧電容 C 容量會隨溫度變化，溫度低時容量小，溫度高時容量大。配諧電容 C 容量變化會影響振蕩回路中心頻率f0，導(dǎo)致噪聲增大。由于文安臺 FHD 質(zhì)子磁力儀噪聲與溫度有關(guān)聯(lián)，我們推測 FHD 噪聲增大可能與儀器配諧電容 C 容量變化有關(guān)。通過調(diào)節(jié)探頭電感量 L 可修正配諧電容容量變化對回路中心頻率f0的影響。經(jīng)多次實驗，精細調(diào)整探頭電感量 L，當選定34.4 mH為探頭電感量值時，發(fā)現(xiàn)地磁 F 分量的噪聲值明顯下降，數(shù)據(jù)曲線變光滑，噪聲干擾排除。

關(guān)鍵詞 FHD-2B 質(zhì)子磁力儀；數(shù)據(jù)噪聲；電感量

中圖分類號： P315.62文獻標識碼： A文章編號：2096-7780（2023）07-0303-06

doi：10.19987/j.dzkxjz.2022-174

Analysis on noise interference ofFHD proton magnetometeratWenan seismic station

Wang Qiugen1），Lin Limin2），Zhang Congcong2），Zhang Qingyang1）

1） Baoding Earthquake Monitoring Center Station， Hebei Baoding 071000， China

2） Linfen Earthquake Monitoring Center Station， Shanxi Linfen 041000， China

AbstractThe FHD-2B proton magnetometer at Wenan station had large F component noise from 2020 to 2021. There is a high correlation between noise interference and temperature，which is generally characterized by high noise inwinter and low noise insummer. Thesurroundingenvironment，powersupply voltage，grounding resistance，powersupply line leakage，probe leakage and power frequency interference of the magnetic room were checked respectively， and no interference sources were found. By studying the principle of FHD proton magnetometer，it is found that theharmonic capacitance capacity of FHD proton magnetometer changes with temperature. The capacity is small at lowtemperature and large at high temperature. The change of the capacity of the harmonic capacitor will affect the centerfrequencyof theoscillationcircuit f0，resultinginincreasednoise. Sincethenoiseinterferenceof theFHDproton magnetometer at Wenan station is highly correlated with temperature，it is speculated that the noise interference may be related to the change in the capacity of the harmonic capacitance. We can adjust the inductance L of the probe to correct the influence of the circuit center frequency of the harmonic capacitance capacity transformation. The probe inductance was adjusted several times through experiments，and 34.4 mH was selected as the most appropriate probe inductance value. After the experiment，the noise value of F component decreased，the data curve became smooth，and the noise interference was solved.

KeywordsFHD-2B proton magnetometer; data noise; probe inductance

引言

FHD 質(zhì)子磁力儀是“九五”期間江蘇省地震局研制的分量質(zhì)子磁力儀，具有自動化、智能化和精度高的特點，該儀器自研發(fā)以來，已在全國多個地磁臺站成功運行[1]。噪聲可以有效監(jiān)測 FHD 質(zhì)子磁力儀的工作狀態(tài)，當儀器在觀測過程中出現(xiàn)問題時，噪聲會有很好的響應(yīng)。彭玉柱等[2]對地磁臺網(wǎng)2009—2014年的 FHD 觀測數(shù)據(jù)進行分析，發(fā)現(xiàn)參考背景噪聲冬季低夏季高，表現(xiàn)出明顯的季節(jié)變化特征。通過與基準站磁通門磁力儀噪聲對比，發(fā)現(xiàn)季節(jié)交替產(chǎn)生的溫度變化是引起儀器噪聲變化的主要原因。溫度變化或儀器老化對儀器配諧電容和探頭電感量產(chǎn)生影響，致使噪聲增大。張聰聰和靳寶萍[3]通過對臨汾中心地震臺 FHD 質(zhì)子磁力儀進行試驗，發(fā)現(xiàn)選擇合適的探頭電感量可以有效降低噪聲。張秀霞等[4]通過對比高郵臺質(zhì)子磁力儀觀測數(shù)據(jù)與當?shù)貧鉁匕l(fā)現(xiàn)，氣溫變化會對質(zhì)子磁力儀觀測數(shù)據(jù)產(chǎn)生影響，在氣溫高于35°和晝夜溫差變化較大時尤為明顯。董海龍等[5]通過分析嘉峪關(guān)臺 FHD 質(zhì)子磁力儀觀測數(shù)據(jù)噪聲，發(fā)現(xiàn)多臺儀器或多天線路間隔太小時，會在儀器觀測數(shù)據(jù)中疊加交流電磁干擾，致使地磁儀器噪聲增大。此外，F(xiàn)HD 質(zhì)子磁力儀觀測數(shù)據(jù)質(zhì)量也會受信號線的屏蔽性和長度的影響。

2020年文安臺 FHD 質(zhì)子磁力儀供電線路經(jīng)改造后，數(shù)據(jù)質(zhì)量較之前有了很大改善，能產(chǎn)出連續(xù)完整的產(chǎn)品數(shù)據(jù)。但每當進入冬季，地磁 F 分量數(shù)據(jù)噪聲就會偏大。為解決這一問題，有效降低噪聲，改善文安臺 FHD 質(zhì)子磁力儀數(shù)據(jù)質(zhì)量，本文從 FHD 數(shù)據(jù)噪聲產(chǎn)生原因入手分析研究。

1 文安臺概述

1.1 臺站地質(zhì)構(gòu)造

文安臺地處河北平原地震構(gòu)造帶北部，冀中拗陷中部，地殼厚度小于35 km，臺基為3000 m 第四系黃土覆蓋層。臺站位于勝芳—柳河—董村活動斷層西側(cè)，周邊有發(fā)育北東向牛東斷裂、任丘斷裂、大城斷裂和北西向徐水斷裂，縱橫交錯。氣候?qū)贉貛Ц蔁犷愋?，四季分明?月最熱，平均氣溫26.7℃;1 月最冷，平均氣溫?4.3℃ ， 年降水約300 mm。臺站附近先后于1967年發(fā)生河間 M6.3地震；1973年發(fā)生里坦 M5.4地震；2006年發(fā)生文安 M5.1地震。

1.2 臺站觀測環(huán)境

文安地震臺位于縣城郊南部，距縣城區(qū)僅700 m，占地約10畝?，F(xiàn)有觀測手段包括：FHD-2B 質(zhì)子磁力儀、氡釷分析儀、極低頻電磁測量儀、電磁擾動儀等。其中 FHD-2B 質(zhì)子磁力儀是 FHD 改進型質(zhì)子磁力儀，儀器安裝在臺站院內(nèi)西南拐角觀測房內(nèi)。臺站西側(cè)主要為農(nóng)田，南側(cè)與養(yǎng)老院為鄰，東側(cè)緊靠毛絨玩具廠，東面距縣氣象站約150 m，其余建筑以民房為主，無大型工礦企業(yè)，臺站 FHD-2B 型質(zhì)子磁力儀觀測環(huán)境良好（圖1）。地磁觀測室為磚混結(jié)構(gòu)，墻體較薄，保溫措施較差，室內(nèi)溫度隨季節(jié)變化顯著。

2 存在問題

文安臺 FHD 質(zhì)子磁力儀數(shù)據(jù)噪聲偏大問題出現(xiàn)在地磁 Z 分量和 F 分量數(shù)據(jù)中。從時間變化軸上看，集中在秋冬天氣偏冷月份，數(shù)據(jù)連續(xù)出現(xiàn)尖峰脈沖，偶然有小臺階。圖2為文安臺 FHD 質(zhì)子磁力儀數(shù)據(jù)噪聲與臺網(wǎng)均值對比圖，2020年春季氣溫低時，地磁 Z 分量和 F 分量數(shù)據(jù)質(zhì)量較差；數(shù)據(jù)質(zhì)量隨氣溫升高逐漸變好，至夏季天氣最熱時，噪聲干擾最??；氣溫降低后，數(shù)據(jù)質(zhì)量變差；2021年1月和2月天氣最冷時，噪聲干擾達到最大；隨著氣溫逐漸升高，噪聲干擾逐漸減??；10月份氣溫降低后，噪聲干擾又隨溫度減低逐漸增大。趨勢上顯示，文安臺 FHD 質(zhì)子磁力儀數(shù)據(jù)噪聲干擾與季節(jié)溫度有較高的相關(guān)度，總體上表現(xiàn)為冬季噪聲干擾大，夏季噪聲干擾小。

3 噪聲干擾分析

我們知道，造成數(shù)據(jù)噪聲增大的原因主要分兩類：儀器內(nèi)部噪聲和環(huán)境條件變化產(chǎn)生的噪聲。儀器噪聲主要為配諧不精確或設(shè)備老化所致。環(huán)境條件變化產(chǎn)生的噪聲主要有：磁性物質(zhì)干擾、磁房溫度變化、墩基不穩(wěn)、供電干擾等因素。

3.1FHD 質(zhì)子磁力儀原理

在地磁場中，具有自旋角動量和自旋磁矩的質(zhì)子會產(chǎn)生拉莫爾旋進運動，地磁場方向為此旋進的軸向，地磁場總強度 F 與旋進頻率f成正比：

式中，F(xiàn) 為地磁場總強度，單位nT；f 為質(zhì)子旋進頻率，單位 Hz；γp 為質(zhì)子旋磁比，γp=2.6751513×108 T · s ， 是推薦常數(shù)；π為圓周率，取值3.141592653。

磁力儀的探頭由線圈和無磁性外殼組成，線圈內(nèi)充滿高純度的煤油，油中含大量氫質(zhì)子，測量時，線圈中質(zhì)子旋進產(chǎn)生與旋進頻率相同的感應(yīng)電動勢。由于線圈中產(chǎn)生的電動勢僅在μV 級，非常微弱，需要進行選頻放大[6]。

質(zhì)子磁力儀使用 LC 并聯(lián)諧振電路選頻放大。中心頻率f0諧振電路關(guān)系式為：

式中，f0為振蕩回路中心頻率；L 為線圈電感量； C 為儀器配諧電容[7]；π與公式（1）相同。

只有感應(yīng)電動勢的旋進頻率 f 在振蕩回路的中心頻率f0附近時，質(zhì)子磁力儀才可選出需要的頻率信號進行放大測量。中心頻率f0越接近質(zhì)子旋進頻率 f，選出的頻率信號越準確，噪聲也越小。在實際電路中存在分布電阻和分布電容，同樣會對中心頻率f0產(chǎn)生影響。

儀器在長期觀測中，受環(huán)境溫度或儀器老化影響，儀器配諧電容 C 發(fā)生變化，導(dǎo)致實際質(zhì)子旋進頻率偏離振蕩回路的中心頻率，使觀測數(shù)據(jù)噪聲增大。由公式（2）可知，當儀器配諧電容 C 發(fā)生變化時，通過調(diào)整探頭電感量 L 可修正振蕩回路的中心頻率f0，使其保持始終接近質(zhì)子旋進頻率f。

3.2 噪聲干擾排查

通常情況下，儀器工作電壓過低、環(huán)境磁性物質(zhì)影響、接地不良、探頭漏油和外界工頻干擾都會影響到 FHD 質(zhì)子磁力儀觀測數(shù)據(jù)質(zhì)量。噪聲偏大可能是以上單個或多個原因疊加造成，如果臺站存在噪聲干擾，需對上述問題現(xiàn)象逐一排查。從磁房周圍環(huán)境進行排查，是否存在干擾源；檢查 FHD質(zhì)子磁力儀主機供電線路和信號線屏蔽在不同線槽內(nèi)，是否存在工頻干擾；測量儀器工作電壓和接地電阻值是否正常，符合要求；檢查供電外線路是否有漏電，儀器探頭是否有漏油情況。

FHD 質(zhì)子磁力儀配諧電容值變化會受環(huán)境溫度影響，冬季比夏季低5%左右。儀器的配諧電容值變化會影響儀器配諧頻率，導(dǎo)致儀器輸出噪聲增大。FHD 質(zhì)子磁力儀配諧電感量 L 是可調(diào)的，通過調(diào)節(jié)電感量 L 來修正配諧電路電容 C 變化量達到精確配諧中心頻率f0的目的。通常，冬季電感量 L 增大約0.2 mH，夏季減小約0.2 mH[8]。文安臺 FHD 質(zhì)子磁力儀噪聲干擾與環(huán)境溫度變化有很高的相關(guān)性，初步判斷噪聲偏大與儀器探頭電感量配選值是否合適有關(guān)。

4 調(diào)整探頭電感量實驗

4.1 實驗時間段選取

地磁日變化分為磁擾日變化和磁靜日變化。磁擾包括地磁脈動、磁暴、灣擾和鉤擾等，發(fā)生磁擾時地磁場變化幅度很大。磁靜日地磁場變化較小。其一般變化規(guī)律為白天變化大，晚上變化小。實驗時段選在磁靜日晚上進行，這個時段 Z 分量變化比較平緩，有利于監(jiān)控調(diào)整過程中 Z 分量的前后變化，避開下午極值變化。本次實驗選在2021年11月26日晚上6點鐘以后進行，該時間段磁場環(huán)境較為穩(wěn)定，具有較好的儀器背景噪聲水平。

4.2 調(diào)整探頭電感量

每次調(diào)整儀器探頭電感量后，即對儀器工作2小時左右測量數(shù)據(jù)進行分析。通過觀察記錄數(shù)據(jù)曲線光滑度和計算數(shù)據(jù)一階差分均方差，粗略判斷每次調(diào)整探頭電感量后觀測數(shù)據(jù)噪聲變化水平（表1和圖3）。一階差分均方差公式為：

式中，D 表示觀測數(shù)據(jù)一階差分均方差；x 表示記錄數(shù)據(jù)一階差分；E 表示測量數(shù)據(jù)的平均值。

實驗前文安臺 FHD 探頭電感量34.1 mH。最初調(diào)節(jié)電感量 L 微調(diào)增加0.2 mH至34.3 mH，計算一階差分均方差值時，發(fā)現(xiàn)儀器輸出噪聲值有所減小；然后，將調(diào)節(jié)電感量 L 減小0.2 mH至33.9 mH，計算一階差分均方差值發(fā)現(xiàn)儀器輸出噪聲有所增加。經(jīng)過反復(fù)多次實驗，得出儀器探頭電感量 L 真實值比實驗前設(shè)定電感量值大些。最后我們將電感量 L 調(diào)至34.5 mH，計算一階差分均方差減小，但儀器輸出噪聲比電感量 L 調(diào)至34.3 mH時要大。儀器探頭電感量 L 應(yīng)在34.2 mH、34.3 mH和34.4 mH其中選擇一個合適值，之后再將電感量 L 分別調(diào)到34.2 mH和 34.4 mH，計算一階差分均方差值差別不大，由于儀器噪聲干擾主要出現(xiàn)在 F 分量，故選定 F 分量噪聲值較小的34.4 mH為文安臺 FHD 儀器探頭電感量值。

經(jīng)調(diào)整電感量 L 實驗記錄數(shù)據(jù)，F(xiàn) 分量一階差分曲線明顯變細，F(xiàn) 分量噪聲從0.32 nT降到0.20 nT； Z 分量噪聲從0.36 nT降到0.27 nT，F(xiàn) 分量和 Z 分量噪聲明顯降低。圖4為文安臺 FHD 質(zhì)子磁力儀修正探頭電感量 L 前后 F 分量數(shù)據(jù)比較曲線，F(xiàn) 分量數(shù)據(jù)突跳幅度明顯降低，曲線變光滑。實驗說明 FHD 質(zhì)子磁力儀通過適當修正儀器探頭電感量 L 可有效減少儀器輸出噪聲水平，改善輸出數(shù)據(jù)質(zhì)量。

5 結(jié)束語

通過對文安臺 FHD 儀器噪聲增大存在的各種可能干擾原因逐一排查和實驗，發(fā)現(xiàn)文安臺 FHD 噪聲增大原因是臺站儀器未能選擇合適儀器探頭電感量 L 造成的。作者通過反復(fù)修正儀器探頭電感量 L 參數(shù)實驗，并對實驗數(shù)據(jù)曲線光滑度和計算數(shù)據(jù)一階差分均方差的分析，發(fā)現(xiàn)合適選擇儀器探頭電感量 L 參數(shù)可有效降低儀器輸出噪聲水平。實驗結(jié)果最終選定儀器探頭電感量為34.4 mH。使得 FHD 儀器 F 分量曲線變光滑、一階差分均方差減小，儀器輸出噪聲明顯降低。噪聲問題在臺站 FHD 質(zhì)子磁力儀觀測運行中經(jīng)常會遇到，本次噪聲問題解決經(jīng)驗可為其他地磁臺站提供借鑒和參考。

參考文獻

［1］何宇飛，夏忠，李軍，等. FHD 質(zhì)子磁力儀改進研究進展[J].國際地震動態(tài)，2019（11）：20-26

He Y F，Xia Z，Li J，et al. Progress of the improved FHD proton magnetometer[J]. Recent Developments in World Seismology，2019（11）：20-26

［2］彭玉柱，何宇飛，夏忠，等.溫度對 FHD 磁力儀觀測數(shù)據(jù)噪聲影響[J].地震地磁觀測與研究，2015，36（2）：103-107

Peng Y Z，He Y F，Xia Z，et al. The temperature influence on data noise of the FHD magnetometer[J]. Seismological and Geomagnetic Observation and Research，2015，36（2）：103-107

［3］張聰聰，靳寶萍.臨汾中心地震臺 FHD 質(zhì)子磁力儀噪聲干擾分析[J].山西地震，2020（3）：30-33

Zhang C C，Jin B P. Noise interference analysis ofFHD proton magnetometer at Linfen central seismic station[J]. Earthquake Research in Shanxi，2020（3）：30-33

［4］張秀霞，蔣延林，趙衛(wèi)紅，等.高郵地震臺質(zhì)子矢量磁力儀觀測與氣溫變化分析[J].山西地震，2017（4）：16-19

Zhang X X，Jiang Y L，Zhao W H，et al. Comparison analysis of proton vector magnetometer observations and temperature changes at Gaoyou seismic station[J]. Earthquake Research in Shanxi，2017（4）：16-19

［5］董海龍，夏忠，閆萬生，等.嘉峪關(guān)臺 FHD 觀測數(shù)據(jù)噪聲的影響因素及對策[J].甘肅科技，2016，32（1）：45-46

Dong H L，Xia Z，Yan W S，et al. Influencing factors and countermeasures of FHD observation data noise at Jiayuguan station[J]. Gansu Science and Technology，2016，32（1）：45-46

［6］居海華，夏忠. FHD-2B 分量質(zhì)子磁力儀維護技術(shù)[J].地震地磁觀測與研究，2012，33（5）：135-141

Ju H H，Xia Z. Maintenance technology for FHD-2B type component proton magnetometer[J]. Seismological and Geomagnetic Observa- tion and Research，2012，33（5）：135-141

［7］盧嘉沁，趙玉紅.都蘭地磁臺 FHD-2B 磁力儀觀測數(shù)據(jù)干擾分析[J].高原地震，2020，32（1）：52-57

Lu J Q，Zhao Y H. Interference analysis ofFHD-2B magnetometer observation data at Dulan geomagnetic station[J]. Plateau Earthquake Research，2020，32（1）：52-57

［8］呂永青，張亮娥，宮靜芝，等.昔陽臺 FHD-2B 質(zhì)子矢量磁力儀背景噪聲超標原因分析及降噪措施[J].山西地震，2020（4）：26-28，32

Lü Y Q，Zhang L E，Gong J Z，et al. Analysis on the cause of excessive background noise of Xiyangtai FHD-2B proton vector mag- netometer and its nosie reduction measures[J]. Earthquake Research in Shanxi，2020（4）：26-28，32