留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于正弦拟合的三轴磁力仪标定方法

贾文抖 林春生 林朋飞 翟国君

贾文抖, 林春生, 林朋飞, 翟国君. 基于正弦拟合的三轴磁力仪标定方法[J]. 武汉大学学报 ( 信息科学版), 2019, 44(9): 1328-1333, 1406. doi: 10.13203/j.whugis20170379
引用本文: 贾文抖, 林春生, 林朋飞, 翟国君. 基于正弦拟合的三轴磁力仪标定方法[J]. 武汉大学学报 ( 信息科学版), 2019, 44(9): 1328-1333, 1406. doi: 10.13203/j.whugis20170379
JIA Wendou, LIN Chunsheng, LIN Pengfei, ZHAI Guojun. Calibration Method of Three-Axis Magnetometer Based on Sine Fitting[J]. Geomatics and Information Science of Wuhan University, 2019, 44(9): 1328-1333, 1406. doi: 10.13203/j.whugis20170379
Citation: JIA Wendou, LIN Chunsheng, LIN Pengfei, ZHAI Guojun. Calibration Method of Three-Axis Magnetometer Based on Sine Fitting[J]. Geomatics and Information Science of Wuhan University, 2019, 44(9): 1328-1333, 1406. doi: 10.13203/j.whugis20170379

基于正弦拟合的三轴磁力仪标定方法

doi: 10.13203/j.whugis20170379
基金项目: 

国家自然科学基金 41374018

国家自然科学基金 41374018

国家重大科研装备研制基金 ZDYZ2012-1

详细信息
    作者简介:

    贾文抖, 博士生, 主要从事磁异常数据解释、磁目标探测的理论与方法研究。jiawd163@163.com

    通讯作者: 林春生, 博士, 教授。lcs_and_zh@163.com
  • 中图分类号: P204

Calibration Method of Three-Axis Magnetometer Based on Sine Fitting

Funds: 

The National Natural Science Foundation of China 41374018

The National Natural Science Foundation of China 41374018

the National Major Scientific Research Equipment Development Fund Project of China ZDYZ2012-1

