一种基于阀门操作和在线用水计量的供水管网漏失高效定位方法与流程

本发明涉及市政工程和城市供水管网领域，具体涉及一种基于阀门操作和在线用水计量的供水管网漏失高效定位方法。

背景技术：

漏失是我国城镇供水管网运行中普遍存在的问题之一。根据《城市供水统计年鉴(2017)》的调查统计，我国主要省市的漏失水量大致占供水量的8％-50％，部分地区的漏失情况相当严重。管网漏失不仅严重浪费了水资源和能源，同时也会引起饮用水安全隐患。因此，供水管网的漏失检测定位一直是供水企业的重要任务之一，对于节约水资源和保障供水安全具有重要意义。供水企业通常采用专业探漏设备来查找漏失，如听漏棒、漏失噪声相关仪、探地雷达等。但是这类设备方法普遍存在成本高、工作量大、效率低的问题，只适用于小范围的漏失精确定位。考虑设备探漏方法的这种特点，一种经济有效的漏失定位实施策略是先通过其他漏失定位方法缩小漏失区域，然后采用设备检测的方法在所确定的漏失区域内进行漏失精确定位。

近年来，随着智能管网技术的快速发展，智能设备已在供水管网中广泛应用，比如智能水表、远控阀门等。其中，智能水表可以实时采集、存储和传输终端用户的水量数据，实现了对管网中用户用水量的在线准确计量。这种基于智能水表的在线用水计量方式可以实现在线水平衡分析，为开发高效漏失定位方法提供了新的机遇。基于此，设计开发高效准确的漏失定位方法以尽量缩小需探测的漏失区域，可以弥补现有漏失定位方法的不足，有效提高漏失定位效率，从而减少供水管网漏失。

国际水协(iwa)推荐采用区域独立计量(dma)方法进行漏失监测定位，即通过将复杂供水管网划分为相对较小的多个dma，并独立计量各dma的流量变化以判断其是否存在漏失。一般来说，当一个dma的最小夜间流量出现明显异常时(如数值超过设定阈值)，说明该dma中产生了新的漏失。然而，该方法只能识别漏失范围至dma水平(通常包含几十公里范围)，而不能准确确定dma中漏失位置。为进一步缩小漏失区域以便于应用设备方法进行精确定位，供水领域的学者和工程师们开发了水力模拟、优化求解、数据分析等多种方法。这类方法在理论上可以有效缩小漏失区域，然而在实际应用中尚存在一些明显问题。一方面，这类方法通常需要建立高精度的水力模型或获取大量准确的实测数据，并且需要专业技术人员进行应用，对供水企业要求过高，因而实用性较差。另一方面，这类方法一般只能定位到可疑的漏失区域(范围相对较大且分散)，准确度难以保证。因此，仍需进一步开发实用高效准确的漏失定位方法。

技术实现要素：

针对上述所提出的漏失定位方法实用性和准确度不足的问题，本发明提出一种基于阀门操作和在线用水计量的供水管网漏失高效定位方法，该方法能够准确高效地确定dma中的漏失区域，且易于实施。

为解决上述技术问题，本发明采用如下技术方案：

一种基于阀门操作和在线用水计量的供水管网漏失高效定位方法，通过优化dma中现存阀门的关阀策略将其划分为两个子区域，并基于智能水表的在线水平衡分析来判别漏失位于哪个子区域内，由此缩小漏失区域；重复实施这种阀门操作和水平衡分析相结合的策略(即持续进行多个阶段)，不断缩小所识别到的漏失区域，直至没有阀门操作可以进一步缩小漏失区域为止。

上述技术方案中，进一步地，所述的定位方法具体包括以下步骤：

1.利用图论方法确定可将dma分解为两个子区域的关阀策略vs＝{v1,v2,...,vn}，形成关阀策略集合v＝{v1,v2,...,vs}。其中，n为dma中现存可操作的阀门总数，s为得到的关阀策略总数；vn∈{0,1}表示阀门n的状态，0代表关闭，1代表打开；关阀策略vs将dma分解为两个子区域可表示为g＝g1(vs)∪g2(vs)，其中g为dma管网图，图中端点代表管网中节点类元素(如用户节点、管段连接点、消火栓等)，边代表管网中线段类元素(如管段、阀门、水泵等)，g1和g2表示被关闭阀门所分解的两个子区域。

