专利名称：基于实时水压信息的城市供水区域动态划分控制方法
技术领域：
本发明涉及的是一种水利监控技术领域的方法，具体是一种基于实时水压信息的城市供水区域动态划分控制方法
背景技术：
随着工业生产能力的提高，城市供水系统的规模在不断扩大，供水管网复杂性也随之提高。但是随之而来的是日益增多的各种各样的爆管事故，很多用户存在着水压不足的现象。根据中心城区管网监测数据，采用聚类方法，将这些测点依据其压力值和压力变化趋势进行了划分，确定平均水压的分区，如水压160kPa为危急低压区；160kPa<水压 < 180kPa为低压预警区；180kPa<水压< 300kPa为正常范围；300kPa <水压< 350kPa为高压预警区；水压> 350kPa为危急高压区。这样，在日常操作中就形成有侧重的监控，可以有效地提高供水效率和有效地减少爆管事故和水压不足等情况的发生。水压区域划分主要是对测压点、泵站、水厂在内的测点数据特性进行分析，采用聚类方法，将这些测点依据其压力值和压力变化趋势进行了划分，传统聚类方法可以分为两种监督式聚类和非监督式聚类。其中监督式聚类的为方法是K-Mean聚类，在这种方法中， 聚类的最终结果与用户初始指定的初始聚类中心有关。用户需要指定聚类数目、初始质心和评判某个成员归属某个聚类集合的阈值大小。这种方法本质上不是根据聚类元素的分布进行自适应聚类，且其聚类精度和效率依赖于初始质心选取等预设参数。非监督聚类中为性的聚类方法有支持向量聚类，是基于支持向量机的聚类方法。 在支持向量聚类中，聚类数据通过高斯核函数从低维空间被映射到高维特征空间，于是可以在这个高维特征空间找到一个最小的超球体将这些数据包围起来(虽然在低维空间无法找到这样的最小超球体将被聚类数据包围起来)，这个超球体接着被映射回原始低维数据空间，被看作是包围被聚类数据点的轮廓，这些轮廓就认为是聚类簇的边界，在同一个轮廓中被包围的数据点认为属于同一个聚类集合，从而完成聚类。但是，支持向量聚类不是基于聚类数据物理分布的，难以反映数据全局分布。近来，一种叫做AP 聚类(Affinity Propagation Clustering)，该算法是 B. J. Frey, Dueck D 于 2007 年发表在 science 上白勺 Clustering by passing messages between data points (基于数据点消息传递的聚类分析，Science, 2007，315 (5814) 972-976)。其基本思想就是通过消息传递，实现数据点的自动聚类。它将数据对象之间的相似度作为输入，而且此相似度可以是非对称的，即数据对象A到数据对象B的相似度可以不等于数据对象B到数据对象A的相似度。实值信息在数据对象之间交换传播直至一组高质量的聚类中心和相应的聚类产生。尽管该文称其聚类效果很好，计算速度也很快，但它也有几个缺点。第一，它也需要事先定义一个相似性度量，从而计算出数据点集之间的相似性矩阵来，这在时间和空间上就需要0 (n2)。第二，迭代次数需要人工设定，而且聚类结果对此也较敏感。第三，获得聚类结果后，不能获得聚类分布的层次性，有时这是不够的。K-Mean算法对于初始聚类中心的选择很关键，因为初始聚类中心选择的好坏直接影响到聚类结果，而且这个算法要求进行聚类时输入聚类数目，这也可以说是对聚类算法的一种限制。不过，这种算法运行速度相对于AP算法要快一些，因此，对于那些小而且数据比较密集的数据集来说，这种聚类算法还是比较好的。然而，在传统的K-Mean聚类技术中， 在某些情况下，不能获得预期的聚类结果，即，接近人类直觉的聚类结果。例如，在本文中， 分类结果不能很好满足相关地理信息。在基于以上各种聚类方法的特点分析，这里采用改进的K-Mean算法，S卩，在一次 K-Mean聚类下，考虑实际限制条件，进行二次聚类。首先，由于数据的充分性和对现实压力等数据的充分了解，可以运用工作经验，对水压区域划分的初始聚类中心进行很好的选择， 这样就充分利用了 K-Mean算法的特点，测压点、泵站、水厂等测点数据是小而且数据比较密集的数据集，因而可以发挥K-Mean算法的快速性特点。考虑到实际水压区域 划分建议既具有理论依据，又完全符合实际要求。本发明的划分结果便于事故分析和监测，一旦发生爆管等事故，影响区域基本可以确定，事故处理更有针对性。
发明内容
本发明针对现有技术存在的上述不足，提供一种基于实时水压信息的城市供水区域动态划分控制方法，通过充分利用测压点获取的数据信息，结合实际地理情况，对城市供水区域进行合理的分区，控制区块规模使之符合要求，便于管理。可广泛应用于大型城市供水调度系统的管理、监测、控制。本发明是通过以下技术方案实现的，本发明利用SCADA系统从供水管网采集测压点、水厂、泵站等节点的压力数据，然后利用数据挖掘手段对数据进行预处理，之后对压力数据样本进行系统聚类分析并根据聚类结果及区块规模进行区域规模指标聚类，最后利用地理连通图进行区块的合并和分割，形成最终划分结果。本发明具体包括以下步骤第一步、通过SCADA(SupervisoryControl And Data Acquisition，数据采集与监视控制)系统对供水管网压力节点进行数据采集，得到采集供水区域内各压力节点的压力数据。所述的数据采集是指含有压力的变化特性的时间域宽度。时间域的选择应当具有广泛性，包含工作日与休息日，尤其以用水量较大的时段为宜，通常跨度范围不少于3个月，每个月采样的天数不少于10天。第二步、将采集到的压力数据进行无效数据剔除处理，然后计算每个压力节点的数据变化值并作为该压力节点的附加属性；所述的无效数据剔除处理是指对所有压力数据进行时间指标归一化，然后依次选取所有压力节点在任一时刻T的压力数据，当在该时刻T没有压力数据时，则进行数据校正和/或坏点剔除。所述的数据校正是指判断当任一时刻T的前后给定误差时间范围At内是否存在压力数据，当存在则选取T+At和T-At时刻的压力数据并求平均后作为T时刻压力节点的压力数据；
所述的坏点剔除是指当连续10个采样时刻无数据时，则判定该时刻的压力数据为坏点，相应删除该压力节点及其对应的压力数据。所述的附加属性通过以下方式得到对校正及坏点剔除后的数据，计算每个压力节点T时刻的压力值与T-I时刻的压力值之差，将所有T时刻与T-I的压力值差值作为该压力节点的附加属性。 第三步、以供水区域内的所有压力节点及其相互之间管线直接联通与否作为元素值建立地理连通图矩阵，具体为对供水区域内的η个压力节点Pi，i = 1，2，. . .，η，η为自然常数建立压力节点管线的地理连通图矩阵Anxn，其中地理连通图矩阵中元素~在压力节点Pi与P」之间存在管线直接联通时取1，否则为0。第四步、每一个压力节点称为一个样本，对所有样本进行聚类计算，得到类间距离平均值及类内距离平均值。所述的聚类计算的步骤如下4. 1)先将η个样本分为η类，计算每一类的中心压力值CPi (i = 1，2，…ncn) (ncn 为当前聚类数):CL 二」-^其中m为一个压力节点类中的压力节点数，Xi (i
=1,2,…m)为每个压力节点的压力特性(包括压力值和压力变化值)；4.2)找出距离最近的两个类并将其合并成一个新类，重新计算新类的中心压力值，并计算该次聚类过程的类间距离平均值和类内距离平均值；4. 3)重复步骤4. 1和步骤4. 2，直到所有的压力节点聚为一类。第五步、根据类间距离平均值和类内距离平均值选择较优聚类数目。所述的较优聚类数目，是指当聚类数从M减为M-I时，反映类内聚合性优劣的类内距离平均值出现超过5%的减小，或反映类间分离性优劣的类间距离平均值出现超过 5%的增加时的聚类数目。第六步、根据聚类的规模，进行区域规模指标聚类，包括步骤如下6. 1)指定类的最大规模，通常根据城市规模、人口密度等设定在5-10内，选取出超过规模的类进行区域规模指标聚类；6.2)单独取出超过规模的类中的压力节点，依照K-Mean聚类法，将其分为两类， 判断新划分的两类是否规模都符合给定要求当符合要求，则用新划分的两类替换原先的一类，完成区域规模指标聚类；当不符合要求，依次增加分类数，重复上述过程，直到所有分类均符合给定规模。第七步、对第三步中得到的地理连通图矩阵进行校正处理，并根据聚类结果，按类对所有压力节点分成若干不同的分区。所述的校正处理是指针对某一个已经确定的类，随机选取其中一个压力节点作为起始节点，以地理连通图作为压力节点之间的连通关系，依照图的遍历方法，判断是否能够遍历该类中所有的节点；当有压力节点无法遍历，则将不能遍历的压力节点取出，组成新类，重复上述过程，直到每一个类中所有压力节点都能够遍历。第八步、利用中心线法确定第七步中的分区界限作为最终分区结果并在城市供水监控、管理系统中加载，实现基于实时水压信息的供水区域动态划分，并实时显示区域划分图，当发生供水问题，比如某压力节点爆管，该节点及所在分区产生警示信号，并通过后台数据库调取该区域面积及人口等信息，作为事故评级和应急调度的评价因素。所述的中心线法是指选择相邻的两个区域，计算每个区域的分类中心，确定它们连线的中点，依次遍历所有两两相邻的区域，确定中点，对所有邻近中点进行连线，得到最终的区域划分结果。 本发明得到的聚类结果，类的规模可以根据实际情况予以调控针对人口密集地区，可以将类的规模控制的较小，管理更为细致；针对人口较分散地区，可以将类的规模设置的较大，节约管理所需的人力、物力和财力。每个类中的压力节点在压力值和变化趋势方面都表现出较强的相关性，便于事故分析和监测。一旦发生爆管等事故，影响区域基本可以确定，事故处理更有针对性。
图1为本发明的实施步骤流程图，图2为本发明实施例中的最终加载结果示意图。
具体实施例方式下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。实施例如图1所示，本实施例实施步骤如下1.通过SCADA系统对供水管网压力节点进行数据采集。某城市东部地区，共有4 大水厂，10条出水管道，因而对应10个压力节点，共有5个泵站，43个测压点。共计58个压力节点。对这58个压力节点进行编号，然后采集各个压力节点从2009年6月1日0时0 分至2009年9月30日23时40分，每20分钟一个压力数据。这样，每一个压力节点就得到了共计8784个压力值数据。在每年的6月到9月，是属于较为炎热的夏季，因此该时间段是用水量较大，压力信息变化较明显的时间段。另一方面，这期间包括了工作日，周末休息日以及暑假，涵盖范围广泛。2.将采集到的压力数据读入计算机并进行预处理，首先搜索第一个压力节点的所有压力值数据，找出其中不满足采集条件，既在某采样时刻附近允许时间误差内未获取压力值的时刻，去除该时刻所有压力节点的压力值数据。然后依次搜索其他节点，直到完成所有节点的数据剔除工作。此时，每个节点符合采集条件的压力数据共计7328个。然后计算每个压力节点压力值的变化数值，共计7327个，将其作为压力节点的附加属性。这样就构成了包含58个压力节点，每个压力节点有14655个压力属性的矩阵。其中每一行的编号为压力节点的编号，每一列为压力节点在某个相同时刻的压力属性。3.依据地理信息，找出每两个压力节点的连通情况，构建地理连通图。4.首先将58个压力节点每个单独作为一类，计算每一类的中心压力序列，以及两两之间的距离。选择距离最近的两个压力节点合并为一个新类并计算该新类的中心压力序列。同时计算该次聚类过程中的类间距离平均值(ado)和类内距离平均值(adi)。之后重复上述过程，直到所有的压力节点聚为一类。聚类信息见表1。表 权利要求
1.一种基于实时水压信息的城市供水区域动态划分控制方法，其特征在于，利用 SCADA系统从供水管网采集测节点的压力数据，然后利用数据挖掘手段对数据进行预处理， 之后对压力数据样本进行系统聚类分析并根据聚类结果及区块规模进行区域规模指标聚类，最后利用地理连通图进行区块的合并和分割，形成最终划分结果。
2.根据权利要求 1所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述划分控制方法具体包括以下步骤第一步、通过数据采集与监视控制系统对供水管网压力节点进行数据采集，得到采集供水区域内各压力节点的压力数据；第二步、将采集到的压力数据进行无效数据剔除处理，然后计算每个压力节点的数据变化值并作为该压力节点的附加属性；第三步、以供水区域内的所有压力节点及其相互之间管线直接联通与否作为元素值建立地理连通图矩阵；第四步、每一个压力节点称为一个样本，对所有样本进行聚类计算，得到类间距离平均值及类内距离平均值；第五步、根据类间距离平均值和类内距离平均值选择较优聚类数目；第六步、根据聚类的规模，进行区域规模指标聚类；第七步、依据第三步中得到的地理连通图矩阵对聚类结果进行校正处理，并根据聚类结果，按类对所有压力节点分成若干不同的分区；第八步、利用中心线法确定第七步中的分区界限作为最终分区结果并在城市供水监控、管理系统中加载，实现基于实时水压信息的供水区域动态划分，并实时显示区域划分图。
3.根据权利要求2所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的无效数据剔除处理是指对所有压力数据进行时间指标归一化，然后依次选取所有压力节点在任一时刻T的压力数据，当在该时刻T没有压力数据时，则进行数据校正和 /或坏点剔除。
4.根据权利要求3所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的数据校正是指判断当任一时刻T的前后给定误差时间范围At内是否存在压力数据，当存在则选取Τ+Δ t和Τ-Δ t时刻的压力数据并求平均后作为T时刻压力节点的压力数据；所述的坏点剔除是指当连续10个采样时刻无数据时，则判定该时刻的压力数据为坏点，相应删除该压力节点及其对应的压力数据。
5.根据权利要求2所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的附加属性通过以下方式得到对校正及坏点剔除后的数据，计算每个压力节点 T时刻的压力值与T-I时刻的压力值之差，将所有T时刻与T-I的压力值差值作为该压力节点的附加属性。
6.根据权利要求2所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的第三步具体为对供水区域内的η个压力节点Pi，i = 1，2，...，η，η为自然常数建立压力节点管线的地理连通图矩阵Anxn，其中地理连通图矩阵中元素％在压力节点 Pi与P」之间存在管线直接联通时取1，否则为0。
7.根据权利要求2所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的聚类计算的步骤如下~4.1)先将η个样本分为η类，计算每一类的中心压力值CPi (i = 1,2, "mcn),其中ncn为当前聚类数:CL =」-^ m为一个压力节点类中的压力节点数，Xi (i = 1，~2，…m)为每个压力节点的压力值和压力变化值；~4. 2)找出距离最近的两个类并将其合并成一个新类，重新计算新类的中心压力值，并计算该次聚类过程的类间距离平均值和类内距离平均值；~4. 3)重复步骤4. 1和步骤4. 2，直到所有的压力节点聚为一类。
8.根据权利要求2所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的较优聚类数目，是指当聚类数从M减为M-I时，反映类内聚合性优劣的类内距离平均值出现超过5%的减小，或反映类间分离性优劣的类间距离平均值出现超过5%的增加时的聚类数目。
9.根据权利要求2所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的第六步具体包括步骤如下~6. 1)指定类的最大规模，通常根据城市规模、人口密度等设定在5-10内，选取出超过规模的类进行区域规模指标聚类；~6. 2)单独取出超过规模的类中的压力节点，依照K-Mean聚类法，将其分为两类，判断新划分的两类是否规模都符合给定要求当符合要求，则用新划分的两类替换原先的一类， 完成区域规模指标聚类；当不符合要求，依次增加分类数，重复上述过程，直到所有分类均符合给定规模。
10.根据权利要求2所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的校正处理是指针对某一个已经确定的类，随机选取其中一个压力节点作为起始节点，以地理连通图作为压力节点之间的连通关系，依照图的遍历方法，判断是否能够遍历该类中所有的节点；当有压力节点无法遍历，则将不能遍历的压力节点取出，组成新类，重复上述过程，直到每一个类中所有压力节点都能够遍历。
11.根据权利要求2所述的基于实时水压信息的城市供水区域动态划分控制方法，其特征是，所述的中心线法是指选择相邻的两个区域，计算每个区域的分类中心，确定它们连线的中点，依次遍历所有两两相邻的区域，确定中点，对所有邻近中点进行连线，得到最终的区域划分结果。
全文摘要
一种水利监控技术领域的基于实时水压信息的城市供水区域动态划分控制方法，利用SCADA系统从供水管网采集测压点、水厂、泵站等节点的压力数据，然后利用数据挖掘手段对数据进行预处理，之后对压力数据样本进行系统聚类分析并根据聚类结果及区块规模进行区域规模指标聚类，最后利用地理连通图进行区块的合并和分割，形成最终划分结果。本发明通过充分利用测压点获取的数据信息，结合实际地理情况，对城市供水区域进行合理的分区，控制区块规模使之符合要求，便于管理。可广泛应用于大型城市供水调度系统的管理、监测、控制。
文档编号G05B13/00GK102385313SQ20111016511
公开日2012年3月21日 申请日期2011年6月17日 优先权日2011年6月17日
发明者仇军, 汪瑞清, 王景成, 葛阳, 赵平伟, 赵金涛 申请人:上海交通大学, 上海市供水调度监测中心