水下辅助导航的重力异常随机场分形插值研究

高伟，鲁国瑞

( 天津城建大学地质与测绘学院, 天津 300384 )

0 引言

地球物理场和惯性导航系统(inertial navigation system，INS)联合的辅助导航技术可有效遏制INS误差累积问题，始终是国内外水下导航领域的研究难点和热点. 水下辅助导航技术主要包括地磁匹配导航、地形匹配导航和重力匹配导航等. 而水下重力匹配导航通过测量重力信息与地球重力场进行匹配实现导航定位，在重力辅助导航系统运行时，采集重力特征显著区域的重力信息与重力基准图对比得到位置信息，实现天空海一体化水下潜器惯性/重力组合导航系统重调[1-4].

地球形状的不规则性和密度的不均匀性导致地球各点的重力场并不相同，表现为空间位置(经度、纬度、高度)的函数[5]. 因此潜航器可在航行过程中，通过重力测量仪器采集航线途经的重力数据，并与预先存储的重力数据相匹配，获得潜航器当前定位信息，进而对INS 累积的位置误差进行校正[6-7]. 重力辅助导航系统在测量重力场数据过程中，潜航器无需露出或接近水面，测量仪器也无需向外部发射信号或接收外部信号，系统可进行无源、隐蔽的导航定位，潜航器在卫星、无线电定位系统失效或遭受破坏的特殊情况下仍可达到自主隐蔽导航的目的[8-10].

随着重力测量仪器及空间测量技术的发展[11]，在全球范围内快速、准确地获取重力数据成为现实，使得重力辅助导航系统具备校正INS 积累位置误差的能力[12-13]. 重力辅助导航系统主要有两种代表性系统，均由Bell Aerospace 研发，分别为于1990 年研制出的重力梯度导航系统[14]，以及于1991 年研制出的重力辅助INS[15]. 重力辅助导航系统由INS、重力测量仪器、数字重力基准图[16]和匹配定位算法等四个主要部分组成，各组成部分的性能差异对重力辅助定位系统的性能有重要影响. 数字重力基准图(重力模型)是重力辅助定位的基础[17-18]，重力基准图描述是否准确，包含的重力特征是否丰富，分辨率是否满足需求，都会影响重力辅助导航系统的性能[19].

为减小线性化误差，对重力异常参考地图进行预处理，选择合适的计算模型对地图进行插值逼近，得到高分辨率、小格网地图，以减小线性化误差对位置观测精度的影响. 在重力异常数据匮乏的地区，随着相邻点间的插值数的增加，拟合函数和Kriging 方法无法满足逼近精度要求. 其中，拟合函数法跟插值的规模、控制点数目、控制点统计特征和地图统计特征差异、地图统计特征相关，不适宜重力异常数据匮乏的区域，而Kriging 方法不适合大尺度条件下的重力插值，因此采用分形重构模拟真实的重力环境，为系统提供良好的仿真环境. 分形通常被用在模拟真实地表和其他自然形态的描述中，它用分维数表达数据场的统计特征. 与Kriging 方法的空间变异函数相似，低分辨率地图不利于长距离的统计特征描述[20-22]. 本文考虑在这种条件下，提出利用分形插值的方法对重力异常随机数据场进行重构，重构模拟将为水下重力辅助INS 提供更接近真实的重力环境[23].

分形就是对不规则和支离破碎的自然形状的描述，分形理论的发展促进了利用它对各种自然形态描述的研究. Goodchild 在1980 年考虑用分形理论分析地文表面，后来Clark 和Schweitzer 定义了具有更好鲁棒性的对地形表面分形估计. 从1991 年开始，Kenichi Arakawa 和Eric Krotkov 对“利用分形重构自然地形”进行了研究，利用分形布朗运动(fractal Brownian motion, FBM)在不规则采样点集上估计自然地形的分形结构，并以此为基础构造出自然地形，Naokazu 等给出基于FBM 模型的插值方法并不断进行完善[24].

本文提出采用FBM 分形插值方法对重力异常随机场进行重构，构造可用于水下重力辅助导航系统的重力异常基准图，通过实验验证分形插值方法的可行性，将为水下重力辅助导航系统的实现提供理论参考.

1 FBM 函数及分维数估计

布朗运动(Brownian Motion，BM)是1982 年英国植物学家R. Brown 发现的，1923 年德国数学家N.Wiener 建立了BM 的数学模型. BM 是一种随机过程B(t)，满足以下两个条件：

式中：C为常数；随机过程的任意两点之间(B(b)-B(a)和B(c)-B(b),a≤b≤c) 的增量是相互独立的.

B(t)被写成

显然，它具有自紧密性.

Mandelbrot 提出FBM 函数，由它产生BM 过程.

Penland 定义了FBM 函数统计自相似性特征

式中：f(t) 为FBM 函数；g(x) 为累计分布函数(g(x)～N(0,σ2))；H为自仿射参数，它与f(t) 的分维数D之间的关系为

n是几何拓朴维数，在重力异常地图中，n=2 .

由式(6)可得

又

可以得到

两边取对数

在对二维随机数据场进行分维数估计时， Δt用地图上两点距离代替，则

2 重力异常地图分维数估计

