CN1261744C - 基于压力信号检测油气管道泄漏的方法 - Google Patents

Abstract

本发明公开了一种基于压力信号检测油气管道泄漏的方法，属于管道泄漏监测技术领域。该方法包括以下过程：采集非平稳状态下的压力变送器输出值序列；采用包括批处理法和递归处理法对非平稳信号加以变换处理；以变换的平稳压力信号序列采用序贯概率比检验法递推公式计算检验参数；判断泄漏是否发生，当检验参数大于报警上限时，则判断管道发生了泄漏，从而可以实现对油气管道泄漏检测的目的。本发明的优点在于：检测精度高，灵敏度好，成本低廉，运行可靠；特别适用于油气管道输送过程中启泵或停泵阶段压力变送器输出值处于非平稳条件下的管道泄漏检测。

Description

基于压力信号检测油气管道泄漏的方法

                            技术领域

本发明涉及一种基于压力信号检测油气管道泄漏的方法，属于管道泄漏检测技术领域。

                            背景技术

长距离输送管道由多个泵站及连接泵站的管道组成，泵站上有多种测量仪表，其中压力变送器是每个泵站都有的仪表，它们被安装在管壁上，直接测量管道内部压力，一般有较高的测量精度，压力变送器的输出值能反映出管道输送状况的变化，因此管道泄漏检测大多采用压力变送器的输出值作为判断的依据。在油气管道输送过程中，由于环境因素如土壤的腐蚀或人为因素如打孔盗油造成管壁破裂，发生油气泄漏，管内压力相应发生变化，使得压力变送器的输出值变化，根据这种变化，利用各种数学算法判断泄漏发生就是管道泄漏检测技术。

利用压力变送器的输出值判断泄漏发生最初是阈值的方法。如果管道正常输送时管内压力是5MPa，选择阈值为0.1MPa，即压力变送器的输出值小于4.9MPa时认为管道上发生了泄漏，系统就给出报警信号，提醒泵站工作人员进行处理。阈值法很容易实施，多数泵站采用此方法进行泄漏的判断，但它也有明显的缺点，现场条件下压力变送器的输出值受到许多干扰，比如管壁的震动、温度的变化，干扰因素使输出值的偏差可能超过阈值而产生误报警，这使阈值法的应用遇到了一定的限制。

针对这种情况，美国壳牌石油公司发明了序贯概率比检验法检测泄漏，并申请了专利。该方法仍然利用压力变送器的输出值来判断泄漏，不同于阈值法的是利用序贯算法对压力变送器的输出值进行了变换，当变换值超过设定的报警限时，系统才认为泄漏发生。举例来说明序贯算法，正常输送时管内压力是5MPa，由于干扰因素的存在，压力变送器的输出值不可能总是5MPa，但干扰因素往往是随机的，即压力变送器的输出值可能是4.9MPa，也可能是5.1MPa，如果将输出值相加后取平均其结果仍是5MPa，而不会超出阈值的范围，因此就减少了误报警率。下面说明现有的基于压力信号检测油气管道泄漏的方法：

(1)采集一段时间内压力变送器的输出值记为z(k)，比如100秒作为一个时间单位，k等于1到100；