2.初始化漏失定位，将dma中所有管道标记为可疑漏失管道，同时设置漏失定位的当前阶段k＝1。本发明以管段表示漏失定位的对象，xm＝1表示可疑漏失管段，xm＝0表示无漏失管段，m＝1,2,...m，m表示dma中管段总数。与管段相连的其他组件处(如节点、消火栓等)的漏失与管段等同考虑。

3.确定当前阶段k的最优阀门操作策略。建立如下的优化目标函数以确定阶段k时的最优阀门操作，

最小化

约束条件：

关阀策略：

将dma分解为两个子区域：g＝g1(v^k)∪g2(v^k)(3)

不可断水用户的压力约束：其中，表示当前阶段k的阀门操作数(相对于上一阶段k-1的阀门状态)，包括阀门开启和关闭操作；表示阀门操作数值归一化。lm表示管段m的长度；表示由关阀策略v^k所分解的两个子区域中可疑漏失管道长度的差异，表示dma中管段总长度；表示可疑漏失管道长度差异数值的归一化。w1和w2表示阀门操作数和可疑漏失管段长度差异这两项的权重，本发明采用w1＝w2＝0.5以等效考虑这两项的权重。

公式(2)表示目标函数的决策变量，即各阀门在阶段k时的操作状态公式(3)表示关阀策略应将dma分解为两个子区域，在步骤1所得到的关阀策略集v即是满足公式(3)约束的可行解集。公式(4)表示管网中不可中断供水的用户节点的压力约束，即阀门操作后这些节点处的压力hu(v^k)不会低于所需最小压力以保证这些用户的用水需求。应用公式(4)时需使用供水管网水力模型进行模拟计算，以判断压力约束是否满足。

公式(1)表示优化目标为最小化关阀策略所产生的阀门操作数和两个子区域中可疑漏失管道长度差异。该目标函数的物理意义是在当前阶段采用最少的阀门操作数分解dma管网，同时尽可能均分可疑漏失管段长度。减少阀门操作数可以降低漏失定位工作量；均分可疑漏失管段长度(即传统二分法)可以减少漏失定位所需的实施阶段数量，提高漏失定位效率。采用第1步所得到的关阀策略集v来评价所定义的优化目标函数值(公式1)，选择满足约束条件(公式4)且目标函数值最小的关阀策略为当前阶段k的最优阀门操作策略。

4.根据最优阀门操作策略关闭和打开相应阀门，将dma分解为两个子区域。在这两个子区域中，有一个子区域应包含至少一个dma入口(即该子区域可正常运行)，而另一个子区域可能正常运行，也可能被关阀所封闭(即中断供水)。由于本发明方法要求在夜间用水时段实施，短暂中断供水并不会产生明显影响。

5.任选一个正常运行的子区域实施在线水平衡分析，以判断漏失是否位于该子区域内。若判定漏失不位于该子区域，则其位于另一个子区域。在线水平衡分析是通过dma的边界流量计和用户端的智能水表实现。进入选定子区域内的流量(入口处流量计在线读数减去出口处流量计在线读数)与选定子区域内用户用水量(智能水表在线读数之和)的差值，即非收益水量(nrw)，可表示所选定子区域内是否存在漏失。如果选定子区域的nrw与漏失定位前dma的nrw数值相近，则说明所选定子区域内存在漏失；否则，漏失位于另一个子区域内。

6.更新当前阶段k的管段标记，将判定不存在漏失的子区域内管段标记为无漏失管段(即x^k＝0)，以缩小漏失区域。

7.判断dma中是否有阀门操作可以进一步缩小漏失区域。若有，则回到第3步进行下一阶段的漏失定位(即k＝k+1)；否则说明已找到最小漏失区域，则停止漏失定位过程。本发明所定位的最小漏失区域为漏失定位方法停止后dma中标记为1的管段及其相连组件。

8.根据漏失定位结果，派遣工作人员至所定位的漏失区域内采用人工听漏法等设备方法精确查找漏点位置，并修复所查找到的漏点。

更进一步地，所述的步骤1具体包括如下步骤：

