一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法技术

本发明专利技术公开了一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法。依据节点灵敏度分布的特点，在备选节点集合中筛除两类不利于进行水力模型校核的节点。主要步骤如下：(1)求解供水管网节点的灵敏度矩阵；(2)计算每个节点的总灵敏度与变异系数；(3)根据总灵敏度与变异系数分别对节点进行排序，在备选节点集合中筛除总灵敏度低与变异系数高的节点。本发明专利技术提出的备选节点筛选方法能有效降低搜索空间，提高监测点布置算法的效率。提高监测点布置算法的效率。提高监测点布置算法的效率。

【技术实现步骤摘要】
一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法


[0001]本专利技术属于城市供水管网类，具体是一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法。

技术介绍

[0002]作为城市公共设施的一个重要组成部分，供水管网系统将饮用水输送到千家万户，以满足人民的生产生活需求。水力模型可以用于提升供水管网运行管理水平的。为了提高供水管网水力模型的精度，就需要以一定频率对其进行校核。而水力模型校核所依托的水力数据来自于管网中布置的监测点。因此，监测点布置方案的好坏对于最终水力模型校核的精度有很大的影响。然而，供水管网模型中节点的数量通常是十分庞大的，如果直接进行穷举搜索，搜索空间十分庞大，计算时间是不可接受的。因此，对节点进行初步筛选，选出备选节点，降低搜索空间是十分必须的。
[0003]鉴于此，本专利技术提出一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法，有效降低搜索空间，提高监测点布置算法的效率。

技术实现思路

