[摘要] """2014年9月山西省地震局與中國(guó)地震局地球物理研究所在山西北部恒山斷裂帶區(qū)域共同建設(shè)了地磁臺(tái)陣，布設(shè)了5臺(tái)GM4型磁通門磁力儀，安裝于代縣、山陰、應(yīng)縣、繁峙和懷仁，以獲取區(qū)域場(chǎng)地磁三分量觀測(cè)資料，臺(tái)陣的5個(gè)測(cè)點(diǎn)位于高壓直流輸電晉南線接地極附近。通過對(duì)比山西地磁臺(tái)陣5個(gè)測(cè)點(diǎn)GM4型磁通門磁力儀與其他臺(tái)站地磁秒數(shù)據(jù)，量算出高壓直流輸電線路（晉北—南京）對(duì)山西地磁臺(tái)陣地磁數(shù)據(jù)干擾幅度，對(duì)其進(jìn)行相關(guān)性和誤差分析，求得擬合公式。結(jié)果表明，山西地磁臺(tái)陣受晉南線高壓直流輸電干擾，Z分量變化幅度隨臺(tái)站位置至?xí)x南線垂直距離增加而減小， H、D分量干擾幅度與Z分量干擾幅度線性相關(guān)，以Z分量變化量為因變量，使用參考擬合公式求取的H、D分量變化量，可應(yīng)用于山西地磁臺(tái)陣受晉南線高壓直流輸電干擾的日常處理工作。

[關(guān)鍵詞] 高壓直流輸電; 線性分析; 干擾幅度" " " " " " " " " " [DOI] 10.19987/j.dzkxjz.2023-163

中圖分類號(hào)：P315.63"""""""""文獻(xiàn)標(biāo)識(shí)碼： A""""""""""文章編號(hào)： 2096-7780（2024）09-0584-06

0 "引言

在高壓直流輸電線路的中段，線路兩邊產(chǎn)生的磁場(chǎng)大小與輸電線路中的不平衡電流的大小成正比，與臺(tái)站距輸電線路的垂直距離成反比，方向相反，三者基本滿足畢奧—薩伐爾定律[1-2]，但是高壓直流輸電接地極附近的地磁臺(tái)站，不僅受高壓直流輸電線路的干擾，還受到接地極走線、接地極位置和周圍敏感建筑物的感應(yīng)場(chǎng)疊加干擾[3-4]，特別是水平分量H、磁偏角D的干擾形態(tài)相對(duì)復(fù)雜、幅度小，判別難度大。截止2023年8月全國(guó)在運(yùn)行高壓直流輸電線36 條，接地極附近少有觀測(cè)點(diǎn)。2017年3月晉南線調(diào)試運(yùn)行，該線路的起始點(diǎn)是朔州（換流站位于平魯區(qū)向陽堡鄉(xiāng)南漢井村，接地極位于忻州市神池縣八角鎮(zhèn)郭家村東約1.5 km），山西地磁臺(tái)陣位于恒山斷裂兩側(cè)，有代縣、繁峙、應(yīng)縣、山陰和懷仁5個(gè)測(cè)點(diǎn)（圖1），安裝有GM4磁通門磁力儀（GM4型磁通門磁力儀，是一種連續(xù)測(cè)量和記錄地磁場(chǎng)D、H、Z時(shí)間變化的數(shù)字化相對(duì)記錄儀，采樣率為1次/s，分辨力為0.1 nT，噪聲為lt;0.1 nT[5]），位于起始點(diǎn)附近，臺(tái)陣5個(gè)測(cè)點(diǎn)觀測(cè)數(shù)據(jù)為研究接地極附近高壓直流輸電干擾特征提供了基礎(chǔ)研究數(shù)據(jù)。

1 "山西地磁臺(tái)陣受晉南線高壓直流輸電干擾幅度資料選取與測(cè)算方法

2013年，地磁臺(tái)網(wǎng)高壓直流輸電判別處理系統(tǒng)的使用對(duì)高壓直流輸電干擾垂直分量Z提供了參考標(biāo)準(zhǔn)[6]。本文首先通過國(guó)家地磁臺(tái)網(wǎng)中心網(wǎng)站獲取晉南線高壓直流輸電對(duì)山西地磁臺(tái)陣臺(tái)站地磁垂直分量Z干擾的大小和方向，再使用對(duì)比法差值檢測(cè)（受干擾臺(tái)站與不受干擾臺(tái)站對(duì)比），測(cè)算出水平分量H、磁偏角D的干擾幅度和方向（GM4型磁通門磁力儀產(chǎn)出D分量、H分量、Z分量數(shù)據(jù)單位均為nT）。選取自2017年以來晉南線（晉北—南京）高壓直流輸電對(duì)山西地磁臺(tái)陣代縣、繁峙、應(yīng)縣、山陰和懷仁5個(gè)地磁臺(tái)各157個(gè)干擾事件，合計(jì)785個(gè)干擾事件進(jìn)行分析（表1）。

2 "山西地磁臺(tái)陣受晉南線高壓直流輸電干擾特征分析

（1）干擾幅度：山西地磁臺(tái)陣至?xí)x南線高壓直流輸電垂直距離P懷仁（68 km）gt;P繁峙（50 km）gt;P應(yīng)縣（40 km）gt;P山陰（15 km）gt;P代縣（lt;1 km），干擾幅度ΔZ（代縣）gt;ΔZ（山陰）gt;ΔZ（應(yīng)縣）gt;ΔZ（繁峙）gt;ΔZ（懷仁），山西地磁臺(tái)陣受晉南線高壓直流輸電干擾Z分量變化幅度隨臺(tái)站位置至?xí)x南線垂直距離增加而減小[7-8]。

（2）干擾方向：代縣、繁峙、應(yīng)縣、懷仁地磁臺(tái)受晉南線高壓直流輸電干擾，H分量、D分量干擾方向一致，但與Z分量干擾方向相反；山陰地磁臺(tái)受晉南線高壓直流輸電干擾，H分量干擾方向與Z分量干擾方向相反，D分量干擾方向部分干擾時(shí)段與Z分量干擾方向相同，部分干擾時(shí)段與Z分量干擾方向相反。

（3）干擾時(shí)段：從統(tǒng)計(jì)的樣本結(jié)果發(fā)現(xiàn)，除懷仁地磁臺(tái)Z分量受晉南線干擾幅度在?0.6～0.3 nT的時(shí)段，H、D分量不受干擾外，其他Z分量受晉南線干擾時(shí)段，代縣、山陰、繁峙、應(yīng)縣和懷仁地磁臺(tái)H、D分量均受到不同程度的干擾。

3 "山西地磁臺(tái)陣受晉南線高壓直流輸電干擾幅度線性擬合

我們定義高壓直流輸電干擾H、D、Z分量的變化量為ΔH、ΔD、ΔZ，通過繪制ΔH、ΔD與ΔZ散點(diǎn)圖，發(fā)現(xiàn)因變量與自變量的變化趨勢(shì)是線性的，因此，可以使用一元線性回歸分析ΔH、ΔD與ΔZ的關(guān)系[9-11]。以ΔH、ΔD為縱坐標(biāo)，ΔZ為橫坐標(biāo)，建立一元回歸方程，則是相應(yīng)于的回歸直線的理論值，而則是相應(yīng)于的實(shí)測(cè)值。本文使用最小二乘法求取a（斜率）、b（截距）、（剩余均方差） 、r（相關(guān)系數(shù)）4個(gè)參數(shù)。當(dāng)rgt;0.8時(shí)稱為高度相關(guān)，當(dāng)rlt;0.3時(shí)稱為低度相關(guān)，其他時(shí)段為中度相關(guān)[12]。本文使用k值評(píng)價(jià)擬合誤差，k越大誤差越大，k越小誤差越小[13]。山西地磁臺(tái)陣受晉南線高壓直流輸電干擾幅度散點(diǎn)圖與擬合曲線見圖2。山西地磁臺(tái)陣受高壓直流輸電干擾幅度參考擬合公式與擬合誤差見表2[14-15]。相關(guān)計(jì)算公式如下：

4 "結(jié)論與討論

通過對(duì)山西地磁臺(tái)陣受晉南線高壓直流輸電干擾分析，可以得出以下結(jié)論：

（1）干擾幅度：山西地磁臺(tái)陣臺(tái)站至?xí)x南線高壓直流輸電垂直距離P懷仁（68 km）gt;P繁峙（50 km）gt;P應(yīng)縣（40 km）gt;P山陰（15 km）gt;P代縣（lt;1 km），干擾幅度ΔZ（代縣）gt;ΔZ（山陰）gt;ΔZ（應(yīng)縣）gt;ΔZ（繁峙）gt;ΔZ（懷仁），山西地磁臺(tái)陣受晉南線高壓直流輸電干擾Z分量變化幅度隨臺(tái)站位置至?xí)x南線垂直距離增加而減小。H分量變化幅度和D分量變化幅度較為復(fù)雜，認(rèn)為與晉南線的走向和臺(tái)站至?xí)x南線垂直距離相關(guān)。

（2）相關(guān)性：ΔZ參考范圍內(nèi)，代縣、繁峙、應(yīng)縣和懷仁的ΔH、ΔD與ΔZ均高度負(fù)相關(guān)；ΔZ參考范圍內(nèi)，山陰ΔH與ΔZ高度負(fù)相關(guān)，ΔD與ΔZ中度正相關(guān)。

（3）擬合誤差：山西地磁臺(tái)陣D分量擬合誤差大于H分量擬合誤差，即k（D）gt; k（H）。

（4）擬合公式：山陰地磁臺(tái)ΔD擬合誤差較大，擬合公式不建議使用，山陰地磁臺(tái)D分量受晉南線高壓直流干擾，變化幅度復(fù)雜，初步認(rèn)為與山陰地磁臺(tái)距換流站和接地極較近有關(guān)，需在今后工作中進(jìn)一步核實(shí)。其他參考擬合公式可用于山西地磁臺(tái)陣受晉南線高壓直流干擾日常處理。由于樣本數(shù)量有限，公式存在偏差，在使用擬合公式時(shí)，應(yīng)嚴(yán)格按照ΔZ參考范圍進(jìn)行擬合計(jì)算。

參考文獻(xiàn)

沈紅會(huì)，周加新，張秀霞，等. 直流輸電對(duì)江蘇地磁觀測(cè)的影響分析[J]. 地震地磁觀測(cè)與研究，2005，26（5）：65-70 """Shen H H，Zhou J X，Zhang X X，et al. The influence analysis of direct current transmit electricity to geomagnetic observation[J]. Seismological and Geomagnetic Observation and Research，2005，26（5）：65-70

