0 引言

地磁场是地球固有的基本特征，由地核场、地壳磁场、变化磁场和感应磁场四部分构成(傅承义等，1985；徐文耀，2009).地磁学是以观测为基础的科学，地磁基础观测数据直接或间接用于研究地球基本磁场及其长期变化(冯丽丽，2015；张素琴等，2016)，用于监测地球外部空间电磁环境(Li et al.，2009；赵旭东等，2010；王源等，2010)，用于震磁关系研究和地震预报并为其他测量手段(如航海、地球物理勘探等)提供通化标准(Gu et al.，2006；陈斌等，2012；李琪等，2015).因此，地磁观测数据不仅具有学术研究意义，也是不可或缺的国家基础战略资源(徐文耀，2003).

截止2014年，我国已基本建成了由160个观测台站组成的覆盖全国数字化地磁前兆观测网络，实现了观测数据的准实时传输和自动汇集，并积累了大量的高精度观测数据.然而，随着城市化进程的加快，铁路、公路、电网等国家基础设施的大规模建设和投入使用，地磁台站观测受到了来自上述干扰源不同程度的影响，表现为各种干扰和数据缺失.数据干扰及缺失的情况严重影响了数字化地磁台站观测数据质量，制约地磁观测工作的开展，同时也影响依赖地磁观测数据的相关科学研究.因此，在数据已经缺失或严重受干扰的情况下尽可能不失真地重构出观测数据，保证数据的可用性对于地磁数据分析而言具有重大意义.

Yao(1988)曾将订正法应用在绝对观测数据的填补研究中，即利用wild公式将多年对应缺数月份的均值作为订正参考值，用于缺失数据大小的衡量.但是对于地磁相对记录而言，其变化主要受空间电流体系影响，磁场活动水平不同，相对记录数据差异较大，无法将多年的数据用于缺失数据填补，因此，订正法不适用于相对记录数据.朱兆才(1989)利用地磁观测数据在空间上较高的相关性，将相邻单台数据与缺失台站数据进行回归分析，提出基于回归方程利用邻台数据进行缺失数据填补的想法.Lockwood等(2013)在重现过去167年地磁活动性与近地行星际关系的研究中，用简单的最小二乘拟合将附近单台数据用于缺失数据的填补，取得了较好的结果.因此，邻台重构是值得尝试的方法.但是该方法对邻台数据依赖性较大，数据质量限制以及日变形态差异等因素必将对数据重构结果产生较大影响，重构精度无法保证.姚休义等(2016)选用了相邻多个台站的数据利用空间加权法对受干扰及缺失数据进行重构，在静日及通日重构中取得了较好的结果，但在强磁扰时段重构精度较差.人工神经网络(Artificial Neural Network，简称ANN)是20世纪80年代中期兴起的一门科学，用来模拟人脑的一些基本特征，利用非线性映射和并行处理方法，通过不断学习和训练，完成输入与输出之间的映射关系，在模式识别、机器学习、决策支持和数据挖掘等领域都得到广泛的应用，也有不少学者将其用于太阳活动和地磁活动预测(Hernandez et al.，1993；Lundstedt and Wintoft，1993；Kugblenu et al.，1999).本文将利用人工神经网络进行受干扰及缺失的地磁相对记录数据重构.

2.需要充分落实工作预案。在对突然发生的事件进行处理时，首先，相关工作人员需要了解清楚事件发生的基本情况，找出突发事件发生的真正原因，然后根据事件发生的时间、地点等对该事件的性质做出明确界定。其次，现场处理人员应当控制好自身的态度和情绪，时刻保持自身头脑清醒，才能够在事件现场做出更为准确高效的判断。

1 数据重构原理与方法

1.1 BP神经网络算法

Rumelhart和McClelland(1986)在研究并行分布式信息处理方法时，提出了利用误差反向传播训练算法构造的一种神经网络Back-Propagation Network，简称BP神经网络(Hecht-Nielsen，1989)，是目前应用最广泛的神经网络之一(沈学利和张纪锁，2010；李翔和朱全银，2013).目前已广泛应用于地学研究中的模式识别和地震预报(蔡煜东等，1993；席道瑛和张涛，1996；董晓娜等，2012；任成冕和熊晶，2012；熊晶等，2013).

研究表明，当两台站间距离小于等于1000 km时，地磁记录有较高的空间相关性，相关系数均大于0.9，随着距离增大，相关性降低(姚休义等，2016).此外，考虑地磁日变的地方时效应(冯志生等，2005)，本文在数据重构过程中以受干扰台站或缺数台站(目标台)为圆心，750 km为半径的范围内选取参考台作为输入，通过BP神经网络的预测值输出作为目标台重构数据.BP神经网络的拓扑结构如图1所示.

