 点击上方“溪流之海洋人生”即可订阅哦ADCP是一种利用声学多普勒频移原理进行矢量流速测量的仪器。目前被广泛用于水运工程建设，水文环境监测，海洋、河口的流场结构调查，港口航道流速和流量测验等。目前国内外的ADCP主要在室内流速检定水槽和室外江河湖海进行流速计量或校验，两种计量测试环境各有利弊。室内流速检定水槽由于受限于建设尺寸，难以计量低频ADCP，难以开展分层流速计量，边界回响噪声影响测试结果,且水槽中的水过于纯净，需要人为添加一定浓度的散射粒子，如石灰、气泡等，添加粒子的方式及保证粒子的均匀分布是较难的课题；室外江河湖海则较为宽广，水体中含有天然均匀的散射粒子，但由于受风浪、紊流、潮汐、声速剖面等自然因素的影响，给测量结果引入诸多不确定度分量。本文依托船闸航道（长180m×宽25m×深10m），宽阔、水深，水面平静无风浪，水体有悬浮物质、可控分层流场等有利ADCP校准的试验环境和新型无人船轻便、高稳定性、超强搭载兼容性的优点，建立测量仪器与计量标准的溯源链，保障ADCP流速量值的准确可靠。
& p* D+ [( \$ ^0 v 一、计量测试系统设计
+ q9 U# p3 K3 Q# g* R2 l& ^ ⒈系统硬件设计
9 w+ m5 t4 n) |6 v 传统ADCP检测方法采用拖车作为载体，整体装置较为复杂且拖车速度受限。以无人船作为搭载试验仪器的载体，其速度能够达到最高3m/s，并且无人船体积较小、更加灵活自由，可按规定航线航行，较为理想。本文采用定制的自动无人船，船身采用高分子聚酯碳纤维材料减少磁性金属对检测设备的影响，并搭载绝对直线技术。校准系统通过转接法兰盘将被测ADCP换能器安装到无人船底部，且将ADCP通讯电缆接入船体电子水密舱。通过法兰盘转接的方式可以适配不同型号的ADCP。将GPS-RTK安装在无人船顶部，并尽量保证二者的轴线在同一直线上，如图1所示。以确保被测ADCP与GPS-RTK所测量的速度为船身相同位置速度值。  # x1 _8 V+ ]4 X# M4 v
图1 测量船示意图 / M4 L+ t/ S9 _: _7 D6 T
⒉试验方法 O% `# f. R) }* W8 A& a9 S" c
通过转接法兰盘将被测ADCP换能器安装到无人船底部，并将GPS-RTK安装在无人船顶部，且保证二者在轴线重合，采样频率均为1Hz。采用全站型电子速测仪（简称全站仪）放样两条180m等长且平行的虚拟直线，一条直线与水池长边平行且通过水池短边中点，以此作为设定航迹线。另一条直线位于近岸，作为标准线便于计量。将船闸分为15m+80m+15m的三段，前15m作为无人船加速段，中间80m作为匀速测量段，后15m作为减速段。确定航迹线起始点与终止点坐标信息，并将其输入到无人船平台系统控制软件进行航线指令编辑。 5 s0 C3 ?6 T) B. U# b
在测量段始末点分别设置安装高速光感应摄像与运动分析系统，作用是采集无人船进入测量段和离开测量段瞬间的运动影像，并逐帧分析无人船特征点进/出测量段的时刻。使用全站仪对测量段水平距离进行准确测量，试验方法示意图如图2所示。  $ g3 H5 Z% b1 @; L4 A' L
图2 试验方法示意图 8 }( v! ~! {' b- u: q: k8 [/ K
无人船根据事先编辑任务指令进入航线并在水上自主航行，同时被测ADCP采集流速数据。根据测量段水平距离和无人船在测量段的运行时间，计算平均速度作为流速参考标准值，与被测ADCP流速测量值比对，计算流速示值误差。 ! w- \! n1 R! O# S
同时收集无人船匀速运动时GPS-RTK测量数据。GPS-RTK所测原始数据为坐标数据，可以利用公式(1)使其转换为船的实时速度值VG，与被测ADCP测得船速进行对比，经计算可以核查无人船的瞬时速度稳定性。 & x- ]- }$ t* V/ I2 a. T
VG＝|(Δx2＋Δy2)-2|/t ⑴ 1 l, r) [- S4 t0 ]
式中，Δx为GPS-RTK在x轴方向相邻两点坐标差，m；Δy为在y轴方向相邻两点坐标差，m；t为两点间测量时间间隔，s。
: h2 {. f/ U4 v. J 二、试验结果 5 D' D" T; Z8 s* f
全站仪所测流速测量段长度L＝79.771m。无人船在测量段内运行时间t＝78s。无人船进入测量段时刻（触发预警线）为15:10:12，无人船离开测量段时刻（触发预警线）为15:11:30。计算平均流速参考标准值V0＝1.02m/s。 . e' Q8 U; v2 W8 \
同一时间内被测ADCP所测数据如图3所示。  . J& e# j/ t# p8 h+ G! L
图3 被测ADCP所测数据 - E) I. K1 J1 _( B2 i
试验过程中，GPS-RTK测得结果与被测ADCP测得流速数据的对比如图4所示。  e& s! {# H: _! Z
（a)  + T: m, X* g7 J" M4 l( B
(b)  fill=%23FFFFFF%3E%3Crect x=249 y=126 width=1 height=1%3E%3C/rect%3E%3C/g%3E%3C/g%3E%3C/svg%3E)
8 X O0 Y1 g/ R; [4 s' D (c)
# M( r& C" z2 C3 M, J. @ 图4 流速数据的对比 4 g) `0 g4 M$ ~+ a
在匀速测量段内记录120组流速数据，测速结果及被测ADCP示值误差如表1所示，表中仅列出间隔10秒的数据。 $ S" k4 I5 I/ t* I: E' ^3 n
( F( w! F2 s6 Y* b 根据HY/T 102-2007《声学多普勒流速剖面仪检测方法》中所描述的数据处理方法，首先利用准则剔除数据中的粗大误差，剩余101组数据，大于100组，符合要求。计算GPS-RTK测量序列的平均值被测ADCP的测量序列的平均值0.982m/s。 ' ?. D" o5 X! M, P$ X+ F; C ?
剔除粗大误差后的101组数据对比示意图如图5所示： % \( N) z* Z3 ^+ g8 _* f6 j/ F% b
图5 101组数据对比图 3 g- Q, U! U& B: q2 X( }& O% I7 K
三、测量结果不确定度分析 5 A4 B( }3 `; W% t$ J x
⒈数学模型及灵敏度系数（略） 1 h6 T6 O/ v2 c; \5 t& \- X9 \. o
⒉不确定度分量计算（略）
6 G$ j; p) A) @9 S ⒊合成不确定度 4 T8 ^! a1 k8 n4 ]8 P3 S- H5 @& h7 p. \
在本次试验过程中，流速测量误差数学模型中的L、t、V及ΔV0各个分量互不相关。根据相关公式可以计算此计量测试系统的合成标准不确定度u(ΔV)＝0.0355m/s，取包含因子k＝2，则测试结果测量不确定度为： . ?2 [0 ]& F2 U/ n& C
u＝K×U(ΔV)＝0.072m/s，k＝2 # R$ O; _# z; v7 r* q
四、结束语
$ {7 ^& H. j7 r. B5 K3 s 声学多普勒流速剖面仪具有直接测出断面流速剖面、不扰动流场、测验历时短、测速范围大等特点，因此被广泛应用于各处港口航道的流速、流量检测，海域水体结构分析，以及海域大范围搜救工作等。本文以无人船为载体设计的ADCP计量测试系统具有体积小、操控灵便的优点；以天津港船闸作为校准场所，综合了室内及室外校准场所的优点，并尽可能规避其干扰。在无人船设定航速1m/s时，被检ADCP流速示值误差-0.029 m/s，评定本测试系统的扩展不确定度U＝0.071m/s，k＝2。
; f4 Z8 i0 L# G/ u2 A, @- B 1 - t* ~ ~5 `' j' S2 o
END : _' h3 ?+ T. T# K! v" H: Z
1 ! G, j% x s& G# i1 G$ ?
【作者简介】文/谢慧 柳义成 张明敏，分别来自天津大学精密测试技术及仪器国家重点实验室和交通运输部天津水运工程科学研究院。第一作者谢慧，女，1995年出生，硕士在读，主要从事海洋声学仪器及软件开发研究；通讯作者柳义成，男，1990年出生，工程师，硕士，主要从事海洋测绘及计量检测研究。本文为基金项目，国家质量基础的共性技术研究与应用重点专项（2018YFF0212202）、中央级公益性科学研究所基本科研业务费专项（TKS160102，TKS170201）。文章来自《仪表技术与传感器》（2020年第10期），参考文献略，编发时有所删节，用于学习与交流，版权归作者及出版社共同拥有，作者已授权转载。 2 ~1 I6 i4 z/ {! n9 H, g
相关阅读推荐 声学技术▏声学多普勒测速仪标校技术研究微课堂▏邓锴：国产声学多普勒流速剖面仪(ADCP) 研发历程及现状海洋论坛▏国产走航式声学海流剖面仪的应用分析海洋技术▏深海船载走航式声学多普勒流速剖面仪海洋技术▏声学深拖系统在海底冷泉调查中的应用海洋论坛▏我国海洋水文调查设备的发展历程 # O! |0 m% ]- y7 D0 m* E7 s/ X
公众号 2 r& C2 D _1 F V$ j% T7 ^4 l
溪流之海洋人生
& M8 E5 |' }5 s! s* @& m: k/ v- ^ 微信号▏xiliu92899 ; u' x# H7 v& ]9 H# Q
用专业精神创造价值
1 @, m: R* S, c& r6 u 用人文关怀引发共鸣 您的关注就是我们前行的动力
. e8 F$ C/ }9 z1 R 投稿邮箱▏12163440@qq.com
) ?& ^8 W- ?8 n% d4 }" X2 }6 h3 j' I% ~ M! _
+ B* U: k+ v; f, |. O( D7 f7 n- p. r$ z0 c+ d
# a5 i/ `0 W9 q0 u. `6 E, n) r | 
