董 博，纪春玲，张环曦，周安聘，茅远哲，唐存智，张尚清

(1. 河北省地震局，石家庄 050021； 2. 河北地质大学， 石家庄 050031)

河北地区近年地磁场变化特征

董 博1,2，纪春玲1*，张环曦1，周安聘1，茅远哲1,2，唐存智2，张尚清2

(1. 河北省地震局，石家庄 050021； 2. 河北地质大学， 石家庄 050031)

通过分析河北及周边10个地磁台2008—2016年的各分量年均值数据和年变率曲线变化形态，发现3个分量均呈单调变化，Z分量和F总场年变率从2008年开始均呈现出上升变化态势；利用克里金插值法对所研究区域地磁场变化特征进行插值计算，分析结果认为：在该研究区域内随着时间的递增，3个分量等值线均逐渐由北向南缓慢的迁移。该研究对于了解近年来河北及附近区域地磁场变化规律有积极指导意义。

地磁场；年变率；克里金插值法；河北

0 引言

地磁学是地球物理学的重要组成部分，地球磁场和地球引力场一样，属地球物理场。按场源位置划分，地磁场由外源场和内源场2部分组成。外源场起源于地表以上的空间电流体系，内源场包括源于地核的基本磁场，源于地壳磁性岩石的地壳场，以及外部变化磁场在地球内部的感应场。地磁场的变化和分布具有空间上的相关性和时间上的延拓性[1]，地磁台站观测的是包括占百分之九十几的内源场与百分之几的外源场在内的整个地磁场。通常来说，磁异常场属于内源场的一部分，也称为岩石圈磁场或地壳磁异常场。而在地磁异常中，有地震引起的地磁异常又占很小的比重，地磁场有各类长期变化和短期变化，只有研究并充分认识他们的规律，才能正确提取异常场[2-4]。

本文所研究的原始数据取自中国地震局地磁台网数据库，通过分析研究河北及周边地区2008—2016年9年间该地区地磁场变化规律，认为台站在所研究区域内分布较合理，台站记录的磁场变化特征基本代表了河北省及附近区域地磁场的时空变化。不仅对了解地磁场变化形态有意义，还有助于识别地磁异常信号，为提供更有价值的地震地磁异常信息做好准备工作。

1 台站概况

河北省现有10个地磁FHD观测台站，考虑到部分台站FHD仪器正式观测时间较晚，且存在大量的数据缺测，故只选取省内资料连续性较好的6个台站数据作为基础研究资料。为了能够更加真实地反映河北地区的地磁场变化特征及避免边界效应的产生，又选取靠近河北边界的周边省份4个地磁台数据作为研究资料(表1)。

表1 台站参数

图1 地磁台站位置分布图

2 资料收集与处理

从中国地震局地磁台网数据库中提取了2008—2016年广平等10个地磁台站各要素(F、H、Z)的原始数据。通过对这10个监测台站9年的地磁资料进行分析，认为数据真实性符合我们对该地区地磁场变化特征分析研究的要求。通过对原始数据进行预处理，突出了地磁数据的趋势变化，并计算得到了河北省及周边地磁台站各要素的年均值。

3 数据分析

3.1 各台数据变化趋势分析

本文所选各地磁台9年的观测资料客观地记录和反映了该地区的地磁场的近期变化特征。对这些台2008—2016年的地磁总场F，垂直分量Z及水平分量H观测资料进行统计、计算后，用年均值做图(图2～4)，可从F、H、Z各自的变化趋势看出，形态变化虽有所不同，但一致性较好；3个分量曲线形态均呈现单调变化，比如各台H分量从2008年开始呈逐年缓慢减小，变化曲线单调递减，近似于线性变化，到2016年约下降了309.2 nT，平均每年下降38.7 nT；Z分量和F总场曲线变化形态则相反，二者均呈现逐年缓慢单调上升现象，各台Z分量和F总场从2008—2016年上升的平均幅度分别为449.8 nT和231.1 nT，平均年变化幅度为56.2 nT和28.9 nT。

表2 各地磁台站水平分量H年均值统计表 nT

图2 2008—2016年各台H分量年均值曲线图

表3 各地磁台站垂直分量Z年均值统计表 nT

图3 2008—2016年各台Z分量年均值曲线图

nT

图4 2008—2016年各台F总场年均值曲线图

3.2 各台站年变率分析

地磁数据分析研究中，通常用1年的长期变化率来表示这一年地磁要素的变化大小[5]，t0年的长期变化率(年变率)的定义是

(1)

其中：A2与A1是t2与t1年某地磁要素的年均值，t0=t2+1/2Δt。为了便于研究，将Δt的时间间隔取为1年，于是公式简化为α=ΔA。令t0=t2。所以本文的年变率是指某年的年均值与前一年年均值之差，通过这样处理，年变率的变化趋势会显得更加突出和直观。图5～7分别给出了所选各台记录的地磁总场F，垂直分量Z和水平分量H的年变率曲线图。

