CN105891755B - 航空器挂载式磁通门磁梯度张量仪的校正方法 - Google Patents
航空器挂载式磁通门磁梯度张量仪的校正方法 Download PDFInfo
- Publication number
- CN105891755B CN105891755B CN201610102276.XA CN201610102276A CN105891755B CN 105891755 B CN105891755 B CN 105891755B CN 201610102276 A CN201610102276 A CN 201610102276A CN 105891755 B CN105891755 B CN 105891755B
- Authority
- CN
- China
- Prior art keywords
- tensor
- fluxgate
- error
- component
- axis
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
- G01R35/005—Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明涉及一种航空器挂载式磁通门磁梯度张量仪的校正方法,一建立单个磁通门误差模型,二建立磁梯度张量分量误差校正模型,三在高空中磁场均匀的区域采集校正数据,四将上述张量分量的误差校正,五在地面磁场不均匀的区域采集校正数据,六求解张量分量的标度因子误差,七结合高空校正和地面校正所获得校正系数,八完整计算校正后的磁梯度张量。本发明的校正算法与现有的航空磁场梯度张量校正算法相比,不仅仅对单个张量分量进行校正,还包含了张量作为一个整体的校正,从而使得张量分量和整体的校正结果更加准确,同时其应用条件更符合实际的航空磁梯度张量探测区域的地质条件。
Description
技术领域:
本发明涉及一种航空物探磁测仪器校正,尤其是航空器挂载式磁通门磁梯度张量仪的校正方法。
背景技术:
磁梯度张量仪测量磁场三分量沿正交坐标系的空间变化率。磁梯度张量仪受地磁影响小,获取信息丰富,特别适合发现浅层潜伏矿,磁性运动目标等,是磁测仪器发展方向之一。在航空磁梯度张量测量中,为了避免张量仪受到航空器如直升机等的磁干扰,通常是把张量仪挂载到离航空器较远的地方,然后在飞行过程中测量目标的磁场梯度,本发明正是针对此类应用时张量仪的校正方法。
由于需要对磁场矢量敏感的传感器,故目前张量仪的实现主要有两种方式:一种是基于超导量子干涉器件的张量仪,如德国的航空超导磁梯度张量仪[R.Stolz et al,Magnetic full-tensor SQUID gradiometer system for geophysical applications,The Leading Edge,2006,25(2):178-180.],澳大利亚的航空超导磁梯度张量仪[P.Schmidt et al,GETMAG-a SQUID magnetic tensor gradiometer for mineral andoil exploration,Exploration Geophysics,2004,35:297-305.]等。超导张量仪虽然具有较高的分辨率,但是价格昂贵,且不能在常温下使用,应用领域有较大的局限性。另一种是基于磁通门传感器构造的张量仪,如美国地质调查局的正四面体形磁通门张量仪[PJ.Brown et al,A case study of magnetic gradient tensor invariants appliedto the UXO problem,U.S.Geological Survey,2004:1-4.]。美国IBM公司研制的磁通门张量仪[RH.Koch et al,Room temperature three sensor magnetic field gradiometer,Review of Scientific Instruments,1996,64(1):230-235.]。
虽然磁通门构造的张量仪具有成本低,温度范围广,分辨率较高等优点,但是磁通门存在三个感应轴标度因子不完全一致造成的误差、三个感应轴不严格正交造成的误差、零点存在偏移、动态误差等问题。这些误差以及各个磁通门误差特性不一致造成磁通门张量仪误差大的问题,从而严重影响了航空磁梯度张量仪的应用效果。针对此问题,目前存在两种校正方法。第一种为实验室标定法,一般是在磁屏蔽环境中,利用亥姆赫兹线圈或麦克斯韦线圈,以及标准电流源等装置进行标定。这类方法除了标定环境难以建立以外,最大问题是标定所处的环境(如屏蔽环境中磁场很小,测量装置处于静止状态)和实际的航空磁梯度张量测量环境(如存在地球磁场,测量装置处于运动状态)相差甚远,直接应用实验室内获得的校正参数并不能保证航空测量应用中校正的正确性。第二种方法是在地面上找一块磁场均匀的区域,然后以多个磁通门测量的总场值或磁场分量值相等为条件求解线性或非线性方程以获得校正参数。这类方法有三个主要的问题:第一,在磁场均匀区域内,单个梯度分量为零是磁梯度张量仪校正中的必要条件,而非充要条件。既不能保证对非零梯度情况响应的正确性,也不能保证各个梯度分量形成张量后作为整体的正确性。第二,在地面难以复现飞行过程中磁通门的动态特性,故难以辨识出多个磁通门动态特性不一致性的问题。而航空磁梯度张量测量始终处于运动状态,动态特性的不一致同样也是形成运动噪声,淹没目标磁探测信号的原因之一。第三,实际航空磁场测量的目标区域中或附近,地面磁性异常往往较大,难以找到磁场均匀的区域进行校正。
航空挂载式磁通门磁梯度张量仪中磁通门的测量误差来源于以下几项:磁通门传感器三个感应轴标度因子不完全一致造成的误差、三个感应轴不严格正交造成的误差、零点存在偏移、由机械安装误差造成的非对准误差、动态特性不一致造成的误差。
发明内容:
本发明的目的是针对上述现有技术的不足,提供一种提高其测量准确性的航空器挂载式磁通门磁梯度张量仪的校正方法。
本发明的目的是通过以下技术方案实现的:
航空器挂载式磁通门磁梯度张量仪的校正方法,包括以下步骤:
步骤一、建立单个磁通门误差模型:
把惯性导航系统和磁通门张量仪进行刚体连接,把各个磁通门传感器坐标系x,y,z轴分别向惯性导航坐标系的x,y,z轴方向对齐,则单个磁通门误差模型为:
式中:F为标度因子误差矩阵;E为非正交误差矩阵;R为零位误差矢量;M为磁通门坐标系与惯性导航坐标系非对准误差矩阵;Dj为磁通门测量值对时间的j阶导数的描述,即动态特性描述矩阵;T为磁通门在自身测量坐标系下的实测值;B为磁通门在惯导坐标系下的真实磁场值;Φ为磁通门的测量噪声;
每个误差矩阵或矢量的具体表达式如下:
(1)标度因子误差矩阵:
其中Sx,Sy,Sz表示磁通门感应轴的标度因子误差,
(2)非正交误差矩阵:
假设实际磁通门坐标系x轴与理想磁通门坐标系x轴的方向一致,ρ为实际磁通门坐标系y轴与理想磁通门坐标系y轴的夹角;为实际磁通门坐标系z轴与理想磁通门坐标系z轴的夹角在磁通门xoz平面上的投影;λ为实际磁通门坐标系z轴与理想磁通门坐标系z轴的夹角在磁通门yoz平面上的投影,
(3)非对准误差矩阵:
式中cα=cosα,cβ=cosβ,cγ=cosγ,sα=sinα,sβ=sinβ,sγ=sinγ,α、β和γ分别为实际磁通门坐标系相对惯导坐标系,绕实际磁通门坐标系x轴,y轴和z轴的非对准角。当非对准角非常小的时候可以将非对准误差矩阵近似表示为:
(4)零位误差矢量:
式中μx,μy,μz分别为磁通门三个轴分别存在的零位误差,
(5)动态特性描述矩阵:
当磁通门由一阶系统描述时取j=1,D1为一阶导数前面的系数矩阵。
将所述三轴磁通门的误差模型改写为校正模型:
式中:C0为T第0阶导数的校正参数矩阵;Cj为Τ第j阶导数的校正参数矩阵;O为零位误差校正矢量;bx,by,bz分别表示磁通门在惯性导航坐标系x,y,z轴上的磁场值,则磁通门中单个磁场分量误差校正模型可表示为:
式中:u表示x,y,z中一个轴,Pu0为矩阵C0中对应bu的行向量校正参数;Quj为矩阵Cj中对应bu的行向量校正参数;Ou为总偏差矢量On中对应bu的校正参数;
步骤二、建立磁梯度张量分量误差校正模型:
磁场沿x,y,z三个方向的二阶导数可构成磁梯度张量G,尽管G有九个分量,由于在不包括场源的域内,磁场的散度和旋度为0,即对称,且迹为0,故只有五个分量是独立的,可表示为:
磁通门张量仪由不同方式构建,主要有十字、四方体和三角形等结构,但都是利用差分近似微分的方式来表示磁梯度张量G各个分量。以十字型结构为例,在惯性导航坐标系下,张量分量表示方式为:
式中:m,n表示磁通门标号,取1,2,3,4;u,v都代表坐标轴的方向,取x,y,z;guv为惯性导航坐标系下某个张量分量;ΔBu为惯性导航坐标系下沿u轴方向的磁场变化量;Δsv表示磁通门n与磁通门m之间的基线距离。
将所述张量分量误差模型改写为:
式中:Ouv为张量分量guv的总偏置误差;
步骤三、在高空中磁场均匀的区域采集校正数据:
将磁通门磁梯度张量仪挂载到航空器如直升机下方较远处。航空器飞行到高空,然后以方形,八字形等多种轨迹进行正常飞行,张量仪采集不同姿态下带误差的磁场测量数据;
步骤四、因为高空中地球磁场的梯度近似为0,设张量分量guv误差模型中磁通门n中沿u轴方向测量值前面的校正参数为puv,将上述张量分量的误差校正模型改写为:
步骤五、在地面磁场不均匀的区域采集校正数据:
在地面上有磁梯度的区域,把张量仪的传感器部分绕其中心点进行旋转,记录此过程的测量值;
步骤六、以不变量进行约束,求解张量分量的标度因子误差:
将高空校正之后所获得的各项校正系数重新带入到公式(13)中,可以获得各个张量分量的相对量值,表达为:
由于式(1)中误差矩阵的数值并不相同,因此每个分量所对应的puv也各不相同,故此时的张量分量可以认为存在张量分量标度因子误差。张量的不变量之一CT,可以由即张量缩并运算获得,其特性是在坐标旋转时,一直保持值不变,当张量分量存在标度因子误差时,旋转时所获得的CT会形成明显的波动,因此旋转过程中CT可以用来校正张量分量的标度因子误差,
把张量分量的校正模型重新表达为:
将公式(16)中的puv都取为1,计算出此时的CT值并进行平均处理。最后,CT的平均值和作为公式(16)中已知量,采用广义最小二乘拟合可估计出puv;
步骤七、结合高空校正获得的各项系数,可以求得每个张量分量的最终校正参数:Pnu0,Pmu0,Qnuj,Qmuj,Ouv;
步骤八、将三轴磁通门磁场测量值和校正参数绝对量输入到磁梯度张量分量校正模型中,计算校正后的磁梯度张量。
有益效果:1.与现有磁梯度张量仪实验室校正方法相比,本发明首先不用在野外构造磁屏蔽环境及应用昂贵的校正设备,其次本校正方法中数据采集的环境与实际航空磁测的环境完全一致,克服了实验室校正方法由于环境不一致造成的校正参数变化的问题。2.与现有磁梯度张量仪野外校正方法相比,本发明提出的校正方法多了在地面利用张量旋转不变量进行校正的过程。此过程是把张量各个分量形成一个整体后进行校正,不但使得单个分量更加准确,也保证了张量作为一个整体的正确性。3.本发明应用中不是在地面寻找磁场均匀区域,而是寻找磁场非均匀区域。由于航空磁测的主要应用是发现铁矿等磁性异常,因此在航空磁测的目标区域中或附近一般难以存在磁场均匀区域。故本发明的使用条件更符合实际的航空磁测探测区域的地质条件,即受环境限制的因素更少。未校正前,均方根误差为785.97,校正之后均方根误差为0.57,改善比为1380.58。其它张量分量的校正效果类似,从而验证了本校正方法的正确性和显著的效果。
附图说明:
图1为航空器挂载式磁通门磁梯度张量仪的校正方法流程图
图2十字全张量仪安装示意图
图3张量分量之一gzz在校正前后与真实值的对比图
图4测线上校正前和校正后gzx和gxz的对比图。
建立“十字”型磁通门张量仪,选取第一个磁通门1坐标系为o0x0y0z06,第二个磁通门2坐标系为o1x1y1z17,第三个磁通门3坐标系为o2x2y2z28,第四个磁通门4坐标系为o3x3y3z39,惯性导航系统(INS)5的坐标系为oxyz 10,磁通门1和磁通门3放置在坐标系10的x轴上关于原点o对称,间距为Δx,采用无磁性材料连接。磁通门2和磁通门4放置在坐标系10的z轴上关于原点o对称,间距为Δz,采用无磁性材料连接。惯性导航系统5放置在y轴上,距离原点o间距为Δy,采用无磁性材料与其它部件连接。安装时尽量保证磁通门1坐标系o0x0y0z06、磁通门2坐标系o1x1y1z17、磁通门3坐标系o2x2y2z28和磁通门4坐标系o3x3y3z39与惯性导航系统5坐标系oxyz 10相同。上述连接构成“十字”型磁通门张量仪。
具体实施方式:
下面结合附图和实施例对本发明做进一步的详细说明。
本发明利用高空中地球磁场的梯度张量近似为零的特点进行校正,然后根据地面磁梯度非零区域中张量仪旋转时不变量的可以指示各个梯度分量标度因子不一致的特点,计算出高空校正后未确定的各个梯度分量的标度因子,最终达到校正磁通门张量仪的目的。
本发明以单个磁通门误差模型为理论基础,推导出磁通门张量仪的误差模型,然后利用高空中地球磁场的梯度张量近似为零的特点,以及张量形成的不变量可以指示各个梯度分量标度因子不一致的特点,针对四个磁通门按“十”字结构排列所构成的磁梯度张量仪,提供一种由航空器挂载后进行测量应用时的校正方法。
航空挂载式磁通门磁梯度张量仪中磁通门的测量误差来源于以下几项:磁通门传感器三个感应轴标度因子不完全一致造成的误差、三个感应轴不严格正交造成的误差、零点存在偏移、由机械安装误差造成的非对准误差、动态特性不一致造成的误差。
在本实施例中以“十字”型磁通门梯度张量仪为例,图2为“十字”型磁通门梯度仪安装结构图,其中包括安装方式(四个磁通门,惯性导航系统,刚性连接)和坐标系。
在具体应用实例中,可以通过仿真进行模拟。仿真中四个磁通门按十字形排列,基线距离为0.1米。地球磁场强度为50000nT,磁偏角为60°,磁倾角为-9°。
航空器挂载式磁通门磁梯度张量仪的校正方法,包括以下步骤:
步骤一、建立单个磁通门误差模型:
把惯性导航系统和磁通门张量仪进行刚体连接,把各个磁通门传感器坐标系x,y,z轴分别向惯性导航坐标系的x,y,z轴方向对齐,则单个磁通门误差模型为:
式中:F为标度因子误差矩阵;E为非正交误差矩阵;R为零位误差矢量;M为磁通门坐标系与惯性导航坐标系非对准误差矩阵;Dj为磁通门测量值对时间的j阶导数的描述,即动态特性描述矩阵;T为磁通门在自身测量坐标系下的实测值;B为磁通门在惯导坐标系下的真实磁场值;Φ为磁通门的测量噪声;
每个误差矩阵或矢量的具体表达式如下:
(1)标度因子误差矩阵:
其中Sx,Sy,Sz表示磁通门感应轴的标度因子误差,
(2)非正交误差矩阵:
假设实际磁通门坐标系x轴与理想磁通门坐标系x轴的方向一致,ρ为实际磁通门坐标系y轴与理想磁通门坐标系y轴的夹角;为实际磁通门坐标系z轴与理想磁通门坐标系z轴的夹角在磁通门xoz平面上的投影;λ为实际磁通门坐标系z轴与理想磁通门坐标系z轴的夹角在磁通门yoz平面上的投影,
(3)非对准误差矩阵:
式中cα=cosα,cβ=cosβ,cγ=cosγ,sα=sinα,sβ=sinβ,sγ=sinγ,α、β和γ分别为实际磁通门坐标系相对惯导坐标系,绕实际磁通门坐标系x轴,y轴和z轴的非对准角。当非对准角非常小的时候可以将非对准误差矩阵近似表示为:
(4)零位误差矢量:
式中μx,μy,μz分别为磁通门三个轴分别存在的零位误差,
(5)动态特性描述矩阵:
当磁通门由一阶系统描述时取j=1,D1为一阶导数前面的系数矩阵。
将所述三轴磁通门的误差模型改写为校正模型:
式中:C0为T第0阶导数的校正参数矩阵;Cj为Τ第j阶导数的校正参数矩阵;O为零位误差校正矢量;bx,by,bz分别表示磁通门在惯性导航坐标系x,y,z轴上的磁场值,则磁通门中单个磁场分量误差校正模型可表示为:
式中:u表示x,y,z中一个轴,Pu0为矩阵C0中对应bu的行向量校正参数;Quj为矩阵Cj中对应bu的行向量校正参数;Ou为总偏差矢量On中对应bu的校正参数;
步骤二、建立磁梯度张量分量误差校正模型:
磁场沿x,y,z三个方向的二阶导数可构成磁梯度张量G,其中只有五个量是独立,利用差分近似微分的方式表示磁梯度张量各个分量。对于“十字”型的磁通门磁梯度张量仪在惯性导航坐标系下磁梯度张量G可表示为:(公式也换了)
张量分量表示方式为:
式中:m,n表示磁通门标号,取1,2,3,4;u,v都代表坐标轴的方向,取x,y,z;guv为惯性导航坐标系下某个张量分量;ΔBu为惯性导航坐标系下沿u轴方向的磁场变化量;Δsv表示磁通门n与磁通门m之间的基线距离。
将所述张量分量误差模型改写为:
式中:Ouv为张量分量guv的总偏置误差;
步骤三、在高空中磁场均匀的区域采集校正数据:
首先,根据地磁场强度、磁偏角和磁倾角获得地磁场三个分量,把地磁场作为四个磁通门在大地坐标系下的真实值。其次,将装载张量仪的吊舱的姿态进行随机变化来模仿高空飞行过程中张量仪的姿态变化,计算四个磁通门相应的空间位置及姿态,结合地磁场值计算四个磁通门在磁通门坐标系下的真实值。然后,按照公式(1)及四个磁通门的误差系数获得每个磁通门测量值。误差系数如表1所示,其中动态特性由一阶导数的系数控制,这是由于实际磁通门可以等效为一个一阶系统。最后,用差分近似微分的方法,可以获得高空中带误差的张量数据,并以此作为高空中采集的所需校正数据。
表1仿真的误差系数
步骤四、因为高空中地球磁场的梯度近似为0,设张量分量guv误差模型中磁通门n中沿u轴方向测量值前面的校正参数为puv,将上述张量分量的误差校正模型改写为:
步骤五、在地面磁场不均匀的区域采集校正数据:
在地面上有磁梯度的区域,把张量仪的传感器部分绕其中心点进行旋转,记录此过程的测量值。地面张量数据模拟方法和高空类似,只有两点区别,一个是把磁偶极子产生磁场(公式17所示)叠加到地磁场上,模仿地面具有磁场梯度的情况。另一个是张量旋转不变量是通过随机改变旋转角后,再利用张量缩并公式计算而得。(此处为了和前面的公式编号保持一致,故编号未做更改)
式中,r为从磁源到某个位置的方位矢量,M为磁矩矢量。
步骤六、以不变量进行约束,求解张量分量的标度因子误差:
将高空校正之后所获得的各项校正系数重新带入到公式(13)中,可以获得各个张量分量的相对量值,表达为:
把张量分量的校正模型重新表达为:
将公式(16)中的puv都取为1,计算出此时的CT值并进行平均处理。最后,CT的平均值和作为公式(16)中已知量,采用广义最小二乘拟合可估计出puv。如表2所示,通过实际和求出的张量分量标度因子误差系数对比,两者的相对误差从0.137%到0.528%,验证不变量校正方法的正确性。
表2实际和求出的张量分量标度因子误差系数对比
pxx | pyx | pzx | pxz | pyz | pzz | |
理论值 | 0.9336 | 1.0201 | 0.9248 | 1.0522 | 0.9888 | 1.0976 |
解算值 | 0.9291 | 1.0187 | 0.9215 | 1.0551 | 0.9934 | 1.0918 |
相对误差 | 0.482% | 0.137% | 0.357% | 0.276% | 0.465% | 0.528% |
步骤七、结合高空校正获得的各项系数,可以求得每个张量分量的最终校正参数:Pnu0,Pmu0,Qnuj,Qmuj,Ouv;
步骤八、将三轴磁通门磁场测量值和校正参数输入到磁梯度张量分量校正模型中,计算校正后的磁梯度张量。为模拟校正前后的数据,首先按航空磁测的方式以地下10米的一个磁偶极子作为探测目标,在目标正上方40米高度,形成一条400米的测线。然后利用公式(18)计算出磁偶极子在测线上形成的理论张量。测量的张量数据由类似高空数据模拟方法生成,只是叠加上公式(17)计算的磁偶极子的磁场。
式中,δij为Kronecker’s delta,i,j=1,2,3表示在笛卡尔坐标系中x,y,z。
本发明方法的具体实施步骤如图1所示,实施结果如图3所示为测线上的理想、未校正和校正后张量分量之一gzz的对比。
Claims (1)
1.一种航空器挂载式磁通门磁梯度张量仪的校正方法,其特征在于,包括以下步骤:
步骤一、建立单个磁通门误差模型:
把惯性导航系统和磁通门张量仪进行刚性连接,并把各个磁通门传感器坐标系x,y,z轴分别向惯性导航坐标系的x,y,z轴方向对齐,则单个磁通门误差模型为:
式中:F为标度因子误差矩阵;E为非正交误差矩阵;R为零位误差矢量;M为磁通门坐标系与惯性导航坐标系非对准误差矩阵;Dj为磁通门测量值对时间的j阶导数的描述,即动态特性描述矩阵;T为磁通门在自身测量坐标系下的实测值;B为磁通门在惯导坐标系下的真实磁场值;Φ为磁通门的测量噪声;
每个误差矩阵或矢量的具体表达式如下:
(1)标度因子误差矩阵:
其中Sx,Sy,Sz表示磁通门感应轴的标度因子误差,
(2)非正交误差矩阵:
假设实际磁通门坐标系x轴与理想磁通门坐标系x轴的方向一致,ρ为实际磁通门坐标系y轴与理想磁通门坐标系y轴的夹角;为实际磁通门坐标系z轴与理想磁通门坐标系z轴的夹角在磁通门xoz平面上的投影;λ为实际磁通门坐标系z轴与理想磁通门坐标系z轴的夹角在磁通门yoz平面上的投影,
(3)非对准误差矩阵:
式中cα=cosα,cβ=cosβ,cγ=cosγ,sα=sinα,sβ=sinβ,sγ=sinγ,α、β和γ分别为实际磁通门坐标系相对惯导坐标系,绕实际磁通门坐标系x轴,y轴和z轴的非对准角,当非对准角非常小的时候将非对准误差矩阵近似表示为:
(4)零位误差矢量:
式中μx,μy,μz分别为磁通门三个轴分别存在的零位误差,
(5)动态特性描述矩阵:
当磁通门由一阶系统描述时取j=1,D1为一阶导数前面的系数矩阵,将上述三轴磁通门的误差模型改写为校正模型:
式中:C0为T第0阶导数的校正参数矩阵;Cj为Τ第j阶导数的校正参数矩阵;O为零位误差校正矢量;bx,by,bz分别表示磁通门在惯性导航坐标系x,y,z轴上的磁场值,则磁通门中单个磁场分量误差校正模型表示为:
式中:u表示x,y,z中一个轴,Pu0为矩阵C0中对应bu的行向量校正参数;Quj为矩阵Cj中对应bu的行向量校正参数;Ou为总偏差矢量On中对应bu的校正参数;
步骤二、建立磁梯度张量分量误差校正模型:
磁场沿x,y,z三个方向的二阶导数构成磁梯度张量G,尽管G有九个分量,由于在不包括场源的域内,磁场的散度和旋度为0,即对称,且迹为0,故只有五个分量是独立的,表示为:
磁通门张量仪由不同方式构建,主要有十字、四方体和三角形结构,但都是利用差分近似微分的方式来表示磁梯度张量G的各个分量,以十字型结构为例,在惯性导航坐标系下,张量分量表示方式为:
式中:m,n表示磁通门标号,取1,2,3,4;u,v都代表坐标轴的方向,取x,y,z;guv为惯性导航坐标系下某个张量分量;ΔBu为惯性导航坐标系下沿u轴方向的磁场变化量;Δsv表示磁通门n与磁通门m之间的基线距离,
将所述张量分量误差模型改写为:
式中:Ouv为张量分量guv的总偏置误差;
步骤三、在高空中磁场均匀的区域采集校正数据:
将磁通门磁梯度张量仪挂载到航空器直升机下方较远处,航空器飞行到高空,然后以方形,八字形多种轨迹进行正常飞行,张量仪采集不同姿态下带误差的磁场测量数据;
步骤四、因为高空中地球磁场的梯度近似为0,设张量分量guv误差模型中磁通门n中沿u轴方向测量值前面的校正参数为puv,将上述张量分量的误差校正模型改写为:
故式(13)中有一个测量值的校正参数为1,这样以此项测量值为已知值,对式(13)列写线性方程组,通过最小二乘拟合方法解算线性方程组,获得式(13)中的各项参数,此时解算出各个磁梯度张量分量误差模型中参数的相对量,
步骤五、在地面磁场不均匀的区域采集校正数据:
在地面上有磁梯度的区域,把张量仪的传感器部分绕其中心点进行旋转,记录此过程的测量值;
步骤六、以不变量进行约束,求解张量分量的标度因子误差:
将高空校正之后所获得的各项校正系数重新带入到公式(13)中,获得各个张量分量的相对量值,表达为:
由于式(1)中误差矩阵的数值并不相同,因此每个分量所对应的puv也各不相同,故此时的张量分量认为存在张量分量标度因子误差,张量的不变量之一CT,由即张量缩并运算获得,其特性是在坐标旋转时,一直保持值不变,当张量分量存在标度因子误差时,旋转时所获得的CT会形成明显的波动,因此旋转过程中CT用来校正张量分量的标度因子误差,
把张量分量的校正模型重新表达为:
则张量不变量校正模型表达为:
将公式(16)中的puv都取为1,计算出此时的CT值并进行平均处理,最后,CT的平均值和作为公式(16)中已知量,采用广义最小二乘拟合估计出puv;
步骤七、结合高空校正获得的各项系数,求得每个张量分量的最终校正参数:Pnu0,Pmu0,Qnuj,Qmuj,Ouv;
步骤八、将三轴磁通门磁场测量值和校正参数输入到磁梯度张量分量校正模型中,计算校正后的磁梯度张量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610102276.XA CN105891755B (zh) | 2016-02-25 | 2016-02-25 | 航空器挂载式磁通门磁梯度张量仪的校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610102276.XA CN105891755B (zh) | 2016-02-25 | 2016-02-25 | 航空器挂载式磁通门磁梯度张量仪的校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105891755A CN105891755A (zh) | 2016-08-24 |
CN105891755B true CN105891755B (zh) | 2018-07-03 |
Family
ID=57014291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610102276.XA Expired - Fee Related CN105891755B (zh) | 2016-02-25 | 2016-02-25 | 航空器挂载式磁通门磁梯度张量仪的校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105891755B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106353824B (zh) * | 2016-09-29 | 2018-03-27 | 吉林大学 | 航空磁通门磁梯度张量仪的系统校正及磁干扰补偿融合方法 |
CN106525079B (zh) * | 2016-11-29 | 2019-03-29 | 北京科技大学 | 一种三轴磁传感器标定方法及装置 |
CN106997035A (zh) * | 2017-03-29 | 2017-08-01 | 吉林大学 | 一种基于磁梯度不变量的磁梯度计校正方法 |
CN108333551B (zh) * | 2018-02-14 | 2021-02-23 | 中国科学院电子学研究所 | 一种磁力计的校正方法 |
CN108828471B (zh) * | 2018-06-11 | 2020-02-07 | 中国科学院地质与地球物理研究所 | 一种多分量海底磁场测量方法及装置 |
CN109615845B (zh) * | 2018-12-24 | 2022-07-29 | 中国船舶重工集团公司第七一九研究所 | 一种声电磁集成探测与通讯一体化缆阵 |
CN109870153B (zh) * | 2019-03-28 | 2020-07-31 | 北京大学 | 一种磁强计正交性标定测试方法及标定测试装置 |
CN110007350B (zh) * | 2019-04-15 | 2020-08-25 | 哈尔滨工业大学 | 一种磁探测方法盲区的分析方法 |
CN110146839B (zh) * | 2019-05-30 | 2021-11-16 | 中国海洋大学 | 一种移动平台磁梯度张量系统校正方法 |
CN112444888A (zh) * | 2019-09-03 | 2021-03-05 | 中国石油天然气集团有限公司 | 地面磁力梯度测量方法、装置及系统 |
NO20220584A1 (en) | 2019-11-14 | 2022-05-16 | Baker Hughes Oilfield Operations Llc | Variation h coils calibration method for triaxial magnetometers |
CN110967041B (zh) * | 2019-12-18 | 2021-09-14 | 自然资源部国土卫星遥感应用中心 | 基于张量不变理论的卫星引力梯度数据精度的验证方法 |
CN111239667B (zh) * | 2020-03-16 | 2021-07-30 | 吉林大学 | 一种各阶磁梯度张量仪的统一校正方法 |
CN112362048A (zh) * | 2020-11-12 | 2021-02-12 | 中国电子科技集团公司第三研究所 | 一种实用的磁梯度张量高精度单点定位方法 |
CN113176529B (zh) * | 2021-03-31 | 2023-10-24 | 西南民族大学 | 一种磁梯度仪校正方法、磁梯度仪校正装置及存储介质 |
CN113156355B (zh) * | 2021-03-31 | 2022-11-08 | 吉林大学 | 一种超导全张量磁梯度测量装置的磁干扰补偿方法 |
CN113253162B (zh) * | 2021-06-18 | 2022-04-26 | 上海交通大学 | 一种微机电系统磁通门地磁张量传感芯片 |
CN113567902B (zh) * | 2021-07-02 | 2024-02-13 | 中国船舶重工集团公司第七一五研究所 | 一种矢量磁力梯度仪转向差校正方法 |
CN114415073B (zh) * | 2022-03-29 | 2022-08-09 | 中国人民解放军海军工程大学 | 一种航磁矢量梯度仪误差模型的地面快速校准方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103245928A (zh) * | 2013-05-23 | 2013-08-14 | 中国科学院上海微系统与信息技术研究所 | 方向可调的均匀磁场和均匀一阶梯度磁场的方法及装置 |
CN103743400A (zh) * | 2013-09-05 | 2014-04-23 | 新乡学院 | 一种磁航向解算装置 |
CN103837845A (zh) * | 2014-01-22 | 2014-06-04 | 哈尔滨工程大学 | 一种飞行器磁干扰场模型参数解算方法 |
CN104199115A (zh) * | 2014-09-18 | 2014-12-10 | 中国人民解放军国防科学技术大学 | 一种基于线性模型的地磁矢量测量误差综合补偿方法 |
CN104237958A (zh) * | 2014-09-29 | 2014-12-24 | 陕西宝成航空仪表有限责任公司 | 基于轨迹约束的地磁场矢量测量误差修正方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3000226B2 (ja) * | 1990-04-20 | 2000-01-17 | セイコーインスツルメンツ株式会社 | 校正用コイル付squid磁束計 |
US8942945B2 (en) * | 2011-04-19 | 2015-01-27 | General Electric Company | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
-
2016
- 2016-02-25 CN CN201610102276.XA patent/CN105891755B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103245928A (zh) * | 2013-05-23 | 2013-08-14 | 中国科学院上海微系统与信息技术研究所 | 方向可调的均匀磁场和均匀一阶梯度磁场的方法及装置 |
CN103743400A (zh) * | 2013-09-05 | 2014-04-23 | 新乡学院 | 一种磁航向解算装置 |
CN103837845A (zh) * | 2014-01-22 | 2014-06-04 | 哈尔滨工程大学 | 一种飞行器磁干扰场模型参数解算方法 |
CN104199115A (zh) * | 2014-09-18 | 2014-12-10 | 中国人民解放军国防科学技术大学 | 一种基于线性模型的地磁矢量测量误差综合补偿方法 |
CN104237958A (zh) * | 2014-09-29 | 2014-12-24 | 陕西宝成航空仪表有限责任公司 | 基于轨迹约束的地磁场矢量测量误差修正方法 |
Non-Patent Citations (2)
Title |
---|
三分量磁通门传感器非正交性误差校正;焦秉刚 等;《现代电子技术》;20110701;第34卷(第13期);全文 * |
基于差分的磁偶极子单点张量定位方法;李光 等;《探测与控制学报》;20121026;第34卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105891755A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105891755B (zh) | 航空器挂载式磁通门磁梯度张量仪的校正方法 | |
CN106353824B (zh) | 航空磁通门磁梯度张量仪的系统校正及磁干扰补偿融合方法 | |
Fang et al. | A novel calibration method of magnetic compass based on ellipsoid fitting | |
CN107272069B (zh) | 基于磁异常梯度的磁性目标追踪方法 | |
CN107544042B (zh) | 一种磁力计阵列校正方法 | |
CN109856689B (zh) | 一种超导航磁梯度张量数据抑噪处理方法和系统 | |
CN106997035A (zh) | 一种基于磁梯度不变量的磁梯度计校正方法 | |
CN110146839A (zh) | 一种移动平台磁梯度张量系统校正方法 | |
Gang et al. | Linear calibration method of magnetic gradient tensor system | |
Sui et al. | Correction of a towed airborne fluxgate magnetic tensor gradiometer | |
Li et al. | Integrated compensation and rotation alignment for three-axis magnetic sensors array | |
CN113156355B (zh) | 一种超导全张量磁梯度测量装置的磁干扰补偿方法 | |
CN107132587A (zh) | 航空超导全张量磁梯度测量系统安装误差标定方法及装置 | |
Sui et al. | Correction and compensation of an airborne fluxgate magnetic tensor gradiometer | |
Pang et al. | Misalignment error suppression between host frame and magnetic sensor array | |
CN109725361A (zh) | 基于磁梯度张量不变量的一种磁性目标定位方法 | |
CN111220932B (zh) | 无人机磁干扰标定方法及分布式磁异常探测系统 | |
CN109633540B (zh) | 一种磁源的实时定位系统及实时定位方法 | |
Pang et al. | Misalignment calibration of geomagnetic vector measurement system using parallelepiped frame rotation method | |
Gang et al. | Integrated calibration of magnetic gradient tensor system | |
Mu et al. | Small Signal Magnetic Compensation Method for UAV-Borne Vector Magnetometer System | |
Huang et al. | Two-step complete calibration of magnetic vector gradiometer based on functional link artificial neural network and least squares | |
Jiang et al. | Scalar calibration of aeromagnetic data using BPANN and LS algorithms based on fixed-wing UAV platform | |
Sui et al. | Analysis and simulation of flight effects on an airborne magnetic gradient tensor measurement system | |
Pang et al. | A new misalignment calibration method of portable geomagnetic field vector measurement system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180703 |
|
CF01 | Termination of patent right due to non-payment of annual fee |