图1 BP神经网络拓扑结构图 Fig.1 Topological structure diagram of BP network

图中X1，X2，…，Xn是BP神经网络的输入值，Y1，Y2，…，Ym是BP神经网络的输出值(预测值)，ωij和ωjk为神经网络权值.BP神经网络可以看成一个非线性函数(Chen，1990)，网络输入值和预测值分别为该函数的自变量和因变量.当输入节点数为n、输出节点数为m时，BP神经网络表达为从n个自变量到m个因变量的函数映射关系.

根据输入变量X、输入层和隐含层之间的连接权值ωij以及隐含层阈值a，根据公式(1)计算隐含层输出H：

总而言之，我国国库集中支付制度推广了十几年，这项制度对高校财务管理造成了较大影响，有利于高校财务管理体系的构建。从当前情况来看，很多高校在财务管理系统方面仍不够完善，在未来的发展中，高校领导要转变观念，充分认识到国库集中支付制度的优点、注重高校内部预算管理的科学化、优化财务管理组织机构、建立完善的财务信息管理系统、拓宽办学资金来源，以高校国库集中支付给财务管理制度带来的变化为出发点，强化自身财务管理体系的构建，实现高校的正常运行和教育目的。

) j=1,2,…,l

(1)

式中，l为隐含层节点数，本文数据重构中取为5，f为隐含层响应函数，本文选则S形响应函数(公式2).根据隐含层输出H，连接权值ωjk和阈值b，计算BP神经网络的预测输出O.

(2)

(3)

根据网络输出O和期望输出Y，计算网络预测误差e，误差函数为

ek=Yk-Ok k=1,2,…,m

(4)

根据网络预测误差e更新网络连接权值ωij和ωjk：

为了计算BP神经网络重构法的重构精度并检验该方法的适用性，本文进行了数据仿真.以数据质量较好的河北红山台(LYH)为目标台，假设该台一整天数据缺失，以距其750 km为半径的范围选取参考台(图2).考虑观测数据质量及相位延迟互补，选取了泰安(TAA)、榆林(YUL)、蒙城(MCH)、乾陵(QIX)四台为参考台进行数据重构，参与重构数据为各台站磁通门秒数据，采样率1 sps，分辨力0.1 nT.各台与LYH的距离分别为253.1 km、454.4 km、478.5 km、662.7 km.选用的BP神经网络类型为4-5-1，即输入层有4个节点，隐含层有5个节点，输出层有1个节点.输入层节点传递函数选用对数S型传递函数(logsig)，隐含层节点传递函数选用线性传递函数(purelin)，训练函数选用动量反传的梯度下降BP算法训练(traingdm)，学习速率为0.9，每次训练的最大迭代次数为1000次.BP神经网络训练时选用缺数日期前一天(或后一天)各台站数据进行训练.

j=1,2,…,l

(5)

ωjk= ωjk+ηHj·ek j=1,2,…,l;

钙硫比的增加对所获得的粉煤灰脱硫活性起到促进作用，但CaCO3的增加必然对炉内煤粉燃烧特性产生影响，若其加入严重影响煤粉燃烧特性，那么对上述研究就失去了意义。为此，笔者进一步针对钙硫比对煤粉燃烧特性的影响进行研究。

k=1,2,…,m

生态文明是工业文明之后人类文明发展的新阶段，它主张人类社会尊重自然、顺应自然、保护自然、合理利用自然，反对漠视自然、糟践自然、滥用自然和盲目干预自然，倡导人与自然和谐共处[1]。党的十八大首次将生态文明建设摆在突出位置，把生态文明纳入中国特色社会主义事业“五位一体”总体布局，并明确“建设生态文明是关系人民福祉、关乎民族未来的长远大计”；党的十九大进一步指出“建设生态文明是关系中华民族永续发展的千年大计”；在2018年5月召开的全国生态环境保护大会上，习近平总书记再次将生态文明建设的高度提升至“关系中华民族永续发展的根本大计”。这一系列论述凸显生态文明建设在社会主义建设事业中的地位愈发重要。

(6)