(2)计算该时间段内均值记为μ，Δμ是压力变化的阈值，由现场情况决定其取值，μ的计算公式为μ＝(∑z(k)/n。其中∑为求和符号，即计算100秒内压力变送器输出值的总和，n等于100；

(3)定义L(k)为序贯算法的输出值，z(k)是其输入值，记L(k)的初始值L(0)为0，L(k)的计算公式为

              L(k)＝L(k-1)+z(k)·Δμ-μ·Δμ-(Δμ·Δμ)/2       (1)

(4)作出管道是否发生泄漏的判断：给出L(k)的报警限，如果计算出的L(k)大于报警限，认为该段时间内有泄漏发生；如果L(k)始终小于报警限，认为该段时间内无泄漏发生；

该方法无须对原有的SCADA系统进行任何改动，即可以进行在线检测，检测成本低，精度高，其泄漏检测灵敏度下限可以达到仪表精度的一半，因此这是一项很有前景的技术。

以上方法基于压力平稳值情况下适用，即均值是不随时间变化的μ＝a(a是常数)，当管道输送处于平稳状态时，这种情况是容易满足的，因此可以利用以上方法进行计算。然而实际上管道的压力并非平稳输出，特别在泵的启动和停泵过程，其压力输出有线性增加和线性减少的趋势，此时计算公式(1)的μ就不是一个常数，而是μ＝a+b·k其中a为截距，b为斜率，k为采样时刻，如果仍然按照公式(1)计算，就会产生线性误差，甚至无法应用，因此必须将压力检测信号加以处理，才能进一步判断泄漏。

                        发明内容

本发明的目的在于提供一种基于压力信号检测油气管道泄漏的方法，该方法特别适用于输油过程中启泵或停泵阶段压力变送器输出值处于非平稳条件下的泄漏检测，具有检测精度高，系统容易实现，运行可靠等特点。

本发明是通过下述技术方案加以实现的，该方法包括以下过程：

1、采集非平稳状态下的压力变送器输出值序列z(k)；

2、对非平稳信号加以变换处理：

当采样全部完成后集中处理时，采用1)批处理的方法，采用一元线性回归法得

              z(k)＝a+b·k+ε(k)                                    (2)

式中z(k)为采集的压力信号序列，ε(k)为不带线性趋势的压力信号序列，a，b是线性系统未知参数，也就是待求的量。由最小二乘法求得未知参数a，b的值，得到一条拟合的直线，截距a，斜率为b，公式为

              f(k)＝a+b·k                                          (3)

由公式(2)和(3)的差得到ε(k)，然后转入第3步再处理。

或当每进行一次采样计算一次ε(k)值时，采用2)递归的处理方法，采用卡尔曼滤波器法，系统状态方程为

              x(k)＝x(k-1)                                          (4)

             z(k)＝H(k)x(k)+v(k)                            (5)

其中x(k)为过程状态变量，其值为(a，b)；z(k)为采集的压力信号序列；H(k)取值为(1，k)；v(k)是不带线性趋势的压力信号序列，同ε(k)。将公式(4)(5)递推得到公式(3)，进而得到ε(k)，转入第3步再处理。

3、以变换的平稳压力信号序列采用序贯概率比检验法递推公式计算检验参数L(k)：

              L(k)＝L(k-1)+ε(k)·Δμ-(Δμ·Δμ)/2       (6)

√公式(6)中Δμ是压力变化的阈值；

4、判断泄漏是否发生：当检验参数L(k)大于报警上限时，则判断管道发生了泄漏，从而可以实现对油气管道泄漏检测的目的。

本发明的特点主要体现在：

1、本方法特别适用于油气管道输送过程中启泵或停泵阶段压力变送器输出值处于非平稳条件下的管道泄漏检测。

2、利用序贯概率比检验法本质上是检验均值的变化，当均值的变化大于仪表精度的一半时就可以检测出泄漏，本发明检测精度高，灵敏度好。

3、当运行参数发生变化时，该方法能很好地跟踪变化，因此适用性较强；本方法可以直接嵌入SCADA系统中，不需要另外的资源，成本低廉，运行可靠。

                        附图说明

图1：数据采集器采集的压力变送器输出值曲线。

图2：图1中压力变送器的输出值曲线去除线性趋势后的结果图。

图3：由图2截取的带有泄漏信号的平稳压力信号图。

图4：由图3经序贯概率比检验法变换后的输出信号图。

图中，横坐标均表示测量压力信号点的个数。

                        具体实施方式

下面结合附图和具体实施方式对本发明作进一步详细的说明：

