摘要
面向深空磁场探测应用的矢量原子磁力仪,不仅需要在弱场下具备较高的灵敏度,还需要具有长期准确测量能力。故需要对矢量原子磁力仪在近零场下的灵敏度和线性度,尤其是准确度方面的性能指标进行地面测试标定。由于目前国家标准计量机构暂时没有可对微弱磁场强度下的磁力仪进行直接测试标定的装置,故提出了一种近零场下间接的测试标定装置和方法。在磁屏蔽筒内利用经过计量标定的CPT原子磁力仪对磁场发生线圈进行标定,通过标定过的线圈产生微弱人工磁场,结合零点旋转测试标定方法,实现对磁力仪准确度的测试标定。完成了近零场高精度原子磁力仪传感器的测试标定和不确定度分析,在±5 nT量程范围内线性度为0.03%,零点合成不确定度为12.90 pT,测量值合成不确定度为15.80 pT。
1 引 言
矢量原子磁力仪是一种典型的磁场探测仪器,可应用于生物医
国内目前已有的弱磁场下的磁力仪测试标定方法基本分为2种,首先都需要屏蔽地磁场的影响,为磁力仪提供1个近零磁场的空
对于面向深空磁场探测应用的矢量磁力仪,需要在地面完成对其灵敏度和准确度等指标的测试标定,故本文提出了1种间接测试标定高精度矢量原子磁力仪的方法。首先利用多层磁屏蔽筒屏蔽地球磁场与外界扰动,在筒内近零磁环境下利用经过中国计量科学研究院计量标定的CPT原子磁力
2 试验装置与方法
2.1 被测矢量原子磁力仪装置
被测矢量原子磁力仪(vector atomic magnetometer,VAM)系统示意图如

图1 被测矢量原子磁力仪系统示意图
Fig.1 Schematic diagram of the measured vector atomic magnetometer system
直流共振光吸收信号噪声较大,利用线圈驱动模块控制产生交流调制磁场,经过调制的共振光信号通过光纤传输到光电检测模块。通过基于现场可编程逻辑门阵列(field programmable gate array,FPGA)的控制系统中的锁相放大器对其进行调制解调产生高信噪比的共振色散信号,与锁定信号零点比较后进入闭环控制器,驱动线圈产生补偿磁场,通过对补偿磁场的精确测量即可得到被测外界矢量磁场的大小。
2.2 测试标定装置
对矢量原子磁力仪测试标定的过程在磁屏蔽筒内完成。如

