影响海洋磁力测量的误差可以分为实时测量误差和非实时测量误差,实时测量误差包括GPS、磁力仪系统产生的测量误差。非实时测量误差包括海洋磁力测量实施过程中的误差、拖鱼定位误差和正常场校正误差,还包括固有的地磁日变改正、船磁改正、拖鱼位置归算误差的影响以及数据整理、处理过程产生的误差。海洋磁力测量准确度以主、联络测线相交点的测量值差的均方根值作为衡量依据,各种误差最终影响均方根的值,在不同精度的海上磁力测量中对这些误差进行比较分析,建立了海洋磁力测量误差数据处理模型,发现误差影响的分配规律变化。 The errors affecting the marine magnetic force measurement can be divided into real-time meas-urement errors and non-real-time measurement errors. The real-time measurement errors include measurement errors generated by GPS and magnetometer systems. Non-real-time measurement errors include errors in the implementation of marine magnetic measurements, dragging fish positioning errors, and normal field correction errors, as well as intrinsic daily changes in geomagnetic field, ship magnetic corrections, influences of dragging fish position reduction errors, and data collation and error generated by the process. Marine magnetic measurement accuracy is based on the root mean square value of the measured value difference between the intersections of the main and contact survey lines. Errors ultimately affect the RMS value. These errors are compared in marine magnetic surveys of different precisions. Based on the analysis, a data pro-cessing model of ocean magnetic measurement error was established, and the distribution rule of the influence of error was found.
杨子群1,张学贤2
1山东科技大学地球科学与工程学院,山东 青岛
2海洋沉积与环境地质国家海洋局重点实验室,国家海洋局第一海洋研究所,山东 青岛
收稿日期:2018年5月9日;录用日期:2018年5月30日;发布日期:2018年6月6日
影响海洋磁力测量的误差可以分为实时测量误差和非实时测量误差,实时测量误差包括GPS、磁力仪系统产生的测量误差。非实时测量误差包括海洋磁力测量实施过程中的误差、拖鱼定位误差和正常场校正误差,还包括固有的地磁日变改正、船磁改正、拖鱼位置归算误差的影响以及数据整理、处理过程产生的误差。海洋磁力测量准确度以主、联络测线相交点的测量值差的均方根值作为衡量依据,各种误差最终影响均方根的值,在不同精度的海上磁力测量中对这些误差进行比较分析,建立了海洋磁力测量误差数据处理模型,发现误差影响的分配规律变化。
关键词 :海洋磁力测量,误差数据处理模型,误差分配规律
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
影响海洋磁力测量的误差来源较多,是各种因素综合的影响 [
海洋磁力测量中,实时测量误差主要指由仪器系统产生的误差,包括GPS、磁力仪系统产生的测量误差,这类误差一般成系统性。
海洋磁力测量中,海洋磁力测量定位数据的采集由船载GPS定位系统完成,GPS的选择应该尽量满足海洋磁力测量的精度要求,并应该充分考虑测区的地理条件和经济合理性。现阶段的海洋磁力测量定位系统大多采用差分GPS设备,其定位精度可达到厘米级,可以满足各种精度要求的磁力测量,很大程度上消除了定位系统产生的误差。
现代海洋磁力仪按工作原理可以分为质子旋进式、欧弗豪塞(Overhauser)式和光泵式磁力梯度仪等3种不同类型 [
针对现今世界市场上海洋磁力仪产品型号较多、技术指标复杂、应用范围不尽相同,产品的选用存在一定困难等问题,标准质子旋进式海洋磁力仪灵敏度较低,存在进向误差,但无死区,价格最为低廉,适合于对灵敏度要求不高的工程和科研地球物理调查。Overhauser海洋磁力仪的灵敏度高,无进向误差,无死区,价格便宜,适合于大多数工程和科研地球物理调查。光泵式海洋磁力仪灵敏度和采样率最高,梯度容忍度最大,但存在死区和进向误差问题,适用于高精度的海洋磁力梯度调查和航空磁力调查。标准质子旋进式海洋磁力仪由于其灵敏度较低,在国内海洋工程和科研地球物理调查中已经较少使用。总体来说,Overhauser海洋磁力仪和光泵式海洋磁力仪各有优势,Overhauser磁力仪一般小巧轻便,价格适宜,而光泵磁力仪采样速度更高一些,梯度容忍度更大一些,因此具体选用应以具体情况而定来减小误差。
测量仪器的误差还存在于晶体振荡器的调试、传感器配谐的调试、静态下仪器信噪比测定与要求、仪器系统的稳定性和仪器系统的稳定性中。进行稳定性测定,若频率漂移,应调试或者更换。测定传感器配谐和选频放大器的每一档中心频率,不满足技术指标时,应更换软件。仪器调试后测量其信号电平值与噪声电平值,信噪比应不小于50,负责要继续调试。调试后的仪器系统进行持续工作状态试验,试验持续时间2d~3d,仪器系统应保持稳定状态 [
在海洋磁力测量中,为了得到测点的磁异常,还会有各种非实时测量误差,包括由仪器产生的非实时测量误差、拖鱼位置误差、正常场校正误差、船磁校正误差、日变校正误差、数据处理整理和算法因素所产生的误差和海洋磁力测量实施过程中的误差。
位置误差即磁力仪观测位置并非实际位置,而产生的误差。海洋磁力测量存在明显的动态性,受到风、流、涌等因素的影响测船定位。从而产生定位误差,本文将其称为位置误差的直接影响。而磁力测量中,为了减小船磁效应,海洋磁力仪通常采取拖曳式作业。这样,使得定位中心和磁测中心不一致,因此,海洋磁力测量数据处理时还必须对拖鱼的位置进行归算,本文将拖鱼位置归算时造成的误差称为位置误差的间接影响。
为了减小位置误差的影响,在施测过程中。应严格控制船速,使测船和拖鱼沿计划航线匀速、直线航行,当实际航线偏离计划航线较大时应该及时修正,以减小位置误差的影响。对于高精度的海洋磁力测量,还可以通过船姿传感器来藏小位置误差的直接影响。为了对拖鱼位置进行准确归算。在数据处理时,可先对测船的航迹线进行复原,即进行测线拟合。清除一些随机干扰的影响。按目前海上调查实施现状,如没有安装探头位置定位器来确定拖曳探头的具体位置,则通过调查船航向、船舷向和航速等记录数据来进行磁测探头的位置改正应是一种较为有效的解决方案。实际海洋磁测中,往往导航定位天线位置与磁测电缆固定位置不一致,这样就需要先利用简单改正推算电缆固定位置后,再结合航向、船舷向及航速等参数计算探头实际位置。同时由于海流的影响,船舷向在测量过程中是不断变化的,但其变化幅度一般不大,这样为简化计算并便于编程实现,可以考虑整条测线采用夹角 的平均值计算,对于较长的测线,或船舷向变化较大的情况,可以考虑分段平均的方法求取不同时段的夹角 平均值计算。合理有效地解决海洋磁测探头位置推算问题,将会提高磁测精度,同时也有利于进一步的磁测资料处理与解释工作。
为了得到测点的磁异常,必须对磁测资料进行正常场校正。正常场校正值是由模拟地磁场时空分布与时间变化的正常场模型提供的。正常场模型有模拟全球地磁场分布的球谐模型,也有模拟区域地磁场分布的多项式模型、球冠模型和矩谐模型 [
当选去IGRF模型作为正常场模型时,国际参考场模型在研究全球地磁场分布时效果较好,各国绘制的不同地区的地磁异常图便于比较。但是由于数据和计算能力的限制,它的分辨能力是有限的,不适宜处理某一地区地磁场或描述空间尺度较小的磁异常差异 [
地磁场是变化的磁场,地磁场随时闻的变化可分为长期变化和短期变化。由于海洋磁力测量的作业时间一般较短,因此影响海洋磁力测量的主要是短期变化,而短期变化影响中起主导作用的是太阳日变化。其变化量有时甚至可达40~60 nT [
由于海水层对地磁日变的衰减作用,在一定水深范围内,地磁日变随着水深的增加而减小,但并非线性关系,在实际中应避免使用深水海底日变数据作地磁日变改正。D. R. Auld等比较沿岸日变站与测区纬度相差不大的海底日变站同步数据发现,经度相差7˚,地磁日变化的差异不超过3 nT [
测量船一般都是用强磁材料建造,因此在船体周围,除地磁场之外,还存在一个来源于船体的附加磁场,简称船磁。船磁对海洋磁力测量的影响是随着测船的航向、测船与拖鱼的距离和所处的磁场面变化的。不同测船在不同方位、不同距离以及不同测区,船磁的影响是不同的。为了有效地减小船磁效应,海上磁力测量时,一般将融力传感器探头拖曳于测量船之后一定距离(《规范》要求海洋磁力仪拖缆长度不得小于3倍船长。在开始测量前后,必须在测区进行船磁试验。船磁试验包括,船磁方位影响的试验和拖缆拖曳距离确定试验,获得必要的船磁改正方位曲线,确定拖缆最佳拖曳距离。
为了减小船磁校正误差,船磁方位测量应进行2或3次。航向误差不得大于10度 [
海洋磁力测量船大都是强磁性材料建造,强磁材料磁性一旦形成很难消失,这就是船磁的固有磁部分。同是,随测量船所处的地磁场变化以及测量船相对地磁场的空间方位变化,船磁也在不断地变化,这部分瞬时而变的附加磁场称为船磁的感应磁部分,因此测量船的此话磁场包括固有磁和感应磁两部分。
同时,传统船磁模型的常数项或改正基准的选取问题一直没有解决好,使得不同航磁、区块的测线数据存在很大的系统差,文章《海洋磁力测量中船磁方位改正的基准确定》注意到了这个问题,提出取垂直于地磁场方位且两两方向相反的多次测量平均值作为改正基准。而这个问题的根本在于传统船磁八方位实验测量中无法消除试验点的地磁异常。
海洋磁力测量资料处理过程中,必然会造成误差的出现,原始数据中的异常数据剔除不完整,异常数据超过所有数据的10%,测线编辑,整理不符合《海洋磁力测量要求》的规定,测点的位置计算不正确,不满足《海洋磁力测量要求》定位精度的规定,正常场模型选择不合理。测区范围较大时,未进行多台站日变改正,或多站日变改正机制归算方法不合理,日变改正值计算方法不得当,结果不准确。船磁校正值计算方法不合理,结果不准确。磁测总精度评价不合理,参加评价的起点书大于总交叉点数的3%,结果不符合《海洋磁力测量要求》精度的规定。数据网格化方法不合理,等值线勾绘较杂乱,土建整饰及绘制精度不符合《海洋磁力测量要求》的规定。等值线间隔应大于磁测精度10倍;地磁异常剖面图中正、负异常的面积分配不均匀。这些都是海洋磁力测量后期数据整理、处理过程中出现误差的方面。
海上磁测方法不合理,定位和磁测数据不同步,定位和磁测采样间隔大于10 s,数据记录较正确、连续,但不够完整;缺失数据的测线里程不超过测线里程的1%,并且每处超过5 km;主测线、检查线方向和测线间隔符合《海洋磁力测量要求》要求。测量过程中的船速和航空控制不够好,航速变化大于0.2 km,测线偏航距离保持在1/5 km。施测中拖鱼状态稳定,因拖鱼不稳定造成的漏测区域超过测线总里程的1% 。未对重点异常进行加密探测等对海洋磁力测量都会产生误差影响。
影响海洋磁力测量的误差是各项误差的综合,根据误差传播定律,
各项误差总影响可由(1)式表示:
T i 2 = T 0 2 + T 1 2 + T 2 2 + T 3 2 + T 4 2 + T 5 2 + T 6 2 (1)
T海洋磁力测量误差总影响;T0为由仪器产生的实时测量误差;T1为拖鱼位置误差;T2为正常场校正误差;T3为船磁校正误差;T4为日变校正误差;M为均方根误差;T5为数据整理、处理过程产生的误差;T6为海洋磁力测量实施过程中的误差。
利用主测线与检查测线上的交叉点的磁测不符值对海洋磁力测量精度进行评定,其公式为:
M = ± ∑ i = 1 n ∑ j = 1 m Δ T i , j 2 N (2)
式中, Δ T i , j 为交差点磁测不符值,N为主测线与检查线实际交叉点的总数。海洋磁力测量中,大部分较大的强磁异常区的交叉点不符值大,可以不参加精度评估,但弃点数不能超过总交点数的3%。
在《海洋调查规范·海洋地质地球物理调查》中海洋地磁测量误差分配与比例尺要求相关,不同的调查比例尺的测量准确度要求如表1 [
海洋磁力测量准确度以主、联络测线相交点的测量值差的均方根值作为衡量依据(表2) [
M 2 ∝ T 2 = T 0 2 + T 1 2 + T 2 2 + T 3 2 + T 4 2 + T 5 2 + T 6 2 (3)
由于数据整理、处理过程产生的误差和海洋磁力测量实施过程中的误差不随海洋磁力测量比例尺的改变而有较大的改变,在求极限时可以视为常数α。
由上表可推出:
1) 当比例尺为 1 : 100 × 10 4 时:
lim T 0 → 2 ( T 0 ) 2 + lim T 1 → 2 ( T 1 ) 2 + lim T 2 → 2 ( T 2 ) 2 + lim T 3 → 2 ( T 3 ) 2 + lim T 4 → 2 ( T 4 ) 2 + α = T 2 ( T i ≥ 0 ) (4)
2) 当比例尺为 1 : 50 × 10 4 时:
lim T 0 → 1 ( T 0 ) 2 + lim T 1 → 2 ( T 1 ) 2 + lim T 2 → 1 ( T 2 ) 2 + lim T 3 → 1 ( T 3 ) 2 + lim T 4 → 2 ( T 4 ) 2 + α = T 2 ( T i ≥ 0 ) (5)
3) 当比例尺为 1 : 20 × 10 4 和 1 : 10 × 10 4 时:
lim T 0 → 1 ( T 0 ) 2 + lim T 1 → 1 ( T 1 ) 2 + lim T 2 → 0.5 ( T 2 ) 2 + lim T 3 → 1 ( T 3 ) 2 + lim T 4 → 1 ( T 4 ) 2 + α = T 2 ( T i ≥ 0 ) (6)
比例尺 | 仪器T0 | 拖鱼位置T1 | 正常场校正T2 | 船磁影响T3 | 日变校正T4 | 均方根差M |
---|---|---|---|---|---|---|
1:100 × 104 | ≤2 | ≤2 | ≤1 | ≤1 | ≤2 | ≤4 |
1:50 × 104 | ≤1 | ≤2 | ≤1 | ≤1 | ≤2 | ≤4 |
1:20 × 104 | ≤1 | ≤1 | ≤0.5 | ≤1 | ≤1 | ≤2 |
1:10 × 104 | ≤1 | ≤1 | ≤0.5 | ≤1 | ≤1 | ≤2 |
表1. 不同的调查比例尺的测量精度范围表
比例尺 | 仪器T0 | 拖鱼位置T1 | 正常场校正T2 | 船磁影响T3 | 日变校正T4 | 总误差 T i 2 |
---|---|---|---|---|---|---|
1:100 × 104 | ( ≤ 2 2 ) + ( ≤ 2 2 ) + ( ≤ 1 2 ) + ( ≤ 1 2 ) + ( ≤ 2 2 ) + T 5 2 + T 6 2 = ( ≤ T i 2 ) | |||||
1:50 × 104 | ( ≤ 1 2 ) + ( ≤ 2 2 ) + ( ≤ 1 2 ) + ( ≤ 1 2 ) + ( ≤ 2 2 ) + T 5 2 + T 6 2 = ( ≤ T i 2 ) | |||||
1:20 × 104 | ( ≤ 1 2 ) + ( ≤ 1 2 ) + ( ≤ 0.5 2 ) + ( ≤ 1 2 ) + ( ≤ 1 2 ) + T 5 2 + T 6 2 = ( ≤ T i 2 ) | |||||
1:10 × 104 | ( ≤ 1 2 ) + ( ≤ 1 2 ) + ( ≤ 0.5 2 ) + ( ≤ 1 2 ) + ( ≤ 1 2 ) + T 5 2 + T 6 2 = ( ≤ T i 2 ) |
表2. 不同调查比例尺的测量精度约束范围表
当Ti取极大值时,不同海洋磁力测量比例尺误差分配情况如图1。
误差全定义域范围内平均分配,可得出关系如图2。
当Ti取中值时,不同海洋磁力测量比例尺误差分配情况如图3。
可得出关系如图4。
Ti在取中值与极大值时具有相似的折线走势,说明在海洋磁力测量中误差的取值具有相似的误差分配特点,在实际的海洋磁力测量过程中,各个误差的大小在中值与极大值之间,比较Ti取中值与极大值得折线统计图,具有一定的说服性。
在实际的海洋磁力测量过程中,各个误差的大小在中值与极大值之间,比较Ti取中值与极大值得折线统计图,可以发现误差随着海洋磁力测量的精度具有一定的分配规律:
1) 不同的Ti取值却有相似的误差分布规律,在小比例尺 1 : 100 × 10 4 时,按照《规范》的准确度容忍度,实测仪器误差、拖鱼位置误差,日变校正误差的对总误差影响较大,正常场校正误差与船磁校正误
图1. Ti取极大值不同比例尺下不同误差范围所占总体误差饼状图。(a) Ti取极大值比例尺为1:100 × 104时,不同误差范围所占总体误差饼状图;(b) Ti取极大值比例尺为1:50 × 104时,不同误差范围所占总体误差饼状图;(c) Ti取极大值比例尺为1:20 × 104和1:10 × 104时,不同误差范围所占总体误差饼状图
图2. Ti取极大值时不同比例尺下误差的变化折线图
图3. Ti取中值时不同比例尺下不同误差范围所占总体误差饼状图。(a) Ti取中值比例尺为1:100 × 104时,不同误差范围所占总体误差饼状图;(b) Ti取极大值比例尺为1:50 × 104时,不同误差范围所占总体误差饼状图;(c) Ti取极大值比例尺为1:20 × 104和1:10 × 104时,不同误差范围所占总体误差饼状图
图4. Ti取中值时不同比例尺下误差的变化折线图
差对总误差影响较小;数据整理、处理过程产生的误差和海洋磁力测量实施过程中的误差介于两者中间,在 1 : 50 × 10 4 比例尺时,拖鱼位置误差与日变校正误差为主要的误差来源,实时仪器测量误差、正常场校正误差和船磁校正误差相似且相对较小,数据处理整理误差和算法因素所产生的误差和海洋磁力测量实施过程中的误差略有增加,已经是不可忽视的误差影响;在 1 : 20 × 10 4 和 1 : 10 × 10 4 比例尺时,正常场误差所占比例最小,实测误差、拖鱼位置误差、船磁校正误差、和日变校正误差均等且在不同Ti取值时所占比例相似,而数据整理、处理过程产生的误差和海洋磁力测量实施过程中的误差在不同Ti取值时,所占比例有较大变化,随着Ti的取值的增大数据整理、处理过程产生的误差和海洋磁力测量实施过程中的误差所占比例增大,但任然是最主要的误差来源。
2) 随着比例尺的增大,实时测量误差即仪器误差的影响所占比例先缩小后变大,在实际海洋磁力测量过程中,要根据比例尺的实际情况做出相应的调整;拖鱼位置误差和日变校正误差基本一致,随着比例尺的增大,先增大后减小,并且在 1 : 50 × 10 4 时所占比例很大,建议在 1 : 50 × 10 4 度时,可以适当提高拖鱼位置和日变校正的准确度减小两者的误差,从而对减小总误差有明显的作用;正常场校正误差所占比例变化不大且较小;船磁校正误差随着比例尺精度增大而增大,在实际测量中随着精度要求的增大,要削弱船磁校正误差的影响;数据整理、处理过程产生的误差和海洋磁力测量实施过程中的误差在小比例尺测量时,影响程度不大,但是在大比例尺,海洋磁力测量精确度高的作业中,数据整理、处理过程产生的误差和海洋磁力测量实施过程中的误差成为主要误差。
3) 海洋磁力测量的误差是多项因素的综合误差,因此,对于其对海洋磁力测量的影响应综合考虑。本文建立了海洋磁力测量的数据误差处理模型,对影响海洋磁力测量的各项误差进行了分析,得出各误差在不同比例尺的海洋磁力测量中误差的变化规律,给出了海洋磁力测量误差分配方案,在确定海洋磁力测量的比例尺精确度要求后,一方面应该在测量实施时对各种误差的影响进行控制,另一方面应对海洋磁力测量的影响综合考虑,根据各项误差分配方案采取措施着重消除突出误差部分。
杨子群,张学贤. 海洋磁力测量过程中的误差分配规律分析 Error Source and Error Compensation Method in Marine Magnetic Survey[J]. 海洋科学前沿, 2018, 05(02): 57-66. https://doi.org/10.12677/AMS.2018.52007