(1)移除管网图g中所有代表阀门的边，识别所产生新图中的所有连通子图。这些连通子图代表管网中通过关闭阀门可封闭的最小区域(包括节点、管段、消火栓等)，也是本发明方法所能定位到的最小漏失区域。

(2)将管网图g转换为新的替代图gv。以管网中阀门为边，上一步中所识别到的所有连通子图为端点，根据阀门与连通子图之间的拓扑连接关系构建新的替代图gv。

(3)确定替代图gv的割集，每个割集表示一种关阀策略vs，割集中的边代表关阀策略中需关闭的阀门(即vn＝0)；将每个割集转换为关阀策略即可形成关阀策略集v。

更进一步地，所述的步骤(3)中，采用karger算法中的contraction算法来确定替代图gv的割集：随机选择图中的一条边将其移除，合并该边的两个端点，保留端点的其他连接边；重复执行该过程直至图中只剩下两个端点，此时这两个剩余端点的连接边即是图的一个割集。重复执行contraction算法至t⁴logt次可得到足够数量的割集(t为图gv中端点总数)。

本发明可在保证漏失定位准确度的情况下明显缩小漏失区域，提高漏失定位的效率。另外，本发明的漏失定位方法要求在夜间用水时段实施，比如02:00～05:00，以减少对dma区域内用户用水的干扰。

本发明所具有的有益效果在于：

1.该方法可采用最少的阀门操作数确定最小的漏失区域(关闭阀门可封闭的最小区域，通常小于5.0km)，定位效率高。相对于传统人工听漏方法，本发明方法可提高漏失定位效率达15-20倍。

2.相对于很多漏失定位方法的准确度不足的问题，该方法通过优化阀门操作和在线水平衡分析逐步缩小漏失区域，可确保漏失定位的准确性。

3.该方法只需操作阀门和进行水平衡分析(普通技术人员即可掌握)，而无需依赖水力模型或复杂数据分析，易于在实际管网中实施。

附图说明

图1为本发明方法的实施路线图。

图2为本发明所采用图论方法的示意图。

图3为本发明的多阶段漏失定位过程示意图。

图4为实施例dma的管网拓扑结构图。

图5为实施例中漏失定位过程中各阶段阀门操作图。

图6为实施例中漏失定位过程中水平衡分析结果图。

具体实施方式

为使本发明所提漏失定位方法易于明白，以下结合附图和具体实施例对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

如图2(a)表示一个简单的dma管网，包含2个入口、1个出口、5个阀门(v1，v2，v3，v4和v5)及若干管段和节点。本发明采用如下具体实施步骤对其进行漏失定位，参见图1。

步骤1：将dma管网视为管网图g(管段和阀门为边，节点为端点)，移除图g中的阀门边，识别图中的所有连通子图。对于图2(a)的示例，识别到5个连通子图，分别以s1、s2、s3、s4和s5表示，如图2(b)所示。每个连通子图代表关阀可分割的管网组件集合。

步骤2：将管网图g转换为替代图gv(阀门为边，连通子图为端点)。对于图2(a)的示例，生成的替代图如图2(c)所示。

步骤3：运行contraction算法t⁴logt次以确定替代图gv的割集(t为图的端点总数)，形成关阀策略集v。对于图2(c)所示替代图，图2(d-f)表示运行contraction算法一次的结果。首先随机选择一条边移除(图2d中的v5)；然后合并所选边的两端节点s4和s5，生成新的端点s45(图2e)；再次随机选择边进行移除和端点合并操作，直至剩下两个端点，即可得到替代图gv的一个割集(图2f所示割集为{v3,v4})。重复运行该算法1005次后，删除割集中的重复值，可得到图gv的割集集合。将割集集合中边转换为关闭的阀门，得到关阀策略集合v。

步骤4：初始化漏失定位程序，设置当前阶段k＝1，并将所有管段设置为可疑漏失管段(即x^k＝1)。图3(a)表示初始化结果，图中标红管段为可疑漏失管段。

步骤5：确定当前阶段的最优阀门操作策略。利用步骤1～3所确定的关阀策略集v评价优化目标函数f^k(公式1)，选择满足约束条件(公式4)且目标函数值f^k最小的关阀策略为最优解。结合当前系统中阀门状态，确定当前阶段的最优阀门操作策略。如图3(b)所示，当前阶段k＝1的最优阀门操作策略为关闭v3和v4。