由图5可以看出，各台H分量的年变率从2008年开始整体呈现缓慢下降趋势，至2012年降至最低，而后发生转折，2012—2014年开始变化为缓慢上升，而后在2015年年变率又开始变为下降，之后2016年开始转折上升。由于所选资料时间有限，至于近年H分量年变率变化形态是否具有短周期性，仍需后期长时间数据进行分析验证。由图6可以看出，各台Z分量年变率呈现逐渐上升的趋势，变化速度逐渐减小，至2014和2015年年变率变化速度趋近于零，2016年开始呈现下降趋势。由图7可以看出，总场F的年变率在2008—2011年变化速度较缓慢，但仍然呈现上升态势，从2011年之后年变率变化速度开始迅速增加，至2013年左右开始变速减小接近于0，至2016年又开始缓慢上升。从图6～7可以看出，Z分量和F总场年变率从2008年开始均呈现上升变化态势，这与张小涛等[6]在分析前20年河北省内的绝对观测地磁数据的年变率时所呈现出的正弦形态变化相符合。按照其分析结果来看，2007后年变率变化曲线应该呈现出上升的态势，且上升时间大概为7 a左右。从年变率曲线能够看出一个共同的特点，在研究区域内地磁场年变化速率南部大于北部，东部大于西部，有自东南向西北逐渐减弱的规律,这与谢用明等[7]的研究成果基本吻合。

图5 各台H分量年均值年变率曲线图

图6 各台Z分量年均值年变率曲线图

图7 各台F总场年均值年变率曲线图

3.3 利用插值法分析区域地磁场变化特征

地磁数据的观测和采样往往在不规则且离散的测网(或测线)进行。为合理并正确解释地球物理异常，在保持原有场的特征前提下，有必要进行均匀且较密集规则网格的数据插值，形成直观的二维等值线，为进一步异常分析提供依据。

常用格网化插值方法包括径向基函数法，克里金法，最小曲率法，三角网线性内插法等。局部地磁场中的地球主磁场，地磁异常场及扰动磁场与克里金法的3种成分假设具有一一对应关系，因此克里金法能较为准确地反映局部地磁场的分布空间相关特性规律，可用于局部磁场空间插值[8]。

在插值进行前，由于所选10个地磁台站海拔高程相差较大，因此需对地磁观测数据进行高度改正。具体思路为，选定一个高度作为标准高度，选择相应经纬度作为计算点得出该地正常场，计算不同高度修正值。若台站海拔高于标准点，则实测值需要加上改正值，反之，减去改正值[9]。应用克里金法对研究区域2008—2016年地磁资料进行插值计算，插值间隔为0.075×0.075,使得x分成96个网格，而y分为100个网格，得到地磁场各分量的二维等值线和三维等值线。因2008—2016年各分量变化趋势差别不大，在此各分量三维等值线图只给出中间年份2011年的图作为参考。

图8(1)～(9)：2008—2016年H分量二维等值线图；图8(10)～(18)：2008—2016年Z分量二维等值线图；图8(19)～(27)：2008—2016年F分量二维等值线图 图8 2008—2016年河北及周边H、Z、F分量二维等值线图

图9 2011年河北及周边H、Z、F分量三维等值线图

由图8(1)～(9)和图9a可以看出，随着纬度的的增大，H值逐渐减小。在(113.5°～118°E，37°～40°N)和(118.5°～120.5°E，40°～40.8°N)范围内变化率较大，其他区域H值变化率较小；(114.5°E，35.7°N)处的H值较周围区域都大，为该研究区域的极大值。2008—2016年该研究区域H分量变化趋势较平稳，如果动态地观察这9 a的等值线图，可以形象地看出在该研究区域内，H分量等值线随着时间的递增，逐渐由北向南缓慢地迁移。

由图8(10)～(18)和图9b可知，随着纬度的增大，Z值逐渐增大。在(113.5°～118°E，36.5°～40.1°N)和(118.2°～120.5°E，40°～40.8°N)范围内变化率较大，其他区域Z值变化率较小；(115.5°E，35.7°N)处的Z值较周围区域都小，为该研究区域的极小值。2008—2016年该研究区域Z分量变化趋势较平稳，如果动态地观察这9年的等值线图，亦可以看出在该研究区域内，Z分量等值线随着时间的递增，逐渐由北向南缓慢地迁移。

由图8(19)～(27)和图9c可知，随着纬度的增大，F值逐渐增大。在(113.5°～118°E，37°～40°N)和(117.5°～120.5°E，40.1°～40.7°N)范围内变化率较大，其他区域F值变化率较小；(116.5°E，35.6°N)处的F值较周围区域都小，为该研究区域的极小值。2008—2016年该研究区域F总场变化趋势较平稳，且与Z分量的变化趋势相一致。总体来看，F总场等值线也是随时间逐渐由北向南缓慢地迁移。