More Information
    Author Bio:

    JIA Wendou, PhD candidate, specializes in the theories and methods of magnetic anomaly data interpretation and the magnetic target detection. E-mail:jiawd163@163.com

    Corresponding author: LIN Chunsheng, PhD, professor. E-mail: lcs_and_zh@163.com
  • 摘要: 三轴磁力仪是一种常用的磁场测量工具,在地磁导航、海洋磁测、地磁勘探等领域应用较为广泛。作为一种磁场矢量测量工具,三轴磁力仪的磁测误差主要来源于三轴非正交、敏感轴灵敏度不一致及零偏误差等方面,磁力仪的标定优劣影响磁力仪的测量性能。为解决三轴磁力仪的标定问题,建立了磁力仪误差补偿模型,提出了基于正弦拟合的三轴磁力仪标定方法。将磁力仪沿平面旋转时的三轴输出数据拟合成正弦曲线,利用正弦曲线的幅值、初相角和平移量等信息计算出磁力仪的标定参数。仿真结果表明:该方法无需非线性优化过程即可实现标定参数的求解,得到的标定参数与设定值吻合度高;磁测噪声对三轴灵敏度系数及非正交角的计算结果影响较小。
  • 图  1  磁力仪三轴指向示意图

    Figure  1.  Diagram of Magnetometer Three-Axis Directions

    图  2  OXmYm面的旋转示意图

    Figure  2.  Diagram of OXmYm Plane Rotating

    图  3  OXmZm面的旋转示意图

    Figure  3.  Diagram of OXmZm Plane Rotating

    图  4  图 2测量模型中磁力仪三轴测量数据和拟合曲线

    Figure  4.  Three-Axis Measurement Data and Fitting Curves of Magnetometer Measurement Model in Fig. 2

    图  5  图 3测量模型中磁力仪x轴和z轴测量数据和拟合曲线

    Figure  5.  Measurement Data and Fitting Curves of x-Axis and z-Axis of Magnetometer Measurement Model in Fig. 3

    图  6  补偿前后的磁场强度曲线

    Figure  6.  Curves of Magnetic Intensity Before and After the Compensation

    图  7  测量数据中含有不同白噪声均值时的磁场强度补偿曲线

    Figure  7.  Magnetic Intensity Compensation Curves with Different White Noises in the Magnetic Measurement Data

    表  1  白噪声及其对应的标定参数

    Table  1.   White Noise and Its Corresponding Calibration Parameters

    噪声均值/nT α/(°) β/(°) γ/(°) kx ky kz Bx0/nT By0/nT Bz0/nT
    5 1.349 32.69 1.071 0.909 1.146 0 1.103 147.77 102.28 128.87
    15 1.346 32.64 1.070 0.909 1.146 0 1.103 158.16 112.43 138.99
    25 1.361 32.85 1.078 0.909 1.146 1 1.103 170.11 123.46 149.17
    40 1.371 32.74 1.082 0.909 1.145 8 1.103 287.23 236.19 266.54
    下载: 导出CSV
  • [1] 尹刚, 张英堂, 石志勇, 等.基于磁异常反演的磁航向误差实时补偿方法[J].武汉大学学报·信息科学版, 2016, 41(7):978-982 http://ch.whu.edu.cn/CN/abstract/abstract5494.shtml

    Yin Gang, Zhang Yingtang, Shi Zhiyong, et al. Real-Time Compensation Method of Magnetic Heading Perturbations Based on Magnetic Anomaly Inversion[J]. Geomatics and Information Science of Wuhan University, 2016, 41(7):978-982 http://ch.whu.edu.cn/CN/abstract/abstract5494.shtml
    [2] 郑晖, 王勇, 王虎彪, 等.地球重磁位场辅助水下潜艇导航仿真研究[J].武汉大学学报·信息科学版, 2012, 37(10):1198-1202 http://ch.whu.edu.cn/CN/abstract/abstract338.shtml

    Zheng Hui, Wang Yong, Wang Hubiao, et al. Simulation Research of Earth's Gravity and Geomagnetism Potential Field Aided Underwater[J]. Geomatics and Information Science of Wuhan University, 2012, 37(10):1198-1202 http://ch.whu.edu.cn/CN/abstract/abstract338.shtml
    [3] Birsan M. Recursive Bayesian Method for Magnetic Dipole Tracking with a Tensor Gradiometer[J]. IEEE Transactions on Magnetics, 2011, 47(2):409-415 doi:  10.1109/TMAG.2010.2091964
    [4] 董大群, 李斌.地磁测量三轴磁通门磁强计自适应校准[J].地球物理学报, 1992, 35(5):272-279

    Dong Daqun, Li Bin.Adaptive Calibration of Tri-axial Fluxgate Magnetometer for Geomagnetic Measurement[J].Chinese Journal of Geophysics, 1992, 35(5):272-279
    [5] 李勇, 刘文怡, 李杰, 等.基于椭球拟合的三轴磁传感器误差补偿方法[J].传感技术学报, 2012, 25(7):917-920 doi:  10.3969/j.issn.1004-1699.2012.07.010

    Li Yong, Liu Wenyi, Li Jie, et al.Error Compensation Method for Three-Axis Magnetic Sensor Based on Ellipsoid Fitting[J].Chinese Journal of Sensors and Actuators, 2012, 25(7):917-920 doi:  10.3969/j.issn.1004-1699.2012.07.010
    [6] Fang Jiancheng, Sun Hongwei, Cao Juanjuan, et al. A Novel Calibration Method of Magnetic Compass Based on Ellipsoid Fitting[J]. IEEE Transactions on Instrumentation & Measurement, 2011, 60(6):2053-2061
    [7] 龙达峰, 刘俊, 张晓明, 等.基于椭球拟合的三轴陀螺仪快速标定方法[J].仪器仪表学报, 2013, 34(6):100-106 http://d.old.wanfangdata.com.cn/Periodical/yqyb201306015

    Long Dafeng, Liu Jun, Zhang Xiaoming, et al. Triaxial Gyroscope Fast Calibration Method Based on Ellipsoid Fitting[J]. Chinese Journal of Scientific Instrument, 2013, 34(6):100-106 http://d.old.wanfangdata.com.cn/Periodical/yqyb201306015
    [8] Kok M, Schon T B. Maximum Likelihood Calibration of a Magnetometer Using Inertial Sensors[J]. Institute of Technology, 2014, 47(3):92-97 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CC0214921642
    [9] Gebre-Egziabher D, Elkaim G H, Powell J D, et al. A Non-linear, Two-step Estimation Algorithm for Calibrating Solid-state Strapdown Magnetometers[C]. The 8th International Conference on Navigation Systems, St. Petersburg, Russia, 2001
    [10] Foster C C, Elkaim G H. Extension of a Two Step Calibration Methodology to Include Nonorthogonal Sensor Axes[J]. IEEE Transactions on Aerospace and Electronic Systems, 2008, 44(3):1070-1078 doi:  10.1109/TAES.2008.4655364
    [11] 庞鸿锋, 潘孟春, 王伟, 等.基于高斯牛顿迭代算法的三轴磁强计校正[J].仪器仪表学报, 2013, 34(7):1506-1511 doi:  10.3969/j.issn.0254-3087.2013.07.010

    Pang Hongfeng, Pan Mengchun, Wang Wei, et al.Error Calibration of Three-Axis Magnetometer Based on Gauss-Newton Iteration Algorithm[J].Chinese Journal of Scientific Instrument, 2013, 34(7):1506-1511 doi:  10.3969/j.issn.0254-3087.2013.07.010
    [12] 吴德会, 黄松岭, 赵伟.基于FLANN的三轴磁强计误差校正研究[J].仪器仪表学报, 2009, 30(3):449-453 doi:  10.3321/j.issn:0254-3087.2009.03.001

    Wu Dehui, Huang Songling, Zhao Wei.Research on Correction of Tri-axial Magnetometer Based on FLANN[J].Chinese Journal of Scientific Instrument, 2009, 30(3):449-453 doi:  10.3321/j.issn:0254-3087.2009.03.001
    [13] Koo W, Sung S, Lee Y J. Error Calibration of Magnetometer Using Nonlinear Integrated Filter Model with Inertial Sensors[J]. IEEE Transactions on Magnetics, 2009, 45(6):2740-2743 doi:  10.1109/TMAG.2009.2020542
    [14] Pang Hongfeng, Li Ji, Chen Dixiang, et al. Calibration of Three-Axis Fluxgate Magnetometers with Nonlinear Least Square Method[J]. Measurement, 2013, 46(4):1600-1606 doi:  10.1016/j.measurement.2012.11.001
    [15] 张滢, 杨任农, 李明阳, 等.基于截断总体最小二乘算法的车载三轴磁力仪标定[J].兵工学报, 2015, 36(3):427-432 doi:  10.3969/j.issn.1000-1093.2015.03.007

    Zhang Ying, Yang Rennong, Li Mingyang, et al. Calibration of Vehicular Three-Axis Magnetometer via Truncated Total Least Squares Algorithm[J]. Acta Armamentarii, 2015, 36(3):427-432 doi:  10.3969/j.issn.1000-1093.2015.03.007
  • [1] 李宗春, 路志勇, 郭迎钢, 张冠宇, 何华, 冯其强, 陈少卿, 王俊威.  基于俯仰轴相交和球拟合的天线旋转中心测量方法 . 武汉大学学报 ( 信息科学版), 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064
    [2] 薛武, 张永生, 王涛, 于英, 曹彬才.  无人直升机面阵影像高精度对地定位 . 武汉大学学报 ( 信息科学版), 2019, 44(2): 246-253. doi: 10.13203/j.whugis20170075
    [3] 唐新明, 陈继溢, 李国元, 高小明, 张文君.  资源三号02星激光测高误差分析与指向角粗标定 . 武汉大学学报 ( 信息科学版), 2018, 43(11): 1611-1619. doi: 10.13203/j.whugis20170058
    [4] 潘国荣, 郭巍, 周跃寅.  一种基于先验误差分解定权的精密测量拟合算法 . 武汉大学学报 ( 信息科学版), 2015, 40(10): 1339-1343. doi: 10.13203/j.whugis20130768
    [5] 谌一夫, 刘璐, 张春玲.  ZY 3卫星在轨几何标定方法 . 武汉大学学报 ( 信息科学版), 2013, 38(5): 557-560.
    [6] 李清泉, 陈小宇, 毛庆洲, 李必军.  光纤陀螺仪关键参数快速标定方法 . 武汉大学学报 ( 信息科学版), 2013, 38(11): 1369-1373.
    [7] 卞光浪, 翟国君, 黄谟涛, 欧阳永忠.  利用剖面测线磁场数据反演三度磁性体参数 . 武汉大学学报 ( 信息科学版), 2012, 37(1): 91-95.
    [8] 于波, 黄谟涛, 翟国君, 刘雁春.  海洋磁力测量空间归算阈值条件确定及其应用 . 武汉大学学报 ( 信息科学版), 2010, 35(2): 172-175.
    [9] 邓洪涛.  地球内核与磁场 . 武汉大学学报 ( 信息科学版), 2010, 35(7): 854-856.
    [10] 李清泉, 雷波, 毛庆洲, 傅智能.  利用激光三角法进行快速车辙检测 . 武汉大学学报 ( 信息科学版), 2010, 35(3): 302-307.
    [11] 詹总谦, 张祖勋, 张剑清.  基于LCD的相机标定精度及其误差分析 . 武汉大学学报 ( 信息科学版), 2008, 33(11): 1142-1145.
    [12] 于波, 刘雁春, 边刚, 肖付民.  海洋磁力测量垂直空间改正研究 . 武汉大学学报 ( 信息科学版), 2008, 33(9): 926-929.
    [13] 张卡, 盛业华, 叶春, 梁诚.  车载三维数据采集系统的绝对标定及精度分析 . 武汉大学学报 ( 信息科学版), 2008, 33(1): 55-59.
    [14] 王解先.  工业测量中一种二次曲面的拟合方法 . 武汉大学学报 ( 信息科学版), 2007, 32(1): 47-50.
    [15] 于波, 刘雁春, 肖付民, 边刚.  海洋磁力测量中确定最佳拖曳距离的一种方法 . 武汉大学学报 ( 信息科学版), 2006, 31(11): 999-1002.
    [16] 黄谟涛, 翟国君, 欧阳永忠, 任来平.  海洋磁力测量误差补偿技术研究 . 武汉大学学报 ( 信息科学版), 2006, 31(7): 603-606.
    [17] 周云耀, 蔡亚先.  甚宽频带地震计频域标定的相关处理技术——用正交相关法提高大噪声条件下正弦波标定的精度 . 武汉大学学报 ( 信息科学版), 2005, 30(7): 632-635.
    [18] 冯文灏, 李欣.  基于三旋转自由度激光经纬仪的结构光工程测量原理与应用 . 武汉大学学报 ( 信息科学版), 1998, 23(4): 345-350,393.
    [19] 李德仁, 王新华.  CCD阵列相机的几何标定 . 武汉大学学报 ( 信息科学版), 1997, 22(4): 308-313,317.
    [20] 杨丹, 吴世明, 吴小玲.  XH-100多功能三轴仪与花岗岩残积土动力特性研究 . 武汉大学学报 ( 信息科学版), 1991, 16(3): 53-58.
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  1212
  • HTML全文浏览量:  158
  • PDF下载量:  84
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-06-05
  • 刊出日期:  2019-09-05

