SAG-2M型与KSS31M型海洋重力仪比测结果分析

2020-07-26张菲菲孙建伟韩波杜润林王万银

物探与化探 2020年4期

张菲菲，孙建伟，韩波，杜润林，王万银

(1.自然资源部天然气水合物重点实验室青岛海洋地质研究所，山东青岛 266071； 2.海洋国家实验室海洋矿产资源评价与探测技术功能实验室，山东青岛 266071； 3.长安大学重磁方法技术研究所，陕西西安 710054；4.长安大学地质工程与测绘学院，陕西西安 710054；5.长安大学西部矿产资源与地质工程教育部重点实验室，陕西西安 710054)

0 引言

海洋重力测量是海洋地球物理测量的主要技术手段之一，重力数据在现代国防领域、深海资源勘探领域以及地球科学研究领域中有广泛的应用[1-5]，是国家的重要战略资源。受研究目标的需求，获得精度更高、分辨率更好的重力场数据是未来海洋重力测量的一项重要任务，而海洋重力数据的精度依赖于采集设备的稳定性能和测量精度。目前，国内外海洋重力调查中使用的主流海洋重力测量设备主要是德国的KSS31型海洋重力仪、美国的L&R型海洋重力仪、美国的BGM型海洋重力仪以及俄罗斯与加拿大合作生产的GT型海洋重力仪[1-2,6-7]，仪器测量精度均可达到1×10-5m/s2。国内海洋重力仪研制起步比较晚，最早是从20世纪60年代开始开展海洋重力仪的研制工作，但均未形成成熟的商业化产品。近年来，随着科技的发展，国内针对海洋重力仪进行了积极的研制工作，其中比较有代表性的是中国航天科技集团公司第九研究院第十三研究所(航天十三所)自主研发的海洋重力仪SAG-2M型海洋重力仪[7]，该重力仪为捷联惯导式重力仪，测试内符合精度为1×10-5m/s2，并实现了小批量生产。

2018年青岛海洋地质研究所拟引进一套SAG-2M重力仪，根据现行的海洋重力测量规范[8-9]，新重力仪投入使用前需明确测量仪器的指标要求，对技术指标进行验证和评估是测量作业规程中非常重要的组成部分。海洋重力仪的技术性能指标一般包括测量范围(量程)、抗干扰能力(动态范围)、动态重复测量精度、零点漂移特性、采样率和工作温度等内容[10-11]，其中动态重复测量精度是海洋重力仪最关键的技术指标。动态重复测量精度可以通过实验室动态重复观测和海洋实际工作两种方式来进行检验[10]，海洋实际工作的测量精度更能体现测量环境时空变化对重力测量成果的影响。海洋作业验证新重力仪的测量精度可以用外部符合方法对重力仪的实际测量精度进行检验，即采用两台或多台重力仪同平台观测，通过对比不同仪器在同一时刻的观测结果获取参试重力仪测量精度的评估信息，该方法被国内外多家海洋调查和科研机构所广泛使用，并取得了较好的比测结果。欧阳永忠等[12]对GT-1A、TAGS、L&R SⅡ、SGA-WZ01和GDP-1五套航空重力仪进行了同机测试，对所获取的重复线和测网成果数据进行了对比分析，GT-1A型航空重力仪综合性能技术指标最佳，国产SGA-WZ01捷联航空重力仪具有最佳的重复线测量精度，此项工作掌握了多类型商业重力仪的技术性能指标，验证国产重力仪数据可靠性；张向宇等[13]对GT-2M、KSS31M和ZLS三种不同类型的重力仪进行了同船比测工作，测试结果表明，相同的测量条件下，3种型号重力仪的测量结果一致性较好，GT-2M型重力仪数据精度略好；张振波[14]对GT-1M和KSS31M海洋重力仪进行了比对，以验证GT-1M重力仪的可靠性，测试结果表明2台重力仪测量结果吻合，但是该工作中KSS31M重力数据使用的是网格数据，而不是同船比测数据。国外也有过GT与美国L&R重力仪同船比对测量的实例[15]，但是针对SAG-2M型海洋重力仪的同船测试结果较少，研发单位曾开展过与GT-1A型重力仪的测试比对，比对精度小于1×10-5m/s2，但未见详细测试结果。

为系统了解并检测SAG-2M型海洋重力仪的技术性能和工作特点，青岛海洋地质研究所于2018年8月搭载SAG-2M和KSS31M重力仪进行了同船比测工作，并对所获取的重力数据进行了数据处理及对比分析。这2种型号的海洋重力仪传感器和工作原理有所不同，但仪器的技术性能指标相近，通过同船比测及数据对比分析可以深入了解不同型号海洋重力仪的工作特点，验证新重力仪数据可靠性，为开展相关技术研究与未来装备选型提供重要参考依据。

