一、泄漏监测与定位的原理
随着管道的建设,各种监测技术也在不断发展,目前应用于管道泄漏监测的主要有压力梯度法、负压力波法、流量平衡法、超声波检测法等物理方法和一些化学方法,这些方法的特点和应用场合各不相同。国外虽然有比较成熟的技术,但是不能够适应中国原油输送的特点,而且价格昂贵。在国内,多家大学和油田均做过一定研究,但到目前为止,只有天津大学研制的系统成功地应用于现场,取得了显著的经济效益。
我国近年来的原油泄漏主要是因为人为破坏造成的,特点是持续时间短、泄漏量较大,属于突发性事故,我们采用负压力波检测法。1997年我们为华东石油管理局濮阳-滑县的管线设计了一套泄漏监测系统,是在DOS操作系统下用C语言实现的,其主要缺点是无法实现多任务,数据采集和通信不能同时进行。从2000年初至今,我们使用LabVIEW平台,设计并不断完善了基于虚拟仪器的原油管道泄漏监测系统,已经成功地安装于临盘-济南和沧州-临邑的两条管线上。

负压力波法是一种声学方法,所谓压力波实际是在管输介质中传播的声波。当管道发生泄漏时,由于管道内外的压差,泄漏点的流体迅速流失,压力下降。泄漏点两边的液体由于压差而向泄漏点处补充。这一过程依次向上下游传递,相当于泄漏点处产生了以一定速度传播的负压力波。根据泄漏产生的负压波传播到上下游的时间差和管内压力波的传播速度就可以计算出泄漏点的位置。定位的原理如图1所示,L为管道长度,X为泄漏点,t1,t2为负压波传播到上下游的时间。图2为采集到的负压波波形,纵坐标单位为MPa。

图2 负压波波形

常规的负压波法定位公式为：

其中K为液体的体积弹性系数,r为液体的密度,E为管材的弹性模量,D为管道直径,e为管壁厚度,C1为与管道约束条件有关的修正系数。其中体积弹性系数K和密度r都是温度的函数。所以采集压力的同时还要采集管道的温度用来校正压力波的传播速度。
流量在监测中也是一个重要的参数。因为我国很多原油管道都不是密闭输送,收油端直接接入大罐,本身压力就很低,再加之大罐的滤波作用,压力的变化可能会淹没在一片噪声之中。若同时监测流量的变化,也可以及时发现泄漏,弥补压力灵敏度低的缺点。另外,除了泄漏,泵站内部的操作,如调泵、调阀,也会产生负压波。我国的现状要求只能在站内安装一台压力变送器,所以仅靠压力就不能判断负压波是站内产生的还是站外产生的。而流量变化对于泄漏和站内操作具有不同性质,通过分析流量变化就可以甄别负压波的来源。
基于以上考虑,我们采取压力流量联合判断的方法,分别在管道的出口和进口加装合适量程的压力变送器、温度变送器和流量变送器,不间断采集原油的压力、温度和流量,监视管道的运行状况。压力变送器和温度变送器输出标准的4~20mA电流模拟信号,流量变送器输出的是脉冲。由于模拟量的采集速率不是很高,大约为几百Hz,脉冲的频率约为700多Hz,我们选用低价的多功能卡PCI-6023E,利用其两路模拟量单端输入和一个通用24位计数器。从功能、价格和开发周期等因素来考虑,这款采集卡都非常适合我们的要求。对于微弱信号的提取效果明显,可以观察到万分之几MPa的微小变化。在临-济线和沧-临线上共使用了10块PCI-6023E,连续工作的稳定性很好。
为了节约存储空间,我们采用二进制格式保存文件。压力和温度每一小时一个文件,流量因为是累积量,每一天一个文件,供泄漏定位和查阅历史数据之用。

二、精确捕捉压力波传播到上下游的时间差
精确获得泄漏引发的压力波传播到上下游传感器的时间差,需要准确地捕捉到泄漏负压力波信号序列的对应特征点。由于不可避免的工业现场的电磁干扰、输油泵的振动等,采集到的压力波信号序列附加了大量噪声,如何从噪声当中准确地提取出信号的特征点是定位的关键。
首先,因为不同条件（管径、温度、压力、输送工艺、环境噪声等）下的信号具有不同的特征,我们利用LabVIEW中丰富的信号处理函数对信号的特征作了深入的研究分析和预处理,使得系统能够针对不同的信号做出相应的处理。例如,在沧州首站,启泵的时候信号中含有比较有规律的尖锐噪声,虽然幅度并不大（只有约千分之一兆帕）,却因为压力比较低而非常明显,如图3。针对这样的情况,我们使用了中值滤波器。中值滤波器能够平滑尖锐噪声,又能保持信号的边沿,而且计算速度非常快,可以用于在线处理。对比下面两图,可以看到滤波的效果。

图4 上图的信号经过中值滤波器之后的波形

对信号作了一定预处理之后,我们采用Signal Processing软件包里提供的小波分析的动态库来捕捉压力波变化的特征点。将含有负压波信号的数据通过两通道滤波器,分为低频概貌和高频细节两路输出,负压波的位置就是高频细节当中最小值对应的索引值。图5为图2所示信号的高频结果。因为我们仅分析信号,不作重构,所选择参数设计更为灵活的双正交的滤波器组,而不必考虑对偶滤波器组是否收敛。利用Wavelet and Filter Bank Design Tookit提供的Design Panel 和1D Test Panel,我们在现场对于不同的信号可以方便地调整滤波器参数,观察滤波效果,大大提高了应用的效率。

图5 图2的信号通过两通道滤波器后的高频结果

三,系统的软件设计
以沧州-临邑长输管线的泄漏监测系统为例介绍一下整个系统的设计。
系统的组成如图6,由安装于沧州、东光、德州和临邑的四套装置和沧州处调度室的中心计算机组成。各站的装置包括各种传感器、计算机、信号调理装置、信号采集卡和调制解调器等。
因为管道全线长180公里,两座泵站之间的距离一般是60公里左右,所以只能通过电话网络实现数据通信。LabVIEW的函数库里没有现成的函数可供调用,我们综合运用Serial Port VIs,In Port, Out Port 和Call Library Function调用自己编写的动态链接库,控制调制解调器实现远程的数据通信,各个子站的数据都可以实时传到中心站,线路中断能够自动重新连接。主程序每一分钟调用一次泄漏判断子程序,该子程序综合运用负压波法、压力梯度法和流量差法分析采集到的工况数据,判断是否有泄漏发生。
沧州站和沧州处调度室仅有几百米的距离,两台计算机之间建立了以太网,使用TCP/IP实现数据传输。现场还有其他PLC系统,我们用DDE实现LabVIEW和Excel之间的数据交换,作为与其他系统的接口。
用户对系统的操作都设计为菜单操作,包括参数设置、文件