在泵的启动阶段，采集压力变送器的输出值，得到图1的上升曲线，其中压力输出值从0.094MPa到0.101Mpa，压力差达到0.007Mpa带有明显的上升趋势。应用序贯法检测泄漏必须由公式L(k)＝L(k-1)+z(k)·Δμ-μ·Δμ-Δμ·Δμ/2计算泄漏检验参数，其中μ应该是常数，即μ＝a(a是常数)。图1中μ值显然不是常数，而是一条截距为a，斜率为b的直线，用公式表示为μ＝a+b·k。因此必须算出直线中a，b的值，代入式中才可以利用序贯算法计算泄漏。本发明就是针对这种情况进行了改进，克服了序贯概率比方法应用中的一大障碍。

计算a，b的值，其一是批处理的方法，即利用最小二乘法对一组数据进行处理，拟合出直线方程，得到a，b的值。图1的数据按本方法计算后得a，b的值分别为0.0943和9.4076e-008；其二是递归的方法，即利用卡尔曼滤波器法，每采集一次压力信号，判断一次拟合曲线的截距与斜率，逐点连接后得到拟合直线方程，同样得到a，b的值。图1的数据按照本方法计算后得到a，b的值分别为0.0943和9.4076e-008。批处理的方法是对一批数据进行综合处理，其结果准确率高，但其判断方法不具有实时性，而递归的方法克服了其不足，适用于在线实时监测。将a，b代入公式(2)中就可以得到不带线性趋势的值ε(k)，也就是图2中所示曲线。

由图2中可以看出在40000点之前，信号均值在0值附近，50000点之后，信号均值在-10附近，均值的这种变化说明管道中有泄漏发生，可以利用序贯法来判断。为了更清楚的表示序贯算法的处理过程，我们取泄漏前后共1800个采样点来说明，如图3所示，取第一个值作为横坐标的原点，共1800个坐标点，纵坐标表示变换后的压力值，其并不是稳定在一个值上，由于干扰的存在，其值在均值上下波动，如果利用阈值法，就会把一些大的波动认为是泄漏发生，则会造成误报警。图3中的值作为输入值，经公式(6)计算后得到输出值L(k)，就是图4所示曲线，其曲线是光滑的，不存在大的波动，避免了阈值法中误报警的情况，图4中给出报警上限为0.08e-004，如图中水平直线所示，在1600点之后L(k)曲线值超过上限，系统产生泄漏报警。

Claims (1)

1.一种基于压力信号检测油气管道泄漏的方法，其特征在于包括以下过程：

1)采集非平稳状态下的压力变送器输出值序列z(k)；

2)对非平稳信号加以变换处理：

当采样全部完成后集中处理时，采用(1)批处理的方法，采用一元线性回归法得

               z(k)＝a+b·k+ε(k)                          (2)式中z(k)为采集的压力信号序列，ε(k)为不带线性趋势的压力信号序列，a，b是线性系统未知参数，也就是待求的量；由最小二乘法求得未知参数a，b的值，得到一条拟合的直线，截距a，斜率为b，公式为

               f(k)＝a+b·k                                (3)由公式(2)和(3)的差得到ε(k)，然后转入第3)步再处理；

或当每进行一次采样计算一次ε(k)值时，采用(2)递归的处理方法，采用卡尔曼滤波器法，系统状态方程为

               x(k)＝x(k-1)                                (4)

               z(k)＝H(k)x(k)+v(k)                         (5)其中x(k)为过程状态变量，其值为(a，b)；z(k)为采集的压力信号序列；H(k)取值为(1，k)；v(k)是不带线性趋势的压力信号序列，同ε(k)；将公式(4)(5)递推得到公式(3)，进而得到ε(k)，转入第3)步再处理；

3)以变换的平稳压力信号序列采用序贯概率比检验法递推公式计算检验参数L(k)：

              L(k)＝L(k-1)+ε(k)·Δμ-(Δμ·Δμ)/2      (6)

其中Δμ是压力变化的阈值；

4)判断泄漏是否发生：当检验参数L(k)大于报警上限时，则判断管道发生了泄漏，从而可以实现对油气管道泄漏检测的目的。

Images