仿真结果如表1所示，与空间加权法一样，BP神经网络重构结果同样呈现出夏季重构残差最大，春秋季与冬季残差较小，磁静日重构残差低于磁扰日，磁场活动越剧烈，重构残差越大的特点.BP神经网络磁静日平均重构残差约为0.12 nT，与空间加权法重构精度(0.12 nT)及仪器分辨力(0.1 nT)相当，低于台站背景噪声水平(0.2 nT)(姚休义，2015)，重构精度较高.重构残差在静日低于扰日及其在季节变化上呈现出夏季最大、冬季略低于春秋季特点的原因均与磁场活动水平有关.地磁活动随季节变化呈现出夏季活动最强，春秋季次之，冬季活动最弱的特点(Ohtani et al.，2005；袁亚红等，2012)，因此，重构结果也随磁场活动水平差异而有所不同.此外，垂直分量Z与地下结构联系最为紧密，最能反映地下结构或观测环境的变化(丁鉴海等，1994；李琪等，2006)，对外源电流体系引起的变化弱于D和H，因此Z分量的重构残差比D和H略大.

(7)

bk=bk+ek k=1,2,…,m.

图4磁扰日数据重构结果，与磁静日有所不同，磁扰日的重构效果略低于磁静日，尤其是在磁扰剧烈时段.如图所示，在06∶11时刻，磁场活动剧烈，变化较快，DZ三要素均在该时段出现了较大的重构残差，尤其是H和Z分量，残差均高于2.0 nT.而其余时段的重构残差均在零值附近波动.计算当天Z平均重构残差Error分别为0.12 nT、0.14 nT、0.14 nT.因此，利用BP神经网络进行磁扰日数据重构时，对于磁场变化迅速的时刻会存在重构残差较大的问题，本文将对磁扰时段重构进行单独讨论.

(8)

根据误差需求，判断算法迭代是否结束，若没有结束，返回隐含层输出计算，再次迭代，经过多次迭代和修正直到输出满足需求误差的权值ωij、ωjk和阈值a、b.权值和阈值确定后，可根据预测样本(参考台数据)利用(1)和(3)式计算预测输出(目标台数据).

n=1,2,…,86400.

1.2 数据预处理

避免因输入数据之间数量级差别较大而造成网络预测误差较大，尤其是地磁观测数据由于仪器补偿不同导致的差异，在进行BP神经网络重构前要将数据进行归一化.本文采用最大最小归一化方法(式9)进行归一化处理.BP神经网络预测的数据再反归一化便可得到实际预测数据.

ASM通常被认为具有保护气管结构、调节气管管径等功能，但Mitzner等[6]认为其像“肺的阑尾”，是人类已退化的结构，并没有重要的生理作用。因此当其发生病变时，通过某种方式除去即可能治愈疾病。

(9)

2 长时间尺度数据仿真

2.1 数据仿真结果

近日，证监会核准了无锡上机数控股份有限公司（以下简称“上机数控”或“公司”）的首发申请，登陆资本市场指日可待。股市动态分析周刊记者发现，上机数控专注于精密机床领域，重视研发，连续两年营收和净利润增速超过100%，拿到IPO批文实至名归。

为了分析地磁场活动水平对重构结果的影响，对2013年每月的最静日和最扰日数据进行了仿真，利用绝对残差(Error)计算其重构精度，并与空间加权重构结果(姚休义等，2016)进行对比，Error按照公式(10)进行计算.式中，f(n)为原始地磁记录数据，为重构数据，N为缺失数据个数：

在变压器检修过程中，按照变压器检修相关标准的要求，在进行变压器排油检修，工作结束后，再次注入变压器油之前，都需要对变压器进行抽真空，在进行主变油箱抽真空期间，多次发生套管底部瓷套破裂导致中性点法兰密封面渗油的现象。

(10)

式中，η为学习速率.利用网络预测误差e更新网络节点阈值a，b.

图2 参考台选取 Fig.2 Distribution of reference stations

表1 不同季节静扰日重构残差 Table 1 Error of reconstruction in different seasons both in quiet and disturbed days