4 结论与讨论

综上资料分析结果，河北地区近年地磁场变化具有如下特征：

1)从F，H，Z各自的变化趋势来看，形态变化虽有所不同，但一致性较好，且3个分量曲线形态均呈现单调变化，各台H分量逐年单调递减，近似于线性变化；Z分量和F总场曲线变化形态则相反，二者均呈现逐年单调上升现象；F、H、Z3个分量各台逐年平均变化量分别为27.6 nT、39.8 nT和56.2 nT，总体变化速度Flt;Hlt;Z。

2)Z分量和F总场年变率从2008年开始均呈现出上升变化态势，这与张小涛等[5]在分析前20 a河北省内的地磁数据的年变率时所呈现出的正弦形态变化相符合。按照其分析结果来看，2007年后，年变率变化曲线应该呈现出上升的态势，且上升时间大概为7 a左右。H分量年变率是否存在短周期变化需后期积累大量的测量数据进行验证。

3)随着纬度的增大，H值逐渐减小，F和Z值逐渐增大。3个分量在(113.3°～117°E，36.5°～39.7°N)和(118.5°～120.5°E，40°～40.8°N)范围内变化率均较大，其他区域变化率相对较小；(114.5°～116.5°E，35.6°～35.7°N)小区域范围内均为F、H、Z极值所在之处，H为极大值，F、Z为极小值。动态地观察这9 a的等值线图，可以形象地看出，在该研究区域内，3个分量等值线均逐渐由北向南缓慢地迁移。

4)F、H、Z等值线分布或许与该区内台站分布情况有关，由于研究区内分布台站较少等原因，分析精度和深度有待加强，积累更多资料后进一步分析和研究。

[1] 薛革, 陶九庆, 张继红. 山东地区地磁背景场的变化特征[J]. 高原地震, 1995, 7(4): 50-55.

[2] 张建国, 乔子云, 吕凤章, 等. 河北省地磁场变化规律研究[J]. 地震地磁观测与研究, 2007, 28(3): 8-12.

[3] 张建国, 姚丽, 刘晓灿, 等. 地震电离层VLF电磁场频谱特征研究[J]. 大地测量与地球动力学, 2012, 32(3): 110-115.

[4] 张建国, 焦立果, 刘晓灿, 等. 汶川MS8.0级地震前后ULF电磁辐射频谱特征研究[J]. 地球物理学报, 2013, 56(4): 1253-1261.

[5] 北京大学, 中国科学技术大学地球物理教研室. 地磁学教程[M]. 北京: 地震出版社, 1983: 56.

[6] 张小涛, 王莉森, 韩丽萍. 河北省近年内地磁场变化规律的初步分析[J]. 地震地磁观测与研究, 2008, 29(6): 17-21.

[7] 解用明. 河北省近年地磁Z份量长期变特征[J]. 山西地震, 1996(2): 47-50.

[8] 刘强, 陶钧, 刘旭, 等. 基于克里金插值的磁法数据格网化研究[J]. 河南科学, 2013, 31(7): 1039-1044.

[9] 董博, 纪春玲, 张环曦, 等. 采用克里金插值法分析河北地区地磁场变化特征[J]. 地震地磁观测与研究, 2016, 37(3): 147-154.

GeomagneticFieldVarivationCharacteristicsinHebeiAreainRecentYears

DONG Bo1,2, JI Chun-ling1*, ZHANG Huan-xi1, ZHOU An-pin1,MAO Yuan-zhe1,2,TANG Cun-zhi2, ZHANG Shang-qing2

(1. Hebei Earthquake Agency, Shijiazhuang 050021, China; 2. Hebei GEO University, Shijiazhuang 050031, China)

In the paper, the annual mean data and the annual variation curve of the 10 geomagnetic stations in Hebei and the surrounding areas from 2008 to 2016 are analyzed. The analysis shows that: the three components show monotonic change; the total annual change rate of the Z component and F component from the beginning of 2008 show an upward trend. Then, the characteristics of magnetic field are studied using Kriging interpolation and the results indicated that: in the study area, the three components of the isoline are slowly moving from north to south as time goes on. The study has great significance to understand the geomagnetic field variation in Hebei and its adjacent areas in recent years.

Hebei; geomagnetic field; annual variation rate; Kriging

董博，纪春玲，张环曦，等. 河北地区近年地磁场变化特征[J]．华北地震科学，2017，35(4):22-28.

2017-03-10

2017年河北省地震局地震科技星火计划项目(DZ20170510065)

董 博(1986—)，男，河北鹿泉人，工程师，主要从事地震监测与分析预报工作.E-mail：dongbo0002@126.com

*

纪春玲(1987—)，女，河北晋州人，助理工程师，主要从事地震监测与分析预报工作.E-mail：279694644@qq.com

P315.7

A

1003-1375(2017)04-0022-07

10.3969/j.issn.1003-1375.2017.04.004