ADCP河流流量测量原理和方法

n = (实测换能器离河底高度 × λ– 盲区 - Δ) / 单元长度 Dc
(10)
参数λ取决于“旁瓣”区的厚度。如果“旁瓣”区的厚度大约为水深的6%，λ=0.94。对于 窄带型ADCP,Δ=0。对于宽带型ADCP，Δ由下式确定：
5
∆
=
1 2
(Dp
+
D0
−
Dc )
(11)
式中： Ｄp ＝脉冲长度 Ｄ0 ＝脉冲延迟长度
RDI公司开发出它的第三代ADCP：“骏马”系列ADCP。第一代和第二代ADCP尺寸和重量都很
大。价格昂贵。第三代ADCP尺寸和重量都大大减小。价格降低了一半左右。第三代ADCP的诞
生，标志着ADCP技术已达到成熟和普及的阶段。至2001年底，大约有3500千多台RDI公司生
产的ADCP应用于世界各地的海洋和河流。
本文简要介绍ADCP和ADCP河流流量测量原理和方法。
2.0 ADCP简介
第一台商品ADCP诞生于七十年代中期。它是由声学多普勒计程仪演变出来的，因此是船载式
ADCP。1982年，美国RDI公司开发出它的第一代ADCP：窄带ADCP；1991年，RDI公司开发出它
的第二代ADCP：宽带ADCP
（关于窄带和宽带ADCP的解释可参见附录2）；1995年，
顾名思义，ADCP是一种利用声学多普勒原理测量水流速度剖面的仪器。与机械式或电磁式流
1
速仪不同，ADCP是一种遥测仪器。即ADCP流速测量点不在ADCP所处的位置。而且，ADCP可以 测量距其一定范围内的许多点的流速。也就是说，一台ADCP可以替代许多台机械式或电磁式 流速仪。
为了测量三维流速，ADCP一般装备有四个（或三个）声换能器。换能器总是安装成与ADCP轴 线成一定倾角。每个换能器即是发射器又是接收器。换能器发射的声波尽可能集中于较窄的 声束范围内。每一个换能器对应于一个声束。换能器发射某一固定频率的声波，然后聆听被 水体中颗粒物散射回来的声波。假定水体中颗粒物与水体流速相同。当颗粒物的移动方向是 接近换能器时，换能器聆听到的回波频率比发射波频率高。当颗粒物的移动方向是背离换能 器时，换能器聆听到的回波频率比发射波频率低。声学多普勒频移，即发射声波频率与回波 频率之差由下式确定：
有两种方法可以测量船速。第一种方法称为“底跟踪”。“底跟踪”是指ADCP通过接收和处 理由河底或海底的回波信号跟踪河底或海底的运动。如果河底或海底没有移动的悬移质，“ 底跟踪”所测量的速度即是船速。第二种方法是利用GPS测量船速。由航迹上任意两点的GPS 座标值可以得到两点间的位移，再除以相对应的时间即得到船速。一般来说GPS方法测船速 不如“底跟踪”方法精确。
迹，流量等。实时测量结果还可以列表显示。图
2
2
为ADCP系统软件屏幕显示。在进行流量测量作业中，ADCP实时测出水体相对于作业船的速度 和作业船相对于河底的速度 (即船速)。水体的真实流速则由相对速度减去船速(矢量差)来得到。ADCP同时还测出水深( 类似于声纳测深)。这些数据(包括流速,船速,水深)由电脑在系统操作软件控制下实时采集 ，处理，并实时计算每一微断面的流量，当作业船沿某断面从河一侧驶至另一侧时,即给出 河流流量。
ADCP通过跟踪水体中颗粒物的运动（称为“水跟踪”）所测量的速度是水流相对于ADCP（也 即ADCP安装平台）的速度。当ADCP安装在固定平台上，“水跟踪”测量的流速即为水流的绝 对速度。当ADCP安装在船上（移动平台），“水跟踪”测量的相对速度中扣除船速（平台的 移动速度）后才得到水流的绝对速度。
图 1 用ADCP进行河流流量测量
图 2 ADCP系统软件屏幕显示 (WinRiver) 3.0 ADCP流量测量原理 从理论上讲，ADCP流量测量原理与传统的人工船测，桥测，缆道测量，和涉水测量的基本原 理是一样的，都是在测流断面上布设多条垂线。在每条垂线处测量水深并测量多点的流速从 而得到垂线平均流速。只是ADCP所测的垂线可以很多，每条垂线上的测点也很多。然而，与 传统方法相比，ADCP方法有如下几点不同： （1） 传统方法是静态方法。即无论是人工船测，桥测，缆道测量，还是涉水测量，仪器总
中部流量由下式给出:
mn
m
∑ ∑ ∑ QM =
f j Dc∆t = [(VxMVby − VyMVbx )]i ⋅ (Z 2 − Z1 )i ∆t
(16)
式中: VxT=表层平均流速(x向分量) VxM=中部平均流速(x向分量) VxB=底层平均流速(x向分量)
中部平均流速和流量
中部平均流速由ADCP直接测出,其值为所有有效单元所测流速之平均。x方向分量由下式算出 (y方向分量类似):
∑ VxM
=
1 n
n
uxj
j=1
(17)
式中: uxj =单元 j 中所测的x向流速分量
2002年4月修改稿
ADCP河流流量测量原理和方法
黄河宁
美国亚迪仪器（RDI）公司
1.0 引言 2.0 ADCP简介 3.0 ADCP流量测量原理 4.0 岸边流量估算 5.0 ADCP流量测量基本要求 6.0 ADCP流量测量中的几个问题 7.0 参考文献
附录A ADCP名词解释 附录B ADCP分类 附录C ADCP流量测量精度分析 附录D “骏马”系列“瑞江”牌河流ADCP简介
1.0 引言
河流流量测量是水文工作者的重要任务之一。传统的河流流量测量方法包括人工船测，桥测 ，缆道测量，和涉水测量等。其基本原理是在测流断面上布设多条垂线。在每条垂线处测量 水深并用流速仪测量一至几个点的流速从而得到垂线平均流速。进而得到断面面积和断面平 均流速。流量则由断面面积和断面平均流速的乘积得到。这种传统方法费工费时，效率低。
ADCP基于如下的公式计算流量:
Q = ∫∫u •ξ ds
(２)
S
式中: Q=流量 S=河流某断面面积 u=河流断面某点处流速矢量 ξ=作业船航迹上的单位法线矢量 ds=河流断面上微元面积
ds由下式确定:
ds = Vb ⋅ dz ⋅ dt
(3)
式中: dz=垂向微元长度
dt=时间微元
Vb=船速矢量 |Vb|=作业船速度(沿航迹) z=垂向坐标,z=0为河底,z=H为水面(H=水深)
Fd
= 2Fs
V C
(1)
式中：Fd
是声学多普勒频移；Fs
是回波频率；V是颗粒物沿声束方向的移动速度（即沿声束方向的水流速度）；C是声波在水
中的传播速度。
ADCP每个换能器轴线即为一个声束座标。每个换能器测量的流速是水流沿其声束座标方向的 速度。任意三个换能器轴线即组成一组相互独立的空间声束座标系。另外，ADCP自身定义有 直角座标系：X-YZ。Z方向与ADCP轴线方向一致。ADCP首先测出沿每一声束座标的流速。然后经过座标转换为 X-Y-Z座标系下的三维流速。然而X-YZ是局部座标系。利用罗盘和倾斜计提供的方向和倾斜数据，X-YZ座标下的流速可转换为地球座标系下的三维流速。
|Vb | 由下式确定:
4
Vb = Vbx 2 + Vby 2
(4)
式中: Vbx 和 Vby 分别为x和y方向的船速。 (x, y) 为水平面上的总体座标系。 公式(２)可以重新写为:
∫ ∫ ∫ ∫ Q
=
T
H
u
⋅
dz
•ξ
| Vb
| ⋅dt
=
T
H
(u ×Vb ) • k ⋅ dz dt
ADCP用于河流流量测量已有十年的历史了。RDI公司的第一代，第二代，和第三代ADCP都有
源自文库
专门为河流流量测量设计和制造的河流型ADCP。图1显示了“骏马”系列“瑞江”牌河流型A
DCP图片。整个测量系统包括三个主要部分:(1)ADCP,（2)电脑，(3)数据采集软件。软件主
要功能包括:实时数据采集和显示，以及流量计算。实时显示的数据有流速矢量，水深，航
设Df为第一个有效单元(靠近换能器)的深度，则：
Df = 换能器深度 + 仪器盲区 + 单元长度 Dc + Δ
(12)
设Dl为最后一个有效单元(靠近河底且未受河底干扰的单元)的深度，则：
Dl = Df + n Dc
(13)
设Z1为河底至最后一个有效单元的高度,Z2为河底至第一个有效单元的高度,则:
1 Z1 = H − Dl − 2 Dc
(14)
1 Z2 = H − D f + 2 Dc
(15)
深度平均流速的x方向分量由下式算出(y方向分量类似):
∫ ∫ ∫ ∫ Vx
=
1 H
H
ux
0
⋅ dz
=
1 H
Z1
Z2
H
[ uxdz + uxdz + uxdz] =
0
Z1
Z2
1 H
[Z1VxB
+ (Z2
− Z1 )VxM
+ (H
− Z2 )VxT ]
是固定于所测垂线上进行测量。ADCP方法是动态方法。ADCP在随测量船运动过程中进 行测量。 （2） 传统方法的测流断面通常要求垂直于河岸。ADCP方法不要求测流断面垂直于河岸。船
3
航行的轨迹可以是斜线或曲线。
为了计算流量，ADCP在走航测量中测量如下数据：
（1） 水的相对速度（相对于船的速度，由“水跟踪”测出） （2） 船速(由“底跟踪”测出,或由GPS算出) （3） 水深（由河底回波强度测出，类似于回声测深仪） （4） 船的航行轨迹（由船速和计时数据算出，或由GPS算出）
用ADCP进行河流流量测量是近十年才发展和应用的新的流量测量方法。ADCP，是英文Acoust ic Doppler Current Profiler的缩写。中文通常译为声学多普勒流速剖面仪。ADCP流量测量方法的发明被认为是 河流流量测量领域的一次革命。当装备有ADCP的测量船从河流某断面一侧航行至另一侧时， ADCP即刻测出河流流量。它比传统的河流流量测量方法提高效率几十倍。它标志着河流流量 测量的现代化。
前面已提到，从相对速度中扣除船速得到水流的绝对速度。
对于某一测量断面，ADCP将测出该断面中部区域内的流速。这个区域称为ADCP实测区。然而 ，在四个边缘小区域内ADCP不能提供测量数据或有效测量数据。第一个区域靠近水面，其厚 度大约为ADCP换能器入水深度，ADCP盲区,以及单元尺寸一半之和。第二个区域靠近河底， 称为“旁瓣”区（河底对声束的干扰区）。其厚度取决于ADCP换能器轴线倾角（与ADCP轴线 的夹角）。例如RDI公司“瑞江”牌ADCP换能器轴线倾角为20度，相应的“旁瓣”区厚度大 约为水深的6%。第三和第四个区域为靠近河岸两侧的区域，其水深较浅，测量船不能靠近或 者不能保证至少一个或两个有效测量单元。这四个小区域称为非实测区。其流速需通过实测 区数据外延来估算。
3.1 流量计算方法
ADCP流量计算方法最早由Christensen 和 Herrick (1982), 以及 Simpson and Oltmann (1993)提出。 他们采用单元流速与船速矢量叉乘积（称为 f 值）进行计算公式的推导。 这里，我们同时采用矢量叉乘积和微断面的深度平均流速矢量来进行公式推导。
ADCP基本上有两种安装方式和应用，一种是安装在固定平台上（如河底或海底），进行定点 垂线（或水平）流速分层测量（在海洋水文中称为流速剖面测量）。例如RDI公司生产的“ 骏马”系列“哨兵”型ADCP的应用。另一种是安装在活动或移动平台上（如调查船上），进 行走航流速剖面测量。例如RDI公司生产的“骏马”系列“瑞江”牌ADCP的应用。
y可以是任意固定坐标。
将沿航迹的断面离散为 m 个微断面,则：
∑ ∫ ∑ ∑ Q = m [ Hi i=1 0
m
m
fi ⋅ dz] ∆t = [(V ×Vb ) • k ]i ⋅ Hi∆t = [VxVby − VyVbx ]i ⋅ Hi∆t
i=1
i=1
(7)
式中:
Hi=微断面 i 处水深 ∆t=相应于微断面的测量时间平均步长 m=断面内总的微断面数目
V=相应于每一微断面的深度平均流速矢量:
每一微断面的深度平均流速矢量V由下式计算:
∫ V
=
1
H
u ⋅ dz
H0
(8)
式中 H由下式确定:
H = 换能器入水深度 + ADCP实测深度
(９)
ADCP在进行测量时,不但在水平方向上将沿航迹的断面划分为微断面(共m个),同时在垂向将 水体划分为小单元。有效单元的数目n由下式确定:
(5)
0 0

00
式中: T=跨断面航行时间 k=垂向单位矢量
定义 f 为点流速流速与船速矢量叉乘积 (Simpson and Oltmann, 1993):
f = (u ×Vb ) • k = uxVby − uyVbx
(6)
式中角标x和y分别表示x和y方向的速度矢量。注意到叉乘积的结果与坐标无关,即:x-