季节日期静/扰日空间加权重构残差(nT)BP重构残差(nT)DHZDHZ春秋季2013-03-08静日0.130.150.120.100.190.122013-03-17扰日0.350.340.360.360.320.362013-04-19静日0.100.130.110.090.150.122013-04-24扰日0.210.260.190.200.220.212013-09-28静日0.080.080.100.080.080.102013-09-19扰日0.170.120.180.170.120.182013-10-05静日0.100.090.110.070.080.122013-10-02扰日0.320.290.290.310.220.40夏季2013-05-29静日0.140.190.160.120.150.152013-05-01扰日0.250.190.280.250.190.282013-06-26静日0.140.140.130.140.140.132013-06-01扰日0.370.290.260.340.270.362013-07-02静日0.170.150.170.130.140.152013-07-06扰日0.180.210.220.190.170.212013-08-29静日0.110.120.120.090.100.122013-08-16扰日0.290.210.210.210.220.24冬季2013-1-22静日0.080.100.100.080.100.102013-01-26扰日0.210.210.190.150.170.172013-02-6静日0.090.130.120.080.110.122013-02-14扰日0.120.140.140.120.140.142013-11-25静日0.090.110.120.080.090.122013-11-11扰日0.180.170.160.160.130.202013-12-2静日0.080.100.100.080.100.102013-12-8扰日0.200.290.270.230.180.28静日均值静日0.110.120.120.100.120.12扰日均值扰日0.240.230.230.240.230.23总均值0.170.170.170.160.160.19

图3和图4分别以2013年12月2日和2013年2月14日为例，分析了静、扰日重构数据与原始数据的差异.A图分别为Z三分量的原始数据和重构数据曲线，其中蓝色实线为原始数据，红色实线为重构数据.整体看来，磁静日重构曲线与原始曲线几乎完全重合，各分量重构数据与原始数据之间相关系数高达1，重构结果较好.B图分别为Z三分量原始数据和重构数据的残差，各分量残差基本在0附近波动，平均重构残差Error分别为0.08 nT、0.10 nT、0.10 nT，与空间加权法结果相当，均不高于仪器分辨力0.1 nT.说明与空间加权法一样，利用BP神经网络进行磁静日地磁数据重构是可行的.

研究表明，在缺镁的土壤上施用镁肥不仅可以提高土壤中镁的含量,而且能提高作物中镁的含量，提高作物产量且改善品质。镁肥对作物的效应受多种因素制约，如土壤酸碱性、质地和交换性镁含量等，还与镁肥种类、作物种类和生育期等有关。

2.2 重构检验

虽然重构数据和原始数据曲线在形态上基本一致，且平均重构残差较小，在一定程度上保证了重构数据的精度.但是数据重构过程中是否遗漏或丢失了某些重要信息，重构数据的可用性是数据重构需要验证的一个问题.因此本文分别计算了原始数据和重构数据的Welch功率谱(Welch，1967).从频率域对原始相对记录数据和重构数据进行了分析.

当竖井掘进机偏斜角度存在时，控制系统按照纠偏算法进行合理纠偏。若偏斜角度过大，例如超过限值，则系统应发出警报，提醒操作人员注意机体偏斜。

图3 磁静日(20131202)重构结果(A)及残差(B) Fig.3 Reconstruction result of quiet day

图4 磁扰日(20130214)重构结果(A)及残差(B) Fig.4 Reconstruction result of disturbed day

图5和图6分别为静、扰日Z三分量原始数据与重构数据的功率谱密度(PSD).二者在频率范围内均表现出低频能量高于高频的特点.整体来说，重构数据及原始数据的功率谱密度曲线在形态及数值上极为一致，尤其是频率低于0.2 Hz(f<0.2 Hz)的低频部分.在个别能量变化迅速及高频部分有细微差别.计算了原始数据功率谱密度与重构数据功率谱密度之间的相关系数，磁静扰日各分量的功率谱密度相关性高达1.0，呈显著相关.高频部分多为仪器及环境噪声成分，本文在进行BP神经网络重构时以参考台数据作为输入，因此输出噪声水平应为各参考台噪声的综合水平，和目标台站噪声水平必然有所差异，因此高频部分的功率谱密度必然也有所不同.总的说来，当对数据高频部分要求不高时，利用BP神经网络方法进行数据重构是可行的，能保障重构数据的可用性.

图5 磁静日(20131202)原始数据(A)与重构数据(B)的功率谱密度对比 Fig.5 Comparison of PSD between original and reconstructed data in quiet day

图6 磁扰日(20130214)原始数据(A)与重构数据(B)的功率谱密度对比 Fig.6 Comparison of PSD between original and reconstructed data in disturbed day

3 磁扰时段数据仿真

对于全天数据缺失的情况，BP神经网络法利用相邻一天的数据作为神经网络训练输入进行缺失数据仿真，当相邻两天磁场活动水平差异较大时，对重构结果影响也较大，而空间加权法由于其利用当天各台数据进行加权重构，整体重构效果尤其是静日重构效果略好.但可以发现在磁扰日磁场活动特别剧烈的时段，BP神经网络及空间加权法重构残差均较大.因此，可考虑将磁扰剧烈时段单独截出，基于同一天的未缺失数据，利用BP神经网络方法进行短时重构.