图2 测试标定装置图
Fig.2 Diagram of testing and calibration system device
首先需要通过CPT原子磁力仪在2 000 nT量程下对磁场线圈系数进行标定。所使用的CPT原子磁力仪已经在中国计量科学院完成了测试标定,标定过程中以准确度为0.1 nT的Overhauser磁力仪作为计量基准装置。综合CPT原子磁力仪的线性度偏差、标准比对不确定度和重复测量不确定度数据得到其总不确定度为0.11 nT。
电流源驱动线圈产生磁场,CPT原子磁力仪探头放置在载物台中心处进行测量。通过对CPT原子磁力仪探测的磁场数据和数字多用表测得的电流值进行线性拟合,即可得到人工磁场线圈的系数。
然后利用标定过的磁场线圈产生±5 nT小量程内的测量磁场,用于测试标定实验室研制的弱场矢量原子磁力仪系统。
3 不确定度评估
3.1 校准场的不确定度
首先利用标定过的CPT原子磁力仪、电流源和数字多用表在大量程下(2 000 nT)对人工磁场线圈的系数进行标定,由
(1) |
拟合计算得到=25.354 2 nT/mA,对大量程下的线圈系数,其不确定度来源主要包括:
1) CPT原子磁力仪的不确定度
根据中国计量科学研究院测试标定结果,CPT原子磁力仪的不确定度为0.11 nT,故=0.11 nT。
2) 电流值的不确定度
线圈系数标定过程中读取电流值的数字多用表在100 mA档位下工作,该档位下示值不确定度为U=0.006 0 mA,线性度为60×1
根据不确定度的传播律,的不确定度可由
(2) |
即
=0.024 nT=24 pT
然后利用标定过的人工磁场线圈在电流源的驱动下产生小量程(±5 nT)范围内的校准磁场,如
= | (3) |
其中。对于小量程下线圈系数,其不确定度来源主要包括:
1) 大量程下磁场系数引入的不确定度
根据上述大量程线圈系数的不确定度24 pT,下放到小量程下得到=0.06 pT。
2) 线圈的线性度
人工磁场线圈系数的最大非线性偏差为-0.934 nT,故线性度为468×1
3) 磁屏蔽筒开合引起的不确定度
测试标定过程中磁屏蔽筒不免进行多次开合,开合过程中机械震动,外界扰动等原因也会在一定程度上对人工线圈产生的磁场有一定的影响。利用磁通门磁力仪在线圈产生的一定磁场环境下进行多次开合实验,记录每次开合关闭屏蔽筒盖后测得的磁场数据如
测试序次 | 磁场数据x |
---|---|
1 | 2 006.7 |
2 | 2 006.6 |
3 | 2 006.6 |
4 | 2 006.6 |
5 | 2 006.5 |
6 | 2 006.5 |
7 | 2 006.6 |
8 | 2 006.6 |
9 | 2 006.5 |
10 | 2 006.5 |
多次测量结果的算数平均值为=2 006.57 nT。由贝塞尔
s | (4) |
计算得到s=0.004 556 nT=4.556 pT
故在2 000 nT下多次屏蔽筒盖开合产生的磁场不确定度为 4.556/ pT=1.44 pT。若将其考虑为对人工磁场线圈系数的影响,相应地在小量程下屏蔽筒开合产生的磁场不确定度为0.004 pT。
则小量程下线圈系数的不确定度可由
= | (5) |
即
= pT=2.34 pT
对于校准磁场,其不确定度来源主要有:
1) 线圈产生磁场带来的不确定度
小量程下线圈系数的不确定度。由上述计算可知=2.34 pT。
电流值的不确定度。在小量程范围内线圈产生校准磁场,数字多用表在小档位1 mA下工作,该档位下的不确定度为U=0.000 049 mA,线性度为49×1
则利用
(6) |
计算得到=6.636 pT。
2) 磁场线圈梯度引起的不确定度
磁力仪探头放置位置不同测得的线圈产生的磁场强度略有不同,磁场线圈在X轴向上引起的非均匀性为0.4%@80 mm。传感器探头敏感单元大小为15 mm×25 mm,则线圈非均匀性在小量程范围内偏移25mm引起的磁场绝对误差为6 pT,考虑到偏移距离的概率分布,磁场梯度所引起的磁场不确定度= pT=3.464 pT。
3) 筒内剩磁不确定度
可由矢量原子磁力仪的灵敏度表征,=1 pT。
综上所述,的合成不确定度利用
(7) |
即
=7.55 pT
对上述的不确定度来源总结如
不确定度来源 | 不确定度 |
---|---|
小量程下线圈系数的不确定度 | 2.34 |
电流示值不确定度 | 0.245 |
磁场梯度引起的不确定度 | 3.464 |
剩磁的不确定度 | 1 |
的合成不确定度 | 7.55 |
除此之外,屏蔽筒内存在的剩余磁场大小会影响磁力仪测得的实际校准磁场大小,但对矢量原子磁力仪进行指标测试时会首先进行零点的旋转标定,以消除筒内剩余磁场偏置的影响。
3.2 待测矢量原子磁力仪零点偏差引起的不确定度
实验室已搭建完成一套完整的弱场矢量原子磁力仪系统,可在磁屏蔽筒内实现磁场矢量测量。矢量磁力仪的测量结果包括传感器零点位置误差,所以首先需要对磁力仪的零点进行测试标定。标定前关闭磁屏蔽筒盖,传感器探头放置在载物台的中央,通过外部驱动完成指定角度旋转,无磁转台旋转位置精度为27.5″,整个实验过程中不再需要进行磁屏蔽筒的开合,可保持该次测试过程中磁场环境稳定。
弱场矢量原子磁力仪传感器探头经过微调旋转,通过极值确定磁轴位置,该初始位置记为0°。筒内存在一定的剩余磁场,此时探头测量外部磁场=,矢量磁力仪测量磁场时本身存在零点偏差。在0°位置测得磁场值如
=[()+()]/2= | (8) |