重力异常地图的分维数的估计也是以式(11)为基础进行的.

在计算时需要得到在不同 ΔX下的E[|f(X+ΔX)-f(X)|]，显然在重力异常地图阵列的基础上，可以得到的 ΔX是离散的，它的取值与地图网格单位距离和地图阵列的规模有关，如图1 所示.

图1 地图网格示意图

能取的距离为ΔX

式中，i=1,2,···,mx；j=1,2,···,my；这里仅讨论ux=uy=u的情形.

显然，ΔX增加，在计算时，可参与统计点数目会急剧减少. 基于该思路，设定距离容差，可以更好地反映地图的统计特征.

即，设 εx为距离容差，当地图上两点距离满足

那么，两点重力异常增量可以计入与ΔX对应的计算中.其中，ΔXP是两点之间的实际网格距离.

按此方法分别对2 个地块区域的重力异常地图进行了分维数计算，结果如图2 所示.

图2 重力异常地图分形特征曲线

用最小二乘法拟合直线，得到的结果如表1 所示.

表1 分形特征计算结果

3 FBM 分形插值方法

用FBM 分形插值，生成新的小尺度重力异常地图的基本步骤如图3 所示.

图3 分形插值方法基本流程

分形特征一个有限的尺度范围内才能充分体现，这个范围称为无标度区间 [ΔXmin,ΔXmax] .

由地图得到的分形特征无标度区间的下限由地图间隔决定

其上限为分形特征曲线线性拟合区域的上限.

分形插值后得到地图的最小尺度uc小于u，因此分形插值所能得到的最小尺度不能由原始地图确定.本文所讨论的重力异常地图分形插值只能用于重力异常地图的粗插值.

分形插值采用随机中点位移法，即插入点附近区域的点理解为一个符合 (0,σ2) 的随机数据与区域的平均值的和，它的实现如图4 所示.

图4 中点位移插值法

在图4 中，填充点为原始地图的控制点，虚线是插值点；插值点分为两类，标号为1 的插值点 (x,y) ，由其周围的4 个原地图控制点为参照进行插值

标号为2 的插值点，所用到的参照点由两个原地图控制点和标号为1 的两个插值点构成

其中， Gauss 是符合均值为0，标准差为1 的随机信号.

4 实验结果与分析

4.1 与原始地图比较

针对某块 60×60 的重力异常地图进行插值，在原始地图上进行等间隔抽样得到一块 30×30 的地图，然后用上述分形插值方法进行插值，并对插值地图和原始地图进行比较. 抽样后的地图分形特征如表2所示.

表2 分形特征参数表

插值后的重力异常地图与原始地图的比较如图5所示.

图5 在原始数据上抽样后分形插值结果

插值后的精度分析如表3 所示，采用误差平均值(E(ge))、绝对误差平均值、绝对误差最大值(max(|ge|))、误差均方值和绝对误差标准差(D(ge) )等参数来评定.

表3 分形插值精度分析

4.2 与原重力异常地图比较

在原重力异常地图上进行插值，分析插值后地图的分形特征的变化. 图6 是插值后的等值线图比较.图7 是插值后地图分形特征曲线. 其特征参数见表2第2 行数据.

图6 原始数据直接分形插值结果

图7 插值后地图分形特征曲线

显然，插值后的地图与插值前地图的分形特征变化很小. 分形插值的精度与标准差相关，标准差越大，精度越低；分形插值引入随机信号，因此插值点的重力异常值具有不确定性；采用中点位移法的插值，插值点位置固定，即uc=u/2 .

分形插值作为地形自然表面的重构方法，具有很强的真实感. 重力异常与地形之间的关联性，产生了利用分形对重力异常进行插值的联想，为得到可靠的重力异常地图的分维数，在分维数定义中加入距离容差，改善了分维数在长空间距离下的统计特征描述.分形插值的精度虽然低，但插值后的地图仍保持原重力异常参考地图的统计特征，该方法为仿真提供了很好的环境模拟. 由于条件的限制，虽然没有对插值结果进行实际测量的验证，但从计算机仿真实验结果看，分形插值作为一种粗插值是可行的.

本文是针对水下重力数据匮乏区域的重力场仿真研究，由于数据异常和条件差，分形插值的实验结果仅与原始数据的计算结果进行了对比研究. 由于重力数据的匮乏，采用其他插值方法的效果不太理想，限于篇幅和数据条件，未进行详细分析，但后续仍需要通过获取大量实测数据来进行反复验证与探讨.

5 结束语

本文针对重力异常数据匮乏的地区，提出了采用FBM 分形插值方法对重力异常随机场进行重构，实现可用于水下重力辅助导航的重力异常基准图生成. 实验结果表明，分形插值后的地图与插值前地图的分形特征变化很小. 分形插值的精度与标准差相关，标准差越大，精度越低；分形插值引入随机信号，因此插值点的重力异常值具有不确定性；采用中点位移法的插值，插值点位置可固定. 分形插值作为地形自然表面的重构方法，具有很强的真实感. 从数值模拟结果来看，分形插值方法作为一种粗插值是可行的和有效的，可为水下重力辅助INS 提供更接近真实的重力环境，为水下重力辅助导航系统的实现提供理参考.