[0004]本专利技术旨在提出一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法。为实现以上目的，本专利技术采取以下技术方案：
[0005]一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法，其特征在于包括如下步骤：
[0006]步骤(一)、求解供水管网节点的灵敏度矩阵；
[0007]步骤(二)、计算每个节点的总灵敏度与变异系数；
[0008]步骤(三)、根据总灵敏度与变异系数分别对节点进行排序，在备选节点集合中筛除总灵敏度低与变异系数高的节点。
[0009]在采用上述技术方案的基础上，本专利技术可以采用以下进一步的技术方案，或对这些进一步的技术方案组合使用：
[0010]所述步骤(一)，具体为：
[0011]假定供水管网模型中有n个节点，已知水厂的总供水量，依据管长，将需水量分配到n个节点，得到节点的基准需水量X＝[x1，x2，...，x
n
]T
，其中X
i
为第i个节点的基准需水量；将节点的基准需水量代入到模型中，对模型进行管网水力平差计算，得到所有节点的基准压力H＝[H1，H2，...，H
n
]T
，其中，H
i
为第i个节点的基准压力值；i＝1,2,3,4
……
n；
[0012]对于节点i，对其需水量进行扰动Δx
i
，得到扰动后的需水量同时保持其它节点需水量不变，即扰动后的n个节点的需水量为将扰动后的节点需水量替换管网模型所有节点的需水量，然后进行水力平差计算，得到扰动后
的节点压力
[0013]通过公式(1)计算节点j压力对节点i压力的灵敏度值S
ij
，
[0014][0015]其中，H
i
与H
j
分别为第i个与第j个节点的基准压力值，与分别为第i个与第j个节点扰动后的节点压力。
[0016]一个节点对自身压力的灵敏度定义为1，即，
[0017]S
ii
＝1(2)
[0018]遍历所有的节点i＝1，2，...，n,对于每一个节点都分别扰动需水量，并根据公式(1)与公式(2)计算该节点对所有节点的灵敏度值，得到，
[0019][0020]进而得到管网中所有节点的灵敏度矩阵S；
[0021][0022]S中共有n个列向量，对应n个节点，每个列向量有n个元素，代表着每个节点对管网中所有n个节点压力的灵敏度。
[0023]所述步骤(二)，具体为：
[0024]采用公式(5)计算每个节点对其他所有节点的总灵敏度TS：
[0025][0026]式中，n为管网模型中节点的总数，S
j，i
为节点i的压力值对节点j的节点需水量的灵敏度；
[0027]采用公式(6)计算每个节点灵敏度的变异系数C
v
：
[0028][0029]式中，σ为数据集的标准差，μ为数据集的平均值。
[0030]所述步骤(三)，具体为：
[0031]得到所有节点的TS值后，依据该数值从大到小对所有节点进行排序。计α1为“灵敏度数值筛除因子”，将所有节点中总灵敏度排在后α1的节点从备选节点集合中移除；对每个节点的灵敏度向量计算变异系数C
v
，同样依据该数值从大到小对所有节点进行排序。定义α2为“变异系数筛除因子”，则将所有节点中变异系数排在前α2的节点从备选节点集合中移除。
[0032]为了阐述步骤(一)至(三)的原理，以下结合图1
‑
3对筛除的基本原理与目的进行说明。绘制节点的灵敏度向量的概率密度分布图，可以大致看出该节点对管网中所有节点
的灵敏度分布情况，一个H市管网中某个位于管网中游传感器灵敏度的累积概率密度分布如图1所示。通过此曲线可以看出，该传感器在不同区域的节点需水量发生变化时，自身压力的变化也是不尽相同的。依此做出推测，比较适宜布置压力传感器的节点，其灵敏度的分布最好满足两个特征：(1)对其他节点灵敏度的数值不能过低。如果监测点所在节点对其他节点的灵敏度过低，在压力传感器存在一定数值的监测噪声误差时，管网中其它节点的需水量发生变化所引起监测点的变化会被传感器的监测误差掩盖，无法对水力模型校核的精度提高有实质性的作用；(2)灵敏度的分布应当较为均匀。为了达到更高的校核精度，监测点的监测数值变化要能对区域中各个节点的需水量变化有更细化的感知。如果一个节点的灵敏度分布是倒L形(即对少数节点的灵敏度极高)，则在校核时只能降低少数节点的偏差，而很难对整个管网其他位置的节点需水量进行校准。上述两个特征的反例分别如图2与图3所示，图2为H市管网中某个位于水厂附近主干管上节点灵敏度的累积概率密度分布图，图3为H市管网中某个位于管道末梢节点灵敏度的概率密度分布图。
[0033]正是基于以上原理，本专利技术通过计算节点的总灵敏度TS以及变异系数Cv，为了降低监测点布置的搜索空间，对这两类不利于节点需水量校核的布置点进行筛除，有效降低搜索空间，提高监测点布置算法的效率。
附图说明
[0034]图1为H市管网中某个位于管网中游传感器灵敏度的累积概率密度分布图；
[0035]图2为H市管网中某个位于水厂附近主干管上节点灵敏度的累积概率密度分布图；
[0036]图3为H市管网中某个位于管道末梢节点灵敏度的概率密度分布图；
[0037]图4为H市供水管网图；
[0038]图5为H市供水管网模型中第一类不利节点位置示意图；
[0039]图6为H市供水管网模型中第二类不利节点位置示意图。
具体实施方式
[0040]下面结合附图和实例，对本专利技术的实现方式进一步详细叙述，具体步骤如下：
[0041]步骤(一).获取供水管网节点的灵敏度矩阵。
[0042]如图4所示，H市供水管网有水厂3个，需水节点4242个(n本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法，其特征在于，包括如下步骤：步骤(一)、求解供水管网节点的灵敏度矩阵；步骤(二)、计算每个节点的总灵敏度与变异系数；步骤(三)、根据总灵敏度与变异系数分别对节点进行排序，在备选节点集合中筛除总灵敏度低与变异系数高的节点。2.如权利要求1所述的一种服务于供水管网水力模型校核的监测点布置备选节点筛选方法，其特征在于，所述步骤(一)，具体为：假定供水管网模型中有n个节点，已知水厂的总供水量，依据管长，将需水量分配到n个节点，得到节点的基准需水量X＝[x1，x2，...，x
n
]
T
，其中X
i
为第i个节点的基准需水量；将节点的基准需水量代入到模型中，对模型进行管网水力平差计算，得到所有节点的基准压力H＝[H1，H2，...，H
n
]
T
，其中，H
i
为第i个节点的基准压力值；i＝1，2，3，4......n；对于节点i，对其需水量进行扰动Δx
i
，得到扰动后的需水量同时保持其它节点需水量不变，即扰动后的n个节点的需水量为将扰动后的节点需水量替换管网模型所有节点的需水量，然后进行水力平差计算，得到扰动后的节点压力通过公式(1)计算节点j压力对节点i压力的灵敏度值S
ij
，其中，H
i
与H
j
分别为第i个与第j个节点的基准压力值，与...

【专利技术属性】
技术研发人员：周华，俞亭超，林犇，申屠华斌，邵煜，翁晓丹，江鑫，
申请(专利权)人：浙江大学，
类型：发明
国别省市：