基于瞬变流的爆管定位方法与流程

1.本发明涉及一种爆管定位方法，属于工程学领域，具体讲,涉及基于瞬变流的爆管定位方法。


背景技术：

2.供水管网是城市基础设施的重要组成部分，被称为城市的生命线。由于不可避免的老化、侵蚀、人为破坏等原因，爆管时有发生，且具有很强的随机性和不确定性，很难被发现。因此，这给我们的生产和生活带来了严重的影响。因此，实现爆管定位具有重要意义。
3.近年来，scada系统被广泛引入供水管网中，对流量和压力进行实时监测。压力的scada系统主要分为两类：(1)低频scada系统：只能监测压力的整体变化；(2)高频scada系统：既能监测压力的瞬时变化—瞬变流异常，又能监测压力的整体变化—系统异常。基于这两类scada系统的爆管定位方法都有大量研究。基于低频scada系统的爆管定位方法只能将爆管定位到某个区域而不能定位到点；基于高频scada系统的爆管定位方法只适用于简单管网，对于复杂管网无法定位。因此研究一种即适用于复杂管网又能将爆管定位到点的爆管定位方法是有必要的。
4.在高频scada系统中，对爆管前后的压力数据进行处理，通过对两类异常的判别，对爆管事件进行识别；然后在管网中每个管段设置不少于3个爆管点，对监测点瞬变流进行数值模拟，通过与实际监测的瞬变流进行对比确定爆管管段；最后在爆管管段等间距设置多个爆管点，再一次对监测点瞬变流进行数值模拟及与监测点瞬变流进行对比，对爆管进行定位。


技术实现要素：

