大同站FHD-2B地磁观测数据典型干扰分析_郝雪景.pdf

郝雪景，韩胜，高龙飞，等.大同站 FHD-2B 地磁观测数据典型干扰分析 J.地震科学进展,2023,53(6):251-260.doi:10.19987/j.dzkxjz.2022-157HaoXJ,HanS,GaoLF,etal.TypicalinterferenceanalysisofFHD-2BobservationdataatDatongseismicstationJ.ProgressinEarthquakeSciences,2023,53(6):251-260.doi:10.19987/j.dzkxjz.2022-157大同站 FHD-2B 地磁观测数据典型干扰分析*郝雪景1,2)韩胜1,2)高龙飞1,2)刘炜1,2)1)山西省地震局，山西太原0300212)太原大陆裂谷动力学国家野外科学观测研究站，山西太原030025摘要山西省大同地震监测中心站 FHD-2B 地磁观测数据质量主要受仪器自身观测系统和周围环境干扰影响，由于仪器出厂后性能接近稳定，因此要提高观测数据质量，可以从改善环境干扰入手。本文从连续率、完整率和背景噪声 3 个方面对 20172021 年的观测数据进行质量评价，发现自 2019 年 6 月仪器更新改造后，仪器背景噪声水平明显降低。通过分析场地施工、车辆运行、高压直流输电、断电等典型干扰事件下的数据曲线形态，发现施工干扰一般表现为各分量受干扰时段一致，干扰幅度不同；高压直流输电干扰常表现为 Z 分量干扰严重，可以通过国家地磁台网中心给出的干扰时段，结合周边地磁台站数据进行处理；车辆干扰对数据各分量均产生干扰，对应于车辆活动时段；断电干扰表现为 D、H、Z 分量观测数据曲线变粗，毛刺增多，甚至出现错误数据，通过修改仪器参数解决。总结不同干扰曲线特征，便于工作人员正确使用 FHD-2B 地磁数据进行地震研究。关键词FHD-2B；地磁；质量评价；干扰；大同中图分类号：P315.63文献标识码：A文章编号:2096-7780(2023)06-0251-10doi：10.19987/j.dzkxjz.2022-157Typical interference analysis of FHD-2B observation data at Datong seismic stationHao Xuejing1,2),Han Sheng1,2),Gao Longfei1,2),Liu Wei1,2)1)ShanxiEarthquakeAgency,ShanxiTaiyuan030021,China2)NationalContinentalRiftValleyDynamicsObservatoryofTaiyuan,ShanxiTaiyuan030025,ChinaAbstractThegeomagneticobservationdataofFHD-2BatDatongseismicstationinShanxiismainlyaffectedbytheinterferenceoftheinstrumentsownobservationsystemandthesurroundingenvironment.Sincetheperformanceofthe instrument is nearly stable after delivery，we can improve the quality of observation data by improving theenvironment.In this paper，the quality of observation data from 2017 to 2021 was evaluated from three aspects：continuityrate，integrityrateandbackgroundnoise，anditisfoundthatthebackgroundnoiseleveloftheinstrumenthasbeensignificantlyreducedsincetheinstrumentwasupdatedinJune2019.Byanalyzingthedatacurvesofdifferentinterferencessuchasconstructioninterference，vehicleinterference，high-voltagedirectcurrenttransmissioninterferenceandpowerfailure，itisfoundthattheconstructioninterferenceisgenerallymanifestedinthesameperiodoftimewhen*收稿日期：2022-10-20；采用日期：2022-11-30。基金项目：山西省地震局科研项目（SBK-2009）和中国地震局地质研究所国家野外科学观测研究站课题（NORSTY20-05）联合资助。通信作者：郝雪景（1994-），女，助理工程师，主要从事地震监测工作。E-mail：。第53卷第6期地震科学进展Vol.53No.62023年6月ProgressinEarthquakeSciencesJun.,2023eachcomponentisinterfered，andtheinterferenceamplitudeisdifferent；theinterferenceofhigh-voltagedirectcurrenttransmissionisusuallyseriousinZcomponent，whichcanbeprocessedthroughtheinterferenceperiodgivenbytheNational Geomagnetic Network Center and combined with the data of surrounding geomagnetic stations；vehicleinterference interferes with each component of data，corresponding to vehicle activity period；the power failureinterferenceshowsthattheobserveddatacurvesofD，HandZcomponentsbecomethicker，burrsincrease，andevenwrongdataappear，whichcanbesolvedbymodifyingtheparametersoftheinstrumenthomepage.Thecharacteristicsofdifferentinterferencecurvesaresummarized，whichisconvenientforworkerstouseFHD-2Bgeomagneticdatacorrectlyforseismicresearch.KeywordsFHD-2B;geomagnetism;qualityevaluation;interference;Datong 引言山西断陷盆地历史上地震频发，其主要原因是断陷盆地内的全新世活动断裂强烈活动引起的1-3。地磁场变化中存在与孕育地震有关的地壳磁场变化4，通过地磁观测可以得到地磁场变化信息，全国目前有逾百个地震台站具有地磁观测设备。许多学者对FHD 地磁观测数据进行过分析研究。彭洪军5对FHD 磁力仪使用过程中的干扰因素进行分析并制订合理的预处理措施；李庆武等6从多方面探寻降低丰宁台 FHD 仪背景噪声的方法；陈贤等7、郭灏明8和彭澎等9分析了信阳台、盐城台和宿迁台地磁观测中的各类干扰识别与处理方法；何宇飞等10根据近年来 FHD 仪观测中出现的问题对仪器进行升级改进。山西省大同地震监测中心站（下文简称大同站）地磁观测配备 1 套 FHD-2B 质子矢量磁力仪、2 套GM4-XL 磁通门磁力仪和 1 套 OVERHAUSER 磁力仪，其中 FHD-2B 磁力仪于 2008 年 1 月正式投入运行，经过 10 多年的观测，积累了大量的数据资料。本文结合国家地磁台网中心的评议结果对大同站20172021 年 FHD-2B 观测数据进行质量评价，分析影响其观测水平的因素，通过查阅观测日志、对比原始与预处理数据曲线，总结不同干扰的曲线变化特征，便于工作人员正确识别干扰并处理数据，提高观测资料质量。1台站及仪器信息1.1 台站概况大同站上皇庄观测站位于大同市马军营乡上皇庄村，距市区约 4km，地理位置如图 1 所示。台站处于大同盆地边缘两大地貌单元的分界处，台站西北部为山区，山地中沟谷发育，切割较深，东南部为洪积倾斜平原，口泉断裂在台站附近通过11。地磁观测场地远离居民住宅区、公路、采石场、电磁线路，测区地形平坦、开阔。N100 m113 13 E113 18113 2340 11 N 40 0539 59图1大同站地理位置图Fig.1LocationmapofDatongstation1.2 仪器介绍FHD-2B 质子矢量磁力仪（下文简称 FHD 仪）观测系统包括主机、探头、分量线圈及连接线路，利用质子旋进原理连续观测地磁场的总强度 F、水平分量 H、磁偏角 D。通过 F、H 可以计算得到垂直分量Z，由数据处理软件计算转换产出。该仪器采用数字化、网络化控制，工作人员通过访问仪器主页，可以下载观测数据，控制仪器参数，监视运行状态等。仪器主要技术指标为：分辨率 0.1nT，总场 F 精度 0.3nT，垂直分量 Z 精度 0.6nT，水平分量 H 精度 0.6nT，磁偏角 D 精度 0.1，测量范围 2000070000nT，采样率为 1min。仪器主机、探头及分量线圈如图 2 所示。大同站 FHD 仪于 2008 年 1 月正式运行，仪器探头和线圈独立放置在探头室的观测墩上，主机放在探头室东北向约 20m 远的地磁房内，使用沿地埋线管布设的双绞线传输。252地震科学进展2023年2质量评价结合国家地磁台网中心 FHD 观测评比结果可知，大同站主要受数据连续率、完整率和背景噪声影响，偶有数据预处理错误等。随着观测系统的不断完善、工作人员处理数据水平的不断进步，除一些不可抗拒的干扰因素外，大同站 FHD 仪观测质量在逐步改善。2.1 连续率和完整率观测数据的连续率反映原始数据的完整性，完整率反映预处理数据的完整性，二者可以直接表征仪器的运行状态及经过处理后数据的保留程度。影响连续率的主要因素是数采漏采、仪器故障、仪器标定、停电等导致的数据缺失，影响完整率的主要因素是场地施工、人为干扰等错误数据的删除。统计大同站 20172021 年的连续率和完整率（表 1）。2019 年连续率和完整率均低于 99%，原因是2019 年 6 月大同站 FHD 仪仪器更新改造，更换了新的主机、探头、线圈和信号线，期间存在数据缺测和施工干扰，并且 12 月观测场地内存在施工干扰（安装保温棚）。2020 年仪器在 4 月和 12 月出现两次主机死机故障，导致数据缺测；受 9 月探头房周边土地修整施工和 10 月办公机房改造干扰的影响导致完整率较低。2021 年仪器运行稳定，连续率和完整率较好。2.2 背景噪声地磁观测参考背景噪声是指固定台站某套仪器对磁场变化响应的灵敏程度，是数据质量评价的一项重要内容12。图 3 为 20172021 年大同站 FHD 仪各测项分量背景噪声值与国家地磁台网中心噪声平均值的对比图。00.40.81.22017201820192020202100.20.40.60.81.02017201820192020202100.20.40.60.82017201820192020202100.20.40.60.820172018201920202021D 分量噪声值/nTH 分量噪声值/nTZ 分量噪声值/nTF 分量噪声值/nT年份年份年份年份台网值大同站更新仪器更新仪器更新仪器更新仪器图3大同站与台网背景噪声值对比Fig.3ThecomparisonofbackgroundnoisebetweenDatongstationandnetworkcenter2019 年 6 月，大同站 FHD 仪进行了仪器更新改造（图 3 箭头所示），更换了新的主机、探头和线圈。可以看出，在此之前大同站背景噪声一直较大，D 分量分布在 0.71.2nT 之间，H 分量分布在 0.30.9nT表120172021 年大同站 FHD 仪数据连续率、完整率Table1The continuity rate and integrity rate of FHD datafrom2017to2021年份连续率/%完整率/%HDFHDF201799.9499.9499.9499.5299.4199.68201899.7399.7399.7399.6799.6699.67201998.4498.4498.4498.2598.2598.24202099.7999.7999.7999.2699.2699.61202199.9999.9999.9999.9999.9999.99（a）仪器主机（a）The mainframe of instrument（b）仪器探头及分量线圈（b）The probe and component coil图2FHD-2B 磁力仪结构图Fig.2ThestructurediagramofFHD-2B第6期郝雪景等：大同站 FHD-2B 地磁观测数据典型干扰分析253之间，Z 分量分布在 0.40.7nT 之间，F 分量分布在0.30.7nT 之间，多数情况下都超出台网噪声平均值，是制约观测数据质量的主要因素。6 月之后，D、Z、F 分量噪声值明显下降，H 分量噪声水平略微下降，普遍低于台网噪声平均值，其中 H、Z、F 三个分量噪声值全年较稳定，D 分量变化幅度不超过 0.2nT，很大程度改善了背景噪声超标的情况。3典型干扰分析大同站 FHD 仪主要受到场地施工、车辆运行、高压直流输电、断电等干扰。在数据处理中能够准确识别各种干扰并采取正确的预处理方式，对保证观测数据质量至关重要。3.1 场地施工干扰当仪器观测场地周围施工时，使用的电磁性工具、建筑材料等不可避免地会对地磁观测产生影响。图 4 为 2019 年 12 月 10 日施工干扰事件的数据曲线。当日大同站工作人员在 FHD 仪探头室东南侧约20m 处安装保温棚，两个施工时段都对观测数据产生干扰（图中方框所示），压制了正常曲线形态，当删除掉干扰数据后，数据正常的日变趋势得以显示。75311原始数据预处理数据281402816028180282002822028240477904781047830478504787055475554955551555535555550123456789 10 11 12 13 14 15 16 17 18 19 20 21 22 230123456789 10 11 12 13 14 15 16 17 18 19 20 21 22 230123456789 10 11 12 13 14 15 16 17 18 19 20 21 22 230123456789 10 11 12 13 14 15 16 17 18 19 20 21 22 23D 分量/H 分量/nTZ 分量/nTF 分量/nTt/ht/ht/ht/h图4大同站 FHD 仪 2019 年 12 月 10 日施工干扰数据曲线Fig.4ConstructioninterferencetoFHDatDatongstationonDecember10，2019施工干扰一般表现为各分量受干扰时段一致，干扰幅度不同。对于场地施工产生的干扰数据，多采取缺数处理，删除相应的错误数据。3.2 高压直流输电干扰随着“西电东输”工程的不断发展，高压直流输电线路也逐渐增多，当输电线路出现故障或调试时，线路中的不平衡电流会对一定范围内的地磁观测产生干扰13。地磁台网高压直流输电判别处理系统能够给出干扰线路、受影响台站、干扰时间和幅度，为工作人员数据预处理提供参考。山西省地磁观测数据受到宁绍线、宁东线、哈郑线、晋南线、上临线、昌宣线、陕湖线和锡泰线等输电线路的影响，各线路与台站的空间位置如图 5 所示。大同站地磁数据主要受到晋南线、宁东线、上临线和锡泰线 4 条输电线路干扰，宁东线、上临线和锡泰线最大干扰幅度均不超过 1.5nT，晋南线最大干扰幅度约 3.3nT。这是由于大同站距离晋南线送端换流站约 110km，受其干扰较其他线路严重，且晋南线对山西省内地磁观测台站均有影响，故本文主要分析晋南线干扰情况。晋南线输电线路以山西朔州为254地震科学进展2023年送端，江苏盱眙为受端，途径山西、河北、山东、河南、安徽和江苏 6 省，于 2016 年 4 月投入运行，额定功率 8000MW，电流 5000A，电压800kV，线路长度1118.5km，避让距离 20km。2019 年 4 月 28 日晋南线输电线路 17:4320:35发生干扰事件，选取当日 1524 时大同站和太原站FHD 仪观测数据曲线进行对比（图 6），干扰时段为图中方框所示。可以看出，大同站 FHD 仪 Z 分量受上海庙昌宣线哈郑线宁绍线宁东线上临线锡泰线晋南线陕湖线郑州哈密盱眙朔州青岛绍兴灵武武汉陕北临沂泰州宣城昌吉锡林浩特承德北京固原榆林蒙城南京昔阳代县大同定襄临汾兰州太原嘉峪关德令哈乌鲁木齐12011010090 E45 N403530250300600 km150输电线路换流站台站省界山西省图5山西省受高压直流输电干扰线路图Fig.5LinediagramofShanxiProvinceinterferedbyHVDCtransmissionD 分量/2818028190282002821047808478124781647820478245549855503555085551355518151617181920212223151617181920212223151617181920212223151617181920212223H 分量/nTZ 分量/nTF 分量/nTD 分量/H 分量/nTZ 分量/nTF 分量/nT1.40.90.40.10.61.1原始数据预处理数据t/h（a）大同站（a）Datong station（b）太原站（b）Taiyuan station151617181920212223151617181920212223151617181920212223151617181920212223t/ht/ht/ht/ht/ht/ht/h1.40.90.40.10.61.129530295402955029560453504535545360453654537045375541205412554130541355414054145图62019 年 4 月 28 日大同站、太原站 FHD 仪受晋南线输电干扰曲线Fig.6JinnantransmissioninterferencetoFHDatDatongstationonApril28，2019第6期郝雪景等：大同站 FHD-2B 地磁观测数据典型干扰分析255干扰，干扰幅度为 2.1nT；太原站 FHD 仪各分量均受干扰，Z、F 分量较严重，最大干扰幅度为 4.8nT，这可能与台站和输电线路的相对位置有关。大同站Z 分量受干扰台阶方向向下，太原站 Z 分量台阶向上，这是由于大同站和太原站位于输电线路两侧，导致 Z 分量干扰变化方向相反。ZIrs蒋延林等14推导了高压直流输电干扰时产生的附加磁场值（nT）与线路中不平衡电流（A）、台站与线路视距离（km）的关系式：Z=0.2Irs(1)由上式可知，对于同一条输电线路的同一次干扰，不平衡电流一致，由于台站和输电线路位置固Zrs定，那么两个台站的比值就等于 的反比，也就是一个恒定值，即关系式：Z1Z2=0.2Irs10.2Irs2=rs2rs1(2)因此，当知道一个台站的干扰幅度后，可以利用上式计算出另一个台站的干扰值。统计 20182020 年大同站、太原站、代县站和定襄站 4 个台站受晋南线输电干扰事件，以大同站Z 分量干扰幅度为基准，将其他 3 个台站的数据与大同站做比值，结果如表 2 所示，将大同站干扰幅度值作为 X 轴，其他台站干扰幅度作为 Y 轴，绘制干扰拟表220182020 年晋南线干扰统计表Table2StatisticsofJinnantransmissioninterferencefrom2018to2020干扰时段大同太原代县定襄幅度/nT幅度/nT比值幅度/nT比值幅度/nT比值2018-02-2116:1316:180.71.62.28617.424.857912.8572018-02-2116:1916:210.71.62.28617.424.857912.8572018-06-0912:1012:120.51.22.40013.126.2006.813.6002018-06-0912:3112:330.51.22.40013.126.2006.813.6002018-06-1507:2907:311.842.22243.424.11122.512.5002018-06-1511:2611:291.842.22243.424.11122.512.5002019-03-2516:3416:362.14.72.23851.324.42926.612.6672019-03-2519:0219:052.14.72.23851.324.42926.612.6672019-04-2817:4317:542.14.82.28652.324.90527.212.9522019-04-2820:3220:352.14.82.28652.424.95227.212.9522019-05-2217:5318:082.14.82.28652.124.81027.112.9052019-05-2218:4518:482.14.82.28652.224.85727.112.9052019-05-2917:2617:292.14.82.28652.124.81027.112.9052019-05-2919:4019:432.14.82.28652.224.85727.112.9052019-05-2920:4320:571.43.22.28634.824.85718.112.9292019-05-2920:5720:590.71.62.2861724.2868.812.5712019-05-2921:1821:202.14.82.28651.824.66726.912.8102020-03-0906:2606:362.14.72.23850.724.14326.312.5242020-03-0908:2108:242.14.72.23850.724.14326.312.5242020-05-2311:5411:591.43.12.21433.824.14317.512.5002020-05-2312:4412:461.43.12.21433.824.14317.512.5002020-06-0914:5114:582.14.82.28651.924.71426.912.8102020-06-0914:5915:012.14.82.28651.924.71426.912.8102020-06-0915:0115:062.14.82.28652.324.90527.112.9052020-06-0915:0715:092.14.82.28652.324.90527.112.905256地震科学进展2023年合曲线，如图 7 所示。由于代县站距离晋南线最近，受干扰的幅度最大，大同站和代县站受干扰方向一致，太原站和定襄站一致。经过分析可以看出，各台站干扰幅度比值拟合效果较好，都接近某一固定值，由于不可避免的受到观测噪声影响，比值也会出现波动。太原站与大同站比值近似为2.261，代县站与大同站比值近似为 24.56，定襄站与大同站比值近似为12.743。60402002040602.5 2.0 1.5 1.0 0.500.51.01.52.02.5太原台代县台定襄台其他台站 Z 分量干扰幅度/nT大同站 Z 分量干扰幅度/nTK=24.56K=2.261K=12.743图7太原站、代县站、定襄站和大同站干扰幅度拟合曲线Fig.7Fitting curve of interference amplitude of Taiyuan，Daixian，DingxiangandDatong高压直流输电干扰的原则是干扰上升总量等于下降总量，数据曲线常表现为尖峰、缓变、台阶，Z 分量受干扰幅度明显大于其他分量，且各台站同时出现干扰，利用地磁台网中心给出的干扰时段，结合周边地磁台站数据，可以有效识别高压直流输电干扰。对于高压直流输电干扰，可以根据曲线形态，结合台阶改正、去除缓变、剔除尖峰进行数据处理。3.3 车辆干扰当有车辆在地磁观测场地附近活动时，往往会对数据造成干扰。图 8 是 2020 年 9 月 20 日大同站FHD 仪数据受车辆干扰曲线。经现场工作人员记录，当日（世界时）06:3007:30 和 10:1510:45 有重达 10t 的大型金属卡车两次进入地磁观测场地，由当日地磁数据曲线可以看出，两次车辆活动对FHD 仪各分量数据都产生了干扰。车辆 06:30 第一次进入场地时，车辆驶入后停靠在探头房东北侧约50m 处，停放 1 小时未活动后驶出场地，对应数据曲线是台阶型干扰，D 分量最大变化 2.3，H 分量变化 13nT，Z 分量变化 4.5nT，F 分量变化 3.5nT。当车辆 10:15 第二次驶入场地时，车辆围绕场地行驶一圈，中途不时停靠，距探头房最近时仅 2m，D 分量最大变化 13.8，H 分量变化 130.9nT，Z 分量变化70.8nT，F 分量变化 62.4nT，曲线形态复杂，突跳和台阶混合，干扰幅度也明显大于第一次车辆活动。当地磁房附近有车辆过往，一般会引起数据突跳，有车辆停放时会产生台阶，干扰在观测数据各19149416原始数据预处理数据D 分量/H 分量/nTZ 分量/nTF 分量/nTt/ht/ht/ht/h28100281502820028250283002835047790478204785047880479104794055520555505558055610556400123456789 10 11 12 13 14 15 16 17 18 19 20 21 22 230123456789 10 11 12 13 14 15 16 17 18 19 20 21 22 230123456789 10 11 12 13 14 15 16 17 18 19 20 21 22 230123456789 10 11 12 13 14 15 16 17 18 19 20 21 22 23车辆驶入车辆驶入车辆驶入车辆驶入车辆驶出车辆驶出车辆驶出车辆驶出车辆活动车辆活动车辆活动车辆活动图82020 年 9 月 20 日大同站 FHD 仪受车辆干扰曲线Fig.8VehicleinterferencetoFHDatDatongstationonSeptember20，2020第6期郝雪景等：大同站 FHD-2B 地磁观测数据典型干扰分析257分量上均有体现，受干扰时段对应于车辆活动时段。工作人员可以结合台站监控设备对观测场地周边环境进行实时监控，利于识别车辆干扰。对于车辆停靠引起的台阶型干扰，可利用台阶改正和剔除尖峰进行处理；当车辆通过引起突跳时可利用剔除尖峰处理；当干扰复杂，在无法保证曲线形态的情况下，可采用缺数处理。3.4 断电干扰2020 年 10 月大同站进行办公机房改造，期间数次断电，恢复供电后发现 FHD 仪观测数据出现了不规律干扰，且每次断电后都会出现这种干扰情况，图 9 为 10 月 13 日和 14 日的两次断电干扰数据曲线，断电时数据缺测，随之 D、H、Z 分量观测曲线变粗，出现了错误数据，而 F 分量正常。经工作人员检查，观测环境正常，仪器接线完好，最后发现 FHD仪主页的仪器参数错误，实时参数与原始参数出现较大偏差，调谐值偏离超过 100nT。当仪器正常工作时，仪器的实时参数与原始参数状态较为接近，调谐值最大偏离约 10nT，仪器参数界面如图 10 所示。原始参数是仪器在初始仪器调试或标定时形成的工作参数，实时工作参数是仪器当前的工作参数。当调谐值严重偏离时，会影响到观测数据的精度，使观测数据曲线线条变粗，毛刺增多，甚至出现乱数现象，仪器工作参数不能实时跟踪变化，在选择“恢复缺省参数”功能后，实时参数恢复为原始参数，仪器正常工作，数据曲线也恢复正常。30201001020原始数据预处理数据D 分量/H 分量/nTZ 分量/nTF 分量/nT281002830028500287002890047500476004770047800479004800048100555455555555565555755558502468 10 12 14 16 18 20 22 02468 10 12 14 16 18 20 2210-1310-1402468 10 12 14 16 18 20 22 02468 10 12 14 16 18 20 2210-1310-1402468 10 12 14 16 18 20 22 02468 10 12 14 16 18 20 2210-1310-1402468 10 12 14 16 18 20 22 02468 10 12 14 16 18 20 2210-1310-14t/ht/ht/ht/h图92020 年 10 月 1314 日大同站 FHD 仪断电干扰曲线Fig.9PoweroffinterferencetoFHDatDatongstationfromOctober13to14，2020对于此种断电导致实时参数变化的干扰，需要工作人员监控通电情况，在恢复供电后及时检查仪器主页参数情况，尽快处理，保证数据准确性。4结论（1）制约大同站 FHD 仪观测水平的主要因素是背景噪声超标，在 2019 年 6 月仪器更新改造后得到明显改善。（2）通过分析场地施工、车辆运行、高压直流输电、断电等典型干扰事件下的数据曲线形态，总结出以下特征：施工干扰会产生不规律的错误数据，各分量受干扰时段一致，干扰幅度不同；车辆干扰表现为观测场地附近有车辆过往时引起数据突跳，有车辆停放时产生台阶，数据各分量均受干扰，对应于258地震科学进展2023年车辆活动时段；高压直流输电干扰常表现为 Z 分量干扰严重，且各台站干扰具有同步性，可以通过国家地磁台网中心给出的干扰时段，结合周边地磁台站数据进行处理；断电干扰表现为 D、H、Z 分量观测数据曲线变粗，毛刺增多，甚至出现错误数据，可以通过修改仪器参数解决。当数据曲线可以通过台阶改正、去除缓变和剔除尖峰处理时，应尽量保证数据的完整性；当干扰复杂无规律时，才采用缺数处理。（3）工作人员在进行数据预处理时，应结合当天干扰情况和数据形态，正确进行数据处理。同时，对于无法判定的数据异常变化，要结合其他地磁台站数据对比分析，利用相对差值检测等方式及时发现数据变化原因，进行合理的预处理，保证数据可靠性。致谢本文地磁观测数据质量评价是基于国家地磁台网中心的评议结果，高压直流输电干扰时段也由地磁台网中心给出，在此表示诚挚的谢意。参考文献闫小兵，周永胜，李自红，等.1695 年临汾 7级地震发震构造研究 J.地震地质，2018，40(4)：883-902YanXB，ZhouYS，LiZH，etal.AstudyontheseismogenicstructureofLinfenM7earthquakein1695J.SeismologyandGeology，2018，40(4)：883-9021闫小兵，周永胜，李自红，等.山西浮山断裂的晚第四纪活动与位移速率 J.地震地质，2022，44(1)：35-45YanXB，ZhouYS，LiZH，etal.TheLateQuaternaryactivityanddisplacementrateofFushanfaultinShanxiJ.SeismologyandGeology，2022，44(1)：35-452曾金艳，李自红，陈文，等.运城盆地盐湖南岸断层晚第四纪活动特征研究 J.第四纪研究，2020，40(1)：124-131ZengJY，LiZH，ChenW，etal.StudyontheactivitycharacteristicsofthesouthbankfaultofYunchengsaltlakeinYunchengbasinsincetheLateQuaternaryJ.QuaternarySciences，2020，40(1)：124-1313郭明瑞，胡久常，王锡娇，等.剑河 M5.5 地震前地磁异常特征分析 J.地震科学进展，2020，50(5)：12-17GuoMR，HuJC，WangXJ，etal.AnalysisofgeomagneticanomaliesbeforetheJianheM5.5earthquakeJ.ProgressinEarthquakeSciences，2020，50(5)：12-174首页 技术指标 仪器参数 数据下载 当天数据 仪器状态 帐户管理 十五规程仪器参数测量参数参数名称探头参数(L)分布电容(C)水平补偿电流(AH)偏置补偿电流(AE)总场调谐值(F)水平磁场调谐值(H)正向合成磁场(R+)反向合成磁场(R)偏置磁场(C)F 修正值(F)H 修正值(H)R+修正值(R+)R 修正值(R)D 修正值(D)时间原始参数33.3007.144.52090.0855617.028114.735154.535101.821060.9+00.0+00.8+00.0+00.00002.042022-09-30 19:0433.3007.144.52090.0855621.428113.835155.135112.321070.3+00.0+00.8+00.0+00.00002.042022-11-26 21:20实时参数恢复缺省参数图10大同站 FHD 仪仪器参数界面Fig.10InstrumentparameterinterfaceofFHDinDatongstation第6期郝雪景等：大同站 FHD-2B 地磁观测数据典型干扰分析259彭洪军.FHD 质子矢量磁力仪资料预处理分析 J.科技与创新，2019(18)：145-146PengHJ.DatapreprocessinganalysisofFHDprotonvectormagnetometerJ.ScienceandTechnology&Innovation，2019(18)：145-1465李庆武，胡秀娟，张玉林，等.丰宁台 FHD 仪背景噪声减小措施分析 J.山西地震，2019(4)：48-51LiQW，HuXJ，ZhangYL，etal.AnalysisofbackgroundnoisereductionmeasuresofFHDinstrumentinFengningStationJ.EarthquakeResearchinShanxi，2019(4)：48-516陈贤，黄恩贤，成万里，等.信阳台 FHD 地磁仪观测数据典型干扰识别及数据处理 J.地下水，2020，42(2)：101-103ChenX，HuangEX，ChengWL，etal.TypicalinterferenceidentificationanddataprocessingofobservationdataofFHDmagneto-meterinXinyangstationJ.Ground