1、第 卷 第期 年月大 地 测 量 与 地 球 动 力 学 ，收稿日期：项目来源：国家自然科学基金（，）。第一作者简介：于文强，硕士生，主要从事地磁场模型研究，：。通讯作者：李厚朴，博士，教授，博士生导师，主要从事海洋大地测量研究，：。：文章编号：（）混沌 神经网络的地磁变化场预测模型于文强李厚朴秦清亮宋立忠王志远海军工程大学电气工程学院，武汉市解放大道 号，海洋环境专项办公室，北京市，部队，青岛市，摘要：提出一种融合混沌理论和 （）神经网络的地磁变化单站预测模型。分析磁场数据的混沌特性，求取关键参数嵌入维数和时间延迟，据此对初始数据进行相空间重构，并将经混沌理论优化的样本集作为神经网络的训练集
2、和测试集进行仿真实验。结果表明，经混沌理论改进后的 神经网络模型可以较为准确地预测地球磁场的变化趋势，对我国地磁场的适用性较好，具有一定的泛化能力。关键词：混沌理论；相空间重构；神经网络；地磁变化中图分类号：文献标识码：地磁变化场是地磁学和空间物理学的研究重点。叠加在地球稳定磁场上的变化磁场可以分为类：）随时间缓慢变化的长期变化地磁场，其变化周期一般为几十年甚至更长；）受到固体地球外部各种空间电流体系影响产生的短期变化地磁场。短期变化地磁场又可分为平静变化和扰动变化类：地磁平静变化比较平缓，有一定的周期性，分为太阳静日变化和太阴日变化；扰动变化的发生较为偶然，持续一段时间后消失，变化短暂且不规
3、律。变化磁场在地球总磁场中占比仅左右，但却对地磁导航影响巨大。通常情况下，地磁平均日变幅可达几十 ，因此对于地磁导航而言，地磁图除了需要高精度的稳定磁场外，还需要高时空分辨率的变化磁场，而现有的地磁变化场模型的时间分辨率不足以满足地磁导航的需求。针对变化磁场，传统的建模方法是根据其物理机制、观测数据、形成猜想以及物理方程建立空间电流体系模型，但其物理过程极为复杂、理论化难度较大，无法达到相应的导航精度。随着计算机技术和人工智能的发展，在变化磁场建模的过程中可以依靠神经网络等方法避开复杂的物理机制，仅在数据驱动下开展研究。首次提出将人工智能应用于地磁日变化模型的开发；等建立地磁变化场多时间尺度组
4、合模型，该方法结合可变地磁场的物理特征，融合多个模型，预测 时 间 跨 度 较 长；等提 出 一 种 神经网络算法，并对和 等地磁活动指数进行预测；程文凯等提出基于 的机器学习方法重构地磁日变数据，相较于反向传播神经网络，其精度更高。本文提出融合混沌理论和 神经网络的预测模 型，在 变 化 磁 场 的 重 构 相 空 间 中，利 用 神经网络逼近相点的演化规律，对变化磁场进行预测。理论与方法 混沌理论混沌是非线性动力系统的固有特性，广泛存在于自然系统和社会系统中。混沌系统是指在一个确定的系统里，存在看似随机发展的不规则运动，这种运动无周期性、不收敛，具有不确定性和不可重复性，并且对于初始值极为
5、敏感。混沌理论是非线性系统在一定的参数条件下，发生分叉、周期与非周期运动相互纠缠而产生的某种非周期有序运动的理论。相空间重构是混沌时间序列预测的重要基础理论，可以从单一时间序列中提取和恢复出原系统的规律。假设变化磁场时间序列为，用时间延迟的方法构造第 卷第期于文强等：混沌 神经网络的地磁变化场预测模型（）个维相空间矢量：（）（，（）（）重构轨道为：，（）式中，和分别为嵌入维数和延迟时间，为维矩阵，由个分量组成，个相点在维相空间中构成一个相型，各相点连线即为描述动力系统在维相空间里的演化轨迹。在相空间重构过程中，合适的和对于地磁变化时间序列中的信息挖掘具有重要意义。地磁变化场起源于磁层和电流层体
6、系，是一个极为复杂的系统，难以用理论模型进行精准建模。但变化磁场数据是一个非线性、非平稳的时间序列，具有混沌特性，因此本文引入混沌理论来挖掘数据深层特征，利用数据驱动寻找隐藏规律，从而完成预测。神经网络 神经网络是一种使用径向基函数作为激活函数的人工神经网络，其输出是输入的径向基函数和神经元参数的线性组合。神经网络由层结构组成，网络拓扑结构见图。第层是输入层，其节点个数等于输入向量维数；第层是隐含层，该层是网络的核心部分，决定网络性能，其变换函数是局部分布的非负非线性函数，对中心点径向对称且逐渐减弱；第层是输出层，该层是隐含层输出的线性加权求和。神经网络有较强的线性拟合能力，收敛速度较快，有较
7、好的泛化能力，可用于构建变化磁场预测模型。图 神经网络拓扑结构 地磁变化预测模型构建 地磁变化场混沌特性检验首先使用李雅普诺夫特性指数（，）判断地磁变化时间序列是否具有混沌特性。假设某一段时间内采集的地磁变化场时间序列为 ，具体判断步骤如下：）通过 变换计算地磁变化时间序列的延迟时间和平均周期。）根据的原则确定嵌入维数，根据和重构相空间：（）（，（）（）式中，（）。将新产生的序列定义为 ，。）根据最小值原理找出离中心点最近的邻点（?），计算最短距离：（）（?）（）（?）（）式中，?，且?。）计算中心点（）的（）：（）（）（?），（，?）（）计算平均（）：（）（）（）式中，为（）的数量。由最小二
8、乘法计算式（）中（）的偏微分，从而求得上式斜率，即最大李雅普诺夫指数。当时，系统有稳定的不动点；当时，对应分叉点或系统的周期解，即系统出现周期现象；当时，系统具有混沌特性。伪最近邻域法确定嵌入维数混沌时间序列是高维相空间混沌运动的轨迹在低维空间上的投影，在投影过程中，混沌轨迹会发生一定程度的扭曲，可能导致高维空间中不相邻的个相点投影到低维空间后成为相邻相点，即伪最近邻点，导致混沌时间序列无规律。伪最近邻域法确定嵌入维数的方法如下：在地磁变化时间序列重构后的维相空间中，每个相点（）（），（），（）都存在某个距离内的最近邻点，采用欧氏 度 量 计 算 点 间 距（）（）（）。将维相空间扩展到（）维
9、上，相点间距会随之发生变化，变为（）（）（）（）。若（）与（）差别较大，则可以认为是高维混沌吸引子个点在投影到低维空间上时变为伪最近邻点。令（）（）（）（）若，则判定（）是（）的虚假最近邻点。为阈值，可在，之间选择。从最小嵌入维数开始计算伪最近邻点比值，逐渐增大嵌入维数，当增大到伪最近邻大 地 测 量 与 地 球 动 力 学 年月点比值小于或伪最近邻点不随嵌入维数的增加而减少时，可以认为奇异吸引子已经完全展开，此时的即为最佳嵌入维数。自相关函数法确定延迟时间对于地磁变化序列 ，其时间跨度为 的自相关函数定义为：（）（）（）（）固定，绘制自相关函数关于时间的函数图像，当自相关函数下降到初始值的倍
10、时，所得时间即为延迟时间。混沌 神经网络组合预测在 神经网络的基础上，融合混沌理论对样本集进行优化，得到更适用于预测地磁变化的混沌 神经网络模型，具体步骤如下。）将已知的地磁变化垂向分量数据作为预测的初始时间序列，根据混沌理论中的伪最近邻域法和自相关函数法计算最优嵌入维数和延迟时间，重构相空间，得到经混沌参数优化的数据集：，（）（）（，（）（）将经相空间重构后的地磁变化垂向分量数据作为 神经网络输入数据，一部分数据作为训练样本集，另一部分数据作为测试样本集。）利用优化后的训练样本对神经网络进行训练，并通过调试参数构建合理网络，最终通过测试样本得到地磁变化垂向分量的预测数据。实验验证使用个台站的
11、数据进行实验验证，分别是北京十三陵（）站、四川郫县（）站、黑龙江漠河（）站、西藏拉萨（）站、广州肇庆（）站，所选台站分布较分散，可以代表我国大部分区域的地磁变化特征（图）。使用磁通门磁力仪测得秒采样数据，时间跨度为 。在实际应用中，匹配分量的时间分辨率越高越好，但是过高的分辨率会导致计算量剧增、预测时间缩短。目前常用的地磁匹配分量为地磁垂向分量和磁场总强度。综合考虑后，采用的 均值作为实验数据，并以样本集前 的数据训练预测模型，选取周的测试数据作为预测检验样本集。站建模预测由于混沌动力系统中噪声信号的不可预测性和高破坏性会掩盖系统的内在动力学特性，影响图台站分布 后续的预测精度以及对系统的准确
12、描述，因此首先要对地磁数据进行去噪和平滑处理。采用较为简单的拉依达准则剔除粗大误差，以 站为例进行建模及预测。首先使用 计算得到地磁变化分量数据的最大李雅普诺夫指数 ，可知地磁变化分量具有混沌特性。然后采用自相关函数法和伪最近邻域法确定延迟时间和嵌入维数，结果见图、。图自相关函数法计算时间延迟 图伪最近邻域法计算嵌入维数 图中红色直线为，由图可知，当自相关函数 下降到初始值的倍时，分量的延迟时间。由图可见，伪邻近第 卷第期于文强等：混沌 神经网络的地磁变化场预测模型率随的增大逐渐减小，并在 时成为常数值，不再随发生变化，说明此时已经能够还原混沌吸引子的结构，即为模型的最佳嵌入维数。当 神经网络
13、每层神经元数目与混沌时间序列重构相空间的嵌入维数相等时，预测效果最佳。因此，本文选取的神经元数目为。将经过平滑和归一化处理后的研究时段前 共 个分量数据的前 个作为训练样本集，经相空间重构后分为训练输入样本集和训练期望输出样本集；将后 个样本作为需要预测的未知样本，经相空间重构后分为测试输入样本集和测试期望输出样本集。选取训练输入样本集中最后（）个数据作为预测输入样本集进行递推预测，将预测结果与传统 神经网络算法的预测结果进行比较（图）。其中，测试输入样本集不参与预测，期望输出样本集作为模型预测结果的参照值。图算法改进前后预测结果及误差绝对值 由图可见，种模型均可以跟踪预测分量的趋势变化，但与
14、传统 神经网络算法相比，混沌 神经网络模型的预测结果明显更接近于实际观测值，这种优势在地磁平缓变化时更加显著。在地磁场发生较大变化时，预测拟合程度是评价模型性能的重要指标。可以发现，相比于 神经网络算法，混沌 神经网络模型对曲线突变的反映效果更好，但是误差仍偏大，这可能是因为受到太阳活动等因素的影响发生地磁扰动时，地磁变化场复杂性增大、可预测性降低。从图还可以看出，神经网络模型误差的绝对值大约在 以内，且离散程度较大；本文改进模型误差的绝对值大约在 以内，且分布较为集中，精度提高约。模型泛用性评估为了评估改进模型在我国其他区域的适用性，对个台站进行建模预测，并计算模型结果的最大误差、平 均 绝
15、对误 差 和均方根误差，结果见表。表各台站误差结果 台站最大误差 由表可见，预测结果的 均小于，均小于，误差较小，说明基于混沌 神经网络的地磁变化预测模型能较好地描述变化磁场的特性。模型对 和 站的预测效果较好，小于，小于，但对 和 站的预测结果误差偏大。为了更准确地研究模型的预测精度，计算各台站误差范围所占比例，结果见表（单位）。表绝对误差统计 由表可见，所有台站的预测结果绝对误差小于 的占比均超过，小于 的占比均超过，除 站以外，其他台站小于 的占比均超过；站的预测效果最好，绝对误差小于 的占比超过。综合分析可知，本文模型对中高纬度地区的预测精度略高于低纬度地区。结语）混沌 神经网络地磁变
16、化单站预测模型能够有效预测地磁变化垂向分量的变化趋势，相较于传统的预测模型，该模型的预测效果更好，精度提高约。）模型预测精度对纬度变化较为敏感，中高纬度地区精度较高，不超过；低纬度地区精度相对较低，约为。参考文献牛超，李夕海，魏一苇区域地磁变化场分析与建模关键技术研究西安：西安电子科技大学出版社，（，大 地 测 量 与 地 球 动 力 学 年月 ：，）管志宁地磁场与磁力勘探北京：地质出版社，（：，）王仕成地磁变化场对地磁匹配制导影响分析西安：西安地图出版社，（：，）牛超，齐玮，李夕海，等地磁场静日变化 研究概述第届国 家 安 全 地 球 物 理 专 题 研 讨 会，北 京，（，），（）：，（），程文凯，杜劲松，陈超，等基于 机器学习的地磁日变重构方法研究地震学报，（）：（，（）：）韩敏，任伟杰，李柏松，等混沌时间序列分析与预测研究综述 信 息 与 控 制，（）：（，（）：）李松，刘力军混沌时间序列智能预测方法及其应用北京：科 学 出 版 社，（，：，）王赤，陈金波，王水地球变化磁场的分形和混沌特征地球物理学报，（）：（，（）：），：，（）韩敏混沌时间序列预测理论与方法北京：中国水利水电
