[发明专利]一种磁体的二阶欧拉反演法

摘要

本发明提供了一种磁体的二阶欧拉反演法，将13个标量磁强计构成标量磁强计阵列，并安装GPS或高精度惯性导航系统，确定磁体的搜索区域，在每航迹段计算磁体的位置坐标和构造指数的估计值及总估计值的均值，再计算均值的最邻近整数值，得到磁体构造指数的反演值，对更新后的第K航迹段上的磁体位置坐标估计值进行解中心计算与迭代筛选，得到磁体位置坐标的反演值。本发明便于磁强计阵列装配和平面单元替换，实现了由多个测量点处的一阶和二阶磁梯度值同时反演出磁体的位置坐标与构造指数，消减了磁梯度的粗大测量误差对反演精度的影响，也免除了构造指数的主观经验选取和梯度测量数据的中心化处理。

主权项

1.一种磁体的二阶欧拉反演法，其特征在于，具体的步骤为：步骤1.将13个标量磁强计按照如下方式构成标量磁强计阵列，在三维空间坐标内，第一磁强计(0)位于坐标(0,0,0)点，即为测量点，第二磁强计(1)位于坐标(‑0.3kl,‑0.4kl,0)，第三磁强计(2)位于坐标(‑0.5kl,‑0.2kl,0)，第四磁强计(3)位于坐标(0.5kl,0.2kl,0)，第五磁强计(4)位于坐标(0.3kl,0.4kl,0)，第六磁强计(5)位于坐标(0,‑0.3kl,‑0.4kl)，第七磁强计(6)位于坐标(0,‑0.5kl,‑0.2kl)，第八磁强计(7)位于坐标(0,0.5kl,0.2kl)，第九磁强计(8)位于坐标(0,0.3kl,0.4kl)，第十磁强计(9)位于坐标(‑0.3kl,0,‑0.4kl)，第十一磁强计(10)位于坐标(‑0.5kl,0,‑0.2kl)，第十二磁强计(11)位于坐标(0.5kl,0,0.2kl)，第十三磁强计(12)位于坐标(0.3kl,0,0.4kl)，并将其捷联于运载体上，运载体安装GPS或高精度惯性导航系统；步骤2.由磁体的磁异常大小和标量磁强计阵列敏感的最小梯度确定磁体的搜索区域，按“己”字形设置运载体的航迹；步骤3.每隔NΣ个测量点将运载体的航迹划分为K段，在每一个航迹段由标量磁强计阵列对地磁总场的测量值及GPS或高精度惯性导航系统对运载体位置坐标的测量值，计算一次磁体的位置坐标和构造指数的估计值；步骤4.计算磁体构造指数的所有估计值的均值，再计算均值的最邻近整数值，得到磁体构造指数的反演值；步骤5.由磁体构造指数的反演值，建立关于磁体位置坐标的另一超定方程组，再使用最小二乘法重新计算第K航迹段上的磁体位置坐标的估计值；步骤6.对更新后的第K航迹段上的磁体位置坐标估计值进行解中心计算与迭代筛选，得到磁体位置坐标的反演值。