图3 传感器探头0°和180°位置磁场示意图
Fig.3 Schematic diagram of magnetic field at 0 ° and 180 ° positions of sensor head
随后在筒内剩余磁场的基础上,利用人工磁场线圈产生磁场,此时磁力仪探头测量外部磁场=+。在磁力仪量程范围内完成7组实验,测得在不同外部磁场环境下,磁力仪的零点偏差如
/nT | /nT | /nT | /pT |
---|---|---|---|
1.062 | -1.081 | -0.009 5 | 10.90 |
2.041 | -2.034 | 0.003 5 | |
-2.064 | 2.093 | 0.014 5 | |
3.027 | -3.046 | -0.009 5 | |
-3.059 | 3.042 | -0.008 5 | |
4.021 | -4.038 | -0.008 5 | |
-4.019 | 3.983 | -0.018 0 |
=10.90 pT
每组实验中测得的与的测量值为一段时间内测量数据的平均值,其标准差分别记为,,视为对应磁场数据的不确定度,如
/nT | /nT | /nT | /nT | /nT | /pT |
---|---|---|---|---|---|
1.062 | 0.003 6 | -1.081 | 0.003 8 | 0.002 6 | 6.90 |
2.041 | 0.002 8 | -2.034 | 0.003 9 | 0.002 4 | |
-2.064 | 0.003 9 | 2.093 | 0.003 6 | 0.002 7 | |
3.027 | 0.004 0 | -3.046 | 0.003 7 | 0.002 7 | |
-3.059 | 0.003 3 | 3.042 | 0.003 6 | 0.002 4 | |
4.021 | 0.004 2 | -4.038 | 0.003 7 | 0.002 8 | |
-4.019 | 0.003 4 | 3.983 | 0.003 9 | 0.002 6 |
对于每组实验下的A类不确定度可通过
(9) |
综合7组实验数据的不确定度,计算得到磁力仪零点测量的A类不确定度为:
=6.90 pT
综上所述,得到磁力仪零点测量不确定度为12.90 pT。
磁屏蔽筒内磁场具有不均匀性,弱场矢量原子磁力仪的敏感单元大小为15 mm25 mm。在磁力仪探头旋转过程中筒内磁场的梯度会导致探头旋转前后-B与B的绝对值不相等,最终影响零点偏差的大小。在上述旋转实验得到的零点偏差中包含了该筒内磁场梯度产生的影响,此影响被计算入磁力仪系统的零点不确定度中,故实际上近零场高精度矢量原子磁力仪的零点不确定度会更好。
3.3 待测矢量原子磁力仪非线性引起的不确定度
利用标定过的人工磁场线圈在±5 nT量程范围内产生多组测量磁场,使用矢量原子磁力仪(VAM)测量该磁场,标定磁力仪系统的线性度,如

图4 矢量原子磁力仪线性度
Fig.4 Linearity of vector atomic magnetometer

图5 线性度偏差
Fig.5 Linearity deviation
3.4 待测矢量原子磁力仪长期稳定性引起的不确定度
将传感器探头放置在无磁转台载物台中央,对磁力仪系统磁场测量的长期稳定度进行测试标定。为尽量保证磁场环境的稳定性,实验选取后半夜进行,记录磁力仪系统在后半夜的6 h连续采样(20 Hz)的磁场测量数据。如