1 SAG-2M和KSS31M海洋重力仪参数对比

1.1 SAG-2M海洋重力仪基本参数

SAG-2M型海洋重力仪是航天十三所自主研发的海洋重力仪，该重力仪采用高精度石英加速度计式重力传感器，将重力传感器集成于捷联惯性/GNSS组合系统中，通过捷联惯性/GNSS组合系统数学解算计算出重力传感器敏感轴指向，然后经过坐标分解得出重力传感器测量值的垂向分量。

SAG-2M型海洋重力仪系统主要的测量部件为重力仪主机(图1)，主机内部集成了石英加速度计式重力传感器、高精度光纤陀螺仪、二次电源转换电路、计算机电路、重力传感器高精度信号转换电路、GNSS OEM板卡、高精度温控电路等。重力仪主要配套部件包括显控记录装置、UPS电源单元和减震支架。该海洋重力仪结构简单，体积小巧，与传统海洋重力仪相比无CC效应，传感器动态范围宽，对载体平台动态环境适应性好，主要技术参数见表1。

图1 SAG-2M型海洋重力仪主机Fig.1 The main frame of SAG-2 Mmarine gravimeter

1.2 KSS31M型海洋重力仪基本参数

KSS31M型海洋重力仪是德国Bodensee公司生产的海洋重力仪[16-19]，其平台抗干扰能力强，仪器测量精度高，被国内外多家单位和科研机构使用，工作性能得到业内认可。该重力仪采用直立弹簧式重力传感器，重力采集系统主要由两部分组成(图2)，分别是重力探头及陀螺平台稳定系统(KT31)和数据采集控制系统(DHS)，其主要技术参数见表1所示。该型号海洋重力仪抗干扰能力强，仪器测量精度高，通常小于1×10-5m/s2。

表1 SAG-2M与KSS31M海洋重力仪技术参数对比

2 数据采集及处理

2.1 数据采集

2018年8月青岛海洋地质研究所组织SAG-2M及KSS31M型海洋重力仪同船搭载“海大号”综合科考船执行调查采集任务，共采集测线44条，其中主测线25条，联络测线19条，测线布设见图3。数据采集过程严格按照仪器操作和海洋地质调查规范执行，测量期间仪器运行正常，在同一泊位进行重力基点比对测量，月漂移量均小于2.0×10-5m/s2。

图3 测线布设(底图为测区地形)Fig.3 Distribution map of survey lines (residual topographic map)

2.2 数据处理

本次数据处理依据新版海洋地质调查规范[1-2,9-10,20]，利用CGG公司LCT重磁处理软件进行了数据处理，处理步骤及参数如下。

1) 数据读取。SAG-2M和KSS31M重力仪采样间隔为 1 s，对原始采集数据进行了整理，按照时间提取了各测线相关数据项：年、儒列日、时、分、秒、重力仪读数。为了进行位置校准，从同船导航数据中提取了经度、纬度数据。

2) 延迟校正。SAG-2M和KSS31M海洋重力仪测量的重力读数与实时的GPS定位数据之间存在滤波延迟，数据处理时需要根据仪器所选档位的滤波延迟时间对重力读数进行延迟校正。本次测量过程中SAG-2M型海洋重力仪滤波延迟时间为189 s，KSS31M型海洋重力仪滤波延迟时间为76 s。

3) 零点漂移校正。重力仪零点漂移均按照线性变化处理，根据出航及返航在港口码头重力基点处重力读数，计算零漂量，按照时间分配到重力仪读数上。零点漂移计算公式为[1-2]：

(1)

其中：t1为航次起始时基点读数时间；t2为航次结束时基点读数时间；ti为采样点观测时读数时间；g零漂为零点漂移量，单位10-5m/s2；δR为零点漂移改正值，单位10-5m/s2。

4) 厄特渥斯校正。厄特渥斯效应是走航式重力测量过程中，科氏力对于安装在测量船上的重力仪所施加的影响，导致测量重力值与实际重力值不符，该效应与航向、航速和船只所在的地理纬度有关。为了消除厄特渥斯效应，利用定位数据计算了船只的航速、航向，对所有测线数据进行了厄特渥斯改正，厄特渥斯改正公式为[1-2]：

δge=7.499×V×sinA·cosφ+0.004V2，

(2)

其中，δge为厄特渥斯改正值，单位10-5m/s2；V为航速，单位m/h；A为航向角；φ为测点的地理纬度。

5) 正常场校正。采用CGCS2000正常重力公式计算正常重力场，其计算公式为[1-2,21]：

(3)

其中：γ0为正常重力场值，单位10-5m/s2，φ为测点的地理纬度。

6) 自由空间重力异常计算。测点自由空间重力异常计算公式为：

g=g0+CΔs+δR+δge-γ0，

(4)

