一种供水管网压力检测点的布置方法
【专利摘要】本发明提出了一种供水管网压力检测点的布置方法，包括以下两个步骤：基于节点压力对管网进行分区，将位置上相邻的并且压力值相近的节点划分为同一个供水子区；和对每个供水子区分别运用灵敏度系数法对管网压力检测点进行选址，以从所述供水子区内的节点中选出压力检测点。本发明的供水管网压力检测点的布置方法，可以满足大型城市供水管网的施工需求，压力监测点的布置位置的计算过程简单，并可以简化计算量。
【专利说明】
一种供水管网压力检测点的布置方法
技术领域
[0001] 本发明涉及给排水工程技术领域，尤其涉及一种供水管网压力检测点的布置方 法。
【背景技术】
[0002] 众所周知，城市供水管网系统中的压力监测点一般指利用有线或者无线进行数据 实时传输的压力记录仪。在管网中设置压力监测点是为了实时的监测城市供水管网的运行 状态。在满足用户用水压力的情况下，根据监测点处的压力值，合理的调节栗站水栗的运行 方式，以降低管网漏水减少栗站的耗能。因此，在城市供水管网系统中的合适位置布置压力 监测点，能够科学合理的调整管网的运行状态，满足居民的用水需求，进而对提高供水的安 全性和可靠性有重大的意义。
[0003] 目前，国内外对有关供水管网压力检测点的选址方法做了一些富有成效的研究， 比较具有代表性的有灵敏度系数法和基于遗传算法的目标函数优化方法。然而，基于灵敏 度系数法的压力监测点布置方法，需要首先求出每个节点相对于其它节点的灵敏度系数， 然后利用聚类方法将灵敏度值相近的节点归为一类，分别选择每一类中最大灵敏度值所对 应的节点作为压力监测点。然而上述方法的缺点是计算量太大，例如，对于节点数目为10万 的大型城市供水管网，需要进行10万次平差计算以得到10万X 10万的灵敏度系数矩阵。如 此巨大的计算量及数据存储量使得该方法难以应用。而遗传算法也有缺陷，当供水管网的 规模较大时，采用遗传算法计算，计算时间较长难以进行实际应用。
[0004] 因此，需要一种新型的供水管网压力检测点的布置方法，以满足大型城市供水管 网的施工需求。
【发明内容】
[0005] 为了解决上述部分或全部技术问题，本发明提供了一种供水管网压力检测点的布 置方法，包括步骤:基于节点压力对管网进行分区，将位置上相邻的并且压力值相近的节点 划分为同一个供水子区；和对每个供水子区分别运用灵敏度系数法对管网压力检测点进行 选址，以从所述供水子区内的节点中选出压力检测点。
[0006] 进一步地，对管网进行分区包括数据准备步骤（1):A，按照管道的编号顺序，分别 计算每根管道相邻节点的压力差Ah1P并将计算的结果作为权值赋予该管道;B，设定分区 阀值S;C，制定分区结果数据格式，将分区计算结果存放在矩阵中，所述矩阵中的每一行表 示一个子区；以及D，初始化矩阵，将矩阵内的所有元素设为零。
[0007] 进一步地，节点压差Ahij为：|hi-hj I，当节点i与节点j相连接时;或0,当节点i与 节点j为同一个节点时;或00，当节点i与节点j不相连接时。
[0008] 进一步地，对管网进行分区还包括计算步骤(2)至(5): (2)按照管道编号顺序从管 道k=l开始搜索；（3)判断该管道k的起始节点i是否存在于矩阵中，若节点i已经存在于矩 阵中，则转到下一个步骤(4);若节点i不存在于矩阵中，则将节点i赋给矩阵中的一个新行； (4)判断矩阵中是否已经存在管道k的终止节点j;若存在，则继续判断所述分区阈值δ是否 大于节点j所在行中所有节点之间的的Ahu，若满足，则合并节点i与节点j所在的行，若不 满足，则转到步骤(5);若不存在，则继续判断阈值δ是否大于节点i与节点j所在行中所有节 点的Ah 1P若满足，则将节点j加入到节点i所在的行中，若不满足，则将节点赋给矩阵中的 一个新行；以及(5)对下一个管道k = k+l重复步骤(3)和步骤(4)直到遍历所有的管道为止。
[0009]进一步地，对管网压力检测点进行选址的过程包括:定义节点之间的压差相似系 数矩阵并对矩阵进行标准化;对所述标准化后的压差相似系数矩阵，利用欧式距离法计算 节点间的相似程度系数，得到模糊相似矩阵；以及分别计算节点与其余节点的平均欧式距 离，选取最小平均欧式距离的节点为压力检测点。
[0010]进一步地，压差相似系数矩阵为由叫组成矩阵Xnxn;
[0011]
[0012]其中Hl、Hj分别为正常状态下，节点i与j的压力值;Hmax与H min分别为正常状态下，管 网中节点压力的最大与最小值;η为管网中节点的总数。
[0013] 进一步地，对矩阵进行标准化包括:对矩阵X进行标准差变换得到矩阵X7 ;和对矩 阵V进行极差变换得到矩阵X"。
[0014] 进一步地，矩阵f由Xi1/组成；
[0015:
[0016]其中，&表示矩阵X中对第k列元素取平均值，
；Sk表示矩阵X中第k 列元素的标准：
η为管网中节点的总数。
[0017] 进一步地，矩阵Χ〃由Xik〃组成；
[0018
[0019] 其中：X 'kmir^示矩阵X'中对第k列元素取最小值;X '表示矩阵X'中对第k列元素 取最大值;〇 < X" 1; η为管网中节点的总数。
[0020] 模糊相似矩阵为由相似程度系数rij组成的矩阵Rnxn;其中rij计算公式为：
[0021]
[0022] 其中x"lk，x" 分别为第i行、第j行各元素;η为管网中节点的总数。
[0023] 进一步地，压力监测点位置的确定需要分别计算节点与其余节点的平均欧式距离
选取最小平均欧式距离的节点为测压点。
[0024] 本发明的供水管网压力检测点的布置方法，可以满足大型城市供水管网的施工需 求，压力监测点的布置位置的计算过程简单，计算量小。
【附图说明】
[0025] 在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中：
[0026] 图1为本发明实施例的供水管网分区计算过程的流程图；
[0027] 图2为本发明实施例的供水管网模型图；
[0028] 图3为本发明实施例的供水管网模型图的分区图；以及
[0029] 图4为本发明实施例的供水管网模型图的压力监测点布置结果图。
[0030] 在附图中，相同的部件使用相同的附图标记。附图并未按照实际的比例。
【具体实施方式】
[0031] 下面将结合附图对本发明作进一步说明。
[0032] 在本发明的一个具体的实施例中，如图2所示，为一个实际供水管网的管网模型 图。该管网模型由6个用水节点1-1、1-2、1-3、1-4、1-5以及1-6，一个水源水库和一个水栗组 成。利用管网平差软件EPANET对管网模型进行计算，计算得到的每个节点的压力，如表所示 1〇
[0033]表1节点压力计算结果
[0035]在该实施例中，首先对管网按照压力采用如图1所示的计算流程进行分区，具体的 计算过程如下：
[0036] (1)数据准备。按照管道的编号顺序，分别计算每根管道相邻节点的压力差Ahlj, 并将计算的结果作为权值赋予该管道;设定分区阈值δ = 2.0，初始化矩阵RESULT。相应的计 算结果如表2所示。
[0037] 弄2管道相郃书点压力差ALi+笪结里
[0039]其中，初始化后的矩阵RESULT为：
[0040]
[00411 (2)按照管道顺序，从管道Pl开始搜索，若Pl起始节点J-I没有在矩阵RESULT中，则 将J-I赋给RESULT-个新行;若Pl终止节点J-2没有在矩阵RESULT中，J-I与J-2之间的压力 差Δ hi」=1.21，小于允许误差δ，因此将J-2加入到J-3所在的行，此时RESULT为：
[0042]
[0043]同理，对于管道?2，？3，？4，？5，？6，？7，按照上述图1中的流程进行计算，最终得到分 区计算结果RESULT，此时RESULT为：
[0044] V u V
^ v v v/
[0045] 从结果矩阵RESULT中可以看出，管网被分为2个区域。分区结果如图3所示。
[0046] 接下来对每个区域设定一个压力监测点，以分区1为例，计算过程如下：
[0047] 由干压差相M赛翁钽陈为由Yd日成钽陈XnXn;
[0048]
[0049] 其中H1A分别为正常状态下，节点i与j的压力值;Hmax与!1_分别为正常状态下，管 网中节点压力的最大与最小值;η为管网中节点的总数。
[0050] (1)计算两两节点的压差相似系数，将表1中的数据带入压差相似系数技术公式 中：
[0051；
[0052；
[0053]
[0054]
[0055]
[0056]
[0057]
[0058]
[0059]
[0060]
[0061]
[0062]
[0063]
[0064]
[0065]
[0066] 同理，对矩阵X的第2列、第3列分别进行计算，得到标准化后的两两节点的压差相 似系数矩阵：
[0067]
[0068]
[0069]
[0070]同理，矩阵R中其余元素按照上述过程进行计算，计算结果为：
[0071：
[0072:
[0073:
[0074:
[0075:
[0076]最后，选取最小平均欧式距离的节点为测压点，即选择节点I，即J-I作为区域1的 压力监测点。
[0077]同样地，可以计算区域2的压力监测点。压力监测点的计算结果如图4所示。
[0078]虽然已经参考优选实施例对本发明进行了描述，但在不脱离本发明的范围的情况 下，可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是，只要不存在结构冲 突，各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文 中公开的特定实施例，而是包括落入权利要求的范围内的所有技术方案。
【主权项】
1. 一种供水管网压力检测点的布置方法，其特征在于，包括如下步骤： 基于节点压力对管网进行分区，将位置上相邻的并且压力值相近的节点划分为同一个 供水子区；和 对每个供水子区分别运用灵敏度系数法对管网压力检测点进行选址，W从所述供水子 区内的节点中选出压力检测点。2. 根据权利要求1所述的压力检测点的布置方法，其特征在于，所述对管网进行分区包 括数据准备步骤(1): A，按照管道的编号顺序，分别计算每根管道相邻节点的压力差Ahu，并将计算的结果 作为权值赋予该管道； B，设定分区阀值δ; C，制定分区结果数据格式，将分区计算结果存放在矩阵中，所述矩阵中的每一行表示 一个子区；W及 D，初始化矩阵，将矩阵内的所有元素设为零。3. 根据权利要求2所述的压力检测点的布置方法，其特征在于，所述节点压差Δ hu为： llu-lij I,当节点i与节点j相连接时；或 0，当节点i与节点j为同一个节点时;或 °°，当节点i与节点j不相连接时。4. 根据权利要求2所述的压力检测点的布置方法，其特征在于，所述对管网进行分区还 包括计算步骤(2)至(5): (2) 按照管道编号顺序从管道k= 1开始捜索； (3) 判断该管道k的起始节点i是否存在于矩阵中，若节点i已经存在于矩阵中，则转到 下一个步骤(4); 若节点i不存在于矩阵中，则将节点i赋给矩阵中的一个新行； (4) 判断矩阵中是否已经存在管道k的终止节点j; 若存在，则继续判断所述分区阔值δ是否大于节点j所在行中所有节点之间的Ahu，若 满足，则合并节点i与节点j所在的行，若不满足，则转到步骤(5); 若不存在，则继续判断阔值S是否大于节点i与节点j所在行中所有节点的Ahu，若满 足，则将节点j加入到节点i所在的行中，若不满足，则将节点赋给矩阵中的一个新行；W及 (5) 对下一个管道k = k+l重复步骤(3)和步骤(4)直到遍历所有的管道为止。5. 根据权利要求1所述的压力检测点的布置方法，其特征在于，所述对管网压力检测点 进行选址的过程包括： 定义节点之间的压差相似系数矩阵并对矩阵进行标准化； 对所述标准化后的压差相似系数矩阵，利用欧式距离法计算节点间的相似程度系数， 得到模糊相似矩阵；W及 分别计算节点与其余节点的平均欧式距离，选取最小平均欧式距离的节点为压力检测 点。6. 根据权利要求5所述的压力检测点的布置方法，其特征在于，所述压差相似系数矩阵 为由：?^组成矩阵ΧηΧη ;其中出、田分别为正常状态下，节点i与j的压力值;Hmax与Hmin分别为正常状态下，管网中 节点压力的最大与最小值;η为管网中节点的总数。7. 根据权利要求5所述的压力检测点的布置方法，其特征在于，所述对矩阵进行标准化 包括： 对矩阵X进行标准差变换得到矩阵公；和 对矩阵公进行极差变换得到矩阵X"。8. 根据权利要求7所述的压力检测点的布置方法，其特征在于，所述矩阵r由xii/组成；其中，而表示矩阵X中对第k列元素取平均值：;Sk表示矩阵X中第k列元 素的标准差，η为管网中节点的总数。9. 根据权利要求7所述的压力检测点的布置方法，其特征在于，所述矩阵X"由xik"组 成；其中：X'k min表示矩阵X'中对第k列元素取最小值;X'k max表示矩阵X'中对第k列元素取 最大值;〇《义"11<《1;11为管网中节点的总数。10. 根据权利要求9所述的压力检测点的布置方法，其特征在于，所述模糊相似矩阵为 由相似程度系数ru组成的矩阵Rnxn;其中ru计算公式为：其中A分别为第i行、第j行各元素;η为管网中节点的总数。
【文档编号】G06F17/50GK105938505SQ201610228129
【公开日】2016年9月14日
【申请日】2016年4月12日
【发明人】李国斌, 曾宪灵, 曾锋, 周扬, 陈圣鹏, 冯炎辉
【申请人】广州京维智能科技有限公司