在人才培养方面，校企合作取得了丰硕成果。2016年学校与九典制药和晶易医药联合，组建新型药物制剂研发湖南省重点实验室培育基地，共同建设一个省级科研平台，真正做到了产学研的深度融合。目前，三方共同申请并合作科研课题已有二十余项，发表高水平论文几十篇。与九典制药共同研发新药荣获“长沙市产学研与科技成果转化奖”。总之，多方面的合作成效，为校企合作展现了无限的可能。

4 实际数据重构

数据仿真结果显示，BP神经网络对于数据重构尤其是短时数据重构有精度高，残差较小的特点，能较大程度地达到完整重构，遗漏或错构数据较少.为了进一步检验该方法的适用性，利用BP神经网络方法对实际地磁观测数据中受干扰情况进行了重构.

数据仿真中假设缺失数据为红山台2012年2月14日06∶00∶00—06∶59∶59一小时的秒数据，参考台同样为泰安(TAN)、榆林(YUL)、蒙城(MCH)、乾陵(QIX)四台.为了减小日变差异及地磁活动水平差异的影响，本次BP神经网络仿真训练数据选择的是同一天中各参考台和目标台未缺失的全部数据，仿真结果如图7所示.虽然相对于其他时段，该时段的重构残差依然略大，但与全天缺失情况相比(图4)，已有较好改进.各分量尤其是偏角D和水平分量H的重构残差基本可以控制在±1 nT以内，该时段的平均重构残差由0.4 nT降低到0.2 nT.而垂直分量Z的重构残差较大主要是由于该时段泰安、榆林两台与蒙城、乾陵两台反向导致, 其反向原因可能与该地区地下结构变化导致电导率异常有关(徐文耀等，1978；姚休义等，2016).因此，BP神经网络在短时重构方面有更高的重构精度.计算该时段内原始数据与重构数据的功率谱密度(图8)，二者在变化形态上基本一致.跟全天缺失情况不同的是，二者不仅在低频部分完全吻合，即使在高频部分即噪声成分也基本一致.因此，BP神经网络更适用于短时间内的数据重构.

以红山台2013年10月4日车辆干扰为例.干扰时间较短，在连续的两小时之内，受干扰时段内磁场活动较平静.选取泰安(TAN)、榆林(YUL)、蒙城(MCH)、乾陵(QIX)四个台站观测数据作为输入进行4-5-1型BP神经构建.输入层节点传递函数选用对数S型传递函数，隐含层节点传递函数选用线性传递函数，训练函数选用动量反传的梯度下降BP算法训练，学习速率为0.9，每次训练的最大迭代次数为1000次.因仅有垂直分量Z受到干扰，因此只对Z进行重构.重构结果如图9所示，图a中灰色实线为原始受干扰数据，图b中黑色实线为重构数据，重构数据不仅避开了干扰并且完整保留了地磁场正常的日变化特征.

图7 磁扰日(20130214)短时重构结果 (a) 图黑色实线为原始数据，灰色实线为重构数据； (B) 重构数据与原始数据的残差，其中虚线部分为残差显著减小部分. Fig.7 Reconstruction result of short time Black solid line on the left indicates original data gray solid line represents reconstructed data；black solid line on the right indicates residual error between original and reconstructed data，and dashed part represent the period in that residual error decreases significantly.

5 结论

本文将BP神经网络应用到地磁观测数据的重构研究中.通过数据仿真和检验，认为BP神经网络在地磁数据重构尤其是磁扰时段的数据重构中有较高的重构精度.

图8 2013-01-24日06∶00∶00—06∶59∶59原始数据(A)与重构数据(B)功率谱密度对比 Fig.8 Comparison of PSD between original and reconstructed data from 06∶00∶00 to 06∶59∶59

图9 受干扰数据重构结果 (a) 原始数据; (b) 重构结果，其中黑色实线为重构数据. Fig.9 Construction result of disturbed data (a) Original data; (b) Construction result, and black solid line indicates reconstructed data.

(1) BP神经网络重构精度受磁场活动水平影响.不同季节及静扰日下重构残差的变化表明，磁场活动越平静，重构残差越小，重构精度越高.磁静日平均重构残差为0.11 nT，磁扰日平均重构残差为0.23 nT.

(2) D和H分量的重构残差小于Z.

