作者简介:赵塔(1987-),男,博士研究生.研究方向:弱磁信号检测,地磁导航.E-mail:zhao_ta_ok@126.com
为了实现水下载体对地磁场的精确测量,提高地磁辅助导航的精度,提出了一种基于多传感器的空间差分地磁场三分量测量方法。首先,分析了水下载体周围磁场的种类及空间分布特性,并建立了磁矢量传感器新的测量模型。然后,利用干扰磁场自身特性以及和地磁场之间的耦合关系构建方程组,推导出空间差分模型表达式。最后,通过仿真分析和地面实验分别将本文方法与传统测量方法进行了对比。结果表明:当外界干扰磁场为1000 nT时,差分测量结果能使测量误差控制在150 nT以内,很好地抑制了干扰源对地磁场测量造成的影响,提高了水下载体地磁场三分量测量的抗干扰能力。
In order to accurately measure the geomagnetic field around an underwater vehicle and improve the precision of geomagnetic aided navigation, a spatial difference measuring method is proposed based on multi-sensor. First, a new measurement model of magnetic vector sensors is established by analyzing the magnetic field types and the spatial distribution characteristics around the vehicle. Then, the spatial difference model expression is deduced according to the characteristics of interference magnetic field and its coupling relationship with the magnetic field. The proposed method was compared with traditional method by simulation and ground experiment. Results show that the proposed method can obviously restrain the impact of the interference magnetic field. When the intensity of the outside interference magnetic field is 1000 nT, the measurement error can be controlled under 150 nT by the spatial difference measuring method, which is able to improve the ability of anti-jamming in geomagnetic field measurement.
进入21世纪以来, 世界上已有10多个国家的1000余艘各种用途的水下载体投入到军用或民用的领域, 其导航方式也根据应用场合的不同逐渐从单一向组合导航系统发展[1]。地磁/惯性组合导航系统能够有效解决惯性导航系统误差随时间累积问题, 既不要求水下载体上浮, 又不向外辐射能量, 具有极好的隐蔽性, 大大提高了载体在水下的生存和应用能力[2, 3]。然而构建足够精确的水下地磁数据库或地磁图是实现地磁/惯性组合导航技术的关键问题。因此, 为了详细了解水下环境磁场特性, 获得更加准确的地磁数据, 有必要进行水下地磁场的实地测量。
由于单一总场存在地磁等值线的问题, 而地磁场的三分量包含更多的地球物理信息。因此, 地磁场的矢量测量必然成为导航技术发展的趋势。日本最早提出一种地磁场三分量的船载式测量方法(STCM), 主要通过传统测量模型对载体自身干扰磁场进行补偿, 地磁场三分量测量精度为(50± 25) nT[4]。国内的闫辉、肖昌汉等[5, 6]也针对海洋环境地磁场三分量的测量进行了研究, 并采用“ 两步法” 验证了测量方法的正确性和可行性。但是上述方法主要采取单传感器测量, 抗干扰能力差, 而且没有考虑外界磁源产生的干扰磁场对测量系统造成的影响, 这也是地磁场实际测量造成误差的原因之一。
本文在STCM的基础之上, 充分考虑了干扰磁场的空间分布特性和来源情况, 通过多传感器测量组建的方程组推导出空间差分模型, 最后通过仿真和地面实验验证了本文方法的正确性和可行性。
目前, 水下大部分载体由铁磁性材料建造而成, 可以将其分为硬磁材料和软磁材料, 它们在地磁场的作用下会被磁化而向外辐射磁场。
同时, 水下环境复杂多变, 载体周围会随时出现一些其他磁源, 这些未知的磁源将产生一些杂散磁场, 从而对测量造成影响。因此, 如图1所示, 上述多种混合磁场的矢量叠加构成了载体周围干扰磁场的主要成分。
硬磁材料具有较高的矫顽力和剩余磁场值, 由于载体在建造过程中或长期停靠时期受到相同方向地磁场的磁化, 产生的这部分剩磁可以保留较长时间且不随载体的航向和位置发生变化[7, 8], 因此也称之为固定磁场, 其矢量在载体坐标系中表示为
软磁材料的矫顽力较低, 且具有较窄的磁滞回线, 能够在较弱的磁场中被磁化而产生感应磁场, 它会随着外磁场的改变而改变, 也会随着载体形状、姿态、位置的不同而发生变化。根据泊松方程, 感应磁场矢量和地磁场矢量之间满足如下关系:
式中:
式中:由
式中:HM 为载体坐标系下传感器的矢量测量值; E为单位矩阵。
在水下环境中, 载体周围还可能存在着很多未知的其他磁源, 如水面上的船只, 海底雷区或沉船、水下的其他运行载体等, 这些铁磁性的物体也会被磁化而向外辐射磁场, 产生的杂散磁场用∑ HN表示。∑ HN是多种磁场的矢量叠加, 具有很强的突发性和不确定性, 既无法测量其大小也无法判断其方向, 但它确实是测量中不可避免而又容易忽视的一种干扰磁场成分。由于在传统模型方法中采用单传感器测量, 模型中无法去除∑ HN而忽略了这部分磁场, 但这部分磁场一旦出现会直接叠加到传感器的测量值中, 对地磁场矢量的最终计算造成较大误差。如何合理地去除∑ HN, 提高地磁场矢量测量精度成为了研究的重点问题。因此, 载体坐标系下传感器更加完整的新测量模型为:
根据磁场强度在水中衰减的经验公式[9]:
式中:
磁场强度随距离衰减曲线如图2所示, 在前20 m的距离磁场梯度较大, 磁场强度衰减了约104倍, 称该区域为快速变化区。而随着距离的增加磁场曲线的斜率也趋于平缓, 即磁场梯度越来越小, 称该区域为缓慢变化区。
通过上面的分析, 将干扰磁场分为自身干扰磁场和外界干扰磁场。自身干扰磁场包括
图3为水下载体坐标系及传感器安装图, x轴在载体左右对称的纵剖面上, 指向载体前进方向为正; y轴在载体中心平面上, 指向右侧为正; z轴垂直xy平面, 向下为正。在载体上中轴线两侧非对称位置捷联安装两个三分量传感器, 根据载体尺寸和传感器测量精度合理选择安装距离, 一般控制在2~5 m以内。由于载体自身产生的干扰磁场在两传感器处产生的磁场梯度较大, 属于图2中快速变化区磁场, 在差分计算时具有明显的差异性, 而∑ HN是由较远距离的单个或多个磁化体产生, 一般距离大于20 m, 属于图2中缓慢变化区磁场, 在载体处产生的磁场梯度很小, 在两个传感器处的响应值大致相同, 因此在相同坐标系下可以通过差分处理将其去掉。
地磁场是一种天然的物理场, 当不存在磁异常点时它的各个分量自身梯度很小, 中国从南到北地磁场水平分量强度主要随纬度变化[10], 水平分量强度每纬度最大变化了10 nT, 垂直分量强度每纬度最大变化了7 nT。而根据两个传感器的实际安装距离, 在载体坐标系下传感器测量的地磁场三分量值与载体位置的地磁场三分量值基本相同, 满足下面关系:
根据式(4)对两个传感器新的测量模型构建如下方程组:
通过上面的分析可知, 外界杂散干扰磁场在两个传感器处满足:
对式(7)中的两个方程做差, 得到:
式(9)为地磁场三分量的差分测量模型, 感应系数矩阵
式中:
为了验证本文提出的空间差分测量方法的可行性, 利用国际标准地磁场IGRF-11/2010模型中提供的地磁场数据, 按照所述的原理应用大型有限元仿真软件Comsol Multiphysics针对本文方法进行仿真验证。
载体模型建立条件如下:采用单椭球模型来构建载体, 模型为长度50 m, 宽度和高度均为10 m的椭球形空腔, 腔体厚度为0.05 m, 后端带有尾翼, 材料的相对磁导率为100。在上表面中轴线两侧捷联式安装两个传感器, 左侧传感器离载体表面高度为0.2 m, 右侧传感器离载体表面高度为1.2 m, 两个传感器的直线距离为2.0 m。在长、宽、高均为100 m的正方体求解域中加入
在理想环境(∑ HN=0, HG变化平缓且已知)中进行感应系数矩阵和固定磁场的仿真求解。首先设置两个传感器处的固定磁场三分量分别为675.0、-575.0、342.0 nT和485.0、-185.0、852.0 nT。然后根据IGRF-11/2010模型选取中国南海北纬5° 东经110° 坐标点的地磁三分量数据40858.8、161.8、-3908.8 nT作为背景地磁场, 使模型进行水平360° 旋转, 并伴有横摇和纵倾动作, 且横摇和纵倾角不超过5° , 角度产生的最大测量误差为0.001° 。
然后在北纬5° 东经110° 到北纬6° 东经111° 海域均匀选取20个地磁场三分量数据作为背景地磁场, 假设仿真过程中载体坐标系与地理坐标系相同。为了更接近真实水下载体周围磁场环境(∑ HN≠ 0, HG未知), 进一步验证地磁场三分量差分测量方法的有效性, 在传感器1的仿真测量数据中加入最大幅值为± 25 nT随机误差来模拟实际环境中的外界杂散磁场, 设
采用最小二乘法拟合出每个传感器处感应系数矩阵和固定磁场的结果, 如表1所示。
表1中
通过图5可以看出, 本文提出的空间差分测量方法对地磁场三分量的测量误差明显低于传统STCM测量方法的误差。由于仿真中设置
而传统STCM方法受到外界杂散磁场的干扰比较大, 在每一个分量上产生的误差基本相同, 最大绝对值都在25 nT左右。说明空间差分方法能够更好地适应水下复杂磁场环境下地磁场分量的测量, 具有更好的抗干扰能力, 可为地磁辅助导航提供更加准确的地磁场信息。
为了方便实验, 载体模型采用长度80 cm、直径20 cm、管壁厚度1 cm的钢管制作(比仿真所用模型缩小了约50倍), 表面涂有灰色油漆, 重45 kg。两侧梯形尾翼对称焊接, 上底边长为10 cm, 下底边长为20 cm, 高度为20 cm, 厚度为1 cm。实验平台用木板制作, 长度为100 cm, 宽度为80 cm, 厚度为3 cm, 平台上固定有载体安装槽、4个Mag-03三分量传感器安装桩和数显倾角仪安装孔等, 此实验中使用的是对角的两个传感器安装桩, 高度差约为5 cm, 全部仪器采用无磁塑料螺丝安装和固定, 具体如图6所示。
主要进行了3个实验的测量:①载体实验平台差分参数估计实验; ②载体实验平台差分方法的测线验证实验; ③载体实验平台差分方法与STCM方法对比实验。
首先无载体时测量点
然后让载体实验平台沿着一条200 m测线进行直线运动, 共记录下11组三分量测量数据, 地磁场
最后保持载体实验平台不动, 让一辆金杯面包车模拟周围突然出现的外界杂散干扰磁场源, 从实验平台左侧3 m处驶过, 地磁场3个分量的测量结果如图8所示。
图8中地磁场分量真值通过无载体时传感器测量得到。通过图7可以看出当背景环境中磁场分布比较均匀, 无杂散干扰源时, 由于航向角的测量会引入误差使STCM传统测量方法和差分测量方法的地磁场测量结果均产生了一定的误差, 但该误差均小于100 nT。当外界出现未知杂散干扰磁源时, 如图8中(a)(b)(c)所示, STCM传统测量方法受到很强烈的干扰, 导致地磁场三分量测量误差随干扰磁源的出现明显增大, 最大值约为1000 nT, 而差分测量方法的测量误差仍然很平稳, 保持在150 nT以内。
水下地磁场的实地测量是实现高精度地磁辅助导航的关键技术。本文针对STCM传统测量方法使用单磁传感器在地磁场测量方面存在的抗干扰能力差、测量精度低的问题, 提出了基于多传感器的空间差分测量方法。能够有效去除载体自身干扰磁场以及外界磁异常体对地磁场测量造成的影响, 增强了测量过程中抵抗外界干扰磁场的能力。仿真结果表明:当外界干扰磁源在载体处产生的最大干扰磁场为25 nT, 且梯度为1 nT/m时, 空间差分测量方法的测量误差为(10± 5) nT, 比传统STCM测量方法的误差明显降低了1/2左右。地面实验结果表明:STCM传统方法产生的误差会随外界未知干扰磁场变化而变化, 而空间差分方法几乎不受外界干扰磁源的影响, 测量误差控制在自身稳定的范围内。因此, 通过上面的研究验证了空间差分方法的合理性和可行性, 为水下地磁辅助导航技术在实际中的应用奠定了基础。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|