5.为克服现有技术的不足，针对复杂供水管网中爆管难以被精确定位的问题，本发明旨在提出一种复杂供水管网爆管精确定位的方法，即通过高频scada系统监测的压力异常信号，对爆管事件进行识别，然后通过在监测点模拟瞬变流与实际监测的瞬变流进行对比，确定爆管管段和实现爆管定位。
6.为此，本发明采用的技术方案如下：
7.一种基于瞬变流的爆管定位方法，步骤如下：
8.s1、根据监测到的高频压力异常特征，识别爆管事件；
9.s2、根据转动阀门的实际操作，数值模拟爆管激发的瞬变流；
10.s3、根据瞬变流传播特征，数值模拟监测点瞬变流；
11.s4、应用监测点瞬变流数值模拟方法，判别爆管管段；
12.s5、再次应用监测点瞬变流数值模拟方法，实现爆管定位。
13.进一步地，在步骤s1中，爆管事件识别的具体方法为：
14.对在高频scada系统中监测到的压力异常数据进行滤波处理，分析其系统异常(压力的整体变化)及瞬变流异常(压力的瞬时变化)，根据其异常特征对爆管事件进行判别。
15.进一步地，在步骤s2中，数值模拟爆管激发瞬变流的具体方法为：
16.在实验中，通过快速打开球阀模拟爆管事件，根据动量定理及打开球阀的实际操作数值模拟爆管激发的瞬变流，其振幅表示式如下：
[0017][0018]
式(1)中，ρ为水密度,a为瞬变流传播速度,v(t)为水流速度,a
c
爆管管段的横截面积，h为爆管处水压,c
d
为泄漏系数,g为重力加速度，μ为爆管系数，a(t)为球阀打开面积(爆管面积)，其表达式如下：
[0019][0020]
式(2)中，r为爆管面半径，α(t)为球阀转杆旋转的角度，随时间的变化关系如下式：
[0021][0022]
式(3)中，c1和c2分别是球阀转动的两个阶段的角加速度，t为时间。
[0023]
进一步地，在步骤s3中，数值模拟监测点瞬变流的具体方法为：
[0024]
激励信号从爆管位置向监测点传播，传播过程存在多个路径，且每个路径中均存在瞬变流的局部损失及摩阻损失，故每条路径中传播至监测点的瞬变流振幅表达式为：
[0025][0026]
式(4)中，k为瞬变流的传播路径序号，t
k
为瞬变流在第k条路径的传播时间，a
ki
第i个节点的透射或反射系数，r
kj
第j段管段单位长度的摩阻因子，e
‑
rkj
第j段管段单位长度的摩阻损失，l
kj
为第j段管段的长度。
[0027]
将来自每条路径瞬变流进行叠加，即监测点的瞬变流，表达式为：
[0028][0029]
运用上述理论，在爆管点输入s2步骤获得的爆管激发瞬变流，在监测点位置可计算出监测点瞬变流的幅值，即完成瞬变流的数值模拟。
[0030]
进一步地，在步骤s4中，判别爆管管段的具体方法为：
[0031]
在管网中每个管段设置不少于3个爆管点，应用s3步骤中的监测点瞬变流数值模拟方法，对监测点瞬变流进行数值模拟，然后计算模拟信号与实际监测信号的误差，则最小误差所在管段被认定为爆管管段。
[0032]
进一步地，在步骤s5中，爆管定位的具体方法为：
[0033]
在s4步骤中判定的爆管管段中等间距设置多个爆管点，再一次应用s3步骤中的监测点瞬变流数值模拟方法，对监测点瞬变流进行数值模拟，并计算拟信号与实际监测信号的误差，则误差最小位置即为爆管位置。
附图说明
[0034]
图1为实验管网示意图
[0035]
图2为爆管事件前后监测点的压力异常曲线图；
[0036]
(a)为爆管事件前后监测点p1的压力异常曲线图
[0037]
(b)为(a)图异常部分在横轴方向上的放大显示
[0038]
(c)为爆管事件前后监测点p2的压力异常曲线图
[0039]
(d)为(c)图异常部分在横轴方向上的放大显示
[0040]
(e)为爆管事件前后监测点p3的压力异常曲线图
[0041]
(f)为(e)图异常部分在横轴方向上的放大显示
[0042]
图3为实验管网各管段3个位置爆管模拟监测信号与实际监测信号的误差分布图；
[0043]
图4为爆管管段不同位置爆管模拟监测信号与实际监测信号的误差分布图；
[0044]
图5为本发明流程图。
具体实施方式
[0045]
下面结合附图及具体实施例对本发明做进一步说明，但下述实施例绝非对本发明有任何限制。
[0046]
以下通过在实验管网中实现爆管定位，具体操作步骤如下：
[0047]
s1、爆管事件识别：
[0048]
实验设备主要包括：水泵、水箱、压力表、高频压力传感器(3个)、dn100、dn50、dn20钢管、dn50pvc管、水阀，各管段组合成具有2个环的管网，高频压力传感器采集频率为10000hz。用一个字母“b”和一个数字对管网内各节点进行编号，例如：b1、b5、b10。各管段编号用一个字母“t”和管段两端节点数字编号进行编号，例如：t2
‑
3、t10
‑
13(图1)。在实验管网，通过快速打开水阀，制作一次爆管事件，监测事件前后的压力数据。
[0049]
实现爆管事件识别，具体如下：
[0050]
对监测到的压力数据进行滤波处理后绘制成曲线(图2)。从图2(a)、图2(c)和图2(e)中可以看出，在30s附近，压力曲线有明显的下降系统异常，且伴有明显的负向突跳异常，经横轴方向放大显示(图2(b)、图2(d)、图2(f)，该异常为瞬变流异常。结合以往总结的爆管事件压力曲线异常特征，判定该事件为爆管事件。
[0051]
s2、数值模拟爆管激发的瞬变流：
[0052]
根据动量定理及上述打开球阀的实际操作数值模拟爆管激发的瞬变流，其振幅表示式如下：
[0053][0054]
式(1)中，ρ为水密度,a为瞬变流传播速度,v(t)为水流速度,a
c
爆管管段的横截面积，h为爆管处水压,c
d
为泄漏系数,g为重力加速度，μ为爆管系数，a(t)为球阀打开面积(爆管面积)，其表达式如下：
[0055][0056]
式(2)中，r为爆管面半径，α(t)为球阀转杆旋转的角度，随时间的变化关系如下
式：
[0057][0058]
式(3)中，c1和c2分别是球阀转动的两个阶段的角加速度，t为时间。
[0059]
s3、数值模拟监测点瞬变流：
[0060]
激励信号从爆管位置向监测点传播，传播过程存在多个路径，且每个路径中均存在瞬变流的局部损失及摩阻损失，故每条路径中传播至监测点的瞬变流振幅表达式为：
[0061][0062]
式(4)中，k为瞬变流的传播路径序号，t
k
为瞬变流在第k条路径的传播时间，a
ki
第i个节点的透射或反射系数，r
kj
第j段管段单位长度的摩阻因子，e
‑
rkj
第j段管段单位长度的摩阻损失，l
kj
为第j段管段的长度。
[0063]
将来自每条路径瞬变流进行叠加，即监测点的瞬变流，表达式为：
[0064][0065]
运用上述理论，在爆管点输入s2步骤获得的爆管激发瞬变流，在监测点位置可计算出监测点瞬变流的幅值，即完成瞬变流的数值模拟。
[0066]
在实验管网的管段上设置爆管点，将式(1)～式(3)获得的爆管激发的瞬变流代入式(4)，获得每一个路径的监测点瞬变流，再将其代入式(5)，其结果则为模拟监测点瞬变流。
[0067]
s4、判别爆管管段：
[0068]
在实验管网中每个管段设置不少于3个爆管点，应用s3步骤中的监测点瞬变流数值模拟方法，对监测点瞬变流进行数值模拟，然后计算模拟信号与实际监测信号的误差(图3)，从图3中可以看出t9
‑
10管段存在最小误差，故认为该管段为爆管管段。
[0069]
s5、爆管定位：
[0070]
在s4步骤中判定的爆管管段中等间距设置多个爆管点，再一次应用s3步骤中的监测点瞬变流数值模拟方法，对监测点瞬变流进行数值模拟，并计算拟信号与实际监测信号的误差(图4)，从图4中可以看出距离b9节点2.3m处误差最小，故认为该位置为爆管位置(表1)。
[0071]
表1：本发明和负压波方法爆管定位结果与实际位置对比表
[0072][0073]
从表1中可以看出，本发明爆管定位位置存在0.2m的绝对误差和2.91％的相对误差；与负压波定位方法相比，本发明定位精度绝对值提高0.98m，相对值提高14.27％；综上，本发明爆管定位的绝对误差和相对误差都较小，解决了复杂管网的精确定位问题。