楊興悅，董海龍，王燕，等. 高壓直流輸電線對(duì)甘肅地磁觀測(cè)的影響分析及應(yīng)對(duì)措施[J]. 地震工程學(xué)報(bào)，2020，42（4）：899-906 """Yang X Y，Dong H L，Wang Y，et al. Impact and countermeasures of HVDC transmission on geomagnetic observation in Gansu Province[J]. China Earthquake Engineering Journal，2020，42（4）：899-906

高熹微，張秀霞，郭宇鑫，等. 地磁H、D分量受高壓直流輸電干擾形態(tài)分析[J]. 地震，2023，43（2）：131-144 """Gao X W，Zhang X X，Guo Y X，et al. Analysis of the interference form of geomagnetic H"and D"components by HVDC transmission[J]. Earthquake，2023，43（2）：131-144

張建濤，王明貴，陳星星，等. 高壓直流輸電對(duì)湖北省地磁觀測(cè)的影響[J]. 地震地磁觀測(cè)與研究，2023，44（3）：95-101 """Zhang J T，Wang M G，Chen X X，et al. Influence of HVDC transmission on geomagnetic observation in Hubei Province[J]. Seismological and Geomagnetic Observation and Research，2023，44（3）：95-101

王曉美，滕云田，王喜珍，等. GM4型磁通門磁力儀性能檢驗(yàn)方法[J]. 地震地磁觀測(cè)與研究，2008，29（6）：88-94 """Wang X M，Teng Y T，Wang X Z，et al. The test method for the performance of fluxgate magnetometer[J]. Seismological and Geomagnetic Observation and Research，2008，29（6）：88-94