(3) 重构数据与原始数据在频率域吻合程度高.原始记录数据与重构数据的功率谱密度相关性高达1，表明BP神经网络重构技术能对信号在频率域内达到较大程度的完整重构，遗漏或错构数据较少.

(4) BP神经网络短时段数据重构效果优于长时段.对长时间段(1 day)数据重构时通常选取邻近一天数据作为训练输入，而对短时间段(1 hr)数据重构时训练输入为当天的未缺失数据.通常情况下，相邻两天的磁场活动水平并不一致，因此长时间段数据重构残差必然大于短时间段重构残差.

6 讨论

磁扰日重构残差较磁静日大，可能的原因是，磁静日磁场活动弱，造成影响的电流体系单一，主要以Sq电流体系为主，因此各参考台站与目标台站之间的关系也相对简单，拟合结果较好.而磁扰日磁场活动剧烈，情况复杂时在一天之内会同时存在多种磁扰事件如磁暴、湾扰或钩扰、脉动等，造成影响的电流体系复杂，各种电流体系在空间上的分布特征不一致，如果选取的参考台区域内还存在地下结构差异导致的电导率异常，则产生的感应磁场也会有所差异，如文中提到的泰安、榆林两台与蒙城、乾陵两台垂直分量反向.这种情况下参考台站与目标台站之间的关系相对复杂，拟合结果必定不如磁静日好，所以重构残差也较大.

对于磁场活动水平较为剧烈的短时间段数据，BP神经网络重构效果明显优于空间加权法，而对于长时间(全天)缺数情况，二者的重构精度相当，均在0.11 nT左右.空间加权法中各参考台的权重只与参考台和目标台之间的距离有关，因此计算较为简单，重构速度较快(姚休义，2016)，而BP神经网络需要经过多次迭代多层训练才能达到理想输出，计算较为繁琐，重构速率较慢.因此，对于缺失数据时间较长且缺数时段内磁场活动水平较弱的情况，建议用空间加权法进行重构，而对于缺失数据时间较短且磁场活动剧烈的情况，建议采用BP神经网络进行重构.需要注意的是，由于环境和仪器噪声差异，两种方法均无法完整重构出信号高频部分，因此，对重构后的数据不建议使用其高频部分.此外，BP神经网络虽然在地磁数据重构中显示了其非线性拟合能力强、重构精度高、信号遗漏少的优点，但重构数据始终是非真实数据，可用于解决数据处理中因数据受干扰或数据缺失造成的数据处理问题，也可用于进行地磁场背景特征研究，但不能作为异常进行提取，尤其在震磁关系研究中，地震前兆异常的提取应避开数据重构时段.

References

Cai Y D, Gan J R, Yao L S. 1993. Artificial neural network in the application of the comprehensive earthquake prediction. Acta Seismologica Sinica (in Chinese), 15(2): 257-260.

Chen B, Gu Z W, Gao J T, et al. 2012. Geomagnetic secular variation in China during 2005—2010 described by IGRF-11 and its error analysis. Progress in Geophysics (in Chinese), 27(2): 512-521.

Chen F C. 1990. Back-propagation neural networks for nonlinear self-tuning adaptive control. IEEE Control Systems Magazine, 10(3): 44-48.

Ding J H, Lu Z Y, Huang X X. 1994. Earthquake Electromagnetics (in Chinese). Beijing: Seismological Press.

Dong X N, Su D L, Li X L, et al. 2012. Application of BP neural network based on rough set in the earthquake case. Journal of Seismological Research (in Chinese), 35(2): 251-259.

Feng L L. 2015. Harmonic spline modeling of regional geomagnetic field [Ph. D. thesis] (in Chinese). Beijing: Institute of Geophysics, China Earthquake Administration.

Feng Z S, Mei W P, Zhang S P, et al. 2005. Preliminary application of the daily-variation spatial correlation method of vertical component′s minutely value of FHD magnetometer. South China Journal of Seismology (in Chinese), 25(3): 1-7.

Fu C Y, Chen Y T, Qi G Z. 1985. Fundamentals of Geophysics (in Chinese). Beijing: Science Press.

Gu Z W, An Z C, Gao J T, et al. 2006. Computation and analysis of the geomagnetic field model in China and its adjacent area for 2003. Acta Seismologica Sinica, 19(2): 145-154.

Hecht-Nielsen R. 1989. Theory of the backpropagation neural network. ∥IJCNN, International Joint Conference on Neural Networks. Washington, DC: IEEE, 1: 593-605.

Hernandez J V, Tajima T, Horton W. 1993. Neural net forecasting for geomagnetic activity. Geophysical Research Letters, 20(23): 2707-2710.

