CN114239282B 一种在极小降压空间下实施供水管网dma分区的方法（太原理工大学）

(19)国家知识产权局(12)发明专利(10)授权公告号CN1142392(65)同一申请的已公布的文献号号李尚明张翀(74)专利代理机构北京恒创益佳知识产权代理事务所(普通合伙)11556专利代理师付金豹(54)发明名称检素与5相邻且水激方向激%的节点集SNO(57)摘要本发明公开了一种在极小降压空间下实施供水管网DMA分区的方法，该方法首先定义节点间相似度计算公式，接着通过MATLAB中的函数gamultiobj优化相似度公式中的系数从而确定分区边界管段。然后设置不同的最小服务水压并将它们分别作为约束条件，以分区后节点平均水龄和分区成本为目标，运用函数gamultiobj得到了多个边界管段上设备布置方案，通过三个原则从中筛选出最优布置方案。最后用模拟退火算法找出最佳的管段更换方案，使得分区后管网的水压可以满足规定的最小服务水压的要求。本发明提供的DMA分区方法能在顺利完成分区的基础上还不会对分区后管网的整体水质造成太大的影(57)摘要据。2A3:以谱聚类算法进行分区后的边界管段的数量及其平均流量、管径和长度为目标函A4:设置一系列不同的最小服务水压并将它们作为约束条件，以分区后节点平均水龄首先找到安装设备后管网中节点压力不满足规定最小服务压力的节点集S₁,然后找到和它们相邻且水流流向它们的节点集S₂,这些节点构成的管段集P就是拟更换管径的第一流流向它们的节点集S₃,然后更新管段集2.如权利要求1所述的一种在极小降压空间下实施供水管网DMA分区的方法，其特征在3.如权利要求2所述的一种在极小降压空间下实施供水管网DMA分4.如权利要求3所述的一种在极小降压空间下实施供水管网DMA分区的方法，其特征在3借助MATLAB中自带的求解多目标优化问题的函数gamultiobj确定三个参数的值；编码方式采用实数编码，决策变量即为相似度计算公式中的三个参数，在优化的过程区间，由此得到的规范化Laplacian矩阵在MATLAB计算特征值和特征向量时不出现虚数；管网不同，c值的选取也不相同，但对于同一管网的不同用水状态，c值的确定可参考上一用水状态，两者的差别很小或者保持不变。5.如权利要求1所述的一种在极小降压空间下实施供水管网DMA分区的方法，其特征在式中：M为管网中除水源外节点的总数目；t;为管网中节点i的节点水龄，h;T是边界管段上安装阀门的数量，个；Cvave,是边界管段上第v个阀门管段流量的列向量；Q是节点流量的列向量；L是管网的回路矩阵；h是管段水头损失的列向4决策变量为在边界管段上安装阀门或者是安装流量计，其维度等于边界管段的数量，6.如权利要求5所述的一种在极小降压空间下实施供水管网DMA分区的方法，其特征在7.如权利要求1所述的一种在极小降压空间下实施供水管网DMA分区的方法，其特征在5技术领域[0001]本申请属于城市供水管网设计技术领域，具体涉及一种对供水管网进行DMA分区的方法，尤其是一种在极小降压空间下实施供水管网DMA分区的方法。背景技术[0002]DMA分区是在边界管段关闭阀门和加装流量计将整个管网划分为若干个相对封闭的独立计量区域，通过对每个区域入口和出口流量的监控及在入口处进行压力调控，达到控制管网漏损和控制管理系统压力的目的。由于分区时要关闭边界阀门，不可避免的会造成管网水头损失的增大，从而降低了管网中各节点压力的冗余度，对于管网末梢降压空间极小的节点来说，这会导致它们分区后的水压很难满足最小服务水压，从而影响管网的正常运行。[0003]目前在DMA分区的研究中对管网末梢节点压力处理存在缺陷。周立典等在进行设备优化布置时以节点压力波动范围在10%以内为约束条件，这可能会使管网末梢的节点压力不满足规定。曾翰等将分区前水压超过16m的节点约束其水压在分区后不低于16m,在分区前低于16m的节点其水压在分区后不低于0m,显然这种约束方式给管网创造了较大的降压空间，很可能导致原本水压就很低的节点分区后水压更低。周中健等将《城镇供水服务》规范中规定的管网末梢压力不应低于14m应用于意大利Modena管网中，增大了管网的降压空间进而完成了DMA分区。发明内容函数gamultiobj结合，确定了理想的分区边界管段。然后通过设置一系列不同的最小服务水压，并将它们分别作为约束条件，运用函数gamultiobj得到边界管段上设备的最优布置方案。最后运用模拟退火算法找出需要更换管径的管段，使得分区后管网中所有节点的压力均可满足最小服务压力。[0005]本发明实现上述目的的技术方案如下。[0008]A2:指定分区数目，定义管网两节点间相似度的计算公式；[0009]A3:以谱聚类算法进行分区后的边界管段的数量及其平均流量、管径和长度为目标函数，通过函数gamultiobj优化相似度计算公式中的参数并确定分区边界管段；[0010]A4:设置一系列不同的最小服务水压并将它们作为约束条件，以分区后节点平均水龄和分区成本为目标函数，经gamultiobj优化计算得到Pareto最优解集，通过三个原则筛选出最优设备布置方案；[0011]A5:以更换管段的成本最小为目标，运用模拟退火算法找出最佳的管段更换方案。6分区结果。7[0023]对于管网末梢降压空间极小的节点来说，一旦管网中新关闭阀门，其节点压力就会很容易低于规定的最小服务水压，因此在优化边界管段上设备布置时，若是将规定的最小服务水压作为约束，很有可能会出现无解或解的质量很差的情况。因此设置一系列不同的最小服务水压作为优化设备布置的约束条件，观察各个压力约束下解的情况，确定最优布置方案。[0024]在进行边界管段上设备的优化布置时应考虑分区后管网的水质安全和分区的经济性，同时还必须满足节点连续性方程，能量守恒方程以及节点的水压约束。[0025]边界管段上优化设备布置的目标函数及约束条件为[0026]式中式中：M为管网中除水源外节点的总数目；t;为管网中节点i的节点水龄，h;T段上安装流量计的数量，台；Cmeter,是边界管段上第m台流量计的价格，元.A是管网的衔接矩阵；q是管段流量的列向量；Q是节点流量的列向量；L是管网的回路矩阵；h是管段水头损[0028]决策变量为在边界管段上安装阀门或者是安装流量计，其维度等于边界管段的数上安装阀门。[0029]根据以下三个原则从Pareto最优解集中选择出最优的布置方案。[0030]原则①每个DMA的入口最多为2个。[0031]原则②布置设备后的节点平均水龄不应高于分区前的水龄。[0032]原则③在满足每个DMA入口数量要求和节点平均水龄限制的基础上，选择运行后能使低于管网规定最小服务水压的节点数量达到最少的设备布置方案。[0033]确定设备布置方案后，建立在该方案下供水管网的水力模型(采用EPANET2.2软件)。[0034]对于不满足管网规定的最小服务水压的节点来说，可以放大流向该节点管段的管径以减少水流在转输过程中的能量损失，从而增大该节点的节点压力，使节点压力满足管网规定的最小服务压力。[0035]在更换过程中需要控制更换管段的数量、长度和管径以节省成本，更换后仍需满足管网的两大基本方程，且在更换完成后管网的节点平均压力不应大于更换前的平均压力，因为压力增大会导致管网中漏损量的增大。8[0040]首先找到安装设备后管网中节点压力不满足规定最小服务压力的节点集S₁,然后找到和它们相邻且水流流向它们的节点集S₂,这些节点构成的管段集P就是拟更换管径的第一批管段，采用模拟退火算法寻找使得目标函数最小且满足约束条件的管径更换方案，且水流流向它们的节点集S₃,然后更新管段集P得到拟更换管径的第二[0041]当首次找到管径更换方案时便可认为已经得到全局最优解是因为每一次继续向DMA分区的质量。附图说明[0047]图3为运用本发明的分区方法与其他研究提出的分区方法对案例管网进行分区后9[0049]图5为本发明中案例管网最终分区结果图。具体实施方式[0051]如图1所示，本发明的一种在极小降压空间下实施供水管网DMA分区的方法，包括如下步骤：[0053]步骤2,指定分区数目，定义管网两节点间相似度的计算公式；[0054]步骤3,以谱聚类算法进行分区后的边界管段的数量及其平均流量、管径和长度为目标函数，通过函数gamultiobj优化相似度计算公式中的参数并确定分区边界管段；[0055]步骤4,设置一系列不同的最小服务水压并将它们作为约束条件，以分区后节点平均水龄和分区成本为目标函数，经gamultiobj优化计算得到Pareto最优解集，通过三个原则筛选出最优设备布置方案；[0056]步骤5,以更换管段的成本最小为目标，运用模拟退火算法找出最佳的管段更换方[0057]下面以一具体实例为例进行说明。[0058]案例管网的拓扑结构如图2所示，该管网共有272个节点(包含4个水源节点)和317条管段，属于中型城镇管网。管网中最大管径为400mm,最小管径为100mm,每个节点有最低压力20m和最高压力32.67～44.11m的限制.在单工况条件下运行时，管网中节点压力最低为20.09m,且有13个节点压力没有超过20.5m,[0059]其他研究中也完成了对案例管网的DMA分区，在综合考虑了经济性和管理等因素后，确定分区数量都为4个。将其他研究得到的分区方案分别记作方案B和方案C。本发明对案例管网实施DMA分区采用的分区数目也设置为4个。方案B来源于现有技术“基于节点自然管网DMA分区方法研究[J].给水排水，2019,55(07):118-123节点能量冗余度的城市供水管网DMA自动分区研究”(参见：蒋浩.基于节点能量冗余度的城市供水管网DMA自动分区研究[D].广州：广东工业大学，2017.)。[0060]采用函数gamultiobj确定分区边界管段，参数设置如下：最优前端个体系数为0.3,种群大小为100,最大进化代数为600,停止代数也为600,适应度函数值偏差为0.01,其他参数为默认值。经试验得到c的取值可以是3.5,因此三个参数的取值限定在(0,3.5)之间。经MATLAB计算后得到Pareto最优解集，在此基础上确定了相似度计算公式中三个参数的值分别为α=0.198、β=2.41、γ=0.674并以此得到了边界管段和分区结果，将其记为方[0061]图3对三种分区方案的相关指标进行了对比，边界管段的流量之和、管径之和以及长度之和分别用fta、fta和f来表示。[0062]由图3可知，方案A中边界管段的平均流量和总流量均是其余两种方案的50%以下，而总的管径长度和平均管径长度均处于另外两种方案的55%～86%之间，边界管段的数量也是三者之中的最低值。[0063]确定分区边界管段后，进行边界管段上设备的优化布置，即在边界管段上安装阀门和流量计。[0064]设置不同的最小服务水压为20m,19m,18m和17m将它们分别作为约束条件，在MATLAB中导入不同管径规格下阀门和流量计的费用，通过gamultiobj函数优化分区边界管段上设备的布置方案。参数设置如下：最优前端个体系数0.1,种群大小100,最大进化代数为50,停止代数也为50,适应度函数值偏差为0.01,其他参数为默认值。[0065]经MATLAB计算后，在不同压力约束下得到的Pareto最优解个数分别为3,4,3,5,其中在压力约束为17m和18m时有2个解重合。表1罗列了15个Pareto最优解的两个目标函数值以及分区所需的流量计台数、布置设备后水压低于20m的节点个数n。和分区后各区域进水口数量。[0066]表1不同压力约束下的Pareto前沿及其他相关信息Tm/台原管网一0一03512034129033128533118331175321173211732116311152111622116221152111[0069]注：*表示该解在最小服务水压约束为17m时重复出现[0070]将不同压力约束下求得的Pareto最优解集按照以下三个原则进行筛选，1)单个DMA的进水口数量不得超过两个；2)管网的节点平均水龄不得高于0.72h;3)运行后水压低11于20m的节点数目最少。[0071]解s-11符合上述原则，其所对应的方案即为最优边界管段上设备布置方案。因此案例管网安装阀门8个，流量计6台，设备安装成本93928元，安装设备后节点的平均水龄为0.70h,单工况条件下运行后管网节点压力最低为17.61m,共有13个节点的压力低于20m。[0072]案例管网规定的节点最小服务水压为20m,分区后13个节点不满足要求，因此需要对管网中部分管段加大管径，减少管网水头损失，使所有节点的压力均可满足要求。[0073]在MATLAB中执行水力模拟，找到拟更换管径的管段集，并对更换管段的管径大小进行优化，通过模拟退火算法得出最优解。算法参数设置初始温度200℃,结束温度0.001℃,降温速率0.98,Mapkob链长50。编码的数字是从0～9的连续自然数，0代表该管段管径不[0074]在拟更换的管段集合更新两次后，计算停止。模拟退火算法运算过程如图4所示，在降温次数E为290次时算法收敛，即得到了全局最优方案。[0075]由计算结果可知需要更换的管段有10条，仅占管网管段总数的3.2%,更换管段的成本为847350元。将在边界管段上布置好设备的案例管网更换管径后，单工况条件下运行时所有节点的压力均大于规定的最小服务水压。最终分区的结果如图5所示。[0076]表2是Modena管网实施DMA分区前后各项运行性能指标的对比。表中，g为管网的综合水龄指数，它能合理反映整个管网水质情况，尤其是大流量用户和管网末梢区域的水表2DMA分区前后管网运行性能指标管网状态分区前分区后分区后的节点水龄相较于分区前没有发生变化，说明案例管网实施DMA分区后未对管网的水质造成太大影响。分区后管网的综合水龄指数下降，从这一角度表明分区后管网的水质尤其是管网末梢和大流量用