基于正弦拟合的三轴磁力仪标定方法

doi: 10.13203/j.whugis20170379
    基金项目:

    国家自然科学基金 41374018

    国家自然科学基金 41374018

    国家重大科研装备研制基金 ZDYZ2012-1

    作者简介:

    贾文抖, 博士生, 主要从事磁异常数据解释、磁目标探测的理论与方法研究。jiawd163@163.com

    通讯作者: 林春生, 博士, 教授。lcs_and_zh@163.com
  • 中图分类号: P204

摘要: 三轴磁力仪是一种常用的磁场测量工具,在地磁导航、海洋磁测、地磁勘探等领域应用较为广泛。作为一种磁场矢量测量工具,三轴磁力仪的磁测误差主要来源于三轴非正交、敏感轴灵敏度不一致及零偏误差等方面,磁力仪的标定优劣影响磁力仪的测量性能。为解决三轴磁力仪的标定问题,建立了磁力仪误差补偿模型,提出了基于正弦拟合的三轴磁力仪标定方法。将磁力仪沿平面旋转时的三轴输出数据拟合成正弦曲线,利用正弦曲线的幅值、初相角和平移量等信息计算出磁力仪的标定参数。仿真结果表明:该方法无需非线性优化过程即可实现标定参数的求解,得到的标定参数与设定值吻合度高;磁测噪声对三轴灵敏度系数及非正交角的计算结果影响较小。