Kugblenu S,Taguchi S, Okuzawa T. 1999. Prediction of the geomagnetic storm associated Dst index using an artificial neural network algorithm. Earth, Planets and Space, 51(4): 307-313.

Li Q, Lin Y F, Zeng X P. 2006. Wavelet analysis as a tool for revealing geomagnetic precursors of the Zhangbei earthquake. Chinese Journal of Geophysics (in Chinese), 49(3): 855-863.

Li Q, Gao Y F, Wang J J, et al. 2009. Local differences in great magnetic storms observed at middle and Low Latitudes. Earth, Planets and Space, 61(8): 995-1001.

Li Q, Yang X, Cai S P. 2015. Case study of applying polarization method to geomagnetic array data. Technology for Earthquake Disaster Prevention (in Chinese), 10(2): 412-417.

Li X, Zhu Q Y. 2013. Prediction of improved BP neural network by Adaboost algorithm. Computer Engineering & Science, 35(8): 96-102.

Lockwood M, Barnard L, Nevanlinna H, et al. 2013. Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr-Part1: a new geomagnetic data composite. Annales Geophysicae, 31(11): 1957-1977.

Lundstedt H, Wintoft P. 1993. Prediction of geomagnetic storms from solar wind data with the use of a neural network. Annales Geophysicae, 12(1): 19-24.

Ohtani S, Ueno G, Higuchi T, et al. 2005. Annual and semiannual variations of the location and intensity of large-scale field-aligned currents. Journal of Geophysical Research: Space Physics, 110(A1): A01216, doi:10.1029/2004JA010634.

Ren C M, Xiong J. 2012. Ionospheric electron density anomalies detected by BP artificial neural network before 2008 M8.6 Indonesia earthquake. Journal of Geodesy and Geodynamics (in Chinese), 32(S1): 28-31.

Rumelhart D E, McClelland J L. 1986. Parallel Distributed Processing: Explorations in the Microstructure of Cognition. Cambridge, MA: MIT Press.

Shen X L, Zhang J S. 2010. Research of intrusion detection based on the BP networks and the improved PSO algorithm. Computer Engineering & Science (in Chinese), 32(6): 34-36, 73.

Wang Y, Hong M H, Chen G X, et al. 2010. Nightside field-aligned current during the April 6, 2000 superstorm. Chinese Science Bulletin, 55, doi: 10.1007/s11434-010-3260-y.

Welch P. 1967. The use of fast Fourier transform for the estimation of power spectra: a method based on time averaging over short, modified periodograms. IEEE Transactions on Audio and Electroacoustics, 15(2): 70-73.

Wu W Y, Qi K, Wang S M. 1978. On the short period geomagnetic variation anomaly of the eastern Gansu province. Acta Geophysica Sinica (in Chinese), 21(3): 218-224.

Xi D Y, Zhang T. 1996. An improvement on BP neural network and the application of the theory of simulated annealing to geosciences. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese), 18(3): 219-223.

Xiong J, Wu Y, Lin J. 2013. Ionospheric electron density anomalies detected by BP artificial neural network before Wenchuan earthquake. Journal of Geodesy and Geodynamics (in Chinese), 33(1): 13-16.

Xu W Y. 2003. Geomagnetism (in Chinese). Beijing:Seismological Press.

Xu W Y. 2009. Yesterday, today and tomorrow of geomagnetic indices. Progress in Geophysics (in Chinese), 24(3): 830-841.

Yao F Z. 1988. The application of correction method in geomagnetic absolute observation. Seismological and Geomagnetic Observation and Research (in Chinese), (6): 35-43.

Yao X Y. 2015. Method of artificial electromagnetic disturbances and data reconstruction [Ph. D. thesis] (in Chinese). Beijing: Institute of Geophysics China Earthquake Administration, 47-54.

Yao X Y, Teng Y T, Yang D M, et al. 2016. On reconstruction technique of geomagnetic observation data. Acta Seismologica Sinica (in Chinese), 38(6): 878-888.

Yuan Y H, Yang D M, Chen H R, et al. 2012. Analysis of the temporal-spatial distribution characteristics of geomagnetic activity index Vr. Chinese Journal of Geophysics (in Chinese), 55(3): 960-969, doi: 10.6038/j.issn.0001-5733.2012.03.026.

Zhang S Q, Fu C H, Zhao X D. 2016. Study of regional geomagnetic model of Fujian and adjacent areas based on 3D Taylor Polynomial model. Chinese Journal of Geophysics (in Chinese), 59(6): 1948-1956, doi: 10.6038/cjg20160602.