步骤6：执行最优阀门操作策略，将dma管网分为两个子区域。如图3(b)所示，关闭v3和v4将示例dma分解为两个子区域，其中一个包含s1、s2和s3，另一个包含s4和s5。

步骤7：任选一个包含dma入口的子区域实施在线水平衡分析以判断漏失位于哪个子区域内。对于图3(b)中两个子区域，任选其一进行水平衡分析，即通过dma入口和出口流量计和用户端智能水表分别计量进入所选子区域的流量和所选子区域内用户用水量，计算其非收益水量nrw。若所选子区域的nrw与dma的nrw相近，则说明漏失位于所选子区域内；否则，漏失位于另一个子区域内。为使漏失定位程序继续进行，假定水平衡分析结果判断出漏失位于包含s4和s5的子区域内。

步骤8：将判定不存在漏失的子区域内管道标记为无漏失(即x^k＝0)。如图3(b)所示，s1、s2和s3内管段被标记为无漏失管段。

步骤9：判断dma中是否有阀门操作可以进一步缩小漏失区域。若有，则回到步骤5继续进行下一阶段的漏失定位(即k＝k+1)；否则说明已找到最小漏失区域，则停止漏失定位程序。对于图3(b)所示情况，可以进一步缩小所示区域，回到步骤5继续进行阶段k＝2的漏失定位过程(最优阀门操作为打开v4，关闭v5)，直至找到最小漏失区域(即图3c中的s4)。

步骤10：根据漏失定位结果，派遣工作人员至所定位的漏失区域内采用人工听漏法等设备方法精确查找漏点位置，并修复所查找到的漏点。

下面结合实际应用场景来说明本发明的实施步骤和应用效果。以下应用场景仅用于说明本发明在具体实践中如何运用，但不用来限制本发明的范围。

图4为某城市供水管网中的一个dma区域，包含2个入口、1个出口、51个可操作阀门和总长为58.7km的管线。在该dma区域的入口和出口处均安装有远传流量计在线监测进入和流出dma的流量，同时各用户处也安装了智能水表以在线记录用户用水量。另外，在该dma区域内存在一个重要工业用户，要求不可中断供水(节点压力不能小于15.0m)。当地水务公司通过对该dma的夜间流量的持续监测分析进行漏失管理。当dma的最小夜间流量明显增加时，说明该区域内产生了新的漏失。下面以一个观测到的漏失事件(最小夜间流量由76.5l/s增加至85.2l/s)为例，应用本发明方法对其进行定位。

按照上述具体实施方式对实例dma中所产生的漏失进行定位，其具体的优化阀门操作和水平衡分析结果分别如图5和图6所示。首先，应用图论方法确定可将dma分解为两个子区域的关阀策略集(步骤1～3)。该dma管网图可由阀门分割为40个连通子图，运行contraction算法9,400,000次可到4940个管网图割集，从而形成具有4940个关阀策略的集合用于后续的优化目标函数求解。应用步骤4～9可逐步缩小漏失区域直至得到最小的漏失区域(图5d中红色管段，总长4.3km)。比如，在阶段1(时间t0～t1)，通过3次阀门操作(关闭3个阀门)将dma管网分为两个子区域(图5a)，对包含出口的子区域进行在线水平衡分析，结果显示其非收益水量nrw为8.5l/s左右(图6)。漏失定位前该dma的nrw约为9.0l/s，这说明漏失位于包含出口的子区域内，由此标记另一个子区域的管段为无漏失管段。类似步骤应用于后续阶段2～4(时间t1～t4)，水平衡分析均是在包含出口的子区域内进行。

最终，本发明通过应用4个阶段共计14次优化的阀门操作和相应水平衡分析可将漏失区域由58.7km缩小至4.3km，漏失定位效率提升约13倍。在确定最小漏失区域后，水务公司派遣听漏工作人员至该区域附近进行精确漏失定位，确认漏失位于所确定的区域内(如图5d所示)。因此，实例应用说明本发明所提方法显著提升了漏失定位的效率，并保证了漏失定位的准确性。

以上实施例只是用于帮助理解本发明的方法及其核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。