海洋重力仪是以舰船为载体,使用重力传感器测定海面重力加速度的地球物理勘测仪器。重力加速度的准确测量,对国防、计量、测绘、地质、地震与资源勘探等领域均有十分重要的意义。随着机械制造工艺日益完善,广大海洋研究人员对海洋重力仪测量精度的要求越来越高。目前在全球海洋科考船上安装应用较为广泛的海洋重力仪型号有:美国Micro-gLaCoste公司生产的MGS6、SeaⅡ,德国Bodenseewerk公司生产的KSS,俄罗斯GravimetricTechnologies公司生产的GT-2M等。与此同时,国内一些科研机构在海洋重力仪研制方面也取得了重大进展,如:中国科学院测量与地球物理研究所研制的CH-Ⅱ型弹簧重力仪、中国船舶重工集团公司第707研究所研制的ZLⅡ型海洋重力仪、中国航天科技集团第十三研究所研制的SAG-2M捷联式重力仪等均日趋成熟,并逐步推向市场。 国内外许多学者针对不同类型海洋重力仪数据处理流程及技术性能等进行了详细的分析研究,如:2018年Yuan等利用“向阳红06”号科考船搭载俄罗斯GT、美国SeaⅡ、中国SAG、ZLⅡ、CH-Ⅱ、SAG多套不同型号的海洋重力仪,进行了多条重复线测量,并对测量数据精度进行对比;Kovrizhnykh等在Caspian海域同船搭载La Coste & RombergZLS、ChekanAM和GT-2M,对上述三套设备测量精度及处理方法进行了比对;张向宇等对GT-2M型海洋重力仪记录的原始数据进行了分析,总结了GT-2M型海洋重力仪数据处理流程;张涛等利用相关分析法对LaCosteSeaⅡ海洋重力仪数据重新进行了交叉耦合分析与改正计算;张向宇等对ZLⅡ型海洋重力仪记录的原始数据进行预处理分析,提出了该型重力仪新的数据处理方法。 2020年中国新建的极地科考破冰船——“雪龙2”号引进了Micro-gLaCoste公司最新生产的SeaⅢ型海洋重力仪,该套设备也是SeaⅢ型在中国的首套应用。“雪龙2”号是无限航区航行破冰科考船,而SeaⅢ型海洋重力仪测量量程大于等于20000×10−5m·s2,满足其全球测量的使用需求。作为该型重力仪的国内首套应用用户,为了测试该型仪器的工作状态、检查其静态工作性能(静态月漂移)及动态精度,分别于2019年12月和2020年6月组织对其进行了实验室内静态测试及海上动态精度测试,为后续其他用户提供相关参考依据。
一、SeaⅢ型海洋重力仪简介 ⒈SeaⅢ型海洋重力仪结构组成及工作原理 Micro-gLaCosteSeaⅢ型海洋重力仪(外形见图1)由工作终端、常平架、重力仪主机、减震阻尼、UPS电源及控制系统等组成,系统结构紧凑、布局合理,便于设备操作与维护。SeaⅢ型海洋重力仪采用第三代高可靠性固态光纤陀螺、固态加速度计,及高度集成的数据控制系统,以提升仪器的可靠性与稳定性。
图1 “雪龙2”号SeaⅢ型海洋重力仪安装图 SeaⅢ型海洋重力仪设计上沿用了SeaⅡ型摆杆-斜拉零长弹簧原理。传感器主体为一个由零长弹簧支撑的铰链摆(图2),摆前部设有上下两个空气阻尼器,空气阻尼器对摆的垂直运动产生较强的阻尼作用;摆后部为由摆和两块固定于传感器内的金属板构成的一组电容器。当重力变化或有垂直干扰加速度作用在摆上时,摆的摆动引起电容器中电容量的变化,电容量的变化率直接反映出摆的摆动速率。电容变化信号通过电容位置指示器转换为直流电压信号输出,经模数转换和系统软件处理即可计算出重力值。
图2 SeaⅢ型海洋重力原理图 ⒉eaⅢ型海洋重力仪性能特点 相比目前国内应用较多的上一代SeaⅡ型海洋重力仪,SeaⅢ型海洋重力仪传感器和常平架体积更小,系统集成度更高;系统重量更轻,约降低30%,易于拆装转移,为不同科考船间的相互共享提供了良好条件;采用的常平架滑环技术和传感器恒温恒压控制,使得系统平台更加稳定可靠;常平架可机械锁止,不需单独拆卸,仪器拆卸和养护更为便捷。表1为SeaⅢ型海洋重力仪精度指标统计表。 表1 SeaⅢ型海洋重力仪精度指标统计表
⒊SeaⅢ型海洋重力仪重力值计算 SeaⅢ型海洋重力仪是通过测定摆在强空气阻尼中的摆动速度来表示重力的微小变化,通过弹簧张力和摆的运动速率换算重力读数,SeaⅢ型海洋重力仪具有在运动状态中实时感应重力值的特性。SeaⅢ型海洋重力仪测量基本方程式为: G=aST+Eotvos+TC⑴ 式中,G为重力读数;a为弹簧张力和重力的转换系数;ST为弹簧张力;Eotvos的值由公式(2)计算;TC为总改正,其计算方式见公式(3)。 Eotvos=7.503×v×cosλ×cosα+0.004154v2⑵ 式中,v为船速(节);λ为纬度(度);α为航向(度)。 TC=CC+BVC+XG+LG⑶ 式中,CC为交叉耦合改正,即仪器在垂向加速度和水平加速度共同作用下引起的附加重力扰动,由公式(4)计算;BVC为摆速改正,由公式(5)计算获得;XG为横向陀螺改正;LG为纵向陀螺改正。
式中,βVCC、βVE、βAL、βAX、βAX2为CC改正系数,仪器出厂时厂家标定后给出;VCC(y〞z′)为纵向加速度与垂向速度的耦合;VE(z〞)2为垂向加速度耦合的平方;AL(y〞z〞)为纵向加速度和垂向加速度的耦合;AX(x〞z〞)为横向加速度和垂向加速度的耦合;AX2[(x〞)2z〞]为横向加速度的平方和垂向加速度的耦合。
式中,dB/dt(i)为摆时间导数的有限差分,由公式(6)计算;B为摆位,从重力数据记录文件中获得;p0、p1、p2、q1、q2仪器自身标定系数,参数值从配置文件中获得。
二、SeaⅢ型海洋重力仪静态数据分析 ⒈月漂移 重力仪月漂移指标是仪器稳定性能的重要评判指标。2019年12月1日至2019年12月30日,在北京某实验室内对SeaⅢ型海洋重力仪进行了为期一个月的静态试验。根据SeaⅢ型海洋重力仪一个月的实验室内静态观测数据计算重力仪读数随时间变化的线性回归直线方程为:y=3.29×10−7x+512(x的单位为秒,y的单位为10−5m·s2),得出该台重力仪的静态掉格呈线性(图3,横坐标单位为秒,纵坐标单位为10−5m·s2),月漂移为0.85×10−5m·s2,优于国家海洋调查规范中对海洋重力仪月漂移不超过3×10−5m·s2的要求。表明“雪龙2”号配备的SeaⅢ型海洋重力仪静态月漂移较小,抗干扰能力强。
图3 SeaⅢ型海洋重力仪一个月静态数据变化 ⒉固体潮影响 固体潮是指在日月引潮力的影响下,引起地球固体部分岩石圈的形态变化,从而引起重力周期性变化的现象。固体潮值由于具有周期性,可以通过模型计算。选取SeaⅢ型海洋重力仪连续三天观测的数据(图4,横坐标单位为秒,纵坐标单位为10−5m·s2),其中锯齿状曲线为重力仪静态数据观测值,光滑曲线为固体潮模型计算值,横坐标为时间,纵坐标(左)为重力仪读数,纵坐标(右)为固体潮模型计算值。如图4所示,重力仪测量值呈现明显周期性变化,低值与固体潮低值相对应,局部高值也相互对应,曲线形态较为一致。说明重力测量值读数变化的周期性是受当地固体潮影响所致,表明“雪龙2”号配备的SeaⅢ型海洋重力仪具有较高的分辨率。
图4 SeaⅢ型海洋重力仪三天静态数据变化 连续3天的重力仪观测读数值显示,观测值叠加了振幅为0.2×10−5m·s2的扰动(图4)。严格意义的实验室内静态测试需要重力仪底座与基岩直接固定连接[10],但国内仅少数专业实验室能够满足上述条件。在普通实验室内深挖数十米到达基岩,并建立基岩与重力仪底座的直接连接是不现实的,所以图中重力仪观测值叠加了振幅高达0.2×10−5m·s2的扰动可能与平台受地基震动影响有关,但也不排除电子信号线路可能存在高频白噪声。 三、“雪龙2”号SeaⅢ型海洋重力仪动态精度测量 海洋重力仪动态精度分为交叉点精度和重复线精度。其测量一般通过布设纵横交错的网格状测线,在主、检测线相交点处进行重复测量,并对测量数据进行符合度解算,由此评估海洋重力仪自身重复性精度和测量作业精度,即为交叉点精度。以重复测线不符值、测线网交叉点不符值为标准,对海洋重力测量的内符合精度进行分析评估。受海况、定位精度影响,科考船海上航行测量完全的重复线较为困难,交叉测线更容易实现,所以一般采用以交叉点精度为主、重复线精度为辅的评价方法进行仪器动态精度衡量。 2020年6月15日至2020年6月28日,“雪龙2”号在中国南海进行了SeaⅢ型海洋重力仪的海上试验,重点对SeaⅢ型海洋重力仪的动态精度进行衡量。试验海区布设东西方向和南北方向各6条测线(图5),每条测线长10海里,测线间隔0.8海里,其中3、9号测线包括方向相反的重复线。试验海区水深介于450~560m之间,海底地形由西北至东南缓慢变深。
图5 SeaⅢ型海洋重力仪南海测网 ⒈交叉点精度 交叉点内符合方程式为:
式中,a为重力观测读数内符合中误差;d为主、检测线交差点重力不符值或重复线观测重力值互差值;n为测线交叉点或重复点个数。统计计算交叉测线交叉点误差值(未做平差),交叉点中误差为0.230×10–5m·s2,见下表2。 表2交叉点误差值(单位:×10−5m·s2)
统计交差点精度分布如图6所示,交叉点误差呈正态分布状分布,误差趋向于0,表明仪器动态性能优良。
图6 交叉点精度分布图 ⒉重复线精度 由于海洋重力仪采样过程中采样数据在空间上是点状离散的分布,任何两条重复测线采样点在空间上都不一定完全重合,为了评价测线重复精度,要求对测线相同位置进行重新插值采样,计算测线相同位置处的重力差异。 重复线精度判定内符合中误差计算公式为:
式中,θRMS为重力观测读数内符合中误差,δ为相邻点重力不符值或重复线观测重力值互差值,n为测线交叉点或重复点个数。 南北向重复线及误差值曲线如图7所示,东西向重复线及误差值曲线如图8所示。
图7 南北向重复线对比
图8 东西向重复线对比 对SeaⅢ型海洋重力仪南北向、东西向重复线进行插值采样计算其差值、平均值、中误差结果如表3所示。 表3SeaⅢ型海洋重力仪重复线分析结果(单位:×10−5m·s2) 东西向、南北向重力仪重复线对比分析显示,重复线最大误差值均分布在测线边缘两端,最小误差值位于测线中间位置处。表明重力仪在刚上线处还未趋于完全相对稳定造成误差值偏大,而当位于测线中间位置重力仪稳定后误差值相对较小,后续任务中为了获取高精度测量值,可适当延长测线。重力仪重复线整体中误差值为0.27×10−5m·s2,表明SeaⅢ型海洋重力仪具有较高的动态精度。 ⒊SeaⅢ型海洋重力仪动态精度分析 ⑴“雪龙2”号船船舶稳定性 “雪龙2”号船长122.5m,宽19m,排水量1.4万吨,尺寸和排水量与国内其他科考船相比均有较大优势,在海洋中航行也更为平稳。“雪龙2”号专门设置了一个重力实验室用于海洋重力观测,该实验室位于船舶稳心位置,有利于走航工况下的重力观测。同时,重力仪在该实验室内的安装位置与“雪龙2”号的中轴线基本重合,重力传感器纵轴与船的艏艉线方向一致,减小了船横摇、纵摇对重力仪观测的影响。且重力实验室舱室面积较小,仅供重力观测使用,无人员频繁进出的扰动,便于建立和保持满足重力仪使用的恒温、恒压的环境条件。“雪龙2”号采用吊舱系统推进,船体内振动噪音较小。SeaⅢ型海洋重力仪避开了船舱机器震动源,避免了船体振动噪声对重力仪测量的影响。 ⑵SeaⅢ型海洋重力仪器重力参数滤波窗函数 重力仪数据采集时,为了消除波浪等垂直加速度的影响,采集软件对弹簧张力ST、总改正TC及Eotovs改正都进行了300s低通滤波。滤波器是利用Exact Blankman窗函数设计的低通滤波器,比起Blankman窗函数,Exact Blankman窗主瓣宽度窄、旁瓣能级底、主瓣能量更加集中,滤波窗函数采用余弦窗裁断信号,可以获取较高的计算精度,公式为: w(n)=a0-a1cos(Ø)+a2cos(2Ø)⑼ 式中,w(n)为窗函数序列,n=0,1,2,…,N,N为滤波长度;a0=7938/18608;a1=9240/18608;a2=1430/18608;Ø=2πn/(N-1)。 四、小结与展望 综上,Micro-gLa Coste SeaⅢ型海洋重力仪零点漂移小,具备较好的静态线性漂移率,月漂移约为0.85×10−5m·s2,满足海洋调查规范要求。该型重力仪分辨率高,重力仪观测值能够清晰反映当地固体潮变化特征。SeaⅢ型海洋重力仪在“雪龙2”号测得的动态精度远优于国家海洋调查规范的要求,在近年来新配的海洋重力仪中精度也处于较优水平,统计计算海上测网交叉点不符值,动态内符合精度为0.23×10−5m·s2。海上动态精度优良,测量可信度高,优良的动态精度得益于“雪龙2”号船舶的稳定性、振动噪声低及重力仪自身采集软件较高的计算精度滤波算法。 目前SeaⅢ型海洋重力仪已经搭载“雪龙2”号上参与相关极地航次调查,该仪器运行稳定,抗干扰能力强。后续将结合极地现场应用情况,进一步分析探讨该型重力仪极区使用性能。 - Z9 `5 _$ M" s, B! i) a6 e
|