Zhao X D, Du A M, Chen H R, et al. 2010. Inversion of the Sq current system and the geomagnetic diurnal variation model. Progress in Geophysics (in Chinese), 25(6): 1959-1967, doi:10.3969/j.issn.1004-2903.2010.06.010.

Zhu Z C. 1989. Spatial correlation analysis in the application of the geomagnetic observational studies. Seismological and Geomagnetic Observation and Research (in Chinese), (5): 44-49.

附中文参考文献

蔡煜东, 甘俊人, 姚林声. 1993. 人工神经网络在地震综合预报中的应用. 地震学报, 15(2): 257-260.

陈斌, 顾左文, 高金田等. 2012. IGRF-11描述的2005—2010年中国地区地磁长期变化及其误差分析. 地球物理学进展, 27(2): 512-521.

丁鉴海, 卢振业, 黄雪香. 1994. 地震电磁学. 北京: 地震出版社.

董晓娜, 苏道磊, 李希亮等. 2012. 基于粗糙集的BP神经网络在震例中的应用研究. 地震研究, 35(2): 251-259.

冯丽丽. 2015. 区域地磁场谐和样条模型. 北京: 中国地震局地球物理研究所.

冯志生, 梅卫萍, 张苏平等. 2005. FHD磁力仪Z分量分钟值日变化空间相关性的初步应用. 华南地震, 25(3): 1-7.

傅承义, 陈运泰, 祁贵仲. 1985. 地球物理学基础. 北京: 科学出版社.

李琪, 林云芳, 曾小苹. 2006. 应用小波变换提取张北地震的震磁效应. 地球物理学报, 49(3): 855-863.

李琪, 杨星, 蔡绍平. 2015. 极化方法应用于地磁台阵的震例分析. 震灾防御技术, 10(2): 412-417.

李翔, 朱全银. 2013. Adaboost算法改进BP神经网络预测研究. 计算机工程与科学, 35(8): 96-102.

任成冕, 熊晶. 2012. 利用BP神经网络探测2012年印尼8.6级地震前的电离层异常扰动. 大地测量与地球动力学, 32(S1): 28-31.

沈学利, 张纪锁. 2010. 基于BP网络与改进的PSO算法的入侵检测研究. 计算机工程与科学, 32(6): 34-36, 73.

王源, 洪明华, 陈耿雄等. 2010. 2000年4月6日磁暴夜侧场向电流变化特征. 科学通报, 55(14): 1409-1415.

徐文耀, 祁暌, 王仕明. 1978. 甘肃省东部地区短周期地磁变化异常及其与地震的关系. 地球物理学报, 21(3): 218-224.

席道瑛, 张涛. 1996. BP神经网络的改进及模拟退火神经网络在地学中应用. 物探化探计算技术, 18(3): 219-223.

熊晶, 吴云, 林剑. 2013. 用BP神经网络技术探测汶川地震前电离层NmF2异常扰动. 大地测量与地球动力学, 33(1): 13-16.

徐文耀. 2003. 地磁学. 北京: 地震出版社.

徐文耀. 2009. 地磁活动指数的过去、现在和未来. 地球物理学进展, 24(3): 830-841.

姚法章. 1988. 订正法在地磁绝对观测中的应用. 地震地磁观测与研究. 6:35-43.

姚休义. 2015. 地磁台站观测异常识别与数据重构技术研究[博士

论文]. 北京: 中国地震局地球物理研究所, 47-54.

姚休义, 滕云田, 杨冬梅等. 2016. 地磁观测数据重构技术研究. 地震学报, 38(6): 878-888.

袁亚红, 杨冬梅, 陈化然等. 2012. 地磁活动指数Vr的时空分布特征分析. 地球物理学报, 55(3): 960-969,doi: 10.6038/j.issn.0001-5733.2012.03.026.

张素琴, 付长华, 赵旭东. 2016. 基于三维Taylor模型的福建及邻近地区地磁模型研究. 地球物理学报, 59(6): 1948-1956, doi: 10.6038/cjg20160602.

赵旭东, 杜爱民, 陈化然等. 2010. Sq电流体系的反演与地磁日变模型的建立. 地球物理学进展, 25(6): 1959-1967, doi:10.3969/j.issn.1004-2903.2010.06.010.

朱兆才. 1989. 空间相关性分析在地磁观测研究中的应用. 地震地磁观测与研究, (5): 44-49.