陳俊，蔣延林，張秀霞，等. 地磁臺(tái)網(wǎng)高壓直流輸電判別處理系統(tǒng)設(shè)計(jì)[J]. 地震地磁觀測(cè)與研究，2014，35（增刊2）：270-274 """Chen J，Jiang Y L，Zhang X X，et al. The design of HVDC discrimination and processing system for geomagnetic network[J]. Seismological and Geomagnetic Observation and Research，2014，35（S2）：270-274

羅詞建，李媛媛，王秋寧，等. 高壓直流輸電對(duì)陜西地磁觀測(cè)的影響[J]. 地震地磁觀測(cè)與研究，2019，40（4）：74-80 """Luo C J，Li Y Y，Wang Q N，et al. Influence of HVDC transmission on observation of geomagnetic in Shaanxi[J]. Seismological and Geomagnetic Observation and Research，2019，40（4）：74-80

李霞，趙玉紅，張朋濤，等. 高壓直流輸電干擾對(duì)青海地磁觀測(cè)數(shù)據(jù)的影響特征[J]. 地震地磁觀測(cè)與研究，2022，43（增刊1）：229-231 """Li X，Zhao Y H，Zhang P T，et al. Characteristic analysis of the influence of HVDC interference on geomagnetic observation in Qinghai Province[J]. Seismological and Geomagnetic Observation and Research，2022，43（S1）：229-231

胡曉霞，譚麗. 線性回歸分析在煤質(zhì)檢測(cè)結(jié)果審核中的應(yīng)用[J]. 煤炭加工與綜合利用，2023（7）：63-67 """Hu X X，Tan L. Application of linear regression in the examination of coal quality test results[J]. Coal Processing amp; Comprehensive Utilization，2023（7）：63-67

