关于地面高精度磁测精确改正方法探讨

关于地面高精度磁测精确改正方法探讨摘要：地面高精度磁测的精度要求已经达到15nT，传统的日变改正、经纬度改正、高度改正、总基点改正方法速度慢、精度低，已经无法满足生产需要。作者在实际工作中总结出了一套高效高精度的解析改正方法并编制了相应的软件。关键词：地面高精度磁测；高度改正；经纬度改正；误差中图分类号：P624献标识码：文章编号：2095-2104（2012）12-0020-021.引言众所周知，磁法测量数据需要进行多项改正：日变改正，消除日变影响；有多个基站时，需要进行总基点改正。测区范围较大时，需要进行经纬度改正；测区高差较大时，需要进行高度改正。各项改正后得到的结果，才能得到真实反映地下磁性体引起的磁异常T。上世纪六十年代前后，磁法工作以寻找地表或埋藏较浅的磁铁矿为主，对磁测数据的精度要求不高（一般20nT左右）。二十世纪末，逐步进入第二找矿空间，矿体埋深较大，在地表能够观测到的磁异常微弱，这就要求磁测数据有更高的精度，精细反映弱磁异常。表1是地面高精度磁测技术规程要求的磁测总精度及误差分配表。有诸多因素影响着磁测总精度，要提高磁测总精度，需要降低每一项的精度。表1磁测误差分配表随着科学技术的进步，质子磁力仪、光泵磁力仪等先进的仪器已经广泛应用。质子磁力仪的灵敏度可达0.01nT，精度可达0.1nT，是目前野外生产的主要机型。仪器一致性误差和噪声误差方面完全可以达到规范要求；操作及点位误差具有较大的随机性，主要依靠外业工作质量保障体系保证；各项改正的精度成为影响磁测总精度的重要因素。传统的改正方法已经不适应高精度磁测工作。2.日变改正和总基点改正在地球的变化磁场中，除大家熟悉的长期变化和静日变化外，还有源于空间磁场的扰动变化，包括周期从0.21000s的快速振动，其幅度般在020nT。磁法测量必需对地磁场的这种短周期变化进行改正，消除其影响。因此需要进行日变观测和日变改正。地磁场短周期变化受地磁感应场的影响特别明显(各种电流系统形成的地球空间磁场中，约包含30的感应场)。而且周期愈短，所受影响愈大。因此必须充分重视地电结构不同对地磁短周期变化的巨大影响，这关系到日变站的控制范围。根据对地电结构相同与不同地区所作一系列日变对比试验结果可知，当磁测均方误差为25nT时，日变站的控制范围不应超过50km，当要求磁测精度优于2nT时，在地电结构基本相同的情况下，日变站的控制范围不应超过30km。因此当测区范围较大、地电结构复杂时，需要设立多个日变观测站，并进行总基点改正。2.1传统的日变改正和总基点改正方法需要求出每个日变观测站的T0值。根据相关规范要求，日变观测站T0值的求取方法为：做较长时间的日变观测，读数间隔不大于20s。选择地磁场变化平稳段，即2h内地磁场平均值变比不超过2nT的时间段，求取Ti的平均值，即为该处的T0值。(1)式中：n参与统计的地磁场Ti值总数，n100。日变改正方法：用日变站的T0值减去日变观测值等于日变改正值：TrT0Tr（2）测点观测值加日变改正值等于日改后测点观测值：TcrTcTr（3）基点改正方法：以上已经计算出了测点的日改后磁场值Tcr，这个数据已经消除了日变影响，用这个值减去总基点的T0值即为测点的磁异常值：Tzj=TcrTzj0（4）公式中Tzj0为为总基点的T0值。由公式（2）和（4）可见，改正时用到了两个T0值：日变站T0和总基点TO。当测区只有一个日变站时，两个T0值相同；当有两个以上日变站时，如果T0值求的不准，就会产生系统误差。事实上，准确求取T0值是很困难的，不同日期、不同时间求取的T0值存在一定的误差。下面是某地日变站不同日期求取的T0值：表2不同日期求取的T0值不同日期求取的T0值误差可达4.13nT。可见，传统的改正方法需要求取T0值，改正速度低，当分基点较多时，这种改正方法必然产生较大的改正误差。2.2精确改正方法为了消除基点T0值求取不准产生的误差，改正分两步进行。第一步：日变改正和分基点改正一次性完成。根据（2）、（3）式，可知：TFjTcrT0TcTrT0TcT0TrT0TcTr（5）式中Tc是观测点磁场值，Tr是测点观测时刻的日变站磁场值，根据日变曲线用内插法求取，这个工作可由专门软件完成，同步时间可精确到秒，高精度质子磁力仪的日变观测时间间隔最小5s，因此，改正精度很高。第二步：进行总基点改正改正公式为：TZJTFjT0（6）式中T0总基点和分基点之间的基本磁场的差值。改正精度取决于总基点和分基点之间的基本磁场的差值。事实证明，总基点与各个分基点进行同步日变观测，求取总基点与分基点之间的基本磁场差值，精度非常高。为避免两台仪器可能存在系统误差带来的影响，联测时两台磁力仪还要互换进行观测，取两次的平均值作为分基点和总基点之间的基本磁场的差值，可以有效消除仪器系统误差带来的影响。表3是某矿调项目基点联测成果的标准离差统计成果。表3基点联测成果标准离差统计表表中14列由加拿大GSM-19T质子磁力仪观测求取；56列由捷克MPG-1质子磁力仪观测求取,为同一基点不同日期联测成果。可见，联测精度非常高，标准离差最高达0.05nT。在某矿调项目，磁测工作区面积较大，总共有5个分基点，表4列出了T0值法和基点联测法求取的总基点和分基点差值成果。表4T0值法和基点联测法成果对比表可见，传统基点改正的误差最高达4.615nT。3.经纬度改正和高度改正3.1传统经纬度改正方法：最传统的方法是：从地磁图查出测区地磁场值，以1nT的间距绘制T0等值线图，从总基点算起，向南改正值为正，向北改正值为负。线性分析法：为了提高改正精度，地球物理工作者又采用了线性分析法，用国际地磁参考场IGRF模型提供的高斯系数，用电子计算机算出测区内1km1km节点地磁场T0值。而后以1nT的间距绘制T0等值线图。用此图作经纬度改正，其作法是以通过总基点的等恒线为零线，向北每过一条等值线减少1nT，向南每过一条等值线增加1nT，依此类推。3.2常用的高程改正方法：高度改