Comparison and Analysis on Conductivity-Temperature-Depth System (CTD) Data Quality
物理海洋学的研究离不开大量的实测海洋数据, 海洋调查是获取实时海洋数据最直接且最为有效的方法, 集成了各种性能优异的传感器仪器设备是海洋调查的基础。海洋的温盐深观测是现今海洋调查的基本内容之一, 从20世纪60年代开始温盐深测量仪(CTD)已经被广泛应用于海洋调查中(陈上及和马继瑞, 1991; 张爱军, 1992), 它的发展随着科学技术的进步而同步前进。CTD能适应如走航实时观测、定点自容观测、抛弃式探头观测等多种观测方式, 并且获取调查海区水体现场实时的高精度物理环境参数, 是目前物理海洋调查中使用最为广泛的仪器设备之一。传感器技术是调查设备的基础, 其各方面性能是衡量仪器设备好坏的关键, 同时也是调查数据质量的保证, 各种数据订正方案应运而生(Fofonoff et al., 1974; Giles and McDougall, 1986; Grose, 1981; 任强等, 2014), 但是在长期的观测中, 传感器的稳定性、漂移、准确度等指标依然是最重要的部分。目前, 世界上已有许多设备生产厂家推出了高精度的CTD, 但是各个厂家采用不同类型的传感器, 因而传感器性能等方面会存在一定的差异性。本文将对3种不同型号的CTD所测数据进行对比分析。
1 传感器性能及资料来源
中国科学院海洋研究所于2014年6月在黄海中部布放坐底潜标, 其上搭载了3种不同厂家生产的温盐深测量仪, 型号分别为美国Sea-Bird公司的SBE 37 Coastal(以下简称37Coastal)、德国Sea & Sun公司的CTD 48, 以及意大利Idronaut公司的304Plus, 3种设备适用于近海应用并且都能够添加溶解氧等附加传感器。表 1为3种CTD的各个传感器参数, 从表 1中可知, 各CTD的压力、温度和电导率传感器的精度和分辨率差别较小, 其传感器指标能满足近海海洋的长期观测, 并且能够获得高精度和高分辨率的压力、温度及盐度等数据。该潜标布放时间为4个月左右, 3种设备在潜标的位置如图 1所示。本文主要针对3种设备所获取的温盐深数据质量进行对比分析, 讨论传感器精度及稳定性等指标参数。
表 1 传感器技术指标表
Tab. 1 Technical index of sensor
型号 |
压力 |
温度 |
电导率 |
范围(m) |
分辨率FS(%) |
精度FS(%) |
范围(℃) |
分辨率(℃) |
精度(℃) |
范围(mS/cm) |
分辨率(mS/cm) |
精度(mS/cm) |
CTD48 |
0~500 |
0.002 |
±0.01 |
–2~36 |
0.0006 |
±0.001 |
0~70 |
0.001 |
±0.001 |
304Plus |
0~500 |
0.0015 |
±0.05 |
–5~35 |
0.0001 |
±0.002 |
0~90 |
0.0003 |
±0.003 |
37Coastal |
0~350 |
0.002 |
±0.1 |
–5~45 |
0.0001 |
±0.002 |
0~70 |
0.0001 |
±0.003 |
2 数据对比分析
2.1 压力数据对比
图 2为3种传感器获取的压力数据, 蓝色曲线代表 37Coastal, 红色代表304Plus, 黑色代表CTD48。从图 2中可以看出, 三者所获得的压力呈周期波动, 其波动周期均为15d, 与大小潮周期一致。37Coastal的压力曲线与CTD48基本重合, 而304Plus的压力曲线与其他二者存在微小偏差。进而, 对三者压力数据进行两两相减并求各自方差(图 3), 结果显示, CTD48与37Coastal的压力差最小, 平均压力差值为0.02m左右, 两者方差为0.0796m; 304Plus比其他两者压力差稍微偏大, 偏差值约为0.2m。由于海底地形复杂, 受地势及海底沉积物的影响, 潜标在海底的真实状态可能是倾斜的, 所以会造成安装设备出现绝对高度差, 导致304Plus测量的压力与其他两者存在偏差。
3种CTD的采样间隔均为半小时, 在海洋中某个点上的要素值在半小时内的正常变化值非常微小, 其对应的传感器获取的数据波动范围也应该很小, 因此能够利用数据梯度值(两两相减)判断获取的数据稳定程度, 同时数据的稳定程度能够反映传感器的稳定性能。根据三者的压力数据求压力梯度▽P的绝对值(图 4), 平均绝对值小则表示数据波动小且数据稳定。计算结果如下: 37Coastal的平均绝对值▽P和▽P标准差分别为0.008m和0.004m; 304Plus的平均绝对值▽P和▽P标准差分别为0.1275m和0.1562m; CTD48的平均绝对值▽P和▽P标准差分别为0.1232m和0.1492m。数值结果表明, 37Coastal平均绝对值▽P和▽P标准差均为三者中最小, 而CTD48和304Plus的数值结果几乎相等。同时, 三者的压力方差分析P值均大于0.05, 说明三者之间不存在显著的差异。
图 4压力梯度▽P与图 3的压力差dP显示, 在相同的时间节点上出现了数据振幅异常增大的现象。而后对2014年6~10月西太平洋台风过境该潜标点的时间进行了统计(表 2), 从时间上可以看出, 压力变化剧烈期与台风过境时间相一致, 说明台风的过境造成水体大幅扰动, 短时间内急剧变化的水体扰动对传感器有一定影响, 使自身测量受到一定的干扰。并且, 不同的传感器对剧烈变化的环境参数响应时间不一致, 使同一时间测量到的压力值不相等, 因此在水团环境剧烈变化时期内, 不同传感器压力测量值有一定偏差。从以上压力数据分析表明, 三种压力传感器性能方面都表现良好, 在不考虑安装角度及潜标姿态的情况下对于本次实验来说CTD48与37Coastal压力数据相差最小, 如果把潜标的姿态造成的高度差进行修正过后, 其三者的压力数据都几乎相等, 因此在数据准确度方面基本相同, 且37Coastal压力传感器在稳定性方面具有优势。
表 2 台风过境时间表
Tab. 2 Transit time of typhoon
台风 |
浣熊 |
麦德姆 |
夏浪 |
娜基莉 |
台风编号 |
201408 |
201410 |
201411 |
201412 |
过境时间 |
2014-07-09 |
2014-07-25 |
2014-08-08 |
2014-08-02 |
JulianDay (d) |
189 |
205 |
218 |
212 |
2.2 温度数据对比
图 5为不同CTD采集的温度数据, 可以看出三者基本一致, 任意两者的相关系数均在0.99以上。对温度进行两两相减(图 6)可以看出, 37Coastal和CTD48的温差最小且数据波动小, 平均温差为–0.005℃, 其标准差为0.0014℃。而其他两组温度差均呈现出一个变大的趋势, 最大温度差为0.055℃, 且标准差大于第一组(37Coastal和CTD48)。三者方差分析结果均为P > 0.05, 说明温度之间不存在显著性差异。
同处理压力数据类似, 对3组温度数据分别求时间序列上的温度梯度▽T (图 7)。结果显示, 三者▽T波动情况一致, 3组数据任意两者相关系数在0.98以上, 说明三者温度测量值偏差非常小。37Coastal的平均绝对值▽T和▽T标准差分别为0.0050℃和0.0084℃, CTD48的平均绝对值▽T和▽T标准差分别为0.0050℃和0.0085℃, 304Plus的平均绝对值▽T和▽T标准差分别为0.0059℃和0.086℃, 数值显示三者温度传感器所获取的数据波动小, 其传感器稳定性能几乎一致。
2.3 盐度数据对比
在长时间序列的测量过程中, 传感器本身的异常电信号干扰等会产生少量的错误异常值, 应先删除部分异常错误数据, 图 8为去掉异常数据后的盐度曲线。数据显示, 三者在160(开始测量时间)~265d(图 8中盐度突变时间点)的数据吻合程度较好, 而在第265d以后37Coastal的数据与其他两者的数据出现了较大偏差, 其他两组数据出现突变且高度一致趋势。
第293d为该潜标回收时间点, 且在潜标回收后对该点进行了SBE 9型CTD的投放, 该CTD是国内外应用非常广泛且具有高数据质量的温盐深测量仪器。在该时间点上37Coastal、304Plus与CTD48的盐度值(PSU)分别为32.9、31.8、31.8, 三者盐度最大相差1.1。图 9为SBE9型CTD所获取的该点盐度数据, 与潜标同一深度的盐度为32.9, 与37Coastal测得的盐度相同。因此, 可以判断304Plus和CTD48盐度数据异常, 同时能够说明304Plus和CTD48在第265d以后所测得的盐度异常。3种CTD基本在同一深度位置, 304Plus和CTD48盐度值同时出现异常突变, 考虑到不同设备的传感器探头同时出现同样故障的概率极低, 导致该现象出现的可能是二者的电导率传感器同时被某生物附着, 电导率测量管堵塞, 测量值出现异常。
排除生物附着等外力因素影响, 对第265d以后的盐度数据进行截断, 取160~265d盐度数据(三者相对稳定区间)进行对比(图 10)。图 10中显示, 三者数据变化趋势整体一致, 其中37Coastal与304Plus数据的相关系数较其他两组数据的相关系数高, 达到了0.97, 37Coastal与CTD48数据的相关系数为0.9。对3组数据求相互盐度差值(图 11), 可以看出第一组盐度差(37Coastal与304Plus)波动最小, 其平均值和标准差分别为0.016和0.03, 均小于其他两组数据。
同处理压力和温度数据类似, 分别对3组盐度数据求时间序列上的盐度梯度▽S (图 12)。计算出37Coastal的平均绝对值▽S和▽S标准差分别为0.003和0.0047, CTD48的平均绝对值▽S和▽S标准差分别为0.003和0.0067, 304Plus的平均绝对值▽S和▽S标准差分别为0.020和0.029。说明37Coastal和CTD48在数据稳定性方面一致, 但是从前文可知, 37Coastal和304Plus数据准确度更高。三种设备中只有37Coastal带有泵结构, 通过泵的作用使三种传感器能对同一微水团进行测量, 从最大程度上保证了数据的准确性, 对于复杂海况条件下及剖面式投放的数据测量具有优势。
3 结语
本文对搭载在坐底潜标平台上的3种不同型号的自容式温盐深测量设备(37Coastal、CTD48和304Plus)所获取的长时间序列的压力、温度和盐度数据进行对比分析。分析结果表明, 三者传感器的数据准确度及稳定性等方面存在微小差异。37Coastal专门被设计用于复杂的近海海况, 因而其在变化剧烈的近海环境下所获取的温度和盐度数据质量均有优异的表现, 而CTD48和304Plus的主要适应对象为深海探测, 在以后的工作中将会通过其在大洋中的数据对比来比对各自传感器的性能。目前, 从集成度、操作性及经济性等方面考虑, CTD48和304Plus两种设备要优于37Coastal。所以, 3种设备对海洋真实的物理环境参数变化观测在其准确度及稳定性等方面都有良好的保障, 对不同的探测环境需要认真选择不同的探测传感器, 以发挥仪器本身最大的应用价值。对于如何选择不同种类的温盐深传感器要结合自身的科研需求以及使用环境等, 根据需求而后详细评估目前各个型号的温盐深传感器的精度、准确度、量程及漂移性等方面指标, 并且在布放的成本以及难度系数等方面也需要认真考虑, 综合各个方面的因素确定最终适合自身科学研究需求的方案。