唐亞寧，周丙常，謝文賢. 一元線性回歸方程的三種顯著性檢驗(yàn)的等價(jià)性分析[J]. 高等數(shù)學(xué)研究，2022，25（6）：76-78 """Tang Y N，Zhou B C，Xie W X. Equivalence analysis of three significant tests of univariate linear regression equation[J]. Studies in College Mathematics，2022，25（6）：76-78

張雪鋒，李忠，陳誠(chéng)，等. 基于Origin最小二乘法和解析法的土工液塑限聯(lián)合測(cè)定試驗(yàn)數(shù)據(jù)處理[J]. 土工基礎(chǔ)，2023，37（5）： 825-829，840 """Zhang X F，Li Z，Chen C，et al. Atterberg limits tests data processing by Origin based least square and analytical methods[J]. Soil Engineering and Foundation，2023，37（5）：825-829，840

萬永革. 數(shù)字信號(hào)處理的MATLAB實(shí)現(xiàn)[M]. 2版. 北京：科學(xué)出版社，2012：4-7 """Wan Y G. MATLAB implementation of digital signal processing[M]. 2nd ed. Beijing：Science Press，2012：4-7

范磊，魏越超，郭宇，等. 太原地震監(jiān)測(cè)中心站地磁數(shù)據(jù)受高壓直流輸電干擾的識(shí)別處理[J]. 山西地震，2022（4）：9-12，37 """Fan L，Wei Y C，Guo Y，et al. Identification and processing of geomagnetic data disturbed by HVDC transmission in Taiyuan Earthquake Monitoring Center Station[J]. Earthquake Research in Shanxi，2022（4）：9-12，37

羅杰，賈靜. 不同線性擬合方法對(duì)直剪試驗(yàn)結(jié)果影響的分析[J]. 水科學(xué)與工程技術(shù)，2023（6）：83-86 """Luo J，Jia J. The influence of different linear fitting methods on the results of direct shear tests[J]. Water Sciences and Engineering Technology，2023（6）：83-86

徐崢. 基于加權(quán)總體最小平方中值法的直線擬合[J]. 鐵道建筑技術(shù)，2023（12）：49-52，135 """Xu Z. Line fitting based on the weighted total least median of squares[J]. Railway Construction Technology，2023（12）：49-52，135

[Abstract] """"In September 2014， Shanxi Earthquake Agency and the Institute of Geophysics， China Earthquake Administration jointly completed the construction of the geomagnetic station array in the Hengshan fault zone in the north of Shanxi. The geomagnetic station array is composed of 5 GM4 type fluxgate magnetometers， which are respectively installed in Daixian， Shanyin， Yingxian， Fanchi and Huairen to obtain the three-component magnetic field data in the observation area. The 5 measurement points of the geomagnetic station array are located near the grounding electrode of the Jinnan line of high-voltage direct current transmission. By comparing the GM4 type fluxgate magnetometer at the 5 measurement points of the geomagnetic station array with the geomagnetic second data of other stations， the interference amplitude of the high-voltage direct current transmission line （Jinbei-Nanjing） on the geomagnetic data of the Shanxi geomagnetic station array was measured， correlation and error analysis also were carried out， and the fitting formula was obtained. The results show that the variation amplitude of the Z"component disturbed by the Jinnan line high-voltage direct current transmission in the Shanxi geomagnetic station array decreases as the vertical distance from the station location to the Jinnan line increases， and the interference amplitudes of the H"and D"components are linearly correlated with the interference amplitude of the Z"component. Taking the variation of the Z"component as the dependent variable， using the reference fitting formula to obtain the variation of the H"and D"components can be applied to the daily processing work of the magnetic field interference of the Shanxi geomagnetic station array under the influence of the Jinnan line high-voltage direct current transmission.

[Keywords] HVDC transmission linear; relationship influence; interference amplitude

*通訊作者："劉素珍（1980-），女，高級(jí)工程師，主要從事地震電磁學(xué)和地震監(jiān)測(cè)預(yù)報(bào)研究。E-mail：liusuzhen@sxmu.edu.cn

作者簡(jiǎn)介："范磊（1986-），男，工程師，主要從事地球物理場(chǎng)觀測(cè)與研究工作。E-mail：471557087@qq.com

基金項(xiàng)目："山西省地震局科研項(xiàng)目（SBK-2330-02）資助。