供水管道泄漏检测及相关仪的原理与使用
（南通市自来水公司 徐少童）
摘 要 介绍了相关仪的基本原理，使用方法等
关键词 相关 数字滤波 噪声
引言
随着我国的经济建设的发展，水资源短缺越来越成为限制我们发展的瓶颈之一，如何解决这个问题已经被逐步提到了战略高度，因此，合理利用水资源，降低漏损就成了我们水利工作者的重中之重。
减少漏损就要有相应的方法，目前我国大部份地区的检漏手段还停留在几十年前的水平，而国外在近二三十年则有了很大的发展，我们要做好这项工作就必须了解他们的技术，并能够最终掌握。
当前，简陋技术最先进的设备当属相关仪了，国外已有普通相关仪，多探头相关记录仪等多种产品，但究其根本，原理都是一样的，本人经过多方学习以及查阅相关资料，对其原理有了进一步的认识，下面就先从相关仪的基本原理说起。
一． 相关仪的基本原理
当管道发生泄漏时，能够产生比普通水声频率高较多的声压波沿管道传播，泄漏噪声频率高低主要取决于泄漏点的大小，泄漏噪声传播速度主要取决于管道直径和管材；通过放置在管道两端（泄漏点包围在中间）的振动传感器或声发射传感器测量泄漏信号，由于泄漏点可能位于管道不同位置，因此泄漏声传播到达两个传感器的时间不同，利用两列信号的互相关分析，一般即可确定泄漏噪声到达两个传感器的时间差。根据该时间差，通过两个传感器间的距离和声波在该管材中的传播速度，即可计算出泄漏点距两个传感器的距离。
设)(),(t y t x 为所测量的两列信号，则其相关函数计算公式如下：
)()()(1
lim )(0τττ-=-=⎰∞→yx T
T xy R dt t x t y T R
若信号为周期信号或一段信号可以反映信号全部特征，则可以采用一个共同周期或一段信号内的均值代替整个历程的平均值。对于泄漏声波信号，只要采集的两列信号均覆盖了在500m 以内泄漏声传播的全过程即可，不必无限制采集。这样，互相关函数计算公式可如下近似：
)()()(1
)(max 0max τττ-=+=⎰yx T xy R dt t y t x T R
相关系数计算公式为：
y x y x xy y x T y
x T xy R dt t y t x T σσμμτσσμμττρ-=-+=⎰∞→)()()(1
lim )(0
其中：τ为延迟时间，y x μμ,分别为)(),(t y t x 两列信号的均值，y x σσ,分别
为)(),(t y t x 两列信号的方差，1)(≤τρxy
。 由于相关系数在[-1，1]范围内，显然，采用相关系数计算比较方便，容易判断。通过互相关系数的计算，位于最大相关系数的τ值，即为泄漏到达两个传感器的时间差。
由于事先无法知道泄漏点距哪一个传感器近，所以两个传感器的两列信号中泄漏噪声到达时间或早或晚于另一个，为了准确确定延迟时间，每次必须计算)()(τρτρ-xy xy 和，由上面可知：
y x y x yx y x y
x xy xy R R σσμμτσσμμττρ-=--=-)()()( 为了区分两个振动传感器，将其标定为红色传感器和蓝色传感器，假定)(t x 为红色传感器测的信号，)(t y 为兰色传感器测的信号，以红色传感器测定的信号)(t x 为标准，若)(τρxy 波形内找到相关系数最大的点，则表明)(t x 信号中泄漏到
达的时间早于)(t y ，即泄漏距红色传感器近些；若)(τρ-xy 波形内找到相关系数
最大的点，则表明)(t x 信号中泄漏到达的时间晚于)(t y ，即泄漏距红色传感器远些，距蓝色传感器近些；若相关分析后计算的延迟时间为零，则表明中间相关，或为噪声干扰或表明泄漏位于两个传感器正中，应移动一个传感器位置后重新相关分析。
泄漏位置的计算：
对于不同管材、不同外径的管道，声音传播速度不同，根据事先给定的速度表，可以知道，两个传感器间的距离也可测的。
假设泄漏点距较近的传感器距离为L ，两个传感器间的距离为D ，声波的速度为V ，泄漏声波到达两个传感器的时间差为Td ，则泄漏点距较近的传感器距离计算公式如下：
)(2d T V L D ⨯+=
)(5.0d T V D L ⨯-=
当某管道只有一个管件可利用时，并且可能的泄漏点位于此管件与管道端部
或关闭的阀门之间时，这时只能放置一个振动传感器，只能采集一列声波，当存在泄漏点时，泄漏声一方面直接经管道传递到传感器，另一方面该波形经管端或关闭的阀门反射后再沿管道传递到传感器，因此该传感器所采集的信号在某一点获得泄漏噪声信号，且一定延迟后，该信号被叠加，通过对该信号进行自相关分析，即可确定该延迟时间。通过该管件声波的传播速度、传感器距管头或关闭的阀门的距离，可以计算出泄漏点距传感器的距离。
设x t ()为所测量的两列信号，则其自相关函数计算公式如下：
)
()()(1lim )()(1lim )()(1lim )(000τττττ-=-=-=+=⎰⎰⎰∞→∞→∞→x T T T T T T x R dt t x t x T
dt t x t x T dt t x t x T
R 自相关系数计算公式为：
202max 22max )()(1
)()(x T x x x
x x dt t x t x T R σμτσμττρ⎰-+≈-==
由于经管头或关闭的阀门反射的泄漏声波一定滞后于直接传递的泄漏声波，
所以只需计算该列信号的相关系数)(τρx ，相关系数最大对应的延迟时间即为泄
漏声直接传递和反射后传递的到达时间差，设该时间差为Td ，传感器距管头或关闭的阀门距离为D ，声波的速度为V ，则泄漏声距传感器的距离L 为：
D T V L d =⨯+2
d T V D L ⋅⨯-=5.0
二． 现场应用中需要解决的问题
在实际应用中，需要克服各种干扰，如管道流水的声音，三通，四通处的分流，汽车经过的电磁干扰和震动噪声，附近工地等其他噪声以及由于没有适合的放置传感器的位置而造成传感器距离太远，有效噪声太微弱，信噪比低等都会干扰正常的噪声信号。另外，管道的材质，粗细，管壁的厚薄，液体压力的大小，漏点的大小，管线所在区域的土质等也将对结果产生重要影响。综上所述，要想获得正确可信的结果，仪器必须具备完善的抗干扰措施，操作者也需要有一定的
