CN114076906B - 一种高温超导全张量磁梯度探头的非正交误差校正方法 - Google Patents
一种高温超导全张量磁梯度探头的非正交误差校正方法 Download PDFInfo
- Publication number
- CN114076906B CN114076906B CN202111353021.8A CN202111353021A CN114076906B CN 114076906 B CN114076906 B CN 114076906B CN 202111353021 A CN202111353021 A CN 202111353021A CN 114076906 B CN114076906 B CN 114076906B
- Authority
- CN
- China
- Prior art keywords
- axis
- squid
- axes
- sensor
- sensitive
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/022—Measuring gradient
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/035—Measuring direction or magnitude of magnetic fields or magnetic flux using superconductive devices
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明为一种高温超导全张量磁梯度探头的非正交误差校正方法,在含有8个SQUID磁强计的高温超导全张量磁梯度探头中,对三个轴向布置的SQUID磁强计传感器之间的非正交误差角度进行测量并对磁梯度分量测量结果进行误差校正。利用三轴亥姆霍兹线圈产生已知参数的磁场与SQUID传感器实际测量值对比的方式从而分别求得8个SQUID磁强计传感器的拾取环敏感轴方向的矢量转换到高温超导全张量磁梯度探头正交直角坐标系的三个坐标轴上所需要知道的角度参数,完成了高温超导全张量磁梯度探头非正交误差的校正,本发明提高了非正交误差的校正精度,由于不需要像传统方法那样在运动过程中进行误差参数的求解从而避免了因为运动产生的磁干扰从而对误差校正结果造成的影响。
Description
技术领域
本发明属于高温超导全张量磁梯度探头领域,具体地而言为一种高温超导全张量磁梯度探头的非正交误差校正方法。
背景技术
超导量子干涉器(SQUID)是一种基于超导约瑟夫森效应和磁通量子化原理的磁通传感器,当SQUID传感器中的约瑟夫森结流入一定的偏置电流时,约瑟夫森结两端的电压变化会呈现出一种宏观量子干涉现象,基于零磁通锁定原理的SQUID读出电路可以将SQUID传感器中的磁通变化转变为线性的电路输出电压变化,研究人员基于此将SQUID传感器应用到生物微弱磁信号探测、仪器无损探伤检测和矿产资源勘探的***开发中。
为了获得磁场的梯度数据,研究人员最开始采用SQUID梯度计作为传感器敏感单元进行***的搭建,但是SQUID梯度计中两个拾取环中间的基线距离由于加工工艺的限制,磁场梯度灵敏度的极值受限,所以为了增加基线距离,研究人员利用如图1所示的SQUID磁强计通过空间上的特殊布局构建成全张量磁梯度探头从而进行磁梯度场全张量9个分量中5个独立分量Bxx、Bxz、Byx、Byz、Bzz的测量。但是该方案存在的问题是由于探头上支撑SQUID传感器的支架本身的加工精度的限制以及SQUID传感器在支架上的安装误差共同导致了三个轴向布置的SQUID磁强计传感器的敏感轴互相之间不垂直,同时同轴向布置的SQUID磁强计传感器的敏感轴互相之间不平行,这两个问题可以统称为高温超导全张量磁梯度测量***中的探头非正交误差。因此,求出非正交误差参数并对磁梯度数据进行校正就显得尤为必要。
发明内容
本发明所要解决的技术问题在于提供一种高温超导全张量磁梯度探头的非正交误差校正方法,解决在含有8个SQUID磁强计的高温超导全张量磁梯度测量***中三个轴向布置的SQUID磁强计传感器的敏感轴互相之间不垂直,同时同轴向布置的SQUID磁强计传感器的敏感轴互相之间不平行的问题,对高温超导全张量磁梯度探头的非正交误差角度进行测量求解并对磁梯度数据进行校正。
本发明是这样实现的,
一种高温超导全张量磁梯度探头的非正交误差校正方法,该方法包括:
在含有8个SQUID磁强计的高温超导全张量磁梯度探头中,对高温超导全张量磁梯度探头坐标系三个轴向布置的SQUID磁强计传感器之间的非正交误差角度进行测量;
三轴亥姆霍兹线圈产生已知参数的磁场与SQUID传感器实际测量值对比,求得8个SQUID磁强计传感器的拾取环敏感轴方向的矢量转换到高温超导全张量磁梯度探头正交直角坐标系的三个坐标轴上所需要的角度参数,完成坐标转化;
采用转换后的温超导全张量磁梯度探头正交直角坐标系的三个坐标轴上所需要的角度参数,对磁梯度全张量的独立分量进行非正交误差的校正。
进一步地,具体包括:
将含有8个SQUID磁强计传感器的高温超导全张量磁梯度探头安装在三轴亥姆霍兹线圈的中心位置,其中在高温超导全张量磁梯度探头坐标系o-xSySzS下,S1号、S2号、S3号SQUID磁强计传感器的拾取环敏感轴法线方向与xS轴平行,S4号、S5号、S6号SQUID磁强计传感器的拾取环敏感轴法线方向与yS轴平行,S7号、S8号SQUID磁强计传感器的拾取环敏感轴法线方向在zS轴上;
利用电流源分别向三轴亥姆霍兹线圈三个轴通入正弦电流,求取8个SQUID传感器的拾取环敏感轴方向与三轴亥姆霍兹线圈三个轴之间的夹角;
建立三轴亥姆霍兹线圈的坐标系为o-xyz,在坐标系o-xyz中,8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈三个轴之间的夹角,得到8个SQUID传感器拾取环敏感轴的法向量坐标,求得S1号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角α17、S2号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角α27、S3号SQUID传感器与S7号SQUID传感器拾取环敏感轴之间的夹角α37、S4号SQUID与S7号SQUID传感器的拾取环敏感轴之间的夹角β47、S5号SQUID与S7号SQUID传感器拾取环敏感轴之间的夹角β57、S6号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角β67、S1号SQUID传感器与S4号SQUID传感器的拾取环敏感轴之间的夹角γ14、S1号SQUID传感器与S5号SQUID传感器拾取环敏感轴之间的夹角γ15、S1号SQUID传感器与S6号SQUID传感器拾取环敏感轴之间的夹角γ16、S2号SQUID传感器与S4号SQUID传感器拾取环敏感轴之间的夹角γ24、S3号SQUID传感器与S4号SQUID传感器拾取环敏感轴之间的夹角γ34;S7号SQUID传感器拾取环敏感轴z7轴与zS轴重合,xS轴在平面x1oz7内,其中x1为S7号SQUID传感器拾取环敏感轴;在高温超导全张量磁梯度探头正交直角坐标系o-xSySzS中,在x1轴和y4轴上分别取A点和B点,其中y4为第4个SQUID传感器拾取环敏感轴,使OA=OB,得出A、B两点的坐标为:
A(OA·sinα17,0,OA·cosα17),B(OA·sinβ47·sinθ,OA·sinβ47·cosθ,OA·cosβ47),
则A、B两点之间的距离为:
在三角形OAB中根据余弦定理得:
联立两式得:
2OA2-2OA2·sinα17·sinβ47·sinθ14-2OA2·cosα17·cosβ47=OA2+OB2-2·OA·OB·cosγ14
整理得:
其中θ14代表S4号SQUID传感器拾取环敏感轴y4轴在平面xSoyS上的投影与yS轴之间的夹角;
同理,在求取θ14过程中,将y4替换为y5和y6,其中y5代表S5号SQUID传感器拾取环敏感轴,y6代表S6号SQUID传感器拾取环敏感轴,求出S5号SQUID传感器拾取环敏感轴y5轴在平面xSoyS上的投影与yS轴之间的夹角θ15和S6号SQUID传感器拾取环敏感轴y6轴在平面xSoyS上的投影与yS轴之间的夹角θ16:
在高温超导全张量磁梯度探头正交直角坐标系o-xSySzS中,x2为第2个SQUID传感器拾取环敏感轴,y4为S4号SQUID传感器拾取环敏感轴,在x2轴和y4轴上分别取C点和D点,使OC=OD,得出C、D两点的坐标为:
C(OC·sinα27·cosθ24,OC·sinα27·sinθ24,OC·cosα27)D(OD·sinβ47·sinθ14,OD·sinβ47·cosθ14,OD·cosβ47)
则C、D两点之间的距离为:
在三角形OCD中根据余弦定理联立上式可得:
又因为sin2θ24+cos2θ24=1
联立两式得:
其中其中θ24代表S2号SQUID传感器拾取环敏感轴x2轴在平面xSoyS上的投影与xS轴之间的夹角;
同理,在求取θ24过程中,将x2替换为x3,其中x3代表S3号SQUID传感器拾取环敏感轴,求出S3号SQUID传感器拾取环敏感轴x3轴在平面xSoyS上的投影与xS轴之间的夹角θ34:
其中
在高温超导全张量磁梯度探头正交直角坐标系o-xSySzS中,z7代表第7个SQUID传感器拾取环敏感轴,z8代表S8号SQUID传感器拾取环敏感轴,z7轴与z8轴之间的夹角
进一步地,对磁梯度全张量的独立分量进行非正交误差的校正包括:8个SQUID传感器输出分别为B1、B2、B3、B4、B5、B6、B7、B8,则经过非正交误差校正后的磁梯度全张量的5个独立分量为:
其中到/>分别指第1个到第8个SQUID传感器输出的零点偏移误差。
进一步地,求取8个SQUID传感器的拾取环敏感轴方向与三轴亥姆霍兹线圈三个轴之间的夹角包括:
利用电流源向三轴亥姆霍兹线圈的X轴通入正弦电流从而使三轴亥姆霍兹线圈的X轴在线圈中心的磁场均匀区内产生BX=AXsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为代表传感器编号;
分别建立8个目标函数,代表传感器编号,利用最小二乘法分别求出8个SQUID传感器输出的正弦信号的峰值/>和8个SQUID传感器的零点偏移误差/>
通过和AX分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈X轴之间的夹角/>
停止向X轴线圈中通入电流,利用电流源向三轴亥姆霍兹线圈的Y轴通入正弦电流从而使三轴亥姆霍兹线圈的Y轴在线圈中心的磁场均匀区内产生BY=AYsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为
分别建立8个目标函数,利用最小二乘法分别求出8个SQUID传感器输出的正弦信号的峰值/>
通过和AY分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈Y轴之间的夹角/>
停止向Y轴线圈中通入电流,利用电流源向三轴亥姆霍兹线圈的Z轴通入正弦电流从而使三轴亥姆霍兹线圈的Z轴在线圈中心的磁场均匀区内产生BZ=AZsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为
分别建立8个目标函数,利用最小二乘法分别求出8个SQUID传感器输出的正弦信号的峰值/>通过/>和AZ分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈Z轴之间的夹角/>
本发明与现有技术相比,有益效果在于:
利用三轴亥姆霍兹线圈产生已知参数的磁场与SQUID传感器实际测量值对比的方式从而分别求得了8个SQUID磁强计传感器的拾取环敏感轴方向的矢量转换到高温超导全张量磁梯度探头正交直角坐标系o-xSySzS的三个坐标轴上所需要知道的角度参数,从而完成了高温超导全张量磁梯度探头非正交误差的校正,该方法相对于传统方法的好处是三轴亥姆霍兹线圈产生的磁场相对于国际地磁场模型求得的磁场更加准确,因此提高了非正交误差的校正精度,另一方面,由于不需要像传统方法那样在运动过程中进行误差参数的求解从而避免了因为运动产生的磁干扰从而对误差校正结果造成的影响。
附图说明
图1为高温超导全张量磁梯度探头中8个SQUID磁强计传感器的安装位置示意图;
图2为高温超导全张量磁梯度探头固定在三轴亥姆霍兹线圈中的磁场均匀区示意图;
图3为SQUID传感器敏感轴所在轴向在三轴亥姆霍兹线圈正交直角坐标系o-xyz中与x、y、z轴夹角示意图;
图4为根据x1、y4、z7轴建立高温超导全张量磁梯度探头正交直角坐标系o-xSySzS并求取y4轴在xSoyS平面上的投影与yS轴的夹角θ14示意图;
图5为求取y5轴在xSoyS平面上的投影与yS轴的夹角θ15示意图;
图6为求取y6轴在xSoyS平面上的投影与yS轴的夹角θ16示意图;
图7为求取x2轴在xSoyS平面上的投影与xS轴的夹角θ24示意图;
图8为求取x3轴在xSoyS平面上的投影与xS轴的夹角θ34示意图;
图9为求取z8轴与z7轴的夹角ε78示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
一种高温超导全张量磁梯度探头的非正交误差校正方法,该方法包括:
在含有8个SQUID磁强计的高温超导全张量磁梯度探头中,对高温超导全张量磁梯度探头坐标系三个轴向布置的SQUID磁强计传感器之间的非正交误差角度进行测量;
三轴亥姆霍兹线圈产生已知参数的磁场与SQUID传感器实际测量值对比,求得8个SQUID磁强计传感器的拾取环敏感轴方向的矢量转换到高温超导全张量磁梯度探头正交直角坐标系的三个坐标轴上所需要的角度参数,完成坐标转化;
采用转换后的温超导全张量磁梯度探头正交直角坐标系的三个坐标轴上所需要的角度参数,对磁梯度全张量的独立分量进行非正交误差的校正。
如图1所示为高温超导全张量磁梯度探头中8个SQUID磁强计传感器的安装位置示意图,SQUID传感器的编号为Si(i=1,2……8),箭头方向为SQUID传感器拾取环敏感轴方向,如图2所示将高温超导全张量磁梯度测量***调试到工作状态并将探头安装在三轴亥姆霍兹线圈的中心位置;
利用电流源向三轴亥姆霍兹线圈的X轴通入正弦电流从而使三轴亥姆霍兹线圈的X轴在线圈中心的磁场均匀区内产生BX=AXsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为(i=1,2……8代表传感器编号);
分别建立8个目标函数,(i=1,2……8代表传感器编号),利用最小二乘法可以分别求出8个SQUID传感器输出的正弦信号的峰值/>和8个SQUID传感器的零点偏移误差/>
通过和AX可以分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈X轴之间的夹角/>(i=1,2……8代表传感器编号);
停止向X轴线圈中通入电流,利用电流源向三轴亥姆霍兹线圈的Y轴通入正弦电流从而使三轴亥姆霍兹线圈的Y轴在线圈中心的磁场均匀区内产生BY=AYsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为(i=1,2……8代表传感器编号);
分别建立8个目标函数,(i=1,2……8代表传感器编号),利用最小二乘法可以分别求出8个SQUID传感器输出的正弦信号的峰值/>
通过和AY可以分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈Y轴之间的夹角/>(i=1,2……8代表传感器编号);
停止向Y轴线圈中通入电流,利用电流源向三轴亥姆霍兹线圈的Z轴通入正弦电流从而使三轴亥姆霍兹线圈的Z轴在线圈中心的磁场均匀区内产生BZ=AZsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为(i=1,2……8代表传感器编号);
分别建立8个目标函数,(i=1,2……8代表传感器编号),利用最小二乘法可以分别求出8个SQUID传感器输出的正弦信号的峰值/>通过/>和AZ可以分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈Z轴之间的夹角(i=1,2……8代表传感器编号);
分别定义S1、S2、S3、S4、S5、S6、S7和S8传感器的敏感轴为x1轴、x2轴、x3轴、y4轴、y5轴、y6轴、z7轴、z8轴;如图3所示,三轴亥姆霍兹线圈的坐标系为o-xyz,在坐标系o-xyz中,根据x1轴、x2轴、x3轴、y4轴、y5轴、y6轴、z7轴、z8轴与x轴、y轴、z轴之间的夹角αi、βi、γi(i=1,2……8代表传感器编号)可以分别得到x1轴、x2轴、x3轴、y4轴、y5轴、y6轴、z7轴、z8轴的法向量坐标,根据定理:空间直角坐标系中两个向量a=(x1,y1,z1)和b=(x2,y2,z2)之间的夹角可以分别求得如图4、5、6、7、8中所示的夹角α17、α27、α37、β47、β57、β67、γ14、γ15、γ16、γ24、γ34;S1号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角α17、S2号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角α27、S3号SQUID传感器与S7号SQUID传感器拾取环敏感轴之间的夹角α37、S4号SQUID与S7号SQUID传感器的拾取环敏感轴之间的夹角β47、S5号SQUID与S7号SQUID传感器拾取环敏感轴之间的夹角β57、S6号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角β67、S1号SQUID传感器与S4号SQUID传感器的拾取环敏感轴之间的夹角γ14、S1号SQUID传感器与S5号SQUID传感器拾取环敏感轴之间的夹角γ15、S1号SQUID传感器与S6号SQUID传感器拾取环敏感轴之间的夹角γ16、S2号SQUID传感器与S4号SQUID传感器拾取环敏感轴之间的夹角γ24、S3号SQUID传感器与S4号SQUID传感器拾取环敏感轴之间的夹角γ34;
建立如图4所示的高温超导全张量磁梯度探头正交直角坐标系o-xSySzS,其中z7轴与zS轴重合,xS轴在平面x1oz7内;在高温超导全张量磁梯度探头正交直角坐标系o-xSySzS中,在x1轴和y4轴上分别取A点和B点,使OA=OB,则可以得出A、B两点的坐标为:
A(OA·sinα17,0,OA·cosα17),B(OA·sinβ47·sinθ,OA·sinβ47·cosθ,OA·cosβ47),
则A、B两点之间的距离为:
在三角形OAB中根据余弦定理可得:
联立两式可得:
2OA2-2OA2·sinα17·sinβ47·sinθ14-2OA2·cosα17·cosβ47=OA2+OB2-2·OA·OB·cosγ14
整理得:
同理可得:
在图7的高温超导全张量磁梯度探头正交直角坐标系o-xSySzS中,在x2轴和y4轴上分别取C点和D点,使OC=OD,则可以得出C、D两点的坐标为:
C(OC·sinα27·cosθ24,OC·sinα27·sinθ24,OC·cosα27)D(OD·sinβ47·sinθ14,OD·sinβ47·cosθ14,OD·cosβ47)
则C、D两点之间的距离为:
在三角形OCD中根据余弦定理联立上式可得:
又因为sin2θ24+cos2θ24=1
联立两式可得:
其中
同理可得:
其中
如图9所示z7轴与z8轴之间的夹角
设传感器S1、S2、S3、S4、S5、S6、S7、S8的输出分别为B1、B2、B3、B4、B5、B6、B7、B8,则经过非正交误差校正后的磁梯度全张量的5个独立分量为:
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种高温超导全张量磁梯度探头的非正交误差校正方法,其特征在于,该方法包括:
在含有8个SQUID磁强计的高温超导全张量磁梯度探头中,对高温超导全张量磁梯度探头坐标系三个轴向布置的SQUID磁强计传感器之间的非正交误差角度进行测量;
三轴亥姆霍兹线圈产生已知参数的磁场与SQUID传感器实际测量值对比,求得8个SQUID磁强计传感器的拾取环敏感轴方向的矢量转换到高温超导全张量磁梯度探头正交直角坐标系的三个坐标轴上所需要的角度参数,完成坐标转化;
采用转换后的温超导全张量磁梯度探头正交直角坐标系的三个坐标轴上所需要的角度参数,对磁梯度全张量的独立分量进行非正交误差的校正;
具体包括:
将含有8个SQUID磁强计传感器的高温超导全张量磁梯度探头安装在三轴亥姆霍兹线圈的中心位置,其中在高温超导全张量磁梯度探头坐标系o-xSySzS下,S1号、S2号、S3号SQUID磁强计传感器的拾取环敏感轴法线方向与xS轴平行,S4号、S5号、S6号SQUID磁强计传感器的拾取环敏感轴法线方向与yS轴平行,S7号、S8号SQUID磁强计传感器的拾取环敏感轴法线方向在zS轴上;
利用电流源分别向三轴亥姆霍兹线圈三个轴通入正弦电流,求取8个SQUID传感器的拾取环敏感轴方向与三轴亥姆霍兹线圈三个轴之间的夹角;
建立三轴亥姆霍兹线圈的坐标系为o-xyz,在坐标系o-xyz中,8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈三个轴之间的夹角,得到8个SQUID传感器拾取环敏感轴的法向量坐标,求得S1号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角α17、S2号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角α27、S3号SQUID传感器与S7号SQUID传感器拾取环敏感轴之间的夹角α37、S4号SQUID与S7号SQUID传感器的拾取环敏感轴之间的夹角β47、S5号SQUID与S7号SQUID传感器拾取环敏感轴之间的夹角β57、S6号SQUID传感器与S7号SQUID传感器的拾取环敏感轴之间的夹角β67、S1号SQUID传感器与S4号SQUID传感器的拾取环敏感轴之间的夹角γ14、S1号SQUID传感器与S5号SQUID传感器拾取环敏感轴之间的夹角γ15、S1号SQUID传感器与S6号SQUID传感器拾取环敏感轴之间的夹角γ16、S2号SQUID传感器与S4号SQUID传感器拾取环敏感轴之间的夹角γ24、S3号SQUID传感器与S4号SQUID传感器拾取环敏感轴之间的夹角γ34;
S7号SQUID传感器拾取环敏感轴z7轴与zS轴重合,xS轴在平面x1oz7内,其中x1为S7号SQUID传感器拾取环敏感轴;在高温超导全张量磁梯度探头正交直角坐标系o-xSySzS中,在x1轴和y4轴上分别取A点和B点,其中y4为S4号SQUID传感器拾取环敏感轴,使OA=OB,得出A、B两点的坐标为:
A(OA·sinα17,0,OA·cosα17),B(OA·sinβ47·sinθ,OA·sinβ47·cosθ,OA·cosβ47),
则A、B两点之间的距离为:
在三角形OAB中根据余弦定理得:
联立两式得:
2OA2-2OA2·sinα17·sinβ47·sinθ14-2OA2·cosα17·cosβ47=OA2+OB2-2·OA·OB·cosγ14
整理得:
其中θ14代表S4号SQUID传感器拾取环敏感轴y4轴在平面xSoyS上的投影与yS轴之间的夹角;
同理,在求取θ14过程中,将y4替换为y5和y6,其中y5代表S5号SQUID传感器拾取环敏感轴,y6代表S6号SQUID传感器拾取环敏感轴,求出S5号SQUID传感器拾取环敏感轴y5轴在平面xSoyS上的投影与yS轴之间的夹角θ15和S6号SQUID传感器拾取环敏感轴y6轴在平面xSoyS上的投影与yS轴之间的夹角θ16:
在高温超导全张量磁梯度探头正交直角坐标系o-xSySzS中,x2为S2号SQUID传感器拾取环敏感轴,y4为S4号SQUID传感器拾取环敏感轴,在x2轴和y4轴上分别取C点和D点,使OC=OD,得出C、D两点的坐标为:
C(OC·sinα27·cosθ24,OC·sinα27·sinθ24,OC·cosα27)D(OD·sinβ47·sinθ14,OD·sinβ47·cosθ14,OD·cosβ47)
则C、D两点之间的距离为:
在三角形OCD中根据余弦定理联立上式可得:
又因为sin2θ24+cos2θ24=1
联立两式得:
其中其中θ24代表S2号SQUID传感器拾取环敏感轴x2轴在平面xSoyS上的投影与xS轴之间的夹角;
同理,在求取θ24过程中,将x2替换为x3,其中x3代表S3号SQUID传感器拾取环敏感轴,求出S3号SQUID传感器拾取环敏感轴x3轴在平面xSoyS上的投影与xS轴之间的夹角θ34:
其中
在高温超导全张量磁梯度探头正交直角坐标系o-xSySzS中,z7代表S7号SQUID传感器拾取环敏感轴,z8代表S8号SQUID传感器拾取环敏感轴,z7轴与z8轴之间的夹角
2.按照权利要求1所述的一种高温超导全张量磁梯度探头的非正交误差校正方法,其特征在于,对磁梯度全张量的独立分量进行非正交误差的校正包括:8个SQUID传感器输出分别为B1、B2、B3、B4、B5、B6、B7、B8,则经过非正交误差校正后的磁梯度全张量的5个独立分量为:
其中到/>分别指第1个到第8个SQUID传感器输出的零点偏移误差。
3.按照权利要求1所述的一种高温超导全张量磁梯度探头的非正交误差校正方法,其特征在于,求取8个SQUID传感器的拾取环敏感轴方向与三轴亥姆霍兹线圈三个轴之间的夹角包括:
利用电流源向三轴亥姆霍兹线圈的X轴通入正弦电流从而使三轴亥姆霍兹线圈的X轴在线圈中心的磁场均匀区内产生BX=AXsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为代表传感器编号;
分别建立8个目标函数,代表传感器编号,利用最小二乘法分别求出8个SQUID传感器输出的正弦信号的峰值/>和8个SQUID传感器的零点偏移误差/>
通过和AX分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈X轴之间的夹角/>
停止向X轴线圈中通入电流,利用电流源向三轴亥姆霍兹线圈的Y轴通入正弦电流从而使三轴亥姆霍兹线圈的Y轴在线圈中心的磁场均匀区内产生BY=AYsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为
分别建立8个目标函数,利用最小二乘法分别求出8个SQUID传感器输出的正弦信号的峰值/>
通过和AY分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈Y轴之间的夹角/>
停止向Y轴线圈中通入电流,利用电流源向三轴亥姆霍兹线圈的Z轴通入正弦电流从而使三轴亥姆霍兹线圈的Z轴在线圈中心的磁场均匀区内产生BZ=AZsinx随时间变化的磁场,与电流源提供的电流信号同步采集8个SQUID传感器的输出为
分别建立8个目标函数,利用最小二乘法分别求出8个SQUID传感器输出的正弦信号的峰值/>通过/>和AZ分别求出8个SQUID传感器拾取环敏感轴方向与三轴亥姆霍兹线圈Z轴之间的夹角/>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111353021.8A CN114076906B (zh) | 2021-11-16 | 2021-11-16 | 一种高温超导全张量磁梯度探头的非正交误差校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111353021.8A CN114076906B (zh) | 2021-11-16 | 2021-11-16 | 一种高温超导全张量磁梯度探头的非正交误差校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114076906A CN114076906A (zh) | 2022-02-22 |
CN114076906B true CN114076906B (zh) | 2023-10-17 |
Family
ID=80283802
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111353021.8A Active CN114076906B (zh) | 2021-11-16 | 2021-11-16 | 一种高温超导全张量磁梯度探头的非正交误差校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114076906B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20030014561A (ko) * | 2001-08-11 | 2003-02-19 | 한국표준과학연구원 | 2 주파수 레이저 간섭계에서의 길이 측정을 위한 위상각측정방법 및 비선형 오차 보상방법 그리고 이를 이용한길이 측정을 위한 위상각 측정방법 및 측정 시스템 |
JP2004028955A (ja) * | 2002-06-28 | 2004-01-29 | Sumitomo Denko Hightecs Kk | 食製品検査装置 |
CN103885020A (zh) * | 2014-03-28 | 2014-06-25 | 卢兆兴 | 一种基于自适应遗传算法的三轴磁强计误差校正技术 |
CN104345348A (zh) * | 2014-11-07 | 2015-02-11 | 吉林大学 | 航空超导全张量磁梯度测量***相关参数获取装置及方法 |
CN104569884A (zh) * | 2013-10-18 | 2015-04-29 | 中国科学院上海微***与信息技术研究所 | 一种超导量子干涉器件三轴磁强计的标定装置及方法 |
WO2016031674A1 (ja) * | 2014-08-26 | 2016-03-03 | 旭化成エレクトロニクス株式会社 | 誤差補正装置、回転角検出装置、回転角センサ、誤差補正方法、およびプログラム |
CN106997035A (zh) * | 2017-03-29 | 2017-08-01 | 吉林大学 | 一种基于磁梯度不变量的磁梯度计校正方法 |
CN107229021A (zh) * | 2017-05-22 | 2017-10-03 | 中国科学院上海微***与信息技术研究所 | 三维磁场测量组件及制备方法 |
CN108387952A (zh) * | 2018-02-06 | 2018-08-10 | 吉林大学 | 三分量squid传感器正交度校正装置及校正方法 |
CN112346149A (zh) * | 2020-10-10 | 2021-02-09 | 成都理工大学 | 一种三分量井中磁力仪测量值校正*** |
CN112698258A (zh) * | 2021-01-20 | 2021-04-23 | 中国人民解放军海军工程大学 | 一种三轴磁强计的一体化误差校正方法 |
CN113156355A (zh) * | 2021-03-31 | 2021-07-23 | 吉林大学 | 一种超导全张量磁梯度测量装置的磁干扰补偿方法 |
CN113156354A (zh) * | 2021-03-29 | 2021-07-23 | 吉林大学 | 一种无屏蔽环境下的三轴亥姆霍兹线圈非正交误差测量方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10288698B2 (en) * | 2016-06-13 | 2019-05-14 | Allegro Microsystems, Llc | Magnetic field sensor having alignment error correction |
-
2021
- 2021-11-16 CN CN202111353021.8A patent/CN114076906B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20030014561A (ko) * | 2001-08-11 | 2003-02-19 | 한국표준과학연구원 | 2 주파수 레이저 간섭계에서의 길이 측정을 위한 위상각측정방법 및 비선형 오차 보상방법 그리고 이를 이용한길이 측정을 위한 위상각 측정방법 및 측정 시스템 |
JP2004028955A (ja) * | 2002-06-28 | 2004-01-29 | Sumitomo Denko Hightecs Kk | 食製品検査装置 |
CN104569884A (zh) * | 2013-10-18 | 2015-04-29 | 中国科学院上海微***与信息技术研究所 | 一种超导量子干涉器件三轴磁强计的标定装置及方法 |
CN103885020A (zh) * | 2014-03-28 | 2014-06-25 | 卢兆兴 | 一种基于自适应遗传算法的三轴磁强计误差校正技术 |
WO2016031674A1 (ja) * | 2014-08-26 | 2016-03-03 | 旭化成エレクトロニクス株式会社 | 誤差補正装置、回転角検出装置、回転角センサ、誤差補正方法、およびプログラム |
CN104345348A (zh) * | 2014-11-07 | 2015-02-11 | 吉林大学 | 航空超导全张量磁梯度测量***相关参数获取装置及方法 |
CN106997035A (zh) * | 2017-03-29 | 2017-08-01 | 吉林大学 | 一种基于磁梯度不变量的磁梯度计校正方法 |
CN107229021A (zh) * | 2017-05-22 | 2017-10-03 | 中国科学院上海微***与信息技术研究所 | 三维磁场测量组件及制备方法 |
CN108387952A (zh) * | 2018-02-06 | 2018-08-10 | 吉林大学 | 三分量squid传感器正交度校正装置及校正方法 |
CN112346149A (zh) * | 2020-10-10 | 2021-02-09 | 成都理工大学 | 一种三分量井中磁力仪测量值校正*** |
CN112698258A (zh) * | 2021-01-20 | 2021-04-23 | 中国人民解放军海军工程大学 | 一种三轴磁强计的一体化误差校正方法 |
CN113156354A (zh) * | 2021-03-29 | 2021-07-23 | 吉林大学 | 一种无屏蔽环境下的三轴亥姆霍兹线圈非正交误差测量方法 |
CN113156355A (zh) * | 2021-03-31 | 2021-07-23 | 吉林大学 | 一种超导全张量磁梯度测量装置的磁干扰补偿方法 |
Non-Patent Citations (13)
Title |
---|
Adaptive sampling of potential-field data: a direct approach to compressive inversion;Foks L等;75th EAGE Conference and Exhibition incorporating;第79卷(第1期);1-9 * |
Axial high-temperature RF SQUID gradiometer systyem for geomagnetic prospecting;Jing Zhao等;2009 9th International Conference on Electronic Measurement & Instruments;第1卷(第1期);704-707 * |
Development and evaluation of an airborne superconducting;T. Jeffrey Gamey;Strategic environmental Reasearch and Development Program;第1卷(第1期);1-44 * |
Multiple cross-gradient joint inversion for geospectral imaging;Gopal B.V 等;Reviews of Geophysics;第34卷(第19期);2473-2488 * |
一种电子磁罗盘航向误差的自适应补偿方法;冯毅博等;仪器仪表学报;第35卷(第11期);2607-2614 * |
三分量磁通门传感器的三轴正交校正与测量;王一等;山东农业大学学报;第46卷(第2期);232-237 * |
五棱台是全张量磁梯度探头侧面倾角优化方法;申茂东等;吉林大学学报;第46卷(第5期);1732-1738 * |
地球物理矢量场磁测技术的研究进展;林君等;中国学科;第62卷(第23期);2606-2618 * |
基于线性误差模型的磁张量***校正;张光等;吉林大学学报(工学版);第45卷(第3期);1012-1016 * |
磁梯度张量***发展及其误差校正研究现状;李青竹等;装甲兵工程学院学报;第31卷(第6期);72-81 * |
航磁三分量探测***研究;王一等;吉林大学;第1卷(第1期);1-112 * |
量子磁力仪再评说;董浩斌等;工程地球物理学报;第7卷(第4期);460-470 * |
高温超导磁梯度仪测控***的设计与实现;赵静等;计算机测量与控制;第19卷(第5期);1052-1075 * |
Also Published As
Publication number | Publication date |
---|---|
CN114076906A (zh) | 2022-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102853760B (zh) | 一种三轴磁传感器磁轴垂直度的标定方法 | |
WO2018006020A1 (en) | Calibration and monitoring for 3-axis magnetometer arrays of arbitrary geometry | |
CN110118948B (zh) | 一种基于超导量子干涉仪的总场测量方法及装置 | |
JP6525336B2 (ja) | 3軸デジタルコンパス | |
CN113325353B (zh) | 一种磁强计空间姿态标定方法及*** | |
CN111413750B (zh) | 一种用于电缆定位的磁通门传感器的误差校正方法 | |
CN107121655B (zh) | 非屏蔽serf原子磁力仪磁场抵消线圈非正交角测量装置及测量方法 | |
JP6664568B1 (ja) | 磁場計測装置、磁場計測方法、磁場計測プログラム | |
CN112130217B (zh) | 线圈矢量磁力仪几何轴与磁轴夹角电学检测***及方法 | |
CN109633491B (zh) | 全张量磁梯度测量***安装误差的标定装置及标定方法 | |
JP2015075465A (ja) | 3次元磁界測定装置及び3次元磁界測定方法 | |
CN108508390A (zh) | 矢量磁传感器的标定方法及其简便标定方法 | |
JP2019215322A (ja) | 磁場計測装置、磁場計測方法、磁場計測プログラム | |
CN107490775B (zh) | 一种三轴线圈常数与非正交角一体化测量方法 | |
JP2004151056A (ja) | 弱磁界発生装置および磁界センサの検査方法 | |
CN114076906B (zh) | 一种高温超导全张量磁梯度探头的非正交误差校正方法 | |
JP4151785B2 (ja) | 方位角計測装置及び方位角計測プログラム、並びに方位角計測方法 | |
Yue et al. | Error Calibration for Full Tensor Magnetic Gradiometer Probe Based on Coordinate Transformation Method | |
CN114779144B (zh) | 一种测量三轴磁强计安装矩阵的方法、芯片和装置 | |
CN111624531B (zh) | 一种用于tmr三轴集成磁传感器的分量解算方法 | |
Zhang et al. | A novel auto-calibration method of the vector magnetometer | |
CN107966670B (zh) | 一种超导全张量探测装置及超导全张量探测方法 | |
WO2020145344A1 (ja) | 磁場校正装置及びこれを用いた磁気計測装置のキャリブレーション方法 | |
CN110702102B (zh) | 一种用于通航飞机的磁导航***及其导航方法 | |
CN110568387B (zh) | 一种基于磁梯度张量的航天器磁矩测试方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |