一种线圈式矢量磁力仪地磁日变误差校正方法与流程

1.本发明涉及磁场测量领域，尤其涉及一种线圈式矢量磁力仪地磁日变误差校正方法。


背景技术：

2.地磁场是矢量场，定义地磁七要素：地磁总场f、水平分量h、北向分量x、东向分量y、垂直分量z、磁倾角i、磁偏角d来表述地磁场，其中水平分量和垂直分量组成的平面为磁子午面。对地磁要素的检测在地质灾害预测、地球科学以及空间科学等领域都有重要应用。因此高精度地磁要素数据对资源探索，自然灾害预测以及全球地磁场模型建立意义重大。
3.地磁要素测量通过磁力仪来实现，按照测量方式可分为总场磁力仪、分量磁力仪、磁向测量磁力仪、多参量磁力仪。总场磁力仪检测精度高，但是只能测量地磁场标量值，无法检测其他要素；分量磁力仪是矢量磁力仪，能够测量地磁场分量，以磁通门磁力仪为例，其可以测量地磁场的三分量，但是存在温漂和正交性误差，不能进行绝对测量；磁向测量磁力仪是di仪，通过磁通门和无磁经纬仪组合来测量磁偏角和磁倾角，但是di仪测量周期长、不能连续观测以及需要专门操作人员，观测人员不同也会引入误差的不足；多参量磁力仪是通过总场传感器和线圈组合成线圈式磁力仪实现对地磁要素测量，主要有fhd磁力仪和didd磁力仪。受太阳周期性活动的影响，地磁矢量场的强度和方向一直连续变化，这种以24小时为周期的地磁场变化现象称之为地磁日变化，而日变干扰是高精度地磁测量时的主要干扰，因为日变干扰是持续存在的，而线圈式磁力仪的检测原理是基于每个测量周期中所有量检测完成之后才能够解算相应的值，因此对检测过程中日变干扰的影响消除是提高系统检测精度的重要工作。
4.目前常用的日变干扰的抑制方法是通过建立日变观测站，用高精度的磁力仪同时进行测量，然后把两者做差，得到的差值就是日变量，得到日变之后再进行相关消除工作，但是该方法无法对磁力仪测量过程中的日变干扰抑制，并且目前很少有针对线圈式磁力仪检测过程中日变误差的分析及抑制方法研究。


技术实现要素：

5.为了解决上述问题，本发明提供了一种线圈式矢量磁力仪地磁日变误差校正方法，主要包括以下步骤：
6.s1：架设线圈式磁力仪和三分量磁力仪，然后记录线圈式磁力仪施加大小相等方向相反的偏置场后的叠加场，记为f
1
和f2‑
，以及不加偏置场时的地磁总场强度，记为f3；
7.s2：结合三分量磁力仪检测的地磁场分量信息，利用向量变换关系将所有时刻得到的叠加场都校准至测量f3的时刻，记为f
3
′
和f3‑
′
；
8.s3：根据校正后的叠加场f
3
′
和f3‑
′
，进而得到日变校正之后的地磁参量信息。
9.进一步地，若f
1
为第一个时刻测量得到的叠加场，f2‑
为第二个时刻测量得到的叠加场，f3为第三个时刻测量得到的叠加场，具体的校准过程如下：
10.记向量f
1
校正后为f
3
′
，得f
3
′
的计算如式(5)：
11.f
3
′
＝a

f3ꢀꢀ
(5)
12.代入f3与f1的关系式f3′
＝f1 δf
13
后，式(5)变换为式(6)：
13.f
3
′
＝f
1
δf
13
ꢀꢀ
(6)
14.将式(6)中等号两边开方可得式(7)
15.f
3
′
＝(f
1 2
2f
1
δf
13
δf
132
)
0.5
ꢀꢀ
(7)
16.其中，f
1
＝a

f3‑
δf
13
，将f
1
代入式(7)，化简后可得式(8)
17.f
3
′
＝[f
21
2δf
13
(f3 a

)
‑
|δf
13
|2]
0.5
ꢀꢀ
(8)
[0018]
同理可得f3‑
′
的计算如式(9)
[0019]
f3‑
′
＝[f
22
‑
2δf
23
(f3 a
‑
)
‑
|δf
23
|2]
0.5
ꢀꢀ
(9)
[0020]
其中，f
3
′
和f3‑
′
分别为f
1
和f2‑
校准后对应到第三时刻的叠加场。
[0021]
进一步地，日变校正之后的磁倾角δi
′
及磁偏角δd
′
的计算公式如下：
[0022][0023]
δi
′
＝θ
′ꢀꢀ
(11)
[0024][0025]
其中，f
3
′
和f3‑
′
分别为f
1
和f2‑
校准后对应到第三时刻的叠加场，f3为第三个时刻测量得到的地磁总场。
[0026]
本发明提供的技术方案带来的有益效果是：使用地磁分量以及施加特定偏置场后叠加场的数据就可以实现日变效应误差校正，无需其他参数；对线圈式磁力仪测量周期内地磁日变效应进行抑制，提高仪器检测精度；地磁日变误差抑制方法是基于线圈式磁力仪本身的操作流程，无需额外的操作步骤。
附图说明
[0027]
下面将结合附图及实施例对本发明作进一步说明，附图中：
[0028]
图1是本发明实施例中线圈式磁力仪测量原理图。
[0029]
图2是本发明实施例中线圈式磁力仪日变干扰机理图。
[0030]
图3是本发明实施例中线圈式磁力仪日变校正方法的流程图。
[0031]
图4是本发明实施例中磁倾角测量示意图。
[0032]
图5是本发明实施例中磁偏角测量示意图。
[0033]
图6是本发明实施例中日变校正前测量结果误差仿真图。
[0034]
图7是本发明实施例中日变校正后测量结果误差仿真图。
具体实施方式
[0035]
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
[0036]
线圈式磁力仪的测量原理如图1所示，在每个测量周期内，认为地磁场不变，然后分别施加大小相等方向相反的偏置场a

和a
‑
后，测量偏置场与地磁场叠加之后的磁场f
3
和
f3‑
，不施加偏置场时测得地磁总场f3，然后通过以下解算方法可以计算得到θ。图1中的f0是仪器开始测量时的地磁场量。
[0037]
在f
3
、a

、f3和f3‑
、a
‑
、f3，组成的两个三角形中，分别由余弦定理可得式(1)和式(2)：
[0038]
f
3 2
＝f
32
a
2
‑
2f
3
a

cos(α)
ꢀꢀ
(1)
[0039]
f3‑2＝f
32
a
‑2‑
2f3‑
a
‑
cos(π
‑
α)
ꢀꢀ
(2)
[0040]
对式(1)和式(2)做差，同时代入a＝a
3
＝a3‑
，α为一个辅助角，通过可得θ(日变校正前记为θ，经过日变校正后计算得到的值记为θ
′
)的计算公式如式(3)：
[0041][0042]
其中，a的表达式为式(4)：
[0043][0044]
对于磁倾角的测量结果为δi＝θ，对于磁偏角的测量则为
[0045]
然而在实际测量时，每一次施加偏置场时地磁场在日变干扰下发生变化，就出现图2中所示的情况，此时检测得到的叠加磁场为f
1
和f2‑
，与图1中的叠加场f
3
和f3‑
相比发生了偏差，从而导致最终解算的磁倾角δi及磁偏角δd存在误差。
[0046]
为解决上述仪器检测过程中的日变干扰，本发明的实施例提供了一种线圈式矢量磁力仪地磁日变误差校正方法，利用三分量磁力仪检测的地磁场三分量信息，通过适量关系转换将所有测量叠加场时刻的值都校准至同一时刻然后再进行地磁参量信息解算，实现对日变误差的抑制，校正示意图如图3所示。假设实际检测时f
1
为第一个时刻测量得到的叠加场，f2‑
为第二个时刻测量得到的叠加场，f3为第三个时刻测量得到的叠加场，需要做的是将f
1
和f2‑
都校准至第三时刻，校准后如图2中f
3
′
和f3‑
′
所示。具体实施步骤如下：
[0047]
记向量f
1
校正后为f
3
′
，f
3
′
的计算如式(5)：
[0048]
f
3
′
＝a

