S?170?(Bi?1ni?B)2n?1?????????(A4)
1n式中：B??Bi，n为参加统计的观测值数。
ni?1A1．4 三种测定噪声方法的比较
为了进行比较，将3台微机质子磁力仪按上述要求在同一地点同时测定仪器的噪声水平(达到秒一级同步观测)，而后用3种方法计算各台仪器的噪声水平如表A1： 4阶差分包络值的3台仪器同时测定法 磁力仪类型及号码 S值 与平均值的 nT 百分误差％ OMNI—4 336 MP—4 316 MP—4 0．16 0．17 0．24 14 11 11 S值 nT 0．13 0．18 0．30 与平均值的 百分误差，％ 7 5 11 7点滑动平均法 1/70法 S值 nT 0．12 0．21 0．27 与平均值的 百分误差，％ 14 11 0 由表可见3种方法算出的噪声水平是很接近的，其差异在统计涨落范围之内，而且3种方法表示仪器噪声水平的次序不变。
A2 测定磁力仪的一致性 A2．1 测定探头一致性
经验表明，制作探头与夹固探头的各种材料的“磁清洁”程度有差异，是造成一致性误差的主要因素。因此，在每测区开工前，要对所有探头的一致性进行测定，方法如下：
首先将成套仪器所配探头(一般是5个)编上号，然后用两台仪器作秒级同步日变观测。其中台站型仪器及一个探头固定不变，即以此为准进行比较。另一台仪器分别轮换同其余4个探头相联结，并注意换探头时主机不能关机，各探头位置应尽量一致，调谐场值预先选好保持不变。每个探头读数三十余次以上，而后分别求出相应与台站仪器读数的差值，并计算各差值数组的算术平均值，比较这4个平均值，即可判断探头的一致性。如以某次测定为例，得出5个探头的一致性如表A2： 探头编号 平均值，nT 1 0．4 2 1．3 3 0 4 0．6 5 1．1 由表A2可见l，3，4号探头一致性较好，而2号与5号探头一致性较好，可以配对使用。使用中不能随意调换探头。以免引入系统误差。 A2．2 校验主机的一致性
从原理上说，质子磁力仪的主机就是一个用来测定核子旋进频率的测频器，而当前测频精度是很高的，主机能以0．00425Hz的分辨率来精确测量频率，出
D-15
厂时都用精确度更高的讯号发生器进行校准，并保证绝对准确度50000nT时为±1nT。所以，一般情况下，主机的一致性都能符合要求。
为校验主机的—致性，可使用同一探头，用不同主机轮换作日变观测，使每台主机读数20～30次，将整个测量段的日变曲线绘出，察看曲线变化趋势是否有脱节现象。若曲线“圆滑”，即表明主机的一致性良好。
附 录 B
高精度磁测的各项改正与基站To值测定工作
(补充件)
B1 正常梯度改正与高度改正
当进行大面积高精度磁测工作时，需要进行正常梯度改正。此时若仍沿用查全国地磁图的办法作正常场梯度改正。就不能满足精度要求了。此时要用国际地磁参考场IGRF1990．0模型提供的高斯系数，用电子计算机算出测区内1km×1km节点地磁场To值。而后以1nT的间距绘制To等值线图。用此图作正常场梯度改正，其作法是以通过总基点的等恒线为零线，向北每过一条等值线减少1nT，向南每过一条等值线增加1nT，以此类推。 B1．1 地磁场各分量的球谐表达式：
mmX???(gncosm??hnsinm?)?n?1m?0n?ndmPn(cos?) d?B1 Y????nmmm(gnsinm??hncosm?)?Pnm(cos?)
n?1m?0sin??mmZ????(n?1)(gncosm??hnsinm?)?Pnm(cos?)
n?1m?0mm 式中：gn、hn——高斯系数，可查表求得；
s)——n次m阶缔合勒让德多项式； Pnm(co??——各节点的经度值；
?——各节点的余纬度。
可见将高斯系数与各节点的坐标值代入(B1)式，即可求出各节点的正常磁场值。
B1．2 国际地磁参考场(IGRF)1990．0模型
1968年10月在华盛顿召开的“地球基本磁场的描述”会议上，通过了1965．0年代国际地磁参考场(IGRF)。后经国际地磁学和高空物理学会执行委员会及世界地磁测量部的同意，作为世界通用的主要磁场标准，其后对国际地磁参考场又作了三次修改补充。DGRF表示确定的地磁参考场，其高斯系数今后不再修改，而每五年改变一次模型，即通过年变率的调整取得。如IGRFl990．0包括1990年代的地磁场模型及一个预测的长期变化模型，即主要磁场模型(m＝n＝10)前80个高斯系数的年变率，用来调整1990．0～1995．0年期间的地磁场模型。如表B1所示。
B 1．3 计算磁场总强度To及梯度值：
D-16
To?X2?Y2?Z2?H2?Z2?????????(B2)
在一级近似的情况下，沿南北向的磁场梯度均半径，R＝6371000 m。
沿垂向的磁场梯度：
?To3ZH?，式中R为地球平?X2RTo?To3T??o?????????(B3) ?RRB1．4 高度改正
由(B3)式可知，当To＝50000nT时，地磁场垂向梯度为－0．024nT／m。高差30 m时，地磁场垂向变化可达－0．72nT。因此本标准要求测点都要测出或从地形图上读出高程，高程误差依本标准第5．5表2执行。
高度改正从总基点高程起算，以To＝50 000nT为例，约每42m高差改正1nT，比总基点高42m时加1nT，比总基点低42m时减1nT。 B2 日变改正
在地球的变化磁场中，除大家熟悉的长期变化和静日变化Sq外，还有源于空间磁场的扰动变化，包括周期从0．2～1000 s的快速振动，其幅度—般为0．01～5nT(有时可达10～20nT)。高精度磁测必需对地磁场的这种短周期变化进行改正，因此对日变观测有较高的要求。再者。目前使用的质子磁力仪可测定绝对地磁场强度，无需使用基点网进行地磁场值的传递，但需消除地磁日变的影响，求出各测站之间地磁场的真正差值。由此可见，高精度磁测的日变观测有两个作用；其一是对测网观测值进行日变改正，其二是在各站之间进行同步日变观测，以选择地球扰动磁场振幅最小时间来确定测站之间地球基本磁场的差值。所以在高精度磁测工作中．主基点的作用被日变站取代了。 B2．1 日变站址的选择：
由于高精度磁测的日变站兼有基点的作用．需要参照对主基点的要求来选择日变站地址。具体要求按本标准第7．2．1．2条各条执行。
B2．2 日变站所用磁力仪的精度应与测网观测磁力仪的精度相同或更高。对微机质子磁力仪的读时准确性要达到秒级，读数间隔5～20 s，日变站与测网观测的同时性亦是秒级，对普通质子磁力仪应准确到0．5min，读数间隔0．5～1min。 B2．3 由于地磁场短周期变化的振幅与微机质子磁力仪噪声均方根值是近于同一数量级的，而噪声是随机的，地磁场脉动变化是有规律的。因此可对日变观测先作5点或7点滑动平均，压低噪声水平之后，再对测网观测进行日变改正，即可提高日变改正精度。尤其当要求日变改正精度优于0．5nT时，必须这样作。 B2．4 地磁场短周期变化受地磁感应场的影响特别明显(各种电流系统形成的地球空间磁场中，约包含30％的感应场)。而且周期愈短，所受影响愈大。因此必须充分重视地电结构不同对地磁短周期变化的巨大影响，这关系到日变站的控制范围。根据对地电结构相同与不同地区所作一系列日变对比试验结果可知，当磁测均方误差为2～5nT时，日变站的控制范围不应超过50 km，当要求磁测精度优于2nT时，在地电结构基本相同的情况下，日变站的控制范围不应超过30 km。当地电结构显著不同时，要增设日变站，其控制范围要经试验确定。
D-17
B 3 基站To值测定工作
用质子磁力仪在总基点与日变站测出的地磁场绝对值Ti是时间的函数，如下式：
Ti?To??(t)?????????(B4)
式中：To——该处地球基本磁场值，不随时间变化；
?(t) ——该处地球变化磁场，是时间的函数。
为准确求出总基点与日变站的To值。需作较长时间的日变观测，读数间隔不大于20 s，观测时间不短于2h。要选择地磁场变化平稳段，即2h内地磁场平均值变比不超过2nT的时间段，求取Ti的平均值Ti，Ti即为该处的To值。
1nTi?To???(t)?????????(B5)
ni?1式中：n——参与统计的地磁场Ti值总数，n＞100。
因为上式中?(t)具有随机性，在n＞l00的情况下，其统计平均值应趋于零。所以Ti?To。
顺便说明。先用分基To值作日变改正后再作总基改正，与直接使用总基To值作日变改正的结果是等效的。
D-18
D-19