图6 矢量原子磁力仪长期磁场数据
Fig.6 Long⁃term magnetic field data of the vector atomic magnetometer
3.5 综合准确度计算分析
上述标定系统中人工磁场线圈产生校准磁场时,小量程下其合成不确定度为7.55 pT。利用该测试标定装置对实验室研制的近零场矢量原子磁力仪系统测试,得到磁力仪零点不确定度为12.90 pT,6 h内磁力仪长期测量磁场数据的不确定度为5.01 pT。综合线圈产生测量磁场的合成不确定度,小量程内待测矢量原子磁力仪的标准非线性偏差,零点偏差引起的不确定度和长期稳定性引起的测量不确定度,计算得到矢量原子磁力仪在±5 nT范围内测试标定的合成不确定度为15.80 pT,可用于表征磁力仪的准确度指标,如
不确定度来源 | 不确定度 |
---|---|
的合成不确定度 | 7.55 |
标准非线性偏差 | 1.11 |
零点不确定度 | 12.90 |
6 h长期测量不确定度 | 5.01 |
准确度 | 15.80 |
4 结 论
本文提出一种间接测试标定近零场下高精度矢量原子磁力仪的方法,利用磁屏蔽筒实现对地球磁场和外界扰动的屏蔽,通过已完成计量标定的CPT原子磁力仪在大量程下标定一组人工磁场线圈,利用线圈产生小量程范围内的弱磁场进而完成对磁力仪的间接测试标定。
利用本文所述测试标定方法,分析测试标定系统中线圈产生微弱测量磁场时的不确定度来源与大小,测试标定了高精度矢量原子磁力仪在近零场下工作时的零点不确定度、线性度和长期测量不确定度。最终得到矢量原子磁力仪的准确度指标为15.80 pT。
参 考 文 献
BOTO E, SHAH V, HILL R M, et al. Triaxial detection of the neuromagnetic field using optically⁃pumped magnetometry: feasibility and application in children[J]. NeuroImage, 2022, 252: 119027. [百度学术]
FENICI R, BRISINDA D, VENUTI A, et al. Thirty years of clinical magnetocardiography at the catholic university of rome: Diagnostic value and new perspectives for the treatment of cardiac arrhythmias[J]. International Journal of Cardiology, 2013, 168(5): 5113-5115 [百度学术]
WITTEVRONGEL B, HOLMES N, BOTO E, et al. Optically pumped magnetometers for practical MEG⁃based brain⁃computer interfacing[J]. Brain⁃Computer Interface Research: A State⁃of⁃the⁃Art Summary 10, doi: 10.1007/978-3-030-79287-9_4. [百度学术]
BUDKER D. Optical magnetometry[M]. New York: Cambridge University Press, 2013. [百度学术]
SCHULTZ G, MHASKAR R, PROUTY M, et al. Integration of micro⁃fabricated atomic magnetometers on military systems[C]// SPIE. Detection and Sensing of Mines, Explosive Objects, and Obscured Targets XXI.2016: 358-367. [百度学术]
BERTRAND F,JAGER T,BONESS A,et al. A
FRATTER I, LÉGER J M, BERTRAND F, et al. Swarm absolute scalar magnetometers first in‑orbit results[J]. Acta Astronautica, 2016, 121: 76-87. [百度学术]
王赤,李晖,郭孝城,等.太阳系边际探测项目的科学问题[J].深空探测学报(中英文),2020,7(6):517-524. [百度学术]
WANG C, LI H, GUO X C, et al. Scientific Objectives for the Exploration of the Boundary of Solar System[J].Journal of Deep Space Exploration, 2020, 7(6): 517-524. [百度学术]
BEHANNON K W, ACUNA M H, BURLAGA L F, et al. Magnetic field experiment for Voyagers 1 and 2[J]. Space Science Reviews, 1977, 21: 235-257. [百度学术]
SMITH E J, CONNOR B V, FOSTER G T. Measuring the magnetic fields of jupiter and the outer solar system[J]. IEEE Transactions on Magnetics, 1975, 11(4): 962-980. [百度学术]
阎正孝, 李大明. 弱磁场定标装置[J]. 电测与仪表, 1986(12): 8-11. [百度学术]
YAN X Z, LI D M. A weak magnetic field calibration device[J]. Electrical Measurement & Instrumentation, 1986(12): 8-11. [百度学术]
程华富.磁通门磁强计校准装置研究[J].宇航计测技术, 2018, 38(5): 55-59. [百度学术]
CHENG H F. Study on a Calibration Device of Fluxgate Magnetometer[J]. Journal of Astronautic Metrology and Measurement, 2018, 38(5): 55-59. [百度学术]
PAN D, LI J, JIN C, et al. A new calibration method for triaxial fluxgate magnetometer based on magnetic shielding room[J]. IEEE Transactions on Industrial Electronics, 2019, 67(5): 4183-4192. [百度学术]
薛洪波, 王劲东, 宋伟, 等. 具有全向性和绝对准确度的航空CPT原子磁力仪[C]// 中国地球物理学会, 中国地球科学联合学术年会. 中国,北京, 2018: 48-50. [百度学术]
ZHANG X F, CHEN W, ZHU X C, et al. Space advanced technology demonstration satellite[J]. Science China Technological Sciences, 2024, 67(1): 240-258. [百度学术]
XUE H B. Prospect of Vector Atomic Magnetometer in the Exploration of Outer Heliosphere and Nearby Interstellar Medium[J]. ISSI-BJ Magazine, TAIKONG, 2020(20): 38-39. [百度学术]
陈怀艳,曹芸,韩洁.测量不确定度的发展和应用研究[J].宇航计测技术, 2014, 34(5): 65-70. [百度学术]
CHEN H Y, CAO Y, HAN J. Research on Development and Application of Measurement Uncertainty[J]. Journal of Astronautic Metrology and Measurement, 2014, 34(5): 65-70. [百度学术]
彭翔, 郭弘. 光泵原子磁力仪技术[J]. 导航与控制, 2022, 21(5-6): 101-121. [百度学术]
PENG X, GUO H. Techniques in Optically-pumped Atomic Magnetometer[J]. Navigation and Control, 2022, 21(5-6): 101-121. [百度学术]
BOTO E, SHAH V, HILL R M, et al. Triaxial detection of the neuromagnetic field using optically-pumped magnetometry: feasibility and application in children[J]. NeuroImage, 2022, 252: 119027. [百度学术]
钱政,王中宇.误差理论与数据处理[M].北京: 科学出版社.2022: 73-115. [百度学术]