f3ꢀꢀ
(5)
[0049]
代入f3′
与f1的关系式f3′
＝f1 δf
13
后，式(5)变换为式(6)：
[0050]
f
3
′
＝f
1
δf
13
ꢀꢀ
(6)
[0051]
将式(6)中等号两边开方可得式(7)
[0052]
f
3
′
＝(f
1 2
2f
1
δf
13
δf
132
)
0.5
ꢀꢀ
(7)
[0053]
其中，f
1
＝a
3
f3‑
δf
13
，将f
1
代入式(7)，化简后可得式(8)
[0054]
f
3
′
＝[f
21
2δf
13
(f3 a

)
‑
|δf
13
|2]
0.5
ꢀꢀ
(8)
[0055]
同理可得f3‑
′
的计算如式(9)
[0056]
f3‑′
＝[f
22
‑
2δf
23
(f3 a
‑
)
‑
|δf
23
|2]
0.5
ꢀꢀ
(9)
[0057]
此时可以得到日变校正之后的磁倾角δi
′
及磁偏角δd
′
，计算过程如下：
[0058]
计算公式(10)：
[0059][0060]
δi
′
＝θ
′ꢀꢀ
(11)
[0061][0062]
通过式(8)和式(9)可知，只要得到δf
jk
、a(所施加偏置场矢量)和f3的向量表示，然后进行运算即可解算出f
3
′
和f3‑
′
，从而实现日变误差的校正。
[0063]
向量δf
jk
可以表示为式(13)，式(13)中的x
j
，y
j
，z
j
和x
k
，y
k
，z
k
分别为f
j
和f
k
时刻地磁场的分量值，可以利用地磁分量测量传感器同步检测得到。
[0064]
δf
jk
＝(x
k
‑
x
j
,y
k
‑
y
j
,z
k
‑
z
j
)
ꢀꢀ
(13)
[0065]
其中，δf
jk
用来表示式(8)或式(9)中δf
13
和δf
23
。
[0066]
在仪器完成初始定向后，仪器位置保持不变，因此所施加偏置场向量方向不变，偏置场坐标需要在初始时刻进行解算，磁倾角测量时，初始时刻各个向量的分布如图4所示。
[0067]
此时可以解算得到a
i
和a
i
‑
在xyz坐标系中的向量表示，分别如式(14)、(15)所示：
[0068]
a
i
＝(
‑
a
i
sin(i0)cos(d0),
‑
a
i
sin(i0)sin(d0),a
i
cos(i0))
ꢀꢀ
(14)
[0069]
a
i
‑
＝(a
i
sin(i0)cos(d0),a
i
sin(i0)sin(d0),
‑
a
i
cos(i0))
ꢀꢀ
(15)
[0070]
其中，a
i
和a
i
‑
分别表示测倾角时施加的等大反向的偏置场，计算θ
′
时公式(8)中的a在测量磁倾角时a＝a
i
＝a
‑
，测量磁倾角时a
i
就对应a，a
i
为c
i
线圈施加的偏置场的标量值，i0和d0分别为初始时刻的磁倾角和磁偏角基值。
[0071]
同理，对磁偏角的测量矢量如图5所示：
[0072]
此时偏置场a
d
与水平面平行，根据几何关系解算得到a
d
和a
d
‑
在xyz坐标系中的向量表示，分别如式(16)、(17)所示：
[0073]
a
d
＝(
‑
a
d sin(d0),a
d cos(d0),0)
ꢀꢀ
(16)
[0074]
a
d
‑
＝(a
d sin(d0),
‑
a
d cos(d0),0)
ꢀꢀ
(17)
[0075]
其中，a
d
和a
d
‑
表示测量磁偏角时所加等大反向的偏置场。
[0076]
由此得到所有计算所需量，然后代入日变校正方法中，完成对线圈式磁力仪地磁日变干扰误差校正。
[0077]
最终对本专利所提出的线圈式磁力仪地磁日变误差校正方法的校正效果进行仿真：设定最大日变干扰为0.45nt，此时计算日变效应对线圈式磁力仪测量结果引入的最大误差超过5"，然后加入日变校正方法后再计算线圈式磁力仪测量结果误差，其误差为0"，仿真图分别如图6和图7所示。
[0078]
本发明的有益效果是：仅用地磁分量信息以及线圈式磁力仪的测量数据就可以实现线圈式磁力仪测量周期中日变效应引入误差的校正，无需其他参数；对线圈式磁力仪测量周期内地磁日变效应进行抑制，提高了仪器检测精度；地磁日变误差抑制方法是基于线圈式磁力仪本身的操作流程，无需额外的操作步骤。
[0079]
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。