其中：g为测点的自由空间重力异常值，单位10-5m/s2；g0为基点绝对重力值，单位10-5m/s2；C为重力仪格值常数；Δs为测点与基点之间的重力仪读数差，单位10-5m/s2。

3 数据比对分析

为方便讨论测量精度问题，笔者使用了3个衡量测量精度的指标参数[11,22-23]，分别为平均绝对偏差(mad)、均方根(rms)以及标准偏差(std)，使用相关系数(rab)来衡量2组数据之间的相关性，计算公式分别为：

(5)

(6)

(7)

(8)

为详细评估2台重力仪所获得数据的测量精度，依据上述指标参数，对2台重力仪自由空间重力数据的交点差、同测线数据、全区网格数据进行了对比分析。

3.1 数据交点差

利用LCT重磁处理软件计算了2台仪器自身及仪器间的交点差，仪器自身是计算主测线与联络测线交点差，仪器间是计算主测线与另一台仪器的联络测线交点差，表2给出了2台仪器自身和仪器之间交点差平均绝对偏差、均方根及标准偏差的统计结果。主测线与联络测线之间交点共有480个，交点差统计结果可以看出2台仪器在本次的测量中精度较高，优于1×10-5m/s2，两者的自身交点差相近，表明2台仪器处理后数据精度相近；2台仪器之间交点差相近，这说明2台仪器处理后数据不存在系统差值。对比所有的交点差结果，SAG-2M相关的交点差相对小，均小于KSS31M自身的交点差，从交点差统计的角度来看，在本次同船比测中，SAG-2M型海洋重力仪处理后数据精度略优于KSS31M型海洋重力仪。

表2 仪器交点差统计 10-5 m/s2

3.2 测线数据比对

选取不同测量位置的8条主测线进行了对比分析，将8条测线处理后自由空间重力异常以平面剖面图展示(图4)，由图4可知，2台重力仪得到的自由空间重力异常总体趋势一致，异常特征高度相似。但2台重力仪在此次测量中也各自存在问题，KSS31M重力仪在Z01、Z07、Z28测线存在跳点数据(图5中绿色圆圈所示)，而SAG-2M重力仪不存在跳点数据；SAG-2M重力仪在Z46、Z55测线存在锯齿状数据(图5中紫色矩形框所示)，即数据抖动的情况，而KSS31M重力仪不存在该情况。

图5 网格数据偏差分布Fig.5 The distribution map of grid data deviation

图4 部分测线自由空间重力异常平面剖面示意Fig.4 The free air gravity anomaly profile map of parts of survey lines

为了分析两者之间的偏差，计算了2台重力仪在同一测线的自由空间重力数据的相关系数，两者偏差的平均绝对偏差、均方根以及标准偏差，计算结果见表3。由2台仪器处理后同测线数据偏差统计结果可得，SAG-2M与KSS31M2台重力仪在8条测线上的测量结果高度线性相关，相关系数均>0.99；各测线之间的偏差相当，平均绝对偏差范围(0.59～1.028)×10-5m/s2；2台仪器在同一测线上的标准偏差较小，均<0.8×10-5m/s2，同测线偏差统计结果略大于仪器间交点差统计结果。通过对测线自由空间重力异常曲线及数据偏差比对，可知SAG-2M重力仪与KSS31M重力仪在本次测量中数据处理结果吻合，这也表明SAG-2M在此次测量中得到的结果是可靠的。

表3 同测线2台重力仪测量偏差统计

3.3 网格数据比对

为了更全面地对比2台重力仪测量数据之间的关系，对全区网格化数据进行了对比分析。将2台重力仪得到的自由空间重力数据进行了网格化，网格间距为2 km，搜索半径为20 km，共得到网格节点8 976个。

对2台重力仪全区网格化数据进行了偏差统计(表4)，两者之间相关系数0.994 5，为高度线性相关，网格数据偏差的绝对平均偏差、均方根以及标准偏差与测线比对结果相近。图5为网格数据偏差，两套数据间偏差多分布于(-2～2)×10-5m/s2之间，约占总网格节点的97.3%。对数据偏差超过3倍标准偏差的点进行了统计(图5紫色圆点)，偏差较大的网格点共有248个点，占总网格节点的 2.76%。选取偏差较大网格点分布相对较多的测线Z01和Z40，对两台仪器的自由空间重力异常曲线进行对比(图6)，Z01线偏差较大点位于测线左侧，两台仪器的自由空间重力异常趋势相似，但是KSS31M数据有跳动，造成了两者之间偏差较大；Z40线偏差较大点位于测线中央位置，两台仪器的自由空间重力异常趋势高度相似，但是幅值略有偏差。分析其原因，由于两台仪器的测量原理有所不同，在实际测量中会因为海况或仪器本身信号处理方式不同造成数据之间有所偏差，但不影响整体数据质量。