English Abstract

贾文抖, 林春生, 林朋飞, 翟国君. 基于正弦拟合的三轴磁力仪标定方法[J]. 武汉大学学报 ( 信息科学版), 2019, 44(9): 1328-1333, 1406. doi: 10.13203/j.whugis20170379
引用本文: 贾文抖, 林春生, 林朋飞, 翟国君. 基于正弦拟合的三轴磁力仪标定方法[J]. 武汉大学学报 ( 信息科学版), 2019, 44(9): 1328-1333, 1406. doi: 10.13203/j.whugis20170379
JIA Wendou, LIN Chunsheng, LIN Pengfei, ZHAI Guojun. Calibration Method of Three-Axis Magnetometer Based on Sine Fitting[J]. Geomatics and Information Science of Wuhan University, 2019, 44(9): 1328-1333, 1406. doi: 10.13203/j.whugis20170379
Citation: JIA Wendou, LIN Chunsheng, LIN Pengfei, ZHAI Guojun. Calibration Method of Three-Axis Magnetometer Based on Sine Fitting[J]. Geomatics and Information Science of Wuhan University, 2019, 44(9): 1328-1333, 1406. doi: 10.13203/j.whugis20170379
  • 三轴磁力仪常被用在地磁导航、磁目标追踪、未爆物探测等方面[1-3],但受制造安装工艺、敏感轴电气性不一致及零点偏移等因素影响,磁力仪本身存在三轴非正交、敏感轴灵敏度不一致及零偏误差等性能缺陷,使得磁测的准确性受到较大影响。例如,在地磁场环境中,如果磁力仪三轴正交误差达到1°,则由正交误差所引起的磁场测量误差将达数百nT之多[4]。为获得准确的磁场数据,在使用之前有必要对三轴磁力仪进行相应的标定。

    椭球模型[5-7]是人们开展三轴磁力仪标定工作时常用的一种数学模型,建立该模型的理论依据是在均匀的磁场环境中,理想三轴磁力仪在各姿态下所测得的总磁场强度对应到空间方位角时是呈球面分布的,但由于受到磁力仪三轴非正交、敏感轴灵敏度不一致及零偏误差的影响,未标定过的磁力仪测得的总磁场强度对应到空间方位角时是呈椭球面分布的。若能得到椭球面的位置和形状参数,即可实现对磁力仪标定参数的求解,从而实现磁测性能的校正。此外,还可以借助光泵等高精度磁测工具测得当地磁场值作为参考量,利用非线性优化方法实现对标定参数的求解。在以上两种不同的求解模型的基础上,人们研究了很多参数求解方法,例如最大似然估计法[8]、两步估计法[9-10]、高斯-牛顿法[11]、神经网络法[12]、卡尔曼滤波法[13]、非线性最小二乘法[14-15]等参数优化方法均可实现对磁力仪标定参数的求解。但这些优化方法多需要通过迭代计算来完成参数的非线性优化求解,计算过程较为复杂。

    为避免求解磁力仪标定参数时进行非线性优化过程,本文提出了一种基于正弦拟合的三轴磁力仪标定方法,将三轴磁力仪沿平面旋转时的三轴输出数据拟合成正弦曲线,利用正弦曲线的幅值、初相角和平移量直接计算标定参数,从而实现对磁力仪的标定和三轴磁测数据的补偿。

    • 图 1所示,建立直角坐标系O-XYZ,未经校准的三轴磁力仪的Xm轴与坐标系的X轴重合,Ym轴在OXY平面内,Ym轴与Y轴的夹角为αZm轴与XXm)轴的夹角记为φZm轴与Z轴的夹角记为γZm轴在OXY面上的投影与X轴的夹角记为β

      图  1  磁力仪三轴指向示意图

      Figure 1.  Diagram of Magnetometer Three-Axis Directions

      在磁场强度为B的环境中,理想情况下,三轴相互正交的磁力仪输出为$B=\left( {{B}_{x}},{{B}_{y}},{{B}_{z}} \right)$,由于磁力仪三轴非正交、各轴灵敏度不一致以及各轴零偏误差的存在,三轴磁测输出可表示为:

      $$ \left[ \begin{matrix} {{B}_{xm}} \\ {{B}_{ym}} \\ {{B}_{zm}} \\ \end{matrix} \right]=\left[ \begin{matrix} {{k}_{x}} & 0 & 0 \\ {{k}_{y}}\text{sin}\alpha & {{k}_{y}}\text{cos}\alpha & 0 \\ {{k}_{z}}\text{cos}\beta \text{sin}\gamma & {{k}_{z}}\text{sin}\beta \text{sin}\gamma & {{k}_{z}}\text{cos}\gamma \\ \end{matrix} \right]\left[ \begin{matrix} {{B}_{x}} \\ {{B}_{y}} \\ {{B}_{z}} \\ \end{matrix} \right]+\left[ \begin{matrix} {{B}_{x0}} \\ {{B}_{y0}} \\ {{B}_{z0}} \\ \end{matrix} \right] $$ (1)

      式中,αβγ是用来描述磁力仪三轴非正交角的参数;kxkykz为描述磁力仪三轴敏感度的灵敏度系数;Bx0By0Bz0为磁力仪三轴零偏误差。根据式(1)有:

      $$ \left[ \begin{matrix} {{B}_{x}} \\ {{B}_{y}} \\ {{B}_{z}} \\ \end{matrix} \right]=\left[ \begin{matrix} {{k}_{11}} & {{k}_{12}} & {{k}_{13}} \\ {{k}_{21}} & {{k}_{22}} & {{k}_{23}} \\ {{k}_{31}} & {{k}_{32}} & {{k}_{33}} \\ \end{matrix} \right]\left[ \begin{matrix} {{B}_{xm}}-{{B}_{x0}} \\ {{B}_{ym}}-{{B}_{y0}} \\ {{B}_{zm}}-{{B}_{z0}} \\ \end{matrix} \right] $$ (2)

      式中,系数kiji, j=1, 2, 3)为中间过程计算参数,可由非正交角和灵敏度系数联合计算得到。计算得到式(1)中的αβγkxkykzBx0By0Bz0这9个标定参数,即可以实现对三轴磁力仪的标定及磁测数据的补偿。

    • 在均匀的磁场环境中,将三轴磁力仪固定在无磁旋转平台上,使面$O{{X}_{m}}{{Y}_{m}}$平行于水平面(见图 2),等角度间隔转动旋转平台,对应不同的旋转角取一组磁测数据$\left( {{B}_{xmi}},{{B}_{ymi}},{{B}_{zmi}} \right)$,记磁力仪三轴输出数据为$\left[ {{B}_{xm1}},{{B}_{xm2}}\cdots \right]$、$\left[ {{B}_{ym1}},{{B}_{ym2}}\cdots \right]$、$\left[ {{B}_{zm1}},{{B}_{zm2}}\cdots \right]$。

      图  2  OXmYm面的旋转示意图

      Figure 2.  Diagram of OXmYm Plane Rotating

      B在水平面上的投影为${{B}_{xy}}=B\text{cos}I$(I是地磁倾角),理想情况下,XmYm轴的输出应满足曲线:

      $$ \left\{ \begin{matrix} {{B}_{xm}}={{B}_{xy}}\text{sin}\left( \theta +{{\theta }_{x}} \right) \\ {{B}_{ym}}={{B}_{xy}}\text{sin}\left( \theta +{{\theta }_{y}} \right) \\ \end{matrix} \right. $$ (3)

      式中,Bxy为磁场BOXY平面上的投影;θ为转动角;对相位角θxθy有|θx-θy |=90°。

      受磁力仪三轴非正交、灵敏度不一致及零偏误差的影响,实际上XmYm轴的输出分别对应曲线:

      $$ \left\{ \begin{array}{*{35}{l}} {{B}_{xm}}={{B}_{xmy}}\text{sin}\left( \theta +{{\theta }_{x0}} \right)+{{B}_{xm0}} \\ {{B}_{ym}}={{B}_{ymx}}\text{sin}\left( \theta +{{\theta }_{y0}} \right)+{{B}_{ym0}} \\ \end{array} \right. $$ (4)

      式中,BxymBymx分别为XmYm轴的输出曲线幅值;θx0θy0为曲线的初相位;Bxm0Bym0为曲线的平移量。相位角θx0θy0的差值即是磁力仪Xm轴和Ym轴之间的夹角,曲线位置的平移量Bxm0Bym0即是磁力仪Xm轴和Ym轴的零偏误差。Xm轴和Ym轴之间的非正交角为:

      $$ \alpha =\left| {{\theta }_{y0}}-{{\theta }_{x0}}-90{}^\circ \right| $$ (5)

      Xm轴和Ym轴的零偏误差分别为:

      $$ \left\{ \begin{array}{*{35}{l}} {{B}_{x0}}={{B}_{xm0}} \\ {{B}_{y0}}={{B}_{ym0}} \\ \end{array} \right. $$ (6)

      对于Xm轴和Ym轴的灵敏度系数kxky,其比例关系为:

      $$ {{k}_{yx}}=\frac{{{k}_{y}}{{B}_{xy}}}{{{k}_{x}}{{B}_{xy}}}=\frac{{{B}_{xmy}}}{{{B}_{ymx}}} $$ (7)

      图 1中的夹角β可根据Xm轴、Zm轴对应拟合曲线的初始相位角θz0θx0计算得到。若Zm轴的输出Bzm保持不变,则Zm轴是垂直于OXm Ym平面的,此时β=γ=0;若Bzm的大小随转动发生变化,根据式(1)有:

      $$ \begin{matrix} {{B}_{zm}}={{k}_{z}}\left( {{B}_{x}}\text{cos}\beta +{{B}_{y}}\text{sin}\beta \right)\text{sin}\gamma +{{k}_{z}}{{B}_{z}}\text{cos}\gamma +{{B}_{z0}}= \\ {{k}_{z}}\left( {{B}_{xy}}\text{sin}\left( \theta +{{\theta }_{x0}} \right)\text{cos}\beta +{{B}_{xy}}\text{cos}\left( \theta +{{\theta }_{x0}} \right)\text{sin}\beta \right)\text{sin}\gamma + \\ {{k}_{z}}{{B}_{z}}\text{cos}\gamma +{{B}_{z0}}= \\ {{k}_{z}}{{B}_{xy}}\text{sin}\left( \theta +{{\theta }_{x0}}+\beta \right)\text{sin}\gamma +{{k}_{z}}{{B}_{z}}\text{cos}\gamma +{{B}_{z0}}= \\ A\text{sin}\left( \theta +{{\theta }_{x0}}+\beta \right)+{{C}_{0}}= \\ A\text{sin}\left( \theta +{{\theta }_{z0}} \right)+{{C}_{0}} \\ \end{matrix} $$ (8)

      Zm轴的输出Bzm在OXm Ym平面上的投影为Bzmxy=Bzm sinγ,相当于只对Bzm的大小作压缩变化,不改变其初相位角θz0的大小,则Zm轴在OXm Ym面上的投影与Xm轴的夹角为:

      $$ \beta =\left| {{\theta }_{z0}}-{{\theta }_{x0}} \right| $$ (9)

      图 3所示,翻转磁力仪,将OXmZm平面水平固定在无磁旋转平台上,按照同样的方法旋转记录磁测数据,记此时的三轴输出为(Bxm2, Bym2, Bzm2)。类似地,拟合得到Xm轴、Zm轴的输出曲线为:

      图  3  OXmZm面的旋转示意图

      Figure 3.  Diagram of OXmZm Plane Rotating

      $$ \left\{ \begin{array}{*{35}{l}} {{B}_{xm2}}={{B}_{xmz}}\text{sin}\left( \theta +\theta {{\text{ }\!\!'\!\!\text{ }}_{x0}} \right)+{{B}_{xm0}} \\ {{B}_{zm2}}={{B}_{zmx}}\text{sin}\left( \theta +\theta {{\text{ }\!\!'\!\!\text{ }}_{z0}} \right)+{{B}_{zm0}} \\ \end{array} \right. $$ (10)

      式中,BxmzBzmx分别为Xm轴、Zm轴的输出曲线的幅值;θ'x0θ'z0为曲线初相位;Bxm0Bzm0为曲线平移量。则Zm轴的零偏误差为:

      $$ {{B}_{z0}}={{B}_{zm0}} $$ (11)

      Xm轴与Zm轴的夹角为:

      $$ \varphi =\left| \theta {{\text{ }\!\!'\!\!\text{ }}_{x0}}-\theta {{\text{ }\!\!'\!\!\text{ }}_{z0}} \right| $$ (12)

      结合图 1,有cosφ=cosβsinγ,则:

      $$ \gamma =\text{arcsin}\left( \text{cos}\varphi /\text{cos}\beta \right) $$ (13)

      Xm轴、Zm轴的灵敏度系数kxkz之间的比例关系为:

      $$ {{k}_{zx}}=\frac{{{k}_{z}}{{B}_{xz}}}{{{k}_{x}}{{B}_{xz}}}=\frac{{{B}_{zmx}}}{{{B}_{xmz}}} $$ (14)

      由式(1)可得:

      $$ \left\{ \begin{array}{*{35}{l}} {{k}_{x}}{{B}_{x}}={{B}_{xm}}-{{B}_{x0}} \\ {{k}_{y}}{{B}_{y}}=\left( {{B}_{ym}}-{{B}_{y0}}-{{k}_{yx}}{{k}_{x}}{{B}_{x}}\text{sin}\alpha \right)/\text{cos}\alpha \\ {{k}_{z}}{{B}_{z}}=\begin{array}{*{35}{l}} ({{B}_{zm}}-{{B}_{z0}}-{{k}_{zx}}{{k}_{x}}{{B}_{x}}\text{cos}\beta \text{sin}\gamma - \\ \left( {{k}_{zx}}/{{k}_{yx}} \right){{k}_{y}}{{B}_{y}}\text{sin}\beta \text{sin}\gamma )/\text{cos}\gamma \\ \end{array} \\ \end{array} \right. $$ (15)

      此时,式(15)等式右侧的数据已经能够计算得到,再联立式(7)、式(14)和式(15),可得Xm轴的灵敏度系数为:

      $$ {{k}_{x}}=\sqrt{\left( {{\left( {{k}_{x}}{{B}_{x}} \right)}^{2}}+{{\left( \frac{{{k}_{y}}{{B}_{y}}}{{{k}_{yx}}} \right)}^{2}}+{{\left( \frac{{{k}_{z}}{{B}_{z}}}{{{k}_{zx}}} \right)}^{2}} \right)/{{B}^{2}}} $$ (16)

      kx代入式(7)和式(14)即可计算出kykz。此时,所需的9个标定参数已全部计算得到,标定后的磁力仪三轴输出可修正为:

      $$ \left\{ \begin{array}{*{35}{l}} {{B}_{x}}=\frac{{{B}_{xm}}-{{B}_{x0}}}{{{k}_{x}}} \\ {{B}_{y}}=\frac{{{B}_{ym}}-{{B}_{y0}}-{{k}_{y}}{{B}_{x}}\text{sin}\alpha }{{{k}_{y}}\text{cos}\alpha } \\ {{B}_{z}}=\frac{{{B}_{zm}}-{{B}_{z0}}-{{k}_{z}}{{B}_{x}}\text{cos}\beta \text{sin}\gamma -{{k}_{z}}{{B}_{y}}\text{sin}\beta \text{sin}\gamma }{{{k}_{z}}\text{cos}\gamma } \\ \end{array} \right. $$ (17)
    • 设磁力仪三轴非正交角分别为α=1.35°、β=32.7°,γ=1.07°,三轴灵敏度系数分别为kx=0.909,ky=1.146,kz=1.103,三轴零偏误差分别为Bx0=143 nT,By0=97.4 nT,Bz0=123.7 nT,地磁场强度B=50 000 nT,地磁倾角I=43°。

    • 将磁力仪的OXmYm面水平固定在旋转平面上,转动旋转平台记录磁力仪的三轴输出数据,并将三轴磁测数据拟合成Asin(θ+C)+D(其中,A为曲线幅值,θ为转动角,C为初相位,D为曲线上下偏移量)的曲线形式,拟合后的三轴输出曲线如图 4所示。

      图  4  图 2测量模型中磁力仪三轴测量数据和拟合曲线

      Figure 4.  Three-Axis Measurement Data and Fitting Curves of Magnetometer Measurement Model in Fig. 2

      图 4可以看出,测量数据点基本都分布在拟合曲线上,两者的吻合度较好。对应的XmYmZm三轴输出数据的拟合曲线描述方程为:

      $$ \left\{ \begin{align} & {{B}_{xm1}}=33\ 240.0\text{sin}\left( \theta +90{}^\circ \right)+143.0 \\ & {{B}_{ym1}}=41\ 906.6\text{sin}\left( \theta +178.65{}^\circ \right)+97.4 \\ & {{B}_{zm1}}=753.2\sin \left( \theta +122.7{}^\circ \right)+37\ 729.4 \\ \end{align} \right. $$

      根据式(5)~(7)、式(9)可计算得到磁力仪的标定参数分别为α=1.35°,β=32.7°,Bx0=143.0 nT,By0=97.4 nT,kyx =1.260 7。

      改变三轴磁力仪的放置姿态,使OXmZm面水平固定于旋转平台上,转动平台,记录磁测数据,并将Xm轴和Zm轴的输出数据拟合成相应的正弦曲线,如图 5所示。

      图  5  图 3测量模型中磁力仪x轴和z轴测量数据和拟合曲线

      Figure 5.  Measurement Data and Fitting Curves of x-Axis and z-Axis of Magnetometer Measurement Model in Fig. 3

      图 5可以看出,拟合得到的曲线与原始测量数据的吻合度较好,测量数据基本都分布在拟合曲线上。对应的Xm轴、Zm轴输出数据的拟合曲线描述方程为:

      $$ \left\{ \begin{align} & {{B}_{xm2}}=33\ 240.0\text{sin}\left( \theta +90{}^\circ \right)+143.0 \\ & {{B}_{zm2}}=40\ 334.2\text{sin}\left( \theta +0.900\ 4{}^\circ \right)+123.7 \\ \end{align} \right. $$

      根据式(11)~(14),可计算得到磁力仪的标定参数分别为φ=89.999 6°,γ=1.07°,Bz0=123.7 nT,kzx=1.213 4。

      结合磁场强度B,再联立式(15)、式(16)可计算得到灵敏度系数kx=0.909,ky=1.146,kz=1.103。此时,所需的9个标定参数已全部得到。

    • 在得到磁力仪的标定参数后,利用式(17)对测得的磁场三分量进行补偿,补偿后的磁场三分量为BxByBz。以下分析测量数据中不含和含有测量噪声时的标定参数求解结果及磁场强度补偿结果变化情况。

    • 利用§3.1中的仿真数据进行计算,在§3.1中,三轴磁力仪的标定参数已计算得到,利用补偿后的三轴输出数据进行磁场强度的计算,补偿前后的计算结果如图 6所示。

      图  6  补偿前后的磁场强度曲线

      Figure 6.  Curves of Magnetic Intensity Before and After the Compensation

      图 6可知,在均匀的磁场环境中,当标定前的磁力仪沿平面转动时,根据测量数据直接计算得到的磁场强度变化幅度很大,磁场强度曲线有较大的波动起伏,偏离真实值可达到5 650 nT;经补偿后,利用三轴输出数据计算得到的磁场强度输出值较为稳定,补偿后的磁场强度曲线基本呈一条直线,补偿后的磁场强度与真实的磁场强度偏差很小。

    • 在仿真数据中添加白噪声用来模拟测量数据,当测量数据中白噪声均值分别为5 nT、15 nT、25 nT、40 nT时,利用上述方法进行标定参数的求解,计算得到的标定参数见表 1。当仿真数据中含有不同均值的白噪声时,补偿后的磁场强度曲线见图 7

      表 1  白噪声及其对应的标定参数

      Table 1.  White Noise and Its Corresponding Calibration Parameters

      噪声均值/nT α/(°) β/(°) γ/(°) kx ky kz Bx0/nT By0/nT Bz0/nT
      5 1.349 32.69 1.071 0.909 1.146 0 1.103 147.77 102.28 128.87
      15 1.346 32.64 1.070 0.909 1.146 0 1.103 158.16 112.43 138.99
      25 1.361 32.85 1.078 0.909 1.146 1 1.103 170.11 123.46 149.17
      40 1.371 32.74 1.082 0.909 1.145 8 1.103 287.23 236.19 266.54

      图  7  测量数据中含有不同白噪声均值时的磁场强度补偿曲线

      Figure 7.  Magnetic Intensity Compensation Curves with Different White Noises in the Magnetic Measurement Data

      表 1可以看出,3种补偿系数的求解均受测量数据中白噪声的影响。仿真测量数据中含有不同均值的白噪声时,随着测量数据中白噪声均值的增加,磁力仪三轴非正交角、灵敏度系数和零偏误差与真实值之间的偏差均有不同程度的增大。其中,三轴非正交角和灵敏度系数受测量白噪声均值的影响相对较小,零偏误差随白噪声均值的增加而变大的幅度较为明显。

      图 7可以看出,当测量数据中含有白噪声并进行磁力仪的标定参数计算时,利用计算得到的标定参数进行磁场数据的补偿,补偿后的磁场强度与磁场真实值之间有一定的偏差,补偿后的磁场曲线是一条在磁场真实值附近做上下浮动的曲线,白噪声均值不同,补偿后的磁场强度曲线的波动程度也不同。随着白噪声均值由5 nT增加到40 nT,曲线浮动变化的幅值也在增大。

    • 本文提出了一种基于正弦拟合的三轴磁力仪标定方法,避免了标定参数求解时复杂的非线性优化过程。利用三轴输出的正弦曲线幅值、初相角和上下平移量等信息对磁力仪的非正交角、灵敏度系数和零偏误差等标定参数进行求解,实现了对磁力仪的标定及三轴输出数据的补偿。仿真结果表明,利用该方法可有效实现磁力仪标定参数的求解,标定后的磁力仪磁测性能改善明显,磁测噪声对三轴灵敏度系数及非正交角的反演结果影响较小。

参考文献 (15)

目录

    /

    返回文章
    返回